Elastic wave propagation in inhomogeneous anisotropic media

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｅｄｉａｉｎｅａｒｔｈａｒｅｑｕｉｔｅｃｏｍｐｌｅｘ．Ｔｈｅｒｅｅｘｉｓｔｓｅｖｅｒａｌｕｎｃｏｎｔｉｎｕｏｕｓｐｌａｉｎｓ．Ｎｏｒｍａｌｙ,ｉｔｉｓｕｓｅｄｔｏａｐｐｒｏｘｉｍａｔｅｒｅａｌｍｅｄｉｕｍｗｉｔｈｌａｙ...     

Love wave propagation in multilayered anisotropic inhomogeneous media

Summary Love wave propagation in a finite set of anisotropic inhomogeneous layers lying between two anisotropic homogeneous half spaces is considered. Generalized frequency equation is obtained by using the Thomson-Haskell matrix method. The usefulness of the general analytical result for discussing more special cases of interest in seismology is brought out in the end.     

Elastic wave propagation in an irregularly layered medium

The indirect boundary element method (IBEM) is used to simulate wave propagation in two-dimensional irregularly layered elastic media for internal line sources. The method is based on the integral representation for scattered elastic waves using single layer boundary sources. Fulfillment of the boundary conditions leads to a system of integral equations. Results are obtained in the frequency domain and seismograins are computed through Fourier synthesis. In order to test and validate the method we present various comparisons between our results and the time series obtained analytically for a buried line source in a half-space and by using the recently developed spectral element method (SEM).     

横向各向同性介质中地震波场谱元法数值模拟

横向各向同性介质是地球内部广泛存在的一种各向异性介质,因此为了能够更好地认识地震波在这种介质中的传播特征,用数值方法进行地震波模拟显得十分必要.本文采用谱元法对横向同性介质中的地震波进行模拟,该方法基于弹性力学方程弱形式基础之上,具有有限元适应任意复杂介质模型的韧性和伪谱法的精度.文中阐述了基于Legendre多项式的谱元法的理论和推导过程,该方法可以形成全局对角质量矩阵,在时间域使用显式的差分算法,提高运算效率,最后通过横向各向同性介质的数值计算,模拟结果表明该方法是一种有效的数值模拟方法.     

Elastic wave scattering by inhomogeneous and anisotropic bodies in a half space

Scattering of elastic waves by inhomogeneous and anisotropic bodies in a half space is considered. The integral equation method is formulated by using the fundamental solution of a homogeneous isotropic body in elastostatics and regarding the resulting inhomogeneous terms as equivalent body forces. Numerical examples are presented for the wave scattering by inhomogeneous and/or anisotropic alluvial valleys and for the dynamic analysis of an inhomogeneous dam. The effect of inhomogeneities and anisotropy on the dynamic behaviour of alluvial valley and dam is discussed.     

弱各向异性介质弹性波的准各向同性近似正演模拟

准各向同性(QI)近似可用于弱各向异性介质的正演模拟.本文通过运用QI方法的零阶和一阶近似,计算了VTI介质模型的地震记录.得出的地震记录与标准各向同性射线理论(IRT)和基于伪谱法的三维地震正演模拟得出的地震记录作了比较,可以认为是精确的合成地震记录.     

SH-wave propagation in anisotropic inhomogeneous crustal layer

Summary In this short contribution,SH-wave propagation in an anisotropic inhomogeneous crustal layer lying on an yielding and rigid isotropic half space is considered. The inhomogeneity is assumed to be present in the directional rigidities and the density. The type of variation in elastic parameters considered herein, is such that the velocities ofSH-waves in horizontal and vertical directions are constant. The frequency equations, governing propagation ofSH-waves are derived for both the cases, i.e. (I) when the lower medium is yielding half space and (II) when it is rigid. Dispersion curves for these cases are also presented.     

Love waves in inhomogeneous and anisotropic earth — I

