Time-harmonic response of transversely isotropic multilayered half-space in a cylindrical coordinate system

With the aid of the analytical layer-element method, a comprehensive analytical derivation of the response of transversely isotropic multilayered half-space subjected to time-harmonic excitations is presented in a cylindrical coordinate system. Starting with the governing equations of motion and the constitutive equations of transversely isotropic elastic body, and based on the Fourier expansion, Hankel and Laplace integral transform, analytical layer-elements for a finite layer and a half-space are derived. Considering the continuity conditions on adjacent layers׳ interfaces and the boundary conditions, the global stiffness matrix equations for multilayered half-space are assembled and solved. Finally, some numerical examples are given to make a comparison with the existing solution and to demonstrate the influence of parameters on the dynamic response of the medium.     

Scattering of shear waves by an elliptical cavity in a radially inhomogeneous isotropic medium

Complex function and general conformal mapping methods are used to investigate the scattering of elastic shear waves by an elliptical cylindrical cavity in a radially inhomogeneous medium. The conformal mappings are introduced to solve scattering by an arbitrary cavity for the Helmholtz equation with variable coefficient through the transformed standard Helmholtz equation with a circular cavity. The medium density depends on the distance from the origin with a power-law variation and the shear elastic modulus is constant. The complex-value displacements and stresses of the inhomogeneous medium are explicitly obtained and the distributions of the dynamic stress for the case of an elliptical cavity are discussed. The accuracy of the present approach is verified by comparing the present solution results with the available published data. Numerical results demonstrate that the wave number, inhomogeneous parameters and different values of aspect ratio have significant influence on the dynamic stress concentration factors around the elliptical cavity.     

柱面波在半空间中洞室周围的散射

采用一种间接边界积分方程法求解了柱面波在半空间中洞室周围的散射问题。通过边界条件的验算以及退化解答与精确解的比较,验证了本文方法的计算精度。文中主要分析了入射波频率、波源与洞室距离等参数对洞室附近地表位移响应和洞室动应力集中的影响,得出了一些有益的结论。     

Pure acoustic wave propagation in transversely isotropic media by the pseudospectral method

Seismic wave propagation in transversely isotropic (TI) media is commonly described by a set of coupled partial differential equations, derived from the acoustic approximation. These equations produce pure P‐wave responses in elliptically anisotropic media but generate undesired shear‐wave components for more general TI anisotropy. Furthermore, these equations suffer from instabilities when the anisotropy parameter ε is less than δ. One solution to both problems is to use pure acoustic anisotropic wave equations, which can produce pure P‐waves without any shear‐wave contaminations in both elliptical and anelliptical TI media. In this paper, we propose a new pure acoustic transversely isotropic wave equation, which can be conveniently solved using the pseudospectral method. Like most other pure acoustic anisotropic wave equations, our equation involves complicated pseudo‐differential operators in space which are difficult to handle using the finite difference method. The advantage of our equation is that all of its model parameters are separable from the spatial differential and pseudo‐differential operators; therefore, the pseudospectral method can be directly applied. We use phase velocity analysis to show that our equation, expressed in a summation form, can be properly truncated to achieve the desired accuracy according to anisotropy strength. This flexibility allows us to save computational time by choosing the right number of summation terms for a given model. We use numerical examples to demonstrate that this new pure acoustic wave equation can produce highly accurate results, completely free from shear‐wave artefacts. This equation can be straightforwardly generalized to tilted TI media.     

横观各向同性饱和地基中埋置荷载的非轴对称瞬态响应

基于孔隙介质的Biot理论,首先利用Laplace变换,给出圆柱坐标系下横观各向同性饱和弹性多孔介质在变换域上的波动方程;将波动方程解耦后,根据方位角的Fourier展开和径向Hankel变换,求解了Biot波动方程,得到以土骨架位移、孔隙水压力和土介质总应力分量的积分形式的一般解;借助一般解,建立了有限厚度饱和土层和饱和半空间的精确动力刚度矩阵,并由土层的层间界面连续条件建立三维非轴对称层状饱和地基的总刚度方程;在此基础上,系统研究了横观各向同性饱和半空间体在内部集中荷载激励下的动力响应,并给出了问题的瞬态解答.该研究为运用边界元法求解饱和地基动力响应奠定了理论基础.     

