区域土壤地震液化预测简化方法

利用可大范围获取的空间参数给出区域土壤液化的评估方法,可预估震前各地区土壤液化可能性,可快速评估震后震区土壤液化情况,对预防减轻地震灾害以及地震快速响应救援都具有重要意义的.本文以我国1976唐山大地震液化调查资料为背景,提出了区域土壤液化预测的四参数简化评估方法,并用2011年新西兰基督城地震进行了检验.选取场地平均剪切波速V_s30、复合地形指数CTI、场地到河流的距离DR以及地面峰值加速度PGA,分别代表土壤密实程度、饱水状态、地质年代和遭遇的地震作用大小,并基于我国唐山地震液化区域调查资料生成的样本,运用经典二元Logistic回归方法,建立了区域土壤液化预测模型和评估公式,回归检验结果良好,整体回判成功率为77%,非液化和液化区域成功率分别为73%和78%.将公式应用于评估2011年新西兰基督城地震液化情况,区域液化和非液化成功率分别为82%和88%,总体成功率为84%.以上结果表明本文模型可靠,方法正确,可用于区域土壤液化震前预估以及土壤液化震后快速评估.

     

The use of liquefaction features in paleoseismology: Lessons learned in the New Madrid seismic zone, central United States

A recent study in the New Madrid seismic zone demonstrates that large uncertainties, often involved but rarely expressed, in paleoliquefaction studies can be reduced by conducting detailed investigations at the most promising sites for dating liquefaction features. During the site investigations, care must be taken to collect samples that will provide close maximum and minimum dates for liquefaction features. It is advisable to use two-sigma calibrated dates, rather than one-sigma calibrated dates or radiocarbon ages, when estimating ages of liquefaction features. Well-constrained ages of individual liquefaction features should provide the basis for estimating the timing of paleoearthquakes and correlating features across a region. As uncertainty in ages of liquefaction features decreases, confidence in estimates of timing, source areas and magnitudes of paleoearthquakes increases. The New Madrid study also shows that modern or historic earthquakes that induced liquefaction in the same region and whose locations and magnitudes are fairly well know can serve as calibration events for paleoearthquakes. Future efforts that could further improve the usefulness of liquefaction feature in paleoseismology include (1) the development of new techniques for dating liquefaction features directly, (2) case studies of modern earthquakes that focus on the site and spatial distributions of liquefaction feature as well as geotechnical properties of liquefaction sites and (3) more rigorous quantification of uncertainties associated with estimates of timing, source areas and magnitudes of paleoearthquakes.     

Anomalous radon emanation at local and regional distances preceding earthquakes in the New Madrid seismic zone and adjacent areas of the central mid-continent of North America, 1981–84

Anomalous soil-radon activity, including several spike-like surges over periods of 5, 2, and 2 1/2 months, and a year-long declining trend, preceded the most significant earthquakes of the central mid-continental region of North America during 1981 and early 1984. The 5-month period of erratic weekly radon activity, February–June 1981, was followed by a tremor of magnitudeM=4.0, 40 km from the monitoring site in the New Madrid seismic zone. An unusual earthquake swarm in central Arkansas, 160 km from the New Madrid seismic zone and 230 km from the monitoring site in the seismic zone, began in January 1982, shortly after a year-long declining trend in anomalous radon emanation. Earthquakes of magnitudes 4.5, 4.1, and 4.0 occurred at the swarm's outset in early 1982, fitting a pattern anticipated for the New Madrid seismic zone on the basis of the radon activity. Two periods of coincident peak radon emanation have since been observed in the Arkansas and New Madrid seismic regions, as have synchronous seismic pulses for the two separate areas.Two more recent periods of highly erratic soil-radon emanation, March–May 1983 and November–January 1984, were followed by a 4.3 earthquake in southwestern Illinois on 15 may 1983, and 3.5–3.6 tremors and swarm activity in the New Madrid seismic zone in late January and mid-February 1984. Prior to the 4.3 event, radon peaked at three widely separated monitoring sites 1–3 weeks before the tremor at distances of 120, 225, and 320 km from the epicentral region, the station at 225 km, in the New Madrid seismic zone, recording the longest period of anomalous radon activity. As for the recent 3.5–3.6 tremors of 1984, seismic activity of this magnitude had been anticipated for January or February on the basis of the amplitude of the November radon anomaly.These observations provide further evidence of (a) the existence of soil-radon anomalies precursory to the large earthquakes in this intraplate region, (b) the utility of such anomalies in anticipating events of small to moderate magnitudes for the region, and (c) the occurrence of regional-scale strain events prior to some of the larger mid-continental earthquakes.A very recent radon anomaly, the strongest yet to be detected in the seven years of monitoring in the mid-continental region, occurred in the New Madrid seismic zone from mid-February through mid-June 1984. A 4.0 earthquake occurred one month after a peak in the radon activity. The amplitude and duration of the anomaly suggest that a significant change in the state of stress or strain may have occurred in the mid-continental region during 1984.     

