重磁实验报告（地大） - 图文

重磁资料采集与处理实习

一、实习目的

（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:(Model 1+Model 3) 输出文件为muti-model_2.grd

Muti-Model_3:(Model 1+Model 2+Model 3) 输出文件为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:

2、海底地形校正，得到布格重力异常，并绘制等值线图

根据布格校正公式BA = FA - 0.0419 * ( 1.64 * Dw)（其中FA为自由空间重力异常，BA为布格重力异常，Dw为海底深度，向上为正），利用surfer中Math代入公式，输出文件BA.grd。然后绘制等值线图。

3、布格重力异常解析延拓

首先要对BA.grd做一个转换处理，变成sufer6 Text.grd ,命名为convert BA.grd。将布格重力异常导入Matlab，分别求1km，2km，5km，10km，20km，50km，并保存计算结果。 将延拓后前后的图像进行对比，第一幅为延拓前，第二幅为延拓后。 BA_UC_1.grd:

BA_UC_2.grd:

BA_UC_5.grd:

BA_UC_10.grd:

BA_UC_20.grd:

BA_UC_50.grd:

4、利用向上延拓方法，实现布格重力异常位场分离

通过3中图像的对比，可以将向上延拓5km和10km得到重力异常作为区域场，然后再求取局部异常，从而实现重力异常位场分离。下图第一个图为原图像，第二个图为通过延拓得到的背景场，第三个为局部异常。

延拓5km作为背景场：

本文来源：https://www.bwwdw.com/article/vmfv.html