首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 328 毫秒
1.
We present a theory for the radiation of high-frequency waves by earthquake faults. We model the fault as a planar region in which the stress drops to the kinematic friction during slip. This model is entirely equivalent to a shear crack. For two-dimensional fault models we show that the high frequencies originate from the stress and slip velocity concentrations in the vicinity of the fault's edges. These stress concentrations radiate when the crack expands with accelerated motion. The most efficient generation of high-frequency waves occurs when the rupture velocity changes abruptly. In this case, the displacement spectrum has an ω-2 behaviour at high frequencies. The excitation is proportional to the intensity of the stress concentration near the crack tips and to the change in the focusing factor due to rupture velocity. We extend these two-dimensional results to more general three-dimensional fault models in the case when the rupture velocity changes simultaneously on the rupture front. Results are similar to those described for two-dimensional faults. We apply the theory to the case of a circular fault that grows at constant velocity and stops suddenly. The present theory is in excellent agreement with a numerical solution of the same problem.
Our results provide upper bounds to the high-frequency radiation from more realistic models in which rupture velocity does not change suddenly. The ω-2 is the minimum possible decay at high frequencies for any crack model of the source.  相似文献   

2.
Summary. A possible mechanism for the occurrence of slow earthquakes is investigated by calculating numerical solutions for the dynamical rupture process on a quasi-three-dimensional fault with heterogeneous frictional strengths. Experimental friction laws for the dependence of sliding frictional stress on slip velocity, which are based on the cohesive properties of fault asperities, are taken into considerations.
It is found that the applied stress does not drop very rapidly with time and the rupture velocity remarkably decreases as the dependence on slip-velocity becomes smaller. These deceleration effects for the rupture propagation are greatly enhanced with increasing heterogeneities in the distribution of frictional strength and as the initial shear stress has lower levels with respect to the average strength. For these cases, the growth of rupture is extremely slow in a nucleus region with the dimension as large as 10 times the initial rupture length, and gains a terminal velocity dependent on the above factors. The displacement-time function becomes noticeably extended in these cases, and indicates a stick—slip-like phenomena in the extended time interval for a strongly heterogeneous fault.
It seems that these results could explain the characteristic features of slow earthquakes.  相似文献   

3.
Summary. We consider the problem of the bilateral extension of a two-dimensional anti-plane shear crack that initiates spontaneously at a point and extends under the influence of cohesive forces at the edges. An approximation to the stresses in the regions beyond the edges of the crack has been found that simplifies the calculation. The exact stresses in these regions are also found iteratively. In the cases of uniformly propagating cracks, the estimates of cohesive forces obtained from the approximation are close to the exact values for high crack speeds but are significantly different for low crack speeds. It is also found that if healing is initiated due to the encounter of one end of a uniformly propagating crack with an unbreakable barrier, the static stress drop in the torn region is constant but may either overshoot or undershoot the dynamical stress drop. In these cases, the final static slip distribution is obtained by freezing the dynamic solution along a characteristic line through the location of the barrier. We find that the crack length cannot be unambiguously derived from the far field spectral properties.  相似文献   

4.
In order to better understand the development of thrust fault‐related folds, a 3D forward numerical model has been developed to investigate the effects that lateral slip distribution and propagation rate have on the fold geometry of pre‐ and syn‐tectonic strata. We consider a fault‐propagation fold in which the fault propagates upwards from a basal decollement and along‐strike normal to transport direction. Over a 1 Ma runtime, the fault reaches a maximum length of 10 km and accumulates a maximum displacement of 1 km. Deformation ahead of the propagating fault tip is modelled using trishear kinematics while backlimb deformation is modelled using kink‐band migration. The applicability of two different lateral slip distributions, namely linear‐taper and block‐taper, are firstly tested using a constant lateral propagation rate. A block‐taper slip distribution replicates the geometry of natural fold‐thrusts better and is then used to test the sensitivity of thrust‐fold morphology to varied propagation rates in a set of fault‐propagation folds that have identical final displacement to length (Dmax/Lmax) ratios. Two stratigraphic settings are considered: a model in which background sedimentation rates are high and no topography develops, and a model in which a topographic high develops above the growing fold and local erosion, transport and deposition occur. If the lateral propagation rate is rapid (or geologically instantaneous), the fault tips quickly become pinned as the fault reaches its maximum lateral extent (10 km), after which displacement accumulates. In both stratigraphic settings, this leads to strike‐parallel rotation of the syn‐tectonic strata near the fault tips; high sedimentation rates relative to rates of uplift result in along‐strike thinning over the structural high, while low sedimentation rates result in pinchout against it. In contrast, slower lateral propagation rates (i.e. up to one order of magnitude greater than slip rate) lead to the development of along‐strike growth triangles when sedimentation rates are high, whereas when sedimentation rates are low, offflap geometries result. Overall we find that the most rapid lateral propagation rates produce the most realistic geometries. In both settings, time‐equivalent units display both nongrowth and growth stratal geometries along‐strike and the transition from growth to nongrowth has the potential to delineate the time of fault/fold growth at a given location. This work highlights the importance of lateral fault‐propagation and fault tip pinning on fault and fold growth in three dimensions and the complex syn‐tectonic geometries that can result.  相似文献   

