摘要 含硫天然氣開發(fā)風(fēng)險主要來自于井噴或地面設(shè)施失效時硫化氫泄漏擴散引起的中毒。為了開展復(fù)雜地形條件下含硫天然氣開發(fā)定量風(fēng)險評價，在含硫天然氣泄漏擴散模擬中引入計算流體力學(xué)，運用復(fù)雜地形數(shù)字高程模型和網(wǎng)格劃分技術(shù)，模擬事故狀態(tài)下硫化氫的時空分布，并基于定量風(fēng)險評價基本原理集成了復(fù)雜地形條件下的定量風(fēng)險評價技術(shù)。以位于復(fù)雜地形下的某含硫氣田發(fā)生天然氣井噴泄漏事故為例，運用“含硫天然氣泄漏擴散中毒定量風(fēng)險評價軟件”進行了個人風(fēng)險值和社會風(fēng)險值的計算。結(jié)果表明，基于數(shù)值模擬的定量風(fēng)險評價技術(shù)較常用的拉格朗日煙團模型計算出的風(fēng)險等值線，更能夠客觀真實地反映復(fù)雜地形條件下的風(fēng)險水平分布。

關(guān)鍵詞  川渝地區(qū)  復(fù)雜地形  定量風(fēng)險評價  計算流體力學(xué)  事故概率  風(fēng)險可接受標準

從20世紀70年代開始，部分國家的核工業(yè)及化工行業(yè)開始應(yīng)用定量風(fēng)險技術(shù)(QRA)并在20世紀90年代逐步規(guī)范化，在英國、荷蘭、新加坡、香港及澳大利亞城市地區(qū)的人口密集區(qū)，已將定量風(fēng)險評價作為高風(fēng)險行業(yè)風(fēng)險評估與控制的法律手段。我國在2011年發(fā)布的安監(jiān)總局40號令《危險化學(xué)品重大危險源監(jiān)督管理暫行規(guī)定》也要求采用定量風(fēng)險評價技術(shù)對重大危險源實施評估和監(jiān)管。

硫化氫是一種劇毒、可燃氣體，含硫天然氣特別是高含硫天然氣開發(fā)主要風(fēng)險來自事故導(dǎo)致硫化氫釋放的中毒風(fēng)險。在天然氣勘探開發(fā)過程中一旦發(fā)生事故，大量含硫天然氣泄漏擴散會導(dǎo)致人員中毒死亡。據(jù)不完全統(tǒng)計，近l0年來我國石油天然氣行業(yè)超過70％的傷亡事故是因為硫化氫中毒所致，因此，針對含硫天然氣勘探開發(fā)開展定量風(fēng)險評價顯得尤為重要。

1 川渝地區(qū)開展QRA的特殊性

1.1 復(fù)雜地形的影響

復(fù)雜地形條件下的大氣邊界層特征和污染物擴散規(guī)律一直是邊界層氣象學(xué)研究的難點與熱點。四川盆地四周為海拔l 000～4 000 m的山地所環(huán)抱，四周叢山阻隔，地形閉塞，冬季北來的冷空氣難以浸入，夏季無臺風(fēng)影響，素有“死水區(qū)”之稱^[1]。盡管盆地大氣混合層雖高，但其復(fù)雜的地形特征形成了特殊的氣象條件——靜風(fēng)頻率高、污染系數(shù)大，是全國大氣凈化能力最弱的地區(qū)之一。

氣象和地形條件對污染物的擴散過程有著決定性的影響。由于山區(qū)起伏地形和不均勻下墊面的強烈影響，低層大氣在水平和垂直兩個方向上形成特殊的風(fēng)場和溫度場，如靜風(fēng)、逆溫等，直接影響大氣污染物在近地面的分布、擴散、沉降等情況，導(dǎo)致污染物不易擴散。

目前，川渝地區(qū)分布多處高含硫天然氣區(qū)塊，這些區(qū)塊主要為山地丘陵地形，而且人口分布稠密。由于復(fù)雜地形對意外釋放含硫天然氣的擴散有阻擋作用，局地風(fēng)場跟地形的耦合度強，使位于丘陵山區(qū)的氣體擴散過程變得非常復(fù)雜。筆者運用復(fù)雜地形數(shù)字高程模型和網(wǎng)格劃分技術(shù)，引入流體力學(xué)計算方法，在高性能計算機系統(tǒng)上并行運算，能夠較真實地模擬受復(fù)雜地形和局地風(fēng)場影響的含硫天然氣云團擴散過程。圖1為復(fù)雜地形數(shù)字高程圖。

