Static Liquefaction Analysis Considering Principal Stress Directions and Anisotropy

Static liquefaction failure of a sloping ground occurs when the shear stress applied by a monotonic triggering load exceeds the undrained yield (peak) shear strength of the saturated liquefiable cohesionless soil. Current practices for determining the in-situ undrained yield strength for ground subject to static shear stress either rely on a suite of costly laboratory tests on undisturbed field samples or empirical correlations based on in-situ penetration tests which do not account for the effects of anisotropic consolidation, intermediate principal stress, and mode of shear on the degree of strain-softening and brittleness of cohesionless soils. This study investigates the effects of variations in the direction and relative magnitudes of principal stresses associated with different modes of shear and ground slopes on static liquefaction failure of cohesionless soils. Empirical relationships are developed between soil brittleness index and maximum excess pore water pressure ratio to characterize soil shearing behavior observed in a database of 271 undrained laboratory shear tests collected from the past literature. The application of these relationships for estimating the static liquefaction triggering strength of cohesionless soils under sloping grounds is described for plane-strain boundary conditions and the results are compared with those back-calculated for several cases of static liquefaction flow failures. The proposed procedure incorporates variations in mode of shear and initial stress anisotropy in an empirical formulation based on in-situ penetration tests.     

Liquefaction analysis and soil improvement in Beydag dam

Beydag dam is under construction on Kucukmenderes River for irrigation purposes. Due to the scarcity of core material and liquefaction of alluvium at the dam site, the original design was changed to Roller Compacted Concete (RCC) from rockfill dam with claycore. Although the new design was safer, it nearly doubled the cost of the dam, so the owner, State Hydraulic Works of Turkey, (DSI) set out to find more economical but equally safe alternative. Since jet-grouting is a cheap ground improvement tool in Turkey, such an alternative was developed for the ground improvement against liquefaction together with concrete face rockfill dam sitting on top of improved ground. This paper presents a detailed discussion of how the new alternative was developed and evaluated: it discusses the determination of jet grouting pattern, the placement of jet grouted blocks, and the assesment of liquefaction. On one hand the soil cement strength of jetgrout columns, internal friction angle of alluvium and rockfill were important in determining the dimensions of the blocks, on the other hand the location of the blocks were highly affected by the areas where liquefaction occurred. One of the most important parameter that has a considerable influence in delineating the boundary betweeen liquefaction and non-liquefaction was the value of stress reduction coefficient (r _d), being primarily sensitive to the weight of overburden, which is calculated by the height from the face of dam to the depth where the calculation was made. This approach is justified by two-dimensional ground response analysis. Most importantly, this paper shows that there exists an alternative solution for building dams on the liquefaction prone sites without removing alluvium by using a well known jet grouting technique for improving ground at only selected places.     

Probabilistic methodology for the analysis of paleoliquefaction features

Paleoliquefaction features can be used to estimate lower bounds on the magnitude and ground motion associated with the earthquake that caused the liquefaction feature. The engineering back-analysis of paleoliquefaction features is usually conducted using state of the practice liquefaction-triggering analysis methodologies. Recent studies have shown that these methodologies are associated with variable probabilities of liquefaction depending on the soil parameters. This would imply that estimates of magnitude and ground motion intensity obtained from these methodologies would not be consistent for all soil sites. Moreover, these estimates could be unconservative. In this paper, the use of a probabilistic methodology for the back-analysis of paleoliquefaction features is proposed. The proposed methodology permits the incorporation of model and parameter uncertainty into the analysis and results in more robust estimates of past magnitude and a measure of the uncertainty associated with these predictions. Previously published paleoliquefaction data are used to demonstrate the applicability of the proposed method. Magnitude estimates obtained with the proposed method do not differ significantly from those obtained using deterministic methodologies, but the proposed methodology permits a quantification of the uncertainty associated with magnitude estimates.     

以标贯试验为依据的砂土液化确定性及概率判别法

