Wave propagation in layered media by time domain BEM

A time domain boundary element in a cylindrical co-ordinate system is developed for the analysis of wave propagation in a layered half-space. The field quantities (displacements and tractions) are expressed as products of Fourier series in the circumferential direction and as linear polynomials in the other spatial directions. An integral equation is written for each layer as an independent domain, and these equations are then assembled into a general equation by virtue of compatibility and equilibrium conditions between the interfaces. Examples of three-dimensional wave propagation in the layered half-spaces due to various forms of surface and inner-domain excitations are reported to demonstrate the accuracy and versatility of the method.     

Wave propagation in semi-infinite media with topographical irregularities using Decoupled Equations Method

In this paper, a novel semi-analytical method, called Decoupled Equations Method (DEM), is presented for modeling of elastic wave propagation in the semi-infinite two-dimensional (2D) media which are involved surface topography. In the DEM, only the boundaries of the problem are discretized by specific subparametric elements, in which special shape functions as well as higher-order Chebyshev mapping functions are implemented. For the shape functions, Kronecker Delta property is satisfied for displacement function. Moreover, the first derivatives of displacement function with respect to the tangential coordinates on the boundaries are assigned to zero at any given node. Employing the weighted residual method and using Clenshaw–Curtis numerical integration, coefficient matrices of the system of equations are transformed into diagonal ones, which leads to a set of decoupled partial differential equations. To evaluate the accuracy of the DEM in the solution of scattering problem of plane waves, cylindrical topographical features of arbitrary shapes are solved. The obtained results present excellent agreement with the analytical solutions and the results from other numerical methods.     

Wave propagation in a seven-story reinforced concrete building: I. Theoretical models

For transient, high frequency, and pulse like excitation of structures in the near field of strong earthquakes, the classical design approach based on relative response spectrum and mode superposition may not be conservative. For such excitations, it is more natural to use wave propagation methods. In this paper (Part I), we review several two-dimensional wave propagation models of buildings and show results for theoretical dispersion curves computed for these models. We also estimate the parameters of these models that would correspond to a seven-story reinforced concrete building in Van Nuys, California. Ambient vibration tests data for this building imply vertical shear wave velocity β_z=112 m/s and anisotropy factor β_x/β_z=0.55 for NS vibrations, and β_z=88 m/s and β_x/β_z=1 for EW vibrations. The velocity of shear waves propagating through the slabs is estimated to be about 2000 m/s. In the companion paper (Part II), we estimate phase velocities of vertically and horizontally propagating waves between seven pairs of recording points in the building using recorded response to four earthquakes.     

Wave propagation in a uniformly rotating fluid

Summary An unsteady flow generated by a harmonically oscillating pressure distribution of frequency acting on the paraboloidal free surface of an inviscid, incompressible fluid rotating with uniform angular velocity has been investigated. It is shown that case (i), >2 , corresponds to the usual surface waves, and case (ii), <2 , in contrast to the surface waves, corresponds to the inertial waves which are originated entirely due to rotation and have no counterpart in a non-rotating fluid motion. An explicit solution of the problem related to the above cases are obtained by the joint Laplace and Hankel transforms treatment in conjunction with asymptotic methods. The significant effects of the Coriolis force and the curvature of the free surface on the wave motions have been investigated. A comparison is made between the solutions of the problems with the horizontal and the paraboloidal free surface curvature. The analysis is concluded by exihibiting the characteristic features of the wave motions.     

Elastic wave propagation in inhomogeneous anisotropic media

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

Wave propagation in a rotating fluid of spherical configuration

Abstract

An asymptotic approximation to the solution of the time-dependent linearized equations governing the motion of an incompressible, inviscid rotating fluid of spherical configuration having uniform density, variable depth and a free upper surface is obtained using the ray method without a shallow water assumption. This result is then modified to obtain a ray approximation to the solution of the time-reduced problem and the free oscillations of the fluid are studied. Axisymmetric modes covering the whole sphere and asymmetric modes trapped in both equatorial and non-equatorial regions are discovered, and all these modes are shown to have countably many resonance frequencies. A shallow water limit is defined and this limit of the time-reduced approximation is obtained. Most of the modes of free oscillation are lost in this limit and the limiting axisymmetric modes are shown to be trapped in the equatorial region and are singular at the wave region boundaries. The limiting approximation is compared to previous results obtained under a shallow water assumption.     

Earthquake processes in a micromorphic continuum

Summary The paper introduces a new model earthquake process based on the theory of micromorphic continua. The processes in a focal region are described by deformations of microstructure in time. It is assumed that the fracturing processes as well as phase transformation of metamorphic phenomena have caused in the past certain non-reversible changes which determine the microstructure of focal region. These internal microstructural elements form the attaching points around which the couple stresses arise. The properties of focal region are determined by the constitutive equations. The micromorphic mechanics considers the existence of body couples as determined by a regional stresses and looks after a response field of stresses, stress moments and strains in the focal region. Further, it is explained how microdislocation field is connected with microdeformations and micromorphic structure. In the considered earthquake structure model a microanisotropy is assumed through the tensor of microinertia. This tensor describes a distribution of microelements. Simple solutions of wave processes in a focal region are presented. The dispersion of waves is discussed.     

Kelvin粘弹性VTI介质中地震波的传播

In this article, under the assumption of weak anisotropy and weak attenuation, we present approximate solutions of anisotropic complex velocities and quality-factors for Kelvin visco-elastic transverse isotropy (KEL-VTI) media, based on the complex physical parameter matrix. Also, combined with the KEL-VTI media model, the propagation characteristics of the qP-, qSV-, and qSH-wave phases and energies are discussed. Further, we build a typical KEL-VTI media model of the Huainan coal mine to model the wave propagation. The numerical simulation results show that the PP- and PSV-wave theoretical wave-fields are close to the wave-fields of three-component P- and converted-waves acquired in the work area. This result proves that the KEL-VTI media model gives a good approximation to this typical coalfield seismic-geologic conditions and is helpful to the study of attenuation compensation of multi-component seismic data.     

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.     

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].     

