Elastic wave propagation in an irregularly layered medium

The indirect boundary element method (IBEM) is used to simulate wave propagation in two-dimensional irregularly layered elastic media for internal line sources. The method is based on the integral representation for scattered elastic waves using single layer boundary sources. Fulfillment of the boundary conditions leads to a system of integral equations. Results are obtained in the frequency domain and seismograins are computed through Fourier synthesis. In order to test and validate the method we present various comparisons between our results and the time series obtained analytically for a buried line source in a half-space and by using the recently developed spectral element method (SEM).     

高铁震源地震信号的稀疏化建模

我国每天有数千趟高铁列车运行在两万多公里的高铁线路上,不但会引起高铁路基的振动,还会激发出地震波.地震检波器所接收到的数据中不仅包含窄带分立谱特性的高铁震源地震信号,还包含宽频带特性的背景信号.如何实现从检波器所接收到的高铁震源地震数据中分离出高铁震源地震信号和宽频带背景信号是准确利用该类信号的关键.考虑到高铁震源地震信号与宽频带信号在频率域明显的形态特征差异,本文首次将形态成分分析这种信号分离手段引入到高铁震源地震信号处理中,实现高铁震源地震信号的稀疏化建模并进而实现从接收数据中分离出高铁震源地震信号以及宽频带背景信号.对北京大学在中国南方某高铁沿线采集到的大量高铁震源地震数据进行处理,结果表明:采用形态成分分析并结合分块坐标松弛算法,能够实现实际采集高铁震源地震数据中的高铁震源地震信号和宽频带信号的分离.     

基于黏弹性边界的P波斜入射波场分解方法精度比较

黏弹性人工边界以其概念清晰,稳定性高,精度良好的优点,在地下结构抗震的研究中已经得到了较为广泛的应用。在使用黏弹性边界进行地震动的输入时,需要把地震波转化为模型网格节点上的等效节点力,这个过程就要涉及到波场分解的问题。波场分解的方式可以不同,但不同的波场分解方式会得出精度不同的计算结果,而学界目前尚未就此达成一致。本文在前人研究的基础上,总结了3种地震波斜入射的波场分解方法,编制了相应的Python脚本来实现黏弹性边界和等效节点力的批量施加,并且以不同角度的地震波斜入射为例在同一有限元模型中计算以进行比较研究。结果表明计算精度与波场分解的方式密切相关:在模型边界上考虑的自由场波动叠加越完整,计算结果越精确。     

城市地下变电站非线性地震反应分析

为研究城市地下变电站结构及其内部电气设备的动力特性及地震反应规律,以西安市快速轨道交通二号线张家堡地下变电站为工程背景,依据黏弹性人工边界理论及有限单元法,建立考虑土体-结构-设备动力相互作用的三维有限元计算模型,首先通过半空间自由场模型验证了地震动输入方法的有效性,继而开展地下变电站系统的模态分析和动力时程分析,得到了该类结构的动力特性及地震响应规律。研究结果表明,结构上覆土体及电气设备在地震作用下产生的惯性力是结构构件应力增长的主要原因,8度罕遇地震作用下地下变电站内电气设备的加速度放大系数明显超出规范要求,需要采取相应保护措施。     

DEFORMATION OF THE BRITTLE-PLASTIC TRANSITION ZONE AT THE POST-SEISMIC RELAXATION PERIOD: A CASE STUDY OF THE RED RIVER FAULT

The transition from microscopic brittle deformation to microscopic plastic deformation is called brittle-plastic transition, which is considered as a key layer for determining the limit of lower continental crust seismicity. The depth and deformation mechanism of the brittle-plastic transition zone is controlled mainly by temperature. Besides, the strain rate and fluid pore pressure also affect the transition during the different deformation stages at the seismic cycle. In this paper, microstructure observation of catalcastic samples collected from the Red River Fault was carried out using optical polarized microscopy and scanning electron microscopy. The morphology, microstructures of deformation characteristics, mineral composition, water-rock reaction, pressure solution, exsolution, crack healing in the samples were systematically observed. The mineral components quantitative analyses were examined using the EDS. Water-rock reaction and pressure solution were systematically observed under SEM. The fabric of the main minerals in the samples was measured using electron backscattered diffraction(EBSD). Based on these analyses, the deformation mode was setup for the brittle-plastic transition zone of the fault during the post-seismic relaxation period. Both brittle deformation and plastic deformation were developed in the cataclastic samples. EBSD data shows that the c axial fabrics of quartz present low-temperature plastic deformation characteristics. The feldspar deformed as cataclastic rock, and the micro-fracture in feldspar was healed by static recrystallized quartz and calcite veins. The calcite vein underwent plastic deformation, which represents the post-seismic relaxation deformation. Based on the analysis of deformation mechanism of cataclastic samples in brittle-plastic transition zone of the Red River Fault, and combined with previous studies, we concluded that the brittle fracture and fracture healing is the main deformation mode at brittle-plastic transition zone in the post-seismic relaxation. High stress and high strain rate at post-seismic relaxation lead to brittle fracture of high-strength minerals such as feldspar in rocks. Plastic deformation occurs in low-strength minerals such as quartz and mica. Under the fluid condition, micro-fractures were healed by quartz and calcite. The minerals such as quartz and calcite in the fracture transformed from static recrystallization to dynamic recrystallization with stress gradually accumulating. With fracture healing and stress accumulation, the fault strength gradually increases which could accumulate energy for the next earthquake.     

