P-wave Anisotropy Research Articles

Theoretical modelling of seismic dispersion, attenuation and frequency-dependent anisotropy in a fluid-saturated porous rock with intersecting fractures

SUMMARYDetection of intersecting open fractures is an important task in many earth science domains. To quantify the seismic responses in the fluid-saturated porous rock with intersecting open fractures, we develop a theoretical model based on Biot's equations of dynamic poroelasticity. The seismic dispersion, attenuation and frequency-dependent anisotropy due to joint effects of fracture–background wave-induced fluid flow (FB-WIFF), and fracture–fracture wave-induced fluid flow (FF-WIFF), as well as elastic scattering are investigated. The numerical results on a fluid-saturated porous and fractured sandstone show that the characteristic frequency of FF-WIFF is controlled by fracture connectivity and fluid viscosity. Variations of fracture connectivity and fluid viscosity may result in the coupling of FF-WIFF with FB-WIFF or elastic scattering. When the fracture connectivity tends to zero, the FF-WIFF vanishes and FB-WIFF becomes most significant. Besides fracture connectivity and fluid viscosity, the fracture geometry and fluid bulk modulus also affect the magnitudes of these three mechanisms and their interplay. Due to effects of these three mechanisms, the P-wave anisotropy varies greatly with frequencies. Furthermore, the fracture intersection angle also influences the P-wave anisotropy significantly. Our model agrees well with previous models in the frequency limits and for the special case with parallel fractures. Since our model incorporates the effects of FF-WIFF, it has a great potential to be applied in the detection for effective fracture networks for fluid flow.

Seismic monitoring of the STIMTEC hydraulic stimulation experiment in anisotropic metamorphic gneiss

Abstract. In 2018 and 2019, we performed STIMulation tests with characterising periodic pumping tests and high-resolution seismic monitoring for improving prognosis models and real-time monitoring TEChnologies for the creation of hydraulic conduits in crystalline rocks (STIMTEC). The STIMTEC underground research laboratory is located at 130 m depth in the Reiche Zeche mine in Freiberg, Germany. The experiment was designed to investigate the rock damage resulting from hydraulic stimulation and to link seismic activity and enhancement of hydraulic properties in strongly foliated metamorphic gneiss. We present results from active and passive seismic monitoring prior to and during hydraulic stimulations. We characterise the structural anisotropy and heterogeneity of the reservoir rocks at the STIMTEC site and the induced high-frequency (>1 kHz) acoustic emission (AE) activity, associated with brittle deformation at the centimetre-to-decimetre scale. We derived the best velocity model per recording station from over 300 active ultrasonic transmission measurements for high-accuracy AE event location. The average P-wave anisotropy is 12 %, in agreement with values derived from laboratory tests on core material. We use a 16-station seismic monitoring network comprising AE sensors, accelerometers, one broadband sensor and one AE hydrophone. All instrumentation was removable, providing us with the flexibility to use existing boreholes for multiple purposes. This approach also allowed for optimising the (near)-real-time passive monitoring system during the experiment. To locate AE events, we tested the effect of different velocity models and inferred their location accuracy. Based on the known active ultrasonic transmission measurement points, we obtained an average relocation error of 0.26±0.06 m where the AE events occurred using a transverse isotropic velocity model per station. The uncertainty resulting from using a simplified velocity model increased to 0.5–2.6 m, depending on whether anisotropy was considered or not. Structural heterogeneity overprints anisotropy of the host rock and has a significant influence on velocity and attenuation, with up to 4 % and up to 50 % decrease on velocity and wave amplitude, respectively. Significant variations in seismic responses to stimulation were observed ranging from abundant AE events (several thousand per stimulated interval) to no activity with breakdown pressure values ranging between 6.4 and 15.6 MPa. Low-frequency seismic signals with varying amplitudes were observed for all stimulated intervals that are more correlated with the injection flow rate rather than the pressure curve. We discuss the observations from STIMTEC in context of similar experiments performed in underground research facilities to highlight the effect of small-scale rock, stress and structural heterogeneity and/or anisotropy observed at the decametre scale. The reservoir complexity at this scale supports our conclusion that field-scale experiments benefit from high-sensitivity, wide-bandwidth instrumentation and flexible monitoring approaches to adapt to unexpected challenges during all stages of the experiment.

