Search citation statements
Paper Sections
Citation Types
Year Published
Publication Types
Relationship
Authors
Journals
The reduction from three‐ to two‐dimensional analysis of the permeability of a fractured rock mass introduces errors in both the magnitude and direction of principal permeabilities. This error is numerically quantified for porous rock by comparing the equivalent permeability of three‐dimensional fracture networks with the values computed on arbitrarily extracted planar trace maps. A method to compute the full permeability tensor of three‐dimensional discrete fracture and matrix models is described. The method is based on the element‐wise averaging of pressure and flux, obtained from a finite element solution to the Laplace problem, and is validated against analytical expressions for periodic anisotropic porous media. For isotropic networks of power law size‐distributed fractures with length‐correlated aperture, two‐dimensional cut planes are shown to underestimate the magnitude of permeability by up to 3 orders of magnitude near the percolation threshold, approaching an average factor of deviation of 3 with increasing fracture density. At low‐fracture densities, percolation may occur in three dimensions but not in any of the two‐dimensional cut planes. Anisotropy of the equivalent permeability tensor varies accordingly and is more pronounced in two‐dimensional extractions. These results confirm that two‐dimensional analysis cannot be directly used as an approximation of three‐dimensional equivalent permeability. However, an alternative expression of the excluded area relates trace map fracture density to an equivalent three‐dimensional fracture density, yielding comparable minimum and maximum permeability. This formulation can be used to approximate three‐dimensional flow properties in cases where only two‐dimensional analysis is available.
The reduction from three‐ to two‐dimensional analysis of the permeability of a fractured rock mass introduces errors in both the magnitude and direction of principal permeabilities. This error is numerically quantified for porous rock by comparing the equivalent permeability of three‐dimensional fracture networks with the values computed on arbitrarily extracted planar trace maps. A method to compute the full permeability tensor of three‐dimensional discrete fracture and matrix models is described. The method is based on the element‐wise averaging of pressure and flux, obtained from a finite element solution to the Laplace problem, and is validated against analytical expressions for periodic anisotropic porous media. For isotropic networks of power law size‐distributed fractures with length‐correlated aperture, two‐dimensional cut planes are shown to underestimate the magnitude of permeability by up to 3 orders of magnitude near the percolation threshold, approaching an average factor of deviation of 3 with increasing fracture density. At low‐fracture densities, percolation may occur in three dimensions but not in any of the two‐dimensional cut planes. Anisotropy of the equivalent permeability tensor varies accordingly and is more pronounced in two‐dimensional extractions. These results confirm that two‐dimensional analysis cannot be directly used as an approximation of three‐dimensional equivalent permeability. However, an alternative expression of the excluded area relates trace map fracture density to an equivalent three‐dimensional fracture density, yielding comparable minimum and maximum permeability. This formulation can be used to approximate three‐dimensional flow properties in cases where only two‐dimensional analysis is available.
SUMMARYThe problem of finite element simulation of incompressible fluid flow in porous medium is considered. The porous medium is characterized by the X‐ray microtomography technique in three dimensions. The finite calculus‐based stabilization technique is reviewed to implement the equal order finite element interpolation functions for both velocity and pressure. A noble preconditioner, the nodal block diagonal preconditioner, is considered whose performance is thoroughly investigated. Combining this preconditioner with a standard iterative solver during the computational homogenization procedure, it is possible to carry out the large‐scale fluid flow simulation for estimating permeability of the porous medium with reasonable accuracy and reliability. Copyright © 2013 John Wiley & Sons, Ltd.
Upscaling methods that need to solve local problems subject to boundary conditions are addressed in this article. We define a new upscaling method based on optimization problems, which can take into account general boundary conditions applied to local problems. The determination of upscaled permeability leads to minimizing the difference of dissipated energies (or averaged velocity) at fine and large scale. Using optimal control techniques, we obtain an effective computing algorithm that allows us to recover, with classical boundary conditions, the well-known results. The uniqueness issue is tackled for the optimization problems introduced in our approach. We show that the method is stable with respect to G-convergence, a property that establishes a link with homogenization theory, and finally, 2D numerical experiments are presented.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.