反演原理及公式介绍.docVIP

  1. 1、本文档共15页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第一章 反演理论 第一节 基本概念 反演和正演 反演 反演是一个很广的概念,根据地震波场、地球自由振荡、交变电磁场、重力场以及热学等地球物理观测数据去推测地球内部的结构形态及物质成分,来定量计算各种有关的物理参数,这些都可以归结为反演问题。在地震勘探中,反演的一个重要应用就是由地震记录得到波阻抗。 有反演,还有正演。要正确理解反演问题,还要知道正演的概念。 正演 正演和反演相反,它是对一个假设的地质模型,给定某些参数(如速度、层数、厚度)用理论关系式(数学模型)推导出某种可测量的量(如地震波)。在地震勘探中,正演的一个重要应用就是制作合成地震记录。 例子 考虑地球内部的温度分布,假定地球内部的温度随深度线性增加,其关系式可表示成: T(z)=a+bz 正演:给定a和b,求不同深度z的对应温度T(z) 反演:已经在不同点z测得T(z),求a和b。 二.反演问题描述和公式表达的几个重要问题 1.应用哪种参数化方式——离散的还是连续的? 2.地球物理数据的性质是什么?观测中的误差是什么? 3.问题能不能作为数学问题提出,如果能够,它是不是适定的? 4.对问题有无物理约束? 5.能获得什么类型的解,达到什么精度?要求得到近似解、解的范围、还是精确解? 6.问题是线性的还是非线性的? 7.问题是欠定的、超定的、还是适定的? 8.什么是问题的最好解法? 9.解的置信界限是什么?能否用其它方法来评价? 第二节 反演的数学基础 一.解超定线性反问题 1.简单线性回归 可利用最小平方法确定参数a、b使误差的平方和最小。 (1-2-1) 拟合公式为: (1-2-2) 该方法的公式原来只适用于解超定问题,但同样适用于欠定问题,当我们有多个参数时,称为多元回归,在地球物理领域广泛采用这种方法。此过程用矩阵形式表示,则称为广义最小平方法矩阵方演。 2.非约束最小平方法反演——广义矩阵方法 由前面讨论可知,参数估计的最小平方方法用矩阵公式表示,所得到的算法等价于一个或多个模型参数的一个或多个数据集反演,步骤为: 问题定义→矩阵公式→最小平方解 线性问题采用广义矩阵形式 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) 约束线性最小平方反演 为了得到最合适的解,通常可在方程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.思路 在实际工作中许多问题都是非线性的,而非线性问题求解通常比较复杂,这样就产生这样一个问题,给定一些非线性问题,而它们又不服从简单的线性变换,那么能否用通用的方法使我们可以用一些线性反演的方法来估算未知模型参数,并最终求得问题的解决呢?答案是肯定的。 2.初始模型和线性化 对于非线性问题 di=fi(m1,m2,…mp)=fi(m), i=1,2,…n (1-2-10) 设m0为初始模型,则其响应为 (1-2-11) 现假定f(m)在m0附近是线性的,从而关于m0的模型响应的微小摄动可以用Taylor级数展开为 或简记为 实际情况要考虑噪声 d=f(m)+e (1-2-12) 令y=d-f(m0),,则有 e=d-=y-Ax (1-2-13) e=y-Ax 这样,非线性问题转化成线性问题,我们可以用线性的方法求出问题的解。 四、无约

文档评论(0)

185****7617 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档