Integrated Multi-Scale Pore Characterization of Carbonate Rocks in the Barra Velha Formation, Santos Basin, Brazil

Carbonate rocks feature heterogeneous porous systems that span multiple scales, from pore level to the reservoir scale. The complexity and diversity of carbonate reservoirs demand a consistent approach to their characterization. The efficient integration of multiscale imaging data and petrophysical data is increasingly important to address the challenges associated with these complex carbonate reservoirs. A crucial step in overcoming these scale gaps in reservoir modeling and simulation involves enhancing the characterization of reservoir flow units and their associations with geological and petrophysical heterogeneities at varying scales. In this study, we focus on the classification of pore types using digital rock analysis and petrophysical evaluation of pre-salt lacustrine carbonates from the Barra Velha Formation (BVF) in the Santos Basin using computerized tomography (CT), core samples description, and petrography. Eight types of pores were identified at the core scale: interparticle, stratiform-vuggy, growth framework, vuggy, vuggy-fracture, fracture, interclast, and intraclast. The distribution and characteristics of these pore types were analyzed at different scales, including thin-sections and micro-CT, and nuclear magnetic resonance (NMR), which highlights the diversity in the porous system and the impact of different pore types on porosity and permeability. NMR analyses illustrated the pore size heterogeneity to provide distinction between tight and porous samples. Hydraulic rock units (HRUs) were defined based on flow zone indicator (FZI) using the probability plot approach. Seven HRUs were defined: HRU1 and HRU2 represent samples with the highest FZI and rock quality index (RQI) values, whereas HRU3 and HRU4 denote intermediate values. HRU5, HRU6, and HRU7 represent units with the lowest values. HRU1 and HRU2 were predominantly associated with vuggy, growth framework, and interparticle porosities, which are often enhanced by dissolution processes. Conversely, HRUs with reduced reservoir qualities (5, 6, and 7), characterized by the lowest permeability values, are more prevalent in intervals with higher silicification and silica and dolomite cementation, presenting a variety of pore types at a macroscale. The integration of multiscale imaging techniques and petrophysical data underscores the complexity of pore systems, providing crucial insights into their reservoir characteristics.

The influence of major faults and fractures on the development of non-matrix porosity system in a pre-salt carbonate reservoir, Santos Basin – Brazil

Faults and fractures are central for characterizing the permeability distribution in carbonate reservoirs since they act as pathways for diagenetic fluids that often favor intense rock dissolution and permeability. Usually, high permeability zones and fractures are not easily recognized in seismic data due to limited resolution and they are often associated with higher concentrations of hydrocarbons or even significant fluid losses during drilling, thus creating a challenge for hydrocarbon exploration. In the Santos Basin, southeast Brazil, the pre-salt carbonate reservoirs from the Barra Velha Formation (BVE) are the main hydrocarbon producers in Brazilian Atlantic margin and well-known for being extremely heterogeneous, exhibiting complex dual-porosity systems. In this study, we built a conceptual model of these fracture zones and non-matrix porosity formation that helped narrowing the understanding of these complex systems. The characterization of faults and fractures was carried out using seismic, well-logs, and borehole image data to understand the influence of these structures in the porosity formation along the Barra Velha Formation. In the study area, three fault sets were defined (F1, F2, and F3) from seismic data. F1 represents to the larger faults, while the F3 fault set represents the smaller faults related to the reactivation of F1; both sets being oriented NE-SW. The F2 fault set is associated with the rift formation and is oriented to NNE-SSW. These three fault sets compartmentalized the studied area into different domains, each exhibiting distinct fracture sets. The natural open fractures were formed during the reactivation of rift faults and are oriented mainly NW, NNE-NNW, NE, and ENE and were identified across the entire study area, but with different intensity values. The fracture intensity closely relates to the distance from major faults where the wells with the highest fracture intensity occurs located 150–590 m from the larger F1 fault set. Scan-lines were conducted throughout the area to determine the fault width, and an average value of 1.2 km was established. Seven non-matrix porosity classes were characterized and classified into stratigraphically concordant and discordant non-matrix pore types at well scale through borehole image interpretation. The Barra Velha Formation exhibit higher occurrence of stratigraphically discordant non-matrix porosity related to fractured zones while stratigraphically concordant non-matrix porosity is mainly controlled by the paleotopography of the study area. Overall, non-matrix porosity formation tends to follow an orientation that suggests a preferential dissolution flow towards NE and ENE directions. Intervals with higher silica content shows a positive correlation with both fracture intensity and stratigraphically discordant non-matrix porosities. This work provides a conceptual model about the fractures and non-matrix porosity distribution in pre-salt carbonate rocks that can help address some of associated structural and stratigraphic uncertainties during field appraisal and development.

