We study a family of numerical schemes applied to a class of multiscale systems of stochastic differential equations. When the time scale separation parameter vanishes, a well-known homogenization or Wong--Zakai diffusion approximation result states that the slow component of the considered system converges to the solution of a stochastic differential equation driven by a real-valued Wiener process, with Stratonovich interpretation of the noise. We propose and analyse schemes for effective approximation of the slow component. Such schemes satisfy an asymptotic preserving property and generalize the methods proposed in a recent article. We fill a gap in the analysis of these schemes and prove strong error estimates, which are uniform with respect to the time scale separation parameter.