循環(huán)載荷下煤樣不同方向滲透特性試驗(yàn)研究*

賈立鋒，董 擎，梁 冰，孫維吉

(1.遼寧工程技術(shù)大學(xué)，礦業(yè)學(xué)院，遼寧 阜新 123000；2. 安陽工學(xué)院，機(jī)械工程學(xué)院，河南 安陽 455000；3. 遼寧工程技術(shù)大學(xué)，力學(xué)與工程學(xué)院，遼寧 阜新 123000)

0 引言

在煤層開采過程中，部分煤體受到采動(dòng)應(yīng)力的反復(fù)作用，例如分層開采的下伏煤層，起支撐和隔離作用的煤柱、反復(fù)水力壓裂的煤體等。煤體內(nèi)含有大量的孔隙裂隙，循環(huán)載荷作用下煤體內(nèi)部結(jié)構(gòu)會(huì)發(fā)生變化，影響煤體的滲透特性。煤是一種典型的沉積巖，具有各向異性的特征[1-3]，循環(huán)載荷作用下煤體不同方向的滲透情況均會(huì)發(fā)生變化。

目前循環(huán)載荷作用下煤體滲透的研究有很多[4-7]。潘榮錕等[5]對(duì)含層理原煤開展了滲透實(shí)驗(yàn)，結(jié)果表明隨著有效應(yīng)力的增大，滲透率降低，加載過程中煤體層理裂隙變形對(duì)裂隙面造成永久性損傷；魏建平等[6]開展了不同含水率條件下2次卸載圍壓的三軸滲流實(shí)驗(yàn)，第2次加載過程中滲透率較第1次平緩，2次卸載過程中，滲透率有一定程度升高，但恢復(fù)不到初始值；郭軍杰等[7]以原煤為研究對(duì)象，進(jìn)行等幅循環(huán)加卸載的滲透率測試，得到在滲透率降低階段，滲透率與循環(huán)次數(shù)呈冪指數(shù)函數(shù)關(guān)系，在滲透率升高階段，滲透率與循環(huán)次數(shù)呈指數(shù)函數(shù)關(guān)系。在煤滲透各向異性的研究方面[8-11]，王登科等[8]開展了含瓦斯煤的各向異性滲流規(guī)律的研究，結(jié)果表明，煤體瓦斯流動(dòng)具有非常明顯的各向異性特征，其滲透率與有效應(yīng)力之間符合負(fù)指數(shù)函數(shù)變化規(guī)律；田坤云等[9]對(duì)加卸載過程平行、斜交及垂直層理方向的原煤試件進(jìn)行了滲透實(shí)驗(yàn)，研究表明，加載過程中煤樣的滲透率與有效應(yīng)力成正比，卸載過程中滲透率與有效應(yīng)力之間成反比；趙宇等[11]在不同圍壓和氣體壓力下對(duì)煤樣的面割理、端割理、垂直層理方向上的滲透率進(jìn)行測試分析，研究認(rèn)為不同氣體壓力下滲透率隨著圍壓增大而減小，氣體壓力對(duì)端割理方向和垂直層理方向滲透率影響不大。

上述研究僅考慮了煤的各向異性滲透特征或者循環(huán)加載條件下單個(gè)方向的滲透特征，未將二者結(jié)合起來。本文以平頂山十二礦己15煤層煤樣為研究對(duì)象，利用自行研制的應(yīng)力-滲流-解吸煤體變形試驗(yàn)裝置，研究了循環(huán)圍壓加載下煤樣不同方向滲透特征。

1 試驗(yàn)方法

1.1 試樣的制備

煤樣取自于河南省平頂山十二礦己15煤層，取樣地點(diǎn)在31030工作面進(jìn)風(fēng)巷，屈服強(qiáng)度2.44～ 2.62 MPa，抗壓強(qiáng)度為4.86～6.37 MPa。將煤樣采集回來后，判斷煤的層理方向和裂隙發(fā)育情況，用巖石切割機(jī)制作尺寸為50 mm×50 mm×50 mm的層理方向正方體試樣，如圖1所示。盡可能避免試樣內(nèi)含有裂隙，最后用ZDM-200磨石機(jī)進(jìn)一步加工煤塊，保證試樣的邊長為50±1 mm。

圖1 試樣及其層理結(jié)構(gòu)特征Fig.1 Specimens and bedding structure characteristics