近地表强衰减介质中的地震波传播模型

建立符合油气储层近地表复杂介质强衰减性质的数学-物理模型是油气勘探和开发的重要课题.本文针对D'Alembert黏弹性介质模型存在的量纲不统一和不能充分刻画孔隙介质结构的不足,通过引入孔隙度和渗透率等参数修改耗散项,改进了原D'Alembert模型,获得了能精细刻画具有强吸收衰减特征的近地表复杂介质模型,即改进的近地表黏弹性模型.基于这种改进的模型,推导了波频散和衰减的表达式,并研究了孔隙度、流体黏度等物理参数对波频散和衰减的影响,获得了相应的规律性认识.为了验证新模型预测近地表介质中波衰减的有效性,本文将新模型应用于胜利油田YX工区的近地表实测数据,同时与弹性Biot模型和BISQ模型以及黏弹性BISQ模型进行了比较.结果表明,与其它三个模型比较,改进的近地表黏弹性模型能够很好地刻画近地表介质的强衰减性,而且新模型所涉及到的物理参数明显少于其它几种模型,有利于新模型在油气储层近地表复杂介质地震勘探的实际应用.

     

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.     

PML solution of longitudinal wave propagation in heterogeneous media

This paper describes the development of a model for unbounded heterogeneous domains with radiation damping produced by an unphysical wave absorbing layer. The Perfectly Matched Layer (PML) approach is used along with a displacement-based finite element. The heterogeneous model is validated using the closed-form solution of a benchmark problem: a free rod with two-part modulus subjected to a specified time history. Both elastically supported and unsupported semi-infinite rods with different degrees of inhomogeneity and loading are considered. Numerical results illustrate the effects of inhomogeneity on the response and are compared with those for equivalent homogeneous domains. The effects of characteristic features of the inhomogeneous problem, presence of local maxima and cut-off frequency are determined. A degenerate case of a homogeneous semi-infinite rod on elastic foundations is produced by tending the magnitude of the foundation stiffness to zero. The response of the latter is compared with that of a free rod. The importance of proper selection of the PML parameters to highly accurate and efficient results is demonstrated by example problems.     

Wave propagation in marine sediments and water saturated soils

Summary A unification of the theories of Biot and Weiskopf has been made to form the suitable equations of motion for porous water saturated soils and marine sediments. It has been shown that the velocities of the body waves depend on the direction of propagation. In the vertical direction there are three, one distortional and two dilatational waves. In the horizontal direction there are two dilatational and two distortional waves. Finally, propagation of Love waves and Rayleigh waves have been discussed. Suitable potential functions have been derived to find the frequency equation for Rayleigh waves.     

Wave propagation in two layered medium under initial stresses

Summary In this paper, the frequency equation for phase velocity of waves propagated in a laminated medium consisting of two eleastic layers of finite thickness under initial stresses, has been obtained. It has been shown that when wave length becomes very small compared to the thickness of each layer, the wave approaches two Rayleigh waves at the two outer surfaces with the possibility of Stoneley waves at the interface. The propagation ofSH-waves in the composite medium under initial stresses has also been discussed. A particular case has been taken to find the velocity of Love wave in the homogeneous half space under initial compressive stresses.Biot's incremental deformation theory has been used.     

Wave propagation through a medium containing anisotropic permeability

Summary Reflection and refraction for an incident plane polarised electromagnetic wave have been studied for a gyrotropic medium of tensorial permeability as a boundary value problem. The evaluated field variables reveal the electromagnetic phenomena for such a medium.     

Finite difference modeling of seismic wave propagation in monoclinic media

Reported in the present paper are the results of the study of propagation of SH waves in the plane of mirror symmetry of a monoclinic multilayered medium with displacement normal to the plane. Dispersion equation has been obtained analytically ussing Haskell’s matrix method, while the finite-difference method has been employed to model the SH-wave propagation to study its phase and group velocities. The stability analysis has been carried out to minimize the exponential growth of the error of finite difference approximation in order to make the finite difference method stable and convergent. Further, variations of phase velocity with respect to both wave number and dispersion parameter for different stability ratios in monoclinic media have been examined and shown graphically. The effect of change of stability ratio on the group velocity of the wave propagation has been also investigated. Likewise, the effects of change of dispersion parameter on phase velocity and the variation of frequency with increase of wave number have been graphically represented and discussed.     

EDA介质中地震波的传播特征

对Christoffel公式进行Bond变换得到EDA介质的Christoffel方程,并由其非零解推导出EDA介质中视横波(qSV)、横波(SH)、视纵波(qP)的相速度、群速度、偏振向量(质点的振动方向)的三维计算公式.通过模型计算分析了具有水平对称轴的各向异性(HTI)介质和EDA介质中介质对称轴的极角和方位角对相速度、群速度及偏振向量的影响,对其随极角、方位角的变化特征进行了分析,并采用Matlab进行了数值计算,对其特征采用三维显示.通过取极角或方位角为零简化得到HTI介质和具有垂直对称轴的各向异性(VTI)介质中地震波的相速度、群速度,对EDA介质中的三维计算结果进行退化验证. 通过数值计算进一步验证了地震波相速度与EDA介质对称轴的相互关系. 结果表明,通过广角地震勘探可探明地下介质的裂隙走向及密度,从而确定灾害体产状.