使用被动地震信号表征一个区域的下方地层的方法以及相应的系统与流程

本发明涉及表征一个区域的地表下面的一般领域，这特别是通过研究被动地震信号(passiveseismicsignal)，并且特别是通过研究低频被动地震信号。

背景技术：

低频被动地震信号是在0.1赫兹(hz)到10hz或者实际上是0.1hz到4hz或5hz范围内的频率处自然出现的代表地层移动的信号。在这些频率范围内，地层的移动可以与产生在海平面上方的陆地内部传播的波的海洋中的波的移动相关联。高于1hz的频率处的波还可以由人类活动(道路、工业)生成。

应当理解的是，在本文中使用的术语“被动”意味着并未由用户或者由用户控制的工具生成地震信号。

还应当观察到，前面所提到的频率范围不同于高频地震信号所涉及的频率范围，其中高频地震信号可以处在10hz到150hz的范围内，并且特别与观测诱发地震有关。

文献“phenomenologyoftremor-likesignalsobservedoverhydrocarbonreservoirs(在油气储层上观测到的类震动信号的现象学)”(s.dangel等人，journalofvolcanologyandgeothermalresearch)描述了在包含油气(因此包含多个流体或者相同流体的多个相)的地下储层上方将观测到处在信号频谱中的1.5hz到4hz范围内的频率峰值。其结果是，观测到这样的信号使得有可能确定在一个区域内的地表下面是否存在多个流体或者单一流体的多个相。

在这样的情况下，通过研究低频被动地震信号，有可能跟踪被用于储存油气(例如天然气)或蒸汽或各种类型的气体(例如co2、h2)的储层中的变化，或者实际上有可能在地热田中实施勘探油气或者包含水和蒸汽的蓄水层的操作。

在现有技术中，以下文献是已知的：

文献us2008/0021656描述了一种处理地震数据的方法，其中采集被动地震信号，获得这些信号的频谱，计算频谱的垂直分量与频谱的水平分量之间的比值，并且对所述比值进行积分。为了实施在该文献中描述的方法，必须选择将在其中实施处理的频率范围。该解决方案不够灵活，并且利用该方法获得的结果也不够准确。

文献wo2010/080366公开了一种通过使用被动地震数据与某种其他类型(例如主动地震数据)的地球物理数据相组合来检测油气的方法。

文献wo2009/027822提出一种确定油气储层的位置的方法，其中获得地震数据，从信号获得频谱，并且确定这些频谱的极大值以便获得地图。

文献ep2030046也描述了一种其中对频谱的幅度之间的比值进行研究的方法。该文献还描述了对频谱进行平滑。

文献ep1960812和文献ep2030046描述了这样的方法，其中获得地震数据，从地震数据获得频谱，并且对频谱进行平滑。

文献wo2009/081210描述了这样一种方法，其中确定通过测量地表下方的移动而获得的频谱的一个频带内的能量。

文献us2014/0254319描述了采集被动地震信号。该文献提出对于0到5hz范围内的频率使用低通滤波器，因为正是在这些频率处基于被动地震信号观测到现象。在该文献中，信号被变换到频率-波数域。

文献wo2014/108843描述了这样一种方法，其中采集微地震信号并且对其进行卷积，随后应用滤波器。

在这些文献中描述的解决方案并不令人满意。这些解决方案没有使得以足够准确的方式确定在一个区域的地表下面是否存在不同的流体或相成为可能。具体来说，这些方法对于人为噪声非常敏感，并且可能难以在其中存在流体的区域与其他区域之间准确地作出区分。

某些方法需要事先确定将要观测的频率范围。该解决方案并不令人满意，因为其没有考虑到出现在将要观测的频率范围之外的异常。

本发明特别试图减轻这些缺点。

技术实现要素：

本发明通过提出一种表征一个区域的地表下面的方法而满足了前述需求，所述方法包括以下步骤：

-准备多个频谱，所述频谱示出在所述区域的地表附近的所述区域的至少一点处获得的被动地震信号的频谱密度，其中在该点处对被动地震信号进行记录，每一个频谱是从示出移动的信号准备的；

-对于出现在每一个频谱中的每一个频谱确定至少一个频谱属性，从而获得与各条记录并且与各个频率相关联的频谱属性的集合；

-把所述属性集合组织在矩阵中，其中每一行与一条记录相关联；以及

-对所述矩阵应用主成分分析方法以便确定主成分，从而推断出所述地表下面的特性。

