Numerical simulation of the propagation of P waves in fractured media

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 ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ 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 ) ⁻¹ 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.     

Observation and modelling of fault-zone fracture seismic anisotropy - I. P, SV and SH travel times

Summary. Three-component VSP borehole seismograms taken in the vicinity of an active normal fault in California show strong systematic shear-wave splitting that increases with proximity to the fault. Using Červený's method of characteristics for ray tracing in anisotropic heterogeneous media and Hudson's formulation of elastic constants for media-bearing aligned fractures, we have fitted a suite of P, SV and SH hanging-wall and foot-wall travel times with a simple model of aligned fractures flanking the fault zone. The dominant fracture set is best modelled as parallel to the fault plane and increasing in density with approach to the fault. The increase in fracture density is non-uniform (power law or Gaussian) with respect to distance to the fault. Although the hanging-wall and the foot-wall rock are petrologically the same unit, the fracture halo is more intense and extensive in the hanging wall than in the foot wall. Upon approach to the fault plane, the fracture density or fracture-density gradient becomes too great for the seismic response to be computed by Hudson–Červený procedures (the maximum fracture density that can be modelled is about 0.08). Within this 25 m fracture domain it appears more useful to model the fault and near field fractures as a low-velocity waveguide. We observe production of trapped waves within the confines of the intense fracture interval.     

Seismic waves in rocks with fluids and fractures

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, CO₂, 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.     

Inferring the relative measure of principal stress components

The phase velocity and the attenuation coefficient of compressional seismic waves, propagating in poroelastic, fluid-saturated, laminated sediments, are computed analytically from first principles. The wavefield is found to be strongly affected by the medium heterogeneity. Impedance fluctuations lead to poroelastic scattering; variations of the layer compressibilities cause inter-layer flow (a 1-D macroscopic local flow). These effects result in significant attenuation and dispersion of the seismic wavefield, even in the surface seismic frequency range, 10–100 Hz. The various attenuation mechanisms are found to be approximately additive, dominated by inter-layer flow at very low frequencies. Elastic scattering is important over a broad frequency range from seismic to sonic frequencies. Biot's global flow (the relative displacement of solid frame and fluid) contributes mainly in the range of ultrasonic frequencies. From the seismic frequency range up to ultrasonic frequencies, attenuation due to heterogeneity is strongly enhanced compared to homogeneous Biot models. Simple analytical expressions for the P -wave phase velocity and attenuation coefficient are presented as functions of frequency and of statistical medium parameters (correlation lengths, variances). These results automatically include different asymptotic approximations, such as poroelastic Backus averaging in the quasi-static and the no-flow limits, geometrical optics, and intermediate frequency ranges.     

Effective anisotropic elastic constants for wave propagation through cracked solids

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.     

Seismic wave properties in time-dependent porosity homogeneous media

It is quantified the properties of seismic waves in fully saturated homogeneous porous media within the framework of Sahay's modified and reformulated poroelastic theory. The computational results comprise amplitude attenuation, velocity dispersion and seismic waveforms. They show that the behaviour of all four waves modelled as a function of offset, frequency, porosity, fluid viscosity and source bandwidth depicts realistic dissipation within the sonic–ultrasonic band. Therefore, it appears that there is no need to include material heterogeneity to model attenuation. By inference it is concluded that the fluid viscosity effects may be enhanced by dynamic porosity.     

Attenuation behaviour of tuffaceous sandstone and granite during microfracturing

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.     

Scattering of waves in elastic media containing passive and active cracks

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.     

Dynamical rupture process on a three-dimensional fault with non-uniform frictions and near-field seismic waves

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.     

High-frequency radiation from crack (stress drop) models of earthquake faulting

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.     

What can be learned from rotational motions excited by earthquakes?

One answer to the question posed in the title is that we will have more accurate data for arrival times of SH waves, because the rotational component around the vertical axis is sensitive to SH waves although not to P-SV waves. Importantly, there is another answer related to seismic sources, which will be discussed in this paper.
Generally, not only dislocations commonly used in earthquake models but also other kind of defects could contribute to producing seismic waves. In particular, rotational strains at earthquake sources directly generate rotational components in seismic waves. Employing the geometrical theory of defects, we obtain a general expression for the rotational motion of seismic waves as a function of the parameters of source defects.
Using this expression, together with one for translational motion, we can estimate the rotational strain tensor and the spatial variation of slip velocity in the source area of earthquakes. These quantities will be large at the edges of a fault plane due to spatially rapid changes of slip on the fault and/or a formation of tensile fractures.     

Modelling the compliance of crustal rock—II. Response to temporal changes before earthquakes

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.     

Reflected-diffracted waves in fracture zone models

Summary. The paper presents the results of modelling of diffracted and reflected-diffracted waves in fracture zones. The Berryhill method was used and the calculations were made for a profile perpendicular to the diffracting edge. Several homogeneous models of the Earth's crust, characterized by different values of crustal thickness, velocity and horizontal distance between shot point and diffracting edge were considered. A dependence of the relative amplitude of diffracted waves on the location of the diffracting edge is given. The pattern of the seismic wavefield depends upon the dimensions of the fracture zone. Amplitude curves of reflected-diffracted waves are presented for a series of models of fracture zones. The possibility of applying the amplitudes of reflected-diffracted wave trains to the interpretation of the structure of fracture zones in the Earth's crust is andysed for different types of fracture zones.     

Wave propagation in fluid-saturated media: waveform and spectral analysis

