Characterisation of ground motion recording stations in the Groningen gas field

The seismic hazard and risk analysis for the onshore Groningen gas field requires information about local soil properties, in particular shear-wave velocity (V_S). A fieldwork campaign was conducted at 18 surface accelerograph stations of the monitoring network. The subsurface in the region consists of unconsolidated sediments and is heterogeneous in composition and properties. A range of different methods was applied to acquire in situ V_S values to a target depth of at least 30 m. The techniques include seismic cone penetration tests (SCPT) with varying source offsets, multichannel analysis of surface waves (MASW) on Rayleigh waves with different processing approaches, microtremor array, cross-hole tomography and suspension P-S logging. The offset SCPT, cross-hole tomography and common midpoint cross-correlation (CMPcc) processing of MASW data all revealed lateral variations on length scales of several to tens of metres in this geological setting. SCPTs resulted in very detailed V_S profiles with depth, but represent point measurements in a heterogeneous environment. The MASW results represent V_S information on a larger spatial scale and smooth some of the heterogeneity encountered at the sites. The combination of MASW and SCPT proved to be a powerful and cost-effective approach in determining representative V_S profiles at the accelerograph station sites. The measured V_S profiles correspond well with the modelled profiles and they significantly enhance the ground motion model derivation. The similarity between the theoretical transfer function from the V_S profile and the observed amplification from vertical array stations is also excellent.     

Shear wave velocities of Mississippi embayment soils from low frequency surface wave measurements

Deep unconsolidated sediments in the Mississippi embayment will influence ground motions from earthquakes in the New Madrid seismic zone. Shear wave velocity profiles of these sediments are important input parameters for modeling wave propagation and site response in this region. Low-frequency, active-source surface wave velocity measurements were performed to develop small-strain shear wave velocity (V_S) profiles at eleven deep soil sites in the Mississippi embayment, from north of New Madrid, Missouri to Memphis, Tennessee. A servo-hydraulic, low-frequency source was used to excite surface wave energy to wavelengths of 600 m, resulting in V_S profiles to depths of over 200 m. The average V_S profile calculated from the eleven sites is in good agreement with common reference V_S profiles that have been used in seismic hazard studies of this region. The variability in V_S profiles is shown to be associated with changes in formation depth and thickness from site-to-site. Using lithologic information at each site, average formation velocities were developed and compared to previous studies. We found average V_S values of about 193 m/s for alluvial deposits, 400 m/s for the Upper Claiborne formations, and 685 m/s for the Memphis Sand formation.     

Three‐dimensional VS profiling using microtremors in Kushiro,Japan

A practical method is presented for determining three‐dimensional S‐wave velocity (V_S) profile from microtremor measurements. Frequency–wave number (f–k) spectral analyses of microtremor array records are combined, for this purpose, with microtremor horizontal‐to‐vertical (H/V) spectral ratio techniques. To demonstrate the effectiveness of the proposed method, microtremor measurements using arrays of sensors were conducted at six sites in the city of Kushiro, Japan. The spectral analyses of the array records yield dispersion characteristics of Rayleigh waves and H/V spectra of surface waves, and joint inversion of these data results in V_S profiles down to bedrock at the sites. Conventional microtremor measurements were performed at 230 stations within Kushiro city, resulting in the H/V spectra within the city. Three‐dimensional V_S structure is then estimated from inversion of the H/V spectra with the V_S values determined from the microtremor array data. This reveals three‐dimensional V_S profile of Kushiro city, together with an unknown hidden valley that crosses the central part of the city. The estimated V_S profile is consistent with available velocity logs and results of subsequent borings, indicating the effectiveness of the proposed method. Copyright © 2008 John Wiley & Sons, Ltd.     

Determination of S-wave velocity structure using microtremors and spac method applied in Thessaloniki (Greece)

