孫盼科, 徐懷民, 黃婭, 姜貽偉

(1.中國石油大學(xué)地球科學(xué)學(xué)院, 北京 102249; 2.中國石油集團(tuán)測井有限公司油氣評價中心, 陜西 西安 710077;3.中國石化普光氣田分公司, 四川 達(dá)州 635000)

0 引 言

地震巖石物理研究[1-3]通過描述巖石彈性參數(shù)與巖石固體基質(zhì)、巖石骨架和孔隙流體這3種物理特征之間的關(guān)系建立儲集層特征參數(shù)變化與地震響應(yīng)特征(主要是速度)之間的定量解釋關(guān)系。Wyllie時間平均方程、Gassmann方程、Biot理論和K-T模型是巖石物理研究中重要的理論模型,常被用于不同巖石物理條件下定量關(guān)系的研究。巖石物理模型的構(gòu)建可簡單歸納為3步:①巖石基質(zhì)等效彈性模量的計算;②干巖石等效彈性模量計算;③飽和巖石的等效彈性模量的計算[4-6]。在這個過程中,干巖石等效彈性模量的求取一直是一個未完全解決的問題,雖然確定方法眾多,但每種方法都有其局限性[7-8]。

碳酸鹽巖儲集層巖石物理表征是通過建立地震速度與碳酸鹽巖儲集層巖石結(jié)構(gòu)、孔隙度和孔隙結(jié)構(gòu)類型等地質(zhì)參數(shù)之間關(guān)系進(jìn)行儲集層表征的一個過程。然而,碳酸鹽巖儲集層強(qiáng)的非均質(zhì)性使得這一實現(xiàn)過程又面臨各種挑戰(zhàn)。由不同沉積環(huán)境和復(fù)雜成巖過程形成的多樣孔隙類型及其對應(yīng)的孔隙結(jié)構(gòu)是導(dǎo)致碳酸鹽巖儲集層強(qiáng)非均質(zhì)性的最主要原因,例如,對中東碳酸鹽巖儲集層的研究發(fā)現(xiàn),在相同孔隙度(25%)情況下,由孔隙結(jié)構(gòu)差異導(dǎo)致的滲透率差異可達(dá)到4個數(shù)量級。此外,由孔隙結(jié)構(gòu)差異引起的地震縱波速度差異達(dá)到了2.5 km/s,甚至更多[9]。因此,孔隙結(jié)構(gòu)是碳酸鹽巖巖石物理研究的核心,有效表征這一參數(shù)的變化以更好地進(jìn)行碳酸鹽巖儲集層孔隙度滲透率表征是國內(nèi)外學(xué)者一直致力解決的問題[10]。另外,流體和巖性也是造成聲波差異的原因之一,但相對而言,孔隙結(jié)構(gòu)是不容忽視的重要因素。據(jù)此,本文綜合巖心、測井和地震資料,充分考慮巖石物理表征過程中干巖石模量計算的復(fù)雜性和碳酸鹽巖的特殊性,選用合適的巖石物理模型,對碳酸鹽巖儲集層的孔隙結(jié)構(gòu)、孔隙度和滲透率進(jìn)行表征。在此基礎(chǔ)上,探討了基于巖石物理的疊前地震儲集層表征的可能性。

1 孔隙結(jié)構(gòu)表征的方法原理

碳酸鹽巖儲集層巖石基質(zhì)礦物主要由方解石、白云石及少量泥質(zhì)組成;孔隙空間主要為孔洞(或稱為硬孔隙或球形孔隙)、粒間孔隙和裂縫,并隨機(jī)分布在巖石骨架中;孔隙流體由油、氣、水組成(見圖1)。以理論巖石物理模型為基礎(chǔ),通過三步法即可建立飽和巖石縱橫波速度與巖石彈性參數(shù)及儲集層參數(shù)的關(guān)系(見圖1)。在這一過程中,孔隙結(jié)構(gòu)類型的表征和干巖石彈性模量的計算仍是亟待解決的問題。

