CN102667528B - 结合来自传感器的地震数据以衰减噪声 - Google Patents

Abstract

为了针对地震勘探进行噪声衰减，将传感器组件部署在地表上，其中传感器组件具有测量经地下结构传播的地震波的地震传感器和包括测量噪声的压力传感器的发散传感器。从地震传感器接收第一数据，从发散传感器接收第二数据。将第一数据和第二数据结合以衰减第一数据中的噪声。

Description

结合来自传感器的地震数据以衰减噪声

背景技术

地震勘探用来验明例如油气储层、淡水含水层、注气区等的地下元素。在地震勘探中，震源被放置在地表或海底上的不同位置，启动震源来产生指向地下结构的地震波。

震源产生的地震波传播到地下结构中，一部分地震波被反射回地面以便被地震接收器(例如地震检波器、加速度仪等)接收。这些地震接收器产生表示检测到的地震波的信号。对来自地震接收器的信号进行处理以生成关于地下结构的内容和特性的信息。

典型的基于陆地的地震勘探布置包括将地震接收器阵列部署在地面上，其中地震接收器以近似网格形式设置。这些地震接收器可以是能够在三个正交方向(竖直z、水平纵测线x和水平横测线y)上测量引入的波场的多分量地震检波器。

对于基于陆地的地震勘探来说，可能存在多种不想要的波场，包括地滚波噪声，例如瑞利或洛夫地面波。这些不想要的波场会污染由地震接收器采集的地震数据。尽管存在多种从地震数据去除不想要的波场的传统技术，但这种技术比较复杂并且可能成本不菲。

发明内容

总的来说，根据一种实施方式，一种用于地震勘探的噪声衰减方法包括在陆地表面上部署传感器组件，其中传感器组件具有测量经地下结构传播的地震波的地震传感器和包括测量噪声的压力传感器的发散传感器。将从地震传感器接收的第一数据和从发散传感器接收的第二数据结合以衰减第一数据中的噪声。

其他或替代的特征通过下面的描述、附图和权利要求将变得更加明白。

附图说明

图1表示根据一种实施方式可以部署成执行基于陆地的地震勘探的传感器组件的实例性布置的示意图；

图2表示根据一种实施方式可以在图1的布置中采用的传感器组件；

图3A-3B是表示根据一种实施方式可以由传感器组件检测的波场的传播的图解；

图4A-4C是表示根据一些实施方式由传感器组件采集的时间偏移域中的数据的图解；以及

图5是根据一种实施方式的执行地震勘探的过程的流程图。

具体实施方式

在下面的描述中，阐述了很多细节来提供对本发明的理解。然而，本领域技术人员将理解，本发明可以在没有这些细节的情况下实践，且来自所描述的实施方式的很多变型或修改是可行的。

根据一些实施方式，为了衰减基于陆地的勘探布置中的地震数据中的噪声，采用均具有至少一个地震传感器和至少一个发散传感器(用于测量噪声)的传感器组件。在一些实施方式中，发散传感器使用填充有其中设置压力传感器(例如水听器)的物质的容器形成。这种布置的压力传感器能够记录主要的噪声，使得来自传感器组件中的压力传感器的数据可用来开发用于净化由地震传感器采集的地震数据的噪声参考模型。其中浸没有压力传感器的物质可以是液体、胶体或例如沙子或塑料的固体。

一种噪声是地滚波噪声。地滚波噪声指的是震源产生的沿着地面朝着地震接收器基本上水平地传播的地震波。这些水平传播的地震波(例如瑞利波或洛夫波)是会污染地震数据的不希望的分量。通常，“噪声”指的是地震数据中不希望的任何信号分量(例如代表从地下元素反射的信号的数据)。其它类型的噪声包括在例如水体或永久冻土体的冰冻表面上采集的数据中存在的弯曲波、和由因诸如风、雨或像交通的人类活动、空气冲击、耀斑噪声或其它工业过程的环境造成的空气噪声。

图1是用于基于陆地的地震勘探的传感器组件100的布置的示意图。传感器组件100部署在地面108(以排或阵列的形式)上。传感器组件100在地面“上”的意思是传感器组件100设置在地面上面或上方或者埋置(完全或部分地)在地面下使得传感器组件100在地面附近10米内。地面108在含有感兴趣的至少一种地下元素106(例如油气储层、淡水含水层、注气区等)的地下结构102的上方。一个或多个震源104(可以是振动器、气枪、爆炸装置等)部署在其中定位有传感器组件100的勘探场中。

