Geospatial contour mapping of shear wave velocity for Mumbai city

Shear wave velocity is one of the most important input parameter in the analysis of geotechnical earthquake engineering problems, particularly to estimate site-specific amplification factor and ground response study. Dynamic in situ tests such as spectral analysis of surface waves (SASW) or multichannel analysis of surface waves (MASW) are very expensive. Also due to lack of specialized personnel, these tests are generally avoided in many soil investigation programs. Worldwide, several researchers have developed correlations between the SPT ‘N’ value and shear wave velocity ‘V _s’, which are useful for determining the dynamic soil properties. In the present study, more than 400 numbers of soil borehole data were collected from various geotechnical investigation agencies, government engineering institutes and geotechnical laboratories from different parts of Mumbai city, which is financial capital of India with highest population density. In this paper, an attempt has been made to develop the correlation between the SPT ‘N’ value and shear wave velocity ‘V _s’ for various soil profile of Mumbai city and compared with other existing correlations for different cities in India. Using Geographical Information System (GIS), a geospatial contour map of shear wave velocity profile for Mumbai city is prepared with contour intervals of 25 and 50 m/s. The scarcity of database or maps of shear wave velocity profile for Mumbai city will make the present geospatial contour maps extremely useful and beneficial to the designer, practitioners for seismic hazard study involved in geotechnical earthquake engineering.     

Use of Surface Waves in Statistical Correlations of Shear Wave Velocity and Penetration Resistance of Chennai Soils

Shear wave velocity (V _s) is one of the most important input parameter to represent the stiffness of the soil layers. It is preferable to measure V _s by in situ wave propagation tests, however it is often not economically feasible to perform the tests at all locations. Hence, a reliable correlation between V _s and standard penetration test blow counts (SPT-N) would be a considerable advantage. This paper presents the development of empirical correlations between V _s and SPT-N value for different categories of soil in Chennai city characterized by complex variation of soil conditions. The extensive shear wave velocity measurement was carried out using Multichannel Analysis of Surface Waves (MASW) technique at the sites where the SPT-N values are available. The bender element test is performed to compare the field MASW test results for clayey soils. The correlations between shear wave velocity and SPT-N with and without energy corrections were developed for three categories of soil: all soils, sand and clay. The proposed correlations between uncorrected and energy corrected SPT-N were compared with regression equations proposed by various other investigators and found that the developed correlations exhibit good prediction performance. The proposed uncorrected and energy corrected SPT-N relationships show a slight variation in the statistical analysis indicating that both the uncorrected and energy corrected correlations can predict shear wave velocity with equal accuracy. It is also found that the soil type has a little effect on these correlations below SPT-N value of about 10.     

Utilization of sample disturbance for dating a marl deposit

Summary Values of the low-amplitude shear modulus, G _o, of a marl underlying the city of Kalamata, Greece, obtained by in situ crosshole and laboratory resonant column techniques are presented in the paper. It was found that the laboratory values were lower than their field counterparts and the difference was considered to be the result of (i) the temporary stress release following sampling and (ii) the mechanical disturbance associated with specimen preparation. By utilizing previous experimental results relating the field versus laboratory difference in the value of G _o to the age of the clay deposit, it became possible to estimate the age of the particular marl which was found to agree with the geologically expected age of the deposit. The approach taken in this study may, upon further improvement, become a useful technique for dating geological materials, especially those possessing cohesion.     

Experimental evaluation of shear strength behaviour of plastic and non-plastic silts

