Markov Chain Monte Carlo Simulation to Estimate Material Properties of a Layered Half-space

Jin Ho Lee; Hieu Van Nguyen; Se Hyeok Lee

doi:10.7734/COSEIK.2023.36.3.203

All Issue

2023 Vol.36, Issue 3 Preview Page

Research Paper

Markov Chain Monte Carlo Simulation to Estimate Material Properties of a Layered Half-space 층상 반무한 지반의 물성치 추정을 위한 마르코프 연쇄 몬테카를로 모사 기법

30 June 2023. pp. 203-211

PDF XML

Abstract

A Markov chain Monte Carlo (MCMC) simulation is proposed for probabilistic full waveform inversion (FWI) in a layered half-space. Dynamic responses on the half-space surface are estimated using the thin-layer method when a harmonic vertical force is applied. Subsequently, a posterior probability distribution function and the corresponding objective function are formulated to minimize the difference between estimations and observed data as well as that of model parameters from prior information. Based on the gradient of the objective function, a proposal distribution and an acceptance probability for MCMC samples are proposed. The proposed MCMC simulation is applied to several layered half-space examples. It is demonstrated that the proposed MCMC simulation for probabilistic FWI can estimate probabilistic material properties such as the shear-wave velocities of a layered half-space.

Keywords

markov chain monte carlo simulation

full waveform inversion

inverse problem

layered half-space

thin-layer method

층상 반무한체에서의 확률론적 완전파형역산을 위한 Markov chain Monte Carlo (MCMC) 모사 기법을 정식화한다. Thin-layer method를 사용하여 조화 수직 하중이 작용하는 층상 반무한체의 지표면에서 추정된 동적 응답과 관측 데이터와의 차이 및 모델 변수의 사전 정보와의 차이를 최소화하도록 목적함수와 모델 변수의 사후 확률밀도함수를 정의한다. 목적함수의 기울기에 기반하여 MCMC 표본을 제안하기 위한 분포함수와 이를 수락 또는 거절할지 결정하는 수락함수를 결정한다. 기본 진동모드 뿐만이 아니라 고차 진동모드가 우세한 경우를 포함하여 다양한 층상 반무한체의 전단파 속도 추정에 제안된 MCMC 모사 기법을 적용하고 그 정확성을 검증한다. 제안된 확률론적 완전파형역산을 위한 MCMC 모사 기법은 층상 반무한체의 전단파 속도와 같은 재료 특성의 확률적 특성을 추정하는 데 적합함을 확인할 수 있다.

키워드

마르코프 연쇄 몬테카를로 모사

완전파형역산

역해석

층상 반무한체

thin-layer method

References

Astaneh, A.V., Guddati, M.N. (2016a) Improved Inversion Algorithms for Near-Surface Characterization, Geophys. J. Int., 206(2), pp.1410~1423. 10.1093/gji/ggw192

Astaneh, A.V., Guddati, M.N. (2016b) Efficient Computation of Dispersion Curves for Multilayered Waveguides and Half- spaces, Comput. Methods Appl. Mech. & Eng., 300, pp. 27~46. 10.1016/j.cma.2015.11.019

Eslaminia, M, Elmeliegy, A.M., Guddati, M.N. (2022) Full Waveform Inversion through Double-Sweeping Solver, J. Comput. Phys., 453, p.110914. 10.1016/j.jcp.2021.110914

Guddati, M.N., Druskin, V., Astaneh, A.V. (2016) Exponential Convergence through Linear Finite Element Discretization of Stratified Subdomains, J. Comput. Phys., 322, pp.492~447. 10.1016/j.jcp.2016.06.045

Kausel, E. (1981) An Explicit Solution for the Green's Functions for Dynamic Loads in Layered Media, MIT Research Report R81-13, Department of Civil Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA.

Kausel, E. (1996) Discussion: Dynamic Response of a Multi- Layered Poroelastic Medium by Rajapakse RKND, Senjuntichai T. Earthq. Eng. & Struct. Dyn., 25(10), pp.1165~1167. 10.1002/(SICI)1096-9845(199610)25:10<1165::AID-EQE600>3.0.CO;2-I

Lee, J.H., Lee, S.-H. (2021a) Full Waveform Inversion to Estimate the Material Properties of a Layered Half-Space, Soil Dyn. & Earthq. Eng., 151, p.106956. 10.1016/j.soildyn.2021.106956

Lee, J.H., Lee, S.-H. (2021b) Estimation of Shear-Wave Velocities of a Layered Half-space Using a Full Waveform Inversion with a Genetic Algorithm, Computational Methods in Applied Mechanics and Engineering, 34(4), pp.221~230. 10.7734/COSEIK.2021.34.4.221

Lee, J.H., Tassoulas, J.L. (2011) Consistent Transmitting Boundary with Continued-Fraction Absorbing Boundary Conditions for Analysis of Soil-Structure Interaction in a Layered Half- Space, Comput. Methods Appl. Mech. & Eng., 200, pp. 1509~1525. 10.1016/j.cma.2011.01.004

Lee, S.-H., Han, S.-W., Lee, J.H. (2022) Probabilistic Full- Waveform Inversion Considering the Discontinuities of the Material Properties of a Layered Half-Space, Int. J. Geomech., 22(11), p.04022207. 10.1061/(ASCE)GM.1943-5622.0002583

Martin, J., Wilcox, L.C., Burstedd, C., Ghattas, O. (2012) A Stochastic Newton MCMC Method for Large-scale Statistical Inverse Problems with Application to Seismic Inversion, SIAM J. Sci. Comput., 34(3), pp.A1460~A1487. 10.1137/110845598

Mashayekh, H., Kallivokas, L.F., Tassoulas, J.L. (2018) Parameter Estimation in Layered Media Using Dispersion-Constrained Inversion, J. Eng. Mech., 144(11), p.04018099. 10.1061/(ASCE)EM.1943-7889.0001527

Nguyen, V.H., Lee, S.-H., Lee, J.H. (2022) Application of the Unscented Kalman Filter to Estimate the Material Properties of a Layered Half-Space, Int. J. Appl. Mech., 15(2), p.2350005. 10.1142/S1758825123500059

Qi, Y., Minka, T.P. (2002) Hessian-based Markov Chain Monte- Carlo Algorithms, Proceedings of the First Cape Cod Workshop on Monte Carlo Methods, Massachusetts, USA: Cape Cod

Tarantola, A. (1987) Inverse Problem Theory: Method for Data Fitting and Model Parameter Estimation, Society for Industrial and Applied Mathematics, Philadelphia, USA.

Virieux, J., Operto, S. (2009) An Overview of Full-Waveform Inversion in Exploration Geophysics, Geophys., 74(6), pp. WCC127~WCC152. 10.1190/1.3238367

Zhao, Z., Sen, M.K. (2021) A Gradient-Based Markov Chain Monte Carlo Method for Full-Waveform Inversion and Uncertainty Analysis, Geophys., 86(1), pp.R15~R30. 10.1190/geo2019-0585.1

Information

Publisher :Computational Structural Engineering Institute of Korea
Publisher(Ko) :한국전산구조공학회
Journal Title :Journal of the Computational Structural Engineering Institute of Korea
Journal Title(Ko) :한국전산구조공학회 논문집
Volume : 36
No :3
Pages :203-211
Received Date : 2023-05-03
Revised Date : 2023-06-11
Accepted Date : 2023-06-12
DOI :https://doi.org/10.7734/COSEIK.2023.36.3.203

Journal of the Computational Structural Engineering Institute of Korea ISSN:1229-3059(Print) 2287-2302(Online) 한국전산구조공학회 논문집

All Issue