weixin_42262783的博客

对称模型，走时场，射线路径，波场等波现象均有良好的对称性，容易现象。对称模型，但走时场，射线路径等波现象均无对称性，比各向同性复杂很多。

在阅读文章的基础上，很多朋友对此研究很感兴趣，想进一步深入学习，在此，提供了敏感核分析的代码，下面是Thomsen参数化的代码，其他参数化的代码一并打包上传。（好吧，发现不能打包上传，谁有办法，评论区告诉我！（需要代码的私信我就好））说明：做该项研究当时用好几种办法推导过 结果都很相似 验证了我的推导是对的。

弹性波方程是通过应力（stress）应变（strain）关系（广义胡克定律，或本构关系）+运动微分方程（动量守恒线性方程 Linearized equation of momentum conservation Navier 方程或牛顿第二定律）导出，描述了波的所有动力学和运动学传播特性。而Christoffel 方程是各向异性弹性波动方程的平面波（单相波，对应相速度）解形式，可以用来独立描述三种类型的波的极化和相速度特征，即运动学（速度）和偏振（振幅以及方向）。3.几大方程的物理意义是什么？

本人做了复现，虽然是复现，但是也花了很久的时间。其中的一些问题，想了很久才解决。当然实现了以后，就觉得所有的问题原来都不是什么特别的难题。2.沿x轴和z轴群速度的计算，可以采用Shooting方法（Qian，2012），或者Fomel（2004）近似公式。1.因果律的物理含义，代表群速度（射线方向）要和走时更新的方向一致，其实代表了射线的传播方向，即群速度方向。上面这俩式子很重要，又能用来推导沿x轴和z轴群速度，还能用来做因果律的判定。直接上结果了，代码有需要的私聊我。

逐渐改进，从各向同性到VTI各向异性，再到TTI介质。模拟效果可以在走时计算的文章里看到，有对比。不懂的欢迎随时交流。

参考各向同性的程序，进一步修改为VTI介质，不懂的随时交流。

涉及为未发表成果，代码暂时不公开，欢迎私下交流讨论，互相学习。