可以观察到，移动可以被感测为移动速度信号或者实际上被感测为移动加速度信号。优选的是使用移动速度信号。还可以观察到，移动可以在三个方向上发生(一个垂直方向和两个水平方向)，并且可以准备仅与垂直移动相关联的频谱以及/或者仅与水平移动相关联的频谱。

所得到的矩阵的规格是记录的数目乘以对应于每一个属性的频率的数目(每个属性的频率数目对于不同的属性可以是不同的)。

发明人观察到，通过使用主成分分析方法，获得了更容易地表现出对应于不同记录的频谱之间的差异的主成分。具体来说，主成分空间是用于表示频谱之间的差异的最佳空间。其结果是，该空间是在其中推断出所述区域的地表下面的特性的良好空间。

如果对多个点进行研究，并且通过使用主成分，有可能获得区域的图形表示，所述表示显示出下方流体的存在。如果只研究一个点，则通过在不同时刻采集的多个信号，只要流体的数量随着时间变化，就有可能观测到随着时间的变化。

本领域技术人员通常用其缩写提及主成分分析：pca。

因此可以理解的是，所获得的频谱是样本，并且其给出有限数目的频率。这些频率可以被选择在较宽的频率范围内，例如从0.1hz到4hz或5hz的频率范围。对于实施本发明不需要定义更窄的频率范围，但是在文献us2008/0021656所描述的方法中则有此必要。

还可以观察到，所述属性可以是表征频谱的频率的任何参数，并且当所述区域的点处在包含流体的区块的垂直上方时发生变化。因此，通过使用多个属性，本发明与文献us2008/0021656中所描述的解决方案的不同之处在于，所使用的仅有的属性是垂直与水平移动之间的比值。发明人观察到，通过使用多个属性，可以有效地检测流体的存在。

此外，在本发明中，可以通过地震仪获得被动地震信号，比如由加拿大供应商nanometrics在“t-40”商标名下销售的装置。这样的装置可以被埋设在区域地表的附近，例如埋设在大约五十厘米的深度。或者其可以被放置在地表上，前提是地表良好地耦合到地层。

在一种特定的实现方式中，所述移动是垂直移动和/或水平移动，并且对应于每一个频率的所述频谱属性属于从由以下各项构成的一组当中选择的类型：对应于垂直地震移动的频谱密度与对应于水平地震移动的频谱密度之间的比值，作为对应于水平地震移动的频率的函数的频谱密度的导数，以及作为对应于垂直地震移动的频率的函数的频谱密度的导数。

发明人观察到，通过组合多个这些频谱属性，有可能获得对于流体存在的有效检测。

在一种特定的实现方式中，通过围绕所选数目的频谱点应用线性回归来计算作为对应于水平地震移动的频率的函数的频谱密度的所述导数以及/或者作为对应于垂直地震移动的频率的函数的频谱密度的所述导数。

所选择的点可以通过将频率轴划分成0.5hz的跨距而获得。

在一种特定的实现方式中，从信号准备每一个频谱包括：

-把信号划分成全部具有相同的持续时间的多个相继的子信号；

-为每一个子信号准备一个频谱密度子频谱；

-对于各个子频谱的每一个频率，从对应于每一个子频谱中的该频率的频谱密度数值确定频谱密度的统计属性(例如中值)；以及

-获得将从所有频率的所有统计属性准备的所述频谱。

该特定实现方式用来获得频谱的良好平滑，这是因为使用了例如中值之类的统计属性。

在一种特定的实现方式中，每一个子信号在前一个子信号的至少一段非零持续时间上与前一个子信号重叠，例如前一个信号的50％的持续时间。

在一种特定的实现方式中，所述记录(或者可选的是所有所述记录)是在区域的不同点处进行的，所述频谱属性集合当中的每一个属性也与一点相关联。换句话说，属性集合当中的每一个属性与一个频率并且与一点相关联，这是因为记录本身与对于不同记录可能是不同的一点相关联。

在一种特定的实现方式中，所述记录是在预定的持续时间期间并且从预定的时间进行的。

优选的是选择处于半夜的起始时间(例如在午夜开始)和大约四个小时的持续时间；这样使得有可能在其中人为噪声处于其最低水平的时段期间采集信号。

在一种特定的实现方式中，所述记录是在同时进行的各组记录中进行的，其中每一组对应于在其间进行该组的记录的一天，各条记录是在所述区域的不同点处并且/或者从不同的时刻进行的。