Array measurements of microtremors at 16 sites in the city of Thessaloniki were performed to estimate the V_s velocity of soil formations for site effect analysis. The spatial autocorrelation method was used to determine phase velocity dispersion curves in the frequency range from 0.8–1.5 to 6–7 Hz. A Rayleigh wave inversion technique (stochastic method) was subsequently applied to determine the Vs profiles at all the examined sites. The determination of V_s profiles reached a depth of 320 m. Comparisons with V_s values from cross-hole tests at the same sites proved the reliability of the SPAC method. The accuracy of the V_s profiles, the ability to reach large penetration depths in densely populated urban areas and its low cost compared to conventional geophysical prospecting, make Mictrotremor Exploration Method very attractive and useful for microzonation and site effects studies. An example of its application for the site characterization in Thessaloniki is presented herein.     

Comparing shear-wave velocity profiles inverted from multichannel surface wave with borehole measurements

Recent field tests illustrate the accuracy and consistency of calculating near-surface shear (S)-wave velocities using multichannel analysis of surface waves (MASW). S-wave velocity profiles (S-wave velocity vs. depth) derived from MASW compared favorably to direct borehole measurements at sites in Kansas, British Columbia, and Wyoming. Effects of changing the total number of recording channels, sampling interval, source offset, and receiver spacing on the inverted S-wave velocity were studied at a test site in Lawrence, Kansas. On the average, the difference between MASW calculated V_s and borehole measured V_s in eight wells along the Fraser River in Vancouver, Canada was less than 15%. One of the eight wells was a blind test well with the calculated overall difference between MASW and borehole measurements less than 9%. No systematic differences were observed in derived V_s values from any of the eight test sites. Surface wave analysis performed on surface data from Wyoming provided S-wave velocities in near-surface materials. Velocity profiles from MASW were confirmed by measurements based on suspension log analysis.     

An overview of surface wave methods and a reliability study of a simplified inversion technique

Following a brief overview of the history and the development of the Surface Wave Method—with a focus on techniques for processing and inverting field data—a Simplified Inversion Method (SIM) is described, which constitutes an improvement of the Satoh et al. (1991) [1] method. The SIM is a direct inversion method of surface wave dispersion data, making use of a penetration depth coefficient, a_R, whose value is a function of Poisson's ratio and the overall shape of the dispersion curve. In the present study the coefficient a_R has been evaluated using data from (a) an extensive database compiled from the technical literature and containing results of inverted surface wave measurements and nearby cross-hole/down-hole measurements, (b) results of side by side surface wave and cross-hole measurements, performed at five sites in the course of this study, (c) finite element analyses simulating the performance of surface wave measurements and thus providing “virtual” data, and (d) applying a current advanced inversion code, available on the Web. Based on all the above data, optimum values of a_R (and of the corresponding uncertainty of the derived V_so vs. depth profile) have been estimated. These values were found to be independent of depth from ground surface. The results of all analyses and comparisons indicate that for the majority of realistic soil profiles (including cases of normal and inverse dispersion conditions) the proposed SIM provides very reliable V_so vs. depth profiles when a value of a_R=0.63–0.67 is used in the inversion process. It is concluded that the SIM can be used with confidence as a direct inversion method of surface wave dispersion data.     

An overview of surface wave methods and a reliability study of a simplified inversion technique

Following a brief overview of the history and the development of the Surface Wave Method—with a focus on techniques for processing and inverting field data—a Simplified Inversion Method (SIM) is described, which constitutes an improvement of the Satoh et al. (1991) [1] method. The SIM is a direct inversion method of surface wave dispersion data, making use of a penetration depth coefficient, a_R, whose value is a function of Poisson's ratio and the overall shape of the dispersion curve. In the present study the coefficient a_R has been evaluated using data from (a) an extensive database compiled from the technical literature and containing results of inverted surface wave measurements and nearby cross-hole/down-hole measurements, (b) results of side by side surface wave and cross-hole measurements, performed at five sites in the course of this study, (c) finite element analyses simulating the performance of surface wave measurements and thus providing “virtual” data, and (d) applying a current advanced inversion code, available on the Web. Based on all the above data, optimum values of a_R (and of the corresponding uncertainty of the derived V_so vs. depth profile) have been estimated. These values were found to be independent of depth from ground surface. The results of all analyses and comparisons indicate that for the majority of realistic soil profiles (including cases of normal and inverse dispersion conditions) the proposed SIM provides very reliable V_so vs. depth profiles when a value of a_R=0.63–0.67 is used in the inversion process. It is concluded that the SIM can be used with confidence as a direct inversion method of surface wave dispersion data.     

