关于布格重力异常计算及资料处理与反演和解释的报告

关于

布格重力异常计算及资料处理与反演和

解释的报告

姓名：***

班级：061084-27

学号：**********

指导老师：***

日期：2011.4.14

目录

前言 (2)

目的 (2)

任务要求 (2)

工作过程 (2)

成果 (2)

工作内容及步骤 (3)

§1-布格重力异常计算 (3)

§2-布格重力异常处理 (3)

1.绘制平面等值线图 (3)

2.异常处理（分离区域异常和局部异常） (6)

§3-布格重力异常反演——特征点法反演 (11)

§4-布格重力异常的解释 (13)

评述与结论 (13)

评述 (13)

结论 (14)

关于布格重力异常计算及资料处理与反演和解释的报告 前言

目的：熟悉并掌握布格重力异常计算及资料处理与反演和解释 任务要求：

根据在一个地区重力测量的结果，计算出布格重力异常，并根据异常进行资料处理和解释，并完成一份工作报告。 工作过程：

(1)利用实测的相对重力值、相对高程值和X,Y 坐标值，计算各种校正（地形校正除外），纬度校正用 计算，自由空间（或高度）校正用

计算，中间层校正用 计算，已知地表物质密度为2.50g/cm 3，起算点纬度为45°；

(2)获得各点处的布格重力异常值后，绘出平面等值线图，等值线距为0.5mGal ；

(3)根据异常（平面或剖面）特征，选用适当的方法进行处理（如压制干扰、消除区域场等）进行处理，并对处理效果进行描述； (4)将处理后的异常进行反演；

(5)写出全部过程和所采用的处理与反演方法之应用理由。

成果：根据布格重力异常数据计算及资料处理与反演初步结果判断，该异

常应由地区下一球体引起，球体埋深98.8m ，剩余质量t 6

103.07⨯，球体中

心在地面的投影点坐标为（248.8,248.8）m 。

):,()2sin(814.0mk

X mGal X g ∆⋅-=ϕδϕ):,(3086.0m h mGal h g f ∆⋅=δ):/:,(0419.03m h cm g mGal h g ρρδσ∆⋅-=

工作内容及步骤

§1-布格重力异常计算

用excel先算出各项校正（除地形校正外），需要注意的是在纬度校正中为测点到总基点间纬向距离，由于测点都位于总基点以北，故取正值；为总基点纬度即45°; 单位要划为km。中间层校正中取值为2.50 。

§2-布格重力异常处理

1.绘制平面等值线图

获得各点处的布格重力异常值后，用surfer绘图软件中的自然邻点网格法绘出平面等值线图，等值线距为0.5mGal，旁边为经过九点平滑法处理后的等值线图如下：

图例：等值线图例：等值线X

∆

ϕX∆ρ

3

/cm

g

—5mGal—

—5mGal—

可以看到经过平滑后的等值线光滑了很多，相当于虑去了高频成分。平滑处理的目的是消除异常数据中由观测引起的偶然误差以及由地表附近密度分布不均匀引起的杂乱无章的重力效应，以得到有意义的地质体引起的异常。

为了更好的了解该区域重力异常特性，再用surfer 做出其一阶导平面等值线图和二阶导平面等值线图：

一阶导在（y=250）线上的剖面图 二阶导在（y=250）线上的剖面图

100

200

300

400

500

-0.04

-0.02

0.02

0.04

0.06

V xz (mGal/m)

X(m)

100

200

300

400

500

-0.1

0.1

0.2

0.3

0.4

V zzz (mGal/m 2)

X(m)

xz V zzz

V

由上两图跟以下球形正演得出的平面等值线图和剖面图相比较，可以基本判定地下引起重力异常的物体形态似球，且根据二阶导平面等值线图还可以推断下面可能是由一个大球和四个小球叠加而产生的重力异常。这为后面的反演过程选用什么模型提供了依据。

下面是球体重力异常一阶及二阶导正演剖面图：

2.异常处理（分离区域异常和局部异常）

根据异常特征，我们知道该地区异常比较稳定，高频成分较少，在前面已摆出经过九点二次平滑法压制干扰后的等值线图，前后基本上变化不大，只是等值线更加光滑滑而已；而我们看到异常的导数等值线图把局部异常突出的比较好，因此我们可以用高阶导数法来提取局部重力异常。

又我们很容易看出该地区区域重力异常呈线性递增，因此我们还可选用趋势分析（最小二乘多项式拟合）的方法进行处理消除区域场等，最后对处理效果进行描述。

鉴于对编程能力的需要，高阶导数法难度较大，短期内无法实现，借助软件的话，我们不知道其演算过程，反演引起的误差也无法估计，因此我暂时选用趋势分析法（最小二乘多项式拟合法）。

趋势分析法是用多项式拟合区域性背景场。其原理如下：

选用一个曲面函数n阶多项式，用于拟合一个区域内的异常；通过若干点上异常值代入多项式可以得到一个方程组，求解方程组可得到多项式的系数，即得到描述区域异常的多项式函数；由于多项式的阶次一般视情况取1~4阶，去拟合相对简单的异常背景，以实现区域异常与局部异常的分离。

eg.下面我先尝试着用matlab针对NE剖面（y=x）编了一个函数来做最小二乘多项式拟合法进行实验，下面是源代码：

clear,clf,clc

x=[0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20];

y=[0.125527 0.7676005 1.2901735 1.789785 2.529205 3.127242 4.11251 5.200932 6.7870845 8.761621 9.9646955 9.618232 8.793035 8.1096465 8.027642 8.1744065 8.4864405 8.804667 9.265239 9.661196 10.209];

cc=polyfit(x,y,1) %求出多项式前的各项系数，改变数字分别