Abstract
The consolidation of clay layers is of great significance for groundwater environmental protection, groundwater storage utilization, and land subsidence. In this study, the governing equation for the excess pore water pressure during the non-linear consolidation process of clay layers under load conditions is obtained based on the one-dimensional non-linear consolidation theory. Analytical solutions are then derived for clay layers with single or double drainage caused by the dissipation of the excess pore water pressure. With these analytical solutions, the groundwater dynamics and deformation of the clay layer are analyzed. Correspondingly, a type curve method is proposed to calculate the hydraulic parameters of the clay layer through laboratory experiments, which verifies the reliability of the analytical solutions. The study results show that the deformation of the clay layer predicted by the non-linear consolidation theory is smaller than that predicted by the linear consolidation theory. The deformation of the clay layer increases with the increase in the thickness of the clay layer, the compressive index, and the overburden load, while it decreases with the increase in the initial void ratio and the initial effective stress. The stable time, at which the consolidation of the clay layer is completed, increases with the increase in the compression index and the thickness of the clay layer, while it decreases with the increase in the initial void ratio, the initial effective stress, and the initial hydraulic conductivity. It does not vary with the load pressure. Conclusively, the deformation prediction based on the non-linear consolidation theory is more accurate and applicable to further load pressures.
Highlights
• Analytical solutions of deformation for clay layer(s) undergoing nonlinear consolidation are derived.
• A type curve method is proposed to determine the hydraulic parameters of clay layers.
• The deformation of the clay layer increases with the increase in the compression index, thickness, and overburden load, and with the decrease in the initial void ratio and initial effective stress.
1 Introduction
With global population growth and economic development, groundwater pollution and geological disasters induced by human activities are becoming increasingly serious, threatening the sustainable development of human society (; ; ; ; ). Clayey soils widely cover the surface of basins, plains, seabed, and groundwater aquifer systems, such as the Yangtze Delta in China and the Virginia coastal plain aquifer system in the U.S. (; ). Clay and silt are the main components of aquitards, and the hydraulic characteristics of aquitards are closely related to groundwater environmental protection, groundwater storage utilization, and land subsidence (; ; ). In addition, clay soil is a significant component of the liner systems installed at hazardous waste disposal sites (e.g., municipal solid waste landfills and nuclear waste disposal sites), thus it has a positive effect on the protection of the water and soil environment around the site, and is an important component of building foundations, dam seepage control systems, mine tailings dams, etc. (; ; ; ; ; ). Therefore, it is important to determine the deformation of the clay layer and its hydraulic parameters for groundwater environmental protection and geological disaster prevention.
Many studies have been conducted on the consolidation and drainage of aquitards (; ). derived an analytical solution to describe the change of excess pore water pressure in an aquitard based on Terzaghi’s one-dimensional consolidation theory. The drawdown was considered to be a constant in an adjacent aquifer and was verified by laboratory tests. proposed a method for estimating the water release from an aquitard during the entire mining history or a limited period, then applied the method to the aquifer system in the Su-Xi-Chang area, but the groundwater level data were insufficient in the study. In addition, proposed a governing equation of drawdown variation to describe the one-dimensional large-strain consolidation behavior of an aquitard without creep effects and derived an analytical solution for the deformation of a large-strain aquitard under a sudden drawdown in an adjacent aquifer.
In addition, many scholars in soil mechanics and hydrogeology have studied the consolidation and deformation of clay layers (; ; ). derived a completely explicit analytical solution for one-dimensional large-strain consolidation in thick and thin soil layers and compared it with the classical small-strain theory. They also extensively studied the non-linear consolidation properties of clays. pointed out that Terzaghi’s one-dimensional linear consolidation theory cannot accurately predict the dissipation rate of pore pressure. proposed a consistent theory for finite-strain consolidation of thick and uniform clay layers under load and found that the traditional consolidation theories may seriously underestimate excess pore water pressure in soft soil layers. Considering the complex non-linearity of materials and the geometric characteristics of the soil consolidation process, derived an analytical solution and obtained the variation law of soil effective stress and void ratio under different conditions through finite element numerical analysis. established a settlement calculation model for soft soil foundations, developed a method to determine the deformation modulus of the soil before and after damage under load, and proposed a non-linear settlement calculation method considering the structural characteristics of soft clay. The effect of temperature and chemical properties of pore water on the seepage and compression of cohesive soils was also analyzed (; ; ).
Moreover, research on soil permeability is also in full swing. and proposed a novel double threshold segmentation algorithm to segment cracks, pores, and grains, and pore-scale variables are defined and extracted from these X-ray CT images to study the geometric characteristics of microstructures of porous geomaterials. , , ; ; successfully predicted the permeability of soil by establishing a model through X-ray CT imaging test, analyzed the internal microstructures of FC and their effects on the compressive strength and permeability, and studied the hydraulic transport properties and fluid flow in single fractures from different projection directions. established a digital quantum mechanism-based neural network (DQNN) to study the permeability using digital porosity, coordination number, and pore network size.
The prediction accuracy of consolidation and deformation processes for an aquitard is obviously affected by the uncertainties and errors in the values of hydraulic properties of the clay layer. Based on laboratory testing data, proposed several methods to calculate the compressive index and precompression stress from soil compressive test data. An analytical method was proposed to determine the hydrogeological parameters using the data of deformation for a large-strain aquitard (). Given that previous analytical solutions were limited by the assumption of a constant head of an overlying unpressurized aquifer, proposed a type curve method for linear drawdown and deformation cases to estimate the hydrogeological parameters of an aquitard in an aquifer system.
Terzaghi’s one-dimensional linear consolidation () problem is a classic problem in soil mechanics and is the most generally used theory for the consolidation of saturated soils. However, saturated clay soils have obvious non-linear trends during consolidation (; ; ). Under the condition that the drawdown of adjacent aquifers is constant, derived an analytical solution to calculate the drawdown and water release of an aquitard undergoing non-linear consolidation described by and proposed a type curve approach to determine the hydraulic parameters of the aquitard. Moreover, many researchers have studied the non-linear consolidation properties of clay by means of analytical solutions or combined analytical solutions with laboratory experiments under the premise that the reduction of permeability is proportional to the reduction of compressibility () (; ; ; ). To the authors’ knowledge, no analytical method has been proposed to determine the hydraulic parameters of clay layers undergoing non-linear consolidation under an overburden load condition. Such a new solution would be useful to analyze the non-linear consolidation phenomenon of clay in practical engineering problems, such as soft soil building foundations and landfill clay liners.
Based on the one-dimensional non-linear consolidation theory developed by (; ), this study first derived the analytical solutions of the excess pore water pressure and the deformation in accordance with Darcy’s law under constant overburden load. Next, the deformation of clay layers with single or double drainage, i.e., singly (SDCL) or doubly draining clay layer (DDCL), was analyzed, and the factors for the deformation were explored by analyzing the sensitivity of multiple parameters. In addition, through a laboratory single and double drainage consolidation test, a type curve approach based on the analytical solutions of non-linear consolidation theory was developed to determine the hydraulic parameters of clay layers and to verify the theoretical research presented in this study.
2 Mathematical model and analytical solutions
2.1 Governing equation for the consolidation of clay layers
The schematic diagram of the one-dimensional non-linear consolidation of a clay layer under an overlying load condition is shown in Figure 1, in which the thickness of the clay layer is . Using the Lagrangian coordinate system, (L) is positive vertically downward, and the origin is located at the top of the clay layer.
FIGURE 1
; gave a governing equation for the non-linear consolidation of large-deformed soil layers (one-dimensional non-linear consolidation theory) governed by the void ratio , while ignoring the self-weight effect of the soil.where is the void ratio (dimensionless), is the initial void ratio (dimensionless), is the hydraulic conductivity (LT-1), which is the function of ; is the effective stress (ML-1T-2), is the time (T), and is the unit weight of water (ML-2T-2).
The one-dimensional non-linear consolidation theory is assumed here (): 1) the water flow in the clay layer obeys Darcy’s law while one-dimensional consolidation of the saturated clay layer occurs; 2) the compression of the clay layer equals the cumulative water release from the clay layer per unit area, because the compressibility of soil particles and pore water is significantly less than that of soil skeleton, and the soil particles and pore water are assumed to be incompressible; 3) and of the clay layer corresponding to the are in accordance with the compression and permeability models (Eqs 2, 3), respectively ():where is the compression index (dimensionless), is the initial effective stress (ML-1T-2), is the penetration index (dimensionless), and is the initial hydraulic conductivity; 4) the clay layer undergoes normal consolidation, and the non-linear variations of the soil compressibility during the consolidation period are assumed to obey Eq. 4 (), and the creep effect is not considered here.where is the final void ratio, which can be calculated by Eq. 2 using the ultimate effective stress at the end of clay layer consolidation, and is equal to the sum of and additional effective stress caused by the overlying load above the clay layer.
Furthermore, the principle of effective stress is observed in the clay layer,where is the load due to some overlying structure or layers (ML-1T-2), and is the excess pore water pressure at position in the clay layer at time (ML-1T-2).
It is assumed that the volumetric compressibility and hydraulic conductivity of the clay layer decrease proportionally (i.e., ) during consolidation and the corresponding consolidation process (; ; ; ). The initial effective stress (, geostatic stress) is distributed uniformly through the entire thickness and has little effect on a thin clay layer, and the calculation error of the water released from a thick clay layer will be greater (; ). Substituting Eqs 2–5 into Eq. 1, the governing equation for the dissipation of excess pore water pressure through a clay layer undergoing one-dimensional non-linear consolidation can be obtained, namely,
Eq. 6 uses excess pore water pressure as the primary variable, andwhich is closely related to the coefficient of consolidation () in Terzaghi’s one-dimensional linear consolidation theory (L2T-1). In addition, relies on the values of , , , and .
2.2 Analytical solutions
Under the conditions of SDCL and DDCL, in the clay layer can be obtained through the governing equation (Eq. 6) with the following initial and boundary conditions (Eqs 8–10). Here, the overlying load is assumed to be constant.
The variations of excess pore water pressure in SDCL and DDCL can be obtained using the variable transformation and characteristic function method, respectively, and expressed as follows:where , the subscripts 1 and 2 of the symbols correspond to the case of single-sided drainage and double-sided drainage, respectively, and the detailed derivations are provided in Appendix A, B, respectively.
The void ratio variation for SDCL and DDCL can be obtained by substituting Eqs 9, 10, respectively, with the assumption proposed by .
Similarly, the hydraulic conductivity variation for SDCL and DDCL can be derived, respectively.
Based on Darcy’s law of seepage, the dimensionless flow velocity for SDCL and DDCL can be obtained, respectively.where is the flow velocity at position and time in SDCL and DDCL (L3T−1), is the dimensionless form of , and is the dimensionless time.
According to the water balance principle, the dimensionless rate of water release for the SDCL and DDCL per unit horizontal area , respectively, where and represent the flow velocity at the upper and lower clay layer interfaces, respectively.
By integrating Eqs 17, 18 over time, the cumulative water fluxes per unit horizontal area of SDCL and DDCL can be obtained, respectively.
Owing to the above assumption about the incompressibility of soil particles, the accumulated released water per unit area of SDCL and DDCL , which is equal to the deformation of the clay layer.where is the deformation of SDCL and DDCL (L). The analytical solutions proposed above are compared with the solution based on one-dimensional linear consolidation theory, where the former considers the non-linear variation of compressibility and permeability during the consolidation of the clay layer (). Therefore, the proposed solution is more practical and accurate in estimating deformation and consolidation problems than the conventional linear theory, especially for newly formed soft sedimentary soils.
2.3 Dimensionless analysis
Figure 2A shows that the evolution of with . is always equal to 0, which leads to for SDCL; while the flow values at the upper and lower clay layer interfaces for DDCL are equal, but in the opposite direction. It is obvious that both SDCL and DDCL have the same variation law: 1) In the early stage of the consolidation process, due to the abrupt overlying load at the initial time, in the clay layer instantaneously reaches the maximum value, and decreases significantly over time. The value of of DDCL is greater than that of SDCL. 2) The value of gradually decreases and tends to 0 in the second stage of the consolidation process, while gradually dissipates. When approaches 0, the consolidation of the clay layer is completed.
FIGURE 2
Figure 2B shows the deformation in SDCL and DDCL. Obviously, the curve also has the same variation characteristics as the curve. The of SDCL is equal to , and the of DDCL is equal to the sum of and . The deformation of the clay layer increases significantly in the early stage of the consolidation process, and then gradually approaches a stable value. At the end of consolidation, the dimensionless deformation of DDCL is the same as that of SDCL. The time to complete consolidation for SDCL is four times longer than that for DDCL.
3 Results and discussion
In order to evaluate the above analytical solutions, the values of , , and of SDCL and DDCL were calculated by case analysis. Then, the values of these parameters based on the one-dimensional linear consolidation theory and the analytical solutions proposed in this study (based on the one-dimensional non-linear consolidation theory) were compared. According to previous studies (; ), the main parameters of the clay layer used in the case study are as follows: = = 0.0532, = 51 kPa, = 1.101, = 2.10 × 10−5 m/day, = 100 kPa, and = 10 m. In the one-dimensional linear consolidation theory, the parameters of the clay layer were considered to be constant, such as and , and the coefficient of volume compressibility () of the clay layer was obtained by .
The values of at different times ( = 1, 3, 5, 10 years) and positions ( = 1, 3, 5, 7, 9 m) were calculated by the one-dimensional linear consolidation theory, and it was calculated by Eqs 11, 12 for SDCL and DDCL, respectively. The results are shown in Figure 3. Obviously, when = 0, is equal to the overlying load . From Figure 3A and Figure 3C, the value of in SDCL decreases with the increase in and . From Figure 3B and Figure 3D, the value of in DDCL decreases with the increase in . Meanwhile, the curve of is symmetric around =5 m, which means that the clay layer in the middle is the greatest. It was also found that the value of predicted by Eqs 11, 12 was greater than that by the linear consolidation theory, which indicated that the dissipation of predicted by the non-linear consolidation theory was lower. The stable time predicted by the two theories for the completion of the consolidation of the clay layer is the same. From Figure 3, is about 21 years for DDCL, which is close to a quarter of that of SDCL. This is because the length of the drainage path in SDCL is , which is approximately twice that of DDCL.
FIGURE 3
The values of at different times and positions in SDCL and DDCL were calculated, and the results are shown in Figure 4. The results show that in SDCL, decreases with the increase in , while in DDCL, it is symmetric around = 5 m. In both SDCL and DDCL, eventually decreases to a stable value while the consolidation of the clay layer completes. The values of predicted by linear and non-linear consolidation theories are 1.076 and 1.056, respectively, for the whole section of the clay layer. Moreover, the value of predicted by non-linear consolidation theory is greater than that predicted by linear consolidation theory. The shapes of the graphs in Figure 4 are similar to those in Figure 3, that is, the variation of predicted by non-linear consolidation theory is similar to that of in the consolidation process of the clay layer.
FIGURE 4
The spatial and temporal distributions of in SDCL and DDCL were calculated by Eqs 15, 16, respectively. During the consolidation process of the clay layer, in SDCL, decreases with the increase in and (Figure 5A; Figure 5C); while in DDCL, has a maximum value in the middle of the clay layer, and the upper and lower parts are symmetrical around = 5 m (Figure 5B; Figure 5D). For both conditions, of the clay layer decreases by 66.2% from the initial value of 7.7 × 10−3 m/year to 2.6 × 10−3 m/year when the consolidation of the clay layer completes. That is, the variation of corresponds to the changes in the void ratio in the consolidation process of the clay layer.
FIGURE 5
The values of in SDCL and DDCL were calculated based on linear and non-linear consolidation theories, and the results are shown in Figure 6. During the consolidation process of the clay layer, the deformation gradually increases to a stable value using both theories. In both SDCL and DDCL, the deformation predicted by each theory ends up being equal. Among them, the deformation predicted by Eqs 23, 24 is 0.116 m, which is significantly smaller than the deformation predicted by the one-dimensional linear consolidation theory (0.216 m). The value of predicted by both theories are the same, 21 years for SDCL and 84 years for DDCL.
FIGURE 6
4 Parametric sensitivity analysis
Under overlying load conditions, the deformation [] of a clay layer depends on the hydraulic [] and soil mechanical parameters [, , , , and ] of the clay layer. Since the values were the same for SDCL and DDCL, we analyzed the parametric sensitivity of for DDCL by varying , , , , , and values taken in the above case study one by one through non-linear and linear consolidation theories.
Figure 7 shows that the value of predicted by Eq. 24 is smaller than that predicted by the one-dimensional linear consolidation theory. The difference in the final stable deformation calculated by linear and non-linear consolidation theories increases with increasing (Figure 7A), (Figure 7D), and (Figure 7F) and decreasing (Figure 7B) and (Figure 7C), and it does not depend on the value of (Figure 7E). The value of has a great effect on the value of the clay layer (Figure 7A). The value of predicted by Eq. 24 increases non-linearly with the increase in , while that predicted by the one-dimensional linear consolidation theory increases linearly with the increase in . When = 0.1, 0.2, 0.3, 0.4, and 0.5, the corresponding values predicted by Eq. 24 are 0.213, 0.421, 0.623, 0.823, and 1.016 m, respectively; while the corresponding values predicted by the one-dimensional linear consolidation theory are 0.405, 0.811, 1.216, 1.621, and 2.027 m, respectively. The predicted value of increases with the increase in . When = 0.1, 0.2, 0.3, 0.4, and 0.5, the predicted values of are 24.7, 54.6, 87.5, 140.4, and 142.1 years, respectively.
FIGURE 7
The impact of on is shown in Figure 7B. From the figure, the values of the clay layer predicted by Eq. 24 and the linear consolidation theory decrease non-linearly with the increase in . When is 0.5, 1, 1.5, and 2, the calculated values by Eq. 24 are 0.159, 0.120, 0.096, and 0.080 m, respectively, and the values from the one-dimensional linear consolidation theory are 0.302, 0.226, 0.181, and 0.151 m, respectively. The values of increase with the increase in . When is 0.5, 1, 1.5, and 2, the obtained values are 33.1, 15.8, 14.4, and 13.2 years in DDCL, respectively.
The effect of on is shown in Figure 7C. When is 10, 20, 50, and 100 kPa, the calculated values by Eq. 24 are 0.250, 0.187, 0.115, and 0.073 m, respectively, and the calculated values by the one-dimensional linear consolidation theory are 1.10, 0.55, 0.22, and 0.11 m, respectively. The value of of the clay layer predicted by Eq.24 decreases non-linearly with the increase in , and that predicted by the one-dimensional linear consolidation theory decreases linearly with the increase in . The predicted decreases with the increase of . When is 10, 20, 50, and 100 kPa, the obtained values are 92.8, 33.1, 13.7, and 6.5 years, respectively.
The impact of on is shown in Figure 7D. The values predicted by both Eq. 24 and the linear theory increase linearly with the increase in . When is 5, 10, 15, and 20 m, the values predicted by Eq. 24 are 0.057, 0.114, 0.171, and 0.228 m, respectively, and that calculated by the one-dimensional linear consolidation theory are 0.108, 0.216, 0.323, and 0.431 m, respectively. The predicted value of increases with the increase in . When is 5, 10, 15, and 20 m, the values are 4.2, 22.8, 33.0, and 67.6 years, respectively.
The impact of on is shown in Figure 7E. The values obtained by Eq. 24 and the linear consolidation theory for all cases are 0.114 and 0.216 m, respectively. The values of decreases linearly with the increase in . When values are 3.15 × 10−3, 6.31× 10−3, 1.58 × 10−2, and 3.15 × 10−2 m/year, the values are 41.9, 20.8, 8.4, and 4.1 years, respectively.
The impact of on is shown in Figure 7F. When is 5, 10, 15, and 20 m, the values predicted by Eq. 24 are 0.072, 0.114, 0.144, and 0.167 m, respectively, and those by the one-dimensional linear consolidation theory are 0.108, 0.216, 0.323, and 0.431 m, respectively. The predicted value of by Eq.24 increases non-linearly with the increase in , and that by the one-dimensional linear consolidation theory increase linearly with the increase in . The value of is not affected by , and has the constant value of about 21 years by both Eq. 24 and the linear consolidation theory.
5 Laboratory experimental studies
In order to validate the analytical solutions proposed in this study, two laboratory experimental studies were performed for SDCL and DDCL. Besides, a type curve approach was developed to determine the clay layer parameters based on the experimental results.
5.1 Determination of clay layer parameters through type curve analysis
According to Eqs 23, 24, the deformation of the clay layer can be expressed as:where is the dimensionless deformation of the clay layer.
By changing Eq. 25 and dimensionless time ( = ) into logarithmic forms, the following equations can be obtained:
Eq.(26) and (27) are both decomposed into modes, where the corresponding variable and constant terms are added. Therefore, in the double logarithmic coordinate system, the shapes of the curves presented by the deformation and the dimensionless deformation are analogous. However, there is an offset in the position of the two curves in the coordinate system ( in the horizontal direction; in the vertical direction), but the order of magnitude of the X and Y-axes corresponds to the same. Therefore, the two curves were translated to enable the data on the two curves to overlap as much as possible, and one overlapping point was selected as the matching point. Substituting the coordinates of the matching point (; , ,; ) into Eqs 28–30, the following parameters can be determined.
5.2 Results of the SDCL experiment
The single-sided drainage experiment carried out by was designed to test the one-dimensional consolidation of the aquitard under load. The thickness of the clay layer was 20 cm, the initial void ratio was 1.28, the load stage was 0–50–100–200 kPa, and the deformations of the clay layer corresponding to the three stages were 6.9, 6.2, and 6.6 mm, respectively.
In this study, the second-level load condition (50–100 kPa) was taken as an example to determine the hydraulic parameters of SDCL with = 50 kPa; = 19.31 cm at this load stage. The deformation of SDCL and the dimensionless deformation curve with time are plotted in Figure 8. The type curve of was then superposed onto this figure, by keeping the axes of the two graphs parallel to each other. The matching point refers to the intersection point of the type curve with the data curve of . The coordinates of the matching point were 6.220 mm; 0.980; 48,303 s,; 9.105 in both systems. Substituting these values into Eqs 28–30, the values of ; ,; were calculated to be 7.03×10−6 m2/s, 0.25, and 6.89×10−8 m/s, respectively. Using the one-dimensional linear consolidation theory; was the same as , volume compressibility () and hydraulic conductivity were calculated to be 0.66 MPa-1 and 4.62×10−8 m/s (), respectively. By the analytical solutions of , the volume compressibility and hydraulic conductivity were calculated to be 0.66 MPa-1 and 4.60×10−8 m/s, respectively. The calculated and values from the analytical solutions of are close to those obtained by linear consolidation theory, and the value predicted by Eq. 29 is larger than that determined by the other two methods.
FIGURE 8
Using the hydraulic parameters determined by the proposed type curve approach, the predicted and measured deformations of SDCL were compared under two load stages (50–100–200 kPa), as shown in Figure 9. From Figure 9, when the load is in the range of 50–100 kPa, the values calculated by the analytical solutions in this study, the one-dimensional linear consolidation theory, and the analytical solutions by are consistent. Under this load stage, the calculated values by the analytical solutions in this study, the one-dimensional linear consolidation theory, and the analytical solutions of are 6.27×10−3 m, 6.35×10−3 m, and 6.22×10−3 m, respectively. In contrast, the measured value of from the experiment is 6.23×10−3 m. During the load stage of 100–200 kPa, the calculated by the analytical solutions in this study, the one-dimensional linear consolidation theory, and the analytical solutions of are 6.50×10−3 m, 1.27×10−2 m, and 1.15×10−2 m, respectively, while the measured value obtained in the laboratory test is 6.57×10−3 m.
FIGURE 9
The value predicted by the analytical solutions in this study agrees well with the experimental results, while those predicted by the one-dimensional linear consolidation theory and the analytical solutions of are significantly larger than the experimental results. Among them, the value calculated by the analytical solutions of is smaller than that by the one-dimensional linear consolidation theory.
5.3 Results of the doubly draining clay layer experiment
The consolidation deformation characteristics of DDCL under overlying load conditions were studied by laboratory tests. The steps of the test were as follows: (a) The inner walls of the consolidation vessel were cleaned. The cohesive soil was air-dried, crushed, screened, milled, and dried. After cooling to room temperature, the soil sample was weighed with an electronic balance and put into a consolidation container (20 mm high). The soil sample was compacted in three layers with the corresponding heights of 6.67mm, 13.33mm, and 20mm, respectively (Figures 10A–C). (b) The sample (Figure 10D) was put into the rock vacuum saturated water test device (DP SJA, institute for Beijing Yaou Depeng Technology Co., Ltd., Beijing, China), the saturated water pressure was set to 20 kPa, the saturated pressure time was set to 300 min, the saturated water time was set to 240 min, and the immersion time was set to 2,440 min. After the soil saturation test was complete, the container of soil samples was removed and the soil samples were taken out (Figure 10E). (c) The double-sided drainage consolidation test was performed on the saturated soil samples. A single-lever arm oedometer (Institute for Nanjing Ningxi Soil Instrument Co., Ltd., Nanjing, China), which met the requirements of the standard consolidation test, was used in the test (Figure 10F). The load in the consolidation vessel was gradually increased through the following pressure sequence (e.g., 50, 100, and 200 kPa). The dial indicator reading was recorded during the test until the deflection stabilized. The initial thickness of the test soil sample was 2 cm, and the initial void ratio was 0.63. In order to avoid the influence of the stress state in the early stage caused by the difference in the previous sample preparation, the first load stage (0–50 kPa) was used as the pre-load stage in this test. Therefore, the characteristics of the clay layer, , , and at the end of the consolidation under the load stage of 0–50 kPa are 1.82 cm, 50 kPa, and 0.54, respectively.
FIGURE 10
The deformation of DDCL and the dimensionless deformation curve with time are plotted in Figure 11. The coordinates of the matching point for DDCL are = 0.632 mm, = 0.980, = 3,600 s, and = 11.218 in both systems. According to the analytical solution proposed in this study, = 1.03×10−6 m2/s, = 1.06×10−8 m/s, and = 0.18. In contrast, according to the one-dimensional linear consolidation theory, was the same as , while and were calculated to be 0.71 MPa-1 and 7.32×10−9 m/s, respectively. By the analytical solutions of , the volume compressibility and hydraulic conductivity were calculated to be 0.71 MPa-1 and 7.3×10−9 m/s, respectively. The values of and calculated by the analytical solutions of are close to those obtained by the one-dimensional linear consolidation theory, and the predicted by the analytical solutions in this study has a greater value than those predicted by the other two methods.
FIGURE 11
Figure 12 compares the deformations of DDCL under two load stages (50–100–200 kPa) by using the hydraulic parameters determined by the type curve approach. Similar to the performance of the SDCL, when the load stage is 50–100 kPa, the calculated values by the analytical solutions in this study, the one-dimensional linear consolidation theory, and the analytical solutions of agree with the experimental results. Under this load stage, the calculated by the analytical solutions in this study, the one-dimensional linear consolidation theory, and the analytical solutions of are 6.03 × 10−3 m, 6.45 × 10−3 m, and 6.32 × 10−3 m, respectively, while the measured value of the experiment is 6.10 × 10−3 m. During the load stage of 100–200 kPa, the value predicted by the analytical solutions in this study agrees well with the experimental results, while the values predicted by the one-dimensional linear consolidation theory and the analytical solutions of are significantly larger than that predicted by the analytical solutions in this study. Under the load stage of 100–200 kPa, the calculated values by the analytical solutions in this study, the one-dimensional linear consolidation theory, and the analytical solutions of are 6.23 × 10−3 m, 1.29 × 10−2 m, and 1.20 × 10−2 m, respectively, while the measured values obtained in the laboratory test is 6.25 × 10−4 m. Therefore, based on these experimental results, the analytical solutions in this study predict the deformation of the clay layer more accurately than the one-dimensional linear consolidation theory and the analytical solutions of .
FIGURE 12
6 Conclusion
In this study, analytical solutions (ignoring creep effects) for SDCL and DDCL undergoing non-linear consolidation, excess pore water pressure, and deformation under overlying load conditions were derived. Combined with the case analysis based on the analytical solution presented in this paper, it is found that under the condition of an overlying load, the deformation in SDCL is equal to that in DDCL. There is a lag in for SDCL compared to DDCL. The drainage path of the consolidated drainage of SDCL is twice as large as that of DDCL, resulting in a lag that is about four times that of DDCL. In addition, there is a difference in under the two theories, i.e., the value by the non-linear consolidation theory is larger than that by the linear consolidation theory. The deformation predicted by the non-linear consolidation theory is smaller than that by the linear consolidation theory. Larger , , and and smaller and lead to a larger difference in the amount of deformation under the two theories. The predicted by the two theories are the same. Through the sensitivity analysis of the parameters, it is obvious that the deformation in SDCL and DDCL is controlled by , , , , and . However, the variation of does not affect the deformation in SDCL and DDCL, and does not depend on the overlying load of the clay layer. Besides, a type curve approach is proposed in this study based on the non-linear consolidation theory, and it agrees well with the measured data from the test. Compared with the linear consolidation theory and the analytical solution of , the analytical solution based on the non-linear consolidation theory in this study is more accurate and suitable for predicting the consolidation and deformation of the clay layer under further load conditions.
7 Open research
Parameter values in Sections 2, 3, and 4 are available through (https://doi.org/10.16285/j.rsm.2010.01.045) and (https://doi.org/10.1002/hyp.11373). The test data in Section 5.2 are available through (https://doi.org/10.1016/j.jhydrol.2022.127994). The test data in Section 5.3 are available through the repository of Dryad (https://datadryad.org/stash/share/frvva1Jd1ETIo-VSViK4-ChOTT6hW9i9EjP2Foj1ExA) and the supporting information files.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
RW is responsible for manuscript preparation and data collection; ZL is responsible for the preparation, data correction and approval of manuscripts, as well as providing fund support; MX, QZ, WI, and HL are responsible for the correction and approval of manuscripts.
Funding
This study was financed by the Open Fund of State Key Laboratory of Geohazard Prevention and Geoenvironment Protection (Chengdu University of Technology) (SKLGP2021K023), Sichuan Provincial Science and Technology Department Youth Fund (2023NSFSC0802), the National Natural Science Foundation of China (41702253) and the China Scholarship Council for working at the University of Waterloo as a visiting scholar. WA. I appreciates the support of the Discovery Grant from the Natural Sciences and Engineering Research Council of Canada (NSERC), which made this collaboration possible.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2023.1131128/full#supplementary-material
References
1
AbdullahW. S.AlshibliK. A.Al-Zou'biM. S. (1999). Influence of pore water chemistry on the swelling behavior of compacted clays. Appl. Clay Sci.15 (5-6), 447–462. 10.1016/S0169-1317(99)00034-4
2
ChenP. S.LiJ. C.HuangM. H.LiD. J. (2021). Consolidation of viscoelastic soil with vertical drains for continuous drainage boundary conditions incorporating a fractional derivative model. Front. Mater.8, 670150. 10.3389/fmats.2021.670150
3
DavisE. H.RaymondG. P. (1965). A non-linear theory of consolidation. Geotechnique15 (2), 161–173. 10.1680/geot.1965.15.2.161
4
DengY. F.YueX. B.CuiY. J.ShaoG. H.LiuS. Y.ZhangD. W. (2014). Effect of pore water chemistry on the hydro-mechanical behaviour of Lianyungang soft marine clay. Appl. Clay Sci.95, 167–175. 10.1016/j.clay.2014.04.007
5
GibsonR. E.EnglandG. L.HusseyM. J. L. (1967). The theory of one-dimensional consolidation of saturated clays. Geotechnique17 (3), 261–273. 10.1680/geot.1967.17.3.261
6
GibsonR. E.SchiffmanR. L.CargillK. W. (1981). The theory of one‐dimensional consolidation of saturated clays. II. Finite nonlinear consolidation of thick homogeneous layers. Geotechnique18, 280–293. 10.1139/t81-030
7
GorelickS. M.ZhengC. (2015). Global change and the groundwater management challenge. Water Resour. Res.51 (5), 3031–3051. 10.1002/2014wr016825
8
GregoryA. S.WhalleyW. R.WattsC. W.BirdN. R. A.HallettP. D.WhitmoreA. P. (2006). Calculation of the compression index and precompression stress from soil compression test data. Soil Tillage Res.89 (1), 45–57. 10.1016/j.still.2005.06.012
9
GuoH. P.HaoA. B.LiW. P.ZangX. S.WangY. L.ZhuJ. Y.et al (2022). Land subsidence and its affecting factors in cangzhou, north China plain. Front. Environ. Sci.10. 10.3389/fenvs.2022.1053362
10
GuoQ. N.LiH. L. (2015). Terrestrial-originated submarine groundwater discharge through deep multilayered aquifer systems beneath the seafloor. Hydrol. Process.29 (4), 295–309. 10.1002/hyp.10163
11
JiaX. Y.HouD. Y.WangL. W.O'CommorD.LuoJ. (2020). The development of groundwater research in the past 40 years: A burgeoning trend in groundwater depletion and sustainable management. J. Hydrology587, 125006. 10.1016/j.jhydrol.2020.125006
12
KonikowL. F.NeuzilC. E. (2007). A method to estimate groundwater depletion from confining layers. Water Resour. Res.43 (7), W07417. 10.1029/2006wr005597
13
LiZ. F.ZhouZ. F.ChenZ.LiuG. Q.ZhouC. Y. (2016). An analytical method to estimate groundwater depletion from a confining layer. Nat. Hazards85 (2), 887–901. 10.1007/s11069-016-2609-1
14
LiZ. F.ZhouZ. F.DaiY. F.DaiB. B. (2019). Contaminant transport in a largely-deformed aquitard affected by delayed drainage. J. Contam. Hydrology221, 118–126. 10.1016/j.jconhyd.2019.02.002
15
LiZ. F.ZhouZ. F.DaiY. F.HuangY.GuoQ. N.ZhouC. Y. (2017). Contaminant transport in a small deformation aquitard affected by the delayed drainage phenomenon. Hydrol. Process.31 (24), 4429–4439. 10.1002/hyp.11373
16
LiZ. F.ZhouZ. F.LiM. Y.ZhangB. R.DaiB. B. (2018). Delayed drainage of a largely deformed aquitard due to abrupt water head decline in adjacent aquifer. Geofluids2018, 1–12. 10.1155/2018/2326491
17
LiuC.ShiB.GuK.ZhangT. S.TangC. S.WangY.et al (2021). Negative pore water pressure in aquitard enhances land subsidence: Field, laboratory, and numerical evidence. Water Resour. Res.58 (1), e2021WR030085. 10.1029/2021wr030085
18
LuoH.LiZ. F.LuoM. L.ZhangQ.IllmanW. A. (2020). An analytical method to calculate groundwater released from an aquitard undergoing nonlinear consolidation. Water Resour. Res.56 (9), e2020WR027320. 10.1029/2020wr027320
19
MaB. H.LiZ.CaiK.HuZ. Y.ZhaoM. H.HeC. B.et al (2021). An improved nonlinear settlement calculation method for soft clay considering structural characteristics. Geofluids2021, 1–7. 10.1155/2021/8837889
20
MitchellJ. K. (1993). Fundamentals of soil behavior. New York: Wiley.
21
NeuzilC. E. (1986). Groundwater flow in low‐permeability environments. Water Resour. Res.22 (8), 1163–1195. 10.1029/wr022i008p01163
22
NichollsR. J.BeavenR. P.StringfellowA.MonfortD.Le CozannetG.WahlT.et al (2021). Coastal landfills and rising sea levels: A challenge for the 21st century. Front. Mar. Sci.8. 10.3389/fmars.2021.710342
23
RomeroE.GensA.LloretA. (2001). Temperature effects on the hydraulic behaviour of an unsaturated clay. Geotechnical Geol. Eng.19, 311–332. 10.1023/A:1013133809333
24
ShuB.ChenY.Amani-BeniM.ZhangR. Z. (2022). Spatial distribution and influencing factors of mountainous geological disasters in southwest China: A fine-scale multi-type assessment. Front. Environ. Sci.10, 1049333. 10.3389/fenvs.2022.1049333
25
TerzaghiK. (1943). Theoretical soil mechanics. New York: Wiley. 10.1002/9780470172766
26
ToufighM. M.OuriaA. (2009). Consolidation of inelastic clays under rectangular cyclic loading. Soil Dyn. Earthq. Eng.29 (2), 356–363. 10.1016/j.soildyn.2008.03.006
27
WangF. Y.HuangH. W.YinZ. Y.HuangQ. (2021). Probabilistic characteristics analysis for the time-dependent deformation of clay soils due to spatial variability. Eur. J. Environ. Civ. Eng.26, 6096–6114. 10.1080/19648189.2021.1933604
28
WuJ.XieX. Y.ZhuX. R. (2010). Study of properties of 1-D complex nonlinear consolidation of saturated soils (in Chinese). Rock Soil Mech.31 (1), 81–86. 10.16285/j.rsm.2010.01.045
29
XieK. H.LeoC. J. (2004). Analytical solutions of one-dimensional large strain consolidation of saturated and homogeneous clays. Comput. Geotechnics31 (4), 301–314. 10.1016/j.compgeo.2004.02.006
30
XieK. H.XieX. Y.JiangW. (2002). A study on one-dimensional nonlinear consolidation of double-layered soil. Comput. Geotechnics29 (2), 151–168. 10.1016/s0266-352x(01)00017-9
31
XueY. Q.ZhangY.YeS. J.WuJ. C.LiQ. F. (2005). Land subsidence in China. Environ. Geol.48 (6), 713–720. 10.1007/s00254-005-0010-6
32
YanH. X.XieH. J.WuJ. W.DingH.QiuZ. H.SunZ. L. (2021). Analytical model for transient coupled consolidation and contaminant transport in landfill liner system. Comput. Geotechnics138, 104345. 10.1016/j.compgeo.2021.104345
33
ZhangY. H.HuoX. F.LuoY. (2023). Prediction of groundwater pollution diffusion path based on multi-source data fusion. Front. Environ. Sci.10, 1116309. 10.3389/fenvs.2022.1116309
34
ZhaoZ.ZhouX. P.ChenJ. W. (2022). Pore-scale hydraulic properties of virtual sandstone microstructures: Spatial variations and voxel scale effects. Archives Civ. Mech. Eng.23 (1), 22. 10.1007/s43452-022-00566-7
35
ZhaoZ.ZhouX. P. (2022). Digital analysis for pore-scale compressive strength and permeability of foamed cement with realistic microstructures by X-ray-μCT imaging. Constr. Build. Mater.346, 128456. 10.1016/j.conbuildmat.2022.128456
36
ZhaoZ.ZhouX. P. (2020a). Digital voxel-based fracture extraction: Insights to characterization of single fracture flow and anisotropy permeability. J. Nat. Gas Sci. Eng.84, 103635. 10.1016/j.jngse.2020.103635
37
ZhaoZ.ZhouX. P. (2020b). Pore-scale effect on the hydrate variation and flow behaviors in microstructures using X-ray CT imaging. J. Hydrology584, 124678. 10.1016/j.jhydrol.2020.124678
38
ZhaoZ.ZhouX. P.QianQ. H. (2021). Dqnn: Pore-scale variables-based digital permeability assessment of carbonates using quantum mechanism-based machine-learning. Sci. China Technol. Sci.65 (2), 458–469. 10.1007/s11431-021-1906-1
39
ZhaoZ.ZhouX. P.QianQ. H. (2020). Fracture characterization and permeability prediction by pore scale variables extracted from X-ray CT images of porous geomaterials. Sci. China Technol. Sci.63 (5), 755–767. 10.1007/s11431-019-1449-4
40
ZhouX. P.ZhaoZ. (2019). Digital evaluation of nanoscale-pore shale fractal dimension with microstructural insights into shale permeability. J. Nat. Gas Sci. Eng.75, 103137. 10.1016/j.jngse.2019.103137
41
ZhouX. P.ZhaoZ.LiZ. (2020). Cracking behaviors and hydraulic properties evaluation based on fractural microstructure models in geomaterials. Int. J. Rock Mech. Min. Sci.130, 104304. 10.1016/j.ijrmms.2020.104304
42
ZhouZ. F.GuoQ. N.DouZ. (2013). Delayed drainage of aquitard in response to sudden change in groundwater level in adjacent confined aquifer: Analytical and experimental studies. Chin. Sci. Bull.58 (25), 3060–3069. 10.1007/s11434-013-5730-5
43
ZhuW.YangX.HeJ.WangX.LuR.ZhangZ. (2021). Investigation and systematic risk assessment in a typical contaminated site of hazardous waste treatment and disposal. Front. Public Health9, 764788. 10.3389/fpubh.2021.764788
44
ZhuangC.YanL.ZhouZ. F.WangJ. G.DouZ. (2021). Estimation of aquitard hydraulic conductivity and skeletal specific storage considering non-Darcy flow. Water Sci. Eng.14 (4), 269–276. 10.1016/j.wse.2021.09.003
45
ZhuangC.ZhouZ. F.IllmanW. A.DouZ.WangJ. G. (2020). Parameter estimation of an overconsolidated aquitard subjected to periodic hydraulic head variations within adjacent aquifers. J. Hydrology583, 124555. 10.1016/j.jhydrol.2020.124555
46
ZhuangC.ZhouZ. F.ZhanH. B.WangG. Y. (2015). A new type curve method for estimating aquitard hydraulic parameters in a multi-layered aquifer system. J. Hydrology527, 212–220. 10.1016/j.jhydrol.2015.04.062
47
ZhuoY.WangJ. G.ZhengS. Y.LiW. J. (2022). Experimental study of one-dimensional nonlinear consolidation of an aquitard under multistage loading. J. Hydrology611, 127994. 10.1016/j.jhydrol.2022.127994
Appendix A Derivation of Eq.11
Using the substitution method:
Substituting Eq. A-1 into –Eq.10, can simplify the original equation to as the independent variable:
By applying the boundary subordination processing to Eq. A-5:
The governing equations and boundary conditions for as the independent variable can be obtained as:
Using the separation of variables method, the solution to Eq. A-8–Eq. A-11 can be derived as:
Combining Eq. A-7 with Eq. A-12, the solution to Eq. A-2 can be given,
Finally, the solution to Eq. 11 is obtained by combining Eq. A-1 and Eq. A-13.
Appendix B Derivation of Eq. 12
As the variable passing method in Appendix A, the independent variable is replaced by .
The corresponding solution procedure is the same as given in Appendix A, and the solution to Eq. B-1–Eq. B-4 can be obtained as:
Finally, the solution to Eq. 12 is obtained by combining Eq. A-1 and Eq. B-5.
Summary
Keywords
non-linear consolidation clay layer, analytical solution, deformation, hydraulic parameters, type curve method
Citation
Wang R, Li Z, Xu M, Zhang Q, Illman WA and Li H (2023) Determination of hydraulic parameters of non-linear consolidation clay layers by type curve method. Front. Earth Sci. 11:1131128. doi: 10.3389/feart.2023.1131128
Received
24 December 2022
Accepted
13 February 2023
Published
07 March 2023
Volume
11 - 2023
Edited by
Xiaoping Zhou, Chongqing University, China
Reviewed by
Nan Xiao, Changsha University of Science and Technology, China
Jinguo Wang, Hohai University, China
Yundong Shou, Wuhan University, China
Updates
Copyright
© 2023 Wang, Li, Xu, Zhang, Illman and Li.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhaofeng Li, lizhaofeng17@cdut.edu.cn
This article was submitted to Geohazards and Georisks, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.