CN113126172B - 静位移校正方法及装置 - Google Patents
静位移校正方法及装置 Download PDFInfo
- Publication number
- CN113126172B CN113126172B CN202010045597.7A CN202010045597A CN113126172B CN 113126172 B CN113126172 B CN 113126172B CN 202010045597 A CN202010045597 A CN 202010045597A CN 113126172 B CN113126172 B CN 113126172B
- Authority
- CN
- China
- Prior art keywords
- apparent resistivity
- original
- resistivity
- apparent
- preset
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012937 correction Methods 0.000 title claims abstract description 185
- 230000003068 static effect Effects 0.000 title claims abstract description 176
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 165
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000001914 filtration Methods 0.000 claims abstract description 134
- 230000000007 visual effect Effects 0.000 claims abstract description 62
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 abstract description 37
- 238000012545 processing Methods 0.000 abstract description 18
- 238000010586 diagram Methods 0.000 description 20
- 238000005259 measurement Methods 0.000 description 13
- 238000012935 Averaging Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000005684 electric field Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 235000012149 noodles Nutrition 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开了一种静位移校正方法及装置,该方法包括:对原始视电阻率进行滤波得到第二视电阻率;利用相对误差不小于预设相对误差阈值的第二视电阻率中的视电阻率替换原始视电阻率的视电阻率得到第三视电阻率数据体;确定第三视电阻率数据体中第四视电阻率;对第四视电阻率进行滤波得到滤波结果第六视电阻率;根据第四视电阻率与对应第六视电阻率,第二视电阻率数据体,确定静校正后的第四视电阻率数据体。本发明能够消除视电阻率中的畸变点,提高静位移校正效果;同时基于相对误差消除畸变点,无需其它处理即可消除视电阻率中的畸变点,提高静位移校正效率。
Description
技术领域
本发明涉及地球物理勘探技术领域,尤其涉及静位移校正方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
在大地电磁探测勘探中,当近地表存在局部电性不均匀体时,在外电场的作用下电性不均匀体表面会产生积累电荷,由此产生一个附加电场,导致观测电场产生畸变,该畸变几乎与频率无关。造成实测视电阻率发生畸变,主要表现为在双对数坐标系中单点曲线整体沿电阻率轴向上或者向下平移了一个固定值,而在剖面等值线中具体表现为横向间隔等值线密集分布的挂面条现象,而相位数据不受影响,这种现象通常称为静态效应或者静位移。
静位移现象在大地电磁勘探中普遍存在,并将会对数据处理、反演及解释结果的精度或者准确性产生严重影响。为了消除这种影响必须进行静位移校正处理。目前国内外静位移校正主要是有曲线平面法、空间滤波法及阻抗张量分解法、瞬变电磁校正发及相位校正法等等。然而,这些方法主要应用于二维大地电磁勘探中,在三维大地电磁勘探中较少应用。而应用于三维大地电磁的静位移校正方法,不仅存在校正效果差的问题,同时还存在校正效率低的问题。
因此,现有的静位移校正存在校正效果差及校正效率低的问题。
发明内容
本发明实施例提供一种静位移校正方法,用以提高静位移校正效果及静位移校正效率,该方法包括:
确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率;原始视电阻率为原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率;
对每个原始视电阻率进行滤波,得到每个原始视电阻率在原始位置的每个第二视电阻率;
根据每个原始视电阻率与每个第二视电阻率中每个对应位置的视电阻率的相对误差,利用相对误差不小于预设相对误差阈值的每个第二视电阻率中的视电阻率替换每个原始视电阻率中对应位置的视电阻率,得到每个第三视电阻率;视电阻率的相对误差反映视电阻率的畸变程度;
确定第三视电阻率数据体中每个第四视电阻率;第三视电阻率数据体为每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值形成的第二视电阻率数据体,与原始视电阻率数据体中每个对应测点的视电阻率相加得到的;每个第四视电阻率为第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率的平均视电阻率;
对每个第四视电阻率进行滤波,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率;
根据每个第四视电阻率与对应每个第六视电阻率,及第二视电阻率数据体,确定对原始视电阻率数据体静校正后的第四视电阻率数据体。
本发明实施例还提供一种静位移校正装置,用以提高静位移校正效果及静位移校正效率,该装置包括:
原始视电阻率确定模块,用于确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率;原始视电阻率为原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率;
第一滤波模块,用于对每个原始视电阻率进行滤波,得到每个原始视电阻率在原始位置的每个第二视电阻率;
替换模块,用于根据每个原始视电阻率与每个第二视电阻率中每个对应位置的视电阻率的相对误差,利用相对误差不小于预设相对误差阈值的每个第二视电阻率中的视电阻率替换每个原始视电阻率中对应位置的视电阻率,得到每个第三视电阻率;视电阻率的相对误差反映视电阻率的畸变程度;
第四视电阻率确定模块,用于确定第三视电阻率数据体中每个第四视电阻率;第三视电阻率数据体为每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值形成的第二视电阻率数据体,与原始视电阻率数据体中每个对应测点的视电阻率相加得到的;每个第四视电阻率为第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率的平均视电阻率;
第二滤波模块,用于对每个第四视电阻率进行滤波,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率;
校正结果获得模块,用于根据每个第四视电阻率与对应每个第六视电阻率,及第二视电阻率数据体,确定对原始视电阻率数据体静校正后的第四视电阻率数据体。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述静位移校正方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述静位移校正方法的计算机程序。
本发明实施例中,利用相对误差不小于预设相对误差阈值的第二视电阻率中的视电阻率替换原始视电阻率中对应位置的视电阻率,能够基于相对误差消除视电阻率的畸变点,从而提高对视电阻率的静位移校正效果。另外,本发明实施例只需根据视电阻率的相对误差,即可消除视电阻率中的畸变点,无需其它额外的处理即可快速消除畸变点,能够提高静位移校正的效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例提供的静位移校正方法的实现流程图;
图2为本发明实施例提供的静位移校正方法中步骤101的实现流程图;
图3为本发明实施例提供的静位移校正方法中步骤102的实现流程图;
图4为本发明实施例提供的静位移校正方法中步骤103的实现流程图;
图5为本发明实施例提供的静位移校正方法中步骤104的实现流程图;
图6为本发明实施例提供的静位移校正方法中步骤105的实现流程图;
图7为本发明实施例提供的静位移校正方法中步骤106的实现流程图;
图8为本发明实施例提供的静位移校正方法的另一实现流程图;
图9为本发明实施例提供的静位移校正装置的功能模块图;
图10为本发明实施例提供的静位移校正装置中原始视电阻率确定模块901的结构框图;
图11为本发明实施例提供的静位移校正装置中第一滤波模块902的结构框图;
图12为本发明实施例提供的静位移校正装置中替换模块903的结构框图;
图13为本发明实施例提供的静位移校正装置中第四视电阻率确定模块904的结构框图;
图14为本发明实施例提供的静位移校正装置中第二滤波模块905的结构框图;
图15为本发明实施例提供的静位移校正装置中校正结果获得模块906的结构框图;
图16为本发明实施例提供的静位移校正装置的另一模块结构图;
图17为本发明实施例提供的静位移校正前XY方向上1.035Hz频率的视电阻率平面等值线示意图;
图18为本发明实施例提供的静位移校正后XY方向上1.035Hz频率的视电阻率平面等值线示意图;
图19为本发明实施例提供的静位移校正前XY方向上某号测线的视电阻率剖面等值线示意图;
图20为本发明实施例提供的静位移校正后XY方向上某号测线的视电阻率剖面等值线示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
图1示出了本发明实施例提供的静位移校正方法的实现流程,为便于描述,仅示出了与本发明实施例相关的部分,详述如下:
如图1所示,静位移校正方法,其包括:
步骤101,确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率;原始视电阻率为原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率;
步骤102,对每个原始视电阻率进行滤波,得到每个原始视电阻率在原始位置的每个第二视电阻率;
步骤103,根据每个原始视电阻率与每个第二视电阻率中每个对应位置的视电阻率的相对误差,利用相对误差不小于预设相对误差阈值的每个第二视电阻率中的视电阻率替换每个原始视电阻率中对应位置的视电阻率,得到每个第三视电阻率;视电阻率的相对误差反映视电阻率的畸变程度;
步骤104,确定第三视电阻率数据体中每个第四视电阻率;第三视电阻率数据体为每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值形成的第二视电阻率数据体,与原始视电阻率数据体中每个对应测点的视电阻率相加得到的;每个第四视电阻率为第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率的平均视电阻率;
步骤105,对每个第四视电阻率进行滤波,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率;
步骤106,根据每个第四视电阻率与对应每个第六视电阻率,及第二视电阻率数据体,确定对原始视电阻率数据体静校正后的第四视电阻率数据体。
三维大地电磁静位移校正处理的原始视电阻率数据体,主要包括XY方向上的视电阻率和YX方向上的视电阻率。这两个方向数据的静位移校正的处理方法是一样的。其中,预设方向为预先设定的方向,例如可以仅是XY方向,也可是仅是YX方向,本领域技术人员可以理解的是,预设方向还可以是XY方向和YX方向,本发明实施例对此不作特别的限制。原始视电阻率数据体,是由预设方向上的所有视电阻率数据组成的。其中,原始视电阻率数据体为三维数据体。
其中,原始视电阻率数据体包括多个测点,每个测点包括多个频率对应的视电阻率。为描述方便,此处我们把每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率称为原始视电阻率。可以理解的是,原始视电阻率数据体对应于多个原始视电阻率。每个测点对应的原始视电阻率均为该测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均值。可以理解的是,每个原始视电阻率均为平面视电阻率数据。
其中,该预设频率为预先设定的频率,例如预先设定该预设频率为1Hz,本领域技术人员可以理解的是,还可以预先设定该预设频率为除上述1Hz之外的其它频率,例如预先设定该预设频率为0.5Hz,1.5Hz或者2Hz等等,本发明实施例对此不作特别的限制。在较优的一实施例中,该预设频率为1Hz。
其中,与预设频率邻近的预设多个中频段频率,为预先设定的与预设频率邻近的多个中频段中的频率。例如,可以预先设定该预设多个中频段频率分别为0.9Hz、1Hz及1.1Hz,本领域技术人员可以理解的是,还可以预先设定该预设多个中频段频率分别为除上述0.9Hz、1Hz及1.1Hz之外的其它中频段的频率,例如预先设定该预设多个中频段频率分别为0.8Hz、1Hz及1.2Hz,本发明实施例对此不作特别的限制。
假设该预设频率为1Hz,某个测点中与1Hz邻近的多个中频段频率分别为频率A、频率B及频率C,则原始视电阻率为频率A的视电阻率、频率B的视电阻率及频率C的视电阻率三者的平均视电阻率。至此,原始视电阻率数据体中每个不同的测点对应有每个不同的原始视电阻率。
在确定每个测点对应的原始视电阻率后,对每个测点对应的每个原始视电阻率进行滤波,以得到每个原始视电阻率在原始位置的滤波结果,即第二视电阻率。可以理解的是,每个测点对应每个原始视电阻率,每个原始视电阻率对应每个第二视电阻率。其中,第二视电阻率即为原始视电阻率在原始位置的滤波结果,即第二视电阻率中的数据点的位置与原始视电阻率中的数据点的位置是一一对应的。此处对原始视电阻率进行滤波所采用的滤波方法,可以为现有技术中常用的滤波方法,此处不再详细赘述。可以理解的是,每个第二视电阻率均为平面视电阻率数据。
在得到每个原始视电阻率的滤波结果每个第二视电阻率后,可以利用第二视电阻率对原始视电阻率中的畸变点进行消除。具体的,首先确定原始视电阻率与第二视电阻率中每个对应位置的数据点的相对误差,该视电阻率的相对误差反映了数据点(视电阻率)的畸变程度,若原始视电阻率中数据点的视电阻率与第二视电阻率中对应数据点的视电阻率的相对误差不小于预设相对误差阈值,说明原始视电阻率中数据点的视电阻率的畸变程度较大,则将原始视电阻率中数据点认为是畸变数据点。
因此,为了消除原始视电阻率中产生畸变的数据点,可以利用第二视电阻率中与原始视电阻率中该畸变数据点对应位置的数据点,替换原始视电阻率中的该畸变数据点,从而达到消除原始视电阻率中畸变数据点的目的。而对于相对误差小于预设相对误差阈值的原始视电阻率中的数据点,说明该部分数据点产生的畸变较小,即认为相对误差小于预设相对误差阈值的原始视电阻率中的数据点没有发生畸变,则保留该部分数据点,据此形成消除畸变后的第三视电阻率。
可以理解的是,鉴于第二视电阻率中数据点的位置与原始视电阻率中数据点的位置是对应的,而第三视电阻率是利用第二视电阻率中对应位置的数据点,替换原始视电阻率对应位置的畸变数据点而来的,由此可以得知,第三视电阻率中数据点的位置与原始视电阻率数据点的位置也是对应的。可以理解的是,每个第三视电阻率均为平面视电阻率数据。
其中,该预设相对误差阈值为预设相对误差阈值区间中的任一相对误差。其中,该预设相对误差阈值区间为预先设定的相对误差阈值区间,例如可以设定该预设相对误差阈值区间为0.01至0.5,本领域技术人员可以理解的是,还可以预先设定该预设相对误差阈值区间为除上述0.01至0.5之外的其它相对误差阈值区间,例如0.02至0.48,或者0.05至0.45等等,本领域技术人员对此不作特别的限制。在较优的一实施例中,该预设相对误差区间为0.01至0.5。
其中,预设相对误差阈值为预先设定的相对误差阈值,且为预设相对误差阈值区间中的任一相对误差。例如,在预设相对误差阈值区间为0.01至0.5时,可以预先设定该预设相对误差阈值为0.08,本领域技术人员可以理解的是,还可以预先设定该预设相对误差阈值为除上述0.08之外的其它相对误差值,例如0.06,或者0.10等等,本领域技术人员对此不作特别的限制。在较优的一实施例中,该预设相对误差阈值为0.08。
至此,针对每个测点对应的每个原始视电阻率及每个第二视电阻率,均可以利用上述相对误差,确定对应的每个第三视电阻率。可以理解的是,该第三视电阻率为平面视电阻率数据。
在得到每个测点对应的每个第三视电阻率后,确定每个原始视电阻率与每个第三视电阻率中每个对应测点之间的差值,所有的原始视电阻率与对应的第三视电阻率中每个对应测点之间的差值,形成第二视电阻率数据体。可以理解的是,该第二视电阻率数据体为三维视电阻率数据体。在得到该第二视电阻率数据体后,进而将该第二视电阻率数据体与原始视电阻率数据体中每个对应测点之间的视电阻率相加,即可得到第三视电阻率数据体。可以理解的是,该第三视电阻率数据体为三维视电阻率数据。
在本发明实施例中,第三视电阻率数据体包括多个测点,每个测点对应每个第四视电阻率。该每个第四视电阻率为第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率的平均视电阻率。
其中,该多个预设高频段频率号为预先设定的多个高频段的频率号,例如预先设定该预设高频段频率号分别为高频段的频率号为2、3及4的频率号,本领域技术人员可以理解的是,还可以预先设定该预设高频段频率号为除上述频率号2、频率号3及频率号4的频率号之外的其它频率号,本领域技术人员对此不作特别的限制。在较优的一实施例中,该多个预设高频段频率号为高频段中频率号为2、3及4的频率号。
据此,假设该多个预设高频段频率号为高频段中频率号为A、B及C的频率号,即频率号A、频率号B及频率号C。则每个第四视电阻率为每个测点中频率号A、频率号B及频率号C的视电阻率的平均视电阻率,以此确定第三视电阻率数据体中每个测点对应的每个第四视电阻率。
在得到第三视电阻率数据体中每个测点对应的每个第四视电阻率后,分别对每个第四视电阻率进行滤波,得到每个第四视电阻率在原始位置的滤波结果,即对应的每个第六视电阻率。可以理解的是,每个第四视电阻率及对应的每个第六视电阻率均为平面视电阻率数据。此处对第四视电阻率进行滤波所采用的滤波方法,可以为现有技术中常用的滤波方法,此处不再详细赘述。该每个第六视电阻率中的数据点的位置与对应每个第四视电阻率中数据点的位置是对应的。
在得到每个第四视电阻率对应的每个滤波结果,即每个第六视电阻率后,进而根据该每个第四视电阻率,对应的每个第六视电阻率,及第二视电阻率数据体,确定对原始视电阻率数据体静校正后的第四视电阻率数据体。可以理解的是,该第四视电阻率数据体为三维视电阻率数据,该第四视电阻率数据体中的视电阻率数据的位置与原始视电阻率数据体中的视电阻率数据的位置是对应的。
在本发明实施例中,利用相对误差不小于预设相对误差阈值的第二视电阻率中的视电阻率替换原始视电阻率中对应位置的视电阻率,能够基于相对误差消除视电阻率的畸变点,从而提高对视电阻率的静位移校正效果。另外,本发明实施例只需根据视电阻率的相对误差,即可确定视电阻率中的畸变点,无需其它额外的处理即可快速消除畸变点,能够提高静位移校正的效率。
图2示出了本发明实施例提供的静位移校正方法中步骤101的实现流程,为便于描述,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效率,如图2所示,步骤101,确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率,包括:
步骤201,获取原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率;
步骤202,将每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率作为每个原始视电阻率。
为了提高静位移校正的效率,只需确定每个测点中与预设频率邻近的预设多个中频段频率的视电阻率,进而计算该几个中频段频率的视电阻率的平均值,即为每个测点对应的的原始视电阻率。
具体的,本发明实施例中该预设频率为1Hz。首先获取每个测点中与1Hz邻近的预设多个中频段频率。假设获取的与1Hz邻近的预设多个中频段频率分别为频率A、频率B及频率C。进而计算频率A的视电阻率,频率B的视电阻率,及频率C的视电阻率三者的平均视电阻率,该三者的平均视电阻率即为原始视电阻率。据此可以得到每个测点对应的每个原始视电阻率。
在本发明实施例中,获取原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率,在确定原始视电阻率时只需要计算每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率,即可快速确定原始视电阻率,因而能够进一步提高静位移校正的效率。
图3示出了本发明实施例提供的静位移校正方法中步骤102的实现流程,为便于描述,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效果,如图3所示,步骤102,对每个原始视电阻率进行滤波,得到每个原始视电阻率在原始位置的每个第二视电阻率,包括:
步骤301,采用规则插值网格对每个原始视电阻率进行克里金插值,得到每个第一插值视电阻率;规则插值网格的网格间隔与原始视电阻率的数据间隔的比值为预设倍数区间内的任一倍数;
步骤302,采用第一滤波窗口对每个第一插值视电阻率进行滑动平均滤波,得到每个第一视电阻率;第一滤波窗口的中心位置数据不参与滤波;
步骤303,对每个第一视电阻率进行反插值,得到每个原始视电阻率在原始位置的滤波结果每个第二视电阻率。
在对原始视电阻率进行滤波处理时,可以依次对原始视电阻率进行规则插值网格的克里金插值、滑动平均滤波及反插值处理。
具体的,首先利用规则插值网格对每个原始视电阻率进行克里金插值,从而确定插值后的每个第一插值视电阻率。其中,在本发明的一实施例中,克里金插值所采用的规则插值网格的网格间隔不大于每个原始视电阻率的数据点间隔,或者说该规则插值网格的网格间隔与每个原始视电阻率的数据点间隔的比值为预设比值。鉴于规则插值网格的网格间隔不大于原始视电阻率的数据点间隔,因此该预设比值为不大于1的值。
例如,该规则插值网格的网格间隔与每个原始视电阻率的数据点间隔相同,或者该规则插值网格的网格间隔与每个原始视电阻率的数据点间隔的比值为0.95或者0.9,即规则插值网格的网格间隔略小于每个原始视电阻率的数据点间隔。本领域技术人员可以理解的是,该规则插值网格的网格间隔与每个原始视电阻率的数据点间隔的比值,还可以除上述0.95或者0.9之外的其它不大于1的值,例如1或者0.9,又或者0.98等等,本发明实施例对此不作特别的限制。
另一方面,该规则插值网格的节点密度为每个原始视电阻率的节点密度的预设倍数。该预设倍数为预设倍数区间中的任一倍数。其中,该预设倍数区间为预先设定的倍数区间,例如可以设定该预设倍数区间为1至2倍,本领域技术人员可以理解的是,还可以预先设定该预设倍数区间为除上述1至2倍之外的其它倍数区间,例如1.2至1.8倍,或者1.5至1.75倍等等,本领域技术人员对此不作特别的限制。在较优的一实施例中,该设定该预设倍数区间为2至3倍。
其中,预设倍数为预先设定的倍数,且为预设倍数区间中的任一倍数。例如,在预设倍数区间为1至2倍时,可以预先设定该预设倍数为1.5倍,本领域技术人员可以理解的是,还可以预先设定该预设倍数为除上述1.5倍之外的其它相对误差值,例如1.4,或者1.6等等,本领域技术人员对此不作特别的限制。
在采用规则插值网格对每个原始视电阻率进行克里金插值,得到对应的每个第一插值视电阻率后,紧接着对得到的对应每个第一插值视电阻率进行滑动平均滤波。
具体的,采用第一滤波窗口对每个第一插值视电阻率进行滑动平均滤波,得到滑动平均滤波后的第一视电阻率。在本发明的一实施例中,该滑动平均滤波的滤波系数为第一预设滤波系数区间中任一滤波系数。另外,第一滤波窗口的中心位置数据不参与滤波。
其中,该第一预设滤波系数为第一预设滤波系数区间中的任一滤波系数。其中,该第一预设滤波系数区间为预先设定的滤波系数区间,例如可以设定该第一预设滤波系数区间为1.01至2,本领域技术人员可以理解的是,还可以预先设定该第一预设滤波系数区间为除上述1.01至2之外的其它滤波系数区间,例如1.05至1.95,或者1.2至1.90等等,本领域技术人员对此不作特别的限制。在较优的一实施例中,该第一预设滤波系数区间为1.01至2。
其中,第一预设滤波系数为预先设定的滤波系数,且为第一预设滤波系数区间中的任一滤波系数。例如,在第一预设滤波系数区间为1.01至2时,可以预先设定该第一预设滤波系数为1.5,本领域技术人员可以理解的是,还可以预先设定该第一预设滤波系数为除上述1.5之外的其它相对误差值,例如1.25,或者1.75等等,本领域技术人员对此不作特别的限制。
另外,在本发明的一实施例中,第一滤波窗口在各个方向上的节点数相等且均为奇数,第一滤波窗口在各个方向上的节点数不超过预设节点数阈值。具体的,假设第一滤波窗口为Wx×Wy×Wz,其中,Wx为第一滤波窗口在x方向上的滤波窗口节点数,Wy为第一滤波窗口在y方向上的滤波窗口节点数,Wz为第一滤波窗口在z方向上的滤波窗口节点数。其中,在本发明的一实施例中,x方向上的滤波窗口节点数Wx、y方向上的滤波窗口节点数Wy,以及z方向上的滤波窗口节点数Wz均为奇数,例如均为3、5或7等奇数值,且各个方向上的节点数不超过预设节点数阈值。另外,在本发明的一实施例中,滤波窗口节点数Wx、Wy及Wz均相等。据此,在较优的一实施例中,第一滤波窗口Wx×Wy×Wz=3×3×3。
其中,该预设节点数阈值为预先设定的节点数阈值,例如预先设定该预设节点数阈值为7,本领域技术人员可以理解的是,还可以预先设定该预设节点数阈值为除上述7之外的其它节点数,例如预先设定该预设节点数阈值为5或者9,本发明实施例对此不作特别的限制。一般来讲,该第一滤波窗口在各个方向上的节点数不宜超过7。
据此,在采用第一滤波窗口对每个第一插值视电阻率进行滑动平均滤波,得到第一视电阻率后,继而对每个第一视电阻率进行反插值,得到对应每个第二视电阻率。该对应每个第二视电阻率即为每个原始视电阻率在原始位置的滤波结果,该每个第二视电阻率中数据点的位置与原始视电阻率中数据点的位置是对应的。
鉴于三维插值后的数据是规则的,相邻数据点之间在沿某个坐标轴方向上的分布是等间距的,空间上相邻最近的四个数据点可以组成一个四面体。具体的,在对每个第一视电阻率进行反插值时,根据每个第一视电阻率中某个散点数据的位置,搜索该散点数据所在的四面体,再根据四面体四个顶点上散点数据的数值,推算出该散点数据的数值,进而通过该种反插值方式计算得到所有散点位置的散点数据。本领域技术人员可以理解的是,此处对每个第一视电阻率进行反插值,还可以采用除上述方法之外的其它反插值算法,本发明实施例对此不作特别的限制。
在本发明实施例中,利用规则插值网格对每个原始视电阻率进行克里金插值得到对应每个第一插值视电阻率,进而采用第一滤波窗口对每个第一插值视电阻率进行滑动平均滤波,得到对应每个第一视电阻率,最后对每个第一视电阻率进行反插值得到原始视电阻率在原始位置的滤波结果每个第二视电阻率,依次对每个原始视电阻率进行规则插值网格的克里金插值、余弦滤波及反插值等处理,能够进一步提高滤波效果,进而能够进一步提高静位移校正的效果。
图4示出了本发明实施例提供的静位移校正方法中步骤103的实现流程,为便于描述,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了在进一步提高静位移校正的效果的同时,进一步提高静位移校正的效率,如图4所示,步骤103,根据每个原始视电阻率与每个第二视电阻率中每个对应位置的视电阻率的相对误差,利用相对误差不小于预设相对误差阈值的每个第二视电阻率中的视电阻率替换每个原始视电阻率中对应位置的视电阻率,得到每个第三视电阻率,包括:
步骤401,确定每个原始视电阻率与对应每个第二视电阻率中每个对应位置的视电阻率的相对误差;
步骤402,在每个对应位置的视电阻率的相对误差小于预设相对误差阈值时,将每个原始视电阻率中对应位置的视电阻率作为每个第三视电阻率对应位置的视电阻率;
步骤403,在每个对应位置的视电阻率的相对误差不小于预设相对误差阈值时,将每个第二视电阻率中对应位置的视电阻率作为每个第三视电阻率对应位置的视电阻率。
为了在消除每个原始视电阻率中畸变点,进一步提交静位移校正的效果的同时,提高消除畸变点的效率,进而提高静位移校正的效率,可以利用相对误差不小于预设相对误差阈值的每个第二视电阻率中的数据点替换每个原始视电阻率中对应位置的数据点,得到消除畸变点后的每个第三视电阻率。
鉴于每个原始视电阻率中的数据点的位置与每个第二视电阻率中的数据点的位置是对应的,首先可以确定每个原始视电阻率与每个第二视电阻率中每个对应位置的数据点的相对误差,相对误差反映了数据点的畸变程度。
具体的,在本发明的一实施例中,可以通过如下公式确定每个原始视电阻率中数据点A1与每个第二视电阻率中与该数据点A1位置对应的数据点A2之间的相对误差:
其中,E表示每个原始视电阻率中数据点A1与第二视电阻率中数据点A2之间的相对误差。
在本发明的一实施例中,以每个原始视电阻率为基础,在每个对应位置的数据点的相对误差小于预设相对误差阈值时,认为每个原始视电阻率的该数据点未产生畸变,则保留每个原始视电阻率的该数据点。在每个对应位置的数据点的相对误差不小于预设相对误差阈值时,认为每个原始视电阻率的该数据点产生了畸变,此时利用相对误差不小于预设相对误差阈值的每个第二视电阻率中的数据点替换每个原始视电阻率中对应位置的产生了畸变的数据点,以此形成消除了畸变的第三视电阻率,达到消除每个原始视电阻率中的畸变点的目的,从而提高静位移校正的效果。
在本发明的一实施例中,以每个原始视电阻率为基础,在每个对应位置的数据点的相对误差小于预设相对误差阈值时,认为每个原始视电阻率的该数据点未产生畸变,将每个原始视电阻率中对应位置的数据点作为每个第三视电阻率对应位置的数据点;在每个对应位置的数据点的相对误差不小于预设相对误差阈值时,认为每个原始视电阻率的该数据点产生了畸变,此时将每个第二视电阻率中与每个原始视电阻率的该数据点对应位置的数据点作为每个第三视电阻率对应位置的数据点,以此形成消除了畸变的每个第三视电阻率,达到消除每个原始视电阻率中的畸变点的目的,从而提高静位移校正的效果。
在本发明实施例中,确定每个原始视电阻率与第二视电阻率中每个对应位置的数据点的相对误差,进而在每个对应位置的数据点的相对误差小于预设相对误差阈值时,将每个原始视电阻率中对应位置的数据点作为每个第三视电阻率对应位置的数据点,在每个对应位置的数据点的相对误差不小于预设相对误差阈值时,将每个第二视电阻率中对应位置的数据点作为每个第三视电阻率对应位置的数据点,基于每个数据点的相对误差消除原始视电阻率中的畸变点,不仅能够进一步提高静位移校正的效果,还能够在进一步提高静位移校正的效果的同时,进一步提高静位移校正的效率。
图5示出了本发明实施例提供的静位移校正方法中步骤104的实现流程,为便于描述,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效率,如图5所示,步骤104,确定第三视电阻率数据体中每个第四视电阻率,包括:
步骤501,确定每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值,形成第二视电阻率数据体;
步骤502,将原始视电阻率数据体与第二视电阻率数据体中每个对应测点的视电阻率相加,得到第三视电阻率数据体;
步骤503,确定第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率;
步骤504,将每个测点中多个预设高频段频率号的视电阻率的平均视电阻率作为每个第四视电阻率。
为了进一步提高静位移校正的效率,在确定第四视电阻率时可以先基于每个原始视电阻率与对应每个第三视电阻率,及第二视电阻率数据体的视电阻率,确定第三视电阻率数据体,进而确定第三视电阻率数据体中的每个第四视电阻率。
具体的,首先确定每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值,获得所有原始视电阻率与对应所有第三视电阻率之间的差值,进而形成第二视电阻率数据体。进而,将原始视电阻率数据体与第二视电阻率数据体中每个对应测点的视电阻率相加,即可得到第三视电阻率数据体。可以理解的是,该第二视电阻率数据体及第三视电阻率数据体均为三维视电阻率数据。
至此,在确定第三视电阻率数据体后,获取第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率,继而计算该每个测点中多个预设高频段频率号的视电阻率的平均视电阻率,该平均视电阻率即为第三视电阻率数据体中每个测点对应的第四视电阻率。可以理解的是,该第四视电阻率为平面视电阻率数据。
其中,该多个预设高频段频率号为预先设定的多个高频段的频率号,例如预先设定该预设高频段频率号分别为高频段的频率号为2、3及4的频率号,本领域技术人员可以理解的是,还可以预先设定该预设高频段频率号为除上述频率号2、频率号3及频率号4的频率号之外的其它频率号,本领域技术人员对此不作特别的限制。在较优的一实施例中,该多个预设高频段频率号为高频段中频率号为2、3及4的频率号。
据此,假设该多个预设高频段频率号为高频段中频率号为A、B及C的频率号,即频率号A、频率号B及频率号C。则每个第四视电阻率为每个测点中频率号A、频率号B及频率号C的视电阻率的平均视电阻率,以此确定第三视电阻率数据体中每个测点对应的每个第四视电阻率。
在本发明实施例中,基于每个原始视电阻率与对应每个第三视电阻率,及第二视电阻率数据体的视电阻率,确定第三视电阻率数据体,进而将第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率的平均视电阻率作为第四视电阻率,基于多个高频段频率号的视电阻率的平均值,能够快速确定第四视电阻率,从而进一步提高静位移校正的效率。
图6示出了本发明实施例提供的静位移校正方法中步骤105的实现流程,为便于描述,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效果,如图6所示,步骤105,对每个第四视电阻率进行滤波,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率,包括:
步骤601,采用规则插值网格对每个第四视电阻率进行克里金插值,得到每个第二插值视电阻率;规则插值网格的网格间隔与第四视电阻率的数据间隔的比值为预设倍数区间内的任一倍数;
步骤602,采用第二滤波窗口对每个第二插值视电阻率进行高斯低通滤波,得到每个第五视电阻率;
步骤603,对每个第五视电阻率进行反插值,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率。
为了进一步提高静位移校正的效果,在得到每个第四视电阻率后,进一步依次对每个第四视电阻率进行规则插值网格的克里金插值,得到每个第四视电阻率在原始位置的滤波结果,即每个第五视电阻率。
具体的,首先依然采用规则插值网格对每个第三视电阻率进行克里金插值,得到插值后的每个第二插值视电阻率。此处规则插值网格与图2及图2对应实施例中的规则插值网格相同,具体可参照上述对应实施例中的描述,此处不在详细赘述。
本发明实施例在得到每个第二插值视电阻率后,对每个第二插值视电阻率进行高斯低通滤波,得到对应的每个第四视电阻率。在本发明的一实施例中,该高斯低通滤波的第二滤波系数为第二预设滤波系数区间中任一滤波系数。其中,该第二预设滤波系数区间为预先设定的滤波系数区间,例如可以设定该第二预设滤波系数区间为0.01至10,本领域技术人员可以理解的是,还可以预先设定该第二预设滤波系数区间为除上述0.01至10之外的其它滤波系数区间,例如0.05至9.95,或者0.2至9.90等等,本领域技术人员对此不作特别的限制。在较优的一实施例中,该第二预设滤波系数区间为0.01至10。
其中,第二预设滤波系数为预先设定的滤波系数,且为第二预设滤波系数区间中的任一滤波系数。例如,在第二预设滤波系数区间为0.01至10时,可以预先设定该第二预设滤波系数为4.5,本领域技术人员可以理解的是,还可以预先设定该第二预设滤波系数为除上述4.5之外的其它相对误差值,例如4.1,或者5.5等等,本领域技术人员对此不作特别的限制。
另外,在本发明的一实施例中,高斯低通滤波在各个方向上的节点数相等且均为奇数,高斯低通滤波在各个方向上的节点数不超过预设节点数阈值。具体的,假设高斯低通滤波为Wx×Wy×Wz,其中,Wx为高斯低通滤波在x方向上的滤波节点数,Wy为高斯低通滤波在y方向上的滤波节点数,Wz为高斯低通滤波在z方向上的滤波节点数。其中,在本发明的一实施例中,x方向上的滤波节点数Wx、y方向上的滤波节点数Wy,以及z方向上的滤波节点数Wz均为奇数,例如均为3、5或7等奇数值,且各个方向上的节点数不超过预设节点数阈值。另外,在本发明的一实施例中,滤波节点数Wx、Wy及Wz均相等。据此,在较优的一实施例中,高斯低通滤波为Wx×Wy×Wz=3×3×3。
其中,该预设节点数阈值为预先设定的节点数阈值,例如预先设定该预设节点数阈值为7,本领域技术人员可以理解的是,还可以预先设定该预设节点数阈值为除上述7之外的其它节点数,例如预先设定该预设节点数阈值为5或者9,本发明实施例对此不作特别的限制。一般来讲,该高斯低通滤波在各个方向上的节点数不宜超过7。
另外,在本发明的一实施例中,通过如下公式确定高斯低通滤波中每个数据的权系数:
W=exp(-factor×(fxi×fxi+fyj×fyj+fzk×fzk));
fxi=i/Wx,i=-Wx/2~Wx/2;
fyj=j/Wy,j=-Wy/2~Wy/2;
fzk=k/Wz,k=-Wz/2~Wz/2;
其中,W表示高斯低通滤波中数据(i,j,k)的权系数,i表示数据(i,j,k)的在x方向上的数据节点编号,j表示数据(i,j,k)在y方向上的数据节点编号,k表示数据(i,j,k)在z方向上的数据节点编号,Wx表示高斯低通滤波在x方向上的节点数,Wy表示高斯低通滤波在y方向上的节点数,Wz表示高斯低通滤波在z方向上的节点数,fxi表示在x方向上的数据节点i的权系数,fyj表示在y方向上的数据节点j的权系数,fzk表示在z方向上的数据节点k的权系数,factor表示高斯低通滤波的滤波系数。
据此,在对每个第二插值视电阻率进行高斯低通滤波,得到对应的每个第五视电阻率后,继而对每个第五视电阻率进行反插值,得到每个第六视电阻率。该每个第六视电阻率即为第四视电阻率在原始位置的滤波结果,该第六视电阻率中数据点的位置与第四视电阻率中数据点的位置是对应的。可以理解的是,该每个第四视电阻率、每个第五视电阻率及每个第六视电阻率均为平面视电阻率数据。
鉴于三维插值后的数据是规则的,相邻数据点之间在沿某个坐标轴方向上的分布是等间距的,空间上相邻最近的四个数据点可以组成一个四面体。具体的,在对每个第五视电阻率进行反插值时,根据每个第五视电阻率中某个散点数据的位置,搜索该散点数据所在的四面体,再根据四面体四个顶点上散点数据的数值,推算出该散点数据的数值,进而通过该种反插值方式计算得到所有散点位置的散点数据,得到每个第六视电阻率。本领域技术人员可以理解的是,此处对第四视电阻率进行反插值,还可以采用除上述方法之外的其它反插值算法,本发明实施例对此不作特别的限制。
在本发明实施例中,利用规则插值网格对每个第四视电阻率进行克里金插值,得到插值后的每个第二插值视电阻率,进而对每个第二插值视电阻率进行高斯低通滤波,得到对应的每个第五视电阻率,最后对每个第五视电阻率进行反插值处理,得到对应的每个第五视电阻率在原始位置的滤波结果第六视电阻率,依次对每个第四视电阻率进行规则插值网格的克里金插值、高斯低通滤波及反插值处理,能够进一步提高滤波效果,进而进一步提高静位移校正的效果。
图7示出了本发明实施例提供的静位移校正方法中步骤106的实现流程,为便于描述,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效率,如图7所示,步骤106,根据每个第四视电阻率与对应每个第六视电阻率,及第二视电阻率数据体,确定对原始视电阻率数据体静校正后的第四视电阻率数据体,包括:
步骤701,确定每个第四视电阻率与对应每个第六视电阻率中每个对应测点之间的差值,形成第三视电阻率数据体;
步骤702,将第二视电阻率数据体与第三视电阻率数据体中每个对应测点的视电阻率相加,得到第四视电阻率数据体。
在确定每个第四视电阻率对应的每个第六视电阻率,即可根据每个第四视电阻率,对应的每个第六视电阻率,及第二视电阻率数据体,确定静位移校正后的第四视电阻率数据体。
具体的,在确定每个第四视电阻率对应的每个第六视电阻率后,首先根据每个第四视电阻率及对应的每个第六视电阻率,确定每个第四视电阻率与对应每个第六视电阻率中每个对应测点之间的差值,所有的第四视电阻率与对应第六视电阻率中每个对应测点之间的差值形成第三视电阻率数据体。在确定第三视电阻率数据体后,将第二视电阻率数据体与第三视电阻率数据体中每个对应测点的视电阻率相加,得到第四视电阻率数据体。该第四视电阻率数据体即为对原始视电阻率数据体进行静位移校正后的视电阻率数据。可以理解的是,该第三视电阻率数据体及第四视电阻率数据体均为三维视电阻率数据。
在本发明实施例中,确定每个第四视电阻率与对应每个第六视电阻率中每个对应测点之间的差值,形成第三视电阻率数据体,进而将第二视电阻率数据体与第三视电阻率数据体中每个对应测点的视电阻率相加,得到静位移校正后的第四视电阻率数据体,能够快速确定静位移校正后的视电阻率数据,因而能够进一步提高静位移校正的效率。
图8示出了本发明实施例提供的静位移校正方法的另一实现流程,为便于描述,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效率,如图8所示,在上述方法步骤的基础上,在步骤101之前,静位移校正方法,还包括:
步骤801,对原始视电阻率数据体取以预设数值为底的对数,得到第一视电阻率数据体;
相应的,步骤102,确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率,包括:
步骤802,确定包括预设方向上所有视电阻率数据的第一视电阻率数据体中每个原始视电阻率。
鉴于在静位移校正处理中采用对数坐标,能够提高静位移校正的效率。因此,为了进一步提高静位移校正的效率,在步骤101之前首先对原始视电阻率数据体进行对数计算,即对原始视电阻率数据体取以预设数值为底的对数,得到对数计算后的第一视电阻率数据体。可以理解的是,该第一视电阻率数据体为三维视电阻率数据。
其中,该预设数值为预先设定的数值,例如预先设定该预设数值为10,本领域技术人员可以理解的是,还可以预先设定该预设数值为除上述10之外的其它数值,例如预先设定该预设数值为5或者8,又或者12,本发明实施例对此不作特别的限制。在较优的一实施例中,该预设数值为10。
在对原始视电阻率数据体进行上述对数处理后,相应的则步骤101中确定原始视电阻率数据体中的每个原始视电阻率,具体为确定第一视电阻率数据体中的每个原始视电阻率。
在本发明实施例中,对原始视电阻率数据体取以预设数值为底的对数,得到第一视电阻率数据体,在对数坐标中进行静位移校正处理,能够进一步提高静位移校正的效率。
在本发明的一实施例中,为了进一步提高静位移校正的效率,在上述方法步骤的基础上,在步骤106之后,静位移校正方法,还包括:
步骤803,对第四视电阻率数据体进行反对数运算,确定原始视电阻率数据体校正后的第五视电阻率数据体。
在进行静位移校正之前,对原始视电阻率数据体进行对数运算,转化为对数坐标中进行静位移校正处理。那么,对应的在得到原始视电阻率数据体校正后的第四视电阻率数据体,该第四视电阻率数据体为对数坐标中的静位移校正结果,因此需要将对数形式表示的第五视电阻率数据体中的视电阻率数据,进行与对数运算对应的反对数运算,以形成非对数坐标形式的、静位移校正后的第五视电阻率数据体。
在本发明实施例中,对第四视电阻率数据体进行反对数运算,确定原始视电阻率数据体校正后的第五视电阻率数据体,能够进一步提高静位移校正的效率。
本发明实施例还提供一种静位移校正装置,如下面的实施例所述。由于这些装置解决问题的原理与静位移校正方法相似,因此这些装置的实施可以参见方法的实施,重复之处不再赘述。
图9示出了本发明实施例提供的静位移校正装置的功能模块,为便于说明,仅示出了与本发明实施例相关的部分,详述如下:
参考图9,所述静位移校正装置所包含的各个模块用于执行图1对应实施例中的各个步骤,具体请参阅图1以及图1对应实施例中的相关描述,此处不再赘述。本发明实施例中,所述静位移校正装置包括原始视电阻率确定模块901、第一滤波模块902、替换模块903、第四视电阻率确定模块904、第二滤波模块905及校正结果获得模块906。
原始视电阻率确定模块901,用于确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率;原始视电阻率为原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率。
第一滤波模块902,用于对每个原始视电阻率进行滤波,得到每个原始视电阻率在原始位置的每个第二视电阻率。
替换模块903,用于根据每个原始视电阻率与每个第二视电阻率中每个对应位置的视电阻率的相对误差,利用相对误差不小于预设相对误差阈值的每个第二视电阻率中的视电阻率替换每个原始视电阻率中对应位置的视电阻率,得到每个第三视电阻率;视电阻率的相对误差反映视电阻率的畸变程度。
第四视电阻率确定模块904,用于确定第三视电阻率数据体中每个第四视电阻率;第三视电阻率数据体为每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值形成的第二视电阻率数据体,与原始视电阻率数据体中每个对应测点的视电阻率相加得到的;每个第四视电阻率为第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率的平均视电阻率。
第二滤波模块905,用于对每个第四视电阻率进行滤波,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率。
校正结果获得模块906,用于根据每个第四视电阻率与对应每个第六视电阻率,及第二视电阻率数据体,确定对原始视电阻率数据体静校正后的第四视电阻率数据体。
在本发明实施例中,替换模块903利用相对误差不小于预设相对误差阈值的第二视电阻率中的视电阻率替换原始视电阻率中对应位置的视电阻率,能够基于相对误差消除视电阻率的畸变点,从而提高对视电阻率的静位移校正效果。另外,本发明实施例只需根据视电阻率的相对误差,即可消除视电阻率中的畸变点,无需其它额外的处理即可快速消除畸变点,能够提高静位移校正的效率。
图10示出了本发明实施例提供的静位移校正装置中原始视电阻率确定模块901的结构示意,为便于说明,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效率,参考图10,所述原始视电阻率确定模块901所包含的各个单元用于执行图2对应实施例中的各个步骤,具体请参阅图2以及图2对应实施例中的相关描述,此处不再赘述。本发明实施例中,所述原始视电阻率确定模块901包括视电阻率获取单元1001及第一平均单元1002。
视电阻率获取单元1001,用于获取原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率。
第一平均单元1002,用于将每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率作为每个原始视电阻率。
在本发明实施例中,视电阻率获取单元1001获取原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率,在确定原始视电阻率时第一平均单元1002只需要计算每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率,即可快速确定原始视电阻率,因而能够进一步提高静位移校正的效率。
图11示出了本发明实施例提供的静位移校正装置中第一滤波模块902的结构示意,为便于说明,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效果,参考图11,所述第一滤波模块902所包含的各个单元用于执行图3对应实施例中的各个步骤,具体请参阅图3以及图3对应实施例中的相关描述,此处不再赘述。本发明实施例中,所述第一滤波模块902包括第一插值单元1101、滑动平均滤波单元1102及第一反插值单元1103。
第一插值单元1101,用于采用规则插值网格对每个原始视电阻率进行克里金插值,得到每个第一插值视电阻率;规则插值网格的网格间隔与原始视电阻率的数据间隔的比值为预设倍数区间内的任一倍数。
滑动平均滤波单元1102,用于采用第一滤波窗口对每个第一插值视电阻率进行滑动平均滤波,得到每个第一视电阻率;第一滤波窗口的中心位置数据不参与滤波。
第一反插值单元1103,用于对每个第一视电阻率进行反插值,得到每个原始视电阻率在原始位置的滤波结果每个第二视电阻率。
在本发明实施例中,第一插值单元1101利用规则插值网格对每个原始视电阻率进行克里金插值得到对应每个第一插值视电阻率,进而滑动平均滤波单元1102采用第一滤波窗口对每个第一插值视电阻率进行滑动平均滤波,得到对应每个第一视电阻率,最后第一反插值单元1103对每个第一视电阻率进行反插值得到原始视电阻率在原始位置的滤波结果每个第二视电阻率,依次对每个原始视电阻率进行规则插值网格的克里金插值、余弦滤波及反插值等处理,能够进一步提高滤波效果,进而能够进一步提高静位移校正的效果。
图12示出了本发明实施例提供的静位移校正装置中替换模块903的结构示意,为便于说明,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了在进一步提高静位移校正的效果的同时,进一步提高静位移校正的效率,参考图12,所述替换模块903所包含的各个单元用于执行图4对应实施例中的各个步骤,具体请参阅图4以及图4对应实施例中的相关描述,此处不再赘述。本发明实施例中,所述替换模块903包括相对误差确定单元1201、保留单元1202及替换单元1203。
相对误差确定单元1201,用于确定每个原始视电阻率与对应每个第二视电阻率中每个对应位置的视电阻率的相对误差。
保留单元1202,用于在每个对应位置的视电阻率的相对误差小于预设相对误差阈值时,将每个原始视电阻率中对应位置的视电阻率作为每个第三视电阻率对应位置的视电阻率。
替换单元1203,用于在每个对应位置的视电阻率的相对误差不小于预设相对误差阈值时,将每个第二视电阻率中对应位置的视电阻率作为每个第三视电阻率对应位置的视电阻率。
在本发明实施例中,相对误差确定单元1201确定每个原始视电阻率与第二视电阻率中每个对应位置的数据点的相对误差,进而保留单元1202在每个对应位置的数据点的相对误差小于预设相对误差阈值时,将每个原始视电阻率中对应位置的数据点作为每个第三视电阻率对应位置的数据点,替换单元1203在每个对应位置的数据点的相对误差不小于预设相对误差阈值时,将每个第二视电阻率中对应位置的数据点作为每个第三视电阻率对应位置的数据点,基于每个数据点的相对误差消除原始视电阻率中的畸变点,不仅能够进一步提高静位移校正的效果,还能够在进一步提高静位移校正的效果的同时,进一步提高静位移校正的效率。
图13示出了本发明实施例提供的静位移校正装置中第四视电阻率确定模块904的结构示意,为便于说明,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效率,参考图6,所述第四视电阻率确定模块904所包含的各个单元用于执行图5对应实施例中的各个步骤,具体请参阅图5以及图5对应实施例中的相关描述,此处不再赘述。本发明实施例中,所述第四视电阻率确定模块904包括第一差值确定单元1301、第一相加单元1302、视电阻率确定单元1303及第二平均单元1304。
第一差值确定单元1301,用于确定每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值,形成第二视电阻率数据体。
第一相加单元1302,用于将原始视电阻率数据体与第二视电阻率数据体中每个对应测点的视电阻率相加,得到第三视电阻率数据体。
视电阻率确定单元1303,用于确定第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率。
第二平均单元1304,用于将每个测点中多个预设高频段频率号的视电阻率的平均视电阻率作为每个第四视电阻率。
在本发明实施例中,第一差值确定单元1301及第一相加单元1302基于每个原始视电阻率与对应每个第三视电阻率,及第二视电阻率数据体的视电阻率,确定第三视电阻率数据体,进而第二平均单元1304将第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率的平均视电阻率作为第四视电阻率,基于多个高频段频率号的视电阻率的平均值,能够快速确定第四视电阻率,从而进一步提高静位移校正的效率。
图14示出了本发明实施例提供的静位移校正装置中第二滤波模块905的结构示意,为便于说明,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效果,参考图14,所述第二滤波模块905所包含的各个单元用于执行图6对应实施例中的各个步骤,具体请参阅图6以及图6对应实施例中的相关描述,此处不再赘述。本发明实施例中,所述第二滤波模块905包括第二插值单元1401、低通滤波单元1402及第二反插值单元1403。
第二插值单元1401,用于采用规则插值网格对每个第四视电阻率进行克里金插值,得到每个第二插值视电阻率;规则插值网格的网格间隔与第四视电阻率的数据间隔的比值为预设倍数区间内的任一倍数。
低通滤波单元1402,用于采用第二滤波窗口对每个第二插值视电阻率进行高斯低通滤波,得到每个第五视电阻率。
第二反插值单元1403,用于对每个第五视电阻率进行反插值,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率。
在本发明实施例中,第二插值单元1401利用规则插值网格对每个第四视电阻率进行克里金插值,得到插值后的每个第二插值视电阻率,进而低通滤波单元1402对每个第二插值视电阻率进行高斯低通滤波,得到对应的每个第五视电阻率,最后第二反插值单元1403对每个第五视电阻率进行反插值处理,得到对应的每个第五视电阻率在原始位置的滤波结果第六视电阻率,依次对每个第四视电阻率进行规则插值网格的克里金插值、高斯低通滤波及反插值处理,能够进一步提高滤波效果,进而进一步提高静位移校正的效果。
图15示出了本发明实施例提供的静位移校正装置中校正结果获得模块906的结构示意,为便于说明,仅示出了与本发明实施例相关的部分,详述如下:
在本发明的一实施例中,为了进一步提高静位移校正的效率,参考图15,所述校正结果获得模块906所包含的各个单元用于执行图7对应实施例中的各个步骤,具体请参阅图7以及图7对应实施例中的相关描述,此处不再赘述。本发明实施例中,所述校正结果获得模块906包括第二差值确定单元1501及第二相加单元1502。
第二差值确定单元1501,用于确定每个第四视电阻率与对应每个第六视电阻率中每个对应测点之间的差值,形成第三视电阻率数据体。
第二相加单元1502,用于将第二视电阻率数据体与第三视电阻率数据体中每个对应测点的视电阻率相加,得到第四视电阻率数据体。
在本发明实施例中,第二差值确定单元1501确定每个第四视电阻率与对应每个第六视电阻率中每个对应测点之间的差值,形成第三视电阻率数据体,进而第二相加单元1502将第二视电阻率数据体与第三视电阻率数据体中每个对应测点的视电阻率相加,得到静位移校正后的第四视电阻率数据体,能够快速确定静位移校正后的视电阻率数据,因而能够进一步提高静位移校正的效率。
图16示出了本发明实施例提供的静位移校正装置的另一功能模块,为便于说明,仅示出了与本发明实施例相关的部分,详述如下:
参考图16,所述静位移校正装置所包含的各个模块用于执行图8对应实施例中的各个步骤,具体请参阅图8以及图8对应实施例中的相关描述,此处不再赘述。本发明实施例中,为了进一步提高静位移校正的效率,在上述模块结构的基础上,所述静位移校正装置,还包括对数模块1601。
对数模块1601,用于对原始视电阻率数据体取以预设数值为底的对数,得到第一视电阻率数据体。
相应的,原始视电阻率确定模块901,具体用于确定包括预设方向上所有视电阻率数据的第一视电阻率数据体中每个原始视电阻率。
在本发明实施例中,对数模块1601对原始视电阻率数据体取以预设数值为底的对数,得到第一视电阻率数据体,在对数坐标中进行静位移校正处理,能够进一步提高静位移校正的效率。
在本发明的一实施例中,为了进一步提高静位移校正的效率,如图16所示,在上述模块结构的基础上,静位移校正装置,还包括反对数模块1603。
反对数模块1603,用于对第四视电阻率数据体进行反对数运算,确定原始视电阻率数据体校正后的第五视电阻率数据体。
在本发明实施例中,反对数模块1603对第四视电阻率数据体进行反对数运算,确定原始视电阻率数据体校正后的第五视电阻率数据体,能够进一步提高静位移校正的效率。
图17示出了本发明实施例的静位移校正前XY方向上1.035Hz频率的视电阻率平面等值线示意,图18示出了本发明实施例的静位移校正前XY方向上1.035Hz频率的视电阻率平面等值线示意,为便于说明,仅示出了与本发明实施例相关的部分,详述如下:
如图17所示,该图为静位移校正前XY方向上1.035Hz频率的视电阻率平面等值线,图中横坐标为x(m)为地理南北坐标,纵坐标y(m)为地理东西坐标。从图17中可以看出,图中存在大量的由静态位移现象引起的畸变点,表现为高频的高异常或者低异常,图中的畸变点会对整体数据产生不可控的不良影响。
如图18所示,该图为利用本发明提供的静位移校正方法对原始视电阻率数据体进行静位移校正后,XY方向上1.035Hz频率的视电阻率平面等值线示意,从图18中可以看出,原始视电阻率中的畸变点被有效消除,等值线平滑。同时保留了数据的整体趋势,符合预期。
对比图17及图18,可以得出本发明提供的静位移校正方法能够有效消除三维视电阻率中的畸变点,提高静位移校正效果,保留数据的整体趋势,恢复了数据的正常分布,对于弱异常也有很好的保留效果,静位移校正效果优异。
图19示出了本发明实施例的静位移校正前XY方向上某号测线的的视电阻率剖面等值线示意,图20示出了本发明实施例的静位移校正后XY方向上某号测线的视电阻率剖面等值线示意,为便于说明,仅示出了与本发明实施例相关的部分,详述如下:
如图19所示,图中横坐标为偏移距,纵坐标为频率。从图中可以看出,图19中所示静位移较正前某号测线的的视电阻率剖面上存在大量的由静态位移现象引起的“挂面条”状异常,表现为整体背景中的横向排列的、条带状高频高异常或低异常,使得整体数据上较为凌乱,大量“挂面条”状异常存在会对数据产生不可控的影响。
如图20所示,该图为经过本发明实施例提供静位移校正方法处理后的某号测线的视电阻率平面等值线,图中横坐标为偏移距,纵坐标为频率。从图20中可以看出,静位移校正后的视电阻率平面等值线整体上过渡平滑,数据趋势明显,“挂面条”状异常被有效消除,符合预期。
对比图19及图20,可以得出本发明实施例提供的静位移校正方法,能够有效消除视电阻率中的“挂面条”状异常,保留数据的整体趋势,恢复数据的正常分布,对于弱异常也有很好的保留效果,静位移校正效果优异。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述静位移校正方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上静位移校正方法的计算机程序。
综上所述,本发明实施例中,利用相对误差不小于预设相对误差阈值的第二视电阻率中的视电阻率替换原始视电阻率中对应位置的视电阻率,能够基于相对误差消除视电阻率的畸变点,从而提高对视电阻率的静位移校正效果。另外,本发明实施例只需根据视电阻率的相对误差,即可消除视电阻率中的畸变点,无需其它额外的处理即可快速消除畸变点,能够提高静位移校正效率。
本领域内的技术人员应明白,本发明的实施例可提供为方法、***、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (16)
1.一种静位移校正方法,其特征在于,包括:
确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率;原始视电阻率为原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率;
对每个原始视电阻率进行滤波,得到每个原始视电阻率在原始位置的每个第二视电阻率;
根据每个原始视电阻率与每个第二视电阻率中每个对应位置的视电阻率的相对误差,利用相对误差不小于预设相对误差阈值的每个第二视电阻率中的视电阻率替换每个原始视电阻率中对应位置的视电阻率,得到每个第三视电阻率;视电阻率的相对误差反映视电阻率的畸变程度;
确定第三视电阻率数据体中每个第四视电阻率;第三视电阻率数据体为每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值形成的第二视电阻率数据体,与原始视电阻率数据体中每个对应测点的视电阻率相加得到的;每个第四视电阻率为第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率的平均视电阻率;
对每个第四视电阻率进行滤波,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率;
根据每个第四视电阻率与对应每个第六视电阻率,及第二视电阻率数据体,确定对原始视电阻率数据体静校正后的第四视电阻率数据体;
其中,对每个原始视电阻率进行滤波,得到每个原始视电阻率在原始位置的每个第二视电阻率,包括:
采用规则插值网格对每个原始视电阻率进行克里金插值,得到每个第一插值视电阻率;规则插值网格的网格间隔与原始视电阻率的数据间隔的比值为预设倍数区间内的任一倍数;
采用第一滤波窗口对每个第一插值视电阻率进行滑动平均滤波,得到每个第一视电阻率;第一滤波窗口的中心位置数据不参与滤波;
对每个第一视电阻率进行反插值,得到每个原始视电阻率在原始位置的滤波结果每个第二视电阻率;
确定第三视电阻率数据体中每个第四视电阻率,包括:
确定每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值,形成第二视电阻率数据体;
将原始视电阻率数据体与第二视电阻率数据体中每个对应测点的视电阻率相加,得到第三视电阻率数据体;
确定第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率;
将每个测点中多个预设高频段频率号的视电阻率的平均视电阻率作为每个第四视电阻率。
2.如权利要求1所述的方法,其特征在于,确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率,包括:
获取原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率;
将每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率作为每个原始视电阻率。
3.如权利要求1所述的方法,其特征在于,根据每个原始视电阻率与每个第二视电阻率中每个对应位置的视电阻率的相对误差,利用相对误差不小于预设相对误差阈值的每个第二视电阻率中的视电阻率替换每个原始视电阻率中对应位置的视电阻率,得到每个第三视电阻率,包括:
确定每个原始视电阻率与对应每个第二视电阻率中每个对应位置的视电阻率的相对误差;
在每个对应位置的视电阻率的相对误差小于预设相对误差阈值时,将每个原始视电阻率中对应位置的视电阻率作为每个第三视电阻率对应位置的视电阻率;
在每个对应位置的视电阻率的相对误差不小于预设相对误差阈值时,将每个第二视电阻率中对应位置的视电阻率作为每个第三视电阻率对应位置的视电阻率。
4.如权利要求1所述的方法,其特征在于,对每个第四视电阻率进行滤波,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率,包括:
采用规则插值网格对每个第四视电阻率进行克里金插值,得到每个第二插值视电阻率;规则插值网格的网格间隔与第四视电阻率的数据间隔的比值为预设倍数区间内的任一倍数;
采用第二滤波窗口对每个第二插值视电阻率进行高斯低通滤波,得到每个第五视电阻率;
对每个第五视电阻率进行反插值,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率。
5.如权利要求1所述的方法,其特征在于,根据每个第四视电阻率与对应每个第六视电阻率,及第二视电阻率数据体,确定对原始视电阻率数据体静校正后的第四视电阻率数据体,包括:
确定每个第四视电阻率与对应每个第六视电阻率中每个对应测点之间的差值,形成第三视电阻率数据体;
将第二视电阻率数据体与第三视电阻率数据体中每个对应测点的视电阻率相加,得到第四视电阻率数据体。
6.如权利要求1所述的方法,其特征在于,在确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率之前,还包括:
对原始视电阻率数据体取以预设数值为底的对数,得到第一视电阻率数据体;
相应的,确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率,包括:
确定包括预设方向上所有视电阻率数据的第一视电阻率数据体中每个原始视电阻率。
7.如权利要求6所述的方法,其特征在于,在根据每个第四视电阻率与对应每个第六视电阻率,及第二视电阻率数据体,确定对原始视电阻率数据体静校正后的第四视电阻率数据体之后,还包括:
对第四视电阻率数据体进行反对数运算,确定原始视电阻率数据体校正后的第五视电阻率数据体。
8.一种静位移校正装置,其特征在于,包括:
原始视电阻率确定模块,用于确定包括预设方向上所有视电阻率数据的原始视电阻率数据体中每个原始视电阻率;原始视电阻率为原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率;
第一滤波模块,用于对每个原始视电阻率进行滤波,得到每个原始视电阻率在原始位置的每个第二视电阻率;
替换模块,用于根据每个原始视电阻率与每个第二视电阻率中每个对应位置的视电阻率的相对误差,利用相对误差不小于预设相对误差阈值的每个第二视电阻率中的视电阻率替换每个原始视电阻率中对应位置的视电阻率,得到每个第三视电阻率;视电阻率的相对误差反映视电阻率的畸变程度;
第四视电阻率确定模块,用于确定第三视电阻率数据体中每个第四视电阻率;第三视电阻率数据体为每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值形成的第二视电阻率数据体,与原始视电阻率数据体中每个对应测点的视电阻率相加得到的;每个第四视电阻率为第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率的平均视电阻率;
第二滤波模块,用于对每个第四视电阻率进行滤波,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率;
校正结果获得模块,用于根据每个第四视电阻率与对应每个第六视电阻率,及第二视电阻率数据体,确定对原始视电阻率数据体静校正后的第四视电阻率数据体;
其中,第一滤波模块包括:
第一插值单元,用于采用规则插值网格对每个原始视电阻率进行克里金插值,得到每个第一插值视电阻率;规则插值网格的网格间隔与原始视电阻率的数据间隔的比值为预设倍数区间内的任一倍数;
滑动平均滤波单元,用于采用第一滤波窗口对每个第一插值视电阻率进行滑动平均滤波,得到每个第一视电阻率;第一滤波窗口的中心位置数据不参与滤波;
第一反插值单元,用于对每个第一视电阻率进行反插值,得到每个原始视电阻率在原始位置的滤波结果每个第二视电阻率;
第四视电阻率确定模块包括:
第一差值确定单元,用于确定每个原始视电阻率与对应每个第三视电阻率中每个对应测点之间的差值,形成第二视电阻率数据体;
第一相加单元,用于将原始视电阻率数据体与第二视电阻率数据体中每个对应测点的视电阻率相加,得到第三视电阻率数据体;
视电阻率确定单元,用于确定第三视电阻率数据体中每个测点中多个预设高频段频率号的视电阻率;
第二平均单元,用于将每个测点中多个预设高频段频率号的视电阻率的平均视电阻率作为每个第四视电阻率。
9.如权利要求8所述的静位移校正装置,其特征在于,原始视电阻率确定模块包括:
视电阻率获取单元,用于获取原始视电阻率数据体中每个测点中与预设频率邻近的预设多个中频段频率的视电阻率;
第一平均单元,用于将每个测点中与预设频率邻近的预设多个中频段频率的视电阻率的平均视电阻率作为每个原始视电阻率。
10.如权利要求8所述的静位移校正装置,其特征在于,替换模块包括:
相对误差确定单元,用于确定每个原始视电阻率与对应每个第二视电阻率中每个对应位置的视电阻率的相对误差;
保留单元,用于在每个对应位置的视电阻率的相对误差小于预设相对误差阈值时,将每个原始视电阻率中对应位置的视电阻率作为每个第三视电阻率对应位置的视电阻率;
替换单元,用于在每个对应位置的视电阻率的相对误差不小于预设相对误差阈值时,将每个第二视电阻率中对应位置的视电阻率作为每个第三视电阻率对应位置的视电阻率。
11.如权利要求8所述的静位移校正装置,其特征在于,第二滤波模块包括:
第二插值单元,用于采用规则插值网格对每个第四视电阻率进行克里金插值,得到每个第二插值视电阻率;规则插值网格的网格间隔与第四视电阻率的数据间隔的比值为预设倍数区间内的任一倍数;
低通滤波单元,用于采用第二滤波窗口对每个第二插值视电阻率进行高斯低通滤波,得到每个第五视电阻率;
第二反插值单元,用于对每个第五视电阻率进行反插值,得到每个第四视电阻率在原始位置的滤波结果每个第六视电阻率。
12.如权利要求8所述的静位移校正装置,其特征在于,校正结果获得模块包括:
第二差值确定单元,用于确定每个第四视电阻率与对应每个第六视电阻率中每个对应测点之间的差值,形成第三视电阻率数据体;
第二相加单元,用于将第二视电阻率数据体与第三视电阻率数据体中每个对应测点的视电阻率相加,得到第四视电阻率数据体。
13.如权利要求8所述的静位移校正装置,其特征在于,还包括:
对数模块,用于对原始视电阻率数据体取以预设数值为底的对数,得到第一视电阻率数据体;
相应的,原始视电阻率确定模块,具体用于确定包括预设方向上所有视电阻率数据的第一视电阻率数据体中每个原始视电阻率。
14.如权利要求13所述的静位移校正装置,其特征在于,还包括:
反对数模块,用于对第四视电阻率数据体进行反对数运算,确定原始视电阻率数据体校正后的第五视电阻率数据体。
15.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7任一所述静位移校正方法。
16.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现权利要求1至7任一所述静位移校正方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010045597.7A CN113126172B (zh) | 2020-01-16 | 2020-01-16 | 静位移校正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010045597.7A CN113126172B (zh) | 2020-01-16 | 2020-01-16 | 静位移校正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113126172A CN113126172A (zh) | 2021-07-16 |
CN113126172B true CN113126172B (zh) | 2024-01-30 |
Family
ID=76772055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010045597.7A Active CN113126172B (zh) | 2020-01-16 | 2020-01-16 | 静位移校正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113126172B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113484920B (zh) * | 2021-08-17 | 2023-05-19 | 成都理工大学 | 一种频域电磁测深资料二维结构化反演方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009124446A1 (zh) * | 2008-04-10 | 2009-10-15 | 中国石油集团东方地球物理勘探有限责任公司 | 三维小面元电磁连续阵列数据采集方法 |
CN104267443A (zh) * | 2014-05-22 | 2015-01-07 | 中国地质科学院矿产资源研究所 | 基于反演模型的大地电磁场静位移校正方法 |
CN104903750A (zh) * | 2014-07-23 | 2015-09-09 | 王雅苹 | 一种电磁数据的保幅静校正方法 |
CN106443189A (zh) * | 2016-12-05 | 2017-02-22 | 国家电网公司 | 一种接地极极址及其周边土壤电阻率三维探测方法和*** |
CN106646666A (zh) * | 2017-01-16 | 2017-05-10 | 中南大学 | 基于平面波电磁测深的静态效应校正方法 |
-
2020
- 2020-01-16 CN CN202010045597.7A patent/CN113126172B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009124446A1 (zh) * | 2008-04-10 | 2009-10-15 | 中国石油集团东方地球物理勘探有限责任公司 | 三维小面元电磁连续阵列数据采集方法 |
CN104267443A (zh) * | 2014-05-22 | 2015-01-07 | 中国地质科学院矿产资源研究所 | 基于反演模型的大地电磁场静位移校正方法 |
CN104903750A (zh) * | 2014-07-23 | 2015-09-09 | 王雅苹 | 一种电磁数据的保幅静校正方法 |
CN106443189A (zh) * | 2016-12-05 | 2017-02-22 | 国家电网公司 | 一种接地极极址及其周边土壤电阻率三维探测方法和*** |
CN106646666A (zh) * | 2017-01-16 | 2017-05-10 | 中南大学 | 基于平面波电磁测深的静态效应校正方法 |
Non-Patent Citations (4)
Title |
---|
CSAMT静态效应校正方法研究与应用;杨妮妮;王志宏;赵晓鸣;;江西科学(第02期);全文 * |
大地电磁三维静态位移校正方法及其应用效果;吴炳良, 邵敏;勘探地球物理进展(第03期);全文 * |
平面聚类静态校正法;许建荣;石油地球物理勘探(第06期);全文 * |
杨妮妮 ; 王志宏 ; 赵晓鸣 ; .CSAMT静态效应校正方法研究与应用.江西科学.2009,(第02期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN113126172A (zh) | 2021-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | A new finite element gradient recovery method: superconvergence property | |
Bortolotti et al. | Uniform Penalty inversion of two-dimensional NMR relaxation data | |
CN106443587B (zh) | 一种高分辨率的快速反卷积声源成像算法 | |
Amjad et al. | Assessment of three-dimensional density measurements from tomographic background-oriented schlieren (BOS) | |
CN113126172B (zh) | 静位移校正方法及装置 | |
CN113945969B (zh) | 不连续地质体检测方法、装置及电子设备 | |
CN107102359B (zh) | 地震数据保幅重建方法和*** | |
Boyle et al. | Methods for calculating the electrode position Jacobian for impedance imaging | |
CN110631499B (zh) | 基于无网格方法的数字图像变形测量方法及电子设备 | |
Hyvönen et al. | Compensation for geometric modeling errors by positioning of electrodes in electrical impedance tomography | |
CN107102056B (zh) | 缺陷漏磁信号的单元伸缩构建方法 | |
CN111291316A (zh) | 一种基于小波变换的多尺度电阻率反演方法及*** | |
Harrach | Interpolation of missing electrode data in electrical impedance tomography | |
CN113627027B (zh) | 非平凡各向异性介质电磁场数值模拟方法及*** | |
CN113504568B (zh) | 一种基于小生境差分进化算法的中值滤波方法 | |
CN112797917B (zh) | 一种高精度的数字散斑干涉相位定量测量方法 | |
Alotaibi et al. | Restoration of Blurred Images Corrupted by Impulse Noise via Median Filters and ℓ p-ℓ q Minimization | |
Nara | Algebraic reconstruction of the general-order poles of a meromorphic function | |
CN117078791B (zh) | 一种ct环形伪影校正方法、装置、电子设备及存储介质 | |
CN113050187A (zh) | 滤波方法及装置、计算机设备及计算机可读存储介质 | |
US20230103493A1 (en) | Method and system for creating dipole moment model | |
CN112230279B (zh) | 增强纵波地震数据品质的方法及装置 | |
CN112799126B (zh) | 沿起伏地表的地震数据重构方法、装置、介质及设备 | |
Li et al. | Determination of three-dimensional atomic positions from tomographic reconstruction using ensemble empirical mode decomposition | |
CN115393547A (zh) | 一种月球卫星重力异常数据的全方向滤波方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |