Optimal First-Arrival Travel-Time Tomography with Multi-Information Constraints for Complex Near-Surface Scheme
PAN Yiming,SANG Yunyun,YAO Xuefeng,ZHANG Kai,LI Zhenchun,XU Xin
DOI: https://doi.org/10.3969/j.issn.1000-1441.2023.02.008
2023-01-01
Abstract:First-arrival travel-time tomography is an effective method for accurately calculating complex near-surface velocity structures.However, it is difficult to apply the traditional first-arrival travel-time tomography algorithm to seismic data prospecting with large data volumes, owing to the need for excessive memory and calculation time.The accuracy of near-surface velocity modeling of super-large offset data and complex piedmont data in western China is insufficient, which seriously affects the quality of subsequent mid-deep modeling and migration imaging.Therefore, it is very important to study first-arrival travel-time tomography for complex near-surface schemes with high computational efficiency and accuracy.In this study, we present an optimal first-arrival travel-time tomography method with multi-information constraints.The underlying idea is to introduce prior information into the objective function through the regularization method and improve the calculation efficiency through the optimization of forward and inversion algorithms.We use the multi-stencil fast marching method（MSFM） to solve the Eikonal equation to calculate the travel-time of grid elements and use an iterative inversion method to linearize the nonlinear inverse problem.The optimal algorithm can improve the accuracy of the first-arrival travel-time calculation of complex near-surface structures and avoid large-scale Frechet matrix calculation and storage, which greatly improves the computational efficiency.In terms of objective function constraints, we adopt Tikhonov regularization to solve the inverse problem of minimizing data difference and model smoothness and increase prior information such as apparent slowness to constrain the solution space of the model.The traditional objective function only fits the travel time of the first-arrival wave; therefore, it is difficult to distinguish the receiving points with the same travel-time and different positions.When dealing with complex near-surface problems, it is easy to observe the phenomenon of "bull’s eye（high-speed abnormal body）".The objective function combined with apparent slowness information can be used to fit the first-arrival travel-time and travel-time offset curves, which can effectively distinguish the receiving points with the same travel-time and different positions and achieve better results for near-surface problems.Adding slowness to the regularization operator to normalize the model can also avoid the problem of decreasing the resolution of inversion in the high-speed region.These methods reduce the ill-conditioning of the target equation and improve the inversion accuracy.In addition, the initial velocity model based on micrologging information can make the inversion result closer to the real velocity structure.The undulating surface model and actual data processing show that this method can effectively solve the problem of "high-speed abnormal bodies in the first-arrival travel-time tomography of 3D seismic data.Compared with traditional first-arrival travel-time tomography, our surface travel-time tomography improves inversion stability and calculation accuracy.