重磁异常剖面反演实现
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
摘 要 介绍 了采用截面为任 意多边形 的有限延伸棱柱体为模型 ， 通过曲线拟合 实现重磁剖面联合反演的方法 ， 复杂地质模型 可由不 同
形态 、 密度 、 磁化强度的多个任意多边形组合而成 。
关 键 词 重磁异常 剖面反演 多边形截 面 有限延伸棱柱体
１ 引 言
在物探 资料 的处理 解释 工作 中 ，除 了需要 根据 地
空间上任意一点引起 的磁异常 △Ｔ 表达式为：
△１ Ｈ ＯＳ 。Ｃ ｕ Ｃ １ ＯＳＤ ｏ七 ｎ Ｃ Ｓ ｓ ｎ Ｕ ＋ ￡ ｓｎ １ Ｏ １ ｉ ｉ 。
质图、 物性资料 、 重磁资料对异常源做出定性判断外 ， 还需要运用各种定量反演技术求取异常源的几何参数 和物性参数 ，笔者采用截面为任意多边形的有限延伸 棱柱体为模型 ， 利用 Ｑ ｅ ｇｅ这一跨平 台图形工具 ｔ ｓ ｎｒ Ｄ ｉ 包 ，在 Ｗｉｄ ｗ 和 Ｌｎ ｘ ｎｏ ｓ ｉｕ 系统 上实 现 了重 磁 剖面 ２５ ．Ｄ
ｃａｓ ｏｅ２ ｐ ｂｉＧ ｏ ｇｎ／ 型体类 ｌ ｄｌＤ： ｕｌ Ｐｌ ｏ／模 ｓＭ ｃ ｙ
｛
ｐ ｂｉ： ｕ ｌｃ
３ 功能设置
程序主要功能包括曲线显示和模型显示修改计算 功能两部分 ， 如图 ２ 。
Ｍ ｄｌＤ ； 默认 构造 函数 ｏｅ２ ０／ ／ Ｍ ｄｌＤ ｃｎｔ ｄｌＤ＆ ｔｅＭｏｅ）／ 贝构 造 ｏｅ （ｏｓＭｏｅ ２ ２ ｏｈｒ ｄ１； 拷 ／ 函数
５５ 对 磁性体 而 言 ．
［］ ５ 吴信才 ．ａ ＧＳ地理信息系统【 ． ： Ｍｐ Ｉ Ｍ］ 北京 电子工业 出版 ［ ］ 立新 ， ６吴 史文中． ： 北京 地理信息系统原理 与算法． 科学
出版社 ．０ ３ ２０．
用 △Ｔ ｚ 来判断倾向都是 比较复杂的，这时 或 ａ ［ ］ ｒ ｕ ｒ ｉ Ｄ ／ｔ ７ Ａ ｔ ｒ ｉ ｔ Ｋ Ｅ２ 编程宝典． ｈ Ｇｆｈ Ｑ 北京： 电子工 业出版 社， ０ ． ２ ２ ０ 可 以利 用总磁 场 Ｔ 来 判 断。根 据规格 化公式 得 ： ａ
Ｑ ｔ ｉ＜ ｏｅ Ｄ ｍ ｐ ｏｅ ｌｔ 模型体 Ｐ Ｌｓ Ｍ ｄｌ ＞＊ ＿ Ｍ ｄｌｉ； ｒ ｔ ２ ＿ ｓ／ ／
Ｑ ｔ ｉ＜ ｏ ｔｄ ｍ ｐ ｅ ｈ ｌｔ／ Ｐｒ ｓ Ｐｉ ２＞＊ Ｈ ｉ ｔｉ； 高程 Ｌｔ ｎ ｇ ＿ ｓ／
立初始模型并输人物性参数一重磁数值计算一对 比拟 合曲线与实测观测异常曲线 与实测观测异常曲线残差最小－ 输出 Ｍ ｐ Ｉ 结果。 ＋ ａＧＳ
量计算 ， 如异常形态近似圆形 ， 可选择 １ 条剖面 ， 异常 重磁异常剖面反演是一种可同时运行在 Ｗｉ ｏ ｓ ｎ ｗ ｄ 形态近似椭圆形可按照长轴、 短轴方向各选择一条； 选 和 Ｌｎｘ系统上 的重 磁 异常解 释 工具 ，已在 资源 潜力 ｉｕ 取剖面 的长度应 取 到异常 两侧正 常场 。 评价 等项 目上得 到 了广 泛 的应用 ，程序 除 了上述 基本
设置 模型指 针
ｖ ｉ ｓｔ ｅ ｈ（ Ｐｒｉ＜ ｏ ｔｄ ｐ ｅｇｔ∥设 置 ｏ Ｈ ｉ ｔＱ ｔ ｓ Ｐ ｉ ２ ＞＊ Ｈ ｉ ） ｄｅ ｇ Ｌｔ ｎ ｈ
删除模型 ： 已有的模型体删除； 将
模 型上 删点 ： 已有 的模 型节 点删除 ； 将
高程
移动模型及节点 ： 移动模型体或节点 ； 修改模型属性 ： 修改模型的物性参数。
ｖ ｕｌｏｌｄ Ｎ ｄ（Ｐｉｔｖ ｓ；增 加节点 ｉ ａｂｏａｄ ｏｅＱ ｏｎａ ｏ ）／ ｔ ｒ ／
ｖ ｕｌｏｌ ｅ ｏｅ Ｐｉ ＆ ｏ ） 删除节点 ｉ ａｂｏｄｌ ｄ（ ｏ ｔ ｐｓ； ｔ ｒ Ｎ Ｑ ｎ ， ，
ｖｔ ｌｏ ａ （Ｐｉ ｅ ＊ａ ｔ ／ ｉｕ ｉ ｄ ｗＱ ａ ｔ ｐｉ ） 绘制模型体 ｒ ａｖ ｄ ｒ ｎ ｒ ｎＪ
ｄｕｌ ｃ ｃｌｅ ｅａ０／ ｏｂ ｌ ａ Ｄ ｌＴ ； 计算 △Ｔ ｅａｕｔ ｔ ／
ｄｕ ｌ ｃ ｃ ｌ ｅ ｒｖｔ ； 计 算重力 ｏｂｅ ａ ｕａ Ｇ ａｉ ０／ ｌ ｔ ｙ ／
ｐｒｖ ｔ ： ｉ ａｅ
４ 程序实现
一
个完整的反演过程为 ：输入重磁剖面异常 建 改模型直到拟合曲线
） 嚣 －ｑ ＝ ｃ 酬ｎ ｓ（ 一啪 ｉ一 ） ｎ ￣
＝ｓｎ￣ｇｎ ， ｙ ＋ＣＯ ｄ ｉ Ｒ ＋ Ｉ Ｓｐ
＋
￣
（ａ 一组 ａｎ ａｎ ｒ ｒ嚣） ｃ ｃ ｃ
∞
ｎ
图 １ 多边 形 棱柱 体 模 型 示 意 图
上述模 型 的重力 异常 经简化后 ，则 在空 间上任 意
一
＿
清除图片： 清除模型区的图片；
网格设置按 钮 ， 置模 型 区坐标 网格 ； 设源自文库
｝ ；
图截异常： Ｇ Ｄ数据上截取剖面作为实测数据； 从 Ｒ 新建模型： 在模型区左键点击 ３ 个点或 ３ 个点以上， 右键结束， 可形成封闭的模型区， 入 输 模型的里端延伸长 度、 外端延伸长度 ， 磁法反演需输人地质体的磁化强度、 磁倾角、 磁偏角， 重力反演需要输 入 地质体的密度。 模 型加点： 在已有的模型上增加节点 ；
５３ 物性 参数 选取 ． 对 于磁性 模 型体而 言磁化 强度 应为岩 矿石 的感应
参考 文献
［］ １ 管志宁． 地磁场与磁力勘探 ． 北京 ： 地质出版社 ， ０ ． ２ ５ ０
［］ ２ 刘天佑． 地球物理勘探概论． 北京 ： 质出版社． ０ ． 地 ２ ７ ０
磁化强度与剩余磁化强度 的矢量合成；磁倾角为岩矿 石的剩余磁化方向与地磁场方向的矢量合成，对大多 数岩矿石来说 ， 由于剩余磁化强度很弱 ， 可以忽略剩余 磁化方向的影响， 直接用地磁场方向作为其磁化倾角。
５２ 异 常边界确 定 ．
功能外还有保存现场 、 图截异常、 体积法储量计算 、 生 异常边界可 以用重磁异常在等值线图上的梯度陡 成模 型产状 、 物性参 数 等实用 功 能 ， 能有 效 的推 断 目标
变带来确定 ；对于规模较小的地质体可用垂向一阶导 地 质体 空 间展 布形 态及 物性参 数 ，提高 重磁异 常 的反 数零值线来圈定 ，对于异常规模较大的地质体可采用 演 效率 和解 释效果 。 垂 向二 阶导 数零值 线来 圈定 。
ｐ ｉ ａ ｅ ｒｖ ｔ ：
图 ２ 重 磁 联 合 反 演界 面
Ｑ ｔ ｉ＜ ｎｍ ｌＮ ｄ２ ５ ｍ ｐ ｏｅｌｔ／ 常 Ｐｒ ｓ Ａ ｏ ａ ｏｅＤ ＞＊ Ｎ ｄ＿ｉ； 异 Ｌｔ ｙ ｓ／
剖面曲线和模型操作 的主要功能如下 ： 打开模型 ： 装入保存的模型体文件；
体节点列表
ｄｏ ｂ ｅｍ ｕ ｌ
＝ 一
磁 场三 分量为 ：
＝ 一
莩 ￣ ） ｓｉ ／ ｉ ＋ ＋ ｎ （
Ｎ
反演程序的开发 ， 程序操作简单 、 功能强 ， 生成 的结果 与常用软件（ Ｍ ｐ Ｉ） 如 ａ ＧＳ完全兼容 ， 大大地提高了重磁
异常 的反演效 率和解 释效 果 。
Ｎ
一
ｓ
ｚ ｏ ｆ￣ ｓ ｃ
ｌ ８
段新力 唐群英
毕武
袁小龙
黄显义
彭仲秋
向诗强 ： 重磁异常剖面反演实现
增刊 ｌ
以上各式 中符号意义如下 ： Ｇ为引力常数 ，为棱 在上 述流程 中对模 型 的操 作 ， ｉ 即模 型建立 、 修改和重 磁 柱体角点标号 ， Ｎ为棱柱 体得边数 ， 、 ０ １ Ｄ 为地磁场 的 异常计算是拟合工作的重要环节，程序定义的模型和 ０ 倾角 、 偏角 ，、 ＩＤ为磁化强度方向的倾角、 偏角 ， 当场源 异常计算类如下： 是由多个多边形组合而成时， 异常可叠加而得 。
２ １ 年 ０１
新 疆 有 色 金 属
１ ９
５ 几点体会
５１ 剖面选 取 ．
体， 不对称时，ａ Ｔ 曲线缓的一侧为磁性体倾 向方向。
５６ 重 视残 余异 常 ．
根 据 已知 矿体 的空 间展 布形态 、矿体 和 围岩 的物 作为定量解释的剖面数据最好使用接近异常中心 性参数正演得到矿体和围岩引起的重磁异常 ，将实测 的原始 测量 剖面数 据 ，当原始 测量 的剖 面方 向与重 磁 重磁异常与其的差 ， 作为残余异常 ， 对残余异常再进行 异常走 向的夹角较小时（ ７。 ） ＜０ ， 应直接从等值线图垂 反演计算 ，可 以判断已知矿体深部或者附近是否还存 直于异常走向的位置截取通过异常 中心的剖面数据 ； 在 矿体 。 对有一定走向的地质体进行计算时 ，选取的剖面不少 于 ２条 ， 对没 有 明确走 向的地 质 体 可选 １２条 进 行定 ６ 结束语 —
Ｍｏ ｅ Ｄ＆ ｐｒｔｒ ｃｎｔ ｄ ｌＤ＆ｏ ｅＭｏ ｅ） ｄ ｌ ｏ ｅａ ＝（ ｓ Ｍｏｅ２ ２ ｏ ｏ ｔ ｒ ｄ １； ｈ
，等号运 算符 ／
ｖ ｔ ｌｏ ｅｔＱ ｔ ｉ＜ ｏ ｔ ＞＊Ｎ ｄｓ， 创 ｉｕ ｉ ｃ ａ （Ｐ Ｌ ｔ ｉ２ ｐ ｏｅ ） ｒ ａｖ ｄ ｒ ｅ ｒ ｓ Ｐ ｎ ｄ － ／ ／ 建模型体 ｖｔ ｌｏ ｏｅ Ｐ ｉ ｏ ） 移动模型体 ｉｕ ｉ ｍ ｖ（ ｏ ｔｖ ｓ； ｒ ａｖ ｄ Ｑ ｎａ ／ ／
Ｍ ｓ １ ｃ ｓ ３］ （ｉ ｎ ｆ ｏ ｆ — ）
２ 算法模型
本程序采用的截面为任意多边形的有限延伸棱柱
体模 型如 图 １ 示 ： 所
＝
（
＋ １ ＋ ，ｆ 肘ｙ２ ３ ｉ ）
其 中：
１ ＝ （ ） （ ） ２＝ （ 一 ， ， （ 一 （ ） ｙ 一 ， 。 ） （ １ ＝ ） ， ， ２ １ ） ３
＿
ｒｏ ／ ｈ ； 密度 ／ Ｍ； 磁 化强度 ／ ／
保存数据 ： 保存模型及反演结果 ； 保存工程 ： 保存工作现场 ；
打 开 图片 ： 装入 地质 剖面 图并 校正 ；
ｄ ｕ ｌ ｏ ｂ ｅｍ
＿
ｄ ｕ ｌ ｏ ｂ ｅｒ ｅ
ｄ ｕｂ ｅｍ ｏ ｌ
＿
Ｉ 磁倾角 ； ／ ／
Ｄ／磁偏 角 ； ／
ｃ ｓ ｎｍ ｙｒ ｔ ／ ｌ ｏ ａ Ｃｅ ｏ／ ｓ ａ Ａ ｌ ａ ｒ 重磁异常计算类
｛
ｐ ｂｌ ： ｕ ｉ ｃ
。 。
Ａ ｏ ａ Ｃｅｔ ０ ／ ｎｍ ｙ ｒ ｏ ； 默认构造函数 ｌ ａｒ ／
ｖｉ ｓｔ ｄｌｉ （ Ｐｒｉ ＜ ｄ ｌＤ ｐ ｏ ｅ ） ｏｄ ｅ Ｍｏ ｅ ｓ Ｑ ｔ ｓ Ｍｏ ｅ ＞ Ｍ ｄｌ ∥ Ｌｔ Ｌｔ ２ ｓ
５４ 模型体 的初 始几何 参数 ．
得。
［］ ３ 张胜业， 潘玉玲． 用地球物理学 原理． 应 武汉 ： 中国地 质
大学出版社 ，０ ４ ２０ ．
［ ］中华人 民共 和国国土资源部 Ｄ ／０ ８ — ０ ６ ４ ＺＴ ０ ２ ２０ ．区域重
力调查规范， ０ ． ２ ６ ０
顶部埋深、中心位置、宽体可通过切线法计算获 社 ， ０４ ２０ ．
，ｌ ＋ ＝ ｆ１ ｏ ｂ ＋ ｃ ｓｑ ＋ ｚ ｉ ｂ ＨＩ ｎ￣ ｓ
ｓｑ ｗｌ ｓｎ ＋ｚｃｏ ｂ ＝一 ｉ
Ｉ
，
点引起的重力异常 △ （表达式为 ： ｇ） ｒ
ｒ＝（ ＋ ） “ “ 。 “ ＋ ＝（ ） “
Ａｇ ＝ ｃｓ２ＩＹ ，＋Ｏ￣） （ ｃａ Ｏ ｏ￣ ［ ２）Ｉｚ，１ Ｉ （ ｉ ｉ
２ １ 矩 ０１
新 疆 有 色 金 属
１ ７
重 磁 异 常 剖 面 反 演 实 现
段新力① 唐群英② 毕 武① 袁小龙③ 黄显义③ 彭仲秋③ 向诗强③
（ 新疆地矿局物化探大队 昌吉 ８０９ ② 新疆矿产实验研究所 乌鲁木齐 ８００ ① ３０ １ ３００
③ 乌鲁木 齐金 维 图文信 息科 技 有限公 司 乌鲁 木 齐 ８０ １） ３０ １
： ＋ ） ｙ＋ “ 其 中， Ｙｎ ｌ 十 ＿ｎ ｕ１ ＋ ｌ 一 ．旦 ｕｎ １
Ｍ ：Ｍ ｃ ｓ ｃ ｓ ｏ Ｉ ｏ Ｄ
Ｒ ＝ Ｊ ｙ ＋ ） ｉ （２ 嵋Ｉ ｌ＋ ＋ ”
Ｍ ＝Ｍ ｃ ｓ ｉ ｏ ｌｓ ｎＤ
－ 警 ｗ ｉ
刮ｇ ） 警
Ｍ ． ＝Ｍ ｎＩ ｓ ｉ
形态 、 密度 、 磁化强度的多个任意多边形组合而成 。
关 键 词 重磁异常 剖面反演 多边形截 面 有限延伸棱柱体
１ 引 言
在物探 资料 的处理 解释 工作 中 ，除 了需要 根据 地
空间上任意一点引起 的磁异常 △Ｔ 表达式为：
△１ Ｈ ＯＳ 。Ｃ ｕ Ｃ １ ＯＳＤ ｏ七 ｎ Ｃ Ｓ ｓ ｎ Ｕ ＋ ￡ ｓｎ １ Ｏ １ ｉ ｉ 。
质图、 物性资料 、 重磁资料对异常源做出定性判断外 ， 还需要运用各种定量反演技术求取异常源的几何参数 和物性参数 ，笔者采用截面为任意多边形的有限延伸 棱柱体为模型 ， 利用 Ｑ ｅ ｇｅ这一跨平 台图形工具 ｔ ｓ ｎｒ Ｄ ｉ 包 ，在 Ｗｉｄ ｗ 和 Ｌｎ ｘ ｎｏ ｓ ｉｕ 系统 上实 现 了重 磁 剖面 ２５ ．Ｄ
ｃａｓ ｏｅ２ ｐ ｂｉＧ ｏ ｇｎ／ 型体类 ｌ ｄｌＤ： ｕｌ Ｐｌ ｏ／模 ｓＭ ｃ ｙ
｛
ｐ ｂｉ： ｕ ｌｃ
３ 功能设置
程序主要功能包括曲线显示和模型显示修改计算 功能两部分 ， 如图 ２ 。
Ｍ ｄｌＤ ； 默认 构造 函数 ｏｅ２ ０／ ／ Ｍ ｄｌＤ ｃｎｔ ｄｌＤ＆ ｔｅＭｏｅ）／ 贝构 造 ｏｅ （ｏｓＭｏｅ ２ ２ ｏｈｒ ｄ１； 拷 ／ 函数
５５ 对 磁性体 而 言 ．
［］ ５ 吴信才 ．ａ ＧＳ地理信息系统【 ． ： Ｍｐ Ｉ Ｍ］ 北京 电子工业 出版 ［ ］ 立新 ， ６吴 史文中． ： 北京 地理信息系统原理 与算法． 科学
出版社 ．０ ３ ２０．
用 △Ｔ ｚ 来判断倾向都是 比较复杂的，这时 或 ａ ［ ］ ｒ ｕ ｒ ｉ Ｄ ／ｔ ７ Ａ ｔ ｒ ｉ ｔ Ｋ Ｅ２ 编程宝典． ｈ Ｇｆｈ Ｑ 北京： 电子工 业出版 社， ０ ． ２ ２ ０ 可 以利 用总磁 场 Ｔ 来 判 断。根 据规格 化公式 得 ： ａ
Ｑ ｔ ｉ＜ ｏｅ Ｄ ｍ ｐ ｏｅ ｌｔ 模型体 Ｐ Ｌｓ Ｍ ｄｌ ＞＊ ＿ Ｍ ｄｌｉ； ｒ ｔ ２ ＿ ｓ／ ／
Ｑ ｔ ｉ＜ ｏ ｔｄ ｍ ｐ ｅ ｈ ｌｔ／ Ｐｒ ｓ Ｐｉ ２＞＊ Ｈ ｉ ｔｉ； 高程 Ｌｔ ｎ ｇ ＿ ｓ／
立初始模型并输人物性参数一重磁数值计算一对 比拟 合曲线与实测观测异常曲线 与实测观测异常曲线残差最小－ 输出 Ｍ ｐ Ｉ 结果。 ＋ ａＧＳ
量计算 ， 如异常形态近似圆形 ， 可选择 １ 条剖面 ， 异常 重磁异常剖面反演是一种可同时运行在 Ｗｉ ｏ ｓ ｎ ｗ ｄ 形态近似椭圆形可按照长轴、 短轴方向各选择一条； 选 和 Ｌｎｘ系统上 的重 磁 异常解 释 工具 ，已在 资源 潜力 ｉｕ 取剖面 的长度应 取 到异常 两侧正 常场 。 评价 等项 目上得 到 了广 泛 的应用 ，程序 除 了上述 基本
设置 模型指 针
ｖ ｉ ｓｔ ｅ ｈ（ Ｐｒｉ＜ ｏ ｔｄ ｐ ｅｇｔ∥设 置 ｏ Ｈ ｉ ｔＱ ｔ ｓ Ｐ ｉ ２ ＞＊ Ｈ ｉ ） ｄｅ ｇ Ｌｔ ｎ ｈ
删除模型 ： 已有的模型体删除； 将
模 型上 删点 ： 已有 的模 型节 点删除 ； 将
高程
移动模型及节点 ： 移动模型体或节点 ； 修改模型属性 ： 修改模型的物性参数。
ｖ ｕｌｏｌｄ Ｎ ｄ（Ｐｉｔｖ ｓ；增 加节点 ｉ ａｂｏａｄ ｏｅＱ ｏｎａ ｏ ）／ ｔ ｒ ／
ｖ ｕｌｏｌ ｅ ｏｅ Ｐｉ ＆ ｏ ） 删除节点 ｉ ａｂｏｄｌ ｄ（ ｏ ｔ ｐｓ； ｔ ｒ Ｎ Ｑ ｎ ， ，
ｖｔ ｌｏ ａ （Ｐｉ ｅ ＊ａ ｔ ／ ｉｕ ｉ ｄ ｗＱ ａ ｔ ｐｉ ） 绘制模型体 ｒ ａｖ ｄ ｒ ｎ ｒ ｎＪ
ｄｕｌ ｃ ｃｌｅ ｅａ０／ ｏｂ ｌ ａ Ｄ ｌＴ ； 计算 △Ｔ ｅａｕｔ ｔ ／
ｄｕ ｌ ｃ ｃ ｌ ｅ ｒｖｔ ； 计 算重力 ｏｂｅ ａ ｕａ Ｇ ａｉ ０／ ｌ ｔ ｙ ／
ｐｒｖ ｔ ： ｉ ａｅ
４ 程序实现
一
个完整的反演过程为 ：输入重磁剖面异常 建 改模型直到拟合曲线
） 嚣 －ｑ ＝ ｃ 酬ｎ ｓ（ 一啪 ｉ一 ） ｎ ￣
＝ｓｎ￣ｇｎ ， ｙ ＋ＣＯ ｄ ｉ Ｒ ＋ Ｉ Ｓｐ
＋
￣
（ａ 一组 ａｎ ａｎ ｒ ｒ嚣） ｃ ｃ ｃ
∞
ｎ
图 １ 多边 形 棱柱 体 模 型 示 意 图
上述模 型 的重力 异常 经简化后 ，则 在空 间上任 意
一
＿
清除图片： 清除模型区的图片；
网格设置按 钮 ， 置模 型 区坐标 网格 ； 设源自文库
｝ ；
图截异常： Ｇ Ｄ数据上截取剖面作为实测数据； 从 Ｒ 新建模型： 在模型区左键点击 ３ 个点或 ３ 个点以上， 右键结束， 可形成封闭的模型区， 入 输 模型的里端延伸长 度、 外端延伸长度 ， 磁法反演需输人地质体的磁化强度、 磁倾角、 磁偏角， 重力反演需要输 入 地质体的密度。 模 型加点： 在已有的模型上增加节点 ；
５３ 物性 参数 选取 ． 对 于磁性 模 型体而 言磁化 强度 应为岩 矿石 的感应
参考 文献
［］ １ 管志宁． 地磁场与磁力勘探 ． 北京 ： 地质出版社 ， ０ ． ２ ５ ０
［］ ２ 刘天佑． 地球物理勘探概论． 北京 ： 质出版社． ０ ． 地 ２ ７ ０
磁化强度与剩余磁化强度 的矢量合成；磁倾角为岩矿 石的剩余磁化方向与地磁场方向的矢量合成，对大多 数岩矿石来说 ， 由于剩余磁化强度很弱 ， 可以忽略剩余 磁化方向的影响， 直接用地磁场方向作为其磁化倾角。
５２ 异 常边界确 定 ．
功能外还有保存现场 、 图截异常、 体积法储量计算 、 生 异常边界可 以用重磁异常在等值线图上的梯度陡 成模 型产状 、 物性参 数 等实用 功 能 ， 能有 效 的推 断 目标
变带来确定 ；对于规模较小的地质体可用垂向一阶导 地 质体 空 间展 布形 态及 物性参 数 ，提高 重磁异 常 的反 数零值线来圈定 ，对于异常规模较大的地质体可采用 演 效率 和解 释效果 。 垂 向二 阶导 数零值 线来 圈定 。
ｐ ｉ ａ ｅ ｒｖ ｔ ：
图 ２ 重 磁 联 合 反 演界 面
Ｑ ｔ ｉ＜ ｎｍ ｌＮ ｄ２ ５ ｍ ｐ ｏｅｌｔ／ 常 Ｐｒ ｓ Ａ ｏ ａ ｏｅＤ ＞＊ Ｎ ｄ＿ｉ； 异 Ｌｔ ｙ ｓ／
剖面曲线和模型操作 的主要功能如下 ： 打开模型 ： 装入保存的模型体文件；
体节点列表
ｄｏ ｂ ｅｍ ｕ ｌ
＝ 一
磁 场三 分量为 ：
＝ 一
莩 ￣ ） ｓｉ ／ ｉ ＋ ＋ ｎ （
Ｎ
反演程序的开发 ， 程序操作简单 、 功能强 ， 生成 的结果 与常用软件（ Ｍ ｐ Ｉ） 如 ａ ＧＳ完全兼容 ， 大大地提高了重磁
异常 的反演效 率和解 释效 果 。
Ｎ
一
ｓ
ｚ ｏ ｆ￣ ｓ ｃ
ｌ ８
段新力 唐群英
毕武
袁小龙
黄显义
彭仲秋
向诗强 ： 重磁异常剖面反演实现
增刊 ｌ
以上各式 中符号意义如下 ： Ｇ为引力常数 ，为棱 在上 述流程 中对模 型 的操 作 ， ｉ 即模 型建立 、 修改和重 磁 柱体角点标号 ， Ｎ为棱柱 体得边数 ， 、 ０ １ Ｄ 为地磁场 的 异常计算是拟合工作的重要环节，程序定义的模型和 ０ 倾角 、 偏角 ，、 ＩＤ为磁化强度方向的倾角、 偏角 ， 当场源 异常计算类如下： 是由多个多边形组合而成时， 异常可叠加而得 。
２ １ 年 ０１
新 疆 有 色 金 属
１ ９
５ 几点体会
５１ 剖面选 取 ．
体， 不对称时，ａ Ｔ 曲线缓的一侧为磁性体倾 向方向。
５６ 重 视残 余异 常 ．
根 据 已知 矿体 的空 间展 布形态 、矿体 和 围岩 的物 作为定量解释的剖面数据最好使用接近异常中心 性参数正演得到矿体和围岩引起的重磁异常 ，将实测 的原始 测量 剖面数 据 ，当原始 测量 的剖 面方 向与重 磁 重磁异常与其的差 ， 作为残余异常 ， 对残余异常再进行 异常走 向的夹角较小时（ ７。 ） ＜０ ， 应直接从等值线图垂 反演计算 ，可 以判断已知矿体深部或者附近是否还存 直于异常走向的位置截取通过异常 中心的剖面数据 ； 在 矿体 。 对有一定走向的地质体进行计算时 ，选取的剖面不少 于 ２条 ， 对没 有 明确走 向的地 质 体 可选 １２条 进 行定 ６ 结束语 —
Ｍｏ ｅ Ｄ＆ ｐｒｔｒ ｃｎｔ ｄ ｌＤ＆ｏ ｅＭｏ ｅ） ｄ ｌ ｏ ｅａ ＝（ ｓ Ｍｏｅ２ ２ ｏ ｏ ｔ ｒ ｄ １； ｈ
，等号运 算符 ／
ｖ ｔ ｌｏ ｅｔＱ ｔ ｉ＜ ｏ ｔ ＞＊Ｎ ｄｓ， 创 ｉｕ ｉ ｃ ａ （Ｐ Ｌ ｔ ｉ２ ｐ ｏｅ ） ｒ ａｖ ｄ ｒ ｅ ｒ ｓ Ｐ ｎ ｄ － ／ ／ 建模型体 ｖｔ ｌｏ ｏｅ Ｐ ｉ ｏ ） 移动模型体 ｉｕ ｉ ｍ ｖ（ ｏ ｔｖ ｓ； ｒ ａｖ ｄ Ｑ ｎａ ／ ／
Ｍ ｓ １ ｃ ｓ ３］ （ｉ ｎ ｆ ｏ ｆ — ）
２ 算法模型
本程序采用的截面为任意多边形的有限延伸棱柱
体模 型如 图 １ 示 ： 所
＝
（
＋ １ ＋ ，ｆ 肘ｙ２ ３ ｉ ）
其 中：
１ ＝ （ ） （ ） ２＝ （ 一 ， ， （ 一 （ ） ｙ 一 ， 。 ） （ １ ＝ ） ， ， ２ １ ） ３
＿
ｒｏ ／ ｈ ； 密度 ／ Ｍ； 磁 化强度 ／ ／
保存数据 ： 保存模型及反演结果 ； 保存工程 ： 保存工作现场 ；
打 开 图片 ： 装入 地质 剖面 图并 校正 ；
ｄ ｕ ｌ ｏ ｂ ｅｍ
＿
ｄ ｕ ｌ ｏ ｂ ｅｒ ｅ
ｄ ｕｂ ｅｍ ｏ ｌ
＿
Ｉ 磁倾角 ； ／ ／
Ｄ／磁偏 角 ； ／
ｃ ｓ ｎｍ ｙｒ ｔ ／ ｌ ｏ ａ Ｃｅ ｏ／ ｓ ａ Ａ ｌ ａ ｒ 重磁异常计算类
｛
ｐ ｂｌ ： ｕ ｉ ｃ
。 。
Ａ ｏ ａ Ｃｅｔ ０ ／ ｎｍ ｙ ｒ ｏ ； 默认构造函数 ｌ ａｒ ／
ｖｉ ｓｔ ｄｌｉ （ Ｐｒｉ ＜ ｄ ｌＤ ｐ ｏ ｅ ） ｏｄ ｅ Ｍｏ ｅ ｓ Ｑ ｔ ｓ Ｍｏ ｅ ＞ Ｍ ｄｌ ∥ Ｌｔ Ｌｔ ２ ｓ
５４ 模型体 的初 始几何 参数 ．
得。
［］ ３ 张胜业， 潘玉玲． 用地球物理学 原理． 应 武汉 ： 中国地 质
大学出版社 ，０ ４ ２０ ．
［ ］中华人 民共 和国国土资源部 Ｄ ／０ ８ — ０ ６ ４ ＺＴ ０ ２ ２０ ．区域重
力调查规范， ０ ． ２ ６ ０
顶部埋深、中心位置、宽体可通过切线法计算获 社 ， ０４ ２０ ．
，ｌ ＋ ＝ ｆ１ ｏ ｂ ＋ ｃ ｓｑ ＋ ｚ ｉ ｂ ＨＩ ｎ￣ ｓ
ｓｑ ｗｌ ｓｎ ＋ｚｃｏ ｂ ＝一 ｉ
Ｉ
，
点引起的重力异常 △ （表达式为 ： ｇ） ｒ
ｒ＝（ ＋ ） “ “ 。 “ ＋ ＝（ ） “
Ａｇ ＝ ｃｓ２ＩＹ ，＋Ｏ￣） （ ｃａ Ｏ ｏ￣ ［ ２）Ｉｚ，１ Ｉ （ ｉ ｉ
２ １ 矩 ０１
新 疆 有 色 金 属
１ ７
重 磁 异 常 剖 面 反 演 实 现
段新力① 唐群英② 毕 武① 袁小龙③ 黄显义③ 彭仲秋③ 向诗强③
（ 新疆地矿局物化探大队 昌吉 ８０９ ② 新疆矿产实验研究所 乌鲁木齐 ８００ ① ３０ １ ３００
③ 乌鲁木 齐金 维 图文信 息科 技 有限公 司 乌鲁 木 齐 ８０ １） ３０ １
： ＋ ） ｙ＋ “ 其 中， Ｙｎ ｌ 十 ＿ｎ ｕ１ ＋ ｌ 一 ．旦 ｕｎ １
Ｍ ：Ｍ ｃ ｓ ｃ ｓ ｏ Ｉ ｏ Ｄ
Ｒ ＝ Ｊ ｙ ＋ ） ｉ （２ 嵋Ｉ ｌ＋ ＋ ”
Ｍ ＝Ｍ ｃ ｓ ｉ ｏ ｌｓ ｎＤ
－ 警 ｗ ｉ
刮ｇ ） 警
Ｍ ． ＝Ｍ ｎＩ ｓ ｉ