Next Article in Journal
Impact of Laser Intensity Noise on Dual-Comb Absolute Ranging Precision
Previous Article in Journal
MsWH: A Multi-Sensory Hardware Platform for Capturing and Analyzing Physiological Emotional Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

College of Instrumentation & Electrical Engineering, Jilin University, Changchun 130012, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(15), 5777; https://doi.org/10.3390/s22155777
Submission received: 24 June 2022 / Revised: 25 July 2022 / Accepted: 30 July 2022 / Published: 2 August 2022
(This article belongs to the Section Physical Sensors)

Abstract

:
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.

1. 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 methods for designing grounding resistance measurement devices. The fall-of-potential method is a basic method widely used in grounding impedance measurements [7]. The 0.618 method and the 29° method are simple applications of the fall-of-potential method. The fall-of-potential method needs to arrange a voltage electrode and a current electrode during measurement. However, there are mutual inductance interference phenomena between the leads of the voltage electrode and the current electrode. Dr. Jozsef Ladanyi introduced a method of measuring grounding resistance by the fall-of-potential method and its error in field measurement [3]. J. Ma analyzed the inductive coupling between the current lead and the voltage lead under a typical lead separation distance [8]. M. A. Salam measured the soil resistivity near the test site. Based on the fall-of-potential method to measure the grounding resistance near the substation, an empirical relationship between the grounding resistance and the distance between the grid and the current electrode was proposed [9]. The heterofrequency measurement method uses a unique hardware and software anti-interference method. Under the condition that the test current frequency is close to the system power frequency resistance, stable and reliable measurement results can be obtained [10]. References [11,12] proposed a theoretical analysis and measurement system for grounding impedance based on short current leads and swept frequency alternating current (AC) sources. The disadvantage of the heterofrequency measurement method is that it still needs to place two leads of voltage and current, and it is difficult to completely avoid the mutual inductance between the measuring leads. The MM [13] and the FEM [14,15] encounter problems when dealing with large-scale grounding grid problems, such as difficulties in subdivision, large amounts of calculation, and difficulties in fully considering various factors [16]. The main disadvantage is that when the size of the grounding grid is too large, the capacity of the computer is limited, especially when the ratio of the size of the grounding grid to the size of the grounding electrode is large. Pseudorandom signals are easy to generate and implement, and their statistical characteristics are similar to those of band-limited white noise. Within a certain frequency range, the power spectrum is almost uniformly distributed [17]. The autocorrelation function of a pseudorandom signal is a triangular pulse, which is similar to the delta (δ) function. The advantage of pseudorandom signals is strong anti-interference that can be easily separated from interference noise. As long as this does not affect the normal operation of the system, the impulse response of the system can be identified. This approach is mostly used in applications such as substation [17,18] and bioimpedance measurements [19].
With an m-sequence pseudorandom signal as the excitation signal and system identification principle as the theoretical basis, the correlation identification technique can not only collect the resistance information but also mitigate the impact of noise and natural potential [20,21]. The m-sequence exhibits desirable autocorrelation characteristics. As the period of m-sequence is extended and the time width decreases, the overall autocorrelation function approaches the response function of the system increasingly, thus achieving the purpose of resistance measurement. On the basis of in-depth research into the m-sequence correlation identification method and correct simulation calculation, this paper proposes a method of denoising by applying the m-sequence correlation identification technology to the measurement of substation resistance. According to this idea, the process of denoising correlation identification was simulated, and its outcomes were analyzed [22]. Compared with the methods of numerical calculation that involve complex processing, for example, segmenting areas and establishing boundary conditions, the pseudorandom-signal-based correlation identification method purposed to measure grounding resistance is more effective in improving the efficiency and accuracy of measurement.
The main purpose of this paper is to simulate the grounding resistance model of a grounding grid based on LabVIEW software [23]. The simulation used a pseudorandom signal or sinusoidal signal as the excitation signal. An excitation signal superimposed with Gaussian white noise was injected into the grounding resistance model to obtain a response signal superimposed with Gaussian white noise. Correlation operations were performed on the excitation signal and response signal, and then FFT was performed on the correlation function [24]. The grounding resistance value could be obtained through the power spectrum of the correlation function. The output results were compared when different excitation signals acted on the system. Finally, the practicability of the designed measuring device and the feasibility of the method used are verified by experiments.

2. 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.

