鲜水河断裂带的地震模拟

本文运用由弹簧—滑块模型来模拟走滑断层的地震活动,讨论了模型中断层强度、断层滑块间弹簧强度、滑块与驱动盘之间片簧强度、滑块摩擦力的影响,以及初始条件敏感性。并以我国最为典型的一条走滑断层———鲜水河断裂带为例,运用试错法调整模型参数以使发生的地震最大限度地逼近实际发生的地震,并利用最终结果讨论了鲜水河断裂带今后的活动性。     

中国大陆强震轮回活动的计算机模型研究

通过对中国大陆7级以上强震的研究发现,地震活动呈现为时间上轮回和空间上迁移的动态图像。为探索这种图像的物理本质,本文应用非线性动力学模型来模拟在统一构造应力场作用下含有多个活动地震带的地震构造块体的地震活动。文中用若干列并联排列的基本元件(由Maxwell体和刚性滑块组成)来模拟地震活动带,并用耦合元件来模拟地震活动带之间和带内各段之间的相互耦合作用。本文的研究结果有助于理解强震的轮回活动以及高潮中成组强震间的相互影响特征。     

大陆地震的动力学模型研究

应用非线性动力学模型，以大陆构造构体为参照对象，对大陆地震的成组孕育和发生过程进行模拟研究，通过对６×８弹簧－滑块－阻尼器组合单元非线性动力学模型的理论计算及其结果的分析，得到了一系列与实际地震较类似的人工地震图像，如模型中的地震活动在时间上有轮回怀，在地点上有条带性，并在不同轮回中有条过移的现象。模型中应力变化十分复杂，但仍能为预测地震提供概略性的信息。     

地震复杂性成因的数值模拟

以最简单的弹簧滑块模型为理论基础,由完全确定性的动力学公式出发,采用细胞自动机技术,研究模拟地震序列时、空、强复杂性的成因。当中公有非线性动力学因素时,模拟序列的时空强特征十分复杂,当引入弹性结构的不均匀性,模拟大事件的空间和强度特征明显受控于断层的不均匀性。由模拟结果推测认为地震的复杂性是由于地球的不均匀性和断层破裂的非线性动力学相互作用而引起的。     

基于大陆地震活动特点建立的简化动力学模型--细胞自动机模型

基于大陆地震活动的特点和已建立的由弹簧－滑块－阻尼器组成的平面网络非线性动力学模型研究结果,根据细胞自动机的原理,提出了对这类模型的一种简化模型－细胞自动机模型。其演化规则是根据大陆地震活动特点和从上述非线性模型中总结出来的,介质参数采用分维分布设定。将新建的模型应用于３０×４０单元组成的较大系统中,探讨了该模型显示的一些地震活动特点。     

中国大陆活动地块边界带主要断层的强震震间晚期综合判定

由于中国大陆强震主要分布在活动地块边界带上,所以活动地块边界带主要断层成为我国大陆型强震研究的重要目标,各强震孕育阶段的判定是大陆型强震原地复发的动力学过程研究主要内容,而目标断层是否处于震间晚期也是强震时间预测的重要研究基础.虽然地震短临预测仍存在诸多科学难题,但最近20年来全球若干强震相关研究表明,如果放宽预测时间尺度的要求,有些方法也可用于强震震间晚期的判定.本文以中国大陆活动地块边界带的391条断层段为研究目标,利用地震地质的强震破裂空段、大地测量的断层运动闭锁段、地震活动的中小地震稀疏段、数值模拟的库仑应力增强显著段等方法,综合判定中国大陆活动地块边界带可能处于震间晚期的主要断层段.本文结果仅是初步结果,该结果的可靠程度有赖于监测条件,其科学性有赖于大陆型强震孕育发生动力学过程的认识水平,虽然本文尝试给出中国大陆活动地块边界带主要断层的震间晚期判定结果,但其结果可靠程度、精细程度等均存在巨大的改善空间.最后,从断层孕震阶段判定需求的角度,本文尝试给出大陆型震源物理模型的具体基础模型,期望起到抛砖引玉的作用,也期望更多地震学家关注大陆型强震的物理机制及其预测基础研究.     

几何弯曲断层活动性的模拟

运用有限元法建立大尺度几何弯曲断层（长约750 km）的二维模型,并利用接触单元技术模拟断层间的作用,模拟了几千年时间尺度内走滑断层的活动,探讨了具有一定几何形态断层对断层系统活动的影响.几何弯曲的断层导致了应力的集中,而且在断层的地震事件中起到了抑制作用,但是也为孕育大震提供了条件.和平直断层的模拟结果比较得知,断层的几何弯曲不仅影响地震破裂长度、改变了地震滑移量,也很大程度影响了能量释放的空间位置,改变了地震轮回周期.此外,弯曲断层应力的不均匀也使主震前往往有前震发生.另外通过模拟可以看出,弯曲断层活动表现出明显的特征地震,使得地震-频次曲线偏离G-R公式,说明了对单个断层而言,特征地震可能更能体现断层的地震活动,这对我们进行地震预测和危险性分析有着很重要的参考意义.     

