史全黨,孔令業(yè),吳超,丁艷雪,劉澤民,于雪,王江

(1. 中國石油新疆油田公司采氣一廠,新疆 克拉瑪依 834000;2. 中國石油大慶油田有限責(zé)任公司 勘探開發(fā)研究院,黑龍江 大慶 163712)

0 引言

隨著含油氣盆地精細(xì)構(gòu)造解釋、精細(xì)儲層預(yù)測和精細(xì)油藏描述工作的深入,對地震資料分辨率的需求越來越高[1-4],精細(xì)儲層預(yù)測的手段也一直以地震波阻抗反演為主[5-8]。根據(jù)目標(biāo)區(qū)的構(gòu)造背景、沉積環(huán)境制定合理的技術(shù)對策和利用波阻抗反演方法來減小波阻抗反演結(jié)果的多解性、提高儲層預(yù)測能力一直是地震波阻抗反演和儲層精細(xì)預(yù)測的攻關(guān)方向[9-10]。姚逢昌等[11]通過對基于模型的波阻抗反演方法進(jìn)行分析,認(rèn)為多解性是基于測井約束性框架模型波阻抗反演方法所固有的特性,取決于初始模型和實際地下構(gòu)造特點(diǎn)、沉積特征的符合程度。現(xiàn)有的測井約束波阻抗反演方法由于各自不同的方法原理,均存在一定局限性[12-13],地震高頻分量主要依賴于初始波阻抗模型從測井曲線中獲取,初始建模主要是在解釋層位控制下的井間插值,沒有考慮斷層在初始建模中的作用,橫向上也沒有充分利用地震信息。當(dāng)目標(biāo)區(qū)構(gòu)造條件復(fù)雜、斷層發(fā)育、沉積環(huán)境和巖性橫向變化較大時,利用鉆井?dāng)?shù)據(jù)井間插值或外推得到的初始模型低頻背景較好,而反映巖性橫向井間局部變化的高頻成分獲取較難,距離插值井點(diǎn)越遠(yuǎn),產(chǎn)生的誤差越大,減小了波阻抗反演方法對初始模型高頻成分的應(yīng)用程度,降低了波阻抗反演儲層預(yù)測精度。

基于小波邊緣分析與建模的波阻抗反演(acoustic impedance inversion based on wavelet edge analysis and modelling)技術(shù)(簡稱AIW),利用小波邊緣分析方法從地震記錄中直接提取地震屬性特征參數(shù),在斷層參與的高密度層控構(gòu)造框架約束下,同測井聲波阻抗數(shù)據(jù)一起建立初始模型,相互約束、相互補(bǔ)充,充分利用地震數(shù)據(jù)橫向分辨率高、縱向控制層位密度大的優(yōu)勢,避免了常規(guī)波阻抗反演過程中初始模型建立不準(zhǔn)而產(chǎn)生的地質(zhì)影響。許崇寶等[14]、崔永福等[15]、謝裕江等[16]應(yīng)用AIW波阻抗反演技術(shù)對準(zhǔn)東煤田和新疆紅沙泉區(qū)塊侏羅系西山窯組含煤地層的煤層厚度、塔里木盆地英買力地區(qū)志留系瀝青砂巖以及中國MOU氣田盒8段河流相心灘砂巖進(jìn)行了預(yù)測,彌補(bǔ)了常規(guī)波阻抗反演建模過程中井間高頻成分缺失,提高了儲層反演精度,但應(yīng)用的區(qū)塊大多構(gòu)造簡單、地層橫向變化不大,很少針對構(gòu)造復(fù)雜、斷裂發(fā)育、沉積環(huán)境多變的斷陷盆地以及低分辨率、低信噪比目標(biāo)區(qū)進(jìn)行研究。

