断层错动引起的上覆土体破裂演化规律研究

断层引起的地面永久大变形是工程特别是生命线工程地震破坏的重要原因之一, 而研究断层错动下上覆土体变形和破裂的发展演化规律, 则是预测地面永久变形状态和分析断层危害性的基础． 本文采用有限元方法对垂直断层错动引起的上覆土体破裂演化规律进行了研究, 建立了垂直断层作用下上覆土体模拟的有限元模型, 对断层错动作用下上覆土体的破裂发展过程进行了模拟分析, 并分析了加载速率、 土体特性及断层倾角等参数对上覆土体的破裂演化发展过程的影响． 结果表明: ① 断层倾角越陡, 地表出现破裂时需要增加的垂直位移越大; ② 由于惯性力的影响, 断层加载速率对地表破裂所需施加位移和土层破裂角产生一定的影响; ③ 断层类型对土层地表破裂角与膨胀角、 摩擦角之间的关系有很大影响． 该分析结果可为新建工程的抗震设计和已建工程结构的抗震加固等工作提供依据．      

Centrifuge modeling of interaction between reverse faulting and tunnel

In this study, a series of centrifuge tests, modeling reverse fault rupture with 60° dip angle, were conducted in a dry sandy soil with a tunnel embedded in the soil layer. The test results showed that the tunnel and soil responses depended on the tunnel position, soil relative density and tunnel rigidity. Tunnels appeared be able to deviate the fault rupture path, while this deviation may be associated with significant rotation and displacement of the tunnel. However, a deeper tunnel was able to diffuse the shear deformation to a wider zone with an unsmooth surface displacement which may cause severe damage to the surface structures. Finally, the tunnel rotation, the location of the fault outcropping, the vertical displacement of the ground surface, the effect of tunnel rigidity on fault rupture path and surface displacement and the effect of soil relative density on fault–tunnel interaction were reported and discussed in this study.     

Influence factors on the ground surface rupture zone induced by buried normal fault dislocation*

The seismic disaster presents a zonal distribution along the fault strike. In this paper, rupture zone of ground surface soil caused by the uniform dislocation, inclined dislocation and warped dislocation of buried normal fault are studied by constituting a three-dimensional finite element model in Automatic Dynamic Incremental Nonlinear Analysis (ADINA). According to the critical value of surface rupture, the variational features and influencing factors of width and starting position of the “avoiding zone” in engineering construction are analyzed by using 96 model calculations. The main results are as follows: (1) Since the rupture zone of the ground surface soil from the point of mechanics is different from the “avoidance zone” from the point of engineering safety, the equivalent plastic strain and the total displacement ratio should be considered to evaluate the effect of the seismic ground movement on buildings. (2) During fault dislocation, plastic failure firstly occurred on the ground surface soil of the footwall side, and then the larger deformation gradually moved to the side of the hanging wall of the fault with the increase of fault displacement. (3) When the vertical displacement of buried fault reaches 3 m, the width of “avoiding zone” in engineering construction varies within the range of 10–90 m, which is most affected by the thickness of overlying soil and the dip angle of the fault.     

震源破裂方式和断层性质对近场强地震动特征的影响

基于数值格林函数法的近场强地震动数值模拟方法,以1994年Northridge地震断层面上位错量的不均匀分布模式和该地区的地层剪切波速度结构为震源模型和计算模型,做了两个方面的模拟研究:(1)直立走滑断层(断层倾角为90°)情况下,模拟分析了有限断层单侧破裂模式和双侧破裂模式对强地震动特征——破裂方向性和上盘效应的影响;(2)对于倾斜断层(倾角为45°),模拟分析了正断层、逆断层和走滑断层情况下,单侧破裂模式对其强地震动主要特征——破裂方向性效应和上盘效应的影响.结果表明:断层的破裂方式直接影响着地表地震动峰值和矢量分布;在近场区无论直立断层还是倾斜断层,其地表地震动峰值分布所表达的破裂方向性效应显著,位于破裂传播前方的地震动强度大,反映了波前被压缩的趋势,破裂后方地震强度明显变小;倾斜断层引起的上盘效应明显,NS向分量和竖向分量的地表地震动峰值的最大值出现在上盘靠近断层迹线处,EW向分量的峰值在断层迹线两侧呈不对称分布,且逆断层引起的地震动峰值最大,走滑断层的次之,正断层的最小.     

