Hydraulic Fracturing at Clair: Unlocking the Potential of Europe's Largest Reservoir

The Greater Clair area, Europe's largest oilfield, has two existing platforms, Clair Phase 1 and Clair Ridge, on production with future potential for a third platform targeting undeveloped Lower Clair Group to the South of Ph1. Clair Phase 1 was the initial development of Clair, targeting Lower Clair Group (LCG) reservoir consisting of a complex Devonian sandstone in six units. Most Phase 1 wells penetrated relatively good quality reservoir enhanced by natural fractures, while more recently Clair Ridge wells took a similar approach targeting natural fractures, however that strategy is continually being evaluated. In some areas however low matrix quality and lack of natural fractures were the dominant characteristics resulting in lower production rates. A brief comparison of the range of production outcomes will be presented, including potential downsides of reliance on natural fractures. Given the large oil volumes in areas of known poorer rock quality, alongside variable production results, a hydraulic fracturing trial was initiated in 2017. Well 206/08-A23 (A23) targeted previously under-developed, poor-quality Unit VI within the Phase 1 Graben area where natural fractures are absent. A pre-frac production test established baseline production of 900BOPD in December 2018. The A23 objectives included subsequent hydraulically fracturing the well to test this techniques ability to unlock production from tight, matrix dominated formation. Detailed analysis of core, log and limited vertical well fracturing data (from initial fracturing trials of 1980's vintage), yielded robust designs. Key challenges included overcoming very low KV/KH ratios with fracture heights exceeding 300ft. The resulting detailed designs provided consistent and predictable hydraulic fracturing execution in A23 in 2019, including placement of four planned 500klbs treatments combined with coil clean-outs after each stage to unload solids and fluids from the well. Initial fracture designs were conservative in terms of pad and proppant scheduling which, alongside learnings around operational logistics offshore West of Shetlands and completion design, offer significant optimisations for future hydraulic fractures. Post frac A23 became the highest non-natural fractured producer across Clair. Initially a six-fold production increase was observed with monitoring of transient production ongoing. Tracer analysis confirmed production contribution from all zones. Proving fracturing technology brings opportunities to unlock poorer Phase 1 and Ridge reservoir areas. Additionally, significant portions of the undeveloped Lower Clair Group to the South of Ph1 comprises lower permeability reservoir with higher viscosity oil and reduced natural fracture presence. Hydraulic fracturing is therefore a crucial completion technique for developing this lower quality reservoir and brings significant value enhancement to the project. Efficient delivery of numerous large fractures in a harsh offshore environment West of Shetlands presents significant challenges. The influence of how the A23 fracturing results and learnings are guiding future hydraulic fracturing concept are detailed, including optimising platform engineering design to facilitate efficient fracturing operations while maintaining the required productivity in this challenging scenario.

Numerical Simulation Study of Proppant Transport in Cross Fractures

Hydraulic fracturing technology is an important means to stimulate unconventional reservoirs, and the placement morphology of proppant in cross fractures is a key factor affecting the effect of hydraulic fracturing. It is very important to study the proppant transport law in cross fractures. In order to study the proppant transportation law in cross fractures, based on the CFD-DEM method, a proppant transport model in cross fractures was established. From the two aspects of the flow field in the fractures and the morphology of the proppant dune, the influence of the natural fracture approach angle, the fracturing fluid viscosity and injection rate on the proppant transport is studied. Based on the principle of hydropower similarity, the conductivity of proppant dune under different conditions is quantitatively studied. The results show that the natural fracture approach angle affects the distribution of proppant and fracturing fluid in natural fractures, and further affects the proppant placement morphology in hydraulic fractures and natural fractures. When the fracturing fluid viscosity is low and the displacement is small, the proppant forms a "high and narrow" dune at the entrance of the fracture. With the increase of the fracturing fluid viscosity and injection rate, the proppant settles to form a "short and wide" placement morphology. Compared with the natural fracture approach angle, the fracturing fluid viscosity and injection rate have a more significant impact on the conductivity of proppant dune. This paper investigated the proppant transportation in cross fractures, and quantitatively analyzes the conductivity of proppant dunes with different placement morphology. The results of this study can provide theoretical guidance for the design of hydraulic fracturing.

