考慮時效影響的深部煤層瓦斯運移特性

張 雷，周宏偉，王向宇，榮騰龍，陳超凡

(1.中國礦業(yè)大學(xué)(北京) 力學(xué)與建筑工程學(xué)院，北京 100083; 2.中國礦業(yè)大學(xué)(北京) 能源與礦業(yè)學(xué)院，北京 100083; 3.中國礦業(yè)大學(xué)(北京) 煤炭資源與安全開采國家重點實驗室，北京 100083; 4.德國亥姆霍茲環(huán)境研究中心(UFZ)，德國 萊比錫 04318)

隨著煤炭資源開采向深部發(fā)展，資源地質(zhì)賦存條件更加的復(fù)雜，煤巖與瓦斯的賦存環(huán)境處于高地應(yīng)力、高地溫、高滲透壓[1]的狀態(tài)，導(dǎo)致煤與瓦斯突出礦井和突出煤層數(shù)量不斷增加，瓦斯突出災(zāi)害的發(fā)生次數(shù)和突出強度相比于淺部煤層呈現(xiàn)上升趨勢。瓦斯同時又是一種高效的清潔能源[2]，若能得到合理開發(fā)利用，可以有效地改善能源結(jié)構(gòu)。研究深部煤層瓦斯抽采運移特征對煤氣資源的安全高效開采具有重要意義[3-4]。

前人基于不同假設(shè)條件，針對煤巖體的孔隙率、滲透率模型進行了深入研究。SEIDLE等[5]提出了一個動態(tài)的孔隙率演化模型，該模型考慮了煤體的膨脹收縮效應(yīng)。PALMER等[6]考慮了孔隙壓力和基質(zhì)膨脹收縮對割理裂隙的孔隙率演化影響，提出了一個基于單軸應(yīng)變和恒定垂向應(yīng)力假設(shè)的滲透率模型。SHI等[7]提出的滲透率模型考慮了吸附解吸引起的煤體變形，并提出了考慮水平有效應(yīng)力的滲透率模型。CUI等[8]使用線彈性多孔彈性介質(zhì)理論，考慮吸附解吸效應(yīng)推導(dǎo)出了與平均有效應(yīng)力有關(guān)的滲透率模型。滲透率隨著有效應(yīng)力的變化歸根結(jié)底是孔隙率隨著有效應(yīng)力的變化[9]。

上述研究均未考慮時間效應(yīng)對煤巖體滲透性的影響，而時間效應(yīng)對于煤巖體的強度和變形有較大影響，同時也會造成滲透率的改變。在煤巖體時效特性的研究方面，周輝等[10]通過分析巖石破裂過程中強度的演化規(guī)律和機制，建立了巖石強度時效性演化模型，定量描述了恒定應(yīng)力條件下巖石強度的弱化過程，從理論和機制上解釋了巖石強度的時間效應(yīng)問題。周宏偉等[11]從分?jǐn)?shù)階導(dǎo)數(shù)出發(fā),在常黏性系數(shù)Abel黏壺基礎(chǔ)上提出了一種新的變黏性系數(shù)的Abel黏壺元件。利用兩種分?jǐn)?shù)階Abel黏壺代替經(jīng)典西原模型中Newton黏壺的方法,建立了基于分?jǐn)?shù)階導(dǎo)數(shù)的鹽巖流變本構(gòu)模型,并給出了該模型的解析解，較好地反映了鹽巖流變的三階段尤其是加速流變階段。尹光志等[12]對煤層頂?shù)装鍘r層試件進行卸圍壓蠕變及周期加載蠕變實驗，并采用Burgers 體模型對采動影響下的巖體進行蠕變理論分析，建立了采動下煤巖蠕變損傷模型，較好地反映巖體蠕變曲線特征。周長冰等[13]分析了高溫三軸應(yīng)力下氣煤蠕變特征，通過對不同溫度下氣煤的滲透率和孔隙率的分析,可以初步判斷氣煤蠕變特征發(fā)生變化的臨界溫度。許江等[14]進行了考慮時間效應(yīng)的蠕變對含瓦斯煤滲透率影響的試驗分析，發(fā)現(xiàn)在相同溫度和有效應(yīng)力條件下，經(jīng)歷蠕變后煤的滲透率會降低。郝富昌等[15]建立了鉆孔周圍煤體黏彈塑性模型和蠕變-滲流耦合作用下瓦斯運移模型，認為不同埋深鉆孔均會隨時間產(chǎn)生縮孔現(xiàn)象，進而確定了不同埋藏深度鉆孔的有效抽采半徑。

