P-wave Anisotropy Research Articles

New insight into the velocity and anisotropy structures of the subduction zone in northern Sumatra

In this study, we conducted seismic tomographic inversions to investigate the velocity and anisotropy structures of northern Sumatra, using 9774 P-wave and 8405 S-wave arrivals from regional earthquakes. Isotropic P-wave velocity, isotropic S-wave velocity, P-wave azimuthal anisotropy, and P-wave radial anisotropy models were generated using eikonal equation-based traveltime tomography methods. The study identified low-velocity zones beneath the Toba and Sinabung volcanoes, potentially indicating the presence of magma reservoirs. Furthermore, low-velocity anomalies above the subduction slab were detected, which are likely caused by the dehydration of the slab and interpreted as channels of upwelling flow. The tomographic results revealed a trench-parallel high-velocity belt in the uppermost mantle, representing the subducting slab of the India-Australian plate. The trench-parallel fast velocity directions in the slab suggested that the subducted oceanic slab retains its frozen-in anisotropy formed at the mid-ocean ridge, or that the anisotropy is induced by the lattice-preferred orientation of the B-type olivine. Negative radial anisotropy in the mantle wedge was observed, reflecting hot upwelling flows and transitions of olivine fabrics in the presence of water due to slab dehydration. The results also indicated a multilevel magma plumbing system beneath the Toba Caldera. In summary, the results of this study provide new insights into the structure and dynamic processes of the northern Sumatra subduction zone.

An iterative weighted least-square fitting method for crustal anisotropy using receiver functions

SUMMARY The harmonic variation of the P-to-S converted phases (i.e. Pms) observed from receiver functions (RFs) includes information on crustal azimuthal anisotropy. However, this harmonic analysis is easily influenced by low-quality RF traces, and the measurements may be misleading. Here, we propose an improved method, named the iterative weighted least-square method (IWLS), to extract the splitting parameters of the crust and simultaneously retrieve the two- and four-lobed components of backazimuthal variation. The quality and weights of different RF traces are estimated properly in the IWLS method. The weight function is related to the sharpness of the Pms phase and the smearing of other signals. We conduct many synthetic tests, and the IWLS method provides stable measurements for poor backazimuthal coverage, strong noise, weak P-wave azimuthal anisotropy and multiple anisotropic layers. We apply the IWLS method to observational data from two temporary stations on the southeastern Tibetan Plateau and North China Craton, respectively. The measurements are comparable to previous results and provide insight into crustal deformation.

DEM Analysis of the P-Wave Anisotropy Response to the Microstructure of Sedimentary Rock Under Biaxial Compression

DEM Analysis of the P-Wave Anisotropy Response to the Microstructure of Sedimentary Rock Under Biaxial Compression

An anisotropic rock physics modeling for the coalbed methane reservoirs and its applications in anisotropy parameter prediction

The elastic properties and anisotropy of coalbed methane reservoirs are still poorly understood due to the complicated pore structure and the existing state of the coalbed methane. Therefore, an anisotropic rock physics model for the coalbed methane reservoirs is proposed based on the related effective medium theories, which focuses on the equivalent calculation of the elastic properties for adsorbed coalbed methane and the characterization of the multiple pore structure. Many factors related to the elastic properties of coal are considered, including the composition, pore structure, adsorbed gas content and pore fluid. The measured ultrasonic and well-logging data demonstrate that the proposed physics model is effective in predicting the elastic constants and velocity anisotropy parameters. The results suggest that the elastic constants increase with the increase of adsorbed gas content to some extent, and the effect of porosity on the elastic constants is much more obvious than that of adsorbed gas content. The velocity anisotropy parameters increase with the aligned fracture porosity and the fracture aspect ratio, and the P-wave anisotropy is more sensitive to these two factors. With increasing organic matter and clay content and porosity, Young's modulus and brittleness index decrease while the Poisson's ratio increases, and the impact of the porosity on the brittleness index is more significant than that of organic matter and clay content. The results are helpful in establishing the relationship between reservoir physics parameters and seismic response and can provide a basis for coalbed methane reservoir prediction and hydraulic fracturing.