Silt is available in many parts of the world in combination with sands and clays. However, due to lack of clear understanding of its engineering behaviour, most of the time it is interpreted in terms of either sands or clays. Structures that are usually built on silty soils are designed to take into account design procedures developed for sandy or clayey soils. Presence of silts in combination with varying amount of sand and clays produces silt that is either plastic or non-plastic in nature. Silt is available in and around the Delhi region, in a majority mixture along with fine sands, which is non-plastic in nature. On the other hand silty deposits found in offshore Bombay High region are found in abundance along with significant amounts of clays and are termed as plastic silts. In this paper a comparison of the stress-strain behaviour of plastic and non-plastic silts is carried out under triaxial compression loading during both drained and undrained conditions. Two representative samples each from Delhi and Bombay High regions were considered for this comparison and results of stress-strain under four sets of confining pressure are discussed in detail. It is observed from this study that behaviour of silts is mainly dependent on the composition and structure of the resultant soil matrix. It is concluded from the results that shear strength parameters as well as volume change/pore pressure response of silty soils is dominated by the constituent soil present along with the silt. It is seen from the comparative behaviour of non-plastic and plastic silts that the presence of sand and clays has a governing effect on pore pressure development and the resultant friction angle. The study also corroborated that the nature of silt is transitional both in the case of plastic and non-plastic forms.     

Seismic piezocone interpretation for shear wave velocity (<Emphasis Type="Italic">V</Emphasis><Subscript>S</Subscript>) determination in the Persian Gulf

Shear wave velocity is one of the important factors representing the dynamic characteristics of soil layers. Hence, many researchers have focused their studies on determining shear wave velocity by direct field measurements or expressions developed by other soil parameters. The shear module and damping ratio of the soil layers also play a similar role in the majority of dynamic soil response analyses. Nevertheless, since they have to be measured in the laboratory by resonant column or cyclic triaxial tests on undisturbed samples, the possibility of preparing such samples and the reliability of the obtained results are of great concerns. In the present study, great effort has been made to determine the above dynamic factors by means of field data obtained from a versatile instrument, namely the seismic piezocone (SPCT_U), and to derive expressions correlating them with some parameters obtainable by much simpler instruments. The reliability of laboratory measurements on undisturbed samples is also evaluated. The seismic piezocone test apparatus has been employed to evaluate the soil properties at 1-m depth intervals by means of measuring tip resistance, sleeve resistance, pore pressure and shear wave velocity. The shear module and the damping ratio are calculated using field data. Meanwhile, in order to assess the laboratory measurements of these parameters, some resonant column tests and cyclic triaxial tests on undisturbed samples of the same soil layers have been carried out. In order to compare the field results of shear modulus and damping ratios with those obtained from laboratory tests, the influences of the soil nature and sample disturbance on the conventional laboratory methods are evaluated and discussed. The shear wave velocity is correlated to overburden pressure and the corrected tip resistance for two groups of fine soils, namely silty clays and carbonate clayey silts, which mainly cover the areas under study in this project, are located in southern parts of Iran near the Persian Gulf. According to the results of the present study, there are narrow limits of shear modulus regarding soils for which the laboratory tests and the field measurements yield approximately the same shear modulus. This limit of shear modulus is about 30–50(MPa) for clay deposits and 70–100 (MPa) for sandy deposits. Also the shear wave velocity can be calculated by a simple expression from total overburden pressure and the tip resistance of simple cone penetration test results conventionally available in many soil explorations prior to engineering practices. However, if the pore pressure inside the saturated soil deposits can be measured by a piezocone apparatus, the shear wave velocity may be calculated using another suggested equation in terms of effective overburden pressure in the present study. Regarding the shear module and the damping ratio, due to the disturbances of the stiff deposits in the sampling process and great deviations of laboratory results from field results, the laboratory measurements of these parameters out of the above limits are not recommended.     

干砂弹性参数测定的弯曲-伸展元试验

利用安装在共振柱测试系统中的弯曲-伸展元,开展了干砂中P波（压缩波）和S波（剪切波）的室内试验,详细地分析了干砂中P波和S波的信号特征,研究了输入频率、土体密实度和有效围压对输出信号的影响。对比各种信号分析方法,并参考共振柱试验结果确定了S波的传播时间。根据实测波速和波动理论,确定了土体的弹性参数,包括剪切模量,侧限模量和泊松比。研究结果表明,P波和S波的输出信号频率在一定程度上随输入信号频率、土体密实度和有效围压的增加而增加,且P波信号比S波信号更容易确定波的传播时间;土体的弹性模量随土体密实度和有效围压的增加而增加,但剪切模量增长比侧限模量快;土体的泊松比并非一个常数,随着土体密实度和有效围压的增加而下降。初步探讨了利用剪切模量估算泊松比,以方便实际工程应用。     

