周宏偉，榮騰龍，牟瑞勇，王路軍，任偉光

(1.中國礦業(yè)大學(xué)(北京) 能源與礦業(yè)學(xué)院，北京 100083; 2.中國礦業(yè)大學(xué)(北京) 力學(xué)與建筑工程學(xué)院，北京 100083; 3.中國礦業(yè)大學(xué)(北京) 煤炭資源與安全開采國家重點(diǎn)實(shí)驗(yàn)室，北京 100083)

煤炭是我國的基礎(chǔ)能源，也是重要的工業(yè)原料，到2030年，煤炭科學(xué)開采預(yù)測(cè)產(chǎn)能總量達(dá)30～35億t[1]。所以，未來一段時(shí)間內(nèi)煤炭生產(chǎn)仍將持續(xù)進(jìn)行，煤礦的安全生產(chǎn)仍需加強(qiáng)保障。煤炭形成過程中伴生有大量的瓦斯氣體，瓦斯既是瓦斯爆炸、煤與瓦斯突出災(zāi)害的誘因，又是一種高效的清潔能源[2-4]。掌握煤層中瓦斯流動(dòng)規(guī)律是減少瓦斯事故、擴(kuò)大瓦斯采集利用的前提，而煤體的滲透率是決定瓦斯流動(dòng)規(guī)律的關(guān)鍵參數(shù)[5]，所以對(duì)煤體滲透率進(jìn)行研究是十分必要的。我國大部分礦區(qū)的煤層滲透率較低，需通過采動(dòng)卸壓、水力壓裂、深孔爆破等手段進(jìn)行人工增透[6-8]，此外煤與瓦斯突出事故多發(fā)于掘進(jìn)工作面[9-11]。因此，建立工作面前方采動(dòng)應(yīng)力下的煤體滲透率模型對(duì)提高瓦斯抽采和減少瓦斯事故具有重要意義。

以往的煤體滲透率模型較多以煤層氣開采為背景[12-14]，同時(shí)基于單軸應(yīng)變條件和上覆載荷不變的假設(shè)建立，而在煤礦采掘活動(dòng)中，煤層的受力特點(diǎn)有別于煤層氣開采。工作面推進(jìn)過程中采煤工作面前方煤體經(jīng)受采動(dòng)應(yīng)力路徑的變化過程，其豎直方向的支承壓力逐漸升高后下降、水平方向的應(yīng)力持續(xù)降低，煤層受力呈現(xiàn)典型的三向應(yīng)力變化特征(圖1(a))。隨著支承壓力的升高以及人工增透活動(dòng)的進(jìn)行，煤體內(nèi)部繼發(fā)出現(xiàn)損傷和破裂現(xiàn)象，煤層內(nèi)裂隙擴(kuò)展貫通造成瓦斯壓力下降，同時(shí)伴隨瓦斯解吸擴(kuò)散的發(fā)生，最終煤體滲透率急劇增長。采動(dòng)應(yīng)力決定了煤體在由原巖應(yīng)力狀態(tài)漸近為側(cè)面臨空過程中的受力狀態(tài)、物理力學(xué)特征、溫度梯度及滲透率變化有別于煤層氣的抽采過程，以往的煤體滲透率模型并不能完全適用于采動(dòng)應(yīng)力條件。

筆者首先概述了煤體滲透率經(jīng)典模型，然后結(jié)合采動(dòng)過程中煤體內(nèi)的力學(xué)變化機(jī)制及滲透率的控制因素提出了采動(dòng)應(yīng)力下煤體滲透模型構(gòu)建中的3個(gè)關(guān)鍵問題，即采動(dòng)煤體各向異性特征、采動(dòng)煤體損傷破裂特征和煤體吸附解吸特征，并將這3個(gè)方面的國內(nèi)外研究成果進(jìn)行了總結(jié)和分析。在此基礎(chǔ)上，將各向異性特征的煤體滲透率模型歸納為由有效應(yīng)力變化表征和幾何參數(shù)變化表征的兩類，又將損傷破裂特征的煤體滲透率模型劃分為本構(gòu)方程中含損傷變量和滲透率表達(dá)式中含損傷變量的兩類。最后，對(duì)采動(dòng)應(yīng)力下煤體滲透率模型的研究進(jìn)展進(jìn)行了介紹并對(duì)當(dāng)前采動(dòng)應(yīng)力下煤體滲透率模型研究的不足和未來的發(fā)展方向進(jìn)行了探討。

1 煤體滲透率經(jīng)典模型概述

為了對(duì)煤體的滲透率演化過程進(jìn)行定量描述，國內(nèi)外學(xué)者建立了大量的滲透率模型，這些研究成果主要集中在煤層氣開采方面，其中一些模型具有較高的代表性，這些經(jīng)典的模型成果在煤體滲透率模型的研究領(lǐng)域中發(fā)揮了重要作用。

GRAY[15]認(rèn)為煤儲(chǔ)層的滲透率與垂直于割理結(jié)構(gòu)的有效水平應(yīng)力大小有直接關(guān)系，較早建立了煤體滲透率演化模型。SEIDLE等[16]將儲(chǔ)層煤體簡(jiǎn)化為“捆綁的火柴棍”幾何結(jié)構(gòu)，得到了煤體滲透率與有效靜水壓力的關(guān)系，從理論角度推導(dǎo)出滲透率變化的指數(shù)函數(shù)形式。PALMER和MANSOORI[17]假設(shè)儲(chǔ)層為單軸應(yīng)變條件，考慮煤層割理壓縮性的變化，基于滲透率k與孔隙率φ的立方關(guān)系，結(jié)合非等溫多孔彈性理論建立了煤層滲透率模型。SHI和DURUCAN[18]從各向同性的彈性本構(gòu)方程出發(fā)，將煤儲(chǔ)層視為連續(xù)介質(zhì)，基于單軸應(yīng)變的假設(shè)條件求出煤儲(chǔ)層水平方向有效應(yīng)力的變化，建立了指數(shù)函數(shù)形式的滲透率模型。CUI和BUSTIN[19]考慮整體平均有效應(yīng)力對(duì)煤體孔隙度和滲透率的影響，利用多孔彈性本構(gòu)方程進(jìn)行推導(dǎo)得到解吸過程中煤體滲透率模型。ZHANG等[20]將煤體視為固體顆粒和孔隙空間組成，假設(shè)煤的孔隙空間同樣存在吸附應(yīng)變，建立了一個(gè)廣義孔隙度的滲透率模型。LIU等[21]認(rèn)為煤體滲透率由基質(zhì)滲透率和裂隙滲透率2部分組成，考慮基質(zhì)和裂隙的相互作用得到合成的滲透率表達(dá)式。經(jīng)典的煤體滲透率模型匯總見表1。

