Gridding and Upscaling*
Single-phase flow-based upscaling on Cartesian grids
single coarse Cartesian gridblock of size DX DY DZ, Darcy’s Law directly
gives the permeability in each direction Kx, Ky, and Kz (using a unit
viscosity single phase incompressible fluid) as Qx DX/( DY DZ DPx), Qy DY/(DX
DZ DPy), and Qz DZ/(DX DY DPz), respectively. The Q/DP rate/pressure
difference ratios come from separate solutions in each direction on the
corresponding (nx,ny,nz) fine grid system, and they may be considered as
boundary condition/solution or as solution/boundary condition – these
definitions and boundary condition values make no difference. In that fine
grid system, the equivalent boundary conditions are effectively applied
using planar wells on the edges, for an exactly equivalent physical system.
method is not [has not been] generally applicable to non-orthogonal grids such as corner
point, SKUA, or unstructured, since Darcy's Law does not give the coarse
block solution for any discretization other than Cartesian.
So, how are the grids in use today upscaled? We
suspect that usually they are not.
It seems that most users today simply generate their grids at the desired
scale without being able to determine the optimal one. Grid size sensitivity
is usually skipped, but if done consists only of comparing results obtained
on different scale grids, which are almost guaranteed to be different
without any of the conditioning of upscaling that makes it effective. That
leads to the wrong conclusion that simulation at the geological scale is
required, or to severely increased uncertainty, in which upscaled rock properties are essentially
unknown, severely inhibiting our abilities.
* taken from an SPE Reservoir Simulation Technical Interest Group
discussion, October 28, 2011 (deleted by SPE, 2014)
New (April 2016)
We have solved the
single-phase flow-based upscaling problem for non-orthogonal grids,
including corner point and SKUA, by adapting our Cartesian upscaling
program. Significant usability and workflow improvements are provided
by integration with Sensor. Application is restricted to unfaulted,
single-porosity reservoirs, and to the matrix portion of dual
permeability reservoirs. Non-orthogonal grids require input
of grid block pore volumes and transmissibilities that are calculated by the
For fractured dual
porosity/permeability reservoirs, upscale as a single porosity model
to the scale indicated by your manual dual porosity/permeability (dp)
upscaling - comparing discrete sp vs. dp results - then add DUAL option and
See or ask about the new Sensor UPSCALE
feature. Options include stopping the run following upscaling and
writing the coarse grid file, or continuing the run with either the upscaled
coarse grid or the original fine grid.
are obvious, not only for the solution of the (single-phase) upscaling
problem for non-orthogonal systems, but for potentially achieving the goal
of quantifying uncertainty in probabilistic forecasts and optimizations, even for field-scale models.
Examples are coming soon.
No matter what
hardware systems you have, when simulation is required to answer the
question(s) and where Sensor is applicable, nothing is more efficient or
effective than running the
optimal number of tuned serial Sensor realizations in parallel (see
Parallel Reservoir Simulation?) That
optimal number is usually equal to the number of effective cores.