In ground-based astronomical spectroscopic observations, there are many telluric absorption lines that are laid on the spectra of celestial objects. To study the physical properties of the celestial objects with these contaminated spectra, the telluric lines should be removed. A conventional method for removing the telluric lines is using the standard stellar spectrum as telluric line. In this paper, we introduce a technique to calculate synthetic telluric spectra and use them to remove telluric lines from a spectrum of a celestial object. We used Line-by-Line Radiative Transfer Model (LBLRTM) for calculating a synthetic spectrum and selected Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) model as atmospheric model. We apply our method to some spectra obtained by Bohyunsan Observatory Echelle Spectrograph (BOES) to show that the telluric lines are well removed from the observed spectra by our model within an accuracy of 2% which is close to the 1-sigma rms of the original spectra.