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

•Fluid dynamics of WAG in fractured systems is systematically studied through a simulation work.•The three-phase hysteresis effects in WAG are mostly due to the gas trapping in the larger pores.•Fracture inclination angle does not considerably vary ultimate recovery factor.•Connecting a vertical fracture to a horizontal fracture enhances matrix-fracture communication in WAG.•The maximum recovery factor is achieved at a WAG ratio of one for the present test case. 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 °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 vertica