重磁数据处理大报告-陈亮

中国地质大学（武汉）地空学院

姓名：陈亮

班级： 061132

学号： 480

指导老师：杨宇山

目录

一、地质任务3

二、工区概况3

三、数据整理4

一、重力资料数据整理4

二、磁场资料数据整理6

四、材料图4

五、研究区重磁异常分析10

六、重磁资料数据处理13

1、重力场延拓13

2、磁场化极处理 16

3、重力场的分离 17

4、磁场的分离18

5、重磁资料导数换算处理20

七、局部重磁异常分析25

八、学习总结25

一、地质任务

（1）将布格重力异常Δg和磁异常ΔT整理出来，计算布格重力异常和磁异常的总精度。

（2）利用surfer绘制测点点位图（即实际材料图），布格重力异常平面图，磁异常ΔT平面图。

（3）根据密度统计表分析研究区的物性特征。

（4）分析研究区重磁异常特征。

（5）对重磁资料进行处理（化极、延拓、导数换算等并绘制结果图件），并进行断裂构造分析。

（6）提取与矿有关的局部重磁异常（绘制结果图件），并进行对应分析，区分矿与非矿异常、磁铁矿与磁铁矿的可能分布范围。

（7）撰写报告。

二、工区概况

研究区位于我国中东部地区，地理坐标为东经°—°，北纬°—°，处在我国非常重要的铁多金属矿成矿带西段。在以往地质、物探工作基础上，2015年3月人们在研究区中部完成了面积为5km²（×2km，线距50m，点距20m，测向方位角0度）的1:5000地面重磁扫面工作。

此次重力施工设计精度为50μGal，磁测施工设计精度为5nT，共完成了3116个测点，检查点159个，重力观测误差为μGal，磁测观测误差为；重力近区地改范围0~20m，在野外完成，采用差分GPS（RTK）进行8方位方形域测量，检查点59个，误差为μGal。点位测量采用RTK差分GPS进行测量，检查313个点，高程测量误差为，平面位置测量误差为。

研究区铁矿赋存于燕山期早的中酸性岩与三叠系地层的接触部位，研究区经历了后期的构造变动，断裂构造发育，浅表磁铁矿经历了风化和淋滤作用后，形

成了50%以上的高品位赤铁矿。

工区地层

工区全为第四系覆盖。在研究区西南方零星有岩体岀露，侏罗系、三叠系、石炭系等地层在工区东北方有小面积岀露。

本次施工共采集了352块各类岩石标本，经过测量得到如下表1物性参数。

表1 岩石物性标本测定结果统计表

三、数据整理

一、重力资料数据整理

利用实测的相对重力值、相对高程值和X,Y 坐标值，计算各种校正。 纬度校正用):,()2sin(814.0mk X mGal X

g ∆⋅-=ϕδϕ计算，

自由空间（或高度）校正用):,(3086.0m h mGal h g f ∆⋅=δ， 中间层校正用):/:,(0419.03m h cm g mGal h g ρρδσ∆⋅-= 地形校正分近区地改和中远区地改：

近区地改：12/ d i i d ∆⨯=∆ρ，△d1:表示密度为2g/cm 3时的近区地改值，原数据中已给出。

中区：20-500米，所选网格为20；密度为2g/cm 3。

远区：500-10000米，所选网格200。密度为2g/cm 3。

其他密度中区地改：22/ d i i d ∆⨯=∆ρ，△d2:表示密度为2g/cm 3时的中区地改值，由软件计算。

其他密度远区地改：32/ d i i d ∆⨯=∆ρ，△d3:表示密度为2g/cm 3

时的远区地

改值，由软件计算。

布格重力异常=相对基点总重力值+所有校正

密度

(kg/cm3)

相关系数

2

3

当密度为cm3时，布格

异常与高程的相关性为零，所以中间层和中远区地改采用的密度为cm3。由于本次实习采用的模型问题，最终还是采用cm3的密度作为改正密度

根据布格异常计算公式

在excel 中计算得到布格异常。

现已知重力施工设计精度要求为50μGal ，重力观测误差为μGal ，重力近区地改误差为μGal ；高程测量误差为，平面位置测量误差为。

高程测量误差引起的重力误差为：{Δgh}.={h}m

σϕδδδδg g g g g f B +++∆=

计算得到误差为**10^-4=μGal

平面位置测量误差引起的重力误差为：{Δgh}.=*sin(2θ)*{ Δx}m 计算得到误差为*sin(2*°)*= μGal 所以，中远区地改的误差范围是：0~μGal

在中区地改中，计算原测区的改正值，将测区的点分别移动半个网格后，在计算其地改值便可以得到新的改正值，用原改正值与新改正值计算中区地改的误差，公式如下：

计算得到中区地改的精度为：μGal 同理，远区地改的精度为： μGal 综上，布格异常的总精度为： ++ +++=μGal ＜50μGal 故满足施工要求。

二、磁场资料数据整理

利用实测的磁测总场、相对高程值和X,Y 坐标值，计算各种校正。 日变校正：由软件完成0δ。

正常场随高程变化校正：1000/9.221h δδ⨯=。 正常场随维度变化校正：1000/)(75.520Y Y -⨯-=δ。 总的磁异常：210δδδδ++=T 。 处理后的磁测数据

四、材料图

实际地形图

实际测点图重力

重力野外资料

布格重力异常磁场

磁场野外数据

磁异常图

五、研究区重磁异常分析

重力异常