核电厂址非基岩场地的地基液化问题是核电厂选址的关键问题,亟需建立核电厂址地基液化判别方法。回顾了以标贯试验和地表峰值加速度为依据的砂土液化判别方法的演化历史,依据Idriss-Boulanger确定液化临界曲线的基本方法,提出了确定液化临界曲线的基本原则,分别依据美国液化数据库、中国抗震规范液化判别式所用的液化数据及综合两者的液化数据资料,给出了相应的液化临界曲线,验证了液化临界曲线的位置对不同的细粒含量、有效上覆压力、现场试验方法的液化数据的合理性,分析了测量或估计土层循环应力比和修正标贯击数各种因素的不确定性对液化临界曲线的敏感性,结果表明：所提的液化临界曲线不易受各种因素的影响。利用Monte Carlo模拟、加权最大似然法和加权经验概率法,给出了液化临界曲线的名义抗液化安全系数与液化概率的经验关系式及概率等值线,并对核电厂Ⅰ类、Ⅱ类和Ⅲ类抗震物项地基,给出了相应的液化临界曲线。     

Liquefaction evaluation discrepancies in tropical lagoonal soils

Field penetration tests and shear wave velocity measurements are both established and accepted methods for evaluating liquefaction potential in soils. The results produced by the two methods are generally well correlated. However, recent studies have shown that when investigating tropical lagoonal deposits, the same accepted methods for evaluating liquefaction potential often produce significant discrepancies in results. This discrepancy is most apparent in saturated lagoonal deposits of calcareous gravelly sand (or sandy gravel), which tend to exhibit low penetration resistance values but relatively high shear wave velocities. These disparate test results can suggest different soil classifications under current building codes. Ambiguity in the code may allow for a potentially unconservative classification, which may in turn allow for the use and construction of less costly, lighter weight foundation systems than warranted. Equally as important, the potential for unconservative design as related to liquefaction appears to be high when shear wave velocity measurements are used as a basis for evaluation in these types of lagoonal deposits. Because of this, it is strongly recommended that caution should be excercised when determining seismic design parameters in these types of geologic environments. A hypothesis to explain the discrepancies in the results of evaluation methods and a suggested design protocol is proposed.     

Assessment and Prediction of Liquefaction Potential Using Different Artificial Neural Network Models: A Case Study

Soil liquefaction as a transformation of granular material from solid to liquid state is a type of ground failure commonly associated with moderate to large earthquakes and refers to the loss of strength in saturated, cohesionless soils due to the build-up of pore water pressures and reduction of the effective stress during dynamic loading. In this paper, assessment and prediction of liquefaction potential of soils subjected to earthquake using two different artificial neural network models based on mechanical and geotechnical related parameters (model A) and earthquake related parameters (model B) have been proposed. In model A the depth, unit weight, SPT-N value, shear wave velocity, soil type and fine contents and in model B the depth, stress reduction factor, cyclic stress ratio, cyclic resistance ratio, pore pressure, total and effective vertical stress were considered as network inputs. Among the numerous tested models, the 6-4-4-2-1 structure correspond to model A and 7-5-4-6-1 for model B due to minimum network root mean square errors were selected as optimized network architecture models in this study. The performance of the network models were controlled approved and evaluated using several statistical criteria, regression analysis as well as detailed comparison with known accepted procedures. The results represented that the model A satisfied almost all the employed criteria and showed better performance than model B. The sensitivity analysis in this study showed that depth, shear wave velocity and SPT-N value for model A and cyclic resistance ratio, cyclic stress ratio and effective vertical stress for model B are the three most effective parameters on liquefaction potential analysis. Moreover, the calculated absolute error for model A represented better performance than model B. The reasonable agreement of network output in comparison with the results from previously accepted methods indicate satisfactory network performance for prediction of liquefaction potential analysis.     

Liquefaction potential of level ground: Deterministic and probabilistic assessment

Two important aspects of the liquefaction of level ground are discussed; the liquefaction potential of a site due to specified earthquake motions and the total probability of liquefaction at a site due to the entire seismic environment surrounding the site. Dynamic effective stress analysis is used to explore two important case histories of liquefaction due to specified ground motions. The probability of liquefaction is assessed by a method which estimates the joint probability of occurence of earthquake magnitude and acceleration amplitude sufficient to cause liquefaction.     

Nonlinear soil response to ground motions during different earthquakes in Nepal,to arrive at surface response spectra