该特定实现方式使用有限数目的测量装置对一个区域进行研究。在该特定实现方式中，每一天移动所述装置以便能够以良好的分辨率覆盖一个区域。

在一种特定的实现方式中，与给定的属性相关联的所述矩阵的各列全部是邻近的。该特定实现方式通过使用主成分分析使得某个属性的影响与另一个属性相比是更加可观测的。

在一种特定的实现方式中，每一组记录与矩阵的一个行组相关联，并且对于每一个行组，属性的数值被归一化。

还可以观察到，在矩阵中，可以通过将其相继置于共同的一组中而把这些行组分组在一起。

因此，由于每一组与一天相关联，归一化用来获得处于相同的数值范围内的属性数值，即使改变传感器的位置会在记录中采集到的信号的幅度中表现出改变。所述归一化可以是约化居中归一化。

在一种特定的实现方式中，所述主成分是投影器(projector)，并且所述矩阵被投影到每一个投影器上，从而对于每一个投影器获得所述区域的图形表示，其中显示出对应于每一点的矩阵投影的结果。

在一种特定的实现方式中，从所述投影器当中确定一定数目的k个投影器。换句话说，根据至少一条标准来选择投影器的数目。举例来说，该标准可以用来确定投影器给出关于区域的地表下面或者地表下面的随着时间的变化的良好图形表示。本领域技术人员知道如何根据通过除了构成本发明的一部分的那些手段之外的其他手段获得的关于地表下面的术语来评估该标准，或者确切来说知道如何根据关于区域中的地层结构的知识来评估该标准。

举例来说，如果图形表示是区域的灰度图，则可以通过确定在预期会有变化的位置处是否出现异常来确定k个投影器。作为指示，异常可以对应于地质陷阱(geologicaltrap)的结构形状。

该特定实现方式特别适合于研究从通过其他手段(井等等)实施的初始分析获知其范围的储层中的变化。

在一种特定的实现方式中，所述方法还包括：

-把所述矩阵投影到所述k个投影器上，从而对于所述矩阵的每一行获得一个长度为k的矢量；

-从所述各个长度为k的矢量获得第二矩阵；

-对第二矩阵应用被组织在一维或二维中的分类方法，以便获得n个行类；

-为矩阵的每一行分配至少一个类数(classnumber)从而表示区域的地表下面的异常的量值；

-为每一个行类准备类头(classhead)；

-从所述各个类头获得维度为n乘k的第三矩阵；以及

-对所述第三矩阵应用伪求逆方法以便获得第四矩阵，所述第四矩阵的维度是n乘以出现在频谱属性的初始矩阵的每一行中的频率的数目。

所述有组织的分类方法可以是本领域技术人员所熟知的方法，比如自组织映射(som)或生成式拓扑映射(gtm)。

n被选择成严格大于1，并且可以是大约十或几十的量级。

对于一维分类，为每一行分配一个类数，以便表示异常(也就是频谱中的变化)的量值。这些类数可以被表示在灰度上。不同的行可以具有相同的类数。

如果使用二维分类，则可以使用其他表示方法，例如通过组合两个标度的互补色。

可以观察到，对于每一个类数，有可能把类头定义为已为之分配类数的所有行的重心。或者，类头可以是最靠近重心的行，作为举例，这可以通过欧式度量确定。

或者，所述方法还包括：

-把所述矩阵投影到所述k个投影器上，从而对于所述矩阵的每一行获得一个长度为k的矢量；

-从所述各个长度为k的矢量获得第二矩阵；

-对所述第二矩阵应用伪求逆方法以便获得第三矩阵，所述第三矩阵具有与所述频谱属性的初始矩阵相同的维度；

-对第三矩阵应用被组织在一维或二维中的分类方法，以便获得n个行类；

-为矩阵的每一行分配至少一个类数从而表示区域的地表下面的异常的量值；

-为每一个行类准备类头；以及

-从所述各个类头获得第四矩阵，所述第四矩阵的维度是n乘以出现在频谱属性的初始矩阵的每一行中的频率的数目。

在该变型中，在应用分类方法之前实施伪求逆步骤。

本发明还提供一种用于表征一个区域的地表下面的系统，所述系统包括：

-用于准备多个频谱的模块，所述频谱表示在所述区域的地表附近的所述区域的至少一点处获得的被动地震信号的频谱密度，其中在该点处对被动地震信号进行记录，每一个频谱是从表示移动的信号准备的；

-用于对于出现在每一个频谱中的每一个频谱确定至少一个频谱属性的模块，从而适合于获得与各条记录并且与各个频率相关联的频谱属性的集合；

