PermeabilityCoefficient ofUnsaturatedLoess and ItsGaseous and Liquid Water Migration Modeling

Watermigration changes the water distribution of loess and significantly affects the engineering properties of loess. However, because of the characteristics of loess such as water sensitivity and special structure, it is still difficult to understand the water migration trend of loess. ,e objectives of this study are to investigate the water diffusion coefficient and soil-water characteristic curve (SWCC) of loess, focusing on the effect of temperature and density, and establish a mixed migration equation for gaseous and liquid water to develop a new calculation method of loess water migration.,e results show that the density of loess with larger liquid water content has a more significant effect on the diffusion coefficient than that of less liquid water content. Furthermore, the density significantly affects the matric suction of unsaturated loess, while the temperature change at normal temperature has slight effect. Based on the data obtained in this study, the soil water permeability coefficient was obtained. Moreover, considering the characteristics of mixed migration of gas and liquid water in unsaturated loess, an equation was derived for the gas-liquid water transfer in unsaturated loess. ,e calculation results are consistent with the experimental results of water migration experiment. Based on the calculation results of water migration, the trend of water migration and proportion of gaseous water during migration were analyzed.


Introduction
Permeability coefficient is important to evaluate the permeability of unsaturated soil. e analysis and calculation of water movement require the measurement of permeability coefficient of unsaturated soil. e permeability coefficient of unsaturated soil has been extensively studied. Jeong et al. [1] studied the effect of rainfall-induced volumetric liquid water content on matric suction. Lemoubou et al. [2] developed an estimation method for unsaturated soil permeability parameters considering the effect of volumetric liquid water content. Zhang et al. [3] studied the effect of instantaneous volumetric liquid water content and saturation on the soil suction of thickened oil sand tailings. Kaewsong et al. [4] studied the effect of suction change caused by the change in saturation degree on the small strain stiffness of unsaturated soil. Tian et al. [5] studied the effect of volumetric liquid water content on the permeability coefficient of unsaturated loess under steady state. Trabelsi et al. [6] studied the suction of clays with different saturations.
Permeability parameters vary with the change in soil types [7][8][9]. Unsaturated loess, one of the most widely distributed soil types, is a special soil with water sensitivity and macropore structure [10][11][12]. Zhang et al. [13] reported that the permeability coefficient of unsaturated loess is affected by pore characteristics, and the pore characteristics of loess with different densities and water contents under different pressures were studied. Based on the study of macropore characteristics of unsaturated loess and effect of consolidation pressure on permeability coefficient, Wei et al. [14] assumed that the pore characteristics significantly affect the permeability coefficient, and the size and distribution of pore determine the permeability coefficient. By conducting constant-head permeability experiments and mercury injection porosimetry on loess specimens of different dry densities, Wang et al. [15] assumed that there is a threshold hydraulic gradient below which no flow is observed, and there is a critical hydraulic gradient below which the relationship between hydraulic gradient and seepage velocity is nonlinear, which is closely related to pore size distribution.
rough theoretical analysis, field test, and software simulation, Wang et al. [16] studied the effect of unsaturated loess permeability coefficient on the grouting effect of collapsible loess stratum. Liu et al. [17] used microstructure models to study air permeability.
Many factors, especially temperature, significantly affect the permeability of unsaturated loess [18][19][20][21]. Owing to a large daily variation of temperature in the loess plateau, the temperature difference between the upper and lower layers of shallow soil is large. With climate change, the temperature change of shallow soil is also relatively large. e change in loess matrix suction caused by temperature affects the mechanical properties of loess and water migration. In addition, because water flows through the pore space of unsaturated loess, the scale of pore space determines its permeability, i.e., the permeability coefficient of unsaturated loess is not constant. At present, some studies have been conducted on the seepage parameters of unsaturated loess, but it is still difficult to understand the trend of water permeability. It is essential to study the water permeability parameters of unsaturated loess. e objectives of this study are to investigate the water diffusion coefficient and soil-water characteristic curve (SWCC) of loess, focusing on the effect of temperature and density, and derive a mixed migration equation for gaseous and liquid water, thus developing a new calculation method for loess water migration. e water diffusion coefficient of unsaturated loess considering density effect was studied at various densities (1.2, 1.4, 1.5, and 1.6 g/cm 3 ). Furthermore, the density and temperature effects of SWCC were evaluated by measuring the matric suction at various densities (1.2, 1.4, 1.5, and 1.6 g/cm 3 ) and temperatures (5, 15, 25, and 35°C). Based on the data of water diffusion coefficient and SWCC, the soil water permeability coefficient was obtained. Moreover, considering the characteristics of mixed migration of gas and liquid water in unsaturated loess, an equation was derived for gas-liquid water transfer in unsaturated loess. e data of water diffusivity and SWCC of loess were substituted into the equation for calculation.

