(完整word版)重力数据处理过程

数据处理与异常推断解释

一、数据处理方法的选择

实测的重力异常是地下由浅至深各类地质体的物性差异在地面综合叠加效

应，其中包括界面起伏，岩性不均匀等诸多地质因素在内。为了从实测异常中提取和强化有用信息，压抑干扰噪声，提高重力勘探综合地质解释的能力，故需对

实测资料进行数据处理和综合分析。

1、数据处理目的

通过不同的数据处理手段，达到突出区域重力场信息、突出与强化断裂带异常信息、突出局部重力异常信息，有效地克服或压制不同干扰异常。顺利达到完成区域重力场特征分析、提取剩余异常、断裂构造划分与分析，圈定钾矿成矿有利部位等地质任务。

2、常用的数据处理方法

数据处理采用中国地质调查局发展研究中心推广的多元信息处理系统软件—GeoExpl及中国地质大学MAGS软件进行数据处理。数据处理的目的是在消除各类误差的基础上从叠加场中分离或突出某些目标物的场，并使其信息形式（或信息结构）更易于识别和定量解释。

常用的处理方法有：各种滤波、趋势分析、解析延拓（上延和下延）、导数转换（水平和垂直导数）、圆滑（圆环法和窗口法）、多次切割、差值场法、小波多尺度分析法等方法。

（1）、数据网格化

为空间分析模块及其它数据处理提供数据源。本次采用克里格法，200米×200米，搜索半径1500米。

（2）、异常分离

采用不同滤波因子的正则化滤波、差值场法、小波多尺度分析法、向上延拓等，可分别求取“区域场”和“局部场”，达到异常分离目的。

（3）、延拓处理

向上延拓：压制了浅部小的地质体场的干扰，了解重力异常衰减规律，随着上延高度增加，突出了深部大的地质体的场。区域场反映了测区深部地质环境和

地质构造特征的差异性，为测区地质构造分区划分提供了重要信息；本次向上延拓自100 m、200 m、500 m、1000 m、2000 m，共5个高度。

向下延拓：利用向下延拓可以分离水平叠加异常。密度体埋深大，异常显得宽缓。越接近密度体，异常的范围越接近其边界。本次向下延拓自100 m、200 m、300m、500 m四个高度。

（4）、水平方向导数及水平总梯度

为了准确划分断裂构造，可求取不同方向的水平方向导数、水平总梯度，以及必要时进行“线性增强”处理。

△gu＝(Vxz2+Vyz2)1/2。其中Vxz是重力异常沿X方向的一阶导数，Vyz是重力异常沿Y方向的一阶导数。水平总梯度与水平方向导数结合，可以更加准确划分和解释断裂构造。

（5）、垂向导数

垂向导数不仅在局部异常分析中起重要作用，主要突出浅源异常，而且垂向二阶导数的0值区（线）与岩体边界关系密切。

（6）、小波多尺度分析法

把小波多尺度分析方法应用于重磁测资料处理，野外观测值ΔG经一阶小

波分解，得到局部场ΔG

局1和区域场ΔG

区1

，把ΔG

区1

作二阶小波分解得ΔG

局2

到和ΔG

区2，再把ΔG

区2

作三阶小波分解可得ΔG

局3

和ΔG

区3

，…,还可以继续分

解下。分解阶数视异常的特征和地质情况来决定，解释时赋于小波逼近部分和各阶的细节明确的地质意义。

根据小波多辩分析的原理，及小波细节的微分特征，实现对位场的多尺度分解及断裂分析。

根据本次1：2.5万重力调查工作的目的任务，重点在于提取可靠的局部重力低值异常，因此，在异常分离上采用多方法进行处理，对比选择抗干扰能力强的方法提取弱局部重力异常。

二、重力异常定性解释

重力异常的解释必需以地层岩石物性资料为基础，注重平面与剖面相结合，定性解释与定量解释相结合，正演与反演相结合。人们对客观事物的认识过程是一个不断实践—认识—再实践的反复过程。同样，对重力资料的处理解释亦是如

此，需要对成果进行不断修正和完善，多次反馈，以便获得一个合理的解释成果。