尊敬的各位老师,上午好!
我叫___,是__级__班的学生,我的论文题目是粘声波正演模拟方法研究,论文是在__导师的悉心指点下完成的,在这里我向我的导师表示深深的谢意,向各位老师不辞辛苦参加我的论文答辩表示衷心的感谢,下面我将本论文设计的目的和主要内容向各位老师作一汇报,恳请各位老师批评指导。
首先,我想谈谈这个毕业论文设计的目的及意义
在地震勘探中,我们通常把地下介质为看成弹性各向同性介质。但是在实际中,地下介质为非完全弹性各向异性介质,存在吸收衰减现象。因此更趋向于粘弹性。粘弹介质正演模拟方法的研究对于还原地下介质的真实情况具有较强的理论意义和实际意义。由于粘弹性介质的正演模拟计算成本较高,计算繁琐,为了降低计算成本,本文进行了声学近似下的粘声波正演模拟方法研究。
下面是关于正演模拟方法的简单分类:
地震波正演模拟方法基于理论基础和表达形式可以分为三种:1>积分方程法(基于惠更斯原理);2>地震波方程数值解法(基于波动方程);3>射线追踪法(基于射线理论)
根据求解方法,地震波方程数值解法可分为有限元法、伪谱法、有限差分法。
有限元法算法比较复杂,计算较慢,而伪谱法不适合物性剧烈变化的复杂模型,相比于这两种方法,有限差分法可以适应剧烈变化的地下介质,算法速度较快。因此,此处选择有限差分法进行数值拟。
第二部分是国内外的研究现状。
从19世纪四十年代第一次提出粘弹性介质理论,粘性介质理论不断发展,到现在,已经逐步发展成熟。国内的相关研究也逐渐增多。
第三部分为正演模拟的过程。
流程图如图所示……
如图所示为粘声波方程,具体的推导过程在论文10-12页有详细的推到步骤。
由于涉及到了网格的选取,网格是用来实现介质模型离散化的一种常用手段,所以此处对网格进行了简单的介绍。由于计算量相同时,交错网格数值模比常规网格更精确,数值计算更稳定,收敛速度更快,压制数值频散更彻底。因此我使用交错网格进行模拟。
如图为推导出的差分系数计算矩阵与交错网格有限差分格式,推导过程详见论文13-16页。
对于震源的选择,此处,我选择的是雷克子波作为加载的震源,来模拟激发地震波。表达式如图所示。为了检测震源是否准确,我选用频率为30hz,做出雷克子波的特征图:
……
由于实际介质更加倾向于半无限空间,所以要进行边界条件的选择,边界条件分为衰减边界条件、吸收边界条件,本文应用的是是PML吸收边界条件,如图所示。然后为了检验边界条件的作用,我用
简单的单层介质模型进行了简单的测试,结果如图所示,从图中可以看出,当没有边界条件的时候,在界面处会产生明显的反射波。
有限差分法最大的局限性就是数值频散问题。造成数值频散现象的原因有很多,例如震源精度,采样间隔,波形畸变等因素的影响。可以用这三种方法来压制数值频散。……
然后我用选取不同阶数的正常声波常规网格有限差分算子进行简单的试算,如图所示:……
第四部分为模型的试算,为了探究地震波的传播规律,首先我建立了一个简单的正序三层分布模型,模型的相关参数以及震源位置如图所示,来模拟地震波在地下介质中的传播。……
然后我有建立了两个三层分布模型,一个含高速夹层,一个含低速夹层,震源位置如图所示。
这是他们在同一时刻的波场快照。……
之后,为了探究了粘声波的衰减规律,我建立了一个四层分布模型,模型参数如图所示,右图分别为正常声波跟粘声波沿测线方向的炮记录,为了对比明显,我选取第30道的声波记录跟粘声记录,放到一个表格中进行对比,如图所示,进行局部放大,可以看出粘声波的地震记录在深层的反射波的振幅明与普通声波的相比,有明显的衰减现象,同相轴相对教弱;在深层,粘声波的记录的波形变化比较严重。也就是说,粘性介质考虑到了地层对能量的衰减作用,反射波的能量比完全弹性的介质减弱很多,更加的符合实际的底层情况。
………
最后,再次谢谢各位老师的聆听,欢迎各位老师批评指正。
因篇幅问题不能全部显示,请点此查看更多更全内容