張大海

(中國石油化工股份有限公司上海海洋油氣分公司研究院,上海 200120)

基于地震速度構(gòu)建的低頻響應(yīng)在波阻抗反演中的應(yīng)用

張大海

(中國石油化工股份有限公司上海海洋油氣分公司研究院,上海 200120)

以V工區(qū)北部為例,應(yīng)用井?dāng)?shù)據(jù)、速度數(shù)據(jù)和地震屬性共同恢復(fù)低頻響應(yīng)數(shù)據(jù),同時使用層速度、層位信息(層的深度和厚度)以及反射振幅的信息進(jìn)行反演,結(jié)果顯示可以改善預(yù)測縱波阻抗的結(jié)果,同時與通常情況下的孔隙度預(yù)測作了對比,結(jié)果表明孔隙度預(yù)測的可信度得到了一定程度的改善。

低頻響應(yīng);地震速度;孔隙度;波阻抗反演;V工區(qū)

1 背景及分析

聲波阻抗作為一種巖石的物性,可以通過聲波時差和密度的測井曲線計算得到,在一定的約束條件下可以由地震數(shù)據(jù)推導(dǎo)得到,同時為其它儲層參數(shù)(孔隙度、巖性、流體性質(zhì)的量等)的深入分析提供重要基礎(chǔ)。地震反演是一個將地震數(shù)據(jù)轉(zhuǎn)化為相對聲波阻抗的過程,結(jié)果是對遠(yuǎn)離井眼處的相對聲波阻抗的估計。假定相對聲波阻抗中缺失的低頻成分可以通過計算得到補(bǔ)償?shù)脑?就可以計算出絕對聲波阻抗進(jìn)而對其它遠(yuǎn)離井眼位置的地區(qū)其它物性參數(shù)進(jìn)行預(yù)測[1]。

與井上的聲波阻抗曲線對比,反演的相對波阻抗在頻譜的高頻和低頻段均有一定的限制,但也沒有限制相對波阻抗是一個有用的地震屬性,有助于更好的定性理解儲層物性,尤其對于碎屑巖儲層,經(jīng)常被用于其他儲層物性的定量估計。碳酸巖的巖相比較均一,巖石骨架特征和孔隙之間的關(guān)系清晰,是一個典型的可以定量預(yù)測的介質(zhì)。但是,碳酸巖的定量儲層預(yù)測受低頻信息缺失的影響很大,聲波阻抗和孔隙度之間的強(qiáng)烈相關(guān)受到低頻趨勢影響很大。

反演的波阻抗數(shù)據(jù)中缺失的低頻成分可以通過多井位置上的波阻抗曲線橫向插值模擬得到補(bǔ)充。通常情況下,井之間的插值受距離影響很大,容易導(dǎo)致假象和非地質(zhì)結(jié)果?；诰嚯x的井資料橫向插值可以通過使用地震速度作為背景資料加以改善,盡管速度的使用是一個很大的改善,但還是受到隨深度增加分辨率和精度大大降低的影響。速度數(shù)據(jù)僅在低頻范圍部分提供缺失的信息。

圖1顯示了一條地震曲線和它的頻譜。同時對波阻抗曲線作了60 Hz以上高切,分辨率明顯降低。

圖1 地震道(左)、及其振幅譜(中)、波阻抗曲線以及其高切濾波曲線(右)Fig.1 Seismic(left)and its amplitude spectrum(middle),impedance log and high-cut filtering curve(right)

圖2是V1井的縱波阻抗曲線,右側(cè)的曲線是縱波阻抗測井曲線和低頻響應(yīng),左側(cè)的曲線是對應(yīng)的相對縱波阻抗曲線。相對縱波阻抗的曲線是由測量得到的縱波阻抗曲線中濾掉低頻響應(yīng)曲線得到的。在縱波阻抗曲線上可以看到深度2 400～2 600 m時阻抗值的下降趨勢以及2 700 m的跳躍特征,這些特征在相對阻抗曲線上是看不到的。顯然對于絕對縱波阻抗這些特征的識別,低頻響應(yīng)起著很重要的作用。

圖2 井上相對波阻抗曲線(左)和波阻抗測井曲線(右)示意Fig.2 Sketch map of relative impedance log(left)and impedance log(right)

2 低頻響應(yīng)的重要性

一般井?dāng)?shù)據(jù)的測量分辨率在0.5英寸,因此測井?dāng)?shù)據(jù)具有很寬的空間頻率響應(yīng),在時間域上的響應(yīng)對應(yīng)的范圍從接近于0一直到幾千赫茲。通常的地震反射數(shù)據(jù)是4 ms采樣率其帶寬可達(dá)到125 Hz的數(shù)據(jù),但是有用的信號帶寬通常不超過6～60 Hz的范圍。實際的頻帶寬度受到多種因素的影響,如接收器的幾何排列,采樣率,脈沖源(強(qiáng)度和頻率),以及地下介質(zhì)的彈性和阻尼情況等。和井?dāng)?shù)據(jù)相比,地震數(shù)據(jù)缺失低頻和高頻成分[2]。

