程　昊

(西北大學(xué)地質(zhì)學(xué)系/大陸動力學(xué)國家重點實驗室，陜西 西安710069)

?

基于測井資料的陸相頁巖層系巖性識別方法

程昊

(西北大學(xué)地質(zhì)學(xué)系/大陸動力學(xué)國家重點實驗室，陜西 西安710069)

[摘要]陸相頁巖層系巖性復(fù)雜，使用常規(guī)的測井巖性識別方法難以有效區(qū)分。因此尋找一種有效的針對不同巖性的識別方法尤為重要。通過分析陸相頁巖層系不同巖性類型測井響應(yīng)特征，統(tǒng)計出不同巖性的測井響應(yīng)分布區(qū)間，研究不同巖性測井響應(yīng)的差異，在明確不同巖性響應(yīng)特征基礎(chǔ)上，采用因子分析法，降低數(shù)據(jù)分析維度，提取對巖性敏感的測井參數(shù)，使用測井反演的方法矯正薄層及圍巖效應(yīng)對測井真實值的影響，運用K-Means聚類分析法對頁巖層系的巖性進行識別。與實際巖心對比表明，該方法能夠有效識別頁巖層系復(fù)雜巖性類型，準確性較好。

[關(guān)鍵詞]陸相；頁巖層系；巖性；測井識別

以美國為代表的北美頁巖氣勘探區(qū)主要為海相地層，沉積環(huán)境穩(wěn)定，巖性相對均一，但從已公開發(fā)表的文獻來看；其中也不乏的砂質(zhì)夾層、砂質(zhì)紋層的發(fā)育；而研究者又恰恰沒有重視這方面的現(xiàn)象，而把過多的精力集中于巖相的劃分[1-5]。鄂爾多斯盆地三疊系延長組長7段陸相沉積，沉積環(huán)境多變，頁巖層系發(fā)育有多種巖性，具有很強的非均質(zhì)性，除了存在有頁巖之外，還發(fā)育著粉砂巖、泥質(zhì)粉砂巖、砂質(zhì)頁巖。這些不同的巖性，厚度不等，巖性變化快，不同程度地增加了儲層巖性的復(fù)雜程度，使得陸相頁巖層系的巖性遠復(fù)雜于海相[6-11]。巖性不同，意味著其礦物成分、有機質(zhì)含量、脆性指數(shù)等有關(guān)參數(shù)的差異，這些參數(shù)對后期頁巖氣有利區(qū)預(yù)測及資源量的計算都有至關(guān)重要的作用。因此準確識別巖性，明確頁巖層系內(nèi)不同巖性的空間展布，對于頁巖氣勘探意義重大。

1陸相頁巖層系的測井響應(yīng)特征

1.1自然伽馬測井響應(yīng)

巖石的放射性主要是由于泥質(zhì)吸附放射性物質(zhì)而產(chǎn)生，還受有機質(zhì)及放射性物質(zhì)的影響。頁巖除了含有泥質(zhì)吸附產(chǎn)生放射性外，由于其含有大量的有機質(zhì)，有機質(zhì)對放射性物質(zhì)的吸附使得放射性進一步增大。因此，頁巖通常具有較高的自然伽馬，而粉砂巖、泥質(zhì)粉砂巖、砂質(zhì)頁巖由于泥質(zhì)含量的減少自然伽馬放射性逐漸降低，其中頁巖的自然伽馬值絕大多數(shù)在140 ApI以上，其它巖性相比于頁巖均有所降低，粉砂巖自然伽馬值最低，砂質(zhì)頁巖、泥質(zhì)粉砂巖的響應(yīng)值介于上述巖性之間(圖1A)。

1.2聲波測井響應(yīng)特征

