首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到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.
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.  相似文献   

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

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

5.
Ultrasonic detection technology is of great significance in the detection and evaluation of physical and mechanical properties of frozen soil, but wave propagation characteristics in frozen soil are unclear. Based on the three-phase composition of frozen saturated soil and the mixture theory, considering Bishop's effective stress formula, the wave propagation equations are establish for frozen saturated soil. In wave propagation, an entropy inequality was introduced to describe the coupling of different phases. The analytic expressions of propagation velocity and attenuation law of waves in frozen soil are obtained, and wave propagation characteristics in frozen saturated soil are discussed. Results show that four types of waves(i.e., P1, P2, P3 and S) are found in frozen saturated soil and all four wave types are dissipative waves, in which the attenuation of P3 is the maximum. The velocity of four waves increases sharply at the excitation frequency range of 10~3–10~9 Hz,but the wave velocity at high-frequency and low-frequency is almost constant. When volume ice content increases, the wave propagation velocity of P1 and S decreases dramatically, and the velocity of P2 increases gradually, but P3 velocity increases first and then decreases to zero with increasing saturation. The attenuation coefficients of P1 and S waves begins to increase gradually when the volume ice content is about 0.4, P2 increases first and then decreases with an increase of volume ice content and P3 increases with the volume ice content and decreases rapidly from extreme to zero.  相似文献   

6.
A large data set of amplitude measurements of minor and major arc Rayleigh waves in the period range 73–171 s is collected. By comparing these amplitudes with the amplitudes of synthetic waveforms calculated by mode summation, maps of lateral variations in the apparent attenuation structure of the Earth are constructed. An existing formalism for predicting the effects of focusing is employed to calculate amplitude perturbations for the same data set. These perturbations are used to construct 'pseudo‐attenuation' maps and these results are compared with the apparent attenuation maps calculated from the data. It is shown that variations in Rayleigh wave amplitude perturbations in the Earth are dominated by attenuation at long wavelengths (below about degree 8) and by elastic structure at shorter wavelengths. It is also shown that the linear approximation for focusing is successful at predicting Rayleigh wave amplitudes using existing phase velocity maps. These results indicate that future attempts to model the velocity structure of the Earth would be assisted by incorporating amplitude data and by jointly inverting for Q structure.  相似文献   

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

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

9.
Summary. A fluid-saturated cubic packing of like elastic spheres is taken to be in equilibrium under the effect of gravity and the effects of a superimposed low-frequency elastic wave are considered. In the first place, expressions for the wave velocity, dispersion and attenuation are derived for the dry packing. This dynamic theory leads to the result that, for very low frequencies, the wave velocity is proportional to the third root of the depth and not the sixth root as is obtained by using the effective elastostatic modulus of the packing. For the fluid-saturated packing, two waves, termed respectively the 'solid wave' and the 'fluid wave', are found to propagate. The 'solid wave' has the characteristics of a wave propagating within a dry packing whose parameters differ in a specified way from those of the original packing, whereas the 'fluid wave' has those of a wave within a homogeneous fluid with similarly modified parameters.  相似文献   

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

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

12.
Summary. It is well-known that the chemical environment and the thermodynamical conditions play a fundamental role in the physics of the fracture properties of solids which in turn appear relevant in the understanding of the earthquake mechanism and related precursory phenomena. We designed and built an experimental apparatus capable of measuring the fracture velocities with control of the physical and chemical environment as well as of the applied stress. The apparatus consists of a chamber where both the total pressure and the partial pressure of gases can be controlled. The stress is applied in mode I (tension) configuration. The crack is detected optically on the surface of the sample with an ordinary microscope and the velocity of propagation is directly measured by the increase in length in the elapsed time. The very good stability of the system over very long periods allows us to measure crack velocities down to 10−12 ms−1. The propagation under such low-velocity regimes appears interesting from a geological standpoint since it implies extremely low stress values, which are generally assumed to have no effect on crustal rocks. The first results obtained will be discussed, with particular regard to the apparent dramatic influence of the partial pressure of water vapour in the crack propagation velocity.  相似文献   

13.
The blockage of the L g wave by crustal barriers such as continental margins and graben structures has long been recognized as providing a very useful tool for mapping large-scale lateral crustal variations along the propagation path. Numerical simulation of L g -wave propagation in complex anelastic media using the pseudospectral method provides insight into the nature of the propagation process using both snapshots of the wavefield and synthetic seismograms. A variety of 2-D structures have been investigated, including the influence of sediments, crustal thickness and attenuation.
Thick sedimentary basins covering a graben structure can have a major influence, since they remove L g energy by generating P conversion and scattering–the principal mechanisms for strong L g attenuation across a graben. The reduction of the L g energy is reinforced by anelastic attenuation in the sediments as well as the influence of the gradually thinning crustal waveguide associated with an elevated Moho.
The extinction of L g in a sequence of explosions fired across the central graben of the North Sea can be simulated by numerical calculations for the structure derived from refraction experiments.  相似文献   