-用于把所述属性集合组织在矩阵中的模块，其中每一行与一条记录相关联；以及

-用于对所述矩阵应用主成分分析方法以便确定主成分的模块，从而推断出所述地表下面的特性。

所述系统可以被配置成实施所述方法的所有前面描述的实现方式。

本发明还提供一种包括指令的计算机程序，当所述程序由计算机执行时，所述指令用于执行如前面所定义的方法的各个步骤。

可以观察到，在本公开内容中所提到的计算机程序可以使用任何编程语言，并且可以具有源代码、对象代码或者介于源代码和对象代码之间的代码的形式，比如具有部分编译形式或者具有任何其他合乎期望的形式。

本发明还提供一种存储计算机程序的计算机可读数据介质，所述计算机程序包括用于执行如前面所定义的方法的各个步骤的指令。

在本公开内容中所提到的数据(或记录)介质可以是能够存储程序的任何实体或设备。举例来说，所述介质可以包括比如只读存储器(rom)之类的存储装置，例如紧致盘(cd)rom或微电子电路rom，或者实际上是磁性记录装置，例如软盘或硬盘。

此外，所述数据介质可以对应于可传输介质，比如适合于通过电缆或光缆、通过无线电或者通过其他手段传送的电信号或光学信号。本发明的程序可以特别从因特网类型的网络下载。

或者，所述数据介质可以对应于程序被合并在其中的集成电路，所述电路适于执行或者被使用在所讨论的方法的执行过程中。

附图说明

本发明的其他特性和优点将从后面参照附图进行的描述而变得显而易见，附图示出了一个不具限制性的实例。

在附图中：

-图1是示出本发明的一种实现方式中的方法的各个步骤的图示；

-图2是示出本发明的一个实施例中的系统的图示；

-图3是一个区域的地表下面的剖面图；

-图4是示出如何从信号获得频谱的图示；

-图5示出了属性被组织在其中的矩阵；

-图6示出了将矩阵投影到投影器上；

-图7是示出重投影矢量的分类的图表；以及

-图8示出了在把参数与每一个重投影矢量相关联之后所获得的图形表示。

具体实施方式

下面将描述根据本发明的一个具体实例的用于表征一个区域的地表下面的方法和系统。

图1是示出表征一个区域的地表下面的一种方法的各个步骤的图示。

所述方法可以被实施来确定在一个区域的地表下面是否存在多个流体或者单一流体的多个相。

使用这种方法的典型应用例如涉及监测包含油气(例如天然气)、蒸汽和各种类型的气体(例如co2、h2)的储层，勘探油气，在地热田中进行勘探。

在所述方法的第一步骤e01中，准备多个频谱，所述频谱示出在所述区域的地表附近的所述区域的至少一点处获得的被动地震信号的频谱密度，其中在该点处对被动地震信号进行记录，每一个频谱是从示出移动的信号准备的。

换句话说，事先实施采集示出水平移动(可能是不同方向上的两个信号)和/或垂直移动的信号的步骤。可以通过使用地震仪来采集这些信号，比如由加拿大供应商nanometrics在“t-40”商标名下销售的装置。这样的装置可以被规则地安排在一个区域的地表附近或者被安排在所述区域的地表上，正如后面参照图3所描述的那样，并且所述装置优选地在夜间被使用以便减少人为噪声。所有信号与对应的时刻和/或所研究的区域的对应的位置或点相关联。

可以观察到，前面提到的该类装置提供表示移动的移动速度信号。

每一个频谱可以通过确定信号的功率谱密度(psd)而从与之对应的信号获得。

还有可能实施试图平滑所获得的频谱的处理，正如后面参照图4所描述的那样。

在步骤e01之后获得频谱，或者实际上是包括与水平移动相关联的频谱并且还有与垂直移动相关联的频谱全部二者的频谱，这些频谱与记录相关联并且因此与其性质相关联，所述性质由区域的点以及采集的时刻或日期构成。

可以观察到，可以从表示两个不同方向上的移动的两个信号获得对应于水平移动的单一频谱，这是通过计算对应于每一个方向的频谱的几何平均。具体来说：

其中f是频率，fh(f)是对应于水平移动的频谱，psde是对应于第一方向上(具体来说是向东)的水平移动的频谱，并且psdn是对应于第二方向上(具体来说是向北)的水平移动的频谱。

因此，所述频谱被采样并且与被包含在先前所选择的较宽范围内的有限数目的频率相关。