启动震源104引起地震波传播到地下结构102中。替代地，代替使用上面所述的受控震源来提供受控源或有源勘探，一些实施方式也可在无源勘探的环境下使用。无源勘探使用传感器组件100来执行下面的一项或多项：(微)地震监控；水压裂监控，其中由于主动注射到地下(例如油气储层)的流体造成的岩石失效而观察到微地震；等等。从地下结构102反射(以及从感兴趣的地下元素106反射)的地震波向上朝着传感器组件100传播。相应传感器组件100中的地震传感器112(例如水听器、加速度仪等)测量从地下结构102反射的地震波。而且，传感器组件100还包括发散传感器114，其被设计为测量噪声，例如地滚波噪声或其它类型的噪声。来自发散传感器114的数据可用来开发噪声参考模型以便衰减测量的地震信号中的噪声。

在一种实施方式中，传感器组件100通过电缆线110互连到控制器116。替代地，代替通过电缆线110连接传感器组件100，传感器组件100可以与控制器116无线连通。在一些实施中，可以将中间路由器或集线器设置在传感器组件100的网络的中间点处，以使得传感器组件100和控制器116之间连通。

图1所示的控制器116还包括能够在处理器122上执行的处理软件120。处理器122连接到存储介质124(例如一个或多个基于盘的存储设备和/或一个或多个存储器设备)。在图1的例子中，存储介质124用来存储从传感器组件100的地震传感器112连通到控制器116的地震传感器数据126，并且存储从传感器组件100的发散传感器114连通的发散数据128。

在操作中，软件120用来处理地震传感器数据126和水听器数据128。使用下面进一步讨论的技术，将水听器数据128与地震传感器数据126结合来衰减地震传感器数据126中的噪声(以产生净化形式的地震传感器数据)。软件120接着可产生输出以基于净化的地震传感器数据126描述地下结构102的特征。

根据一些实施方式的传感器组件100在图2中更详细地描述。传感器组件中的地震传感器112可以是用于测量由地下结构102中的地震波引起的粒子速度的地震检波器，或者替代地，地震传感器112可以是用于测量由经过地下结构102传播的地震波引起的加速度的加速度仪。

在一些实施方式中，地震传感器112是用于测量竖直方向(图1中的轴线z表示的)上的地震波的竖直分量地震传感器。在替代实施方式中，传感器组件100可另外或替代地包括用于检测基本水平方向、例如与地面108基本平行的x或y方向上的地震波的地震传感器。

也是传感器组件100的一部分的发散传感器114(在传感器组件100的壳体101内)用于测量与地震传感器112测量的经过地下结构102传播的地震波不同的输入(例如噪声)。在替代实施方式中，传感器组件100的发散传感器114可以在物理上与地震传感器112间隔开一些预定距离。

发散传感器114具有密封的封闭容器200。容器200容纳容器200内的一定体积的液体202(或其它物质，例如胶体或像沙子或塑料的固体)。而且，容器200容纳浸没在液体202(或其它物质)中的水听器204(或者其它类型的压力传感器)。压力传感器浸没在物质中的意思是压力传感器被所述物质包围或者以其他方式附接到物质或与物质接触。在下面的讨论中，以浸没在液体202中的水听器204为参考，注意在替代实施方式中，其他类型的压力传感器可以浸没在其它类型的物质中。中立且悬浮地浸没在液体202中的水听器204在机械上与容器200的壁断开。因此，水听器204只对经容器200的壁感应到液体202中的声波敏感。为了保持固定的位置，水听器204通过联接机构206附接，该联接机构206抑制声波经其传播。

液体202的例子包括：煤油、矿物油、植物油、硅油和水。在其它实施方式中，可以采用其他类型的液体。可以使用具有更高粘性的液体来改变对不同类型的波的敏感度，包括P(压缩)波、S(剪切)波、瑞利波和洛夫波。而且，设置在发散传感器114的容器200中的液体202的量决定水听器204的敏感度。只部分填充液体的容器200记录较弱的信号。在一些实施方式中，容器200可部分填充液体以在容器200内提供扩张体积。液体202的扩张，例如由于液体202的温度升高，可以在扩张体积(填充有气体)中得到适应。

