Pore fluid viscosity effects on P‐ and S‐wave anisotropy in synthetic silica‐cemented sandstone with aligned fractures

Ultrasonic (500 kHz) P‐ and S‐wave velocity and attenuation anisotropy were measured in the laboratory on synthetic, octagonal‐shaped, silica‐cemented sandstone samples with aligned penny‐shaped voids as a function of pore fluid viscosity. One control (blank) sample was manufactured without fractures, another sample with a known fracture density (measured from X‐ray CT images). Velocity and attenuation were measured in four directions relative to the bedding fabric (introduced during packing of successive layers of sand grains during sample construction) and the coincident penny‐shaped voids (fractures). Both samples were measured when saturated with air, water (viscosity 1 cP) and glycerin (100 cP) to reveal poro‐visco‐elastic effects on velocity and attenuation, and their anisotropy. The blank sample was used to estimate the background anisotropy of the host rock in the fractured sample; the bedding fabric was found to show transverse isotropy with shear wave splitting (SWS) of 1.45 ± 1.18% (i.e. for S‐wave propagation along the bedding planes). In the fractured rock, maximum velocity and minimum attenuation of P‐waves was seen at 90° to the fracture normal. After correction for the background anisotropy, the fractured sample velocity anisotropy was expressed in terms of Thomsen's weak anisotropy parameters ε, γ & δ. A theory of frequency‐dependent seismic anisotropy in porous, fractured, media was able to predict the observed effect of viscosity and bulk modulus on ε and δ in water‐ and glycerin‐saturated samples, and the higher ε and δ values in air‐saturated samples. Theoretical predictions of fluid independent γ are also in agreement with the laboratory observations. We also observed the predicted polarisation cross‐over in shear‐wave splitting for wave propagation at 45° to the fracture normal as fluid viscosity and bulk modulus increases.     

Observations of fluid‐dependent shear‐wave splitting in synthetic porous rocks with aligned penny‐shaped fractures‡

P‐ and S‐wave velocity and attenuation coefficients (accurate to ±0.3% and ±0.2 dB/cm, respectively) were measured in synthetic porous rocks with aligned, penny‐shaped fractures using the laboratory ultrasonic pulse‐echo method. Shear‐wave splitting was observed by rotating the S‐wave transducer and noting the maximum and minimum velocities relative to the fracture direction. A block of synthetic porous rock of fracture density 0.0201 ± 0.0068 and fracture size 3.6 ± 0.38 mm (measured from image analysis of X‐ray CT scans) was sub‐sampled into three 20–30 mm long, 50 mm diameter core plugs oriented at 0°, 45° and 90° to the fracture normal (transversely isotropic symmetry axis). Full waveform data were collected over the frequency range 500–1000 kHz for both water and glycerin saturated cores to observe the effect of pore fluid viscosity at 1 cP and 100 cP, respectively. The shear‐wave splitting observed in the 90° core was 2.15 ± 0.02% for water saturated and 2.39 ± 0.02% for glycerin saturated, in agreement with the theory that suggests that the percentage splitting should be 100 times the fracture density and independent of the saturating fluid. In the 45° core, by contrast, splitting was 0.00 ± 0.02% for water saturation and ?0.77 ± 0.02% for glycerin saturation. This dependence on fracture orientation and pore fluid viscosity is consistent with the poro‐visco‐elastic theory for aligned, meso‐scale fractures in porous rocks. The results suggest the possible use of shear‐ or converted‐wave data to discriminate between fluids on the basis of viscosity variations.     

Fluid-dependent shear-wave splitting in a poroelastic medium with conjugate fracture sets

