大深度激电法中电磁耦合效应分析与去除方法
Analysis and removal of electromagnetic coupling effects in the large-depth induced polarization method
通讯作者: 陈儒军(1973-),副教授,博士,主要从事地球物理电法勘探仪器与信号处理工作。Email:chrujun12358@gmail.com
第一作者:
收稿日期: 2024-12-26 修回日期: 2025-03-25
| 基金资助: |
|
Received: 2024-12-26 Revised: 2025-03-25
随着浅层矿产资源的逐渐枯竭,深部矿产勘探成为矿业发展的必然趋势,增大勘探深度最直接、最有效的方法是增加激电装置的供电极距。然而,长供电极距会引发强烈的电磁耦合效应,从而对激电信号造成严重干扰。针对这一问题,本文基于解析方法计算了不同测量装置在均匀半空间和层状介质中的电磁耦合效应,通过对比,分析了测量装置类型、布线方式、供电极距、大地电阻率及频率等因素对电磁耦合强度的影响。最后,根据激电效应和电磁耦合效应在频率域内相位的差异,推导了相对相位谱的计算公式,并针对不同的应用场景对其去耦效果进行理论分析。研究结果表明,供电极距的增大、大地电阻率的降低以及工作频率的升高均会显著增强电磁耦合的干扰强度。在相同条件及探测深度下,四极装置受到的电磁耦合干扰通常高于三极装置。相较于传统的激电相位谱,相对相位谱的引入使三极装置和四极装置的最大工作频率分别提高了4倍和10.6倍,表明相对相位法在大深度激电法勘探中具有一定的去耦效果。本研究对大深度激电勘探的现场施工具有重要的指导意义。
关键词:
With the gradual depletion of shallow mineral resources, deep mineral exploration has emerged as an essential development trend in the mining industry. Increasing the distance between receiver electrodes of induced polarization (IP) devices is the most direct and effective approach to enhance the exploration depth. However, a long distance can cause strong electromagnetic (EM) coupling effects, severely interfering with IP signals. Addressing this challenge, this study calculated the EM coupling effects of various measuring devices in homogeneous half-space and layered media using analytical methods. Furthermore, this study comparatively analyzed the impacts of various factors, including measuring device type, wiring layout, distance between receiver electrodes, earth resistivity, and frequency, on the EM coupling intensity. Based on the phase differences between IP and EM coupling effects in the frequency domain, this study derived the calculation equation of the relative phase spectrum, followed by a theoretical analysis of the decoupling effects in application scenarios. The results indicate that increasing distance between receiver electrodes, decreasing earth resistivity, and raising working frequency all significantly intensified the EM coupling interference. Under consistent conditions and detection depths, the Schlumberger array suffered from higher EM coupling interference compared to the pole-dipole array. Compared to the traditional IP phase spectrum, the relative phase spectrum enhanced the maximal working frequency of the pole-dipole and Schlumberger arrays by four and 10.6 times, respectively, suggesting the decoupling capability of the relative phase method in large-depth IP exploration. Overall, this study provides significant guidance for the field implementation of large-depth IP exploration.
Keywords:
本文引用格式
冀振兴, 秦浩杰, 陈儒军, 王权功.
JI Zhen-Xing, QIN Hao-Jie, CHEN Ru-Jun, WANG Quan-Gong.
0 引言
上述问题可以从硬件采集端和数据处理两个方面来解决。前者可通过增大发射电流提高数据采集信噪比,再经选频接收、相干检测等手段增强接收机抗干扰能力[10
第一类是假设频率域激电效应的振幅随频率对数近似线性增加,而电磁耦合的振幅则随频率呈指数增加,从而可以通过两频校正方法去除电磁耦合效应[14-
上述研究成果对激电信号的高质量提取有着重要的意义,但他们的研究对象均为浅层矿体,没有考虑在深部勘探中电磁耦合对激电信号的影响。在大深度勘探中,由于供电极距的增加导致电磁耦合强度增大,测量信号也会变得更加复杂,上述方法假设的有效性需要进一步的验证。为此,本文基于Thomas等[22]发布的CR1Dmod开源程序对大深度激电中的电磁耦合效应进行研究,分析装置类型、布线方式、大地电阻率阻值和频率等因素对电磁耦合强度的影响。同时从相对相位谱的计算原理出发,分析其在处理大深度激电数据的适用性,并定量考察其在三极装置和四极装置中的去耦效果,得到的结论对大深度激电找矿的现场施工有一定的指导意义。
大深度激电法中的电磁耦合是利用激电法进行矿产勘探时,供电线路与测量线路之间的电性感应,主要包括电容耦合和电磁耦合。电容耦合是指供电线路与测量导线以及大地之间因存在分布电容而导致的假激电异常,Madden等[23]通过研究证明,只要供电导线和测量导线之间的距离足够大,电容耦合的影响就可以忽略不计。因此,本文仅考虑外在因素对电磁耦合强度的影响。
1 电磁耦合计算原理
电磁耦合指的是发射导线和接收导线间随频率变化的互感应阻抗。Sunde[24]在1968年指出地表任意两条接地导线之间的电磁耦合响应可以通过对ds和dS这两个偶极子的互阻抗沿各自的导线分别进行积分得到,并给出了其理论解,如式(1)所示:
式中:ds和dS分别是测量导线MN和供电导线AB上面的线元;r为两个线元之间的距离;θ是两个线元之间的夹角。为了简单起见,本文仅考虑供电线路和测量线路之间共线或平行的情况,即θ=0。
式(1)中的P(r)只与供电导线、测量导线和大地相互间的电磁耦合效应有关,因此称为感应函数,Q(r)只与供电电极和测量电极的相互位置有关,因此称为接地函数。Dey等[25]在1973年进一步给出了P(r)和Q(r)的表达式:
式中:u0=
在均匀半空间中,P(r)和Q(r)的表达式可以通过推导进一步简化为式(3)的形式[26]:
接地相Q(r)只与电极位置有关,其表达式可以表示为:
因此在均匀半空间地层中,供电线路和测量线路之间总阻抗的表达式为:
式中:K为装置系数,只与电极位置有关。对于层状介质中的电磁耦合,本文利用发布的CR1Dmod开源程序进行计算,其正确性已经在多篇文献中得到了验证[27-
常规的激电法勘探中常用测量阻抗的相位来表示测量信号,其计算方法如式(6)所示:
式中:Re和Im分别代表取阻抗的实部和虚部;ϕ为相位,主要有激电信号ϕIP和电磁耦合信号ϕEM两部分组成。
2 影响电磁耦合强度的因素分析
图1
图2
图2
不同电极装置中激电相位和电磁耦合相位
Fig.2
The phase of IP and EM coupling in different electrode arrays
图2展示的是不同装置参数下的纯激电相位谱,可以看到当频率极低或极高时相位较小,在中间的某个频段取得极大值,这是由不同频率下的探测深度差异引起的。当频率极低时,探测深度较大,测量信号只受第二层电阻率的影响(近似为均匀半空间),因此相位几乎为0。同样,当频率极高时,测量信号只受第一层电阻率的影响,相位值也相对较低。此外,图2中4条激电信号曲线完全重合,说明四极装置和三极装置的探测深度都只与AO的长度有关,从而验证了它们理论探测深度的统一性(即三极装置的探测深度d1≈AO/3,对称四级的装置的探测深度d2≈AB/6=AO/3)。从图2b中可以看出不同装置参数下的电磁耦合大小是不同的,四极装置受电磁耦合干扰程度明显大于三极装置,且三极装置中h越大,电磁耦合效应越小。因此在实际测量时,三极装置中的B电极应放置在尽可能远的位置。
为进一步测试电磁耦合效应随供电极距和频率的变化规律,计算了不同供电极距下的三极装置在均匀半空间中电磁耦合效应。由于在三极装置中,BO>5AO时即可认为电极B在无穷远的位置[29],因此在下述实验中,BO长度均被设置为15 000 m。激电法在大深度测量时,工作频率通常小于1 Hz,因此共设置了1/16 Hz、2/16 Hz、4/16 Hz、8/16 Hz和1 Hz 5个频点,具体计算结果如图3所示。其中图3是电阻率为100 Ω·m时的计算结果,图3b是电阻率为1 000 Ω·m时的计算结果。通过图3和3b对比发现,无论在低阻地层中还是在高阻地层中,电磁耦合效应随供电极距和频率的变化规律是完全一致的,即电磁耦合效应总是随着供电极距的增加而增大,随着频率的增高而变强。此外能发现在同样的供电极距和频率下,高阻层中的电磁耦合强度总是比低阻层更弱,说明电磁耦合效应随大地电阻率的升高而降低。同时也在斯伦贝谢装置中进行了同样的实验,得到的结论与三极装置中完全一致,此处不再赘述。
图3
图3
三极装置中电磁耦合效应随供电极距的变化规律
Fig.3
The variation of the EM coupling with the spacing of the current-electrodes in the pole-dipole array
在图1所示的旁侧中梯装置中,常通过改变旁侧距离来调整测量点位置和探测深度,因此有必要研究旁侧距离大小对电磁耦合效应的影响。图4展示了在不同大地电阻率和频率下,电磁耦合强度随旁侧距离的变化规律。其中,x和y的负值分别代表测量导线MN中点O1在供电导线中点O的左侧和下方,反之则代表O1在O的右侧和上方。从图中可以看出,在不同电阻率和频率下,电磁耦合强度的数值差别很大,但每一幅子图中电磁耦合强度随旁侧距离的变化规律完全一致,即旁侧距离越大,电磁耦合效应越小。此外,每幅子图均呈现中心对称性,这说明电磁耦合强度只与旁侧距离的大小有关,而与O及O1之间的相对位置无关。不同子图中电磁耦合强度的差异再次验证了上述实验结论,即频率越高,大地电阻率越小,电磁耦合强度越大。
图4
图4
中梯装置中旁侧距离对电磁耦合强度的影响
Fig.4
Impact of side distance on EM coupling in intermediate gradient array
3 相对相位谱在不同测量装置中的去耦效果分析
常规的激电相位谱作为激电法的测量信号,在大深度勘探时极易受到电磁耦合效应的影响。前人已经验证了相对相位谱在常规深度勘探中的良好去耦效果,本章将进一步从计算原理出发,分析其适用场景,并定量考察其在大深度激电勘探中的应用效果。
3.1 相对相位谱计算原理
相对相位谱的计算如式(7)所示:
式中:ϕr(f)和ϕ(f)分别是频率为f时的相位和相对相位;k为相对相位法中的频比,一般来说,k越大,激电相位谱和相对相位谱的差异越小[19]。
相位谱由激电信号和电磁耦合效应两部分组成,二者的频率特性有很大差异,在一定频率范围内近似满足式(8):
将式(8)代入式(7)可将相对相位谱的计算公式进一步推导为式(9),可以发现最终得到的相对相位ϕr与激电相位ϕIP近似相等,从而消除了电磁耦合效应ϕEM的影响。
从上述公式的推导可以发现,相对相位法的去耦效果主要取决于式(8)中“≈”两侧数值之间的差异。当两侧数值完全相等时,电磁耦合的影响能够完全消除;反之,二者之间的差异越大,相对相位谱计算过程中引入的误差也越大,从而影响去耦效果。为了定性的分析相对相位法的应用场景,计算了斯伦贝谢装置在不同电阻率下的频率特性,结果如图5所示。可以看到图5中由于供电极距较小,几乎所有曲线都呈现良好的线性关系,说明相对相位谱在常规勘探深度中具有良好的去耦效果。然而,随着供电极距(勘探深度)的增大,低阻和高频段的曲线形态不再满足式(8)所述的线性关系,且供电极距越大,电磁耦合效应与频率之间的线性关系越差,表明相对相位谱在长极距、高频和低阻条件下的去耦效果较差。
图5
图5
斯伦贝谢装置中电磁耦合响应频率特性
Fig.5
The frequency response characteristics of EM coupling in Schlumberger array
3.2 相对相位法在不同装置类型中的去耦效果
在长极距、低阻和高频情况下,相对相位法难以完全去除电磁耦合效应,此时式(8)中的ϕr≈ϕIP不再成立,而是也包含了一部分电磁耦合干扰,即ϕr=ϕIP+ϕEM。因此,可以通过计算不同装置类型、频率和大地电阻率下的ϕEM分量,从而定量分析相对相位法在大深度(主要针对1 000 m的勘探深度)激电法中的应用效果。
图6
图6
相对相位法中频比对去耦效果的影响
Fig.6
The impact of frequency ratio on decoupling effect in the relative phase method
为分析相对相位法在不同测量装置中的去耦效果,首先基于均匀半空间(Rt=350 Ω·m)地层计算了三极装置(AO=3 000 m)和四极装置(AB=6 000 m)常规激电相位谱中的电磁耦合分量,如图7中的两条实线所示。然后将频比设置为2,计算相对相位谱中的电磁耦合分量,如图7中的两条虚线所示。对同一测量装置而言,在低频段激电相位谱中的电磁耦合效应总是强于相对相位谱,说明相对相位法能够去除一部分耦合干扰。然而,随着频率的增高,激电相位谱和相对相位谱之间的差异逐渐减小,说明相对相位法的去耦效果随着频率的升高而变差,尤其是当频率增加到10 Hz时,相位和相对相位几乎相等,说明此时由于式(8)中“≈”两侧数值之间的差异过大,导致相对相位法不再具有去耦能力。此外,在同一频率下,三极装置测量的激电相位或相对相位中的电磁耦合分量均小于四极装置,再次验证了在相同的勘探深度下,三极装置受电磁耦合干扰程度更小。
图7
图7
不同测量装置的激电相位谱和相对相位谱对比
Fig.7
Comparison of IP phase and relative phase in different electrode arrays
针对1 000 m的探测深度,将三极装置的供电极距AO设置为3 000 m,斯伦贝谢装置的供电极距AB设置为6 000 m,分别计算它们在不同电阻率下的最大工作频率,结果如图8所示。可以看到,随着大地电阻率的增加,由于电磁耦合干扰强度的降低,常规的激电相位谱法和相对相位谱法的最大工作频率都在线性变大。利用激电相位谱测量时,三极装置的最大工作频率明显大于四极装置,而采用相对相位谱进行测量时,虽然三极装置的最大工作频率依然大于四极装置,但二者之间的差距变得更小,说明相对相位法在不同测量装置中的去耦效果是不同的。图8b进一步展示了相对相位谱的去耦效果,可以看到三极装置在采用相对相位谱后,最大工作频率提升了4倍,而四极装置则提升了近11倍,且最大工作频率提升的倍数均为常数,不随大地电阻率的变化而变化。通过将激电相位谱和相对相位谱的最大工作频率作对比,证明相对相位法在大深度激电法中也有一定的去耦效果。
图8
图8
激电相位谱和相对相位谱的最大工作频率对比
Fig.8
Comparison of the maximum working frequency of IP phase and relative phase
3.3 相对相位法去耦案例
为了验证相对相位法去除电磁耦合干扰的有效性,本文对阿拉善呼和沙拉某矿区的扩频激电测量数据进行了分析。数据采集时使用三极装置,供电电极(A)数量为66个,分布在-3 600~3 600 m的测线上,测量电极共有51对,均匀分布在-1 000~1 000 m的位置上。通过对所有测量数据进行处理,形成如图9所示的两个断面,分别为激电相位谱断面(图9)和相对相位谱断面(图9b)。从图9中可以看出常规的激电相位谱受到了严重的电磁耦合干扰,导致相位谱图中出现了严重的噪声和伪影,无法区分出激电异常信号,导致无法识别出潜在的异常区域。采用相对相位法后,图9b中的电磁耦合干扰显著减少,激电异常信号的识别变得更加清晰,在测线-700~-500 m、300~500 m的区域内出现了明显的激电异常。通过对比两幅图像可以发现,相对相位法能够有效去除电磁耦合的影响,显著提高了激电异常识别的准确性。因此,基于相对相位法的去耦技术在扩频激电测量中的应用具有重要的实际意义。
图9
图9
相对相位法去耦效果(以呼和沙拉某矿区为例)
Fig.9
The decoupling effect of the relative phase method (a case study of the Huhe-shala mining area)
4 结论
1)三极装置和对称四级装置的理论探测深度只与AO的长度有关,在AO长度相同时,两种测量装置测得的纯激电信号是完全一样的,而主要区别在于在相同的探测条件下,四极装置受电磁耦合干扰程度总是强于三极装置。
2)三极装置和四极装置中的电磁耦合强度主要受到工作频率、供电极距和大地电阻率的影响,其中大地电阻率越小、电极距越长、工作频率越高,电磁耦合强度越大。此外,旁侧中梯装置的电磁耦合强度还受到旁侧距离的影响,旁侧距离越小,电磁耦合强度越大,当垂直旁侧距离和水平旁侧距离均为0时,电磁耦合强度最大。三极装置中的电磁耦合强度还受到OB距离的影响,即B点放置的越远,电磁耦合效应越小,因此在实际施工时,电极B应放置的尽可能远。
3)由于相对相位法中假设的条件,即k·ϕEM(f)≈ϕEM(k·f)并不总是满足,特别是在长极距、高频、低阻情况下,因此其去耦效果及应用场景是有限的。针对1 000 m的探测深度,相对相位谱的引入可将对称四极装置和三极装置的最大工作频率分别提高10.6倍和4倍。
参考文献
金属矿产资源的深部找矿、勘探与成矿的深层动力过程
[J].
exploration and exploitation for metal mineral resocrces and its deep dynamical process of formation
[J].
深部金属矿产资源地球物理勘查与应用
[J].
Geophysical exploration and application of deep metallic ore resources
[J].
深部地质资源地球物理探测技术研究发展
[J].
Research progress of geophysical exploration technology for deep geological resources
[J].
Electrical-prospecting method for hydrocarbon search using the induced-polarization effect
[J].
Benefits of the induced polarization geoelectric method to hydrocarbon exploration
[J].
Investigation of a Bronze Age plankway by spectral induced polarization
[J].DOI:10.1002/arp.v9:4 URL
Application of time domain induced polarization to the mapping of lithotypes in a landfill site
[J].
DOI:10.5194/hess-16-1793-2012
URL
[本文引用: 1]
. A direct current (DC) resistivity and time domain induced polarization (TDIP) survey was undertaken at a decommissioned landfill site situated in Hørløkke, Denmark, for the purpose of mapping the waste deposits and to discriminate important geological units that control the hydrology of the surrounding area. It is known that both waste deposits and clay have clear signatures in TDIP data, making it possible to enhance the resolution of geological structures compared to DC surveys alone. Four DC/TDIP profiles were carried out crossing the landfill, and another seven profiles in the surroundings provide a sufficiently dense coverage of the entire area. The whole dataset was inverted using a 1-D laterally constrained inversion scheme, recently implemented for TDIP data, in order to use the entire decay curves for reconstructing the electrical parameters of the soil in terms of the Cole-Cole polarization model. Results show that it is possible to resolve both the geometry of the buried waste body and key geological structures. In particular, it was possible to find a silt/clay lens at depth that correlates with the flow direction of the pollution plume spreading out from the landfill and to map a shallow sandy layer rich in clay that likely has a strong influence on the hydrology of the site. This interpretation of the geophysical findings was constrained by borehole data, in terms of geology and gamma ray logging. The results of this study are important for the impact of the resolved geological units on the hydrology of the area, making it possible to construct more realistic scenarios of the variation of the pollution plume as a function of the climate change.
Electromagnetic coupling in frequency-domain induced polarization data:A method for removal
[J].DOI:10.1111/gji.2001.145.issue-1 URL [本文引用: 1]
Fundamental characteristics of an approximate correction method for electromagnetic coupling in frequency-domain induced polarization
[J].
DOI:10.1190/1.1441913
URL
[本文引用: 1]
Our studies have determined that both gradient and dipole-dipole arrays may produce a rather strong electromagnetic (EM) coupling effect under certain conditions. Variations in these apparent effects are closely correlated with changes of surface current density and also influenced by topographic relief. The dipole-dipole array can provide EM coupling anomalies with the same sign and shape as bona fide induced-polarization (IP) anomalies.It has been shown by theoretical calculations and model studies that response curves of frequency effects of EM coupling in logarithmic coordinates have distinctive features of sectional linearization. This can provide a basic method to correct approximately EM coupling. Compared with methods being used inside and outside China, the suggested correction method has the advantage of wide application. This has been already proven by a few field examples.
高精度多频电法数据采集系统
[J].
The data acquisition system in the high-precision multi frequency electric method
[J].
相位激电法在强干扰区的应用试验
[J].
Phase-IP experimental effects under the condition of strong interference and difficult grounding
[J].
频率域激电发送端去耦建模与仿真分析
[J].
Decoupling model and simulation analysis from signal transmitters in frequency domain IP
[J].
组合激电测深装置的应用与研究
[J].
The application and study of induced polarization group device
[J].
频率域激电法的电磁耦合校正应用研究
[J].
Applied research of electromagnetic coupling correction in frequency domain IP
[J].
变频法电磁耦合效应校正原理与实例
[J].
The correction principle and examples of electromagnetic coupling effect in frequency conversion method
[J].
多频去耦在双频激电数据去耦中的应用
[J].
The application of multiple frequency decoupling method to the removal of dual-frequency IP coupling effect
[J].
中梯装置相位激电电磁耦合研究
[J].
A study on electromagnetic coupling of central gradient array phase-IP
[J].
相位激电电磁耦合两频校正技术
[J].
The phase-IP electromagnetic coupling correction technique of two frequencies
[J].
多频激电相对相位谱研究
[J].
The study of relative phase spectrum in multi-frequency induced polarization
[J].
相位激电法电磁耦合效应的校正
[J].
The Phase-IP electromagnetic coupling correction
[J].
利用指数函数拟合去除扩频激电信号中的电磁耦合感应
[J].
A new method for removing EM-coupling in spread spectrum induced polarization based on exponential function fitting
[J].
CRID mod:A matlab program to model ID complex resisivity effects in electrical and electromagnetic surveys
[J].DOI:10.1016/j.cageo.2006.01.001 URL [本文引用: 1]
Induced polarization,a review
[C]//
Electromagnetic coupling in frequency and time-domain induced-polarization surveys over a multilayered earth
[J].
DOI:10.1190/1.1440348
URL
[本文引用: 1]
Electromagnetic coupling responses in frequency and time-domain induced-polarization measurements over a multilayered earth are evaluated. For collinear dipole-dipole and pole-dipole configurations over a dissipative layered subsurface, the percent frequency effects of electromagnetic coupling are seen to be as high as 60 percent for large L2(σ1f) values, where L is the length of the receiving dipole, σ1 is the conductivity of the top layer of the half-space, and f is the higher frequency of excitation used. In both frequency and time-domain analyses, the distinctive effects of layering compared to that of a homogeneous half-pace response are shown for different electrode configurations, layer geometry, and electrical parameters of the subsurface. The pole-dipole configuration of electrodes, in general, exhibits higher coupling compared to the dipole-dipole configuration. In time-domain measurements, the late off-time transient decays reflect almost entirely the normal polarizability of the layered subsurface, in that the coupling responses are significant only during the early off-time of the transient.
CR1Dinv:A Matlab program to invert 1D spectral induced polarization data for the Cole-Cole model including electromagnetic effects
[J].DOI:10.1016/j.cageo.2008.06.001 URL
/
| 〈 |
|
〉 |