Deformation microstructures of blueschists in Alpine Corsica, France, and implications for seismic anisotropy and the low-velocity layer in subducting oceanic crust

Understanding the deformation microstructures and seismic properties of blueschists is important for interpreting the seismic anisotropy at the top of subducting oceanic crust. The deformation microstructures and seismic properties of epidote blueschist, epidote-lawsonite blueschist, and lawsonite blueschist from Alpine Corsica, France, were studied using electron backscatter diffraction mapping and transmission electron microscopy. The crystal-preferred orientations (CPOs) of glaucophane show a strong maximum of the [001] axes aligned sub-parallel to the lineation. The maximum of the (010) poles of epidote is aligned sub-parallel to the lineation. The [001] axes of lawsonite are strongly aligned sub-normal to the foliation. Intracrystalline misorientation, dislocation, and occasional fragmentation of glaucophane indicate that glaucophane CPO was developed mainly by dislocation creep. The intracrystalline deformation features of epidote, such as subgrain boundaries and dislocations, indicate that the CPO of the epidote was formed by dislocation creep. In contrast, our data indicate that lawsonite CPO was developed by a combination of rigid-body rotation and dislocation creep with a minor effect of twinning. The P-wave anisotropy (AVp) and maximum S-wave anisotropy (max.AVs) of epidote blueschists (AVp = 12.3–16.9% and max.AVs = 6.53–11.75%) are higher than those of lawsonite blueschists. The seismic velocities of the blueschists are lower than those of mantle peridotites. The coexistence of epidote and lawsonite in the blueschist results in variations in seismic anisotropy and P-wave velocity in the whole rock, depending on the modal content of epidote and lawsonite. A high content ratio of epidote to lawsonite in blueschist produces strong P- and S-wave anisotropies of blueschist when the glaucophane content is more than ∼50%. The P-wave velocity of blueschist decreases proportionally with increasing epidote content. Our results indicate that the seismic properties of deformed blueschists contribute to the seismic anisotropy and the low-velocity layer in subducting oceanic crusts in subduction zones.

Deformation Behavior and Seismic Characteristics of Sandy Facies Opalinus Clay During Triaxial Deformation Under Dry and Wet Conditions

Unconsolidated, undrained triaxial deformation tests were performed on sandy facies Opalinus Clay at 50 MPa confining pressure to characterize the effect of water and microfabric orientation on the deformation behavior, mechanical properties, and P-wave velocity evolution. Dry and wet (≈ 8 and > 95% initial water saturation, respectively) samples with 12.6 ± 0.4 vol% porosity were deformed parallel and perpendicular to the bedding direction at a constant strain rate of 5 × 10–6 s−1. Dry samples revealed semi-brittle behavior and exhibited strain localization at failure, while deformation was more ductile at saturated conditions, promoting stable, slow faulting. Peak strength, Young’s modulus, and number of cumulative acoustic emissions decreased significantly for wet samples compared to dry samples; the opposite was observed for Poisson’s ratio. P-wave velocity anisotropy was significantly altered by differential stress, primarily due to the interplay between pore and fracture closure and stress-induced microcrack formation. For samples that were deformed perpendicular to bedding, we observed a reduction and reversal of P-wave velocity anisotropy with increasing differential stress, whereas anisotropy of parallel samples increased. The results suggest that water saturation reduces the pressure at the brittle-ductile transition and that the elastic properties and anisotropy of sandy facies Opalinus Clay can be significantly altered in an anisotropic stress field, e.g., adjacent to fault zones or tunnel excavations. Changes in elastic anisotropy are primarily controlled by the orientation between the pre-existing microfabric and the maximum principal stress direction, stress magnitude, and the degree of water saturation.