Integration between experimental investigation and numerical simulation of alkaline surfactant foam flooding in carbonate reservoirs

In Brazil, pre-salt carbonate reservoirs are largely responsible for the current increase in oil production. However, due to its peculiar characteristics, increasing oil recovery by water injection is not enough. Therefore, we seek to evaluate the recovery potential using chemical methods (cEOR). Among these, the Alkali Surfactant Foam (ASF) method appears with high potential, a variant of Alkali Surfactant Polymers (ASP) without the problems presented by it. Therefore, this work presents an innovative methodology, which seeks to evaluate the potential for recovery with ASF in carbonate reservoirs by integrating experimental characterization and recovery prediction using reservoir simulation. For this, phase behavior and adsorption analyses were carried out. The experimental results provided key parameters for the simulation, such as optimal salinity, surfactant adsorption, foam mobility reduction factors. The results are from two case studies of AS and ASF flooding, using a section of UNISIM-II benchmark, using a one-quarter of five-spot model. Having the modelling for these cEOR methods defined, an optimization process for each method was applied, allowing a reliable comparison among the methods and over a base case of water injection, seeking the maximization of the net present value (NPV). As a result, in the experimental part, a low interfacial tension (IFT) value of 0.003 mN/m was achieved with a surfactant adsorption reduction of 17.9% for an optimal setting among brine (NaCl), alkali (NaBO2.4H2O), and surfactant (BIO-TERGE AS 40). In the reservoir simulation part, using a fast genetic algorithm in the optimization process, a NPV of US$ 14.43 million higher than the base case (water injection) and a 4.5% increase in cumulative oil production for the ASF injection case were obtained. Considering the analyses of production curves (cumulative oil production and oil rate) and oil saturation maps, a considerable oil production anticipation was observed, which was the main reason for NPV improvement, proving the high potential for application of the ASF method in carbonate reservoirs.

Paleokarst features in the Aptian carbonates of the Barra Velha Formation, Santos Basin, Brazil

Seismic data, borehole image logs, and conventional well logs were used to investigate the distribution and characteristics of paleokarst features in the Aptian carbonates of the Barra Velha Formation in a pilot area of the Santos Basin, Brazil. Multiple seismic attributes were used to enhance details on the seismic data and highlight key seismic parameters including strata deformation and geometry, continuity of seismic events, and fault patterns. The study found that karst structures are controlled by faults and fractures along structural highs, which served as a conduit for the flow of dynamic fluids that dissolved the carbonate materials. Several closed, circular depressions and bright spots identified in the northeastern portion of the study area represent possible sinkhole structures. Epigenic and hypogenic processes due to the action of meteoric water, hydrothermal activity, and intra-formation acidity along regional unconformities in the Barremian-Aptian may have been responsible for the dissolution. Limitations of this study are related to the difficulty of integrating multiple datasets with various scales. However, the higher confidence for the occurrence of the karst features is provided by borehole images at the sub-seismic scale. The findings of this study hold significant relevance for the strategic planning of energy development and carbon sequestration initiatives in the Brazilian continental margins, thereby aiding in informed decision-making.

Comparing WAG-CO 2 injection with continuous water and gas injection in separate wells for the development and management of a CO2-rich light oil fractured carbonate reservoir subject to full gas recycling

