全波形反演中三种优化方法对比研究

发布时间：2020-05-28 09:13

【摘要】：随着勘探要求的不断提高,普通地震勘探的成像精度已经难以满足人们勘探的需要。全波形反演技术(Full waveform inversion)因具有刻画精细构造的潜力,而受到当今地球物理学家的追捧。虽然全波形反演在理论上几乎完美,但是在具体的实现过程中却困难重重。本文实现了基于一阶声波方程的全波形反演方法,研究过程中论文对以下部分进行了详细讨论:首先,论文对全波形反演中的正演部分进行了研究,详细推导了一阶声波方程的交错网格差分格式,讨论了PML吸收边界条件和CPML吸收边界条件的优劣点,以及交错网格有限差分的稳定性条件和频散条件。通过编程实现了简单模型和Marmousi模型的正演模拟,模拟结果表明在满足稳定性条件和频散条件下的时间二阶空间十六阶的交错网格有限差分正演具有较高的精度;通过CPML吸收边界条件和PML吸收边界条件的吸收效果对比,可以明显得出CPML吸收边界条件的吸收效果要优于PML吸收边界条件。其次,论文详细推导了基于伴随状态法的梯度求取公式,并对迭代步长的选取、多尺度反演策略和初始模型的选择进行讨论。通过编程实现了全波形反演中梯度的求取,验证了伴随状态法求取梯度的有效性。通过在异常体模型和Marmousi模型的反演结果,表明了本文所采取的方法和技术是可行的,多尺度反演策略可以提高反演精度,初始模型的选取对反演结果有很强的影响。此外,针对梯度优化问题,论文对L-BFGS算法进行了详细的公式推导,对比了最速下降法、共轭梯度方法和L-BFGS算法的全波形反演结果,可以明显得出基于L-BFGS算法的全波形反演结果具有最好的效果,基于共轭梯度法的反演结果次之,最速下降法相对前两者表现最差;针对全波形反演编程困难以及计算成本巨大的问题,本文讨论了新一代的Fortran语言特性和OpenMP并行技术在全波形反演中的应用效果,结果表明运用新一代的Fortran语言的特性具有易于编程,使代码便于维护的优点,对全波形反演的代码编写十分有利,利用OpenMP并行技术可以在很大的程度上减少计算成本。
【学位授予单位】：长安大学
【学位级别】：硕士
【学位授予年份】：2019
【分类号】：P631.4

【参考文献】

相关期刊论文 前10条

1 宋鹏;解闯;李金山;谭军;刘伟;谭惠文;;基于MPI+OpenMP的三维声波方程正演模拟[J];中国海洋大学学报(自然科学版);2015年09期

2 苗永康;;基于L-BFGS算法的时间域全波形反演[J];石油地球物理勘探;2015年03期

3 王义;董良国;;L-BFGS法时间域全波形反演中初始矩阵的选择方法[J];石油物探;2014年05期

4 张猛;王华忠;任浩然;冯波;隋志强;王延光;;基于CPU/GPU异构平台的全波形反演及其实用化分析[J];石油物探;2014年04期

5 张生强;刘春成;韩立国;杨小椿;;基于L-BFGS算法和同时激发震源的频率多尺度全波形反演[J];吉林大学学报(地球科学版);2013年03期

6 严红勇;刘洋;;黏弹TTI介质中旋转交错网格高阶有限差分数值模拟[J];地球物理学报;2012年04期

7 龙桂华;赵宇波;李小凡;高琴;王周;;三维交错网格有限差分地震波模拟的GPU集群实现[J];地球物理学进展;2011年06期

8 张军华;臧胜涛;单联瑜;石林光;;高性能计算的发展现状及趋势[J];石油地球物理勘探;2010年06期

9 卞爱飞;於文辉;周华伟;;频率域全波形反演方法研究进展[J];地球物理学进展;2010年03期

10 皮红梅;蒋先艺;刘财;王成祥;姜绍辉;;波动方程数值模拟的三种方法及对比[J];地球物理学进展;2009年02期

，

本文编号：2685026