Estimation of shearwave velocity in drifts using multichannel analysis of surface wave (MASW) technique — A case study from Jammu & Kashmir,India

Multichannel analysis of surface waves (MASW) is a non-destructive seismic prospecting method utilizing Rayleigh waves for imaging and characterizing shallow sub-surface structure. Multichannel analysis of surface waves (MASW) studies were conducted in drift areas of two bridge sites in the hilly terrain of J&K for imaging and characterizing shallow sub-surface structure. The purpose of the present study is to estimate the shear wave velocity (V_S) and subsurface structure in four drifts made in a hilly terrain for construction of two bridges. Rayleigh waves are having dispersive properties, travelling along or near the ground surface and are usually characterized by relatively low velocity, low frequency, and high amplitude. The study area comprises of Tertiary group of rocks which are underlain by Siwalik group. The main rock type in the study area is dolomite which has undergone various geological processes like weathering, jointing, fracturing and shearing. MASW data was collected inside four drifts in the mountainous terrain of J&K state which are located on either sides of Chenab river. The data was analyzed by relevant processing software using dispersion and inversion technique. Shear wave velocities were estimated up to 30 m depth. Average shear wave velocity (V_S ³⁰) up to top 30m was also computed. It is observed that, V_S in the range 400–800 m/s upto 10–15 m corresponding to weathered rock, followed by compact dolomite rock up to the depth of about 30 m with V_S in the range 1200–1600 m/s. Some low velocity zones are also identified from these sections which represent shear zones.     

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.     

Application of MASW in road failure investigation

Multichannel analysis of surface waves (MASW) survey was conducted to measure shear wave velocities in order to ascertain the likely causes of road failure along LASU-IBA expressway in Alimosho local government area, Lagos, Nigeria. MASW data were acquired along the express road. The acquired dataset was processed and transformed into two-dimensional structure reflective of depth and surface wave velocity distribution within a depth of investigation using SurfSeis software. The MASW shear wave velocity data were compared to geophysical data that was acquired along the same profile. The comparison was also done with geotechnical data that had been acquired prior to the study some meters away from the study area. The correlation between N values to measured shear wave velocity using MASW was generated. The comparison illustrates the accuracy and consistency of MASW-derived shear wave velocity profiles. We concluded that (1) the low-velocity region that varies between 100 and 250 m/s at surface down to 4 m beneath the surface is characterized by loose/peat materials and may have been responsible for the road failure within the study area; this region depicts a very loose compaction area. (2) The MASW technique is a time–cost-effective tool for obtaining reliable shear wave velocity profiles, and (3) the MASW is particularly attractive in areas that cannot be readily assessed by other geophysical and geotechnical tools.     

Surface soil assessment in the Ubhur area,north of Jeddah,western Saudi Arabia,using a multichannel analysis of surface waves method

Ten boreholes drilled in Ubhur area up to the depth of bedrock indicted the shallow depth of bedrock where the average depth ranges between 10 and 15 m. The standard penetration test N-values of these boreholes were measured and averaged. Based on N-values to the depth of bedrock, Ubhur area can be classified as site class C and D. Multichannel analysis of surface waves technique has been applied along seventy six profiles using 24-channel geophone array and 4.5Hz vertical geophones with 1m geophone spacing and sledgehammer and/or weight drop as seismic energy sources. Values of shear wave velocity to 30 m are calculated and then averaged (Vs30) where it ranges between 310.08 m/s and 1139.8 m/s. Therefore, Ubhur area can be classified into site class B, C and D based on site classification of the national earthquake hazards reduction program (NEHRP) recommendations. Accordingly, the greatest part of the study area falls in site class C while class B and D covered limited areas in the western and the eastern parts respectively. Depending on the shallow depth of bedrock in the study area, the Vs30 parameter is not applicable in the study area so the average values of Vs for the soil thickness, excluding the bedrock, have been calculated and mapped for site class C and D only. So Vs30 approach is not applicable for areas with shallow depth of bedrock which gives higher classification.     

