**Pre-conference publication**
22 Apr 2020

**Pre-conference publication** | 22 Apr 2020

# Influence of pore water pressure change on consolidation behavior of saturated poroelastic medium

Chao-Lung Yeh Wei-Cheng Lo Cheng-Wei Lin and Chung-Feng Ding

^{1},

^{2},

^{1},

^{1}

**Chao-Lung Yeh et al.**Chao-Lung Yeh Wei-Cheng Lo Cheng-Wei Lin and Chung-Feng Ding

^{1},

^{2},

^{1},

^{1}

^{1}Tainan Hydraulics Laboratory, National Cheng Kung University, Tainan 701, Taiwan^{2}Department of Hydraulic and Ocean Engineering, National Cheng Kung University, Tainan 701, Taiwan

^{1}Tainan Hydraulics Laboratory, National Cheng Kung University, Tainan 701, Taiwan^{2}Department of Hydraulic and Ocean Engineering, National Cheng Kung University, Tainan 701, Taiwan

**Correspondence**: Chao-Lung Yeh (clyeh@thl.ncku.edu.tw)

**Correspondence**: Chao-Lung Yeh (clyeh@thl.ncku.edu.tw)

There are many factors causing land subsidence, and groundwater extraction is one of the most important causes of subsidence. A set of coupled partial differential equations are derived in this study by using the poro-elasticity theory and linear stress-strain constitutive relation to describe the one-dimensional consolidation in a saturated porous medium subjected to pore water pressure change due to groundwater table depression. Simultaneously, the closed-form analytical solutions for excess pore water pressure and total settlement are obtained. To illustrate the consolidation behavior of the poroelastic medium, the saturated layer of clay sandwiched between two sand layers is simulated, and the dimensionless pore water pressure changes with depths and the dimensionless total settlement as function of time in the clay layer are examined. The results show that the greater the water level change in the upper and lower sand layers, the greater the pore water pressure change and the total settlement of the clay layer, and the more time it takes to reach the steady state. If the amount of groundwater replenishment is increased, the soil layer will rebound.

The consolidation of soil induced by excessive withdrawal of groundwater is a worldwide problem and usually known as land subsidence which has occurred in many cities, for example, Bangkok, Venice, Mexico, Shanghai and so on. Therefore, analyzing and predicting soil consolidation due to groundwater withdrawal in a multiaquifer system is important for achieving the control of land subsidence and the sustainable use of groundwater resource. The most widely used theories of consolidation behavior are Terzaghi (1925) and Biot (1941). The study of consolidation was firstly proposed by Terzaghi (1925) considering the unidirectional consolidation rate of saturated clay layers. Despite the significant advance of this theory in capturing the essence of consolidation process, the coupling between fluid flow and deformation in a fully-saturated porous medium is not truly laid on a firm theoretical base. This conceptual breakthrough was achieved later by Biot (1941), who brought the role of the solid and fluid constituents on equal footing to formulate a pair of coupled equilibrium equations of motion using the displacement vector of solid and fluid pore pressure as dependent variables, now known as poroelasticity. Tuncay and Corapcioglu (1996) used volume-average theory to simulate the consolidation problem of poroelastic media saturated by two immiscible Newtonian fluids. In addition, Luo and Zeng (2011) constructed a coupled three-dimensional model of consolidation and groundwater extraction based on the viscoelastic-plastic constitutive relation and rheological theory of Biot's poroelasticity, and used finite element method to simulate and compare groundwater level changes and soil compression. Lo et al. (2014) have recently developed the poroelasticity theory of three-dimensional consolidation in unsaturated soils, which features the displacement vector of the solid phase together with the excess pore water and air pressures as dependent variables.

Tsai et al. (2006) and Tseng et al. (2008) applied the linear poro-elasticity theory to investigate the soil consolidation of multi-aquifer system caused by groundwater level decline due to pumping, but they didn't consider the viscous coupling caused by viscous drag. Considering the coupling between fluid flow and deformation, we analyze the soil consolidation due to groundwater table depression in a multiaquifer system based on the coupled poroelasticity equations formulated by Lo et al. (2014). After the closed-form analytical solutions for excess pore water pressure together with total settlement are obtained, numerical studies are undertaken for highly permeability porous sand (aquifers) surrounding an impervious clay (aquitard).

According to the poroelasticity theory of consolidation developed by Lo et al. (2014), the momentum balance equations for a deformable homogeneous porous medium bearing one compressible, viscous fluid, in the absence of body force and inertial force, can be expressed as

where the subscripts “s” and “f” represent the solid phase and fluid phase,
respectively; ∅ represents the porosity; *p*_{f}
indicates the gauge pressure (excess pressure); *u*_{f} and
*u*_{s} designate the displacements of each phase;
${R}_{\mathrm{22}}(=-{\mathrm{\varnothing}}^{\mathrm{2}}{\mathit{\mu}}_{\mathrm{f}}/{k}_{\mathrm{s}})$ signifies the
coefficient of viscous drag due to the velocity difference between solid
phase and fluid phase, where *μ*_{f} is the dynamic shear viscosity of
the fluid phase, *k*_{s} is the intrinsic permeability. $\stackrel{\mathrm{\u203e}}{\stackrel{\mathrm{\u203e}}{\mathit{\sigma}}}$ denotes the total stress tensor for the porous medium, using the
sign convention that tension is positive. On the basis of Biot's theory
(1962), the total stress components $\stackrel{\mathrm{\u203e}}{\stackrel{\mathrm{\u203e}}{\mathit{\sigma}}}$ of the
bulk material is defined as

where ${\stackrel{\mathrm{\u203e}}{\stackrel{\mathrm{\u203e}}{t}}}_{\mathrm{s}}$ is the stress tensor for the solid phase, $\stackrel{\mathrm{\u203e}}{\stackrel{\mathrm{\u203e}}{\mathit{\delta}}}$ is the unit tensor.

## 2.1 Linear stress-strain relations

According to the mass balance equation, the linear stress-strain relations for a porous medium containing single fluid

where $\stackrel{\mathrm{\u203e}}{\stackrel{\mathrm{\u203e}}{e}}=\frac{\mathrm{1}}{\mathrm{2}}(\mathbf{\nabla}{\mathit{u}}_{\mathrm{s}}+\mathbf{\nabla}{\mathit{u}}_{\mathrm{s}}^{\mathrm{T}})$ is the solid
phase strain tensor, the superscript “T” representing the transpose; *G* is
the shear modulus of porous medium frame; $e=\mathbf{\nabla}\cdot {\mathit{u}}_{\mathrm{s}}$ is the dilatation of solid phase;
$\mathit{\epsilon}=\mathbf{\nabla}\cdot {\mathit{u}}_{\mathrm{f}}$ is the
dilatation of fluid phase; The elastic coefficients *A*, *Q*, and *R* can be
expressed as function of four elastic moduli (see the detail from Lo et al., 2005). After introduction of volume-averaged form and the linear
stress-strain relation, Eq. (2) becomes

According to Eq. (3b), the fluid dilatation *ε* can be expressed
in terms of the fluid pressure *p*_{f} and the solid dilatation *e* as

In Eq. (5), the parameter *d*_{1} is dimensionless, while
*d*_{2} has the unit of inverse stress. It is conventional to use
stress variables as dependent variables rather than strain variables.
Therefore, we take divergence of both sides of Eq. (1a), and then
substitute Eq. (5) into the result and obtain:

In the same spirit, the linear stress-strain relation can be written as:

For a one-dimensional problem with a time-invariant total compaction stress, Eqs. (6) and (7) can be reduced to

where *w* is the component of the displacement of vector of the solid phase
in the *z* direction. Equation (8b) can be rearranged as

where ${K}^{\left(u\right)}=[\mathrm{2}G+A+Q+(Q+R\left){d}_{\mathrm{1}}\right]$. Thus, the term $\frac{\partial w}{\partial z}$ in Eq. (8a) can be replaced by Eq. (9) to yield the following coupled partial differential equation

where the parameters *q*_{1} and *b*_{1} are defined as

Equations (10) represent a one-dimensional poroelastic problem of consolidation caused by pore water pressure changes. Coupling occurs only in the first-order time-derivative term of the pore water pressure. The diffusivity coefficient ${c}_{v}={b}_{\mathrm{1}}/{q}_{\mathrm{1}}$ denotes the coefficient of consolidation.