5.
A theoretical approach to the propagation of interacting cracks   总被引：2，自引：0，他引：2  
We propose a scheme to compute interaction effects between two randomly oriented cracks under compressive stresses and we discuss the role crack interactions play in the crack coalescence process. Stress intensity factors are computed by using an iterative technique based on the method of successive approximations. Once crack propagation occurs, curved wing cracks grow from the initial crack tips. The stress intensity factors at the wing crack tips are calculated as the sum of two terms: a component for a single wing crack subjected to both the applied stresses and the interaction effect, and a component due to the sliding of the initial crack. We have applied our procedure to various crack geometries. Our results show that interaction effects act on the crack propagation path. For cracks under tension, our approach correctly predicts the curving, hook-shaped paths of interacting cracks that have been observed in various materials. For en echelon compressive cracks, interaction effects depend on the geometry of stepping. For right-stepping cracks, no mode I crack coalescence occurs. A mixedmode propagation criterion may be introduced to check whether coalescing secondary shear fractures initiate. For left-stepping cracks, depending on whether or not there is overlapping, crack coalescence is achieved by tension wing cracks at the inner crack tips. Without overlapping, the growing wing cracks delimit a region where a tensile secondary fracture may develop and lead to coalescence. These results are consistent with previous work and show that our procedure may be now extended to a population of cracks.  相似文献   

6.
The dynamic coalescence of two mode II cracks on a planar fault is simulated here using the elastodynamic boundary integral equation method. We focus on the complexity of the resultant slip rate and seismic radiation in the crack coalescence model (CCM) and on the reconstruction of a single crack model (SCM) that can reproduce the CCM waveforms from heterogeneous source parameters rather than coalescence. Simulation results reveal that localized higher slip rates are generated by coalescence as a result of stress interaction between the approaching crack tips. The synthesized seismic radiation exhibits a distinct coalescence phase that has striking similarities to stopping phases in the radiation and propagation properties. The corresponding SCM yields a singular increase in the stress drop distribution, which is accompanied by a sudden decrease in it across the point of coalescence in the CCM. This implies that the generation of high-frequency radiation is more efficient from coalescence than from stopping, although both phenomena exhibit the same strong  ω−2  -type displacement spectra.  相似文献   

7.
Source models such as the k -squared stochastic source model with k -dependent rise time are able to reproduce source complexity commonly observed in earthquake slip inversions. An analysis of the dynamic stress field associated with the slip history prescribed in these kinematic models can indicate possible inconsistencies with physics of faulting. The static stress drop, the strength excess, the breakdown stress drop and critical slip weakening distance D c distributions are determined in this study for the kinematic k -squared source model with k -dependent rise time. Several studied k -squared models are found to be consistent with the slip weakening friction law along a substantial part of the fault. A new quantity, the stress delay, is introduced to map areas where the yielding criterion of the slip weakening friction is violated. Hisada's slip velocity function is found to be more consistent with the source dynamics than Boxcar, Brune's and Dirac's slip velocity functions. Constant rupture velocities close to the Rayleigh velocity are inconsistent with the k -squared model, because they break the yielding criterion of the slip weakening friction law. The bimodal character of D c / D tot frequency–magnitude distribution was found. D c approaches the final slip D tot near the edge of both the fault and asperity. We emphasize that both filtering and smoothing routinely applied in slip inversions may have a strong effect on the space–time pattern of the inferred stress field, leading potentially to an oversimplified view of earthquake source dynamics.  相似文献   

