A Single Equation to Depict Bottomhole Pressure Behavior for a Uniform Flux Hydraulic Fractured Well

Pressure transient analysis for a vertically hydraulically fractured well is evaluated using two different equations, which cater for linear flow at the early stage and radial flow in the later stage. However, there are three different stages that take place for an analysis of pressure transient, namely linear, transition and pseudo-radial flow. The transition flow regime is usually studied by numerical, inclusive methods or approximated analytically, for which no specific equation has been built, using the linear and radial equations. Neither of the approaches are fully analytical. The numerical, inclusive approach results in separate calculations for the different flow regimes because the equation cannot cater for all of the regimes, while the analytical approach results in a difficult inversion process to compute well test-derived properties such as permeability. There are two types of flow patterns in the fracture, which are uniform and non-uniform, called infinite conductivity in a high conductivity fracture. The study was conducted by utilizing an analogous study of linear flow equations. Instead of using the conventional error function, the exponential integral with an infinite number of wells was used. The results obtained from the developed analytical solution matched the numerical results, which proved that the equation was representative of the case. In conclusion, the generated analytical equation can be directly used as a substitute for current methods of analyzing uniform flow in a hydraulically fractured well.

Evaluating the Performance of Horizontal Multi-Frac Wells in a Depleted Gas Condensate Reservoir in Sultanate of Oman

The Saih Rawl gas condensate field has been producing for 20 years from multiple fractured vertical wells covering a very thick gross interval with varying reservoir permeability. After many years of production, the remaining reserves are mainly in the lowest permeability upper units. A pilot program using horizontal multi-frac wells was started in 2015, and five wells were drilled, stimulated and tested over a four-year period. The number of stages per horizontal well ranged from 6 to 14, but in all cases production was much less than expected based on the number of stages and the production from offset vertical wells producing from the same reservoir units with a single fracture. The scope of this paper is to describe the work that was performed to understand the reason for the lower than expected performance of the horizontal wells, how to improve the performance, and the implementation of those ideas in two additional horizontal wells completed in 2020. The study workflow was to perform an integrated analysis of fracturing, production and well test data, in order to history match all available data with a consistent reservoir description (permeability and fracture properties). Fracturing data included diagnostic injections (breakdown, step-rate test and minifrac) and main fracture treatments, where net pressure matching was performed. After closure analysis (ACA) was not possible in most cases due to low reservoir pressure and absence of downhole gauges. Post-fracture well test and production matching was performed using 3D reservoir simulation models including local grid refinement to capture fracture dimensions and conductivity. Based on simulation results, the effective propped fracture half-length seen in the post-frac production was extremely small, on the order of tens of meters, in some of the wells. In other wells, the effective fracture half-length was consistent with the created propped half-length, but the fracture conductivity was extremely small (finite conductivity fracture). The problems with the propped fractures appear to be related to a combination of poor proppant pack cleanup, low proppant concentration and small proppant diameter, compounded by low reservoir pressure which has a negative impact on proppant regained permeability after fracturing with crosslinked gel. Key conclusions from this study are that 1) using the same fracture design in a horizontal well with transverse fractures will not give the same result as in a vertical well in the same reservoir, 2) the effect of depletion on proppant pack cleanup in high temperature tight gas reservoirs appears to be very strong, requiring an adjustment in fracture design and proppant selection to achieve reasonable fracture conductivity, and 3) achieving sufficient effective propped length and height is key to economic production.

No Need for Radio Active Tracer When Mother Nature Helps Detecting the Effective Hydraulic Fracture Height

The tight gas reservoirs of Haima Supergroup provide the majority of gas production in the Sultanate of Oman. The paper discusses a possibility of using the anomalies from natural radioactivity to evaluate the fracture height for complex tight gas in mature fields of Oman. The standard industry practice is adding radioactive isotopes to the proppant. Spectral Gamma Ray log is used to determine near wellbore traced proppant placement. Spectral Noise log in combination with Production logs helps to identify the active fractures contributing to production. These methods complement each other, but they are obviously associated with costs. Hence, majority of wells are fracced without tracers or any other fracture height diagnostics. However, in several brown fields, an alternative approach to identify fracture height has been developed which provides fit-for-purpose results. It is based on the analysis of naturally occurring radioactive minerals (NORM) precipitation. The anomalies were observed in the many gas reservoirs even in cases when tracers were not used. At certain conditions, these anomalies can be used to characterize fracture propagation and optimize future wells hydraulic Fracture design. A high number of PLTs and well test information were analyzed. Since tight formations normally don't produce without fracturing, radioactive anomalies flag the contributing intervals and hence fracture propagation. The main element of analysis procedure is related to that fact that if no tracers applied, the discrepancy between normalized Open Hole Gamma Ray and Gamma Ray taken during PLT after 6-12 months of production can be used instead to establish fracture height. This method cannot be applied for immediate interpretation of fracture propagation because time is required to precipitate NORM and using the anomalies concept. The advantage of this method is that it can be used in some fields to estimate the frac effectiveness of wells without artificial tracers. It is normally assumed that the Natural radioactivity anomalies appear mainly due to co-production of the formation water. However, in the fields of interest the anomalies appear in wells producing only gas and condensate. This observation provides an opportunity for active fracture height determination at minimum cost.