如图2进一步所示，传感器组件100还包括与地震传感器112和发散传感器114都电联接的电路208。电路208可包括存储元件、处理元件和用于在电缆线110上将由地震传感器112和发散传感器114采集的数据连通到控制器116的通信元件(图1)。

如图2所示，地震传感器112定位在发散传感器114的容器200的上方和外部。替代地，地震传感器112可具有相对于发散传感器114的一些其他布置。发散传感器114的至少一部分在地面108下，使得水听器204可以在地面108处或者在地面108下方，而不是在地面108上方。当安置时，传感器组件100的发散传感器114与地面108下的土地介质稳固地接触，这改善了通过发散传感器114中的水听器204采集的信号的数据质量。

在采用缆线110的实施方式中，经缆线110从远程电源(例如位于控制器116处的电源)向传感器组件100传输功率。在采用无线连通而不使用缆线110的实施方式中，传感器组件100可设有提供本地功率的电池。

在基于陆地的地震勘探中，粒子位移(或速度或加速度)通过自由表面(地面108)正下方的地震传感器测量。因此，观测到的信号分量不仅包含即将到来的冲击地震波，而且包括在固体-空气界面(地面108)处向下反射/转变的波的额外贡献。图3A表示了引入的压缩波或P波的例子，图3B表示引入的剪切波或S波的例子。P波和S波从地下结构102引入(图1)。P波在地震波的传播方向上延伸，而S波在基本垂直于地震波传播方向的方向上延伸。

在图3A和3B的每个中，z方向代表竖直方向，而x方向代表纵测线水平方向。空气-固体界面对应于图1中所示的地面108。即将到来的P波表示为P_up(图3A)，而即将到来的S波表示为S_up(图3B)。地震传感器112不仅记录即将到来的P波或S波，而且记录在固体-空气界面处向下反射/转变的波的额外贡献，包括响应于P_up波反射的R_up和R_ps波(图3A)。类似地，从S_up波反射的反射/转变的波表示为图3B中的R_sp和R_ss。

R_pp波是从P_up波反射的P波，而R_ps波是从P_up波反射的S波。类似地，R_sp波是从S_up波反射的P波，而R_ss波是从S_up波反射的S波。

与地震传感器112相反，图2所示的传感器114中的填充有液体202的容器200内的水听器204对波的传播方向和角度不敏感。结果，这导致近乎竖直入射角(小慢度)处事件的破坏性总和，并且因此，水听器204的测量值(表示为U_H)记录主要的表面噪声(其具有大慢度)。慢度与视速度的倒数成比例：小慢度由高视速度造成，而大慢度由小视速度造成。因此，U_H(来自图2的发散传感器114中的水听器204的数据)比U_X(由地震传感器测量的地震数据的x分量)提供更好的噪声表示。因此，基于自适应减法和/或极化，U_H可用于在U_Z上去除噪声。U_Z指的是在z方向上测量的地震波。在下面的讨论中，假设U_X、U_Z和U_H测量自由表面正下方的粒子位移或速度或加速度或压力。

U_Z和U_X根据公式1和2表示如下：

U_z＝(-q_αα+R_ppq_αα-R_pspβ)P_up+(pβ-R_sspβ+R_spq_aa)S_up    (公式1)

U_x＝(pα+R_pppα+R_psq_ββ)P_up+(q_ββ+R_ssq_ββ+R_sppα)S_up   (公式2)

其中P_up和S_up是入射P波和S波(分别如图3A和3B所示)，α和β是接近地面的P波速度和S波速度，且P＝sini/α＝sinj/β是水平慢度。对于P波和S波的竖直慢度分别是q_α＝(α^-2-p²)^0.5和q_β＝(β^-2-p²)^0.5。R_ij项是入射i波在传感器正上方的固体-空气界面处向后反射/转变成j波的反射/转变系数。公式1和2各自的左手侧涉及考虑到图3A所示的入射P波的自由表面效应，而公式1和2各自的右手侧涉及考虑图3B所示的入射S波的自由表面效应。注意，地震传感器(地震检波器或者加速度仪)信号分量是矢量测量-例如，如果即将到来的P波情形被认为在z分量上，则公式1和2各自的左手侧中的+和-涉及传播方向，即上行(-)或下行(+)，而项q_αα(cosi)或Pβ(sinj)涉及传播角度，其中i和j分别是P波和S波(相对于竖向)的传播角度。

