Teatini, and 5 co-authors, reported in 2015 that the stress variation induced by aquifer overdraft in sedimentary basins with shallow bedrock may cause rupture in the form of pre-existing fault activation or earth fissure generation. The process is causing major detrimental effects on a many areas in China and Mexico. Ruptures yield discontinuity in both displacement and stress field that classic continuous finite element (FE) models cannot address. Interface finite elements (IE), typically used in contact mechanics, may be of great help and are implemented herein to simulate the fault geomechanical behavior. Two main approaches, i.e. Penalty and Lagrangian, are developed to enforce the contact condition on the element interface.
The incorporation of IE incorporation into a three-dimensional (3-D) FE geomechanical simulator shows that the Lagrangian approach is numerically more robust and stable than the Penalty, thus providing more reliable solutions. Furthermore, the use of a Newton-Raphson scheme to deal with the non-linear elasto-plastic fault behavior allows for quadratic convergence. The FE – IE model is applied to investigate the likely ground rupture in realistic 3-D geologic settings. The case studies are representative of the City of Wuxi in the Jiangsu Province (China), and of the City of Queretaro, Mexico, where significant land subsidence has been accompanied by the generation of several earth fissures jeopardizing the stability and integrity of the overland structures and infrastructure.
For a similar characterization of over-pumping groundwater and subsidence within unconsolidated to partly consolidated sediments in and around Houston, Texas, see:
Campbell, et al., 2018: http://www.i2mconsultants.com/downloads/corp/JGG-1-015.pdf