首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 328 毫秒
1.
Seismic anisotropy — the state of the art: II   总被引：1，自引：0，他引：1  
Summary. The theory, causes, observations, and possible applications of seismic anisotropy in the Earth have developed considerably since the previous state of the art paper was published in 1977. The behaviour of waves in layered anisotropic media is now much better understood and the evidence for seismic anisotropy indicates that anisotropy is likely to be present throughout much of the crust and upper mantle. The top few hundred kilometres of the mantle appears to be anisotropic with the orientations aligned by the present or palaeo stress-field. The upper part of the crust is frequently anisotropic, probably due to cracks differentially aligned by the non-lithostatic stresses. The possibility of being able to monitor crack geometry by seismic techniques opens a wide range of applications in currently important activities.  相似文献   

2.
The basis for earthquake prediction   总被引：3，自引：0，他引：3  
Summary. Recent advances in understanding the behaviour of shear waves propagating in the crust make the routine prediction of earthquakes seem practicable. Accumulating evidence suggests that most of the Earth's crust is pervaded by distributions of fluid-filled cracks and microcracks that are aligned by the contemporary stress-field so that the cracked rockmass is effectively anisotropic to seismic waves. This causes shear-waves to split, and shear-wave splitting is observed whenever shear-waves propagating along suitable raypaths in the crust are recorded by three-component instruments. These distributions of cracks are known as extensive-dilatancy anisotropy or EDA. Many characteristics of the crack- and stress-geometry can be monitored by analyzing shear-waves propagating through the cracked rockmass. Observations of temporal variations of the behaviour of shear-wave splitting in seismic gaps confirm these hypotheses, and suggest that stress changes before earthquakes may be monitored by analyzing shear-waves. In particular, monitoring earthquake preparation zones with three-component shear-wave vertical-seismic-profiles could lead to techniques for the routine prediction of earthquakes.  相似文献   

3.
Summary. The polarizations of shear waves recorded by networks of digital three-component seismometers immediately above small earthquakes near the North Anatolian Fault in Turkey display shear-wave splitting on almost all shear-wave seismograms recorded within the shear-wave window. This splitting is incompatible with source radiation-patterns propagating through simple isotropic structures but is compatible with effective anisotropy of the internal structure of the rock along the ray paths. This paper interprets the phenomena in terms of widespread crack-induced anisotropy. Distributions of stress-induced cracks model many features of the observations, and synthetic polarization diagrams calculated for propagation through simulated cracked rock are similar to the observed patterns. This evidence for widespread crack-induced anisotropy lends strong support to the hypothesis of extensive-dilatancy anisotropy (EDA) suggested by laboratory experiments in subcritical crack-growth. The crucial evidence confirming some form of EDA would be observations of temporal changes in shear-wave splitting as the stress field alters the crack density and crack geometry. There is some weak evidence for such temporal changes at one site, but further analysis of suitable digital three-component seismometer networks in seismic areas is required to confirm EDA.  相似文献   

4.
Summary. Polarization anomalies in seismic shear wavetrains, diagnostic of propagation through anisotropic media, have now been observed in dilatancy zones in seismic regions. Stress-induced dilatancy will open cracks with preferred orientations, which will be effectively anisotropic to short-period seismic waves. The polarization anomalies are due to the shear waves splitting, in propagation through anisotropic media, into components with different polarizations and different velocities. This writes characteristic signatures into the shear wavetrains. The paper examines ways in which the differential shear-wave anisotropy (the delay between the split shear-waves) varies with direction by plotting stereograms of the relative delays, and their polarizations, for possible dilatancy symmetry-systems. It seems likely, that if sufficient observations of these anomalies can be obtained at each stage of the dilatancy episode, it will be possible to estimate the symmetry directions of the dilatancy and the geometry of the stress-induced crack-system, as well as monitoring the progress of the dilatancy episode.  相似文献   

