Wave Velocity Anisotropy Research Articles

Insight into wave propagation in polyimide films and resistive grid sandwich structures towards a hybrid monitoring of hypervelocity impact

Micro-Meteoroid and Orbital Debris pose a significant threat to the safe operation of orbiting spacecraft, potentially leading to mission failure in space exploration. Quantitative characterization of hypervelocity impact (HVI) is crucial to ensure the safety and successful completion of on-orbit missions. Firstly, this study designed a three-layer sandwich structure of polyimide film with orthogonally laid resistive wires, combined with piezoelectric and resistive wire sensors, for the simultaneous acquisition of acoustic emission (AE) signals generated by HVI and measurement of perforation dimensions. Secondly, a semi-analytical finite element (SAFE) analysis of wave dispersion properties in the periodic sandwich structure is conducted with Bloch’s theorem, together with a hybrid model based on three-dimensional smoothed particle hydrodynamics and finite element methods (SPH-FEM) to comprehensively understand the AE waves and damage characteristics induced by HVI. The resulting anisotropic wave propagation characteristics with SAFE and SPH-FEM are closely matched. Thirdly, a time delay-multiplication (TDM) imaging algorithm considering wave velocity anisotropy is proposed for accurate real-time “visualization” of HVI locations. Lastly, correlations are established between projectile and perforation dimensions. The proposed algorithm for HVI multi-parameter quantification and damage detection helps evaluate the space HVI environment and HVI-induced damage to spacecraft.

Guided wave multi-frequency damage localization method in variable-thickness structures by one pair of sensors based on frequency-dependent velocity anisotropy

Variable thickness structures are prevalent in aircraft, ships, and other machines, necessitating numerous sensors for health monitoring to reduce safety hazards. This paper presents a guided wave multi-frequency localization method based on frequency-dependent velocity anisotropy. This method achieves damage localization in variable-thickness structures with a pair of sensors and can effectively reduce the number of sensors used for monitoring. Variations in structural thickness cause a gradient in guided wave velocity that bends the propagation path. Different thickness variations with different directions cause wave velocity anisotropy. As a result, variations in thickness cause possible damage loci determined by echo time to deviate from an elliptical shape. Because the velocity anisotropy is frequency-dependent, damage loci at different frequencies are close but do not overlap and intersect only at the damage location. So, the multi-frequency method can increase the damage information acquired by a single pair of sensors, enabling damage localization. Experimental validation was conducted on a steel plate with linearly varying thicknesses. The feasibility of the multi-frequency localization method was verified by successfully locating the damage at three different locations using a pair of receiver-excitation sensors. In addition, the experiments demonstrated the capability of this multi-frequency method in improving the localization accuracy of sensor networks. The method has potential applications in monitoring systems lightweight, phased arrays, and imaging enhancement.

3D shear wave velocity and azimuthal anisotropy structure in the shallow crust of Binchuan Basin in Yunnan, Southwest China, from ambient noise tomography

The Binchuan Basin in northwest Yunnan, southwest China, is a rift basin developed at the intersection of the Red River Fault and Chenghai Fault, where historical earthquakes have occurred. Understanding the fine velocity structure of the shallow crust in this region can help improve earthquake location accuracy and our understanding of the relationship between fault zone structures and fault slip behaviors. Using the continuous waveform data recorded by 381 dense array stations in 2017, we obtained 7 915 Rayleigh-wave phase velocity dispersion curves in the period band of 0.2–6 ​s from ambient noise cross-correlation functions after rigorous data processing and quality control. We determined 3D isotropic and azimuthally anisotropic shear wave velocity models at depths above 6 ​km in the shallow crust based on the direct surface wave azimuthal anisotropic tomography method. The isotropic model reveals a strong correspondence between the S-wave velocity structure at depths of 0–1 ​km and the regional topography and lithology. The Binchuan depocenter, Zhoucheng depocenter, Xiangyun Basin, and Xihai Rift Basin are primarily composed of Quaternary deposits, which show low-velocity anomalies, while the regions with the Paleozoic shale, limestone, and basalt exhibit high-velocity anomalies. The nearly N–S orientation of fast directions from azimuthal anisotropy models are mainly controlled by the active Binchuan Fault with N–S strike as well as the NNW-oriented primary compressive stress.

Orientation-dependent multi-spall performance of monocrystalline NiTi alloys under shock compression

