一种变网格差分的快速行进法

doi:10.11720/wtyht.2019.1026

Abstract
Figure/Table
References
Related Citation (15)

Download: PDF (5782 KB) HTML (0 KB)
Export: BibTeX | EndNote (RIS)

Abstract

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.

Key words： grid-various double grid finite-difference fast marching method traveltime forward modeling

Received: 17 January 2018 Published: 20 February 2019

P631.4

	Service

	E-mail this article
	Add to my bookshelf
	Add to citation manager
	E-mail Alert
	RSS
	Articles by authors
	You-Juan HE
	Yu-Lei QIAO
	Li-Juan HOU
	Jun ZHU
	Gang GAO
	Peng WANG

Cite this article:

You-Juan HE,Yu-Lei QIAO,Li-Juan HOU, et al. A grid-variable finite-difference fast marching method[J]. Geophysical and Geochemical Exploration, 2019, 43(1): 199-208.

URL:

https://www.wutanyuhuatan.com/EN/10.11720/wtyht.2019.1026 OR https://www.wutanyuhuatan.com/EN/Y2019/V43/I1/199

Split mode of variable grid FMM

Schematic diagram of narrowband technology

Fig.4~Fig.7 the same
">

Relative error calculated by conventional first-order FMM traveltime
a—no angular point;b—add angular point;Fig.4~Fig.7 the same

Relative error calculated by conventional second-order FMM traveltime

Relative error calculated by grid-variable first-order FMM traveltime

Relative error calculated by grid-variable second-order FMM traveltime

Relative error calculated by double grid first-order FMM traveltime

Number of Grid points corresponding to different grid spacing

Number of grid points corresponding to different grid spacing

Time comparison about several calculations

Forward modeling with high-speed anomalous bodies

Three methods isochronal line for traveltime
a—the whole area;b—fractionated gain

Expanding wavefronts law

The abnormal body speed is 2300 m/s

The abnormal body speed is 2500 m/s

The abnormal body speed is 3000 m/s

Traveltime isochron of Marmousi model

[1]	王乾龙 . 基于快速推进法的三维随机介质地震波走时计算[D]. 长春:吉林大学, 2017.
[1]	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.
[5]	王飞, 曲昕馨, 刘四新 , 等. 一种新的基于多模板快速推进算法和最速下降法的射线追踪方法[J]. 石油地球物理勘探, 2014,49(6):1106-1114.
[5]	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.
[7]	张双杰, 朱培民, 赵仁基 , 等. 快速推进法计算精度分析及改进[J]. 工程地球物理学报, 2009,6(3):254-265.
[7]	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.
[9]	张风雪, 吴庆举, 李永华 , 等. FMM射线追踪方法在地震学正演和反演中的应用[J]. 地球物理学进展, 2010,25(4):1197-1205.
[9]	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.
[11]	李永博, 李庆春, 吴琼 . 快速行进法射线追踪提高旅行时计算精度和效率的改进措施[J]. 石油地球物理勘探, 2016,51(3):467-473.
[11]	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.
[12]	孙章庆, 孙建国, 韩复兴 . 针对复杂地形的三种地震波走时算法及对比[J]. 地球物理学报, 2012,55(2):560-568.
[12]	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.
[13]	朱生旺, 曲寿利, 魏修成 , 等. 变网格有限差分弹性波方程数值模拟方法[J]. 石油地球物理勘探, 2012,55(2):560-568.
[13]	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.
[14]	张志禹, 侯文婷, 苗永康 . 步长自适应的有限差分复杂地表波场数值模拟[J]. 地球物理学进展, 2017,32(3):1321-1330.
[14]	Zhang Z Y, Hou W T, Miao Y K . Seismic wave simulation method based on variable grid step[J]. Progress in Geophysics, 2017,32(3):1321-1330.
[15]	李振春, 张慧, 张华 . 一阶弹性波方程的变网格高阶有限差分数值模拟[J]. 石油地球物理勘探, 2008,43(6):711-716.
[15]	Li Z C, Zhang H, Zhang H . Variable-grid high-order finite-difference numeric simulation of first-order elastic wave equation[J]. Oil Geophysical Prospecting, 2008,43(6):711-716.
[16]	孙章庆 . 起伏地表条件下的地震波走时与射线路径计算[D]. 长春:吉林大学, 2011.
[16]	Sun Z Q . The seismic traveltimes and raypath computation under undulating earth's surface condition[D]. Changchun:Jilin University, 2011.