CO2相變爆破致裂增透裂紋擴(kuò)展數(shù)值模擬研究

王曉東

(華陽新材料科技集團(tuán)有限公司技術(shù)中心，山西 陽泉 045000)

1 前言

煤層氣作為近年來不斷崛起的優(yōu)質(zhì)能源逐漸走進(jìn)人們的視野。煤層氣又稱“瓦斯”，是一種儲存在煤層內(nèi)以甲烷為主要?dú)怏w、部分游離于煤孔隙中或溶解于煤層水中的烴類氣體。我國煤層氣資源儲量極其豐富，其總儲量可排世界前三位。煤層氣作為一種“綠色能源”“可持續(xù)發(fā)展型能源”，既是煤礦井下災(zāi)害事故的致因因素之一，又是造成溫室效應(yīng)的一種有害氣體源[1]。因此，對于煤層氣資源的開發(fā)利用可有效增加新能源在能源結(jié)構(gòu)中的占比、改善保障煤礦生產(chǎn)安全、減少環(huán)境污染。受地質(zhì)條件以及煤炭沉降條件的影響，我國大部分煤層的透氣性較低，實施瓦斯抽采相對困難。對于單一煤層而言，現(xiàn)有的瓦斯抽采解決辦法有設(shè)置密集鉆孔或者對單一煤層進(jìn)行長時間抽采，目前所采取的這些方法不僅提高了抽采成本，而且降低了抽采效率。因此，增加煤層的透氣性，提高煤層瓦斯抽采效率以及降低瓦斯抽采的成本是目前急需解決的問題。

液態(tài)CO2相變爆破煤層致裂增透技術(shù)作為一種新型提高瓦斯抽采效率的技術(shù)，實施方便簡單，成本低廉，整體安全性高，受到了眾多學(xué)者以及煤炭行業(yè)從業(yè)者的青睞。張健[2]利用二氧化碳相變爆破技術(shù)在煤礦開展現(xiàn)場爆破試驗，極大地提高了煤層的瓦斯抽采效率，并且也有效地解決壓力集中導(dǎo)致的煤和瓦斯突出和掘進(jìn)工作面應(yīng)力集中的問題。洪林[3]利用FLAC3D與Comsol模擬軟件模擬了二氧化碳相變爆破，研究認(rèn)為液態(tài)二氧化碳多孔連續(xù)爆破能夠達(dá)到最理想的爆破增透效果，且控制孔的半徑越大，增透效果越明顯，瓦斯抽采效率越高。黃曉實[4]將二氧化碳爆破致裂器應(yīng)用在巖石深孔預(yù)裂爆破中，分析了致裂器噴嘴形狀對于爆破致裂效果的影響。

前人對于二氧化碳相變爆破煤層致裂的研究多側(cè)重于現(xiàn)場試驗中[5-10]，較少利用實驗室數(shù)值模擬研究其致裂的效果。因此，將二氧化碳相變爆破煤層致裂現(xiàn)場試驗與LS-DYNA數(shù)值模擬中煤層爆破裂紋擴(kuò)展表征致裂效果相結(jié)合，對于提高二氧化碳相變爆破致裂效果，提高爆破致裂效率有著良好的指導(dǎo)意義。

2 CO2相變爆破致裂煤層增透技術(shù)及原理分析

CO2相變爆破致裂煤層增透技術(shù)主要是通過液態(tài)CO2瞬間受熱由液態(tài)轉(zhuǎn)化為氣態(tài)，體積迅速增大膨脹產(chǎn)生的巨大壓力，壓力產(chǎn)生的應(yīng)力波和高壓氣體共同作用于周圍煤體域，促使原煤層裂隙擴(kuò)展的同時產(chǎn)生新的裂隙 ，增大煤層透氣性的同時提高煤層內(nèi)瓦斯抽采效率[11]。

