复杂介质中时间域航空电磁数据仿真技术研究

来源 :吉林大学 | 被引量 : 0次 | 上传用户:freecase
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
时间域航空电磁法采用机载移动平台,可以对高山、沙漠、沼泽等难以开展地面工作的地区进行快速高效的勘探,现已被广泛应用于资源勘查、环境工程等领域。然而,在数值模拟方面,现有的三维时间域航空电磁正演方法存在稳定性差、计算效率低、无法处理复杂地电结构等诸多问题,严重制约了三维时间域航空电磁数据处理和解释技术的发展。为了解决上述问题,本文采用基于非结构网格的时间域矢量有限元算法开展复杂介质条件下三维时间域航空电磁数据仿真研究,并对电磁响应特征及起伏地形和各向异性的影响规律进行分析。本文首先从时间域电场扩散方程出发,基于非结构四面体网格的灵活性对复杂几何模型进行精细剖分,准确拟合地下电性结构。然后,采用自动满足电场切向连续和无散条件的矢量插值基函数进行空间离散,并利用伽辽金方法建立有限元控制方程。通过引入后推欧拉格式进行时间域离散,获得无条件稳定的隐式方程,放宽对时间步长的限制,根据电流变化特征和电磁扩散规律灵活设定时间步长。同时,应用并行技术计算单元矩阵,并对系数矩阵进行压缩存储以提高计算效率,减少内存消耗。最后,利用多波前并行求解器MUMPS实现有限元方程的快速求解。当时间步长不变时,只需对系数矩阵进行一次分解,通过不断回代右端项即可实现航空电磁多时间道正演模拟。为解决多源问题,本文基于时间域航空电磁系统的有限影响范围(footprint),采取局部网格策略,将观测点分组并分别设计独立网格。通过与一维解析解进行对比,验证了本文算法的准确性。随后以三维地质体模型为例,分析局部网格策略的加速效果。为模拟复杂场源任意发射波形的时间域航空电磁响应,本文提出基于偶极子离散和瞬时电流脉冲的场源处理技术。首先将任意姿态、形状的发射线圈离散成若干段首尾相连的导线,并近似成电偶极子,实现场源的精确模拟,进而根据数值实验分析发射源姿态和形状变化对电磁场的影响规律。然后,本文将瞬时电流脉冲直接加载到源项中,通过改变各个时刻的脉冲强度模拟任意发射波形航空电磁响应。由于该方法无需计算电流二阶导数,数值稳定性得到有效提高。通过计算三维模型实际发射电流波形的航空电磁响应,并与传统褶积算法结果进行比较,检验本文算法模拟复杂波形航空电磁响应的优越性。鉴于网格剖分对三维数值模拟精度影响较大,本文将面向目标的自适应网格加密技术应用于时间域航空电磁三维正演模拟中。首先,基于法向电流密度连续性条件进行后验误差估计;然后通过在观测点处施加虚拟点源计算影响函数,调整各个单元后验误差的影响权重;最后根据加权误差进行局部网格加密,改善网格剖分质量,进而达到提高数值模拟精度的目的。通过对均匀半空间和起伏地表模型上航空电磁响应进行模拟,检验本文自适应算法的可靠性。基于上述研究,本文首先采用时间域有限元算法对球体、直立/水平/倾斜板状体、基底隆起和断层等典型地电模型航空电磁响应进行正演模拟,并分析响应曲线特征。结果显示,不同模型的电磁响应差异较大,根据响应曲线形态和峰值位置可以对目标体的产状和水平位置进行判断,为异常拾取和识别提供参考。其次,本文研究山峰、山谷等起伏地表模型的航空电磁响应特征,分析地形影响规律,验证本文算法处理复杂三维地质结构的有效性。数值实验结果表明,地形对航空电磁响应影响严重,且主要集中在早期,并随时间逐渐减弱,响应形态与地形几何形状呈镜像关系。起伏地表与地下异常体相互耦合使得观测信号变得十分复杂,给数据解释造成困难。最后,本文模拟任意各向异性大地复杂地电结构模型的时间域航空电磁响应,并分析各向异性影响特征。对各向异性异常体、各向异性围岩以及起伏地形各向异性大地三种典型模型的研究结果表明各向异性对航空电磁响应的强度和分布形态均产生严重影响,根据异常分布特征可以判断各向异性特征(电导率主轴方向和旋转角度)。本文研究的时间域航空电磁仿真技术可以为航空电磁系统研发和施工设计提供理论数据,同时为复杂介质条件下三维航空电磁数据反演提供正演模块。本文的研究成果对提高航空电磁数据解释水平具有一定的理论和实用价值。
其他文献
借助开放式创新理论、组织学习理论和动态理论,探讨并厘清参与者惯例变更、贡献一致治理变量,在合作创新过程中的作用路径和治理效果。对368份来自67个参与合作创新的企业数
本文对Chapman阳光掠射函数[Ch(z_p,χ)]进行了数值积分,求得了其在实际模式大气中随观测高度z_p及天顶角χ的变化.计算并讨论了低热层高温度梯度、分子与湍流扩散、重力场及
本文对中国东部地磁台站链地磁静日变化东向分量进行了分析,有结果表明,在冬、夏季清晨,南北半球之间存在着电离层发电机驱动的场向电流,其方向是从夏季半球到冬季半球;在黄
本文讨论宽频带静电噪声对等离子体片边界层处粒子的随机加速问题.计算表明:宽频带静电噪声会使离子轨道在相空间Poincare截面上一定范围内无规行走,该范围与扰动电场的振幅
Iss-b卫星探测结果表明,|电离层离子丰度呈现经度变化.本文选取较真实的电离层模型参数,对130°E和114°E两个子午面内传播的低纬哨声进行了射线跟踪.我们发现,由于[
精密检调焦技术是空间相机必不可少的重要组成部分,它直接影响相机的成像质量.本文介绍了一种自准法精密检调焦系统的方案,并对测量数据的处理方法进行了分析.实验结果表明,
对ISEE-3人造卫星在1980年5月—1981年8月中,观测到的48个X射线耀斑进行了分析,发现其中有1/3是在6个活动区中重复爆发的.研究这部分X射线耀斑的物理性质与所在活动区的黑子
本文用相空间重构方法,把电离层参数的时变序列重构一个m维相空间,研究相点之间的距离关系,提取关联维数和Lyapunov特征指数,从而揭示电离层参数的混沌特性,及参数的最大可预报时间尺度。
甭哥金矿床产于青藏高原东缘玉树-义敦弧南端。甭哥矿区发育超钾质正长岩有别于与俯冲相关的钙碱性中酸性岩。义敦弧岩浆带多发育斑岩型Cu矿床,甭哥金矿作为少数已知的独立Au
本文从铝合金铸件气孔类别分析入手,阐述了铝合金铸件气孔的成因以及其对铸件性能的影响,并针对铝合金铸件气孔形成的主要因素提出了相应的预防措施,并扼要介绍了铝合金铸件