張 飛，尚瑋煒，羅華貴，王 寧，賈進章，李 斌

(1.山西晉煤集團技術研究院有限責任公司 通風安全技術服務分公司，山西 晉城 048006；2.遼寧工程技術大學 安全科學與工程學院，遼寧 阜新 123000；3.礦山熱動力災害與防治教育部重點實驗室，遼寧 葫蘆島 125105)

目前，我國有3/5以上的煤礦都是高瓦斯礦井且具有突出危險性，這些礦井不僅威脅井下工人的個人安全，還嚴重地影響到我國煤炭資源的開采效率[1-2]。針對煤層瓦斯治理問題，當采用通風的方式無法有效降低瓦斯在巷道內(nèi)的濃度時，應進行煤層瓦斯抽采治理，通過煤層鉆孔抽采，快速地降低煤層中賦存的瓦斯含量[3-4]。煤層鉆孔瓦斯抽采有效半徑對瓦斯抽采效果的影響較大，準確預測煤層鉆孔瓦斯抽采有效半徑是提高瓦斯抽采率和減少瓦斯災害的重要技術手段[5]。

針對低滲透高瓦斯煤層抽采有效半徑測定的問題，科研人員通過理論分析、數(shù)值模擬和現(xiàn)場試驗開展了大量的研究工作，也取得了許多重要成果。鄒士超等[6]理論分析了徑向流場中瓦斯壓力的分布規(guī)律，并以0.5 MPa瓦斯壓力為臨界點，利用模擬方法得到了四季春煤礦6號煤層瓦斯抽采有效半徑范圍，現(xiàn)場驗證誤差小于10%；張明杰等[7]在瓦斯流動理論基礎上，研究了鉆孔內(nèi)瓦斯含量、涌出量和殘存量三者之間的關系，提出了基于鉆孔瓦斯自然涌出規(guī)律來測定抽采有效半徑的方法，并進行現(xiàn)場應用，測試結(jié)果誤差在3%左右；王冬[8]根據(jù)井下現(xiàn)場實際情況，利用模擬軟件分析了瓦斯抽采有效半徑與抽采周期之間的關系，為青崗坪煤礦鉆孔布置提供了指導；李潤芝等[9]通過測定鉆孔周圍煤體特定時刻的瓦斯含量來測算煤層鉆孔瓦斯抽采有效半徑，經(jīng)現(xiàn)場驗證，準確率較高；王兆豐等[10]提出了利用非線性滲流方程來測定瓦斯抽采半徑，工程試驗表明此方法能夠優(yōu)化布孔間距，提高抽采效率，節(jié)約工程成本；舒才等[11]建立了利用瓦斯抽采量確定抽采半徑的數(shù)學模型，并將該模型應用在興隆煤礦進行驗證，現(xiàn)場實測結(jié)果與模型預測結(jié)果基本一致。

上述研究成果雖然推動了瓦斯抽采技術的應用和發(fā)展，但由于影響抽采半徑的因素較多，各因素之間關系又模糊不清，現(xiàn)有計算瓦斯抽采半徑的方法各異，導致相關研究缺乏統(tǒng)一的評價抽采半徑標準。鑒于此，筆者擬研究瓦斯抽采有效半徑影響因素的主次順序，厘清各影響因素對瓦斯抽采有效半徑的影響程度。借助COMSOL軟件，模擬瓦斯壓力、煤層滲透率、抽采負壓、孔徑及抽采周期對瓦斯抽采有效半徑的影響；應用灰色關聯(lián)分析法對這5個因素進行分析，確定多種影響因素中的主控因素，以期為煤層鉆孔瓦斯抽采技術的應用研究提供理論參考與借鑒。

1 煤層瓦斯抽采數(shù)學模型的構建

煤層預抽是降低瓦斯?jié)舛鹊谋匾胧?，然而進行瓦斯抽采數(shù)值模擬時，往往受到多方面的影響，為減少這些因素的影響，給出如下基本假設：假設煤層頂板及底板密封不透氣；煤層透氣性系數(shù)不受瓦斯壓力變化影響；瓦斯為理想氣體且運移規(guī)律符合達西定律；瓦斯吸附規(guī)律適用朗格繆爾方程[11-12]?；谏鲜鰲l件，構建瓦斯抽采數(shù)學模型如下：

