岩盐断层带摩擦滑动的速度依赖性转换及其地震学意义

为了深入理解断层带摩擦滑动速度依赖性转换及其机制,利用双轴摩擦实验对干燥及含水条件下岩盐断层带摩擦的速度依赖性进行了实验研究,并观测了摩擦滑动过程中的声发射,分析了断层带的微观结构.实验结果表明,干燥岩盐断层带在0.1~100 μm/s的速度范围内表现为速度弱化,增大σ₂会使断层带向速度强化转变;含水条件下岩盐断层带在1~100 μm/s的速度范围内表现为速度弱化,而在0.1~0.01 μm/s的速度范围内表现为速度强化,速度依赖性转换出现在0.1~1 μm/s,其中断层表现为振荡或应力释放时间较长的黏滑事件;岩盐断层带在干燥条件下表现出很强的声发射活动,每个黏滑均对应一丛声发射事件,而在含水条件下一次黏滑只对应一个声发射事件.显微观察表明,局部化的脆性破裂是速度弱化域的主要变形机制,分布式的碎裂流动是干燥岩盐断层带在速度强化域的变形机制,颗粒边界迁移以及压溶作用的塑性变形是含水条件下岩盐断层带在速度强化域的主要变形机制,而脆性破裂和塑性变形共同控制着速度依赖性转换域断层带的变形.水的存在促进岩盐发生塑性变形,进而导致断层带从速度弱化向速度强化转换.上述结果有助于理解断层带上地震活动的特征和慢地震的机制.     

白云母对岩盐断层带摩擦速度依赖性影响的实验研究

为更好地理解层状硅酸盐对断层强度、滑动速度依赖性及地震活动特征的影响,利用双轴摩擦实验对含白云母岩盐断层带在干燥及含水条件下摩擦的速度依赖性进行了实验研究,并观测了摩擦滑动过程中的声发射,分析了断层带的微观结构.实验结果表明,干燥条件下含白云母岩盐断层带在0.1 ～ 100μm/s的速度范围内表现为黏滑和速度弱化,增大σ2会使断层带从速度弱化向速度强化转化,速度依赖性转换出现在0.1 μm/s,其中断层滑动表现为稳滑或应力释放时间较长的黏滑事件;含水条件下含白云母岩盐断层带在0.05 ～0.01μm/s的速度范围内表现为速度强化,0.1 ～10μm/s的速度范围内表现为速度弱化,50～100μm/s的速度范围内又转换为速度强化行为.含白云母岩盐断层带在干燥条件下一次黏滑伴随一个或一丛声发射事件,而在含水条件下与稳滑相对应,滑动过程中并未记录到声发射事件.显微结构观察表明,速度弱化域的主要变形机制是岩盐颗粒的脆性破裂和局部化的滑动;干燥条件下,速度强化域的主要变形机制是岩盐颗粒的均匀破裂;含水条件下2个速度强化域对应不同的微观机制,高速域的速度强化受控于岩盐颗粒在白云母相互连结形成的网状结构上的滑动及其均匀碎裂作用,而低速域的速度强化还受岩盐的压溶作用控制.通过与岩盐断层带摩擦实验结果对比可知,白云母的存在对于燥岩盐断层带摩擦滑动方式和速度依赖性没有显著影响,而在含水条件下白云母的存在使得岩盐断层带滑动趋于稳定.实验结果为分析含层状硅酸盐断层的强度和稳定性提供了依据.此外,在速度依赖性转换域上观察到的应力缓慢释放的现象进一步证实了在岩盐断层带摩擦滑动过程中观察到的现象,这对慢地震机制研究具有参考意义.     

流体对石灰岩断层摩擦滑动影响的实验研究

在气体介质三轴高温岩石力学实验仪器上,采用意大利Scaglia Bianca石灰岩,在温度50～300℃、围压150MPa,含50MPa孔隙压、无孔隙压含饱和水和完全干燥三种条件下,开展摩擦滑动实验.实验力学数据和显微结构表明,完全干燥样品在120℃时出现慢滑移,实验样品中没有出现溶解与沉淀.无孔隙压含饱和水条件下,100℃、120℃、150℃条件下出现典型的慢滑移,实验样品中含有微弱的溶解与沉淀;300℃条件下出现黏滑,实验样品中出现沉淀.在含50MPa孔隙压条件下,50℃时的实验表现为典型的稳滑,实验样品中以溶解为主;在100～150℃时,出现慢滑移,实验样品中以溶解为主,沉淀为辅;在200～300℃时,出现典型黏滑,实验样品中以沉淀为主.实验结果表明,石灰岩断层泥摩擦滑动稳定性随温度变化,受流体中碳酸钙的溶解和沉淀作用控制,因此,流体中矿物质的饱和度这一化学性质对断层带的摩擦强度和摩擦滑动稳定性具有显著影响.     

