Finite-difference modeling of surface waves in poroelastic media and stress mirror conditions

During seismic wave propagation on a free surface, a strong material contrast boundary develops in response to interference by P- and S- waves to create a surface-wave phenomenon. To accurately determine the effects of this interface on surface-wave propagation, the boundary conditions must be accurately modeled. In this paper, we present a numerical approach based on the dynamic poroelasticity for a space–time-domain staggered-grid finite-difference simulation in porous media that contain a free-surface boundary. We propose a generalized stess mirror formulation of the free-surface boundary for solids and fluids in porous media for the grid mesh on which lays the free-surface plane. Its analog is that used for elastic media, which is suitable for precise and stable Rayleigh-type surface-wave modeling. The results of our analysis of first kind of Rayleigh （R1） waves obtained by this model demonstrate that the discretization of the mesh in a similar way to that for elastic media can realize stable numerical solutions with acceptable precision. We present numerical examples demonstrating the efficiency and accuracy of our proposed method.

Dynamic stiffness matrix of a poroelastic multi-layered site and its Green's functions

Few studies of wave propagation in layered saturated soils have been reported in the literature.In this paper,a general solution of the equation of wave motion in saturated soils,based on one kind of practical Blot's equation, was deduced by introducing wave potentials.Then exact dynamic-stiffness matrices for a poroelastic soil layer and half- space were derived,which extended Wolf's theory for an elastic layered site to the case of poroelasticity,thus resolving a fundamental problem in the field of wave propagation and soil-structure interaction in a poroelastic layered soil site.By using the integral transform method,Green's functions of horizontal and vertical uniformly distributed loads in a poroelastic layered soil site were given.Finally,the theory was verified by numerical examples and dynamic responses by comparing three different soil sites.This study has the following advantages:all parameters in the dynamic-stiffness matrices have explicitly physical meanings and the thickness of the sub-layers does not affect the precision of the calculation which is very convenient for engineering applications.The present theory can degenerate into Wolf's theory and yields numerical results approaching those for an ideal elastic layered site when porosity tends to zero.

Diffraction of plane SV waves by a cavity in poroelastic half-space

This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green＇s functions of compressive and shear wave sources are derived based on Biot＇s theory. The scattered waves are constructed using fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by： （1） checking the satisfaction extent of the boundary conditions, （2） comparing the degenerated solutions of a single-phased case with well- known solutions, and （3） examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.

Screening of plane S waves by an array of rigid piles in poroelastic soil

An array of rigid piles used as a screening barrier for plane shear (S) waves is investigated in a homogeneous unbounded space. The dynamic poroelastic theory of Biot is employed, under the assumption of an incompressible solid grain. Using Fourier-Bessel series, the problem of multiple scattering is solved by imposing continuity conditions and equilibrium conditions at the soil-pile interfaces with the translational addition theorem. A parametric analysis is conducted to investigate the influence of the permeability of poroelastic soil, separation between piles, number of piles and frequency of incident waves on screening effectiveness of the barrier, and the results are compared with those in an elastic soil medium. Computed results show that the intrinsic permeability of the soil medium displays an apparent effect on the screening of plane S waves.

Green's functions for uniformly distributed loads acting on an inclined line in a poroelastic layered site

Based on one type of practical Biot＇s equation and the dynamic-stiffness matrices ofa poroelastic soil layer and half-space, Green＇s functions were derived for unitformly distributed loads acting on an inclined line in a poroelastie layered site. This analysis overcomes significant problems in wave scattering due to local soil conditions and dynamic soil-structure interaction. The Green＇s functions can be reduced to the case of an elastic layered site developed by Wolf in 1985. Parametric studies are then carried out through two example problems.

Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅰ): Formulation

This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green＇s functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot＇s theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.

Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅱ): Numerical results and discussion

This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated （either drained or undrained） cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.

Green's function of multi-layered poroelastic half-space for models of ground vibration due to railway traffic

