Seismic ray tracing in models with complex surface topography using linear traveltime interpolation on a hybrid grid

jianguo sun,zhangqing sun,fuxing han

DOI: https://doi.org/10.1190/segam2013-0232.1

2013-01-01

SEG Technical Program Expanded Abstracts

Abstract:PreviousNext No AccessSEG Technical Program Expanded Abstracts 2013Seismic ray tracing in models with complex surface topography using linear traveltime interpolation on a hybrid gridAuthors: Jianguo SunZhangqing SunFuxing HanJianguo SunCollege of Geoexploration Science and Technology, Department of Geophysics, Jilin USearch for more papers by this author, Zhangqing SunCollege of Geoexploration Science and Technology, Department of Geophysics, Jilin USearch for more papers by this author, and Fuxing HanCollege of Geoexploration Science and Technology, Department of Geophysics, Jilin USearch for more papers by this authorhttps://doi.org/10.1190/segam2013-0232.1 SectionsSupplemental MaterialAboutPDF/ePub ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions ShareFacebookTwitterLinked InRedditEmail Abstract If the earth's surface has a complex geometry, the implementation of the grid method for solving the eikonal equation will become difficult. In fact, in this case we must first solve the following six problems, namely (1) how to describe the curved earth's surface, (2) how to divide the surface into discrete pieces, (3) how to describe the medium directly near the surface, (4) how to divide the medium directly near the surface into discrete elements, (5) how to compute the local traveltime, and (6) how to capture and to propagate the local wavefront near the curved earth's surface or interface. Previously, for solving the problems (1) to (5), we used a second order upwind finite difference scheme that uses nonuniform grid spacing in the regions near the earth's surface or the interface (Sun et al., 2011). For solving the problem (6), we modified the fast marching method (FMM) based on the narrow band technique by introducing new point types, namely the surface point, the point above the surface, the interface point, and the point under the interface (Sun et al., 2011). Here, for solving the problems listed above we present a method that uses a hybrid grid and that combines the linear interpolation method and the FMM. Specifically, the surface is approximated by broken lines (2D cases) and the medium directly near the surface is divided by an irregular grid that is composed of rectangular and triangular grid elements. For the local travetime computation and the local wavefront propagation, we use the linear interpolation method and the modified FMM, respectively. Furthermore, for constructing rays from the computed traveltimes, we search the point from which a ray comes from in reversed direction, i.e., from the bottom of the model to the top. Numerical results show that the method presented here not only can flexibly approximate the traveltimes, but also can give accurate ray trajectories in the models with strong heterogeneities. Permalink: https://doi.org/10.1190/segam2013-0232.1FiguresReferencesRelatedDetailsCited byA 3D reflection ray-tracing method based on linear traveltime perturbation interpolationTongyu Li, Jie Liu, and Jianzhong Zhang26 April 2019 | GEOPHYSICS, Vol. 84, No. 4Multiples ray tracing using multistage group marching upwind hybrid method in regional multiscale random media with topographyZhangqing Sun, Jianguo Sun, Fuxing Han, Xueqiu Wang, Yiteng Tian, and Peng Qi17 August 2017Diffraction Imaging and Modeling and Seismic Modeling Complete Session17 August 2017 SEG Technical Program Expanded Abstracts 2013ISSN (print):1052-3812 ISSN (online):1949-4645Copyright: 2013 Pages: 5258 Publisher:Society of Exploration Geophysicists HistoryPublished Online: 19 Aug 2013 CITATION INFORMATION Jianguo Sun, Zhangqing Sun, and Fuxing Han, (2013), "Seismic ray tracing in models with complex surface topography using linear traveltime interpolation on a hybrid grid," SEG Technical Program Expanded Abstracts : 3465-3469. https://doi.org/10.1190/segam2013-0232.1 Plain-Language Summary PDF DownloadLoading ...