Direct estimation of rupture depths of earthquake faults from coseismic surface deformation

The rupture dimensions of earthquake faults are important parameters for characterizing earthquake ruptures and ground motions. Two key parameters to be determined are the rupture depth and dip angle of earthquake faults. Dislocation theory in an elastic half space indicates that if a seismic rupture directly runs up to the ground surface, there exist zero points of horizontal strain in the surface deformation, which correspond to the rupture depths, except for pure strike-slip faults. In this study, we use numerical simulations to investigate the possibility of inferring rupture depths from zero-strain points for cases of buried faults and heterogeneous media. The results show that the correspondence of zero-strain points to the rupture depths can be influenced by the heterogeneity of the underground media and the stress field. For buried faults, the correspondence relationship is approximately valid when the fault depth is <1 km. In addition, the range of earthquake fault dip angles can be estimated by horizontal displacements on the ground. We also study how to determine the rupture depths of faults from InSAR data after large earthquakes, and successfully apply the method to the 2008 Wenchuan earthquake. The method proposed here, which determines the parameters of fault geometry according to surface deformation, is simple and easy to perform. With independent of aftershocks, it can provide valuable constraints to kinematic inversions.     

隐伏正断层错动引发上覆土体破裂过程的三维数值模拟

通过建立三维数值模型,对隐伏正断层在均匀错动和倾斜错动方式下土体的破裂过程进行研究。利用应力罗德参数和等效塑性应变分别对断层错动过程中上覆土体的应力状态和破坏形式进行分析,并提出土体破裂的判别方法。通过对数值模拟结果的分析得到以下结论:① 在断层错动过程中,下盘一侧受断层错动影响的上覆土体的应力状态经压剪→纯剪→拉剪逐渐变化,而上盘一侧上覆土体的应力状态变化较为复杂,经压剪→纯剪→拉剪→纯剪→压剪重复变化;② 在断层均匀错动过程中,断层下盘一侧土体的破裂率先出现在地表拉剪区内,随错动量的增大,破裂带向两侧、向深部扩展;同时,下盘一侧土体的底部产生破坏,并斜向上扩展,逐渐与顶部破裂相连;③ 在断层倾斜错动过程中,地表破裂出现的位置和上覆土体的厚度有关。对于厚度较大的土体,正断层倾斜错动能够在地表形成与断层走向有一定夹角、且与断层长度相比长度很短的地表破裂或地裂缝,而数值模拟可对正断层错动导致的地表破裂的模式加以补充,为研究地裂缝的形成机理和分布形式提供依据。      

基于有限断层动力学模型的逆断层—三维层状沉积盆地地震动模拟

近断层效应使得沉积盆地对地震动放大效应更为复杂。本文针对逆断层发震下三维层状沉积盆地地震反应,基于波动谱元法,采用有限断层动力学模型,模拟断层动力破裂、地壳层地震波传播和层状沉积盆地对地震波散射全过程。在此基础上,对比分析了层状和均质沉积盆地对近断层地震动放大效应的影响,讨论了不同断层倾角下层状沉积盆地地震动加速度特性。结果表明：层状沉积盆地PGA空间分布与均质沉积盆地存在较大差异,由于近断层效应和盆地效应,层状沉积盆地地表局部范围竖向PGA大于水平向PGA;90°断层倾角下层状沉积盆地地表地震动放大范围与60°断层倾角结果明显不同,主要集中在盆地中心区域和断层附近,且幅值远小于60°断层倾角下结果;沿断层走向,盆地内地表地震动加速度峰值对应时刻较盆地外延后。     

隐伏正断层错动致地表破裂变形特征的研究