The dependence of shear‐wave splitting in fractured reservoirs on the properties of the filling fluid may provide a useful attribute for identifying reservoir fluids. If the direction of wave propagation is not perpendicular or parallel to the plane of fracturing, the wave polarized in the plane perpendicular to the fractures is a quasi‐shear mode, and therefore the shear‐wave splitting will be sensitive to the fluid bulk modulus. The magnitude of this sensitivity depends upon the extent to which fluid pressure can equilibrate between pores and fractures during the period of the deformation. In this paper, we use the anisotropic Gassmann equations and existing formulations for the excess compliance due to fracturing to estimate the splitting of vertically propagating shear waves as a function of the fluid modulus for a porous medium with a single set of dipping fractures and with two conjugate fracture sets, dipping with opposite dips to the vertical. This is achieved using two alternative approaches. In the first approach, it is assumed that the deformation taking place is quasi‐static: that is, the frequency of the elastic disturbance is low enough to allow enough time for fluid to flow between both the fractures and the pore space throughout the medium. In the second approach, we assume that the frequency is low enough to allow fluid flow between a fracture set and the surrounding pore space, but high enough so that there is not enough time during the period of the elastic disturbance for fluid flow between different fracture sets to occur. It is found that the second approach yields a much stronger dependence of shear‐wave splitting on the fluid modulus than the first approach. This is a consequence of the fact that at higher wave frequencies there is not enough time for fluid pressure to equilibrate and therefore the elastic properties of the fluid have a greater effect on the magnitude of the shear‐wave splitting.     

Seismic anisotropy in a hydrocarbon field estimated from microseismic data

The study of seismic anisotropy in exploration seismology is gaining interest as it provides valuable information about reservoir properties and stress directions. In this study we estimate anisotropy in a petroleum field in Oman using observations of shear‐wave splitting from microseismic data. The data set was recorded by arrays of borehole geophones deployed in five wells. We analyse nearly 3400 microearthquakes, yielding around 8500 shear‐wave splitting measurements. Stringent quality control reduces the number of reliable measurements to 325. Shear‐wave splitting modelling in a range of rock models is then used to guide the interpretation. The difference between the fast and slow shear‐wave velocities along the raypath in the field ranges between 0–10% and it is controlled both by lithology and proximity to the NE‐SW trending graben fault system that cuts the field formations. The anisotropy is interpreted in terms of aligned fractures or cracks superimposed on an intrinsic vertical transversely isotropic (VTI) rock fabric. The highest magnitudes of anisotropy are within the highly fractured uppermost unit of the Natih carbonate reservoir. Anisotropy decreases with depth, with the lowest magnitudes found in the deep part of the Natih carbonate formation. Moderate amounts of anisotropy are found in the shale cap rock. Anisotropy also varies laterally with the highest anisotropy occurring either side of the south‐eastern graben fault. The predominant fracture strikes, inferred from the fast shear‐wave polarizations, are consistent with the trends of the main faults (NE‐SW and NW‐SE). The majority of observations indicate subvertical fracture dip (>70° ). Cumulatively, these observations show how studies of shear‐wave splitting using microseismic data can be used to characterize fractures, important information for the exploitation of many reservoirs.     

A strategy for automated analysis of passive microseismic data to image seismic anisotropy and fracture characteristics

Monitoring of induced seismicity is gaining importance in a broad range of industrial operations from hydrocarbon reservoirs to mining to geothermal fields. Such passive seismic monitoring mainly aims at identifying fractures, which is of special interest for safety and productivity reasons. By analysing shear‐wave splitting it is possible to determine the anisotropy of the rock, which may be caused by sedimentary layering and/or aligned fractures, which in turn offers insight into the state of stress in the reservoir. We present a workflow strategy for automatic and effective processing of passive microseismic data sets, which are ever increasing in size. The automation provides an objective quality control of the shear‐wave splitting measurements and is based on characteristic differences between the two independent eigenvalue and cross‐correlation splitting techniques. These differences are summarized in a quality index for each measurement, allowing identification of an appropriate quality threshold. Measurements above this threshold are considered to be of good quality and are used in further interpretation. We suggest an automated inversion scheme using rock physics theory to test for best correlation of the data with various combinations of fracture density, its strike and the background anisotropy. This fully automatic workflow is then tested on a synthetic and a real microseismic data set.     

Fluid‐dependent shear‐wave splitting in fractured media

