Finite-difference calculation of traveltimes in three dimensions
نویسنده
چکیده
The traveltimes of first arriving seismic rays through most velocity structures can be computed rapidly on a three-dimensional numerical grid by finite-difference extrapolation. Head waves are properly treated and shadow zones are filled by the appropriate diffractions. Differences of less than 0.11 percent are found between the results of this technique and ray tracing for a complex but smooth model. This scheme has proven useful for earthquake location and shows promise as an inexpensive, well-behaved substitute for ray tracing in forward-modeling and Kirchhoff inversion applications.
منابع مشابه
Shortest-path calculation of first arrival traveltimes by expanding wavefronts
A new approach to computing traveltimes and ray paths by solving the shortest path problem is presented. The technique is based on a partitioning of the shortest path optimization problem into smaller problems. We recursively evaluate the solution on expanding wavefronts instead of finding the global shortest paths from the source. To solve the local minimization, we apply a modified version of...
متن کاملA Fast Adaptive Vortex Method in Three Dimensions *
The method of local corrections (MLC) developed by Anderson for two spatial dimensions is a particle-particle panicle-mesh method. in which the calculation of the velocity field induced by a collection of vortices is split into two parts: (i) a finite difference velocity field calculation using a fast Poisson solver, the results of which are used to represent the velocity field induced by vorti...
متن کاملAn adaptive finite-difference method for traveltimes and amplitudes
The point-source traveltime field has an upwind singularity at the source point. Consequently, all formally high-order, finite-difference eikonal solvers exhibit firstorder convergence and relatively large errors. Adaptive upwind finite-difference methods based on high-order Weighted Essentially NonOscillatory (WENO) RungeKutta difference schemes for the paraxial eikonal equation overcome this ...
متن کاملAn efficient finite difference time domain algorithm for band structure calculations of Phononic crystal
In this paper, a new algorithm for studying elastic wave propagation in the phononic crystals is presented. At first, the displacement-based forms of elastic wave equations are derived and then the forms are discretized using finite difference method. So the new algorithm is called the displacement-based finite difference time domain (DBFDTD). Three numerical examples are computed with this met...
متن کاملAn Amplitude and Traveltime Calculation Using a Higher-order Parabolic Equation
A higher-order parabolic equation is used to compute the traveltime (phase) and the amplitude in constant density acoustic media. This approach is in the frequency domain, thereby avoiding the high frequency approximation inherent in the Eikonal equation. Intrinsic attenuation can be naturally incorporated into the calculation. The error at large angles of propagation caused by the expansion of...
متن کامل