摘要：煤層滲透率是煤層氣開采的重要因素，為了有效的計(jì)算沁水盆地柿莊地區(qū)滲透率，根據(jù)該地區(qū)煤層裂隙的地質(zhì)特點(diǎn)和煤層瓦斯?jié)B流原理，采用達(dá)西定律推導(dǎo)的裂縫滲透率計(jì)算公式并增加調(diào)節(jié)系數(shù)做改進(jìn)，建立了分別適用于3#和15#煤層的滲透率計(jì)算公式。利用該地區(qū)井3#和15#煤層測井?dāng)?shù)據(jù)驗(yàn)證公式的有效性，計(jì)算結(jié)果與原始滲透率的相對誤差總體在30%左右。實(shí)驗(yàn)結(jié)果表明：改進(jìn)的達(dá)西方程對于柿莊地區(qū)煤層滲透率計(jì)算精度高，具有一定的地區(qū)實(shí)用性。

關(guān)鍵詞：滲透率 達(dá)西定律 柿莊地區(qū) 煤層氣

中圖分類號：TD31 文獻(xiàn)標(biāo)識碼：A 文章編號：1674-098X（2015）08（b）-0056-03

Abstract：Coal reservoir permeability is an important factor of coal-bed gas exploitation.In order to permeability is calculated effectively in ShiZhuang area of QinShui basin. According to the geological characteristics of coal fracture and seepage theory of mashgas in that area，which using the formula of fracture permeability derived from the Darcy’s law and added adjustment coefficient for improvement.Calculation formula of permeability is established by suitable for 3# and 15# coal seam respectively.Used logging data of 3# and 15# coal seam to verify the availability of the formula，which the relative error calculated general at around 30%.The experimental results show that coal reservoir permeability of Shi Zhuang area is calculated by the improved Darcy equation that had precision more higher and regional practicability.

Key Word：permeability；Darcy equation；Shi Zhuang Area；Coal-Seam gas

沁水盆地柿莊地區(qū)煤層氣的合理開采具有提高煤礦生產(chǎn)安全和增加能源產(chǎn)量的雙重意義。該地區(qū)煤層滲透率是煤層氣勘探和解吸的重要影響因素，目前尚未有統(tǒng)一的計(jì)算方法。國內(nèi)學(xué)者對于煤層滲透率計(jì)算方法主要有多元回歸，支持向量機(jī)，神經(jīng)網(wǎng)絡(luò)等統(tǒng)計(jì)方法[1-3]進(jìn)行預(yù)測分析，但柿莊地區(qū)地質(zhì)結(jié)構(gòu)復(fù)雜[4]，煤層滲透率普遍偏小，規(guī)律性波動較大，導(dǎo)致統(tǒng)計(jì)方法預(yù)測精度不高，并且需根據(jù)大量數(shù)據(jù)建模分析容易導(dǎo)致計(jì)算量大，過學(xué)習(xí)和欠學(xué)習(xí)等問題。煤層氣滲流力學(xué)隨著理論不斷深化，應(yīng)用范圍越來越廣，受到研究者們的廣泛關(guān)注。煤層中存在相互流通的孔隙和裂隙，我國周世寧[5-6]院士指出煤層氣在大孔隙和裂隙中的線性滲流原理基本符合達(dá)西定律，之后孫培德[7-9]又完善了線性流動數(shù)學(xué)模型。該理論應(yīng)用在柿莊地區(qū)需要參數(shù)變型和根據(jù)地質(zhì)環(huán)境特征進(jìn)行系數(shù)調(diào)整。該文是利用改進(jìn)的達(dá)西方程對柿莊地區(qū)煤層滲透率實(shí)際值進(jìn)行計(jì)算分析，可以避免統(tǒng)計(jì)方法帶來的問題，提高了預(yù)測精度，為柿莊地區(qū)煤層氣的開采提供良好的基礎(chǔ)。

1 沁水盆地區(qū)域地質(zhì)特征

沁水盆地柿莊地區(qū)地質(zhì)結(jié)構(gòu)復(fù)雜，其中山西組（3#）煤層由薄層石灰?guī)r，泥巖夾煤層以及濱海相砂巖組成，煤為黑色塊狀，具垂直節(jié)理較多，內(nèi)生裂隙較發(fā)育，組厚39～111 m；太原組（15#）由海陸交替石灰?guī)r、相砂巖、泥巖和煤層組成，煤層為階梯狀橫向斷口，內(nèi)生裂隙發(fā)育，組厚51～149米[10]。15#煤層位于3#煤層下部，依據(jù)地質(zhì)結(jié)構(gòu)可以判斷3#煤層垂直裂縫較多，15#煤層水平裂縫較多，是煤層氣開采的主力煤層。

2 煤層氣滲流特點(diǎn)

