Adopting the quasi-three-dimensional(Quasi-3D) numerical method to optimize the anti-freeze design parameters of an underground pipeline usually involves heavy numerical calculations. Here, the fitting formulae between the safe conveyance distance(SCD) of a water pipeline and six influencing factors are established based on the lowest water temperature(LWT) along the pipeline axis direction. With reference to the current widely used anti-freeze design approaches for underground pipelines in seasonally frozen areas, this paper first analyzes the feasibility of applying the maximum frozen penetration(MFP) instead of the mean annual ground surface temperature(MAGST) and soil water content(SWC) to calculate the SCD. The results show that the SCD depends on the buried depth if the MFP is fixed and the variation of the MAGST and SWC combination does not significantly change the SCD. A comprehensive formula for the SCD is established based on the relationships between the SCD and several primary influencing factors and the interaction among them. This formula involves five easy-to-access parameters: the MFP, buried depth, pipeline diameter, flow velocity, and inlet water temperature. A comparison between the analytical method and the numerical results based on the Quasi-3D method indicates that the two methods are in good agreement overall. The analytic method can be used to optimize the anti-freeze design parameters of underground water pipelines in seasonally frozen areas under the condition of a 1.5 safety coefficient. Three-dimensional (Quasi-3D) numerical method to optimize the anti-freeze design parameters of an underground pipeline usually referred to heavy numerical calculations. Here, the fitting formulae between the safe conveyance distance (SCD) of a water pipeline and six influencing factors are established based on the pipeline water direction (LWT) along the pipeline axis direction. With reference to the current wide used anti-freeze design approaches for underground pipelines in seasonally frozen areas, this paper first analyzes the feasibility of applying the maximum instead of the mean annual ground surface temperature (MAGST) and soil water content (SWC) to calculate the SCD. The results show that the SCD depends on the buried depth if the MFP is fixed and the variation of the MAGST and SWC combination does not significantly change the SCD. A comprehensive formula for the SCD is established based on the relationships between the SCD and several prim This formula involves five easy-to-access parameters: the MFP, buried depth, pipeline diameter, flow velocity, and inlet water temperature. A comparison between the analytical method and the numerical results based on the Quasi-3D method indicates that the two methods are in good agreement overall. The analytic method can be used to optimize the anti-freeze design parameters of underground water pipelines in seasonally frozen areas under the condition of a 1.5 safety coefficient.