P-wave Tomography in Inhomogeneous Orthorhombic Media

— A P-wave tomographic method for 3-D complex media (3-D distribution of elastic parameters and curved interfaces) with orthorhombic symmetry is presented in this paper. The technique uses an iterative linear approach to the nonlinear travel-time inversion problem. The hypothesis of orthorhombic anisotropy and 3-D inhomogeneity increases the set of parameters describing the model dramatically compared to the isotropic case. Assuming a Factorized Anisotropic Inhomogeneous (FAI) medium and weak anisotropy, we solve the forward problem by a perturbation approach. We use a finite element approach in which the FAI medium is divided into a set of elements with polynomial elastic parameter distributions. Inside each element, analytical expressions for rays and travel times, valid to first-order, are given for P waves in orthorhombic inhomogeneous media. More complex media can be modeled by introducing interfaces separating FAI media with different elastic properties. Simple formulae are given for the Fréchet derivatives of the travel time with respect to the elastic parameters and the interface parameters. In the weak anisotropy hypothesis the P-wave travel times are sensitive only to a subset of the orthorhombic parameters: the six P-wave elastic parameters and the three Euler angles defining the orientation of the mirror planes of symmetry. The P-wave travel times are inverted by minimizing in terms of least-squares the misfit between the observed and calculated travel times. The solution is approached using a Singular Value Decomposition (SVD). The stability of the inversion is ensured by making use of suitable a priori information and/or by applying regularization. The technique is applied to two synthetic data sets, simulating simple Vertical Seismic Profile (VSP) experiments. The examples demonstrate the necessity of good 3-D ray coverage when considering complex anisotropic symmetry.     

忽略TTI介质对称轴倾角的可行性

假设横向各向同性(TI)介质的对称轴是垂直的(VTI)或者水平的(HTI)能给实际资料处理带来便利,然而实际TI介质的对称轴往往是倾斜的(TTI),忽略对称轴倾角可能给各向异性参数提取和成像带来偏差,因此需要研究是否能、以及什么条件下能忽略TTI介质对称轴倾角.本文通过理论研究和数值分析研究了与TTI介质弹性性质最接近的VTI介质(OAVTI)的弹性常数和各向异性参数与原TTI介质的弹性常数和各向异性参数之间的联系与差别.结果表明:OAVTI介质各向异性参数与原TTI介质各向异性参数之间的差别可统一表示成F(α₀/β₀,ε,δ,γ)ξ²的形式,其中F(α₀/β₀,ε,δ,γ)是无量纲各向异性参数(ε, δ, γ)的线性函数,ξ是对称轴倾角;ξ的大小对各参数的误差起主导作用,一般不建议忽略20°~25°以上的对称轴倾角;当ξ较小时,即使是对强各向异性的TTI介质作VTI近似,引起的P波各向异性参数误差也很小,因此在纵波资料处理中忽略TTI介质对称轴倾角通常是可行的;即使在小ξ条件下,倾斜对称轴对SV波也有显著影响,因此在转换波资料处理中,不建议忽略TTI介质的对称轴倾角.本文的研究为分析忽略TTI介质对称轴倾角的可行性提供了理论依据和简便的判据.     

Seismic, Structural and Petrological Models of the Subcrustal Lithosphere in Southern Germany: A Quantitative Revaluation

