用脊回法反演重力异常的多层密度及其界面-地震工程学报.PDF

第$@ 卷 第+ 期 西 北 地 震 学 报 O/8 $@ B/8 + $##+ $##+ 年$ 月 BCD9EFGH9GDB HGIHJCKCLIM(K .CNDB(K P4Q8 ， 用脊回法反演重力异常的多层密度及其界面 杨长福 （浙江大学地球科学系，浙江杭州 !##$% ） 摘 要：将研究区域划分成具有固定宽度的矩形网格，以网格密度和厚度作为模型参数，在此基 础上形成重力异常的反演目标函数，计算出对模型参数的偏导数矩阵，然后采用脊回归法对重力 异常进行反演而同时得到密度及其界面。以此方法对理论模型进行了反演试验。 关键词：脊回归法；重力异常；反演 中图分类号：’!$ 文献标识码：( 文章编号：### ) #*++（$##+ ）#+ ) #$,! ) #- # 引 言 反演总会因实际资料的观测误差及其数量的有限，而表现为非唯一性。反演的目标就是寻找一些模 型既能解释观测资料又能满足物理条件的限制。人们从最小二乘法发展了很多线性反演方法来解决这类 问题，但这些方法的迭代过程中都需要好的初始模型，以便迭代收敛。但即便如此，也未必保证能收敛于 真实模型［］ ［$ ］ 。脊回归法 就常被用来解决这类地球物理反演问题。这种方法的反演效果很大程度上取 决于所采用的阻尼过程，通常是在开始时使阻尼因子取很大的正值，使收敛过程趋近于最速下降法；然后 在每次迭代后乘以小于 的因子，使迭代过程接近最小二乘法。然而，./012342 不用最小特征值，并在迭 代发散时，用最大特征值来代替它［! ］。脊回归法仍带有上述两种方法中调整阻尼因子的特征，但采用了 ［+ ］ 两步最小化过程，使得效率提高 。 ［- ］ 密度界面的反演多局限于单层界面的反演 ，而多层密度 界面的反演比较困难。目前反演多层密度界面的方法主要有两 种：一是把多层密度界面的重力场，通过分离场的方法逐层分 离，逐层反演［@ ］ ［% ］ ；二是直接反演，如王一新的正则化反演法 和 ［* ］ 陈胜早的! 层变密度反演法 。这些方法都很复杂，且要求密 度分布已知，使用起来很不方便。目前同时进行多层密度及其 界面反演的还不多。本文将以单元密度和网格厚度为模型参 数，进行重力异常的多层密度及其界面的同时反演。 正演模型及计算方法 如图 所示，将二维地球介质划分成 个密度均匀分布的矩 形网格。为简便起见，网格划成宽度不等的 列，各列含有数目 # 相同，厚度不等，总数为$ 个的矩形网格，按行从左至右从上到 图 二维密度模型 下依次排序。图 中第% 个网格’() 四顶点的坐标分别为（#% ， 5678 904 :;/ ) 6=4236/21 4236:? =/48 $ ），（# ，$ ），（# ，$ ），（# ，$ ），其剩余密度为 ，则在坐标原点* 处（当计算点不在原点，可将坐标系 % %$ %$ %! %! %+ %+ !% 沿水平方向平移到该计算点，计算公式不变）的重力异常可由下式计算： 收稿日期：$##+A#+A#, 作者简介：杨