陳培元，段曉夢，郭麗娜，倪軍娥

（中海油研究總院，北京100028）

?

伊拉克B油田白堊系Mishrif組碳酸鹽巖儲層分類

陳培元，段曉夢，郭麗娜，倪軍娥

（中海油研究總院，北京100028）

伊拉克B油田白堊系Mishrif組為厚層礁灘相碳酸鹽巖，受成巖作用影響儲層非均質(zhì)性極強，采用常規(guī)的儲層分類方法難以準確的刻畫儲層孔隙結(jié)構(gòu)之間的差異性。根據(jù)取心井實測孔隙度和滲透率資料，采用W inland R35將研究區(qū)儲層劃分為五種類型，每種類型儲層孔隙度和滲透率之間相關(guān)性較好。結(jié)合毛管壓力和相滲曲線，研究發(fā)現(xiàn)不同儲層類型對應(yīng)的毛管壓力特征和相滲特征差異較大，充分體現(xiàn)了碳酸鹽巖儲層受成巖作用影響非均質(zhì)性強的特征。儲層類型的劃分對于構(gòu)建精細的地質(zhì)模型和提高數(shù)值模擬的精度具有重要的意義。

儲層分類；Winland R35；碳酸鹽巖；伊拉克

碳酸鹽巖儲層在世界油氣田分布中占有十分重要的地位，據(jù)統(tǒng)計在碳酸鹽巖儲層中發(fā)現(xiàn)的油氣儲量已接近全球總儲量的一半，油氣產(chǎn)量則占全球總產(chǎn)量的60%以上［1］。世界重要產(chǎn)油氣區(qū)的儲集層都是以碳酸鹽巖儲層為主，如墨西哥環(huán)礁帶、波斯灣盆地、利比亞的錫爾特盆地、北美地臺區(qū)的密執(zhí)安盆地、二疊盆地、伊利諾斯盆地、西部盆地和辛辛那提隆起以及阿爾伯達地區(qū)等［2，3］。在國內(nèi)，隨著許多大中型油氣田，如塔里木盆地奧陶系油氣藏、四川盆地油氣區(qū)、鄂爾多斯盆地奧陶系油氣藏等［4-6］相繼發(fā)現(xiàn)和開發(fā)，對該類油氣藏開展的研究也越來越多。與碎屑巖儲層相比，碳酸鹽巖地層沉積時代久遠，經(jīng)歷的成巖作用期次多、類型廣，致使其各向異性變大、非均質(zhì)性變強，許多成功的應(yīng)用在碎屑巖儲層中的評價方法在碳酸鹽巖儲層中卻難以適用。因此，碳酸鹽巖儲層評價成為一直以來困擾該類油氣藏勘探開發(fā)的主要技術(shù)難題之一。本文以伊拉克B油田白堊系碳酸鹽巖儲層為例，通過采用目前碳酸鹽巖油藏最常用的方法Winland R35對研究區(qū)儲層進行了分類評價，為油藏地質(zhì)建模和數(shù)值模擬提供了依據(jù)。

1　地質(zhì)概況

B油田位于伊拉克東南部，毗鄰伊朗邊界，南距巴士達約175 km，西北距巴格達約350 km。構(gòu)造上，位于Zagros山前陸褶皺帶中部，為一NW-SE向的狹長背斜［7］。區(qū)域沉積表明，受早阿爾卑斯運動的影響，美索不達米亞地區(qū)在白堊紀時期主要以淺海碳酸鹽巖沉積為主，只有在Nahr Umr組形成了廣泛的碎屑巖沉積［8，9］。白堊系Mishrif組是研究區(qū)主要的生產(chǎn)層，與上覆Khasib組呈不整合接觸，與下伏Rumaila組呈整合接觸［10］，地層厚度320 m～370 m，自上而下Mishrif組可以劃分為MA1、MA2、MB1、MB2、MC1及MC2六個油組及對應(yīng)的六個亞油組，其中MB21亞油組是研究區(qū)的主力產(chǎn)層段。區(qū)域沉積表明，Mishrif組沉積期，阿拉伯地區(qū)發(fā)生區(qū)域構(gòu)造活動，形成地貌高地和低地，大部分高地成為有機體建造發(fā)育部位，其中研究區(qū)主要發(fā)育緩坡碳酸鹽巖臺地相［7］。取心井巖心及鏡下薄片統(tǒng)計分析表明，Mishrif組儲層段以顆粒灰?guī)r為主，次為泥晶灰?guī)r。

圖1　汞飽和度35%時孔喉半徑示意圖及研究區(qū)巖心孔喉分布曲線

2　Winland R35評價方法

Winland R35評價方法又稱孔喉幾何形狀法，是由Alden等提出的一種利用孔喉半徑R35來劃分和評價巖石物理流動單元的方法［11］。Winland R35方法描述了地層巖石孔喉尺寸與孔隙度和滲透率的關(guān)系，該方法不考慮地層中巖石的沉積過程，直接通過巖石孔喉結(jié)構(gòu)劃分巖石類型，反映了不同巖石類型地層在目前狀態(tài)下的滲流能力［12］。Winland研究表明壓汞曲線上進汞達35%時的孔喉半徑R35大小可以反映該巖石中的流體流動和開發(fā)動態(tài)（見圖1a）?？梢詫悠返膲汗€分析求取R35值，但是樣品的壓汞分析資料畢竟有限，不易在全區(qū)內(nèi)統(tǒng)計，因此Winland提出利用測井曲線進行R35值的求取，Winland方程：