表1煤體滲透率經(jīng)典模型
Table1Classicalpermeabilitymodelsofcoals

文獻(xiàn)來源模型表達(dá)式參數(shù)釋義GRAY[15]σeh-σeh0=-ν1-ν(p-p0)+E(1-ν)ΔεsΔpsΔpsk=1.013×10-0.31σehσeh為水平有效應(yīng)力;ν為泊松比;E為彈性模量;p為氣體壓力;Δεs/Δps為單位等效吸附壓力改變下的煤體應(yīng)變;Δps為等效吸附壓力變化增量;k為煤的滲透率;下標(biāo)0為初始狀態(tài)SEIDLE等[16]kk0=exp[-3Cf(σehy-σehy0)]Cf為割理壓縮系數(shù);σehy為有效的靜水壓力PALMER和MANSOORI[17]φφ0=1+cmφ0(p-p0)+clφ0KM-1 Bp1+Bp-Bp01+Bp0 kk0=φφ0 3M=E(1-ν)(1+ν)(1-2ν),K=E3(1-2ν),cm=1M-KM+G-1 γs,其中γs為顆粒壓縮系數(shù);φ為孔隙率;cl為Langmuir吸附應(yīng)變常數(shù);B為Langmuir吸附平衡常數(shù);G為比例系數(shù)SHI和DURUCAN[18]σeh-σeh0=-ν1-ν(p-p0)+E3(1-ν)εLpp+pL-p0p0+pL kk0=exp[-3Cf(σeh-σeh0)]εL為Langmuir吸附體積應(yīng)變常數(shù);pL為Langmuir吸附壓力常數(shù)CUI和BUSTIN[19]kk0=exp-3Kp[(σm-σm0)-(p-p0)] Kp為孔隙體積模量;σm為平均法向應(yīng)力ZHANG等[20]kk0=11+S(1+S0)+αbφ0(S-S0) ? 3S=εV+(p/Ks)-εs,S0=(p0/Ks)-εLp0/(p0+pL),αb=1-K/Ks,Ks為煤固體顆粒體積模量;K為煤的體積模量;εV為煤的總體積應(yīng)變;εs為煤的吸附應(yīng)變LIU等[21]kk0=km0km0+kf01+Rmpmφm0K 3+kf0km0+kf01+1-Rmφf0(ΔεV-Δεs) ?3km為基質(zhì)滲透率;kf為裂隙滲透率;φm為基質(zhì)孔隙率;φf為裂隙孔隙率;Rm為整體彈性模量損失率;pm為基質(zhì)上的氣體壓力

2 采動(dòng)應(yīng)力下煤體滲透率模型構(gòu)建的關(guān)鍵問題

2.1 采動(dòng)煤體各向異性特征

煤屬于沉積巖，在沉積方向具有明顯的層理結(jié)構(gòu)，而在垂直層理的方向又存在空間狹窄的裂隙結(jié)構(gòu)，這種結(jié)構(gòu)被稱為割理結(jié)構(gòu)。割理可分為面割理和端割理結(jié)構(gòu)，通常將兩種割理結(jié)構(gòu)簡(jiǎn)化為相互正交的形式。因此，煤體具有顯著的各向異性特征(圖1(b))，進(jìn)而煤體在各個(gè)方向上存在不同的滲透率變化趨勢(shì)。除此之外，采動(dòng)應(yīng)力下煤體滲透率模型的構(gòu)建必須把握三向受力的特點(diǎn)，建立可表征豎向加卸載、水平卸載過程的滲透率理論模型。因此，無論是煤體內(nèi)部結(jié)構(gòu)特征還是外部受力條件均要求使用各向異性的滲透率模型對(duì)采動(dòng)應(yīng)力下煤體的滲透率演化過程進(jìn)行描述。筆者根據(jù)滲透率模型的建模角度將現(xiàn)有的各向異性煤體滲透率模型成果歸納為有效應(yīng)力變化和幾何參數(shù)變化兩類進(jìn)行介紹。各向異性特征的煤體滲透率模型匯總見表2。

表2各向異性滲透率模型
Table2Anisotropicpermeabilitymodels