8.
Modelling dynamic rupture for complex geometrical fault structures is performed through a finite volume method. After transformations for building up the partial differential system following explicit conservative law, we design an unstructured bi-dimensional time-domain numerical formulation of the crack problem. As a result, arbitrary non-planar faults can be explicitly represented without extra computational cost. On these complex surfaces, boundary conditions are set on stress fluxes and not on stress values. Prescribed rupture velocity gives accurate solutions with respect to analytical ones depending on the mesh refinement, while solutions for spontaneous propagation are analysed through numerical means. An example of non-planar spontaneous fault growth in heterogeneous media demonstrates the good behaviour of the proposed algorithm as well as specific difficulties of such numerical modelling.  相似文献   

9.
We perform analytical and numerical studies of scaling relations of earthquakes and partition of elastic strain energy between seismic and aseismic components using a thermodynamically based continuum damage model. Brittle instabilities occur in the model at critical damage level associated with loss of convexity of the strain energy function. A new procedure is developed for calculating stress drop and plastic strain in regions sustaining brittle instabilities. The formulation connects the damage rheology parameters with dynamic friction of simpler frameworks, and the plastic strain accumulation is governed by a procedure that is equivalent to Drucker–Prager plasticity. The numerical simulations use variable boundary forces proportional to the slip-deficit between the assumed far field plate motion and displacement of the boundary nodes. These boundary conditions account for the evolution of elastic properties and plastic strain in the model region. 3-D simulations of earthquakes in a model with a large strike-slip fault produce scaling relations between the scalar seismic potency, rupture area, and stress drop values that are in good agreement with observations and other theoretical studies. The area and potency of the simulated earthquakes generally follow a linear log–log relation with a slope of 2/3, and are associated with stress drop values between 1 and 10 MPa. A parameter-space study shows that the area-potency scaling is shifted to higher stress drops in simulations with parameters corresponding to lower dynamic friction, more efficient healing, and higher degree of seismic coupling.  相似文献   

10.
A crack model in antiplane shear configuration is shown representing creep processes interpreted in terms of 'viscous' deformation of a narrow plastic layer, characterized by inhomogeneous rheological properties, embedded within a homogeneous elastic medium. The evolution in time of slip and stress over the crack plane is studied through a truncated expansion in Chebyshev polynomials, and convergence is proved to be fast in the simple examples considered. Finite-stress solutions are found which are compatible with constitutive relations of elasto-plastic materials and furthermore these allow us to simulate creep propagation and stress transfer between locked and unlocked fault segments. This model provides a simple interpretation of the shallow depth of the seismogenic layer observed in several areas of the world and lends itself to modelling creep processes during either post-seismic rebound or pre-seismic stress buildup. Stress transfer is accomplished mostly by the slow extension of the creeping section. During a seismic cycle it is envisaged that different regimes dominate over deep, intermediate and shallow sections of faults: (i) slow pre-seismic stress build-up accompanied by creep and stress migration toward intermediate depths; (ii) brittle fracture over shallow and intermediate sections of faults; (iii) post-seismic rebound over intermediate and deep sections of faults. The present crack model, while providing finite-stress solutions, allows a better understanding of how stress may accommodate at different depths over a fault plane during a seismic cycle.  相似文献   

11.
Summary. Current models of the structure of an active fault zone recognize two important subdivisions – an upper zone, extending to mid-crustal depths, in which processes associated with brittle fracture and friction dominate the fault behaviour, and a lower zone, extending into the mantle, within which stresses may be relieved by ductile flow. Anisotropy directly or indirectly induced by stress might occur throughout the fault zone, especially if caused by some form of stress-induced crack alignment. Dilatancy associated with high stresses is likely to be a very localized phenomenon in the vicinity of high strength regions (asperities), but alignments caused by subcritical crack growth at low stress and strain rate ( extensive-dilatancy anisotropy ) could give rise to anisotropy throughout the fault region.  相似文献   

12.
We develop a Galerkin finite element boundary integral equation method (GaBIEM) for spontaneous rupture propagation problems for a planar fault embedded in a homogeneous full 2-D space. A 2-D antiplane rupture propagation problem, with a slip-weakening friction law, is simulated by the GaBIEM. This method allows one to eliminate the strong singularities from the integral representation of the traction, and to separate explicitly the expression for the traction into an instantaneous component; static and time-dependent components with weakly (logarithmic) singular kernels; and a dynamic component and a quasi-static component, with continuous, bounded, kernels. Simulated results throw light into the performance of the GaBIEM and highlight differences with respect to that of the traditional, collocation, boundary integral equation method (BIEM). Both methods converge with a power law with respect to grid size, with different exponents. There is no restriction on the CFL stability number for the GaBIEM since an implicit, unconditionally stable method is used for the time integration. The error of the approximation increases with the time step, as expected, and it can remain below that of the BIEM.  相似文献   

