Boundary element simulation of scattering of elastic waves by 3-D cracks

Numerical modelling techniques are now becoming common for understanding the complicated nature of seismic wave propagation in fractured rock. Here the Indirect Boundary Element Method (IBEM) is applied to study scattering of elastic waves by cracks. The problem addressed in this paper is the diffraction of P and S waves by open 3-D cracks of arbitrary shape embedded in a homogeneous isotropic medium. The IBEM yields the value of the jump of displacements between opposite surfaces of the crack, often called Crack Opening Displacement (COD). This is used to evaluate the solution away from the crack. We use a multi-regional approach which consists of splitting a surface S into two identical surfaces S⁺ and S⁻ chosen such that the crack lies at the interface. The resulting integral equations are not hyper-singular and wave propagation within media that contain open cracks can be rigorously solved. In order to validate the method, we compare results of displacements of a penny-shaped crack for a vertical incident P-wave with the classic results by Mal (1970) obtaining excellent agreement. This comparison gives us confidence to study cases where no analytic solutions exist. Some examples of incidence of P or S waves upon cracks with various shapes are depicted and the salient aspects of the method are also discussed. Both frequency and time-domain results are included.     

Scholte waves on fluid-solid interfaces by means of an integral formulation

The present work shows the propagation of Scholte interface waves at the boundary of a fluid in contact with an elastic solid, for a broad range of solid materials. It has been demonstrated that by an analysis of diffracted waves in a fluid it is possible to infer the mechanical properties of the elastic solid medium, specifically, its propagation velocities. For this purpose, the diffracted wave field of pressures and displacements, due to an initial wave of pressure in the fluid, are expressed using boundary integral representations, which satisfy the equation of motion. The source in the fluid is represented by a Hankel’s function of second kind and zero order. The solution to this wave propagation problem is obtained by means of the Indirect Boundary Element Method, which is equivalent to the well-known Somigliana representation theorem. The validation of the results is carried out by using the Discrete Wave Number Method and the Spectral Element Method. Firstly, we show spectra of pressures that illustrate the behavior of the fluid for each solid material considered, then, we apply the Fast Fourier Transform to show results in time domain. Snapshots to exemplify the emergence of Scholte’s waves are also included.     

Rayleigh波入射下相邻山体地震响应边界元法模拟

采用间接边界元法（IBEM）,对Rayleigh波入射下两邻近山体的地震响应进行了定量分析。结果表明:与单个山体在Rayleigh波入射下的地震响应相比,两山地形的地震反应规律更为复杂,反应特征受控于入射波频率和山体间距等因素。总体上看:入射Rayleigh波频率较低时,两山间距对双山地形地震反应影响较大。低频波入射,受邻近山体影响,迎波面山体水平位移峰值约为入射波水平位移的6.3倍,放大效应可达单山作用的1.5倍。且山体山脚处竖向位移的频谱振荡更为剧烈,特定频率下山脚的竖向位移反应可达单山放大效应的1.6倍。较高频波入射时,迎波面一侧山体受邻近山体影响较小,且对波表现出明显的屏障效应,背波面山体地震反应强度被削弱。     

The discrete wave number formulation of boundary integral equations and boundary element methods: A review with applications to the simulation of seismic wave propagation in complex geological structures

We review the application of the discrete wave number method to problems of scattering of seismic waves formulated in terms of boundary integral equation and boundary element methods. The approach is based on the representation of the diffracting surfaces and interfaces of the medium by surface distributions of sources or by boundary source elements, the radiation from which is equivalent to the scattered wave field produced by the diffracting boundaries. The Green's functions are evaluated by the discrete wave number method, and the boundary conditions yield a linear system of equations. The inversion of this system allows the calculation of the full wave field in the medium. We investigate the accuracy of the method and we present applications to the simulation of surface seismic surveys, to the diffraction of elastic waves by fractures, to regional crustal wave propagation and to topographic scattering.     

Visualization of whole mantle propagation of seismic shear energy using normal mode summation