The primary objective of this work was committed to large-scale non-equilibrium molecular dynamics modeling for the investigation on multiple spall performance, microstructural evolutions, and melting state of NiTi alloys along three typical crystal orientations. We discovered that there was a single-wave structure present in the [100] orientation, whereas elastic-plastic two-wave structures were found in [110] and [111], exhibiting a significant anisotropy in elastic-plastic wave velocity. Compared to [110] and [111] orientations, multi-spallation was more likely to be generated in the [100] direction. During accumulative dynamic damage, there was also a strong sensitivity of void evolution to orientation. For [100], voids tended to be densely distributed in the material and exhibited a strip-shaped manner. In contrast, a more scattered void distribution was observed for the other two orientations, showing a flaky pattern. Interestingly, dislocations in the spall region were almost annihilated before approaching void nucleation. According to void statistics, the peak number of voids, void nucleation rate, and spall region size were identified to be in a good uniformity along [110], [100], and [111]. For all the cases, the samples experienced partial melting in the released state and were followed by a recovery of structural order degree. However, there was an insignificant anisotropy in the melting state for each crystal orientation during multi-spallation.

Shear wave anisotropy response characteristics of laminated shale and its correction methods

Laminated shale exhibits strong anisotropic characteristics in horizontal wells, significantly affecting the application of shear wave travel time data. Using equivalent medium theory and the elastic wave equation, we constructed a physical model for laminated shale, simulating the impact of lamination angle and lamination density on shear wave velocity and shear wave anisotropy coefficient. We established the relationship between the shear wave anisotropy coefficient and the cosine of the lamination angle. The simulation results indicate that shear wave velocity decreases with an increase in lamination angle and density, while the shear wave anisotropy coefficient gradually decreases with an increase in lamination angle and initially increases and then decreases with an increase in lamination density. Moreover, it exhibits a power exponential relationship with the cosine of the lamination angle. After correction, the minimum horizontal principal stress and fracture pressure calculated from shear wave data show an average relative error reduction of 6.72% and 7.77%, respectively, compared to measured data. This demonstrates that the correction method effectively eliminates the impact of laminated shale anisotropy on shear waves.

Phase diagram and physical properties anisotropy of strontianite

The physical properties and anisotropy of carbonates at high pressure are fundamental to understanding the deep carbon cycle and storage. Here, we redraw the phase diagram of strontianite (SrCO3) by first-principles calculations combined with lattice dynamics and study its physical properties and anisotropy at high pressure. The results show that there are only Pnma and Pmmn phases in SrCO3 in the entire Earth's mantle, and there is no R3m phase. The SrCO3 polymorphs have a large elastic anisotropy. SrCO3 polymorphs are a kind of carbonate with high density, low wave velocity, and low wave velocity anisotropy. The linear thermal expansion coefficient of SrCO3-Pmmn was obtained for the first time by using the linear elastic compression coefficient combined with the volume thermal expansion coefficient. The results are of great significance for understanding the high-pressure behavior of Sr-carbonates.

Crustal shear wave velocity and radial anisotropy beneath the Mississippi embayment from ambient noise tomography

Imaging of shallow multiple basins and middle to lower crustal structure beneath the Mississippi embayment is still limited. No tomography studies have delineated multiple basins as effectively as those from reflection and drilling experiments. Several studies have constructed different velocity models for the middle to lower crust. We take advantage of the Northern Embayment Lithosphere Experiment dataset and long-term deployment of 277 3-component broadband stations installed from 1990 to 2018 to image crustal shear wave velocity and radial anisotropy in the Mississippi embayment. We first construct 5–40s group velocity maps for Rayleigh wave and 5–30s ones for Love wave. In the 5 s group velocity map, we observe that low velocity features correlate spatially with the sediment thickness from the Crust 1.0 model. High velocity anomalies are associated with the Ozark uplift and the Nashville Dome. In the 32 s group velocity map, high-low velocity boundaries align with 40 km crustal thickness contours from the Crust 1.0 model. We then invert Rayleigh and Love group velocity dispersion for shear wave velocity and radial anisotropy at each grid point and construct a 3D model accordingly. At the depth of 0–5 km of the model, we find that low velocity features characterize the Mississippi Valley Graben, Rough Creek Graben, Black Warrior Basin, and coastal flood plain. High velocity features mark thinly sedimented regions, like the Ozark uplift and the Nashville Dome. Regions with high gravity anomalies are characterized with negative anisotropies. At the depth of 5–12 km, we observe low velocity features in the Missouri Batholith, centering around 10 km. This region is characterized with negative radial anisotropy (Vsh < Vsv), inferring possible intrusions. High-low velocity boundaries correlate with the Appalachian-Ouachita thrust front. At the depth of 12–24 km, the Mississippi Valley Graben is filled with high velocity material, which may be related to the formation of the graben due to isostatic equilibrium. At the depth of 24–43 km, high velocity anomalies characterize the lower crust of the Mississippi embayment, consistent with high density features from the inversion of free air anomaly. The high velocity anomaly might be responsible for the formation of the embayment due to isostatic adjustment. The average velocity in the lower crust is approximately 4.3 km/s. Mafic Pikwitonei granulite is a possible material, indicating relative strong lower crust. This work provides useful information to explain earthquake occurrence in the New Madrid Seismic Zone and regional tectonic settings.

