Simulation and Experiment of Substation Grounding Resistance Measurement Based on Pseudorandom Signal

The grounding resistance of a substation is an important parameter that should be designed within a reasonable range to prevent operational accidents from damaging electrical equipment due to overvoltage and ensure the safe operation of an electrical system. However, simply and accurately measuring the grounding resistance of a substation has been a difficult problem faced by engineers and technicians for a long time. This paper proposes a method of denoising by applying the m-sequence correlation identification technology to the measurement of substation resistance. We established a grounding resistance model of a grounding grid and used LabVIEW to simulate it. Based on system identification and correlation function theory, pseudorandom signals or sinusoidal signals were used as excitation signals. The output results of the system were compared when pseudorandom signals and sinusoidal signals were used as excitation signals. It was verified that the grounding resistance value measured by a pseudorandom signal was closer to the actual value, which met the design requirements. Laboratory test results verify that the method of calculating grounding resistance based on the correlation analysis method is feasible.


Introduction
Grounding resistance is the main technical index of a grounding system for power plants and substations. The effectiveness and safety of a grounding system are judged by grounding resistance. However, the measurement of grounding resistance is a complicated problem because it can be affected by many factors. It is not only related to the size, shape, and ground resistivity of the grounding device [1] but also affected by the surrounding electromagnetic field, metal substances in the soil, the uniformity of earth resistivity, the measurement method, and the electrode arrangement [2,3]. In practice, substation accidents are more common in the substation system environment. An investigation showed that China has had many direct or indirect accidents caused by grounding resistance not meeting requirements or excessive grounding resistance caused by aging grounding grids, accounting for approximately 30% of grid operation accidents [4]. Therefore, before each physical substation is put into operation, it is necessary to carry out systematic and scientific power grid testing, specifically to master the grounding resistance design of the grounding grid to meet the operational requirements of the substation [5]. Moreover, after a substation is put into operation, it is necessary to regularly check the grounding resistance parameters of the grounding grid of the substation and to adjust the substation system equipment in real time to meet the actual operation requirements. If a device is unqualified, it must be replaced or repaired to ensure the safe operation of electrical equipment [6].
The measurement methods currently studied mainly include the fall-of-potential method, the heterofrequency measurement method, the moment method (MM), and the finite element method (FEM). At present, the IEEE standards and China's power industry standards are still based on the theory of the fall-of-potential method and its derivative

Materials and Methods
In this section, we introduce the commonly used wiring methods when measuring grounding resistance. The concept of related functions is briefly described. The principle of system identification based on the cross-correlation method is analyzed. The theoretical knowledge of pseudorandom signals and their autocorrelation function are mentioned.

Grounding Resistance Measurement Method
Theoretically, the grounding resistance is equal to R = V/I, where V is the ground potential rise (GPR) of the grounding system relative to the remote soil, and I is the current injected into the grounding grid. Note that the injection current is assumed to return from infinity to the current electrode (return electrode). However, in reality, it is impossible to place the current electrode at infinity when injecting current into the grounding grid. Similarly, when measuring the GPR of a grounding system, it is impossible to place the voltage electrode at infinity [25]. In the measurement of the fall-of-potential method [26][27][28][29], the position of the voltage electrode needs to be continuously adjusted. The position distribution of the voltage electrode and the current electrode in the test is shown in Figure 1, which shows the commonly used arrangement when using the fall-of-potential method to measure the grounding resistance.
power spectrum of the correlation function. The output results were compa ferent excitation signals acted on the system. Finally, the practicability of measuring device and the feasibility of the method used are verified by exp

Materials and Methods
In this section, we introduce the commonly used wiring methods wh grounding resistance. The concept of related functions is briefly described. of system identification based on the cross-correlation method is analyzed. T knowledge of pseudorandom signals and their autocorrelation function are

