拐折断层黏滑过程的实验研究

采用位移、应变和声发射等测量手段,研究了预置5°拐折断层的房山花岗岩样品的黏滑过程,分析了不同加载速率下5°拐折断层失稳的黏滑特征及相关物理场的演化过程。实验结果表明: 1)5°拐折断层的黏滑周期与加载速率在数值上呈负对数相关关系; 2)在不同的加载速率下,大多数的5°拐折断层失稳是双震事件, 2次子事件的间隔时间大多在100～200ms之间; 3)采用不同的观测手段,即使采样速率一致,其临震响应也存在差异性,如断层失稳前沿断层的应变测量结果呈现明显的应变弱化,断层位移则未见明显的变化; 4)断层黏滑过程中的声发射事件呈现明显的沿断层迁移的特征。认识强震的发生机理和余震特征需要进一步研究断层失稳过程的动力学信息。     

5°拐折断层在黏滑过程中物理场演化与交替活动的实验研究

在实验室利用96通道应变记录采集系统和分布式多通道瞬态信号采集系统,观测了预切5.拐折断层的标本在变形失稳过程中应变场和声发射事件的时空演化.实验在双轴伺服加载系统上进行.在Y方向按位移控制方式加载,位移速率先后取0.5μm/s、1μm/s、0.5μm/s和0.1μm/s.观测得到:1)标本沿断层发生周期性的黏滑失稳,...     

Spatio-temporal characteristics of acoustic emission during the deformation of rock samples with compressional and extensional en-echelon faults

Introduction By the experimental techniques of acoustic emission (AE) and surface strain observation, a lot of achievements on the study of deformation and fracture process of samples with en-echelon faults has been obtained (Bomblakis, 1973; LIU, et al, 1986; DU, et al, 1989; LI, et al, 1989; ZHANG, FANG, 1989; ZHANG, LI, 1989; MA, et al, 1995a, b; LI, et al, 2000). However, there still exist many problems. For example, most of discussion on the spatio-temporal evolution of microf…     

雁列式断层组合变形过程中的声发射活动特征

着重研究拉张或挤压型雁列式断层标本变形过程中声发射(简称AE，下同)的时空演化特征. 结果表明，预置构造对AE空间分布格局具有较强的控制作用，随着差应力的增加，AE首先在两个端点附近丛集，之后向两端点连线附近集中，出现明显的破裂局部化现象，较大事件还通常集中于某一端点附近反复发生. 前期微破裂丛集图象指示后期宏观破裂的扩展方向及扩展范围. 拉张和挤压型雁列区宏观破裂方向分别与轴向应力方向垂直和平行. 雁列式断层标本变形过程中，破裂前的弱化阶段相对较长、弱化过程明显. 微破裂事件累积频次指数增长可能是系统失稳前的典型征兆之一，而雁列区宏观破裂之后，AE数量逐渐减少、应变释放相对减弱. 摩擦滑动过程中，大的粘滑失稳前未见有AE活动前兆性的增强过程.雁列区的b值变化在失稳前显示趋势性降低——快速回升这一典型的变化过程，b值降低一般发生在差应力增强过程之中，并有可能延续至弱化阶段，而快速回升则一般发生在破裂失稳前的短时期内. 对比研究表明，构造差异所导致的b值差异远大于b值随差应力的增加而产生的变化量，而对同一构造标本，力学状态的改变会导致AE序列时序特征的急剧变化，较高的加载速率对应较高的应变能释放及明显的低b值. 先期破碎带由于较低的破坏强度，其对差应力的微小变化具有特殊敏感性，从而成为源、兆分离、窗口或敏感点效应等地震活动性前兆现象的一种可能的原因.      

Analysis of slip-softening instability on a fault with asperities