Experimental Study on Water Diffusion Coefficient
To measure the water diffusivity of unsaturated loess, a water diffusivity experiment was carried out. Disturbed loess was taken from an engineering site in Xi'an at a depth of 5 m below the surface, and the property parameters of loess are shown in Table 1. First, compacted loess samples of different densities were prepared indoors. en, a horizontal soil column infiltration was used to test the water diffusivity of unsaturated loess. Figure 1 shows the experimental apparatus for water diffusivity. Markov bottle can supply water for the loess samples. e sleeve of columnar loess simple, a plexiglass cylinder with a length of 5 cm and diameter of 4.2 cm, can be removed from each other to facilitate easy disassembly at the end of experiment to determine the water content. Baffle plates were added at the head and end of soil column and fixed with screws. e total length of soil column is 75 cm. At the inlet boundary, the water chamber with an exhaust hole can keep the water level constant and the soil saturated. A filtering layer with low bubble pressure and high conductivity was installed between the inlet chamber and soil column to ensure that the only driving force in the soil is matrix suction.
First, enough samples (air-dried soil) were prepared, and loess with uniform bulk density was loaded in each cylinder. Second, the screw was tightened, and the soil column was horizontally placed after the soil column was installed. en, the inlet chamber was instantaneously filled with water, and the exhaust hole was kept open until the air in the inlet chamber was completely exhausted. Subsequently, we start to record the elapsed time and the initial water level reading. Finally, the water supply was stopped, and the fastening screw was loosened after the experiment time.
e mass water content of each section of loess was measured by weighing and converted to volumetric water content, and the experiment duration and total water inflow were recorded.
According to the time-dependent water content in the soil column along the length direction, the water diffusion coefficients were obtained using the following equation: where D(θ) is the soil water diffusion coefficient (cm 2 /min); λ � xt − (1/2) is the transformation parameter; x is the distance from a point in the soil trough to the interface between soil and water tank (cm); t is the test duration (min); θ 0 is the uniform initial volume water content of soil column (cm 3 /cm 3 ); θ is the volume water content of soil column during the diffusion (cm 3 /cm 3 ).  2 Advances in Civil Engineering Figure 2 expresses the difference in the magnitude of permeability coefficient, and semilog plot was used in the curve, making the original smooth curve wavy. In the logarithmic coordinate, the change in the curve of water diffusion coefficient with volumetric liquid water content is an approximate straight line. us, the water diffusion coefficient increases with the increase in volumetric liquid water content. e higher the volumetric liquid water content, the faster the rate of increase of water diffusion coefficient. With the increase in volumetric liquid water content, the water diffusion coefficient of unsaturated loess increases by several orders of magnitude. Besides, the water diffusion coefficient of unsaturated loess increases with the decrease in dry density, and it is more obvious when the dry density is low. erefore, the effect of density should be fully considered when determining the water diffusion coefficient of loess.

Experimental Study on SWCC of Unsaturated Loess
To measure the SWCC of unsaturated loess, a matric suction experiment was carried out. e physical properties of soil samples used in the experiment are the same as those of loess used in the water diffusion experiment. e SWCC was measured using a high-speed centrifuge (Hitachi cr21g PF). First, four ring cutters were loaded with the same day density of compacted loess, and those samples were saturated with distilled water. en, the saturated sample was placed in a centrifuge and dehydrated at a specific speed. Because of the centrifugal force, the water in soil sample will be thrown out into the container below. e matric-suction-led adsorption of water in soil could be counteracted by centrifugal force. is behavior causes the numerical equality between centrifuge force and suction of soil sample after centrifuge dewatering; that is, the centrifuge force of soil can be regarded as matric suction. us, matric suction can be represented by a suction head and calculated from the speed and radius of centrifuge equation (2). e measurement accuracy of matric suction in this study is 1 kPa.
By measuring the volume of water thrown out, the volume of residual water in soil sample can be calculated, and the volumetric liquid water content of soil sample can be obtained after the experiment. Finally, by changing the speed of centrifuge, the water content of loess under different suctions can be measured. e final volumetric liquid water content and matric suction are the average values of four samples, and the SWCC of loess with a certain density can be obtained. e above experiments were repeated with loess of other densities. e relationship between volumetric liquid water content and suction head of unsaturated loess with different bulk densities can be obtained. e suction head h is where r is the centrifugal radius (cm), rpm is the rotational speed (r/min), and h is the suction head (cm). Figure 3 shows that temperature has no significant effect on matric suction. For the soil sample with a certain volume moisture content, the suction change caused by temperature change has a certain value. is is small at a high water content and can be ignored but is larger at a low water content. e larger the temperature difference, the greater the change of suction. When studying the water migration of unsaturated soil, the difference in moisture content at different temperatures is not significant when the suction is the same. erefore, the effect of temperature on suction can be ignored. Figure 4 shows that when the soil density ratio is small, the change in suction with the density of unsaturated loess at different volumetric liquid water contents is small. When the soil density is large, the change in suction with the density of unsaturated loess at different volumetric liquid water contents is relatively large. e effect of density on matric suction is very significant, mainly caused by the change in density, pore size of soil, and curvature radius of membrane water change.

