频率域波动方程高精度有限差分格式及并行模拟算法研究

发布时间：2020-07-12 00:22

【摘要】：随着易采矿藏的逐步衰竭和科技的进步,勘探技术对效率和精度的要求越来越高。现有先进地震成像技术大多需要一种高效精确的正演模拟算法,而频率域波动方程求解则是其中之一。波动方程可由有限差分法、有限元法、不连续Galerkin法或谱元法等来进行离散。有限差分法以其简单及计算效率高在地震勘探领域得到了广泛的应用。设计一种高精度的有限差分方法可以降低空间采样精细度或提高模拟频率,这对提高成像精度大有裨益。此外,由于离散形式的频率域正演模拟等同于求解一个大型稀疏线性方程组,因此一个快速稳定的解法对许多频率域内地震成像技术(如全波形反演和逆时偏移成像等)的实现也十分重要。考虑到地下介质高度的非均匀性,各向异性以及其他复杂的性质,由波动方程离散化导出的线性方程组系数矩阵不定且高度病态,求解此类方程组是一个巨大的挑战。由此,本文基于以上两点做了如下工作:首先,针对频率域声波正演模拟问题探讨了二阶差分格式,混合网格差分格式和四阶差分格式,在理论上分析了这几种格式的收敛阶,并通过数值实验证实了相应的收敛阶。将数值解与精确解作对比,进一步考察了不同格式的绝对误差。通过对比求解所需计算时间,分析了三种格式各自的执行效率。结果显示混合网格格式是其他两种格式的一个折中,在保证数值解较高精度的同时仍能将计算量维持在与二阶精度相近的水平上。其次,对于二维频率域弹性波场模拟,从二阶交错网格格式出发,通过与时间域地震记录的对比验证了格式的精确性。在四阶交错网格格式基础上引入反集中质量技巧,在频散分析的框架下通过极小化归一化相速度与单位一之间的偏差,给出了最优的四阶交错网格格式。该格式的频散较经典四阶格式有了大幅的减少,因而可以在单位波长内采取更少的网格点而同时保证数值解的精度,这将会大大降低波场模拟的计算量。通过在复杂的Marmousi2模型中进行模拟并将结果与时间域地震记录作对比,验证了该格式的精确性。再次,为了得到三维情况下最优的四阶交错网格差分格式,建立三维(粘)弹性波动方程频散分析理论框架,分别应用基于梯度的Levenberg-Marquardt算法和结合了下山单纯形法的全局收敛模拟退火算法来求解归一化相速度和单位一之间偏差泛函的极小化问题,确定了最优的质量加权平均系数和最优四阶差分系数。频散曲线分析指出该最优格式受Poisson比影响较小,在群速度误差上限设置为1%时,单位波长内仅需3.7个网格点即可。定量地分析了最优格式和经典格式在应用直接法和迭代法时内存耗用的差异,指出了最优格式可大大降低波场模拟的内存耗用量。通过与声波、弹性波方程精确解的对比,在数值上验证了最优格式的高精度,即使在单位波长内仅有3.3个网格点的情况下,最优格式仍能较准确地匹配精确解。最后,在二、三维弹性波场模拟中应用高度并行的CARP-CG方法以应对大规模计算带来的挑战。通过数值实验研究了Poisson比、自由表面边界条件和地震衰减效应对该迭代法收敛性的影响。针对并行计算问题,设计了多频不同规模的波场模拟以考察CARP-CG方法在多核计算时的可扩展性,讨论了该方法在大规模波场模拟中的可行性。通过分析计算复杂度,探讨了CARP-CG方法应用于频率域波场模拟相较于直接在时间域进行模拟的优势。将CARP-CG方法与GMRES,CGNR和BiCGSTAB方法做对比,指出了CARP-CG方法对于不同的介质参数和实验设定均能收敛,展现了极好的稳定性。由此,CARP-CG方法可为地震成像技术提供一个高效的大规模正演模拟算法。
【学位授予单位】：哈尔滨工业大学
【学位级别】：博士
【学位授予年份】：2016
【分类号】：O241.8

【相似文献】

相关期刊论文 前10条

1 王秀琴;徐琛梅;;一类非线性偏微分方程的有限差分格式的稳定性研究(英文)[J];数学季刊;2009年03期

2 余跃玉;胡兵;闵心畅;;非线性Sobolev-Galpern方程的一种有限差分格式[J];四川大学学报(自然科学版);2011年01期

3 潘书勤;冉政;;有限差分格式的群分类[J];上海大学学报(自然科学版);2012年01期

4 周朦,方国洪;二维长波方程的一个无条件稳定有限差分格式[J];海洋与湖沼;1988年02期

5 杨云,王元明;求解非线性反应扩散方程的有限差分格式(英文)[J];华东师范大学学报(自然科学版);2002年04期