2004年苏门答腊地震粘滞性松弛效应对华南地区地壳水平活动的影响

2004年苏门答腊地震后, 不同学者根据不同观测数据(地震波、 GPS), 得到了此次地震的断层滑动模型。 反演过程中使用半无限空间模型时, 无法利用远场观测数据进行约束, 势必影响远场形变的解释。 基于Hoechner等使用的断层几何模型和GPS同震位移数据, 本研究利用球体位错理论反演方法反演了2004年苏门答腊地震断层滑动模型, 得到的矩震级为9.24, 最大滑移量为30.4 m, 由于考虑了曲率的效应, 该模型在远场同震位移的计算结果与GPS数据吻合较好。 然后, 选取了2001—2004年和2004—2007年两期的GPS水平位移速度场, 研究2004年苏门答腊地震对华南地区地壳水平活动的影响, 从两期的GPS水平位移速度场差异可以看出地震后华南块体有向西南方向的运动趋势, 华南块体受到此次地震明显的震后影响。 最后, 基于反演得到的断层模型, 利用Tanaka等提出的粘弹性球体位错理论对华南块体两期GPS水平位移速度场差异进行模拟, 得到华南块体内部粘滞性系数为2×10¹⁹ Pa·s, 当考虑地幔粘滞性松弛效应后, 两期的速度场差异的均方根值由3.2 mm减少为1.9 mm。 可见在研究2004年前后中国大陆GPS水平位移速度场时, 若继续以华南块体为基准, 需考虑此次地震的地幔粘滞性松弛效应。     

A prediction model for frequency spectrum of blast‐induced seismic wave in viscoelastic medium

A prediction model for frequency spectrum of blast‐induced seismic waves is established. The effect of explosive sources is considered in this model. Our model implies that the frequency spectrum of blast‐induced seismic wave is mainly influenced by the initial pressure and the adiabatic exponent of explosives. The dominant frequency increases with the decreasing of initial pressure or the increasing of adiabatic exponent. In addition, this prediction model is verified by the experiment. The error of the dominant frequency is 4%–6%. It is indicated that the proposed model in this paper can reasonably predict the frequency spectrum of blast‐induced seismic waves, and then, we can provide a better frequency spectrum by optimizing the explosion source.     

2011年日本MW9.0地震震后形变机制与震源区总体构造特征

本文利用较为完备的球体位错理论,结合4.5年的震后位移数据,优化了2011年日本M_W9.0地震震源区岩石圈弹性层厚度与地幔黏滞性因子,更新了该强震断层余滑时空演化过程.首先,基于日本列岛215个均匀分布的GPS连续观测站震前2年与震后4.5年的观测数据,提取了2011年日本M_W9.0地震引起的震后位移时空变化;接着,依据断层余滑衰减相对较快的特点,利用黏弹性球体位错理论对震后3~4.5年的GPS观测数据进行反复拟合,确定2011年日本M_W9.0地震震源区地幔黏滞性系数和岩石圈弹性层厚度的最优解分别为6×10¹⁸ Pa·s和30 km;然后,从震后3年内GPS观测数据中剔除地幔黏滞性松弛效应,获取了断层余滑对应的震后位移场;最后,利用基于球体位错理论的反演算法,反演了2011年日本M_W9.0地震断层余滑的时空演化过程.结果表明,2011年日本M_W9.0地震引起的断层余滑在震后半年内变化显著,震后2年主震区域余滑基本停止,断层的两端存在一定的余滑效应,断层余滑的累计矩震级达到8.59;地震后4年,地幔黏滞性松弛效应对震后位移场的贡献在总体上超过断层余滑的贡献.     

Seismic performance evaluation of typical dampers designed by Chinese building code

Adding dampers is a commonly adopted seismic risk mitigation strategy for modern buildings, and the corresponding design procedure of dampers has been well established by the Chinese Building Code. Even though all types of dampers are designed by the same procedure, actual seismic performance of the building may differ from one to the others. In this study, a nine-story benchmark steel building is established, and three different and typical types of dampers are designed according to the Chinese Building Code to realize structural vibration control under strong earthquake excitation. The seismic response of the prototype building equipped with a viscoelastic damper, viscous damper and buckling-restrained brace(BRB) subjected to 10 earthquake records are calculated, and Incremental Dynamic Analysis(IDA) is performed to describe progressive damage of the structure under increasing earthquake intensity. In the perspective of fragility, it shows that the viscoelastic damper has the highest collapse margin ratio(CMR), and the viscous damper provides the best drift control. Both the BRB and viscoelastic dampers can effectively reduce the floor acceleration responses in the mid-rise building.