The main aim of designing a safe grounding system is to provide a low impedance path to the flow of the earth fault currents without exceeding the operational constraints and equipment limits that will ensure electrical continuity. One of the most significant and well-known parameters to design a safe grounding system in power systems is the exact determination of the maximum earth fault current division factor. This paper presents a simplified and accurate method to calculate the earth fault current division factor in different states of the earth fault occurrence within and the vicinity of the substation under study. In the proposed method, a hybrid overhead-cable line, the impact of the frozen soil, mutual coupling between phase conductors and guard wires, grounding system resistance of adjacent substations, different tower footing resistances, the impact of the phase conductor impedance, and different spans in transmission lines can be considered. In addition, a closed-form formulation is also developed for the occurrence of the earth fault on the transmission line. Finally, details of the analysis results of this study have been compared with other methods in the literature. The validity and accuracy of the proposed approach also have been assayed and confirmed in details.