官學(xué)源

中國石油天然氣管道工程有限公司沈陽(yáng)分公司

摘要：針對埋地含硫天然氣管道泄漏問(wèn)題，采用有限體積法，建立泄漏擴散模型和燃燒模型，考慮土壤對泄漏氣體的影響，在風(fēng)場(chǎng)影響下，對比分析埋地含硫天然氣管道泄漏時(shí)有毒氣體的擴散情況和發(fā)生燃燒后有毒氣體的擴散情況，得出泄漏氣體燃燒擴散比單一的氣體泄漏擴散更加危險，燃燒產(chǎn)物二氧化硫為主要危害氣體，同時(shí)，存在大量的二氧化碳和水蒸氣以及未發(fā)生反應的硫化氫，危害性極大，需及時(shí)采取防范和急救措施。

關(guān)鍵詞：含硫天然氣管道；燃燒；數值模擬；泄漏；多孔介質(zhì)

天然氣管道泄漏事故不僅影響管道安全生產(chǎn)，而且其易燃易爆的危險也會(huì )給人們生活和環(huán)境構成威脅。針對管道泄漏擴散問(wèn)題，我國學(xué)者進(jìn)行了大量的數值研究[1-5]，而針對管道泄漏著(zhù)火的研究相對較少。本文通過(guò)FLUENT軟件，基于有限體積法，建立埋地含硫天然氣管道持續泄漏擴散和燃燒模型，對比分析泄漏時(shí)有毒氣體的擴散情況和發(fā)生燃燒后有毒氣體的擴散情況。為埋地含硫天然氣管道的風(fēng)險評估、安全運行和應急管理提供參考。

1  數學(xué)模型

1.1  組分輸運守恒方程 [式（1）]

式中：Ri為第 i 組分化學(xué)反應速率，kgmol/（m3·s）；Si為擴散相加用戶(hù)定義的源項；Jj為組分 j 的擴散通量，mol/（m2·s）；Yi為第 i 種物質(zhì)的質(zhì)量分數，% ；u為流體流速，m/s；ρ為密度，kg/m3；t為時(shí)間，s 。

湍流中的質(zhì)量擴散方程 [式（2）]

式中：Sct為湍流施密特數；Di,m為混合物中第 i 種物質(zhì)的擴散系數，μt為湍流黏度，Pa·s。

1.2  渦耗散模型 [式（3）]

式中：Ri,r為反應 r 中物質(zhì) i 的產(chǎn)生/分解速度，kgmol/（m3·s）；Mw,i為第 i 種物質(zhì)的質(zhì)量分數，%；ε/k 為大尺度渦混合時(shí)間尺度，s；N為系統化學(xué)物質(zhì)數目；YP，YR為某種產(chǎn)物的質(zhì)量組分；A，B為經(jīng)驗常數，這里分別取0.4和0.5； ν'i,r為反應 r 中反應物 i 的化學(xué)計量系數。

1.3  完全燃燒的化學(xué)反應：

2  數值模擬與結果分析

2.1  問(wèn)題描述

某埋地含硫天然氣管道，取計算管長(cháng)間距為5 km， 管道外徑650 mm，埋深1.5 m。管道上層土壤為多孔介質(zhì)，孔隙度為0.267，密度為2650 kg/m3，導熱系數1.512 W/m·K。高含硫天然氣持續泄漏著(zhù)火的模擬區域范圍為5000×10011.5 m：泄漏口中心在x=0、y=﹣1.5處，泄漏口直徑為0.05 m，泄漏速度為252 m/s，泄漏方向垂直向上，火源在泄漏口附近。天然氣中甲烷（CH4）體積分數為95%，硫化氫（H2S）體積分數為5%；空氣組分中甲烷的爆炸體積濃度區間為5%～15%，硫化氫的爆炸體積濃度區間為4.3%～46.0%；硫化氫的中毒體積濃度下線(xiàn)為0.001%，二氧化硫（SO2）的中毒體積濃度下限為0.04%；泄漏氣體溫度相同為300 K，環(huán)境溫度與天然氣溫度相同；大氣壓力為101325 Pa；風(fēng)速為4 m/s；空氣中氮氣組分占78%，氧氣占22%。

2.2  FLUENT軟件模擬結果

采用二維單精度SILPLE方法進(jìn)行隱式耦合求解，采用矩形網(wǎng)格在泄漏口處網(wǎng)格加密，模擬結果見(jiàn)圖 1～5。

圖 1 30 min內硫化氫泄漏擴散圖示

圖 2 泄漏趨于穩定時(shí)硫化氫擴散云圖

圖 3 燃燒云團

圖 4 燃燒產(chǎn)物擴散云圖

