馬平華，邵先杰，霍夢穎，褚慶忠，霍春亮，梁武斌(.燕山大學(xué) 石油工程系，河北 秦皇島 066004； .中海石油有限公司 天津分公司，天津 30045)

近20年來，我國煤層氣開發(fā)經(jīng)歷了礦場試驗到規(guī)模化開發(fā)，獲得了突破性進展。根據(jù)統(tǒng)計，截止到2017年底，經(jīng)過地質(zhì)與石油工作者的不懈努力，累計探明煤層氣地質(zhì)儲量超過了30×1012m3，已完鉆煤層氣井超過14 000口，年產(chǎn)量達到了70.2×108m3，取得了巨大成績。但有相當大的一部分井為低產(chǎn)井，穩(wěn)產(chǎn)時間短，影響整體開發(fā)效益，目前迫切需要提高煤層氣井的產(chǎn)能[1]。

國家依托“973”項目和科技重大專項，開展了大量的基礎(chǔ)性研究。針對中國的中高煤階，在煤儲層基本地質(zhì)特征、微觀孔隙結(jié)構(gòu)、煤層氣富集、解吸和滲流等基礎(chǔ)理論方面取得了很多成果；在鉆井、采氣和壓裂等工藝方面形成了相對完善的礦場技術(shù)。但是在煤儲層地質(zhì)建模方面明顯滯后于礦場生產(chǎn)，甚至在這方面是一項空白。地質(zhì)建模是氣藏數(shù)值模擬和開發(fā)動態(tài)分析的基礎(chǔ)，沒有可靠、準確、能真實反映煤儲層特征的地質(zhì)模型根本就得不出可信的數(shù)值模擬結(jié)果，嚴重影響了開發(fā)方案編制和開發(fā)措施的實施，從而影響到了煤層氣井的產(chǎn)量。

煤儲層地質(zhì)建模存在的主要困難是：①煤巖樣易碎，很難通過室內(nèi)實驗獲得煤層真實的滲透率；②煤儲層具有雙孔隙介質(zhì)的特性，儲層物性受煤巖組構(gòu)、割理和微裂隙等因素影響[2-4]，存在高度的非均質(zhì)性，即使做出個別樣品的實驗結(jié)果，也不能完全代表儲層的物性；③煤巖具有彈性模量低，泊松比高的特點，對應(yīng)力十分敏感，地面條件下測試的滲透率不能反映地下的情況；④煤巖滲透率極低，與孔隙度之間相關(guān)性不明顯，砂巖儲層中常用的測井物性解釋方法在煤儲層滲透率解釋中很難應(yīng)用。

國外煤儲層厚度大[5]、滲透率高、相對均質(zhì)，基本上是采用砂巖儲層的建模方法。國內(nèi)的礦場生產(chǎn)實踐證明，砂巖儲層的建模方法對于我國中高煤階和超低滲煤儲層不適應(yīng)。因此，迫切需要建立一套完善的、在煤巖組構(gòu)和地應(yīng)力約束下，能夠反映煤儲層特性的地質(zhì)建模理論與方法，為煤層氣的高效開發(fā)奠定基礎(chǔ)。

1 煤儲層地質(zhì)建模思路

煤儲層滲透率是成煤物質(zhì)在成巖作用下經(jīng)壓實、失水、收縮、有機質(zhì)演化和破裂等地質(zhì)過程綜合作用的結(jié)果，當前地質(zhì)條件下的地應(yīng)力對煤巖的滲透率影響很大，因為當?shù)貞?yīng)力高時，煤巖基質(zhì)顆粒堆積緊密，割理和微裂隙閉合，滲透率往往比較低，反之則高。

在長期的科研生產(chǎn)實踐中，我們探討了煤儲層新的地質(zhì)建模思路和方法(圖1)。

1) 利用試井資料確定滲透率。韓城礦區(qū)部分井投產(chǎn)前進行了試井，通過試井資料的分析，計算出滲透率，這一結(jié)果基本反映了煤層地下的真實滲透率。

2) 依據(jù)壓裂資料確定地應(yīng)力。煤層氣井(直井)投產(chǎn)前都要實施壓裂，利用壓裂資料計算地應(yīng)力的大小。

3) 煤巖類型測井解釋。不同類型煤巖的巖石力學(xué)參數(shù)不同，密度不同，測井響應(yīng)不同。通過煤巖測井響應(yīng)關(guān)系的研究，建立煤巖識別圖板，定性識別出各類煤巖，為煤巖約束滲透率奠定基礎(chǔ)。