5.
Summary. Particle-motion plots of shear waves have been studied for the section FG of the FENNOLORA seismic experiment. Shear-wave splitting is observed on some records and the polarization of the first arriving shear waves show two peaks at about N35°W and N65°E. These results can be interpreted as being due to crack-induced anisotropy with the crack direction dominated by a (dominant) horizontal stress around N35°W. This is consistent with in situ stress measurements and focal mechanism studies in Scandinavia. the results show that seismic refraction experiments may be useful in providing evidence of crack-induced anisotropy in the stable continental crust.  相似文献   

6.
Summary. Almost all shear-waves from local earthquakes recorded on closely-spaced three-component seismometer networks deployed near the North Anatolian Fault, Turkey, in two experiments in 1979 and 1980, display shear-wave splitting. The observations are consistent with the presence of EDA (extensive-dilatancy anisotropy), distributions of fluid-filled cracks and microcracks aligned by the regional stress field. Temporal changes in the stress-field, which may occur before an earthquake, may modify the geometry and possibly the orientation of the EDA-microcracks, and lead to corresponding changes in the behaviour of the split shear-waves. A third experiment was undertaken in 1984 to investigate EDA further and to search for possible temporal variations of the polarization of the leading split shear-wave and the time delay between split shear-waves. Observations indicate that the polarization alignments, which are parallel to the strike of the parallel vertical EDA-cracks, are unaltered between 1979 and 1984, implying that the direction of the regional stress field has not changed significantly. Temporal changes in the stress field are more likely to cause changes in the crack density and/or aspect ratio, which would result in a corresponding change in time delay between the split shear-waves. We examine observations of time delay in relation to their propagation path with respect to the crack geometry since it is then possible to separate the effects of changes in crack density and changes in aspect ratio. With this procedure, a small temporal variation of time delays is found between 1979 and 1984, consistent with a decrease in crack density, and consequently a relaxation of stress, in this time period. No evidence was found for any observable variation of time delay over a six month observation span in 1984. We suggest that analysis of repeated shear-wave VSPs offers a technique for monitoring stress changes before earthquakes.  相似文献   

7.
There have been several claims that seismic shear waves respond to changes in stress before earthquakes. The companion paper develops a stress-sensitive model (APE) for the behaviour of low-porosity low-permeability crystalline rocks containing pervasive distributions of fluid-filled intergranular microcracks, and this paper uses APE to model the behaviour before earthquakes. Modelling with APE shows that the microgeometry and statistics of distributions of such fluid-filled microcracks respond almost immediately to changes in stress, and that the behaviour can be monitored by analysing seismic shear-wave splitting. The physical reasons for the coupling between shear-wave splitting and differential stress are discussed.
In this paper, we extend the model by using percolation theory to show that large crack densities are limited at the grain-scale level by the percolation threshold at which interacting crack clusters lead to pronounced increases in rock-matrix permeability. In the simplest formulation, the modelling is dimensionless and almost entirely constrained without free parameters. Nevertheless, APE modelling of the evolution of fluid-saturated rocks under stress reproduces the observed fracture criticality and the narrow range of shear-wave azimuthal anisotropy in crustal rocks. It also reproduces the behaviour of temporal variations in shear-wave splitting observed before and after the 1986, M = 6, North Palm Springs earthquake, Southern California, and several other smaller earthquakes.
The agreement of APE modelling with a wide range of observations confirms that fluid-saturated crystalline rocks are stress-sensitive and respond to changes in stress by critical fluid-rock interactions at the microscale level. This means that the effects of changes in stress and other parameters can be numerically modelled and monitored by appropriate observations of seismic shear waves.  相似文献   

