The Inﬂuence of Permeability on the Propagation Characteristics of the Waves in Different Saturated Soils

: The permeability of saturated soils has great inﬂuence on the velocities and attenuation characteristics of fast compressional wave P1, low compressional wave P2, and shear wave S in saturated soils, respectively. In three different cases, namely zero, ﬁnite, and inﬁnite permeability, the wave equations and theoretical velocities of P1, P2, and S wave in saturated soils are given based on the u-w- p equation, respectively. According to the solutions of the wave equations, the real velocities and attenuation coefﬁcients of three waves are redeﬁned, respectively. In different saturated soils, the inﬂuences of the permeability and the loading frequency on the wave velocities and attenuation are discussed, respectively. Moreover, the suitable application scope of the u - p equation is discussed based on different permeabilities and loading frequencies.


Introduction
Saturated soils are a saturated porous medium consisting of a solid skeleton and pore water. Because of the different materials and the fluid-solid couplings, the propagation characteristics of waves in saturated soils are different from those of single-phase soils. For engineering involving the dynamic problems of saturated soils, such as petroleum engineering, geophysics, earthquake engineering and the seabed acoustics etc. [1,2], the analysis of the wave theory of saturated soils is employed to provide guidance and reference.
Biot first constructed the wave theory of saturated porous medium in the 1950s, namely the Biot theory [3,4]. However, the Biot theory is uncomprehensive due to the abstract concepts of several parameters in the controlling equation, and the parameters are difficult to assign; therefore, many scholars further developed the Biot theory [5][6][7][8]. Zienkiewicz proposed the dynamic u-w-p equation, which has a clear physical concept by ignoring the added mass density and introducing the dynamic of permeability to replace fluid viscosity. Moreover, the u-p formulation [9,10] is the simplified form of the u-w-p equation and is widely applied due to its simple equation and high computational efficiency. However, because of the simplified assumption of ignoring the acceleration of the fluid phase relative to the solid phase, the u-p equation is limited in solving dynamic problems under low or medium frequency loadings.
Biot [3] successfully predicted that there are three body waves (two compressional waves, the P1 and P2 waves, and one shear wave, the S wave) in saturated soils, and acquired the velocity of each wave. He stated that the interaction between solid and fluid phases in saturated soils consists of inertia, mechanical, and viscous couplings [3,4,11,12]. The inertia and mechanical couplings are the main reasons for the generation of the three body waves, and the propagation of waves in saturated soils is attenuated and dispersed due to the viscous coupling. Moreover, the existences of the three waves in saturated soils were confirmed by many experiments. Plona [13] used a model of glass balls submerged in

Equations of Saturated Soils
Saturated soils are regarded as a two-phase medium consisting of a solid skeleton and the pore fluid in the voids of the solid skeleton. The pore fluid and the interaction between the two phases make the mechanical properties of saturated soils very different from one-phase soils. Therefore, two governing equations describing the mechanical properties of saturated soils are determined.

u-w-p Formulation
The u-w-p formulation is a complete equation of saturated porous media. It consists of the basic governing equations, shown as follows: (1) Mass conservation equation of the pore fluid: (2) Dynamic equilibrium equations: (3) Effective stress principle: where the definitions of the symbols used in Equations (1)-(4) are listed in Table 1.
Equations (1)-(4) consist of the basic governing equations of saturated porous media.
where, K T , K S and K f are the bulk modulus of soil skeleton, soil particles and pore fluid respectively. For the material of soil, soil particles are nearly incompressible relative to the soil skeleton, thus they have the relations of K S >> K D , α ≈ 1.

Q b
Compressibility coefficient of fluid Equations (1) and (4) are substituted into Equations (2) and (3) and the u-w-p formulation is expressed by displacement and the pore pressure variables can be written as follows:

