震源动力学破裂过程数值模拟研究

首先阐述了震源动力学过程研究的重要意义, 在此基础上, 研究了复杂的断层几何形态及介质模型对动力学破裂过程的影响, 并对常用的有限元方法、 离散元方法、 有限差分方法和边界积分方程方法等进行了相应介绍. 讨论了这些数值模拟方法各自的优缺点, 建议在方法的选择上应视具体问题及计算的精度而定. 最后对动力学数值模拟的关键部分, 滑动摩擦准则进行了论述. 常用的滑动摩擦准则有滑动弱化准则、 速率弱化准则和速率-状态依赖摩擦准则. 在单纯考虑某个地震的动力学破裂传播过程时, 滑动弱化准则较为常用, 其中滑动弱化距离的选取至关重要. 但若考虑整个地震循环, 速率-状态依赖摩擦准则更为合适.      

地震结构爆裂理论与短临预测

在地震结构爆裂动力学理论基础上,揭示了临震前兆波的客观机制。作为储能体的地震包体,其流体压力达到覆盖岩层的断裂强度和重力之和时,在引潮力的作用下,包体发生局部破裂,裂纹发生扩张,产生破裂流和破裂波,交替发生冲击加载和垮塌卸载,同潮汐波叠加,形成地震前兆共振波,以相速度和群速度向围岩传播,其振动频率和周期表征包体的固有频率和周期。波速和周期均可以被灵敏度高和抗干扰性能强的PS100HRT地电仪在地震爆发前捕捉到,从而可用于确定地震的震中和震级。依据共振波出现后振幅随时间增加的规律,月潮规则,可以确定发震时间。地震和临震理论以及诊断技术均已经成熟,中国有望在世界上首先实现地震预报。查清地震包体,地震能可以如同油气那样供人类利用。     

地震破裂模式的模糊模式识别

把几种可能的物理机制作为识别对象,定义了计算各识别对象与物理现象相符合的证据量．从这种思想出发,提出了判定强震破裂模式的模糊模式识别准则,只要把各识别对象的“证据量”带人这个识别准则,就可识别出某指定强震的地震破裂机制该方法被具体地应用于炉霍地震（1973,Ms＝7.9）破裂机制的识别,结果表明炉霍地震符合barrier破裂机制这个结论与有关炉霍地震破裂机制的地震学研究结果是完全一致的该方法比已有的其他方法更为简便,因此有推广价值.     

Discrete modelling of time‐dependent rockfill behaviour

This modelling study deals with the time‐dependent behaviour of rockfill media, which is of particular interest during the life of rockfill dams. Breakage of rock blocks and crack propagation are the main processes responsible for rockfill creep and collapse. The modelling procedure presented here is performed on two scales: on the rock block scale, where the grain is taken to be an assembly of rigid particles initially endowed with cohesive bonds, and on the rockfill scale, which is taken to involve a set of breakable grains interacting via contact and friction processes. The grain breakage process is described in term of a thermodynamically consistent damage interface model, where the damage is a gradual delayed process. This model was implemented in a non‐smooth contact dynamics code. The effects of the main parameters involved were analysed by performing numerical studies. The ability of the model to predict the creep behaviour of rockfill media is confirmed by presenting several simulations. Copyright © 2008 John Wiley & Sons, Ltd.     

Dynamic nucleation process of shallow earthquake faulting in a fault zone

Two distinct phases are commonly observed at the initial part of seismograms of large shallow earthquakes: low-frequency and low-amplitude waves following the onset of a P wave ( P ₁) are interrupted by the arrival of the second impulsive phase P₂ enriched with high-frequency components. This observation suggests that a large shallow earthquake involves two qualitatively different stages of rupture at its nucleation.
We propose a theoretical model that can naturally explain the above nucleation behaviour. The model is 2-D and the deformation is assumed to be anti-plane. A key clement in our model is the assumption of a zone in which numbers of pre-existing cracks are densely distributed; this cracked zone is a model for the fault zone. Dynamic crack growth nucleated in such a zone is intensely affected by the crack interactions, which exert two conflicting effects: one tends to accelerate the crack growth, and the other tends to decelerate it. The accelerating and decelerating effects are generally ascribable to coplanar and non-coplanar crack interactions, respectively. We rigorously treat the multiple interactions among the cracks, using the boundary integral equation method (BIEM), and assume the critical stress fracture criterion for the analysis of spontaneous crack propagation.
Our analysis shows that a dynamic rupture nucleated in the cracked zone begins to grow slowly due to the relative predominance of non-coplanar interactions. This process radiates the P₁ phase. If the crack continues to grow, coalescence with adjacent coplanar cracks occurs after a short time. Then, coplanar interactions suddenly begin to prevail and crack growth is accelerated; the P₂ phase is emitted in this process. It is interpreted that the two distinct phases appear in the process of the transition from non-coplanar to coplanar interaction predominance.     

