弹性波方程的建立弹性波方程的解析解数值模拟求解方法概览数值.PPT

弹性波方程及其数值模拟求解 弹性波方程的解析解 二阶三维声波方程 三、平面谐波解： 平面谐波解的简单形式 平面谐波解的具体形式 或 弹性波方程及其数值模拟求解 弹性波方程的解析解 平面谐波解的实质：将弹性波看做简谐波的叠加 简谐波 弹性波方程及其数值模拟求解 弹性波方程的解析解 平面谐波解的性质之一：若函数A，B均具有平面谐波解的形式，则对A，B的线性运算相当于对A，B中的简谐波分别对应进行线性运算，再将简谐波运算得到的结果叠加。最终得到的结果C仍然具有平面谐波解的形式。 弹性波方程及其数值模拟求解 弹性波方程的解析解 平面谐波解的性质之二：对于每一个简谐波来说，若有点A以及过A点的平面，该平面垂直于简谐波的传播方向，则A点和该平面上所有点的振动情况都相同。 弹性波方程及其数值模拟求解 弹性波方程的解析解 二阶三维声波方程 平面简谐波的性质之三：平面波以相速度V沿k的方向传播。 平面简谐波解的具体形式 角频率 角波数 相速度 弹性波方程及其数值模拟求解 弹性波方程的建立 弹性波方程的解析解 数值模拟求解方法概览 数值模拟算法的频散分析 数值模拟算法的稳定性条件 总结 弹性波方程及其数值模拟求解 数值模拟求解方法概览 解析解的局限性： 1. 大多数解析解包含了前提假设，实际情况下很难与前提假设相符，或仅在一定条件下可以满足前提假设。例如达朗贝尔解是在无限长弦的条件下得到的；平面波解隐含了空间无限广以及均匀介质的假设。 2.解析解中包含了复杂的数学计算，很多情况下难以求解，或者即使能够求解，过程和结果都太过复杂。例如要求得具体的达朗贝尔解，需要知道震源函数或其他定解条件；克西霍夫积分解需要求取曲面积分。 3.求取解析解要求信号是连续的函数，而实际生产中，地震信号是离散的。因而许多解析解无法应用，或做离散近似之后误差太大。 弹性波方程及其数值模拟求解 数值模拟求解方法概览 一阶声波方程 二阶声波方程 数值模拟方法：在已知一部分离散的波场值的前提下，通过求取离散的波场值的时间导数和空间导数，进而求取全部的离散波场值。 弹性波方程及其数值模拟求解 数值模拟求解方法概览 一阶声波方程 二阶声波方程 第一类：时间导数在时间域求取，空间导数在空间域求取： 有限差分法，有限单元法。 弹性波方程及其数值模拟求解 数值模拟求解方法概览 二阶声波方程 第二类：时间导数在频率域求取，空间导数在空间域求取： 频率-空间域正演模拟。 弹性波方程及其数值模拟求解 数值模拟求解方法概览 二阶声波方程 第三类：时间导数在时间域求取，空间导数在波数域求取： 伪谱法正演模拟。 本节参考 刘喜武 著《弹性波场论基础》 参见数学物理方程中的弦振动方程 参见 伊尔玛资 著 《地震资料分析》 参见 伊尔玛资 著 《地震资料分析》 参见 牟永光著 《三维复杂介质地震数值模拟》 参见 参见《accuracy_of_finite_difference_method_of _acoustic_wave_equations》 参见《任意偶数阶精度有限差分法数值模拟_刘洋》 作者：谷谷谷 弹性波方程及其数值模拟求解 弹性波方程的建立 弹性波方程的解析解 数值模拟求解方法概览 数值模拟算法的频散分析 数值模拟算法的稳定性条件 总结 弹性波方程及其数值模拟求解 弹性波方程是描述了弹性介质中质点振动的方程。 弹性波方程由几何方程，虎克定律，运动微分方程三套方程构建的。 弹性波方程的建立 位移矢量 应力张量 应变张量 弹性波方程及其数值模拟求解 几何方程：描述了位移和应变之间的关系。 弹性波方程的建立 弹性波方程及其数值模拟求解 虎克定律：描述了应变和应力之间的关系。 弹性波方程的建立 弹性波方程及其数值模拟求解 虎克定律：描述了应变和应力之间的关系。 弹性波方程的建立 弹性系数矩阵 各向同性介质的弹性矩阵 弹性波方程及其数值模拟求解 运动微分方程：描述了位移和应力之间的关系。 弹性波方程的建立 弹性波方程及其数值模拟求解 弹性波方程的建立 二阶位移方程的建立 几何方程的简化形式 微分算子矩阵 弹性波方程及其数值模拟求解 弹性波方程的建立 二阶位移方程的建立 微分算子矩阵 运动微分方程的简化形式 弹性波方程及其数值模拟求解 弹性波方程的建立 二阶位移方程的建立 弹性波方程及其数值模拟求解 弹性波方程的建立 二阶位移方程的建立 完成！ 二阶弹性波方程是描述了弹性介质中质点振动的幅度的方程。 略去外力项 弹性波方程及其数值模拟求解 弹性波方程的建立 二阶位移方程的优缺点 二阶位移方程的优点： 形式简单，只含三个方程和三个变量。 在均匀介质的前提下，容易对方程进行各种变形。 二阶位移方程的缺点： 无法利用能量通量密度矢量（Poynting矢量）求取波场传播方向。 无法