声波方程正演模拟.pptVIP

相对于上述几种方法，有限差分法是一种更为快速有效的方法。虽然其精度比不上有限元法，但因其具有计算速度快，占用内存较小的优点，在地震学界受到广泛的重视与应用。内容提纲*声波方程的有限差分法数值模拟04数值算法类型及其局限性03波动方程类型及其局限性02地震勘探基本原理01四、声波方程的有限差分法数值模拟*对于二维速度-深度模型，地下介质中地震波的传播规律可以近似地用声波方程描述：01是介质在点（x,z）处的纵波速度，02为描述速度位或者压力的波场，03为震源函数。04(4-1)05空间模型网格化（如图4-1所示）：*图4-1差分网格划分示意图网格间隔长度*表示(i,j)点k时刻的波场值时间采样步长时间二阶、空间二阶差分格式推导如下：*(4-3)在(i,j)点k时刻用Taylor展式展开：将(4-2)将上两式相加，略去高阶小量，整理得(i,j)点k时刻的二阶时间微商为：0102(4-4)同理可得(i,j)点k时刻的二阶空间微商分别为：*1(4-5)2(4-6)(4-7)是空间离散步长式中为时间采样率为介质速度的空间离散值这就实现了用网个点波场值的差商代替了偏微分方程的微商，将上三个式子代入(4-1)式中得：同理可得空间四阶精度的差分格式为：同理可得时间二阶、空间四阶精度的声波方程差分格式为：(4-8)前言前言前言前言前言前言前言前言*声波方程数值模拟

――地球物理场论

基础Ⅰ期末作业（1）任课教师：宋鹏内容提纲*地震勘探基本原理波动方程类型及其局限性数值算法类型及其局限性声波方程的有限差分法数值模拟一、地震勘探基本原理*＊▽▽▽▽▽xt同相轴为双曲线，即反射波的时距曲线为双曲线，反射波一个同相轴可带来一个地层的信息。实际地下介质非常复杂，所得到的炮集记录也包含更多的地下信息。实际的炮集记录见图1-1和1-2。图1-1陆上某区实际地震记录*图1-2海上某区实际地震记录*广义的地震反演即是从地震炮集记录出发，经过复杂的去噪、速度分析以及偏移成像处理等手段得到反映地下的地质结构的地震剖面。实际的地震剖面见图1-3和1-4。图1-3 陆上某区地震剖面*图1-4海上某区地震剖面*地震波场模拟即地震正演，是指已知模型结构，通过物理或数值计算的方法模拟该地质结构下的地震波的传播，最终合成地震记录，也可以认为其是野外数据采集过程的室内再现。物理模拟花费昂贵，人们一般采用比较经济的数值模拟技术。地震波场数值模拟是在给定数学模型（如弹性波方程，声波方程等）、震源和地下几何界面、物性参数（岩层密度、速度等）情况下，研究弹性波或声波的传播规律。内容提纲*地震勘探基本原理01波动方程类型及其局限性02数值算法类型及其局限性03声波方程的有限差分法数值模拟04二、波动方程类型及其局限性*一阶压力-速度方程组：二阶标量声波方程：声波方程：12波动方程类型及其局限性能够描述且只能描述纵波的传播规律，包括直达波、反射波、透射波、折射波等，但不能描述转换波传播规律。需要的已知条件包括：震源函数2）地层速度/密度3）边界条件2、弹性波方程：*21能够描述纵、横波的传播规律，包括直达波、反射波、透射波、折射波以及转换波等。需要的已知条件包括：震源函数2）地层速度或根据方程的类型需要提供的地层的其它弹性参数3）边界条件粘声波/弹性波方程前面讨论的是理想弹性介质，波在其中传播时，没有能量的损耗，介质中应力和应变关系严格遵循胡克定律（这种理想介质称虎克固体），但波在实际介质中传播时，是有能量损耗的，这就是所谓的弹性波吸收。波在传播过程中，实际介质的不同部位之间会出现某种摩擦力，称为内摩擦力或粘滞力。这种力导致机械能向其他形式能量转换，最终转化为热能消耗掉。但是到目前为止，在地震资料反演处理中应用最多的还是声波方程，弹性波以及粘弹性波方程的应用还只是停留在模拟层次上。在地震勘探中，地震波传播的实际介质是十分复杂的。在一定条件下，即震源作用时间短，作用力微小，地球介质可以看作完全弹性模型，但随着地震勘探技术的发展，勘探精度要求提高，面临复杂地质目标时，要求地震勘探采用更加符合实际的介质模型进行研究。粘弹性介质模型更符合实际。内容提纲*数值算法类型及其局限性声波方程的有限差分法数值模拟波动方程类型及其局限性地震勘探基本原理三、数值算法类型及其优缺点*地震波波动方程数值模拟方法主要