在巖石骨架一致的情況下，巖石的聲波時差主要受到流體成分、有機質(zhì)含量等影響。頁巖中有機質(zhì)含量高、又含有天然氣，使得頁巖聲波測井表現(xiàn)出高值。相比于頁巖，粉砂巖、泥質(zhì)粉砂巖、砂質(zhì)頁巖則主要受礦物成分和巖石本身的致密程度表現(xiàn)出不同的聲波時差值。在響應(yīng)值上，頁巖的聲波時差響應(yīng)值主要在240～280 us/m之間，遠高于其它巖性，砂質(zhì)頁巖、泥質(zhì)粉砂巖、粉砂巖則呈現(xiàn)出依次降低的趨勢(圖1B)。

1.3中子測井響應(yīng)

富含有機質(zhì)的地層中的含氫量高，因此對快中子減速能力強，因此在頁巖層系內(nèi)頁巖顯示為高中子孔隙度。受地層巖性、流體性質(zhì)影響較大發(fā)育在頁巖中的砂質(zhì)頁巖繼承了頁巖的中子響應(yīng)特征；粉砂巖、泥質(zhì)粉砂巖中子測井響應(yīng)值整體上較小，由巖性由泥質(zhì)粉砂巖向粉砂巖變化，中子響應(yīng)值并無明顯的變化，從其中子測井響應(yīng)值上可以看出，兩者幾乎位于同一數(shù)值變化范圍內(nèi)(圖1C)。

1.4密度測井響應(yīng)特征

由于其礦物組成、流體含量有機質(zhì)等差異，不同巖性在密度測井曲線上表現(xiàn)為不同的響應(yīng)特征。頁巖中富含有機質(zhì)，其密度低于其他巖層，在密度測井曲線上顯示為低異常；而粉砂質(zhì)頁巖、泥質(zhì)粉砂巖、粉砂巖等石英長石含量較多，密度表現(xiàn)出較高值的特征。從統(tǒng)計的不同巖性密度測井響應(yīng)值分布區(qū)間來看頁巖密度較低，粉砂巖密度最大，泥質(zhì)粉沙巖、砂質(zhì)頁巖的密度響應(yīng)值介于兩者之間，密度響應(yīng)值存在著較大的重疊(圖1D)。

1.5電阻率測井響應(yīng)

電阻率測井曲線的變化是多種因素綜合的結(jié)果，而且由于其所含的礦物成分、孔隙大小、流體性質(zhì)以及含油氣多少的差異，電阻率表現(xiàn)出不同的變化特征。頁巖有機質(zhì)含量高，而有機質(zhì)導(dǎo)電性差，因此頁巖的電阻率呈現(xiàn)高值的特征。相對于頁巖，粉砂質(zhì)頁巖、泥質(zhì)粉砂巖等電阻率整體較小。電阻率測井曲線能夠分辨出頁巖外，對于砂質(zhì)頁巖、泥質(zhì)粉砂巖以及粉砂巖，在電阻率曲線上幾乎無法區(qū)分(圖1E)。

1.6光電吸收截面指數(shù)

光電吸收界面指數(shù)與構(gòu)成地層的礦物有著密切的關(guān)系。石英、長石等PE值較小，伊利石、綠泥石等粘土礦物具有較高的PE值。頁巖層系的主要礦物成分為石英長石以及粘土礦物；由于頁巖中粘土礦物占有主要優(yōu)勢，因此PE值為高值，而砂質(zhì)頁巖略微有所降低；泥質(zhì)粉砂巖、粉砂巖由于粘土礦物含量減小以及石英長石等的增加，PE值逐相對較小。光電吸收截面指數(shù)值分布特征表明，不同巖性的光電吸收截面指數(shù)分布區(qū)間重疊程度較大(圖1F)。

圖1　長7段頁巖層系不同巖性測井響應(yīng)值統(tǒng)計

2頁巖層系巖性識別方法

由頁巖到粉砂巖，隨著巖性的演變，不同測井變量并沒有呈現(xiàn)出與之對應(yīng)線性變化關(guān)系。因此，簡單的使用常規(guī)方法(交會圖、重疊圖)進行巖性幾乎不可能同時區(qū)分出如此復(fù)雜巖性[12-13]。而且受巖性厚度變化以及圍巖測井響應(yīng)值的影響，薄層的測井響應(yīng)值并不是其真實值。為此，基于巖石測井響應(yīng)為基礎(chǔ)，采用因子分析方法，降低數(shù)據(jù)分析維度，優(yōu)選出對巖性敏感的測井參數(shù)，應(yīng)用非線性測井反演方法，矯正圍巖效應(yīng)的影響，使用K-means聚類分析的方法，進而實現(xiàn)對巖性的快速準確識別。

