热油管道防腐层大修期间热力参数的数值模拟的研究.docVIP

热油管道防腐层大修期间热力参数的数值模拟研究 李伟1(宇波1 张劲军1 张争伟1 黄金萍2 （1．中国石油大学（北京）城市油气输配技术北京市重点实验室，北京昌平102249；2．中油管道分公司沈阳调度中心，辽宁沈阳，110031） 摘要：基于对防腐层大修期间施工进度与输油工况的合理简化，建立了有开挖段站间热油管道正常运行数学模型，并采用有限容积法对该模型进行数值离散。通过对该模型的数值求解，可以模拟防腐层大修期间不同环境条件、不同运行参数和开挖回填不同阶段管道沿线原油和土壤温度变化。根据温度模拟结果，可以确定大修期间开挖站间的最大允许开挖长度。利用所编计算程序对2006年铁秦输油管道葫芦岛－绥中站间防腐层大修期间运行工况的模拟结果显示，进站油温Figure1: sketch map of buried pipeline 油流的换热方程： （1） 结蜡层、管壁和防腐层的导热方程： （2） 土壤导热方程： （3） 边界条件： 当时 （4） 当时， （5） 当或时， （6） 当，；或，时， （7） 其中，为原油密度，kg/m3； 为时间，s；为油流平均速度，m/s； 为重力加速度，m/s2； 为油流截面平均压力，Pa；为达西摩阻系数；为原油定压比热容，J/(kg·℃)；为管道内直径，m；为原油在单位管壁面积上单位时间的散热量，w/m2；为原油温度，℃；为原油膨胀系数，℃-1；为第层（结蜡层、管壁和防腐层）的密度，kg/m3；为第层（结蜡层、管壁和防腐层）的比热容，J/(kg·℃)；为第层（结蜡层、管壁和防腐层）的温度，℃；为第层（结蜡层、管壁和防腐层）的导热系数，W/(m·℃)；为径向位置，m；为环向弧度；为土壤密度，kg/m3；为土壤比热容，J/(kg·℃)；为土壤温度，℃；为土壤导热系数，W/(m·℃)；为垂直于轴向的水平位置，m；为深度，m ；为油流对管内壁的放热系数，W/(m2·℃)；为管内壁温度，℃；为地表向大气的放热系数，W/(m2·℃)；为大气温度，℃；为自然热力状态下的土壤温度，℃；为热影响区域的边界至管中心的水平距离，m；为热影响区域的边界至地表的垂直距离，m；为管道埋深，m；为管道最外层半径，m。 对站间开挖段，管外介质变为空气，其边界条件处理为： 在没有天气突变、泡水及回填的情况下，认为开挖段管道直接和空气进行对流换热，根据钢管的导热量和管外壁与大气对流换热量的平衡关系列方程。 回填后，假定管道周围回填土壤温度均匀分布，初始温度为当时土壤自然温度。 ⑶ 在有泡水的情况下，根据不同的泡水深度进行处理: ① 当泡水完全淹没开挖段输油管道时，对泡水区域进行网格划分，用水的物性参数求解泡水带走的热量。 ② 当泡水段液面在管道下方，没有触及管道时，按照管道外壁和大气对流换热来处理。 ③ 当泡水淹没管道一部分，并没有完全淹没的时候，根据上述两种方法按照泡水区域所占的比例来进行处理。 2．数学模型的数值离散 采用有限容积法对模型中二维非稳态导热方程进行离散化，首先对管道周围矩形热力影响区域进行区域离散。图2为求解区域内网格划分示意图，图3为该区域内任意节点的控制容积示意图。针对任意的控制容积，导热方程（3）的积分形式如下： （8） 其中V为控制容积的体积（对二维导热问题为控制容积的面积），A为控制容积界面的面积矢量，其正方向与外法线单位矢量一致，如图2所示，符号“”表示两个矢量的内积。将上式应用于如图3中所示三角形控制容积，可得 （9） 式中为重心为P0的三角形的面积，和分别为时间间隔的当前时层和上一时层P0点的温度值。是界面1，2，3上的平均温度梯度。界面上的平均温度梯度可以通过节点上的温度梯度线形插值得到： （10） 其中，和为插值因子。 为了确定离散方程，采用最小二乘法来确定温度梯度如下： （11） 其中，表示P0节点的温度梯度在i坐标轴上的分量，为从P0到Pj的有向线段。