張　波，謝雄剛，劉洋成

(1.貴州大學 礦業(yè)學院，貴州 貴陽 550025；2.喀斯特地區(qū)優(yōu)勢礦產(chǎn)資源高效利用國家地方聯(lián)合工程實驗室，貴州 貴陽 550025；3.貴州省非金屬礦產(chǎn)資源綜合利用重點實驗室，貴州 貴陽 550025)

0　引言

我國煤層地質(zhì)條件復雜多變，煤層瓦斯含量大，據(jù)統(tǒng)計，全國有60%以上的煤礦屬于高瓦斯或煤與突出礦井[1]。隨著煤層開采深度的增加，溫度、地應(yīng)力和瓦斯含量增大，工作面前方瓦斯涌出量較多，容易導致煤與瓦斯突出事故[2]。且我國煤礦低滲透氣性煤層較多，鉆孔抽采瓦斯可以降低瓦斯含量，但較低的滲透率導致瓦斯抽采效果不佳，增加了煤礦開采的難度。因此，研究低滲透氣性煤層在溫度和地應(yīng)力條件下瓦斯?jié)B流規(guī)律顯得尤為重要。

在研究煤層中瓦斯的運移規(guī)律時，學者們分析了煤與瓦斯的作用機理，認為煤體有效應(yīng)力改變、瓦斯壓力變化和溫度升降會造成煤體變形，導致煤層孔隙率和滲透率發(fā)生改變，分別提出了各種數(shù)學模型來適應(yīng)不同的地質(zhì)條件。尹光志等[3-5]通過引入瓦斯吸附膨脹應(yīng)力，建立了能描述含瓦斯煤流固耦合作用下的骨架可變形性和氣體可壓縮性的流固耦合模型，得出鉆孔抽采瓦斯的有效半徑會隨著時間推移逐漸變大，最后趨于一個定值；鮮學福，周軍平等[6-9]基于熱力學方程建立了含瓦斯煤由吸附瓦斯引起的煤巖膨脹應(yīng)變的理論模型，得到了吸附瓦斯引起的吸附膨脹變形；許江，李波波等[10]對不同溫度條件下，型煤在瓦斯壓力升降過程中瓦斯?jié)B流特性進行了研究，得出在升降壓作用下，瓦斯?jié)B透率的變化規(guī)律；陶云奇[11]在流固耦合研究的基礎(chǔ)上，提出了溫度、地應(yīng)力和瓦斯壓力共同作用下的熱流固三場耦合，并建立了相對應(yīng)的數(shù)學模型，指出在深部開采時，溫度不可忽略；范超軍，李勝等[12]分析了地下水對瓦斯?jié)B透率的影響，提出了瓦斯和水的相對滲透率來表征瓦斯-水兩相流滲透率相互作用關(guān)系，建立了深部煤層瓦斯-水兩相熱-流-固耦合模型。然而，以上研究忽略了受采動影響下工作面前方溫度、瓦斯壓力和地應(yīng)力的相互作用。本文通過模擬熱流固耦合作用下瓦斯?jié)B流過程，建立了工作面前方瓦斯運移規(guī)律數(shù)學模型。

1　深部瓦斯?jié)B流模型

1.1　瓦斯運移模型

煤層內(nèi)的瓦斯存在吸附和游離兩種狀態(tài)，吸附瓦斯占比80%～90%左右[13]，工作面前方受采動影響，在卸壓區(qū)形成較多新鮮裂隙，大量吸附瓦斯轉(zhuǎn)化為游離狀態(tài)；由于在卸壓區(qū)和應(yīng)力集中區(qū)存在較大瓦斯壓力梯度，煤層內(nèi)部瓦斯?jié)B流到工作面，充滿整個工作面前方，造成工作面瓦斯涌出量劇增，在應(yīng)力集中區(qū)容易發(fā)生瓦斯突出事故，如圖1所示。

圖1　工作面前方瓦斯運移模型Fig.1　Methane transport model in front of working face