Applying Slim Dipole Sonic Tool to Evaluate Hydraulic Fracture Geometry and Orientation Using the Differential Sonic Method

Successful hydraulic fracturing is critical for hydrocarbon recovery from tight reservoirs. Fracture geometry is one essential quality indicator of the created fracture. The geometry provides information about the size of the created fracture and containment and verifies the pre-job modeling. Different techniques are applied to determine fracture geometry, and each has its own advantages and limitations. Due to its simplicity, the radioactive tracer log is commonly used to determine fracture placement and fracture height. Its main drawbacks include shallow depth of investigation, time dependency, and the requirement for multiple interventions for multistage fracturing operations. The crosswell microseismic technique probes a larger volume and it is potentially capable of providing fracture height, length, and orientation. Operational complexity and long processing turnaround time are the main challenges of this technique. Time-lapse shear slowness anisotropy analysis is an effective method to determine hydraulic facture height and orientation. In this technique, the shear slowness anisotropy is recorded before and after the fracture is created. The observed shear anisotropy difference indicates the intervals where the fractures were created, allowing these intervals lengths to be measured. Combining this analysis with gyroscopic data allows determining the fracture orientations. Compared to a tracer log, the differential casedhole sonic anisotropy (DCHSA) has a deeper depth of investigation, and it is time independent. Thus, the repeated log can be acquired at the end of the multistage fracturing operations. Compared to the microseismic technique, this new technique provides more precise fracture height and orientation. The new generation slim dipole sonic technology of 2.125-in. diameter extends the applicability of the DCHSA technique to smaller casing sizes. The shear differential method was applied to a vertical well that targeted the Athel formation in the south of the Sultanate of Oman. This formation is made of silicilyte and is characterized by very low permeability of about 0.01 md on average. Thus, hydraulic fracturing plays a critical role for the economic oil recovery in this reservoir. Aiming to achieve a better zonal contribution, the stimulation design was changed from a limited number of large fractures to an extensive multistage fracturing design in the subject well. Sixteen hydraulic fracturing stages were planned. The DCHSA was applied to provide accurate and efficient fracture geometry evaluation. The DCHSA accurately identified fracture intervals and their corresponding heights and orientations. This enabled effectively determining the created fracture quality and helped explain the responses of the production logs that were recorded during the well test. This study provided a foundation for the placement and completion design of the future wells in the subject reservoir. It particularly revealed adequate fracturing intervals and the optimum number of stages required to achieve optimum reservoir coverage and avoid vertical overlapping.

Integrated Hydraulic Fracturing and Well-Test Data Analytics using R

1. Abstract Various datasets are generated during hydraulic fracturing, flowback- and well-testing operations, which require consistent integration to lead to high-quality well performance interpretations. An automated digital workflow has been created to integrate and analyze the data in a consistent manner using the open-source programming language R. This paper describes the workflow, and it explains how it automatically generates well performance models and how it analyzes raw diagnostic fracture injection test (DFIT) data using numerical algorithms and Machine Learning. This workflow is successfully applied in a concession area located in the center of the Sultanate of Oman, where to date a total of 25+ tight gas wells are drilled, hydraulically fractured and well-tested. It resulted in an automated and standardized way of working, which enabled identifying trends leading to improved hydraulic fracturing and well-testing practices.

Hydraulic Fracturing to Successfully Exploit Depleted Gas Reserves: A Case History from the North Sea