為了考慮煤巖體的時效性特征對瓦斯在煤層中運移的影響以及準(zhǔn)確描述煤巖的孔隙率和滲透率演化過程，筆者以煤層割理裂隙的滲透性受三向應(yīng)力影響的角度出發(fā)，借鑒CUI等[8]引入平均有效應(yīng)力變化量來描述煤層滲透情況的思想，利用CUI等[8]建立的孔隙率、滲透率與平均有效應(yīng)力變化量之間的指數(shù)關(guān)系，并結(jié)合蠕變本構(gòu)建立了考慮時效性的孔隙率和滲透率模型，對比和分析了在模擬深部瓦斯運移過程中的作用。

1 考慮時效性的煤體本構(gòu)方程

隨著埋深的增加和瓦斯抽采時間的延長，含瓦斯煤的時效性特征越來越明顯。 NISHIHARA[16]提出了準(zhǔn)靜態(tài)流變模型——西原模型，用來描述時間效應(yīng)影響下巖石的應(yīng)力應(yīng)變關(guān)系，如圖1所示。該模型中包括胡克體，黏彈性體和黏塑性體。

圖1 經(jīng)典西原模型Fig.1 Classical Nishihara model

根據(jù)不同元件的串聯(lián)方式，整個模型的總應(yīng)變可表示為

ε=εe+εve+εvp

(1)

式中，εe,εve和εvp分別為胡克體、黏彈性體和黏塑性體的應(yīng)變。

本文只考慮煤層所受應(yīng)力小于其屈服應(yīng)力σs的情況，則西原模型的蠕變本構(gòu)關(guān)系可表示為

(2)

式中，Ee為煤體的彈性模量,GPa;Eve為煤體的黏彈性模量,GPa;ηve為煤體的黏塑性模量,GPa。

不考慮時間效應(yīng)的應(yīng)力應(yīng)變本構(gòu)方程由JAEGER[17]提出:

(3)

則考慮煤體基質(zhì)蠕變效應(yīng)、瓦斯吸附解吸的有效應(yīng)力應(yīng)變本構(gòu)關(guān)系可由式(3)改寫成

(4)

式(4)中系數(shù)A,B，C表達式為

(5)

(6)

式中，εL為Langmuir體應(yīng)變常數(shù);PL為Langmuir壓力常數(shù),MPa;p為孔隙壓力,MPa。

對于煤儲層，可認為其受三向應(yīng)力狀態(tài)，受力示意如圖2所示。同時假設(shè)在鉆孔周圍很大范圍內(nèi)的煤層邊界在x,y方向應(yīng)變較小，可忽略不計，在z方向存在應(yīng)變[15]。

圖2 煤層受力示意Fig.2 Distribution of coal seam stress

在上述應(yīng)變假設(shè)下，有效應(yīng)力應(yīng)變本構(gòu)關(guān)系可由式(4)簡化為

(7)

考慮上覆巖層總載荷τ不變，即Δτ=0，利用有效應(yīng)力原理導(dǎo)出z方向的有效應(yīng)力變化量，具體形式為

Δσz=Δ(τ-αp)=Δτ-Δαp=-Δαp

(8)

式中，α為Biot系數(shù)。

