Penetration level of solar photovoltaic (PV) energy in the utility network is steadily increasing. This changes the fault level and causes protection problems. Furthermore, multi-tapped structure of distribution network deployed to integrate solar PV energy to the grid and supplying loads at the same time also raised the protection challenges. Hence, this manuscript is aimed at introducing an algorithm to identify and classify the faults incident on the network of utilities where penetration level of the solar PV energy is high. This fault recognition algorithm is implemented in four steps: (1) calculation of Stockwell transform-based fault index (STFI) (2) calculation of Wigner distribution function-based fault index (WDFI) (3) calculation of combined fault index (CFI) by multiplying STFI and WDFI (4) calculation of index for ground fault (IGF) used to recognize the involvement of ground in a fault event. The STFI has the merits that its performance is least affected by the noise associated with the current signals and it is effective in identification of the waveform distortions. The WDFI employs energy density of the current signals for estimation of the faults and takes care of the current magnitude. Hence, CFI has the merit that it considers the current magnitude as well as waveform distortion for recognition of the faults. The classification of faults is achieved using the number of faulty phases. An index for ground fault (IGF) based on currents of zero sequence is proposed to classify the two phase faults with and without the ground engagement. Investigated faults include phase to ground, two phases fault without involving ground, two phases fault involving ground and three phase fault. Fault recognition algorithm is tested for fault recognition with the presence of noise, various angles of fault incidence, different impedances involved during faulty event, hybrid lines consisting of overhead line (OHL) and underground cable (UGC) sections, and location of faults on all nodes of the test grid. Fault recognition algorithm is also tested to discriminate the transients due to switching operations of feeders, loads and capacitor banks from the faulty transients. Performance of the fault recognition algorithm is compared with the algorithms based on discrete wavelet transform (DWT), Stockwell transform (ST) and hybrid combination of alienation coefficient and Wigner distribution function (WDF). Effectiveness of the fault recognition algorithm is established using a detailed study on the IEEE-13 nodes test feeder modified to incorporate solar PV plant of capacity 1 MW in MATLAB/Simulink. Algorithm is also validated on practical utility grid of Rajasthan State of India.