基于区间B样条小波有限元法的频率域电磁法三维多分辨率正演研究

来源 :吉林大学 | 被引量 : 0次 | 上传用户:xyw6623
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
频率域电磁法是一种高效的地球物理勘探手段,能够服务于矿产与油气资源勘探、环境与工程勘查等领域,从而推动国家基础设施建设和经济发展。开展频率域电磁法的三维正演研究能够促进三维反演技术的发展,同时为野外实际工作提供理论支撑。为了改善数值模拟算法的性能,本文引入基于区间B样条小波的小波有限元法开展频率域航空电磁法(Airborne EM)与海洋电磁法(Marine controlled-source EM)的三维正演模拟研究。小波有限元是一种将小波理论中的尺度函数或小波函数作为有限元基函数的数值模拟方法。相比于常规的多项式插值,小波具有更加优异的数值性质,因此能够改善有限元法的性能。小波有限元能够采用不同尺度的基函数近似单元中未知的场分布,实现多分辨率的求解。同时,由于小波的稀疏性,小波有限元方程的系数矩阵通常更加稀疏,有利于加快求解速度。区间B样条小波定义在参考闭区间[0,1]内的节点序列上,并由节点间的分段B样条多项式组合而成。区间B样条小波具有显式的表达式。同时,区别于传统的定义于实数轴上的小波,区间B样条小波定义在闭区间内,且对边界处的尺度函数与小波函数进行了特殊的定义,确保了其在边界处的数值精度。本文由Maxwell方程出发,定义磁矢量势与电标量势,并基于库仑规范推导得到频率域耦合势方程。为了克服场源奇异性,减少网格数,本文采用场分离算法,以背景电场构造右端项并对二次耦合势方程进行求解。对于航空电磁正演问题,本文采用磁偶极子在空气全空间中激发的响应作为背景场。由于其响应存在解析解,因此可以直接求得。对于海洋电磁正演问题,本文采用电偶极源在层状介质中激发的响应作为背景场。由于其响应不存在解析解,本文通过一维正演程序计算其半解析解。进而,我们利用区间B样条小波的三维张量积构造有限元基函数对单元内未知的二次势进行近似。通过规则六面体网格对计算区域进行剖分,并通过伽辽金法离散控制方程,形成大型线性方程组。对于航空电磁问题,本文采用直接求解器MUMPS对多源问题进行快速回代求解。对于海洋电磁问题,本文采用稳定化双共轭梯度法以及代数多重网格预条件子进行求解以减少内存需求,实现大尺度问题的正演模拟。最终,通过移动最小二乘插值对二次势进行处理得到测点处的电磁场响应。对于航空电磁正演问题,为验证本文算法有效性,本文采用不同的网格尺寸与基函数尺度,基于均匀半空间模型对算法精度进行验证。同时,与常规的节点、矢量有限元法的对比表明本文算法可用较少的时间获得精度更高的正演结果。对三维低阻异常体的航空电磁正演模拟结果表明:本文算法与谱元法结果吻合情况良好。进而,本文对典型模型的航空电磁正演响应进行了模拟与分析,并对响应特征进行了总结。对于海洋电磁正演问题,本文首先模拟了海底存在高阻层的水平层状模型的正演响应,通过与一维半解析解的对比验证了本文算法精度。而后,本文分析了网格尺寸与基函数尺度对海洋电磁响应和正演效率的影响。通过与Jacobi预条件子、SOR预条件子进行对比验证了代数多重网格预条件子的有效性。最终,本文分析了海水深度和异常体埋深对海洋电磁响应结果的影响。数值算例表明:基于区间B样条小波的小波有限元法具有较高的数值模拟精度,能够以不同的尺度对正演问题的解进行近似,同时可以采取细化网格尺寸和提高基函数尺度两种策略灵活地提升正演计算精度。
其他文献
学位
学位
Sawan气田开发生产区位于Thar沙漠,在地质上形成了中印度河盆地的台地斜坡区。构造上它位于Jacobabad古高地的东南边缘。Sawan气田被两条主要的走滑断层划分为南、中、北三个部分。本文所研究的三口井都位于Sawan气田的中部。本文利用油气物理测井的测井相解释对Sawan气田储层质量进行了评价和分析,以全面了解Jacobabad高压南侧储层沉积学的地下特征。岩石物理参数由页岩体积、密度、中
学位
学位
学位
学位
学位
学位
地球物理正反演的核心任务之一是通过求解麦克斯韦方程组模拟准静态区域的电磁场。正演模拟作为电磁数据反演的基础,是电磁数据三维反演技术达到实用化的关键一环。然而,地下目标体形态复杂,多维数值模拟难以开展。此时,为获得精确的模拟结果,需要对地下半空间进行精细的离散化,导致大型线性系统的求解难题。由此三维反演技术无法做为一种常规的数据处理手段,用于实际地球探测问题。因此,实现对大尺度模型的高效正演是电磁三