2D inversion of seismic first-arrival traveltime based on FCM clustering model constraint
LIU Jiacheng1, ZHANG Zhiyong1, ZHOU Qinyuan2, LI Man1, LI Hongli3
1. School of Geophysics and Measurement-control Technology, East China University of Technology, Nanchang, Jiangxi 330013, China; 2. Hunan Transportation Planning, Survey and Design Institute Co., LTD., Changsha, Hunan 410008, China; 3. Guangxi University of Science and Technology, Liuzhou, Guangxi 545006, China
Abstract:The 2D inversion of seismic first-arrival traveltime with minimum structure model constraint regularization suffers from the issue of insufficient delineation in model boundaries, particularly when dealing with sparse ray distribution within geological structures, leading to unsatisfactory inversion results. To address this challenge, this paper introduces a model constraint function based on fuzzy C-means (FCM) clustering, so as to improve the accuracy of the inversion results in delineating the model boundaries. This constraint incorporates prior information as reference cluster centers and employs an iterative process to repeatedly modify the cluster centers and the membership degrees of each grid cell to the cluster centers, enabling automatic classification of velocities. On this basis, a multi-grid inversion strategy guided by model sensitivity information is adopted to enhance the stability and effectiveness of the inversion. Parameter selection schemes for FCM clustering model constraint weights and weights of the prior information guidance term are discussed using simplified models. A comparison is made between the inversion results of unsupervised learning and prior information-supervised learning approaches. The latter approach successfully addresses the issue of blurring in depicting the boundaries of the velocity model during inversion, effectively improving the resolution of the inversion results. Finally, the inversion of real measured data verifies the applicability and effectiveness of the approach in practical applications.
JULIAN B R, GUBBINS D. Three-dimensional seismic ray tracing[J]. Journal of Geophysics, 1977, 43(1):95-113.
[5]
ČERVENY V, PŠENČÍK I. Gaussian beams and paraxial ray approximation in three-dimensional elastic inhomogeneous media[J]. Journal of Geophysics, 1983, 53(1):1-15.
[6]
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.
[7]
VIDALE J.Finite-difference calculation of travel times[J]. Bulletin of the Seismological Society of America, 1988, 78(6):2062-2076.
[8]
SETHIAN J A, POPOVICI A M.3-D traveltime computation using the fast marching method[J].Geophysics, 1999, 64(2):516-523.
[9]
王迪, 吴翔, 白超英, 等. 倾斜正交各向异性介质中多次反射qP波射线追踪[J]. 石油地球物理勘探, 2023, 58(3):632-640.WANG Di, WU Xiang, BAI Chaoying, et al.Seismic ray tracing for multiple reflected qP wave in tilted orthorhombic anisotropic media[J]. Oil Geophysical Prospecting, 2023, 58 (3):632-640.
[10]
LELIÈVRE P G, FARQUHARSON C G, HURICH C A. Computing first-arrival seismic traveltimes on unstructured 3-D tetrahedral grids using the fast marching method[J]. Geophysical Journal International, 2011, 184(2):885-896.
[11]
WHITE D J. Two-dimensional seismic refraction tomography[J]. Geophysical Journal International, 1989, 97(2):223-245.
[12]
PAIGE C C, SAUNDERS M A.LSQR:An algorithm for sparse linear equations and sparse least squares[J]. ACM Transactions on Mathematical Software, 1982, 8(1):43-71.
[13]
KENNETT B L N, SAMBRIDGE M S, WILLIAMSON P R. Subspace methods for large inverse problems with multiple parameter classes[J]. Geophysical Journal International, 1988, 94(2):237-247.
[14]
HIRAHARA K.Detection of three-dimensional velocity anisotropy[J].Physics of the Earth and Planetary Interiors, 1988, 51(1-3):71-85.
[15]
刘玉柱, 董良国, 夏建军.初至波走时层析成像中的正则化方法[J].石油地球物理勘探, 2007, 42(6):682-685.LIU Yuzhu, DONG Liangguo, XIA Jianjun.Normalized approach in tomographic imaging of first breaks travel time[J].Oil Geophysical Prospecting, 2007, 42(6):682-685.
[16]
LELIEVRE P G, FARQUHARSON C G, HURICH C A.Inversion of first-arrival seismic traveltimes without rays, implemented on unstructured grids[J].Geophysical Journal International, 2011, 185(2):749-763.
[17]
MAAG E, LI Y. Discrete-valued gravity inversion using the guided fuzzy C-means clustering technique[J].Geophysics, 2018, 83(4):G59-G77.
[18]
谭家炜, 李静, 李飞达, 等.基于超虚干涉法约束的模糊C均值聚类地震初至自动提取[J].石油地球物理勘探, 2020, 55(5):979-990.TAN Jiawei, LI Jing, LI Feida, et al.Automatic pick-up of seismic P-wave first arrivals via fuzzy C-means method constrained by super-virtual interferometry[J]. Oil Geophysical Prospecting, 2019, 55(5):979-990.
[19]
易柯, 张志勇, 李曼, 等. 基于FCM聚类算法的二维RMT电阻率与介电常数联合反演[J]. 地球物理学报, 2022, 65(6):2340-2350.YI Ke, ZHANG Zhiyong, LI Man, et al.Joint inversion of resistivity and permittivity for two dimensional RMT data based on FCM clustering[J]. Chinese Journal of Geophysics, 2022, 65(6):2340-2350.
[20]
SUN J, LI Y.Joint inversion of multiple geophysical data using guided fuzzy C-means clustering[J].Geophysics, 2016, 81(3):ID37-ID57.
[21]
李曼, 于鹏, 张志勇, 等.二维直流电阻率与音频大地电磁自适应渐进联合反演[J]. 同济大学学报(自然科学版), 2020, 48(1):114-122.LI Man, YU Peng, ZHANG Zhiyong, et al. Study of 2-D joint inversion of direct current resistivity and audio magnetotelluric data using adaptive progressive mesh refinement strategy[J]. Journal of Tongji University (Natural Science), 2020, 48(1):114-122.
[22]
程三, 张志勇, 周峰, 等. 非结构化网格渐进自适应正则化反演算法[J]. 石油地球物理勘探, 2022, 57(2):467-477.CHENG San, ZHANG Zhiyong, ZHOU Feng, et al. Study on step-by-step regularization inversion based on adaptive unstructured mesh[J]. Oil Geophysical Prospecting, 2022, 57(2):467-477.
[23]
ZHDANOV M S. Geophysical Inverse Theory and Regularization Problems[M]. Elsevier, Amsterdam, 2002.
[24]
LELIEVRE P G, FARQUHARSON C G. Gradient and smoothness regularization operators for geophysical inversion on unstructured meshes[J]. Geophysical Journal International, 2013, 195(1):330-341.
[25]
李曼, 林文东. 基于最小支持的2.5维直流电阻率正则化反演研究[J]. 东华理工大学学报(自然科学版), 2014, 37(3):292-298.LI Man, LIN Wendong. The study of regularization inversion for 2.5 dimension DC resistivity based on minimum support stabilizing factor[J]. Journal of East China University of Technology (Natural Science Edition), 2014, 37(3):292-298.