共查询到20条相似文献，搜索用时 62 毫秒
1.
Nucleation of earthquakes and its implication to precursors 总被引：1，自引：0，他引：1
IntroductionIn the middle of 1990s, the argument about nucleation phase occurred. Ellsworth and Beroza(1 995) studied the wave records of 30 earthquake events, which initiated slowly. Ellsworth, et al calledthem as nucleation phase. They discussed the relationship between their moment, duration time andthe relationship with the main shock behind. Scherbaum (1997a, b) proposed his different view andsuggested that some of the recording of the long period initial of seismic wave was misinterpret… 相似文献
2.
A two-stage model of fracture of rocks 总被引：7，自引：0，他引：7
In this paper we propose a two-stage model of rock fracture. In the first stage, cracks or local regions of failure are uncorrelated and occur randomly throughout the rock in response to loading of pre-existing flaws. As damage accumulates in the rock, there is a gradual increase in the probability that large clusters of closely spaced cracks or local failure sites will develop. Based on statistical arguments, a critical density of damage will occur where clusters of flaws become large enough to lead to larger-scale failure of the rock (stage two). While crack interaction and cooperative failure is expected to occur within clusters of closely spaced cracks, the initial development of clusters is predicted based on the random variation in pre-existing flaw populations. Thus the onset of the unstable second stage in the model can be computed from the generation of random, uncorrelated damage. The proposed model incorporates notions of the kinetic (and therefore time-dependent) nature of the strength of solids as well as the discrete hierarchic structure of rocks and the flaw populations that lead to damage accumulation. The advantage offered by this model is that its salient features are valid for fracture processes occurring over a wide range of scales including earthquake processes. A notion of the rank of fracture (fracture size) is introduced, and criteria are presented for both fracture nucleation and the transition of the failure process from one scale to another. 相似文献
3.
Mitiyasu Ohnaka 《Pure and Applied Geophysics》2004,161(9-10):1915-1929
It is shown that the rupture nucleation length increases up to the critical length with time according to a power law, and that the accelerating phase of nucleation leading up to the critical point is scaled in the framework of fracture mechanics based on slip-dependent constitutive formulation. Geometric irregularity of the rupturing surfaces plays a fundamental role in scaling the accelerating phase of nucleation up to the critical point. A power-law scaling relation between the rupture growth length and the nucleation time to the critical point is derived from theoretical consideration based on laboratory data. This power-law scaling relation has no singularity, and hence it may be useful for the predictive purpose of an imminent, large earthquake. 相似文献
4.
-- We have simulated a rupture transition from quasi-static growth to dynamic propagation using the boundary integral equation method. In order to make a physically reasonable model of earthquake cycle, we have to evaluate the dynamic rupture propagation in the context of quasi-static simulation. We used a snapshot of the stress distribution just before the earthquake in the quasi-static simulation. The resultant stress will be fed back to the quasi-static simulation. Since the quasi-static simulation used the slip-and time-dependent constitutive relation, the friction law itself evolves with time. Thus, we used the slip-weakening constitutive relation for dynamic rupture propagation consistent with that used for the quasi-static simulation. We modeled a San Andreas type strike-slip fault, in which two different size asperities existed. 相似文献
5.
James R. Rice 《Pure and Applied Geophysics》1983,121(3):443-475
Constitutive relations for fault slip are described and adopted as a basis for analyzing slip motion and its instability in
the form of earthquakes on crustal faults. The constitutive relations discussed include simple rate-independent slip-weakening
models, in which shear strength degrades with ongoing slip to a residual frictional strength, and also more realistic but
as yet less extensively applied slip-rate and surface-state-dependent relations. For the latter the state of the surface is
characterized by one or more variables that evolve with ongoing slip, seeking values consistent with the current slip rate.
Models of crustal faults range from simple, single-degree-of-freedom spring-slider systems to more complex continuous systems
that incorporate nonuniform slip and locked patches on faults of depth-dependent constitutive properties within elastic lithospheric
plates that may be coupled to a viscoelastic asthenosphere.
Most progress for the rate and state-dependent constitutive relations is at present limited to single-degree-of-freedom systems.
Results for stable and unstable slip with the various constitutive models are summarized. Instability conditions are compared
for spatially uniform versus nonuniform slip, including the elastic — brittle crack limit of the nonuniform mode. Inferences
of constitutive and fracture parameters are discussed, based on earthquake data for large ruptures that begin with slip at
depth, concentrating stress on locked regions within a brittle upper crust. Results of nonlinear stability theory, including
regimes of complex sustained stress and slip rate oscillations, are outlined for rate and state-dependent constitutive relations,
and the manner in which these allow phenomena like time-dependent failure, restrengthening in nearly stationary contact, and
weakening in rapidly accelerated slip, is discussed. 相似文献
6.
We numerically study the dynamic interaction of propagating cracks. It is assumed that propagating cracks can nucleate and drive subsidiary cracks because of shear strain enhancement near the propagating crack tips. The critical strain fracture criterion is assumed in the analysis. Intense interaction is expected to occur among the cracks. All the cracks are assumed to be parallel and antiplane strain deformation is assumed in the computation.In the interaction of two non-coplanar cracks, a strain shadow is formed in the neighborhood of each crack because of the strain release by the introduction of the crack. The growth of each crack is accelerated when the propagating tips of each crack are outside of the strain shadow of the other crack. In general, the crack tips enter the strain shadow, and the crack tips decelerate. The calculation shows that only one of the two cracks can continue to grow, and the other's growth is decelerated and arrested. If we can assume that the suite of cracks interact in a pairwise manner only, then this may suggest that only a limited number of cracks can continue to grow during the final stage of the rupture process. Hence the crack interaction causes complexity in dynamic earthquake faulting. The concepts of barrier and asperity have been employed by many researchers for the interpretation of complex seismic wave data. However, the physical realities of such concepts are obscure. Our calculations show that dynamic crack interactions can produce barriers and asperities in some cases; the crack tip deceleration or arrest due to the interactions among non-coplanar cracks can be interpreted as being due to a barrier. The dynamic coalescence among the coplanar cracks can be regarded as an asperity.Umeda found a localized area that strongly radiates high-frequency seismic waves in the epicentral areas of some large shallow earthquakes. He defined this as an earthquake bright spot. Our analysis implies that only a limited number of cracks continue to grow when many interactive cracks nucleate, and that all other cracks stop extending soon after nucleation. Hence, if the nucleation and termination of several cracks occur in a localized area, it will be observed seismologically as an earthquake bright spot. This is because it is theoretically known that the sudden termination of crack growth and dynamic crack coalescence efficiently emits high-frequency elastic waves. 相似文献
7.
在以前的工作中,考虑直立走滑型断层地震,假设断层面微元破裂强度遵循Weibull概率分布,由细观力学方法推导出断层面的宏观本构关系是一个非线性函数,表现为弹性-软化塑性特征,在此基础上用稳定性理论研究了地震稳定性问题.而实际断层大多是倾斜的,为此,本文首先建立了由围岩和倾斜断层构成的平面地震力学模型,采用宏观的断层载荷-变形的全过程曲线,详细讨论了倾斜断层地震的不稳定性问题.结果表明,远场一旦施加位移,断层也同时错动,这可能与实际情况不符合.为了更好的模拟断层的初始能量累计过程,进一步对断层本构模型进行改进.考虑断层面破裂强度,采用Coulomb破裂准则,则断层表现为刚塑性本构关系,只有当断层面剪应力达到一个临界值时,断层才开始错动.研究表明,对于倾斜断层地震,与直立走滑型断层地震一样,系统刚度比β（围岩切线刚度与断层刚度最大值之比）是决定地震失稳的重要参数,只有当β<1时才会出现地震失稳,且伴随应力突跳和围岩应变能释放.当β≥1时,仅仅是断层无震滑动,不会发生地震.在远场应施以位移形式边界条件,以致地震失稳发生在平衡路径的位移转向点并伴有应力突跳. 相似文献
8.
天然地震发生后,地震波及区域内的地下岩层渗透率常常会发生显著改变,其变化曲线显示出独有的特征,造成这一现象的机理较为复杂,传统渗流理论尚不能给出合理解释.针对这一问题,从震后渗透率变化规律入手,深入分析了地下岩层裂缝体系对渗透率的影响,给出了裂缝结构参数与渗透率之间的定量关系.结合岩层黏弹特性以及天然地震所产生的地下岩层体应变特征,基于裂缝体系分维度正比于外部应力的实验事实,将黏弹体应力松弛机制引入该体系,对裂缝分形渗透率模型进行了含时推广,建立起震后地下岩层渗透率的时间演化模型,理论预测曲线与实验曲线吻合较好.在此基础上提出‘分形裂缝渗透率松弛效应’这一全新概念.本研究为震控流体运移研究提供了新思路,对于揭示震后断层恢复机制,探讨断层活动与孕震的关联有一定的理论价值和现实意义. 相似文献
9.
简述了最近20年来国内外岩石高速摩擦实验研究领域的进展和动态:岩石高速摩擦实验技术的发展实现了对高滑动速率、大位移的地震过程的实验模拟;其结果揭示了岩石和断层泥在地震滑动速率下的力学性状,深化了对断层滑动弱化机制、临界滑动距离、以及地震发生过程的认识和理解;实验在假玄武玻璃成因方面取得了重要进展,并提出了断层发生地震滑动可能留下的其它地质证据,可望为研究断层滑动性状与地震物理过程提供新的思路和信息.岩石高速摩擦实验今后的发展方向主要包括:发展具有加温系统和孔隙压系统的岩石高速摩擦实验装置,研究水热作用下岩石和断层泥的高速摩擦性状;室内实验和地震资料分析相结合研究断层滑动和地震机制;室内实验和野外地质调查相结合探索断层发生地震错动的地质证据等等. 相似文献
10.
基于断层面强度的非均匀性,将断层面的宏观破裂过程看作是断面局部微元的破裂累积过程,假设断层局部微元强度遵循Weibull概率分布,从统计力学角度推导出了宏观的断层载荷.变形的全过程本构关系.采用一维地震力学模型,以远场位移a为控制变量,用稳定性理论研究了地震不稳定性问题.研究表明,系统刚度比(围岩刚度与断层刚度之比)β是影响地震发生的重要参数,只有当β<1时才会出现地震失稳,且应力突跳发生在平衡路径曲线的位移转向点,并给出了地震过程的3个重要参数(地震后断层半错距、地震应力降和释放的弹性能)的表达式.当β≥1时,不会发生地震,仅是缓慢的断层滑动,属于无震滑动. 相似文献
11.
12.
利用FLAC,模拟了双轴压缩岩样的破坏过程及剪应力异常.在峰值强度之前及之后,岩石的本构模型分别取为线弹性及莫尔库仑剪破坏与拉破坏复合的应变软化模型.本文仅分析了3个具有较大应力降的周期的剪应力异常及塑性区分布.在第1及第2周期,剪破坏分别发生在试样的两侧及背景空区内部,前兆明显,地震属于断错型.在第3周期,背景空区之外的新剪切破裂带引起断错型地震;空区内部的老破裂带引起走滑型地震,未观测到前兆.在塑性区边缘,剪应力梯度较高.破裂带位置剪应力值反而较低.尽管在应力-应变曲线的软化阶段之前,不同单元的剪应力表现出千姿百态的复杂形态,但尚有规律可循.若几个单元都位于某条破裂带上,则剪应力的变化可能是同步的,或演变规律类似.鉴于原地复发地震类型可能不同,前兆将有差别. 相似文献
13.
W. Müller 《Pure and Applied Geophysics》1986,124(4-5):693-709
Fracture phenomena in rocks are associated with mainly mode I crack growth, sometimes superposed by shear or torsion. The present paper contributes to a fracture mechanics analysis of mode I and mixed mode crack propagation, by presenting reliable fracture toughness data for some rocks which include the effect of induced crack propagation rate, and the influence of effective pressure, and by numerical calculations on fracture propagation in layered rock formations. Empirical relations between fracture toughness,K
Ic' and induced crack opening displacement rate, as well as effective pressure, are given. The observedK
Ic pressure relation supports a theoretical model which takes into account the existence of microcracks in the crack tip region. Finite element calculations of fracture propagation in layered rock formations demonstrate the important effect of mixed mode crack growth. The numerical approach is particularly applied to single crack growth in hydraulic fracturing and in three point bending tests on layered single edge crack specimens. 相似文献
14.
A conceptual model for describing effective saturation in fractured hard rock is presented. The fracture network and the rock matrix are considered as an equivalent continuum medium where each fracture is conceptualized as a porous medium of granular structure and the rock matrix is assumed to be impermeable. The proposed model is based on the representation of a rough‐walled fracture by an equivalent porous medium, which is described using classical constitutive models. A simple closed‐form equation for the effective saturation is obtained when the van Genuchten model is used to describe saturation inside fractures and fractal laws are assumed for both aperture and number of fractures. The relative hydraulic conductivity for the fractured rock is predicted from a simple relation derived by Liu and Bodvarsson. The proposed constitutive model contains three independent parameters, which may be obtained by fitting the proposed effective saturation curve to experimental data. Two of the model parameters have physical meaning and can be identified with the reciprocal of the air entry pressure values in the fractures of minimum and maximum apertures. Effective saturation and relative hydraulic conductivity curves match fairly well the simulated constitutive relations obtained by Liu and Bodvarsson. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
15.
—Forced torsional oscillation techniques have been used to explore the seismic-frequency shear mode viscoelasticity of specimens of two crustal rocks (Cape Sorell quartzite and Delegate aplite), cycled between room temperature and 700°C under conditions of moderate confining pressure. The anisotropy and intergranular inhomogeneity of thermal expansivity in these materials give rise to large deviatoric stresses, resulting in thermal cracking at temperatures above a pressure-dependent threshold temperature, associated with the onset of very pronounced temperature sensitivity of the shear modulus, in general accord with the predictions of fracture mechanics models. For Delegate aplite in particular, the shear modulus behaves reproducibly during multiple thermal cycles at different confining pressures, consistent with the notion that the thermal cracks are of low aspect ratio (minimum/maximum dimension), and are therefore readily closed by the prevailing confining pressure once the thermal stresses are removed. Marked frequency-dependent dissipation of shear strain energy is observed on heating each rock to temperatures ≥ 500°C, although the attenuation varies significantly with prior thermal history, probably as a result of progressive dehydration and relaxation of deviatoric stresses. Temperature and pressure dependent crack densities for Delegate aplite have been estimated by comparison of the observed shear moduli with those expected for a crack-free aggregate. In parallel with the forced oscillation tests, measurements have been made of the rate at which (argon) pore pressure equilibrium is re-established following a perturbation. Combination of these results, which provide a proxy for permeability, with the inferred crack densities indicates that the variation of permeability with crack density is well described by a percolation model with a threshold crack density of ～0.2. 相似文献
16.
Chen Yongming 《中国地震研究》2004,18(4):337-347
The investigation on damages to frozen soil sites during the West Kunlun Mountains Pass earthquake with Ms 8.1 in 2001 shows that the frozen soil in the seismic area is composed mainly of moraine, alluvial deposit, diluvial deposit and lacustrine deposit with the depth varying greatly along the earthquake rupture zone. The deformation and rupture of frozen soil sites are mainly in the form of coseismic fracture zones caused by tectonic motion and fissures,liquefaction, seismic subsidence and collapse resulting from ground motion. The earthquake fracture zones on the surface are main brittle deformations, which, under the effect of sinlstral strike-slip movement, are represented by shear fissures, tensional cracks and compressive bulges. The distribution and configuration patterns of deformation and rupture such as fissures, liquefaction, seismic subsidence and landslides are all related to the ambient rock and soil conditions of the earthquake area. The distribution of earthquake damage is characterized by large-scale rupture zones, rapid intensity attenuation along the Qinghai-Xizang (Tibet) Highway, where buildings distribute and predominant effect of rock and soil conditions. 相似文献
17.
An understanding of the frictional sliding on faults that can lead to earthquakes requires a knowledge of both constitutive behavior of the sliding surfaces and its mechanical interaction with the loading system. We have determined the constitutive parameters for frictional sliding of initially bare surfaces of Westerly granite, using a recently developed high pressure rotary shear apparatus that allows long distances of sliding and therefore a greater assurance of attaining steady state behavior. From experiments conducted at room temperature and normal stresses of 27–84 MPa several important results have been found. (1) A gouge layer 100 to 200 m thick was developed from the initially bare rock surfaces after 18 to 70 mm of sliding. (2) The steady state frictional resistance, attained after about 10 mm of sliding, is proportional to the negative of the logarithm of the sliding velocity. (3) Abrup changes in the velocity of sliding result in initial changes in the frictional resistance, which have the same sign as the velocity change, and are followed by a gradual decay to a new steady state value over a characteristic distance of sliding. This velocity weakening behavior is essentially identical with that found by several previous workers on the same material at lower normal stress. (4) Our results are well described by a two state variable constitutive law. The values of the constitutive parameters are quite similar to those found previously at low normal stress, but the characteristic distance is about an order of magnitude smaller than that found at 10 MPa normal stress with thicker layers of coarser gouge. (5) We have approximated our results with a one state variable constitutive law and compared the results with the predictions of existing nonlinear stability analysis; in addition, we have extended the stability analysis to systems possessing two state variables. With such formulations good agreement is found between the experimentally observed and theoretically predicted transitions between stable and unstable sliding. These results allow a better understanding of the instabilities that lead to earthquakes. 相似文献
18.
The development of a powerful numerical model to simulate the fracture behavior of concrete material has long been one of the dominant research areas in earthquake engineering.A reliable model should be able to adequately represent the discontinuous characteristics of cracks and simulate various failure behaviors under complicated loading conditions.In this paper,a numerical formulation,which incorporates a sophisticated rigid-plastic interface constitutive model coupling cohesion softening,contact,friction and shear dilatation into the XFEM,is proposed to describe various crack behaviors of concrete material.An effective numerical integration scheme for accurately assembling the contribution to the weak form on both sides of the discontinuity is introduced.The effectiveness of the proposed method has been assessed by simulating several well-known experimental tests.It is concluded that the numerical method can successfully capture the crack paths and accurately predict the fracture behavior of concrete structures.The influence of mode-Ⅱ parameters on the mixed-mode fracture behavior is further investigated to better determine these parameters. 相似文献
19.