Coupling Production Data with the DFN to Unravel Natural Fracture Networks in a Tight Naturally Fractured Jurassic Reservoir, Kuwait

2021 ◽  
Author(s):  
Raphael Altman ◽  
Mariela Pichardi ◽  
Pratik Sangani ◽  
Tahani Al Rashidi ◽  
Girija Shankar Padhy ◽  
...  

Abstract The Jurassic Najmah-Sargelu of west Kuwait can be thought of as a "hybrid" between a conventional and an unconventional reservoir. These systems form an increasingly important resource for operators, but their performance is unpredictable because matrix permeability is in the micro-Darcy range and production depends on natural fractures. Success depends on how well the static models are aligned to the dynamic production, and the effectiveness of a fit-for-purpose multistage completion on project economics. In this work we present our lessons learnt in production modelling these reservoirs and the coupling between reservoir simulation and the discrete fracture network (DFN). Our reservoir models were constructed using a highly integrated approach incorporating data from all scales and disciplines (drilling, geophysical, geological, reservoir and production) and the production simulations were run using dual porosity and black oil models. As expected, the DFN played a key part of this effort. An iterative approach was used to adjust the DFN so that it was consistent with production observations. However, in all cases care was made to ensure the new DFN honoured the seismic, geological, well log and drilling data from which it was generated. Final, smaller adjustments were made to the simulation model at the log scale to match PLT data. We used uncertainty analysis to run hundreds of simulation cases and found that the character of the natural fractures is quite well imprinted in the observed production data, particularly pressure buildup data. This gave us a better understanding of whether the natural fractures are diffuse and laterally extensive away from the wellbore or if they are localized close to the wellbore. Where reservoir simulation history matches inferred laterally extensive natural fractures, an good correlation was obtained with the natural fracturing from the DFN. This correlation was poor where natural fracturing was confined to a smaller depth interval (as observed from PLT), and is a result of the limitation in seismic resolution to resolve these natural fractures. The lessons learnt from our work helps towards improved understanding of production mechanisms of these reservoirs and their natural fracture networks. This, together with higher resolution azimuthal seismic, advanced wellbore characterization data and multistage completions are the desired key ingredients for technically enhancing production in these reservoirs.

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.


2014 ◽  
Vol 852 ◽  
pp. 831-834
Author(s):  
Bin Dou ◽  
Hui Gao ◽  
Gang Zhou ◽  
Lei Ren

This paper presents an advanced computational modeling of natural fracture networks in hot dry rock (HDR) reservoirs. A mathematical model is developed for describing the heat energy extracted from a HDR in a multi-well system. The model stochastically simulates discrete properties of natural fractures, utilizing multi-set orientation and fractal mathematics. The simulated fracture networks are essential for further stimulation and fluid flow studies. The results show that the heat extraction effectiveness is affected significantly by the well spacing, well radius, reservoir thickness, and pumped flow rate in a multi-well system. The water temperature decreases with increasing pumping rate and increases with the well spacing, well radius, and reservoir thickness. This paper also examines the detrimental effects of the simulated natural fracture network on the stimulated fluid flow capacity. The effective permeability enhancement (due to hydraulic stimulation) is found almost proportional to the density of the reservoir natural fractures.


Author(s):  
Hannes Hofmann ◽  
Tayfun Babadagli ◽  
Günter Zimmermann

The creation of large complex fracture networks by hydraulic fracturing is imperative for enhanced oil recovery from tight sand or shale reservoirs, tight gas extraction, and Hot-Dry-Rock (HDR) geothermal systems to improve the contact area to the rock matrix. Although conventional fracturing treatments may result in bi-wing fractures, there is evidence by microseismic mapping that fracture networks can develop in many unconventional reservoirs, especially when natural fracture systems are present and the differences between the principle stresses are low. However, not much insight is gained about fracture development as well as fluid and proppant transport in naturally fractured tight formations. In order to clarify the relationship between rock and treatment parameters, and resulting fracture properties, numerical simulations were performed using a commercial Discrete Fracture Network (DFN) simulator. A comprehensive sensitivity analysis is presented to identify typical fracture network patterns resulting from massive water fracturing treatments in different geological conditions. It is shown how the treatment parameters influence the fracture development and what type of fracture patterns may result from different treatment designs. The focus of this study is on complex fracture network development in different natural fracture systems. Additionally, the applicability of the DFN simulator for modeling shale gas stimulation and HDR stimulation is critically discussed. The approach stated above gives an insight into the relationships between rock properties (specifically matrix properties and characteristics of natural fracture systems) and the properties of developed fracture networks. Various simulated scenarios show typical conditions under which different complex fracture patterns can develop and prescribe efficient treatment designs to generate these fracture systems. Hydraulic stimulation is essential for the production of oil, gas, or heat from ultratight formations like shales and basement rocks (mainly granite). If natural fracture systems are present, the fracturing process becomes more complex to simulate. Our simulation results reveal valuable information about main parameters influencing fracture network properties, major factors leading to complex fracture network development, and differences between HDR and shale gas/oil shale stimulations.


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>


