大尺度雁列断层粘滑运动数值模拟

通过ANSYS软件建立二维雁列断层模型,采用粘弹性力学参数,建立摩擦系数接触单元,模拟10万年时间大尺度雁列走滑断层活动（地震相关活动）,分析挤压条件下的断层活动特征,研究结果表明：①断层时空区域不同,所受应力也不同,接近贯通区域重合地区所受应力较大,粘滑运动剧烈,由贯通区域向断层两侧逐渐延伸,所受应力依次减小;②摩擦系数越小,地震发生频率越高,地震运动周期越短,则位错量相对较小,易发生震级较小地震;摩擦系数越大,地震发生频率越低,地震运动周期越长,则位错量相对较大,易发生震级较大地震;③给定的边界挤压速率越大,地震发生频率越高,地震运动周期越短,则位错量增大,易发生震级较大地震。     

应力途径与岩石的摩擦滑动

用双剪法对济南辉长岩和点苍山大理岩进行了摩擦滑动实验.实验中的应力变化方式有两种.A 型实验:先使断层面上的正应力增加到一定值,然后保持正应力不变,并增加剪应力使断层发生粘滑;B 型实验:先使断层面上的正应力增加到一定的值,保持正应力不变,并增加剪应力到断层发生粘滑前的某一应力状态,再保持剪应力不变,减小正应力直到粘滑发生.实验表明,B 型实验中岩石的摩擦强度高于 A 型实验.A 型实验中粘滑发生前有声发射率增加的前兆,B 型实验中粘滑发生前看不到声发射率的明显增加.由实验得到一个启示,即闭锁断层的开锁可能采取两种形式:一种是冲开闭锁,即以剪应力的增加使断层发生错动;另一种是解开闭锁,即以正应力的减小使断层发生错动.断层的粘滑采取哪种形式,由断层带的应力变化途径决定.      

Experimental study on nucleation process of stick-slip instability on homogeneous and non-homogeneous faults

The nucleation process of stick-slip instability was analyzed based on the experimental measurements of strain and fault slip on homogeneous and non-homogeneous faults. The results show that the nucleation process of stick-slip on the homogeneous fault is of weak slip-weakening behavior under constant loading point velocity. The existence of a short “weak segment” on the fault makes slip-weakening phenomenon in nucleation process more obvious, while the existence of a long “weak segment” on the fault makes the nucleation process changed. The nucleation is characterized by accelerating slip in a local region and rapid increase of shear stress along the fault in this case, which is more coincident with the rate and state friction law. During the period when fault is locked, increasing of shear stress causes lateral elastic dilation near the fault, and the rebound of the dilation at the time of instability causes an instantaneous increase of normal stress in the fault plane, which is an important factor making fault be rapidly locked and its strength recovered.     

Experimental study on nucleation process of stick-slip instability on homogeneous and non-homogeneous faults

The nucleation process of stick-slip instability was analyzed based on the experimental measurements of strain and fault slip on homogeneous and non-homogeneous faults. The results show that the nucleation process of stick-slip on the homogeneous fault is of weak slip-weakening behavior under constant loading point velocity. The existence of a short "weak segment" on the fault makes slip-weakening phenomenon in nucleation process more obvious, while the existence of a long "weak segment" on the fault makes the nucleation process changed. The nucleation is characterized by accelerating slip in a local region and rapid increase of shear stress along the fault in this case, which is more coincident with the rate and state friction law. During the period when fault is locked, increasing of shear stress causes lateral elastic dilation near the fault, and the rebound of the dilation at the time of instability causes an instantaneous increase of normal stress in the fault plane, which is an important factor making fault be rapidly locked and its strength recovered.     

An experimental study of triggered stick-slip