本文應(yīng)用小波邊緣分析與井—震聯(lián)合建模的波阻抗反演技術(shù),考慮準(zhǔn)噶爾盆地陸梁隆起帶3D地震的波組反射特性以及地下地震波場的分布特點(diǎn),通過小波邊緣分析提取反映巖性局部變化的地震特征參數(shù),在斷層參與的高密度層控構(gòu)造框架約束下,參與初始模型的建立和波阻抗的反演迭代,確定波阻抗反演的關(guān)鍵環(huán)節(jié)、優(yōu)化克拉美麗氣田DX14井區(qū)儲層預(yù)測流程,不但精細(xì)刻畫了克拉美麗氣田DX14井區(qū)扇三角洲前緣砂體的分布,也提高了有利儲層的預(yù)測精度。

1 區(qū)域地質(zhì)概況

克拉美麗氣田位于準(zhǔn)噶爾盆地陸梁隆起東南部的滴南凸起中段,滴南凸起整體為一大型西傾的鼻狀構(gòu)造,受準(zhǔn)噶爾盆地多期的構(gòu)造建造和改造影響,斷裂發(fā)育、構(gòu)造復(fù)雜。滴南凸起中段是準(zhǔn)噶爾盆地腹部最早展開勘探的地區(qū)之一,至今已發(fā)現(xiàn)了石炭系、二疊系、侏羅系和白堊系多個油氣藏,其中二疊系梧桐溝組油氣藏埋深3 500～4 000 m,而且構(gòu)造破碎、斷裂發(fā)育,多發(fā)育為斷層控制的小斷塊和小斷鼻構(gòu)造,氣藏類型主要為構(gòu)造—巖性氣藏。主要含油氣層為梧桐溝組一段,巖性為細(xì)砂巖和粉砂巖,儲集空間類型主要為剩余粒間孔及粒內(nèi)溶孔。根據(jù)電性及儲層特征,自下而上可劃分為兩段(P3wt2、P3wt1),油層主要發(fā)育在P3wt1,砂巖儲層不但低孔、特低滲、非均質(zhì)性強(qiáng),而且儲層橫向連續(xù)性差,縱向變化大,使得滴南凸起中段DX14井區(qū)的地震資料不但構(gòu)造成像效果差、分辨率低,而且砂體的地震響應(yīng)特征不明顯。

通過滴南凸起中段DX14井區(qū)原始地震資料分析,在目的層段梧桐溝組存在弱—空白以及低分辨率的地震反射,地震資料主頻為30 Hz,如果地層平均速度按3 900 m/s計算,滴南凸起中段DX14井區(qū)3D地震資料可識別砂體的最小厚度為32.5 m,而滴西DX14井區(qū)含氣砂體厚度一般為2.0～15 m,所以常規(guī)地震資料的分辨能力是無法滿足精細(xì)構(gòu)造解釋和精細(xì)儲層預(yù)測的地質(zhì)需求,因此,應(yīng)在精細(xì)井—震一體化初始模型建立的基礎(chǔ)上開展高精度的波阻抗反演和精細(xì)儲層預(yù)測。

2 技術(shù)方法

2.1 多尺度小波邊緣分析與地震特征參數(shù)建模

小波邊緣分析是一個理想的地震特征參數(shù)提取方法,利用小波邊緣分析從地震數(shù)據(jù)本身提取地震特征參數(shù)得到地震數(shù)據(jù)局部的變化信息,通過地震數(shù)據(jù)的小波變換函數(shù)的一階微分極大值和二階微分零值進(jìn)行處理,對處理后的梯度點(diǎn)進(jìn)行判斷,確定其為邊緣點(diǎn)。

對地震數(shù)據(jù)d進(jìn)行小波變換得到Ws(d),計算公式為

(1)

對式(1)兩側(cè)求偏導(dǎo),且令

(2)

(3)

式中:d為地震數(shù)據(jù);Ws(d)為地震數(shù)據(jù)的小波變換;s為小波尺度;x為平移量,s;ψ為母小波函數(shù);t 為時間,s;b為時間延遲x秒的地震數(shù)據(jù),小波尺度s控制小波函數(shù)的伸縮,平移量x控制小波函數(shù)的平移。

聯(lián)合式(1)、式(2)、式(3),求解此方程組,可得到地震記錄中不同尺度、不同級別的地震特征參數(shù)邊緣點(diǎn)(斷層和巖性體的邊緣)。