煤層結(jié)構(gòu)相對于巖層結(jié)構(gòu)較疏松，內(nèi)部充滿了大孔隙以及裂隙。 煤層氣以吸附態(tài)和游離態(tài)存在于煤層孔隙和裂隙中。吸附態(tài)的煤層氣會受到游離態(tài)的游離態(tài)瓦斯的影響而擴(kuò)散到煤粒表面，繼而進(jìn)入到連通的大孔隙以及裂隙中去。游離態(tài)的煤層氣體在大孔隙裂隙中向低壓處流動或者溶入地層水中而隨之流動。煤層氣在煤層中的滲流為主要存在狀態(tài)，流動過程不僅受到壓力影響，同時(shí)也受到煤壁的粘著力影響。

3 改進(jìn)的達(dá)西方程基本原理

煤層中水和氣是以層流狀態(tài)流動，在煤體孔隙中滲透時(shí)，由于滲透阻力的作用，沿程必然伴隨著能量的損失，而滲透能量損失與滲流速度之間的相互關(guān)系即為達(dá)西定律。所以根據(jù)煤層中水氣流動滿足達(dá)西定律煤巖流量計(jì)算公式可推導(dǎo)出裂縫滲透率計(jì)算公式[11-12]：

（1）

其中為有效孔隙度，為煤層裂縫。煤層層間裂縫計(jì)算公式由斯倫貝謝公司的Faivre和Sibbit兩位資深學(xué)者根據(jù)二維有限元方法進(jìn)行數(shù)值研究創(chuàng)立的雙側(cè)向測井解釋計(jì)算方法[13-14]，煤層裂縫公式如下：

垂直裂縫：

（2）

水平裂縫：

（3）

其中，為分別為淺/深側(cè)向電阻率，為泥漿電導(dǎo)率。

，為泥漿電阻率。

沁水盆地柿莊地區(qū)煤層分為山西組（3#）和太原組（15#）。根據(jù)該地區(qū)山西組（3#）煤層結(jié)構(gòu)選擇將式（2）垂直裂縫表達(dá)式帶入式（1），利用測井?dāng)?shù)據(jù)分析對方程進(jìn)行調(diào)整并加入調(diào)節(jié)系數(shù)0.5723，整理出柿莊地區(qū)山西組（3#）煤層滲透率計(jì)算公式：

（4）

同理根據(jù)該地區(qū)太原組（15#）煤層結(jié)構(gòu)選擇將式（3）水平裂縫表達(dá)式帶入式（1），利用測井?dāng)?shù)據(jù)分析對方程進(jìn)行調(diào)整并加入調(diào)節(jié)系數(shù)0.2741，整理出柿莊地區(qū)太原組（15#）煤層滲透率計(jì)算公式：

（5）

4 滲透率計(jì)算實(shí)驗(yàn)

4.1 研究區(qū)3#煤層滲透率計(jì)算

利用柿莊地區(qū)SX-005，SX-009，SX-011等10口井3#煤層的測井?dāng)?shù)據(jù)進(jìn)行（4）式改進(jìn)達(dá)西滲透率方程的驗(yàn)證。每組數(shù)據(jù)包括滲透率PERM，有效孔隙度POR，泥漿電阻率Rm，深側(cè)向電阻率Rt，淺側(cè)向電阻率Rs，其中滲透率PERM作為因變量，其余參數(shù)作為自變量。將自變量帶入方程中求出滲透率計(jì)算值并與原始滲透率進(jìn)行對比。表1列出計(jì)算值與真實(shí)值部分?jǐn)?shù)據(jù)的誤差分析結(jié)果，其中EA為絕對誤差，EB/%為相對誤差。由表1可知，相對誤差平均值在0.015左右，絕對誤差平均值在12.5%左右。

圖1是柿莊地區(qū)井3#煤層測井?dāng)?shù)據(jù)中部分樣本的達(dá)西方程計(jì)算值和真實(shí)值對比圖。由圖1顯示這10口井的計(jì)算結(jié)果與對應(yīng)的滲透率真實(shí)值基本保持一致。

4.2 研究區(qū)15#煤層滲透率計(jì)算

利用柿莊地區(qū)SX-008，SX-010，SX-015等11口井15#煤層的真實(shí)測井?dāng)?shù)據(jù)對（5）式改進(jìn)的達(dá)西滲透率方程做驗(yàn)證。表2列出計(jì)算值與真實(shí)值部分?jǐn)?shù)據(jù)的誤差分析結(jié)果，由表2可知，相對誤差平均值在0.0242左右，絕對誤差17.7818%左右。

圖2是柿莊地區(qū)井15#煤層測井?dāng)?shù)據(jù)中部分樣本的滲透率達(dá)西方程計(jì)算值和真實(shí)值對比圖。由圖2顯示雖然有個(gè)別點(diǎn)存在偏差較大的情況，但總體這11口井的滲透率計(jì)算結(jié)果與真實(shí)值相接近。通過對3#和15#煤層驗(yàn)證結(jié)果表明：改進(jìn)的達(dá)西方程對于該地區(qū)煤層滲透率計(jì)算結(jié)果較理想，具有一定的合理性。