在第二步骤e02中，确定频谱属性。这些属性可以从由以下各项构成的一组当中选择：对应于垂直地震移动的频谱密度与对应于水平地震移动的频谱密度之间的比值，作为对应于水平地震移动的频率的函数的频谱密度的导数，以及作为对应于垂直地震移动的频率的函数的频谱密度的导数。

这三个属性可以被如下写出：

以及

其中fv(f)是对应于垂直移动的频谱。

可以观察到，可以通过围绕所选数目的频谱点应用线性回归来计算作为对应于水平地震移动的频率的函数的频谱密度的导数以及/或者作为垂直地震移动的频率的函数的频谱密度的所述导数。

作为指示，所选择的点可以通过将频率轴划分成0.5hz的跨距而获得。

步骤e02用来获得一个属性集合，所述属性集合与频率相关联(所述频率在不同类型的属性之间可以是不同的)，与记录相关联并且因此与区域的点相关联，并且还与信号被采集的时刻和/或日期相关联。

在步骤e03中，该属性集合被组织到矩阵中，其中每一行与一条记录相关联(也就是说与区域的一个点并且与进行记录时的一个时刻和/或一个日期相关联)。

将参照图5更加详细地描述这一组织。

在步骤e04中，对所述成分应用主成分分析方法，以便确定主成分并且从而推断出所述地表下面的特性。

图2是示出适合于实施参照图1所描述的步骤e01到e04的系统1的图示。

系统1可以是计算机系统，并且其包括处理器2和存储器3。

计算机程序4的指令被存储在存储器3中。计算机程序4包括用于实施步骤e01的指令41，用于实施步骤e02的指令42，用于实施步骤e03的指令43，以及用于实施步骤e04的指令44。

指令41到44和处理器一起组合形成分别适于实施步骤e01到e04的系统1的各个模块。

图3是希望通过实施本发明的方法来表征的一个区域的地表下面的剖面图。

为此目的，地震仪100被埋设在区域的地表附近，并且地震仪100形成可以在所述剖面的平面中看到的一个群组101的一部分。作为举例，地震仪100被埋设在大约五十厘米的深度。这样的安装对于技术人员是特别简单的。

或者，地震仪可以被放置在地表上，前提是配置允许获得与地层的良好耦合。本领域技术人员知道如何放置地震仪以便获得良好的耦合。

在该例中，区域的地表下面包括区块200和区块300，其中区块200包含气体并且区块300包含水。该区域可能是储层。这两个具有不同相的流体的存在使得有可能实施本发明的方法。

在所示出的实例中，在群组101中，被安排在图的中部的地震仪100所给出的频谱将不同于被安排在右侧和左侧的地震仪100的频谱，这是因为仅有中部的地震仪处在储层的垂直上方。

为了允许用很少的装置在一个区域内进行记录，有可能分组进行测量。

举例来说，在第一天从午夜到早晨4点，地震仪100被安排形成群组101并且采集数据。在第二天从午夜到早晨4点，地震仪100被安排形成群组102并且采集数据。在第三天从午夜到早晨4点，地震仪100被安排形成群组103并且采集数据。在第四天从午夜到早晨4点，地震仪100被安排形成群组104并且采集数据。

图4是示出如何从信号获得频谱的图示，所述信号例如是通过参照图3描述的地震仪100所获得的信号。

在该图中，可以看到表示移动的信号sig，其特别表示沿着一个方向的移动速度。所述信号是在从午夜开始的四小时长的采集过程中所采集的：这是用来减少人为噪声的出现。

信号sig可以被划分成全部具有相同的持续时间的多个子信号，各个子信号是相继的，并且该例中的每一个子信号在前一个子信号的至少一半的持续时间上与前一个子信号重叠(这一重叠并不是至关重要的)。在图中，子信号由信号sig下方的弯括号表示。

可以观察到，如果其呈现出过多噪声，则随后可以省略其中一些子信号并且不作处理。作为指示，有可能消除具有超出一定阈值的绝对值的子信号。举例来说，有可能排除其绝对值位于针对整个信号sig定义的99％分位点之上的子信号。

随后，为每一个子信号准备一个子频谱。在图中示出了三个频谱密度子频谱：psd_1、psd_2和psd_3。

对于子频谱的每一个频率，根据对应于每一个子频谱中的该频率的频谱密度数值确定对应于频谱密度的中值。随后从所有各个中值获得频谱psd_m。换句话说，该频谱由这些中值构成。

