Aregular two- parameter perturbation analysis based upon the boundary layer approximation is presented here to study the radiative effects of both firs- and second- order resistances due to a solid matrix on the natural convection flows in porous media