王欣，薛林福，遇運良，白燁，石玉江，張海濤

（1.吉林大學(xué)地球科學(xué)學(xué)院，吉林長春130062；2.中國石油長慶油田公司勘探開發(fā)研究院，陜西西安710018）

壓實法石英砂巖成巖相識別及其應(yīng)用

王欣1，薛林福1，遇運良1，白燁1，石玉江2，張海濤2

（1.吉林大學(xué)地球科學(xué)學(xué)院，吉林長春130062；2.中國石油長慶油田公司勘探開發(fā)研究院，陜西西安710018）

提出利用測井資料進行石英砂巖成巖相有效劃分的壓實法。在地層埋藏史恢復(fù)的基礎(chǔ)上，通過回歸分析方法確定僅有機械壓實作用情況下的理論孔隙度－深度曲線與實測孔隙度－深度曲線進行對比，配合顯微鏡下鑄體薄片觀察及掃描電鏡觀察，X－衍射和陰極發(fā)光分析，確定石英砂巖成巖相劃分標(biāo)準(zhǔn)。根據(jù)成巖相劃分標(biāo)準(zhǔn)，利用測井資料進行實際孔隙度計算，并進行流體校正，對未取心井段石英砂巖成巖作用類型進行識別。該方法用于蘇里格地區(qū)石英砂巖成巖相劃分，把蘇里格地區(qū)石英砂巖劃分為強、中、弱溶蝕相和強、中、弱膠結(jié)相等6類成巖相類型。研究結(jié)果表明，該方法提高了石英砂巖測井成巖相識別能力和定量評價精度，具有較好的應(yīng)用前景。

測井解釋；石英砂巖；壓實法；成巖相；成巖作用

0 引 言

隨著勘探開發(fā)程度的提高，成巖作用研究在儲層描述、油藏分析和后期的勘探開發(fā)等方面的重要性日益突出［1］。成巖相是在成巖與構(gòu)造等作用下沉積物經(jīng)歷一定成巖作用和演化階段的產(chǎn)物［2］。國內(nèi)外關(guān)于碳酸鹽成巖相的研究相對較多，對石英砂巖成巖相的研究相對較少，而石英砂巖作為主要的砂巖儲層，它的成巖相越來越受到研究者的關(guān)注。常規(guī)的分析方法如鑄體薄片、掃描電鏡、X－射線衍射、流體包裹體等成本高，只能進行定性判別，不能滿足儲層評價的需要。作為對地層特征響應(yīng)的測井資料包含了豐富地質(zhì)信息，具有全井段采集的優(yōu)點。有效挖掘測井信息，探討測井資料用于進行石英砂巖成巖相類型的定量判別具有可行性與重要性。研究作為對分析及化驗資料的擴充和補充，可以滿足勘探開發(fā)的要求。

1 壓實法成巖相判別

石英砂巖的成巖作用主要為機械壓實作用、溶蝕作用、膠結(jié)作用。機械壓實構(gòu)成了孔隙度－深度關(guān)系基準(zhǔn)，溶蝕和膠結(jié)作用造成了巖石在孔隙度－深度基準(zhǔn)線左右的變化。要確定石英砂巖發(fā)生了溶蝕作用還是膠結(jié)作用，需要建立一個只考慮機械壓實作用的孔隙度－深度關(guān)系的理論曲線作為基準(zhǔn)，與實際測井獲得的孔隙度－深度曲線對比，配合巖心薄片對比觀察，確定石英砂巖成巖相劃分標(biāo)準(zhǔn)，在此基礎(chǔ)上進行成巖相定量的類型劃分。壓實法是通過恢復(fù)僅有壓實作用情況下的理論壓實曲線與現(xiàn)今孔隙度進行比較來確定成巖作用的一種方法。如果石英砂巖在埋藏壓實作用基礎(chǔ)上經(jīng)歷了膠結(jié)，則巖石的現(xiàn)今孔隙度小于理論孔隙度，如果現(xiàn)今巖石孔隙度大于理論孔隙度，則巖石可能經(jīng)歷了溶蝕過程。