4) 滲透率與地應(yīng)力定量關(guān)系的建立。煤儲層對地應(yīng)力比較敏感，通過其關(guān)系的定量研究，為地應(yīng)力約束煤巖滲透率提供依據(jù)。

5) 煤儲層空間幾何分布模型建立。在鉆井和沉積相約束下采用隨機建模技術(shù)建立煤層空間幾何分布模型。

6) 煤巖相(類型)模型建立。利用測井煤巖類型識別技術(shù)和圖板，逐井逐層識別，在井點約束下采用示性點過程模擬方法建立煤巖相模型。

7) 地應(yīng)力模型建立。在構(gòu)造背景約束下，通過井點的控制，采用隨機建模技術(shù)建立地應(yīng)力模型。

圖1 多參數(shù)約束下的煤儲層隨機地質(zhì)建模技術(shù)路線Fig.1 Technology roadmap for stochastic geology modeling of coalbed under the constraints of multi-parameters

8) 滲透率模型建立。有試井數(shù)據(jù)的井點，滲透率以試井數(shù)據(jù)控制，沒有試井但有壓裂的井，利用地應(yīng)力計算的滲透率控制井點滲透率?？臻g上在沉積相、煤巖相模型和地應(yīng)力模型的約束下，采用隨機建模技術(shù)建立煤儲層滲透率模型。

2 地質(zhì)概況

2.1 地質(zhì)背景

韓城地區(qū)位于鄂爾多斯盆地東南緣，南鄰渭河斷陷，東接汾河地塹，西與鄂爾多斯盆地西緣逆沖構(gòu)造帶相毗鄰，面積約為1 100 km2。由于受韓城和渭河兩條邊界斷層的影響，礦區(qū)邊緣地帶和淺部，構(gòu)造復(fù)雜，斷裂發(fā)育；中深部地層則較為平緩，僅存在少量寬緩褶皺起伏，總體為單斜構(gòu)造形態(tài)，傾向北西，淺部傾角較大，深部較平緩，一般為5°～10°。礦區(qū)的北部主要受擠壓作用的控制，南部受拉伸作用控制。煤層氣的主產(chǎn)層位為山西組發(fā)育的3#煤層和太原組的5#和11#煤層[6]。

3#，5#和11#煤層埋藏深度一般為400～1 000 m，地層壓力為2.42～2.8 0 MPa，壓力系數(shù)為0.6～0.8。變質(zhì)程度較高，以瘦煤和貧煤為主。實驗樣品的滲透率為0.01×10-3～2.5×10-3μm2。含氣量為9.78～11.23 m3/t，平均資源豐度為2.91×108m3/km2。

2.2 沉積環(huán)境與沉積相

韓城礦區(qū)的3#煤層發(fā)育相對較薄，厚度在0.8～4.4 m；5#煤層厚度中等，為0.5～9.1 m，中部夾一比較穩(wěn)定的夾矸；11#煤層厚度大，并且相對穩(wěn)定，為4.0～11.6 m。

依據(jù)研究區(qū)地質(zhì)背景資料，通過巖心觀察和對測井資料、巖性特征、巖性組合以及沉積構(gòu)造的分析，認為研究區(qū)太原組屬于潟湖-潮坪-障壁島沉積體系，可進一步劃分為潟湖亞相、沼澤亞相、潮坪亞相和障壁島亞相。5#和11#煤層屬濱海沼澤環(huán)境形成的泥炭堆積，環(huán)境相對穩(wěn)定，厚度較大(圖2)，總體趨勢是東北方向厚度大，西南方向厚度減??；煤巖顯微組分中硫化鐵類礦物和碳酸鹽類礦物含量偏高，平均值分別為0.85%和5.75%。山西組為河控三角洲沉積體系[7]，包括三角洲平原分流河道、分流間、間灣沼澤和決口扇等微相類型。3#煤層屬于三角洲平原間灣沼澤環(huán)境堆積，由于環(huán)境不穩(wěn)定，煤層厚度較薄，總體上呈不規(guī)則的條帶狀分布；煤巖顯微組分中硫化鐵類礦物和碳酸鹽類礦物含量較其它兩層低，平均分別為0.2%和3.2%。

3 地應(yīng)力分布規(guī)律

3.1 地應(yīng)力計算

地應(yīng)力的測試方法很多，主要有差應(yīng)變法、Kaiser效應(yīng)法、井眼崩落法和地球物理法等，但對于煤層氣井來講，幾乎每口井都壓裂投產(chǎn)，壓裂資料豐富，利用壓裂資料計算地應(yīng)力是最便捷有效的方法。

