A New Theory of Love Waves in Multi-layered Media with Irregular Interfaces

In this article, we have derived a new and more general formulation of Love waves in arbitrarily irregular multi-layered media by using the global generalized reflection/transmission (abbreviated to R/T thereafter) matrices method developed earlier by Chen [17~20]. From the basic principle that the modal solutions are the non-trivial solutions of the free elastodynamic equation under appropriate boundary conditions, we naturally derived the characteristic frequencies and the corresponding distorted modes of Love wave in irregular multi-layered media. Moreover, we have derived the corresponding excitation formulation of Love waves in such laterally heterogeneous media by using the general solution of elastodynamic equation [17~20]. Similar to the result for laterally homogeneous layered structure, the Love waves radiated from a point source in irregular multi-layered media can be expressed as a superposition of distorted modes. Since the structure model used here is quite arbitrary, it can be used for so     

Seismic wavefield modeling in media with fluid-filled fractures and surface topography

We present a finite difference (FD) method for the simulation of seismic wave fields in fractured medium with an irregular (non-flat) free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations (in a curvilinear coordinate system) in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry (HTI). Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures.     

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second-order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite-element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.     

Vector-scattering of elastic waves by directional structural space gradients

A nonstochastic and noniterative theory of vector scattering in inhomogeneous media is presented. The elastodynamic vector wave-equation for 3D inhomogeneous media is solved for a weak heterogeneity at the high-frequency region. It is shown that there exists a forward scattered field which decays slowly along the source-receiver path. Its rate of attenuation depends on the azimuth of the path relative to the direction of the inhomogeneity, but is independent of frequency. The Green's tensor for the above regime is derived in closed form and leads to the quantification of fields of dipolar sources in weak inhomogeneous media. The inhomogeneity at the source creates a source-induced scattering (in addition to path-scattering) having a radiation-pattern that bears the signature of the source. The availability of the analytic Green's tensor, in conjunction with the Huygens-Kirchhoff-Helmholtz formalism, opens new ways to calculate the scattered fields due to various structural inhomogeneities applicable to exploration and earthquake seismology. The theoretical results of this study point to the conclusion that the scalar wave approximation may not always be valid for the propagation of seismic waves in the earth's lithosphere.     

A fundamental solution of a multilayered half-space due to an impulsive ring source

The fundamental solutions of axisymmetric elastodynamic problem for the multilayered half-space due to an impulsive ring source acting within a layered elastic media are derived in time domain with the aid of Laplace–Hankel mixed transform and transfer matrix techniques. In addition, an effective numerical procedure, which utilizes the fast Hankel transform algorithm, is also proposed to calculate these solutions. Illustrative examples have been given to demonstrate that the fundamental solutions can be readily evaluated and the numerical results are of high accuracy. The present solutions can be directly applied to determine the transient wave fields caused by a seismic source and show the potential application to the elastodynamic problems solved by the boundary element method.     

非均匀介质全局广义R/T递推传播矩阵法二维SH地震波模拟研究

地球深部圈层及沉积盆地是一种分区非均匀介质系统,其中不规则地层边界(含起伏地表)对地震波的主要特征有显著影响,而地层的随机非均匀性则主要影响地震波的散射和衰减特征.为了精确刻画不规则地层边界对地震波的反射、透射效应以及非均质体散射引起的地震波衰减效应,全局广义R/T递推传播矩阵法(GGRTM)被提出并逐步发展成为继有限...     

Least square approach to simulate wave propagation in irregular profiles using the indirect boundary element method

In earthquake engineering and seismology it is of interest to know the surface motion at a given site due to the incoming and scattered seismic waves by surface geology. This can be formulated in terms of diffraction of elastic waves and then the indirect boundary element method (IBEM) for dynamic elasticity is used. It is based on the explicit construction of diffracted waves at the boundaries from which they radiate. This provides the analyst with insight on the physics of diffraction. The IBEM has been applied to study the amplification of elastic waves in irregular soil profiles. From the strong or weak satisfaction of boundary conditions and a simple analytical discretization scheme a linear system of equations for the boundary sources is obtained. Here, we explore the use of a weak discretization strategy with more collocation points than force densities. The least squares enforcement of boundary conditions leads to a system with reduced number of unknowns. This approach naturally allows one to use both coarser and finer boundary discretizations for smooth and rapidly varying profiles, respectively. A well studied semicircular canyon under incident P or SV in-plane waves is used to calibrate this method. Several benefits are obtained using mixed meshing that leads to the least squares condensation of the IBEM.     

离散介质中地震波传播的数值模拟

伪谱法应用付里叶变换和有限差分求解波动方程，本利用伪谱法模拟了三维离散介质中地震波的传播。结果表明：(1)该方法的计算精度高；(2)即使在均匀、各向同性的离散介质中，地震波的传播也表现出方位各向异性的特点。     

