无限弹性介质中的弹性波.ppt

相对于上述几种方法，有限差分法是一种更为快速有效的方法。虽然其精度比不上有限元法，但因其具有计算速度快，占用内存较小的优点，在地震学界受到广泛的重视与应用。第157页,共180页，星期六，2024年，5月 一、地震勘探基本原理

二、波动方程类型及其局限性

三、数值算法类型及其局限性

四、声波方程的有限差分法数值模拟内容提纲第158页,共180页，星期六，2024年，5月四、声波方程的有限差分法数值模拟对于二维速度-深度模型，地下介质中地震波的传播规律可以近似地用声波方程描述：是介质在点（x,z）处的纵波速度，为描述速度位或者压力的波场，为震源函数。(4-1)第159页,共180页，星期六，2024年，5月空间模型网格化（如图4-1所示）：图4-1差分网格划分示意图第160页,共180页，星期六，2024年，5月网格间隔长度，时间采样步长表示(i,j)点k时刻的波场值第161页,共180页，星期六，2024年，5月 时间二阶、空间二阶差分格式推导如下：(4-2)(4-3)第162页,共180页，星期六，2024年，5月将上两式相加，略去高阶小量，整理得(i,j)点k时刻的二阶时间微商为：(4-4)第163页,共180页，星期六，2024年，5月同理可得(i,j)点k时刻的二阶空间微商分别为：(4-5)(4-6)第164页,共180页，星期六，2024年，5月这就实现了用网个点波场值的差商代替了偏微分方程的微商，将上三个式子代入(4-1)式中得：(4-7)式中第165页,共180页，星期六，2024年，5月同理可得空间四阶精度的差分格式为：第166页,共180页，星期六，2024年，5月无限弹性介质中球面空腔产生的弹性波胀缩点震源引起的球面波取由中心向四周扩散的球面波的解，即为则位移场即第125页,共180页，星期六，2024年，5月无限弹性介质中球面空腔产生的弹性波胀缩点震源引起的球面波这样，在震源附近，质点的位移基本上重复震源强度函数变化规律，称近震源场，远离震源时，质点的位移是震源强度函数的导数。这说明球面波在其传播过程中波形逐渐改变，这也是区别平面波的一个重要特点。第126页,共180页，星期六，2024年，5月谐波波有着多种多样的形式，从波动方程的通解表达形式知波函数与激发子波或震源函数有关，但最简单而最重要的莫过于谐波（Harmonicwave），即简谐振动在弹性介质中的传播而形成的波。由Fourier级数理论级可知，任何复杂的振动都可以看成若干各不相同简谐振动的合成；因此任何复杂的弹性波也总可以看成若干个具有各不相同的频率，相位和振幅的谐波合成起来的。于是谐波在弹性介质中传播时，任何一点处的扰动都不仅随时间而且随空间作余弦（或正弦）型变化。第127页,共180页，星期六，2024年，5月谐波--平面谐波所谓的平面谐波是由若干个作简谐振动的质点发出的扰动所形成的几何包络面为平面的谐波。如下图所示。第128页,共180页，星期六，2024年，5月谐波--平面谐波沿x轴方向传播的平面谐波对于沿着x轴正向传播的平面谐波的扰动函数可以简单写成：（a）也可写成：第129页,共180页，星期六，2024年，5月谐波--平面谐波若设振动周期为T，圆频率而波数波传播速度为第130页,共180页，星期六，2024年，5月谐波--平面谐波图6?8振动图第131页,共180页，星期六，2024年，5月谐波--平面谐波图6?9波剖面图第132页,共180页，星期六，2024年，5月谐波--平面谐波相速度：指一定的相位移动的速度。群速度：指一定的振幅包络线移动的速度。相速度和群速度皆为波速，如果各个简谐波相速度相同，则群速度与相速度相等；若各个简谐波相速度不同，群速度与相速度不等，此时即发生频散。频散：由各种不同谐波成分组成的波，虽然受同一起始扰动下，但各自以不同的速度传播，并且起始扰动的形状在传播中将产生变化。扰动经传播以后将扩展成为一更长的波列，这种现象我们称之为频散。第133页,共180页，星期六，2024年，5月谐波--球面简谐纵波所谓球面谐波是由若干作简谐振动质点发出的扰动所形成的波阵面为球面的谐波，如下图所示，对于球对称情况，则为球面简谐纵波。波函数（扰动函数）可以写成：写成复指数形式设为位移场的标量位，则：当r波长时，第一项起主要作用，称为远震源场。第134页,共180页，星期六，2024年，5月广义胡克定律第135页,