趙晉斌

（山西晉城無(wú)煙煤礦業(yè)集團(tuán)有限責(zé)任公司,山西 晉城 048000)

1672-5050（2017)05-0073-04

10.3919/j.cnki.issn1672-5050sxmt.2017.10.021

2017-07-13

趙晉斌（1987-),男,山西晉城人,碩士,從事煤層氣開(kāi)發(fā)研究工作。

基于趨勢(shì)面理論的沙爾湖凹陷煤層氣儲(chǔ)量評(píng)估

趙晉斌

（山西晉城無(wú)煙煤礦業(yè)集團(tuán)有限責(zé)任公司,山西 晉城 048000)

以新疆吐哈盆地為依托,通過(guò)運(yùn)用趨勢(shì)面分析理論的數(shù)學(xué)地質(zhì)方法,結(jié)合相關(guān)的軟件,對(duì)吐哈地區(qū)沙爾湖凹陷賦存的巨厚煤進(jìn)行模擬。研究表明,該地區(qū)地質(zhì)條件較為復(fù)雜,煤層厚度起伏較大,但煤炭資源可觀,煤層氣資源量豐富,具有相當(dāng)好的開(kāi)采前景。通過(guò)建立一整套適合該區(qū)域儲(chǔ)量計(jì)算的流程,為后續(xù)的勘探開(kāi)發(fā)提供了可行的依據(jù)。

煤層氣;趨勢(shì)面分析;儲(chǔ)量估計(jì)

吐哈盆地煤層氣資源儲(chǔ)量巨大,是未來(lái)吐哈油氣資源以外重要的接替資源。該地區(qū)地質(zhì)條件較為復(fù)雜,可采煤層較多,厚度和埋深起伏較大,以沙爾湖為例,煤層埋深300 m ～1 200 m,其中凹陷區(qū)為深部含煤區(qū),巨厚煤層,最大單煤層厚度217 m,煤炭和煤層氣資源量豐富[1]。本文以沙爾湖凹陷為例,以趨勢(shì)面分析為基礎(chǔ)理論,對(duì)該地區(qū)的煤層氣儲(chǔ)量進(jìn)行評(píng)估,為今后的該地區(qū)的研究提供借鑒。

1 沙爾湖凹陷煤層氣成因

沙爾湖凹陷位于新疆吐哈盆地南緣中段,在構(gòu)造位置上處于沙爾湖背斜和覺(jué)羅塔格背斜之間,北鄰沙爾湖隆起帶,南以F1斷層與覺(jué)羅塔格復(fù)背斜相鄰,區(qū)域構(gòu)造綱要圖，見(jiàn)圖1。

圖1 區(qū)域構(gòu)造綱要圖Fig.1 Regional tectonics outline

沙爾湖凹陷的基底地層為二疊系火山巖和陸相碎屑巖,呈平行或微角度不整合于二疊系火山巖上的中侏羅統(tǒng)西山窯組是沙爾湖煤田的含煤地層,其總體特征為煤層層數(shù)多,厚度變化大。西山窯組的巖性為河湖相-泥炭沼澤相沉積的泥巖、砂巖、礫巖和煤層[2]。

根據(jù)煤層生氣的階段性,在Ro<0.4%時(shí),煤層以生成生物氣為主,Ro>0.4%時(shí),煤層進(jìn)入熱解生氣階段。沙爾湖凹陷煤階處于低變質(zhì)褐煤-長(zhǎng)焰煤階段,Ro介于0.31%～0.7%,煤層氣具有熱演化和生物氣兩種成因。

沙爾湖凹陷主煤層具有罕見(jiàn)的低灰-特低灰特點(diǎn),灰分含量只有3%～8%,非常有利于有機(jī)組分對(duì)甲烷的吸附,也有利于煤層割理及裂縫的發(fā)育,同時(shí),灰分含量低,意味著煤層孔隙很少被泥質(zhì)充填,煤層孔隙發(fā)育,滲透性好。值得指出的是該區(qū)煤層的相對(duì)密度明顯偏低,一般只有1.28～1.36,原因除了煤中礦物質(zhì)含量低以外,地層壓實(shí)程度、煤的演化程度也很低,煤基質(zhì)層間縫、粒間孔隙、植物組織孔、氣孔等孔隙都非常發(fā)育。

凹陷北部沙試1井鉆探證實(shí),西山窯組主煤層具有較好的含氣性。在500 m ～600 m 井段,原煤含氣量1.5 m3/t ～3.09 m3/t,含氣飽和度80.5%～91.0%,對(duì)于低煤階褐煤而言,煤層含氣量及含氣飽和度都比較高。實(shí)驗(yàn)數(shù)據(jù)表明隨著埋深增加含氣量、含氣飽和度有明顯增加。

2 趨勢(shì)面方程的原理與生成

趨勢(shì)面分析,是一種基于數(shù)學(xué)來(lái)模擬曲面的研究分析,是元素在空間上展布所反映研究對(duì)象變化趨勢(shì)的一種數(shù)學(xué)方法。作為回歸分析的一種重要類(lèi)型,本質(zhì)上就是運(yùn)用最小二乘法來(lái)擬合的一個(gè)二維乃至多維的非線性函數(shù)[3]。