通過水力壓裂施工曲線以及壓裂施工過程中的壓力變化等參數(shù)，可直接讀取出煤巖的破裂壓力、裂縫延伸壓力和瞬時停泵壓力[8]。而閉合壓力可通過壓裂施工過程中壓力變化數(shù)據(jù)求取。突然停止向井中泵入液體，此時已經(jīng)壓開的裂縫由于地應(yīng)力的作用而開始閉合，根據(jù)閉合壓力可以計算出最小水平主應(yīng)力[7]。再根據(jù)巖石的抗張強度與破裂壓力，可以進一步求出最大水平主應(yīng)力[9]。其關(guān)系式如下：

式中：σh為最小水平主應(yīng)力，MPa；σH為最大水平主應(yīng)力，MPa；pw為井筒液柱壓力，MPa；pc為裂縫閉合壓力，MPa；pm為井筒摩阻，MPa；pf為煤巖石破裂壓力，MPa；po為煤層孔隙壓力，MPa；St為煤巖的抗張強度，MPa。

根據(jù)上述方法計算出了韓城礦區(qū)所有試井層位的地應(yīng)力值，共計118個井層。

3.2 地應(yīng)力分布

根據(jù)計算結(jié)果，不同井以及不同層之間的差別很大，最小水平主應(yīng)力最大的井是H3-041井的3#煤層，38.92 MPa；最小值是H3-037井11#煤層，7.1 MPa。同一口井不同層位間也存在比較大的差別，比如H3-032井，3#,5#和11#煤層的最小水平主應(yīng)力分別為13.67,16.73和28.55 MPa?？傮w趨勢是隨著深度增加地應(yīng)力呈增大的趨勢，但具體到某一構(gòu)造部位的井并不全是這樣的規(guī)律。

根據(jù)繪制的最小水平主應(yīng)力圖分析，3#煤層高應(yīng)力區(qū)分布在西南和東北部，東南部為低地應(yīng)力區(qū)，西北部為中應(yīng)力區(qū)(圖3)；5#煤層的高應(yīng)力區(qū)分布在北部，南部為地應(yīng)力區(qū)；11#煤層的高應(yīng)力區(qū)分布在南部，北部為低地應(yīng)力區(qū)。11#煤層深度大，斷層少，整體特征是最小水平應(yīng)力值分布比較均勻，多在22～32 MPa范圍內(nèi)，3#和5#煤層埋藏淺，斷層發(fā)育，平面上地應(yīng)力變化大。

圖3 韓城礦區(qū)3#煤層最小水平主應(yīng)力分布Fig.3 Distribution of minimum horizontal principal stress in the coalbed 3# of the Hancheng mining area

4 煤儲層滲透率約束因素

4.1 試井滲透率計算

通過試井壓力恢復(fù)曲線數(shù)可以得到儲層的有效滲透率[10-12]，其分析方法有MDH法、Horner法、Ramey圖板法及Boudet圖板法等[10-11]。本次研究采用Horner分析法求取煤儲層的滲透率。

(3)

在直角坐標系中，Pws與lg[(tP+Δt)/Δt]呈線性關(guān)系，若設(shè)其斜率為m，則滲透率為：

(4)

式中：pws為井底壓力，MPa；pi為原始地層壓力，MPa；tp為開井生產(chǎn)流動時間，h；Δt為開井生產(chǎn)流動時間，h；K為滲透率，10-3μm2；q為地面的日產(chǎn)氣量，m3；μ為氣體粘度，mPa·s；B為氣體體積系數(shù)；h為有效厚度，m；m為Horner直線段的斜率。

根據(jù)計算結(jié)果，滲透率分布區(qū)間為0.001 3×10-3～0.071 0×10-3μm2，平均為0.022 0×10-3μm2。各井層間滲透率差別比較大。3#煤層的滲透率為0.002 6×10-3～0.044 0×10-3μm2，平均為0.018 3×10-3μm2；5#煤層的滲透率為0.006 0×10-3～0.071 0×10-3μm2，平均為0.023 2×10-3μm2；11#煤層的滲透率為0.001 3×10-3～0.051 7×10-3μm2，平均為0.012 8×10-3μm2；總體情況是5#煤層的滲透率稍高，3#煤層中等，11#煤層最小。

4.2 滲透率影響因素

煤儲層滲透率的影響因素有很多，比如煤巖變質(zhì)程度、煤層厚度、煤體結(jié)構(gòu)、煤巖類型、割理、微裂隙以及地應(yīng)力等因素。其中，最直觀、影響最明顯的因素是煤巖類型、割理(包括微裂隙)以及地應(yīng)力。