Numerical simulation of liquefaction effects on seismic SSI

The present paper deals with the influence of soil non-linearity, introduced by soil liquefaction, on the soil–foundation–structure interaction phenomena. The objective is to reveal the beneficial or unfavourable effects of the non-linear SSI on both structural drift and settlement of a given structure. Factors such as the signal modification due to liquefaction, and ratios of fundamental frequencies of soil, structure and signal may play an important role on the damage of the structure. The importance of each of these factors is evaluated through a significant parametric study. A 2D coupled finite element modelling is carried out using an elastoplastic multi-mechanism model to represent the soil behaviour. This paper presents the research work we did in the framework of the European Community project NEMISREF (New methods of mitigation of seismic risk on existing foundations, GRDI-40457), to study possible retrofitting measures using GEFDYN computational tools.     

Liquefaction analysis from seismic velocities and determination of lagoon limits Kumluca/Antalya example

This study analyzes liquefaction in the Kumluca/Antalya residential area and surroundings, using seismic velocities of soil deposits and the predominant period of the earthquake wave. The liquefaction analysis calculates shear–stress ratio, shear–resistance ratio and safety factor. Shear wave velocity used in liquefaction analysis was determined through surface waves. Moreover, the dynamic parameters of the ground were calculated through seismic velocities. Distributions of groundwater, shear wave velocity, adjusted shear wave velocity, predominant period of vibration, soil amplification and ground acceleration of the research area were mapped. In addition, the liquefied and non-liquefied areas as a result of liquefaction analysis in Kumluca were determined and presented as maps. Examining these maps, among all these maps, the limits of the lagoon sandbar and the old lake area were determined using only the liquefaction map.     

Threshold seismic energy and liquefaction distance limit during the 2008 Wenchuan earthquake

Estimating the possible region of liquefaction occurrence during a strong earthquake is highly valuable for economy loss estimation, reconnaissance efforts and site investigations after the event. This study identified and compiled a large amount of liquefaction case histories from the 2008 Wenchuan earthquake, China, to investigate the relationship between the attenuation of seismic wave energy and liquefaction distance limit during this earthquake. Firstly, we introduced the concept of energy absorption ratio, which is defined as the absorbed energy of soil divided by the imparted energy of seismic waves at a given site, and the relationship between the energy absorption ratio and the material damping ratio was established based on shear stress–strain loop of soil element and the seismic wave propagation process from the source to the site. Secondly, the threshold imparted seismic energy of liquefaction was obtained based on existing researches of absorbed energy required to trigger liquefaction of sandy soils and the ground motion attenuation characteristics of the 2008 Wenchuan earthquake, and the liquefaction distance limit of this earthquake was estimated according to the proposed magnitude–energy–distance relationship. Finally, the field liquefaction database of 209 sites of the 2008 Wenchuan earthquake was used to validate such an estimation, and the field observed threshold imparted seismic energy to cause liquefaction in recent major earthquakes worldwide was back-analyzed to check the predictability of the present method, and several possible mechanisms were discussed to explain the discrepancy between the field observations and the theoretical predictions. This study indicates that seismic energy attenuation and liquefaction distance limit are regional specific and earthquake dependent, and 382 J/m³ is the average level of threshold imparted seismic energy to cause liquefaction for loose saturated sandy soils, and the corresponding liquefaction distance limit is approximately 87.4 km in fault distance for a M_w?=?7.9 event in the Chengdu Plain. The proposed regional energy attenuation model and threshold imparted seismic energy may be considered as an approximate tool in evaluating the liquefaction hazard during potential earthquakes in this area.     

