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

有效

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

1.一種基于高維小波變換的地震資料不連續性檢測方法，其包括以下步驟：

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))]

其中，imag(F(x，y，t))表示對F(x，y，t)復地震道取虛部；real(F(x，y，t))表示對F(x，y，t)取實部；

進而獲得三維相位余弦數據體PhC(x)：

PhC(x，y，t)＝cos[Phase(x，y，t)]

3)對三維相位余弦數據體PhC(x)進行分塊，記第n塊三維分塊相位余弦數據體為PhCP(n，x)；

4)分別對每個三維相位余弦數據體的子塊PhCP(n，x)進行三維連續小波變換；

其中，x(x＝(x，y，t))表示平移因子；ψ(x)為三維母小波函數；a表示三維連續小波變換的尺度因子；θ表示地質目標選擇傾角；表示地質目標選擇方位角；表示與傾角θ及方位角有關的旋轉算子；PhCP_k(n，k)是PhCP(n，x)的波數域形式；是ψ(x)的波數域形式，其中k為波數；b表示臨時變量；

5)利用每個三維相位余弦數據體子塊的三維連續小波變換系數來進行三維地震資料的不連續性檢測，首先利用三維連續小波變換系數來確定地震資料不連續性的方向：

即將三維連續小波模最大值對應的傾角θ_dis(n，a，x)及方位角定義為x點的不連續性的方向，那么x點的不連續性度量為：

6)將生成的各個子塊的不連續性數據體Discon(n，a，x)進行拼接，得到整個地震資料的不連續性數據體，完成檢測。

2.如權利要求1所述的一種基于高維小波變換的地震資料不連續性檢測方法，其特征在于：所述步驟3)中，對三維相位余弦數據體PhC(x)進行分塊時，子塊各個方向上的點數盡量選擇2的整數次冪，為64*64*64或者128*128*128大的子塊。