液氨泄漏事故擴散模擬

武偉男++楊平

摘 要：系統對比了高斯多煙團模式與SLAB模型模擬液氨儲罐泄漏后的氨氣擴散特征。結果表明，兩種模型的模擬結果存在較為明顯差異。在模擬設定條件下，事故發生點下風向60～2000 m范圍內，SLAB模型得到的最高濃度高于多煙團模式，前者是后者的1.01～35.2倍，且差別隨距離增大而增大。事故發生點下風向600 m以內，SLAB模型模擬得到的橫向影響距離大于多煙團模式；而在下風向600 m以外，多煙團模式模擬得到的橫向距離大于SLAB模型，差距隨下風向距離增加而增大。下風向同一地點，SLAB模型得到的氨氣最高濃度出現時間較多煙團模式較早，SLAB模型計算得到的氨氣煙團出現到消散時間也較多煙團模式更短。上述結果可為化學品泄漏導致突發環境事件的預防和應急中模型選擇提供參考。

關鍵詞：液氨 泄漏 擴散模擬 多煙團模型 SLAB模型

中圖分類號：X937 文獻標識碼：A 文章編號：1674-098X（2017）03（b）-0024-05

Diffusion Simulation of Liquid Ammonia Leakage

Comparison of the Multi-puff Model and SLAB Model

Wu Weinan1 Yang Ping2

（1.Solid waste Management Center in Liaoning Provine， Shenyang Liaoning， 110161， China； 2.Panjin Liaoning Fried Dough Sticks as for as sludge Treatment and Utillzation co.， LTD，Panjing Liaoing，124218，China）

Abstract：Simulation results of diffusion after liquid ammonia leakage calculated by the Gaussian multi-puff model and SLAB model were systematically compared. Results showed that there were obvious differences between the two models. Under the setting conditions， the round maximum ammonia concentrations simulated by the SLAB model were higher than those by the multi-puff model within 60 to 2000 m downstream the resource. And the former was 1.01 to 35.2 times that of the latter， and the difference increased with increasing distance. Higher cross-affected distances were found by SLAB model within 600 m downstream the resource， while cross-affected distances simulated by the multi-puff model were higher outside 600 m downstream， and the differences between the two models increases with the distances. In the same location downwind， the highest concentration of ammonia came earlier in SLAB model， while the time period from appearance and dissipation was shorter in multi-puff model. These results may provide a reference on diffusion model selection for prevention and response of environmental emergencies caused by chemical releases.

Key Words：Liquid ammonia；Leakage；Diffusion simulation；Multi-plume model；SLAB model

近年來，突發性環境事件頻發。以液氨等有毒氣體泄漏為代表的突發性環境事件往往導致嚴重后果，易形成大面積的危險區域，對周圍的環境和人員造成嚴重的危害。液氨是一種易燃易爆、有毒有害的化工原料，有腐蝕性并極易揮發。低濃度氨對粘膜有刺激作用，高濃度可造成組織溶解壞死[1]。氨氣泄漏和爆炸事故往往會導致眾多人員中毒或死亡，給公眾的生命健康和環境安全造成非常嚴重的影響。2013年6月3日，吉林省德惠市寶源豐禽業有限公司因氨氣泄漏爆炸，導致121人死亡，76人受傷，直接經濟損失1.82億元；2013年8月31日，上海市寶山區上海翁牌冷藏實業有限公司發生液氨泄漏，造成15人死亡，25人不同程度受傷[2]。因此，對有毒氣體發生泄漏后的擴散范圍、泄漏物質空氣中含量的時空分布、對人體造成危害的區域進行模擬預測和環境風險分析，對于突發環境事故預防和應急均具有重要意義。