通过建立三维计算模型,对隐伏正断层在均匀错动、倾斜错动和翘倾错动方式下地表土体的应力路径、破裂和变形特征进行了研究。根据地表破裂临界值,分析了工程建设“避让带”的宽度和起始位置的变化特征。根据行业规范,提出工程建设“关注带”的确定方法,分析了“关注带”的宽度和起始位置的变化特征,得到以下主要结论:①在断层错动过程中,位于两侧的地表土体应力路径变化明显不同,下盘一侧和上盘一侧分别以三轴拉伸和三轴压缩为主;②地表强变形带与地表破裂带的分布并不一致,需要综合考虑等效塑性应变和总位移比2个指标来评价同震地表错动对建筑物的影响;③当隐伏断层错动的垂直位移达到3m时,工程建设“避让带”的宽度在10—90m范围内变化,受上覆土体厚度和断层倾角的影响最大,而工程建设“关注带”的宽度在150—400m范围内变化,受上覆土体的性质影响最大。     

Response analysis of buried pipelines crossing fault due to overlying soil rupture

A 3-D soil-pipe nonlinear finite element model with contact element is suggested and the influences of the rupture mode, thickness and rigidity of overlying soil on the response of buried pipeline are analyzed. The numerical results show that the soil rupture mode determines the location of the large deformation or failure of the pipeline, and the plastic deformation of the pipeline occurs at the zone where the plastic deformation or rupture of the overlying soil appears. When the fault dip angle on bedrock...     

覆盖土层场地地震断裂反应分析方法

基于有限元和拟静力的根本原理，提出了一种基于断层位错的覆盖土层场地地震断裂反应弹塑性有限元分析方法，计算中考虑了土层厚度、断层种类和倾角、软弱夹层以及土的非线性等因素，通过编制的程序计算了地表断裂位移，模拟了土层的破裂过程。     

1999年台湾集集大地震的地表断层破裂特征

野外考察结果表明 ,1999年 9月 2 1日台湾集集大地震是由车笼埔断层发生逆冲作用造成的。地震产生的地表破裂长约 80km ,具有明显的挤压逆断层特征 ,其活动方式为具左旋性质的逆倾滑动。实测逆冲断层以 30°～ 50°的角向西北逆冲而上。断层的垂直位移量 ,南段约 2～ 3m ,北段约 3～ 8m ;断层的水平位移量 ,南段 0～ 3m ,北段 3～ 5m ;垂直断层的水平缩短量 ,南段 2～ 3m ,北段 3～ 6m。从台湾西部麓山带的地质构造剖面分析 ,地震震源恰好位于台湾西部麓山带中生代基底与其上的沉积盖层的界面的深度 ,而西部麓山带第三纪地层和其下的基底的分离面为一滑动面 (decolle ment)。在菲律宾海板块的挤压作用下 ,沿该区中生代基底之上滑动面的错动导致了地震的发生     

应用动态复合震源模型模拟汶川Mw7.9地震强地面运动

2008年5月12日中国汶川地区发生Mw7.9地震,震中位置103.4°E,31.0°N.主要发震断层空间展布长达300多公里,由南西方向到北东方向呈现明显的分段性,汶川一映秀段逆冲为主兼有少量的右旋走滑分量;安县一北川段为逆冲一右旋走滑的断层错动;青川段以右旋走滑为主兼有少量逆冲分量.采用改进后的复合震源强地面运动预...     

利用断层围陷波资料研究汶川M_S8.0地震构造特征

利用汶川地震区不同地段的断层围陷波记录,分析了该地震断层的分段性特征。对断层北东段的关庄测线分析研究结果表明：地壳内破碎带的宽度大约160~180m,地下破碎带的中间与地表破裂的位置对应,并且地下破碎带在断层的两盘边缘较均匀地分布,反映了北东段的断层倾角较陡,近似直立断层。对断层南西段的虹口测线研究结果表明：地壳内破碎带的宽度大约180~200m,地下破碎带主要分布在地表断层陡坎上盘所对应的地壳内,反映了南西段断层倾角比北东段断层倾角小。本文的研究结果可以为汶川8.0级地震的构造背景研究提供依据。     

剪切破裂源及其近场影响