## 2.2 Boundary and initial conditions

A stratum of clay sandwiched between sandy strata that are highly permeable
and much stiffer than the clay is shown in Fig. 1, water level depression
*h*_{1} and *h*_{2} due to pumping take place in the sandy
strata above and below the clay. Due to significant differences of
permeability and compressibility between sand and clay, excessive pore
pressure only exists in the clay as consolidation proceeds. And nearly all
of the consolidation happens due to the volume change within the clay, while
the sandy strata may be considered rigid media as compared with the clay.
Furthermore, the steady-state solutions, which are usually applied in
engineering practice, will be considered.

The initial conditions can be expressed as ${p}_{\mathrm{f}}(z,\phantom{\rule{0.125em}{0ex}}\mathrm{0})=\mathrm{0}$. With respect to the boundary conditions, the top and bottom surfaces of the porous layers are considered to be permeable so that the pore fluids can drain at both sides of the sample, which leads to ${p}_{\mathrm{f}}\left(\mathrm{0},\phantom{\rule{0.125em}{0ex}}t\right)=-{\mathit{\rho}}_{\mathrm{f}}g{h}_{\mathrm{2}}$ and ${p}_{\mathrm{f}}\left(L,\phantom{\rule{0.125em}{0ex}}t\right)=-{\mathit{\rho}}_{\mathrm{f}}g{h}_{\mathrm{1}}$. Based on Tsai et al. (2006), combined flow equation and Darcy's law and assuming an isotropic and homogeneous porous media and one-dimensional steady-state consolidation, the pore water pressure in the clay stratum can be represented as function of depth as ${p}_{\mathrm{f}}=-{\mathit{\rho}}_{\mathrm{f}}g[{h}_{\mathrm{2}}+z({h}_{\mathrm{1}}-{h}_{\mathrm{2}}\phantom{\rule{0.125em}{0ex}})/L]$. By combining of initial and boundary equations into Eq. (10), we can get

The time-dependent total settling *s*(*t*) can be
consequently deduced

Equation (13) implies that, as for *t*→∞, the final
total settlement Δ*s* is equal to

In order to investigate quantitatively the deformation of porous media due to water table depression, the governing equations have been implemented in MATLAB and a textbook case is simulated. The ground parameters used in the study are summarized in Table 1, which also lists the numerical value of elastic coefficients and hydraulic parameters necessary for the numerical simulations performed in the present study. Before the consolidation of soil subjected to water level depression is examined, the implementation of the poroelastic consolidation is verified. Considering water level depressions ${h}_{\mathrm{1}}/L={h}_{\mathrm{2}}/L=\mathrm{0.01}$ and elapsed dimensionless time $T={c}_{v}t/{L}^{\mathrm{2}}$, the excess dimensionless pore water pressure ${P}^{\ast}(={p}_{\mathrm{f}}(zt)/({\mathit{\rho}}_{\mathrm{f}}gL\left)\right)$ with respect to the dimensionless depth ${Z}^{\ast}(=z/L)$ are shown in Fig. 2. The pore water pressure increases with time at the inside of the clay stratum. Because the boundary of clay stratum is subjected to a constant water table depression, the dimensionless pore water pressure at top and bottom of clay stratum remains equal to the value caused by the water table depression. The dimensionless pore water pressure at different depth for various dimensionless elapsed times is also plotted in Fig. 3. The dimensionless pore water pressure near the boundary is greater than that near the center of stratum. The change of pore water pressure in the clay layer is directly related to pressures at the upper and lower boundaries resulting from pumping, and the change of pore water pressure is from the boundary to the medium of the layer.

Figure 4 shows the dimensionless settlement ${s}^{\ast}(=\mathrm{s}(t)/\mathrm{\Delta}s)$ with respect to the dimensionless elapsed time. We can easily calculate the quantity of settlement at any time from Fig. 4. The dissipation of pore water pressure is significantly influenced by permeability (Fig. 5). When the intrinsic permeability is higher, the pore water pressure reaching stabilize is faster.