目前，廣泛應用的氣體擴散模型包括高斯模型[3]、SLAB模型[4]、SUTTON模型[5]、ALOHA模型[6]等。國內外學者針對上述模型開展了應用性研究。莫秀忠等[7]基于MATLAB計算平臺高斯煙團模型建立了液氨泄漏后的濃度分布模型。鄒旭東等利用SLAB模型模擬了氯氣泄漏后氯氣擴散的時間、范圍和對周圍環境的危害[8]。王爽和王志榮以某化工廠的氯化氫泄漏事故為背景，利用ALOHA重氣擴散模型對該事故進行模擬，分析了敏感點濃度和人體接觸劑量隨時間的變化[6]。

已有研究主要集中于不同模型的應用、濃度模擬、影響條件分析、風險區域劃分等方面，缺少不同模型的結果的橫向對比。該研究以液氨泄漏事故為例，對比分析高斯多煙團模式和SLAB模型模擬結果的差異，以期為環境應急管理過程中擴散模型選擇提供參考。

1 模型及模擬條件

1.1 多煙團模式

多煙團模式基于高斯模型，是我國《環境風險評價技術導則》（HJ/T-2004）的推薦模型[3]，適用于瞬時泄漏擴散。該模式把風險源煙團輸送時間分割為若干時段，假定每個時段發射一個煙團，計算每個煙團在各時刻對關心點的貢獻[9]。

第i個煙團在時刻、在點（x，y，0）產生的濃度為：

（1）

式中，為煙團排放量，mg，；為釋放率，mg/s；為時段長度，s；，，為煙團在w時段沿x、y、z方向的等效擴散參數，m，可按照式（2）、（3）估算；、分別為第w時段結束時第i煙團質心的x、y坐標，按照式（4）、（5）計算。

（j = x， y， z） （2）

（3）

（4）

（5）

各煙團對某關心點t時刻的貢獻濃度，按照式（6）計算。

（6）

式中，n為需要跟蹤的煙團數，由式（7）確定。

（7）

式中，f為小于1的系數，根據計算要求確定。

1.2 SLAB模型

SLAB模型由美國能源部的Lawrence Livermore國家實驗室開發，用于比空氣稠密氣體泄放的大氣擴散模擬[4]，是美國EPA推薦危險化學品意外泄放事故模擬的應急模型。SLAB模型的模擬源可以是持續的，有限持續時間的，或者瞬間泄放。持續和有限持續時間泄放應用于蒸發池、水平射流和垂直射流泄放源。瞬間泄放則假設為瞬間體源進行模擬。模型可以處理的泄放類型包括：地平面蒸發池、有高度的水平射流或垂直射流、瞬間泄放以及液體溢漏。SLAB模型可用于稠密氣體釋放或者液體溢漏而蒸發出的稠密氣體擴散。盡管氨氣密度低于空氣，但由于氨氣經常以液氨形式儲存，液氨泄漏后因氣化時大量吸熱而具有重氣體的特點，屬于SLAB適用類型的后者。

SLAB模型假設事件發生在沒有障礙物的平坦區域，模型沒有考慮有坡度的地形條件。泄放物質的大氣擴散由守恒方程來計算，包括質量、動量、能量和組分守恒。持續泄放作為穩態煙羽處理。有限持續時間泄放的起始煙云擴散用穩態煙羽模式來解釋，一直持續到泄放源停止泄放。當泄放源被切斷時，煙云作為煙團處理，其隨后的擴散使用瞬變煙團模式進行計算。在預測濃度隨時間變化方面，SLAB模型在穩定、中度穩定及不穩定的大氣環境下均能得到較好的預測結果[8， 10]。

1.3 模擬條件

研究模擬一存有6 000 kg液氨儲罐泄漏事故后的液氨擴散情景。儲罐內壓力為250 kPa，裂口面積為0.0004 m2，裂口之上液位高度為2 m，持續時間為10 min，環境溫度為25℃。經計算得到泄漏持續時間內，泄漏量為4.43 kg/s。假設事故發生時大氣穩定度為D，風速為2.0 m/s。對于多煙團模式，假定10 s一個煙團；對于SLAB模型，選用水平射流模式。模擬時間為20 min，模擬范圍為事故源下風向縱向2 000 m，橫向1000 m。