The offset‐midpoint travel‐time pyramid for P‐wave in 2D transversely isotropic media with a tilted symmetry axis

For pre‐stack phase‐shift migration in homogeneous isotropic media, the offset‐midpoint travel time is represented by the double‐square‐root equation. The travel time as a function of offset and midpoint resembles the shape of Cheops’ pyramid. This is also valid for transversely isotropic media with a vertical symmetry axis. In this study, we extend the offset‐midpoint travel‐time pyramid to the case of 2D transversely isotropic media with a tilted symmetry axis. The P‐wave analytical travel‐time pyramid is derived under the assumption of weak anelliptical property of the tilted transverse isotropy media. The travel‐time equation for the dip‐constrained transversely isotropic model is obtained from the depth‐domain travel‐time pyramid. The potential applications of the derived offset‐midpoint travel‐time equation include pre‐stack Kirchhoff migration, anisotropic parameter estimation, and travel‐time calculation in transversely isotropic media with a tilted symmetry axis.     

Vertical and horizontal vibrations of a rigid disc on a multilayered transversely isotropic half-space

A half-space containing horizontally multilayered regions of different transversely isotropic elastic materials as well as a homogeneous half-space as the lowest layer is considered such that the axes of material symmetries of different layers and the lowest half-space to be as depth-wise. A rigid circular disc rested on the free surface of the whole half-space is considered to be under a forced either vertical or horizontal vibration of constant amplitudes. Because of the involved integral transforms, the mixed boundary value problems due to mixed condition at the surface of the half-space are changed to some dual integral equations, which are reduced to Fredholm integral equations of second kind. With the help of contour integration, the governing Fredholm integral equations are numerically solved. Some numerical evaluations are given for different combinations of transversely isotropic layers to show the effect of degree of anisotropy of different layers on the response of the inhomogeneous half-space.     

A tilted transversely isotropic slowness surface approximation

The relation between vertical and horizontal slownesses, better known as the dispersion relation, for transversely isotropic media with a tilted symmetry axis (TTI) requires solving a quartic polynomial equation, which does not admit a practical explicit solution to be used, for example, in downward continuation. Using a combination of the perturbation theory with respect to the anelliptic parameter and Shanks transform to improve the accuracy of the expansion, we develop an explicit formula for the vertical slowness that is highly accurate for all practical purposes. It also reveals some insights into the anisotropy parameter dependency of the dispersion relation including the low impact that the anelliptic parameter has on the vertical placement of reflectors for a small tilt in the symmetry angle.     

Forced vertical vibration of rigid circular disc buried in an arbitrary depth of a transversely isotropic half space

This paper is concerned with the investigation of the vertical vibration of a rigid circular disc buried at an arbitrary depth in a transversely isotropic half space in such a way the axis of material symmetry of the half space is normal to the surface of it and parallel to the vibration direction. By using the Hankel integral transforms, the mixed boundary-value problem is transformed to a pair of integral equations called dual integral equations, which generally can be reduced to a Fredholm integral equation of the second kind. With the aid of complex variable or contour integration, the governing integral equation is numerically solved in the general dynamic case. Two degenerated cases (i) the disc is buried in a transversely isotropic full space, and (ii) rigid circular disc is attached on the surface of the half space are discussed. The reduced static case of the dual integral equations is solved analytically and the vertical displacement, the contact pressure and the static impedance/compliance function are explicitly found. It is shown that the vertical pressure and the compliance function reduced for isotropic half space are identical to the previous solutions reported in the literature. The dynamic contact pressure under the disc and the impedance function are numerically evaluated in general dynamic case and graphically shown that the singularity exists in the contact pressure at the edge of the disc is the same as the static case. In addition, the impedance functions evaluated here for the isotropic domain are collapsed on the solution given by Luco and Mita. To show the effect of different material anisotropy, the numerical evaluations are given for some different transversely isotropic materials and compared.     

横观各向同性层状半空间中的弹性位错

本文在柱向量函数系下,利用传播矩阵法求解了层状横观各向同性半空间由内部点源位锚引起的变形;对六个基本点源位错,以等价体力法推出了横观各向同性情形下的点源函数,并且给出了内部任意剪切位错源引起的地表位移的积分表达式。为研究地球的层状结构,特别是其上部的横观各向同性对地表的地震位移、应变以及倾斜场的影响提供了计算公式。     

