A stochastic mathematical model and finite element framework by combination ofKarhunen-Lo eve expansion and polynomial chaos expansion of uncertainty quantification ofnatural convection in random media were developed based on stochastic theory and heattransfer theory in porous media.In this approach， the input stochastic field(porosity) is re-presented by the Karhunen-Lo eve expansion and the output fields(e.g.， velocity， pressure，temperature， and streamfunction) are expressed by the polynomial chaos expansion.Usingthe spectral decomposition， the stochastic problem in random porous media is reformulatedto a set of deterministic problems to be solved for each polynomial chaos.A stochastic pro-jection method is implemented to obtain the chaos coefficients in the corresponding deter-ministic governing equations and the statistics of the numerical solution are obtained.Theprediction results are compared against results obtained using a MonteCarlo method.Excel-lent agreement between these results indicates the efficiency and accuracy of the proposedmethod.