u-p Formulation
Zienkiewicz [10] assumed that the relative acceleration of the fluid phase can be ignored as the relative motion between the fluid and solid phases are slow. By ignoring the terms of the fluid acceleration in the u-w-p formulation, the simplified u-p formulation can be written as: where the definitions of the symbols in Equation (7) are listed in Table 1. The simplified u-p formulation has simpler equations and a higher computational efficiency relative to the u-w-p formulation. However, the assumption [10] limits the scope of application of the u-p formulation, which will be discussed in Section 6.

Wave Equations and Theoretical Wave Velocities
The pore fluid makes the propagation properties of the waves in saturated soils different from those of common soils, especially under dynamic loadings. The propagation procedures of the waves in saturated soils are closely related to the viscous coupling, so the seepage procedure of the pore fluid affects the velocities and the attenuation of the waves in saturated soils. Therefore, analysis of the propagation theory of saturated soils is necessary. In this section, for three different viscous couplings (finite, infinite, and zero permeability), the u-w-p formulation is transformed into three wave equations for the P1, P2 and S waves, and the corresponding theoretical velocities are obtained.

Finite Permeability
Using Fourier transform, the form of Equation (5) in the frequency domain can be written as: where ρ m is introduced to make the formulation similar to the wave equations and has the relation of: ρ m = ρ f n − 1 k i ω . Displacement u and pore pressure p are decomposed by Helmholtz decomposition: where ϕ s and ϕ f are the irrotational potential function of the solid and fluid phases, respectively. ψ s is the rotational potential function of the solid phase. Substituting Equations (10) and (11) into Equations (8) and (9), and using the relation of ∇ · ψ s e = 0, the following equations can be obtained: ϕ s , ϕ f and ψ s in Equation (12) should satisfy the following relations: Equations (13)- (15) can be expressed as three scalar equations and is written into the matrix forms shown below: where,

Compressional Wave
The characteristic equation of the Equation (16) is: The dispersion equation can be obtained by unfolding Equation (19): The solution of Equation (20) is: where, Equation (21) is two complex solutions of Equation (20), respectively, where wave numbers l P1,2 and wave velocities C P1,2 are the theoretical wave numbers and velocities of two compressional waves in saturated soils, respectively. P1 and P2 represent the fast and low compressional waves in saturated soils, respectively.
In addition, the simplifier form can be written as below, by substituting Equation (14) into Equation (13) The potential functions of the solid and fluid displacements are, respectively, expressed by the potential functions of the P1 and P2 waves [34] where, Substituting Equation (24) into Equation (23) and using Equation (20), we can obtain the following equations: Equation (26) is obtained as follows using the Fourier inverse transform: In Equation (27), the equations of the potential functions of the P1 and P2 waves have the forms of wave equations.

Shear Wave
Similarly, Equation (17) of the rotational potential function ψ s is transformed into the time domain: where C S is the theoretical velocity of the shear wave (S wave). Equation (28) is the wave equation of the S wave and its characteristic equation is: The solution of Equation (29) is expressed as: where l S is the theoretical wave number of the S wave.

Infinite Permeability
For the case of infinite permeability, the pore fluid flow through the solid porosity without any obstacle and there is no viscous coupling between the pore fluid and the solid skeleton in saturated soils. The properties of saturated soils are similar to the single-phase soils. In Equations (8) and (9), the permeability is set to k = ∞ and the relation ρ m ≈ ρ f /n is satisfied. Substituting the relation of k = ∞, Equations (8) and (9) is transformed into: Using the similar derivational process of the case of the finite permeability, the wave numbers and velocities of P1, P2, and S waves in infinite permeability cases, respectively, are: where, From Equation (33), we can see that there are two compressional waves and one shear wave in saturated soils in the case of infinite permeability. The forms of the wave equations of P1, P2 and S waves are same with Equations (27) and (28).