Martin等根據(jù)孔喉半徑R35值的分布范圍將儲層分為四種類型［13］（見圖1b）：（1）巨孔喉儲層：R35＞10μm；（2）大孔喉儲層：2.0μm＜R35≤10.0μm；（3）中孔喉儲層：0.5μm＜R35≤2.0μm；（4）微孔喉儲層：R35≤0.5μm。

孔喉的大小是控制流體流動的主要因素之一，特別是較粗孔喉對流體流動起著主導(dǎo)作用。因此可以利用Winland方程計算求得的孔喉半徑R35來劃分儲層類型。但是在應(yīng)用Winland方程時必須要注意其使用條件，如果儲層中具有連通的溶洞和裂縫等次生孔隙時，計算得到的R35會過大，失去其真實性［13］。

3　儲層類型劃分

通過對研究區(qū)取心井實測物性數(shù)據(jù)分析發(fā)現(xiàn)，研究區(qū)孔隙度和滲透率之間受成巖作用影響相關(guān)性差。采用R35方法對取心井1 928塊巖心樣品進行分析，結(jié)果（見圖2）。整體上研究區(qū)Mishrif組巖心樣品的孔喉分布在0.5μm～10.0μm，僅有少量的樣品孔喉半徑大于10μm。

圖2　研究區(qū)Mishrif組巖心孔滲分布圖

圖3　基于R35方法的孔滲分類圖

圖4　不同R35對應(yīng)的儲層鏡下特征

根據(jù)Alden的劃分方案，將研究區(qū)儲層劃分為四類，每種儲層類型的孔滲擬合情況（見圖3a），從圖上可以看出每種儲層類型的擬合結(jié)果，整體上采用R35分類擬合的誤差較小，但對于大孔喉儲層（2.0μm＜R35≤10.0μm）孔滲之間相關(guān)性較差，方差僅為0.375 6。進一步通過鏡下觀察發(fā)現(xiàn)（見圖4），2.0μm＜R35與5.0μm＜R35的儲層盡管滲透率差異不大，但孔隙結(jié)構(gòu)存在著明顯的差異，主要體現(xiàn)在5.0μm＜R35的儲層孔喉要顯著大于2.0μm＜R35的儲層?；诖耍瑢?.0μm＜R35≤10.0μm的儲層進一步細分為2.0μm＜R35≤5.0μm與5.0μm＜R35≤10.0μm兩種類型，重新對兩種類型的儲層孔滲關(guān)系進行擬合，發(fā)現(xiàn)兩者之間的相關(guān)性明顯得到了提高（見圖3b）。

為了更好的驗證采用R35方法進行分類的可行性，用擬合的孔滲關(guān)系式計算了取心井的巖心滲透率，并與實際分析得到的滲透率數(shù)據(jù)進行了對比。從對比圖上可以看出（見圖5），采用R35預(yù)測的滲透率與實測滲透率數(shù)據(jù)沿斜率為1的直線分布，兩者的吻合率較高。

圖5　實測滲透率與R35預(yù)測滲透率交匯圖

4　儲層類型的應(yīng)用

碳酸鹽巖沉積物與碎屑巖具有極大的差別，在成巖過程中極易受成巖流體和構(gòu)造作用的影響，形成大量的溶蝕孔、洞和縫，導(dǎo)致碳酸鹽巖儲層可以在較小的范圍內(nèi)展現(xiàn)出較大的儲層物性（孔隙度，滲透率）差異［14］，難以采用常規(guī)的碎屑巖“相控”思路進行建模。此外，碳酸鹽巖儲層孔隙結(jié)構(gòu)類型復(fù)雜，使得儲存流體和供流體通過的孔隙空間結(jié)構(gòu)比較復(fù)雜，從而導(dǎo)致碳酸鹽巖儲層滲透率的定量計算難以準確表達。

而儲層類型的劃分，可以從根本上解決以上兩個問題。在地質(zhì)建模過程中，可以將不同的儲層類型定義為不同的“儲集相”。由于不同的儲層類型可以反映不同的孔隙結(jié)構(gòu)特征，對應(yīng)的毛管壓力曲線和相滲曲線也不盡相同（見圖6），這樣就為油藏地質(zhì)建模和數(shù)值模擬提供了依據(jù)。通過模擬不同儲層類型的“儲集相”和“孔隙度”模型，根據(jù)擬合的孔滲關(guān)系式可以準確的計算滲透率，從而構(gòu)建滲透率模型。這樣一方面可以提高模型的精度，另一方面可以反映碳酸鹽巖儲層的強非均質(zhì)性。

圖6　不同類型儲層對應(yīng)的毛管壓力曲線及相滲曲線

5　結(jié)論

