It is well known that the stratospheric quasi-biennial oscillation (QBO) is forced by equatorial waves with different horizontal/vertical scales, including Kelvin waves, mixed Rossby–gravity (MRG) waves, inertial gravity waves (GWs), and mesoscale GWs, but the relative contribution of each wave is currently not very clear. Proper representation of these waves is critical to the simulation of the QBO in general circulation models (GCMs). In this study, the vertical resolution in the Beijing Climate Center Atmospheric General Circulation Model (BCC-AGCM) is increased to better represent large-scale waves, and a mesoscale GW parameterization scheme, which is coupled to the convective sources, is implemented to provide unresolved wave forcing of the QBO. Results show that BCC-AGCM can spontaneously generate the QBO with realistic periods, amplitudes, and asymmetric features between westerly and easterly phases. There are significant spatiotemporal variations of parameterized convective GWs, largely contributing to a great degree of variability in the simulated QBO. In the eastward wind shear of the QBO at 20 hPa, forcing provided by resolved waves is 0.1–0.2 m s−1 day−1 and forcing provided by parameterized GWs is ~0.15 m s−1 day−1. On the other hand, westward forcings by resolved waves and parameterized GWs are ~0.1 and 0.4–0.5 m s−1 day−1, respectively. It is inferred that the eastward forcing of the QBO is provided by both Kelvin waves and mesoscale convective GWs, whereas the westward forcing is largely provided by mesoscale GWs. MRG waves barely contribute to the formation of the QBO in the model.