煤是一種典型的沉積巖，具有橫觀各向異性的特征，可以認為煤層在x,y方向的物理力學(xué)性質(zhì)相同，而垂直于xy平面的z方向物理力學(xué)性質(zhì)不同，則在x,y方向有效應(yīng)力的變化量可通過式(7)變換成

(9)

則考慮時間效應(yīng)的平均有效應(yīng)力改變量可由式(8)和式(9)表示為

(10)

式(10)中第1項、第2項表示彈性和黏彈性割理裂隙壓縮性，而第3項代表彈性和黏彈性基質(zhì)收縮。其中

(11)

2 考慮時效性的煤體孔隙率和滲透率模型

為了計算時間效應(yīng)和瓦斯吸附解吸對煤體滲透率演化的影響。筆者借助CUI-BUSTIN模型(C-B模型)的思想[8,18]，在描述孔隙率與滲透率動態(tài)變化過程中引入平均有效應(yīng)力變化量。CUI等[18]指出孔隙率與平均有效應(yīng)力關(guān)系表示為

(12)

式中，φ0為初始孔隙率;σe為平均有效應(yīng)力改變量,MPa;Kp為孔隙體積模量,GPa。

此外，針對孔隙體積模量Kp與煤體基質(zhì)的彈性模量K之間關(guān)系可表示為

(13)

其中，E為煤體的彈性模量,GPa;ν為泊松比。將式(13)代入式(12)，則孔隙率演化模型可表示為

(14)

此外，滲透率與孔隙率之間可認為滿足如下冪函數(shù)關(guān)系式:

(15)

則滲透率演化可表示為

(16)

其中，cf=1/Kp為割理裂隙壓縮系數(shù)。根據(jù)所述孔隙率、滲透率與平均有效應(yīng)力變化量之間的關(guān)系，將式(10)和(11)代入到式(14)和(16)，導(dǎo)出基于時間效應(yīng)考慮的孔隙率與滲透率模型:

(17)

(18)

3 瓦斯運移方程

鑒于瓦斯在煤體中的流動過程較為復(fù)雜，控制方程基于如下假設(shè)構(gòu)建:① 主要考慮割理裂隙滲透作用;② 瓦斯在裂隙割理中的運動規(guī)律滿足達西定律，忽略氣體的Klinkenberg效應(yīng);③ 利用Langmuir等溫吸附方程描述氣體解吸過程。

氣體在流動過程中的連續(xù)性方程可表示為

(19)

式中，H為自由項和解吸項的氣體含量,kg/m3;v為氣體的達西速度,m/s;ρ為氣體密度,kg/m3;Qs為氣體源匯項,kg/(m3·s)。

根據(jù)CUI和BUSTIN[8]在文獻中表述，H可表示為

(20)

其中,ρa為標(biāo)準(zhǔn)大氣壓下的氣體密度,kg/m3;ρc為煤的密度,kg/m3;VL是Langmuir體積常數(shù),m3/kg;PL是Langmuir壓力常數(shù),MPa;p為氣體壓力,MPa。其中，氣體密度可認為符合理想氣體狀態(tài)方程，形式為

(21)

式中，M為瓦斯氣體的摩爾質(zhì)量,g/mol;R為理想氣體常數(shù),J/(mol·K);T為溫度,K。

一般氣體滲流過程忽略重力的影響，因此，在不考慮重力影響的條件下，根據(jù)達西定律，氣體流速為

(22)

式中，μ為瓦斯氣體動力黏度系數(shù),Pa·s。

根據(jù)式(19)～(22)，新的氣體連續(xù)性方程可表示為

(23)

將考慮時間效應(yīng)的滲透率和孔隙率演化模型代入式(23)，通過COMSOL Multiphysics軟件進行數(shù)值求解，從而得到同時考慮時間效應(yīng)和解吸效應(yīng)的瓦斯運移規(guī)律。

4 瓦斯運移特性模擬分析

4.1 計算模型及參數(shù)