Zero Permeability
For the case of zero permeability, there is no relative motion between the pore fluid and the solid skeleton. The strong viscous coupling makes the P2 wave attenuate to zero rapidly. Because w is equal to zero, the fluid relative velocity . w and acceleration ..
w are equal to zero, and the terms of permeability in Equations (5) and (6) are ignored. Equations (5) and (6) is simplified as: Referring to the method in the literature [35], Equation (36) is substituted into Equation (35). Two uncoupled equations about only the displacements are shown below: In two dimensional cylindrical coordinates, the displacement potential functions, respectively, are expressed by the radial and tangential components: Equation (38) is substituted into Equation (37), and the wave equations of the displacement potential functions are obtained: where the velocities of P1 and S waves, respectively, are: Equation (39) and Equation (40), respectively, are the wave equations of P1 and S waves under the zero-permeability condition.

Solutions of Wave Equations
In saturated soils with different permeabilities, all of the equations to describe the dynamic behaviors of the three body waves have the forms of wave equations. In this section, the solutions of the wave equations are determined in the forms of plane, cylindrical, and spherical waves.

Plane Wave
For the one-dimensional plane wave, the solutions of the potential functions of three waves are expressed as below [36]: where A 1 , A 2 and A 3 are the amplitudes of the P1, P2, and S waves, respectively. ω is circular frequency, x is the direction of wave propagation, and t represents time.

Cylindrical and Spherical Waves
For cylindrical and spherical waves, the solutions of the potential functions of the three waves are approximately expressed in Equations (43) and (44) [37], respectively where f is random function of outgoing waves. P1 and P2 waves have the same wave functions, f. r is the radial direction of the wave propagation. However, different from the one-dimensional plane wave, the cylindrical and spherical waves have geometric attenuation during the diffusion process.

Discussions of Wave Velocities, Dispersion and Attenuation Characteristics
In the cases of infinite and zero permeability, the theoretical velocities in Equations (33) and (41) of the three body waves are all real numbers and independent of circular frequency. This indicates that the three waves have constant velocities and the propagation procedures are not attenuated and dispersed. For zero permeability, the density of the mixture medium is used to compute the velocity of the shear wave. Both of the solid and fluid phases undertake the shear force. For infinite permeability, only the density of the solid phase is used to compute the velocity of the shear wave. So, the saturated soil can be treated as a one-phase soil and the shear force is only undertaken by the solid skeleton.
For the finite permeability, the relation of (21) and (30). The theoretical velocities of three waves are complex numbers, which are functions of the permeability and the circular frequency. Therefore, the propagation properties of the three waves in saturated soils are analyzed from the three aspects of dispersion, velocity, and attenuation characteristics, respectively.

Dispersion and Real Wave Velocity
Dispersion is a procedure where waves are gradually separated with different frequencies. Thus, the wave velocity, initial waveform, and vibration mode are changed during the propagation process.
In the case of the finite permeability, theoretical velocities of the three waves are uniformly expressed by the complex number C. Thus, the velocity can be written in the following form: where j is the root of the dispersion equation. J = 1, 2, 3. C r and C i represent the real and imaginary parts of the root, respectively. In Equation (36), real and imaginary parts of the velocity are expressed as a function of the circle frequency, which reflects the dispersion of the waves. The one-dimensional plane wave is employed to simplify the analysis progress. Equation (42) is uniformly written as the form shown below: where, Equation (47) is substituted into Equation (46), and the following equation is obtained: where, For the problems of a one-dimensional plane wave, the potential functions of three waves in a saturated porous medium can be expressed in a unified form of Equation (48). In Equation (48), the term A j exp ω C ij (ω) x describes the attenuation of the wave propagation, and the term exp iω 1 C rj (ω) x + t reflects the propagation of the waves in a simple harmonic mode. Therefore, the real part of Equation (50) is the real wave velocity for the finite permeability case in a saturated porous medium.

