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

Abstract

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

Description

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

技术领域

本发明涉及三维地震观测系统分析评价技术，尤其涉及一种三维地震观测系统多频率偏移噪声获取方法及装置。

背景技术

三维地震勘探是石油与天然气勘探的主要工具。三维地震勘探的主要过程包括：(1)数据资料采集。在陆上地震勘探数据采集现场一般进行下述三个工作：地震观测系统设计、根据设计方案在野外布设震源和检波器，地震波的激发和接收。首先在室内进行地震观测系统设计，以确定震源点和检波点的最佳摆放位置。然后根据设计方案，在野外布设震源和检波器，陆上地震勘探的震源点一般采用炸药震源，并沿地震测线等间距布置多个检波器来接收地震波信号，现代地震勘探中检波器的数量多在1000或10000个。震源点在爆炸后产生出地震波，地震波遇岩层界面反射回来被检波器接收并传到仪器车，仪器车将检波器传来的信号记录下来，这就获得了用以研究地下油气埋藏情况的地震记录。(2)地震数据处理。地震数据处理是把第一步骤采集到的地震数据资料输入专用电子计算机，按不同要求用一系列功能不同的程序进行处理运算，把数据进行归类编排，突出有效的，除去无效和干扰的，最后把经过各种处理的数据进行叠加和偏移，最终得到二维或三维地震数据体文件。(3)资料解释。资料解释是把经过处理的地震信息变成地质成果的过程，包括运用波动理论和地质知识，综合地质、钻井、测井等各项资料，作出构造解释、地层解释、岩性和烃类检测解释及综合解释，绘出有关成果图件，对勘探区域作出含油气评价，提出钻探井位置等。

在整个地震勘探的处理流程中，地震叠前偏移方法是最重要的地震资料处理流程之一，其效果直接决定了地震资料处理的成败甚至整个地震勘探的成败。前面的介绍曾提到，三维地震观测系统设计是整个地震勘探的第一个环节，因此，必须从三维地震观测系统设计阶段考虑地震叠前偏移的效果(即偏移噪声)，以选择合适的地震观测系统设计方案来完成地震勘探任务。目前直接采用地震叠前偏移理论计算三维地震观测系统偏移噪声的方法主要有两种：(1)基于Beylkin理论，结合检波点与震源点的分析方法；(2)Berkhout共聚焦分析方法。目前这两种方法均采用地震波射线追踪或地震波Kirchhoff算子实现，因此仅适用于均匀介质或层状介质情况，无法适用于复杂的地下介质条件，且计算精度较低。

发明内容

本发明实施例提供一种三维地震观测系统多频率偏移噪声获取方法及装置，通过仅仅进行一次单频率偏移噪声运算，利用快速偏移噪声插值算法将其推广至多频率条件下，在保证计算精度的同时不增加算法的整体计算量，从而实现三维地震观测系统多频率偏移噪声的快速获取。

一方面，本发明实施例提供了一种三维地震观测系统多频率偏移噪声获取方法，所述三维地震观测系统多频率偏移噪声获取方法包括：获取三维地震观测系统数据；根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声；对所述三维地震观测系统的单频率偏移噪声进行插值，获取所述三维地震观测系统的单频率偏移噪声之外的一种或多种频率下偏移噪声；对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声。

另一方面，本发明实施例提供了一种三维地震观测系统多频率偏移噪声获取装置，所述三维地震观测系统多频率偏移噪声获取装置包括：数据获取单元，用于获取三维地震观测系统数据；单频率偏移噪声单元，用于根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声；插值单元，用于对所述三维地震观测系统的单频率偏移噪声进行插值，获取所述三维地震观测系统的单频率偏移噪声之外的一种或多种频率下偏移噪声；多频率偏移噪声单元，用于对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声。