切比雪夫伪谱法模拟地震波场

介绍了切比雪夫伪谱法以及快速傅立叶算法在其中的应用,并用切比雪夫伪谱法模拟二维有限区域弹性介质地震波场。分别计算了兰姆问题,均匀介质中心爆破源问题,介质内部的速度异常体问题以及各种分界面情况下地震首波传播问题。     

2d and 3D elastic wave propagation by a pseudo-spectral domain decomposition method

A new numerical method is presented for propagating elastic waves in heterogeneous earth media, based on spectral approximations of the wavefield combined with domain decomposition techniques. The flexibility of finite element techniques in dealing with irregular geologic structures is preserved, together with the high accuracy of spectral methods. High computational efficiency can be achieved especially in 3D calculations, where the commonly used finite-difference approaches are limited both in the frequency range and in handling strongly irregular geometries. The treatment of the seismic source, introduced via a moment tensor distribution, is thoroughly discussed together with the aspects associated with its numerical implementation. The numerical results of the present method are successfully compared with analytical and numerical solutions, both in 2D and 3D.     

Calculating complex‐valued P‐wave first‐arrival traveltimes in attenuative vertical transversely isotropic media with an irregular surface

The complex‐valued first‐arrival traveltime can be used to describe the properties of both velocity and attenuation as seismic waves propagate in attenuative elastic media. The real part of the complex‐valued traveltime corresponds to phase arrival and the imaginary part is associated with the amplitude decay due to energy absorption. The eikonal equation for attenuative vertical transversely isotropic media discretized with rectangular grids has been proven effective and precise to calculate the complex‐valued traveltime, but less accurate and efficient for irregular models. By using the perturbation method, the complex‐valued eikonal equation can be decomposed into two real‐valued equations, namely the zeroth‐ and first‐order traveltime governing equations. Here, we first present the topography‐dependent zeroth‐ and first‐order governing equations for attenuative VTI media, which are obtained by using the coordinate transformation from the Cartesian coordinates to the curvilinear coordinates. Then, we apply the Lax–Friedrichs sweeping method for solving the topography‐dependent traveltime governing equations in order to approximate the viscosity solutions, namely the real and imaginary parts of the complex‐valued traveltime. Several numerical tests demonstrate that the proposed scheme is efficient and accurate in calculating the complex‐valued P‐wave first‐arrival traveltime in attenuative VTI media with an irregular surface.     

基于精细积分法的三维弹性波数值模拟(英文)

波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。     

Consistency of the spatial autocorrelation method with seismic interferometry and its consequence

We have cross‐checked the conventional theory of the spatial autocorrelation method and the consequence of seismic interferometry: the retrieval of the elastodynamic Green's function. Their mutual consistency is almost complete. The basic formulas of the conventional spatial autocorrelation theory can be derived by an alternative approach based on the retrieval of the elastodynamic Green's function. The only discrepancy is found with the average of the complex coherence function over azimuth in a wavefield dependent on azimuth. It is hypothesized, in discussion, that this discrepancy is due to the way of representing the wavefield in the background theory of seismic interferometry that can produce only wavefields moderately dependent on azimuth and that the mentioned consequence of seismic interferometry can also only make sense in a wavefield moderately dependent on azimuth. Our field experiment with a wavefield dependent on azimuth showed that the consequence of seismic interferometry in the logical framework of the conventional spatial autocorrelation theory is appropriate under such degrees of approximation as the measure proposed in this study, i.e., the deviation of the total dispersion curves is between about 10 and 16 per cent at the maximum from those averaged over azimuth. The acceptance of the retrieval of Green's function gives a proper physical meaning to the complex coherence function: the real part of the elastodynamic Green's function normalized by its zero‐offset version. This makes it possible to take a deterministic approach rather than the statistical one on which the conventional spatial autocorrelation method is based and gives fruitful new aspects and perspectives. For example, the formula for the multi‐mode case is given and the possibility of exploration of two or three dimensional velocity structures is suggested.     

各向异性介质qP波传播描述II:分离纯模式标量波

在各向异性地震波场中,qP波与qS波常常是耦合在一起的.多分量地震数据处理中一个关键环节就是波型分离(即模式解耦),以纵波成分为主的常规单分量地震数据的成像则需要合理描述标量qP波的传播算子.本文作者曾构建了在运动学上同弹性波动方程等价,动力学上突出标量qP波的伪纯模式波动方程.为了彻底消除qS波残余,本文根据波矢量与qP波偏振矢量之间的偏差,提出从伪纯模式波场提取纯模式标量qP波的方法.数值分析展示了投影偏差算子在波数域和空间域的特征.基于不同复杂程度理论模型的试验结果表明,联合"伪纯模式传播算子"与"投影偏差校正"可为各向异性介质分离模式波场传播过程提供一种简便的描述工具.     

