Methods for upscaling digital rock modeling data are described. Core-plug samples for pore-scale modeling are strategically chosen using whole-core minipermeability grids and conventional CT (Computed Tomography) scans. Pore models or pore-network models are used for flow modeling. Computed numerical SCAL (Special Core AnaLysis) properties are validated using laboratory-derived data, then they are used to populate borehole-scale models. Borehole-scale models use MPS (Multi-Point Statistics) to combine minipermeability grids and conventional CTscans of whole core with electrical borehole images to create 3D numerical pseudocores for each RRT (Reservoir Rock Type). SCAL properties determined from pore-scale models are distributed for each petrophysical facies in numerical pseudocores. Effective SCAL properties computed from various MPS borehole-scale realizations or models are used to populate interwell-scale models for each RRT. At the interwell scale, seismic attributes and variogram statistics from LWD (logging while drilling) data are used to populate digital rock models. Effective properties computed from flow simulations for interwell volumes are used to populate full-field scale models. At the full-field scale, outcrop analogs, sequence stratigraphy, forward stratigraphic models, diagenetic models, and basin-scale models are combined using MPS to improve flow simulations. At every stage, REVs (representative element volumes) are computed to be certain rock heterogeneities have been captured.