Modeling fracture propagation, seismic identification, and gas production in shale reservoirs

Author(s):  
Naddia D. Arenas ◽  
Gabriela B. Savioli ◽  
Patricia M. Gauzellino ◽  
Juan E. Santos
2015 ◽  
Vol 04 (02) ◽  
pp. 23-32 ◽  
Author(s):  
Gary Feast ◽  
S. Sina Hosseini Boosari ◽  
Kim Wu ◽  
John Walton ◽  
Zufang Cheng ◽  
...  

2019 ◽  
Vol 8 (2S11) ◽  
pp. 2726-2737

Unconventional gas reservoirs are now the targets for meeting the demand for gas. These reservoirs are at the depth of more than 10,000 ft (even over 15000 depth as well) and are difficult to be exploited by conventional methods. For the last decades hydraulic fracturing has become the tool to develop these resources. Mathematical models (2D and pseudo-3D) have been developed for fracture geometry, which should be realistically created at the depth by surface controllable treatment parameters. If the reservoir rock is sandstone, then proppant fracturing is suitable and if the rock is carbonates, then acid fracturing is applicable. In both cases, proper design of controllable treatment parameters within constraints is essential. This needs proper optimization model which gives real controllable parametric vales. The model needs the most important analyses from geomechanical study and linear elastic fracture mechanics of rock containing unconventional gas so that fracture geometry makes maximum contact with the reservoirs for maximum recovery. Currently available software may lack proper optimization scheme containing geomechanical stress model, fracture geometry, natural fracture interactions, real field constraints and proper reservoir engineering model of unconventional gas resources, that is, production model from hydraulically fractured well (vertical and horizontal). An optimization algorithm has been developed to integrate all the modules, as mentioned above, controllable parameters, field constraints and production model with an objective function of maximum production (with or without minimization of treatment cost). Optimization is basically developed based on Direct Search Genetic and Polytope algorithm, which can handle dual objective function, non-differentiable equations, discontinuity and non-linearity. A dual objective function will meet operator’s economic requirements and investigate conflict between two objectives. The integrated model can be applied to a vertical or horizontal well in tight gas or ultra-tight shale gas deeper than over 10,000 ft. A simulation (with industrial simulators) was conducted to investigate and analyse fracture propagation behavior, under varying parameters with respect to the fracture design process, for tight gas reservoirs. Results indicate that hydraulic fracture propagation behavior is not uninhibited in deep reservoirs as some may believe that minor variations of variables such as in-situ stress, fluid properties etc. are often detrimental to fracture propagation in some conditions. Application of this model to a hypothetical tight and ultra-tight unconventional gas formations indicates a significant gas production at lower treatment cost; whereas the resources do not flow without any stimulation (hydraulic fracturing).


2016 ◽  
Vol 19 (03) ◽  
pp. 415-428 ◽  
Author(s):  
Najeeb S. Alharthy ◽  
Tadesse W. Teklu ◽  
Thanh N. Nguyen ◽  
Hossein Kazemi ◽  
Ramona M. Graves

Summary Understanding the mechanism of multicomponent mass transport in the nanopores of unconventional reservoirs, such as Eagle Ford, Niobrara, Woodford, and Bakken, is of great interest because it influences long-term economic development of such reservoirs. Thus, we began to examine the phase behavior and flow characteristics of multicomponent flow in primary production in nanoporous reservoirs. Besides primary recovery, our long-term objectives included enhanced oil production from such reservoirs. The first step was to evaluate the phase behavior in nanopores on the basis of pore-size distribution. This was motivated because the physical properties of hydrocarbon components are affected by wall proximity in nanopores as a result of van der Waals molecular interactions with the pore walls. For instance, critical pressure and temperature of hydrocarbon components shift to lower values as the nanopore walls become closer. In our research, we applied this kind of critical property shift to the hydrocarbon components of two Eagle Ford fluid samples. Then, we used the shifted phase characteristics in dual-porosity compositional modeling to determine the pore-to-pore flow characteristics, and, eventually, the flow behavior of hydrocarbons to the wells. In the simulation, we assigned three levels of phase behavior in the matrix and fracture pore spaces. In addition, the flow hierarchy included flow from matrix (nano-, meso-, and macropores) to macrofractures, from macrofractures to a hydraulic fracture (HF), and through the HF to the production well. From the simulation study, we determined why hydrocarbon fluids flow so effectively in ultralow-permeability shale reservoirs. The simulation also gave credence to the intuitive notion that favorable phase behavior (phase split) in the nanopores is one of the major reasons for production of commercial quantities of light oil and gas from shale reservoirs. It was determined that the implementation of confined-pore and midconfined-pore phase behavior lowers the bubblepoint pressure, and this, in turn, leads to a slightly higher oil recovery and lesser gas recovery. Also it was determined that the implementation of midconfined-pore and confined-pore phase-behavior shift reduces the retrograde liquid-condensation region, which in turn, leads to lower liquid yield while maintaining the same gas-production quantity. Finally, the important reason that we are able to produce shale reservoirs economically is “rubblizing” the reservoir matrix near HFs, which creates favorable permeability pathways to improve reservoir drainage. This is why multistage hydraulic fracturing is so critical for successful development of shale reservoirs.