Correct representation of seismic waveforms propagating through the mantle from a 600 km deep earthquake is presented using graphic interpolation between synthetic seismograms computed across a grid of mantle depths and distances. All torsional normal modes with periods above 12 s are summed to create 72,846 seismograms at depths between the surface and the core-mantle boundary. The resulting time snapshots show the manner by which seismic shear energy propagates around the core away from the source.     

The interaction of elastic waves with a solid-liquid interface,with applications to the core-mantle boundary

The paper discusses basic ideas and principles underlying methods, which have proved useful in the interpretation of diffraction and scattering phenomena by a smooth or slightly rough solid-liquid interface. Generally, the wave interaction may be formulated as an excitation problem; it involves (1) finding an equivalent dislocation or source distribution on the interface, and (2) evaluating the excited wave field. These steps are taken through perturbation theory and/or adopting the appropriate generalization of ray theory. In this context an explicit form of Green's function is also given. The methods have been applied to the core-mantle boundary, with a view toward the interpretation of recent data: (a) diffracted P and S waves around the core (in particular, their attenuation); (b) scattered short-period core phases (in particular, precursors to PKP and PKKP). Other types of wave interaction, and implications for models of the core-mantle boundary structure, are briefly mentioned.Paper presented at the EGS/ESC workshop on Generation and propagation of seismic waves in Neustadt, Federal Republic of Germany, August 1978.     

Propagation of seismic waves through liquefied soils

To predict the earthquake response of saturated porous media it is essential to correctly simulate the generation, redistribution, and dissipation of excess pore water pressure during and after earthquake shaking. To this end, a reliable numerical tool requires a dynamic, fully coupled formulation for solid–fluid interaction and a versatile constitutive model. Presented in this paper is a 3D finite element framework that has been developed and utilized for this purpose. The framework employs fully coupled dynamic field equations with a u–p–U formulation for simulation of pore fluid and solid skeleton interaction and a SANISAND constitutive model for response of solid skeleton. After a detailed verification and validation of the formulation and implementation of the developed numerical tool, it is employed in the seismic response of saturated porous media. The study includes examination of the mechanism of propagation of the earthquake-induced shear waves and liquefaction phenomenon in uniform and layered profiles of saturated sand deposits.     

Scattering of elastic waves in cracked media using a finite difference method

This paper presents a simple, flexible way of introducing stress-free boundary conditions for including cracks and cavities in 2D elastic media by a finite difference method (FDM). The surfaces of cracks and cavities are discretized in a staircase on a rectangular grid scheme. When zero-stress is applied to free surfaces, the resulting finite difference schemes require a set of adjacent fictitious points. These points are classified based on the geometry of the free surface and their displacement is computed as a prior step to later calculation of motion on the crack surface. The use of this extra line of points does not involve a significant drain on computational resources. However, it does provide explicit finite difference schemes and the construction of displacement on the free surfaces by using the correct physical boundary conditions. An accuracy analysis compares the results to an analytical solution. This quantitative analysis uses envelope and phase misfits. It estimates the minimum number of points per wavelength necessary to achieve suitable results. Finally, the method is employed to compute displacement in various models with cavities in the P-SV formulation. The results show suitable construction of the reflected P and S waves from the free surface as well as diffraction produced by these cavities.     

Reflection of waves from the boundary of a random elastic semi-infinite medium

In this paper the smooth perturbation technique is employed to investigate the problem of reflection of waves incident on the plane boundary of a semi-infinite elastic medium with randomly varying inhomogeneities. Amplitude ratios have been obtained for various types of incident and reflected waves. It has been shown that an incidentSH orSV type of wave gives rise to reflectedSH, P andSV waves, the main components beingSH andP, SV in the respective cases. The reflected amplitudes have been calculated depending upon the randomness of the medium to the square of the small quantity , where measures the deviation of the medium from homogeneity. An incidentP-type wave produces mainly aP component and also a weakSH component to the order of ². The reflected amplitudes obtainable for elastic media are also altered by terms of the same order. The direction of the reflected wave is influenced by randomness in some cases.     

Analytical expressions for the borehole seismic signatures excited by an explosive source buried in an elastic half‐space

