首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The effect of cracks on the elastic properties of an isotropic elastic solid is studied when the cracks are saturated with a soft fluid. A polynomial equation in effective Poisson's ratio is obtained, whose coefficients are functions of Poisson's ratio of the uncracked solid, crack density and saturating fluid parameter. Elastic and dynamical constants used in Blot's theory of wave propagation in poroelastic solids are modified for the introduction of cracks. The effects of cracks on the velocities of three types of waves are observed numerically. The frequency equation is derived for the propagation of Rayleigh-type surface waves in a saturated poroelastic half-space lying under a uniform layer of liquid. Dispersion curves for a particular model of oceanic crust containing cracks are plotted. The effects of variations in crack density and saturation on the phase and group velocity are also analysed.  相似文献   

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

3.
Summary. The Oblique Seismic Experiment (OSE) has been proposed to increase the usefulness of the IPOD crustal borehole as a means of investigating layer 2 of oceanic crust. Specific objectives are: to determine the lateral extent of the structure intersected by the borehole, to analyse the role of cracks in the velocity structure of layer 2, to look for anisotropy which may be caused by large cracks with a preferred orientation and, finally, to measure attenuation in oceanic crust.
The first successful Oblique Seismic Experiment in oceanic crust was carried out in 1977 March in a hole 400 miles north of Puerto Rico. An adequate study of lateral velocity variations was impossible because the hole was not deep enough, the hole was inadequately logged, and the small scale basement topography was not known. In general both P - and S -wave velocity profiles suggest that the crack density decreases with depth in layer 2. Velocities at the bottom of layer 2 are the same as matrix velocities for basalt, implying that crack density may be negligible at this depth. No convincing evidence for anisotropy in either layer 2 or 3 is found from travel time analysis. The hole was not deep enough to measure attenuation from normal incidence shots and amplitudes were not consistent enough to obtain a measure of attenuation from long range shots.  相似文献   

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

5.
Hydraulic crack propagation in a porous medium   总被引：2，自引：0，他引：2  
We develop a model for the propagation of a fluid-filled crack in a porous medium. The problem is motivated by the mechanism whereby drainage networks may form in partially molten rock below the Earth's lithosphere. Other applications include the propagation of hydraulic fractures in jointed rocks and in oil drilling operations, and the formation of dessication cracks in soils. Motivated by the lithospheric problem, we study a situation in which gravity acts in the direction of crack propagation. The model couples the elastoliydrodynamic equations of crack propagation with a pore pressure field in the porous rock, which drives the fluid flow which supplies the crack. The effect of the pore flow is to include a diffusional term in the evolution equation for the crack width, thus allowing a crack initiated at the base of the lithosphere to propagate down into the asthenosphere. Asymptotic and numerical solutions are presented for this crack evolution. However, the predicted drainage of melt into this crack is tiny compared with the upward percolative melt migration, and the predicted width of cracks (millimetres) is much too small to allow propagation of melt into the lithosphere without freezing. As a mechanism to explain magma fracturing in the lithosphere, the process described here therefore requires further refinement.  相似文献   

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.
Summary. Following the classic work of Eshelby, the slip and stress distributions due to an elliptical plane shear crack are evaluated. The relation between average (or maximum) slip on the crack and the (constant) static stress drop, for faults of different aspect ratios, is found. The slip vector is not parallel to the applied stress but makes a small angle to it, except when the stress is applied along the major or minor axis of the ellipse. The stress -distribution around the crack shows that in addition to the expected stress concentration along the crack edge, there are broad regions of stress increase off the crack plane for circular and elliptical cracks, similar to those known to exist for in-plane but not for antiplane shear cracks. Whether the stress- intensity factor at the end of one axis is greater or less than that at the end of the other axis ( ka ≶ kb ), depends on the condition: √ b/a ≶ (1 − v ) where a and b are the semi-axes of the ellipse, ka and kb are the stress-intensity factors at the end of the a- and b -axes and v is Poisson's ratio. The total stress-intensity factor varies smoothly along the edge of the ellipse from one axis to the other and it is found that this variation is rather small.  相似文献   

8.
Laboratory measurements of ultrasonic wave propagation in tuffaceous sandstone (Kimachi, Japan) and granite (Iidate, Japan) were performed during increasing fracturing of the samples. The fracturing was achieved by unconfined uniaxial compression up to and beyond the point of macrofracture of the specimen using a constant low strain rate. The observed variation of wave velocity (up to 40 per cent) due to the development of micro- and macrofractures in the rock is interpreted by rock models relating velocity changes to damage and crack density. The calculated density of the newly formed cracks reaches higher values for the sandstone than for the granite. Using the estimated crack densities, the attenuation behaviour is interpreted in terms of different attenuation mechanisms; that is, friction and scattering. Rayleigh scattering as described by the model of Hudson (1981 ) may explain the attenuation qualitatively if the largest plausible crack dimensions are assumed in modelling.  相似文献   

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

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