這樣,通過小波變換把地震信號分解為多個不同尺度的分量,利用其良好的局部時頻分辨力提取不同精度、不同奇異度的邊緣點(diǎn),最后連接起來就形成綜合邊緣。提取的地震特征參數(shù)的邊緣點(diǎn)不但反映了原始剖面中不易識別的地質(zhì)特征和地下巖性縱、橫向的局部變化,將綜合邊緣點(diǎn)控制的地震特征參數(shù)同測井曲線一起參與波阻抗反演建模,克服了井間插值反演建模過程中無法充分利用地震信息的不足,獲取了反映巖性橫向井間局部變化的高頻成分,從而也補(bǔ)充了波阻抗反演過程中高、低頻地震分量。

2.2 AIW波阻抗反演算法

AIW波阻抗反演采用非線性全局優(yōu)化的非?？焖倌M退火法(VFSA)作為基本反演方法,該算法通過模型反復(fù)擾動修改不但可以避免線性反演算法強(qiáng)烈依賴初始模型而落入局部極值的弊端,且可以高精度求得反演問題的全局最優(yōu)解[17],使得AIW波阻抗反演結(jié)果與地下地質(zhì)情況更加符合,其算法如下:

1)將歸一化相似系數(shù)定義為AIW波阻抗反演的目標(biāo)函數(shù)E[18],來計算初始框架模型g0的目標(biāo)函數(shù)值E(g0),計算公式為

2)聯(lián)合地震特征參數(shù)和測井?dāng)?shù)據(jù)一起來約束模型擾動中的模型變量的變化區(qū)間[Ai,Bi],修改初始模型g0獲取新模型g,并根據(jù)關(guān)系式(4)計算新模型g的目標(biāo)函數(shù)值E(g)以及新模型目標(biāo)函數(shù)E(g)與當(dāng)前模型目標(biāo)函數(shù)E(g0)差ΔE,即ΔE=E(g)-E(g0)。

3)當(dāng)ΔE<0時,確定新模型g為最終的初始模型,如果確定了新模型g,置g0=g,E(g0)=E(g)。

4)當(dāng)ΔE>0,則按依賴于溫度的似Cauchy分布產(chǎn)生新模型,計算公式為:

g′i=gi+yi(Bi-Ai),

(5)

yi=Tsgn(u-0.5)[(1+1/T)|2u-1|-1] ,

(6)

式中:g′i為修改后的第i模型,g′i∈[Ai,Bi];gi為當(dāng)前的第i模型,gi∈[Ai,Bi];yi為擾動因子;T為溫度,K;sgn為符號函數(shù);u為均勻分布的隨機(jī)量,且u∈[0,1]。

5)同一溫度下,多次重復(fù)步驟3)和4),緩慢降低溫度,重復(fù)步驟3)～5),直至收斂條件ΔE<0,反復(fù)迭代反演,獲取高精度的AIW波阻抗反演數(shù)據(jù)體。

2.3 AIW波阻抗反演的關(guān)鍵環(huán)節(jié)

AIW波阻抗反演的關(guān)鍵在于利用小波邊緣分析方法從地震記錄中直接提取地震屬性特征參數(shù),由于提取的地震屬性特征參數(shù)不但反映地下巖性橫向局部變化,而且直接參與了波阻抗反演過程中初始模型的建立和反復(fù)迭代反演,彌補(bǔ)了井間插值建模過程中井間高頻成分的缺失和井間巖性的局部變化,避免了常規(guī)波阻抗反演過程中初始模型建立不準(zhǔn)確而產(chǎn)生的影響。主要關(guān)鍵環(huán)節(jié)為:

1)利用測井和地震資料開展井震精細(xì)對比,確定不同波組的地質(zhì)表征,通過縱向細(xì)劃地層對比單元、橫向精細(xì)對比,對三維地震資料開展精細(xì)的構(gòu)造解釋,得到精細(xì)的構(gòu)造和沉積格架,聯(lián)合測井曲線建立波阻抗低頻模型;