Catastrophic damages reported during an earthquake include building damages, excessive ground shaking, uneven settlements and liquefaction. While most of the seismic hazard studies map the probable level of ground shaking at the bedrock level, their use in assessing the above damages is very limited until the response of the local soil is also taken into account. Determination of the local soil response needs regionally recorded ground motions, dynamic soil properties, in situ geotechnical details, etc., which most of the time are not readily available for the region under study. In the present work, the response of local soil for Nepal has been studied indirectly taking into account the surface level of ground shaking during various past as well as recent EQs observed at various locations. Based on the present analysis, a low value of amplification factor for high peak horizontal acceleration and vice versa is observed in central, western as well as southern parts of Nepal. These observations suggest nonlinear soil behavior and are in accordance with the available literature. Further, the ground motion records during 2015 Nepal EQ show maximum soil response at 0.3 s which is exactly matching with the site class C obtained from in situ data for the above locations. Based on the above observations, various correlations between the high peak horizontal acceleration and the surface spectral acceleration are proposed to obtained site specific surface response spectrum for Nepal.     

Laboratory Investigation into the Effects of Silty Fines on Liquefaction Susceptibility of Chlef (Algeria) Sandy Soils

In a number of recent case studies, the liquefaction of silty sands has been reported. To investigate the undrained shear and deformation behaviour of Chlef sand–silt mixtures, a series of monotonic and stress-controlled cyclic triaxial tests were conducted on sand encountered at the site. The aim of this laboratory investigation is to study the influence of silt contents, expressed by means of the equivalent void ratio on undrained residual shear strength of loose, medium dense and dense sand–silt mixtures under monotonic loading and liquefaction potential under cyclic loading. After an earthquake event, the prediction of the post-liquefaction strength is becoming a challenging task in order to ensure the stability of different types of earth structures. Thus, the choice of the appropriate undrained residual shear strength of silty sandy soils that are prone to liquefaction to be used in engineering practice design should be established. To achieve this, a series of undrained triaxial tests were conducted on reconstituted saturated silty sand samples with different fines contents ranging from 0 to 40 %. In all tests, the confining pressure was held constant at 100 kPa. From the experimental results obtained, it is clear that the global void ratio cannot be used as a state parameter and may not characterize the actual behaviour of the soil as well. The equivalent void ratio expressing the fine particles participation in soil strength is then introduced. A linear relationship between the undrained shear residual shear strength and the equivalent void ratio has been obtained for the studied range of the fines contents. Cyclic test results confirm that the increase in the equivalent void ratio and the fines content accelerates the liquefaction phenomenon for the studied stress ratio and the liquefaction resistance decreases with the increase in either the equivalent void ratio or the loading amplitude level. These cyclic tests results confirm the obtained monotonic tests results.     

Liquefaction analysis and damage evaluation of embankment-type structures

The increasing importance of performance-based earthquake engineering analysis points out the necessity to assess quantitatively the risk of liquefaction of embankment-type structures. In this extreme scenario of soil liquefaction, devastating consequences are observed, e.g., excessive settlements, lateral spreading and slope instability. The present work discusses the global dynamic response and interaction of an earth structure-foundation system, so as to determine quantitatively the collapse mechanism due to foundation’s soil liquefaction. A levee-foundation system is simulated, and the influence of characteristics of input ground motion, as well as of the position of liquefied layer on the liquefaction-induced failure, is evaluated. For the current levee model, its induced damage level (i.e., induced crest settlements) is strongly related to both liquefaction apparition and dissipation of excess pore water pressure on the foundation. The respective role of input ground motion characteristics is a key component for soil liquefaction apparition, as long duration of mainshock can lead to important nonlinearity and extended soil liquefaction. A circular collapse surface is generated inside the liquefied region and extends toward the crest in both sides of the levee. Even so, when the liquefied layer is situated in depth, no significant effect on the levee response is found. This research work provides a reference case study for seismic assessment of embankment-type structures subjected to earthquake and proposes a high-performance computational framework accessible to engineers.     

1998 Adana-Ceyhan (Turkey) Earthquake and a Preliminary Microzonation Based on Liquefaction Potential for Ceyhan Town

