CN114488326A - 一种提高探测航磁数据地质体边界能力的方法及系统 - Google Patents

一种提高探测航磁数据地质体边界能力的方法及系统 Download PDF

Info

Publication number
CN114488326A
CN114488326A CN202210136957.3A CN202210136957A CN114488326A CN 114488326 A CN114488326 A CN 114488326A CN 202210136957 A CN202210136957 A CN 202210136957A CN 114488326 A CN114488326 A CN 114488326A
Authority
CN
China
Prior art keywords
boundary
matrix
data
aeromagnetic
magnetic gradient
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.)
Granted
Application number
CN202210136957.3A
Other languages
English (en)
Other versions
CN114488326B (zh
Inventor
王林飞
王明
熊盛青
林晓星
屈进红
薛典军
李芳�
何辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Aero Geophysical Survey and Remote Sensing Center for Natural Resources
Original Assignee
China Aero Geophysical Survey and Remote Sensing Center for Natural Resources
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Aero Geophysical Survey and Remote Sensing Center for Natural Resources filed Critical China Aero Geophysical Survey and Remote Sensing Center for Natural Resources
Priority to CN202210136957.3A priority Critical patent/CN114488326B/zh
Publication of CN114488326A publication Critical patent/CN114488326A/zh
Application granted granted Critical
Publication of CN114488326B publication Critical patent/CN114488326B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/38—Processing data, e.g. for analysis, for interpretation, for correction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Computing Systems (AREA)
  • Geophysics (AREA)
  • Remote Sensing (AREA)
  • Algebra (AREA)
  • Geology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Measuring Magnetic Variables (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Harvesting Machines For Root Crops (AREA)

Abstract

本发明属于航空磁测领域，具体提供一种提高探测航磁数据地质体边界能力的方法及系统，其中方法包括S1，获取实测或计算的全磁力梯度张量数据矩阵M；S2，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3；S3，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；S4，建立边界检测函数E，并获取E的最大值max(|E|)；E＝λ1·λ2·λ3·A；S5，计算δ·max(|E|)进而得到新的边界检测均衡滤波器iLP，其中，δ为均衡深浅部异常的调节系数。该方案具有较高的分辨率、较强的抗倾斜磁化能力和抗噪声干扰性，不仅能够均衡不同振幅大小航磁数据目标地质体的边界信息，而且还能避免产生多余的虚假边界异常。提高了计算稳定性、增强信噪比，还能提高横向分辨力。

Description

一种提高探测航磁数据地质体边界能力的方法及系统
技术领域
本发明涉及航空磁测领域，更具体地，涉及一种提高探测航磁数据地质体边界能力的方法及系统。
背景技术
航空磁场测量数据是不同深度、不同形态、不同规模的磁性地质体磁场信息在观测面上的综合反映；由于磁场的叠加效应，使得某些具有一定地质意义的异常变得复杂，在原始图件上很难识别，给地质解释工作带来了难度。随着工程技术与磁梯度张量探测仪器研发技术的不断地发展与成熟，应用磁力张量数据在分析、处理上述问题也得到了相应的发展。磁力张量数据是磁场矢量分量的梯度，包含了的磁场信息，能够反映目标体的矢量磁矩信息，张量数据具有高精度、高分辨率、多参量的优点，可用于描述场源体的磁化方向和几何形态，提高对目标地质体的分辨率。
虽然局部相位法(Local Phase，LP)可以探测航磁数据目标地质体的边界，突出更多有用信息，但是该方法缺点是数值稳定性较弱和易产生虚假异常边界，且当分母
Figure BDA0003504752240000011
或接近0时，局部相位法存在“解析奇点”，会使得计算结果不稳定，影响实际应用效果。
发明内容
本发明针对现有技术中存在的局部相位法探测航磁数据目标地质体的边界时，数值稳定性较弱且易产生虚假异常边界的技术问题。
本发明提供了一种提高探测航磁数据地质体边界能力的方法，包括以下步骤：
S1，获取实测或计算的全磁力梯度张量数据矩阵M；
S2，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3；
S3，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；
S4，建立边界检测函数E，并获取E的最大值max(|E|)；
E＝λ1·λ2·λ3·A；
S5，计算δ·max(|E|)进而得到新的边界检测均衡滤波器iLP，其中，δ为均衡深浅部异常的调节系数。
优选地，所述全磁力梯度张量数据矩阵的特征值λ1、λ2、λ3均对应着航磁数据目标地质体的边界；全磁力梯度张量总模值A包含了全部9个张量元素的信息，其最大值对应着地质体的边界。
优选地，所述δ的默认值为0.001。
优选地，所述全磁力梯度张量数据矩阵M的矩阵表达形式为：
Figure BDA0003504752240000021
根据三维直角坐标系确定磁场数据在x,y,z三个方向的磁场分量以及每个磁场分量分别在x,y,z方向的一阶梯度分量，得到12个分量数据，其中9个一阶梯度分量(Mxx、Mxy、Mxz、Myx、Myy、Myz、Mzx、Mzy、Mzz)构成全张量磁梯度数据。
优选地，所述新的均衡边界识别滤波器iLP，表达式为：
Figure BDA0003504752240000031
其中，δ是均衡深浅部异常的调节系数，取值为0到1之间任意值。
优选地，所述总模值计算公式为：
Figure BDA0003504752240000032
本发明还提供了一种提高探测航磁数据地质体边界能力的系统，包括：
全磁力梯度张量数据获取模块，用于获取实测或计算的全磁力梯度张量数据矩阵M；
计算模块，用于根据全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3，然后计算得到该矩阵的总模值A；建立边界检测函数E，并获取E的最大值max(|E|)；E＝λ1·λ2·λ3·A；计算δ·max(|E|)进而得到新的边界检测均衡滤波器iLP，其中，δ为均衡深浅部异常的调节系数。
本发明还提供了一种电子设备，包括存储器、处理器，所述处理器用于执行存储器中存储的计算机管理类程序时实现提高探测航磁数据地质体边界能力的方法的步骤。
本发明还提供了一种计算机可读存储介质，其上存储有计算机管理类程序，所述计算机管理类程序被处理器执行时实现提高探测航磁数据地质体边界能力的方法的步骤。
有益效果：本发明提供的一种提高探测航磁数据地质体边界能力的方法及系统，其中方法包括S1，获取实测或计算的全磁力梯度张量数据矩阵M；S2，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3；S3，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；S4，建立边界检测函数E，并获取E的最大值max(|E|)；E＝λ1·λ2·λ3·A；S5，计算δ·max(|E|)进而得到新的边界检测均衡滤波器iLP，其中，δ为均衡深浅部异常的调节系数。该方案具有较高的分辨率、较强的抗倾斜磁化能力和抗噪声干扰性，不仅能够均衡不同振幅大小航磁数据目标地质体的边界信息，而且还能避免产生多余的虚假边界异常，解决了现有局部相位法存在“解析奇点”问题，提高计算稳定性、增强信噪比，还能提高横向分辨力。
附图说明
图1为本发明提供的一种提高探测航磁数据地质体边界能力的方法流程图；
图2为本发明提供的一种可能的电子设备的硬件结构示意图；
图3为本发明提供的一种可能的计算机可读存储介质的硬件结构示意图；
图4为本发明提供的一种提高探测航磁数据地质体边界能力的方法的结果图；
图5为现有LP法计算结果图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
传统方案在实际航空磁测数据处理工作中利用局部相位法检测目标地质体边界时，存在其存在精度低、稳定性和抗噪能力差、不能有效的均衡深浅异常的振幅，易产生虚假的航磁异常边界结果的问题。基于此，参考图1，为本发明提供的一种提高探测航磁数据地质体边界能力的方法及系统，其中方法包括S1，获取实测或计算的全磁力梯度张量数据矩阵M；S2，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3；S3，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；S4，建立边界检测函数E，并获取E的最大值max(|E|)；E＝λ1·λ2·λ3·A；S5，计算δ·max(|E|)进而得到新的边界检测均衡滤波器iLP，其中，δ为均衡深浅部异常的调节系数。该方案具有较高的分辨率、较强的抗倾斜磁化能力和抗噪声干扰性，不仅能够均衡不同振幅大小航磁数据目标地质体的边界信息，而且还能避免产生多余的虚假边界异常，解决了现有局部相位法存在“解析奇点”问题，提高计算稳定性、增强信噪比，还能提高横向分辨力。
在一个具体的实施场景中，提高局部相位法探测航磁数据地质体边界能力的方法技术，包括以下步骤：
(1)获取实测或计算的全磁力梯度张量数据矩阵M；
(2)根据步骤(1)中的全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3，全磁力梯度张量数据矩阵的特征值λ1、λ2、λ3均对应着航磁数据目标地质体的边界。
(3)根据步骤(1)中的全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；全磁力梯度张量总模值A包含了全部9个张量元素的信息，其最大值对应着地质体的边界。
(4)根据步骤(2)中全磁力梯度张量数据矩阵的特征值λ1、λ2、λ3和步骤(3)中全磁力梯度张量数据矩阵的总模值A，建立边界检测函数E；因E同时具备了张量特征值和总模值的属性，能够提高对浅部目标地质体的识别精度，但检测深部目标地质体的分辨率低，所以需要进一步提高检测深部目标地质体的边界的能力。
E＝λ1·λ2·λ3·A。
(5)根据步骤(4)和实测或计算获得的磁梯度张量数据值M，在现有局部相位法(LP)基础上，重新定义并构造一个新的均衡边界识别滤波器iLP，表达式为：
Figure BDA0003504752240000061
其中，δ是均衡深浅部异常的调节系数，一般取值为0-1。
与现有的局部相位法(LP)表达式对比：
Figure BDA0003504752240000062
其中，T代表是航磁总场值；
Figure BDA0003504752240000063
分别是T在x、y、z三个方向一阶导数(或一阶梯度)。
进一步的，如上所述的一种新的提高局部相位法探测航磁数据地质体边界能力的方法技术，步骤(1)中，全磁力梯度张量数据矩阵M是航空磁场在不同方向分量的变化率，矩阵表达形式为：
Figure BDA0003504752240000071
根据三维直角坐标系确定磁场数据在x,y,z三个方向的磁场分量以及每个磁场分量分别在x,y,z方向的一阶梯度分量，得到12个分量数据，其中9个一阶梯度分量构成全张量磁梯度数据。
进一步的方案，如上所述的一种新的提高局部相位法探测航磁数据地质体边界能力的方法技术，步骤(2)和步骤(3)中，计算全磁力梯度张量数据矩阵M的特征值λ1、λ2、λ3和总模值A。
Figure BDA0003504752240000072
进一步的，如上所述的一种新的提高局部相位法探测航磁数据地质体边界能力的方法技术，步骤(4)中，建立一个新的边界检测函数E。
E＝λ1·λ2·λ3·A。
进一步的，如上所述的一种新的提高局部相位法探测航磁数据地质体边界能力的方法技术，步骤(5)中，δ是均衡深浅部异常的调节系数，取值为0-1。
本发明与现有技术相比，有益效果和优势在于：
本发明基于全磁力梯度张量数据，突破了现有局部相位法定义的思想，新构建了一个合理的均衡航磁数据目标地质体边界检测方法，能更好的探测出埋深不同的多源场物体的边界，使边界识别结果更收敛，且该方法有效地避免了磁化方向和噪声对结果的干扰，解决了现有局部相位法存在“解析奇点”问题，提高计算稳定性和消除产生虚假航磁数据地质体边界干扰、增强信噪比，以及提高了不同埋深的地质体的边界位置增强和提取能力，具有更高的分辨率和精度。
为了验证方法的有效性和技术优势，对其进行了实测数据计算，并与现有局部相位法对比。如图4和图5所示，从图中可以看出，本方案提高了不同埋深的地质体的边界位置增强和提取能力，具有更高的分辨率和精度。
可以理解的是，本发明提供的一种提高探测航磁数据地质体边界能力的系统，与前述各实施例提供的提高探测航磁数据地质体边界能力的相关技术特征一样，在此不再赘述。
请参阅图2为本发明实施例提供的电子设备的实施例示意图。如图2所示，本发明实施例提了一种电子设备，包括存储器1310、处理器1320及存储在存储器1310上并可在处理器1320上运行的计算机程序1311，处理器1320执行计算机程序1311时实现以下步骤：S1，获取实测或计算的全磁力梯度张量数据矩阵M；
S2，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3；
S3，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；
S4，建立边界检测函数E，并获取E的最大值max(|E|)；
E＝λ1·λ2·λ3·A；
S5，计算δ·max(|E|)进而得到新的边界检测均衡滤波器iLP，其中，δ为均衡深浅部异常的调节系数。
请参阅图3为本发明提供的一种计算机可读存储介质的实施例示意图。如图3所示，本实施例提供了一种计算机可读存储介质1400，其上存储有计算机程序1411，该计算机程序1411被处理器执行时实现如下步骤：S1，获取实测或计算的全磁力梯度张量数据矩阵M；
S2，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3；
S3，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；
S4，建立边界检测函数E，并获取E的最大值max(|E|)；
E＝λ1·λ2·λ3·A；
S5，计算δ·max(|E|)进而得到新的边界检测均衡滤波器iLP，其中，δ为均衡深浅部异常的调节系数。
需要说明的是，在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详细描述的部分，可以参见其它实施例的相关描述。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式计算机或者其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包括这些改动和变型在内。

Claims (9)

1.一种提高探测航磁数据地质体边界能力的方法，其特征在于，包括以下步骤：
S1，获取实测或计算的全磁力梯度张量数据矩阵M；
S2，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3；
S3，根据步骤S1中的全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；
S4，建立边界检测函数E，并获取E的最大值max(|E|)；
E＝λ1·λ2·λ3·A；
S5，计算δ·max(|E|)进而得到新的边界检测均衡滤波器iLP，其中，δ为均衡深浅部异常的调节系数。
2.根据权利要求1所述的提高探测航磁数据地质体边界能力的方法，其特征在于，所述全磁力梯度张量数据矩阵的特征值λ1、λ2、λ3均对应着航磁数据目标地质体的边界；全磁力梯度张量总模值A包含了全部9个张量元素的信息，其最大值对应着地质体的边界。
3.根据权利要求1所述的提高探测航磁数据地质体边界能力的方法，其特征在于，所述δ的默认值为0.001。
4.根据权利要求1所述的提高探测航磁数据地质体边界能力的方法，其特征在于，所述全磁力梯度张量数据矩阵M的矩阵表达形式为：
Figure FDA0003504752230000021
根据三维直角坐标系确定磁场数据在x,y,z三个方向的磁场分量以及每个磁场分量分别在x,y,z方向的一阶梯度分量，得到12个分量数据，其中9个一阶梯度分量(Mxx、Mxy、Mxz、Myx、Myy、Myz、Mzx、Mzy、Mzz)构成全张量磁梯度数据。
5.根据权利要求4所述的提高探测航磁数据地质体边界能力的方法，其特征在于，所述新的均衡边界识别滤波器iLP，表达式为：
Figure FDA0003504752230000022
其中，δ是均衡深浅部异常的调节系数，取值为0到1之间任意值。
6.根据权利要求1所述的提高探测航磁数据地质体边界能力的方法，其特征在于，所述总模值计算公式为：
Figure FDA0003504752230000023
7.一种提高探测航磁数据地质体边界能力的系统，其特征在于，包括：
全磁力梯度张量数据获取模块，用于获取实测或计算的全磁力梯度张量数据矩阵M；
计算模块，用于根据全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3，然后计算得到该矩阵的总模值A；建立边界检测函数E，并获取E的最大值max(|E|)；E＝λ1·λ2·λ3·A；计算δ·max(|E|)进而得到新的边界检测均衡滤波器iLP，其中，δ为均衡深浅部异常的调节系数。
8.一种电子设备，其特征在于，包括存储器、处理器，所述处理器用于执行存储器中存储的计算机管理类程序时实现如权利要求1-7任一项所述的提高探测航磁数据地质体边界能力的方法的步骤。
9.一种计算机可读存储介质，其特征在于，其上存储有计算机管理类程序，所述计算机管理类程序被处理器执行时实现如权利要求1-7任一项所述的提高探测航磁数据地质体边界能力的方法的步骤。
CN202210136957.3A 2022-02-15 2022-02-15 一种提高探测航磁数据地质体边界能力的方法及系统 Active CN114488326B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210136957.3A CN114488326B (zh) 2022-02-15 2022-02-15 一种提高探测航磁数据地质体边界能力的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210136957.3A CN114488326B (zh) 2022-02-15 2022-02-15 一种提高探测航磁数据地质体边界能力的方法及系统

Publications (2)

Publication Number Publication Date
CN114488326A true CN114488326A (zh) 2022-05-13
CN114488326B CN114488326B (zh) 2023-07-28

Family

ID=81481136

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210136957.3A Active CN114488326B (zh) 2022-02-15 2022-02-15 一种提高探测航磁数据地质体边界能力的方法及系统

Country Status (1)

Country Link
CN (1) CN114488326B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115236755A (zh) * 2022-07-25 2022-10-25 中国自然资源航空物探遥感中心 基于张量特征值的航磁异常边界检测方法、装置
CN116360000A (zh) * 2023-02-17 2023-06-30 中国自然资源航空物探遥感中心 一种地壳物质结构的航空磁法探测技术

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0818363D0 (en) * 2008-10-07 2008-11-12 Gravitec Instr Ltd Gradiometer for measuring gravitational and magnetic field gradients with improved sensor
GB201118371D0 (en) * 2007-08-02 2011-12-07 Arkex Ltd Geophysical data processing systems
CN112749493A (zh) * 2021-01-25 2021-05-04 中国自然资源航空物探遥感中心 基于全磁力梯度张量特征值的地质体边界检测方法及系统
CN113886753A (zh) * 2021-10-09 2022-01-04 中国自然资源航空物探遥感中心 基于张量特征值的Tilt法航磁边界检测方法、装置
CN113917544A (zh) * 2021-10-08 2022-01-11 中国科学院空天信息创新研究院 基于磁梯度张量特征值的近地表目标位置快速圈定方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB201118371D0 (en) * 2007-08-02 2011-12-07 Arkex Ltd Geophysical data processing systems
GB0818363D0 (en) * 2008-10-07 2008-11-12 Gravitec Instr Ltd Gradiometer for measuring gravitational and magnetic field gradients with improved sensor
CN112749493A (zh) * 2021-01-25 2021-05-04 中国自然资源航空物探遥感中心 基于全磁力梯度张量特征值的地质体边界检测方法及系统
CN113917544A (zh) * 2021-10-08 2022-01-11 中国科学院空天信息创新研究院 基于磁梯度张量特征值的近地表目标位置快速圈定方法
CN113886753A (zh) * 2021-10-09 2022-01-04 中国自然资源航空物探遥感中心 基于张量特征值的Tilt法航磁边界检测方法、装置

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
GUO,CW;XING,Z;WANG.LF;MA,Y;HUAN,HF: "Three-directional analytic signal analysis and interpretation of magnetic gradient tensor.", APPLIED GEOPHYSICS *
吴招才;刘天佑;: "磁力梯度张量测量及应用", 地质科技情报, no. 03 *
姚刚;董向欣;李丽丽;马国庆;孟令顺;: "东海陆架盆地构造划分的高精度重磁解释技术", 吉林大学学报(地球科学版), no. 02 *
田招招: "基于重磁及梯度张量数据的边界识别与位置反演研究", 中国优秀硕士学位论文全文数据库 *
郑强: "基于磁力梯度全张量的均衡边界识别方法研究", 中国优秀硕士学位论文全文数据库 *
郑强;郭华;张贵宾;韩松;王明;刘浩军;: "基于磁力梯度全张量特征值的均衡边界识别方法", 石油地球物理勘探, no. 02 *
陈安国;周涛发;刘东甲;张舒;: "基于Theta法改进的均衡滤波器在位场边界识别中的应用――以庐枞矿集区为例", 地球物理学报, no. 02 *
魏永强;倪卫冲;房江奇;: "边界增强法在航磁数据处理中的应用", 物探与化探, no. 01 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115236755A (zh) * 2022-07-25 2022-10-25 中国自然资源航空物探遥感中心 基于张量特征值的航磁异常边界检测方法、装置
CN115236755B (zh) * 2022-07-25 2023-10-03 中国自然资源航空物探遥感中心 基于张量特征值的航磁异常边界检测方法、装置
CN116360000A (zh) * 2023-02-17 2023-06-30 中国自然资源航空物探遥感中心 一种地壳物质结构的航空磁法探测技术
CN116360000B (zh) * 2023-02-17 2023-09-05 中国自然资源航空物探遥感中心 一种地壳物质结构的航空磁法探测方法

Also Published As

Publication number Publication date
CN114488326B (zh) 2023-07-28

Similar Documents

Publication Publication Date Title
CN110095738B (zh) 磁梯度仪参数设计方法
CN114488326B (zh) 一种提高探测航磁数据地质体边界能力的方法及系统
CA2686215C (en) Determining borehole corrected formation properties
CN113886753B (zh) 基于张量特征值的Tilt法航磁边界检测方法、装置
CN108333628A (zh) 基于正则化约束的弹性波最小二乘逆时偏移方法
Yu et al. Acoustic-elastic coupled equation for ocean bottom seismic data elastic reverse time migration
CN112749493A (zh) 基于全磁力梯度张量特征值的地质体边界检测方法及系统
CN110133716B (zh) 基于组合模型加权函数的磁异常数据三维反演方法
US10877175B2 (en) Seismic acquisition geometry full-waveform inversion
CN111399066B (zh) 一种基于正交基函数处理标量磁异常梯度信号的方法
CN110187393B (zh) 一种基于广义回归神经网络的航磁补偿方法
CN110361792B (zh) 一种地球物理数据融合及成像方法、介质与设备
CN107479091B (zh) 一种提取逆时偏移角道集的方法
CN112668146B (zh) 一种基于欧拉反褶积法实用性改进的场源位置估算方法
US20210396901A1 (en) Acoustic dispersion curve identification based on reciprocal condition number
CN115236755B (zh) 基于张量特征值的航磁异常边界检测方法、装置
CN113886754B (zh) 基于张量特征值的Theta Map法航磁边界检测方法、装置
CN112773396A (zh) 基于全波形反演的医学成像方法、计算机设备及存储介质
CN112731520A (zh) 一种基于结构张量扩散滤波的全波形反演方法及系统
Zhang et al. First break of the seismic signals in oil exploration based on information theory
CN115826038A (zh) 一种vti介质伴随状态法走时多参数层析成像方法及系统
CN112666615B (zh) 基于Tilt法获取精确航空磁测数据的方法、装置
Fonte-Boa et al. Interpolation artifacts as a result of spatial aliasing: A case study of the airborne magnetic data set of southeastern Minas Gerais, Brazil
CN109242770B (zh) 图像导引的地震速度插值方法及计算机可读存储介质
WO2017082912A1 (en) Two-dimensional imaging with multi-stage processing

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