4.2.1 煤巖類型對滲透率的影響

光亮煤和半亮煤脆度大，機械強度小，易破碎，內(nèi)生裂隙發(fā)育，灰分含量低，密度小，含氣量高，滲透率相對高。半暗煤光澤較弱，較堅硬，密度較大，內(nèi)生裂隙欠發(fā)育，灰分含量較高，滲透率較低。暗淡煤光澤暗淡，質(zhì)地堅硬，密度大，內(nèi)生裂隙不發(fā)育，礦物質(zhì)含量高，灰分含量高，滲透率低。

根據(jù)實際資料統(tǒng)計[13]，不同的煤巖類型，割理發(fā)育程度存在很大差別(圖4)，從光亮煤到暗淡煤，割理密度逐漸降低，光亮煤中割理密度比暗淡煤割理密度高出1倍以上，滲透率高出8倍以上。由此可見,煤巖類型對割理的分布和滲透率大小的影響是至關(guān)重要的因素。

4.2.2 割理(微裂隙)對滲透率的影響

煤巖雖然屬于雙孔隙介質(zhì)結(jié)構(gòu)，但基質(zhì)孔隙屬于納米級孔隙，對滲透率貢獻很小，幾乎可以忽略，在生產(chǎn)過程中，對滲流起貢獻作用的主要是割理和微裂隙[11]。

根據(jù)實驗數(shù)據(jù)統(tǒng)計[13](圖5)，滲透率與割理(微裂隙)具有正相關(guān)關(guān)系：

K=0.161 4lnM-0.305 4,R=0.865 2

(5)

式中：K為滲透率，10-3μm2；M為割理(微裂隙)密度×開啟率，條/(5 cm)。

圖4 韓城礦區(qū)不同煤巖類型的微裂隙和滲透率Fig.4 Various microcracks of different coal-rocks and their permeability in the Hancheng mining area

圖5 韓城礦區(qū)滲透率與割理發(fā)育程度的關(guān)系Fig.5 Relationship between permeability and cleat frequency of occurrence in the Hancheng mining area

4.2.3 地應(yīng)力對滲透率的影響

因為煤巖對應(yīng)力比較敏感，煤儲層滲透性與地應(yīng)力關(guān)系密切，一般情況下，地應(yīng)力小的地區(qū)，煤儲層滲透率高，煤層氣井產(chǎn)氣量也高；反之，煤儲層滲透率低。但煤巖類型不同，巖石力學(xué)參數(shù)不同，對地應(yīng)力的敏感程度不同，因此，研究過程中分煤巖類型統(tǒng)計滲透率與地應(yīng)力的關(guān)系相關(guān)程度會提高(圖6)，四種煤巖類型的滲透率都隨最小水平主應(yīng)力值的增大呈指數(shù)下降的趨勢，且相關(guān)性較好。

光亮煤滲透率與最小主應(yīng)力關(guān)系：

K=0.100 7e-0.093 7σh,R=0.846 4

(6)

半亮煤滲透率與最小主應(yīng)力關(guān)系：

K=0.101 8e-0.073σh,R=0.832 7

(7)

半暗煤滲透率與最小主應(yīng)力關(guān)系：

K=0.108 9e-0.148 3σh,R=0.832 8

(8)

暗淡煤滲透率與最小主應(yīng)力關(guān)系：

K=0.059 1e-0.125σh,R=0.866 1

(9)

5 多因素約束下的煤儲層地質(zhì)建模

由于煤層滲透率參數(shù)的獲得一直比較困難，因此需要從多個方面研究確定。有試井的井點其滲透率利用試井滲透率控制，沒有試井但有壓裂資料的井點根據(jù)試井計算的地應(yīng)力與滲透率關(guān)系模型，在煤巖類型約束下確定井點的滲透率，其它區(qū)域的滲透率在煤巖相和地應(yīng)力模型的約束下采用序貫高斯法確定?？偨Y(jié)為四步建模法：首先建立煤儲層空間格架模型，然后在格架模型上分別建立煤巖相(類型)模型和地應(yīng)力模型，最后在三者的約束下完成滲透率模型。

圖6 韓城礦區(qū)煤滲透率與最小水平主應(yīng)力關(guān)系Fig.6 Relationship between permeability and minimum horizontal principal stress in the Hancheng mining areaa.光亮煤滲透率與最小水平主應(yīng)力關(guān)系;b.半亮煤滲透率與最小水平主應(yīng)力的關(guān)系;c.半暗煤滲透率與最小水平主應(yīng)力的關(guān)系;d.暗淡煤滲透率與最小水平主應(yīng)力的關(guān)系

