Wave propagation characteristics in frozen saturated soil

ChengCheng Du , DongQing Li , Feng Ming , YuHang Liu , XiangYang Shi

1. State Key Laboratory of Frozen Soil Engineering, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, Gansu 730000, China

2. University of Chinese Academy of Sciences, Beijing 100049, China

1 Introduction

According to the classical wave theory, wave velocity is inextricably interwoven with physical and mechanical parameters of the medium (Li and Li,1989). The study of wave propagation in different materials has always been an important research subject.It not only has important theoretical and academic value, but also has an important guiding significance for engineering applications. In recent years, with the development of engineering construction in permafrost regions, dynamic characteristics of frozen soils has garnered increased attention for academics and engineers. A porous medium is a material containing pores, with the skeletal portion of the material called the "matrix" or "frame". The pores fill with a fluid(liquid or gas). Many natural substances such as rocks and soil, zeolites, biological tissues, and man-made materials such as cements and ceramics can be considered as porous media. Frozen soil is a porous medium composed of soil particles, pore ice, pore water and gas. Therefore, the propagation theory in porous media is an effective way to study wave propagation characteristics in frozen soil.

Mixture theory and continuum method can used to deal with the problems of multiphase porous media with irregular and complex internal structures (Yang and Yu, 2000). Therefore, Vardoulakis and Beskos(1986) established dynamic equations of unsaturated porous media based on mixture theory, and discussed the wave propagation characteristics of nearly saturated soil in low frequency ranges. Based on the mixture theory, Lu and Hanyga (2004, 2005) established the three-phase dynamic control equations of unsaturated soil, and analyzed the velocity and attenuation characteristics in unsaturated soil by numerical examples. Albers (2006, 2009) established a macro linear model based on the simplified mixture model and focused on analysis of the effects of saturation and excitation frequency on wave velocity and attenuation coefficient in unsaturated porous media. Assuming that the skeleton displacement, water head and temperature are the basic variables, Lu (2001) established two coupling seepage–elastic models,thermal–elastic and thermal–seepage. In addition, the propagation characteristics of saturated soil and unsaturated soil is studied in different angles, and the application of the mixture theory in wave problems is promoted to some extent (Zhouet al., 2008; Chenet al., 2012; Liuet al., 2014, 2017; Heet al., 2016). On the base of existing research of frozen soil, Miaoet al.(1995) introduced the mixture theory into the study of the thermodynamic properties of frozen soil to establish a complete constitutive relation.

Based on the Continuum Theory of Mixture, Zhou and Lai (2011) established the elastic wave dispersion equation of porous medium for frozen saturated soil and derived analytical expressions for propagation velocity and attenuation of elastic waves in frozen saturated soil. This theory applied the Bishop effective stress principle and ignored the coupling between components.

Based on the aforementioned results, mixture theory and continuum method play an important role in the study of wave propagation in unsaturated porous media, but most of the coupling between porous media is considered inadequate (Laiet al., 2003; Xuet al., 2011). With a porous medium consisting of soil skeleton, pore water, pore ice and gas, the equilibrium, deformation compatibility and energy balance conditions of frozen soil is influenced by the coupling of each component. Thus, the coupling between the components cannot be neglected. In this paper, by introducing entropy inequality to describe the coupling between the three-phase saturated soil material and using continuum porous medium method and mixture theory, we derived an analytical expression of wave velocity and attenuation and discussed the wave propagation characteristics in frozen saturated soils.

2 Wave equation for wave propagation in frozen saturated soil

In this paper, frozen saturated soils is considered as a three-phase porous media composed of soil skeleton, pore water and pore ice, designated by the superscripts "s", "w", and "i", respectively. The volume fraction of three-phase composition can be expressed as:

2.1 Continuity equation

According to the assumption of small deformation and no moment of momentum transfer between the phases, the continuity equation requires that the mass of any matter domain is a constant, so that the continuity equation of frozen saturated soil is:

whereβis used to denote an individual phase,β=w,i;

2.2 Motion equation

2.3 Entropy inequality

Considering the coupling between three phases of frozen saturated soil, entropy is introduced into the second law of thermodynamics in order to express the irreversibility of deformation process. For linear dynamic process of frozen saturated soil, the entropy inequality assumes the following form:

2.4 Constitutive relations

According to previous studies (Wei and Muraleetharan, 2006), Helmholtz free energy of frozen saturated soil can be written as:

For simplicity, the free energy equations of each component are solved with the following polynomials by the sum decomposition method:

The material derivative of pore water and pore ice are rewritten as:

Substituting Equations(4)–(6)into Equation(3),we obtained the control inequalities of thermodynamics:

To make Equation(7)absolutely valid, the following hypothesis needs to be satisfied:

Neglecting the free energy produced by strain deviator, the increment of specific free energy of each phase medium can be described as:

The state of frozen saturated soil before and after the treatment of external force (static equilibrium state

Assuming that the frozen soil deformation is in the small strain range, Equation(1)can be formulated as a linear mass conservation equation:

Using Equation(8)–(10), the following stressstrain relations is obtained:

2.5 Wave equation

The Bishop's effective stress formula for frozen saturated soil is expressed as:

In addition, to make Equation(7)absolutely established, capillary pressure can be expressed as:

If the volume ice content of frozen saturated soil is given, the relation between capillary pressure and pore water can be described as (Van Genuchten,1980):

We assume that the incremental form of capillary pressure can be described as:

By combining Equation(2), with(10)–(16), we can obtain wave equation of wave propagation in frozen saturated soil:

3 Elastic wave solution

The following general form for the dilatations can be assumed:

wherekpandksare complex number of P and S waves;ωis the excitation frequency.

Substituting Equation(19)into Equation(18)gives:

Assuming that Equation(20)has nonzero solutions, the characteristic equations of compressional wave and shear wave can be obtained:

Rearranging Equation(21)gives:

As is shown in Equation(22), three types of compressional (P1, P2 and P3) waves and one kind of shear wave (S) are found in frozen saturated soil. If the excitation frequencyωis given, velocity and attenuation coefficient can be obtained:

where Re and Im denote real and imaginary parts of the complex number, respectively.

4 Numerical results and analysis

Sections 2 and 3 indicate that the propagation velocity in frozen saturated soil is related to excitation frequency and volume ice content separately. In order to evaluate the influence of excitation frequency and volume ice content on propagation characteristics of frozen saturated soil, the physical parameters are listed in Table 1. Meanwhile, the following parameters are assumed: porosityns=0.2, excitation frequency range 0–1016Hz, volume ice content range 0.001–0.999, respectively.

Table 1 Physical parameters

Figure 1 Change law of P1, P2, P3 and S wave propagation velocity under different excitation frequency

Figure 2 shows the variation curves of P1, P2, P3 and S wave attenuation with excitation frequency in frozen saturated soil. One can see that with increasing excitation frequency, all wave attenuations of the four wave types increases, in which the attenuation of P2 and P3 is larger than that of P1 and S. The attenuation law of P1 and S waves is similar to that of magnitude, and it decays rapidly in intermediate-frequency. Take the condition of the excitation fre-

Figure 2 Change law of P1, P2, P3 and S wave attenuation under different excitation frequency

Figure 3 Change law of P1, P2, P3 and S wave propagation velocity under different volume ice content

Figure 4 Change law of P1, P2, P3 and S wave attenuation under different volume ice content

5 Conclusions

This paper is based on continuum porous medium method and mixture theory and considers soil skeleton, pore water and pore ice as fundamental quantities.By introducing entropy inequality to describe the coupling between the three-phase saturated soil material, we derived the analytical expressions of wave velocity and attenuation and discussed the propagation characteristics of waves in frozen saturated soils.

(1) There are three types of compressional waves(P1, P2, P3) and shear wave (S) in frozen saturated soil. The velocity of P1 and S is larger, attenuation is smaller, followed by P2, and P3 wave decays rapidly to zero.

(2) Under increasing excitation frequency, all velocities of four wave types increases and increases sharply in the excitation frequency range of 103–109Hz, but the wave velocity is almost unchanged outside this excitation frequency range. The attenuation of P2 and P3 waves is larger than that of P1 and S waves. The attenuation law of P1 and S waves is similar to that of magnitude, and the attenuation rate is greatly influenced by volume ice content.

(3) Under increasing volume ice content, the wave propagation velocity of P1 and S waves decreases dramatically when the volume ice content reached 0.6, and the velocity of P2 is increasing gradually. But the P3 wave velocity increases first and then decreases to zero with increasing saturation. The attenuation coefficient of P1 and S waves begins to increase gradually when the volume ice content is about 0.4. The attenuation coefficient of P2 wave increases first and then decreases with the increase of volume ice content. The attenuation coefficient of P3 wave increases with the volume ice content and decreases rapidly from extreme to zero.

(4) Due to the existence of ice particles in frozen soil, the wave propagation characteristics of frozen soil are different from those of unfrozen saturated soil. The wave propagation velocity in frozen saturated soil increases with volume ice content because waves travel faster in ice than in water. Under the one volume ice content, the phenomenon of the existence of an inflection point in the velocity curve of S wave is related to the characteristics of shear wave. As a result, in frozen saturated soils, volume ice content has an intimate relationship with the wave propagation characteristics.

This work is supported by the National Natural Science Foundation of China (No. 41271080 and No.41701060) and the funding of the State Key Labor-atory of Frozen Soil Engineering (No. SKLFSE-ZT-17).

Albers B, 2006. On results of the surface wave analyses in poroelastic media by means of the Simple Mixture Model and the Biot model.Soil Dynamics and Earthquake Engineering, 26(6–7): 537–547.DOI: 10.1016/j.soildyn.2006.01.007.

Albers B, 2009. Analysis of the propagation of sound waves in partially saturated soils by means of a macroscopic linear poroelastic model. Transport in Porous Media, 80(1): 173–192. DOI:10.1007/s11242-009-9360-y.

Chen WY, Xia TD, Chen W,et al., 2012. Propagation of plane P-waves at the interface between an elastic solid and an unsaturated poroelastic medium. Applied Mathematics and Mechanics, 33(7):781–795. DOI: 10.3879/j.issn.1000-0887.2012.07.001.

Hanyga A, 2004. Two-fluid porous flow in a single temperature approximation. International Journal of Engineering Science,42(13–14): 1521–1545. DOI: 10.1016/j.ijengsci.2004.04.001.

He PF, Xia TD, Liu ZJ,et al., 2016. Reflection of P1 wave from free surfaces of double-porosity media. Rock and Soil Mechanics,37(6): 1753–1761. DOI: 10.16285/j.rsm.2016.06.028.

Lai YM, Zhang LX, Xu WZ,et al., 2003. Temperature features of broken rock mass embankment in the Qinghai-Tibetan railway.Journal of Glaciology and Geocryology, 25(3): 291–296. DOI:10.3969/j.issn.1000-0240.2003.03.009.

Li DQ, Wu ZW, Fang JH,et al., 1998. Heat stability analysis of embankment on the degrading permafrost district in the East of the Tibetan Plateau, China. Cold Regions Science and Technology,28(3): 183–188. DOI: 10.1016/S0165-232X(98)00018-4.

Li XW, Li XY, 1989. Wave propagation with mass-coupling effect in fluid-saturated porous media. Applied Mathematics and Mechanics, 10(4): 321–327. DOI: 10.1007/BF02017772.

Liu B, Su Q, Zhao WH,et al., 2017. Reflection and transmission of P waves at interfaces of saturated sandwiched coarse granular structure. Journal of Southwest Jiaotong University, 52(2): 280–287.DOI: 10.3969/j.issn.0258-2724.2017.02.010.

Liu ZJ, Xia TD, Zhang QF,et al., 2014. Parametric studies of propagation characteristics of bulk waves in two-phase porous media. Rock and Soil Mechanics, 35(12): 3443–3450. DOI: 10.16285/j.rsm.2014.12.001.

Lu HL, 2001. Continuum theory of mixture for freezing and thawing of water-saturated porous media. Journal of Southwest Jiaotong University, 36(6): 599–603. DOI: 10.3969/j.issn.0258-2724.2001.06.011.

Lu JF, Hanyga A, 2005. Linear dynamic model for porous media saturated by two immiscible fluids. International Journal of Solids and Structures, 42(9–10): 2689–2709. DOI: 10.1016/j.ijsolstr.2004.09.032.

Miao TD, Zhu JJ, Ding BY, 1995. Essay on constitutive relation of wave propagation in saturated porous media. Chinese Journal of Theoretical and Applied Mechanics, 27(5): 536–543. DOI:10.6052/0459-1879-1995-5-1995-464.

Ming F, Zhang Y, Li DQ, 2016. Experimental and theoretical investigations into the formation of ice lenses in deformable porous media.Geosciences Journal, 20(5): 667–679. DOI: 10.1007/s12303-016-0005-1.

Van Genuchten, 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal, 44(44): 892–898.

Vardoulakis I, Beskos DE, 1986. Dynamic behavior of nearly saturated porous media. Mechanics of Materials, 5(1): 87–108. DOI:10.1016/0167-6636(86)90017-7.

Wei CF, Muraleetharan KK, 2006. Acoustical characterization of fluidsaturated porous media with local heterogeneities: theory and application. International Journal of Solids and Structures, 43(5):982–1008. DOI: 10.1016/j.ijsolstr.2005.06.008.

Xu XT, Lai YM, Liu F,et al., 2011. A study of mechanical test methods of frozen soil. Journal of Glaciology and Geocryology, 33(5):1132–1138.

Yang SY, Yu MH, 2000. Constitutive descriptions of multiphase porous media. Acta Mechanica Sinica, 32(1): 11–24. DOI:10.3321/j.issn:0459-1879.2000.01.002.

Zhou FX, Lai YM, 2011. Propagation characteristics of elastic wave in saturated frozen soil. Rock and Soil Mechanics, 32(9): 2669–2674.DOI: 10.3969/j.issn.1000-7598.2011.09.018.

Zhou XM, Sun MY, Xia TD,et al., 2008. Effect of the fluid on propagation characteristics of compressible waves in saturated soils. Journal of Harbin Institute of Technology, 40(6): 974–977.DOI: 10.3321/j.issn:0367-6234.2008.06.031.