Spatial distribution of mid-lower crustal flow in the SE Tibetan Plateau revealed by P-wave velocity and azimuthal anisotropy beneath the Lijiang–Xiaojinhe fault and its vicinity

SUMMARY The Lijiang–Xiaojinhe fault (LXF) and its vicinity are located in the transition zone among the Tibetan Plateau (TP), the South China block and the Indochina block. Researchers believe that this area has acted as a key tectonic zone during the evolution of the TP. Owing to the continuous growth and SE-ward expansion of the TP, the LXF and its vicinity have experienced intense deformation. Although different models, such as the rigid block extrusion and mid-lower crustal flow models, have been proposed to explain this intense deformation, a consensus has not yet been achieved. To better understand the deformation of the LXF and its vicinity, a high-resolution image of the subsurface structure must be constructed. In this study, we construct images of P-wave velocity and azimuthal anisotropy structures by using an eikonal equation-based traveltime tomography method. We collect high-quality seismic data from 276 broad-band seismic stations and manually pick a total of 48 037 first arrivals for the tomography study. Our tomographic results reveal a strong low-velocity body below the LXF and its vicinity. In addition, a strong azimuthal anisotropy structure with an N–S-oriented fast velocity direction is distributed along the low-velocity body. These features indicate the occurrence of mid-lower crustal flow, that penetrates across the LXF and extends to the Dianzhong block (DZB). In addition, we find obvious low-velocity perturbations in the mid-lower crust and uppermost mantle beneath the DZB. The low velocities may be attributed to the upwelling of hot materials from the upper mantle. We consider the limited distribution of mid-lower crustal flow on the margin of the SE TP, and mid-lower crustal flow may not play a significant role in the expansion of the TP.

Anisotropy of Consolidation Yield Stress and Comparison with Anisotropy of Elastic Wave Velocity and Tensile Strength in Sedimentary Soft Rocks

Anisotropy of Consolidation Yield Stress and Comparison with Anisotropy of Elastic Wave Velocity and Tensile Strength in Sedimentary Soft Rocks

3-D crustal shear wave velocity and azimuthally anisotropic structure of the Pearl River Delta onshore-offshore area and its tectonic implications

The Pearl River Delta (PRD) onshore-offshore area lies in the transition zone between South China and northern continental margin of the South China Sea. To better understand the tectonic and deformation patterns at depth, we developed a 3-D high-resolution crustal shear wave velocity (Vs) and azimuthal anisotropy model by ambient noise tomography in this region. The data sets comprise ∼30 days of seismic ambient noise recorded by 124 land seismic stations and 28 Ocean Bottom Seismometers (OBS). Our results reveal a significant lateral contrast in azimuthal anisotropy in the crust. The fast axes of the azimuthal anisotropy in land show an approximately NW-SE direction, which may be mainly controlled by the regional stress field. In the offshore PRD, however, the fast axes of anisotropy strike nearly in the NE-SW direction, likely dominated by the shearing of the littoral fault zone. Combined with other previous results, the anisotropic model reveals three major crustal deformation processes in this region since the Yanshanian. Furthermore, an upper-middle crustal frozen mafic magma chamber with high Vs (∼3.85–4.0 km/s) is imaged beneath Shenzhen, Zhuhai, and Hong Kong. It is interpreted as a source of intense magmatic activities in and around Hong Kong during the late Mesozoic. Our new model provides important constraints for understanding the crustal architecture and the multistage deformation processes in and around the PRD onshore-offshore area.