饱和砂土地震液化判别的分形插值模型

借鉴分形基本理论,提出了基于分形插值模型的饱和砂土地震液化判别方法.该方法首先选取影响饱和砂土地震液化判别的7个主要因素,根据分类标准,采用在每级标准中随机内插的方法,得到40个标准样本,用于构建饱和砂土地震液化判别的分形插值模型;其次根据最大似然分类原则确定每个饱和砂土地震液化判别指标的评价分维数;然后利用加权求和法计算样本的综合评价值,并根据样本综合评价值与经验等级之间的关系建立分形插值评价模型;最后,进行了实例分析结果表明:该模型的评价结果合理、客观,计算得到的每个样本具体得分值,即使对属于同一级的样本也可以给出其地震液化程度的顺序,为饱和砂土地震液化评价工作提供了一种新的研究方法与思路.     

Probabilistic seismic performance and loss assessment of a bridge–foundation–soil system

This paper presents the probabilistic seismic performance and loss assessment of an actual bridge–foundation–soil system, the Fitzgerald Avenue twin bridges in Christchurch, New Zealand. A two-dimensional finite element model of the longitudinal direction of the system is modelled using advanced soil and structural constitutive models. Ground motions at multiple levels of intensity are selected based on the seismic hazard deaggregation at the site. Based on rigorous examination of several deterministic analyses, engineering demand parameters (EDP's), which capture the global and local demand, and consequent damage to the bridge and foundation are determined. A probabilistic seismic loss assessment of the structure considering both direct repair and loss of functionality consequences was performed to holistically assess the seismic risk of the system.It was found that the non-horizontal stratification of the soils, liquefaction, and soil–structure interaction had pronounced effects on the seismic demand distribution of the bridge components, of which the north abutment piles and central pier were critical in the systems seismic performance. The consequences due to loss of functionality of the bridge during repair were significantly larger than the direct repair costs, with over a 2% in 50 year probability of the total loss exceeding twice the book-value of the structure.     

Assessment of variability and uncertainties effects on the seismic response of a liquefiable soil profile

The present work deals with the evaluation of the effect of the randomness of both soil parameters and input motion on the seismic response of a sandy soil profile. Special attention is given to estimate the relative contribution of model input parameters on (i) the probability of liquefaction apparition and (ii) the surface seismic response. The Monte Carlo simulations were used for that purpose. This analysis shows that, for the considered case, the choice of the earthquake input signal remains the most subtle parameter in order to define the liquefaction probability. Otherwise, spatial variability of both soil properties and soil model's parameters have a weak impact on the response of the soil profile.     

加固地基增强地基的抗震性能

加固地基是增强地基抗震性能的一种行之有效的方法。加固地基可提高地基抗液化能力，改善地基土的动力特性，减轻砂土地基的液化势；通过对地基土产生预震效应，提高了加固地基整体的强度、刚度和稳定性。对震害严重的厚冲积层或厚软土层，用高压喷射注浆法进行加固处理有其独特的优越性。     

Bounding surface plasticity model for the seismic liquefaction analysis of geostructures

This paper presents the constitutive relations and the simulative potential of a new plasticity model developed mainly for the seismic liquefaction analysis of geostructures. The model incorporates the framework of critical state soil mechanics, while it relies on bounding surface plasticity with a vanished elastic region to simulate the non-linear soil response. Key constitutive ingredients of the new model are: (a) the inter-dependence of the critical state, the bounding and the dilatancy (open cone) surfaces on the basis of the state parameter ψ, (b) a (Ramberg–Osgood type) non-linear hysteretic formulation for the “elastic” strain rate, (c) a discontinuously relocatable stress projection center related to the “last” load reversal point, which is used for mapping the current stress point on model surfaces and as a reference point for introducing non-linearity in the “elastic” strain rate and finally (d) an empirical index of the directional effect of sand fabric evolution during shearing, which scales the plastic modulus. In addition, the paper outlines the calibration procedure for the model constants, and exhibits its accuracy on the basis of a large number of laboratory element tests on Nevada sand. More importantly, the paper explores the potential of the new model by presenting simulations of the VELACS centrifuge tests of Models No 1 and 12, which refer to the free-field liquefaction response of Nevada sand and the seismic response of a rigid foundation on the same sand, respectively. These simulations show that the new model can be used successfully for the analysis of widely different boundary value problems involving earthquake soil liquefaction, with the same set of model constants calibrated on the basis of laboratory element tests.     

A surface seismic approach to liquefaction

The liquefaction potential of soils is traditionally assessed through geotechnical approaches based on the calculation of the cyclical stress ratio (CSR) induced by the expected earthquake and the ‘resistance’ provided by the soil, which is quantified through standard penetration (SPT), cone penetration (CPT), or similar tests. In more recent years, attempts to assess the liquefaction potential have also been made through measurement of shear wave velocity (V_S) in boreholes or from the surface. The latter approach has the advantage of being non-invasive and low cost and of surveying lines rather than single points. However, the resolution of seismic surface techniques is lower than that of borehole techniques and it is still debated whether it is sufficient to assess the liquefaction potential.In this paper we focus our attention on surface seismic techniques (specifically the popular passive and active seismic techniques based on the correlation of surface waves such as ReMi^TM, MASW, ESAC, SSAP, etc.) and explore their performance in assessing the liquefaction susceptibility of soils. The experimental dataset is provided by the two main seismic events of M_L=5.9 and 5.8 (M_W=6.1, M_W=6.0) that struck the Emilia-Romagna region (Northern Italy) on May 20 and 29, 2012, after which extensive liquefaction phenomena were documented in an area of 1200 km².The CPT and drillings available in the area allow us to classify the soils into four classes: A) shallow liquefied sandy soils, B) shallow non-liquefied sandy soils, C) deep non-liquefied sandy soils, and D) clayey–silty soils, and to determine that on average class A soils presented a higher sand content at the depth of 5–8 m compared to class B soils, where sand was dominant in the upper 5 m. Surface wave active–passive surveys were performed at 84 sites, and it was found that they were capable of discriminating among only three soil classes, since class A and B soils showed exactly the same V_S distribution, and it is possible to show both experimentally and theoretically that they appear not to have sufficient resolution to address the seismic liquefaction issue.As a last step, we applied the state-of-the art CSR–V_S method to assess the liquefaction potential of sandy deposits and we found that it failed in the studied area. This might be due to the insufficient resolution of the surface wave methods in assessing the Vs of thin layers and to the fact that Vs scales with the square root of the shear modulus, which implies an intrinsic lower sensitivity of Vs to the shear resistance of the soil compared to parameters traditionally measured with the penetration tests. However, it also emerged that the pure observation of the surface wave dispersion curves at their simplest level (i.e. in the frequency domain, with no inversion) is still potentially informative and can be used to identify the sites where more detailed surveys to assess the liquefaction potential are recommended.     