13.
The spontaneous growth of a dynamic in-plane shear crack is simulated using a newly developed method of analysis in which no a priori constraint is required for the crack tip path, unlike in other classical studies. We formulate the problem in terms of boundary integral equations; the hypersingularities of the integration kernels are removed by taking the finite parts. Our analysis shows that dynamic crack growth is spontaneously arrested soon after the bending of the crack tips, even in a uniformly stressed medium with homogeneously distributed fracture strengths. This shows that the dynamics of crack growth has a significant effect on forming the non-planar crack shape, and consequently plays an essential role in the arrest of earthquake rupturing.  相似文献   

14.
Summary. Dynamical rupture process on the fault is investigated in a quasi-three-dimensional faulting model with non-uniform distributions of static frictions or the fracture strength under a finite shearing pre-stress. The displacement and stress time functions on the fault are obtained by solving numerically the equations of motion with a finite stress—fracture criterion, using the finite difference method.
If static frictions are homogeneous or weakly non-uniform, the rupture propagates nearly elliptically with a velocity close to that of P waves along the direction of pre-stress and with a nearly S wave velocity in the direction perpendicular to it. The rise time of the source function and the final displacements are larger around the centre of the fault. In the case when the static frictions are heavily non-uniform and depend on the location, the rupture propagation becomes quite irregular with appreciably decreased velocities, indicating remarkable stick-slip phenomena. In some cases, there remain unruptured regions where fault slip does not take place, and high stresses remain concentrated up to the final stage. These regions could be the source of aftershocks at a next stage.
The stick—slip faulting and irregular rupture propagation radiate high-frequency seismic waves, and the near-field spectral amplitudes tend to show an inversely linear frequency dependence over high frequencies for heavily non-uniform frictional faults.  相似文献   

15.
The Middle Durance fault system, southeastern France, is a slow active fault that produced moderate-size historical seismic events and shows evidence of at least one   M w ≳ 6.5  event in the last 29 000 yr. Based on dynamic rupture simulation, we propose earthquake scenarios that are constrained by knowledge of both the tectonic stress field and of the 3-D geometry of the Durance fault system. We simulate dynamic rupture interaction among several fault segmentations of different strikes, dips and rakes, using a 3-D boundary integral equation method. 50 combinations of reasonable stress field orientations, stress field amplitudes and hypocentre locations are tested. The probability of different rupture evolutions is then computed. Each segment ruptures mainly as a single event (44 per cent of the 50 simulations test in this paper). However, the probability that an event triggers simultaneously along three segments is high (26 per cent), leading to a potential rupture length of 45 km. Finally, 2 per cent of the simulations occur along four adjacent segments, producing the greatest total rupture length of 55 km. The simulation results show that the southernmost segment is most easily ruptured (40 per cent), because of its favourable orientation with respect to the tectonic stress and of its favourable location for interaction with the other segments. South-bound unilateral propagation is slightly preferable (41 per cent), compared to north-bound unilateral and bilateral propagation modes. Although, these rupture scenarios cannot be directly translated into probabilities of occurrence, they do provide a better insight as to which rupture scenarios are more likely, an important element to better estimate near-field strong ground motion and seismic hazard.  相似文献   

16.
We seek to understand how the stress interactions and the slip-weakening process combine within a non-coplanar, normal fault network to allow a slip instability to develop, and shape the final slip distribution on the system. In a first part, we perform a non-linear spectral analysis to investigate the conditions of stability and the process of slip initiation in an antiplane non-coplanar fault system subject to a slip-dependent friction law. That numerical model allows determining the zones that are able to slip within a fault network, as well as the location of the stress singularities. The resulting slip profiles on the faults show only a few different shapes, some of them with long, linear sections. This leads to formulate a general classification of slip profiles that can be used to infer the degree of fault interaction within any non-coplanar system. In a second part of work, we use our modelling to try reproducing the cumulative slip profiles measured on three real normal interacting faults forming a large-scale en echelon system. For that, we assume that cumulative slip profiles can be compared to the first static modal solution of our conceptual model. We succeed reproducing the profiles quite well using a variable weakening along the faults. Overall, the weakening rate decreases in the direction of propagation of the fault system. Yet, modelling the slip along the propagating, isolated termination segment of the system requires an unlikely distribution of weakening. This suggests that factors not considered in our analysis may contribute to slip profile shaping on isolated, propagating faults.  相似文献   

