An extension of the replica Ornstein-Zernike (ROZ) equations for partly quenched polydisperse systems is presented. Explicit calculations have been performed for a hard sphere fluid confined by a polydisperse hard sphere disordered matrix. A numerical solution of the ROZ equations makes use of the orthonormal polynomials with the weight function corresponding to the distribution function of the diameters of matrix species.