高瓦斯低透氣性煤層CO2相變致裂增透技術(shù)研究*

洪紫杰，王　成,2，熊祖強,2

(1.河南理工大學 能源科學與工程學院，河南 焦作 454003；2.煤炭安全生產(chǎn)河南省協(xié)同創(chuàng)新中心，河南 焦作 454003)

0　引言

我國90%以上的礦井屬于高瓦斯低透氣性礦井，瓦斯治理問題日益突出。面對瓦斯抽采濃度低、衰減速率快等難題，國內(nèi)外學者提出多種抽采治理瓦斯的措施，其中增透技術(shù)受到煤炭行業(yè)普遍青睞[1-3]。為致力于提高瓦斯抽采效果，我國煤炭行業(yè)研究者提出了開采保護層、水力掏槽[4-5]、深孔預裂爆破、水力沖孔、水力壓裂、水力割縫等[6-7]多種增透措施，雖然上述技術(shù)措施取得了一定效果，但由于工藝及技術(shù)措施不完善，沒有取得預期的理想效果[8]。上世紀伴隨著CO2相變致裂技術(shù)的迅速發(fā)展，在煤炭企業(yè)中取得了顯著成績，CO2相變致裂技術(shù)克服了炸藥爆破增透方法的高危險性及破壞性[9-11]，為煤層增透、工作面放頂及巷道開挖提供了新思路。但是CO2相變致裂機理、致裂能量分析及關(guān)鍵部件釋放管幾何參數(shù)的研究仍停留在發(fā)展的初期階段[12-15]，尤其對釋放管幾何參數(shù)研究更是缺少系統(tǒng)的理論與實踐，以上原因使CO2相變致裂技術(shù)處于不成熟階段。

為了更加科學合理有效的利用CO2相變致裂技術(shù)增透促抽瓦斯，通過理論分析計算，F(xiàn)luent數(shù)值模擬[16]等技術(shù)措施，詳細分析了CO2致裂原理，驗算了CO2致裂能量，選取了釋放管幾何最優(yōu)參數(shù)[17-18]，為進一步推廣致裂技術(shù)和提高瓦斯抽采效果提供依據(jù)。

1　CO2相變致裂增透技術(shù)原理及煤層損傷分析

1.1　CO2相變致裂原理

CO2致裂技術(shù)原理是預先在瓦斯抽放孔內(nèi)裝入CO2儲液管，通過接通裝置電流啟動儲液管發(fā)熱片，致使管內(nèi)CO2在20 ms內(nèi)受熱體積瞬間膨脹600倍以上；儲液管內(nèi)CO2也集聚吸熱發(fā)生汽化反應，由液態(tài)轉(zhuǎn)化為氣態(tài)，此時儲液管氣態(tài)CO2壓力達到120～276 MPa；氣態(tài)CO2在如此高壓下直接沖破儲液管中的定壓剪切片裝置，由釋放管排放孔瞬間噴發(fā)。這種充分利用氣體高壓沖擊、破碎煤巖體的相變致裂技術(shù)，進一步實現(xiàn)了煤巖體的致裂效果，增加了煤巖層的透氣性。

1.2　致裂后煤體損傷方程模型分析

CO2致裂是高壓氣體瞬間沖擊、破碎煤體，應力-應變值達到損傷閾值時，煤體微觀單元即會發(fā)生破壞損傷，煤體由于受到高壓氣體作用破碎變形，同時伴隨大量裂隙的產(chǎn)生，煤層透氣性也明顯增大，損傷單元彈性模量[19]如下式所示：

E=(1-D)E0

(1)

式中:E為損傷單元彈性模量， MPa；E0為無損傷單元彈性模量 MPa；D為損傷系數(shù)。

煤體微觀單元透氣系數(shù)損傷耦合方程遵循2種破壞規(guī)律，若煤巖體微觀單元拉伸應力超過抗拉強度σb損傷閾值，拉應力損傷系數(shù)D表示如下：

σt<σb

(2)

(3)

式中：εt1為抗拉極限應變；εt0為彈性拉伸極限應變；σbr為抗拉殘余強度， MPa。

對應的單元透氣性系數(shù)為：

(4)

式中：θ為瓦斯壓力系數(shù)；β為應力影響系數(shù)；ξ為透氣系數(shù)的增大系數(shù)；ξ0為單元損傷后透氣系數(shù)的增大系數(shù)；λ0為初始透氣系數(shù)；p0為瓦斯壓力， MPa；