Multichannel analysis of surface waves (MASW) for seismic site characterization using 2D genetic algorithm at Bahrah area,Wadi Fatima,Saudi Arabia

The present study deals with the seismic site classification of Bahrah area, Wadi Fatima, to characterize the local site conditions. The dynamic behavior of sediments was studied by the application of surface wave inversion. The multichannel analysis of surface waves (MASW) shallow geophysical technique was utilized for site classification. MASW survey was carried out at 66 sites along with 13 seismic refraction profiles at suitable localities. MASW and seismic refraction profiles were processed and compared with the available borehole data. The integration of MASW techniques with seismic refraction and borehole data progressively enhanced the subsurface visualization and reliability of the shear wave velocity estimation in the subsurface in the study area. The subsurface shear-wave velocity model was achieved by the solution of an inverse problem-based dispersion of surface waves and propagation in a vertically heterogeneous medium. The 2D genetic algorithm was employed for the inversion of dispersion curves to obtain velocity and thickness of subsurface layers. The depth to engineering bedrock and velocity of shear waves in the first 30 m was deciphered and mapped. The depth of bedrock in study area varies from 4 to 30 m, and V _{S 30} ranges from 320 to 800 m/s. The most of study area falls in B and C class categories in addition to few sites of D class according to the NEHRP guidelines.     

Shear Wave Velocities in the Estimation of Earthquake Hazard Over Alluvium in Seismically Active Region

It is well known that the potential hazard during an earthquake is mainly in alluvium or alluvium filled basins; shear wave velocity plays a significant role in estimating the possible hazard during an earth quake in such an area. This paper presents shear wave velocity profile from Jabalpur, Central India mainly dominated by alluvial soil that was moderately affected by an earthquake of magnitude 6.5 in May, 1997. The acquired shear wave velocity by Multichannel Analysis of Surface Wave (MASW) in as many as 36 sites over alluvial soil ranges from 200 m/sec to 400 m/sec and in a few sites marginally less than 200 m/sec corresponding to a sub-surface depth of 30–35 m. Further, the computed N values vary as low as near zero to less than 25. The study is substantiated by the estimation of P-wave velocity by refraction seismic method at the same locations of MASW which ranges from 350 m/sec to 2200 m/sec. The results suggest that the damage during an earth quake appears to be highly unlikely in view of the marginally high Vs up to depth of 30 m. This study on seismic hazard is substantiated by the estimation of frequency of the ground as well as amplification which is found to be a maximum of about 2.5 in the frequency band of 2–6 Hz in west and north western portion of the study area.     

Noise Reduction Technique Applied to the Multichannel Analysis of Surface Waves

This paper shows the presence of noises and technique to reduce these noises during the surface wave analysis. The frequency-dependent properties of Rayleigh-type surface waves can be used for imaging and characterizing the shallow subsurface. Interference by coherent source-generated noise inhibits the reliability of shear-wave velocities determined through inversion of the phase velocities of Rayleigh waves. Among these interferences by non-planar, non-fundamental mode Rayleigh waves (noise) are body waves, scattered and non-source-generated surface waves, and higher-mode surface waves. For the reduction of noise, the filtering technique is implemented in this paper for the multichannel analysis of surface wave method (MASW). With the de-noising technique during the MASW method, more robust and reliable outcome is achieved. The significance of this paper is to obtain pre-awareness about noises during surface wave analysis and take better outcomes with de-noising performance in near surface soil investigations.     

Liquefaction Hazard Assessment of Earth Quake Prone Area: a Study Based on Shear Wave Velocity by Multichannel Analysis of Surface Waves (MASW)

