• Previous Article
    3D Electrical Impedance Tomography reconstructions from simulated electrode data using direct inversion $ \mathbf{t}^{\rm{{\textbf{exp}}}} $ and Calderón methods
  • IPI Home
  • This Issue
  • Next Article
    Where Bayes tweaks Gauss: Conditionally Gaussian priors for stable multi-dipole estimation
October  2021, 15(5): 1121-1134. doi: 10.3934/ipi.2021031

Velocity modeling based on Rayleigh wave dispersion curve and sparse optimization inversion

1. 

Shandong University of Science and Technology, College of Earth Science and Engineering, Qingdao 266590, China

2. 

Key Laboratory of Petroleum Resources Research, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China

3. 

University of Chinese Academy of Sciences, Beijing 100049, China

4. 

Innovation Academy for Earth Science, Chinese Academy of Sciences, Beijing 100029, China

* Corresponding author: Y. F. Wang

Received  July 2020 Revised  February 2021 Published  October 2021 Early access  May 2021

Fund Project: This work is sponsored by the National Key R & D Program of The Ministry of Science and Technology of China (Nos. 2018YFC1504203, 2018YFC0603500) and the key programs IGGCAS-2019031 & SZJJ-201901

This paper studies the S wave velocity modeling based on the Rayleigh wave dispersion curve inversion. We first discuss the forward simulation, and present a fast root-finding method with cubic-order of convergence speed to obtain the Rayleigh wave dispersion curve. With the Rayleigh wave dispersion curve as the observation data, and considering the prior geological anomalies structural information, we establish a sparse constraint regularization model, and propose an iterative solution method to solve for the S wave velocity. Experimental tests are performed both on the theoretical models and on the field data. It indicates from the experimental results that our new inversion scheme possesses the characteristics of easy calculation, high computational efficiency and high precision for model characterization.

Citation: Yan Cui, Yanfei Wang. Velocity modeling based on Rayleigh wave dispersion curve and sparse optimization inversion. Inverse Problems and Imaging, 2021, 15 (5) : 1121-1134. doi: 10.3934/ipi.2021031
References:
[1]

X. F. Chen, A systematic and efficient method of computing normal modes for multilayered half-space, Geophysical Journal International, 115 (1993), 391-409.  doi: 10.1111/j.1365-246X.1993.tb01194.x.

[2]

X. F. Chen, Seismogram synthesis in multi-layered half-space Part Ⅰ, Theoretical formulations, Earthquake Research in China, 13 (1999), 149-174. 

[3]

H. A. DiaoP. J. Li and X. K. Yuan, Inverse elastic surface scattering with far-field data, Inverse Problems and Imaging, 13 (2019), 721-744.  doi: 10.3934/ipi.2019033.

[4]

M. V. de HoopG. Nakamura and J. Zhai, Unique recovery of piecewise analytic density and stiffness tensors from the elastic-wave Dirichlet-to-Neumann map, SIAM Journal on Applied Mathematics, 79 (2019), 2359-2384.  doi: 10.1137/18M1232802.

[5]

M. V. de Hoop, A. Iantchenko and R. D. van der Hilst and J. Zhai, Semiclassical inverse spectral problem for seismic surface waves in isotropic media Ⅰ: Love waves, Inverse Problems, 36 (2020), 075015, 27 pp. doi: 10.1088/1361-6420/ab8446.

[6]

M. V. de Hoop, A. Iantchenko and R. D. van der Hilst and J. Zhai, Semiclassical inverse spectral problem for seismic surface waves in isotropic media Ⅱ: Rayleigh waves, Inverse Problems, 36 (2020), 075016, 31 pp. doi: 10.1088/1361-6420/ab8447.

[7]

P. GabrielsR. Snieder and G. Nolet, In situ measurements of shear-wave velocity in sediments with higher-mode Rayleigh waves, Geophysical Prospecting, 35 (1987b), 187-196.  doi: 10.1111/j.1365-2478.1987.tb00812.x.

[8]

N. A. Haskell, The dispersion of surface waves on multilayered media, Bulletin of the Seismological Society of America, 43 (1953), 17-34. 

[9]

C. G. LaiG. J. RixS. Foti and V. Roma, Simultaneous measurement and inversion of surface wave dispersion and attenuation curves, Soil Dynamics and Earthquake Engineering, 22 (2002), 923-930.  doi: 10.1016/S0267-7261(02)00116-1.

[10]

G. Nakamura and G. Uhlmann, Inverse problems at the boundary for an elastic medium, SIAM Journal on Mathematical Analysis, 26 (1995), 263-279.  doi: 10.1137/S0036141093247494.

[11]

F. Schwab and L. Knopoff, Surface-wave dispersion computations, Bulletin of the Seismological Society of America, 60 (1970), 321-344. 

[12]

