CN112666615B - 基于Tilt法获取精确航空磁测数据的方法、装置 - Google Patents

基于Tilt法获取精确航空磁测数据的方法、装置 Download PDF

Info

Publication number
CN112666615B
CN112666615B CN202011497559.1A CN202011497559A CN112666615B CN 112666615 B CN112666615 B CN 112666615B CN 202011497559 A CN202011497559 A CN 202011497559A CN 112666615 B CN112666615 B CN 112666615B
Authority
CN
China
Prior art keywords
gradient
data
horizontal gradient
wmean
aeromagnetic
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
CN202011497559.1A
Other languages
English (en)
Other versions
CN112666615A (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 CN202011497559.1A priority Critical patent/CN112666615B/zh
Publication of CN112666615A publication Critical patent/CN112666615A/zh
Application granted granted Critical
Publication of CN112666615B publication Critical patent/CN112666615B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Magnetic Variables (AREA)

Abstract

本发明公开了一种基于Tilt法获取精确航空磁测数据的方法、装置，方法包括：根据航空磁测数据T，计算在x方向水平梯度
Figure DDA0002842646620000011
在y方向水平梯度
Figure DDA0002842646620000012
和在z方向垂向梯度
Figure DDA0002842646620000013
根据x方向水平梯度
Figure DDA0002842646620000014
以及在y方向水平梯度
Figure DDA0002842646620000015
计算总水平梯度THDR，根据所述z方向垂向梯度
Figure DDA0002842646620000016
和总水平梯度THDR计算z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)根据z方向垂向梯度平均值VDR、总水平梯度THDR、z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)，获取精确航空磁测数据。本发明提供的获取精确航空磁测数据的方法、装置及存储介质，可以解决Tilt法存在“解析奇点”问题，提高计算稳定性，消除产生虚假航磁数据地质体的边界干扰，增强信噪比，以及提高航磁数据处理转换的质量和地质体边界识别效果。

Description

