第28卷 第3期2008年3月
煤气与热力
G AS &HE AT
Vol . 28No . 3Mar . 2008
燃气管道泄漏流量的计算
黄小美, 彭世尼, 徐海东, 杨茂华
1
1
2
3
(1. 重庆大学城市建设与环境工程学院, 重庆400045; 2. ,
山东潍坊261061; 3. )
摘 要: , 联合两种模型计算任何泄漏孔口直径, 进行了实例计算并对比了不同模型的计算结果。 关键词:泄漏流量计算; 管道模型; 小孔模型; 管道小孔综合模型; 流量限制
中图分类号:T U996 文献标识码:A 文章编号:1000-4416(2008) 03-0B11-06
Ca lcul a ti on of L eakage Ra te fro m Ga s P i peli n e
HUANG Xiao 2mei, PENG Shi 2ni, XU Hai 2dong, Y ANG M ao 2hua
Abstract: According t o the basic equati ons of one 2di m ensi on gas fl ow, a hole model f or calcula 2ti on of hole leakage rate in adiabatic p r ocess and a p i peline model for calculati on of leakage rate in varia 2ble p r ocess suited t o full rup ture of p i peline are deducted . These t w o kinds of models are combined t o calculate the leakage rate fr om leakage holes with different dia meters . The li m itati on of the maxi m um gas leakage rate is discussed, the exa mp le calculati on is carried out, and the calculati on results of different models are compared .
Key words: calculati on of leakage rate; p i peline model; hole model; combined model of p i peline model and hole model; li m itati on of fl o w rate
概述1
在燃气管道事故定量风险评价、事故抢险预案制定和漏气损失评估时, 首先要计算泄漏流量。燃气管道在事故破损时, 燃气可通过两种途径进入到大气中, 一种是燃气直接泄漏到大气环境中, 另一种是泄漏到土壤中, 通过土壤渗透进入大气环境。前者可以通过理论推导得出泄漏流量的计算公式, 后者理论计算比较复杂且不确定性很大。本文主要分析和讨论前一种情况下的泄漏流量计算。第三方破坏是城市燃气管道泄漏的主要原因之一, 其主要表现是挖掘机器、钻孔机器破坏管道, 在这种情况下, 燃气通常直接泄漏到大气中。此外, 架空管道泄漏也是直接泄漏到大气中。
口当作孔径很小的小孔, 从而建立泄漏流量计算的模型
。
图1 管道泄漏示意图
图中 点1———管道起点
点2———泄漏口入口点点3———泄漏口出口截面上的点点4———点2上游附近的某点L ———泄漏点至管道起点的距离, m
3
q V , U ———泄漏点上游管道体积流量, m /h
q V ———泄漏体积流量, m /h
3
小孔模型的推导2
管道泄漏示意图见图1。小孔模型是将泄漏孔
图1中, 点1通常为该管道上游的调压器出口,
第28卷 第3期
煤气与热力
www . watergasheat . com
2
2
κκ
其压力通常保持不变。假设点4的断面流量及其平均流速方向不受泄漏影响, 而点4下游至泄漏口处的任何点管道断面平均流速由于受到泄漏影响而不再沿管道轴线方向, 点4至点2的距离非常小, 可以忽略不计, 因而点4的压力近似等于点2的压力。
小孔模型假设管内燃气全部从该小孔泄漏, 即管道上游无支管或支管燃气流量为0, 这样假设是为了保证从小孔泄漏的燃气流量是最大值; 由于泄漏小孔孔径较小, 泄漏流量有限, 因而忽略管道沿程阻力, , 即:
p 2=p 1
式中 p 2———图12的绝对压力, Pa
p 1———图1中点1的绝对压力, Pa
v 3
2
-
v 2, L
2
κ
R T 1-κ-1con 2
p 3p 2
(5)
式中 v 3———点3的燃气断面平均流速, 即为燃气
的泄漏出口流速, m /s
v 2, L ———点2沿管道泄漏口轴线方向的流速,
m /s
T 2———点2, K p 3点33(62]
:
(6) (7)
3(1)
(p a >p c ) p c (p a ≤p c )
p p 2
=β
κ
κ-1
在泄漏孔处, 燃气流速一般较快, 燃气没有足够的时间与环境进行热量交换, 因此燃气泄漏过程, 即从点2到点3的燃气流动过程可被视为可压缩气体绝热流动过程, 可见泄漏孔口与喷嘴相似。孔口泄漏瞬间的流动可以看作是一维流动, 气体的一元流动欧拉运动微分方程为+v d v =0
ρ
[1]
β=
κ+(8)
式中 p a ———环境压力, Pa
p c ———临界压力, Pa
β———临界压力比
在小孔模型下, v 2, L 一般很小或近似为0, 因而
v 2, L 远远小于v 3, 可以忽略不计, 因此式(5) 可简化
2
2
:
(2)
为: v 3=
R T 1--1con 2
p 3p 2
κ
式中 p ———燃气绝对压力, Pa
3
ρ———燃气密度, kg/m
v ———燃气断面平均流速, m /s
(9)
泄漏孔出口压力与入口压力的比值等于临界压
[2]
因为泄漏过程为绝热过程, 所以有κ=C 1ρ
κ—式中 ——燃气的等熵指数
C 1———常量:
(3)
力比时的泄漏出口燃气流速称为临界流速, 将式
(7) 、(8) 代入式(9) 有: v c =
R T -1con 2
(10)
式中 v c ———泄漏出口燃气临界流速, m /s
当p 3达到临界压力p c 后, 点3的流速v 3就等于临界流速v c 。
根据孔口出流的质量流量公式, 可得出泄漏质量流量为:
πd ρ q m =0. 25μ3v 3
式中 q m ———泄漏质量流量, kg/s
μ———流量系数, 可取0. 90~0. 98
d ———泄漏孔口当量直径, m
3ρ——点3处的燃气密度, kg/m3—
2
κ是温度的函数, 在常温下理想气体的κ可近似当作定值, 对于天然气等由多原子分子组成的
气体, κ取1. 29。
因为城市燃气压力大多在1. 6MPa 以下, 燃气温度为常温, 所以燃气可以看作理想气体有:
=R con T
ρ
式中 R con ———燃气的气体常数, J /(kg ・K )
T ———燃气的温度, K [2]
[2]
, 因此
(11)
(4)
不规则孔口当量直径按下式计算: d =
A j
将式(3) 代入式(2) 并在小孔入口点2和小孔出口点3进行积分, 然后将式(4) 代入整理可得:
www . watergasheat . com
2
黄小美, 等:燃气管道泄漏流量的计算 第28卷 第3期
式中 A j ———泄漏孔口面积, m
(6) 、(9) 、(10) 代入式(11) 并整理可将式(4) 、
得:
当
p a
>β时:p 2
μπd 2・ q m =0. 25当
p a p 2
p 2con T 2
-p p κ
-
p p κ
(12)
β时:≤
2
q m =0. con ++κ-(13)
在小孔模型下, 认为式(1) 成立, 且由于管内流速小, 管内流动可视为等温过程, 因而对于任一管道, 只要知道管道起点压力、管内燃气温度和泄漏孔口当量直径, 便可按照式(12) 和(13) 计算泄漏质量流量。
式中 q m , U ———泄漏点上游管道内燃气的质量流
量, kg/s
D ———管道内径, m λ———摩擦阻力系数, 可按文献[4]或[5]提
供的公式计算
式(15) 中n =1时, 表示管内流动为等温过程, , 管; 当n =κ时, 表示, ; 实, n 的值在1和κ之间, 为了简化计算, 通常在管内燃气流速较小, 管道较长时, n 取1, 而在管道流速很大或者管道很短时, n 取κ。在管道模型下, 点2为管道末端且暴露于大气环境中, 因而泄漏口处的压力等于大气环境压力, 即:
(16) p 2=p a
此时只要知道管径、泄漏处至管道起点的距离、
管道起点的压力, 选取适当的值, 便可按照式(15) 计算管道的质量流量q m , U 。
管道模型的推导3
管道模型适合于管道完全断裂的情形, 即泄漏当量直径等于管道内径, 点2即为管道末端, 点2和点3重合, 管道泄漏流量就等于管输流量, 此时可按管道水力计算公式来计算管道流量。
燃气泄漏过程中, 管道上游阀门关闭之前, 管内燃气流动可以视为稳定流动。在燃气管道完全断裂泄漏情况下, 燃气流速较大, 管内燃气可能没有足够的时间和周边环境进行充分的热交换, 管内流动也不能看作等温过程, 而只能看作多变过程, 因而对于管内燃气流动过程有:n =C 2
ρ
[3]
管道-小孔综合模型4
当燃气管道完全断裂时, 按管道模型计算泄漏流量比较准确, 当燃气管道只有孔径很小的破损孔时, 按小孔模型计算泄漏流量比较准确。但在实际情况下, 特别是由于施工开挖导致的断裂, 泄漏口既不是小孔, 也不是完全断裂, 因此用这两种模型都不准确。管道-小孔综合模型就是将管道模型和小孔
[6]
模型结合起来进行泄漏流量计算。
事实上, 在小孔模型下, 只要有燃气泄漏, 管道内燃气就会流动, 由于摩擦的存在, 管道内必然会有沿程阻力, 因此式(1) 并不成立, 实际的p 2还需要根据式(15) 计算。前面已经假设管输燃气全部从小孔泄漏, 这种假设对于风险评价和事故应急抢险来说是保守的, 因而有:
(17) q m =q m , U
式(17) 中q m 可根据式(12) 和(13) 计算, 式(12) 和(13) 中的T 2可由式(4) 和(14) 联立求得,
(12) 、q m , U 可以根据式(15) 计算, 因此联立式(4) 、
(13) 、(14) 、(15) 、(17) 便可计算出泄漏质量流量q m 、泄漏处的管内压力p 2及管内温度T 2。该模型对于任何泄漏口尺寸的稳定泄漏过程都比较准确, 计
(14)
式中 n ———多变指数
C 2———常量
仍然假设从管道起点至泄漏点之间的管段上没有支管, 或所有支管的流量都为0, 同时考虑到管内
[4]
燃气流速较大时, 不能忽略对流项, 将文献[4]提供的管内流动基本微分方程组联合式(14) 可以推导得到式(15) 。
2
16(n +1) q m , U n n n
p 1-p 2p 1・
π2D 4n R con T 1
p 1λL ln 2D n p (15)
第28卷 第3期
煤气与热力
www . watergasheat . com
算时需要利用计算机语言编程求解。
为了符合工程习惯, 通常需要将泄漏质量流量转换成泄漏体积流量: q V =3600
q m
M a
(25)
ρ0
3
(18)
又根据管内稳定流动质量流量方程有:
πD 2ρ(26) 0. 25v =q m , U
(23) 、(26) 代入式(22) , 整理可得:将式(14) 、 M a =
2ρ式中 ——标准状态下燃气密度, kg/m0—
2
最大泄漏流量限制的讨论5
燃气在管道内的最大流量有限制, 并不是按照
管道模型或者管道-小孔模型计算得到的结果就一定是实际流量, 校验:界流量的限制。
① 城市燃气管道都与上游调压器相连, 管道内的流量和泄漏流量都不能大于上游调压器的最大通过能力q V , max , 调压器说明书通常都给出了该调压器在各进口压力下的q V , max 值, 因此无论用哪种模型计算, 泄漏流量的计算结果最大只能取q V , max 。
② 管内燃气临界流量的限制
与泄漏孔处的绝热流动相似, 燃气管内流动也存在临界流量问题。由于沿程摩擦阻力对燃气运动参数的影响, 实际气体一元运动微分方程为:
λv 2(19) +v d v =0
p 2D
将式(14) 代入式(19) , 参照文献[1]推导可得:
2
λ(20) -=2
p n -κ2D M a
2
κλd L M a
(21) 2v n -κ2D M a
式中 M a ———马赫数
马赫数按式(22) 计算: M a c
C 2q n -1
m , U
v 2
(27)
, (27) 中参数除断, v 增大而, 4。有: M a 4
(28)
式中 M a 4———点4处的马赫数
(14) 、(27) 代入式(28) , 整理可得:将式(4) 、 q m , U
2
(0. 25πD 2) 2p 2p 1
R con T 1
p n
(29)
式中 p 4———点4处的燃气绝对压力, Pa
T 1———点1处的燃气温度, K
根据图1中的假设有:
(30) p 4=p 2
将式(30) 代入式(29) 并转化成体积流量, 则有:
q V , U
2n
(31)
令管内临界体积流量q V , c, P 为: q V , c, P =900πD 2p 1
p 2n
p con T 10
3
(32)
(22) (23)
式中 q V , c, P ———管内燃气临界体积流量, m /h
无论按照何种模型计算所得的体积流量必须校验是否满足式(31) , 若计算得到的泄漏体积流量大于q V , c, P , 则只能取q V , c, P 。
6 实例计算
c =ρ
式中 c ———当地声速, m /s
管道起始端流速一般较小, 马赫数远小于1, 因而有:
2
(24) n -κM a >0
由式(20) 和(21) 可知, 随着流动距离的增长, 压力不断减小, 速度不断增加, 但压力不可能无限降
2
低, 流速不可能无限增大, 即n -κM a 不可能为0, 因此式(25) 恒成立:
某管道受施工破坏而连续泄漏, 泄漏口近似为圆形。该管道上游调压器出口绝对压力为0. 5MPa, 调压器最大通过能力为35000m /h, 泄漏点
3
距离调压器1k m , 管道内径为200mm 。已知天然
3
气密度为0. 76kg/m, 管道周边环境温度为288K 。分析在不同泄漏当量直径下的泄漏流量。
① 最大流量限制在管道完全断裂时泄漏流量最大, 此时p 2=101325Pa 。由于燃气流速较大, 假设管内流动为绝
www . watergasheat . com 黄小美, 等:燃气管道泄漏流量的计算 第28卷 第3期
热过程, n 取1. 29, 将p 2代入式(32) 计算得到管道内燃气临界体积流量为:
3
q V , c, P =53184m /h
而该管道上游调压器最大通过能力为35000m /h, 因此最大泄漏体积流量限制为35000m /h 。
3
3
② 管道模型计算
(14) 、(15) 、管道完全断裂情况下, 根据式(4) 、
(16) 编程计算, 可得燃气泄漏的体积流量为213143
m /h, 小于最大泄漏体积流量限制值, 因此该值为管道在此位置的最大可能泄漏体积流量③ 管道-采用管道, 4) 、(12) 、(13) 、(14) ) 、() 编程计算, 其中n 按式(33) 确定:
n =1+(κ-1) 2
D
2
(33)
式(33) 假设多变指数与泄漏孔口面积成线性关系, 当管道完全断裂时, 管内流动为绝热过程, 而当泄漏孔口非常小的时候, 为等温过程。
经计算, 泄漏孔口当量直径从1mm 增加到200mm , 管道-小孔综合模型泄漏体积流量q V 、泄漏处的管内压力p 2、泄漏处的管内温度T 2随泄漏孔口当量直径的变化规律分别见图2、3、4中的曲线1, 小孔模型计算结果则分别见图2、3、4中的曲线2
。
算结果比后者低, 曲线1前段为直线是由于考虑了
最大流量的限制而形成的
。
图5 泄漏点位置与该点最大泄漏体积流量的关系
结论及建议7
图2 泄漏体积流量与泄漏孔口当量直径的关系
④ 不同泄漏位置最大泄漏体积流量计算本文推导的管道模型与文献[4]或[5]中提供的管道水力计算公式相比, 考虑了管内流动的温度变化和对流项的影响。为了比较这两个公式计算结果的差别, 在本例中假设泄漏当量直径等于管道内径, 分别用两个公式计算管道不同泄漏位置的泄漏体积流量, 结果对比见图5。当泄漏处至管道起点
较远时, 两者计算结果基本一致(前者计算结果比后者略高) , 而当泄漏处距管道起点较近时,
前者计
燃气管道泄漏流量计算是燃气管道风险评价的重要环节, 也是事故应急抢险的重要依据, 利用合理假设推导的管道-小孔综合模型计算泄漏流量, 同时考虑泄漏流量的限制条件, 能够得到比较精确的结果, 能够满足风险评价和应急抢险的需要。
本文的模型没有考虑非稳定过程, 实际工程中, 在管道开始泄漏瞬间和当泄漏管道的上游阀门被关闭后, 燃气流动是非稳定流动, 但在燃气持续泄漏期间, 燃气流动可看作稳定流动, 在风险评价和事故应急抢险过程中, 往往关心的是这种持续泄漏的影响,
第28卷 第3期
煤气与热力
www . watergasheat . com
因此可不必进行复杂的非稳定流动计算。若需要更为精确的计算结果, 则需要建立管内燃气与周边环境之间的传热关系, 建立非稳定流动和传热方程组, 确定边界条件, 采用计算流体力学方法计算。参考文献:
[1] 周谟仁. 流体力学 泵与风机(第3版) [M].北京:
[4] 段常贵. 燃气输配(第3版) [M].北京:中国建筑工
业出版社, 2001.
[5] G B 50028—2006, 城镇燃气设计规范[S].
[6] Montiel H, V ilchez A, Casal J, et al . Mathe matical model 2
ling of accidental gas releases[J ].Journal of Hazardous Materials, 1998, (59) :211-233.
中国建筑工业出版社, 1994.
[2] 邱信立, 廉乐明, 李力能. 工程热力学(第3版) [M].
作者简介:黄小美(1980-) , 湖南安仁人,
、节能及。
:E -ma il:hx m1980@126. com
北京:中国建筑工业出版社, 1992.
[3] Dong Y, Gao H, Zhou J, et al . of
gas release m gineering (:-241.
收稿日期:2007-01-12; 修回日期:2008-01-09
・标准规范简讯・
(G B 50028—2006) 勘误表《城镇燃气设计规范》
城镇燃气设计规范国家标准管理组将“《城镇燃气设计规范》G B 50028—2006勘误表”作为管理组“国燃规标字[2007]第31号函”的附件上报给建设部标准定额司。2007年9月7日, 建设部标准定额司以“建标标便[2007]28号函”批准了此勘误表。待中国建筑工业出版社重印《城镇燃气设计规范》时, 将按照此勘误表对相关内容进行更正。
(G B 50028—2006) 勘误表《城镇燃气设计规范》
序号
1234567
页号-行数
79-3143-倒10149-9158-8、291-倒11302-14
误
用阴极保护。
2……()
正
用阴极保护。
2……()
……现场人工开启型, , , ……现场人工开启型。, 城市
气管道, 在避雷保护范围以外时, ……停电时紧急自动切断阀必须处于关闭状态;
137-倒7、137-倒6
建筑沿外墙架设的燃气管道, ……
8
150-12、(150-13
) ;
(城镇燃气设计规范国家标准管理组 供稿)
第28卷 第3期2008年3月
煤气与热力
G AS &HE AT
Vol . 28No . 3Mar . 2008
燃气管道泄漏流量的计算
黄小美, 彭世尼, 徐海东, 杨茂华
1
1
2
3
(1. 重庆大学城市建设与环境工程学院, 重庆400045; 2. ,
山东潍坊261061; 3. )
摘 要: , 联合两种模型计算任何泄漏孔口直径, 进行了实例计算并对比了不同模型的计算结果。 关键词:泄漏流量计算; 管道模型; 小孔模型; 管道小孔综合模型; 流量限制
中图分类号:T U996 文献标识码:A 文章编号:1000-4416(2008) 03-0B11-06
Ca lcul a ti on of L eakage Ra te fro m Ga s P i peli n e
HUANG Xiao 2mei, PENG Shi 2ni, XU Hai 2dong, Y ANG M ao 2hua
Abstract: According t o the basic equati ons of one 2di m ensi on gas fl ow, a hole model f or calcula 2ti on of hole leakage rate in adiabatic p r ocess and a p i peline model for calculati on of leakage rate in varia 2ble p r ocess suited t o full rup ture of p i peline are deducted . These t w o kinds of models are combined t o calculate the leakage rate fr om leakage holes with different dia meters . The li m itati on of the maxi m um gas leakage rate is discussed, the exa mp le calculati on is carried out, and the calculati on results of different models are compared .
Key words: calculati on of leakage rate; p i peline model; hole model; combined model of p i peline model and hole model; li m itati on of fl o w rate
概述1
在燃气管道事故定量风险评价、事故抢险预案制定和漏气损失评估时, 首先要计算泄漏流量。燃气管道在事故破损时, 燃气可通过两种途径进入到大气中, 一种是燃气直接泄漏到大气环境中, 另一种是泄漏到土壤中, 通过土壤渗透进入大气环境。前者可以通过理论推导得出泄漏流量的计算公式, 后者理论计算比较复杂且不确定性很大。本文主要分析和讨论前一种情况下的泄漏流量计算。第三方破坏是城市燃气管道泄漏的主要原因之一, 其主要表现是挖掘机器、钻孔机器破坏管道, 在这种情况下, 燃气通常直接泄漏到大气中。此外, 架空管道泄漏也是直接泄漏到大气中。
口当作孔径很小的小孔, 从而建立泄漏流量计算的模型
。
图1 管道泄漏示意图
图中 点1———管道起点
点2———泄漏口入口点点3———泄漏口出口截面上的点点4———点2上游附近的某点L ———泄漏点至管道起点的距离, m
3
q V , U ———泄漏点上游管道体积流量, m /h
q V ———泄漏体积流量, m /h
3
小孔模型的推导2
管道泄漏示意图见图1。小孔模型是将泄漏孔
图1中, 点1通常为该管道上游的调压器出口,
第28卷 第3期
煤气与热力
www . watergasheat . com
2
2
κκ
其压力通常保持不变。假设点4的断面流量及其平均流速方向不受泄漏影响, 而点4下游至泄漏口处的任何点管道断面平均流速由于受到泄漏影响而不再沿管道轴线方向, 点4至点2的距离非常小, 可以忽略不计, 因而点4的压力近似等于点2的压力。
小孔模型假设管内燃气全部从该小孔泄漏, 即管道上游无支管或支管燃气流量为0, 这样假设是为了保证从小孔泄漏的燃气流量是最大值; 由于泄漏小孔孔径较小, 泄漏流量有限, 因而忽略管道沿程阻力, , 即:
p 2=p 1
式中 p 2———图12的绝对压力, Pa
p 1———图1中点1的绝对压力, Pa
v 3
2
-
v 2, L
2
κ
R T 1-κ-1con 2
p 3p 2
(5)
式中 v 3———点3的燃气断面平均流速, 即为燃气
的泄漏出口流速, m /s
v 2, L ———点2沿管道泄漏口轴线方向的流速,
m /s
T 2———点2, K p 3点33(62]
:
(6) (7)
3(1)
(p a >p c ) p c (p a ≤p c )
p p 2
=β
κ
κ-1
在泄漏孔处, 燃气流速一般较快, 燃气没有足够的时间与环境进行热量交换, 因此燃气泄漏过程, 即从点2到点3的燃气流动过程可被视为可压缩气体绝热流动过程, 可见泄漏孔口与喷嘴相似。孔口泄漏瞬间的流动可以看作是一维流动, 气体的一元流动欧拉运动微分方程为+v d v =0
ρ
[1]
β=
κ+(8)
式中 p a ———环境压力, Pa
p c ———临界压力, Pa
β———临界压力比
在小孔模型下, v 2, L 一般很小或近似为0, 因而
v 2, L 远远小于v 3, 可以忽略不计, 因此式(5) 可简化
2
2
:
(2)
为: v 3=
R T 1--1con 2
p 3p 2
κ
式中 p ———燃气绝对压力, Pa
3
ρ———燃气密度, kg/m
v ———燃气断面平均流速, m /s
(9)
泄漏孔出口压力与入口压力的比值等于临界压
[2]
因为泄漏过程为绝热过程, 所以有κ=C 1ρ
κ—式中 ——燃气的等熵指数
C 1———常量:
(3)
力比时的泄漏出口燃气流速称为临界流速, 将式
(7) 、(8) 代入式(9) 有: v c =
R T -1con 2
(10)
式中 v c ———泄漏出口燃气临界流速, m /s
当p 3达到临界压力p c 后, 点3的流速v 3就等于临界流速v c 。
根据孔口出流的质量流量公式, 可得出泄漏质量流量为:
πd ρ q m =0. 25μ3v 3
式中 q m ———泄漏质量流量, kg/s
μ———流量系数, 可取0. 90~0. 98
d ———泄漏孔口当量直径, m
3ρ——点3处的燃气密度, kg/m3—
2
κ是温度的函数, 在常温下理想气体的κ可近似当作定值, 对于天然气等由多原子分子组成的
气体, κ取1. 29。
因为城市燃气压力大多在1. 6MPa 以下, 燃气温度为常温, 所以燃气可以看作理想气体有:
=R con T
ρ
式中 R con ———燃气的气体常数, J /(kg ・K )
T ———燃气的温度, K [2]
[2]
, 因此
(11)
(4)
不规则孔口当量直径按下式计算: d =
A j
将式(3) 代入式(2) 并在小孔入口点2和小孔出口点3进行积分, 然后将式(4) 代入整理可得:
www . watergasheat . com
2
黄小美, 等:燃气管道泄漏流量的计算 第28卷 第3期
式中 A j ———泄漏孔口面积, m
(6) 、(9) 、(10) 代入式(11) 并整理可将式(4) 、
得:
当
p a
>β时:p 2
μπd 2・ q m =0. 25当
p a p 2
p 2con T 2
-p p κ
-
p p κ
(12)
β时:≤
2
q m =0. con ++κ-(13)
在小孔模型下, 认为式(1) 成立, 且由于管内流速小, 管内流动可视为等温过程, 因而对于任一管道, 只要知道管道起点压力、管内燃气温度和泄漏孔口当量直径, 便可按照式(12) 和(13) 计算泄漏质量流量。
式中 q m , U ———泄漏点上游管道内燃气的质量流
量, kg/s
D ———管道内径, m λ———摩擦阻力系数, 可按文献[4]或[5]提
供的公式计算
式(15) 中n =1时, 表示管内流动为等温过程, , 管; 当n =κ时, 表示, ; 实, n 的值在1和κ之间, 为了简化计算, 通常在管内燃气流速较小, 管道较长时, n 取1, 而在管道流速很大或者管道很短时, n 取κ。在管道模型下, 点2为管道末端且暴露于大气环境中, 因而泄漏口处的压力等于大气环境压力, 即:
(16) p 2=p a
此时只要知道管径、泄漏处至管道起点的距离、
管道起点的压力, 选取适当的值, 便可按照式(15) 计算管道的质量流量q m , U 。
管道模型的推导3
管道模型适合于管道完全断裂的情形, 即泄漏当量直径等于管道内径, 点2即为管道末端, 点2和点3重合, 管道泄漏流量就等于管输流量, 此时可按管道水力计算公式来计算管道流量。
燃气泄漏过程中, 管道上游阀门关闭之前, 管内燃气流动可以视为稳定流动。在燃气管道完全断裂泄漏情况下, 燃气流速较大, 管内燃气可能没有足够的时间和周边环境进行充分的热交换, 管内流动也不能看作等温过程, 而只能看作多变过程, 因而对于管内燃气流动过程有:n =C 2
ρ
[3]
管道-小孔综合模型4
当燃气管道完全断裂时, 按管道模型计算泄漏流量比较准确, 当燃气管道只有孔径很小的破损孔时, 按小孔模型计算泄漏流量比较准确。但在实际情况下, 特别是由于施工开挖导致的断裂, 泄漏口既不是小孔, 也不是完全断裂, 因此用这两种模型都不准确。管道-小孔综合模型就是将管道模型和小孔
[6]
模型结合起来进行泄漏流量计算。
事实上, 在小孔模型下, 只要有燃气泄漏, 管道内燃气就会流动, 由于摩擦的存在, 管道内必然会有沿程阻力, 因此式(1) 并不成立, 实际的p 2还需要根据式(15) 计算。前面已经假设管输燃气全部从小孔泄漏, 这种假设对于风险评价和事故应急抢险来说是保守的, 因而有:
(17) q m =q m , U
式(17) 中q m 可根据式(12) 和(13) 计算, 式(12) 和(13) 中的T 2可由式(4) 和(14) 联立求得,
(12) 、q m , U 可以根据式(15) 计算, 因此联立式(4) 、
(13) 、(14) 、(15) 、(17) 便可计算出泄漏质量流量q m 、泄漏处的管内压力p 2及管内温度T 2。该模型对于任何泄漏口尺寸的稳定泄漏过程都比较准确, 计
(14)
式中 n ———多变指数
C 2———常量
仍然假设从管道起点至泄漏点之间的管段上没有支管, 或所有支管的流量都为0, 同时考虑到管内
[4]
燃气流速较大时, 不能忽略对流项, 将文献[4]提供的管内流动基本微分方程组联合式(14) 可以推导得到式(15) 。
2
16(n +1) q m , U n n n
p 1-p 2p 1・
π2D 4n R con T 1
p 1λL ln 2D n p (15)
第28卷 第3期
煤气与热力
www . watergasheat . com
算时需要利用计算机语言编程求解。
为了符合工程习惯, 通常需要将泄漏质量流量转换成泄漏体积流量: q V =3600
q m
M a
(25)
ρ0
3
(18)
又根据管内稳定流动质量流量方程有:
πD 2ρ(26) 0. 25v =q m , U
(23) 、(26) 代入式(22) , 整理可得:将式(14) 、 M a =
2ρ式中 ——标准状态下燃气密度, kg/m0—
2
最大泄漏流量限制的讨论5
燃气在管道内的最大流量有限制, 并不是按照
管道模型或者管道-小孔模型计算得到的结果就一定是实际流量, 校验:界流量的限制。
① 城市燃气管道都与上游调压器相连, 管道内的流量和泄漏流量都不能大于上游调压器的最大通过能力q V , max , 调压器说明书通常都给出了该调压器在各进口压力下的q V , max 值, 因此无论用哪种模型计算, 泄漏流量的计算结果最大只能取q V , max 。
② 管内燃气临界流量的限制
与泄漏孔处的绝热流动相似, 燃气管内流动也存在临界流量问题。由于沿程摩擦阻力对燃气运动参数的影响, 实际气体一元运动微分方程为:
λv 2(19) +v d v =0
p 2D
将式(14) 代入式(19) , 参照文献[1]推导可得:
2
λ(20) -=2
p n -κ2D M a
2
κλd L M a
(21) 2v n -κ2D M a
式中 M a ———马赫数
马赫数按式(22) 计算: M a c
C 2q n -1
m , U
v 2
(27)
, (27) 中参数除断, v 增大而, 4。有: M a 4
(28)
式中 M a 4———点4处的马赫数
(14) 、(27) 代入式(28) , 整理可得:将式(4) 、 q m , U
2
(0. 25πD 2) 2p 2p 1
R con T 1
p n
(29)
式中 p 4———点4处的燃气绝对压力, Pa
T 1———点1处的燃气温度, K
根据图1中的假设有:
(30) p 4=p 2
将式(30) 代入式(29) 并转化成体积流量, 则有:
q V , U
2n
(31)
令管内临界体积流量q V , c, P 为: q V , c, P =900πD 2p 1
p 2n
p con T 10
3
(32)
(22) (23)
式中 q V , c, P ———管内燃气临界体积流量, m /h
无论按照何种模型计算所得的体积流量必须校验是否满足式(31) , 若计算得到的泄漏体积流量大于q V , c, P , 则只能取q V , c, P 。
6 实例计算
c =ρ
式中 c ———当地声速, m /s
管道起始端流速一般较小, 马赫数远小于1, 因而有:
2
(24) n -κM a >0
由式(20) 和(21) 可知, 随着流动距离的增长, 压力不断减小, 速度不断增加, 但压力不可能无限降
2
低, 流速不可能无限增大, 即n -κM a 不可能为0, 因此式(25) 恒成立:
某管道受施工破坏而连续泄漏, 泄漏口近似为圆形。该管道上游调压器出口绝对压力为0. 5MPa, 调压器最大通过能力为35000m /h, 泄漏点
3
距离调压器1k m , 管道内径为200mm 。已知天然
3
气密度为0. 76kg/m, 管道周边环境温度为288K 。分析在不同泄漏当量直径下的泄漏流量。
① 最大流量限制在管道完全断裂时泄漏流量最大, 此时p 2=101325Pa 。由于燃气流速较大, 假设管内流动为绝
www . watergasheat . com 黄小美, 等:燃气管道泄漏流量的计算 第28卷 第3期
热过程, n 取1. 29, 将p 2代入式(32) 计算得到管道内燃气临界体积流量为:
3
q V , c, P =53184m /h
而该管道上游调压器最大通过能力为35000m /h, 因此最大泄漏体积流量限制为35000m /h 。
3
3
② 管道模型计算
(14) 、(15) 、管道完全断裂情况下, 根据式(4) 、
(16) 编程计算, 可得燃气泄漏的体积流量为213143
m /h, 小于最大泄漏体积流量限制值, 因此该值为管道在此位置的最大可能泄漏体积流量③ 管道-采用管道, 4) 、(12) 、(13) 、(14) ) 、() 编程计算, 其中n 按式(33) 确定:
n =1+(κ-1) 2
D
2
(33)
式(33) 假设多变指数与泄漏孔口面积成线性关系, 当管道完全断裂时, 管内流动为绝热过程, 而当泄漏孔口非常小的时候, 为等温过程。
经计算, 泄漏孔口当量直径从1mm 增加到200mm , 管道-小孔综合模型泄漏体积流量q V 、泄漏处的管内压力p 2、泄漏处的管内温度T 2随泄漏孔口当量直径的变化规律分别见图2、3、4中的曲线1, 小孔模型计算结果则分别见图2、3、4中的曲线2
。
算结果比后者低, 曲线1前段为直线是由于考虑了
最大流量的限制而形成的
。
图5 泄漏点位置与该点最大泄漏体积流量的关系
结论及建议7
图2 泄漏体积流量与泄漏孔口当量直径的关系
④ 不同泄漏位置最大泄漏体积流量计算本文推导的管道模型与文献[4]或[5]中提供的管道水力计算公式相比, 考虑了管内流动的温度变化和对流项的影响。为了比较这两个公式计算结果的差别, 在本例中假设泄漏当量直径等于管道内径, 分别用两个公式计算管道不同泄漏位置的泄漏体积流量, 结果对比见图5。当泄漏处至管道起点
较远时, 两者计算结果基本一致(前者计算结果比后者略高) , 而当泄漏处距管道起点较近时,
前者计
燃气管道泄漏流量计算是燃气管道风险评价的重要环节, 也是事故应急抢险的重要依据, 利用合理假设推导的管道-小孔综合模型计算泄漏流量, 同时考虑泄漏流量的限制条件, 能够得到比较精确的结果, 能够满足风险评价和应急抢险的需要。
本文的模型没有考虑非稳定过程, 实际工程中, 在管道开始泄漏瞬间和当泄漏管道的上游阀门被关闭后, 燃气流动是非稳定流动, 但在燃气持续泄漏期间, 燃气流动可看作稳定流动, 在风险评价和事故应急抢险过程中, 往往关心的是这种持续泄漏的影响,
第28卷 第3期
煤气与热力
www . watergasheat . com
因此可不必进行复杂的非稳定流动计算。若需要更为精确的计算结果, 则需要建立管内燃气与周边环境之间的传热关系, 建立非稳定流动和传热方程组, 确定边界条件, 采用计算流体力学方法计算。参考文献:
[1] 周谟仁. 流体力学 泵与风机(第3版) [M].北京:
[4] 段常贵. 燃气输配(第3版) [M].北京:中国建筑工
业出版社, 2001.
[5] G B 50028—2006, 城镇燃气设计规范[S].
[6] Montiel H, V ilchez A, Casal J, et al . Mathe matical model 2
ling of accidental gas releases[J ].Journal of Hazardous Materials, 1998, (59) :211-233.
中国建筑工业出版社, 1994.
[2] 邱信立, 廉乐明, 李力能. 工程热力学(第3版) [M].
作者简介:黄小美(1980-) , 湖南安仁人,
、节能及。
:E -ma il:hx m1980@126. com
北京:中国建筑工业出版社, 1992.
[3] Dong Y, Gao H, Zhou J, et al . of
gas release m gineering (:-241.
收稿日期:2007-01-12; 修回日期:2008-01-09
・标准规范简讯・
(G B 50028—2006) 勘误表《城镇燃气设计规范》
城镇燃气设计规范国家标准管理组将“《城镇燃气设计规范》G B 50028—2006勘误表”作为管理组“国燃规标字[2007]第31号函”的附件上报给建设部标准定额司。2007年9月7日, 建设部标准定额司以“建标标便[2007]28号函”批准了此勘误表。待中国建筑工业出版社重印《城镇燃气设计规范》时, 将按照此勘误表对相关内容进行更正。
(G B 50028—2006) 勘误表《城镇燃气设计规范》
序号
1234567
页号-行数
79-3143-倒10149-9158-8、291-倒11302-14
误
用阴极保护。
2……()
正
用阴极保护。
2……()
……现场人工开启型, , , ……现场人工开启型。, 城市
气管道, 在避雷保护范围以外时, ……停电时紧急自动切断阀必须处于关闭状态;
137-倒7、137-倒6
建筑沿外墙架设的燃气管道, ……
8
150-12、(150-13
) ;
(城镇燃气设计规范国家标准管理组 供稿)