Many projects in the Brazilian pre-salt assume the use of water alternating gas (WAG-CO2) injection as an ecologically safe carbon storage strategy, with improved hydrocarbon recovery. However, studies that compare these advantages with a simpler management plan are not common. The objective of this work is to compare WAG-CO2 injection with continuous injection of water and gas (CIWG) rich in CO2 in separate wells for the development and management of a light-oil fractured carbonate reservoir subject to full gas recycling. We employed the UNISIM-II benchmark model, a naturally fractured carbonate reservoir with Brazilian pre-salt characteristics, which enables an application in controlled environment where the reference response is known (UNISIM-II-R). We used a model-based decision analysis for production strategy selection, hierarchical optimization of the decision variables and algorithms to maximize the objective function. Representative models (RM) are selected from the ensemble of models and used to incorporate the effects of geological, reservoir, and operational uncertainties into the optimization process. The net present value is the objective function during the nominal optimization of candidate strategies of each RM and the expected monetary value and risk analysis are considered to select the final production strategy considering uncertainties. The risk analysis was quantified based on downside risk and upside potential relation to a benchmark return. We optimized two alternative development plans (one considering WAG-CO2 injection and the other continuous injection of water and gas in separate wells) and compared their performance indicators and decision variables, including design variables (number, type and placement of well, and size of production facilities) and life-cycle control rules (management of equipment over time). We then applied a cross-simulation, where the best strategy optimized for one recovery method was applied to the other and the injection strategy was optimized again. We were therefore able to assess the need to pre-define the recovery method before defining design variables to validate the flexibility of each strategy for possible future changes in the recovery mechanism. Finally, we repeated the study for different reservoir scenarios to compare the alternatives considering typical uncertainties of the Brazilian pre-salt and validated the final strategies in the reference model to quantify the real value in decision making. The strategies reached a full gas recycling in both recovery methods and allowed a comparison of their advantages and disadvantages. The operations of WAG-CO2 injection can be more complex and the equipment more expensive. The novelty of this work is the consideration of continuous injection of water and gas in separate wells as a simpler alternative to the development and management of pre-salt oil fields, since this method may also meet operators’ and environmental demands, bearing simpler operating challenges and promoting good recovery and profitability.

Construction of Single-Porosity and Single-Permeability Models as Low-Fidelity Alternative to Represent Fractured Carbonate Reservoirs Subject to WAG-CO2 Injection Under Uncertainty

Fractured carbonate reservoirs are typically modeled in a system of dual-porosity and dual-permeability (DP/DP), where fractures, vugs, karsts and rock matrix are represented in different domains. The DP/DP modeling allows for a more accurate reservoir description but implies a higher computational cost than
the single-porosity and single-permeability (SP/SP) approach. The time may be a limitation for cases that require many simulations, such as production optimization under uncertainty. This computational cost is more challenging when we couple DPDP models with compositional fluid models, such as in the case of fractured light-oil reservoirs where the production strategy accounts for water-alternating-gas (WAG) injection. In this context, low fidelity models (LFM) can be an interesting alternative for initial studies. This work shows the potential of compositional single-porosity and single-permeability models based on pseudo-properties (SP/SP-P) as LFM applied to a fractured benchmark carbonate reservoir, subject to WAG- CO2 injection and gas recycle. Two workflows are proposed to assist the construction of SP-P models for studies based on (i) nominal approach and (ii) probabilistic approach of reservoir properties. Both workflows begin with a parametrization step, in which the pseudo-properties are optimized for a base case in order to minimize the mismatch between forecasts of the SP/SP-P and DP/DP models. The new parametrization methods proposed in this work showed to be viable for the construction of the SP/SP-P models. For studies under uncertainties, the workflow proposes obtaining pseudo-properties by robust optimizations based on representative models from a DP/DP ensemble, which proved to be an effective method. The case study is the benchmark UNISIM-II-D-CO with an ensemble of 197 DP/DP models and two different production strategies. The risk curves for production, injection and economic indicators obtained from DP/DP and SP/SP-P ensembles showed good match and the computational time spent on simulations of the SP/SP-P ensemble was 81% faster than DP/DP models, on average. Finally, the responses obtained from both ensembles were validated in a reference model (UNISIM-II-R) that represents the true response and is not part of the ensemble. The results indicate the SP/SP-P modeling as a good LFM for preliminary assessments of highly time-consuming studies. Besides, the workflows proposed in this work can be very useful for assisting the construction of SP/SP-P models for different case studies. However, we recommend the use of the high-fidelity models to support the final decision.