Variations in shear wave velocity and soil site class in Kolkata city using regression and sensitivity analysis

The detrimental effects of an earthquake are strongly influenced by the response of soils subjected to dynamic loading. The behavior of soils under dynamic loading is governed by the dynamic soil properties such as shear wave velocity, damping characteristics and shear modulus. Worldwide, it is a common practice to obtain shear wave velocity (V _s in m/s) using the correlation with field standard penetration test (SPT) N values in the absence of sophisticated dynamic field test data. In this paper, a similar but modified advanced approach has been proposed for a major metro city of eastern India, i.e., Kolkata city (latitudes 22°20′N–23°00′N and longitudes 88°04′E–88°33′E), to obtain shear wave velocity profile and soil site classification using regression and sensitivity analyses. Extensive geotechnical borehole data from 434 boreholes located across 75 sites in the city area of 185 km² and laboratory test data providing information on the thickness of subsoil strata, SPT N values, consistency indices and percentage of fines are collected and analyzed thoroughly. A correlation between shear wave velocity (V _s) and SPT N value for various soil profiles of Kolkata city has been established by using power model of nonlinear regression analysis and compared with existing correlations for other Indian cities. The present correlations, having regression coefficients (R ²) in excess of 0.96, indicated good prediction capability. Sensitivity analysis predicts that significant influence of soil type exists in determining V _s values, for example, typical silty sand shows 30.4 % increase in magnitude of V _s as compared to silt of Kolkata city. Moreover, the soil site classification shows Class D and Class E category of soil that exists typically in Kolkata city as per NEHRP (Recommended provisions for seismic regulations for new buildings and other structures—Part 1: Provisions. Prepared by the Building Seismic Safety Council for the Federal Emergency Management Agency (Report FEMA 450), Washington, DC, 2003) guidelines and thereby highlighting the seismic vulnerability of the city. The results presented in this study can be utilized for seismic microzonation, ground response analysis and hazard assessment for Kolkata city.

     

A non‐coaxial critical state soil model and its application to simple shear simulations

The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd.     

基于SCPTU的软土最大剪切模量测试分析研究

最大剪切模量是土的基本力学特性参数,对土动力特性分析和岩土工程抗震设计有着重要的意义。目前最大剪切模量的确定主要依赖于室内试验,土样的采取和室内试验既耗时又不能保证精度。以江苏北部里下河古泻湖相软土地区高速公路建设为工程背景,采用地震波孔压静力触探（SCPTU）和下孔法（DHT）进行了土层剪切波速的测试,基于SCPTU和DHT剪切波速测试资料建立了最大剪切模量Gmax与SCPTU测试参数之间的相关关系。结果表明：采用SCPTU测试的锥尖阻力和孔压参数能够用于对软土的最大剪切模量的初步评价。     

Stratigraphic delineation by three-dimensional clustering of piezocone data

Piezocone soundings are a fast and economical approach for geotechnical site characterization, providing three separate and continuous channels of data with depth, including: tip resistance q _T, porewater pressure u ₂ and sleeve friction f _s. Literally hundreds to thousands of data points are collected by a single sounding. Since these readings are functions of both soil type and soil behaviour, they can be used for the delineation of soil stratigraphy.

One way to process large amounts of data involves clustering. Cluster analysis is an efficient statistical way to analyse the stratigraphic vertical profiling of geomaterials and means to detect the inherent similarity between data sets and group them together. Clustering in previous geotechnical research was based on only two channels of piezocone data (q _T and u₂). The method works well for soils that are under the groundwater table and was applied to soundings in clay deposits.

