CN113886754B - 基于张量特征值的Theta Map法航磁边界检测方法、装置 - Google Patents

基于张量特征值的Theta Map法航磁边界检测方法、装置 Download PDF

Info

Publication number
CN113886754B
CN113886754B CN202111177603.5A CN202111177603A CN113886754B CN 113886754 B CN113886754 B CN 113886754B CN 202111177603 A CN202111177603 A CN 202111177603A CN 113886754 B CN113886754 B CN 113886754B
Authority
CN
China
Prior art keywords
gradient
boundary detection
data matrix
total
tensor data
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
Application number
CN202111177603.5A
Other languages
English (en)
Other versions
CN113886754A (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 CN202111177603.5A priority Critical patent/CN113886754B/zh
Publication of CN113886754A publication Critical patent/CN113886754A/zh
Application granted granted Critical
Publication of CN113886754B publication Critical patent/CN113886754B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于张量特征值的Theta Map法航磁边界检测方法、装置及存储介质，方法包括：获取全磁力梯度张量数据矩阵M；根据M建立边界检测函数E以及深度分辨力增益因子Mz；分别计算E在x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez；根据Ex、Ey和Ez计算总水平梯度THDR以及总梯度模ASM；根据THDR、ASM、Mz、M，实现航磁边界检测。本发明提供的基于张量特征值的Theta Map法航磁边界检测方法、装置及存储介质可以解决Theta Map法存在“解析奇点”问题，提高计算稳定性，消除虚假干扰信息，增强信噪比，以及提高航磁数据处理转换的质量和地质体边界识别效果。

Description

基于张量特征值的Theta Map法航磁边界检测方法、装置
技术领域
本发明是关于航空磁测技术领域，特别是关于一种基于张量特征值的Theta Map法航磁边界检测方法、装置。
背景技术
航空磁力测量是将航空磁力仪(例如光泵式、核旋式和磁通门式)系统安装在飞行器，中通过观测地磁场参数(例如地磁场总强度T或总磁场异常△T或其梯度)，来寻找磁性或与磁性有关的矿体，以了解地质构造、进行磁性填图、解决城市和工程稳定性和考古等问题。
航空磁力测量数据是不同深度、不同形态、不同规模的磁性地质体磁场信息在观测面上的综合反映。但是，由于测量数据的误差或磁场的叠加，使得测量数据难以区分，给地质解释工作带来了难度。
工程技术与磁梯度张量探测仪器研发技术的不断地发展与成熟，应用磁力张量数据在分析、处理上述问题也得到了相应的发展。磁力张量数据是磁场矢量分量的梯度，包含了的磁场信息，能够反映目标体的矢量磁矩信息，张量数据具有高精度、高分辨率、多参量的优点，可用于描述场源体的磁化方向和几何形态，提高对目标地质体的分辨率。
目前，可以通过构造边界识别滤波器Theta Map的方法，提高对航磁异常的分辨能力。
基于此，本申请的发明人发现，Theta Map现有的Theta Map法在实际应用中存在一定程度上受磁异常分量与磁化方向影响、分辨率低和易产生虚假异常边界，且当分母
Figure BDA0003295911390000021
或接近0时，Theta Map法存在“解析奇点”，会使得计算结果不稳定，影响实际应用效果。
公开于该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域一般技术人员所公知的现有技术。
发明内容
本发明的目的在于提供一种基于张量特征值的Theta Map法航磁边界检测方法、装置，其能够解决Theta Map法存在“解析奇点”问题，提高计算稳定性，消除虚假干扰信息。
为实现上述目的，本发明提供了一种基于张量特征值的Theta Map法航磁边界检测方法，包括：
获取全磁力梯度张量数据矩阵M，其中，所述全磁力梯度张量数据矩阵M包括三维直角坐标系下，x,y,z三个方向的磁场分量分别在x,y,z方向的9个一阶梯度分量；
根据所述全磁力梯度张量数据矩阵M建立边界检测函数E以及深度分辨力增益因子Mz；
分别计算E在x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez；
根据所述x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez计算总水平梯度THDR以及总梯度模ASM；
根据总水平梯度THDR、总梯度模ASM、深度分辨力增益因子Mz、全磁力梯度张量数据矩阵M，实现航磁边界检测。
在本发明的一实施方式中，所述根据总水平梯度THDR、总梯度模ASM、深度分辨力增益因子Mz、全磁力梯度张量数据矩阵M，实现航磁边界检测包括：
根据以下公式实现航磁边界检测，所述公式包括：
Figure BDA0003295911390000031
其中，
Figure BDA0003295911390000032
Figure BDA0003295911390000033
δ是取值为0-1的调节系数，max(|E|)为边界检测函数E的最大值。
在本发明的一实施方式中，获取全磁力梯度张量数据矩阵M，包括：
获取实测的航磁数据；
根据三维直角坐标系确定航磁数据在x,y,z三个方向的磁场分量以及每个磁场分量分别在x,y,z方向的一阶梯度分量，9个所述一阶梯度分量构成全张量磁梯度数据。
在本发明的一实施方式中，所述
Figure BDA0003295911390000034
在本发明的一实施方式中，所述根据所述全磁力梯度张量数据矩阵M建立边界检测函数E包括：
根据所述全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3；
根据所述全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；
根据所述全磁力梯度张量数据矩阵的特征值λ1、λ2、λ3和所述全磁力梯度张量数据矩阵的总模值A，建立边界检测函数E，其中，E＝λ1·λ2·λ3·A。
在本发明的一实施方式中，所述全磁力梯度张量总模值A包含了全部9个张量元素的信息，其最大值对应着地质体的边界，
Figure BDA0003295911390000041
在本发明的一实施方式中，所述根据所述全磁力梯度张量数据矩阵M建立深度分辨力增益因子Mz包括：
根据以下公式，计算深度分辨力增益因子Mz，所述公式为：
Figure BDA0003295911390000042
在本发明的一实施方式中，根据以下公式计算均衡深浅部异常的调节系数δ，所述计算公式为：
Figure BDA0003295911390000043
其中，max(Mz)表示Mz最大值；min(Mz)表示Mz最小值。
为实现上述目的，本发明提供了一种基于张量特征值的Theta Map法航磁边界检测的装置，包括：
获取模块，用于获取全磁力梯度张量数据矩阵M，其中，所述全磁力梯度张量数据矩阵M包括三维直角坐标系下，x,y,z三个方向的磁场分量分别在x,y,z方向的9个一阶梯度分量；
建立模块，用于根据所述全磁力梯度张量数据矩阵M建立边界检测函数E以及深度分辨力增益因子Mz；
梯度计算模块，用于分别计算E在x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez；
总水平梯度THDR以及梯度模ASM计算模块，用于根据x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez计算总水平梯度THDR以及总梯度模ASM；
检测模块，用于根据总水平梯度THDR、总梯度模ASM、深度分辨力增益因子Mz、全磁力梯度张量数据矩阵M，实现航磁边界检测。
为实现上述目的，本发明还提供了一种存储介质，所述存储介质存储有计算机可执行指令，所述计算机可执行指令用于执行上述的基于张量特征值的Theta Map法航磁边界检测方法。
与现有技术相比，根据本发明的基于张量特征值的Theta Map法航磁边界检测方法、装置及存储介质，能更好的探测出埋深不同的多源场物体的边界，使边界识别结果更收敛，且该方法有效地避免了磁化方向和噪声对结果的干扰，解决了现有Theta Map法存在“解析奇点”问题，提高计算稳定性和消除产生虚假航磁数据地质体边界干扰、增强信噪比，以及提高了不同埋深的地质体的边界位置增强和提取能力，具有更高的分辨率和精度。
附图说明
图1是根据本发明一实施方式的基于张量特征值的Theta Map法航磁边界检测方法的流程图；
图2是根据本发明一实施方式的基于张量特征值的Theta Map法航磁边界检测装置的结构示意图。
主要附图标记说明：
1-获取模块，2-建立模块，3-梯度计算模块，4-总水平梯度THDR以及总梯度模ASM计算模块，5-检测模块。
具体实施方式
下面结合附图，对本发明的具体实施方式进行详细描述，但应当理解本发明的保护范围并不受具体实施方式的限制。
除非另有其它明确表示，否则在整个说明书和权利要求书中，术语“包括”或其变换如“包含”或“包括有”等等将被理解为包括所陈述的元件或组成部分，而并未排除其它元件或其它组成部分。
本发明实施例提供一种基于张量特征值的Theta Map法航磁边界检测方法的，参见图1，其为基于张量特征值的Theta Map法航磁边界检测方法的流程图，包括：步骤S1-步骤S4。
步骤1，获取全磁力梯度张量数据矩阵M，其中，所述全磁力梯度张量数据矩阵M包括三维直角坐标系下，x,y,z三个方向的磁场分量分别在x,y,z方向的9个一阶梯度分量。
在一种实现方式中，全磁力梯度张量数据矩阵M可以通过实测全磁力梯度张量数据来获取，或者，可以通过以下步骤获取：
获取实测的航磁数据；
根据三维直角坐标系确定磁场数据在x,y,z三个方向的磁场分量以及每个磁场分量分别在x,y,z方向的一阶梯度分量，9个所述一阶梯度分量构成全张量磁梯度数据。具体的，航磁数据为地下磁性地质体引起的磁异常信号。
其中，
Figure BDA0003295911390000071
步骤2，根据所述全磁力梯度张量数据矩阵M建立边界检测函数E以及深度分辨力增益因子Mz。
在一种实现方式中，步骤2中的根据所述全磁力梯度张量数据矩阵M建立边界检测函数E可以包括：
步骤201，根据所述全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3，其中，所述全磁力梯度张量数据矩阵的特征值λ1、λ2、λ3均对应着航磁数据目标地质体的边界。
步骤202，根据所述全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A。
具体的，全磁力梯度张量总模值A包含了全部9个张量元素的信息，其最大值对应着地质体的边界，其中，
Figure BDA0003295911390000072
步骤203，根据所述全磁力梯度张量数据矩阵的特征值λ1、λ2、λ3和所述全磁力梯度张量数据矩阵的总模值A，建立边界检测函数E。
E同时具备了张量特征值和总模值的属性，能够提高对浅部目标地质体的识别精度，但检测深部目标地质体的分辨率低，所以需要进一步提高检测深部目标地质体的边界的能力。
其中，E＝λ1·λ2·λ3·A。
在一种实现方式中，步骤2中的根据所述全磁力梯度张量数据矩阵M建立深度分辨力增益因子Mz，可以包括：
根据以下公式计算深度分辨力增益因子Mz，所述公式为：
Figure BDA0003295911390000081
由此，可以提高垂向探测能力，即提高检测深部目标地质体的边界的能力。
步骤3，分别计算E在x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez，其中
Figure BDA0003295911390000082
步骤4，根据所述x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez计算总水平梯度THDR以及总梯度模ASM；
其中，
Figure BDA0003295911390000083
Figure BDA0003295911390000084
步骤5，根据总水平梯度THDR、总梯度模ASM、深度分辨力增益因子Mz、全磁力梯度张量数据矩阵M，实现航磁边界检测。
所述根据总水平梯度THDR、总梯度模ASM、深度分辨力增益因子Mz、全磁力梯度张量数据矩阵M，实现航磁边界检测包括：
根据以下的改进后的Theta Map(iTheta)公式获取精确航空磁测数据，所述公式包括：
Figure BDA0003295911390000085
即：
Figure BDA0003295911390000091
其中，
Figure BDA0003295911390000092
Figure BDA0003295911390000093
δ是均衡深浅部异常的调节系数，一般取值为0-1，max(|E|)为边界检测函数E的最大值。
在一种实现方式中，可以根据以下公式计算均衡深浅部异常的调节系数δ，所述计算公式为：
Figure BDA0003295911390000094
max(Mz)表示Mz最大值；min(Mz)表示Mz最小值。
本实施例提供与现有的Theta Map表达式对比：
Figure BDA0003295911390000095
其中，T为航空磁测数据；
Figure BDA0003295911390000096
分别是T在x、y、z三个方向一阶导数(或一阶梯度)。航空磁测数据T即航空磁力异常场是地壳中的含铁磁性地质体在地磁场作用下所产生的附加磁场。
在根据上述步骤得构建Theta Map滤波器之后，可以获取精确航空磁测数据，进一步准确推断构造体场源的边界、深度、产状、规模及场分布规律和物理性质等，对划分大地构造单元、进行构造分区、确定断裂构造带的位、区分不同岩性与地层的分布及进行物性填图等问题有重要意义。
通过本实施例提供的基于张量特征值的Theta Map法航磁边界检测方法，新构建了一个合理的均衡航磁数据目标地质体边界检测方法，能更好的探测出埋深不同的多源场物体的边界，使边界识别结果更收敛，且该方法有效地避免了磁化方向和噪声对结果的干扰，解决了现有Theta Map法存在“解析奇点”问题，提高计算稳定性和消除产生虚假航磁数据地质体边界干扰、增强信噪比，以及提高了不同埋深的地质体的边界位置增强和提取能力，具有更高的分辨率和精度。
本发明实施例还提供了一种基于Theta Map法获取精确航空磁测数据的装置，请参阅图2，其为基于张量特征值的Theta Map法航磁边界检测装置的结构示意图，包括：获取模块1、建立模块2、梯度计算模块3，总水平梯度THDR以及总梯度模ASM计算模块4，以及检测模块5。
获取模块1用于获取全磁力梯度张量数据矩阵M，其中，所述全磁力梯度张量数据矩阵M包括三维直角坐标系下，x,y,z三个方向的磁场分量分别在x,y,z方向的9个一阶梯度分量。
建立模块2用于根据所述全磁力梯度张量数据矩阵M建立边界检测函数E以及深度分辨力增益因子Mz。
梯度计算模块3用于分别计算E在x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez。
总水平梯度THDR以及总梯度模ASM计算模块4用于根根据所述x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez计算总水平梯度THDR以及总梯度模ASM。
检测模块5用于根据总水平梯度THDR、总梯度模ASM、深度分辨力增益因子Mz、全磁力梯度张量数据矩阵M，实现航磁边界检测。
本发明实施例还提供了一种存储介质，所述存储介质存储有计算机可执行指令，其包含用于执行上述基于张量特征值的Theta Map法航磁边界检测方法的程序，该计算机可执行指令可执行上述任意方法实施例中的方法。
其中，所述存储介质可以是计算机能够存取的任何可用介质或数据存储设备，包括但不限于磁性存储器(例如软盘、硬盘、磁带、磁光盘(MO)等)、光学存储器(例如CD、DVD、BD、HVD等)、以及半导体存储器(例如ROM、EPROM、EEPROM、非易失性存储器(NAND FLASH)、固态硬盘(SSD))等。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
前述对本发明的具体示例性实施方案的描述是为了说明和例证的目的。这些描述并非想将本发明限定为所公开的精确形式，并且很显然，根据上述教导，可以进行很多改变和变化。对示例性实施例进行选择和描述的目的在于解释本发明的特定原理及其实际应用，从而使得本领域的技术人员能够实现并利用本发明的各种不同的示例性实施方案以及各种不同的选择和改变。本发明的范围意在由权利要求书及其等同形式所限定。

Claims (7)

1.一种基于张量特征值的Theta Map法航磁边界检测方法，其特征在于，包括：
获取全磁力梯度张量数据矩阵M，其中，所述全磁力梯度张量数据矩阵M包括三维直角坐标系下，x,y,z三个方向的磁场分量分别在x,y,z方向的9个一阶梯度分量；
根据所述全磁力梯度张量数据矩阵M建立边界检测函数E以及深度分辨力增益因子Mz；
分别计算E在x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez；
根据所述x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez计算总水平梯度THDR以及总梯度模ASM；
根据总水平梯度THDR、总梯度模ASM、深度分辨力增益因子Mz、全磁力梯度张量数据矩阵M，实现航磁边界检测；
所述根据总水平梯度THDR、总梯度模ASM、深度分辨力增益因子Mz、全磁力梯度张量数据矩阵M，实现航磁边界检测包括：
根据以下公式实现航磁边界检测，所述公式包括：
Figure FDA0003520199400000011
其中，
Figure FDA0003520199400000012
Figure FDA0003520199400000013
δ是取值为0-1的调节系数，max(|E|)为边界检测函数E的最大值；
所述根据所述全磁力梯度张量数据矩阵M建立边界检测函数E包括：
根据所述全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3；
根据所述全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；
根据所述全磁力梯度张量数据矩阵的特征值λ1、λ2、λ3和所述全磁力梯度张量数据矩阵的总模值A，建立边界检测函数E，其中，E＝λ1·λ2·λ3·A；
所述根据所述全磁力梯度张量数据矩阵M建立深度分辨力增益因子Mz包括：
根据以下公式，计算深度分辨力增益因子Mz，公式为：
Figure FDA0003520199400000021
2.根据权利要求1所述的基于张量特征值的Theta Map法航磁边界检测方法，其特征在于，获取全磁力梯度张量数据矩阵M，包括：
获取实测的航磁数据；
根据三维直角坐标系确定航磁数据在x,y,z三个方向的磁场分量以及每个磁场分量分别在x,y,z方向的一阶梯度分量，9个所述一阶梯度分量构成全张量磁梯度数据。
3.根据权利要求1或2所述的基于张量特征值的Theta Map法航磁边界检测方法，其特征在于，所述
Figure FDA0003520199400000031
4.根据权利要求2所述的基于张量特征值的Theta Map法航磁边界检测方法，其特征在于，所述全磁力梯度张量数据矩阵的总模值A包含了全部9个张量元素的信息，其最大值对应着地质体的边界，
Figure FDA0003520199400000032
5.根据权利要求2所述的基于张量特征值的Theta Map法航磁边界检测方法，其特征在于，
根据以下公式计算均衡深浅部异常的调节系数δ，计算公式为：
Figure FDA0003520199400000033
其中，max(Mz)表示Mz最大值；min(Mz)表示Mz最小值。
6.一种基于张量特征值的Theta Map法航磁边界检测的装置，其特征在于，包括：
获取模块，用于获取全磁力梯度张量数据矩阵M，其中，所述全磁力梯度张量数据矩阵M包括三维直角坐标系下，x,y,z三个方向的磁场分量分别在x,y,z方向的9个一阶梯度分量；
建立模块，用于根据所述全磁力梯度张量数据矩阵M建立边界检测函数E以及深度分辨力增益因子Mz；
梯度计算模块，用于分别计算E在x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez；
总水平梯度THDR以及梯度模ASM计算模块，用于根据x方向水平梯度Ex、y方向水平梯度Ey和z方向垂向梯度Ez计算总水平梯度THDR以及总梯度模ASM；
检测模块，用于根据总水平梯度THDR、总梯度模ASM、深度分辨力增益因子Mz、全磁力梯度张量数据矩阵M，实现航磁边界检测；
所述根据总水平梯度THDR、总梯度模ASM、深度分辨力增益因子Mz、全磁力梯度张量数据矩阵M，实现航磁边界检测包括：
根据以下公式实现航磁边界检测，所述公式包括：
Figure FDA0003520199400000041
其中，
Figure FDA0003520199400000042
Figure FDA0003520199400000043
δ是取值为0-1的调节系数，max(|E|)为边界检测函数E的最大值；
所述根据所述全磁力梯度张量数据矩阵M建立边界检测函数E包括：
根据所述全磁力梯度张量数据矩阵M计算得到该矩阵的三个特征值λ1、λ2、λ3；
根据所述全磁力梯度张量数据矩阵M计算得到该矩阵的总模值A；
根据所述全磁力梯度张量数据矩阵的特征值λ1、λ2、λ3和所述全磁力梯度张量数据矩阵的总模值A，建立边界检测函数E，其中，E＝λ1·λ2·λ3·A；
所述根据所述全磁力梯度张量数据矩阵M建立深度分辨力增益因子Mz包括：
根据以下公式，计算深度分辨力增益因子Mz，公式为：
Figure FDA0003520199400000051
7.一种存储介质，其特征在于，所述存储介质存储有计算机可执行指令，所述计算机可执行指令用于执行权利要求1-5任意一项所述的基于张量特征值的Theta Map法航磁边界检测方法。
CN202111177603.5A 2021-10-09 2021-10-09 基于张量特征值的Theta Map法航磁边界检测方法、装置 Active CN113886754B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111177603.5A CN113886754B (zh) 2021-10-09 2021-10-09 基于张量特征值的Theta Map法航磁边界检测方法、装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111177603.5A CN113886754B (zh) 2021-10-09 2021-10-09 基于张量特征值的Theta Map法航磁边界检测方法、装置

Publications (2)

Publication Number Publication Date
CN113886754A CN113886754A (zh) 2022-01-04
CN113886754B true CN113886754B (zh) 2022-04-15

Family

ID=79005649

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111177603.5A Active CN113886754B (zh) 2021-10-09 2021-10-09 基于张量特征值的Theta Map法航磁边界检测方法、装置

Country Status (1)

Country Link
CN (1) CN113886754B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115236755B (zh) * 2022-07-25 2023-10-03 中国自然资源航空物探遥感中心 基于张量特征值的航磁异常边界检测方法、装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108508490A (zh) * 2018-03-07 2018-09-07 吉林大学 一种基于解析信号的磁张量梯度数据均衡边界识别方法
CN111007571A (zh) * 2019-11-28 2020-04-14 吉林大学 一种基于三维构造张量的航磁数据地质体边界识别方法
CN112749493A (zh) * 2021-01-25 2021-05-04 中国自然资源航空物探遥感中心 基于全磁力梯度张量特征值的地质体边界检测方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110133749B (zh) * 2019-05-30 2021-06-22 中国地质科学院 一种地质资源勘探中重磁数据处理方法及其系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108508490A (zh) * 2018-03-07 2018-09-07 吉林大学 一种基于解析信号的磁张量梯度数据均衡边界识别方法
CN111007571A (zh) * 2019-11-28 2020-04-14 吉林大学 一种基于三维构造张量的航磁数据地质体边界识别方法
CN112749493A (zh) * 2021-01-25 2021-05-04 中国自然资源航空物探遥感中心 基于全磁力梯度张量特征值的地质体边界检测方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于磁力梯度全张量特征值的均衡边界识别方法;郑强 等;《石油地球物理勘探》;20200430;第55卷(第4期);第454-464页 *
重力张量数据的目标体边缘检测方法探讨;朱自强等;《石油地球物理勘探》;20110615;第46卷(第03期);第482-488页 *

Also Published As

Publication number Publication date
CN113886754A (zh) 2022-01-04

Similar Documents

Publication Publication Date Title
CN113886753B (zh) 基于张量特征值的Tilt法航磁边界检测方法、装置
CN110095738B (zh) 磁梯度仪参数设计方法
CN107272069B (zh) 基于磁异常梯度的磁性目标追踪方法
CN103153790B (zh) 使用运动传感器和附接至装置的磁力计的测量数据估计该装置在重力参照系中的偏航角的设备和方法
Dannemiller et al. A new method for determination of magnetization direction
Gang et al. Detection of ferromagnetic target based on mobile magnetic gradient tensor system
CN110187393B (zh) 一种基于广义回归神经网络的航磁补偿方法
CN112749493A (zh) 基于全磁力梯度张量特征值的地质体边界检测方法及系统
CN110133716A (zh) 基于组合模型加权函数的磁异常数据三维反演方法
CN114488326B (zh) 一种提高探测航磁数据地质体边界能力的方法及系统
CN113886754B (zh) 基于张量特征值的Theta Map法航磁边界检测方法、装置
Yin et al. Magnetic dipole localization and magnetic moment estimation method based on normalized source strength
Schneider et al. Inversion of geo-magnetic full-tensor gradiometer data
Ma et al. Interpretation of magnetic anomalies by horizontal and vertical derivatives of the analytic signal
CN115236755B (zh) 基于张量特征值的航磁异常边界检测方法、装置
Li et al. Application of Helbig integrals to magnetic gradient tensor multi-target detection
CN113917544A (zh) 基于磁梯度张量特征值的近地表目标位置快速圈定方法
Liu et al. The IDQ curve: A tool for evaluating the direction of remanent magnetization from magnetic anomalies
CN112666619B (zh) 基于标准偏差法取精确航空磁测数据的方法、装置
CN112666615B (zh) 基于Tilt法获取精确航空磁测数据的方法、装置
Zhang et al. A component compensation method for magnetic interferential field
Wang et al. The stability optimization algorithm of second-order magnetic gradient tensor
Li et al. Preprocessed method and application of magnetic gradient tensor data
Guo et al. Three-directional analytic signal analysis and interpretation of magnetic gradient tensor
Mendonça Inversion of gravity-field inclination to map the basement relief of sedimentary basins

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