Remnants of shifting early Cenozoic Pacific lower mantle flow imaged beneath the Philippine Sea Plate

Seismic anisotropy could provide vital information about the evolution and internal convection of the deep Earth interior. Although previous seismological studies have revealed a wide distribution of seismic anisotropy in the upper portion of the lower mantle beneath many subduction zones, the existence of anisotropy at these depths away from subducted slabs remains debated. Here we use P-wave azimuthal anisotropy tomography to image the crust and mantle down to 1,600-km depth. We find prominent anisotropic patterns in the upper portion of the lower mantle beneath the Philippine Sea Plate. Substantial azimuthal anisotropy with N–S fast-velocity directions occurs at 700–900-km depths. We interpret this azimuthal anisotropy as a remnant of the Pacific lower mantle flow field about 50 million years ago. Two isolated high-velocity anomalies at 700–1,600-km depths may be vestigial pieces of the subducted Izanagi slab with seismic velocity features suggesting a shift in the Pacific lower mantle flow field by about 40 million years ago. Our findings provide seismic evidence for the existence of complex lower mantle flows and deformation mechanisms away from subduction zones.

Variation in olivine crystal-fabrics and their seismic anisotropies in the Horoman Peridotite Complex, Hokkaido, Japan

We examined the microstructures and crystal-fabrics of peridotites within a large area (6×5km) of the Horoman Peridotite Complex in the Hidaka metamorphic belt of Hokkaido, Japan. Thirteen peridotite samples were analyzed for olivine and orthopyroxene grain sizes, fabric strength (J-index), and crystallographic preferred orientations (CPOs). Mean grain sizes of olivine and orthopyroxene were ranged in 295–497µm and in 257–537µm, respectively. The olivine fabric strength values decreased from the lower to the upper part of the complex, whereas the orthopyroxene fabric strength values showed no systematic trends. The peridotites contained three different olivine CPOs, previously known as A, E, and AG types. Combined with a previous study, we found that olivine CPOs showed a transitional distribution from E to A to AG type from south to north. E type peridotites occur at the basement of the complex in the south, suggesting that local water infiltration might occur at the basement of the complex. The A type peridotites occurred mainly in the middle of the studied area and subsequently the AG type peridotites occurred towards the north. Moreover, we calculated the seismic properties of peridotite as olivine 100% aggregates and mixed (olivine and orthopyroxene) aggregates. It showed that orthopyroxene CPOs reduce P-wave anisotropies of peridotite (0.2–2.2%) without modification of the P-wave propagation patterns.

Experimental Investigation of Stress Sensitivity of Elastic Wave Velocities for Anisotropic Shale in Wufeng–Longmaxi Formation

The shale of the Wufeng–Longmaxi formation in the Sichuan Basin is the preferred layer for shale gas exploration in China, and its petrophysical characteristics are the key to geological and engineering sweet spot prediction. However, the characteristics and impact mechanisms of its acoustic wave velocity and elastic anisotropy are currently unclear. In this paper, the Wufeng–Longmaxi shale is taken as the research object, and the P-wave and S-wave velocities of the samples are tested under the loading and unloading processes of confining pressure. The stress sensitivity variations in parameters such as wave velocity, wave velocity ratio, and anisotropy are discussed. P-wave and S-wave anisotropy parameters are correlated under different pressure conditions. X-ray diffraction, casting thin sections, scanning electron microscopy, micron CT scanning, and other analytical techniques are used to explore the mechanisms of stress sensitivity of elastic parameters. The research results indicate that: (1) the acoustic velocities of samples from different angles are V90° > V45° > V0°, and there is a positive correlation between the wave velocity and the confining pressure. After unloading the confining pressure, irreversible plastic deformation occurs due to the closure of some microfractures in the rock core, causing the wave velocity to be higher than the initial value. (2) The stress sensitivity coefficient of the P-wave (The mean is 3.00 m·s−1·MPa−1) is higher than that of the S-wave (the mean is 1.23 m·s−1·MPa−1), and the stress sensitivity coefficient of the compacted stage (the mean is 3.02 m·s−1·MPa−1) is higher than that of the elastic stage (the mean is 1.21 m·s−1·MPa−1). (3) The anisotropy of the P-wave and S-wave is negatively correlated with the confining pressure. When the confining pressure is loaded to 65 MPa, the change rate of the P-wave anisotropy coefficient is 23%, and its stress sensitivity is higher than that of S-wave anisotropy coefficient (the change rate is 13.7%). After unloading the confining pressure, the degree of anisotropy is reduced due to the closure of some microfractures. The empirical formula of P-wave and S-wave anisotropy parameters under different pressures is established through linear regression, which can provide a reference for mutual predictions. (4) The variation in wave velocity anisotropy with stress can be divided into stress and material anisotropy, which are related to the directional arrangement of microfractures and clay minerals, respectively. The quantitative characterization of shale anisotropy can be realized by evaluating the development degree of reservoir fractures and mineral components, providing a reference for logging interpretations, sweet spot prediction, and fracturing construction of shale gas reservoirs.