文獻(xiàn)來源模型表達(dá)式參數(shù)釋義CHEN等[22]ki=ki0exp[3Cfi(σei-σei0)] i=x,yCfi為不同方向的割理壓縮系數(shù);σei為不同方向的有效應(yīng)力LIU等[23]ki=ki0exp[-3Cfi(σem-σem0)] i=x,y,zCfi為不同方向的平均割理壓縮系數(shù);σem為有效的平均法向應(yīng)力ZHANG等[24]ki=ki0exp[-αiCpσ(σi-σi0)] i=x,y,zCpσ為孔隙壓縮系數(shù);αi為冪函數(shù)系數(shù);σi為不同方向的應(yīng)力PAN和CONNELL[25]ki=ki0exp[3Cfi(σei-σei0)] i=x,y—GU和CHALATURNYK[26]kiki0=1+ajbjΔεfj 31-Δεfj i=x,yaj為j方向的基質(zhì)單元寬度;bj為j方向的裂隙單元寬度;Δεfj為j方向的裂隙應(yīng)變?cè)隽狂R天然等[27]ki=∑3j=1(bj0+Δbj)312(aj0+Δaj) i=x,y,z—LIU等[28]kiki0=∑i≠j121+3(1-Rm)φf0(Δεtj-13Δεs) ?Δεtj為j方向煤體的總應(yīng)變YANG等[29-31]kxkx0=121+3(1-R2)φ0 ?(Δεt2-13Δεs) 3+ 121+3(1-R3)φ0 ?(Δεt3-13Δεs) 3kyky0=121+3(1-R1)φ0 ?(Δεt1-13Δεs) 3+ 121+3(1-R3)φ0 ?(Δεt3-13Δεs) 3kzkz0=121+3(1-R1)φ0 ?(Δεt1-13Δεs) 3+ 121+3(1-R2)φ0 ?(Δεt2-13Δεs) 3Ri為i方向的模量損失率;Δεti為i方向的煤體總應(yīng)變WANG等[32]kiki0=1+2(1-Rm)φf0Δεej ?3 i,j=x,y,z;i≠jΔεej為j方向的煤體有效體積應(yīng)變WU等[33]kmkm0=1(1+S)φm0[(1+S0)φm0+αb(S-S0)] 3kiki0=∑i≠j121-1φf0+3KfK13αTΔT+13Δεs-1KΔσei ?3Kf為裂隙剛度;αT為熱膨脹系數(shù);T為溫度WANG等[34]ki=ki0 1+1φi0Δσej-vbjkΔσek-vbijΔσeiEbj-Δσej-vmjkΔσek-vmijΔσeiEmj+Δσek-vbkiΔσei-vbjkΔσejEbk-Δσek-vmkiΔσei-vmjkΔσeiEmk-εLjFljppLj+p-Fljp0pLj+p0 -εLkFlkppLk+p-Flk0p0pLk+p0 ? 3i≠j≠kEbi為i方向的煤體彈性模量;Emi為i方向的基質(zhì)彈性模量;εLi為i方向的Langmuir吸附體積應(yīng)變常數(shù);pLi為i方向的Langmuir吸附壓力常數(shù);Fli為i方向的內(nèi)膨脹應(yīng)變系數(shù)

2.1.1 有效應(yīng)力變化表征的煤體滲透率模型

有效應(yīng)力變化表征的煤體滲透率模型都以煤體所受有效應(yīng)力為自變量，在建模過程中較多將煤基質(zhì)和裂隙視為一個(gè)整體，利用彈性本構(gòu)方程對(duì)煤體在坐標(biāo)軸3個(gè)方向的應(yīng)力和應(yīng)變關(guān)系進(jìn)行推導(dǎo)，建立滲透率和有效應(yīng)力之間的關(guān)系。此類模型[22-25]均可以歸納為以下形式

(1)

(2)

2.1.2 幾何參數(shù)變化表征的煤體滲透率模型

幾何參數(shù)變化表征的煤體滲透率模型大多以煤體內(nèi)孔裂隙尺寸、裂隙應(yīng)變和孔隙度為自變量，建模過程基本是將煤體劃分為基質(zhì)和裂隙兩部分，通過計(jì)算應(yīng)力、吸附作用下裂隙尺寸或孔隙度變化，建立煤體裂隙開度、孔隙度和滲透率的關(guān)系。為了考慮煤體各向異性的結(jié)構(gòu)特征(圖1(b))，這些滲透率模型[26-34]構(gòu)建過程中大多將煤體抽象為立方體幾何結(jié)構(gòu)模型(圖1(c))。

一些學(xué)者將氣體在煤層裂隙中的流動(dòng)簡(jiǎn)化為平板間流體流動(dòng)過程，通過POISEUILLE定律得到滲透率與裂隙開度之間的立方關(guān)系(式(3))，并以此開始推導(dǎo),即

(3)

式中，a，b分別為基質(zhì)單元和裂隙單元的寬度。

GU和CHALATURNYK[26]將煤體簡(jiǎn)化為火柴棍模型，而馬天然等[27]將煤體簡(jiǎn)化為立方體模型，分別對(duì)式(4)進(jìn)行變換直接得到了滲透率和裂隙開度的3次方關(guān)系式。LIU等[28]將煤體抽象為立方體幾何結(jié)構(gòu)，認(rèn)為煤體為各向同性的線彈性材料，提出了模量損失率來表征各方向彈性模量的變化過程，建立了不同邊界條件下的各向異性滲透率模型。YANG等[29-31]認(rèn)為工程分析中應(yīng)忽略基質(zhì)滲透率的影響，借此對(duì)LIU等[28]模型加以修正，構(gòu)建出橫觀各向同性的滲透率模型。WANG等[32]從式(4)出發(fā)對(duì)壓縮過程中裂隙孔徑的變化進(jìn)行推導(dǎo)，得出一個(gè)與LIU等[28]模型形式類似的表達(dá)式，同時(shí)定義了基質(zhì)與整體單元應(yīng)變比的新參數(shù)，建立了一個(gè)2維各向異性的滲透率模型。

另外，一些學(xué)者定量表示出三向有效應(yīng)力變化過程中煤體內(nèi)部孔隙度的變化，并將孔隙度變化結(jié)果代入滲透率和孔隙度的立方關(guān)系[17](式(4))建立了相應(yīng)的各向異性滲透率模型。WU等[33]將煤體的物理力學(xué)性質(zhì)視為各向同性，假設(shè)坐標(biāo)軸3個(gè)方向下熱膨脹和吸附作用產(chǎn)生的煤體應(yīng)變相等而有效應(yīng)力不相等，建立了雙孔隙度的煤體各向異性滲透率模型。WANG等[34]基于立方體模型提出了內(nèi)膨脹應(yīng)變系數(shù)以對(duì)吸附過程中煤體基質(zhì)的內(nèi)膨脹量進(jìn)行修正，建立了一個(gè)以三向有效應(yīng)力和氣體壓力為自變量的滲透率模型。

(4)

式中，k為滲透率;φ為孔隙率。

從建模過程上來說，有效應(yīng)力變化和幾何參數(shù)變化的兩個(gè)建模角度可認(rèn)為是互逆的相反過程。有效應(yīng)力的角度是將煤體所受有效應(yīng)力通過孔隙度、應(yīng)變、模量、泊松比、吸附常數(shù)等參數(shù)綜合表示出來，而幾何參數(shù)的角度是將裂隙開度或孔隙度用應(yīng)力、模量、泊松比、吸附常數(shù)等參數(shù)進(jìn)行表達(dá)，進(jìn)而各自代入相應(yīng)的框架基礎(chǔ)公式。此外，根據(jù)各自建模所依據(jù)的框架基礎(chǔ)，有效應(yīng)力角度建模結(jié)果基本為指數(shù)型函數(shù)、幾何參數(shù)角度建模的結(jié)果多數(shù)為3次方的冪函數(shù)。

