A grid-variable finite-difference fast marching method
You-Juan HE1, Yu-Lei QIAO2, Li-Juan HOU3, Jun ZHU4, Gang GAO1, Peng WANG1
1. Key Laboratory of Exploration Technologies for Oil and Gas Resources(Yangtze University),Ministry of Education,Wuhan 430100,China 2. Shengli Oil Field Exploration and Development Research Institute,SINOPEC,Dongying 257000,China 3. No.1 Oil Production Plant,Qinghai Oil Field,PetroChina,Haixi 816400,China 4. No.1 Oil Production Plant,Northwest Oilfield Company,SINOPEC,Luntai 841600,China
The accuracy of seismic wave traveltimes directly affects the reliability of research results in such fields as seismic inversion,pre-stack migration imaging and tomography.Therefore,it is of great significance to study the improvement of the accuracy of seismic wave traveltimes.Based on the double grid technology,this paper comes up with a fast marching method (FMM) based on the grid-various finite-difference scheme to calculate the traveltimes of seismic wave.It analyzes the advantages and applicability of the grid-various finite-difference FMM by the forward simulation of uniform model as well as the existence of high-speed anomalous body model,Marmousi model.The results show that the corner points need to be included in the calculation when the traveltimes are calculated by using Eikonal equation so as to reduce the error.Under the background of uniform model,the grid-variable finite-difference FMM has the same advantages as the double grid FMM.Nevertheless,under the background of the existence of high-speed anomalous body model,the double grid FMM may violate the law of wavefront expansion to cause a greater error.The grid-variable finite-difference FMM does not have such a problem,and its advantage is remarkable.Therefore,this method is an effective way to improve the accuracy and efficiency of traveltime calculation,which not only enhances the applicability of the FMM but also expands the application range of grid-various technology.
Wang Q L . The calculation of seismic traveltime based on FMM in 3-D random medium[D]. Changchun:Jilin University, 2017.
[2]
Cervený V. Seismic ray theory[M]. Cambridge: Cambridge University Press, 2001.
[3]
Zhao A H, Zhang Z J, Teng J W .Minimum traveltime tree algorithm for seismic ray tracing:improvement in efficieney[J].J. Geophys. Eng., 2004(1):245-251.
[4]
Vidale J E . Finite-difference calculations of traveltimes[J]. Bull Seism Soc Am, 1988,78:2062-2076.
Wang F, Qu X X, Liu S X , et al. A new ray tracing approach based on both multistencils fast marching and the steepest descent[J]. Oil Geophysical Prospecting, 2014,49(6):1106-1114.
[6]
Sethian J A, Popovici A M . 3-D traveltime computation using the fast marching method[J]. Geophysics, 1999,64(2):516-523.
doi: 10.1190/1.1444558
Zhang S J, Zhu P M, Zhao R J , et al. Analysis on calculation accuracy of fast marching method and its improvement[J]. Chinese Journal of Engineering Geophysics, 2009,6(3):254-265.
[8]
Rawlinson N, Sambridge M . Multiple reflection and transmission phases in complex layered media using a multistage fast marching method[J]. Geophysics, 2004,69:1338-1350.
doi: 10.1190/1.1801950
Zhang F X, Wu Q J, Li Y H , et al. Application of FMM ray tracing to forward and inverse problems of seismology[J]. Progress in Geophysics, 2010,25(4):1197-1205.
[10]
孙章庆 . 复杂地表条件下地震波走时计算方法研究[D]. 长春:吉林大学, 2007.
[10]
Sun Z Q . Study on the computation method of seismic traveltimes including surface topography[D]. Changchun:Jilin University, 2007.
Li Y B, Li Q C, Wu Q . Improved fast marching method for higher calculation accuracy and efficiency of traveltime[J]. Oil Geophysical Prospecting, 2016,51(3):467-473.
Sun Z Q, Sun J G, Han F X . The comparison of three schemes for computing seismic wave traveltimes in complex topographical conditions[J]. Chinese Journal of Geophysics, 2012,55(2):560-568.
Zhu S W, Qu S L, Wei X C , et al. Numeric simulation by grid various finite difference elastic wave equation[J]. Oil Geophysical Prospecting, 2007,42(6):634-639.