Attenuation Coefficient
In the case of finite permeability, the velocity has an imaginary part, which illustrates that waves are attenuated in the propagation procedure.
Similarly, the potential function in Equation (48) can be expressed as an equivalent function of wave number l.
where the l r and l i are the real and imaginary parts of the wave numbers, respectively. The wave number is l j (ω) = l rj (ω) + l ij (ω)i. The term of exp[(l r x + ωt)i] represents the mode of wave in the propagation procedure, and the amplitude of A exp(−l i x) reflects that the amplitude of the wave attenuates in the exponential form.
The attenuation of the velocity can be described by the attenuation coefficient β, which is defined that the ratio of the reduction of the amplitude after a wavelength and the amplitude in a wavelength. The detailed formulation is shown as Equation (52) where L is the wavelength in a period. The amplitude attenuation coefficient is a value between zero and one (0 < β < 1). The calculation of wave number l refers to Equations (21) and (30).

Comparison of Velocities and Frequency Dispersion in Different Saturated Soils
Four saturated porous medium, saturated gravel soil, sand soil, clay silt, and silt, were selected to analyze propagation characteristics. The mechanical parameters are listed in detail in Table 2.   Figure 1 shows the relation of the frequency, permeability, and velocities of the P1, P2, and S waves in the case of finite permeability. From Figure 1, as the permeability k < 1 × 10 −6 m/s, the velocities of the three waves do not change with the increase in permeability and loading frequency. The values always maintain the minimum value, which is close to the minimum velocity under the condition of zero permeability. As the permeability k > 1 × 10 −6 m/s, with the increase in permeability, the velocities of the three waves increase gradually to their maximum velocities, which are close to the velocities in the condition of infinite permeability. As the permeability is in the range of k < 1 × 10 −6 m/s, the velocities of the three waves increase with the increase in the loading frequency. This indicates that the three waves are dispersed in the propagation procedure. Moreover, the permeability is proportional to the velocities, and the influence of loading frequency on the velocities is less than the permeability. to replace the actual wave velocities in this condition. For the saturated soils with a permeability, such as sand and gravel, the velocities fall in between the maximum minimum velocities. The velocities of the three waves should be calculated by the velocity equations in this condition.  Figure 2 shows the relations of permeability, loading frequency, and the attenu coefficients of the P1, P2, and S waves in the case of finite permeability. For P1 waves, as the permeability reaches k < 1 × 10 −4 m/s the attenuation coefficients alwa In Figure 1, by contrasting the relations of permeability, loading frequency, and the velocities in the four saturated soils, we can see that different saturated soils do not change the shapes of the velocity curves, but vary the amplitudes of the curves. The harder the soil, the faster the propagation wave is. The permeability coefficients of the grey planes in Figure 1 represent the actual permeability coefficients of the four saturated soils. They are k = 3.24 × 10 −3 , 4.6 × 10 −4 , 7.9 × 10 −6 and 6 × 10 −10 m/s, respectively. For the saturated soils with low permeability, such as clay silt and silt, the permeability coefficients are less than 1 × 10 −6 m/s. The corresponding velocities are close to the minimum velocities and do not change with the loading frequency. Therefore, the minimum wave velocities can be used to replace the actual wave velocities in this condition. For the saturated soils with a high permeability, such as sand and gravel, the velocities fall in between the maximum and minimum velocities. The velocities of the three waves should be calculated by the actual velocity equations in this condition. Figure 2 shows the relations of permeability, loading frequency, and the attenuation coefficients of the P1, P2, and S waves in the case of finite permeability. For P1 and S waves, as the permeability reaches k < 1 × 10 −4 m/s the attenuation coefficients always remain at zero, which illustrates that the two waves do not attenuate in the propagation procedure. With the increase in the permeability, a wave crest appears in the range of k > 1 × 10 −4 − 1 m/s and the maximum values of the attenuation coefficients appear as permeability is about 1 × 10 −2 m/s for the P1 and S waves. This indicates that two waves show slight attenuation. For the P2 wave, as permeability k > 1 × 10 −4 m/s, the attenuation coefficients change from 1.0 to 0.0 with the increase of the permeability. The P2 wave gradually attenuates to zero. Contrasted with the loading frequency, the permeability has more influence on the attenuation coefficient. Comparing the relations of permeability, loading frequency, and the attenuation efficient in the four saturated soils, we can find that different saturated soils do not cha the shapes of the attenuation coefficient curves, but do change the amplitudes of th curves. Especially when the saturated soils changed from the gravel to the silt, the am tudes of the attenuation coefficients of the S wave obviously shows an increased tende The S wave attenuates more significantly in soft soils.  Comparing the relations of permeability, loading frequency, and the attenuation coefficient in the four saturated soils, we can find that different saturated soils do not change the shapes of the attenuation coefficient curves, but do change the amplitudes of these curves. Especially when the saturated soils changed from the gravel to the silt, the amplitudes of the attenuation coefficients of the S wave obviously shows an increased tendency. The S wave attenuates more significantly in soft soils.