相對縱波阻抗盡管本身存在一定的局限性,但是它還是一個很有用的屬性,因為它的高分辨率可以用于詳細(xì)的地層解釋,并且可以提供儲層性質(zhì)的額外信息。只要能在合適的頻段通過井上交會分析找到相對縱波阻抗和其他儲層性質(zhì)之間的關(guān)系,相對縱波阻抗就可以作為巖性、孔隙度以及含烴指示的量。在比較純凈的石灰?guī)r中由于孔隙和相對縱波阻抗存在很強(qiáng)的相關(guān),所以可以被用于孔隙預(yù)測,但是當(dāng)存在含烴流體時它的敏感度受到一定的限制。當(dāng)儲層中含有天然氣或高水平的溶解氣的時候都會對阻抗值在水平上有一定的影響。所以為了從縱波阻抗推導(dǎo)出孔隙度,恢復(fù)低頻響應(yīng)和建立絕對縱波阻抗顯得尤為重要。

圖3是縱波阻抗與孔隙度曲線的交會分析,其中顏色代表含水飽和度,暖色調(diào)為高含水飽和度。左面的圖為全帶寬縱波阻抗與孔隙度交會圖,右面圖為將全帶寬縱波阻抗曲線通過帶寬濾波得到的帶寬縱波阻抗曲線與孔隙度的交會分析?？梢钥吹阶竺鎴D的相關(guān)性明顯較好,且低含水飽和度與低的絕對縱波阻抗相關(guān),明顯與含烴有關(guān)。同樣含烴的影響在相對波阻抗與孔隙度的交會圖中很難被識別。

圖3 相對波阻抗、絕對波阻抗與孔隙度交會圖Fig.3 Crossplot of relative impedance,absolute impedance and porosity

因此,為了使用反演的波阻抗進(jìn)行孔隙度的預(yù)測,包括全頻段的信息是很重要的,也就是說必須彌補(bǔ)地震數(shù)據(jù)體中缺失的低頻的成分。

3 應(yīng)用實例

本文以V工區(qū)北部為例,使用層速度、層位信息(層的深度和厚度)以及反射振幅的信息,結(jié)果顯示可以明顯改善預(yù)測低頻響應(yīng)的精度,同時與通常情況下的孔隙度預(yù)測作了對比。這樣地震速度(0～2 Hz帶寬)、井?dāng)?shù)據(jù)(2～8 Hz)和地震數(shù)據(jù)(高于8 Hz)共同提供了頻率豐富的資料基礎(chǔ)。

圖4為V工區(qū)北部Dg-I層局部構(gòu)造圖,該地區(qū)為西高東低,地勢落差較大,該構(gòu)造為斷背斜構(gòu)造,構(gòu)造高部位與東部凹陷區(qū)時間深度相差大約一秒多,凹陷處的縱波阻抗應(yīng)明顯高于構(gòu)造高部位的縱波阻抗,但由于該工區(qū)斷裂系統(tǒng)發(fā)育,且已知井較少,構(gòu)造高部位的V1井很難控制東部凹陷區(qū),因此在由疊后的波阻抗反演數(shù)據(jù)體得到的該層平面展布圖(圖5左)中可以看出東部凹陷區(qū)的縱波阻抗并沒有升高,這與構(gòu)造解釋相矛盾。當(dāng)我們增加了由地震速度構(gòu)建的低頻響應(yīng)到地震反演的波阻抗體中,得到的平面展布圖(圖5右)更加合理,便于解釋,進(jìn)而又從井上的縱波阻抗和孔隙度之間的關(guān)系得到了孔隙度的預(yù)測平面圖。

圖4 V工區(qū)北部Dg-I層局部構(gòu)造圖Fig.4 Local structure map in the north of area V

4 小結(jié)

由于地震頻帶寬度的限制,很難得到全頻帶地震反演數(shù)據(jù)體,這對儲層物性的預(yù)測帶來了很大的限制,尤其在井比較少的情況下。通過井上分析可知,低頻成分對物性之間關(guān)系的建立起著很重要的作用。

圖5 疊后波阻抗反演(左)和應(yīng)用地震速度建模(右)的反演結(jié)果Fig.5 Poststack impedance inversion result without(left)and with(right)seismic velocity modeling

本文以V工區(qū)北部為例,在已有低頻響應(yīng)模型基礎(chǔ)上,通過應(yīng)用地震速度構(gòu)建的低頻響應(yīng)模型進(jìn)行反演得到的縱波阻抗數(shù)據(jù)體,改善了波阻抗反演結(jié)果的精度,進(jìn)而增加了孔隙度預(yù)測的可信度。

[1]Jason Geoscience Workbench Manual[M].2001.

[2]Hampson D P,Schuelke J S,QuireinJ A.Use of multiattribute transforms to predict log properties from seismic data[J].Geophysics,2000,66(1):220-236.

The application of low-frequency response model based on seismic velocity to acoustic impedance inversion

Zhang Dahai
(Institute of SINOPEC Shanghai Offshore Oil&Gas Company,Shanghai200120)

Low-frequency response model based on well data,velocity data and seismic attributes was applied to acoustic impedance inversion.For the northern part of area V,inversion was done using interval velocity,horizon information and reflection amplitude information,and the results showed P-impedance prediction and porosity prediction were both improved.

low-frequency response;seismic velocity;porosity;acoustic impedance inversion;area V

book=3,ebook=143

P631.4+2

A

10.3969/j.issn.1008-2336.2010.03.015

1008-2336(2010)03-0015-04

2010-04-12;改回日期:2010-06-10

張大海,男,高級工程師,2001年畢業(yè)于浙江大學(xué)地球科學(xué)系,獲理學(xué)博士學(xué)位,現(xiàn)從事儲層及油氣預(yù)測工作。E-mail:zhangdh @shopc.com.cn。