This study proposes a Green＇s function, an essential representation of water-saturated ground under moving excitation, to simulate ground borne vibration from trains. First, general solutions to the governing equations of poroelastic medium are derived by means of integral transform. Secondly, the transmission and reflection matrix approach is used to formulate the relationship between displacement and stress of the stratified ground, which results in the matrix of the Green＇s function. Then the Green＇s function is combined into a train-track-ground model, and is verified by typical examples and a field test. Additional simulations show that the computed ground vibration attenuates faster in the immediate vicinity of the track than in the surrounding area. The wavelength of wheel-rail unevenness has a notable effect on computed displacement and pore pressure. The variation of vibration intensity with the depth of ground is significantly influenced by the layering of the strata soil. When the train speed is equal to the velocity of the Rayleigh wave, the Mach cone appears in the simulated wave field. The proposed Green＇s function is an appropriate representation for a layered ground with shallow ground water table, and will be helpful to understand the dynamic responses of the ground to complicated moving excitation.

THE NON-AXISYMMETRICAL DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED POROELASTIC MEDIA

The Blot's wave equations of transversely isotropic saturated poroelastic media excited hy non-axisymmetrical harmonic source were solved by means of Fourier expansion and Hankel transform. Then the components of total stress in porous media are expressed with the solutions of Biot's wave equations. The method of research on non-axisymmetrical dynamic response of saturated porous media is discussed, and a numerical result is presented.

Soil pressure and pore pressure for seismic design of tunnels revisited: considering water-saturated, poroelastic half-space

This paper describes a systematic study on the fundamental features of seismic soil pressure on underground tunnels, in terms of its magnitude and distribution, and further identifi es the dominant factors that signifi cantly infl uence the seismic soil pressure. A tunnel embedded in water-saturated poroelastic half-space is considered, with a large variety of model and excitation parameters. The primary features of both the total soil pressure and the pore pressure are investigated. Taking a circular tunnel as an example, the results are presented using a fi nite element-indirect boundary element(FE-IBE) method, which can account for dynamic soil-tunnel interaction and solid frame-pore water coupling. The effects of tunnel stiffness, tunnel buried depth and input motions on the seismic soil pressure and pore pressure are also examined. It is shown that the most crucial factors that dominate the magnitude and distribution of the soil pressure are the tunnel stiffness and dynamic soil-tunnel interaction. Moreover, the solid frame-pore water coupling has a prominent infl uence on the magnitude of the pore pressure. The fi ndings are benefi cial to obtain insight into the seismic soil pressure on underground tunnels, thus facilitating more accurate estimation of the seismic soil pressure.

Amplification of in-plane seismic ground motion by group cavities in layered half-space (Ⅱ): with saturated poroelastic soil layers

As the continuation study on amplification of in-plane seismic ground motion by underground group cavities in layered half-space, this study extends to the case of poroelastic half-space with dry poroelastic and saturated poroelastic soil layers. The influence of poroelastic layers on the amplification of seismic ground motion is studied both in frequency domain and time domain using indirect boundary element method （IBEM）. It is shown that for the example of a saturated poroelastic site in Tianjin under the excitation of Taft wave and E1 Centro wave, the amplification of seismic ground motion in poroelastic case is slightly smaller than that in the elastic case, and the amplification of PGA （peak ground acceleration） and its PRS （peak response spectrum）.. can be increased up to 38.8% and 64.6%; the predominant period of response spectra in poroelastic case becomes shorter to some extent compared with that in the elastic case. It is suggested that the effect of underground group cavities in poroelastic half-space on design seismic ground motion should be considered.

Effect of Rotation on the Propagation of Waves in Hollow Poroelastic Circular Cylinder with Magnetic Field

Employing Biot’s theory of wave propagation in liquid saturated porous media,the effect of rotation and magnetic field on wave propagation in a hollow poroelastic circular of infinite extent is investigated.An exact closed form solution is presented.General frequency equations for propagation of poroelastic cylinder are obtained when the boundaries are stress free.The frequencies are calculated for poroelastic cylinder for different values of magnetic field and rotation.Numerical results are given and illustrated graphically.The results indicate that the effect of rotation,and magnetic field are very pronounced.Such a model would be useful in large-scale parametric studies of mechanical response.

