Al-Maqtari, Ameen N. (SAFER E&D Operations Company) | Saleh, Ahmed A. (SAFER E&D Operations Company) | Al-Haygana, Adel (SAFER E&D Operations Company) | Al-Adashi, Jaber (SAFER E&D Operations Company) | Alogily, Abdulkhalek (SAFER E&D Operations Company) | Warren, Cassandra (Schlumberger) | Mavridou, Evangelia (Schlumberger) | Schoellkopf, Noelle (Schlumberger) | Sheyh Husein, Sami (Schlumberger) | Ahmad, Ammar (Schlumberger) | Baig, Zeeshan (Schlumberger) | Teumahji, Nimuno Achu (Schlumberger) | Thiakalingam, Surenthar (Schlumberger) | Khan, Waqar (Schlumberger) | Masurek, Nicole (Schlumberger) | Andres Sanchez Torres, Carlos (Schlumberger)
A 3D petroleum systems model (PSM) of Block 18 in the Sab'atayn basin, onshore western Yemen, was constructed to evaluate the untapped oil and gas potential of the Upper Jurassic Madbi formation. 3D PSM techniques were used to analyze petroleum generation for conventional reservoirs and the petroleum saturations retained in the source rock for the unconventional system. Block 18 has several proven petroleum systems and producing oil and gas fields. The principal source rocks are within the Madbi Formation, which comprises two units, the Lam and the Meem members. Both contain transgressive organically rich "hot" shales with total organic carbon (TOC) of 8 to 10%; these are located stratigraphically at the base of each member. Additional organic-rich intervals within the Lam and Meem are less-effective source rocks, with lower TOC values.
The PSM consisted of 17 depositional events and 2 hiatuses. To accurately replicate geochemical and stratigraphic variations, the Lam and Meem members were further divided into sublayers. The model was calibrated to present-day porosity, permeability, and pressure data, and it incorporated vertical and lateral lithofacies and organic facies variations. Further calibrations used observed maturities (vitrinite reflectance and pyrolysis Tmax) and present-day temperatures and considered laterally variable heat flow from the Early Jurassic to the Late Miocene. Finally, petrophysical analyses from wells provided calculated hydrocarbon saturations, which were used to calibrate the saturation output from the model. The model satisfactorily reproduces the distribution of the main gas and oil fields and discoveries in the study area and is aligned with well test data.
Maturity results indicate that the upper Lam intervals currently sit within the main to early oil window but are immature at the edges of Block 18 (based on the Sweeney and Burnham Easy R0% kinetics). The lowest Lam unit enters the wet gas window in the center of the block. The underlying Meem member ranges from wet gas to early oil window maturity. Like the Lam, the Meem remains immature along the edges of Block 18. However, in the south of the block, the richest source rocks within the Meem are mainly in the oil window. The degree of transformation of the Meem and Lam varies throughout the members. The model predicts that, at present, the lowest part of the Meem, containing the greatest TOC, has 90% of its kerogen transformed into hydrocarbons.
The model confirms that the Madbi formation is a promising unconventional shale reservoir with a high quantity of hydrocarbons retained within it. Despite the higher quantity of hydrocarbons retained in the upper Meem, in terms of liquid and vapor hydrocarbons predicted in this model, the lower Lam is the most-prospective conventional tight sand reservoir, and the Meem has very small potential as tight sand reservoirs. This study provided a novel application of 3D PSM technology to assess new unconventional as well as conventional plays in this frontier area.
Hydrodynamics and geothermics are important tools for understanding the complex distribution of reservoir fluids in the Montney Formation in Alberta and British Columbia, Canada. The Montney comprises a conventional system in the east and an unconventional, Deep Basin-style hydrocarbon system in the west, where an underpressured, oil-dominated fairway just west and downdip of the conventional system grades further downdip into overpressured liquids and gas fairways.
The first part of this study addresses how these systems can be mapped from a pressure and temperature perspective. The Montney hydrodynamics system is explained using pressure versus elevation graphs. Key contours are taken from pressure-depth ratio maps to define the general boundaries between systems, noting that these boundaries change with depth. Geothermal gradient mapping is used to identify areas of prominent high or low geothermal gradients, which can have a significant effect on the positioning of gas liquids fairways. Key current day isotherms are also identified to represent the current phase windows by relating present-day formation temperatures to Tmax data.
To evaluate how pressure and temperature affect liquids production within the Montney, liquids production trends need to be considered. The second half of the paper discusses how mapping gas composition, particularly C2+ Wet Gas Index (WGI), may serve as a good proxy for liquids yields.
While the authors appreciate the complexities of phase behavior and the various factors influencing liquids production, the objective of this paper is to link trends that can be observed in liquids production to trends in pressure, temperature and gas composition. Ultimately, this paper examines ways in which hydrodynamics and geothermics can be used to help predict spatial variations in observed liquids production. By analyzing the co-relationships of the pressure, temperature and WGI data, the Montney segregates into two distinct domains which we term the Northern (British Columbia) Play and the Southern (Alberta) Play. This analysis can be tied in with other data sets for a better understanding of the reservoir such as: isotope geochemistry to gain insights into hydrocarbon migration; Special Core Analysis (SCAL) data to gain insights into fluid mobility; vapour-liquid equilibrium data to examine hydrocarbon fractionation during production; and completions data to provide a more complete picture of reservoir deliverability.
Tavassoli, Shayan (The University of Texas at Austin) | Krishnamurthy, Prasanna (The University of Texas at Austin) | Beckham, Emily (The University of Texas at Austin) | Meckel, Tip (The University of Texas at Austin) | Sepehrnoori, Kamy (The University of Texas at Austin)
Storage of large amounts of CO2 within deep underground aquifers has great potential for long-term mitigation of climate change. The U.S. Gulf Coast is an attractive target for CO2 storage because of the favorable formation properties for injection and containment of CO2. Deltaic formations are one of the primary targeted depositional environments in the Gulf Coast. This paper investigates CO2 storage in deltaic saline aquifers through a combination of geological modeling and flow simulation.
The geological model in our study is developed based on a laboratory-scale 3D flume experiment replicating the formation of a delta structure and populated with geologic properties according to Miocene Gulf of Mexico natural analogues. We used invasion percolation simulations to understand the gravity- driven flow and the relationship between architecture, stratigraphy, and fluid migration pathways. The results were used to develop an upscaled model for compositional simulation with the key features of the original geological model and to determine injection schemes that maximize the injection capacity and minimize the amount of mobile CO2 in the formation. In order to achieve this, we used compositional reservoir simulations to study the pressure-driven flow and phase behavior.
The results of invasion percolation simulations were used to identify the key stratigraphic units affecting CO2 migration. The realistic geometries and high resolution of the model facilitate the transfer of results from synthetic to subsurface data. The results allow for the analysis of deltaic depositional environments, important stratigraphic surfaces, and their impact on CO2 storage. The reservoir simulation model and phase behavior were validated against available field and lab data. The results of reservoir simulations were used to investigate the effects of main mechanisms, such as gas trapping and solubilization, on storage capacity. We compared our simulation results on the basis of invasion percolation (gravity driven) and reservoir simulation (pressure driven). The comparison is helpful to understand the strengths and weaknesses of each approach and determine best practices to evaluate CO2 migration within similar formations.
The unique and extremely well characterized deltaic model allows for unprecedented representation of the depositional aquifer architecture. This research combines geologic modeling, flow simulation, and application for CO2 storage. The integrated conclusions will constrain predictions of actual subsurface flow performance and CO2 storage capacity in deltaic systems, while identifying potential risks and primary stratigraphic migration pathways. This research gives insights on prediction of CO2 storage performance and characterization of prospective saline aquifers.
Diatto, Paolo (Eni S.p.A.) | Cerioli Regondi, Anita (Eni S.p.A.) | Doering, Sascha (Eni S.p.A.) | Italiano, Domenico (Eni S.p.A.) | Maffeis, Ivan (Eni S.p.A.) | Marchesini, Marco (Eni S.p.A.) | Martin, Marco (Eni S.p.A.)
With the aim of improving the understanding of production behaviour in a multi-discovery asset and the evaluation of near-field exploration opportunities, an integrated study has been carried out involving three different disciplines: Fluid Thermodynamics (PVT), Organic Geochemistry and Petroleum Systems Modelling (PSM). The synergistic workflow has been undertaken starting from an accurate quality check of the initial dataset related to fluid samples and lab tests. By merging PVT and geochemical data, it was possible to carry out a robust statistical survey and explore correlations across different parameters and features; in this way, strict connection among many physical parameters and some oil maturity and biodegradation indices were identified. In the following step, after geo-referencing the fluid samples in the framework of the Petroleum Systems Model and tracking the locations of the source rocks, a reliable interpretation of the oil expulsion and migration history became possible over the whole reservoir fluid system. Finally, taking into account the simulated fluid phase envelopes, further insights were drawn in terms of the fluid phase behavior in different areas, contributing to reduce uncertainty and exploration risk for future activity in nearby prospects.
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.
Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. In the context of fluids, migration is the movement of fluids, generally away from a source rock through permeable layers towards a trap or vent. In a seismic context, migration is a computation applied to seismic data that returns reflection events (signals) to their origin in the subsurface.
Although variables that affect kick-killing do not necessitate a change in the basic procedural structure, they may cause unexpected behaviors that can mislead an operator into choosing the wrong procedure. The one-circulation method is used to demonstrate the effect of these variables. The influx type entering the wellbore plays a key role in casing-pressure behavior. The influx can range from heavy oil to fresh water. The most common is gas or salt water; each has a pronounced casing pressure curve and different downhole effects.
In recent years, deformation of the reservoir host rocks has become a subject of great interest, prompted in part by the dramatic subsidence observed at Ekofisk platforms in the North Sea. One method of monitoring deformation is by passive seismic monitoring. It is called "passive" because the geopysicist does not activate a seismic source, but rather uses existing geophones to monitor ongoing changes in the rocks due to downhole conditions. Deformation is an important aspect of reservoir production, even without a significant compaction drive in many cases. Previous studies have been published in the scientific and earthquake literature relating earthquakes to oil/gas production and to injection practices.
The ability of seismic reflection technology to image subsurface targets is possible largely through the geometry of sources and receivers. A method similar to triangulation is used to place reflections in their correct locations with (more-or-less) correct amplitudes, which can then be interpreted. The amplitudes are indicative of relative changes in impedance, and the seismic volume can be processed to yield impedances between the reflecting boundaries. The geometry of sources and receivers in a typical reflection seismic survey yields a number of seismic traces with common midpoints or central bins for stacking. These traces were recorded at different offset distances, and the travel times for seismic waves traveling to and from a given reflecting horizon varies with that distance (Figure 1).
Despite the now-routine use of prestack depth migration (PSDM) for unconventionals, confusion abounds on the topic of how to best incorporate near-surface velocity estimates into the PSDM shallow-model-building process. The present work seeks to eliminate the confusion via a carefully-controlled synthetic experiment in which the (known) near-surface velocity distribution mimics typical Permian Basin shallow geology. In this experiment, various methods for near-surface model building are tested, ranging from simplistic to sophisticated, and PSDM results are compared against the ideal image. These tests clearly demonstrate that gather flattening improves dramatically with application of the more sophisticated shallow model building approaches. In the case of the most primitive approaches (e.g, migration-from-flat-datum or migration from topography where the shallow velocity cells are flooded with a spatially uniform “replacement” velocity), the migrated gathers exhibit significant residual moveout, and applying a tomographic velocity update to improve flattening leads to a significant error in event depth location (i.e, “depthing”). This depthing error suggests that downstream anisotropic parameter estimation will be compromised unless a more sophisticated shallow model building approach is employed. The concept of differential statics is introduced and is demonstrated to be a useful tool which can provide good gather flattening, accurate event depthing, and also improved lateral continuity of events in the common case where the near-surface velocity estimate from refraction statics analysis is not suitable for verbatim insertion into the shallow PSDM model. Key findings from the synthetic experiments are corroborated by analogous observations on real data, suggesting that the experiments are indeed capturing realistic effects.
It is well known that near-surface heterogeneity can cause significant traveltime distortion of reflected signals, and, furthermore, that such distortion poses a major challenge in land seismic imaging. Addressing this challenge is particularly important in unconventional plays, where accurate depthing of subtle features is crucial for applications such as landing and steering optimization. Recently, some notable advances have been made, including the use of novel refraction statics techniques (Diggins et al., 2016), application of full-waveform inversion (e.g., Roy et al., 2017), and incorporation of gravity/EM data (Colombo et al., 2012), all of which seek to better estimate the near-surface velocity field. At the same time as these advances are unfolding, prestack depth migration is beginning to see widespread use in many North American unconventional shale plays (Rauch-Davies et al., 2018).