建立孔隙度－深度關(guān)系理論方程

式中，φ0為初始孔隙度；c為壓實系數(shù)，反映孔隙度隨不同深度的變化規(guī)律；z為砂巖埋藏深度；φ（z）為沉積物在深度z處的孔隙度。

壓實法成巖相判別步驟見圖1。

圖1 壓實法成巖相判別流程圖

2 壓實法在蘇里格地區(qū)的應(yīng)用

2.1 蘇里格地區(qū)地質(zhì)概況

蘇里格地區(qū)位于鄂爾多斯盆地伊陜斜坡西部，勘探范圍西起內(nèi)蒙古自治區(qū)鄂托克前旗，東到陜西省靖邊縣，北起內(nèi)蒙古自治區(qū)鄂托克后旗的敖包加汗，南至陜西省吳旗縣，勘探范圍［3］約4×104m2（見圖2），是發(fā)育于上古生界碎屑巖系中的大型砂巖巖性圈閉氣藏［4］。上古生界自下而上可劃分為石炭系本溪組、二疊系太原組、山西組、下石盒子組、上石盒子組和石千峰組，其中，上、下石盒子組自上而下可劃分為盒1至盒8等8個含氣層段［5］；山西組劃分為山1、山2段。本文研究的主要為二疊系下石盒子組盒8段和山西組山1段的石英砂巖。

圖2 蘇里格地區(qū)構(gòu)造區(qū)劃圖

2.2 蘇里格地區(qū)石英砂巖成巖作用類型及對孔隙的影響

通過顯微鏡下鑄體薄片觀察及掃描電鏡觀察，結(jié)合X－衍射分析、陰極發(fā)光分析，對蘇里格地區(qū)上古生界盒8段及山1段地層石英砂巖的成巖作用類型及其表現(xiàn)、成巖作用對石英砂巖的影響進行了系統(tǒng)研究。結(jié)果表明，研究層段石英砂巖主要發(fā)育機械壓實、溶蝕和膠結(jié)3種成巖作用。

2.2.1 機械壓實作用

蘇里格地區(qū)壓實作用主要表現(xiàn)為儲層中原生粒間孔隙的縮小與減少［6－7］。隨著沉積物埋藏深度增大，壓實作用強度逐漸加大，雜基含量逐漸減少，膠結(jié)類型呈現(xiàn)出由基底式→孔隙式→接觸式→無膠結(jié)物式的變化趨勢；相應(yīng)地顆粒接觸方式由飄浮狀→點接觸→線接觸→凹凸接觸逐漸向縫合線接觸演化。

2.2.2 溶蝕作用

石英顆粒不易溶解，根據(jù)薄片觀察資料，蘇里格地區(qū)石英砂巖溶解作用主要表現(xiàn)為易溶填充物、凝灰質(zhì)、巖屑和雜基的溶蝕。在溶蝕作用強烈情況下可發(fā)生少數(shù)石英巖巖屑和碎屑的溶蝕。溶蝕作用對孔隙度的演化具有建設(shè)性的作用［8－10］。

2.2.3 膠結(jié)作用

蘇里格地區(qū)石英砂巖主要的膠結(jié)作用類型為硅質(zhì)膠結(jié)、方解石和黏土充填膠結(jié)等作用，硅質(zhì)膠結(jié)主要以次生石英和粒間充填自生微晶粒狀石英2種方式產(chǎn)出。膠結(jié)類型以孔隙式膠結(jié)為主，其次為再生－孔隙式膠結(jié)。黏土礦物膠結(jié)作用主要以高嶺石膠結(jié)為主［11－12］。

2.3 考慮機械壓實作用的孔隙度－深度關(guān)系理論曲線的確定

2.3.1 聲波時差法剝蝕量恢復(fù)

石英砂巖的粒度、分選和磨圓、埋藏歷史和后期成巖作用影響著石英砂巖的理論孔隙度－深度曲線，所以要建立理論孔隙度－深度關(guān)系曲線，必須進行埋藏史的恢復(fù)。剝蝕量影響著理論孔隙度的計算，巖石的理論孔隙度應(yīng)等于最大埋深時的孔隙度。通過剝蝕量恢復(fù)，確定巖石最大埋藏深度，可以計算出巖石的理論孔隙度。