Seismic AVO statistical inversion incorporating poroelasticity

Seismic amplitude variation with offset(AVO) inversion is an important approach for quantitative prediction of rock elasticity,lithology and fluid properties.With Biot-Gassmann's poroelasticity,an improved statistical AVO inversion approach is proposed.To distinguish the influence of rock porosity and pore fluid modulus on AVO reflection coefficients,the AVO equation of reflection coefficients parameterized by porosity,rock-matrix moduli,density and fluid modulus is initially derived from Gassmann equation and critical porosity model.From the analysis of the influences of model parameters on the proposed AVO equation,rock porosity has the greatest influences,followed by rock-matrix moduli and density,and fluid modulus has the least influences among these model parameters.Furthermore,a statistical AVO stepwise inversion method is implemented to the simultaneous estimation of rock porosity,rock-matrix modulus,density and fluid modulus.Besides,the Laplace probability model and differential evolution,Markov chain Monte Carlo algorithm is utilized for the stochastic simulation within Bayesian framework.Models and field data examples demonstrate that the simultaneous optimizations of multiple Markov chains can achieve the efficient simulation of the posterior probability density distribution of model parameters,which is helpful for the uncertainty analysis of the inversion and sets a theoretical fundament for reservoir characterization and fluid discrimination.

Diffraction of plane P waves around an alluvial valley in poroelastic half-space

Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green＇s fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.

Teleseismic waves reveal anisotropic poroelastic response of wastewater disposal reservoir

Connecting earthquake nucleation in basement rock to fluid injection in basal,sedimentary reservoirs,depends heavily on choices related to the poroelastic properties of the fluid-rock system,thermo-chemical effects notwithstanding.Direct constraints on these parameters outside of laboratory settings are rare,and it is commonly assumed that the rock layers are isotropic.With the Arbuckle wastewater disposal reservoir in Osage County,Oklahoma,high-frequency formation pressure changes and collocated broadband ground velocities measured during the passing of large teleseismic waves show a poroelastic response of the reservoir that is both azimuthally variable and anisotropic;this includes evidence of static shifts in pressure that presumably relate to changes in local permeability.The azimuthal dependence in both the static response and shear coupling appears related to tectonic stress and strain indicators such as the orientations of the maximum horizontal stress and faults and fractures.Using dynamic strains from a nearby borehole strainmeter,we show that the ratio of shear to volumetric strain coupling is~0.41 which implies a mean Skempton's coefficient of A=0.24 over the plausible range of the undrained Poisson's ratio.Since these observations are made at relatively low confining pressure and differential stress,we suggest that the hydraulically conductive fracture network is a primary control on the coupling between pore pressure diffusion and elastic stresses in response to natural or anthropogenic sources.

Effect of magnetic field on poroelastic bone model for internal remodeling

This paper studies the effects of the magnetic field and the porosity on a poroelastic bone model for internal remodeling. The solution of the internal bone remodeling process induced by a magnetic field is presented. The bone is treated as a poroelastic material by Biot＇s formulation. Based on the theory of small strain adaptive elasticity, a theoretical approach for the internal remodeling is proposed. The components of the stresses, the displacements, and the rate of internal remodeling are obtained in analytical forms, and the numerical results are represented graphically. The results indicate that the effects of the magnetic field and the porosity on the rate of internal remodeling in bone are very pronounced.

Propagation of plane P-waves at interface between elastic solid and unsaturated poroelastic medium

