Impact of Complex Fracture Networks on Rate Transient Behavior of Wells in Unconventional Reservoirs Based on Embedded Discrete Fracture Model

2021 ◽  
pp. 1-12
Author(s):  
Jiazheng Qin ◽  
Yingjie Xu ◽  
Yong Tang ◽  
Rui Liang ◽  
Qianhu Zhong ◽  
...  

Abstract It has recently been demonstrated that complex fracture networks (CFN) especially activated natural fractures (ANF) play an important role in unconventional reservoir development. However, traditional rate transient analysis (RTA) methods barely investigate the impact of CFN or ANF. Furthermore, the influence of CFN on flow regime is still ambiguous. Failure to consider these effects could lead to misdiagnosis of flow regimes and underestimation of original oil in place (OOIP). A novel numerical RTA method is therefore presented herein to improve the quality of reserves assessment. A new methodology is introduced. Propagating hydraulic fractures (HF) can generate different stress perturbations to allow natural fractures (NF) to fail, forming various ANF pattern. An embedded discrete fracture model (EDFM) of ANF is stochastically generated instead of local grid refinement (LGR) method to overcome the time-intensive computation time. These models are coupled with reservoir models using non-neighboring connections (NNCs). Results show that except for simplified models used in previous studies subjected to traditional concept of stimulated reservoir volume (SRV), in our study, the ANF region has been discussed to emphasis the impact of NF on simulation results. Henceforth, ANF could be only concentrated around the near-wellbore region, and it may also cover the whole simulation area. Obvious distinctions could be viewed for different kinds of ANF on diagnostic plots. Instead of SRV-dominated flow mentioned in previous studies, ANF-dominated flow developed in this work is shown to be more reasonable. Also, new flow regimes such as interference flow inside and outside activated natural fracture flow region (ANFR) are found. In summary, better evaluation of reservoir properties and reserves assessment such as OOIP are achieved based on our proposed model compared with conventional models. The novel RTA method considering CFN presented herein is an easy-to-apply numerical RTA technique that can be applied for reservoir and fracture characterization as well as OOIP assessment.

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.


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.


SPE Journal ◽  
2020 ◽  
Vol 25 (04) ◽  
pp. 2143-2161 ◽  
Author(s):  
Olufemi Olorode ◽  
Bin Wang ◽  
Harun Ur Rashid

Summary Most unconventional oil and gas reservoirs are known to have several natural fractures in different orientations, which are consistent with the prevailing stresses when they were created. The accurate and efficient modeling of natural and hydraulic fractures presents a significant computational challenge. In this work, we show the limitations of the embedded discrete-fracture model (EDFM) and present the first 3D projection-based EDFM (pEDFM) algorithm and compositional simulation studies with realistic fracture networks in a fully 3D space. The simulation results from this work indicate that the pEDFM presented can model realistic fractured unconventional reservoirs accurately and efficiently. To validate the model, we present some simplistic fracture cases that can be meshed and modeled easily using explicit-fracture modeling in commercial-reservoir simulators. From the cases studied, we observe that using progressively finer grids near the hydraulic-fracture surfaces helps to improve model accuracy because this allows us to capture the sharp pressure drops expected near these fracture surfaces. The simulation results show that, unlike EDFM, the robust pEDFM algorithm presented here is accurate even at the low fracture-conductivity values expected in many of these ubiquitous natural fractures. In this paper, we present the first full 3D compositional modeling with pEDFM. We demonstrate that our model can accurately and efficiently model multiply fractured horizontal wells in unconventional reservoirs, which have complex networks of thousands of fractures at various orientations.


2021 ◽  
Author(s):  
Wei Yu ◽  
Anuj Gupta ◽  
Ravimadhav N. Vaidya ◽  
Kamy Sepehrnoori

Abstract The complexity of dynamic modeling for naturally fractured reservoirs has increased in recent years to incorporate more data and physics, as well as to handle advanced completion designs and development scenarios. While these complex models can provide more insight to difficult problems, they come with higher computational costs. Such a limitation prohibits an asset team from working with a large number of well/fracture scenarios that correctly represent geological uncertainty. This study presents a powerful non-intrusive Embedded Discrete Fracture Model (EDFM) method to efficiently handle millions of natural and hydraulic fractures with hundreds of horizontal wells, which has never been modeled in the literature. Specifically, we built a 3D geological model using a black oil reservoir simulator with 100 square miles in the horizontal area and 11 layers of 165 ft thickness. The total number of matrix cells without considering fractures is over 3 million. In total, 400 horizontal wells with well length of 6000 ft were modeled in two target layers. Each layer contains 200 wells. Each well has 112 hydraulic fractures with cluster spacing of 50 ft. The total number of hydraulic fractures is 44,800. In addition, we generated three cases with 10K, 100K and 1 million 3D natural fractures with dip angle from 70 to 90 degrees. For the case with 1 million natural fractures, the total number of cells is over 42 million. Well performance for the field example, with and without natural fractures, was investigated. This work adds significant value to the well and fracture spacing optimization process during field development planning. The non-intrusive EDFM method has been proven to be an efficient fracture modeling tool for simulating million-level complex hydraulic/natural fractures, which significantly improves accuracy and reduces computational time.


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.