Azimuthal anisotropy in rocks can result from the presence of one or more sets of partially aligned fractures with orientations determined by the stress history of the rock. A shear wave propagating in an azimuthally anisotropic medium splits into two components with different polarizations if the source polarization is not aligned with the principal axes of the medium. For vertical propagation of shear waves in a horizontally layered medium containing vertical fractures, the shear‐wave splitting depends on the shear compliance of the fractures, but is independent of their normal compliance. If the fractures are not perfectly vertical, the shear‐wave splitting also depends on the normal compliance of the fractures. The normal compliance of a fluid‐filled fracture decreases with increasing fluid bulk modulus. For dipping fractures, this results in a decrease in shear‐wave splitting and an increase in shear‐wave velocity with increasing fluid bulk modulus. The sensitivity of the shear‐wave splitting to fluid bulk modulus depends on the interconnectivity of the fracture network, the permeability of the background medium and on whether the fracture is fully or partially saturated.     

Distribution of gas hydrate in fractured reservoirs: Insights from anisotropic seismic measurements

Seismic properties of hydrate-bearing reservoirs that are affected significantly by the hydrate distribution are key for quantitative assessment of the reservoir. The knowledge of hydrate distribution in fractured reservoirs remains poorly understood. To obtain such knowledge, we measured and analyzed five anisotropic velocities needed to fully characterize the seismic anisotropy in an artificial sandstone with aligned fractures during hydrate formation associated with varying distribution. We showed that while the formation of hydrate improved the velocities, the improvement was more significant for hydrate saturation above 10%. We also showed that the increasing trends varied among the anisotropic velocities when hydrate saturation was above 10%. Specifically, the compressional wave velocity travelling vertical to the bedding plane and the shear wave velocity with polarization perpendicular to the bedding plane increased more rapidly than the other compressional and shear wave velocities, respectively. Interpretation of the anisotropic seismic results suggested that the hydrate tends to bind to the grains in the fractures at low hydrate saturation, and becomes to bridge the fracture surfaces when the hydrate saturation exceeds 10%. The results have provided new insights into the hydrate distribution and its resulting anisotropic seismic properties in fractured reservoirs. This will pave the way for the successful assessment of hydrate in fractured reservoirs.     

Anomalous electrical resistivity anisotropy in clean reservoir sandstones

We report novel laboratory measurements of the full electrical resistivity tensor in reservoir analogue quartzose sandstones with clay contents less than 1.5%. We show that clean, homogeneous, visually uniform sandstone samples typically display between 15% and 25% resistivity anisotropy with minimum resistivity normal to the bedding plane. Thin‐section petrography, analysis of fabric anisotropy, and comparison to finite‐element simulations of grain pack compaction show that the observed anisotropy symmetries and magnitudes can be explained by syn‐depositional and post‐depositional compaction processes. Our findings suggest that: electrical resistivity anisotropy is likely to be present in most clastic rocks as a consequence of ballistic deposition and compaction; compaction may be deduced from measurements of electrical anisotropy; and the anisotropy observed at larger scales in well logging and controlled‐source electromagnetic data, with maximum resistivity normal to bedding, is most likely the result of meso‐scale (10^?1 m–10¹ m) periodic layering of electrically dissimilar lithologies.     

Finite-difference modelling of S-wave splitting in anisotropic media

