CN112666619B - 基于标准偏差法取精确航空磁测数据的方法、装置 - Google Patents

基于标准偏差法取精确航空磁测数据的方法、装置 Download PDF

Info

Publication number
CN112666619B
CN112666619B CN202011503984.7A CN202011503984A CN112666619B CN 112666619 B CN112666619 B CN 112666619B CN 202011503984 A CN202011503984 A CN 202011503984A CN 112666619 B CN112666619 B CN 112666619B
Authority
CN
China
Prior art keywords
standard deviation
gradient
data
horizontal gradient
wmean
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
CN202011503984.7A
Other languages
English (en)
Other versions
CN112666619A (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 CN202011503984.7A priority Critical patent/CN112666619B/zh
Publication of CN112666619A publication Critical patent/CN112666619A/zh
Application granted granted Critical
Publication of CN112666619B publication Critical patent/CN112666619B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Magnetic Variables (AREA)

Abstract

本发明公开了基于标准偏差法取精确航空磁测数据的方法、装置，方法包括：根据航空磁测数据T，计算航空磁测数据T在x方向水平梯度
Figure DDA0002844333020000011
在y方向水平梯度
Figure DDA0002844333020000012
和在z方向垂向梯度
Figure DDA0002844333020000013
分别计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0；根据z方向垂向梯度的标准偏差Sz计算z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及根据x方向、y方向、z方向的标准偏差之和S0计算标准偏差S0平均值Wmean(S0)；根据上述计算获取精确航空磁测数据。本发明提供的基于标准偏差法取精确航空磁测数据的方法、装置，可以解决NSTD法存在“解析奇点”问题，提高计算稳定性，增强信噪比，以及提高航磁数据处理转换的质量和地质体边界识别效果。

Description

基于标准偏差法取精确航空磁测数据的方法、装置
技术领域
本发明是涉及航空磁测技术领域，特别是涉及一种基于标准偏差法取精确航空磁测数据的方法、装置。
背景技术
航空磁力测量是将航空磁力仪(例如光泵式、核旋式和磁通门式)系统安装在飞行器，中通过观测地磁场参数(例如地磁场总强度T或总磁场异常△T或其梯度)，来寻找磁性或与磁性有关的矿体，以了解地质构造、进行磁性填图、解决城市和工程稳定性和考古等问题。
航空磁力测量数据是不同深度、不同形态、不同规模的磁性地质体磁场信息在观测面上的综合反映。但是，由于测量数据的误差或磁场的叠加，使得测量数据难以区分，给地质解释工作带来了难度。
目前，可以通过边界识别滤波器归一化标准偏差法(NSTD)的方法，提高对航磁异常的分辨能力。现有的归一化标准偏差法表达式为：
Figure BDA0002844331000000011
其中，
Figure BDA0002844331000000012
S0＝Sx+Sy+Sz，Wmean()表示滑动窗口平均值，T为实测的航空磁测数据，x、y为空间坐标的两个方向，Δ表示正则化调节因子，不同的取值对边界识别结果有影响；Δ的取值0-1，其通过深度学习样本库的自动获取方法，其受不同地质条件产生变化，根据地质解释结果的优劣确定最佳参数值；Δ取值过小会产生虚假多余边界。
基于此，本申请的发明人发现，现有归一化标准偏差法在实际应用中存在分辨力低、精度差和易产生虚假异常边界，且当分母
Figure BDA0002844331000000021
或接近0时，归一化标准偏差法存在“解析奇点”，会使得计算结果不稳定，影响实际应用效果。
公开于该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域一般技术人员所公知的现有技术。
发明内容
为解决上述问题，本发明实施例的目的在于提供基于标准偏差法获取精确航空磁测数据的方法、装置。
为实现上述目的，本发明提供了基于标准偏差法获取精确航空磁测数据的方法，包括：根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000022
所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000023
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000024
根据所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000025
以及所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000026
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000027
分别计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0；根据所述z方向垂向梯度的标准偏差Sz计算z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及根据所述x方向、y方向、z方向的标准偏差之和S0计算标准偏差S0平均值Wmean(S0)；根据x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz、z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)和标准偏差S0平均值Wmean(S0)，获取精确航空磁测数据。
在一优选的实施方式中，所述根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000031
所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000032
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000033
之前，包括：获取实测的航空磁测数据T，并进行预处理。
在一优选的实施方式中，所述预处理包括坐标转换、正常场校正、日变校正、滞后校正以及磁场水平调整中的一种或多种。
在一优选的实施方式中，所述根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000034
所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000035
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000036
包括：根据预处理后的航空磁测数据T进行频率域内二维傅里叶变换；根据二维傅里叶变换与希尔伯特变换之间的关系进行计算，得到频率域希尔伯特变换；根据所述频率域希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000037
预处理后的航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000038
和z方向垂向梯度
Figure BDA0002844331000000039
在一优选的实施方式中，所述根据所述希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure BDA00028443310000000310
预处理后的航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000041
和z方向垂向梯度
Figure BDA0002844331000000042
包括：二维希尔伯特变换在x方向上的分量对应为
Figure BDA0002844331000000043
二维希尔伯特变换在y方向上的分量对应为
Figure BDA0002844331000000044
预处理后的T对应为
Figure BDA0002844331000000045
在一优选的实施方式中，根据所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000046
以及所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000047
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000048
分别计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0包括：提供大小为m×n的滑动窗口；在所述滑动窗口内，根据以下公式计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0，
所述公式包括：
Figure BDA0002844331000000049
Figure BDA00028443310000000410
Figure BDA00028443310000000411
Figure BDA00028443310000000412
其中，μx为滑动窗口内x方向航空磁测数据的算术平均值，μy为滑动窗口内x方向航空磁测数据的算术平均值，μz为滑动窗口内x方向航空磁测数据的算术平均值。
在一优选的实施方式中，根据所述z方向垂向梯度的标准偏差Sz计算z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及根据所述x方向、y方向、z方向的标准偏差之和S0计算标准偏差S0平均值Wmean(S0)包括：
在所述滑动窗口内计算窗口内所有元素在z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及在所述滑动窗口内计算窗口内所有元素的标准偏差S0平均值Wmean(S0)。
在一优选的实施方式中，所述根据x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz、z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)和标准偏差S0平均值Wmean(S0)，获取精确航空磁测数据包括：
根据以下公式获取精确航空磁测数据，所述公式包括：
Figure BDA0002844331000000051
其中，其中，
Figure BDA0002844331000000052
S0＝Sx+Sy+Sz，Wmean()表示滑动窗口均值，T为实测的航空磁测数据，x、y为空间坐标的两个方向，Δ为正则化调节因子。
在一优选的实施方式中，所述正则化调节因子为0.001。
为实现上述目的，本发明还提供了一种基于标准偏差法获取精确航空磁测数据的装置，包括：梯度计算模块，用于根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000061
所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000062
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000063
标准偏差计算模块，用于根据所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000064
以及所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000065
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000066
分别计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0；均值计算模块，用于根据所述z方向垂向梯度的标准偏差Sz计算z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及根据所述x方向、y方向、z方向的标准偏差之和S0计算标准偏差S0平均值Wmean(S0)；精确航空磁测数据获取模块，用于根据x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz、z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)和标准偏差S0平均值Wmean(S0)，获取精确航空磁测数据。
为实现上述目的，本发明还提供了一种存储介质，所述存储介质存储有计算机可执行指令，所述计算机可执行指令用于执行上述任意一项所述的基于标准偏差法获取精确航空磁测数据的方法。
与现有技术相比，根据本发明的基于标准偏差法获取精确航空磁测数据的方法、装置，可以解决NSTD法存在“解析奇点”问题，提高计算稳定性，消除产生虚假航磁数据地质体的边界干扰，增强信噪比，以及提高航磁数据处理转换的质量和地质体边界识别效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出了本发明实施例所提供的基于标准偏差法获取精确航空磁测数据的方法的流程图；
图2示出了本发明实施例所提供的基于标准偏差法获取精确航空磁测数据的装置的结构示意图。
具体实施方式
下面结合附图，对本发明的具体实施方式进行详细描述，但应当理解本发明的保护范围并不受具体实施方式的限制。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
除非另有其它明确表示，否则在整个说明书和权利要求书中，术语“包括”或其变换如“包含”或“包括有”等等将被理解为包括所陈述的元件或组成部分，而并未排除其它元件或其它组成部分。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。
本发明实施例提供的基于标准偏差法获取精确航空磁测数据的方法的流程图，参见图1，包括：步骤S1-步骤S4。
步骤1，根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000081
所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000082
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000083
其中，航空磁测数据T即航空磁力异常场是地壳中的含铁磁性地质体在地磁场作用下所产生的附加磁场。
在一种实现方式中，步骤S1之前，还可以包括：获取实测的航空磁测数据T，并进行预处理。所述预处理包括坐标转换、正常场校正、日变校正、滞后校正以及磁场水平调整中的一种或多种。通过预处理，可以消除在航空磁场测量数据中，因多种因素会造成数据存在噪声和误差。
在一种实现方式中，步骤S1进一步可以包括：
步骤S101，根据预处理后的航空磁测数据T进行频率域内二维傅里叶变换；
步骤S102，根据二维傅里叶变换与希尔伯特变换之间的关系进行计算，得到频率域希尔伯特变换；
步骤S103，根据所述频率域希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000084
预处理后的航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000085
和z方向垂向梯度
Figure BDA0002844331000000086
具体的，根据所述希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000091
预处理后的航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000092
和z方向垂向梯度
Figure BDA0002844331000000093
包括：二维希尔伯特变换在x方向上的分量对应为
Figure BDA0002844331000000094
二维希尔伯特变换在y方向上的分量对应为
Figure BDA0002844331000000095
预处理后的T对应为
Figure BDA0002844331000000096
二维希尔伯特变换可以在频率中直接计算，其傅里叶域变换因子为：
DH(u,v)＝-i·sign(u,v)，
Figure BDA0002844331000000097
其中，i2＝-1，u、v是x和y方向上的圆波数。希尔伯特变换在x和y方向上的分量可以表示为：
Figure BDA0002844331000000098
其中，F[]、F-1[]分别表示傅里叶正变换和逆变换；Hx与Hy分别代表二维希尔伯特变换在x和y方向上的分量。
因此，x方向水平梯度
Figure BDA0002844331000000099
y方向水平梯度
Figure BDA00028443310000000910
和z方向垂向梯度
Figure BDA00028443310000000911
可以利用希尔伯特变换计算得到，二者具有等价关系：
Figure BDA00028443310000000912
应用上述等价关系优点在于直接希尔伯特变换不会放大航空磁测数据中的噪声干扰。
步骤2，根据所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000101
以及所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000102
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000103
分别计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0。
具体的，提供大小为m×n的滑动窗口；在所述滑动窗口内，根据以下公式计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0，所述公式包括：
Figure BDA0002844331000000104
Figure BDA0002844331000000105
Figure BDA0002844331000000106
Figure BDA0002844331000000107
其中，μx为滑动窗口内x方向航空磁测数据的算术平均值，μy为滑动窗口内x方向航空磁测数据的算术平均值，μz为滑动窗口内x方向航空磁测数据的算术平均值。
步骤3，根据所述z方向垂向梯度的标准偏差Sz计算z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及根据所述x方向、y方向、z方向的标准偏差之和S0计算标准偏差S0平均值Wmean(S0)。
具体的，步骤3可以包括：
在所述滑动窗口内计算窗口内所有元素在z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及在所述滑动窗口内计算窗口内所有元素的标准偏差S0平均值Wmean(S0)。
具体的，可以通过以下公式进行计算：
Figure BDA0002844331000000111
Figure BDA0002844331000000112
其中，i表示m×n大小的滑动窗口内的第几个元素。
步骤4，根据x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz、z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)和标准偏差S0平均值Wmean(S0)，获取精确航空磁测数据。
步骤4可以根据以下公式获取精确航空磁测数据，所述公式包括：
Figure BDA0002844331000000113
其中，
Figure BDA0002844331000000114
S0＝Sx+Sy+Sz，Wmean()表示滑动窗口均值，T为预处理后的航空磁测数据，x、y为空间坐标的两个方向，Δ为正则化调节因子，Δ的取值范围为0-1。其中，正则化调节因子可以为0.001。
不同的正则化调节因子的取值对边界识别结果有影响。正则化调节因子是通过深度学习样本库的自动获取方法获得的，得到Δ的取值范围为0-1。Δ取值过小会产生虚假多余边界，Δ取值过大会使结果逐渐趋近与总水平导数其受不同地质条件产生变化，根据地质解释结果的优劣确定最佳参数值。作为最优的一种实现方式，正则化调节因子的取值可以为0.001。由此，可以提高边界的识别效果。
在根据x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz、z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)和标准偏差S0平均值Wmean(S0)构建NSTD滤波器之后，可以获取精确航空磁测数据，进一步准确推断构造体场源的边界、深度、产状、规模及场分布规律和物理性质等，对划分大地构造单元、进行构造分区、确定断裂构造带的位、区分不同岩性与地层的分布及进行物性填图等问题有重要意义。
通过引入正则化调节因子Δ和希尔伯特变换，建立了一个合理的均衡滤波器，重新定义了归一化标准偏差法，能更好的探测出埋深不同的多源场物体的边界，使边界识别结果更收敛，具有更高的分辨率和精度。
由此，解决NSTD法存在“解析奇点”问题，提高计算稳定性，消除产生虚假航磁数据地质体的边界干扰，增强信噪比，以及提高航磁数据处理转换的质量和地质体边界识别效果。
本发明实施例还提供了一种获取精确航空磁测数据的装置，包括：梯度计算模块1、标准偏差计算模块2、均值计算模块3以及精确航空磁测数据获取模块4。
梯度计算模块1用于根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000121
所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000122
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000123
标准偏差计算模块2用于根据所述航空磁测数据T在x方向水平梯度
Figure BDA0002844331000000131
以及所述航空磁测数据T在y方向水平梯度
Figure BDA0002844331000000132
和所述航空磁测数据T在z方向垂向梯度
Figure BDA0002844331000000133
分别计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0。
均值计算模块3用于根据所述z方向垂向梯度的标准偏差Sz计算z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及根据所述x方向、y方向、z方向的标准偏差之和S0计算标准偏差S0平均值Wmean(S0)。
精确航空磁测数据获取模块4用于根据x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz、z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)和标准偏差S0平均值Wmean(S0)，获取精确航空磁测数据。
本发明实施例还提供了一种存储介质，所述存储介质存储有计算机可执行指令，其包含用于执行上述基于标准偏差法获取精确航空磁测数据的方法的程序，该计算机可执行指令可执行上述任意方法实施例中的方法。
其中，所述存储介质可以是计算机能够存取的任何可用介质或数据存储设备，包括但不限于磁性存储器(例如软盘、硬盘、磁带、磁光盘(MO)等)、光学存储器(例如CD、DVD、BD、HVD等)、以及半导体存储器(例如ROM、EPROM、EEPROM、非易失性存储器(NAND FLASH)、固态硬盘(SSD))等。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

