A semi-analytical boundary element method for scattering of waves in a half space

A boundary element method using finite strip discretization is proposed to investigate the three-dimensional scattering behaviour of waves in a half space. The scatterers examined are an irregular surface and a cavity in a half space. The numerical examples demonstrate the remarkably good accuracy of the results and validate the application of the present method to solve three-dimensional scattering problems.     

A local time-domain transmitting boundary for simulating cylindrical elastic wave propagation in infinite media

Finite element simulation of the time-dependent wave propagation in infinite media requires enforcing the transmitting boundary to replace the truncated far-field infinite domain so as to model the effect of the wave radiation towards infinity. This paper proposed a novel local time-domain transmitting boundary for simulating the cylindrical elastic wave radiation problem. This boundary is a mechanical model consisting of the spring, dashpot and mass elements, with the auxiliary degrees of freedom introduced, which is dynamically stable and easily implemented into the commercial finite element codes. Numerical analysis of the cylindrical elastic wave radiation problem indicates that the proposed transmitting boundaries with the order N=3 for cylindrical P and SV waves and with the order N=4 for cylindrical SH wave have very high accuracy, even when the artificial boundary at wave source. The proposed transmitting boundary with order N=0 can be applied approximately to the general two-dimensional infinite elastic wave problems that contain the more complex outgoing wave fields at artificial boundary than the cylindrical waves. The plane-strain Lamb problem is analyzed with the acceptable engineering accuracy achieved. On the other hand, the proposed transmitting boundary with higher order can be a tool to localize the temporal convolution that appears in an exact time-domain transmitting boundary for the general infinite wave problems. This potential applicability is mentioned.     

Surface wave propagation in fiber-reinforced anisotropic elastic layer between liquid saturated porous half space and uniform liquid layer

Surface wave propagation in fiber-reinforced anisotropic elastic layer between a liquid saturated porous half space and a uniform liquid layer is considered. Equation of motion and suitable boundary conditions give rise to a dispersion equation in the form of a ninth order determinant. Phase velocity and group velocity of a particular model have been studied.     

Combined paraxial-consistent boundary conditions finite element model for simulating wave propagation in elastic half-space media

In this paper, a finite element model of a soil island is coupled to both a consistent transmitting boundary and a paraxial boundary, which are then used to model the propagation of waves in semi-infinite elastic layered media. The formulation is carried out in the frequency domain while assuming plane strain conditions. It is known that a discrete model of this type, while providing excellent results for a wide range of physical parameters in the context of a half-space problem, may deteriorate rapidly at low frequencies of excitation. This is so because at low frequencies the various waves in the model eventually attain characteristic wavelengths which exceed the distance of the bottom boundary, which then causes that boundary to fail. Also, the paraxial boundaries themselves break down at very low frequencies. In this paper, this difficulty is overcome and the model׳s performance is improved upon dramatically by incorporating an artificial buffer layer sandwiched between the bottom of the soil medium and the underlying elastic half-space. Applications dealing with rigid foundations resting on homogenous or layered half-space media are shown to exhibit significant improvement. Following extensive simulations, clear guidelines are provided on the performance of the coupled model and an interpretation is given on the engineering significance of the findings. Finally, clear recommendations are provided for the practical use of the proposed modelling strategy.     

Transient analysis of wave propagation in non‐homogeneous elastic unbounded domains by using the scaled boundary finite‐element method

The scaled boundary finite‐element method has been developed for the dynamic analysis of unbounded domains. In this method only the boundary is discretized resulting in a reduction of the spatial dimension by one. Like the finite‐element method no fundamental solution is required. This paper extends the scaled boundary finite‐element method to simulate the transient response of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. To reduce the number of degrees of freedom and the computational cost, the technique of reduced set of base functions is applied. The scaled boundary finite‐element equation for an unbounded domain is reformulated in generalized coordinates. The resulting acceleration unit‐impulse response matrix is obtained and assembled with the equation of motion of standard finite elements. Numerical examples of non‐homogeneous isotropic and transversely isotropic unbounded domains demonstrate the accuracy of the scaled boundary finite‐element method. Copyright © 2006 John Wiley & Sons, Ltd.     

弹性波数值模拟的延迟边界方法

在地震波场的波动方程数值模拟中,由于计算量的限制,必须加入人为的边界,使模拟计算可以在一定的空间范围内进行. 由于边界节点上的波场值不能像模拟区域内部的节点一样使用中心差分来计算,使其计算精度大大降低,从而产生边界反射. 为了消除边界反射,本文提出了延迟边界方法,根据弹性波在传播方向上等距离质点的等相位延迟性质和振幅衰减特性,由内部波场的时空分布,推算出边界波场的相位延迟的大小和振幅衰减系数,从而提高边界节点上的波场值计算精度,消除边界反射的产生.     

A viscous-spring transmitting boundary for cylindrical wave propagation in saturated poroelastic media