圖1 飛仙關(guān)組碳酸鹽巖儲集層地震巖石物理建模流程

通過對Biot理論模型的延伸,Sun和Goldberg基于裂縫型多孔介質(zhì)的動力學(xué)理論建立了巖石物理模型[11-12]。在這個模型中,反映裂縫型多孔介質(zhì)內(nèi)部構(gòu)造特征的拓?fù)鋮?shù)骨架柔韌性因子(γ)被引入模型,用來定量表征碳酸鹽巖復(fù)雜孔隙結(jié)構(gòu)類型及其復(fù)雜連通關(guān)系對巖石聲波速度的影響。骨架柔韌性因子相對獨立于孔隙度,主要與孔隙結(jié)構(gòu)類型相關(guān),還受孔隙連通性、巖石結(jié)構(gòu)及礦物粒度大小的影響。此外,在骨架柔韌性因子的計算過程中并不需要考慮干巖石的彈性模量,求解過程也更為簡便、可靠,因此被廣泛應(yīng)用于定量刻畫碳酸鹽巖儲集層孔隙結(jié)構(gòu)類型研究中[13-15]。

對于碳酸鹽巖儲集層孔隙結(jié)構(gòu)類型通常采用孔隙縱橫比進(jìn)行表征(見圖2)。球形孔隙相比粒間孔隙或者裂縫連通性差,在外界應(yīng)力作用下,流體無法在三維空間內(nèi)流動,基本是靜止在原孔隙位置。因此,其產(chǎn)生強(qiáng)的反作用力,使巖石相對“硬”。在球形孔隙相對發(fā)育的儲集層中,地震波更容易沿著巖石骨架基質(zhì)傳播,速度快。對于發(fā)育粒間孔隙或裂縫的儲集層,孔隙之間的連通性要比球形孔隙高很多倍。在外界應(yīng)力作用下,孔隙內(nèi)的流體可在三維空間內(nèi)流動,這樣對外界應(yīng)力的反作用力也小,使巖石表現(xiàn)為相對“軟”。好的孔隙連通性使得這類儲集層骨架基質(zhì)的連通性要比球形孔隙發(fā)育的儲集層差很多,所以,地震波在這類儲集層中傳播時必然要經(jīng)過含有流體的孔隙,流體和巖石間的磨擦使地震波傳播能量減弱,進(jìn)而引起地震波傳播速度的降低。物理上,骨架柔韌性因子與孔隙縱橫比成負(fù)相關(guān)關(guān)系[12](見圖2)。為此,同一孔隙度條件下,低骨架柔韌性因子的樣品代表了具有高孔隙縱橫比和聲波傳播速度的球形孔隙發(fā)育的儲集層,而相對高骨架柔韌性因子的樣品則代表了粒間孔隙或裂縫發(fā)育的儲集層。

圖2 碳酸鹽巖孔隙結(jié)構(gòu)類型示意圖

基于巖石物理分析的骨架柔韌性因子計算公式為

γ=1+lnfln (1-φ)

(1)

其中f=1-KfKs+1-KfKsφFk(1-φ)1-KfKsFk

(2)

Fk=Ks-Kφ(Ks-Kf)

(3)

(4)

式中,vp、vs分別為縱、橫波速度,km/s;ρ為體積密度,g/cm3;K為體積模量,GPa;Ks、Kf分別為巖石固體基質(zhì)體積模量和流體體積模量,GPa;φ為孔隙度,%;Fk、f為中間變量。

為準(zhǔn)確估算骨架柔韌性因子(γ),還必須準(zhǔn)確計

算組成巖石礦物的含量、孔隙度、巖石固體基質(zhì)和混合流體的體積模量。研究采用測井復(fù)雜巖性分析方法(CRA)確定組成巖石不同礦物的含量和孔隙度。利用Viogt-Reuss-Hill模型和Wood模型分別對巖石固體基質(zhì)的體積模量和混合流體的體積模量進(jìn)行計算,有

