In order to model and evaluate the reliability of long-life high-voltage relays with small-sample fault data characteristics, a reliability analysis method integrating average rank, the minimum mean square distance empirical distribution function, and total least squares estimation is proposed. In the random truncation experiment, considering the influence of random truncation data, the average rank method is used to correct the rank of small-sample fault data; then, the optimal empirical distribution function for small-sample fault data is obtained through the minimum average square distance, which can overcome the impact of small-sample fault data randomness. Under the assumption of the Weibull distribution model, the total least squares estimation method is used for reliability model parameter estimations, and the linear correlation coefficient and d-test method are used for model hypothesis testing. If two or more distribution models pass the linear correlation coefficient test and the d-test simultaneously, the root mean square error and relative root mean square error are applied to determine the optimal reliability model. The effectiveness of this method is verified by comparing it with the maximum likelihood estimation method.