SH波地表散射的局域边界元模拟

本文提出了一种计算不规则起伏地形中SH波散射的有效方法——局域边界元法.本方法基于传统边界元法,为计算复杂地表散射问题提供了一种更加高效的解决方案.根据地震波满足的边界积分方程中牵引力格林函数的特性,我们将自由边界分解成水平部分和起伏部分.通过公式推导,可将水平部分的位移由起伏部分的位移通过格林函数线性叠加表示,因此只需对起伏部分的位移进行直接求解,从而极大地减少了待求解的未知数个数,显著提高了计算效率.通过与半圆形山谷SH波平面波入射的解析解比较,验证了方法的正确性.数值模型比较显示,局域边界元模拟结果与传统边界元数值解完全吻合,但是大幅提高了计算效率.因此,局域边界元法可以作为模拟不规则地形中地震波散射的有效工具.     

Simulation of improved pure P-wave equation in transversely isotropic media with a horizontal symmetry axis

Characterizing the expressions of seismic waves in elastic anisotropic media depends on multiparameters. To reduce the complexity, decomposing the P-mode wave from elastic seismic data is an effective way to describe the considerably accurate kinematics with fewer parameters. The acoustic approximation for transversely isotropic media is widely used to obtain P-mode wave by setting the axial S-wave phase velocity to zero. However, the separated pure P-wave of this approach is coupled with undesired S-wave in anisotropic media called S-wave artefacts. To eliminate the S-wave artefacts in acoustic waves for anisotropic media, we set the vertical S-wave phase velocity as a function related to propagation directions. Then, we derive a pure P-wave equation in transversely isotropic media with a horizontal symmetry axis by introducing the expression of vertical S-wave phase velocity. The differential form of new expression for pure P-wave is reduced to second-order by inserting the expression of S-wave phase velocity as an auxiliary operator. The results of numerical simulation examples by finite difference illustrate the stability and accuracy of the derived pure P-wave equation.     

Modified image algorithm to simulate seismic channel waves in 3D tunnel model with rugged free surfaces

The existence of rugged free‐surface three‐dimensional tunnel conditions in the coal seams, caused either by geological or mining processes, will inevitably influence wave propagation characteristics when the seismic waves go through the coal mines. Thus, a modified image algorithm has been developed to account for seismic channel waves propagating through this complicated topography with irregular free surfaces. Moreover, the seismic channel waves commonly exhibit damped and dispersive signatures, which is not only because of their own unique sandwich geometry of rock–coal–rock but also because of the viscoelastic behavior of coal. Considering the complexity of programming in three‐dimensional tunnel models with rugged free surfaces, an optimized vacuum grid search algorithm, enabling to model highly irregular topography and to compute efficiently, is also proposed when using high‐order staggered finite‐difference scheme to simulate seismic channel wave propagations in viscoelastic media. The numerical simulations are implemented to investigate the accuracy and stability of the method and the impact of coal's viscoelastic behavior on seismic channel wave propagation characteristics. The results indicate that the automatic vacuum grid search algorithm can be easily merged into high‐order staggered finite‐difference scheme, which can efficiently be applied to calculate three‐dimensional tunnel models with rugged free surfaces in the viscoelastic media. The simulation also suggests that the occurrence of a three‐dimensional tunnel with free surfaces has a remarkable influence on the seismic channel wave propagation characteristics and elastic energy distribution.     

3D高阶抛物Radon变换在不规则地震数据保幅重建中的应用

3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换（NHOPRT）地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而获得最接近完整规则数据的高阶抛物Radon变换域系数.在用SVD求解反演方程过程中,利用高阶抛物Radon变换算子在频率域为指数函数,具有线性可分解特性,将二维空间的高阶抛物Radon变换算子分解为两个独立的一维空间变换算子,减小了变换算子的矩阵大小,从而很大程度地提高了计算效率.理论模型和实际地震数据重建测试证明了本文方法的有效性以及实用性.     

The plurality effect of topographical irregularities on site seismic response

Topography can have significant effects on seismic ground response during an earthquake because topographic irregularities cause considerable differences between the seismic waves emitted by the source and the waves reaching the ground surface. When a seismic motion happens in a topographically irregular area, seismic waves are trapped and reflected between the topographic features. Therefore, the interaction between topographies can amplify seismic ground response. In order to reveal how interaction between topographies influences seismic response, several numerical finite element studies have been performed by using the ABAQUS program. The results show that topographic features a greater distance between the seismic source and the site would cause greater seismic motion amplification and is perceptible for the hills far away from the source and the ridges. Also, site acceleration response is impacted by surrounding topography further than site velocity and displacement response.