Slip-softening instability on a vertical strike-slip fault with asperities has been analysed. The fault strength is uniform in depth, but the strength is nonuniform in the strike direction, i.e., there are asperities on the fault. These asperities and other segments of the fault have the same type of constitutive law but different peak stresses. The material surrounding the fault is represented by elastic plates, of which the top and bottom surfaces are stress-free.We use a finite element method to study the evolution of theoretical displacement, stress and strain field with a growing displacement applied at the remote plate ends. The slip and frictional stress are obtained as part of the solution. We have compared the difference of theoretical displacement, strain field and the distribution of frictional stress on the fault between unstable and stable slip. In addition, we have studied the effect of size and strength of asperities on instability, and the softening behaviour of asperities before instability.We find that (1) the failure of the fault zone may be due to either dynamic instability or rapid quasistable slip. A general characteristic of unstable mode is that slippage, on some parts of asperities increases indefinitely for a small finite increase in remote imposed displacement until, immediately before the unstable slip; (2) the size and peak strength of asperities have a large effect on instability. Reducing the size and peak strength of asperities tends to replace inertially unstable deformation with stable deformation; (3) the location with maximum acceleration during unstable slip, as the plausible nucleating seismic source, is in asperities; (4) the shapes of the changes in theoretical stress and strain at a given location, caused by the nonlinear constitutive property of the fault, are all similar whether instability, happens or not. This fact suggests that the changes of peak type or bend type in crustal deformation are not required for earthquake instability.     

Granite deformation and behavior of acoustic emission sequence under the temperature and pressure condition at different crust depths

IntroductionFocal depths of shallow strong intraplate earthquakes are mostly distributed in the highstrength range of lithosphere controlled by the rheology property of granite and diorite (Brace,Kohlstedt, 1980; Sibson, 1982; Meissner, Strehlau, 1982; CHEN, Molnap 1983). This is obviouslyrelated to the change of rock deformation characteriStics at different crust depths. So, for earthquake stUdy, both of the rock failure types (fractUre or rock flow) and its mechanical instabilityforms (…     

Evolution of Wear and Friction Along Experimental Faults

We investigate the evolution of wear and friction along experimental faults composed of solid rock blocks. This evolution is analyzed through shear experiments along five rock types, and the experiments were conducted in a rotary apparatus at slip velocities of 0.002–0.97 m/s, slip distances from a few millimeters to tens of meters, and normal stress of 0.25–6.9 MPa. The wear and friction measurements and fault surface observations revealed three evolution phases: A) An initial stage (slip distances <50 mm) of wear by failure of isolated asperities associated with roughening of the fault surface; B) a running-in stage of slip distances of 1–3 m with intense wear-rate, failure of many asperities, and simultaneous reduction of the friction coefficient and wear-rate; and C) a steady-state stage that initiates when the fault surface is covered by a gouge layer, and during which both wear-rate and friction coefficient maintain quasi-constant, low levels. While these evolution stages are clearly recognizable for experimental faults made from bare rock blocks, our analysis suggests that natural faults “bypass” the first two stages and slip at gouge-controlled steady-state conditions.     

不同温压条件下声发射应变能释放特征——加速模型参数物理含义的初步讨论

利用不同温、压条件下的花岗岩变形实验数据,研究声发射（AE）事件应变释放特征,探讨加速模型参数m值与温压环境的关系.常温条件下,声发射应变显示一定的加速释放特征,但m值随围压增加未显示出趋势性的变化,表明常温条件下m值与岩石强度关系不密切.围压固定时,m值随温度升高逐渐变大,声发射应变从加速释放逐渐过渡到匀速释放,这意味着不同温度条件下岩石变形过程中内部微破裂形式的差异,可能导致应变释放类型的较大差异（即m值的较大差异）.在浅表地层的温压条件下,岩石破坏前显示一定的加速释放特征,m小于1;在渐进式破坏区段,应变释放呈逐渐减弱的减速释放态势,m明显大于1;在深部温压条件下,应变释放加速特征明显,m值明显较低.此外,完整岩样破裂前声发射应变加速释放特征显著,而宏观剪裂面的黏滑之前,声发射应变基本上匀速释放.     

利用模板匹配技术检测米尺度岩石断层黏滑实验中的声发射事件