Claims (9)

1.一种基于标准偏差法获取精确航空磁测数据的方法，其特征在于，包括：
根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure FDA0003150450790000011
所述航空磁测数据T在y方向水平梯度
Figure FDA0003150450790000012
和所述航空磁测数据T在z方向垂向梯度
Figure FDA0003150450790000013
根据所述航空磁测数据T在x方向水平梯度
Figure FDA0003150450790000014
以及所述航空磁测数据T在y方向水平梯度
Figure FDA0003150450790000015
和所述航空磁测数据T在z方向垂向梯度
Figure FDA0003150450790000016
分别计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0；
根据所述z方向垂向梯度的标准偏差Sz计算z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及根据所述x方向、y方向、z方向的标准偏差之和S0计算标准偏差S0平均值Wmean(S0)；
根据x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz、z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)和标准偏差S0平均值Wmean(S0)，获取精确航空磁测数据；
所述根据x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz、z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)和标准偏差S0平均值Wmean(S0)，获取精确航空磁测数据包括：
根据以下公式获取精确航空磁测数据，所述公式包括：
Figure FDA0003150450790000021
其中，其中，
Figure FDA0003150450790000022
S0＝Sx+Sy+Sz，Wmean()表示滑动窗口均值，T为实测的航空磁测数据，x、y为空间坐标的两个方向，Δ为正则化调节因子，Δ的取值范围为0-1。
2.根据权利要求1所述的方法，其特征在于，所述根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure FDA0003150450790000023
所述航空磁测数据T在y方向水平梯度
Figure FDA0003150450790000024
和所述航空磁测数据T在z方向垂向梯度
Figure FDA0003150450790000025
之前，包括：
获取实测的航空磁测数据T，并进行预处理。
3.根据权利要求2所述的方法，其特征在于，所述预处理包括坐标转换、正常场校正、日变校正、滞后校正以及磁场水平调整中的一种或多种。
4.根据权利要求2所述的方法，其特征在于，所述根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure FDA0003150450790000026
所述航空磁测数据T在y方向水平梯度
Figure FDA0003150450790000027
和所述航空磁测数据T在z方向垂向梯度
Figure FDA0003150450790000028
包括：
根据预处理后的航空磁测数据T进行频率域内二维傅里叶变换；
根据二维傅里叶变换与希尔伯特变换之间的关系进行计算，得到频率域希尔伯特变换；
根据所述频率域希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure FDA0003150450790000031
预处理后的航空磁测数据T在y方向水平梯度
Figure FDA0003150450790000032
和z方向垂向梯度
Figure FDA0003150450790000033
5.根据权利要求4所述的方法，其特征在于，所述根据所述希尔伯特变换的结果确定预处理后的航空磁测数据T在x方向水平梯度
Figure FDA0003150450790000034
预处理后的航空磁测数据T在y方向水平梯度
Figure FDA0003150450790000035
和z方向垂向梯度
Figure FDA0003150450790000036
包括：
二维希尔伯特变换在x方向上的分量对应为
Figure FDA0003150450790000037
二维希尔伯特变换在y方向上的分量对应为
Figure FDA0003150450790000038
预处理后的T对应为
Figure FDA0003150450790000039
6.根据权利要求4所述的方法，其特征在于，根据所述航空磁测数据T在x方向水平梯度
Figure FDA00031504507900000310
以及所述航空磁测数据T在y方向水平梯度
Figure FDA00031504507900000311
和所述航空磁测数据T在z方向垂向梯度
Figure FDA00031504507900000312
分别计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0包括：
提供大小为m×n的滑动窗口；
在所述滑动窗口内，根据以下公式计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0，所述公式包括：
Figure FDA0003150450790000041
Figure FDA0003150450790000042
Figure FDA0003150450790000043
Figure FDA0003150450790000044
其中，μx为滑动窗口内x方向航空磁测数据的算术平均值，μy为滑动窗口内y方向航空磁测数据的算术平均值，μz为滑动窗口内z方向航空磁测数据的算术平均值。
7.根据权利要求6所述的方法，其特征在于，根据所述z方向垂向梯度的标准偏差Sz计算z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及根据所述x方向、y方向、z方向的标准偏差之和S0计算标准偏差S0平均值Wmean(S0)包括：
在所述滑动窗口内计算窗口内所有元素在z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及在所述滑动窗口内计算窗口内所有元素的标准偏差S0平均值Wmean(S0)。
8.根据权利要求1所述的方法，其特征在于，
所述正则化调节因子为0.001。
9.一种基于标准偏差法获取精确航空磁测数据的装置，其特征在于，包括：
梯度计算模块，用于根据航空磁测数据T，计算所述航空磁测数据T在x方向水平梯度
Figure FDA0003150450790000051
所述航空磁测数据T在y方向水平梯度
Figure FDA0003150450790000052
和所述航空磁测数据T在z方向垂向梯度
Figure FDA0003150450790000053
标准偏差计算模块，用于根据所述航空磁测数据T在x方向水平梯度
Figure FDA0003150450790000054
以及所述航空磁测数据T在y方向水平梯度
Figure FDA0003150450790000055
和所述航空磁测数据T在z方向垂向梯度
Figure FDA0003150450790000056
分别计算x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz及x方向、y方向、z方向的标准偏差之和S0；
均值计算模块，用于根据所述z方向垂向梯度的标准偏差Sz计算z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)，以及根据所述x方向、y方向、z方向的标准偏差之和S0计算标准偏差S0平均值Wmean(S0)；
精确航空磁测数据获取模块，用于根据x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz、z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)和标准偏差S0平均值Wmean(S0)，获取精确航空磁测数据；
所述根据x方向水平梯度的标准偏差Sx、y方向水平梯度的标准偏差Sy和z方向垂向梯度的标准偏差Sz、z方向垂向梯度的标准偏差Sz平均值Wmean(Sz)和标准偏差S0平均值Wmean(S0)，获取精确航空磁测数据包括：
根据以下公式获取精确航空磁测数据，所述公式包括：
Figure FDA0003150450790000061
其中，其中，
Figure FDA0003150450790000062
S0＝Sx+Sy+Sz，Wmean()表示滑动窗口均值，T为实测的航空磁测数据，x、y为空间坐标的两个方向，Δ为正则化调节因子，Δ的取值范围为0-1。
CN202011503984.7A 2020-12-17 2020-12-17 基于标准偏差法取精确航空磁测数据的方法、装置 Active CN112666619B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011503984.7A CN112666619B (zh) 2020-12-17 2020-12-17 基于标准偏差法取精确航空磁测数据的方法、装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011503984.7A CN112666619B (zh) 2020-12-17 2020-12-17 基于标准偏差法取精确航空磁测数据的方法、装置