假設(shè)某煤礦含煤區(qū)域的鉆孔實(shí)際觀測(cè)數(shù)據(jù)為:（Xi,Yi,Zi)（i=1,2,…,n).

則有趨勢(shì)擬合值（xi,yi,zi),且有如下關(guān)系:

（Xi,Yi,Zi)=（xi,yi,zi)+ei.

其中,ei為剩余值。顯然,當(dāng)元素（Xi,Yi,Zi)在空間上變動(dòng)時(shí),上式就形象的描述了元素的實(shí)際布分曲面、趨勢(shì)面以及剩余面之間的互動(dòng)關(guān)系。趨勢(shì)面的核心就是從實(shí)際觀測(cè)值出發(fā),利用回歸分析,使得剩余值平方和區(qū)域最小,即最小二乘法下的趨勢(shì)面擬合。表達(dá)式如下:

用來(lái)計(jì)算趨勢(shì)面的數(shù)學(xué)方程一般有多項(xiàng)式函數(shù)和傅立葉級(jí)數(shù),其中最常用的就是以多項(xiàng)式函數(shù)的形式計(jì)算的[4]。一次、二次、三次趨勢(shì)面方程如下:

z1=a0+a1x+a2y,

z2=b0+b1x+b2y+b3x2+b4xy+b5y2,

z3=c0+c1x+c2y+c3x2+c4xy+c5y2

本報(bào)訊 11月13日，在經(jīng)濟(jì)活動(dòng)分析會(huì)上，湖北化肥公司財(cái)務(wù)相關(guān)負(fù)責(zé)人用圖表、趨勢(shì)圖對(duì)標(biāo)分析，深入淺出介紹了1-10月物資庫(kù)存情況。

+c6x3+c7x2y+c8xy2+c9y3.

綜上,趨勢(shì)面方程所求的實(shí)質(zhì)就是根據(jù)觀測(cè)值Xi,Yi,Zi（i=1,2,…,n),來(lái)確定多項(xiàng)式的系數(shù)a0、a1、a2、b0,…,b5、c0,…,c9,…,使剩余值平方和最小。

2.1趨勢(shì)面方程建立

沙爾湖區(qū)域煤厚起伏較大,范圍較大,為了能最好將所有區(qū)域覆蓋到,選取了包括沙1、2、3、4鉆孔在內(nèi)的94組觀測(cè)點(diǎn)。其中,鉆孔坐標(biāo)（x,y)作為自變量,z（煤厚)作為應(yīng)變量,構(gòu)成了趨勢(shì)面分析的基本數(shù)據(jù)。一般來(lái)說(shuō),回歸分析次數(shù)越高,計(jì)算越復(fù)雜,并且擬合關(guān)鍵的數(shù)據(jù)點(diǎn)時(shí),會(huì)造成很大的誤差。所以一般每次的數(shù)據(jù)變量只做1到3次的趨勢(shì)面分析。

由于趨勢(shì)面分析設(shè)計(jì)到十分復(fù)雜的矩陣運(yùn)算,為了便于編譯計(jì)算,選用了MATLAB作為程序的編譯工具。本文根據(jù)趨勢(shì)面方程回歸時(shí)的擬合程度的局限性,對(duì)趨勢(shì)方程做到三次的擬合運(yùn)算。通過(guò)程序編譯后載入后輸出分為兩部分,分別為一次二次三次趨勢(shì)方程擬合后與實(shí)際z值（煤厚)的剩余值（和各趨勢(shì)方程函數(shù)的系數(shù)值,即a0、a1、a2、b0，…b5、c0，…，c9的值。趨勢(shì)方程整理后如下:

z1=-9.588 8×103+0.000 6x+0.009 9y,

7.001 8×10-6xy-9.747 6×10-6y2,

z3=-7.657×10-7x2+1.625 6×10-6xy-

4.306 6×10-10x3+2.780 9×10-1x2y-

5.981 8×10-9xy2+4.286 3×10-9y3.

2.2趨勢(shì)面方程輸出及該地區(qū)沉積特點(diǎn)

根據(jù)MATLAB所求的方程,重新擬合數(shù)據(jù)點(diǎn),結(jié)合繪圖軟件SURFER,繪制出該區(qū)域一次、二次、三次趨勢(shì)面圖及對(duì)應(yīng)的剩余圖（偏差圖)如圖2-圖7。其中三次趨勢(shì)方程的F值遠(yuǎn)遠(yuǎn)超過(guò)F檢驗(yàn)的臨界值F（a),趨勢(shì)面方程顯著,即該模型從整體上是可用的。從趨勢(shì)圖、偏差圖、等值線圖（圖8)以及該地區(qū)的地質(zhì)情況可以清晰的觀測(cè)到沙爾湖地區(qū)的煤層空間分布,在該區(qū)域南部以及斷層西側(cè)的觀測(cè)點(diǎn)有較大的剩余值,因而其厚度較大,而周邊范圍剩余值較小,煤厚變化幅度小。反應(yīng)在沉積當(dāng)時(shí),沉積基底沉降速度由四周向斷層區(qū)域和沙爾湖東部中心區(qū)逐漸變大,沉積物的堆積速度逐漸變快。

