一种基于高维小波变换的地震资料不连续性检测方法

摘要

本发明涉及一种基于高维小波变换的地震资料不连续性检测方法，其包括以下步骤：1)采集地震资料，预处理后得到三维地震数据体；2)对三维地震数据体每道进行Hilbert变换，以此得到三维相位余弦数据体；3)对三维相位余弦数据体进行分块；同时根据不连续的尺度和地质目标选择尺度、倾角和方位角的搜索范围；4)对每个三维相位余弦数据体子块进行三维连续小波变换；5)利用每个三维相位余弦数据体子块的三维连续小波变换系数模给出地震资料的不连续性度量，得到分块不连续性数据体；6)将生成的分块不连续数据体进行拼接，得到整个地震资料的不连续性数据体。本发明方法运算量少，能够多尺度的检测地震资料的不连续性，因此，可广泛用于地震资料的检测过程中。

主权项

1.一种基于高维小波变换的地震资料不连续性检测方法，其包括以下步骤：1)采集地震资料，预处理后得到三维地震数据体f(x，y，t)，其中，x，y表示位置坐标，t表示时间轴；2)对三维地震数据体f(x)每道进行Hilbert变换，得到复地震道F(x)：$F(x,y,t)=f(x,y,t)+i\frac{1}{\pi }{\int }_{-\infty }^{+\infty }f(x,y,b)\frac{1}{t-b}\mathrm{ab}$其中，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)进行拼接，得到整个地震资料的不连续性数据体，完成检测。