Grounding Resistance Measurement Method
Theoretically, the grounding resistance is equal to R = V/I, where V potential rise (GPR) of the grounding system relative to the remote soil, and I injected into the grounding grid. Note that the injection current is assumed t infinity to the current electrode (return electrode). However, in reality, it is place the current electrode at infinity when injecting current into the ground ilarly, when measuring the GPR of a grounding system, it is impossible to p age electrode at infinity [25]. In the measurement of the fall-of-potential me the position of the voltage electrode needs to be continuously adjusted. The tribution of the voltage electrode and the current electrode in the test is sho 1, which shows the commonly used arrangement when using the fall-of-pote to measure the grounding resistance. In Figure 1, E is the ground electrode. P is the potential electrode. C electrode. r is the distance between P and E. θ is the angle between P and C. this represents the 0.618 method derived from the fall-of-potential method diagram of the arrangement of the voltage electrode and the current electrod   In Figure 1, E is the ground electrode. P is the potential electrode. C is the current electrode. r is the distance between P and E. θ is the angle between P and C. When θ = 0 • , this represents the 0.618 method derived from the fall-of-potential method. A schematic diagram of the arrangement of the voltage electrode and the current electrode is shown in Figure 2.

Materials and Methods
In this section, we introduce the commonly used wiring methods wh grounding resistance. The concept of related functions is briefly described of system identification based on the cross-correlation method is analyzed. T knowledge of pseudorandom signals and their autocorrelation function are

Grounding Resistance Measurement Method
Theoretically, the grounding resistance is equal to R = V/I, where V potential rise (GPR) of the grounding system relative to the remote soil, and injected into the grounding grid. Note that the injection current is assumed infinity to the current electrode (return electrode). However, in reality, it is place the current electrode at infinity when injecting current into the ground ilarly, when measuring the GPR of a grounding system, it is impossible to age electrode at infinity [25]. In the measurement of the fall-of-potential m the position of the voltage electrode needs to be continuously adjusted. Th tribution of the voltage electrode and the current electrode in the test is sh 1, which shows the commonly used arrangement when using the fall-of-po to measure the grounding resistance. In Figure 1, E is the ground electrode. P is the potential electrode. C electrode. r is the distance between P and E. θ is the angle between P and C this represents the 0.618 method derived from the fall-of-potential method diagram of the arrangement of the voltage electrode and the current electro Figure 2.   To facilitate the calculation, the grounding grid is equivalent to a hemisphere. According to the electromagnetic field calculation method, the voltage U(d PE ) at the voltage electrode can be obtained as After Equation (1), I is the current injected into the grounding grid, U(d PE ) is the voltage between the two points P and E, ρ is the resistivity, d PE is the distance between the ground electrode and the voltage electrode, d CE is the distance between the ground electrode and the current electrode, and d CP is the distance between the voltage electrode and the current electrode. If R is the grounding resistance of the grounding grid, the voltage is where U is the voltage at the reference point. Then, the potential difference is where V is the ground potential rise (GPR) of the grounding system. Let d PE be a variable; the solution can be d PE = 0.618d CE [30,31]. That is, when the voltage electrode is located between the ground electrode and the current electrode and is 0.618d CE away from the ground electrode, an accurate voltage value can be obtained. In fact, the grounding grid structure falls between the disk electrode and the ring electrode. Therefore, when the 0.618 method is used for grounding resistance detection, usually d PE = (4 − 5)D [32,33], d PE = (0.5 − 0.6)d CE . If the size of the ground electrode is much smaller than d PE and d CE , the above relationship is still valid.

