Reservoir

Estimating Residual Oil Saturation in Carbonate Rocks: A Combined Approach of Direct Simulation and Data-Driven Analysis

This paper presents a novel modeling framework for predicting residual oil saturation in carbonate rocks. The proposed framework uses supervised machine learning models trained on data generated by pore-scale simulations and aims to supplement conventional coreflooding tests or serve as a tool for rapid residual oil saturation evaluation of a reservoir.

221498_hero.jpg
Source: Paper SPE 221498

Estimating residual oil saturation (Sor) post-waterflooding is critical for selecting enhanced oil recovery strategies, further field development, and production prediction. We established a data-driven workflow for evaluating Sor in carbonate samples using microcomputed tomography (μ-CT) images.

The two-phase lattice Boltzmann method facilitated the flooding simulation on 7,192 μ-CT samples. Petrophysical parameters (features) obtained from pore network modeling and feature extraction from μ-CT images were utilized to develop tree-based regression models for predicting Sor. Petrophysical features include porosity, absolute permeability, initial water saturation, pore size distribution, throat size distributions, and surface roughness distribution.

Our method excludes vugs and macro/nanoporosity, which complicates multiscale simulations—a recognized challenge in modeling carbonate rocks. When subdividing the image into numerous subvolumes, certain subvolumes may contain vugs exceeding the dimensions of the subvolume itself. Hence, these vugs were omitted given the entirety of the image constitutes a vug. Conversely, vugs with dimensions smaller than those of the subvolume were not excluded. Despite scale limitations, our subsampling, supported by substantial data volume, ensures our microscale porosity predictions are statistically reliable, setting a foundation for future studies on vugs and nanoporosity’s impact on simulations.

The results show that features obtained from dry-sample images can be used for data-driven Sor prediction. We tested three regression models: gradient boosting (GB), random forest, and extreme gradient boosting. Among these, the optimized GB-based model demonstrated the highest predictive capacity for Sor prediction (R2 = 0.87, mean absolute error = 1.87%, mean squared error = 0.12%).

Increasing the data set size is anticipated to enhance the models’ ability to capture a broader spectrum of rock properties, thereby improving their prediction accuracy. The proposed predictive modeling framework for estimating Sor in heterogeneous carbonate formations aims to supplement conventional coreflooding tests or serve as a tool for rapid Sor evaluation of the reservoir.


This abstract is taken from paper SPE 221498 by A. S. Rizk, M. Tembely, W. AlAmeri, E. W. Al-Shalabi, R. Farmanov, and S. Markovic, Khalifa University of Science and Technology. The paper has been peer reviewed and is available as Open Access in SPE Journal on OnePetro.