1.2 泄漏概率的確定

QRA方法的核心是對事故發(fā)生概率和事故后果進行擬合^[2]。目前，國內(nèi)外在事故后果方面的研究較多，認識也比較成熟，而在事故發(fā)生概率方面的研究則相對較少。氣田集氣場站、天然氣凈化廠、集氣管道等設(shè)備失效概率可參考的數(shù)據(jù)庫包括英國HSE烴類泄漏數(shù)據(jù)庫(UK HSE)、美國化學(xué)流程工業(yè)設(shè)備可靠性數(shù)據(jù)庫(CCPS)、基于風(fēng)險的檢驗導(dǎo)則(API 581)、英國陸上管道作業(yè)者協(xié)會(UKOPA)、歐洲氣體管道事故數(shù)據(jù)庫(EGIG)等。中國石油西南油氣田公司(以下簡稱西南油氣田公司)從20世紀60年代開始在含硫天然氣開發(fā)方面積累了一定的經(jīng)驗，也收集了相關(guān)設(shè)備失效數(shù)據(jù)，形成了不同氣質(zhì)條件下的管道失效數(shù)據(jù)庫。筆者所述事故概率通過借鑒并修正采用國外失效數(shù)據(jù)庫與西南油氣田公司同類設(shè)備歷史統(tǒng)計數(shù)據(jù)相結(jié)合的方式來確定。借鑒國外數(shù)據(jù)時考慮擬評價對象的工藝、介質(zhì)、自然和社會環(huán)境狀況，分析數(shù)據(jù)來源及其可靠性，應(yīng)用事故樹模型進行裝置事故發(fā)生概率計算，然后結(jié)合項目設(shè)備及風(fēng)險管理實際情況進行修正采用^[3]。

1.3 風(fēng)險可接受標準的確定

目前國內(nèi)除《危險化學(xué)品重大危險源監(jiān)督管理暫行規(guī)定》確定了重大危險源對人口密集場所、學(xué)校等敏感目標的個體風(fēng)險標準值以及重大危險源的社會風(fēng)險F-N曲線外，還沒有其他適用于偏遠地區(qū)、人居分散條件下含硫天然氣開發(fā)的個人風(fēng)險標準值可供參考。為此，研究英國、荷蘭等國家風(fēng)險標準值制訂的歷史、背景等資料后^[4]，分析我國近10年采掘業(yè)事故傷亡率、公眾意外死亡風(fēng)險以及石油天然氣行業(yè)傷亡情況，提出了員工個人風(fēng)險(1.5×10^-3)、公眾個人風(fēng)險(0.8×10^-4)以及社會風(fēng)險的基準值l.4×10^-3 (N=1)，斜率為-1。風(fēng)險可接受標準的最終確定應(yīng)由相應(yīng)的風(fēng)險承擔(dān)方與風(fēng)險制造方在客觀風(fēng)險基礎(chǔ)上平等協(xié)商決定。

2 基于計算流體力學(xué)的后果計算

計算流體動力學(xué)(Computational Fluid Dynamics，簡稱CFD)是把原來在時間域及空間域上連續(xù)的物理量的場，如速度場和壓力場，用一系列有限個離散點上的變量值集合來代替，通過一定的原則和方式建立起關(guān)于這些離散點上各場變量之間關(guān)系的代數(shù)方程組，然后求解代數(shù)方程組獲得各場變量的近似值^[5]。利用基于CFD原理的有限元軟件來完成復(fù)雜地形條件下含硫天然氣井噴或地面設(shè)施泄漏所致的含硫天然氣三維擴散數(shù)值模擬。

CFD可以看作是在流動基本方程(質(zhì)量守恒方程、動量守恒方程、能量守恒方程)控制下對流動數(shù)值的模擬。通過這種數(shù)值模擬，可以得到極其復(fù)雜問題的流場內(nèi)各個位置上基本物理量(如速度、壓力、溫度、濃度等)的分布，以及這些物理量隨時間的變化情況。由于含硫天然氣是混合物質(zhì)，因此，控制方程還應(yīng)包括組分守恒方程。4個方程的表達式如下。