横观各向同性层状半空间中的弹性位错

本文在柱向量函数系下,利用传播矩阵法求解了层状横观各向同性半空间由内部点源位锚引起的变形;对六个基本点源位错,以等价体力法推出了横观各向同性情形下的点源函数,并且给出了内部任意剪切位错源引起的地表位移的积分表达式。为研究地球的层状结构,特别是其上部的横观各向同性对地表的地震位移、应变以及倾斜场的影响提供了计算公式。     

横观各向同性岩土材料的应变局部化分析

为分析横观各向同性岩土材料的应变局部化行为,推导对应于Lade横观各向同性屈服准则向前欧拉算法的迭代格式及本构矩阵,通过有限元软件ABAQUS的UMAT接口实现相应的程序代码。数值算例调查分析材料主方向对横观各向同性岩土结构的极限承载力和变形局部化模式的影响。结果表明该模型能较好地模拟横观各向同性岩土结构的应变局部化行为。     

Scattering pattern analysis and true-amplitude generalized Radon transform migration for acoustic transversely isotropic media with a vertical axis of symmetry

In multi-parameter ray-based anisotropic migration/inversion, it is essential that we have an understanding of the scattering mechanism corresponding to parameter perturbations. Because the complex nonlinearity in the anisotropic inversion problem is intractable, the construction of true-amplitude linearized migration/inversion procedures is needed and important. By using the acoustic medium assumption for transversely isotropic media with a vertical axis of symmetry and representing the anisotropy with P-wave normal moveout velocity, Thomsen parameter δ and anelliptic parameter η, we formalize the linearized inverse scattering problem for three-dimensional pseudo-acoustic equations. Deploying the single-scattering approximation and an elliptically anisotropic background introduces a new linear integral operator that connects the discontinuous perturbation parameters with the multi-shot/multi-offset P-wave scattered data. We further apply the high-frequency asymptotic Green's function and its derivatives to the integral operator, and then the scattering pattern of each perturbation parameter can be explicitly presented. By naturally establishing a connection to generalized Radon transform, the pseudo-inverse of the integral operator can be solved by the generalized Radon transform inversion. In consideration of the structure of this pseudo-inverse operator, the computational implementation is done pointwise by shooting a fan of rays from the target imaging area towards the acquisition system. Results from two-dimensional numerical tests show amplitude-preserving images with high quality.     

Scanning anisotropy parameters in horizontal transversely isotropic media

The horizontal transversely isotropic model, with arbitrary symmetry axis orientation, is the simplest effective representative that explains the azimuthal behaviour of seismic data. Estimating the anisotropy parameters of this model is important in reservoir characterisation, specifically in terms of fracture delineation. We propose a travel‐time‐based approach to estimate the anellipticity parameter η and the symmetry axis azimuth ? of a horizontal transversely isotropic medium, given an inhomogeneous elliptic background model (which might be obtained from velocity analysis and well velocities). This is accomplished through a Taylor's series expansion of the travel‐time solution (of the eikonal equation) as a function of parameter η and azimuth angle ?. The accuracy of the travel time expansion is enhanced by the use of Shanks transform. This results in an accurate approximation of the solution of the non‐linear eikonal equation and provides a mechanism to scan simultaneously for the best fitting effective parameters η and ?, without the need for repetitive modelling of travel times. The analysis of the travel time sensitivity to parameters η and ? reveals that travel times are more sensitive to η than to the symmetry axis azimuth ?. Thus, η is better constrained from travel times than the azimuth. Moreover, the two‐parameter scan in the homogeneous case shows that errors in the background model affect the estimation of η and ? differently. While a gradual increase in errors in the background model leads to increasing errors in η, inaccuracies in ?, on the other hand, depend on the background model errors. We also propose a layer‐stripping method valid for a stack of arbitrary oriented symmetry axis horizontal transversely isotropic layers to convert the effective parameters to the interval layer values.     

Three-dimensional Green's functions for transversely isotropic thermoporoelastic bimaterials

