SUMMARYTo be able to use a limited number of relatively large grid cells in numerical oil reservoir simulators and groundwater models, upscaling of the absolute permeability is frequently applied. The spatial ÿne-scale permeability distribution, which is generally obtained from geological and geostatistical models, is incorporated in the relatively large grid cells of the numerical model. If the porous medium may be approximated as a periodic medium, upscaling can be performed by the homogenization method. Numerical homogenization gives rise to an approximation error. The complementarity between the conformal-nodal ÿnite element method and the mixed-hybrid ÿnite element method has been used to quantify this error. The two methods yield, respectively, upper and lower bounds for the eigenvalues of the coarse-scale permeability tensor. Results of numerical experiments obtained using tetrahedral meshes are shown both in the far ÿeld and in the near well region.