8.
Summary. The section of the North Anatolian Fault lying near the city of Izmit, at the east of the Marmara Sea, has been identified as a seismic gap and the possible site of a future major earthquake. Previously published studies of records from an earthquake swarm within the gap (TDP1 and TDP2) provided the first evidence that shear-wave splitting occurs in earthquake source regions, a conclusion since verified by many studies at other locations. A third field study (TDP3) was mounted in the Izmit region during the summer of 1984. Observations were made over an eight-month period and included geomagnetic and geoelectric measurements in addition to a series of observations utilising dense arrays of three-component seismometers. Earthquake activity in the principal study area was monitored over a period of eight months. Records showed features similar to those observed in the earlier studies. In particular: (1) almost all shear waves emerging within the shear-wave window displayed shear-wave splitting; and (2) the polarizations of the first arriving (faster) split shear-waves showed sub-parallel alignments, characteristic of propagation through a distribution of parallel vertical cracks striking perpendicular to the minimum compressional stress.
These and other observations support the conclusion of earlier studies – that the upper crust is pervaded by distributions of micro-cracks aligned by stress, known as extensive-dilatancy anisotropy. A search for time dependence in shear-wave phenomena has revealed temporal variations in the delays between the split shear-waves throughout the course of the TDP3 study, but as yet this has not been correlated with specific earthquake activity.  相似文献   

9.
Seismic wave propagation through the earth is often strongly affected by the presence of fractures. When these fractures are filled with fluids (oil, gas, water, CO2, etc.), the type and state of the fluid (liquid or gas) can make a large difference in the response of the seismic waves. This paper summarizes recent work on methods of deconstructing the effects of fractures, and any fluids within these fractures, on seismic wave propagation as observed in reflection seismic data. One method explored here is Thomsen's weak anisotropy approximation for wave moveout (since fractures often induce elastic anisotropy due to non-uniform crack-orientation statistics). Another method makes use of some very convenient crack/fracture parameters introduced previously that permit a relatively simple deconstruction of the elastic and wave propagation behaviour in terms of a small number of crack-influence parameters (whenever this is appropriate, as is certainly the case for small crack densities). Then, the quantitative effects of fluids on these crack-influence parameters are shown to be directly related to Skempton's coefficient B of undrained poroelasticity (where B typically ranges from 0 to 1). In particular, the rigorous result obtained for the low crack density limit is that the crack-influence parameters are multiplied by a factor  (1 − B )  for undrained systems. It is also shown how fracture anisotropy affects Rayleigh wave speed, and how measured Rayleigh wave speeds can be used to infer shear wave speed of the fractured medium in some cases. Higher crack density results are also presented by incorporating recent simulation data on such cracked systems.  相似文献   

10.
Wave speeds and attenuation of elastic waves in material containing cracks   总被引：38，自引：0，他引：38  
Summary. Expressions now exist from which may be calculated the propagation constants of elastic waves travelling through material containing a distribution of cracks. The cracks are randomly distributed in position and may be randomly orientated. The wavelengths involved are assumed to be large compared with the size of the cracks and with their separation distances so that the formulae, based on the mean taken over a statistical ensemble, may reasonably be used to predict the properties of a single sample. The results are valid only for small concentrations of cracks.
Explicit expressions, correct to lowest order in the ratio of the crack size to a wavelength, are derived here for the overall elastic parameters and the overall wave speeds and attenuation of elastic waves in cracked materials where the mean crack is circular, and the cracks are either aligned or randomly orientated. The cracks may be empty or filled with solid or fluid material. These results are achieved on the basis of simply the static solution for an ellipsoidal inclusion under stress.
The extension to different distributions of orientation or to mixtures of different types of crack is quite straightforward.  相似文献   

11.
Summary. Theoretical developments of Hudson demonstrate how to calculate the variations of velocity and attenuation of seismic waves propagating through solids containing aligned cracks. The analysis can handle a wide variety of crack configurations and crack geometries. Hudson associates the velocity variations with effective elastic constants. In this paper we associate the variation of attenuation with the imaginary parts of complex effective elastic constants. These complex elastic constants permit the simulation of wave propagation through two-phase materials by the calculation of wave propagation through homogeneous anisotropic solids.  相似文献   

12.
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.  相似文献   

13.
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.  相似文献   

14.
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.  相似文献   

15.
P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.  相似文献   