Comparison of wave-propagation simulations in fractured domains using discrete fractures and equivalent media

SUMMARY Fractures largely control reservoir permeability and, therefore, it is of immediate importance to know the geometrical parameters of fracture sets and their effects on seismic data. To understand the interaction between the fractures and seismic wavefield, we performed numerical simulations of elastic wave propagation in fractured digital rocks (FDRs) using the discontinuous Galerkin method, the linear-slip model and sets of randomly distributed fractures. We compared the results with those obtained using Hudson’s equivalent media theory (HEM) and observed that, when the fracture density is 0.08 or less, there is no statistical difference between the FDR and HEM results; however, when the fracture density is higher than 0.08, the results of HEM diverge from those of FDR. Furthermore, HEM accuracy depends not only on the fracture density but also on the P- to S-wave velocity ratio. The P-wave anisotropy induced by the fractures is observed as a delay, which can be due to fracture density, length or a mixture of both. The Pwave is delayed in all directions, but mostly in the direction perpendicular to the fractures’ orientation.

Mantle structure and flow beneath the central-western US: Constraints from anisotropic tomography

To investigate lateral and depth variations of seismic anisotropy beneath the central-western United States, we determined a detailed 3-D model of P-wave anisotropic tomography by inverting a large number of arrival-time data of local and teleseismic events. Our results reveal significant azimuthal anisotropies in the crust and lithosphere, which are associated with ancient orogenic collisional and magmatic activities. As depth increases, the fast-velocity direction (FVD) pattern becomes gradually trended and small features fade away. There is a boundary in the FVD distribution, which separates the tectonically active region in the west from the stable cratonic region in the east. Frozen-in anisotropy with a NW-SE FVD is preserved in the thick Wyoming cratonic lithosphere that exhibits as a high-velocity (high-V) anomaly to a depth of ~250 km. In the asthenosphere beneath the western thin lithosphere, FVDs are generally parallel with the absolute motion direction of the North American plate due to shearing between the plate and the asthenosphere. In the deeper areas, the subducted and fragmented slab exhibiting as high-V anomalies leads to slab-related mantle flows. These results indicate that seismic anisotropies exist in both the lithosphere and asthenosphere with different geodynamic mechanisms and it is feasible to link the P-wave azimuthal anisotropy to lithospheric deformations, fossil anisotropy in the lithosphere, and flows in the asthenosphere.

Core-Derived Velocity Systematics, Mississippian Meramec Formation, Oklahoma

Summary Elastic wave velocities are commonly used to predict porosity, mineralogy, and lithology from formation properties. When only P-wave sonics are available in historical wells, systematics for predicting shear velocities are useful for developing elastic models. Although much research has been done on conventional reservoir velocity systematics, the equivalency for unconventional formations is still a work in progress. There has also been a limited number of research studies with laboratory measures published. Using laboratory pulse transmission ultrasonic data, we created a Vp-Vs systematic for the Meramec Formation in this study. The effects of porosity and mineralogy on velocities are explored, as well as a comparison of Meramec velocity systematics with well-established literature systematics. Vp and Vs measurements were taken on 385 dodecane-saturated core samples from seven Meramec wells (106 vertical and 279 horizontal plugs). S-wave and P-wave anisotropy in Meramec Formation samples used in this study are typically less than 10%. Each sample was also tested for porosity and mineralogy. We find that velocities are more sensitive to porosity than mineralogy by a factor of 10. Below are our equations for predicting Vp and Vs (in km/s), when only clay content and porosity are known. In these equations, φ is the volume fraction pores, and Clays is the weight fraction of clay. These equations are for those samples in which there is low P-wave and S-wave anisotropies:(1)Vp=6.4−1.2*Clays−15.4*φ(R2=0.5),(2)Vs=3.6−0.5*Clays−5.2*φ(R2=0.4). We suggest two methods for calculating Vs from Vp: Ignoring anisotropy, we combined both Vp and Vs measurements from all vertical plugs and low anisotropy horizontal plugs to create a single shear wave predictor; and considering anisotropy, Vp measurements from horizontal plugs were corrected using Thomsen’s compressional wave anisotropy parameter, after which a shear velocity predictor was generated. The shear wave predictors for dodecane-saturated measurements are as follows (all velocities are km/s):(3)Method 1: Vs= 0.90 + 0.42*Vp (R2=0.7),(4)Method 2: Vs= 0.80 + 0.45*Vp (R2=0.6). The residual and estimated error in Eq. 3 is slightly less than in Eq. 4. Even though there is a significant variance in measurement frequency, the Meramec velocity systematic shows good agreement with dipole wireline measurements using the first equation. The Meramec velocity systematics differ significantly from previously published systematics, such as the trend line by Greenberg and Castagna (1992) and the shale trend line by Vernik et al. (2018). Using the correlations by Greenberg and Castagna (1992) for limestone or dolomite, the shear velocities of the samples in this study cannot be predicted. These data have yielded shear wave systematics, which can be used in wireline and seismic investigations. The results suggest that the method of ignoring anisotropy yields a better Vs estimate than the one that takes anisotropy into account. Using well-established shear wave velocity systematics from the published literature can result in an estimated inaccuracy of greater than 16%. It is important to calibrate velocity systematics to the target formation.

