Combining genetic and linearized algorithms for a two-step joint inversion of Rayleigh wave dispersion and H/V spectral ratio curves

SUMMARY The joint inversion of Rayleigh wave dispersion and H/V curves from environmental noise measurements allows the retrieval of S-wave velocity profiles for the shallow subsoil. For this purpose, genetic and linearized algorithm have been combined in a two-step inversion procedure, that allows the principal drawbacks typical of the application of each algorithm separately to be overcome. In the first step, a genetic algorithm procedure is used to constrain the subvolume of the parameter space where the absolute minimum of the misfit function is located. In the second step, a linearized inversion algorithm, having as an initial guess the minimum misfit model deduced from the first step, is applied to force the inversion towards the optimal solution. To evaluate the feasibility and effectiveness of this approach, seismic noise recordings at a test site in the Po river valley (North Italy) have been analysed. Here, detailed geophysical and geological information is available along with earthquake recordings, which allow a well constrained definition of both the local shear wave profile and transfer function. Comparisons between theoretical and experimental S-wave velocity profiles and, above all, between the theoretical and experimental site response functions shows that this combination of inversion procedures can very efficiently to manage the extreme non-linearity of the problem.

[1]  Stefano Parolai,et al.  Joint inversion of phase velocity dispersion and H/V ratio curves from seismic noise recordings using a genetic algorithm, considering higher modes , 2005 .

[2]  Hiroaki Yamanaka,et al.  Application of genetic algorithms to an inversion of surface-wave dispersion data , 1996, Bulletin of the Seismological Society of America.

[3]  Stratos Zacharopoulos,et al.  Guidelines for the implementation of the H/V spectral ratio technique on ambient vibrations measurements, processing and interpretation , 2004 .

[4]  M. Sambridge.  Geophysical inversion with a neighbourhood algorithm—I. Searching a parameter space , 1999 .

[5]  Kohji Tokimatsu,et al.  Effects of Multiple Modes on Rayleigh Wave Dispersion Characteristics , 1992 .

[6]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[7]  M. Sambridge,et al.  Monte Carlo analysis of inverse problems , 2002 .

[8]  Kohji Tokimatsu,et al.  S-Wave Velocity Profiling by Inversion of Microtremor H/V Spectrum , 2004 .

[9]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[10]  R. Wiggins,et al.  The general linear inverse problem - Implication of surface waves and free oscillations for earth structure. , 1972 .

[11]  S. Kramer.  Geotechnical Earthquake Engineering , 1996 .

[12]  Sesame Team,et al.  Sesame project - Deliverable D23.12 - WP12 - GUIDELINES FOR THE IMPLEMENTATION OF THE H/V SPECTRAL RATIO TECHNIQUE ON AMBIENT VIBRATIONS MEASUREMENTS, PROCESSING AND INTERPRETATION , 2004 .

[13]  David M. Boore,et al.  Rayleigh wave particle motion and crustal structure , 1969 .

[14]  Stefano Parolai,et al.  Joint inversion of H/V ratios and dispersion curves from seismic noise: Estimating the S‐wave velocity of bedrock , 2005 .

[15]  Donat Fäh,et al.  A theoretical investigation of average H/V ratios , 2001 .

[16]  R. D. Woods,et al.  Use of Rayleigh Modes in Interpretation of SASW Test , 1991 .

[17]  S. Parolai,et al.  S-wave Velocity Profiles for Earthquake Engineering Purposes for the Cologne Area (Germany) , 2006 .

[18]  A Robust Two-Step Inversion of Complex Magnetotelluric Apparent Resistivity Data , 2005 .

[19]  Frank Scherbaum,et al.  Determination of shallow shear wave velocity profiles in the Cologne, Germany area using ambient vibrations , 2003 .

[20]  M. Ohori,et al.  A Comparison of ESAC and FK Methods of Estimating Phase Velocity Using Arbitrarily Shaped Microtremor Arrays , 2002 .

[21]  T. Ohmachi,et al.  Ground Motion Characteristics Estimated from Spectral Ratio between Horizontal and Verticcl Components of Mietremors. , 1997 .

[22]  H. Pedersen,et al.  Constraints on absolute S velocities beneath the Aegean Sea from surface wave analysis , 2005 .

[23]  Jianghai Xia,et al.  Estimation of near‐surface shear‐wave velocity by inversion of Rayleigh waves , 1999 .

[24]  D. Albarello,et al.  QUICK ESTIMATES OF SOFT SEDIMENT THICKNESSES FROM AMBIENT NOISE HORIZONTAL TO VERTICAL SPECTRAL RATIOS: A CASE STUDY IN SOUTHERN ITALY , 2004 .

[25]  W. Menke.  Geophysical data analysis : discrete inverse theory , 1984 .

[26]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[27]  Kohji Tokimatsu,et al.  S-Wave Velocity Profiling by Joint Inversion of Microtremor Dispersion Curve and Horizontal-to-Vertical (H/V) Spectrum , 2005 .

[28]  David M. Boore,et al.  Comparing Shear-Wave Velocity Profiles From Inversion of Surface-Wave Phase Velocities with Downhole Measurements: Systematic Differences Between the CXW Method and Downhole Measurements at Six USC Strong-Motion Sites , 1998 .

[29]  Frank Scherbaum,et al.  How good are shear wave velocity models obtained from inversion of ambient vibrations in the Lower Rhine Embayment (N.W. Germany) , 2004 .

[30]  David E. Goldberg,et al.  Genetic Algorithms in Search Optimization and Machine Learning , 1988 .

[31]  Shear-wave velocity structure of sediments from the inversion of explosion-induced Rayleigh waves: Comparison with cross-hole measurements , 1997 .

[32]  Michel Campillo,et al.  Contribution of radar interferometry to a two-step inversion of the kinematic process of the 1992 Landers earthquake , 1999 .