In the present paper, a new cluster analysis approach is developed based on all three channels of data, thus extending the method to soils above the water table and applicable to sands, silts, and clays. Example soil profiles derived by three-channel cluster analysis are presented herein and compared with conventional soil boring and sampling data.     

Dynamic soil properties for microzonation of Delhi,India

Delhi, the capital of India, has experienced mild seismic shaking during several earthquakes in the past. The large variations of depth to bedrock and ground water table coupled with different soil types at different locations of Delhi necessitate a seismic microzonation study. Dynamic soil properties such as shear wave velocity, modulus reduction and damping characteristics of local soils are the basic and essential input parameters for conducting even a preliminary ground response analysis which is an essential input in microzonation studies. Shear wave velocity is not measured routinely due to its high cost and lack of the required expertise. Several researchers in the past developed correlations between shear wave velocity (V _s ) and routinely measured N values. In the present study, shear wave velocity profiles measured in the field at more than 80 borehole locations to a depth of about 20 to 32m using Spectral Analysis of Surface Waves (SASW) are presented and correlations between shear wave velocity and N values are also presented for use by engineers and designers. Results of strain and stress controlled cyclic triaxial tests on remoulded samples of sand-silt mixtures in the high strain range are used for generating the modulus reduction and damping curves and are compared with the well-known curves in the literature. The results presented in this article can be used for microzonation studies as well as site specific ground response analyses at Delhi.     

Undrained Shear Strength in Cohesive Soils Estimated by Directional Modes of In-Situ Shear Wave Velocity

The estimated undrained shear strength (s_u) is often not a unique value because it can be evaluated by various test types and/or procedures, such as different failure modes, shear strain rates, and boundary conditions. This study explores (1) the relationship between reference undrained shear strength and in situ shear wave velocity in terms of the effective overburden stress, and (2) the independent relationships to evaluate the undrained shear strength with special consideration of different directional and polarization modes (VH, HV, HH shear waves), which has not been reported. This evaluation is done via a worldwide database compiled from 43 well-documented geotechnical test sites associated with soft ground. Finally, new correlation models are proposed to estimate the undrained shear strength based on the in situ shear wave velocity as well as the plasticity index or the overconsolidation ratio. The application of the shear wave velocity–undrained shear strength relation is illustrated through two independent case studies. The proposed relationships are expected to contribute to reasonable estimates of undrained shear strength as well as offer practical guidance on even extrapolation beyond the data that is available to geotechnical engineers.     

Formulation of a unified constitutive model for clays and sands

This paper presents a new generalized effective stress model, referred to as MIT-S1, which is capable of predicting the rate independent, effective stress–strain–strength behaviour of uncemented soils over a wide range of confining pressures and densities. Freshly deposited sand specimens compressed from different initial formation densities approach a unique condition at high stress levels, referred to as the limiting compression curve (LCC), which is linear in a double logarithmic void ratio, e, mean effective stress space, p′. The model describes irrecoverable, plastic strains which develop throughout first loading using a simple four-parameter elasto-plastic model. The shear stiffness and strength properties of sands in the LCC regime can be normalized by the effective confining pressure and hence can be unified qualitatively, with the well-known behaviour of clays that are normally consolidated from a slurry condition along the virgin consolidation line (VCL). At lower confining pressures, the model characterizes the effects of formation density and fabric on the shear behaviour of sands through a number of key features: (a) void ratio is treated as a separate state variable in the incrementally linearized elasto-plastic formulation: (b) kinematic hardening describing the evolution of anisotropic stress–strain properties: (c) an aperture hardening function controls dilation as a function of ‘formation density’; and (d) the use of a single lemniscate-shaped yield surface with non-associated flow. These features enable the model to describe characteristic transitions from dilative to contractive shear response of sands as the confining pressure increases. This paper summarizes the procedures used to select input parameters for clays and sands, while a companion paper compares model predictions with measured data to illustrate the model capability for describing the shear behaviour of clays and sands. Copyright © 1999 John Wiley & Sons, Ltd.     