1.2　基本假設(shè)

依據(jù)瓦斯賦存條件，提出以下假設(shè)：①煤是一種孔隙-裂隙雙重介質(zhì)的彈性體；②煤體被單項瓦斯所飽和；③游離瓦斯的慣性力、瓦斯的體積力忽略不計；④飽和孔隙-裂隙介質(zhì)的體積變形包括骨架變形和孔隙、裂隙變形[14];⑤瓦斯在煤層中的滲流規(guī)律符合達西定律；⑥煤體中吸附游離瓦斯分別服從修正Langmuir方程和真實氣體狀態(tài)方程[15]；⑦煤體變形是微量的，處于線彈性變形階段，服從廣義虎克定律；⑧應(yīng)力應(yīng)變符號與彈性力學相同，壓應(yīng)力為負，拉應(yīng)力為正，位移沿坐標軸正方向為正，反方向為負。

1.3　瓦斯運移模型控制方程

1.3.1應(yīng)力場方程

煤體應(yīng)變由3部分組成，包括瓦斯吸附膨脹(解吸收縮)引起的煤基質(zhì)變形[16]，地應(yīng)力對煤體的壓迫變形，溫度變化導致的熱膨脹變形。三者共同作用相互耦合引起煤體發(fā)生變形，具體作用效果如圖2所示。

圖2　煤體受力變形Fig.2　Coal body force deformation

從宏觀上看，煤顆粒排列疏松，在一定外應(yīng)力作用下，一方面使煤顆粒排列更加緊密，另一方面壓縮煤顆粒改變顆粒體積形狀，使顆粒體積變小，進而導致孔隙率和滲透率變??；微觀形態(tài)上，煤基質(zhì)體積改變是由地應(yīng)力壓縮導致的地應(yīng)力應(yīng)變，瓦斯吸附膨脹(解吸收縮)引起的吸附膨脹(解吸收縮)應(yīng)變，溫度改變造成的熱膨脹應(yīng)變[17]。結(jié)合彈性力學中平衡微分方程和廣義虎克定律，可以推出煤體受力變形的應(yīng)力應(yīng)變本構(gòu)方程。

瓦斯壓力引起的吸附膨脹應(yīng)變：

(1)

式中：KY為體積壓縮系數(shù)，MPa-1；p為瓦斯壓力，Pa；p0為初始狀態(tài)下瓦斯壓力，Pa；ρ為煤的視密度，t/m3；R為普適氣體常數(shù)，取8.3143；T為絕對溫度，K；Vm=22.4×10-3m3/mol，為氣體摩爾體積；a，b為煤的吸附常數(shù)。

熱膨脹引起的應(yīng)變：

(2)

地應(yīng)力引起的應(yīng)變：

(3)

式中：G為拉梅系數(shù)；υ為煤的泊松比；α為煤有效體積變形。

含瓦斯煤的總應(yīng)變?yōu)?

ε=εp+εT+εW

(4)

THM耦合本構(gòu)方程：

(5)

通過本構(gòu)方程推出應(yīng)力場方程的張量形式為：

(6)

1.3.2滲流場方程

瓦斯在煤中的流動規(guī)律，在較大裂隙中符合菲克擴散定律，在小孔隙中遵循達西滲流定律，本文只研究煤層氣在煤體中的滲流特性，因此，選用達西定律來計算滲透率，即：

(7)

(8)

(9)

式中：q為瓦斯?jié)B流速度，m/s；▽p為瓦斯壓力梯度，Pa/m;μ為氣體運動粘度，Pa·s。

修正的langmuir方程和真實氣體狀態(tài)方程[18]為：

(10)

(11)

綜合以上公式可以得出滲流場方程為：

(12)

式中：e為煤總體積形變；m為滑脫因子。

1.3.3溫度場方程

在瓦斯流動過程中，煤層中瓦斯的吸附解吸和滲流都具有熱效應(yīng)，是一個非等溫過程,含瓦斯煤在外力作用下發(fā)生的變形同樣產(chǎn)生熱效應(yīng)[19]，因此滲流場、應(yīng)力場和溫度場存在著相互耦合作用。