Stick-slips have been studied in the laboratory on granite, labrodorite and sandstone samples of two different sizes. Different roughness was achieved on the sawcut surfaces by finishing them with different grinding compounds ranging from grit 40 to grit 1000. Stick-slips occurred as a result of 1) slowly increasing the shear and normal stresses, 2) superimposing a sinusoidal stress modulation (0.1 and 10 Hz) on the slowly increasing stresses, 3) triggering by a stress impulse when the shear stress was well below the levels where stick-slips occurred without the impulse, and 4) foreshocks.Stick-slips triggered by impulses or foreshocks occurred long after the beginning of the triggering events, i.e., long in comparison with elastic wave travel times through the sample. All triggered events were very rich in high frequencies (corner frequency of 100 kHz). The untriggered stick-slips did not contain much energy at the high frequencies (corner frequency of 10 kHz). The dynamic friction coefficients for the triggered stick-slips were smaller than for the untriggered events.The long delay between the onset of the trigger and the stick-slip, and the high frequencies may be a consequence of corrosion of asperities. The ultimate triggering and the rate of corrosion are likely related to the interplay of the normal and shear stresses as they load and unload the fault surface. The consistent shape of the high frequency spectra is probably due to sample resonances which are excited rather than being characteristic of the details of the stick-slips. If these laboratory observations are directly applicable to earthquake seismology, the spectra of earthquakes which were triggered by other earthquakes should be anomalously rich in high frequencies.     

Dynamic Rupture in a 3-D Particle-based Simulation of a Rough Planar Fault

An appreciation of the physical mechanisms which cause observed seismicity complexity is fundamental to the understanding of the temporal behaviour of faults and single slip events. Numerical simulation of fault slip can provide insights into fault processes by allowing exploration of parameter spaces which influence microscopic and macroscopic physics of processes which may lead towards an answer to those questions. Particle-based models such as the Lattice Solid Model have been used previously for the simulation of stick-slip dynamics of faults, although mainly in two dimensions. Recent increases in the power of computers and the ability to use the power of parallel computer systems have made it possible to extend particle-based fault simulations to three dimensions. In this paper a particle-based numerical model of a rough planar fault embedded between two elastic blocks in three dimensions is presented. A very simple friction law without any rate dependency and no spatial heterogeneity in the intrinsic coefficient of friction is used in the model. To simulate earthquake dynamics the model is sheared in a direction parallel to the fault plane with a constant velocity at the driving edges. Spontaneous slip occurs on the fault when the shear stress is large enough to overcome the frictional forces on the fault. Slip events with a wide range of event sizes are observed. Investigation of the temporal evolution and spatial distribution of slip during each event shows a high degree of variability between the events. In some of the larger events highly complex slip patterns are observed.     

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

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

Frictional heat generation and seismic radiation in a foam rubber model of earthquakes

Results from a study of stick-slip particle motion at the interface between two stressed foam rubber blocks indicate that normal vibrations and interface separation are an important part of the stick-slip process in foam rubber. The dimension of the dynamic slip pulse is small compared to the dimension of the model (approximately 10 cm vs. 200 cm) consistent with the abrupt-locking slip pulse model ofBrune (1970, 1976), andHeaton (1990). A comparison of frictional heat generation between stable-sliding and stick-slip foam rubber models indicates a linear relation between the temperature increase on the fault surface (for a given distance of slip) and the driving shear force for the stable-sliding model, while for the stick-slip model there is essentially no variation in frictional heat generation with an increase in shear stress. We performed experiments to investigate the ratio of normal motion to shear motion at different levels of normal stress in the stick-slip foam rubber model. Preliminary result indicate that the normal component of the particle motion increases more rapidly with increasing normal stress than the shear component. The phenomenon of interface separation and normal vibrations may thus explain some of the most frustrating problems in earthquake mechanics, e.g., the heat flow paradox, the long-term weakness of major active faults, and anomalousP-wave radiation.     

Effect of fault jogs on frictional behavior: An experimental study