Calculation Method for Permeability Coefficient of Unsaturated Loess
First, the parameters of SWCC model were analyzed and determined. Soil density significantly affects the matric suction of unsaturated loess; therefore, the density factor must be considered when fitting the equation. Although temperature can hardly affect the matric suction of unsaturated loess, the temperature factor can be ignored in the fitting equation. e models used to represent SWCC are usually the Gardner model [22], VG model [23], and Fredlund model [24]. In this study, the SWCC is represented by VG model: where θ r is the residual volumetric liquid water content; θ s is the quasi-saturated volumetric liquid water content; and α, m, and n are coefficients; m � 1 − (1/n). e specific water volume of unsaturated loess is − (dθ/dh), obtained by deriving θ(h) and represented by C(θ) in the following equation: e experimental results of four types of density of unsaturated loess were fitted and analyzed. θ r , θ s , α, n, and m can be written as follows: e specific water capacity of unsaturated loess can be obtained by substituting θ r , θ s , α, n and m into equation (4). en, regression analysis was carried out on the experimental values of water diffusivity D(θ) of unsaturated loess with four densities. e regression equation is where p 1 , p 2 , p 3 are the parameters varying with density and θ is the volumetric liquid water content. Regression analysis was carried out on the change of p 1 , p 2 , and p 3 with density. e regression equations are shown in the three following equations: By substituting C(θ) and D(θ), the permeability coefficient of unsaturated loess considering the effect of density was obtained using the following equation: To verify the equation, the parallel method was used in the water diffusion experiments and matrix suction experiments.
e calculated results were compared with the measurement results obtained by parallel experiments. e correlation coefficients are shown in Table 2, indicating that equation (10) can be used to fit the experiment results appropriately.

Analysis of Permeability Coefficient Equation
Equation (10) was obtained based on the results. It was applied to calculate the permeability coefficient of loess with different densities and water contents to explore the change in permeability coefficient of unsaturated loess with density and water content. e calculation results are shown in Figures 5 and 6. Figure 5 shows the variation of permeability coefficient with the water content of unsaturated loess with different densities. Figure 6 shows the change in permeability coefficient with the density of unsaturated loess at different water contents. e semilog plot was used. Figure 5 shows that the change in water content of unsaturated loess can lead to several orders of magnitude changes in its permeability coefficient; under the condition of constant dry density, the relationship between permeability coefficient and water content in semilog plot is approximately linear, indicating that the value of permeability coefficient increases monotonously with the increase in water content, similar to the increase in water diffusion coefficient with water content. From the absolute change, when the water content is large, the change in permeability coefficient with water content is relatively large, and when the water content is small, the change in permeability coefficient with water content is relatively small. Because of the use of order coordinate, a wave-type curve appeared in the calculation curve, which reflects that the increased rate of permeability coefficient increases with the increase in water content; i.e., the higher the soil moisture content, the faster the increase in speed.  Figure 6 shows that when the soil density is small, the change in permeability coefficient with the density of unsaturated loess at different water contents is small. When the soil density is relatively large, the change in permeability coefficient with the density of unsaturated loess at different water contents is relatively large. e permeability coefficient of dense loess is sensitive to the change in density.