Experimental Validation and Calibration of the Galvin Model with Artificial Tight Sandstones with Controlled Fractures

The study of fractures in the subsurface is very important in unconventional reservoirs since they are the main conduits for hydrocarbon flow. For this reason, a variety of equivalent medium theories have been proposed for the estimation of fracture and fluid properties within reservoir rocks. Recently, the Galvin model has been put forward to model the frequency-dependent elastic moduli in fractured porous rocks and has been widely used to research seismic wave propagation in fractured rocks. We experimentally investigated the feasibility of applying the Galvin model in fractured tight stones. For this proposal, three artificial fractured tight sandstone samples with the same background porosity (11.7% ± 1.2%) but different fracture densities of 0.00, 0.0312, and 0.0624 were manufactured. The fracture thickness was 0.06 mm and the fracture diameter was 3 mm in all the fractured samples. Ultrasonic P- and S-wave velocities were measured at 0.5 MHz in a laboratory setting in dry and water-saturated conditions in directions at 0°, 45°, and 90° to the fracture normal. The results were compared with theoretical predictions of the Galvin model. The comparison showed that model predictions significantly underestimated P- and S- wave velocities as well as P-wave anisotropy in water-saturated conditions, but overestimated P-wave anisotropy in dry conditions. By analyzing the differences between the measured results and theoretical predictions, we modified the Galvin model by adding the squirt flow mechanism to it and used the Thomsen model to obtain the elastic moduli in high- and low-frequency limits. The modified model predictions showed good fits with the measured results. To the best of our knowledge, this is the first study to validate and calibrate the frequency-dependent equivalent medium theories in tight fractured rocks experimentally.

Crustal Anisotropy from the Birefringence of P-to-S Converted Waves: Bias Associated with P-Wave Anisotropy

Many researchers have used the birefringence of P‑to‑S converted waves from the Moho discontinuity to constrain the anisotropy of Earth’s crust. However, this practice ignores the substantial influence that anisotropy has on the initial amplitude of the converted wave, which adds to the splitting acquired during its propagation from Moho to the seismometer. We find that large variations in Ps birefringence estimates with back-azimuth occur theoretically in the presence of P‑wave anisotropy, which normally accompanies S‑wave anisotropy. The variations are largest for crustal anisotropy with a tilted axis of symmetry, a geometry that is often neglected in birefringence interpretations, but is commonly found in Earth’s crust. We simulated globally-distributed P‑coda datasets for 36 distinct 4‑layer crustal models with combinations of elliptical shear anisotropy or compressional anisotropy, and also incorporated the higher-order anisotropic Backus parameter C. We tested both horizontal and tilted symmetry-axis geometries and tested the birefringence tradeoff associated with Ps converted phases at the top and bottom of a thin high‑ or low‑velocity basal layer. We computed composite receiver functions (RFs) with harmonic regression over back azimuth, using multipletaper correlation with moveout corrections for the epicentral distances of 471 events, to simulate a realistic data set. We estimate Ps birefringence from the radial and transverse RFs, a strategy that is similar to previous studies. We find that Ps splitting can be a useful indicator of bulk crustal anisotropy only under restricted circumstance, either in media with no compressional anisotropy, or if the symmetry axis is horizontal throughout. In other, more-realistic cases, the inferred fast polarization of Ps birefringence estimated from synthetic RFs tends either to drift with back-azimuth, form weak penalty-function minima, or return splitting times that depend on the thickness of an anisotropic layer, rather than the birefringence accumulated within it.