現(xiàn)場試驗中所用CO2致裂器如圖1所示。其操作過程分為4個步驟：CO2致裂器的拼裝、向致裂器內(nèi)沖入適當(dāng)?shù)囊簯B(tài)CO2、將制備好的致裂器安裝在井下和致裂起爆。液態(tài)CO2充裝設(shè)備由液態(tài)CO2鋼瓶、壓力泵、空氣壓縮機(jī)等設(shè)備組成，利用壓力泵將液態(tài)CO2充入特制的儲液管內(nèi)，并將壓力保持在固定范圍。組裝其他部件，將制備好的致裂器安裝在井下準(zhǔn)備起爆，將封孔器用封孔高壓管路與致裂器連接，并將發(fā)爆母線固定于頂桿之上。將頂桿推進(jìn)到符合設(shè)計要求的深度。檢查孔內(nèi)設(shè)備的連接情況，確認(rèn)可靠后，連接壓力泵及測壓表。用壓力泵向封孔器內(nèi)注水至一定壓力后，調(diào)節(jié)頂桿，固定致裂器。安裝工藝如圖2所示。

1—主管；2—充氣閥；3—泄能閥；4—發(fā)熱裝置；5—電極引出體

1—CO2致裂器；2—連接件；3—封孔器；4—封孔引出桿；5—固定器

3 CO2相變爆破致裂增透裂紋數(shù)值模擬及結(jié)果分析

3.1 有限元模型的建立

有限元模型的建立以現(xiàn)場試驗相關(guān)數(shù)據(jù)為基礎(chǔ)。建立40 m×20 m的二維模型，鉆孔直徑0.094 m。根據(jù)文獻(xiàn)[12]中的致裂半徑公式，代入表1中的具體參數(shù)進(jìn)行計算可得致裂半徑約為8.4 m。因現(xiàn)場試驗鉆孔布置數(shù)量為16個，為簡化模型及減少計算量，有限元模型中爆破鉆孔數(shù)量布置為3個，通過LS- PrePost軟件設(shè)置點火時間一致控制三個鉆孔同時起裂，分別建立鉆孔相距10 m、5 m、3 m三個數(shù)值模擬計算模型。在所建立的有限元模型中，炸藥采用*MAT_HIGH_EXPLOSIVE_BURN。空氣模型采用*MAT_NULL，煤體采用 *MAT_PLASTIC_KINEMATIC[13-16]。以鉆孔間距10 m的數(shù)值模型為例，建立的有限元模型如圖3所示，模型中主要參數(shù)具體數(shù)值見表1。

圖3 有限元模型(鉆孔間距10 m)

表1 有限元模型參數(shù)值設(shè)置

根據(jù)上表建立有限元模型，由于該模型模擬的是煤體爆炸區(qū)域的一部分，為保證模擬的準(zhǔn)確性及提高計算效率，在煤體域邊界添加無反射邊界條件，計算時長設(shè)置為0.002 s。

3.2 數(shù)值模擬結(jié)果分析

本次數(shù)值模擬設(shè)置了爆破鉆孔相距10 m、5 m、3 m三種不同條件下二氧化碳相變爆破致裂模型，通過數(shù)值模擬結(jié)果裂紋發(fā)展對比分析致裂效果。

1)鉆孔相距10 m煤體爆破裂紋效果

圖4所示為鉆孔間距為10 m二氧化碳相變爆破致裂裂紋擴(kuò)展隨時間變化的壓力云圖。

圖4 鉆孔間距10 m爆破致裂裂紋擴(kuò)展及壓力云圖

根據(jù)圖4可以看出，在t=0.000 2 s時刻，二氧化碳相變爆破沖擊波傳播至煤體域并逐漸在煤體域內(nèi)形成爆破致裂空腔區(qū)。在t=0.000 4 s時，沖擊波繼續(xù)向四周傳播，逐漸形成煤體域的壓碎區(qū)。t=0.001 6 s時，爆破致裂裂紋擴(kuò)展發(fā)育情況逐漸穩(wěn)定，煤層二氧化碳相變爆破致裂區(qū)域明顯，通過裂紋拓展區(qū)域可判斷煤層致裂范圍。