Publications (2)

Publication Number Publication Date
CN112666619A CN112666619A (zh) 2021-04-16
CN112666619B true CN112666619B (zh) 2021-08-24

Family

ID=75406430

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011503984.7A Active CN112666619B (zh) 2020-12-17 2020-12-17 基于标准偏差法取精确航空磁测数据的方法、装置

Country Status (1)

Country Link
CN (1) CN112666619B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113866836B (zh) * 2021-11-23 2023-10-20 重庆大学 一种基于归一化磁异常导数标准差的多目标边界识别方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107797148A (zh) * 2017-10-18 2018-03-13 中国国土资源航空物探遥感中心 一种基于三维地质建模的航磁异常场分离方法及系统
CN109407161A (zh) * 2018-09-21 2019-03-01 中国国土资源航空物探遥感中心 用于提取地球物理磁异常场边界的磁场刻痕分析方法
CN109814163A (zh) * 2019-02-28 2019-05-28 中国科学院遥感与数字地球研究所 一种基于扩展补偿模型的航磁张量数据抑噪方法及系统
CN109856690A (zh) * 2019-02-28 2019-06-07 中国科学院遥感与数字地球研究所 基于混合范数拟合的航磁梯度张量数据抑噪方法及系统
CN110632674A (zh) * 2019-10-28 2019-12-31 中国自然资源航空物探遥感中心 一种航空重力测量数据的弱信息提取方法
CN111007571A (zh) * 2019-11-28 2020-04-14 吉林大学 一种基于三维构造张量的航磁数据地质体边界识别方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7705599B2 (en) * 2007-07-09 2010-04-27 Kjt Enterprises, Inc. Buoy-based marine electromagnetic signal acquisition system
US9933540B2 (en) * 2012-09-28 2018-04-03 Cgg Data Services Ag Multiple receivers for airborne electromagnetic surveying
US9335434B2 (en) * 2012-11-02 2016-05-10 Pgs Geophysical As Method and system for processing data acquired in an electromagnetic survey
KR101403296B1 (ko) * 2013-12-09 2014-06-03 한국지질자원연구원 3차원 항공 자력 탐사 시스템 및 이를 이용한 3차원 항공 자력 탐사 방법
CN107037492A (zh) * 2017-05-26 2017-08-11 贵州省地质矿产勘查开发局0三地质大队 一种地质数据分析建模方法
CN109975879B (zh) * 2019-03-29 2020-06-26 中国科学院电子学研究所 一种基于磁传感器阵列的磁偶极子目标跟踪方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107797148A (zh) * 2017-10-18 2018-03-13 中国国土资源航空物探遥感中心 一种基于三维地质建模的航磁异常场分离方法及系统
CN109407161A (zh) * 2018-09-21 2019-03-01 中国国土资源航空物探遥感中心 用于提取地球物理磁异常场边界的磁场刻痕分析方法
CN109814163A (zh) * 2019-02-28 2019-05-28 中国科学院遥感与数字地球研究所 一种基于扩展补偿模型的航磁张量数据抑噪方法及系统
CN109856690A (zh) * 2019-02-28 2019-06-07 中国科学院遥感与数字地球研究所 基于混合范数拟合的航磁梯度张量数据抑噪方法及系统
CN110632674A (zh) * 2019-10-28 2019-12-31 中国自然资源航空物探遥感中心 一种航空重力测量数据的弱信息提取方法
CN111007571A (zh) * 2019-11-28 2020-04-14 吉林大学 一种基于三维构造张量的航磁数据地质体边界识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Tilt-Euler 方法在位场数据处理及解释中的应用";王明 等;《物探与化探》;20120229;第36卷(第1期);第126-132页 *
"航磁软补偿质量评价方法及软件实现";王林飞 等;《物探与化探》;20131231;第37卷(第6期);第1027-1030页 *
"重磁场二维希尔伯特变换— 直接解析信号解释方法";骆遥 等;《地球物理学报》;20110731;第54卷(第7期);第1912-1920页 *

