ALBERT

Description: Summary The initial water saturation in a reservoir is important for both hydrocarbon volume estimation and distribution of multiphase flow properties such as relative permeability. Often, a practical reservoir engineering approach is to relate relative permeability to flow property regions by binning of the initial water saturation. The rationale behind this approach is that initial water saturation is related to both the pore-throat radius distribution and the wettability of the rock, both of which affect relative permeability. However, pore-throat radius and wettability are usually not explicitly included in geomodel property modeling. Therefore, the saturation height model should not only capture an average hydrocarbon pore volume but also reflect the underlying mechanisms from hydrocarbon migration history and its impact on initial water saturation distribution. This work introduces and describes a new term, excess water, for more precise classification of saturation height model scenarios in reservoirs in which multiple mechanisms have interacted and caused a complex water saturation distribution. An example of the presence of transition zones related to drained local perched aquifers (excess water) in oil-down-to (ODT) wells is shown using a limited data set from a North Sea reservoir. The physical basis for drainage and imbibition transition zones connected to both regional and perched aquifers is given. The distribution of initial water saturation in reservoirs containing excess water is demonstrated through numerical modeling of oil migration over millions of years. Highly permeable reservoirs are more likely to have locally trapped water because of lower capillary forces. A static situation occurs in areas where the capillary forces cannot maintain a high enough water saturation for further water drainage. On the other hand, both high- and low-permeability reservoirs may have significant excess water because of ongoing dynamic effects. In both cases, long distances for water to drain laterally to a regional aquifer enhance the possibility for a dynamic excess water situation.

Description: Summary The use of multilateral wells started in the mid-1990s in particular in Canada, and they have since been used in many countries. However, few papers on multilateral wells focus on their production performance. Thus, what can be expected from such wells in terms of production is not clear, and this paper will attempt to address that gap. Taking advantage of public data, the production performance of multilateral wells in various Western Canadian fields has been studied. In the cases reviewed in this paper, these wells always target a single formation; they have been used in a variety of fields and reservoirs, mostly for primary production but also for polymer flooding in some cases. Multiple examples will be provided, mostly in heavy oil reservoirs, and production performance will be compared with nearby horizontal wells whenever possible. From the more classical dual and trilateral, to more complex architectures with seven or eight laterals, and the more exotic with laterals drilled from laterals, the paper will present the architecture and performance of these complex wells and of some fields that have been developed almost exclusively with multilateral wells. Interestingly, multilateral wells have not been used much for secondary or tertiary recovery, probably because of the difficulty of controlling water production after breakthrough. However, field results suggest that this may not be such a difficult proposition. One of the most remarkable wells producing a 1,250-cp oil under polymer flood has achieved a cumulative production of more than 3 million bbl, which puts it among the top producers in Canada. Although multilateral wells have been in use for more than 25 years, very few papers have been devoted to the description of their production performance. This paper will bring some clarity to these aspects. It will also attempt to address when multilateral wells can be used and to compare their performance to that of horizontal wells in the same fields. It is hoped that this paper will encourage operators to reconsider the use of multilateral wells in their fields.

Description: Summary Forecasting production from hydrocarbon fields is challenging because of the large number of uncertain model parameters and the multitude of observed data that are measured. The large number of model parameters leads to uncertainty in the production forecast from hydrocarbon fields. Changing operating conditions [e.g., implementation of improved oil recovery or enhanced oil recovery (EOR)] results in model parameters becoming sensitive in the forecast that were not sensitive during the production history. Hence, simulation approaches need to be able to address uncertainty in model parameters as well as conditioning numerical models to a multitude of different observed data. Sampling from distributions of various geological and dynamic parameters allows for the generation of an ensemble of numerical models that could be falsified using principal-component analysis (PCA) for different observed data. If the numerical models are not falsified, machine-learning (ML) approaches can be used to generate a large set of parameter combinations that can be conditioned to the different observed data. The data conditioning is followed by a final step ensuring that parameter interactions are covered. The methodology was applied to a sandstone oil reservoir with more than 70 years of production history containing dozens of wells. The resulting ensemble of numerical models is conditioned to all observed data. Furthermore, the resulting posterior-model parameter distributions are only modified from the prior-model parameter distributions if the observed data are informative for the model parameters. Hence, changes in operating conditions can be forecast under uncertainty, which is essential if nonsensitive parameters in the history are sensitive in the forecast.