2 結果與討論

設定條件下，多煙團模式和SLAB模型得到的模擬結果分別如圖1、圖2所示。2中模型都能得到氨氣地面濃度隨時間、位置的變化。可以看出，多煙團模式模擬煙團移動較SLAB模型更快。本設定液氨的泄漏時間為10 min，多煙團模式認為泄漏停止后，氨氣煙團將立即離開事故發生點，向下風向遷移；而SLAB模型假設液氨泄漏至地面，形成液池，盡管在液氨儲罐停止泄漏后，液池中的液氨將繼續揮發并持續一段時間。因此，SLAB模型模擬事故發生13 min后，仍有氨氣從事故發生點揮發。這種假設上的差異也導致了SLAB模型模擬的煙團跨度較多煙團模式更大。

2種模型得到的氨氣軸線（沿x軸方向）地面最大濃度與不同下風向距離關系如圖3所示。SLAB模型結果表明距儲罐1～2 m范圍內的最高氨氣體積百分比達到100%，而多煙團模式結果表明距儲罐13～14 m最高氨氣體積百分比均到100%。60～2 000 m范圍，SLAB模型得到的最高濃度是多煙團模式得到的最高濃度的1.01～35.2倍；且隨著下風向距離增加，SLAB與多煙團模式最高濃度差距加大。若以最高氨氣濃度達到1 390 mg/m3（半致死濃度）的地點為疏散區域，多煙團模式模擬得到的疏散半徑約為360 m，而SLAB模型模擬得到的疏散半徑約為550 m。

在與下風向垂直的橫向（y軸方向）方面，分析下風向縱向2 000 m、橫向1 000 m范圍內（即x=0～2 000 m，y=-1 000～1 000 m），0～20 min時刻內，氨氣在各點的最高濃度，下風向各點橫向位置濃度與軸線濃度之比超過10%的范圍定義為橫向影響距離。橫向影響距離與下風向距離關系如圖4所示。在設定的模擬條件下，下風向600 m以內，SLAB模型模擬得到的橫向影響距離大于多煙團模式；而在下風向600 m以外，多煙團模式模擬得到的橫向距離大于SLAB模型，并隨下風向距離增加，差距增大。從前述分析可知，2種模型得到的疏散半徑均在600 m以內，因而考慮到橫向的影響范圍，SLAB模型模擬得到的疏散面積也大于多煙團模式，分別約為6.63萬m2和1.72萬m2。

分析下風向某處在事故發生后氨氣濃度隨時間變化，多煙團模式和SLAB模型模擬結果如圖5所示。可以看出，越靠近事故源的地點2種模型得到的模擬結果相似度越高。下風向100 m以后，2種模型結果差異逐漸明顯。第一，如前所述，SLAB模型給出的最高濃度高于多煙團模式；第二，SLAB模型給出的氨氣最高濃度出現時間較高斯模式較早，且差異隨下風向距離增加而增大；第三，同一地點SLAB模型計算得到的氨氣煙團出現到消散時間也較多煙團模式更短。

孫召賓[11]使用Burro現場實驗數據，對國內多煙團模式、SLAB模型和ALOHA模型計算的液化氣泄漏模擬結果可靠性進行了驗證，結果表明，高斯模型較不適用于稠重氣云擴散的數值模擬，SLAB模型、ALOHA模型和高斯模型模擬結果的可靠性排序為SLAB模型>ALOHA模型>高斯模型。高凌[12]將SLAB模型與液氯模擬泄漏試驗結果進行了對比，發現事故下風向300～1 000 m范圍內，SLAB模型模擬結果是實測結果的約4～6倍。國內尚缺少針對液氨泄漏的模型模擬與實測結果的對比研究，因而無法判斷“導則”推薦的多煙團模式和SLAB模型更為準確。該研究結果表明液氨泄漏后，SLAB模型模擬的地面氨氣濃度高于多煙團模式，得到的致死區域也高于多煙團模式。

