Topological analysis of fracture networks integrated with flow simulation models for equivalent fracture permeability estimation

2021 ◽  
Vol 147 ◽  
pp. 104338
Author(s):  
Jefferson Pedro Silva ◽  
Igor Fernandes Gomes ◽  
Rafael Fernandes Vieira Correia Santos ◽  
Tiago Siqueira de Miranda ◽  
Ricardo Pereira Guedes ◽  
...  
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.


2006 ◽  
Vol 9 (05) ◽  
pp. 513-520 ◽  
Author(s):  
Yu Ding ◽  
Remy Basquet ◽  
Bernard Bourbiaux

Summary One difficulty in fracture upscaling for field-scale dual-porosity reservoir simulation is the determination of equivalent gridblock fracture permeability, which depends on the type of boundary conditions imposed on the discrete-fracture-network (DFN) simulation. Actually, classical upscaling procedures usually are based on linearly varying pressure boundary conditions, which cannot capture the near-well flow behavior. As a result, the well productivity calculated by a dual-porosity flow simulator can be very different from that calculated on a DFN model. This paper proposes a near-well fracture-upscaling procedure based on the geological DFN model to improve the accuracy of well productivity in fractured-reservoir simulators. This procedure enables us to represent the actual flow through the fractures and the exchanges between matrix and fractures in the well vicinity. On the basis of the computed near-well flow pattern, equivalent fracture transmissibilities as well as numerical well indices are determined and assigned to the gridblocks of the dual-porosity reservoir simulator. The reliability and necessity of using the near-well upscaling procedure are demonstrated by examples. Introduction Advanced characterization methodologies are now able to provide realistic models of geological fracture networks (Cacas et al. 2001). In addition, production logging and transient well tests can be simulated with DFN models to validate the geological fracture-network geometry and calibrate the hydraulic properties of fractures (Sarda et al. 2002). However, because of computational limitations, the complex geological DFN model cannot be used straightforwardly to simulate a multiphase-flow production scenario at field scale (Bourbiaux et al. 2002). For such simulations, a dual-porosity reservoir simulator is typically used. The dual-porosity reservoir model, using large gridblocks to discretize the whole reservoir, is a conceptual representation of the actual geology of the fractured medium. The flow properties of the fracture network are then homogenized on gridblocks through upscaling procedures. The upscaling of fracture properties is the problem of translating the geological and hydraulic description of fracture networks into reservoir-simulation parameters. The dual-porosity model requires the determination of equivalent fracture permeability and equivalent matrix-block dimensions or shape factors (Bourbiaux et al. 1997; Sarda et al. 1997). This paper discusses methodologies for upscaling the permeability of a fracture network, especially in the vicinity of the well. Upscaling of fracture permeability has been studied extensively. The commonly used method is numerical, based on flow simulation on a model of the actual fracture network with specific boundary conditions to compute an equivalent gridblock permeability (Sarda et al. 1997). Other methods were also developed; for example, Oda (1985) proposed an analytical equation to calculate the fracture-permeability tensor, and Lough et al. (1997) presented an approach using the boundary-element method, which integrates the contribution of matrix in the equivalent permeability of the fractured medium. When using a numerical approach to determine the equivalent permeability of a fracture network, the upscaled result depends on the type of boundary conditions imposed in the flow simulation. Actually, classical upscaling procedures are usually based on flow simulation in a parallelepipedic model with linear-type pressure boundary conditions, which cannot capture the near-well flow behavior. As a result, the well productivity calculated by a dual-porosity flow simulator can be very different from that calculated on a near-wellbore DFN model.


2005 ◽  
Vol 8 (04) ◽  
pp. 300-309 ◽  
Author(s):  
Zeno G. Philip ◽  
James W. Jennings ◽  
Jon E. Olson ◽  
Stephen E. Laubach ◽  
Jon Holder

