第1章 绪论 1
1.1大地电磁法数值模拟国内外发展现状 1
1.1.1积分方程法 2
1.1.2有限差分法 3
1.1.3有限单元法 4
1.1.4矢量有限单元法 6
1.1.5自适应有限单元法 7
1.2大地电磁法反演国内外发展现状 8
第2章 大地电磁法三维矢量有限元模拟 13
2.1 MT的三维边值问题 13
2.1.1 Maxwell方程组 13
2.1.2边界条件 14
2.1.3准静态极限下的电场矢量波动方程 15
2.1.4 MT三维边值问题 15
2.2三维MT的伽辽金有限元公式 17
2.2.1伽辽金（Galerkin）法 17
2.2.2三维MT边值问题的Galerkin有限元方程 19
2.2.3截断边界条件的处理 20
2.3三维MT矢量有限元分析 21
2.3.1网格剖分 21
2.3.2三维矢量形函数 22
2.3.3三维矢量单元分析 24
2.4数值模拟算例 27
2.4.1均匀大地模型 27
2.4.2水平三层模型 31
2.4.3单个异常体模拟 34
2.5本章小结 37
第3章 非结构化三维MT矢量有限元模拟 38
3.1 D网格及优化 38
3.1.1映射法 39
3.1.2 AFT方法 40
3.1.3栅格法 41
3.1.4 Delaunay方法 42
3.1.5非结构化网格局部加密技术 4
3.2四面体的矢量有限元单元分析 48
3.3坐标变换及矢量形状函数 49
3.3.1坐标变换 49
3.3.2矢量形状函数 50
3.4单元矩阵计算及合成 52
3.5大型稀疏复线性方程组求解 54
3.5.1大型稀疏复线性方程组求解方法 54
3.5.2 Krylov子空间方法 56
3.5.3 GMRES（m）求解器 56
3.6数值计算与分析 58
3.6.1均匀大地模型 59
3.6.2水平四层模型 61
3.6.3 COMMEMI3D-1模型 62
3.6.4 COMMEMI3D-2模型 63
3.6.5单个高阻体模型 65
3.6.6组合异常体模拟 66
3.7本章小结 75
第4章 三维MT自适应矢量有限元模拟 76
4.1自适应有限元模拟的数学基础 76
4.1.1 Sobolev函数向量空间 76
4.1.2 Hemlholtz空间分解 77
4.2 MT模型的残差后验误差估计 78
4.3自适应策略 82
4.4数值模拟算例 84
4.4.1均匀大地模型 84
4.4.2 COMMEMI3D-1MT模型 88
4.4.3两个高阻异常体模型 88
4.5本章小结 91
第5章 被动源电磁测深双模联合反演 94
5.1基于最小构造的双模联合反演 94
5.1.1双模联合反演的目标函数 94
5.1.2目标函数的约束机制 95
5.1.3最小构造模型约束 96
5.1.4迭代求解策略 99
5.2提高正反演精度与速度的策略 101
5.2.1基于二次场的有限元正演模拟 101
5.2.2提高反演速度的双重网格方案 105
5.3模型反演算例与分析 106
5.4本章小结 160
第6章 反演灵敏度分析及灵敏度矩阵快速计算 162
6.1灵敏度矩阵求取方法 162
6.2音频大地电磁法灵敏度分析 165
6.2.1一维灵敏度分析 166
6.2.2二维灵敏度分析 169
6.3灵敏度矩阵快速计算 173
6.3.1灵敏度矩阵快速计算方案 173
6.3.2模型试验及分析 183
6.4本章小结 191
第7章 基于背景分离的被动源电磁测深双模联合反演 192
7.1多分辨小波分解理论 192
7.1.1小波变换的定义 192
7.1.2多分辨（多尺度）分析 193
7.2被动源电磁测深背景分离方法 194
7.3基于背景分离的被动源电磁测深双模联合反演方法 202
7.3.1背景分离反演的思想 202
7.3.2模型计算及结果分析 202
7.4实测数据反演与解释 209
7.4.1测区地质情况 209
7.4.2音频大地电磁数据一维分析 209
7.4.3音频大地测深数据二维反演与解释 210
7.5本章小结 219
参考文献 220