基于谱元法的大地电磁二维数值模拟

基于谱元法的大地电磁二维数值模拟

论文摘要

大地电磁法(又称为大地电磁测深法,简称MT)是以天然交变电磁场为场源的一种探测地下电性结构的地球物理勘探方法,在能源勘探、矿产普查、地质调查、地震预报和工程地质等领域有着广泛的应用。近年来随着数值算法的进步、数据采集系统的提高和高性能计算机的迅速发展,大大促进了电磁法数值模拟技术的发展,出现了多种高效的数值算法,有效提高了数值模拟的速度和精度。其中谱元法作为一种新的数值方法在近十几年来得到快速发展,为了提高大地电磁数值模拟的精度和效率,本文将谱元法引入到大地电磁测深领域当中。谱元法是有限元和谱方法结合的一种数值模拟方法,具备有限元法处理复杂结构模型的灵活性及谱方法的高精度和指数收敛性,使得谱元法在流体力学、动力学、声学和地震波场等领域得到广泛的应用。近几年来,谱元法已实现了在海洋可控源电磁的数值模拟,但还未有学者将谱元法应用到大地电磁领域的数值模拟,因此本文将谱元法应用于大地电磁二维正演模拟当中。本文介绍了谱元法的理论基础,从频率域的麦克斯韦方程组出发,详细推导了MT二维正演的所满足的偏微分方程和二维介质大地电磁场所满足的边值问题,基于Galerkin加权余量法,在求解域内采用矩形单元进行剖分,单元内对未知量(电磁场值)在插值节点处插值基函数选用Gauss-Lobatto-Legendre(GLL)正交多项式,采用GLL数值积分最后形成大型线性方程组,用Intel的Pardiso求解器直接求解线性方程组得到电磁场分布,实现了大地电磁二维谱元法模拟。为了检验本文算法和程序的可靠性,在FORTRAN平台上开发了二维大地电磁谱元法正演程序,并给出层状模型和国际标准COMMEMI 2D-1模型进行对比验证,数值模拟结果表明了本文编写谱元算法程序的可靠性和稳定性。此外还计算了两个低阻异常模型和起伏地形模型,并讨论了地形模型对MT二维正演模拟的影响。

