{"id":276,"date":"2022-01-13T23:17:39","date_gmt":"2022-01-13T23:17:39","guid":{"rendered":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/chapter\/box-3-3d-poroelasticity-equations\/"},"modified":"2022-01-17T19:09:16","modified_gmt":"2022-01-17T19:09:16","slug":"box-3-3d-poroelasticity-equations","status":"publish","type":"chapter","link":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/chapter\/box-3-3d-poroelasticity-equations\/","title":{"raw":"Box 3 3\u2011D Poroelasticity Equations","rendered":"Box 3 3\u2011D Poroelasticity Equations"},"content":{"raw":"
\r\n

Theoretically, land subsidence is best analyzed according to the theory of consolidation (Biot, 1941), which holds that consolidation itself represents the response of a compressible porous medium to changes in the flow field operating within it. A complete analysis of land subsidence requires determination of the 3\u2011D deformation field accompanying the 3\u2011D flow field, and must be accomplished in a complex multi\u2011aquifer system. A few basic principles underlie the consolidation process. As outlined above, the first principle advanced by Terzaghi (1923) states that the total stress \u03c3<\/em>tot<\/em><\/sub> at any point of the porous medium is equal to the sum of the effective intergranular \u03c3<\/em>eff<\/em><\/sub> and the neutral pore pressure p<\/em>:<\/p>\r\n

\u03c3<\/em>tot<\/em><\/sub> = \u03c3<\/em>eff<\/em><\/sub> + p<\/em><\/p>\r\n

Deformation of the porous body is controlled exclusively by variation of the effective stress \u03c3<\/em>eff<\/em><\/sub>. If we consider changes relative to an initial undisturbed state of equilibrium, the Cauchy equations of equilibrium are cast in terms of incremental effective stress and pore pressure as shown in Equation\u00a0Box\u00a03\u20111.<\/p>\r\n\r\n\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\frac{\\partial \\sigma _{xx}}{\\partial x}+\\frac{\\partial \\tau _{xy}}{\\partial y}+\\frac{\\partial \\tau _{xz}}{\\partial z}=\\frac{\\partial p}{\\partial x}[\/latex]<\/td>\r\n<\/td>\r\n<\/tr>\r\n
<\/td>\r\n[latex]\\displaystyle \\frac{\\partial \\tau _{yx}}{\\partial x}+\\frac{\\partial \\sigma _{yy}}{\\partial y}+\\frac{\\partial \\tau _{yz}}{\\partial z}=\\frac{\\partial p}{\\partial y}[\/latex]<\/td>\r\n(Box 3-1)<\/td>\r\n<\/tr>\r\n
<\/td>\r\n[latex]\\displaystyle \\frac{\\partial \\tau _{zx}}{\\partial x}+\\frac{\\partial \\tau _{zy}}{\\partial y}+\\frac{\\partial \\sigma _{zz}}{\\partial z}=\\frac{\\partial p}{\\partial z}[\/latex]<\/td>\r\n<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nwhere:\r\n\r\n\r\n\r\n\r\n\r\n
\u03c3<\/em>xx<\/em><\/sub><\/td>\r\n=<\/td>\r\nincremental normal effective stress in x direction<\/td>\r\n<\/tr>\r\n
\u03c3<\/em>yy<\/em><\/sub><\/td>\r\n=<\/td>\r\nincremental normal effective stress in y direction<\/td>\r\n<\/tr>\r\n
\u03c3<\/em>zz<\/em><\/sub><\/td>\r\n=<\/td>\r\nincremental normal effective stress in z direction<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nand, the incremental shear stresses:\r\n\r\n\r\n\r\n\r\n\r\n
\u03c4<\/em>xy<\/em><\/sub><\/td>\r\n=<\/td>\r\n\u03c4<\/em>yx<\/em><\/sub><\/td>\r\n<\/tr>\r\n
\u03c4<\/em>xz<\/em><\/sub><\/td>\r\n=<\/td>\r\n\u03c4<\/em>zx<\/em><\/sub><\/td>\r\n<\/tr>\r\n
\u03c4<\/em>yz<\/em><\/sub><\/td>\r\n=<\/td>\r\n\u03c4<\/em>zy<\/em><\/sub><\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

The relationships between the incremental effective stress tensor \u03c3<\/em> and the incremental strain tensor \u03b5<\/em> for a geomechanical isotropic medium are shown in Equation\u00a0Box\u00a03\u20112.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\begin{bmatrix}\r\n\\sigma _{xx} \\\\\r\n\\sigma _{yy} \\\\\r\n\\sigma _{zz} \\\\\r\n\\tau _{xy} \\\\\r\n\\tau _{xz} \\\\\r\n\\tau _{yz}\\end{bmatrix}=D^{-1}\\begin{bmatrix}\r\n\\varepsilon _{xx} \\\\\r\n\\varepsilon _{yy} \\\\\r\n\\varepsilon _{zz} \\\\\r\n\\varepsilon _{xy} \\\\\r\n\\varepsilon _{xz} \\\\\r\n\\varepsilon _{yz}\\end{bmatrix}[\/latex]<\/td>\r\n(Box 3-2)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

with matrix D<\/em>-1<\/sup> as shown in Equation\u00a0Box\u00a03\u20113.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle D^{-1}=\\frac{E}{(1+\\nu )(1-2\\nu )} \\begin{bmatrix}\r\n1-\\nu & \\nu & \\nu & 0 & 0 & 0 \\\\\r\n\\nu & 1-\\nu & \\nu & 0 & 0 & 0 \\\\\r\n\\nu & \\nu & 1-\\nu & 0 & 0 & 0 \\\\\r\n0 & 0 & 0 & \\frac{1-2\\nu }{2} & 0 & 0 \\\\\r\n0 & 0 & 0 & 0 & \\frac{1-2\\nu }{2} & 0 \\\\\r\n0 & 0 & 0 & 0 & 0 & \\frac{1-2\\nu }{2} \\\\\r\n\\end{bmatrix}[\/latex]<\/td>\r\n(Box 3-3)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nwhere:\r\n\r\n\r\n\r\n\r\n
E<\/em><\/td>\r\n=<\/td>\r\nYoung\u2019s modulus (stiffness of a solid as the ratio of its tensile stress and axial strain, ML-1<\/sup>T-2<\/sup>)<\/td>\r\n<\/tr>\r\n
\u03bd<\/em><\/td>\r\n=<\/td>\r\nPoisson\u2019s ratio (reflects deformation of a solid in the direction perpendicular to loading, the negative ratio of transverse strain to axial strain, dimensionless)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

Typically, in layered aquifer systems laid down in a depositional environment, the geomechanical properties along the vertical direction (v<\/em>) are different from those in a horizontal direction (h<\/em>). The geomechanical properties of a transversally isotropic porous medium are fully described by five independent parameters E<\/em>v,<\/em><\/sub> E<\/em>h,<\/em><\/sub> v<\/em>v<\/em><\/sub>,<\/em><\/sub> E<\/em>h,<\/em><\/sub> G<\/em>v,<\/em><\/sub> with G<\/em> the shear modulus. G<\/em>h<\/em><\/sub> is dependent on E<\/em>h<\/em><\/sub> and v<\/em>h<\/em><\/sub> through Equation\u00a0Box\u00a03\u20114.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle G_{h}=\\frac{E_{h}}{2(1-\\nu _{h})}[\/latex]<\/td>\r\n(Box 3-4)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

Thermodynamic consistency requires the positive definiteness of matrix C <\/em>\u22121<\/sup> relating the stress tensor to the strain tensor, which implies (Ferronato et al., 2013):<\/p>\r\n

[latex]\\displaystyle 1-\\nu {_{h}}^{2}>0[\/latex] \u00a0\u00a0 and \u00a0\u00a0 [latex]\\displaystyle 1-\\nu {_{h}}-2\\nu {_{h}}^{2}\\frac{E_{v}}{E_{h}}>0[\/latex]<\/p>\r\n

then, as shown in Equations Box\u00a03\u20115, setting,<\/p>\r\n\r\n\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\vartheta =\\frac{E_{h}}{E_{v}}[\/latex]<\/td>\r\n<\/td>\r\n<\/tr>\r\n
<\/td>\r\n[latex]\\displaystyle \\eta =\\frac{E_{h}}{2(1+\\nu _{h})G_{v}}[\/latex]<\/td>\r\n(Box 3-5)<\/td>\r\n<\/tr>\r\n
<\/td>\r\n[latex]\\displaystyle \\alpha =\\frac{1}{E_{v}}\\left ( 1-\\frac{2\\nu {_{v}}^{2}}{1-\\nu _{h}}\\frac{E_{v}}{E_{h}} \\right )[\/latex]<\/td>\r\n<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

the constitutive matrix C<\/em>\u00a0<\/em>-<\/sup>1<\/sup> (equivalent to D<\/em>-1<\/sup> for a transversally isotropic medium) is Equation\u00a0Box\u00a03\u20116 (Ferronato et al., 2013).<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle C^{-1}=\\frac{1}{(1-\\nu {_{h}}^{2})\\alpha }\\begin{bmatrix}\r\nC{_{1}}^{-1} & 0 \\\\\r\n0 & C{_{2}}^{-1} \\\\\r\n\\end{bmatrix}[\/latex]<\/td>\r\n(Box 3-6)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

where:<\/p>\r\n

[latex]\\displaystyle C{_{1}}^{-1}=\\begin{bmatrix}\r\n\\vartheta -\\nu {_{v}}^{2} & \\nu {_{v}}^{2}+\\vartheta \\nu _{h} & \\nu _{v}(1+\\nu _{h}) \\\\\r\n\\nu {_{v}}^{2}+\\vartheta \\nu _{h} & \\vartheta -\\nu {_{v}}^{2} & \\nu _{v}(1+\\nu _{h}) \\\\\r\n\\nu _{v}(1+\\nu _{h}) & \\nu _{v}(1+\\nu _{h}) & 1-\\nu {_{v}}^{2} \\\\\r\n\\end{bmatrix}[\/latex]<\/p>\r\n

[latex]\\displaystyle C{_{2}}^{-1}=\\frac{\\vartheta (1-\\nu _{h})-2\\nu {_{v}}^{2}}{2}\\begin{bmatrix}\r\n1 & 0 & 0 \\\\\r\n0 & 1\/\\eta & 0 \\\\\r\n0 & 0 & 1\/\\eta \\\\\r\n\\end{bmatrix}[\/latex]<\/p>\r\n

The coefficient \u03b1 <\/span><\/em>provided in Equation\u00a0Box\u00a03\u20115 is the vertical oedometric compressibility of the medium prevented from expanding laterally (Gambolati et al., 1984). Setting v<\/em>h<\/em><\/sub>\u00a0<\/em><\/sub>=\u00a0v<\/em>v<\/em><\/sub>, E<\/em>h<\/em><\/sub>\u00a0=\u00a0E<\/em>v<\/em><\/sub>, and G<\/em>h<\/em><\/sub>\u00a0=\u00a0G<\/em>v<\/em><\/sub>, Equation\u00a0Box\u00a03\u20116 becomes Equation\u00a0Box\u00a03\u20112 and \u03b1<\/em> becomes Equation\u00a0Box\u00a03\u20117.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\alpha =\\frac{(1+\\nu )(1-2\\nu )}{(1-\\nu )E}[\/latex]<\/td>\r\n(Box 3-7)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

Equation\u00a0Box\u00a03\u20117 is the vertical compressibility of an isotropic soil. If we replace the relations between the effective stress and the strain above into the Cauchy equations, we obtain the equilibrium equations for a porous medium subject to internal pore pressure variations, p<\/em>, written in terms of displacements (isotropic medium) as in Equation\u00a0Box\u00a03\u20118.<\/p>\r\n\r\n\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle G\\triangledown ^{2}u+(\\lambda +G)\\frac{\\partial \\varepsilon }{\\partial x}=\\frac{\\partial p}{\\partial x}[\/latex]<\/td>\r\n<\/td>\r\n<\/tr>\r\n
<\/td>\r\n[latex]\\displaystyle G\\triangledown ^{2}u+(\\lambda +G)\\frac{\\partial \\varepsilon }{\\partial y}=\\frac{\\partial p}{\\partial y}[\/latex]<\/td>\r\n(Box 3-8)<\/td>\r\n<\/tr>\r\n
<\/td>\r\n[latex]\\displaystyle G\\triangledown ^{2}u+(\\lambda +G)\\frac{\\partial \\varepsilon }{\\partial z}=\\frac{\\partial p}{\\partial z}[\/latex]<\/td>\r\n<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nwhere:\r\n\r\n\r\n\r\n\r\n\r\n\r\n
u<\/em>, v<\/em>, w<\/em><\/td>\r\n=<\/td>\r\ncomponents of the incremental position vector along the coordinate axes x<\/em>, y<\/em>, and z<\/em>, respectively (L)<\/td>\r\n<\/tr>\r\n
\u22072<\/sup><\/td>\r\n=<\/td>\r\nLaplace operator<\/td>\r\n<\/tr>\r\n
\u03bb<\/em><\/td>\r\n=<\/td>\r\nLam\u00e8 constant equal to \u03bdE<\/em>\/[(1\u00a0-\u00a02\u03bd<\/em>)\u00a0(1\u00a0+\u00a0\u03bd<\/em>)]<\/td>\r\n<\/tr>\r\n
\u03b5<\/em><\/td>\r\n=<\/td>\r\n\u03b5<\/em>xx<\/em><\/sub> + \u03b5<\/em>yy<\/em><\/sub> + \u03b5<\/em>zz<\/em><\/sub>, volume strain or dilatation<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

Similar equations hold for a transversally isotropic medium, not given here, however, because of their greater complexity. There are three equations with four unknowns: u<\/em>, v<\/em>, w<\/em>, and p<\/em>. The additional equation needed to close the system is provided by the groundwater flow equation that controls subsurface flow within the aquifer.<\/p>\r\n

The flow equation is based on the principle of mass conservation for both solid grains and water. Thus, Darcy\u2019s law must be cast in terms of the relative velocity of fluid to grains. Cooper (1966) and Gambolati (1973a) derived the flow equations by assuming a grain velocity different from zero, and worked with material derivatives (total derivatives and substantial derivatives) in the appropriate places in the development. Gambolati (1973b) showed that the grain velocity can be discarded, that is, assumed to be zero, as long as the final soil settlement does not exceed 5 percent of the original aquifer thickness, a condition reached in nearly all applications. DeWiest (1966) took into consideration the dependence of the hydraulic conductivity on the water\u2019s specific weight, \u03b3<\/em>, via the intrinsic permeability and the dependence of \u03b3 <\/span><\/em>on the incremental pressure variation. Gambolati (1973b) again showed that the influence of the dependence of \u03b3<\/em> on the hydraulic conductivity is slight, and can safely be neglected. Later, within this framework, the groundwater flow equation as originally developed by Biot (1941, 1955) was elegantly and clearly derived by Verruijt (1969), and thus the fourth equation to be added to the above Equation\u00a0Box\u00a03\u20118 is Equation\u00a0Box\u00a03\u20119.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\frac{1}{\\gamma }\\triangledown \\cdot (K_{ij}\\triangledown p)=\\phi \\beta \\frac{\\partial p}{\\partial t}+\\frac{\\partial \\varepsilon }{\\partial t}[\/latex]<\/td>\r\n(Box 3-9)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nwhere:\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
\u2207<\/td>\r\n=<\/td>\r\n\u2202<\/em>\/\u2202<\/em>x<\/em> + \u2202<\/em>\/\u2202<\/em>y<\/em> + \u2202<\/em>\/\u2202<\/em>z<\/em><\/td>\r\n<\/tr>\r\n
K<\/em>ij<\/em><\/sub><\/td>\r\n=<\/td>\r\nk<\/em>ij<\/em><\/sub>\u03b3<\/em>\/\u03bc<\/em>, hydraulic conductivity tensor with principal components K<\/em>xx<\/em><\/sub>, K<\/em>yy<\/em><\/sub>, and K<\/em>zz<\/em><\/sub><\/td>\r\n<\/tr>\r\n
k<\/em>ij<\/em><\/sub><\/td>\r\n=<\/td>\r\nintrinsic permeability tensor<\/td>\r\n<\/tr>\r\n
\u03bc<\/em><\/td>\r\n=<\/td>\r\nviscosity of water<\/td>\r\n<\/tr>\r\n
\u03d5<\/em><\/td>\r\n=<\/td>\r\nmedium porosity<\/td>\r\n<\/tr>\r\n
\u03b2<\/em><\/td>\r\n=<\/td>\r\ncompressibility of water<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

Equation\u00a0Box\u00a03\u20118 together with Equation\u00a0Box\u00a03\u20119 form the mathematical basis of the so\u2011called \"coupled\" (or Biot) formulation of flow and stress in an isotropic porous medium experiencing groundwater flow. It is the most sophisticated theoretical approach to the simulation of land subsidence in the area of linear elasticity. Gambolati (1974) showed that at any point P<\/em> of the porous medium, the deformation may be expressed as the sum of two contributing factors: (1) the pointwise deformation caused by the incremental pore pressure acting at P<\/em> and (2) the deformation caused by the pressure p<\/em> acting outside P<\/em>, namely in the remainder of the medium. Gambolati (1974) called the second factor the \"three\u2011dimensional effect\": it vanishes, of course, in one\u2011dimensional media. The first factor is expressed as:<\/p>\r\n

[latex]\\displaystyle \\varepsilon =\\frac{1}{E_{v}}\\left ( 1-\\frac{2\\nu {_{v}}^{2}}{1-\\nu _{h}}\\frac{E_{v}}{E_{h}} \\right )p=\\alpha p[\/latex]<\/p>\r\n

in a geomechanical transversally isotropic medium, and<\/p>\r\n

[latex]\\displaystyle \\varepsilon =\\frac{(1+\\nu )(1-2\\nu )}{(1-\\nu )E}p=\\alpha p[\/latex]<\/p>\r\n

in a geomechanical isotropic medium, with \u03b1<\/em> the vertical compressibility previously defined. Replace the above expression for \u03b5<\/em> in the flow Equation\u00a0Box\u00a03\u20119) and you obtain the so\u2011called \"uncoupled\" formulation of flow and stress. In the uncoupled formulation the flow equation is solved for p<\/em> independently of the stress equation, with the gradient of the pore pressure variations later integrated into the equilibrium equations\u00a0(Equation\u00a0Box\u00a03\u20118) as a known external source of strength. The uncoupled flow equation is thus Equation\u00a0Box\u00a03\u201110.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\triangledown \\cdot (K_{ij}\\triangledown p)=\\gamma (\\alpha +\\phi \\beta )\\frac{\\partial p}{\\partial t}[\/latex]<\/td>\r\n(Box 3-10)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

Assuming the medium to be transversally isotropic as far as the hydraulic conductivity is concerned as well, having axes coincident with the principal directions of anisotropy, Equation\u00a0Box\u00a03\u201110 becomes Equation\u00a0Box\u00a03\u201111.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\frac{\\partial }{\\partial x}\\left ( K_{xx}\\frac{\\partial p}{\\partial x} \\right ) +\\frac{\\partial }{\\partial y}\\left ( K_{yy}\\frac{\\partial p}{\\partial y} \\right ) +\\frac{\\partial }{\\partial z}\\left ( K_{zz}\\frac{\\partial p}{\\partial z} \\right ) =S_{s}\\frac{\\partial p}{\\partial t}[\/latex]<\/td>\r\n(Box 3-11)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

The coefficient S<\/em>s<\/em><\/sub>\u00a0<\/em><\/sub>=\u00a0\u03b3<\/em>(\u03b1<\/em>\u00a0+\u00a0\u03d5\u03b2<\/em>) is the specific elastic storage coefficient referred to previously. The uncoupled equation has been the basis of classical groundwater hydrology from the very beginning of quantitative hydrogeology\u2019s development (for example, Theis, 1935; Jacob, 1940; Todd, 1960; and Bear, 1972), and is still universally used today. The superiority of the coupled approach in predicting land subsidence due to groundwater pumping has been disputed by Gambolati et al. (2000), who showed that the uncoupled pressure solution can be safely used in predicting land subsidence in compacting sedimentary basins, the coupled and uncoupled solutions being virtually indistinguishable at any time of practical interest.<\/p>\r\n

It may also be of interest to mention some basic definitions of oedometer vertical soil compressibility, which is the main rock parameter controlling land subsidence. The definition of \u03b1<\/em> given above is the one derived from the classical theory of elasticity assuming reversible elastic properties of the porous medium. The problem of defining various rock compressibilities is thoroughly discussed by Zimmerman (1991). In the present analysis, we restrict our discussion to the comparison between \u03b1<\/em> as defined above, and the compressibility c<\/em>b<\/em><\/sub> as is typically defined in geotechnics by Equation\u00a0Box\u00a03\u20114). Assume a 1\u2011D soil sample with initial length \u0394z<\/em> experiencing a vertical (oedometer) deformation \u03b4<\/em>(\u0394z<\/em>). In the classical elastic theory, the vertical compressibility \u03b1<\/em> is defined as Equation\u00a0Box\u00a03\u201112:<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\alpha =\\frac{\\delta (\\Delta z)}{\\Delta z}\\frac{1}{p}=\\frac{\\varepsilon}{p}[\/latex]<\/td>\r\n(Box 3-12)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

where<\/p>\r\n\r\n

p<\/span><\/div>\r\n, equal and opposite to the incremental effective stress, is negative in the sample compaction \u03b4<\/em>(\u0394z<\/em>). Using the void ratio, we can write Equation\u00a0Box\u00a03\u201113:\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\delta (\\Delta z)=[\\Delta z+\\delta (\\Delta z)]\\frac{e}{1+e}-\\Delta z\\frac{e_{0}}{1+e_{0}}[\/latex]<\/td>\r\n(Box 3-13)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nwhere:\r\n\r\n\r\n\r\n
e<\/em>0<\/sub><\/td>\r\n=<\/td>\r\ninitial void ratio prior to compaction (Figure 2<\/a>)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

Equation\u00a0Box\u00a03\u201113 assumes that the individual soil grains are incompressible, so that the sample volume \u03b4<\/em>(\u0394z<\/em>) is equal to the variation of the porous volume (Figure\u00a02). By dividing both sides of Equation\u00a0Box\u00a03\u201113 by \u0394z<\/em> and rearranging, we obtain Equation\u00a0Box\u00a03\u201114.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\varepsilon =\\frac{\\delta (\\Delta z)}{\\Delta z}=\\frac{e-e_{0}}{1+e_{0}}[\/latex]<\/td>\r\n(Box 3-14)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

also,<\/p>\r\n

[latex]\\displaystyle \\alpha =\\frac{\\varepsilon }{p}=\\frac{e-e_{0}}{p(1+e_{0})}[\/latex]<\/p>\r\n

and if \u03b1<\/em> does not depend on p<\/em> then Equation\u00a0Box\u00a03\u201115 can be written.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\frac{de}{dp}=\\alpha (1+e_{0})[\/latex]<\/td>\r\n(Box 3-15)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

That is, the void ratio is proportional to the incremental pressure p<\/em> (for any given initial e<\/em>0<\/sub>). Substitution of Equation\u00a0Box\u00a03\u201115 into Equation\u00a05<\/a> of the main portion of this book, with dp<\/em>\u00a0= \u2212d<\/em>\u03c3<\/em>z<\/em><\/sub>, leads to:<\/p>\r\n

[latex]\\displaystyle c_{b}=\\alpha \\frac{1+e_{0}}{1+e} =\\alpha \\frac{1+e_{0}}{1+e_{0}+\\alpha p(1+e_{0})} =\\frac{\\alpha }{1+\\alpha p}[\/latex]<\/p>\r\n

Only when the incremental pressure p<\/em> approaches 0, do \u03b1<\/em> and c<\/em>b<\/em><\/sub> coincide. In general, the two compressibilities \u03b1<\/em> and c<\/em>b <\/em><\/sub><\/span>are not equal and cannot be considered simultaneously constant. The expression of c<\/em>b<\/em><\/sub> versus \u03b5<\/em> is (using Equation\u00a0Box\u00a03\u201114) is as shown in Equation\u00a0Box\u00a03\u201116.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle c_{b}=\\frac{1+e_{0}}{1+e}\\frac{d\\varepsilon }{dp} =\\frac{1}{1+\\varepsilon }\\frac{d\\varepsilon }{dp}[\/latex]<\/td>\r\n(Box 3-16)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

If \u03b1<\/em> is constant and d<\/em>\u03b5<\/em>\/dp<\/em>\u00a0=\u00a0\u03b1<\/em> we have Equation\u00a0Box\u00a03\u201117.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle c_{b}=\\frac{\\alpha }{1+\\varepsilon }=\\frac{\\alpha }{1+\\alpha p}[\/latex]<\/td>\r\n(Box 3-17)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

Gambolati (1973b) has shown that the assumption of constant \u03b1<\/em> can be easily removed to give the general correct relationship between \u03b1<\/em> and c<\/em>b<\/em><\/sub> as in Equation\u00a0Box\u00a03\u201118.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle c_{b}=\\frac{p\\frac{d\\alpha }{dp}+\\alpha }{1+\\alpha p }[\/latex]<\/td>\r\n(Box 3-18)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

If c<\/em>b<\/em><\/sub> is constant, Equation\u00a0Box\u00a03\u201118 can be integrated to provide \u03b1<\/em> as expressed in Equation\u00a0Box\u00a03\u201119.<\/p>\r\n\r\n\r\n\r\n\r\n
<\/td>\r\n[latex]\\displaystyle \\alpha =\\frac{\\textup{exp}(pc_{b})-1}{p}[\/latex]<\/td>\r\n(Box 3-19)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n

The assumption that the individual grains are incompressible is fully warranted by the fact that the compressibility of any aquifer system is orders\u2011of\u2011magnitude greater than the compressibility of the single grain. Geerstma (1973) provides the value of \u03b1<\/em>\u00a0=\u00a01.610\u2011<\/sup>6<\/sup> bar\u2011<\/sup>1<\/sup> for grains of silicate. In contrast, the compressibility of aquifer systems is, typically, orders\u2011of\u2011magnitude larger than the compressibility of single grains as shown in Figure\u00a0Box\u00a03\u20111. Figure\u00a0Box\u00a03\u20111 provides an example of the compressibility of an aquifer system in terms of the relationship of \u03b1<\/em> versus depth and vertical effective intergranular stress \u03c3<\/em>z<\/em><\/sub>z<\/em><\/sub> in the sedimentary basin of the river Po plain, Italy (Gambolati et al., 1991, 1999; and Comerlati et al., 2004). However, as long as the ultimate relative compaction \u03b1<\/em>p<\/em> does not exceed 5 percent of the compacting unit (which is typically the case in geologic formations, particularly in shallow formations), the difference between \u03b1<\/em> and c<\/em>b<\/em><\/sub> does not exceed 2\u20113 percent (Gambolati, 1973b, Figure\u00a014) and for practical applications the two definitions are interchangeable.<\/p>\r\n\"Graph\r\n

Figure<\/strong>\u00a0<\/strong>Box\u00a0<\/strong>3<\/strong>\u20111\u00a0<\/strong>\u2011\u00a0<\/strong>Uniaxial vertical compressibility, \u03b1<\/em>, versus effective stress \u03c3<\/em>zz<\/em><\/sub> and depth z<\/em> in the Po river plain, Italy (after Comerlati et al., 2004).<\/p>\r\n

Finally, it is worth mentioning that, when comprehensive in situ and lab soil characterizations are available, more realistic constitutive formulations taking into account plastic or viscoplastic behavior may be developed and used for the simulation and prediction of land subsidence in soft under\u2011consolidated alluvial basins (for example, Ye et al., 2012).<\/p>\r\n

Return to where text linked to <\/span>Box\u00a0<\/span>3<\/span><\/a><\/p>\r\n

<\/p>\r\n\r\n<\/div>","rendered":"

\n

Theoretically, land subsidence is best analyzed according to the theory of consolidation (Biot, 1941), which holds that consolidation itself represents the response of a compressible porous medium to changes in the flow field operating within it. A complete analysis of land subsidence requires determination of the 3\u2011D deformation field accompanying the 3\u2011D flow field, and must be accomplished in a complex multi\u2011aquifer system. A few basic principles underlie the consolidation process. As outlined above, the first principle advanced by Terzaghi (1923) states that the total stress \u03c3<\/em>tot<\/em><\/sub> at any point of the porous medium is equal to the sum of the effective intergranular \u03c3<\/em>eff<\/em><\/sub> and the neutral pore pressure p<\/em>:<\/p>\n

\u03c3<\/em>tot<\/em><\/sub> = \u03c3<\/em>eff<\/em><\/sub> + p<\/em><\/p>\n

Deformation of the porous body is controlled exclusively by variation of the effective stress \u03c3<\/em>eff<\/em><\/sub>. If we consider changes relative to an initial undisturbed state of equilibrium, the Cauchy equations of equilibrium are cast in terms of incremental effective stress and pore pressure as shown in Equation\u00a0Box\u00a03\u20111.<\/p>\n\n\n\n\n\n
<\/td>\n\"\displaystyle \frac{\partial \sigma _{xx}}{\partial x}+\frac{\partial \tau _{xy}}{\partial y}+\frac{\partial \tau _{xz}}{\partial z}=\frac{\partial p}{\partial x}\"<\/td>\n<\/td>\n<\/tr>\n
<\/td>\n\"\displaystyle \frac{\partial \tau _{yx}}{\partial x}+\frac{\partial \sigma _{yy}}{\partial y}+\frac{\partial \tau _{yz}}{\partial z}=\frac{\partial p}{\partial y}\"<\/td>\n(Box 3-1)<\/td>\n<\/tr>\n
<\/td>\n\"\displaystyle \frac{\partial \tau _{zx}}{\partial x}+\frac{\partial \tau _{zy}}{\partial y}+\frac{\partial \sigma _{zz}}{\partial z}=\frac{\partial p}{\partial z}\"<\/td>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

where:<\/p>\n\n\n\n\n\n
\u03c3<\/em>xx<\/em><\/sub><\/td>\n=<\/td>\nincremental normal effective stress in x direction<\/td>\n<\/tr>\n
\u03c3<\/em>yy<\/em><\/sub><\/td>\n=<\/td>\nincremental normal effective stress in y direction<\/td>\n<\/tr>\n
\u03c3<\/em>zz<\/em><\/sub><\/td>\n=<\/td>\nincremental normal effective stress in z direction<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

and, the incremental shear stresses:<\/p>\n\n\n\n\n\n
\u03c4<\/em>xy<\/em><\/sub><\/td>\n=<\/td>\n\u03c4<\/em>yx<\/em><\/sub><\/td>\n<\/tr>\n
\u03c4<\/em>xz<\/em><\/sub><\/td>\n=<\/td>\n\u03c4<\/em>zx<\/em><\/sub><\/td>\n<\/tr>\n
\u03c4<\/em>yz<\/em><\/sub><\/td>\n=<\/td>\n\u03c4<\/em>zy<\/em><\/sub><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

The relationships between the incremental effective stress tensor \u03c3<\/em> and the incremental strain tensor \u03b5<\/em> for a geomechanical isotropic medium are shown in Equation\u00a0Box\u00a03\u20112.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle \begin{bmatrix}<\/td>\n(Box 3-2)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

with matrix D<\/em>-1<\/sup> as shown in Equation\u00a0Box\u00a03\u20113.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle D^{-1}=\frac{E}{(1+\nu )(1-2\nu )} \begin{bmatrix}<\/td>\n(Box 3-3)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

where:<\/p>\n\n\n\n\n
E<\/em><\/td>\n=<\/td>\nYoung\u2019s modulus (stiffness of a solid as the ratio of its tensile stress and axial strain, ML-1<\/sup>T-2<\/sup>)<\/td>\n<\/tr>\n
\u03bd<\/em><\/td>\n=<\/td>\nPoisson\u2019s ratio (reflects deformation of a solid in the direction perpendicular to loading, the negative ratio of transverse strain to axial strain, dimensionless)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Typically, in layered aquifer systems laid down in a depositional environment, the geomechanical properties along the vertical direction (v<\/em>) are different from those in a horizontal direction (h<\/em>). The geomechanical properties of a transversally isotropic porous medium are fully described by five independent parameters E<\/em>v,<\/em><\/sub> E<\/em>h,<\/em><\/sub> v<\/em>v<\/em><\/sub>,<\/em><\/sub> E<\/em>h,<\/em><\/sub> G<\/em>v,<\/em><\/sub> with G<\/em> the shear modulus. G<\/em>h<\/em><\/sub> is dependent on E<\/em>h<\/em><\/sub> and v<\/em>h<\/em><\/sub> through Equation\u00a0Box\u00a03\u20114.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle G_{h}=\frac{E_{h}}{2(1-\nu _{h})}\"<\/td>\n(Box 3-4)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Thermodynamic consistency requires the positive definiteness of matrix C <\/em>\u22121<\/sup> relating the stress tensor to the strain tensor, which implies (Ferronato et al., 2013):<\/p>\n

\"\displaystyle 1-\nu {_{h}}^{2}>0\" \u00a0\u00a0 and \u00a0\u00a0 \"\displaystyle 1-\nu {_{h}}-2\nu {_{h}}^{2}\frac{E_{v}}{E_{h}}>0\"<\/p>\n

then, as shown in Equations Box\u00a03\u20115, setting,<\/p>\n\n\n\n\n\n
<\/td>\n\"\displaystyle \vartheta =\frac{E_{h}}{E_{v}}\"<\/td>\n<\/td>\n<\/tr>\n
<\/td>\n\"\displaystyle \eta =\frac{E_{h}}{2(1+\nu _{h})G_{v}}\"<\/td>\n(Box 3-5)<\/td>\n<\/tr>\n
<\/td>\n\"\displaystyle \alpha =\frac{1}{E_{v}}\left ( 1-\frac{2\nu {_{v}}^{2}}{1-\nu _{h}}\frac{E_{v}}{E_{h}} \right )\"<\/td>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

the constitutive matrix C<\/em>\u00a0<\/em>–<\/sup>1<\/sup> (equivalent to D<\/em>-1<\/sup> for a transversally isotropic medium) is Equation\u00a0Box\u00a03\u20116 (Ferronato et al., 2013).<\/p>\n\n\n\n
<\/td>\n\"\displaystyle C^{-1}=\frac{1}{(1-\nu {_{h}}^{2})\alpha }\begin{bmatrix}<\/td>\n(Box 3-6)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

where:<\/p>\n

\"\displaystyle C{_{1}}^{-1}=\begin{bmatrix}<\/p>\n

\"\displaystyle C{_{2}}^{-1}=\frac{\vartheta (1-\nu _{h})-2\nu {_{v}}^{2}}{2}\begin{bmatrix}<\/p>\n

The coefficient \u03b1 <\/span><\/em>provided in Equation\u00a0Box\u00a03\u20115 is the vertical oedometric compressibility of the medium prevented from expanding laterally (Gambolati et al., 1984). Setting v<\/em>h<\/em><\/sub>\u00a0<\/em><\/sub>=\u00a0v<\/em>v<\/em><\/sub>, E<\/em>h<\/em><\/sub>\u00a0=\u00a0E<\/em>v<\/em><\/sub>, and G<\/em>h<\/em><\/sub>\u00a0=\u00a0G<\/em>v<\/em><\/sub>, Equation\u00a0Box\u00a03\u20116 becomes Equation\u00a0Box\u00a03\u20112 and \u03b1<\/em> becomes Equation\u00a0Box\u00a03\u20117.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle \alpha =\frac{(1+\nu )(1-2\nu )}{(1-\nu )E}\"<\/td>\n(Box 3-7)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Equation\u00a0Box\u00a03\u20117 is the vertical compressibility of an isotropic soil. If we replace the relations between the effective stress and the strain above into the Cauchy equations, we obtain the equilibrium equations for a porous medium subject to internal pore pressure variations, p<\/em>, written in terms of displacements (isotropic medium) as in Equation\u00a0Box\u00a03\u20118.<\/p>\n\n\n\n\n\n
<\/td>\n\"\displaystyle G\triangledown ^{2}u+(\lambda +G)\frac{\partial \varepsilon }{\partial x}=\frac{\partial p}{\partial x}\"<\/td>\n<\/td>\n<\/tr>\n
<\/td>\n\"\displaystyle G\triangledown ^{2}u+(\lambda +G)\frac{\partial \varepsilon }{\partial y}=\frac{\partial p}{\partial y}\"<\/td>\n(Box 3-8)<\/td>\n<\/tr>\n
<\/td>\n\"\displaystyle G\triangledown ^{2}u+(\lambda +G)\frac{\partial \varepsilon }{\partial z}=\frac{\partial p}{\partial z}\"<\/td>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

where:<\/p>\n\n\n\n\n\n\n
u<\/em>, v<\/em>, w<\/em><\/td>\n=<\/td>\ncomponents of the incremental position vector along the coordinate axes x<\/em>, y<\/em>, and z<\/em>, respectively (L)<\/td>\n<\/tr>\n
\u22072<\/sup><\/td>\n=<\/td>\nLaplace operator<\/td>\n<\/tr>\n
\u03bb<\/em><\/td>\n=<\/td>\nLam\u00e8 constant equal to \u03bdE<\/em>\/[(1\u00a0–\u00a02\u03bd<\/em>)\u00a0(1\u00a0+\u00a0\u03bd<\/em>)]<\/td>\n<\/tr>\n
\u03b5<\/em><\/td>\n=<\/td>\n\u03b5<\/em>xx<\/em><\/sub> + \u03b5<\/em>yy<\/em><\/sub> + \u03b5<\/em>zz<\/em><\/sub>, volume strain or dilatation<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Similar equations hold for a transversally isotropic medium, not given here, however, because of their greater complexity. There are three equations with four unknowns: u<\/em>, v<\/em>, w<\/em>, and p<\/em>. The additional equation needed to close the system is provided by the groundwater flow equation that controls subsurface flow within the aquifer.<\/p>\n

The flow equation is based on the principle of mass conservation for both solid grains and water. Thus, Darcy\u2019s law must be cast in terms of the relative velocity of fluid to grains. Cooper (1966) and Gambolati (1973a) derived the flow equations by assuming a grain velocity different from zero, and worked with material derivatives (total derivatives and substantial derivatives) in the appropriate places in the development. Gambolati (1973b) showed that the grain velocity can be discarded, that is, assumed to be zero, as long as the final soil settlement does not exceed 5 percent of the original aquifer thickness, a condition reached in nearly all applications. DeWiest (1966) took into consideration the dependence of the hydraulic conductivity on the water\u2019s specific weight, \u03b3<\/em>, via the intrinsic permeability and the dependence of \u03b3 <\/span><\/em>on the incremental pressure variation. Gambolati (1973b) again showed that the influence of the dependence of \u03b3<\/em> on the hydraulic conductivity is slight, and can safely be neglected. Later, within this framework, the groundwater flow equation as originally developed by Biot (1941, 1955) was elegantly and clearly derived by Verruijt (1969), and thus the fourth equation to be added to the above Equation\u00a0Box\u00a03\u20118 is Equation\u00a0Box\u00a03\u20119.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle \frac{1}{\gamma }\triangledown \cdot (K_{ij}\triangledown p)=\phi \beta \frac{\partial p}{\partial t}+\frac{\partial \varepsilon }{\partial t}\"<\/td>\n(Box 3-9)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

where:<\/p>\n\n\n\n\n\n\n\n\n
\u2207<\/td>\n=<\/td>\n\u2202<\/em>\/\u2202<\/em>x<\/em> + \u2202<\/em>\/\u2202<\/em>y<\/em> + \u2202<\/em>\/\u2202<\/em>z<\/em><\/td>\n<\/tr>\n
K<\/em>ij<\/em><\/sub><\/td>\n=<\/td>\nk<\/em>ij<\/em><\/sub>\u03b3<\/em>\/\u03bc<\/em>, hydraulic conductivity tensor with principal components K<\/em>xx<\/em><\/sub>, K<\/em>yy<\/em><\/sub>, and K<\/em>zz<\/em><\/sub><\/td>\n<\/tr>\n
k<\/em>ij<\/em><\/sub><\/td>\n=<\/td>\nintrinsic permeability tensor<\/td>\n<\/tr>\n
\u03bc<\/em><\/td>\n=<\/td>\nviscosity of water<\/td>\n<\/tr>\n
\u03d5<\/em><\/td>\n=<\/td>\nmedium porosity<\/td>\n<\/tr>\n
\u03b2<\/em><\/td>\n=<\/td>\ncompressibility of water<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Equation\u00a0Box\u00a03\u20118 together with Equation\u00a0Box\u00a03\u20119 form the mathematical basis of the so\u2011called “coupled” (or Biot) formulation of flow and stress in an isotropic porous medium experiencing groundwater flow. It is the most sophisticated theoretical approach to the simulation of land subsidence in the area of linear elasticity. Gambolati (1974) showed that at any point P<\/em> of the porous medium, the deformation may be expressed as the sum of two contributing factors: (1) the pointwise deformation caused by the incremental pore pressure acting at P<\/em> and (2) the deformation caused by the pressure p<\/em> acting outside P<\/em>, namely in the remainder of the medium. Gambolati (1974) called the second factor the “three\u2011dimensional effect”: it vanishes, of course, in one\u2011dimensional media. The first factor is expressed as:<\/p>\n

\"\displaystyle \varepsilon =\frac{1}{E_{v}}\left ( 1-\frac{2\nu {_{v}}^{2}}{1-\nu _{h}}\frac{E_{v}}{E_{h}} \right )p=\alpha p\"<\/p>\n

in a geomechanical transversally isotropic medium, and<\/p>\n

\"\displaystyle \varepsilon =\frac{(1+\nu )(1-2\nu )}{(1-\nu )E}p=\alpha p\"<\/p>\n

in a geomechanical isotropic medium, with \u03b1<\/em> the vertical compressibility previously defined. Replace the above expression for \u03b5<\/em> in the flow Equation\u00a0Box\u00a03\u20119) and you obtain the so\u2011called “uncoupled” formulation of flow and stress. In the uncoupled formulation the flow equation is solved for p<\/em> independently of the stress equation, with the gradient of the pore pressure variations later integrated into the equilibrium equations\u00a0(Equation\u00a0Box\u00a03\u20118) as a known external source of strength. The uncoupled flow equation is thus Equation\u00a0Box\u00a03\u201110.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle \triangledown \cdot (K_{ij}\triangledown p)=\gamma (\alpha +\phi \beta )\frac{\partial p}{\partial t}\"<\/td>\n(Box 3-10)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Assuming the medium to be transversally isotropic as far as the hydraulic conductivity is concerned as well, having axes coincident with the principal directions of anisotropy, Equation\u00a0Box\u00a03\u201110 becomes Equation\u00a0Box\u00a03\u201111.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle \frac{\partial }{\partial x}\left ( K_{xx}\frac{\partial p}{\partial x} \right ) +\frac{\partial }{\partial y}\left ( K_{yy}\frac{\partial p}{\partial y} \right ) +\frac{\partial }{\partial z}\left ( K_{zz}\frac{\partial p}{\partial z} \right ) =S_{s}\frac{\partial p}{\partial t}\"<\/td>\n(Box 3-11)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

