重力异常 matlab,球体重力异常

《工程与环境地球物理》

实验四

一、 实验目的和要求

1. 计算球体重力异常，并分析其重力异常的特点，加深理解，通过绘图了解影响重力

异常大小的因素；

2. 掌握用Matlab 绘制重力异常曲线的方法，提高简单程序的编辑能力。

二、 实验软件

1. Matlab 编程软件

三、 实验内容

1. 计算球体的重力异常，给出平面图、主剖面图和任意一条侧剖面图；

计算内容包括重力异常、重力异常的二阶导数以及三阶垂直导数；

2. 在埋深固定情况下，计算异常最大值与球体大小关系曲线；

3. 在球体大小固定情况下，计算异常最大值与理深的关系曲线；

四、 结果分析

通过实验，得出成果图，进行分析：

由以上程序得到计算结果图，结合理论知识可知：球体为三度体，在主剖面图中重力异常在球体正上方最大，向两侧逐渐减小，侧剖面上，形态大体相同；以过球体正中心的测线垂线方向处的重力异常值最大，两侧慢慢减小。平面图上表现为以球体在地面上的投影点为中心的一系列的同心圆。

等等

五、 提交报告

纸质打印版，内容包括程序、结果图、以及结果图分析，交到地质宫315，曾老师 ps ：

球体重力异常的计算公式为：

3332222222

252222

52222

2225222222272222

4()

3[()()()]()3()3()2()(233)()r Gh Gm z h g x y z h x y h hx Vxz Gm x y h hy Vyz Gm x y h h x y Vzz Gm

x y h h x y h

Vzzz Gm x y h πρξη--?==-+-+-++-=++-=++--=++--=++剩余