基于任意六面体谱元法频率/时间域航空电磁三维正演模拟研究

来源 :吉林大学 | 被引量 : 0次 | 上传用户:yykk110
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
近年来,航空电磁勘探技术在近地表资源勘查、水文、工程与环境勘查等领域发挥越来越重要的作用。由于该技术具有采集数据量大、对地下浅层介质分辨率高等特点,其正反演需要高效率和高精度算法作为保障,这为航空电磁实际数据处理和正反演带来严峻挑战。目前,航空电磁一维正反演算法较为成熟,却不足以满足复杂地电模型的模拟与反演解释,已经开发的三维正演算法受计算精度与效率的双重限制,制约着实测数据三维反演和数据解释顺利进行。为同时改善航空电磁三维正演算法的效率与精度,本论文提出以Gauss-Lobatto-Legendre多项式为基函数的谱元法开展频率域和时间域航空电磁复杂三维地电模型正演模拟研究,具体分析不同物理网格与不同阶数基函数对航空电磁正演模拟精度与效率的改善程度,探究谱元法对各向异性地电模型、复杂形态异常体、起伏地形等条件下航空电磁响应模拟能力。为实现基于谱元法航空电磁正演,本论文首先采用规则六面体对地电模型进行剖分,分别从矢量Helmholtz方程和电磁场矢量波动方程边值问题出发,将全张量电阻率引入到正演模拟中,推导各向异性条件下频率域和时间域航空电磁正演边值问题,并对各向异性边值问题进行简化,给出各向同性条件下正演边值问题。为避免航空电磁场在发射源附近解的强烈奇异性,本论文在正演过程中利用场分离算法。针对频率域正演问题,以空气为背景的全空间解析解为一次场。针对时间域问题,以均匀大地半空间半解析解为背景场。为求解频率域问题,本论文以Gauss-Lobatto-Legendre多项式为基函数对计算空间进行离散,根据Galerkin加权余量法推导频率域问题的谱元控制方程,通过离散空间物理域与基函数空间参考域的映射关系完成单元矩阵分析过程,结合基函数在配置点处δ函数属性,采用降阶数值积分,以增加总体矩阵稀疏性。考虑到航空电磁多源问题,本论文采用直接求解器并行求解多源航空电磁正演问题,以提高多源问题求解效率。为求解时间域问题,本论文提出无条件稳定的后推欧拉差分格式与谱元法相结合的时间域谱元法,实现稳定而不受时间步长限制的时间域正演问题的直接求解,根据Galerkin加权余量法和后推欧拉差分格式推导了求解任意时间道航空电磁响应的时域谱元法控制方程;利用降阶数值积分形成的对角质量矩阵,简化总体矩阵与右端项计算。引入直接求解器并行求解多源多时间道时间域问题,在时间步长不变时,通过一次分解多次回代方式求解,避免了多时间道多次耗时分解,实现多源多时间道问题的快速求解。为检验频率域正演算法,本论文利用一维半解析解证实了基于规则六面体网格谱元法各向同性与各向异性航空电磁正演的精确性。并与有限元、有限差分算法对比,验证本论文算法高效性,和其在单一物理网格中对“非线性”变化的航空电磁响应刻画能力。本论文以频率域正演问题为主探究谱元法模拟电磁响应的综合能力,讨论不同物理网格剖分形式与不同基函数阶数对精度与效率的影响,根据谱元法自身特点,探究基于谱元法航空电磁正演计算的“二元”优化策略(即物理网格剖分与谱元基函数阶数的最佳结合方式)。针对典型各向异性模型,分析航空电磁各向异性响应特征,并给出航空电磁响应的各向异性效应的识别方法。为验证时间域正演问题算法的精度,本论文讨论了在不同物理网格剖分和不同基函数阶数条件下时间域正演计算的有效性,分析谱元基函数阶数和物理网格剖分对计算结果精度和效率的影响,讨论典型三维地电模型多时间道多源时间域航空电磁响应特征。为模拟复杂地电条件航空电磁响应,探究谱元法模拟的灵活性,本论文采用形变六面体单元离散地电模型,建立基于形变六面体网格谱元法频率域和时间域航空电磁正演方法。通过引入形变六面体单元形函数表示Jacobian矩阵,建立地电模型对应的物理坐标系与基函数对应的参考坐标系映射关系,推导了包含Jacobian矩阵相关函数的单元矩阵形式,并结合数值积分完成了谱元矩阵的分析。为优化复杂地电模型矩阵分析过程,本论文提出形变六面体网格与规则六面体网格相结合的方式建模,以不规则六面体单元拟合任意的物性界面,在物性不变区域采用规则六面体单元离散,以加速矩阵分析过程。通过与一维半解析解、有限元法、基于规则六面体谱元法对比分析,验证了基于形变六面体网格谱元法在频率域和时间域航空电磁正演中的有效性。为测试基于形变六面体网格谱元法的灵活性,本论文采用形变六面体网格对多种复杂地电模型(形变异常体、倾斜异常体、和球形异常体)的频率域和时间域航空电磁响应特征进行分析,探究异常体形变和倾斜程度对航空电磁响应的影响特征。为进一步探究算法的灵活性,本论文以频率域航空电磁正演为例,测试了基于形变六面体网格谱元法模拟起伏地形地电模型航空电磁响应的能力,并讨论起伏地形地电模型的航空电磁响应特征。为探究简单物理网格离散实现复杂地电模型航空电磁正演求解的可行性,本论文提出通过基于耦合无单元Galerkin法思想的谱元法进行航空电磁正演模拟。参考无单元Galerkin法节点离散物理模型的方式,本论文以数值积分节点离散物理模型。采用规则六面体单元对物理域进行剖分,从Galerkin加权余量法出发推导了频率域航空电磁正演控制方程。通过与一维半解析解和常规谱元法对比,验证改进谱元法的有效性,总结了基函数阶数与数值积分阶数对模拟结果的影响规律,并对复杂模型和实际地质体航空电磁响应进行模拟,以探究算法的灵活性。本论文系统地建立了规则和形变六面体网格谱元法的频率域和时间域航空电磁三维正演理论,并开发了相应算法,实现了基于谱元法任意各向异性复杂地电模型航空电磁三维正演模拟,为频率域与时间域航空电磁大数据量高精度和快速反演奠定基础。
其他文献
促性腺激素的异质性及其在家畜繁殖生理中的作用(中)D.J.Cooke等著李震傅国栋摘译4β-亚单位由于β-亚单位决定着促性腺激素的生物学特异性(Pierce等,1981),所以经常把β-亚单位叫做促性腺激素特异性亚单位
刘禹锡的<竹枝词>是唐诗宝库中一颗璀璨的明珠.作者以诗意的眼光看生活,努力去发现生活之美,以民歌的手法写生活.
<正>《环球军事》2012年2月上半月版刊登了一篇《老火炮打出新风采》的文章,写的是中国人民解放军沈阳军区某炮兵团的事迹。该文开篇颇具气势:"‘东北望,
期刊
为解决用户在线课程选择过程中忽略课程学习的进阶性,进而导致知识间出现断层的问题,借助“中国大学MOOC”平台,提出了一种基于规则库的用户选课健康状态感知及推荐方法。该
<正>360多年前,徐霞客曾四次到衢州进行旅行探险。在乡村野外、崇山峻岭中,跋山涉水,风餐露宿,以命相搏,孜孜不倦地追求科学真理,用实际行动进行乡村旅游的实践,矗立起"游圣"
会议
由于精神病人的心理活动处于严重紊乱和不稳定之中,常常发生兴奋、冲动、攻击和暴力行为.据调查,精神病患者的攻击行为发生率约为一般人群的10倍[1].
目的 对照比较成人近视不同验光方法间存在的差异。方法 选取2016年1月~2017年1月我院眼科收治的行散瞳电脑验光、视网膜检影的成人近视眼患者124例(248眼)为研究对象,研究球镜
利用三维荧光光谱研究了大沽河流域地下水中溶解有机物(DOM)的荧光组分类型及其空间分布与变化,并与地表河水进行了比较。基于平行因子分析模型在大沽河流域识别出2个类腐殖
<正>氮氧化物(NOx)是主要的大气污染物之一,氨选择性催化还原(NH3-SCR)被认为是消除NOx最有效的技术,催化剂是SCR技术的核心。在本工作中,我们开发了一种新型的介孔TiO2球负载的M