Summary The frequency equation is derived for the propagation of Love waves in the earth's crust, composed of transversely isotropic layers and overlying anisotropic and inhomogeneous mantle. The exact boundary value problem is solved for a single layer and extended to multilayered media by generalizing theHaskell's technique. In fact the problem of deriving the frequency equation has been reduced to finding out the solution of the equation of motion subject to the appropriate boundary conditions. To illustrate the method, the author has derived frequency equations of Love waves for linear, exponential and generalized power law variation of vertical shear wave velocity with depth in the half space overlain by transversely isotropic inhomogeneous stratum.     

On the wave propagation in inhomogeneous non-isotropic media

Summary FollowingEason, we have discussed here the propagation of elastic waves in non-homogeneous spheres and cylinders when the curved surface is given a uniform normal loading. The material is assumed to be transversely isotropic with respect to a direction of symmetry, the stress and displacement components within the body may be assumed to depend on one space co-ordinate and time alone. The particular case in which the elastic parameters are proportional to (radius) ⁿ has been considered.     

Love wave dispersion: A comparison of results for a semi-infinite medium with inhomogeneous layers and for its approximation by homogeneous layers

Love wave dispersion in various semi-infinite media consisting of inhomogeneous layers is discussed. The phase and group velocities are computed when shear wave velocity and density in each inhomogeneous layer are varying exponentially with depth. At the beginning one or two inhomogeneous layers over a homogeneous semi-infinite medium are considered. The dispersion results for these structures are compared with those for their approximations with homogeneous layers. Comparisons show that differences of phase and group velocities for the original models from those for their approximated models (i) increase with the increase of wave number and (ii) are larger for group velocity than for phase velocity. The difference is approximately proportional to the rate of change of parameters in the layers. Finally, dispersion curves are obtained for model IP3MC, which consists of many inhomogeneous and homogeneous layers over a homogeneous semi-infinite medium. The results are compared with the observed group velocity data across the Indian Peninsula.     

Rayleigh wave propagation in a two-layer anisotropic media

Summary The possibility of propagation of Rayleigh waves in an incompressible crust of constant density and rigidity varying exponentially with depth lying on a semi-infinite transversely isotropic base has been discussed in this paper. Frequency equation has been derived and numerical calculations are made. The result obtained in this case is compared with that ofNewlands [3]²) andDutta [4].     

Exact solutions ofSH wave equation in transversely isotropic inhomogeneous elastic media

Summary TheSH wave equation in a transversely isotropic inhomogeneous elastic medium, where the elastic parameters and density are functions of vertical coordinate, is considered. A general procedure is given for finding the inhomogeneities for which the equation can be solved in terms of hypergeometric, Whittaker, Bessel and exponential functions. A few simple inhomogeneities and the corresponding solutions in terms of these transcendental functions are presented.     

热弹性介质中波传播特征

深层-超深层油气地震勘探涉及高温介质地震波传播问题,热弹介质参数对地震波传播有重要影响.含弛豫时间修正项的Lord-Shulman双曲型耦合热弹波动方程从理论上预测了热弹性介质中存在快纵波、慢纵波(一种准静态慢纵波,简称热波)和横波的传播,两个纵波为热耗散衰减波而横波不受介质热特性的影响.本文结合平面波频散分析和格林函...     

裂缝诱导各向异性双孔隙介质波场传播特征

基于裂缝诱导各向异性和双相介质理论,对裂缝诱导的具有水平对称轴的横向各向同性（HTI）双孔隙介质的本构关系进行了研究,与等效连续介质模型相结合,综合考虑裂缝系统和基质孔隙系统的两种孔隙度和两种渗透率参数,得到裂缝诱导HTI双孔隙介质的等效孔隙度和等效渗透率,进而得到介质的运动平衡方程;并进一步推导出介质的一阶速度-应力方程.采用交错网格高阶有限差分法对模型进行了数值模拟,结果揭示了介质中两套系统的存在对其波场传播特征的影响,为进一步研究实际地球介质的波场特征奠定了基础.     

An equivalent medium model for wave simulation in fractured porous rocks

