An advanced discrete fracture model for variably saturated flow in fractured porous media
KAUST DepartmentAli I. Al-Naimi Petroleum Engineering Research Center (ANPERC)
Energy Resources and Petroleum Engineering
Energy Resources and Petroleum Engineering Program
Physical Science and Engineering (PSE) Division
Online Publication Date2020-05-06
Print Publication Date2020-06
Embargo End Date2021-05-11
Permanent link to this recordhttp://hdl.handle.net/10754/662841
MetadataShow full item record
AbstractAccurate modeling of variably saturated flow (VSF) in fractured porous media with the discrete fracture-matrix (DFM) model is a computationally challenging problem. The applicability of DFM model to VSF in real field studies at large space and time scales is often limited, not only because it requires detailed fracture characterization, but also as it involves excessive computational efforts. We develop an efficient numerical scheme to solve the Richards equation in discretely fractured porous media. This scheme combines the mixed hybrid finite element method for space discretization with the method of lines for time integration. The fractures are modeled as lower-dimensional interfaces (1D), within the 2D porous matrix. We develop a new mass-lumping (ML) technique for the fractures to eliminate unphysical oscillations and convergence issues in the solution, which significantly improves efficiency, enabling larger field applications. The proposed new scheme is validated against a commercial simulator for problems involving water table recharge at the laboratory scale. The computational efficiency of the developed scheme is examined on a challenging problem for water infiltration in fractured dry soil, and compared with standard numerical techniques. We show that the ML technique is crucial to improve robustness and efficiency, which outperforms the commonly used methods that we tested. The applicability of our method is then demonstrated in a study concerning the effect of climate change on groundwater resources in a karst aquifer/spring system in El Assal, Lebanon. Simulations, including recharge predictions under climate change scenarios, are carried out for about 80 years, up to 2099. This study demonstrates the applicability of our proposed scheme to deal with real field cases involving large time and space scales with high variable recharge. Our results indicate that the water-table level is sensitive to the presence of fractures, where neglecting fractures leads to an overestimation of the available groundwater amount. The proposed numerical approach is generic for DFM model and can be extended to different 2D and 3D finite-element frameworks.
CitationKoohbor, B., Fahs, M., Hoteit, H., Doummar, J., Younes, A., & Belfort, B. (2020). An advanced discrete fracture model for variably saturated flow in fractured porous media. Advances in Water Resources, 140, 103602. doi:10.1016/j.advwatres.2020.103602
SponsorsMarwan Fahs would acknowledge the support from the national school of water and environmental engineering of Strasbourg through the research project PORO6100. The characterization of the case study, including data collection was funded in the framework of a USAID PEER Science project Award number (102881; Cycle 3). The data used in this work are available on the GitHub repository: https://github.com/fahs-LHYGES.
JournalAdvances in Water Resources