賀振華,賈義蓉,蔣煉,黃德濟(jì)

(1.成都理工大學(xué) 油氣藏地質(zhì)及開發(fā)工程國家重點實驗室,成都 610059;2.成都理工大學(xué) 信息工程學(xué)院,成都 610059)

碳酸鹽巖礁灘油氣儲層地震預(yù)測方法探討

賀振華1,2,賈義蓉2,蔣煉2,黃德濟(jì)2

(1.成都理工大學(xué) 油氣藏地質(zhì)及開發(fā)工程國家重點實驗室,成都 610059;2.成都理工大學(xué) 信息工程學(xué)院,成都 610059)

儲層結(jié)構(gòu)和孔隙流體預(yù)測是目前碳酸鹽巖礁灘油氣儲層地震預(yù)測的重點和難點。這里從流體敏感性參數(shù)的選擇和基于儲層結(jié)構(gòu)模擬的孔隙度預(yù)測二方面,研究了碳酸鹽巖礁灘油氣儲層的流體識別問題。其中,基于測井資料統(tǒng)計分析和/或巖石物理巖樣測試的參數(shù)交會圖的制作,是優(yōu)選烴類敏感參數(shù)的基礎(chǔ)。在單一敏感參數(shù)的基礎(chǔ)上,構(gòu)組復(fù)合型的流體識別因子,能獲得更好的流體識別效果。儲層孔隙度預(yù)測是這樣實現(xiàn)的:首先,對Gassm an流體替換方程通過引入近似關(guān)系βp-βs≈βp進(jìn)行簡化;然后引入Eshelby-W alsh儲層結(jié)構(gòu)參數(shù)以獲得直接計算孔隙度的表達(dá)式;最后,根據(jù)彈性反演得到的縱波、橫波阻抗(或縱波、橫波速度)等參數(shù),計算得到孔隙度及孔隙流體預(yù)測剖面。經(jīng)實際地震資料的流體預(yù)測結(jié)果顯示,新方法比常規(guī)方法預(yù)測的精度高。

礁灘儲層;流體預(yù)測;地震孔隙度反演;孔隙結(jié)構(gòu)

0 概述

碳酸鹽巖油氣儲層的預(yù)測,主要包括縫洞型和礁灘型,或二者的復(fù)合型儲層的預(yù)測。在這里,我們主要涉及的是生物礁灘型油氣儲層的識別和預(yù)測問題。

生物礁是碳酸鹽巖隆的一種類型,由各種抗浪耐沖擊的,具有極強生命力的造礁生物組成的特殊碳酸鹽巖建造。礁、灘巖體相對圍巖具有較高的孔隙度和滲透率,是石油和天然氣的良好儲存場所,具有油氣豐度大、產(chǎn)能高的明顯特征。世界上曾經(jīng)達(dá)到日產(chǎn)萬噸級油氣的井總共只有九口井,八口井在碳酸鹽巖儲層中,而其中四口井在生物礁儲層中。由此可見,生物礁油氣儲層在油氣勘探開發(fā)中占有十分重要的地位。

我國石油工業(yè)部門在四川盆地東北部、塔里木盆地和南海珠江口等盆地,先后發(fā)現(xiàn)了不同類型和特征的大型礁灘油氣田,有巨大的潛在資源量。因此,加強生物礁、灘儲層的勘探開發(fā),對緩解我國能源供求緊張的局面,具有現(xiàn)實而又長遠(yuǎn)的戰(zhàn)略意義。

然而,生物礁灘儲層屬于隱蔽型巖性油氣藏或者巖性～構(gòu)造復(fù)合油氣藏,有相當(dāng)大的預(yù)測難度。在以二維地震資料為主的勘查階段,預(yù)測生物礁灘儲層的鉆探成功率很低,據(jù)劉劃一等人[1、2]的資料顯示,在八十年代初投入了較多的礁灘勘探工作量以后,至1997年底在已知的三十口礁井中,屬于正確預(yù)測的只有10%。目前,隨著三維地震資料的廣泛應(yīng)用,新的預(yù)測方法不斷涌現(xiàn),以及勘探精度的逐步提高,預(yù)測生物礁灘儲層的鉆探成功率已有大大提高[3～5],有的鉆探成功率已達(dá)90%以上,但這只是指鉆遇生物礁體本身而言。而生物礁灘儲層的內(nèi)部結(jié)構(gòu)和含流體性質(zhì)的預(yù)測依然十分困難,并且鉆遇干井和水井的風(fēng)險很大。主要原因在于:

