which is utilized to connect the electrical equipment bodies to a neutral reference node. Consequently, it could decrease the touch and step voltages during or after a contingency, such as a fault occurrence or lightning strokes inside the substation or its proximity. In the second perspective, the principle is to maintain insulation and to ensure the safety of electrical equipment. Then, in order to protect the equipment insulation and the movement of the traveling waves due to the lightning strike and the switching in the transmission system of power, there must be a parallel path to a reliable grounding system, to discharge these currents with very low impedance [5]. Since the ground impedance of the grounding grid is considerable, it will prevent rapid discharge and cause a hypothetically dangerous local potential at the point of lightning evacuation and increase ground potential rise (GPR) [6; 7]. The transient behavior of grounding systems directly and significantly affects the performance of the electrical system under SC fault conditions or lightning strike collisions [8]. Due to the fact that most power network modeling to analyze the electromagnetic transient situations in software packages such as electro-magnetic transient program (EMTP) are in lumped circuits, the lack of a proper method with the precision that can provide this criterion has led to modelling the grounding grids with simple passive RLC parameters. There are some investigations to use the models presented in [9; 10; 11; 12], such as circuit theory and distributed transmission lines. In recent years, the researchers have been trying to provide a suitable model for grounding systems that also behave appropriately in electromagnetic transient states. However, these studies have usually succeeded in providing models that are either usable at low frequencies or should be considered a constant classification for all grounding grids [13; 14]. Thoroughly, various modeling approaches have been presented so far, which, based on numerical methods, investigate the behavior of the grounding systems at high frequencies. For instance, to implement and simulate the grounding system into an alternative transient program (ATP), a simple circuit based model is proposed by [PERSON] et al. in [15]. In this technique, after dividing the earth into \(n\) pieces, each component is modeled with an equivalent circuit with lumped elements, including self-inductance (M) and mutual inductances (L) and earth leakage conductance (G). The analysis based on electromagnetic fields is the most complicated method to model the high-frequency behavior of the grounding grids, because [PERSON]'s equations should be solved with the least assumptions and initializations [16; 17; 18; 19]. Therefore, by solving the [PERSON] equations, it is possible to calculate the electromagnetic fields at any points of the electrodes, and then, the values of voltage and current can be calculated without approximation. [PERSON] in [18; 19; 20] has presented this strategy for the first time, in which the potential in different locations will be measurable on the surface of the conductors into the Earth. Then, the inverse fast Fourier transform (IFFT) tool was used to calculate the grounding system parameters in the time domain. However, the IFFT tool is approximated due to the dependence on the appropriate window selection for a sequential sampling of the signal (in the frequency domain). Moreover, the observance of the distance between sampling steps is necessary, in accordance with the Nyquist law, results in a large volume of specimens in the frequency domain samples that are required to achieve the time domain parameters [21]. Another problem is that it needs to be corrected at low sampling frequencies, and it cannot be used with an inverse Fourier transform directly. Since the theory of the electromagnetic field is based on the method of moments (MoM), and solving [PERSON]'s equations and has minimum assumptions, it has a high accuracy to determine the grounding grid variables.