This invention relates to a method of estimating fluxes for the processes of matter and field transport through fluid-saturated or gas-saturated porous solid. The method comprises obtaining three-dimensional porous solid images by, but not limited, X-ray microtomography, 3D NMR imaging, 3D reconstruction from petrographic thin-section analysis etc., digital processing and morphological analysis of the 3D core images by consecutive application of the image filtering, segmentation and multiple property recognition for obtaining digital 3D models of porous solid samples and performing a set of morphological and geometrical statistical property analysis. For the above mentioned 3D model (models) heat, mass, chemical and electric fluxes are modeled (separately or in combination) under given boundary conditions by means of numerical solver. The new models, which are statistically equivalent to the abovementioned model (models) are generated by means of random field and stochastic geometry theory; heat, mass, chemical and electric fluxes are simulated for new models. The obtained fluxes are averaged over realizations to be used in macroscopic calculations.