## Abstract

Here processes of wave propagation in a two-component Biot’s medium are considered, which are generated by arbitrary forces actions. By using Fourier transformation of generalized functions, a fundamental solution, Green tensor, of motion equations of this medium has been constructed in a non-stationary case and in the case of stationary harmonic oscillation. These tensors describe the processes of wave propagation (in spaces of dimensions 1, 2, 3) under an action of power sources concentrated at coordinates origin, which are described by a singular delta-function. Based on them, generalized solutions of these equations are constructed under the action of various sources of periodic and non-stationary perturbations, which are described by both regular and singular generalized functions. For regular acting forces, integral representations of solutions are given that can be used to calculate the stress-strain state of a porous water-saturated medium.

### Keywords

- Biot’s medium
- solid and liquid components
- Green tensor
- Fourier transformation
- regularization

## 1. Introduction

Various mathematical models of deformable solid mechanics are used to study the seismic processes of earth’s crust. The processes of wave propagation are most studied in elastic media. But these models do not take into account many real properties of an ambient array. These are, for example, the presence of groundwater, which affects the magnitude and distribution of stresses. Models, which take into account the water saturation of earth’s crust structures, presence of gas bubbles, etc., are multicomponent media. A variety of multicomponent media, complexity of processes associated with their deformation, lead to a large difference in methods of analysis and modeling used in solving such problems.

Porous medium saturated with liquid or gas, from the point of view of continuum mechanics, is essentially a two-phase continuous medium, one phase of which is particles of liquid (gas) and other solid particles are its elastic skeleton. There are various mathematical models of such media, developed by various authors. The most famous of them are the models of Biot, Nikolaevsky, and Horoshun [1, 2, 3, 4, 5]. However, the class of solved tasks to them is very limited and mainly associated with the construction and study of particular solutions of these equations based on methods of full and partial separation of variables and theory of special functions in the works of Rakhmatullin, Saatov, Filippov, Artykov [6, 7], Erzhanov, Ataliev, Alexeyeva, Shershnev [8, 9], etc. In this regard, it is important to develop effective methods of solution of boundary value problems for such media with the use of modern mathematical methods.

Periodic on time processes are very widespread in practice. By this cause, here we consider also processes of wave propagation in Biot’s medium, posed by the periodic forces of different types. Based on the Fourier transformation of generalized functions, we constructed fundamental solutions of oscillation equations of Biot’s medium. It is Green tensor, which describes the process of propagation of harmonic waves at a fixed frequency in the space–time of dimension N = 1, 2, 3, under the action concentrated at the coordinates origin. By using this tensor, we construct generalized solutions of these equations for arbitrary sources of periodic disturbances, which can be described as both regular and singular distributions. They can be used to calculate the stress-strain state of a porous water-saturated medium by seismic wave propagation.

## 2. The parameters and motion equations of a two-component Biot’s medium

The equations of motion of a homogeneous isotropic two-component Biot’s medium are described by the following system of second-order hyperbolic equations [1, 2, 3]:

Here *N* is the dimension of the space. At a plane deformation *N* = 2, the total spatial deformation corresponds to *N* = 3, at *N* = 1 the equations describe the dynamics of a porous liquid-saturated rod.

We denote *N*).

Constants

where *m* is a porosity of the medium. The constant of attached density *Q, R* characterize an interaction of a skeleton with a liquid on the basis of.

### 2.1 Biot’s law for stresses

Here

Further we use the next notations for partial derivatives:

There are three sound speeds in this medium:

where the next constants were introduced as:

The first two speeds *dilatational waves*. The second slower dilatation wave is called *repackaging wave*. A third velocity *shear waves* and at

We introduce also two velocities of propagation of dilatational waves in corresponding elastic body and in an ideal compressible fluid:

## 3. Problems of periodic oscillations of Biot’s medium

Construction of motion equation solutions by periodic oscillations is very important for practice since existing power sources of disturbances are often periodic in time and therefore can be decomposed into a finite or infinite Fourier series in the form:

where periods of oscillation of each harmonic *T*. Therefore, it is enough to consider the case of stationary oscillations, when the acting forces are periodic on time with an oscillation frequency

The solution of Eq. (1) can be represented in the similar form:

where complex amplitudes of displacements

which give us the solution of problem for forces (4).

We get equations for complex amplitudes by stationary oscillations, substituting (6) into the system (1):

To construct the solutions of this system for different forces, we define Green tensor of it.

## 4. Green tensor of Biot’s equations by stationary oscillations

Let us construct

Here *x* = 0. The upper index of this tensor (…^{[k]}) fixes the current concentrated force and its direction. The lower index corresponds to component of movement of a skeleton and a fluid, respectively,

Their complex amplitudes

Since fundamental solutions are not unique, we’ll construct such, which tend to zero at infinity:

and satisfy the radiation condition of type of Sommerfeld radiation conditions [10].

Matrix of such fundamental equations is named *Green tensor* of Eq. (8).

## 5. Fourier transform of fundamental solutions

To construct

where

Let us apply Fourier transformation to Eq. (10) and use property of Fourier transform of derivatives [10]:

Then we get the system of 2 *N* linear algebraic equations for Fourier components of this tensor:

By using gradient divergence method, this system has been solved by us. For this the next basic function were introduced:

and the next theorem has been proved [11, 12].

**Theorem 1.***Components of Fourier transform of fundamental solutions have the form:*

*where the next constants have been introduced as:*