上述技术方案具有如下有益效果：因为采用所述三维地震观测系统多频率偏移噪声获取方法包括：获取三维地震观测系统数据；根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声；对所述三维地震观测系统的单频率偏移噪声进行插值，获取所述三维地震观测系统的单频率偏移噪声之外的一种或多种频率下偏移噪声；对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声的技术手段，所以达到了如下的技术效果：克服了现有三维地震观测系统多频率偏移噪声获取方法运算速度较慢的缺陷，提出了一种可以适用于简单或复杂模型的三维地震观测系统多频率偏移噪声的快速获取方法。可将仅仅进行一次单频率偏移噪声运算，利用快速偏移噪声插值算法将其推广至多频率条件下，在保证计算精度的同时不增加算法的整体计算量，从而实现三维地震观测系统多频率偏移噪声的快速获取。本发明实施例大幅度提升多频率偏移噪声计算的速度，有力地促进了偏移噪声分析技术在实际地震勘探中的应用。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例一种三维地震观测系统多频率偏移噪声获取方法流程图；

图2为本发明实施例一种三维地震观测系统多频率偏移噪声获取装置结构示意图；

图3为本发明应用实例多频率偏移噪声计算过程的示意图；

图4为本发明应用实例某三维地震观测系统多频率偏移噪声分析结果示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

如图1所示，为本发明实施例一种三维地震观测系统多频率偏移噪声获取方法流程图，所述三维地震观测系统多频率偏移噪声获取方法包括：

101、获取三维地震观测系统数据；

优选的，所述三维地震观测系统数据包括：震源点的三维空间坐标、检波点三维空间坐标、震源点与检波点之间的对应关系；其中，所述震源点的三维空间坐标包括：高程；所述检波点三维空间坐标包括：高程。

102、根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声；

优选的，所述根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声，包括：

对于均匀介质，利用如下公式计算单频率偏移噪声：

$B_D\left(f\right)\≈\frac{1}{4{\mathrm{\π}}^2}\∫\∫\exp \left[i(k_x\mathrm{\Δx}+k_y\mathrm{\Δy})\right]D^{\′}{\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)为单频率偏移噪声。

103、对所述三维地震观测系统的单频率偏移噪声进行插值，获取所述三维地震观测系统的单频率偏移噪声之外的一种或多种频率下偏移噪声；

104、对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声。

优选的，所述对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声，包括：

利用如下公式计算所述三维地震观测系统的多频率偏移噪声：

$B\≈\sqrt{\mathrm{\Δf}\underset{j=1}{\overset{N}{\mathrm{\Σ}}}\left[j^2{\left|w\left(\mathrm{j\Δf}\right)B(\frac{\mathrm{j\Δx}}{N},\frac{\mathrm{j\Δy}}{N},\mathrm{N\Δf})\right|}^2\right]},$

其中，B为多频率偏移噪声，N是频率采样点数，Δf为频率采样间隔，w(f)为雷克子波频率域表达式 $w\left(\mathrm{j\Δf}\right)=\frac{2}{\sqrt{\mathrm{\π}}}\frac{{\left(\mathrm{j\Δf}\right)}^2}{f_p^3}\exp \left[\frac{-{\left(\mathrm{j\Δf}\right)}^2}{f_p^2}\right],$ f_p为地震波主频。

对应于上述方法实施例，如图2所示，为本发明实施例一种三维地震观测系统多频率偏移噪声获取装置结构示意图，所述三维地震观测系统多频率偏移噪声获取装置包括：

数据获取单元21，用于获取三维地震观测系统数据；

单频率偏移噪声单元22，用于根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声；

插值单元23，用于对所述三维地震观测系统的单频率偏移噪声进行插值，获取所述三维地震观测系统的单频率偏移噪声之外的一种或多种频率下偏移噪声；

多频率偏移噪声单元24，用于对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声。

优选的，所述三维地震观测系统数据包括：震源点的三维空间坐标、检波点三维空间坐标、震源点与检波点之间的对应关系；其中，所述震源点的三维空间坐标包括：高程；所述检波点三维空间坐标包括：高程。