System Identification and Correlation Functions
For a signal x(t) with a period of T, the average characteristic of its own correlation degree is called the autocorrelation function R x (τ), which represents the average correlation degree between two different moments separated by τ, namely, If a linear system exists, the input signal is x(t), the output signal is y(t), and the input value x(t) at time t has an impact on the output value y(t + τ). The degree of this impact is the cross-correlation function R xy (τ) to describe, namely, The principle of using the cross-correlation method to identify the dynamic characteristics of the system is shown in Figure 3.  R (τ) is the cross-correlation function of the input signal x(t) and the output signal z(t), R (τ) is the autocorrelation function of the test signal x(t), S (ω) is the power spectrum of the cross-correlation function R (τ), and S (ω) is the power spectrum of the autocorrelation function R (τ). The impulse response function is h(t), all internal noises are modeled as noise n(t) superimposed on the response signal y(t), and the ob-  R xz (τ) is the cross-correlation function of the input signal x(t) and the output signal z(t), R x (τ) is the autocorrelation function of the test signal x(t), S xz (ω) is the power spectrum of the cross-correlation function R xz (τ), and S x (ω) is the power spectrum of the autocorrelation function R x (τ). The impulse response function is h(t), all internal noises are modeled as noise n(t) superimposed on the response signal y(t), and the observed output signal z(t) is The test signal x(t) is broadband noise, which is superimposed on the normal operating signal f(t), and the signal input to the identified system is where u(t) is the sum of the signals x(t) and f(t). Generally, f(t) and x(t) are not related to each other. For a physically achievable linear system, the input and output signals meet the following convolution relationship: By performing cross-correlation processing on the test signal x(t) and the observed output signal z(t), we obtain where R xz (τ) is the cross-correlation function of the input signal x(t) and the output signal z(t). Substituting (9) into (10) obtains After Equation (11), R xn (τ) is the cross-correlation function of the input signal x(t) and the noise n(t). R xf (τ − α) is the cross-correlation function of the input signal x(t) and the normal operating signal f(t) at time (τ − α). R x (−α) is the autocorrelation function of the test signal x(t) at time α. The test signal x(t) is not related to the noise n(t) and the production signal f(t), so Bringing this solution into (11) obtains where R x (τ) is the autocorrelation function of the test signal x(t) at time τ. The relationship expressed by (13) is called the Wiener-Hopf formula. Note that the existence of any signals other than x(t) in the system does not affect the cross-correlation result of (13). Regardless of whether other signals are applied externally or generated inside the system, as long as they are not related to the test signal, the above conclusion is valid. If the test signal x(t) is white noise, then its autocorrelation function R x (τ) is a δ function, and its power is set to K, which can be obtained by (13): The impulse response function of h(t) the system can be directly obtained from the cross-correlation function, and the frequency response function of the identified system can be obtained by the FFT of (14): where S xz (ω) is the power spectrum of the cross-correlation function R xz (τ). H(ω) is the frequency response function. From the above analysis, we can see that when the input signal is random white noise, the system impulse response can be obtained by the correlation between the input signal and the output response. The amplitude of this impulse response becomes K times the original amplitude. The phase has not changed, and the shape is consistent with the real impulse response. Therefore, using white noise as the input Sensors 2022, 22, 5777 6 of 14 signal, the system impulse response can be obtained through the cross-correlation method. However, in practice, white noise cannot be directly generated. Therefore, researchers often generate sequences with white noise properties to give them the correlation properties of random signals. In general, the m-sequence pseudorandom signal has periodicity, but as long as its period is much longer than the width of the autocorrelation function, it can be approximately considered to have the characteristics of white noise. The m-sequence is the earliest widely used pseudorandom sequence and has good randomness and balance. The m-sequence is called the maximum length shift register sequence (MLS), which is a positive-and negative-level signal. The n-stage shift register can be used to generate an m-sequence with length N = 2 n − 1. The m-sequence with an amplitude of ±1 generated by the 5-stage shift register is shown in Figure 4.
the frequency response function. From the above analysis, we can see that when the inpu signal is random white noise, the system impulse response can be obtained by the corre lation between the input signal and the output response. The amplitude of this impulse response becomes K times the original amplitude. The phase has not changed, and the shape is consistent with the real impulse response. Therefore, using white noise as the input signal, the system impulse response can be obtained through the cross-correlation method. However, in practice, white noise cannot be directly generated. Therefore, re searchers often generate sequences with white noise properties to give them the correla tion properties of random signals. In general, the m-sequence pseudorandom signal has periodicity, but as long as its period is much longer than the width of the autocorrelation function, it can be approximately considered to have the characteristics of white noise.
The m-sequence is the earliest widely used pseudorandom sequence and has good randomness and balance. The m-sequence is called the maximum length shift register se quence (MLS), which is a positive-and negative-level signal. The n-stage shift register can be used to generate an m-sequence with length N = 2 n − 1. The m-sequence with an am plitude of ±1 generated by the 5-stage shift register is shown in Figure 4.  The autocorrelation function R x (τ) of the bipolar m-sequence pseudorandom signal x(t) with a level of ±a, a clock period of ∆t and a period of T = N∆t is a triangular wave similar to the δ function. The autocorrelation function of the pseudorandom signal x(t) is When N is large enough, in the interval of 0 ∼ N∆t, R x (τ) can be approximated as a δ function [34].
The test signal parameters used for system identification are the clock period ∆t, signal length N, repetition period q, and signal amplitude ±a. Only by correctly selecting each parameter can an accurate impulse response be obtained. The following is the selection of signal length N and signal amplitude ±a. When the clock period ∆t is determined, the signal length N is selected to satisfy N∆t > T. The rule for selecting N is N∆t = (1.25 − 1.5)T s . If the adjustment time T s of the system is unknown and the maximum time constant T max of the system is known, the signal length N is selected as N = (1.25 − 1.5) * (4 − 5)T max . The amplitude a of the pseudorandom signal should be as large as possible without affecting the normal operation of the system to improve the signal-to-noise ratio.
If x(t) is a sinusoidal function s(t) superimposed with uncorrelated noise n(t), then In the formula, A is the signal amplitude, ω 0 is the signal angular frequency, ϕ is the initial phase angle of the signal, and its autocorrelation function is where R s (τ) is the autocorrelation function of the sinusoidal function s(t). R n (τ) is the autocorrelation function of the noise n(t). The long-term integration result of the first cosine function in the integrand is zero, and the second term is not a function of time t, so If n(t) is broadband noise, R n (τ) is concentrated around τ = 0. When τ is large, the amplitude A and frequency f of the signal s(t) can be measured by R x (τ). In this way, after autocorrelation processing, the amplitude and frequency of the sinusoidal signal are extracted from the noise.

