Learning multi-phase flow and transport in fractured porous media with auto-regressive and recurrent graph neural networks

Journal Article (2026)
Author(s)

Mohammed Al Kobaisi (TU Delft - Electrical Engineering, Mathematics and Computer Science, Khalifa University of Science and Technology)

Wenjuan Zhang (Qingdao University of Technology)

Waleed Diab (Khalifa University of Science and Technology)

Hadi Hajibeygi (TU Delft - Civil Engineering & Geosciences)

Research Group
Numerical Analysis
DOI related publication
https://doi.org/10.1016/j.ees.2025.11.002 Final published version
More Info
expand_more
Publication Year
2026
Language
English
Research Group
Numerical Analysis
Journal title
Earth Energy Science
Issue number
1
Volume number
2
Article number
100035
Downloads counter
43
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Abstract

In the past three decades, a wide array of computational methodologies and simulation frameworks have emerged to address the complexities of modeling flow and transport processes in fractured porous media. The conformal mesh approaches which explicitly align the computational grid with fracture surfaces are considered by many to be the most accurate. However, such methods require excessive fine-scale meshing, rendering them impractical for large or complex fracture networks. The Embedded Discrete Fracture Model (EDFM) offers a good balance between accuracy and efficiency and has gained a lot of traction in recent years. Nonetheless, it is not free of drawbacks; EDFM can, and often will, generate fracture cells that have orders of magnitudes smaller volumes than the matrix cells, which significantly impacts the convergence performance of nonlinear solvers. In this work, we propose to learn the complex flow and transport dynamics in fractured porous media with graph neural networks (GNN). GNNs are well suited for this task due to the unstructured topology of the computation grid resulting from EDFM discretization. We propose two deep learning architectures, a GNN and a recurrent GNN. Both networks follow a two-stage training strategy: an autoregressive one step roll-out, followed by a fine-tuning step where the model is supervised using the whole ground-truth sequence. We demonstrate that the two-stage training approach is effective in mitigating error accumulation during autoregressive model rollouts in the testing phase. Our findings indicate that both GNNs generalize well to unseen fracture realizations. While the second stage of training proved to be beneficial for the GNN model, its impact on the recurrent GNN model was less pronounced. Finally, the performance of both GNNs for temporal extrapolation is tested. The recurrent GNN significantly outperformed the GNN in terms of accuracy, thereby underscoring its superior capability in predicting long sequences.