Proactive Decision-Making Through Real-Time Geomechanical Support Leveraging Drilling a Long Horizontal Section Through a Tight Unconventional Reservoir of an Oil and Gas Field: Middle East

This paper discusses the Integrated Role of Geomechanics and Drilling Fluids Design for drilling a well oriented towards the minimum horizontal stress direction in a depleted, yet highly stressed and complex clastic reservoir. There are multiple challenges related to such a well that need to be addressed during the planning phase. In this case, the well needs to be drilled towards the minimum horizontal stress direction (Shmin) to benefit multi-stage hydraulic fracturing. At the same time, the most prominent challenge is that this well orientation is more prone to wellbore failure and requires a maximum mud weight, due to the present strike slip stress environment. Well planning challenges in such an environment include (a) the determination of formation characteristics and rock properties, (b) the anticipation of higher formation collapse pressure during the course of drilling the lateral section within the reservoir, (c) the determination of the upper bound mud weight to prevent lost circulation due to a low fracture gradient against depleted sections, or due to the presence of pre-existing natural fractures, d) mitigating the higher risk of differential sticking against depleted porous layers, and determining appropriate bridging in the drilling fluids, (e) recognizing the prolonged exposure time of the formation due to the length of the lateral and the lower rate of penetration against the tight highly dense formations. For successful drilling, and to mitigate the above risks, the first step is to prepare a predrill GeoMechanical model along with adequate fluid design and drillers action plans to be considered during drilling. Offset well petrophysical logs and core data are considered for the preparation of the predrill GeoMechanical model, along with the drilling experiences in the offset locations. Based on the above, a predrill GeoMechanical model is prepared, a risk matrix is being established, and a representative mud weight window is recommended (Wellbore Stability Analysis). In most cases, the offset well locations considered are vertical- or inclined-, or lateral wells of different trajectory azimuth than the target well location and the predrill GeoMechanical model can incorporate such variations easily; however, any Geology uncertainty, leading to a different rock property- and stress set-up (or even different pore pressure than expected), at the actual well location will be part of the uncertainty of the predrill GeoMechanical model and Wellbore Stability Analysis. This is where the real time monitoring is playing out its full potential: giving an updated model and wellbore stability analysis during drilling. While drilling the lateral section, the wellbore condition is being monitored using LWD (logging while drilling) tools, e.g. Gamma Ray, Density, Neutron, Acoustic Caliper, Azimuthal density image and ECD (equivalent circulating density). While gamma ray helps in determining the lithology, density logs help to understand the formation hardness, and they can be used to generate a calibrated pseudo acoustic log. Based on this pseudo acoustic log, the rock strength and other rock mechanical properties of the pre- GeoMechanical model can be updated as soon as they become available. This gives insight into the model differences and helps to understand model variations and adjust Wellbore Stability recommendations accordingly. While the neutron log helps to determine the zones of high porosity, and thus potential risk zones for differential sticking, the azimuthal density image clearly indicates the breakout zones caused by the shear failure of the wellbore. The presence of wellbore failure (breakout) is further confirmed by acoustic caliper data, and accordingly wellbore stability related recommendations are communicated to the operator, for an increase in the specific gravity of the mud, and thus, to balance the wellbore. From a mud rheology perspective, high performance OBM (oil-based mud) parameters are maintained consistent with the formation properties, to minimize fluid loss, optimize wellbore strengthening characteristics and minimize at the same time solids concentrations in order to avoid excessive ECD (equivalent circulating density) which may open pre-existing natural fractures resulting in downhole losses and in consequence might lead to differential sticking. In the case study presented herein, the proactive implementation of GeoMechanics and its Wellbore Stability application as well as the integration of drilling fluids services, resulted in the smooth and successful drilling of the lateral section, and also in the delivery of an in gauge hole necessary for multi-stage fracturing (MSF) completion optimization.

Real-Time Wellbore Stability and Hole Quality Evaluation Using LWD Azimuthal Photoelectric Measurements