計算模型如圖3所示，二維均勻介質(zhì)長寬均為10 m，鉆孔半徑為0.15 m，模型中初始瓦斯氣壓為6 MPa，在xy平面內(nèi)，假設(shè)初始孔隙率為0.005、初始滲透率為1×10-15m2，并且在各個方向均相等。鉆孔邊界處壓力為0.1 MPa，煤層內(nèi)與鉆孔外存在氣體壓力差，導(dǎo)致氣體的釋放及瓦斯解吸。根據(jù)文獻[19]的假設(shè)，煤體在水平面內(nèi)(x，y方向)不存在應(yīng)變，且在z方向上整體荷載恒定。由于模型尺寸較大，邊界不受抽采影響，因此，所有邊界均為不可滲邊界，鉆孔邊界除外。模型中具有深部特征的參數(shù)來自文獻[19-21]，見表1。

圖3 數(shù)值模擬模型Fig.3 Numerical simulation model

表1 數(shù)值模擬參數(shù)Table 1 Numerical simulation parameters

(24)

通過式(24)的簡化形式，認為長時間抽采過程中x,z方向的彈性系數(shù)可表示為黏彈性模量和彈性模量的線性疊加，并將式(24)代入式(17),(18)進行計算。

4.2 滲透率模型有效性驗證

選取圖3中距鉆孔中心0.3 m的位置為監(jiān)測點。為了驗證本文提出的考慮時間效應(yīng)模型的有效性，將本文模型、DANESH-CHEN模型[19](簡稱D-C模型)、現(xiàn)場實測數(shù)據(jù)[21]3者滲透率比率(k/k0)隨瓦斯壓力變化的曲線進行對比，結(jié)果如圖4所示。

圖4 滲透率比與壓力關(guān)系Fig.4 Relationship between permeability ratio and gas pressure

由圖4可知，抽采初期瓦斯壓力從6 MPa開始下降，本文考慮時間效應(yīng)模型、D-C模型與現(xiàn)場數(shù)據(jù)擬合較好，說明在該壓力區(qū)間本文模型有較好的適用性;抽采中后期即瓦斯壓力小于4 MPa階段，隨著瓦斯壓力逐漸下降，本文模型、D-C模型分別與現(xiàn)場數(shù)據(jù)存在不同程度的誤差，其中D-C模型與現(xiàn)場數(shù)據(jù)的差別逐漸增大，而本文模型卻與現(xiàn)場數(shù)據(jù)逐漸接近，尤其在壓力為1 MPa時已經(jīng)近似等于現(xiàn)場數(shù)據(jù)的滲透率比率。針對模擬長期瓦斯抽采，通過不同瓦斯壓力階段的對比，本文提出考慮時效性的模型優(yōu)于D-C模型，能較好地模擬孔隙率、滲透率在抽采過程中動態(tài)變化情況。

本文模型在長期抽采的有效性已經(jīng)驗證，現(xiàn)討論短期瞬態(tài)過程有效性?？紤]時間效應(yīng)的孔隙率、滲透率模型與其他瞬態(tài)模型存在的區(qū)別在于時間項t,現(xiàn)將時間t=0代入式(10)，從而使式(17)和式(18)退化成未考慮時間效應(yīng)的孔隙率與滲透率的演化模型，可表示為

(25)

(26)

則瞬態(tài)有效性可采用未考慮時間效應(yīng)的滲透率模型式(26)驗證，并與未考慮時間效應(yīng)的SHI-DURUCAN模型[7](簡稱S-D模型)進行對比，具體趨勢如圖5所示。

圖5 滲透率與壓力關(guān)系Fig.5 Relationship between permeability and gas pressure

