We present a model-dependent method to estimate the red shifts of three TeV BL Lac objects (BL Lacs) through fitting their (quasi-) simultaneous multi-waveband spectral energy distributions (SEDs) with a one-zone leptonic synchrotron self-Compton model. Considering the impact of electron energy distributions (EEDs) on the results, we use three types of EEDs to fit the SEDs: a power-law EED with exponential cut-off (PLC), a log-parabola (PLLP) EED and the broken power-law (BPL) EED. We also use a parameter alpha to describe the uncertainties of the extragalactic background light models, as in Abdo et al. We then use a Markov chain Monte Carlo method to explore the multi-dimensional parameter space and obtain the uncertainties of the model parameters based on the observational data. We apply our method to obtain the red shifts of three TeV BL Lac objects in the marginalized 68 per cent confidence, and find that the PLC EED does not fit the SEDs. For 3C66A, the red shift is 0.14-0.31 and 0.16-0.32 in the BPL and PLLP EEDs. For PKS1424+240, the red shift is 0.55-0.68 and 0.55-0.67 in the BPL and PLLP EEDs. For PG1553+113, the red shift is 0.22-0.48 and 0.22-0.39 in the BPL and PLLP EEDs. We also estimate the red shift of PKS1424+240 in the high stage to be 0.46-0.67 in the PLLP EED, roughly consistent with that in the low stage.