In this paper, we first prove the local existence of strong solutions to the 3D Boussinesq equations in a bounded domain with Navier boundary conditions. Then we show the global stability of strong large solutions under a suitable integral condition.