Wave tracing: Ray tracing for the propagation of band-limited signals: Part 2 — Applications

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE385-VE393 ◽  
Author(s):  
John K. Washbourne ◽  
Kenneth P. Bube ◽  
Pedro Carillo ◽  
Carl Addington

Modeling seismic propagation is critically important to our work; unfortunately, we often must trade simulation accuracy for reduced computational expense. We present a new seismic-modeling method that is as simple and computationally efficient as Snell’s law ray tracing but provides propagation paths and arrival times more consistent with finite-bandwidth data. We refer to this modeling method as wave tracing and apply it to nonlinear traveltime tomography and depth imaging. By replacing Snell’s law ray tracing with wave tracing, we get better ray coverage, more robust and faster ray bending (fewer iterations), and a much more robust and faster algorithm for nonlinear tomography (fewer iterations, too). A very significant benefit is increased stability and robustness of tomographic inversion with respect to small changes in model parameterization and regularization. A related benefit is the increased stability of depth images with respect to small changes in velocity, which can increase confidence in interpretation. The velocity models that result from wave tracing match picked arrival times in band-limited data better and generate improved depth images. These advantages of wave tracing relative to conventional Snell’s law ray tracing have been tested on both synthetic and real data examples for crosswell seismic geometry.

2004 ◽  
Vol 04 (02) ◽  
pp. 325-339
Author(s):  
KAZUHIRO KOYAMA ◽  
YOSHIAKI TOMIZAWA ◽  
MINORU OKADA

In this paper we propose two improvements to the beam tracing method, one is to reduce the final redundant-drawing of pixel fragments into the frame buffer, and the other is to use Snell's law instead of the tangent law. It is well known that the ray tracing scheme is one of the most effective methods for high quality CG image synthesis. The beam tracing method was introduced because traditional ray tracing algorithms have a significant calculation cost. In the improved method, the projection screen is recursively divided into non-overlapping coherent region segments based on the coherence of the ray-trace-tree structure of a given scene, and an image can be synthesized after the segmentation process by tracing one or several rays in each region segment. We introduced Snell's law for the refractions of the rays. Therefore, we can reduce the computational cost of ray tracing, eliminate the extra calculations caused by overlapping polygons on a screen, and increase the quality of CG images generated by our proposed method.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE377-VE384 ◽  
Author(s):  
Kenneth P. Bube ◽  
John K. Washbourne

Many seismic imaging techniques require computing traveltimes and travel paths. Methods to compute raypaths are usually based on high-frequency approximations. In situations such as head waves, these raypaths minimize traveltime but are not paths along which most of the energy travels. We have developed a new approach to computing raypaths, using a modification of ray bending that we call wave tracing; it computes raypaths and traveltimes that are more consistent with the paths and times for the band-limited signals in real data than the paths and times obtained using high-frequency approximations. Wave tracing shortens the raypath while keeping the raypath within the Fresnel zone for a characteristic frequency of the signal.


Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. S251-S259 ◽  
Author(s):  
Claudia Vanelle ◽  
Miriam Spinner ◽  
Thomas Hertweck ◽  
Christoph Jäger ◽  
Dirk Gajewski

True-amplitude Kirchhoff migration (TAKM) is an important tool in seismic-reflection imaging. In addition to a structural image, it leads to reflectivity maps of the subsurface. TAKM is carried out in terms of a weighted diffraction stack where the weight functions are computed with dynamic ray tracing (DRT) in addition to the diffraction traveltimes. DRT, however, is time-consuming and imposes restrictions on the velocity models, which are not always acceptable. An alternative approach to TAKM is proposed in which the weight functions are directly determined from the diffraction traveltimes. Because other methods exist for the generation of traveltimes, this approach is not limited by the requirements for DRT. Applications to a complex synthetic model and real data demonstrate that the image quality and accuracy of the reconstructed amplitudes are equivalent to those obtained from TAKM with DRT-generated weight functions.


Author(s):  
Julie Bentley ◽  
Craig Olson

Geophysics ◽  
1993 ◽  
Vol 58 (7) ◽  
pp. 987-996 ◽  
Author(s):  
Robert Fischer ◽  
Jonathan M. Lees

A technique for improving the efficiency of shortest path ray tracing (SPR) is presented. We analyze situations where SPR fails and provide quantitative measures to assess the performance of SPR ray tracing with varying numbers of nodes. Our improvements include perturbing the ray at interfaces according to Snell’s Law, and a method to find correct rays efficiently in regions of low velocity contrast. This approach allows the investigator to use fewer nodes in the calculation, thereby increasing the computational efficiency. In two‐dimensional (2-D) cross‐borehole experiments we find that with our improvements, we need only use 2/3 as many nodes, saving up to 60 percent in time. Savings should be even greater in three dimensions. These improvements make SPR more attractive for tomographic applications in three dimensions.


Geophysics ◽  
2007 ◽  
Vol 72 (4) ◽  
pp. T27-T36 ◽  
Author(s):  
Chao-ying Bai ◽  
Stewart Greenhalgh ◽  
Bing Zhou

We present an accurate 3D ray-tracing algorithm based on a modified (more flexible and economical) shortest-path method (SPM). Unlike the regular SPM in the 3D case, which uses only primary nodes at the corners of each cell and whose accuracy depends on actual cell size, the new method can work with much larger cell sizes by introducing secondary nodes along all bounding surfaces of the cell. This increases the ray angular coverage and permits detailed specification of the velocity field. The modified SPM simultaneously calculates first-arrival times and gradually locates the related raypaths on all grid nodes as the wave field evolves. Its advantages over the regular SPM are its ability to handle high-contrast velocity models more easily, lower memory requirements and less CPU time, and the capability to calculate a relatively large 3D model without losing accuracy. The maximum relative error bound in the computed traveltimes of the modified SPM is established for a uniform velocity field, which may be considered an upper error bound for the whole model in real problems. The modified method in this study is compared with the regular SPM theoretically and on two specific velocity models. The Marmousi model is used to further test the performance of the new approach for both accuracy and flexibility in a complex velocity field. The study shows that the modified SPM is preferable to regular SPM for real 3D problems.


Sign in / Sign up

Export Citation Format

Share Document