(1)儲層結(jié)構(gòu)和流體預(yù)測本身是國內(nèi)、外的高難度研究課題,探索性強。

(2)碳酸鹽巖礁灘儲層內(nèi)部巖性和孔隙流體結(jié)構(gòu)的非均質(zhì)性較其它儲層更強,與圍巖的差異更小,使得地震信號雜亂,甚至空白。

(3)目前儲層識別和流體預(yù)測的一些有效方法和統(tǒng)計關(guān)系,大多是針對碎屑巖儲層的,直接利用有一定困難。

(4)碳酸鹽巖儲層往往埋藏深,這一方面限制了多波和大出射角地震資料和相關(guān)先進(jìn)方法的利用;另一方面使地震資料的分辨率降低。

因此,發(fā)展有針對性的預(yù)測方法是十分必要的。這里,我們根據(jù)近年來在生物礁灘型油氣儲層識別和流體預(yù)測方面的實踐,探討一些具體問題,主要包括流體敏感性預(yù)測參數(shù)的選擇和基于儲層結(jié)構(gòu)模擬的孔隙流體預(yù)測方法。

1 流體敏感性預(yù)測參數(shù)的選擇

基于地震資料儲層流體預(yù)測的一個重要前提是選擇對流體比較敏感的物性參數(shù)和地震響應(yīng)屬性。所謂敏感參數(shù)系是指能將儲層和非儲層,油氣層和干層,油氣層和水層有效區(qū)分的參數(shù)。其主要方法有:

(1)基于測井資料和/或巖石物理樣品測試分析的流體敏感參數(shù)的選擇。

(2)應(yīng)用物理和數(shù)學(xué)分析方法,構(gòu)組高靈敏度流體識別因子,提高地震反演預(yù)測流體的可靠性。

1.1 基于測井資料或巖石物理樣品測試分析的流

體敏感參數(shù)選擇

首先,要根據(jù)任務(wù)要求,并利用各種參數(shù)的交會圖對有關(guān)參數(shù)逐一分析,最后確定有利于流體識別的敏感參數(shù)或參數(shù)組合。由于生物礁灘儲層的復(fù)雜性,不同地區(qū)的儲層,其流體敏感參數(shù)可能很不相同;即使在同一地區(qū),同一目的層,各個井的敏感參數(shù)也會存在較大差異,需要全面分析,仔細(xì)選擇。

圖1為ZJ地區(qū)礁灘儲層有效孔隙度和密度的測井參數(shù)交會圖。由圖1可見,在相同孔隙度情況下,利用密度可區(qū)分碳酸鹽巖和碎屑巖層,但很難識別油層和水層。

圖2(見下頁)為ZJ地區(qū)礁灘儲層中有效孔隙度與橫波速度、縱橫波阻抗差、縱橫波阻抗平方差和自然伽瑪?shù)慕粫D。其中,在相同孔隙度且孔隙度小于20%的情況下,利用縱橫波阻抗差和縱橫波阻抗平方差有可能區(qū)分油～水層。而橫波速度和自然伽瑪不能可靠區(qū)分。

圖1 ZJ地區(qū)礁灘儲層有效孔隙度和密度的測井參數(shù)交會圖Fig.1 C rossp lotof effective porosity and density for carbonate reef-shoal reservoirwell-log data in ZJ area

1.2 復(fù)合型流體識別因子的構(gòu)組

通過基本的彈性參數(shù)λ(拉梅系數(shù))、μ(剪切模量)、ρ(密度)和縱橫波阻抗Ip、Is可以構(gòu)組更多、更敏感的復(fù)合型流體識別因子,以獲得更好的流體識別效果,如公式(1)～公式(5)所示。

下頁圖3～后面圖5展示了部份流體識別因子對氣～水層的潛在識別能力。

(1)泊松比。

(2)Goodw ay(1997)等提出的識別因子。

(3)Russell(2003)等提出的流體屬性。

式中 f為流體因子;c為調(diào)節(jié)參數(shù)。

(4)賀振華等提出的高靈敏度流體識別因子

式中 B為調(diào)節(jié)參數(shù)。當(dāng)B=2 Ip/Is時,得