Description: Summary Volumetric maximum curvature attribute computed from 3D ocean bottom cable (OBC) seismic data, production logging tool (PLT), inorganic chemical tracer data, and fractures observed from core and full-bore formation microimager (FMI) logs were integrated to characterize fractured carbonate reservoirs of an offshore oil field in Abu Dhabi, United Arab Emirates (UAE). The extracted maximum curvature anomalies are predominantly orientated in NNE-SSW and NE-SW, a trend perpendicular to the dominant fault direction in the oil field and similar to the dominant strike directions of fractures measured from core data and FMI logs. Because the fracture strike directions of well data mimic the strike directions of curvature anomalies at corresponding reservoir levels, we interpreted the maximum curvature anomalies to represent dilatational fractured zones or fracture corridors. Integration of dynamic data, such as PLT and chemical tracers, and maximum curvature anomalies demonstrate that the inferred fracture zones can determine water breakthroughs as well as inter- and intrareservoir communications. As a result, this study highlights possible fracture zones and their internal architecture, as well as their potential flow capabilities. These results play a key role in reservoir management and monitoring of water movement through structural pathways.

Description: Summary Geoscientific and engineering experiments in petrophysics, rock physics, and rock mechanics depend on multiple, costly, and sometimes rare samples used to characterize the properties of natural rocks. Testing these samples helps in modeling various hydrocarbon recovery and stimulation scenarios, as well as understanding the fluid-rock interactions in the subsurface under various pressure and temperature conditions. Over the last decade, 3D printing has matured to become a more commonly available tool to enable repeatable experiments with controllable materials and pore system geometries to investigate petrophysical, geomechanical, and geophysical properties of porous rocks. This review introduces the development, characteristics, and capabilities of 3D printing technology that are specifically used in research. Applications in the realm of petrophysics highlight the issues of replicating the pore network geometry and subsurface physics, aiming at understanding fluid flow in porous media problems. Using 3D-printed models in rock mechanics experiments focuses on generating comparable geomechanical properties and reproducing fractures, joint surfaces, and other rock structures, whereas in rock physics, geophysical forward modeling is highlighted to take advantage of 3D printing technology. By summarizing the recent advances in 3D printing as applied to petrophysics, rock physics, and rock mechanics, this review paper presents the current state of the art and the challenges in scale, cost, time, and materials, as well as the directions for advancing this frontier discipline to answer various fundamental questions regarding porous media research using 3D printing technology.

Description: Summary Carbon dioxide (CO2) enhanced oil recovery (EOR) has long been practiced in the US as an efficient mean for enhancing oil production. Many of the US CO2-EOR developments have been designed horizontally. This is because of a viscous-dominated CO2 flow regime that is prevalent in these developments driven by thin and low-permeability reservoirs. Reservoirs and fluid properties are different in the North Sea. Pays are usually thicker with better petrophysical properties. Lighter oils can also be found in North Sea reservoirs. This suggests that a dissimilar flow regime might prevail CO2 displacements in the North Sea developments, which could favor a dissimilar CO2-EOR process design. This study thus compares CO2 flow regimes between several North Sea and US reservoirs. We use scaling analysis to characterize and compare CO2 flow regimes between these two classes of reservoirs. Scaling analysis characterizes CO2 displacement in each reservoir system using three dimensionless numbers: gravity, effective aspect ratio, and mobility ratio. Displacement experiments conducted in stochastically generated permeability fields, under exactly matched magnitudes of the derived dimensionless numbers, reveal the prevailing CO2 flow regime in each reservoir system. Results of scaling analysis indicate that CO2 flooding in the North Sea reservoirs can be generally characterized with a larger gravity number, smaller effective aspect ratio, and smaller mobility ratio than the average US CO2 flooded reservoirs. Flow regime analysis indicates that unlike the majority of the US CO2 flooded reservoirs, CO2 flow regimes tend to be more gravity-dominated in the North Sea class of reservoirs. CO2 flow regimes in the North Sea systems are expected to suffer from a higher degree of instability because of thicker North Sea pays, which limit effective crossflow. Understanding the differences and characteristics of CO2 flow regimes in the North Sea prospects can help operators design their CO2 flooding more efficiently, which could increase the recovery factor (RF) as well as address CO2 storage requirements, a necessary consideration for CO2-EOR deployment in the North Sea.