17.
Two distinct phases are commonly observed at the initial part of seismograms of large shallow earthquakes: low-frequency and low-amplitude waves following the onset of a P wave ( P 1) are interrupted by the arrival of the second impulsive phase P2 enriched with high-frequency components. This observation suggests that a large shallow earthquake involves two qualitatively different stages of rupture at its nucleation.
We propose a theoretical model that can naturally explain the above nucleation behaviour. The model is 2-D and the deformation is assumed to be anti-plane. A key clement in our model is the assumption of a zone in which numbers of pre-existing cracks are densely distributed; this cracked zone is a model for the fault zone. Dynamic crack growth nucleated in such a zone is intensely affected by the crack interactions, which exert two conflicting effects: one tends to accelerate the crack growth, and the other tends to decelerate it. The accelerating and decelerating effects are generally ascribable to coplanar and non-coplanar crack interactions, respectively. We rigorously treat the multiple interactions among the cracks, using the boundary integral equation method (BIEM), and assume the critical stress fracture criterion for the analysis of spontaneous crack propagation.
Our analysis shows that a dynamic rupture nucleated in the cracked zone begins to grow slowly due to the relative predominance of non-coplanar interactions. This process radiates the P1 phase. If the crack continues to grow, coalescence with adjacent coplanar cracks occurs after a short time. Then, coplanar interactions suddenly begin to prevail and crack growth is accelerated; the P2 phase is emitted in this process. It is interpreted that the two distinct phases appear in the process of the transition from non-coplanar to coplanar interaction predominance.  相似文献   

18.
Summary. We consider the problem of the unilateral extension of a two-dimensional anti-plane crack that initiates spontaneously at a point. The crack extends under the influence of cohesive resistance at the edge and dynamical friction along the crack walls. The stresses in the region beyond the edge of the crack are approximated so that they are exactly equal to the cohesive stresses near the edge of the crack, and are zero on the wavefront. An exact method of solving such problems is also given and can be used to determine the validity of the approximation. We find that the crack will not grow if the cohesion exceeds some critical value; this is consistent with an earlier result obtained by Knopoff, Mouton & Burridge for a similar one-dimensional model of crack propagation.  相似文献   

19.
Summary. The method proposed by Mendiguren to determine the source parameters from free oscillation data is applied to the 1970 July 31 deep Colombian earthquake. The results indicate a source propagating horizontally for about 150 km along the lithosphere and cutting across its width. The slab behaves as a guide for source propagation. The horizontal propagation velocity is determined as 3.8 km/s. The intensity of the source grew proportionately to the second power of the propagation distance. This rate of source intensity growth may be interpreted either by a fan-shaped fault model or by a cone-shaped volume source. The average slip and stress drop are estimated as 360 cm and 300 bar for the fault model. For the volume source model the transformational shear strain and stress are estimated as 11 × 10−5 and 160 bar. There is no evidence of a double couple radiation preceding the P origin time. It is shown that the isotropic and deviatoric components of the moment tensor cannot be uniquely resolved when only observations of a single mode are available. It is observed that, from a statistical basis, the available 0 Sn data for Colombian shock can be equally well explained by a pure deviatoric source model or by a source model including an isotropic component. Numerical experiments indicate that the inclusion of higher mode data does not change this situation. But, on the other hand, numerical experiments show that the available data and the scheme used for the inversion would not result in a solution including an artificial implosive component if the actual source were pure deviatoric. If the departure from a pure deviatoric source is produced by noise, it has to be non-random, as it could be produced by lateral heterogeneities not included in the inversion scheme.  相似文献   

20.
We study the effects of structural inhomogeneity on the quasi-static growth of strike-slip faults. A layered medium is considered, made up of an upper layer bounded by a free surface and welded to a lower half-space with different elastic property. Mode III crack is employed as a mathematical model of strike-slip fault, which is nucleated in the lower half-space and then propagates towards the interface. We adopt FEM-β, newly proposed analysis method for failure, to simulate the quasi-statistic crack growth governed by the stress distribution in layered media. Our results show that along planar traces across interfaces a compliant upper layer has significant effects on promoting/suppressing crack growth before/after its extension into the layer and vice versa for a rigid one. This proposes a possibility that surface breaks due to strike-slip faulting could be arrested by deposit layers at the topmost part of the Earth's crust.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号