二维声波方程全波形反演方法研究

发布时间：2020-11-06 06:22

　　 全波形反演(Full Waveform Inversion,简称FWI)具有高精度、高分辨率的特点,其结果可以为岩性判断、油气藏识别等勘探问题提供重要支撑。本文给出了一种新的二维时间域声波方程的全波形反演方法。我们从正演模拟方法、迭代方向和搜索步长这三个方面对全波形反演方法进行了研究:(1)使用 ONADM(Optimal Nearly Analytic Discrete Method)方法对层状模型和Marmousi模型进行正演模拟,从模拟结果的地表地震记录可知:两个模型的直达波和反射波明显、整个地震记录清晰且无可见数值频散出现,可见ONADM方法可以用于全波形反演中;(2)选择CG-DESCENT方法作为迭代方向的选取方法;(3)选择改进后退法作为步长搜索方法。在此基础上,本文提出了一种新的全波形反演方法一CDIB(CG-DESCENT method and improved back-off method)全波形反演方法。使用CDIB全波形反演方法分别对层状模型和Marmousi模型进行反演试算。从层状模型的反演结果可知:在满足收敛条件的前提下,CDIB全波形反演方法迭代了 50次即可得到较好的反演结果。在对Marmousi模型反演过程中,我们同时使用L-BFGS算法联合Wolfe-Powell非精确一维搜索方法(简称LBWP全波形反演方法)反演。通过两种方法的反演结果作对比,我们发现:(1)在满足收敛条件的前提下,CDIB全波形反演方法迭代了 68次,而LBWP全波形反演方法迭代了 96次,且LBWP全波形反演方法所需计算时间是CDIB全波形反演方法的4.61倍,故可知CDIB全波形反演方法具有收敛速度快、计算时间少、计算效率高的优点;(2)CDIB全波形反演方法浅层反演较好,而LBWP全波形反演方法深部反演较好。数值实验说明了本文提出的CDIB全波形反演方法具有一定的有效性。
【学位单位】：西南石油大学
【学位级别】：硕士
【学位年份】：2018
【中图分类】：O241.8
【部分图文】：

图２－２层状模型１．Ｉｓ的地表地震记录图??注：图屮横十标为水平位置，中．位为ｋｍ：纵肀标为传播时间，ｉ＇ｆｉ．位为Ｓ

图２－４?Ｍａｒｍｏｕｓｉ模型５ｓ的地表地震记录??注：图屮横坐标为水平位置，单位为ｋｎｖ，纵唯标为传播时间，单位为ｓ

图３－２最速下降法和共轭梯度法对比图??注：图中（ｆｌ）最速下降法迭代过程变量变化示意图，呈“锯齿形＼?共轭梯度法迭代过程变量变化示意??图
【参考文献】

相关期刊论文 前8条

1 冯英杰;杨长春;吴萍;;地震波有限差分模拟综述[J];地球物理学进展;2007年02期

2 丁继才;常旭;刘伊克;汪长永;;反射地震数据的逐层波形反演[J];地球物理学报;2007年02期

3 佘德平;波场数值模拟技术[J];勘探地球物理进展;2004年01期

4 张美根,王妙月,李小凡,杨晓春;时间域全波场各向异性弹性参数反演[J];地球物理学报;2003年01期

5 杨顶辉;双相各向异性介质中弹性波方程的有限元解法及波场模拟[J];地球物理学报;2002年04期

6 孟鸿鹰,刘贵忠;小波变换多尺度地震波形反演[J];地球物理学报;1999年02期

7 杨顶辉，滕吉文，张中杰;三分量地震波场的近似解析离散模拟技术[J];地球物理学报;1996年S1期

8 符力耘，牟永光;弹性波边界元法正演模拟[J];地球物理学报;1994年04期

相关博士学位论文 前1条

1 马啸;波动方程保辛近似解析离散化算法研究[D];清华大学;2013年

相关硕士学位论文 前4条

1 崔海娟;改进共轭梯度法求解无约束优化问题[D];渤海大学;2014年

2 陈永芮;时间域全波形反演方法研究[D];中国石油大学(华东);2014年

3 倪长宽;复杂介质地震波场正演模拟方法研究[D];中国石油大学;2008年

4 宋国杰;三维弹性波方程的改进近似解析离散化方法及波场模拟[D];清华大学;2008年

本文编号：2872771