Water Migration of Unsaturated Loess
To apply the obtained permeability parameters to the calculation of water migration, a one-dimensional (1D) water transfer model of unsaturated loess was established. e obtained permeability parameters were substituted into the model to obtain the calculated values of water migration. According to the measured results of water migration of unsaturated loess, the calculated values of water migration  based on the permeability parameters of loess were compared with the measurement values of water migration experiments. e comparison results can be used to verify the diffusion coefficient and SWCC of loess measured above.

Experimental Setup for Water Migration.
Experiments were conducted on the water migration of unsaturated loess, and the results of water migration were obtained. First, Arid soil (with a small water content) and wet soil (with a large water content) were prepared and then placed in a PVC pipe with a length of 10 cm and sealed at one end. en, they were docked in pairs. Figure 7 shows the experiment apparatus of water migration, and Table 3 shows the physical properties of loess samples.

eoretical Model for Gaseous Water
Migration. e gaseous water present in loess is mainly transported by the diffusion of saturated gaseous water. Matrix suction makes pore water form a concave meniscus, affecting the saturated gaseous water pressure of gaseous water migration and diffusion. According to the law of conservation of mass and Fick's first law, a 1D diffusion equation of gaseous water in unsaturated loess equations (11) and (12) can be obtained.
where w a is the moisture content of gaseous water; J x is the water air flow rate in x direction; D x is the diffusion coefficient of gaseous water in x direction; S is the saturation of soil; n is the porosity; R v is the gaseous water constant with a value of 0.4615 (KJ/(kg · K)); T is the absolute temperature; P v0 is the gaseous water pressure in the pore in the saturated state, where the value is related to the temperature and matric suction of soil. Considering the relationship between matric suction and P v0 , P v0 can be expressed using Kelvin equation in the two following equations: where P v is the saturated gaseous water pressure on free surface; e is the surface tension of curved liquid surface, which is determined according to the SWCC obtained in this study; r is the radius of curvature; ρ v0 is the density of gaseous water; ρ w is the density of water. In unsaturated soils, the matric suction u a − u w causes the interface between gaseous water and liquid water to bend, and it can be obtained using the following equation: By combining equations (13) and (15), P v0 can be represented as follows: P v is also affected by temperature and can be represented as follows:   Advances in Civil Engineering where B and F are constants with values of 2.315 × 10 8 kPa and 2489.3 (kJ/kg), respectively; u a − u w is the matric suction of unsaturated soil. By combining equations (16) and (17), P v0 can be represented as follows: D x can be obtained from the following formula: where D v is the gaseous water diffusion coefficient (cm 2 /s) with the reference system of aerated pores in soil (D x taking the soil as the reference system); D a is the gaseous water diffusion coefficient in still air (cm 2 /s), and the value of D a is 2.12 × 10 − 5 (T/273.15) 2 , a constant in this study; a is the curvature, 0.66 in this study; θ v is the volume moisture content (cm 3 /cm 3 ); η is the porosity (%); and ρ is the saturated density of gaseous water, which is related to temperature.

eoretical Model for Liquid Water Migration.
Liquid water migration conforms to Darcy's law. e liquid water migration in a loess area is unsteady, and the water diffusion rate and matrix suction affect the permeability coefficient. Based on Darcy's law and law of conservation of mass, a 1D equation of liquid water transport in unsaturated loess can be obtained.
where θ L is the volumetric liquid water content of liquid water, which should be converted into the mass water content of liquid water in calculation. K is the permeability coefficient of liquid water in the direction of seepage, which has been measured in this study; and ψ m is the matrix potential of liquid water, which was obtained from the SWCC measured in this study. By combining equations (11) and (20), a governing equation for gaseous and liquid water flow can be obtained using the following mass conservation equation: where θ is the mixed volumetric water content, namely, the sum of liquid and gaseous water contents.

Boundary and Initial
Conditions. e boundary conditions are presented using the following equation: e initial conditions were obtained using the following equation: where w is the mass water content (w l + w a ) (%) and l is the distance from the arid end (cm).

Calculation of Water Migration.
e finite element equations of liquid water and gaseous water migration were transformed using a backward difference scheme. e boundary conditions are consistent with the boundary conditions of water transfer experiment. Combined with the model parameters obtained by the author, the distribution laws of gaseous and liquid water at different migration times were iterated. Considering the interaction of gaseous water migration and liquid water migration, the migration amounts of gaseous water and liquid water in the i-th time step were calculated simultaneously. After the iteration converges to the optimal value, the gaseous water content and liquid water content were added to obtain the total water content. Using the total water content of the i-th time step, the migration amounts of gaseous water and liquid water in the i+1-th time step can be calculated. In this manner, the results of gas-liquid mixed water migration can be obtained under a predetermined water transport time. e calculation process is shown in Figure 8. To determine the accuracy of calculation results, dt and ds were set as 2.4 s and 0.001 m, respectively, by trial calculation. e calculation program was compiled by FORTRAN.