Summary In conventional reservoir simulations, gridblock permeabilities are frequently assigned values larger than those observed in core measurements to obtain reasonable history matches. Even then, accuracy with regard to some aspects of the performance such as water or gas cuts, breakthrough times, and sweep efficiencies may be inadequate. In some cases, this could be caused by the presence of substantial flow through natural fractures unaccounted for in the simulation. In this paper, we present a numerical investigation into the effects of coupled fracture-matrix fluid flow on equivalent permeability. A fracture-mechanics-based crack-growth simulator, rather than a purely stochastic method, was used to generate fracture networks with realistic clustering, spacing, and fracture lengths dependent on Young's modulus, the subcritical crack index, the bed thickness, and the tectonic strain. Coupled fracture-matrix fluid-flow simulations of the resulting fracture patterns were performed with a finite-difference simulator to obtain equivalent permeabilities that can be used in a coarse-scale flow simulation. The effects of diagenetic cements completely filling smaller aperture fractures and partially filling larger aperture fractures were also studied. Fractures were represented in finite-difference simulations both explicitly as grid cells and implicitly using nonneighbor connections (NNCs) between grid cells. The results indicate that even though fracture permeability is highly sensitive to fracture aperture, the computed equivalent permeabilities are more sensitive to fracture patterns and connectivity. Introduction High-permeability fracture networks in a matrix system can create high-conductivity channels for the flow of fluids through a reservoir, producing larger flow rates and, therefore, larger apparent permeabilities. The presence of fractures can also cause early breakthrough of the displacing fluid and lead to poorer sweep efficiencies in displacement processes. A better understanding of reservoir performance in such cases may be obtained by including the details of the fluid flow in fractures in a coupled fracture-matrix reservoir flow model. It is very difficult to directly measure interwell fracture-network geometry in sufficient detail to model its effect on reservoir behavior. Thus, most modeling approaches have been statistical, using data from outcrop and wellbore observations to determine distributions of fracture attributes such as fracture length, spacing, and aperture to randomly populate a field. In this paper, we use a mechanistic approach to generate the fracture patterns. Attributes of the fracture network depend on the applied boundary conditions and material properties.


2008 ◽  
Vol 11 (06) ◽  
pp. 1071-1081 ◽  
Author(s):  
Amy Whitaker ◽  
C. Shah Kabir ◽  
Wayne Narr

Summary The extent to which fractures affect fluid pathways is a vital component of understanding and modeling fluid flow in any reservoir. We examined the Wafra Ratawi grainstone for which production extending for 50 years, including recent horizontal drilling, has provided some clues about fractures, but their exact locations, intensity, and overall effect have been elusive. In this study, we find that a limited number of total fractures affect production characteristics of the Ratawi reservoir. Although fractures occur throughout the Wafra field, fracture-influenced reservoir behavior is confined to the periphery of the field where the matrix permeability is low. This work suggests that for the largest part of the field, explicit fractures are not necessary in the next-generation Earth and flow-simulation models. The geologic fracture assessment included seismic fault mapping and fracture interpretation of image logs and cores. Fracture trends are in the northeast and southwest quadrants, and fractures are mineralized toward the south and west of the field. Pressure-falloff tests on some peripheral injectors indicate partial barriers, and most of these wells lie on seismic-scale faults in the reservoir, suggesting partial sealing. A few wells show fractured-reservoir production characteristics, and rate-transient analysis on a few producers indicates localized dual-porosity behavior. Producers proximal to dual-porosity wells display single-porosity behavior, however, to attest to the notion of localized fracture response. The spatially restricted fracture-flow characteristics appear to correlate with fracture or vug zones in a low-permeability reservoir. Presence of fracture-flow behavior was tested by constructing the so-called flow-capacity index (FCI), the ratio of khwell (well test-derived value) to khmatrix (core-derived property). Data from 80 wells showed khmatrix to be consistently higher than khwell, a relationship that suggests insignificant fracture production in these wells. Introduction The Wafra field is in the Partitioned Neutral Zone (PNZ) between Kuwait and Saudi Arabia, as shown in Fig. 1. The field has been producing since the 1950s and has seen renewed drilling activity since the late 1990s, including horizontal drilling and implementation of peripheral water injection (Davis and Habib 1999). The Lower Cretaceous Ratawi formation contains the most reserves of the producing intervals at Wafra. The Ratawi oolite (a misnomer--it is a grainstone) reservoir has variable porosity (5 to 35%) and permeability that ranges from tens to hundreds of md (Longacre and Ginger 1988). The main Wafra structure is a gentle (i.e., interlimb angle >170°), doubly plunging anticline trending north-northwest to south-southeast, which culminates near its northern end. The East Wafra spur is a north-trending branch that extends from the center of the main Wafra structure. As seen in Fig. 1, relief on the Main Wafra structure exceeds that on East Wafra. The Ratawi oolite in the Wafra field has been studied at length, and various authors have reported geologic and engineering elements, leading to reservoir characterization and understanding of reservoir performance. Geologic studies are those of Waite et al. (2000) and Sibley et al. (1997). In contrast, Davis and Habib (1999) presented implementation of peripheral water injection, whereas Chawathé et al. (2006) discussed realignment of injection pattern owing to lack of pressure support in the reservoir interior. Previous studies considered the reservoir to behave like a single-porosity system. But recent image-log fracture interpretations indicate high fracture densities, suggesting that the implementation of a dual-porosity model may be necessary because the high impact of fractures during field development has been recognized in some Middle East reservoirs for more than 50 years (Daniel 1954). Static and dynamic data are required to characterize fracture reservoir behavior accurately (Narr et al. 2006). Geologic description of the fracture system, by use of cores, borehole images, seismic data, and well logs, does not in itself determine whether fractures affect reservoir behavior. While seismic and some image logs were available to locate fractures in the Wafra Ratawi reservoir, no dynamic testing with the specific objective of understanding fracture impact has occurred. So, to determine whether fractures influence oil productivity significantly, we used diagnostic analyses of production data and well tests of available injectors. The assessment of fracture effects in the Ratawi reservoir will be used to guide the next generation of geologic and flow-simulation models. Dynamic data involving pressure and rate have the potential to reveal the influence of open fractures in production performance. Unfortunately, pressure-transient testing on single wells does not always provide conclusive evidence about the presence of fractures with the characteristic dual-porosity dip on the pressure-derivative signature (Bourdet et al. 1989). That is because a correct mixture of matrix/fracture storativity must be present for the characteristic signature to appear (Serra et al. 1983). In practice, interference testing (Beliveau 1989) between wells appears to provide more-definitive clues about interwell connectivity, leading to inference about fractures. In contrast to pressure-transient testing, rate-transient analysis offers the potential to provide the same information without dedicated testing. In this field, all wells are currently on submersible pumps. Consequently, the pump-intake pressure and measured rate provided the necessary data for pressure/rate convolution or rate-transient analysis. We provide the Ratawi-reservoir case study primarily as an example of the integration of diverse geologic and engineering data to develop an assessment of fracture influence on reservoir behavior. It illustrates the use of production-data diagnostic tests to determine fracture influence in the absence of targeted fracture-analysis testing. The workflow can be applied to similar static/dynamic problems, such as fault-transmissivity determination. Secondly, this analysis illustrates the process of deciding that fractures, although present throughout the reservoir, may not lead to widespread fractured-reservoir characteristics (e.g., Allan and Sun 2003).


