Journal of Jilin University(Earth Science Edition) ›› 2019, Vol. 49 ›› Issue (4): 1160-1168.doi: 10.13278/j.cnki.jjuese.20180198

Previous Articles     Next Articles

Travel-Time Computation in VTI Media Based on Fast Marching Method

Xiao Han, Wang Deli   

  1. College of GeoExploration Sicence and Technology, Jilin University, Changchun 130026, China
  • Received:2018-07-23 Online:2019-07-26 Published:2019-07-26
  • Supported by:
    Supported by National Nature Science Foundation of China (2016ZX05026-002-003) and National Natural Science Foundation of China (41374108)

Abstract: Full-wave seismic inversion is a trend of geological structure interpretation. Seismic first-arrival travel time is an important parameter in processes of providing initial models like tomography. An efficient, accurate and stable method for travel time computation is very important for anisotropic modeling. Firstly, we use the perturbation theory and Taylor's theorem to expand the eikonal equation of the VTI media, by which we achieve the analytical solution of the homogeneous media (Alkhalifa2012); next, we introduce a fast marching method of isotropic media, and use the upwind difference to obtain the single point travel time. Combining with narrowband propulsion technique, we set a new travel-time computation method in VTI media based on the fast matching method. The calculated result is evaluated by comparing with the analytical values of the weak anisotropic homogenous model. The effectiveness and stability of the method is verified through calculating the travel time of the reflected waves in the layered medium model. The theoretical analysis and model trial show that this method is efficient and precise, and can be applied to anisotropic tomography and full wave field migration.

Key words: anisotropy, fast marching method, narrow band technique, travel time calculation

CLC Number: 

  • P631.4