K. T. Tran and D. R. Hiltunen, One-dimensional inversion of full waveforms using a genetic algorithm, Journal of Environmental & Engineering Geophysics, 17 (2012), 197-213.  doi: 10.2113/JEEG17.4.197.

[13]

K. T. Tran and D. R. Hiltunen, One-dimensional inversion of full waveforms using a genetic algorithm, Journal of Geotechnical & Geoenvironmental Engineering, 138 (2012b), 1075-1090. 

[14]

Y. F. Wang and T.-Y. Xiao, Fast realization algorithms for determining regularization parameters in linear inverse problems, Inverse Problems, 17 (2001), 281-291.  doi: 10.1088/0266-5611/17/2/308.

[15] Y. F. WangI. E. StepnovaV. N. Titarenko and A. G. Yagola, Inverse Problems in Geophysics and Solution Methods, Higher Education Press, Beijing, 2011a. 
[16]

Y. F. Wang, A. G. Yagola and C. C. Yang (eds.), Optimization and Regularization for Computational Inverse Problems and Applications, Springer, Berlin, 2011b.

[17]

J. H. XiaR. D. Miller and C. B. Park, Estimation of near-surface shear-wave velocity by inversion of Rayleigh wave, Geophysics, 64 (1999), 691-700.  doi: 10.1190/1.1444578.

[18] J. H. Xia, High-Frequency Surface-wave Methods, China University of Geosciences Press, Wuhan, 2015. 
[19]

W. C. Yang, Theory and Methods of Geophysical Inversion, Geological Publishing House, Beijing, 1997.

show all references

References:
[1]

X. F. Chen, A systematic and efficient method of computing normal modes for multilayered half-space, Geophysical Journal International, 115 (1993), 391-409.  doi: 10.1111/j.1365-246X.1993.tb01194.x.

[2]

X. F. Chen, Seismogram synthesis in multi-layered half-space Part Ⅰ, Theoretical formulations, Earthquake Research in China, 13 (1999), 149-174. 

[3]

H. A. DiaoP. J. Li and X. K. Yuan, Inverse elastic surface scattering with far-field data, Inverse Problems and Imaging, 13 (2019), 721-744.  doi: 10.3934/ipi.2019033.

[4]

M. V. de HoopG. Nakamura and J. Zhai, Unique recovery of piecewise analytic density and stiffness tensors from the elastic-wave Dirichlet-to-Neumann map, SIAM Journal on Applied Mathematics, 79 (2019), 2359-2384.  doi: 10.1137/18M1232802.

[5]

M. V. de Hoop, A. Iantchenko and R. D. van der Hilst and J. Zhai, Semiclassical inverse spectral problem for seismic surface waves in isotropic media Ⅰ: Love waves, Inverse Problems, 36 (2020), 075015, 27 pp. doi: 10.1088/1361-6420/ab8446.

[6]

M. V. de Hoop, A. Iantchenko and R. D. van der Hilst and J. Zhai, Semiclassical inverse spectral problem for seismic surface waves in isotropic media Ⅱ: Rayleigh waves, Inverse Problems, 36 (2020), 075016, 31 pp. doi: 10.1088/1361-6420/ab8447.

[7]

P. GabrielsR. Snieder and G. Nolet, In situ measurements of shear-wave velocity in sediments with higher-mode Rayleigh waves, Geophysical Prospecting, 35 (1987b), 187-196.  doi: 10.1111/j.1365-2478.1987.tb00812.x.

[8]

N. A. Haskell, The dispersion of surface waves on multilayered media, Bulletin of the Seismological Society of America, 43 (1953), 17-34. 

[9]

C. G. LaiG. J. RixS. Foti and V. Roma, Simultaneous measurement and inversion of surface wave dispersion and attenuation curves, Soil Dynamics and Earthquake Engineering, 22 (2002), 923-930.  doi: 10.1016/S0267-7261(02)00116-1.

[10]

G. Nakamura and G. Uhlmann, Inverse problems at the boundary for an elastic medium, SIAM Journal on Mathematical Analysis, 26 (1995), 263-279.  doi: 10.1137/S0036141093247494.

[11]

F. Schwab and L. Knopoff, Surface-wave dispersion computations, Bulletin of the Seismological Society of America, 60 (1970), 321-344. 

[12]

K. T. Tran and D. R. Hiltunen, One-dimensional inversion of full waveforms using a genetic algorithm, Journal of Environmental & Engineering Geophysics, 17 (2012), 197-213.  doi: 10.2113/JEEG17.4.197.

[13]

K. T. Tran and D. R. Hiltunen, One-dimensional inversion of full waveforms using a genetic algorithm, Journal of Geotechnical & Geoenvironmental Engineering, 138 (2012b), 1075-1090. 

[14]