Comparing liquefaction evaluation methods using penetration-VS relationships

Three methods that follow the general format of the Seed-Idriss simplified procedure for evaluating liquefaction resistance of soils are compared in this paper. They are compared by constructing relationships between penetration resistance and small-strain shear–wave velocity (V_S) implied from cyclic resistance ratio (CRR) curves for the three methods, and by plotting penetration-V_S data pairs. The penetration-V_S data pairs are from 43 Holocene-age sand layers in California, South Carolina, Canada, and Japan. It is shown that the V_S-based CRR curve is more conservative than CRR curves based on the Standard Penetration Test (SPT) and Cone Penetration Test (CPT), for the compiled Holocene data. This result agrees with the findings of a recent probability study where the SPT-, CPT-, and V_S-based CRR curves were characterized as curves with average probability of liquefaction of 31, 50, and 26%, respectively. New SPT- and CPT-based CRR equations are proposed that provide more consistent assessments of liquefaction potential for the Holocene sand layers considered.     

Pre-loading effect on dynamic soil properties: Seismic methods and their efficiency in geotechnical aspects

There are many publications on the investigation of soil properties using seismic prospecting. Among these properties, special attention has been given to shear wave velocity V_S, using more than two different methods for soil and site characterization. In this study, the in-hole, non-invasive refraction and surface wave inversion methods to evaluate soil improvement are investigated. The investigation was conducted on the new Egnatia highway (Northern Greece). Wave velocity profiles have been measured before and after preloading for the construction of an embankment at a soft soil site. The purpose is to quantify the dynamic properties and to evaluate the efficiency of the applied tools in detecting their variation. Among others, an emphasis was given to the observed improvement at particular layers of high sand content.     

Characterization of site conditions (soil class,V<Subscript>S30</Subscript>, velocity profiles) for 33 stations from the French permanent accelerometric network (RAP) using surface-wave methods

Data provided by accelerometric networks are important for seismic hazard assessment. The correct use of accelerometric signals is conditioned by the station site metadata quality (i.e., soil class, V_S30, velocity profiles, and other relevant information that can help to quantify site effects). In France, the permanent accelerometric network consists of about 150 stations. Thirty-three of these stations in the southern half of France have been characterized, using surface-wave-based methods that allow derivation of velocity profiles from dispersion curves of surface waves. The computation of dispersion curves and their subsequent inversion in terms of shear-wave velocity profiles has allowed estimation of V_S30 values and designation of soil classes, which include the corresponding uncertainties. From a methodological point of view, this survey leads to the following recommendations: (1) perform both active (multi-analysis surface waves) and passive (ambient vibration arrays) measurements to derive dispersion curves in a broadband frequency range; (2) perform active acquisitions for both vertical (Rayleigh wave) and horizontal (Love wave) polarities. Even when the logistic contexts are sometimes difficult, the use of surface-wave-based methods is suitable for station-site characterization, even on rock sites. In comparison with previous studies that have mainly estimated V_S30 indirectly, the new values here are globally lower, but the EC8-A class sites remain numerous. However, even on rock sites, high frequency amplifications may affect accelerometric records, due to the shallow relatively softer layers.     

四川、甘肃地区VS30经验估计研究

