瞬变电磁扩散场到虚拟波场的精细积分变换算法

现有的瞬变电磁解释与反演方法主要对电阻率参数进行解释,较难给出准确的地质构造信息.为了能够较多获得地质构造信息,可将瞬变电磁扩散场信号转换为虚拟波场信号,提高瞬变电磁法的分辨率.本文采用精细积分法,实现了瞬变电磁扩散场到虚拟波场的转换.将一个高度病态的线性方程组转换成求积分的过程,大大降低了解决病态问题的难度.积分步长以2的指数增加,在保证精度的前提下,具有极高的计算效率;设置合理的终止迭代条件,进一步提高本文方法的适用性.对典型地电模型的计算,表明本文方法具有较高的精度;减小两个波峰之间的距离,证明了本文方法具有较好的分辨能力;对含有噪声信号的扩散场进行波场反变换,证明了本文方法具有较好的抗噪性.最后计算三维模型和实测数据的虚拟波场,可以很好地反映出三维地质体的界面信息,证明本文方法稳定、可靠,可以对地下地质目标体的界面进行有效识别.本文方法相比预条件正则化共轭梯度法(PRCG),在计算效率方面提升了 4倍左右;本文方法的相对误差在5％以内,而PRCG的最大相对误差可达40％～50％,精度最高可提升10倍;在分辨率与稳定性方面,也表现出较好的效果.     

瞬变电磁虚拟波场的三维曲面延拓成像

根据瞬变电磁场扩散方程与波动方程间存在的数学对应形式,将已知瞬变电磁场数据转换为虚拟波场数据,在此基础上,提出了用克希霍夫积分法将所转换的波场从地面向地下反向延拓的方法.在延拓的数值计算中,为提高计算效率,采用三维边界元技术,把边界积分分解为诸三角单元积分的积分之和, 建立了曲面延拓方程式;通过地面测点的波场值,求出地下某一点的波场值,实现了地表为曲面的向下延拓成像计算.通过对理论模型计算和实际资料处理,证明了该方法可以增强瞬变电磁法识别地下电性分界面的能力,使瞬变电磁法对地下目标体的三维精细探测成为可能.     

考虑关断时间的回线源激发TEM三维时域有限差分正演

从麦克斯韦旋度方程出发可以直接导出瞬变电磁场扩散方程,然而扩散方程不含电场对时间的一阶导数,不能构成显式的时域有限差分方程,借鉴du Fort-Frankel有限差分离散方法引入虚拟位移电流项构建显式时域有限差分方程.对Wang和Hohmann的经典时域算法进行了两点改进:第一,通过将矩形回线源电流密度加入麦克斯韦方程组的安培环路定理方程,实现回线源瞬变电磁激发源加入;第二,在计算中考虑关断时间.第一点改进使时域有限差分方程考虑了一次场的计算,并且源的计算不再依赖均匀半空间模型响应作为初始条件,使算法能够适应表层电阻率不均匀时的三维复杂模型.由于实际观测中不可能出现阶跃电流的关断形式,第二点改进可以方便设置发射电流下降沿.采用改进的三维时域有限差分正演算法对均匀半空间模型、四类三层模型、均匀半空间中含有低阻块体模型进行了计算并分别与解析解、线性数字滤波解、积分方程解和Wang的三维时域有限差分解进行了对比验证.以H模型为例,采用建立的三维时域有限差分正演算法计算了不同关断时间的斜阶跃脉冲回线源瞬变电磁中心点感应电动势衰减曲线.以实际地质资料为基础,构建包含两层采空区的三维复杂模型,以1 μs的极短关断时间进行了复杂模型定回线源瞬变电磁响应计算,并计算了该复杂模型的视电阻率曲线.     

隧道瞬变电磁克希霍夫偏移成像与地质灾害探测

隧道绿色施工不仅需要预报工作面前方突水突泥等低阻异常,更需要查明灾害体的地质构造,为灾害预警与防治提供保障。传统的瞬变电磁解释方法仅能够提供单一电阻率信息不能满足隧道地质灾害治理要求。因此,将虚拟波场成像技术引入到瞬变电磁隧道探测技术中,期望实现灾害体电性与结构综合解释。首先,根据波场变换原理将瞬变电磁数据转换到虚拟波场;然后,利用电阻率成像算法计算工作面前方电阻率建立虚拟波场速度模型,利用波动方程克希霍夫积分解进行波场延拓实现瞬变电磁虚拟波场偏移成像;最后,综合电阻率与偏移成像结果对工作面前方地质体进行解释,判断地质灾害体的电性与构造特征。分别采用充水溶洞、充水断层两类常见灾害模型对算法进行验证。充水断层模型电阻率成像结果可以识别工作面前方低阻异常,但是异常分布范围略有增加,且不易判断倾斜角度;而偏移成像结果有效地识别异常的边界位置且对应准确,更易容判断异常倾斜方向。充水溶洞模型电阻率成像结果可识别溶洞的低阻异常与位置,但是异常形态与实际模型略有偏差;偏移成像结果圈定了模型的前后边界且与模型吻合较好。实测数据偏移成像结果有效圈定了低阻异常位置与裂隙形态,预测结果与已知出水孔情况相符。理论模型与实测数据结果表明,虚拟波场克希霍夫偏移成像结果既包含地质灾害体的电阻率分布,又能识别灾害体电性分界面,为地质灾害预警与防治提供了丰富的地质信息。      

