[发明专利]一种大型原油浮顶储罐非稳态传热过程数值模拟方法

说明书

本发明涉及一种大型原油浮顶储罐非稳态传热过程数值模拟方法，解决了现有数值模拟研究无法全面反映相关因素对大型浮顶油罐温度场的影响规律的问题。具体采用网格生成技术建立的大型浮顶原油储罐内油品传热过程的数学模型，能够最大限度地保证数值模拟过程与实际工况相吻合；并运用迭代的方法对该模型进行求解，其具有简单、灵活以及通用性强的特点，容易在计算机上实现，既能保证计算效率，又能节约时间。

技术领域

本发明涉及一种数值模拟方法，特别涉及一种大型原油浮顶储罐非稳态传热过程数值模拟方法。

背景技术

随着原油储备需求的增加，大型浮顶油罐因其技术和经济上的优势已成为大规模原油库的首选储油设施。油罐在生产运行中，当罐内原油温度低于析蜡点时，含蜡原油会因析蜡在浮顶、罐壁和罐底所包围的内沿形成一定厚度和强度的凝油层，严重时可能发生凝罐等事故。因此，为确保储罐安全生产，必须准确掌握罐内原油温度场的变化规律。国内外针对储罐油品温度场的研究方法主要集中在数值模拟、现场测试两个方面。采用现场测试，对特定的工程问题能够得到较可靠的结果，但其只能在条件相同的情况下使用，具有较大的局限性，而应用传热学理论对油罐进行数值模拟，则可以较好地弥补这一不足。

应用数值模拟方法研究油罐温度场的变化规律，由传统的经验公式计算到现今的软件模拟求解，建立的物理数学模型与实际油罐越来越贴近，相应地数值求解过程也越来越复杂。采用经验公式法计算罐内原油温度场变化规律时，对边界传热系数的设定往往采取用经验值的方法，将罐内原油温度看成一个整体进行计算，计算过程虽简单便捷，但将实际情况进行了过度简化，导致计算结果精度有限，在实际工程中只能作为参考使用。现有商业软件，如Fluent等CFD数值模拟软件可以用于原油储罐非稳态传热过程的数值模拟，其具有丰富的物理模型、先进的数值方法和强大的前后处理功能，但是计算耗时较多，仅能反映对象较短一段时间内的温降特点，不易应用于工程实践，并且这类软件不能得到整体温度场的变化规律。

综上所述，目前针对大型原油浮顶储罐非稳态传热过程数值模拟研究，均具有一定的局限性，无法全面反映相关因素对大型浮顶油罐温度场的影响规律。

发明内容

为了克服现有技术的上述缺点，本发明提供了一种大型原油浮顶储罐非稳态传热过程数值模拟方法。

本发明解决其技术问题所采用的技术方案是：一种大型原油浮顶储罐非稳态传热过程数值模拟方法，包括如下内容：

步骤一：以原油储罐过轴线的剖面为研究对象，建立储罐二维非稳态传热数值模型；

步骤二：运用有限差分法对数值模型区域进行离散，得到外部角点、边界节点以及内部节点的非稳态离散方程，并根据离散方程的稳定性，确定空间步长以及时间步长；

步骤三：根据罐顶、罐壁以及罐底的边界条件，建立原油物性参数、边界传热系数与温度的对应关系，所述边界传热系数包括罐顶传热系数、罐壁传热系数及罐底传热系数，具体按照以下步骤进行：

(1)假定罐顶温度，与油品温度之间的平均值为定性温度，然后计算出在此定性温度下油品的物性参数：相对密度、粘度、导热系数、比热容以及体膨胀系数；

(2)根据格拉晓夫准则和普朗特准则，确定中间系数和值，然后计算出内部放热系数、外部放热系数以及辐射放热系数，如果是双盘浮顶罐，还需要根据浮舱定性温度计算出舱内放热系数；

(3)计算得出罐顶传热系数，根据热平衡原理进行验算；如果满足精度，则认为假定的平均温度是准确的；如不满足，应重新假定平均温度，再计算罐顶传热系数，同理可得罐壁传热系数及罐底传热系数；

步骤四：对罐内原油初始温度进行赋值，运用插值法将相应温度下的原油物性参数及边界传热系数，代入到外部角点、边界节点以及内部节点离散方程中，然后迭代进行下一个时间步长的计算，得出罐内原油温度场分布，直到整个模拟过程计算结束，得出罐内原油温度场的变化规律。