Alkandari, Dana K. (Australian College of Kuwait) | AlTheferi, Ghaneima M. (Australian College of Kuwait) | Almutawaa, Hawra'a M. (Australian College of Kuwait) | Almutairi, Maryam (Australian College of Kuwait) | Alhindi, Nora (Australian College of Kuwait) | Al-Rashid, Sherifa M. (Australian College of Kuwait) | Al-Bazzaz, Waleed H. A. (Kuwait Institute For Scientific Research)
Formation damage is the impairment of permeability of rocks inside a petroleum reservoir. This occurs during drilling, production, stimulation and enhanced oil recovery operations, by various mechanisms such as chemical, mechanical, biological and thermal. Near wellbore formation damages have a great impact on productivity of the damaged formation. Acidizing is a stimulation method to remove the effect of near wellbore damage through reacting with damaging materials or the formation rocks (carbonate or sandstone rocks) to restore or improve permeability around the wellbore. Several experiments are conducted to study the effect of temperature and acid concentration combined on the efficiency of matrix acidizing. Three different concentrations scenarios of hydrochloric acid (3%, 15%, and 28%) and 4 different temperatures scenarios (25 °C, 35 °C, 70 °C, and 100 °C) were tested to investigate pore-enlargement success effect on permeability. The purpose of this experiment is to introduce the concept of optimized temperature augmented with optimized acid concentration in carbonate matrix acidulation. Morphology of pore geometry and area measurement software is used. New Advancement in imaging that captured pore area enlargement as big-data necessarily for artificial intelligence modeling. Captured pores before treatment and captured pores after thermal-HCL acid treatment have demonstrated that image processing of the actual acidized rock data can select the optimized recipe concentration of acid that will increase permeability, hence recovery. The results show that matrix acidizing is an effective method to improve permeability and enhance production, as it demonstrates that using less acid concentration with the optimized temperature can result in a favorable and satisfying outcomes.
Acid-tunneling is an acid jetting method for stimulating carbonate reservoirs. Several case histories from around the world were presented in the past showing optimistic post-stimulation production increases in open-hole wells, comparing to conventional coiled tubing (CT) acid jetting, matrix acidizing, and acid fracturing. However, many questions about the actual tunnel creation and tunneling efficiency are still not answered. In this paper, the results of an innovative full-scale research program involving water and acid jetting are reported for the first time.
The tunnels are constructed through chemical reaction and mechanical erosion by pumping hydrochloric (HCl) acid through conventional CT and a bottom-hole assembly (BHA) with jetting nozzles and two pressure-activated bending joints that control the tunnel initiation directions. If the jetting speed is too high and the acid is not consumed in front of the BHA during the main tunneling process, then unspent acid flows toward the back of the BHA and creates main wellbore and tunnel enlargement with potential wormholes as fluid leaks off, lowering the tunneling length efficiency.
Full-scale water and acid jetting tests were performed on Indiana limestone cores with 2-4 mD permeability and 12-14% porosity. Many field-realistic combinations of nozzle sizes, jetting speeds, and back pressures were included in the testing program. The cores were 3.75-in. in diameter by 6-in. in length for the water tests, and 12-in. in diameter by 18-in. in length for the tests with 15-wt% HCl acid. The jetting BHA was moved as the tunnels were constructed, at constant force on the nozzle mole, to minimize the nozzle stand-off distance. Six acid tests were performed at the ambient temperature of 46F and two at 97F. The results from the acid tests show that the acid tunneling efficiency can be optimized by reducing the nozzle size and pump rate. The results from the water and acid tests with exactly the same parameters to match the actual CT operations in the field show that the tunnels are constructed mostly by chemical reaction and not by mechanical erosion. The acid tunneling efficiencies obtained from the full-scale acid tests are superior to the average tunneling efficiency of more than 500 actual tunnels constructed during more than 100 acid tunneling operations performed to date worldwide.
The paper describes the full-scale water and acid jetting tests on Indiana limestone cores. The major novelty of this test program consists of performing all measurements with back pressure, unlike all previous water and acid jetting studies reported in literature, to more accurately mimic the downhole well conditions. The novel understanding of the combined effect of the nozzle size, pump rate, and back pressure significantly improves the actual acid-tunneling efficiency.
We demonstrate how key geological uncertainties in a giant onshore carbonate reservoir in the Middle East, most notably fracture permeability and saturation distribution, impact the quality of the history match and change the performance forecasts of a planned Miscible Water Alternating Gas (MWAG) injection process. To achieve this, we used a history matching and multi-objective optimisation (MOO) workflow that was tightly integrated with an innovative reservoir modelling workflow that paid particular attention to the fracture and saturation modelling.
Different geological models for the reservoir were designed by integrating static and dynamic data. These data indicated the need to consider fault-related fractures and to update the saturation distribution in the reservoir model. The effective medium theory was therefore used to estimate effective permeability in order to capture the presence of low-intensity fault-controlled fractures in the reservoir. The integration of Special Core Analysis (SCAL) and log-derived J-functions allowed us to build alternative saturation models that honoured well data with great accuracy. The resulting history matched models therefore accounted for the key geological uncertainties present in the reservoir. Afterwards, MOO was applied for each history matched model to identify well controls that optimally balanced the need to maximise the time on the plateau rate while adhering to the field's gas production constraints.
Our results clearly show that including low-intensity fault-controlled fractures in the reservoir model improved the quality of the history match for the gas oil ratio (GOR), bottom hole pressure (BHP) and water cut. This is especially true for wells located near faults, which were difficult to match in the past. Moreover, our results further show that the updated saturation model improved the quality of the history match for the water cut, particularly for wells located in the transition zone. These different history matched models yielded different production forecasts, where the time at which the reservoir can be produced at the plateau rate varied by up to ten years.
Applying MOO for each history matched model then allowed us to identify well controls for the MWAG injection that could extend the time at which the reservoir would be produced at the plateau rate for up to nine years and the risk of losing production plateau down to two years, while always adhering to the current field operational constraints.
We demonstrate how the integration of MOO with an innovative workflow for fracture and saturation modelling impacts the prediction of a planned MWAG injection in a giant onshore carbonate reservoir. Our work clearly illustrates the potential of integrating MOO with new reservoir characterisation methods to improve the quantification of uncertainties in reservoir performance predictions in carbonate reservoirs.
Numerous carbonate reservoir discoveries were made in Indonesia (
The process involves multiple cycles—from formation evaluation (e.g., geomechanics analysis, design of an effective fracturing method, and production forecasting) through the economic impact to the operator. During the early phase of this integrated study, the uncertainties of all static and dynamic parameters (i.e., geological complexity, rock physics, and stress profile) were considered for fracturing design. Production performances from multiple fracturing stimulation scenarios were then modeled and compared to select the plan that optimizes production for the Berai Formation.
Results demonstrated an effective multidiscipline approach toward a comprehensive strategy to meet the ultimate objective in optimizing production. This project leveraged formation evaluation and fracturing design to deliver integrated solutions from exploration to accurate production forecast. The well stimulations were performed by carefully selecting fluid characteristics based on geological-petrophysical properties, pressure, and stress profiles within the area. Results yielded excellent production gains—for the best case, up to 50% with an average of 40% in comparison with initial production by using an acid that provides optimum fracture geometry and permeability.
This opportunity demonstrated the importance of understanding formation behavior and the parameters that aid the selection of an appropriate fracturing design for a low porosity/permeability carbonate reservoir.
Wei, Chenji (Research Institute of Petroleum Exploration and Development, CNPC) | Zheng, Jie (Research Institute of Petroleum Exploration and Development, CNPC) | Ouyang, Xiaohu (China Petroleum Pipeline Engineering Co., Ltd, CNPC) | Ding, Yutao (China National Oil and Gas Exploration and Development Company Ltd. CNPC) | Ding, Mingming (China National Oil and Gas Exploration and Development Company Ltd. CNPC) | Lin, Shiyao (China National Oil and Gas Exploration and Development Company Ltd. CNPC) | Song, Hongqing (University of Science and Technology Beijing)
Understanding the heterogeneity is critical for a successful water injection in a carbonate reservoir. Thief zone is one of the most obvious forms of heterogeneity, which indicates the thin layer with higher permeability compared to the average reservoir permeability. The existence of thief zone results in earlier water breakthrough and faster water cut increase, which then lead to lower sweep efficiency and smaller recovery factor. Therefore, determining the distribution of thief zone and its impact towards production, and proposing a corresponding development plan are very important.
In this paper, a novel method is established to determine the thief zone distribution based on dynamic surveillance data. A new index is proposed as the relative contribution index to characterize the relative contribution of a certain layer, which is fundamental for thief zone determination. In addition, effect on water flooding development of thief zone's location is studied by experimental and theoretical analysis. The changes of water cut and production rate are analyzed under different conditions such as location of the thief zone, injection rate, and variogram. Finally, optimized development strategy is proposed to deal with the existence of thief zone.
Distribution of thief zone is characterized based on the proposed method, which indicates that thief zone development has intimate relationship with depositional facies and diagenesis. Experimental and theoretical analysis results show that the present model considering stratified water-flood is consistent with the experimental results. The water displacement effect is the best when the thief zone is located in the upper reservoir. This paper also points out the optimal adjustment period for water shutoff and profile control of the reservoir with thief zones. In addition, the greater the injection rate, the faster the water cut increase. Furthermore, the smaller the variogram, the slower the water cut increase, and the later the water breakthrough time.
This study provides a method to characterize thief zone, which can be used as a reference for similar oilfield development. In addition, it provides a quick and reasonable guide in the later adjustment of water flooding development of carbonate reservoirs with thief zones.
A particular challenge inherent to carbonate reservoirs is reservoir rock typing which impacts model initialisation and saturation distributions and hence STOIIP, phase mobilities, and flow behaviours. We explore how flow diagnostics can be used best to detect subtle differences in reservoir dynamics arising from different model initialisations by comparing flow diagnostics simulations with full-physics simulations.
Flow diagnostics are applied to two reservoirs, a synthetic but realistic model representing an analogue for the Arab-D formation and a giant carbonate reservoir from the Middle East. Saturation modelling and reservoir rock typing is based on uniform and heterogeneous Pc and kr distributions, and further employs a state-of-the-art software that integrates of SCAL data and log-derived saturations. Sweep efficiency and dynamic Lorenz coefficients are then derived from the flow diagnostics results to quantify and compare the dynamic behaviour of the reservoir models. The full-physics simulations, which are used to validate the flow diagnostics results, are carried out with a commercial Black Oil simulator.
The flow diagnostics results can clearly distinguish between different homogenous and heterogeneous rock-type distributions, wettability trends, as well as novel saturation modelling approaches that use dedicated software tools. Flow diagnostics capture the same trends in recovery predictions as the full-physics simulations. Importantly though, the total CPU time for a single flow diagnostics calculation including model loading is on the order of seconds, compared to minutes and hours for a single full-physics simulation. These observation give confidence that flow diagnostics can be used effectively to compare and contrast the impact of reservoir rock typing, saturation modelling, and model initialisation on reservoir performance before running full-physics simulations. Flow diagnostic hence allow us to reduce the number of reservoir models from a model ensemble and select a small number of diverse yet realistic reservoir models that capture the full range of geological uncertainties which are then subjected to more detailed reservoir simulation studies.
Flow diagnostics are particularly well suited for complex carbonate reservoirs which are geologically more complex than clastic reservoirs and often exhibit significant uncertainties. Giant carbonate reservoirs are also challenging to simulate using full-physics simulators due to their size, so the impact of geological uncertainty on the predicted reservoir performance is often underexplored. Flow diagnostics are hence an effective complement to quantify uncertainty in state-of-the-art reservoir modelling, history matching and optimisation workflows, particularly for giant carbonate reservoirs.
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.
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.
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.
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