This paper presents a hybrid numerical model for computation of the phase-to-ground fault current distribution between soil and neutral conductors of grounding system. Hybrid model is combination of electromagnetic model and transmission line model. In numerical realization, the finite element technique is used. Using electromagnetic model, conductive coupling between all buried non-insulated conductive parts of grounding system is taken into account. Using transmission line model, inductive coupling between phase conductors and ground wires of parallel overhead power lines is used. Similar transmission line approach is used for parallel power cable lines and cable line grounding wires. In computer programs developed for ground fault current distribution, the conductive coupling between all buried non-insulated conductive parts of grounding system can be taken into account or neglected. From the numerous numerical examples, it can be deduced that one can introduce a significant error by neglecting the conductive coupling. In presented numerical example, error caused by neglecting the conductive coupling between all buried non-insulated conductive parts of grounding system can be observed in rms values and phase angles of currents and potentials.