A linear viscoporoelastic model is developed to describe the problem of reflection and transmission of an obliquely incident plane P-wave at the interface between an elastic solid and an unsaturated poroelastic medium, in which the solid matrix is filled with two weakly coupled fluids （liquid and gas）. The expressions for the amplitude reflection coefficients and the amplitude transmission coefficients are derived by using the potential method. The present derivation is subsequently applied to study the energy conversions among the incident, reflected, and transmitted wave modes. It is found that the reflection and transmission coefficients in the forms of amplitude ratios and energy ratios are functions of the incident angle, the liquid saturation, the frequency of the incident wave, and the elastic constants of the upper and lower media. Numerical results are presented graphically. The effects of the incident angle, the frequency, and the liquid saturation on the amplitude and the energy reflection and transmission coefficients are discussed. It is verified that in the transmission process, there is no energy dissipation at the interface.

Normal compression wave scattering by a permeable crack in a fluid-saturated poroelastic solid

A mathematical formulation is presented for the dynamic stress intensity factor (mode I) of a finite permeable crack subjected to a time-harmonic propagating longitudinal wave in an infinite poroelastic solid. In particular, the effect of the wave-induced fluid flow due to the presence of a liquid-saturated crack on the dynamic stress intensity factor is analyzed. Fourier sine and cosine integral transforms in conjunction with Helmholtz potential theory are used to formulate the mixed boundary-value problem as dual integral equations in the frequency domain. The dual integral equations are reduced to a Fredholm integral equation of the second kind. It is found that the stress intensity factor monotonically decreases with increasing frequency, decreasing the fastest when the crack width and the slow wave wavelength are of the same order. The characteristic frequency at which the stress intensity factor decays the fastest shifts to higher frequency values when the crack width decreases.

Mathematical osteon model for examining poroelastic behaviors

An extended and reasonable stress boundary condition at an osteon exte- rior wall is presented to solve the model proposed by Remond and Naili. The obtained pressure and fluid velocity solutions are used to investigate the osteonal poroelastic behaviors. The following results are obtained. （i） Both the fluid pressure and the velocity amplitudes are proportional to the strain amplitude and the loading frequency. （ii） In the physiological loading state, the key role governing the poroelastic behaviors of the osteon is the strain rate. （iii） At the osteon scale, the pressure is strongly affected by the permeability variations, whereas the fluid velocity is not.

Poroelastic finite-difference modeling for ultrasonic waves in digital porous cores

Scattering attenuation in short wavelengths has long been interesting to geophysicists. Ultrasonic coda waves, observed as the tail portion of ultrasonic wavetrains in laboratory ultrasonic measurements, are important for such studies where ultrasonic waves interact with smallscale random heterogeneities on a scale of micrometers, but often ignored as noises because of the contamination of boundary reflections from the side ends of a sample core. Numerical simulations with accurate absorbing boundary can provide insight into the effect of boundary reflections on coda waves in laboratory experiments. The simulation of wave propagation in digital and heterogeneous porous cores really challenges numerical techniques by digital image of poroelastic properties, numerical dispersion at high frequency and strong heterogeneity, and accurate absorbing boundary schemes at grazing incidence. To overcome these difficulties, we present a staggered-grid high-order finite-difference （FD） method of Biot＇s poroelastic equations, with an arbitrary even-order （2L） accuracy to simulate ultrasonic wave propagation in digital porous cores with strong heterogeneity. An unsplit convolutional perfectly matched layer （CPML） absorbing boundary, which improves conventional PML methods at grazing incidence with less memory and better computational efficiency, is employed in the simulation to investigate the influence of boundary reflections on ultra- sonic coda waves. Numerical experiments with saturated poroelastic media demonstrate that the 2L FD scheme with the CPML for ultrasonic wave propagation significantly improves stability conditions at strong heterogeneity and absorbing performance at grazing incidence. The boundary reflections from the artificial boundary surrounding the digital core decay fast with the increase of CPML thick- nesses, almost disappearing at the CPML thickness of 15 grids. Comparisons of the resulting ultrasonic coda Qsc values between the numerical and experimental ultrasonic S waveforms for a cylindrical rock sample demonstrate that the boundary reflection may contribute around one-third of the ultrasonic coda attenuation observed in laboratory experiments.