Description: Summary The knowledge of the effects of instability and heterogeneity on displacements, primarily enhanced oil recovery, and carbon dioxide storage are well known, although they remain difficult to predict. The usual recourse to modeling these effects is through numerical simulation. Simulation remains the gold standard for prediction; however, its results lack generality, being case-specific. There are also several analytic models for displacements that are usually more informative than simulation results. However, these methods apply to steady-state, incompressible flow. Carbon dioxide injection for storage uses compressible fluids and, in the absence of producers, will not approach steady-state flow (Wu et al. 2017). Consequently, it is unlikely that storage will be in reservoirs of open boundaries (steady-state flow). Flow of compressible fluid necessitates the use of closed or partially sealed boundaries, a factor that is consistent with compressible flow. This work deals with the conditions that cause the onset of incipient viscous fingering or Saffman-Taylor (ST) instability. The actual growth and propagation of fingers, a subject of much recent literature, is not discussed here. The original ST formalism of M 〉 1 for gravity-free flow is highly restrictive: it is for linear flow of nonmixing incompressible fluids in steady-state flow. In this work, we relax the incompressible flow restriction and thereby broaden the ST criterion to media that have sealing and/or partially sealing outer boundaries. We use the nonlinear partial differential equation for linear flow and developed analytic solutions for a tracer flow analog and also for a two-fluid compressible flow. The analysis is restricted to stabilized flow and to constant compressibility fluids, but we are not restricted to small compressibility fluids. There is no transition (mixing) zone between displacing and displaced fluids; the displacement is piston-like. The absence of a transition zone means that the results apply to both miscible and immiscible displacements, absent dispersion, or local capillary pressure. The assumption of a sharp interface is to focus on the combined effect of mobility ratio and compressibility. We use the product of the fluid compressibility and pressure drop (cfΔP) to differentiate the compressibility groups (Dake 1978; Dranchuk and Quon 1967), where ΔP is defined as the pressure drop within the specific fluid region. The results will be based on proposed analytical solutions compared to numerical simulation. The proposed formulation is less restrictive than the original ST formalism of M 〉 1 and allows evaluation of viscous fingering initiation or ST stability criterion in the presence of different boundary conditions (open vs. closed boundaries) with compressible fluids under the stated assumptions, which is the scope of this work. The key contribution here is the effect of external boundaries, which consequently makes necessary the use of compressible fluids. Absent compressibility, the necessary condition for the growth of a viscous finger is simply the mobility ratio, M 〉 1. It is the objective of this work to study how the ST criterion is affected by the presence of sealing and partially sealing outer boundaries with the consequent inclusion of compressible flows as in carbon dioxide storage and enhanced oil recovery by gas injection. The results show that adding compressibility always makes displacements more unstable for stabilized background flow, even for a favorable mobility ratio (M  5×10−3 1/psi). For a sealed external boundary (no production or leakage), displacements will become more stable as a front approaches an external boundary for all mobility ratios (M) investigated.

Description: Summary An approach to the analysis of production data from waterflooded oil fields is proposed in this paper. The method builds on the established techniques of rate-transient analysis (RTA) and extends the analysis period to include the transient- and steady-state effects caused by a water-injection well. This includes the initial rate transient during primary production, the depletion period of boundary-dominated flow (BDF), a transient period after injection starts and diffuses across the reservoir, and the steady-state production that follows. RTA will be applied to immiscible displacement using a graph that can be used to ascertain reservoir properties and evaluate performance aspects of the waterflood. The developed solutions can also be used for accurate and rapid forecasting of all production transience and boundary-dominated behavior at all stages of field life. Rigorous solutions are derived for the transient unit mobility displacement of a reservoir fluid, and for both constant-rate-injection and constant-pressure-injection after a period of reservoir depletion. A simple treatment of two-phase flow is given to extend this to the water/oil-displacement problem. The solutions are analytical and are validated using reservoir simulation and applied to field cases. Individual wells or total fields can be studied with this technique; several examples of both will be given. Practical cases are given for use of the new theory. The equations can be applied to production-data interpretation, production forecasting, injection-water allocation, and for the diagnosis of waterflood-performanceproblems. Correction Note: The y-axis of Fig. 8d was corrected to "Dimensionless Decline Rate Integral, qDdi". No other content was changed.