硬石膏断层带摩擦稳定性转换与微破裂特征的实验研究

断层带摩擦稳定性转换及其对应的微破裂特征对于地震成核条件和慢地震机理研究具有重要的意义.本文利用双轴实验装置研究了硬石膏断层带摩擦稳定性的转换及其对应的应变变化、微破裂特征,并分析了实验标本的微观结构.实验结果表明,σ₂和加载点速度对断层滑动稳定性具有显著影响.在低σ₂条件下,硬石膏断层带出现不稳定滑动,变形以局部化的脆性破裂和摩擦为主;随σ₂的增加,断层由不稳定滑动向稳定滑动转换,断层带变形方式逐渐转变为分布式的破裂.在低σ₂条件下,硬石膏断层带在较低加载点速度下表现为速度强化且滑动稳定,在中等加载点速度下表现为速度弱化并伴有准周期性的黏滑,在较高加载点速度下又有转向速度强化的趋势,σ₂的提高使得速度弱化的范围逐渐减少,滑动趋于稳定.上述两次转换对应不同的微破裂特征,在较高速度下从速度弱化转换为速度强化时,断层滑动伴有能量较小但频度很高的微破裂活动,而在较低速度下从速度弱化转换为速度强化时,断层滑动伴有间歇性的微破裂,这与断层带的微观结构特征有较好的对应关系,表明其转换机制是不同的.     

水热条件下富层状硅酸盐矿物糜棱岩的摩擦特性实验研究

为探索断层岩石摩擦特性对于断层力学性质的影响,我们采集了龙门山汶茂断裂韧性剪切带中的富含层状硅酸盐矿物的糜棱岩样品进行了水热条件下摩擦滑动实验研究.实验在三轴压机之上完成,实验温度为100~600℃,有效正压力100MPa,孔隙水压分别为30MPa和130MPa.为获得摩擦滑动的稳定性参数(a-b),剪切滑移速率在1.22μm·s-1,0.244μm·s-1和0.0488μm·s-1之间切换.实验发现在200~500℃的温度范围内,摩擦系数随着温度的增加而显著增大(约0.56~0.72).在200~300℃范围内,随温度的升高糜棱岩的摩擦滑动表现出由稳定的速度强化向不稳定速度弱化转变的趋势.在有效正压力不变的情况下,孔隙水压的增大会促进糜棱岩的摩擦滑动在500~600℃温度范围内由不稳定的速度弱化向稳定的速度强化的转变.实验给出的断层在原地深度处的脆性和塑性变形机制的转变,有助于理解断层深部的地震成核机制以及成核的温压条件.     

蒙脱石的脱水作用对断层摩擦本构行为的影响

利用高温双轴摩擦装置，研究了含蒙脱石的断层带在不同温度下摩擦滑动的速度依赖性，以期了解脱水作用对摩擦行为的影响。结果表明，断层带摩擦强度随温度而升高，而速度依赖性较为复杂，以１．４ｕ／ｓ为界，室温和１００℃时，低滑动速率下表现为微弱的速度弱化，高滑动速率下则表现为速度强化；２００℃时均为速度强化；３００℃时高滑动速率下仍为速度强化，但低滑动速率下转变为速度弱化；４００℃以上，均为明显的速度弱化。摩擦行为的变化与脱水过程及相应的断层物质变形方式的变化密切相关     

脆塑性转化带Carrara大理岩断层稳定性的实验研究

