High-precision reflection traveltime tomography for velocity structure with an irregular surface
Zhang Xin-Yan,Xu Tao,Bai Zhi-Ming,Gao Rui,Li Qiu-Sheng,Liu You-Shan,Zhang Zhi
DOI: https://doi.org/10.6038/cjg20170209
2017-01-01
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.