2.1. 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.
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.
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
U d P E = ρ I 2 π d P E ρ I 2 π ( d C E d P E )
After Equation (1),   I is the current injected into the grounding grid, U d P E is the voltage between the two points P and E, ρ is the resistivity,   d P E is the distance between the ground electrode and the voltage electrode, d C E is the distance between the ground electrode and the current electrode, and d C P is the distance between the voltage electrode and the current electrode. If R is the grounding resistance of the grounding grid, the voltage is
U = U 0 = RI ρ I 2 π d C E
where U is the voltage at the reference point. Then, the potential difference is
V = U U d PE
V = RI ρ I 2 π 1 d P E 1 d C P + 1 d C E
where V is the ground potential rise (GPR) of the grounding system.
Let d P E be a variable; the solution can be d P E = 0.618 d C E [30,31]. That is, when the voltage electrode is located between the ground electrode and the current electrode and is 0.618 d C E 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 P E = 4 5 D [32,33], d P E = 0.5 0.6 d CE . If the size of the ground electrode is much smaller than d P E and d C E , the above relationship is still valid.

2.2. 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,
R x τ = lim T 1 T 0 T x t x t + τ dt
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,
R xy τ = lim T 1 T 0 T x t y t + τ dt
The principle of using the cross-correlation method to identify the dynamic characteristics of the system is shown in Figure 3.
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
z t = y t + n t
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
u t = x t + f t
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:
y t = h t u t = 0 h α u t α d α
By performing cross-correlation processing on the test signal x t and the observed output signal z t , we obtain
R xz τ = E x τ y t + n t
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
R xz τ = R xn τ + 0 h α R x α + R xf τ α d α
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
R xn τ = R xf τ = R xf τ α = 0
Bringing this solution into (11) obtains
R xz τ = h τ R x τ
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):
h t = 1 K R xz t
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):
H ω = 1 K S xz ω
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 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 = 2n − 1. The m-sequence with an amplitude 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
R x τ = a 2 1 τ Δ t N + 1 N ,           Δ t τ Δ t         a 2 N ,           Δ t τ N 1 Δ t
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
x t = s t + n t = Asin ω 0 t + φ + n t
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
R x τ = R s τ + R n τ = lim T 1 2 T T T A 2 2 cos 2 ω 0 t ω 0 τ + 2 φ cos ω 0 τ dt + R n τ
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
R x τ = A 2 2 cos ω 0 τ + R n τ
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.

3. 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.

3.1. Simulation Model Design and Parameter Design

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,
R = S xz ω S x ω
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.
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.

3.2. 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 = 2n − 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.
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 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.

3.3. 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.
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.

3.4. 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.

4. 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.

Author Contributions