Development of a New Empirical Relation to Assess P-wave Velocity Anisotropy of Rocks

Several physical parameters and anisotropy related to rock textural arrangements, schistosity and weakness planes such as cracks and joints affect the P-wave velocity (VP). First, VP anisotropy of quartz-mica schist as a common type of widespread metamorphic rock was compared with VP anisotropy of jointed homogeneous limestone specimens to clarify effect of these two different types of anisotropies. The results showed that the VP anisotropy of quartz-mica schist texture is stronger than the VP anisotropy of jointed limestone, because all body of quartz-mica schist specimens have VP anisotropy behavior. Many rocks are anisotropic and degree of anisotropy varies from one to another. Various investigations have been carried out on VP anisotropy but there is not a unique comprehensive relation to represent the influence of different degrees of anisotropy on the VP for different rocks. The relation between VP and angle (θ) between the axis of symmetry (perpendicular to weakness planes) with the wave propagation direction was analyzed for a wide range of anisotropy degrees using the results of nine different types of rocks including: Angouran quartz-mica schist, Golgohar mica schist, amphibole schist, mica-quart schist, Marcellus shale, Withby shale WUK47B, WUK70 and WUK2, and Veroia-Polymylos gneiss. A new simple empirical relation fitted to all groups of results was obtained to assess VP for different degrees of anisotropies with a good correlation of determination (R2 = 0.937), low RMSE (RMSE = 320 m/s) and low CV (CV = 7.0%). P wave velocity anisotropy can simply be predicted by the developed relation using only two parameters of VP0 and VP90 or ε that is the percentage change of VP0 with respect to VP90. A VP anisotropy classification diagram was also developed based on the different values of ε.

Feasibility study of fracture interpretation using multicomponent seismic data: SEAM Phase II Barrett model

Several P-wave azimuthal anisotropy studies have been conducted for the SEG Advanced Modeling Program Phase II Barrett model data. However, these analyses provide fracture property estimation that is inconsistent with the actual model properties. Therefore, we perform a feasibility study to understand the influence of the overburden and reservoir properties, and processing and inversion steps, which together determine the success of the fracture interpretation from seismic data. The 1D model properties (orthorhombic for overburden and reservoir) are first extracted from the actual Barrett model properties at two locations. Anisotropic prestack reflectivity modeling exposes the true orthorhombic response of the 1D medium in the form of common-offset and common-azimuth gathers. The true anisotropic response is obscured in the Barrett data (generated by finite-element modeling) due to the mild lateral velocity variations and orthorhombic anisotropy in the overburden. We then expose the reservoir anisotropic response using an isotropic overburden in the reflectivity modeling. This indicates that the P-wave velocity variation with azimuth (VVAZ) responses generated by the reservoir itself are weak, which leads to an unstable VVAZ inversion to estimate the interval normal moveout (NMO) velocity anisotropy. The reservoir thickness (125 m or 65 ms of two-way traveltime) or NMO velocity anisotropy (6%–7%) needs to be at least doubled to obtain a stable VVAZ inversion. Anisotropic geometric-spreading correction improves the amplitude variation with azimuth inversion results when reflectivity modeling models the orthorhombic overburden. The converted wave (C-wave) has a stronger VVAZ response compared to that of the P-wave. We suggest that the C-wave data could be useful to constrain fracture interpretation in the Barrett model. We conclude that the results of previous studies are due to the combination of the residual influence of overburden after processing and imaging and the weak anisotropy responses from the reservoir.