2)通過小波邊緣分析提取反映地下巖性橫向局部變化的地震屬性特征參數(shù),結(jié)合精細(xì)標(biāo)定的測井曲線以及應(yīng)用解釋層位、解釋斷層和測井曲線建立的波阻抗低頻模型重新建立新的AIW波阻抗模型,進(jìn)一步獲取波阻抗高、低頻成分;

3)提取的地震屬性特征參數(shù)同鉆井聲波測井曲線中速度或聲波阻抗數(shù)據(jù)一起進(jìn)行模型擾動和修改,利用非線性全局優(yōu)化的非?？焖倌M退火法,反復(fù)迭代反演,獲取高精度的AIW波阻抗反演數(shù)據(jù)體。

3 應(yīng)用實例

3.1 井—震精細(xì)標(biāo)定與儲層敏感性分析

圖1 波阻抗重構(gòu)前(a)、后(b)的電阻率與波阻抗交匯Fig.1 Intersection of resistivity and wave impedance before(a) and after(b) wave impedance reconstruction

3.2 構(gòu)造框架模型

以往波阻抗反演的構(gòu)造框架模型建立主要是依賴于地震解釋的區(qū)域性組、段級標(biāo)準(zhǔn)反射層位,而很少應(yīng)用一些層間的砂組和油層反射層位,同時也完全不考慮斷層對反演結(jié)果的影響,在構(gòu)造簡單、地層橫向變化不大的凹陷性盆地這種影響較小,但在構(gòu)造復(fù)雜、斷層發(fā)育、地層橫向變化較大的斷陷盆地就會產(chǎn)生相當(dāng)大的影響:

1)由于斷陷盆地地層橫向變化,在大的區(qū)域性標(biāo)準(zhǔn)反射層控下內(nèi)插的井間測井信息就會發(fā)生構(gòu)造誤差,為了減小內(nèi)插測井信息的構(gòu)造誤差,開展層序地層研究,縱向細(xì)化地層單元,開展砂組級層序?qū)Ρ?以砂組級的地震反射層面來實現(xiàn)層控構(gòu)造框架建模,減小區(qū)域性標(biāo)準(zhǔn)反射層控下井間內(nèi)插測井信息發(fā)生的構(gòu)造誤差;

2)由于斷陷盆地復(fù)雜的斷裂結(jié)構(gòu),井間內(nèi)插時也沒有考慮斷層兩側(cè)內(nèi)插測井?dāng)?shù)據(jù)的突變,不但影響斷層兩側(cè)的巖性接觸關(guān)系,降低斷層兩側(cè)波阻抗反演精度,而且降低了波阻抗反演剖面對斷層的識別能力,為了解決這一問題,將斷層轉(zhuǎn)化為地震層位,明確斷層與縱向細(xì)化地層單元間的接觸關(guān)系,同地震反射層位一起控制測井信息的井間內(nèi)插,克服了斷層兩側(cè)由于井間測井?dāng)?shù)據(jù)內(nèi)插而引起的測井?dāng)?shù)據(jù)、地震屬性特征參數(shù)誤差,實現(xiàn)了斷層兩側(cè)內(nèi)插測井?dāng)?shù)據(jù)以及地震屬性特征參數(shù)突變。

3.3 小波邊緣分析

小波邊緣分析具有多尺度分辨的特點(diǎn),可以有效地區(qū)分地震剖面中不同尺度、不同級別的地震特征參數(shù)邊緣點(diǎn),利用小波邊緣分析獲取的地震特征參數(shù)中的邊緣點(diǎn)可以是斷層、巖性體的邊緣。

