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.
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.
This paper discusses the use of a novel data-driven method for automated facies classification and characterization of carbonate reservoirs. The approach makes an extensive use of wireline and while drilling electrical borehole image logs and provides a direct and fast recognition of the main geological features at multi-scale level, together with secondary porosity estimation. This embodies an unbiased and valuable key-driver for rock typing, dynamic behavior understanding and reservoir modeling purposes in these puzzling scenarios.
The implemented methodology takes advantage of a non-conventional approach to the analysis and interpretation of image logs, based upon image processing and automatic classification techniques applied in a structural and petrophysical framework. In particular, the Multi-Resolution Graph-based Clustering (MRGC) algorithm that is able to automatically shed light on the significant patterns hidden in a given image log dataset. This allows the system to perform an objective multi-well analysis within a time-efficient template. A further characterization of the facies can be established by means of the Watershed Transform (WT) approach, based on digital image segmentation processes and which is mainly aimed at quantitative porosity partition (primary and secondary).
The added value from this data-driven image log analysis is demonstrated through selected case studies coming from vertical and sub-horizontal wells in carbonate reservoirs characterized by high heterogeneity. First, the MRGC has been carried out in order to obtain an alternative log-facies classification with an inherent textural meaning. Next, the WT-based algorithm provided a robust quantification of the secondary porosity contribution to total porosity, in terms of connected vugs, isolated vugs, fractures and matrix contribution rates. Finally, image log-facies classification and quantitative porosity partition have been integrated with production logs and pressure transient analyses to reconcile the obtained carbonate rock types with the effective fluid flows and the associated dynamic behavior at well scale.
The presented novel methodology is deemed able to perform an automatic, objective and advanced interpretation of field-scale image log datasets, avoiding time-consuming conventional processes and inefficient standard analyses when the number of wells to be handled is large and/or in harsh circumstances. Moreover, secondary porosity can be proficiently identified, evaluated and also characterized from the dynamic standpoint, hence representing a valuable information for any 3D reservoir models.
This paper presents a diagnostic workflow to understand and implement rock and fluid modeling in a diagenetically heterogeneous and hydrodynamically pressured Middle East carbonate field. The workflow allows interactive field data integration, provides guidance for reservoir property distribution and fluid contact generation in order to improve reserves and forecasting estimation. The workflow is useful to a reservoir modeler in QA/QC role and in this case it proves particularly applicable in an organization with constrained resources during the farm-in process. The workflow runs on numerical methods within the static model to avoid database discrepancy during the diagnostic process. Using the core (CCAL, SCAL), log and pressure database, the geoscientist can assess subsurface modeling outputs from the simplest to more complex deterministic scenarios. The process aims to minimize the discrepancy between data input and model output while continuously honoring the data, maintaining realistic correlations (e.g. between static permeability and water saturation) and respecting inherent uncertainty.
Using a data-rich Middle East carbonate reservoir, the pre- and post-diagnostic comparison of 3D modeled reservoir properties to the input data are demonstrated. Diagnostic steps have helped to understand potential subsurface scenarios and thus minimize the discrepancy post exercise. The value of the workflow is its ability to pinpoint the key uncertainties in rock and fluid modeling from the field’s vast dataset in a shorter diagnostic time. The application of the workflow in this carbonate reservoir case study increases the importance of geological and property driven rock type classification and its 3D distribution in matching the water saturation profile. This proved particularly challenging in this case study due to the field’s compartmentalization - fluid contact scenario.
Hydrocarbon in place volumes are often inaccurate as a result of poor representation of the reservoir structure (by means of a 3D grid), that in combination with the use of traditional saturation calculation methods, lead to erroneous hydrocarbon volumes and poor investment decisions.
Traditionally a reservoir model is represented with a 3D grid, in a complex setting such as fault intersections and stacked reservoirs. A corner point grid is often used, which has limitations to represent this complexity. Further, the hydrocarbon saturations are then derived on a cell by cell basis on that 3D grid using simple averaging techniques of saturation height functions. The poor structure representation on the pillar grid in addition to the simplistic averaging methods lead to inaccuracies of the in place volumes especially where a prominent transition zone is present.
This paper presents new advanced saturation averaging methods (volume and height weighted) using saturation height functions on 3D grids. The new advanced saturation averaging methods are used on different reservoir models to compare the saturation distribution and volumetric differences against the traditional saturation calculation methods. A 4-way dip closure reservoir model with a tilted free water level (typical example of a carbonate reservoir in the Middle East), and a faulted S-grid model of the F3-FA field (North Sea) are used.
For the 4-way dip closure reservoir model, when comparing the advanced ‘volume weighted’ and traditional ‘by center of the part of the cell’ saturation averaging methods, a significant difference in the water saturations is observed which leads to about 5% difference in the calculation of in place hydrocarbon volumes. Further, it is observed that changing the thickness and orientation of the 3D grid cells can result in even larger differences of 5-10%.
The faulted F3 model shows that the difference between the hydrocarbon saturation values is largest where it matters most, that is, around the fluid contacts and in the transition zone. The new advanced saturation averaging methods give accurate hydrocarbon saturations irrespective of the size or complexity of the 3D grid and without any discretization effects.
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.
In this paper we describe how incorporating inter-well flow of gas tracers into a numerical simulation model of an oil producing field allowed us to improve the reservoir characterization. Gas tracers evidenced the lateral and vertical reservoir connectivity, identified preferential flow paths and eventually provided an additional tool for the dynamic history match.
The field produces light oil with the support of miscible gas injection, from a reservoir composed of two stacked fluvial sandstones units. To improve the reservoir characterization four inter-well gas tracer campaigns involving a total of 13 injector wells have been successfully completed. All these tracer injections have then been modeled at our full field numerical simulation model with the purpose of challenging the reservoir description in it. Input from this exercise have been later-on used during the construction of geomodels, benefiting from our improved reservoir knowledge.
The reservoir is composed of two units, Middle (M) and Upper (U), both deposited as laterally amalgamated fluvial channels. Both units appear as vertically separated by an impervious shale interval, only absent in 2 wells out of 30.
The shale interval was originally considered as a possible flow barrier by the earliest geological models. However, numerical simulation models were only able to replicate observed tracers arrivals when specific vertical connections existed between both units, indicating the shale interval was no as laterally continuous as formerly suspected.
History matching the tracers arrivals in the western field area was also helpful to reveal a fast gas breakthrough between wells which were aligned perpendicularly to the main channel orientation (NW-SE). This finding confirmed that this prevailing channel orientation was not the only responsible for a good reservoir communication, but also the lateral amalgamation (SW-NE) of channels was exerting a significant control in low sinuosity fluvial systems, as well as secondary flow directions in high sinuosity systems.
The improved reservoir characterization have been reflected in subsequent reservoir geological models and numerical simulators, avoiding misleading history matching solutions. Also, it is worth to note that this have had a direct impact on the gas injection strategy followed in the field.
Quite commonly interwell tracer flow is not fully incorporated into numerical simulators. This is a time consuming process, which, in addition to the conventional model uncertainties, requires sensitivities to the associated tracer parameters. This paper demonstrate through a real case how valuable this additional effort may result, and how it may improve the geological and dynamic understanding of our field.
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.
For oil and gas companies, accelerating the time to first hydrocarbon is a strategic objective. Special core analysis programs for tight gas reservoirs may take many months because of the long equilibration times involved in the tests. This represents a bottleneck for achieving the goal of reducing the time-to- market. Both log interpretation and reservoir modelling activities are impacted by the long SCAL durations. In order to face the challenge, a suite of fast methods have been developed. They are fast because they operate under non-equilibrium conditions. The methods give the m&n parameters for electric log interpretation, the endpoint gas relative permeability and the relationship linking initial gas saturation, trapped gas saturation and endpoint water relative permeability in a couple of days.
Reservoir A is being developed in early and interim phases in order to gather static & dynamic data to minimize the risk associated to subsurface uncertainties. In early and interim phases, only production is taking places. During full field, water injection scheme will be implemented using mainly 5-spot pattern. It is very crucial to measure the subsurface uncertainties and their impact on the reservoir development. For this purpose, the uncertainty parameters are identified and their ranges are selected based on the current well performances during probabilistic History matching (PHM) phase. In full field runs, the uncertain subsurface parameters are quantified to prioritize the future reservoir monitoring and data gathering plans. Note that wells are equipped with the permanent downhole pressure gauges.
Reservoir A is one of the major reservoirs of a green-field located offshore Abu Dhabi and is being developed with a 5-spot water injection pattern. The producers and water injectors are horizontal wells which are drilled across different flow unit within the reservoir. The reservoir properties are variable across all the flow units, which may results in a non-uniform water front. Being a green field, there are more uncertainties as compared to the brown field. More than three years production & pressure data is available which is used in this uncertainty study. This production data is mainly used to achieve the probabilistic History match on well-wise basis. In this uncertainty study, previous HM parameters are removed. However, based on previous history matching learnings, the subsurface uncertain parameters ranges are selected for this probabilistic History match phase. The criteria for filtering the valid runs during this phase are set to be ±150 Psi compared to the actual downhole pressure readings. In case of decreasing this filtering range to 75 Psi, results in reduction in the reserve range in P90 to P10. Based on ±150 Psi principle, the subsurface parameter ranges are furthered reformed for full field uncertainty study/run. The industry standard workflow is followed to quantify the subsurface parameters during this phase. In this study, we used the Permeability modifiers based on RRT, Faults transmisibilities, Relative Perm curves (based on SCAL data), Kv/Kh ratio (from PTA), etc. as uncertain parameters. The impact of each parameter is measured and quantified with respect to plateau and total reserves.