目前我国建筑工程抗震设计规范中对于工程场地条件的判断依据主要是地表以下20m深度范围内土层的等效剪切波速,简称VS20。相比之下,国外应用较广的是地表以下30m深度范围内的等效剪切波速,简称VS30。这种差别导致国内科研工作者在应用国外的地震工程、工程抗震模型时经常遇到对场地条件描述不准确的困难。为了解决这个问题,本文根据147个四川、甘肃地区国家强震动台站20m左右深度的钻孔剪切波速数据,利用延拓方法、场地分类统计方法以及基于地形特征的VS30估计方法研究各台站VS30与VS20的经验关系,对比发现基于速度梯度延拓的结果最为可取。参考国际上通用的Geomatrix Classification场地分类标准,最终得到四川、甘肃地区各类场地的平均VS30,此结果可以为缺乏钻孔数据的工程场地的VS30估计提供参考。     

Evaluation of shear wave velocity profile using SPT based uphole method

The uphole method is a field seismic test which uses receivers on the ground surface and an underground source. A modified form of the uphole method is introduced in order to obtain efficiently the shear wave velocity (V_S) profile of a site. This method is called the standard penetration test (SPT)-uphole method because it uses the impact energy of the split spoon sampler in the SPT test as a source. Since the SPT-uphole method can be performed simultaneously with the SPT test it is economical and not labor intensive compared to the original uphole methods which use small explosives or a mechanical source. Field testing and interpretation procedures for the proposed method are described. To obtain reliable travel time information of the shear wave, the first peak point of the shear wave using two component geophones is recommended. Through a numerical study using the finite element method (FEM), the procedure of the proposed method was verified. Finally, the SPT-uphole method was performed at several sites, and the field applicability of the proposed method was verified by comparing the V_S profiles determined by the SPT-uphole method with the profiles determined by the downhole, the spectral analysis of surfaces waves (SASW) method and from the SPT-N values.     

Upper lithospheric structure in the central Fennoscandian shield: Constraints from P- and S-wave velocity models and V<Subscript>P</Subscript>/V<Subscript>S</Subscript> ratio distribution of the BALTIC wide-angle seismic profile

The paper presents an analysis of the crust and upper mantle structure in the central Fennoscandian shield based on new P- and S-wave 2D velocity models of the BALTIC wide-angle reflection and refraction profiles. Using reprocessing of the old data, new P- and S-wave velocity models and V _P /V _S ratio distribution were developed. Moving from SW to NE, the thickness of the crust varies strongly, from ∼36 km to extremely thick, 58–64 km, crossing Wiborg rapakivi massif, Saimaa and Outokumpu areas, and Eastern Finland complex. Based on the lateral variations of V _P , V _P /V _S and thickness of the crust, three main blocks of the crust and upper mantle were distinguished from SW to NE: southwestern, associated with Wiborg rapakivi massif; the central, having the highest thickness of the crust; and the northeastern, not well documented, with Archaean basement.     

Simultaneous inversion of velocity structures and hypocentral locations: Application to the Friuli seismic area NE Italy

A layeredP- andS-wave velocity model is obtained for the Friuli seismic area using the arrival time data ofP- andS-waves from local earthquakes. A damped least-squares method is applied in the inversion.The data used are 994P-wave arrival times for 177 events which have epicenters in the region covered by the Friuli seismic network operated by Osservatorio Geofisico sperimentale (OGS) di Trieste, which are jointly inverted for the earthquake hypocenters andP-wave velocity model. TheS-wave velocity model is estimated on the basis of 978S-wave arrival times and the hypocenters obtained from theP-wave arrival time inversion. We also applied an approach thatP- andS-wave arrival time data are jointly used in the inversion (Roecker, 1982). The results show thatS-wave velocity structures obtained from the two methods are quite consistent, butP-wave velocity structures have obvious differences. This is apparent becauseP-waves are more sensitive to the hypocentral location thanS-waves, and the reading errors ofS-wave arrival times, which are much larger than those ofP-waves, bring large location errors in the joint inversion ofP- andS-wave arrival time. The synthetic data tests indicated that when the reading errors ofS-wave arrivals are larger than four times that ofP-wave arrivals, the method proposed in this paper seems more valid thanP- andS-wave data joint inversion. Most of the relocated events occurred in the depth range between 7 and 11 km, just above the biggest jump in velocity. This jump might be related to the detachment line hypothesized byCarulli et al. (1982). From the invertedP- andS-wave velocities, we obtain an average value 1.82 forV _p /V _s in the first 16 km depth.     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.     