圖5呈現(xiàn)了兩種未考慮時間效應(yīng)的滲透率模型變化情況，隨著瓦斯壓力逐漸降低，滲透率呈現(xiàn)明顯的上升趨勢。抽采初期，兩種模型滲透率數(shù)值差別較小，抽采中后期，兩種模型滲透率數(shù)值差別明顯增大。通過與圖4中數(shù)據(jù)對比，發(fā)現(xiàn)兩種未考慮時間效應(yīng)模型的滲透率在數(shù)值上大于考慮時間效應(yīng)下的滲透率，對比結(jié)果與許江等[14]室內(nèi)蠕變滲流實驗結(jié)果相同，未經(jīng)歷蠕變的滲透率較大。同時，圖4,5呈現(xiàn)的曲線表明，吸附解吸作用對上述的每個滲透率模型有重要的影響，在初始儲層壓力6 MPa至鉆孔壁標(biāo)準(zhǔn)大氣壓整個壓力降范圍內(nèi)，使得滲透率隨著壓力降低呈指數(shù)型增長。通過對比考慮時效性模型、D-C模型在模擬瓦斯抽采時分別與現(xiàn)場數(shù)據(jù)的擬合程度，驗證了本文提出的模型在模擬短期、長期瓦斯抽采的有效性。

圖6 孔隙率、滲透率和壓力隨時間變化趨勢Fig.6 Evolution of porosity,permeability and pressure with time

4.3 考慮時效影響的瓦斯運移參數(shù)變化分析

圖6(a)為監(jiān)測點孔隙率計算結(jié)果，孔隙率隨著抽采時間增長呈現(xiàn)增長趨勢?？紤]時間效應(yīng)與未考慮時間效應(yīng)影響模型的孔隙率變化會受到解吸效應(yīng)的影響，氣體解吸導(dǎo)致基質(zhì)收縮和裂隙內(nèi)瓦斯壓力增大，裂隙寬度變大，孔隙率隨之增大。而考慮時間效應(yīng)的孔隙率明顯小于未考慮時間效應(yīng)模型的孔隙率，說明隨著時間增長，一方面是蠕變本構(gòu)中黏彈性元件的存在使得煤體更為致密，孔隙率變小，另一方面是抽采中應(yīng)力重分布導(dǎo)致煤層2次壓實，孔隙率變小。而圖6(b)中滲透率變化趨勢取決于孔隙率的變化，則通過圖6(a)可知，滲透率會隨著抽采時間增長呈現(xiàn)上升趨勢。隨著抽采時間增長，時間效應(yīng)越發(fā)明顯，考慮時間效應(yīng)模型與未考慮時間效應(yīng)模型的滲透率差值逐漸增大。

瓦斯抽采過程中壓力隨時間的變化情況如圖6(c)所示，瓦斯壓力隨時間增長呈下降趨勢，瓦斯儲層初始壓力為6 MPa，隨著抽采進行，瓦斯壓力下降趨勢明顯，且在考慮時間效應(yīng)的瓦斯壓力下降幅度略低于未考慮時效性模型，在90 000 s時下降到0.735 MPa。由于瓦斯解吸效應(yīng)的存在，瓦斯壓力受解吸氣體含量和游離態(tài)氣體含量控制，而瓦斯解吸量與游離態(tài)瓦斯氣體壓力密切相關(guān)。因此，時間效應(yīng)對鉆孔周圍的瓦斯氣體壓力分布影響較小。

考慮時效性的滲透率隨瓦斯壓力變化的趨勢如圖7所示，通過對比發(fā)現(xiàn)在滿足煤體所受應(yīng)力小于其屈服強度和瓦斯壓力降相同時，考慮時間效應(yīng)模型和未考慮時間效應(yīng)模型中滲透率隨瓦斯壓力變化差別明顯，在同樣氣體壓力下，前者小于后者的滲透率。究其原因，蠕變本構(gòu)條件下黏彈性元件的存在，使煤體更致密，且氣體解吸效應(yīng)影響被弱化，隨著氣體壓力的減小，時間效應(yīng)會更為明顯。

圖7 滲透率與壓力關(guān)系Fig.7 Relationship between permeability and gas pressure

4.4 考慮時效影響的瓦斯運移特征分析

考慮模擬抽采鉆孔周圍的瓦斯運移情況，為實際瓦斯抽采過程中預(yù)估鉆孔周圍的滲透率與壓力的變化情況提供數(shù)值參考。