Design and Result Analysis
In this section, we built a pure resistance model and reasonably selected simulation parameters. The output results of the system when the m-sequence pseudorandom signal and sinusoidal signal were used as excitation signals were compared. The simulation results of grounding resistance when the signal was superimposed with or without noise were analyzed. Finally, laboratory test verify that is feasible. Figure 5 shows the grounding grid simulation model. An m-sequence pseudorandom signal or sinusoidal signal is used as the excitation signal x(t) and injected into the grounding grid. The grounding resistance of the grounding grid is set to h(t) to obtain the response signal y(t). The response signal is superimposed with external interference noise n(t) to obtain the output signal z(t). The autocorrelation function of the excitation signal and the cross-correlation function of the excitation signal and the output signal are calculated. FFT is performed to obtain the power spectrum of the autocorrelation function and the cross-correlation function. The cross-correlation power spectrum is divided by the autocorrelation power spectrum to obtain the grounding resistance value, namely,

Simulation Model Design and Parameter Design
where S xz (ω) is the power spectrum of the cross-correlation function R xz (τ), and S x (ω) is the power spectrum of the autocorrelation function R x (τ). R is the grounding resistance. and the cross-correlation function. The cross-correlation power spectrum is divided by the autocorrelation power spectrum to obtain the grounding resistance value, namely, where S (ω) is the power spectrum of the cross-correlation function R (τ), and S (ω) is the power spectrum of the autocorrelation function R (τ). R is the grounding resistance. Then, take the average and calculate the final result, considering that in actual measurements, there is inevitable noise when using analog-to-digital conversion (ADC) to collect excitation signals and response signals. Therefore, in the simulation process, an interference noise signal n(t) is also superimposed on the input signal x(t). The interference noise signals used in the simulation are all Gaussian white noise signals.
In the simulation model, the sampling frequency of the excitation signal x(t) was F = 1 kHz, the number of sampling points N of the signal was 1000, and the frequency resolution was F = = 1. The number of FFT sampling points N and N was usually an integer power of 2, so N = 2 = 1024 was selected in the simulation experiment. The standard deviation of the superimposed Gaussian white noise was 0.5. In practice, Then, take the average and calculate the final result, considering that in actual measurements, there is inevitable noise when using analog-to-digital conversion (ADC) to collect excitation signals and response signals. Therefore, in the simulation process, an interference noise signal n(t) is also superimposed on the input signal x(t). The interference noise signals used in the simulation are all Gaussian white noise signals. In the simulation model, the sampling frequency of the excitation signal x(t) was F s = 1 kHz, the number of sampling points N 0 of the signal was 1000, and the frequency resolution was F = F s N 0 = 1. The number of FFT sampling points N ≥ F s F and N was usually an integer power of 2, so N = 2 10 = 1024 was selected in the simulation experiment. The standard deviation of the superimposed Gaussian white noise was 0.5. In practice, the noise superimposed by the collected signal is not overly large to better verify the characteristics of pseudorandom signals that are easy to separate from noise. We chose a larger intensity of noise.

