重磁工作介绍_2010_1_17_简化

常的影响被削弱了年, ‘ 我们用上述方法编制了 , 币机算法语言程序 , 并作了模型试验 , , 后来又处理了实际资料闭有一定效果。

原始数据中存在有局部异常成分也是造成虚假异常的原因之一数据点随机分布叠加异常的划分趋势分析方法去处理规则格网数据吼 , 所以可以用不规则测网。

‘’ 至今仍是一个没有很好解决的问题、 , 还有许多工作要做二异常的反演 , 异常的反演一多模型最优化选择法、即根据观测异常求地质体的位置。

大小、产状和物性参数数据处理和解释工作中的一个不可缺少的环节欲反演的目标大体上有三种—是重磁资料求地质体用规则几何形体近似 , 的几何参数和物求一定范围物。

性参数性的空间分布的实质在于 , 。

求一定范围的物性分界面 , 反映某一地质层位的起伏目前采用的反演方法可分为两类直接法 , 。

根据观测整理的异常直接曲线拟合一选择法选择法 , 计算地质体的某些参数多用于解释简单的异常一次完成。

将实测异常曲线与一系列已知形状模型体产生的理论异常曲线进行比较当实测曲线与某一理论曲线符合为实际地质体的近似结果。

在给定的误差范围内选择法需迭代完成。

时 , 就将该理论曲线所对应的模型体作这里先介绍多模型最优化选择法在一个矿区、的应用情况 , 然后讨论此法的应用条件 , 年 , 武汉地质学院磁法组应用长方体组合模型采用了十五个模型 , 改进的马奎特法。

〕冀东对。

‘ 一区的地磁异常进行了反演川得到了各模型体的参数的理论曲线。

图表示反演得到的模型体的平面位置和由它算出理论曲线与实测曲线有些模型体的参数也。

一拟合得较好 , , 滋儡节静之火浓姆丫一之文侧气二么爹又一夕‘ 之毛‘ ‘’ 币—、、比较接近附近钻孔中的见矿情况月又卜女例如间见到体 , 孔位于第块之间 , 块和第米剩图一汤火在井深入米赤铁矿和 , 米磁性矿米磁块和第按计算结果是该处应有又如在第一、 , 性矿体块之间的孔 , 按反演推米磁性矿米到米米磁铁算在体 , 米左右应有实际钻探结果在当然此外 , 之间见到了矿区。

重磁资料采集与处理实习一、实习目的（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.grdModel 2：X1 = 120; %长方体X方向起点坐标X2 = 180; %长方体X方向终点坐标Y1 = 120; %长方体Y方向起点坐标Y2 = 180; %长方体Y方向终点坐标Z1 = 1; %长方体Z方向起点坐标Z2 = 20; %长方体Z方向终点坐标经过matlab运行后导出mod_2.grdModel 3：X1 = -75; %长方体X方向起点坐标X2 = -125; %长方体X方向终点坐标Y1 = -75; %长方体Y方向起点坐标Y2 = -125; %长方体Y方向终点坐标Z1 = 1; %长方体Z方向起点坐标Z2 = 20; %长方体Z方向终点坐标经过matlab运行后导出mod_3.grd2、利用Sufer软件绘制重力异常平面等值线图Model 1:Model 2：Model 3：3、合并生成多个长方体组合模型的重力异常利用surfer中grid中的math进行组合。

中国地质大学（武汉）地空学院姓名：陈亮班级： 061132学号： 480指导老师：杨宇山目录一、地质任务3二、工区概况3三、数据整理4一、重力资料数据整理4二、磁场资料数据整理6四、材料图4五、研究区重磁异常分析10六、重磁资料数据处理131、重力场延拓132、磁场化极处理 163、重力场的分离 174、磁场的分离185、重磁资料导数换算处理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%以上的高品位赤铁矿。

重磁勘探课程教学实践与探索一、重磁勘探课程设置与教学内容重磁勘探课程一般作为地质类专业的选修课或专业方向课程，在地球物理学、地球科学、石油与天然气勘探等相关专业中均有涉及。

该课程主要围绕重磁勘探的基本原理、方法与技术进行教学，培养学生对地球重磁场的认识与应用能力。

课程内容主要包括：1. 重磁场基本原理：介绍地球重磁场的形成机制、特点及其在地球物理勘探中的作用；2. 重磁数据采集与处理：介绍重磁数据的采集方法、处理技术及数据解释与分析；3. 重磁勘探仪器与设备：介绍常用的重磁勘探仪器与设备的原理、结构及使用方法；4. 重磁勘探实例分析：通过实际案例分析，展示重磁勘探在地质勘探中的应用与效果。

以上内容旨在为学生建立对重磁勘探基础知识的认识和理解，为将来从事相关工作奠定基础。

1. 教学方法多样化在重磁勘探课程的教学中，应采用多种教学方法，如理论授课、案例分析、实践操作等，以提高学生的学习兴趣和参与度。

重磁勘探这门课程是一门理论联系实际的课程，因此教学内容的理论知识与实践技能应该相互结合，通过实地勘探、数据处理软件的使用等方式，激发学生的学习潜能，增强他们的实践操作能力。

2. 实验室教学与实地实习相结合重磁勘探课程的实践教学应该注重实验室教学与实地实习的结合。

通过实验室教学，学生可以了解重磁勘探仪器的使用方法、数据采集与处理技术等基础知识，而通过实地实习，学生可以亲身感受到重磁勘探在实际地质勘探中的应用效果，加深对课程内容的理解和应用能力。

1. 数据处理与解释难度较大重磁勘探课程中，数据处理与解释是一个比较困难的部分。