—Anisotropy in the subcontinental lithosphere becomes increasingly important, because it is observed in many seismic studies especially for P _n -waves. Typical rocks of the uppermost mantle are peridotites, which predominantly exhibit a pronounced elastic anisotropy. This anisotropy is mainly caused by the anisotropic elastic properties and the lattice preferred orientation (here referred to as texture) of olivine. To evaluate the elastic anisotropy of peridotites from the subcontinental lithosphere, specimens of the Northern Hessian Depression (Germany) and the Balmuccia Ultramafic Massif (Northern Italy) have been used. They comprise four olivine texture types, which are characteristic for olivine textures observed worldwide. The bulk rock elastic properties have been calculated using olivine and orthopyroxene textures, their single-crystal elastic constants at ambient pressure/temperature conditions and their volume fraction. Clinopyroxene and spinel are assumed to be randomly distributed. The effect of four different orientations of the foliation within the uppermost mantle has been evaluated, since this orientation is usually unknown.¶Two of the olivine textures have a pronounced azimuthal dependence of compressional waves when a horizontal foliation within the uppermost mantle is presumed. These variations cause significant azimuthal variations of the P-wave reflections coefficients at the Moho. Primarily, we predict a significant azimuthal dependence of the critical points where the reflected amplitude increases from approximately 15% to 95%. Possibly, these azimuthal variations can be detected by seismic reflection measurements carried out at earth surface.¶The remaining two texture types only manifest a small directional dependence. When anisotropy of compressional waves is observed in seismic studies, these latter types can only be of subordinate importance. However, all of the peridotites investigated are able to explain the seismically observed azimuthal variations of compressional waves when a vertical foliation is proposed. This ambiguity can be substantially reduced when shear waves (S-waves) are considered. The directional distribution of S-wave velocities and of the S-wave splitting exhibits characteristic patterns for the different olivine texture types. This could be used to discriminate between different texture types and orientations of the foliation within the uppermost mantle. A fundamental requirement for a more comprehensive interpretation is the availability of detailed S-wave observations. The maximum S-wave splitting in the peridotites investigated coincides with the maximum of the faster (leading) S-wave. This may be of importance to detect S-wave splitting in future seismic studies.     

Elastic anisotropy of ferromagnesian post-perovskite in Earth's D” layer

We report experimental observation of a sizable elastic anisotropy in a polycrystalline sample of ferromagnesian silicate in post-perovskite (ppv) structure. Using a novel composite X-ray transparent gasket to contain and synthesize ppv in a panoramic diamond-anvil cell along with oblique X-ray diffraction geometry, we observed the anisotropic lattice strain and {1 0 0} or {1 1 0} slip-plane texture in the sample at 140 GPa. We deduced the elasticity tensor (c_ij), orientation-dependent compressional wave velocities, polarization-dependent shear-wave velocities, and the velocity anisotropy of the silicate ppv. Our results are consistent with calculations and indicate that with sufficient preferred orientation, the elastic anisotropy of this phase can produce large shear-wave splitting.     

Search for orthorhombic or higher symmetry in rocks

Summary The properties of quadrics of stress and strain, which were derived from the eigenvectors of the matrix of elastic constants, are analysed for a medium with orthorhombic or higher symmetry. It was found that the orientation of the axes of the quadrics can be employed to determine the crystallographic axes of the medium. This formalism derived for minerals is used for finding the axes of the internal symmetry of a rock for which we assume a distribution close to that of the orthorhombic system with a view to the spatial distribution of P-wave velocities.     

Classes of Anisotropic Media: A Tutorial

The purpose of the present article is to give a precise definition and analysis from first principals of anisotropy, as the term applies to elastic media, taking care to avoid unnecessary assumptions. Two fundamental concepts, material invariance and symmetry group of a material, are defined purely in terms of the stress-strain relation. The implications of material symmetry, or in other words, of anisotropy, for the structure of the stiffness tensor are then investigated. Using the reduced notation of Voigt, these results are presented as the well-known simplifications in the form taken by the six-by-six stiffness matrix that represents the material's stiffness tensor. A new, simple proof is given for the remarkable fact that an elastic medium cannot have rotational symmetry by an angle of less than 90° without being transversely isotropic. In addition, the mutual relation that the notions of elastic symmetry and crystal symmetry have with respect to the so-called orthogonal group is sketched. Despite the historical association between anisotropic elastic materials and the study of crystals, the given presentation shows that conceptually the notion of anisotropy in elastic media is entirely independent of that of crystal symmetry.     

On Choosing Effective Elasticity Tensors Using a Monte-Carlo Method

A generally anisotropic elasticity tensor can be related to its closest counterparts in various symmetry classes. We refer to these counterparts as effective tensors in these classes. In finding effective tensors, we do not assume a priori orientations of their symmetry planes and axes. Knowledge of orientations of Hookean solids allows us to infer properties of materials represented by these solids. Obtaining orientations and parameter values of effective tensors is a highly nonlinear process involving finding absolute minima for orthogonal projections under all three-dimensional rotations. Given the standard deviations of the components of a generally anisotropic tensor, we examine the influence of measurement errors on the properties of effective tensors. We use a global optimization method to generate thousands of realizations of a generally anisotropic tensor, subject to errors. Using this optimization, we perform a Monte Carlo analysis of distances between that tensor and its counterparts in different symmetry classes, as well as of their orientations and elasticity parameters.     

