Impact of Hydraulic Fracture Fairway Development in Multi-Stage Horizontal Laterals - A Production Flow Simulation Study

2021 ◽  
Author(s):  
Abdul Muqtadir Khan ◽  
Jon Olson

Abstract The vast shale gas and tight oil reservoirs cannot be economically developed without multi-stage hydraulic fracture treatments. Owing to the disparity in the density of natural fractures and the different in-situ stress conditions in these formations, micro-seismic fracture mapping has shown that hydraulic fracture treatments develop a range of large-scale fracture networks. The effect of these various fracture geometries on production is a subject matter in question. The fracture networks approximated with micro-seismic mapping are integrated with a commercial numerical production simulator that discretely models different network structures. Two fracture geometries have been broadly proposed, i.e., orthogonal and transverse. The orthogonal pattern represents a network with cross-cutting fractures orthogonal to each other, whereas transverse profile maps uninterrupted fractures achieving maximum depth of penetration into the reservoir. The response for a single stage is further investigated by comparing the propagation of each stage to be dendritic versus planar. A dendritic propagation is a bifurcation of the induced hydraulic fracture due to the intersection with the natural fracture (failure along the plane of weakness). For the same injected fracture treatment volume, the transverse network attains a higher penetration into the reservoir, achieves a higher stimulated reservoir volume (SRV), and produces around 40-65% more than the orthogonal network over a timespan of 10 years. The SRV will largely dictate the drainage area in a tight environment. The cumulative production rises until the pressure drawdown reaches the extent of the fracture fairway. For the orthogonal network, the unstimulated reservoir boundary is reached at a sooner time than the transverse network. It is found that by increasing the fracture spacing in both the networks from 100 ft to 150 ft, the relative production was enhanced in the orthogonal network by 41%, but when it was further increased to 200 ft- there was a minor drop (not increase) in the relative production (4.5%). For an infinite conductivity fracture, the width of the fracture has minimal effect on oil and gas production. For the dendritic pattern, the size of the SRV created due to the interaction between the induced and natural fractures largely depends on the length of natural fractures and the point of interaction (center, off-center, or extremity). Effect of length, distance of natural fracture from wellbore, and the point of interaction is evaluated. A novel approach for reservoir simulation is used, where porosity (instead of permeability) is used as a scaling parameter for the fracture width. The forward modeling effort, including the comparative fracture geometries setup, induced, and natural fracture interaction parametric study, is unique.

2016 ◽  
Vol 56 (1) ◽  
pp. 225 ◽  
Author(s):  
Kunakorn Pokalai ◽  
David Kulikowski ◽  
Raymond L. Johnson ◽  
Manouchehr Haghighi ◽  
Dennis Cooke

Hydraulic fracturing in tight gas reservoirs has been performed in the Cooper Basin for decades in reservoirs containing high stress and pre-existing natural fractures, especially near faults. The hydraulic fracture is affected by factors such as tortuosity, high entry pressures, and the rock fabric including natural fractures. These factors cause fracture plane rotation and complexities, leading to fracture disconnection or reduced proppant placement during the treatment. In this paper, rock properties are estimated for a targeted formation using well logs to create a geomechanical model. Natural fracture and stress azimuths within the interval were interpreted from borehole image logs. The image log interpretations inferred that fissures are oriented 30–60° relative to the maximum horizontal stress. Next, diagnostic fracture injection test (DFIT) data was used with the poro-elastic stress equations to predict tectonic strains. Finally, the geomechanical model was history-matched with a planar 3D hydraulic fracturing simulator, and gave more insight into fracture propagation in an environment of pre-existing natural fractures. The natural fracture azimuths and calibrated geomechanical model are input into a framework to evaluate varying scenarios that might result based on a vertical or inclined well design. A well design is proposed based on the natural fracture orientation relative to the hydraulic fracture that minimises complexity to optimise proppant placement. In addition, further models and diagnostics are proposed to aid predicting the hydraulically induced fracture geometry, its impact on gas production, and optimising wellbore trajectory to positively interact with pre-existing natural fractures.


2015 ◽  
Author(s):  
Hisanao Ouchi ◽  
Amit Katiyar ◽  
John T. Foster ◽  
Mukul M. Sharma