The shear wave velocity (V_S) profile based on the dispersive characteristics of fundamental mode of Rayleigh type surface waves indicate underground stiffness change with depth as well as near surface stiffness. The most important utility of shear wave velocity (V_S) is to estimate the liquefaction hazard potential of an area particularly in seismically active region. Rayleigh type surface waves were utilized to estimate the velocity (V_S) of shallow subsurface covering a depth range of 30–50 m employing multichannel analysis of surface waves. The liquefaction hazard map predicts an approximate percentage of an area that will have surface manifestation of liquefaction during an earth quake. The surface wave data acquired in an earth quake prone region of Jabalpur (Seismic zone III), India, yields a velocity (V_S) range of 200–750 m/s corresponding to the subsurface depth of 30–35 m. The results were analyzed for possible liquefaction hazard in the study area and presented here besides the N values.     

土层剪切波速不确定性对场地刚性判断的影响

岩土工程中的不确定性是工程风险的重要来源。在我国7个主要地区的20个工程场地上开展了国内外首次大范围现场剪切波速不确定性专项试验,共由47家单位用11种工程常用仪器采用单孔法完成,获取了600组土层剪切波速Vs实测数据。依据实测结果拟合了场地时间平均剪切波速Vs, z的变异系数COV与计算深度Z的对应公式,揭示了土层剪切波速不确定性对场地刚性判断的影响,进而分析了我国场地分类的可能潜在误判区域。结论如下：在目前我国现场剪切波速测试水平下,Vs, z的变异系数与深度负相关——在近地表处达到极大值15%,并随着深度增加明显减小;地表至10 m深度区间内,Vs, z的变异系数迅速减小;10～20 m深度区间内,Vs, z的变异系数减小趋势变缓;超过20 m深度后,Vs, z的变异系数小于5%且基本不变;我国III、IV类场地分类采用的等效剪切波速Vse和欧美场地分类指标Vs, 30的变异系数基本相等;依据我国建筑抗震设计规范GB50011－2010,对I类场地以及覆盖层厚度小于5 m或等效剪切波速接近500 m/s的II类场地分类时需考虑剪切波速不确定性的潜在影响,其他类型的II类场地、III类和IV类场地在分类时基本可忽略其影响。     

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.     

Research of reservoir watershed fine zoning and flood forecasting method

In the present study, a set of correlation relations between standard penetrations tests (SPT-N values) and shear wave velocity (Vs) for different categories of soils is developed for Dholera region, Gujarat state, Western India. Shear wave velocities were measured using multichannel analysis of surface wave (MASW at 42 sites) and by PS logging (at 16 sites). SPT-N values were taken in total 87 geotechnical boreholes at depth interval of 3 m. Seismic site classification is done as per National Earthquake Hazards Reduction Program of the USA. Total 336 pairs of SPT-N values and Vs data at different depths are used to develop the regression correlation between uncorrected SPT-N value and Vs and compared with the results of other workers from India and worldwide for all soils type. It is found that regression correlations developed for clay and all soils are almost similar to each other, for sand, the coefficient value is less and for silt, it is higher. The new regression equation gives good prediction performance. The present correlation can be used for the seismic hazard study for the study area and also for the other areas having similar soil strata using a process of validation.     

On the use of rayleigh wave group velocities for the analysis of continental margins

Rayleigh wave group velocities provide a low-cost means for a quick assessment of averaged local properties of the Earth's crust in continental margin regions of the Atlantic type. Sufficiently accurate measurements (with a standard error of 0.3 km/s or less) of group velocities in continental shelf areas at periods between 5 and 30 seconds provide important information about structural parameters. They may resolve the Moho depth to within 4 or 5 km, depending on crustal thickness and give useful estimates of the average velocities in the upper part of the crust. The group velocities of Rayleigh waves in this period range are influenced most by the shear velocity at all depths and the compressional velocity and the density near the surface. For continental rise regions, the dominating influence of the water layer limits the effectiveness of the method.     

Rayleigh面波各阶模式频散曲线对横波速度和层厚的敏感性探讨