本試驗(yàn)共加工出2塊試樣，分別編號(hào)為PDS-1B和PDS-2，2塊試樣的層理面平行于面3，垂直于面1和面2，表面無明顯裂紋，而PDS-2試樣面1和面2交叉處下方有缺角。

1.2 試驗(yàn)裝置

試驗(yàn)使用的是自行研制的應(yīng)力-滲流-解吸煤體變形試驗(yàn)裝置(見圖2)，該裝置詳細(xì)介紹見文獻(xiàn)[12-13]，本試驗(yàn)未使用裝置中的變形測量系統(tǒng)。為了能夠測試試樣不同方向的滲透率，給煤樣夾持器添加了1個(gè)正方體壓頭。改進(jìn)前只能測試尺寸為50 mm×50 mm×100 mm的長方體試樣，改進(jìn)后可以對(duì)截面積50 mm×50 mm×50 mm的試樣進(jìn)行測試(見圖3)。滲透試驗(yàn)采用的氣體為純度99.99%的N2。

圖2 應(yīng)力-滲流-解吸煤體變形試驗(yàn)裝置Fig.2 Experimental device of coal matrix deformation caused by stress-seepage-desorption

圖3 夾具中的試樣Fig.3 A specimen in the gripper

1.3 試驗(yàn)原理及步驟

滲透率是表征多孔介質(zhì)滲透特性的重要參數(shù)。根據(jù)穩(wěn)態(tài)法測試煤巖樣滲透率原理，滲透率可以通過公式(1)計(jì)算：

(1)

式中：k為煤體的滲透率，m2；pout為煤體出口壓力，Pa；pin為煤體進(jìn)氣口壓力，Pa；Qout為試樣出口穩(wěn)定流量，m3/s；A為流體通過的截面積，m2；L為煤樣的長度，m；μ為氮?dú)獾膭?dòng)力黏度，17.58×10-6Pa·s。

試驗(yàn)研究循環(huán)載荷作用下煤樣不同方向上滲透特性，需要對(duì)試樣的3個(gè)面進(jìn)行測試，測試順序?yàn)榇怪庇诿?方向、垂直于面2方向、垂直于面3方向，每個(gè)面測試完畢后依次需要卸載軸壓和圍壓、拆卸試件更換測試面、封裝試件、施加軸壓、圍壓。

用PDS-1B試樣測試了不同方向滲透的初始值，為了避免圍壓和軸壓對(duì)滲透率造成干擾，圍壓和軸壓取值為2 MPa，小于煤體的屈服強(qiáng)度。進(jìn)氣口壓力以0.3 MPa為初始值，以0.3 MPa/次的增量加載至1.5 MPa結(jié)束。依次測試垂直于面1、面2、面3方向的滲透率。

用PDS-2試樣開展循環(huán)載荷下煤體不同方向滲透率試驗(yàn)。具體試驗(yàn)步驟如下：1)將試驗(yàn)裝置按圖2連接，檢查試驗(yàn)裝置的氣密性；2)將PDS-2煤試樣依圖3方式安裝，安裝時(shí)保證垂直于面1的方向?yàn)闈B透方向；3)將軸壓加載至3 MPa，圍壓加載至4 MPa，待軸壓和圍壓穩(wěn)定后加載孔隙壓力。進(jìn)氣口壓力(相對(duì)壓力值)以0.5 MPa為初始值，以0.5 MPa/次的增量加載至3 MPa，測試滲流穩(wěn)定時(shí)單位時(shí)間內(nèi)排出氣體的體積；4)將圍壓以1 MPa/次的增量由4 MPa加載至8 MPa，每個(gè)圍壓下均重復(fù)步驟3)中進(jìn)氣口壓力的加載路徑，并測滲流穩(wěn)定時(shí)單位時(shí)間內(nèi)排出氣體的體積；5)面1測試完畢后，重復(fù)步驟3)～5)對(duì)試樣的垂直于面2和面3方向的滲透率進(jìn)行測試。

2 試驗(yàn)結(jié)果分析

假設(shè)氣體壓力在試樣內(nèi)均勻分布，試樣沿滲透方向的孔隙壓力是不斷變化的，可以采用進(jìn)口壓力和出口壓力之和的平均值來壓計(jì)算孔隙壓力[14]，本文以平均孔隙壓力pm為指標(biāo)，分析其與滲透率的關(guān)系，如式(2)所示：

(2)

2.1 層理特征對(duì)煤體滲透率的影響

對(duì)PDS-1B不同方向滲透率測試結(jié)果分析可以得到平均孔隙壓力與滲透率的關(guān)系，如圖4所示。

