Modeling permeability in imperfectly layered porous media. I. Derivation of block-scale permeability tensor for thin grid-blocks

Practical expressions are given for the nine components of the block-scale permeability tensor of a thin block. These expressions are derived from the local-scale continuity equation and Darcy's law in an anisotropic layered porous medium. The flow problem is separated in a bottom-flux problem and a top-flux problem, both of which can be solved in essentially the same way. The bottom-flux problem has been worked out in detail, and has been separated in two parts: a vertical potential difference and a horizontal potential difference part. Each is solved with a different approach specially designed for it. Depth-averaged expressions are obtained first, after which block-scale expressions are obtained by assuming a constant depth-averaged flux. In the zeroth order, this results in the well-known Dupuit approximation in geohydrology, and the vertical equilibrium (VE) approximation in petroleum reservoir engineering. The novelty of the theory presented here stems from the application of a perturbation technique to obtain first-order corrections to these well-known results. The local-scale laws are applied in the coordinate system coinciding with the principal axes of the local-scale permeability tensor. Only in this coordinate system the local-scale permeability tensor has zero off-diagonal components. However, since the porous medium is imperfectly layered, the first-order corrections show that the off-diagonal components of the block-scale permeability tensor are not zero. Furthermore, the block-scale permeability tensor is generally nonsymmetric, which implies that a coordinate system in which the off-diagonal terms disappear does not exist.