与以矢量方式记录P波和S波两者的地震传感器112(图2)相比，液体单元中的水听器(图2中的204)只记录P波场因为S波传播不受液体支持导致的压力波动。另外，假定该传感器(发散传感器114)对传播方向不敏感，因为压力是标量。因此，与公式1和2表示的地震传感器数据相比，对于完全各向同性的水听器传感器，角度-方向有关的项消失，且水听器数据U_H可表示为：

U_H＝C₁(P_UP+R_PPP_up+R_SPS_up)            (公式3)

其中C₁是常数，包括：(1)补偿水听器204(测量压力)和地震传感器112(测量位移或速度或加速度)之间的敏感度差异的校正系数，和(2)从土地到发散传感器114的容器的传输系数。水听器204不直接记录S波，但注意由于自由表面处的S-到-P转变(公式3的最后项)，S波场有关的能量将存在于水听器数据中。最后，通过将作为U_Z和U_X的函数的P_up和S_up表示和用于反射系数的表示包括在公式3中，可以获得：

U_H＝C₂PU_x        (公式4)

其中C₂取决于C₁和接收器附近的近地面属性。

在慢度域中，U_H仅仅是U_X的依赖P(P表示慢度)的换算形式。这意味着，在小慢度(小P)处，或者等同地在小入射角处，甚至与U_X相比，U_H分量上的波幅是极其小的。对于S波以及P波来说也是如此。U_H上的这种自然P波衰减源自于上行和下行反射事件之间的破坏性总和，在竖直入射(p＝0)处R_pp等于-1。总而言之，通过比大慢度(慢表观波，通常为地滚波噪声)多得多地衰减小慢度(像反射那样的高视速度)，水听器204用作自然速度过滤器。

通过将公式4(慢度域中的)改写成传统时间偏移域(P＝δx/δt)中的，获得：

$\frac{\mathrm{\δ}{U}_H(t,x)}{\mathrm{\δt}}=C_2\frac{\mathrm{\δ}{U}_x(t,x)}{\mathrm{\δ\χ}},$ (公式5)

上面的纵测线(2D)情况可以延伸到3维情况：

$\frac{\mathrm{\δ}{U}_H}{\mathrm{\δt}}\∝(\frac{{\mathrm{\δU}}_x}{\mathrm{\δx}}+\frac{\mathrm{\δ}{U}_y}{\mathrm{\δy}}),$ (公式6)

公式6表示的水听器信号分量U_H的时间导数与波场(自由表面正下方)的散度成比例。

图4A-4C表示时间偏移域中的模拟数据，其中偏移指每个传感器组件100和震源之间的距离。图4A表示针对U_Z的时间偏移模拟数据，图4B表示针对U_X的时间偏移模拟数据，图4C表示针对U_H的时间偏移模拟数据。图4B-4C中的每个中的大体锥形的图样402代表噪声。注意，噪声出现在每个U_Z、U_X和U_H中。然而，实际地震数据(通过图4A-4C中的弯曲结构表示)在U_H中特别在小偏移处衰减。因此，很明显，输出U_H的水听器204提供自然速度过滤。

与地滚波锥形402相比，体波以及特别是在小偏移处反射的P波在U_H上比在U_X上具有甚至更小的幅度。因此，针对U_Z上的地滚波衰减，U_H比U_X提供更好的噪声参考模型。另外，U_H是比U_X(或者更一般地与地震传感器记录的地震数据正交或接近正交的分量)更好的噪声模型，因为U_H还记录离线发散事件(与U_X相比，存在于U_Z和U_y)。

U_H和U_Z之间的全局地滚波噪声关联比U_X和U_Z之间的好，因为U_H上的S波和洛夫波比U_X上的弱。结果，基于自适应减法或极化过滤，U_H可用来提供U_Z上的更好噪声衰减。通过采用U_H代替U_X，第二信号保留应当被改善，因为U_H含有较少的体波能量(特别是在地滚波噪声占主导的小慢度偏移处)。

加权的HZ总和(U_H和U_Z信号分量的总和)能够使得噪声从U_Z去除。基本的自适应减法可表示为：

$U_z^{\mathrm{clean}}(t,x)=U_z(t,x)-F(t,x).U_H(t,x),$ (公式7)

其中F可以是例如通过在变化的时间偏移窗口中匹配U_H与U_Z(并且例如在至少二次方意义上)估计的标量或频率依赖算子。在上面的公式7中，表示噪声分量被去除后的地震数据。