Seismic wave propagation in reservoir rocks is often strongly affected by fractures and micropores. Elastic properties of fractured reservoirs are studied using a fractured porous rock model, in which fractures are considered to be embedded in a homogeneous porous background. The paper presents an equivalent media model for fractured porous rocks. Fractures are described in a stress‐strain relationship in terms of fracture‐induced anisotropy. The equations of poroelasticity are used to describe the background porous matrix and the contents of the fractures are inserted into a matrix. Based on the fractured equivalent‐medium theory and Biot's equations of poroelasticity, two sets of porosity are considered in a constitutive equation. The porous matrix permeability and fracture permeability are analysed by using the continuum media seepage theory in equations of motion. We then design a fractured porous equivalent medium and derive the modified effective constants for low‐frequency elastic constants due to the presence of fractures. The expressions of elastic constants are concise and are directly related to the properties of the main porous matrix, the inserted fractures and the pore fluid. The phase velocity and attenuation of the fractured porous equivalent media are investigated based on this model. Numerical simulations are performed. We show that the fractures and pores strongly influence wave propagation, induce anisotropy and cause poroelastic behaviour in the wavefields. We observe that the presence of fractures gives rise to changes in phase velocity and attenuation, especially for the slow P‐wave in the direction parallel to the fracture plane.     

Elastic constants computed for a medium composed of this anisotropic layers

Summary The problem of the elastic properties of a medium composed of thin anisotropic layers is treated. The study is based on the conditions of stress and strain and on Hooke's Law under the assumption of close contact between the layers. The algorithm described is suitable for a computer.Dedicated to RNDr. Jan Pícha, CSc., on his 60th Birthday     

Love waves in a sedimentary layer lying over an anisotropic inhomogeneous half space

Summary The effects of anisotropy and inhomogeneity on the propagation of Love waves in a sedimentary layer, overlying the inhomogeneous and transversely isotropic half space, are studied in this paper. The results of numerical analysis show an appreciable variation of phase- and group-velocity of Love waves in low frequency region compared to high frequency region due to the presence of transverse isotropy and inhomogeneity in the half space. The higher values for phase velocity are found for the increasing values of anisotropy factor as well as for the greater power of density variation. However, the presence of higher anisotropy factor and inhomogeneity in the half space reduce group velocity considerably in the lower frequency region.     

Elastic wave propagation and scattering in prestressed porous rocks

Poro-acoustoelastic theory has made a great progress in both theoretical and experimental aspects, but with no publications on the joint research from theoretical analyses, experimental measurements, and numerical validations. Several key issues challenge the joint research with comparisons of experimental and numerical results, such as digital imaging of heterogeneous poroelastic properties, estimation of acoustoelastic constants, numerical dispersion at high frequencies and strong heterogeneities, elastic nonlinearity due to compliant pores, and contamination by boundary reflections. Conventional poroacoustoelastic theory, valid for the linear elastic deformation of rock grains and stiff pores, is modified by incorporating a dualporosity model to account for elastic nonlinearity due to compliant pores subject to high-magnitude loading stresses. A modified finite-element method is employed to simulate the subtle effect of microstructures on wave propagation in prestressed digital cores. We measure the heterogeneity of samples by extracting the autocorrelation length of digital cores for a rough estimation of scattering intensity. We conductexperimental measurements with a fluid-saturated sandstone sample under a constant confining pressure of 65 MPa and increasing pore pressures from 5 to 60 MPa. Numerical simulations for ultrasound propagation in the prestressed fluid-saturated digital core of the sample are followed based on the proposed poro-acoustoelastic model with compliant pores. The results demonstrate a general agreement between experimental and numerical waveforms for different stresses, validating the performance of the presented modeling scheme. The excellent agreement between experimental and numerical coda quality factors demonstrates the applicability for the numerical investigation of the stress-associated scattering attenuation in prestressed porous rocks.     

A note on wave propagation in granular medium

Summary The body wave propagation in granular medium is discussed in two sections. In Section I, the solutions of general equations of granular medium involving displacement vector and rotation vector are shown to be dependent on the solutions of four different equations. In Section II, the plane wave propagation is studied, showing the dispersion and dissipation of body waves in the granular medium.