The coefficient S<\/em>s<\/em><\/sub>\u00a0<\/em><\/sub>=\u00a0\u03b3<\/em>(\u03b1<\/em>\u00a0+\u00a0\u03d5\u03b2<\/em>) is the specific elastic storage coefficient referred to previously. The uncoupled equation has been the basis of classical groundwater hydrology from the very beginning of quantitative hydrogeology\u2019s development (for example, Theis, 1935; Jacob, 1940; Todd, 1960; and Bear, 1972), and is still universally used today. The superiority of the coupled approach in predicting land subsidence due to groundwater pumping has been disputed by Gambolati et al. (2000), who showed that the uncoupled pressure solution can be safely used in predicting land subsidence in compacting sedimentary basins, the coupled and uncoupled solutions being virtually indistinguishable at any time of practical interest.<\/p>\n

It may also be of interest to mention some basic definitions of oedometer vertical soil compressibility, which is the main rock parameter controlling land subsidence. The definition of \u03b1<\/em> given above is the one derived from the classical theory of elasticity assuming reversible elastic properties of the porous medium. The problem of defining various rock compressibilities is thoroughly discussed by Zimmerman (1991). In the present analysis, we restrict our discussion to the comparison between \u03b1<\/em> as defined above, and the compressibility c<\/em>b<\/em><\/sub> as is typically defined in geotechnics by Equation\u00a0Box\u00a03\u20114). Assume a 1\u2011D soil sample with initial length \u0394z<\/em> experiencing a vertical (oedometer) deformation \u03b4<\/em>(\u0394z<\/em>). In the classical elastic theory, the vertical compressibility \u03b1<\/em> is defined as Equation\u00a0Box\u00a03\u201112:<\/p>\n\n\n\n
<\/td>\n\"\displaystyle \alpha =\frac{\delta (\Delta z)}{\Delta z}\frac{1}{p}=\frac{\varepsilon}{p}\"<\/td>\n(Box 3-12)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