Influence of the microstructure on stress-dependent P-wave anisotropy in sandstone

SUMMARY To understand the factors that affect the stress-dependent P-wave velocity anisotropy, a method is proposed to simulate anisotropic microcracks and minerals based on the discrete element method (DEM). Laboratory triaxial tests and numerical simulations were performed on sandstone samples with bedding orientations parallel and perpendicular to the maximum principal stress. The ellipse fitting method was applied to analyse the variation in P-wave anisotropy. The micromechanism of stress-dependent P-wave anisotropy was revealed. The evolution of microcracks is the main reason for the change in P-wave anisotropy under compression. As the confining pressure increases, the magnitude of the P-wave anisotropy is reduced. The weakening of the P-wave anisotropy results from the decrease in the number of open microcracks. Under deviatoric stress loading, the P-wave anisotropy of the bedding-parallel sample in the axial direction is strengthened. Anisotropy reversal occurs in the bedding-normal sample. The microcrack behaviour depends on the direction of maximum principal stress. The variation in microcrack anisotropy induced by stress controls the evolution of P-wave velocity anisotropy. The stress at which anisotropic reversal occurs depends on the preferred orientation mineral. The DEM model offers the unique ability to directly examine the variation in microstructure anisotropy that causes the change in P-wave anisotropy.

Effects of melt-percolation, refertilization and deformation on upper mantle seismic anisotropy: constraints from peridotite xenoliths, Marie Byrd Land, West Antarctica

Abstract We report on the petrology, microstructure and seismic properties of 44 peridotite xenoliths extracted from the upper mantle beneath Marie Byrd Land (MBL), West Antarctica. The aim of this work is to understand how melt-rock reaction, refertilization, and deformation affected the seismic properties (velocities, anisotropy) of the West Antarctic upper mantle, in the context of MBL tectonic evolution and West Antarctic Rift System formation. Modal compositions, mineral major element compositions, microstructures, and crystallographic preferred orientations (CPOs) provide evidence for diachronous reactive melt percolation and refertilization. Olivine shows three main CPO patterns, the A-type, axial-[010], and axial-[100] texture types. Average seismic properties of the MBL mantle lithosphere are mainly controlled by the strength of olivine crystallographic texture. Reactive melt percolation and refertilization likely modified seismic velocities and anisotropy, as is suggested by a systematic decrease in maximum P-wave and S-wave anisotropies with increasing modal abundance of pyroxene. At larger spatial scales, the seismic properties of the MBL mantle xenoliths are dominated by the anisotropy resulting from the A-type olivine CPO. Variations between individual volcanic centres, however, attest to spatial variations in the mantle structure, potentially related to 3D deformation and the prolonged tectonic history of MBL.

Prediction of sweet spots in tight sandstone reservoirs based on anisotropic frequency-dependent AVO inversion

