Application Investigation of Light Proppant in Hydraulic Fracturing of Marcellus Shale

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.

2021 ◽  
Author(s):  
Mohamed El Sgher ◽  
Kashy Aminian ◽  
Ameri Samuel

Abstract The objective of this study was to investigate the impact of the hydraulic fracturing treatment design, including cluster spacing and fracturing fluid volume on the hydraulic fracture properties and consequently, the productivity of a horizontal Marcellus Shale well with multi-stage fractures. The availability of a significant amount of advanced technical information from the Marcellus Shale Energy and Environment Laboratory (MSEEL) provided an opportunity to perform an integrated analysis to gain valuable insight into optimizing fracturing treatment and the gas recovery from Marcellus shale. The available technical information from a horizontal well at MSEEL includes well logs, image logs (both vertical and lateral), diagnostic fracture injection test (DFIT), fracturing treatment data, microseismic recording during the fracturing treatment, production logging data, and production data. The analysis of core data, image logs, and DFIT provided the necessary data for accurate prediction of the hydraulic fracture properties and confirmed the presence and distribution of natural fractures (fissures) in the formation. Furthermore, the results of the microseismic interpretation were utilized to adjust the stress conditions in the adjacent layers. The predicted hydraulic fracture properties were then imported into a reservoir simulation model, developed based on the Marcellus Shale properties, to predict the production performance of the well. Marcellus Shale properties, including porosity, permeability, adsorption characteristics, were obtained from the measurements on the core plugs and the well log data. The Quanta Geo borehole image log from the lateral section of the well was utilized to estimate the fissure distribution s in the shale. The measured and published data were utilized to develop the geomechnical factors to account for the hydraulic fracture conductivity and the formation (matrix and fissure) permeability impairments caused by the reservoir pressure depletion during the production. Stress shadowing and the geomechanical factors were found to play major roles in production performance. Their inclusion in the reservoir model provided a close agreement with the actual production performance of the well. The impact of stress shadowing is significant for Marcellus shale because of the low in-situ stress contrast between the pay zone and the adjacent zones. Stress shadowing appears to have a significant impact on hydraulic fracture properties and as result on the production during the early stages. The geomechanical factors, caused by the net stress changes have a more significant impact on the production during later stages. The cumulative gas production was found to increase as the cluster spacing was decreased (larger number of clusters). At the same time, the stress shadowing caused by the closer cluster spacing resulted in a lower fracture conductivity which in turn diminished the increase in gas production. However, the total fracture volume has more of an impact than the fracture conductivity on gas recovery. The analysis provided valuable insight for optimizing the cluster spacing and the gas recovery from Marcellus shale.


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.


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-18
Author(s):  
Jaeyoung Park ◽  
Candra Janova

This paper introduces a flow simulation-based reservoir modeling study of a two-well pad with long production history and identical completion parameters in the Midland Basin. The study includes building geologic model, history matching, well performance prediction, and finding optimum lateral well spacing in terms of oil volume and economic metrics. The reservoir model was constructed based on a geologic model, integrating well logs, and core data near the target area. Next, a sensitivity analysis was performed on the reservoir simulation model to better understand influential parameters on simulation results. The following history matching was conducted with the satisfactory quality, less than 10% of global error, and after the model calibration ranges of history matching parameters have substantially reduced. The population-based history matching algorithm provides the ensemble of the history-matched model, and the top 50 history-matched models were selected to predict the range of Estimate Ultimate Recovery (EUR), showing that P50 of oil EUR is within the acceptable range of the deterministic EUR estimates. With the best history-matched model, we investigated lateral well spacing sensitivity of the pad in terms of the maximum recovery volume and economic benefit. The results show that, given the current completion design, the well spacing tighter than the current practice in the area is less effective regarding the oil volume recovery. However, economic metrics suggest that the additional monetary value can be realized with 150% of current development assumption. The presented workflow provides a systematic approach to find the optimum lateral well spacing in terms of volume and economic metrics per one section given economic assumptions, and the workflow can be readily repeated to evaluate spacing optimization in other acreage.


Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-21
Author(s):  
Hongfei Ma ◽  
Wenqi Zhao ◽  
Meng Sun ◽  
Xiaodong Wang ◽  
Lun Zhao ◽  
...  

The volume fracturing technique has been widely used to improve the productivity of ultralow-permeability reservoirs. This paper presents a new semianalytical model to simulate the pressure transient and production behaviour of finite conductivity vertical fractured wells with stimulated reservoir volume (SRV) in heterogeneous reservoirs. The model is based on the five-linear flow model, the Warren-Root model, and fracture conductivity influence function. The model is validated by comparing its results with a numerical model. One novelty of this model is its consideration of three different kinds of production prediction models. Constant rate, constant pressure, and compound working systems are taken into account. This paper illustrates the effects of the SRV size and shape, mobility ratio, initial flow rate, limiting wellbore pressure, and hydraulic fracture parameters under different working systems. Results show that the SRV and parameters of fractures have a significant influence on long-term well performance. Moreover, the initial rate can extend the constant rate period by 418%, and limiting wellbore pressure can effectively improve the cumulative recovery rate by 23%. Therefore, this model can predict long-term wells’ behaviour and provide practical guiding significance for hydraulic fracturing design.