式中ρ為密度；t為時間；u為速度矢量；μ為動力黏度；p為壓力；T為溫度；k為傳熱系數(shù)；C_p為比熱容；C為體積濃度；D為擴散系數(shù)；S_u為源項。

對于含硫天然氣擴散引起的風(fēng)險，需要計算特定位置硫化氫含量的變化情況，即各網(wǎng)格點或敏感目標處硫化氫含量隨時間變化情況。分別以某井場井噴和某天然氣凈化裝置事故為例，得到事故發(fā)生時硫化氫的時空分布情況(圖2～5)。

3 定量風(fēng)險計算實例

以某高含硫氣田在鉆井過程中發(fā)生井噴失控事故為例，利用具有自主知識產(chǎn)權(quán)的定量風(fēng)險計算軟件(QRA-AYY)，通過建立事故模型，導(dǎo)入項目地形圖，輸入項目所在地的大氣穩(wěn)定度、風(fēng)向及風(fēng)速聯(lián)合頻率，模擬計算泄漏出來的含硫天然氣在不同氣象條件下的擴散影響區(qū)域范圍，然后輸入人口分布及關(guān)心點的地理坐標，可計算出關(guān)心點的風(fēng)險值、個人風(fēng)險等值線以及表示社會風(fēng)險水平的F-N曲線。

源項數(shù)據(jù)：泄漏點離地面1 m高，垂直方向，井噴持續(xù)時間為15 min^[6]，釋放速率(最大絕對無阻流量)為1 260×10⁴ m³／d，H₂S含量為16％(體積分數(shù))，天然氣密度為0.9 kg／m³，溫度為25℃，井噴失控事故概率為4.5×10^-4／a。計算結(jié)果如圖6、7所示。

個人風(fēng)險等值線的變化與井場所在的地形、局地風(fēng)場密切相關(guān)。若分別按1×10^-4、1×10^-5、1×10^-6作為風(fēng)險呵接受標準進行評價，井口外個人風(fēng)險等值線最遠處分別為東南350 m、430 m和700 m，最近處距井口 l00-150 m，相應(yīng)的不規(guī)則區(qū)域范圍外人員風(fēng)險處于可接受水平。社會風(fēng)險F-N曲線說明井噴的風(fēng)險處于可接受區(qū)域，但還是應(yīng)采取技術(shù)合理可行、成本盡可能低的控制措施進一步降低風(fēng)險。

4 結(jié)論

1)基于計算流體力學(xué)(CFD)的有毒氣體擴散模擬方法，建立了適用于復(fù)雜地形條件下含硫天然氣井噴或地面設(shè)施泄漏所致的硫化氫擴散仿真模型，在此基礎(chǔ)上進行定量風(fēng)險計算，其個人風(fēng)險等值線能與地形、氣象條件很好地吻合。

2)事故概率通過借鑒并修正采用國外失效數(shù)據(jù)庫與西南油氣田公司同類設(shè)備歷史統(tǒng)計數(shù)據(jù)相結(jié)合的方式來確定，符合我國含硫天然氣開發(fā)技術(shù)現(xiàn)狀及社會、自然環(huán)境特征，據(jù)此完成的定量風(fēng)險計算結(jié)果能有效地為含硫天然氣開發(fā)場所風(fēng)險排序、風(fēng)險削減措施的取舍等提供更科學(xué)的依據(jù)。

參考文獻

[1] 王漢臣.大氣保護與能源利用[M].北京：中國環(huán)境科學(xué)出版社，1 992.

[2] UIJT de HAAG P A M，ALE B J M.CPRl8E Guideline for quantitative risk assessment[S].Hague：VROM，2005.

[3] APLAPI 581 Risk based inspection[S].Washington：American Petroleum Institute，2008.

[4] Center for Chemical Process Safety(CCPS).Guidelines for developing quantitative safety risk criteria[S].New York：John Wiley＆Sons Inc.，2009.

[5] 王福軍.計算流體動力學(xué)分析——CFD軟件原理與應(yīng)用[M].北京：清華大學(xué)出版社，2004.

[6] 國家安全生產(chǎn)監(jiān)督管理總局.AQ201 6含硫化氫天然氣井失控井口點火時間規(guī)定[S].北京：中國標準出版社，2008.

本文作者：翁幫華張林霞 胥云麗 劉坤 周東 李宇

作者單位：中國石油西南油氣田公司安全環(huán)保與技術(shù)監(jiān)督研究院