Abstract The frequency-dependent amplitude-versus-offset (FAVO) method has become a practical method for fluid detection in sand reservoirs. At present, most FAVO inversions are based on the assumption that reservoirs are isotropy, but the application effect is not satisfactory for fractured reservoirs. Hence, we analyse the frequency variation characteristics of anisotropy parameters in tight sandstone reservoirs based on a new petrophysical model, and propose a stepwise anisotropic FAVO inversion method to extract frequency-dependent attributes from prestack seismic field data. First, we combine the improved Brie's law with the fine-fracture model to analyse frequency-dependent characteristics of velocities and Thomsen anisotropy parameters at different gas saturations and fracture densities. Then, we derive an anisotropic FAVO inversion algorithm based on Rüger's approximation formula and propose a stepwise anisotropic FAVO inversion method to obtain the dispersions of anisotropy parameters. Finally, we propose a method that combines the inversion spectral decomposition with the stepwise anisotropy FAVO inversion and apply it to tight sand reservoirs in the Xinchang area. We use P-wave velocity dispersion and anisotropy parameter ε dispersion to optimise favourable areas. Numerical analysis results show that velocity dispersion of the P-wave is sensitive to fracture density, which can be used for fracture prediction in fractured reservoirs. In contrast, anisotropic parameter dispersion is sensitive to gas saturation and can be used for fluid detection. The seismic data inversion results show that velocity dispersion of the P-wave and anisotropic parameter dispersion are sensitive to fractured reservoirs in the second member of Xujiahe Group, which is consistent with logging interpretation results.

Seismic properties across an amphibolite- to greenschist-facies strain gradient (Neves area, eastern Alps, Italy): New considerations for shear zone imaging

Detecting the presence and geometry of crustal shear zones by geophysical methods relies on our understanding of the intrinsic parameters controlling the seismic properties of these deformed rocks, over the range of pressure-temperature conditions expected in the Earth's crust. To this end, we aimed to track changes in P-wave propagation velocity (VP), anisotropy (AVP) and S-wave splitting (AVS) across a natural shear zone using experimental and electron backscatter diffraction methods. Five samples were collected across a meter-scale shear zone developed under amphibolite- to greenschist-facies in the Neves metagranodiorite (Tauern Window, Eastern Alps). With increasing strain, seismic properties show a non-linear evolution resulting from interactions among structural, textural and mineralogical changes. Geometry and magnitude of seismic anisotropy are controlled by the interference of strongly anisotropic phyllosilicates with less anisotropic but more abundant quartz and plagioclase. The amphibolite-facies shear zone margin (> 500 °C) is free of major mineral transformations, and rather characterized by the development of S-C structure, inducing an orthorhombic symmetry to the VP distribution and up to 12 % AVP, with the fast axis parallel to lineation. AVS is maximized (6 %) in the two directions parallel to the C-planes. Further deformation under greenschist-facies (450–500 °C) induces significant crystallization of muscovite, chlorite, and albite. A dominant hexagonal symmetry develops, with a fast VP plane and a fast S-wave polarization plane along the C-planes. AVP and AVS exceed 10 %. The computed elastic tensor of the whole shear zone emphasizes the importance of the strain gradient on the seismic signature of the overall structure. In our model, the bulk seismic properties are dominated by the orthorhombic symmetry of the protomylonitic zone rather than by the ultramylonite core with hexagonal symmetry (C-plane), which may explain why seismic studies of shear zones commonly involve a fast velocity plane oblique to the layer boundaries.

Pn tomography and anisotropic study of the Indian shield and the adjacent regions

High-resolution P-wave velocity and anisotropy structure of the hitherto elusive uppermost mantle beneath the Indian shield and its surrounding regions are presented to unravel the tectonic imprints in the lithosphere. We inverted high quality 19,500 regional Pn phases from 172 seismological stations for 4780 earthquakes at a distance range of 2° to 15° with a mean apparent Pn velocity of 8.22 km/s. The results suggest that the Pn velocity anomalies with fast anisotropic directions are consistent with the collision environments in the Himalaya, Tibetan Plateau, Tarim Basin, and Burmese arc regions. The higher Pn anomalies along the Himalayan arc explicate the subducting cold Indian lithosphere. The cratonic upper mantle of the Indian shield is characterized by Pn velocity of 8.12–8.42 km/s, while the large part of the central Indian shield has higher mantle-lid velocity of ~8.42 km/s with dominant anisotropic value of 0.2–0.3 km/s (~7.5%) suggesting the presence of mafic ‘lava pillow’ related to the Deccan volcanism. The impressions of the rifts and the mobile belts are conspicuous in the velocity anomaly image indicating their deep seated origin. The Pn anisotropy in the Indian shield exhibits a complex pattern and deviates from the absolute plate motion directions derived from the SKS study, demonstrating the presence of frozen anisotropy in the Indian lithospheric uppermost mantle, due to the large scale tectonic deformation after its breakup from the Gondwanaland. Whereas, Pn and SKS anisotropic observations are well consistent in Tarim basin, Tibetan regions, eastern Himalayan syntaxis and the Burmese arc. The modeled anisotropic Pn clearly manifests a lower velocity anomaly bounded by 85°E and 90°E ridges in the southern Bay of Bengal. Further, 85°E ridge spatially separates the BoB lithosphere into faster and slower regions consistent with the body wave tomography and free-air gravity observation.