Comparative Analysis of Sinusoidal Signal and Pseudorandom Signal
The main LabVIEW program was mainly composed of an excitation signal generation module, a signal superimposed noise module, an autocorrelation and cross-correlation signal module, an FFT module, and a result calculation module. First, set the standard value of resistance R1; second, inject the m-sequence pseudorandom signal or sinusoidal signal superimposed with noise into the pure resistance model as the excitation signal to obtain the current signal and voltage signal at both ends of the resistance. Then, the data were processed based on the correlation identification method to get the measurement result R2. By comparing the values of R1 and R2, the accuracy of R2 was verified.
The amplitude of the sinusoidal signal was A = 4, and the frequency was 10 Hz; the simulation waveform is shown in Figure 6. In the simulation experiment, the number of signal sampling points was selected to be 1000. If the signal is periodic, then its autocorrelation function is also periodic. The autocorrelation function of the sine signal is the cosine function of the energy attenuation. When τ = 1000, R x (1000) = A 2 2 = 8, which conforms to the theoretical value. If the number of pseudorandom signal levels was n = 5, the signal length was N = 2 n − 1 = 31, and the amplitude was a = 4. The simulation waveform is shown in Figure 7. When τ = 1000, R x (1000) = a 2 = 16. The other parameters of the two excitation signals and other parameters of the program were the same.    Figure 8 shows a curve diagram of the grounding resistance simulation results. Figure 8a shows that the resistance is a constant value when there is no interference noise, and the curve is a straight line. Figure 8b is the resistance curve after the interference noise is superimposed. Because of the presence of noise interference, the simulation result curve has a spike, and the final result can be approximated to the actual value by averaging.    Figure 8 shows a curve diagram of the grounding resistance simulation results. Figure 8a shows that the resistance is a constant value when there is no interference noise, and the curve is a straight line. Figure 8b is the resistance curve after the interference noise is superimposed. Because of the presence of noise interference, the simulation result curve has a spike, and the final result can be approximated to the actual value by averaging.   Figure 8a shows that the resistance is a constant value when there is no interference noise, and the curve is a straight line. Figure 8b is the resistance curve after the interference noise   Figure 8 shows a curve diagram of the grounding resistance simulation results. Figure 8a shows that the resistance is a constant value when there is no interference noise, and the curve is a straight line. Figure 8b is the resistance curve after the interference noise is superimposed. Because of the presence of noise interference, the simulation result curve has a spike, and the final result can be approximated to the actual value by averaging.  The grounding resistance in the substation grounding grid is usually less than 4 Ω, so a resistance value between 0-10 Ω was randomly selected for simulation. Table 1 shows the randomly selected resistance values and the grounding resistance values obtained by the simulation of the two excitation signals and their error percentages. In order to avoid stochastic uncertainty, we measured three groups of data for each measured point and calculated the average value of them. Then the error percentage between the average value and the standard value was calculated. Figure 9 shows the error percentage curve.  It can be seen in the above figures that when a sinusoidal signal is used as an excitation signal, the grounding resistance value obtained differs greatly from the actual value. Moreover, the data vary widely and are extremely unstable. The error percentage can reach 2.83%. The pseudorandom signal can be well separated from the signal. The grounding resistance value obtained by the simulation measurement is close to the standard resistance value, and the error percentage is also very small, basically below 0.29%, which is allowed in the actual device design. The grounding resistance measurement result is stable. Moreover, this result is obtained when the noise is relatively large. In practice, the superimposed noise of the signal is relatively small, which verifies that the pseudorandom signal as the excitation signal can identify the system well. There will be much interference in the surrounding environment of the substation, including harmonic interference and Gaussian interference. The pseudorandom signal is easy to separate from the noise, which shows its strong anti-interference ability.
In addition, the pseudorandom signal is a binary signal, which can be realized by turning off the switch. However, if a sinusoidal signal is used as the power supply current It can be seen in the above figures that when a sinusoidal signal is used as an excitation signal, the grounding resistance value obtained differs greatly from the actual value. Moreover, the data vary widely and are extremely unstable. The error percentage can reach 2.83%. The pseudorandom signal can be well separated from the signal. The grounding resistance value obtained by the simulation measurement is close to the standard resistance value, and the error percentage is also very small, basically below 0.29%, which is allowed in the actual device design. The grounding resistance measurement result is stable. Moreover, this result is obtained when the noise is relatively large. In practice, the superimposed noise of the signal is relatively small, which verifies that the pseudorandom signal as the excitation signal can identify the system well. There will be much interference in the surrounding environment of the substation, including harmonic interference and Gaussian interference. The pseudorandom signal is easy to separate from the noise, which shows its strong anti-interference ability.
In addition, the pseudorandom signal is a binary signal, which can be realized by turning off the switch. However, if a sinusoidal signal is used as the power supply current of the measuring device, its power supply voltage is usually not very low. The distortion will increase as the output amplitude increases. The greater the distortion is, the more inaccurate the measurement results are. This proves that pseudorandom signals can be more easily generated and accurately used for the measurement of substation grounding resistance.

