CN103455724B - 一种含硫天然气管道的泄漏中毒潜在影响区域计算方法 - Google Patents
一种含硫天然气管道的泄漏中毒潜在影响区域计算方法 Download PDFInfo
- Publication number
- CN103455724B CN103455724B CN201310403544.8A CN201310403544A CN103455724B CN 103455724 B CN103455724 B CN 103455724B CN 201310403544 A CN201310403544 A CN 201310403544A CN 103455724 B CN103455724 B CN 103455724B
- Authority
- CN
- China
- Prior art keywords
- pipeline
- natural gas
- sulfur
- gas pipeline
- potential impact
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Pipeline Systems (AREA)
Abstract
本发明公开了一种含硫天然气管道的泄漏中毒潜在影响区域计算方法,包括以下步骤:步骤一,根据输气管道设计资料和现场踏勘,收集管道的管径、壁厚、运行压力、天然气含硫量、所处地区的地表环境等基础数据;步骤二,根据基础数据对计算潜在影响区所需的相关参数进行计算;步骤三,将计算得到的参数值代入到含硫天然气管道的泄漏中毒潜在影响区域计算公式中,最终确定潜在影响区域的范围。本发明的计算方法简便,计算结果能表达潜在影响区的具体形状和位置,完善了含硫天然气管道泄漏潜在影响区的计算理论和计算结果,为含硫天然气管道的完整性管理提供了管治依据。
Description
技术领域
本发明属于天然气管道完整性管理领域,尤其涉及一种含硫天然气管道的泄漏中毒潜在影响区域的计算方法。
背景技术
天然气是一种清洁、高效的能源,而管道输送是众多天然气运输方式中最为经济的方式之一。目前,我国建成的天然气管道总长已超过8万公里。随着经济社会的发展,我国对天然气资源的需求也持续增加,众多已探明的含硫天然气矿也得以开采和利用。但是,含硫天然气对输送管道的腐蚀能力远大于净化天然气,在管道的运营过程中,因腐蚀所产生的管道泄漏事故也时有发生。因此,输送含硫天然气管道的完整性管理越来越得到国内外管道运营者的重视。
管道完整性管理的基础是对管道高后果区的划分,而管道风险的潜在影响区域(Potential ImpactArea)则是管道高后果区的划分依据。受国际公认的完整性管理标准ASME B31.8S《Managing Integrity System ofGas Pipeline》和API1160《ManagingSystem Integrity forHazardous Liquid Pipelines》,以及我国的行业内标准SY/T6621《输气管道系统完整性管理》均定义管道风险的潜在影响区域是指管道失效后果的量化指标,用于表征受失效后果危害影响的区域范围。目前,以上三种标准均采用相同的潜在影响区域计算公式,其表述如下:
式中r——潜在影响区半径,m;
d——天然气管道的外直径,m;
p——天然气管道的输送压力,kPa。
据文献《Amodel for sizing high consequence areas associated withnatural gas pipelines.》(Mark J.Stephens.C-FER,Report99068.Edmonton,Canada,2000.)论述,该公式建立的核心假设是:从管道中泄漏的天然气在泄漏点位置被立即点燃;该公式建立的理论基础是:火灾所带来的热辐射导致人、物的损害。
分析天然气管道泄漏事件不难发现,“立即点燃”并非该事件的唯一后果。当泄漏点附近不存在点火源时,例如泄漏点位于水稻田附近,泄漏的天然气将不会被立即点燃,而是在扩散作用下,于泄漏点下风方向形成一个可燃气体覆盖区域。另一方面,对于含硫天然气而言,火灾热辐射也并非造成人员伤亡的唯一原因。由于H2S具有极强的毒性,人体接触极微量的H2S就可能中毒身亡。因此,含硫天然气在未点燃条件下所形成的扩散云区域内,还有可能引发人员中毒的后果。显然,要表达这一中毒潜在影响区域应基于大气扩散理论,以及H2S的人员中毒剂量。由此可见,目前标准中所采用的潜在影响区域计算公式,从建立的假设到理论基础都不能表达含硫天然气在延迟点燃条件下所引发的人员中毒后果,其计算结果存在缺陷。
发明内容
本发明提供一种含硫天然气管道的泄漏中毒潜在影响区域的计算方法,以解决现有技术中的各种缺点不足,本发明采取的技术方案如下:
提供一种含硫天然气管道的泄漏中毒潜在影响区域计算方法,该方法包括以下步骤:
步骤一,根据输气管道设计资料和现场踏勘,收集管道的管径、壁厚、运行压力、天然气中的含硫量、所处地区的地表环境、地表平均风速、天气条件、云量和太阳辐射强度基础数据;
步骤二,根据基础数据对计算潜在影响区所需的管道所处地区的大气稳定程度、管道所处地区的地表粗糙度、管道所处地区的扩散参数进行计算;
步骤三,将计算得到的参数值代入到含硫天然气管道的泄漏中毒潜在影响区域计算公式中,其中,公式的约束条件为:最终确定潜在影响区域的范围。
本发明提供的含硫天然气管道泄漏扩散潜在影响区域的计算方法,是以大气扩散理论为基础,以H2S的短期接触中毒阈值为边界条件所推导出的新计算公式,解决了泄漏的含硫天然气在未立即点燃的条件下,在大气中扩散形成的人员中毒潜在影响区域的计算问题,涵盖了原标准方法未涉及的人员中毒后果,完善了含硫天然气管道泄漏潜在影响区的计算理论和计算结果。本发明的计算方法简便,计算结果能表达潜在影响区的具体形状和位置,可以为含硫天然气管道的完整性管理提供重要的管治依据。
附图说明
图1是本发明实施例提供的天然气管道泄漏扩散潜在影响区域的计算方法的流程图。
图2是本发明方法计算得到的潜在影响区域与原标准计算的影响区域对比。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明提供一种含硫天然气管道的泄漏中毒潜在影响区域计算方法,该方法包括以下步骤:
步骤一,根据输气管道设计资料和现场踏勘,收集管道的管径、壁厚、运行压力、天然气含硫量、所处地区的地表环境等基础数据;
步骤二,根据基础数据对计算潜在影响区所需的相关参数进行计算;
步骤三,将计算得到的参数值代入到含硫天然气管道的泄漏中毒潜在影响区域计算公式中,最终确定潜在影响区域的范围。
在步骤一中,管道基础参数收集包括管道的管径、壁厚、运行压力、天然气中的含硫量、所处地区的地表环境、地表平均风速、天气条件、云量和太阳辐射强度。
在步骤二中,计算潜在影响区域的相关参数值包括:确定管道所处地区的大气稳定程度;确定管道所处地区的地表粗糙度;计算管道所处地区的扩散参数。其中,管道所处地区的大气稳定程度按表1查询确定。
表1大气稳定度级别划分表
*(1)夜晚的定义为日落前一小时至次日日出后一小时;(2)无论天气条件如何,夜晚前后一小时内的大气稳定度均为D级(中性);(3)强日照是指仲夏晴天正午的日照强度,弱日照是指寒冬晴天正午的日照强度。
其中,管道所处地区的地表粗糙度按表2查询确定。
表2地表有效粗糙度
地面类型 | 地表粗糙度Z0/m |
草原、平坦开阔地 | <0.1 |
农作物地区 | 0.1—0.3 |
村落、分散的树林 | 0.3—1 |
分散的高、矮建筑物 | 1—4 |
密集的高、矮建筑物 | 4 |
其中,管道所处地区的扩散参数计算公式为:
σy=σy0fy
σz=σz0fz
式中x,y,z——直角坐标系中的三个主轴,x轴与风向平行,其正方向表示风向,y轴与风速垂直,z轴为竖直方向,与重力加速度方向平行;
σy,σz——天然气在y、z方向上的扩散参数,m;
σy0、σz0——初始扩散参数,按表3查询计算公式,公式中的x表示计算位置与天然气泄漏点间的距离在x轴方向上的投影,m;
fy、fz——扩散参数修正系数,若地表粗糙度小于0.1则扩散修正系数均为1,若地表粗糙度大于0.1,则有:
fy(Z0)=1+a0Z0
a0,b0,c0,d0,e0,f0,g0——计算常数,按表4查询。
表3初始扩散参数计算表
大气稳定度 | σy0(m) | σz0(m) |
A | 0.22x(1+0.0001x)-1/2 | 0.2x |
B | 0.16x(1+0.0001x)-1/2 | 0.12x |
C | 0.11x(1+0.0001x)-1/2 | 0.08x(1+0.0002x)-1/2 |
D | 0.08x(1+0.0001x)-1/2 | 0.06x(1+0.0015x)-1/2 |
E | 0.06x(1+0.0001x)-1/2 | 0.03x(1+0.0003x)-1 |
F | 0.04x(1+0.0001x)-1/2 | 0.016x(1+0.0003x)-1 |
表4计算常数表
进一步,在步骤三中,计算扩散潜在影响区域的公式为:
其中,公式的约束条件为:
式中x,y——天然气管道泄漏扩散潜在影响区的边界点坐标;
C(x,y)——扩散天然气中H2S在坐标位置(x,y)处的浓度,在计算中取为H2S的短期接触中毒阈值,C(x,y)=10×10-3kg/m3。
η——天然气中H2S的含量,kg/m3;
d——管道直径,m;
δ——管道的壁厚,m
p——管道运行压力,kPa;
σy,σz——天然气在y、z方向上的扩散参数,m;
u——地表平均风速,m/s。
其中,潜在影响区域计算公式迭代求解的方法是:设定增长步长Δx,以天然气泄漏点为原点,依次将坐标(x+n·Δx)(n=1,2,3,…)代入计算公式中,求得特定浓度下的对应坐标y,至到(x+n·Δx)值不再符合边界条件公式,据此可得到一系列潜在影响区域的边界点坐标(xn,yn)。(n=1,2,3,…),将这些坐标用曲线连接起来就可以得到天然气管道泄漏扩散的潜在影响区域。
下面结合具体实例对本发明的应用原理作进一步描述。
例:一条输送天然气的管道,管径为273mm,壁厚为10mm,工作压力为7.36MPa,天然气中的含硫量为3.36×10-3kg/m3,管道泄漏点地表为农作物地区,平均风速为3m/s,以本发明所给出的方法,计算晴朗白天时,中等太阳辐射条件下的泄漏潜在影响区。
第一步:收集管道基础参数:管道的管径为273mm,壁厚为10mm,运行压力7.36MPa,所处地区的地表为农作物地区、地表平均风速为3m/s、天气条件为晴朗的白天和太阳辐射强度为中等。
第二步:根据基础参数对计算潜在影响区所需的相关参数进行计算。
(1)确定管道所处地区的大气稳定程度
根据表1查得,管道所处地区的大气稳定程度为B级。
(2)确定管道所处地区的地表粗糙度
根据表2查得,管道所处地区的地表粗糙度Z0=0.3。
(3)计算管道所处地区的扩散参数
根据表3查得,管道所处地区的初始扩散参数计算公式为:
σy0=0.16x(1+0.0001x)-1/2
σz0=0.12x
由于粗糙度Z0>0.1,因此需要对初始扩散参数进行修正,则有:
σy=σy0(1+a0Z0)=0.16x(1+0.0001x)-1/2(1+a0Z0)
根据表3查得,上式中的常数为:a0=0.115,b0=1.5,c0=0.045,d0=0.853,e0=0.0128,f0=0.156,g0=0.0136。
第三步,将得到的各项参数带入到潜在影响区域的计算公式中进行迭代计算。
(1)设定泄漏点为坐标原点,设定迭代步长增量Δx=8m;
(2)将x1=x0+Δx=8m带入约束条件,判定得到成立;
(3)将x1=8m以及第一、二步中得到的各项参数带入:
中,计算得到y1=±4.36m;
(4)令xi+1=xi+Δx,重复步骤(2)–(3),直至不成立,退出循环,整理(xi,yi)可得计算结果如表5所示。
表5应用实例的计算结果
X轴坐标 | Y轴正向 | Y轴负向 | X轴坐标 | Y轴正向 | Y轴负向 |
0 | 0 | 0 | 128 | 34.06 | -34.06 |
8 | 4.36 | -4.36 | 144 | 34.73 | -34.73 |
16 | 8.78 | -8.78 | 152 | 34.79 | -34.79 |
24 | 12.11 | -12.11 | 160 | 34.66 | -34.66 |
32 | 15.09 | -15.09 | 176 | 33.75 | -33.75 |
48 | 20.21 | -20.21 | 184 | 32.94 | -32.94 |
56 | 22.43 | -22.43 | 192 | 31.86 | -31.86 |
64 | 24.44 | -24.44 | 208 | 28.73 | -28.73 |
80 | 27.89 | -27.89 | 216 | 26.54 | -26.54 |
88 | 29.35 | -29.35 | 224 | 23.77 | -23.77 |
96 | 30.63 | -30.63 | 240 | 15.22 | -15.22 |
112 | 32.68 | -32.68 | 248 | 6.02 | -6.02 |
120 | 33.46 | -33.46 | 254 | 0 | 0 |
根据原有标准规定,该管道的潜在影响区域为圆形,其半径为:
将两种方法计算得到的潜在影响区域进行作图对比,其结果如图2所示。从图2中可以看出,本发明方法计算得到的潜在影响区域与原标准计算的结果有较大不同。含硫天然气泄漏后的中毒潜在影响区域近似为椭圆形,且在下风方向上的扩散距离远大于其扩散宽度。使用原标准方法,不能表征这一区域,说明本发明确有其必要性和实用性。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种含硫天然气管道的泄漏中毒潜在影响区域计算方法,包括以下步骤:
步骤一,根据输气管道设计资料和现场踏勘,收集管道的管径、壁厚、运行压力、天然气中的含硫量、所处地区的地表环境、地表平均风速、天气条件、云量和太阳辐射强度基础数据;
步骤二,根据基础数据对计算潜在影响区所需的管道所处地区的大气稳定程度、管道所处地区的地表粗糙度、管道所处地区的扩散参数进行计算;
步骤三,将计算得到的参数值代入到含硫天然气管道的泄漏中毒潜在影响区域计算公式中,其中,公式的约束条件为:最终确定潜在影响区域的范围;
其中式中x,y——天然气管道泄漏扩散潜在影响区的边界点坐标;
C(x,y)——扩散天然气中H2S在坐标位置(x,y)处的浓度,在计算中取为H2S的短期接触中毒阈值,C(x,y)=10×10-3kg/m3;
η——天然气中H2S的含量,kg/m3;
d——管道直径,m;
δ——管道的壁厚,m;
p——管道运行压力,kPa;
σy,σz——天然气在y、z方向上的扩散参数,m;
u——地表平均风速,m/s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310403544.8A CN103455724B (zh) | 2013-09-06 | 2013-09-06 | 一种含硫天然气管道的泄漏中毒潜在影响区域计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310403544.8A CN103455724B (zh) | 2013-09-06 | 2013-09-06 | 一种含硫天然气管道的泄漏中毒潜在影响区域计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103455724A CN103455724A (zh) | 2013-12-18 |
CN103455724B true CN103455724B (zh) | 2017-02-15 |
Family
ID=49738078
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310403544.8A Expired - Fee Related CN103455724B (zh) | 2013-09-06 | 2013-09-06 | 一种含硫天然气管道的泄漏中毒潜在影响区域计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103455724B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112115623B (zh) * | 2020-10-20 | 2022-03-15 | 西南石油大学 | 一种泄漏工况下输气管道阀室压降速率计算方法 |
CN114441513B (zh) * | 2020-11-03 | 2024-02-09 | 中国石油化工股份有限公司 | 一种含硫化氢设备泄漏扩散距离的估算方法 |
CN113343498A (zh) * | 2021-07-05 | 2021-09-03 | 西南石油大学 | 一种含h2s/co2环境管道钢腐蚀速率计算方法 |
CN117437329B (zh) * | 2023-12-12 | 2024-03-01 | 江西云眼视界科技股份有限公司 | 高后果区域绘制方法、系统、可读存储介质及计算机设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102486261A (zh) * | 2010-12-03 | 2012-06-06 | 中国石油天然气股份有限公司 | 输气管道调度方案综合评价方法 |
-
2013
- 2013-09-06 CN CN201310403544.8A patent/CN103455724B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102486261A (zh) * | 2010-12-03 | 2012-06-06 | 中国石油天然气股份有限公司 | 输气管道调度方案综合评价方法 |
Non-Patent Citations (6)
Title |
---|
A model for sizing high consequence areas associated with natural gas pipelines;Stephens MJ;《2002 4th international pipleine conference》;20021231;759-767 * |
天然气管道泄漏H_2S扩散及影响区域的数值模拟与分析;罗志云;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20080815;B019-40 * |
输气管道稳态计算公式改进研究;李长俊等;《西南石油大学学报(自然科学版)》;20081031;第30卷(第5期);155-157 * |
静置油罐的小呼吸损耗量在大气环境中的扩散;段贵话;《石油库与加油站》;20120630;第21卷(第3期);16-19 * |
高含硫天然气管线泄漏毒害研究;吴昊;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20110415;B019-322 * |
高含硫天然气管道泄漏毒害评价模型分析研究;梁瑞等;《化学工程与装备》;20100630(第6期);34-37 * |
Also Published As
Publication number | Publication date |
---|---|
CN103455724A (zh) | 2013-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103455724B (zh) | 一种含硫天然气管道的泄漏中毒潜在影响区域计算方法 | |
Nazir et al. | Impacts of renewable energy atlas: Reaping the benefits of renewables and biodiversity threats | |
Li et al. | Rapid urbanization effect on local climate: intercomparison of climate trends in Shenzhen and Hong Kong, 1968-2013 | |
Li et al. | Response of energy consumption for building heating to climatic change and variability in Tianjin City, China | |
CN104298851A (zh) | 一种用于预报强降水天气的数据处理方法 | |
Zhang et al. | The source apportionment of primary PM2. 5 in an aerosol pollution event over Beijing-Tianjin-Hebei region using WRF-Chem, China | |
Chuai et al. | Carbon neutrality check in spatial and the response to land use analysis in China | |
Wang et al. | Research on the wind pressure coefficient in natural wind calculations for extra-long highway tunnels with shafts | |
Park et al. | Air conditioning system design to reduce condensation in an underground utility tunnel using CFD | |
Yuan et al. | Aerosol vertical mass flux measurements during heavy aerosol pollution episodes at a rural site and an urban site in the Beijing area of the North China Plain | |
CN108831116A (zh) | 输电线路山火蔓延行为多级校正方法、系统及存储介质 | |
CN110955744A (zh) | 一种基于多源数据的输电线路山火监测预警方法及系统 | |
Popovicheva et al. | Black carbon seasonal trends and regional sources on Bely Island (Arctic) | |
CN115618718A (zh) | 一种基于表面监测和遗传算法优化的垃圾填埋场甲烷排放估算方法 | |
Wong et al. | Development of an integrated Micro-environmental monitoring system for construction sites | |
Li et al. | Ventilation and environmental control of underground spaces: a short review | |
CN107014743A (zh) | 一种新型的埋地钢管腐蚀速率的获取方法 | |
Tidblad et al. | Mapping of acid deposition effects and calculation of corrosion costs on zinc in China | |
Wang et al. | A method for calculating the moisture load in the tailrace surge tank of underground hydropower station based on field tests | |
Yıldırım et al. | Evaluation of cities in the context of energy efficient urban planning approach | |
Liu et al. | Using insulation in china’s buildings: potential for significant energy savings and carbon emission reductions | |
An et al. | Building wind corridor network using roughness length | |
Xu et al. | Modeling and simulation of haze process based on Gaussian model | |
Techarat et al. | Coal dust dispersion from a coal storage pile | |
Charloth et al. | Analysis of particulate dispersion in Suralaya coal-fired power plant |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170215 Termination date: 20190906 |
|
CF01 | Termination of patent right due to non-payment of annual fee |