一种三维地震观测系统多频率偏移噪声获取方法及装置

摘要

本发明实施例提供一种三维地震观测系统多频率偏移噪声获取方法及装置，所述方法包括：获取三维地震观测系统数据；根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声；对所述三维地震观测系统的单频率偏移噪声进行插值，获取所述三维地震观测系统的单频率偏移噪声之外的一种或多种频率下偏移噪声；对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声。本发明实施例仅仅进行一次单频率偏移噪声运算，在保证计算精度的同时不增加算法的整体计算量，从而实现三维地震观测系统多频率偏移噪声的快速获取。

主权项

一种三维地震观测系统多频率偏移噪声获取方法，其特征在于，所述三维地震观测系统多频率偏移噪声获取方法包括：获取三维地震观测系统数据；根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声，包括：对于均匀介质，利用如下公式计算单频率偏移噪声：$B_D\left(f\right)\approx \frac{1}{4{\pi }^2}\int \int \exp \left[i(k_x\Delta x+k_y\Delta y)\right]D^{\prime }{\mathrm{dk}}_x{\mathrm{dk}}_y,$其中，i为虚数单位，k_x和k_y分别为水平x方向和水平y方向波数，Δx和Δy分别为检波点到地下目标反射点x方向和y方向的水平距离分量，D′为检波点矩阵D的波数域变换，B_D(f)为单频率偏移噪声；对于非均匀介质，利用如下公式计算单频率偏移噪声：B_D(f)≈f²∫∫exp[2πif(p_xΔx+p_yΔy)]D″dp_xdp_y，其中，i为虚数单位，f为频率，p_x和p_y分别为水平x方向和水平y方向慢度，Δx和Δy分别为检波点到地下目标反射点x方向和y方向的水平距离分量，D″为检波点矩阵D的Radon域变换，B_D(f)为单频率偏移噪声；对所述三维地震观测系统的单频率偏移噪声进行插值，获取所述三维地震观测系统的单频率偏移噪声之外的一种或多种频率下偏移噪声；对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声，包括：利用如下公式计算所述三维地震观测系统的多频率偏移噪声：$B\approx \sqrt{\Delta f\underset{j=1}{\overset{N}{\Sigma }}\left[j^2{\left|w\left(j\Delta f\right)B(\frac{j\Delta x}{N},\frac{j\Delta y}{N},N\Delta f)\right|}^2\right]},$其中，B为多频率偏移噪声，N是频率采样点数，Δf为频率采样间隔，w(f)为雷克子波频率域表达式f_p为地震波主频。