圖2和圖3分別為DX14井區(qū)三維L1036地震剖面以及經(jīng)小波邊緣分析得到的地震特征參數(shù)邊緣檢測特征點(diǎn)剖面。從圖3可見,經(jīng)過小波多尺度邊緣分析處理后,縱向地層有效分開,相鄰道的不連續(xù)性明顯,微斷裂的斷點(diǎn)被突出、地層微小錯動特征清晰,反映的地層接觸關(guān)系更準(zhǔn)確,提取的地震特征參數(shù)的邊緣點(diǎn)可以反映原始剖面中不易識別的地質(zhì)特征,將邊緣點(diǎn)控制的地震特征參數(shù)同測井曲線一起參與波阻抗反演建模,不但補(bǔ)充了波阻抗高、低頻成分,而且把地層橫向間的接觸關(guān)系完全考慮在波阻抗反演的整個過程中,并且在迭代過程中應(yīng)用此參數(shù)對其進(jìn)行擾動限制,使AIW波阻抗反演模型反映的地下地質(zhì)情況更加真實,減少了波阻抗反演的多解性。

圖2 DX14井區(qū)三維L1036地震剖面Fig.2 L1036 seismic profile in DX14 well area

3.4 井—震一體化初始模型建立

在波阻抗反演框架模型的建立過程中,低頻分量主要來源于兩種途徑:一是將測井?dāng)?shù)據(jù)中的高頻成分濾去;二是直接應(yīng)用地震處理的速度譜資料,目前大多數(shù)波阻抗反演的低頻成分主要是濾去聲波測井曲線中的高頻成分獲得。

同常規(guī)波阻抗反演不同,DX14井區(qū)AIW波阻抗反演初始模型建立過程中突出井—震一體化、構(gòu)造—巖性統(tǒng)一:

1)應(yīng)用鉆井?dāng)M聲波測井曲線獲得速度或者阻抗的低頻成分,然后在斷層參與的高密度層控構(gòu)造框架控制下橫向沿著縱向細(xì)化的層序界面外推。為了縱向細(xì)化低頻模型,保證低頻模型井間構(gòu)造精度,縱向選擇梧桐溝組一段3個砂組和梧桐溝組二段1個砂組的底界以及百口泉組底界5個地震反射層做為主要的控制層,當(dāng)遇到斷層時,按該斷層轉(zhuǎn)換后的層位與縱向細(xì)化的層序界面間的接觸關(guān)系處理,斷層上盤層序界面與斷層關(guān)系按超覆處理,下盤層序界面與斷層關(guān)系按剝蝕處理,建立準(zhǔn)確的DX14井區(qū)AIW波阻抗反演的低頻阻抗模型。

2)通過井—震結(jié)合精細(xì)標(biāo)定鉆井有利儲層和縱向細(xì)化的層序界面,提取擬波阻抗曲線的低頻阻抗分量,在建立的斷層參與的高密度層控構(gòu)造框架控制下內(nèi)插和外推低頻阻抗建立低頻阻抗模型后,根據(jù)AIW 波阻抗反演基本原理,依據(jù)小波邊緣分析技術(shù)從經(jīng)過拓頻處理的提頻三維地震資料中提取反映巖性局部變化的地震屬性特征參數(shù),在斷層參與的高密度層控構(gòu)造框架控制下,用小波邊緣檢測方法得到的邊緣信息結(jié)合測井沿層位外推的波阻抗來估算反演的初始波阻抗模型,同擬聲波阻抗的低頻阻抗模型一起建立AIW波阻抗反演的最終初始模型,并且在迭代過程中應(yīng)用此參數(shù)對其進(jìn)行擾動限制,會使反演結(jié)果收獲意想不到的效果。

由于斷層參與的高密度層控構(gòu)造框架減小了層間測井信息內(nèi)插引起的層間構(gòu)造誤差、實現(xiàn)了斷層兩側(cè)內(nèi)插測井?dāng)?shù)據(jù)的突變,通過井—震精細(xì)標(biāo)定,依據(jù)小波邊緣分析技術(shù)從雙向拓頻地震數(shù)據(jù)中提取地震屬性特征參數(shù)不但補(bǔ)充了井間內(nèi)插造成的高頻成分缺失,也加入了反映巖性變化信息的地震屬性參數(shù),實現(xiàn)井—震一體化、構(gòu)造—巖性統(tǒng)一的波阻抗反演初始模型。