Numerical simulation on seismic liquefaction by adaptive mesh refinement due to two recovered fields in error estimation

The finite element method is widely employed in numerical analysis of seismic liquefaction at present. However, due to the mesh distortion in strong dynamic action, the results are generally inaccurate and can even stop the calculation. Therefore, an adaptive mesh refinement (AMR) scheme is introduced to improve the numerical results based on a finite element and finite difference coupled dynamic method (FEM–FDM). The changing law about the errors due to the linear and the bilinear recovered fields of coarse meshes, AMR meshes and fine meshes is discussed. The seismic response on soil is compared between different AMR meshes. The AMR strategy can better simulate seismic liquefaction and offer more accurate results than the normal finite element simulation.     

Liquefaction in silty soils—screening and remediation issues

Current techniques for liquefaction screening, ground modification for liquefaction mitigation, and post-improvement verification rely on knowledge gained from extensive research on clean sands, field observations of liquefied ground, and judicial correlation of normalized penetration resistance [(N₁)₆₀,q_c1N] or shear wave velocity (v_s1) data with field liquefaction observations. Uncertainties prevail on the direct extrapolation of such techniques for silty soil sites. This paper examines laboratory data on liquefaction resistance, strength, and v_s1 of sands and silty soils using grain contact density as the basis. Effect of silt content on cyclic resistance, strength, m_v, and c_v is examined in this light. Rational insights on effects of silt content on the current screening techniques based on (N₁)₆₀, q_c1N, and v_s1 to silty soils are offered. Recent advances and modifications to the traditional densification, drainage, and permeation grouting techniques to make them viable for silty soils are discussed.     