Mapping of stress and structure controlled upper crustal anisotropy in Kumaon-Gharwal Himalaya

Shear wave splitting analysis of local earthquakes provides valuable insights into the structure and stress controlled upper crustal anisotropy signatures. We examine these signatures at the Kumaon-Garhwal Himalaya and provides an excellent example of anisotropic crustal structure in the ongoing continent-continent collision settings between Indian and the Eurasian plates. A total of 256 local earthquakes were selected for the analysis (1.0 ≤ M ≤ 5.4) between January 2017 and February 2021, recorded at a dense network of 51 broadband stations. The result indicates the observed crustal anisotropy is parallel to stress-aligned micro-cracks far from the major Himalayan fault zones and structure parallel near the fault zones. The dominant fast polarization directions (FPD) in the Inner Lesser Himalaya (ILH) and Higher Himalaya, are consistent with the maximum compressive horizontal stress directions (SHmax), essentially influenced by the local stress field. Whereas in the Outer Lesser Himalaya (OLH) and Sub Himalaya, fast directions are sub-parallel to the structural trends, suggesting that the anisotropy is associated with the shear fabric from recent deformation episodes related to structural induced anisotropy. The computed normalized delay times show a mean value of 1 ms/km, within the OLH and Sub-Himalaya, while in ILH and Higher Himalaya, this value increases up to 5.3 ms/km. The associated crack densities are 0.0038 and 0.0207, with shear wave velocity anisotropy of 0.38% and 2.07% respectively. Significant scatter within the depth range of 10–15 km is observed in normalized delay times, suggesting the source of anisotropy within the upper crust.

Evolution law of ultrasonic characteristics and its relationship with coal-measure sandstone mechanical properties during saturation and desaturation

Underground buildings are easily destabilised by water and there is a strong link between their mechanical properties and ultrasonic characteristics. To grasp the above features and the relationship among them, the mechanism of sandstone degradation under water erosion was briefly reviewed. Using the coal-measure sandstone as the object, experimental studies on the dynamic evolution of water content, ultrasonic wave velocity, waveform, uniaxial compressive strength (UCS), elastic modulus (E), and damage mode during saturation-desaturation were carried out. The findings demonstrate that coal-measure sandstones usually exhibit more mechanical parameter deterioration following saturation than non-coal-measure sandstones. During saturation, the water absorption velocity slows down; the P-wave velocity (VP) first fluctuates down, then fluctuates up and finally decreases slightly (VP increases after saturation), and the wave velocity anisotropy coefficient (δi) decreases; the waveform shows a trend of unequal time and unequal amplitude → equal time and equal amplitude → unequal time and equal amplitude, and the first wave reaches a progressively shorter time. During desaturation, the dehydration rate slows down, VP first drops rapidly and then fluctuates slightly and rises slightly near drying (VP decreases after desaturation), δi increases, and the waveform changes in the opposite trend to saturation. Saturation treatment significantly weakened the mechanical properties, while desaturation restored some of the strength lost due to immersion. The sandstone mechanical properties and ultrasonic features are more closely related during saturation than desaturation.

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° &gt; V45° &gt; 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.

Lamina Influences on Tensile Strength of Shallow Marine Shales from Upper Ordovician, Western Ordos Basin.

To investigate the influence of the lamina effect on the tensile strength of shallow marine shales and improve the shortcomings of the existing Brazilian standard disc splitting method, the Lamina shale of the Upper Ordovician Wulalik Fm in the western margin of the Ordos Basin was selected as the experimental object. Based on the radial wave velocity anisotropy test, the direction of crustal stress was determined, and standard cores were drilled. The Brazilian standard disc splitting experiment on Lamina shale with different loading angles was designed and carried out. The influence of lamina on the tensile strength of shale was summarized, and an improved calculation method of tensile strength was proposed. The experimental results indicate that the presence of lamina makes the tensile strength of shallow marine shale exhibit significant anisotropy, and the fracture surface morphology of standard discs under different loading angles varies greatly. The overall failure characteristics can be classified into two types: linear and curved. When the loading angle is 0° or 90°, the fracture surface of the disc belongs to tensile failure (linear type), and the traditional splitting method has good applicability. When the loading angle is greater than 0° and less than 90°, the fracture surface of the disc belongs to tensile shear failure (curve type), and traditional splitting methods are not applicable. There is a difference in tensile strength between vertical and horizontal wells, and vertical wells should consider the comprehensive tensile strength of the rock matrix and lamina at a 90° loading angle. Horizontal wells should consider the tensile strength of the weak lamina plane with a loading angle of 0°. The improved Brazilian splitting method solves the problem of the traditional method, calculating lower tensile strength values when the loading angle is greater than 0° and less than 90°. This provides important basic data support for wellbore stability evaluation and reservoir stimulation transformation.