Prediction Analysis of Induced Casing Deformation of Horizontal Well in Fractured Shale Reservoir: Case Study of Weiyuan Gasfield

Casing deformation (CD) is a major challenge for shale gas development in Weiyuan gasfield, natural fracture (NF) slippage is one of the main causes of CD in Weiyuan gas filed. In order to study the mechanism and regularity of NF slippage induced CD, a wellbore shear stress calculation model and a CD degree prediction model are established. And results show that, the approach angle and ground principal stress difference have significant influence on wellbore shear stress, high wellbore shear stress occurs when wellbore orientation is perpendicular to the NF trend. Wellbore shear stress increases with the increase of fracture fluid pressure and NF area, improving casing strength or cementing quality has limited effect on reducing the risk of CD. The smaller the young's modulus, the higher the CD degree, Poisson's ratio has limited effect on CD degree. NF approach and fracture fluid pressure determines the value of CD degree. Field case shows that reasonable fracturing technology to control fracture net pressure and wellbore position arrangement are helpful for reducing CD risk, and the model proposed in this paper can be used to predict CD risk and calculate the CD degree.

A Case Study Using Integrated Multi Disciplinary Approach to Model Microseismic Events During Stimulation

For tight reservoirs where hydraulic fracturing is required to enable sufficient fluid mobility for economic production, it is critical to understand the placement of induced fractures, their connectivity, extent, and interaction with natural fractures within the system. Hydraulic fracture initiation and propagation mechanisms are greatly influenced by the effect of the stress state, rock fabric and pre-existing features (e.g. natural fractures, faults, weak bedding/laminations). A pre-existing natural fracture system can dictate the mode, orientation and size of the hydraulic fracture network. A better understanding of the fracture growth phenomena will enhance productivity and also reduce the environmental footprint as less fractures can be created in a much more efficient way. Assessing the role of natural fractures and their interaction with hydraulic fractures in order to account for them in the hydraulic fracture model is achieved by leveraging microseismicity. In this study, we have used a combination of borehole and surface microseismic monitoring to get high vertical resolution locations and source mechanisms. 3D numerical modelling of hydraulic fracturing in complex geological conditions to predict fracture propagation is essential. 3D hydraulic fracturing simulation includes modelling capabilities of stimulation parameters, true 3D fracture propagation with near wellbore 3D complexity including a coupled DFN and the associated microseismic event generation capability. A 3D hydraulic fracture model was developed and validated by matching model predictions to microseismic observations. Microseismic source mechanisms are leveraged to determine the location and geometry of pre-existing features. In this study, we simulate a DFN based on the recorded seismicity of multi stage hydraulic fractures in a horizontal well. The advanced 3D hydraulic fracture modelling software can integrate effectively and efficiently data from a variety of multi-disciplinary sources and scales to create a subsurface characterization of the unconventional reservoir. By incorporating data from 3D seismic, LWD/wireline, core, completion/stimulation monitoring, and production, the software generates a holistic reservoir model embedded in a modular, multi-physics software platform of coupled numerical solvers that capture the fundamental physics of the processes being modelled. This study illustrates the importance of a powerful software tool that captures the necessary physics of stimulation to predict the effects of various completion designs and thereby ensure the most accurate representation of an unconventional reservoir response to a stimulation treatment.

Applying Boundary Element Simulation to Generate Stress-Based Fracture Drivers: A Case Study for the Montney Unconventional Oil/Gas Play in Western Canada Basin