滲透率模型中如果能將煤體的邊界條件和結(jié)構(gòu)參數(shù)均考慮為各向異性且持續(xù)變化的特征是最為理想和貼近實(shí)際的情況，然而受目前試驗(yàn)條件和測(cè)試手段所限，大量學(xué)者也只能將模型中的部分參數(shù)做出近似和假設(shè)處理。從上述各向異性的模型成果中可以看出，外部邊界的三向應(yīng)力和三向應(yīng)變參數(shù)較易跟蹤和獲取，所以三向應(yīng)力的特征在模型中容易實(shí)現(xiàn)[22-25]。結(jié)構(gòu)異性參數(shù)大多較難準(zhǔn)確把握，比如3個(gè)方向的割理壓縮系數(shù)、彈性模量和泊松比等，因此一些模型中將結(jié)構(gòu)參數(shù)視為各向同性[26-33]。

2.2 采動(dòng)煤體損傷破裂特征

受開采擾動(dòng)的影響，支承壓力的空間分布具有分區(qū)化特征，煤體由原巖應(yīng)力逐漸升高至應(yīng)力峰值后迅速卸壓(圖1(a))，在此過程中煤體由原生初始狀態(tài)漸進(jìn)為損傷破裂狀態(tài)，伴隨煤體內(nèi)微裂紋的擴(kuò)展貫通，滲透率將會(huì)迅速增大，可達(dá)初始值的數(shù)十倍，乃至百倍[35-36]。因此，對(duì)于采場(chǎng)前方煤體必須考慮開采擾動(dòng)引起的損傷破裂對(duì)滲透率的影響。

目前，國內(nèi)外學(xué)者已經(jīng)開始探究損傷破裂影響下滲透率的變化規(guī)律。開采擾動(dòng)下煤體中瓦斯流動(dòng)及滲透率的變化是一個(gè)典型的氣固耦合問題。一種解決此問題的手段就是將應(yīng)力場(chǎng)和滲流場(chǎng)的控制方程進(jìn)行耦合求解。根據(jù)損傷力學(xué)理論可以定義多種形式的損傷變量，損傷變量對(duì)滲透率的控制作用通過耦合方程的形式進(jìn)行傳遞。這種調(diào)控方式下，損傷變量處于固體力學(xué)的本構(gòu)方程之中，而滲透率表達(dá)式內(nèi)并不含損傷變量。另一種處理此問題的方式就是將定義出的損傷變量直接嵌入滲透率模型的表達(dá)式之中。筆者依據(jù)損傷變量所處的位置不同將損傷特征的實(shí)現(xiàn)方式歸納為兩種:一種為損傷變量存在于本構(gòu)方程中，另一種為損傷變量存在于滲透率表達(dá)式中。考慮損傷破裂特征的滲透率模型成果匯總見表3。

表3考慮損傷破裂特征的滲透率模型
Table3Damage-basedpermeabilitymodels

文獻(xiàn)來源模型表達(dá)式參數(shù)釋義YANG等[48]k=k0e-β1(σ1-α1p) D=0ξk0e-β1(σ1-α1p) D>0 k=k0e-β1(σ3-α1p) D=0ξk0e-β1(σ3-α1p) 0

2.2.1 本構(gòu)方程中含損傷變量的煤體滲透率模型

一些學(xué)者基于宏觀唯象的連續(xù)損傷力學(xué)理論建立了煤體損傷本構(gòu)并將其代入固體力學(xué)控制方程，與瓦斯?jié)B流方程聯(lián)立進(jìn)行求解?；贚EMAITRE[37]提出的應(yīng)變等效假設(shè)(式(5))，魏明堯等[38]建立了彈塑性損傷本構(gòu);HU等[39]定義了雙重孔隙結(jié)構(gòu)下煤體的等效彈性損傷柔度張量，構(gòu)建了雙重孔隙介質(zhì)損傷的本構(gòu)模型;WANG等[40-41]將HU[39]定義的等效彈性損傷柔度張量進(jìn)行簡(jiǎn)化后應(yīng)用在單一孔隙率的滲流方程之中,即

(5)

另外，一部分學(xué)者通過細(xì)觀損傷力學(xué)的方法建立了煤體的損傷本構(gòu)，將煤體視為大量代表性體積單元(RVE)組成，通過研究每個(gè)RVE的力學(xué)破壞特征來反映煤體內(nèi)損傷演化過程。ZHU和TANG[42-43]考慮巖石的非均質(zhì)性，假設(shè)巖石參數(shù)符合Weibull分布(式(6))，將最大拉應(yīng)力和Mohr-Coulomb準(zhǔn)則作為判斷煤體RVE破壞的準(zhǔn)則(式(7))，建立了煤體的彈性損傷本構(gòu)。沿用這種研究思路，一些學(xué)者建立了不同的氣固耦合模型用以探究煤層中損傷-應(yīng)力-滲流的問題，間接搭建了損傷變量和滲透率之間的關(guān)系[44-52]。除了根據(jù)滲透率與孔隙度立方關(guān)系(式(4))推導(dǎo)滲透率表達(dá)式的常見方法[44-45]之外，YANG等[48]利用分段函數(shù)依次建立了煤體單元在壓縮和拉伸條件下的滲透率函數(shù)式;劉力源等[51-52]將煤體視為雙重孔隙介質(zhì)，結(jié)合煤體吸附CO2出現(xiàn)的損傷演化過程分別得到煤基質(zhì)和裂隙的滲透率計(jì)算式。

(6)

式中，C1，C2為兩個(gè)損傷破裂閾值函數(shù);u為單元參數(shù);u0為單元參數(shù)的平均值;m為函數(shù)的形狀參數(shù);σ1,σ3分別為單元的最大主應(yīng)力和最小主應(yīng)力;ft0為單軸抗拉強(qiáng)度;fc0為單軸抗壓強(qiáng)度;φ′為內(nèi)摩擦角。