where<\/p>\n

p<\/span><\/div>\n

, equal and opposite to the incremental effective stress, is negative in the sample compaction \u03b4<\/em>(\u0394z<\/em>). Using the void ratio, we can write Equation\u00a0Box\u00a03\u201113:<\/p>\n\n\n\n
<\/td>\n\"\displaystyle \delta (\Delta z)=[\Delta z+\delta (\Delta z)]\frac{e}{1+e}-\Delta z\frac{e_{0}}{1+e_{0}}\"<\/td>\n(Box 3-13)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

where:<\/p>\n\n\n\n
e<\/em>0<\/sub><\/td>\n=<\/td>\ninitial void ratio prior to compaction (Figure 2<\/a>)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Equation\u00a0Box\u00a03\u201113 assumes that the individual soil grains are incompressible, so that the sample volume \u03b4<\/em>(\u0394z<\/em>) is equal to the variation of the porous volume (Figure\u00a02). By dividing both sides of Equation\u00a0Box\u00a03\u201113 by \u0394z<\/em> and rearranging, we obtain Equation\u00a0Box\u00a03\u201114.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle \varepsilon =\frac{\delta (\Delta z)}{\Delta z}=\frac{e-e_{0}}{1+e_{0}}\"<\/td>\n(Box 3-14)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