由圖8,9變化趨勢可知，在同一抽采時刻，隨著距鉆孔中心距離的增大，滲透率呈現(xiàn)降低的趨勢，瓦斯壓力與滲透率的變化趨勢相反，距中心距離越大，壓力越高;而當(dāng)臨近鉆孔中心時,滲透率上升幅度和瓦斯壓力下降的幅度都比較大，這是由于越接近鉆孔中心,煤層的卸壓越明顯,煤層滲透率也相應(yīng)變大,瓦斯壓力下降也就越明顯。

圖8 不同時刻滲透率演化Fig.8 Evolution of permeability at different times

圖9 不同時刻瓦斯壓力演化Fig.9 Evolution of gas pressure at different times

圖10 不同時刻瓦斯壓力和滲透率隨演化趨勢Fig.10 Evolution of gas pressure and permeability at different times

在考慮時間效應(yīng)的情況下，隨著抽放時間的推移,煤層中的滲透率上升程度和瓦斯壓力下降程度逐漸減小。這類趨勢表明:在鉆孔抽采瓦斯的初始階段,受時間效應(yīng)影響較小，瓦斯抽放量較大;而一段時間后受瓦斯壓力降低和時間效應(yīng)的影響瓦斯抽放量將趨于穩(wěn)定。當(dāng)模擬瓦斯抽采時間為1,10,30 d時，滲透率隨著抽采時間的增長逐漸增大，各時刻模型內(nèi)最小滲透率分別為0.4×10-14,1.6×10-14,1.81×10-14m2;而瓦斯壓力逐漸減小，各時刻模型內(nèi)最大瓦斯壓力分別為1.95,0.21,0.107 MPa。

不同時刻滲透率和瓦斯壓力的演化趨勢如圖10所示。當(dāng)模擬抽采時間逐漸增長時，相同位置處的滲透率逐漸增大，瓦斯壓力逐漸減小，原因是受煤層內(nèi)瓦斯壓力降低和時間效應(yīng)的影響。抽采1 d時，距鉆孔中心距離最遠處滲透率與瓦斯壓力數(shù)值分別為3.67×10-15m2,1.95 MPa，臨近鉆孔中心處滲透率與瓦斯壓力數(shù)值分別為13 ×10-15m2,0.387 MPa。而抽采30 d時，距鉆孔中心距離最遠處滲透率與瓦斯壓力數(shù)值分別為18.1 ×10-15m2,0.107 MPa，臨近鉆孔中心處滲透率與瓦斯壓力數(shù)值分別為18.2×10-15m2,0.1 MPa，整個模型近似達到一種內(nèi)外壓力平衡狀態(tài)。計算結(jié)果與實際深部煤層瓦斯從滲流通道運移至鉆孔的過程中參數(shù)的變化趨勢相符。

5 結(jié) 論

(1)結(jié)合平均有效應(yīng)力的影響，推導(dǎo)出考慮時效性影響的孔隙率和滲透率模型，并利用現(xiàn)場數(shù)據(jù)對模型有效性進行了驗證，確認模型能較好的模擬深部瓦斯長期運移規(guī)律。

(2)通過與未考慮時間效應(yīng)模型的計算結(jié)果對比可知，隨著抽采過程的進行，考慮時間效應(yīng)下的孔隙率、滲透率數(shù)值明顯減小，抽采時間越長，時間效應(yīng)越明顯，考慮時間效應(yīng)與否的兩種模型差值越大。

(3)運用考慮時間效應(yīng)的滲透率模型分析了瓦斯抽采過程中鉆孔周圍的滲流場，在本文模擬條件下，當(dāng)抽采時間為1 d時，臨近鉆孔中心處滲透率較大、瓦斯壓力較小;當(dāng)抽采時間達到30 d時，模型內(nèi)滲透率與瓦斯壓力的演化趨于平衡狀態(tài)，臨近與遠離鉆孔中心處的數(shù)值差別較小。