[1] Um J, Thurber C. A Fast Algorithm for Two-Point Seismic Ray Tracing[J]. Bulletin of the Seismological Society of America, 1987, 77(3):972-986.
[2] Moser T J, Nolet G, Snieder A R. Ray Bending Revisited[J]. Bulletin of the Seismological Society of America, 1992, 82(1):259-288.
[3] Langan R T. Tracing of Rays Through Heterogeneous Media:An Accurate and Efficient Procedure[J]. Geophysics, 1985, 50(9):1456-1465.
[4] Dijkstra E W. A Note on Two Problems in Connection with Graphs[J]. Numerische Mathematics, 1959, 1(1):269-271.
[5] 孙章庆,孙建国,韩复兴.复杂地表条件下基于线性插值和窄带技术的地震波走时计算[J]. 地球物理学报, 2009, 52(11):2846-2853. Sun Zhangqing, Sun Jianguo, Han Fuxing. Traveltimes Computation Using Linear Interpolation and Narrow Band Technique Under Complex Topographical Conditions[J]. Chinese Journal of Geophysics, 2009, 52(11):2846-2853.
[6] 孙章庆, 孙建国, 岳玉波,等. 基于快速推进迎风双线性插值法的三维地震波走时计算[J]. 地球物理学报, 2015,58(6):2011-2023. Sun Zhangqing, Sun Jianguo, Yue Yubo, et al. 3D Traveltime Computation Using Fast Marching Upwind Bilinear Interpolation Method[J]. Chinese Journal of Geophysics, 2015, 58(6):2011-2023.
[7] 李培明, 梅胜全, 马青坡. 一种改进的双线性插值射线追踪方法[J]. 石油地球物理勘探, 2013, 48(4):553-558. Li Peiming, Mei Shengquan, Ma Qingpo. An Improved Bilinear Interpolation Traveltime Ray Tracing Method[J]. Oil Geophysical Prospecting, 2013, 48(4):553-558.
[8] 梅胜全,邓飞,钟本善,等.基于改进的双线性走时插值的三维射线追踪[J]. 物探化探计算技术,2010,32(2);152-157. Mei Shengquan, Deng Fei, Zhong Benshan, et al.The 3D Ray Tracing Method Base on the Improved Bilinear Traveltime Interpolation[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2010, 32(2):152-157.
[9] Alkhalifah T. Traveltime Computation with the Linearized Eikonal Equation for Anisotropic Media[J]. Geophysical Prospecting, 2000, 50(4):373-382.
[10] Vidale J. Finite-Difference Calculation of Travel Times[J]. Bulletin of the Seismological Society of America, 1988, 78(6):2062-2076.
[11] 孙章庆, 孙建国, 韩复兴. 三维起伏地表条件下地震波走时计算的不等距迎风差分法[J]. 地球物理学报, 2012, 55(7):2441-2449. Sun Zhangqing, Sun Jianguo, Han Fuxing. Traveltime Computation Using the Upwind Finite Difference Method with Nonuniform Grid Spacing in a 3D Undulating Surface Condition[J]. Chinese Journal of Geophysics, 2012, 55(7):2441-2449.
[12] 李兴旺, 白超英, 李晓玲. 四面体单元剖分下三维各向异性TI介质中多次波射线追踪[J]. 石油地球物理勘探, 2017, 52(1):48-55. Li Xingwang, Bai Chaoying, Li Xiaoling. Multiples Raytracing in 3D Anisotropic TI Media with a Tetrahedron Cell Model[J]. Oil Geophysical Prospecting, 2017, 52(1):48-55.
[13] 李晓玲, 白超英, 胡光义. 起伏层状TI介质中多次波射线追踪[J]. 石油地球物理勘探, 2013, 48(6):924-931. Li Xiaoling, Bai Chaoying, Hu Guangyi. Multiple Ray Tracing in an Undulated Layered TI Media[J]. Oil Geophysical Prospecting, 2013, 48(6):924-931.
[14] 高亮, 李幼铭, 陈旭荣,等. 地震射线辛几何算法初探[J]. 地球物理学报,2000, 43(3):402-410. Gao Liang, Li Youming, Chen Xurong,et al. An Attempt to Seismic Ray Tracing with Symplectic Algorithm[J]. Chinese Journal of Geophysics, 2000, 43(3):402-410.
[15] 秦孟兆, 陈景波. Maslov渐近理论与辛几何算法[J]. 地球物理学报, 2000, 43(4):522-533. Qin Mengzhao, Chen Jingbo. Maslov Asymptotic Theory and Symplectic Algorithm[J]. Chinese Journal of Geophysics, 2000, 43(4):522-533.
[16] Cameron M, Fomel S, Sethian J. Time-to-Depth Conversion and Seismic Velocity Estimation Using Time-Migration Velocity[J]. Geophysics, 2008, 73(5):205-210.
[17] Rawlinson N, Sambridge M. Multiple Reflection and Transmission Phases in Complex Layered Media Using a Multistage Fast Marching Method[J]. Geophysics, 2004, 69(5):1338-1350.
[18] Rawlinson N, Sambridge M, Saygin E. A Dynamic Objective Function Technique for Generating Multiple Solution Models in Seismic Tomography[J]. Geophysical Journal International, 2008, 174(1):295-308.
[19] Hao Q, Alkhalifah T. An Acoustic Eikonal Equation for Attenuating Transversely Isotropic Media with a Vertical Symmetry Axis[J]. Geophysics, 2017, 82(1):C9-C20.
[20] Sethian J A. AFast-Marching Level Set Method for Monotonically Advancing Fronts[J]. Proceedings of the National Academy of Sciences of the United States of America, 1996, 93(4):1591-1595.
[21] Sethian J A. Level Set Methods and Fast Marching Method[M]//Sethian J A. Level Set Methods and Fast Marching Methods. London:Cambridge University Press, 1999:400.
[22] Alkhalifah T, Fomel S. Implementing the Fast Marching Eikonal Solver:Spherical Versus Cartesian Coordinates[J]. Geophysical Prospecting, 2010, 49(2):165-178.
[23] Popovici A M, Sethian J A. 3-D Imaging Using Higher Order Fast Marching Traveltimes[J]. Geophysics, 2002, 67(10):604-609.
[24] Rawlinson N,Sambridge M. Wave Front Evolution in Strongly Heterogeneous Layered Media Using the Fast Marching Method[J]. Geophysical Journal of the Royal Astronomical Society, 2004, 156(3):631-647.
[25] Rawlinson N,Sambridge M, Rawlinson N, et al. The Fast Marching Method:An Effective Tool for Tomographic Imaging and Tracking Multiple Phases in Complex Layered Media[J]. Exploration Geophysics, 2005, 36(4):341-350.
[26] De Kool M, Rawlinson N,Sambridge M. A Practical Grid-Based Method for Tracking Multiple Refraction and Reflection Phases in Three-Dimensional Heterogeneous Media[J]. Geophysical Journal of the Royal Astronomical Society, 2006, 167(1):253-270.
[27] 张文生, 何樵登, 朱建伟,等. 横向各向同性介质中群速度的计算[J]. 物探化探计算技术, 1997,19(2):97-102. Zhang Wensheng, He Qiaodeng, Zhu Jianwei, et al. Calculation of Group Velocity in Transversely Isotropic Media[J]. Geophysical and Geochemical Exploration Technology, 1997, 19(2):97-102.
[28] 董良国, 魏建新. 横向各向同性介质中弹性波的物理模拟[J]. 石油物探, 1999, 38(1):76-84. Dong Liangguo, Wei Jianxin. Physical Simulation of Elastic Waves in Transversely Isotropic Media[J]. Petroleum Geophysical Exploration, 1999, 38(1):76-84.
[29] 孟庆生, 何樵登, 王德利. 均匀横向各向同性介质中P波及SV波的射线规律[J]. 吉林大学学报(地球科学版), 2002, 32(4):378-381. Meng Qingsheng, He Qiaodeng, Wang Deli. Ray Law of P Wave and SV Wave in Uniform Transversely Isotropic Medium[J]. Journal of Jilin University (Earth Science Edition), 2002, 32(4):378-381.
[30] 张建中. 三维TI介质中P波NMO速度及VSP走时联合反演[D]. 北京:中国地震局地质研究所, 2005. Zhang Jianzhong. P-Wave NMO Velocity and VSP Travel Time Joint Inversion in 3D TI Media[D]. Beijing:Institute of Geology, China Earthquake Administration, 2005.
[31] 郝奇. VTI介质速度和各向异性参数建模研究[D]. 长春:吉林大学, 2010. Hao Qi. Modeling of Velocity and Anisotropy Parameters of VTI Media[D]. Changchun:Jilin University, 2010.
[32] 白海军, 孙赞东, 王学军. 基于波前构建法的TTI介质射线追踪[J]. 石油地球物理勘探, 2011,46(a01):1-6. Bai Haijun, Sun Zandong, Wang Xuejun. TTI Media Ray Tracing Based on Wavefront Construction[J]. Oil Geophysical Prospecting, 2011, 46(a01):1-6.
[33] Alkhalifah T. Acoustic Approximations for Seismic Processing in Transversely Isotropic Media[J]. Geophysics, 1998, 63(2):623-631. doi:10.1190/1.1444361.
[34] Bender C M, Orszag S A. Advanced Mathematical Methods for Scientists and Engineers[M]. London:McGraw-Hill,1999.
[35] Alkhalifah T. Scanning Anisotropy Parameters in Complex Media[J]. Geophysics, 2011, 76(2):U13-U22. doi:10.1190/1.3553015.
[1] Yang Guoxin, Pang Yumao, , Ma Ruishan , Guo Xingwei, , Cao Hui , Cai Laixing. Influence of Porosity and Fluid Medium on Thermal Conductivity of Clastic Rocks [J]. Journal of Jilin University(Earth Science Edition), 2025, 55(5): 1728-1741.
[2] Long Gang, Shen Jinsong, Su Zhaoyang, Ran Shang. Time Domain Electromagnetic Response Simulation and Characteristic Analysis of Electroanisotropic Seafloor Massive Sulfide Deposits [J]. Journal of Jilin University(Earth Science Edition), 2025, 55(1): 274-288.
[3] Cao Changxin, Sun Hongyue. Optimization of Siphon Drainage Hole Spacing on Slope in Loess Region [J]. Journal of Jilin University(Earth Science Edition), 2021, 51(4): 1152-1159.
[4] Chen Xiangzhong, Wang Bin. AVOA Inversion Methodof Fractured Reservoir Based on Petrophysical Model [J]. Journal of Jilin University(Earth Science Edition), 2021, 51(1): 266-276.
[5] Feng Kai, Qin Ce. Anisotropy Forward Modeling of Magnetotelluric(MT) Adaptive Finite Element [J]. Journal of Jilin University(Earth Science Edition), 2020, 50(6): 1887-1896.
[6] Luo Teng, Feng Xuan, Guo Zhiqi, Liu Cai, Liu Xiwu. Seismic Inversion of Anisotropy Parameters of Fractured Reservoirs by Simulated Annealing and Particle Swarm Optimization [J]. Journal of Jilin University(Earth Science Edition), 2019, 49(5): 1466-1476.
[7] Hong Bo, Li Xi'an, Wang Li, Li Lincui. Permeability Anisotropy and Microstructure of Yan'an Q3 Loess [J]. Journal of Jilin University(Earth Science Edition), 2019, 49(5): 1389-1397.
[8] Zhang Wenqiang, Yin Changchun, Liu Yunhe, Zhang Bo, Ren Xiuyan. Field Continuation for MCSEM Forward Modeling and Identification of Anisotropic Media [J]. Journal of Jilin University(Earth Science Edition), 2019, 49(2): 578-590.
[9] Zhou Yang, Su Shengrui, Li Peng, Ma Hongsheng, Zhang Xiaodong. Microstructure and Mechanical Properties of Broken Phyllite [J]. Journal of Jilin University(Earth Science Edition), 2019, 49(2): 504-513.
[10] Yang Zhilong, Yin Changchun, Zhang Bo, Liu Yunhe, Ren Xiuyan, Hui Zhejian. 3-D Arbitrarily Anisotropic Modeling for Towed Marine DC Resistivity Method in Deep Ocean [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(6): 1845-1853.
[11] Deng Xinhui, Liu Cai, Guo Zhiqi, Liu Xiwu, Liu Yuwei. Full Wave Field Seismic Response Simulation and Analysis of Anisotropic Shale Reservoir in Luojia Area of Jiyang Depression [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(4): 1231-1243.
[12] Zhang Bing, Guo Zhiqi, Xu Cong, Liu Cai, Liu Xiwu, Liu Yuwei. Fracture Properties and Anisotropic Parameters Inversion of Shales Based on Rock Physics Model [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(4): 1244-1252.
[13] Hu Xinlei, Lü Yanfang, Sun Yonghe, Sun Tongwen. Comprehensive Quantitative Evaluation of Vertical Sealing Ability of Faults in Caprock: An Example of Ed2 Mudstone Caprock in Nanpu Sag [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(3): 705-718.
[14] Yin Changchun, Yang Zhilong, Liu Yunhe, Zhang Bo, Qi Yanfu, Cao Xiaoyue, Qiu Changkai, Cai Jing. Characteristics of 3D DC Resistivity Response for Arbitrary Anisotropic Models Using Circular Scanning Measurement [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(3): 872-880.
[15] Sun Jianguo, Li Yilong, Sun Zhangqing, Miao He. Computation of Seismic Traveltimes and Raypath Based on Model Parameterization [J]. Journal of Jilin University(Earth Science Edition), 2018, 48(2): 343-349.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] CHU Feng-you, SUN Guo-sheng,LI Xiao-min,MA Wei-lin, ZHAO Hong-qiao. The Growth Habit and Controlling Factors of the CobaltRich Crusts in Seamount of the Central Pacific[J]. J4, 2005, 35(03): 320 -0325 .
[2] ZHOU Li-ping, SHEN Xiang-dong, LI Xue-bin, BAI Zhong-qiang. Experiment Study of Mechanical Properties of Natual Pumice Powder Cement-Soil[J]. J4, 2009, 39(3): 492 -497 .
[3] HUANG Yu-long, WANG Pu-jun, SHAO Rui. Porosity and Permeability of Pyroclastic Rocks of the Yingcheng Formation in Songliao Basin[J]. J4, 2010, 40(2): 227 -236 .
[4] PAN Dian-qi,ZHANG Zu-pei,PAN Dian-cai,CHEN Yi-min,XU Rui. A Test Research on Longitudinal Wave Velocity of Artificial Frozen Clay Under Different Temperature and Moisture Conditions[J]. J4, 2006, 36(04): 588 -591 .
[5] FU Zhe,ZHOU Yun-xuan, LIU Dian-wei, LIU Wan-song. Design of the FeatureBased ObjectOriented Data Model and Implementation of Prototype System for a Virtual GIS[J]. J4, 2006, 36(04): 647 -652 .
[6] LU Jin-cai, LI Yu-hong, WEI Xian-yang, WEI Jian-she. Research on the Depositional Environment and Resources Potential of the Oil Shale in the Chang7 Member,Triassic Yanchang Formation in the Ordos Basin[J]. J4, 2006, 36(6): 928 -0932 .
[7] HUANG Ji-guo,HUANG Guo-xin,NIE Guang-zheng,WEI Hai-juan,LV Ai-min,WANG Xue-song. Experimental Study on the Treatment of Landfill Leachate in a TwoPhase(UBF-UASB) Anaerobic System at Mid-Temperature[J]. J4, 2007, 37(1): 144 -0147 .
[8] SUN Cai-zhi, LIU Yu-lan, YANG Jun. Research on the Ecological and Sustainable Groundwater Table Regulation in the Lower Liaohe River Plain[J]. J4, 2007, 37(2): 249 -254 .
[9] SHU Ping, QU Yan-ming, WANG Guo-jun, DING Ri-xin, AI Xing-bo,JI Xue-yan, TANG Hua-feng, BIAN Wei-hua, WANG Pu-jun. Geological features of the volcanic reservoirs of the Songliao Basin and their geophysical detection[J]. J4, 2007, 37(4): 726 -0733 .
[10] ZHAO Yu-yan, HAO Li-bo, ZHANG Zhi-li,LU Ji-long, SUN Guang-rui. Design and Realization of Metal Deposit Exploration Information System[J]. J4, 2008, 38(1): 161 -0166 .