非線性反演是美國學(xué)者Heidari等2011年提出并成功用于地層評價的一種有效的方法。該方法通過分析地層中的骨架成分的組成，采用數(shù)值模擬的方法，分別計算地層在不同厚度條件下的測井響應(yīng)值，分析不同厚度圍巖、薄層之間的相互影響，通過模擬不同條件下實測值與模擬值之間的差異，以矯正由于圍巖效應(yīng)對薄層的影響。

因子分析假設(shè)一組多元數(shù)據(jù)與控制因素的關(guān)系，通過量化各個測井參數(shù)之間的協(xié)方差，分析可能被同一個因素影響的參數(shù)，給出不同參數(shù)之間的相關(guān)關(guān)系。受同一因素影響的變量，在因子分析變量圖中，分布在相近或者相同的區(qū)域。通過預(yù)測各測井參數(shù)重疊性，找出這些參數(shù)相互之間的重疊部分，分析這些重疊部分所反映巖石的共同特征，達到減少分析數(shù)據(jù)的維度及數(shù)量。

聚類分析是基于統(tǒng)計理念數(shù)據(jù)集的分類方法，尤其是K-means聚類分析法廣泛應(yīng)用于具有相似屬性特征分類，該方法旨在定量化獲得大量數(shù)據(jù)點N維空間特征，通過探索大量數(shù)據(jù)之間的結(jié)構(gòu)從而避免了研究多維變量。在實際應(yīng)用具體類型的劃分類應(yīng)基于研究區(qū)目的層段地質(zhì)背景。主要原理是基于具有相同或者相似屬性的樣品會落入同一片區(qū)域，然后根據(jù)樣品點之間的距離，最終劃分出不同的類型。

采用因子分析方法，優(yōu)選出巖性識別的敏感曲線，應(yīng)用測井反演的方法，矯正薄層及圍巖效應(yīng)對測井曲線真實值的影響，使用K-Means聚類分析對B井進行了巖性識別。從變量在因子分析結(jié)果圖(圖2)的分布區(qū)域來看，主要分為三組。自然伽馬和電阻率反映地層有機質(zhì)、泥質(zhì)含量；聲波時差、補償中子反映巖石孔隙空間；光電吸收截面指數(shù)、密度反映巖石骨架成分。根據(jù)變量對巖性的分辨程度(見圖1)，最終挑選出聲波時差、自然伽馬、密度作為巖性識別的參數(shù)。

圖2　A井因子分析結(jié)果

巖心觀察可知，該層段主要發(fā)育有四種巖性，因此，聚類分析中的類型設(shè)置為4類。圖3為應(yīng)用因子分析法優(yōu)選出的對巖性區(qū)分敏感的測井參數(shù)，經(jīng)過測井反演之后，使用K-Means聚類分析法識別的不同巖性的測井響應(yīng)值分布區(qū)間。其中，黃色樣品點具有非常高的自然伽馬、聲波時差以及低密度的特征，與頁巖測井響應(yīng)特征對應(yīng)；藍色樣品點的自然伽馬、聲波時差略小于黃色樣品點，巖性類型為砂質(zhì)頁巖；紅色樣品點具有非常高的密度值，聲波時差、自然伽馬明顯低于其他類型，為粉砂巖；綠色樣品點則為泥質(zhì)粉砂巖。將上述識別結(jié)果與實測巖心進行對照分析表明(圖4)，采用該方法識別出的巖性類型與實際巖心巖性信息較為一致。

圖3　A井聚類分析后各類型的聲波、自然伽馬、密度交會圖

圖4　A井巖性識別結(jié)果與巖心對比圖

3結(jié)語