Y. F. Wang and T.-Y. Xiao, Fast realization algorithms for determining regularization parameters in linear inverse problems, Inverse Problems, 17 (2001), 281-291.  doi: 10.1088/0266-5611/17/2/308.

[15] Y. F. WangI. E. StepnovaV. N. Titarenko and A. G. Yagola, Inverse Problems in Geophysics and Solution Methods, Higher Education Press, Beijing, 2011a. 
[16]

Y. F. Wang, A. G. Yagola and C. C. Yang (eds.), Optimization and Regularization for Computational Inverse Problems and Applications, Springer, Berlin, 2011b.

[17]

J. H. XiaR. D. Miller and C. B. Park, Estimation of near-surface shear-wave velocity by inversion of Rayleigh wave, Geophysics, 64 (1999), 691-700.  doi: 10.1190/1.1444578.

[18] J. H. Xia, High-Frequency Surface-wave Methods, China University of Geosciences Press, Wuhan, 2015. 
[19]

W. C. Yang, Theory and Methods of Geophysical Inversion, Geological Publishing House, Beijing, 1997.

Figure 1.  Configuration and coordinate system of a multilayered half space
Figure 2.  Schematic diagram of the generalized reflection-transmission coefficient method. $ {{\rho }^{\left( j \right)}},{{\lambda }^{\left( j \right)}},{{\mu }^{\left( j \right)}} $ and $ {{h}^{\left( j \right)}} $ are the density, the Lamé constants and the layer thickness of the $ j $-th layer; $ {{\mathbf{C}}^{\left( j \right)}} = {{\left[ \mathbf{C}_{d}^{\left( j \right)},\mathbf{C}_{u}^{\left( j \right)} \right]}^{T}} $ are the corresponding coefficients of these downgoing and upgoing waves of the $ j $-th layer; $ {{\mathbf{R}}^{\left( j \right)}} = {{[\mathbf{R}_{ud}^{\left( j \right)},\mathbf{R}_{du}^{\left( j \right)}]}^{T}} $ and $ {{\mathbf{T}}^{\left( j \right)}} = {{\left[ \mathbf{T}_{u}^{\left( j \right)},\mathbf{T}_{d}^{\left( j \right)} \right]}^{T}} $ are the reflection coefficients and the transmission coefficients of the $ j $-th layer
Figure 3.  Inversion results using SIRW and LSM. (a) Data fitting for noiseless data: observed data (blue asterisks), dispersion curves corresponding to the initial model (brown dashed line), the models inverted from the noise-free data using SIRW (red solid line) and LSM (green solid line); (b) Inverted results with noiseless data: the true S wave velocity model (blue solid line), the initial S wave velocity model (blue dashed line), the S wave velocity models inverted by the SIRW (red solid line) and the LSM (green solid line) for the noise-free data; (c) Data fitting for noisy data: observed data (blue asterisks), dispersion curves corresponding to the initial model (brown dashed line), the models inverted from the noisy data using SIRW (red solid line) and LSM (green solid line); (d) Inverted results with noisy data: the true S wave velocity model (blue solid line), the initial S wave velocity model (blue dashed line), the S wave velocity models inverted by the SIRW (red solid line) and the LSM (green solid line) for the noisy data
Figure 4.  Inversion results for the field data (a) Seismic record; (b) Normalized power spectrum and the extracted fundamental dispersion curve (red solid line). (c) Initial S wave velocity model (green dashed line), S wave velocity structure inverted by the SIRW (red solid line) and S wave logging results (blue dots). (d) A geologic profile ([13])
Table 1.  Parameters of the five-layer synthetic model
Layer S wave velocity P wave velocity Density Thickness
numbers (m/s) (m/s) (g/cm3) (m)
1 300 560 1.7 2
2 340 640 1.8 5
3 400 750 1.8 3
4 360 680 1.9 4
5 500 900 2.0 $ \infty $
Layer S wave velocity P wave velocity Density Thickness
numbers (m/s) (m/s) (g/cm3) (m)
1 300 560 1.7 2
2 340 640 1.8 5
3 400 750 1.8 3
4 360 680 1.9 4
5 500 900 2.0 $ \infty $
[1]

Darryl D. Holm, Ruiao Hu. Nonlinear dispersion in wave-current interactions. Journal of Geometric Mechanics, 2022  doi: 10.3934/jgm.2022004

[2]

Antti Lipponen, Aku Seppänen, Jari Hämäläinen, Jari P. Kaipio. Nonstationary inversion of convection-diffusion problems - recovery from unknown nonstationary velocity fields. Inverse Problems and Imaging, 2010, 4 (3) : 463-483. doi: 10.3934/ipi.2010.4.463

[3]

Christina Surulescu, Nicolae Surulescu. Modeling and simulation of some cell dispersion problems by a nonparametric method. Mathematical Biosciences & Engineering, 2011, 8 (2) : 263-277. doi: 10.3934/mbe.2011.8.263

