燃气管道泄漏扩散模型及数值模拟.docx

城镇燃气管道泄漏扩散模型及数值模拟韩克顺 天津城市建设学院 天津摘要：城镇燃气管道的分布区域人口及建筑众多，燃气管道一旦发生泄漏，将有可 能造成重大的财产损失甚至人员伤亡因此 ，为了量化城镇燃气泄露危害，针对管道不同的破坏情况及气源建立了燃气泄漏各种源模型以及扩散模型 ，并且建立了燃气管线动态泄漏扩散模型及伤害性危险范围对第三方破坏所造成的城镇燃气管道泄漏模型进 行了模拟，采用CFD技术对管道泄漏燃气的扩散进行模拟研究 ，获得了泄漏气体的扩散数值模拟结果，为城镇燃气管道安全运行提供了理论依据关键词：燃气管道；泄漏；泄漏模型；扩散模型；数值模拟Abstract : The areas that city gas pipelines distributed in always crowd with people and buildings, and once gas releases through damaged pipes, accidents would happen involving substantial economic losses and even victims amongst the population. Therefore, in order to quantify the hazards of city gas leak, different leakage and diffusion models were established according the characters of gas and the damage. The diffusion model of unsteady leakage of gas pipeline and the damaging and destroying areas were established. Based on the computational fluid dynamics (CFD) technique the diffusion range of leaked gas during accident of underground gas pipeline that caused by the third party damage was studied. According to the results, the dispersion of the leaked gas was obtained and providing the theorybasis for safety operation of city gas pipelines.Key words: Gas pipelines; Leakage; Leakage model; Diffusion model; Numerical simulation由于管道老化、腐蚀、管材和焊缝缺陷等原因 ，尤其是随着市政建设的发展，城镇内违章施工挖断、压裂燃气管道的事故屡见不鲜 ，泄漏事故频发，引发火灾及爆炸事故，造成人员伤亡及环境污染的恶劣后果。

燃气管道泄漏扩散的危险性主要在于两个方面 ：一是一定区域内燃气浓度到达阈限值，使人员中毒；二是一定区域内燃气浓度到达爆炸 极限,范围内发生火灾或爆炸事故针对这种情况 ，本文根据不同条件建立了相应适用的管道泄漏率模型，并应用FLUENTT法进行了城燃管道泄漏模拟1模型分析1.1 燃气泄漏源模型泄漏速率确实定是分析泄漏扩散及评价事故后果的基础 ，其依据的主要模型根据泄漏孔径大小分别为小孔泄漏模型、管道上I型和大孔泄漏模型总的来说 ，腐蚀穿孔容易发生小孔径泄漏事故，而对于外力破坏(诸如施工、自然灾害等)大孔径泄漏事故发生的 可能性较大1.1.1 孔口泄漏 小孔泄漏模型适用于穿孔泄漏的情形 ，穿孔泄漏是指管道或设备由于腐蚀等原因形成小孔，燃气从小孔泄漏对于理想气体等嫡膨胀任意点处的质量流率2M k ./科产+“启Ip J ⑴式中,Q为气体通过孔洞的质量流量，kg/s;自当裂口形状为圆形时取 1.00,三角形时取0.95,长方形时取0.90,孔口为内 层腐蚀形成的渐缩孔（钝角入口）取0.9

从安全工作的角度考虑，都需要通过孔洞流出蒸汽的最大流量引起最大流速的压 力比为:⑵（临界压力），其泄漏量为⑶一岛）一塞压邑是导致孔洞或管道流动最大流量的下游最大压力—qA飞况时1.1.2 管道泄漏 管道泄漏模型适用于开裂泄漏的情形开裂泄漏的原因通常是由于 外力干扰或超压破裂，属于大面积泄漏，泄漏口面积通常为管道截面积的 80%〜100%开裂泄漏瞬时泄漏量大，导致管道中的压力明显降低对于理想气体绝热管道流 ，声速和非声速情况下的质量通量都可以用 Darcy公式表示：m 版HA-P』G= = Oj 一一 （4）式中,G为质量通量；m为气体的质量流率，kg/s;A为孔面积，金1;黑为气体膨胀系数；、为重力常数；p为上游气体密度,kg//；R为上游气体压力，Pa;匚为下游气体压力,Pa;为差损失项，可用4fL/d替换，f为摩擦系数，d为内径,L为管长,m1.2 燃气泄漏及扩散模型城镇燃气一般包括天然气、液化石油气和人工煤气其中天然一及人工煤气密度 】气|,泄漏后为浮性气体被动扩散;液化石油气的气态密度比空气重,泄漏后属于重 气扩散对于不同性质的气体需要应用不同的扩散模型1.2.1 中性浮力扩散 中性浮力扩散模型，用于估算释放发生后释放气体与空气混合 ，并导致混合气云具有中性浮力后下风向各处的浓度 ，对于燃气管道泄漏一般用烟羽模型来描述。

