声波方程正演模拟省公开课一等奖全国示范课微课金奖课件.pptx

1声波方程数值模拟

――地球物理场论

基础Ⅰ期末作业（1）任课教师：宋鹏第1页

2一、地震勘探基本原理

二、波动方程类型及其不足

三、数值算法类型及其不足

四、声波方程有限差分法数值模拟内容提纲第2页

3一、地震勘探基本原理＊▽▽▽▽▽xt第3页

4 同相轴为双曲线，即反射波时距曲线为双曲线，反射波一个同相轴可带来一个地层信息。实际地下介质非常复杂，所得到炮集统计也包含更多地下信息。实际炮集统计见图1-1和1-2。第4页

5图1-1陆上某区实际地震统计第5页

6图1-2海上某区实际地震统计第6页

7广义地震反演即是从地震炮集统计出发，经过复杂去噪、速度分析以及偏移成像处理等伎俩得到反应地下地质结构地震剖面。实际地震剖面见图1-3和1-4。第7页

8图1-3 陆上某区地震剖面第8页

9图1-4海上某区地震剖面第9页

10地震波场模拟即地震正演，是指已知模型结构，经过物理或数值计算方法模拟该地质结构下地震波传输，最终合成地震统计，也能够认为其是野外数据采集过程室内再现。物理模拟花费昂贵，人们普通采取比较经济数值模拟技术。地震波场数值模拟是在给定数学模型（如弹性波方程，声波方程等）、震源和地下几何界面、物性参数（岩层密度、速度等）情况下，研究弹性波或声波传输规律。第10页

11 一、地震勘探基本原理

二、波动方程类型及其不足

三、数值算法类型及其不足

四、声波方程有限差分法数值模拟内容提纲第11页

12二、波动方程类型及其不足1、声波方程：一阶压力-速度方程组：二阶标量声波方程：第12页

13二、波动方程类型及其不足能够描述且只能描述纵波传输规律，包含直达波、反射波、透射波、折射波等，但不能描述转换波传输规律。需要已知条件包含：1）震源函数2）地层速度/密度3）边界条件第13页

142、弹性波方程：第14页

15能够描述纵、横波传输规律，包含直达波、反射波、透射波、折射波以及转换波等。需要已知条件包含：1）震源函数2）地层速度或依据方程类型需要提供地层其它弹性参数3）边界条件第15页

163、粘声波/弹性波方程 前面讨论是理想弹性介质，波在其中传输时，没有能量损耗，介质中应力和应变关系严格遵照胡克定律（这种理想介质称虎克固体），但波在实际介质中传输时，是有能量损耗，这就是所谓弹性波吸收。波在传输过程中，实际介质不一样部位之间会出现某种摩擦力，称为内摩擦力或粘滞力。这种力造成机械能向其它形式能量转换，最终转化为热能消耗掉。第16页

17 在地震勘探中，地震波传输实际介质是十分复杂。在一定条件下，即震源作用时间短，作用力微小，地球介质能够看作完全弹性模型，但伴随地震勘探技术发展，勘探精度要求提升，面临复杂地质目标时，要求地震勘探采取愈加符合实际介质模型进行研究。粘弹性介质模型更符合实际。 不过到当前为止，在地震资料反演处理中应用最多还是声波方程，弹性波以及粘弹性波方程应用还只是停留在模拟层次上。第17页

18 一、地震勘探基本原理

二、波动方程类型及其不足

三、数值算法类型及其不足

四、声波方程有限差分法数值模拟内容提纲第18页

19三、数值算法类型及其优缺点地震波波动方程数值模拟方法主要包含克希霍夫积分法、傅里叶变换法、有限元法和有限差分法等。克希霍夫积分法引入射线追踪过程，本质上是波动方程积分解一个数值计算，在某种程度上相当于绕射叠加。该方法计算速度较快，但因为射线追踪中存在着诸如焦散、多重路径等问题，故其普通只能适合于较简单模型，难以模拟复杂地层波场信息。第19页

20傅里叶变换法是利用空间全部信息对波场函数进行三角函数插值，能愈加准确地模拟地震波传输规律，同时，利用快速傅里叶变换(FFT)进行计算，还能够提升运算效率，其主要优点是精度高，占用内存小，但缺点是计算速度较慢，对模型适用性差，尤其是不适应于速度横向改变猛烈模型。第20页

21波动方程有限元法做法是：将变分法用于单元分析，得到单元矩阵，然后将单元矩阵总体求和得到总体矩阵，最终求解总体矩阵得到波动方程数值解；其主要优点是理论上可适宜于任意地质体形态模型，确保复杂地层形态模拟逼真性，到达很高计算精度，但有限元法主要问题是占用内存和运算量均较大，不适合用于大规模模拟，所以该方法在地震波勘探中还未得到广泛地应用。第21页

22相对于上述几个方法，有限差分法是一个更为快速有效方法。即使其精度比不上有限元法，但因其含有计算速度快，占用内存较小优点，在地震学界受到广泛重视与应用。