In this paper, the stochastic parameters of the effective thermal conductivity of multilayer composites are considered. The examined specimens of composites were built with a different number of layers and each had a different saturation density of a composite matrix with fibers. For each case of laminate built with a prescribed number of layers and assumed saturation density, 10,000 tests of its effective thermal conductivity were carried out using numerical experiments. It was assumed that the fibers located in each layer were rectilinear, had a circular cross-section and that they could take random positions in their repeatable volume elements (RVEs). In view of the mentioned assumptions, the heat flux passing throughout a cross-section of a composite sample, perpendicular to the fibers’ direction, was considered. The probability density functions were fitted to the obtained data and then the chosen stochastic parameters of the effective thermal conductivity coefficients were determined.