Efficient Modeling of Depletion Induced Fracture Deformation in Unconventional Reservoirs

2021 ◽  
Author(s):  
Yanli Pei ◽  
Kamy Sepehrnoori

Abstract The change of fracture conductivity during reservoir depletion significantly affects the well performance and stress evolution in unconventional formations. A common practice is to model fracture deformation using the traditional finite element method with very dense unstructured grids representing complex fracture geometries. However, the associated computational cost is high, so previous studies mainly use empirical correlations to catch the fracture conductivity loss or neglect fracture deformation during the production period. This work proposes a novel coupled flow and geomechanics model with embedded fracture methods to capture the fracture deformation accurately yet efficiently in unconventional reservoirs. Under a single set of structured grids, an embedded discrete fracture model (EDFM) is employed to characterize fluid flow through discrete fractures by introducing non-neighboring connections, and an extended finite element method (XFEM) is applied to simulate discontinuities over fracture walls by adding phantom nodes. In addition, a modified proppant model is incorporated to represent interactions between proppants and hydraulic surfaces, and an iterative coupling scheme is implemented to link the fracture-related fluid flow and solid mechanics. Being validated against the classical benchmark problem, the coupled model is used to investigate the impacts of proppant strength, closure stress, and bottomhole pressure on fracture deformation, well production, and in-situ stresses. Numerical results indicate that weaker proppant support induces more fracture aperture and production losses, resulting in greater stress changes and higher residual pressure in the depletion region. In comparison, the fracture deformation for a well-propped scenario is modest and barely affects the well performance and stress redistribution. Less stressed formation corresponds to lower closure stress on fracture walls, which triggers limited fracture closure and stabilizes well production. Moreover, a moderate bottomhole pressure decline rate avoids significant fracture closure while preserves relatively high initial production rates. The coupled flow and geomechanics model with embedded fracture methods resolves computational difficulties in modeling complex fracture deformations and delivers more insights on production forecast and stress changes crucial to refracturing and infill operations.

Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-13
Author(s):  
Kailei Liu ◽  
Boyue Xu ◽  
Changjea Kim ◽  
Jing Fu

Horizontal well fracturing technology is widely used in unconventional reservoirs such as tight or shale oil and gas reservoirs. Meanwhile, the potential of enhanced oil recovery (EOR) methods including huff-n-puff miscible gas injection are used to further increase oil recovery in unconventional reservoirs. The complexities of hydraulic fracture properties and multiphase flow make it difficult and time-consuming to understand the well performance (i.e., well production) in fractured shale reservoirs, especially when using conventional numerical methods. Therefore, in this paper, two methods are developed to bridge this gap by using the machine learning technique to forecast well production performance in unconventional reservoirs, especially on the EOR pilot projects. The first method is the artificial neural network, through which we can analyze the big data from unconventional reservoirs to understand the underlying patterns and relationships. A bunch of factors is contained such as hydraulic fracture parameters, well completion, and production data. Then, feature selection is performed to determine the key factors. Finally, the artificial neural network is used to determine the relationship between key factors and well production performance. The second is time series analysis. Since the properties of the unconventional reservoir are the function of time such as fluid properties and reservoir pressure, it is quite suitable to apply the time series analysis to understand the well production performance. Training and test data are from over 10000 wells in different fractured shale reservoirs, including Bakken, Eagle Ford, and Barnett. The results demonstrate that there is a good match between the available and predicated well performance data. The overall R values of the artificial neural network and time series analysis are both above 0.8, indicating that both methods can provide reliable results for the prediction of well performance in fractured shale reservoirs. Especially, when dealing with the EOR field cases, such as huff-n-puff miscible gas injection, Time series analysis can provide more accurate results than the artificial neural network. This paper presents a thorough analysis of the feasibility of machine learning in multiple fractured shale reservoirs. Instead of using the time-consuming numerical methods, it also provides a more robust way and meaningful reference for the evaluation of the well performance.