Ks=12∑NiKi+1∑NiKi

(5)

式中,Ni為組成巖石的第i中礦物的體積百分?jǐn)?shù);Ki為對應(yīng)礦物的體積模量。

Kf=(So/Ko+Sg/Kg+Sw/Kw)-1

(6)

式中,Ko、Kg和Kw分別為油、氣和水的體積模量;So、Sg和Sw分別為油、氣和水的飽和度,并且So+Sg+Sw=1。

2 實例分析

普光氣田位于川東北斷褶帶東北段的雙石廟-普光構(gòu)造帶上,是一個北北東向的構(gòu)造-巖性復(fù)合型氣藏,形成于燕山期,定型于喜山期。下三疊統(tǒng)飛仙關(guān)組臺地邊緣相鮞粒灘儲集層是普光氣田最主要的產(chǎn)氣位置,儲集層巖性主要以顆粒白云巖和細(xì)-中晶白云巖為主,受早期沉積作用和后期同生成巖、淺-中埋藏、深埋藏和后期變形4個成巖階段的影響,發(fā)育有粒內(nèi)溶孔,鑄?？?粒間溶孔,晶間孔和裂縫等不同骨架柔韌性因子的儲集空間。

2.1 基于巖石物理分析的速度—孔隙度關(guān)系

普光氣田白云巖儲集層白云石含量大部分在90%以上,方解石4%～8%,含極少量泥質(zhì);流體類型為高含H2S和CO2的天然氣。表1是利用CRA方法進(jìn)行計算時所用的參數(shù)值,圖3為計算得到的A井取心段巖石不同礦物組成的含量和孔隙度值。對比巖心分析孔隙度和計算得到的孔隙度可以發(fā)現(xiàn)兩者具有較高的一致性,說明了CRA方法的適用性。取心井段的巖心薄片資料展示了垂向上孔隙結(jié)構(gòu)的變化,5 602.34 m和5 642.24 m處的粒內(nèi)溶孔較其他位置的儲集層具有較大的孔隙度和較高的速度。研究區(qū)裂縫主要發(fā)育于飛仙關(guān)組中上部的非有效灰?guī)r儲集層中,為此,研究中剔除了裂縫對骨架柔韌性因子的影響。

表1 礦物含量和孔隙度計算時所用的參數(shù)值

圖3 A井取心段測井曲線組合及其巖石物理特征

圖4 A井取心段密度孔隙度與縱波速度交會圖

圖4為A井取心段巖心對應(yīng)密度孔隙度和縱波速度之間的交會圖,色標(biāo)指示了骨架柔韌性因子的變化。從圖4中可以看到,縱波速度和孔隙度具有明顯的負(fù)相關(guān)關(guān)系。在相同的孔隙度條件下,由于孔隙結(jié)構(gòu)的差異,速度具有很強(qiáng)的發(fā)散性。當(dāng)孔隙度等于12.5%時,速度差異可以達(dá)到1 000 m/s。根據(jù)骨架柔韌性因子的差別,可以將縱波速度—孔隙度關(guān)系劃分為3個不同的趨勢區(qū)域:橙色趨勢區(qū)域的骨架柔韌性因子小于4,藍(lán)色趨勢區(qū)域的骨架柔韌性因子介于4和6之間,綠色趨勢區(qū)域的骨架柔韌性因子大于6。高孔隙度低速度的樣本點表現(xiàn)為低的骨架柔韌性因子,而低孔隙度高速度的樣本點則對應(yīng)于高的骨架柔韌性因子。

