管道泄漏及放空计算教材.docx

根据一元气体流动基本方程式，推导了孔口泄漏在绝热过程下泄漏 流量计算的小孔模型和适合管道完全断裂的多变过程泄漏流量计算 的管道模型，联合两种模型计算任何泄漏孔口直径下的泄漏流量， 讨论了燃气最大泄漏流量的限制，进行了实例计算并对比了不同模 型的计算结果关键词：泄漏流量计算；管道模型；小孔模型；管道小孔综合模型; 流量限制Calculation of Leakage Rate from Gas PipelineHUANG Xiao-mei， PENG Shini， XU Hai-dong， YANG Mao-huaAbstract: According to the basic equations of one-dimensio n gas flow， a hole model for calculation of hole leakage r ate in adiabatic process and a pipeline model for calculat ion of leakage rate in variable process suited to full rup ture of pipeline are deducted. These two kinds of models a re combined to calculate the leakage rate from leakage hol es with different diameters. The limitation of the maximum gas leakage rate is discussed，the example calculation is carried out， and the calculation results of different mode ls are compared.Key words: calculation of leakage rate； pipeline model； ho le model； combined model of pipeline model and hole model： limitation of flow rate1概述在燃气管道事故定量风险评价、事故抢险预案制定和漏气损 失评估时，首先要计算泄漏流量。

燃气管道在事故破损时，燃气可通过两种途径进入到大气中，一种是燃气直接泄漏到大气环境中， 另一种是泄漏到土壤中，通过土壤渗透进入大气环境前者可以通 过理论推导得出泄漏流量的计算公式，后者理论计算比较复杂且不 确定性很大本文主要分析和讨论前一种情况下的泄漏流量计算第三方破坏是城市燃气管道泄漏的主要原因之一，其主要表现是挖 掘机器、钻孔机器破坏管道，在这种情况下，燃气通常直接泄漏到 大气中此外，架空管道泄漏也是直接泄漏到大气中2小孔模型的推导管道泄漏示意图见图1小孔模型是将泄漏孔口当作孔径很 小的小孔，从而建立泄漏流量计算的模型图中点1——管道起点点2——泄漏口入口点点3——泄漏口出口截面上的点泄漏点至管道起点的距离，mqvu 泄漏点上游管道体积流量，ms/hqv——泄漏体积流量，ms/h图1中，点1通常为该管道上游的调压器出口，其压力通常保 持不变假设点4的断面流量及其平均流速方向不受泄漏影响，而 点4下游至泄漏口处的任何点管道断面平均流速由于受到泄漏影响 而不再沿管道轴线方向，点4至点2的距离非常小，可以忽略不计， 因而点4的压力近似等于点2的压力小孔模型假设管内燃气全部从该小孔泄漏，即管道上游无支 管或支管燃气流量为0,这样假设是为了保证从小孔泄漏的燃气流 量是最大值；由于泄漏小孔孔径较小，泄漏流量有限，因而忽略管 道沿程阻力，认为泄漏处的管内压力等于管道起点压力，即：p2=pi (1)式中p2——图1中点2的绝对压力，PaPi——图1中点1的绝对压力，Pa在泄漏孔处，燃气流速一般较快，燃气没有足够的时间与环 境进行热量交换，因此燃气泄漏过程，即从点2到点3的燃气流动 过程可被视为可压缩气体绝热流动过程，可见泄漏孔口与喷嘴相似。

孔口泄漏瞬间的流动可以看作是一维流动，气体的一元流动欧拉运 动微分方程为［1］：业 + vdv = 0 (2)式中p——燃气绝对压力，Pa燃气密度，kg/ms燃气断面平均流速，m/s因为泄漏过程为绝热过程，所以有⑵：f = G (3)式中K ——燃气的等熵指数C1——常量K是温度的函数，在常温下理想气体的K可近似当作定值 ⑵，对于天然气等由多原子分子组成的气体，K取1.29因为城市燃气压力大多在1.6MPa以下，燃气温度为常温， 所以燃气可以看作理想气体⑵，因此有：,5 (4)式中Rcon——燃气的气体常数，J/(kg - K)T——燃气的温度，K将式⑶代入式⑵并在小孔入口点2和小孔出口点3进行积 分，然后将式(4 )代入整理可得：式中V3——点3的燃气断面平均流速，即为燃气的泄漏出口流速，m /sv2l——点2沿管道泄漏口轴线方向的流速，m/sT2——点2处的燃气温度，Kp3——点3的绝对压力，PaP3可按式(6〜8)计算⑵：式中Pa——环境压力，PaPc 临界压力，PaP ——临界压力比在小孔模型下，V2L 一般很小或近似为0,因而V2 l2远远小于V32；, 可以忽略不计，因此式(5)可简化为：泄漏孔出口压力与入口压力的比值等于临界压力比时的泄漏出口燃气流速称为临界流速，将式(7)、(8)代入式(9)有：式中vc——泄漏出口燃气临界流速，m/s当p3达到临界压力pc后，点3的流速V3就等于临界流速秽，。

