Estimation of Fracture Porosity of Naturally Fractured Reservoirs with No Matrix Porosity Using Fractal Discrete Fracture Networks

Author(s):  
Tae Hyung Kim ◽  
David Stuart Schechter
2009 ◽  
Vol 12 (02) ◽  
pp. 232-242 ◽  
Author(s):  
Tae H. Kim ◽  
David S. Schechter

Summary Matrix porosity is relatively easy to measure and estimate compared to fracture porosity. On the other hand, fracture porosity is highly heterogeneous and very difficult to measure and estimate. When matrix porosity of naturally fractured reservoirs (NFRs) is negligible, it is very important to know fracture porosity to evaluate reservoir performance. Because fracture porosity is highly uncertain, fractal discrete fracture network (FDFN) generation codes were developed to estimate fracture porosity. To reflect scale-dependent characteristics of fracture networks, fractal theories are adopted. FDFN modeling technique enables the systematic use of data obtained from image log and core analysis for estimating fracture porosity. As a result, each fracture has its own fracture aperture distribution, so that generated FDFN are similar to actual fracture systems. The results of this research will contribute to properly evaluating the fracture porosity of NFR where matrix porosity is negligible.


Fractals ◽  
2019 ◽  
Vol 27 (01) ◽  
pp. 1940015 ◽  
Author(s):  
WEIFENG LV ◽  
GUOLIANG YAN ◽  
YONGDONG LIU ◽  
XUEFENG LIU ◽  
DONGXING DU ◽  
...  

The fracture has great impact on the flow behavior in fractured reservoirs. Fracture traces are usually self-similar and scale-independent, which makes the fractal theory become a powerful tool to characterize fracture. To obtain three-dimensional (3D) digital rocks reflecting the properties of fractured reservoirs, we first generate discrete fracture networks by stochastic modeling based on the fractal theory. These fracture networks are then added to the existing digital rocks of rock matrixes. We combine two low-permeable cores as rock matrixes with a group of discrete fracture networks with fractal characteristics. Various types of fractured digital rocks are obtained by adjusting different fracture parameters. Pore network models are extracted from the 3D fractured digital rock. Then the permeability is predicted by Darcy law to investigate the impacts of fracture properties to the absolute permeability. The permeability of fractured rock is subject to exponential increases with fracture aperture. The relationship between the permeability and the fractal dimension of fracture centers is exponential, as well as the relationship between permeability and the fractal dimension of fracture lengths.


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.


SPE Journal ◽  
2017 ◽  
Vol 22 (04) ◽  
pp. 1064-1081 ◽  
Author(s):  
Sanbai Li ◽  
Dongxiao Zhang ◽  
Xiang Li

Summary A fully coupled thermal/hydromechanical (THM) model for hydraulic-fracturing treatments is developed in this study. In this model, the mixed finite-volume/finite-element method is used to solve the coupled system, in which the multipoint flux approximation L-method is used to calculate interelement fluid and heat flux. The Gu et al. (2011) crossing criterion is extended to a 3D scenario to delineate the crossing behaviors as hydraulic fractures meet inclined natural fractures. Moreover, the modified Barton et al. (1985) model proposed by Asadollahi et al. (2010) is used to estimate the fracture aperture and model the shear-dilation effect. After being (partially) verified by means of comparison with results from the literature, the developed model is used to investigate complex-fracture-network propagation in naturally fractured reservoirs. Numerical experiments show that the key factors controlling the complexity of the induced-fracture networks include stress anisotropy, injection rate, natural-fracture distribution (fracture-dip angle, strike angle, spacing, density, and length), fracture-filling properties (the degree of cementation and permeability), fracture-surface properties (cohesion and friction angle), and tensile strength of intact rock. It is found that the smaller the stress anisotropy and/or the lower the injection rate, the more complex the fracture network; a high rock tensile strength could increase the possibility of the occurrence of shear fractures; and under conditions of large permeability of fracture filling combined with small cohesive strength and friction coefficient, shear slip could become the dominant mechanism for generating complex-fracture networks. The model developed and the results presented can be used to understand the propagation of complex-fracture networks and aid in the design and optimization of hydraulic-fracturing treatments.


SPE Journal ◽  
2013 ◽  
Vol 19 (02) ◽  
pp. 289-303 ◽  
Author(s):  
Ali Moinfar ◽  
Abdoljalil Varavei ◽  
Kamy Sepehrnoori ◽  
Russell T. Johns

Summary Many naturally fractured reservoirs around the world have depleted significantly, and improved-oil-recovery (IOR) processes are necessary for further development. Hence, the modeling of fractured reservoirs has received increased attention recently. Accurate modeling and simulation of naturally fractured reservoirs (NFRs) is still challenging because of permeability anisotropies and contrasts. Nonphysical abstractions inherent in conventional dual-porosity and dual-permeability models make them inadequate for solving different fluid-flow problems in fractured reservoirs. Also, recent technologies for discrete fracture modeling may suffer from large simulation run times, and the industry has not used such approaches widely, even though they give more-accurate representations of fractured reservoirs than dual-continuum models. We developed an embedded discrete fracture model (DFM) for an in-house compositional reservoir simulator that borrows the dual-medium concept from conventional dual-continuum models and also incorporates the effect of each fracture explicitly. The model is compatible with existing finite-difference reservoir simulators. In contrast to dual-continuum models, fractures have arbitrary orientations and can be oblique or vertical, honoring the complexity of a typical NFR. The accuracy of the embedded DFM is confirmed by comparing the results with the fine-grid, explicit-fracture simulations for a case study including orthogonal fractures and a case with a nonaligned fracture. We also perform a grid-sensitivity study to show the convergence of the method as the grid is refined. Our simulations indicate that to achieve accurate results, the embedded discrete fracture model may only require moderate mesh refinement around the fractures and hence offers a computationally efficient approach. Furthermore, examples of waterflooding, gas injection, and primary depletion are presented to demonstrate the performance and applicability of the developed method for simulating fluid flow in NFRs.


Sign in / Sign up

Export Citation Format

Share Document