声发射观测是岩石摩擦实验中研究强震孕育机制的重要观测手段之一.传统的声发射观测以触发记录为主,难以捕捉并提取微弱的声发射信号.近年来天然地震学中发展的模板匹配技术可以在连续波形记录中识别出微弱的事件,大大提升事件检测能力.本文发展了适用于多通道连续记录的声发射信号的模板匹配技术,并将该技术应用于分析米尺度岩石断层黏滑失稳过程.研究结果显示:模板匹配技术识别出的声发射数量约为传统方法识别的5倍,相对完备震级降低约0.3.临近失稳,检测目录完备震级以上声发射率表现出幂律增加的特征,且声发射事件聚集在失稳破裂起始点附近.失稳后,检测目录完备震级以上声发射率表现出幂律衰减,其中失稳后早期表现出较慢的衰减速度.另外,失稳后早期声发射事件表现出随对数时间沿断层迁移的趋势.最后,通过与天然地震学观测的结果对比表明,声发射模板匹配识别技术是研究大尺度岩石摩擦滑动失稳过程中前震和余震时空演化规律的有效手段,该技术助于在方法上为室内实验和野外观测进行融合研究提供思路.     

Spatial and temporal characteristics of acoustic emission activity during deformation of rock samples with inhomogeneous fault under biaxial compression

The spatial and temporal distributions of acoustic emission （AE） during the deformation of samples containing an inhomogeneous fault have been studied under biaxial compression. The results show that the fault strength and the duration from loading to failure increase and the failure mode changes from abrupt instability to gradual failure with increase of lateral stress σ2. The pre-setting fault and its heterogeneity play an important role in controlling AE spatial distribution during the deformation. The basic pattern of AE spatial distribution is controlled by the pre-setting fault, especially the parts with inhomogeneous strength and the strong segment, and the localization of fracturing starts from the positions with inhomogeneous strength, With increase of σ5, the dense AE distribution area spreads from the positions with inhomogeneous strength to the whole strong segment gradually. AE temporal sequence is significantly affected by σ2. The fault shows abrupt instability in final failure, which occurs on a background of ＂enhancement-quiescence＂ in microfracturing activity at lower σ2. At higher σ5, the fault shows gradual failure, and AE occurrence rate increases continuously and AER increase exponentially before and after the failure, The effect of σ2 on b-value is also remarkable, b-value shows precursory decrease in the weakening stage when fault behavior is abrupt instability, but it shows balanced change when fault behavior is gradual failure,     

含障碍体平直断层标本变形过程中群体微破裂事件的时空演化特征

设计了含圆柱形障碍体的平直断层物理模型，在概念上可模 拟断层面上与地震孕育相关的障碍体等强固区域. 开展了中等尺度岩石标本的双轴压缩实验 ，研究标本在变形过程中群体微破裂事件的时空演化. 结果表明：(1)由于标本预置断层两 盘与障碍体之间、以及断层两盘彼此之间复杂的相互作用，使得交替活动成为含障碍体平直 断层标本变形过程中最为显著的特征； (2)含障碍体平直断层标本能够发生动力学失稳的前 提条件是障碍体在外加载荷作用下首先发生破坏，微破裂群体在障碍体破坏前后显示不同的 时空演化特征；(3)含障碍体平直断层标本失稳前的总变形量较大，弹性变形阶段微破裂累 积频次呈指数增长. 而在其后的位移弱化阶段，在障碍体及围绕障碍体的围岩区域内微破裂 事件稀少，突发失稳发生在相对平静的背景之上；(4)不同构造组合变形过程中的微破裂群 体具有明显不同的b值，b值与构造相关的区域性差异是明显的.     

汶川地震和芦山地震的孕震机理及震前中长期地震危险性研究