Analysis of the u-p Equation
Based on the u-p formulation, wave velocities and attenuated coefficients are deduced using the same procedure as the u-w-p formulation. Moreover, the scope of application of the u-p formulation is discussed. Equation (7) is transformed into the frequency domain: In the same way, using the Helmholtz decomposition, displacement u and pore pressure p, respectively, are substituted into Equations (53) and (54). We can obtain the equations expressed by potential functions: Equations (55)-(57) is written in a matrix form, shown as follows: where, The characteristic equation of Equations (55)-(57) is: The solutions of Equation (59) are shown as follows: where, From Equation (62), we can see that the theoretical velocities of the P1 and P2 waves have the imaginary parts, but the shear wave has a constant velocity. Thus, the two compressional waves are attenuated and the shear wave is not attenuated during the propagation procedure.
In addition, Equation (60) is expanded as: The Equation (56) is substituted into Equation (55) and the following equation can be obtained Moreover, the eigenvectors corresponding to the eigenvalue in Equation (62) are expressed as follows: where, Using Equation (67), the potential functions of the solid and fluid displacements can be expressed by the potential function of the P1 and P2 waves. The detailed relations are shown below: The following two compressional wave equations can be obtained by substituting Equation (68) into Equation (66) and using Equation (65).
From Equation (69), the potential functions of displacements of P1 and P2 waves have the forms of wave equations. The wave equation of the S wave is the same as Equation (28). The velocities of the three waves refer to Equations (62) and (63).
To discuss the application of the u-p formulation, the product of k·ω of the permeability and the frequency is used as a variation. To explain the accuracy of the simplified u-p formulation, Figure 3 shows the differences in the velocities of P1, P2, and S waves, obtained from the simplified u-p formulation and the completed u-w-p formulation, respectively. In addition, comparisons of the velocities are conducted in saturated gravel soil, sand soil, clay silt and silt, respectively. In Figure 3, blue shade fields represent the actual range of k·ω. By discuss the relative errors of the velocities of the P1, P2 and S wave in the actual range, the scope of the application of the u-p formulation can be gotten preliminarily. soil, sand soil, clay silt and silt, respectively. In Figure 3, blue shade fields represent the actual range of k·ω. By discuss the relative errors of the velocities of the P1, P2 and S wave in the actual range, the scope of the application of the u-p formulation can be gotten preliminarily.
In clay silt and silt, value of k·ω is less than 10 −3 and 10 −6 , respectively. In this scope, three velocities obtained from u-p formulation are same with them from the u-w-p formulation. The relative errors of the velocities between the two formulations are zeros. In gravel and sand soils, the values of k·ω are 10 −5 -10 2 and 10 −7 -1, respectively. In these scopes, both of the relative errors of the velocities between the two formulations are beyond the relative error limit of 3% in the two soils. Therefore, the application of the u-p formulation would be unsuitable in soils with high permeability, like sand and gravel. We mainly discuss the influence of the permeability and frequency on the wave velocities gained from the u-p formulation as the value of k·ω is between 10 −7 -10 2 . In this range, the primary range of the loading frequency is 10 −2 -10 2 Hz, so the range of the permeability is 10 −5 − 1 m/s.    Figure 4 shows the relations of the permeability, frequency, and the velocities in gravel. By comparing the relative errors of the velocities of P1, P2, and S waves, obtained from the u-p and u-w-p formulations, the variation of scope of the u-p formulation with increasing permeability and frequency are further discussed.
In Figure 4, the relative errors of the P2 wave velocities are the largest of the three waves. Therefore, taking the relative errors of the P2 wave velocities in sand soils as an example (in Figure 4), the relative errors of the three wave velocities obtained, respectively, from the u-p formulation and the u-w-p formulation are less than 3% as the permeability k is less than 10 −2 m/s and frequencyω is less than 1 Hz. In addition, the permeability k and frequencyω satisfy the relation of ω = -100k + 101. Within this field, the wave velocities are mostly closed to the real velocities, and the u-p formulation can be safely used. In clay silt and silt, value of k·ω is less than 10 −3 and 10 −6 , respectively. In this scope, three velocities obtained from u-p formulation are same with them from the u-w-p formulation. The relative errors of the velocities between the two formulations are zeros. In gravel and sand soils, the values of k·ω are 10 −5 -10 2 and 10 −7 -1, respectively. In these scopes, both of the relative errors of the velocities between the two formulations are beyond the relative error limit of 3% in the two soils. Therefore, the application of the u-p formulation would be unsuitable in soils with high permeability, like sand and gravel. We mainly discuss the influence of the permeability and frequency on the wave velocities gained from the u-p formulation as the value of k·ω is between 10 −7 -10 2 . In this range, the primary range of the loading frequency is 10 −2 -10 2 Hz, so the range of the permeability is 10 −5 − 1 m/s. Figure 4 shows the relations of the permeability, frequency, and the velocities in gravel. By comparing the relative errors of the velocities of P1, P2, and S waves, obtained from the u-p and u-w-p formulations, the variation of scope of the u-p formulation with increasing permeability and frequency are further discussed.
In Figure 4, the relative errors of the P2 wave velocities are the largest of the three waves. Therefore, taking the relative errors of the P2 wave velocities in sand soils as an example (in Figure 4), the relative errors of the three wave velocities obtained, respectively, from the u-p formulation and the u-w-p formulation are less than 3% as the permeability k is less than 10 −2 m/s and frequency ω is less than 1 Hz. In addition, the permeability k and frequency ω satisfy the relation of ω = −100k + 101. Within this field, the wave velocities are mostly closed to the real velocities, and the u-p formulation can be safely used.

Conclusions
The propagation properties of waves in saturated soil are the theoretical basis for analyzing the dynamic problems of saturated soils. The waves and their propagation properties are also different in saturated soils with different permeabilities. In this paper, the following research works are carried out for saturated soils with different permeabilities, and the corresponding conclusions are obtained: (1) Based on the u-w-p formulation of saturated porous medium, the wave equations and the corresponding velocities of P1, P2, and S waves are acquired under zero, finite, and infinite permeabilities, respectively. (2) The differences in the permeability and loading frequency of saturated soils have different influences on propagation properties, such as the velocities, dispersion and attenuation characteristics.

Conclusions
The propagation properties of waves in saturated soil are the theoretical basis for analyzing the dynamic problems of saturated soils. The waves and their propagation properties are also different in saturated soils with different permeabilities. In this paper, the following research works are carried out for saturated soils with different permeabilities, and the corresponding conclusions are obtained: (1) Based on the u-w-p formulation of saturated porous medium, the wave equations and the corresponding velocities of P1, P2, and S waves are acquired under zero, finite, and infinite permeabilities, respectively. (2) The differences in the permeability and loading frequency of saturated soils have different influences on propagation properties, such as the velocities, dispersion and attenuation characteristics.