CN108508490A - 一种基于解析信号的磁张量梯度数据均衡边界识别方法 - Google Patents
一种基于解析信号的磁张量梯度数据均衡边界识别方法 Download PDFInfo
- Publication number
- CN108508490A CN108508490A CN201810185428.6A CN201810185428A CN108508490A CN 108508490 A CN108508490 A CN 108508490A CN 201810185428 A CN201810185428 A CN 201810185428A CN 108508490 A CN108508490 A CN 108508490A
- Authority
- CN
- China
- Prior art keywords
- analytic signal
- boundary recognition
- ratio
- magnetic
- boundary
- 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
Links
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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于解析信号的磁张量梯度数据均衡边界识别方法,包括如下步骤:S1、基于磁张量数据方向解析信号的特征,建立合理的解析信号比值函数来实现边界识别工作;S2、利用比值函数可均衡不同深度地质体的效应,通过方向解析信号不同阶导数比值的均衡边界识别技术,获得收敛的边界结果;S3、通过基于解析信号比值的均衡边界识别滤波器,获得地质体的范围;S4、基于解析信号导数比值的阶跃均衡边界识别滤波器增强边界识别结果的分辨率,本发明可准确和清晰地给定区域矿产的分布,且在降低倾斜磁化干扰的同时提高了对较深地质体的分辨率,对于深部矿产资源勘探具有较大的应用前景。
Description
技术领域
本发明属于地球科学技术领域,更具体地说,尤其涉及一种基于解析信号的磁张量梯度数据均衡边界识别方法。
背景技术
磁张量梯度是磁矢量在不同方向的变化率,能更好地完成场源体几何形态的描述。磁张量梯度数据具有高精度、高分辨率、多参数的优点,能更准确地描述磁性体的分布,在矿产资源勘探方面具有重要的用途。但是,磁张量梯度数据依旧会受到倾斜磁化的干扰而使异常形态不规则,为解释带来了困难。解析信号具有受倾斜磁化干扰小的优点,但是解析信号随着深度的增加衰减严重,因此单独采用不能很好地识别较深地质体边界。为此提出基于解析信号比值的均衡边界识别技术进行磁张量数据的解释,均衡边界识别法是利用不同方向解析信号的比值进行地质体的边界识别,其能有效地均衡不同深度地质体的响应,在降低倾斜磁化干扰的同时提高了对较深地质体的分辨率。
发明内容
本发明的目的是为了解决现有技术中存在的缺点,而提出的一种基于解析信号的磁张量梯度数据均衡边界识别方法。
为实现上述目的,本发明提供如下技术方案:一种基于解析信号的磁张量梯度数据均衡边界识别方法,包括如下步骤:
S1、基于磁张量数据方向解析信号的特征,建立解析信号比值函数来实现边界识别工作;
S2、利用比值函数均衡不同深度地质体的效应,通过方向解析信号不同阶导数比值的均衡边界识别技术,获得收敛的边界结果;
S3、通过基于解析信号比值的均衡边界识别滤波器,获得地质体的范围;
S4、基于解析信号导数比值的阶跃均衡边界识别滤波器增强边界识别结果的分辨率。
优选的,所述步骤S1中磁张量梯度数据是磁场M不同方向分量的变化率,矩阵形式为:
优选的,所述磁场M矩阵M每一行的平方和为解析信号,称作x、y、z方向解析信号,x、y、z方向解析信号的方程式为:
优选的,所述步骤S3中基于解析信号比值的均衡边界识别滤波器的公式为:
优选的,所述步骤S4中基于解析信号导数比值的阶跃均衡边界识别滤波器的公式为:
本发明的技术效果和优点:本发明提供的一种基于解析信号的磁张量梯度数据均衡边界识别方法,基于磁张量数据方向解析信号的特征,建立合理的解析信号比值函数来实现边界识别工作,该方法有效地避免了磁化方向对结果的干扰,利用比值函数可均衡不同深度地质体的效应,从而能很清晰地给定较深磁源体的分布,此外,为了获得更加清晰地地质体边界,给出了方向解析信号不同阶导数比值的均衡边界识别技术,其能获得更加收敛的边界结果,本发明可准确和清晰地给定区域矿产的分布,且在降低倾斜磁化干扰的同时提高了对较深地质体的分辨率,对于深部矿产资源勘探具有较大的应用前景。
附图说明
图1为本发明一种基于解析信号的磁张量梯度数据均衡边界识别方法工艺流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参阅图1,一种基于解析信号的磁张量梯度数据均衡边界识别方法,包括如下步骤:
S1、基于磁张量数据方向解析信号的特征,建立合理的解析信号比值函数来实现边界识别工作,可以有效地避免了磁化方向对结果的干扰;
S2、利用比值函数可均衡不同深度地质体的效应,能很清晰地给定较深磁源体的分布,对于深部矿产资源勘探具有较大的应用前景;为了获得更加清晰地地质体边界,通过方向解析信号不同阶导数比值的均衡边界识别技术,获得收敛的边界结果;
S3、通过基于解析信号比值的均衡边界识别滤波器,获得地质体的范围,采用均衡边界识别滤波器BAS能同时显示不同深度地质体的位置,能更加清晰和准确地获得地质体的边界,提高了较深地质体的分辨率;
S4、基于解析信号导数比值的阶跃均衡边界识别滤波器增强边界识别结果的分辨率,采用阶跃边界识别滤波器,可以使边界识别结果更加收敛,易于完成数据的解释工作。
具体的,所述步骤S1中磁张量梯度数据是磁场M不同方向分量的变化率,矩阵形式为:
具体的,所述磁场M矩阵M每一行的平方和为解析信号,称作x、y、z方向解析信号,x、y、z方向解析信号的方程式为:
具体的,所述步骤S3中基于解析信号比值的均衡边界识别滤波器的公式为:
具体的,所述步骤S4中基于解析信号导数比值的阶跃均衡边界识别滤波器的公式为:
传统的是利用水平方向解析信号和水平方向解析信号垂直导数进行张量数据的边界识别,其表达式为:
传统的两种方法的极值对应于地质体的边界。通过埋深不同的棱柱体模型所产生的磁张量异常来试验不同方法的应用效果,采用方向解析信号进行地质体边界的识别,能清晰地显示较浅地质体的水平范围,且解析信号垂直导数的识别结果更加收敛,但两种方法对于较深地质体的范围均难以识别。
综上所述:本发明提供的一种基于解析信号的磁张量梯度数据均衡边界识别方法,基于磁张量数据方向解析信号的特征,建立合理的解析信号比值函数来实现边界识别工作,该方法有效地避免了磁化方向对结果的干扰,利用比值函数可均衡不同深度地质体的效应,从而能很清晰地给定较深磁源体的分布,此外,为了获得更加清晰地地质体边界,给出了方向解析信号不同阶导数比值的均衡边界识别技术,其能获得更加收敛的边界结果,本发明可准确和清晰地给定区域矿产的分布,且在降低倾斜磁化干扰的同时提高了对较深地质体的分辨率,对于深部矿产资源勘探具有较大的应用前景。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于解析信号的磁张量梯度数据均衡边界识别方法,其特征在于,包括如下步骤:
S1、基于磁张量数据方向解析信号的特征,建立解析信号比值函数来实现边界识别工作;
S2、利用比值函数均衡不同深度地质体的效应,通过方向解析信号不同阶导数比值的均衡边界识别技术,获得收敛的边界结果;
S3、通过基于解析信号比值的均衡边界识别滤波器,获得地质体的范围;
S4、基于解析信号导数比值的阶跃均衡边界识别滤波器增强边界识别结果的分辨率。
2.根据权利要求1所述的一种基于解析信号的磁张量梯度数据均衡边界识别方法,其特征在于:所述步骤S1中磁张量梯度数据是磁场M不同方向分量的变化率,矩阵形式为:
3.根据权利要求2所述的一种基于解析信号的磁张量梯度数据均衡边界识别方法,其特征在于:所述磁场M矩阵M每一行的平方和为解析信号,称作x、y、z方向解析信号,x、y、z方向解析信号的方程式为:
4.根据权利要求1所述的一种基于解析信号的磁张量梯度数据均衡边界识别方法,其特征在于:所述步骤S3中基于解析信号比值的均衡边界识别滤波器的公式为:
5.根据权利要求1所述的一种基于解析信号的磁张量梯度数据均衡边界识别方法,其特征在于:所述步骤S4中基于解析信号导数比值的阶跃均衡边界识别滤波器的公式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810185428.6A CN108508490B (zh) | 2018-03-07 | 2018-03-07 | 一种基于解析信号的磁张量梯度数据均衡边界识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810185428.6A CN108508490B (zh) | 2018-03-07 | 2018-03-07 | 一种基于解析信号的磁张量梯度数据均衡边界识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108508490A true CN108508490A (zh) | 2018-09-07 |
CN108508490B CN108508490B (zh) | 2019-12-27 |
Family
ID=63377021
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810185428.6A Expired - Fee Related CN108508490B (zh) | 2018-03-07 | 2018-03-07 | 一种基于解析信号的磁张量梯度数据均衡边界识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108508490B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109707375A (zh) * | 2018-12-26 | 2019-05-03 | 吉林大学 | 一种通过井下磁梯度张量数据辨识地下磁层信息的方法 |
CN110133735A (zh) * | 2019-05-16 | 2019-08-16 | 山东大学 | 深部异常体边界识别与分步增强方法、***、介质及设备 |
CN110687610A (zh) * | 2019-09-19 | 2020-01-14 | 长安大学 | 一种基于重磁数据相关性分析的场源定位及属性识别方法 |
CN111007571A (zh) * | 2019-11-28 | 2020-04-14 | 吉林大学 | 一种基于三维构造张量的航磁数据地质体边界识别方法 |
CN112327383A (zh) * | 2020-11-25 | 2021-02-05 | 吉林大学 | 基于梯度张量特征值的位场数据地质构造位置圈定方法 |
CN112749493A (zh) * | 2021-01-25 | 2021-05-04 | 中国自然资源航空物探遥感中心 | 基于全磁力梯度张量特征值的地质体边界检测方法及*** |
CN113886753A (zh) * | 2021-10-09 | 2022-01-04 | 中国自然资源航空物探遥感中心 | 基于张量特征值的Tilt法航磁边界检测方法、装置 |
CN113886754A (zh) * | 2021-10-09 | 2022-01-04 | 中国自然资源航空物探遥感中心 | 基于张量特征值的Theta Map法航磁边界检测方法、装置 |
CN116774303A (zh) * | 2023-06-26 | 2023-09-19 | 中国自然资源航空物探遥感中心 | 一种场源边界定位方法、装置及计算机可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090179646A1 (en) * | 2008-01-16 | 2009-07-16 | Sascha Fath | Gradient coil and method to manufacture a gradient coil |
CN101904745A (zh) * | 2009-06-05 | 2010-12-08 | 复旦大学 | 一种基于磁共振弥散张量成像确定目标边界的方法 |
CN106019394A (zh) * | 2016-04-27 | 2016-10-12 | 中国地质科学院矿产资源研究所 | 海洋大地电磁场非线性共轭梯度三维并行反演方法 |
CN106405658A (zh) * | 2016-08-30 | 2017-02-15 | 中国人民解放军海军工程大学 | 一种基于矢量磁梯度计的运动式磁性目标定位方法 |
-
2018
- 2018-03-07 CN CN201810185428.6A patent/CN108508490B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090179646A1 (en) * | 2008-01-16 | 2009-07-16 | Sascha Fath | Gradient coil and method to manufacture a gradient coil |
CN101904745A (zh) * | 2009-06-05 | 2010-12-08 | 复旦大学 | 一种基于磁共振弥散张量成像确定目标边界的方法 |
CN106019394A (zh) * | 2016-04-27 | 2016-10-12 | 中国地质科学院矿产资源研究所 | 海洋大地电磁场非线性共轭梯度三维并行反演方法 |
CN106405658A (zh) * | 2016-08-30 | 2017-02-15 | 中国人民解放军海军工程大学 | 一种基于矢量磁梯度计的运动式磁性目标定位方法 |
Non-Patent Citations (1)
Title |
---|
马国庆 等: ""磁张量数据的边界识别和解释方法"", 《石油地球物理勘探》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109707375A (zh) * | 2018-12-26 | 2019-05-03 | 吉林大学 | 一种通过井下磁梯度张量数据辨识地下磁层信息的方法 |
CN110133735A (zh) * | 2019-05-16 | 2019-08-16 | 山东大学 | 深部异常体边界识别与分步增强方法、***、介质及设备 |
CN110687610A (zh) * | 2019-09-19 | 2020-01-14 | 长安大学 | 一种基于重磁数据相关性分析的场源定位及属性识别方法 |
CN110687610B (zh) * | 2019-09-19 | 2021-03-19 | 长安大学 | 一种基于重磁数据相关性分析的场源定位及属性识别方法 |
CN111007571B (zh) * | 2019-11-28 | 2022-05-13 | 吉林大学 | 一种基于三维构造张量的航磁数据地质体边界识别方法 |
CN111007571A (zh) * | 2019-11-28 | 2020-04-14 | 吉林大学 | 一种基于三维构造张量的航磁数据地质体边界识别方法 |
CN112327383A (zh) * | 2020-11-25 | 2021-02-05 | 吉林大学 | 基于梯度张量特征值的位场数据地质构造位置圈定方法 |
CN112327383B (zh) * | 2020-11-25 | 2021-12-31 | 吉林大学 | 基于梯度张量特征值的位场数据地质构造位置圈定方法 |
CN112749493A (zh) * | 2021-01-25 | 2021-05-04 | 中国自然资源航空物探遥感中心 | 基于全磁力梯度张量特征值的地质体边界检测方法及*** |
CN113886754A (zh) * | 2021-10-09 | 2022-01-04 | 中国自然资源航空物探遥感中心 | 基于张量特征值的Theta Map法航磁边界检测方法、装置 |
CN113886754B (zh) * | 2021-10-09 | 2022-04-15 | 中国自然资源航空物探遥感中心 | 基于张量特征值的Theta Map法航磁边界检测方法、装置 |
CN113886753A (zh) * | 2021-10-09 | 2022-01-04 | 中国自然资源航空物探遥感中心 | 基于张量特征值的Tilt法航磁边界检测方法、装置 |
CN116774303A (zh) * | 2023-06-26 | 2023-09-19 | 中国自然资源航空物探遥感中心 | 一种场源边界定位方法、装置及计算机可读存储介质 |
CN116774303B (zh) * | 2023-06-26 | 2024-05-07 | 中国自然资源航空物探遥感中心 | 一种场源边界定位方法、装置及计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108508490B (zh) | 2019-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108508490A (zh) | 一种基于解析信号的磁张量梯度数据均衡边界识别方法 | |
US20200326389A1 (en) | Method for Making a Magnetic Gradiometer with High Detection Accuracy and Success Rate | |
CN104374385B (zh) | 一种海底磁传感器阵列目标定位的新方法 | |
CN106405664B (zh) | 一种磁异常化极方法 | |
Yuan et al. | Edge detection of potential field data with improved structure tensor methods | |
CN105259571B (zh) | 一种地层倾角检测方法 | |
CN105589108A (zh) | 基于不同约束条件的瞬变电磁快速三维反演方法 | |
Wang et al. | Improved curvature gravity gradient tensor with principal component analysis and its application in edge detection of gravity data | |
Li et al. | An approach for estimating the magnetization direction of magnetic anomalies | |
CN109425906A (zh) | 一种磁异常探测矢量磁目标识别方法 | |
CN108072906A (zh) | 一种分布式磁探测磁目标识别方法 | |
Luo et al. | An ensemble-based framework for proactive geosteering | |
CN109298464A (zh) | 致密砂岩储层成岩相测井识别方法及装置 | |
Wang et al. | New improvements for lineaments study of gravity data with improved Euler inversion and phase congruency of the field data | |
ATE422245T1 (de) | Seismisches aufnahmesystem | |
CN107918153A (zh) | 一种地震信号相干性高精度检测方法 | |
Wang et al. | PS-wave processing and S-wave velocity inversion of OBS data from Northern South China Sea | |
CN104570104B (zh) | 一种基于两步法avf的纵横波地震品质因子提取方法 | |
CN109521470A (zh) | 分析地质构造对地震反演裂缝密度影响的方法 | |
CN113900141B (zh) | 油气分布预测方法及装置 | |
CN106353801B (zh) | 三维Laplace域声波方程数值模拟方法及装置 | |
Mushayandebvu et al. | Magnetic imaging using extended Euler deconvolution | |
Li et al. | Research on RTP aeromagnetic gradient data and its applicability in different latitudes | |
Fullagar et al. | 3D magnetic modelling and inversion incorporating self-demagnetisation and interactions | |
Zhou et al. | Magnetic anomaly depth and structural index estimation using different height analytic signals 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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191227 |
|
CF01 | Termination of patent right due to non-payment of annual fee |