First-principles study on the high-pressure physical properties of orthocarbonate Ca2CO4

Orthorhombic Ca2CO4 is a recently discovered orthocarbonate whose high-pressure physical properties are critical for understanding the deep carbon cycle. Here, we study the structure, elastic and seismic properties of Ca2CO4-Pnma at 20–140 GPa using first-principles calculations, and compare them with the results of CaCO3 polymorphs. The results show that the structural parameters of Ca2CO4-Pnma are in good agreement with the experimental results. It could be the potential host of carbon in the Earth's mantle subduction slab, and its low wave velocity and small anisotropy may be the reason why it cannot be detected in seismic observation. The thermodynamic properties of Ca2CO4-Pnma at high temperature and high pressure are obtained using the quasi-harmonic approximation method. This study is helpful in understanding the behavior of Ca-carbonate in the Earth’s lower mantle conditions.

Analytic formula for the group velocity and its derivatives with respect to elastic moduli for a general anisotropic medium

Studying the anisotropy of seismic wave velocity can provide a theoretical basis for understanding the characteristics of seismic wave propagation in anisotropic media, such as forward modeling of seismic wavefield and simultaneous inversion of anisotropic parameters. It is an important work to calculate the group velocities of three kinds of seismic waves (qP, qS1, and qS2) and the analytical expressions of the partial derivatives of the group velocities with respect to 21 elastic parameters in general anisotropic media. Therefore, a new relationship between the phase and group velocity in the general anisotropic medium is introduced first, and then the analytical expressions of the group velocity of three kinds of waves are deduced. Then, the analytical formulas of partial derivatives of the group (or phase) velocity with respect to 21 elastic parameters are derived. Finally, the distribution of partial derivatives of group slowness with respect to 21 elastic parameters with varied ray angles is analyzed and discussed. The final analytical expressions of group velocity and its partial derivatives are relatively concise and can deal with the singular points and multivalued problems of qS waves. By drawing comparison to the numerical solution of the eigenvector method, the accuracy of the analytical formula of group velocity and its partial derivatives with respect to elastic parameters are verified. The numerical simulation results indicate that the algorithm can accurately obtain the partial derivatives of group slowness with respect to 21 elastic parameters and would provide a theoretical basis for seismic ray tracing and traveltime inversion in general anisotropic media.

Evaluation and engineering applications of the in situ stress state of deep tight sandstone reservoirs in the Xujiahe Formation of the Yingshan-Longgang in the central Sichuan Basin

Introduction: The Xujiahe Formation in the Yingshan-Longgang area of the Sichuan Basin has abundant tight oil and gas resources. However, the complex geological conditions and insufficient understanding of the magnitude and direction of in-situ stresses severely limit exploration and development of oil and gas resources in the region. This has made it difficult to carry out hydraulic fracturing and horizontal well deployment in the region.Methods: Based on experimental analyses, including acoustic emission experiments, wave velocity anisotropy experiments, and paleomagnetic experiments, this study evaluated the magnitude and direction of in-situ stresses, while also combining electrical imaging logging and acoustic wave logging data. This provides a basis for hydraulic fracturing in the research area.Results and Discussion: The experimental results indicate that the three-way stress state of the Xujiahe Formation is a strike-slip stress state, and the three-way stress value gradually increases with increasing burial depth. The Longgang area is particularly affected by the depth. The average azimuth angle of the maximum horizontal principal stress in the study area is N100°E, while tectonic deformation may cause a counterclockwise deviation of the maximum horizontal principal stress azimuth in the Yingshan area. Considering the influence of changes in the magnitude of in situ stress in the longitudinal direction on hydraulic fracturing extension mode, it is recommended to select the middle and upper sections of the Xu II for fracturing and transformation. Similarly, considering the distribution of natural fractures and the direction of the in situ stresses, it is suggested that the deployment azimuth angle for horizontal wells be between N20°–30°E.