Reference transversely isotropic medium approximating a given generally anisotropic medium

For a given stiffness tensor (tensor of elastic moduli) of a generally anisotropic medium, we can estimate the extent to which the medium is transversely isotropic, and determine the direction of its reference symmetry axis. In this paper, we rotate the given stiffness tensor about this reference symmetry axis, and determine the reference transversely isotropic (uniaxial) stiffness tensor as the average of the rotated stiffness tensor over all angles of rotation. The obtained reference transversely isotropic (uniaxial) stiffness tensor represents an analytically differentiable approximation of the given generally anisotropic stiffness tensor. The proposed analytic method is compared with a previous numerical method in two numerical examples.     

Simulation of full-waveform log in saturated cracked formations using Hudson's approach

We study the propagation of elastic waves that are generated in a fluid‐filled borehole surrounded by a cracked transversely isotropic medium. In the model studied the anisotropy and borehole axes coincide. To obtain the effective elastic moduli of a cracked medium we have applied Hudson's theory that enables the determination of the overall properties as a function of the crack orientation in relation to the symmetry axis of the anisotropic medium. This theory takes into account the hydrodynamic mechanism of the elastic‐wave attenuation caused by fluid filtration from the cracks into a porous matrix. We have simulated the full waveforms generated by an impulse source of finite length placed on the borehole axis. The kinematic and dynamic parameters of the compressional, shear and Stoneley waves as functions of the matrix permeability, crack orientation and porosity were studied. The modelling results demonstrated the influence of the crack‐system parameters (orientation and porosity) on the velocities and amplitudes of all wave types. The horizontally orientated cracks result in maximal decrease of the elastic‐wave parameters (velocities and amplitudes). Based on the fact that the shear‐ and Stoneley‐wave velocities in a transversely isotropic medium are determined by different shear moduli, we demonstrate the feasibility of the acoustic log to identify formations with close to horizontal crack orientations.     

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

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

The wave field patterns of the propagation of longitudinal and transverse elastic waves in grain-oriented rocks

In the work under consideration, on the basis of data analysis on the velocities of quasi-longitudinal Vp and quasi-transverse Vs waves—measured for rock samples at high pressures and temperatures, and analogous velocities, calculated on the basis of the grain orientation distribution function, reconstructed from the neutron-diffraction textural experiment—the indicative inconsistencies between the experimental and model characteristics were inferred. The theoretical analysis of the wave field patterns of the propagation of longitudinal and transverse elastic waves in the anisotropic media is carried out. It is established that, in the general case, in the anisotropic inhomogeneous media the velocities of Vp and Vs propagation, measured experimentally and obtained from the modeling, cannot coincide due to the existence of the physical coupling between the vibrations of two types: transverse and longitudinal vibrations.     

Converted PS-wave Reflection Coefficients in Weakly Anisotropic Media

—?I derive converted P S-wave reflection coefficients at a planar weak-contrast interface separating two weakly anisotropic half-spaces using first-order perturbation theory. The general expressions are further specified for the interface separating any of the two following media: isotropic, transversely isotropic with a vertical symmetry axis (VTI), transversely isotropic with a horizontal symmetry axis (HTI) and orthorhombic. Relatively simple forms of small-angle reflection coefficients are also obtained. The coefficients are expressed as functions of Thomsen-type medium parameters and incidence and azimuthal phase angles. Derived expressions, as well as their application, are more complicated than the corresponding expressions for P P-wave reflection coefficients. General characteristics and pitfalls are discussed. Numerical tests reveal a good agreement between exact and approximate coefficients for most models presented.     

Determination of the reference symmetry axis of a generally anisotropic medium which is approximately transversely isotropic

For a given stiffness tensor (tensor of elastic moduli) of a generally anisotropic medium, we estimate to what extent the medium is transversely isotropic (uniaxial) and determine the direction of its reference symmetry axis expressed in terms of the unit reference symmetry vector. If the medium is exactly transversely isotropic (exactly uniaxial), we obtain the direction of its symmetry axis. We can also calculate the first–order and second–order spatial derivatives of the reference symmetry vector which may be useful in tracing the reference rays for the coupling ray theory. The proposed method is tested using various transversely isotropic (uniaxial) and approximately transversely isotropic (approximately uniaxial) media.     