5.1 煤儲層宏觀空間分布模型建立

根據(jù)沉積相演化規(guī)律和煤厚度變化趨勢分析，煤層分布明顯具有分帶性。3#煤層大致上呈北東向條帶狀分布，受分流河道的方向控制，通過厚度變差函數(shù)的分析，優(yōu)勢方向確定為50°，變程設(shè)定為260 m，次方向變程設(shè)定為120 m。5#和11#煤層厚帶基本上都是分布在東北部，即向陸方向厚度增大，向海方向厚度減薄，等厚線基本上是平行岸線。把優(yōu)勢方向確定為150°(即平行于岸線方向)，變程設(shè)定為280 m，次方向的變程確定為100 m。煤層厚度雖然在空間上變化比較大，但以漸變?yōu)橹?，橫向上具有連續(xù)性和過渡性，因此采用序貫高斯模擬法，依據(jù)井點數(shù)據(jù)，在沉積相背景的約束下，通過變程控制完成空間分布模型。

5.2 煤巖相(類型)模型建立

前面分析可以看出，煤巖類型不同，其化學(xué)組分、顯微組分不同，基質(zhì)顆粒的大小不同，巖石力學(xué)性質(zhì)也不同，對應(yīng)力敏感程度不同，直接影響到滲透率。煤巖類型是滲透率的重要約束因素，所以識別出不同的煤巖類型，建立煤巖類型空間分布模型是建立滲透率模型的基礎(chǔ)。

5.2.1 煤巖類型的測井解釋

根據(jù)對韓城煤層氣井巖電關(guān)系的研究，從光亮煤到暗淡煤，補償密度測井值逐漸增大，補償聲波時差值逐漸減小，擴徑率逐漸降低[13]。光亮煤的補償密度測井一般小于1.50 g/cm3，補償聲波時差測井大于425 μs/m井徑擴徑率大于50%；半亮煤的密度測井為1.48～1.60 g/cm3，補償聲波時差測井是400～430 μs/m，井徑擴徑率20%～50%；半暗煤的密度測井為1.58～1.67 g/cm3，補償聲波時差測井是390～420 μs/m，井徑擴徑率是10%～20%；暗淡煤的補償密度測井為1.65～1.75 g/cm3，補償聲波時差測井為365～400 μs/m，井徑擴徑率小于10%；夾矸的補償密度測井大于1.75 g/cm3，補償聲波時差測井小于370 μs/m，幾乎不擴徑。

利用煤巖測井響應(yīng)關(guān)系建立煤巖測井解釋圖板(圖7)，依據(jù)圖板完成全區(qū)各井層煤巖類型的解釋。

5.2.2 煤巖相建模

煤巖相的分布存在隨機性，各類型之間在縱向和橫向上不存在過度性，界限通常是突變的，因此選擇了示性點過程的模擬方法。示性點過程是基于目標的隨機模擬方法[14]，根據(jù)點過程的概率產(chǎn)生模擬對象在空間上的分布，然后將模擬對象的屬性標注于網(wǎng)格節(jié)點上。從地質(zhì)統(tǒng)計學(xué)的角度看，示性點過程主要是模擬對象及其屬性在三維空間上的聯(lián)合分布[15-18]。

圖7 韓城礦區(qū)煤巖類型解釋圖板Fig.7 Interpretation chart of coal lithotypes in the Hancheng mining area