2021 ◽  
Author(s):  
Dimitry Chuprakov ◽  
Ludmila Belyakova ◽  
Ivan Glaznev ◽  
Aleksandra Peshcherenko

Abstract We developed a high-resolution fracture productivity calculator to enable fast and accurate evaluation of hydraulic fractures modeled using a fine-scale 2D simulation of material placement. Using an example of channel fracturing treatments, we show how the productivity index, effective fracture conductivity, and skin factor are sensitive to variations in pumping schedule design and pulsing strategy. We perform fracturing simulations using an advanced high-resolution multiphysics model that includes coupled 2D hydrodynamics with geomechanics (pseudo-3D, or P3D, model), 2D transport of materials with tracking temperature exposure history, in-situ kinetics, and a hindered settling model, which includes the effect of fibers. For all simulated fracturing treatments, we accurately solve a problem of 3D planar fracture closure on heterogenous spatial distribution of solids, estimate 2D profiles of fracture width and stresses applied to proppants, and, as a result, obtain the complex and heterogenous shape of fracture conductivity with highly conductive cells owing to the presence of channels. Then, we also evaluate reservoir fluid inflows from a reservoir to fracture walls and further along a fracture to limited-size wellbore perforations. Solution of a productivity problem at the finest scale allows us to accurately evaluate key productivity characteristics: productivity index, dimensional and dimensionless effective conductivity, skin factor, and folds of increase, as well as the total production rate at any day and for any pressure drawdown in a well during well production life. We develop a workflow to understand how productivity of a fracture depends on variation of the pumping schedule and facilitate taking appropriate decisions about the best job design. The presented workflow gives insight into how new computationally efficient methods can enable fast, convenient, and accurate evaluation of the material placement design for maximum production with cost-saving channel fracturing technology.


2021 ◽  
Author(s):  
Christian Windisch

Abstract This paper presents a holistic approach to modern oilfield and well surveillance through the inclusion of state-of-the-art edge computing applications in combination with a novel type of data transmission technology and algorithms developed in-house for automatic condition monitoring of SRP systems. The objective is to enable the responsible specialist staff to focus on the most important decisions regarding oilfield management, rather than wasting time with data collection and preparation. An own operated data communication system, based on LPWAN-technology transfers the dyno-cards, generated by an electric load cell, into the in-house developed production assistance software platform. Suitable programmed AI-algorithms enable automatic condition detection of the incoming dyno cards, including conversion and analysis of the corresponding subsurface dynamograms. A smart alarming system informs about occurring failure conditions and specifies whether an incident of rod rupture, pump-off condition, gas lock or paraffin precipitation occurred in the well. A surface mounted measuring device delivers liquid level and bottomhole pressure information automatically into the software. Based on these diverse data, the operations team plans the subsequent activities. The holistic application approach is illustrated using the case study of an SPR-operated well in an Austrian brownfield.


2018 ◽  
Vol 18 (3) ◽  
pp. 323-337
Author(s):  
Nguyen Huu Truong

Kinh Ngu Trang oilfield is of the block 09-2/09 offshore Vietnam, which is located in the Cuu Long basin, the distance from that field to Port of Vung Tau is around 140 km and it is about 14 km from the north of Rang Dong oilfield of the block 15.2, and around 50 km from the east of White Tiger in the block 09.1. That block accounts for total area of 992 km2 with the average water depth of around 50 m to 70 m. The characteristic of Oligocene E reservoir is tight oil in sandstone, very complicated with complex structure. Therefore, the big challenges in this reservoir are the low permeability and the low porosity of around 0.2 md to less than 1 md and 1% to less than 13%, respectively, leading to very low fracture conductivity among the fractures. Through the Minifrac test for reservoir with reservoir depth from 3,501 mMD to 3,525 mMD, the total leak-off coefficient and fracture closure pressure were determined as 0.005 ft/min0.5 and 9,100 psi, respectively. To create new fracture dimensions, hydraulic fracturing stimulation has been used to stimulate this reservoir, including proppant selection and fluid selection, pump power requirement. In this article, the authors present optimisation of hydraulic fracturing design using unified fracture design, the results show that optimum fracture dimensions include fracture half-length, fracture width and fracture height of 216 m, 0.34 inches and 31 m, respectively when using proppant mass of 150,000 lbs of 20/40 ISP Carbolite Ceramic proppant.