14.
Offset-dependent characteristics of seismic scattering are useful for characterizing fractured reservoirs. We use two models that have different background medium properties and different azimuthal AVO responses to study elastic wave propagation and scattering in gas-saturated, heterogeneously fractured reservoirs. Heterogeneous fracture density distributions are built through stochastic modelling. Synthetic seismograms are generated by 3-D finite difference modelling, and waveforms along crack-normal and strike directions are considered in this paper. The multiple signal classification (MUSIC) frequency estimator is used in waveform estimation to provide frequency-domain attributes related to seismic wave scattering by fracture heterogeneity. Our results indicate that the strength of the scattering field is a function of the background medium. The strength also increases with increasing fracture scatterer density and with decreasing correlation length of spatial variations of fracture density. The scattering field is weak at the top of the fractured reservoir. The first-order results are dominated by velocity anisotropy of the mean fracture density field. However, the base of the fractured reservoir corresponds to a strong scattering field on which fracture heterogeneity has a larger effect and is characterized by the loss of coherence.  相似文献   

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

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

17.
Summary. Asymptotic ray theory is applied to surface waves in a medium where the lateral variations of structure are very smooth. Using ray-centred coordinates, parabolic equations are obtained for lateral variations while vertical structural variations at a given point are specified by eigenfunctions of normal mode theory as for the laterally homogeneous case. Final results on wavefields close to a ray can be expressed by formulations similar to those for elastic body waves in 2-D laterally heterogeneous media, except that the vertical dependence is described by eigenfunctions of 'local' Love or Rayleigh waves. The transport equation is written in terms of geometrical-ray spreading, group velocity and an energy integral. For the horizontal components there are both principal and additional components to describe the curvature of rays along the surface, as in the case of elastic body waves. The vertical component is decoupled from the horizontal components. With complex parameters the solutions for the dynamic ray tracing system correspond to Gaussian beams: the amplitude distribution is bell-shaped along the direction perpendicular to the ray and the solution is regular everywhere, even at caustics. Most of the characteristics of Gaussian beams for 2-D elastic body waves are also applicable to the surface wave case. At each frequency the solution may be regarded as a set of eigenfunctions propagating over a 2-D surface according to the phase velocity mapping.  相似文献   

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

19.
A quadrangle-grid velocity–stress finite difference method, based on a first-order hyperbolic system that is equivalent to Biot's equations, is developed for the simulation of wave propagation in 2-D heterogeneous porous media. In this method the velocity components of the solid material and of the pore fluid relative to that of the solid, and the stress components of three solid stresses and one fluid pressure are defined at different nodes for a staggered non-rectangular grid. The scheme uses non-orthogonal grids, allowing surface topography and curved interfaces to be easily modelled in the numerical simulation of seismic responses of poroelastic reservoirs. The free-surface conditions of complex geometry are achieved by using integral equilibrium equations on the surface, and the source implementations are simple. The algorithm is an extension of the quadrangle-grid finite difference method used for elastic wave equations.  相似文献   

20.
Effects of fractures on seismic-wave velocity and attenuation   总被引:1,自引:0,他引:1  
The effects of fractures on the seismic velocity and attenuation of a rock are investigated using theoretical results and experimental data. Fractures in a rock mass influence the traveltimes and amplitudes of seismic waves that have propagated through them. The displacement discontinuity model, recently employed in fracture investigations, is modified to describe the effect of fractures on seismic-wave velocity and attenuation. This new model, the modified displacement discontinuity model (MDD), is formulated in a way analogous to transmission-line analysis. The fractures are treated as transmission lines for the passage of seismic waves. The MDD takes into consideration realistic fracture parameters which include the fracture length, the fractional area of a fracture surface in contact, and the nature of the infilling material. A single fracture of varying geometric and material properties is shown to affect dramatically the transmission properties of a propagating waveform, and hence the seismic velocity and attenuation. These effects have been shown to result in a frequency-dependent velocity and attenuation. The sensitivity of the fracture parameters to seismic-wave velocity and attenuation was investigated and interesting results were obtained. Fracture parameters used in designing experimental models consisting of synthetically manufactured cracks were fed into the MDD and a well-known crack model, Hudson's model, for comparison. Velocities as a function of the incident-wave angle were obtained from both numerical models and were compared with the results from the experimental modelling. For P waves, the MDD model results show better agreement with those of the experimental model for all crack densities investigated than those from Hudson's model.  相似文献   

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

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