Rayleigh面波勘探的目的在于有效利用频散曲线反演地层厚度及横波速度,而不同模式的频散曲线对横波速度和层厚的敏感性不同。通过求取介质参数变化10%后与参数不变化时的二组频散曲线的差值,得到各阶模式的频率～相速度差曲线,分析了Rayleigh面波各模式频散曲线对横波速度、层厚的敏感性。试验结果表明,基阶模式对于浅层的横波速度和层厚比较敏感,敏感区域主要集中在较窄的频带范围内。而高阶模式对于相对较深层的横波速度和层厚比较敏感,且频率范围分布较大,敏感性强的频段分布比较分散。研究结果可以为Rayleigh面波多模式联合反演提供理论依据。     

Advanced performance in geotechnical engineering using tomography analysis

This paper presents the technique to demonstrate reliable two-dimensional (2-D) tomography of near-surface soil through multichannel analysis of surface wave (MASW) method. The MASW method is used as seismic method for determining the shear wave velocity profile of near-surface soil with better performance. Normally, this method represents one-dimensional (1-D) soil layer profile to delineate anomalous subsurface materials and detect soil characteristics in geotechnical investigations. Conventionally, the MASW method is able to represent depth corresponding shear wave velocity through 1-D representation. In this paper, this method is improved through representation of depth and distance versus shear wave velocity profile using 2-D tomography analysis in geotechnical investigations. The outcomes of the newly developed method are more reliable and informative in comparison to the outcome of the conventional MASW method. The significance of this research is incorporating advanced tomography technique with MASW method to obtain 2-D tomography of geotechnical characteristics with consistency.     

基于自动安置排列的多道面波分析方法

横波速度由于更能直接反映岩石的物理性质,从而在许多工程、环境和地下水勘查中有重要作用.多道面波分析方法(MASW)是一种新技术,通过求得的一维近地表横波速度剖面构建一个二维横波速度场,从而准确探测近地表介质情况.针对该方法数据的采集,提出了一种检波器排列自动安置的技术(自动安置排列),即使用牵引车辆的液压驱动装置在几秒内自动安置几十个检波器,并将其与传统安置排列的方法做了比较,结果表明,自动安置排列技术可用于多道面波分析方法,而且可以大大减少多道面波分析方法的采集时间和费用.     

Crustal and upper mantle S-wave velocity structures across the Taiwan Strait from ambient seismic noise and teleseismic Rayleigh wave analyses

Although orogeny tapers off in western Taiwan large and small earthquakes do occur in the Taiwan Strait, a region largely untouched in previous studies owing mostly to logistical reasons. But the overall crustal structure of this region is of particular interest as it may provide a hint of the proto-Taiwan before the orogeny.By combining time domain empirical Green’s function (TDEGF) from ambient seismic noise using station-pairs and traditional surface wave two-station method (TS) we are able to construct Rayleigh wave phase velocity dispersion curves between 5 and 120 s. Using Broadband Array in Taiwan for Seismology (BATS) stations in Taiwan and in and across the Strait we are able to derive average 1-D Vs structures in different parts of this region. The results show significant shear velocity differences in the upper 15 km crust as expected. In general, the highest Vs in the upper crust observed in the coastal area of Mainland China and the lowest Vs appears along the southwest offshore of the Taiwan Island; they differ by about 0.6–1.1 km/s. For different parts of the Strait, the upper crust Vs structures are lower in the middle by about 0.1–0.2 km/s relative to those in the northern and southern parts. The upper mantle Vs structure (Moho – 150 km) beneath the Taiwan Strait is about 0.1–0.3 km/s lower than the AK135 model. The overall crustal thickness is approximately 30 km, much thinner and less variable than under the Taiwan Island. The inversion of seismic velocity structures using shorter period band dispersion data in the sea areas with water depth deeper than 1000 m should take water layer into consideration except for the continental shelves.     

NEHRP soil classification and estimation of 1-D site effect of Dehradun fan deposits using shear wave velocity