式(5)結(jié)合了Ip/Is和的優(yōu)點。

2 基于儲層孔隙結(jié)構(gòu)模擬的孔隙流體預(yù)測方法探討

我們將Gassm an流體替換方程和由Eshelby 1957年提出,后經(jīng)W alsh簡化的方程[6,7]相結(jié)合,以便在考慮儲層結(jié)構(gòu)的同時,預(yù)測孔隙流體。因為孔隙結(jié)構(gòu)和孔隙度,對儲層的影響都很重要。

圖2 ZJ地區(qū)礁灘儲層流體敏感參數(shù)的測井統(tǒng)計分析Fig.2 Crossp lotsof effective porosity and shearwave velocity(a)difference between P and S impedance(b)square difference of P and S impedance(c)Gamm a curves from ZJwell-log data

圖3 L地區(qū)礁灘儲層高靈敏度流體識別因子(HSFIF)與縱波阻抗的交會圖Fig.3 Crossp lotof high sensitive fluid factor(HSFIF)and Pwave's impedance from carbonate reef-shoal reservoir in L area

圖4 L地區(qū)礁灘儲層泊松比與縱波阻抗的交會圖Fig.4 Crossp lotof high sensitive fluid factor and poisson'sratio from the area sam e as Fig.3

2.1 方法原理

Gassm an流體替換方程[8～12]如下:

式中 βs、βD、βp、ˉβ和η分別為含流體雙相介質(zhì)的基質(zhì)(骨架)壓縮系數(shù)、干燥(含空氣)巖石的壓縮系數(shù)、孔隙流體的壓縮系數(shù)、含孔隙流體儲層的有效壓縮系數(shù)和儲層的孔隙度。考慮到孔隙流體(油、氣、水)的壓縮系數(shù),一般大于巖石基質(zhì)(骨架)的壓縮系數(shù)達(dá)一個數(shù)量級以上,可得βp-βs≈βp的近似關(guān)系(對碳酸鹽巖來說,其基質(zhì)壓縮系數(shù)比碎屑巖的小一倍以上,可使該項近似的可行性更好!)。于是,式(6)可簡化為:

圖5 L地區(qū)礁灘儲層縱波、橫波阻抗的交會圖Fig.5 Crossp lotof P-wave and S-wave impedances in L area

由式(7),得

由Eshelby和W alsh的干燥巖石橢球包體近似公式[5、6]:

式中 m為干燥巖石橢球包體的結(jié)構(gòu)參數(shù)之一;α是另一結(jié)構(gòu)參數(shù),表示橢球孔隙或裂隙的縱橫(或長短軸)比;βD在式(8)和式(9)中具有相同的涵義。將式(9)代入式(8),經(jīng)化簡后得式(10)。

由A、B的表達(dá)式,很易得到:

ˉβ和η為飽和流體的壓縮系數(shù)和孔隙度,可由巖石物理測試分析或測井參數(shù)統(tǒng)計獲得。βp的參考值可用公式計算[16、17],βs為巖石基質(zhì)壓縮系數(shù),一般較難得到。但通過式(10)的統(tǒng)計分析和線形擬合,可求出A、B,再按式(11),很易求得βs。上述關(guān)系式對地震彈性正演、反演和流體識別,有著重要參考價值。

2.2 應(yīng)用效果

圖6 基于礁灘儲層內(nèi)部孔隙結(jié)構(gòu)模擬的孔隙度預(yù)測剖面(紅黃色的中低孔隙度指示含油層)Fig.6 Comparison between new porosity inversion(a)and traditional Pwave velocity inversion(b).The resultof porosity inversion is better than the one in(b)

通過地震疊前彈性參數(shù)反演,首先獲得縱橫波阻抗或縱橫波速度,然后利用上述公式獲得考慮了孔隙結(jié)構(gòu)的孔隙度和孔隙流體反演剖面。圖6為某區(qū)通過Ky1井的孔隙度地震反演剖面。圖6(b)為同一資料的常規(guī)疊前縱波速度反演剖面。剖面中的紅黃色標(biāo)代表中低孔隙度和中低縱波速度,是油層存在的標(biāo)志。二者相比,孔隙度反演的結(jié)果與Ky1井資料的吻合度比縱波速度反演結(jié)果的吻合程度更好。

3 結(jié)論

