重磁实验报告(地大).doc

重磁资料采集与处理实习 实习目的 （1）通过本次实习，加深对理论知识的认识和理解。 （2）熟悉Grapher和sufer以及matlab软件的使用，会进行基本的操作和数据处理。 二、实习内容 （1）重磁数据的光滑、拟合、插值和网格化 利用Grapher软件实现磁异常曲线的光滑、拟合与去噪 上图红线代表线性光滑后的结果，可见磁异常在局部呈锯齿状，很可能地下分布有基性的喷出岩；蓝线代表10阶多项式拟合后的结果，可以反映区域场的变化情况。 将原始曲线改为散点图，可看出光滑后的效果。 利用Surfer软件实现磁异常数据的网格化与显示 测区内测点分布图如下： 打开sufer，点击Grid中出现Data，然后选中目标文件进行网格化，将网格化的文件在sufer中显示如下： 组合长方体重力异常计算与分析 计算出多个长方体的重力异常，并将结果导出为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 利用Sufer软件绘制重力异常平面等值线图 Model 1: Model 2： Model 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 利用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) 上述三种组合模型，第一个图像为原始磁异常，第二个为滑动平均后的结果，第三个为得到的局部磁异常图像，从图上可以看出滑动平均的不足，即会产生虚假异常。 组合长方体模型重力异常边界分析 计算出多个长方体的重力异常，并绘制其平面等值线图 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 利用Surfer软件求水平一阶方向导数 Muti-Model_1:(Model 1+Model 2) SN first derivative: EW first derivative: NE first derivative: SE first derivative: Muti-Model_2:(Model 1+Model