In this paper, we study an inverse problem consisting in the determination of the mechanical properties of a layered solid elastic medium in contact with a fluid medium by measuring the variation of the pressure in the fluid while propagating a seismic/acoustic wave. The estimation of mechanical parameters of the solid is obtained from the simulation of a seismic wave propagation model governed by a system of partial differential equations.
Two stochastic methods, Markov chain Monte Carlo with an accelerated version and simultaneous perturbation stochastic approximation, are implemented and compared with respect to cost and accuracy.
Bishop, T.N. , Bube, K.P., Cutler, R.T., Langan, R.T., Love, P.L., Resnick, J.R., Shuey, R.T., Spindler, D.A., and Wyld, H.W., Tomographic determination of velocity and depth in laterally varying media . Geophysics, 50, 903-923 (1985).
2.
Tarantola, A.Inverse Problem Theory and Methods for Model Parameter Estimation , SIAM, Philadelphia, PA, 2005
3.
Cerveny, V.Seismic Ray Theory, Cambridge University Press , Cambridge, 2001.
4.
Ma, C.-C. , Liu, S.-W., and Chang, C.-M.Inverse calculation of material parameters for a thin-layer system using transient elastic waves . Acoustical Society of America Journal, 112, 811 (2002 )
5.
Tsogka, C.Modélisation Mathématique et Numérique de la Propagation des Ondes Élastiques Tridimensionnelles dans les Milieux Fissurés. Ph.D thesis, Chapters 1 and 2, Université Paris Dauphine-Paris IX, 1999.
6.
Diaz, J. and Joly, P.Robust High Order Non Conforming Finite Element Formulation for Time Domain Fluid-Structure. Journal of Computational Acoustics , 13, 403-431 (2005).
7.
Plessix, R.-E.A review of the adjoint-state method for computing the gradient of a functional with geophysical applications. Geophysical Journal International, 167, 495-503 (2006)
8.
Niederreiter, H.Random Number Generation and Quasi-Monte Carlo Methods, CBMS-NSF Regional Conference Series in Applied Mathematics , Vol. 63, vi, p.241, Society of Industrial Applied Mathematics, Philadelphia, PA, 1992.
9.
Andres Christen, J. and Fox, C.MCMC Using Approximation. Journal of Computational and Graphical Statistics, 14, 795-810 (2005).
10.
Maryak, J.L. and Chin, D.C.Global Random Optimization by Simultaneous Perturbation Stochastic Approximation , in Proceedings of the American Control Conference , Arlington, VA, pp. 756-762, 2001.
11.
Spall, J.C.Implementation of the Simultaneous Perturbation Algorithm for Stochastic Optimization. IEEE Transactions of Aerospace and Electronic Systems, 34, 817-823 (1998).
12.
Spall, J.C.A stochastic approximation algorithm for large-dimensional systems in the Kiefer-Wolfowitz setting, in Proceedings of IEEE Conference on Decision and Control, pp. 1544-1548, 1988.
13.
Spall, J.C.Multivariate stochastic approximation using a simultaneous perturbation gradient approximation. IEEE Transactions on Automated Control, 37, 332-341 (1992).
14.
Spall, J.C.Introduction to Stochastic Search and Optimization: Estimation, Simulation, and Control, Chapter 7, Wiley, Hoboken, NJ, 2003.
15.
Bjorck, A.Numerical Methods for Least Squares Problems, SIAM, Philadelphia, PA, 1996.
16.
Robert, C.Méthodes de Monte Carlo par chaines de Markov, Economica , Paris, 1996.
17.
Meyer, S. , Christensen, N. and Nicholls, G.Fast Bayesian Reconstruction of Chaotic Dynamical Systems Via Extended Kalman Filtering . Physical Review, 3535-3542 (2001 ).
18.
Tan, M.S. , Fox, C. and Nicholls, G.Physics 707SC Course Materials, Chapters 7-9, Auckland University.
19.
Assous, F. and Collino, F.A numerical method for the exploration of sensitivity: the case of the identification of the 2D stratified elastic medium . Inverse Problems, 6, 487-513 (1990)
20.
Al Khoury, P., Chavant, G., Clément , F. and Hervé, P.Inversion of spectroscopic data, application on CO2 radiation of flame combustion. Inverse Problems in Science and Engineering, 13, 219-240 (2005)
21.
Gerecér, L. and Vágó, Z.The Mathematics of Noise-Free SPSA. Proceedings of the IEEE Conference on Decision and Control, pp. 4400-4405 (2001 ).