优选的，所述单频率偏移噪声单元22对于均匀介质，利用如下公式计算单频率偏移噪声：

$B_D\left(f\right)\≈\frac{1}{4{\mathrm{\π}}^2}\∫\∫\exp \left[i(k_x\mathrm{\Δx}+k_y\mathrm{\Δy})\right]D^{\′}{\mathrm{dk}}_x{\mathrm{dk}}_y,$

其中，i为虚数单位，k_x和k_y分别为水平x方向和水平y方向波数，Δx和Δy分别为检波点到地下目标反射点x方向和y方向的水平距离分量，D′为检波点矩阵D的波数域变换，B_D(f)为单频率偏移噪声。

优选的，所述单频率偏移噪声单元22对于非均匀介质，利用如下公式计算单频率偏移噪声：

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)为单频率偏移噪声。

优选的，所述多频率偏移噪声单元24利用如下公式计算所述三维地震观测系统的多频率偏移噪声：

$B\≈\sqrt{\mathrm{\Δf}\underset{j=1}{\overset{N}{\mathrm{\Σ}}}\left[j^2{\left|w\left(\mathrm{j\Δf}\right)B(\frac{\mathrm{j\Δx}}{N},\frac{\mathrm{j\Δy}}{N},\mathrm{N\Δf})\right|}^2\right]},$

其中，B为多频率偏移噪声，N是频率采样点数，Δf为频率采样间隔，w(f)为雷克子波频率域表达式 $w\left(\mathrm{j\Δf}\right)=\frac{2}{\sqrt{\mathrm{\π}}}\frac{{\left(\mathrm{j\Δf}\right)}^2}{f_p^3}\exp \left[\frac{-{\left(\mathrm{j\Δf}\right)}^2}{f_p^2}\right],$ f_p为地震波主频。

本发明实施例上述技术方案具有如下有益效果：因为采用所述三维地震观测系统多频率偏移噪声获取方法包括：获取三维地震观测系统数据；根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声；对所述三维地震观测系统的单频率偏移噪声进行插值，获取所述三维地震观测系统的单频率偏移噪声之外的一种或多种频率下偏移噪声；对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声的技术手段，所以达到了如下的技术效果：克服了现有三维地震观测系统多频率偏移噪声获取方法运算速度较慢的缺陷，提出了一种可以适用于简单或复杂模型的三维地震观测系统多频率偏移噪声的快速获取方法。可将仅仅进行一次单频率偏移噪声运算，利用快速偏移噪声插值算法将其推广至多频率条件下，在保证计算精度的同时不增加算法的整体计算量，从而实现三维地震观测系统多频率偏移噪声的快速获取。本发明实施例大幅度提升多频率偏移噪声计算的速度，有力地促进了偏移噪声分析技术在实际地震勘探中的应用。

以下结合应用实例对本发明实施例上述技术方案进行详细说明：

通常偏移噪声分析均在单频率下进行,然而实际的地震波频率总是会有一定的带宽，即总是多频率的。为了解决上述矛盾，必须将单频率偏移噪声分析扩展至多频率偏移噪声分析。如图3所示，为本发明应用实例多频率偏移噪声计算过程的示意图。多频率偏移噪声结果可用下式表示：

$B=\sqrt{\mathrm{\Δf}\underset{j=1}{\overset{N}{\mathrm{\Σ}}}{\left|w\left(\mathrm{j\Δf}\right)B\left(\mathrm{j\Δf}\right)\right|}^2}---\left(1\right)$

方程(1)表示所有频率波场的能量和,其中,B为多频率偏移噪声，N是频率采样点数，Δf为频率采样间隔，B(jΔf)为单频率偏移噪声，w(f)为雷克子波频率域表达式 $w\left(\mathrm{j\Δf}\right)=\frac{2}{\sqrt{\mathrm{\π}}}\frac{{\left(\mathrm{j\Δf}\right)}^2}{f_p^3}\exp \left[\frac{-{\left(\mathrm{j\Δf}\right)}^2}{f_p^2}\right],$ f_p为地震波主频。