溫度場耦合方程為：

(13)

將方程(6)，(12)和(13)聯(lián)立，即得到了深部煤層瓦斯?jié)B流的THM模型，利用COMSOL Multiphysics的PDE方程模塊、固體力學模塊和熱方程模塊進行數(shù)值模擬。

2　工作面前方瓦斯?jié)B流數(shù)值模擬

2.1　幾何模型和邊界條件

該模型選取貴州省盤州市松和煤礦15#煤層12150采煤工作面為工程實例，為了計算方便，沿回采面走向取單位厚度剖面作為計算的平面模型，模型長度和厚度分別為20 m和3 m。瓦斯?jié)B流場邊界條件：瓦斯在煤層中流動的初始瓦斯壓力為p=1.27 MPa，工作面瓦斯壓力為p0=0.1 MPa，煤層頂?shù)装鍨閹r層，瓦斯流量為0。應(yīng)力場邊界條件：在模型上下邊界和左右邊界位移為0，煤層受上覆巖層重力作用，上邊界有均勻載荷，同時煤層自身重力作用一并考慮。溫度場邊界條件：煤層各邊界設(shè)置溫度為293 K，由于溫度場的存在，煤層內(nèi)部存在溫度、瓦斯壓力和地應(yīng)力耦合作用，假設(shè)模型邊界和外部不存在熱量傳遞。計算過程中的基本參數(shù)如表1所示。

表1　數(shù)值模擬基本參數(shù)

2.2　工作面前方瓦斯壓力變化規(guī)律

由于采動應(yīng)力影響，工作面前方卸壓區(qū)瓦斯裂隙較為發(fā)育，工作面煤壁瓦斯壓力為0.1 MPa，與初始瓦斯壓力1.27 MPa存在很大瓦斯壓力梯度，應(yīng)力的卸除，促進煤體內(nèi)部瓦斯解吸釋放，瓦斯從煤層內(nèi)部運移到工作面，導致工作面瓦斯涌出量不斷增加，煤層內(nèi)部瓦斯壓力逐漸降低。根據(jù)數(shù)值模擬得到不同時間點(5,10,15,20,25,30 h)的瓦斯壓力分布特征(見圖3～4)，工作面前方8 m內(nèi)，從卸壓區(qū)到應(yīng)力集中區(qū)的過渡段，煤層瓦斯壓力下降梯度較大，煤體解吸瓦斯沿裂隙快速釋放，垂直應(yīng)力的解除，釋放很大能量，解釋了該過渡段容易發(fā)生瓦斯突出的原因，也說明了應(yīng)力集中區(qū)易發(fā)生煤與瓦斯突出。

圖3　各時間段瓦斯壓力分布Fig.3　Distribution of gas pressure distribution in each time period

圖4　工作面前方瓦斯壓力變化Fig.4　Gas pressure changes in front of work surface

如圖4可知，工作面前方20 m內(nèi)，在0～10 m處瓦斯壓力急劇升高，新鮮裂隙較多，瓦斯壓力下降明顯；從10～20 m 處瓦斯壓力下降緩慢，位于應(yīng)力集中區(qū)；20 m以后瓦斯壓力為原始瓦斯壓力。工作面前方5 m處，開采時間為5，10，15，20，25，30 h 的瓦斯壓力分別為1.23，1.0，0.91，0.85，0.8，0.78 MPa。從應(yīng)力集中區(qū)到卸壓區(qū)，不同時間節(jié)點下瓦斯壓力下降趨勢隨開采時間的增加，瓦斯壓力梯度逐漸減小，吸附瓦斯含量降低，瓦斯總含量降低，瓦斯運移速率相對減慢；在工作面前方同一位置，瓦斯壓力隨著時間的增加逐漸降低，但降低速率減慢，最后，吸附瓦斯和游離瓦斯達到新的動態(tài)平衡，瓦斯壓力趨于一個定值。圖5為工作面前方1 m處不同時間節(jié)點的瓦斯壓力分布，也符合以上分析規(guī)律。

