A 3D Computational Study of Effective Medium Methods Applied to Fractured Media

• Published in: Transport in porous media 2013-10, Vol.100 (1), p.115-142
• Main Authors: Sævik, Pål Næverlid, Berre, Inga, Jakobsen, Morten, Lien, Martha
• Format: Article
• Language: English
• Subjects: Apertures; Asymmetry; Asymptotic methods; Asymptotic properties; Civil Engineering; Classical and Continuum Physics; Computer simulation; Curve fitting; Earth and Environmental Science; Earth Sciences; Earth, ocean, space; Effective medium theory; Ellipsoids; Engineering and environment geology. Geothermics; Exact sciences and technology; Finite element method; Fractures; Geotechnical Engineering & Applied Earth Sciences; Hydrocarbons; Hydrogeology; Hydrology. Hydrogeology; Hydrology/Water Resources; Inclusions; Industrial Chemistry/Chemical Engineering; Inhomogeneity; Methods; Numerical stability; Percolation theory; Permeability; Pollution, environment geology; Predictions; Sedimentary rocks; Stability analysis
• ISSN: 0169-3913; 1573-1634
• DOI: 10.1007/s11242-013-0208-0

This work evaluates and improves upon existing effective medium methods for permeability upscaling in fractured media. Specifically, we are concerned with the asymmetric self-consistent, symmetric self-consistent, and differential methods. In effective medium theory, inhomogeneity is modeled as ellipsoidal inclusions embedded in the rock matrix. Fractured media correspond to the limiting case of flat ellipsoids, for which we derive a novel set of simplified formulas. The new formulas have improved numerical stability properties, and require a smaller number of input parameters. To assess their accuracy, we compare the analytical permeability predictions with three-dimensional finite-element simulations. We also compare the results with a semi-analytical method based on percolation theory and curve-fitting, which represents an alternative upscaling approach. A large number of cases is considered, with varying fracture aperture, density, matrix/fracture permeability contrast, orientation, shape, and number of fracture sets. The differential method is seen to be the best choice for sealed fractures and thin open fractures. For high-permeable, connected fractures, the semi-analytical method provides the best fit to the numerical data, whereas the differential method breaks down. The two self-consistent methods can be used for both unconnected and connected fractures, although the asymmetric method is somewhat unreliable for sealed fractures. For open fractures, the symmetric method is generally the more accurate for moderate fracture densities, but only the asymmetric method is seen to have correct asymptotic behavior. The asymmetric method is also surprisingly accurate at predicting percolation thresholds.