EXPERIMENTAL STUDY OF TEMPERATURE EVOLUTION AND IDENTIFICATION OF INSTABILITY POSITION OF PLANAR STRIKE-SLIP FAULT DURING PROCESS OF STICK-SLIP

Stick-slip of fault in laboratory accompanies change of temperature. Temperature change is not only concerned with sliding friction, but also with the stress state of the sample. In this article, we use infra-red thermal imaging system as wide-range observation means to study the temperature variation of different stages during the deformation of sample. The rock sample for the experiment is made of granodiorite from Fangshan County with a size of 300mm×300mm×50mm. It is cut obliquely at an angle of 45°, forming a planar fault. Two-direction servo-control system was used to apply load on the sample. The load in both directions was forced to 5MPa and maintained constant (5MPa) in the X direction, then the load in the Y direction was applied by a displacement rate of 0.5μm/s, 0.1μm/s and 0.05μm/s successively. The left and below lateral of the sample were fixed, and the right and top lateral of the sample were slidable when loaded. The experiment results show not only the temperature change from increase to decrease caused by conversion of stress accumulation to relaxation before and after the peak stress, but also opposite variation of temperature increase on fault and temperature decrease in rock during instability stage. Most important of all, we have found the temperature precursor identifying the position of instability through the temperature variation with time along the fault. It shows that rate of temperature increase of instability position keeps relative high value since the stage of strongly off-linear stage, and accelerates in stage of meta-instability. After separating the effect of friction and stress, we found that temperature increase occurs in the rock near the fault instead of on the fault, which means the mechanism of temperature increase is stress accumulation. Temperature of fault at the instability position does not increase, which means the position is locked. We speculate that the position of locked area on fault with high stress accumulation near the fault may be the future instability position. It is of significance of studying temperature variation during stick-slip to the monitoring of earthquake precursors. Heat caused by friction of earthquake needs long time to transfer to the surface and could not be detected as a precursor. While the stress of surface rock near the fault would change as the stress of interior rock changes, which could cause detectable temperature variations. The research purpose of this article is to find special change positions before instability. As the temperature variations are caused by stress and slip of fault, the results are also meaningful to analysis of stress and displacement data related to earthquake precursors.     

Fracture in Westerly Granite under AE Feedback and Constant Strain Rate Loading: Nucleation, Quasi-static Propagation, and the Transition to Unstable Fracture Propagation

