地震层析成像doc
摘要:地震层析成像LTI-SIRT算法研究与分析摘要:本文研究反演二维近地表速度结构的初至波层析方法,正演采用LTI算法,反演采用SIRT算法。经过长时间的实践过程,本文总结出LTI-SIRT算法具体实现方法与步骤,并对LTI算法和SIRT算法进行了改进和发展。为了得到更好的成像结果,还结合平滑,插值,加权的处理手段,给出实验结果,分析实验误差,对算法作出评价。关键词:初至波;层析反演;LTI;SIRT;射线追踪1引言层析反演一般分两步:一是正演,即计算旅行时及射线追踪;二是反演求解,修正模型。不断地重复这两步,直到模型修正到满足一定的要求。射线追踪的算法有多种,比如最短路径法,有限差分方程法,旅行时线性插值法(LTI,LinearTraveltimeInterpolation)等.反演求解的方法有反投影法(BPT),代数重建法(ART),同时迭代重建法(SIRT),奇异值分解法(SVD),最小二乘QR分解法(LSQR)等多种算法。各种算法都有自己的优缺点,经综合考虑,本文采用的是射线追踪精度高的LTI法以及比较容易实现的SIRT法。以下给出本文利用LTI-SIRT进行层析反演的实现方法与步骤。2用LTI-SIRT算法进行层析反演的实现方法与具体步骤第一步:初至拾取;第二步:建立初始模型;第三步:射线追踪,找出地震波从炮点到接收点的传播路径,并算出理论走时;第四步:根据理论走时与实际拾取走时的差值来反演求解,对模型进行修正;第五步:对修正后的模型进行插值处理;第六步:对修正后的模型进行平滑;第七步:重复上述四个步骤,直到模型修正满足一定的要求;下面对以上步骤
温馨提示:当前文档最多只能预览
5 页,若文档总页数超出了
5 页,请下载原文档以浏览全部内容。
本文档由 匿名用户 于 2022-07-05 23:26:41上传分享