Abstract
Porosity change is a common characteristic of natural soils in fluid-solid interaction problems, which can lead to an obvious change of the soil-water retention curve (SWRC). The influence of porosity on soil water retention phenomena is investigated by a theoretical model and an experimental test in this study. A model expressing the change in suction with porosity and effective saturation is put forward theoretically. The model is based on an idealization of three-phase porous materials, the pore structures of which are homogeneous and isotropic. It accounts for the porosity effect on soil water retention, using four parameters with clear physical meanings. The presented model can obtain the SWRC at any porosity, which will reduce the test number required in characterizing the hydraulic behavior of soil. A laboratory experiment for loamy sand with different porosities is performed. The test results show that suction has a significant variation with changes in porosity and decreases with the increase of porosity. The formulation is verified by both the test data and the literature data for FEBEX bentonite and Boom clay. The very good agreements between measured and predicted results show that the SWRC model is reliable and feasible for various soils.
1. Introduction
The soil-water retention curve (SWRC) is defined as the relationship between matric suction and the degree of saturation in unsaturated soils. The SWRC relation is a key hydraulic property to describe fluid flow phenomena in unsaturated soils, which can be affected by the porosity or the density of the soil [1, 2]. A variation in soil porosity is common in fluid-solid interaction problems [3, 4], and this change will result in an obvious change of SWRC [5]. Sometimes, it needs to predict the SWRCs for soil with different porosities. However, many traditional SWRC models [6–11] do not consider the effect of porosity on SWRC, which may cause an inaccurate result. Therefore, it needs to establish the relationship between porosity and SWRC in agricultural engineering and geotechnical engineering, especially involving water-induced landslides [12–14] and water-rock coupling engineering [15–17].
Various SWRC models have been proposed to take the influence of porosity into account by relating suction and/or the pore size distribution index to void ratio or density. The first group of the models linked the current void ratio to the initial void ratio using a volume change law which describes the relation between suction and void ratio [18–22]. It is difficult to establish a direct formulation between volume and suction, which usually involves the hydraulic and mechanical processes of soil. Some models used an empirical volume function to quantify the influence of initial porosity on SWRC based on experimental findings [21]. Even for the simplest models [23], it still requires a lot of additional parameters to describe the variation of soil volume. The other approaches describe the porosity dependence by shifting SWRCs with a porosity-dependent air entry value or pore size distribution index [3, 24–28]. Such approaches can be used in SWRC laws which contain the item of air entry value, such as the Brooks and Corey model [7]. But unsaturated soils are three-phase mixtures of solid, water, and air. Water and air flow in structured soils depend not only on soil texture but also on pore shape and pore distribution [29]. Some important intrinsic relations for the pore systems may not be satisfied. It is necessary to consider the influence of pore distribution, shape, and size on SWRC from the mixture theory.
In this paper, a SWRC model accounting for the effect of porosity is presented, assuming the unsaturated soil as a continuous, isotropic, and homogeneous three-phase pore mixture. This equation is expressed by a law of effective saturation, suction, and porosity, with only four parameters. An experimental test is also conducted for loamy sand to study the porosity effect on water retention. The proposed model is verified by the test data for loamy sand with different porosities, as well as the literature data for FEBEX bentonite and Boom clay.
2. Derivation of the Porosity-Dependent SWRC Model
Suction, , can be expressed as the relation between pore gas pressure () and pore water pressure ():
in soil can be set as the atmospheric pressure, for gas flow is relatively free to the water flow ([30, 31]). Therefore, it only needs to determine .
2.1. Pore Water Pressure
Although the soil pore structure is composed of nonuniform pore size and the pore distribution is quite complex, most soils can be assumed to be a continuous, isotropic, and homogeneous porous medium, in which the distributions of pore water and air are also isotropic and homogeneous at a macroscopic scale [29]. Therefore, on per unit area of two-dimensional (2D) cross-section cutting through the soil for any angle, the total pore area , pore gas area , and pore water area are constant. The macroscopic porosity and saturation can be described by where , , and are the total volume of pore, water, and soil, respectively.
The water in the pore mainly includes the interfaces of water-solid and water-gas and the water between the two interfaces, as shown in Figure 1(a). In order to calculate on any cross-section, the pore shape can be arbitrary. Here, the pores are idealized as ellipses for their shape can be round or flat. An elliptical pore randomly distributes on arbitrary cross-section. Each pore contains an identical elliptic air bubble (Figure 1(b)) [29]. The total normal force of a given cross-section is expressed as where represents water surface tension, is the angle between the tangent plane and the section plane at point of the water-gas interface (Figure 1(b)), and represents the differential length of . is the area occupied by water but not including the water-gas interface. The differential area of is represented by represents the static water pressure in the area of .
(a)
(b)
According to the assumption of isotropy and homogeneity of medium, the angle at any point at the water-air interface always has . The part, , in Eq. (3) can be expressed as
Another part, , in Eq. (3) is rewritten as
According to the theory of surface science, the surface tension of liquid can be regarded as a constant when temperature is constant [32]. can be described as where ; the unit for is N/m.
The shapes of elliptic gas bubbles and pores are represented by two semiaxes, and for pores () and and for gas bubbles (). A pore contains gas bubble, so there are and . Defining and , where and , the macroscopic porosity and saturation can be expressed as
The equivalent pore number per area for the cross-section is represented by , which is expressed as . According to Eqs. (7) and (8), one can obtain and .
The perimeter of an elliptical air bubble can be written as . The total perimeter length of the air bubbles on the cross-section:
Eq. (6) can be rewritten as
2.2. Formulations of the New SWRC Model
Substituting Eq. (10) into Eq. (1) yields
The equation of suction is composed of the functions of pore density, saturation, shape of air bubbles, and water stress. The influence of porosity on suction is reflected through the parameters (porosity), (pore density), and (shape coefficient of air bubble). (pore density) is related to the pore structure of soil and can be written as a function of porosity. A simple function is used to express the relation. is related to saturation, and its effect can be expressed by saturation. So, Eq. (11) can be rewritten as where , is in units of N/m2, and the other parameters are dimensionless. and are parameters changed with .
Defining yields
can be obtained from the conditions that and :
Hence, one has
The equation is not a new one and has been used by many people (including Brooks and Corey and van Genuchten). describes the ratio of current saturation to the maximum. Noting that represents the saturation when and represents the saturation when .
The function represents the effect of the shape of air bubbles. It can be approximately expressed by a simple power law function through a good curve fitting method (Figure 2):
The air bubble shape, represented by , should vary with , , (temperature), and some other factors. For the difficulty in measuring directly in practice or developing analytical function form of , a function for is used here for simplicity, , where and are parameters related to the concerned soil. Eq. (13) can be simplified as
The proposed model involves four parameters, . is related to the surface tension. The unit of is N/m2, and the other parameters are dimensionless. is connected to pore density and related to porosity. represents the combined effect of temperature, hysteresis effects, and other factors, and is only related to porosity.
3. Validation
3.1. Sample Preparation and Test Procedure
In order to verify Eq. (17), a pressure-plate system is used to determine SWRCs at different porosities. The setup of the experimental system is shown in Figure 3. A cylindrical soil sample (70 mm in diameter and 20 mm in height) enclosed tightly in a cutting ring is placed on a ceramic disk, and the ceramic disk (with air entry value 500 kPa) connects well with the bottom. With the application of air pressure, water flows out from the soil sample to a chamber and the mass of water is measured by a balance.
A kind of loamy sand (83.92% sand, 14.19% silt, and 1.89% clay) extracted from the Three Gorges Region, China, was used as test material. According to the specification of soil test, three samples were prepared and labeled as L1, L2, and L3 (Figure 4), and the corresponding porosities were 0.4666, 0.4346, and 0.4097, respectively. The test started with saturated soil and was carried out in room temperature (20°C). The degree of saturation corresponding to each suction step was calculated from weighing the water outflow, the soil sample over the test, and the dry sample after drying in an oven.
3.2. Validation against the Experimental Results
The test data for loamy sand with all porosities was adopted to verify Eq. (17). There are four parameters, , in the proposed SWRC model. The parameters are calibrated by using a nonlinear least squares (NLS) fitting method: , , , and , respectively.
Figure 5 shows the comparison between the model-predicted suctions and the test data for loamy sand. The predicted suctions are in good agreement with the test data, with most discrepancies less than 10%. The measured and predicted SWRCs using the proposed model at different porosities are shown in Figure 6. It can be seen that the predicted SWRCs by the proposed model match well with the experimental data, which indicates that Eq. (17) successfully reproduces the influence of porosity on the SWRC. Figure 6 also shows that both the calculated and measured suctions decrease with the increase of porosity, and the shape and position of SWRCs change obviously with different porosities.
3.3. Validation against Literature Data
In order to further examine the applicability of the proposed model, the experimental data for two different soils found in literature are used.
The first is FEBEX bentonite (clay: 92%; silt and sand: 8%), which is a kind of engineered clay [29]. Three dry densities were chosen here; they are g/cm3, g/cm3, and g/cm3, respectively. The test temperature was 20°C. The porosity was calculated as , where g/cm3. So, the corresponding porosities of the three samples are , , and , respectively.
The second is Boom clay [33]. This material was obtained by compacting natural Boom clay. The Boom clay (10–20% smectite, 20–30% illite, and 20–30% kaolinite) has a plastic limit of 29%, a liquid limit of 56%, and a specific gravity of 2.7, and half of the particles are smaller than 2 μm. The tested samples were prepared at dry densities of 1.37 g/cm3 and 1.67 g/cm3, and the corresponding void ratios are 0.932 and 0.591, respectively. So, the porosities of the two samples are 0.4824 and 0.3715, respectively. The test temperature was kept at 22°C.
The four parameters, , in Eq. (17) are calculated by using the NLS method and listed in Table 1. In Figure 7, the experimental data for the two different soils are fitted with Eq. (17). It reveals very good agreements between the model predictions and experimental data, with negligible discrepancies.
(a)
(b)
4. Conclusions
A theoretical model is developed to describe the influence of porosity on SWRC. The model is derived for a three-phase mixture with an idealized pore structure. The idealized air bubble shape and pore system meet the equivalence of macroscopic physical effects. On this basis, a SWRC model is proposed, where suction is expressed as a function of porosity and effective saturation. There are four characterization parameters in the model, , with clear physical meaning. is related to the surface tension; is related to porosity; represents the combined effect of temperature, hysteresis effects, and other factors; and is only related to porosity. An experimental verification of this model was carried out on SWRCs of loamy sand. The model predictions are close to the test results, which also show that suction decreases with the increase of porosity, and the shape and position of SWRCs change obviously with different porosities. The predictions of the proposed model were also compared with the experimental results for FEBEX bentonite and Boom clay in published articles. The good agreements show that the SWRC model is reliable and feasible for wide soils.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by the National Key Research and Development Program of China (grant number 2017YFC1501100), the National Nature Science Foundation of China (grant numbers 51939004 and 51279090), and the Open Fund of Key Laboratory of Geological Hazards on Three Gorges Reservoir Area, Ministry of Education, China Three Gorges University (grant number 2020KDZ11).