Soil characterization of Tınaztepe region (İzmir/Turkey) using surface wave methods and nakamura (HVSR) technique

To determine the shear wave velocity structure and predominant period features of T?naztepe in ?zmir, Turkey, where new building sites have been planned, active–passive surface wave methods and single-station microtremor measurements are used, as well as surface acquisition techniques, including the multichannel analysis of surface waves (MASW), refraction microtremor (ReMi), and the spatial autocorrelation method (SPAC), to pinpoint shallow and deep shear wave velocity. For engineering bedrock (V _s > 760 m/s) conditions at a depth of 30 m, an average seismic shear wave velocity in the upper 30 m of soil (AV_s30) is not only accepted as an important parameter for defining ground behavior during earthquakes, but a primary parameter in the geotechnical analysis for areas to be classified by V _s30 according to the National Earthquake Hazards Reduction Program (NEHRP). It is also determined that Z1.0, which represents a depth to V _s = 1000 m/s, is used for ground motion prediction and changed from 0 to 54 m. The sediment–engineering bedrock structure for T?naztepe that was obtained shows engineering bedrock no deeper than 30 m. When compared, the depth of engineering bedrock and dominant period map and geology are generally compatible.     

Dynamic Site Characterization and Correlation of Shear Wave Velocity with Standard Penetration Test ‘N’ Values for the City of Agartala,Tripura State,India

Seismic site characterization is the basic requirement for seismic microzonation and site response studies of an area. Site characterization helps to gauge the average dynamic properties of soil deposits and thus helps to evaluate the surface level response. This paper presents a seismic site characterization of Agartala city, the capital of Tripura state, in the northeast of India. Seismically, Agartala city is situated in the Bengal Basin zone which is classified as a highly active seismic zone, assigned by Indian seismic code BIS-1893, Indian Standard Criteria for Earthquake Resistant Design of Structures, Part-1 General Provisions and Buildings. According to the Bureau of Indian Standards, New Delhi (2002), it is the highest seismic level (zone-V) in the country. The city is very close to the Sylhet fault (Bangladesh) where two major earthquakes (M _w > 7) have occurred in the past and affected severely this city and the whole of northeast India. In order to perform site response evaluation, a series of geophysical tests at 27 locations were conducted using the multichannel analysis of surface waves (MASW) technique, which is an advanced method for obtaining shear wave velocity (V _s) profiles from in situ measurements. Similarly, standard penetration test (SPT-N) bore log data sets have been obtained from the Urban Development Department, Govt. of Tripura. In the collected data sets, out of 50 bore logs, 27 were selected which are close to the MASW test locations and used for further study. Both the data sets (V _s profiles with depth and SPT-N bore log profiles) have been used to calculate the average shear wave velocity (V _s30) and average SPT-N values for the upper 30 m depth of the subsurface soil profiles. These were used for site classification of the study area recommended by the National Earthquake Hazard Reduction Program (NEHRP) manual. The average V _s30 and SPT-N classified the study area as seismic site class D and E categories, indicating that the city is susceptible to site effects and liquefaction. Further, the different data set combinations between V _s and SPT-N (corrected and uncorrected) values have been used to develop site-specific correlation equations by statistical regression, as ‘V _s’ is a function of SPT-N value (corrected and uncorrected), considered with or without depth. However, after considering the data set pairs, a probabilistic approach has also been presented to develop a correlation using a quantile–quantile (Q–Q) plot. A comparison has also been made with the well known published correlations (for all soils) available in the literature. The present correlations closely agree with the other equations, but, comparatively, the correlation of shear wave velocity with the variation of depth and uncorrected SPT-N values provides a more suitable predicting model. Also the Q–Q plot agrees with all the other equations. In the absence of in situ measurements, the present correlations could be used to measure V _s profiles of the study area for site response studies.     

Analysis of Short-Period Rayleigh Waves Recorded in the Bohemian Massif Area During Celebration 2000 Experiment