压剪应力场中组合裂纹非共面破裂特征的实验研究

通过实验力学方法研究了几种典型组合裂纹非共面破裂特征．实验结果表明，在双轴加载条件下，脆性材料内部的初始不闭合裂纹在受力闭合之前出现新裂纹，扩展裂纹中不仅有张性裂纹，而且有剪切裂纹，部分甚至还出现逆向剪切裂纹．对单轴加载和双轴加载两种条件下裂纹的破坏规律作了比较．     

Simulation of the spontaneous growth of a dynamic crack without constraints on the crack tip path

The spontaneous growth of a dynamic in-plane shear crack is simulated using a newly developed method of analysis in which no a priori constraint is required for the crack tip path, unlike in other classical studies. We formulate the problem in terms of boundary integral equations; the hypersingularities of the integration kernels are removed by taking the finite parts. Our analysis shows that dynamic crack growth is spontaneously arrested soon after the bending of the crack tips, even in a uniformly stressed medium with homogeneously distributed fracture strengths. This shows that the dynamics of crack growth has a significant effect on forming the non-planar crack shape, and consequently plays an essential role in the arrest of earthquake rupturing.     

Rupture model of the 2013 M_W 6.6 Lushan (China) earthquake constrained by a new GPS data set and its effects on potential seismic hazard

Vertical records are critically important when determining the rupture model of an earthquake, especially a thrust earthquake. Due to the relatively low fitness level of near-field vertical displacements, the precision of previous rupture models is relatively low, and the seismic hazard evaluated thereafter should be further updated. In this study, we applied three-component displacement records from GPS stations in and around the source region of the 2013 MW6.6 Lushan earthquake to re-investigate the rupture model.To improve the resolution of the rupture model, records from both continuous and campaign GPS stations were gathered, and secular deformations of the GPS movements were removed from the records of the campaign stations to ensure their reliability. The rupture model was derived by the steepest descent method(SDM), which is based on a layered velocity structure. The peak slip value was about 0.75 m, with a seismic moment release of 9.89 × 10~(18) N·m, which was equivalent to an M_W6.6 event. The inferred fault geometry coincided well with the aftershock distribution of the Lushan earthquake. Unlike previous rupture models, a secondary slip asperity existed at a shallow depth and even touched the ground surface. Based on the distribution of the co-seismic ruptures of the Lushan and Wenchuan earthquakes, post-seismic relaxation of the Wenchuan earthquake, and tectonic loading process, we proposed that the seismic hazard is quite high and still needs special attention in the seismic gap between the two earthquakes.     

2004年苏门答腊大地震激发的T-波

T-波是由海底地震或者海陆边界俯冲带附近地震激发,并在海洋低速层中传播的声波.2004年12月26日,在印度洋东部印尼苏门答腊岛附近发生MW=9.3级大地震,其产生的能量在印度洋中激发了巨大的海啸,造成了严重的人员伤亡和财产损失,受到了世界科学家们极大的关注.本文从台站(PALK)及台站(DGAR)记录到的地震的信号中,提取出了清晰的高频T-波,并在频率域内分析,最终得到了T-波的频谱已及频率随时间变化图像.另外,通过对大地震时间域和频率域内T-波信号的分析,了解到此次大地震断层破裂过程持续的时间大致为500 s,其间伴随有两次明显的能量释放过程.分析数据表明两次能量释放过程的间隔大致为80~100 s.T-波分析将为推断海洋地震以及海陆边界俯冲带附近地震的特征,提供一种独立的研究手段和方法.     

地震破裂过程的几何学与运动学特征的模拟

本应用形变破裂和激光全息光弹实验，结合断裂力学的观点研究了地震破裂过程的三个方面：１．地震破裂的力学机制；２．地震破裂的应变特征与运动过程；３．地震破裂过程的应力场分布特征。并分析了地震破裂过程的几何学与运动学特征以及探讨了它们与地震发生。前兆和余震迁移的密切关系。此项研究为活动断层分段研究和地震预报提供了实验证据。