Computational fluid dynamic simulation of multi-phase flow in fractured porous media during water-alternating-gas injection process

Naturally fractured reservoirs (NFRs) are among major oil and gas producing reservoirs that are commonly targeted for various enhanced oil recovery (EOR) techniques, particularly water-alternating-gas (WAG) injection. Production from NFRs is complicated due to the flow communication between the matrix and fractures in such porous media. The implementation of WAG injection in NFRs features inherent complexities not only related to the three-phase flow, saturation history, and cycle-dependent hysteresis of the individual phases, but also the fracture-matrix communication, fingering, and early breakthrough of the injected phases particularly during gas injection processes. This paper provides the key results on the computational fluid dynamics (CFD) simulation of WAG injection in a fractured system. We evaluate the impacts of hysteresis, fracture characteristics (aperture, orientation, and fracture density in the network), and the three-phase relative permeability of available phases during the WAG injection. The CFD model simulates an immiscible WAG injection, and the CFD results are compared to the experimental data in a strongly water-wet sand-pack. Similar to the experiments, we simulate Maroon crude as the oil phase, and synthetic brine, and pure CO2 at 100 degrees C and atmospheric outlet pressure. The CFD results are in excellent agreement with the experimental data. The absolute relative error is less than 12% for predicting the ultimate oil recovery factors (RFs) in water flooding (WF) and gas injection (GI) cycles. Including the three-phase hysteresis significantly increases the accuracy of the WAG process simulation, while excluding hysteresis underestimates the instantaneous RFs at each cycle (especially GI cycles) and the ultimate RF by 4%. We also analyze the fracture pattern and configuration; adding fractures to the system increases the system effective permeability, leading to more contact between the injecting fluid and trapped oil, and consequently higher oil recovery. Connecting a vertical fracture to the horizontal fracture enhances the recovery through strong connections between the vertical and horizontal blocks. An increase in the fracture aperture from 0.5 to 3 mm increases the RF from 50% to 59%. The fracture inclination angle from 30 degrees to 90 degrees increases the ultimate RF by only 2%. Including gravity forces in vertical model orientation results in overall improvement in RF through engaging both matrix and fracture media in all cycles. As the permeability contrast between the matrix and fracture media decreases, the flow communication between the two regions increases, which improves the recovery performance of the WAG process. Our simulation results can help to further understand the effect of various parameters and operational conditions on WAG injection performance in fractured reservoirs.