基于Tilt法获取精确航空磁测数据的方法、装置
技术领域
本发明是涉及航空磁测技术领域，特别是涉及一种基于Tilt法获取精确航空磁测数据的方法、装置。
背景技术
航空磁力测量是将航空磁力仪(例如光泵式、核旋式和磁通门式)系统安装在飞行器，中通过观测地磁场参数(例如地磁场总强度T或总磁场异常△T或其梯度)，来寻找磁性或与磁性有关的矿体，以了解地质构造、进行磁性填图、解决城市和工程稳定性和考古等问题。
航空磁力测量数据是不同深度、不同形态、不同规模的磁性地质体磁场信息在观测面上的综合反映。但是，由于测量数据的误差或磁场的叠加，使得测量数据难以区分，给地质解释工作带来了难度。
目前，可以通过构造边界识别滤波器Tilt的方法，提高对航磁异常的分辨能力。Tilt法是通过一阶导数的比值平衡高幅值和低幅值异常达到航磁数据边界增强。Tilt法表达式：
Figure BDA0002842646600000011
其中，
Figure BDA0002842646600000012
T为实测的航空磁测数据，x、y为空间坐标的两个方向。
基于此，本申请的发明人发现，Tilt法在实际应用中存在横向分辨力低、精度差和易产生虚假异常边界，且当分母
Figure BDA0002842646600000013
或接近0时，Tilt法存在“解析奇点”，会使得计算结果不稳定，影响实际应用效果。
公开于该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域一般技术人员所公知的现有技术。
发明内容
为解决上述问题，本发明实施例的目的在于提供一种基于Tilt法获取精确航空磁测数据的方法、装置。
为实现上述目的，本发明提供了一种基于Tilt法获取精确航空磁测数据的方法，包括：根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000021
所述航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000022
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002842646600000023
根据所述航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000024
以及所述航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000025
计算总水平梯度THDR；根据所述z方向垂向梯度
Figure BDA0002842646600000026
和总水平梯度THDR计算z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)；根据z方向垂向梯度
Figure BDA0002842646600000027
总水平梯度THDR、z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)，获取精确航空磁测数据。
在一优选的实施方式中，所述根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000028
所述航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000029
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002842646600000031
之前，包括：获取实测的航空磁测数据T，并进行预处理。
在一优选的实施方式中，所述预处理包括坐标转换、正常场校正、日变校正、滞后校正以及磁场水平调整。
在一优选的实施方式中，所述根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000032
所述航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000033
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002842646600000034
包括：根据预处理后的航空磁测数据T进行频率域内二维傅里叶变换；根据二维傅里叶变换与希尔伯特变换之间的关系进行计算，得到频率域希尔伯特变换；根据所述频率域希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000035
预处理后的航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000036
和z方向垂向梯度
Figure BDA0002842646600000037
在一优选的实施方式中，所述根据所述希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000038
预处理后的航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000039
和z方向垂向梯度
Figure BDA00028426466000000310
包括：二维希尔伯特变换在x方向上的分量对应为
Figure BDA00028426466000000311
二维希尔伯特变换在y方向上的分量对应为
Figure BDA00028426466000000312
预处理后的T对应为
Figure BDA00028426466000000313
在一优选的实施方式中，根据所述在z方向垂向梯度
Figure BDA00028426466000000314
和总水平梯度THDR计算z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)包括：提供大小为m×n的滑动窗口；在所述滑动窗口内计算窗口内所有元素在z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)。
在一优选的实施方式中，所述根据z方向垂向梯度
Figure BDA0002842646600000041
总水平梯度THDR、z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)，获取精确航空磁测数据包括：
根据以下公式获取精确航空磁测数据，所述公式包括：
Figure BDA0002842646600000042
其中，
Figure BDA0002842646600000043
Wmean()表示滑动窗口均值，T为实测的航空磁测数据，x、y为空间坐标的两个方向，Δ为正则化调节因子。
在一优选的实施方式中，所述正则化调节因子为0.001。
为实现上述目的，本发明还提供了一种基于Tilt法获取精确航空磁测数据的装置，包括：梯度计算模块，用于根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000044
所述航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000045
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002842646600000046
总水平梯度计算模块，用于根据所述航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000047
以及所述航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000051
计算总水平梯度THDR；均值计算模块，用于根据所述z方向垂向梯度
Figure BDA0002842646600000052
和总水平梯度THDR计算z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)；精确航空磁测数据获取模块，用于根据z方向垂向梯度
Figure BDA0002842646600000053
总水平梯度THDR、z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)，获取精确航空磁测数据。
为实现上述目的，本发明还提供了一种存储介质，所述存储介质存储有计算机可执行指令，所述计算机可执行指令用于执行上述任意一项所述的基于Tilt法获取精确航空磁测数据的方法。
与现有技术相比，根据本发明的基于Tilt法的获取精确航空磁测数据的方法、装置及存储介质，可以解决Tilt法存在“解析奇点”问题，提高计算稳定性，消除产生虚假航磁数据地质体的边界干扰，增强信噪比，以及提高航磁数据处理转换的质量和地质体边界识别效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出了本发明实施例所提供的基于Tilt法获取精确航空磁测数据的方法的流程图；
图2示出了本发明实施例所提供的基于Tilt法获取精确航空磁测数据的装置的结构示意图。
具体实施方式
下面结合附图，对本发明的具体实施方式进行详细描述，但应当理解本发明的保护范围并不受具体实施方式的限制。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
除非另有其它明确表示，否则在整个说明书和权利要求书中，术语“包括”或其变换如“包含”或“包括有”等等将被理解为包括所陈述的元件或组成部分，而并未排除其它元件或其它组成部分。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
本发明实施例提供一种基于Tilt法获取精确航空磁测数据的方法的流程图，参见图1，包括：步骤S1-步骤S4。
步骤1，根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000061
所述航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000062
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002842646600000063
航空磁测数据T即航空磁力异常场是地壳中的含铁磁性地质体在地磁场作用下所产生的附加磁场。
在一种实现方式中，步骤S1之前，还可以包括：获取实测的航空磁测数据T，并进行预处理。所述预处理包括坐标转换、正常场校正、日变校正、滞后校正以及磁场水平调整中的一种或多种。通过预处理，可以消除在航空磁场测量数据中，因多种因素会造成数据存在噪声和误差。
在一种实现方式中，步骤S1进一步可以包括：
步骤S101，根据预处理后的航空磁测数据T进行频率域内二维傅里叶变换；
步骤S102，根据二维傅里叶变换与希尔伯特变换之间的关系进行计算，得到频率域希尔伯特变换；
步骤S103，根据所述频率域希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000071
预处理后的航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000072
和z方向垂向梯度
Figure BDA0002842646600000073
具体的，根据所述希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000074
预处理后的航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000075
和z方向垂向梯度
Figure BDA0002842646600000076
包括：二维希尔伯特变换在x方向上的分量对应为
Figure BDA0002842646600000077
二维希尔伯特变换在y方向上的分量对应为
Figure BDA0002842646600000078
预处理后的T对应为
Figure BDA0002842646600000079
二维希尔伯特变换，其傅里叶域变换因子为：
DH(u,v)＝-i·sign(u,v)，
Figure BDA00028426466000000710
其中，i2＝-1，u、v是x和y方向上的圆波数。希尔伯特变换在x和y方向上的分量可以表示为：
Figure BDA0002842646600000081
其中，F[]、F-1[]分别表示傅里叶正变换和逆变换；Hx与Hy分别代表二维希尔伯特变换在x和y方向上的分量。
因此，x方向水平梯度
Figure BDA0002842646600000082
y方向水平梯度
Figure BDA0002842646600000083
和z方向垂向梯度
Figure BDA0002842646600000084
可以利用希尔伯特变换计算得到，二者具有等价关系：
Figure BDA0002842646600000085
应用上述等价关系优点在于直接希尔伯特变换不会放大航空磁测数据中的噪声干扰。
步骤2，根据所述航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000086
以及所述航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000087
计算总水平梯度THDR。
具体的，可以通过以下公式总水平梯度THDR，所述公式包括：
Figure BDA0002842646600000088
其中，
Figure BDA0002842646600000091
为航空磁测数据T在x方向水平梯度，
Figure BDA0002842646600000092
为航空磁测数据T在在y方向水平梯度。
步骤3，根据所述z方向垂向梯度
Figure BDA0002842646600000093
和总水平梯度THDR计算z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)。
具体的，步骤3可以包括：提供大小为m×n的滑动窗口；在所述滑动窗口内计算窗口内所有元素在z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)。
具体的，可以通过以下公式进行计算：
Figure BDA0002842646600000094
Figure BDA0002842646600000095
其中，i表示m×n大小的滑动窗口内的第几个元素。
步骤4，根据z方向垂向梯度
Figure BDA0002842646600000096
总水平梯度THDR、z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)，获取精确航空磁测数据。
步骤4可以根据以下公式获取精确航空磁测数据，所述公式包括：
Figure BDA0002842646600000097
其中，
Figure BDA0002842646600000098
Wmean()表示滑动窗口均值，T为预处理后的航空磁测数据，x、y为空间坐标的两个方向，Δ为正则化调节因子，Δ的取值范围为0-1。
不同的正则化调节因子的取值对边界识别结果有影响。正则化调节因子是通过深度学习样本库的自动获取方法获得的，得到Δ的取值范围为0-1。Δ取值过小会产生虚假多余边界，Δ取值过大会使结果逐渐趋近与总水平导数其受不同地质条件产生变化，根据地质解释结果的优劣确定最佳参数值。作为最优的一种实现方式，正则化调节因子的取值可以为0.001。由此，可以提高边界的识别效果。
在根据z方向垂向梯度
Figure BDA0002842646600000101
总水平梯度THDR、z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)构建Tilt滤波器之后，可以获取精确航空磁测数据，进一步准确推断构造体场源的边界、深度、产状、规模及场分布规律和物理性质等，对划分大地构造单元、进行构造分区、确定断裂构造带的位、区分不同岩性与地层的分布及进行物性填图等问题有重要意义。
由此，解决Tilt法存在“解析奇点”问题，提高计算稳定性，消除产生虚假航磁数据地质体的边界干扰，增强信噪比，以及提高航磁数据处理转换的质量和地质体边界识别效果。
通过本实施例提供的方法，可以使航磁异常输出振幅趋于均衡化，可以增强深层源弱异常，准确圈定较小较窄的异常体边界，为后期资料处理解释提供丰富的地质体信息；其对浅源和深源响应均等，并且能处理同一水平上幅度动态范围大的场源，由反正切函数的特性所决定，不管总水平梯度的绝对值或垂直梯度的振幅如何变化，始终使计算结果控制
Figure BDA0002842646600000102
在之间，有效的平衡了高幅值异常和低幅值异常；对于深部场源，即使它的垂向梯度和水平梯度都很小，两者的比值仍然会很大，零值处对应场源边界的位置，不受埋深的影响；因此，能更好的探测出埋深不同的多源场物体的边界。
本发明实施例还提供了一种基于Tilt法获取精确航空磁测数据的装置，包括：梯度计算模块1、总水平梯度计算模块2、均值计算模块3以及精确航空磁测数据获取模块4。
梯度计算模块1用于根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000111
所述航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000112
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002842646600000113
总水平梯度计算模块2用于根据所述航空磁测数据T在x方向水平梯度
Figure BDA0002842646600000114
以及所述航空磁测数据T在y方向水平梯度
Figure BDA0002842646600000115
计算总水平梯度THDR。
均值计算模块3用于根据所述z方向垂向梯度
Figure BDA0002842646600000116
和总水平梯度THDR计算z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)。
精确航空磁测数据获取模块4用于根据z方向垂向梯度
Figure BDA0002842646600000117
总水平梯度THDR、z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)，获取精确航空磁测数据。
本发明实施例还提供了一种存储介质，所述存储介质存储有计算机可执行指令，其包含用于执行上述基于Tilt法获取精确航空磁测数据的方法的程序，该计算机可执行指令可执行上述任意方法实施例中的方法。
其中，所述存储介质可以是计算机能够存取的任何可用介质或数据存储设备，包括但不限于磁性存储器(例如软盘、硬盘、磁带、磁光盘(MO)等)、光学存储器(例如CD、DVD、BD、HVD等)、以及半导体存储器(例如ROM、EPROM、EEPROM、非易失性存储器(NAND FLASH)、固态硬盘(SSD))等。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

