时间域航空电磁法2.5维有限元模拟

采用三角网格剖分的有限元法,研究了2.5维航空瞬变电磁法正演模拟问题。利用时频变换数值方法将时间域电磁场转换到拉氏域,再利用傅里叶变换将三维问题降维变为2.5维问题,然后由有限元法求解得到拉氏域二维电磁场,逆拉氏变换后得到时间域航空瞬变响应。为了回避正演模拟中总感应磁场在场源处的奇异性问题,采用异常场算法,场源响应通过在微分方程中施加背景电磁场实现。由于瞬变电磁信号具有较大的动态范围,而且需要经过两次正、逆拉氏变换和傅里叶变换,每个环节的计算精度和速度要严格控制在较高的水平上,否则积累误差会非常大。模型计算表明均匀大地和层状大地模型解析解与数值解吻合很好。这证明该算法是正确可行的,可作为研究二维复杂地质体的方法手段。     

定回线源瞬变电磁三维异常特征反演

分析了实现定回线源瞬变电磁三维异常特征反演的可行性,简要地介绍了改进的广义逆反演方法基本原理,给出了定回线源瞬变电磁三维异常特征反演的流程。对正演模拟数据进行反演的实例表明,定回线源瞬变电磁三维异常特征反演已得到初步实现。     

基于非结构化网格的瞬变电磁2．5维有限元正演模拟

利用Delaunay三角化这种网格非结构化方法。通过编程实现了二维模型的非结构化三角形网格剖分,并编写了中心回线法瞬变电磁2．5维有限元正演程序。与前人计算结果对比,在取得相同计算精度的情况下,与结构化网格相比,非结构化网格所需网格和节点数量大大减少,计算效率更高。通过将非结构化网格法引入到瞬变电磁2．5维正演模拟中,实现了对复杂二维地电模型的有限元数值模拟,提高了现有有限元算法的应用范围。     

ANSYS在瞬变电磁三维正演中的应用

瞬变电磁正演是研究瞬变电磁响应规律的有效途径,瞬变电磁的高精度、快速三维正演技术,是当前瞬变电磁正演研究的热点.ANSYS软件中的Multiphysics模块也适用于瞬变电磁法的正演模拟问题.本文通过一些模型算例介绍了ANSYS软件在瞬变电磁场三维正演中的应用方法,包括建模、加载、求解设置、后处理的一般过程.数值计算结...     

线源频率电磁测深二维正演(一)

随着频率电磁测深(包括可控源声频大地电磁测深CSAMT)的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础      

线源频率电磁测深二维正演(一)

随着频率电磁测深(包括可控源声频大地电磁测深 CSAMT)的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础。     

用线性数字滤波法计算大回线源在地下形成的瞬变电磁场

本文讨论了地面大回线源在地下形成瞬变电磁场的正演计算方法。其正演计算是将频率域中相同模型、相同装置条件下的谐变场结果转换到时间域来。计算方法上采用了安德森等人提出的线性数字滤波法。计算过程中采用了三次样条函数插值法，提高了计算速度。通过均匀半无限空间解析解与数值解的计算结果比较表明，两者吻合很好，说明计算方法是可信的。     

线源频率电磁测深二维正演(二)

随着频率电磁测深（包括可控源声频大地电磁测深ＣＳＡＭＴ）的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础。     

任意角度瞬变电磁测深正演问题计算方法

针对目前比较普遍的瞬变测深晚期场计算精度较低的情况,本文依据频率域电磁场实虚分量的奇偶性及付氏变换的导数性质,利用余弦变换实现了由频率域到时间域电磁场的转换。计算表明,本算法计算晚期瞬变电磁及视电阻率响应具有较高的精度。文中重点计算了单一偶极供电,多点测量的任意角度瞬变测深的正演问题,并对其规律性作了总结。     

MT有限元模拟中截断边界的影响

利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。     

有限单元法求解水中二维体问题正演研究

潜艇在现代海军中占有举足轻重的地位,所以反潜技术也越来越重要,相对于传统的声波、磁法等反潜技术,电磁法探测是近几年才开始起步的。这里研究了电磁法探测水下二维体的问题,基于天然电磁场的有限单元正演,同时考虑电阻率与磁导率,采用三角剖分,通过模型试算验证了天然电磁场的有限元算法在水下二维体探测中,具有较高的精度和适用性,模拟结果 TM模式优于TE模式,可以为电磁法探测水下潜艇打下坚实基础。     

电磁波测井的实验技术及结果分析

为了对多频电磁波测井有更深入的了解,结合实验研究工作,针对油田的地质情况,研制了模拟实际测井环境的油水层模型,归纳出一套较为系统的实验方法和测量技术,并对测量结果进行了分析和讨论,取得了较好的实验效果.     

多辐射场源地空瞬变电磁一维反演方法研究

地空瞬变电磁法融合了地面与航空瞬变电磁的方法原理与观测模式,将发射源置于地表,并利用航空器在空中接收响应信号。因此,地空瞬变电磁法同时具有地面装置大测深、高信噪的优势以及航空装置中采集高效的优势。先有研究偏重于单一辐射源,忽略了多辐射场源在勘探深度与信噪比上的优势。因此需要研究适用于多辐射场源地空装置的数据解释方法。在解决多辐射场源地空瞬变电磁正演问题之后,以Occam反演理论为基础,对多辐射场源地空瞬变电磁一维反演进行研究。首先要构建目标函数,本文通过引入拉格朗日乘子,将数据拟合差及模型粗糙度组合起来进而构建目标函数;其次解决灵敏度矩阵的构建问题。模型计算结果表明本文多辐射场源地空瞬变电磁法一维反演算法的有效性,可为三维反演与精细解释提供良好的基础。通过对典型模型的反演计算,从其理论模型及反演结果的对比可以看出,Occam反演方法可以应用到多辐射场源地空瞬变电磁系统中,为多辐射场源地空瞬变电磁解释增加了新途径。     

