[发明专利]一种时移重力增量的去噪方法和装置

说明书

技术领域

本发明涉及石油勘探技术领域，特别涉及一种时移重力增量的去噪方法和装置。

背景技术

在油气生产过程中，地下的油气物质会重新分布，从而导致地下的重力场产生局部的变化。针对地下重力场变化所进行的测量称为时移微重力测量，对地球微重力场产生的局部变化进行重复性测量所得到的数据称为时移微重力数据。因此，进行时移微重力测量时需要对同一工区进行多次观测，对相同测点进行多次观测多次读数，从而工区内的测点均存在时移微重力增量。

目前，对微重力勘探数据进行去噪处理，即，对勘探数据中的个别突跳点进行平滑滤波所采用的方法主要有：向上延拓的去噪处理方法以及低通滤波的去噪处理方法。

然而，上述两种方法并没有考虑进行时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系，因而并不适用于时移微重力增量的去噪处理情况。

针对上述问题，目前尚未提出有效的解决方案。

发明内容

本发明实施例提供了一种时移重力增量的去噪方法和装置，以解决现有技术中采用常规的去噪处理方法对时移重力数据进行去噪处理时，无法考虑时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系的问题。

本发明实施例提供了一种时移重力数据的去噪方法，包括：获取工区内各个测点的时移重力增量，其中，所述时移重力增量为测点在第一时间段内多次测量的重力平均值和在第二时间段内多次测量的重力平均值之间的差值；获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差；对所述各个测点中每个测点采用以下操作得到所述各个测点中每个测点去噪后的时移重力增量：选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。

在一个实施例中，在得到所述各个测点中每个测点去噪后的时移重力增量之后，还包括：将所述各个测点中每个测点去噪后的时移重力增量作为未去噪的时移重力增量，对所述各个测点中每个测点重复执行以下去噪操作直至达到预设的重复次数：选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。

在一个实施例中，所述预设的重复次数不小于10次。

在一个实施例中，按照以下公式计算得到去噪后的所述当前测点的时移重力增量：

当A<(D/2)和B<(D/2)时，M＝N；

当(D/2)<A<(3D/4)和(D/2)<B<(3D/4)时，M＝N×a+P×b；

当(D/2)<A<(3D/4)和B>(3D/4)、(D/2)<B<(3D/4)和A>(3D/4)满足二者之一时， M＝N×c+P×d；

当A>(3D/4)且B>(3D/4)时，M＝N×e+P×f；

其中，A表示以所述预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差，B表示在第二时间段内多次测量的重力方差，D表示以所述预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差中的最大值，M表示去噪后的所述当前测点的时移重力增量，N表示所述时移重力增量，P表示以所述预设去噪区间作为边界所限定的范围内的所有测点的时移重力增量的平均值，a、b、c、d、e、f 均表示在0％至100％之间的加权系数。

在一个实施例中，获取工区内各个测点的时移重力增量，包括：对所述时移动重力增量进行网格化处理；获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差，包括：对所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差进行网格化处理。

在一个实施例中，所述网格化处理方法包括：克里金插值法。

在一个实施例中，所述预设去噪区间包括：n×n的矩阵，其中，所述n表示所述预设去噪区间中的测点个数。