任意空间取向TI弹性张量解析表述

本文理论给出任意空间取向TI(ATI)四阶弹性张量的解析表述,其以VTI弹性常数及其简单组合为系数,包括各向同性项、TI对称轴方向矢量分量的二次项和四次项,其中TI对称轴方向矢量可以在固定坐标系定义, 也可以相对三维倾斜界面甚至相对波传播方向.相比四阶张量变换法和Bond变换法,ATI弹性张量能简洁而透明地为本构关系和波动方程提供四阶张量的所有元素. ATI弹性张量能为诸多方面的理论研究提供支撑.     

P-SH Conversions in Layered Media with Hexagonally Symmetric Anisotropy: A CookBook

—Reflectivity synthetic seismograms demonstrate that the type, layering and orientation of 1-D anisotropy influences strongly the coda of teleseismic P waves at periods T > 1 sec, particularly P-SH converted waves. We assume the simplest form of anisotropy described by an elastic tensor with a symmetry axis ? of arbitrary orientation. The resulting phase velocities vary as cos 2ξ with respect to that axis. Using three families of simple crustal models, we compare the effects of an anisotropic surface layer with reverberations caused by both "thick" and "thin" layers of anisotropy at depth. If anisotropy in the surface layer is significant, the polarization of direct P can be distorted to generate a transverse component, followed by Ps and a prominent shear reverberation converted from direct P at the free surface. If the anisotropic layer is buried, the first, and often the most prominent, arrival on the transverse component is the P-to-SH conversion at its upper surface. If the anisotropic layer is sufficiently thin, P-to-SH conversions from its boundaries interfere to form a derivative pulse shape on the transverse component, which could be mistaken as the signature of shear-wave splitting. If ? is horizontal, compressional (P) and shear (S) anisotropy both produce similar waveform perturbations with four-lobed azimuthal patterns, suggesting that a weighted stack of P coda from different back-azimuths would improve signal-to-noise. For ? tilted between the horizontal and vertical, however, the effects of P- and S-anisotropy differ greatly. The influence of P-anisotropy on P-to-S conversion is greatest for a symmetry axis tilted at 45° to the vertical, where its azimuthal pattern has two lobes, rather than four. Combinations of P- and S-anisotropy typically lead to a composite azimuthal dependence in the P-coda reverberations.     

Space-time variations in the parameters of split S waves excited by local earthquakes in southern Kamchatka

An analysis of the distribution of split S-wave parameters due to small local earthquakes occurring down to a depth of 200 km beneath the PET IRIS station (in the town of Petropavlovsk-Kamchatskii) for the period 1993–2002 revealed an inhomogeneous distribution of anisotropic properties and time-dependent changes in the state of stress and strain in the sinking slab. Three time intervals have been identified from the behavior of the split-wave parameters: 1993–1995, 1996–1998, and 1999–2002. The most orderly orientation of fast-wave polarization azimuths at all depths and lower stresses in the earth was observed in 1996–1998. Regions of relatively stable and unstable parameter behavior consistent with increased and decreased P and S velocities have been identified. The central block down to a depth of 120 km and the lower part of the Benioff zone are classified as stable regions dominated by E-SE fast-wave azimuths. Unstable behavior of fast-wave azimuths was noticed for depths of 60 to 90 km and around the top of the Benioff zone. The highest values of S-wave delay time and regions where the delay times were migrating mostly occur at contacts between more rigid and weakened blocks.     

Seismic Velocities and Anisotropy of the Lower Continental Crust: A Review

—Seismic anisotropy is often neglected in seismic studies of the earth’s crust. Since anisotropy is a common property of many typically deep crustal rocks, its potential contribution to solving questions of the deep crust is evaluated. The anisotropic seismic velocities obtained from laboratory measurements can be verified by computations based on the elastic constants and on numerical data pertaining to the texture of rock-forming minerals. For typical lower crustal rocks the influence of layering is significantly less important than the influence of rock texture. Surprisingly, most natural lower crustal rocks show a hexagonal type of anisotropy. Maximum anisotropy is observed for rocks with a high content of aligned mica. It seems possible to distinguish between layered intrusives and metasediments on the basis of in situ measurements of anisotropy, which can thus be used to validate different scenarios of crustal evolution.     

