Rayleigh wave dispersion curve inversion via genetic algorithms and Marginal Posterior Probability Density estimation |
| |
| Institution: | 1. Institute of Rock Structure and Mechanics, Academy of Sciences of the Czech Republic, Prague, Czech Republic;2. Studio Tecnogeo2000, Udine, Italy;3. Datamind — Scientific and Technological Park, Udine, Italy;1. Departamento de Física, Ingeniería de Sistemas y Teoría de la Señal, Universidad de Alicante, Ap. Correos 99, 03080, Alicante, Spain;2. Instituto Universitario de Física Aplicada a las Ciencias y las Tecnologías, Universidad de Alicante, Ap. Correos 99, 03080, Alicante, Spain;3. Departamento de Ciencias de la Tierra y del Medio Ambiente, Universidad de Alicante, Ap. Correos 99, 03080 Alicante, Spain;4. Departamento de Ingeniería Civil, Universidad de Granada, Campus Fuentenuevax, Ap. Correos/n, 18071 Granada, Spain;5. Departamento de Física, Campus Las Lagunillas, Universidad de Jaén, 23071, Jaén, Spain;6. Dpt. Scienze della Terra. “Sapienza”, Università di Roma, P.le Aldo Moro 5, 00185 Roma, Italy;1. Institute of Rock Structure and Mechanics, Academy of Sciences of the Czech Republic, Prague, Czech Republic;2. Faculdade de Ciências da Universidade do Porto (DGAOT), Porto, Portugal;3. Geology and Geophysics Department, Faculty of Sciences, King Saud University, Riyadh, Saudi Arabia;4. Seismology Dept., National Research Institute of Astronomy and Geophysics, Cairo, Egypt;1. Hubei Subsurface Multi-scale Imaging Lab (SMIL), Institute of Geophysics and Geomatics, China University of Geosciences, 388 Lumo Road, Wuhan, Hubei 430074, China;2. Subsurface Imaging and Sensing Laboratory, Institute of Geophysics and Geomatics, China University of Geosciences, 388 Lumo Road, Wuhan, Hubei 430074, China;1. State Key Laboratory of Coal Resources and Safe Mining, China University of Mining & Technology (Beijing), Beijing 100083, China;2. College of Geoscience and Surveying Engineering, China University of Mining & Technology (Beijing), Beijing 100083, China;1. Hubei Subsurface Multi-scale Imaging Lab (SMIL), Institute of Geophysics and Geomatics, China University of Geosciences, 388 Lumo Road, Wuhan, Hubei 430074, China;2. Subsurface Imaging and Sensing Laboratory, Institute of Geophysics and Geomatics, China University of Geosciences, 388 Lumo Road, Wuhan, Hubei 430074, China |
| |
| Abstract: | Surface wave dispersion curve inversion is a challenging problem for linear inversion procedures due to its highly non-linear nature and to the large numbers of local minima and maxima of the objective function (multi-modality). In order to improve the reliability of the inversion results, we implemented and tested a two-step inversion scheme based on Genetic Algorithms (GAs). The proposed scheme performs several preliminary “parallel” runs (first step) and a final global run using the previously-determined fittest models as starting population.In this work we focus on the inversion of shear-wave velocity and layer thickness while fixing compressional-wave velocity and density according to user-defined Poisson's ratios and velocity–density relationship respectively. The procedure can nonetheless perform the inversion under different degrees of regularization, depending on the a priori information and the desired degree of freedom of the system.Thanks to the large number of considered models, in addition to the fittest model, a mean model and its accuracy are evaluated by means of a statistical approach based on the estimation of the Marginal Posterior Probability Density (MPPD).We tested the proposed GA-based inversion scheme on three synthetic models reproducing a complex structure with low-to-moderate velocity cover (also including a low-velocity channel) lying over hard bedrock. For all the considered cases the bedrock velocity and depth were properly identified, and velocity inversion was reconstructed with minor uncertainties.The performed tests also investigate the influence of the first higher mode, the reduction of the frequency range of the considered dispersion curve as well as the use of different number of strata. While a limited frequency range of the dispersion curve (maximum frequency reduced from 80 to 40 Hz) does not seem to significantly limit the accuracy of the retrieved model, the adoption of the correct number of strata and the addition of the first higher mode help better focus the final solution.In conclusion, the proposed approach represents an improvement of a purely GA-based optimization scheme and the MPPD-based mean model typically offers a more significant and precise solution than the fittest one.Results of the inversion performed on a field data set were validated by borehole stratigraphy. |
| |
| Keywords: | |
| 本文献已被 ScienceDirect 等数据库收录！ |
|