运用非连续变形分析法与三维有限元法相结合的方法，以GPS资料作为位移速率和震源机制的约束条件，通过数值模拟研究了青藏高原及其东侧邻区构造地块的运动、变形、相互作用及其与近30年来发生于该区的大地震之间的关系。研究中引入了以应力与摩擦强度的比值定义的断层“失稳危险度”，通过数值模拟计算得到了研究区地壳块体边界断层的失稳危险度分布。结果表明，失稳危险度高的地段与近期该区发生的M_S≥7.0地震所在的位置基本一致，其中龙门山断裂带上包括汶川和芦山大地震的发震断层均为失稳危险度最高值地区。计算得到的应变率强度分布图显示，青藏高原东部边缘整条地带均为应变率强度的陡变带，特别是以龙门山断裂带上的陡变最为明显，其西侧应变率强度为东侧的近4倍，而且，这个带位于宽度相同、走向与龙门山断裂带走向相一致的高应变能密度带中，表明这两次大地震前，作为其发震断层的龙门山断裂带已积累了相当高的应变能，失稳危险度高，处于力学上的不稳定状态。模拟计算得到在上地壳层中，2001年昆仑山口西M_S8.1地震引起汶川、芦山地震发震断层的库仑破裂应力增加约0.016 MPa，相当于龙门山断裂带约两年的应力积累，也就是说，使汶川、芦山地震发震断层的失稳破裂提前了约两年。 此外，关于2008年汶川M_S8.0地震的模拟计算表明，汶川地震的发生也使包括芦山地震发震断层的龙门山断裂带西南段和东昆仑断裂带东南端的库仑破裂应力增大，应变能积累增强，这说明汶川M_S8.0地震的发生对已处于失稳危险度较高状态的2013年芦山地震和2017年九寨沟地震发震断层的提前失稳破裂起到了促进作用。      

A Damage Model of Acoustic Emission in Rock and Its Initial Application in Seismological Research

The model of the relationship of AE parameter,damage variable and strain is derived by applying the damage theory and micro-element statistical strength theory.The relation between AE characteristics during rock failure and machine stiffness is analyzed under uniaxial compression with the above model.Based on the above analysis,the internal connection among AE activity law and seisraogenic process and earthquake activity is discussed.     

岩石声发射的损伤模式及其在地震研究中的初步应用

本文运用损伤理论和微元强度统计理论导出了声发射参数与损伤参量及应变的关系模型。依据这一模型,分析了单轴压缩时岩石破坏阶段声发射特性与压机刚度的关系。在此基础上,探讨了声发射活动规律与孕震过程及地震活动的内在本质联系。     

Characterization of Stress Drops on Asperities Estimated from the Heterogeneous Kinematic Slip Model for Strong Motion Prediction for Inland Crustal Earthquakes in Japan

Dense strong motion observation networks provided us with valuable data for studying strong motion generation from large earthquakes. From kinematic waveform inversion of seismic data, the slip distribution on the fault surface of large earthquakes is known to be spatially heterogeneous. Because heterogeneities in the slip and stress drop distributions control the generation of near-source ground motion, it is important to characterize these heterogeneities for past earthquakes in constructing a source model for reliable prediction of strong ground motion. The stress changes during large earthquakes on the faults recently occurring in Japan are estimated from the detailed slip models obtained by the kinematic waveform inversion. The stress drops on and off asperities are summarized on the basis of the stress change distributions obtained here. In this paper, we define the asperity to be a rectangular area whose slip is 1.5 or more times larger than the average slip over the fault according to the previous study for inland crustal earthquakes. The average static stress drops on the asperities of the earthquakes studied here are in the range 6?C23?MPa, whereas those off the asperities are below 3?MPa. We compiled the stress drop on the asperities together with a data set from previous studies of other inland earthquakes in Japan and elsewhere. The static stress drop on the asperity depends on its depth, and we obtained an empirical relationship between the static stress drop and the asperity??s depth. Moreover, surface-breaking asperities seemed to have smaller stress drops than buried asperities. Simple ground motion simulations using the characterized asperity source models reveal that deep asperities generate larger ground motion than shallow asperities, because of the different stress drops of the asperities. These characteristics can be used for advanced source modeling in strong ground motion prediction for inland crustal earthquakes.     

4010声发射仪的改装及其二维定位功能的初步实验