選取圖4中10個鑄體薄片對應(yīng)的樣本點,A1～A4、B1～B3、C1～C3分別位于不同趨勢區(qū)域中,其中A3和B1、A4和B2、B3和C1分別對應(yīng)于孔隙度5.5%、7%和14%。對比樣本點B3和C1可以發(fā)現(xiàn)兩者在相同孔隙度條件下,速度差異達(dá)到了400 m/s。進(jìn)一步分析發(fā)現(xiàn),B3和C1具有相同的巖性,但巖石結(jié)構(gòu)差異明顯。B3是結(jié)晶白云巖,而C1是顆粒白云巖。此外,B3具有晶間孔和部分粒內(nèi)溶孔混合孔隙類型,C1則發(fā)育相對簡單的粒內(nèi)溶孔或鑄模孔(見圖4)。樣本點A3、B1、A4和B2都是結(jié)晶灰?guī)r,但是B1和B2具有更多的粒內(nèi)溶孔。因此,在相同的孔隙度條件下,B1和B2具有相對強(qiáng)的巖石剛度和大的縱波傳播速度。樣本點A3和A4以晶間孔為主,相比其他孔隙結(jié)構(gòu)類型的儲集層而言,縱波在這類儲集層中傳播時,由于流體和巖石之間的作用會引起速度顯著的降低(見圖4)。

圖5 B井飛仙關(guān)組密度孔隙度與縱波速度交會圖

圖5為B井飛仙關(guān)組全井段密度孔隙度與縱波速度的交會圖。對比發(fā)現(xiàn),圖5和圖4具有類似的趨勢,發(fā)育粒內(nèi)溶孔或鑄?？椎臉颖军c沿黑色趨勢線分布,如圖5中C1～C3,其骨架柔韌性因子小于5;發(fā)育粒內(nèi)溶孔、晶間溶孔混合孔隙類型的樣本點沿著紅色趨勢線分布,如圖5中的B1～B3,其骨架柔韌性因子介于5～8之間;發(fā)育晶間溶孔和微孔隙的樣本點沿著綠色的趨勢線分布,如圖5中的A1～A4,其骨架柔韌性因子大于8。骨架柔韌性因子將單一的速度—孔隙度關(guān)系劃分為具有不同孔隙結(jié)構(gòu)特征的3個區(qū)域,從而利用3個速度—孔隙度關(guān)系表征具有不同孔隙結(jié)構(gòu)特征的儲集層。該方法充分考慮了儲集層孔隙結(jié)構(gòu)變化對速度的影響,充分利用數(shù)據(jù)點,用3個高精度的速度—孔隙度關(guān)系表征儲集層,提高儲集層表征的精度。此外,該方法還能有效避免忽略掉儲集層表征過程中高孔隙度低速度端鑄模孔或粒內(nèi)溶孔發(fā)育的有效儲集層和低孔隙度高速度端晶間溶孔發(fā)育的有效儲集層。

圖6 普光氣田飛仙關(guān)組儲集層孔隙度與滲透率交會圖*非法定計量單位,1 mD=9.87×10-4 μm2,下同

2.2 基于巖石物理分析的孔隙度—滲透率關(guān)系

圖7 不同孔隙結(jié)構(gòu)儲集層壓汞曲線特征

圖8 骨架柔韌性因子與流動帶指數(shù)交會圖