基于龙门山断裂带中段平行逆断层格局 和动力学背景的特征地震数值模拟实验

本文构建一种应用有限元开展特征地震数值模拟的新方法, 并以龙门山断裂带中段的浅层构造和动力学机制为背景, 研究了平行逆冲断层分布格局对区域地震活动性的影响. 结果表明, 从断层活动相互影响的角度看, 包含3条平行逆断层的断裂带的整体地震活动性并不适用严格周期的特征地震模型, 当断层间距在20 km以下时, 随着断层间距的缩短, 对单条断层应用特征地震模型的适用性会逐渐降低. 龙门山断裂带中段的模拟计算结果显示, 后山断裂的地震活动相对独立, 区域活动性和中央断裂的断层活动很可能不适用严格周期的特征地震模型.      

人工地震目录研究进展

基于地震目录的统计分析及由此衍生的地震复发周期与地震预测概率分析一直是地震活动性研究的重要部分,而人类有仪器记录的地震活动历史极其有限,因此采用其他方法建立人工地震目录进而研究地震活动特征与规律具有一定的理论意义.本文系统阐述了弹簧-滑块模型、细胞自动机模型、断层系统地震活动模型、有限元-滑块耦合模型、有限元模型等几种人工地震目录研究的理论模型与方法进展;归纳总结了不同模型在实现地震孕育、发生时不可忽视的关键问题:摩擦本构关系,并详细阐述了摩擦本构关系的进展及其区别;最后对有限元数值模拟人工地震目录研究时存在的问题及需要的进一步改进做了探讨.     

2011年M_W9.0东日本大地震动力学破裂过程的数值模拟

力学上,地震可以看作在应力场作用下由于断层带介质的突然损伤或软化导致的断层带失稳事件.本文基于这个地震动力学模型,利用一种可以模拟断层大位错的有限元方法,研究了2011年MW9.0东日本大地震(Tohoku-Oki)的动力学破裂过程.比较了无障碍体和具有不同刚度障碍体的断层带模型产生的断层位移、位错和应力降.主要结果表明,障碍体的存在并不明显地改变障碍体区域的初始构造应力场.对有障碍体情形,准静态结果显示断层上盘最大逆冲位移和最大剪切位错分别为51m和58m,均发生在海底表面海沟处,与无障碍体的结果(最大剪切位错约55m)相比差别不大;下盘最大倾向位移(-10m)并不与上盘最大值出现在同一位置,而是在障碍体处.障碍体处剪应力降(约11 MPa)大于周围非障碍体区域.障碍体处正应力降的最大值约为3 MPa.模拟结果似乎不支持海山是导致本次地震异乎寻常大位错的原因,而倾向于断层带剪切刚度在地震过程中极度损伤或软化.     

A Computer Model for Cyclic Activity of Strong Earthquakes in Continental Seismic Zones

Based on the study of activity of earthquakes with M**********>7.0 in mainland China,we have found a dynamic pattern,i.e.,the cyclic characteristics in time and migration from one seismic zone to another in space.In order to understand the physical mechanism of this pattern,we use a nonlinear dynamic model to simulate the seismic activity in fault zones under a unified tectonic stress field.The basic elements in our model consist of a Maxwell body and a rigid sliding block.Basic elements in a column represent a fault,and coupling elements connecting different columns simulate the interaction among faults and fault segments.The results provide insights to the cyclic activity of strong earthquakes and to the feature of mutual influence between strong earthquakes in groups in the climax of seismic activity.     

Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region, China