(1)陸相頁巖層系巖性類型復(fù)雜，發(fā)育有頁巖、砂質(zhì)頁巖、泥質(zhì)粉砂巖以及粉砂巖，頁巖表現(xiàn)出高自然伽馬、高聲波時差、高電阻率、光電吸收截面指數(shù)、高中子以及低密度測井響應(yīng)特征，砂質(zhì)頁巖、泥質(zhì)粉砂巖、粉砂巖的測井響應(yīng)幅度相對較小，不同巖性測井響應(yīng)特征存在著較大的重疊。

(2)采用因子分析法，優(yōu)選出對巖性敏感的測井參數(shù)，運用測井反演的方法矯正薄層及圍巖效應(yīng)對地層真實響應(yīng)值的影響，使用K-Means聚類分析法，對頁巖層系巖性進行了識別，與實際巖心對比表明，該方法在頁巖層系巖性識別中效果較好。

參考文獻

[1]Loucks R G，Stephen C. Ruppel.Mississippian Barnett Shale： Lithofacies and depositional setting of a deep-water shale-gas succession in the Fort Worth Basin, Texas[J]. AAPG Bulletin，2007， 91(4 )：579-601.

[2]Ursula Hammes，Hamlin H S， Thomas E E. Geologic analysis of the Upper Jurassic Haynesville Shale in east Texas and west Louisiana[J]. AAPG Bulletin，2011, 95(10 )： 1643-1666.

[3]Guochang Wang， Timothy R. Carr. Organic-rich Marcellus Shale lithofacies modeling and distribution pattern analysis in the Appalachian Basin[J]. AAPG Bulletin，2013，97(4 )：2173-2205.

[4]James J，Hickey，Bo Henk. Lithofacies summary of the Mississippian Barnett Shale, Mitchell 2 T.P. Sims well,Wise County, Texas[J]. AAPG Bulletin，2007，91(4 )：437-443.

[5]Mohamed O. Abouelresh and Roger M. Slatt.Lithofacies and sequence stratigraphy of the Barnett Shale in east-central Fort Worth Basin, Texas. AAPG Bulletin[J].2012, 96(1)：1-22.

[6]徐士林，包書景. 鄂爾多斯盆地三疊系延長組頁巖氣形成條件及有利發(fā)育區(qū)預(yù)測[J]. 天然氣地球科學(xué).2009，20(3)：460-465.

[7]王社教，李登華，李建忠，等. 鄂爾多斯盆地頁巖氣勘探潛力分析[J]. 天然氣工業(yè).2011，31(12)：40-46.

[8]董大忠，鄒才能，楊樺，等. 中國頁巖氣勘探開發(fā)進展與發(fā)展前景[J]. 石油學(xué)報.2012，33(1)：107-114.

[9]姜呈馥，王香增，張麗霞，等. 鄂爾多斯盆地東南部延長組長7段陸相頁巖氣地質(zhì)特征及勘探潛力評價[J]. 中國地質(zhì).2013，40(6)：1880-1888.

[10]王香增，高勝利，高潮. 鄂爾多斯盆地南部中生界陸相頁巖氣地質(zhì)特征[J]. 石油勘探與開發(fā).2014，41(3)：294-304.

[11]張燁毓，周文，唐瑜，等. 鄂爾多斯盆地三疊系長7油層組頁巖儲層特征[J].成都理工大學(xué)學(xué)報(自然科學(xué)版).2013，40(6)：671-676.

[12]嚴偉，王建波，劉帥，等. 四川盆地焦石壩地區(qū)龍馬溪組泥頁巖儲層測井識別[J]. 天然氣工業(yè).2014，06：30-36.

[13]齊寶權(quán)，楊小兵，張樹東，等. 應(yīng)用測井資料評價四川盆地南部頁巖氣儲層[J]. 天然氣工業(yè).2011，31(4)：44-47.

[中圖分類號]P619.22+7

[文獻標識碼]B

[文章編號]1004-1184(2016)02-0111-03

[作者簡介]程昊(1991-)，男，陜西商州人，在讀碩士研究生，主攻方向：測井資料處理與解釋。

[收稿日期]2015-03-14