2021 ◽  
Author(s):  
Ahmed Farid Ibrahim ◽  
Mazher Ibrahim ◽  
Matt Sinkey ◽  
Thomas Johnston ◽  
Wes Johnson

Abstract Multistage hydraulic fracturing is the common stimulation technique for shale formations. The treatment design, formation in-situ stress, and reservoir heterogeneity govern the fracture network propagation. Different techniques have been used to evaluate the fracture geometry and the completion efficiency including Chemical Tracers, Microseismic, Fiber Optics, and Production Logs. Most of these methods are post-fracture as well as time and cost intensive processes. The current study presents the use of fall-off data during and after stage fracturing to characterize producing surface area, permeability, and fracture conductivity. Shut-in data (15-30 minutes) was collected after each stage was completed. The fall-off data was processed first to remove the noise and water hammer effects. Log-Log derivative diagnostic plots were used to define the flow regime and the data were then matched with an analytical model to calculate producing surface area, permeability, and fracture conductivity. Diagnostic plots showed a unique signature of flow regimes. A long period of a spherical flow regime with negative half-slope was observed as an indication for limited entry flow either vertically or horizontally. A positive half-slope derivative represents a linear flow regime in an infinitely conductive tensile fracture. The quarter-slope derivative was observed in a bilinear flow regime that represents a finite conductivity fracture system. An extended radial flow regime was observed with zero slope derivative which represents a highly shear fractured network around the wellbore. For a long fall-off period, formation recharge may appear with a slope between unit and 1.5 slopes derivative, especially in over-pressured dry gas reservoirs. Analyzing fall-off data after stages are completed provides a free and real-time investigation method to estimate the fracture geometry and a measure of completion efficiency. Knowing the stage properties allows the reservoir engineer to build a simulation model to forecast the well performance and improve the well spacing.


2022 ◽  
Author(s):  
John E. Busteed ◽  
Jesus Arroyo ◽  
Francisco Morales ◽  
Mohammed Omer ◽  
Francisco E. Fragachan

Abstract Uniformly distributing proppant inside fractures with low damage on fracture conductivity is the most important index of successful fracturing fluids. However, due to very low proppant suspension capacity of slickwater and friction reducers fracturing fluids and longer fracture closure time in nano & pico darcies formations, proppants settles quickly and accumulates near wellbore resulting in worse-than-expected well performance, as the fracture full capacity is not open and contributing to production. Traditionally, cross-linked polymer fluid systems are capable to suspend and transport high loading of proppants into a hydraulically generated fracture. Nevertheless, amount of unbroken cross-linked polymers is usually left in fractures causing damage to fracture proppant conductivity, depending on polymer loading. To mitigate these challenges, a low viscosity-engineered-fluid with excellent proppantcarrying capacity and suspension-in excess of 30 hours at static formation temperature conditions - has been designed, enhancing proppant placement and distribution within developed fractures, with a 98% plus retained conductivity. In this work experimental and numerical tests are presented together with the path followed in developing a network of packed structures from polymer associations providing low viscosity and maximum proppant suspension. Challenges encountered during field injection with friction are discussed together with the problem understanding characterized via extensive friction loop tests. Suspension tests performed with up to 8-10 PPA of proppant concentration at temperature conditions are shared, together with slot tests performed. Physics-based model results from a 3D Discrete Fracture Network simulator that computes viscosity, and elastic parameters based on shear rate, allows to estimate pressure losses along the flow path from surface lines, tubular goods, perforations, and fracture. This work will demonstrate the advanced capabilities and performance of the engineered fluid over conventional fracturing fluids and its benefits. Additionally, this paper will present field injection pressure analysis performed during the development of this fluid, together with a field case including production results after 8 months of treatment. The field case production decline observed after fracture treatment demonstrates the value of this system in sustaining well production and adding additional reserves.