Based on the concept of "active blocks" and spatial distribution of historical earthquakes with surface ruptures as well as major and subordinate active faults. The Sichuan-Yunnan region can be divided into four first-order blocks. They are the Markam block (I), the Sichuan-Yunnan rhombic block (II), Baoshan-Pu'er block (III), and Mizhina-Ximeng block (IV). Cut by sub-ordinate NE-trending active faults, the Sichuan-Yunnan rhombic block (II) can be further divided into two sub-blocks: the northwestern Sichuan sub-block (II1) and the middle Yunnan sub-block (II2), while the Baoshan- Pu'er block (III) can be further divided into three sub-blocks: Baoshan sub-block (III1), Jinggu sub-block (III2), and Mengla sub-block (III3). A quantitative study of offset landforms is carried out and the basic types of active faults and their long-term slip rates along the major boundaries of active blocks of different orders in the Sichuan-Yunnan region are determined, through slip vector analysis, the motion states of the active blocks are clarified and the deformation coordination on the block margins is discussed. It is suggested that the tectonic motion of the blocks in this region is a complex or superimposition of three basic types of motions: southeastward sliding, rotating on vertical axis, and uplifting. The Markam block (I), the northwestern Sichuan sub-block (II1), and middle Yunnan sub-block (II2) have a southeastward horizontal sliding rate of 1-5 mm/a, clockwise rotating angular rate of 1.4-4(/Ma, and uplifting rate of about 1 mm/a. The Baoshan-Pu'er (III) and Mizhina-Ximeng (IV) blocks have also been extensively clockwise rotated. This pattern of motion is a strain response to the collision between the Indian and Eurasian plates and the localized deformation and differential slip on the block margins associated with the northward motion of the Indian Plate. Because a set of transverse thrusts between the blocks absorbs and transforms some components of eastward or southeastward sliding motion, the eastward escape or extrusion of the Tibetan Plateau is limited as "imbricated thrusting transformation-limited extrusion model".     

Linearization of rigid body dynamics on frictional interfaces under harmonic loading

The non-linear dynamic response of the rigid block is linearized by means of a friction model that implicitly considers block response through an experimental parameter obtained from shaking table experiments. In view of the great difficulty in carrying out shaking table experiments, in this technical note some recommendations to estimate friction model parameters are given. The selection of parameters considers the sliding response mode of the block–plane-excitation system: stick–slip or continuous sliding. Once all the friction parameters and block response mode were estimated, a methodology was proposed to compute rigid block dynamic response. The numerical results were then compared to actual experimental data for a rigid block sliding on a geotextile–wood interface, along an inclined plane subjected to base harmonic acceleration. Experiments were carried out for both the stick–slip and the continuous sliding modes. Computed and measured responses for both cases showed good agreement, thus indicating that the methodology developed in this research is adequate to capture the physics (of non-linear nature) of rigid blocks sliding on frictional interfaces subjected to complex harmonic loading. The findings encourage the extension of the linearization technique to the more general seismic loading case.     

Nonplanar Faults: Mechanics of Slip and Off-fault Damage

Stress interactions and sliding characteristics of faults with random fractal waviness in a purely elastic medium differ both qualitatively and quantitatively from those of faults with planar surfaces. With nonplanar fault models, solutions for slip diverge as resolution of the fractal features increases, and the scaling of fault slip with fault rupture dimension becomes nonlinear. We show that the nonlinear scaling of slip and divergence of solutions arise because stresses from geometric interactions at irregularities along nonplanar faults grow with increasing slip and produce backstresses that progressively impede slip. However, in real materials with finite strength, yielding will halt the growth of the interaction stresses, which will profoundly affect slip of nonplanar faults. We infer that in the brittle seismogenic portion of the Earth’s crust, off-fault yielding occurs on pervasive secondary faults. Predicted rates of stress relaxation with distance from major faults with random fractal roughness follow a power-law relationship that is consistent with reported clustering of background seismicity up to 15 kilometers from faults.     

Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region,China

Based on the concept of “active blocks” and spatial distribution of historical earthquakes with surface ruptures as well as major and subordinate active faults. The Sichuan-Yunnan region can be divided into four first-order blocks. They are the Markam block (I), the Sichuan-Yunnan rhombic block (II), Baoshan-Pu’er block (III), and Mizhina-Ximeng block (IV). Cut by sub-ordinate NE-trending active faults, the Sichuan-Yunnan rhombic block (II) can be further divided into two sub-blocks: the northwestern Sichuan sub-block (ll₁) and the middle Yunnan sub-block (ll₂), while the Baoshan- Pu’er block (III) can be further divided into three sub-blocks: Baoshan sub-block (Ill₁), Jinggu sub-block (lll₂), and Mengla sub-block (lll₃). A quantitative study of offset landforms is carried out and the basic types of active faults and their long-term slip rates along the major boundaries of active blocks of different orders in the Sichuan-Yunnan region are determined, through slip vector analysis, the motion states of the active blocks are clarified and the deformation coordination on the block margins is discussed. It is suggested that the tectonic motion of the blocks in this region is a complex or superimposition of three basic types of motions: southeastward sliding, ro-tating on vertical axis, and uplifting. The Markam block (I), the northwestern Sichuan sub-block (ll₁), and middle Yunnan sub-block (ll₂) have a southeastward horizontal sliding rate of 1-5 mm/a, clockwise rotating angular rate of 1.4-4°/Ma, and uplifting rate of about 1 mm/a. The Baoshan-Pu’er (III) and Mizhina-Ximeng (IV) blocks have also been extensively clockwise rotated. This pattern of motion is a strain response to the collision between the Indian and Eurasian plates and the localized deformation and differential slip on the block margins associated with the northward motion of the Indian Plate. Because a set of transverse thrusts between the blocks absorbs and transforms some components of eastward or southeastward sliding motion, the eastward escape or extrusion of the Tibetan Plateau is limited as “imbricated thrusting transformation-limited extrusion model”.     

Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region,China

Based on the concept of “active blocks” and spatial distribution of historical earthquakes with surface ruptures as well as major and subordinate active faults. The Sichuan-Yunnan region can be divided into four first-order blocks. They are the Markam block (I), the Sichuan-Yunnan rhombic block (II), Baoshan-Pu’er block (III), and Mizhina-Ximeng block (IV). Cut by sub-ordinate NE-trending active faults, the Sichuan-Yunnan rhombic block (II) can be further divided into two sub-blocks: the northwestern Sichuan sub-block (ll₁) and the middle Yunnan sub-block (ll₂), while the Baoshan- Pu’er block (III) can be further divided into three sub-blocks: Baoshan sub-block (Ill₁), Jinggu sub-block (lll₂), and Mengla sub-block (lll₃). A quantitative study of offset landforms is carried out and the basic types of active faults and their long-term slip rates along the major boundaries of active blocks of different orders in the Sichuan-Yunnan region are determined, through slip vector analysis, the motion states of the active blocks are clarified and the deformation coordination on the block margins is discussed. It is suggested that the tectonic motion of the blocks in this region is a complex or superimposition of three basic types of motions: southeastward sliding, ro-tating on vertical axis, and uplifting. The Markam block (I), the northwestern Sichuan sub-block (ll₁), and middle Yunnan sub-block (ll₂) have a southeastward horizontal sliding rate of 1-5 mm/a, clockwise rotating angular rate of 1.4-4°/Ma, and uplifting rate of about 1 mm/a. The Baoshan-Pu’er (III) and Mizhina-Ximeng (IV) blocks have also been extensively clockwise rotated. This pattern of motion is a strain response to the collision between the Indian and Eurasian plates and the localized deformation and differential slip on the block margins associated with the northward motion of the Indian Plate. Because a set of transverse thrusts between the blocks absorbs and transforms some components of eastward or southeastward sliding motion, the eastward escape or extrusion of the Tibetan Plateau is limited as “imbricated thrusting transformation-limited extrusion model”.

     

Dynamic analysis of stacked rigid blocks

The dynamic behavior of a structural model of two stacked rigid blocks subjected to ground excitation is examined. Assuming no sliding, the rocking response of the system standing free on a rigid foundation is investigated. The derivation of the equations of motion accounts for the consecutive transition from one pattern of motion to another, each being governed by a set of highly nonlinear differential equations. The system behavior is described in terms of four possible patterns of response and impact between either the two blocks or the base block and the ground. The equations governing the rocking response of the system to horizontal and vertical ground accelerations are derived for each pattern, and an impact model is developed by conservation of angular momentum considerations. Numerical results are obtained by developing an ad hoc computational scheme that is capable of determining the response of the system under an arbitrary base excitation. This feature is demonstrated by using accelerograms from the Northridge, CA, 1994, earthquake. It is hoped that the two-blocks model used herein can facilitate the development of more sophisticated multi-block structural models.     

Seismic assessment of the rigid sliding displacements caused by pulse motions

This paper presents a comprehensive study on the rigid block sliding displacement of slopes subjected to ground motions with large velocity pulses. A comparison of the performance of various existing empirical displacement models is provided through analyses of the displacement residuals of slopes subject to pulse-like motions. Except for the PGA- and PGV-based Saygili and Rathje model (2008, referred to as SR08), positive medians of residuals are observed for selected models, indicating an under-estimation. There is a negative constant shift in the total residuals for the SR08 model, which can be easily fixed by changing the constant term in the predictive equation. The residuals from the SR08 model also have the smallest standard deviation compared to the other models. A modified SR08 model is developed for predicting rigid block sliding displacements for pulse-like motions. The modified predictive model is used in probabilistic seismic displacement analyses of slopes in a hypothetical near-fault region.     

抗震销的模拟及其对异型板桥抗震性能影响分析

抗震销在上海莘庄立交工程的设计中,由于构造的特殊性,其力学特性难以用常规方法表达。基于位移和多级设防的抗震思想,以莘庄立交工程中的一异型板桥为例,首次以两个包含极端边界条件的结构分析模型考虑抗震销这一特殊构造措施对结构的影响,并按照实际结构构造充分考虑材料非线性及连接构件的非线性等因素,分别以弹塑性纤维梁柱单元、空间滑动支座单元和接触单元模拟桥墩屈服、支座滑动、相邻结构之间碰撞,建立了空间非线性时程反应分析模型,并进行了设计和罕遇地震水平两个概率水准下的地震反应分析和抗震性能评估,并根据分析结果的对比,提出强化抗震销刚度的建议。