利用普光氣田取心段實測的覆壓孔隙度滲透率資料,建立研究區(qū)孔隙度和滲透率之間的關(guān)系(見圖6)。由復(fù)雜孔隙結(jié)構(gòu)特征導(dǎo)致的碳酸鹽巖儲集層強(qiáng)的非均質(zhì)性使得研究區(qū)滲透率離散分布。利用骨架柔韌性因子,可以將研究區(qū)離散的孔隙度和滲透率關(guān)系劃分為3種不同的趨勢(見圖6)。橙色樣本點為粒內(nèi)溶孔或鑄模孔發(fā)育的儲集層(見圖6中C1),其骨架柔韌性因子小于4。該類儲集層孔隙結(jié)構(gòu)以大孔微喉、中孔微喉為主,孔隙之間連通性差。圖7中①、②毛細(xì)管壓力曲線表現(xiàn)為細(xì)歪度,分選中等,儲集性較好。藍(lán)色樣本點為晶間溶孔和粒內(nèi)溶孔混合孔隙發(fā)育的儲集層(見圖6中B2、B3),其骨架柔韌性因子介于4～6之間。該類儲集層孔隙結(jié)構(gòu)以大孔中喉、中孔細(xì)喉為主,孔隙之間連通性相對較好,主要取決于晶間孔隙的發(fā)育程度。圖7中③、④毛細(xì)管壓力曲線表現(xiàn)為中-細(xì)歪度,雙平臺的特征,分選中等,儲集性較好。綠色樣本點為晶間溶孔發(fā)育的儲集層(見圖6中A4),其骨架柔韌性因子大于6。該類儲集層孔隙結(jié)構(gòu)以中孔中喉、中孔細(xì)喉為主,孔隙之間連通性最好。圖7中⑤、⑥毛細(xì)管壓力曲線表現(xiàn)為粗-中歪度,分選較好,儲集性好。在相同孔隙度條件下,骨架柔韌性因子大于6的儲集層由于孔喉連通性最好而使其滲透率要比其他兩類儲集層高幾個數(shù)量級。據(jù)圖6可知,研究區(qū)內(nèi)滲透率最高的儲集層是晶間溶孔最發(fā)育的儲集層而不是孔隙度最大的儲集層。根據(jù)骨架柔韌性因子,研究區(qū)存在3類孔隙度和滲透率關(guān)系,相比圖6中紅色曲線所示的單一孔隙度和滲透率關(guān)系,其相關(guān)系數(shù)明顯提高,儲集層滲透率的表征精度也隨之得到顯著提高。

2.3 流動帶指數(shù)對比分析

巖石物理相為沉積、成巖和構(gòu)造的綜合效應(yīng),決定了碳酸鹽巖儲集層中孔隙、溶孔及裂縫的發(fā)育程度和共生關(guān)系。流動帶指數(shù)是進(jìn)行儲集層巖石物理相劃分的重要指標(biāo),是表征儲集層微觀孔隙結(jié)構(gòu)、反映宏觀滲流特征的重要參數(shù)[16]。在地震巖石物理研究過程中,骨架柔韌性因子是分析地震波在不同孔隙結(jié)構(gòu)類型儲集層中傳播規(guī)律,建立巖性參數(shù)、物性參數(shù)和地震速度、密度等彈性參數(shù)之間關(guān)系的橋梁。基于骨架柔韌性因子的巖石物理表征和基于流動帶指數(shù)的巖石物理相分類都體現(xiàn)了儲集層巖性、物性和孔隙結(jié)構(gòu)對儲集層質(zhì)量的控制作用,兩者之間必然存在內(nèi)在關(guān)系。圖8為A井取心段骨架柔韌性因子和流動帶指數(shù)的交會圖,圖8中可以看到兩者之間具有很好的正相關(guān)關(guān)系。在碳酸鹽巖儲集層中,流動帶指數(shù)高值區(qū)通常代表具有好的孔喉連通性、晶間溶孔相對發(fā)育的儲集層,其對應(yīng)于高γ值。流動帶指數(shù)低值區(qū)代表了孔喉連通性差、粒內(nèi)溶孔或鑄模孔發(fā)育的儲集層,其對應(yīng)于低γ值,在這類儲集層中,縱波以相對較快的速度沿骨架傳播。混合孔隙發(fā)育的儲集層,孔喉連通性介于以上2類儲集層之間,故對應(yīng)于中等的流動帶指數(shù)和γ值。骨架柔韌性因子和流動帶指標(biāo)之間高的相關(guān)性從巖石物理相(沉積相、成巖相和裂縫相)的角度進(jìn)一步驗證了骨架柔韌性因子在具有強(qiáng)非均質(zhì)性的碳酸鹽巖含氣儲集層孔滲表征過程中的適用性。