Claims (9)

1.一种基于Tilt法获取精确航空磁测数据的方法，其特征在于，包括：
根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure FDA0003172277420000011
所述航空磁测数据T在y方向水平梯度
Figure FDA0003172277420000012
和所述航空磁测数据T在z方向垂向梯度
Figure FDA0003172277420000013
根据所述航空磁测数据T在x方向水平梯度
Figure FDA0003172277420000014
以及所述航空磁测数据T在y方向水平梯度
Figure FDA0003172277420000015
计算总水平梯度THDR；
根据所述z方向垂向梯度
Figure FDA0003172277420000016
和所述总水平梯度THDR计算z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)；
根据所述z方向垂向梯度
Figure FDA0003172277420000017
所述总水平梯度THDR、所述z方向垂向梯度平均值Wmean(VDR)和所述总水平梯度平均值Wmean(THDR)，获取精确航空磁测数据；
所述根据z方向垂向梯度
Figure FDA0003172277420000018
总水平梯度THDR、z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)，获取精确航空磁测数据包括：
根据以下公式获取精确航空磁测数据，所述公式包括：
Figure FDA0003172277420000021
其中，
Figure FDA0003172277420000022
Wmean()表示滑动窗口均值，T为实测的航空磁测数据，x、y为空间坐标的两个方向，Δ为正则化调节因子。
2.根据权利要求1所述的方法，其特征在于，所述根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure FDA0003172277420000023
所述航空磁测数据T在y方向水平梯度
Figure FDA0003172277420000024
和所述航空磁测数据T在z方向垂向梯度
Figure FDA0003172277420000025
之前，包括：
获取实测的航空磁测数据T，并进行预处理。
3.根据权利要求2所述的方法，其特征在于，所述预处理包括坐标转换、正常场校正、日变校正、滞后校正以及磁场水平调整中的一种或多种。
4.根据权利要求2所述的方法，其特征在于，所述根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure FDA0003172277420000026
所述航空磁测数据T在y方向水平梯度
Figure FDA0003172277420000027
和所述航空磁测数据T在z方向垂向梯度
Figure FDA0003172277420000028
包括：
根据预处理后的航空磁测数据T进行频率域内二维傅里叶变换；
根据二维傅里叶变换与希尔伯特变换之间的关系进行计算，得到频率域希尔伯特变换；
根据所述频率域希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure FDA0003172277420000031
预处理后的航空磁测数据T在y方向水平梯度
Figure FDA0003172277420000032
和z方向垂向梯度
Figure FDA0003172277420000033
5.根据权利要求4所述的方法，其特征在于，所述根据所述希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure FDA0003172277420000034
预处理后的航空磁测数据T在y方向水平梯度
Figure FDA0003172277420000035
和预处理后的航空磁测数据T在z方向垂向梯度
Figure FDA0003172277420000036
包括：
二维希尔伯特变换在x方向上的分量对应为
Figure FDA0003172277420000037
二维希尔伯特变换在y方向上的分量对应为
Figure FDA0003172277420000038
二维希尔伯特变换在z方向上的分量对应为
Figure FDA0003172277420000039
6.根据权利要求4所述的方法，其特征在于，根据所述在z方向垂向梯度
Figure FDA00031722774200000310
和所述总水平梯度THDR计算z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)包括：
提供大小为m×n的滑动窗口；
在所述滑动窗口内计算窗口内所有元素在z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)。
7.根据权利要求1所述的方法，其特征在于，
所述正则化调节因子为0.001。
8.一种基于Tilt法获取精确航空磁测数据的装置，其特征在于，包括：
梯度计算模块，用于根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure FDA0003172277420000041
所述航空磁测数据T在y方向水平梯度
Figure FDA0003172277420000042
和所述航空磁测数据T在z方向垂向梯度
Figure FDA0003172277420000043
总水平梯度计算模块，用于根据所述航空磁测数据T在x方向水平梯度
Figure FDA0003172277420000044
以及所述航空磁测数据T在y方向水平梯度
Figure FDA0003172277420000045
计算总水平梯度THDR；
均值计算模块，用于根据所述z方向垂向梯度
Figure FDA0003172277420000046
和所述总水平梯度THDR计算z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)；
精确航空磁测数据获取模块，用于根据所述z方向垂向梯度
Figure FDA0003172277420000047
所述总水平梯度THDR、所述z方向垂向梯度平均值Wmean(VDR)和所述总水平梯度平均值Wmean(THDR)，获取精确航空磁测数据；
所述根据z方向垂向梯度
Figure FDA0003172277420000048
总水平梯度THDR、z方向垂向梯度平均值Wmean(VDR)和总水平梯度平均值Wmean(THDR)，获取精确航空磁测数据包括：
根据以下公式获取精确航空磁测数据，所述公式包括：
Figure FDA0003172277420000051
其中，
Figure FDA0003172277420000052
Wmean()表示滑动窗口均值，T为实测的航空磁测数据，x、y为空间坐标的两个方向，Δ为正则化调节因子。
9.一种存储介质，其特征在于，所述存储介质存储有计算机可执行指令，所述计算机可执行指令用于执行权利要求1-7任意一项所述的基于Tilt法获取精确航空磁测数据的方法。
CN202011497559.1A 2020-12-17 2020-12-17 基于Tilt法获取精确航空磁测数据的方法、装置 Active CN112666615B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011497559.1A CN112666615B (zh) 2020-12-17 2020-12-17 基于Tilt法获取精确航空磁测数据的方法、装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011497559.1A CN112666615B (zh) 2020-12-17 2020-12-17 基于Tilt法获取精确航空磁测数据的方法、装置