2015 ◽  
Author(s):  
Mark W. McClure ◽  
Mohsen Babazadeh ◽  
Sogo Shiozawa ◽  
Jian Huang

Abstract We developed a hydraulic fracturing simulator that implicitly couples fluid flow with the stresses induced by fracture deformation in large, complex, three-dimensional discrete fracture networks. The simulator can describe propagation of hydraulic fractures and opening and shear stimulation of natural fractures. Fracture elements can open or slide, depending on their stress state, fluid pressure, and mechanical properties. Fracture sliding occurs in the direction of maximum resolved shear stress. Nonlinear empirical relations are used to relate normal stress, fracture opening, and fracture sliding to fracture aperture and transmissivity. Fluid leakoff is treated with a semianalytical one-dimensional leakoff model that accounts for changing pressure in the fracture over time. Fracture propagation is treated with linear elastic fracture mechanics. Non-Darcy pressure drop in the fractures due to high flow rate is simulated using Forchheimer's equation. A crossing criterion is implemented that predicts whether propagating hydraulic fractures will cross natural fractures or terminate against them, depending on orientation and stress anisotropy. Height containment of propagating hydraulic fractures between bedding layers can be modeled with a vertically heterogeneous stress field or by explicitly imposing hydraulic fracture height containment as a model assumption. The code is efficient enough to perform field-scale simulations of hydraulic fracturing with a discrete fracture network containing thousands of fractures, using only a single compute node. Limitations of the model are that all fractures must be vertical, the mechanical calculations assume a linearly elastic and homogeneous medium, proppant transport is not included, and the locations of potentially forming hydraulic fractures must be specified in advance. Simulations were performed of a single propagating hydraulic fracture with and without leakoff to validate the code against classical analytical solutions. Field-scale simulations were performed of hydraulic fracturing in a densely naturally fractured formation. The simulations demonstrate how interaction with natural fractures in the formation can help explain the high net pressures, relatively short fracture lengths, and broad regions of microseismicity that are often observed in the field during stimulation in low permeability formations, and which are not predicted by classical hydraulic fracturing models. Depending on input parameters, our simulations predicted a variety of stimulation behaviors, from long hydraulic fractures with minimal leakoff into surrounding fractures to broad regions of dense fracturing with a branching network of many natural and newly formed fractures.


SPE Journal ◽  
2016 ◽  
Vol 21 (04) ◽  
pp. 1302-1320 ◽  
Author(s):  
Mark W. McClure ◽  
Mohsen Babazadeh ◽  
Sogo Shiozawa ◽  
Jian Huang

Summary We developed a hydraulic-fracturing simulator that implicitly couples fluid flow with the stresses induced by fracture deformation in large, complex, 3D discrete-fracture networks (DFNs). The code is efficient enough to perform field-scale simulations of hydraulic fracturing in DFNs containing thousands of fractures, without relying on distributed-memory parallelization. The simulator can describe propagation of hydraulic fractures and opening and shear stimulation of natural fractures. Fracture elements can open or slide, depending on their stress state, fluid pressure, and mechanical properties. Fracture sliding occurs in the direction of maximum resolved shear stress. Nonlinear empirical equations are used to relate normal stress, fracture opening, and fracture sliding to fracture aperture and transmissivity. Fluid leakoff is treated with a semianalytical 1D leakoff model that accounts for changing pressure in the fracture over time. Fracture propagation is modeled with linear-elastic fracture mechanics. The Forchheimer equation (Forchheimer 1901) is used to simulate non-Darcy pressure drop in the fractures because of high flow rate. A crossing criterion is implemented that predicts whether propagating hydraulic fractures will cross natural fractures or terminate against them, depending on orientation and stress anisotropy. Height containment of propagating hydraulic fractures between bedding layers can be modeled with a vertically heterogeneous stress field or by explicitly imposing hydraulic-fracture-height containment as a model assumption. Limitations of the model are that all fractures must be vertical; the mechanical calculations assume a linearly elastic and homogeneous medium; proppant transport is not included; and the locations of potentially forming hydraulic fractures must be specified in advance. Simulations were performed of a single propagating hydraulic fracture with and without leakoff to validate the code against classical analytical solutions. Field-scale simulations were performed of hydraulic fracturing in a densely naturally fractured formation. The simulations demonstrate how interaction with natural fractures in the formation can help explain the high net pressures, relatively short fracture lengths, and broad regions of microseismicity that are often observed in the field during stimulation in low-permeability formations, and that are not predicted by classical hydraulic-fracturing models. Depending on input parameters, our simulations predicted a variety of stimulation behaviors, from long hydraulic fractures with minimal leakoff into surrounding fractures to broad regions of dense fracturing with a branching network of many natural and newly formed fractures.


Sign in / Sign up

Export Citation Format

Share Document