王正科，卢琳，刘花婷

(陕西地矿第二综合物探大队有限公司,陕西 西安 710016)

0 引言

实测重力异常包含了从地表(观测面)到深部所有密度不均匀引起的重力效应，即是观测面以下所有不均匀密度体的叠加异常。在解译某个(类)目标地质体所产生的重力异常时，必须从实测重力异常中分离出单纯由该地质体引起的异常，再对这个异常进行反演获取目标地质体形体、规模等参数。然而，要把单个地质体所引起的重力异常从实测重力异常中分离出来还是非常困难的，因此实测重力异常的分离就是地质解释工作中的重中之重。

对于实测重力异常的复杂性[1]，划分区域(场)异常和局部异常的方法多种多样，诸如：小波多尺度分析法、补偿圆滑滤波法、正则化位场分离法、圆周平均法、滑动窗口平均法、插值切割法、滑动趋势分析法等。这些方法均是对已经给定的场进行分离，每一种方法都有其局限性和不足之处，其中滑动窗口平均法的不足之处为：① 正方形窗口过大达不到分离局部异常的目的；过小则区域场中保留许多局部异常；② “虚假异常”是方法的固有原因，不可避免，使用此方法需要选择不同参数而做出多种相关图件，并需要在解释中相互参照对比和补充。

滑动趋势分析法的不足之处则是多项式的次数选择问题，次数高则会导致拟合异常的震荡效应加剧，叠加的局部异常趋势不明显。一般选取“1、2、3”次多项式，按照“最小二乘法”进行回归分析，取“最佳”(窗口多项式拟合值与节点异常值之差值最小)作为网格节点的局部异常值。……