Behavior of stratified sand–silt–gravel composites under seismic liquefaction conditions

Laboratory seismic liquefaction studies have generally dealt with homogeneous soil conditions only, although stratified soils exist for various soil deposits. The main objective of this research project was to compare the behavior of stratified and homogeneous sand–silt–gravel composites during seismic liquefaction conditions for various silt and gravel contents. An experimental program was undertaken in which a total of eighty stress-controlled undrained cyclic triaxial tests were performed. Two methods of sample preparation were used for each soil type. These methods included moist tamping (representing homogeneous soil conditions) and sedimentation (representing layered soil conditions). The silt contents ranged from 0 to 50%, and soils with 10 and 30% gravel contents were tested. The confining pressure in all test series was 100 kPa. The results indicate that the liquefaction resistances of layered and uniform soils are not significantly different, despite the fact that the soil fabric produced by the two methods of sample preparation is totally different. This finding justifies applying the laboratory tests results to the field conditions for the range of variable studied.     

汶川地震砾性土液化场地特征解析

通过成都平原砾性土场地勘察测试,研究汶川地震中大量砾性土液化场地的基本特性,找出一般规律,对砾性土场地液化发生主客观原因提出解释,并修正以往若干认识偏差.分析表明：汶川地震液化砾性土层粒径范围宽,含砾量5%~85%甚至更大,同时其实测剪切波速140~270 m·s^-1,修正剪切波速160~314 m·s^-1,都远超历史记录;液化砾性土场地1/2集中在Ⅷ度区内,表明如砂土层液化一样,砾性土场地大规模液化需要较强地震动触发,但超过触发强度后液化规模增长均有限;成都平原浅表地层二元基本结构是汶川地震中出现大量砾性土场地的客观条件之一,该结构可使饱和砾性土层处于封闭状态,构成了砾性土液化的基本条件;虽然液化砾性土层剪切波速很高,但实际上大多松散状态,是此次地震大量砾性土场地发生液化的客观条件之二;地震中地表（井中）喷出物与地下实际液化土类大相径庭,且液化层埋深大多小于6.0 m,以往以地表喷出物反推地下液化层土性类型的做法不再成立;认为砾性土层波速大、透水性好而不会液化的传统认识也不再成立,但砾性土层液化条件与砂土层液化条件不同,前者要求更高.

     

An empirical method to estimate shear wave velocity of soils in the New Madrid seismic zone

In this study, a set of charts are developed to estimate shear wave velocity of soils in the New Madrid seismic zone (NMSZ), using the standard penetration test (SPT) N values and soil depths. Laboratory dynamic test results of soil samples collected from the NMSZ showed that the shear wave velocity of soils is related to the void ratio and the effective confining pressure applied to the soils. The void ratio of soils can be estimated from the SPT N values and the effective confining pressure depends on the depth of soils. Therefore, the shear wave velocity of soils can be estimated from the SPT N value and the soil depth. To make the methodology practical, two corrections should be made. One is that field SPT N values of soils must be adjusted to an unified SPT N′ value to account the effects of overburden pressure and equipment. The second is that the effect of water table to effective overburden pressure of soils must be considered. To verify the methodology, shear wave velocities of five sites in the NMSZ are estimated and compared with those obtained from field measurements. The comparison shows that our approach and the field tests are consistent with an error of less than of 15%. Thus, the method developed in this study is useful for dynamic study and practical designs in the NMSZ region.     

