变网格有限差分波场数值模拟的研究

发布时间：2020-03-23 02:01

【摘要】：地震波场数值模拟是勘探地震学的重要研究内容之一,高精度、高效率的数值模拟是研究波动现象的重要手段,在油气藏勘探开发中发挥着重要作用。随着油气藏勘探的逐步发展,地震勘探的重心已从浅中层构造、简单粗糙构造向复杂深层非均质构造偏移,在这种情况下既能对这些精细(低速)岩层准确刻画,又能保证整个模拟的计算效率就显得尤为重要。传统的常规网格有限差分法都是基于相同大小的网格对计算区域进行划分,即整个模拟区域的采样步长是固定的。在模拟这种精细岩层时,为了减小数值频散和提高模拟精度,整个模型需要用较小步长的网格来划分,这样在计算波场时每个波长上的采样点数就会增加,导致计算量增大而降低模拟效率。因此对于复杂地质构造,寻求一种更加灵活适应性更强的波数值模拟方法是非常必要的。论文在研究基于声波方程的常规(交错)网格的高阶有限差分算法的基础上,分析了有限差分正演模拟的稳定性和数值频散问题。针对常规网格有限差分法进行波场数值模拟时存在频散和模拟效率相矛盾的问题,提出采用变步长交错网格有限差分算法进行数值模拟,该算法首先根据介质速度的变化采用不同步长的网格进行模拟区域的划分,然后利用交错网格算法计算波场值,针对计算时出现的网格过渡区域,通过特定的边界区域网格步长的变化推导出波场差分公式,在一个循环内完成对所有过渡带网格点的波场计算,实现模拟区域从大(小)网格到小(大)网格的平滑稳定过渡。针对计算机数值模拟引入的人工边界问题,通过多次透射边界条件,Clayton吸收边界条件以及完全匹配层(PML)边界条件的对比分析,发现进行波场数值模拟时采用PML边界条件计算效率高且对边界反射吸收效果好。最后,针对地震波场建立不同的介质模型进行波场数值模拟,从模拟精度、计算效率和内存需求三个方面与常规(交错)网格算法进行对比,实验分析结果表明,变步长交错网格算法在计算效率、模拟精度等方面优于常规(交错)网格算法。论文的研究内容为后续地震勘探工作的顺利进行及地震反演的研究奠定了坚实的理论和技术基础。
【图文】：

西安工程大学硕士学位论文最常用的是矩形网格；第二是将微分方程转网格点的函数值。差分离散的途径主要有两接转化为时间和空间的二阶中心差分离散求方法计算简单，，程序易于实现，但差分精度双曲线方程组用应力和速度的交错网格格式度，目前，主流的差分算法就是基于交错网

图 2-2 速度-应力交错网格示意时间上 2 N阶的差分近似，应力和位移速度分别是在 tTt 和2 时刻差度时间差分近似式。力：)2((21)!1222121121,21,ttunuunnNnnkijkij 度：)2((21)!12221211,1,nnNnnkijkijttvnv v =2 时，时间二阶精度的差分算子为：uuukijkijk 21,21,
【学位授予单位】：西安工程大学
【学位级别】：硕士
【学位授予年份】：2018
【分类号】：O241.3

【参考文献】

相关期刊论文 前10条

1 孟凡顺;张杰;李洋森;;复杂介质多尺度变网格地震波数值模拟[J];地球物理学进展;2015年01期

2 孟凡顺;李洋森;王璐;邓志勇;张亮;;基于可变交错网格的地震波数值模拟[J];地球物理学进展;2013年06期

3 王维红;柯璇;裴江云;;完全匹配层吸收边界条件应用研究[J];地球物理学进展;2013年05期

4 胡昊;刘伊克;常旭;王一博;杜向东;杨仁虎;;逆时偏移计算中的边界处理分析及应用[J];地球物理学报;2013年06期

5 孙成禹;丁玉才;;波动方程有限差分双变网格算法的精度分析[J];石油地球物理勘探;2012年04期

6 李国平;姚逢昌;石玉梅;黄文锋;程利敏;;有限差分法地震波数值模拟的几个关键问题[J];地球物理学进展;2011年02期

7 黄超;董良国;;可变网格与局部时间步长的高阶差分地震波数值模拟[J];地球物理学报;2009年01期

8 孙成禹;李胜军;倪长宽;张玉华;;波动方程变网格步长有限差分数值模拟[J];石油物探;2008年02期

9 朱生旺,魏修成;波动方程非规则网格任意阶精度差分法正演[J];石油地球物理勘探;2005年02期

10 裴正林,牟永光;地震波传播数值模拟[J];地球物理学进展;2004年04期

相关博士学位论文 前1条

1 杨宇;场地地震波动模拟中透射边界稳定性问题研究[D];中国地震局地球物理研究所;2011年

相关硕士学位论文 前3条

1 王馨;2D各向同性和各向异性介质中高阶交错网格有限差分波场模拟[D];长安大学;2014年

2 李洋森;复杂介质可变网格地震波数值模拟研究[D];中国海洋大学;2013年

3 武英婷;交错网格高阶有限差分法的弹性波波场数值模拟[D];西安理工大学;2010年

本文编号：2595991