圖5　工作面前方1m處不同時間節(jié)點瓦斯壓力分布Fig.5　Distribution of gas pressure at different time points at 1m ahead of working face

3　壓力變化對滲透率的影響

孔隙率和滲透率由于有效壓力、瓦斯壓力和溫度作用，從煤層內(nèi)部到工作面，先降低后逐漸增大，最后趨于穩(wěn)定。有效應(yīng)力較高，壓縮煤體，導致煤顆粒排列緊密，造成孔隙率相對較低，進而瓦斯?jié)B透率較低，隨著采煤活動的進行，內(nèi)部裂隙逐漸增大孔隙率變大，瓦斯?jié)B透率逐漸升高。

圖6為滲透率隨工作面長度變化規(guī)律。如圖6所示，煤層深部未受采動影響的原始應(yīng)力區(qū)，瓦斯?jié)B透率為8×10-15m2，在應(yīng)力集中區(qū)，瓦斯壓力和有效應(yīng)力較高，壓縮煤體，導致煤顆粒排列緊密，滲流通道變窄，造成孔隙率相對較低，滲透率降低。在卸壓區(qū)，煤體體積形變逐漸變大，產(chǎn)生了很多新裂隙，發(fā)生擴容，滲流通道貫通，導致滲透率急劇增加，因此在應(yīng)力最大處形成了煤層滲透率最低點，并且隨著時間的推移，滲透率最低點逐步遠離工作面。工作面前方20 m處煤層滲透率趨于穩(wěn)定為7.8×10-15m2，比初始滲透率偏低，這是由于深部煤層溫度較高，使煤體顆粒變大，壓縮孔隙和裂隙空間，導致煤層滲透率降低。

圖6　滲透率隨工作面長度變化規(guī)律Fig.6　Permeability with the change of working face length

4　溫度變化對瓦斯?jié)B流的影響

隨著煤層深度加深，煤體溫度有所上升，煤體溫度升高時，一方面，瓦斯分子熱運動加劇，吸附瓦斯獲取一部分內(nèi)能后，脫離煤基質(zhì)表面成為游離瓦斯；另一方面，煤體受熱后向內(nèi)膨脹，使?jié)B流通道變窄，減小滲流空間，導致滲透率減小。

圖7為滲透率隨溫度變化規(guī)律。由圖7可以看出隨著溫度從293～296 K變化中，滲透率從7.8×10-15降低到7.74×10-15，煤層滲透率逐漸減小，距離工作面越遠處滲透率減小幅度越大。由此可知，在采煤工作面前方，煤層滲透率受溫度影響，雖然溫度升高后瓦斯熱運動加劇，有促進瓦斯?jié)B透率的趨勢，但由于工作面前方外應(yīng)力較大，煤體受熱膨脹應(yīng)力小于有效應(yīng)力，導致煤體內(nèi)膨脹，滲流空間減小，造成滲透率降低。綜上可知，采煤工作面前方瓦斯?jié)B透率受瓦斯壓力梯度、煤層溫度和有效應(yīng)力相互影響，開采初期，瓦斯壓力梯度較大，由于采動影響，產(chǎn)生大量裂隙和通道，瓦斯?jié)B透率較大。隨著開采進行，當溫度升高時，煤體受熱膨脹，但由于熱膨脹應(yīng)力小于有效應(yīng)力，煤體內(nèi)膨脹，滲流通道減小，滲透率降低。綜合效果看，在同一位置，隨著時間推移滲透率是逐漸減小的。

圖7　滲透率隨溫度變化規(guī)律Fig.7　Permeability with temperature variation

5　結(jié)論