Laboratory Test
According to the method proposed, a device was designed for the purpose of measuring the grounding resistance. The device consists of an m-sequence pseudorandom signal generation module, a solid-state relay module, a current sensor module, an ADC acquisition module, a serial communication module, and an upper computer LabVIEW module. STM32 single-chip microcomputer (STMicroelectronics, Geneva, Switzerland) is intended to generate the m-sequence pseudorandom signal treated as an excitation signal. The solid-state relay is equivalent to a non-contact switch purposed to separate the control end from the load end. After capturing the excitation signal, the current sensor converts it into a voltage signal, which is conducive to processing the signal for the single-chip microcomputer. As the core of the ADC acquisition circuit, the ADS1256 chip achieves the conversion between analog signals and digital signals. The upper computer LabVIEW module is responsible mainly for autocorrelation and cross-correlation operations, FFT calculation, and the display of results. A grounding resistance model was constructed in the laboratory, with the connection diagram shown in Figure 10.
In the model, the potential electrode, the current electrode, and the ground electrode were connected, with the power switched on for testing. The upper computer LabVIEW was used to display the m-sequence pseudorandom signal waveform, the voltage signal waveform, and the grounding resistance value, as shown in Figure 11. We conducted test on different resistance values, with the measurement results and their errors listed in Table 2. end from the load end. After capturing the excitation signal, the current sensor converts it into a voltage signal, which is conducive to processing the signal for the single-chip microcomputer. As the core of the ADC acquisition circuit, the ADS1256 chip achieves the conversion between analog signals and digital signals. The upper computer LabVIEW module is responsible mainly for autocorrelation and cross-correlation operations, FFT calculation, and the display of results. A grounding resistance model was constructed in the laboratory, with the connection diagram shown in Figure 10. In the model, the potential electrode, the current electrode, and the ground electrode were connected, with the power switched on for testing. The upper computer LabVIEW was used to display the m-sequence pseudorandom signal waveform, the voltage signal waveform, and the grounding resistance value, as shown in Figure 11. We conducted test on different resistance values, with the measurement results and their errors listed in Table 2.  It can be seen from the above figures and table that the error between the measurement value and the actual value is less than 4%, which falls within the allowable range. In the meantime, it was verified that the method of calculating the grounding resistance based on the correlation analysis method is feasible. In addition, the method demonstrates high precision and robustness to interference.  It can be seen from the above figures and table that the error between the measurement value and the actual value is less than 4%, which falls within the allowable range. In the meantime, it was verified that the method of calculating the grounding resistance based on the correlation analysis method is feasible. In addition, the method demonstrates high precision and robustness to interference.

