反演原理及公式介绍.docxVIP

第一章反演理论 第一节基本概念 反演和正演 反演 反演是一个很广的概念，根据地震波场、地球自由振荡、交变电磁场、重力场以及热学等地球物理观测数据去推测地球内部的结构形态及物质成分，来定量计算各种有关的物理参数，这些都可以归结为反演问题。在地震勘探中，反演的一个重要应用就是由地震记录得到波阻抗。 有反演，还有正演。要正确理解反演问题，还要知道正演的概念。 正演 正演和反演相反，它是对一个假设的地质模型，给定某些参数(如速度、层数、厚度)用理论关系式(数学模型)推导出某种可测量的量(如地震波】在地震勘探中，正演的一个重要应用就是制作合成地震记录。 例子 考虑地球内部的温度分布，假定地球内部的温度随深度线性增加，其关系式可表示成： T(z)=a+bz 正演：给定a和b，求不同深度z的对应温度T(z) 反演：已经在不同,曲测得T(z)，求a和b。 反演问题描述和公式表达的几个重要问题 应用哪种参数化方式一一离散的还是连续的？ 地球物理数据的性质是什么？观测中的误差是什么？ 问题能不能作为数学问题提出，如果能够，它是不是适定的？ 对问题有无物理约束？ 能获得什么类型的解，达到什么精度？要求得到近似解、解的范围、还是精确解？ 问题是线性的还是非线性的？ 问题是欠定的、超定的、还是适定的？ 什么是问题的最好解法？ 解的置信界限是什么？能否用其它方法来评价？ 第二节反演的数学基础 解超定线性反问题 1.简单线性回归 可利用最小平方法确定参麴、b使误差的平方和最小。 _Zy—b￡x__La一一y—bxn(1-2-1)s，一一一、b=(nE-ZxA)nEx2-(Zx)2 (1-2-1) 拟合公式为： y=a+bx(1-2-2) 该方法的公式原来只适用于解超定问题但同样适用于欠定问题，当我们有多个参数时，称为多元回归，在地球物理领域广泛采用这种方法。此过程用矩阵形式表示，则称为广义最小平方法矩阵方演。 非约束最小平方法反演一一广义矩阵方法 由前面讨论可知，参数估计的最小平方方法用矩阵公式表示，所得到的算法等价于一个或多个模型参数的一个或多个数据集反演，步骤为： 问题定义一矩阵公式一最小平方解 线性问题采用广义矩阵形式d=Gm(1-2-3) 对于精确的数据模型，参数n为m=G-1d(1-2-4) 但是由于试验误差，实际数据将不能精确拟合获得，故采用最小平方法求解。解的矩阵表示式为m=[GtG]-1Grd(1-2-5) 上式具体计算时可用奇异值分解方法G=UAVT 最后，得人，、，，、 m=(GtG)-1GTd=VA-1UTd(1-2-6) 约束线性最小平方反演 为了得到最合适的解，通常可在方fS=Gm中加先验信息，进行约束反演。约束方程为Dm=h(1-2-7) D一般为只有对角线有值的矩阵，我们希望朝着偏置m使得中最小。jj中=(d-Gm)^(d-Gm)+B2(Dm-h)r(Dm-h)(1-2-8) 如果D是单位矩阵，可以得到约束解m=(GTG+B2I)-1(GTd+B2h)(1-2-9)c 式中，B称为Lagrange乘子。 解非线性反演问题 思路 在实际工作中许多问题都是非线性的，而非线性问题求解通常比较复杂，这样就产生这样一个问题，给定一些非线性问题，而它们又不服从简单的线性变换，那么能否用通用的方法使我们可以用一些线性反演的方法来估算未知模型参数，并最终求得问题的解决呢？答案是肯定的。 初始模型和线性化 对于非线性问题d=f(m，m，…皿)=f(m)，i=1，2，???n(1-2-10)ii12pi 设mo为初始模型，则其响应为d0=/(m0)(1-2-11) 现假定f(m)在mo附近是线性的，从而关于mo的模型响应的微小摄动可以用Taylor级数展开为f(m)=f(mo+8m,mo+8m,mo+8m,…mo+8m) =f(mo)+直^8m+直^8m+直^8m++f?8m+高次项TOC\o1-5\h\z1dm问题的公式化目标函数：q=ee=(d—f(m))r(d—f(m))(1-2-14)利用前述结果，上式改写为q=Te=(y-Ax)T(y-Ax)(1-2-15)dm问题的解法：Gauss-Newton 问题的公式化目标函数：q=ee=(d—f(m))r(d—f(m))(1-2-14)利用前述结果，上式改写为q=Te=(y-Ax)T(y-Ax)(1-2-15) 问题的解法：Gauss-Newton法 对参数摄动的最小平方解X=(AtA)-1Aiy(1-2-16)将摄动(x=5m)应用于起始模型n。，迭代公式如下： 或简记为df帝)、f(m)=f(m)+[2f__)I8m[+0(118mll2)oIdmm=moj V.J=1.〕j实际情况要考虑噪声d=f(m)+e(1-2-12)df帝)、e=d-f(m)=d-f(mo)-]2■(m)I8m.