为探究脆塑性转化带断层的力学性质和滑动稳定性,本文采用干燥的Carrara大理岩预切断层(saw-cut)样品,在气体介质三轴高温岩石力学实验仪上开展了摩擦实验研究,实验温度70～400℃,围压30～100 MPa,位移速率在0.08μm·s^-1, 0.4μm·s^-1, 2μm·s^-1之间切换.实验力学数据揭示,不同围压下Carrara大理岩断层摩擦系数随温度变化规律不同：低围压(30 MPa)下,摩擦系数随温度升高先增大后减小,中高围压(≥70 MPa)下摩擦系数则表现为随温度先减小后增大.断层摩擦滑动行为在100～300℃的范围内表现出由稳定的速度强化转化为不稳定的速度弱化,且在400℃左右重新转变为稳定的速度强化.实验后断层滑动面形貌和微观结构分析表明,稳定滑动断层面为高反射镜面,擦痕清晰;黏滑断层面为有光泽的凹凸不平的表面;最高围压下蠕滑的断层面粗糙无光泽,擦痕不可辨别.本文认为受温度激活的塑性变形过程逐步主导了岩石变形,对断层激活发生不稳定滑动至关重要,而高围压则会抑制断层的不稳定滑动.本研究结果不仅为识别野外...     

热水条件下黑云母断层泥的摩擦强度与稳定性

黑云母是自然界常见的层状硅酸盐矿物,其摩擦系数不高且化学稳定性好,对其摩擦性质的关注可能会对弱断层的研究有所帮助.本次工作选取的实验温度条件对应于典型地壳强度模型中脆塑性转化带的范围,为300 ℃和400 ℃.有效正应力为200 MPa,孔隙水压包括10 MPa和30 MPa,在此条件下对黑云母模拟断层泥进行摩擦实验研究.实验得出黑云母的摩擦系数平均在0.36左右.速度依赖性随温度升高速度弱化的程度增强,表现为300 ℃为十分微弱的速度弱化,而在400 ℃出现了黏滑行为,代表了更强的速度弱化.显微结构中同时出现了脆性剪切变形和塑性扭折变形,但决定宏观力学性质的显然是脆性剪切变形.在黑云母存在的情况下,本研究的实验结果有助于理解大陆地壳脆塑性转化带中地震的可能性和弱断层深部的变形机制、宏观力学行为以及地震活动.     

基岩区断层黏滑与蠕滑的地质标志 和岩石力学实验证据

文中总结了基岩断层带黏滑与蠕滑的地质标志与岩石力学实验证据,分析了控制黏滑与蠕滑的物理机制。断层带内的矿物组成、矿物变形机制、流体作用和断层带变形方式等是控制黏滑与蠕滑的主要因素。富含黏土矿物的断层泥具有速度强化型摩擦滑动,控制着断层蠕滑,而以方解石、石英、长石及辉石等造岩矿物为主的断层泥在大陆浅源地震的震源深度条件下具备黏滑条件。脆性破裂伴随的扩容过程是断层黏滑的必要条件,而压实、碎裂和塑性剪切变形形成的叶理和小褶皱对应于蠕滑。在流体作用下,压溶使孔隙和微裂隙愈合,有利于断层强度的恢复和断层闭锁,既是断层发生不稳定滑动的根源,也是断层带局部存在高压流体的条件,而在流体作用下的退变质反应与水解反应生成黏土矿物和层状及环状硅酸盐矿物,不仅降低了断层带的强度,还导致断层向蠕滑转变。断层带内均匀分布多个剪切面和较宽的变形带对应于蠕滑,局部化的R剪切及Y剪切、窄变形带和摩擦镜面对应于黏滑。     

几种模拟断层泥摩擦滑动速度依赖性转换的实验研究

基于速率-状态摩擦定律,速度弱化是断层失稳的必要条件,速度弱化向速度强化的转换控制着断层从不稳定滑动向稳定滑动的转换.因此,断层摩擦滑动的速度依赖性转换是涉及断层带上地震活动特征、地震成核深度等的重要问题.     

Simulation of the Influence of Rate- and State-dependent Friction on the Macroscopic Behavior of Complex Fault Zones with the Lattice Solid Model

