三维复杂山地多级次群推进迎风混合法多波型走时计算

doi:10.6038/cjg20170521

摘要
图/表
参考文献
相关文章 (7)

全文: PDF (3403 KB) HTML (1 KB)
输出: BibTeX | EndNote (RIS)

摘要三维复杂山地条件下的各种地震波型的走时计算技术，可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时，提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题，利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题，利用群推进法模拟三维复杂地表条件下地震波前的扩展问题，利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明：新方法具有很好的计算精度与效率，且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	孙章庆
	孙建国
	王雪秋
	高正辉
	江兆南

关键词 ：三维复杂山地, 多级次算法, 群推进法, 迎风混合法, 多波型, 走时计算

Abstract：The techniques for computing traveltime of different types of seismic waves in mountainous areas with 3D complex topographic conditions can be directly used for analyzing the kinematical characteristics of seismic waves, designing the seismic data acquisition system, and developing some seismic data processing techniques directly for the 3D complex earth's surface. To compute traveltime of different types of seismic waves with high accuracy, great flexibility, and unconditional stability in 3D complex topographic conditions, we propose a multistage group marching upwind hybrid method. In this new method, we use the upwind finite-difference method with non-uniform grid spacing to compute the local seismic traveltime on or nearby the earth's surface with unconditional stability and concisely. The local seismic traveltime on the other uniform grids is computed by using the upwind bilinear interpolation method which has high accuracy. We also use the group marching method to simulate flexibly the extended processing of the seismic wave-front in 3D complex topographic conditions. Finally, the computational problem of multiple seismic traveltime is solved flexibly by the multistage algorithm. The contrast analysis to computational accuracy and the numerical example show that the new method presented in this paper has high accuracy and the flexible adaptability to solve the problem of multiple seismic traveltime computation in mountainous areas with complicated media and 3D complex surface topography.

Key words： 3D complex topography Multistage algorithm Group marching method Upwind hybrid method Multiple seismic wave Traveltime computation

收稿日期: 2016-04-28

PACS:

P631

基金资助:国家自然科学基金项目（41404085、41274120）、中国博士后科学基金第6批特别资助和第52批面上资助项目（2013T60329、2012M520680）、吉林大学优秀青年教师培养计划（419080500337）联合资助.

作者简介: 孙章庆,男,1982年生,博士,副教授,主要从事地震波传播理论与成像技术、直流电场正演模拟、计算地球物理等方面的研究工作.E-mail:sun_zhangq@jlu.edu.cn

链接本文:

http://www.geophy.cn/CN/10.6038/cjg20170521 或 http://www.geophy.cn/CN/Y2017/V60/I5/1861