根据孔口出流的质量流量公式，可得出泄漏质量流量为：q =0.25p n d2p 3v3 (11)式中qm 泄漏质量流量，kg/s流量系数，可取0.90〜0.98泄漏孔口当量直径，m气——点3处的燃气密度，kg/m3不规则孔口当量直径按下式计算:d=—\ 7T式中气一一泄漏孔口面积，m2将式（4）、（6）、（9）、（10）代入式（11）并整理可得:当"时：E•心匡tiiETi (⑵当广项时:Pi吝(13)在小孔模型下，认为式（1）成立，且由于管内流速小，管内 流动可视为等温过程，因而对于任一管道，只要知道管道起点压力、 管内燃气温度和泄漏孔口当量直径，便可按照式（12）和（13）计算泄 漏质量流量3管道模型的推导管道模型适合于管道完全断裂的情形⑶，即泄漏当量直径等 于管道内径，点2即为管道末端，点2和点3重合，管道泄漏流量 就等于管输流量，此时可按管道水力计算公式来计算管道流量燃气泄漏过程中，管道上游阀门关闭之前，管内燃气流动可 以视为稳定流动在燃气管道完全断裂泄漏情况下，燃气流速较大， 管内燃气可能没有足够的时间和周边环境进行充分的热交换，管内 流动也不能看作等温过程，而只能看作多变过程，因而对于管内燃 气流动过程有：式中n——多变指数C2——常量仍然假设从管道起点至泄漏点之间的管段上没有支管，或所 有支管的流量都为0，同时考虑到管内燃气流速较大时，不能忽略 对流项（九），将文献［4］提供的管内流动基本微分方程组联合式（14） 可以推导得到式（15）。

式中qmu——泄漏点上游管道内燃气的质量流量，kg/sD——管道内径，m入一一摩擦阻力系数，可按文献［4］或［5］提供的公式计算式（15）中n=1时，表示管内流动为等温过程，此时认为管内 燃气与周边环境有充分的热交换，管内温度等于环境温度且保持不 变；当n=K时，表示管内流动为绝热过程，此时认为管内流速太快 或者管道太短，管内燃气完全没有和环境进行热交换；实际上这两 种理想状况都不存在，n的值在1和K之间，为了简化计算，通常 在管内燃气流速较小，管道较长时，n取1,而在管道流速很大或者 管道很短时，n取k在管道模型下，点2为管道末端且暴露于大气环境中，因而 泄漏口处的压力等于大气环境压力，即：P2=Pa （16）此时只要知道管径、泄漏处至管道起点的距离、管道起点的 压力，选取适当的值，便可按照式（15）计算管道的质量流量楫广4管道-小孔综合模型当燃气管道完全断裂时，按管道模型计算泄漏流量比较准 确，当燃气管道只有孔径很小的破损孔时，按小孔模型计算泄漏流 量比较准确但在实际情况下，特别是由于施工开挖导致的断裂， 泄漏口既不是小孔，也不是完全断裂，因此用这两种模型都不准确 管道一孔综合模型就是将管道模型和小孔模型结合起来进行泄漏流 量计算（九）。

事实上，在小孔模型下，只要有燃气泄漏，管道内燃气就会 流动，由于摩擦的存在，管道内必然会有沿程阻力，因此式（1）并不 成立，实际的p2还需要根据式（15）计算前面已经假设管输燃气全 部从小孔泄漏，这种假设对于风险评价和事故应急抢险来说是保守 的，因而有：qm=qmu （17）式（17）中qm可根据式（12）和（13）计算，式（12）和（13）中的 咒可由式（4）和（14）联立求得，qmu可以根据式（15）计算，因此联立 式（4）、（12）、（13）、（14）、（15）、（17）便可计算出泄漏质量流量q m、泄漏处的管内压力p2及管内温度T2该模型对于任何泄漏口尺寸 的稳定泄漏过程都比较准确，计算时需要利用计算机语言编程求解为了符合工程习惯，通常需要将泄漏质量流量转换成泄漏体 积流量：(h = 3 600— (18)P()式中Po——标准状态下燃气密度，kg/m35最大泄漏流量限制的讨论燃气在管道内的最大流量有限制，并不是按照管道模型或者 管道一小孔模型计算得到的结果就一定是实际流量，而必须对计算 结果进行两个方面的校验：一是调压器最大通过能力的限制，二是 管内临界流量的限制① 调压器最大通过能力的限制城市燃气管道都与上游调压器相连，管道内的流量和泄漏流 量都不能大于上游调压器的最大通过能力qvmax，调压器说明书通常 都给出了该调压器在各进口压力下的qvmax值，因此无论用哪种模型 计算，泄漏流量的计算结果最大只能取气呻。

② 管内燃气临界流量的限制与泄漏孔处的绝热流动相似，燃气管内流动也存在临界流量 问题由于沿程摩擦阻力对燃气运动参数的影响，实际气体一元运 动微分方程为：业+血+啰=0 (19)将式(14)代入式(19)，参照文献［1 ］推导可得:p n - kMo2 2Ddv kM(l2 A(1Av n - kMu2 2式中Ma——马赫数马赫数按式(22)计算:Ma =——P式中c 当地声速，m/s管道起始端流速一般较小，马赫数远小于1,因而有:n-K Ma2>0(24)由式(20)和(21)可知，随着流动距离的增长，压力不断减小， 速度不断增加，但压力不可能无限降低，流速不可能无限增大，即 n-K Ma2不可能为0，因此式(25)恒成立：又根据管内稳定流动质量流量方程有:0.25p n d2p v=qm,U(26)将式(14)、(23)、(26)代入式(22)，整理可得:在等内径稳定管道流动中，式(27)中参数除断面平均流速移外 均为常数，因而马赫数随刨增大而增大，图1中点4处马赫数最大有：(28)式中Ma4——点4处的马赫数2 (0.25”少)2伏守(29)将式(4)、(14)、(27)代入式(28)，整理可得:Pi式中p4—一点4处的燃气绝对压力，PaT1——点1处的燃气温度，K根据图1中的。