若煤巖體微觀單元壓應力超過損傷閾值，損傷系數(shù)D按下式表示:

(5)

式中：σcr為抗壓殘余強度；εr為殘余應變；σc0為最大壓應變；

對應的單元透氣系數(shù):

(6)

式中：σc為壓應力， MPa。

煤體損傷后裂紋的擴張可用有效應力強度因子衡量，表達式如下：

(7)

式中：δ0為應力強度因子；δ1為有效應力強度因子；D為損傷系數(shù)。

由(7)可知，致裂作用時，煤體損傷系數(shù)越大，有效應力強度因子也越大，增透致裂效果越好。

2　CO2致裂當量計算分析

選用了常用的SD390型定壓泄能片、儲液管及相應的配套設(shè)置為例[20]，進行計算CO2相變致裂的TNT當量。參考與CO2儲罐爆炸相應的研究文獻，CO2儲罐爆炸當量計算方法選用壓縮氣體與飽和水蒸氣容器爆破模型的能量計算方法。公式如下：

(8)

式中:Eg為氣體爆炸能量，kJ；v為容器容積，2.16 m3；p為容器內(nèi)氣體絕對壓力，276 MPa；k為氣體的絕熱指數(shù)，取1.29；p1為大氣壓。

液態(tài)CO2相變致裂裝置的近似TNT當量WTNT參考下式計算：

(9)

式中：ETNT為1 kgTNT爆炸能量，4 250 kJ/kg。

由式(8)和式(9)可得,高壓CO2氣體爆炸能量約為1 688 kJ,當量WTNT約為398 g，因此，1 kg液態(tài)CO2相變致裂時，其爆炸的當量相當于398 gTNT炸藥，其致裂能量足以滿足現(xiàn)場對煤層致裂效果的要求。

3　CO2釋放管幾何參數(shù)優(yōu)化研究

由CO2相變致裂的基本原理可知，釋放管是整個相變致裂裝置的終端元件，同時也是形成高壓CO2氣體最關(guān)鍵部分，因此釋放管整個結(jié)構(gòu)及對應的幾何參數(shù)將會從根本上影響煤層的致裂效果。在整個致裂設(shè)備工作過程中，CO2集聚吸熱發(fā)生汽化反應后，通過釋放管直接作用于煤巖體。鑒于此，優(yōu)化研究CO2釋放管的幾何參數(shù)，對于增強其致裂效果有著至關(guān)重要的意義。論文選用Fluent數(shù)值模擬方法優(yōu)化釋放管的幾何參數(shù)。

圖2　4種結(jié)構(gòu)噴嘴動壓力Fig.2　Four nozzle structure dynamic stress nephogram

3.1　釋放管噴嘴的形式及幾何參數(shù)

針對目前釋放管噴嘴的研究情況，釋放管噴嘴主要分為以下4種結(jié)構(gòu)形式，見圖1。

圖1　FLUENT 釋放管噴嘴類型Fig.1　Fluent software release pipe nozzle type

釋放管噴嘴的幾何參數(shù)[6]主要包括：入口直徑D1，出口直徑D2，噴嘴長度L，錐形噴嘴和錐直型噴嘴收縮角θ。

3.2　釋放管噴嘴模擬結(jié)果及分析

3.2.1噴嘴幾何類型模擬研究與分析

根據(jù)理論分析，CO2氣體在高壓力條件下噴出時，其對煤巖體致裂效果相對較好，因此首先選擇現(xiàn)有的最高壓力276 MPa進行數(shù)值模擬分析，最后再進行驗證不同壓力的致裂效果。

由圖2對4種結(jié)構(gòu)的壓力云圖分析可知，通過模擬CO2氣體在高壓276 MPa條件下噴出的內(nèi)、外部流場得知，不同的噴嘴結(jié)構(gòu)，CO2氣體由噴嘴噴出后的動壓力分布差別很大。綜合比較，(b)，(c)，(d) 3種結(jié)構(gòu)在噴嘴處所產(chǎn)生的噴射流場比較細，(a)結(jié)構(gòu)在噴嘴處所產(chǎn)生的噴射流場相對較大，相當于后3種結(jié)構(gòu)的2～4倍；(a)結(jié)構(gòu)射流段致裂長度明顯長于后3種結(jié)構(gòu)，但最大動壓力量綱并無很大差異；另外，由(a)結(jié)構(gòu)噴出的CO2氣體擴散范圍也比較大，而且處于較高動壓力狀態(tài)下，有利于煤層內(nèi)裂隙的擴展，能進一步增透煤層的透氣性。鑒于此，(a)結(jié)構(gòu)空心圓柱結(jié)構(gòu)噴嘴性能效果顯著，有利于提高瓦斯抽放效果。

