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.

A fast amplitude preserving three-parameter 3D parabolic Radon transform and its application on multiple attenuation

A 3D constrained hybrid linear-pseudo hyperbolic Radon transform

High-precision Seismic Data Reconstruction with Multi-Domain Sparsity Constraints Based on Curvelet and High-Resolution Radon Transforms

Seismic Data Two-Step Recovery Approach Combining Sparsity-Promoting and Hyperbolic Radon Transform Methods

Amplitude-Preserving Nonlinear Adaptive Multiple Attenuation Using the High-Order Sparse Radon Transform

An Efficient Amplitude-Preserving Radon Transform With Frequency-Dependent Curvature for Multiple Attenuation

Amplitude-preserving Iterative Deblending of Simultaneous Source Seismic Data Using High-Order Radon Transform

Fast High-Resolution Hyperbolic Radon Transform

2D and 3D Prestack Seismic Data Regularization Using an Accelerated Sparse Time-Invariant Radon Transform

An interpolating scaling functions method with low-storage five-stage fourth-order explicit Runge-Kutta schemes for 3D ground penetrating radar simulation

Three-Dimensional Reverse Time Migration of Ground-Penetrating Radar Signals

3D Seismic Data Reconstruction based on Weighted Fast Iterative Shrinkage Thresholding algorithm

Three-dimensional passive seismic waveform imaging around the SAFOD site, California, using the generalized Radon transform

Automatic Velocity Analysis Using High-Resolution Hyperbolic Radon Transform

Robust multi-dimensional reconstruction via Group Sparsity with Radon operators

Improving Receiver Function Imaging with High-Resolution Radon Transform

Plane-wave Orthogonal Polynomial Transform for Amplitude-Preserving Noise Attenuation

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

3-D Tunnel Space Elastic Reverse Time Migration in Cylindrical Coordinates Using Reconstructed Time-Domain Wavefields

An efficient 3D wave-equation pre-stack time migration for high-density and wide-azimuth data

Seismic data extrapolation based on multi-scale dynamic time warping