SBI与RRI结合的MT二维反演方法研究.pdfVIP

The10mChinaInternational 2011 Geo-Electromagnetic Workshop SBI反演从模型构建上着手，把地质结构看着由有限数量具有均匀电导率的电性单元组成，而这些单元的深 度在横向上是可变的，模型参数用各单元的电阻率值及其下底边界深度表示，对于一个包含f层每层含r／个最小 单元的半空间模型，模型参数m可描述为【101： m=[Iog(p．)，Iog(p：)，．．．，log(p，)，in(d．。)’ ，l、 In(d．：)，．．．，ln(d，．)，．．．，ln(d。)，ln(d，：)，．．．，In(d。)】 式中肛为第f层的电阻率·以为第z层第／／个单元下底界面的深度，由于各参数的变化范围较大，为了反 演的稳定性，采用它们的对数值进行反演。 反演目标函数为 u=IlSnl|12+／z一01w(d一以m))112一右】 (2) 通过求导得到模型参数的迭代系列为 (3) mt+。=Lu(s7s)+(wJ。)7(wJ。)】．I(wJ。)7Wd。 对于给定的初始模型，可以通过这个迭代系列寻找极小可能构造模型12】。 2．2RRI反演基本原理 RRI反演通过解与一维相近的反演问题，来计算在每个测点位置下面的电阻率扰动，把二维反演问题转化 为一系列一维反演问题。 L． E一 出 对于TE,Ⅲ模式分别定义变量矿=亡孕oz=毗卺’c，=芒H警=鲁H作扰动分析’建立数据扰动和模型 参数扰动之间的线性积分方程： ％=志肌暗粉胁皿 (4) 帆=志肌培粉b啦 (5) 其中d丸和甜体分别为TE和TM模式下观测数据与理论数据的差值，ao(Z)为模型改变前的电导率值， 值，Ey,o(儿，z)和ex．o∽，z)是初始模型或者本次迭代前模型在第f个测点下某深度的理论电场。 综合考虑模型横向和垂向的不均匀性，构造如下目标函数： 蜘J1帮吲z，刳，。斋卜 ‘6’ 这是一个在各测点上的标度Laplace范数。式中，函数．，。(z)可以控制标度尺的长度，是用来度量不同深度的 构造，取／0)=m(z+％)，常数Zo通常是选取模型表层电阻率值和最高频率情况下的趋肤深度。 m=tn(cr)=-re(p)。gQ)是起控制水平方向构造的惩罚因子【’61。 2．3RR!与SBI结合反演 从上面原理可以看到，由于两种反演的模型参数不同，所以不能把两种方法以加权的形式放在一起加以计算。 因此我们先对数据作RRI反演，从半空间模型开始，得到一个较光滑的模型结果，记该结果为ml。将ml作为 SBI反演的初始模型进行SBI反演，得到反演结果m2，将m2作为最终反演结果。由于ml是一个以不同测点不 同深度电阻率表示的结果，而SBI反演的初始模型为块状结构模型，因此需要把ml转化为块状结构模型。我们 采用最新的模糊分割方法将ml转化为块状结构模型。模糊分割方法的原理如下： 将数据集分成L组·则可把数据集表示成：K(f=1，…，N；j=1，…，M)，在构建反演初始模型时，N表 示电阻率的个数，M通常等于3(电阻率的值、Y、Z的坐标)，L表示要构建初始模型的层数；那么数据集第，个 238 201 m十届叶】目目％地球电m≠}H*2 矢量在第f类的成姗j数目表示为U=[％]删(u也称分割矩阵)特数据集模糊分割成L娄的函数满足 以下条件： 0U。