论文目录

  • 摘要
  • Abstract
  • 第1章 绪论
  •   1.1 选题背景及研究意义
  •   1.2 国内外研究现状
  •     1.2.1 大地电磁法正演问题研究进展
  •     1.2.2 谱元法在电磁领域的研究进展
  •   1.3 本文研究的方法及结构纲要
  • 第2章 谱元法理论基础
  •   2.1 正交多项式
  •   2.2 Chebyshev正交多项式
  •   2.3 Legendre正交多项式
  •   2.4 Fourier展开
  • 第3章 基于谱元法的大地电磁法的正演理论
  •   3.1 大地电磁正演理论
  •     3.1.1 控制方程及边界条件
  •     3.1.2 Galerkin法求解MT边值问题
  •   3.2 构造谱元空间及基函数
  •     3.2.1 网格剖分及映射关系
  •     3.2.2 构造二维空间基函数
  •   3.3 微分方程的离散形式求解
  •   3.4 边界条件的施加
  •   3.5 视电阻率及相位的计算
  •   3.6 线性方程组的求解技术
  •     3.6.1 迭代法
  •     3.6.2 直接法
  •   3.7 压缩存储技术
  •   3.8 谱元法程序在MT的实现过程
  • 第4章 数值算例
  •   4.1 一维层状介质模型
  •   4.2 国际标准模型COMMEMI2D-1
  •   4.3 两个低阻异常体模型
  •   4.4 起伏地形模型
  •     4.4.1 地垒模型
  •     4.4.2 地堑模型
  • 第5章 结论及建议
  • 参考文献
  • 附录
  • 致谢
  • 文章来源

    类型: 硕士论文

    作者: 方小姣

    导师: 李长伟

    关键词: 谱元法,大地电磁,正交多项式,二维正演

    来源: 桂林理工大学

    年度: 2019

    分类: 基础科学,工程科技Ⅰ辑

    专业: 地质学,地球物理学,矿业工程

    单位: 桂林理工大学

    分类号: P631.325

    DOI: 10.27050/d.cnki.gglgc.2019.000038

    总页数: 63

    文件大小: 1511K

    下载量: 21

    相关论文文献

    • [1].谱元法在求解刚架结构动力学问题中的应用[J]. 动力学与控制学报 2012(01)
    • [2].谱元法叠前逆时偏移研究[J]. 地球物理学进展 2008(03)
    • [3].考虑剪切变形的复杂刚架结构动力学特性分析[J]. 中国机械工程 2014(14)
    • [4].基于谱元法的板结构中导波传播机理与损伤识别[J]. 机械工程学报 2012(21)
    • [5].基于谱元法的空间刚架动力学特性分析[J]. 计算力学学报 2016(05)
    • [6].谱元法与透射边界的配合使用及其稳定性研究[J]. 工程力学 2015(11)
    • [7].卫星周期桁架结构振动特性应用研究[J]. 噪声与振动控制 2016(02)
    • [8].谱消去粘性谱元方法求解对流扩散方程[J]. 华电技术 2017(10)
    • [9].透射边界与谱元法的结合及对波动模拟精度的改进[J]. 振动与冲击 2017(02)
    • [10].谱元法应用于涡声传播问题的研究[J]. 西安交通大学学报 2016(11)
    • [11].谱元法模拟地裂缝场地地表地震波传播[J]. 地震工程与工程振动 2014(S1)
    • [12].复杂构造网格化及高精度地震波场谱元法数值模拟[J]. 石油物探 2014(04)
    • [13].内源输入时盆地特征参数对盆地地震效应的影响[J]. 苏州科技大学学报(工程技术版) 2018(03)
    • [14].最优控制问题弱间断解的一个自适应算法[J]. 控制与决策 2013(02)
    • [15].Galbrun方程谱元法数值计算[J]. 九江学院学报(自然科学版) 2010(02)
    • [16].核密度估计法在板件概率损伤识别中的应用[J]. 振动.测试与诊断 2015(06)
    • [17].汶川大地震地震波传播的谱元法数值模拟研究[J]. 中国科学(D辑:地球科学) 2009(04)
    • [18].基于谱元法的二维大地电磁数值模拟[J]. 矿产与地质 2019(05)
    • [19].谱元方法求解对流扩散方程及其稳定性分析[J]. 西安交通大学学报 2015(01)
    • [20].基于多辛结构谱元法的保结构地震波场模拟[J]. 地球物理学报 2012(10)
    • [21].人工源宽角地震实验的谱元法模拟和逆时偏移成像方法研究[J]. 中国地震 2020(02)
    • [22].透射边界条件在波动谱元模拟中的实现:二维波动[J]. 力学学报 2017(04)
    • [23].极坐标系下Fourier-Legendre谱元方法与有限差分法数值扩散的比较[J]. 计算力学学报 2013(03)
    • [24].基于谱元法的梁结构中Lamb波传播特性研究[J]. 噪声与振动控制 2009(06)
    • [25].多网格谱元法及其在高性能计算中的应用[J]. 应用声学 2018(01)
    • [26].聚集单元谱元法在承受冲击载荷结构动态分析中的应用[J]. 噪声与振动控制 2016(06)
    • [27].弹性水平成层场地稳态波动的谱单元模拟[J]. 南京工业大学学报(自然科学版) 2013(01)
    • [28].基于谱元法的各向异性黏弹性介质地震波场模拟和分析[J]. 世界地质 2017(01)
    • [29].基于谱元法的玉田震害异常研究[J]. 国际地震动态 2015(04)
    • [30].界面剥离钢-混凝土组合结构应力波传播谱元法模拟研究[J]. 工程力学 2017(02)

    标签:;  ;  ;  ;  

    基于谱元法的大地电磁二维数值模拟
    下载Doc文档

    猜你喜欢