also,<\/p>\n

\"\displaystyle \alpha =\frac{\varepsilon }{p}=\frac{e-e_{0}}{p(1+e_{0})}\"<\/p>\n

and if \u03b1<\/em> does not depend on p<\/em> then Equation\u00a0Box\u00a03\u201115 can be written.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle \frac{de}{dp}=\alpha (1+e_{0})\"<\/td>\n(Box 3-15)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

That is, the void ratio is proportional to the incremental pressure p<\/em> (for any given initial e<\/em>0<\/sub>). Substitution of Equation\u00a0Box\u00a03\u201115 into Equation\u00a05<\/a> of the main portion of this book, with dp<\/em>\u00a0= \u2212d<\/em>\u03c3<\/em>z<\/em><\/sub>, leads to:<\/p>\n

\"\displaystyle c_{b}=\alpha \frac{1+e_{0}}{1+e} =\alpha \frac{1+e_{0}}{1+e_{0}+\alpha p(1+e_{0})} =\frac{\alpha }{1+\alpha p}\"<\/p>\n

Only when the incremental pressure p<\/em> approaches 0, do \u03b1<\/em> and c<\/em>b<\/em><\/sub> coincide. In general, the two compressibilities \u03b1<\/em> and c<\/em>b <\/em><\/sub><\/span>are not equal and cannot be considered simultaneously constant. The expression of c<\/em>b<\/em><\/sub> versus \u03b5<\/em> is (using Equation\u00a0Box\u00a03\u201114) is as shown in Equation\u00a0Box\u00a03\u201116.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle c_{b}=\frac{1+e_{0}}{1+e}\frac{d\varepsilon }{dp} =\frac{1}{1+\varepsilon }\frac{d\varepsilon }{dp}\"<\/td>\n(Box 3-16)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