圖4和圖5分別為DX14井區(qū)3D地震L1036線測井聲波阻抗內(nèi)插外推得到的低頻模型以及小波邊緣分析得到的地震屬性特征參數(shù)聯(lián)合測井聲波阻抗、構(gòu)造框架建立的井—震統(tǒng)一的AIW波阻抗反演模型。對比圖4和圖5可見,井—震統(tǒng)一的AIW波阻抗反演模型不但包含測井聲波阻抗內(nèi)插外推的低頻模型井間橫向內(nèi)插的低頻成分,同時也引入了反映井間巖性變化的地震屬性數(shù)據(jù)的小波邊緣分析結(jié)果,在圖4中,由于井間聲波阻抗橫向連續(xù)內(nèi)插,無法體現(xiàn)井間巖性的變化和地層間的疊置關(guān)系,而在井—震統(tǒng)一的AIW波阻抗反演模型圖5中,地層縱向上巖性變化可以通過圖4中的低頻成分來區(qū)分,而橫向上巖性的變化、地層間接觸關(guān)系以及斷層、斷點(diǎn)則通過地震屬性特征參數(shù)來表征,同測井聲波阻抗內(nèi)插外推的低頻模型相比,井—震統(tǒng)一的AIW波阻抗反演模型更能反映地下的實際地質(zhì)情況。

圖4 L1036線測井阻抗低頻模型Fig.4 Low frequency logging impedance model of L1036 line

圖5 井震統(tǒng)一的AIW波阻抗反演模型Fig.5 AIW wave impedance inversion model for well seismology unification

4 應(yīng)用效果

在滴南凸起DX14井區(qū),梧桐溝組一段是梧桐溝組的主力含氣層系,上部砂巖發(fā)育,儲層較厚,下部泥巖增多。梧桐溝組一段沉積時期,構(gòu)造活動相對增強(qiáng),地貌高差構(gòu)造變化較大,受東部克拉美麗山物源控制,在DX14井區(qū)發(fā)育一系列扇三角洲前緣砂體,砂體頂部和底部主要發(fā)育泥巖、粉砂質(zhì)泥巖沉積,儲層結(jié)構(gòu)具有典型的“泥包砂”特點(diǎn),單層砂巖厚度一般2.0～7.0 m,前緣砂體累計厚度一般大于為5.0～18.0 m。應(yīng)用AIW波阻抗反演方法,充分考慮滴南凸起DX14井區(qū)復(fù)雜的構(gòu)造特點(diǎn)、多變的沉積環(huán)境形成的“小而薄”扇三角洲前緣砂體以及三維地震資料波組反射特征、地下全方位地震波場的分布特征與變化因素,依據(jù)小波邊緣分析技術(shù)從經(jīng)過拓頻處理的三維地震資料中提取反映巖性局部變化的地震屬性特征參數(shù),在斷層參與的高密度層控構(gòu)造框架控制的基礎(chǔ)上同測井?dāng)?shù)據(jù)一起參與初始模型建立,并反復(fù)迭代修改反演道的波阻抗模型,最終得到滴南凸起DX14井區(qū)高精度的地震反演波阻抗資料。

圖6為DX14井區(qū)在同一個構(gòu)造地質(zhì)框架模型基礎(chǔ)上,分別采用AIW波阻抗反演和常規(guī)基于模型波阻抗反演方法得到的過DX14井波阻抗反演剖面,從圖6可見,AIW波阻抗反演剖面的分辨率高于常規(guī)基于模型波阻抗反演結(jié)果和常規(guī)疊后地震剖面,AIW波阻抗反演結(jié)果反映的砂層間相互接觸關(guān)系不但清晰,而且明顯反映出了地層巖相、巖性的橫向變化,地層巖性信息也更加豐富。通過過DX14井的波阻抗反演剖面也可以看出梧桐溝組一段的砂體橫向分布范圍較小、不連續(xù),呈現(xiàn)出扇三角洲前緣砂體“小而薄”的特點(diǎn),符合在DX14井區(qū)發(fā)育扇三角洲前緣砂體的地質(zhì)認(rèn)識。縱向上準(zhǔn)確預(yù)測了DX14井梧桐溝組一段發(fā)育的5套砂體,其中3 532.00～3 536.00 m井段厚度為4.0 m的含氣砂體同3 542.00～3 550.00 m井段厚度為8.0 m的含油砂體呈明顯的疊置關(guān)系,而3 522.00～3 524.00 m井段厚度為2.0 m的含氣砂體、3 562.00～3 572.00 m井段厚度為10.0 m的含氣砂體則呈現(xiàn)明顯的透鏡狀砂體,向上傾和下傾方向尖滅,反演結(jié)果并與實鉆結(jié)果高度吻合,而常規(guī)基于模型波阻抗反演結(jié)果不但儲層分辨率不高,而且無法反映砂體間的疊置關(guān)系。

