《重磁勘探》物探1203张小涵
利用欧拉反褶积法确定水下磁性体的位置
背景场变化不大时, 可以在一定程度上减小背景 场的影响。该方法的原理如下所述。
设观测平面上 z = 0, 并假设观测值为 h, 其中 包含了背景值 b, 即 T = h - b , 将 T 代入( 5) 式, 并整理得:
x0
h x
+
y0
h y
+
z0
h z
-
Nh +
Nb
=x
h x
+
y
h y
( 6)
假设在每个矩形数据窗口内, b 是常量, 对窗
( 1. D ep ar tment of Ear th Sciences , Zhej iang Univer sity , H angz hou 310027, China; 2. CSI C N o . 715 R esearch I nstitute , H ang z hou 310012, China; 3. Geo- Ex p lo ration Science and T echnology I nstitute of J ilin Univer sity , Chang chun 130026, China)
( 1. 浙江大学 地球科学系, 杭州 310027; 2. 中国船舶重工集团公司 715 所, 杭州 310012; 3. 吉林大学 地球探测科学与技术学院, 长春 130021)
《重磁勘探》读书报告整理
《重磁勘探》——读书报告专业所在院（系、部）核工程技术学院研究生姓名郭猛猛学号 2010070807专业名称固体地球物理学日期 2011年6月30日上了重磁勘探这门课以后，我觉得弄清楚每个定义是很有必要的，所以我把这些常用到的定义总结下来，作为我的读书报告。
重力勘探:：根据地球重力场研究地球构造及寻找矿产资源的一门地球物理学科或地球物理方法。
等位面：凡位之值相等的各点所构成的曲面称为等位面。
重力位：场中任意一点的重力位等于将一单位质量从无限远处移至该点是重力所做的功。
重力等位面：重力位值相等的各点所构成的曲面。
地球的形状：两极稍平，赤道略鼓的不规则球体。
这个形状的一级近似可视为平均半径为6376km的正球面，二级近似是一个两极半径略小于赤道半径的二轴托球面。
重力的变化：包括空间因素与时间因素。
空间因素是1.地球不是一个正球体，近似于两极压缩的扁球体，地表面又是起伏不平的。
2.地球绕一定的轴旋转。
3.地下物质密度分布不均匀。
时间因素是1.太阳、月亮等天体引力引起的重力的变化，它表现有一定的周期性，也称为潮汐变化。
2.地球形状的变化和地下物质运动引起的变化为非周期性的，也称为非潮汐变化。
重力异常：在重力勘探中，将由于地下岩石、矿物密度分布不均匀所引起的重力变化，或地质体与围岩密度的差异引起的重力变化。
自由空气重力异常：重力观测值经过高度校正或自由空气校正得到的重力值称为自由空气重力异常。
布格重力异常：经过高度校正及中间层校正得到的重力值称为布格重力异常。
仪器灵敏度：单位重力的变化所能引起相对重力测量仪器的平衡体偏角的大小，如果偏角越大则表示仪器越灵敏。
基点：用重力仪在测点上进行观测是，需要有一些精度跟高、重力值一直的点来控制，这些点称为基点。
基点网：重力基点在观测是都要联成封闭的网络，这些网叫做基点网。
基点网的作用：控制重力普通电的观测精度，避免误差的积累；检查重力仪在某一段工作时间内的零点漂移，确定零点漂移校正系数；推算全区重力测点上的相对重力值或绝对重力值。
重磁力测量在石油勘探的应用及进展 ppt课件
重磁任务 磁法反演，磁性基底 的深度，垂直断层识 别，梯度分析，区域 沉积中心及沉积路径 的突出 沉积中心及沉积路径 的突出，综合盆地模 拟，密度反演，磁性 源确定，沉积磁性分 析 ，古磁分析，综合 速度分析（二维和三 维）
ppt课件
综合资料 综合钻探及露头资料， 地形资料，遥感资料， 地震资料，层序地层 分析，地震地层学分 析 综合地震资料，岩性 资料（露头及钻探）， 层序地层分析，生物 地层资料
ppt课件
19
在墨西哥湾，为了描述圈闭的翼部，在一个盐丘上方采 集了三维地震资料。处理地震资料之前首先处理和解释 了与三维地震测量同时采集的高精度重力数据。重力资 料解释预测在盐丘一侧有一个大的突出部位 （overhang），而盐丘另一侧面接近垂直。由重力资 料解释得到的这个盐丘模型，指导了地震处理人员选择 适当的叠前深度偏移速度模型，指导三维地震解释者尽 可能靠近盐丘面进行勘探。 在墨西哥湾，根据在海面上方250ft处的高分辨率航磁 测量，发现了由盐丘构造引起的少有的磁异常极小值。 同时磁法模拟确定了抗磁性盐丘的棱边及总体厚度。
ppt课件 4
一、重磁力测量所解决的油气地质问题
近年来，随着重磁力原始数据的采集速度 和精度（含各项外部改正的精度）的成倍 提高，处理解释技术的飞快进步，大比例 （1/1万～ 1/5万），高密度点的重磁力勘 探日趋广泛，过去被忽视的微弱异常研究 愈越来越受到重视，并在沉积体系划分、 油气的重力资料检测方面获得突破，正得 到广泛地运用，取得了很好的经济效益， 在油气田的开发过程中，也有人运用微重 力勘探对区块评价并配合测井等资料进行 储量计算等。
ppt课件
16
在墨西哥湾的深海区，根据卫星重力资料对局 部大地水准面所约束的海洋环流模型作出预测， 钻探者作好了充分准备，当涡流（Eddy）通过 后马上恢复钻探操作，节省了时间。
重磁勘探方法在石油勘察中的应用研究
143在自然环境当中，地质结构十分复杂，无论应用哪项地球物理方法，都会受到一定的制约和限制，例如一些重要的地质界面当中，并不一定是弹性波阻抗界面，所以在相应的地震勘探资料当中没有清晰的显示。
当前，最明显的问题为：其一，如果盖层当中具有高波阻抗界面，如存在玄武岩覆盖时，地震波产生的穿透力会有所缩减，其下各地震反射界面便不能清晰在资料当中展示。
此外，盖层下存在着并不均一的基层岩性，并且有剧烈的顶面起伏，当表面深度存在的起伏特征以及弹性波速特征，在地震勘探资料当中的显示十分混乱并且不清晰时，便出现了难以判断的问题。
为了对这一问题进行解决。
需要利用综合资料进行研究和解释，这便需要应用重、磁、地震以及相应的地球物理资料给予解释。
目前，对于油气田当中的各种物化探资料实施了综合考量，使得资料地质解释的可靠性有了非常大的提升。
对于重、磁方法的应用，在石油勘探当中有广泛应用，有益于提升对资料进行解释的水平。
1 重、磁勘探技术的进展分析1.1 提升重、磁勘探仪器的发展和精度当前，对于重、磁仪器制造的工艺水准，有了非常大的改善和提升，仪器当中的技术指标，借助科学力量的发展，达到了全新的水平，尤其是电子技术对其产生的带动作用，已完成了自动化以及电子化的发展，并向智能化逐步迈进，精度有了非常大的提升。
因为仪器的探测精度越来越强，如没有相关干扰，会对非常多的异常细节以及微弱的异常特征信息进行有效获取[1]。
1.2 广泛应用计算机以及计算机技术计算机科技力量的发展，为重、磁勘探提供了更加有力的技术支撑以及保障，是在石油勘探过程中广泛应用重、磁方法的基础前提[2]。
在80年代之后，计算机的发展速度非常迅猛，如速度更快、容量更大等等，使得重、磁勘测方法在具体应用时开拓了更多的功能，发展前景广阔。
目前，重、磁资料当中的人机交互解释方法，以及图形图像快速处理技术越来越完善，相关技术内容也更加成熟，所以在使用时十分便捷和高效。
1.3 更新和完善解释方法重、磁资料解释属于重、磁探测环节当中的关键性内容，尤其是在近40年的研究当中，已经从之前处理和解释单一空间领域逐步发展到空间域以及波数域两大处理系统；从以往形态简单、物性均匀的物体正反演逐步发展到了任意形状、不均匀物质的正反演，并且之前的解释方法非常单一，目前已经发展到了重、磁联合反演，重震联合反演，并可以综合处理和解释物探资料[3]。
重磁电法在金属矿山勘探中的应用及优化策略研究
181地质勘探G eological prospecting重磁电法在金属矿山勘探中的应用及优化策略研究谢宇飞甘肃省地质调查院，甘肃 兰州 ７３００００摘 要：随着全球矿产资源需求的日益增长，金属矿山勘探成为地球科学领域的关键任务，直接影响着资源开发的可持续性和效率。
重磁电法融合重力法、磁力法和电法，能够提供多维度的地下结构信息，尤其在复杂地质环境中表现出其优越性。
对此，本文首先对重磁电法的技术原理进行介绍，然后对重磁电法的应用流程进行分析，并结合实例对金属矿山勘探中重磁电法的应用要点以及优化策略展开探究，旨在提高勘探效率和准确性。
关键词：金属矿山；勘探；重磁电法中图分类号：Ｐ６１８．２ 文献标识码：Ａ 文章编号：１００２－５０６５（２０２４）０１－０１８１－３Research on the Application and Optimization Strategy of Gravity Magnetoelectric Methodin Metal Mine ExplorationXIE Yu-feiＧｅｏｌｏｇｉｃａｌ Ｓｕｒｖｅｙ ｏｆ Ｇａｎｓｕ Ｐｒｏｖｉｎｃｅ，Ｌａｎｚｈｏｕ ７３００００，ＣｈｉｎａAbstract: Ｗｉｔｈ ｔｈｅ ｉｎｃｒｅａｓｉｎｇ ｄｅｍａｎｄ ｆｏｒ ｇｌｏｂａｌ ｍｉｎｅｒａｌ ｒｅｓｏｕｒｃｅｓ， ｍｅｔａｌ ｍｉｎｉｎｇ ｅｘｐｌｏｒａｔｉｏｎ ｈａｓ ｂｅｃｏｍｅ ａ ｋｅｙ ｔａｓｋ ｉｎ ｔｈｅ ｆｉｅｌｄ ｏｆ ｅａｒｔｈ ｓｃｉｅｎｃｅ， ｄｉｒｅｃｔｌｙ ａｆｆｅｃｔｉｎｇ ｔｈｅ ｓｕｓｔａｉｎａｂｉｌｉｔｙ ａｎｄ ｅｆｆｉｃｉｅｎｃｙ ｏｆ ｒｅｓｏｕｒｃｅ ｄｅｖｅｌｏｐｍｅｎｔ． Ｔｈｅ ｉｎｔｅｇｒａｔｉｏｎ ｏｆ ｇｒａｖｉｔｙ， ｍａｇｎｅｔｉｃ， ａｎｄ ｅｌｅｃｔｒｉｃａｌ ｍｅｔｈｏｄｓ ｃａｎ ｐｒｏｖｉｄｅ ｍｕｌｔｉ－ｄｉｍｅｎｓｉｏｎａｌ ｕｎｄｅｒｇｒｏｕｎｄ ｓｔｒｕｃｔｕｒａｌ ｉｎｆｏｒｍａｔｉｏｎ， ｅｓｐｅｃｉａｌｌｙ ｉｎ ｃｏｍｐｌｅｘ ｇｅｏｌｏｇｉｃａｌ ｅｎｖｉｒｏｎｍｅｎｔｓ， ｄｅｍｏｎｓｔｒａｔｉｎｇ ｉｔｓ ｓｕｐｅｒｉｏｒｉｔｙ． Ｉｎ ｔｈｉｓ ｒｅｇａｒｄ， ｔｈｉｓ ａｒｔｉｃｌｅ ｆｉｒｓｔ ｉｎｔｒｏｄｕｃｅｓ ｔｈｅ ｔｅｃｈｎｉｃａｌ ｐｒｉｎｃｉｐｌｅｓ ｏｆ ｔｈｅ ｇｒａｖｉｔｙ， ｍａｇｎｅｔｉｃ ａｎｄ ｅｌｅｃｔｒｉｃａｌ ｍｅｔｈｏｄ， ａｎａｌｙｚｅｓ ｔｈｅ ａｐｐｌｉｃａｔｉｏｎ ｐｒｏｃｅｓｓ ｏｆ ｔｈｅ ｇｒａｖｉｔｙ， ｍａｇｎｅｔｉｃ ａｎｄ ｅｌｅｃｔｒｉｃａｌ ｍｅｔｈｏｄ， ａｎｄ ｅｘｐｌｏｒｅｓ ｔｈｅ ａｐｐｌｉｃａｔｉｏｎ ｐｏｉｎｔｓ ａｎｄ ｏｐｔｉｍｉｚａｔｉｏｎ ｓｔｒａｔｅｇｉｅｓ ｏｆ ｔｈｅ ｇｒａｖｉｔｙ， ｍａｇｎｅｔｉｃ ａｎｄ ｅｌｅｃｔｒｉｃａｌ ｍｅｔｈｏｄ ｉｎ ｍｅｔａｌ ｍｉｎｅ ｅｘｐｌｏｒａｔｉｏｎ ｔｈｒｏｕｇｈ ｅｘａｍｐｌｅｓ， ａｉｍｉｎｇ ｔｏ ｉｍｐｒｏｖｅ ｅｘｐｌｏｒａｔｉｏｎ ｅｆｆｉｃｉｅｎｃｙ ａｎｄ ａｃｃｕｒａｃｙ．Keywords: ｍｅｔａｌ ｍｉｎｅｓ； Ｅｘｐｌｏｒａｔｉｏｎ； Ｇｒａｖｉｔｙ ｍａｇｎｅｔｏ ｅｌｅｃｔｒｉｃ ｍｅｔｈｏｄ收稿日期：２０２３－１１作者简介：谢宇飞，男，生于１９９２年，汉族，甘肃天水人，本科，工程师，研究方向：地球物理重磁电固体矿产勘查。
1重磁勘探-绪论
地 质
物 理
第27页，共46页。
5、地球物理场（数据）
地球物理场：是指存在于地球内部及其周围的、具有物理作
用的物质空间。
例如：1）地球内部及其周围具有重力作用的物质空间，称为重力 场；
2）地球的磁场以及地壳内磁性体产生的磁场。
地球的重力场﹑地磁场都是天然存在的地球物理场。
• 地球物理场：各种地球物理方法在地表或地表附近测量的各 种物理现象的信息可以统称为地球物理场的信息。
•物探方法的基础是建立在什么上面的？ 第8页，共46页。
2、应用地球物理学学科分类
△ 按工作原理分类 重力勘探、磁法勘探、电法勘探、
地震勘探、放射性勘探、地热勘探。
Radioactive Heat
第9页，共46页。
基本实质：利用岩矿石的六种主要物理性质或物性参数，建 立相应的六种应用地球物理方法。
重磁勘探-绪论
第1页，共46页。
教材与主要参考书
• 1、《煤田重力与磁法勘探》，高培泽等编. 煤炭工业出版社 ，1990。 • 2、应用地球物理教程（重力勘探、磁法勘探 ），罗孝宽等编.北京：
地质出版社，1991 。 • 3、重力场与重力勘探，曾华霖.北京：地质出版社， 2005 。 • 4、地磁场与磁力勘探，管志宁.北京：地质出版社， 2005 。 • 5、应用地球物理学教程，傅良魁等编.北京：地质出版社，1991 。 • 6、勘探重力学与地磁学，孟令顺编.北京：地质出版社，2008
◎ 地质工作各阶段 ◎ 全球与区域构造研究研究
◎ 工程无损检测 ◎ 考古调查
第37页，共46页。
7、地球物理勘探面临的任务、问题和发
展趋势
地球物理勘探面临的任务
1.基础地质调查
《重磁勘探》物探1203张小涵复习进程
《重磁勘探》物探1203张小涵《重磁勘探》课程设计报告专业：勘查技术与工程班级：物探 1203班姓名：张小涵学号： 201211020321 指导教师：张春灌二〇一四年十二月二日《重磁勘探》课程设计任务书目录第1章引言 (1)第2章模型设计 (2)第3章程序实现 (4)第4章结果分析 (6)第5章总结 (10)参考文献 (11)第一章引言一、课程设计目的和意义本次课程设计是《重磁勘探》学习的延续，目的是巩固所学的重力和磁法勘探的理论知识，加深对基本原理的理解，会用所学程序设计语言（如Matlab，Surfer）等软件完成课程设计题目的程序编写、数据计算，利用现有绘图软件完成数据成图。
二、课程设计内容两个不同空间位置的球形矿体在地面引起的重力异常。
三、课程设计要求1.设计两个不同空间位置的球形矿体，并画出其平面和主剖面空间位置示意图(假设地面水平)。
2.推导两个不同空间位置的球形矿体在地面引起重力异常的计算公式，利用推导出的公式，用Matlab语言编程实现计算重力异常。
3.根据所提供的课程设计报告格式编写报告，报告中应附上课程设计任务书，报告内容应包括所推导的公式、所编写的程序、结果图件等。
第二章 模型设计一、 计算公式23222323222)(34)(h Y X h R G h Y X GMh g ++=++=∆σπ球若常数G 取6.67*10-3，剩余密度σ的单位取g/cm 3，半径R 的单位取m ，中心埋深h 的单位取m ，则重力异常Δg 的单位为mGal 。
二、模型参数三、模型图件图2-1 球体平面空间位置示意图图2-2 球体主剖面空间位置示意图第三章程序实现一、 Matlab程序语句a=500;b=500;clf;GA=6.67e-003;DPI=3.1415926;P1=0.5;P2=1.0;R1=50.0;R2=75.0;H1=100.0;H2=175.0;M1=(4.0/3.0)*DPI*(R1^3)*P1;M2=(4.0/3.0)*DPI*(R2^3)*P2;x=0:5:a;y=0:5:b;for i=1:length(y)for j=1:length(x)z(i,j)=(GA*M1*H1)/(((x(j)-100)^2+(y(i)-100)^2+H1^2)^1.5)+(GA*M2*H2)/(((x(j)-325)^2+(y(i)-275)^2+H2^2)^1.5);endendaxis([0,a,0,b,min(min(z)),max(max(z))]);colormap(flipud(winter));surf(x,y,z),xlabel('X 轴 (m)'),ylabel('Y 轴 (m)'),title('两个不同空间位置的球形矿体在地面引起的重力异常');save chenwentao.dat z -ascii;save chenwentao z;二、 Matlab软件所显示的结果图件图3-1 利用MATLAB软件实现球形矿体重力异常图图3-2利用MATLAB显示X-Z方向球形矿体重力异常图第四章结果分析一、制图流程1.1 先利用给定的模型参数在Surfer软件上设计出给定地下地质球体，分别绘制出给定球形矿体的平面及主剖面空间位置图；1.2 根据给定模型参数及Sufer软件绘制图件写出Matlab程序语言，分别求出不同条件下的球形矿体重力异常数据；1.3 在Matlab软件上输入程序语言，同时得出球形地质体在地面产生异常的三维立体图；1.4 利用异常数据在Surfer软件上显示不同条件下的球形矿体重力异常的平面等值线图、矢量图、影像图等一系列图件，并对其结果进行理论分析。
高精度重磁技术在矿山采区探测中的应用
地质勘探G eological prospecting高精度重磁技术在矿山采区探测中的应用薛玲香（中石化石油工程地球物理有限公司地理地质信息勘查分公司，湖北 武汉 ４３０２２３）摘 要：在矿山地质勘察中， 采区探测是一个较新的研究课题。
在常规的勘探和电磁法研究基础上，针对矿山采区的成矿特征和分布规律， 提出高精度重磁技术在矿山采区探测中的应用。
在对磁异常进行计算的基础上、进行磁异常弱信息提取与增强，之后根据重磁异常推断矿区边界，得出有效结果。
试验结果表明，本文方法在矿山采区探测中，与传统方法相比，其探测精度有较大提高，具有较强的有效性和实用性。
关键词：高精度重磁技术；矿山；采区；探测中图分类号：Ｐ６３１ 文献标识码：Ａ 文章编号：１００２－５０６５（２０２０）０４－０１３５－２Application of High Precision Heavy Magnetic Technology in Mining Area ExplorationXUE Ling-xiang（Ｓｉｎｏｐｅｃ Ｇｅｏｐｈｙｓｉｃａｌ Ｃｏｒｐｏｒａｔｉｏｎ，Ｇｅｏｌｏｇｉｃａｌ Ｇｅｏｌｏｇｉｃａｌ Ｉｎｆｏｒｍａｔｉｏｎ Ｅｘｐｌｏｒａｔｉｏｎ Ｂｒａｎｃｈ，Ｗｕｈａｎ ４３０２２３，Ｃｈｉｎａ）Abstract: Ｉｎ ｍｉｎｅ ｇｅｏｌｏｇｉｃａｌ ｓｕｒｖｅｙ， ｍｉｎｉｎｇ ａｒｅａ ｅｘｐｌｏｒａｔｉｏｎ ｉｓ ａ ｒｅｌａｔｉｖｅｌｙ ｎｅｗ ｒｅｓｅａｒｃｈ ｔｏｐｉｃ． Ｏｎ ｔｈｅ ｂａｓｉｓ ｏｆ ｃｏｎｖｅｎｔｉｏｎａｌ ｅｘｐｌｏｒａｔｉｏｎ ａｎｄ ｅｌｅｃｔｒｏｍａｇｎｅｔｉｃ ｍｅｔｈｏｄ， ｔｈｅ ａｐｐｌｉｃａｔｉｏｎ ｏｆ ｈｉｇｈ ｐｒｅｃｉｓｉｏｎ ｒｅｍａｇｎｅｔｉｓｍ ｔｅｃｈｎｏｌｏｇｙ ｉｎ ｍｉｎｉｎｇ ｅｘｐｌｏｒａｔｉｏｎ ｉｓ ｐｕｔ ｆｏｒｗａｒｄ ｉｎ ｖｉｅｗ ｏｆ ｔｈｅ ｍｅｔａｌｌｏｇｅｎｉｃ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ａｎｄ ｄｉｓｔｒｉｂｕｔｉｏｎ ｌａｗ ｏｆ ｍｉｎｉｎｇ ａｒｅａ． Ｂａｓｅｄ ｏｎ ｔｈｅ ｃａｌｃｕｌａｔｉｏｎ ｏｆ ｍａｇｎｅｔｉｃ ａｎｏｍａｌｙ， ｔｈｅ ｗｅａｋ ｉｎｆｏｒｍａｔｉｏｎ ｏｆ ｍａｇｎｅｔｉｃ ａｎｏｍａｌｙ ｉｓ ｅｘｔｒａｃｔｅｄ ａｎｄ ｅｎｈａｎｃｅｄ， ａｎｄ ｔｈｅｎ ｔｈｅ ｂｏｕｎｄａｒｙ ｏｆ ｍｉｎｉｎｇ ａｒｅａ ｉｓ ｉｎｆｅｒｒｅｄ ａｃｃｏｒｄｉｎｇ ｔｏ ｈｅａｖｙ ｍａｇｎｅｔｉｃ ａｎｏｍａｌｙ， ａｎｄ ｔｈｅ ｅｆｆｅｃｔｉｖｅ ｒｅｓｕｌｔ ｉｓ ｏｂｔａｉｎｅｄ． Ｔｈｅ ｔｅｓｔ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ ｄｅｔｅｃｔｉｏｎ ａｃｃｕｒａｃｙ ｏｆ ｔｈｉｓ ｍｅｔｈｏｄ ｉｓ ｇｒｅａｔｌｙ ｉｍｐｒｏｖｅｄ ｃｏｍｐａｒｅｄ ｗｉｔｈ ｔｈｅ ｔｒａｄｉｔｉｏｎａｌ ｍｅｔｈｏｄ， ａｎｄ ｉｔ ｈａｓ ｓｔｒｏｎｇ ｖａｌｉｄｉｔｙ ａｎｄ ｐｒａｃｔｉｃａｂｉｌｉｔｙ．Keywords: ｈｉｇｈ ｐｒｅｃｉｓｉｏｎ ｒｅｍａｇｎｅｔｉｓｍ ｔｅｃｈｎｏｌｏｇｙ； ｍｉｎｅ； ｍｉｎｉｎｇ ａｒｅａ； ｄｅｔｅｃｔｉｏｎ通过重磁异常对目标体边界进行探测是地质物探技术解释的一项极其重要的研究工作。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《重磁勘探》课程设计报告
专业：勘查技术与工程班级：物探1203班姓名：张小涵
学号：201211020321 指导教师：张春灌
二〇一四年十二月二日
《重磁勘探》课程设计任务书
目录
第1章引言 (1)
第2章模型设计 (2)
第3章程序实现 (4)
第4章结果分析 (6)
第5章总结 (10)
参考文献 (11)
第一章引言
一、课程设计目的和意义
本次课程设计是《重磁勘探》学习的延续，目的是巩固所学的重力和磁法勘探的理论知识，加深对基本原理的理解，会用所学程序设计语言（如Matlab，Surfer）等软件完成课程设计题目的程序编写、数据计算，利用现有绘图软件完成数据成图。
二、课程设计内容
两个不同空间位置的球形矿体在地面引起的重力异常。
三、课程设计要求
1.设计两个不同空间位置的球形矿体，并画出其平面和主剖面空间位置示意图(假设地面水平)。
2.推导两个不同空间位置的球形矿体在地面引起重力异常的计算公式，利用推导出的公式，用Matlab语言编程实现计算重力异常。
3.根据所提供的课程设计报告格式编写报告，报告中应附上课程设计任务书，报告内容应包括所推导的公式、所编写的程序、结果图件等。
第二章 模型设计
一、 计算公式
2
322
2
32
322
2
)
(34
)
(h Y X h R G h Y X GMh g ++=
++=
∆σπ球
若常数G 取6.67*10-3，剩余密度σ的单位取g/cm 3，半径R 的单位取m ，中心埋深h 的单
位取m ，则重力异常Δg 的单位为mGal 。
二、模型参数
三、模型图件
图2-1 球体平面空间位置示意图
图2-2 球体主剖面空间位置示意图
第三章程序实现
一、 Matlab程序语句
a=500;
b=500;
clf;
GA=6.67e-003;
DPI=3.1415926;
P1=0.5;
P2=1.0;
R1=50.0;
R2=75.0;
H1=100.0;
H2=175.0;
M1=(4.0/3.0)*DPI*(R1^3)*P1;
M2=(4.0/3.0)*DPI*(R2^3)*P2;
x=0:5:a;
y=0:5:b;
for i=1:length(y)
for j=1:length(x)
z(i,j)=(GA*M1*H1)/(((x(j)-100)^2+(y(i)-100)^2+H1^2)^1.5)+(GA*M2*H2)/(((x(j)-3 25)^2+(y(i)-275)^2+H2^2)^1.5);
end
end
axis([0,a,0,b,min(min(z)),max(max(z))]);
colormap(flipud(winter));
surf(x,y,z),xlabel('X 轴 (m)'),ylabel('Y 轴 (m)'),title('两个不同空间位置的球形矿体在地面引起的重力异常');
save chenwentao.dat z -ascii;
save chenwentao z;
二、 Matlab软件所显示的结果图件
图3-1 利用MATLAB软件实现球形矿体重力异常图
图3-2利用MATLAB显示X-Z方向球形矿体重力异常图
第四章结果分析
一、制图流程
1.1 先利用给定的模型参数在Surfer软件上设计出给定地下地质球体，分别绘制出给定球形矿体的平面及主剖面空间位置图；
1.2 根据给定模型参数及Sufer软件绘制图件写出Matlab程序语言，分别求出不同条件下的球形矿体重力异常数据；
1.3 在Matlab软件上输入程序语言，同时得出球形地质体在地面产生异常的三维立体图；
1.4 利用异常数据在Surfer软件上显示不同条件下的球形矿体重力异常的平面等值线图、矢量图、影像图等一系列图件，并对其结果进行理论分析。
二、Sufer软件部分图件显示
图4-1 利用Sufer软件显示不同球形矿体重力异常表面图
图4-2 利用Sufer软件显示不同球形矿体重力异常矢量图
图4-3 利用Sufer软件显示不同球形矿体重力异常渐变地形图
图4-4 利用Sufer软件显示不同球形矿体重力异常线框图
图4-5 利用Sufer软件显示不同球形矿体重力异常影像图
图4-6 利用Sufer软件显示不同球形矿体重力异常平面等值线图
三、结果分析
通过观察分析上述不同球形矿体重力异常图件，我们可以得出以下结论：
（1）不同的地下球形地质体在地面产生的重力异常异常明显不同；
(2) 半径越大的地下球形地质体产生的重力异常在地面的表现越突出；
(3) 每个球形地质体产生的异常曲线是围绕其异常中心半径不等的同心圆；
(4) 地下球体的存在所产生的的重力异常矢量方向都是沿其异常中心向外呈发散状；
(5) 剩余密度的不同对于地下球形矿体所产生的重力异常的影响也是很大的，所以在我们实际的勘探工作中对所遇到的重力异常要进行综合分析，注意多因素的存在对数据的影
响，同时还应结合工区内的地质资料进行重力异常数据的反演，减少多解性的影响和干扰。
第五章总结
当老师跟我们通知这个《重磁勘探》这门课程之后的课程设计可能需要一周的时候，我觉得应该不会那么复杂的吧，所以就简单的看了老师给的课程设计指导书，以为很短的时间就可以完成这个任务，但是当自己真正实践操作的时候就暴露出了自己在实践方面的不足和使用常用专业软件技能方面的欠缺，因为之前的课程实验中基本接触并简单使用过Sufer软件，所以实际操作起来难度并不是很大，但还是发现自己对于具体操作中的一些不足，例如对于球体平面及主剖面位置绘制的时候并不知道如何确定数据的范围，如何使得表格的主分割线不影响我们图件的正常显示，这些都在老师的悉心指导下得以完成。
还有对于Matlab软件的使用更显生疏，因为是第一次接触，之前也未在课程设计前做好基本的预习工作，在编写Matlab程序语言的时候根本无从下手，所以只能借鉴已有的程序去分析每一条程序的具体含义，之后才得以顺利写出程序。
我们利用Matlab软件将所需数据计算出来后还涉及到利用该数据在Sufer软件下做出各种具体的球体重力异常图件，在这一方面涉及的数据的改写和文件格式的转化我基本上完全不会，最后还是在老师的指导下顺利完成本次课程设计。
从早上9:00开始一直到晚上6:00，我整个一天的时间基本都在实验室完成本次课程设计，自己原本以为就几个小时可以完成的任务用了9 个小时左右，一方面是由于自己前期知识的欠缺，并未对本次课程设计所需的软件有所了解，另一方面也在这具体的实验过程中更加强化对这一知识的理解，随着学习的不断深入，我觉得对常用软件的熟练操作是我们必备的技能之一，通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，从而提高自己的实际动手能力和独立思考的能力。
在设计的过程中遇到问题，通过努力最终都得到了解决。
在这里，由衷的感谢老师和同学们的帮助。
此次设计也让我明白了思路即出路，有什么不懂不明白的地方要及时请教或上网查询，只要认真钻研，动脑思考，动手实践，就没有弄不懂的知识。
参考文献
[1] 袁炳强，张春灌编著，重磁勘探，西安：西安石油大学, 2014年
[2] 焦新华，重力与磁法勘探，北京：地质出版社，2009年。