GeoArabia ◽  
2001 ◽  
Vol 6 (1) ◽  
pp. 27-42
Author(s):  
Stephen J. Bourne ◽  
Lex Rijkels ◽  
Ben J. Stephenson ◽  
Emanuel J.M. Willemse

ABSTRACT To optimise recovery in naturally fractured reservoirs, the field-scale distribution of fracture properties must be understood and quantified. We present a method to systematically predict the spatial distribution of natural fractures related to faulting and their effect on flow simulations. This approach yields field-scale models for the geometry and permeability of connected fracture networks. These are calibrated by geological, well test and field production data to constrain the distributions of fractures within the inter-well space. First, we calculate the stress distribution at the time of fracturing using the present-day structural reservoir geometry. This calculation is based on a geomechanical model of rock deformation that represents faults as frictionless surfaces within an isotropic homogeneous linear elastic medium. Second, the calculated stress field is used to govern the simulated growth of fracture networks. Finally, the fractures are upscaled dynamically by simulating flow through the discrete fracture network per grid block, enabling field-scale multi-phase reservoir simulation. Uncertainties associated with these predictions are considerably reduced as the model is constrained and validated by seismic, borehole, well test and production data. This approach is able to predict physically and geologically realistic fracture networks. Its successful application to outcrops and reservoirs demonstrates that there is a high degree of predictability in the properties of natural fracture networks. In cases of limited data, field-wide heterogeneity in fracture permeability can be modelled without the need for field-wide well coverage.


Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-11
Author(s):  
Gou Feifei ◽  
Liu Chuanxi ◽  
Ren Zongxiao ◽  
Qu Zhan ◽  
Wang Sukai ◽  
...  

Unconventional resources have been successfully exploited with technological advancements in horizontal-drilling and multistage hydraulic-fracturing, especially in North America. Due to preexisting natural fractures and the presence of stress isotropy, several complex fracture networks can be generated during fracturing operations in unconventional reservoirs. Using the DVS method, a semianalytical model was created to analyze the transient pressure behavior of a complex fracture network in which hydraulic and natural fractures interconnect with inclined angles. In this model, the complex fracture network can be divided into a proper number of segments. With this approach, we are able to focus on a detailed description of the network properties, such as the complex geometry and varying conductivity of the fracture. The accuracy of the new model was demonstrated by ECLIPSE. Using this method, we defined six flow patterns: linear flow, fracture interference flow, transitional flow, biradial flow, pseudoradial flow, and boundary response flow. A sensitivity analysis was conducted to analyze each of these flow regimes. This work provides a useful tool for reservoir engineers for fracture designing as well as estimating the performance of a complex fracture network.


2020 ◽  
Vol 54 ◽  
pp. 149-156
Author(s):  
Ajay K. Sahu ◽  
Ankur Roy

Abstract. It is well known that fracture networks display self-similarity in many cases and the connectivity and flow behavior of such networks are influenced by their respective fractal dimensions. In the past, the concept of lacunarity, a parameter that quantifies spatial clustering, has been implemented by one of the authors in order to demonstrate that a set of seven nested natural fracture maps belonging to a single fractal system, but of different visual appearances, have different clustering attributes. Any scale-dependency in the clustering of fractures will also likely have significant implications for flow processes that depend on fracture connectivity. It is therefore important to address the question as to whether the fractal dimension alone serves as a reasonable proxy for the connectivity of a fractal-fracture network and hence, its flow response or, if it is the lacunarity, a measure of scale-dependent clustering, that may be used instead. The present study attempts to address this issue by exploring possible relationships between the fractal dimension, lacunarity and connectivity of fractal-fracture networks. It also endeavors to study the relationship between lacunarity and fluid flow in such fractal-fracture networks. A set of deterministic fractal-fracture models generated at different iterations and, that have the same theoretical fractal dimension are used for this purpose. The results indicate that such deterministic synthetic fractal-fracture networks with the same theoretical fractal dimension have differences in their connectivity and that the latter is fairly correlated with lacunarity. Additionally, the flow simulation results imply that lacunarity influences flow patterns in fracture networks. Therefore, it may be concluded that at least in synthetic fractal-fracture networks, rather than fractal dimension, it is the lacunarity or scale-dependent clustering attribute that controls the connectivity and hence the flow behavior.


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.


Sign in / Sign up

Export Citation Format

Share Document