2.2.2 滲透率表達(dá)式中含損傷變量的煤體滲透率模型

通過理論推導(dǎo)構(gòu)建損傷變量并將其嵌入滲透率表達(dá)式的研究主要限于連續(xù)損傷力學(xué)的框架內(nèi)進(jìn)行，大多數(shù)學(xué)者基于應(yīng)變等效假設(shè)展開研究，而損傷變量的構(gòu)建方法多樣。XUE等[53]通過孔隙度和有效應(yīng)力定義了損傷變量;ZHANG等[54]采用煤體的體積應(yīng)變、吸附應(yīng)變和初始孔隙率構(gòu)建了損傷變量;LI等[55]以軸向應(yīng)變構(gòu)建了損傷變量。另外，也有學(xué)者采用回歸分析的方法，通過擬合系數(shù)來確定滲透率和不同損傷程度下應(yīng)力狀態(tài)之間的關(guān)系。CHEN等[56]為了描述煤體從彈性狀態(tài)到峰后破壞的滲透率變化全過程，提出了修正系數(shù)L，通過回歸函數(shù)得到有效偏應(yīng)力與修正系數(shù)L的“S”形關(guān)系。ZHANG等[57]基于達(dá)西定律和裂隙的本構(gòu)方程，通過擬牛頓算法的回歸分析確定了以應(yīng)力表示的滲透率函數(shù)。為了搭建滲透率和損傷變量之間的“橋梁”，另外一些學(xué)者立足于試驗(yàn)結(jié)果反映出的變化規(guī)律，在前人建立的滲透率模型基礎(chǔ)之上進(jìn)行改造，直接添加反映損傷影響的作用項(xiàng)或修正系數(shù)對(duì)滲透率變化過程進(jìn)行推演，也取得了一些研究成果。ZHU等[58]在滲透率和孔隙度的立方關(guān)系(式(4))上進(jìn)行改造，添加了指數(shù)形式的損傷作用項(xiàng);ZHENG等[59]以CUI和BUSTIN[19]模型為基礎(chǔ)，在該模型中添加了損傷作用項(xiàng)，假設(shè)了滲透率和損傷變量的指數(shù)關(guān)系;XUE等[60]通過室內(nèi)試驗(yàn)的滲透率和聲發(fā)射數(shù)據(jù)反分析提出在峰后階段引入一個(gè)作用項(xiàng)來反映損傷影響，并假設(shè)峰后階段滲透率和損傷變量呈線性關(guān)系;XIE等[61]假設(shè)損傷作用下滲透率的變化符合式(8)的關(guān)系,即

(8)

式中，ζ為與材料有關(guān)的系數(shù)。

可以看出，考慮損傷破裂特征的滲透率模型研究在近十年才逐漸展開，該領(lǐng)域的研究尚處于起始階段，模型中損傷效應(yīng)的研究深度尚淺，損傷變量的構(gòu)建方式也較為分散。在模型適用性上，損傷變量存在本構(gòu)方程中的實(shí)現(xiàn)方式可以適用較多邊界條件情況，具有更廣的適用范圍，但是該方式下求解煤體滲透率需要依次計(jì)算應(yīng)力場(chǎng)和滲流場(chǎng)控制方程，計(jì)算工作量較大，所以該方式更適合采用數(shù)值計(jì)算。損傷變量存在滲透率表達(dá)式的實(shí)現(xiàn)方式使得研究滲透率與各因素之間的數(shù)量關(guān)系更加便捷，但其大多基于已有模型成果進(jìn)行改進(jìn)，受原模型構(gòu)建過程中的假設(shè)條件限制，諸如單軸應(yīng)變條件、上覆載荷不變等。損傷變量所處不同位置的兩種處理方式都可以達(dá)到在滲透率計(jì)算過程中考慮損傷破裂影響的目的。究其本質(zhì)而言，損傷變量存在于本構(gòu)方程的情況可視為損傷變量存在滲透率表達(dá)式情況的隱式函數(shù)形式。

2.3 煤體吸附解吸特征的表征方法

不同于其他氣體在巖石中的滲透行為，瓦斯在煤層中運(yùn)移時(shí)存在吸附解吸的特殊現(xiàn)象。煤層采掘過程中，煤體損傷破裂伴隨著瓦斯?jié)B流通道的形成。與此同時(shí)，開采擾動(dòng)打破了煤層內(nèi)原有的瓦斯壓力平衡狀態(tài)造成瓦斯壓力持續(xù)下降，煤體內(nèi)吸附態(tài)瓦斯釋放引起解吸現(xiàn)象的發(fā)生。吸附解吸作用直接影響煤基質(zhì)的脹縮變形，在滲透率模型構(gòu)建的過程中，大多通過吸附應(yīng)變這一中間變量間接反映煤體所受有效應(yīng)力的大小或煤體內(nèi)部滲流通道幾何尺寸的變化[62]，因此關(guān)于吸附應(yīng)變的計(jì)算方法就成為了研究吸附解吸影響下滲透率變化的焦點(diǎn)。

在早期的煤體滲透率模型中，通常將吸附應(yīng)變與氣體壓力的增量視為線性關(guān)系進(jìn)行處理。GRAY[15]通過提出等效吸附壓力的概念建立了吸附應(yīng)變與氣體壓力之間的線性關(guān)系。LEVINE[63]利用基質(zhì)收縮系數(shù)確定了基質(zhì)吸附應(yīng)變和氣體壓力增量之間的線性函數(shù)。國內(nèi)外學(xué)者普遍將煤對(duì)瓦斯的吸附認(rèn)為是物理吸附過程且為單分子層吸附?；诖?，Langmuir等溫吸附模型[64]可以較好的用于煤對(duì)瓦斯氣體吸附量的計(jì)算(式(9))。HARPALANI和CHEN[65-66]認(rèn)為在氣體吸附過程中吸附應(yīng)變和氣體吸附量成正比，采用Langmuir理論公式計(jì)算瓦斯氣體吸附量，建立了煤體吸附應(yīng)變和瓦斯氣體吸附量的線性關(guān)系(式(10))。鑒于Langmuir等溫吸附模型可以較好的描述氣體壓力改變下煤對(duì)瓦斯的吸附氣體量，LEVINE[63]對(duì)Langmuir模型作出進(jìn)一步改進(jìn)，采用一個(gè)和Langmuir吸附模型具有相同數(shù)學(xué)函數(shù)形式的方程(式(11))來描述煤體吸附應(yīng)變。