If \u03b1<\/em> is constant and d<\/em>\u03b5<\/em>\/dp<\/em>\u00a0=\u00a0\u03b1<\/em> we have Equation\u00a0Box\u00a03\u201117.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle c_{b}=\frac{\alpha }{1+\varepsilon }=\frac{\alpha }{1+\alpha p}\"<\/td>\n(Box 3-17)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Gambolati (1973b) has shown that the assumption of constant \u03b1<\/em> can be easily removed to give the general correct relationship between \u03b1<\/em> and c<\/em>b<\/em><\/sub> as in Equation\u00a0Box\u00a03\u201118.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle c_{b}=\frac{p\frac{d\alpha }{dp}+\alpha }{1+\alpha p }\"<\/td>\n(Box 3-18)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

If c<\/em>b<\/em><\/sub> is constant, Equation\u00a0Box\u00a03\u201118 can be integrated to provide \u03b1<\/em> as expressed in Equation\u00a0Box\u00a03\u201119.<\/p>\n\n\n\n
<\/td>\n\"\displaystyle \alpha =\frac{\textup{exp}(pc_{b})-1}{p}\"<\/td>\n(Box 3-19)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

The assumption that the individual grains are incompressible is fully warranted by the fact that the compressibility of any aquifer system is orders\u2011of\u2011magnitude greater than the compressibility of the single grain. Geerstma (1973) provides the value of \u03b1<\/em>\u00a0=\u00a01.610\u2011<\/sup>6<\/sup> bar\u2011<\/sup>1<\/sup> for grains of silicate. In contrast, the compressibility of aquifer systems is, typically, orders\u2011of\u2011magnitude larger than the compressibility of single grains as shown in Figure\u00a0Box\u00a03\u20111. Figure\u00a0Box\u00a03\u20111 provides an example of the compressibility of an aquifer system in terms of the relationship of \u03b1<\/em> versus depth and vertical effective intergranular stress \u03c3<\/em>z<\/em><\/sub>z<\/em><\/sub> in the sedimentary basin of the river Po plain, Italy (Gambolati et al., 1991, 1999; and Comerlati et al., 2004). However, as long as the ultimate relative compaction \u03b1<\/em>p<\/em> does not exceed 5 percent of the compacting unit (which is typically the case in geologic formations, particularly in shallow formations), the difference between \u03b1<\/em> and c<\/em>b<\/em><\/sub> does not exceed 2\u20113 percent (Gambolati, 1973b, Figure\u00a014) and for practical applications the two definitions are interchangeable.<\/p>\n