A new numerical approach to the solution of waves propagating in a fluid-saturated medium, using Biot's theory as a foundation, has important implications for oil reservoir management and earthquake prediction. A numerical scheme is developed using an exponential transformation that explicitly treats the petrophysical and fluid properties of the medium within the framework of a generalized model. The scheme accounts for wave dissipation and velocity modifications. The numerical solution is used to perform numerical experiments to study the dynamic behaviour of waves in a fluid-saturated medium at well-logging frequencies (15 kHz). The results from the numerical experiments indicate that the degree of saturation by a high-viscosity fluid (HVF) such as oil, the temperature and the porosity of a medium strongly influence the spectral power distribution, frequency content and the velocity of waves propagating through the medium. An increase in HVF saturation causes enhanced attenuation of the low-frequency components, and increases the seismic velocity. An increase in porosity, however, enriches the low-frequency components and decreases the seismic velocity. A spectral quantification procedure is suggested and used to obtain information about the petrophysical and fluid properties of the medium from the spectral characteristics of the transmitted waveform. The procedure involves segmentation of the energy or power distribution of the transmitted waveforms into specified energy bands. The energy or power in these bands is then estimated. The extracted quantification variables are found to have strong correlations with the degree of HVF saturation, and the temperature and the porosity of the medium.     

Seismic radiation from dynamic coalescence, and the reconstruction of dynamic source parameters on a planar fault

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  ω⁻²  -type displacement spectra.     

Acoustoelasticity in cracked solids

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.     

Evaluation of Engineered Geothermal Systems as a Heat Source for Oil Sands Production in Northern Alberta

This paper evaluates the application of geothermal energy by numerically modeling the heat extraction that would result from the injection of cold water into an artificially fractured hot dry rock (HDR). The HDR that would be utilized in Alberta is expected to be granite with a network of pre-existing natural fractures. However, to ensure a continued flow of injected water from the reservoir to the production wells, creation of additional fractures is required. Thus, the properties of these fractures are of prime importance to the efficiency of geothermal energy production. The fracture networks for the simulations were created using a numerical code and were converted into a grid format to be used in a commercial thermal simulator. A new approach to embed a complex fracture system into the numerical model was applied. Various properties of the fractures such as aperture, length, and spacing were changed and their absolute and relative effects on energy production were quantified and the results are presented in this paper. This modeling technique was also verified by comparison with the conventional dual porosity model and by performing a history match with real field data obtained from literature. The applicability of this approach to provide heat for oil sands extraction was investigated using the volumes of water currently needed in northern Alberta. Based on these constraints, numerical simulations were run to evaluate the optimum well spacing that would be required using a three-well configuration. In this simulation, the fracture parameters (density and aperture) were kept fixed assuming that they are not affected by cold water injection. The results of this study suggest that geothermal energy has a potential to be a sustainable form of thermal energy for oil sands extraction in northern Alberta.     

Compressional and shear velocities of dry and saturated jointed rock: a laboratory study

Summary. An experimental and theoretical study was made to investigate the effects of large, through-going fractures on the seismic velocities of dry and saturated rocks. Fractures were simulated using cut and ground surfaces and the velocities were measured in a variety of igneous and metamorphic rocks at ultrasonic frequencies under confining pressures to 200 MPa. The 'fractures' caused a decrease in compressional and shear velocities. The effect was greater at lower pressures, with higher average fracture frequency, for rougher surfaces, in dry rather than saturated rocks, and in rocks containing a lower microcrack porosity. By modelling the surfaces as two rough surfaces in contact, deforming elastically under stress, equations are developed which appear to account adequately for the measurements, at least where the assumption of elasticity is valid. Where non-elastic deformation occurs, the theory no longer applies. Application of the theory to field data from the literature indicates that the model has some validity and that extraction of fracture parameters from seismic measurements is complicated by the large number of variables involved, including fracture frequency or spacing, stress state, fracture roughness, and degree of saturation.     

Inversion of seismic attributes for velocity and attenuation structure

We have developed an inversion formuialion for velocity and attenuation structure using seismic attributes, including envelope amplitude, instantaneous frequency and arrival times of selected seismic phases. We refer to this approach as AFT inversion for amplitude, (instantaneous) frequency and time. Complex trace analysis is used to extract the different seismic attributes. The instantaneous frequency data are converted to t * using a matching procedure that approximately removes the effects of the source spectra. To invert for structure, ray-perturbation methods are used to compute the sensitivity of the seismic attributes to variations in the model. An iterative inversion procedure is then performed from smooth to less smooth models that progressively incorporates the shorter-wavelength components of the model. To illustrate the method, seismic attributes are extracted from seismic-refraction data of the Ouachita PASSCAL experiment and used to invert for shallow crustal velocity and attenuation structure. Although amplitude data are sensitive to model roughness, the inverted velocity and attenuation models were required by the data to maintain a relatively smooth character. The amplitude and t * data were needed, along with the traveltimes, at each step of the inversion in order to fit all the seismic attributes at the final iteration.     

Teleseismic Sn: a guided wave in the mantle

The short-period seismic phase Sn has been interpreted by Stephens & Isacks as a lid wave' in which the seismic energy is constrained to the uppermost few tens of kilometres of the mantle. We have extended their normal-mode interpretation for structures both with and without low-velocity zones (LVZ) in the upper mantle. We have used spherical, anelastic models of the Earth. For a model of an oceanic mantle with a LVZ, we agree that Sn is a lid wave for sources above 200–250 km, if only the onset of Sn is considered. The later portions of the Sn wave train sample the structure as deeply as the 420-km discontinuity. For deeper foci, the pseudo-lid wave does not appear to be generated; even the onset of Sn samples the deeper mantle structure. For a model of a continental mantle without a LVZ, in general, sources at all depths above the 420-km discontinuity appear to generate teleseismic Sn which samples the entire mantle as deeply as the discontinuity and which travels with a velocity significantly greater than the lid velocity. Thus the velocity of Sn may be an important diagnostic to determine whether or not a LVZ exists in the upper mantle.