重力异常正演实验.pdf

吉林大学地球探测科学与技术学院

实验报告

重力与磁法勘探

课程名称

成绩评定：

报告人王浩然9

班号

(姓名、学号)

2015/11/2地质宫应用地教师：

实验时间地点

球物理实验室

简单条件下规则圆柱体重力异常的正演计算

实验题目年月日

要求学生熟悉计算机常用的编程语言，能够编制简单的计算程序。学习

和掌握简单条件下水平圆柱体的重力异常计算方法，并能使用常用的绘

实验目的

图软件对所计算的结果绘制剖面和平面图，以便加深对圆柱体重力异常

及要求

特征的认识。

实验内容(包括：实验内容、实验步骤、问题讨论、结论及建议)：

实验步骤：

1.教师将本班学生分为4组。

2.使用水平圆柱体进行编程计算；

3.将计算结果绘制成相应的平面、剖面图；

4.讨论实验中遇到的以及对水平圆柱体重力异常理解过程中产生的问题；

5.对实验的认识、讨论和建议。

实验内容：

实验的模型和参数：

无限长水平圆柱体，线密度λ=1，埋深10m，20m

程序原代码：

clear

x=-105:1:105;

y=-105:1:105;

[X,Y]=meshgrid(x,y);

G=6.67;

lam=1;

h=10;

g=2*G*lam*h./(X.^2+h.^2);

vxz=4*G*lam.*h.*X./((X.^2+h.^2).^2);

vxz=-4*G*lam.*h.*X./((X.^2+h.^2).^2);

vzz=2*G*lam.*(h.^2-X.^2)./((X.^2+h.^2).^2);

vzzz=4*G*lam.*(h.^2-3.*X.^2)./((X.^2+h.^2).^3);

figure(1);

contour(g,20);title(水平圆柱体δg等值线图);

figure(2);

subplot(2,2,1),plot(x,g(106,:)),axis([-105,105,0,1.5]),title(水平圆柱体δg);

subplot(2,2,2),plot(x,vxz(106,:)),axis([-105,105,-0.1,0.1]),title(水平圆柱体vxz);

subplot(2,2,3),plot(x,vzz(106,:));,axis([-105,105,-0.05,0.15]),title(水平圆柱体vzz);

1

subplot(2,2,4),plot(x,vzzz(106,:));,axis([-105,105,-0.001,0.003]),title(水平圆柱体vzzz);

figure(3);

surf(g);title(δg);

G=6.67;

lam=1;

h=20;

g1=2*G*lam*h./(X.^2+h.^2);

figure(4);

plot(x,g(106,:),g),axis([-105,105,0,1.5]),title(水平圆柱体δg（当h变化的时候）),

Legend(h=10);

holdon;

plot(x,g1(106,:),r),legend(h=20)

;

2

3