The impressive features of the fast marching method, stability and efficiency, are achieved at the expense of accuracy. The approach is based on a first-order approximation of the traveltime derivatives with respect to position. This low-order approximation can result in relatively large traveltime errors when the wavefront curvature is large and the wavefront propagation is diagonal to the grid orientation. In Cartesian coordinates, spherical wavefronts can spend a long time satisfying these two destructive conditions. As a result, implementation of the fast marching method in the Cartesian coordinate, unless the gird is very fine, can result in relatively large errors. Recently, () suggested a spherical coordinate implementation of the fast marching method. In spherical coordinates, wavefronts emanating from a source typically spend less time traveling diagonally with respect to the coordinate system than the case in Cartesian coordinates. As a result, spherical coordinate implementation achieves more accurate traveltimes even for complex models like the Marmousi model, without compromising stability or efficiency.
The SEG/EAGE 3-D velocity model was used to generate 3-D prestack synthetic datasets. These datasets are commonly used to test migration algorithms, especially 3-D prestack Kirchhoff ones. At the center of a 3-D Kirchhoff-based migration is 3-D traveltime computation, and the quality of the migration depends heavily on the accuracy of the traveltime calculation. Thus, the 3-D slat-dome velocity model is useful for testing 3-D traveltime calculation methods, including the fast marching method, which is the main focus here.
In this paper, I will highlight the benefits of executing the fast marching method of traveltime calculation in spherical coordinates. To show these benefits, I will use the SEG/EAGE 3-D salt-dome velocity model. The complexity of the 3-D salt structure provides us with a good measure of the accuracy of spherical coordinate implementation. We will also look into the role that head-waves and other low, but fast, energy arrivals play in traveltime maps.