西藏沙让特大型钼矿的发现与瞬变电磁法探查

近10年来,西藏冈底斯斑岩成矿带找矿取得重要研究成果,成为研究热点地区之一.西藏沙让地表为粘土化低温矿化蚀变、总体品位不够.穿插有小而富的脉状辉钼矿,规模有限,深部成矿状况不明.利用大回线源瞬变电磁法新技术进行工作,成功改装瞬变电磁接收装置;定义大回线源瞬变电磁晚期视电率计算公式;采用修正后的视电阻率计算公式进行数据处...     

瞬变电磁法三维拟地震成像信息提取技术

在瞬变电磁法中,为提高对地下电性界面的分辨率,可移植地震勘探中的诸多技术对界面直接成像,以克服反演中的不适定性.首先为使得波形在高低不同的电性基底所产生的虚拟反射波振幅有可比性,在将原有瞬变电磁场逆变换成虚拟波场数据的基础之上,需要求取均衡系数;接着通过对虚拟波场的离散数据求取反褶积,可基本抵消波场变换的展宽作用和大地滤波影响,以达到对电性界面精确定位的效果;然后,对波形做低通滤波处理后运用克希霍夫积分公式,将虚拟地震波场由地面向地下反向逆时偏移成像.把波场分析原理用于对瞬变电磁场的解释,就形成了瞬变电磁偏移方法.从波场角度看,它拓展和丰富了瞬变电磁场的内涵.使用本方法,可更加容易得到从前用常规解释方法处理无法提取的信息.在延拓成像中可更精确的确定界面深度,并为在瞬变电磁法中合理引入拟地震处理解释方法提供重要的理论依据和实验参照.     

瞬变电磁虚拟波场二阶Born近似成像算法

依据瞬变电磁扩散场与虚拟波动场之间的数学表达式,通过波场反变换得到虚拟波场,进而借用地震勘探上成熟的成像方法定位地质界面的位置以及形态.本文首先将时间域波动方程转换为频率域波动方程,将总波动场分为背景场和散射场、总波速分为背景速度和扰动速度,然后运用Green定理得到了散射场关于地质界面速度扰动量的表达式,并使用Bor...     

双轴各向异性介质中回线源瞬变电磁三维拟态有限体积正演算法

采用模拟离散的有限体积法实现了双轴各向异性地层回线源瞬变电磁三维正演.首先引入内积定义,采用自然边界条件,将瞬变电磁法的控制方程转化为弱形式表示.将计算区域划分为一系列的控制体积单元,采用交错网格对控制方程进行模拟有限体积空间离散,包括旋度算子离散和空间内积离散.基于斯托克斯定理的旋度积分定义公式实现旋度算子离散.中点平均实现电导率双轴各向异性的空间内积离散,从而得到离散化的控制方程.时间步迭代采用无条件稳定的欧拉后向差分格式.并通过均匀全空间中稳定电流回线源的磁场解析表达式得到回线源初始时刻的电磁场分布.为了同时保证计算精度和效率,本文采用分段等间隔的时间步迭代,利用直接法求解器PARDISO实现其快速求解.最后通过对比层状模型和各向异性半空间模型的正演计算结果,验证了本文算法的计算精度和计算效率;计算三维双轴各向异性模型的正演响应可知,水平方向电导率变化对电磁响应产生显著影响,而垂直方向的电导率变化对电磁响应几乎没有影响.产生这一现象的主要原因是回线源产生的感应电流主要是水平方向的,因此响应主要受到水平方向电导率的影响,垂直方向的电导率影响很小.     

考虑关断时间的地面瞬变电磁三维带地形反演

起伏地形和关断时间对地面瞬变电磁响应影响严重,这给传统基于水平地表模型和理论阶跃波形的瞬变电磁数据解释技术带来很大困难.为此,本文开展考虑起伏地形和关断时间的地面瞬变电磁三维反演算法研究.正演采用基于非结构网格和后退欧拉隐式时间离散格式的时间域有限元算法,快速模拟起伏地表模型瞬变电磁响应.反演采用L-BFGS算法,减少每次反演迭代的计算量.针对发射线圈随地表起伏变化的特点,利用基于偶极子离散的场源处理技术模拟发射源的实际形状,采用瞬时电流脉冲技术实现考虑关断时间的地面瞬变电磁三维正演模拟.我们首先将本文开发的三维反演算法应用于理论模型的反演计算中,检验本文算法的可靠性,并分析地形和关断时间对反演结果的影响特征.在此基础上,进一步将本文算法应用于实测数据反演,验证本文算法的实用性.     

瞬变电磁虚拟波场高分辨成像技术研究

本文提出了从瞬变电磁场到波场的优化算法.目的是通过对现有渡场变换过程进行插值的办法提高波场变换的分辨率,从而为后期数据处理的质量打好坚实的基础.在进行波场反变换的过程中,将正则化算法用于逆变换过程的计算中,通过采用偏差原理和Newton迭代格式选出最优的正则化参数,使得反变换所得到的波场稳定、可靠.解决了积分系数过多而...