Description: Summary The extensive depletion of the development target triggers the demand for infill drilling in the upside target of multilayer unconventional reservoirs. However, such an infill scheme in the field practice still heavily relies on empirical knowledge or pressure responses, and the geomechanics consequences have not been fully understood. Backed by the data set from the Permian Basin, in this work we present a novel integrated reservoir-geomechanics-fracture model to simulate the spatiotemporal stress evolution and locate the optimal development strategy in the upside target of the Bone Spring Formation. An embedded discrete fracture model (EDFM) is deployed in our fluid-flow simulation to characterize complex fractures, and the stress-dependent matrix permeability and fracture conductivity are included through the compaction/dilation option. After calibrating reservoir and fracture properties by history matching of an actual well in the development target (i.e., third Bone Spring), we run the finite element method (FEM)-based geomechanics simulation to model the 3D stress state evolution. Then a displacement discontinuity method (DDM) hydraulic fracture model is applied to simulate the multicluster fracture propagation under an updated heterogeneous stress field in the upside target (i.e., second Bone Spring). Numerical results indicate that stress field redistribution associated with parent-well production indeed vertically propagates to the upside target. The extent of stress reorientation at the infill location mainly depends on the parent-child horizontal offset, whereas the stress depletion is under the combined impact of horizontal offset, vertical offset, and infill time. A smaller parent-child horizontal offset aggravates the overlap of the stimulated reservoir volume (SRV), resulting in more substantial interwell interference and less desirable oil and gas production. The same trend is observed by varying the parent-child vertical offset. Moreover, the efficacy of an infill operation at an earlier time is less affected by parent-well depletion because of the less-disturbed stress state. The candidate infill-well locations at various infill timings are suggested based on the parent-well and child-well production cosimulation. Being able to incorporate both pressure and stress responses, the reservoir-geomechanics-fracture model delivers a more comprehensive understanding and a more integral solution of infill-well design in multilayer unconventional reservoirs. The conclusions provide practical guidelines for the subsequent development in the Permian Basin.

Description: Summary Achieving effective results using conventional thermal recovery technology is challenging in the deep undisturbed reservoir with extra-heavy oil in the LKQ oil field. Therefore, in this study, a novel approach based on in-situ combustion huff-and-puff technology is proposed. Through physical and numerical simulations of the reservoir, the oil recovery mechanism and key injection and production parameters of early-stage ultraheavy oil were investigated, and a series of key engineering supporting technologies were developed that were confirmed to be feasible via a pilot test. The results revealed that the ultraheavy oil in the LKQ oil field could achieve oxidation combustion under a high ignition temperature of greater than 450°C, where in-situ cracking and upgrading could occur, leading to greatly decreased viscosity of ultraheavy oil and significantly improved mobility. Moreover, it could achieve higher extra-heavy-oil production combined with the energy supplement of flue gas injection. The reasonable cycles of in-situ combustion huff and puff were five cycles, with the first cycle of gas injection of 300 000 m3 and the gas injection volume per cycle increasing in turn. It was predicted that the incremental oil production of a single well would be 500 t in one cycle. In addition, the supporting technologies were developed, such as a coiled-tubing electric ignition system, an integrated temperature and pressure monitoring system in coiled tubing, anticorrosion cementing and completion technology with high-temperature and high-pressure thermal recovery, and anticorrosion injection-production integrated lifting technology. The proposed method was applied to a pilot test in the YS3 well in the LKQ oil field. The high-pressure ignition was achieved in the 2200-m-deep well using the coiled-tubing electric igniter. The maximum temperature tolerance of the integrated monitoring system in coiled tubing reached up to 1200°C, which provided the functions of distributed temperature and multipoint pressure measurement in the entire wellbore. The combination of 13Cr-P110 casing and titanium alloy tubing effectively reduced the high-temperature and high-pressure oxygen corrosion of the wellbore. The successful field test of the comprehensive supporting engineering technologies presents a new approach for effective production in deep extra-heavy-oil reservoirs.