Hydraulic fracturing has been well established in the Southern North Sea (SNS) since the mid-1980s; however, it has typically been conducted as the final phase of development in new gas fields. One of these fields is Chiswick located in the Greater Markham area 90 miles offshore UK in 130 ft of water. Following an unsuccessful well repair of the multi-fractured horizontal well C4, it was decided to cost-effectively and expediently exploit the remaining pressure-depleted reserves near the toe via a single large fracture initiated from a deviated sidetrack wellbore designated C6. A deviated wellbore was chosen versus the original near-horizontal well to reduce well risk and costs and ultimately deliver an economic well. Several key challenges were identified, and mitigating measures were put in place. Modular formation dynamics tester data from the sidetrack open hole indicated the reservoir pressure gradient had depleted to 0.23 to 0.25 psi/ft, raising concerns about the ability of the well to unload the fluid volumes associated with a large fracture treatment. Wellbore deviation and azimuth with the associated potential for near-wellbore tortuosity would drive a typically short perforation interval (i.e., 3 ft). However, a compromise to mitigate convergent pressure loss in depletion was required, and the perforation interval was therefore set at 14 ft with provision made for robust step-down tests (SDT) and multi-mesh sand slugs. To further offset any near-well convergence pressure drop during cleanup, an aggressive tip screenout (TSO) proppant schedule, including a high concentration tail-in (12 PPA) with an aggressive breaker schedule, was executed to fully develop propped hydraulic width. Following formation breakdown and SDT to 40 bbl/min, the well went on near-instantaneous vacuum. Clearly, an extremely conductive feature had been created or contacted. However, upon use of a robust crosslinked gel formulation and 100-mesh sand, the bottomhole and positive surface pressure data allowed a suitable fracture design to be refined and placed with a large width, as evidenced by the extreme 2,309-psi net pressure development over that of the pad stage while placing 500,500 lbm of 16/30 resin-coated (RC) intermediate strength proppant (ISP) to 12 PPA. Although a lengthy nitrogen lift by coiled tubing (CT) was planned, the well cleanup response in fact allowed unaided hydrocarbon gas flow to surface within a short period. The well was then further beaned-up under well test conditions to a flow rate of approximately 26 MMscf/D under critical flowing conditions with a higher bottomhole flowing pressure than that of the original C4 well. Given the last producing rate of the original multiple fractured horizontal wellbore was 27 MMscf/D at a drawdown of 1,050 psi through two separate hydraulic fractures, then the outcome of this well was judged to be highly successful and at the limit of predrill expectations. This case history explains and details the rationale, methods, and techniques employed in well C6 to address the challenge of successful hydraulic fracture stimulation in a depleted formation. Challenges were addressed by combining a number of techniques, coupled with field experience, resulting in a highly productive well despite the relatively low reservoir pressure coupled with a limited time frame to plan and execute. These techniques are transferrable to other offshore gas fields in the region where reservoir depletion makes economic recovery difficult or indeed prohibitive.

An Integrated Approach of Numerical Well Test for Well Intersecting Fractures Based on FMI Image

Fluid flow is strongly affected by fractures in unconventional reservoirs. It is essential to deeply understand the flow characteristics with fractures for improving the production and efficiency of unconventional reservoir exploitation. The purpose of this work is to develop an accurate numerical model to evaluate the transient-pressure response for well intersecting fractures. The meshes generated from Fullbore Formation Micro-Imager (FMI) images ensure an efficient numerical description of the geometries for fractures and interlayers. The numerical simulation is implemented by an inhouse finite element method-based code and benchmarked with drill stem test (DST) data. The results show that three flow regimes appear in the reservoir with fractures within the test period: wellbore afterflow, pseudolinear flow, and radial flow. In contrast, only the wellbore afterflow and radial flow appear for the wells without fractures. The results also reveal that fractures dominate the flow near the wellbore. Verification and application of the model show the practicability of the integrated approach for investigating the transient-pressure behaviors in the unconventional reservoir.

Determination of Safe Mud Weight Window in Rumaila Oilfield, Southern Iraq

The oil and gas industry, wellbore instability plays an important role in financial losses and stops the operations while the drilling which leads to extra time known as non-productive time. In this work, a comprehensive study was carried out to realize the nature of the instability problems of the wellbore in Rumaila oilfield to improve the well design. The study goal is to develop a geomechanical model in one dimension by utilizing Schlumberger Techlog (Version 2015) software. Open hole wireline measurements were needed to develop the model. The model calibrating and validating with core laboratory tests (triaxial test), well test (Mini-frac test), repeated formation test. Mohr-Coulomb, Mogi-Coulomb, and Modified Lade are the three failure criteria which utilized to analyze the borehole breakouts and to determine the minimum mud weight needed for a stable wellbore wall. For more accuracy of the geomechanical model, the predicted profile of the borehole instability is compared with the actual failure of the borehole that is recorded by caliper log. The results of the analysis showed that the Mogi-Coulomb criteria are closer to the true well failure compared with the other two criteria and considered as the better criteria in predicting the rock failure in the Rumaila oilfield. The wellbore instability analysis revealed that the vertical and low deviated wells (less than 40º) is safer and more stable. While, the horizontal and directional wells should be drilled longitudinally to the direction of the minimum horizontal stresses at a range between 140º–150º North West-South East and the mud weight recommended is increased to 10.5 ppg to avoid most of instabilities problems. The results contribute in development plan of the wells nearby the studied area and decreasing NPT and cost.