F(t，x)是算子，一般情况下它可以是几点长度的微波，或者在滤波长度＝1的特殊情况下简单地为一个数(标量)。在后一情况下，数F简单地通过在选择的时间偏移窗口中用U_Z除U_H得到：F(t，x)＝Z(x，t)/H(x，t)。F不依赖于频率，因为它仅仅是U_H和U_Z之间的换算系数(同样是在具体的时间偏移窗口中的)。

更一般地，当考虑一般情况时，F(t，x)可以称为维纳滤波或者U_H和U_Z之间的转换函数，并且F(t，x)变为依赖频率。在这种情况下，F(t，x)评估基于Z(x，t)和H(x，t)的相互关联和自关联的使用。很多其他算法可用来评估F，例如在L1范数或L2范数和最小二乘方中求解F。其他设计和应用滤波算子F的合适方法包括极化、自适应、静态、确定、多渠道和多维过滤技术。

尽管已经参考了净化U_Z地震信号(z方向上的地震信号分量)，注意发散传感器输出U_H也可以用来净化水平方向上的地震信号分量U_X或U_y。

通过采用根据一些实施方式的噪声衰减技术，可以增加传感器组件之间的间距。例如，相邻传感器组件之间的间距可以设成使得传感器组件彼此以大于由发散传感器记录的噪声的最短波长的一半的距离间隔开。结果，需要部署较为稀疏的传感器组件阵列，这降低了设备成本，并降低了与传感器组件在场中的部署相关的人工成本。

图2表示具有基本上立方体容器200的发散传感器114。容器200的形状可改为其他形状，例如下列形状：平行六面体形状、金字塔形状、四边棱锥形、双棱锥形、椭圆形和球形。改变容器200的形状引入了各向异性，这可产生波场分解机会(例如分解成P和S波)。

代替将水听器204(图2)浸没在液体中，水听器204可以替代地浸没在更具粘性的流体(例如硅油)、胶体或者甚至是例如沙子或塑料的固体中，以允许检测S波。这可允许波场分解，例如P波与S波或者瑞利波与洛夫波。

图5表示根据一种实施方式用于执行地震勘探的总过程。部署传感器组件100(502)以便进行地震勘探。接着启动一个或多个震源(例如图1中的104)(504)。响应于震源的启动，地震波传播到地下结构102(图1)。反射的地震波通过传感器组件100检测。

地震数据(由地震传感器112测量的)和发散数据(由发散传感器114测量的)通过控制器116接收(506)。控制器116接着将地震数据和发散数据结合(508)以例如根据公式7净化地震数据。被净化的地震数据接着通过控制器116处理(510)以描述地下结构102的特征。注意地震测量和发散测量的结合可以在从传感器接收数据后通过中央控制单元执行，或者在图2中的传感器单元中的电子器件内侧或者在任何中间位置本地地进行。所述结合可以在数据处理中心实时地、在每次放炮后、完成勘探后或稍后进行。

根据一些实施方式的噪声衰减技术或系统可在具有障碍的勘探布置中采用，这些障碍例如阻止传感器组件以规则间距设置的建筑物或自然障碍。换句话说，所述障碍干扰传感器组件的规则形式。通过使用根据基于发散传感器使用的一些实施方式的噪声衰减技术，可以得到更好的结果。

上面描述的软件(包括图1的软件120)的指令装载成在处理器(诸如图1中的处理器122)上执行。处理器包括微处理器、微控制器、处理器模块或子系统(包括一个或多个微处理器或微控制器)、或其他控制或计算设备。“处理器”可指单个部件或者多个部件(例如，一个CPU或多个CPU)。

数据和指令(软件的)存储在相应的存储设备中，作为一个或多个计算机可读或计算机可使用的存储介质实施。存储介质包括不同形式的存储器，包括：例如动态或静态随机存取存储器(DRAM或SRAM)、可擦除和可编程只读存储器(EPROM)、可电擦除和编程只读存储器(EEPROM)和闪存的半导体存储设备；例如固定、软盘和可移动盘的磁盘；包括磁带的其他磁性介质；以及诸如光盘(CD)或数字视频光盘(DVD)的光学介质。