圖4 PDS-1B不同方向滲透率隨平均孔隙壓力變化Fig.4 Variation of permeability in different directions with average pore pressure

由圖4可以看出，垂直于面1方向滲透率的范圍為0.96～1.22 mD，垂直于面2方向滲透率范圍為0.91～1.19 mD，垂直于面1和面2滲透方向?yàn)槠叫杏趯永砻娣较颍谄骄紫秹毫ο嗤瑫r(shí)滲透率相差不大。垂直于面3滲透的方向?yàn)榇怪庇趯永砻娣较颍瑵B透率范圍為0.54～0.74 mD，小于平行層理面的滲透率。在相同的軸壓、圍壓和平均孔隙壓力下，煤在垂直層理面方向的滲透率大于平行層理面，平行層理面內(nèi)的滲透率相差不大。

2.2 循環(huán)載荷下滲透壓差與流量的關(guān)系

滲透壓差是試樣進(jìn)氣口和出氣口的壓力差Δp，由于壓力表的讀數(shù)為相對(duì)壓力值，出口壓力為大氣壓，因此滲透壓差數(shù)值上等于進(jìn)口壓力表值。滲透壓力與流量的關(guān)系如圖5所示。

圖5 流量與滲透壓差之間的關(guān)系Fig.5 Relationship between rate of flow and osmotic pressure difference

由圖5可以看出，隨著滲透壓差的增大，通過煤樣的流量增加，且二者之間的關(guān)系可以用函數(shù)Q=D1(Δp2+0.2Δp)表示，其中D1為擬合系數(shù)，二次函數(shù)擬合出的相關(guān)性系數(shù)R2均大于0.999。如圖5(a)所示，垂直于面1滲透時(shí)，圍壓為4 MPa，滲透壓差為3 MPa時(shí)流量為26.49 mL/s；圍壓為5，6，7和8 MPa，滲透壓差為3 MPa時(shí)流量分別為12.41，8.35，5.30和4.11 mL/s，在相同滲透壓差下，通過試樣的流量隨著圍壓的增大而減小。如圖5(b)所示，垂直于面2滲透時(shí)，該方向流量隨壓差的變化規(guī)律與垂直于面1方向相同，但在相同圍壓和滲透壓差下，垂直于面2方向上通過試樣的流量大于垂于面1方向。垂直于面3滲透時(shí)，在相同滲透壓差下，圍壓對(duì)透過試樣流量影響比較復(fù)雜，在滲透壓差為1 MPa時(shí)，Qσ3=6 MPa>Qσ3=8 MPa>Qσ3=4 MPa；滲透壓差為2 MPa時(shí)，Qσ3=6 MPa>Qσ3=4 MPa>Qσ3=8 MPa。其原因可能是圍壓為6 MPa時(shí)，試樣內(nèi)有新裂隙產(chǎn)生，流量增大，而圍壓由6 MPa增加至8 MPa時(shí)裂隙閉合。

2.3 循環(huán)載荷下平均孔隙壓力與滲透率的關(guān)系

滲透率與平均孔隙壓力之間的關(guān)系如圖6所示。

圖6 滲透率與平均孔隙壓力之間的關(guān)系Fig.6 Relationship between permeability and average pore pressure

由圖6(a)可以看出，垂直于面1滲透時(shí)，圍壓為4 MPa時(shí)煤樣的滲透率隨著平均孔隙壓力的增加呈現(xiàn)先減小后增加的趨勢，平均孔隙壓力為1.6 MPa時(shí)的滲透率達(dá)到最大值0.196 mD，而圍壓為5，6，7和8 MPa時(shí)煤樣的滲透率隨平均孔隙壓力變化不大；從圖6(b)可以看出，垂直于面2滲透時(shí)，滲透率隨著平均孔隙壓力的增大而上下波動(dòng)，隨著圍壓的增大滲透率逐漸減小，但在相同的平均孔隙壓力和圍壓下，垂直于面2方向的滲透率明顯高于垂直于面1方向上的滲透率。垂直于面1和面2方向上滲透率隨平均孔隙壓力變化不大，滲透率只與多孔介質(zhì)本身的結(jié)構(gòu)特征有關(guān)，與流體特征無關(guān)[15]。

從圖6(c)可以看出，垂直于面3滲透時(shí)的變化規(guī)律與垂直于面1和面2時(shí)不同，首先滲透率隨平均孔隙壓力的增大而減小，其規(guī)律可以用考慮克林伯格效應(yīng)的滲透率公式表示：