Borehole guided waves that are excited by explosive sources outside of the borehole are important for interpreting borehole seismic surveys and for rock property inversion workflows. Borehole seismograms are typically modelled using numerical methods of wave propagation. In order to benchmark such numerical algorithms and partially to interpret the results of modelling, an analytical methodology is presented here to compute synthetic seismograms. The specific setup is a wavefield emanating from a monopole point source embedded within a homogeneous elastic medium that interacts with a fluid‐filled borehole and a free surface. The methodology assumes that the wavelength of the seismic signal is much larger than the borehole radius. In this paper, it is supposed that there is no poroelastic coupling between the formation and the borehole. The total wavefield solution consists of P, PP, and PS body waves; the surface Rayleigh wave; and the low‐frequency guided Stoneley wave (often referred as the tube wave) within the borehole. In its turn, the tube wave consists of the partial responses generated by the incident P‐wave and the reflected PP and PS body waves at the borehole mouth and by the Rayleigh wave, as well as the Stoneley wave eigenmode. The Mach tube wave, which is a conic tube wave, additionally appears in the Mach cone in a slow formation with the tube‐wave velocity greater than the shear one. The conditions of appearance of the Mach wave in a slow formation are formulated. It is shown that the amplitude of the Mach tube wave strongly depends on Poisson's ratio of the slow surrounding formation. The amplitude of the Mach tube wave exponentially decreases when the source depth grows for weakly compressible elastic media with Poisson's ratio close to 0.5 (i.e., saturated clays and saturated clay soils). Asymptotic expressions are also provided to compute the wavefield amplitudes for different combinations of source depth and source‐well offset. These expressions allow an approximate solution of the wavefield to be computed much faster (within several seconds) than directly computing the implicit integrals arising from the analytical formulation.     

Experimental observations of elastic wave radiation characteristics from tensile cracks and pre-existing shear faults

Elastic waves were recorded from tensile cracks and prepared shear faults in the laboratory. The velocity of tensile crack propagation was measured. The duration of crack propagation was correlated with the signal pulse duration. Experiments with prepared shear faults showed that obstacles on the fault change the pulse form and the relation between wave amplitude (magnitude) and seismic moment.     

Scattering of seismic waves by cracks in multi-layered geological regions: I. Mechanical model

In this work, a hybrid boundary integral equation method (BIEM) is developed, based on both displacement and hypersingular traction formulations, for the analysis of time-harmonic seismic waves propagating through cracked, multi-layered geological regions with surface topography and under plane strain conditions. Specifically, the displacement-based BIEM is used for a multi-layered deposit with interface cracks, while the regularized, traction-based BIEM is used when internal cracks are present within the layers. The standard uni-dimensional boundary element with parabolic shape functions is employed for discretizing the free surface and the layer interfaces, while special discontinuous boundary elements are placed near the crack tips to model the asymptotic behaviour of both displacements and tractions. This formulation yields displacement amplitudes and phase angles on the free surface of a geological deposit, as well as stress intensity factors near the tips of the cracks. Finally, in the companion paper, numerical results are presented which show that both scattered wave and stress concentration fields are sensitive to the incidence seismic wave parameters and to specific site conditions such as surface topography, layering, the presence of cracks and crack interaction.     

Propagation of partially coherent non-stationary random waves in a viscoelastic layered half-space

This paper presents a highly accurate method based on the precise integration method (PIM) and on the pseudo excitation method (PEM). The method computes the propagation behaviour of partially coherent non-stationary random waves in a viscoelastic, transversely isotropic solid, which consists of a multi-layered soil resting on a homogeneous semi-infinite space. The excitation source is a local rupture between two layers, which causes a partially coherent non-stationary random field. The analysis of non-stationary random wave propagation is transformed into that for deterministic waves by using PEM. The resulting governing equations in the frequency-wavenumber domain are linear ordinary differential equations, which are solved very precisely by using PIM. The evolutionary power spectral densities of the ground level responses are investigated and some typical earthquake phenomena are explained.     

Numerical modelling of wave propagation in anisotropic soil using a displacement unit-impulse-response-based formulation of the scaled boundary finite element method

