Inversion of borehole weak motion records observed in Istanbul (Turkey)

SUMMARY The estimation of shear wave velocity and attenuation in near-surface geology is of primary importance in engineering seismology. In fact, such knowledge is essential for site response studies when preparing improved seismic hazard scenarios. In this study, we propose a linear inversion of the spectra of a deconvolved wavefield collected by a borehole array in Istanbul, Turkey. The spectra are calculated using as a reference the recordings collected by a sensor at the surface. This allows us to minimize the effect of the deconvolution filter on the peaks’ amplitudes. The feasibility of the proposed inversion scheme and the parametrization of the velocity and attenuation models were tested and assessed using synthetic data. The real data inversion is carried out using observation from three weak motion events. Several starting models are used for each event, changing the values of the quality factor Qs, whereas the starting Vs profiles are fixed to what was obtained by the analysis of the deconvolved wavefield in the time domain. Our results showed that fairly well-constrained Vs models can be obtained while the variability of the Qs profile, both between different inversions of the same data set, and between the inversions of different data sets, might be large. This indicates that Qs variations observed in models derived by strong motion data inversions should be considered with care. However, the high agreement of the estimated soil profiles in the uppermost layer, both when using different events and different starting models, suggests that any variation in the S-wave velocity that might be estimated when using strong motion data should be reliable.

[1]  J. Sheng,et al.  Interferometric/daylight seismic imaging , 2004 .

[2]  安芸 敬一 Space and time spectra of stationary stochastic waves, with special reference to microtremors , 1959 .

[3]  Richard L. Weaver,et al.  On the emergence of the Green's function in the correlations of a diffuse field: pulse-echo using thermal phonons. , 2001, Ultrasonics.

[4]  Atilla Ansal,et al.  The Ataköy vertical array (Turkey): Insights into seismic wave propagation in the shallow-most crustal layers by waveform deconvolution , 2009 .

[5]  Walter H. F. Smith,et al.  Free software helps map and display data , 1991 .

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

[7]  Wei Li,et al.  Site Amplification and Attenuation via Downhole Array Seismogram Inversion: A Comparative Study of the 2003 Miyagi-Oki Aftershock Sequence , 2008 .

[8]  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 .

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

[10]  David Halliday,et al.  Seismic interferometry, surface waves and source distribution , 2008 .

[11]  Jamison H. Steidl,et al.  Attenuation and Velocity Structure for Site Response Analyses via Downhole Seismogram Inversion , 2006 .

[12]  R. Weaver,et al.  On the emergence of the Green's function in the correlations of a diffuse field: pulse-echo using thermal phonons. , 2001, Ultrasonics.

[13]  Atilla Ansal,et al.  Seismic Input Motion Determined from a Surface–Downhole Pair of Sensors: A Constrained Deconvolution Approach , 2010 .

[14]  Atilla Ansal,et al.  Determination of shallow S-wave attenuation by down-hole waveform deconvolution: a case study in Istanbul (Turkey) , 2010 .

[15]  Roel Snieder,et al.  Extraction of near-surface properties for a lossy layered medium using the propagator matrix , 2007 .

[16]  Mario Bertero,et al.  Introduction to Inverse Problems in Imaging , 1998 .

[17]  H. Kawase,et al.  Evaluation of local site effects and their removal from borehole records observed in the Sendai region, Japan , 1995, Bulletin of the Seismological Society of America.

[18]  H. Okada The Microtremor Survey Method , 2003 .

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

[20]  Roel Snieder,et al.  Equivalence of the virtual-source method and wave-field deconvolution in seismic interferometry. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  Dominic Assimaki,et al.  A Wavelet-based Seismogram Inversion Algorithm for the In Situ Characterization of Nonlinear Soil Behavior , 2009 .

[22]  Michel Campillo,et al.  Emergence of broadband Rayleigh waves from correlations of the ambient seismic noise , 2004 .

[23]  Roel Snieder,et al.  Downhole Receiver Function: a Case Study , 2007 .

[24]  F. D. Martin,et al.  Nonlinear Soil Response of a Borehole Station Based on One-Dimensional Inversion during the 2005 Fukuoka Prefecture Western Offshore Earthquake , 2010 .

[25]  T. Satoh Inversion of QS of Deep Sediments from Surface-to-Borehole Spectral Ratios Considering Obliquely Incident SH and SV Waves , 2006 .

[26]  Rongjiang Wang,et al.  A simple orthonormalization method for stable and efficient computation of Green's functions , 1999, Bulletin of the Seismological Society of America.