Model-Based Petroleum Field Management in Three Stages: Life-Cycle, Short-Term, and Real-Time

The objective of this work is to present a new practical methodology to manage petroleum fields considering three stages (life-cycle, short-term, and real-time) that can run alongside different model fidelities and characteristics. The model-based field management process follows the general methodology proposed by Schiozer et al. (2019) with four activities: (1) fit-for-purpose models construction, (2) data assimilation for uncertainty reduction, (3) life-cycle production optimization and (4) short-term optimization for real-time implementation. The selection of the production strategy for field management comprehends the last two activities. Life-cycle optimization is the first stage of the process and generates control setpoints for short-term analysis. Short-term optimization is then used to improve the quality of the solutions considering the control parameters of the next cycle (considering a closed-loop procedure). Real-time solution is then implemented considering operational disturbances from real operations. The methodology was applied to a benchmark case (UNISIM-IV-2026) which is a case based on a typical carbonate field from the Brazilian Pre-salt, with light oil and submitted to Water-Alternate-Gas injection with CO2 (WAG-CO2). The results show that the methodology is applicable to real and complex fields. As the three stages can run simultaneously, one can (1) use different model fidelities to improve the quality of the solutions and (2) use model-based solutions for real-time implementation. Life-cycle optimization using complex simulation models and long-term objectives can run in the background to generate control setpoints for short-term analysis in which lower fidelity models and simplified solutions can be used for the control and field revitalization parameters of each closed-loop cycle. Real-time solutions can be implemented considering operational problems and disturbances. This work presents a novel procedure to integrate three stages for production optimization that can run in parallel, allowing the integration of life-cycle and real-time solutions. The methodology (1) allows the use of complex reservoir simulation models from the life-cycle production strategy optimization, (2) focuses short-term control parameters that improve the quality of the short-term solution, and (3) guides real-time implementation, so it can be the basis to a digital field management.

Numerical Study on the Impact of Advanced Phenomena in a Fractured Carbonate Reservoir Subjected to WAG-CO2 Injection

Advanced phenomena related to water-alternating-gas (WAG) injection are usually neglected in numerical simulations. This work evaluates the impact of different physical phenomena on field indicators, considering a typical pre-salt carbonate reservoir (UNISIM-II-D-CO, a dual-por dual-perm compositional case) subjected to WAG-CO2 injection. Additionally, the computational cost incurred by each of these phenomena is evaluated, since it represents a great challenge in optimization and probabilistic studies. The following phenomena are evaluated considering a nominal base case: (i) matrix-fracture transfer calculation, (ii) relative permeability hysteresis, (iii) CO2 and CH4 solubilities in aqueous phase, (iv) diffusion, (v) numerical dispersion control models, and (vi) velocity-dependent dispersion. CO2 and CH4 solubilities in the aqueous phase, as well as molecular diffusion, did not have a significant impact on field indicators, but they increased simulation runtime more than two times. Matrix-fracture transfer modeling was the most impactful factor, followed by hysteresis and velocity-dependent dispersion. Therefore, the impact of these phenomena was also investigated in a probabilistic approach, considering an ensemble of 197 geostatistical scenarios under uncertainty. Risk curves revealed that the advanced matrix-fracture transfer models improve sweep efficiency. This effect is mainly due to gravity force which acts as a driving mechanism for the oil moving from the matrix to fractures. The capillary effect, in turn, was small compared to gravity. The impact of dispersion and hysteresis on risk curves were smaller than the effect of matrix-fracture transfer modelling. However, these phenomena are particularly interesting in UNISIM-II-D-CO due to the presence of Super-K facies. Hysteresis, when applied to low and high permeability layers, reduced gas mobility and, consequently, the gas produced, contributing to the NPV for most models under uncertainty. On the other hand, the velocity-dependent dispersion mainly affected fluid flows in the regions adjacent to Super-K layers, promoting better oil recovery. The inclusion of advanced phenomena related to WAG-CO2 injection can hold importance when modeling fractured carbonate fields, like those found in the Pre-Salt in Brazil. Nevertheless, computational costs might make their inclusion impractical in full-field simulation models employed for optimization and probabilistic studies. In such cases, it is recommended to assess low-fidelity models or alternatives to accelerate simulations, focusing mainly on the most impactful phenomena related to WAG-CO2 injection.

MODEL-BASED ANALYSIS TO EVALUATE THE EFFECT OF POLYMER PROPERTIES AND PHYSICAL PHENOMENA ON POLYMER FLOODING OPERATIONS

Polymer flooding, known for enhancing heavy oil displacement, may encounter efficiency-limiting physical phenomena. This work assesses the impact of key factors like viscosity, shear rate, and adsorption in field applications using a model-based approach on the EPIC001 case, a real offshore Brazilian heavy oil reservoir. In simulation results, increasing the displacing fluid viscosity using apparent or zero-shear functions, oil production and economic returns show improvements over waterflooding. The shear rate effect slightly increases oil production and enhances injectivity loss due to shear-thinning in polymer flow. However, modeling it increases computation costs, as it extends simulation run time from minutes to days, making it impractical for intensive processes like production optimization. An analysis of the method’s effectiveness shows that it varies based on the adsorption level considered. At its highest value, even with a higher oil recovery factor, the economic return was lower than using waterflooding. Combining shear rate and adsorption has a minimal impact on field indicators when compared to adsorption alone. This work enhances the comprehension of physical phenomena and non-Newtonian behavior in tertiary polymer flooding on heavy oil reservoirs and its impact on the production forecast. It also highlights important considerations for modeling-based procedures.