Green's functions for transversely isotropic thermoporoelastic bimaterials are established in the paper. We first express the compact general solutions of transversely isotropic thermoporoelastic material in terms of harmonic functions and introduce eight new harmonic functions. The three-dimensional Green's function having a concentrated liquid source or a concentrated heat source in steady state is completely solved using these new harmonic functions. The analytical results show some new phenomena of the pore fluid pressure increment, thermal increment and stress distributions at the interface. In the two materials, the pore fluid pressure has the same distribution because of the common fluid permeability, but the situation is different for the thermal increment. Shear failure is most likely at the two sources due to the highly degenerated direction of shear stress contours.     

On the stresses produced by a moving stress discontinuity at the boundary between two transversely isotropic media

Summary This paper deals with the stresses produced by a moving stress discontinuity between two semi-infinite transversely isotropic media. The solutions for the stresses have been obtained in simple closed forms for certain types of moving stress discontinuity. Discontinuities in the normal stress as well as in the shear stress has been considered. The case in which one or both of the media are isotropic may be deduced as a special case from the results obtained here.     

An efficient wave extrapolation method for anisotropic media with tilt

Wavefield extrapolation operators for elliptically anisotropic media offer significant cost reduction compared with that for the transversely isotropic case, particularly when the axis of symmetry exhibits tilt (from the vertical). However, elliptical anisotropy does not provide accurate wavefield representation or imaging for transversely isotropic media. Therefore, we propose effective elliptically anisotropic models that correctly capture the kinematic behaviour of wavefields for transversely isotropic media. Specifically, we compute source‐dependent effective velocities for the elliptic medium using kinematic high‐frequency representation of the transversely isotropic wavefield. The effective model allows us to use cheaper elliptic wave extrapolation operators. Despite the fact that the effective models are obtained by matching kinematics using high‐frequency asymptotic, the resulting wavefield contains most of the critical wavefield components, including frequency dependency and caustics, if present, with reasonable accuracy. The methodology developed here offers a much better cost versus accuracy trade‐off for wavefield computations in transversely isotropic media, particularly for media of low to moderate complexity. In addition, the wavefield solution is free from shear‐wave artefacts as opposed to the conventional finite‐difference‐based transversely isotropic wave extrapolation scheme. We demonstrate these assertions through numerical tests on synthetic tilted transversely isotropic models.     

Dynamic stress concentration in pre-stressed poroelastic media due to moving punch influenced by shear wave

During the occurrence of earthquake, the shear wave propagates in the rocks present inside/at the Earth’s crust. The propagation of shear wave may lead to the progression of punch present inside the rock medium. As a result of this, substantial stress accumulated at the vicinity of propagating punch inside rock medium which significantly affects the stability of various geological and human-made structure and, hence, may cause failure of structure. Therefore, the analysis of stress concentration at the vicinity of punch moving due to shear wave propagation has become prominent in the area of seismology. In the present paper, an analytical perspective has been employed to discuss the influence of velocity of moving punch associated with the propagation of shear wave on developed dynamic stress concentration (DSC) in three types of pre-stressed vertical transversely isotropic (VTI) poroelastic media viz. granite (an igneous rock); sandstone (a sedimentary rock); and marble (a metamorphic rock). The closed form expression of DSC for the force of constant intensity has been derived with the aid of Weiner-Hopf technique along with Galilean and two-sided Fourier integral transformations. The noticeable influence of different affecting parameters (viz. velocity of moving punch associated with the shear wave propagation, horizontal compressive/tensile initial stresses, vertical compressive/tensile initial stress, porosity, and anisotropy parameter) on dynamic stress concentration has also been reported. Numerical computation and graphical illustrations have been carried out for the aforementioned three different types of porous rocks to investigate the profound impact of affecting parameters on DSC. Moreover, some noteworthy peculiarities have also been derived from the obtained expression of dynamic stress concentration.     

曲面上的位场转换（二维和三维情况）

在文献[1]中我们利用复变函数论的方法只解决了二维的曲面上的位场转换问题.在本文中,我们利用单层位势理论将它化为一个奇异积分方程求解的问题,并且提出了用奇点分离法来解此奇异积分方程,这对二维和三维情况都适用.通过理论模型的计算,结果的精度是令人满意的.