Comparison of Calculated and Experimental Results.
e experiment soil samples with a length of 200 mm were divided into 100 sections. e step length dx is 2 mm, and the time step dt is 2.4 s. e experiment was carried out under different conditions. Figure 9 shows that the    Advances in Civil Engineering content is close to the experimental value, but the change in calculated water content is still greater. is is mainly because the model calculation ignores the hysteresis effect of water migration, and the contact area of arid and wet soil samples has a moisture content gradient, where the water migration is intense, and the hysteresis effect of water migration is obvious. e hysteresis effect is most obvious in the initial stage of water migration, but, with the increase in migration time, the hysteresis effect gradually weakens, and the fitting value and measured value converge. On the whole, the water content distribution obtained by the coupling calculation of water transfer between gas and liquid with the measured permeability parameters in this study is consistent with the obtained water content distribution in the experiment. Figure 9 shows that water migration flowing from wet soil to arid soil results in significant redistribution of water between arid and wet soils. e change in water content can be divided into three stages according to the time: In the first stage (days 0-8), the water content of loess near the arid wet contact surface changes significantly, but the change in soil moisture content far away from the contact surface is not obvious. In the second stage (days [8][9][10][11][12][13][14][15][16], the water content of loess near the arid wet contact surface remains stable, and the water content far away from the contact surface changes. In the third stage (after day 16), the stable stage of water migration is reached, and the water content in different places does not change. In addition, the change in water between arid and wet soils is antisymmetric from the arid wet contact surface. To facilitate the analysis, the arid soil (water diversion end) was taken for analysis. In the first stage (days 0-8), the water content of right side of arid soil increases rapidly. e closer the distance between the arid and wet contact surfaces, the greater the increase in water content. However, the water content of left side of arid soil (far away from the arid wet contact surface) changes slightly.
is is mainly because, in the early stage of water migration, only the water content gradient near the arid wet interface is used as the driving force, while the moisture content in the soil is uniform, and no moisture gradient acts as the driving force. erefore, the moisture of wet soil flows into the right end of arid soil from the contact surface and accumulates at the right end. With water accumulation in the right end, a moisture content gradient is formed on the left and right sides of arid soil, and the water migration enters the second stage (days [8][9][10][11][12][13][14][15][16]. Driven by the gradient, the water present on the right side moves to the left side, and the water content of the left side soil increases rapidly. However, because of the water supply from the wet soil on the right side, the water content of right side remains unchanged and does not decrease. With the increase in water content in the right end, the water contents at the left and right ends gradually become equal. With the disappearance of driving force, the water migration enters a stable stage, and the water migration enters the third stage (after day 16). Figure 10 shows that gaseous water migration is the main water migration in total water migration. According to the water migration experiment on the 8th day, the gaseous water migration accounts for about 70% of total migration. On the 12th day of water migration experiment, the proportion of gaseous water migration decreased but still accounted for 65% of total migration. On the 16th day of water migration experiment, the proportion of gaseous water migration further decreased, accounting for about 50% of total migration. us, the migration of gaseous water was faster than that of liquid water. With the increase in time, the migration of gaseous water weakened, and the migration of liquid water increased. However, in general, the Advances in Civil Engineering migration of gaseous water accounts for the main part of total migration.

Conclusion
(1) e relationship between water diffusivity and water content of loess with different densities was obtained using horizontal soil column infiltration method, and the SWCC of loess with different densities was obtained using a high-speed centrifuge method. e experiment results show that when the water content is large, the effect of density on the diffusion rate and matric suction of loess is very significant.
(2) A method was developed to determine the permeability coefficient of unsaturated loess considering the effect of density by regression analysis. e variation trend of permeability coefficient of unsaturated loess with density and water content was analyzed.
e permeability coefficient increases monotonically with the increase in water content, and the permeability coefficient of dense loess is sensitive to the change in density. unsaturated loess. e water diffusivity and SWCC of loess obtained by experiment were substituted into the equation for calculation. e calculated results are consistent with the experimental results of water migration experiment. is shows that the model of water migration established in this study can be used to calculate the water migration of unsaturated loess, and it is reasonable to apply the measured water diffusivity and SWCC in this study to the calculation of water migration.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.