[發明專利]一種基于高維小波變換的地震資料不連續性檢測方法

有效

• 申請號: 201110134352.2
• 申請日: 2011-05-23
• 公開（公告）號: CN102298155A
• 公開（公告）日: 2011-12-28
• 發明（設計）人: 朱偉林;趙偉;高靜懷;王曉凱;陳文超;姜秀娣;宋亮;翁斌;安琪;劉永江;李麗霞;王清振;糜芳
• 申請（專利權）人: 中國海洋石油總公司;中海石油研究中心;西安交通大學
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京紀凱知識產權代理有限公司 11245
• 代理人: 徐寧;關暢

技術領域

本發明涉及一種油氣勘探中地震資料的檢測方法，特別是關于一種基于高維小波變換的地震資料不連續性檢測方法。

背景技術

隨著地震勘探技術的深入和發展，油氣勘探已經由常規的構造性油氣藏勘探向隱蔽油氣藏勘探轉變。利用地震資料檢測地層傾角、地下斷層、裂縫、巖性變化以及地質體的邊緣等結構對隱蔽油氣藏勘探有著重要的意義。

目前業界通常采用三維地震勘探技術，其產生的數據量巨大，即便是疊后數據，其數據量一般也為幾十GB(千兆)。因此，在處理如此龐大的數據時，所采用技術的高效快速性尤為重要。油氣勘探中所采集的海量地震數據中往往夾雜著各種噪聲，而各種噪音對后期處理將產生嚴重的影響。故而所采用技術必須具有良好的抗噪性能。

常規的地震資料橫向不連續性檢測技術包括第一代相干體算法、基于相似性度量的第二代相干體算法、基于特征結構分解的第三代相干體算法以及局部結構熵算法。第一代相干體算法采用相鄰道之間的相關來度量橫向不連續性，對噪聲敏感；第二代相干體算法采用分析窗內各道的相似性來度量橫向不連續性，分辨率受限；第三代相干體算法對分析窗內各道構成的協方差矩陣進行特征分解，利用最大特征值與所有特征值之比來度量地震資料的橫向不連續性。上述各算法的時間分析窗通常要大于等于一個子波的持續期，且其運算量隨著分析窗的增大而急劇增大。

由于2DCWT(2-Dimensional?Continuous?Wavelet?Transform，兩維連續小波變換)可以通過二維FFT(Fast?Fourier?Transform，快速傅里葉變換)來快速實現，并且可以對地震信號進行多分辨分析，即可以在受噪聲干擾較小的尺度下分析地震資料，故而2DCWT可以解決常用技術的運算量大和抗噪性能差的問題。Bouchereau于1997年首次將2DCWT用于檢測地震資料的不連續結構，所使用的小波為二維Morlet小波，但二維Morlet小波在低波數時其空間局部化性質較差，因此不能大規模用于地震資料處理。另外，目前業內經常使用三維地震資料，二維算法并不能充分利用三維地震資料的空間相關性。

發明內容

針對上述問題，本發明的目的是提供一種基于高維連續小波變換的地震資料不連續性檢測方法，該方法可充分利用三維地震資料的空間相關性，大規模用于地震資料處理。

為實現上述目的，本發明采取以下技術方案：一種基于高維小波變換的地震資料不連續性檢測方法，其包括以下步驟：1)采集地震資料，預處理后得到三維地震數據體f(x，y，t)，其中，x，y表示位置坐標，t表示時間軸；2)對三維地震數據體f(x)每道進行Hilbert變換，得到復地震道F(x)：

F(x,y,t)=f(x,y,t)+i1π∫-∞+∞f(x,y,b)1t-bab]]>

其中，x＝(x，y，t)；b表示臨時變量；在此基礎上提取瞬時相位得到三維瞬時相位數據體Phase(x)：

Phase(x，y，t)＝tan^-1[imag(F(x，y，t))/real(F(x，y，t))]