Horizontal and high-inclination deep wells are routinely drilled to enhance hydrocarbon recovery. To sustain production rates, these wells are generally designed to be drilled in the direction of minimum horizontal stress in strike slip stress regime to facilitate transverse fracture growth during fracturing operations. These wells can also cause wellbore instability challenges due to high stress concentration due to compressional or strike-slip stress regimes. Hence, apart from pre-drill wellbore stability analysis for an optimum mud weight design, it is important to continuously monitor wellbore instability indicators during drilling. With the advancements of logging-while-drilling (LWD) techniques, it is now possible to better assess wellbore stability during drilling and, if required, to take timely decisions and adjust mud weight to help mitigate drilling problems. The workflow for safely drilling deep horizontal wells starts with analyzing the subsurface stress regime using data from offset wells. Through a series of steps, data is integrated to develop a geomechanics model to select an optimum drilling-fluid density to maintain wellbore stability while minimizing the risks of differential sticking and mud losses. Due to potential lateral subsurface heterogeneity, continuous monitoring of drilling events and LWD measurements is required, to update and calibrate the pre-well model. LWD measurements have long been used primarily for petrophysical analysis and well placement in real time. The use of azimuthal measurements for real-time wellbore stability evaluation applications is a more recent innovation. Shallow formation density readings using azimuthal LWD measurements provide a 360° coverage of wellbore geometry, which can be effectively used to identify magnitude and orientation of borehole breakout at the wellbore wall. Conventional LWD tools also provide auxiliary azimuthal measurements, such as photoelectric (Pe) measurement, derived from the near detector of typical LWD density sensors. The Pe measurement, with a very shallow depth of investigation (DOI), is more sensitive to small changes in borehole shape compared with other measurements from the same sensor, particularly where a high contrast exists between drilling mud and formation Pe values. Having azimuthal measurements of both Pe and formation density while drilling facilitates better control on assess wellbore stability assessment in real time and make decisions on changes in mud density or drilling parameters to keep wellbore stable and avoid drilling problems. Time dependency of borehole breakout can also be evaluated using time-lapse data to enhance analysis and reduce uncertainty. Analyzing LWD density and Pe azimuthal data in real time has guided real-time decisions to optimize drilling fluid density while drilling. The fluid density indicated by the initial geo-mechanical analysis has been significantly adjusted, enabling safe drilling of deep horizontal wells by minimizing wellbore breakouts. Breakouts identified by LWD density and photoelectric measurements has been further verified using wireline six-arm caliper logs after drilling. Contrary to routinely used density image, this paper presents use of Pe image for evaluating wellbore stability and quality in real time, thereby improving drilling safety and completion of deep horizontal wells drilled in the minimum horizontal stress direction.

Integrated Mechanical Earth Modeling for Predicting Sand Production: A Case Study

Summary Sand production is a serious problem in oil and gas wells, and one of the main concerns of production engineers. This problem can damage downhole equipment and surface production facilities. This study presents a sand production case and quantifies sanding risks for an oil field in Iraq. The study applies an integrated workflow of constructing 1D Mechanical Earth Modeling (MEM) and predicting the sand production with multiple criteria such as shear failure during drilling, B index, and critical bottomhole pressure (CBHP) or critical drawdown pressure (CDDP). Wireline log data were used to estimate the mechanical properties of the formations in the field. The predicted sand production propensity was validated based on the sand production history in the field. The interpretation results of some wells anticipated in this study showed that when a shear failure occurs during drilling, the B index is around 2 × 104 MPa or less and the CBHP is equal to the formation pore pressure. For this case, sand control shall be carried out in the initial stage of production. On the other hand, when the shear failure does not exist, the B index is always greater than 2 × 104 MPa, and the CBHP is mostly less than the formation pore pressure. In this case, implementing sand control methods could be postponed as the reservoir pressure undergoes depletion. However, for the anticipated field, sand control is recommended to be carried out in the initial stage of well production even when the CBHP is less than the formation pore pressure since sanding will be inevitable when the reservoir pressure depletes to values close to the initial reservoir pressure. The tentative evaluation of the stress regime showed that a normal fault could be the stress regime for the formations. For a normal fault stress regime, the study explained that when the reservoir permeability is isotropic, an openhole vertical wellbore has less propensity for sand production than a horizontal wellbore. Moreover, when the wellbore azimuth is in the direction of the minimum horizontal stress, the CBHP will be lower than in any other azimuth, and sanding will take place at higher wellbore inclination angles. For the anticipated field, because of the casedhole well completion and the anisotropic reservoir permeability, a horizontal well drilled in the direction of minimum horizontal stress with oriented perforation in the direction of maximum horizontal stress is an alternative method for controlling sand production.