Natural fracture distribution is critical to the hydrocarbon production from the Early Triassic Montney unconventional oil and gas play. The formation underwent several tectonic events, creating a unique natural fracture system. Identifying tectonic events and their stress field evolution is an import component in fracture system modeling and prediction. The objective of this paper is to identify the evolution of paleo-stress domains, to establish related tectonic models, and to generate the drivers for fracture network modeling which will aid in reservoir understanding and overall play development. Compared with other geomechanical approaches, the boundary element method (BEM) is better suited for the structural characteristics in the study area. Hence, the corresponding boundary element simulation (BES) was applied for the evolution of the paleo-stress domains. The methodology is a combination of 3D BEM and Monte Carlo simulations. The inputs include seismic interpreted faults and natural fractures from Formation Microimager logs. After applying the methodology, several best fit realizations were calculated, and the admissible paleo-stress domains were characterized by the tectonic models which are consistent with the regional tectonic evolution of the formation. The study area is about 400 km2 located at northeast British Columbia in the Western Canada Basin. The main structural features are the thrust and back-thrust faults, forming different fault blocks without any significant deformation structures. The Montney formation within the study area underwent several tectonic events: (1) regime of terrane collision, indentation and lateral escape during end of Middle Jurassic to Middle Cretaceous; (2) regime of left-lateral transpression dominated by strike-slip during end of Late Cretaceous and Paleocene; and (3) regime of right-lateral transtension dominated by strike-slip during end of Early and Middle Eocene which is maintained till present day. Three major stress domains were identified in the study area by the application of the BES method, one reverse event and two strike-slip events, representing paleo and present-day stress domains. These stress domains are consistent with the regional tectonic evolution history of the foreland basin. The stress field parameters, such as stress ratio and maximum horizontal stress azimuth, are consistent. The derived tectonic models are shown to be reliable drivers for the subsequent fracture modeling and geomechanics study.

Hydraulic Fracturing Design Developments Towards Improved Effectiveness in North Kuwait Reservoirs

The North Kuwait Jurassic Gas (NKJG) reservoirs pose productivity challenges due to their geological heterogeneity, complex tectonic settings, high stress anisotropy, high pore pressure, and high bottom-hole temperature. Additionally, high natural fracture intensity in clustered areas play an important role in the wells hydrocarbon deliverability. These challenges are significant in field development starting from well design and stimulation up to production stages. The Gas Field Development Group (GFDG) are introducing for the first time in Kuwait new completion designs at high fracturing intensity; open-hole Multi Stage Completions (MSC), 4.5" Monobores and hybrid completions along with customized and efficient stimulation methods. This development strategy designed to overcome reservoir difficulties and enhance the well performance during initial testing and long-term production phases. At early stages of production, most of the wells were stimulated with simple matrix acidizing jobs and this method was sufficient to reach commercial production in conventional reservoirs. However, the reservoir depletion trend has negatively affected stimulation effectiveness and the wells performance in the recent years; hence, short and long-term solutions introduced to manage the sub-hydrostatic reservoir pressure. Our current focus is on the short-term stimulation solutions as they are relatively easier to apply compared to the long-term solutions that require additional resources, which are not available in the country. As a result, the stimulation methods, specifically the hydraulic fracturing treatments, increased production dramatically compared to previous years and it applied across North Kuwait Fields in conventional and unconventional reservoirs to reach the production targets of 2020-2021. The hydraulic fracturing treatment designs improved during the 2020-2021 fiscal year. The number of operations tripled compared to before and alternative chemical treatments with new fracturing designs implemented. In addition, these treatments executed across different well completions and reservoir properties. The objectives behind each fracturing treatment were different; for example: discovering new areas, re-stimulating under-performing wells, fracturing unconventional reservoirs, etc. Some promising wells did not flow as per expectation after matrix acidizing treatments despite the logs showing good reservoir quality similar to offset wells with good production. After re-stimulating with acid fracturing, the wells performed much better and one of them set a benchmark as the best producer amongst the offset wells. This paper evaluates the gaps in developing NKJG reservoirs, including fracturing treatments and highlights of the pros/cons for each operation, which in future supports the improvement of stimulation job designs. Moreover, it reveals the future requirements that control the operation success and how to reduce the well cleaning time post-fracturing in the event of low reservoir pressure. Finally, it describes how the outcome of the analyses directly assists reaching the production targets; since NKJG's production mainly depends on stimulation techniques.

Fracture Characterization in Deep Gas Reservoirs to Identify Fracture Enhanced Flow Units, Offshore Abu Dhabi