Publications (2)

Publication Number Publication Date
CN112666615A CN112666615A (zh) 2021-04-16
CN112666615B true CN112666615B (zh) 2021-09-14

Family

ID=75404885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011497559.1A Active CN112666615B (zh) 2020-12-17 2020-12-17 基于Tilt法获取精确航空磁测数据的方法、装置

Country Status (1)

Country Link
CN (1) CN112666615B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113886753B (zh) * 2021-10-09 2022-04-12 中国自然资源航空物探遥感中心 基于张量特征值的Tilt法航磁边界检测方法、装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110133735A (zh) * 2019-05-16 2019-08-16 山东大学 深部异常体边界识别与分步增强方法、系统、介质及设备
CN111007571A (zh) * 2019-11-28 2020-04-14 吉林大学 一种基于三维构造张量的航磁数据地质体边界识别方法
CN111045089A (zh) * 2019-12-20 2020-04-21 核工业北京地质研究院 一种增强重磁场边界识别方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9383475B2 (en) * 2008-06-09 2016-07-05 Rock Solid Images, Inc. Geophysical surveying

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110133735A (zh) * 2019-05-16 2019-08-16 山东大学 深部异常体边界识别与分步增强方法、系统、介质及设备
CN111007571A (zh) * 2019-11-28 2020-04-14 吉林大学 一种基于三维构造张量的航磁数据地质体边界识别方法
CN111045089A (zh) * 2019-12-20 2020-04-21 核工业北京地质研究院 一种增强重磁场边界识别方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Enhancement of the total horizontal gradient of magnetic anomalies using the tilt angle;Francisco J. F. Ferreira et al.;《GEOPHYSICS》;20130630;第78卷(第3期);第J33-J41页 *
NAV-Edge: Edge detection of potential-field sources using normalized anisotropy variance;Heng Lei Zhang et al.;《GEOPHYSICS》;20140630;第79卷(第3期);第J43-J53页 *
Tilt-Euler方法在位场数据处理及解释中的应用;王明等;《物探与化探》;20120229;第36卷(第1期);第126-131页 *
重磁场二维希尔伯特变换-直接解析信号解释方法;骆遥等;《地球物理学报》;20110731;第54卷(第7期);第1912-1919页 *

