重力异常物性正则化反演研究

周宇轩 雷宛

经典位场物性反演是建立在数与模型之间的线性关系,利用离散正演公式将每个剖分单元相对观测点的核函数进行计算,以寻求观测数据与核函数有最好的相关性。该类方法在计算过程最大的问题就是在解巨大超定方程组时,导致求解极不稳定。在前人的研究基础上,提出一种基于先验密度模型约束正则化条件下的重力异常视密度求解方法,通过L曲线法选取最佳的参数,使求解过程中有较高的稳定性和求解精度。这里以一理论模型作为观测数据,并引入噪声作为观测误差,通过算法的对比表明,该方法有很高的计算精度和稳定性,对于需要求解大型超定矩阵的三维模型,在给定先验初始模型的条件下,该方法仍然能够取得很好的效果。