Based on the u–U formulation of Biot equation and the assumption of zero permeability coefficient, a viscous-spring transmitting boundary which is frequency independent is derived to simulate the cylindrical elastic wave propagation in unbounded saturated porous media. By this viscous-spring boundary the effective stress and pore fluid pressure on the truncated boundary of the numerical model are replaced by a set of spring, dashpot and mass elements, and its simplified form is also given. A u–U formulation FEA program is compiled and the proposed transmitting boundaries are incorporated therein. Numerical examples show that the proposed viscous-spring boundary and its simplified form can provide accurate results for cylindrical elastic wave propagation problems with low or intermediate values of permeability or frequency content. For general two dimensional wave propagation problems, spuriously reflected waves can be greatly suppressed and acceptable accuracy can still be achieved by placing the simplified boundary at relatively large distance from the wave source.     

多人工波速优化透射边界在谱元法地震波动模拟中的应用

将作者最近发展的多人工波速优化透射边界(记为ca j-MTF)应用于高精度谱元法的地震波动模拟中,并与经典的廖氏透射(MTF)边界、完美匹配层(PML)边界、黏弹性边界以及一阶旁轴近似边界进行了比较分析.结果显示:①ca j-MTF边界与MTF边界在形式上非常接近,它继承了后者公式简单、易于实现、精度可控、计算量低以及...     

饱和多孔介质近场波动分析的一种黏弹性人工边界

采用有限模型数值方法求解能量开放系统中的波动问题时,虚拟人工边界的处理方式对计算结果的准确性和精度具有重要的影响.本文针对无限域饱和多孔介质中波传播问题的人工边界处理方式进行了研究,提出了饱和多孔介质近场波动分析的一种黏弹性人工边界处理方法.在考虑多孔介质中固相和液相的相互作用的情况下,通过在人工边界处分别施加反映固相和液相介质波传播效应的弹簧及阻尼来模拟饱和多孔介质中波的能量辐射效应影响.算例表明,本文建议的黏弹性人工边界具有较好的模拟效果.     

A two-dimensional horizontal wave propagation and mud mass transport model

The present study offers a two-dimensional horizontal wave propagation and morphodynamic model for muddy coasts. The model can be applied on a general three-dimensional bathymetry of a soft muddy coast to calculate wave damping, fluid mud mass transport and resulting bathymetry change under wave actions. The wave propagation model is based on time-dependent mild slope equations including the wave energy dissipation due to the wave-mud interaction of bottom mud layers as well as the combined effects of the wave refraction, diffraction and breaking. The constitutive equations of the visco-elastic–plastic model are adopted for the rheological behavior of fluid mud. The mass transport velocity within the fluid mud layer is calculated combining the Stokes’ drift, the mean Eulerian velocity and the gravity-driven mud flow. The results of the numerical model are compared against a series of conducted wave basin experiments, wave flume experiments and field observations. Comparisons between the computed results with both the field and laboratory data reveal the capability of the proposed model to predict the wave transformation and mud mass transport.     

Hybrid absorbing boundary condition for three-dimensional elastic wave modeling

Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer (PML) boundary condition requires special treatment for the absorbing zone, and in three-dimensional (3D) modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition (HABC) has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional (2D) modeling into 3D modeling based on the 1st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.     

A discontinuous Galerkin method for seismic wave propagation in coupled elastic and poroelastic media

Numerical simulation in coupled elastic and poroelastic media is important in oil and gas exploration. However, the interface between elastic and poroelastic media is a challenge to handle. In order to deal with the coupled model, the first-order velocity–stress wave equations are used to unify the elastic and poroelastic wave equations. In addition, an arbitrary high-order discontinuous Galerkin method is used to simulate the wave propagation in coupled elastic–poroelastic media, which achieves same order accuracy in time and space domain simultaneously. The interfaces between the two media are explicitly tackled by the Godunov numerical flux. The proposed forms of numerical flux can be used efficiently and conveniently to simulate the wave propagation at the interfaces of the coupled model and handle the absorbing boundary conditions properly. Numerical results on coupled elastic–poroelastic media with straight and curved interfaces are compared with those from a software that is based on finite element method and the interfaces are handled by boundary conditions, demonstrating the feasibility of the proposed scheme in dealing with coupled elastic–poroelastic media. In addition, the proposed method is used to simulate a more complex coupled model. The numerical results show that the proposed method is feasible to simulate the wave propagation in such a media and is easy to implement.     

Measuring wave propagation characteristics in artificial sand-clay mixtures

Ultrasonic compressional (V _p ) and shear (V _s ) velocities have been measured on artificial sand-clay mixtures. The measurements were carried out in a drained triaxial load cell using a pulse transition method. The measuring device was equiped with a waveform storage facility. The investigated mixtures consisted mainly of kaolinite and quartz sand. Some mixtures also contained Na-montmorillonite, illites or quartz-flour. The acoustic behaviour was observed during a pressure increase up to 72 MPa vertical and 36 MPa horizontal pressure. At a given pressure,V _p andV _s in pure sand turned out to be similar to those in pure kaolinite. As predicted by the sand-clay model of Marion (1990), a velocity maximum corresponds to a minimum in total porosity. This porosity minimum marks the transition from a clayey sand to a sandy clay. It is not only reflected in bothV _p andV _s , but also in the quality of the received pulse. The effective tension of the received signal during 20µs after the first arrival, was used as an indication for P-wave pulse attenuation. This apparent attenuation decreases with increasing clay content and increases with increasing porosity. It is shown that clay mineralogy does not measurably affect wave velocities in clayey sands.     

