李旭東，張園園，張巨偉，奚譽(yù)寧，韓明桓

（遼寧石油化工大學(xué)環(huán)境與安全工程學(xué)院，遼寧撫順 113001）

硫化氫是一種重要的化工物料，具有易燃易爆、有毒等特點(diǎn)。管道閥門(mén)或法蘭密封不良、儲(chǔ)罐溫度或壓力過(guò)大、液位過(guò)高等因素都可能導(dǎo)致硫化氫泄漏，遇火源可能發(fā)生爆炸；若未遇到火源則可能在空氣中散逸致人員中毒[1]。不同劑量硫化氫導(dǎo)致人體損害程度不同。分析硫化氫擴(kuò)散的影響因素，特別是分析隨著不同影響因素量的變化導(dǎo)致硫化氫擴(kuò)散后果的變化情況，對(duì)于制定應(yīng)急救援預(yù)案，安排人員疏散具有重要意義。

1 理論簡(jiǎn)介

硫化氫擴(kuò)散后果以湍流擴(kuò)散微分方程為基礎(chǔ)進(jìn)行分析[2]。根據(jù)泄漏原因有連續(xù)泄漏和瞬時(shí)泄漏，擴(kuò)散模型根據(jù)泄漏類(lèi)型分為煙團(tuán)模型和煙羽模型。

1.1 瞬時(shí)泄漏分析

如果石油化工裝置爆破片破裂、安全閥錯(cuò)誤開(kāi)啟或者管道發(fā)生斷裂、石油化工裝置發(fā)生沖料、放空閥門(mén)異常開(kāi)啟等導(dǎo)致大量介質(zhì)瞬間泄漏，擴(kuò)散后果用煙團(tuán)模型來(lái)描述，煙團(tuán)模型分為無(wú)風(fēng)煙團(tuán)模型和有風(fēng)煙團(tuán)模型[3]。

1.1.1 無(wú)風(fēng)情況

無(wú)風(fēng)情況下泄漏擴(kuò)散后果模型[4-5]見(jiàn)式（1）。

式中：C為氣云濃度，kg/m3；Q為泄漏源強(qiáng)，kg/s；Kx、Ky、Kz為x、y、z湍流擴(kuò)散系數(shù)；x，y，z為x，y，z軸；t為時(shí)間，s。

1.1.2 有風(fēng)情況

有風(fēng)情況下泄漏擴(kuò)散后果模型[6-8]見(jiàn)式（2）。

式中：u為環(huán)境風(fēng)速，m/s。

1.2 連續(xù)泄漏分析

連續(xù)泄漏一般發(fā)生在石油化工裝置出現(xiàn)穿孔、裂紋等，導(dǎo)致介質(zhì)長(zhǎng)時(shí)間緩慢持續(xù)泄漏。介質(zhì)擴(kuò)散過(guò)程類(lèi)似羽毛一樣緩慢飄散稱(chēng)作煙羽模型[9-11]。煙羽模型需要考慮無(wú)風(fēng)和有風(fēng)兩種狀況。

1.2.1 無(wú)風(fēng)情況

無(wú)風(fēng)情況下泄漏擴(kuò)散后果模型[12-13]具體見(jiàn)式（3）。

1.2.2 有風(fēng)情況

有風(fēng)情況下泄漏擴(kuò)散后果模型[14-16]具體見(jiàn)式（4）。

2 MATLAB 數(shù)值模擬及影響因素分析

2.1 數(shù)值模擬條件

硫化氫的擴(kuò)散受到風(fēng)速、大氣穩(wěn)定度、源高、源強(qiáng)等因素影響，該文采用MATLAB軟件模擬多因素對(duì)硫化氫擴(kuò)散后果的影響。模擬過(guò)程中只讓某一種因素的量值發(fā)生變化，其它因素保持固定，分析擴(kuò)散后果的變化情況。模擬的參數(shù)為：地形平坦，有風(fēng)，硫化氫連續(xù)泄漏，大氣穩(wěn)定度為中性。

通過(guò)模擬軟件對(duì)有風(fēng)連續(xù)煙羽模型下的6 種不同情況進(jìn)行擴(kuò)散后果模擬。具體模擬參數(shù)如表1 所示，模擬結(jié)果如圖1所示。

表1 擴(kuò)散后果模擬參數(shù)

2.2 模擬結(jié)果分析

條件（1）模擬結(jié)果如圖1（a）所示。在橫軸x等于41 m 處硫化氫達(dá)到峰值濃度，峰值濃度為374.26 mg/m3。條件（2）源強(qiáng)變?yōu)樵瓉?lái)的2倍，模擬結(jié)果如圖1（2）所示，峰值濃度變?yōu)?48.52 mg/m3為原來(lái)2倍，仍出現(xiàn)在x等于41 m處。