Considering the current situation, the groundwater table may be a time series change. The upper water level changes with time was assumed as ${h}_{\mathrm{1}}=\mathrm{0.1}\times \mathrm{sin}\left(\frac{t\mathit{\pi}}{\mathrm{12}}\right)m$ (the time interval is 1 h, so the water level changes for 12 h as a cycle), lower water level changes do not change with time ${h}_{\mathrm{2}}/L=\mathrm{0.02}$. It can be clearly seen from the results in Fig. 6 that because of the greater fluctuation of the fixed water level at the lower boundary, the trend of the consolidation with time is similar to the above-mentioned fixed water level fluctuation, but at the same time, it is affected by the fluctuation of the upper water level, and the periodic fluctuation appears. When the upper groundwater table raises, the settlement significantly decreases, which is regarded as the phenomenon of soil swelling. It can also be seen that the amount of rebound is less than that of compression. Therefore, increasing the groundwater replenishment from land subsidence prevention, the soil layer will partially rebound.

A set of coupled partial differential equations are presented in this study by using the poro-elasticity theory and linear stress-strain constitutive relation to describe the one-dimensional consolidation in a saturated porous medium subjected to pore water pressure change due to groundwater table depression. Simultaneously, the closed-form analytical solutions for excess pore water pressure together with total settlement are written. To quantitatively analyze the consolidation behavior of the poroelastic medium, a saturated layer of clay sandwiched between two sand layers is taken as example, and the dimensionless pore water pressure changes with depth and the dimensionless total settlement as a function of time in the clay layer are examined. The results show that the greater the water level change in the upper and lower sand layers, the greater the pore water pressure change and the total settlement of the clay layer are, and the more time it takes to reach the steady state. The change of pore water pressure in the clay layer is directly related to the pumping in the upper and lower boundaries, and the change of pore water pressure is from the boundary to the medium of the layer. When the groundwater table rises again, the phenomenon of soil swelling appears. If the amount of groundwater replenishment is increased, the soil layer will rebound partially.

The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.

CLY designed a study concept and made a general text writing and editing, WCL help to develop the concept of the approach. CWL and CFD provided assistance for data analyses.

The authors declare that they have no conflict of interest.

This article is part of the special issue “TISOLS: the Tenth International Symposium On Land Subsidence – living with subsidence”. It is a result of the Tenth International Symposium on Land Subsidence, Delft, the Netherlands, 17–21 May 2021.

Biot, M. A.: General theory of three-dimensional consolidation, J. Appl. Phys., 12, 155–164, 1941.

Lo, W. C., Sposito, G., and Majer, E.: Wave propagation through elastic porous media containing two immiscible fluids, Water Resour. Res., 41, W02025, https://doi.org/10.1029/2004WR003162, 2005.

Lo, W. C., Yeh, C. L., and Tsai, C. T.: Effect of soil texture on the propagation and attenuation of acoustic wave at unsaturated conditions, J. Hydrol., 338, 273–284, 2007.

Lo, W. C., Sposito, G., and Chu, H.: Poroelastic theory of consolidation in unsaturated soils, Vadose Zone J., 13, 5, https://doi.org/10.2136/vzj2013.07.0117, 2014.

Luo, Z. J. and Zeng, F.: Finite element numerical simulation of land subsidence and groundwater exploitation based on visco-elastic-plastic Biot's consolidation theory, J. Hydrodynam., 23, 615–624, 2011.

Rawls, W. J., Ahuja, J. R., and Brakensiek, D. L.: Estimating soil hydraulic properties from soils data. Proceedings of Workshop on Indirect Methods for Estimating the Hydraulic Properties of Unsaturated Soils, Riverside, CA, 329–341, 1992.

Terzaghi, K.: Erdbaumechanik auf bodenphysikalisher Grundlage: Vienna, Austria, Deuticke, 399 pp., 1925.

Tsai, T. L., Chang, K. C., and Huang, L. H.: Body force effect on consolidation of porous elastic media due to pumping, J. Chin. Inst. Eng., 29, 75–82, 2006.

Tseng, C. M., Tsai, T. L., and Huang, L. H.: Effects of body force on transient poroelastic consolidation due to groundwater pumping, Environ. Geol., 54, 1507–1516, 2008.

Tuncay, K. and Corapcioglu, M. Y.: Consolidation of elastic porous media saturated by two immiscible fluids, J. Eng. Mech., 122, 1077–1085, 1996.