圖9 3種不同孔隙結(jié)構(gòu)類型碳酸鹽巖儲集層反射系數(shù)隨入射角變化圖

2.4 成因機(jī)制簡析

普光氣田不同孔隙結(jié)構(gòu)儲集層的發(fā)育在宏觀上受沉積相的控制,微觀上受成巖作用的控制。臺地邊緣鮞粒灘亞相是有利儲集層集中發(fā)育位置,其他相帶儲集條件差。粒內(nèi)溶孔或鑄模孔發(fā)育的儲集層主要發(fā)育于鮞粒灘灘核的位置,同生成巖階段的大氣淡水選擇性溶蝕、膠結(jié)和混合白云巖化作用一方面使儲集層發(fā)育豐富的粒內(nèi)溶孔、鑄模孔;另一方面形成強(qiáng)烈的粒間白云石膠結(jié)使形成的孔隙得到很好的保存。后期成巖改造對該類儲集層影響較弱,孔隙多呈孤立狀,連通性較差,使該類儲集層總體上表現(xiàn)為相對高孔隙度、低滲透率的特征。晶間溶孔發(fā)育的儲集層主要發(fā)育于鮞粒灘灘緣的位置,早期成巖作用較弱,受后期埋藏白云巖化和重結(jié)晶作用的影響,儲集層原始顆粒結(jié)構(gòu)基本被破壞,重結(jié)晶成細(xì)-中晶白云巖,局部可見殘余鮞粒結(jié)構(gòu)。白云巖儲集層受后期非選擇性溶蝕作用的影響,晶間孔繼續(xù)擴(kuò)大形成孔喉連通較好的晶間溶孔,使該類儲集層總體上表現(xiàn)為相對低孔隙度、高滲透率的特征。混合孔隙結(jié)構(gòu)類型發(fā)育的儲集層主要形成于灘核-灘緣的過渡位置,受整個成巖階段的影響,儲集層形成孔喉結(jié)構(gòu)復(fù)雜的混合孔隙組合類型,孔隙度滲透率分布介于上述兩類儲集層之間[17-18]。

3 疊前地震儲層表征的探討

碳酸鹽巖儲集層的空間表征是在儲集層測井評價的基礎(chǔ)上,基于地震資料進(jìn)行的。相比疊后地震資料,疊前地震數(shù)據(jù)在碳酸鹽巖孔滲表征過程中可以提供更多、更準(zhǔn)確的關(guān)于孔隙結(jié)構(gòu)特征的信息,從而更好地預(yù)測碳酸鹽巖儲集層的非均質(zhì)性[19-21]。大量理論研究表明,不同孔隙結(jié)構(gòu)類型的碳酸鹽巖儲集層在給定孔隙度條件下具有不同的AVO響應(yīng)特征。為此,研究中構(gòu)建了以飛仙關(guān)組儲蓋組合為原型的雙層介質(zhì)模型,以探討疊前地震儲層表征的可行性。模型參數(shù)來自于A井中具有相似孔隙度、不同孔隙結(jié)構(gòu)類型的3個樣本點(見表2)。通過求解Zoeppritz方程可以得到具有不同孔隙結(jié)構(gòu)類型的3類儲集層縱、橫波反射系數(shù)隨入射角變化而變化情況。圖9中,不同孔隙結(jié)構(gòu)類型的儲集層具有不同的臨界角,并且從粒內(nèi)溶孔到晶間溶孔儲集層,縱波和轉(zhuǎn)換橫波的臨界角角度逐漸增加。就縱波反射系數(shù)而言,不同孔隙結(jié)構(gòu)類型儲集層反射系數(shù)的差異隨著入射角的增大而增大。例如,在入射角小于30°的情況下反射系數(shù)的差值近似等于0.05,當(dāng)入射角大于30°時反射系數(shù)的差異迅速增大,不同孔隙結(jié)構(gòu)類型儲集層的地震響應(yīng)特征差異越來越明顯。對于轉(zhuǎn)換橫波的反射系數(shù),其具有和縱波反射系數(shù)類似的特征,但反射系數(shù)發(fā)生明顯變化對應(yīng)的入射角要大于40°。上述AVO模擬結(jié)果表明,孔隙結(jié)構(gòu)類型的變化的確對地震響應(yīng)特征有重要的影響,相比小偏移距地震數(shù)據(jù)而言,大偏移距地震數(shù)據(jù)能更好地區(qū)分碳酸鹽巖儲集層中孔隙結(jié)構(gòu)類型的變化以便更好地進(jìn)行儲集層孔滲表征。以上結(jié)果進(jìn)一步印證了普光氣田飛仙關(guān)組儲集層疊前地震表征的可能性,利用疊前地震數(shù)據(jù)體,結(jié)合骨架柔韌性因子的分析,可以按不同孔隙結(jié)構(gòu)類型進(jìn)行井間孔隙度滲透率的表征。

