!天然气管道大孔泄漏速率的建模与分析_陈平

收稿日期：2013-09-14；定稿日期：2013-11-01基金项目：中央高校基本科研业务费专项资金资助项目（FRF-TP-12-067A）作者简介：陈平（1973-），女，北京人，副教授，博士。主要研究方向为材料摩擦磨损、工程图学。E-mail：[e-mail]目前世界上95%以上的天然气都采用管道输送方式[1]，但是输气管道在运行过程中可能会受到各种因素的干扰，致使管道发生泄漏事故[2-3]。根据泄漏孔径的大小，通常将管道泄漏模型分为小孔模型、管道模型和大孔模型3种类型[4-5]。国内外学者针对上述3种不同失效模式下的气体泄漏率的计算模型展开了一系列研究，并取得了可喜的成果[3-4,6-12]。但是对于大孔泄漏模型的研究仍然不太成熟，比如不同压力下的状态分析与计算模型等，本文结合上述文献，对上面的3种模型的适用范围进行了研究，对大孔泄漏模型泄漏速率的变化进行了理论分析，并且对大孔泄漏过程中存在的3种状态在高压和低中压情况下进行详细的分析讨论，得出更加符合实际天然气管道大孔泄漏速率的建模与分析陈平，梅华锋，赵海林（北京科技大学机械工程学院，北京100083）摘要：通过对天然气管道泄漏过程的分析，建立了气体的大孔泄漏模型。对大孔泄漏过程存在的3种状态进行讨论分析，验证了在高压情况下只存在管内亚临界流、孔口临界流和管内孔口都是临界流两种情况，在低中压情况下只存在管内亚临界流、孔口临界流和管内孔口均为亚临界流两种情况，并且分析了不同泄漏点和不同压力情况下发生泄漏的速率变化情况。采用本文的计算模型可以较好的计算出不同泄漏孔径的泄漏速率，为长输管道的定量风险评价提供依据。关键词：管道；泄漏速率；大孔泄漏；建模；分析中图分类号：TE873文献标识码：A文章编号：2095-302X(2014)03-0486-04ModelingandAnalysisfortheBig-holeLeakageRateofGasPipelineChenPing,MeiHuafeng,ZhaoHailin(SchoolofMechanicalEngineering,UniversityofScienceandTechnologyBeijing,Beijing100083,China)Abstract:Throughtheanalysisontheprogressoftheleakageofgaspipeline，thecalculatingmodelsofbig-holeleakageareestablished.Basedontheleakageprocessofbig-hole，thethreekindsofstatearediscussedandanalyzed.Itistestifiedthat,underthesituationsofhigh-pressure,thereonlyexisttwopossiblestates,whicharesub-criticalflowinpipeandcriticalflowinorifice,orcriticalflowinpipeandorifice.Whileinthecaseoflowandmediumpressure,thereonlyexistssub-criticalflowinpipe,criticalflowinorificeandsub-criticalflowinpipeandorifice.Therateofleakageatdifferentpointsisalsoanalyzed.Theleakyrateofdifferentholediameterscanbewellcalculatedbythemodelsinthepaper,anditprovidesbasisforquantitativeriskassessmentoflonggaspipeline.Keywords:pipeline;leakyrate;big-holeleakage;modeling;analysis2014年6月图学学报June2014第35卷第3期JOURNALOFGRAPHICSVol.35No.3情况的泄漏计算模型，并且通过实例验证其正确性，可为长输管道的定量风险评价提供依据。1基本方程图1为管道气体泄漏示意图[4]。如图1所示，管道泄漏点处距离输气起始点的间距为Le，1为起始点处中心位置，2为与泄漏孔在同一截面上的管内某点，3为泄漏口处位置，Pi、ρi、Ti、ωi（i=1、2、3）分别表示1、2、3点处的气体压力、密度、温度、气体流速，Pa、ρa、Ta分别表示外界环境的压力、密度和温度。图1管道气体泄漏示意图泄漏孔口处泄漏率的大小取决于泄漏孔处的流动状态，其判别标准由临界压力比CPR值[13]来确定，即：1221aPCPRP(1)式中，Pa：大气环境的压力，Pa（绝）；P2：气体在点2处的压力，Pa（绝）；γ：泊松比。（1）当2/aPPCPR时，孔口气体为临界泄漏，气体泄漏率m为：112221gdhWmCAPZRT(2)式中，m：气体泄漏率(kg/s)；Cd：流量系数，在雷诺数Re30000的非临界流条件下取0.61，其他情况取1.0；Ah：泄漏孔口面积(m2)；Wg：气体的摩尔质量(kg·mol-1)；Z：压缩因子；R：理想气体常数，为8.314(J·mol-1·K-1)；T2：气体在点2处的温度(K)。（2）当2/aPPCPR≥时，孔口气体为亚临界泄漏，气体泄漏率m为：21222221gaadhWPPmCAPZRTPP(3)2大孔泄漏模型大孔泄漏模型是一种孔径介于小孔和管道之间的模型，Montiel等[7]认为在泄漏过程中存在3种可能的情况：管内亚临界流、泄漏孔口临界流状态；管内、泄漏孔口亚临界流状态；管内、泄漏孔口临界流状态。董玉华等[3-4]认为第三种情况不存在。而王兆芹等[14-15]则认为在高压和低压情况下管内、泄漏孔口的状态会有所不同。结合上述文献，下面分别对3种情况进行详细的论述。（1）管内亚临界流、泄漏孔口临界流状态：当气体的运动满足式(4)的条件时，气体为管内亚临界流、泄漏孔口临界流状态：212112211aMaPPMaPCPRP(4)气体泄漏率的计算公式由式(2)来表示。管内点1和点2各参数间的关系式(5)表示[6,16]：2122122121212222212121212121212121MaTTMaMaPMaPMaMaMaMaMaMa(5)式中，Mai：点i处的马赫数(i=1,2)；ρi：点i处的气体密度kg/m³(i=1,2)；Ti：点i处的气体温度(K)(i=1,2)。（2）管内、泄漏孔口亚临界流状态：当气体的运动满足式(6)的条件时，气体在管内、泄漏孔口均为亚临界流状态。212112211aMaPPMaPCPRP≥(6)气体泄漏率的计算公式由式(3)来表示。管内点1和点2各参数间的关系由式(5)表示。（3）管内、泄漏孔口临界流状态：当气体的运动满足式(7)的条件时，气体在管内和泄漏孔口均为临界流状态。212112211aMaPPMaPCPRP(7)第3期陈平等：天然气管道大孔泄漏速率的建模与分析487气体泄漏率的计算公式由式(8)来表示：1111121gdhWmCAMaPZRT(8)管内点1和点2各参数间的关系由式(5)表示。3计算案例与分析某天然气长输管道长度为872km，管道内径为660mm，管道内壁粗糙度为46µm，管线的起始点压力为5MPa，起始点气体温度为293K。天然气的摩尔质量为17.1kg/kmol，压缩因子为0.9，黏度为1.01×10-5Pa·s，理想气体常数为8.314J·mol-1·K-1，取泊松比为1.3，管道在距离起始点等效长度Le=126km处发生泄漏。假设气体在泄漏过程中，管道起始点的压力，即点1处的压力P1为定值。根据文献中的公式计算小孔泄漏、管道泄漏情况下的泄漏速率以及采用上述讨论的大孔泄漏模型得出的计算结果见图2~6。泄漏速率（kg/s）图2不同模型的气体泄漏速率与孔径比的关系曲线图3小孔模型和大孔模型泄漏速率比值与孔径比的关系曲线图4大孔模型和管道模型泄漏速率比值与孔径比的关系曲线从图2~4可知，当管道的泄漏孔径较小，即在其孔径比（d/D）小于0.15时，小孔模型和大孔模型计算所得的结果相近，相对误差小于5%，可以用小孔模型代替大孔模型来计算。随着泄漏孔径的不断增大，在孔径比（d/D）大于0.9时，大孔泄漏模型所得计算结果与管道模型相近，相对误差小于5%，可以用管道模型代替大孔模型来计算。图5不同起始压力下的气体泄漏速率与孔径比的关系曲线（Le=126km）从图5可以看出，随着起始压力的增大，同一孔径比所对应的气体泄漏速率也逐渐增大，在相同的起始压力下，气体的泄漏速率随着孔径比的增大而逐渐增大，最终趋于平缓。在P1≤0.184MPa时，泄漏孔口处的泄漏状态一直为亚临界流状态；在0.184MPa＜P1＜7MPa时，当孔径比增大到某一点时，泄漏孔口处的泄漏状态将由临界流状态转变为亚临界流状态，如P1=5MPa时，在孔径比为0.87（泄漏孔径大约为0.57m）左右时气体的泄漏状态发生了变化；在P1≥7MPa时，当孔径比从0逐渐增大到1的过程中，气体的泄漏一直处于临界流状态。从而验证了上面对大孔模型的论述的正确性，即在高压时只有管内为亚临界流状态、泄漏孔口为临界流状态和488应用与交流2014年管内、泄漏孔口均为临界流状态两种情况，低中压时只有管内、泄漏孔口均为亚临界流状态和管内为亚临界流状态、泄漏孔口为临界流状态两种情况。0.00.10.20.30.40.50.60.70.80.91.0孔径比(d/D)10080604020050km800km126km200km400km(a)P1=5MPa0.00.10.20.30.40.50.60.70.80.91.0孔径比(d/D)14012010080604020050km800km126km200km400km(b)P1=7MPa图6不同泄漏点距离下的气体泄漏率与孔径比的关系曲线从图6可知，在压力一定的情况下，泄漏点处离起始点距离越远，同一孔径比下的泄漏速率就越小，这主要是管道内壁摩擦致使气体压力下降所造成的。从图中还可以看出当P1=5MPa时，虽然泄漏点距离不同，但是气体在刚开始时均处于临界流状态，随着孔径比的增大，泄漏状态都转变为亚临界流状态；当P1=7MPa时，无论泄漏距离为多少，气体都一直处于临界流状态。这也从另一方面验证了前面大孔模型讨论的正确性。4结论（1）当管道的泄漏孔径较小，即在其相对孔径比（d/D）小于0.15时，小孔模型和大孔模型计算所得的结果相近，误差小于5%，可用小孔模型来代替大孔模型计算气体的泄漏速率；当泄漏孔径接近管径，即在其相对孔径比（d/D）大于0.90时，大孔模型计算所得结果与管道模型趋于相等，误差小于5%，可用管道模型来代替大孔模型；在处于中间范围时，则只能采用大孔模型来计算气体的泄漏速率。（2）对于大孔模型，在高压情况下，只存在管内亚临界流、孔口临界流和管内、孔口均为临界流两种情况；在低中压情况下，只存在管内、孔口亚临界流和管内亚临界流、孔口临界流两种情况。（3）在压力一定的情况下，泄漏点处离起始点距离越远，同一孔径比下的泄漏速率就越小，而且不同孔径下的泄漏速率变化不同。另外在不同的压力情况下，气体泄漏所处的状态也会不一样。参考文献[1]霍桂霞.管线球阀的设计特点研究[J].宁夏机械,2010,1(2):19-21.[2]FearnehoughGA.Thecontrolofriskin