-- In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by Dieterich and Kilgore (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D_c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.     

八宝山断层的变形行为与降雨及地下水的关系

利用北京丰台大灰厂观测台站1970——2003年的长期综合观测资料,系统分析了降雨和地下水对北京八宝山断层变形行为的影响. 研究结果表明:当降雨量持续稳定周期性变化,则断层孔隙压也呈稳定周期性变化,断层变形行为也表现出明显的规则周期性变化;当降雨量明显偏低或偏离正常周期性变化,则断层孔隙压周期性消失,断层的变形行为方式也发生改变. 降雨通过改变断层带孔隙压力的变化而影响着八宝山断层的变形行为. 结合该断层已有研究结果,认为降雨和地下水有可能通过改变断层变形行为方式而影响着区域构造应力/应变场的调整. 这一结果将对研究地球浅部流体与固体相互作用提供直接的观测证据.      

FRICTIONAL PROPERTIES OF A NATURAL FAULT GOUGE FROM DRILLED CORES IN THE LONGMENSHAN FAULT ZONE CUTTING GRANITIC ROCK

In this paper, we report friction experiments performed on natural fault gouge samples embedded in granitic rock from drilled core by a project entitled "the Longmenshan Fault Shallow Drilling(LMFD)". Compared with other natural fault gouge, this yellow-greenish gouge(YGG)is dominantly chlorite-rich. The maximum content of chlorite reaches 47%in the YGG. To understand the frictional properties of the YGG sample, experiments were performed at constant confining pressure of 130MPa, with constant pore pressure of 50MPa and at different temperatures from 25℃ to 150℃. The experiments aim to address the frictional behavior of the YGG under shallow, upper crustal pressure, and temperature conditions. Compared with previous studies of natural gouge, our results show that the YGG is stronger and shows a steady state friction coefficient of 0.47~0.51. Comparison with previous studies of natural gouge with similar content of clay minerals indicates a sequence of strengths of different clay minerals:chlorite > illite > smectite. At temperatures up to 150℃ hence depths up to~8km in the Longmenshan region, the YGG shows stable velocity-strengthening behavior at shallow crustal conditions. Combined with the fact of strong direct velocity effect, i.e., (a-b)/a>0.5, faults cutting the present clastic lithology up to~8km depth in the Longmenshan fault zone(LFZ)are likely to offer stable sliding resistance, damping co-seismic rupture propagating from below at not-too-high slip rates. However, as the fault gouge generally has low permeability, co-seismic weakening through thermal pressurization may occur at high slip rates(>0.05m/s), leading to additional hazards.     

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.     

Discrete Element Modeling of Stick-Slip Instability and Induced Microseismicity

Using Particle Flow Code, a discrete element model is presented in this paper that allows direct modeling of stick-slip behavior in pre-existing weak planes such as joints, beddings, and faults. The model is used to simulate a biaxial sliding experiment from literature on a saw-cut specimen of Sierra granite with a single fault. The fault is represented by the smooth-joint contact model. Also, an algorithm is developed to record the stick-slip induced microseismic events along the fault. Once the results compared well with laboratory data, a parametric study was conducted to investigate the evolution of the model’s behavior due to varying factors such as resolution of the model, particle elasticity, fault coefficient of friction, fault stiffness, and normal stress. The results show a decrease in shear strength of the fault in the models with smaller particles, smaller coefficient of friction of the fault, harder fault surroundings, softer faults, and smaller normal stress on the fault. Also, a higher rate of displacement was observed for conditions resulting in smaller shear strength. An increase in b-values was observed by increasing the resolution or decreasing the normal stress on the fault, while b-values were not sensitive to changes in elasticity of the fault or its surrounding region. A larger number of recorded events were observed for the models with finer particles, smaller coefficient of friction of the fault, harder fault surroundings, harder fault, and smaller normal stress on the fault. The results suggest that it is possible for the two ends of a fault to be still while there are patches along the fault undergoing stick-slips. Such local stick-slips seem to provide a softer surrounding for their neighbor patches facilitating their subsequent stick-slips.     

Frictional Properties of a Low-Angle Normal Fault Under In Situ Conditions: Thermally-Activated Velocity Weakening

The Zuccale fault is a regional, low-angle, normal fault, exposed on the Isle of Elba in central Italy that accommodated a total shear displacement of 6–8 km. The fault zone structure and fault rocks formed at <8 km crustal depth. The present-day fault structure is the final product of several deformation processes superposed during the fault history. In this study, we report results from a series of rotary shear experiments performed on 1-mm thick powdered gouges made from several fault rock types obtained from the Zuccale fault. The tests were done under conditions ranging from room temperature to in situ conditions (i.e., at temperatures up to 300 °C, applied normal stresses up to 150 MPa, and fluid-saturated.) The ratio of fluid pressure to normal stress was held constant at either λ = 0.4 or λ = 0.8 to simulate an overpressurized fault. The samples were sheared at a constant sliding velocity of 10 μm/s for at least 5 mm, after which a velocity-stepping sequence from 1 to 300 μm/s was started to determine the velocity dependence of friction. This can be represented by the rate-and-state parameter (a–b), which was determined by an inversion of the data to the rate-and-state equations. Friction of the various fault rocks varies between 0.3 and 0.8, similar to values obtained in previous studies, and decreases with increasing phyllosilicate content. Friction decreases mildly with temperature, whereas normal stress and fluid pressure do not affect friction values systematically. All samples exhibited velocity strengthening, conditionally stable behavior under room temperature conditions and (a–b) increased with increasing sliding velocity. In contrast, velocity weakening, accompanied by stick–slips, was observed for the strongest samples at 300 °C and sliding velocities below 10 μm/s. An increase in fluid pressure under these conditions led to a further reduction in (a–b) for all samples, so that they exhibited unstable, stick–slip behavior at low sliding velocity. The results suggest that phyllosilicate-bearing fault rocks can deform by stable, aseismic creep at low, resolved shear stress and low shear rate conditions. An increase in fluid pressure or loading of stronger portions could lead to a runaway instability. The runaway instability might be limited in size because of (1) the fault heterogeneity, (2) the observed strengthening at higher sliding velocities, and (3) a co-seismic drop in pore-fluid pressure.     

Spatio-temporal Slip, and Stress Level on the Faults within the Western Foothills of Taiwan: Implications for Fault Frictional Properties

We use preseismic, coseismic, and postseismic GPS data of the 1999 Chi-Chi earthquake to infer spatio-temporal variation of fault slip and frictional behavior on the Chelungpu fault. The geodetic data shows that coseismic slip during the Chi-Chi earthquake occurred within a patch that was locked in the period preceding the earthquake, and that afterslip occurred dominantly downdip from the ruptured area. To first-order, the observed pattern and the temporal evolution of afterslip is consistent with models of the seismic cycle based on rate-and-state friction. Comparison with the distribution of temperature on the fault derived from thermo-kinematic modeling shows that aseismic slip becomes dominant where temperature is estimated to exceed 200° at depth. This inference is consistent with the temperature induced transition from velocity-weakening to velocity-strengthening friction that is observed in laboratory experiments on quartzo-feldspathic rocks. The time evolution of afterslip is consistent with afterslip being governed by velocity-strengthening frictional sliding. The dependency of friction, μ, on the sliding velocity, V, is estimated to be ${{\partial \mu }/{\partial \, {\rm ln}\, V}} = 8 \times 10^{ - 3}$ . We report an azimuthal difference of about 10–20° between preseismic and postseismic GPS velocities, which we interpret to reflect the very low shear stress on the creeping portion of the décollement beneath the Central Range, of the order of 1–3 MPa, implying a very low friction of about 0.01. This study highlights the importance of temperature and pore pressure in determining fault frictional sliding.     

SIMULATION OF SEISMIC RISK IN THE DALIANGSHAN SUB-BLOCK AND ADJACENT AREAS USING THE NONLINEAR FRICTION FEM METHOD

Most earthquakes result from fault activity under heterogeneous loading and complex physical properties, also affected by fault structure and interaction between faults. Such a complicated mechanism makes often failures of the "seismic gap" theory in the effort of medium-and long-term earthquake prediction. This study attempts to address this issue using the finite element method(FEM).The friction behavior of faults can be used to simulate the non-uniformity of rupture processes of the seismogenic structure. So we use the FEM containing non-linear friction to simulate fault ruptures in the Daliangshan sub-block and adjacent areas, and compare the results with time-space evolution of historical M_S ≥ 7 earthquakes since 1840 in this region. In the simulation, the sequence of large-batch fault contact nodes change from "stick state" to "slip state" in short time, which mimics the sudden fault slip and the occurrence of major earthquakes. The results show that the fault breaking lengths from simulation are largely consistent with the magnitudes of historical earthquakes in the study area, such as the 1850 Puge-Xichang M_S7.5, and 1887 Shiping M_S7.0 earthquakes. The simulation also shows the development of seismic gaps and "gap breaks" by major earthquakes on the Xianshuihe fault, such as 1955 Kangding M_S7.5 earthquake. Especially, the results illustrated the very long time of the seismogenic process of the 2008 Wenchuan M_S8.0 earthquake, and the corresponding sudden big rupture along the Longmenshan Fault, which is very similar to the observed surface rupture and very long incubation time and sudden co-seismic process. Then, this simulation is further applied to long-term earthquake prediction for the study area by calculation on a much longer time. The simulation results suggest that the Xiaojiang fault and the Zemuhe fault have relatively higher seismic risk, while moderate-sized earthquakes might occur on the Daliangshan fault and the Aninghe fault, and major earthquakes might rupture the northern segment of the Xianshuihe fault in a much longer time.