2021 ◽  
Author(s):  
Mokhles Mezghani ◽  
Mustafa AlIbrahim ◽  
Majdi Baddourah

Abstract Reservoir simulation is a key tool for predicting the dynamic behavior of the reservoir and optimizing its development. Fine scale CPU demanding simulation grids are necessary to improve the accuracy of the simulation results. We propose a hybrid modeling approach to minimize the weight of the full physics model by dynamically building and updating an artificial intelligence (AI) based model. The AI model can be used to quickly mimic the full physics (FP) model. The methodology that we propose consists of starting with running the FP model, an associated AI model is systematically updated using the newly performed FP runs. Once the mismatch between the two models is below a predefined cutoff the FP model is switch off and only the AI model is used. The FP model is switched on at the end of the exercise either to confirm the AI model decision and stop the study or to reject this decision (high mismatch between FP and AI model) and upgrade the AI model. The proposed workflow was applied to a synthetic reservoir model, where the objective is to match the average reservoir pressure. For this study, to better account for reservoir heterogeneity, fine scale simulation grid (approximately 50 million cells) is necessary to improve the accuracy of the reservoir simulation results. Reservoir simulation using FP model and 1024 CPUs requires approximately 14 hours. During this history matching exercise, six parameters have been selected to be part of the optimization loop. Therefore, a Latin Hypercube Sampling (LHS) using seven FP runs is used to initiate the hybrid approach and build the first AI model. During history matching, only the AI model is used. At the convergence of the optimization loop, a final FP model run is performed either to confirm the convergence for the FP model or to re iterate the same approach starting from the LHS around the converged solution. The following AI model will be updated using all the FP simulations done in the study. This approach allows the achievement of the history matching with very acceptable quality match, however with much less computational resources and CPU time. CPU intensive, multimillion-cell simulation models are commonly utilized in reservoir development. Completing a reservoir study in acceptable timeframe is a real challenge for such a situation. The development of new concepts/techniques is a real need to successfully complete a reservoir study. The hybrid approach that we are proposing is showing very promising results to handle such a challenge.


Author(s):  
Shovan Lal Chattoraj ◽  
Prashant K. Champati ray ◽  
Sudhakar Pardeshi ◽  
Vikram Gupta ◽  
Yateesh Ketholia