Energies ◽  
2019 ◽  
Vol 12 (18) ◽  
pp. 3405 ◽  
Author(s):  
Jamiu M. Ekundayo ◽  
Reza Rezaee

The true contribution of gas desorption to shale gas production is often overshadowed by the use of adsorption isotherms for desorbed gas calculations on the assumption that both processes are identical under high pressure, high temperature conditions. In this study, three shale samples were used to study the adsorption and desorption isotherms of methane at a temperature of 80 °C, using volumetric method. The resulting isotherms were modeled using the Langmuir model, following the conversion of measured excess amounts to absolute values. All three samples exhibited significant hysteresis between the sorption processes and the desorption isotherms gave lower Langmuir parameters than the corresponding adsorption isotherms. Langmuir volume showed positive correlation with total organic carbon (TOC) content for both sorption processes. A compositional three-dimensional (3D), dual-porosity model was then developed in GEM® (a product of the Computer Modelling Group (CMG) Ltd., Calgary, AB, Canada) to test the effect of the observed hysteresis on shale gas production. For each sample, a base scenario, corresponding to a “no-sorption” case was compared against two other cases; one with adsorption Langmuir parameters (adsorption case) and the other with desorption Langmuir parameters (desorption case). The simulation results showed that while gas production can be significantly under-predicted if gas sorption is not considered, the use of adsorption isotherms in lieu of desorption can lead to over-prediction of gas production performances.


Energies ◽  
2020 ◽  
Vol 14 (1) ◽  
pp. 2
Author(s):  
Liang Gong ◽  
Yuan Zhang ◽  
Na Li ◽  
Ze-Kai Gu ◽  
Bin Ding ◽  
...  

The rapid growth in energy consumption and environmental pollution have greatly stimulated the exploration and utilization of shale gas. The injection of gases such as CO2, N2, and their mixture is currently regarded as one of the most effective ways to enhance gas recovery from shale reservoirs. In this study, molecular simulations were conducted on a kaolinite–kerogen IID composite shale matrix to explore the displacement characteristics of CH4 using different injection gases, including CO2, N2, and their mixture. The results show that when the injection pressure was lower than 10 MPa, increasing the injection pressure improved the displacement capacity of CH4 by CO2. Correspondingly, an increase of formation temperature also increased the displacement efficiency of CH4, but an increase of pore size slightly increased this displacement efficiency. Moreover, it was found that when the proportion of CO2 and N2 was 1:1, the displacement efficiency of CH4 was the highest, which proved that the simultaneous injection of CO2 and N2 had a synergistic effect on shale gas production. The results of this paper will provide guidance and reference for the displacement exploitation of shale gas by injection gases.


2015 ◽  
Author(s):  
Wu Kan ◽  
Jon E. Olson

Abstract Complex fracture networks have become more evident in shale reservoirs due to the interaction between pre-existing natural and hydraulic fractures. Accurate characterization of fracture complexity plays an important role in optimizing fracturing design, especially for shale reservoirs with high-density natural fractures. In this study, we simulated simultaneous multiple fracture propagation within a single fracturing stage using a complex hydraulic fracture development model. The model was developed to simulate complex fracture propagation by coupling rock mechanics and fluid mechanics. A simplified three-dimensional displacement discontinuity method was implemented to more accurately calculate fracture displacements and fracture-induced dynamic stress changes than our previously developed pseudo-3d model. The effects of perforation cluster spacing, differential stress (SHmax - Shmin) and various geometry natural fracture patterns on injection pressure and fracture complexity were investigated. The single stage simulation results shown that (1) higher differential stress suppresses fracture length and increases injection pressure; (2) there is an optimal choice for the number of fractures per stage to maximize effective fracture surface area, beyond which increasing the number of fractures actually decreases effective fracture area; and (3) fracture complexity is a function of natural fracture patterns (various regular pattern geometries were investigated). Natural fractures with small relative angle to hydraulic fractures are more likely to control fracture propagation path. Also, natural fracture patterns with more long fractures tend to increase the likelihood to dominate the preferential fracture trend of fracture trajectory. Our numerical model can provide a physics-based complex fracture network that can be imported into reservoir simulation models for production analysis. The overall sensitivity results presented should serve as guidelines for fracture complexity analysis.


Sign in / Sign up

Export Citation Format

Share Document