Three-dimensional P- and S-wave velocity profiling of geotechnical sites using full-waveform inversion driven by field data

Abstract We discuss the application of a recently developed full-waveform-inversion-based technique to the imaging of geotechnical sites using field-collected data. Specifically, we address the profiling of arbitrarily heterogeneous sites in terms of P- and S-wave velocities in three dimensions, using elastic waves as probing agents. We cast the problem of finding the spacial distribution of the elastic soil properties as an inverse medium problem, directly in the time domain, and use perfectly-matched-layers (PMLs) to account for the semi-infinite extent of the site under investigation. After briefly reviewing the theoretical and computational aspects of the employed technique, we focus on the characterization of the George E. Brown Jr. Network for Earthquake Engineering Simulation site in Garner Valley, California (NEES@UCSB). We compare the profiles obtained from our full-waveform-inversion-based methodology against the profiling obtained from the Spectral-Analysis-of-Surface-Waves (SASW) method, and report agreement. In an attempt to validate our methodology, we also compare the recorded field data at select control sensors that were not used for the full-waveform inversion, against the response at the same sensors, computed based on the full-waveform-inverted profiles. We report very good agreement at the control sensors, which is a strong indicator of the correctness of the inverted profiles. Overall, the systematic framework discussed herein seems robust, general, practical, and promising for three-dimensional site characterization purposes.

[1]  G. Biros,et al.  An inverse problem formulation for parameter estimation of a reaction–diffusion model of low grade gliomas , 2014, Journal of mathematical biology.

[2]  Jeffrey C Bamber,et al.  Evaluation of the adjoint equation based algorithm for elasticity imaging , 2004, Physics in medicine and biology.

[3]  Richard D. Miller,et al.  Multichannel analysis of surface waves , 1999 .

[4]  Roland Martin,et al.  Improved forward wave propagation and adjoint-based sensitivity kernel calculations using a numerically stable finite-element PML , 2014 .

[5]  Omar Ghattas,et al.  Site characterization using full waveform inversion , 2013 .

[6]  Hiroaki Yamanaka,et al.  Two-dimensional shallow soil profiling using time-domain waveform inversion , 2015 .

[7]  S. Kucukcoban The inverse medium problem in PML-truncated elastic media , 2010 .

[8]  A. Fathi,et al.  Full-waveform inversion in three-dimensional PML-truncated elastic media , 2015, 1504.08340.

[9]  Dimitri Komatitsch,et al.  A hybrid method to compute short-period synthetic seismograms of teleseismic body waves in a 3-D regional model , 2013 .

[10]  Brittany D. Froese,et al.  Application of the Wasserstein metric to seismic signals , 2013, 1311.4581.

[11]  Michael McVay,et al.  Site characterization using Gauss–Newton inversion of 2-D full seismic waveform in the time domain , 2012 .

[12]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[13]  N. Petra,et al.  Model Variational Inverse Problems Governed by Partial Differential Equations , 2011 .

[14]  P. Mora Nonlinear two-dimensional elastic inversion of multioffset seismic data , 1987 .

[15]  Stephen J. Wright,et al.  Springer Series in Operations Research , 1999 .

[16]  Marc Bonnet,et al.  A Modified Error in Constitutive Equation Approach for Frequency-Domain Viscoelasticity Imaging Using Interior Data. , 2015, Computer methods in applied mechanics and engineering.

[17]  Craig M. Newtson,et al.  A Gauss–Newton full-waveform inversion for material profile reconstruction in viscoelastic semi-infinite solid media , 2016 .

[18]  Yi Wang,et al.  Three-dimensional full waveform inversion of short-period teleseismic wavefields based upon the SEM–DSM hybrid method , 2015 .

[19]  Loukas F. Kallivokas,et al.  Time‐domain hybrid formulations for wave simulations in three‐dimensional PML‐truncated heterogeneous media , 2015 .

[20]  Loukas F. Kallivokas,et al.  The inverse medium problem in heterogeneous PML-truncated domains using scalar probing waves , 2011 .

[21]  A. Oberai,et al.  Approaches to accommodate noisy data in the direct solution of inverse problems in incompressible plane strain elasticity , 2014, Inverse problems in science and engineering.

[22]  Loukas F. Kallivokas,et al.  Direct time-domain soil profile reconstruction for one-dimensional semi-infinite domains , 2009 .

[23]  F. Tröltzsch Optimal Control of Partial Differential Equations: Theory, Methods and Applications , 2010 .

[24]  John G. Proakis,et al.  Digital Signal Processing 4th Edition , 2006 .

[25]  Murthy N. Guddati,et al.  A novel wave equation solver to increase the efficiency of full waveform inversion , 2014 .

[26]  Andreas Fichtner,et al.  Multiscale full waveform inversion , 2013 .

[27]  Dimitri Komatitsch,et al.  High-resolution seismic array imaging based on an SEM-FK hybrid method , 2014 .

[28]  Assad A. Oberai,et al.  INVERSE PROBLEMS PII: S0266-5611(03)54272-1 Solution of inverse problems in elasticity imaging using the adjoint method , 2003 .

[29]  O. Ghattas,et al.  A Newton-CG method for large-scale three-dimensional elastic full-waveform seismic inversion , 2008 .

[30]  A Tikhonov,et al.  Solution of Incorrectly Formulated Problems and the Regularization Method , 1963 .

[31]  J. M. Roesset,et al.  Characterization of geotechical sites by SASW method , 1994 .

[32]  Ushnish Basu,et al.  Explicit finite element perfectly matched layer for transient three‐dimensional elastic waves , 2009 .

[33]  Arash Fathi Full-waveform inversion in three-dimensional PML-truncated elastic media : theory, computations, and field experiments , 2015 .

[34]  J. Tromp,et al.  Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements , 2011 .

[35]  Marc Bonnet,et al.  Large Scale Parameter Estimation Problems in Frequency-Domain Elastodynamics Using an Error in Constitutive Equation Functional. , 2013, Computer methods in applied mechanics and engineering.

[36]  C. Bunks,et al.  Multiscale seismic waveform inversion , 1995 .