磁异常ΔT三维相关成像
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
［４］ 积法的实用性 ． 在常规欧拉反褶积 ２ ０ ０ ７年 Ｓ ａ ｌ ｅ ｍ１
犕狇 ０ ２ ） 犜 狓， 狕）＝ μ · ５ 犃·（ 狓 Δ 狔， 狇（ 狇 －狓 ４ π 狉 ２ ） ） ） ·（ 狓 狕 ＋犅·（ ＋犆·（ 狔 狇 －狔 狇 －狓 狇 －狕
［
） ） ） ·（ 狕 狓 ＋犈·（ ＋ 犇·（ 狔 狇 －狔 狇 －狕 狇 －狓
犜ｄ ａ ｔ ａｏ ｆｓ ｎ ｔ ｈ ｅ ｓ ｉ ｚ ｅ ｄｍ ｕ ｌ ｔ ｉ ｌ ｅｒ ｅ ｃ ｔ ａ ｎ ｕ ｌ ａ ｒｐ ｒ ｉ ｓ ｍ ｓｍ ｏ ｄ ｅ ｌａ ｎ ｄｔ ｈ ｅ ｔ ｅ ｓ ｔｏ ｎｔ ｈ ｅｍ ａ ｎ ｅ ｔ ｉ ｃａ ｎ ｏ ｍ ａ ｌ ｙ ｐ ｇ ｇ ｙΔ ， ｒ ａ ｗｍ ａ ｎ ｅ ｔ ｉ ｃｄ ａ ｔ ａ ｓ ｈ ｏ ｗｔ ｈ ａ ｔ ｔ ｈ ｅｍ ｅ ｔ ｈ ｏ ｄｃ ａ ｎ ｉ ｍ ａ ｅ ｔ ｈ ｅ ｓ ａ ｃ ｅ ｅ ｘ ｔ ｅ ｎ ｄ ｅ ｄｓ ｔ ａ ｔ ｅ ｓ ａ ｎ ｄ ｔ ｈ ｅｄ ｉ ｓ ｔ ｒ ｉ ｂ ｕ ｔ ｉ ｏ ｎ ｇ ｇ ｐ
３ 犇犮 狅 狉 狉 犲 犾 犪 狋 犻 狅 狀 犻 犿 犪 犻 狀 犳 狅 狉犿 犪 狀 犲 狋 犻 犮犪 狀 狅 犿 犪 犾 犜犱 犪 狋 犪 犵 犵 犵 狔Δ
１２ １２ １２ ，ＭＥ ＧＵＯ Ｌ ｉ ａ ｎ Ｈ ｕ ｉ ＮＧ Ｘ ｉ ａ ｏ Ｈ ｏ ｎ Ｓ Ｈ ＩＬ ｅ ｉ ｇ ｇ ， ， ， ，
利
、 Ｐ ｉ ｌ ｋ ｉ ｎ ｔ ｏ ｎ 等提出了一些改进方法和措 ｇ 在一定程度上改善了反演效果和效率 ． 施，
］ １ ３ 磁异常欧拉 反 褶 积 方 法 ［ 以欧拉齐次方程为
［ １ ０， １ １］
［ ］ １ ２
运用磁异常 、 空间导数以及各种地质体具有特 基础 ， 定的 “ 构造指数 ” 来 确 定 异 常 场 源 的 位 置． 它无需场 计算简便快速 ， 适合于大面积磁 源物性的先验知识 ， 但是实际场源构造指数的未知 测资料的分析解释 ． 性、 反褶积解的发散 性 等 实 际 问 题 制 约 着 欧 拉 反 褶
４ ３ ６
地 球 物 理 学 报（ ） Ｃ ｈ ｉ ｎ ｅ ｓ ｅＪ ． Ｇ ｅ ｏ ｈ ｓ ． ｐ ｙ
５ ３卷
实际地质情况的物性分布和几何形态 ． 几十年来 ， 国 内外众多学者
［ １～７］
在基准面上 ， 假设测区地下任意 狕 取垂直向下为正 ， 坐标为 （ 、 体积为 狏 狓 狕 狔 狇， 狇， 狇） 狇 的第狇 个均匀磁化 球 体的磁化强度为犑狇 ， 磁矩 犕狇 ＝ 犑 地磁场倾角为 狏 狇 狇， 偏角为 犃 磁化 强度 倾 角 为 犐、 偏角为 犃 由于 犐 ′ ′． ０、 ０， 均匀磁化球体外部磁场与位于球心的偶极子磁场等 因此 ， 这 里 把 均 匀 磁 化 球 体 视 为 偶 极 子． 则 第狇 效， 个偶极子在测区上任意点 （ 狓， 狕）处的 总强 度磁 异 狔，
·犔 －狀·犖 － 犿·犕 ， 犾 犃 ＝ ２· ·犔 －狀·犖 ， 犅 ＝ ２·犿·犕 －犾 ， ·犖 ＋狀·犔） 犾 犆 ＝ ３·（ ， 犿·犖 ＋狀·犕 ） 犇 ＝ ３·（ ， ·犕 ＋ 犿·犔） 犾 犈 ＝ ３·（ ·犔． 狀·犖 － 犿·犕 －犾 犉 ＝ ２· 定义测区上实测 Δ 犜 异常与第狇 个 偶极 子测 区 ） 互相关系数为 上的 Δ 犜 Ｐ ｅ ａ ｒ ｓ ｏ ｎ 狇 异常的皮氏 （
第５ ３卷 第２期 ２ ０ １ ０年２月
地 球 物 理 学 报
ＣＨ Ｉ Ｎ Ｅ Ｓ Ｅ Ｊ ＯＵＲ ＮＡ Ｌ Ｏ Ｆ Ｇ Ｅ Ｏ Ｐ ＨＹ Ｓ Ｉ Ｃ Ｓ
Ｖ ｏ ｌ ． ５ ３，Ｎ ｏ ． ２ ， Ｆ ｅ ｂ ． ２ ０ １ ０
郭良辉 ， 孟小红 ， 石 磊． 磁异常 Δ 地 球 物 理 学 报， ： ： ／ 犜 三 维 相 关 成 像． ２ ０ １ ０， ５ ３（ ２） ４ ３ ５～４ ４ １， Ｄ Ｏ Ｉ １ ０． ３ ９ ６ ９ ． ｉ ｓ ｓ ｎ ． ０ ０ ０ １ ５ ７ ３ ３． ｊ ２ ０ １ ０． ０ ２． ０ ２ ２ ） ， Ｇ ｕ ｏＬ Ｈ，Ｍ ｅ ｎ Ｓ ｈ ｉ Ｌ． ３ Ｄｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｍ ａ ｉ ｎ ｆ ｏ ｒｍ ａ ｎ ｅ ｔ ｉ ｃ ａ ｎ ｏ ｍ ａ ｌ 犜ｄ ａ ｔ ａ ． 犆 犺 犻 狀 犲 狊 犲 犑． 犌 犲 狅 犺 狊 ．（ ｉ ｎＣ ｈ ｉ ｎ ｅ ｓ ｅ ２ ０ １ ０， ｇＸ Ｈ， ｇ ｇ ｇ ｙΔ 狆 狔 （ ） ： ： ／ ５ ３ ２ ４ ３ ５～４ ４ １， Ｄ Ｏ Ｉ １ ０． ３ ９ ６ ９ ． ｉ ｓ ｓ ｎ ． ０ ０ ０ １ ５ ７ ３ ３． ２ ０ １ ０． ０ ２． ０ ２ ２ ｊ
法的基础上提出了 磁 倾 角 欧 拉 反 褶 积 方 法 ， 它不需 要已知场 源 的 构 造 指 数 就 能 估 计 出 场 源 的 位 置 分 进一步可估计出场源的构造指数 ． 这将促进欧拉 布， 反褶积方法的发展和应用 ．
［５］ 概率 成 像 方 法 是 Ｐ 首次提出用于自然 ａ ｔ ｅ ｌ ｌ ａ１ ［６］ 电位异常的解释 ， 将这种方 法推 广到了 Ｍ ａ ｕ ｒ ｉ ｅ ｌ ｌ ｏ１ ［７］ 在任意 标 量 和 矢 量 地 球 物 理 场 重力领域 ． Ｉ ｕ ｌ ｉ ａ ｎ ｏ１
２ ） ） ， 狕 ＋犉·（ × （ 狔 狇 －狔 狇 －狕
］
（ ） １
其中 ， ０ 为真空磁导率 ， μ
２ ２ ２ （ ） ） ） ， 狉＝ 槡 狓 狕 ＋（ ＋（ 狔 狇 －狓 狇 －狔 狇 －狕 犾＝ ｃ ｏ ｓ 犐 ｃ ｏ ｓ 犃 ′ 犿 ＝ｃ ｏ ｓ 犐 ｓ ｉ ｎ 犃 ′ 狀 ＝ｓ ｉ ｎ 犐 ０ ０， ０ ０， ０， ｏ ｓ 犐 ｃ ｏ ｓ 犃 ′， 犕 ＝ｃ ｏ ｓ 犐 ｓ ｉ ｎ 犃 ′， 犖 ＝ｓ ｉ ｎ 犐 ′， 犔 ＝ ｃ
，犕 １犓 犲 犪 犫 狅 狉 犪 狋 狅 狉 犲 狅 犇 犲 狋 犲 犮 狋 犻 狅 狀（ 犆 犺 犻 狀 犪犝 狀 犻 狏 犲 狉 狊 犻 狋 犲 狅 狊 犮 犻 犲 狀 犮 犲 狊， 犅 犲 犻 犻 狀 犻 狀 犻 狊 狋 狉 犱 狌 犮 犪 狋 犻 狅 狀， 犅 犲 犻 犻 狀 ０ ０ ０ ８ ３， 犆 犺 犻 狀 犪 狔犔 狔狅 犳犌 狔狅 犳犌 犼 犵） 狔狅 犳犈 犼 犵１ ２犛 犮 犺 狅 狅 犾 狅 犲 狅 犺 狊 犻 犮 狊犪 狀 犱犐 狀 狅 狉 犿 犪 狋 犻 狅 狀犜 犲 犮 犺 狀 狅 犾 狅 犆 犺 犻 狀 犪犝 狀 犻 狏 犲 狉 狊 犻 狋 犲 狅 狊 犮 犻 犲 狀 犮 犲 狊， 犅 犲 犻 犻 狀 ０ ０ ０ ８ ３， 犆 犺 犻 狀 犪 犳犌 狆 狔 犳 犵 狔， 狔狅 犳犌 犼 犵１
基金项目 国家 自 然 科 学 基 金 青 年 科 学 基 金 项 目 （ 、 国家８ 和国家专项课题 ４ ０ ９ ０ Hale Waihona Puke Baidu ０ ３ ３） ６ ３计划课题（ ２ ０ ０ ６ ＡＡ ０ ６ Ｚ １ １ １， ２ ０ ０ ６ ＡＡ ０ ６ Ａ ２ ０ １ ３） （ ） 联合资助 ． Ｓ ｉ ｎ ｏ Ｐ ｒ ｏ ｂ ｅ ０ １ ０ ５ 作者简介 郭良辉 ， 男， 讲师 ， 博士生 ， 从事勘探地球物理学领域的教学和科研工作 ． ： ＿ １ ９ ８ ０ 年生 ， Ｅ ｍ ａ ｉ ｌ ｕ ｏ ｌ ｉ ａ ｎ ｈ ｕ ｉ ６ ３． ｃ ｏ ｍ ＠１ ｇ ｇ
ｏ ｆ ｅ ｕ ｉ ｖ ａ ｌ ｅ ｎ ｔｍ ａ ｎ ｅ ｔ ｉ ｓ ｍｏ ｆｇ ｅ ｏ ｌ ｏ ｉ ｃ ａ ｌａ ｎ ｏ ｍ ａ ｌ ｏ ｕ ｓｂ ｏ ｄ ｉ ｅ ｓｗ ｉ ｔ ｈｈ ｉ ｈｒ ｅ ｓ ｏ ｌ ｕ ｔ ｉ ｏ ｎｂ ｏ ｔ ｈｌ ｏ ｎ ｉ ｔ ｕ ｄ ｉ ｎ ａ ｌ ｌ ｑ ｇ ｇ ｇ ｇ ｙ ａ ｎ ｄｔ ｒ ａ ｎ ｓ ｖ ｅ ｒ ｓ ｅ ｌ ． ｙ ， ， ， 犓 犲 狑 狅 狉 犱 狊 ａ ｎ ｅ ｔ ｉ ｃａ ｎ ｏ ｍ ａ ｌ ３ Ｄｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｍ ａ ｉ ｎ Ｅ ｎ ｔ ｒ ｏ ｉ ｌ ｔ ｅ ｒ Ａ ｎ ｏ ｍ ａ ｌ ｅ ａ ｒ ａ ｔ ｉ ｏ ｎ Ｍ ｇ ｙ ｇ ｇ ｐ ｙｆ ｙｓ ｐ 狔 实测磁异 常 的 反 演 解 释 则 是 磁 力 勘 探 解 发挥作用 ．
［９ ］ 可表示为 常Δ 犜 狓， 狕）１ 狔， 狇（
对磁异常反演方法和算 法 作了 较
多的研究与应用 ． 由于这种反演方法在数学上呈病 态， 具有固有的多 解 性 ， 需 引 入 先 验 约 束 条 件， 计算 维数巨大 、 计算量庞大 ， 这些问题都制约着反演效果
［］ ［］ 近 年 来， 姚长 和实 际 应 用 ． Ｐ ｏ ｒ ｔ ｎ ｉ ａ ｕ ｉ ｎ ｅ８ 、 Ｌ ｉ９ 、 ｇ
磁异常 Δ 犜 三维相关成像
２ ２ ２ 郭良辉１， ， 孟小红１， ， 石 磊１， 中国地质大学 ，北京 ） ，北京 １ １ 地下信息探测技术与仪器教育部重点实验室 （ ０ ０ ０ ８ ３ 北京 ） 地球物理与信息技术学院 ，北京 １ ２ 中国地质大学 （ ０ ０ ０ ８ ３
摘 要 本文将重力和重力梯度数据三维相关成像方法推广到磁 力 勘 探 领 域 ， 推导并建立了磁异常 Δ 犜 的三维相 关成像方法 ， 同时提出了基于熵滤波分离异 常 的 三 维 相 关 成 像 算 法 来 提 高 成 像 分 辨 率 ． 组合模型磁异常 Δ 犜 数据 和实际磁测资料试验分析表明 ， 本文方法能成像出地 下 地 质 体 的 空 间 赋 存 状 态 和 等 效 磁 性 分 布 ， 具有良好的横向 和纵向分辨率 ． 关键词 磁异常 ， 三维相关成像 ， 熵滤波 ， 异常分离 ： ／ ． ｉ ｓ ｓ ｎ ． ０ ０ ０ １ ５ ７ ３ ３ ． ２ ０ １ ０ ． ０ ２ ． ０ ２ ２ Ｄ Ｏ Ｉ １ ０ ． ３ ９ ６ ９ ｊ 中图分类号 Ｐ ６ ３ １ 收稿日期 ２ ， ０ ０ ９ １ １ ０ ３ ２ ０ １ ０ ０ １ １ ８收修定稿
犃 犫 狊 狋 狉 犪 犮 狋 ｅｅ ｘ ａ ｎ ｄｔ ｈ ｅ３ Ｄｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｍ ａ ｉ ｎ ｅ ｔ ｈ ｏ ｄｏ ｆｇ ｒ ａ ｖ ｉ ｔ ｎ ｄｇ ｒ ａ ｖ ｉ ｔ ｒ ａ ｄ ｉ ｏ ｍ ｅ ｔ ｒ ａ ｔ ａ Ｗ ｐ ｇ ｇｍ ｙａ ｙｇ ｙｄ ｉ ｎ ｔ ｏｍ ａ ｎ ｅ ｔ ｉ ｃ ｅ ｘ ｌ ｏ ｒ ａ ｔ ｉ ｏ ｎ ．Ｗ ｅｄ ｅ ｒ ｉ ｖ ｅ ａ ｎ ｄｐ ｒ ｏ ｏ ｓ ｅ ｔ ｈ ｅ ３ Ｄｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｍ ａ ｉ ｎ ｅ ｔ ｈ ｏ ｄ ｆ ｏ ｒｍ ａ ｎ ｅ ｔ ｉ ｃ ｇ ｐ ｐ ｇ ｇｍ ｇ ， ａ ｎ ｏ ｍ ａ ｌ 犜ｄ ａ ｔ ａａ ｎ ｄａ ｌ ｓ ｏｐ ｒ ｏ ｏ ｓ ｅ ｔ ｈ ｅａ ｌ ｏ ｒ ｉ ｔ ｈ ｍｏ ｆ ｔ ｈ ｅ３ Ｄｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｍ ａ ｉ ｎ ａ ｓ ｅ ｄｏ ｎｅ ｎ ｔ ｒ ｏ ｙΔ ｐ ｇ ｇ ｇｂ ｐ ｙ ，ｗ ｆ ｉ ｌ ｔ ｅ ｒ ｆ ｏ ｒａ ｎ ｏ ｍ ａ ｌ ｅ ａ ｒ ａ ｔ ｉ ｏ ｎ ｈ ｉ ｃ ｈｃ ａ ｎｂ ｅｕ ｓ ｅ ｄｔ ｏｉ ｍ ｒ ｏ ｖ ｅｔ ｈ ｅｒ ｅ ｓ ｏ ｌ ｕ ｔ ｉ ｏ ｎｏ ｆ ｔ ｈ ｅｉ ｍ ａ ｉ ｎ ．Ｔ ｈ ｅ ｐ ｇ ｇ ｙｓ ｐ
１ 引 言
磁力 勘 探 是 所 有 地 球 物 理 方 法 中 发 展 最 早 、 应 用最广泛的一种方法 ， 在区域地质调查 、 矿产资源和 油气勘探 、 工程地质 和 考 古 等 方 面 都 可 不 同 程 度 地
常用的磁异常反 决各种实际地质问 题 的 重 要 环 节 ． 演方法有物性反演法和欧拉反褶积法 ． 磁异 常 物 性 反 演 方 法 是 基 于 反 演 理 论 、 在最小 二乘意义下使目标函数达到极小的线性或非线性反 演． 在引入足够约束条件的情况下 ， 能够给出接近于
犕狇 ０ ２ ） 犜 狓， 狕）＝ μ · ５ 犃·（ 狓 Δ 狔， 狇（ 狇 －狓 ４ π 狉 ２ ） ） ） ·（ 狓 狕 ＋犅·（ ＋犆·（ 狔 狇 －狔 狇 －狓 狇 －狕
［
） ） ） ·（ 狕 狓 ＋犈·（ ＋ 犇·（ 狔 狇 －狔 狇 －狕 狇 －狓
犜ｄ ａ ｔ ａｏ ｆｓ ｎ ｔ ｈ ｅ ｓ ｉ ｚ ｅ ｄｍ ｕ ｌ ｔ ｉ ｌ ｅｒ ｅ ｃ ｔ ａ ｎ ｕ ｌ ａ ｒｐ ｒ ｉ ｓ ｍ ｓｍ ｏ ｄ ｅ ｌａ ｎ ｄｔ ｈ ｅ ｔ ｅ ｓ ｔｏ ｎｔ ｈ ｅｍ ａ ｎ ｅ ｔ ｉ ｃａ ｎ ｏ ｍ ａ ｌ ｙ ｐ ｇ ｇ ｙΔ ， ｒ ａ ｗｍ ａ ｎ ｅ ｔ ｉ ｃｄ ａ ｔ ａ ｓ ｈ ｏ ｗｔ ｈ ａ ｔ ｔ ｈ ｅｍ ｅ ｔ ｈ ｏ ｄｃ ａ ｎ ｉ ｍ ａ ｅ ｔ ｈ ｅ ｓ ａ ｃ ｅ ｅ ｘ ｔ ｅ ｎ ｄ ｅ ｄｓ ｔ ａ ｔ ｅ ｓ ａ ｎ ｄ ｔ ｈ ｅｄ ｉ ｓ ｔ ｒ ｉ ｂ ｕ ｔ ｉ ｏ ｎ ｇ ｇ ｐ
３ 犇犮 狅 狉 狉 犲 犾 犪 狋 犻 狅 狀 犻 犿 犪 犻 狀 犳 狅 狉犿 犪 狀 犲 狋 犻 犮犪 狀 狅 犿 犪 犾 犜犱 犪 狋 犪 犵 犵 犵 狔Δ
１２ １２ １２ ，ＭＥ ＧＵＯ Ｌ ｉ ａ ｎ Ｈ ｕ ｉ ＮＧ Ｘ ｉ ａ ｏ Ｈ ｏ ｎ Ｓ Ｈ ＩＬ ｅ ｉ ｇ ｇ ， ， ， ，
利
、 Ｐ ｉ ｌ ｋ ｉ ｎ ｔ ｏ ｎ 等提出了一些改进方法和措 ｇ 在一定程度上改善了反演效果和效率 ． 施，
］ １ ３ 磁异常欧拉 反 褶 积 方 法 ［ 以欧拉齐次方程为
［ １ ０， １ １］
［ ］ １ ２
运用磁异常 、 空间导数以及各种地质体具有特 基础 ， 定的 “ 构造指数 ” 来 确 定 异 常 场 源 的 位 置． 它无需场 计算简便快速 ， 适合于大面积磁 源物性的先验知识 ， 但是实际场源构造指数的未知 测资料的分析解释 ． 性、 反褶积解的发散 性 等 实 际 问 题 制 约 着 欧 拉 反 褶
４ ３ ６
地 球 物 理 学 报（ ） Ｃ ｈ ｉ ｎ ｅ ｓ ｅＪ ． Ｇ ｅ ｏ ｈ ｓ ． ｐ ｙ
５ ３卷
实际地质情况的物性分布和几何形态 ． 几十年来 ， 国 内外众多学者
［ １～７］
在基准面上 ， 假设测区地下任意 狕 取垂直向下为正 ， 坐标为 （ 、 体积为 狏 狓 狕 狔 狇， 狇， 狇） 狇 的第狇 个均匀磁化 球 体的磁化强度为犑狇 ， 磁矩 犕狇 ＝ 犑 地磁场倾角为 狏 狇 狇， 偏角为 犃 磁化 强度 倾 角 为 犐、 偏角为 犃 由于 犐 ′ ′． ０、 ０， 均匀磁化球体外部磁场与位于球心的偶极子磁场等 因此 ， 这 里 把 均 匀 磁 化 球 体 视 为 偶 极 子． 则 第狇 效， 个偶极子在测区上任意点 （ 狓， 狕）处的 总强 度磁 异 狔，
·犔 －狀·犖 － 犿·犕 ， 犾 犃 ＝ ２· ·犔 －狀·犖 ， 犅 ＝ ２·犿·犕 －犾 ， ·犖 ＋狀·犔） 犾 犆 ＝ ３·（ ， 犿·犖 ＋狀·犕 ） 犇 ＝ ３·（ ， ·犕 ＋ 犿·犔） 犾 犈 ＝ ３·（ ·犔． 狀·犖 － 犿·犕 －犾 犉 ＝ ２· 定义测区上实测 Δ 犜 异常与第狇 个 偶极 子测 区 ） 互相关系数为 上的 Δ 犜 Ｐ ｅ ａ ｒ ｓ ｏ ｎ 狇 异常的皮氏 （
第５ ３卷 第２期 ２ ０ １ ０年２月
地 球 物 理 学 报
ＣＨ Ｉ Ｎ Ｅ Ｓ Ｅ Ｊ ＯＵＲ ＮＡ Ｌ Ｏ Ｆ Ｇ Ｅ Ｏ Ｐ ＨＹ Ｓ Ｉ Ｃ Ｓ
Ｖ ｏ ｌ ． ５ ３，Ｎ ｏ ． ２ ， Ｆ ｅ ｂ ． ２ ０ １ ０
郭良辉 ， 孟小红 ， 石 磊． 磁异常 Δ 地 球 物 理 学 报， ： ： ／ 犜 三 维 相 关 成 像． ２ ０ １ ０， ５ ３（ ２） ４ ３ ５～４ ４ １， Ｄ Ｏ Ｉ １ ０． ３ ９ ６ ９ ． ｉ ｓ ｓ ｎ ． ０ ０ ０ １ ５ ７ ３ ３． ｊ ２ ０ １ ０． ０ ２． ０ ２ ２ ） ， Ｇ ｕ ｏＬ Ｈ，Ｍ ｅ ｎ Ｓ ｈ ｉ Ｌ． ３ Ｄｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｍ ａ ｉ ｎ ｆ ｏ ｒｍ ａ ｎ ｅ ｔ ｉ ｃ ａ ｎ ｏ ｍ ａ ｌ 犜ｄ ａ ｔ ａ ． 犆 犺 犻 狀 犲 狊 犲 犑． 犌 犲 狅 犺 狊 ．（ ｉ ｎＣ ｈ ｉ ｎ ｅ ｓ ｅ ２ ０ １ ０， ｇＸ Ｈ， ｇ ｇ ｇ ｙΔ 狆 狔 （ ） ： ： ／ ５ ３ ２ ４ ３ ５～４ ４ １， Ｄ Ｏ Ｉ １ ０． ３ ９ ６ ９ ． ｉ ｓ ｓ ｎ ． ０ ０ ０ １ ５ ７ ３ ３． ２ ０ １ ０． ０ ２． ０ ２ ２ ｊ
法的基础上提出了 磁 倾 角 欧 拉 反 褶 积 方 法 ， 它不需 要已知场 源 的 构 造 指 数 就 能 估 计 出 场 源 的 位 置 分 进一步可估计出场源的构造指数 ． 这将促进欧拉 布， 反褶积方法的发展和应用 ．
［５］ 概率 成 像 方 法 是 Ｐ 首次提出用于自然 ａ ｔ ｅ ｌ ｌ ａ１ ［６］ 电位异常的解释 ， 将这种方 法推 广到了 Ｍ ａ ｕ ｒ ｉ ｅ ｌ ｌ ｏ１ ［７］ 在任意 标 量 和 矢 量 地 球 物 理 场 重力领域 ． Ｉ ｕ ｌ ｉ ａ ｎ ｏ１
２ ） ） ， 狕 ＋犉·（ × （ 狔 狇 －狔 狇 －狕
］
（ ） １
其中 ， ０ 为真空磁导率 ， μ
２ ２ ２ （ ） ） ） ， 狉＝ 槡 狓 狕 ＋（ ＋（ 狔 狇 －狓 狇 －狔 狇 －狕 犾＝ ｃ ｏ ｓ 犐 ｃ ｏ ｓ 犃 ′ 犿 ＝ｃ ｏ ｓ 犐 ｓ ｉ ｎ 犃 ′ 狀 ＝ｓ ｉ ｎ 犐 ０ ０， ０ ０， ０， ｏ ｓ 犐 ｃ ｏ ｓ 犃 ′， 犕 ＝ｃ ｏ ｓ 犐 ｓ ｉ ｎ 犃 ′， 犖 ＝ｓ ｉ ｎ 犐 ′， 犔 ＝ ｃ
，犕 １犓 犲 犪 犫 狅 狉 犪 狋 狅 狉 犲 狅 犇 犲 狋 犲 犮 狋 犻 狅 狀（ 犆 犺 犻 狀 犪犝 狀 犻 狏 犲 狉 狊 犻 狋 犲 狅 狊 犮 犻 犲 狀 犮 犲 狊， 犅 犲 犻 犻 狀 犻 狀 犻 狊 狋 狉 犱 狌 犮 犪 狋 犻 狅 狀， 犅 犲 犻 犻 狀 ０ ０ ０ ８ ３， 犆 犺 犻 狀 犪 狔犔 狔狅 犳犌 狔狅 犳犌 犼 犵） 狔狅 犳犈 犼 犵１ ２犛 犮 犺 狅 狅 犾 狅 犲 狅 犺 狊 犻 犮 狊犪 狀 犱犐 狀 狅 狉 犿 犪 狋 犻 狅 狀犜 犲 犮 犺 狀 狅 犾 狅 犆 犺 犻 狀 犪犝 狀 犻 狏 犲 狉 狊 犻 狋 犲 狅 狊 犮 犻 犲 狀 犮 犲 狊， 犅 犲 犻 犻 狀 ０ ０ ０ ８ ３， 犆 犺 犻 狀 犪 犳犌 狆 狔 犳 犵 狔， 狔狅 犳犌 犼 犵１
基金项目 国家 自 然 科 学 基 金 青 年 科 学 基 金 项 目 （ 、 国家８ 和国家专项课题 ４ ０ ９ ０ Hale Waihona Puke Baidu ０ ３ ３） ６ ３计划课题（ ２ ０ ０ ６ ＡＡ ０ ６ Ｚ １ １ １， ２ ０ ０ ６ ＡＡ ０ ６ Ａ ２ ０ １ ３） （ ） 联合资助 ． Ｓ ｉ ｎ ｏ Ｐ ｒ ｏ ｂ ｅ ０ １ ０ ５ 作者简介 郭良辉 ， 男， 讲师 ， 博士生 ， 从事勘探地球物理学领域的教学和科研工作 ． ： ＿ １ ９ ８ ０ 年生 ， Ｅ ｍ ａ ｉ ｌ ｕ ｏ ｌ ｉ ａ ｎ ｈ ｕ ｉ ６ ３． ｃ ｏ ｍ ＠１ ｇ ｇ
ｏ ｆ ｅ ｕ ｉ ｖ ａ ｌ ｅ ｎ ｔｍ ａ ｎ ｅ ｔ ｉ ｓ ｍｏ ｆｇ ｅ ｏ ｌ ｏ ｉ ｃ ａ ｌａ ｎ ｏ ｍ ａ ｌ ｏ ｕ ｓｂ ｏ ｄ ｉ ｅ ｓｗ ｉ ｔ ｈｈ ｉ ｈｒ ｅ ｓ ｏ ｌ ｕ ｔ ｉ ｏ ｎｂ ｏ ｔ ｈｌ ｏ ｎ ｉ ｔ ｕ ｄ ｉ ｎ ａ ｌ ｌ ｑ ｇ ｇ ｇ ｇ ｙ ａ ｎ ｄｔ ｒ ａ ｎ ｓ ｖ ｅ ｒ ｓ ｅ ｌ ． ｙ ， ， ， 犓 犲 狑 狅 狉 犱 狊 ａ ｎ ｅ ｔ ｉ ｃａ ｎ ｏ ｍ ａ ｌ ３ Ｄｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｍ ａ ｉ ｎ Ｅ ｎ ｔ ｒ ｏ ｉ ｌ ｔ ｅ ｒ Ａ ｎ ｏ ｍ ａ ｌ ｅ ａ ｒ ａ ｔ ｉ ｏ ｎ Ｍ ｇ ｙ ｇ ｇ ｐ ｙｆ ｙｓ ｐ 狔 实测磁异 常 的 反 演 解 释 则 是 磁 力 勘 探 解 发挥作用 ．
［９ ］ 可表示为 常Δ 犜 狓， 狕）１ 狔， 狇（
对磁异常反演方法和算 法 作了 较
多的研究与应用 ． 由于这种反演方法在数学上呈病 态， 具有固有的多 解 性 ， 需 引 入 先 验 约 束 条 件， 计算 维数巨大 、 计算量庞大 ， 这些问题都制约着反演效果
［］ ［］ 近 年 来， 姚长 和实 际 应 用 ． Ｐ ｏ ｒ ｔ ｎ ｉ ａ ｕ ｉ ｎ ｅ８ 、 Ｌ ｉ９ 、 ｇ
磁异常 Δ 犜 三维相关成像
２ ２ ２ 郭良辉１， ， 孟小红１， ， 石 磊１， 中国地质大学 ，北京 ） ，北京 １ １ 地下信息探测技术与仪器教育部重点实验室 （ ０ ０ ０ ８ ３ 北京 ） 地球物理与信息技术学院 ，北京 １ ２ 中国地质大学 （ ０ ０ ０ ８ ３
摘 要 本文将重力和重力梯度数据三维相关成像方法推广到磁 力 勘 探 领 域 ， 推导并建立了磁异常 Δ 犜 的三维相 关成像方法 ， 同时提出了基于熵滤波分离异 常 的 三 维 相 关 成 像 算 法 来 提 高 成 像 分 辨 率 ． 组合模型磁异常 Δ 犜 数据 和实际磁测资料试验分析表明 ， 本文方法能成像出地 下 地 质 体 的 空 间 赋 存 状 态 和 等 效 磁 性 分 布 ， 具有良好的横向 和纵向分辨率 ． 关键词 磁异常 ， 三维相关成像 ， 熵滤波 ， 异常分离 ： ／ ． ｉ ｓ ｓ ｎ ． ０ ０ ０ １ ５ ７ ３ ３ ． ２ ０ １ ０ ． ０ ２ ． ０ ２ ２ Ｄ Ｏ Ｉ １ ０ ． ３ ９ ６ ９ ｊ 中图分类号 Ｐ ６ ３ １ 收稿日期 ２ ， ０ ０ ９ １ １ ０ ３ ２ ０ １ ０ ０ １ １ ８收修定稿
犃 犫 狊 狋 狉 犪 犮 狋 ｅｅ ｘ ａ ｎ ｄｔ ｈ ｅ３ Ｄｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｍ ａ ｉ ｎ ｅ ｔ ｈ ｏ ｄｏ ｆｇ ｒ ａ ｖ ｉ ｔ ｎ ｄｇ ｒ ａ ｖ ｉ ｔ ｒ ａ ｄ ｉ ｏ ｍ ｅ ｔ ｒ ａ ｔ ａ Ｗ ｐ ｇ ｇｍ ｙａ ｙｇ ｙｄ ｉ ｎ ｔ ｏｍ ａ ｎ ｅ ｔ ｉ ｃ ｅ ｘ ｌ ｏ ｒ ａ ｔ ｉ ｏ ｎ ．Ｗ ｅｄ ｅ ｒ ｉ ｖ ｅ ａ ｎ ｄｐ ｒ ｏ ｏ ｓ ｅ ｔ ｈ ｅ ３ Ｄｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｍ ａ ｉ ｎ ｅ ｔ ｈ ｏ ｄ ｆ ｏ ｒｍ ａ ｎ ｅ ｔ ｉ ｃ ｇ ｐ ｐ ｇ ｇｍ ｇ ， ａ ｎ ｏ ｍ ａ ｌ 犜ｄ ａ ｔ ａａ ｎ ｄａ ｌ ｓ ｏｐ ｒ ｏ ｏ ｓ ｅ ｔ ｈ ｅａ ｌ ｏ ｒ ｉ ｔ ｈ ｍｏ ｆ ｔ ｈ ｅ３ Ｄｃ ｏ ｒ ｒ ｅ ｌ ａ ｔ ｉ ｏ ｎ ｉ ｍ ａ ｉ ｎ ａ ｓ ｅ ｄｏ ｎｅ ｎ ｔ ｒ ｏ ｙΔ ｐ ｇ ｇ ｇｂ ｐ ｙ ，ｗ ｆ ｉ ｌ ｔ ｅ ｒ ｆ ｏ ｒａ ｎ ｏ ｍ ａ ｌ ｅ ａ ｒ ａ ｔ ｉ ｏ ｎ ｈ ｉ ｃ ｈｃ ａ ｎｂ ｅｕ ｓ ｅ ｄｔ ｏｉ ｍ ｒ ｏ ｖ ｅｔ ｈ ｅｒ ｅ ｓ ｏ ｌ ｕ ｔ ｉ ｏ ｎｏ ｆ ｔ ｈ ｅｉ ｍ ａ ｉ ｎ ．Ｔ ｈ ｅ ｐ ｇ ｇ ｙｓ ｐ
１ 引 言
磁力 勘 探 是 所 有 地 球 物 理 方 法 中 发 展 最 早 、 应 用最广泛的一种方法 ， 在区域地质调查 、 矿产资源和 油气勘探 、 工程地质 和 考 古 等 方 面 都 可 不 同 程 度 地
常用的磁异常反 决各种实际地质问 题 的 重 要 环 节 ． 演方法有物性反演法和欧拉反褶积法 ． 磁异 常 物 性 反 演 方 法 是 基 于 反 演 理 论 、 在最小 二乘意义下使目标函数达到极小的线性或非线性反 演． 在引入足够约束条件的情况下 ， 能够给出接近于