重磁实验报告(地大)

重磁资料采集与处理实习

一、实习目的

（1）通过本次实习，加深对理论知识的认识和理解。

（2）熟悉Grapher和sufer以及matlab软件的使用，会进行基本的操作和数据处理。

二、实习内容

（1）重磁数据的光滑、拟合、插值和网格化

1、利用Grapher软件实现磁异常曲线的光滑、拟合与去噪

上图红线代表线性光滑后的结果，可见磁异常在局部呈锯齿状，很可能地下分布有基性的喷出岩；蓝线代表10阶多项式拟合后的结果，可以反映区域场的变化情况。

将原始曲线改为散点图，可看出光滑后的效果。

2、利用Surfer软件实现磁异常数据的网格化与显示

测区内测点分布图如下：

打开sufer，点击Grid中出现Data，然后选中目标文件进行网格化，将网格化的文件在sufer中显示如下：

（2）组合长方体重力异常计算与分析

1、计算出多个长方体的重力异常，并将结果导出为GRD格式Model 1：

X1 = -100; %长方体X方向起点坐标

X2 = 100; %长方体X方向终点坐标

Y1 = -100; %长方体Y方向起点坐标

Y2 = 100; %长方体Y方向终点坐标

Z1 = 10; %长方体Z方向起点坐标

Z2 = 55; %长方体Z方向终点坐标

经过matlab运行后导出mod_1.grd

Model 2：

X1 = 120; %长方体X方向起点坐标

X2 = 180; %长方体X方向终点坐标

Y1 = 120; %长方体Y方向起点坐标

Y2 = 180; %长方体Y方向终点坐标

Z1 = 1; %长方体Z方向起点坐标

Z2 = 20; %长方体Z方向终点坐标

经过matlab运行后导出mod_2.grd

Model 3：

X1 = -75; %长方体X方向起点坐标

X2 = -125; %长方体X方向终点坐标

Y1 = -75; %长方体Y方向起点坐标

Y2 = -125; %长方体Y方向终点坐标

Z1 = 1; %长方体Z方向起点坐标

Z2 = 20; %长方体Z方向终点坐标

经过matlab运行后导出mod_3.grd

2、利用Sufer软件绘制重力异常平面等值线图Model 1:

Model 2：

Model 3：

3、合并生成多个长方体组合模型的重力异常

利用surfer中grid中的math进行组合。

Muti-Model_1:(Model 1+Model 2)

输出文件为muti-model_1.grd

输出文件为muti-model_2.grd

输出文件为muti-model_3.grd

4、利用Surfer滑动平均方法，分离大小长方体异常

首先点击grid，然后选择filtering中的moving-average，用滑动平均后的结果减去输出结果，即可得到。

Muti-Model_1:(Model 1+Model 2)

Muti-Model_2:(Model 1+Model 3)

Muti-Model_3:(Model 1+Model 2+Model 3)

上述三种组合模型，第一个图像为原始磁异常，第二个为滑动平均后的结果，第三个为得到的局部磁异常图像，从图上可以看出滑动平均的不足，即会产生虚假异常。

（3）组合长方体模型重力异常边界分析

1、计算出多个长方体的重力异常，并绘制其平面等值线图

Muti-Model_1:(Model 1+Model 2)

输出文件为muti-model_1.grd，平面等值线图如下：

Muti-Model_2:(Model 1+Model 3)

输出文件为muti-model_2.grd

Muti-Model_3:(Model 1+Model 2+Model 3)

输出文件为muti-model_3.grd

2、利用Surfer软件求水平一阶方向导数

Muti-Model_1:(Model 1+Model 2)

SN first derivative:

EW first derivative:

NE first derivative:

SE first derivative:

SN first derivative:

EW first derivative:

NE first derivative:

SE first derivative:

SN first derivative:

EW first derivative:

NE first derivative:

SE first derivative:

Muti-Model_1:(Model 1+Model 2)

Muti-Model_2:(Model 1+Model 3)

Muti-Model_3:(Model 1+Model 2+Model 3)

（4）组合球体磁异常分析

假设组合球体模型位于武汉地区（磁化倾角45°，磁偏角0°），地面水平，且不考虑剩磁

假设有这样三个磁性球体：

---磁性球体1

--X点坐标0 m

--Y点坐标0 m

-- 中心埋深50 m

-- 半径20 m

-- 磁化率0.2

---磁性球体2

--X点坐标-30 m

--Y点坐标-30 m

-- 中心埋深10 m

-- 半径 3 m

-- 磁化率0.2