3.2.2噴嘴直徑模擬研究與分析

通過模擬噴嘴長度對致裂效果的影響，得出在其他條件相同的條件下，長度在18 mm，20 mm，24 mm時噴嘴長度對射流段致裂長度及煤層內(nèi)裂隙的擴展影響不大，綜合考慮裝置的便捷性及經(jīng)濟性，確定18 mm為噴嘴最優(yōu)長度。在噴嘴長度一定的情況下，模擬噴嘴直徑分別為18 mm，20 mm，24 mm而其他條件均相同的空心圓柱結(jié)構(gòu)噴嘴在276 MPa條件下內(nèi)、外部流場變化。

圖4　釋放管動壓力Fig.4　Dynamic pressure of the pipe

由圖3對3種不同直徑的壓力云圖分析可知，噴嘴直徑在18 mm和20 mm時，噴嘴處的噴射流暢及致裂長度相差不大；當噴嘴直徑達到24 mm時，噴嘴處的噴射流暢及致裂長度都明顯增大，且CO2氣體擴散范圍也比較大。總之，隨噴嘴直徑的增大，噴嘴口處的噴射流場也隨之變粗，且致裂段長度也增加。綜合分析得出，確定24 mm為釋放管噴嘴直徑的最優(yōu)尺寸。

圖3　不同直徑(a)空心圓柱結(jié)構(gòu)噴嘴動壓力Fig.3　Different diameter of a hollow cylindrical nozzle structure dynamic pressure nephogram

通過上述模擬分析，確定了釋放管噴嘴的最優(yōu)結(jié)構(gòu)和幾何參數(shù)，即最優(yōu)結(jié)構(gòu)(a)空心圓柱結(jié)構(gòu)，噴嘴最優(yōu)長度18 mm，最優(yōu)直徑24 mm。為進一步研究不同壓力對致裂效果的影響，分別模擬最優(yōu)噴嘴在126 MPa，190 MPa，236 MPa及276 MPa時不同壓力下內(nèi)外部流場變化，通過結(jié)果分析比較可知，隨壓力增大噴嘴處噴射段的長度和粗細程度都無明顯變化，但流場最大動壓力增大，因此確定最優(yōu)壓力為276 MPa。

在上述最優(yōu)模擬計算條件下，對最優(yōu)釋放管的流場進行模擬，其動壓力云圖如圖4所示。

從釋放管動壓力云圖得出：各釋放管由于釋放孔數(shù)目的不同，所產(chǎn)生的壓力等值線均不相同。由圖中可以看出單組釋放孔型釋放管產(chǎn)生的流場比較連續(xù)均勻，影響范圍較大，更有利于致裂煤層。而其他幾種釋放管由于釋放孔間產(chǎn)生流場的交互影響，流場分布不均，且影響范圍也較小。

4　CO2相變致裂增透技術(shù)現(xiàn)場應用

4.1　相變致裂增透技術(shù)方案的確定

相變致裂點數(shù)即在鉆孔內(nèi)布置1根或數(shù)根儲液管對煤層致裂，采用不同數(shù)量的液態(tài)CO2儲液管，由于CO2用量不同，煤層致裂效果各異。CO2致裂啟動后，將形成以鉆孔為中心的裂隙區(qū)，且裂隙區(qū)域的大小與液態(tài)CO2用量關(guān)系密切。為了更科學合理的指導液態(tài)CO2致裂技術(shù)，研究確定CO2用量與瓦斯抽放影響半徑的關(guān)系(圖5)，至關(guān)重要，通過對現(xiàn)場實測數(shù)據(jù)記錄分析，采用數(shù)值擬合的研究方法，初步獲得了致裂管數(shù)與鉆孔瓦斯抽放影響半徑的函數(shù)關(guān)系：

R=0.19N3-1.46N2+3.72N+25.3

式中：R為瓦斯抽放影響半徑，m；N為液態(tài)CO2儲液管數(shù)，個。