Geofluids ◽  
2018 ◽  
Vol 2018 ◽  
pp. 1-12 ◽  
Author(s):  
Jing Huang ◽  
Lan Ren ◽  
Jinzhou Zhao ◽  
Zhiqiang Li ◽  
Junli Wang

Refracturing is an encouraging way to uplift gas flow rate and ultimate gas recovery from shale gas wells. A numerical model, considering the stimulated reservoir volume and multiscale gas transport, is applied to simulate the gas production from a refractured shale gas well. The model is verified against field data from a shale gas reservoir in Sichuan Basin. Two refracturing scenarios: refracturing through existing perforation clusters and refracturing through new perforation zones, are included in the simulation work. Three years after production is determined to be the optimum time for refracturing based on the evolution analysis of reservoir pressure, effective stress, fracture permeability, and gas recovery. The role that the hydraulic fracture conductivity and hydraulic fracture half-length play in gas production for different refracturing cases is explored. Pumping parameters of the refracturing job in Sichuan Basin are discussed combining with sensitivity analysis, and suggestions for pumping parameters optimization are proposed.


Energies ◽  
2021 ◽  
Vol 14 (19) ◽  
pp. 6007 ◽  
Author(s):  
Christopher R. Clarkson ◽  
Zhenzihao Zhang ◽  
Farshad Tabasinejad ◽  
Daniela Becerra ◽  
Amin Ghanizadeh

The current practice for multi-fractured horizontal well development in low-permeability reservoirs is to complete the full length of the well with evenly spaced fracture stages. Given methods to evaluate along-well variability in reservoir quality and to predict stage-by-stage performance, it may be possible to reduce the number of stages completed in a well without a significant sacrifice in well performance. Provision and demonstration of these methods is the goal of the current two-part study. In Part 1 of this study, reservoir and completion quality were evaluated along the length of a horizontal well in the Montney Formation in western Canada. In the current (Part 2) study, the along-well reservoir property estimates are first used to forecast per-stage production variability, and then used to evaluate production performance of the well when fewer stages are completed in higher quality reservoir. A rigorous and fast semi-analytical model was used for forecasting, with constraints on fracture geometry obtained from numerical model history matching of the studied Montney well flowback data. It is concluded that a significant reduction in the number of stages from 50 (what was implemented) to less than 40 could have yielded most of the oil production obtained over the forecast period.


SPE Journal ◽  
2016 ◽  
Vol 21 (03) ◽  
pp. 1028-1040 ◽  
Author(s):  
Wei Yu ◽  
Kan Wu ◽  
Kamy Sepehrnoori

Summary Two key technologies such as horizontal drilling and hydraulic fracturing have led to the economic production of unconventional resources such as shale gas and tight oil. In reality, a nonplanar hydraulic-fracture geometry with varying fracture width and fracture permeability is created during the hydraulic-fracturing process. However, it is challenging to simulate well performance from the nonplanar fracture geometry. For the sake of simplicity, the nonplanar fracture geometry is often represented by ideal planar fracture geometry with constant fracture width, which one can easily handle analytically, semianalytically, and numerically. However, such ideal fracture geometry is inadequate to capture the physics of the transient flow behavior of the nonplanar fracture geometry. Although significant efforts were made in recent years to numerically model well performance from the complex fracture geometry, these approaches are still challenging to model the nonplanar fracture geometry with varying width because of a large considerable fracture-gridding issues and an expensive computation cost. In addition, the effect of nonplanar fracture geometry on well productivity and transient flow behavior was not reported in the literature. Hence, a model to handle the nonplanar fracture geometry by considering varying fracture width and fracture permeability is still lacking in the petroleum industry. Zhou et al. (2014) proposed a semianalytical model to handle the complex fracture geometry with constant fracture width. However, the semianalytical model did not consider the effects of stress-dependent fracture conductivity and the nonplanar fracture geometry as well as planar fracture geometry with varying fracture width along the fracture. In this work, we extended the semianalytical model to simulate production from the nonplanar fracture geometry as well as planar fracture geometry with varying width. In addition, the effect of stress-dependent fracture conductivity was implemented in the model. We verified the semianalytical model against a numerical reservoir simulator for single planar fracture with constant width. We performed two case studies. The first case contains a comparison of two planar fractures, one with constant fracture width and another with varying fracture width. In the second study, we compared two fractures with different fracture geometries such as planar fracture geometry and nonplanar fracture geometry, which were generated from the fracture-propagation model. In addition, transient flow regimes were investigated on the basis of a log-log graph of the dimensionless pressure drop and pressure-drop derivative vs. the dimensionless time. This work can provide critical insights into understanding the well performance from tight oil reservoirs with the nonplanar hydraulic-fracture geometry.