在煤儲層空間幾何框架模型上，井點由測井煤巖相圖板識別結(jié)果約束，井間根據(jù)點過程的概率產(chǎn)生模擬對象在空間的分布。由于不同層間煤巖相的變化比較大，因此不同層給出不同的煤巖相的幾何約束參數(shù)。3#煤層，厚度為0.8～4.4 m，無光亮煤分布；半亮煤呈不規(guī)則的條帶狀分布，面積為0.566 3～3.664 9 km2，長軸方向是30°～90°，長度為651.3～523 2.1 m，短軸方向是0°～315°，寬度為268.6～810.1 m；半暗煤平面上連片分布，面積為2.243 2～7.534 1 km2，長軸方向是10°～45°，長度為4 182.3～7 395.6 m，短軸方向是100°～135°，寬度為1 000～1 900 m；暗淡煤平面上呈不規(guī)則圓形分布，半徑為149.9～495.1 m，面積為0.138 0～0.647 5 km2。5#煤層，厚度是0.5～9.1 m，光亮煤平面上呈不規(guī)則的橢圓狀分布，面積是0.236 5～0.914 3 km2，長軸方向為45°左右，長度是728.3～1 610.6 m，短軸方向為315°，寬度為378.7～756.0 m；半亮平面上連片分布，面積是1.022 4～9.328 4 km2，長軸方向為25°～50°，長度是4 235.7～7 253.4 m，短軸方向為300°，寬度是411.8～806.0 m；半暗煤平面上呈不規(guī)則的條帶狀、圓狀分布，面積為0.133 6～2.375 6 km2，長軸方向為40°，長度為488.3～4 009.0 m，短軸方向為310°，寬度為400.3～668.9 m；暗淡煤平面上呈不規(guī)則的圓狀分布，半徑為325.0～510.0 m。11#煤層，煤厚度為4.0～11.6 m，無光亮煤分布；半亮平面上呈條帶狀和橢圓狀分布，面積為0.415 1～2.566 0 km2，長軸方向是90°左右，長度為792.3～2 406.4 m，短軸方向0°左右，寬度100.4～1 174.9 m；半暗煤呈不規(guī)則的片狀分布，面積1.024 5～2.338 0 km2，無方向性；暗淡煤呈條帶狀分布，面積為0.231 0～2.698 2 km2，長軸方向為36°～130°，長度為408.3～5 043.6 m，短軸方向為40°～306°，寬度為590.4～922.9 m。

最終完成了3個主力煤層的煤巖相空間模型(圖8是其中的3#煤巖模型)。

5.3 地應(yīng)力模型建立

地應(yīng)力受區(qū)域構(gòu)造背景和局部構(gòu)造控制，在空間上一般具有連續(xù)性和過渡關(guān)系。根據(jù)井點數(shù)據(jù)的分析，確定不同方向的變差函數(shù)，優(yōu)勢方向是0°，變程確定為260 m，次方向變程為150 m。在模擬方法上，選用地質(zhì)統(tǒng)計學(xué)中適用于連續(xù)變量模擬的序貫高斯模擬算法，采用井點數(shù)據(jù)約束，變程控制，井間插值的模擬過程，確保模擬的結(jié)果與區(qū)域背景一致，并能體現(xiàn)局部特征，同時要與手工分析的地應(yīng)力趨勢一致。

5.4 煤儲層屬性模型的隨機建立

根據(jù)相控建模原則，在地應(yīng)力模型約束下，確定不同煤巖相的滲透率變差函數(shù)，采用協(xié)同模擬方法建立各模擬單元中每個煤巖相的滲透率模型。

具體步驟為：①有試井的井點滲透率由試井滲透率控制；②無試井但有壓裂的井點，利用煤巖約束下的最小水平主應(yīng)力與滲透率關(guān)系式，由計算出的滲透率數(shù)值約束；③以井點滲透率數(shù)據(jù)為依據(jù)，煤巖相模型和地應(yīng)力模型為約束條件，結(jié)合煤巖割理和裂隙發(fā)育的優(yōu)勢方向，通過變差函數(shù)分析，采用序貫高斯法模擬滲透率屬性三維模型，根據(jù)長寬比以及井距，采用球形函數(shù)作為變差函數(shù)。3#煤層，主方向為115°，變程為240 m，次方向變呈為120 m，垂向變程1 m；5#煤層，主方向為120°，變程為260 m，次方向變呈為130 m，垂向變程1 m；11#煤層，主方向為130°，變程為275 m，次方向變呈為135 m，垂向變程1 m。以此完成了三個煤層的地質(zhì)模型(圖9是其中的3#煤層滲透率模型)。

5.5 模型檢驗與評價

采用礦場實際生產(chǎn)數(shù)據(jù)對模型進行可靠性驗證，模型與井的生產(chǎn)現(xiàn)狀基本吻合，高滲透率區(qū)多為中高產(chǎn)井，低滲透率區(qū)以中低產(chǎn)井為主。按照厚度對韓城礦區(qū)內(nèi)各井各煤層的滲透率數(shù)值進行加權(quán)平均計算(因為區(qū)內(nèi)多為合采井)，與穩(wěn)產(chǎn)期的平均日產(chǎn)氣量進行相關(guān)性分析，結(jié)果顯示二者具有比較好的的正相關(guān)關(guān)系(圖10)。

模型中高滲透率區(qū)域與地應(yīng)力的分布趨勢也基本吻合，同時也受煤巖類型控制。如圖9是3#煤層的滲透率模型，東南方向為高滲透率分布區(qū)域，與低地應(yīng)力區(qū)基本一致(圖3)。南部也出現(xiàn)了一個相對高的滲透率分布區(qū)，主要原因是該區(qū)以半亮煤為主，使得物性比較好。

圖10 韓城礦區(qū)加權(quán)滲透率與穩(wěn)產(chǎn)期平均日產(chǎn)氣量關(guān)系Fig.10 Relationship between weighted permeability and average daily gas production in stable stages in the Hancheng mining area

6 結(jié)論

1) 煤巖類型不同，其巖石力學(xué)特征不同，割理發(fā)育程度不同，滲透率不同。光亮煤和半亮煤密度小，脆度大，機械強度小，易破碎，內(nèi)生裂隙發(fā)育，滲透率相對高。半暗煤較堅硬，密度較大，內(nèi)生裂隙欠發(fā)育，滲透率較低。暗淡煤質(zhì)地堅硬，密度大，礦物質(zhì)含量高，內(nèi)生裂隙不發(fā)育，滲透率低。煤巖類型是滲透率的重要約束條件。