Conceptualization, Q.D. and B.Z.; Data curation, B.Z.; Formal analysis, B.Z. and H.Z.; Methodology, B.Z. and S.L.; Project administration, Q.D.; Resources, B.Z.; Software, B.Z., Q.D. and H.Z.; Supervision, Q.D.; Visualization, H.Z.; Writing—original draft, S.L. and B.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Graduate Innovation Fund of Jilin University under Grant 101832020CX188.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pratama, R.A.; Hermawan; Facta, M. Analysis of Grounding System in 150 kV Kudus Substation. In Proceedings of the 2018 5th International Conference on Information Technology, Computer, and Electrical Engineering (ICITACEE), Semarang, Indonesia, 27–28 September 2018; pp. 245–250. [Google Scholar]
  2. Unde, M.G.; Kushare, B.E. Impact of seasonal variation of soil resistivity on safety of substation grounding system. In Proceedings of the Fifth International Conference on Advances in Recent Technologies in Communication and Computing (ARTCom 2013), Bangalore, India, 20–21 September 2013; pp. 173–182. [Google Scholar]
  3. Ladanyi, J.; Smohai, B. Influence of auxiliary electrode arrangements on earth resistance measurement using the Fall-of-Potential method. In Proceedings of the 2013 4th International Youth Conference on Energy (IYCE), Siofok, Hungary, 6–8 June 2013; pp. 1–6. [Google Scholar]
  4. Li, Y. Research on Short-Distance Measurement Method of Grounding Resistance of Large-Scale Substation Grounding Grid. Master’s Thesis, Zhengzhou University, Zhengzhou, China, 2017. (In Chinese). [Google Scholar]
  5. Phayomhom, A.; Sirisumrannukul, S.; Kasirawat, T.; Puttarach, A. Safety design planning of ground grid for outdoor substations in MEA’s power distribution system. In Proceedings of the ECTI-CON2010, Chiang Mai, Thailand, 19–21 May 2010; pp. 298–302. [Google Scholar]
  6. Colak, I.; Garip, I.; Issi, F. An application of maintaining constant grounding resistance of renewable energy sources by using a dspic. In Proceedings of the 2014 Ninth International Conference on Ecological Vehicles and Renewable Energies (EVER), Monte Carlo, Monaco, 25–27 March 2014; pp. 1–6. [Google Scholar]
  7. Nassereddine, M.; Rizk, J.; Nagrial, M.; Hellany, A. Substation earth grid measurement using the fall of potential method (FOP) for a limited test area. In Proceedings of the 2014 Australasian Universities Power Engineering Conference (AUPEC), Perth, WA, Australia, 28 September–1 October 2014; pp. 1–6. [Google Scholar]
  8. Ma, J.; Dawalibi, F.P. Effects of inductive coupling between leads in ground impedance measurements using the fall-of-potential method. In Proceedings of the IEEE T&D Conference, New Orleans, LA, USA, 11–16 April 1999; Volume 1, pp. 266–271. [Google Scholar]
  9. Salam, M.A.; Noh, M. Measurement of grounding resistance with square grid and rods near substations. In Proceedings of the IEEE EPEC, London, ON, Canada, 10–12 October 2012; pp. 123–127. [Google Scholar]
  10. Liu, Y. Measurement method and application analysis of substation grounding resistance. J. Pop. Electric. 2019, 34, 21–22. (In Chinese) [Google Scholar]
  11. Zhang, B.; Zeng, R.; He, J.; Sun, W.; Yao, J.; Su, Q.; Han, S.W. Novel measurement system for grounding impedance of substations and power plants. In Proceedings of the International Conference on PowerCon 2000, Perth, WA, Australia, 4–7 December 2000; Volume 3, pp. 1389–1394. [Google Scholar]
  12. Lima, A.B.; Paulino, J.O.S.; Boaventura, W.C.; Guimarães, M.F. Transient ground impedance measurement using a very short current lead. In Proceedings of the 2013 International Symposium on Lightning Protection (XII SIPDA), Belo Horizonte, Brazil, 7–11 October 2013; pp. 177–183. [Google Scholar]
  13. Vieira, P.H.; Moura, R.A.; Schroeder, M.A.; Lima, A.C. Symmetry exploitation to reduce impedance evaluations in grounding grids. Int. J. Electr. Power Energy Syst. 2020, 123, 106268. [Google Scholar] [CrossRef]
  14. Guemes, J.A.; Hernando, F.E. Method for calculating the ground resistance of grounding grids using FEM. IEEE Trans. Power Deliv. 2004, 19, 595–600. [Google Scholar] [CrossRef]
  15. Hong, T.P.; Van, Q.D.; Viet, T.V. Grounding resistance calculation using FEM and reduced scale model. In Proceedings of the IEEE CEIDP, Virginia Beach, VA, USA, 18–21 October 2009; pp. 278–281. [Google Scholar]
  16. Ma, J.; Dawalibi, F.P. Analysis of grounding systems in soils with finite volumes of different resistivities. IEEE Trans. Power Deliv. 2002, 17, 596–602. [Google Scholar] [CrossRef]
  17. Li, F.; Li, Z.; Li, H.; Yang, Z.; Gao, P.; Zhou, W. Application of pseudo random correlation identification method based on LabVIEW in electromagnetic prospecting. In Proceedings of the 37th CCC, Wuhan, China, 25–27 July 2018; pp. 10289–10293. (In Chinese). [Google Scholar]
  18. Liu, Y.; Dong, H.; Liu, X.; He, Z. A preliminary study on the application of m-sequence in electrical exploration. J. Eng. Geophys. 2010, 7, 159–163. (In Chinese) [Google Scholar]
  19. Gao, Y.; Tang, R.; Liang, J.; Zhao, T.; Lu, B.; Jiang, W. Identification of transmission characteristics of airborne noise in car based on m pseudo-random sequence. J. Jilin Univ. Eng. Technol. Ed. 2012, 42, 1384–1389. (In Chinese) [Google Scholar]
  20. Song, X.; Wang, X.; Dong, Z.; Zhao, X.; Feng, X. Analysis of Pseudo-Random Sequence Correlation Identification Parameters and Anti-Noise Performance. Energies 2018, 11, 2586. [Google Scholar] [CrossRef]
  21. Cheng, L. Research on Suppression of Controllable Source Electromagnetic Sounding Interference Based on M-Sequence Correlation Identification. Master’s Thesis, China University of Geosciences, Beijing, China, 2015. (In Chinese). [Google Scholar]
  22. Liu, J. Research on Correlation Identification of M-Sequence and Its Application in Induced Polarization Method and Denoising Analysis. Master’s Thesis, China University of Geosciences, Beijing, China, 2010. (In Chinese). [Google Scholar]
  23. Wang, H. Research and Implementation of Frequency Domain Correlation Identification Based on LabVIEW. Master’s Thesis, China University of Geosciences, Beijing, China, 2013. (In Chinese). [Google Scholar]
  24. Zhi, Q.; Deng, X.; Wu, J.; Chen, X.; Zhang, J.; Wang, X.; Yang, Y. Research on the anti-jamming mechanism of electromagnetic method correlation identification and numerical simulation. J. Prog. Geophys. 2018, 33, 1755–1766. [Google Scholar]
  25. Ma, J.; Dawalibi, F.P. Extended analysis of ground impedance measurement using the fall-of-potential method. IEEE Trans. Power Deliv. 2002, 17, 881–885. [Google Scholar] [CrossRef]
  26. Korasli, C. Ground resistance measurement with alternative fall-of-potential method. IEEE Trans. Power Deliv. 2005, 20, 1657–1661. [Google Scholar] [CrossRef]
  27. Colella, P.; Pons, E.; Tommasini, R.; Silvestre, M.L.D.; Sanseverino, E.R.; Zizzo, G. Fall of Potential Measurement of the Earth Resistance in Urban Environments: Accuracy Evaluation. IEEE Trans. Ind. Appl. 2019, 55, 2337–2346. [Google Scholar] [CrossRef]
  28. Dimcev, V.; Handjiski, B.; Sekerinska, R. Alternative fall-of-potential method for grounding grids impedance measurements and inductive coupling between leads. In Proceedings of the 2003 IEEE International Symposium on Electromagnetic Compatibility, EMC ’03., Istanbul, Turkey, 11–16 May 2003; pp. 74–77. [Google Scholar]
  29. Dimcev, V.; Handjiski, B.; Vrangalov, P.; Sekerinska, R. Impedance measurement of grounding systems with alternative fall-of-potential method. In Proceedings of the Conference Record of the 2000 IEEE Industry Applications Conference, Thirty-Fifth IAS Annual Meeting and World Conference on Industrial Applications of Electrical Energy (Cat. No.00CH37129), Rome, Italy, 8–12 October 2000; pp. 3215–3219. [Google Scholar]
  30. Alcantara, F.R. Simulation of Measurements of Resistances of Grounding Systems by a Simple Hemispheric Model. In Proceedings of the 2019 IEEE XXVI International Conference on Electronics, Electrical Engineering and Computing (INTERCON), Lima, Peru, 12–14 August 2019; pp. 1–4. [Google Scholar]
  31. Alcantara, F.R. An Approximated Procedure to Find the Correct Measurement Point in the Fall-of-Potential Method. In Proceedings of the 2018 IEEE PES Transmission & Distribution Conference and Exhibition—Latin America (T&D-LA), Lima, Peru, 18–21 September 2018; pp. 1–5. [Google Scholar]
  32. Yu, Z. The application research and improvement of the fall of potential method in the grounding impedance test of the grounding grid. J. Eng. Constr. Des. 2018, 17, 84–86. (In Chinese) [Google Scholar]
  33. Tan, H.; Wu, C.; Zhang, J. Research on grounding impedance test of substation grounding grid. J. Guangxi Acad. Sci. 2020, 36, 213–217. [Google Scholar]
  34. Li, B.-N. Pseudo-Random Signal and Related Identification. Beijing, China. 1987, pp. 26–27, [Online]. Available online: http://book.sciencereading.cn/shop/book/Booksimple/onlineRead.do?id=BBB2FA6CD52C648C3B9551D48E541CE55000&readMark=1 (accessed on 8 August 2021). (In Chinese).