2)鉆孔相距5 m煤體爆破裂紋效果

圖5所示為鉆孔間距為5 m二氧化碳相變爆破致裂裂紋擴(kuò)展隨時間變化的壓力云圖。

圖5 鉆孔間距5 m爆破致裂裂紋擴(kuò)展及壓力云圖

根據(jù)圖5爆破鉆孔間距5 m條件下不同時刻爆破致裂煤體域壓力云圖變化，可以看出在t=0.000 2 s時刻，二氧化碳相變爆破沖擊波傳播至煤體域并逐漸在煤體域內(nèi)形成爆破致裂空腔區(qū)，該階段致裂耗時與鉆孔間距10 m模擬數(shù)值模擬結(jié)果一致。根據(jù)上圖可知，未到t=0.000 4 s時刻，即已形成爆破致裂壓碎區(qū)，該階段耗時較短。t=0.001 3 s時，煤體域裂紋區(qū)發(fā)育擴(kuò)展趨于穩(wěn)定，通過裂紋可判斷爆破致裂半徑范圍。

3)鉆孔相距3 m煤體爆破裂紋效果

圖6所示為鉆孔間距為3 m二氧化碳相變爆破致裂裂紋擴(kuò)展隨時間變化的壓力云圖。

圖6 鉆孔間距3 m爆破致裂裂紋擴(kuò)展及壓力云圖

根據(jù)圖6可知，在鉆孔間距為3 m條件下t=0.000 2 s時爆破沖擊波傳播至煤體域，沖擊波逐漸向煤體四周傳播并逐漸形成裂紋。t=0.001 8 s時裂隙區(qū)發(fā)育完成。

由圖4、5、6爆破裂紋裂隙拓展及壓力云圖變化情況可知，鉆孔間距設(shè)置對于同等裝藥量條件下爆破致裂效果有著顯著影響，隨著鉆孔間距設(shè)置的增大，二氧化碳爆破致裂效果更佳。結(jié)合二氧化碳相變爆破致裂現(xiàn)場試驗結(jié)果與數(shù)值模擬計算結(jié)果可知，在二氧化碳相變爆破實際煤體增透過程中，爆破鉆孔間距設(shè)置為10 m情況下可獲得理想的最佳爆破致裂效果。

4 CO2相變爆破致裂增透現(xiàn)場試驗分析

4.1 試驗設(shè)計

根據(jù)以往抽采情況，在貴州省安順市宏發(fā)煤礦1903工作面回風(fēng)巷進(jìn)行普通抽采鉆孔抽采瓦斯效果不太理想，而CO2致裂是提高瓦斯抽采效率的有效手段，且現(xiàn)場試驗中的煤體的煤層堅固性系數(shù)f值為0.54(煤層堅固性系數(shù)大于突出臨界值 0.5，表明煤層煤質(zhì)較硬)，符合CO2致裂的應(yīng)用條件。在貴州省安順市宏發(fā)煤礦1903工作面回風(fēng)巷開展CO2相變爆破致裂增透現(xiàn)場試驗，根據(jù)1903工作面回風(fēng)巷具體施工特點，實際在1903工作面回風(fēng)巷下幫布置14個致裂孔，上幫布置2個致裂孔，致裂孔布置示意圖如圖7所示。

圖7 鉆孔布置示意圖

根據(jù)鉆孔的具體施工深度確定致裂器的數(shù)量，致裂前檢查測量作業(yè)地點附近的瓦斯、CO2濃度，待瓦斯?jié)舛鹊陀?.8%、CO2濃度低于0.5%時進(jìn)行致裂工作。

4.2 致裂試驗數(shù)據(jù)分析