We have implemented a 3D finite‐difference scheme to simulate wave propagation in arbitrary anisotropic media. The anisotropic media up to orthorhombic symmetry were modelled using a standard staggered grid scheme and beyond (monoclinic and triclinic) using a rotated staggered grid scheme. The rationale of not using rotated staggered grid for all types of anisotropic media is that the rotated staggered grid schemes are more expensive than standard staggered grid schemes. For a 1D azimuthally anistropic medium, we show a comparison between the seismic data generated by our finite‐difference code and by the reflectivity algorithm; they are in excellent agreement. We conducted a study on zero‐offset shear‐wave splitting using the finite‐difference modelling algorithm using the rotated staggered grid scheme. Our S‐wave splitting study is mainly focused on fractured media. On the scale of seismic wavelenghts, small aligned fractures behave as an equivalent anisotropic medium. We computed the equivalent elastic properties of the fractures and the background in which the fractures were embedded, using low‐frequency equivalent media theories. Wave propagation was simulated for both rotationally invariant and corrugated fractures embedded in an isotropic background for one, or more than one, set of fluid‐filled and dry fractures. S‐wave splitting was studied for dipping fractures, two vertical non‐orthogonal fractures and corrugated fractures. Our modelling results confirm that S‐wave splitting can reveal the fracture infill in the case of dipping fractures. S‐wave splitting has the potential to reveal the angle between the two vertical fractures. We also notice that in the case of vertical corrugated fractures, S‐wave splitting is sensitive to the fracture infill.     

Anisotropic azimuthal surface-wave inversion for a vertically fractured and transverse isotropy medium: Theory and examples from a controlled field experiment

Fractures in elastic media add compliance to a rock in the direction normal to the fracture strike. Therefore, elastic wave velocities in a fractured rock will vary as a function of the energy propagation direction relative to the orientation of the aligned fracture set. Anisotropic Thomson–Haskell matrix Rayleigh-wave equations for a vertically transverse isotropic media can be used to model surface-wave dispersion along the principal axes of a vertically fractured and transversely isotropic medium. Furthermore, a workflow combining first-break analysis and azimuthal anisotropic Rayleigh-wave inversion can be used to estimate P-wave and S-wave velocities, Thomsen's ε, and Thomsen's δ along the principal axes of the orthorhombic symmetry. In this work, linear slip theory is used to map our inversion results to the equivalent vertically fractured and transversely isotropic medium coefficients. We carried out this inversion on a synthetic example and a field example. The synthetic data example results show that joint estimation of S-wave velocities with Thomsen's parameters ε and δ along normal and parallel to the vertical fracture set is reliable and, when mapped to the corresponding vertically fractured and transversely isotropic medium, provides insight into the fracture compliances. When the inversion was carried out on the field data, results indicated that the fractured rock is more compliant in the azimuth normal to the visible fracture set orientation and that the in situ normal fracture compliance to tangential fracture compliance ratio is less than half, which implies some cementation may have occurred along the fractures. Such an observation has significant implications when modelling the transport properties of the rock and its strength. Both synthetic and field examples show the potential of azimuthal anisotropic Rayleigh-wave inversion as the method can be further expanded to a more general case where the vertical fracture set orientation is not known a priori.     

AVO investigations of shallow marine sediments

Amplitude‐variation‐with‐offset (AVO) analysis is based on the Zoeppritz equations, which enable the computation of reflection and transmission coefficients as a function of offset or angle of incidence. High‐frequency (up to 700 Hz) AVO studies, presented here, have been used to determine the physical properties of sediments in a shallow marine environment (20 m water depth). The properties that can be constrained are P‐ and S‐wave velocities, bulk density and acoustic attenuation. The use of higher frequencies requires special analysis including careful geometry and source and receiver directivity corrections. In the past, marine sediments have been modelled as elastic materials. However, viscoelastic models which include absorption are more realistic. At angles of incidence greater than 40°, AVO functions derived from viscoelastic models differ from those with purely elastic properties in the absence of a critical angle of incidence. The influence of S‐wave velocity on the reflection coefficient is small (especially for low S‐wave velocities encountered at the sea‐floor). Thus, it is difficult to extract the S‐wave parameter from AVO trends. On the other hand, P‐wave velocity and density show a considerably stronger effect. Attenuation (described by the quality factor Q) influences the reflection coefficient but could not be determined uniquely from the AVO functions. In order to measure the reflection coefficient in a seismogram, the amplitudes of the direct wave and the sea‐floor reflection in a common‐midpoint (CMP) gather are determined and corrected for spherical divergence as well as source and streamer directivity. At CMP locations showing the different AVO characteristics of a mud and a boulder clay, the sediment physical properties are determined by using a sequential‐quadratic‐programming (SQP) inversion technique. The inverted sediment physical properties for the mud are: P‐wave velocity α=1450±25 m/s, S‐wave velocity β=90±35 m/s, density ρ=1220±45 kg/m³, quality factor for P‐wave Q_P=15±200, quality factor for S‐wave Q_S=10±30. The inverted sediment physical properties for the boulder clay are: α=1620±45 m/s,β=360±200 m/s,ρ=1380±85 kg/m³,Q_P=790±660,Q_S=25±10.     