1)建立了工作面前方瓦斯?jié)B流熱流固耦合數(shù)學模型，通過數(shù)值模擬分析了煤層瓦斯壓力和滲透率變化規(guī)律，給出了受工作面前方“三帶”影響下，煤層瓦斯壓力分布情況。從卸壓區(qū)到應(yīng)力集中區(qū)的過渡段，煤層瓦斯壓力下降梯度較大，煤體解吸瓦斯沿裂隙快速釋放，垂直應(yīng)力的解除，釋放很大能量，解釋了該過渡段容易發(fā)生瓦斯突出的原因。

2)在采動影響下，在應(yīng)力集中區(qū)，瓦斯壓力和有效應(yīng)力較高，壓縮煤體，導致煤顆粒排列緊密，滲流通道變窄，造成孔隙率相對較低，滲透率降低。在卸壓區(qū)，煤體體積形變逐漸變大，產(chǎn)生了很多新裂隙，發(fā)生擴容，滲流通道貫通，導致滲透率急劇增加，因此在應(yīng)力最大處形成了煤層滲透率最低點為7.64×10-15m2。并且隨著時間的推移，滲透率最低點逐步遠離工作面。

3)在采煤工作面前方，雖然溫度升高后瓦斯熱運動加劇，有促進瓦斯?jié)B透率的趨勢，但由于工作面前方有效應(yīng)力較大，煤體受熱膨脹應(yīng)力小于有效應(yīng)力，導致煤體內(nèi)膨脹，滲流空間減小，造成滲透率降低。

[1]張子敏.瓦斯地質(zhì)學[M].徐州：中國礦業(yè)大學出版社，2009：2-3.

[2]陳亮.工作面前方煤體變形破壞和滲透率演化及其應(yīng)用研究[D].北京：中國礦業(yè)大學，2016.

[3]尹光志，李銘輝，李生舟，等.基于含瓦斯煤巖固氣耦合模型的鉆孔抽采瓦斯三維數(shù)值模擬[J].煤炭學報，2013，38(4)：469-472.

YIN Guangzhi,LI Minghui,LI Shengzhou,et al.Three-dimensional numerical simulation of gas drainage by drilling based on gas-solid coupling model of coal and gas [J].Journal of China Coal Society,2013,38 (4):469-472.

[4]尹光志，李文璞，李銘輝，等.加卸載條件下原煤滲透率與有效應(yīng)力的規(guī)律[J].煤炭學報，2014，39(8)：1497-1503.

YIN Guangzhi,LI Wenpu,LI Minghui,et al.Law of permeability and effective stress of raw coal under loading and unloading [J].Journal of China Coal Society,2014,39(8):1497-1503.

[5]尹光志，李文璞，李銘輝，等.不同加卸載條件下含瓦斯煤力學特性試驗研究[J].巖石力學與工程學報，2013,32(5)：893-899.

YIN Guangzhi,LI Wenpu,LI Minghui,et al.Experimental study on mechanical properties of gas coal under different loading and unloading conditions [J].Chinese Journal of Rock Mechanics and Engineering,2013,32(5):893-899.

[6]周軍平，鮮學福，李曉紅，等.多孔介質(zhì)中不混溶兩相滲流的本構(gòu)關(guān)系[J].四川大學學報(工程科學版)，2010，42(3)：107-111.

ZHOU Junping,XIAN Xuefu,LI Xiaohong,et al.Constitutive relationship of immiscible two-phase seepage in porous media [J].Journal of Sichuan University(Engineering Science Edition),2010,42 (3):107-111.

[7]董賀，王佰順，李春亭.穿層預抽鉆孔傾角與煤層氣抽采效果關(guān)系的研究[J].中國安全生產(chǎn)科學技術(shù)，2015，11(2) :99-104.

DONG He,WANG Baishun,LI Chunting.Research on the relationship between pre-pumping borehole dip angle and coalbed methane drainage effect [J].Journal of Safety Science and Technology,2015,11(2):99-104.

[8]王臣，鮮學福，周軍平，等.含不同氣體煤巖全應(yīng)力應(yīng)變滲透特性試驗研究[J].地下空間與工程學報，2013，9(3)：492-495.