Acoustic and mechanical tests of sandstone-shale composites in Songliao Basin and prediction of uniaxial compressive strength

Rock uniaxial compressive strength is a vital rock mechanics parameter, which accurate assessment is critical for high-efficiency unconventional reservoir development. However, its realization for standard cores is quite problematic due to difficulties and high costs related to the collection and preparation of samples. Unconventional fracturing treatment have evolved from single to composite lithology. This study collected shale and sandstone samples at different bedding angles (0°, 45°, 90°)to prepare the sandstone-shale composite samples at different interlayer ratios. It measured the prepared composites’ longitudinal wave and transverse wave velocities with ultrasonic velocity testing. The uniaxial compressive strength of the sandstone-shale composites with various bedding angles were tested via a servo rock triaxial testing machine. The anisotropic characteristics of the prepared sandstone-shale composites with various bedding angles were analyzed, and relation between the uniaxial compressive strength of the anisotropic sandstone-shale composite and P-wave velocity was derived and experimentally validated. The results show that: the wave velocity test results of the sandstone-shale combination show that with the increasing proportion of shale, the wave velocity of the longitudinal wave and the transverse wave both increase; The wave velocity results of different bedding planes in shale are as follows: Vp(0°)> Vp(45°)> Vp(90°). The average P-wave and S-wave velocities at 45° are 0.92, while the average P-wave and S-wave velocities at 90° are 0.85 compared to 0°; the P-wave and S-wave anisotropy of shale are both 0.14, and the shale shows significant anisotropy. In the sandstone-shale composite, the uniaxial compressive strength gradually decreases and the strain increases with the increase of shale proportion. The value of uniaxial compressive strength is UCS(0°) > UCS(90°) > UCS(45°). Based on the P-wave modulus, a fitting relationship formula for the uniaxial compressive strength of sandstone-shale composites with different bedding angles and interlayer ratios was established and verified. The anisotropic strength predictive formula established based on experimental data is simple, precise, and easy to obtain, with strong practicability. The present research results can provide proper technical support for sandstone-shale composite fracturing optimization and reservoir development.

Lattice-preferred orientation (LPO) of olivine and amphibole in amphibole peridotites and neighboring hornblendites from Gapyeong, South Korea and implications for seismic anisotropy