Abstract A novel fully coupled hydraulic fracturing model based on a nonlocal continuum theory of peridynamics is presented and applied to the fracture propagation problem. It is shown that this modeling approach provides an alternative to finite element and finite volume methods for solving poroelastic and fracture propagation problems and offers some clear advantages. In this paper we specifically investigate the interaction between a hydraulic fracture and natural fractures. Current hydraulic fracturing models remain limited in their ability to simulate the formation of non-planar, complex fracture networks. The peridynamics model presented here overcomes most of the limitations of existing models and provides a novel approach to simulate and understand the interaction between hydraulic fractures and natural fractures. The model predictions in two-dimensions have been validated by reproducing published experimental results where the interaction between a hydraulic fracture and a natural fracture is controlled by the principal stress contrast and the approach angle. A detailed parametric study involving poroelasticity and mechanical properties of the rock is performed to understand why a hydraulic fracture gets arrested or crosses a natural fracture. This analysis reveals that the poroelasticity, resulting from high fracture fluid leak-off, has a dominant influence on the interaction between a hydraulic fracture and a natural fracture. In addition, the fracture toughness of the rock, the toughness of the natural fracture, and the shear strength of the natural fracture also affect the interaction between a hydraulic fracture and a natural fracture. Finally, we investigate the interaction of multiple completing fractures with natural fractures in two-dimensions and demonstrate the applicability of the approach to simulate complex fracture networks on a field scale.


Processes ◽  
2018 ◽  
Vol 6 (8) ◽  
pp. 113 ◽  
Author(s):  
Shen Wang ◽  
Huamin Li ◽  
Dongyin Li

To investigate the mechanism of hydraulic fracture propagation in coal seams with discontinuous natural fractures, an innovative finite element meshing scheme for modeling hydraulic fracturing was proposed. Hydraulic fracture propagation and interaction with discontinuous natural fracture networks in coal seams were modeled based on the cohesive element method. The hydraulic fracture network characteristics, the growth process of the secondary hydraulic fractures, the pore pressure distribution and the variation of bottomhole pressure were analyzed. The improved cohesive element method, which considers the leak-off and seepage behaviors of fracturing liquid, is capable of modeling hydraulic fracturing in naturally fractured formations. The results indicate that under high stress difference conditions, the hydraulic fracture network is spindle-shaped, and shows a multi-level branch structure. The ratio of secondary fracture total length to main fracture total length was 2.11~3.62, suggesting that the secondary fractures are an important part of the hydraulic fracture network in coal seams. In deep coal seams, the break pressure of discontinuous natural fractures mainly depends on the in-situ stress field and the direction of natural fractures. The mechanism of hydraulic fracture propagation in deep coal seams is significantly different from that in hard and tight rock layers.


Author(s):  
Yunsuk Hwang ◽  
Jiajing Lin ◽  
David Schechter ◽  
Ding Zhu

Multiple hydraulic fracture treatments in reservoirs with natural fractures create complex fracture networks. Predicting well performance in such a complex fracture network system is an extreme challenge. The statistical nature of natural fracture networks changes the flow characteristics from that of a single linear fracture. Simply using single linear fracture models for individual fractures, and then summing the flow from each fracture as the total flow rate for the network could introduce significant error. In this paper we present a semi-analytical model by a source method to estimate well performance in a complex fracture network system. The method simulates complex fracture systems in a more reasonable approach. The natural fracture system we used is fractal discrete fracture network model. We then added multiple dominating hydraulic fractures to the natural fracture system. Each of the hydraulic fractures is connected to the horizontal wellbore, and some of the natural fractures are connected to the hydraulic fractures through the network description. Each fracture, natural or hydraulically induced, is treated as a series of slab sources. The analytical solution of superposed slab sources provides the base of the approach, and the overall flow from each fracture and the effect between the fractures are modeled by applying the superposition principle to all of the fractures. The fluid inside the natural fractures flows into the hydraulic fractures, and the fluid of the hydraulic fracture from both the reservoir and the natural fractures flows to the wellbore. This paper also shows that non-Darcy flow effects have an impact on the performance of fractured horizontal wells. In hydraulic fracture calculation, non-Darcy flow can be treated as the reduction of permeability in the fracture to a considerably smaller effective permeability. The reduction is about 2% to 20%, due to non-Darcy flow that can result in a low rate. The semi-analytical solution presented can be used to efficiently calculate the flow rate of multistage-fractured wells. Examples are used to illustrate the application of the model to evaluate well performance in reservoirs that contain complex fracture networks.


2020 ◽  
Vol 10 (8) ◽  
pp. 3333-3345
Author(s):  
Ali Al-Rubaie ◽  
Hisham Khaled Ben Mahmud