对仅具有线定位功能的4010声发射仪进行了改装,使之能进行面定位并能实时输出。增加同步时间控制器,使两组线定位器能保证对同一事件进行两个轴向同时定位;改单探头为多探头并联,并且在探头和样品间加声导板,减小定位误差,使非线性定位造成的误差与仪器的分辨精度一致,综合误差为样品尺度的1.5%;用Z-80单板机程控声发射仪和宽行打印机,将两个轴向的定位坐标和样品声发射率(或事件累集数,或能量值)实时打印。用改装后的4010声发射定位仪,在单轴压缩下对板状样品(30015013mm)进行直至破裂的定位实验,结果表明仪器的定位功能理想。同时,AE空间分布表明:于样品破坏前在崩裂型破裂区周围形成明显空区;AE密集区均被宏观张性或压性破裂贯穿;空间分布图象具有明显的阶段性。      

Dynamic faulting under rate-dependent friction

We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law.     

雁列式断层变形过程中物理场演化的实验研究（一）

在双轴压缩条件下对挤压型和拉张型雁列式断层变形破坏过程中声发射、断层位移及应变场的时空演化进行了实验研究。研究表明，两类雁列式断层具有类似的变形破坏过程，即前期以雁列区的破裂贯通为主，后期以沿断层的滑动为主。但挤压型雁列区内可积累较高的应变能，有较强的应变释放和声发射活动，雁列区对滑动始终起着阻碍作用；而拉张型雁列区内难以产生快速的应变释放和较强的声发射活动，雁列区对后期的滑动也无明显的阻碍作用。两类雁列区贯通所引起的失稳均有明显的前兆，而在随后滑动过程中失稳前兆的出现，取决于断层不同部位或段落之间是否有应变传递     

障碍体和凹凸体对地震矩的影响

本文利用具有非均匀剪切应力分布的圆盘状裂缝问题的解,定量计算了障碍体和凹凸体对地震矩的影响.我们发现,除非障碍体的尺寸很小或很接近于整个断层的尺寸,含有障碍体的断层的地震矩,大致为具有均匀应力降的等面积断层的地震矩的40%.由此可以推知,在障碍体的这种尺寸范围内,视应力降大致也是△Af=SA-f的40%.这里SA是远场应力,f是残余摩擦应力.另外,已存在滑动带对于凹凸体破坏的地震矩有明显的放大作用,这种作用在凹凸体尺寸较小时特别显著.例如,当凹凸体半径是整个断层半径的2/10时,凹凸体上的矩被放大7倍多;当其半径为整个断层半径的1/10时,凹凸体的矩将被放大约30倍. 文中还比较了三维效应和二维效应的差异.我们发现,Mx0/M0c在三维情形对障碍体尺寸的变化并不敏感,这是与二维情形很不相同的.这里,Mx0是具有障碍体的断层破裂所产生的地震矩,M0c是相应的均匀断层破裂的矩.所以,二维分析不适用于具有障碍体的三维断层.然而,在分析已存在滑动带对凹凸体破裂的矩的影响时,平面应变结果与三维结果差异并不很大.      

Hypocenter distribution and focal mechanism of AE events during two stress stage creep in Yugawara andesite

More than a thousand acoustic-emission (AE) hypocenters were determined in a cylindrical andesite specimen under two-stage uniaxial creep at stresses of 204 and 214 MPa. Strains were monitored for 6 peripheral points at the middle part of the cylindrical specimen's wall. The strain data indicate gradual increase of nonuniform deformation during steady creep and strong intensification of the nonuniformity during acceleration creep and, therefore, biased stress distribution within the specimen. The correlation between dilatant strain and AE hypocenters was investigated for whether or not tensile cracks emit AE. The region with high AE activity shows only a small dilatant strain. This negative correlation between AE and the dilatant region may eliminate tensile cracks as possible AE sources. A composite focal-mechanism solution of local AE events, covering a wide solid angle of the focal hemisphere, indicates that shear fractures emit AE waves. The direction of the compressional axis in this solution shows a significant deviation from that inferred from the applied external force, suggesting that the local stress field is governed by preexisting weak zones that are, presumably, produced by tensile cracks within the specimen. AE hypocenters tended to form clusters during steady creep under the constant compressional stress. During acceleration creep caused by a small step increase of the external stress, the preceding clusters disappeared while a new cluster appeared in an incipient fault plane. This suggests that changes in seismicity pattern such as migrations or quiescences of swarm—important clues for earthquake predictions—may be caused by an instantaneous change in the tectonic-stress levels.