一种时移重力增量的去噪方法和装置的制造方法
【专利摘要】本发明提供了一种时移重力增量的去噪方法和装置，其中，该方法包括：通过对工区内所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量。在本发明实施例中，采用加权处理的方法有效解决了现有技术中进行去噪处理时没有考虑进行时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系的问题，去除了时移重力增量的观测误差，即观测噪声，使得最终得到的数据与实际的地质情况吻合度更高。
【专利说明】
_种时移重力増量的去噪方法和装置
技术领域
[0001]本发明涉及石油勘探技术领域，特别涉及一种时移重力增量的去噪方法和装置。
【背景技术】
[0002]在油气生产过程中，地下的油气物质会重新分布，从而导致地下的重力场产生局部的变化。针对地下重力场变化所进行的测量称为时移微重力测量，对地球微重力场产生的局部变化进行重复性测量所得到的数据称为时移微重力数据。因此，进行时移微重力测量时需要对同一工区进行多次观测，对相同测点进行多次观测多次读数，从而工区内的测点均存在时移微重力增量。
[0003]目前，对微重力勘探数据进行去噪处理，S卩，对勘探数据中的个别突跳点进行平滑滤波所采用的方法主要有:向上延拓的去噪处理方法以及低通滤波的去噪处理方法。
[0004]然而上述两种方法并没有考虑进行时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系，因而并不适用于时移微重力增量的去噪处理情况。
[0005]针对上述问题，目前尚未提出有效的解决方案。
【发明内容】
[0006]本发明实施例提供了一种时移重力增量的去噪方法和装置，以解决现有技术中采用常规的去噪处理方法对时移重力数据进行去噪处理时，无法考虑时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系的问题。
[0007]本发明实施例提供了一种时移重力数据的去噪方法，包括:获取工区内各个测点的时移重力增量，其中，所述时移重力增量为测点在第一时间段内多次测量的重力平均值和在第二时间段内多次测量的重力平均值之间的差值;获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差;对所述各个测点中每个测点采用以下操作得到所述各个测点中每个测点去噪后的时移重力增量:选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。
[0008]在一个实施例中，在得到所述各个测点中每个测点去噪后的时移重力增量之后，还包括:将所述各个测点中每个测点去噪后的时移重力增量作为未去噪的时移重力增量，对所述各个测点中每个测点重复执行以下去噪操作直至达到预设的重复次数:选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。
[0009]在一个实施例中，所述预设的重复次数不小于10次。
[0010]在一个实施例中，按照以下公式计算得到去噪后的所述当前测点的时移重力增量:
[0011]当A〈(D/2)和B〈(D/2)时，M = N;
[0012]当(D/2)〈A〈(3D/4)和(D/2)〈B〈(3D/4)时，M= NXa+PXb;
[0013]当(0/2)〈厶〈(30/4)和8>(30/4)、(0/2)〈8〈(30/4)和厶>(30/4)满足二者之一时，]?=N X c+P X d ；
[0014]当A> (3D/4)且 B〉(3D/4)时，M = NXe+PXf;
[0015]其中，A表示以所述预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差，B表示在第二时间段内多次测量的重力方差，D表示以所述预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差中的最大值，M表示去噪后的所述当前测点的时移重力增量，N表示所述时移重力增量，P表示以所述预设去噪区间作为边界所限定的范围内的所有测点的时移重力增量的平均值，a、b、c、d、e、f均表示在O %至100 %之间的加权系数。
[0016]在一个实施例中，获取工区内各个测点的时移重力增量，包括:对所述时移动重力增量进行网格化处理;获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差，包括:对所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差进行网格化处理。
[0017]在一个实施例中，所述网格化处理方法包括:克里金插值法。
[0018]在一个实施例中，所述预设去噪区间包括:nXn的矩阵，其中，所述η表示所述预设去噪区间中的测点个数。
[0019]本发明实施例还提供了一种时移重力增量的去噪装置，包括:平均值获取模块，用于获取工区内各个测点的时移重力增量，其中，所述时移重力增量为测点在第一时间段内多次测量的重力平均值和在第二时间段内多次测量的重力平均值之间的差值;方差获取模块，用于获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差;加权处理模块，用于对所述各个测点中每个测点采用以下操作得到所述各个测点中每个测点去噪后的时移重力增量:选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。
[0020]在一个实施例中，所述时移重力增量的去噪装置还包括:重复模块，用于在得到所述各个测点中每个测点去噪后的时移重力增量之后，将所述各个测点中每个测点去噪后的时移重力增量作为未去噪的时移重力增量，对所述各个测点中每个测点重复执行以下去噪操作直至达到预设的重复次数:单次重复模块，用于选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。
[0021 ]在一个实施例中，所述预设的重复次数不小于10次。
[0022]在上述实施例中，在对时移重力增量进行去噪的过程中，通过对所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量，即采用加权处理的方法有效解决了现有技术中进行去噪处理时，没有考虑进行时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系的问题，去除了时移重力增量的观测误差，即观测噪声，使得最终得到的数据与实际的地质情况吻合度更高。
【附图说明】
[0023]此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中:
[0024]图1是本发明实施例的一种时移重力增量的去噪方法的流程图；
[0025]图2是本发明实施例的一种时移重力增量的去噪装置的一种结构框图。
【具体实施方式】
[0026]为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
[0027]针对采用常规的去噪处理方法对时移重力数据进行去噪处理时所出现的问题，发明人提出了一种时移重力增量的去噪方法，以解决采用向上延拓的去噪处理方法以及低通滤波的去噪处理方法进行去噪处理时，并没有考虑进行时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系的问题。具体地，该一种时移重力增量的去噪方法可以如图1所示，包括以下步骤:
[0028]步骤101:获取工区内各个测点的时移重力增量，其中，所述时移重力增量为测点在第一时间段内多次测量的重力平均值和在第二时间段内多次测量的重力平均值之间的差值；
[0029]步骤102:获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差；
[0030]其中，本实施例中和重力相关的数据均特指微重力。微重力测量具体指的是，微弱重力异常的测量技术与方法。微重力也称为微伽级重力，毫伽是重力测量的常用单位，微伽是毫伽的千分之一。
[0031]在实际应用的过程中进行重力测量时，需要使用多台仪器对同一个测点进行多次测量。其中，重力平均值指的是多台仪器多次重力测量结果的平均值。同样地，重力方差指的是多台仪器多次重力测量结果的方差值。
[0032]其中，上述第一时间段和第二时间段为两个独立的时间段，二者之间没有重叠。
[0033]分别对上述观测得到的各个测点的时移重力增量、第一重力方差以及第二重力方差进行网格化处理。在本实施例中，可以采用克里金插值法分别对上述观测得到的各个测点的时移重力增量、第一重力方差以及第二重力方差进行网格化处理，从而可以得到时移重力增量网格文件、第一重力方差网格文件以及第二重力方差网格文件。
[0034]步骤103:对所述各个测点中每个测点采用以下操作得到所述各个测点中每个测点去噪后的时移重力增量:选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量。
[0035]进一步地，在得到各个测点中每个测点去噪后的时移重力增量之后，可以将各个测点中每个测点去噪后的时移重力增量作为未去噪的时移重力增量，对各个测点中每个测点重复执行以下去噪操作直至达到预设的重复次数:选择当前测点作为中心点，将以该中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量。
[0036]其中，上述预设去噪区间可以包括但不限于:nXn的矩阵，其中，所述η表示预设去噪区间中的测点个数;同时，上述预设的重复次数可以设定为不小于10次。
[0037]在本实施例中，可以按照以下公式计算得到去噪后的当前测点的时移重力增量:
[0038]当A〈(D/2)和B〈(D/2)时，M = N;
[0039 ]当(D/2) <A〈( 3D/4)和(D/2) <B〈( 3D/4)时，M = N X a+P X b;
[0040]当(0/2)〈厶〈(30/4)和8>(30/4)、(0/2)〈8〈(30/4)和厶>(30/4)满足二者之一时，]? =N X c+P X d ；
[0041 ]当 A> (3D/4)且 B〉(3D/4)时，M = NXe+PXf;
[0042]其中，A表示以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差，B表示在第二时间段内多次测量的重力方差，D表示以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差中的最大值，M表示去噪后的当前测点的时移重力增量，N表示时移重力增量，P表示以预设去噪区间作为边界所限定的范围内的所有测点的时移重力增量的平均值，a、b、c、d、e、f均表示在O %至100 %之间的加权系数。
[0043]例如:a可以设置为10%，此时，相应的，b可以设置为90% ;c可以设置为21%，此时，相应的，d可以设置为79 % ; e可以设置为41 %，此时，相应的，f可以设置为59 %。值得注意的是，上述对a、b、c、d、e、f的取值只是为了更好的说明本申请，也可以取别的数值，本申请对此不作限定。但是，8、13、0、(1、6、;1^需要满足以下要求:3+&=100% ,c+d= 100%，e+f =100%。
[0044]下面以对某一测点进行时移重力增量去噪处理为例说明本实施例的去噪过程，可以包括:首先，确定进行去噪处理的当前测点并确定去噪处理的重复次数;再以该测点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量;然后，将得到的去噪后的当前测点的时移重力增量作为未去噪的时移重力增量，重复上述去噪处理的过程，直至达到上述确定的去噪处理的重复次数;最后，得到去噪后的当前测点的时移重力增量。对所测定的时移重力增量网格文件、第一重力方差网格文件以及第二重力方差网格文件中的所有测点均执行上述去噪过程，可以得到去噪后的时移重力增量网格文件。
[0045]下面以某地区的时移重力增量网格文件、第一重力方差网格文件以及第二重力方差网格文件进行去噪处理为例对本实施例所提供的一种时移重力增量的去噪方法进行说明，然而值得注意的是，该具体实施例仅是为了更好地说明本发明，并不构成对本发明的不当限定。
[0046]在本实施例中，可以选择5X5的矩阵，S卩n = 5，作为进行时移重力增量去噪运算的区间。当然，也可以选择3X3的矩阵作为进行时移重力增量去噪运算的区间，本发明对此不作限定。
[0047]进一步地，设置时移重力增量去噪处理的重复次数为50次。在本实施例中，选择50次作为去噪处理的重复次数是因为进行50去噪后所得到的时移重力增量网格文件的去噪效果最佳。同样地，进行时移重力去噪处理时也可以进行60次去噪处理，本发明对此不作限定。
[0048]相应地，进行去噪处理时，对选定的以待去噪测点区间，可以按照以下规则进行加权处理:
[0049]当A〈4微伽并且B〈4微伽时，M= N;
[0050]当4微伽<A〈6微伽并且4微伽<B〈6微伽时，M = NX80%+PX20% ；
[0051 ]当4微伽<A〈6微伽并且B>6微伽、或者4微伽<B〈6微伽并且A>6微伽，这两个条件满足二者之一时，M=NX60%+PX40% ；
[0052]当A>6微伽并且B>6微伽时，M= NX30%+PX70% ；
[0053]其中，A表示以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差，B表示在第二时间段内多次测量的重力方差，D表示以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差中的最大值，M表示去噪后的当前测点的时移重力增量，N表示时移重力增量，P表示以预设去噪区间作为边界所限定的范围内的所有测点的时移重力增量的平均值。
[0054]对某地区的时移重力增量网格文件、第一重力方差网格文件以及第二重力方差网格文件进行去噪处理的过程中。首先，可以确定进行去噪处理的当前测点;再以该测点为中心，以5X5的矩阵区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量;然后，将得到的去噪后的当前测点的时移重力增量作为未去噪的时移重力增量，重复上述去噪处理的过程，直至重复次数达到50次后，可以得到去噪后的当前测点的时移重力增量。对某地区所测定的时移重力增量网格文件、第一重力方差网格文件以及第二重力方差网格文件中的所有测点均执行上述去噪过程，可以得到去噪后的时移重力增量网格文件。
[0055]基于同一发明构思，本发明实施例中还提供了一种时移重力增量的去噪装置，如下面的实施例所述。由于时移重力增量的去噪装置解决问题的原理与时移重力增量的去噪方法相似，因此时移重力增量的去噪装置的实施可以参见时移重力增量的去噪方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。图2是本发明实施例的时移重力增量的去噪装置的一种结构框图，如图2所示，包括:平均值获取模块201、方差获取模块202和加权处理模块203，下面对该结构进行说明。
[0056]平均值获取模块201，用于获取工区内各个测点的时移重力增量，其中，所述时移重力增量为测点在第一时间段内多次测量的重力平均值和在第二时间段内多次测量的重力平均值之间的差值；
[0057]方差获取模块202，用于获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差；
[0058]加权处理模块203，用于对所述各个测点中每个测点采用以下操作得到所述各个测点中每个测点去噪后的时移重力增量:选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。
[0059]在一个实施例中，所述装置还可以包括:重复模块，用于在得到所述各个测点中每个测点去噪后的时移重力增量之后，将所述各个测点中每个测点去噪后的时移重力增量作为未去噪的时移重力增量，对所述各个测点中每个测点重复执行以下去噪操作直至达到预设的重复次数:选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。
[0060]在一个实施例中，所述预设的重复次数不小于10次。
[0061]在一个实施例中，所述加权处理模块可以按照以下公式计算得到去噪后的所述当前测点的时移重力增量:
[0062]当A〈(D/2)和B〈(D/2)时，M = N;
[0063 ]当(D/2) <A〈( 3D/4)和(D/2) <B〈( 3D/4)时，M = N X a+P X b;
[0064]当(0/2)〈厶〈(30/4)和8>(30/4)、(0/2)〈8〈(30/4)和厶>(30/4)满足二者之一时，]?=N X c+P X d ；
[0065]当A>(3D/4)且 B>(3D/4)时，M = NXe+PXf;
[0066]其中，A表示以所述预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差，B表示在第二时间段内多次测量的重力方差，D表示以所述预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差中的最大值，M表示去噪后的所述当前测点的时移重力增量，N表示所述时移重力增量，P表示以所述预设去噪区间作为边界所限定的范围内的所有测点的时移重力增量的平均值，a、b、c、d、e、f均表示在O %至100 %之间的加权系数。
[0067]在一个实施例中，所述平均值获取模块用于对所述时移动重力增量进行网格化处理;所述方差获取模块用于对所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差进行网格化处理。
[0068]在一个实施例中，所述网格化处理方法包括:克里金插值法。
[0069]在一个实施例中，所述预设去噪区间包括:nXn的矩阵，其中，所述η表示所述预设去噪区间中的测点个数。
[0070]从以上的描述中，可以看出，本发明实施例实现了如下技术效果:在对时移重力增量进行去噪的过程中，通过对所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量，即采用加权处理的方法有效解决了现有技术中进行去噪处理时没有考虑进行时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系的问题，去除了时移重力增量的观测误差，即观测噪声，使得最终得到的数据与实际的地质情况吻合度更尚O
[0071]显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
[0072]以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1.一种时移重力增量的去噪方法，其特征在于，包括: 获取工区内各个测点的时移重力增量，其中，所述时移重力增量为测点在第一时间段内多次测量的重力平均值和在第二时间段内多次测量的重力平均值之间的差值； 获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差； 对所述各个测点中每个测点采用以下操作得到所述各个测点中每个测点去噪后的时移重力增量:选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量。2.如权利要求1所述的方法，其特征在于，在得到所述各个测点中每个测点去噪后的时移重力增量之后，所述方法还包括: 将所述各个测点中每个测点去噪后的时移重力增量作为未去噪的时移重力增量，对所述各个测点中每个测点重复执行以下去噪操作直至达到预设的重复次数:选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量。3.如权利要求2所述的方法，其特征在于，所述预设的重复次数不小于10次。4.如权利要求1所述的方法，其特征在于，按照以下公式计算得到去噪后的所述当前测点的时移重力增量: 当 A〈(D/2)和 B〈(D/2)时，M=N; 当(D/2) <A< (3D/4)和(D/2) <B< (3D/4)时，M=N X a+P X b; 当(0/2)〈厶〈(30/4)和8>(30/4)、(0/2)〈8〈(30/4)和厶>(30/4)满足二者之一时，]?=~\0+PXd； 当八>(30/4)且8>(30/4)时，]\1=~\6+?\匕 其中，A表示以所述预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差，B表示在第二时间段内多次测量的重力方差，D表示以所述预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差中的最大值，M表示去噪后的所述当前测点的时移重力增量，N表示所述时移重力增量，P表示以所述预设去噪区间作为边界所限定的范围内的所有测点的时移重力增量的平均值，a、b、c、d、e、f均表示在O %至100 %之间的加权系数。5.如权利要求1所述的方法，其特征在于，获取工区内各个测点的时移重力增量，包括:对所述时移动重力增量进行网格化处理； 获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差，包括:对所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差进行网格化处理。6.如权利要求5所述的方法，其特征在于，所述网格化处理方法包括:克里金插值法。7.如权利要求1所述的方法，其特征在于，所述预设去噪区间包括:nXn的矩阵，其中，所述η表示所述预设去噪区间中的测点个数。8.一种时移重力增量的去噪装置，其特征在于，包括: 平均值获取模块，用于获取工区内各个测点的时移重力增量，其中，所述时移重力增量为测点在第一时间段内多次测量的重力平均值和在第二时间段内多次测量的重力平均值之间的差值； 方差获取模块，用于获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差； 加权处理模块，用于对所述各个测点中每个测点采用以下操作得到所述各个测点中每个测点去噪后的时移重力增量:选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。9.如权利要求8所述的装置，其特征在于，所述装置还包括: 重复模块，用于在得到所述各个测点中每个测点去噪后的时移重力增量之后，将所述各个测点中每个测点去噪后的时移重力增量作为未去噪的时移重力增量，对所述各个测点中每个测点重复执行以下去噪操作直至达到预设的重复次数:选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。10.如权利要求9所述的装置，其特征在于，所述预设的重复次数不小于10次。
【文档编号】G01V7/06GK105974488SQ201610517205
【公开日】2016年9月28日
【申请日】2016年7月4日
【发明人】刘云祥, 赵文举, 郑胜贤
【申请人】中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司