A moderate earthquake (M_s = 6.2) occurred in the Cukurova region in the southern part of Turkey, on 27 June 1998. It resulted in loss of 145 lives and significant damage particularly in the settlements close to the epicenter at the south of Ceyhan town. Widespread liquefaction and associated sand boils, ground fissures and ground deformations due to lateral spreading occurred during this earthquake. In this study, main characteristics of the earthquake are presented and liquefaction throughout the site was assessed. An attempt was also made to establish preliminary microzonation maps for Ceyhan using the data from liquefaction susceptibility analyses. The results of the analyses indicated that the data from the liquefied sites were within the empirical bounds suggested by the field-performance evaluation method. Fortunately, most of the riversides were used for agricultural purposes alone, damage to structures from liquefaction and associated ground failures were rather limited. Preliminary assessments indicated that at depths of about 5 m the liquefaction potential of thin sand layers tends to diminish.     

Ground Investigations for Changi East Reclamation Projects

Geotechnical ground investigation is important in all construction works particularly, land reclamation projects. The Changi East reclamation project is a mega project that involves creation of land space as large as 3,000 hectares in multiple phases between 1991 and 2005. Ground investigation plays a critical role in the Changi project which involved substantial amount of ground improvement works after reclamation by sand filling. One important aspect of the investigation was to evaluate the geotechnical characteristics of the underlying compressible soils beneath the original seabed for the selection of suitable methods of ground improvement. The other important aspect of the investigation was to characterize and verify the compactness of the reclaimed sand fill. Many types of in situ tests were extensively used in the project. The field vane shear test provides correlations between the undrained shear strength and depth for the upper and the lower marine clay strata in the Northern part and the Southern Part of the reclaimed site. Combined with laboratory undrained shear strength tests, the field vane strength profiles allow the derivation of the undrained shear strength over effective stress ratios for the upper and the lower marine clays. In addition, results from three other in situ tests, the piezocone cone penetration test, the flat dilatometer test, and the self boring pressuremeter test provided useful verification of these correlations. A field-performance based method of assessing degree of consolidation in the underlying clay was developed by combining field monitoring, laboratory testing and conventional as well as specialized in situ testing. An effective use of in situ testing methods for assessment, the outcome of ground improvement and for fill quality control in the densification of granular soils is illustrated with field observation data collected at the project site.     

黄河口不同强度粉土液化特性的试验研究

研制真空压缩装置以制备不同强度的粉土土样,并进行动三轴试验,测试不同强度粉土的液化特性。建立了土的强度与动剪应力比的关系,确立孔压增长模型及其参数,探讨土的强度和动应力对液化的影响。结果表明,（1）真空压缩装置制备土样过程中,土体超孔压在24 h内基本消散完毕,28 d贯入阻力达300～400 N,不排水抗剪强度达8 kPa,达到黄河口软弱土的强度;（2）土的强度与动剪应力比基本呈线性关系,且土的强度越高,孔压增长曲线越呈现上凸趋势,破坏时的孔压比也越大;（3）指数模型能够较好模拟黄河口粉土孔压增长情况,其中参数a和b分别位于0.77～5.63和0.17～4.65之间;对于孔压比上限,参数a, b分别为0.92和4.65;孔压比下限,参数a, b分别为1.25和0.89。     

基于土层常规参数的剪切波速液化概率计算公式

分析国际上现有液化场地剪切波速473组数据进行参数特征后,以此为基础,采用分区方法,利用成熟的Logistic模型,提出了以地表峰值加速度、剪切波速、地下水位、可液化层埋深等4个参数表达的土体液化概率计算公式和不同概率下液化临界值计算公式,研究了不同概率水平下公式的适用性,并与现有主要方法进行比较。研究表明:地震动强度为液化判别首要影响参数,液化层与非液化层剪切波速区分度不显著,采用以往一个公式构造液化判别公式的方式难以要达到基本要求;现有Andrus公式和石兆吉公式会将很多明显为非液化的场地误判为液化场地,违背了现有认识,达到了不可接受的程度,需要改进。文中公式取50%概率时液化点和非液化点回判成功率符合对等原则,不同烈度下成功率均在70%左右;总体看,公式表现良好,可为中国工程建设提供一个合理、可操作性强的剪切波速液化概率判别方法。     

Assessment of Liquefaction Potential of Guwahati City: A Case Study