一种模拟非均匀介质中弹性波传播新的k-space方法

传统的伪谱（PS）方法,采用傅里叶变换（FT）计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分（FD）近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.

     

Transient response of a homogeneous half space with due regard for displacement currents

The solution for the half-space model is represented directly in the time domain as computationally stable convolution integrals. The influence of the geoelectric parameters of the earth and transmitter current waveform are then investigated for both infinitesimal and finite-dimensional transmitter loops. Simple empirical formulae are derived to account for the finite duration of the transmitter current turn off time.The whole transient process is divided into three essentially different stages: the propagation stage, the intermediate stage and the diffusion stage. The first is characterized by extremely complicated signal behavior. Apparently, interpretation of the field data using any kind of model fitting inversion algorithm is impossible in this stage. The diffusion stage virtually coincides with that used in the quasi-static case and is, therefore, unsuitable for detecting the dielectric properties of the earth. The intermediate stage is, thus, the only possible time range in which the dielectric properties can be detected using the dynamic characteristics of the signal.The duration of each stage is evaluated depending on the geoelectric parameters of the earth for different transmitter current waveforms.     

弹性波传播中由微结构相互作用导致的尺度效应分析

相比于传统弹性波动方程,非对称弹性波动方程增加的独立自由项,包含有介质特征尺度参数.基于非对称弹性波动方程,可以分析弹性波传播中,由介质内微孔缝隙结构相互作用所导致的地震波传播尺度效应.本文从介质应变能密度函数出发,并结合几何方程和平衡方程,给出修正偶应力理论下的非对称弹性波动方程以及对应的非对称SH型横波波动方程的数...     

Dispersion ofSH body waves in a continuously stratified elastic half space

Summary The dispersion ofSH body waves in a continously stratified isotropic medium is studied theoretically. Limiting forms of the displacement are given for zero, small and large values of the horizontal component of the wave number. It has been shown that below a certain characteristic frequency the dependence of the displacement on the distance from the free surface is monotonic instead of periodic.     

A new series solution method for two-dimensional elastic wave scattering along a canyon in half-space

This paper presents a semi-analytical method for studying the two-dimensional problem of elastic wave scattering by surface irregularities in a half-space. The new method makes use of the member of a c-completeness family of wave functions to construct the scattering fields, and then applies equal but opposite tractions to those of the foregoing constructed scattering fields on the horizontal surface of the half-space to produce additional scattering fields. These additional scattering fields are a series of Lamb's solutions. Thus the whole scattering field constructed in the series automatically satisfies the Navier equations, the condition of zero traction on the half-space surface, and the radiation boundary conditions at infinity. Using the traction-free conditions along the canyon surface, the coefficients of the series solutions are determined via a least-squares method. For incident P, SV, and Rayleigh waves, the numerical results are presented for the scattering displacements in the vicinity of a semi-circular canyon in the half-space.     

A high-order time-domain transmitting boundary for cylindrical wave propagation problems in unbounded saturated poroelastic media

Based on the u–p formulation of Biot equation with an assumption of zero permeability coefficient, a high-order transmitting boundary is derived for cylindrical elastic wave propagation in infinite saturated porous media. By this transmitting boundary the total stresses on the truncated boundaries of a numerical model, such as a finite element model, are replaced by a set of spring, dashpot and mass elements, with some additionally introduced auxiliary degrees of freedom. The transmitting boundaries are incorporated into the DIANA SWANDYNE II program and an unconditionally stable implicit time integration algorithm is adopted. Despite the assumption made in the derivation of the transmitting boundary, numerical examples show that it can provide highly accurate results for cylindrical elastic wave propagation problems in infinite saturated porous medium in case the u–p formulation is applicable. Although the direct applications of the proposed transmitting boundary to general two dimensional wave problems in infinite saturated porous media are not highly accurate, acceptable accuracy can still be achieved by placing the transmitting boundary at relatively large distance from the wave source.     

Effects of source and cavity depths on wave fields in homogeneous half spaces

In homogeneous half spaces excited by small buried, spherically symmetric P-wave sources, Rayleigh waves (R-waves) could be generated. The effects of the source depth on the induced wave pattern and propagation behavior of R-waves are analyzed using the thin layer method. When a cavity is present in a homogeneous half space, R-waves could be formed in the scattered wave field. It is found that the energy of R-waves in the incident (direct) surface wave-field is related to the ratio of the source depth to the wavelength of R-waves; R-waves have relatively strong energy when the ratio is less than 1; the buried source induced R-waves approximately travel at the velocity of the planar R-waves in the range of the offset beyond about one wavelength; the energy of R-waves in the back-scattered surface wave-field depends on the ratio of the depth of cavity to the wavelength of R-waves; and for the case of the cavity presented at depths less than one wavelength, R-waves can be clearly observed in the back-scattered wave field. The results are helpful for selecting the source depth and the frequency component in seismic surveys and interpreting both the incident wave and the scattered wave patterns.