A Level Set Based Eulerian Method for Paraxial Multivalued Traveltimes

We apply the level-set methodology to compute multivalued solutions of the paraxial eikonal equation in both isotropic and anisotropic metrics. This paraxial equation is obtained from 2-D stationary eikonal equations by using one of the spatial directions as the arti cial evolution direction. The advection velocity eld used to move level sets is obtained by the method of characteristics; therefore the motion of level sets is de ned in a phase space, and the zero level set yields the location of bicharacteristic strips in the reduced phase space. The multivalued traveltime is obtained from solving another advection equation with a source term. The complexity of the algorithm is O(N LogN) in the worst case and O(N as the typical Lagrangian ray-tracing, where N is the number of the sampling points along one of the spatial directions. Numerical experiments including the well-known Marmousi synthetic model illustrate the accuracy and the eciency of the Eulerian method