大落差天然气管道清管冲击分析-天然气与石油.DOC

基金项目：国家自然科学基金面上项目。 大落差天然气管道清管冲击分析 Dropping impact analysis of hilly natural gas pipeline pigging 李长俊1,2刘浠尧1贾文龙1彭阳1 （1西南石油大学石油与天然气工程学院，成都610500 2 CNPC复杂天然气集输重点研究室，西南石油大学，成都 610500） 摘要：中缅天然气管道沿线地势复杂，形成了很多的大落差管段，最大处可达1000多米，给清管作业带来了相当大风险。针对某大落差管段，采用OLGA软件分析了清管器在管内的运行速度与持液率变化，得出管道承受的冲击载荷，再利用CAESARII软件分析了在此荷载作用下管道应力与位移变化。结果表明，清管器在作业时，管道内部的输送压力对大落差管段的冲击荷载影响明显；管道在土壤中的横向和纵向位移均较小；在积液量较少情况下，清管器速度与持液率的增加对冲击应力影响较小。 关键词：大落差管段；清管；冲击荷载；OLGA；CAESARII 1.引言 从云南瑞丽市进入中国境内的天然气管道，在清管作业时，液弹在清管器推动下经过大落差管段会产生较大的冲击荷载，极易导致管道失效。国内学者多以实验与数值模拟相结合研究清管器在管道中的运行规律[1,2]，以及沿坡敷设输气管道应力[3]。但关于大落差管道在清管作业时液体对管道冲击影响的研究[4]，国内却较少开展。根据中缅天然气管道落差较大，里程较长的特点，分析了某大落差管段的应力与位移影响，以评估清管作业时管道的安全性。 2.基于OLGA大落差管道清管分析 2.1大落差管段模型建立 以云南省内某段中缅大落差埋地管道为例，利用OLGA软件[5]建立的清管模型如图2-1所示，其沿线高程如图2-2所示。 图2-1 中缅天然气某大落差管段清管仿真模型 图2-2中缅天然气某大落差管段沿线里程与高程 根据该管线实际出口压力6470kPa，入口流量1420万方/天，总积液量500m3，可以得到清管器运行速度和管道内持液率的变化。 2.2清管器运行速度变化 设定管道稳定运行1.5h后，从首站投入清管器，模拟清管过程中清管球的运动速度，跟踪清管球的运动位置，结果如图2-3所示。 图2-3清管球运动速度与运行位置分析 由图2-3可以 REF _Ref406793609 \h \* MERGEFORMAT 看出，清管器最大运行速度为8.0m/s。特别地，清管器到达22#阀室后，清管器的速度出现了较大的波动，说明液体主要积聚在22#阀室至末站的大落差管段间。因此，对该管段持液率进一步分析。 2.3管内持液率变化 分别模拟了5种不同时间（12h，25h，26h，31h，32h）从22#阀室至末站管内持液率变化情况，如图2-4~2-8所示。 图2-4 清管开始12h后22#阀室至末站持液率情况（清管器未到22#阀室） 图2-5 清管开始25h后22#阀室至末站持液率情况（清管器已达到22#阀室） 图2-6清管开始26h后22#阀室至末站持液率情况 图2-7清管开始31h后22#阀室至末站持液率情况 图2-8清管开始32h后22#阀室至末站持液率情况（清管器已到达末站） 可知，当清管进行26h后管内的持液率介于0.4~1.0之间；当清管进行31h后，管道持液率最大为1.0；当清管进行32h时，管线内的积液被清出，几乎无积液。 2.4冲击荷载计算 对液体通过弯头，进行受力分析，如图2-9所示。由动量和能量方程得到液体对管道的冲击荷载为： （2-1） 式中Rx，Ry为液柱对管道在轴向和纵向上的推力，N；v为清管球速度，m/s；A为管道的过流面积，m2；α为弯头与水平坐标轴之间的角度；ρ为清管液柱的密度，kg/m3；hL为管道坡底弯头处的持液率。 图2-9 清管过程液体对管道的受力分析图 3.基于CAESARII大落差管段冲击分析 3.1大落差管段模型建立 22#阀室至末站管段全长12.08km，高程差为1016.36m，该管段清管工况分析模型如图3-1所示，管段局部放大模型如图3-2所示。沿线管段参数和土壤参数，见表3-1和表3-2。 表3-1 管道参数 材料 外径/mm 直管壁厚mm 弯管壁厚mm 腐蚀裕量mm 介质密度 kg/m3 保温层厚度/mm 安装温度 ℃ 操作温度 ℃ X70钢 1016 17.5 22.5 0.1 95 0 20 50 表3-2 土壤参数 摩擦系数 土壤密度 （kg/m3） 到管顶埋深/m 摩擦角度/（°） 压实系数 热膨胀系数 0.6 0.001686 1.6 0 3 0.015 图3-1 22#阀室至末站某大落差段模型图 3-2 22#阀室至末站某大落差段局部放