大地电磁测深一维反演解释.docxVIP

大地电磁测深一维反演解释 0 mt数据反演和一维反演 图根于20世纪50年代提出了图根河的地震计算方法。这是一种利用自然交变电阻器研究地球电性结构的地球物理勘探方法。由于它不用人工供电,成本低,工作方便,不受高阻层屏蔽,对低阻层分辨率高,且勘探深度仅与电磁场的频率有关,浅可达几十米,深可达数百公里,因此在许多领域内得到了广泛的应用。MT资料反演解释就是根据地表测得的地球电磁场响应,如视电阻率、阻抗相位、表面阻抗、倾子等,通过一定的优化处理求得一个合理的地电模型。目前二维数据处理技术日趋成熟,比较典型的二维反演解释方法有: OCCAM法、 RRI法、 REBOCC法、 NLCG法和ABIC法。同时,一维反演方法也很多,应用广泛的有:OCCAM反演、一维连续介质反演的曲线对比法、“正演修正法”一维反演、二次函数逼近非线性反演和自适应正则化反演。 另外,一维反演通常是二维反演计算的第一步骤,并且可以采用一维反演方法对地质模型进行拟二维反演解释,往往能对MT数据解释带来有用的效果。拟二维反演方法的精度比不上二维反演,但由于其反演速度快、内存需求少,所以在实际工作中,可以对野外大地电磁测深数据进行初步解释。因此,使用一维反演方法进行拟二维反演解释还是非常有意义的。 作者在本文采用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率,并将Bostick反演结果作为反演计算的初始模型。 1 地表波阻抗zi 假定地电剖面是水平分层均匀的,如果地电剖面共有N层,则共有2N-1个参数,hi(i=1,2,…,N-1)和ρi(i=1,2,…,N)分别代表第i层的厚度和电阻率。 对于上述的一维层状介质模型,计算视电阻率ρa和相位?a的公式如式(1)。 ρa=1ωμ|Z1|2?a=arctanIm[Z1]Re[Z1](1)ρa=1ωμ|Ζ1|2?a=arctanΙm[Ζ1]Re[Ζ1](1) 其中Z1表示第一层地表的波阻抗;μ为导磁率;ω=2πTω=2πΤ为角频率;Z1可用下面的递推公式计算。 Zi=Z0iZ0i(1?e?2kihi)+Zi+1(1+e?2kihi)Z0i(1+e?2kihi)+Zi+1(1?e?2kihi)ZN=ωμkN=Z0N(2)Ζi=Ζ0iΖ0i(1-e-2kihi)+Ζi+1(1+e-2kihi)Ζ0i(1+e-2kihi)+Ζi+1(1-e-2kihi)ΖΝ=ωμkΝ=Ζ0Ν(2) 其中ki=iωμρi???√ki=iωμρi是第i层的复传播系数;Z0i为第i层的特征阻抗;Zi是第i层顶面的波阻抗。 由此可见,对N层地电断面而言,若固定每一层的厚度,则视电阻率ρa和相位?a可表示为关于信号周期及层电阻率参数之间的函数,即 ρa=f(ρ1,ρ2,…,ρN,T) (3) ?a=f(ρ1,ρ2,…,ρN,T) (4) 2 mt方案的二维反演 2.1 模型约束目标函数 大地电磁反演问题,可以抽象地描述为已知观测数据求之于对应模型的过程。设d是观测数据向量,m是模型参数向量, F为把地球模型映射到理想数据的函数,则 d=F(m) (5) 式中 为正演函数。 大地电磁反演问题是不适定的,其反演结果是不稳定的,并且具有非唯一性,这也就意味着不同的地电模型对观测数据的拟合具有同样的精度。在地球物理反演中,为了改善解的稳定性和非唯一性问题,通常是引入Tikhonov的正则化思想,见式(6)。 Pα(m)=?(m)+αs(m) (6) 其中Pα(m)为总目标函数;α为正则化因子;?(m)为观测数据与预测数据之差的平方和(即数据目标函数);s(m)为稳定器(即模型约束目标函数),这里采用基于先验模型的最小模型约束。 因此,大地电磁反演问题的总目标函数表示为式(7)。 Pα(m)=‖Wd[dobs-F(m)]‖2+ α‖Wm(m-mref)‖2(7) 式中Wd为数据权系数矩阵;Wm为模型权系数矩阵;mref为先验模型。 将F(m)用泰勒公式展开为式(8)。 F(mk+△m)=F(mk)+Jk△m+ O‖(△m)2‖ (8) 其中mk为模型的第k次迭代值。 于是可得: dk+1≈dk+Jk△m(9) 这里dk=F(mk);△m=mk+1-mk;Jk是雅克比灵敏度矩阵。 于是有: pα(mk+1)=‖Wd(dobs-dk-Jk△m)‖2+ α‖Wm(mk-mref)‖2(10) 将上述目标函数对△m求导并令其等于0,可得线性方程组(11)。 (JkTWTdWdJk+αWTmWm)△m= JkTWTdWd(dobs-dk)+ αWTmWm(mref-mk) (11) 解方程组式(11),可得到模型修正量△m,将其加到预测模型参数矢量中,可以得到新的模型参数矢量。再重复该过程,直