11.
Scattering of wavefields in a 3-D medium that includes passive and/or active structures, is numerically solved by using the boundary integral equation method (BIEM). The passive structures are velocity anomalies that generate scattered waves upon incidence, and the active structures contain endogenous fracture sources, which are dynamically triggered by the dynamic load due to the incident waves. Simple models are adopted to represent these structures: passive cracks act as scatterers and active cracks as fracture sources. We form cracks using circular boundaries, which consist of many boundary elements. Scattering of elastic waves by the boundaries of passive cracks is treated as an exterior problem in BIEM. In the case of active cracks, both the exterior and interior problems need to be solved, because elastic waves are generated by fracturing with stress drop, and the growing crack boundaries scatter the incident waves from the outside of the cracks. The passive cracks and/or active cracks are randomly distributed in an infinite homogeneous elastic medium. Calculations of the complete waveform considering a single scatter show that the active crack has weak influence on the attenuation of first arrivals but strong influence on the amplitudes of coda waves, as compared with those due to the passive crack. In the active structures, multiple scattering between cracks and the waves triggered by fracturing strongly affect the amplitudes of first arrivals and coda waves. Compared to the case of the passive structures, the attenuation of initial phase is weak and the coda amplitudes decrease slowly.  相似文献   

12.
Non-linear elastic behaviour of damaged rocks   总被引：1，自引：0，他引：1  
The pervasive damage of rocks by microcracks and voids strongly affects their macroscopic elastic properties. To evaluate the damage effects, we derive here the macroscopic stress-strain relations for a 3-D elastic solid with non-interacting cracks embedded inside a homogeneous matrix. The cracks considered are oriented either perpendicular to the maximum tension axis, or perpendicular to the maximum compression axis. In the first case they dilate during loading and in the second case they contract during loading. We derive a solution for the elastic energy of this rock following the self-consistent scheme of Budiansky & O'Connell (1976). The solution describes the stress-strain relations in terms of Λd and μd, which are the modified Lame constants, and an additional parameter Λ. The latter accounts for the non-linear behaviour of the solid and is related to crack density. The solution predicts a non-linear elastic rheology even for an infinitesimal strain of ɛ < 0.001, and abrupt change in the elastic moduli when the loading reverses from uniaxial compression to uniaxial tension.
We use the derived solution to analyse rock-mechanics experiments in which beams of Indiana limestone were deformed under four-point loading. This configuration provides simultaneously the apparent tensile and compressive moduli for small strains. The apparent moduli fit the effective elastic moduli calculated with the present damage model well, including the differences between tensile and compressive moduli.  相似文献   

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

14.
We derive a set of non-hypersingular boundary integral equations, both elastodynamic and elastostatic, for the analysis of arbitrarily shaped 2-D anti-plane and in-plane cracks located in an infinite homogeneous isotropic medium, rendered in a unified nomenclature for all cases. The hypersingularities that appear in the usual formulations for the dynamic cases, existent both at the source point and at the wavefront, are removed by way of a regularization technique based on integration by parts. The equations for the in-plane cases are presented in terms of a local Cartesian coordinate system, one of the axes of which is always held locally tangential to the crack trace. The expressions for the elastic field at any point on the model plane are also given.
Our formulations are shown to yield accurate numerical results, as long as appropriate stabilization measures are taken in the numerical scheme. The numerical applicability of our method to non-planar crack problems is illustrated by simulations of dynamic growth of a hackly crack which has small off-plane side-branches. The results imply that the branching of a crack brings about a significant decrease in the crack-tip stress concentration level and consequently may play an essential role in the arrest of earthquake rupturing.  相似文献   

15.
It is pointed out, in the context of the boundary integral equation method (BIEM), that, in the mechanics of 2-D curved in-plane shear cracks, a smooth curve, along which the crack orientation changes continuously, and an abrupt kink, across which it changes discontinuously, are not equivalent to each other. The discrepancy is illustrated by numerical results, and a set of conceptual models is used to demonstrate analytically how the equations that govern the crack mechanics have inherently distinct forms depending on whether the crack orientation changes continuously or abruptly across a bend, as long as one abides by the principles of linear elasticity theory. This has serious implications for the numerical treatment of a curved crack, which can be modelled as a chain of finite elements that are connected either smoothly or at abrupt kinks, the two methods producing different numerical outcomes. No similar paradox arises in the cases of anti-plane shear or open in-plane cracks.  相似文献   

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

17.
Numerical simulation of the propagation of P waves in fractured media   总被引：1，自引：0，他引：1  
We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q −1 of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q −1 of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) −1 in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.  相似文献   

18.
summary . A new technique is presented for modelling the elastic constants of cracked structures with application to systems with weak concentrations of parallel cracks, and of simple biplanar and triplanar cracks. The velocities and Vp/Vs ratios of these anisotropic structures are used to provide quantitative models for some earthquake precursors. These results indicate the great importance of crack geometry to the behaviour of precursors. The velocities of saturated cracks appear to favour the dilatancy-diffusion model of precursory phenomena. Synthetic seismograms are calculated for propagation through possible dilatancy zones. The seismograms show some characteristic features which may be useful for the investigation of earthquake dilatancy.  相似文献   

19.
A new model that accounts for the stress dependence of the phase velocity of elastodynamic waves propagating in a cracked solid under compression is presented. The phase velocities of longitudinal and shear waves are derived from the effective elastic properties of a cracked solid, which are evaluated within the framework of Kachanov's approach. Following Kachanov, the extra-compliance tensor of the cracked solid is related to the crack compliances, which display a marked non-linear behaviour when subjected to a compressive load. Such non-linear behaviour is shown to be derived from the elastic interaction between the contacting crack faces under compression. This work does not address the effect of mutual interaction among cracks and the generation of higher harmonics due to the medium non-linearity. Numerical examples are presented that illustrate the phase velocity changes occurring in a solid with a random distribution of parallel cracks as a function of an external compressive load. A distinctive feature of the acoustoelastic effect in solids with large parallel fractures and in solids with distributions of aligned microcracks is also illustrated.  相似文献   

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

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

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