表2 用于正演的不同孔隙結(jié)構(gòu)類型碳酸鹽巖儲集層參數(shù)

4 結(jié) 論

(1) 骨架柔韌性因子是用來定量表征碳酸鹽巖儲集層復(fù)雜孔隙結(jié)構(gòu)類型對巖石聲波速度影響的一個參數(shù)。在相同孔隙度條件下,低骨架柔韌性因子代表了具有高孔隙縱橫比和聲波傳播速度的球形孔隙發(fā)育的儲集層,而相對高骨架柔韌性因子則代表了粒間孔隙或裂縫發(fā)育的儲集層。在骨架柔韌性因子的計算過程中不需要考慮干巖石的彈性模量,求解過程更為簡便、可靠。

(2) 骨架柔韌性因子可以量化不同的孔隙結(jié)構(gòu),從而建立骨架柔韌性因子約束下的多個速度—孔隙度—滲透率關(guān)系式表征具有不同孔隙結(jié)構(gòu)特征儲集層的孔隙度滲透率特征,提高儲集層表征的精度。這一方法避免了單一關(guān)系模型未考慮儲集層孔隙結(jié)構(gòu)變化對速度影響和精度不足的問題。

(3) 沉積和成巖控制了碳酸鹽巖儲集層不同孔隙結(jié)構(gòu)的發(fā)育,基于骨架柔韌性因子的巖石物理表征和基于流動帶指標(biāo)參數(shù)的巖石物理相分類之間高的相關(guān)性進(jìn)一步證明了骨架柔韌性因子在具有強(qiáng)非均質(zhì)性的碳酸鹽巖含氣儲集層孔隙度滲透率表征過程中的適用性。

(4) 基于AVO正演模擬的結(jié)果表明,大偏移距的疊前地震數(shù)據(jù)具有評價研究區(qū)碳酸鹽巖儲集層復(fù)雜孔隙結(jié)構(gòu)和進(jìn)行井間孔隙度滲透率表征的可能性。

參考文獻(xiàn):

[1] 鄧?yán)^新, 王歡, 周浩, 等. 龍馬溪組頁巖微觀結(jié)構(gòu)、地震巖石物理特征與建模 [J]. 地球物理學(xué)報, 2015(6): 2123-2136.

[2] 張廣智, 陳懷震, 王琪, 等. 基于碳酸鹽巖裂縫巖石物理模型的橫波速度和各向異性參數(shù)預(yù)測 [J]. 地球物理學(xué)報, 2013, 56(5): 1707-1715.

[3] 張廣智, 陳嬌嬌, 陳懷震, 等. 基于巖石物理模版的碳酸鹽巖含氣儲層定量解釋 [J]. 吉林大學(xué)學(xué)報: 地球科學(xué)版, 2015, 45(2): 630-638.