The liquefaction potential of saturated cohesionless deposits in Guwahati city, Assam, was evaluated. The critical cyclic stress ratio required to cause liquefaction and the cyclic stress ratio induced by an earthquake were obtained using the simplified empirical method developed by Seed and Idriss (J soil Mech Found Eng ASCE 97(SM9):1249–1273, 1971, Ground motions and soil liquefaction during earthquakes. Earthquake Engineering Research Institute, Berkeley, CA, 1982) and Seed et al. (J Geotech Eng ASCE 109(3):458–483, 1983, J Geotech Eng ASCE 111(12):1425–1445, 1985) and the Idriss and Boulanger (2004) method. Critical cyclic stress ratio was based on the empirical relationship between standard penetration resistance and cyclic stress ratio. The liquefaction potential was evaluated by determining factor of safety against liquefaction with depth for areas in the city. A soil database from 200 boreholes covering an area of 262 km² was used for the purpose. A design peak ground acceleration of 0.36 g was used since Guwahati falls in zone V according to the seismic zoning map of India. The results show that 48 sites in Guwahati are vulnerable to liquefaction according to the Seed and Idriss method and 49 sites are vulnerable to liquefaction according to the Idriss and Boulanger method. Results are presented as maps showing zones of levels of risk of liquefaction.     

Arias Intensity Assessment of Liquefaction Test Sites on the East Side of San Francisco Bay Affected by the Loma Prieta,California, Earthquake of 17 October 1989

Uncompacted artificial-fill deposits on the east side of San Francisco Bay suffered severe levels of soil liquefaction during the Loma Prieta earthquake of 17 October 1989. Damaged areas included maritime-port facilities, office buildings, and shoreline transportation arteries, ranging from 65 to 85 km from the north end of the Loma Prieta rupture zone. Typical of all these sites, which represent occurrences of liquefaction-induced damage farthest from the rupture zone, are low cone penetration test and Standard Penetration Test resistances in zones of cohesionless silty and sandy hydraulic fill, and underlying soft cohesive Holocene and Pleistocene sediment that strongly amplified ground motions. Postearthquake investigations at five study sites using standard penetration tests and cone penetration tests provide a basis for evaluation of the Arias intensity-based methodology for assessment of liquefaction susceptibility.     

A study on the identification of liquefaction-induced failures on ground surface based on the data from the 1999 Kocaeli and Chi-Chi earthquakes

One of the major causes of earthquake damage is liquefaction. However, it doesn't result in severe harm unless it leads to ground surface damage or ground failure. Therefore, prediction of potential for ground surface damage due to liquefaction is one of the important issues in microzonation studies for liquefaction-induced damage in areas with high seismicity. In 1985, based on a database compiled from Chinese and Japanese earthquakes, Ishihara considered the influence of the non-liquefied cap soil on the occurrence or non-occurrence of ground failure (mainly sand boiling), and proposed an empirical approach to predict the potential for ground surface damage at sites susceptible to liquefaction. However, some investigators indicated that this approach is not generally valid for sites susceptible to lateral spread or ground oscillation. In this study, a contribution to improve the approach by Ishihara is made. For the purpose, an index called liquefaction severity index (LSI) and data from two devastating earthquakes, which occurred in Turkey and Taiwan in 1999, were employed. The data from liquefied and non-liquefied sites were grouped and then analysed. Based on the observations reported by reconnaissance teams who visited both earthquake sites and the results of the liquefaction potential analyses using the filed-performance data, a chart to assess the potential for ground surface disruption at liquefaction-prone areas was produced. The analyses suggest that the procedure proposed by Ishihara is quite effective particularly for the occurrence of sand boils, while the bounds suggested in this method generally may not be valid for the prediction of liquefaction-induced ground surface disruption at sites susceptible to lateral spreading. The chart proposed in this study shows an improvement over the Ishihara's approach for predicting the liquefaction-induced ground surface damage. The microzonation maps comparing the liquefaction sites observed along the southern shore of Izmit Bay and in Yuanlin, and the surface damage and non-damage zones predicted from the proposed chart can identify accurately the liquefaction (sand boiling and lateral spreading) and no-liquefaction sites.     

Engineering geologic and geotechnical analysis of paleoseismic shaking using liquefaction effects: field examples