圖5　儲液管數(shù)量與瓦斯抽放影響半徑關(guān)系Fig.5　Liquid storage tube quantity relationship with gas drainage radius

由函數(shù)關(guān)系及圖5可知，鉆孔瓦斯抽放影響半徑與致裂管數(shù)呈三次函數(shù)關(guān)系，即隨著儲液管數(shù)量增多，瓦斯抽放孔影響半徑呈現(xiàn)回歸式拋物線趨勢增加，但增長率不大，同時考慮到施工程序及致裂成本，本次現(xiàn)場試驗儲液管數(shù)量定為1個，即單孔一點致裂。

4.2　試驗地點概括

試驗地點選在永錦能源云蓋山二礦23302工作面，該煤層位于山西組下部的二1煤層，煤層厚度0.2～9.7 m，平均厚度5.0 m。23302工作面地面標高+416 m，工作面開采上限+83 m，工作面開采下限+66.5 m，平均埋深342 m，煤層傾角8°～17°，23302工作面走向長度730 m，傾向長度165 m。煤層為Ⅲ類不易自燃煤層，煤塵爆炸指數(shù)17.99%，煤塵具有爆炸危險性，煤層節(jié)理比較發(fā)育，煤層頂?shù)装寰鶠樯百|(zhì)泥巖，透氣性較差。瓦斯含量19.35 m3/t，實測瓦斯壓力4.1 MPa。

4.3　相變致裂煤層透氣性分析

CO2相變致裂技術(shù)針對低透氣性高瓦斯煤層的增透瓦斯預抽效果明顯，在試驗區(qū)掘進巷180 m處沿巷道走向設(shè)置鉆孔，每組施工12個鉆孔，且分兩列布置，間距3 m，在鉆孔中間設(shè)置3個致裂孔，鉆孔布置見圖6。

圖6　鉆孔現(xiàn)場布置Fig.6　Drill site layout

CO2相變致裂后，煤層松動鉆孔內(nèi)裂隙擴展發(fā)育，煤層透氣性顯著增加；然而隨著煤體致裂后，要發(fā)生收縮變形，致裂后的裂隙自動閉合，煤層透氣性又逐漸降低，40 d后又趨于穩(wěn)定。為直觀形象觀測煤層透氣性變化，論文采用徑向流量法測量煤層透氣性系數(shù)，通過對比得知，趨于穩(wěn)定后透氣性系數(shù)提高19.83～24.45倍，增透效果顯著，見圖7。

圖7　煤層透氣性系數(shù)Fig.7　Permeability coefficient of coal seam

4.4　瓦斯抽采效果分析

圖8　致裂前后瓦斯抽采效果Fig.8　Effect of gas drainage before and after fracturing

本次現(xiàn)場試驗以CO2致裂后鉆孔的瓦斯平均濃度及抽采純量為評價指標，致裂前后其隨抽采時間的變化關(guān)系圖8所示。由圖8得出，致裂前10 d瓦斯平均濃度偏低，平均濃度為15.3%，在瓦斯抽采10 d以后，瓦斯平均濃度基本維持在10%以下，且瓦斯?jié)舛人p率也較快，抽采周期較短；而致裂后瓦斯平均濃度明顯提高，致裂后前30 d內(nèi)，瓦斯平均濃度在53.6%左右，在致裂50 d后，瓦斯平均濃度仍維持在45%以上，相對于致裂前提高了4倍左右，同時瓦斯衰減周期明顯增長，比致裂前至少提高2.5倍，且目前判斷，瓦斯抽采濃度依然沒出現(xiàn)衰減趨勢。由圖8得出，瓦斯抽采純量由致裂前平均1.8 m3/d提高到致裂后平均6 m3/d，提升幅度將近達到3.33倍，明顯提高了瓦斯抽采量。

綜上所述可知，低透氣性煤層采用CO2致裂技術(shù)，能夠明顯提高瓦斯的抽采濃度及抽采純量，并延長了衰減周期，取得了比較理想的增透促抽效果。

5　結(jié)論

1)CO2相變致裂技術(shù)利用高壓氣體沖擊破碎煤巖體，實現(xiàn)了煤巖體的裂隙發(fā)育，增大了煤巖層的透氣性，CO2致裂作用時，煤體損傷系數(shù)越大，有效應力強度因子也越大，增透致裂效果越顯著。