High gas flow rates in deep-buried dolomitized reservoir from an offshore field Abu Dhabi cannot be explained by the low matrix permeability. Previous permeability multiplier based on distance to major faults is not a solid geological solution due to over-simplifying reservoir geomechanics, overlooking folding-related fractures, and lack of detailed fault interpretation from poor seismic. Alternatively, to characterize the heterogeneous flow related with natural fractures in this undeveloped reservoir, fracture network is modelled based on core, bore hole imager (BHI), conventional logs, seismic data and test information. Limited by investigation scale, vertical wells record apparent BHI, and raw fracture interpretation cannot represent true 3D percolation reflected on PLT. To overcome this shortfall, correction based on geomechanics and mechanical layer (ML) analysis is performed. Young's modulus (E), Poisson ratio (ν), and brittleness index are calculated from logs, describing reservoir tendency of fracturing. Other than defining MLs, bedding plane intensity from BHI is also used as an indicator of fracture occurrence, since stress tends to release at strata discontinuity and forms bed-bounded fractures observed from cores. Subsequently, a new fracture intensity is generated from combined geomechanics properties and statistics average of BHI-derived fracture occurrence within the ML frame, which improves match with PLT and distinguishes fracture enhance flow intervals consistently in all wells. Seismic discontinuity attributes are used as static fracture footprints to distribute fractures from wells to 3D. The final hybrid DFN comprises large-scale deterministic zone-crossing fractures and small-scale stochastic bed-bounded fractures. Sub-vertical open fractures are dominated by NE-SW wrenching fractures related with Zagros compression and reactive salt upward movement. There is no angle rotation of fractures in different fault blocks. Open fractures in other strikes are supported by partial cements and mismatching fracture walls on computerized tomography (CT) images. ML correlation shows vertical consistence across stratigraphic framework and its intensity indicates fracture potential of vertical zones reflected by tests. Fracture-enhanced flow units are further constrained by a threshold in both combined geomechanics properties and statistics average of raw BHI fracture intensity in ML frame. As a result, final fracture network maps reservoir brittleness and flow potential both vertically and laterally, identifying fracture regions along folding axis not just major faults, evidenced by wells and seismic. According to the upscaling results, the case study reveals a type-III fractured reservoir, where fractures contribute to flow not to volume. Fracture network enhances bed-wise horizontal communication but also opens vertical feeding channels. Fracture permeability is mainly influenced by aperture and intensity, while aspect ratio, fracture length, and proportion of strikes and dips mainly influence permeability distribution rather than absolute values. This study provides a production-oriented characterization workflow of natural fracture heterogeneity based on correction of raw BHI in undeveloped fields.

Coupling Production Data with the DFN to Unravel Natural Fracture Networks in a Tight Naturally Fractured Jurassic Reservoir, Kuwait

The Jurassic Najmah-Sargelu of west Kuwait can be thought of as a "hybrid" between a conventional and an unconventional reservoir. These systems form an increasingly important resource for operators, but their performance is unpredictable because matrix permeability is in the micro-Darcy range and production depends on natural fractures. Success depends on how well the static models are aligned to the dynamic production, and the effectiveness of a fit-for-purpose multistage completion on project economics. In this work we present our lessons learnt in production modelling these reservoirs and the coupling between reservoir simulation and the discrete fracture network (DFN). Our reservoir models were constructed using a highly integrated approach incorporating data from all scales and disciplines (drilling, geophysical, geological, reservoir and production) and the production simulations were run using dual porosity and black oil models. As expected, the DFN played a key part of this effort. An iterative approach was used to adjust the DFN so that it was consistent with production observations. However, in all cases care was made to ensure the new DFN honoured the seismic, geological, well log and drilling data from which it was generated. Final, smaller adjustments were made to the simulation model at the log scale to match PLT data. We used uncertainty analysis to run hundreds of simulation cases and found that the character of the natural fractures is quite well imprinted in the observed production data, particularly pressure buildup data. This gave us a better understanding of whether the natural fractures are diffuse and laterally extensive away from the wellbore or if they are localized close to the wellbore. Where reservoir simulation history matches inferred laterally extensive natural fractures, an good correlation was obtained with the natural fracturing from the DFN. This correlation was poor where natural fracturing was confined to a smaller depth interval (as observed from PLT), and is a result of the limitation in seismic resolution to resolve these natural fractures. The lessons learnt from our work helps towards improved understanding of production mechanisms of these reservoirs and their natural fracture networks. This, together with higher resolution azimuthal seismic, advanced wellbore characterization data and multistage completions are the desired key ingredients for technically enhancing production in these reservoirs.