反演原理及公式介绍

《反演原理及公式介绍》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、第一章反演理论第一节基本概念一．反演和正演1．反演反演是一个很广的概念，根据地震波场、地球自由振荡、交变电磁场、重力场以及热学等地球物理观测数据去推测地球内部的结构形态及物质成分，来定量计算各种有关的物理参数，这些都可以归结为反演问题。在地震勘探中，反演的一个重要应用就是由地震记录得到波阻抗。有反演，还有正演。要正确理解反演问题，还要知道正演的概念。2．正演正演和反演相反，它是对一个假设的地质模型，给定某些参数（如速度、层数、厚度）用理论关系式（数学模型）推导出某种可测量的量（如地震波）。在地震勘探中，正演

2、的一个重要应用就是制作合成地震记录。3．例子考虑地球内部的温度分布，假定地球内部的温度随深度线性增加，其关系式可表示成：T(z)=a+bz正演：给定a和b，求不同深度z的对应温度T(z)反演：已经在不同点z测得T(z)，求a和b。二．反演问题描述和公式表达的几个重要问题1．应用哪种参数化方式——离散的还是连续的？2．地球物理数据的性质是什么？观测中的误差是什么？3．问题能不能作为数学问题提出，如果能够，它是不是适定的？4．对问题有无物理约束？5．能获得什么类型的解，达到什么精度？要求得到近似解、解的范围、还

3、是精确解？6．问题是线性的还是非线性的？7．问题是欠定的、超定的、还是适定的？8．什么是问题的最好解法？9．解的置信界限是什么？能否用其它方法来评价？第二节反演的数学基础一．解超定线性反问题151．简单线性回归可利用最小平方法确定参数a、b使误差的平方和最小。（1-2-1）拟合公式为：（1-2-2）该方法的公式原来只适用于解超定问题，但同样适用于欠定问题，当我们有多个参数时，称为多元回归，在地球物理领域广泛采用这种方法。此过程用矩阵形式表示，则称为广义最小平方法矩阵方演。2．非约束最小平方法反演——广义矩阵

4、方法由前面讨论可知，参数估计的最小平方方法用矩阵公式表示，所得到的算法等价于一个或多个模型参数的一个或多个数据集反演，步骤为：问题定义→矩阵公式→最小平方解线性问题采用广义矩阵形式d=Gm（1-2-3）对于精确的数据模型，参数m为m=G-1d（1-2-4）但是由于试验误差，实际数据将不能精确拟合获得，故采用最小平方法求解。解的矩阵表示式为（1-2-5）上式具体计算时可用奇异值分解方法G=U∧VT最后，得=（GTG）-1GTd=V∧-1UTd（1-2-6）一．约束线性最小平方反演15为了得到最合适的解，通常可

5、在方程d=Gm中加先验信息，进行约束反演。约束方程为Dm=h（1-2-7）D一般为只有对角线有值的矩阵，我们希望朝着偏置使得最小。=（d-Gmd-Gm）+β2（Dm-hDm-h）（1-2-8）如果D是单位矩阵，可以得到约束解=（GTG+β2I）（GTd+β2h）（1-2-9）式中，β称为Lagrange乘子。三．解非线性反演问题1．思路在实际工作中许多问题都是非线性的，而非线性问题求解通常比较复杂，这样就产生这样一个问题，给定一些非线性问题，而它们又不服从简单的线性变换，那么能否用通用的方法使我们可以用一些

6、线性反演的方法来估算未知模型参数，并最终求得问题的解决呢？答案是肯定的。2．初始模型和线性化对于非线性问题di=fi（m1，m2，…mp）=fi（m），i=1，2，…n（1-2-10）设m0为初始模型，则其响应为（1-2-11）现假定f（m）在m0附近是线性的，从而关于m0的模型响应的微小摄动可以用Taylor级数展开为15或简记为实际情况要考虑噪声d=f（m）+e（1-2-12）令y=d-f（m0），，则有e=d-=y-Ax（1-2-13）e=y-Ax这样，非线性问题转化成线性问题，我们可以用线性的方法求

7、出问题的解。四、无约束非线性反演1．问题的公式化目标函数：q=eTe=(d-f(m))T(d-f(m))（1-2-14）利用前述结果，上式改写为q=eTe=(y-Ax)T(y-Ax)（1-2-15）2．问题的解法：Gauss-Newton法对参数摄动的最小平方解（1-2-16）将摄动（x=δm）应用于起始模型m0，迭代公式如下：(1-2-17)其中mk为Jacobian矩阵A的赋值。3．Gauss-Newton法的局限性当AT15A病态（本征值很小或近于0）时，计算的解会大到令人难以置信。因此在实践当中，必

8、须对mk做x的微小校正。4．最速下降（梯度）法初始模型仅在目标函数q的负梯度方向予以校正，即（1-2-18）其中k是合适的常数，进一步推导可得（1-2-19）以上方程中以[ATA]-1取代常数因子2k，将变为方程1-2-16所定义的Gauss-Newton法，k值决定校正步长。但以上方程并不含有任何逆矩阵，因此较Gauss-Newton法具备更好的起始收敛特征。最速下降法当采用最小平方解法时，其收敛速率将下降，因