We prove the convergence of an implicit time Euler scheme for a 2 dimensional Boussinesq model on the torus subject to a random perturbation. The rate of convergence in probability is polynomial and that in $L^2(Omega)$ is logarithmic. The key ingredients are the time regularity of the velocity and temperature, and bounds of various moments of both processes, as well as of their discretized versions, uniformly on the time interval. This is joint work with Hakima Bessaih.