重磁自动反演主要约束技术的分析实验.pdfVIP

·700- 中国地球物理2011 (19)地质调查和矿产勘查地球物理 重磁自动反演主要约束技术的分析实验 Researchandtestsofsomeconstraintson and inversion gravitymagnetic 周岐山木姚长利 郑元满 Zhou Yuanman Qishan’YaoChangliZheng 中国地质大学(北京)地下信息探测技术与仪器教育部重点实验室北京100083 1．引言 自动反演方法是当前勘探业的研究热点。随着现在勘探难度越来越大，测量数据精度越来越高，重、 磁勘探对反演工作提出了更高的技术要求。针对这种现实需要，自动反演方法应运而生。然而实际工作中 发现，自动反演方法具有多解性，直接求解的结果往往与实际情况不符，于是通过数学方法对自动反演施 加约束成为自动反演工作中最重要的一个环节。本文研究了几种比较重要的约束反演方法，分析了这些方 法所各自适用的条件，为自动反演的实用化提供了依据。 2．方法技术 自动反演方法的基本思想是将一定范围的地下空间剖分成许多规则单元，这些规则单元具有简单的几 何形态和均匀的密度或者磁性。反演计算的目的是计算每一个剖分单元的物性，所以这种方法也称为物性 反演方法。为了便于计算，一般都采用规则剖分的方式将地下空间剖分为长方体的组合。由位场理论得知， 密度或磁化率与模型的几何参数(磁性反演中还有磁化方向参数)之间是线性关系。则可推导线性反演目 标函数如下： 伊=(d—Gin)7’(d—Gm)+A缈 式中，d为观测数据向量，G为正演核函数矩阵，m为反演所求的剖分单元物性参数向量。目标函数分为 两项，前一项表示求取最小方差解，后一项是通过正则化方法加入的约束条件。 一个比较合理的反演计算结果需要多种约束条件共同作用，本文研究了下面一系列约束方法。 (1)深度约束和物性范围约束：深度约束是约束反演中最重要的一项约束条件。无约束反演会出现趋 肤效应，这是由于不同深度剖分单元在核函数中权重不同所导致的。深度约束通过对核函数矩阵加权克服 了这一现象。物性范围约束可以通过迭代计算中对超出限制的反演结果强制赋值实现。可以去除一些明显 不合理的结果。在自动反演中，这两条约束条件一般是不可缺少的。 (2)光滑约束、最小模型约束和轴向约束：这三种约束都是针对不同的场源形态提出的，要根据实际情 况斟酌使用。光滑约束通过约束模型梯度最小求最光滑模型，这种约束方法能较好的约束块状场源的形态， 但是对深度不同的多个场源的分辨能力很弱；最小模型约束通过使非零模型的数量最少求最小体积模型，在 剖分较粗糙、场源形态很简单的情况下，该约束方法基本能够完全约束出场源的形状和物性，但是当剖分单 元远远小于实际场源时，该约束的反演结果会使物性分布集中到一个单元内，所以该方法具有较强的局限性； 轴向约束通过约束模型转动惯量最小求沿某些己知的轴线方向展布的模型，对于轴线型的场源有出色的分辨 能力，不足之处在于无法反演块状的场源。这三种约束方法都能够通过正则化方法实现，并且研究证明光滑 约束可以结合轴线约束做多约束反演，这种多约束反演综合了两者的优点，是很实用的反演方法。 3．结论 针对不同情况，灵活运用不同的约束反演方法，毹够较好的克制自动反演中出现的多解性，达到较理 想的反演效果。 参考文献 [I】王家映．地球物理反演理论[M】．北京：高等教育出版社．1998． D．W．3-Dinversion l：394--408． 【2】Li．Y．，Oldenburg ofmagnetic 996，6 eata册．Geophysics,l