多子波地震道分解技術在平湖油氣田的應用

高偉義，張 昊，安 平

(1.上海石油天然氣有限公司，上海 200041;2.杰奧世博(北京）技術有限公司，北京 100016）

多子波地震道分解技術在平湖油氣田的應用

高偉義1，張 昊1，安 平2

(1.上海石油天然氣有限公司，上海 200041;2.杰奧世博(北京）技術有限公司，北京 100016）

多子波地震道分解是一種新的地震道分解技術，它是以分解算法將地震道分解成不同形狀，不同頻率的子波。利用多子波地震道分解等技術，對東海平湖油氣田花港組油藏地震數據體進行多子波地震道分解與重構、頻譜衰減和數據體波形分解等應用研究，找出了與鉆井資料吻合，并含有地質意義的屬性，勾畫出該地區主要油藏的含油氣范圍和有利目標，取得了較好的應用效果。

平湖油氣田;多子波地震道分解與重構;頻譜衰減;數據體波形分解

平湖油氣田花港組H2、H6和H7油藏是多層系的巖性和構造復合油藏，砂體發育，儲層物性好。但由于該油藏為塊狀砂巖，各砂巖體相互疊置且含油氣情況不明，無法按構造油水邊界和常規地震解釋方法確定花港組油藏含油氣分布范圍。

針對這一情況，采用多子波地震道分解與重構技術，通過將地震道數據進行分解，尋找其中最能反映油藏分布的頻段進行重構。在此基礎上進行含油氣研究，預測了各砂體含油氣情況，圈定了含油氣范圍。

1 多子波分解與重構的原理

常規地震道解釋的基本模型是單一子波與地層反射系數系列的褶積，地震信號處理和解釋所用的褶積和反褶積以及地震道反演一般都是基于單一地震子波的假設。但是不同物理特性的地層，如儲層和非儲層，含油氣儲層和不含油氣儲層的地震響應是不同的，地震子波在穿過不同地層時所受到的改造也不同，其形狀會發生相應的變化。基于單一地震子波的褶積模型的一些常規儲層及油氣預測的方法就存在一定的局限性，對此引進了多子波地震道模型的概念。

設R(t）代表時間域的反射系數，Ri(t）代表只含有一個非零的反射系數序列，那么R(t）=(t）;而wj(t）代表一組不同形狀(不同頻率）的子波，則地震道

基于這一概念的多子波地震道分解方法就是將地震道分解成不同形狀、不同頻率雷克子波［1］的響應，如果將這些子波響應重新疊加，就可得到分解前的地震道。該技術突破了常規地震信號處理和解釋中單一子波的假設，地震道分解結果更符合客觀實際。

2 多子波地震道分解技術應用分析

在多子波地震道分解技術的應用中，首先在高品質的疊后地震資料上進行多子波地震道分解，在分解后的數據體上進行層位和斷層解釋，再結合鉆井資料，對不同的解釋目標，找出含油和不含油特征進行多子波地震道重構，并選擇其對應的不同地震窗口，在非重構和重構的地震數據體上進行頻譜、頻譜衰減、數據體波形分解等應用研究，找出與鉆井資料吻合，并且具有地質含義的屬性，進行含油氣預測研究，圈定含油氣范圍。下面就多子波地震道分解與重構技術的三種應用方式結合實例進行分析。

2.1 多子波地震道分解與重構［2］

2.1.1 地震層位解釋

地震資料的干擾訊號通常給解釋帶來困難，常規的頻率濾波只是去掉或減少干擾波所在的頻率，并不能去掉干擾子波。在應用多子波分解法將地震分解后，可對子波進行篩選，去掉干擾子波，提高信噪比。

應用多子波地震道分解與重構進行地震解釋時，首先將本區地震道分解成不同頻率的雷克子波，然后根據解釋層位特征，不斷改變分解后數據體的重構子波范圍，選擇能最佳反映該解釋層位特征，重構子波地震剖面進行層位對比解釋。圖1是本區一條穿過斷層的原始剖面和5～49 Hz的重構子波地震剖面，可見重構后的地震剖面信噪比高，F斷層更加清晰，更易進行層位解釋。

圖1 原始地震剖面(左）和5～49 Hz重構地震剖面(右）Fig.1 Original seismic section(left）and reconstruction with 5 ～49 Hz(right）

2.1.2 儲層和油氣預測

應用多子波地震道分解與重構技術進行儲層和含油氣預測時，假定同一地區和同一層位對同一類型的儲層或含油氣層具有相似的地震響應，根據已知鉆井資料，選取與該儲層或含油儲層變化相關的子波，去掉無直接關系的子波，合成最大限度反映該儲層或含油儲層變化的新的地震道。

圖2是本區花港組H7層位過A4和A1井的原始剖面和20～50 Hz的重構子波地震剖面，A4井鉆遇了H7油層，A1井鉆遇了H7水層，在原始剖面上二口鉆井都存在與儲層有關的紅色負相位同相軸，表明二口鉆井都存在相同儲層。在對子波進行篩選重構后，A1水井的子波幾乎沒有了，而A4油井的子波仍然很強，所以可清楚區別油層和水層，從而進行儲層和含油氣預測。

圖2 過A4和A1井的原始剖面(左）和20～50 Hz重構子波剖面(右）Fig.2 Original section(left）amd reconstruction with 20 ～50 Hz(right）which cross A4 well and A1 well

2.2 頻譜及頻譜衰減

地震波穿過不同物理特性的地層(油層或非油層）時，其頻率的衰減也不同。由于采用的是雷克子波進行分解，雷克子波的頻譜是一解析函數，通過子波就可計算任意大小的頻譜。可在預測層位的上下各取一個窗口，分別計算頻譜，并對上下頻譜進行能量歸一化處理，消除窗口內巖性的影響差，然后計算頻譜差，估算地震訊號穿過該層位時的頻譜衰減，從而進行含油氣預測。

本次在頻譜分析時采用了二種窗口:一是巖性窗口，該窗口主要包括目標層對應的數據段，主要反映該層的巖性特征。二是衰減窗口，該窗口主要對應目標層之下，反映地震波穿過該層后的衰減頻譜。圖3是本區花港組H2層連井線的巖性窗口頻譜、衰減窗口頻譜和H2層上下窗口的頻譜差(頻譜衰減）。

巖性窗口頻譜特征表明:①和③含油區巖性頻譜頻率成份偏高，⑤含水區巖性頻譜頻率偏低。衰減窗口頻譜表明這三個區的衰減頻譜具有明顯不同的特征，①區以10 Hz較強能量為特征，③區在10 Hz附近能量要低于①區，⑤區以20 Hz較強寬緩能量為特征。頻譜衰減過井線表明這三個區有著完全不同的頻譜特征，①區頻譜衰減嚴重，表示地震訊號穿過H2層位后高頻部分減少，低頻部分增加，較強的負值出現在10 Hz附近，③區也存在衰減，但沒有①區明顯，⑤區基本上沒有衰減，寬緩的負值分布在在20～30 Hz。

該研究成果與鉆井結果吻合度較高，①區對應為高含油飽和度的厚油層，儲層物性較好，③區對應為低含油飽和度的油層，儲層物性較差，⑤區為含水儲層。

圖3 巖性窗口頻譜(左）、衰減窗口頻譜(中）和頻譜衰減(右）Fig.3 Frequency spectrums of lithology window(left） ，frequency spectrums of attenuation window(mid）and attenuation of frequency spectrum(right）

2.3 地震體波形分解與重構［3］

地震體波形分解是一種提取線性相似地震波形狀的分解與合成方法，這一方法把地震數據分解成若干分量，第一分量的波形代表輸入地震數據中具有最大共性、最大能量的波形。第二分量代表去掉第一分量后的地震數據中具有最大共性、最大能量的波形，以此類推。因此，地震體波形分解為我們提供一個有效識別不同等級中不同地震響應的地質巖性和流體的有效方法。地質分析認為第一分量可反映該地區最大地質巖性或巖相分類，而相似的振幅可反映同一類型或相似的地質巖相或巖相。第二分量則反映次一級的地質巖性或巖相分類，以此類推。

應用該方法，針對花港組不同層位在不同的重構數據體上進行分解，找出最佳反映該油藏或巖性特征分布的窗口及分量，進行儲層和含油性預測。圖4為全子波重構地震數據體形分解得到的第一分量花港組H7層位相對振幅平面圖和20～55 Hz子波重構數據體形分解的第一分量H7層位相對振幅平面圖以及25～60 Hz子波重構數據體形分解的第一分量H6層位相對振幅平面圖。

與實際鉆井比較，全子波重構地震數據體形分解得到的第一分量更多的是反映H7儲層的分布，而20～55 Hz和25～60 Hz子波重構數據體形分解的第一分量更多的是反映含油儲層的分布。我們在該圖上勾畫出H7和H6層可能的含油范圍及有利區塊，該含油范圍和振幅顯示與鉆井結果吻合度較高，A1和A2井所鉆遇的H7層砂體為含水儲層，在其上覆的H6層砂體為含油層。

3 結論

本次應用多子波分解技術對平湖油氣田花港組的地震數據體進行了多子波分解，在此基礎上對花港組主要油藏進行了重構、頻譜、頻譜衰減和數據體形分解與重構等多項應用研究，得出如下結論:

圖4 全子波重構H7(左）和子波重構H7(中）及H6(右）的第一分量振幅平面圖Fig.4 First component amplitude map of original H7(left）and reconstruction of H7(mid）and H6(right）

(1）多子波地震道分解技術針對基于單一地震子波褶積模型的常規儲層及油氣預測方法存在的局限性，引進了多子波地震道模型概念，地震道分解結果符合客觀實際。

(2）在預測層位的上下分別計算頻譜和頻譜差，估算地震訊號的頻譜衰減，可以進行含油氣預測。

(3）地震體波形分解與重構提供了有效識別不同等級中不同地震響應的地質巖性和流體的方法。

(4）多子波分解和數據體積分解等技術所圈定的含油氣范圍與鉆井結果吻合度較高，具有較好的應用價值，在平湖油氣田油藏挖潛中起到了重要作

用。

［1］Ricker，N.The Form and Laws of Propagation of Seismic Wavelets［J］.Geophysics，1953，18(1）:10－40.

［2］Ping An.Application of multi-wavelet seismic trace decomposition and reconstruction to seismic data interpretation and reservoir characterization［C］//SEG International Exposition ＆ Annual Meeting，New Orleans.2006:973－977.

［3］Ping An.Case studies on stratigraphic interpretation and sand mapping using volume-based seismic waveform decomposition［C］//SEG.SEG International Exposition＆ Annual Meeting，New Orleans:2006:496－499.

Application study of multi-wavelet seismic trace decomposition and reconstruction in Pinghu oil and gas field

GAO Weiyi1，ZHANG Hao1，AN Ping2
(1.Shanghai Petroleum Co.Ltd.，Shanghai 200041，China;2.GeoCyber Solutions Inc.，Beijing 100016，China）

Multi-wavelet seismic trace decomposition and reconstruction is a new kind of seismic trace decomposition technology，which separates the seismic traces into a series of wavelets with different shape and different frequency using the decomposition algorithm.By applying the multi-wavelet seismic trace decomposition and reconstruction，study on frequency spectrum attenuation and volume based waveform decomposition processing on the seismic data which represents the oil reservoirs in Huagang formation was carried on，those attributes which have the distinct geological meaning and good relativities with wells were found out，and hydrocarbon bearing area in different reservoir were clearly portrayed.

Pinghu oil and gas field;multi-wavelet seismic trace decomposition and reconstruction;frequency spectrum attenuation;volume based waveform decomposition

P631.4

A

10.3969/j.issn.1008－2336.2011.04.020

1008－2336(2011）04－0020－04

2011－06－10;改回日期:2011－07－01

高偉義，男，1962年生，高級工程師，碩士，現主要從事石油勘探與開發研究工作。E-mail:gaowy@shpc.com.cn。