尽管已经关于有限数量的实施方式公开了本发明，但受益于此披露的本领域技术人员将认识到对其进行的多种修改和变型。所附的权利要求意在覆盖落入本发明的真正精神和范围内的这种修改和变型。

Claims (23)

1.一种用于陆上地震勘探的噪声衰减方法，包括：

在地表上部署传感器组件，其中所述传感器组件具有测量经地下结构传播的地震波的地震传感器以及包括测量噪声的压力传感器的发散传感器；

从所述地震传感器接收第一数据，从所述发散传感器接收第二数据；以及

结合所述第一数据和所述第二数据以衰减所述第一数据中的噪声。

2.根据权利要求1所述的方法，其中所述地震传感器和所述发散传感器在物理上间隔开预定距离。

3.根据权利要求1所述的方法，其中结合所述第一数据和所述第二数据包括从所述第一数据减去所述第二数据。

4.根据权利要求3所述的方法，其中从所述第一数据减去所述第二数据包括从所述第一数据减去所述第二数据和滤波器算子的乘积。

5.根据权利要求1所述的方法，其中所述发散传感器包括容纳物质和浸没在所述物质中的所述压力传感器的容器。

6.根据权利要求5所述的方法，其中所述物质选自包括液体、胶体和固体的组。

7.根据权利要求1所述的方法，其中来自所述发散传感器的所述第二数据提供比与来自所述地震传感器的所述第一数据正交或者接近正交的分量更好的噪声模型。

8.根据权利要求1所述的方法，其中通过所述发散传感器测量的噪声包括地滚波噪声。

9.根据权利要求1所述的方法，其中通过所述发散传感器测量的噪声包括弯曲波噪声和空气噪声中的一种。

10.根据权利要求1所述的方法，还包括在所述地表上部署另外的传感器组件，其中所述另外的传感器组件中的每个具有测量经地下结构传播的地震波的地震传感器和测量噪声的发散传感器。

11.根据权利要求10所述的方法，其中部署所述传感器组件包括在具有干扰传感器组件的规则形式的一个或多个障碍的环境中部署所述传感器组件，其中所述发散传感器的设置甚至能够在不具有传感器组件的规则形式下实现噪声衰减。

12.根据权利要求11所述的方法，其中部署所述传感器组件包括使得传感器组件彼此以大于由所述发散传感器记录的噪声的最短波长的一半的距离分开设置。

13.一种系统，包括：

传感器组件，其用于部署在地表上以便执行地震勘探；和

控制器，其具有接收由所述传感器组件收集的数据的处理器，其中每个所述传感器组件具有测量经地下结构传播的地震波的地震传感器和包括测量噪声的压力传感器的发散传感器，

其中所述处理器能够将来自所述地震传感器的第一数据与来自所述发散传感器的第二数据结合以衰减来自所述地震传感器的第一数据中的噪声。

14.根据权利要求13所述的系统，其中每个所述发散传感器包括填充物质和浸没在所述物质中的相应压力传感器的容器。

15.根据权利要求14所述的系统，其中所述物质选自包括液体、胶体和固体的组。

16.根据权利要求14所述的系统，其中在每个所述传感器组件中，所述地震传感器设置在所述发散传感器的所述容器的上方和外部。

17.根据权利要求13所述的系统，其中所述控制器能够在缆线上与所述传感器组件连通。

18.根据权利要求13所述的系统，其中所述控制器能够与所述传感器组件无线连通。

19.根据权利要求13所述的系统，其中将所述第一数据和所述第二数据结合包括从所述第一数据减去所述第二数据。

20.根据权利要求19所述的系统，其中从所述第一数据减去所述第二数据包括从所述第一数据减去所述第二数据和滤波器算子的乘积。

21.一种物品，其包括至少一个包含指令的计算机可读存储介质，在执行所述指令时引起处理器：

接收由布置在地表上的传感器组件收集的数据，其中每个所述传感器组件具有测量经地下结构传播的地震波的地震传感器和包括测量噪声的压力传感器的发散传感器；和

将来自所述地震传感器的第一数据与来自所述发散传感器的第二数据结合以衰减来自所述地震传感器的数据中的噪声。

22.根据权利要求21所述的物品，其中从所述发散传感器接收第二数据包括从均包括填充物质和浸没在所述物质中的压力传感器的容器的发散传感器接收数据。

23.根据权利要求21所述的物品，其中将所述第一数据和所述第二数据结合包括从所述第一数据减去所述第二数据。