Abstract. Debris flows, a type of landslides, are not nowadays limited only to the periodic devastation of the geologically fragile Himalaya but also ubiquitous in weathered Deccan Volcanic Province of the cratonic south Indian peninsula. Comprehensive assessment of landslide hazard, pertinently, requires process-based modeling using simulation methods. Development of precipitation triggered debris flow simulation models of real events are still at a young stage in India, albeit, especially in tectonically less disturbed regions. A highly objective simulation technique has therefore been envisaged herein to model the debris flow run-out happened in Malin. This takes cues from a high- resolution DEM and other ancillary ground data including geotechnical and frictional parameters. The algorithm is based on Voellmy frictional (dry and turbulent frictional coefficients, μ and ξ respectively) parameters of debris flow with pre-defined release area identified on high-resolution satellite images like LISS-IV and Cartosat-1. The model provides critical quantitative information on flow 1) Velocity, 2) Height, 3) Momentum, and 4) Pressure along the entrainment path. The simulated velocity of about 16 m/s at mid-way the slide plummeted to 6.2 m/s at the base with intermittently increased and decreased values. The simulated maximum height was 3.9 m which gradually declined to 1.5 m near the bottom. The results can be beneficial in engineering intervention like the construction of check dams to digest the initial thrust of the flow and other remedial measures designed for vulnerable slope protection.


2020 ◽  
Vol 23 (02) ◽  
pp. 518-533
Author(s):  
Manuel Gomes Correia ◽  
João Carlos von Hohendorff Filho ◽  
Denis José Schiozer

2015 ◽  
Vol 18 (04) ◽  
pp. 481-494 ◽  
Author(s):  
Siavash Nejadi ◽  
Juliana Y. Leung ◽  
Japan J. Trivedi ◽  
Claudio Virues

Summary Advancements in horizontal-well drilling and multistage hydraulic fracturing have enabled economically viable gas production from tight formations. Reservoir-simulation models play an important role in the production forecasting and field-development planning. To enhance their predictive capabilities and to capture the uncertainties in model parameters, one should calibrate stochastic reservoir models to both geologic and flow observations. In this paper, a novel approach to characterization and history matching of hydrocarbon production from a hydraulic-fractured shale is presented. This new methodology includes generating multiple discrete-fracture-network (DFN) models, upscaling the models for numerical multiphase-flow simulation, and updating the DFN-model parameters with dynamic-flow responses. First, measurements from hydraulic-fracture treatment, petrophysical interpretation, and in-situ stress data are used to estimate the initial probability distribution of hydraulic-fracture and induced-microfracture parameters, and multiple initial DFN models are generated. Next, the DFN models are upscaled into an equivalent continuum dual-porosity model with analytical techniques. The upscaled models are subjected to the flow simulation, and their production performances are compared with the actual responses. Finally, an assisted-history-matching algorithm is implemented to assess the uncertainties of the DFN-model parameters. Hydraulic-fracture parameters including half-length and transmissivity are updated, and the length, transmissivity, intensity, and spatial distribution of the induced fractures are also estimated. The proposed methodology is applied to facilitate characterization of fracture parameters of a multifractured shale-gas well in the Horn River basin. Fracture parameters and stimulated reservoir volume (SRV) derived from the updated DFN models are in agreement with estimates from microseismic interpretation and rate-transient analysis. The key advantage of this integrated assisted-history-matching approach is that uncertainties in fracture parameters are represented by the multiple equally probable DFN models and their upscaled flow-simulation models, which honor the hard data and match the dynamic production history. This work highlights the significance of uncertainties in SRV and hydraulic-fracture parameters. It also provides insight into the value of microseismic data when integrated into a rigorous production-history-matching work flow.


2005 ◽  
Vol 9 (4) ◽  
pp. 313-321 ◽  
Author(s):  
R. R. Shrestha ◽  
S. Theobald ◽  
F. Nestmann

Abstract. Artificial neural networks (ANNs) provide a quick and flexible means of developing flood flow simulation models. An important criterion for the wider applicability of the ANNs is the ability to generalise the events outside the range of training data sets. With respect to flood flow simulation, the ability to extrapolate beyond the range of calibrated data sets is of crucial importance. This study explores methods for improving generalisation of the ANNs using three different flood events data sets from the Neckar River in Germany. An ANN-based model is formulated to simulate flows at certain locations in the river reach, based on the flows at upstream locations. Network training data sets consist of time series of flows from observation stations. Simulated flows from a one-dimensional hydrodynamic numerical model are integrated for network training and validation, at a river section where no measurements are available. Network structures with different activation functions are considered for improving generalisation. The training algorithm involved backpropagation with the Levenberg-Marquardt approximation. The ability of the trained networks to extrapolate is assessed using flow data beyond the range of the training data sets. The results of this study indicate that the ANN in a suitable configuration can extend forecasting capability to a certain extent beyond the range of calibrated data sets.


2006 ◽  
Vol 10 (2) ◽  
pp. 241-264 ◽  
Author(s):  
D. Rickenmann ◽  
D. Laigle ◽  
B. W. McArdell ◽  
J. Hübl

Sign in / Sign up

Export Citation Format

Share Document