Amphibole peridotites and neighboring hornblendites are often found in subduction zones. To understand the effect of amphibole-rich rocks on seismic anisotropy in subduction zones, we studied the lattice-preferred orientation (LPO) of olivine and amphibole in amphibole peridotites and neighboring hornblendites in Gapyeong, South Korea. The major minerals of amphibole peridotites were olivine (31–51% in volume), amphibole (28–47%), and orthopyroxene (7–16%). Amphibole in amphibole peridotites showed relatively high SiO2 and low TiO2 and Na2O contents (S-type amphibole), indicating that it was formed under supra-subduction conditions. Amphibole in amphibole peridotites showed the type-I, type-II, and type-IV LPO, whereas amphibole in neighboring hornblendites showed the type-III and type-IV LPO. In the case of olivine, most samples showed a mixture of A- and B-type LPO, and one sample showed a mixture of B- and C-type LPO. Many serpentine inclusions were found in olivines. Fourier-transform infrared (FTIR) analysis of the samples showed that the olivines contained a large amount of water (∼29000–45000 ppm H/Si). We also found many dislocations in olivines. These observations indicate that samples showing a mixture of A- and B-type LPO and a mixture of B- and C-type LPO of olivine were deformed under water-rich conditions by dislocation creep. In amphibole peridotites, the P-wave anisotropy of olivine was relatively low (0.9–4.8%), whereas the P-wave anisotropy of amphibole was high (6.5–17.7%). The maximum S-wave anisotropy of olivine was also relatively low (0.87–2.89%), whereas the maximum S-wave anisotropy of amphibole was high (3.81–15.19%). In hornblendites, the P-wave anisotropy and maximum S-wave anisotropy of amphibole were high (6.9–13.6% and 4.27–10.61%, respectively). The P-wave anisotropy and maximum S-wave anisotropy of the amphibole peridotites were in the range of 2.6–8.4% and 1.73–7.30%, respectively. The seismic velocity and anisotropy pattern of amphibole peridotites were more similar to those of amphibole than those of olivine, indicating that the seismic properties of amphibole peridotites were more strongly affected by amphibole than by olivine. Furthermore, the seismic anisotropy of the mixture of amphibole peridotite and hornblendite in subduction zones was also found to be significantly affected by the amphibole LPO in hornblendite.

PS reflection moveout in a homogeneous anisotropic layer of arbitrary symmetry and tilt

Use of unconverted reflected PP waves leads to moveout formulae, which relate the moveout to only a limited number of parameters specifying anisotropy. Use of converted reflected PS waves extends the number of parameters to, generally, a complete set of twenty-one parameters. It leads, on the other hand, to more complicated formulae and phenomena unknown from studies of unconverted PP waves. In anisotropic media one has to deal with the existence of two S waves and their singularities, acoustic approximation cannot be used. We present and test approximate formulae for the reflection moveout of PS waves converted at a horizontal reflector underlying a homogeneous layer of arbitrary anisotropy symmetry and orientation. For their derivation, we use the combination of weak-anisotropy approximation and A-parameters specifying anisotropy. The complete set of twenty-one A-parameters represents an alternative to twenty-one independent elements of the stiffness tensor. In addition to the moveout formulae for separate PS1 and PS2 waves, we also present moveout formula for a P to common S wave whose traveltimes represent an average of traveltimes of PS1 and PS2 waves. Use of the P to common S wave reflection leads to several simplifications. Presented tests indicate that for most offsets and models with P-wave anisotropy up-to to 26% and S-wave anisotropy over 25%, the relative traveltime errors do not exceed 2.5%.

Interpreting the effects of shale rock properties on seismic anisotropy by statistical and machine learning methods

The elastic anisotropy of unconventional shale reservoirs plays a crucial role in their economic hydrocarbon production. The strong anisotropy of organic-rich shale rocks, characterized by clay platelets, elongated pores, and organic material, presents challenges for traditional reservoir characterization methods, which are designed for isotropic or weakly anisotropic conventional sandstone reservoirs. To effectively describe the geomechanical behaviors of organic-rich shales and estimate changes in stress caused by hydraulic fracturing, it is necessary to develop practical methods for predicting anisotropy parameters from conventional rock property logs. This study proposes a statistical approach for predicting independent elastic stiffness coefficients for vertically transverse isotropic (VTI) shale from three shale rock properties: clay volume fraction, kerogen volume fraction, and total porosity. The results of this study show that the stiffness coefficients, C11 and C33, are related to P-wave anisotropy and porosity has a greater impact on P-wave anisotropy than clay and kerogen volume fractions. In contrast, changes in clay volume have a more sensitive effect on C44 and C66, providing S-wave anisotropy information that can help detect lithology changes. Furthermore, kerogen variation has a greater influence on both C12 and C13 than the effects of porosity and clay variations. From these relationships, we develop exponential nonlinear regression (ENLR) models for predicting independent elastic stiffness coefficients from three shale rock properties. These models are validated with the Midland Basin log data, and we improve the accuracy of the C13 estimation model using the deep feedforward neural network (DFNN) method. The results show that the statistical approach can help understand the effect of shale rock components on seismic anisotropy variations and better predict the variations in lithology, fluid, and organic contents in shales from anisotropic seismic data. Furthermore, the methodology can provide robust estimates of the stress changes from hydraulic fracturing, supporting the efficient and economically viable production of shale hydrocarbons.