Abstract All reservoirs are fractured to some degree. Depending on the density, dimension, orientation and the cementation of natural fractures and the location where the hydraulic fracturing is done, preexisting natural fractures can impact hydraulic fracture propagation and the associated flow capacity. Understanding the interactions between hydraulic fracture and natural fractures is crucial in estimating fracture complexity, stimulated reservoir volume, drained reservoir volume and completion efficiency. However, because of the presence of natural fractures with diffuse penetration and different orientations, the operation is complicated in naturally fractured gas reservoirs. For this purpose, two numerical methods are proposed for simulating the hydraulic fracture in a naturally fractured gas reservoir. However, what hydraulic fracture looks like in the subsurface, especially in unconventional reservoirs, remain elusive, and many times, field observations contradict our common beliefs. In this study, the hydraulic fracture model is considered in terms of the state of tensions, on the interaction between the hydraulic fracture and the natural fracture (45°), and the effect of length and height of hydraulic fracture developed and how to distribute induced stress around the well. In order to determine the direction in which the hydraulic fracture is formed strikethrough, the finite difference method and the individual element for numerical solution are used and simulated. The results indicate that the optimum hydraulic fracture time was when the hydraulic fracture is able to connect natural fractures with large streams and connected to the well, and there is a fundamental difference between the tensile and shear opening. The analysis indicates that the growing hydraulic fracture, the tensile and shear stresses applied to the natural fracture.


2015 ◽  
Vol 55 (1) ◽  
pp. 351
Author(s):  
Alireza Keshavarz ◽  
Alexander Badalyan ◽  
Raymond Johnson ◽  
Pavel Bedrikovetski

A method is proposed for enhancing the conductivity of micro-fractures and cleats around the hydraulically induced fractures in coal bed methane reservoirs. In this technique, placing ultra-fine proppant particles in natural fractures and cleats around hydraulically induced fractures at leak-off conditions keeps the coal cleats open during water-gas production, and this consequently increases the efficiency of hydraulic fracturing treatment. Experimental and mathematical studies for the stimulation of a natural cleat system around the main hydraulic fracture are conducted. In the experimental part, core flooding tests are performed to inject a flow of suspended particles inside the natural fractures of a coal sample. By placing different particle sizes and evaluating the concentration of placed particles, an experimental coefficient is found for optimum proppant placement in which the maximum permeability is achieved after proppant placement. In the mathematical modelling study, a laboratory-based mathematical model for graded proppant placement in naturally fractured rocks around a hydraulically induced fracture is proposed. Derivations of the model include an exponential form of the pressure-permeability dependence and accounts for permeability variation in the non-stimulated zone. The explicit formulae are derived for the well productivity index by including the experimentally found coefficient. Particle placement tests resulted in an almost three-times increase in coal permeability. The laboratory-based mathematical modelling, as performed for the field conditions, shows that the proposed method yields around a six-times increase in the productivity index.


2020 ◽  
pp. 014459872096083
Author(s):  
Yulong Liu ◽  
Dazhen Tang ◽  
Hao Xu ◽  
Wei Hou ◽  
Xia Yan

Macrolithotypes control the pore-fracture distribution heterogeneity in coal, which impacts stimulation via hydrofracturing and coalbed methane (CBM) production in the reservoir. Here, the hydraulic fracture was evaluated using the microseismic signal behavior for each macrolithotype with microfracture imaging technology, and the impact of the macrolithotype on hydraulic fracture initiation and propagation was investigated systematically. The result showed that the propagation types of hydraulic fractures are controlled by the macrolithotype. Due to the well-developed natural fracture network, the fracture in the bright coal is more likely to form the “complex fracture network”, and the “simple” case often happens in the dull coal. The hydraulic fracture differences are likely to impact the permeability pathways and the well productivity appears to vary when developing different coal macrolithtypes. Thus, considering the difference of hydraulic fracture and permeability, the CBM productivity characteristics controlled by coal petrology were simulated by numerical simulation software, and the rationality of well pattern optimization factors for each coal macrolithotype was demonstrated. The results showed the square well pattern is more suitable for dull coal and semi-dull coal with undeveloped natural fractures, while diamond and rectangular well pattern is more suitable for semi-bright coal and bright coal with more developed natural fractures and more complex fracturing fracture network; the optimum wells spacing of bright coal and semi-bright coal is 300 m and 250 m, while that of semi-dull coal and dull coal is just 200 m.


