CN115128680B - 一种磁性源多波形组合的瞬变电磁靶向测量方法 - Google Patents
一种磁性源多波形组合的瞬变电磁靶向测量方法 Download PDFInfo
- Publication number
- CN115128680B CN115128680B CN202210752160.6A CN202210752160A CN115128680B CN 115128680 B CN115128680 B CN 115128680B CN 202210752160 A CN202210752160 A CN 202210752160A CN 115128680 B CN115128680 B CN 115128680B
- Authority
- CN
- China
- Prior art keywords
- waveform
- current
- abnormal
- field
- equation
- 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
Links
Images
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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
-
- 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
-
- 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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
- G01V2003/084—Sources
-
- 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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
- G01V2003/085—Receivers
-
- 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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
- G01V2003/086—Processing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种磁性源多波形组合的瞬变电磁靶向测量方法,目的在于对梯形波、三角波、半正弦波及其组合电流波形的探测效果进行分析,进行靶向测量来实现瞬变电磁精细化探测。采用CN‑FDTD方法对电磁波动方程进行离散,将发射电流转化为电流密度,加载到电场控制方程进行含源计算;加源完成后,将电场转化为感应电动势,用磁场控制方程进行无源计算。首先对电流波形频谱进行分析,再采用均匀半空间模型、层状异常体模型和复杂异常体模型对不同电流波形的响应和累积灵敏度分布进行分析,建立靶向激励参数集;根据目标体特征进行精准化靶向测量,对纵向分布多异常体和低阻层下方异常体进行有效测量,解决瞬变电磁法对复杂异常体探测精度低的问题。
Description
技术领域
本发明涉及一种磁性源多波形组合的瞬变电磁靶向测量方法,适用于电磁法地球物理勘探或复杂地质结构精细化探测领域,尤其是为磁性源多波形组合靶向测量进行工作参数优化。
背景技术
磁性源时域电磁探测方法,利用发射机控制发射线圈向地下发射双极性脉冲电流,脉冲电流可以为梯形波、三角波和半正弦波等,在发射电流发射过程及完全关断后,通过接收线圈或超导量子干涉仪(SQUID,Superconducting Quantum Interference Device)与接收机采集一次磁场和二次感应场随时间变化的衰减曲线,获得感应电动势和磁场信息,对数据进行处理解释获得地下介质丰富的电性信息。时域电磁法具有探测深度大、经济便捷的优点,在地下资源勘探、工程地质勘探、城市地下空间浅层探测、地下水资源与环境保护等领域得到了广泛应用。
传统的时域电磁法由于发射装置的限制仅能对梯形波进行激发,且由于发射线圈为感性负载,发射电流关断后呈e指数衰减,在电流尾部还会有过冲和震荡现象,这些原因造成发射电流波形质量较低;传统的接收装置是采用发射线圈获得感应电动势信息,在发射电流关断之后进行采集,测量参数单一;发射系统性能低下和解释参数单一严重限制了瞬变电磁法的进一步发展和应用。
为满足多深度分布复杂异常体精细化探测的需要,传统的双极性梯形波不再适用,双极性三角波脉宽较小,高频成分丰富,其趋肤深度较小,能够对浅层目标体进行有效识别;双极性半正弦波脉宽较大,已有研究成果表明,在发射过程中进行测量,可有效获得近地表地质体信息;采用组合电流波形可对不同深度的多个异常体、异常体上方含有低阻层的异常体进行有效识别,提高测量效率和测量精度,减小电磁法探测盲区。
电磁波动方程是由麦克斯韦基本方程组导出,将电场和磁场分离开来,可以将发射电流波形转化为电流面密度加载到x、y分量的电场控制方程,进行含激励源计算;加源过程完成后,将电场转化为感应电动势,采用磁场垂直分量控制方程继续进行迭代;与传统的麦克斯韦控制方程相比,避免了电场和磁场的交叉迭代,且能够进行含激励源计算,满足了多波形组合发射波形数值模拟和响应分析的需求。
三维时域有限差分(Three-dimensional Finite-difference Time-domain)方法是对时间或空间偏导项进行中心差分近似,通过在时域的迭代模拟电磁波的传播过程,但是时间步长的选取受到数值稳定性(CFL)条件的限制,时间步长选取较小,迭代速度慢且稳定性较差;CN-FDTD离散算法是将当前时刻的空间偏导项表示为前一时刻和后一时刻空间偏导项和的形式,这样能提高迭代方程的稳定性,增大时间步长,提高计算速度;采用矩阵方程形式对控制方程进行求解即可获得空间中电磁场的响应。
为了表征电磁场能量在不同深度、不同时刻的分布情况,参考Adam Smiarowski和Greg Hodges的研究成果,使用地下空间中感应电动势的分布来计算系统的灵敏度,观测并记录电流发射过程中及电流关断后灵敏度的变化趋势。由于感兴趣的是灵敏度的深度分布,而不是侧向分布,因此对于每个深度,我们首先对该深度所有水平位置网格节点上的感应电动势进行求和,以给出每个深度水平的相对灵敏度。累积灵敏度从地面向下计算,对当前深度上方的相对灵敏度求和,再与该时刻的总响应进行归一化,得出累积灵敏度随时间变化的二维平面图;为便于观察,观测累积灵敏度为50%、60%、70%、80%和90%的分布曲线,累积灵敏度大于90%的区域,能量已经比较微弱,地面的测量结果很难反映出这些区域的地质特征,因此我们主要对不同发射电流波形的响应平面图特征和累积灵敏度为90%时的分布曲线进行对比分析,对不同发射电流波形的探测能力进行评价。
孙怀凤等用CN-FDTD方法对麦克斯韦方程组进行离散,增大了时间步长,减小了迭代的次数,并采用均匀半空间模型和层状模型对迭代精度进行了验证,发现该方法计算速度快且相对误差明显降低。Adam Smiarowski和Greg Hodges对半正弦波的响应进行了分析,且引入了电流累积灵敏度表示电磁波能量的分布,结果表明,在半正弦波下降过程中对浅地表地质体的灵敏度最高,可对地下30m左右的目标体进行有效识别。
中国专利CN112285788B公开了一种基于电磁波动方程的CPML边界条件加载方法,通过采用电磁波动方程作为磁场垂直分量的控制方程,显著提高了计算速度和精度,证明了用有限差分算法离散电磁波动方程的有效性。
中国专利CN112698410B公开了电性源双相导电介质感应-极化共生时域电磁探测方法,通过分析感应-极化共生效应特征,构建双可控沿梯形波发射靶向激励关系,快关断用于测量感应场,慢关断用于测量极化场,证明了通过优化激励源参数进行靶向激励对于提高测量精度的有效性。
以上所述方法公布了CN-FDTD方法和电磁波动方程的应用,以及电性源梯形波靶向激励的测量方法。但对于多波形组合发射电流的数值模拟及响应分析,已有数值模拟方法几乎都无法满足需求,无法通过响应分析优化靶向激励参数;绝大多数瞬变电磁测量装置均是发射梯形波,在关断之后采集感应场的衰减曲线,在面对复杂地质体时无法对地下结构进行高精度识别及精细化探测;如何在时域电磁法探测中建立靶向激励参数集、发射靶向激励电流波形、连续低噪声磁场和感应电动势采集,进行多参数数据解释实现对复杂异常体的高精度探测是本领域技术人员迫切解决的一个技术问题。
发明内容
本发明所要解决的技术问题在于提供一种磁性源多波形组合的瞬变电磁靶向测量方法,目的在于对梯形波、三角波、半正弦波及其组合电流波形的探测效果进行分析,进行靶向测量来提高瞬变电磁探测的灵敏度和精度。采用CN-FDTD方法对电磁波动方程进行离散,将发射电流转化为电流面密度,加载到电场控制方程进行含源计算;加源完成后,将电场转化为感应电动势,用磁场控制方程进行无源计算。首先对电流波形频谱进行分析,再采用均匀半空间模型、层状异常体模型和复杂异常体模型对不同电流波形的响应和累积灵敏度分布进行分析,建立靶向激励参数集;根据目标体特征进行精准化靶向测量,对纵向分布多异常体和含低阻层异常体进行有效测量,解决瞬变电磁法对复杂异常体探测精度低的问题。
本发明是这样实现的,一种多波形组合的瞬变电磁靶向测量方法:
1)、采用CN-FDTD方法对电磁波动方程进行离散,克服数值稳定性条件的限制,能够提高电磁场控制方程的稳定性,增大迭代的时间步长,减小内存占用,提高计算速度;
2)、对网格进行全空间不均匀剖分,包含空气层和地下层,将不同波形发射电流激励源转化为电流面密度,加载到电磁波动方程的电场分量控制方程中,进行含激励源迭代,发射电流关断后,将电场转化为感应电动势,采用磁场垂直分量控制方程继续进行迭代,进一步提高计算速度;
3)、采用矩阵方程求解电场或磁场值,获得不同激励源的地下异常体响应及累积灵敏度分布曲线;
4)、根据不同发射电流波形参数对复杂异常体模型的识别效果,确定最佳的单电流波形参数或组合电流波形参数,提高激励源对地下异常体的敏感度,建立靶向激励参数集;
5)、通过应用大功率IGBT组成的H桥路、无源钳位技术和RLC串联谐振技术设计参数可控的梯形波、三角波、半正弦波及其组合波形发射装置,同时采用24位AD转换、差分输入、低通滤波技术设计低噪声、连续和高采样率接收装置,构成磁性源多波形组合的瞬变电磁靶向测量系统;
6)、根据已有地质资料或者期望测量的地层深度及特征选择最佳的单波形或组合波形靶向激励参数,在发射过程中用电流采集装置记录发射波形,采用接收线圈或超导量子干涉仪(SQUID)记录感应电动势或磁场随时间变化的总场;
7)、根据记录的发射电流波形获得一次场信息,用总场减去一次场获得纯二次场响应,对数据进行预处理,再进行多参数联合成像,对地下异常体信息进行解释;
其中步骤1)中,电磁波动方程电场x、y分量和磁场垂直分量表达式为:
FDTD离散方法是将时间偏导项和空间偏导项(以x方向为例)用差分格式表示,离散格式为:
其中Δtn+1=tn+1-tn,Δtn=tn-tn-1;Δt为迭代时间步长,且Δtn+1=Δtn=Δt;Δxi和Δxi+1为沿x方向的网格边长,i、j、k分别为不同方向的网格数。
CN-FDTD离散方法是将当前时刻空间偏导项表示为前一时刻和后一时刻空间偏导项的和,能够显著提高数值模拟稳定性,CN-FDTD离散格式为:
式(5)中的第一项和第二项为有源迭代阶段,第三项为无源迭代阶段,g1、g2分别为不同时刻拉普拉斯展开项的系数,一般情况下取g1=g2=0.5,m是有源阶段的迭代次数,n是无源阶段的迭代次数;
步骤2)中,将发射电流转化为电流面密度,进行含源阶段迭代的电场控制方程为:
其中ΔI/Δt为发射电流随时间的变化率,Δt为迭代时间步长,lx、ly、lz为发射线圈所在位置的网格边长,Jx、Jy分别为x、y方向的电流面密度;
无源迭代阶段将电场转化为感应电动势Hmz,继续进行迭代,控制方程为:
步骤3)中,对网格进行剖分并对节点进行编号,网格节点数量为P=(i×j×k),求解电场或磁场的矩阵方程形式为:
其中矩阵A为方阵,其维度为(i×j×k)×(i×j×k),p为网格节点编号,矩阵中元素a,b通过公式(6)和(7)获得,求解磁场时a,b的表达式如公式(9)所示:
为了表征电磁场能量在不同深度、不同时刻的分布情况,对于每个深度,首先对该深度所有水平位置网格节点上的感应电动势进行求和,以给出每个深度水平的相对灵敏度;累积灵敏度从地面向下计算,对当前深度上方的相对灵敏度求和,再与该时刻的总响应进行归一化,累积灵敏度S(zw,tn)在数学上的表达式为:
其中U、V、W分别为x,y,z方向的网格剖分节点总数,w为当前深度所在的z方向的网格剖分节点数,tn为当前迭代次数n的时刻;
步骤4)中,首先采用三角波、梯形波和半正弦波对不同深度的单个异常体进行响应分析,观察二维平面图、三维切片图、累积灵敏度在异常体所在深度的分布评价对异常体识别的效果,在脉冲宽度大于等于2ms的半正弦波下降过程中对浅地表异常体识别效果更好,在脉冲时间小于等于1ms的三角波关断之后对较浅层异常体识别效果更好,在关断时间小于等于1ms的梯形波关断之后对深层异常体识别效果更好;
接下来采用单电流波形及组合电流波形对不同深度的双异常体进行响应分析,结果表明,三角-半正弦组合电流波形对浅地表-较浅层双异常体识别效果更好,梯形-半正弦组合电流波形对浅地表-深层双异常体识别效果更好,梯形-三角组合电流波形对较浅层-深层双异常体识别效果更好;
采用组合电流波形对异常体上方含有低阻层的模型进行响应分析,结果表明,先用半正弦波或三角波对浅地表或较浅层层状低阻体进行识别,再用三角波或梯形波的总响应减去层状异常体响应,即可对层状低阻体下方的较浅层低阻体或深层低阻体进行识别;
根据对多种类型、多种参数异常体识别的效果,针对异常体特征,确定最佳的单电流波形参数或组合电流波形参数,提高激励源对地下异常体的敏感度,建立靶向激励参数集,作为多波形发射装置的发射参数选择依据。
步骤5)中,接收装置为了适应多波形发射的需求,采用24位AD转换、差分输入、低通滤波技术实现低噪声、连续和高采样率电压采集;多波形发射和接收装置构成了磁性源多波形组合的瞬变电磁靶向测量系统;
步骤6)中,根据已有的粗略测量地质资料或者期望测量的地层深度及特征确定靶向激励参数,用霍尔传感器组成的电流采集装置记录发射波形信息;接收线圈抗干扰能力强,有效面积大,记录感应电动势随时间变化的总场,SQUID分辨率很高,有效信号持续时间长,记录磁场随时间变化的总场,接收线圈与SQUID配合使用,能够进一步提高测量的准确度二
步骤7)中,根据记录的发射电流波形信息获得一次场信息,用总场减去与一次场同相位的信号获得纯二次场响应;对二次场进行基线校正、叠加、滤波预处理,对感应电动势和磁场信号进行多参数联合成像,解释地下异常体信息,完成磁性源多波形组合的瞬变电磁靶向测量,有效提高探测灵敏度和精度。
本发明与现有技术相比,带来的有益效果在于:采用CN-FDTD方法离散电磁波动方程,同时进行全空间剖分和含激励源迭代,能够有效提高数值模拟精度,同时满足了多波形组合发射电流数值模拟的需求;通过对多种类型、多种参数异常体模型的响应和累积灵敏度分布曲线进行分析,能够优化激励源工作参数,建立靶向激励参数集,有效提高探测灵敏度;磁性源多波形组合的瞬变电磁靶向测量系统能够对地下目标体进行靶向激励,同时进行连续、低噪声、高分辨率感应电动势和磁场采集,最大限度减小探测盲区并提高数据的准确性和有效性;剔除一次场后,对二次场纯响应进行预处理和多参数联合成像,完成对地下异常体的解释和识别,促进时域电磁法精细化探测的发展进程。
附图说明
图1是磁性源多波形组合的瞬变电磁靶向测量方法流程图;
图2是本发明瞬变电磁多波形靶向激励系统的整体框图;
图3是本发明连续、低噪声、高分辨率采集装置的整体框图;
图4是三种单电流波形在均匀半空间中的累积灵敏度分布图;
图5是三种单电流波形在较浅层高阻层状体中的累积灵敏度分布图;
图6是电流波形示意图及频谱分布图;
图7是三角-半正弦组合波识别浅层高阻体-较浅层低阻体切片图;
图8是梯形-三角组合波识别较浅层低阻体-深层低阻体切片图;
图9是梯形-半正弦组合波识别较浅层层状低阻体-深层低阻体切片图;
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例
参见图1,一种磁性源多波形组合的瞬变电磁靶向测量方法,包括:
1)、根据CN-FDTD算法离散电磁波动方程,推导电场和磁场迭代方程,计算区域包含空气层和地下层,对网格进行不均匀剖分,整个电磁场计算区域尺寸为6780m×6780m×5340m,最小网格步长为10m,最大网格设置为240m,边界区域设置6层CFS-PML吸收边界,空气层z方向距离为1950m,用于观测的地下中心区域尺寸为3900m×3900m×1950m,网格节点设置为39×39×20个,其中x、y方向上的网格数均为39个,z方向上网格数20个。
2)、在整个计算区域内设置电导率、磁导率、人工介电常数、CPML系数、发射电流波形节点值、异常体位置及电导率等参数,空气层电导率设置为σair=1×10-6S/m,大地电导率设置为σ0=0.01S/m,磁导率设置为真空磁导率。
3)、建立大型稀疏矩阵,构建矩阵各参数,将发射电流转化为电流面密度,进行含源循环迭代,每迭代一次完成后,将电场转化为感应电动势保存下来;
4)、判断含源阶段是否迭代完成,当未完成时,重复步骤3;当含源阶段迭代完成后,将最后两个时刻保存的感应电动势作为无源阶段磁场迭代的初值,继续进行迭代,保存每个时刻的各节点感应电动势的值;
5)、判断无源阶段是否迭代完成,当未完成时,重复步骤4;当无源阶段迭代完成后,停止迭代,获得一个四维矩阵,包含迭代过程各个时刻各个节点的感应电动势。
6)、对结果进行响应分析,获取累积灵敏度随时间的分布图、感应电动势分布平面图和纯响应切片图,根据当前电流波形对异常体的识别效果进行综合评价;改变异常体模型,改变发射电流波形参数,针对异常体模型优化激励参数,建立靶向激励参数集。
7)、首先对电导率为σ0=0.01S/m的均匀半空间模型进行响应分析,三种电流波形累积灵敏度分布曲线如图4所示,在电流波形发射阶段,主要能量分布在浅地表,完全关断之后电磁波迅速扩散到200m深度,浅地表信息丢失;从累积灵敏度分布曲线可以看出,半正弦波下降阶段适合测量浅地表信息,三角波适合测量较浅层信息,梯形波适合测量深层信息。为了进一步确定不同波形对不同深度异常体的灵敏度差异,在距地面90m-270m深度设置层状高阻体,电导率为σl=0.001S/m,如图5所示,三种电流波形的90%累积灵敏度曲线在该深度的分布时间分别为0.2ms、0.3ms和0.4ms,三角波对较浅层高阻体更敏感。
8)、为了分析组合电流波形对电磁探测的有效性,对比了单电流波形及组合电流波形的频谱,如图6所示,组合电流波形频谱分布更加广泛;随后设置了复杂异常体模型,三角-半正弦组合波对浅层高阻-较浅层低阻异常体模型的纯响应切片图如图7所示,在半正弦波下降阶段(12.05ms-12.50ms),能够识别浅层高阻目标体,而在三角波关断很短时间内(0.205ms-0.24ms),浅层高阻信息已经丢失;在三角波及半正弦波关断后0.1ms和0.3ms时刻,三角波对较浅层低阻体的分辨能力更高;与之相似,梯形-三角组合波对较浅层低阻-深层低阻异常体模型的纯响应切片图如图8所示,梯形波(3.50ms)对深层低阻异常体识别更清晰,三角波(10.4ms)对较浅层低阻异常体识别更清晰。为了对实际地形中低阻层状体下方的异常体进行探测,采用梯形-半正弦组合波对浅层层状低阻-深层低阻异常体模型进行激励,其纯响应切片图如图9所示,在半正弦波下降沿阶段(11ms-12ms)可以识别浅层层状低阻体,梯形波(5.50ms)可以识别层状低阻体下方的低阻体。继续对其他模型进行响应分析,确定最佳激励参数,建立多波形靶向激励集。
9)、根据已有地质资料或期望探测范围确定靶向激励参数,采用多波形组合的瞬变电磁测量系统进行激励和采集,获取感应电动势和磁场信息二随后对数据进行预处理和多参数联合成像,对结果进行解释,获得地下异常体信息,实现瞬变电磁精细化探测。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种磁性源多波形组合的瞬变电磁靶向测量方法,其特征在于,包括如下步骤:
1)、采用CN-FDTD方法对电磁波动方程进行离散,克服数值稳定性条件的限制,能够提高电磁场控制方程的稳定性,增大迭代的时间步长,减小内存占用,提高计算速度;
2)、对网格进行全空间不均匀剖分,包含空气层和地下层,将不同波形发射电流激励源转化为电流面密度,加载到电磁波动方程的电场分量控制方程中,进行含激励源迭代,发射电流关断后,将电场转化为感应电动势,采用磁场垂直分量控制方程继续进行迭代,进一步提高计算速度;
3)、采用矩阵方程求解电场或磁场值,获得不同激励源的地下异常体响应及累积灵敏度分布曲线;
4)、根据不同发射电流波形参数对复杂异常体模型的识别效果,确定最佳的单电流波形参数或组合电流波形参数,提高激励源对地下异常体的敏感度,建立靶向激励参数集;
5)、通过应用大功率IGBT组成的H桥路、无源钳位技术和RLC串联谐振技术设计参数可控的梯形波、三角波、半正弦波及其组合波形发射装置,同时采用24位AD转换、差分输入、低通滤波技术设计低噪声、连续和高采样率接收装置,构成磁性源多波形组合的瞬变电磁靶向测量系统;
6)、根据已有地质资料或者期望测量的地层深度及特征选择最佳的单波形或组合波形靶向激励参数,在发射过程中用电流采集装置记录发射波形,采用接收线圈或超导量子干涉仪(SQUID)记录感应电动势或磁场随时间变化的总场;
7)、根据记录的发射电流波形获得一次场信息,用总场减去一次场获得纯二次场响应,对数据进行预处理,再进行多参数联合成像,对地下异常体信息进行解释。
2.根据权利要求1所述的一种磁性源多波形组合的瞬变电磁靶向测量方法,其特征在于:
步骤1)中,电磁波动方程电场x、y分量和磁场垂直分量表达式为:
CN-FDTD离散方法是将当前时刻空间偏导项表示为前一时刻和后一时刻空间偏导项的和,能够显著提高数值模拟稳定性,CN-FDTD离散格式为:
式(2)中的第一项和第二项为有源迭代阶段,第三项为无源迭代阶段,g1、g2分别为不同时刻拉普拉斯展开项的系数,一般情况下取g1=g2=0.5,m是有源阶段的迭代次数,n是无源阶段的迭代次数;
步骤2)中,将发射电流转化为电流面密度,进行含源阶段迭代的电场控制方程为:
其中ΔI/Δt为发射电流随时间的变化率,Δt为迭代时间步长,lx、ly、lz为发射线圈所在位置的网格边长,Jx、Jy分别为x、y方向的电流面密度;
无源迭代阶段将电场转化为感应电动势Hmz,继续进行迭代,控制方程为:
步骤3)中,对网格进行剖分并对节点进行编号,网格节点数量为P=(i×j×k),求解电场或磁场的矩阵方程形式为:
其中矩阵A为方阵,其维度为(i×j×k)×(i×j×k),p为网格节点编号,i、j、k分别为不同方向的网格数,矩阵中元素a,b通过公式(3)和(4)获得,求解磁场时a,b的表达式如公式(6)所示:
为了表征电磁场能量在不同深度、不同时刻的分布情况,对于每个深度,首先对该深度所有水平位置网格节点上的感应电动势进行求和,以给出每个深度水平的相对灵敏度;累积灵敏度从地面向下计算,对当前深度上方的相对灵敏度求和,再与该时刻的总响应进行归一化,累积灵敏度S(zw,tn)在数学上的表达式为:
其中U、V、W分别为x,y,z方向的网格剖分节点总数,w为当前深度所在的z方向的网格剖分节点数,tn为当前迭代次数n的时刻;
步骤4)中,首先采用三角波、梯形波和半正弦波对不同深度的单个异常体进行响应分析,观察二维平面图、三维切片图、累积灵敏度在异常体所在深度的分布评价对异常体识别的效果,在脉冲宽度大于等于2ms的半正弦波下降过程中对浅地表异常体识别效果更好,在脉冲时间小于等于1ms的三角波关断之后对较浅层异常体识别效果更好,在关断时间小于等于1ms的梯形波关断之后对深层异常体识别效果更好;
接下来采用单电流波形及组合电流波形对不同深度的双异常体进行响应分析,结果表明,三角-半正弦组合电流波形对浅地表-较浅层双异常体识别效果更好,梯形-半正弦组合电流波形对浅地表-深层双异常体识别效果更好,梯形-三角组合电流波形对较浅层-深层双异常体识别效果更好;
采用组合电流波形对异常体上方含有低阻层的模型进行响应分析,结果表明,先用半正弦波或三角波对浅地表或较浅层层状低阻体进行识别,再用三角波或梯形波的总响应减去层状异常体响应,即可对层状低阻体下方的较浅层低阻体或深层低阻体进行高精度识别;
根据对多种类型、多种参数异常体识别的效果,针对异常体特征,确定最佳的单电流波形参数或组合电流波形参数,提高激励源对地下异常体的敏感度,建立靶向激励参数集,作为多波形发射装置的发射参数选择依据。
3.根据权利要求1所述的一种磁性源多波形组合的瞬变电磁靶向测量方法,其特征在于:
步骤5)中,接收装置为了适应多波形发射的需求,采用24位AD转换、差分输入、低通滤波技术实现低噪声、连续和高采样率电压采集;多波形发射和接收装置构成了磁性源多波形组合的瞬变电磁靶向测量系统;
步骤6)中,根据已有的粗略测量地质资料或者期望测量的地层深度及特征确定靶向激励参数,用霍尔传感器组成的电流采集装置记录发射波形信息;接收线圈抗干扰能力强,有效面积大,记录感应电动势随时间变化的总场,SQUID分辨率很高,有效信号持续时间长,记录磁场随时间变化的总场,接收线圈与SQUID配合使用,能够进一步提高测量的准确度;
步骤7)中,根据记录的发射电流波形信息获得一次场信息,用总场减去与一次场同相位的信号获得纯二次场响应;对二次场进行基线校正、叠加、滤波预处理,对感应电动势和磁场信号进行多参数联合成像,解释地下异常体信息,完成磁性源多波形组合的瞬变电磁靶向测量,有效提高探测灵敏度和精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210752160.6A CN115128680B (zh) | 2022-06-28 | 2022-06-28 | 一种磁性源多波形组合的瞬变电磁靶向测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210752160.6A CN115128680B (zh) | 2022-06-28 | 2022-06-28 | 一种磁性源多波形组合的瞬变电磁靶向测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115128680A CN115128680A (zh) | 2022-09-30 |
CN115128680B true CN115128680B (zh) | 2023-06-02 |
Family
ID=83379486
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210752160.6A Active CN115128680B (zh) | 2022-06-28 | 2022-06-28 | 一种磁性源多波形组合的瞬变电磁靶向测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115128680B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116626772B (zh) * | 2023-07-25 | 2023-11-28 | 青岛海月辉科技有限公司 | 磁指纹特征分析方法及金属物质主动式探测方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103576205A (zh) * | 2013-11-17 | 2014-02-12 | 成都理工大学 | 基于组合磁性源技术的地空瞬变电磁勘查方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2794532B1 (fr) * | 1999-06-02 | 2001-08-03 | Commissariat Energie Atomique | Procede de detection electromagnetique d'objets conducteurs |
RU2301431C2 (ru) * | 2005-03-24 | 2007-06-20 | Общество с ограниченной ответственностью "Сибирская геофизическая научно-производственная компания" | Способ электроразведки с использованием пространственного дифференцирования поля становления на нескольких разносах |
US10126456B2 (en) * | 2015-01-07 | 2018-11-13 | The Regents Of The University Of California | System and method for groundwater detection and evaluation |
CN105158808B (zh) * | 2015-07-01 | 2017-10-31 | 长安大学 | 一种浅海瞬变电磁海空探测及其解释方法 |
CN107065019A (zh) * | 2016-12-07 | 2017-08-18 | 北京桔灯地球物理勘探股份有限公司 | 应用于道路灾害和塌陷检测的三维电磁成像装置及使用方法 |
CN106772644B (zh) * | 2016-12-14 | 2018-06-26 | 中国矿业大学 | 矿井瞬变电磁三分量探测方法 |
CA3122828C (en) * | 2020-12-16 | 2023-10-31 | Jilin University | Squid-based electromagnetic detection method for induction-polarization symbiotic effect of two-phase coducting medium |
CN112698410B (zh) * | 2020-12-16 | 2022-03-15 | 吉林大学 | 电性源双相导电介质感应-极化共生时域电磁探测方法 |
CN113158517B (zh) * | 2021-03-31 | 2022-08-09 | 深圳大学 | 基于系统矩阵结合的cn-fdtd仿真方法、装置及相关组件 |
CN113866835B (zh) * | 2021-11-10 | 2023-07-14 | 吉林大学 | 一种时域三波形组合的电磁发射系统及控制方法 |
CN114114434B (zh) * | 2021-11-23 | 2024-06-14 | 重庆璀陆探测技术有限公司 | 瞬变电磁发射机组合式双脉冲发射电路的控制方法 |
-
2022
- 2022-06-28 CN CN202210752160.6A patent/CN115128680B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103576205A (zh) * | 2013-11-17 | 2014-02-12 | 成都理工大学 | 基于组合磁性源技术的地空瞬变电磁勘查方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115128680A (zh) | 2022-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Song et al. | Application of particle swarm optimization to interpret Rayleigh wave dispersion curves | |
CN105044793B (zh) | 一种多道瞬变电磁探测数据的反演方法和装置 | |
CN103955003B (zh) | 一种在超导瞬变电磁应用中的噪声抑制方法 | |
CN105785455B (zh) | 一种基于b样条插值的二维地面核磁共振反演方法 | |
CN112068212A (zh) | 一种无人机半航空时间域电磁探测数据分析解释方法 | |
CN112949134B (zh) | 基于非结构有限元方法的地-井瞬变电磁反演方法 | |
CN110133733B (zh) | 一种基于粒子群优化算法的电导-极化率多参数成像方法 | |
CN115128680B (zh) | 一种磁性源多波形组合的瞬变电磁靶向测量方法 | |
CN110348568B (zh) | 一种适用于强电磁干扰地区的深部采空区探测方法 | |
US9335434B2 (en) | Method and system for processing data acquired in an electromagnetic survey | |
Zhou et al. | Migration velocity analysis and prestack migration of common-transmitter GPR data | |
Christiansen et al. | A quantitative appraisal of airborne and ground-based transient electromagnetic (TEM) measurements in Denmark | |
CN114779355A (zh) | 基于发射电流全波形的地面瞬变电磁法反演方法及装置 | |
Gupta et al. | Straightforward inversion of vertical electrical sounding data | |
CN108776357B (zh) | 沉积地层瞬变电磁法电磁干扰的校正方法及装置 | |
CN112287545B (zh) | 一种双相导电介质的时-空分数阶电导率建模及模拟方法 | |
Li et al. | One-dimensional full-waveform inversion for magnetic induction data in ground-based transient electromagnetic methods | |
CN114721059B (zh) | 一种基于瞬变电磁法的地下成像方法 | |
CN116973983A (zh) | 一种磁性源多波形组合精细化探测的数据处理方法 | |
CN115906559B (zh) | 一种基于混合网格的大地电磁自适应有限元正演方法 | |
CN115563791B (zh) | 基于压缩感知重构的大地电磁数据反演方法 | |
Wang et al. | Characterization of excitation effects and data interpretation of combined time-domain multiwaveform transmission currents | |
Baruah et al. | Ground motion parameters of Shillong plateau: one of the most seismically active zones of northeastern India | |
CN111965714B (zh) | 一种基于暂态过程的电磁探测方法、设备及存储设备 | |
CN114675337A (zh) | 一种基于多匝线圈和瞬变电磁法的地下测深方法 |
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 |