Practical Aspects of Upscaling Geocellular Geological Models for Reservoir Fluid Flow Simulations: A Case Study in Integrating Geology, Geophysics, and Petroleum Engineering Multiscale Data from the Hunton Group

Optimal upscaling of a high-resolution static geologic model that reflects the flow performance of the reservoir is important for reasons such as time and calculation efficiency. In this study, we demonstrate that honoring reservoir heterogeneity is critical in predicting accurate production and reducing the time and cost of running reservoir flow simulations for the Hunton Group carbonate. We integrated three-dimensional (3D) seismic data, well logs, thin sections, outcrops, multiscale fracture studies, discrete fracture networks, and geostatistical methods to create a 100 x 150 x 1 ft gridded representative geologic model. We calibrated our gridded porosity and permeability parameters, including the evaluation of fractures, by history matching the simulated production rate and cumulative production volumes from a baseline fine-scale model generated from petrophysical and production data obtained from five wells. We subsequently reperformed the simulations using a suite of coarser grids to validate our property upscaling workflow. Compared to our baseline history matching, increasing the horizontal grid cell sizes (i.e., horizontal upscaling) by factors of 2, 4, 8, and 16 results in cumulative production errors ranging from +0.5% for two time (2x) coarser to +74.22% for 16x coarser. The errors associated with vertical upscaling were significantly less, i.e., ranging from +0.5% for 2x coarser to +10.8% for 16x coarser. We observed higher production history matching errors associated with natural fracture size. Results indicate that greater connectivity provided by natural fracture length has a larger effect on production compared to the higher permeability provided by larger apertures. We also estimated the trade-off between accuracy and run times using two methods: (a) using progressively larger grid cell sizes; (b) applying 1, 5, 10, and 20 parallel processes. Computation time reduction in both scenarios may be described by simple power law equations. Observations made from our case study and upscaling workflow may be applicable to other carbonate reservoirs.