Abstract:Seismic traveltime tomography is a key technique for near-surface velocity model building. Specifically, the wave-equation tomography can correctly handle the first-order diffraction phenomena when the scale of velocity heterogeneity is comparable with the incident wavelength. Although wave-equation-based traveltime inversion can theoretically produce a velocity model with higher resolution compared with ray tomography, the computational cost is much higher which prohibits its implementation for industry-scaled 3-D problems. To address this problem, by introducing the random boundary condition and the discrete Fourier integration, we propose a practical 3-D waveequation first -arrival tomography using the monochromatic gradient inversion scheme. The advantages of the proposed method are: in terms of computational cost and storage space, the proposed method only needs very few memory spaces (i. e., compared with wave-equation forward modeling, only two additional monochromatic wavefields and one image volume are required), and meanwhile, the computational cost is decreased at least by 1/3 (i. e., compared with the wavefield reconstruction scheme, only two wavefield propagations are needed, and no special treatment for numerical boundaries). Besides, the introduction of random boundary conditions further simplifies the modeling algorithm which in turn makes it well suited for GPU (Graphics Processing Unit) acceleration. Theoretically, considering that the major area that may influence the first -order diffraction of seismic waves is usually referred to as the first Fresnel zone, which has a very similar shape to the monochromatic sensitivity kernel calculated with the dominant frequency. On one hand, the monochromatic sensitivity kernel within the first Fresnel zone will be coherently increased because the scattered wavefield caused by random boundaries becomes less spatial correlative. On the other hand, the multi-fold seismic acquisition system makes the sensitivity kernels outside the first Fresnel zone further canceled. Numerical example demonstrates that the smoothed monochromatic gradient contains sufficient low-wavenumber components which are adequate for background velocity inversion. Moreover, a 3-D synthetic data test demonstrates that a high-resolution near-surface velocity model can be inverted using wave-equation traveltime tomography with affordable computer resources. Compared with the traditional wave-equation tomography methods, the proposed method not only reduces the memory requirement remarkably but also significantly increases the computational efficiency. Therefore, it is very promising for dealing with industrial-sized 3-D inverse problems.

Traveltime computation using the upwind finite difference method with nonuniform grid spacing in a 3D undulating surface condition

Modeling Of Rayleigh-Wave Propagation In Three-Dimensional Media With Topography Using Staggered-Grid Finite Difference Scheme

Numerical Study of Turbidity Current over a Three-Dimensional Seafloor

An Efficient Immersed Free Surface Boundary Method for 3-D Scalar Seismic Waves Finite-Difference Modeling in Presence of Topography

A Three-dimensional Immersed Boundary Method for Accurate Simulation of Acoustic Wavefields with Complex Surface Topography

Application of grid successive subdivision algorithm in 3-D seismic ray tracing

Numerical Investigation of Rayleigh-wave Propagation on Topography Surface

A travel-time tomography method for improving horizontal resolution of high-frequency surface-wave exploration

3D frequency-domain modeling of marine controlled source electromagnetic responses with topography using finite volume method

A New Scheme for Fast Calculation of Seismic Traveltimes and Ray Paths in Heterogeneous Media

Simulate the elastic wavefields in media with an irregular surface topography based on staggered grid finite difference

3D Seismic-Wave Modeling with a Topographic Fluid–Solid Interface at the Sea Bottom by the Curvilinear-Grid Finite-Difference Method

Viscoelastic wave finite-difference modeling in the presence of topography with adaptive free-surface boundary condition

3-D Irregular Grids Finite Difference Method for Elastic Wave Propagation in Anisotropic Media

Rugged Topography Near-Surface Velocity Inversion Using 3-D Finite-Frequency Tomography

A Nonbalanced Staggered-Grid FDTD Scheme for the First-Order Elastic-Wave Extrapolation and Reverse-Time Migration

2-D Poroelastic Wave Modelling with a Topographic Free Surface by the Curvilinear Grid Finite-Difference Method

Incorporating topographic effects in surface wave tomography based on shortest-path ray tracing

Elastic reverse time migration using an efficient and accurate adaptive variable grid discretization method

3D Wave-Equation First-Arrival Tomography Using Monochromatic Traveltime Sensitivity Kernel

Seismic Propagation Simulation in Complex Media with Non-Rectangular Irregular-Grid Finite-Difference