2020 ◽  
Author(s):  
Simon Oldfield ◽  
Mikael Lüthje ◽  
Michael Welch ◽  
Florian Smit

<p>Large scale modelling of fractured reservoirs is a persistent problem in representing fluid flow in the subsurface. Considering a geothermal energy prospect beneath the Drenthe Aa area, we demonstrate application of a recently developed approach to efficiently predict fracture network geometry across an area of several square kilometres.</p><p>Using a strain based method to mechanically model fracture nucleation and propagation, we generate a discretely modelled fracture network consisting of individual failure planes, opening parallel and perpendicular to the orientation of maximum and minimum strain. Fracture orientation, length and interactions vary following expected trends, forming a connected fracture network featuring population statistics and size distributions comparable to outcrop examples.</p><p>Modelled fracture networks appear visually similar to natural fracture networks with spatial variation in fracture clustering and the dominance of major and minor fracture trends.</p><p>Using a network topology approach, we demonstrate that the predicted fracture network shares greater geometric similarity with natural networks. Considering fluid flow through the model, we demonstrate that hydraulic conductivity and flow anisotropy are strongly dependent on the geometric connection of fracture sets.</p><p>Modelling fracture evolution mechanically allows improved representation of geometric aspects of fracture networks to which fluid flow is particularly sensitive. This method enables rapid generation of discretely modelled fractures over large areas and extraction of suitable summary statistics for reservoir simulation. Visual similarity of the output models improves our ability to compare between our model and natural analogues to consider model validation.</p>


Author(s):  
Zhaozhong Yang ◽  
Rui He ◽  
Xiaogang Li ◽  
Zhanling Li ◽  
Ziyuan Liu

The tight sandstone gas reservoir in southern Songliao Basin is naturally fractured and is characterized by its low porosity and permeability. Large-scale hydraulic fracturing is the most effective way to develop this tight gas reservoir. Quantitative evaluation of fracability is essential for optimizing a fracturing reservoir. In this study, as many as ten fracability-related factors, particularly mechanical brittleness, mineral brittleness, cohesion, internal friction angle, unconfined compressive strength (UCS), natural fracture, Model-I toughness, Model-II toughness, horizontal stress difference, and fracture barrier were obtained from a series of petrophysical and geomechanical experiments are analyzed. Taking these influencing factors into consideration, a modified comprehensive evaluation model is proposed based on the analytic hierarchy process (AHP). Both a transfer matrix and a fuzzy matrix were introduced into this model. The fracability evaluation of four reservoir intervals in Jinshan gas field was analyzed. Field fracturing tests were conducted to verify the efficiency and accuracy of the proposed evaluation model. Results showed that gas production is higher and more stable in the reservoir interval with better fracability. The field test data coincides with the results of the proposed evaluation model.


SPE Journal ◽  
2018 ◽  
Vol 24 (01) ◽  
pp. 302-318 ◽  
Author(s):  
Jixiang Huang ◽  
Joseph P. Morris ◽  
Pengcheng Fu ◽  
Randolph R. Settgast ◽  
Christopher S. Sherman ◽  
...  

Summary A fully coupled finite-element/finite-volume code is used to model 3D hydraulically driven fractures under the influence of strong vertical variations in closure stress interacting with natural fractures. Previously unknown 3D interaction mechanisms on fracture-height growth are revealed. Slipping of a natural fracture, triggered by elevated fluid pressure from an intersecting hydraulic fracture, can induce both increases and decreases of normal stress in the minimum-horizontal-stress direction, toward the center and tip of the natural fracture, respectively. Consequently, natural fractures are expected to be able to both encourage and inhibit the progress of hydraulic fractures propagating through stress barriers, depending on the relative locations between the intersecting fractures. Once the hydraulic fracture propagates above the stress barrier through the weakened segment near a favorably located natural fracture, a configuration consisting of two opposing fractures cuts the stress barrier from above and below. The fluid pressure required to break the stress barrier under such opposing-fracture configurations is substantially lower than that required by a fracture penetrating the same barrier from one side. Sensitivity studies of geologic conditions and operational parameters have also been performed to explore the feasibility of controlled fracture height. The interactions between hydraulic fractures, natural fractures, and geologic factors such as stress barriers in three dimensions are shown to be much more complex than in two dimensions. Although it is impossible to exhaust all the possible configurations, the ability of a 3D, fully coupled numerical model to naturally capture these processes is well-demonstrated.


Sign in / Sign up

Export Citation Format

Share Document