CN108983296A - 一种瞬变电磁和激发极化效应解耦方法 - Google Patents
一种瞬变电磁和激发极化效应解耦方法 Download PDFInfo
- Publication number
- CN108983296A CN108983296A CN201810629059.5A CN201810629059A CN108983296A CN 108983296 A CN108983296 A CN 108983296A CN 201810629059 A CN201810629059 A CN 201810629059A CN 108983296 A CN108983296 A CN 108983296A
- Authority
- CN
- China
- Prior art keywords
- induced polarization
- transient electromagnetic
- effect
- signals
- polarization effect
- 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
- 230000010287 polarization Effects 0.000 title claims abstract description 67
- 230000001052 transient effect Effects 0.000 title claims abstract description 66
- 230000000694 effects Effects 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000012880 independent component analysis Methods 0.000 claims abstract description 23
- 238000001514 detection method Methods 0.000 claims abstract description 7
- 239000000284 extract Substances 0.000 claims abstract description 3
- 238000000926 separation method Methods 0.000 claims description 29
- 239000011159 matrix material Substances 0.000 claims description 23
- 238000005457 optimization Methods 0.000 claims description 18
- 230000001808 coupling effect Effects 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 description 6
- 230000005284 excitation Effects 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- UCKMPCXJQFINFW-UHFFFAOYSA-N Sulphide Chemical compound [S-2] UCKMPCXJQFINFW-UHFFFAOYSA-N 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
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/02—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current
-
- 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
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种瞬变电磁和激发极化效应解耦方法，将接收信号的分成直流分量与交流分量，以交流分量为混合信号，利用独立分量分析算法实现瞬变电磁和激发极化效应解耦，分别提取单纯的瞬变电磁信号与激发极化信号，然后分别进行三维反演和解释。本发明方法能能够有效的提取地矿探测中低阻地质异常信息、压制激发极化效应，可以有效提高矿体勘探中瞬变电磁与激发极化的综合解释效果。
Description
技术领域
本发明属于地球物理探测技术领域，具体涉及一种瞬变电磁和激发极化效应解耦方法。
背景技术
地矿勘探中，大多数成矿岩石受激发极化效应的影响，电阻率是一个与发射频率、岩矿石电化学性质等因素有关的复数，这就导致了野外测量所得的电磁场信号包含电磁感应效应与激发极化效应双重因素影响；这一现象主要表现为在金属矿或硫化矿产等激发极化效应较强的地区进行瞬变电磁勘探时，瞬变电磁采集的晚期信号出现负值的现象，即瞬变电磁的激发极化效应影响。
理论上若无激发极化效应影响，瞬变电磁响应为一严格衰减的正值信号，但由于耦合效应的存在，瞬变电磁晚期信号往往不满足衰减规律，并伴随大量负值，这使得传统的解释方法完全失效，导致采集数据完全报废。
因此，探索瞬变电磁与激发极化响应的解耦，实现瞬变电磁与激发极化的联合解释成为矿产资源勘探研究的焦点问题，但由于电磁耦合的复杂性，目前尚没有一种具有普适性的方法。
发明内容
针对上述现有技术存在的不足，本发明的目的在于，提供一种瞬变电磁和激发极化效应解耦与联合解释方法，用以实现瞬变电磁与激发极化的联合解释。
为了实现上述任务，本发明采用以下技术方案：
一种瞬变电磁和激发极化效应解耦方法，包括以下步骤：
将采集到的信号分解为直流分量和交流分量，提取的交流分量反映了瞬变电磁和激发极化效应的耦合效果，以此作为独立分量分析的信号源；
采用最小互信息独立分量分析算法进行激发极化效应分离，提取探测区域瞬变电磁信号与激发极化效应信号；
对于提取的激发极化效应信号，基于电阻率法进行激发极化三维反演；对于提取的瞬变电磁信号，采用瞬变电磁解释技术进行瞬变电磁解释。
进一步地，所述的最小互信息独立分量分析算法采用互信息极小化判据、基于自然梯度优化，包括两个方面：优化判据和寻优算法。
进一步地，所述的优化判据采用互信息极小化判据，表示为：
上式中，x是交流分量，y是输出信号，yi是输出信号y中的元素，m是输出信号个数，W是分离矩阵，log|W|为分离矩阵信息度量，H(·)为信息熵。
进一步地，所述的采用最小互信息独立分量分析算法采用自然梯度法进行优化，通过将常规梯度用自然梯度代替，从而简化调节参数时的分离矩阵求逆计算，保证计算的稳定性。
进一步地，通过逐步调节分离矩阵W使目标函数取得极小值完成寻优过程；调节分离矩阵的调节量ΔW表示为：
上式中，W是分离矩阵，μ为调节参数。
一种瞬变电磁和激发极化效应解耦系统，包括依次连接的信号采集分解模块、分离模块以及反演解释模块，其中：
所述的信号采集分解模块用于将采集到的信号分解为直流分量和交流分量，提取的交流分量反映了瞬变电磁和激发极化效应的耦合效果，以此作为独立分量分析的信号源；
所述的分离模块采用最小互信息独立分量分析算法进行激发极化效应分离，提取探测区域瞬变电磁信号与激发极化效应信号；
所述的反演解释模块用于对于提取的激发极化效应信号，基于电阻率法进行激发极化三维反演；对于提取的瞬变电磁信号，采用瞬变电磁解释技术进行瞬变电磁解释。
本发明具有以下技术特点：
1.本发明方法提出利用独立分量分析法对两种信号进行分离，得到单纯的瞬变电磁响应与激发极化响应，然后分别反演进行综合解释，同时给出极化率、电阻率以及深度信息，实现含有激发极化效应的瞬变电磁响应信息的综合利用。
2.本发明利用最小互信息独立分量分析算法进行激发极化效应解耦，实现矿体特征信号提取，提高了瞬变电磁法在矿体勘探中的有效应用，综合利用了瞬变电磁响应与激电特征。
附图说明
图1为本发明方法的示意图。
具体实施方式
参照图1，一种瞬变电磁和激发极化效应解耦方法包括以下步骤：
步骤1，将采集到的信号分解为直流分量和交流分量，提取的交流分量反映了瞬变电磁和激发极化效应的耦合效果，以此作为独立分量分析的信号源；
独立分量分析法要求源信号是零均值的实随机变量且在任意时刻均相互统计独立，为了满足这一条件限制，首先将采集到的信号进行分解，表达为直流分量和交流分量和的形式，而提取的交流分量反映了瞬变电磁和激发极化效应的耦合效果，以此作为独立分量分析的信号源。
步骤2，采用互信息极小化判据、基于自然梯度优化的最小互信息独立分量分析算法进行激发极化效应分离，从而提取探测区域瞬变电磁信号与激发极化效应信号；而所述的瞬变电磁信号即可反映出地质特征。
所述采用互信息极小化判据含义如下：
设n维列向量y的概率密度函数为Py(y)，各分量yi的概率密度函数为Pi(yi)，则互信息极小化判据为互信息I(y，W)取得极小值：
其中，是Py(y)与之间的KL散度(相对熵)，KL值越大，则二者越不相似，当两者相同时，KL值为0。
将交流分量信号作为混合信号输入，利用独立分量分析算法进行激发极化效应分离，实现激发极化效应信号与瞬变电磁信号的分离。独立分量分析算法包括两个方面：优化判据和寻优算法。
所述的优化判据采用互信息极小化判据，表示为：
其中，x是混合信号，即所述的交流分量，y是输出信号，yi是输出信号y中的元素，m是输出信号个数，W是分离矩阵，log|W|为分离矩阵信息度量，H(·)为信息熵。
所述的采用最小互信息独立分量分析算法采用自然梯度法进行优化，通过将常规梯度用自然梯度代替，从而简化调节参数时的分离矩阵求逆计算，保证计算的稳定性。由于目标函数I(y,W)是分离矩阵W的函数，所述的寻优算法为，通过逐步调节分离矩阵W使目标函数取得极小值完成寻优过程，所对应的输出信号y即为盲信号分离后的输出信号。调节分离矩阵W的调节量如式3所示。将式3展开，最后的调节参数涉及到对分离矩阵求逆。矩阵求逆往往计算量大；若矩阵条件数较大，往往计算不稳定。采用自然梯度法进行优化，将常规梯度用自然梯度代替，能够避免对分离矩阵求逆，不仅节省计算量，而且计算稳定。当然，其中调节参数中的参量需要转化成分离信号的高阶统计量来计算。
其中，y为输出信号，μ为调节参数，W是式2中的分离矩阵。
采用互信息极小化判据、基于自然梯度优化的最小互信息独立分量分析算法对输出的各独立成分的排序和幅度比例变化具有不变性。因此可以容易区分瞬变电磁信号与激电信号(激发极化效应信号)。
步骤3，对于提取的激发极化响应(激发极化效应信号)，基于电阻率法进行激发极化三维反演；对于提取的瞬变电磁信号，利用现有的瞬变电磁解释技术进行瞬变电磁解释。
将信号分离后，可分别获得无电磁影响的激电响应与无激电影响的瞬变电磁响应。将两者分别利用现有解释技术进行解释，获得瞬变电磁成果图件与激发极化成果图件，实现二者的综合解释。
综上所述，本实施例提供的一种瞬变电磁和激发极化效应解耦方法，将接收信号的分成直流分量与交流分量，以交流分量为混合信号，利用独立分量分析算法实现瞬变电磁和激发极化效应解耦，分别提取单纯的瞬变电磁信号与激发极化信号，更加准确提取矿体特征信号，可以提高矿体勘探中瞬变电磁与激发极化的综合解释效果。
Claims (6)
1.一种瞬变电磁和激发极化效应解耦方法，其特征在于，包括以下步骤：
将采集到的信号分解为直流分量和交流分量，提取的交流分量反映了瞬变电磁和激发极化效应的耦合效果，以此作为独立分量分析的信号源；
采用最小互信息独立分量分析算法进行激发极化效应分离，提取探测区域瞬变电磁信号与激发极化效应信号；
对于提取的激发极化效应信号，基于电阻率法进行激发极化三维反演；对于提取的瞬变电磁信号，采用瞬变电磁解释技术进行瞬变电磁解释。
2.如权利要求1所述的瞬变电磁和激发极化效应解耦方法，其特征在于，所述的最小互信息独立分量分析算法采用互信息极小化判据、基于自然梯度优化，包括两个方面：优化判据和寻优算法。
3.如权利要求2所述的瞬变电磁和激发极化效应解耦方法，其特征在于，所述的优化判据采用互信息极小化判据，表示为：
上式中，x是交流分量，y是输出信号，yi是输出信号y中的元素，m是输出信号个数，W是分离矩阵，log|W|为分离矩阵信息度量，H(·)为信息熵。
4.如权利要求2所述的瞬变电磁和激发极化效应解耦方法，其特征在于，所述的采用最小互信息独立分量分析算法采用自然梯度法进行优化，通过将常规梯度用自然梯度代替，从而简化调节参数时的分离矩阵求逆计算，保证计算的稳定性。
5.如权利要求3所述的瞬变电磁和激发极化效应解耦方法，其特征在于，通过逐步调节分离矩阵W使目标函数取得极小值完成寻优过程；调节分离矩阵的调节量ΔW表示为：
上式中，W是分离矩阵，μ为调节参数。
6.一种瞬变电磁和激发极化效应解耦系统，其特征在于，包括依次连接的信号采集分解模块、分离模块以及反演解释模块，其中：
所述的信号采集分解模块用于将采集到的信号分解为直流分量和交流分量，提取的交流分量反映了瞬变电磁和激发极化效应的耦合效果，以此作为独立分量分析的信号源；
所述的分离模块采用最小互信息独立分量分析算法进行激发极化效应分离，提取探测区域瞬变电磁信号与激发极化效应信号；
所述的反演解释模块用于对于提取的激发极化效应信号，基于电阻率法进行激发极化三维反演；对于提取的瞬变电磁信号，采用瞬变电磁解释技术进行瞬变电磁解释。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810629059.5A CN108983296B (zh) | 2018-06-19 | 2018-06-19 | 一种瞬变电磁和激发极化效应解耦方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810629059.5A CN108983296B (zh) | 2018-06-19 | 2018-06-19 | 一种瞬变电磁和激发极化效应解耦方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108983296A true CN108983296A (zh) | 2018-12-11 |
| CN108983296B CN108983296B (zh) | 2020-05-22 |
Family
ID=64540605
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810629059.5A Expired - Fee Related CN108983296B (zh) | 2018-06-19 | 2018-06-19 | 一种瞬变电磁和激发极化效应解耦方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108983296B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2021042952A1 (zh) * | 2019-09-05 | 2021-03-11 | 中国科学院地质与地球物理研究所 | 一种接地导线源瞬变电磁响应中ip信息的提取方法 |
| CN114047550A (zh) * | 2021-09-27 | 2022-02-15 | 中国地质科学院地球物理地球化学勘查研究所 | 一种压制慢极化效应的差分波形瞬变电磁测量方法 |
| CN115619786A (zh) * | 2022-12-19 | 2023-01-17 | 中国科学技术大学 | 磁场图像处理方法、装置 |
| CN118348598A (zh) * | 2024-06-18 | 2024-07-16 | 山东大学 | 一种隧道超前地质预报全时域激发极化探测系统及方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102053281A (zh) * | 2009-11-10 | 2011-05-11 | 中国石油化工集团公司 | 一种运用长偏移距瞬变电磁阵列法进行油气检测的方法 |
| CN105676295A (zh) * | 2016-01-30 | 2016-06-15 | 吉林大学 | 基于squid的磁源激发极化-感应的联合探测系统与方法 |
| CN107526109A (zh) * | 2017-08-16 | 2017-12-29 | 桂林电子科技大学 | 基于虚拟波域法的瞬变电磁建模和反演方法 |
-
2018
- 2018-06-19 CN CN201810629059.5A patent/CN108983296B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102053281A (zh) * | 2009-11-10 | 2011-05-11 | 中国石油化工集团公司 | 一种运用长偏移距瞬变电磁阵列法进行油气检测的方法 |
| CN105676295A (zh) * | 2016-01-30 | 2016-06-15 | 吉林大学 | 基于squid的磁源激发极化-感应的联合探测系统与方法 |
| CN107526109A (zh) * | 2017-08-16 | 2017-12-29 | 桂林电子科技大学 | 基于虚拟波域法的瞬变电磁建模和反演方法 |
Non-Patent Citations (9)
| Title |
|---|
| JIANMEI ZHOU ET AL.: "Simulation of 3D marine CSEM response in anisotropic media by coupled potential finite volume method", 《NSAPC 2015》 * |
| 冯兵等: "利用CSAMT电磁场响应提取激电效应的方法初探", 《地球物理学进展》 * |
| 孙怀凤等: "考虑关断时间的回线源激发TEM三维时域有限差分正演", 《地球物理学报》 * |
| 张发启等: "《盲信号处理及应用》", 31 October 2006, 西安：西安电子科技大学出版社 * |
| 徐信等: "基于EMD-ICA的激电数据降噪处理方法", 《计算机应用研究》 * |
| 王伟明: "考虑激电效应的瞬变电磁中心回线一维正反演", 《新疆有色金属》 * |
| 罗倩等: "瞬变电磁信号降噪算法", 《计算机仿真》 * |
| 邱宁: "瞬变电磁信号的预处理和激电效应识别研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 马建仓等: "《盲信号处理》", 30 June 2006, 国防工业出版社 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2021042952A1 (zh) * | 2019-09-05 | 2021-03-11 | 中国科学院地质与地球物理研究所 | 一种接地导线源瞬变电磁响应中ip信息的提取方法 |
| CN114047550A (zh) * | 2021-09-27 | 2022-02-15 | 中国地质科学院地球物理地球化学勘查研究所 | 一种压制慢极化效应的差分波形瞬变电磁测量方法 |
| CN114047550B (zh) * | 2021-09-27 | 2023-07-07 | 中国地质科学院地球物理地球化学勘查研究所 | 一种压制慢极化效应的差分波形瞬变电磁测量方法 |
| CN115619786A (zh) * | 2022-12-19 | 2023-01-17 | 中国科学技术大学 | 磁场图像处理方法、装置 |
| CN118348598A (zh) * | 2024-06-18 | 2024-07-16 | 山东大学 | 一种隧道超前地质预报全时域激发极化探测系统及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108983296B (zh) | 2020-05-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108983296B (zh) | 一种瞬变电磁和激发极化效应解耦方法 | |
| Huang et al. | Mathematical morphological filtering for linear noise attenuation of seismic data | |
| Liu | Noise reduction by vector median filtering | |
| Wu et al. | Removal of multisource noise in airborne electromagnetic data based on deep learning | |
| Fedi et al. | Understanding imaging methods for potential field data | |
| CN110133737B (zh) | 一种隐伏矿的电磁学预测方法 | |
| Ji et al. | Noise reduction of grounded electrical source airborne transient electromagnetic data using an exponential fitting-adaptive Kalman filter | |
| Bouchedda et al. | Sferics noise reduction in time-domain electromagnetic systems: application to MegaTEMII signal enhancement | |
| CN108983300B (zh) | 一种隧道掘进机施工条件下的瞬变电磁隧道超前预报方法 | |
| Lv et al. | Noise removal for semi-airborne data using wavelet threshold and singular value decomposition | |
| Xue et al. | Airborne electromagnetic data denoising based on dictionary learning | |
| Feng et al. | Resistivity-depth imaging with the airborne transient electromagnetic method based on an artificial neural network | |
| Li et al. | Using wavelet packet denoising and a regularized ELM algorithm based on the LOO approach for transient electromagnetic inversion | |
| Sun et al. | A self-constrained inversion of magnetic data based on correlation method | |
| CN104424393B (zh) | 一种基于主成分分析的地震数据储层反射特征加强方法 | |
| CN116090662B (zh) | 基于知识图谱的境外铜矿潜力预测方法、系统及电子设备 | |
| CN115016008A (zh) | 基于神经网络的电性源感应-极化共生效应多参数成像方法 | |
| Zhang et al. | Application of u-net for the recognition of regional features in geophysical inversion results | |
| Chukwu et al. | Unsupervised machine learning and depth clusters of Euler deconvolution of magnetic data: a new approach to imaging geological structures | |
| Zhang et al. | Development and future prospects of quantitative mineral assessment in China | |
| CN110865409B (zh) | 一种基于波阻抗低秩正则化的地震波阻抗反演方法 | |
| CN111366976A (zh) | 基于尺度导向的地震属性自适应中值滤波方法 | |
| Marchant et al. | 3D inversion of electromagnetic logging-while-drilling data | |
| CN110632664B (zh) | 基于多地球物理参数相关性的含油气预测方法与装置 | |
| Deng et al. | Magnetic Object Identification Based on Convolutional Neural Network |
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: 20200522 Termination date: 20210619 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |