局部波数在磁异常解释应用中的方法技术.pdfVIP

　第３５卷第６期 物　探　与　化　探 Ｖｏｌ．３５，Ｎｏ．６　 　２０１１年１２月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ Ｄｅｃ．，２０１１　 局部波数在磁异常解释应用中的方法技术 崔莉１，２，３，４，王万银２，３，４ （１．惠州工业科技学校，广东 惠州　５１６００１；２．长安大学 重磁方法技术研究所，陕西 西安　 ７１００５４；３．长安大学地质工程与测绘学院，陕西 西安　７１００５４；４．长安大学 西部矿产资源与地质 工程教育部重点实验室，陕西 西安　７１００５４） 摘 要：研究了提高频率域中导数计算结果的精度以及局部波数反演结果的稳定性和精度，对比了余弦扩边方法和 最小曲率扩边方法在频率域中计算导数的精度。通过对比表明，最小曲率扩边结果的精度高于余弦扩边结果的精 度，提出了向上延拓某一高度来进行局部波数的反演技术。此技术不但能够解决高频干扰，提高导数计算精度，而 且利用不同延拓高度的反演结果能够得到一个平稳的反演结果，从而提高了局部波数反演的精度。理论模型试算 验证了方法技术的精度，实际资料处理验证了方法技术的稳定性，达到了满意的效果。 关键词：磁异常解释；局部波数；最小曲率；余弦扩边；向上延拓；导数计算 中图分类号：Ｐ６３１．２　　　文献标识码：Ａ　　　文章编号：１０００－８９１８（２０１１）０６－０７７９－０６ 　　磁异常边缘参数反演主要有线性特征法、解析 是以向下无限延深薄板的磁异常理论为基础来求地 信号法、欧拉反褶积法、Ｗｅｒｎｅｒ反褶积法以及局部 质体接触面的水平位置、上顶埋深、磁化强度的大小 ［８］ 波数法。线性特征法是利用具有明显走向和倾向特 及方向。该方法最早是由Ｓ．Ｗｅｒｎｅｒ 提出，当时该 征的抛物面来拟合磁性体对应的磁异常水平梯度的 方法主要用于分析瑞典的矿化岩墙。Ｈａｒｔｍａｎ等［９］ 异常，并求出所有网格节点上的特征值，由最大曲率 发表了应用这种方法的结果，他在该文章中通过模 特征值所分布的平面位置来确定磁性体的边界，再 型试算证明了Ｗｅｒｎｅｒ反褶积方法可以快速的反演 利用总梯度模及其二阶导数项，就可以求出磁性体 磁性体的埋深。 ［１］ ［２］ 顶面的埋深 。解析信号法是 Ｒｏｅｓｔ等人 提出。 ２０世纪９０年代局部波数法得到了广泛的关 解析信号振幅在磁性差异处出现极大值，根据这些 注，近十年来该方法得到了快速的发展。局部波数 极大值可以确定磁源边界位置。然而，除了薄板、垂 法适用于特定的场源模型，即倾斜的断层、倾斜的薄 直接触带等特殊情况外，在一般情况下，磁异常总梯 板和水平圆柱体，利用二维磁测数据的导出参数，可 度模极大值位置与磁源边界位置都存在一定偏差， 以反演异常体的特征点和埋深。该方法利用了局部 ［３］ 所以，黄临平和管志宁 采用理论推导与理论模型 波数的振幅通常不依赖于磁化强度和磁倾角并且不 ［１０］ 相结合的方法，分析了磁异常总梯度模极大值确定 需要先验信息就可求出磁性体的埋深的特点 。 磁源边界产生偏差的原因，并提出了减小误差的办 因此在确定磁性体位置上比磁场值更直接，并在磁 法。欧拉反褶积法是一种快速的反演方法，它不需 测资料解释中快捷简单，对于前期的勘探工作具有 要对地质模型作任何的假设，也无需对磁测数据进 重要意义。 行化极处理，因此其适用