4D Finite element modeling of stress distribution in depleted reservoir of south Iraq oilfield

The harvest of hydrocarbon from the depleted reservoir is crucial during field development. Therefore, drilling operations in the depleted reservoir faced several problems like partial and total lost circulation. Continuing production without an active water drive or water injection to support reservoir pressure will decrease the pore and fracture pressure. Moreover, this depletion will affect the distribution of stress and change the mud weight window. This study focused on vertical stress, maximum and minimum horizontal stress redistributions in the depleted reservoirs due to decreases in pore pressure and, consequently, the effect on the mud weight window. 1D and 4D robust geomechanical models are built based on all available data in a mature oil field. The 1D model was used to estimate all mechanical rock properties, stress, and pore pressure. The minimum and maximum horizontal stress were determined using the poroelastic horizontal strain model. Furthermore, the mechanical properties were calibrated using drained triaxial and uniaxial compression tests. The pore pressure was tested using modular dynamic tester log MDT. The Mohr–Coulomb model was applied in the 4D model to calculate the stress distribution in the depleted reservoir. According to study wells, the target area has been classified into four main groups in Mishrif reservoir based on depletion: highly, moderately, slightly, and no depleted region. Also, the results showed that the units had been classified into three main categories based on depletion state (from above to low depleted): L1.1, L1.2, and M1. The mean average reduction in minimum horizontal stress magnitude was 322 psi for L1.1, 183.86 psi for L1.2, and 115.56 psi for M1. Thus, the lower limit of fracture pressure dropped to a high value in L1.1, which is considered a weak point. As a result of changing horizontal stress, the mud weight window became narrow.

Multiscale Assessment of Caprock Integrity for Geologic Carbon Storage in the Pennsylvanian Farnsworth Unit, Texas, USA

Leakage pathways through caprock lithologies for underground storage of CO2 and/or enhanced oil recovery (EOR) include intrusion into nano-pore mudstones, flow within fractures and faults, and larger-scale sedimentary heterogeneity (e.g., stacked channel deposits). To assess multiscale sealing integrity of the caprock system that overlies the Morrow B sandstone reservoir, Farnsworth Unit (FWU), Texas, USA, we combine pore-to-core observations, laboratory testing, well logging results, and noble gas analysis. A cluster analysis combining gamma ray, compressional slowness, and other logs was combined with caliper responses and triaxial rock mechanics testing to define eleven lithologic classes across the upper Morrow shale and Thirteen Finger limestone caprock units, with estimations of dynamic elastic moduli and fracture breakdown pressures (minimum horizontal stress gradients) for each class. Mercury porosimetry determinations of CO2 column heights in sealing formations yield values exceeding reservoir height. Noble gas profiles provide a “geologic time-integrated” assessment of fluid flow across the reservoir-caprock system, with Morrow B reservoir measurements consistent with decades-long EOR water-flooding, and upper Morrow shale and lower Thirteen Finger limestone values being consistent with long-term geohydrologic isolation. Together, these data suggest an excellent sealing capacity for the FWU and provide limits for injection pressure increases accompanying carbon storage activities.

Diagnostic Plots of Pressure Transient, Rate-Transient, and Diagnostic Fracture-Injection/Falloff Tests

All transient test interpretation methods rely on or utilize diagnostic plots for the identification of wellbore or fracture storage distortion, flow regimes, and other parameters (e.g., minimum horizontal stress). Although all "test" interpretations of interest are transient test data (i.e., those involving an "event"), the associated diagnostic plots are not interchangeable between such tests. The objective of this work is to clearly define the appropriate diagnostic plot(s) for each type of transient test. The work applies the appropriate transient test theory to demonstrate the applicability of each diagnostic plot along with clearly defining the characteristic features that make a given plot "diagnostic." For pressure transient testing, the material is largely a review, but for rate transient tests and diagnostic fracture-injection/falloff tests, new ideas are introduced and documented to justify appropriate diagnostic plots. Data examples are provided for illustration and application. In general, pressure transient test diagnostic plots are not misused, but the same cannot be said for diagnostic fracture-injection/falloff tests (or DFITs) where it is common to ascribe flow regimes and/or draw other erroneous conclusions based on observations from an inappropriately constructed or interpretated diagnostic plot. The examples provided illustrate both the correct diagnostic plot and interpretations, but also illustrate how data can be easily misinterpreted in common practice.

