Reduced-Order Modeling for Subsurface Flow Simulation in Fractured Reservoirs

Ming Gao,Wenyue Sun,Jianchun Xu,Junchao Li
DOI: https://doi.org/10.2118/223949-pa
2024-01-01
SPE Journal
Abstract:Summary Reservoir simulation for fractured reservoirs is often challenging and time-consuming due to the strong heterogeneity and complex flow dynamics introduced by fracture-matrix interactions. In this study, we introduce a novel reduced-order modeling procedure to speed up the flow simulation of fractured reservoirs. The reduced-order model (ROM) is developed based on proper orthogonal decomposition (POD) in conjunction with the embedded discrete fracture model (EDFM) that provides full-order simulation results. With the full-order training simulation, snapshots of reservoir pressure and saturation state at different timesteps are captured and assembled into separate data matrices. Singular value decomposition (SVD) is then applied to these data matrices to obtain a reduced set of orthogonal base vectors for pressure and saturation solutions, respectively. These base vectors enable the projection of high-dimensional linear equations into much lower-dimensional spaces, which significantly accelerates the process of solving nonlinear governing equations under the EDFM approach. The developed reduced-order modeling procedure is implemented in the MATLAB reservoir simulation toolbox (MRST) and tested via multiple cases for both 2D and 3D fractured reservoirs under different boundary and well control scenarios. In certain challenging cases, the use of multiple training simulations is explored and is shown to provide improved predictions. Overall, the proposed ROM approach is able to provide simulation results that are very consistent with those obtained from the full-order simulations while achieving computational speedups of about an order of magnitude for large-scale cases. These observations indicate that the proposed ROM exhibits satisfactory generalization performance, making it suitable for problems that require many flow simulations under different settings, such as production optimization.
What problem does this paper attempt to address?