Surface wave studies for shear wave velocity and bedrock depth estimation over basalts

Shear wave velocity (V _S) estimation is of paramount importance in earthquake hazard assessment and other geotechnical/geo engineering studies. In our study, the shear wave velocity was estimated from ground roll using multichannel analysis of surface wave (MASW) technique making use of dispersive characteristics of Rayleigh type surface waves followed by imaging the shallow subsurface basaltic layers in an earthquake-prone region near Jabalpur, India. The reliability of MASW depends on the accurate determination of phase velocities for horizontally traveling fundamental mode Rayleigh waves. Inversion of data from surface waves resulted in a shear wave velocity (V _S) in the range of 200–1,200 m/s covering the top soil to weathering and up to bedrock corresponding to a depth of 10–30 m. The P-wave velocity (V _P) obtained from refraction seismic studies at these locations found to be comparable with V _S at an assumed specific Poisson’s ratio. A pair of selected set of V _S profiles over basalt which did not result in a hazardous situation in an earthquake of moderate magnitude are presented here as a case study; in other words, the shear wave velocity range of more than 200 m/s indicate that the area is highly unlikely prone to liquefaction during a moderate or strong earthquake. The estimated depth to basalt is found to be 10–12 m in both the cases which is also supported by refraction studies.     

An Experimental Characterization of Shear Wave Velocity (V<Subscript>s</Subscript>) in Clean and Hydrocarbon-Contaminated Sand

The characteristics of hydrocarbon-contaminated soils have been among major concerns of geotechnical engineers due to its significant frequency of event and also its influential consequences on our surroundings from various environmental and engineering viewpoints. Heretofore, the effects of diverse kinds of hydrocarbon contaminants on majority of geotechnical properties of fine- and coarse-grained soils such as grain size, hydraulic conductivity, plasticity, compressibility, internal friction, cohesion, and shear strength have been investigated. However, there has not been a concentrated research study examining shear wave velocity (\({\text{V}}_{\text{s}}\)) of hydrocarbon-contaminated soils as an important geotechnical property of soil due to this fact that, in small/very small strain levels, the maximum shear modulus of soils (\({\text{G}}_{ \hbox{max} }\)) can be determined using shear wave velocity (\({\text{G}}_{ \hbox{max} } =\uprho{\text{V}}_{\text{s}}^{2}\)). This paper aims to investigate effects of hydrocarbon contamination on shear wave velocity of sandy soils by comparing shear wave velocities in identically prepared clean and contaminated samples. To this aim, an Iranian light crude oil, a standard type of silica sand (Ottawa sand), and a bender element apparatus were used to minutely measure shear wave velocity of clean and crude oil contaminated sand samples. Moreover, dry and quasi-moist tamping methods were employed in order to provide comparable clean and contaminated specimens (containing 4, 6, 8, 10, and 12 wt% of crude oil), respectively. Firstly, a comprehensive bender element (BE) and resonant column tests were conducted on the identically prepared clean sand samples at various amounts of frequency (2–20 kHz) and under various confining pressure (50–500 kPa) to find the best methods of accurately determining shear wave travel time in BE tests. Thereafter, BE tests were conducted to examine shear wave velocity in contaminated specimens. Based on the results, it was found that there was a critical value for crude oil content with the maximum shear wave velocity so that shear wave velocity of 4 wt% contaminated sand (V_s-4 wt%) was about 1.2 times higher than clean one (V_s-clean), and contrastingly adding further crude oil up to 6 wt% made a significant reduction in value of shear wave velocity to some extent that V_s-6 wt% was slightly lower than V_s-clean (V_s-6 wt% = 0.95–0.97V_s-clean). Moreover, adding more contaminant (8–12 wt%) into sand had negligible influences on shear wave velocity. In this paper, the effects of crude oil contamination on sand microstructure were also evaluated using scanning electron microscopy.     

Soil void ratio correlation with shear wave velocities and SPT N values for Indo-Gangetic basin

