The performance of ray-based tomography inversion is affected by many factors,such as initial model error and low-velocity interlayer.The conventional tomography method based on first-arrival wave travel time,which constrains or smooths models,destroys the relative relationship between model parameters and rays and affects the inversion stability.By testing the performance of first-arrival wave travel time-based tomography inversion under different initial models,this study proposed a first-arrival wave travel time-based tomography inversion method with surface wave information as constraints.The process of this method is as follows:(1)Given that surface waves feature high energy and frequency dispersion in seismic data,the surface-wave frequency dispersion curves are obtained through the multi-channel analysis of surface waves;(2)Using the damped least squares method,the shallow-surface shear wave (S-wave) velocities are determined through inversion;(3)With the S-wave velocity structure as the constraint,the initial compressional wave (P-wave) model is established,and accordingly,the first-arrival wave travel time-based tomography inversion that considers regularization is achieved.This method improves the accuracy and stability of shallow structure inversion by fully utilizing the surface wave information in seismic data to counteract the inherent defects of tomography inversion.The effectiveness of this method has been verified using actual data.
Li H, Wang H Z, Zhang B. The study of regularization in tomography[J]. Geophysical Prospecting for Petroleum, 2015, 54(5):569-581.
doi: 10.3969/j.issn.1000-1441.2015.05.010
Zhang J H, Lyu N, Tian L Y, et al. An overview of the methods and techniques for seismic data noise attenuation[J]. Progress in Geophysics, 2005, 20(4):1083-1091.
Bi Y Y, Wang J J, Xu X H, et al. Ground roll attenuation based on the combination of discrete curvelet transform dictionary and two-dimensional local discrete cosine transform dictionary[J]. Geophysical Prospecting for Petroleum, 2017, 56(2):222-231.
doi: 10.3969/j.issn.1000-1441.2017.02.009
Li J W, Zang D G, Diao Y B, et al. Adaptive subtraction and Curvelet transforms combine to suppress surface waves[J]. Oil Geophysical Prospecting, 2020, 55(5):1005-1015.
Wu D S, Sun C Y, Lin M Y. High resolution dispersion imaging of surface waves based on multiple signal classification in frequency-velocity domain[J]. Geophysical Prospecting for Petroleum, 2017, 56(1):141-149.
doi: 10.3969/j.issn.1000-1441.2017.01.016
[10]
Stokoe K H, Nazarian S P. Effectiveness of ground improvement from spectral analysis of surface waves[C]// Helsinki:8th European Conference on Soil Mechanics and Foundation Engineering, 1983:91-94.
Liu J P, Hou W S, Xu S F. Adjacent-channel transient Rayleigh wave method and its application in compression strength test of water-tight wall[J]. Yangtze River, 2003, 34(2):34-36,56.
[12]
Xia J, Miller R D, Park C B. Estimation of near-surface shear-wave velocity by inversion of Rayleigh wave[J]. Geophysics, 1999, 64(3):1390-1395.
Wang Z N, Sun C Y, Wu D S. Near-surface 3D shear wave velocity modeling based on surface wave information[C]// Chinese Petroleum Society, 2019:1186-1189.
[14]
Li H X, et al. Multichannel analysis of surface waves based on short array stacked Correlation gather[J]. Soil Dynamics and Earthquake Engineering, 2021:146.
Deng X J, Feng S Y, Zuo Y, et al. Research on shallow structure using the surface wave and primary wave of shallow reflection seismicdata[J]. Journal of Geodesy and Geodynamics, 2019, 39(4):425-431.
Chen M, Wang Z H, Liu Z D, et al. Exploration of urban underground space resources:combined wave and first arrival tomography to detect near surface velocity structure in Jinan spring area[J]. Progress in Geophysics, 2022, 37(2):786-796.
Huang X G, Sun J G, Sun Z Q, et al. Calculation method for multiple travel based on the return function equation and the improved rapid propulsion method[J]. Oil Geophysical Prospecting, 2016, 51(6):1109-1118.
[18]
Li S W, Vladimirsky A, Fomel S. First-break traveltime tomography with the double-square-root eikonal equation[J]. Geophysics, 2013, 78(6):U89-U101.
doi: 10.1190/geo2013-0058.1
[19]
邓乐翔. 瑞雷波场正演模拟及频散曲线的提取[D]. 西安: 长安大学, 2010.
[19]
Deng L X. Rayleigh wave field simulation modeling and frequency dispersion curves of extraction.[D]. Xi’an: Chang’an University, 2010.
[20]
卢建旗. 多道面波分析方法及其应用研究[D]. 哈尔滨: 中国地震局工程力学研究所, 2013.
[20]
Lu J Q. Multi-channel analysis of surface wave and its utility[D]. Harbin: Institute of Engineering Mechanics,China Earthquake Administration, 2013.
[21]
伍敦仕. 基于面波信息的近地表参数反演方法研究[D]. 青岛: 中国石油大学(华东), 2018.
[21]
Wu D S. Study on inversion methods of near-surface parameters based on surface-wave information[D]. Qingdao: China University of Petroleum(East China), 2018.
[22]
Sun C Y, Wang Y Y, Wu D S, et al. Nonlinear Rayleigh wave inversion based on the shuffled frog-leaping algorithm[J]. Applied Geophysics, 2017, 14(4):551-558.
doi: 10.1007/s11770-017-0641-x
Zhao D, Wang G J, Wang X T, et al. The application of genetic algorithm to Rayleigh wave inversion[J]. Geophysical and Geochemical Exploration, 1995, 19(3):178-185.
Shen H Y, Li Q C, Yan Y Y, et al. Phase velocity analysis of multi-channel transient surface wave[J]. Geophysical Prospecting for Petroleum, 2016, 55(5):692-702.
doi: 10.3969/j.issn.1000-1441.2016.05.008
[27]
Fomel S. Shaping regularization in geophysical-estimation problems[J]. Geophysics, 2007, 72(2):R29-R36.
doi: 10.1190/1.2433716