Studying the effect of geometrically irregular bodies on the mechanical behavior of fault activity is of significance in understanding the seismic activity along a fault zone. By using rock mechanics ex- periment with medium-scale samples, we have studied the effect of fault jogs, the most common irregularity along fault zones, on frictional behavior. The research indicates that extensional fault jog can be easily fractured because of its low strength and the fractured jog has no obvious resistance to fault sliding, and the micro-fractures occurring in the jog are indicative of stick-slip along the faults. The fault zone containing extensional jogs is characterized by velocity weakening and can be described by rate and state friction law. Compressional fault jog makes fault sliding more difficult because of its high fracturing strength, but the micro-fractures occurring in the tensile areas around fault ends at higher stress level can provide necessary condition for occurrence of stick-slip along the faults before the jog is fractured and thus act as precursors of fault instability. Compression jog can be taken as a stable indicator of fault segmentation until the jog is completely fractured and two faults are linked.     

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.     

Stick-slip and strain waves in the physics of earthquake rupture: experiments and models

The study of stick-slip is directly connected with the formulation of the consistent concept of seismicity migration and explanation of the nature of strain waves observed at the contacts of blocks and plates. We propose a phenomenological model to describe the initiation of stick-slip at the rough contact of blocks of rocks. The model contains the leading factors of stick-slip process (friction, roughness of contact surfaces, and asperity). The model reproduces the universal profile of slip velocity and displacement, velocity of dynamic rupture observed in experiments. The stick-slip motion at a rough surface of a fault is shown to be a nonlinear process and can be described by the generalized sine-Gordon equation. The results of calculations and a comparison with experimental data testify to the existence of waves of a new type — solitary waves of sliding. The analogies are found between the strain waves generated due to stick-slip at the contact of the blocks of rocks, strain waves in the crustal faults, and those within the lithosphere.     

Instability on a weakening fault

Weakening of a prestressed sawcut in Westerly granite under laboratory condition is accomplished by injecting pressurized fluid into the sawcut. After injection a sequence of stick-slips is observed while the deviatoric stress decreases successively with each stick-slip. On the basis of the experimental observation we develop a model of fault instability due to inhomogeneous and progressive weakening of the fault. According to this model, the fault surface is divided into the slipped and the locked regions, depending on whether or not the local state of stress satisfies the friction criterion. The average shear stress in the slipped region decreases with time and, in order to maintain a quasi-static equilibrium, shear stress in the remaining locked region on the fault surface increases. This situation may last until a critical state of stress on the fault is met, at which a sudden instability (stick0slip) may occur. We suggest that this mechanism of stress transfer may be a viable mechanism of induced seismicity and aftershocks, in addition to the well-known mechanism of a local increase of pore pressure. By comparing the experimental data with model predictions we show that the critical condition for slip instability is when the average shear stress over the locked region becomes equal to the value given by the friction criterion. Thus the friction criterion established for slip on fractures on which the state of stress is macroscopically uniform may also be applicable to fractures on which the stress state is macroscopically heterogeneous.     

利用超声尾波观测1.5 m长岩石断层黏滑实验的波速变化

在实验室内利用超声尾波观测大尺度(1.5 m)岩石断层的黏滑过程.利用基于尾波干涉的观测方法,我们获得了高达10~(-6)的相对波速变化的观测精度,这相当于~10 kPa的应力变化.利用高精度的测量,我们获取两种不同加载速率下(1μm·s~(-1),10μm·s~(-1))黏滑过程三个阶段(恢复、加载和滑动)基于波速变化的特征量.我们更进一步获取了断层失稳阶段波速变化的时空演化过程.最后讨论了该观测方法需要改进的地方.以上研究结果表明作为一种对现有实验观测手段的有益补充,利用超声尾波观测实验室大尺度岩石断层的动力学过程是可行的.     

岩石的双剪摩擦及震源时间函数的实验研究

本文介绍了岩石的双剪摩擦滑动实验。在实验宁测量了不同正应力条件下的粘滑应力降和位移值,求出了静、动摩擦系数;用瞬态波形存储器记录了粘滑过程中位移随时间的变化曲线,得出了震源时间函数,它的形状与布隆提出的震源时间函数极为相似,两者皆以指数形式下降,它的数值与野外实测结果大致相近。为验证实验数据的准确性,本文还对粘滑的理想模型进行了分析和计界。理论、实验和野外观测的结果大体一致。     