εs=αsVs(10)

式中，Vs為吸附的瓦斯氣體含量;VL為Langmuir體積常數(shù);αs為煤的吸附膨脹系數(shù);εL為Langmuir吸附體積應(yīng)變常數(shù);pL為Langmuir吸附壓力常數(shù)。

煤體對(duì)瓦斯吸附解吸的內(nèi)在機(jī)理可通過熱力學(xué)理論進(jìn)行探究，所以除了采用上述Langmuir形式的方程對(duì)吸附應(yīng)變進(jìn)行表征之外，一些學(xué)者通過吸附熱力學(xué)方法建立了吸附應(yīng)變的計(jì)算式[67-72]，詳見表4。

表4吸附熱力學(xué)角度建立的吸附應(yīng)變表達(dá)式
Table4Expressionsofadsorptionstrainbasedonadsorptionthermodynamics

文獻(xiàn)來源吸附應(yīng)變表達(dá)式參數(shù)釋義PAN和CONNELL[67]εs=RTAln(1+Bp)ρmEmf(x,νm)-pEm(1-2νm)R為氣體常數(shù);T為溫度;A為Langmuir吸附常數(shù);ρm為煤基質(zhì)密度;Em為基質(zhì)的彈性模量;νm為基質(zhì)的泊松比LIU和HARPALANI[68]εl=3AρmRTEAVg∫p0B1+Bpdpεl為吸附/解吸引起的線性應(yīng)變;EA為吸附解吸的膨脹模量;Vg為氣體摩爾體積LIU等[69]εs=nsatRTVcEc[ln(1+Bpt)-ln(1+Bpr)]nsat為飽和吸附摩爾數(shù);Vc為煤的體積;Ec為煤的彈性模量;pr為起始?xì)怏w壓力;pt為終止氣體壓力VERMOREL和PIJAUDIER-CABOT[70]εs-εs0=∫pp0dpK1-KKs (1-χ)-1-1 ?χ為約束度HOL等[71]εsij=Vm0ρ(Aij+δij)3Cσεsij為各個(gè)方向的吸附應(yīng)變;Vm0為吸附相的偏摩爾體積;ρ為煤的密度;Aij為各向異性膨脹張量;δij為克羅內(nèi)克符號(hào);Cσ為吸附分子的摩爾數(shù)

采用以氣體壓力增量為自變量的線性函數(shù)對(duì)煤吸附應(yīng)變進(jìn)行描述是最早的研究方法，其優(yōu)點(diǎn)是瓦斯氣體壓力的測(cè)定較為簡(jiǎn)單，但是線性函數(shù)卻有夸大瓦斯壓力對(duì)吸附應(yīng)變影響的缺點(diǎn)。通過對(duì)比試驗(yàn)結(jié)果發(fā)現(xiàn)Langmuir等溫吸附理論可以很好的描述煤對(duì)瓦斯的吸附解吸過程，上述Langmuir形式的吸附應(yīng)變表達(dá)式具有一定的理論基礎(chǔ)且方程形式較為簡(jiǎn)潔，表達(dá)式中的VL,εL,pL等參數(shù)也較易測(cè)試，所以式(10)[73-76]和式(11)[77-80]被后來學(xué)者大量效仿使用。通過熱力學(xué)方法建立的吸附應(yīng)變大多以吸附熱力學(xué)為基礎(chǔ)進(jìn)行推導(dǎo)，由此得出的吸附應(yīng)變表達(dá)式具有更為嚴(yán)密的理論基礎(chǔ)，然而此類方法的模型參數(shù)較多且部分參數(shù)難以準(zhǔn)確測(cè)定，所以此類吸附應(yīng)變表達(dá)式的實(shí)用性有所限制。

3 采動(dòng)應(yīng)力下煤體滲透率模型研究進(jìn)展

現(xiàn)有的煤體滲透率模型中同時(shí)考慮各向異性、損傷破裂、吸附解吸等特征的滲透率模型還較少。近些年，部分學(xué)者已經(jīng)開始關(guān)注這一問題，并建立了相應(yīng)的采動(dòng)應(yīng)力下煤體滲透率模型。筆者將這些模型成果歸納為有效應(yīng)力型、幾何參數(shù)型和系數(shù)擬合型的3種進(jìn)行介紹。

3.1 有效應(yīng)力型的煤體滲透率模型

有效應(yīng)力型的煤體滲透率模型主要以有效應(yīng)力、氣體壓力、彈性模量和泊松比等力學(xué)參數(shù)為基本變量對(duì)滲透率的變化進(jìn)行表征。其中，ZHANG等[81]基于立方體模型結(jié)構(gòu)(圖1(c))，運(yùn)用損傷和斷裂力學(xué)理論提出了一個(gè)可考慮采動(dòng)應(yīng)力引起微裂紋擴(kuò)展、吸附解吸造成煤體變形的有效應(yīng)力型滲透率模型(式(12))，并根據(jù)謝和平等[82]提出的3種典型開采方式下采動(dòng)應(yīng)力路徑(圖2)，對(duì)采動(dòng)煤體滲透率變化規(guī)律進(jìn)行了擬合分析。

圖2(a)～(c)為不同開采方式下煤體滲透率的試驗(yàn)和擬合結(jié)果?？梢钥闯鲈?種開采方式下，隨著采動(dòng)卸壓過程的進(jìn)行，煤體內(nèi)微裂隙不斷萌生和擴(kuò)展，損傷積累不斷增大，滲透率的實(shí)測(cè)和擬合結(jié)果持續(xù)升高。無煤柱和放頂煤開采方式的滲透率在加卸載途中出現(xiàn)陡然上升的情況，但擬合結(jié)果卻沒有出現(xiàn)陡然升高的特征。不過整體來看，式(12)基本可以反映出煤體滲透率隨著采動(dòng)卸壓過程的進(jìn)行而持續(xù)升高的變化規(guī)律。