Fracture Shape and Orientation Contributions to P-Wave Velocity and Anisotropy of Alpine Fault Mylonites

P-wave anisotropy is significant in the mylonitic Alpine Fault shear zone. Mineral- and texture-induced anisotropy are dominant in these rocks but further complicated by the presence of fractures. Electron back-scattered diffraction and synchrotron X-ray microtomography (micro-CT) data are acquired on exhumed schist, protomylonite, mylonite, and ultramylonite samples to quantify mineral phases, crystal preferred orientations, microfractures, and porosity. The samples are composed of quartz, plagioclase, mica and accessory garnet, and contain 3–5% porosity. Based on the micro-CT data, the representative pore shape has an aspect ratio of 5:2:1. Two numerical models are compared to calculate the velocity of fractured rocks: a 2D wave propagation model, and a differential effective medium model (3D). The results from both models have comparable pore-free fast and slow velocities of 6.5 and 5.5 km/s, respectively. Introducing 5% fractures with 5:2:1 aspect ratio, oriented with the longest axes parallel to foliation decreases these velocities to 6.3 and 5.0 km/s, respectively. Adding both randomly oriented and foliation-parallel fractures hinders the anisotropy increase with fracture volume. The anisotropy becomes independent of porosity when 80% of fractures are randomly oriented. Modeled anisotropy in 2D and 3D are different for similar fracture aspect ratios, being 30 and 15%, respectively. This discrepancy is the result of the underlying assumptions and limitations. Our numerical results explain the effects that fracture orientations and shapes have on previously published field- and laboratory-based studies. Through this numerical study, we show how mica-dominated, pore-free P-wave anisotropy compares to that of fracture volume, shape and orientation for protolith and shear zone rocks of the Alpine Fault.

Lattice-Preferred Orientation and Seismic Anisotropy of Minerals in Retrograded Eclogites from Xitieshan, Northwestern China, and Implications for Seismic Reflectance of Rocks in the Subduction Zone

Various rock phases, including those in subducting slabs, impact seismic anisotropy in subduction zones. The seismic velocity and anisotropy of rocks are strongly affected by the lattice-preferred orientation (LPO) of minerals; this was measured in retrograded eclogites from Xitieshan, northwest China, to understand the seismic velocity, anisotropy, and seismic reflectance of the upper part of the subducting slab. For omphacite, an S-type LPO was observed in three samples. For amphibole, the <001> axes were aligned subparallel to the lineation, and the (010) poles were aligned subnormal to foliation. The LPOs of amphibole and omphacite were similar in most samples. The misorientation angle between amphibole and neighboring omphacite was small, and a lack of intracrystalline deformation features was observed in the amphibole. This indicates that the LPO of amphibole was formed by the topotactic growth of amphibole during retrogression of eclogites. The P-wave anisotropy of amphibole in retrograded eclogites was large (approximately 3.7–7.3%). The seismic properties of retrograded eclogites and amphibole were similar, indicating that the seismic properties of retrograded eclogites are strongly affected by the amphibole LPO. The contact boundary between serpentinized peridotites and retrograded eclogites showed a high reflection coefficient, indicating that a reflected seismic wave can be easily detected at this boundary.

Seismic Anisotropy in Subduction Zones: Evaluating the Role of Chloritoid