图5示出了如何将属性组织在一个矩阵m内，该矩阵用于组织对于在不同的点处并且在可能不同的日期采集的信号所获得的属性。

在矩阵m中使用以下标记法：

attr_i：类型i的属性；

x_j：区域中的点j(该点与测量的日期相关联)；

f_k：频谱属性的频率k。

在矩阵m中，每一行与一条记录并且与区域的一点x_j相关联，每一列与属性类型attr_i并且与频谱属性频率f_k相关联。

在所述矩阵中，与给定的属性相关联的矩阵的各列全都是邻近的。因此，对应于各组同时进行的记录(例如在参照图3描述的群组101到104内)的所述矩阵的各行在矩阵中都被分组在一起从而形成各个行组，并且每一组在该例中与一天相关联。

优选的是，对于每一天或者对于每一个行组，属性的数值被归一化。

通过把属性组织在矩阵m中使得有可能实施主成分分析方法，其中每一行是一个个体并且每一列是一个变量。主成分分析用来获得被称作“投影器(projector)”的主成分。因此投影器是被写成p的矢量，其长度l等于不同属性类型的数目乘以对于每一种属性类型所存在的频率的数目。

矩阵m的索引为i(其处在1到m的范围内，其中m是矩阵m中的行的数目)的行到投影器p上的投影被如下计算：

由于该投影的结果对应于区域的一点，因此有可能获得在该处实施了信号采集的区域的所有各点处的矩阵m的投影的图形表示。

这样的图形表示被示出在图6中。该图示出了对应于投影的四个图形表示：prj1、prj2、prj3和prj4。已知的储层res的轮廓在这些图形表示上被突出显示。

呈现出最佳地对应于预期的空间表示的灰度变化的图形表示被认为是良好的投影器。本领域技术人员知道如何解释这些映射。

在区域的每一个图形表示下方，还示出了作为频率的函数的投影器本身。

在所示出的实例中，投影器prj1和prj2被认为是良好的投影器。因此，投影器的数目k被确定为等于两个投影器，这些投影器被写成p1和p2。

随后，可以通过使用下面的公式把矩阵投影到两个投影器p1和p2上：

<m|(p1,p2)>(i)＝(<m|p1>,<m|p2>)

这样就获得了属于二维空间的一定数目的m(矩阵中的行的数目或者记录的数目)个矢量。这些矢量可以被组织形成第二矩阵的各行。

图7在其由轴x1、x2表示的初始参照系中示出了存在于该第二矩阵中的各个个体。每一个个体对应于图中的一个叉号。

因此，通过应用被组织在一维中以便获得n个类的分类方法，有可能对这些个体进行分类，这是通过确定近似各个个体的被写成ξ的曲线，并且随后通过确定由沿着曲线ξ的曲线横坐标上的圆形表示的类，其中曲线ξ对应于类数(在该例中从-1到1变化)。可以观察到，所述圆形的半径对应于其协方差半径。

该图还示出对应于已被保留的两个投影器的轴，所述轴被参照为e1和e2，并且不会允许异常被足够好地表示。

在该例中，类数表示区域的地表下面的异常的量值。

因此，所述圆形的中心被视为类头。

于是有可能从各个类头获得维度是n乘k的第三矩阵。随后有可能实施矩阵的伪求逆以便获得第四矩阵，所述第四矩阵的维度是n乘以出现在频谱属性的初始矩阵的每一行中的频率的数目。

图8示出了通过对应于区域的每一点的类数显示异常的量值的数值而获得的映射，其中有可能在伪求逆之后获得该映射。

在该图中，对于每一个类头还示出了曲线，所述曲线揭示出作为对应于类头的属性中的频率的函数的变化(这是对于两个属性进行的)，作为频率的函数的对应于垂直移动的频谱的导数，以及垂直和水平移动之间的比值。

应当观察到，这些曲线显示出非常不同的行为，包括在低于1hz的频率范围内也是如此。

发明人已观察到，当在储层的垂直上方进行测量时，在频谱中由文献“phenomenologyoftremor-likesignalsobservedoverhydrocarbonreservoirs(在油气储层上观测到的类震动信号的现象学)”(s.dangel等人，journalofvolcanologyandgeothermalresearch)规定的峰值之前实际上存在下降。通过使用较大的频率范围并且通过使用主成分分析，本发明允许比现有技术的解决方案中更清楚地显示出异常。

还可以观察到，在本发明的描述中示出的各种图表是从在已知的储层之上和周围进行的测量所获得的。