[4] Li Jingye, Chen Xiaohong. A Rock-phsical Modeling Method for Carbonate Reservoirs at Seismic Scale [J]. Applied Geophsics, 2013, 10(1): 1-13.

[5] 劉欣欣, 印興耀, 張峰. 一種碳酸鹽巖儲層橫波速度估算方法 [J]. 中國石油大學(xué)學(xué)報: 自然科學(xué)版, 2013, 37(1): 42-49.

[6] 熊曉軍, 呂龑, 林凱, 等. 一種碳酸鹽巖流體替換新方法 [J]. 中國石油大學(xué)學(xué)報: 自然科學(xué)版, 2012, 36(6): 57-61.

[7] 鄔龍, 劉堂晏, 王紅濤, 等. 干巖石模量的計算方法對比分析及其應(yīng)用 [J]. 地球物理學(xué)進(jìn)展, 2014, 29(5): 2272-2277.

[8] 伍向陽, 陳祖安. 干巖石體積模量估計 [J]. 地球物理學(xué)報, 2001, 44: 146-151.

[9] Sun Y F, Berteussen K, Vega S, et al. Effects of Pore Structure on 4D Seismic Signals in Carbonate Reservoirs [J]. SEG Technical Program Expanded Abstracts, 2006(1): 3260.

[10] 姜均偉, 朱宇清, 徐星, 等. 伊拉克H油田碳酸鹽巖儲層的孔隙結(jié)構(gòu)特征及其對電阻的影響 [J]. 地球物理學(xué)進(jìn)展, 2015, 30(1): 203-209.

[11] Sun Y F, Goldberg D. Effects of Aspect Ratio on Wave Velocities in Fractured Rocks [J]. SEG Expanded Abstracts, 1997(1): 925.

[12] Sun Y F. Pore Structure Effects on Elastic Wave Propagation in Rocks: AVO Modelling [J]. Journal of Geophysics & Engineering, 2004, 1(4): 268-276.

[13] Gartner G L B, Wagner P D, Baechle G T, et al. Obtaining Permeability from Seismic Data-A New Breakthrough in Carbonate Reservoir Modelling [J]. International Petroleum Technology Conference, 2005.

[14] Shedid S A, Almehaideb R A. Enhanced Reservoir Description of Heterogeneous Carbonate Reservoirs [J]. Journal of Canadian Petroleum Technology, 2003, 42(7): 23-28.

[15] Weger R J, Eberli G P, Baechle G T, et al. Quantification of Pore Structure and its Effect on Sonic Velocity and Permeability in Carbonates [J]. AAPG Bulletin, 2009, 10(10): 1297-1317.

[16] 王起琮, 趙淑萍, 魏欽廉, 等. 鄂爾多斯盆地中奧陶統(tǒng)馬家溝組海相碳酸鹽巖儲集層特征 [J]. 古地理學(xué)報, 2012, 14(2): 229-242.

[17] 馬永生, 郭彤樓, 趙雪鳳, 等. 普光氣田深部優(yōu)質(zhì)白云巖儲層形成機(jī)制 [J]. 中國科學(xué) D輯: 地球科學(xué), 2007, 37: 43-52.

[18] 王恕一, 蔣小瓊, 管宏林, 等. 川東北普光氣田鮞粒白云巖儲層粒內(nèi)溶孔的成因 [J]. 沉積學(xué)報, 2010, 28(1): 10-16.

[19] 龔洪林, 劉偉方, 張喜梅. 塔北輪古東碳酸鹽巖儲層巖石物理響應(yīng)特征 [J]. 地球物理學(xué)進(jìn)展, 2014, 29(1): 331-338.

[20] 孟憲軍, 劉福平, 王玉梅, 等. 地震反射系數(shù)的相角變化 [J]. 石油地球物理勘探, 2010: 121-124.

[21] 彭通, 徐平. 縱波反射系數(shù)影響因素分析 [J]. 現(xiàn)代礦業(yè), 2009(8): 67-70.