共查询到20条相似文献，搜索用时 31 毫秒
1.
Simulation of the Micro-physics of Rocks Using LSMearth 总被引：4，自引：0，他引：4
-- The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (M0ra and Place, 1994; Place and Mora, 1999). A new modular and flexible LSM approach has been developed that allows different micro-physics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions. 相似文献
2.
-- 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 Dc on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments. 相似文献
3.
Particle Dynamics Simulations of Rate- and State-dependent Frictional Sliding of Granular Fault Gouge 总被引：1，自引：0，他引：1
Julia K. Morgan 《Pure and Applied Geophysics》2004,161(9-10):1877-1891
— Recent simulations using the particle dynamics method (PDM) have successfully captured many features of natural faults zones as illuminated in laboratory studies. However, 2-D simulations conducted on idealized assemblages of particles using simple elastic-frictional contact laws, yield friction values considerably lower than natural materials, and lack time- and velocity-dependent changes in strength that influence dynamic fault slip. Here, preliminary results of new PDM simulations are described, in which particle motions are restricted as a proxy for particle interlocking and out of plane contacts, and time-dependent contact healing is introduced to capture temporal strengthening of granular assemblages. Frictional strength is increased, and in the absence of interparticle rolling, can attain values observed in the laboratory. The resulting mechanical behavior is qualitatively similar to that described by empirically-based rate-state friction laws, providing new physical insight into the discrete mechanics of natural faults. 相似文献
4.
Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model 总被引：4，自引：0，他引：4
— The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes. 相似文献
5.
The microstructural state and evolution of fault gouge has important implications for the mechanical behaviour, and hence
the seismic slip potential of faults. We use 3D discrete element (DEM) simulations to investigate the fragmentation processes
operating in fault gouge during shear. Our granular fault gouge models consist of aggregate grains, each composed of several
thousand spherical particles stuck together with breakable elastic bonds. The aggregate grains are confined between two blocks
of solid material and sheared under a given normal stress. During shear, the grains can fragment in a somewhat realistic way
leading to an evolution of grain size, grain shape and overall texture. The ‘breaking up’ of the fault gouge is driven by
two distinct comminution mechanisms: grain abrasion and grain splitting. The relative importance of the two mechanisms depends
on applied normal stress, boundary wall roughness and accumulated shear strain. If normal stress is sufficiently high, grain
splitting contributes significantly to comminution, particularly in the initial stages of the simulations. In contrast, grain
abrasion is the dominant mechanism operating in simulations carried out at lower normal stress and is also the main fragmentation
mechanism during the later stages of all simulations. Rough boundaries promote relatively more grain splitting whereas smooth
boundaries favor grain abrasion. Grain splitting (plus accompanying abrasion) appears to be an efficient mechanism for reducing
the mean grain size of the gouge debris and leads rapidly to a power law size distribution with an exponent that increases
with strain. Grain abrasion (acting alone) is an effective way to generate excess fine grains and leads to a bimodal distribution
of grain sizes. We suggest that these two distinct mechanisms would operate at different stages of a fault’s history. The
resulting distributions in grain size and grain shape may significantly affect frictional strength and stability. Our results
therefore have implications for the earthquake potential of seismically active faults with accumulations of gouge. They may
also be relevant to the susceptibility of rockslides since non-cohesive basal shear zones will evolve in a similar way and
potentially control the dynamics of the slide. 相似文献
6.
本文对龙门山断裂带金河磷矿浅钻岩芯中的三种断层泥开展了低速到高速摩擦滑动的实验研究,并对实验变形样品开展了BET比表面积研究.摩擦实验在干燥和孔隙水压条件下开展,速率范围涵盖20 μm·s-1~1.4 m·s-1.实验结果显示,三种断层泥在干燥条件下的摩擦性质差别不大,但在孔隙水压条件下,三者的中低速摩擦强度与层状硅酸盐矿物的种类而非总含量紧密相关,蒙脱石和伊利石相比绿泥石更能有效地弱化断层.三种断层泥在孔隙水压条件下存在中低速率域的速度强化,暗示着对断层的加速滑动存在一定的阻碍作用.孔隙水压下,黄绿色和灰绿色断层泥的初始动态弱化非常迅速并伴随断层泥层的瞬时扩容,凹凸体急剧加热导致的局部热压作用可能是造成这种力学行为的物理机制.在经历高速滑动之后,三种断层泥在干、湿条件下的BET比表面积都显著降低,暗示着可能发生了颗粒烧结.中低速域内,孔隙水的存在使得断层泥呈现分散式的剪切变形,BET比表面积的增加因此比干燥条件下更加明显.对表面能的估算表明,颗粒磨碎所消耗的能量至多不超过摩擦力做功的8%,暗示着断层作用中颗粒磨碎所占的能量比例较低. 相似文献
7.
断层岩的粒度分布包含岩石破裂机制、摩擦性质和地震能量分配等重要信息。筛分-称重和激光测量是分析断层岩三维粒度的2种有效方法,但每一种方法的测量范围仅有3个量级,难以全面反映断层岩的粒度分布特征。利用上述2种方法对汶川地震断层滑动带上的断层岩（简称断层岩）的粒度分布进行了测量,粒径测量范围从0.2μm至16mm,跨度达到5个数量级。结果显示：1）存在临界粒径dc（0.95～1.90μm）。粒度大于和小于dc的颗粒满足不同的颗粒数（Nd）-粒径（d）分布关系,表明该断层岩的粒径分布不具有自相似性特征。2）利用粒度大于dc的颗粒计算出的分形维数与断层岩类型有很好的相关性,即断层带边缘的角砾岩的平均分形维数为2.6,核部压碎角砾岩的平均分形维数约为3.0,中心断层泥的分形维数约为3.5。粒径小于dc的颗粒的分形维数为1.7～2.1。分形维数的突变反映出断层破裂机制的复杂性,即在不同的粒度域,岩石的破裂机制不尽相同。3）依据粒度分析结果,估算出汶川地震断层泥的单位破裂能（Es）为0.63MJ/m2。 相似文献
8.
In this study, 3-D Lattice Solid Model (LSMearth or LSM) was extended by introducing particle-scale rotation. In the new model,
for each 3-D particle, we introduce six degrees of freedom: Three for translational motion, and three for orientation. Six
kinds of relative motions are permitted between two neighboring particles, and six interactions are transferred, i.e., radial,
two shearing forces, twisting and two bending torques. By using quaternion algebra, relative rotation between two particles
is decomposed into two sequence-independent rotations such that all interactions due to the relative motions between interactive
rigid bodies can be uniquely decided. After incorporating this mechanism and introducing bond breaking under torsion and bending
into the LSM, several tests on 2-D and 3-D rock failure under uni-axial compression are carried out. Compared with the simulations
without the single particle rotational mechanism, the new simulation results match more closely experimental results of rock
fracture and hence, are encouraging. Since more parameters are introduced, an approach for choosing the new parameters is
presented. 相似文献
9.
The kinematics of gouge deformation 总被引：13，自引：0，他引：13
The methods and first results of a new approach to examining fault gouge are described. Samples of undisturbed fault gouge from the exhumed Lopez fault zone in the San Gabriel Mountains, California were impregnated with low viscosity epoxy resin and sectioned to produce microscope slides. The slides were photographed using optical and electron microscopy with magnifications ranging in factors of 2 from 12.5 to 1600. At all scales, the particles appeared angular with planar faces, suggesting tensile failure. No shear zones were discernable. The particle size distribution was studied. At each magnification the particles were sorted by diameter into four classes, differing in mean diameter by factors of 2. The numbers in each class were then scaled by the characteristic class dimension. The process revealed a remarkable degree of self-similarity. Over the range examined, the fractal dimension was within 5% of 2.60.On the basis of the observations, a new model for the mechanical processes that generate gouge is offered. It is argued that self-similarity results from repeated tensile splitting of grains. Unlike earlier models that consider splitting probability to be either independent of particle size or due to the preexisting distribution of defects, we propose that failure probability depends largely on the relative size of nearest neighbors. If nearest neighbors of the same size are preferentially broken, any initial distribution of particles will tend toward a self-similar distribution having a fractal dimension of 2.58.The model allows us to outline a procedure whereby the observed comminution in a fault zone can be related to the shear strain that the zone has accommodated and propose a theoretical frequency magnitude relation for the seismic energy emitted by the fracture process. 相似文献
10.
A systematical study is made on the characteristics of fault gouge from NW-trending faults in western Shandong.The macroscopic features of the fault gouge zones were investigated.From the data on surface structure of quartz fragments from the fault gouge and its maturity,the relative age,the movement mode,the rupture pattern,and the activities of the NW-trending faults are determined.Moreover,their seismogeological significance and some problems in study of fault gouge are discussed. 相似文献
11.
Internal structures and high-velocity frictional properties of a bedding-parallel carbonate fault at Xiaojiaqiao outcrop activated by the 2008 Wenchuan earthquake 
下载免费PDF全文
下载免费PDF全文 Linfeng Hou Shengli Ma Toshihiko Shimamoto Jianye Chen Lu Yao Xiaosong Yang Yuji Okimura 《地震科学（英文版）》2012,25(3):197-217
This paper reports internal structures of a bedding-parallel fault in Permian limestone at Xiaojiaqiao outcrop that was moved by about 0.5 m during the 2008 MW7.9 Wenchuan earthquake. The fault is located about 3 km to the south from the middle part of Yingxiu-Beichuan fault, a major fault in the Longmenshan fault system that was moved during the earthquake. The outcrop is also located at Anxian transfer zone between the northern and central segments of Yingxiu-Beichuan fault where fault system is complex. Thus the fault is an example of subsidiary faults activated by Wenchuan earthquake. The fault has a strike of 243° or N63°E and a dip of 38°NW and is nearly optimally oriented for thrust motion, in contrast to high-angle coseismic faults at most places. Surface outcrop and two shallow drilling studies reveal that the fault zone is several centimeters wide at most and that the coseismic slip zone during Wenchuan earthquake is about 1 mm thick. Fault zone contains foliated cataclasite, fault breccia, black gouge and yellowish gouge. Many clasts of foliated cataclasite and black gouge contained in fault breccia indicate multiple slip events along this fault. But fossils on both sides of fault do not indicate clear age difference and overall displacement along this fault should not be large. We also report results from high-velocity friction experiments conducted on yellowish gouge from the fault zone using a rotary shear low to high-velocity frictional testing apparatus. Dry experiments at normal stresses of 0.4 to 1.8 MPa and at slip rates of 0.08 to 1.35 m/s reveal dramatic slip weakening from the peak friction coefficient of around 0.6 to very low steady-state friction coefficient of 0.1-0.2. Slip weakening parameters of this carbonate fault zone are similar to those of clayey fault gouge from Yingxiu-Beichuan fault at Hongkou outcrop and from Pingxi fault zone. Our experimental result will provide a condition for triggering movement of subsidiary faults or off-fault damage during a large earthquake. 相似文献
12.
A Parallel Implementation of the Lattice Solid Model for the Simulation of Rock Mechanics and Earthquake Dynamics 总被引：2，自引：0，他引：2
— The Lattice Solid Model has been used successfully as a virtual laboratory to simulate fracturing of rocks, the dynamics of faults, earthquakes and gouge processes. However, results from those simulations show that in order to make the next step towards more realistic experiments it will be necessary to use models containing a significantly larger number of particles than current models. Thus, those simulations will require a greatly increased amount of computational resources. Whereas the computing power provided by single processors can be expected to increase according to “Moore’s law,” i.e., to double every 18–24 months, parallel computers can provide significantly larger computing power today. In order to make this computing power available for the simulation of the microphysics of earthquakes, a parallel version of the Lattice Solid Model has been implemented. Benchmarks using large models with several millions of particles have shown that the parallel implementation of the Lattice Solid Model can achieve a high parallel-efficiency of about 80% for large numbers of processors on different computer architectures. 相似文献
13.
断层泥的再生显微结构特征及其地震地质意义 总被引：3，自引：1，他引：3
三轴剪切摩擦实验后的断层泥与天然断层泥的再生显微结构特征研究表明，断层泥的显微结构特征与断层滑动方式之间有一定的关系，稳滑使断层泥变形均匀，产生低角度剪切（＜１４°）、布丁构造和颗粒碎裂流动。粘滑使断层泥局部发生强烈变形、高角度剪切（＞１４°）和碎粒出现随机裂纹等。断层泥的再生显微结构特征可用作鉴别古地震的存在 相似文献
14.
为更好地理解层状硅酸盐对断层强度、滑动速度依赖性及地震活动特征的影响,利用双轴摩擦实验对含白云母岩盐断层带在干燥及含水条件下摩擦的速度依赖性进行了实验研究,并观测了摩擦滑动过程中的声发射,分析了断层带的微观结构.实验结果表明,干燥条件下含白云母岩盐断层带在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个速度强化域对应不同的微观机制,高速域的速度强化受控于岩盐颗粒在白云母相互连结形成的网状结构上的滑动及其均匀碎裂作用,而低速域的速度强化还受岩盐的压溶作用控制.通过与岩盐断层带摩擦实验结果对比可知,白云母的存在对于燥岩盐断层带摩擦滑动方式和速度依赖性没有显著影响,而在含水条件下白云母的存在使得岩盐断层带滑动趋于稳定.实验结果为分析含层状硅酸盐断层的强度和稳定性提供了依据.此外,在速度依赖性转换域上观察到的应力缓慢释放的现象进一步证实了在岩盐断层带摩擦滑动过程中观察到的现象,这对慢地震机制研究具有参考意义. 相似文献
15.
Internal structures and high-velocity frictional properties of Longmenshan fault zone at Shenxigou activated during the 2008 Wenchuan earthquake 下载免费PDF全文
下载免费PDF全文 Yu Wang Shengli Ma Toshihiko Shimamoto Lu Yao Jianye Chen Xiaosong Yang Honglin He Jiaxiang Dang Linfeng Hou Tetsuhiro Togo 《地震科学（英文版）》2014,27(5):499-528
This paper reports internal structures of a wide fault zone at Shenxigou, Dujiangyan, Sichuan province, China, and high-velocity frictional properties of the fault gouge collected near the coseismic slip zone during the 2008 Wenchuan earthquake. Vertical offset and horizontal displacement at the trench site were 2.8 m (NW side up) and 4.8 m (right-lateral), respectively. The fault zone formed in Triassic sandstone, siltstone, and shale about 500 m away from the Yingxiu-Beichuan fault, a major fault in the Longmenshan fault system. A trench survey across the coseismic fault, and observations of outcrops and drill cores down to a depth of 57 m revealed that the fault zone consists of fault gouge and fault breccia of about 0.5 and 250–300 m in widths, respectively, and that the fault strikes N62°E and dips 68° to NW. Quaternary conglomerates were recovered beneath the fault in the drilling, so that the fault moved at least 55 m along the coseismic slip zone, experiencing about 18 events of similar sizes. The fault core is composed of grayish gouge (GG) and blackish gouge (BG) with very complex slip-zone structures. BG contains low-crystalline graphite of about 30 %. High-velocity friction experiments were conducted at normal stresses of 0.6–2.1 MPa and slip rates of 0.1–2.1 m/s. Both GG and BG exhibit dramatic slip weakening at constant high slip rates that can be described as an exponential decay from peak friction coefficient μ p to steady-state friction coefficient μ ss over a slip-weakening distance D c. Deformation of GG and BG is characterized by overlapped slip-zone structures and development of sharp slickenside surfaces, respectively. Comparison of our data with those reported for other outcrops indicates that the high-velocity frictional properties of the Longmenshan fault zones are quite uniform and the high-velocity weakening must have promoted dynamic rupture propagation during the Wenchuan earthquake. 相似文献
16.
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. 相似文献
17.
Fractal research of fault gouge 总被引：1，自引：0，他引：1
ＦｒａｃｔａｌｒｅｓｅａｒｃｈｏｆｆａｕｌｔｇｏｕｇｅＳＨＵＮ－ＭＥＩＳＨＡＯ（邵顺妹）ａｎｄＪＩＮ－ＣＨＡＮＧＺＯＵ（邹瑾敞）ＥａｒｔｈｑｕａｋｅＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅｏｆＬａｎｚｈｏｕ，ＳｔａｔｅＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒｅａ... 相似文献
18.
19.
R. M. Stesky 《Pure and Applied Geophysics》1978,116(4-5):690-704
This paper reviews many of the mechanical properties of faulted and jointed rock under pressure and temperature and in the presence of water. At low effective confining pressures (below about 1 kilobar), the friction strength is quite variable and depends on the frictional resistance between gouge particles or asperities and on the dilatancy of the fault. At higher pressures the friction strength is nearly independent of mineralogy, temperature, and rate, at least for rocks whose friction strength is less than the failure strength. Water tends to slightly weaken the fault. The type of sliding motion, whether stick-slip or stable sliding, is much more affected by environmental and mineralogical factors. In general, stick-slip is dominant at high pressures and low temperatures, in the presence of strong minerals such as quartz and feldspar, in the absence of gouge, for lower surface roughness, and perhaps in the presence of water. The microscopic deformation mechanisms are poorly understood. At low temperatures, cataclasis dominates in rocks containing mostly quartz or feldspar, and plastic deformation in rocks containing mostly calcite or platy silicates. At high temperature most minerals deform plastically, producing a greater temperature-and rate-dependence of the friction strength. Glass has been found in some sliding surfaces in sandstone. 相似文献
20.
基于速率与状态依赖性摩擦本构关系理论框架,在热水条件下研究了角闪石断层泥的摩擦滑动性质并与闪长岩的另一种主要矿物斜长石的摩擦滑动性质进行了对比.摩擦实验是在三轴实验系统上完成,有效正应力200 MPa,孔隙压力30 MPa,并将加载速率在1.22 μm/s和 0.122 μm/s之间实施了切换.结果表明角闪石的摩擦系数均值为0.70±0.01,随着温度增加没有系统性的变化,整体低于斜长石的摩擦系数(0.75±0.01);计算与实验表明,角闪石和斜长石的摩擦系数的体积分数加权平均值与闪长岩的摩擦系数基本一致;角闪石在实验温度范围内(100～614 ℃)显示速率强化(a-b>0),与斜长石在整个温度范围内的速率弱化(a-b<0)正好相反;角闪石的速率依赖性在整个实验温度范围内无系统性的变化. 相似文献