透氣性系數(shù)對比考察采用劉明舉、何學(xué)秋提出的煤層透氣性系數(shù)方法測定[17]。向煤層打鉆孔，封孔的同時安置瓦斯壓力表，待表示數(shù)穩(wěn)定后將其拆除，并安裝煤氣表采用日測1次的方法測定鉆孔瓦斯自然流量，直至涌出量穩(wěn)定。根據(jù)煤層徑向不穩(wěn)定流動理論結(jié)合測得的煤層瓦斯壓力、瓦斯含量及鉆孔基本參數(shù)綜合計算，最終確定煤層的透氣性系數(shù)。所涉及到的計算公式為

(1)

(2)

式中：A——參數(shù)；

B——參數(shù)；

p——煤層瓦斯壓力，MPa；

p0——絕對瓦斯壓力，通常取0.1 MPa；

r1——鉆孔半徑，m；

λ——煤層透氣性系數(shù)；

由式(1)～(2)可知A、B的數(shù)值，然后任選F值根據(jù)式(3)～(8)反推λ，再根據(jù)式(9)求出F值，判斷是否符合范圍。以此方法反復(fù)推演最終計算得出煤層的透氣性系數(shù)，具體范圍下的選用公式如式(4)～(9)所示。

λ=A1.61·B1/1.64(F=10-2～1)

(3)

λ=A1.39·B1/2.56(F=1～10)

(4)

λ=1.10·A1.25·B1/4(F=10～102)

(5)

λ=1.83·A1.14·B1/7.3(F=102～103)

(6)

λ=2.1·A1.11·B1/9(F=103～105)

(7)

λ=3.14·A1.07·B1/14.4(F=105～107)

(8)

F=B·λ

(9)

式中：T——從排放瓦斯開始到測到瓦斯流量q時的時間間隔，d；

q——煤層單位暴露面積上的瓦斯流量，m3/m2·d，q=Q/(2πr1L)；

Q——T時刻的瓦斯流量，m3/d；

L——鉆孔見煤長度，m，可取煤層厚度；

λ——煤層的透氣性系數(shù)，m2/MPa2·d。

未1、未3孔為未致裂增透區(qū)域的煤層透氣性系數(shù)測定孔，作為原始對照組；6#、7#、8#致裂孔為致裂增透后煤層透氣性系數(shù)的測定孔，作為實驗對照組。對各對照組的瓦斯自然流量進(jìn)行測量，結(jié)果見表2，鉆孔瓦斯自然流量測試變化曲線如圖8所示。

表2 煤層鉆孔瓦斯自然流量測試結(jié)果

結(jié)合測量數(shù)據(jù)及計算公式計算可得各鉆孔的煤層透氣性系數(shù)數(shù)值見表3。

表3 煤體致裂后煤層的透氣性系數(shù)

液態(tài)CO2相變爆破致裂增透前該煤層的透氣性系數(shù)λ為0.126～0.135 m2/MPa2·d。進(jìn)行爆破致裂增透后，煤層的透氣性系數(shù)λ上升為0.363～0.441 m2/MPa2·d。表明液態(tài)CO2相變爆破煤層致裂后，煤層的透氣性顯著提高其數(shù)值增大為原始煤層的3倍，說明液態(tài)CO2相變爆破致裂增透技術(shù)對于煤層的增透和提高低透氣性煤層瓦斯抽采效率均有較好的意義及效果。

5 結(jié)論

(1)利用LS-DYNA數(shù)值模擬二氧化碳相變爆破鉆孔間距設(shè)置為10 m、5 m、3 m三個數(shù)值模擬計算模型，分析其裂紋擴(kuò)展致裂效果及壓力云圖變化情況可知，在爆破鉆孔間距設(shè)置為10 m情況下，爆破致裂裂紋擴(kuò)展發(fā)育情況最佳，可獲得良好的破碎致裂效果。

(2)開展現(xiàn)場爆破致裂增透裂紋擴(kuò)展試驗后，煤層的透氣性系數(shù)λ上升為0.363～0.441 m2/MPa2·d，增大為原始煤層的3倍左右，煤層的透氣性顯著提高。表明該項技術(shù)能有效的提高煤層氣開采的效率對于開采災(zāi)害的降低以及煤礦井下安全事故的預(yù)防均有一定意義。