基于Julia加速的射频大地电磁法正演模拟与位移电流作用研究

射频大地电磁法(Radio-magnetotelluric, RMT)是浅地表电磁勘探的重要手段之一。由于勘探频段为10～300 kHz,其电磁场传播受地下介质介电常数的影响较大。传统准静态条件下的电磁响应严重制约RMT正演模拟精度,并进一步影响反演成像分辨率。针对这一问题,提出了一种基于Julia并行加速的全电流RMT电磁响应数值模拟方法,利用Julia的分布式计算将各频点的计算发送到不同进程进行求解,从而达到加速计算的目的,同时在计算中考虑位移电流的影响,提升正演的模拟精度。通过计算几种典型高阻/高介电模型的RMT响应,分析并总结了位移电流对射频段电磁场视电阻率及相位响应的影响规律。数值模拟结果表明：当浅部存在高阻覆盖层时,基于准静态假设条件计算的RMT视电阻率和相位响应偏高,且频率越高、覆盖层电阻率越大,响应偏差越大;对于煤炭采空区模型,RMT法能有效反映异常体位置,但忽略位移电流会在采空区及其附近引起较大计算误差;起伏地形的算例表明地形会覆盖地下异常体的RMT数值响应,尤其是地形拐角处;2种不同规模的并行算例对比证明了并行算法的高效性,且随着求解问题规模增大,并行算法效率也随之...     

电导率分块线性变化的二维高频大地电磁法有限元正演模拟

为了适应电性参数在水平方向和垂直方向是连续变化的实际需要,这里采用矩形剖分,线性插值的有限元法研究了电导率分块线性变化的高频大地电磁高精度,快速正演模拟。首先给出了二维高频大地电磁测深的变分问题以及电导率分块线性变化时的有限元数值解法,编制了相应的有限元程序。利用该程序对层状模型进行了计算,并与解析解的结果做了对比分析,证明了程序的正确性和有效性,然后对典型电导率随深度线性变化的模型和典型地堑模型进行了正演模拟,结果表明能够有效地解决电导率在水平和垂直方向分块线性变化的高频大地电磁正演问题。     

地一井瞬变电磁三维交错网格有限差分正演及响应特性

为了研究地-井瞬变电磁响应特征,获得不同井位低阻薄板异常体的响应规律:首先采用交错网格有限差分技术离散二次场满足的频率域赫姆霍兹方程;然后结合虚框叠加等效和虚拟界面法将发射源和接收位置扩展到任意层位,解决了全空间背景格林函数的计算问题;之后利用MUMPS求解器求解频率域二次场,再经过余弦变换,获得井中任意位置的瞬变电磁响应;最后采用三维模型对本文算法的正确性进行验证,设计了均匀半空间和嵌入低阻薄板的三维模型,获得垂直磁场的三维分布,分析均匀导电半空间中低阻薄板对地-井瞬变电磁三维响应的影响特点.结果表明:本文三维地-井瞬变电磁计算方案的计算精度与前人基本相同;水平导电薄板的存在主要影响板体附近及穿过薄板的井中瞬变响应,最大的特点是在异常体位置附近的中期时间道响应出现变号现象.本文的研究为定性解释地-井瞬变电磁法异常提供一个技术手段,也为地-井瞬变电磁三维反演奠定基础.     

基于传输线法模拟探地雷达在多层介质模型中的传播

数值模拟是研究探地雷达电磁波在多层介质中传播规律的有效手段,传输线算法是一种基于Huygens原理的时域分析电磁场方法。这里阐述了TLM方法的基本原理及基本关系式,用TLM方法对多层介质进行了正演模拟,将模拟结果与FDTD方法合成结果进行了对比,结果表明了TLM用于GPR正演计算的正确性及有效性。应用该算法对已有模型进行了定性分析,验证了随着介电常数和电导率增大,雷达电磁波衰减加快的结论。     

CSAMT一维全频视电阻率反演

本文采用全频率视电阻率反演,使用2~(-3)~2~(13)Hz频率,并且不做近区和过渡区校正,直接用卡尼亚视电阻率作为反演参数。正演采用虚界面法计算有限长导线在水平层状介质的电磁场分量,反演采用有限内存拟牛顿法。数值模拟采用三层模型和实测数据进行反演,得到了符合地下电性变化的反演模型结果,证实了CSAMT一维全频视电阻率反演具有可实施性。     

基于移动阵列电极的盾构超前探测正演响应分析

针对现有超前探测技术探测空间狭小、电磁干扰严重且难以应用于隧道掘进机（Tunnel Boring Machine,简称TBM）施工隧道等难题,采用一种基于移动阵列电极的盾构超前探测方法。运用Comsol Multiphysics有限元仿真软件建立地电模型,对不同地质情况进行正演响应分析。结果表明:存在保护电极时,电压分布范围明显扩大,表明探测范围增大,为提高探测效果,有必要添加保护电极;异常体的尺寸、形状、位置改变会影响测量电压分布呈现规律性变化,表明这种方法对异常反应的响应良好;通过移动阵列电极盾构超前探测物理模型实验结果对比证明了基于移动阵列电极的盾构超前探测方法的可行性。