5 實(shí)驗(yàn)結(jié)果分析

實(shí)際測井過程中由于地質(zhì)結(jié)構(gòu)變化以及設(shè)備誤差，難免會有一些數(shù)據(jù)點(diǎn)異常，剔除數(shù)據(jù)中的奇異點(diǎn)對實(shí)驗(yàn)結(jié)果的正確性很有必要。該文在實(shí)驗(yàn)過程中對于煤層數(shù)據(jù)進(jìn)行了篩選，對于致密層的數(shù)據(jù)點(diǎn)進(jìn)行了排除。經(jīng)過研究發(fā)現(xiàn)：合理的測井滲透率數(shù)據(jù)越多，越能反映煤儲層的實(shí)際情況。

實(shí)驗(yàn)通常采用絕對誤差和相對誤差兩個(gè)評判標(biāo)準(zhǔn)評價(jià)計(jì)算結(jié)果的準(zhǔn)確度，改進(jìn)的達(dá)西方程對于柿莊地區(qū)共計(jì)21口井的3#和15#煤層滲透率計(jì)算相對誤差小于0.03，絕對誤差基本在25%以內(nèi)。雖然有部分?jǐn)?shù)據(jù)點(diǎn)對于原始滲透率誤差較大，但計(jì)算結(jié)果的總體趨勢與真實(shí)值相吻合，體現(xiàn)了改進(jìn)的達(dá)西方程在該地區(qū)煤層滲透率計(jì)算上的有效性。

6 結(jié)論

改進(jìn)的達(dá)西方程對于柿莊地區(qū)大孔隙和裂縫多的煤層段具有更好的計(jì)算效果。

該文實(shí)驗(yàn)建立的改進(jìn)達(dá)西方程只適合應(yīng)用于柿莊地區(qū)的滲透率計(jì)算，并且這只是該地區(qū)井下的普遍規(guī)律，對于個(gè)別地質(zhì)復(fù)雜井的滲透率計(jì)算，還需要根據(jù)實(shí)際情況增加新的調(diào)節(jié)系數(shù)。

改進(jìn)的達(dá)西方程是根據(jù)測井參數(shù)計(jì)算實(shí)際滲透率的值，為柿莊地區(qū)的煤層氣開采提供依據(jù)。

參考文獻(xiàn)

[1]王相業(yè)，李建武，楊志遠(yuǎn)，等.柳林地區(qū)煤層滲透率逐步回歸分析與預(yù)測[J].煤田地質(zhì)與勘探，2013（3）：18-22.

[2]尹光志，李銘輝，李文璞，等.基于改進(jìn)BP神經(jīng)網(wǎng)絡(luò)的煤體瓦斯?jié)B透率預(yù)測模型[J].煤炭學(xué)報(bào)，2013（7）：1179-1184.

[3]王雷.支持向量機(jī)預(yù)測煤儲層滲透性[D].中國石油大學(xué)，2011.

[4]王存武，柳迎紅，梁建設(shè)，等.沁水盆地南部柿莊北地區(qū)煤層氣勘探潛力研究[J].中國煤層氣，2014（3）：3-6.

[5]周世寧，孫輯正.煤層瓦斯流動理論及其應(yīng)用[J].煤炭學(xué)報(bào)，1965（1）：24-37.

[6]周世寧.瓦斯在煤層中流動的機(jī)理[J].煤炭學(xué)報(bào)，1990（1）：15-24.

[7]孫培德.煤層瓦斯流場流動規(guī)律的研究[J].煤炭學(xué)報(bào)，1987（4）：74-82.

[8]孫培德.煤層瓦斯流動方程補(bǔ)正[J].煤田地質(zhì)與勘探，1993（5）：34-35.

[9]Sun P.COAL GAS DYNAMICS AND ITS APPLICATIONS[J].Scientia Geologica Sinica，1994（1）：83-90.

[10]陸小霞，黃文輝，王佳旗，等.沁水盆地柿莊北深部煤層煤體結(jié)構(gòu)發(fā)育特征[J].煤田地質(zhì)與勘探，2014（3）：8-11，16.

[11]侯俊勝.煤層氣儲層測井評價(jià)方法及其應(yīng)用[M].北京：冶金工業(yè)出版社，2000.

[12]何曉輝.煤層氣測井解釋方法研究與應(yīng)用[D].長江大學(xué)，2012.

[13]桑孝偉.深部煤炭儲層參數(shù)反演預(yù)測[D].中國地質(zhì)大學(xué)（北京），2010.

[14]黃烈林，高純福，布志虹，等.雙側(cè)向測井確定裂縫等效寬度—— 兼論Sibbit公式中的一個(gè)錯(cuò)誤[J].江漢石油學(xué)院學(xué)報(bào)，2002（4）：42-44，4.