对于均匀介质，偏移噪声公式可表示为：

$B_D\left(f\right)\≈\frac{1}{4{\mathrm{\π}}^2}\∫\∫\exp \left[i(k_x\mathrm{\Δx}+k_y\mathrm{\Δy})\right]D^{\′}{\mathrm{dk}}_x{\mathrm{dk}}_y---\left(2\right)$

其中，i为虚数单位，k_x和k_y分别为水平x方向和水平y方向波数，Δx和Δy分别为检波点到地下目标反射点x方向和y方向的水平距离分量，D′为检波点矩阵D的波数域变换。

对于非均匀介质，偏移噪声公式可表示为：

B_D(f)≈f²∫∫exp[2πif(p_xΔx+p_yΔy)]D″dp_xdp_y   (3)

其中，i为虚数单位，f为频率，p_x和p_y分别为水平x方向和水平y方向慢度(和速度相对，是速度的倒数)，Δx和Δy分别为检波点到地下目标反射点x方向和y方向的水平距离分量，D″为检波点矩阵D的Radon域变换。

综合方程(2)和(3),可得：

$B(\mathrm{\Δx},\mathrm{\Δy},\mathrm{j\Δf})\≈\frac{1}{j^2}B(\mathrm{j\Δx},\mathrm{j\Δy},\mathrm{j\Δf})---\left(4\right)$

代入方程(1)，即可得到多频率偏移噪声快速获取算法：

$B\≈\sqrt{\mathrm{\Δf}\underset{j=1}{\overset{N}{\mathrm{\Σ}}}\left[j^2{\left|w\left(\mathrm{j\Δf}\right)B(\frac{\mathrm{j\Δx}}{N},\frac{\mathrm{j\Δy}}{N},\mathrm{N\Δf})\right|}^2\right]}---\left(5\right)$

方程(5)表明,进行一次单频率偏移噪声分析与插值运算,即可近似得到三维地震观测系统多频率偏移噪声结果，如图4所示，为本发明应用实例某三维地震观测系统多频率偏移噪声分析结果示意图。

参照图1，本发明应用实例包括以下步骤：

1)导入复杂地表下某地震观测系统设计方案

观测系统设计方案包括检波点的三维空间坐标(含高程)、震源点三维空间坐标(含高程)，和震源点与检波点之间的对应关系。一个地震观测系统设计方案包含大量的震源点和检波点的三维坐标信息。这里的震源点指地震勘探普遍采用的人工炸药震源。检波点指地震勘探中在地表布设的检波器，一般至少有上千个，用于接收震源发出经地下反射后回到地表的地震波。

2)计算三维地震观测系统单频率偏移噪声

利用方程(2)或(3)计算三维地震观测系统单频率偏移噪声。

3)插值得到其它频率偏移噪声结果

利用方程(4)插值得到其它频率下偏移噪声结果。

4)对所有频率成分求和得到多频率偏移噪声结果

利用方程(5)得到多频率偏移噪声，实现仅利用一次单频率偏移噪声分析与插值运算,即可近似得到多频率下偏移噪声分析结果。

参照图4所示本发明应用实例某三维地震观测系统多频率偏移噪声分析结果，与传统多频率偏移噪声分析不同，该技术仅需要进行一次单频率偏移噪声运算，进而利用快速插值算法得到多频率偏移噪声结果，在保证计算精度的同时不增加算法的整体计算量，从而实现三维地震观测系统多频率偏移噪声的快速获取方法。该发明大幅度提升多频率偏移噪声计算的速度，有力的促进了偏移噪声分析技术在实际地震勘探中的应用。