Shear wave velocity of the near surface soil at nearly 50 sites in the sub Himalayan mountain exit covering Doon fan deposits, was determined using Multi-channel Analysis of Surface Waves (MASW), a seismic reflection technique. Based on the average shear wave velocity of the upper 30 m soil column, sites in the Dehradun fan are predominantly classified as class ‘D’ (180–360 m/s). Similarly, sites located on the northwestern, eastern and southeastern sides of the fan deposit have shear wave velocities (in the upper 30 m soil) greater than 360 m/s, thereby classifying them as class ‘C’ (360–760 m/s) in accordance with NEHRP provisions. Some of the sites towards the southwestern side of the fan deposits had average shear wave velocities less than 180 m/s and could be classified as soil class ‘E’. One dimensional site effects, including amplification and dynamic period were calculated for the majority of the sites. However, some of the representative suite of sites across the north–south profile of Dehradun fan has been discussed here. Although the attenuation is greater on the southwestern side of the Dehradun fan deposits (i.e. thicker, low velocity sediments) and the sites had been classified as class ‘D’ and ‘E’ but the site amplification tends to be greater in the northern and northwestern part of the city due to large impedance contrast with in the near surface soils.     

High-Frequency Rayleigh-Wave Method

High-frequency (≥2 Hz) Rayleigh-wave data acquired with a multichannei recording sys-tem have been utilized to determine shear (S)-wave velocities in near-surface geophysics since the early 1980s. This overview article discusses the main research results of high-frequency surface-wave tech-niques achieved by research groups at the Kansas Geological Survey and China University of Geosciences in the last 15 years. The multichannel analysis of surface wave (MASW) method is a nou-iuvasive acoustic approach to estimate near-surface S-wave velocity. The differences between MASW results and direct borehole measurements are approximately 15% or less and random. Studies show that simultaneous inversion with higher modes and the fundamental mode can increase model resolution and an investigation depth. The other important seismic property, quality factor (Q), can also be estimated with the MASW method by inverting attenuation coefficients of Rayleigh waves. An inverted model (S-wave velocity or Q) obtained using a damped least-squares method can be assessed by an optimal damping vector in a vicinity of the inverted model determined by an objective function, which is the trace of a weighted sum of model-resolution and model-covariance matrices. Current developments include modeling high-frequency Rayleigh-waves in near-surface media, which builds a foundation for shallow seismic or Rayleigh-wave inversion in the time-offset domain; imaging dispersive energy with high resolution in the frequency-velocity domain and possibly with data in an arbitrary acquisition geometry, which opens a door for 3D surface-wave techniques; and successfully separating surface-wave modes, which provides a valuable tool to perform S-wave velocity profiling with high-horizontal resolution.     

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.     

Directional dependence of P- and S-wave propagation and polarization in foliated rocks from the Kola superdeep well: Evidence from laboratory measurements and calculations based on TOF neutron diffraction

We have measured P- and S-wave velocities on two amphibolite and two gneiss samples from the Kola superdeep borehole as a function of pressure (up to 600 MPa) and temperature (up to 600 °C). The velocity measurements include compressional (Vp) and shear wave velocities (Vs1, Vs2) propagating in three orthogonal directions which were in general not parallel to inherent rock symmetry axes or planes. The measurements are accompanied by 3D-velocities calculations based on lattice preferred orientation (LPO) obtained by TOF (Time Of Flight) neutron diffraction analysis which allows the investigation of bulk volumes up to several cubic centimetres due to the high penetration depth of neutrons. The LPO-based numerical velocity calculations give important information on the different contribution of the various rock-forming minerals to bulk elastic anisotropy and on the relations of seismic anisotropy, shear wave splitting, and shear wave polarization to the structural reference frame (foliation and lineation). Comparison with measured velocities obtained for the three propagation directions that were not in accordance with the structural frame of the rocks (foliation and lineation) demonstrate that for shear waves propagating through anisotropic rocks the vibration directions are as important as the propagation directions. The study demonstrates that proper measurement of shear wave splitting by means of two orthogonal polarized sending and receiving shear wave transducers is only possible when their propagation and polarization directions are parallel and normal to foliation and lineation, respectively.