Geochemical features and seismic imaging of the tectonic zone between the Tibetan Plateau and Ordos Block, central northern China

The Tibetan Plateau is growing by both vertical uplift and horizontal extension. It is a continuing debate how the Tibetan Plateau interacts with its surrounding plates and blocks. Due to intense tectonic activity, which produced catastrophic earthquakes, the tectonic zone between the northeast margin of the horizontal extending Tibetan Plateau and the stable Ordos Block has garnered considerable interest. This study investigated the spatial distribution of gas geochemical anomalies (e.g., high flux of CO2 in correspondence of the main faults) at regional scale together with the seismic tomography in correspondence of this tectonic zone with the aim to figure out the domain of convergent boundary between the Ordos block and Tibetan plateau, and trace the tectonic discontinuities which are able to transfer fluids through the crustal layers between the two main geological units. From northwest to southeast, obvious difference of spatial distributions of geochemical and geophysical features in the tectonic zone between the northeast margin of the Tibetan Plateau and the Ordos Block is inferred. The northeast area (Zone A) is dominated by thrust and strike-slip faults with clear velocity boundary underneath, where low crack density (ε), saturation rate (ξ) and Poisson’ ratio (σ) in the middle-lower crust coincided with the low values of heat flow and CO2 emissions, tectonic compression and regional locked-fault can be inducements. The southeast area (Zone C) is dominated by extensional tectonics with roughly E-W fast-velocity direction (FVD) of P-wave azimuthal anisotropy, where high permeability and porosity can be deduced from crustal high ε, ξ and relatively high σ anomalies, resulting in high heat flow, CO2 concentrations and fluxes at the surface, and predominantly crustal-derived gases. The intermediate area (Zone B) also dominated by thrust and strike-slip faults is an extraordinary zone, where intensely locked-fault were clearly revealed, while the predominant anisotropic FVDs in the middle crust changed obviously, more contribution of shallow gas component was detected, and CO2 flux, heat flow, and regional ε, ξ, and σ in the upper crust were higher, compared with those in Zone A, which indicated the regional crushing fragmentation underneath Zone B. The adopted multidisciplinary approach demonstrated that Zone B is the convergent boundary between the Tibetan Plateau and the Ordos Block.

Waveform inversion of large data sets for radially anisotropic Earth structure

SUMMARY The amount of high-quality seismic data is expanding rapidly, and there is a need for algorithms that take advantage of classical methods to achieve high efficiency using widely available computing power. In this study, we develop a novel waveform inversion method to retrieve radially anisotropic Earth models that can be used to investigate deformation and flow in the mantle. Our method is comprised of two parts: (1) extraction and fitting of the fundamental mode and (2) fitting of the full synthetic waveform. The waveform inversion method results in path average model constraints with uniquely determined independent uncertainties. We demonstrate through synthetic testing that the method is able to retrieve radially anisotropic perturbations down to the mantle transition zone, and leakage effects due to ignoring P-wave anisotropy are minimal. We apply the method to ∼16 000 waveforms generated by earthquakes occurring in the East Sea (Sea of Japan) region, and we demonstrate that the subsequent linear inversion of radially anisotropic path constraints produces models that are similar to those resulting from full waveform adjoint tomography methods. We validate our model by predicting waveforms for earthquakes not included in our inversion, and we show that our method is able to extract structural information. Our results indicate low-velocity anomalies and weak radial anisotropy in NE Japan, which may be due to competing influences from ascending fluids and/or melts and horizontal flow in the lower crust and upper mantle. In the southern East Sea, we image low velocities and relatively high radial anisotropy, which may reflect high temperatures, shallow dehydration and olivine LPO in the upper mantle.