In this study attempt has been made to understand in-situ void ratio in Indo-Gangetic basin (IGB) and to form empirical relations between void ratio and shear wave velocity (V_s), N values considering subsoil investigation data. Multichannel analysis of surface wave (MASW) test and standard penetration test was carried out along with soil property measured at 25 locations. The general soil profile varied from silty sand to clay of low compressibility, ground water level fluctuated between 1-27 m, depth of borehole varied from 20-40 m. Regression analysis was conducted on 202 data sets of void ratio and shear wave velocity, 293 data sets of void ratio and SPT- N value, which resulted in inverse correlations between void ratio and V_s, SPT N value. The datas were segregated into fine, coarse grained data based on engineering classification and relations were developed separately. Until now, no studies have related in-situ void ratio to V_s and SPT N. These correlations will be useful to predict void ratio for sites having measured values of V_s and N value. These void ratios can be further used to assess liquefaction susceptibility.     

On Garg's solution of Biot's equations for wave propagation in a one-dimensional fluid-saturated elastic porous solid

Garg's approximate analytical solutions of Biot's equations for wave propagation in a fluid-saturated elastic porous solid of infinite extent subjected to a velocity boundary condition of a Heaviside function at one end are examined for small and large drag. Garg's approximations were apparently introduced to facilitate exact inversion of Laplace transforms of certain quantities. The approximate solutions are compared with carefully evaluated numerical inverses of the Laplace transform solutions for different soils with widely varying properties. It is seen that for most soils (clays, silts and, sands) the error in Garg's approximate solutions in insignificant, and the solutions can be used as benchmarks for verifying numerical analysis procedures such as finite element codes.     

Unified 3D critical state bounding-surface plasticity model for soils incorporating continuous plastic loading under cyclic paths. Part I: Constitutive relations

The unified three-dimensional (3D) critical state bounding-surface plasticity model gUTS enables clays, silts and sands to be treated within a single framework. Furthermore, loose and dense states of a particular soil subjected to a wide range of confinements are viewed as a single material defined by the same set of constants. The model is able to handle both monotonic and complex cyclic paths including those involving a rotation of the principal stress directions. The model incorporates the following features: combined use of radial and deviatoric mapping rules and the use of an apparent normal consolidation line for sands; use of a non-associated flow rule where the ratio of the rates of volumetric plastic strain to deviatoric plastic strain is a function only of the ratio of deviatoric to mean effective stresses and the Lode angle; adoption of a bi-linear critical state line projected onto the plane of the void ratio versus logarithm of mean effective stress; inclusion of a sub-elliptic, or super-elliptic, segment in the plastic dilatancy surface for stress ratios less than critical; use of elliptic segments in the deviatoric planes; movement of the projection centre in the deviatoric mapping region and incorporation of a plastic stiffening effect for cyclic paths which repeatedly load in the same deviatoric direction.     

Evaluation of a constitutive model for clays and sands: Part I – sand behaviour

This paper evaluates the performance of a generalized effective stress soil model for predicting the rate independent behaviour of freshly deposited sands, while a companion paper describes model capabilities for clays and silts. Most material parameters can be obtained from standard laboratory data, including hydrostatic or one‐dimensional compression, drained and undrained triaxial shear testing. A compilation of data on compression behaviour allows for estimation of compression parameters when this type of data is not available. Extensive comparisons of model predictions with measured data from undrained triaxial shear tests shows that the model gives excellent predictions of the transition from dilative to contractive shear response as the confining pressure and/or the initial formation void ratio increases. A parametric study of drained response shows that the model describes realistically the variation of peak friction angle and dilation rate as a function of confining pressure and density when compared with an empirical correlation valid for many sands. The proposed formulation predicts a unique critical state locus for both drained and undrained triaxial testing which is non‐linear over the entire range of stresses and is in excellent agreement with recent experimental data. Overall, the model provides excellent predictions of the stress–strain–strength relationships over a wide range of confining pressures and formation densities. Copyright © 2002 John Wiley & Sons, Ltd.