Subduction zones are often characterized by the presence of strong trench-parallel seismic anisotropy and large delay times. Hydrous minerals, owing to their large elastic anisotropy and strong lattice preferred orientations (LPOs), are often invoked to explain these observations. However, the elasticity and the LPO of chloritoid, which is one of such hydrous phases relevant in subduction zone settings, are poorly understood. In this study, we measured the LPO of polycrystalline chloritoid in natural rock samples, obtained the LPO-induced seismic anisotropy, and evaluated the thermodynamic stability field of chloritoid in subduction zones. The LPO of chloritoid aggregates displayed a strong alignment of the [001] axes subnormal to the rock foliation, with a girdle distribution of the [100] axes and the (010) poles subparallel to the foliation. New elasticity data of single-crystal chloritoid showed a strong elastic anisotropy of chloritoid with 47% for S-waves (VS) and 22% for P-waves (VP), respectively. The combination of the LPO and the elastic anisotropy of the chloritoid aggregates produced a strong S-wave anisotropy with a maximum AVS of 18% and a P-wave anisotropy with an AVP of 10%. The role of chloritoid LPO in seismic anisotropy was evaluated in natural rock samples and a hypothetical blueschist. Our results indicate that the strong LPO of chloritoid along the subduction interface and in subducting slabs can influence the trench-parallel seismic anisotropy in subduction zones with “cold” geotherms.

Deformation Microstructures of Phyllite in Gunsan, Korea, and Implications for Seismic Anisotropy in Continental Crust

Muscovite is a major constituent mineral in the continental crust that exhibits very strong seismic anisotropy. Muscovite alignment in rocks can significantly affect the magnitude and symmetry of seismic anisotropy. In this study, deformation microstructures of muscovite-quartz phyllites from the Geumseongri Formation in Gunsan, Korea, were studied to investigate the relationship between muscovite and chlorite fabrics in strongly deformed rocks and the seismic anisotropy observed in the continental crust. The [001] axes of muscovite and chlorite were strongly aligned subnormal to the foliation, while the [100] and [010] axes were aligned subparallel to the foliation. The distribution of quartz c-axes indicates activation of the basal<a>, rhomb<a> and prism<a> slip systems. For albite, most samples showed (001) or (010) poles aligned subnormal to the foliation. The calculated seismic anisotropies based on the lattice preferred orientation and modal compositions were in the range of 9.0–21.7% for the P-wave anisotropy and 9.6–24.2% for the maximum S-wave anisotropy. Our results indicate that the modal composition and alignment of muscovite and chlorite significantly affect the magnitude and symmetry of seismic anisotropy. It was found that the coexistence of muscovite and chlorite contributes to seismic anisotropy constructively when their [001] axes are aligned in the same direction.

Crystal preferred orientations, deformation mechanisms and seismic properties of high pressure metamorphic rocks from the Central Qiangtang metamorphic belt, Tibetan Plateau

Crystal preferred orientation (CPO) of high pressure metamorphic rocks is an important contributor to seismic anisotropy in the subduction interface and continental crust. However, the mechanisms of deformation and CPO development, as well as the effect of CPO patterns on the seismic anisotropies remain not fully solved. By analyzing the microstructures, CPOs and seismic properties of the high pressure rocks from the Central Qiangtang metamorphic belt in the northern Tibetan Plateau, we found that the mineral CPOs (e.g., omphacite and amphibole) are developed dominantly by non-dislocation based deformation mechanisms in the eclogites and amphibolites, while quartz and phengite CPOs are primarily formed by dislocation creep in the garnet-quartz-mica schists. The variations of CPO patterns are mainly controlled by the strain geometry and/or strain magnitude, which also regulate the first-order symmetry patterns of seismic anisotropies among samples. Besides, amphibole enrichment can remarkably increase the anisotropies and decrease the velocities of seismic waves in the eclogite‒amphibolite suite. A combination of seismic parameters such as Vp, Vs, Vp/Vs ratio, intensities and symmetries of P-wave anisotropy (AVp) and shear wave splitting (AVs), and P- or S-wave reflection coefficients can be used to discriminate eclogites, amphibolites and their garnet-quartz-mica schists country rocks at depth when high resolution seismic methods can be applied. Notably, amphibolite and garnet-quartz-mica schists showing moderate-to-large AVp and max. AVs can make considerable contributions to field-observed crustal anisotropy near the strike-slip continental shear zones or mantle wedge anisotropy in the subduction zones.