Also Published As

Publication number Publication date
CN112666615A (zh) 2021-04-16

Similar Documents

Publication Publication Date Title
Oruç Edge detection and depth estimation using a tilt angle map from gravity gradient data of the Kozaklı-Central Anatolia Region, Turkey
Dannemiller et al. A new method for determination of magnetization direction
CN113886753B (zh) 基于张量特征值的Tilt法航磁边界检测方法、装置
CN106154323A (zh) 基于地震拓频处理的相控随机反演薄储层预测方法
CN109799540B (zh) 基于地质信息约束下的火山岩型铀矿床磁化率反演方法
CN112666615B (zh) 基于Tilt法获取精确航空磁测数据的方法、装置
Zhou et al. Self-constrained Euler deconvolution using potential field data of different altitudes
Williams et al. Comparison of grid Euler deconvolution with and without 2D constraints using a realistic 3D magnetic basement model
Ding et al. Reliability analysis of seismic attribute in the detection of fault-karst
Ma et al. Interpretation of magnetic anomalies by horizontal and vertical derivatives of the analytic signal
Abdelrahman et al. Magnetic interpretation using a least-squares, depth-shape curves method
CN108169802A (zh) 一种粗糙介质模型的时域电磁数据慢扩散成像方法
Zhu et al. Normalized vertical derivatives in the edge enhancement of maximum-edge-recognition methods in potential fields
CN112666619B (zh) 基于标准偏差法取精确航空磁测数据的方法、装置
Schiffler et al. Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data
CN113886754B (zh) 基于张量特征值的Theta Map法航磁边界检测方法、装置
CN115236755B (zh) 基于张量特征值的航磁异常边界检测方法、装置
CN115857028A (zh) 基于非接触观测法的水文地质参数的时空定量表达方法
Guo et al. Three-directional analytic signal analysis and interpretation of magnetic gradient tensor
Guo et al. Apparent magnetization mapping in the presence of strong remanent magnetization: The space-domain inversion approach
Pan et al. Euler deconvolution of the analytic signals of the gravity gradient tensor for the horizontal pipeline of finite length by horizontal cylinder calculation
CN113552623A (zh) 曲流河沉积砂体边界的确定方法和确定装置
Mendonça Inversion of gravity-field inclination to map the basement relief of sedimentary basins
Sun et al. Scratch recognition and analysis of gravity field in Chinese continent
Gerovska et al. Three-dimensional interpretation of magnetic and gravity anomalies using the finite-difference similarity transform

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