1)煤層孔隙率的動態(tài)方程

瓦斯抽采過程中鉆孔周圍的煤體孔隙率處于動態(tài)變化的過程中，假設在煤層中僅僅含有瓦斯，則根據(jù)定義，孔隙率φ可表達為[13]：

(1)

式中：Vp為煤的孔隙體積，m3；Vt為煤的總體積，m3；Vp0為初始孔隙體積，m3；Vt0為初始總體積，m3；ΔVp為孔隙體積變化量，m3；ΔVt為總體積變化量，m3；φ0為初始孔隙率；εV為體積應變；ΔVs為骨架體積變化量，m3；Vs0為初始骨架體積，m3；Kγ為壓縮的體積系數(shù)；Δp為瓦斯壓力的變化量，MPa。

2)煤層瓦斯?jié)B透率的動態(tài)方程

煤體的瓦斯壓力變化能夠引起孔隙率φ和滲透率k發(fā)生動態(tài)變化，二者之間存在的關系可以通過Kozeny-Carman方程來表示[14]：

(2)

式中：Kz為無量綱常數(shù)；As為孔隙表面積，m2。

3)瓦斯流動控制方程

根據(jù)質(zhì)量守恒定律，瓦斯總的質(zhì)量變化量等于單位時間內(nèi)瓦斯流入、流出的質(zhì)量差[15]，可以表示為：

(3)

式中：Q為煤體內(nèi)瓦斯含量，m3/t；ρ為瓦斯密度，kg/m3；v為滲流速度，m/s。

瓦斯含量在單位體積煤中可表達為[15]：

(4)

式中：β為瓦斯壓縮系數(shù)，kg/(m3·Pa)；p為瓦斯壓力，MPa；a為質(zhì)量極限吸附量，m3/t；b為吸附常數(shù)；c為校正系數(shù)；ρn為標準狀態(tài)下瓦斯密度，kg/m3。

瓦斯在煤體中的流動規(guī)律適用達西定律，同時考慮到Klinkenberg效應對瓦斯流動也有一定程度的影響，故瓦斯的滲流速度可以表示為[16]：

(5)

式中：μ為動力黏度系數(shù)，Pa/s；m為Klinkenberg系數(shù)。

由公式(3)、(4)、(5)，可得瓦斯流動控制方程：

(6)

2 煤層抽采數(shù)值模型構建及邊界條件

使用COMSOL數(shù)值模擬軟件進行模擬，以寺河煤礦的6301工作面為基本參數(shù)選取模型參數(shù)。通過井下實測，具體參數(shù)見表1。

表1 煤層基礎參數(shù)

模型尺寸長為40 m，高為 25 m，其中煤層的厚度為5 m，如圖1所示。頂?shù)装鍨椴煌笟鈳r石，模型底部固定，上覆巖層壓力為 13 MPa，鉆孔位于矩形中間，深度為 60 m，鉆孔邊界壓力為抽采負壓。為使模擬結(jié)果更為準確，模型采用自定義劃分網(wǎng)格，網(wǎng)格的最大單元選擇0.2 m，最小單元為0.1 m，網(wǎng)格的曲率因子取0.7。

圖1 單孔抽采鉆孔數(shù)值模型

3 抽采有效半徑影響因素模擬結(jié)果分析

大量研究結(jié)果表明，煤層瓦斯抽采有效半徑的主要影響因素有5個方面：孔徑、抽采周期、初始瓦斯壓力、滲透率、抽采負壓強度[17-18]。利用COMSOL 模擬軟件，分別模擬以上5個因素對煤層瓦斯抽采有效半徑的影響情況。依據(jù)煤礦相關標準[19]，以殘存瓦斯含量Qc=8 m3/t作為有效半徑確定的指標，通過分析以上5個因素對抽采有效半徑的影響情況，得出影響煤層瓦斯抽采鉆孔有效半徑的主控因素。

3.1 初始瓦斯壓力對抽采有效半徑的影響分析

