Bigoni, Francesco (Eni S.p.A) | Pirrone, Marco (Eni S.p.A) | Trombin, Gianluca (Eni S.p.A) | Vinci, Fabio Francesco (Eni S.p.A) | Raimondi Cominesi, Nicola (ZFOD) | Guglielmelli, Andrea (ZFOD) | Ali Hassan, Al Attwi Maher (ZFOD) | Ibrahim Uatouf, Kubbah Salma (ZFOD) | Bazzana, Michele (Eni Iraq BV) | Viviani, Enea (Eni Iraq BV)
The Mishrif Formation is one of the important carbonate reservoirs in middle, southern Iraq and throughout the Middle East. In southern Iraq, the formation provides the reservoir in oilfields such as Rumaila/West Qurna, Tuba and Zubair. The top of the Mishrif Formation is marked by a regional unconformity: a long period of emersion in Turonian (ab. 4.4 My) regionally occurred boosted by a warm humid climate, associated to heavy rainfall. In Zubair Field, within the Upper interval of Mishrif Formation, there are numerous evidences of karst features responsible of important permeability enhancements in low porosity intervals that are critical for production optimization and reservoir management purposes.
In the first phase, the integration of Multi-rate Production logging and Well Test analysis was very useful to evaluate the permeability values and to highlight the enhanced permeability (largely higher than expected Matrix permeability) intervals related to karst features; Image log analysis, on the same wells, allowed to find out a relationship between karst features and vug densities, making possible to extend the karst features identification also in wells lacking of well test and Production logging information. This approach has allowed to obtain a Karst/No Karst Supervised dataset for about 60 wells.
In the second phase different seismic and geological attributes have been considered in order to investigate possible correlations with karst features. In fact there are some parameters that show somehow a correlation with Karst and/or NoKarst wells: the Spectral Decomposition (specially 10 and 40 Hz volumes), the detection of sink-holes at top Mishrif on the Continuity Cube and its related distance, the sub-seismic Lineaments (obtained from Curvature analysis and subordinately from Continuity), distance from Top Mishrif. In the light of these results, the most meaningful parameters have been used as input data for a Neural Net Process ("Supervised Neural Network") utilizing the Supervised dataset both as a Trained dataset (70%) and as a Verification dataset (30%). A probability 3D Volume of Karst features was finally obtained; the comparison with verification dataset points out an error range around 0.2 that is to say that the rate of success of the probability Volume is about 80%.
The final outcomes of the workflow are karst probability maps that are extremely useful to guide new wells location and trajectory. Actually, two proof of concept case histories have demonstrated the reliability of this approach. The newly drilled wells, with optimized paths according to these prediction-maps, have intercepted the desired karst intervals as per the subsequent image log interpretation, which results have been very valuable in the proper perforation strategy including low porous intervals but characterized by high vuggy density (Karst features). Based on these promising results the ongoing drilling campaign has been optimized accordingly.
This paper discusses the re-construction of the long-term development plan for an offshore giantfield located in Abu Dhabi with the aim to mitigate the rising challenges in the maturing field. The primary objective is to understand the reservoir behavior in terms of fluid movement incorporating the learning from the vast history while correlating with the geological features.
The field has been divided into segments based on multiple factors considering the static properties such as facies distribution, diagenesis, faults, and fractures while incorporating the dynamic behaviors including pressure trends and fluid movements.
On further analysis, various trends have been identified relating these static and dynamic behaviors. The production mechanism for each of the reservoirs and the subsequent sub reservoirs were analyzed with the help of Chan plots, Hall plots and Lorentz plots which distinctly revealed trends that further helped to classify the wells into different production categories.
Using the above methodology the field has been categorized in segments and color coded to indicate areas of different ranking. The green zone indicates area of best interest which currently has strong pressure support and wells can be planned immediately. The wells in this area are expected to produce with a low risk of water and gas. The yellow zone indicates areas of caution where special wells including smart wells maybe required to sustain production. This area showed relatively lower pressure support owing the location of the water injectors and the degraded facies quality between the injectors and the producers. The red zone highlights areas which are relatively mature compared to the neighboring zones and will require new development philosophy to improve the recovery. The findings from this study were used as the basis for a reservoir simulation study using a history matched model, to plan future activities and improve the field recovery.
This study involved an in-depth analysis incorporating the latest findings with respect to the static and dynamic properties of the reservoir. This has helped to classify the reservoir based on the development needs and will play a critical role in designing the future strategies in less time.
A reservoir simulation study of different fishbone well designs performance compared to a base development well design of extra-long maximum reservoir contact (MRC) single lateral wells is presented. The objective is to compare different well design concepts in a waterflood recovery scheme to achieve production target rate and maximize resource value for economic development of an undeveloped tight carbonate reservoir. The studied reservoir is located in a giant offshore oil field in the Middle East and was used as a representation of the different tight reservoirs within the field. It is characterized by poor quality rocks with a permeability trending from 2 – 0.5 md in a SE – NW direction.
The study compromises an assessment of the achievable initial maximum oil rate, volumetric reservoir sweep and expected ultimate oil recovery factor for different well design concepts for a base short well spacing utilized for effective pressure support. In addition to that, the impact of fishbone well design on well count reduction potential utilizing twice the base short well spacing compared to single lateral wells development design utilizing the base short well spacing was evaluated.
A sector model with equal producer to injector ratio was used with refined gridding to wells and bulk area gridded with a cell size of 10 m by 10 m in a representative area of the reservoir. The modeled wells incorporated with vertical flow performance tables with gas lift capabilities. The analysis also incorporated generating streamlines for analyzing fishbone well designs areal reservoir sweep and an examination of remaining movable oil areal distribution. An assessment matrix was formulated for comparing extra-long MRC single laterals base development design versus different fishbone well designs. The assessment matrix incorporated in addition to reservoir related flow performance indicators: drilling complexity and well cost, well life cycle activities, etc. for a comprehensive assessment.
The main findings show that fishbone well designs have complicated areal sweep performance, especially with sealed motherbore, that result in a lower oil recovery factor with higher hydrocarbon pore volume injected and water oil ratio compared to extra-long MRC single laterals. Also, fishbone well designs have serious limitations during well life cycle activities compared to extra-long MRC single lateral design in terms of stimulation, well accessibility and well intervention options making the extra-long MRC single laterals the preferred field development concept within tight reservoirs especially with the base short well spacing. Finally, the analysis has shown that Fishbone well designs can’t reduce the well count since base short well spacing is still needed for effective pressure support by water injection in addition to maximizing the oil recovery factor within the field life time and building and sustaining the target plateau.
In this work we discuss the successful application of our previously developed automated scenario reduction approach applied to life-cycle optimization of a real field case. The inherent uncertainty present in the description of reservoir properties motivates the use of an ensemble of model scenarios to achieve an optimized robust reservoir development strategy. In order to accurately span the range of uncertainties it is imperative to build a relatively large ensemble of model scenarios. The size of the ensemble is directly proportional to the computational effort required in robust optimization. For high-dimensional, complex field case models this implies that a large ensemble of model scenarios which albeit accurately captures the inherent uncertainties would be computationally infeasible to be utilized for robust optimization. One of the ways to circumvent this problem is to work with a reduced subset of model scenarios. Methods based on heuristics and ad-hoc rules exist to select this reduced subset. However, in most of the cases, the optimal number of model realizations must be known upfront. Excessively small number of realizations may result in a subset that does not always capture the span of uncertainties present, leading to sub-optimal optimization results. This raises the question on how to effectively select a subset that contains an optimal number of realizations which both is able to capture the uncertainties present and allow for a computationally efficient robust optimization. To answer this question we have developed an automated framework to select the reduced ensemble which has been applied to an original ensemble of 300 equiprobable model scenarios of a real field case. The methodology relies on the fact that, ideally, the distance between the cumulative distribution functions (CDF) of the objective function (OF) of the full and reduced ensembles should be minimal. This allows the method to determine the smallest subset of realizations that both spans the range of uncertainties and provides an OF CDF that is representative of the full ensemble based on a statistical metric. In this real field case application we optimize the injection rates throughout the assets life-cycle with expected cumulative oil production as the OF. The newly developed framework selected a small subset of 17 model scenarios out of the original ensemble which was used for robust optimization. The optimal injection strategy achieved an average increase of 6% in cumulative oil production with a significant reduction, approximately 90%, in the computational effort. Validation of this optimal strategy over the original ensemble lead to very similar improvements in cumulative oil production, highlighting the reliability and accuracy of our framework.
Yonebayashi, Hideharu (INPEX CORPORATION) | Iwama, Hiroki (INPEX CORPORATION) | Takabayashi, Katsumo (INPEX CORPORATION) | Miyagawa, Yoshihiro (INPEX CORPORATION) | Watanabe, Takumi (INPEX CORPORATION)
CO2 injection is one of widely applied enhanced oil recovery (EOR) techniques, moreover, it is expected to contribute to the climate change from a viewpoint of storing CO2 in reservoir. However, CO2 is well known to accelerate precipitating asphaltenes which often deteriorate production. To understand in-situ asphaltene-depositions, unevenly distributed in composite carbonate core during a CO2 flood test under reservoir conditions, were investigated through numerical modelling study.
Tertiary mode CO2 core flood tests were performed. A core holder was vertically placed in an oven to maintain reservoir temperature and to avoid vertical segregation. A composite core consisting of four Ø1.5" × L2.75" plug cores, which had similar porosity range but slightly varied air permeabilities, was retrieved from a core holder after the flooding test. The remaining hydrocarbon was extracted by Dean-stark method, and heptane insoluble materials were extracted from each plug core via IP-143 method to observe distribution of asphaltene deposits. The variation of asphaltene mass in plug cores was investigated to explain its mechanism thermodynamically.
The core flood test was completed to achieve a certain additional oil recovery by 15 pore volume CO2 injection without any unfavorable differential pressure. The remaining asphaltene mass in each plug core revealed a trend in which more asphaltene collected from the inlet-side core. We assumed a scenario to explain the uneven asphaltene distribution by incorporating the vaporized-gas-drive and CO2 condensing mechanism. Namely, asphaltenes deposited immediately when pure CO2 contacted with oil. The contact between more pure CO2 and oil might be more frequently occurred in inlet-side core. To reproduce the scenario, a cubic-plus-association (CPA) model was generated to estimate asphaltene precipitating behavior as injected gas composition varied. In the first plug core, more pure CO2 gas was considered to contact with fresh reservoir oil compared with the downstream cores which might have less pure CO2 because of its condensation. The light-intermediate hydrocarbon gas vaporized by CO2 was also considered to emphasize the trend of more asphaltene deposits in upstream-side cores. The CPA model revealed consistent phenomenon supporting the scenario.
A flow simulation-driven time-lapse seismic feasibility study is performed for the Amberjack field that leverages existing multi-vintage 4D time-lapse seismic data. The focus is a field consisting of stacked shelf and deepwater reservoir sands situated in the Gulf of Mexico in Mississippi Canyon Block 109 in 1,030 ft of water. The solution leverages seismic interpretation, seismic inversion, earth modeling, and reservoir simulation [including embedded petro-elastic modeling (PEM) capabilities] to enable the reconciliation of data across multiple seismic vintages and forecast the optimal future seismic survey acquisition in a closed-loop. The overarching feasibility solution is integrated and simulation-driven involving multi-vintage seismic inversion, spatially constraining the petrophysical property model by seismic inversion, and performing reservoir simulation with the embedded PEM. The PEM is used to compute P-impedance and Vp/Vs dynamically, which enables tuning to both historical production and multi-vintage seismic data. The process considers a hybrid fine-scale 3D geocellular model in which the only upscaling of petrophysical properties occurs when the P-impedance from seismic inversion is blocked to the 3D geocellular grid. This process minimizes resampling errors and promotes direct tuning of the simulator response with registered seismic that has been blocked to a geocellular earth model grid. The results illustrate a three-part simulation-to-seismic calibration procedure that culminates with a prediction step which leads to a simulation-proposed time-lapse seismic acquisition timeline that is consistent with the calibrated reservoir simulation model. The first calibration tunes the model to historical production profiles. The second calibration reconciles the dynamic P-impedance estimate of the simulated shallow reservoir with that of the seismic inversion blocked to the 3D geocellular grid. The combination of these two steps outline a seismic-driven history matching process whereby the simulation model is not only consistent with production data but also the subsurface geologic and fluid saturation description. Large and short wavelength disparities in the P-impedance calibration existing between the simulator response and the time-lapse seismic data are attributed to resampling errors as a result of seismic inversion-derived P-impedance being blocked to the 3D geocelluar grid, as well as sparse well control in the earth model which leads to the obscuring of some asset-specific characteristics. The results of the third calibration step show how the time-lapse seismic feasibility solution accurately confirms prior seismic surveys undertaken in the asset. Given this confirmation, the solution achieves a suitable prediction of seismic-derived rock property response from the reservoir simulator as well as the optimal future time-lapse seismic acquisition time.
There are a vast number of reservoirs with drill cuttings and core images that have classification problems associated with them. This could be due to the images not being classified in the first place, or the images may be available but the interpretation reports could be missing. Another problem is that images from different wells could be interpreted by different wellsite geologists/sedimentologists and hence result in an inconsistent classification scheme. Finally, there could also be the problem of some images being incorrectly classified. Ergo it would be desirable to have an unbiased objective system that could overcome all of these issues. Step in convolutional neural networks. Advances during this decade in using convolutional neural networks for visual recognition of discriminately different objects means that now object recognition can be achieved to a significant extent. Once the network is trained on a representative set of lithological classes, then such a system just needs to be fed the raw drill cuttings or core images that it has not seen before and it will automatically assign a lithological class to each image and an associated probability of the image belonging to that class. In so doing, images below a certain probability threshold can be automatically flagged for further human investigation. The benefit of such a system would be to improve reservoir understanding by having all available images classified in a consistent manner hence keeping the characterization consistent as well. It would further help to reduce the time taken to get human expertise to complete the task, as well as the associated cost.
Abdul Ghani, Mohamad (IFP Energies nouvelles) | Ayache, Simon Victor (IFP Energies nouvelles) | Batôt, Guillaume (IFP Energies nouvelles) | Gasser-Dorado, Julien (IFP Energies nouvelles) | Delamaide, Eric (IFP Technologies Canada Inc)
Although SAGD is a very popular in-situ extraction method in Canada, this thermal process relies on huge energy and water consumption to generate the steam. Irregular growth of the steam-chamber due to heterogeneities further degrades its yield. Contact between the steam chamber and the overburden also leads to heat losses. The objective of this paper is to investigate how Foam Assisted-SAGD could mitigate these technical issues and improve the efficiency of the SAGD process. Compositional thermal reservoir simulations are used to simulate and analyze a Foam Assisted-SAGD pilot. The shear-thinning effect close to the wells is also accounted for. The simulations are run on a homogeneous model mimicking the Foster Creek project in Alberta, Canada. Several type of injection sequences have been analyzed in terms of foam formation, back-produced surfactants and cumulative Steam-Oil-Ratio. Results are compared with the original SAGD performance. In order to propagate the foaming surfactants throughout the steam chamber the injection sequence needs to be properly determined. A simple continuous Foam Assisted-SAGD injection would lead to an accumulation of surfactant between the wells due to gravity segregation, preventing the foam from acting on the upper part of the steam chamber. Furthermore surfactant production occurs after a few weeks due to the proximity of the producer and the injector. A proper injection strategy of the type SAGD/slug/SAGD/slug is found to delay the chemical breakthrough and increase the amount of surfactant retained in the reservoir while allowing the surfactant propagation throughout the steam chamber. After optimization the Foam Assisted-SAGD process appears to be technically promising.
Al-Farisi, Omar (Khalifa University of Science and Technology) | Zhang, Hongtao (Khalifa University of Science and Technology) | Raza, Aikifa (Khalifa University of Science and Technology) | Ozzane, Djamel (ADNOC) | Sassi, Mohamed (Khalifa University of Science and Technology) | Zhang, TieJun (Khalifa University of Science and Technology)
Automated image processing algorithms can improve the quality and speed of classifying the morphology of heterogeneous carbonate rock. Several commercial products have worked to produce petrophysical properties from 2D images and with less extent from 3D images, relying on image processing and flow simulation. Images are mainly micro-computed tomography (μCT), optical images of thin-section, or magnetic resonance images (MRI). However, most of the successful work is from the homogeneous and clastic rocks. In this work, we have demonstrated a Machine Learning assisted Image Recognition (MLIR) approach to determine the porosity and lithology of heterogeneous carbonate rock by analyzing 3D images form μCT and MRI. Our research method consists of two parts: experimental and MLIR. Experimentally, we measured porosity of rock core plug with three different ways: (i) weight difference of dry and saturated rock, (ii) NMR T2 relaxation of saturated rock, and (iii) helium gas injection of rock after cleaning and drying.
We performed MLIR on 3D μCT and MRI images using random forest machine-learning algorithm. Petrophysicist provided a set of training data with classes (i.e., limestone, pyrite, and pore) as expert knowledge of μCT Image intensity correspondence to petrophysical properties. MLIR performed, alone, each task for identifying different lithology types and porosity. Determined volumes have been checked and confirmed with three different experimental datasets. The measured porosity, from three experiment-based approaches, is very close. Similarly, the MLR measured porosity produced excellent results comparatively with three experimental measurements, with an accuracy of 97.1% on the training set and 94.4% on blind test prediction.
Bhushan, Yatindra (ADNOC Onshore) | Ali Al Seiari, Reem (ADNOC Onshore) | Igogo, Arit (ADNOC Onshore) | Hashrat Khan, Sara (ADNOC Onshore) | Al Mazrouei, Suhaila (ADNOC Onshore) | Al Raeesi, Muna (ADNOC Onshore) | Al Tenaiji, Aamna (ADNOC Onshore)
A reservoir simulation study has been performed to assess the enhanced oil recovery benefits for a proposed pilot on Simultaneous Injection of Miscible Gas (CO2) and Polymer (SIMGAP) in a giant carbonate reservoir (B) in Abu Dhabi. The model has been used to carry out uncertainty analysis for various input parameters and analyze their impact on pilot performance. The paper discusses the uncertainty analysis in detail.
Reservoir-B consists of B_Upper and B_Lower layers which are in full hydrodynamic equilibrium. However, in the southern and western parts of the reservoir, the B_Upper layer has permeabilities that are one to two orders of magnitude higher than the B_Lower layer. The reservoir is on plateau production under waterflooding, however, it is observed that there is water override in B_Upper. The B_Upper layer is being waterflooded very efficiently, while the B_Lower layer remains largely unflooded and forms the key target for enhanced oil recovery (EOR).
The proposed SIMGAP pilot plans to inject polymer into the B_Upper layer and CO2 into the B_Lower layer with producers in the B_Lower layer. The pilot will utilize a line drive pattern at 250m spacing using 3000 ft horizontal wells. There will be two central horizontal injectors (one in B_Upper and the other in B_Lower) and two horizontal producers (one on either side of the central injectors).
Pilot uncertainty analysis cases have been run by varying different parameters that could impact the pilot performance. The parameters that have been varied are polymer viscosity, polymer adsorption, residual resistance factor, thermal stability of polymer, residual oil to miscible flooding (Sorm), residual oil to water flooding (Sorw), Krw end point, high perm streaks, fracture possibility and extension to B_Upper or B_Lower layers, three phase oil relative permeability models, maximum trapped gas saturation, dense zone permeability and pore volume uncertainty. In addition, a grid sensitivity study was undertaken to test the sensitivity of the process to varying levels of dispersion. The results suggest that the key uncertainties which have impact on recovery are polymer viscosity, polymer adsorption, residual oil saturation to water and CO2, presence of high perm streaks and maximum trapped gas saturation values. Vertical observation wells located between the injector and producer wells (equivalent to 0.3 to 0.4 PV of CO2 injection in B_Lower), will be used to confirm whether the SIMGAP process has been successful in containing CO2 in the B_Lower layer and thereby suppressing crossflow.