2)得到了1 kg液態(tài)CO2相變致裂爆炸當量，其爆炸的當量相當于398 gTNT炸藥。

3)確定了釋放管噴嘴的最優(yōu)結(jié)構(gòu)和幾何參數(shù)，即最優(yōu)結(jié)構(gòu)為(a)空心圓柱結(jié)構(gòu)，噴嘴最優(yōu)長度18 mm，最優(yōu)直徑24 mm，最優(yōu)壓力276 MPa；單組釋放孔型釋放管產(chǎn)生的流場比較連續(xù)均勻，影響范圍較大，更有利于致裂煤層。

4)得到了鉆孔瓦斯抽放影響半徑與致裂管數(shù)呈三次函數(shù)關(guān)系，提出了單孔一點致裂技術(shù)方案，現(xiàn)場應用取得了顯著的增透促抽效果。

[1]王兆豐，李豪君，陳喜恩.液態(tài)CO2相變致裂煤層增透技術(shù)布孔方式研究[J].中國安全生產(chǎn)科學技術(shù),2015,34(11):12-15.

WANG Zhaofeng, LI Haojun, CHEN Xi′en. Study on hole layout of liquid CO2phase-transforming fracture technology for permeability improvement of coal seam[J].Journal of Safety Science and Technology, 2015,34(11):12-15.

[2]魏國營，郭中海，謝倫榮，等.煤巷掘進水力掏槽防治煤與瓦斯突出技術(shù)[J].煤炭學報，2007，32(2)：172-176.

WEI Guoying，GUO Zhonghai，XIE Lunrong, et al. Hydraulic slotting technology to prevent coal and gas outburst during heading excavation [J].Journal Of China Coal Society. 2007，32(2)：172-176.

[3]周聲才，李棟，張鳳舞.煤層瓦斯抽采爆破卸壓的鉆孔布置優(yōu)化分析及應用[J].巖石力學與工程學報，2013，32(4)：807-808.

ZHOU Shengcai, LI Dong, ZHANG Fengwu. Optimization analysis of drilling layout based on blasting releasing pressure and its application [J]. Chinese Journal of Rock Mechanics and Engineering, 2013，32(4)：807-808.

[4]王兆豐，孫小明，陸庭侃.液態(tài)CO2相變致裂強化瓦斯預抽試驗研究[J].河南理工大學學報(自然科學版)，2015,34(1)：2-5.

WANG Zhaofeng，SUN Xiaoming，LU Tingkai. Experiment research on strengthening gas drainage effect with fracturing technique by liquid CO2phase transition[J].Journal of HeNan Polytechnic University(Natural Science), 2015,34(1)：2-5.

[5]龔敏，黃毅華，王德勝，等.松軟煤層深孔預裂爆破力學特性的數(shù)值分析[J].巖石力學與工程學報，2008，27(8)：1674-1681.

GONG Min, HUANG Yihua, WANG Desheng，et al. Numerical simulation on mechanical characteristics of deep-hole presplitting blasting in soft coal bed[J].Chinese Journal of Rock Mechanics and Engineering, 2008，27(8)：1674-1681.

[6]徐幼平，林柏泉，翟成，等.定向水力壓裂裂隙擴展動態(tài)特征分析及其應用[J].中國安全科學學報，2011，21(7):104-110.

XU Youping，LIN Baiquan，ZHAI Cheng，et al.Analysison dynamic characteristics of cracks extension in directional hydraulic fracturing and its application[J]. China Safety Science Journal，2011，21(7):104-110.

[7]胡國忠，許家林，王宏圖，等.低滲透煤與瓦斯的固-氣動態(tài)耦合模型及數(shù)值模擬[J].中國礦業(yè)大學學報，2011，40(1)：1-6.

HU Guozhong, XU Jialin, WANG Hongtu, et al. Research on a dynamically coupled deformation and gas flow model applied to low-permeability coal[J]. Journal of China University of Mining& Technology,2011，40(1)：1-6.

[8]唐書恒，馬彩霞，葉建平，等.注二氧化碳提高煤層甲烷采收率的實驗模擬[J].中國礦業(yè)大學學報，2006,35(5):606-611.

TANG Shuheng, MA Caixia, YE Jianping, et al.A modeling experiment of enhancing coalbed methane recovery by carbon dioxide injection[J]. Journal of China University of Mining& Technology, 2006,35(5):606-611.