在數(shù)值模擬的過程中，分別取煤層的初始瓦斯壓力p0為1.0、1.2、1.4、1.6、1.8 MPa，模擬結(jié)果如圖2所示，在150 d抽采周期內(nèi)，上述5個不同初始瓦斯壓力值所對應的初始瓦斯含量Q0分別約為8.92、10.36、12.51、13.59、14.94 m3/t。初始瓦斯壓力為 1.0 MPa時，抽采有效半徑約為1.79 m，其余4種情況下，鉆孔的抽采有效半徑在0.30～0.75 m，抽采效果相對較差。

圖2 不同初始瓦斯壓力對瓦斯含量影響曲線

從圖2可以看出，煤層瓦斯抽采有效半徑隨著煤層的初始瓦斯壓力的增大而減小(瓦斯含量變化曲線與Qc交點所對應的橫坐標值為有效半徑，下同)。因此，在井下實際施工過程中，煤層瓦斯抽采鉆孔間距和位置要根據(jù)初始瓦斯壓力進行合理布置，以達到最佳的抽采效果。

抽采鉆孔在不同初始瓦斯壓力條件下，抽采30 d和150 d瓦斯含量變化率曲線如圖3所示。模擬結(jié)果表明，煤層的瓦斯含量變化率隨著煤層初始瓦斯壓力的增大而減小，這說明煤層的初始瓦斯壓力過大時，需要投入更大的瓦斯治理成本，才能將瓦斯含量降至安全范圍之內(nèi)，以避免災害的發(fā)生。

(a)抽采30 d

3.2 滲透率對抽采有效半徑的影響分析

滲透率分別取1.5×10-18、3.5×10-18、6.5×10-18、9.5×10-18m2，煤層瓦斯抽采有效半徑模擬結(jié)果見圖4,測得抽采有效半徑分別約為1.7、3.0、4.5、5.8 m。

圖4 不同滲透率對瓦斯含量影響曲線

由圖4可以看出，煤層的滲透率越大，煤層的抽采有效半徑越大。原因在于煤層滲透率大，有助于游離態(tài)瓦斯在抽采負壓狀態(tài)下的運移，有利于瓦斯抽采。對于滲透率較高的煤層在實施瓦斯抽采時，可適當增加布孔間距來節(jié)約成本。

3.3 抽采負壓對抽采有效半徑的影響分析

鉆孔抽采負壓分別取10、15、20、25 kPa，煤層抽采有效半徑測試結(jié)果如圖5所示，測得抽采有效半徑約為1.66、1.72、1.80、1.89 m。

圖5 不同抽采負壓對瓦斯含量影響曲線

由圖5可以看出，抽采負壓越高，抽采有效半徑越大。原因在于抽采負壓所產(chǎn)生的壓力差是煤體內(nèi)瓦斯運移的主要動力，壓力梯度越大，孔隙中瓦斯流動的速度越快，從而提高瓦斯抽采效率。但抽采負壓不宜過高，否則會加劇鉆孔密封段的漏風。因此，在井下現(xiàn)場瓦斯抽采時，可以適當增加負壓值來提高抽采有效半徑的范圍。

3.4 鉆孔直徑對抽采有效半徑的影響分析

抽采鉆孔直徑分別取94、113、153 mm，煤層抽采有效半徑測試結(jié)果分別約為1.81、1.84、1.95 m，如圖6所示。

圖6 不同鉆孔直徑對瓦斯含量影響曲線

從圖6可以看出，抽采鉆孔的孔徑越大，抽采有效半徑越大。原因在于抽采鉆孔孔徑越大，原煤暴露出的部分越多，施工時，煤體受到的擾動就越強，從而使煤層瓦斯抽采有效半徑增大。因此，適當增大孔徑可以提高煤層的抽采效果，在實際施工過程中，建議抽采孔孔徑設置為94 mm或113 mm，因為孔徑過小會影響瓦斯抽采效果，孔徑過大，施工過程中會出現(xiàn)塌孔現(xiàn)象并且容易誘發(fā)瓦斯突出事故。

3.5 抽采周期對抽采有效半徑的影響分析

抽采周期分別選取30、90、150、210、270、330、360 d，瓦斯抽采有效半徑測試結(jié)果分別約為1.1、1.8、2.4、2.9、3.5、4.0、4.3 m，如圖7所示。

