2O14年10月 第1O期 天然气管道环境风险评价方法探讨 杨媚 (1.河北省环境科学研究院,河北石家庄050051;2.河北正奇环境科技有限公司,河北石家庄050051) 摘要:指出了天然气是国家大力发展的清洁能源,随着天然气管道覆盖范围的扩大,天然气管道事故数量 也呈现上升态势。结合天然气管道环境风险评价实例,分析了天然气管道在风险识别、源项分析、后果计 算、风险计算和评价以及风险管理等各个阶段的工作要点和注意事项,介绍了射流火灾、天然气泄漏爆炸 事故影响预测模式,以及最大可信事故发生概率、风险值的计算方法。同时,结合天然气管道环境风险评 价工作经验,提出了现有评价方法存在的问题和解决问题的建议。 关键词:天然气管道;环境风险;评价方法 中图分类号:X820。4 文献标识码:A 文章编号:1674 9944(2014)10—0214—05 人身安全与环境的影响和损害,进行评估,提出防范、应 1 引言 天然气是一种优质、清洁的能源,具有使用方便、热 急与减缓措施。环境风险评价与安全评价的主要区别 是:环境风险评价关注点是事故对厂(场)界外环境的影 响,应把事故引起厂(场)界外人群的伤害、环境质量的 恶化及对生态系统影响的预测和防护作为评价工作 重点。 值高、无毒、燃烧后不留残渣等优点,是理想的城市气 源 ]。天然气的使用对于改善能源结构,保护生态环 境,提高城市的现代化水平,走可持续发展的道路具有 重要的意义 ]。根据国家发展战略,天然气作为洁净能 源已被列入国家重点基础设施建设项目。据统计,我国 现有海上和陆上天然气、原有和成品油长距离输送管道 根据《建设项目环境风险评价技术导则》(HJ/T1 69 2004),环境风险评价的基本内容包括:风险识别、源 项分析、后果计算、风险计算和评价、风险管理。 2.1 风险识别 已接近8万km口j,未来几年,我国天然气管道的建设仍 将处于高速发展时期。 使用天然气代替燃煤、燃油,有利于缓解我国现在 严重的大气环境污染状况,也能为居民生活提供便利, 风险识别范围包括生产设施风险识别和生产过程 所涉及的物质风险识别。在进行风险识别时,应经全方 面、全过程的风险识别。生产设施风险识别要包括主要 但天然气的燃爆特性也给人们在利用过程中带来种种 安全隐患。我国在役的天然气管道中,已有5O 管道 服役年龄超过20年,即国家标准规定的设计寿命;随着 服役时间的增长,因管道材质问题或施T、腐蚀和外力 等作用造成的损伤,使管道安全状况日益恶化,潜藏的 安全隐患不断显现,大量的管线已处于事故频发阶 段 。原国家环境保护总局2004年12月11日发布的 生产装置、贮运系统、公用T程系统、T程环保设施及辅 助生产设施等。物质风险识别要包括主要原材料及辅 助材料、燃料、中间产品、最终产品以及生产过程排放的 “三废”污染物等。 2.1.1 物质风险识别 天然气的主要成分为甲烷、乙烷、乙烯等,具有燃爆 特性,具体应根据企业提供的《天然气成分分析报告》进 行分析。要特别关注天然气中H。s的含量。H S是剧 毒物质,天然气的毒性主要来自于H。S。同时,经分析 和比较加拿大、西欧的事故统计结果,发现含硫天然气 《建设项目环境风险评价技术导则》(HJ/T169—2004) 要求,对于涉及有毒有害和易燃易爆物质的生产、使用、 贮运等新建、改建和技术改造项目进行环境风险评价。 因此,环境风险评价已经成为天然气管道项目环境影响 评价的重要内容。 管道事故率比净化天然气事故率高;引发含硫天然气 管道发生事故的最主要因素为管道内腐蚀,比例达到了 53 ;引发净化天然气管道发生事故的最主要因素为第 三方破坏,比例达到了52 要意义。 2.1.2 生产设施风险识别 本文将结合张家口某天然气输气管道建设工程的 环境风险评价,就天然气输气管道环境风险评价采用的 技术方法进行研究探讨。 。因此,确定天然气中 H s的含量对识别风险类型、确定最大可信事故具有重 2环境风险评价的内容 建设项目环境风险评价是指对建设项目建设和运 行期间发生的可预测突发性事件或事故(一般不包括人 长输管线]二程主要包括站场、阀室和输气管线,其 主要_[艺设备一一管道、阀门、过滤分离器等均为密闭 高压系统,生产过程中若设备超压、阀门失灵等原因可 为破坏及自然灾害)引起有毒有害、易燃易爆等物质泄 漏,或突发事件产生的新的有毒有害物质,所造成的对 收稿日期:2014 08—27 能引起设备超压爆炸,若设备密封措施失效可引起天然 作者简介:杨214 媚(1983),女,辽宁沈阳人,工程师,主要从事环境保护、环境影响评价及清洁生产方面的研究。 杨媚:天然气管道环境风险评价方法探讨 气泄漏,天然气遇到点火源会发生火灾、爆炸事故。 2.2源项分析 2.2.1 重大危险源辨识 长输管线天然气存贮量为进出站阀门、输送管道内 两截断阀之间的贮存量。进行重大危险源辨识时,除应 根据管道实际输送天然气量及《建设项目环境风险评价 导则》(HJ/Tl69—2004)附录A规定的贮存场所临界 量以外,还应考虑管道压力。根据重大危险源辨识结果 和《建设项目环境风险评价技术导则》(HJ/T169— 2004)中有关规定确定环境风险评价工作等级。 2.2.2 最大可信事故及其概率 应对国内外天然气管线事故进行统计与分析,找出 引起天然气管道事故的主要类型、产生原因、影响因素、 事故后果,为确定最大可信事故提供依据。 长输管线工程跨越区域周边的环境状况存在较大 的差异,因此,在确定最大可信事故时不能简单地按照 管线长短、压力及天然气储量确定,而应该综合分析天 然气储量、周边环境状况、敏感点分布情况及规模、管道 特性、发生事故的类型和概率,并对各事故因素所占的 权重进行统计、分析,最终确定最大可信事故。以张家 口市某天然气管道项目为例,其冯家沟截断阀室一大仓 村截断阀室段管线长30.1km,天然气储量为221.23t, 管线两侧500m范围内有6个居民点,与管线最近的居 民点为158m;王家楼截断阀室一冯家沟截断阀室段管 线长29.9kin,天然气储量为219.76t,管线两侧500m 范围内有9个居民点,与管线最近的居民点为130m。 综合比较可知,虽然冯家沟截断阀室一大仓村截断阀室 段管线天然气储量较大,但王家楼截断阀室 冯家沟截 断阀室段发生天然气泄漏爆炸事故产生的环境风险更 大,因此应将王家楼截断阀室~冯家沟截断阀室段管线 发生泄漏从而导致火灾、爆炸事故确定为最大可信 事故。 最大可信事故发生概率的确定方法,现在环评中常 用的有两种:事故树法和类比法。事故树法是将最大可 信事故作为顶上事件,通过分析导致最大可信事故发生 的基本事件及各类基本事件的最小割集、最小径集组 合,并为不同危险源因素赋予了不同的危险性重要度 值,经计算确定最大可信事故发生概率。该方法针对性 强、定量性高,但是计算过程复杂,如果基本事件分析不 全面,计算结果也会产生较大的偏差。类比法是现在使 用较多的一种方法,其主要是类比目前权威部门公开发 表的关于各类型事故发生概率的统计结果。该方法简 单易行,但是不同项目在不同环境下发生事故的概率不 可能完全相同,因此,此方法只是对最大可信事故发生 概率的一种估算。 2.2.3 天然气泄漏量计算 泄漏事故规模通常划分为小型、中型、大型及特大 型等几个等级。一般环评选取大型泄漏事故作为评价 对象,即施工作业不慎挖破管道,致使管道严重破裂,输 送介质大量泄漏。 发生泄漏事故时天然气泄漏量采用《导则》中的公 式计算,泄漏速率Qr计算如下: 环境及保护 Q;=YC APMK (2 K+1 式中: 一气体泄漏速度,kg/s;P一容器压力, MPa; 一气体泄漏系数;A一裂口面积,rn ;M一分子 量;R一气体常数; 一气体温度,K。 2.3后果计算 天然气泄漏引发的典型火灾爆炸事故主要是射流 火灾事故以及外泄的天然气与空气形成爆炸性混合物 遇明火引发爆炸事故。 2.3.1喷射火灾影响分析 在一定压力下的天然气泄漏时会形成射流,如果在 泄漏裂口处被点燃,则形成喷射火。把整个喷射火看成 是由喷射中心线上的全部点热源组成,每个点热源的热 辐射通量相等,则产生的热辐射通量按下式计算: 口=,7QH( 式中:口一点热源热辐射通量,W; 效率因子; Q’一泄漏速度,kg/s;H 一燃烧热,J/kg。 从理论上讲,喷射火的火焰长度等于从泄漏口到可 燃混合气燃烧下限(LFL)的射流轴线长度。对表面火 焰热通量,则集中在LFL/1.5处。对风险评价分析而 言,点热源数 一般取5即可。 射流轴线上某点热源i到距离该点z处一点的热 辐射强度为: qR 一4丁r 式中:J 一点热源i至目标点 处的热辐射强度, W/m ;q--点热源的辐射通量,W;R一辐射率2; 一点 热源到目标点的距离,rrl。 某一目标点处的入射热辐射强度等于喷射火的全 部点热源对目标的热辐射强度的总和: ,一∑J 式中: 一计算时选取的点热源数。 火灾产生的热辐射强度足够大时,可使周围的物体 燃烧或变形。强烈的热辐射可能烧毁设备甚至造成人员 伤亡等。火灾损失估算建立在辐射通量与损失等级的 相应关系的基础上,表1为不同入射通量造成伤害或损 失的情况。 表1热辐射的不同入射通最所造成的损失 ( 毛 对设备的损害 对人的伤害 s7.5操作设备全部损坏1: 一一 在无火焰、长时间辐射下, 重大烧伤/lOs, ~ 木材燃烧的最小能量 i00 死亡/1rain … 有火焰时,木材燃烧, 1度烧伤11o s, … 塑料融化的最低能量 l 死亡/1min 20s以上感觉疼痛, … 未必起泡 I.6 …一 长期辐射无不舒服感。 以张家口市某天然气管道为例,假设该管道的泄漏 孔半径为5cm,泄漏速度为62.43kg/s,天然气燃烧热为 37.29×10。J/kg,则发生射流火灾时不同入射通量对应 215 2014年1O月 绿 色 科技 第10期 的伤害半径如表2所示。 对照该管道沿线敏感点情况,如发生射流火灾时, 该工程全线无村庄位于射流火灾的死亡半径、重大烧伤 半径、1度烧伤半径和轻伤半径涉及的范围内,但约有9 个村庄的部分区域处于无不舒服感觉半径涉及的范围 内,会不同程度地受到射流火灾的影响。因此,应采取 可靠的风险防范措施,防止射流火灾的发生,降低射流 火灾的影响。 表2射流火灾伤害半径一览 2.3.2爆炸事故影响分析 天然气管道发生泄漏后,除可能发生射流火灾以 外,也可能发生爆炸事故。目前环境影响评价中一般使 用蒸气云爆炸(UVCE)模型预测泄漏的天然气在理想 状态下产生的伤害程度;但笔者认为,应预测计算天然 气管道泄漏后的爆炸极限范围,结合天然气爆炸冲击波 伤害程度分析结果,说明拟建工程天然气管道泄漏后发 生爆炸对环境的影响情况。 (1)蒸气云爆炸影响分析。蒸气云爆炸(UVCE)是 由于气体或易于挥发的液体燃料的大量快速泄漏,与周 围空气混合形成覆盖很大范围的“预混云”,在某一有限 制空间遇点火而导致的爆炸。根据气云爆炸冲击波预 测模型的特点和复杂性,可以将它们分成数值模型、物 理模型和相关模型三种。其中相关模型也就是缩放比 率模型,依靠试验结果而建立起来的,易应用于灾害风 险评价方面,典型的UVCE相关模型包括TNT模型、 TNO模型、ME模型、CAM模型。长期以来,军事上一 直对高能炸药的破坏可能性很感兴趣,积累了很多 TNT量与破坏之间关系的实验数据。因此,使用TNT 当量来描述事故爆炸的威力就比较方便。 ①TNT当量计算。 WTNT一1.8 V ,Q,/QTN1 式中:1.8一地面爆炸系数;W 一蒸气云的TNT 当量,kg; 一蒸气云的TNT当量系数;W,一蒸气云中 燃料的总质量,kg;Q厂燃料的燃烧热,MJ/kg;Q 一 TNT的爆炸热,MJ/kg。 ②死亡半径R 。 R 一13.6(W r/1000)。” ③重伤半径R 。 由下列方程式求解: △P 一0.137Z-。+0.119Z +0.269Z一一0.019 Z—R /(E/P。) 。 △P =17000/P 式中:P。 环境压力,Pa;E爆炸总能量,E一1.8 A ,Q,×3O ,3O 为假定参与爆炸反应的气体百 分比。 ④轻伤半径R。。 轻伤半径R。由下列方程组求解 AP =0.137Z 十0.1l9Z~。+0.269Z一0.019 21 6 z—R f(E,P .E) zkPs一17000/P。 ⑤财产损失半径R 。 对于爆炸性破坏,财产损失半径R (破坏半径)的 计算公式为: R 一KIIWTNT 邝/(1+(3175/WTNT) ) 。 式中:K 一二级破坏(砖砌房外表5O ~70 破 损,墙壁下部危险)系数,K 一4.6。 根据以上公式对天然气管道泄漏爆炸危害进行计 算,可得冲击波影响下的死亡半径、重伤半径、轻伤半径 和财产损失半径。 (2)天然气泄漏地面浓度分析。天然气泄漏后地面 浓度计算采用下述多烟团模式: Cb(… 一面 唧c— H:、 f(x--z ) (y-- ) 1 2 2 Xp{一 ~ 7 式中: ( ,Y,O,£ )一第i个烟团在t 时刻(即第 w时段)在点(22,Y,O)产生的地面浓度;Q 一烟团排放 量(mg),Q 一QzXt;Q为释放率(mg/s),At为时段长度 (s); , , , , 一烟团在w时段沿37、Y、z方向的 等效扩散参数(m),可由下式估算: 2 一∑ 2.^ ( —z,Y, ) 式中: 2m= 2. (t )一 , (t 1);z 和Ywi~第砌时 段结束时间第i烟团质心的z和Y坐标,由下述两式 计算: z 一 …(£一t 一1)+∑tt (t 一t 1) 一U . (t--t 1)+∑ (t^一t 一1) 各烟团对某个关心点t小时的浓度贡献,按下式 计算: C(x,y,O, )一∑G(z,y,O, ) 式中: 一需要跟踪的烟团数,可由下式确定: C + ( ,y,O, )三三三f∑C (1z, ,O, ) 式中:_厂为小于1的系数,可根据计算要求确定。 环境影响评价按不利情况,采用上述模式预测计算 不同时刻下风向天然气最大落地浓度及出现距离,下风 向天然气浓度大于等于爆炸下限区域范围以及下风向 天然气浓度大于等于爆炸上限区域范围。 天然气浓度达到爆炸极限的范围为浓度介于爆炸 上限和爆炸下限之间的环状区域,并且随着泄漏时间的 延长,爆炸极限范围出现的位置变远,覆盖范围变宽。 在该范围内,如果动火或存在明火会发生爆炸事故。 (3)天然气抬升高度分析。由于天然气管道内压力 高于外界大气压,且天然气密度小于空气密度,因此天 然气从管道泄漏出来后会在动力和热力的作用下继续 向上攀升,达到一定高度后才逐渐沿风向变平。 一般烟气抬升公式由动量项和热量项组成,其中动 量项是指烟气在烟囱内向上具有的动量使它离开烟囱 以后继续上升,热量项是指烟气温度比周围气温高,密 杨 媚:天然气管道环境风险评价方法探讨 度较小,因浮力而上升。本节评价对象为高压天然气管 道的泄漏事故,泄漏天然气的热释放率很低,而天然气 喷出速度接近音速,在这种情况下,动量项对烟气抬升 高度所起的作用远大于热流项的作用 。综合分析,可 选取霍兰德公式来分析天然气泄漏情形下的抬升高度: △H 一K v sa E1.5+2.68×10 sPd( 』S )] 式中:△H,一一校正后的天然气抬升高度;K一霍兰 德公式的修正系数,根据美国ASME对其进行的规定, 保守估计取A、B稳定度等级,K一1.15;U 一天然气的 n 最大泄漏速度,U 一— ; 一泄露孔径,mill; 一平均 nzp0 风速,rn;P~操作压力,kPa;Ts一天然气温度,K; 一 大气温度,K。 以张家口市某天然气管道为例,假设该管道的泄漏 孔半径为5cm,泄漏速度为62.43kg/s;管道操作压力为 6300kPa,根据以上计算得到校正后的天然气抬升高度 为12.85m。 可见,一般天然气管道发生泄漏后,天然气的抬升 高度往往可达10m以上,如果管道泄漏事故发生在空 旷的乡村,在10m左右的高空很少会有明火出现,故引 发火灾爆炸的风险较低;如果管道泄漏事故发生在城 市,由于城市建筑较高,故即使是在高空也可能存在明 火。综合比较可知,天然气管道泄漏事故发生在城市 比发生在农村风险值要高得多。 (4)泄漏后爆炸影响分析。要综合考虑由蒸气云爆 炸影响范围和天然气管道泄漏扩散后达到爆炸极限的 覆盖范围,最终确定泄漏后爆炸影响范围。根据管道沿 线居民区分布情况调查,确定全线管道中心线两侧受影 响的敏感点。 2.4风险计算和评价 风险可表述为: 风险值(后果/时间):概率(事故数/单位时间)× 危险程度(后果/每次事故) 可利用2.2节确定的最大可信事故发生概率和2.3 分析的最大可信事故危害程度计算求得。对于社会公 众而言最大可接受风险不应高于常见的风险值。 风险的单位多采用“死亡人数/年”。安全和风险是 相伴而生的,风险事故的发生频率不可能为0。通常事 故危害所致风险水平可分为最大可接受水平和可忽略 水平。表3列出了一些机构和研究者推荐的最大可接 受风险水平和可忽略水平。 在工业和其他活动中,一般而言,环境风险的可接 受程度对有毒有害工业以自然灾害风险值(即10 /a) 为背景值。 2.5风险管理 风险管理包括风险防范措施和应急预案两部分 内容。 根据《建设项目环境风险评价技术导则》(HJ/T169 2004),风险防范措施可从选址、总图布置和建筑安全 防范措施、危险化学品贮运安全防范措施、工艺技术设 环境及保护 表3 最大可接受水平和可忽略水平的推荐值 计安全防范措施、自动控制设计安全防范措施、电气电 讯安全防范措施、消防及火灾报警系统、紧急救援站或 有毒气体防护站设计等方面提出。但是笔者认为,天然 气管道工程在施工期、运营期应采取的安全防范措施是 不同的,如果单单按照以上方法对安全防范措施进行分 类,则提出的措施会显得混乱,因此应按照工程前期设 计阶段、施工期、运营期分别提出风险管理要求,同时, 对于不同时期的安全防范措施,应按照预防、减缓、恢复 三个层次提出,使防范措施一目了然,以利于将来的竣 工环保验收。 应急预案对于降低事故危害具有重要作用。向环 境排放污染物的企业事业单位,生产、贮存、经营、使用、 运输危险物品的企业事业单位,产生、收集、贮存、运输、 利用、处置危险废物的企业事业单位,以及其他可能发 生突发环境事件的企业事业单位,均应当编制独立的环 境应急预案文件,并组织评估小组对本单位编制的环境 应急预案进行评估,根据评估结果,对应急预案草案进 行修改后报所在地环境保护主管部门备案。因此,环评 中只对天然气管线工程如何建立事故应急预案进行概 括性描述,从事故应急预案制定原则、应急预案主要内 容等方面提出原则性的要求,供生产单位在编制事故应 急预案时参考。 3问题与建议 3.1 问题 (1)目前,在进行天然气管道泄漏事故情景预设时, 往往假定天然气输气管道整体断裂(此为极端最不利情 况),但是在实际生活中,管道整体破损情况多是由于在 施工活动中的人为破坏或是自然灾害;埋地管道如焊缝 不严或腐蚀等原因,很难造成管道整体断裂。但是根据 《建设项目环境风险评价技术导则》(HJ/T169—2004) 的有关规定,建设项目环境风险评价一般不包括人为破 坏及自然灾害。人为破坏和自然灾害造成的管道破坏 应属于安全评价范畴。 (2)一般情况下,埋地管道发生泄漏后,泄漏的天然 气会经管道上的覆土慢慢地渗漏到地表,其泄漏速度比 露天容器要小得多。现有风险评价泄漏预测模型为考 虑管道覆土对泄漏气体的阻力,因此并不完全适用埋地 管道。 (3)国内外天然气管道事故统计是天然气管道最大 可信事故概率分析的类比依据,而目前国内缺乏这方面 的统计资料 ]。 217 2Ol4年1O月 绦色科技 Journal of Green Science and Technology 第1O期 环境保护工作中开展污染企业信用评级的探讨 韩庆蓉 (山西省太原市环保局迎泽分局环境监测站,山西太原030012) 摘要:通过对现有环境保护工作中污染企业风险事件管理的现状分析,对污染企业信用评级的概念、内涵、 模式进行了研究,探讨了开展污染企业信用评级的可能性,并就相关技术问题进行了讨论,提出了在污染 企业开展信用评级的政策性建议。 关键词:污染企业;信用评级;建议;技术 中图分类号:F206 文献标识码:A 文章编号:1674—9944(2014)10—0218—03 2.1.1 传统信用评级的定义 1 引言 为进一步推进环境保护工作的发展,2o13年,由环 保部起草的《企业环境信用评价办法(试行)》征求意见 稿征求意见。该征求意见稿显示,我国将对污染物排放 总量大、环境风险高、生态环境影响大的10类企业进行 环境信用评级。其中,未批先建、恶意偷排及构成环境 犯罪等将被直接“一票否决”评为环保不良企业。 信用评级(Credit Rating),又称资信评级,是一种 社会中介服务,为社会提供资信信息,或为单位自身提 供决策参考。 通常意义的信用评级,是指独立的第三方信用评级 中介机构对债务人如期足额偿还债务本息的能力和意 愿进行评价,并用简单的评级符号表示其违约风险和损 失的严重程度。广义的信用评级则是对评级对象履行 相关合同和经济承诺的能力和意愿的总体评价 ]。 与之对应的,环保行业的信用评级研究工作尚未展 开,相关的评级原理、评级标准和方法、评级指标体系、 评级操作系统及评级数据库等等的研究,基本处于空白 状态,原因在于,传统的信用评级是指信用评级机构按 照规范化的程序,采用一定的分析技术,给出关于受评 对象的信用风险的判断意见,并将结果用一个简单的符 号表示出来。与之相对应的,评级机构是以盈利为目的 从概念上讲,信用评级与环境保护没有直接的关 系。但是,随着环境保护工作发展的需要,如何对重点 污染行业,包括重点污染行业中的企业进行信用评级, 为其建立信用档案,通过为社会提供正确的环境保护信 用信息,有效防止环境风险事件的发生,并积极引导企 业通过信用评级进行自我管理和约束,积极推进社会生 态文明的发展,有重要的意义。 2.1.2 污染企业信用评级的含义 的市场参与者,在盈利与声誉之问找平衡,评级机构不 承担评级结果的法律责任,因是咨询机构,不是权威的 鉴证类机构。 污染企业信用评级,指的是由独立的第三方中介机 构,通过信用评级的技术,对评级对象的环境保护状况 进行评价及预测,为社会提供准确的企业环保信用 信息。 如何在环境保护工作中有效开展信用评级研究,并 能将研究成果迅速转化为评级能力,是中国环保工作者 及信用评级工作者值得深思的问题。本文将就环境保 护信用评级的概念、特点、途径、建议等方面展开探讨。 污染企业信用评级主要包括以下几方面含义。 2概念和特点 2.1 概念 收稿日期:2O14一O9—10 (1)环保信用与其他信用关系指向的对象一样,表 现形式为带有约束关系的义务关系,这种关系带有契约 性质,其中,评级对象的责任是向社会提供环境保护公 作者简介:韩庆蓉(1972~),山西沁源人,工程师,主要从事环保工作。 3.2 建议 与环境工程,2006,13(4):63~65. (1)尽快建立我国天然气管道事故数据库,统计各 类天然气管道事故发生概率、事故原因和影响范围,为 风险评价提供技术支撑。 (2)尽快开发适用于埋地天然气管道泄漏、爆炸等 事故的预测模型,使天然气管道事故影响预测更符合 实际。 [2]贺力荃,彭小武.城镇天然气管道输送工程环境风险评价探讨 [J].新疆环境保护,2008,30(4):17~2o. [3]蒲 明.2009年我国油气管道新进展[J].国际石油经济,201l (3). [43雷 达.天然气输气管道泄漏事故风险分析[D].北京:北京化工 大学,2012. [5]罗小兰,向启贵,银小兵,等.关于天然气管道环境风险评价的认识 [J].石油与天然气化工,2008(6). 参考文献: [1]刘 斐,刘 茂.城市输气管线火灾事故的风险定量计算[J].安全 218 [6]吕 奎,孔益平,周 宇,等.天然气管道风险分析及应急管理[J]. 安全与环境工程,2Ol2(2).