摘要：正演数值模拟是反演成像的基础.为了实现磁法勘探精细化反演成像与定量解释,本文利用把一个大问题分解为多个小问题的思路,提出一种基于泊松方程的高效、高精度空间波数混合域二度体磁异常数值模拟方法.该方法利用傅里叶变换把二度体磁位偏微分方程转换为一维常微分方程,并采用基于二次插值的一维有限单元法求解该方程,进而通过反傅里叶变换得到空间域磁异常.在模型算例中,分别设计截面为矩形的常磁化率和变磁化率二度体模型,针对本文算法的计算精度和计算效率进行了验证.模型算例结果表明:该算法计算精度高,相对误差绝对值均小于1%;计算速度快,网格节点剖分2501×2501的模型模拟时间为4.18s;适用于任意复杂地形模型.

摘要：提出一种基于积分解的空间-波数混合域二度体磁异常数值模拟快速算法。该方法将磁异常二维空间域卷积问题,通过傅里叶变换转换为空间-波数混合域垂向一维积分问题,将一个复杂问题分解为多个小问题,不同波数的小问题之间具有高度并行性;保留深度方向为空间域,采用二次插值的形函数计算垂向一维积分,便于浅层网格适当加密,深层网格适当稀疏,兼顾计算精度、计算效率及模拟复杂地形。在此基础上,根据一维形函数积分的特点,提出了一种适用于起伏地形条件下的磁异常快速计算方法,核心思想是对于相同的单元积分进行存储,避免重复计算,进一步提高了计算效率,尤其适用于复杂地形条件下的模拟。模型算例中分别设计了突变介质模型、起伏地形模型和复杂模型,通过数值解与解析解对比结果表明:本算法正确、可靠,且具有数值模拟精度高、计算速度快,适用于任意复杂地形的特点。

摘要：针对任意密度分布的二度体重力异常的计算,提出了一种空间-波数混合域二度体重力异常正演方法。该方法将二度体引力位满足的偏微分方程通过一维傅里叶变换,转化为不同波数相互独立的常微分方程,并采用基于二次插值的一维有限单元法求解该常微分方程。该方法充分利用了傅里叶变换的高效性和垂向网格剖分的灵活性,实现了二度体重力异常的高效、高精度数值模拟。采用截面为矩形的常密度和变密度的二度体模型,对所提算法的计算精度和速度进行检验,结果表明,与传统的二度体正演方法相比,算法在保证计算精度的同时,提高了计算效率。