圖 5 未燃燒的硫化氫擴散云圖

2.3  結果分析

埋地含硫天然氣管道泄漏初期，由于管道上層多孔介質(zhì)對泄漏氣體的影響，使其不能形成高速射流。氣體進(jìn)入土壤后，土壤的毛管壓力、土壤阻力和地表的張力對泄漏氣體的共同作用，使泄漏氣體的湍能大量減少，但持續泄漏使氣體的湍能不斷增大，仍會(huì )通過(guò)土壤迅速且持續不斷地涌出地表，進(jìn)入大氣進(jìn)行射流運動(dòng)和擴散運動(dòng)。土壤作為多孔介質(zhì)，具有一定的孔隙度和獨特的儲存能力。一定面積的土壤飽和后，氣體將會(huì )繼續向土壤橫向擴散，造成更大面積的污染，同時(shí)縱向方面氣體將滲出土壤外地表，消耗氣體動(dòng)能。因此，在短時(shí)間內有大量的高濃度含硫天然氣聚集在泄漏口附近的地表處。隨著(zhù)時(shí)間的推移大量泄漏氣體進(jìn)入環(huán)境中，其主要有毒危害氣體為硫化氫氣體。在低空處，泄漏2 min時(shí)，硫化氫有毒云團距離泄漏口下方向475 m；泄漏8 min時(shí)，距離1637 m；泄漏10 min～30 min內，距離近似呈線(xiàn)性增長(cháng)，在30 min時(shí)，距離達4045 m（ 圖 1）。泄漏趨于穩定后，硫化氫有毒云團距離泄漏口下方向超過(guò)4800 m，高為248 m（圖 2）。

當泄漏氣體涌出地表進(jìn)入大氣，在泄漏口附近遇到火源時(shí)，泄漏氣體組分中的甲烷和硫化氫遇空氣中的氧氣迅速發(fā)生劇烈化學(xué)反應且放出大量熱量，因為風(fēng)場(chǎng)的影響，整個(gè)燃燒區域在泄漏口下方向，最高溫度達2158 K（圖 3）。燃燒反應生成大量的二氧化硫、二氧化碳和水蒸汽。其中二氧化硫有毒濃度云團水平擴散距離超過(guò)4800 m，最大體積濃度高達12.6%，遠高于其中毒濃度下線(xiàn)；水蒸汽和二氧化碳最大體積濃度分別超過(guò)13%和12%，其水平擴散距離超過(guò)2000 m，雖然水蒸汽和二氧化碳氣體無(wú)毒，但是其滯留在空氣中可能發(fā)生化學(xué)反應生成酸，不僅污染大氣，當遇到金屬、建筑物時(shí)，可能對其造成腐蝕破壞（圖 4）。而在泄漏口下方向0～1500 m范圍內聚集著(zhù)大量高濃度的未反應的硫化氫氣體，最大濃度為50.7%（圖 5）。對比二氧化硫氣體和硫化氫氣體的擴散情況可知，二氧化硫危害范圍遠高于硫化氫氣體危害區域，并且完全覆蓋其危險范圍。

3   結束語(yǔ)

對比兩種擴散情況，單一的泄漏擴散主要有毒成分為硫化氫氣體；而發(fā)生燃燒時(shí)情況相對復雜且危害程度更高，主要有毒成分為二氧化硫氣體，同時(shí)還存在高濃度的硫化氫有毒氣體和潛在的危害組分水蒸汽和二氧化碳。應對不同情況采取相應的防范和急救措施。

參考文獻：

[1]程浩力，劉德俊，劉倩倩，等.城燃管道街道峽谷泄漏擴散CFD數值模擬[J].石油化工高等學(xué)校學(xué)報，2011（4）：60-63.

[2]程浩力，劉德俊.城鎮燃氣管道泄漏擴散模型及數值模擬[J].遼寧石油化工大學(xué)學(xué)報，2011，2（31）：27-31.

[3]李朝陽(yáng)，馬貴陽(yáng).高含硫天然氣管道泄漏數值模擬[J].化學(xué)工 程，2011，39（7）：88-92.

[4]李振林，姚孝庭，張永學(xué)，等.基于FLUENT的高含硫天然氣管道泄漏擴散模擬[J].油氣儲運，2008，27（5）：38-41.

[5]付吉強，熊新強，王蕓,等.架空含硫天然氣管道泄漏擴散數值模擬研究[J]．當代化工，2011，40（10）：1079-1083．

作者簡(jiǎn)介： 官學(xué)源 ，工程師，碩士， 2013年畢業(yè)于遼寧石油化工大學(xué) 。聯(lián)系方式： 18040014317，guanxueyuan00@126.com。