Abstract:Global orogenic belts and the regions of middle and western China have strongly varied topography. Crustal structure imaging in such condition is of great interest in both the mineral resource exploration and the geodynamics study. The probing depth of seismic reflection Waves in deep seismic sounding can reach up to the whole crust and even the top of the upper mantle of the earth, which is an effective supplement to the first-arrival refraction waves traveling in shallower depth. Current methods to deal with irregular surface usually take the step shape approximation or use a low velocity layer to overlay the surface. However, these methods may cause accuracy loss or even imaging distortion. This paper aims to introduce a grid-based method for traveltime tomography using the reflection traveltimes, to invert for velocity structure of the whole crust with an irregular surface. To deal with the irregular surface problem, this study is based on an irregular surface flattening scheme. The irregular surface flattening, which involves the transformation between the curvilinear and the Cartesian coordinate systems, can describe rugged terrain with no accuracy loss. Wave-front traveltime is calculated using the fast-sweeping method (FSM) to solve the topography-dependent Eikonal equation, and a multistage technique is then applied to calculate the reflection waves by reinitializing FSM in the incident layer. Ray paths are found by following the steepest traveltime gradient from the receiver to the interface and then to the source. Using the back-projection algorithm, the slowness perturbations are obtained, which have a relationship with the already existing traveltime. The final velocity model is obtained after a few iterations during which the slowness is updated, until the fit of the reflection traveltimes is satisfactory. We take four numerical examples to verify our method. First, the perfect match of traveltimes and ray paths obtained by our method and the shortest path method verifies that our forward modeling of the reflection traveltime tomography is accurate and reliable. Second, we test the dependency of the traveltime inversion on the initial velocity model. All four initial models converge to the true model within fairly good accuracy, which shows that our method is robust and effective with regard to the choice of the initial models. Third, our method is applied to a relatively complex model with strong variation of topography, and the final inverted velocity model gives correct ray paths of the reflections, very low traveltime residual, and high resolution velocity structure. Finally, both the full recovery of checkerboard patterns and the high-precision results inversed with the reflection traveltimes with Gaussian noise added show robustness and effectiveness of the proposed method. In summary, this paper addresses the problem of the irregular surface in seismic traveltime tomography. Based on the body-conforming grid and coordinate transformations, the wave-fronts of the reflection waves beneath the irregular surface are traced with high accuracy using the FWM and the multistage technique. Consequently, the first-arrival traveltime tomography is extended to the reflection traveltime tomography to image the deep structure of the crust in regions with complex surface topography. Numerical experiments show its reliability and effectiveness. This method can thus have a great potential application in mountain regions and basins, where high accuracy and high resolution seismic imaging is crucial to unravel their structure and reveal the tectonic evolution.

Wide‐spectrum Reconstruction of a Velocity Model Based on the Wave Equation Reflection Inversion Method and Its Application

Reflection Waveform Inversion Based on Full-Band Seismic Data Reconstruction for Salt Structure Inversion

Reflection Multi‐scale Envelope Inversion

Multiscale direct-envelope inversion using wavenumber filtering based on reflection tomography for background velocity recovery of large-scale, strong-contrast media

Waveform Imaging Based on Linear Forward Representations for Scalar Wave Seismic Data

Data-domain wave equation reflection traveltime tomography

Reflected wave least squares reverse time migration with angle illumination compensation

Background Velocity Model Building by Wave-Equation Reflection Traveltime Tomography

Characteristic reflector-based wave equation reflection traveltime inversion

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

Velocity modeling for joint tomography of data domain first-arrival traveltime and imaging domain reflection traveltime

Elastic Model Low- to Intermediate-Wavenumber Inversion Using Reflection Traveltime and Waveform of Multicomponent Seismic Data

Wave-equation Traveltime Slope Inversion by Combining Finite Difference and Cross-Correlation Methods

High-precision reflection traveltime tomography for velocity structure with an irregular surface

Full Waveform Inversion with Combined Misfit Functions and Application in Land Seismic Data

Envelope-Based Sparse-Constrained Deconvolution for Velocity Model Building

Converted-wave Reflection Traveltime Inversion with Free-Surface Multiples for Ocean-Bottom-node Data

Elastic Wave-equation Reflection Traveltime Inversion Using Dynamic Warping and Wave Mode Decomposition

Correlation-Based Reflection Full-Waveform Inversion

Regularized Seismic Velocity Inversion Based on Inverse Scattering Theory

Acoustic Wave‐equation Traveltime and Waveform Inversion of Crosshole Seismic Data