2) 根據(jù)試井計算出的煤儲層滲透率和壓裂資料確定的地應(yīng)力資料分析，在分煤巖類型條件下，滲透率受最小水平主應(yīng)力約束明顯，并且具有比較好的相關(guān)性。

3) 在煤儲層空間格架模型、煤巖相(類型)模型和地應(yīng)力模型約束下完成的滲透率模型比較好地刻畫了煤儲層的屬性特征。模型與生產(chǎn)情況吻合比較好，該建模方的法理論依據(jù)充分，實際應(yīng)用中可操作性強。

致謝：國家煤層氣工程中心和中國地質(zhì)大學(xué)(北京)能源與環(huán)境教研室提供了大部分礦區(qū)資料，在此表示感謝！衷心感謝審稿專家提出的寶貴修改意見！

參 考 文 獻

[1] 邵先杰,董新秀,湯達禎,等.韓城礦區(qū)煤層氣中低產(chǎn)井治理技術(shù)與方法[J].天然氣地球科學(xué),2014,25(3):435-443.

Shao Xianjie,Dong Xinxiu,Tang Dazhen,et al.Treatment techonlogy and method of low-to-moderate production coalbed methane wells in Hancheng mining area[J].Natural Gas Geoscience,2014,25(3):435-443.

[2] Chen Yue,Tang Dazhen,Xu Hao,et al.Pore and fracture characteristics of different rank coals in the eastern margin of the Ordos Basin,China[J].Journal of Natural Gas Science and Engineering.2015,26:1264-1277.

[3] 馬東民,張遂安,藺亞兵.煤的等溫吸附—解吸實驗及其精確擬合[J].煤炭學(xué)報,2011,36(3):477-480.

Ma Dongmin,Zhang Suian,Lin Yabing.Isothermal adsorption and desorption experiment of coal and experimental results accuracy fit ting[J].Journal of China Coal Society,2011,36(3):477-480.

[4] 欒偉娜,林玉祥,吳玉琛.沁水盆地煤巖頂板致密砂巖氣充注過程及運移機理研究[J].山東科技大學(xué)學(xué)報(自然科學(xué)版),2015,34(3):38-44.DOI:10.16452/j.cnki.sdkjzk.20150619.010.

Luan Weina,Lin yuxiang,Wu Yuchen.Filling process and migration mechanism of tight sandstone gas on the coal rock roof in qinshui basin[J].Journal of Shandong University of Science and Technology(Natural Science),2015,34(3):38-44.DOI:10.16452/j.cnki.sdkjzk.20150619.010.

[5] 龍勝祥,李辛子,葉麗琴,等.國內(nèi)外煤層氣地質(zhì)對比及其啟示[J].石油與天然氣地質(zhì),2014,35(5):696-703.

Long Shengxiang,Li Xinzi,Ye Liqin,et al.Comparison and enlighte-nment of coalbed methane geology at home and abroad[J].Oil & Gas Geology,2014,35(5) :696-703.

[6] 張玉平.韓城礦區(qū)煤層氣資源評價與開發(fā)分析[J].煤炭技術(shù),2007,26(4):100-102.

Zhang Yuping.Evaluation and exploration analysis ofcoal seam gas in Hancheng Mining Area[J].Coal Technology,2007,26(4):100-102.

[7] 郝蜀民,李良,張威,等.鄂爾多斯盆地北緣石炭系—二疊系大型氣田形成條件[J].石油與天然氣地質(zhì),2016,37(2):149-154.

