基于Parker算法的磁性双界面模式正反演研究.pdfVIP

维普资讯 学兔兔 石油天然气学报 (江汉石油学院学报) 2008年6月 第3o卷 第3期 · 78 · Journal of Oil and Gas Technology(J．j PO Jun．2008 Vo1．30 No．3 基于Parker算法的磁性双界面模式正反演研究 相 鹏，刘 展 (={]国石油大学(华东)地球资源与信息学院，山东东营257061) [摘要]Parker界面反演算法中包含了一个不稳定因子——向下延拓算子，导致其不稳定，通常的做法是 在反演过程中加入一个低通滤波器以限制高频振荡。这种做法在消除高频假波的同时会对有效信号进行 一 定的压制，造成反演精度的降低。重新推导了正反演公式，提出了双磁性界面模式正反演公式和解决 反演发散问题的方法。双界面模式正反演公式在实际应用中更具合理性和普适性，解决反演发散问题的 方法保证了反演过程能够稳定收敛。讨论了算法实现的关键技术和措施，通过模型实验分析算法的精度 和可用性，可以发现算法计算精度很高，误差为0．8 ～4％，且具有一定抗噪性 用该方法处理中国某 地磁测数据，得到较好效果。 [关键词]双界面模式；界面反演；磁法反演；Parker法 [中图分类号]P631．325 [文献标识码]A [文章编号]1000—9752(2008)03—0078—05 随着区域地球物理工作的开展，研究区域及深部地质构造、寻找开发深部矿产资源已E{益为地质、 地球物理学家所重视。磁法勘探的应用领域也在不断扩大。确定磁性基底深度以及居里等温面深度对研 究区域及深部地质构造具有十分重要的意义。 很多学者提出了各种根据已知位场异常计算界面起伏的方法，包括利用等效棱柱体拟合界面进行反 演。该方法先计算每个棱柱体产生的异常，然后求和就能得到相应界面产生的异常，但该方法由于计算 速度慢而没有广泛采用。1973年，Parker在位场界面正演计算中引入快速傅里叶变换l】]，其方法具有 速度快、功能强的特点，成为重磁异常定量解释中的工具，适合于含油气构造区的地球物理资料综合解 释及联合反演。其后0ldenburg在此基础上，提出了确定起伏密度界面的迭代反演法 ]。尽管如此，还 存在3个问题有待讨论：一是所有的磁性体不应是单一界面模型，而应是由上、下界面所围的双界面模 式。尽管双界面模型可由两个单界面模型组合而成，但直接用双界面模型讨论更为合理，也具有普遍 性；二是如何加快正演时级数求和的收敛，Parker对单界面模型正演时的收敛速度作了讨论，但还不 完善，对于其在双界面模型中的应用需重新讨论；三是如何保证反演问题的稳定收敛，这是目前所面临 的最大困难。反演时指数因子的高频放大影响会使计算结果发散剧烈，即使当迭代剩余量很小时也会如 此，以往解决这一问题的唯一办法是添加低通滤波器，但此法并不能保证反演问题稳定收敛，并且低通 滤波器的参数选择难以控制。笔者在前人工作的基础上 ]，重新推导了正、反演公式，并进行了改 进，有效地解决了以上3个问题。 1 磁性界面正演原理 1．1 正演问题基本原理 正确地建立磁性界面的模型是研究界面反演方法的基础，由于界面模型在频率域可以实现快速计 算，故首先研究频率域中建立界面模型问题 (图1)。 设直角坐标系OXYZ的原点位于观测平面卜，Z轴垂直向下，观测点坐标为P(x， ，z)，磁性层点的坐 标为Q( ，叩， )。观测点P的矢量 为 一．r—i+ + ， 在观测平面 的投影为 一 + ；Q点的坐 [收稿日期]2008一o3一o2 [基金项目]国家 “863”计划项目 (2005AA06z204)。 [作者简介]相鹏 (1