Figure 1. Commonly used arrangement when measuring grounding resistance using the fall-of-potential method.
Figure 1. Commonly used arrangement when measuring grounding resistance using the fall-of-potential method.
Sensors 22 05777 g001
Figure 2. The voltage electrode and current electrode are located on the same side of the grounding grid. D is the diagonal length of the grounding grid.
Figure 2. The voltage electrode and current electrode are located on the same side of the grounding grid. D is the diagonal length of the grounding grid.
Sensors 22 05777 g002
Figure 3. Principle of system identification based on the cross-correlation method.
Figure 3. Principle of system identification based on the cross-correlation method.
Sensors 22 05777 g003
Figure 4. Five-level m-sequence pseudorandom signal.
Figure 4. Five-level m-sequence pseudorandom signal.
Sensors 22 05777 g004
Figure 5. Grounding grid simulation model.
Figure 5. Grounding grid simulation model.
Sensors 22 05777 g005
Figure 6. Sine signal autocorrelation function.
Figure 6. Sine signal autocorrelation function.
Sensors 22 05777 g006
Figure 7. Pseudorandom autocorrelation function.
Figure 7. Pseudorandom autocorrelation function.
Sensors 22 05777 g007
Figure 8. Curve of grounding resistance simulation results: (a) curve of grounding resistance without interference noise; (b) curve of grounding resistance with superimposed interference noise.
Figure 8. Curve of grounding resistance simulation results: (a) curve of grounding resistance without interference noise; (b) curve of grounding resistance with superimposed interference noise.
Sensors 22 05777 g008
Figure 9. Percentage error curve of the sine signal and pseudorandom signal measuring resistance.
Figure 9. Percentage error curve of the sine signal and pseudorandom signal measuring resistance.
Sensors 22 05777 g009
Figure 10. Schematic diagram of the grounding resistance measuring device and its connection: Measuring instrument.
Figure 10. Schematic diagram of the grounding resistance measuring device and its connection: Measuring instrument.
Sensors 22 05777 g010
Figure 11. The LabVIEW interface displayed the signal waveforms and the measurement value of the grounding resistance.
Figure 11. The LabVIEW interface displayed the signal waveforms and the measurement value of the grounding resistance.
Sensors 22 05777 g011
Table 1. Standard resistance values and grounding resistance values measured by two excitation source signals.
Table 1. Standard resistance values and grounding resistance values measured by two excitation source signals.
Standard Resistance Value (Ω)Resistance Value Measured by Sinusoidal Signal (Ω)Error Percentage (%)Resistance Value Measured by M-Sequence Pseudorandom Signal (Ω)Error Percentage (%)
Group 1Group 2Group 3Average ValueGroup 1Group 2Group 3Average Value
0.6560.670.670.680.672.640.670.650.650.660.10
1.6981.641.661.651.652.831.691.71.711.70.12
2.2731.972.452.332.251.012.272.272.282.270.01
3.5473.483.63.473.520.863.533.543.543.540.29
4.5884.674.664.654.661.574.574.594.574.580.25
4.7544.944.824.694.821.324.784.744.754.760.06
5.3415.325.375.395.360.365.345.345.345.340.02
6.2176.216.166.196.190.496.216.226.216.210.06
6.5436.446.466.66.50.666.526.566.556.540.01
8.4218.058.378.58.311.368.48.418.428.410.13
9.2349.439.119.149.230.089.249.259.249.240.10
10.56810.6510.6510.6910.660.9010.5510.5610.5710.560.08
Table 2. The grounding resistance measurement value and the error between the measurement value and the actual value.
Table 2. The grounding resistance measurement value and the error between the measurement value and the actual value.
The Actual Resistance Value (Ω)The Measurement Resistance Value (Ω)Error Percentage (%)
0.30.3062.00
0.50.5132.60
1.01.0131.30
1.51.4423.87
1.81.7651.94
21.9383.10
2.32.2422.52
2.52.4760.96
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Duan, Q.; Zou, B.; Li, S.; Zhao, H. Simulation and Experiment of Substation Grounding Resistance Measurement Based on Pseudorandom Signal. Sensors 2022, 22, 5777. https://doi.org/10.3390/s22155777

AMA Style

Duan Q, Zou B, Li S, Zhao H. Simulation and Experiment of Substation Grounding Resistance Measurement Based on Pseudorandom Signal. Sensors. 2022; 22(15):5777. https://doi.org/10.3390/s22155777

Chicago/Turabian Style

Duan, Qingming, Bofeng Zou, Shanshan Li, and Hu Zhao. 2022. "Simulation and Experiment of Substation Grounding Resistance Measurement Based on Pseudorandom Signal" Sensors 22, no. 15: 5777. https://doi.org/10.3390/s22155777

APA Style

Duan, Q., Zou, B., Li, S., & Zhao, H. (2022). Simulation and Experiment of Substation Grounding Resistance Measurement Based on Pseudorandom Signal. Sensors, 22(15), 5777. https://doi.org/10.3390/s22155777

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop