郭良辉，孟小红，石 磊，陈召曦

1 地下信息探测技术与仪器教育部重点实验室（中国地质大学，北京），北京 100083

2 中国地质大学（北京）地球物理与信息技术学院，北京 100083

3 中国地震局地球物理研究所，北京 100081

1 引 言

根据重力观测值整理得到的重力异常，包含了从地表到深部所有密度不均匀体引起的重力效应，信息非常丰富.实测重力异常是所有这些重力效应的叠加.要根据实测重力异常反演解释某个地质体（比如岩矿体、地质构造或物性分界面），必须首先从叠加异常中分离出单纯由这个地质体引起的异常，然后用这个异常进行反演解释.然而，从叠加异常中分离出某个地质体引起的异常，或者把叠加异常分解为几个地质体引起的单一异常，是比较困难的.重力异常分离是目前重力资料处理解释中没有很好解决的难题之一.

当前，重力异常分离方法众多，如最小二乘平滑法、带通滤波法、向上延拓法、匹配滤波法［1］、补偿圆滑滤波法［2］、正则化滤波法［3］、维纳滤波法［4-5］、小波变换法［6-8］、熵滤波法［9-10］、非线性滤波法［11］等.这些方法均有各自的特点和优势，但同时又有各自的局限性.

传统向上延拓法存在的一个主要问题是在进行向上延拓时，浅源短波长信息和深源长波长信息同时延拓，在浅源信息被压制的同时，深源信息也受到压制，将实测异常减去上延后的区域异常所得到的剩余异常中，仍包含有区域成分，异常分离不彻底［12-13］.针对该问题，Pawlowski［12］根据维纳滤波和……