In Situ Stress and Stress Regime in the Onshore Part of the Northeast Java Basin

In situ stress is importance in the petroleum industry because it will significantly enhance our understanding of present-day deformation in a sedimentary basin. The Northeast Java Basin is an example of a tectonically active basin in Indonesia. However, the in situ stress in this basin is still little known. This study attempts to analyze the regional in situ stress (i.e., vertical stress, minimum and maximum horizontal stresses) magnitude and orientation, and stress regime in the onshore part of the Northeast Java Basin based on twelve wells data, consist of density log, direct/indirect pressure test, and leak-off test (LOT) data. The magnitude of vertical ( and minimum horizontal ( stresses were determined using density log and LOT data, respectively. Meanwhile, the orientation of maximum horizontal stress ( was determined using image log data, while its magnitude was determined based on pore pressure, mudweight, and the vertical and minimum horizontal stresses. The stress regime was simply analyzed based on the magnitude of in situ stress using Anderson’s faulting theory. The results show that the vertical stress ( ) in wells that experienced less erosion can be determined using the following equation: , where is in psi, and z is in ft. However, wells that experienced severe erosion have vertical stress gradients higher than one psi/ft ( . The minimum horizontal stress ( ) in the hydrostatic zone can be estimated as, while in the overpressured zone, . The maximum horizontal stress ( ) in the shallow and deep hydrostatic zones can be estimated using equations: and , respectively. While in the overpressured zone, . The orientation of is ~NE-SW, with a strike-slip faulting stress regime.

Rupture Directivity in 3D Inferred From Acoustic Emissions Events in a Mine-Scale Hydraulic Fracturing Experiment

Rupture directivity, implying a predominant earthquake rupture propagation direction, is typically inferred upon the identification of 2D azimuthal patterns of seismic observations for weak to large earthquakes using surface-monitoring networks. However, the recent increase of 3D monitoring networks deployed in the shallow subsurface and underground laboratories toward the monitoring of microseismicity allows to extend the directivity analysis to 3D modeling, beyond the usual range of magnitudes. The high-quality full waveforms recorded for the largest, decimeter-scale acoustic emission (AE) events during a meter-scale hydraulic fracturing experiment in granites at ∼410 m depth allow us to resolve the apparent durations observed at each AE sensor to analyze 3D-directivity effects. Unilateral and (asymmetric) bilateral ruptures are then characterized by the introduction of a parameter κ, representing the angle between the directivity vector and the station vector. While the cloud of AE activity indicates the planes of the hydrofractures, the resolved directivity vectors show off-plane orientations, indicating that rupture planes of microfractures on a scale of centimeters have different geometries. Our results reveal a general alignment of the rupture directivity with the orientation of the minimum horizontal stress, implying that not only the slip direction but also the fracture growth produced by the fluid injections is controlled by the local stress conditions.

In Situ Stress and Stress Regime in the Onshore Part of the Northeast Java Basin

In situ stress is importance in the petroleum industry because it will significantly enhance our understanding of present-day deformation in a sedimentary basin. The Northeast Java Basin is an example of a tectonically active basin in Indonesia. However, the in situ stress in this basin is still little known. This study attempts to analyze the regional in situ stress (i.e., vertical stress, minimum and maximum horizontal stresses) magnitude and orientation, and stress regime in the onshore part of the Northeast Java Basin based on twelve wells data, consist of density log, direct/indirect pressure test, and leak-off test (LOT) data. The magnitude of vertical ( and minimum horizontal ( stresses were determined using density log and LOT data, respectively. Meanwhile, the orientation of maximum horizontal stress ( was determined using image log data, while its magnitude was determined based on pore pressure, mudweight, and the vertical and minimum horizontal stresses. The stress regime was simply analyzed based on the magnitude of in situ stress using Anderson’s faulting theory. The results show that the vertical stress ( ) in wells that experienced less erosion can be determined using the following equation: , where is in psi, and z is in ft. However, wells that experienced severe erosion have vertical stress gradients higher than one psi/ft ( . The minimum horizontal stress ( ) in the hydrostatic zone can be estimated as, while in the overpressured zone, . The maximum horizontal stress ( ) in the shallow and deep hydrostatic zones can be estimated using equations: and , respectively. While in the overpressured zone, . The orientation of is ~NE-SW, with a strike-slip faulting stress regime.