CN113866834A - 一种基于熵滤波的场源中心位置反演方法 - Google Patents
一种基于熵滤波的场源中心位置反演方法 Download PDFInfo
- Publication number
- CN113866834A CN113866834A CN202111078344.0A CN202111078344A CN113866834A CN 113866834 A CN113866834 A CN 113866834A CN 202111078344 A CN202111078344 A CN 202111078344A CN 113866834 A CN113866834 A CN 113866834A
- Authority
- CN
- China
- Prior art keywords
- window
- entropy
- anomaly
- filtering
- magnetic
- 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.)
- Pending
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 238000000926 separation method Methods 0.000 claims abstract description 14
- 238000003384 imaging method Methods 0.000 claims abstract description 11
- 230000002159 abnormal effect Effects 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 8
- 238000012800 visualization Methods 0.000 description 6
- 230000005415 magnetization Effects 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/081—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the magnetic field is produced by the objects or geological structures
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于地质资源勘探领域,是一种基于熵滤波的场源中心位置反演方法,步骤为:S1:获取实测磁异常数据;S2:根据已反演的实测磁异常数据,利用熵滤波异常分离方法分离出对应地下不同深度层位的磁异常,熵滤波异常分离方法是空间域滑动窗口滤波方法,分离出的各层剩余异常局部极大值水平位置即为场源中心的水平位置;S3:对S2得到的各层结果置于相应地下深度层位,进行三维成像。该三维数据局部极值位置即场源的中心位置。有益效果:本发明方法提高了场源中心位置的计算精度,具有更高的分辨率及定位效果,且方法简便易行。
Description
技术领域
本发明涉及地球物理地质资源勘探技术领域,具体涉及一种基于熵滤波的场源中心位置获取方法。
背景技术
磁法勘探是地球物理勘探方法中发展最早,应用较广的一种方法。因其效率高、成本低、工作领域广、应用范围广等特点受到人们关注。在磁异常反演与解释工作中,由于受到纵向分辨率的限制,难以准确定位出场源的中心位置。同时,当地下存在多个场源时,场源之间相互干扰也会影响场源中心的定位精度。除此之外,也会受到噪声干扰等影响。因此,寻找一种能够准确定位出场源中心位置的方法是磁异常反演与解释中的关键问题之一。
目前常用的磁异常反演方法主要有两种:物性反演法,通过将地下半空间剖分为规则的网格单元,反演各个网格单元的物性信息,进而通过各网格的物性分布来推断地下场源的分布情况。由于物性反演法具有多解性,需给出特定的约束条件,这些问题制约着该方法的反演效果。欧拉反褶积法,以欧拉齐次方程为基础,应用不同形态的地质体所特有的构造指数来反演地下场源的分布,计算较为简便,且无需大量的约束条件,但由于在实际工作中往往缺少先验信息,难以确定地下场源的构造指数,因此制约了该方法的反演效果。
当前磁异常分离方法众多,但其中大多数都难以从众多异常中准确地定位出场源的中心位置。异常分离依旧是磁异常数据处理与解释中的一个难题。因此,如何进一步提高成像的分辨率,改善成像效果,进而更加准确地定位出场源的中心位置成为该领域的努力方向。
发明内容
为解决现有技术存在的问题,本发明提出了基于熵滤波的场源中心位置反演方法,应用于获得地质资源勘探中场源的中心位置。
本发明的技术方案:一种基于熵滤波的场源中心位置反演方法,包括以下步骤:
S1:获取实测磁异常数据,所述实测磁异常数据为目标体磁场数据;
S2:根据已获得的实测磁异常数据,利用熵滤波异常分离方法分离出对应地下不同深度层位的磁异常,其中所述熵滤波异常分离方法是一种空间域滑动窗口滤波方法,分离出的各层剩余异常局部极大值水平位置即为场源中心的水平位置;
S3:对S2得到的各层结果置于相应地下深度层位,进行三维成像,三维成像数据局部极值位置即得到场源的中心位置。
上述空间域滑动窗口滤波方法如下:
b、计算窗口内各点磁异常与磁异常平均值之差的绝对值的最大值;
c、对窗口内各点计算各自的熵值,利用各点熵值计算各点熵权系数;
e、滑动窗口至下一待滤波点并重复上述步骤。
所述步骤a中计算窗口内磁异常最大值,最小值及平均值分别为:
所述步骤b中窗口内各点磁异常与磁异常平均值之差的绝对值的最大值为:
所述步骤c中窗口内各点熵值计算为:
可得窗口内各点熵权系数为:
所述步骤d中熵滤波计算的区域异常结果为:
因实施本发明需要的计算量大且需要可视化显示,故需要计算机辅助计算。所需硬件设备包括:一个或多个处理器、用于存储一个或多个计算程序的存储装置,当一个或多个计算程序被所述一个或多个处理器执行,使得所述一个或多个处理器执行本发明方法的有关数据计算。
本发明的有益效果:本方法利用熵滤波分离异常,相对于传统的场源中心反演方法,具有无需约束条件及先验信息等优势,可以快速定位出场源的中心位置;本发明相较于其他传统方法,具有更高的分辨率及定位效果。利用本发明可以进行场源的中心位置高精度计算,有利于地质资源勘探的精度提高。方法简便易行,有利于推广。
附图说明
图1为两个球体模型的示意图;
图2为两个球体模型产生的磁异常平面等值线图;
图3为采用本发明方法对图2的模型磁异常数据进行中心位置计算所得的三维可视化图;
图4为采用本发明方法对图2的模型磁异常数据进行中心位置计算所得的平面可视化图;
图5为两个立方体模型的示意图;
图6为两个立方体模型产生的磁异常平面等值线图;
图7为采用本发明方法对图5的模型磁异常数据进行中心位置计算所得的三维可视化图;
图8为采用本发明方法对图5的模型磁异常数据进行中心位置计算所得的平面可视化图。
具体实施方式
一种基于熵滤波的场源中心位置反演方法,包括以下步骤:
S1:获取实测磁异常数据,所述实测磁异常数据为目标体磁场数据;
S2:根据已获得的实测磁异常数据,利用熵滤波异常分离方法分离出对应地下不同深度层位的磁异常,其中所述熵滤波异常分离方法是在熵的理论基础上提出的一种空间域滑动窗口滤波方法,分离出的各层剩余异常局部极大值水平位置即为场源中心的水平位置;
S3:对S2得到的各层结果置于相应地下深度层位,进行三维成像,三维成像数据局部极值位置即得到场源的中心位置。
上述空间域滑动窗口滤波方法如下:
b、计算窗口内各点磁异常与磁异常平均值之差的绝对值的最大值;
c、对窗口内各点计算各自的熵值,利用各点熵值计算各点熵权系数;
e、滑动窗口至下一待滤波点并重复上述步骤。
所述步骤a中计算窗口内磁异常最大值,最小值及平均值分别为:
所述步骤b中窗口内各点磁异常与磁异常平均值之差的绝对值的最大值为:
所述步骤c中窗口内各点熵值计算为:
可得窗口内各点熵权系数为:
所述步骤d中熵滤波计算的区域异常结果为:
实验例1
如图1所示的两个球体模型的模型示意图,设定磁化倾角为90°,磁化偏角为0°,两个球体模型中心距地表100m,两个球体模型半径为50m,磁化率为0.1SI,产生的磁异常如图2所示。采用本发明方法对原始磁异常数据进行基于熵滤波异常分离的场源中心位置计算,本发明反演的两个球体模型的中心深度均为100m,对计算结果进行可视化得到如图3所示的三维场源中心位置计算结果以及图4所示的平面场源中心切片图。可以看出,本发明可以较准确地反演场源中心位置计算结果,并实现了结果的可视化。
实验例2
如图5所示的两个立方体模型的模型示意图,设定磁化倾角为90°,磁化偏角为0°,两个立方体模型中心距地表200m,两个立方体模型长和宽均为100m,高为200m,磁化率为0.1SI,产生的磁异常如图6所示。采用本发明方法对原始磁异常数据进行基于熵滤波异常分离的场源中心位置计算,本发明反演的两个球体模型的中心深度均为100m,对计算结果进行可视化得到如图7所示的三维场源中心位置计算结果以及图8所示的平面场源中心切片图。可以看出,本发明可以较准确地反演场源中心位置计算结果,并实现了结果的可视化。
Claims (7)
1.一种基于熵滤波的场源中心位置反演方法,其特征在于,包括以下步骤:
S1:获取实测磁异常数据,所述实测磁异常数据为目标体磁场数据;
S2:根据已获得的实测磁异常数据,利用熵滤波异常分离方法分离出对应地下不同深度层位的磁异常,其中所述熵滤波异常分离方法是一种空间域滑动窗口滤波方法,分离出的各层剩余异常局部极大值水平位置即为场源中心的水平位置;
S3:对S2得到的各层结果置于相应地下深度层位,进行三维成像,三维成像数据局部极值位置即得到场源的中心位置。
7.一种基于熵滤波的场源中心位置反演方法的实施装置,其特征在于,含有至少一个处理器、至少一个存储装置,存储装置用于存储计算机程序,处理器用于执行计算机程序,通过计算机实施如下步骤:
S1:输入获取的实测磁异常数据,所述实测磁异常数据为目标体磁场数据;
S2:根据已获得的实测磁异常数据,利用熵滤波异常分离方法分离出对应地下不同深度层位的磁异常,其中所述熵滤波异常分离方法是一种空间域滑动窗口滤波方法,分离出的各层剩余异常局部极大值水平位置即为场源中心的水平位置;
S3:对S2得到的各层结果置于相应地下深度层位,进行三维成像,三维成像数据局部极值位置即得到场源的中心位置;
所述空间域滑动窗口滤波步骤为:
b、计算窗口内各点磁异常与磁异常平均值之差的绝对值的最大值;
c、对窗口内各点计算各自的熵值,利用各点熵值计算各点熵权系数;
e、滑动窗口至下一待滤波点并重复上述步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111078344.0A CN113866834A (zh) | 2021-09-15 | 2021-09-15 | 一种基于熵滤波的场源中心位置反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111078344.0A CN113866834A (zh) | 2021-09-15 | 2021-09-15 | 一种基于熵滤波的场源中心位置反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113866834A true CN113866834A (zh) | 2021-12-31 |
Family
ID=78995936
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111078344.0A Pending CN113866834A (zh) | 2021-09-15 | 2021-09-15 | 一种基于熵滤波的场源中心位置反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113866834A (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040028273A1 (en) * | 2000-12-01 | 2004-02-12 | Shugo Hamahashi | Entropy filter, and area extracting method using the filter |
US20120075440A1 (en) * | 2010-09-28 | 2012-03-29 | Qualcomm Incorporated | Entropy based image separation |
CN110515136A (zh) * | 2019-07-03 | 2019-11-29 | 吉林大学 | 一种基于重磁界面反演的大地热流估计方法 |
CN110687610A (zh) * | 2019-09-19 | 2020-01-14 | 长安大学 | 一种基于重磁数据相关性分析的场源定位及属性识别方法 |
CN110967774A (zh) * | 2019-11-15 | 2020-04-07 | 中国科学院电子学研究所 | 一种基于传感器阵列的磁异常检测方法 |
CN110989005A (zh) * | 2019-12-23 | 2020-04-10 | 西北工业大学青岛研究院 | 基于标量磁力仪阵列的弱磁异常自适应实时探测方法 |
CN111337992A (zh) * | 2020-03-23 | 2020-06-26 | 兰州大学 | 一种基于位场数据向下延拓的场源深度获得方法 |
CN111913227A (zh) * | 2020-08-11 | 2020-11-10 | 哈尔滨工程大学 | 一种磁异常信号频率特征的计算方法和*** |
CN113359192A (zh) * | 2021-06-04 | 2021-09-07 | 北京中科导控科技有限公司 | 一种微弱磁异常目标检测及定位方法 |
-
2021
- 2021-09-15 CN CN202111078344.0A patent/CN113866834A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040028273A1 (en) * | 2000-12-01 | 2004-02-12 | Shugo Hamahashi | Entropy filter, and area extracting method using the filter |
US20120075440A1 (en) * | 2010-09-28 | 2012-03-29 | Qualcomm Incorporated | Entropy based image separation |
CN110515136A (zh) * | 2019-07-03 | 2019-11-29 | 吉林大学 | 一种基于重磁界面反演的大地热流估计方法 |
CN110687610A (zh) * | 2019-09-19 | 2020-01-14 | 长安大学 | 一种基于重磁数据相关性分析的场源定位及属性识别方法 |
CN110967774A (zh) * | 2019-11-15 | 2020-04-07 | 中国科学院电子学研究所 | 一种基于传感器阵列的磁异常检测方法 |
CN110989005A (zh) * | 2019-12-23 | 2020-04-10 | 西北工业大学青岛研究院 | 基于标量磁力仪阵列的弱磁异常自适应实时探测方法 |
CN111337992A (zh) * | 2020-03-23 | 2020-06-26 | 兰州大学 | 一种基于位场数据向下延拓的场源深度获得方法 |
CN111913227A (zh) * | 2020-08-11 | 2020-11-10 | 哈尔滨工程大学 | 一种磁异常信号频率特征的计算方法和*** |
CN113359192A (zh) * | 2021-06-04 | 2021-09-07 | 北京中科导控科技有限公司 | 一种微弱磁异常目标检测及定位方法 |
Non-Patent Citations (1)
Title |
---|
郭良辉: "磁异常∆T三维相关成像", 《地球物理学报》, vol. 53, no. 2, pages 437 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dannemiller et al. | A new method for determination of magnetization direction | |
Fedi et al. | Understanding imaging methods for potential field data | |
EP1877834B1 (en) | 3D fast fault restoration | |
CN106777598B (zh) | 任意磁化率分布复杂磁性体磁场梯度张量数值模拟方法 | |
CN111337992B (zh) | 一种基于位场数据向下延拓的场源深度获得方法 | |
Ren et al. | Gravity gradient tensor of arbitrary 3D polyhedral bodies with up to third-order polynomial horizontal and vertical mass contrasts | |
CN110687610B (zh) | 一种基于重磁数据相关性分析的场源定位及属性识别方法 | |
CN111856599B (zh) | 一种基于pde的磁测数据等效源化极与类型转换方法 | |
CN110346835B (zh) | 大地电磁的正演方法、正演***、存储介质及电子设备 | |
KR20080023946A (ko) | 오일러 해를 이용한 지하공동의 3차원 중력역산 방법 및이를 이용한 3차원 영상화 방법 | |
CN110346834B (zh) | 三维频率域可控源电磁的正演方法、*** | |
CN113962077B (zh) | 三维各向异性强磁场数值模拟方法、装置、设备及介质 | |
Yangsheng et al. | Three-dimensional fractal distribution of the number of rock-mass fracture surfaces and its simulation technology | |
CN109254327B (zh) | 三维强磁性体的勘探方法及勘探*** | |
CN115292973B (zh) | 一种任意采样的空间波数域三维磁场数值模拟方法及*** | |
CN110703347A (zh) | 基于构造背景的重力断裂影像识别方法 | |
CN114004127A (zh) | 二维主轴各向异性强磁场数值模拟方法、装置、设备及介质 | |
CN110660051B (zh) | 一种基于导航金字塔的张量投票处理方法 | |
Li et al. | Fast 3D forward modeling of the magnetic field and gradient tensor on an undulated surface | |
Haber et al. | An adaptive mesh method for electromagnetic inverse problems | |
CN112668146A (zh) | 一种基于欧拉反褶积法实用性改进的场源位置估算方法 | |
CN117034719A (zh) | 重磁数据区域场去除方法 | |
CN113866834A (zh) | 一种基于熵滤波的场源中心位置反演方法 | |
Wang et al. | Interpretation of gravity data for fault distribution near the Mongolia–Hinggan metallogenic belt in the eastern China-Mongolia frontier area | |
Vatankhah et al. | A method for 2-dimensional inversion of gravity data |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211231 |