Governing equations of poroelastodynamics in time and frequency domain are derived. The continuity equation complements the momentum balance equations. After reduction for spherical symmetry (geometry and loading), the governing equations in frequency domain are solved by introducing wave potentials. The wave propagation velocities are obtained as the real parts of the characteristic equation of the coupled ODE system. Time domain solution for Dirac type boundary pressure is obtained through numerical inversion of transformed solutions. The results are compared to the solution in classical elasticity theory found in the literature.

1. Introduction

In wave mechanics, cavity expansion problems (circular and spherical) have received considerable attention since early 1930s [1]. These problems generally are amenable to exact solution, and the analytic solutions cast new light to the nature of wave propagation in solid media. The problem has practical applications in geophysics, seismology, and tunnel and mining engineering, as earthquake sources, underground detonation, and seismic probing. The exact solutions of such simple problems serve to understand more complex wave motions; these solutions can also be used as benchmark problems to assess the accuracy of numerical methods (FEM, BEM, FDM, etc.).

The literature on cavity problems in elastodynamics abounds [26], while corresponding work in poroelasticity has only recently been emerging [7]. An analytical solution for a dynamically loaded poroelastic column [8, 9] is also available.

This work concerns frequency domain analytical solution of dynamic pressurization of a spherical cavity in an infinite poroelastic medium. Time domain solution for Dirac type boundary pressure (pressure suddenly applied and removed) is also investigated using FFT (fast Fourier transform) technique. The developed solutions have parallel applications as the classical elastic cavity problem and clearly show the existence of a second pressure wave, the so-called slow wave.

In the sequel, a short review of Biot’s poroelasticity theory is given. The governing equations are reduced for spherical symmetry. The analytical solution in frequency domain is obtained for a permeable boundary with arbitrary but spherically symmetric pressure inside the cavity. Time domain solution is obtained via a numerical inversion using FFT.

1.1. Biot’s Theory of Poroelasticity

Unlike the classical elasticity, Biot’s theory of poroelasticity is a coupled defromation flow theory of a porous solid matrix with interstitial fluid. Biot introduced his linear quasistatic theory in 1941 [10] and later extended it to cover the dynamic range [11, 12]. An extensive review of quasi-static poroelasticity can be found in [13, 14].

The constitutive equations of linear-isotropic poroelasticity are given as (summation convention applies)where are the strains in the solid , and are components of solid displacement vector and total stress tensor, is the fluid pressure, is the variation of fluid volume per unit reference volume, and is the Kronecker delta. Here, tensile and are positive, while pore pressure is positive when it is compressive. Four elastic constants,: drained Lame’s modulus, as defined in classical elasticity (dimension F/L2);: drained shear modulus, as defined in classical elasticity (dimension F/L2);: Biot’s effective stress coefficient (dimensionless);: Biot deformation modulus, corresponds to the reciprocal of constrained storage coefficient in hydrogeology (dimension F/L2),characterize the poroelastic material.

The balance of linear momentum for the solid and fluid parts gives, respectively, Here, are relative fluid displacements, and are the components of the fluid flux vector, the three densities , , and are the bulk density, fluid density, and added fluid mass density, respectively. is the porosity of the porous matrix. is the coefficient of permeability in units L3T/M. , are the body forces per unit volume of solid and fluid, respectively. The third equation is sometimes unconventionally referred to as generalized Darcy’s law. The last balance law is the continuity equation for the fluid part where is fluid generation rate per unit volume.

In the absence of body forces and fluid generation, Fourier transforms of (1a), (1b), (2), (3), and (4) with respect to time give, when combined, the 3D governing equations in frequency domain It has been assumed that the medium is initially at rest. Here, , represent the Fourier transformed variables, that is, and is a coefficient defined as In dyadic form (5) reads where is the solid displacement vector. Analytical solutions of the 3D governing equations are scarse if not impossible. There are the point load solutions for infinite poroelastic medium, the fundamental solutions, for both quasi-static and dynamic poroelasticity [1518] and poroelastic Lamb’s problem [19] available in the literature. Also a good compendium of anlytical solutions in the quasi-static range can be found in [14].

2. Governing Equations for Spherical Symmetry

Here, a spherical cavity of radius “” in an infinite PE medium (Figure 1) will be considered. By substituting the gradient and divergence in (8), governing equations in spherical coordinates can be obtained.