Seismic dispersion and attenuation in layered porous rocks with fractures of varying orientations

Wave‐induced fluid flow plays an important role in affecting the seismic dispersion and attenuation of fractured porous rocks. While numerous theoretical models have been proposed for the seismic dispersion and attenuation in fractured porous rocks, most of them neglect the wave‐induced fluid flow resulting from the background anisotropy (e.g. the interlayer fluid flow between different layers) that can be normal in real reservoirs. Here, according to the theories of poroelasticity, we present an approach to study the frequency‐dependent seismic properties of more realistic and complicated rocks, i.e. horizontally and periodically layered porous rock with horizontal and randomly orienting fractures, respectively, distributed in one of the two periodical layers. The approach accounts for the dual effects of the wave‐induced fluid flow between the fractures and the background pores and between different layers (the interlayer fluid flow). Because C₃₃ (i.e., the modulus of the normally incident P‐wave) is directly related to the P‐wave velocity widely measured in the seismic exploration, and its comprehensive dispersion and attenuation are found to be most significant, we study mainly the effects of fracture properties and the stiffness contrast between the different layers on the seismic dispersion and attenuation of C₃₃. The results show that the increasing stiffness contrast enhances the interlayer fluid flow of the layered porous rocks with both horizontal and randomly orienting fractures and weakens the wave‐induced fluid flow between the fractures and the background pores, especially for the layered porous rock with horizontal fractures. The modelling results also demonstrate that for the considered rock construction, the increasing fracture density reduces the interlayer fluid flow while improves the dispersion and attenuation in the fracture‐relevant frequency band. Increasing fracture aspect ratio is found to reduce the dispersion and attenuation in the fracture‐relevant frequency band only, especially for the layered porous rock with horizontal fractures.     

Effect of fracture fill on frequency‐dependent anisotropy of fractured porous rocks

In fractured reservoirs, seismic wave velocity and amplitude depend on frequency and incidence angle. Frequency dependence is believed to be principally caused by the wave‐induced flow of pore fluid at the mesoscopic scale. In recent years, two particular phenomena, i.e., patchy saturation and flow between fractures and pores, have been identified as significant mechanisms of wave‐induced flow. However, these two phenomena are studied separately. Recently, a unified model has been proposed for a porous rock with a set of aligned fractures, with pores and fractures filled with two different fluids. Existing models treat waves propagating perpendicular to the fractures. In this paper, we extend the model to all propagation angles by assuming that the flow direction is perpendicular to the layering plane and is independent of the loading direction. We first consider the limiting cases through poroelastic Backus averaging, and then we obtain the five complex and frequency‐dependent stiffness values of the equivalent transversely isotropic medium as a function of the frequency. The numerical results show that, when the bulk modulus of the fracture‐filling fluid is relatively large, the dispersion and attenuation of P‐waves are mainly caused by fractures, and the values decrease as angles increase, almost vanishing when the incidence angle is 90° (propagation parallel to the fracture plane). While the bulk modulus of fluid in fractures is much smaller than that of matrix pores, the attenuation due to the “partial saturation” mechanism makes the fluid flow from pores into fractures, which is almost independent of the incidence angle.     

Non-hyperbolic moveout inversion of wide-azimuth P-wave data for orthorhombic media

