弹性波论文.docVIP

声波方程数值模拟实验报告 实验题目：声波方程数值模拟 实验目的： （1）通过模拟人工合成的地震记录，体会地震勘探基本原理和方法，验证地震波传播能量波形变化趋势。 英特尔? 酷睿? 2、声波方程的有限差分法数值模拟 对于二维速度-深度模型，地下介质中地震波的传播规律可以近似地用声波的方程描： （2） V（x，z）是介质在点（x , z）处的纵波速度，u为描述速度位或者压力的波场，s（t）为震源函数。 为求上式的数值解，必须将此式离散化，用差商代替微商。把空间模型网格化，如下图： 对空间微分采用四阶精度差分格式，时间微分采用二阶精度差分格式，得： （3） 式中为介质速度的空间离散值，是空间离散步长，为时间离散步长，s（t）为震源函数，确定震源位置。 在实际应用中，速度、时间步长、空间步长需满足一定的条件，具体情况如下： 1、稳定性条件： （4） 这里表示的是地下介质的最大波速；若地下介质网格间隔、最小速度、及时间采样间隔不符合(4-8)式时，第推求解（4-7）式，波场值会出现误差（高阶小量）