Discussion
Compared with numerical calculation methods that require complex processing procedures such as dividing areas and establishing boundary conditions, the correlation identification method based on pseudorandom signal is used to measure the grounding resistance to improve efficiency. In this section, a pure resistance model is built. Specifically, the sampling frequency F s of the excitation signal x(t) and the number of sampling points N 0 of the signal in the simulation model were 1 kHz and 1000, respectively. In this paper, the superimposed noise is Gaussian white noise with a standard deviation of 0.5. Thus, the sampling frequency F s of 1 KHz will not be affected by the frequency of 50 Hz. The output results of the system were compared when pseudorandom signals and sinusoidal signals were used as excitation signals. The pseudorandom signal can be well separated from the signal. The grounding resistance value obtained by the simulation measurement is close to the standard resistance value. It verifies that the pseudorandom signal as the excitation signal can well identify the system and the laboratory test results verify that the method of calculating grounding resistance based on the correlation analysis method is feasible. Additionally, there is much interference in the surrounding environment of the substation, including harmonic interference and Gaussian interference. The pseudorandom signal is easy to separate from the noise, suggesting its strong anti-interference ability.

Conclusions
This paper proposed a simulation experiment on the grounding resistance of the substation grounding grid. First, we analyzed the layout of the lead wires when measuring the grounding resistance of the substation. Then, based on the system identification and correlation function theory, a pseudorandom signal was used as the excitation signal, and autocorrelation and cross-correlation operations were performed. Then, FFT was carried out on the correlation function to measure the grounding resistance value in the simulation model. Finally, a comparative analysis with the sinusoidal signal was carried out, which verified that the pseudorandom signal was easy to separate from the noise interference signal and reflected its strong anti-interference characteristics.
Through this simulation experiment, we carried out the design of a substation grounding resistance measurement device. Based on the theoretical methods used in the simulation experiment, considering various interferences around the substation, the hardware circuit was rationally designed. These circuits included signal acquisition circuits, power amplifier circuits, and signal conditioning circuits. Then, the algorithm used in the simulation experiment was used for data processing. Finally, a practical grounding resistance measurement device could be obtained.
In the actual measurement, to avoid the coupling problem between the current lead and the voltage lead, the short arrangement and accurate position of the lead as far as possible are required. Moreover, when arranging the electrode, technicians should avoid the underground metal and other places where electromagnetic interference will occur. Finally, high precision ADC is needed during the design of a grounding resistance measuring instrument. It can be used to realize accurate collection and conversion of the current signal and the voltage signal, thus increasing the accuracy of grounding resistance measurement.