[9]張春會，趙全勝，于永江.非均勻煤巖雙重介質(zhì)滲流-應力耦合模型[J].采礦與安全工程學報，2009,26(4)：481-485.

ZHANG Chunhui, ZHAO Quansheng, YU Yongsheng.Double media seepage-stress coupling model for heterogeneous coal and rock[J]. Journal of Mining敬& Safety Engineering, 2009,26(4)：481-485.

[10]趙陽升，楊棟，胡耀青，等.低滲透煤儲層煤層氣開采有效技術(shù)途徑的研究[J].煤炭學報，2001，26(5):455-458.

ZHAO Yangsheng, YANG Dong, HU Yaoqing, et al. Research on effective technical way of coal bed methane extraction in low permeability coal seam[J].Journal Of China Coal Society, 2001，26(5):455-458.

[11]董慶祥,王兆豐,韓亞北.液態(tài) CO2相變致裂的TNT當量研究[J].中國安全科學學報,2014,24(1):85-87.

DONG Qingxiang, WANG Zhangfeng, HAN Yabei. Research on TNT equivalent of liquid CO2phase-transition fracturing[J].China Safety Science Journal, 2014, 24(1):85-87.

[12]韓亞北.液態(tài)二氧化碳相變致裂增透機理研究[D].焦作: 河南理工大學，2014：7.

[13]趙萬春，艾池，李玉偉，等.基于損傷理論雙重介質(zhì)水力壓裂巖體劣化與孔滲特性變化理論研究[J].巖石力學與工程學報，2009(28)：3490-3496.

ZHAO Wanchun, AI Chi, LI Yuwei, et al. Study on rockmass deterioration and porosity and permeability in double-porosity medium under hydraulic fracture based on damage theory[J].Chinese Journal of Rock Mechanics and Engineering, 2009(28)：3490-3496.

[14]張建英，王兆豐，辛新平.高壓射流割縫卸壓增透技術(shù)在單一低滲透煤層瓦斯抽采中的應用[J].中州煤炭,2015(1):57-60.

ZHANG Jianying, WANG Zhaofeng, XIN Xinping. Application of high pressure water jet cutting technology in single low permeability coal seam gas extraction[J].Zhong Zhou Coal,2015(1):57-60.

[15]王海東.突出煤層掘進工作面CO2可控相變致裂防突技術(shù)[J].煤炭科學技術(shù)， 2016,44(3):70-73.

WANG Haidong.CO2controllable phase transition fracturing and outburst prevention technology of gateway driving face in outburst seam[J].Coal Science and Technology ,2016,44(3):70-73.

[16]張悅，張民波，朱天玲，等低透氣性煤層 CO2增透預裂技術(shù)應用[J].科技導報，2013，31( 23) : 36-39.

ZHANG Yue，ZHANG Minbo，ZHU Tianling，et al. Low permeability coal seam CO2anti-reflection pre-splitting technology application[J].Science & Technology Review，2013，31(23):36-39.

[17]范迎春.二氧化碳爆破增透技術(shù)在低透氣性煤層中的應用研究[J].中州煤炭,2015(5):1-3.

FAN Yingchun .Application research on low permeability seam of carbon dioxide blasting fracturing technology[J].Zhongzhou Coal, 2015(5):1-3.

[18]霍中剛.二氧化碳致裂器深孔預裂爆破煤層增透新技術(shù)[J].煤炭科學技術(shù)，2015，43( 2) : 80-83.

HUO Zhonggang.New technology of carbon dioxide fracture applied to deep borehole pre-cracking blasting for seam permeability improvement[J].Coal Science and Technology，2015，43(2):80-83.

[19]楊天鴻，唐春安，徐濤.巖石破裂過程滲流特性-理論模型與應用[M].北京: 科學出版社，2004：69-70.

YANG Tianhong, TANG Chunan, XU Tao. Seepage characteristics of rock failure process: a theoretical model and its application[M]；Beijing: Science Press, 2004：69-70.

[20]李文煒，狄剛，王瑞欣.船運液態(tài)CO2儲罐爆炸事故的原因分析[J].安全與環(huán)境工程，2010，17(1):95-98.

LI Wenwei, DI Gang, WANG Runxin. Analysis of a liquid CO2tank explosion on a ship[J]. Safety and Environmental Engineering,2010,17(1):95-98.