基于N阶多项式密度函数的重力异常反演方法

摘要

本发明涉及一种基于N阶多项式密度函数的重力异常反演方法，属于地球物理勘探技术领域。其包括如下步骤：第一步,对实测重力异常进行预处理，得到研究目标体产生的重力异常；第二步,单元划分,即将研究区划分为竖向并排的N_r个矩形单元，各个矩形单元纵向长度相同，且大于或等于研究目标区的深度；第三步，设定密度多项式函数的最大阶数，构建重力反演的核函数矩阵；第四步，构建反演目标函数，建立相应的重力反演线性方程组，求解方程，使计算的理论重力异常与实测重力异常残差在允许的范围内，得到相应各矩形单元密度函数的多项式系数；第五步，利用反演的多项式系数计算整个研究区的地下密度分布。本发明实现了对研究区重力异常的反演。