基于Chester-Higgs模型探讨摩擦生热对断层演化进程的影响

速率和状态相依赖的摩擦定律是本文采用的重要定律。结合Chester-Higgs摩擦模型和McKenzie-Brune摩擦生热模型,在一维弹簧-滑块-断层近似模型下,利用四阶变步长的Dormand-Prince算法,研究探讨了断层摩擦生热对断层演化的影响。结果表明:与忽略温度影响的情形相比,摩擦生热造成的温度上升可导致断层滑移时刻的略微提前,并伴随着摩擦系数和状态变量的下降,同时也使得断层的滑移量和应力降略有减小,而滑移速率有所增大;另外,在考虑温度影响时,有效正应力和临界滑移距离也会影响断层的演化过程,断层上的有效正应力越大,断层失稳时刻越提前,温度上升越明显;断层的临界滑移距离越大,断层失稳时刻则越迟,温度上升越显著,但当临界滑移距离超过5 cm时,具有不同临界滑移距离的断层,失稳时的温度则基本保持一致。      

Experimental Study on the Stick-slip Phenomenon of Intact Granite under Solid Confining Pressure

In this paper,compression tests of intact granite samples have been made in a triaxial testing machine with solid confining pressure.From the tests,the influences of confining pressure and loading rate(axial strain rate)on the deformation and fracture process of rock samples,on the stress drop and recurrence interval of stick-slip events,and on the geometric distribution of the main fracture have been studied.The experimental results show that the loading rate influences the stress drop and recurrence interval of stick-slip events greatly.At lower loading rates,the stress drop of stick-slip events is greater,their recurrence interval is longer and shows no regularity in distribution.When the loading rate goes higher,the stress drop will become smaller and the recurrence interval will tend to be constant and stick-slip events show a quasi-periodicity.At lower confining pressures and strain rates,the main fracture may evolve into 2 X-shaped conjugate shear faults; whereas at higher confining pressures and     

正断层黏滑错动作用下隧道抗错断技术模型试验研究

为提高穿越黏滑断层隧道的结构安全性和稳定性,开展了正断层黏滑错动作用下隧道抗错断措施作用效果室内模型试验。通过分析结构纵向应变、接触压力及安全系数,对未采取措施、结构加强、结构加强并施设减错层、结构加强并施设减错缝、结构加强并施设减错层和减错缝5种工况进行了对比分析。研究结果表明:采用减错措施（减错层/缝或减错层+减错缝）后,沿隧道纵向结构纵向应变、接触压力增加倍数由剧烈变化变为较均匀变化;仅采用结构加强措施抵抗正断层黏滑错动对隧道结构的影响,效果有限;通过采用减错措施消减正断层黏滑错动对隧道结构的影响,效果明显,结构安全系数最小值增加8倍以上;通过采用结构加强并施设减错层和减错缝措施抵抗和消减正断层黏滑错动对隧道结构的影响,效果显著,结构安全系数最小值增加25.45倍。研究成果可为高烈度艰险山区黏滑正断层隧道的抗错断设防设计提供参考。     

Effect of frictional heating and thermal advection on pre-seismic sliding: a numerical simulation using a rate-, state- and temperature-dependent friction law

