各向异性TTI介质快速扫描法（Fast Sweeping Method，FSM）走时计算

该方法原创作者为Waheed（2015），文章在地球物理权威期刊Geophysics上发表，向原创致敬。本人做了复现，虽然是复现，但是也花了不少的时间。其中的一些问题，想了很久才解决。当然所有问题都解决以后回头再看，就觉得所有的问题原来都不是什么特别的难题。

方法原理极力推荐阅读参考其原创文章，在这里强调几个重点：

1.因果律的物理含义，代表群速度（射线方向）要和走时更新的方向一致，其实代表了射线的传播方向，即群速度方向。

2.沿x轴和z轴群速度的计算，可以采用Shooting方法（Qian，2012），或者Fomel（2004）近似公式。但更准确的是直接用Hamilton形式的程函方程推导。

3.TTI椭圆各向异性介质程函方程对应的Hamilton可统一表达为

令，,则有

上面这俩式子很重要，既能用来推导沿x轴和z轴群速度，又能用来做因果律的判定。

重点强调完了，直接上结果，有需要代码的直接私聊我，本人非常愿意相互学习共享基础性研究的成果。