(3)

式中：k0為試樣的絕對(duì)滲透率， mD；b為克林伯格系數(shù)。圍壓為4，6，8 MPa時(shí)的曲線擬合相關(guān)性系數(shù)分別為0.864，0.945，0.964，擬合度較高。說明垂直于面3滲透時(shí)，滲透率受克林伯格效應(yīng)的影響。其次，相同滲透壓差下，隨著圍壓的增加，滲透率變化規(guī)律復(fù)雜，和圍壓對(duì)透過試樣流量的影響類似。

2.4 循環(huán)載荷對(duì)不同方向滲透率的影響

將同一圍壓不同平均孔隙壓力下的滲透率視為1組數(shù)據(jù)，求出其算數(shù)平均值。將PDS-1B在圍壓2 MPa時(shí)的算術(shù)平均值作為滲透率初始值，其垂直于面1，面2，面3方向的初始滲透率分別為1.085，1.005和0.627 mD。將PDS-2垂直于面1，面2，面3方向的不同圍壓下滲透率的算術(shù)平均值與初始滲透率值比值(下文簡稱滲透率比值)作為分析指標(biāo)，建立其與圍壓的關(guān)系，如圖7所示。雖然垂直于面3方向滲透率隨平均孔隙壓力變化很大，為了分析循環(huán)圍壓加載對(duì)滲透率的影響，也對(duì)其求算術(shù)平均值。

圖7 不同方向滲透率比值與圍壓之間的關(guān)系Fig.7 Relationship between confining pressure and permeability ratio in different directions

由圖7可知，垂直于面1和面2方向?yàn)槠叫杏趯永矸较颍跏紳B透率在1 mD左右，第1次加載和第2次加載時(shí)滲透率比值均小于0.3，這是因?yàn)閲鷫涸黾樱瑢?dǎo)致煤體內(nèi)裂隙閉合，滲透率變小。在相同的圍壓下垂直于面3的滲透率比值大于垂直于面2同時(shí)也大于垂直于面1的滲透率比值，這是因?yàn)檠h(huán)圍壓大于煤樣的最大屈服強(qiáng)度和單軸抗壓強(qiáng)度，煤樣內(nèi)部垂直于面2經(jīng)歷了1次加載，面3經(jīng)歷了2次加載，裂隙得到擴(kuò)展。說明滲透率比值隨循環(huán)加載次數(shù)的增加而增大。煤的初始滲透率為垂直于面1，面2，面3時(shí)的初始滲透率，分別為1.085，1.005和0.627 mD，大小順序垂直于面1大于垂直于面2，大于垂直于面3，經(jīng)歷3次循環(huán)加載后，滲透率比值垂直于面3方向大于垂直于面2，大于垂直于面1。

隨著圍壓的增加，垂直于面1和面2方向的滲透率比值逐漸減小，這是因?yàn)閲鷫涸黾樱紫读严吨饾u閉合。垂直于面3方向滲透率比值在圍壓6 MPa時(shí)較大，這可能是由于在圍壓6 MPa時(shí)煤體破壞，較4 MPa時(shí)平均滲透率增加，圍壓為8 MPa時(shí)裂隙較6 MPa時(shí)閉合，滲透率比值減小。

3 結(jié)論

1)在相同的軸壓、圍壓和平均孔隙壓力下，平行層理面方向的滲透率大于垂直層理，平行層理面內(nèi)的滲透率相差不大。

2)通過煤樣的流量隨著滲透壓差的增大而增加，且二者之間的關(guān)系可以用Q=D1(Δp2+0.2Δp)表示。

3)PDS-2垂直于面1和面2滲透時(shí)，隨著平均孔隙壓力的增加，滲透率變化幅度不大。垂直于面3滲透時(shí)，滲透率隨平均孔隙壓力的增大而減小，滲透率與平均孔隙壓力之間的關(guān)系可以考慮用克林伯格效應(yīng)的滲透率公式描述。

4)圍壓增加，導(dǎo)致裂隙閉合，滲透率減小，當(dāng)循環(huán)圍壓大于煤屈服強(qiáng)度和抗壓強(qiáng)度時(shí)，裂隙擴(kuò)展，滲透率增加。循環(huán)圍壓加載可以改變煤樣原有不同方向的滲透率大小順序，滲透率與原初始滲透率比值隨循環(huán)加載次數(shù)的增加而增大。