16.
We analyse the relation between rock fabric, expressed by the preferred orientation of rock-forming minerals and microcracks, and elastic anisotropy of crystalline rock from the KTB pilot well. Detailed analyses of mineralogical composition, textures and microcrack fabrics were performed. In addition, ultrasonic velocity measurements of spherical samples in several directions were carried out at various confining pressures, and inverted in terms of the complete set of 21 elastic constants. By comparing the elastic tensors of the rocks at the final confining pressure (at which most of the microcracks are closed) with those at a lower pressure level, it is possible to separate the anisotropy induced by microcracks from that caused by mineral alignment. In contrast to previous work, no a priori knowledge of the type of anisotropy (triclinic, monoclinic, orthotropic etc.), or of the spatial orientation of the symmetry elements (planes, axes) of the cracked rock or of the intact rock is assumed. Furthermore, no restrictive assumptions on the orientation distribution function and the shape of the cracks are needed.
The results show that the elastic anisotropy characteristics, whether they are related to the microcracks or to the rock-forming minerals, are clearly correlated with the directly observed rock fabrics. We show that the symmetry directions of the mineral fabric and of microcrack fabric agree. A further result is that the microcrack-induced anisotropy dominates the other causes of anisotropy at confining pressures smaller than a few tens of megapascals, the situation being reversed at higher pressures. The laboratory data are quantitatively compared with sonic log data from the KTB well, showing the influence of pore fluids, effective pressure and crack density reduction on the anisotropy in situ .  相似文献   

17.
Summary . Expressions are available in the literature for the propagation constants of a 'mean' wave travelling in a material with a random distribution of cracks. These are approximations for small crack density and crack dimensions small compared with a wavelength. The formulae provided to second order in crack density by the method of smoothing are here extended to the case where the cracks consist of two or more sets, aligned in different directions.  相似文献   

18.
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.  相似文献   

19.
Summary. Data from Japanese local seismograph networks suggest that the stresses in double seismic zones are in-plate compression for the upper zone and in-plate tension for the lower zone; the stresses do not necessarily appear to be down-dip. It may therefore be possible to identify other double seismic zones on the basis of data which indicate that events with differing orientations of in-plate stresses occur in a given segment of slab.
A global survey of published focal mechanisms for intermediate depth earthquakes suggests that the stress in the slab is controlled, at least in part, by the age of the slab and the rate of convergence. Old and slow slabs are under in-plate tensile stresses and the amount of in-plate compression in the slab increases with increasing convergence rate or decreasing slab age. Young and fast slabs are an exception to this trend; all such slabs are down-dip tensile. Since these slabs all subduct under continents, they may be bent by continental loading. Double seismic zones are not a feature common to all subduction zones and are only observed in slabs which are not dominated by tensile or compressive stresses.
Unbending of the lithosphere and upper mantle phase changes are unlikely to be the causes of the major features of double zones, although they may contribute to producing some of their characteristics. Sagging or thermal effects, possibly aided by asthenospheric relative motion, may produce the local deviatoric stresses that cause double zones.  相似文献   

20.
In the southwest U.S., cracks in alluvial fan surface clasts have a preferred orientation independent of rock fabric and shape. In this paper, we show that differential insolation of incipient cracks of random orientations predicts a distribution of crack orientations consistent with field observations. In this model, crack growth by hydration and/or thermal weathering is primarily a function of local water content at the crack tip. Crack tips that experience minimal solar insolation maintain a greater average moisture and, hence, weather more rapidly than cracks that experience greater solar insolation. To show this, we used a numerical radiative transfer code to quantify the solar insolation of rectangular cracks at 35° N. latitude with a range of depths and orientations. The amount of solar energy reaching the bottom of each crack was calculated at 5-min intervals over the day for several days of the year to determine hourly, daily, seasonal, and annual energy deposition as a function of crack depth and orientation. By assuming that only crack orientations that effectively shield their interiors and minimize their water loss are able to grow, the pattern of cracks produced by the model is consistent with field observations. The annual average insolation, which controls water retention, is associated with the two primary modes of crack orientation. The effect of daily recharge by summer rains of the North American monsoon system is consistent with the observed deviations from these primary modes. Model results suggest that both the annual average insolation and the daily pattern of rainfall is recorded in the preferred crack orientations of surface clasts in the southwest U.S.  相似文献   

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

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