学生需要掌握大量的数据处理软件、方法与技术，例如：梯度法、模型法、傅里叶变换等。

在解释数据时需要对地质结构有一定的认识与理解，这对学生的地质学基础要求较高。

2. 实地勘探条件限制实地勘探是重磁勘探课程的重要内容之一，但由于实地条件的限制，学校在设置实地勘探实习课程时需要克服地域条件的限制，尽可能提供更多的实习机会与资源，以满足学生的学习需求。

重磁电勘探简介重力勘探一、重力勘探的基本概念1．重力重力的实质是牛顿万有引力和离心力的合力。

万有引力是牛顿总结前人伽里略研究行星运动规律提出来的，认为任何物体相互之间都有吸引力，吸引力的大小与两物体的质量乘积成正比，与两物体之间的距离平方成反比，其相互之间量的关系为122m m R F f R R=⋅(6—1) 式中 m 1，m 2——分别为任意两物体的质量；R——两物体相互间的距离；f——引力常数，其值在CGS制中为6．67×10—8cm3／g·s2。

上式即为牛顿万有引力定律，F力的方向对ml 来说，是由ml指向m2，对m2来说则相反。

地球是有质量的，对地球表面上任一物体来说，都有地球的吸引力。

设地球的质量为M，地面上任一物体的质量为m，则它们之间相互的吸引力F可根据式(6—1)来确定，其方向如图6—1(a)所示。

由于地球近似一个球体，对地面的m物体来说，其引力的方向指向地心。

由于地球在不断地自转，地球表面上任何物体都具有一个离心力P，其大小由下式来决定2P mr ω=(6—2)式中 r ——m 到地轴的垂向距离；ω——地球自转的角速度。

力P 的方向如图6—1(a)所示，径向指向外。

离心力P 随纬度的不同而变化，随着r 向两极减小而减小，从赤道的最大值减小到两极为零。

为了描述重力的空间分布，通常采取直角坐标系，以数学解析式表示，如图6—1(b)所示。

设地心为坐标原点，z 轴与地球的自转轴重合，x ，y 轴在赤道面上。

设任意点A 的坐标为(x ，y ，z)，地球内部某一质量单元dm 的坐标为(,,ξηζ)，A 点到dm 的距离为r ，则dm 对A 点单位质量的引力为2dm rdF f r r =(6—3)式中 r r——A 到dm 方向的单位矢量，其方向是从A 到dm 。

r 对x ，y ，z 三个坐标方向的余弦分别为：,,x y z r r rξηζ---，那么dF 在x ，y ，z 三个坐标方向的引力分量为图6—1 地球引力示意图地球的全部质量对A 点所产生的引力分量为积分号下的V 表示对整个地球进行积分。

中国地质大学（武汉）地空学院姓名：陈亮班级：061132学号：20131004480指导老师：杨宇山目录一、地质任务 (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）撰写报告。

二、工区概况研究区位于我国中东部地区，地理坐标为东经113.98697°—114.00797°，北纬29.508730°—29.522709°，处在我国非常重要的铁多金属矿成矿带西段。

在以往地质、物探工作基础上，2015年3月人们在研究区中部完成了面积为5km ²（2.5km×2km，线距50m，点距20m，测向方位角0度）的1:5000地面重磁扫面工作。

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

重磁电勘探简介欧阳歌谷（2021.02.01）重力勘探一、重力勘探的基本概念1．重力重力的实质是牛顿万有引力和离心力的合力。

万有引力是牛顿总结前人伽里略研究行星运动规律提出来的，认为任何物体相互之间都有吸引力，吸引力的大小与两物体的质量乘积成正比，与两物体之间的距离平方成反比，其相互之间量的关系为122m m R F f R R=⋅(6—1) 式中 m 1，m 2——分别为任意两物体的质量； R ——两物体相互间的距离；f ——引力常数，其值在CGS 制中为6．67×10—8cm 3／g ·s 2。

上式即为牛顿万有引力定律，F 力的方向对m l 来说，是由m l 指向m 2，对m 2来说则相反。

地球是有质量的，对地球表面上任一物体来说，都有地球的吸引力。

设地球的质量为M，地面上任一物体的质量为m，则它们之间相互的吸引力F可根据式(6—1)来确定，其方向如图6—1(a)所示。

由于地球近似一个球体，对地面的m物体来说，其引力的方向指向地心。

由于地球在不断地自转，地球表面上任何物体都具有一个离心力P，其大小由下式来决定2=(6—2)P mrω式中 r——m到地轴的垂向距离；ω——地球自转的角速度。

力P的方向如图6—1(a)所示，径向指向外。

离心力P随纬度的不同而变化，随着r向两极减小而减小，从赤道的最大值减小到两极为零。

为了描述重力的空间分布，通常采取直角坐标系，以数学解析式表示，如图6—1(b)所示。

设地心为坐标原点，z轴与地球的自转轴重合，x，y轴在赤道面上。

设任意点A的坐标为(x，y，z)，地球内部某一质量单元dm的坐标为(,,ξηζ)，A点到dm的距离为r，则dm 对A 点单位质量的引力为2dm rdF f r r =(6—3)式中 r r ——A 到dm 方向的单位矢量，其方向是从A 到dm 。

r 对x ，y ，z 三个坐标方向的余弦分别为：,,x y zr r r ξηζ---，那么dF 在x ，y ，z 三个坐标方向的引力分量为图6—1 地球引力示意图地球的全部质量对A 点所产生的引力分量为 积分号下的V 表示对整个地球进行积分。