The elastic properties and anisotropy of artificial compacted clay samples

Clay minerals are a major component of hydrocarbon reservoir rocks, and they are known to play important roles in the physical and elastic properties of rocks. However, it is difficult to directly measure these properties of single-crystal clays due to their small particle size. Therefore, we have constructed three sets of artificial clay samples with different compaction stresses to investigate the effect of the compaction stress and clay mineralogy on their elastic properties and anisotropy. All of the dry samples are measured by the pulse-transmission method. The results indicate that the compaction stress and clay mineralogy have a significant influence on the physical and elastic properties of the clay samples. The microstructures of clay samples indicate that the clay platelets are aligned almost perpendicularly to the direction of compaction stress, and the ultrasonic velocity analysis validates the assumption of transverse isotropy of our clay samples. The velocities increase with the compaction stress, especially at low stress, which corresponds to the rapid porosity reduction at low stress levels. Velocity anisotropy parameters increase with increasing of compaction stress due to the increase of texture sharpness for clay minerals during the compaction process. The elastic moduli of the clay samples display a significant stress sensitivity and a strong directional dependence, with the Young’s moduli increasing and the Poisson’s ratios decreasing with the compaction stress. A simple theoretical template is used to quantify the orientation distribution functions (ODFs) of clay platelets, and the generalized Legendre coefficients of ODF increase with the increase of compaction stress, especially at low stress. Further, the compressional-wave (P-wave) and shear-wave anisotropy increase with the ODF coefficients [Formula: see text] and [Formula: see text], especially P-wave anisotropy.

Stress-path-dependent effective medium model for granular media — Comparison with experimental data

Modeling velocity changes in response to stress changes plays an important part in understanding seismic responses from the subsurface. One branch of such modeling consists of treating an assemblage of grains as an effective medium and using established grain contact theories to determine the elastic moduli. Such models are commonly limited to hydrostatic or uniaxial strain scenarios, not capable of capturing the anisotropy induced by a general triaxial stress state. A new set of expressions is developed by extending an existing effective medium theory to a general stress state in which the radial (horizontal) stress is different from the axial (vertical) stress. The theory is valid at the limits of slip and no slip. Novel functions to combine the no-slip and slip limits are implemented to match the model to observed laboratory data on glass beads and sand grain assemblages loaded along different stress paths. The new expressions provide a good agreement between the modeled and measured stress and stress path dependence of the compressional wave (P-wave) and shear wave velocities and associated P-wave anisotropy. This stress dependence is of particular interest in rock-physics modeling workflows evaluating time-lapse feasibility for shallow or unconsolidated sand reservoirs or for characterizing burial history.

Unexpected Consequences of Transverse Isotropy

ABSTRACTIn a series of articles, Kawakatsu et al. (2015) and Kawakatsu (2016a,b, 2018) introduced and discussed a new parameter, ηκ, that characterizes the incidence angle dependence (relative to the symmetry axis) of seismic body-wave velocities in a transverse isotropy (TI) system. During the course of these exercises, several nontrivial consequences of TI were realized and summarized as follows: (1) P-wave velocity (anisotropy) strongly influences the conversion efficiency of P-to-S and S-to-P, as much as S-wave velocity perturbation does; (2) Rayleigh-wave phase velocity has substantial sensitivity to P-wave anisotropy near the surface; (3) a trade-off exists between ηκ and the VP/VS ratio if the latter is sought under an assumption of isotropy or the elliptic condition. Among these findings, the first two deserve careful attention in interpretation of results of popular seismic analysis methods, such as receiver function analysis and ambient-noise Rayleigh-wave dispersion analysis. We present simple example cases for such problems to delineate the effect in actual situations, as well as scalings among TI parameters of the crust and mantle materials or models that might help understanding to what extent the effect becomes important.