求得的剝蝕厚度若大于上覆沉積厚度，則剝蝕地層及其下伏的各地層最大埋深大于現(xiàn)有的埋深，應(yīng)將它們分別加上剝蝕厚度，再進行古厚度恢復(fù)；若剝蝕厚度小于上覆沉積厚度，則上覆地層將補償剝蝕影響，在孔隙度－深度關(guān)系上是正常壓實趨勢，此時無需作任何校正。

中生代以來，蘇里格地區(qū)發(fā)生了4期均勻抬升和地層剝蝕事件：三疊系／侏羅系（印支期，200a）、侏羅系延安組／侏羅系直羅組（燕山運動早期，175 a）、侏羅系／白堊系（燕山運動中期，145a）和白堊系／第四系（燕山運動晚期，65a）［13］。經(jīng)過前人研究，蘇里格地區(qū)前3期剝蝕量較小，不需要進行剝蝕量恢復(fù)，因此，只需對第4期白堊紀(jì)／第四系的剝蝕量進行計算。巖石聲波時差隨深度變化存在指數(shù)關(guān)系，該關(guān)系不受地層遭受剝蝕影響，通過現(xiàn)有聲波時差數(shù)據(jù)擬合可以建立標(biāo)準(zhǔn)指數(shù)曲線，曲線上各點所對應(yīng)的埋藏深度值與現(xiàn)有深度值的差就是所經(jīng)歷的剝蝕量值。本文利用聲波時差法進行剝蝕量恢復(fù)［14］。

2.3.2 理論曲線的確定

理論上，僅考慮機械壓實作用的沉積物孔隙度與深度呈負指數(shù)遞減關(guān)系變化［15－22］。

（1）初始孔隙度φ0確定。其中，track為分選系數(shù)，track＝P25／P75，P25和P75由粒度累計曲線上25%和75%所對應(yīng)的顆粒直徑求得。分選系數(shù)越低，砂巖的分選性越好，其原始孔隙度越高。對蘇里格地區(qū)48個砂巖樣品的分選系數(shù)進行計算（見圖3）求得分選系數(shù)在1.45～2.86之間，平均分選系數(shù)為1.85，帶入M.Scherer公式，即

求得初始孔隙度為33.2%。

圖3 分選系數(shù)頻率圖

（2）計算壓實系數(shù)。因條件限制，不能實測到足夠多的樣品，又考慮到聲波時差值反映均勻分布的原始孔隙度，所以使用自行開發(fā)的計算孔隙度的Simple軟件，利用測井聲波時差值對孔隙度進行計算，并進行流體校正，求取不同深度砂巖孔隙度值。

在恢復(fù)研究井剝蝕量的基礎(chǔ)上，對蘇82井、召60井、召61井、召65井、統(tǒng)40井、蘇102井、蘇54井、蘇77井、蘇52井等9口井砂巖孔隙度－深度數(shù)據(jù)進行擬合，求取壓實系數(shù)的平均值為0.000 4，以該值作為蘇里格地區(qū)壓實系數(shù)c值。

2.3.3 理論孔隙度－深度曲線與實測數(shù)據(jù)對比和成巖相劃分標(biāo)準(zhǔn)的確定

把理論曲線和實測孔隙度－深度曲線進行對比。通過巖石薄片描述對比分析，實測數(shù)據(jù)位于理論曲線的左側(cè)，砂巖薄片顯示為膠結(jié)作用，實測數(shù)據(jù)位于理論曲線右側(cè)，砂巖薄片顯示為溶蝕作用見圖4。

根據(jù)偏移程度，建立成巖相劃分標(biāo)準(zhǔn)

（1）石英砂巖成巖相劃分指標(biāo)（簡稱偏離度）：偏離度＝實測孔隙度值－理論孔隙度值

（2）根據(jù)偏離度，石英砂巖的成巖相可以劃分為6種（見表1）。