The aim of this paper is to show the application of short-period surface waves recorded during deep seismic sounding experiment for constraining shallow velocity structure of the crust. Phase velocity of fundamental mode Rayleigh waves, observed along the CELEBRATION 2000 experiment profile CEL09, were obtained by a p-ω method and has been subsequently inverted for one-dimensional shear velocity models for the top 2 km. Multiple filter technique applied to one shot gather was used to carry out a joint inversion of phase and group velocity data and to provide γ_R data to be used for Q_β inversion. Validity of obtained V_S and Q_β models was confirmed by the reflectivity method. Noticeably, no clear dispersive wawes were observed in the Tepla-Barrandian Unit. Quasi-2D model based on the individual 1D V_S models is well correlated with the surface geology. Lower V_S are observed in the Saxothuringian Zone in comparison to the Moldanubian Zone. In the vicinity of the Central Bohemian and Moldanubian Plutons, the near-surface V_S values are relatively low, but below 1 km depth, they are higher than in surrounding areas. We interpret it as the result of the weathering and cracks within the granitoid rocks.     

Anisotropic 4D seismic response inferred from ultrasonic laboratory measurements: A direct comparison with the isotropic response

Pore-pressure depletion causes changes in the triaxial stress state. Pore-pressure depletion in a flat reservoir, for example, can be reasonably approximated as uniaxial compaction, in which the horizontal effective stress change is smaller than the vertical effective stress. Furthermore, the stress sensitivity of velocities can be angle-dependent. Therefore, time-lapse changes in reservoir elastic anisotropy are expected as a consequence of production, which can complicate the interpretation of the 4D seismic response. The anisotropic 4D seismic response caused by pore-pressure depletion was investigated using existing core velocity measurements. To make a direct comparison between the anisotropic 4D seismic response and the isotropic response based only on vertical velocities, pseudoisotropic elastic properties were utilized, and the two responses were compared in terms of a dynamic rock physics template. A comparison of the dynamic rock physics templates indicates that time-lapse changes in reservoir elastic anisotropy have a noticeable impact on the interpretation of 4D seismic data. Changes in anisotropy as a result of pore-pressure depletion cause a time-lapse amplitude variation with offset response as if there is a reduction in V_P/V_S (i.e., pseudoisotropic V_P/V_S decreases), although the vertical V_P/V_S increases. The impact of time-lapse changes in anisotropy on the amplitude variation with offset gradient was also investigated, and the time-lapse anisotropy was found to enhance changes in the amplitude variation with offset gradient for a given case.     

Determination of shale stiffness tensor from standard logs

A technique allowing inversion of the shale stiffness tensor from standard logging data: sonic velocities, density, porosity and clay content is developed. The inversion is based on the effective medium theory. The testing of the technique on laboratory measurements of the elastic wave velocities in shale samples shows that the inversion makes it possible to predict the elastic wave velocities V_P, V_S1 and V_S2 in any direction within an error of a few per cent. The technique has been applied for the stiffness tensor inversion along a well penetrating a shale formation of the Mississippian age altered by thin layers of limestone. It is demonstrated that the symmetry of a stiffness tensor inverted at the sonic frequency (2 kHz) is slightly orthorhombic and taking into account the experimental errors, can be related to the vertical transverse isotropy symmetry. For the productive interval of the shale formation, the Thomsen parameters ?, γ, and δ average, respectively, 0.32, 0.25 and 0.21, which indicate anelliptic behaviour of the velocities in this shale. The coefficients of anisotropy of this shale interval are around 24% and 20% for the compressional and shear waves, respectively. The values of the inverted velocities in the bedding plane for this interval are in good agreement with the laboratory measurements. The technique also allows inversion of the water saturation of the formation (Sw) and the inverted values are in agreement with the Sw values available for this formation. A Backus‐like upscaling of the inverted stiffness tensors is carried out for the lower and upper bounds of the frequency band used in the crosswell tomography (100 Hz and 500 Hz). These results can serve as an initial velocity model for the microearthquake location during hydrofracking of the shale formation.