重力实习报告

由于校园内实测仅限于一条测线，所以我们运用二连浩特的实测数据来实现向上延拓的处理方法，二连浩特原始异常图如下所示：

图7 二连浩特原始异常图

通过matlab处理，将上述异常向上延拓500m得图8和图9所示的异常。

图8 二连浩特原始异常上延500m异常图

- 9 -

图9 二连浩特原始异常上延500m剩余异常图

2、导数换算的原理和方法

不同形状的地址体的重力异常的倒数具有不同的特征，可以通过特征不同对异常进行解释和分类；重力异常的倒数可以突出浅而小的地质体异常而压制区域性深部的地质重力效应，在一定程度上分离深浅地质体，且倒数次数越高，这种分辨能力越强；重力高阶导数可以将几个互相靠近、埋深差不多的相邻地质体引起的叠加异常分离出来。

倒数换算的计算公式 总水平梯度： 22gh?gx?gy

222gx?gy?gz总梯度模： ga?

???1?gz?gt?tanTilt梯度： 22??g?gxy??

同样对二连浩特数据进行处理，运用matlab计算得的导数换算图如图10、图11所示。

图10 导数换算1

- 10 -

图11 导数换算2

3、DEXP法场源位置估计的原理和方法

对位场f进行延拓高度为z的频率域向上延拓，得到延拓后的重力异常f(z)。对f(z)作n阶DEXP变换，得到尺度场。尺度场Wn(z)极值点对应于重力场源深度，因此DEXP法反演结果反映了地下场源大致位置。取不同的延拓高度，重复(1)和(2)步骤，就得到位场的DEXP法反演结果。

运用matlab计算得测线B的场源位置图如图12所示。

图12 测线B的DEXP法场源位置图

4、人机交互正反演拟合的原理和方法

人机交互式反演拟合即通过人为改变地质体的模型和计算机的计算，使计算异常曲线逐步接近测量曲线，从而反演出与实际地质体相近的模型，这种方法与最优化选择法基本相同。其计算过程如下：

（1）建立直角坐标系，横坐标沿测线方向，垂直向上的坐标表示重力异常，垂直向下的坐标表示深度。观测重力异常一般为一些离散分布的点值。

- 11 -

（2）根据重力异常的形态及已知的地质资料在坐标系中，用光笔画一个初始模型，例如一个矩形，这时在上界面上出现了由这个矩形引起的理论异常曲线。

（3）根据理论异常与观测异常的差异，增加多边形截面的角点数，并调整角点的深度及水平位置，使理论曲线与观测异常点吻合得最好。在这种情况下得到的多边形截面的形状就是地质体形状的最佳解答

图13 人机交互反演结果图

- 12 -