Laboratory measurements of ultrasonic wave velocities and anisotropy across a gold-hosting structure: A case study of the Thunderbox Gold Mine, Western Australia

Most lode gold deposits worldwide are associated with structures such as shear zones. Thanks to their capacity to couple resolution and depth of investigation, seismic methods can identify these indirect indicators of mineralization and help extend gold exploration targets to greater depths. Rocks from shear zones are usually seismically anisotropic. Seismic anisotropy is generally related to the intrinsic texture of the rock and the presence of cracks at depth, although the impact of the latter relative to the former fades out with increasing depth. Understanding the seismic impedance contrast between the rocks and determining seismic anisotropy in relation to the texture of the rock, and its evolution with depth (pressure) is therefore necessary to help interpret exploratory seismic surveys. To do so, we characterized in the laboratory 126 core samples representing different stages of alteration and different lithologies from the Thunderbox Gold Mine in Western Australia. Laboratory measurements consist of ultrasonic wave velocity at ambient and in situ pressure, density, mineralogy, texture analysis, and estimation of P-wave anisotropy by inverting the ultrasonic data. These experimental data were used as input in different rock physics models to calculate texture-derived velocities and to model the effects of seismic anisotropy on seismic reflectivity. Except for massive basalt samples that present a remarkable elastic impedance contrast with all lithologies and for the talc schist samples (shear zone samples) that present high seismic anisotropy values, the seismic reflectivity between the lithological units encountered at the Thunderbox Gold Mine is low. According to seismic reflectivity modelling, seismic anisotropy along the shear zone significantly affects the seismic reflectivity, with values of reflectivity reducing more with the increasing angle of seismic reflection than if the shear zone is considered as an isotropic interface. The good agreement between calculated texture-derived velocities with experimental measurements shows that the crystallographic preferred orientation of minerals of the shear zone samples is the main source of seismic anisotropy. This study seeks to improve the understanding of the seismic response across mineral deposits that are structurally controlled by shear zones.

P-wave velocity anisotropy in Hamra Quartzites reservoir, Hassi Messaoud oil field in Algeria

P-wave velocity anisotropy in Hamra Quartzites reservoir, Hassi Messaoud oil field in Algeria

Pressure-dependent elastic anisotropy: A Bakken petroleum system case study

The complexity of the micropore fabric and organic texture in tight unconventional reservoirs often complicates their correlations with physical rock property parameters. Anisotropy, which is a common physical property of tight oil reservoirs, can affect directional deformation and permeability during the reservoir development and pressure depletion period. Thus, accurately measuring seismic anisotropy and its relationships with petrophysical properties can refine our understanding of vertical transverse isotropy rocks and ultimately improve fracturing design as well as production forecasting. We have developed a measurement technique of ultrasonic velocities (P and S) from multiple directions, that is, bedding-parallel, bedding-perpendicular, and at 45° to bedding under simulated in situ reservoir pressure. Following these measurements, we also correlate these findings to pressure-dependent porosity and permeability, in addition to geologic attributes and textural anisotropy. The samples that we investigate are tight carbonaceous siltstone and organic-rich mudrock cores originating from the Bakken petroleum system, which include Lodgepole, upper Bakken shale, middle Bakken, lower Bakken shale, and Three Forks. Depending on the rock fabric and mineral constituents, each formation exhibits distinctive anisotropic behavior in response to pressure changes. The results indicate that the elastic parameters and petrophysical properties are inversely related under pressure for all samples. We also find that P-wave anisotropy and SH-wave anisotropy parameters decrease gradually with increasing confining pressure. The anisotropic behavior can be related to volumetric concentration and directional connectivity of clay layers, kerogen veins, and subparallel microfractures.