Upscaling of Constitutive Relations I Unsaturated Heterogeneous Porous Media

When numerical model are used for modeling field scale flow and transport processes in the subsurface, the problem of ''upscaling'' arises. Typical scales, corresponding to spatial resolutions of subsurface heterogeneity in numerical models, are generally much larger than the measurement scale of the parameters and physical processes involved. The upscaling problems is, then, one of assigning parameters to gridblock scale based on parameter values measured on small scales. The focus of this study is to develop an approach to determine large-scale (upscaled) constitutive relations (relationships among relative permeability, capillary pressure and saturation) from small-scale measurements for porous media for a range of air entry values that are typical for the tuff matrix in the unsaturated zone of Yucca Mountain. For porous media with large air entry values, capillary forces play a key role in determining spatial water distribution at large-scales. Therefore, a relatively uniform capillary pressure approximately exists even for a large gridblock scale under steady state flow conditions. Based on these reasoning, we developed formulations that relate upscaled constitutive relations to ones measured at core-scale. Numerical experiments with stochastically generated heterogeneous porous media were used to evaluate the upscaling formulations