3.2 幾何參數(shù)型的煤體滲透率模型

幾何參數(shù)型的煤體滲透率模型主要以孔隙率和煤體應(yīng)變?yōu)榛緟?shù)來構(gòu)建滲透率模型。ZHENG等[83]將LIU等[28]提出的煤體各向異性滲透率模型和ZHU等[58]定義的損傷破裂特征修正項(xiàng)進(jìn)行聯(lián)結(jié)。通過細(xì)觀損傷力學(xué)的方法構(gòu)建損傷變量以對(duì)滲透率進(jìn)行修正，合成出可表征采動(dòng)應(yīng)力下煤體損傷的幾何參數(shù)型滲透率表達(dá)式:

圖2 不同開采方式下應(yīng)力路徑及煤體滲透率演化[81-82]Fig.2 Stress paths and permeability of coal under different mining layouts[81-82]

式中，εj為j方向的應(yīng)變;εt0,εc0分別為最大拉伸主應(yīng)變和最大壓縮主應(yīng)變;ε1,ε3分別為最大主應(yīng)變和最小主應(yīng)變。

其后，利用式(13)分別對(duì)考慮和不考慮損傷效應(yīng)的保護(hù)層底板下煤巖體滲透率進(jìn)行了模擬分析。由圖3可以看出距離保護(hù)層底板下方80 m處的巖體滲透率為初始狀態(tài);在與底板距離縮進(jìn)過程中煤巖體滲透率逐漸升高;在距離底板12 m之外，考慮損傷與否的滲透率結(jié)果相同;而在12 m范圍之內(nèi)，考慮損傷效應(yīng)的煤巖體滲透率比率出現(xiàn)急劇上升，從初始滲透率比率的180倍陡增至680倍左右，而不考慮損傷的煤巖體滲透率比率僅升高至200倍左右。

圖3 保護(hù)層底板下煤體滲透率分布[83]Fig.3 Permeability distribution of the coal under the first-mined seam floor[83]

3.3 系數(shù)擬合型的煤體滲透率模型

系數(shù)擬合型的煤體滲透率模型主要通過試驗(yàn)數(shù)據(jù)擬合的手段來確定一系列待定系數(shù)的變化與不同區(qū)段滲透率的關(guān)系。XUE等[84]根據(jù)煤體三軸壓縮試驗(yàn)過程中應(yīng)力-應(yīng)變-滲透率關(guān)系的變化特征(圖4)，將全應(yīng)力-應(yīng)變過程劃分為3階段:彈性階段、塑性軟化階段、殘余狀態(tài)階段。同時(shí)認(rèn)為彈性階段沒有微裂紋產(chǎn)生，塑性軟化階段萌生大量裂紋，進(jìn)而結(jié)合巷道開挖后圍巖損傷的3區(qū)特征，將圍巖的漸進(jìn)損傷區(qū)與全應(yīng)力-應(yīng)變過程3階段進(jìn)行了一一對(duì)應(yīng):彈性階段對(duì)應(yīng)無損區(qū)、塑形軟化階段對(duì)應(yīng)開挖損傷區(qū)、殘余狀態(tài)階段對(duì)應(yīng)強(qiáng)開挖損傷區(qū)。由此，得到了巷道開挖后，掘進(jìn)工作面周圍煤體的系數(shù)擬合型滲透率表達(dá)式

(14)

式中，αki,βki,γki(i=1,2,3)分別為不同應(yīng)力和應(yīng)變條件下描述滲透率變化的系數(shù);εip為煤進(jìn)入塑性狀態(tài)的應(yīng)變;εir為煤進(jìn)入殘余狀態(tài)的應(yīng)變。

圖4 煤的應(yīng)力-應(yīng)變-滲透率關(guān)系[84]Fig.4 Diagram of stress-strain-permeability for coal samples[84]

圖5為掘進(jìn)工作面推進(jìn)過程中周圍煤體內(nèi)滲透率分布演化結(jié)果，周圍煤體的滲透率分布同樣符合三區(qū)特征。距離巷道中心較遠(yuǎn)處，煤體滲透率處于初始狀態(tài);隨著與巷道中心距離的縮減，無損區(qū)內(nèi)的煤體滲透率由于初始裂隙的閉合而滲透率逐漸減小(圖5 Ⅰ)，距離至塑性區(qū)的邊界時(shí)，滲透率達(dá)到最小值;開挖損傷區(qū)內(nèi)由于產(chǎn)生大量裂紋而造成滲透率持續(xù)升高(圖5 Ⅱ);強(qiáng)開挖損傷區(qū)內(nèi)煤體破裂程度較為嚴(yán)重，由此出現(xiàn)大量貫通裂隙引起滲透率急劇升高(圖5 Ⅲ)。對(duì)比不同開挖時(shí)間下的滲透率分布結(jié)果可知，隨著時(shí)間的推移，煤體內(nèi)滲透率的變化幅度呈上升趨勢(shì)。

圖5 掘進(jìn)工作面周圍煤體滲透率分布[84]Fig.5 Permeability distribution of the coal around the roadway[84]

開采擾動(dòng)下工作面前方煤體或底板下煤巖體的滲透率均呈現(xiàn)分區(qū)段變化的特征，上述3個(gè)滲透率模型分別從斷裂力學(xué)、細(xì)觀損傷力學(xué)和連續(xù)損傷力學(xué)的角度考慮了采動(dòng)過程中煤體損傷破裂造成的滲透率分區(qū)特征。在此之中，有效應(yīng)力型的模型運(yùn)用斷裂力學(xué)理論，確定了微裂紋的開裂判據(jù)，通過損傷理論構(gòu)建了微裂紋和損傷變量的關(guān)系，建立的模型較好的體現(xiàn)了采動(dòng)應(yīng)力和損傷破裂特征，模型的理論推導(dǎo)過程較為嚴(yán)密，具有較好的基礎(chǔ)理論支撐。幾何參數(shù)型的模型加入了煤體參數(shù)的微觀統(tǒng)計(jì)分布，更加貼近煤體內(nèi)部孔裂隙結(jié)構(gòu)導(dǎo)致的非均勻性特征，同時(shí)模型采用拉伸和剪切破壞準(zhǔn)則作為單元失效判據(jù)更為全面，模型從細(xì)觀損傷本構(gòu)出發(fā)考慮損傷破裂的因素對(duì)煤體變形和滲透的影響，適用范圍較廣，但計(jì)算中將煤體劃分為細(xì)觀單元后，求解計(jì)算的工作量較大，所以較為適合進(jìn)行數(shù)值計(jì)算。系數(shù)擬合型的模型從形式上來看較為簡(jiǎn)潔，可顯著區(qū)分工作面前方煤體的損傷分布特征，可以較好的反映各損傷區(qū)煤體滲透率的變化規(guī)律，此外該模型也易于進(jìn)行理論計(jì)算，尤其是便于解析計(jì)算。