The azimuthally varying non‐hyperbolic moveout of P‐waves in orthorhombic media can provide valuable information for characterization of fractured reservoirs and seismic processing. Here, we present a technique to invert long‐spread, wide‐azimuth P‐wave data for the orientation of the vertical symmetry planes and five key moveout parameters: the symmetry‐plane NMO velocities, V⁽¹⁾_nmo and V⁽²⁾_nmo , and the anellipticity parameters, η⁽¹⁾, η⁽²⁾ and η⁽³⁾ . The inversion algorithm is based on a coherence operator that computes the semblance for the full range of offsets and azimuths using a generalized version of the Alkhalifah–Tsvankin non‐hyperbolic moveout equation. The moveout equation provides a close approximation to the reflection traveltimes in layered anisotropic media with a uniform orientation of the vertical symmetry planes. Numerical tests on noise‐contaminated data for a single orthorhombic layer show that the best‐constrained parameters are the azimuth ? of one of the symmetry planes and the velocities V⁽¹⁾_nmo and V⁽²⁾_nmo , while the resolution in η⁽¹⁾ and η⁽²⁾ is somewhat compromised by the trade‐off between the quadratic and quartic moveout terms. The largest uncertainty is observed in the parameter η⁽³⁾ , which influences only long‐spread moveout in off‐symmetry directions. For stratified orthorhombic models with depth‐dependent symmetry‐plane azimuths, the moveout equation has to be modified by allowing the orientation of the effective NMO ellipse to differ from the principal azimuthal direction of the effective quartic moveout term. The algorithm was successfully tested on wide‐azimuth P‐wave reflections recorded at the Weyburn Field in Canada. Taking azimuthal anisotropy into account increased the semblance values for most long‐offset reflection events in the overburden, which indicates that fracturing is not limited to the reservoir level. The inverted symmetry‐plane directions are close to the azimuths of the off‐trend fracture sets determined from borehole data and shear‐wave splitting analysis. The effective moveout parameters estimated by our algorithm provide input for P‐wave time imaging and geometrical‐spreading correction in layered orthorhombic media.     

Error in shear-wave polarization and time splitting

Shear‐wave polarization and time delay are attributes commonly used for fracture detection and characterization. In time‐lapse analysis these parameters can be used as indicators of changes in the fracture orientation and density. Indeed, changes in fracture characteristics provide key information for increased reservoir characterization and exploitation. However, relative to the data uncertainty, is the comparison of these parameters over time statistically meaningful? We present the uncertainty in shear‐wave polarization and time delay as a function of acquisition uncertainties, such as receiver and source misorientation, miscoupling and band‐limited random noise. This study is applied to a time‐lapse borehole seismic survey, recorded in Vacuum Field, New Mexico. From the estimated uncertainties for each survey, the uncertainty in the difference between the two surveys is 31° for the shear‐wave polarization angle and 4 ms for the shear‐wave time delay. Any changes in these parameters greater than these error estimates can be interpreted with confidence. This analysis can be applied to any time‐lapse measurement to provide an interval of confidence in the interpretation of shear‐wave polarization angles and time splitting.     

Effect of sub‐core scale heterogeneities on acoustic and electrical properties of a reservoir rock: a CO2 flooding experiment of brine saturated sandstone in a computed tomography scanner

The effect of sub‐core scale heterogeneity on fluid distribution pattern, and the electrical and acoustic properties of a typical reservoir rock was studied by performing drainage and imbibition flooding tests with CO₂ and brine in a laboratory. Moderately layered Rothbach sandstone was used as a test specimen. Two core samples were drilled; one perpendicular and the other parallel to the layering to allow injection of fluids along and normal to the bedding plane. During the test 3D images of fluid distribution and saturation levels were mapped by an industrial X‐ray CT‐scanner together with simultaneous measurement of electrical resistivity, ultrasonic velocities as well as amplitudes. The results showed how the layering and the flooding direction influenced the fluid distribution pattern and the saturation level of the fluids. For a given fluid saturation level, the measured changes in the acoustic and electrical parameters were affected by both the fluid distribution pattern and the layering orientation relative to the measurement direction. The P‐wave amplitude and the electrical resistivity were more sensitive to small changes in the fluid distribution patterns than the P‐wave velocity. The change in amplitude was the most affected by the orientation of the layering and the resulting fluid distribution patterns. In some instances the change due to the fluid distribution pattern was higher than the variation caused by the change in CO₂ saturation. As a result the Gassmann relation based on ‘uniform' or ‘patchy' saturation pattern was not suitable to predict the P‐wave velocity variation. Overall, the results demonstrate the importance of core‐imaging to improve our understanding of fluid distribution patterns and the associated effects on measured rock‐physics properties.     

