ALBERT

All Library Books, journals and Electronic Records Telegrafenberg

feed icon rss

Your email was sent successfully. Check your inbox.

An error occurred while sending the email. Please try again.

Proceed reservation?

Export
  • 1
    Electronic Resource
    Electronic Resource
    PO Box 1354, 9600 Garsington Road , Oxford OX4 2XG , UK . : Blackwell Science Ltd
    Geophysical prospecting 53 (2005), S. 0 
    ISSN: 1365-2478
    Source: Blackwell Publishing Journal Backfiles 1879-2005
    Topics: Geosciences , Physics
    Notes: For multivalued traveltime computation on dense grids, we propose a wavefront-orientated ray-tracing (WRT) technique. At the source, we start with a few rays which are propagated stepwise through a smooth two-dimensional (2D) velocity model. The ray field is examined at wavefronts and a new ray might be inserted between two adjacent rays if one of the following criteria is satisfied: (1) the distance between the two rays is larger than a predefined threshold; (2) the difference in wavefront curvature between the rays is larger than a predefined threshold; (3) the adjacent rays intersect. The last two criteria may lead to oversampling by rays in caustic regions. To avoid this oversampling, we do not insert a ray if the distance between adjacent rays is smaller than a predefined threshold. We insert the new ray by tracing it from the source. This approach leads to an improved accuracy compared with the insertion of a new ray by interpolation, which is the method usually applied in wavefront construction. The traveltimes computed along the rays are used for the estimation of traveltimes on a rectangular grid. This estimation is carried out within a region bounded by adjacent wavefronts and rays. As for the insertion criterion, we consider the wavefront curvature and extrapolate the traveltimes, up to the second order, from the intersection points between rays and wavefronts to a gridpoint. The extrapolated values are weighted with respect to the distances to wavefronts and rays. Because dynamic ray tracing is not applied, we approximate the wavefront curvature at a given point using the slowness vector at this point and an adjacent point on the same wavefront. The efficiency of the WRT technique is strongly dependent on the input parameters which control the wavefront and ray densities. On the basis of traveltimes computed in a smoothed Marmousi model, we analyse these dependences and suggest some rules for a correct choice of input parameters. With suitable input parameters, the WRT technique allows an accurate traveltime computation using a small number of rays and wavefronts.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 2
    Electronic Resource
    Electronic Resource
    Oxford, UK : Blackwell Science Ltd
    Geophysical prospecting 50 (2002), S. 0 
    ISSN: 1365-2478
    Source: Blackwell Publishing Journal Backfiles 1879-2005
    Topics: Geosciences , Physics
    Notes: To carry out a 3D prestack migration of the Kirchhoff type is still a task of enormous computational effort. Its efficiency can be significantly enhanced by employing a fast traveltime interpolation algorithm. High accuracy can be achieved if second-order spatial derivatives of traveltimes are included in order to account for the curvature of the wavefront. We suggest a hyperbolic traveltime interpolation scheme that permits the determination of the hyperbolic coefficients directly from traveltimes sampled on a coarse grid, thus reducing the requirements in data storage. This approach is closely related to the paraxial ray approximation and corresponds to an extension of the well-known 〈inlineGraphic alt="inline image" href="urn:x-wiley:00168025:GPR0285:GPR_0285_mu1" location="equation/GPR_0285_mu1.gif"/〉 method to arbitrary heterogeneous and complex media in 3D. Application to various velocity models, including a 3D version of the Marmousi model, confirms the superiority of our method over the popular trilinear interpolation. This is especially true for regions with strong curvature of the local wavefront. In contrast to trilinear interpolation, our method also provides the possibility of interpolating source positions, and it is 5–6 times faster than the calculation of traveltime tables using a fast finite-difference eikonal solver.
    Type of Medium: Electronic Resource
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
  • 3
    Publication Date: 2021-08-11
    Description: Full waveform inversion or adjoint tomography has routinely been performed to image the internal structure of the Earth at high resolution. This is typically done using the Fréchet kernels and the approximate Hessian or the approximate inverse Hessian because of the high-computational cost of computing and storing the full Hessian. Alternatively, the full Hessian kernels can be used to improve inversion resolutions and convergence rates, as well as possibly to mitigate interparameter trade-offs. The storage requirements of the full Hessian kernel calculations can be reduced by compression methods, but often at a price of accuracy depending on the compression factor. Here, we present open-source codes to compute both Fréchet and full Hessian kernels on the fly in the computer random access memory (RAM) through simultaneously solving four wave equations, which we call Quad Spectral-Element Method (QuadSEM). By recomputing two forward fields at the same time that two adjoint fields are calculated during the adjoint simulation, QuadSEM constructs the full Hessian kernels using the exact forward and adjoint fields. In addition, we also implement an alternative approach based on the classical wavefield storage method (WSM), which stores forward wavefields every kth (k≥1) timestep during the forward simulation and reads required fields back into memory during the adjoint simulation for kernel construction. Both Fréchet and full Hessian kernels can be computed simultaneously through the QuadSEM or the WSM code, only doubling the computational cost compared with the computation of Fréchet kernels alone. Compared with WSM, QuadSEM can reduce the disk space and input/output cost by three orders of magnitude in the presented examples that use 15,000 timesteps. Numerical examples are presented to demonstrate the functionality of the methods, and the computer codes are provided with this contribution.
    Print ISSN: 0895-0695
    Electronic ISSN: 1938-2057
    Topics: Geosciences
    Location Call Number Expected Availability
    BibTip Others were also interested in ...
Close ⊗
This website uses cookies and the analysis tool Matomo. More information can be found here...