\"Graph<\/p>\n

Figure<\/strong>\u00a0<\/strong>Box\u00a0<\/strong>3<\/strong>\u20111\u00a0<\/strong>\u2011\u00a0<\/strong>Uniaxial vertical compressibility, \u03b1<\/em>, versus effective stress \u03c3<\/em>zz<\/em><\/sub> and depth z<\/em> in the Po river plain, Italy (after Comerlati et al., 2004).<\/p>\n

Finally, it is worth mentioning that, when comprehensive in situ and lab soil characterizations are available, more realistic constitutive formulations taking into account plastic or viscoplastic behavior may be developed and used for the simulation and prediction of land subsidence in soft under\u2011consolidated alluvial basins (for example, Ye et al., 2012).<\/p>\n

Return to where text linked to <\/span>Box\u00a0<\/span>3<\/span><\/a><\/p>\n

\n<\/div>\n","protected":false},"author":1,"menu_order":31,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-276","chapter","type-chapter","status-publish","hentry"],"part":177,"_links":{"self":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapters\/276","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/wp\/v2\/users\/1"}],"version-history":[{"count":14,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapters\/276\/revisions"}],"predecessor-version":[{"id":466,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapters\/276\/revisions\/466"}],"part":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/parts\/177"}],"metadata":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapters\/276\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/wp\/v2\/media?parent=276"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/pressbooks\/v2\/chapter-type?post=276"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/wp\/v2\/contributor?post=276"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/books.gw-project.org\/land-subsidence-and-its-mitigation\/wp-json\/wp\/v2\/license?post=276"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}