WANG Chen,XIAN Xuefu,ZHOU Junping,et al.Experimental study of total stress and strain permeability of coal with different gas [J].Journal of Underground Space and Engineering,2013,9(3):492-495.

[9]周軍平，鮮學福，姜永東，等.考慮基質(zhì)收縮效應(yīng)的煤層氣應(yīng)力場-滲流場耦合作用分析[J].巖土力學，2010，31(7)：2318-2322.

ZHOU Junping,XIAN Xuefu,JIANG Yongdong,et al.Coupling effect analysis of stress field-seepage field in coalbed methane considering matrix shrinkage effect [J].Rock and Soil Mechanics,2010,31(7):2318-2322.

[10]許江，張敏，李波波，等.不同溫度條件下氣體壓力升降過程中瓦斯運移規(guī)律的試驗研究[J].巖土力學，2016，37(6)：1580-1585.

XU Jiang,ZHANG Min,LI Bobo,et al.Experimental study on gas migration during gas pressure rise and fall under different temperature conditions [J].Rock and Soil Mechanics,2016,37 (6):1580-1585.

[11]陶云奇.含瓦斯煤THM耦合模型及煤與瓦斯突出模型研究[D].重慶：重慶大學，2009.

[12]范超軍,李勝,羅明坤,等.基于流-固-熱耦合的深部煤層氣抽采數(shù)值模擬[J].煤炭學報,2016,41(12):3076-3085.

FAN Chaojun,LI Sheng,LUO Mingkun,et al.Deep CBM extraction numerical simulation based on hydraulic-mechanical-thermal coupledmodel[J].Journal of China Coal Society,2016,41(12) :3076-3085.

[13]肖偉晶，王曉軍，李士超，等.滲透壓力作用下巖石三軸壓縮過程變形特性分析[J].中國安全生產(chǎn)科學技術(shù)，2017，13(12)：38-42.

XIAO Weijing,WANG Xiaojun,LI Shichao,et al.Analysis on deformation characteristics of rock in triaxial compression process with seepage pressure[J].Journal of Safety Science and Technology,2017,13(12):38-42.

[14]趙陽升.煤體-瓦斯耦合數(shù)學模型及數(shù)值解法[J].巖石力學與工程學報，1994，26(3)：229-239.

ZHAO Yangsheng.Coal-gas coupling mathematical model and numerical solution [J].Journal of Rock Mechanics and Engineering,1994,26 (3):229-239.

[15]李祥春，郭勇義，吳世躍.煤吸附膨脹變形與孔隙率、滲透率關(guān)系的分析[J].太原理工大學學報，2005，36(3)：264-266.

LI Xiangchun,GUO Yongyi,WU Shiyue.Analysis of the relationship between coal adsorption and expansion deformation,porosity and permeability [J].Journal of Taiyuan University of Technology,2005,36 (3):264-266.

[16]孫培德.SUN模型及其應(yīng)用[M].北京：科學出版社，1995.

[17]馬忠.基于熱流固耦合的瓦斯抽采數(shù)值模擬研究[D].合肥：安徽理工大學，2014.

[18]題正義，潘進，田臣，等.7.0 m采高綜采面回采巷道變形破壞規(guī)律研究[J].中國安全生產(chǎn)科學技術(shù)，2016，12(4)：57-61.

TI Zhengyi,PAN Jin,TIAN Chen,et al.Research on deformation and failure of mining gateway with height of 7 m in fully mechanized coal mining face [J].Journal of Safety Scienceand Technology,2016,12 (4):57-61.

[19]陳亮，吳兵，許小凱，等.泥、砂巖交互地層綜放開采覆巖破壞高度的確定[J].采礦與安全工程學報，2017，4(3)：431-436.

CHEN Liang,WU Bing,XU Xiaokai,et al.Determine the failure height of overburden strata in fully mechanized caving mining with interaction between mud and sandstone [J].Mining and Safety Engineering Journal,2017,4(3):431-436.