劉 健，馬 赟，張永國

(1.山西能源學(xué)院，山西 太原 030006；2.太原師范學(xué)院，山西 晉中 030619； 3.忻州師范學(xué)院，山西 忻州 034000；4.太原理工大學(xué) 礦業(yè)工程學(xué)院，山西 太原 030024)

近年來，我國的能源供需格局在不斷優(yōu)化，能源結(jié)構(gòu)在多元化發(fā)展[1]。而煤層氣作為新能源的代表在我國能源結(jié)構(gòu)中的地位正在逐步上升。煤層氣的主要成分是CH4，在長期的地質(zhì)作用下煤層氣與煤層是伴生存在的[2]。煤體內(nèi)部的孔洞和裂隙結(jié)構(gòu)為瓦斯儲(chǔ)存和流動(dòng)的通道，自然狀態(tài)下煤體內(nèi)部發(fā)育孔洞裂隙較少，滲透性極低，抽采瓦斯的難度較高[3，4]。而增大煤體的滲透性可以減小瓦斯抽采的難度，提高瓦斯抽采的效率。有研究表明[5，6]，對(duì)地層中賦存煤體滲透特性起主要影響作用的因素為溫度、靜水壓力和孔隙壓力，但針對(duì)各個(gè)因素的影響程度還未做過相應(yīng)研究。根據(jù)影響煤體滲透性的最主要因素可以為煤層氣高效開采方案的制定提供依據(jù)，故本文以煤體溫度、靜水壓力以及孔隙壓力為因素，將煤體滲透率作為研究指標(biāo)，基于正交試驗(yàn)法對(duì)不同條件下煤體的滲透率進(jìn)行了實(shí)驗(yàn)室測試，研究影響煤體滲透率的主控因素和煤體滲透率達(dá)到最大時(shí)三個(gè)參數(shù)的取值，以期為煤層氣開采最佳條件的選擇提供理論基礎(chǔ)。

1 試驗(yàn)介紹

1.1 試驗(yàn)方法

在進(jìn)行煤體的滲透率測試時(shí)，將標(biāo)準(zhǔn)煤樣置于反應(yīng)釜內(nèi)，對(duì)反應(yīng)釜進(jìn)行水浴加熱，通過溫度傳感器可得到水的實(shí)時(shí)溫度；利用驅(qū)替泵可提供煤樣所處的靜水壓力；將氮?dú)庾鳛榱鲃?dòng)氣體來測試煤樣滲透率，通過計(jì)量泵可控制進(jìn)口處氮?dú)鈮毫?，也就是煤樣的孔隙壓力[7]。為了模擬煤層真實(shí)的賦存環(huán)境，設(shè)計(jì)了三因素三水平正交試驗(yàn)，方案見表1。選擇測試溫度分別為20℃、30℃、和40℃，靜水壓力分別為7MPa、9MPa和11MPa，孔隙壓力分別為1MPa、2MPa和3MPa。

表1 不同測試方案下煤體滲透率統(tǒng)計(jì)表

1.2 試驗(yàn)原理

通過氮?dú)鈦頊y試煤體滲透率，主要原因如下：①無論煤體所處的環(huán)境如何改變，煤體均不會(huì)與氮?dú)獍l(fā)生化學(xué)反應(yīng)；②煤體對(duì)氮?dú)獠痪哂形叫?，則氮?dú)庠诿后w內(nèi)部的裂隙中流動(dòng)時(shí)不會(huì)發(fā)生滯留，所測煤體滲透率的結(jié)果更加精確。

基于達(dá)西定律[6]對(duì)煤體的滲透率進(jìn)行計(jì)算，氮?dú)庠诿后w內(nèi)的滲透狀態(tài)可通過公式(1)表示：

其中：Q為氮?dú)獾牧髁?，cc/s；μ為氮?dú)獾酿ざ认禂?shù)，Pa·s；L為煤體長度，cm；k為試樣的滲透系數(shù)；A為煤體截面尺寸，cm2；h0和h1分別為出口處和進(jìn)口處的流體高度，cm。

在式(1)中，有P=ρgh，則可得：

式中，P0為出口處的氮?dú)鈮毫?，Pa；P1為進(jìn)口處的氮?dú)鈮毫Γ彩敲后w的孔隙壓力，Pa。

令煤體的滲透率為K，則K=k/ρg，由此可得：

理想狀態(tài)下，出口處和進(jìn)口處氮?dú)獾牧髁繚M足：

P0Q0=P1Q1=PQ

結(jié)合式(1)—(4)可得煤體滲透率[6]計(jì)算公式：

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

