基于高阶多项式密度函数的重力反演

摘要

重力反演是研究地下密度体空间展布的地球物理手段之一.现有的密度反演方法通过对地下介质的网格剖分,直接反演每个网格单元的密度值.本文用多项式函数近似表示地下复杂的密度变化,提出一种通过反演多项式密度函数的系数确定密度变化的新方法.不同于常规网格剖分方法,该方法无需在垂向上划分单元,在一定程度上缓解了网格疏密、内存占用与反演精细度之间的矛盾.理论模型测试表明多项式系数反演与多种约束相结合,能清晰地突出局部异常体的位置、规模和边界等信息,且效果优于L2模约束反演的结果.该方法成功应用于济阳坳陷区古潜山和凹陷的识别,根据反演的密度差剖面可大致确定不同岩性密度体的边界,弥补了地震剖面反映古潜山展布的不足.