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.
Contreras Perez, David Rafael (OMV E&P GmbH - Abu Dhabi) | Al Zaabi, Ruqaya Abdulla (ADNOC Offshore - GUL) | Viratno, Bernato (OMV E&P GmbH - Abu Dhabi) | Sellar, Christopher (OMV E&P GmbH - Abu Dhabi) | Susanto, Maria Indriaty (OMV E&P GmbH - Abu Dhabi)
This paper summarizes an efficient workflow for building a reliable static model reference case by improving the accuracy of well placement in a hydrocarbon bearing structure. This is beneficial in optimising upcoming well target position and trajectory planning as well as during the dynamic history matching process. In a non-operated venture, the ability to generate an up-to-date static model that maintains pace with operations, provides valuable insight to advise the operator on the upcoming drilling plan and continuously supports the dynamic model for reserves booking, is highly sought after.
The systematic approach described in this paper is applied to a geo-model from a Middle East carbonate reservoir consisting of over 50 wells with good quality PSDM seismic data. The workflow presented begins with seismic mapping, utilizing volume-based modelling techniques, followed by structural element correction using borehole images (e.g. structural formation dip and true stratigraphic thickness estimate) and finally introduces alternative control points, which enable drilled wellbore trajectories to be structurally anchored, based on layer thicknesses and structural trends within the target reservoir.
Using this approach it is possible to generate a consistent structural model that honours geological markers, measured dip ranges and structural trends seen from seismic data and image logs. During the process one learns more about data quality (e.g. scale of data resolution and depth of investigation), associated with specific fields and carbonate reservoirs through the interaction between geological, geophysical and petrophysical disciplines and ensures their correct use. Data are used to improve the raw interpreted seismic horizons by calibrating mapped thickness distribution against the well tops. 2D visualizations are generated on a well-by-well basis, including map views, curtain sections (along each horizontal well), composite cross-sections and 3D visualizations to show inter-well relationships within different geological layers. As a result the well is placed in the correct structural position. Correct well placement, especially of highly deviated/horizontal wells, provides more accurate identification of reservoir sweet spots, leading to improved well target position and trajectory planning for upcoming wells, and a robust baseline to achieve production/well test history match during the dynamic modelling process.
The identification of the fluid fill history is a necessity for the development strategy of any field, in particular in the Middle East where tectonic history is often reported to affect fluid distribution and contacts in many fields. The fluid fill concept for a low permeability carbonate field has been re-evaluated and modified from a tilted contact interpretation with imbibition of the deepest unit to a field-wide flat contact and primary drainage saturation distribution. The oil volumes in the reservoir under study are sensitive to minor changes in the structure and fluid fill due to the relatively low structural dip and low permeability transitional nature of the reservoir. The paper highlights the importance of removing preconceptions in data analysis and ensuring consistency on interpretations between different available data sources. It also demonstrates how data quality could completely change the fluid fill concept.
The three main reservoir units of the Lower Shuaiba A, Lower Shuaiba B and Kharaib have been charged from two oil migration events. Structural changes post the first primary drainage are revealed by regional seismic images of the shallower horizons. Due to the rock low permeability, the water saturations are above irreducible value and the whole interval is in the "transition zone". Kharaib unit was believed to be imbibed by the aquifer after charge and was not developed. Three possible fluid fill scenarios were investigated: a) tilted contact due to structural changes post-charge, b) imbibition of the deeper interval, c) primary drainage with field-wide flat contact related to the second pulse of charge. Each scenario impacts the development of the three units positively or negatively. Water saturation logs vs. True Vertical Depth plots were the main diagnostic tool used to rule out fluid fill scenarios. The plots were used to recognise lateral changes of the saturation profile and investigate imbibition signatures. Production data were also used to cross check the expected fluid fill scenario. The resistivity tools’ types and mud resistivities were examined.
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.
Identification of tidal channels fairways is key for predicting behavior of areas at higher risk to water breakthrough or otherwise have a significant impact on the development and monitoring of reservoir performance. However, tidal channels in carbonates are not often easily characterized using conventional seismic attributes. It is important to decipher the complexity of the carbonate tidal channel architecture with integrated multisource data and a variety of approaches.
In this paper, petrological characteristics and petrographic analysis is conducted on well logs and validated carefully using core data. Then, the second step is to compare the carbonate channel systems with modern analogue in Bahama tidal flat and outcrop scales in Wadi Mi'Aidin (Northern Oman). Thereafter, the supervised probabilistic neural network (PNN) and linear regression method were undertaken to detect an additional channel distribution.
The relationship of high porosity with low acoustic impedance appeared mostly in the channel facies which reflects good reservoir quality grainstone channels. Outside these channels, the rock is heavily mud filled by peritidal carbonates and characterized by a high acoustic impedance anomaly with low quality of porosity distribution. The new observation of PNN porosity volume revealed a lateral distribution of the Mishrif carbonate tidal channels in terms of paleocurrent direction and the connectivity. Additionally, the prior information from core data and the geological knowledge indicate a good consistency with classified lithology. These observations implied that Mishrif channels consist of a wide range of lithology and porotype fluctuations due to the impact of depositional environment.
The work enables us to provide a new insight into the distribution of channel bodies, and petrophysical properties with quantification of their influence on dynamic reservoir behavior of the main producing reservoir. This work will not only provide an important guidance to the development and production of this case study, however also deliver an integrated work path for the similar geological and sedimentary environment in the nearby oil fields of Southern Iraq.
This work presents a new open access carbonate reservoir case study that uniquely considers the major uncertainties inherent to carbonate reservoirs using one of the most prolific aggradational parasequence carbonate formation set in the U.A.E; the Late Barremian Upper Kharaib Mb. as an analogue. The ensemble considers a range of interpretational scenarios and geomodelling techniques to capture the main components of its reservoir architectures, stratal geometries, facies, pore systems, diagenetic overprints and wettability variations across its platform-to-basin profile.
Fully anonymized data from 43 wells across 22 fields in the Bab Basin, U.A.E from different geo-depositional positions and height above FWL’s (specified to capture multiple structural positions) within an area of 36,000 km2 was used. The data comprises of a full suite of open hole logs and core data which has been anonymized, rescaled, repositioned and structurally deformed; FWL’s were normalized and the entire model was placed in a unique coordinate system. Our petrophysical model captures the geological setting and reservoir heterogeneities of selected fields but now at a manageable scale.
The novelty of this work has been to create semi-synthetic
This paper has an objective of identifying the nature of formation fluid from an extreme tight fractured reservoir. A good understanding of petrophysical properties of the reservoir rock as well as the fluid it contains constitutes a real challenge for tight reservoirs, that are the most common unconventional sources of hydrocarbons. The front-line characterization mean is the Wireline logging which comes directly after drilling the well or while drilling, knowing that for low to extreme low porosity-permeability reservoirs any attempt of conventional well testing will not bring any added value not rather than a confirmation of reservoir tightness. A tailored workflow was adopted to design the most appropriate formation testing module, select the best depths for fluid sampling, and distinguish hydrocarbon from water bearing intervals. This workflow involves ultrasonic and Electric Borehole Images in combination with Sonic Scanner for natural fractures detection, localization and characterization, integrating Dielectric recording and processing for petrophysical evaluation, then Formation Testing was carried out for fluid identification and sampling. The use of borehole electric and sonic imager coupled with advanced sonic acquisition helped not only to identify the natural fractures depths, but also the nature of these fractures. This integration was used for selecting the sampling station.
Albertini, Cristian (Eni Spa) | Bigoni, Francesco (Eni Spa) | Francesconi, Arrigo (Eni Spa) | Lazzeri, Riccardo (Eni Spa) | Vercellino, Alberto (Eni Spa) | Borromeo, Ornella (Eni Spa) | Gabellone, Tatyana (Eni Spa) | Consonni, Alberto (Eni Spa) | Geloni, Claudio (Eni Spa)
The reservoir quality of Karachaganak Carbonates Field results significantly affected by diagenetic processes. In particular, the replacive dolomitization affects porosity, permeability and irreducible water saturation while the precipitation of anhydrite reduces both porosity and permeability. Such impacting processes were therefore analysed and described in the reservoir 3D Model following geologically consistent rules that honour well data.
The field scale diagenetic study was performed following five steps:
Core data studies Lithological logs analysis Hydrological processes identification Hydrological processes reactive transport simulations 3D Lithological model building
Core data studies
Lithological logs analysis
Hydrological processes identification
Hydrological processes reactive transport simulations
3D Lithological model building
The dolomite distribution, estimated from the lithological log analysis and cores data, results mainly confined on the flanks of the paleo-high. This distribution was endorsed by the results of 3D field scale reactive transport modelling related to Kohout geothermal convection mechanism acting in the shallow burial of the carbonate paleo-high at each stratigraphic unit. The final lithological 3D Model was built consistently with this hydrological process calibrated with well data used as verification data set in the stochastic simulations.
The anhydrite distribution, estimated from lithological log analysis and cores data, is, generally, present in a few percentage of volume and, mainly, in the upper section of the reservoir (less than 250 m, below the bottom of the overlaying Kungurian evaporites). This anhydrite was related to diffuse downward percolation of the Kungurian brine and, marginally, to dolomitization. The occurrence of higher concentration of anhydrite was also locally observed but generally connected to fracture infill and, sometimes, also in the deeper section of the reservoir. These events were related to brine percolation exploiting a network of syn-depositional fractures, particularly along the flanks of the carbonate bank (Neptunian dykes). Such hydrological processes was endorsed by 2D reactive transport modelling. In fact, the anhydrite infilling fractures may have a significant impact on the reservoir flow path and therefore a workflow for identification of these Neptunian dykes was applied, based on seismic attributes (Continuity and Curvatures) according to the Eni proprietary workflow utilized for the identification of sub-seismic discontinuities (Tfrac-Sibilla).
The so estimated dolomite distribution represents about the 15% of the lithology at field scale but up to the 60% on the flanks of the carbonate build-up, marginal areas investigated by very few wells but impacting on about the 30% of the field total GBV. Accordingly, the petrophysical characteristics of the field flanks result affected, in the 3D Reservoir Model, by the presence of dolomite, i.e. increased porosity, permeability and irreducible water saturation. Moreover, the identification of the sub-seismic discontinuities filled by anhydrite allows a better description of the permeability baffles affecting the 3D model flow paths.
Full field development of the Upper Jurassic carbonates, offshore Abu Dhabi is exceedingly challenging. The heterogeneous texture, complicated pore systems and intensive lithology changes all mark the regressive cycles of sedimentation. Such complicated characteristics obscure formation evaluation of these formations. Advanced well logging tools and interpretation methodologies are implemented to minimize the petrophysical uncertainties to qualify the products as field development critical elements. This case study highlights a newly applied NMR log interpretation approach. The results help to understand the complex pore system in a tight carbonate layer, along a horizontal drain drilled close to the oil-water contact.
NMR log data was acquired in real-time while drilling simultaneously with Gamma Ray, Resistivity and Image Logs. Earlier field studies recommended swapping standard T2 free fluid relaxation cutoff values by actual laboratory NMR measurements for a higher precision suitable for the reservoir texture heterogeneity, the study itself supported the application of higher cutoff values to better discriminate the free fluid in well-connected macro pores from the irreducible which will have a direct impact on the computed permeability.
In this case study, a variable free-fluid T2 cutoff was firstly implemented based on arbitrary estimations to match the computed Coates permeability to the offset core values. Free-fluid, irreducible fluids were sequentially computed. A unique NMR-Gamma Inversion (NMR-GI) workflow is further utilized as a mathematically defined approach to process the raw data using probabilistic functions. The result is a more precise pore size distribution, coherent with the geological variations. NMR Capillary pressure was computed.
The complex formation texture could be accurately tracked for thousands of feet drilled along the horizontal drain. After validation with offset core, the NMR-GI interpretation was combined with, Archie saturation and Image log analysis for a conclusive assessment. Hydraulic flow units were combined. Successful completion design and production zone selection articulated on the defined open hole log interpretation.
NMR while drilling logging and the applied (NMR-GI) methodology prove to be leading tools to assist in resolving carbonate reservoir complexities. Not only that they help to understand the pore system characteristics, but they effectively support well placement, completion and production.
This study examines which is the margin of usability for Artificial Intelligence (AI) algorithms related to the rock properties distribution in static modeling. This novel method shows a forward modeling approach using neural networks and genetic algorithms to optimize correlation patterns among seismic traces of stack volumes and well rock properties. Once a set of nonlinear functions is optimized in the well locations, to correlate seismic traces and rock properties, spatial response is estimated using the seismic volume. This seismic characterization process is directly dependent on the error minimization during the structural seismic interpretation process, as well as, honoring the structural complexity while modeling. Previous points are key elements to obtain an adequate correlation between well data and seismic traces. The joint mechanism of neural networks and genetic algorithms globally optimize the nonlinear functions and its parameters to minimize the cost function. Estimated objective function correlates well rock properties with seismic stack data. This mechanism is applied to real data, within a high structural complexity and several wells. As an output, calibrated petrophysical time volumes in the interval of interest are obtained. Properties are used initially to generate a geological facies model. Subsequently, facies and seismic properties are used for the three-dimensional distribution of petrophysical properties such as: rock type, porosity, clay volume and permeability. Therefore, artificial intelligence algorithms can be widely exploited for uncertainty reduction within the rock property spatial estimation.