Laboratory experiments on simulated faults in rocks clearly show the temperature dependence of dynamic rock friction. Since rocks surrounding faults are permeable, we have developed a numerical method to describe the thermo-mechanical evolution of the pre-seismic sliding phase which takes into account both the rate-, state- and temperature-dependent friction law and the heat advection term in the energy equation. We consider a laminar fluid motion perpendicular to a vertical fault plane and assume that fluids move away from the fault plane. A semi-analytical temperature solution which accounts for the variability of slip velocity and stress on the fault has been found. This solution has been generalized to the case of a time varying fluid velocity and then was used to include the thermal pressurization effect. After discretizing the temperature solution, the evolution of the system is obtained by the solution of a system of first order differential equations which allows us to determine the evolution of slip, slip rate, friction coefficient, effective normal stress, temperature and fluid velocity. The numerical solutions are found using a Runge-Kutta method with an adaptative stepsize control in time. When the thermal pressurization effects can be neglected, the heat advection effect gives rise to a delay, with respect to the purely conductive case, of the earthquake occurrence time. This delay increases with increasing permeability H of the system. When the thermal pressurization effects are taken into account the situation is opposite, i.e. the onset of instability tends to precede that of the purely conductive case. The advance in the time of occurrence of instability increases with increasing coefficient of thermal pressurization. In the small permeability range (H ≤ 10^?18 m²), the seismic moment and nucleation length of the pre-seismic phase are significantly smaller than those predicted by the purely conductive model.     

Seismicity and strain accumulation around Karliova Triple Junction (Turkey)

GPS studies in Turkey date back to the early 1990s, but were mostly focused on the seismically active North Anatolian Fault System (NAFS), or on the more populated Western Anatolia. Relatively few studies were made of the seismically less-active East Anatolian Fault System (EAFS), although it has the potential to produce large earthquakes. In this study, we present the results of a combination of geodetic and seismological data around the Karliova Triple Junction (KTJ), which lies at the intersection of the North- and East Anatolian Fault Systems. In particular, the geodetic slip rates obtained through block modeling of GPS velocities were compared with b-values to assess seismicity in the region. Yedisu segment, one of the best-known seismic gaps in Turkey, was specifically analyzed. The relatively low b-values across Yedisu segment verify the accumulation of seismic energy in this segment, and the GPS-derived geodetic slip rates suggest that it has the potential to produce an earthquake of Mw 7.5 across an 80-km rupture zone.Additionally, analysis of earthquake data reveals that the study area has a ductile or rigid–ductile behavior with respect to its surroundings, characterized by varying b-values. Although, seismic events of moderate- to high magnitudes are confined along the major fault zones, there are also low-seismicity zones along the eastern part of the Bitlis Suture Zone and around Yedisu. Since the high seismicity areas within the region may not accumulate sufficient stress for a large earthquake to occur, it is considered that the deformation in such areas occurs in a ductile manner. On the other hand, the areas characterized by low b-values may have the capacity of stress accumulation, which could lead to brittle deformation.     

2017年米林6.9级地震与1950年察隅8.6级地震的关系及两次地震对周边活动断层的影响

2017年11月18日,我国西藏林芝市米林县发生M6.9地震.其东南220 km距离处,1950年发生了察隅8.6级大地震.二者同处喜马拉雅东构造结大拐弯处,周边分布多条活动断裂带.察隅地震作为20世纪我国内陆最大的地震,它的发生对周边断层活动性的影响目前研究尚少.本文基于分层半无限空间黏弹性地球模型,计算了察隅地震对周围活动断裂带和米林地震的影响,同时利用有限断层破裂模型计算了米林地震对周边活动断层产生的同震库仑应力变化,并分析了两次地震的关系,讨论了下地壳及以下介质的黏滞系数和断层有效摩擦系数对计算结果的影响.结果表明,察隅地震影响范围广、强度大且持续时间长,对喜马拉雅东构造结周边活动断层均产生了较大的影响,库仑应力变化达到数MPa量级.在下地壳黏滞系数为1.0×10²⁰ Pa·s、摩擦系数为0.4情况下,按照青藏高原主要活动断裂每年由正常构造运动应力累积1~4 kPa计算,察隅地震的发生使米林地震提前了相当于2.79~11.15年的时间.米林6.9级地震影响范围和强度有限,只对震源附近的活动断层产生一定影响,库仑应力变化最大为数十kPa量级,对远处的断层影响较小.下地壳及以下介质的黏滞系数在震后黏滞松弛过程中影响逐渐明显,断层有效摩擦系数对同震库仑应力计算影响较大.