3 結語

該研究系統對比了多煙團模式與SLAB模型模擬液氨儲罐泄漏后的氨氣擴散。結果表明，兩種模型均能模擬液氨泄漏后地面氨氣濃度的時空分布，但兩種模型的模擬結果存在明顯差異。在模擬設定條件下，事故發生點下風向60～2 000 m范圍內，SLAB模型得到的最高濃度是多煙團模式得到的最高濃度的1.01～35.2倍。事故發生點下風向600 m以內，SLAB模型模擬得到的橫向影響距離大于多煙團模式；而在下風向600 m以外，多煙團模式模擬得到的橫向距離大于SLAB模型，并隨下風向距離增加，差距增大。以最高氨氣濃度達到半致死濃度為疏散區域，SLAB模型模擬得到的疏散半徑和疏散范圍分別是多煙團模式的1.5和3.9倍。下風向同一地點，SLAB模型得到的氨氣最高濃度出現時間較多煙團模式較早，SLAB模型計算得到的氨氣煙團出現到消散時間也較多煙團模式更短。

化學品泄漏等環境突發事件的風險防控越來越受到關注，我國多地構建了基于擴散模型的環境風險應急支持系統。該研究結果表明，不同模型得到的模擬結果存在顯著差異，可為液氨泄漏突發環境事件預防和應急中模型選擇提供參考。模型給出的濃度過高，將增大救援、疏散的范圍，可能會影響到安全區域內人們的正常生產生活，造成一些不必要的浪費；而模型給出的濃度過低，則將使受影響人群不能及時疏散，造成人員傷亡。因此，在未來，針對環境突發事件發生頻率較高的化學品的擴散開展試驗研究，對模型進行篩選和優化，才能更有效地指導環境風險防控和應急響應工作。

參考文獻

[1] 張杰，趙明.液氨泄漏事故的定量風險評價研究[J].安全與環境工程，2012，19（1）：69-72.

[2] 夏登友，錢新明，黃金印，等.液氨泄漏擴散模擬及危害評估[J].中國安全科學學報，2014，24（3）：22-27.

[3] HJ/T 169-2004，建設項目環境風險評價技術導則[S].

[4] Donald L. Ermak. User's Manual For SLAB：An Atmospheric Dispersion Model For Denser-Than-Air Releases[R].Lawrence Livermore National Laboratory，1990.

[5] 孫莉， 趙穎，曹飛，等.危險化學品泄漏擴散模型的研究現狀分析與比較[J].中國安全科學學報，2011，21（1）：37-42.

[6] 王爽，王志榮.利用ALOHA軟件對一起氯化氫泄漏事故的模擬分析[J].滅火指揮與救援，2010，29（8）：698-700.

[7] 莫秀忠，吳欣甜，謝飛，等.基于MATLAB的液氨瞬時泄漏模擬及應急措施研究[J].南開大學學報：自然科學版，2014，47（4）：1-5.

[8] 鄒旭東，楊洪斌，汪宏，等.SLAB在突發大氣污染事件應急模擬中的應用[J].環境科學與技術，2010，33（12F）： 588-590.

[9] 胡二邦.環境風險評價實用技術、方法和案例[M].北京：中國環境科學出版社，2009.

[10] 瞿子晶，鐘圣俊.WebGIS和SLAB模型的突發性大氣污染事故模擬和應用[J].環境科學與技術，2014，37（5）： 107-111.

[11] 孫召賓.危險化學品泄放事故后果計算模型的研究及應用[D].大連：大連理工大學，2012.

[12] 高凌.SLAB View在化學泄漏事故應急救援中的應用[J].消防科學與技術，2011，30（9）：833-836.