Characterization of dipping fractures in a transverselyisotropic background

Although it is believed that natural fracture sets predominantly have near‐vertical orientation, oblique stresses and some other mechanisms may tilt fractures away from the vertical. Here, we examine an effective medium produced by a single system of obliquely dipping rotationally invariant fractures embedded in a transversely isotropic with a vertical symmetry axis (VTI) background rock. This model is monoclinic with a vertical symmetry plane that coincides with the dip plane of the fractures. Multicomponent seismic data acquired over such a medium possess several distinct features that make it possible to estimate the fracture orientation. For example, the vertically propagating fast shear wave (and the fast converted PS‐wave) is typically polarized in the direction of the fracture strike. The normal‐moveout (NMO) ellipses of horizontal reflection events are co‐orientated with the dip and strike directions of the fractures, which provides an independent estimate of the fracture azimuth. However, the polarization vector of the slow shear wave at vertical incidence does not lie in the horizontal plane – an unusual phenomenon that can be used to evaluate fracture dip. Also, for oblique fractures the shear‐wave splitting coefficient at vertical incidence becomes dependent on fracture infill (saturation). A complete medium‐characterization procedure includes estimating the fracture compliances and orientation (dip and azimuth), as well as the Thomsen parameters of the VTI background. We demonstrate that both the fracture and background parameters can be obtained from multicomponent wide‐azimuth data using the vertical velocities and NMO ellipses of PP‐waves and two split SS‐waves (or the traveltimes of PS‐waves) reflected from horizontal interfaces. Numerical tests corroborate the accuracy and stability of the inversion algorithm based on the exact expressions for the vertical and NMO velocities.     

Shear Wave Splitting Analysis to Estimate Fracture Orientation and Frequency Dependent Anisotropy

Shear wave splitting is a well-known method for indication of orientation, radius, and length of fractures in subsurface layers. In this paper, a three component near offset VSP data acquired from a fractured sandstone reservoir in southern part of Iran was used to analyse shear wave splitting and frequency-dependent anisotropy assessment. Polarization angle obtained by performing rotation on radial and transverse components of VSP data was used to determine the direction of polarization of fast shear wave which corresponds to direction of fractures. It was shown that correct implementation of shear wave splitting analysis can be used for determination of fracture direction. During frequencydependent anisotropy analysis, it was found that the time delays in shearwaves decrease as the frequency increases. It was clearly demonstrated throughout this study that anisotropy may have an inverse relationship with frequency. The analysis presented in this paper complements the studied conducted by other researchers in this field of research.     

Crack-induced anisotropy in the crust from shear wave splitting observed in Tangshan region,North China

Using the cross correlation function analysis method, this paper discusses shear wave splitting and crack-induced anisotropy in the crust beneath Tangshan, North China, by the digital data from Tangshan strong ground motion temporary arrays. Sixteen of twenty-one stations in the arrays recorded earthquake events available for studying from 1982 to 1984. Having calculated 131 available records, we get slower shear wave time delayτ and faster shear wave polarization azimuthPaz in Tangshan region, and the cracks densityε is got further from them. The analysis shows that the stress field is very complicated in Tangshan region and has strongly regional feature. Because of the complicated distribution of faults, different shear wave splitting characteristics are shown in 16 stations, scatteredτ and differentPaz. And they also were observed that theτ andPaz values were diverse within the time scale of hours in more than one station. In Tangshan region the average results ofτ,Paz andε are 0.0071 s · km^?1, northwest-west near to east-west and 0.022 respectively. Meantime, the standard deviations were calculated in this paper.