相对于现有技术，本发明应用实例技术方案不需要对每个频率单独进行一次偏移噪声分析运算，再对所有的频率成分进行求和得到多频率下结果，效率较高，适宜在石油勘探领域的推广和应用。

本领域技术人员还可以了解到本发明实施例列出的各种说明性逻辑块(illustrativelogical block)，单元，和步骤可以通过电子硬件、电脑软件，或两者的结合进行实现。为清楚展示硬件和软件的可替换性(interchangeability)，上述的各种说明性部件(illustrativecomponents)，单元和步骤已经通用地描述了它们的功能。这样的功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用，可以使用各种方法实现所述的功能，但这种实现不应被理解为超出本发明实施例保护的范围。

本发明实施例中所描述的各种说明性的逻辑块，或单元都可以通过通用处理器，数字信号处理器，专用集成电路(ASIC)，现场可编程门阵列或其它可编程逻辑装置，离散门或晶体管逻辑，离散硬件部件，或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器和微处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或任何其它类似的配置来实现。

本发明实施例中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软件模块、或者这两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其它任意形式的存储媒介中。示例性地，存储媒介可以与处理器连接，以使得处理器可以从存储媒介中读取信息，并可以向存储媒介存写信息。可选地，存储媒介还可以集成到处理器中。处理器和存储媒介可以设置于ASIC中，ASIC可以设置于用户终端中。可选地，处理器和存储媒介也可以设置于用户终端中的不同的部件中。

在一个或多个示例性的设计中，本发明实施例所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的媒介上，或以一个或多个指令或代码形式传输于电脑可读的媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其它地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但不限于RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁性存储装置，或其它任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或通用或特殊处理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例如，如果软件是从一个网站站点、服务器或其它远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线方式传输的也被包含在所定义的电脑可读媒介中。所述的碟片(disk)和磁盘(disc)包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述的组合也可以包含在电脑可读媒介中。

以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种三维地震观测系统多频率偏移噪声获取方法，其特征在于，所述三维地震观测系统多频率偏移噪声获取方法包括：

获取三维地震观测系统数据；

根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声；

对所述三维地震观测系统的单频率偏移噪声进行插值，获取所述三维地震观测系统的单频率偏移噪声之外的一种或多种频率下偏移噪声；

对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声。

2.如权利要求1所述三维地震观测系统多频率偏移噪声获取方法，其特征在于，所述三维地震观测系统数据包括：震源点的三维空间坐标、检波点三维空间坐标、震源点与检波点之间的对应关系；其中，所述震源点的三维空间坐标包括：高程；所述检波点三维空间坐标包括：高程。

3.如权利要求1所述三维地震观测系统多频率偏移噪声获取方法，其特征在于，所述根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声，包括：

对于均匀介质，利用如下公式计算单频率偏移噪声：

$B_D\left(f\right)\≈\frac{1}{4{\mathrm{\π}}^2}\∫\∫\exp \left[i(k_x\mathrm{\Δx}+k_y\mathrm{\Δy})\right]D^{\′}{\mathrm{dk}}_x{\mathrm{dk}}_y,$

其中，i为虚数单位，k_x和k_y分别为水平x方向和水平y方向波数，Δx和Δy分别为检波点到地下目标反射点x方向和y方向的水平距离分量，D′为检波点矩阵D的波数域变换，B_D(f)为单频率偏移噪声。

4.如权利要求1所述三维地震观测系统多频率偏移噪声获取方法，其特征在于，所述根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声，包括：

对于非均匀介质，利用如下公式计算单频率偏移噪声：

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)为单频率偏移噪声。

5.如权利要求3或4所述三维地震观测系统多频率偏移噪声获取方法，其特征在于，所述对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声，包括：

利用如下公式计算所述三维地震观测系统的多频率偏移噪声：