本文利用有限差分法对垂直走滑断层和垂直倾滑断层两种剪切破裂源在不同埋深和不同破裂尺度(指垂直于地面方向上断层尺度)情况下的近场理论地震图进行了讨论,并得到了地表不同震中距处的加速度付氏谱。通过对上述两种源所产生的近场效应的对比分析,本文得到了以下初步结果: (1) 在埋深相同,断层面宽度相同的情况下,垂直走滑断层源对近场地表的影响比垂直倾滑断层源大。(2) 在垂直倾滑断层源的埋深、破裂尺度相同的情况下,在地表同一震中距处产生的水平位移要比垂直位移大。(3) 两种破裂源在地表同一震中距处产生的位移的最大振幅随破裂尺度的增大而增大,且近似为线性关系;在极震区内不同震中距处产生的位移的最大振幅随震中距的增大而减小,且随震中距的变化率在减小;在地表不同震中距处产生的加速度的付氏谱的主峰周期随震中距的增大而变长,随破裂的增大,付氏谱中主峰周期也加长,且长周期成份愈多;随着破裂源埋深的增加,地表不同震中距处产生的位移的最大振幅变小,即等震线变疏。(4) 走滑断层源产生的加速度付氏谱的主峰周期比倾滑断层源长。     

Numerical analyses of fault–foundation interaction

Field evidence from recent earthquakes has shown that structures can be designed to survive major surface dislocations. This paper: (i) Describes three different finite element (FE) methods of analysis, that were developed to simulate dip slip fault rupture propagation through soil and its interaction with foundation–structure systems; (ii) Validates the developed FE methodologies against centrifuge model tests that were conducted at the University of Dundee, Scotland; and (iii) Utilises one of these analysis methods to conduct a short parametric study on the interaction of idealised 2- and 5-story residential structures lying on slab foundations subjected to normal fault rupture. The comparison between numerical and centrifuge model test results shows that reliable predictions can be achieved with reasonably sophisticated constitutive soil models that take account of soil softening after failure. A prerequisite is an adequately refined FE mesh, combined with interface elements with tension cut-off between the soil and the structure. The results of the parametric study reveal that the increase of the surcharge load q of the structure leads to larger fault rupture diversion and “smoothing” of the settlement profile, allowing reduction of its stressing. Soil compliance is shown to be beneficial to the stressing of a structure. For a given soil depth H and imposed dislocation h, the rotation Δθ of the structure is shown to be a function of: (a) its location relative to the fault rupture; (b) the surcharge load q; and (c) soil compliance.     

近场强地震动数值模拟的简化计算方法

近场强地震动除受场地条件的影响外,还受到震源破裂面上子源的空间分布特点、子源破裂先后顺序的强烈控制,基于数值格林函数法的近场强地震动数值模拟方法可以综合考虑震源、传播途径及局部场地条件的影响,对计算过程进行合理简化,分2步完成地震动模拟:第1步,在介质均匀区采用矩张量的解析解计算所有子源在盖层底面的位移,形成下一步有限元计算的输入场;第2步,在盖层介质不均匀区,结合局部人工透射边界技术,采用时、空解耦的波动显式有限元方法计算地表强地震动。在有限断层模型中,采用具有9个力偶的等效地震矩张量表达断层产状、滑动方向等的影响,采用Brune模型定义各子源的滑动时间函数,描述滑动的时、空不均匀分布特征,从而细化震源模型。通过对Northridge地震中4个基岩台站地表地震动的模拟结果和强震记录,验证了此简化计算方法的可行性     

基于FLAC3D数值模拟的清徐地裂缝与交城断裂成因关系研究

基于断裂力学理论,应用FLAC3D数值模拟方法对交城断裂活动时附近土体的变形特征进行模拟计算。计算结果表明:断层活动带动上覆土层差异沉降及应力场变化,当断层错动量达到一定程度,断层附近出现拉应力区,土体在拉张应力和自身重力作用下产生垂直差异形变,从而产生地裂缝。将数值模拟结果与方山探槽揭示的地裂缝剖面特征进行对比,二者结果相互印证,说明清徐地裂缝是交城断裂活动的结果,二者具有明显的对应关系,地裂缝为构造成因裂缝。数值模拟结果为研究断裂和地裂缝成因关系提供了佐证。     

Earthquake fault rupture—shallow foundation interaction in undrained soils: a simplified analytical approach

