Abstract
A fully coupled mathematical model of land subsidence caused by groundwater pumping was established based on the mechanics of porous seepage and the theory of fluid–solid interaction. The mathematical model employing the Galerkin finite element method was proposed to simulate the deformation dependencies of hydraulic properties due to the water pressure decrease in aquifers. This model has been verified by comparing with the known analytical solutions in the confined aquifer. Results show that the simulated drawdown and displacements match well with those of analytical solutions. To evaluate the nonlinear effects of hydraulic properties in the seepage and consolidation coupling phenomena, the numerical model is applied to an ideal three layers numerical experiment. The results show that the subsidence rate is faster than conventional groundwater theory when the nonlinear hydraulic properties (NHP) are considered in the coupled model. The reason is that the water level decline due to groundwater withdrawal induces the soil compression and the porosity and permeability decrease. Decrease of permeability in regions adjacent to the pumping well produces hydraulic gradient and seepage forces that may result in accelerated subsidence. Therefore, the effect of the NHP is an interaction process of pore water pressure and soil consolidation deformation and should not be ignored. Further studies of various hydrogeology problems are recommended to consider the heterogeneity concerning different stratigraphic condition in the field.
Similar content being viewed by others
References
Abidin HZ, Gumilar I, Andreas H, Murdohardono D, Fukuda Y (2013) On causes and impacts of land subsidence in Bandung Basin, Indonesia. Environ Earth Sci 68(6):1545–1553
Bai M, Elsworth D (1994) Modeling of subsidence and stress dependent hydraulic conductivity for intact and fractured porous media. Rock Mech Rock Eng 27(4):209–234
Bear J, Corapcioglu MY (1981) Mathematical model for regional land subsidence due to pumping, 2. Integrated aquifer subsidence equations for vertical and horizontal displacements. Water Resour Res 17:947–958
Biot MA (1941) General theory of three-dimensional consolidation. J Appl Phys 12:144–164
Biot MA (1955) Theory of elasticity and consolidation for a porous anisotropic solid. J Appl Phys 26(2):182–185
Budhu M, Adiyaman B (2010) Mechanics of land subsidence due to groundwater pumping. Int J Numer Anal Methods Geomech 34(14):1459–1478
Budhu M, Adiyaman B (2013) The influence of clay zones on land subsidence from groundwater pumping. GroundWater 51(1):51–57
Carman PC (1937) Fluid flow through granular beds. Trans Inst Chem Eng 15:150–166
Gambolati G, Teatini P, Baú D, Ferronato M (2000) Importance of poroelastic coupling in dynamically active aquifers of the Po river basin, Italy. Water Resour Res 36(9):2443–2460
Gatmiri B, Delage P (1997) A formulation of fully coupled thermal–hydraulic–mechanical behaviour of saturated porous media—numerical approach. Int J Numer Anal Methods Geomech 21(3):199–225
Hsieh PC (2006) A viscoelastic model for the dynamic response of soils to periodical surface water disturbance. Int J Numer Anal Methods Geomech 30(12):1201–1212
Huang H, Wattenbarger RC, Gai XL, Brown WP, Hehmeyer OJ, Wang JL, Long TA (2013) Using a fully coupled flow and geomechanical simulator to model injection into heavy oil reservoirs. Int J Numer Methods Fluids 71(6):671–686
Huyakorn PS, Springer EP, Guvanasen V, Wadsworth TD (1986) A three-dimensional finite-element model for simulating water flow in variably saturated porous media. Water Resour Res 22(13):1790–1808
Khan AS, Khan SD, Kakar DM (2013) Land subsidence and declining water resources in Quetta Valley, Pakistan. Environ Earth Sci 70(6):2719–2727
Kim JM (2000) A fully coupled finite element analysis of water-table fluctuation and land deformation in partially saturated soils due to surface loading. Int J Numer Anal Methods Geomech 49:1101–1119
Kim JM, Parizek RR (1999) A mathematical model for the hydraulic properties of deforming porous media. GroundWater 37(4):546–554
Kim JM, Parizek RR (2005) Numerical simulation of the Rhade effect in layered aquifer systems due to groundwater pumping shutoff. Adv Water Resour 28:627–642
Kozeny J (1927) Über kapillare Leitung des Wassers im Boden. Sitzungsber Akad Wiss Wien 136:271–306
Leake SA, Hsieh PA (1997) Simulation of deformation of sediments from decline of groundwater levels in an aquifer underlain by a bedrock step. USGS Open File Report
Lewis RW, Schrefler BA (1978) A fully coupled consolidation model of the subsidence of Venice. Water Resour Res 14(2):223–230
Lewis RW, Schrefler BA (1998) The finite element method in the static and dynamic deformation and consolidation of porous media. Wiley, New York
Li PC, Kong XY, Lu DT (2003) Mathematical modeling of flow in saturated porous media on account of fluid-structure coupling effect. J Hydrodyn (In Chinese) 18(4):419–426
Ouria A, Toufigh MM, Nakhai A (2009) An investigation on the effect of the coupled and uncoupled formulation on transient seepage by the finite element method. Am J Appl Sci 4(12):950–956
Ouria A, Toufigh MM, Nakhai A (2011) Nonlinear analysis of transient seepage by the coupled finite element method. Int J Mech 5(1):35–39
Safai NM, Pinder GF (1979) Vertical and horizontal land deformation in a desaturating porous medium. Adv Water Resour 2:19–25
Schrefler BA, Zhan XY (1993) A fully coupled model for water flow and airflow in deformable porous media. Water Resour Res 29(1):155–167
Terzaghi K (1943) Theoretical Soil Mechanics. Wiley, New York
Thomée V (2006) Galerkin finite element methods for parabolic problems, 2nd edn. Springer, Berlin
Wang HF (2000) Theory of linear poroelasticity with application to geomechanics and hydrogeology. Princeton University Press, Princeton
Xue Q, Feng XT, Liu JJ (2005) Study on multi-field coupled model and numerical simulation of landfill gas transport. J Syst Simul (In Chinese) 17:20–24
Yeh HD, Lu RH, Yeh GT (1996) Finite element modeling for land displacements due to pumping. Int J Numer Anal Methods Geomech 20(2):79–99
Zhang CH, Zhao QS, Yu YJ (2011) Model of coupled gas flow and deformation process in heterogeneous coal seams and its application. J Coal Sci Eng 17:76–80
Zimmerman RW (2000) Coupling in poroelasticity and thermoelasticity. Int J Rock Mech Min Sci 37:79–87
Zimmerman RW, Somerton WH, King MS (1986) Compressibility of porous rock. J Geophys Res 91:12765–12777
Acknowledgments
This work was supported by Beijing Science and Technology Commission project “Research on water resources sustainable utilization of new air area in Beijing” (Z131100005613001).
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Yang, Y., Song, X.F., Zheng, F.D. et al. Simulation of fully coupled finite element analysis of nonlinear hydraulic properties in land subsidence due to groundwater pumping. Environ Earth Sci 73, 4191–4199 (2015). https://doi.org/10.1007/s12665-014-3705-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s12665-014-3705-8