地球探测科学与技术学院
沈阳及其附近地区重磁数据处理与解释
报告
姓名：李雪垒
学号：
班级：四班
专业：勘查技术与工程（应用地球物理）
指导教师：吴燕冈教授
目录
前言 (2)
第一章重磁数据处理基本原理与方法 (3)
一、重力场与磁场的波谱介绍 (3)
二、数据处理的基本方法 (3)
三、Surfer、Grapher简介 (3)
第二章地质概况 (5)
一、东北及其附近地区地质概况 (5)
二、实验区内的地质概况 (5)
第三章区内重磁异常综合解释 (8)
一、重力数据异常处理与解释 (8)
二、磁异常数据异常处理与解释 (12)
三、重磁异常场综合分析 (15)
第四章本次实验的初步结论 (16)
主要参考文献 (16)
前言
重力勘探是测量与围岩有密度差异的地质体在其周围引起的重力异常﹐以确定这些地质体存在的空间位置﹑大小和形状，从而对工作地区的地质构造和矿产分布情况作出判断的一种地球物理勘探方法。磁法勘探是通过观测和分析由岩石、矿石（或其他探测对象）磁性差异所引起的磁异常，进而研究地质构造和矿产资源（或其他探测对象）的分布规律的一种地球物理勘探方法。二者有广泛的应用，如研究地壳深部构造；研究区域地质构造，划分成矿远景区；掩盖区的地质填图，包括圈定断裂﹑断块构造﹑侵入体等；广泛用于普查与勘探可燃性矿床(石油﹑天然气﹑煤)；查明区域构造，确定基底起伏，发现盐丘﹑背斜等局部构造；普查与勘探金属矿床(铁﹑铬﹑铜﹑多金属及其他)，主要用于查明与成矿有关的构造和岩体，进行间接找矿；也常用于寻找大的﹑近地表的高密度矿体，并计算矿体的储量；工程地质调查；如探测岩溶，追索断裂破碎带等。
随著电子技术的发展和微处理机的广泛应用，测量磁场3个分量及其梯度的高精度航空磁力仪已经制成。加上高精度的导航和数据处理，绘图和资料解释推断的自动化，今后航空磁法勘探将代替部分地面磁法勘探，并在工作过程中自动作出解释，绘出磁性体空间分布图。利用这些图件，再结合其他资料，能可靠地对工作地区的地质构造作出推断，供找矿﹑找地下水﹑工程建设和地震预报等方面应用。我国在改革开放以后，随着科学技术的飞速发展，在重磁勘探领域取得了令人瞩目的成就，在测量精度方面大大提高。
由于重磁法勘探应用广泛，成本不高，因此在勘探领域一般是其他勘探方法之前的首选方法。由于地球区域复杂，通常要对所采集的数据进行各种处理，以去除各种无关影响，提取所要的结果。同时根据处理结果对其进行解释，其中解释又分为定性解释与定量解释，其处理方法与解释方法在本次实验中均有所简单涉及。
在本次课程中，我们主要学习了重磁异常的空间域处理与转换，重磁异常的波数域处理，重磁异常的反演方法以及重磁资料的地质解释和在勘探中的应用等。
基于本学期学习的内容及理论知识，结合自己的理解对东北地区的重磁异常做初步的处理及解释。本次实验作图工具使用的为Surfer和Grapher，同时也参考了一些前人的研究成果，以此作为基本出发点，进而得出一些初步的结论。
第一章重磁数据处理基本原理与方法
一、重力场与磁场的波谱介绍
重磁数据的处理分为空间域和波数域，在空间域中重磁数据的处理十分复杂，而经过傅里叶变换到波数域中后就变得十分简单方便。在波数域中对数据所有的处理都是一些因子的乘积；而且由于快速傅里叶变换的出现运算速度明显加快；除此之外波数域处理的另一个优点是没有了边缘损失。鉴于上述优点重磁数据处理都是在波数域中进行的。
在波数域数据处理时要先把数据扩充到2的整数幂倍，然后利用FFT计算原始数据的谱，再将原始异常的谱乘以处理过程的权因子求转换谱，最后用逆FFT 转换回空间域异常。用到泊松公式波数谱，其中重磁异常波数谱由水平尺寸、深度因子、位移因子和磁化因子乘积构成。
二、数据处理的基本方法
1、解析延拓：
根据某观测平面上的实测异常，换算场源以外其它空间位置的异常称为解析延拓。换算平面位于实测平面之上，称为向上延拓，换算平面位于实测平面之下，称为向下延拓。向上延拓可以突出深部异常，压制浅部异常；计算形状参数勾绘断面图。而向下延拓可以突出局部异常，压制深部异常；划分水平叠加异常；评价低缓异常；增加解释推断的信息。延拓的波数谱只与深度因子有关，h0为向上延拓，h0为向下延拓。
下面以无限长直立台阶为模型说明解析延拓的作用：
模型参数：埋深H=10km ，h=8km ，台阶厚度t=H-h=2km ，剩余密度ρ=cm3，剖面线从x从-100km到+100km ，y为1到128km，间隔均为1km 。
图1-1（黑色为原始数据，蓝色向下延拓2km，红色为向上延拓2km）
由图1-1对比可知，向上延拓后重力异常幅值变小了，向下延拓后重力异常幅值
变大了，当在边界时会出现局部振荡。
2、导数异常的计算：
（1）垂向导数的计算：求垂向导数就是要根据实测重磁场T 的分布来计算n n T z ∂∂在该平面上的分布。n 为求导的阶数，常用的是n=1和n=2，即计算垂向
一阶导数或垂向二阶导数。因为n n n T r T z ∂=∂，所以n
r 是计算年n 阶垂向导数的波数响应。即垂向一阶导数的波数响应为r, 垂向二阶导数的波数响应为r 2。垂向导数异常能区分相邻异常体的异常，减少其相互迭加的影响。或者减轻围岩的干扰，分离迭加在背景场中的局部场。从理论上分析用垂向二阶导数的零值线可以圈定异常体韵范围和位置。
（2）水平导数的计算：水平导数即对x 方向的导数T
iuT x
∂=∂和y 方向的导数 T ivT y ∂=∂若s 是实测平面上某一任意方向，它与x 方向的夹角为
，则cos sin (cos sin )T T T i u v T s x y αααα∂∂∂=+=+∂∂∂。方向导数可以用来突出某一方向
的异常特征。x 方向的方向导数对于x 方向的高频成分有放大作用，y 方向的方向导数对于y 方向的高频成分有放大作用，因此实践中常用方向导数来分析区内某一方向的构造线特征。
3、区域场与局部场的分离——匹配滤波：
区域场与局部场的分离是区域重磁场数据处理的一个重要方面。在实际资料解释中也有重要意义。区域场以低频成分为主，局部场则以高频成分为主。采用提取不同波数成分的场就可以完成场的分离。
匹配滤波：当实测场是由区域场与局部场迭加而成时，从径向平均对数功率谱曲线土就会显示两个拟合直线段，它们分别反映了区域场和局部场的主频段。因而采用一般的波数滤波方法就可以进行场的分离。
由于在本学期的课程中对这方面的内容涉及较少，并且在实验中也并未尝试，因此此种方法在本次报告中就不再详细介绍，在数据处理时暂时不做此方法的处理。
4、化磁极：
化极是将斜磁化的Z a 、△T 转换成垂直磁化的垂直磁异常。它消除由于磁化场的倾角和偏角引起的磁异常不对称性的一种滤波技术，以便使异常似乎是在磁极处获得的，从而异常形态简化，利于推断解释。由于在本学期的课程中并未对该种方法做详细的介绍，因此在以下数据的处理中不再应用。
三、Surfer 、Grapher 简介。
Surfer 主要用于画等值线图和三维立体图，画等值线图Countour Map 和三维立体图surface map 的文件必须是grd 文件，画图框Base map 要用bln 文件，