Let , , be the spherical components of solid displacements, because the cavity is spherically symmetric, and we assume a time varying but spherically symmetric pressure inside; waves emanating from such a source will have spherical symmetry, that is, Hence, the governing equations of 3D poroelasticity in FTS, in this case, reduce to the following: Here, denotes the only nonzero displacement for simplicity, that is, radial component .

The stress components are related to the radial displacement as Other .

3. General Solution of the Governing Equations

If one introduces a displacement potential such that then (10) and (11) become where Now, we assume solutions of the form where is a frequency dependent parameter. Substituting (14) in (12), one obtains the eigenvalue problem: Hence, are the roots of the characteristic equation With the following abbreviations: the roots of the characteristic equation are and the solution becomes It is clear that the terms represent waves propagating inwards. In an infinite medium, no waves propagate from infinity (Sommerfeld’s radiation condition [1]); discarding these terms, (19) becomes Moreover, the amplitudes , are related through (15) as introducing (21) in (20) Furthermore, since and Finally, the general solution for radial displacement and pore pressure becomes The integration constants and are to be evaluated from boundary conditions at .

3.1. Boundary Conditions

Among various possible combinations of traction, displacement, pore pressure, and fluid flux, only the permeable boundary condition will be considered here, that is, The stress boundary condition corresponds to a loading of in time domain, where is the strength of loading and is the Dirac delta function (or distribution). When (25) are introduced in (24) one obtains where

3.2. Analytical Solution in Frequency Domain

Thus, the complete analytical solutions in frequency domain for Dirac type pulse are Stresses can be found by inserting (28) in (11). The two terms in each of (28) represent two -phases with different propagation characteristics. The reciprocal of real part of and corresponds to fast and slow (second) longitudinal wave velocities. The second (slow) -wave is highly dispersive and highly damped.

3.3. Solution in the Classical Theory of Elasticity

The elastodynamic solution in frequency domain for suddenly pressurized spherical cavity can be derived following a similar outline and is found in [1, 20]; here we restate the form of the solution, which complies with the Fourier transform definition in (6): where is the displacement, , are the - and -wave velocities in an idealized elastic medium, and other parameters are as defined previously.

4. Results in Frequency and Time Domain

In all of the computations material constants (Table 1) for Berea Sandstone are used [21].

The material parameters for corresponding elastic medium are as in Table 2.

The real and complex parts of displacement and pore pressure solutions in frequency domain are plotted in Figures 2 and 3, respectively, note that the elastic theory doesn’t involve pore pressure.

The following (Figure 4) are the dispersion curves for the fast and slow wave. As seen in the figure, the wave front of fast -wave travels with a speed equal to 3137.2 m/s where as the slow wave front velocity is 1037 m/s. While, both waves display dispersion, the slow wave is rather more dispersive than the fast wave.

4.1. Time Domain Solution Using FFT Technique

The frequency domain analytical solutions (28) are extremely complex (if not impossible) to be inverse transformed (by analytical means) into time domain. Thus, time domain solutions are obtained using numerically inverting (28) by FFT (fast Fourier transform) technique.

Since Dirac impulse is highly transient, a high frequency resolution is required, thus the following parameters are chosen for FFT: seconds (analysis interval), (exponent in FFT),

which give a frequency resolution of 0.628 rad/s and a Nyquist frequency of 41177 rad/sec. The results are plotted in Figure 5.

As seen from Figure 5, the fast primary wave in poroelastic medium travels faster (arrives earlier) than -wave in elastic medium, with the amplitude of elastic wave being slightly larger. The second slow wave in poroelastic wave cannot have been captured, probably because of its small amplitude and high damping and dispersion.

5. Conclusions

Analytical solution of waves propagating from the surface of a suddenly pressurized poroelastic spherical cavity with permeable boundary is obtained in frequency domain. The problem is one of the classical problems in wave propagation in elastic solids [1, 20]. Among many other possibilities of boundary conditions, this paper, which employs permeable boundary conditions, is one of the first attempts to the solution of a variety of other problems related to spherical wave propagation from cavities in poroelastic media.

Time domain solution is obtained with inverse FFT. Since, Dirac type pulse is a highly transient excitation, a high frequency resolution is necessary. The results clearly show that -waves in poroelastic media propagate faster, with slightly lower amplitude, as compared to elastic media with similar material properties. The slow wave could not have been visualized in the graphics, because of small amplitude, high damping, and dispersion characteristics of this phase.