圖4 蘇39井理論曲線和實測孔隙度－深度曲線對比圖

表1 成巖相劃分標(biāo)準(zhǔn)表

2.3.4 壓實法石英砂巖成巖相劃分結(jié)果

采用壓實法對有掃描電鏡資料的井段進行了石英砂巖成巖作用分析。通過理論孔隙度與實測孔隙度的比較可以確定成巖相，其結(jié)果與根據(jù)掃描電鏡資料所確定的砂巖的成巖作用類型基本一致（見表2）。

表2 蘇里格地區(qū)成巖相識別效果對照表

3 結(jié) 論

（1）提出了一種石英砂巖成巖相的有效識別方法壓實法，利用測井資料，在地層埋藏史恢復(fù)的基礎(chǔ)上，通過回歸分析方法對石英砂巖成巖相進行了定量的判別。

（2）壓實法應(yīng)用于蘇里格地區(qū)的石英砂巖成巖相的識別，建立了蘇里格地區(qū)石英砂巖成巖相判別標(biāo)準(zhǔn)。識別出6類成巖相類型為強、中、弱溶蝕相和強、中、弱膠結(jié)相。識別結(jié)果與根據(jù)掃描電鏡資料所確定的砂巖成巖相類型進行對比驗證，證明了該方法識別結(jié)果準(zhǔn)確，具有較好的應(yīng)用前景。

（3）提高了石英砂巖測井成巖相識別能力和定量評價精度，進而為蘇里格地區(qū)的深入勘探與進一步開發(fā)提供了科學(xué)依據(jù)。

［1］ 張小莉，沈英，陳文學(xué).利用測井資料分析成巖作用對儲集層的影響［J］.沉積學(xué)報，2008，18（1）：127－137.

［2］ 段新國，宋榮彩，李國輝，等.四川盆地須二段綜合成巖相特征研究［J］.西南石油大學(xué)學(xué)報：自然科學(xué)版，2011，1（33）：8－16.

［3］ 劉紅岐，彭仕宓，唐洪.蘇里格氣田低阻氣層分布特征及識別方法研究［J］.天然氣工業(yè)，2004，10（24）：41－44.

［4］ 何自新，付金華，席勝利，等.蘇里格大氣田成藏地質(zhì)特征［J］.石油學(xué)報，2003，24（2）：6－12.

［5］ 李清，楊鎖衛(wèi)，喻智勇，等.提高蘇里格氣田鉆井速度的探討［J］.西部探礦工程，2010，4：74－76.

［6］ 王琪，史基安，肖立新，等.石油侵位對碎屑儲集巖成巖序列的影響及其與孔隙演化的關(guān)系：以塔西南坳陷石炭系石英砂巖為例［J］.沉積學(xué)報，1998，16（3）：97－101.

［7］ 丁曉琪，張哨楠.鄂爾多斯盆地西南緣中生界成巖作用及其對儲層物性的影響［J］.油氣地質(zhì)與采收率，2011，1（18）：18－23.

［8］ 寧寧，陳孟晉，劉銳娥，等.鄂爾多斯盆地東部上古生界石英砂巖儲層成巖及孔隙演化［J］.天然氣地球科學(xué)，2007，18（3）：334－338.

［9］ 鄭瑞林.陜甘寧盆地煤系地層中石英砂巖成巖作用及其孔隙演化［J］.石油勘探與開發(fā)，1989，6（3）：31－40.

［10］關(guān)維東，李興訓(xùn).臺3一北10井區(qū)侏羅系貧石英砂巖的成巖作用及儲層特征研究［J］.新疆石油地質(zhì)，1991，12（2）：173－179.

［11］劉偉，竇齊豐，黃述旺，等.成巖作用的定量表征與成巖儲集相研究：以科爾沁油田交2斷塊區(qū)九佛堂組（J3jf）下段為例［J］.中國礦業(yè)大學(xué)學(xué)報，2002，5（31）：399－403.

［12］李壯福，郭英海，李熙哲，等.王府凹陷儲集層成巖作用及對儲集性能的影響［J］.中國礦業(yè)大學(xué)學(xué)報：自然科學(xué)版，2001，1（30）：94－98.