An efficient method for modelling the propagation of elastic waves in unbounded domains is developed. It is applicable to soil–structure interaction problems involving scalar and vector waves, unbounded domains of arbitrary geometry and anisotropic soil. The scaled boundary finite element method is employed to derive a novel equation for the displacement unit-impulse response matrix on the soil–structure interface. The proposed method is based on a piecewise linear approximation of the first derivative of the displacement unit-impulse response matrix and on the introduction of an extrapolation parameter in order to improve the numerical stability. In combination, these two ideas allow for the choice of significantly larger time steps compared to conventional methods, and thus lead to increased efficiency. As the displacement unit-impulse response approaches zero, the convolution integral representing the force–displacement relationship can be truncated. After the truncation the computational effort only increases linearly with time. Thus, a considerable reduction of computational effort is achieved in a time domain analysis. Numerical examples demonstrate the accuracy and high efficiency of the new method for two-dimensional soil–structure interaction problems.     

Modifications of the dispersion relations for surface waves in a layer on a half-space

Summary Dispersion relations for Love and Rayleigh waves in a layer on a half-space are modified by introducing quadratic slownesses instead of velocities. The advantages of this approach are demonstrated on analytical formulae for computing the group velocity.     

3D seismic response of a limited valley via BEM using 2.5D analytical Green's functions for an infinite free-rigid layer

This paper presents analytical solutions for computing the 3D displacements in a flat solid elastic stratum bounded by a rigid base, when it is subjected to spatially sinusoidal harmonic line loads. These functions are also used as Greens functions in a boundary element method code that simulates the seismic wave propagation in a confined or semi-confined 2D valley, avoiding the discretization of the free and rigid horizontal boundaries.The models developed are then used to simulate wave propagation within a rigid stratum and valleys with different dimensions and geometries, when struck by a spatially sinusoidal harmonic vertical line load. Simulations are performed in the frequency domain, for varying spatial wave numbers in the axial direction of the valley. Time results are obtained by means of inverse Fourier transforms, to help understand how the geometry of the valley may affect the variation of the displacement field.     

Analysis of Rayleigh waves in saturated porous elastic media by finite element method

A numerical procedure for the analysis of Rayleigh waves in saturated porous elastic media is proposed by use of the finite element method. The layer stiffness matrix, the layer mass matrix and the layer damping matrix in a layered system are presented for the discretized form of the solid-fluid equilibrium equation proposed by Biot. In order to consider the influence of the permeability coefficient on the behavior of Rayleigh waves, attention is focused on the following states: ‘drained’ state, ‘undrained’ state and the states between two extremes of ‘drained’ and ‘undrained’ states. It is found from computed results that the permeability coefficient exerts a significant effect on dispersion curves and displacement distributions of Rayleigh waves in saturated porous media.     

Transmission fluctuations across an array and heterogeneities in the crust and upper mantle

Adopting the spectral approach, we derive the formulation of angular coherence and transverse coherence of transmission fluctuations. Our derivation and results provide new insight on transmission fluctuation analysis. A review of research work on fluctuation analysis using observations at large seismic arrays such as LASA and NORSAR-follows. We point out that the model of a single-layer Gaussian medium cannot explain the angular coherence of NORSAR data and a more general model of a non-Gaussian, multi-scale, vertically inhomogeneous random media is needed. The model of a two-layer power-law medium proposed by Flatté and Wu is among the simplest of such models.     

The reflection of plane waves in a poroelastic half-space saturated with inviscid fluid

This paper discusses surface displacements, surface strain, rocking, and energy partitioning during reflection-of-plane waves in a fluid-saturated poroelastic half-space. The medium is modeled by Biot's theory, and is assumed to be saturated with inviscid fluid. A linear porosity-modulus relation based on experimental data on sandstones is used to determine the material parameters for Biot's model. Numerical results in terms of angle of incident waves and Poisson's ratio are illustrated for various porosities and degrees of solid frame stiffness. The results show that the amount of solid frame stiffness controls the response of a fluid-saturated porous system. A poroelastic medium with essentially dry-frame stiffness behaves like an elastic medium, and the influence of pore fluid increases as dry-frame stiffness is reduced. The effects of a second P-wave become noticeable in poroelastic media with low dry-frame stiffness.