煤體滲透率測試時(shí)，可得到氮?dú)饬髁亢兔后w孔隙壓力，將這兩個(gè)參數(shù)代入式(5)可得到不同試驗(yàn)方案下煤體的孔隙率。測試結(jié)果見表1。從表1中可以看出，方案3下煤體的滲透率最小，方案7下煤體的滲透率最大，達(dá)到了方案3的將近9倍，可見滲透率得到了質(zhì)的提高。究其原因：煤體是不均勻介質(zhì)，內(nèi)部存在強(qiáng)度不等的結(jié)構(gòu)面[8，9]，隨著溫度的升高，煤體發(fā)生熱膨脹，但煤體內(nèi)部礦物的熱膨脹系數(shù)是不同的，這就導(dǎo)致了煤體的熱損傷，拓寬了原有的裂縫，晶體間的熱破裂也開始形成，煤體的滲透性增強(qiáng)；煤體所受的靜水壓力會(huì)限制煤體的變形，阻礙煤體的熱膨脹，在較高的靜水壓力作用下煤體內(nèi)部的裂隙逐漸發(fā)生閉合，一定程度上堵塞的滲流通道。故當(dāng)煤體的溫度較高且所受的靜水壓力較小時(shí)，煤體的熱膨脹所受的阻力較小，滲透率相對(duì)較高。對(duì)影響煤體滲透率的3個(gè)主要因素進(jìn)行極差分析，得到結(jié)果見表2。

表2 煤體滲透率極差分析結(jié)果統(tǒng)計(jì)表

從表2中可知，煤體滲透率隨著溫度和孔隙壓力的增大而增大，隨著靜水壓力的升高而減小。靜水壓力對(duì)煤體滲透性的影響效果最為顯著，其次為孔隙壓力和溫度，具體表現(xiàn)為：靜水壓力>孔隙壓力≈溫度。

為了定量評(píng)價(jià)三個(gè)因素對(duì)煤體滲透率的影響效果，通過MATLAB進(jìn)行了多元回歸分析，得到：

S=0.3A+0.0045B-0.0443C+0.0478 (6)

式中，S為煤體滲透率，mD；A為煤體溫度，℃；B為靜水壓力，MPa；C為孔隙壓力，MPa。

該公式的相關(guān)系數(shù)為0.9351，與F檢驗(yàn)中對(duì)應(yīng)的概率P=0.0017，說明煤體滲透率與這三個(gè)參數(shù)之間具有很好的線性相關(guān)關(guān)系[10]。在得知煤體上溫度、靜水壓力以及孔隙壓力這三個(gè)參數(shù)的基礎(chǔ)上，可以通過式(6)來較為精確的預(yù)測煤體滲透率，所得結(jié)果與實(shí)際值較為接近。

3 不同孔隙壓下滲透率分布特征

煤層的溫度和所受的應(yīng)力隨著煤層埋深的改變而改變，通過滲透率測試試驗(yàn)得到了溫度和靜水壓力不同時(shí)煤體的滲透率變化情況[11，12]。為了直觀地分析溫度和靜水壓力雙重作用下煤體滲透率的分布情況，通過origin軟件進(jìn)行繪制，得到了不同孔隙壓下滲透率的分布，如圖1所示。

圖1 不同孔隙壓下煤體滲透率分布云圖

在各個(gè)滲透率分布云圖中，不同的滲透率通過不同的顏色反應(yīng)。當(dāng)煤體的孔隙壓不同時(shí)，滲透率分布的形態(tài)基本相似；滲透率從0.03mD增大到0.36mD的過程中，滲透率分布范圍先增大后減小，幾乎呈正態(tài)分布。孔隙壓力為1MPa時(shí)，煤體滲透率主要分布在0.15～0.24mD；孔隙壓力為2MPa時(shí)，煤體滲透率主要分布在0.12～0.27mD；孔隙壓力為2MPa時(shí)，煤體滲透率同樣主要分布在0.12～0.27mD。綜合分析，認(rèn)為當(dāng)孔隙壓力大于2MPa時(shí)，孔隙壓力對(duì)煤體滲透率的影響不明顯。

4 結(jié) 論

1)溫度升高會(huì)引起煤體的熱膨脹，裂縫會(huì)在孔隙壓的作用下發(fā)生擴(kuò)展，而靜水壓力會(huì)限制煤體的變形；滲透率測試結(jié)果表明，當(dāng)煤體溫度為40℃、靜水壓力為11MPa以及孔隙壓力為3MPa時(shí)，煤體的滲透率最大，達(dá)到了0.3617mD。

2)不同因素對(duì)煤體滲透性的影響效果表現(xiàn)為：靜水壓力>孔隙壓力≈溫度，而且煤體滲透率與這三個(gè)參數(shù)之間具有很好的線性相關(guān)關(guān)系。

3)滲透率從0.03mD增大到0.36mD的過程中，滲透率分布范圍先增大后減小，幾乎呈正態(tài)分布；當(dāng)孔隙壓力大于2MPa時(shí)，孔隙壓力對(duì)煤體滲透率的影響不明顯。