基于MPI的频率域航空电磁法有限元二维正演并行计算研究

【类型】学位论文

【作者】李小康 

【关键词】 频率域航空电磁法 有限元二维正演 并行计算 MPI

【摘要】航空电磁法(Airborne Electromagnetic Method)是航空物探中常用的观测方法之一,广泛的应用于地质普查、资源勘探和环境监测等领域。本文针对吊舱式直升机频率域航空电磁系统IMPULSE的各方面特点进行深入分析,提出了二维正反演研究的重要性。频率域航空电磁法二维正演问题实际上是,求解三维场源在二维地电断面中引起的电磁响应,也被称为2.5维问题。利用电磁场沿走向以场源中心呈对称分布的特点,可通过沿走向的傅立叶变换将空间域2.5维问题简化为波数域二维问题进行求解。本文首先按照场分离的思路,对一次场和二次场进行分离,通过将背景电导率分布设为全空气介质,使用解析解计算一次场,而二次场则使用有限元法求解。二次场求解过程中,使用伽辽金加权余量法和抛物线等参单元法进行处理。在求解大型线性方程组时,没有对整体刚度矩阵进行集成,而是选择波前法求解,从而避免的对整体刚度矩阵的存储,提高了计算效率。在完成公式推导后,根据得到的算法设计流程图,使用Fortran语言编写了串行代码。精度检验的结果显示,该频率域航空电磁二维正演算法具有相当的可靠性。随后,对8组不同的地下异常体模型和观测装置组合进行正演计算,并对得到的电磁响应结果进行分析。鉴于频率域航空电磁法多维正反演问题的复杂性,并为将来的反演计算打下基础,本文尝试对提出的正演算法进行并行化。Message Passing Interface(MPI)是用于并行编程的一种主要标准和规范,通过对其深入研究和对提出的频率域航空电磁法有限元正演模拟算法的分析,选择了域分解算法设计和SPMD并行编程思路。通过在小型网络机群的测试,正演计算的效率得到了大幅提升。

【学位名称】博士

【学位授予单位】中国地质大学（北京）

【学位授予年度】2016

【导师姓名】谭捍东

【读秀链接】读秀链接