4 結(jié)論與展望

(1)相比天然氣和煤層氣的開采過程，采動(dòng)應(yīng)力環(huán)境下的煤層在受力方式和結(jié)構(gòu)形式上更為復(fù)雜。工作面推進(jìn)過程中采動(dòng)應(yīng)力的3個(gè)應(yīng)力主方向出現(xiàn)不同的變化趨勢(shì)，采動(dòng)應(yīng)力可視為外部邊界條件，而煤層由于沉積作用出現(xiàn)的層理和割理結(jié)構(gòu)屬于內(nèi)部結(jié)構(gòu)條件。煤炭開采過程中，其內(nèi)部結(jié)構(gòu)和外部邊界條件的主方向并非一致，所以需要考慮外部邊界條件和內(nèi)部滲透率主方向的夾角問題，統(tǒng)一兩者之間的關(guān)系。與此同時(shí)，割理結(jié)構(gòu)的各向異性也造成煤體物理力學(xué)參數(shù)存在異性特征。因此，需要結(jié)合地質(zhì)條件和各向異性的物理力學(xué)參數(shù)建立一個(gè)考慮采動(dòng)應(yīng)力、煤層產(chǎn)狀和滲透率方向的模型，以此確定煤層瓦斯?jié)B透的最優(yōu)化方向，指導(dǎo)煤層增透和瓦斯抽采工作。

(2)采動(dòng)卸壓是一個(gè)打破原始地應(yīng)力后應(yīng)力持續(xù)重分布引起煤層損傷破裂的增透過程?，F(xiàn)有的煤體滲透率模型已經(jīng)開始逐漸關(guān)注到這個(gè)問題，然而在已有成果中完全基于損傷理論推導(dǎo)出的損傷滲透率模型還較少，部分經(jīng)驗(yàn)?zāi)Ｐ偷膶?shí)用性較好，但推導(dǎo)過程存在較多前提假設(shè)。因此，損傷滲透率模型的理論基礎(chǔ)仍需進(jìn)一步夯實(shí)。除此之外，結(jié)合斷裂力學(xué)理論對(duì)煤體內(nèi)裂隙演化過程進(jìn)行定量描述，探究裂隙連通度對(duì)煤體滲透性的影響，確定采場(chǎng)尺度下不同裂隙網(wǎng)絡(luò)空間分布與煤層滲透性的關(guān)系[85-88]，建立一個(gè)多理論交叉、宏細(xì)觀結(jié)合的滲透率模型是十分必要的。

(3)煤體吸附應(yīng)變的研究中，大量的實(shí)驗(yàn)和理論成果都是基于等溫條件獲得，而采面前方煤體內(nèi)的氣體壓力和溫度由表及內(nèi)卻是不斷升高的。此外，在瓦斯吸附和解吸過程中，煤體應(yīng)變嚴(yán)格來講并不具有完全可逆的時(shí)空演化特征[89-91]。因此，加強(qiáng)非等溫條件下瓦斯吸附特性的研究，掌握瓦斯解吸和吸附過程中煤體應(yīng)變的差異性對(duì)采動(dòng)影響下煤層滲透率演化過程的精確描述大有裨益。

總體來說，多因素影響下的采動(dòng)煤體滲透率模型的理論體系并未完全建立，能夠較為系統(tǒng)和完善體現(xiàn)上述特征的采動(dòng)煤體滲透率模型也還鮮有出現(xiàn)，但煤體滲透率模型在各向異性特征、損傷破裂特征和吸附解吸特征的方面均已取得大量研究成果，以這些成果為基礎(chǔ)的采動(dòng)煤體滲透率模型將會(huì)迅速發(fā)展。同時(shí)，以上各方面的研究深度及廣度仍需進(jìn)一步加強(qiáng)。除此之外，采動(dòng)應(yīng)力下煤體滲透率模型的研究還有一些問題需要面對(duì)。如今，煤炭開采深度的迅速增大給煤體滲透率的研究帶來巨大挑戰(zhàn)。一方面，深部煤體的賦存環(huán)境具有“三高”特征，即高溫、高地應(yīng)力和高巖溶水壓[92-94]。深部開采是多場(chǎng)多相環(huán)境共同影響的工程活動(dòng)，采掘和增透過程中“三高”因素的變化相互耦合，任一因素的變化都將嚴(yán)重影響煤層的透氣性效果。因此，深部采動(dòng)煤體滲透率模型需要解決多場(chǎng)多相耦合的問題。另一方面，深部煤體開采過程中具有顯著的強(qiáng)擾動(dòng)和強(qiáng)時(shí)效特征[95-96]。強(qiáng)擾動(dòng)是指采掘活動(dòng)對(duì)深部高應(yīng)力區(qū)巖體形成的強(qiáng)烈卸荷效應(yīng);強(qiáng)時(shí)效是指深部環(huán)境下采動(dòng)巖體具有與采掘活動(dòng)無關(guān)的、明顯的流變效應(yīng)。顯而易見，強(qiáng)擾動(dòng)和強(qiáng)時(shí)效特征將會(huì)強(qiáng)烈影響深部煤體的滲透特性。所以，為了全面反映深部煤體滲透率的控制因素，深部采動(dòng)煤體滲透率模型還需要考慮強(qiáng)擾動(dòng)和強(qiáng)時(shí)效特征的影響。