Implementing the fast marching eikonal solver: spherical versus Cartesian coordinates

Spherical coordinates are a natural orthogonal system for describing wavefronts emanating from a point source. A regular grid distribution in the Cartesian‐coordinate system tends to undersample the wavefront description near the source (at the highest wavefront curvature) and oversample it away from the source. Spherical coordinates, in general, provide a more balanced grid distribution for characterizing point‐source wavefronts. Our numerical implementation confirms that the recently introduced fast marching algorithm is both a highly efficient and an unconditionally stable eikonal solver. However, its first‐order approximation of traveltime derivatives can induce relatively large traveltime errors for waves propagating in a diagonal direction with respect to the coordinate system. Examples, including the IFP Marmousi and the SEG/EAGE 3D salt‐dome models, show that a spherical‐coordinate implementation of the method results in far fewer errors in traveltime calculation than the conventional Cartesian‐coordinate implementation, and with practically no loss in computational advantages.