宋浩然,林柏泉,趙 洋,孔 佳,查 偉
(1.中國(guó)礦業(yè)大學(xué) 煤礦瓦斯與火災(zāi)防治教育部重點(diǎn)實(shí)驗(yàn)室,江蘇 徐州 221116;2.中國(guó)礦業(yè)大學(xué) 安全工程學(xué)院,江蘇 徐州 221116)
我國(guó)深部煤層瓦斯普遍具有高地應(yīng)力、高瓦斯壓力、高瓦斯含量以及低滲透性特征[1-2],各種復(fù)雜因素使得當(dāng)前煤礦開采和瓦斯抽采難度較大,礦山安全生產(chǎn)受到嚴(yán)重威脅[3]。瓦斯是引起礦井動(dòng)力災(zāi)害的主要因素之一,容易導(dǎo)致煤礦發(fā)生爆炸和煤與瓦斯突出危險(xiǎn)[4]。因此,為了提高瓦斯的資源化利用并減少礦井災(zāi)害發(fā)生,提高鉆孔瓦斯抽采效率十分必要[5]。目前,水力割縫卸壓瓦斯抽采技術(shù)是實(shí)現(xiàn)煤礦井下瓦斯資源化的主要手段之一[6],也是防治瓦斯災(zāi)害的重要手段[7]。對(duì)此,相關(guān)學(xué)者已經(jīng)開展了大量試驗(yàn)研究。Candela等進(jìn)行了一系列關(guān)于滲透率增強(qiáng)機(jī)制的試驗(yàn)研究,從微觀角度探索了煤層的卸壓增透機(jī)理[8-9]。尹光志等利用三軸滲透儀對(duì)突出煤進(jìn)行了試驗(yàn)研究,分析了瓦斯壓力對(duì)滲流的影響作用機(jī)制[10]。魯義等開展了單鉆孔和多鉆孔協(xié)同瓦斯抽采試驗(yàn),驗(yàn)證了抽采半徑與布孔間距關(guān)系的正確性[11]。數(shù)值模擬是研究低滲透煤層瓦斯抽采的有效方法,許多學(xué)者考慮過不同物理場(chǎng)的耦合效應(yīng)。Lu等建立了鉆孔水力壓裂后煤層瓦斯流動(dòng)的多物理場(chǎng)耦合模型,分析了煤層瓦斯壓力的變化規(guī)律和滲透率的演化規(guī)律[12]。Zhou 等建立并模擬了深部煤層的氣-固耦合模型,探索了高瓦斯煤礦防止瓦斯突出的有效方法[13]。Wei 等分析了瓦斯抽采過程中有效應(yīng)力的變化以及煤基質(zhì)解吸對(duì)煤體孔隙度影響的作用機(jī)理,建立了滲透率動(dòng)態(tài)演化模型[14]。林柏泉等建立了應(yīng)力場(chǎng)、滲流場(chǎng)和擴(kuò)散場(chǎng)多場(chǎng)耦合模型,研究了瓦斯抽采過程中煤層瓦斯流場(chǎng)的演化規(guī)律,并結(jié)合楊柳煤礦的現(xiàn)場(chǎng)實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證了模型的準(zhǔn)確性。煤層滲透率具有各向異性,即煤層平行層理滲透率與垂直層理滲透率存在差異[15]。Anggara等通過實(shí)驗(yàn)室測(cè)量,證實(shí)了煤層平行層理2個(gè)方向的滲透率基本相同,而各向異性的主要差異在平行層理和垂直層理2個(gè)方向[16]。Tan等研究了煤體的各向異性和非均質(zhì)性對(duì)瓦斯儲(chǔ)存的影響機(jī)制,實(shí)驗(yàn)結(jié)果表明煤的非均質(zhì)性和各向異性對(duì)氣體擴(kuò)散行為有顯著影響[17]。岳高偉等依據(jù)各向異性煤體的滲透率測(cè)試結(jié)果進(jìn)行了氣-固耦合模擬,得出了平行層理方向抽采效果優(yōu)于垂直層理的結(jié)論[18]。Wang等考慮了煤的各向異性,探究了煤層的各向異性在模擬水平滲透率隨壓力變化中的作用[19]。
目前對(duì)于水力割縫煤層瓦斯?jié)B流規(guī)律的研究大都是基于煤體各向同性來開展,且對(duì)于割縫鉆孔周圍瓦斯流場(chǎng)的時(shí)空演化規(guī)律研究較少,因此研究煤體的非均質(zhì)性和各向異性對(duì)水力割縫鉆孔瓦斯抽采影響規(guī)律是非常必要的。文中基于煤體非均質(zhì)性和各向異性,建立了考慮煤體應(yīng)力場(chǎng)和滲流場(chǎng)的單孔與多孔抽采模型,模擬分析了抽采鉆孔周圍瓦斯流場(chǎng)的時(shí)空演化規(guī)律。該研究結(jié)果可為現(xiàn)場(chǎng)抽采鉆孔布置和提高煤層瓦斯的回采率提供依據(jù)。
煤體瓦斯抽采的流固耦合模型在如下的假設(shè)條件下建立。
1)煤體變形處于線彈性變形階段,服從廣義胡克定律;
2)煤層瓦斯為理想氣體,且在裂隙中的流動(dòng)符合達(dá)西定律;
3)煤體中的瓦斯流動(dòng)場(chǎng)穩(wěn)定變化不大,按照等溫流動(dòng)處理;
4)含瓦斯煤體是由煤基質(zhì)與裂隙2部分組成,且具有非均質(zhì)性和各向異性。
煤是由多種礦物顆粒、膠結(jié)物等組成的混合體,煤體經(jīng)過復(fù)雜的地質(zhì)演變和構(gòu)造運(yùn)動(dòng)形成結(jié)構(gòu)十分復(fù)雜的固體材料[20]。由此可知,煤體的組成單元性質(zhì)在空間上的分布往往是非均質(zhì)的。為了表征其非均質(zhì)性,假定煤體的細(xì)觀單元的力學(xué)性質(zhì)服從Weibull分布[21],該分布可按如下分布密度函數(shù)來定義
(1)
式中u為滿足Weibull分布的參數(shù)數(shù)值(如彈性模量、泊松比等);u0為一個(gè)與所有單元參數(shù)平均值有關(guān)的參數(shù);m為Weibull分布的形狀參數(shù)。
圖1為不同均質(zhì)度系數(shù)下細(xì)觀單元力學(xué)性質(zhì)在煤體空間中的分布情況。根據(jù)Weibull分布的基本性質(zhì),m越大,材料單元的均質(zhì)性越好,反之表明單元屬性越離散。因此,u0和m被稱作材料的分布參數(shù)。采用式(1),可在數(shù)值模擬中產(chǎn)生煤體材料的非均質(zhì)性參數(shù),考慮現(xiàn)場(chǎng)煤體的實(shí)際力學(xué)性質(zhì),在前人研究基礎(chǔ)上[20-21],綜合確定m值取2.0.
圖1 基于Weibull分布的煤體細(xì)觀單元力學(xué)性質(zhì)分布Fig.1 Distribution of mechanical properties of meso-cells in coal based on Weibull distribution
含瓦斯煤體變形的控制方程[22]
(2)
式中G為煤體的剪切模量,MPa;ui,kk和uk,ki分別為不同方向的位移分量,m;v為煤體的泊松比;α為Biot系數(shù);pi為瓦斯壓力,MPa;K為煤的體積模量,MPa;εs,j為煤應(yīng)變張量;fi為煤體的體積力,MPa.
實(shí)驗(yàn)驗(yàn)證了煤體吸附和解吸瓦斯的應(yīng)變符合朗格繆爾型曲線[23-24],體積應(yīng)變可用Langmuir方程計(jì)算
(3)
式中εs為體積應(yīng)變;εL為L(zhǎng)angmuir體積應(yīng)變常數(shù);p為瓦斯壓力,MPa;PL為體積應(yīng)變等于0.5εL的孔隙壓力,MPa.
瓦斯?jié)B流符合質(zhì)量守恒方程,其質(zhì)量平衡方程可以定義為
(4)
式中m為單位氣體含量,kg·m-3;t為時(shí)間,s;ρg為瓦斯密度,kg·m-3;qg為達(dá)西速度向量,m·s-1;Qs為氣源,kg·m3·s-1.
單位體積煤體裂隙中的瓦斯含量由Langmuir方程計(jì)算表示
(5)
式中φ為煤體孔隙率;ρga為標(biāo)準(zhǔn)條件下的氣體密度,kg·m-3;ρc為煤體密度,kg·m-3;VL為L(zhǎng)angmuir體積常數(shù);PL為L(zhǎng)angmuir壓力常數(shù)。
根據(jù)達(dá)西定律,煤層瓦斯流動(dòng)速度qg可以表示為
(6)
式中μ為動(dòng)力黏度,Pa·s;k為滲透率,m2.
考慮煤體應(yīng)變,煤體孔隙率φ可以表示為[25]
(7)
由Kozeny-Carman方程可知[26],煤體的滲透率與孔隙率的關(guān)系為
(8)
式中k0為煤體初始滲透率,m2.
聯(lián)立式(4)(5)(6)計(jì)算化簡(jiǎn),可得煤體瓦斯流動(dòng)的控制方程
(9)
式中pa為大氣壓力,101.325 kPa.
在COMSOL Multiphysics軟件中根據(jù)應(yīng)力場(chǎng)、滲流場(chǎng)相互耦合作用建立2種瓦斯抽采模型:?jiǎn)慰淄咚钩椴赡P秃投嗫淄咚钩椴赡P?。模型的尺寸分別為20 m×30 m和90 m×120 m.根據(jù)現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)設(shè)置了邊界條件,如圖2所示,模型底部為固定約束,煤體左右邊界為輥支撐,模型上部巖層壓力10 MPa.煤體原始瓦斯壓力P0=2 MPa,鉆孔抽采負(fù)壓P=15 kPa.由于水力割縫鉆孔周圍存在擾動(dòng)破裂區(qū)域,因此煤層鉆孔半徑設(shè)置為0.5 m.從表1可以看出,煤體的平均彈性模量、泊松比和密度等基本參數(shù)來自現(xiàn)場(chǎng)實(shí)測(cè),其余部分來自前人研究[27]。分別開展了2種情形下的瓦斯流場(chǎng)演化模擬。
情形一:?jiǎn)慰壮椴蓷l件下,考慮煤層的非均質(zhì)性和各向異性,分析抽采過程煤體關(guān)鍵參數(shù)的時(shí)空演化規(guī)律,幾何模型如圖2(a)所示。
情形二:多孔抽采條件下,考慮煤層的非均質(zhì)性和各向異性,分析抽采過程抽采率的變化規(guī)律,幾何模型如圖2(b)所示。
圖2 幾何模型圖Fig.2 Geometric model diagram
表1 模型的基本參數(shù)
煤的形成經(jīng)歷了長(zhǎng)期的地質(zhì)演變和構(gòu)造運(yùn)動(dòng),因而煤層的結(jié)構(gòu)十分復(fù)雜。如圖3所示,煤層的層理結(jié)構(gòu)存在結(jié)構(gòu)異性,這種結(jié)構(gòu)異性特征導(dǎo)致了煤層的滲透率表現(xiàn)出明顯的各向異性,煤層平行層理和垂直層理方向的瓦斯流場(chǎng)存在明顯差異。根據(jù)前人的實(shí)驗(yàn)研究可知[18],平行層理方向滲透率約是垂直層理方向滲透率的2~4倍,文中數(shù)值模擬選用的煤體平行層理方向初始滲透率為垂直層理方向初始滲透率的3倍。圖4是煤層的彈性模量分布圖,煤體彈性模量的分布符合Weibull分布,這與煤層的非均質(zhì)性比較相符。文中數(shù)值模擬中彈性模量的分布范圍為2.0~3.0 GPa.
圖3 煤層的層理結(jié)構(gòu)Fig.3 Layer structure of coal seam
圖4 煤層的彈性模量分布Fig.4 Elastic modulus distribution of coal seam
圖5給出了距鉆孔半徑2 m煤體滲透率變化。從總體上看,整個(gè)極坐標(biāo)中煤體的滲透率曲線呈現(xiàn)為“不規(guī)則橢圓形”。這是因?yàn)?,在煤體各向異性影響下,平行層理方向煤體的滲透率大于垂直層理方向煤體的滲透率。在煤體非均質(zhì)性影響下,滲透率曲線的變化不規(guī)則。抽采初期,鉆孔周圍煤體的各向滲透率較小,隨著抽采時(shí)間增加,煤體的各向滲透率均逐漸增大。這是因?yàn)橥咚箟毫Φ慕档蛯?dǎo)致了煤基質(zhì)的收縮,引起了煤體裂隙寬度的增大,進(jìn)而滲透率逐漸增大。
圖5 鉆孔周圍滲透率變化Fig.5 Permeability change diagram around the borehole
圖6為煤層瓦斯壓力沿著平行層理和垂直層理方向的空間演化規(guī)律,整體上看距離抽采鉆孔越遠(yuǎn),煤層瓦斯壓力越大,最終趨于穩(wěn)定。這是因?yàn)殡S著瓦斯抽采的進(jìn)行,鉆孔周圍瓦斯容易抽出,距離鉆孔較遠(yuǎn)的瓦斯較難抽出。隨著抽采時(shí)間的增加,平行層理方向和垂直層理方向瓦斯壓力均逐漸減小。對(duì)比平行層理方向圖6(a)和垂直層理方向圖6(b),可以看出距離鉆孔相同的位置(以距離鉆孔2 m位置為例),平行層理方向煤層裂隙瓦斯壓力從2.0 MPa降為1.68 MPa,降幅較大,而垂直層理方向的瓦斯壓力從2.0 MPa降至1.91 MPa,降幅較小。這是因?yàn)槠叫袑永矸较虻臐B透率較大,而垂直層理方向的滲透率較小,所以煤層平行層理方向裂隙的瓦斯更容易被抽離。
圖6 煤層瓦斯壓力的空間演化規(guī)律Fig.6 Spatial evolution law of coal seam gas pressure
圖7是煤層滲透率沿著平行層理和垂直層理方向的空間演化規(guī)律,滲透率曲線表現(xiàn)為不規(guī)“鋸齒形”,這是煤的非均質(zhì)性導(dǎo)致的。距離鉆孔越遠(yuǎn),煤體的滲透率越低,并趨于穩(wěn)定。這是由于煤層在地應(yīng)力影響下,鉆孔周圍存在明顯應(yīng)力集中區(qū),水力割縫鉆孔周圍擾動(dòng)區(qū)域煤體破碎損傷嚴(yán)重,而距離鉆孔中心越遠(yuǎn)的煤體損傷越小,所以距離抽采孔越近滲透率越大,距離較遠(yuǎn)的煤體滲透率逐漸減小。抽采時(shí)間越久,煤體的滲透率越大。這是由于煤層瓦斯含量不斷下降,煤基質(zhì)收縮,煤體裂隙寬度變大,滲透率增大。對(duì)比平行層理方向圖7(a)和垂直層理方向圖7(b),可知平行層理方向煤體滲透率較大,垂直層理方向煤體滲透率較小,且平行層理方向的滲透率在距離鉆孔5 m以后趨于穩(wěn)定,而垂直層理方向的滲透率在距離鉆孔3 m以后趨于穩(wěn)定。這是因?yàn)殂@孔周圍存在應(yīng)力集中區(qū),平行層理方向煤體所受應(yīng)力較大,而垂直層理方向煤體所受應(yīng)力較小,所以平行層理方向的抽采影響范圍較大,垂直層理方向的抽采影響范圍較小。
圖7 煤層滲透率的空間演化規(guī)律Fig.7 Spatial evolution law of coal seam permeability
圖8給出了各測(cè)點(diǎn)瓦斯壓力隨時(shí)間的演化規(guī)律。各個(gè)測(cè)點(diǎn)的瓦斯壓力隨著抽采時(shí)間的增加總體上呈現(xiàn)下降趨勢(shì),測(cè)點(diǎn)瓦斯壓力衰減速度距離鉆孔越近衰減速度越快,即P1>P2>P3>P4,P5>P6>P7>P8.以距離鉆孔1 m的P1和P5兩測(cè)點(diǎn)為例,P1測(cè)點(diǎn)抽采120 d后瓦斯壓力下降為0.4 MPa左右,而P5測(cè)點(diǎn)抽采120 d后瓦斯壓力下降為0.6 MPa,相同抽采時(shí)間下P1測(cè)點(diǎn)的瓦斯衰減幅度大于P5,其他各點(diǎn)情況相似,這是因?yàn)樵诔椴蛇^程中煤層內(nèi)瓦斯總含量是下降趨勢(shì),而平行層理方向煤層滲透率大于垂直層理方向的滲透率,距離抽采鉆孔相同的位置,平行層理方向煤體裂隙的瓦斯更容易被抽離出煤層。
圖8 各測(cè)點(diǎn)瓦斯壓力的時(shí)間演化規(guī)律Fig.8 Time evolution law of gas pressure at each measuring point
圖9 各測(cè)點(diǎn)滲透率的時(shí)間演化規(guī)律Fig.9 Time evolution law of permeability of each measuring point
圖9是各個(gè)測(cè)點(diǎn)滲透率隨著時(shí)間的演化規(guī)律,各點(diǎn)的滲透率隨著抽采時(shí)間的增加都有明顯的上升,并且瓦斯抽采初期煤體滲透率上升速度較快,抽采后期煤體滲透率上升速度較為緩慢穩(wěn)定。這是因?yàn)槌椴沙跗诿后w裂隙內(nèi)瓦斯迅速被抽離,煤基質(zhì)收縮,進(jìn)而導(dǎo)致煤體裂隙寬度增大,煤體滲透率快速增加,經(jīng)過一段時(shí)間的抽采,煤層瓦斯含量下降,瓦斯?jié)B流速度減緩,滲透率緩慢增加。將平行層理方向的4個(gè)測(cè)點(diǎn)和垂直層理方向的4個(gè)測(cè)點(diǎn)進(jìn)行對(duì)比,可以發(fā)現(xiàn)平行層理方向測(cè)點(diǎn)的滲透率都大于垂直層理方向的滲透率,并且測(cè)點(diǎn)P3的滲透率略小于測(cè)點(diǎn)P4,而測(cè)點(diǎn)P7的滲透率略大于測(cè)點(diǎn)P8.這是由于煤體具有非均質(zhì)的特性,選取的測(cè)點(diǎn)P3和P4兩點(diǎn)的煤體力學(xué)性質(zhì)不同,因此煤層存在滲透率變化異常點(diǎn)。
根據(jù)單孔瓦斯抽采模擬結(jié)果以及楊柳礦現(xiàn)場(chǎng)實(shí)際情況建立了多孔瓦斯抽采模型,為驗(yàn)證數(shù)學(xué)模型和物理模型的合理性,將模擬中多孔抽采下的煤體抽采率與楊柳礦現(xiàn)場(chǎng)實(shí)驗(yàn)所測(cè)抽采率進(jìn)行對(duì)比分析。
筆者早前對(duì)楊柳礦4#鉆場(chǎng)進(jìn)行了現(xiàn)場(chǎng)實(shí)驗(yàn),4#鉆場(chǎng)為6列7排共42個(gè)水力割縫鉆孔,孔間距為6 m.為了分析抽采鉆孔的卸壓增透效果,監(jiān)測(cè)記錄了4#鉆場(chǎng)瓦斯抽采純流量。根據(jù)楊柳礦4#鉆場(chǎng)的鉆孔數(shù)量和布置,在COMSOL Multiphysics軟件中建立了相同的幾何物理模型,根據(jù)表1和現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)設(shè)置了相同參數(shù),經(jīng)過模擬計(jì)算得到整個(gè)鉆場(chǎng)的抽采率。圖10顯示了多孔抽采瓦斯壓力分布云圖,從圖中可以看出鉆孔群周圍瓦斯壓力隨抽采時(shí)間增加逐漸降低,抽采60d后抽采影響范圍呈現(xiàn)為“橢圓形”,平行層理的瓦斯抽采效果優(yōu)于垂直層理。圖11顯示了數(shù)值模擬與現(xiàn)場(chǎng)實(shí)驗(yàn)抽采率的對(duì)比分析。從圖11可以看出,前十天多孔抽采鉆孔模型模擬的抽采率結(jié)果與現(xiàn)場(chǎng)監(jiān)測(cè)記錄的結(jié)果誤差較大,為10%左右,這是因?yàn)槟M條件較現(xiàn)場(chǎng)條件相對(duì)理想,且現(xiàn)場(chǎng)測(cè)試存在誤差。而十天之后平均誤差減小為6%左右,從整體上看,模型模擬的抽采率結(jié)果能夠預(yù)測(cè)整體規(guī)律,具有良好的擬合性。結(jié)果表明所建立的模型是正確的、合理的,驗(yàn)證了該模型的可靠性以及工程適用性,可以對(duì)今后現(xiàn)場(chǎng)鉆孔的布置和施工提供相應(yīng)的理論指導(dǎo)。
圖10 多孔協(xié)同瓦斯抽采過程中煤體瓦斯壓力分布云圖Fig.10 Gas pressure distribution in coal body during multi-hole gas drainage
圖11 數(shù)值模擬與現(xiàn)場(chǎng)測(cè)試抽采率的對(duì)比分析Fig.11 Comparative analysis of numerical simulation and field test gas production rate
1)受煤體各向異性的影響,等間距測(cè)點(diǎn)的瓦斯壓力降幅不同,平行層理方向測(cè)點(diǎn)的瓦斯壓力下降速度明顯快于垂直層理方向。隨著瓦斯抽采時(shí)間的增加,瓦斯抽采影響范圍逐漸增大,瓦斯壓力云圖呈現(xiàn)“橢圓形”,平行層理方向的瓦斯抽采效果優(yōu)于垂直層理方向。
2)受煤體非均質(zhì)性的影響,等間距測(cè)點(diǎn)的滲透率曲線呈現(xiàn)“鋸齒形”,導(dǎo)致在滲透率的空間演化規(guī)律中會(huì)出現(xiàn)異常點(diǎn)情況。隨著瓦斯抽采時(shí)間的增加,煤層瓦斯含量逐漸降低,煤基質(zhì)收縮,煤體裂隙寬度增加,煤層的滲透率逐漸增大。
3)結(jié)合楊柳礦4#鉆場(chǎng)實(shí)際監(jiān)測(cè)數(shù)據(jù)和多孔協(xié)同瓦斯抽采的數(shù)值模擬結(jié)果,可以得到現(xiàn)場(chǎng)實(shí)測(cè)的抽采率與模擬得到的抽采率相互吻合,從而驗(yàn)證了所建模型的合理性以及工程適用性。