Also Published As

Publication number Publication date
CN112666619A (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
CN113886753B (zh) 基于张量特征值的Tilt法航磁边界检测方法、装置
Dannemiller et al. A new method for determination of magnetization direction
Pilkington et al. Potential field continuation between arbitrary surfaces—Comparing methods
CN103792576A (zh) 一种基于梯度结构张量的储层非均质检测方法及设备
CN110133716A (zh) 基于组合模型加权函数的磁异常数据三维反演方法
CN112666619B (zh) 基于标准偏差法取精确航空磁测数据的方法、装置
CN112749493A (zh) 基于全磁力梯度张量特征值的地质体边界检测方法及系统
Abdelrahman et al. Magnetic interpretation using a least-squares, depth-shape curves method
CN112666615B (zh) 基于Tilt法获取精确航空磁测数据的方法、装置
Ma et al. Interpretation of magnetic anomalies by horizontal and vertical derivatives of the analytic signal
Ding et al. Reliability analysis of seismic attribute in the detection of fault-karst
Zhu et al. Normalized vertical derivatives in the edge enhancement of maximum-edge-recognition methods in potential fields
CN113886754B (zh) 基于张量特征值的Theta Map法航磁边界检测方法、装置
CN112668146B (zh) 一种基于欧拉反褶积法实用性改进的场源位置估算方法
CN112328955B (zh) 重磁数据的处理方法、存储介质及装置
Liu et al. The IDQ curve: A tool for evaluating the direction of remanent magnetization from magnetic anomalies
CN115236755A (zh) 基于张量特征值的航磁异常边界检测方法、装置
CN114876454B (zh) 一种大斜度井复杂岩性水淹层识别方法及系统
CN113818867B (zh) 一种构建伪毛管压力曲线方法、系统、介质、设备及应用
CN105093309A (zh) 基于构造应力分析进行裂缝检测及有效性验证的方法
Guo et al. Apparent magnetization mapping in the presence of strong remanent magnetization: The space-domain inversion approach
Gerovska et al. Three-dimensional interpretation of magnetic and gravity anomalies using the finite-difference similarity transform
Sun et al. Scratch recognition and analysis of gravity field in Chinese continent
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