（1）碳酸鹽巖儲層受成巖作用的影響具有較強的非均質(zhì)性，孔滲之間相關(guān)性差，單純的根據(jù)孔隙度和滲透率難以進行儲層類型的劃分。

（2）采用Winland R35方法可以將研究區(qū)碳酸鹽巖儲層劃分為五種類型，其中不同類型的儲層對應(yīng)的孔喉半徑具有明顯的差別。

（3）根據(jù)儲層類型劃分結(jié)果可以進行碳酸鹽巖“儲集相”建模，進而建立對應(yīng)的滲透率模型，提高了地質(zhì)建模和數(shù)值模擬的精度。

［1］Roehl PO.Carbonate petroleum reservoirs［M］.Berlin，Heidelberg，New York：Springer Verlag，1985：1-15.

［2］Alsharhan A S，Nairn A E M.Sedimentary basins of petroleum geology of the Middle East［M］.Amsterdam：Elsevier Science，1997：1-942.

［3］Wilson J L.Limestone and dolomite reservoir［A］.Hobson G D.Developments in petroleum geology-2［M］.London：App lied Science Publishers Ltd，1980：1-51.

［4］羅平，張靜，劉偉，等.中國海相碳酸鹽巖油氣儲層基本特征［J］.地學前緣，2008，15（1）：36-50.

［5］趙宗舉，范國章，吳興寧，等.中國海相碳酸鹽巖的儲層類型、勘探領(lǐng)域及勘探戰(zhàn)略［J］.海相油氣地質(zhì)，2007，12（1）：1-11.

［6］吳興寧，孫六一，于洲，等.鄂爾多斯盆地西部奧陶紀巖相古地理特征［J］.巖性油氣藏，2015，27（6）：87-96.

［7］Aqrawi A A M.，Goff JC，Horbury A，et al.The petroleum geology of Iraq［M］.Scientific Press，Beaconsfield，UK，2010：424.

［8］Alsharhan A Sand Nairn A EM.A review of the Cretaceous formations in the Arabian Peninsula and Gulf：Part I．Lower Cretaceous（Thamama Group）stratigraphy and paleogeography［J］．Journal of Petroleum Geology，1986，9（4）：365-392.

［9］Alsharhan A S，Nairn A EM.A review of the ctretaceous formations in the Arabian peninsula and gulf：partii.mid-cretaceous（wasia group）stratigraphy and paleogeography［J］. Journal of Petroleum Geology，1988，11（1）：89-112.

［10］Aqrawi A A M，Thehni G A，et al.Mid-cretaceous rudistbearing carbonates of the Mishrif Formation：an important reservoir sequence in the Mesopotamian Basin［J］.Journal of Petroleum Geology，1998，（21）：57-82.

［11］Kolodzie S Jr.Analysis of pore throat size and use of the Waxmann-Smits equation to determine OOIP in Spindle Field，Colorado.Proc.Soc.Petrol.1980，Engn.55th Ann. Tech.Fall Conf［J］.Society of Petroleum Engineers，SPE 9382.

［12］譚學群，廉培慶.碳酸鹽巖油藏巖石分類方法研究［J］.科學技術(shù)與工程，2013，13（14）：3963-3967.

［13］Martin A J，Solomon S T and Hartmann D J.Characterization of petrophysical flow units in carbonate reservoirs［J］. AAPG Bulletin，1997，81：734-759.

［14］劉石磊，鄭榮才，顏文全，等.阿姆河盆地阿蓋雷地區(qū)牛津階碳酸鹽巖儲層特征［J］.巖性油氣藏，2012，24（1）：57-63.

Classification research on carbonate reservoir of the cretaceous M ishrif formation in B oilfield，Iraq

CHEN Peiyuan，DUAN Xiaomeng，GUO Lina，Ni Jun'e
（CNOOC Research Institute，Beijing 100028，China）

The cretaceous Mishrif formation of B oilfield ismainly consist of reef-shoal facies carbonate，and the reservoir heterogeneity is strong caused by the diagenesis.The conventional reservoir classification method is difficult to accurately depict the differences between reservoir pore structure.Porosity and permeability data of coring core samples are used by Winland R35method to divided reservoir types，which can be divided into five types in the study area，and each type of reservoir has good correlation between porosity and permeability.Combining with capillary pressure and relative permeability curves，it found that different reservoir with different capillary pressure and relative permeability characteristics，which can be fully embodies the carbonate reservoir heterogeneity affected by the diagenesis.The division of reservoir types has great importantmeaning to build fine geologicalmodel and improve the accuracy of numerical simulation.

reservoir classification；Winland R35；carbonate；Iraq

TE122.22

A

1673-5285（2016）07-0082-05

10.3969/j.issn.1673-5285.2016.07.020

2016－06-14

國家科技重大專項“亞太及南美地區(qū)復(fù)雜油田開發(fā)地質(zhì)關(guān)鍵技術(shù)研究”，項目編號：2011ZX05030-005。

陳培元，男（1984-），博士，研究方向為油氣田開發(fā)與地質(zhì)，郵箱：chenpy4@cnooc.com.cn。