[4]

Atul Kumar, R. R. Yadav. Analytical approach of one-dimensional solute transport through inhomogeneous semi-infinite porous domain for unsteady flow: Dispersion being proportional to square of velocity. Conference Publications, 2013, 2013 (special) : 457-466. doi: 10.3934/proc.2013.2013.457

[5]

Jihong Zhao, Ting Zhang, Qiao Liu. Global well-posedness for the dissipative system modeling electro-hydrodynamics with large vertical velocity component in critical Besov space. Discrete and Continuous Dynamical Systems, 2015, 35 (1) : 555-582. doi: 10.3934/dcds.2015.35.555

[6]

Jibin Li, Yan Zhou. Bifurcations and exact traveling wave solutions for the nonlinear Schrödinger equation with fourth-order dispersion and dual power law nonlinearity. Discrete and Continuous Dynamical Systems - S, 2020, 13 (11) : 3083-3097. doi: 10.3934/dcdss.2020113

[7]

P. M. Jordan, Louis Fishman. Phase space and path integral approach to wave propagation modeling. Conference Publications, 2001, 2001 (Special) : 199-210. doi: 10.3934/proc.2001.2001.199

[8]

Stéphane Junca, Bruno Lombard. Stability of neutral delay differential equations modeling wave propagation in cracked media. Conference Publications, 2015, 2015 (special) : 678-685. doi: 10.3934/proc.2015.0678

[9]

Evelyn Herberg, Michael Hinze, Henrik Schumacher. Maximal discrete sparsity in parabolic optimal control with measures. Mathematical Control and Related Fields, 2020, 10 (4) : 735-759. doi: 10.3934/mcrf.2020018

[10]

Bruno Sixou, Valentina Davidoiu, Max Langer, Francoise Peyrin. Absorption and phase retrieval with Tikhonov and joint sparsity regularizations. Inverse Problems and Imaging, 2013, 7 (1) : 267-282. doi: 10.3934/ipi.2013.7.267

[11]

Hung-Wen Kuo. The initial layer for Rayleigh problem. Discrete and Continuous Dynamical Systems - B, 2011, 15 (1) : 137-170. doi: 10.3934/dcdsb.2011.15.137

[12]

Jinqiao Duan, Vincent J. Ervin, Daniel Schertzer. Dispersion in flows with obstacles and uncertainty. Conference Publications, 2001, 2001 (Special) : 131-136. doi: 10.3934/proc.2001.2001.131

[13]

Chien Hsun Tseng. Analytical modeling of laminated composite plates using Kirchhoff circuit and wave digital filters. Journal of Industrial and Management Optimization, 2020, 16 (5) : 2213-2252. doi: 10.3934/jimo.2019051

[14]

Brenton LeMesurier. Modeling thermal effects on nonlinear wave motion in biopolymers by a stochastic discrete nonlinear Schrödinger equation with phase damping. Discrete and Continuous Dynamical Systems - S, 2008, 1 (2) : 317-327. doi: 10.3934/dcdss.2008.1.317

[15]

Le Li, Lihong Huang, Jianhong Wu. Flocking and invariance of velocity angles. Mathematical Biosciences & Engineering, 2016, 13 (2) : 369-380. doi: 10.3934/mbe.2015007

[16]

Deren Han, Zehui Jia, Yongzhong Song, David Z. W. Wang. An efficient projection method for nonlinear inverse problems with sparsity constraints. Inverse Problems and Imaging, 2016, 10 (3) : 689-709. doi: 10.3934/ipi.2016017

[17]

Ville Kolehmainen, Matthias J. Ehrhardt, Simon R. Arridge. Incorporating structural prior information and sparsity into EIT using parallel level sets. Inverse Problems and Imaging, 2019, 13 (2) : 285-307. doi: 10.3934/ipi.2019015

[18]

Jian Lu, Lixin Shen, Chen Xu, Yuesheng Xu. Multiplicative noise removal with a sparsity-aware optimization model. Inverse Problems and Imaging, 2017, 11 (6) : 949-974. doi: 10.3934/ipi.2017044

[19]

Linh V. Nguyen. A family of inversion formulas in thermoacoustic tomography. Inverse Problems and Imaging, 2009, 3 (4) : 649-675. doi: 10.3934/ipi.2009.3.649

[20]

Birol Yüceoǧlu, ş. ilker Birbil, özgür Gürbüz. Dispersion with connectivity in wireless mesh networks. Journal of Industrial and Management Optimization, 2018, 14 (2) : 759-784. doi: 10.3934/jimo.2017074

2021 Impact Factor: 1.483

Metrics

  • PDF downloads (371)
  • HTML views (245)
  • Cited by (0)

Other articles
by authors

[Back to Top]