Asakawa E, Kawanaka T. 1993. Seismic ray tracing using linear traveltime interpolation. Geophysical Prospecting, 41(1): 99-111.
Bai C Y, Greenhalgh S, Zhou B. 2007. 3D ray tracing using a modified shortest-path method. Geophysics, 72(4): T27-T36.
Bai C Y, Huang G J, Li Z S. 2011. Simultaneous inversion combining multiple-phase travel times within 3D complex layered media. Chinese J. Geophys. (in Chinese), 54(1): 182-192, doi: 10.3969/j.issn.0001-5733.2011.01.019.
Huang G J, Bai C Y. 2013. Simultaneous inversion of three model parameters with multiple classes of arrival times. Chinese J. Geophys. (in Chinese), 56(12): 4215-4225, doi: 10.6038/cjg20131224.
Huang G J, Bai C Y, Qian W. 2015. Simultaneous inversion of three model parameters with multiple phases of arrival times in spherical coordinates. Chinese J. Geophys. (in Chinese), 58(10): 3627-3638, doi: 10.6038/cjg20151016.
Jing X L, Yang C C, Wang S Q. 2007. A improved seismic reflection tomographic method. Chinese J. Geophys. (in Chinese), 50(6): 1831-1836.
Julian B R, Gubbins D. 1977. Three-dimensional seismic ray tracing. J. Geophys., 43: 95-113.
Kim S, Folie D. 2000. The group marching method: An O(N) level set eikonal solver.//70th Annual International Meeting, SEG, Expanded Abstracts, 2297-2300.
Lan H Q, Zhang Z J. 2013a. A high-order fast-sweeping scheme for calculating first-arrival travel times with an irregular surface. Bulletin of the Seismological Society of America, 103(3): 2070-2082.
Lan H Q, Zhang Z J. 2013b. Topography-dependent eikonal equation and its solver for calculating first-arrival traveltimes with an irregular surface. Geophys. J. Int., 193(2): 1010-1026.
Lan H Q, Chen J Y, Zhang Z J. 2014. A fast sweeping scheme for calculating P wave first-arrival travel times in transversely isotropic media with an irregular surface. Pure Appl. Geophys., 171(9): 2199-2208.
Li P M, Mei S Q, Ma Q B. 2013. An improved bilinear interpolation travel-time ray-tracing method. Oil Geophysical Prospecting (in Chinese), 48(4): 553-558.
Liu F, Zhang D, Yang Y, et al. 2012. A fast numerical solution to minimum equation in 3-D seismic LTI ray tracing. Journal of Wuhan University (Natural Science Edition) (in Chinese), 58(5): 395-400.
Mei S Q, Deng F, Zhong B S, et al. 2010. The 3D ray tracing method base on the improved bilinear traveltime interpolation. Computing Techniques for Geophysical and Geochemical Exploration (in Chinese), 32(2): 152-157.
Nakanishi I, Yamaguchi K. 1986. A numerical experiment on nonlinear image reconstruction from first-arrival times for two-dimensional island arc structure. Journal of Physics of the Earth, 34(2): 195-201.
Qu C, Zhou H L, Zhao D P. 2007. Deep structure beneath the west margin of Philippine Sea Plate and South China Sea from P and S wave travel time tomography. Chinese J. Geophys. (in Chinese), 50(6): 1757-1768.
Rawlinson N, Sambridge M. 2004. Multiple reflection and transmission phases in complex layered media using a multistage fast marching method. Geophysics, 69(5): 1338-1350.
Sethian J A, Popovici A M. 1999. 3-D traveltime computation using the fast marching method. Geophysics, 64(2): 516-523.
Sun J G. 2000. Limited-aperture migration. Geophysics, 65(2): 584-595.
Sun Z Q, Sun J G, Han F X. 2012. Traveltime computation using the upwind finite difference method with nonuniform grid spacing in a 3D undulating surface condition. Chinese J. Geophys. (in Chinese), 55(7): 2441-2449, doi: 10.6038/j.issn.0001-5733.2012.07.028.
Sun Z Q, Sun J G, Yue Y B, et al. 2015. 3D traveltime computation using fast marching upwind bilinear interpolation method. Chinese J. Geophys. (in Chinese), 58(6): 2011-2023, doi: 10.6038/cjg20150616.
Tang X P, Bai C Y. 2009. Multiple ray tracing within 3-D layered media with the shortest path method. Chinese J. Geophys. (in Chinese), 52(10): 2635-2643, doi: 10.3969/j.issn.0001-5733.2009.10.024.
Vidale J. 1988. Finite-difference calculation of travel times. Bulletin of the Seismological Society of America, 78(6): 2062-2076.
Vidale J E. 1990. Finite-difference calculation of traveltimes in three dimensions. Geophysics, 55(5): 521-526.
Vinje V, Iversen E, Gjøystdal H. 1993. Traveltime and amplitude estimation using wavefront construction. Geophysics, 58(8): 1157-1166.
Wang H, Chang X. 2000. 3-D ray tracing method based on graphic structure. Chinese J. Geophys. (in Chinese), 43(4): 534-541.
Wang H Z, Ma Z T. 1999. Traveltime calculation in 3D media with arbitrary velocity distribution.//69th Annual International Meeting, SEG, Expanded Abstracts, 1263-1266.
Zhang D, Fu X R, Yang Y, et al. 2009. 3-D seismic ray tracing algorithm based on LTI and partition of grid interface. Chinese J. Geophys. (in Chinese), 52(9): 2370-2376, doi: 10.3969/j.issn.0001-5733.2009.09.023.
Zhang D, Zhang T T, Qiao Y F, et al. 2013. A 3-D ray tracing method based on B-spline traveltime interpolation. Oil Geophysical Prospecting (in Chinese), 48(4): 559-566.
附中文参考文献
白超英, 黄国娇, 李忠生. 2011. 三维复杂层状介质中多震相走时联合反演成像. 地球物理学报, 54(1): 182-192, doi: 10.3969/j.issn.0001-5733.2011.01.019.
黄国娇, 白超英. 2013. 多震相走时联合三参数同时反演成像. 地球物理学报, 56(12): 4215-4225, doi: 10.6038/cjg20131224.
黄国娇, 白超英, 钱卫. 2015. 球坐标系下多震相走时三参数同时反演成像. 地球物理学报, 58(10): 3627-3638, doi: 10.6038/cjg20151016.
井西利, 杨长春, 王世清. 2007. 一种改进的地震反射层析成像方法. 地球物理学报, 50(6): 1831-1836.
李培明, 梅胜全, 马青坡. 2013. 一种改进的双线性插值射线追踪方法. 石油地球物理勘探, 48(4): 553-558.
刘锋, 张东, 杨艳等. 2012. 三维LTI射线追踪极小值方程的快速数值解法. 武汉大学学报(理学版), 58(5): 395-400.
梅胜全, 邓飞, 钟本善等. 2010. 基于改进的双线性旅行时插值的三维射线追踪. 物探化探计算技术, 32(2): 152-157.
瞿辰, 周蕙兰, 赵大鹏. 2007. 使用纵波和横波走时层析成像研究菲律宾海板块西边缘带和南海地区的深部结构. 地球物理学报, 50(6): 1757-1768.
孙章庆, 孙建国, 韩复兴. 2012. 三维起伏地表条件下地震波走时计算的不等距迎风差分法. 地球物理学报, 55(7): 2441-2449, doi: 10.6038/j.issn.0001-5733.2012.07.028.
孙章庆, 孙建国, 岳玉波等. 2015. 基于快速推进迎风双线性插值法的三维地震波走时计算. 地球物理学报, 58(6): 2011-2023, doi: 10.6038/cjg20150616.
唐小平, 白超英. 2009. 最短路径算法下三维层状介质中多次波追踪. 地球物理学报, 52(10): 2635-2643, doi: 10.3969/j.issn.0001-5733.2009.10.024.
王辉, 常旭. 2000. 基于图形结构的三维射线追踪方法. 地球物理学报, 43(4): 534-541.
张东, 傅相如, 杨艳等. 2009. 基于LTI和网格界面剖分的三维地震射线追踪算法. 地球物理学报, 52(9): 2370-2376, doi: 10.3969/j.issn.0001-5733.2009.09.023.
张东, 张婷婷, 乔友锋等. 2013. 三维旅行时场B样条插值射线追踪方法. 石油地球物理勘探, 48(4): 559-566.