條件（3）源強(qiáng)變?yōu)樵瓉?lái)的4倍，硫化氫峰值濃度模擬結(jié)果為1 497.04 mg/m3（模擬圖省略），峰值濃度亦變?yōu)樵瓉?lái)4 倍。說(shuō)明擴(kuò)散峰值濃度會(huì)隨著泄漏源強(qiáng)的增加呈等倍數(shù)變化，但峰值出現(xiàn)位置不變。

條件（3）風(fēng)速變?yōu)?倍。模擬結(jié)果如圖1（3）所示，硫化氫峰值濃度變?yōu)?87.13 mg/m3，為原來(lái)1/2，仍出現(xiàn)在x等于41 m處；當(dāng)模擬條件為：源高H＝8 m，源強(qiáng)Q＝0.35 kg/s，風(fēng)速u(mài)＝12 m/s，風(fēng)速變?yōu)闂l件（1）4倍，模擬結(jié)果硫化氫峰值濃度為93.56 mg/m3（模擬圖省略），峰值濃度變?yōu)樵瓉?lái)的1/4倍。說(shuō)明泄漏峰值濃度會(huì)隨著風(fēng)速的增加成相反倍數(shù)的變化，但是峰值出現(xiàn)在x軸的位置不變。風(fēng)速對(duì)峰值濃度的出現(xiàn)具有稀釋作用。

條件（4）和（1）相比，源高變?yōu)樵瓉?lái)的2倍。模擬結(jié)果如圖1（4）所示，硫化氫峰值濃度由374.26 mg/m3改變?yōu)?3.76 mg/m3，峰值濃度變?yōu)樵瓉?lái)的1/4倍，峰值濃度出現(xiàn)在x等于82 m處；

條件（5）和（1）相比，源高為原來(lái)4 倍，模擬結(jié)果如圖1（5）所示，硫化氫峰值濃度為23.53 mg/m3，為原來(lái)的1/16 倍，峰值濃度出現(xiàn)在x等于166 m處；

條件（6）和（1）相比，源高為原來(lái)8 倍，模擬結(jié)果如圖1（6）所示，硫化氫峰值濃度為5.92 mg/m3，為原來(lái)的1/63，峰值濃度出現(xiàn)在x等于341 m 處；可看出源高對(duì)擴(kuò)散區(qū)域峰值影響較大，隨著源高增高，擴(kuò)散區(qū)域峰值濃度下降較快。隨著源高的變化，擴(kuò)散峰值所在x軸的位置基本與源高呈等比例變化。

圖1 6種情況模擬結(jié)果

3 模擬結(jié)果的應(yīng)用分析

硫化氫泄漏后最關(guān)鍵的問(wèn)題是組織人員疏散。以上模擬結(jié)果對(duì)于組織人員疏散和制定應(yīng)急救援預(yù)案具有重要的指導(dǎo)意義。

根據(jù)模擬結(jié)果可知：源強(qiáng)和風(fēng)速的變化不會(huì)影響硫化氫擴(kuò)散后峰值濃度出現(xiàn)的位置，但是源高的變化會(huì)影響硫化氫擴(kuò)散后峰值濃度出現(xiàn)的位置，所以管理人員在應(yīng)急救援疏散的時(shí)候可以根據(jù)以上規(guī)律，避開(kāi)硫化氫峰值濃度出現(xiàn)的區(qū)域，組織人員向遠(yuǎn)離峰值濃度出現(xiàn)的區(qū)域疏散。

4 結(jié)論

（1）該文分析了硫化氫連續(xù)泄漏過(guò)程中，源強(qiáng)、平均風(fēng)速、泄漏源高度呈倍數(shù)變化后擴(kuò)散區(qū)域峰值的變化的結(jié)果，擴(kuò)散區(qū)域峰值所處x位置的變化結(jié)果。

（2）模擬結(jié)果表明，源強(qiáng)、平均風(fēng)速、泄漏源高度的變化都會(huì)影響擴(kuò)散區(qū)域硫化氫峰值濃度。源強(qiáng)的增加會(huì)增加擴(kuò)散區(qū)域硫化氫峰值的濃度，泄漏峰值濃度會(huì)隨著泄漏源強(qiáng)的增加等倍數(shù)的變化，但不會(huì)改變硫化氫峰值濃度出現(xiàn)的區(qū)域。平均風(fēng)速的增加會(huì)稀釋擴(kuò)散區(qū)域硫化氫峰值的濃度，泄漏峰值濃度會(huì)隨著平均風(fēng)速的增加呈相反倍數(shù)的變化，不會(huì)改變硫化氫峰值濃度出現(xiàn)的區(qū)域。隨著泄漏源高度的增加，擴(kuò)散區(qū)域硫化氫峰值濃度顯著下降。

（3）石油化工企業(yè)硫化氫中毒事故造成嚴(yán)重后果，以上研究結(jié)果對(duì)企業(yè)安全預(yù)防具有重要意義。