2021 ◽  
Author(s):  
Aymen Alhemdi ◽  
Ming Gu

Abstract Slickwater-sand fracturing design is widely employed in Marcellus shale. The slickwater- sand creates long skinny fractures and maximizes the stimulated reservoir volume (SRV). However, due to the fast settling of sand in the water, lots of upper and deeper areas are not sufficiently propped. Reducing sand size can lead to insufficient fracture conductivity. This study proposes to use three candidate ultra-lightweight proppants ULWPs to enhance the fractured well performance in unconventional reservoirs. In step 1, the current sand pumping design is input into an in-house P3D fracture propagation simulator to estimate the fracture geometry and proppant concentrations. Next, the distribution of proppant concentration converts to conductivity and then to fracture permeability. In the third step, the fracture permeability from the second step is input into a reservoir simulator to predict the cumulative production for history matching and calibration. In step 4, the three ULWPs are used to replace the sand in the frac simulator to get new frac geometry and conductivity distribution and then import them in reservoir model for production evaluation. Before this study, the three ULWPs have already been tested in the lab to obtain their long-term conductivities under in-situ stress conditions. The conductivity distribution and production performance are analyzed and investigated. The induced fracture size and location of the produced layer for the current target well play a fundamental effect on ultra-light proppant productivity. The average conductivity of ULWPs with mesh 40/70 is larger and symmetric along the fracture except for a few places. However, ULWPs with mesh 100 generates low average conductivity and create a peak conductivity in limited areas. The ULW-3 tends to have less cumulative production compared with the other ULWPs. For this Marcellus Shale study, the advantages of ultra-lightweight proppant are restricted and reduced because the upward fracture height growth is enormous. And with the presence of the hydrocarbon layer is at the bottom of the fracture, making a large proportion of ULWPs occupies areas that are not productive places. The current study provides a guidance for operators in Marcellus Shale to determine (1) If the ULWP can benefit the current shale well treated by sand, (2) what type of ULWP should be used, and (3) given a certain type of ULWP, what is the optimum pumping schedule and staging/perforating design to maximize the well productivity. The similar workflow can be expanded to evaluate the economic potential of different ULWPs in any other unconventional field.


SPE Journal ◽  
2020 ◽  
Vol 25 (06) ◽  
pp. 3280-3299
Author(s):  
Hongyang Chu ◽  
Xinwei Liao ◽  
Zhiming Chen ◽  
W. John John Lee

Summary Because of readily available production data, rate-transient analysis (RTA) is an important method to predict productivity and reserves, and for reservoir and completion characterization in unconventional reservoirs. In addition, multihorizontal well pads are a common development method for unconventional reservoirs. Close well spacing between multifractured horizontal wells (MFHWs) in the multiwell pads makes interference from adjacent MFHWs especially significant. For RTA of production data from multihorizontal well pads, the influence of adjacent MFHWs cannot be ignored. In this work, we propose a semianalytic RTA model for the multihorizontal well pad with arbitrary multiple MFHW properties and starting-production times. Combining Laplace transformation and finite-difference analysis, we obtained a general solution of a multiwell mathematical model to use in RTA. Our model is applicable to cases of multiple MFHWs with different bottomhole pressures (BHPs), varying hydraulic-fracture properties, and different starting-production times. In the solutions, we observe bilinear flow, linear flow, transition flow, and multi-MFHW flow. Rate-normalized pressure (RNP) and its derivative are also affected by multi-MFHW flow. Two case studies revealed that the negative effect of interwell interference on the parent-well productivity is closely related to the pressure distribution caused by the production of child wells.


Sign in / Sign up

Export Citation Format

Share Document