(1)流體敏感性預(yù)測參數(shù)的選擇,對碳酸鹽巖礁灘油氣儲層以及其它類型儲層的流體預(yù)測都是十分重要的。其中,根據(jù)測井參數(shù)和巖石物理測試參數(shù)所做的交會圖,是進(jìn)行流體敏感參數(shù)選擇的基礎(chǔ)性圖件。在制作不同參數(shù)的交會圖件時應(yīng)當(dāng)盡可能全面,以便挑選出真正的敏感性參數(shù)或參數(shù)組合。

(2)復(fù)合型流體識別因子較單一的敏感參數(shù)有更好的流體識別效果,而且構(gòu)組方法簡單。讀者可根據(jù)研究任務(wù)的需要,自行構(gòu)組更多、更有效的復(fù)合型流體識別因子。

(3)基于儲層內(nèi)部孔隙結(jié)構(gòu)模擬的孔隙度預(yù)測(反演)方法的優(yōu)點是考慮了孔隙結(jié)構(gòu)這一重要因素,有利于儲層基質(zhì)壓縮系數(shù)的估算,并見到了一定的實際應(yīng)用效果。但仍然存在孔隙結(jié)構(gòu)參數(shù)估算的困難,而且基質(zhì)壓縮系數(shù)的估算值要隨孔隙流體的性質(zhì)變化,這些問題還需要進(jìn)一步認(rèn)識和探索。

[1]劉劃一,王一剛,楊雨,等.川東上二疊統(tǒng)生物礁氣藏多元信息綜合預(yù)測方法研究[J].天然氣工業(yè),1999,19(4):13.

[2]陳太源.九龍山構(gòu)造南側(cè)存在長興組臺緣礁的探討[J].天然氣工業(yè),1989,9(1):6.

[3]賀振華,蒲勇,熊曉軍,等.川東北長興～飛仙關(guān)組礁灘儲層的三維地震識別[J].物探化探計算技術(shù),2009,31(1):1.

[4]賀振華,黃德濟(jì),文曉濤.裂縫油氣藏地球物理預(yù)測[M].成都:四川科學(xué)技術(shù)出版社,2007.

[5]牟書令.中國海相油氣勘探理論技術(shù)與實踐[M].北京:地質(zhì)出版社,2009.

[6]陳颙,黃庭芳.巖石物理學(xué)[M].北京:北京大學(xué)出版社,2001.

[7]陳颙,黃庭芳,劉恩儒.巖石物理學(xué)[M].合肥:中國科學(xué)技術(shù)大學(xué)出版社,2009.

[8]SHER IFFR E.ReservoirGeophysics[M].SEG,1992.

[9]AMOSNUR.雙相介質(zhì)中波的傳播[M].許云,譯.北京:石油工業(yè)出版社,1986.

[10]PER AVSETH,TAPAN MUKERJI,GARY MAVKO.Quantitative Seism ic In terp retation—App lying Rock Physics Too ls to Reduce Interp retation R ick[M].Cambridge University Press,2005.

[11]LEON THOM SEN.Understanding seism ic anisotropy in exp loration and exp loitation[M].SEG,2002.

[12]MAVKO G,MUKERJIT,DVORK IN J.The rock physics handbook-o ls for seism ic in porousm edia[M].Cam bridgeUniversity Press,2003.

[13]RUSSELL B H.The app lication ofm u ltivariate statisticsand neu ralnetwo rks to the p red iction of reservoirparam eters using seism ic attributes.Calgary[M].Departm entof Geology and Geophysics,Canada,2004.

[14]CONNOLLY.Elastic impedance[J].The Leading Edge,1999,18(4):438.

[15]王之敬(W ang Zee).“巖石物理概述”學(xué)術(shù)報告(PPT)[D].2002.

[16]鄒文.基于地震資料的流體識別技術(shù)研究[D].成都:成都理工大學(xué),2008.

[17]劉雯林.油氣田開發(fā)地震技術(shù)[M].北京:石油工業(yè)出版社,1996.

TE 122.2+21

A

1001—1749(2011)01—0001—05

國家自然科學(xué)基金資助(40739907,40774064)

2010-08-23

賀振華(1938-),男,湖北大悟縣人,教授,博士導(dǎo)師,現(xiàn)主要從事油氣地球物理和巖石物理等領(lǐng)域的教學(xué)與研究工作。