The greatest impediments to the widespread acceptance of back-calculated ground motion characteristics from paleoliquefaction studies typically stem from three uncertainties: (1) the significance of changes in the geotechnical properties of post-liquefied sediments (e.g., “aging” and density changes), (2) the selection of appropriate geotechnical soil indices from individual paleoliquefaction sites, and (3) the methodology for integration of back-calculated results of strength of shaking from individual paleoliquefaction sites into a regional assessment of paleoseismic strength of shaking. Presented herein are two case studies that illustrate the methods outlined by Olson et al. [Engineering Geology, this issue] for addressing these uncertainties.

The first case study is for a site near Memphis, Tennessee, wherein cone penetration test data from side-by-side locations, one of liquefaction and the other of no liquefaction, are used to readily discern that the influence of post-liquefaction “aging” and density changes on the measured in situ soil indices is minimal. In the second case study, 12 sites that are at scattered locations in the Wabash Valley and that exhibit paleoliquefaction features are analyzed. The features are first provisionally attributed to the Vincennes Earthquake, which occurred around 6100 years BP, and are used to illustrate our proposed approach for selecting representative soil indices of the liquefied sediments. These indices are used in back-calculating the strength of shaking at the individual sites, the results from which are then incorporated into a regional assessment of the moment magnitude, M, of the Vincennes Earthquake. The regional assessment validated the provisional assumption that the paleoliquefaction features at the scattered sites were induced by the Vincennes Earthquake, in the main, which was determined to have M7.5. The uncertainties and assumptions used in the assessment are discussed in detail.     

A Novel Methodology to Classify Soil Liquefaction Using Deep Learning

In this research, deep learning (DL) model is proposed to classify the soil reliability for liquefaction. The applicability of the DL model is tested in comparison with emotional backpropagation neural network (EmBP). The database encompassing cone penetration test of Chi–Chi earthquake. This study uses cone resistance (q_c) and peck ground acceleration as inputs for prediction of liquefaction susceptibility of soil. The performance of developed models has been assessed by using various parameters (receiver operating characteristic, sensitivity, specificity, Phi correlation coefficient, Precision–Recall F measure). The performance of DL is excellent. Consistent results obtained from the proposed deep learning model, compared to the EmBP, indicate the robustness of the methodology used in this study. In addition, both the developed model was also tested on global earthquake data. During validation on global data, both the models shows good results based on fitness parameters. The developed classification models a simple, but also efficient decision-making tool in engineering design to quantitatively assess the liquefaction potential. The finding of this paper can be further used to capture the relationship between soil and earthquake parameters.

     

Liquefaction hazard mapping at the town of Edessa,Northern Greece

The town of Edessa is located on Northern Greece at a region that is characterized as low seismicity zone due to the fact that few moderate events of M < 6 occurred during the last century. According to the Greek Seismic Code, the expected acceleration having a 10% probability of being exceeded in 50 years is equal to 0.16g. However, an amplification of ground motion is likely to occur due the local geology that is consisted of Holocene fluvio-torrential deposits. The basic aim of this paper is to evaluate the site amplification due to geological conditions and to assess the liquefaction hazard. In order to achieve this, 1-D site response analyses were performed. The data that were employed for the construction of the numerical models have been collected from borings with standard penetrations tests (SPT) that were drilled for construction purposes. Afterward, the liquefaction potential of the subsoil layers was evaluated taking into consideration two seismic scenarios. The first scenario was based on the seismic parameters, earthquake magnitude and PGA, assigned by the Greek Seismic Code. On the second seismic model, we employed the values of acceleration, resulted from the 1-D analyses and the earthquake magnitude as it was defined by the Greek Seismic Code. In order to compile the liquefaction hazard maps, we initially estimated the liquefaction potential index (LPI) of the soil columns using the parameters provided by SPT, for both seismic loadings, and afterward we correlated these values with the proposed classification of the severity of liquefaction-induced deformations. In addition, having computed the value of probability based on the LPI, liquefaction manifestations probability maps were compiled for both scenarios. The result of this study was that liquefaction-induced ground disruptions are likely to occur at the center of the city, among the branches of Voda River, only when the amplified values of acceleration are taken into account to the computation of liquefaction potential.