6 张厚柱;张宇;孙正;;高维波动方程数值模拟的隐式分裂有限差分格式[J];石油物探;2007年06期

7 徐琛梅;杨巍纳;;一类非线性偏微分方程的有限差分格式的稳定性问题[J];河南大学学报(自然科学版);2008年03期

8 刘欣;黄凯美;李福乐;;一类非线性抛物型方程的有限差分格式[J];青岛农业大学学报(自然科学版);2009年02期

9 M.雷里;R.伦高尼;V.彭莱梯;姚敬之;;二维有限差分格式的简洁分析表达式[J];华水科技情报;1985年04期

10 D.L.罗伯特;M.S.赛里姆;陈玉田;;解一维抛物型偏微分方程的六种显式和两种隐式有限差分格式的比较[J];华水科技情报;1985年04期

相关会议论文 前10条

1 申义庆;杨国伟;高智;;高分辨率有限紧致格式[A];中国力学学会学术大会'2005论文摘要集（下）[C];2005年

2 王焕军;林年添;郝继杰;孙剑;刘建龙;;时间—频率域的极化分析及其应用研究[A];中国地球物理学会第二十七届年会论文集[C];2011年

3 曹思远;牟永光;;时间-频率域分析综述[A];1993年中国地球物理学会第九届学术年会论文集[C];1993年

4 王博;雷军;;频率域信号分离设计及其在环境振动研究中的应用[A];中国地球物理2013——第十八专题论文集[C];2013年

5 许琨;;频率域最小二乘反演[A];工程地质力学创新与发展暨工程地质研究室成立50周年学术研讨会论文集[C];2008年

6 阮百尧;村上裕;徐世浙;;三维地形频率域垂直磁偶极子场边界元模拟方法[A];2000年中国地球物理学会年刊——中国地球物理学会第十六届年会论文集[C];2000年

7 阮百尧;;三维地形频率域人工源电磁场的边界元模拟方法[A];中国地球物理学会第二十届年会论文集[C];2004年

8 杨顶辉;滕吉文;张中杰;;二维各向异性介质数值模拟中的TVD方法[A];1995年中国地球物理学会第十一届学术年会论文集[C];1995年

9 高智;李明军;朱力立;;对流扩散方程的变步长摄动有限差分格式[A];计算流体力学研究进展——第十一届全国计算流体力学会议论文集[C];2002年

10 曹书红;陈景波;;声波方程频率域高精度正演的17点格式及数值实现[A];中国科学院地质与地球物理研究所2012年度（第12届）学术论文汇编——油气资源研究室[C];2013年

相关博士学位论文 前10条

1 李扬;频率域波动方程高精度有限差分格式及并行模拟算法研究[D];哈尔滨工业大学;2016年

2 李婷婷;基于高阶有限差分格式的Inverse Lax-Wendroff方法及其稳定性分析[D];中国科学技术大学;2017年

3 吕永强;Landau-Lifshitz方程的有限差分格式与整体正则解[D];中国工程物理研究院;2004年

4 张蕾;几类偏微分方程非标准有限差分格式的研究[D];哈尔滨工业大学;2014年

5 张磊;非标准有限差分法和参数扰动法在微分方程上的应用[D];哈尔滨工业大学;2012年

6 纪晓琳;频率域等效源曲面位场数据处理和转换方法研究[D];长安大学;2015年

7 王瑞文;流体力学及海洋数值模拟基于POD技术的降维方法研究[D];首都师范大学;2007年

8 李文杰;频率域航空电磁数据处理技术研究[D];中国地质大学（北京）;2008年

9 沈俊;生物型体竞争模型的高分辨率高精度方法[D];中国科学技术大学;2007年

10 杜娟;流体力学方程基于POD方法的降维数值解法研究[D];北京交通大学;2011年

相关硕士学位论文 前10条

1 赵小红;一类非线性偏微分方程的若干有限差分格式[D];燕山大学;2008年

2 赵丽梅;Cahn-Hilliard方程的几个有限差分格式的稳定性和收敛性分析[D];南京信息工程大学;2016年

3 李晓旺;两类阻尼波动方程的有限差分格式[D];山西大学;2016年

4 张磊;二维Ginzburg-Landau方程有限差分格式的长时间性态[D];黑龙江大学;2007年

5 王珏;带有多项式非线性项的三维反应扩散方程有限差分格式的长时间行为[D];黑龙江大学;2006年

6 董明华;双曲守恒方程的一类高精度算法[D];哈尔滨工业大学;2011年

7 李福乐;一类带有边界反馈的Timoshenko梁方程的有限差分格式[D];东南大学;2005年

8 常红;Camassa-Holm方程的有限差分格式[D];江苏大学;2009年

9 胡迪;带阻尼项的Sine-Gordon方程有限差分格式的长时间性态[D];黑龙江大学;2010年

10 国萃;某些非线性方程数值解的长时间性态[D];黑龙江大学;2005年

本文编号：2751141