圖7 不同抽采周期下瓦斯含量變化曲線

由圖7可以看出，隨著抽采周期的增加，煤層瓦斯抽采有效半徑越來越大，但相對的變化率卻越來越小。對于低滲透高瓦斯煤層，適當增加抽采周期可以降低井下瓦斯事故的威脅。

4 鉆孔瓦斯抽采有效半徑主控因素確定

根據(jù)上述5個瓦斯抽采參數(shù)對抽采有效半徑影響情況，引入灰色關聯(lián)理論[20]，對這5個抽采參數(shù)與抽采有效半徑之間的關聯(lián)度做進一步主次順序分析。

選取抽采有效半徑作為母因素，子因素分別為初始瓦斯壓力、滲透率、抽采負壓、孔徑和抽采周期。各影響因素的取值如表2所示。

表2 灰色關聯(lián)分析因素及其取值

灰色關聯(lián)度理論數(shù)學模型建立如下[20]：

1)絕對關聯(lián)度ε0i計算公式如下：

(7)

(8)

(9)

(10)

式中：x0(k)為參照數(shù)列;xi(k)為對比數(shù)列。

采用MATLAB軟件計算灰色絕對關聯(lián)度，計算結(jié)果見表3。

表3 各因素與抽采有效半徑灰色絕對關聯(lián)度數(shù)據(jù)

2)在計算灰色相對關聯(lián)度γ0i時，需要進行無量綱化處理，以確保數(shù)據(jù)的統(tǒng)一性。采用均值化方法對抽采有效半徑數(shù)據(jù)(參考數(shù)列)和5個影響因素數(shù)據(jù)(對比數(shù)列)進行無量綱化處理，計算公式如下[20]：

(11)

經(jīng)處理后，將數(shù)據(jù)代入公式(7)～(10)，便可以得到灰色相對關聯(lián)度γ0i。

灰色相對關聯(lián)度計算結(jié)果見表4。

表4 各因素與抽采有效半徑灰色相對關聯(lián)度數(shù)據(jù)

3)灰色綜合關聯(lián)度(δ)涵蓋灰色絕對和相對關聯(lián)度共同的優(yōu)勢，既可以分析數(shù)據(jù)之間的相關性，又可以比較相對于數(shù)據(jù)始點變化趨勢，其計算公式如下[20]：

δ=λε0i+(1-λ)γ0i

(12)

式中λ為分配系數(shù)，取0.5。

灰色綜合關聯(lián)度計算結(jié)果見表5。

表5 各因素與抽采有效半徑灰色綜合關聯(lián)度數(shù)據(jù)

從表5可以看出，初始瓦斯壓力的綜合關聯(lián)度最大，其次是滲透率、抽采周期、孔徑，抽采負壓的綜合關聯(lián)度最小，即這5個因素對煤層瓦斯抽采有效半徑的影響程度排序為：初始瓦斯壓力>滲透率>抽采周期>孔徑>抽采負壓。因此，在工程實踐中進行瓦斯抽采時，應該首先勘測本煤層初始瓦斯壓力和滲透率，當煤層初始瓦斯壓力較高時，應適當增大抽采鉆孔的孔徑、抽采負壓和抽采周期，這樣才能達到更好的煤層瓦斯抽采效果。

5 結(jié)論

1)通過構建瓦斯抽采數(shù)學模型，借助模擬軟件分析得到煤層瓦斯抽采有效半徑與初始瓦斯壓力呈遞減關系，與滲透率、抽采負壓、孔徑和抽采周期呈遞增關系。

2)通過構建灰色關聯(lián)分析數(shù)學模型，得到各因素對煤層瓦斯抽采有效半徑的影響程度綜合排序為：初始瓦斯壓力>滲透率>抽采周期>抽采鉆孔直徑>抽采負壓。

3)在礦井的同一工作面中，煤層的初始瓦斯壓力和滲透率幾乎保持一致，因此，為了提高煤層瓦斯抽采有效半徑范圍，應該適當?shù)卦黾映椴芍芷?、孔徑和抽采負壓?/p>