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.

Anisotropic media hybrid grid tomographic inversion method based on angle domain imaging gathers

An efficient least-squares reverse time migration in image domain

A Guided Filter-Based 3D Hybrid Variational Optical Flow for Accurate Tomographic PIV Measurements

Reverse-time Migration Based on Generalized Synthesis of Seismic Data in F Requency-Space Domain

Study of Travel Time Tomography by Iterative Tikhonov Regularization Inversion

Extracting angle-domain common-image gathers in VTI media using ansiotropic-Helmholtz P/S wave-mode decomposition

Study on Inversion Methods for Travel Time Tomography Based on Generalized Inverse Theory

Travel time tomography by ray tracing using the fast sweeping method

Wave-equation First-Arrival Tomography for 3-D Near-Surface Velocity Inversion

Adjoint-state Traveltime Tomography for Azimuthally Anisotropic Media in Spherical Coordinates

INVERSION OF P-WAVE ANISOTROPIC PARAMETERS FROM CROSSWELL SEISMIC DATA

Elastic reverse time migration based on nested triangular mesh in 2D isotropic media

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

Slope Tomography Dynamically Weighted According to the Locations of the Reflection Points

Two-Dimensional Joint Inversion of Seismic Velocity and Electrical Resistivity Using Seismic Travel Times and Full Channel Electrical Measurements Based on Alternating Cross-Gradient Structural Constraint

The Analysis and the Solution of the Grid Size Problem to Travel-Time Tomography

Multigrid-Based Inversion for Volumetric Radar Imaging With Asteroid Interior Reconstruction as a Potential Application

Direct inversion of surface wave dispersion data with multiple-grid parametrizations and its application to a dense array in Chao Lake, eastern China

Multiple-grid Model Parametrization for Seismic Tomography with Application to the San Jacinto Fault Zone

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

Anisotropic Wave Separation Elastic Reverse Time Migration Based on the Pseudo-Decoupled Wave Equations in VTI Media