在2005年至2015年间，BL Lac天体PG 1553+113多次观测到超高能（VHE；> 100 GeV）伽马射线耀发，这些耀发被切伦科夫望远镜High Energy Stereoscopic System (HESS)、Major Atmospheric Gamma Imaging Cherenkov (MAGIC)以及Very Energetic Radiation Imaging Telescope Array System (VERITAS)探测到。尽管其红移（z）存在不确定性，但多项研究尝试估算这一值。本研究利用四种不同的星系际背景光（EBL）模型结合光致强子框架，分析了PG 1553+113的十七组独立观测到的VHE伽马射线能谱。通过全局$\chi^2$拟合对所有观测数据进行处理，以确定每种EBL模型的最佳拟合红移。此外，计算了在所有EBL模型下红移的置信水平（CL）区间。结果表明，光致强子框架能够有效描述所有观测到的能谱。在各种EBL模型中，Saldana-Lopez等（2021）提出的模型获得了最低的总$\chi^2$值58.82（自由度为58），而Dominguez等（2011）提出的模型对应的最高值为75.67。对于PG 1553+113红移的统计误差，Saldana-Lopez等（2021）模型的95%置信区间为0.500 < z < 0.537，而Dominguez等（2011）模型的置信区间为0.491 < z < 0.527。这些结果突显了光致强子方法在解释观测到的VHE伽马射线能谱方面的一致性。

Between 2005 and 2015, the BL Lacertae object PG 1553+113 exhibited multiple very high-energy (VHE; > 100 GeV) gamma-ray flares, which were detected by the Cherenkov telescopes, High Energy Stereoscopic System (HESS), Major Atmospheric Gamma Imaging Cherenkov (MAGIC), and the Very Energetic Radiation Imaging Telescope Array System (VERITAS). Despite the uncertainty surrounding its redshift (z), various studies have sought to estimate this value. In this study, seventeen independently observed VHE gamma-ray spectra of PG 1553+113 are analyzed using four distinct extragalactic background light (EBL) models alongside the photohadronic framework. A global $\chi^2$ fit is applied to all observational data to determine the best-fitting redshift for each EBL model. Additionally, confidence level (CL) intervals for the redshift are calculated across all EBL models. The findings demonstrate that the photohadronic framework effectively describes all observed spectra. Among the EBL models, the lowest total $\chi^2$ value of 58.82 (with 58 degrees of freedom) was obtained using the model from Saldana-Lopez et al. (2021), while the highest value of 75.67 was associated with the model from Dominguez et al. (2011). The 95 per cent CL intervals for the statistical error of the redshift of PG 1553+113 are 0.500 < z < 0.537 for the Saldana-Lopez et al. (2021) model and 0.491 < z < 0.527 for the Dominguez et al. (2011) model. These results highlight the consistency of the photohadronic approach in interpreting the observed VHE gamma-ray spectra.