New observations of fracture nucleation are presented from three triaxial compression experiments on intact samples of Westerly granite, using Acoustic Emission (AE) monitoring. By conducting the tests under different loading conditions, the fracture process is demonstrated for quasi-static fracture (under AE Feedback load), a slowly developing unstable fracture (loaded at a `slow' constant strain rate of 2.5 × 10<sup>−6</sup> /s) and an unstable fracture that develops near instantaneously (loaded at a `fast' constant strain rate of 5 × 10⁻⁵ /s). By recording a continuous ultrasonic waveform during the critical period of fracture, the entire AE catalogue can be captured and the exact time of fracture defined. Under constant strain loading, three stages are observed: (1) An initial nucleation or stable growth phase at a rate of ~ 1.3 mm/s, (2) a sudden increase to a constant or slowly accelerating propagation speed of ~ 18 mm/s, and (3) unstable, accelerating propagation. In the ~ 100 ms before rupture, the high level of AE activity (as seen on the continuous record) prevented the location of discrete AE events. A lower bound estimate of the average propagation velocity (using the time-to-rupture and the existing fracture length) suggests values of a few m/s. However from a low gain acoustic record, we infer that in the final few ms, the fracture propagation speed increased to 175 m/s. These results demonstrate similarities between fracture nucleation in intact rock and the nucleation of dynamic instabilities in stick slip experiments. It is suggested that the ability to constrain the size of an evolving fracture provides a crucial tool in further understanding the controls on fracture nucleation.     

油气田典型岩石三轴压缩变形破坏与声发射活动特征——四川盆地震旦系白云岩及页岩的破坏过程

随着二氧化碳地质封存、深部地热开采、地下储气库建设、页岩气开发、二次驱油/驱气等工业应用的快速发展,与地下流体注入有关的诱发地震活动呈现一定的增加趋势.利用声发射实验观测油气田典型岩石在三轴压缩条件下变形破坏过程与声发射活动特征,对研究注水诱发地震过程有着重要意义.本文利用四川盆地现场采集的震旦系白云岩及页岩,采用实验室声发射技术观测研究岩石三轴压缩变形破坏过程中地震波速度等物性参数及声发射事件时空分布特征.实验结果表明:震旦系白云岩及页岩在变形破坏过程中均有一定的声发射活动.根据声发射定位结果,声发射主要集中在破坏前后的较短时间内,页岩的层理面为结构弱面,控制最终破坏面的形态及声发射特征.根据应力-应变结果,白云岩在压缩的后期阶段有一定的扩容现象,但页岩在整个压缩阶段均没有明显扩容现象.研究结果表明四川盆地较古老的白云岩及页岩具有脆性破坏特征,地下流体注入容易诱发微震活动,形成裂缝,有利于页岩气压裂开采.微震活动有利于监测裂缝的发生发展,但同时在页岩气开发及二氧化碳地质封存时应采取相应的预防控制措施进行安全合理的储盖层管理,避免灾害性诱发地震的发生.     

从断裂力学观点探讨 b 值的物理实质

岩石试验声发射及天然地震的 b 值到底反映了什么样的物理实质,是一个长期以来引起争论的问题.本文运用断裂力学的观点和方法,研究岩石试验试件内微裂纹系中各裂纹的扩展顺序及声发射能量,由此确定整个声发射序列,并进而确定 b 值.如果假设裂纹长度的分布密度函数为 p()=B-.经过力学分析和数学推导,可得到关系式 b=3/2.由此看出:b 值的物理实质是,它反映了介质的断裂构造状态————包括介质中裂纹系的空间分布及其它影响材料中裂纹扩展的物理量(如断裂韧度、摩擦系数等)的空间分布.这一结论与大多数实验及观测结果相符.我们的研究结果与茂木清夫观点有共同之处.茂木清夫认为:介质的不均匀性是决定 b 值的重要因素.介质的不均匀性一词固然也包括了介质的断裂构造状态在内,但是本文之观点比茂木更进一步说明了问题的实质.      

声发射矩张量反演

地震矩张量反演是获取震源过程的有效方法.岩石变形过程中的声发射与地震类似,均是弹性应变能快速释放.如假设条件得到满足,矩张量反演方法同样可用于了解声发射震源过程.声发射矩张量反演可使用P波位移进行计算.当样品尺寸远大于声发射波长,且声发射由微破裂产生,声发射源的尺度很小时,P波矩张量反演可采用远场近似.本文首先针对声发射的特点,实现根据远场P波反演声发射矩张量的算法,并通过人工声发射实验对算法的正确性和可靠性进行了检验.最后,用声发射矩张量反演方法对花岗岩单轴压缩实验的声发射源特征进行了分析.结果表明：对于纯剪切破裂模式,声发射矩张量反演可得断层面;对于非纯剪切破裂模式,如纵向挤压导致的横向张性劈裂,由于存在多解性不能得到断层面,但可通过矩张量的迹区分破裂模式.     

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值. 先期破碎带由于较低的破坏强度，其对差应力的微小变化具有特殊敏感性，从而成为源、兆分离、窗口或敏感点效应等地震活动性前兆现象的一种可能的原因.      

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,     

剪切破裂与粘滑——浅源强震发震机制的研究

周口店花岗闪长岩的高温高压三轴实验和理论分析表明,剪切破裂和摩擦滑移具有类似的孕育过程和发生机制。剪切破裂贯通强度就是一种摩擦强度。剪切破裂和摩擦滑移各自都有渐进式和突发式之分。突发式摩擦滑移是已有断层的粘滑滑移。突发式剪切破裂则是完整岩石的初始粘滑滑移。考虑到地壳温度随深度增加,完整岩石剪裂强震要求较高的围压,因此,多数浅源强震的发震方式很可能是已有断层的粘滑     

The Assessment of Damage Around Critical Engineering Structures Using Induced Seismicity and Ultrasonic Techniques

—?Two large-diameter boreholes have been excavated vertically from the floor of a tunnel at the Äspö Hard Rock Laboratory, Sweden. The two deposition holes will have simulated high-level radioactive waste canisters installed in them in an experiment undertaken to test the retrievability of waste from a proposed repository. Induced seismicity and other acoustic monitoring techniques have been used to investigate the Excavation Damaged Zone (EDZ) around the two holes. High-frequency acoustic emission (AE) monitoring has been used to delineate regions of stress-induced microfracturing on the millimetre scale. This has been shown to locate in clusters around the perimeter of the deposition hole at azimuths orthogonal to the far-field maximum principal stress. Three-dimensional velocity surveys have been conducted along ray paths that pass through the damaged region and through a stress-disturbed zone around the excavation. Induced microfracturing and stress disturbance have been observed as sharp decreases in velocity as the excavation proceeds through the rock mass. The combination of the high-resolution velocity measurements and the AE source locations has allowed the linking of the velocity measurements to a volume of excavation damaged rock. This has provided a quantitative estimate of the effect of the EDZ on the rock mass.     

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

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

标本非均匀性对岩石变形声发射时空分布的影响及其地震学意义

基于一组实验结果,讨论了岩石标本的非均匀性（预存微裂纹和宏观裂纹）对声发射 时空分布的影响. 对比分析表明,预存微裂纹的存在使得破裂成核前声发射率快速增加、 b值表现出一种短期异常现象,即在下降背景上出现起伏,从而增加了破裂时间的可预报 性. 宏观结构（节理、层面等）的存在对声发射率和b值的影响与微裂纹相同,而且宏 观结构对声发射的空间分布具有控制作用,声发射丛集的宏观构造通常控制着未来的主破裂 . 这意味着宏观构造的存在使得主破裂的时间和位置预报成为可能.      

固体围压下完整花岗岩粘滑现象的实验研究

本文用完整的花岗岩样品在固体围压三轴实验装置上压缩,研究围压和应变速率对岩样变形破裂过程、粘滑应力降、粘滑复发间隔及样品主破裂几何分布的影响。结果表明,加载速率较低时,粘滑应力降较大,复发间隔较长且分布无规律。加载速率越大,粘滑应力降越小,复发间隔也近似相等,粘滑事件表现出准周期性。围压和应变率较低时,岩石的主破裂会演变成两个交叉的共轭断裂面;而围压和应变率较高时,岩石的主破裂则演变成单一断裂面或入字形断裂面。本文结果对认识中国大陆板内地震孕育、发生及重复过程;研究地震重复发生的机制及影响地震复发间隔的主要因素都有重要意义。     

广州地磁台转换函数的长期变化和季节变化

本文全面系统地研究了1960——1987年广州地磁台转换函数的变化.结果表明:(1)季节变化的特点是夏季幅度大,冬季幅度小,具有12和6个月的主要周期成份;长期变化则表现以每年0.0025的速率下降.转换函数 B 的长期变化和季节变化不明显;(2)由转换函数 A,B 求出的帕金森矢量的方向明显地受海洋的影响,反映了海洋效应;(3)转换函数 A,B 的逐月值与地磁活动相关,而年均值与地磁活动无关.A 的年均值无明显的11年或22年周期.长期变化可能主要受地球内部电性结构的影响.      

亚失稳阶段雁列断层热场演化的实验研究

寻找地震前兆是国内外地震预报的难点.在实验室开展岩石变形实验,观测相关物理量的演化特征和规律,是前兆探索中重要的基础工作.本文对压性雁列断层进行变形实验,根据应力时间曲线将具有前兆意义的变形阶段分为强偏离线性阶段、亚失稳阶段和失稳阶段,利用红外热像仪记录各个阶段的热场变化,特别是亚失稳阶段的热场表现.研究结果显示:强偏离线性阶段,岩桥区挤压升温;亚失稳阶段,断层带增温,且高温点扩展连通,断层带内侧卸载降温,断层带上的协同化增温和断层带内侧(包含岩桥区)的协同化降温是失稳前的重要信号;失稳阶段,应力释放,标本整体降温,只有断层带因摩擦升温.总之,亚失稳阶段的热场具有明显特征,利于识别出亚失稳态;不同变形阶段、不同构造部位温度变化不同,在寻找与构造活动相关的热异常时,应考虑变形所处的阶段及异常所处的构造部位.     

三向受压状态下地壳岩体的张性变形和张性破裂

本文叙述了处于三向受压状态的地壳岩体产生张性变形和张性破裂的原理,并以珠江三角洲地区为例说明:①岩体内应力集中点的浅部,可产生张性变形和张性破裂;②地质考察所观察到的张性或张剪性断裂,其断裂面上作用的正应力未必是拉应力,极有可能是压应力,但却一定是拉应变;③张性破裂引起的地震,震级较小,一般在4.5级以下,④地震引起的张性地裂缝的走向,就是震源应力场的主压应力方向。     

水氡异常的水动力学机制

在水氡异常机制研究中，人们以往都强调岩石力学机制，即孕震过程中岩石微破裂产生氡射气。然而，大量的震例表明，水氡异常可出现在孕震早期，可见于距震中较远的井中，那里的与那时的含水层岩石很难处在微破裂状态。因此，笔者提出水氡异常的水动力学机制，即在较低的应力水平下含水岩体的变形与其水动力状态的变化引起水氡异常的机制     

煤岩变形破裂电磁辐射的实验研究

对受载煤体变形破裂产生电磁辐射进行了实验研究及规律分析,并对煤岩变形破裂电磁辐射的机理进行了探讨分析．研究结果表明,煤体变形破裂时能够产生电磁辐射,电磁辐射信号与声发射信号并非严格同步,电磁辐射信号较声发射信号丰富,电磁辐射与煤岩体的变形破裂过程密切相关．煤岩电磁辐射技术在揭示煤岩变形破裂机理,预测、预报地震,矿山煤岩灾害动力现象及岩石混凝土建筑稳定性方面有着非常广泛的应用前景．     

In‐situ stress at a site close proximity to the Gofukuji Fault,central Japan,measured using drilling cores

The stress state near the Gofukuji Fault, central Japan, was estimated by simultaneously applying deformation rate analysis (DRA), the acoustic emission (AE) method, and AE rate analysis (AERA) to drilling core samples recovered from depths of 327 and 333 m at a site close to the fault. The obtained stress state was for a strike‐slip fault. It is predicted from the estimated stresses that a tensile stress of 6.4 ± 1.1 MPa acted on the borehole wall at 333 m depth during drilling. This is comparable to the typical tensile strength of granodiorite. The tensile stress estimated at 327 m depth is less than the tensile strength. This is consistent with observations of drilling‐induced tensile fracture (DITF) at depths between 329 and 334 m. Ratios of the shear stress to normal stress (normalized shear stress) acting on the fault are calculated to be 0.4–1.0, which are comparable to friction coefficients of rocks measured in laboratories. The large normalized shear stress may represent strength recovery of the Gofukuji Fault. The impendency of earthquake occurrence on this fault is evaluated to be high from the recurrence interval of earthquakes on the fault.     

Mechanisms of slip nucleation during earthquakes

Slip nucleation during earthquakes is apparently analogous to rupture nucleation within an intact rock sample subjected to triaxial loading. The observations indicate that both these nucleation processes initiate within a relatively small volume and in both the slip propagates unstably along a quasi-planar surface. In both processes a single, pre-existing, shear fracture cannot nucleate the large-scale slip, and in both a ‘process zone' that includes several interacting fractures in a small volume are required to initiate the unstable slip. Both processes require rupture of intact rocks, generate complex fracture geometry, and are associated with intense energy-release rate during slip. Recent observations and analyses are used to correlate rupture nucleation in laboratory tests with nucleation events of large earthquakes. It is proposed that earthquake nucleation occurs by the interaction among multiple fractures within a small volume that develops into unstable yielding of the healed fault zone.