重磁剖面2.5D反演拟合使用说明
浅谈重磁电剖面研究工作中图件处理心得
浅谈重磁电剖面研究工作中图件处理心得摘要：随着地质调查工作的不断开展，重磁方法在找矿工作中的应用得以发展，其作为指导地质找矿的方法之一，数据采集后如何快速成图，就成为内业整理人员的头等大事，本文为此展开分析，对平时工作中在重磁电剖面图件处理的一些心得进行分项。
希望为地质项目提供了有力的基础地质数据支撑，发挥了巨大的作用。
关键词：重磁，图件处理，应用成果，A brief discussion on the drawing processing experience in the research of heavy magnetic and electric profileLiShao-hua WangPengShaanxi Second Comprehensive Geophysical Exploration GroupCo.,LTD.,Xi'an710016,ChinaAbstract:With the continuous development of geological survey work,the application of gravity and magnetic method in prospecting work has been developed. As one of the methods to guide geological prospecting, how to quickly map after data acquisition has become the top priority of the internal industry collation personnel. This paper analyzes this, and some of the experience of processing the gravity, magnetic and electric profile pieces in the normal work are pided into several items. Hope provides strong basic geological data support for geological projects and plays a huge role.Key words: gravity and magnetism, map processing, application results,1资料处理方法技术重力异常是叠加场，是地下不同密度体的总体反映。
重磁电联合反演方法在天津基岩构造研究中的应用_马为
第37卷第3期地质调查与研究Vol.37No.3 2014年9月GEOLOGICAL SURVEY AND RESEARCH Sep.2014重磁电联合反演方法在天津基岩构造研究中的应用马为，李世斌，徐新学，郑军，赵洪鹏，李华强（天津市地球物理勘探中心，天津300170）摘要:以地层和岩石的综合物性研究为基础，采用重磁电联合反演方法，对横穿天津地区的重力、大地电磁(MT)长剖面进行数据处理，并对天津基岩构造中主要地层、断裂、隐伏岩体特征等进行了综合解释。
结果表明，重磁电联合反演对克服多解性、提高定量解释精度具有较好的效果。
本文还通过应用实例，介绍了重磁电综合分析解释方法的运用，可供类似地区的基岩地质构造研究借鉴和参考。
关键词:重磁电联合反演；构造；地层；定量解释中图分类号：P631.1;P631.2;P631.3文献标识码：A文章编号：1672－4135（2014）03－0212-05各类地球物理方法是根据地质构造及岩石不同的物理参数进行解释推断，而各单一地球物理反演结果各有侧重，不可避免的存在反问题的多解性，影响了地质解释成果的精度。
联合反演方法以同一地质目标体的多种地球物理参数（电性、密度、磁性）为基础，互为约束,可以提高反演结果的分辨率与置信度，适应了地球物理综合解释的需要，因此在基岩构造研究、油气矿产勘查等多种领域内得到了越来越多的重视，其中应用较多的联合反演为重磁电方法的联合[1-7]。
为研究天津地区的基岩地质构造，笔者完成了北东方向的重力、大地电磁（MT）综合长剖面。
该剖面西起武清，东至北塘，与本区主要构造单元走向垂直，为进行重磁电联合反演提供了基础。
鉴于本区具有丰富的地层和岩石物性研究成果和钻孔资料，可以为联合反演提供准确的约束条件，同时联合反演技术可以结合重力划分断裂和大地电磁进行垂向地层分层的优势，提高基岩构造解释工作的精度。
故本文采用重磁电联合反演技术，将地质与物探有机结合，建立最优化的地质－地球物理模型，解决了天津基岩构造研究中地层、断裂、隐伏岩体等问题，划分了“两坳一隆”的格局，取得了较好的效果。
重磁资料数据处理某些方法.
常的影响被削弱了年, ‘ 我们用上述方法编制了 , 币机算法语言程序 , 并作了模型试验 , , 后来又处理了实际资料闭有一定效果。
原始数据中存在有局部异常成分也是造成虚假异常的原因之一数据点随机分布叠加异常的划分趋势分析方法去处理规则格网数据吼 , 所以可以用不规则测网。
‘’ 至今仍是一个没有很好解决的问题、 , 还有许多工作要做二异常的反演 , 异常的反演一多模型最优化选择法、即根据观测异常求地质体的位置。
大小、产状和物性参数数据处理和解释工作中的一个不可缺少的环节欲反演的目标大体上有三种—是重磁资料求地质体用规则几何形体近似 , 的几何参数和物求一定范围物。
性参数性的空间分布的实质在于 , 。
求一定范围的物性分界面 , 反映某一地质层位的起伏目前采用的反演方法可分为两类直接法 , 。
根据观测整理的异常直接曲线拟合一选择法选择法 , 计算地质体的某些参数多用于解释简单的异常一次完成。
将实测异常曲线与一系列已知形状模型体产生的理论异常曲线进行比较当实测曲线与某一理论曲线符合为实际地质体的近似结果。
在给定的误差范围内选择法需迭代完成。
时 , 就将该理论曲线所对应的模型体作这里先介绍多模型最优化选择法在一个矿区、的应用情况 , 然后讨论此法的应用条件 , 年 , 武汉地质学院磁法组应用长方体组合模型采用了十五个模型 , 改进的马奎特法。
〕冀东对。
‘ 一区的地磁异常进行了反演川得到了各模型体的参数的理论曲线。
图表示反演得到的模型体的平面位置和由它算出理论曲线与实测曲线有些模型体的参数也。
一拟合得较好 , , 滋儡节静之火浓姆丫一之文侧气二么爹又一夕‘ 之毛‘ ‘’ 币—、、比较接近附近钻孔中的见矿情况月又卜女例如间见到体 , 孔位于第块之间 , 块和第米剩图一汤火在井深入米赤铁矿和 , 米磁性矿米磁块和第按计算结果是该处应有又如在第一、 , 性矿体块之间的孔 , 按反演推米磁性矿米到米米磁铁算在体 , 米左右应有实际钻探结果在当然此外 , 之间见到了矿区。
人工鱼算法的重磁位场反演方法
优 状态 就替 代 之 。 该 算 法 具 有 良好 的求 取 全 局 极 值 能 力 ， 具有 并 对初值、 参数选择不敏感、 鲁棒性强、 简单易实现等 诸 多优 点 ， 但是 当 问题 规 模 过 大 时 求 解 困难 并 且求 解 初期 收敛 较快 ， 期较 慢 。 后 人工 鱼 群算 法 的改 进研 究， 如在 算法 中增 加 鱼群 的协 调 行 为 ， 人 工鱼 群算 将 法与 模拟 退 火算 法相 融 和 的混合 智 能算 法 等 。人工 鱼 群算 法 目前 已用 于 组 合优 化 、 数估 计 、 Ｉ 控 参 ＰＤ 制器 的参 数 整定 、 神经 网络 优化 等 。 ２ 算 例分 析 下 面 理 论 模 型 数 值 计算 有 两 个 目的 ： 一 是 为 其 了检 验 ＡＦ Ａ 用于 重 磁 位 场 反演 的有 效 性 ； 二 是 Ｓ 其 为 了说 明总磁 异常 对 不 同尺度 磁 源反 演 的效 里 。我 们 设计 如 图 １所 示 的 直 立板 板 组 合 模 型 ， Ｌ个 直 用 立 板拟 合 一 个 二度 异 常 体 ， 各个 板 上 顶 和 下 底 埋 深 分 别 为 ｈ ｌ ｈ ２ ｈ １ ｈ ２ … 、 ｉ、ｉ、 、ｌ、１， ｌ 、 １ 、 ２ 、 ２ 、 ｈ１ｈ２ … ｈｌｈ２ 第 一个 板 中心 的横 坐 标为 Ｘ ， 的 宽度 均 为 ｂ 再 加 Ｏ板 ， 上 磁化 强度 的分 量 Ｊ 、ｚ其 中 ： ｘ Ｊ， Ｌ一５ ｈ 一 ［ ． ． ， ｌ ４ ０６ ０
Ｓｅ ２ 人 工 鱼 从 一 个 状 态 转 移 到 另 一 个 状 态 ： ｔｐ ：
Ｘｉ Ｘｉ Ｒｓ ｄ ｍ （ ｉ ｕ １ ； ＝ ＋ ｎ ｏ ｖｓ ａ ）
Ｓｅ ３ 如 果 适 应 值 ｙ ＜ Ｙｊ则 Ｘｉ 一 Ｘｉ ｔｐ ： ｉ ， Ｉ ＋
重磁资料正反演
成都理工大学工程技术学院勘查技术与工程专业重磁资料处理报告书班级：学号：姓名：指导老师：适用专业：勘察技术与工程提交时间：年月日目录(请使用更新域)目录(请使用更新域) ................................................................................................ - 1 - 1 重磁勘探原理简述. (2)1.1重力勘探原理简述（标题使用黑体三号） (2)1.2 磁法勘探原理简述 (2)2 空间域常规重磁资料处理方法 (3)2.1 原始模型部分 (3)2.2 重磁异常的导数换算 (3)2.2.1 重力一阶导数 (3)2.2.2 重力垂向二阶导数 (5)2.2.3 磁法一阶导数 (5)2.2.4 磁法二阶导数 (6)2.3 重磁异常的空间换算 (8)2.3.1 向上延拓 (8)2.3.2 向下延拓 (9)2.4磁异常分量及不同磁化方向二度磁异常之间的换算 (10)2.5 综合应用 (11)3 频率域中重磁异常的换算...................................................... 错误！未定义书签。
3.1......................................................................................... 错误！未定义书签。
3.2......................................................................................... 错误！未定义书签。
3.3......................................................................................... 错误！未定义书签。
Res2dinv反演软件完成反演后
Res2dinv反演软件完成反演后，可以将反演数据保存为surfer数据格式，以便在surfer软件中重新绘制。
输出时会自动生成4个数据文件，一般生成模型数据格式为：原文件名_modres.dat原文件名_modres.bln原文件名_modres.lvl原文件名_modres_post.datmod是模型的英文简写，res是电阻率的英文简写。
如下图如果是生成带地形的模型数据格式为原文件名_topres.dat原文件名_topres.bln原文件名_topres.lvl原文件名_topres_post.dattop是地形的英文简写。
第一个文件是反演模型电阻率值的surfer格式数据文件。
第二个文件是模型值的surfer边界文件，是对网格文件做白化是用的。
第三个文件是surfer等级数据文件，用颜色来表示模型电阻率的大小。
第四个文件是张贴图数据文件，用于在地表面上表示斜距。
下面我们一步步使用这4个文件来把surfer等值线画出来。
假设原文件为TEST.DAT。
一，首先使用test_topres.dat文件在surfer8中点击“网格”---“数据”—读入“test_topres.dat”文件。
网格化数据时，网格化方法可选择克里格方法或自然邻近点法，克里格法在网格化较小数值时，容易插值出负值，而自然邻近点法则不会出现此问题。
网格线的间距要选小一点，即行数要多一些，这样做白化时，白化的边界不会出现很严重的锯齿状条纹。
当然，如果行数太大，计算量也较大，实际中根据情况灵活选择，满足要求即可。
网格文件创建后，网格化插值算法会将数据文件插值为一个矩形的网格文件(以X，Y方向的最大值插满整个网格)如果直接绘制等值线图，会形成一个矩形的等值线图，这并不是我们需要的，因为反演后的模型数据不一定是矩形的，而且有地形的话，形状还会更复杂一些，这时，就需要使用网格白化的功能将没有实际数据的地方但又插值出的数据去掉。
一，使用test_topres.bln边界文件点击“网格”---“白化”，从弹出的打开网格对话框里选中上一步生成的网格文件“test_topres.grd”，打开后，从打开边界文件对话框中选择“test_topres.bln”边界文件，打开后会出现数据白化后要保存的网格文件名对话框，一般把原来的网格文件覆盖即可。
RGIS教程软件结构和功能
工作空间管理
RGIS2009基本操作
• 系统提供了不同的数据查看方式：
– 地图窗口 – 浏览窗口 – 布局窗口
地图窗口
• 以常规地图的方式表达信息，使数据的 地理模式直观化。 是最常用的地图显示 和操作窗口，包含了常用的编辑工具。
浏览窗口
• 以列表方式表达信息，通过浏览窗口可以全面 查看表中各项数据值和选定数据。该窗口数据 可以实现与地图窗口的联动。
等。
地图图层
• 系统是按图层的方式组织表，一张表构 成一个图层实体。每一个图层包含整个 地图的不同方面。每张表都作为独立的 图层显示。每个图层可以包含不同的地 图对象（点、线、面、文本）。
地图图层
• 例如，第一个图层包含重力异常等值线， 第二个图层是表示异常线的标注，第三 个图层是异常图的面元，把它们层层叠 加就形成一幅完整的地图。在创建每一 个图层时，都要为其建立一张表，系统 就是通过这种方式在表与地图之间建立 了联系。（见下图）
RGISMAP主要功能
• 数据库建设和维护：重力、磁测数据库建设和维护的工具，能够 实现对省级和全国重力和磁测数据库管理和维护。
• 数据预处理：坐标转换、网格化、空白区处理、数据扩边、地磁 要素计算、磁化强度计算、剖面数据准备。
• 重磁数据处理功能：空间域/ 频率域的延拓、导数、曲化平，异常 分离。 • 磁测数据处理：化磁极、低纬度化极、三分量转换、任意方向分 量转换、磁源重力异常计算；剖面数据处理。
表的基本概念
• 系统通过表的形式将数据与地图有机地 结合在一起。当用户在系统中打开数据 文件时，系统将创建一个表。这个表至 少由两个独立的文件组成，一个是包含 数据结构的文件（*.tab)，另一个是包含 原始数据的文件(*.dat)。关于文件，我们 将在下面描述。
宁夏卫宁北山金场子—二人山岩体重磁资料处理解释
宁夏卫宁北山金场子—二人山岩体重磁资料处理解释朱丹;刘天佑;代小强【摘要】卫宁北山是宁夏重要多金属矿成矿区,渗流卤水—表生改造型金矿及铅、银等多金属矿的形成与深部岩体的作用有关.利用磁异常化极、水平总梯度模法对磁异常进行处理,得出金场子—二人山岩体的平面特征及岩体边界.对过异常中心的东西、南北向两条重磁剖面,采用欧拉齐次方程法、二维磁化率成像和2.5D重磁异常人机交互反演进行反演解释,得出金场子—二人山深部岩体形态及特征,该结果与电法、钻探一致.对于金场子—二人山地区进一步寻找金与铅、银等多金属矿具有指导意义.【期刊名称】《工程地球物理学报》【年(卷),期】2015(012)006【总页数】6页(P766-771)【关键词】卫宁北山;金场子—二人山岩体;重磁异常;多金属矿【作者】朱丹;刘天佑;代小强【作者单位】中国地质大学地球物理与空间信息学院,湖北武汉430074;中国地质大学地球物理与空间信息学院,湖北武汉430074;宁夏地球物理地球化学勘察院,宁夏银川750004【正文语种】中文【中图分类】P631卫宁北山—香山矿集区位于宁夏西部，大地构造位于北祁连走廊过渡带东段与阿拉善地块南缘、鄂尔多斯地块西缘的转换交界处[1]，是宁夏金属矿产勘查重要的地区之一。
卫宁北山位于矿集区的北部，通过多年的勘查工作，发现了多处以金、铅、银等为主的金属矿床，如金场子金矿、照璧山铁矿、大通沟铜矿等[2]。
其中金场子金矿、二人山银铅矿以及一些小型的多金属矿点都主要赋存于石炭系海陆交互相沉积所形成的碳酸盐含膏岩建造中[3]，受东西向断裂带(F18)的控制。
矿床的空间产出部位与硅化、黄钾铁钒化、褐铁矿化、高岭土化和少量的绿泥石化等蚀变体关系密切[4]。
金场子金矿的金矿物为微细金球，其形态为圆球状、椭球状等。
球粒均由多金体组成，成色极高，平均为998，这一特点具有表生金的成色特征[5]，为赋存与巨厚氧化带中的“渗流卤水一表生改造型金矿”(图1)。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第一节 软件功能
一、软件简介
1. 方法原理简介
剖面位场正演拟合是采用二度半多边形截面棱柱体重磁正反演公式[]计算磁
性体模型正演理论曲线，然后与实测异常曲线进行对比，使理论曲线拟合实测曲
线。同时采用奇异值分解与阻尼最小二乘法相结合的方法，得到收敛速度快而且
稳定的计算结果，此方法适合于任意起伏地形条件。
2. 功能简介
剖面位场可视化正演拟合软件是在WINDOWS下开发的具有友好界面的高精
度重磁剖面解释软件。所选模型为水平有限长的棱柱体，截面为任意多边形，其
任意组合可以逼近任意形态的地质体。使用者可以根据实际测量的数据情况，进
行圆滑、滤波预处理及化极、延拓等位场转换处理。根据磁场的曲线形态，可在
计算机屏幕上直观地建立模型，动态地修改模型，且能同时看到其重磁场与实测
场的拟合情况。另外还可以快速直观地反演模型的物性。该程序系统功能强，操
作简便，使用者可以把精力集中于要解决的目标问题上，因而极大地提高了异常
的反演效率和解释效果。
二、主窗口功能介绍
图4.1.1主窗口
如图4.1.1主窗口由上到下由四部分组成，即菜单项、工具条、工作区、状态
条。
图4.1.2 菜单条
1、菜单条
1.1 文件
文件菜单项可以新建模型，对数据文件和模型文件进行装入、保存以及打印
等操作。
1.2 查看
查看菜单项的功能有a.查看部分模型的曲线；b.模型的合理性检查；c.工
具条和状态条的显示/隐藏设置；d.数据区和模型区信息显示。
1.3 编辑修改
编辑修改菜单项可以对模型进行剖分，对窗口进行更新。
1.4 设置
设置菜单项用来对正反演系统进行设置，其设置项包括 a.选择要反演的重
磁场类型；b.设置原始数据曲线、计算数据曲线、数据区、模型区以及所选部分
模型的颜色；c.选择模型移动方式；d.设置地磁场参数以及剖面方位角；e.设置
模型角点加/不加标志；f.设置角点的有效范围。
1.5 预处理
预处理菜单项可以选择性地对原始数据进行三点圆滑、非线性滤波、位场上延、
化极以及调整剖面水平等处理。
1.6 反演
反演菜单项可以通过调整物性约束范围，选择模型进行最优化物性反演。
1.7 格式转换
格式转换菜单项可以对以下四种格式的数据进行转换。a.DOC版数据格式；
b.线数据格式；c.HC-90D格式；d.数据库格式。
1.8 输出
输出菜单项可以以文件的形式输出预处理结果和剖面拟合结果。
1.9 帮助
帮助菜单项可以查到主要功能键的介绍。
2.工具条
图4.1.3 工具条
如图4.1.3，工具条中的功能按钮从左到右分别对应：
2.1 建一新模型文件（清除已有的模型）
2.2 打开（装载）一已有数据模型文件
2.3 保存模型数据文件
2.4 将预处理结果作为原始数据
2.5 设置模型的移动状态为一整体
2.6 设置模型的移动状态为分离
2.7 设置模型的移动状态为拉伸
2.8 是否为模型角点作一有效范围的标志
2.9 对模型一分为二
为满足物性不均匀性的需要，则可对已建模型进行剖分。
2.10 反演模型的物性
当需要反演某些模型的物性时，在反演功能按钮上单击鼠标左键一下，然后
在相应的模型里单击左键，则选中这些模型，再单击一下反演功能按钮，则选择
结束，反演开始。因为采取了一系列措施，使得反演达到实时效果。（此后，这
些模型的物性对话框有反演前后物性变化情况。）
2.11 显示单个或某些模型的异常曲线
先在该按钮上单击鼠标左键，然后选择要显示的模型，选择完毕，则再按一下该
按钮，则会显示所选模型的场值曲线。
2.12 垂向比例尺的选择变化
缺省情况下，界面模型区水平与垂直比例关系是1：1。但有时剖面的物理长度
很大时，例如500公里长的剖面，而一般情装况下，地质体的深度不超过二三十
公里，这时有必要调整水平与垂直之间的比例关系。工具条中有两个按钮，分别
放大与缩小水平与垂直之比。放大、缩小皆是两倍关系，最大为8倍及1/8。
2.13 模型合理性检验(检验模型与地形及飞行线是否交叉)。
2.14 预览输出结果。