Hao Shumin,Li Liang,Zhang Wei,et al.Forming conditions of large-scale gas fields in Permo-Carboniferous in the northern Ordos Basin[J].Oil & Gas Geology,2016,37(2):149-154.

[8] 李勇,湯達禎,許浩,等.鄂爾多斯盆地柳林地區(qū)煤儲層地應(yīng)力場特征及其對裂隙的控制作用[J].煤炭學(xué)報,2014,39(1):164-168.

Li Yong,Tang Dazhen,Xu Hao,et al.Characteristic of in-situ stress field in Liulin area,Ordos Basin and its control on coal fractures[J].Journal of China Coal Society,2014,39(S1):164-168.

[9] 趙慶.應(yīng)用壓裂資料計算地應(yīng)力的一種方法[J].河南石油,2005,19(1):43-45.

Zhao Qing.A computing method togeostress by fracturing data[J].Henan Oil,2005,19(1):43-45.

[10] 張楓,王振升,馬立軍,等.一種檢驗試井解釋有效滲透率合理性的方法[J].天然氣地球科學(xué),2010,21(3):367-370.

Zhang Feng,Wang Zhensheng,Ma Lijun,et al.A method of testing reasonability of effective permeability by well test Interpretation[J].Natural Gas Geoscience,2010,21(3):367-370.

[11] 安小平,李相方,程時清,等.不同方法獲取滲透率的對比分析[J].油氣井測試,2005,14(5):14-17.

An Xiaoping,Li Xiangfang,Cheng Shiqing,et al.Comprartive analysis for permeability acquired from different methods[J].Well Testing,2005,14(5):14-17.

[12] 劉峻杉.煤儲層滲透率影響因素探討[J].重慶師范大學(xué)學(xué)報(自然科學(xué)版),2014,31(2):100-104.

Liu Junbin.The influencing factors of coal reservoir permeability [J].Journal of Chongqing Normal University (Natural Science),2014,31(2):100-104.

[13] 湯達禎,王生維,等.煤儲層物性控制機理及有利儲層預(yù)測方法[M].北京:科學(xué)出版社,2010.

Tang Dazhen,Wang Shengwei.et al.Control mechanism of coal prope-rties and prediction technique of favorable reservoir[M].Beijing:Science Press,2010.

[14] 邵先杰,董新秀,湯達禎,等.韓城礦區(qū)煤巖類型測井解釋技術(shù)及產(chǎn)能預(yù)測方法[J].測井技術(shù),2013,37(6):671-675.

Shao Xianjie,Dong Xinxiu,Tang Dazhen,et al.Coal type definition and productivity prediction of Hancheng mining area based on logging data[J].Well Logging Technology,2013,37(6):671-675.

[15] 王家華,張團峰.油氣儲層隨機地質(zhì)建模[M].北京:石油工業(yè)出版社,2001.

Wang Jiahua,Zhang Tuanfeng.Stochasticreservoir modellng[M].Beijing:Petroleum Industry Press,2001

[16] 霍春亮,古莉,趙春明,等.基于地震、測井和地質(zhì)綜合一體化的儲層精細建模[J].石油學(xué)報,2007,28(6):66-71.

Huo Chunliang,Gu Li,Zhao Chunming,et al.Integrated reservoir geo-logical modeling based on seismic,log and geological data[J].Acta Petrolei Sinica,2007,28(6):66-71.

[17] 劉占良,朱新磊,杜支文，等.蘇里格氣田ZT1區(qū)塊儲層地質(zhì)建模[J].山東科技大學(xué)學(xué)報(自然科學(xué)版),2015,34(1):1-8.DOI:10.16452/j.cnki.sdkjzk.2015.01.001.

Liu Zhanliang,Zhu Xinlei,Du zhiwen,et al.Reservoir geological modeling of ZT1 block in sulige gas field[J].Journal of Shandong University of Science and Technology(Natural Science).2015,34(1):1-8.DOI:10.16452/j.cnki.sdkjzk.2015.01.001.

[18] 李百川,周巨標,丁玉盛,等.不同覆壓條件下儲層物性變化特征及水驅(qū)油實驗研究[J].山東科技大學(xué)學(xué)報(自然科學(xué)版),2016,35(2):22-28.DOI:10.16452/j.cnki.sdkjzk.2016.02.014.

Li Baichuan,Zhou Jubiao,Ding Yusheng,et al.Experimental study of reservoir physical property change and water flooding under different overburden pPressures[J].Journal of Shandong University of Science and Technology(Natural Science),2016,35(2):22-28.DOI:10.16452/j.cnki.sdkjzk.2016.02.014.