The interaction between normal or reverse fault ruptures and shallow foundations resting on a homogeneous undrained soil layer is investigated. After performing a thorough set of non‐linear finite element simulations, three simple kinematic mechanisms are proposed to predict the conditions for the fault rupture to be diverted by the foundation. The results obtained by both numerical and analytical approaches are in good agreement, and support the adequacy of the proposed mechanisms. A relationship is established for assessing the minimum foundation bearing load needed for diversion of the fault rupture trace, either normal or reverse. This relationship is proved to be independent of the fault type and dip angle, so that it can be applied easily for engineering applications, even in the absence of specific tectonic information. Copyright © 2006 John Wiley & Sons, Ltd.     

走滑断层地震地表破裂带分布影响因素数值模拟研究——以1973年炉霍MS7.6地震为例

地震后在断层两侧的强变形与破裂带是地震灾害最严重的区域.为系统、定量研究同震地表变形带特征及其影响因素,本研究建立了走滑断层的三维有限元模型,分别探讨了断层位错量、断层倾角、错动方式、上覆松散层厚度、沉积层土性等因素的影响规律.模拟结果显示：走滑断层同震地表变形表现为以断层为中心的近似对称单峰分布,强地表变形集中在断层两侧各50 m宽度范围,地表变形量峰值随位错量增加而增大,破裂带宽度也随位错量增加而增大,但增量逐渐减小,并趋于一个渐近值;断层倾角对地表变形与破裂带宽度影响表现为随倾角减小变形量峰值点向上盘小距离偏移;走滑兼正断位错引起的变形量峰值最大,但地表破裂带宽度最小,走滑兼逆断引起的变形量峰值最小,但地表破裂带宽度最大,直立纯走滑断层的两参量都居中;走滑断层地表变形量峰值随上覆松散层厚度增大而减小,但随厚度减小的速率逐渐变小,松散层厚度从5 m增加到20 m时,破裂带宽度随厚度增加而缓慢增加,但自厚度大于20 m时,破裂带宽度开始随厚度增加而逐渐下降;当不同土性覆盖层（粗砂、粉砂、黏土）厚度相同时,地震引起的地表变形量峰值自粗砂、粉砂、黏土逐次增大,当粗砂厚度为60 m以上时,3.6 m的同震水平位错已不能形成地表破裂,而粉砂的厚度为70 m以上,黏土的厚度则为75 m以上.     

ESTIMATING THE WENCHUAN EARTHQUAKE RUPTURE PROCESS USING NEAR FIELD STRONG MOTION RECORDS AND COSEISMIC DISPLACEMENTS

Based on the extensive near field broadband strong-motion records with uniform azimuthal coverage and coseismic displacements, the rupture process of 2008 Wenchuan earthquake is inversed by the non-negative least square method and multiple-time window technique. The possible rupture sequence among southern Beichuan Fault, Pengguan Fault and Xiaoyudong Fault and the initial rupture time of high dip angle part of southern Beichuan Fault are analyzed from kinetic aspects, which have been seldom focused on. The results indicate that:(1) The near field waveform fitting residuals and the coseismic displacements show that only a bilateral rupture occurs on the intersection between the southern Beichuan Fault and Xiaoyudong Fault can the synthetic records of the stations located near the southwestern end of the Beichuan Fault conform to the observed ones, and meanwhile, the Pengguan fault cannot generate large slips on its southwestern part. The possible rupture sequence is that the earthquake started at the low dip angle part of Beichuan Fault and propagates to the Pengguan Fault in the shallow area, the Xiaoyudong Fault is triggered by the Pengguan Fault, and then producing bilateral rupture on the high dip angle part of Beichuan Fault at the intersection with the Xiaoyudong Fault. (2) Through analysis of the synthetic second packet records of stations at the southwest area of the fault, we obtain the initial rupture time on the high dip angle part of Beichuan Fault may have a 8s stagnation. In terms of timing, there may be rupture sequence between the southern Beichuan Fault and Pengguan Fault which are parallel to each other. The rupture of the southern shallow part of Beichuan Fault with high dip angle may lag behind the Pengguan Fault. At the same time, there may be a multipoint rupture in the southern section of the Beichuan Fault. (3) There is a good correspondence between the area on the fault with larger slip rate and the surrounding stations with larger PGV. In areas where slip rate on the fault plane is large, the stations tend to have larger peak ground velocities.