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.

Velocity inversion method with reflection travel-time based on Gaussian packet

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

Regularized Seismic Velocity Inversion Based on Inverse Scattering Theory

Preconditioned Transmission + Reflection Joint Traveltime Tomography with Adjoint-State Method for Subsurface Velocity Model Building

On the method of multidimensional Q-factor and velocity inversion with VSP data

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

First arrival traveltime tomography based on an improved scattering-integral algorithm

Travel-time inversion for converted waves of ocean-bottom multi-component seismograms considering water-layer reverberations

Reflection-based traveltime and waveform inversion with second-order optimization

Data-domain wave equation reflection traveltime tomography

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

Memory‐efficient Frequency‐domain Gauss–Newton Method for Wave‐equation First‐arrival Traveltime Inversion

First-arrival waveform inversion method and its nmnerical test

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

Wave Equation Reflection Traveltime Inversion Using Gauss–Newton Optimization

Characteristic reflector-based wave equation reflection traveltime inversion

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

Automatic Traveltime Inversion Via Sparse Decomposition of Seismic Data

Computational aspects of finite-frequency traveltime inversion kernels

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

Inversion of Deep-Water Velocity Using the Munk Formula and the Seabed Reflection Traveltime: An Inversion Scheme That Takes the Complex Seabed Topography Into Account.