利用伪谱法模拟横向各向同性介质中的波

介质的弹性常数为三维四阶张量的分量,共有８１个,由于应力张量和应变张量的对称性及能量密度是应变的二次函数,一般各向异常性介质的独立弹性常数可减为２１个,如果介质具有较高的对称性,独立弹性常数的数目会更少。 对于地壳和上地幔,具有５个独立弹性常数的横向各向同性介质是一个非常好的近似,本研究中横向各向同性介质的对称轴方向可以是任意的（即对称轴可以不平等于铅直方向）,在此情况下,需要进行坐标变换,如果已知介质在某一坐标系（其坐标轴平行或垂直于介质的对称轴）中的弹性常数,我们能够容易地利用变换公式得到变换后新坐标系中的弹性常数。 本文提出了一种方案,利用伪谱法既能模拟横向各向同性介质中的平面波,也能模拟点源激发的波场。在勘探地球物理和地震学中,模拟横向各向同性介拮中传播的平面波及区域源产生的波是最重要的研究课题之一。然而在一般各向异性介质中,很难或不可能确定弹性波的相速度和偏振方向,但在横向各向同性介质中,则可以通过坐标变换来实现,这里我们所提出的方法可以用于横向各向同性介质中弹性波的模拟。     

Ambiguous Moment Tensors and Radiation Patterns in Anisotropic Media with Applications to the Modeling of Earthquake Mechanisms in W-Bohemia

Anisotropic material properties are usually neglected during inversions for source parameters of earthquakes. In general anisotropic media, however, moment tensors for pure-shear sources can exhibit significant non-double-couple components. Such effects may be erroneously interpreted as an indication for volumetric changes at the source. Here we investigate effects of anisotropy on seismic moment tensors and radiation patterns for pure-shear and tensile-type sources. Anisotropy can significantly influence the interpretation of the source mechanisms. For example, the orientation of the slip within the fault plane may affect the total seismic moment. Also, moment tensors due to pure-shear and tensile faulting can have similar characteristics depending on the orientation of the elastic tensor. Furthermore, the tensile nature of an earthquake can be obscured by near-source anisotropic properties. As an application, we consider effects of inhomogeneous anisotropic properties on the seismic moment tensor and the radiation patterns of a selected type of micro-earthquakes observed in W-Bohemia. The combined effects of near-source and along-path anisotropy cause characteristic amplitude distortions of the P, S1 and S2 waves. However, the modeling suggests that neither homogeneous nor inhomogeneous anisotropic properties alone can explain the observed large non-double-couple components.The results also indicate that a correct analysis of the source mechanism, in principle, is achievable by application of anisotropic moment tensor inversion.     

Texture and anisotropy analysis of Qusaiba shales

Scanning and transmission electron microscopy, synchrotron X‐ray diffraction, microtomography and ultrasonic velocity measurements were used to characterize microstructures and anisotropy of three deeply buried Qusaiba shales from the Rub’al‐Khali basin, Saudi Arabia. Kaolinite, illite‐smectite, illite‐mica and chlorite show strong preferred orientation with (001) pole figure maxima perpendicular to the bedding plane ranging from 2.4–6.8 multiples of a random distribution (m.r.d.). Quartz, feldspars and pyrite crystals have a random orientation distribution. Elastic properties of the polyphase aggregate are calculated by averaging the single crystal elastic properties over the orientation distribution, assuming a nonporous material. The average calculated bulk P‐wave velocities are 6.2 km/s (maximum) and 5.5 km/s (minimum), resulting in a P‐wave anisotropy of 12%. The calculated velocities are compared with those determined from ultrasonic velocity measurements on a similar sample. In the ultrasonic experiment, which measures the effects of the shale matrix as well as the effects of porosity, velocities are smaller (P‐wave maximum 5.3 km/s and minimum 4.1 km/s). The difference between calculated and measured velocities is attributed to the effects of anisotropic pore structure and to microfractures present in the sample, which have not been taken into account in the matrix averaging.