$B\≈\sqrt{\mathrm{\Δf}\underset{j=1}{\overset{N}{\mathrm{\Σ}}}\left[j^2{\left|w\left(\mathrm{j\Δf}\right)B(\frac{\mathrm{j\Δx}}{N},\frac{\mathrm{j\Δy}}{N},\mathrm{N\Δf})\right|}^2\right]},$

其中，B为多频率偏移噪声，N是频率采样点数，Δf为频率采样间隔，w(f)为雷克子波频率域表达式 $w\left(\mathrm{j\Δf}\right)=\frac{2}{\sqrt{\mathrm{\π}}}\frac{{\left(\mathrm{j\Δf}\right)}^2}{f_p^3}\exp \left[\frac{-{\left(\mathrm{j\Δf}\right)}^2}{f_p^2}\right],$ f_p为地震波主频。

6.一种三维地震观测系统多频率偏移噪声获取装置，其特征在于，所述三维地震观测系统多频率偏移噪声获取装置包括：

数据获取单元，用于获取三维地震观测系统数据；

单频率偏移噪声单元，用于根据所述三维地震观测系统数据，计算所述三维地震观测系统的单频率偏移噪声；

插值单元，用于对所述三维地震观测系统的单频率偏移噪声进行插值，获取所述三维地震观测系统的单频率偏移噪声之外的一种或多种频率下偏移噪声；

多频率偏移噪声单元，用于对所述三维地震观测系统的单频率偏移噪声和所述单频率偏移噪声之外的一种或多种频率下偏移噪声进行求和，获取所述三维地震观测系统的多频率偏移噪声。

7.如权利要求6所述三维地震观测系统多频率偏移噪声获取装置，其特征在于，所述三维地震观测系统数据包括：震源点的三维空间坐标、检波点三维空间坐标、震源点与检波点之间的对应关系；其中，所述震源点的三维空间坐标包括：高程；所述检波点三维空间坐标包括：高程。

8.如权利要求6所述三维地震观测系统多频率偏移噪声获取装置，其特征在于，所述单频率偏移噪声单元对于均匀介质，利用如下公式计算单频率偏移噪声：

$B_D\left(f\right)\≈\frac{1}{4{\mathrm{\π}}^2}\∫\∫\exp \left[i(k_x\mathrm{\Δx}+k_y\mathrm{\Δy})\right]D^{\′}{\mathrm{dk}}_x{\mathrm{dk}}_y,$

其中，i为虚数单位，k_x和k_y分别为水平x方向和水平y方向波数，Δx和Δy分别为检波点到地下目标反射点x方向和y方向的水平距离分量，D′为检波点矩阵D的波数域变换，BD(f)为单频率偏移噪声。

9.如权利要求6所述三维地震观测系统多频率偏移噪声获取装置，其特征在于，所述单频率偏移噪声单元对于非均匀介质，利用如下公式计算单频率偏移噪声：

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)为单频率偏移噪声。

10.如权利要求8或9所述三维地震观测系统多频率偏移噪声获取装置，其特征在于，所述多频率偏移噪声单元利用如下公式计算所述三维地震观测系统的多频率偏移噪声：

$B\≈\sqrt{\mathrm{\Δf}\underset{j=1}{\overset{N}{\mathrm{\Σ}}}\left[j^2{\left|w\left(\mathrm{j\Δf}\right)B(\frac{\mathrm{j\Δx}}{N},\frac{\mathrm{j\Δy}}{N},\mathrm{N\Δf})\right|}^2\right]},$

其中，B为多频率偏移噪声，N是频率采样点数，Δf为频率采样间隔，w(f)为雷克子波频率域表达式 $w\left(\mathrm{j\Δf}\right)=\frac{2}{\sqrt{\mathrm{\π}}}\frac{{\left(\mathrm{j\Δf}\right)}^2}{f_p^3}\exp \left[\frac{-{\left(\mathrm{j\Δf}\right)}^2}{f_p^2}\right],$ f_p为地震波主频。