地球物理MATLAB源代码PFAST:用于传输和反射地震旅行时间层析成像的高性能软件包概括
PFAST是基于并行地震层析成像的并行快速扫描方法的缩写。档为用户提供了快速入门使用该程序的基本教程。PFAST支持在2-D/3-D任意异质各向同性模型中获得传播时间,执行常规的初到行时间层析成像,以及在2-D中进行反射到达断层扫描(版本1.1),并使用第一个和反射到达的二维图像进行联合层析成像(版本1.1)。
PFAST的基本算法基于以下三个要素:
(a) 快速扫描法(FSM),一种基于网格的Eikonal方程求解器;
(b) 惠更斯原理,使用FSM计算反射传播时间;
(c) 伴随方法,用于获得非线性目标函数的梯度,而无需费时地评估Fréchet导数矩阵。有关该理论的详细信息,请参见发表在《国际地球物理杂志》上的论文(Huang和Bellefleur,2012年)及其在具有复杂近地表喀斯特湖的北极多年冻土地区的应用(Huang和Bellefleur,2011年)。目前,PFAST支持2-D各向同性模型中的联合层析成像和3-D各向同性模型。