圖2 沙爾湖含煤盆地一次趨勢(shì)圖Fig.2 Simple trend of coal-bearing basin in Sha’erhu

圖3 沙爾湖含煤盆地一次偏差圖Fig.3 Simple deviation of coal-bearing basin in Sha’erhu

圖4 沙爾湖含煤盆地二次趨勢(shì)圖Fig.4 Quadratic trend of coal-bearing basin in Sha’erhu

圖5 沙爾湖含煤盆地二次偏差圖Fig.5 Quadratic deviation of coal-bearing basin in Sha’erhu

圖6 沙爾湖含煤盆地三次趨勢(shì)圖Fig.6 Cubic trend of coal-bearing basin in Sha’erhu

圖7 沙爾湖含煤盆地三次偏差圖Fig.7 Cubic deviation of coal-bearing basin in Sha’erhu

圖8 沙爾湖煤厚等值線圖Fig.8 Contours of coal thickness in Sha’erhu

3 煤層氣儲(chǔ)量評(píng)估與結(jié)論

根據(jù)按照煤層厚度、含氣量、煤層埋深劃分不同等級(jí)的區(qū)塊,進(jìn)而進(jìn)行區(qū)塊劃分,按照煤層厚度、含氣量加權(quán)平均計(jì)算儲(chǔ)量。具體計(jì)算公式如下:

Gi=AihiDiCi=ViDiCi（i=1,2,…,n).

式中:Gi為第i塊煤層氣儲(chǔ)量,m3;Ai為第i塊含氣面積,km2;hi為第i塊煤層加權(quán)平均厚度,m;Vi為第i塊含煤盆地體積,m3;Di為第i塊煤層的平均密度,t/m3;Ci為第i塊煤層平均含氣量,t/m3。則:

G=G1+G2+G3+…+Gn,即G為煤層氣田的總儲(chǔ)量。

通過(guò)第2部分的驗(yàn)證,選取的三次趨勢(shì)方程是最能擬合沙爾湖地區(qū)煤層展布情況,通過(guò)surfer軟件的輸出可以得到以下結(jié)果和認(rèn)識(shí):

1)由z的最小值和最大值,可以得出模擬后的沙爾湖的煤厚范圍,其中最厚煤層（沙爾湖凹陷中心)模擬后的厚度值為z最大值與最小值的差值,即H=236.66-（-79.07)=315.73,與沙爾湖探測(cè)中部聚煤中心累計(jì)最厚煤層309.27m相接近。

2)根據(jù)surfer中的梯形規(guī)則,得出:V約為613 m3,即為該地區(qū)煤層展布區(qū)域體積。由該地區(qū)煤巖密度r=1.3 t/m3,可估算該區(qū)域的煤炭?jī)?chǔ)約為796.75億t。結(jié)合生氣成因的因素,以及對(duì)各部分含氣面積的計(jì)算、含氣量參數(shù)的確定,得出沙爾湖地區(qū)含氣量約為1 464.47×108m3。

[1] 楊珍祥,李巧梅.吐哈盆地沙爾湖洼陷煤層氣成藏地質(zhì)特征[J].新疆石油地質(zhì),2008,29（6):713-715.

YANG Zhenxiang,LI Qiaomei.Geologic Characteristics of Coal Seam Gas Accumulation in Shaerhu Sub-Sag,Tuha Basin[J].Xinjiang Petroleum Geology,2008,29（6):713-715.

[2] 吳濤,趙文智.吐哈盆地煤系油氣田形成和分布[M].北京:石油工業(yè)出版社,1997.

[3] 李哲,張淑英.基于多項(xiàng)式趨勢(shì)面分析理論的天然氣需求預(yù)測(cè)[J].資源與產(chǎn)業(yè),2008,10（2):105-107.

LI Zhe,ZHANG Shuying.Prediction of Natural Gas Demamds Based on Muliti-trendanalysis[J].Resources & Industries,2008,10（2):105-107.

[4] 康永尚,沈金松,諶卓恒.現(xiàn)代數(shù)學(xué)地質(zhì)[M].北京:石油工業(yè)出版社,2005.

ReservesEstimationofCoalbedMethaneinSha’erhuSagBasedonTrendSurfaceTheory

ZHAOJinbin

（ShanxiJinchengAnthraciteMiningGroupCo.,Ltd.,Jincheng048000,China)

Based on Turpan-Hami Basin, Xinjiang, the mathematical and geological methods of trend surface theory are used to simulate the thick coal in Sha’erhu Sag. In spite of the complexity of the geological situation with large terrain variation, the results show huge reserves and rich coalbed methane resources with a promising future. By establishing a series of processes to calculate the reserves for the area, the study provides a feasible basis for the future exploration.

coalbed methane; trend surface analysis; reserves estimation

P618

A

（編輯:樊 敏)