重力异常正演matlab程序,圆锥体完全均衡下重力异常正演 [MATLAB]

在完全均衡的模型下,若地表有一圆锥体(山峰等),计算跨越山顶的截面上所得到的各种重力异常。

地壳密度 $kg\cdot m^{-3}$

上地幔密度 $g\cdot cm^{-3}$

地表地形圆锥体半径 (km)

地表地形圆锥体高度 (km)

计算莫霍面形变圆锥半径 (km)

计算莫霍面形变圆锥高度 (km)

地壳厚度 (km)

$2.8\times 10^{3}$

$3.5\times 10^{3}$

$2.0$

$2.0$

$2.0$

$8.0$

$30.0$

计算结果如下。横坐标单位:m,纵坐标单位:mGal

6c2ef6efc2130d4ea0998391b234d36a.png

MATLAB代码如下:

% 生成地下体的布格重力异常

syms r a z x h;

f = r / sqrt((z + h)^2 + r^2 + x^2 - 2*r*x*cos(a))^3;

dense = - 700; G = 6.67E-11;

depth = 30000; subheight = 8000; height = 2000;

total_spots = 81;

total_anom = zeros(1, 81);

total_xvec = zeros(1, 81);

spots = 20; from = 50000; to = 2500; interval = -2500;

xvec = from:interval:to;

anom = zeros(1, spots);

for no = 1:spots

rad = from + (no - 1)*interval;

gap = 800;

N = subheight/gap;

anomaly = 0;

for n = 0:N-1

Radius = 2000 - (gap/4)*n;

fc = subs(f, [z, x, h], [depth + gap*n, rad, 0]);

func =

  • 2
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是一个简单的立方体重力异常正演MATLAB代码,供您参考: ```matlab % define constants G = 6.674e-11; % gravitational constant rho = 2670; % density of the cube L = 1000; % length of the cube x = -5000:100:5000; % x-coordinates y = -5000:100:5000; % y-coordinates z = 0:100:10000; % z-coordinates [X,Y,Z] = meshgrid(x,y,z); % create the grid % calculate gravity anomalies gx = zeros(size(X)); gy = zeros(size(Y)); gz = zeros(size(Z)); for i = 1:numel(x) for j = 1:numel(y) for k = 1:numel(z) % calculate the distance between the grid point and the cube r = sqrt((X(i,j,k)^2) + (Y(i,j,k)^2) + (Z(i,j,k)^2)); % calculate the gravitational attraction gx(i,j,k) = gx(i,j,k) + (G*rho*L*(X(i,j,k)/r^3)); gy(i,j,k) = gy(i,j,k) + (G*rho*L*(Y(i,j,k)/r^3)); gz(i,j,k) = gz(i,j,k) + (G*rho*L*(Z(i,j,k)/r^3)); end end end % plot the results figure; subplot(1,3,1); imagesc(x,y,squeeze(gx(:,:,end))); colorbar; title('gx'); subplot(1,3,2); imagesc(x,y,squeeze(gy(:,:,end))); colorbar; title('gy'); subplot(1,3,3); imagesc(x,y,squeeze(gz(:,:,end))); colorbar; title('gz'); ``` 这个代码用了一个立方体模型,以100m为间隔生成x、y、z三个方向的网格,然后计算每个网格点上的重力异常。最后,将计算结果通过三个子图的方式展示出来,分别对应x、y、z三个方向上的重力异常。需要注意的是,这个模型非常简单,只是为了演示立方体重力异常正演的基本过程。如果您需要更复杂的模型或更多的计算参数,可以根据自己的需求进行修改。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值