磁性界面反演方法.docxVIP

第九节磁性界面反演方法 具有一定磁性差异的地质界面，如结晶基底面、大岩体的上顶面等，是找矿 勘探与基础地质研究中常见的地质现象。 磁性界面反演方法是确定这一类地质模 型界面深度的方法。磁性界面反演方法有空间域和频率域两大类。 在磁性界面反 演中，常常把磁性界面划分成大量的离散二度水平棱柱体或三度直立棱柱体组合 模型，由于未知参数太多不能采用直接解法， 往往采用迭代法或其他方法。空间 域方法未知参数多、计算时间长、效率低。 Parker（ 1972, 1973, 1974）采用了 连续模型，得出了频率域重磁位场正反演的理论公式， Oldenburg（ 1974）把它 推广成迭代形式并做了二维计算。由于引入快速傅立叶变换，在相同精度下，频 率域方法比空间域方法反演速度至少要快一个数量级以上。 因此，频率域磁性界 面的反演方法成为界面反演的一种常用方法，用于区域磁测资料解释与油气勘探 中研究基底构造。 一、磁性界面异常的正演 如图7-9-1所示的磁性界面，其上下界面磁化强度差为 M，为简单起见，设 M垂直向下。 图7-9-1磁性界面示意图 若磁化率为常数，考虑n=0时，即泰勒展开式第一项在空间域为常数项，略 去n=0项则有： aZ（u,v） = 01M |e出送（一s） hn （7-9-13） 2 -心 n! - 上式表示，当给定了平均深度H及平均深度上的起伏h ,，取泰勒展开式 ■ n 有限项数n=3~8,就可以计算出h和Z（u,v），禾I」用快速傅立叶变换即可得到空 间域的磁异常值 Z （x,y,0。 二、磁性界面异常的反演 式（7-9-13）是磁性界面正演计算公式，稍作一下变化，就可以当作反演迭 代公式。我们把和式中n=1的项写出并移项得 ▼ 2 e 心旳-n (_s) h Z- (h) (7-9-14) 卩0 Ms n： n! 若令等式右端的h值为第i次近似值，而等式左端的h为第(i+1)次近似值, 即为Parker提出的迭代反演公式： (7-9-15)-i 1 -2 eHs ? _「： n (7-9-15) n!% Msh ~Z_ n! % Ms n z2 h值，而不是 h值，而不是h的修正量Ah，不容 易保证迭代收敛。可以采用如下的迭代方式： 取式(7-9-15)的一级近似值： (1) 2 eHs h Z (7-9-16) k0 Ms …(1)其中h表示h …(1) 其中h表示h的一级近似值 h⑴的频谱。 (1) 把一级近似值h代入式 (7-9-13)得 (1) Z (一s) (h(1)n) (7-9-17) ?⑵h把厶Z⑴按下式求h ?⑵ h ——L - Z(1)](7-9-18) J0 Ms r⑵ 把h反变换求得h⑵，则h(1) - h⑵为h的二次近似值，把它带入式(7-9-13) 求得厶Z的二级近似值AZ⑵。 按下式求h⑶: Hs (3) 2 e ⑵ h (u,v) [ ：Z(u,v)- ：Z (u,v)](7-9-19) % Ms 则h(1) - h⑵-h(3)为h的三次近似值，再继续求⑶ 如此下去，直到h的 第n次近似值与第(n-1)次近似值之差满足所要求的精度为止。这时得到的界 面起伏为 n h 八 h(i) ( 7-9-20) i d