假定风向沿x轴，风速恒定为u, Pasquill-Gifford 模型所描述的位于地面 H高处的连续稳态源的浓度为：〈C〉(x,y,z)= 式中，〈C〉(x,y,z)为平均浓度，g/：J;G为连续泄漏率，kg/s;扩散系数/ % %分别代表下风向、侧风向和垂直方向 (x,y,z)浓度的标准偏差u为风速，m/s;y为风向垂直方向距离，m;z为距里面距离，m;H为泄漏源高度，m下风向地面上最大浓度出现的位置 ，可由下式求得：(6)=H/地面上沿x轴的最大浓度111ax由下式求得:g也」叩1m媪 ⑺1.2.2 重气扩散密度大于其扩散环境空气密度的气体均称为重气 ，现有的重气扩散模型大都是在箱模型(Box-Model)的基础上建立的高斯模型和重气扩散箱式模型有很 大的局限性为了比较精确并贴近实际的模拟重气扩散 ，人们采用计算流体力学(CFD)的方法模拟其扩散的三维非定常态湍流流动过程泄漏燃气视为不可压缩流体 ，其温度视为与环境温度相同，选择k-e双方程模型，则泄漏燃气的运动与扩散规律可由以下控制 方程描述连续性方程：1=0 ⑻动量方程：(9)k- £双方程模型沙洲(10)(11)扩散方程：(12)—+ u— + u— + w— = —fD—J印 瘠 冽 A 第' 索式中各系数的经验数值为：固=0.09, 3=1.44, =1.92, EZj=1.0, 0=1.3采取一定的离散方法将上式离散，并根据具体情况构造边界条件解控制方程 ，便可 得到泄漏场所的燃气浓度分布和速度分布。

在 FLUEN邙软件中建立泄漏场所几何模型 选择动量方程模型、边界条件及其他己知参数，便可自动计算2实例模拟2.1 模型建立本文对燃气泄漏扩散过程进行数值模拟，为了能得到具有指导意义的一般性规律 ，需要从各种实例中抽象出一种有代表意义的模型用三维模型建立燃气管道泄漏这样 的超大空间计算量是巨大的，因此,采用FLUEN啾件建立二维模型模拟分析，模型区域 如图1所示a)模型区域凤必口=>自由出口燃气泄漏口地面(b)计算模型图1泄漏模型2.2 FLUENTT散模拟该管道埋于地下，管道的内径150 7m,燃气及环境温度为300 K, 泄漏口为100 mmx 150 mm,风速为5 m/s,环境压力1 x 1/，燃气泄漏速度为v=330 m/s其中天然气采用甲烷，液化石油气采用丙烷代替2.2.1 天然气泄漏模拟 风速为1 m/s, 3 m/s,10 m/s 时天然气泄漏情况如图 2所示就小知聃皿川1 7 *taRl7 CKta-01 6506-01EgWSgOi 4 gqi4 gMi 3他小 300B-CI12 50=(11 20Cfe-OI 1皿Ml1幅小 务值版 OWl+WCchiIm「窜 erf Fraetjon erf ch4Apr 19 2012FLUENT 6.3 ：2

可以看出，随着泄漏的进行，相对较高的浓度区扩大速度变的 缓慢2.2.1 液化石油气泄漏模拟 风速为1 m/s, 3 m/s,10 m/s 时液化石油气泄漏情况如图3所示9班小 4001411整的001/ SOB 闻 I 7OCta-0le5Qg-dl 触 CJQB4J15珈m 专 gQ0l44 OOUll ^50h-OI25CS-Q1 之钝好＜11 I灿加 loceci^009+00CmIdu巾 of %比bon of c3nAApr 19 曾 12FLUENT 6.3 2d dp邛 spe. ske)(a)风速为1 m/sCmlGXiri at 牖的生 faction of cM.Apr 19 2012FLUENT 6,3 ?2(js dp. pbns,3麻)1TO+00融孙gi 8g9/1 TWb-OI 74411 B4OI bCJQBTJI 1^™写如m专g0l 4441 4O0BO1 3独31 34川I 25Ce-<]1 2g41 1处01 10dl &gg QW-HKI〔c〕风速为10 m/s图3液化石油气扩散浓度分布云图〔b〕风速为3 m/s由图2与图3相比可知，风速对燃气的泄漏扩散情况占主导地位。

风速较大时 体会向下风向扩散，并且在竖直方向的扩散会受到抑制。