The frequency comb which is characterized by equally-spaced frequency lines with high mode coherence has received much attention since its first demonstration in near-infrared and optical frequency range. In the terahertz frequency range, the electrically-pumped terahertz quantum cascade laser (THz QCL) based on semiconductors is an ideal candidate for achieving frequency comb operation in a frequency range between 1 THz and 5 THz. The group velocity dispersion (GVD) is a key factor for the frequency comb. A higher GVD can pull the frequencies from their equidistant values and limit the comb bandwidth. Therefore the laser dispersion needs to be compensated for in order to make the total GVD sufficiently low and flat, such as using a Gires-Tournois interferometer (GTI) or the double chirped mirror (DCM). However, a successful design still depends on the knowledge of the total GVD in the laser. In this paper, we show how to calculate the GVD in metal-metal waveguide THz QCLs by taking into account the dispersions from the GaAs material, the waveguide, and the laser gain, which conduces to the understanding of the frequency comb behavior. The waveguide loss is modelled by the finite element method. The loss due to intersubband absorption is calculated by Fermi's gold rule. All the losses, i.e., waveguide loss, mirror loss, and intersubband absorption loss, are summed up to calculate the clamped gain. The material loss can be calculated by using the reststrahlen band model. Because of these losses and gain, the refractive index needs to be replaced by a complex refractive index. The real part of the complex refractive index is the refractive index, which can be calculated from the Kramers-Kronig relationship that connects the loss or gain with the refractive index. Then the GVD introduced by the material loss, waveguide loss, and clamped gain can be finally calculated. The results show that the total GVD of THz QCL is approximately –8 × 10<sup>5</sup>~8 × 10<sup>5</sup> fs<sup>2</sup>/mm which is strongly determined by the clamped gain. Finally, the developed numerical model is employed to study the dispersion compensation effect of a GTI mirror which is coupled into a QCL gain cavity. The design of the THz QCL based on GTI structure is more flexible and feasible than that of the DCM. The result shows that by carefully designing the geometry of GTI, the dispersion of a THz QCL can be compensated for, thus achieving the broadband terahertz frequency combs.