井中三分量磁测的梯度张量欧拉反褶积及应用
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
［ ２ ， ７ － １ ０ ］
３ 井中三分量磁测的梯度张量反演— — —欧 拉反褶积方法
３ ． １ 方法原理 欧拉反褶积方法源于 １ ８世纪的欧拉齐次方程， １ ９ ６ ５年， Ｈ ｏ ｏ ｄ首先提出了针对磁法的欧拉齐次方 程， 并利用 ２ Ｄ欧拉方程和特征点法得到了场源深 度、 总场极大值、 垂直梯度极大值以及构造指数之间 的关系。 Ｓ ｌ ａ ｃ ｋ 、 Ｂ ａ ｒ ｏ ｎ ｇ ｏ 、 Ｔ ｈ ｏ ｍ ｐ ｓ ｏ 、 Ｒ ｅ ｉ ｄ等进一步发 展了欧拉反褶积方法， 该方法在地面高精度磁测资
２ ２１ ／ ２ · ｚ
；
２ ２１ ／ ２ · ｚ
； 。
Ｓ （ ｕ ， ｖ ， ｚ ） ＝ Ｂ ｚ ２ ｕ ＋ｖ ） Ｓ （ ｕ ， ｖ ， ０ ） ｅ π（ Ｂ
收稿日期： ２ ０ １ ２－ ０ ３－ ３ ０
笔者研究了磁力梯度张量在井中三分量磁测中 的应用。利用井中三分量磁测资料 Ｂ 、 Ｂ 、 Ｂ ｘ ｙ ｚ 换算 磁力梯度张量 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ ， ｘ ｘ ｘ ｙ ｘ ｚ ｙ ｘ ｙ ｙ ｙ ｚ ｚ ｘ ｚ ｙ ｚ ｚ 再根据欧拉反褶积方法的原理， 采用联合反演方法 求解场源体的空间位置。理论模型结果表明， 该方 法不仅具有方便快速、 且不必知道场源体的物性参 数与准确几何形态就能够求解场源体位置的优点， 而且磁力梯度张量比单分量反演结果精度更高。将 该方法用于湖北大冶铁矿 １ ８ ２井三分量磁测资料 的解释， 取得较好的地质效果。
２ 井中磁测三分量梯度张量的计算
目前还没有仪器可直接用来测量井中磁力梯度 、 Ｂ 、 Ｂ 张量。但可以由井中磁测三分量 Ｂ ｘ ｙ ｚ 换算磁 力梯度 张 量 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 。 ｘ ｘ ｘ ｙ ｘ ｚ ｙ ｘ ｙ ｙ ｙ ｚ ｚ ｘ ｚ ｙ ｚ ｚ 具体的换算采用频率域方法或空间域方法， 频率域
第 ３ ７卷第 ４期 ２ ０ １ ３年 ８月
物 探 与 化 探
Ｇ Ｅ Ｏ Ｐ Ｈ Ｙ Ｓ Ｉ Ｃ Ａ Ｌ＆Ｇ Ｅ Ｏ Ｃ Ｈ Ｅ Ｍ Ｉ Ｃ Ａ ＬＥ Ｘ Ｐ Ｌ Ｏ Ｒ Ａ Ｔ Ｉ Ｏ Ｎ
Ｖ ｏ ｌ ． ３ ７ ， Ｎ ｏ ． ４ Ａ ｕ ｇ ． ， ２ ０ １ ３
Ｄ Ｏ Ｉ ： １ ０ ． １ １ ７ ２ ０ ／ ｊ ． ｉ ｓ ｓ ｎ ． １ ０ ０ ０－ ８ ９ １ ８ ． ２ ０ １ ３ ． ４ ． １ ２
３－ ６ ］ １ 磁力梯度张量概念 ［
磁力梯度张量定义为
１ １ ｇ １ ２ ｇ １ ３ ｇ Ｇ＝ ｇ ｇ ｇ １ ２ ２ ２ ３ ＝ ２ ｇ ３ １ ｇ ３ ２ ｇ ３ ３ ２ Ｕ ２ ｘ ２ Ｕ ｙ ｘ ２ Ｕ ｚ ｘ ２ ２ Ｂ Ｂ Ｂ Ｕ Ｕ Ｘ Ｘ Ｘ ｘ ｙ ｘ ｚ ｘ ｙ ｚ ２ ２ Ｂ Ｂ Ｂ Ｕ Ｕ ＝ Ｙ Ｙ Ｙ （ １ ） ， ２ ｙ ｚ ｘ ｙ ｚ ｙ ２ ２ Ｂ Ｂ Ｂ Ｕ Ｕ Ｚ Ｚ Ｚ ２ ｚ ｙ ｘ ｙ ｚ ｚ
wk.baidu.com
欧拉齐次方程可写为 （ ｘ－ｘ ） ｆ ／ ｘ＋（ ｙ－ｙ ） ｆ ／ ｙ＋ ０ ０ （ ｚ －ｚ ） ｆ ／ ｚ＝－Ｎ ·ｆ ， ０ （ ２ ） 式中， ｆ 在地面磁测资料解释中为总磁场强度异常 Ｔ ，（ ｘ ， ｙ ， ｚ ） 为计算点的坐标， （ ｘ ， ｙ ， ｚ ） 为场源体 Δ ０ ０ ０ Ｎ为与场源体类型有关的参数， 称构造指 的位置， １ ， 偶极线源 Ｎ＝ ２ ， 偶极子源 Ｎ 数， 如单极线源 Ｎ＝ ＝ ３ 。该方法通过设置一个滑动窗口， 利用窗口内的 Ｔ观测值， 根据最小二乘法求解线性方程组得到对 Δ 场源水平位置、 深度的估计。通过移动滑动窗口与 改变滑动窗口的大小， 利用统计方法剔除不合理的 解， 最后得到场源的水平位置、 深度值。 若式（ ２ ） 中ｆ 为井中磁测的三个分量 Ｂ 、 Ｂ 、 ｘ ｙ Ｂ ， 根据联合反演的原理， 把单分量的欧拉反褶积方 ｚ 法推广为三分量的联合反演， 可以写出井中三分量 磁测联合反演的方程为 ｗ ［ （ ｘ－ｘ ） Ｂ ｙ－ｙ ） Ｂ ｚ －ｚ ） Ｂ ］＋ １ ０ ｘ ｘ ＋（ ０ ｘ ｙ ＋（ ０ ｘ ｚ ｗ ［ （ ｘ－ｘ ） Ｂ ｙ－ｙ ） Ｂ ｚ －ｚ ） Ｂ ］＋ ２ ０ ｙ ｘ ＋（ ０ ｙ ｙ ＋（ ０ ｙ ｚ ｗ ［ （ ｘ－ｘ ） Ｂ ｙ－ｙ ） Ｂ ｚ －ｚ ） Ｂ ］＝ ３ ０ ｚ ｘ ＋（ ０ ｚ ｙ ＋（ ０ ｚ ｚ －Ｎ ［ ｗ （ Ｂ ）＋ｗ （ Ｂ ）＋ｗ （ Ｂ ） ］ ， １ ｘ －ｂ ｘ ２ ｙ －ｂ ｙ ３ ｚ －ｂ ｚ ３ ） （ 式中， ｗ 、 ｗ 、 ｗ 、 Ｂ 、 Ｂ 可 １ ２ ３ 分别为 Ｂ ｘ ｙ ｚ分量的权系数， 以根据观测精度来定权系数进行三分量的联合反 演， 也可以只对某一分量进行欧拉反褶积。如当 ｗ １ ＝ １ 、 ｗ ０ 、 ｗ ０或 ｗ ０ 、 ｗ １ 、 ｗ ０或 ｗ ２＝ ３＝ １＝ ２＝ ３＝ １＝ ０ 、 ｗ ０ 、 ｗ １时， 式（ ３ ） 就是利用单一分量 Ｂ ２＝ ３＝ ｘ或 Ｂ 、 Ｂ 、 ｂ 、 ｂ ｙ ｚ求场源深度的欧拉反褶积方法。 ｂ ｘ ｙ ｚ为 背景场（ 或称区域场） 。场源点坐标 ｘ 、 ｙ 、 ｚ ０ ０ ０ 与背 、 ｂ 、 ｂ 个待求的未知参数， 而 ９个张量 景场 ｂ ｘ ｙ ｚ为 ６ 元素 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 与三分量 ｘ ｘ ｘ ｙ ｘ ｚ ｙ ｘ ｙ ｙ ｙ ｚ ｚ ｘ ｚ ｙ ｚ ｚ Ｂ 、 Ｂ 、 Ｂ 、 构造指数 Ｎ为已知值。 ｘ ｙ ｚ 由式（ ３ ） 可以看出： 欧拉方程的齐次项， 即三分 量Ｂ 、 Ｂ 、 Ｂ 、 ｙ 、 ｚ 的偏导数 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 ｘ ｙ ｚ对 ｘ ｘ ｘ ｘ ｙ ｘ ｚ ｙ ｘ ｙ ｙ
井中三分量磁测的梯度张量欧拉反褶积及应用
刘天佑１， 高文利２， 冯杰２， 习宇飞１， 欧洋１
（ １ ． 中国地质大学 地球物理与空间信息学院， 湖北 武汉 ４ ３ ０ ０ ７ ４ ； ２ ． 中国地质科学院 地球物理地 球化学勘查研究所， 河北 廊坊 ０ ６ ５ ０ ０ ０ ）
摘 要： ２ ０ ０ ６年 Ｔ ｈ ｅＬ ｅ ａ ｄ ｉ ｎ ｇＥ ｄ ｇ ｅ 第一期详细介绍了磁力梯度张量的理论、 仪器与野外试验结果， 磁力梯度张量技术 已成为磁力勘探新的热点之一。笔者介绍了磁力梯度张量的概念及优点； 利用频率域与空间域方法把井中三分量 磁测资料换算磁力梯度张量； 根据联合反演原理与欧拉反褶积方法， 推导了磁力梯度张量的联合反演方程， 该方程 通过权函数矩阵可以灵活对一个或多个分量反演， 比文献［ １ ４ ］ 方法更具普遍性。理论模型结果表明， 磁力梯度张 量反演方法对井底异常的定位准确， 对于 ３ Ｄ模型， 一口钻井的资料也能较好确定空间位置。将该方法用于湖北大 ８ ２井三分量磁测资料的解释， 得出磁力梯度张量的欧拉解集中在 １ ０ ０～ １ ８ ０ｍ与 ５ ０ ０～ ５ ５ ０ｍ两个深度， 冶铁矿 １ 与钻探结果十分吻合。该结果证实了地质上关于铁矿体分布具两个台阶的推论， 对大冶铁矿的深部找矿具有实际 意义。 关键词： 磁力梯度张量； 井中三分量磁测；联合反演方程； 欧拉反褶积； 大冶铁矿 中图分类号：Ｐ ６ ３ １ ． ８ 文献标识码：Ａ 文章编号： １ ０ ０ ０－ ８ ９ １ ８ （ ２ ０ １ ３ ） ０ ４－ ０ ６ ３ ３－ ０ ７
地磁场是具有方向和幅值的矢量场， 在三维空 间中， 可以用由 ９个（ ３× ３的矩阵）空间梯度 Ｂ 、 ｘ ｘ Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 组成的张量来表示。 ｘ ｙ ｘ ｚ ｙ ｘ ｙ ｙ ｙ ｚ ｚ ｘ ｚ ｙ ｚ ｚ 直接进行磁场三分量的垂直梯度和水平梯度测量， 能获得更多的反映场源特点与细节的信息， 对磁异 常的解释十分重要， 磁力梯度技术也日益得到人们 的重 视。２ ０ ０ ６年 美 国 勘 探 地 球 物 理 协 会 的 刊 物 《 Ｔ ｈ ｅＬ ｅ ａ ｄ ｉ ｎ ｇＥ ｄ ｇ ｅ 》 第一期为关于磁力梯度勘探的 专刊， 详细介绍了磁力梯度张量的理论、 仪器与野外 采集、 试验结果及优点， 标志着磁力梯度张量技术已 １ － ２ ］ 成为磁力勘探中一个新的热点 ［ 。在我国， 地面 磁力测量还处于标量测量阶段， 井中磁测虽然开展 了三分量测量， 但没有进一步对资料进行深入处理。 对于磁力梯度张量的研究， 无论是在仪器的研制还 是梯度张量的理论方面， 目前我国做的工作还很少。 近年来， 井中磁测仪器有了很大发展， 在国内， 重庆、 北京、 上海等地质仪器厂研制生产了各种型号 的 井 中 磁 测 仪 器， 如 Ｊ Ｓ Ｚ 、 Ｊ Ｃ Ｘ １ 、 Ｊ Ｃ Ｘ ２ 、Ｊ Ｃ Ｘ ３ 、 Ｇ Ｊ Ｃ Ｘ １ 、 Ｊ Ｇ Ｓ １ 、 Ｊ Ｃ Ｃ ３ １ 、 Ｊ Ｓ Ｚ ０ ５等， 仪器性能逐渐稳 Ｃ Ｘ ３型井中三分 定， 自动化程度和精度提高。如 Ｊ 量磁力仪， 垂直分量精度大约可达到 ３ ０ｎ Ｔ ， 水平分 量精度大约可达到 ５ ０～ ６ ０ｎ Ｔ 。这使得我们有可能 利用磁力梯度张量的理论对井中三分量磁测资料作 进一步的研究。
７ － ９ ， １ ３ － １ ４ ］ 料Δ Ｔ得到广泛的应用 ［ 。
： 张量分量的等值线图类似于垂直磁
化、 对场源体分辨率高、 易于解释， 能够减少反演的 多解性； 每个张量分量代表一个方向滤波器， 能够突 出场源体相应方向上的结构； 张量分量的联合处理 解释能够给出场源的磁化方向信息； 能够直接得到 “ 点源” 的方位与磁矩， 提高了线源的分辨率等。磁 力梯度张量的这些优点预示了它巨大的应用潜力， 得到人们的普遍关注。
基金项目： 国土资源部“ 十二五” 项目（ １ ２ １ ２ ０ １ １ １ ２ ０ １ ９ ５ ） 、 国土资源部公益性行业科研专项（ ２ ０ ０ ９ １ １ ０ １ ７ ０ ３ ）
·６ ３ ４ ·
物 探 与 化 探
８ ， １ １ － １ ２ ］ Ｖ Ｄ算法 ［ 。
３ ７卷
式中， Ｕ为磁性体的磁位， 磁力梯度张量即磁位对 ｘ 、 ｙ 、 ｚ 的二阶导数， 或井中磁测三分量 Ｂ 、 Ｂ 、 Ｂ ｘ ｙ ｚ对 ｘ 、 ｙ 、 ｚ 的一阶导数。从式（ １ ） 可以看出， 磁力梯度张 ｇ ， 并且， 由于在无源区域 Ｕ遵 量是对称的， 即ｇ ｉ ｊ＝ ｊ ｉ 守Ｌ ａ ｐ ｌ ａ ｃ ｅ ＇ ｓ 方程， 因此 Ｇ 的对角元素之和为零， 即 ｇ ｇ ｇ ０ 。从这两个性质中可以得知梯度张 １ １＋ ２ ２＋ ３ ３＝ 量 ９个元素中只有 ５个元素是独立。 与总场和梯度相比， 磁力梯度张量有很多的优 点， 如
７ ］ ： 方法如下 ［
（ １ ） 把实测的井中磁测三分量异常值 Ｂ （ ｘ ， ｙ ， ０ ） 通过富里叶变换得到频谱 Ｓ （ ｕ ， ｖ ， ｚ ） ； Ｂ （ ２ ） 将Ｂ （ ｘ ， ｙ ， ０ ） 的频谱 Ｓ （ ｕ ， ｖ ， ｚ ） 乘以频率响 Ｂ
２ ２ １ ２ ｉ ｕ ， ２ ｉ ｖ ， ２ （ ｕ ＋ ｖ ）／ ， 得到频谱 应函数 ２ π π π ２ （ ｕ＋ ｖ） π Ｓ （ ｕ ， ｖ ， ｚ ） ＝２ ｉ ｕ Ｓ （ ｕ ， ｖ ， ０ ） ｅ π Ｂ Ｂ ｘ ２ （ ｕ＋ ｖ） π Ｓ （ ｕ ， ｖ ， ｚ ） ＝２ ｉ ｖ Ｓ （ ｕ ， ｖ ， ０ ） ｅ π Ｂ Ｂ ｙ
３ 井中三分量磁测的梯度张量反演— — —欧 拉反褶积方法
３ ． １ 方法原理 欧拉反褶积方法源于 １ ８世纪的欧拉齐次方程， １ ９ ６ ５年， Ｈ ｏ ｏ ｄ首先提出了针对磁法的欧拉齐次方 程， 并利用 ２ Ｄ欧拉方程和特征点法得到了场源深 度、 总场极大值、 垂直梯度极大值以及构造指数之间 的关系。 Ｓ ｌ ａ ｃ ｋ 、 Ｂ ａ ｒ ｏ ｎ ｇ ｏ 、 Ｔ ｈ ｏ ｍ ｐ ｓ ｏ 、 Ｒ ｅ ｉ ｄ等进一步发 展了欧拉反褶积方法， 该方法在地面高精度磁测资
２ ２１ ／ ２ · ｚ
；
２ ２１ ／ ２ · ｚ
； 。
Ｓ （ ｕ ， ｖ ， ｚ ） ＝ Ｂ ｚ ２ ｕ ＋ｖ ） Ｓ （ ｕ ， ｖ ， ０ ） ｅ π（ Ｂ
收稿日期： ２ ０ １ ２－ ０ ３－ ３ ０
笔者研究了磁力梯度张量在井中三分量磁测中 的应用。利用井中三分量磁测资料 Ｂ 、 Ｂ 、 Ｂ ｘ ｙ ｚ 换算 磁力梯度张量 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ ， ｘ ｘ ｘ ｙ ｘ ｚ ｙ ｘ ｙ ｙ ｙ ｚ ｚ ｘ ｚ ｙ ｚ ｚ 再根据欧拉反褶积方法的原理， 采用联合反演方法 求解场源体的空间位置。理论模型结果表明， 该方 法不仅具有方便快速、 且不必知道场源体的物性参 数与准确几何形态就能够求解场源体位置的优点， 而且磁力梯度张量比单分量反演结果精度更高。将 该方法用于湖北大冶铁矿 １ ８ ２井三分量磁测资料 的解释， 取得较好的地质效果。
２ 井中磁测三分量梯度张量的计算
目前还没有仪器可直接用来测量井中磁力梯度 、 Ｂ 、 Ｂ 张量。但可以由井中磁测三分量 Ｂ ｘ ｙ ｚ 换算磁 力梯度 张 量 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 。 ｘ ｘ ｘ ｙ ｘ ｚ ｙ ｘ ｙ ｙ ｙ ｚ ｚ ｘ ｚ ｙ ｚ ｚ 具体的换算采用频率域方法或空间域方法， 频率域
第 ３ ７卷第 ４期 ２ ０ １ ３年 ８月
物 探 与 化 探
Ｇ Ｅ Ｏ Ｐ Ｈ Ｙ Ｓ Ｉ Ｃ Ａ Ｌ＆Ｇ Ｅ Ｏ Ｃ Ｈ Ｅ Ｍ Ｉ Ｃ Ａ ＬＥ Ｘ Ｐ Ｌ Ｏ Ｒ Ａ Ｔ Ｉ Ｏ Ｎ
Ｖ ｏ ｌ ． ３ ７ ， Ｎ ｏ ． ４ Ａ ｕ ｇ ． ， ２ ０ １ ３
Ｄ Ｏ Ｉ ： １ ０ ． １ １ ７ ２ ０ ／ ｊ ． ｉ ｓ ｓ ｎ ． １ ０ ０ ０－ ８ ９ １ ８ ． ２ ０ １ ３ ． ４ ． １ ２
３－ ６ ］ １ 磁力梯度张量概念 ［
磁力梯度张量定义为
１ １ ｇ １ ２ ｇ １ ３ ｇ Ｇ＝ ｇ ｇ ｇ １ ２ ２ ２ ３ ＝ ２ ｇ ３ １ ｇ ３ ２ ｇ ３ ３ ２ Ｕ ２ ｘ ２ Ｕ ｙ ｘ ２ Ｕ ｚ ｘ ２ ２ Ｂ Ｂ Ｂ Ｕ Ｕ Ｘ Ｘ Ｘ ｘ ｙ ｘ ｚ ｘ ｙ ｚ ２ ２ Ｂ Ｂ Ｂ Ｕ Ｕ ＝ Ｙ Ｙ Ｙ （ １ ） ， ２ ｙ ｚ ｘ ｙ ｚ ｙ ２ ２ Ｂ Ｂ Ｂ Ｕ Ｕ Ｚ Ｚ Ｚ ２ ｚ ｙ ｘ ｙ ｚ ｚ
wk.baidu.com
欧拉齐次方程可写为 （ ｘ－ｘ ） ｆ ／ ｘ＋（ ｙ－ｙ ） ｆ ／ ｙ＋ ０ ０ （ ｚ －ｚ ） ｆ ／ ｚ＝－Ｎ ·ｆ ， ０ （ ２ ） 式中， ｆ 在地面磁测资料解释中为总磁场强度异常 Ｔ ，（ ｘ ， ｙ ， ｚ ） 为计算点的坐标， （ ｘ ， ｙ ， ｚ ） 为场源体 Δ ０ ０ ０ Ｎ为与场源体类型有关的参数， 称构造指 的位置， １ ， 偶极线源 Ｎ＝ ２ ， 偶极子源 Ｎ 数， 如单极线源 Ｎ＝ ＝ ３ 。该方法通过设置一个滑动窗口， 利用窗口内的 Ｔ观测值， 根据最小二乘法求解线性方程组得到对 Δ 场源水平位置、 深度的估计。通过移动滑动窗口与 改变滑动窗口的大小， 利用统计方法剔除不合理的 解， 最后得到场源的水平位置、 深度值。 若式（ ２ ） 中ｆ 为井中磁测的三个分量 Ｂ 、 Ｂ 、 ｘ ｙ Ｂ ， 根据联合反演的原理， 把单分量的欧拉反褶积方 ｚ 法推广为三分量的联合反演， 可以写出井中三分量 磁测联合反演的方程为 ｗ ［ （ ｘ－ｘ ） Ｂ ｙ－ｙ ） Ｂ ｚ －ｚ ） Ｂ ］＋ １ ０ ｘ ｘ ＋（ ０ ｘ ｙ ＋（ ０ ｘ ｚ ｗ ［ （ ｘ－ｘ ） Ｂ ｙ－ｙ ） Ｂ ｚ －ｚ ） Ｂ ］＋ ２ ０ ｙ ｘ ＋（ ０ ｙ ｙ ＋（ ０ ｙ ｚ ｗ ［ （ ｘ－ｘ ） Ｂ ｙ－ｙ ） Ｂ ｚ －ｚ ） Ｂ ］＝ ３ ０ ｚ ｘ ＋（ ０ ｚ ｙ ＋（ ０ ｚ ｚ －Ｎ ［ ｗ （ Ｂ ）＋ｗ （ Ｂ ）＋ｗ （ Ｂ ） ］ ， １ ｘ －ｂ ｘ ２ ｙ －ｂ ｙ ３ ｚ －ｂ ｚ ３ ） （ 式中， ｗ 、 ｗ 、 ｗ 、 Ｂ 、 Ｂ 可 １ ２ ３ 分别为 Ｂ ｘ ｙ ｚ分量的权系数， 以根据观测精度来定权系数进行三分量的联合反 演， 也可以只对某一分量进行欧拉反褶积。如当 ｗ １ ＝ １ 、 ｗ ０ 、 ｗ ０或 ｗ ０ 、 ｗ １ 、 ｗ ０或 ｗ ２＝ ３＝ １＝ ２＝ ３＝ １＝ ０ 、 ｗ ０ 、 ｗ １时， 式（ ３ ） 就是利用单一分量 Ｂ ２＝ ３＝ ｘ或 Ｂ 、 Ｂ 、 ｂ 、 ｂ ｙ ｚ求场源深度的欧拉反褶积方法。 ｂ ｘ ｙ ｚ为 背景场（ 或称区域场） 。场源点坐标 ｘ 、 ｙ 、 ｚ ０ ０ ０ 与背 、 ｂ 、 ｂ 个待求的未知参数， 而 ９个张量 景场 ｂ ｘ ｙ ｚ为 ６ 元素 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 与三分量 ｘ ｘ ｘ ｙ ｘ ｚ ｙ ｘ ｙ ｙ ｙ ｚ ｚ ｘ ｚ ｙ ｚ ｚ Ｂ 、 Ｂ 、 Ｂ 、 构造指数 Ｎ为已知值。 ｘ ｙ ｚ 由式（ ３ ） 可以看出： 欧拉方程的齐次项， 即三分 量Ｂ 、 Ｂ 、 Ｂ 、 ｙ 、 ｚ 的偏导数 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 ｘ ｙ ｚ对 ｘ ｘ ｘ ｘ ｙ ｘ ｚ ｙ ｘ ｙ ｙ
井中三分量磁测的梯度张量欧拉反褶积及应用
刘天佑１， 高文利２， 冯杰２， 习宇飞１， 欧洋１
（ １ ． 中国地质大学 地球物理与空间信息学院， 湖北 武汉 ４ ３ ０ ０ ７ ４ ； ２ ． 中国地质科学院 地球物理地 球化学勘查研究所， 河北 廊坊 ０ ６ ５ ０ ０ ０ ）
摘 要： ２ ０ ０ ６年 Ｔ ｈ ｅＬ ｅ ａ ｄ ｉ ｎ ｇＥ ｄ ｇ ｅ 第一期详细介绍了磁力梯度张量的理论、 仪器与野外试验结果， 磁力梯度张量技术 已成为磁力勘探新的热点之一。笔者介绍了磁力梯度张量的概念及优点； 利用频率域与空间域方法把井中三分量 磁测资料换算磁力梯度张量； 根据联合反演原理与欧拉反褶积方法， 推导了磁力梯度张量的联合反演方程， 该方程 通过权函数矩阵可以灵活对一个或多个分量反演， 比文献［ １ ４ ］ 方法更具普遍性。理论模型结果表明， 磁力梯度张 量反演方法对井底异常的定位准确， 对于 ３ Ｄ模型， 一口钻井的资料也能较好确定空间位置。将该方法用于湖北大 ８ ２井三分量磁测资料的解释， 得出磁力梯度张量的欧拉解集中在 １ ０ ０～ １ ８ ０ｍ与 ５ ０ ０～ ５ ５ ０ｍ两个深度， 冶铁矿 １ 与钻探结果十分吻合。该结果证实了地质上关于铁矿体分布具两个台阶的推论， 对大冶铁矿的深部找矿具有实际 意义。 关键词： 磁力梯度张量； 井中三分量磁测；联合反演方程； 欧拉反褶积； 大冶铁矿 中图分类号：Ｐ ６ ３ １ ． ８ 文献标识码：Ａ 文章编号： １ ０ ０ ０－ ８ ９ １ ８ （ ２ ０ １ ３ ） ０ ４－ ０ ６ ３ ３－ ０ ７
地磁场是具有方向和幅值的矢量场， 在三维空 间中， 可以用由 ９个（ ３× ３的矩阵）空间梯度 Ｂ 、 ｘ ｘ Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 、 Ｂ 组成的张量来表示。 ｘ ｙ ｘ ｚ ｙ ｘ ｙ ｙ ｙ ｚ ｚ ｘ ｚ ｙ ｚ ｚ 直接进行磁场三分量的垂直梯度和水平梯度测量， 能获得更多的反映场源特点与细节的信息， 对磁异 常的解释十分重要， 磁力梯度技术也日益得到人们 的重 视。２ ０ ０ ６年 美 国 勘 探 地 球 物 理 协 会 的 刊 物 《 Ｔ ｈ ｅＬ ｅ ａ ｄ ｉ ｎ ｇＥ ｄ ｇ ｅ 》 第一期为关于磁力梯度勘探的 专刊， 详细介绍了磁力梯度张量的理论、 仪器与野外 采集、 试验结果及优点， 标志着磁力梯度张量技术已 １ － ２ ］ 成为磁力勘探中一个新的热点 ［ 。在我国， 地面 磁力测量还处于标量测量阶段， 井中磁测虽然开展 了三分量测量， 但没有进一步对资料进行深入处理。 对于磁力梯度张量的研究， 无论是在仪器的研制还 是梯度张量的理论方面， 目前我国做的工作还很少。 近年来， 井中磁测仪器有了很大发展， 在国内， 重庆、 北京、 上海等地质仪器厂研制生产了各种型号 的 井 中 磁 测 仪 器， 如 Ｊ Ｓ Ｚ 、 Ｊ Ｃ Ｘ １ 、 Ｊ Ｃ Ｘ ２ 、Ｊ Ｃ Ｘ ３ 、 Ｇ Ｊ Ｃ Ｘ １ 、 Ｊ Ｇ Ｓ １ 、 Ｊ Ｃ Ｃ ３ １ 、 Ｊ Ｓ Ｚ ０ ５等， 仪器性能逐渐稳 Ｃ Ｘ ３型井中三分 定， 自动化程度和精度提高。如 Ｊ 量磁力仪， 垂直分量精度大约可达到 ３ ０ｎ Ｔ ， 水平分 量精度大约可达到 ５ ０～ ６ ０ｎ Ｔ 。这使得我们有可能 利用磁力梯度张量的理论对井中三分量磁测资料作 进一步的研究。
７ － ９ ， １ ３ － １ ４ ］ 料Δ Ｔ得到广泛的应用 ［ 。
： 张量分量的等值线图类似于垂直磁
化、 对场源体分辨率高、 易于解释， 能够减少反演的 多解性； 每个张量分量代表一个方向滤波器， 能够突 出场源体相应方向上的结构； 张量分量的联合处理 解释能够给出场源的磁化方向信息； 能够直接得到 “ 点源” 的方位与磁矩， 提高了线源的分辨率等。磁 力梯度张量的这些优点预示了它巨大的应用潜力， 得到人们的普遍关注。
基金项目： 国土资源部“ 十二五” 项目（ １ ２ １ ２ ０ １ １ １ ２ ０ １ ９ ５ ） 、 国土资源部公益性行业科研专项（ ２ ０ ０ ９ １ １ ０ １ ７ ０ ３ ）
·６ ３ ４ ·
物 探 与 化 探
８ ， １ １ － １ ２ ］ Ｖ Ｄ算法 ［ 。
３ ７卷
式中， Ｕ为磁性体的磁位， 磁力梯度张量即磁位对 ｘ 、 ｙ 、 ｚ 的二阶导数， 或井中磁测三分量 Ｂ 、 Ｂ 、 Ｂ ｘ ｙ ｚ对 ｘ 、 ｙ 、 ｚ 的一阶导数。从式（ １ ） 可以看出， 磁力梯度张 ｇ ， 并且， 由于在无源区域 Ｕ遵 量是对称的， 即ｇ ｉ ｊ＝ ｊ ｉ 守Ｌ ａ ｐ ｌ ａ ｃ ｅ ＇ ｓ 方程， 因此 Ｇ 的对角元素之和为零， 即 ｇ ｇ ｇ ０ 。从这两个性质中可以得知梯度张 １ １＋ ２ ２＋ ３ ３＝ 量 ９个元素中只有 ５个元素是独立。 与总场和梯度相比， 磁力梯度张量有很多的优 点， 如
７ ］ ： 方法如下 ［
（ １ ） 把实测的井中磁测三分量异常值 Ｂ （ ｘ ， ｙ ， ０ ） 通过富里叶变换得到频谱 Ｓ （ ｕ ， ｖ ， ｚ ） ； Ｂ （ ２ ） 将Ｂ （ ｘ ， ｙ ， ０ ） 的频谱 Ｓ （ ｕ ， ｖ ， ｚ ） 乘以频率响 Ｂ
２ ２ １ ２ ｉ ｕ ， ２ ｉ ｖ ， ２ （ ｕ ＋ ｖ ）／ ， 得到频谱 应函数 ２ π π π ２ （ ｕ＋ ｖ） π Ｓ （ ｕ ， ｖ ， ｚ ） ＝２ ｉ ｕ Ｓ （ ｕ ， ｖ ， ０ ） ｅ π Ｂ Ｂ ｘ ２ （ ｕ＋ ｖ） π Ｓ （ ｕ ， ｖ ， ｚ ） ＝２ ｉ ｖ Ｓ （ ｕ ， ｖ ， ０ ） ｅ π Ｂ Ｂ ｙ