2021 ◽  
Author(s):  
Fengyuan Zhang ◽  
Hamid Emami-Meybodi

Abstract This study presents a new type-curve method to characterize hydraulic fracture (HF) attributes and dynamics by analyzing two-phase flowback data from multi-fractured horizontal wells (MFHWs) in hydrocarbon reservoirs.The proposed method includes a semianalytical model, as well as a workflow to estimate HF properties (i.e., initial fracture pore-volume and fracture permeability) and HF closure dynamics (through iterating fracture compressibility and permeability modulus).The semianalytical model considers the coupled two-phase flow in the fracture and matrix system, the variable production rate at the well, as well as the pressure-dependent reservoir and fluid properties. By incorporating the contribution of fluid influx from matrix into the fracture effective compressibility, a new set of dimensionless groups is defined to obtain a dimensionless solution for type-curve analysis. The accuracy of the proposed method is tested using the synthetic data generated from six numerical simulation cases for shale gas and oil reservoirs. The numerical validation confirms the unique behavior of type curves during fracture boundary dominated flow and verifies the accuracy of the type-curve analysis in the characterization of fracture properties. For field application, the proposed method is applied to two MFHWs in Marcellus shale gas and Eagle Ford shale oil.The agreement of interpreted results between the proposed method and straight-line analysis not only demonstrates the practicality in field application but also illustrates the superiority of the type-curve method as an easy-to-use technique to analyze two-phase flowback data. The analysis results from both of the field examples reveal the consistency in the estimated fracture properties between the proposed method and long-term history matching.


Energies ◽  
2021 ◽  
Vol 14 (11) ◽  
pp. 3137
Author(s):  
Amine Tadjer ◽  
Reider B. Bratvold ◽  
Remus G. Hanea

Production forecasting is the basis for decision making in the oil and gas industry, and can be quite challenging, especially in terms of complex geological modeling of the subsurface. To help solve this problem, assisted history matching built on ensemble-based analysis such as the ensemble smoother and ensemble Kalman filter is useful in estimating models that preserve geological realism and have predictive capabilities. These methods tend, however, to be computationally demanding, as they require a large ensemble size for stable convergence. In this paper, we propose a novel method of uncertainty quantification and reservoir model calibration with much-reduced computation time. This approach is based on a sequential combination of nonlinear dimensionality reduction techniques: t-distributed stochastic neighbor embedding or the Gaussian process latent variable model and clustering K-means, along with the data assimilation method ensemble smoother with multiple data assimilation. The cluster analysis with t-distributed stochastic neighbor embedding and Gaussian process latent variable model is used to reduce the number of initial geostatistical realizations and select a set of optimal reservoir models that have similar production performance to the reference model. We then apply ensemble smoother with multiple data assimilation for providing reliable assimilation results. Experimental results based on the Brugge field case data verify the efficiency of the proposed approach.


2021 ◽  
Author(s):  
Obinna Somadina Ezeaneche ◽  
Robinson Osita Madu ◽  
Ishioma Bridget Oshilike ◽  
Orrelo Jerry Athoja ◽  
Mike Obi Onyekonwu

Abstract Proper understanding of reservoir producing mechanism forms a backbone for optimal fluid recovery in any reservoir. Such an understanding is usually fostered by a detailed petrophysical evaluation, structural interpretation, geological description and modelling as well as production performance assessment prior to history matching and reservoir simulation. In this study, gravity drainage mechanism was identified as the primary force for production in reservoir X located in Niger Delta province and this required proper model calibration using variation of vertical anisotropic ratio based on identified facies as against a single value method which does not capture heterogeneity properly. Using structural maps generated from interpretation of seismic data, and other petrophysical parameters from available well logs and core data such as porosity, permeability and facies description based on environment of deposition, a geological model capturing the structural dips, facies distribution and well locations was built. Dynamic modeling was conducted on the base case model and also on the low and high case conceptual models to capture different structural dips of the reservoir. The result from history matching of the base case model reveals that variation of vertical anisotropic ratio (i.e. kv/kh) based on identified facies across the system is more effective in capturing heterogeneity than using a deterministic value that is more popular. In addition, gas segregated fastest in the high case model with the steepest dip compared to the base and low case models. An improved dynamic model saturation match was achieved in line with the geological description and the observed reservoir performance. Quick wins scenarios were identified and this led to an additional reserve yield of over 1MMSTB. Therefore, structural control, facies type, reservoir thickness and nature of oil volatility are key forces driving the gravity drainage mechanism.


Sign in / Sign up

Export Citation Format

Share Document