Dynamic properties and liquefaction resistance of two soil materials in an earthfill dam—Laboratory test results

Results are presented of laboratory resonant column and cyclic triaxial tests on specimens of two compacted soils (a sandy–silty clay and a sand–gravel mixture), planned to be used in the core and the shells, respectively, of a proposed earthfill dam. The values of low-amplitude shear modulus of the clayey material were found to increase with increasing confining pressure and decreasing water content, with deviations of ±20% from the predictions of the “Hardin equation”. On the other hand, the low-amplitude damping ratio was found to be at least four times higher than the values corresponding to natural undisturbed cohesive soils. The proposed G/G₀−γ_c curve for the compacted cohesive soil was found to be independent of confining pressure and small variations of the water content on either side of the optimum value and showed a remarkable agreement with recently published similar curves for natural cohesive soils. For the case of the sand–gravel mixture, normalized G/G₀–γ_c and D–γ_c curves are proposed, based on recently published results for gravelly soils and the limited data of the present study. The liquefaction resistance of the saturated sand–gravel mixture was found to be strongly dependent on its relative density, especially for high values of cyclic stress ratio. A relative density of at least 55% was found to be necessary to assure safety against earthquake-induced liquefaction of the material. The results presented herein may be used (directly or as guide) in the seismic analysis of (new or existing) earth dams constructed from similar soil materials and in addition they provide insight into the dynamic behavior of compacted soils.     

Development of regional simulation of seismic ground‐motion and induced liquefaction enhanced by GPU computing

Much research has been conducted for physics‐based ground‐motion simulation to reproduce seismic response of soil and structures precisely and to mitigate damages caused by earthquakes. We aimed at enabling physics‐based ground‐motion simulations of complex three‐dimensional (3D) models with multiple materials, such as a digital twin (high‐fidelity 3D model of the physical world that is constructed in cyberspace). To perform one case of such simulation requires high computational cost and it is necessary to perform a number of simulations for the estimation of parameters or consideration of the uncertainty of underground soil structure data. To overcome this problem, we proposed a fast simulation method using graphics processing unit computing that enables a simulation with small computational resources. We developed a finite‐element‐based method for large‐scale 3D seismic response analysis with small programming effort and high maintainability by using OpenACC, a directive‐based parallel programming model. A lower precision variable format was introduced to achieve further speeding up of the simulation. For an example usage of the developed method, we applied the developed method to soil liquefaction analysis and conducted two sets of simulations that compared the effect of countermeasures against soil liquefaction: grid‐form ground improvement to strengthen the earthquake resistance of existing houses and replacement of liquefiable backfill soil of river wharves for seismic reinforcement of the wharf structure. The developed method accelerates the simulation and enables us to quantitatively estimate the effect of countermeasures using the high‐fidelity 3D soil‐structure models on a small cluster of computers.     

1976年唐山大地震CPT液化数据库重构

1976年唐山大地震导致了范围广、灾害严重的液化震害。地震发生后,已有学者对液化场地进行了2次静力触探测试。本文首先给出基于2次测试数据的唐山大地震CPT液化数据库,利用中国规范方法和NCEER推荐方法对数据库进行液化判别,发现针对第1次测试数据的判别成功率较高,而针对第2次测试数据的判别成功率较低。经30年的时间,绝大多数测试点的液化层强度与埋深均发生了较大变化,第2次测试时土层液化可能性已发生了较大变化,其测试数据对液化判别方法的改进意义较小。利用第2次测试的摩阻比Rf将第1次测试数据ps分解为qc和fs,相比第2次测试数据,分解指标具有更高的可靠性。为此基于第1次测试分解数据构建唐山大地震CPT液化数据库,为液化判别方法改进提供数据支持。