利用迭代反演方法解决变磁性磁界面反演问题.pdfVIP

1 中国地球物理201 ·975· (26)海洋地球物理学 利用迭代反演方法解决变磁性磁界面反演问题 lterativeInversionMethodtoResolvetheInverfionProblem Using of Interface VaryingMagnetisation 张向宇1吴健生2 Wu ZhangXiangyuJiansheng 1．广州海洋地质调查局 广州510760；2同济大学海洋地质国家重点实验室上海200092 近年来，由于地震勘探方法的限制，使得南海基底构造及构造演化推断工作受到影响，因此随着勘探 的深入，重磁勘探方法得到了进一步重视。同时近些年来，也有多项相关重磁研究在南海区域进行，我 们可以通过磁界面的反演计算来揭示基础层底界面，从而加深对这一地区地壳上部结构的认识。 南海地区经历了多次构造运动，使得地层岩石磁性在横向上呈现较大差异。目前常用的磁界面反演方 法如欧拉算法、维纳反褶积法以及希尔伯特变换法等多是针对常磁性磁界面的计算，若直接将其应用于南 海北部磁性界面的反演中必定会带进很大误差，这就给常规磁界面反演方法的应用提出了挑战。因此通过 对常用的界面迭代反演方法进行改进，获得适用于变磁性磁界面的迭代反演方法，探索得到一个可用于解 决此类情况的方法手段，为今后的类似研究提供一个思路。其次磁性界面可能与基底界面有关联，这样通 过反演磁性界面可以得到基底的情况，为进一步的地质一地球物理解释工作提供了依据。 对于一底界面以下磁化强度为零的磁界面，我们可以将其看作由N(N0)个直立长方体组合而成，并近似 认为每一个长方体内的高度和磁化强度为常数，则观测面上任一观测点(工，y，z)处的垂直磁化垂直磁异常值为： ^， ，c- 、， 、 ． ^ ． Za上(x，Y，z)。毪罗彬arctan垒羔三冬掣lm+差r一善r一量(1) 斗刀百 ，．L台一z， 抽★一于．tot一亏二o^一号 长方体的中心坐标；ak、bk、Q分别为第k个长方体的长、宽、高；r 实际进行计算时，已知界面磁异常和磁化强度值，便可通过(1)式，利用矩阵求取界面深度，完成反演工 作。从这个过程不难发现，若要达到可以将每一个长方体的磁化强度值和高度值看做不变的近似条件，则划分的 长方体数N要较大才可，而随着N的增大，矩阵计算的工作量也大大增加。若界面的磁化强度值和界面起伏都 变化不大的时候，我们便可以做近一步近似，取滑动窗口，将窗口处的界面深度值和磁化强度值看做是恒定的， 为窗口内所有长方体的深度和磁化强度的平均值，在每一个窗口处通过解矩阵求得界面深度值，当窗口滑动遍历 整个计算区域，就得到了整个界面的反演结果。这样通过引入滑动窗口的方式，将大大减小反演计算量。 本文采用的变磁性磁界面迭代反演方法就是建立在这个思想上，对常磁性磁界面迭代反演方法改进得 到的。首先输入研究区化极后的历磁异常数据和事先通过研究区地质资料得到的平均界面深度和磁化强 度值：为了减小边界影响，需要将异常下延到平均界面深度处。计算每一个窗口处对应平均界面埋深的磁 异常线性正演矩阵的解，得到每个滑动窗口处对应的界面深度，这样通过滑动窗口遍历整个研究区域就可 以得到整个区域的界面深度：通过对得到的界面深度值进行磁异常正演与输入的实际磁异常进行比较，判 断均方误差是否满足迭代结束条件，若不满足迭代结束条件，则求取界面深度修正值，返回继续进行迭代 计算直到满足迭代精度或迭代次数要求，当迭代结束后，就可得到界面深度反演结果。 为了对方法的可靠性进行验证说明，并且分析得到：①通过球体模型的验证，反演误差均在5％以内， 得到较理想的反演结果，说明该方法可靠性较好。②由于磁异常下延的影响，使该反演方法对地形