圖6 AIW波阻抗反演(a)與常規(guī)基于模型波阻抗反演(b)剖面對比Fig.6 Comparison between AIW wave impedance inversion section and conventional model-based wave impedance inversion section

利用AIW波阻抗反演和常規(guī)基于模型波阻抗反演結(jié)果,分別對梧桐溝組一段下砂組砂體頂、底界面進(jìn)行精細(xì)井震對比、追蹤解釋,構(gòu)造成圖,以梧桐溝組一段下砂組頂、底地震反射界面為層段時窗界限,按照砂巖和泥巖間波阻抗界限6 700 [(kg·m-3)·(m·s-1)]預(yù)測砂體分布結(jié)果如圖7,從中可見常規(guī)基于模型波阻抗反演資料分辨率較低,砂體的地震響應(yīng)特征不明顯,砂體邊界及巖性體間的關(guān)系不清(圖7b)。而利用AIW波阻抗反演結(jié)果可以明顯看出:受東部克拉美麗山物源控制,在DX14井區(qū)發(fā)育一系列扇三角洲前緣砂體,砂體與斷裂相配合在DX14井區(qū)形成斷層—巖性和巖性油氣藏(圖7a)。

圖7 DX14井區(qū)梧桐溝組一段波阻抗反演(a)與常規(guī)基于模型波阻抗反演(b)波阻抗平面Fig.7 P3wt2 1 AIW wave impedance inversion(a) and conventional model-based wave impedance inversion(b) wave impedance plan of the first member of Wutonggou Formation in DX14 well area

5 結(jié)論

1)應(yīng)用小波邊緣分析技術(shù)建立的初始波阻抗模型更加接近真實地下地質(zhì)情況,提高了反演的真實性,利用所得的特征點(diǎn)信息結(jié)合測井波阻抗對反演的約束起到了很好的效果,為儲層的研究提供了指導(dǎo)作用。

2)AIW波阻抗反演技術(shù)將測井波阻抗、地震特征邊緣信息作為約束條件建立的初始模型,彌補(bǔ)了常規(guī)反演過程中缺少井間的局部構(gòu)造、巖性變化的地震信息,增加了波阻抗反演過程中高頻成分,同常規(guī)反演相比,AIW波阻抗反演結(jié)果不但能夠準(zhǔn)確反映井間巖性變化,而且更接近實際地質(zhì)情況。

3)在斷層參與的高密度層控構(gòu)造框架控制下,測井波阻抗和地震屬性特征參數(shù)聯(lián)合建模實現(xiàn)了井—震一體化、構(gòu)造—巖性統(tǒng)一,但應(yīng)用的地震屬性特征參數(shù)比較單一,多屬性聯(lián)合建模仍是波阻抗反演的研究方向。

4)利用基于小波邊緣分析的井—震聯(lián)合建模AIW波阻抗高反演方法,對準(zhǔn)噶爾盆地陸梁隆起DX14井區(qū)梧桐溝組一段含氣砂巖儲層進(jìn)行精細(xì)反演和預(yù)測,刻畫了扇三角洲前緣砂體,預(yù)測砂體厚度的相對誤差均小于2.84%,預(yù)測精度高,完全適用于準(zhǔn)噶爾等復(fù)雜斷陷盆地以及復(fù)雜目標(biāo)區(qū)構(gòu)造—巖性等各種隱蔽性油氣藏的精細(xì)儲層預(yù)測和精細(xì)油藏描述。