Lithospheric deformation in Mongolia revealed by anisotropic Rayleigh wave tomography: Implications for the mechanism of intracontinental orogeny

The Mongolian Plateau is situated in the eastern part of the Central Asian Orogenic Belt that was formed during the Paleozoic era. In the Cenozoic era, the plateau has been reactivated, leading to significant topographic uplift, volcanism, and seismic activity. Several destructive earthquakes, ranging from magnitudes 6 to 8, have occurred in the 20th century. But the mechanism of these neotectonics remains unclear. To better understand these processes, we determined S wave velocity and azimuthal anisotropy in the crust and upper mantle under Mongolia using Rayleigh wave dispersion curves. Our findings indicate that the lithosphere is dominated by NW-SE fast velocity directions (FVDs). However, we observed different patterns in two specific regions. The lithospheric anisotropy under the Hangay Dome is weak, suggesting that it is accompanied by significant mantle upwelling under the dome. On the other hand, the FVDs in the Gobi Desert exhibit a circular pattern that wraps around the Ordos Block, which may have been caused by past and present deformation imposed on the stable block. These results offer new insights into the lithospheric deformation and dynamics of the Mongolian Plateau.

Quantifying fabric anisotropy of granular materials using wave velocity anisotropy: a numerical investigation

Fabric anisotropy is a sought-after micro index to correlate macro mechanical responses of granular materials. In this work, the discrete-element method is utilised to simulate multi-directional bender element tests in granular soils to obtain the evolution of wave velocities during drained conventional triaxial (CT) and true triaxial (TT) tests; the contact normal-based fabric is simultaneously monitored for bridging the fabric anisotropy and wave velocity anisotropy. The results show that stress-normalised wave velocities and microscopic fabric, including contact normal distribution and coordination number, remain nearly constant until a stress ratio threshold is reached. After the threshold value is reached, stress-normalised wave velocities start to decrease, especially in the minor principal stress direction, accompanied by significant adjustment of coordination number and fabric anisotropy. The results also reveal that the normalised wave velocity depends on the contact normal densities in the wave propagation and particle oscillation directions. With the contact normal distribution represented by a density function, a good linear relationship between the microscopic fabric anisotropy and macroscopic wave velocity anisotropy is obtained for both CT and TT tests.

Deformation of the Qinling belt revealed byP-wave velocity and azimuthal anisotropy tomography

SUMMARYThe Qinling belt is a transitional zone lying among three units: the North China block (NCB), the South China block (SCB) and the northeastern Tibetan Plateau (NETP). Owing to the interaction of these units, complex deformation has occurred in the Qinling belt. Although many studies have been conducted to understand the deformation mechanism in the Qinling belt, some key issues are still under debate, such as whether middle-lower crustal flow exists beneath the western Qinling belt (WQB). High-resolution images of subsurface structures are essential to shed light on the deformation mechanism. In this paper, high-resolution images of the velocity structure and azimuthal anisotropy beneath the Qinling belt are obtained by using an eikonal equation-based traveltime tomography method. Our seismic tomography inverts 38 719 high-quality P-wave first arrivals from 1697 regional earthquakes recorded by 387 broad-band seismic stations. In the WQB, our tomography results show low-velocity anomalies but relatively weak anisotropy in the middle-lower crust. These features suggest that middle-lower crustal flow may not exist in this area. In the central Qinling belt (CQB), we find low-velocity anomalies in the middle-lower crust; however, the fast velocity directions no longer trend E–W but vary from NNE–SSW to N–S. These characteristics can be ascribed to the convergence and collision between the NCB and the SCB. In addition, we find strong low-velocity anomalies in the uppermost mantle beneath the CQB, which may indicate delamination of the lower crust. In the southern Qinling belt, we observe significant high-velocity anomalies in the upper crust beneath the Hannan–Micang and Shennong–Huangling domes. These high-velocity anomalies indicate a mechanically strong upper crust, which is responsible for the arc-shaped deformation process of the Dabashan fold. Based on the P-wave velocity and azimuthal anisotropic structures revealed by the inversion of high-quality seismic data, the deformation of the Qinling belt is affected mainly by the convergence between the NCB and the SCB rather than by the middle-lower crustal flow from the Tibetan Plateau.