No abstract

Heterogeneous porous media are generally characterized by many scales of variations. A methodology is required for transporting information from one scale to another. Also, grid refinement techniques require that we are able to assign effective permeabilities on grid-blocks of different sizes. Recently, methods for calculating effective permeabilities on blocks of arbitrary size have been developed for models in groundwater-flow. We have extended these methods to a two-phase model, where a new dispersion coefficient is derived. Based on operator-splitting techniques and adaptive grid-refinement, we have developed a discretization strategy able to resolve local phenomena in space and time. The solution procedure includes the concept of domain decomposition where the global solution is constructed from local computations on sub-domains. This approach allows the resolution and even the entire solution strategy to depend on local properties of the model. The algorithm is well suited for parallel computer architectures as well as integration into existing large scale simulators. In the present paper these techniques are tested for heterogeneous, two-dimensional two-phase models. Numerical experiments based on fine-scale discretization are compared with the results from the calculation based on averaged equations.