This form is very convenient for constructing originals of Green tensor.

## 6. Stationary Green tensor construction: radiation conditions

In this case let us construct the originals of the basic function but only over

which, in accordance with its definition (14), satisfies the equation:

Using property (12) for derivatives from here, we get Helmholtz equation for fundamental solution (accurate within a factor

Fundamental solutions of Helmholtz equation, which satisfy Sommerfeld conditions of radiation:

are well known [10]. They are unique. Using them, we obtain:

where

These functions (subject to factor *x* = 0 to infinity and decay at infinity.

The last property is true only for *N* = 2,3. In the case *N* = 1, all fundamental solutions of Eq. (16):

do not decay at infinity.

From Theorem 1, the next theorem follows.

**Theorem 2.***The components of Green tensor of Biot’s equations at stationary oscillations with frequency ω, which satisfy the radiation conditions, have the form:*

where

*Proof*. By using originals of basic functions, property (12) of derivatives, we can obtain from formulas for

Then formulas of Theorem 2 follow from formulas of Theorem 1.

## 7. Generalized solutions by arbitrary periodic forces

Under the action of arbitrary mass forces with frequency

Note that mass forces may be different from the space of generalized vector function, singular and regular. Since Green tensor is singular and contains delta-functions, this convolution is calculated on the rule of convolution in generalized function space. If a support of acting forces are bounded (contained in a ball of finite radius), then all convolutions exist. If supports are not bounded, then the existence conditions of convolutions in formula (17) requires some limitations on behavior of forces at infinity which depends on a type of mass forces and space dimension.

The obtained solutions allow us to study the dynamics of porous water- and gas-saturated media at the action of periodic sources of disturbances of a sufficiently arbitrary form. In particular, they are applicable in the case of actions of certain forces on surfaces, for example, cracks, in porous media that can be simulated by simple and double layers on the crack surface.

There is another interesting feature of the Green tensor of the Biot’s equations, which contains, as one of the terms, the delta-function that complicates the application of this tensor for solving boundary value problems based on analogues of Green formulas for elliptic systems of equations or the boundary element method. Here, when constructing the model, the viscosity of the liquid is not taken into account, which, apparently, leads to the presence of such terms, and it requires improvement of this model taking into account a viscosity.

## 8. Green tensor of Biot’s equations by non-stationary motion

To construct the non-stationary Green tensor, at first we also construct the originals of the basic functions in an initial space-time:

They are originals of the classic wave equation:

Depending on the dimension of a space, solutions of this wave equation that satisfy the radiation conditions have the following form [10]:

Here *simple layer on the sound cone*

Using regularization of the general function

and the properties of Fourier transform of generalized functions convolution:

It is easy to show that the next lemma is true.

**Lemma.***The originals of the primitives of the basic functions satisfying the radiation conditions are representable in the following form:*

*for N = 3*

for N = 2

for N = 1

Using these functions and the properties of the Fourier transform, we obtain the components of the Green tensor from the formulas of Theorem 1. We formulate the result in the next theorem.

**Theorem 3.***The components of Green tensor of motion equations of two-component Biot’s medium have the following forms:*

Here

for *N* = 1

for *N* = 2

for *N* = 3

## 9. Generalized solutions of Biot’s equations by non-stationary forces

Using the properties of Green tensor, we obtain generalized solutions of non-stationary Biot’s equations under the action of arbitrary mass forces in the Biot’s medium, which satisfy the radiation condition at infinity. They have the form of tensor functional convolution:

It’s taken according to the rules of convolution of generalized functions depending on the type of mass forces [10].

In order to get the classic solution, we must present formulas (28) in regular integral forms. For this, let us present matrix of Green tensor as sum of regular functions and singular functions, which contain delta-function:

Then also write:

Here

The convolution with singular part is equal to:

In 3D space, there are convolutions with simple layers on sound cones (see (27)). To construct their integral presentation, use this rule:

Here the internal integral is taken over sphere with center in the point *x*, and its radius is equal to

If components of acting forces

Substituting the formulas of Theorem 4 into (29), we obtain displacements and stresses of skeleton and liquid in Biot’s medium in spaces of dimension *N* = 1, 2, 3. Calculation of these convolutions by using these formulas essentially depends on the form of acting forces and gives possibility to construct regular presentation of generalized solution for wide class of acting forces, which are the classic solution of Biot’s equation.

## 10. Calculation of the stress state of Biot’s medium

Using Biot’s law (2), we can define the generalized stresses in skeleton and a pressure in a liquid:

These formulas also can be written in integral form by using the same rules. But we can apply here the next lemma, which was proved in [11].

**Lemma.***Fourier transformations of the divergences of Green tensor have the next form:*

From this lemma, we can prove easily the next theorem.

**Theorem 4.***Divergences of elastic and liquid displacement of Green tensor have the next form:*

Substituting these formulas in (30), we define the stresses in the skeleton and the pressure in the liquid of Biot’s medium.

If we paste

## 11. Conclusion

The obtained solutions give possibility to study the dynamics of porous water- and gas-saturated media and rods under actions of disturbance sources of different forms and can be used for solutions of boundary value problems in porous media by using boundary element method.

These solutions can be used for describing wave processes by explosions and earthquakes. In these cases mass forces are described by using singular generalized function, such as multipoles, simple and double layers, and others.

## Acknowledgments

This work was financially supported by the Ministry of Education and Science of the Republic of Kazakhstan (Grant AP05132272).