Commun. Comput. Phys., 10 (2011), pp. 1315-1332. |
On the Importance of the Stokes-Brinkman Equations for Computing Effective Permeability in Karst Reservoirs Marcin Krotkiewski ^{1*}, Ingeborg S. Ligaarden ^{2}, Knut-Andreas Lie ^{3}, Daniel W. Schmid ^{1} 1 Physics of Geological Processes, University of Oslo, P.O. Box 1048 Blindern, NO-0316 Oslo, Norway.2 SINTEF ICT, Department of Applied Mathematics, P.O. Box 124 Blindern, NO-0314 Oslo, Norway. 3 SINTEF ICT, Department of Applied Mathematics, P.O. Box 124 Blindern, NO-0314 Oslo, Norway; and Center of Mathematics for Applications, University of Oslo, P.O. Box 1053 Blindern, NO-0316 Oslo, Norway. Received 29 June 2010; Accepted (in revised version) 2 February 2011 Available online 11 August 2011 doi:10.4208/cicp.290610.020211a Abstract Cavities and fractures significantly affect the flow paths in carbonate reservoirs and should be accurately accounted for in numerical models. Herein, we consider the problem of computing the effective permeability of rock samples based on high-resolution 3D CT scans containing millions of voxels. We use the Stokes-Brinkman equations in the entire domain, covering regions of free flow governed by the Stokes equations, porous Darcy flow, and transitions between them. The presence of different length scales and large (ten orders of magnitude) contrasts in permeability leads to highly ill-conditioned linear systems of equations, which are difficult to solve. To obtain a problem that is computationally tractable, we first analyze the relative importance of the Stokes and Darcy terms for a set of idealized 2D models. We find that, in terms of effective permeability, the Stokes-Brinkman equations are only applicable for a special parameter set where the effective free-flow permeability is less than four orders of magnitude different from the matrix permeability. All other cases can be accurately modeled with either the Stokes or the Darcy end-member flows, depending on if there do or do not exist percolating free-flow regions. The insights obtained are used to perform a direct computation of the effective permeability of a rock sample model with more than 8 million cells. AMS subject classifications: 65N30, 68W10, 86-02, 86-08Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164 Key words: Upscaling, Stokes equation, Brinkman equation, porous media. *Corresponding author. Email: marcink@fys.uio.no (M. Krotkiewski), Ingeborg.Ligaarden@sintef.no (I. S. Ligaarden), Knut-Andreas.Lie@sintef.no (K.-A. Lie), schmid@fys.uio.no (D. W. Schmid) |