［13］張文忠.蘇里格地區(qū)上古生界氣田異常低壓成因研究［D］.北京：中國地質(zhì)大學(xué)，2009.

［14］彭清華，周江羽，揭異新.聲波時差法對南堡凹陷東營組剝蝕量的恢復(fù)［J］，斷塊油氣田，2009，6（16）：50－53.

［15］李傳亮，孔祥言.孔－深關(guān)系曲線方程研究［J］.新疆石油地質(zhì)，2001，2（22）：5－6.

［16］Vidar Storvoll，Knut Bj?rlykke，Nazmul H.Mondol Velocity－depth Rrends in Mesozoic and Cenozoic Sediments from the Norwegian Shelf［J］.AAPG Bulletin，2005，3（89）：359－381.

［17］Roberto Aguilera.Sandstone vs.Carbonate Petroleum Reservoirs：A global Perspective on Porosity－depth and Porosity Permeability Relation－ships：Discussion［J］.AAPG Bulletin，2006，5（90）：807－810.

［18］Paul H Nadeau，Stephen N Ehrenberg.Sandstone vs.Carbonate Petroleum Reservoirs：A Global Perspective on Porosity－depth and Porosity－permeability Relationships：Reply［J］.AAPG Bulletin，2006，5（90）：811－813.

［19］Ehrenberg S N，Nadeau P H，?.Steen.Petroleum Reservoir Porosity Versus Depth：Influence of Geological Age［J］.AAPG Bulletin，2009，10（93）：1281－1296.

［20］Olav Walderhaug.Modeling Quartz Cementation and Porosity in Middle Jurassic Brent Group Sandstones of the Kvitebj?rn Field，Northern North Sea［J］.AAPG Bulletin，2000，9（84）：1325－1339.

［21］Paxton S T，Szabo J O，Ajdukiewicz J M，Klimentidis R E.Construction of an Intergranular Volue Compaction Curve for Evaluating and Predicting Compaction and Porosity Loss in Rigid－grain Sandstone Reservoirs［J］.AAPG Bulletin，2002，12（86）：2047－2067.

［22］?yvind Marcussen，Tom Erik Maast，Nazmul H Mondol，et al.Changes in Physical Properties of a Peservoir Sandstone as a Function of Burial Depth－The Etive Formation，Northern North Sea［J］.Marine and Petroleum Geology，2010：1725－1735.

Recognition of Quartz Sandstone Diagenetic Facies with Compaction Method

WANG Xin1，XUE Linfu1，YU Yunliang1，BAI Ye1，SHI Yujiang2，ZHANG Haitao2
（1.Earth Science College，Jilin University，Changchun，Jilin 130062，China；2.Institute of Exploration and Development，Changqing Oilfield Company，PetroChina，Xi’an，Shaanxi 710018，China）

Proposed is an effective method called compaction method which is used to recognize the diagenetic facies of quartz sandstone with log data.On the basis of denudation quantity recovery，this method creates a porosity－depth relations curve as standard.The curve only reflects mechanical compaction by regression analysis method.Then the values of the theoretical and actual porosity are compared.Combined with core observation，thin section identification and electronic probe spectrum，we set the division standard of the diagenetic facies using the compaction method.According to the division standard，the diagenetic facies of quartz sandstone in no coring wells is identified by calculated porosity after correction of fluid.This method is used in Sulige area.The diagenetic facies of quartz sandstone can be divided into six types：weak／

moderate／strong dissolution，weak／moderate／strong cementation.The above results show that the method improves the identification ability and quantitative evaluation of quartz sandstone diagenetic facies using log data，and so has good prospect.

log interpretation，quartz sandstone，compaction method，diagenetic facies，diagenesis

P539.2

A

2012－05－03 本文編輯 李總南）

1004－1338（2012）04－387－05

國家自然科學(xué)基金（41174096）和國家重大專項基金（2011ZX05009、2011ZX05044）聯(lián)合資助

王欣，女，碩士研究生，從事儲層研究。