一种时移重力增量的去噪方法和装置制造方法及图纸

本发明专利技术提供了一种时移重力增量的去噪方法和装置，其中，该方法包括：通过对工区内所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量。在本发明专利技术实施例中，采用加权处理的方法有效解决了现有技术中进行去噪处理时没有考虑进行时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系的问题，去除了时移重力增量的观测误差，即观测噪声，使得最终得到的数据与实际的地质情况吻合度更高。

【技术实现步骤摘要】

本专利技术涉及石油勘探
，特别涉及一种时移重力增量的去噪方法和装置。
技术介绍
在油气生产过程中，地下的油气物质会重新分布，从而导致地下的重力场产生局部的变化。针对地下重力场变化所进行的测量称为时移微重力测量，对地球微重力场产生的局部变化进行重复性测量所得到的数据称为时移微重力数据。因此，进行时移微重力测量时需要对同一工区进行多次观测，对相同测点进行多次观测多次读数，从而工区内的测点均存在时移微重力增量。目前，对微重力勘探数据进行去噪处理，即，对勘探数据中的个别突跳点进行平滑滤波所采用的方法主要有：向上延拓的去噪处理方法以及低通滤波的去噪处理方法。然而上述两种方法并没有考虑进行时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系，因而并不适用于时移微重力增量的去噪处理情况。针对上述问题，目前尚未提出有效的解决方案。
技术实现思路
本专利技术实施例提供了一种时移重力增量的去噪方法和装置，以解决现有技术中采用常规的去噪处理方法对时移重力数据进行去噪处理时，无法考虑时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系的问题。本专利技术实施例提供了一种时移重力数据的去噪方法，包括：获取工区内各个测点的时移重力增量，其中，所述时移重力增量为测点在第一时间段内多次测量的重力平均值和在第二时间段内多次测量的重力平均值之间的差值；获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差；对所述各个
测点中每个测点采用以下操作得到所述各个测点中每个测点去噪后的时移重力增量：选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。在一个实施例中，在得到所述各个测点中每个测点去噪后的时移重力增量之后，还包括：将所述各个测点中每个测点去噪后的时移重力增量作为未去噪的时移重力增量，对所述各个测点中每个测点重复执行以下去噪操作直至达到预设的重复次数：选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。在一个实施例中，所述预设的重复次数不小于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表示所述预设去噪区间中的测点个数。本专利技术实施例还提供了一种时移重力增量的去噪装置，包括：平均值获取模块，用于获取工区内各个测点的时移重力增量，其中，所述时移重力增量为测点在第一时间段内多次测量的重力平均值和在第二时间段内多次测量的重力平均值之间的差值；方差获取模块，用于获取所述各个测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差；加权处理模块，用于对所述各个测点中每个测点采用以下操作得到所述各个测点中每个测点去噪后的时移重力增量：选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。在一个实施例中，所述时移重力增量的去噪装置还包括：重复模块，用于在得到所述各个测点中每个测点去噪后的时移重力增量之后，将所述各个测点中每个测点去噪后的时移重力增量作为未去噪的时移重力增量，对所述各个测点中每个测点重复执行以下去噪操作直至达到预设的重复次数：单次重复模块，用于选择当前测点作为中心点，将以所述中心点为中心，以预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和未去噪的时移重力增量进行加权处理，得到去噪后的所述当前测点的时移重力增量。在一个实施例中，所述预设的重复次数不小于10次。在上述实施例中，在对时移重力增量进行去噪的过程中，通过对所有测点在第一时间段内多次测量的重力方差、在第二时间段内多次测量的重力方差和时移重力增量进行加权处理，得到去噪后的当前测点的时移重力增量，即采用加权处理的方法有效解决了现有技术中进行去噪处理时，没有考虑进行时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系的问题，去除了时移重力增量的观测误差，即观测噪声，使得最终得到的数据与实际的地质情况吻合度更高。附图说明此处所说明的附图用来提供对本专利技术的进一步理解，构成本申请的一部分，并不构成对本专利技术的限定。在附图中：图1是本专利技术实施例的一种时移重力增量的去噪方法的流程图；图2是本专利技术实施例的一种时移重力增量的去噪装置的一种结构框图。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本专利技术做进一步详细说明。在此，本专利技术的示意性实施方式及其说明用于解释本专利技术，但并不作为对本专利技术的限定。针对采用常规的去噪处理方法对时移重力数据进行去噪处理时所出现的问题，专利技术人提出了一种时移重力增量的去噪方法，以解决采用向上延拓的去噪处理方法以及低通滤波的去噪处理方法进行去噪处理时，并没有考虑进行时移微重力多次观测得到的数据与单次观测多次读数所得到的数据之间的联系的问题。具体地，该一种时移重力增量的去噪方法可以如图1所示，包括以下步骤：步骤101：获取工区内各个测点的时移重力增量，其中，所述时移重力增量为测点在第一时间段内多次测量的重力平均值和在第二时间段内多次本文档来自技高网...

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

【技术特征摘要】
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×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表示以所述预设去噪区间作为边界所限定的范围内的所有测点在第一
\t时间段内多次测量的重力方差，B表示在第二时间段内多次测量的重力方差，D表示以所述预设去噪区间作为边界所限定的范围内的所有测点在第一时间段内多次测量的重力方差以及在第二时间段内多次测量的重力方差中的最大...

【专利技术属性】
技术研发人员：刘云祥，赵文举，郑胜贤，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京;11