Electromagnetic Environment of Transmission Line Based on Full Parameter Online Estimation

Zidan Sun* , Xiaofeng Zhou** , Likai Liang* and Yang Mo*

Abstract

Abstract: The parameters of transmission lines have an influence on the electromagnetic environment surrounding the line. This paper proposes a method based on phasor measurement unit (PMU) and supervisory control and data acquisition (SCADA) to achieve online estimation of transmission line full parameters, such as resistance, reactance and susceptance. The proposed full parameter estimation method is compared with the traditional method of estimating resistance independently based on SCADA system. Then, the electromagnetic environ-ment is analyzed based on the different parameter estimation methods. The example results illustrate that online estimation of transmission line full parameters is more accurate in the analysis of electromagnetic environment, which further confirms its necessity and significance in engineering application.

Keywords: Current Capacity , Electromagnetic Environment , Line Parameter , Power Grid , Transmission Line

1. Introduction

With the development of society, the scale expanding of power grid becomes the general trend [1]. The transmission lines are getting closer to public living areas, resulting in many electromagnetic environmental problems [2]. At the same time, public awareness of environmental protection is increasing, and complaints caused by electromagnetic environment often occur, leading to local power grid construction lagging behind the demand for electricity [3]. It not only affects the residential electricity, but also restricts the development of local economy.

The electromagnetic environment of alternating current transmission lines includes power frequency electric field, power frequency magnetic field, radio interference and audible noise [4]. Power frequency magnetic field has become one of the key issues of current research because of its highest public concern and numerous disputes [5,6]. The system power flow includes the line current and the line voltage [7]. The variation of current will impact the magnetic field of transmission lines, and the variation of voltage will impact the electric field [8]. The accurate calculation of current and voltage is beneficial to the exploration of electromagnetic fields [9]. A lot of literatures analyze the electromagnetic environment of high-voltage transmission lines based on the theoretical parameter values [10,11]. In the actual operation of the line, there is a certain gap between the theoretical value and the actual value of the parameters such as the resistance and the reactance, which will affect the system power flow, such as the voltages and currents of all transmission lines [12].

The actual parameters deviate from the design values due to the long-term operation of the lines [13,14]. Online tracking and prediction of transmission line resistance can be achieved by phasor measurement unit (PMU) [15]. PMU is characterized by the short sampling period, the measurable voltage and the current phase, which can quickly evaluate the performance of power system and achieve rapid response under the unexpected circumstance [16]. However, it is impractical to configure PMU on all the transmission lines [17]; hence, supervisory control and data acquisition (SCADA) systems are applied to track transmission line resistance [18]. SCADA system can provide abundant real-time measurements of the voltage, the current, the active power and the reactive power. Online estimation of transmission line parameters can be achieved. Although the author [18] achieves the tracking and prediction of transmission line resistance, the theoretical value of other parameters is still adopted. To some extent, it affects the accuracy of online tracking results so that the electromagnetic environment surrounding the line can't be calculated accurately. During the long-time operation, the parameters may change because of the aging of conductors or the influence of climate, geography and environment [19,20].

In this paper, considering the above problems, the method of online estimation of transmission line parameters, and the analysis method of parameters impacting the electromagnetic environment are presented. The rest of this paper is organized as follows. Section 2 introduces the online estimation of transmission line parameters based on PMU and SCADA system, and the calculation of electromagnetic environment. The online estimation results of parameters on the proposed method are presented in Section 3. The method of full-parameter online estimation is compared with the traditional method of estimating resistance independently. Based on the online estimation results, the electromagnetic environment of the IEEE 5-bus system under three cases are also analyzed and compared in Section 3. Conclusions are given finally in Section 4.

2. Method and Model

2.1 Online Estimation of Transmission Line Parameters Based on PMU

Online tracking of the full parameters, such as the resistance, reactance and admittance can be achieved by PMU system. The model of [TeX:] $$\Pi -type$$ equivalent line is shown in Fig. 1. [TeX:] $$P_{s}$$ and s are the active power and reactive power at the sending end. [TeX:] $$P_{r}$$ and r are the active power and reactive power at the receiving end. [TeX:] $$R_{l} \text { and } X_{l}$$ is the resistance and the reactance of the transmission line. [TeX:] $$\mathbf{I}_{\mathbf{s}}$$ is the injection current and [TeX:] $$\mathrm{U}_{\mathrm{s}}$$ is the bus voltage at the sending end. [TeX:] $$I_{r}$$ is the injection current and [TeX:] $$U_{r}$$ is the bus voltage at the receiving end. [TeX:] $$\dot{I}$$ is the current carrying of transmission line. [TeX:] $$B_{s} \text { and } B_{r}$$ are the susceptance at the sending end and receiving end respectively.

Ignoring the difference of the susceptance between the sending and receiving end, such as [TeX:] $$B_{s}=B_{r}=B_{l}.$$ According to two-port circuit theory of [TeX:] $$II -type$$ equivalent circuit, Eq. (1) is established.

(1)
[TeX:] $$\begin{array}{l} \dot{I}_{s}=\left(\dot{U}_{s}-\dot{U}_{r}\right) /\left(R_{l}+j X_{l}\right)+\dot{U}_{s}\left(j B_{l}\right) \\ \dot{I}_{r}=\left(\dot{U}_{s}-\dot{U}_{r}\right) /\left(R_{l}+j X_{l}\right)-\dot{U}_{r}\left(j B_{l}\right) \end{array}$$

Fig. 1.
The [TeX:] $$\Pi- type$$ equivalent line.

Matrix form of Eq. (1) is shown in Eq. (2).

(2)
[TeX:] $$\mathrm{I}=\mathrm{UY}$$

(3)
[TeX:] $$I=\left[\begin{array}{l} \dot{I}_{s} \\ \dot{I}_{r} \end{array}\right] \mathrm{U}=\left[\begin{array}{cc} \dot{U}_{s}-\dot{U}_{r} & \dot{U}_{s} \\ \dot{U}_{s}-\dot{U}_{r} & -\dot{U}_{r} \end{array}\right] \quad \mathrm{Y}=\left[\begin{array}{c} 1 /\left(R_{l}+j X_{l}\right) \\ j B_{l} \end{array}\right]$$

When the measured data [TeX:] $$m(m>2),$$ there are overdetermined equations. EM is set as the error phasor, Eq. (4) is established.

(4)
[TeX:] $$\mathrm{E}_{M}=\mathrm{I}_{M}-\mathrm{U}_{M} \mathrm{Y}_{M}$$

where

(5)
[TeX:] $$\mathrm{E}_{\mathrm{M}}=\left[\begin{array}{llll} \mathrm{e}_{1} & \mathrm{e}_{2} & \cdots & \mathrm{e}_{\mathrm{m}} \end{array}\right]^{\mathrm{T}}$$

(6)
[TeX:] $$\mathrm{I}_{\mathrm{M}}=\left[\mathrm{I}_{1}^{\mathrm{T}} \quad \mathrm{I}_{2}^{\mathrm{T}} \quad \cdots \quad \mathrm{I}_{\mathrm{m}}^{\mathrm{T}}\right]^{\mathrm{T}}$$

(7)
[TeX:] $$\mathrm{U}_{\mathrm{M}}=\left[\mathrm{U}_{1}^{\mathrm{T}} \quad \mathrm{U}_{2}^{\mathrm{T}} \quad \cdots \quad \mathrm{U}_{\mathrm{m}}^{\mathrm{T}}\right]^{\mathrm{T}}$$

(8)
[TeX:] $$Y_{M}=Y$$

The square sum of the error component is set as J, which is used as the performance index.

(9)
[TeX:] $$\min \mathrm{J}=\sum_{i=1}^{m} e_{i}^{2}=E_{M}^{T} E_{M}=\left(I_{M}-U_{M} Y_{M}\right)^{T}\left(I_{M}-U_{M} Y_{M}\right)$$

When the minimum value of J is taken as a necessary condition, Eq. (10) is obtained.

(10)
[TeX:] $$\left.\frac{\partial \mathrm{J}}{\partial \mathrm{Y}_{M}}\right|_{\mathrm{Y}_{\mathrm{M}}=\hat{\mathrm{Y}}_{\mathrm{M}}}= -2 \mathrm{U}_{M}^{T} \mathrm{I}_{M}+2 \mathrm{U}_{M}^{T} \mathrm{U}_{M} \hat{\mathrm{Y}}_{M}$$

Therefore,

(11)
[TeX:] $$\hat{\mathrm{Y}}_{M}=\left(\mathrm{U}_{M}^{T} \mathrm{U}_{M}\right)^{-1} \mathrm{U}_{M}^{T} \mathrm{I}_{M}$$

Online estimation of [TeX:] $$\Pi -type$$ equivalent line can be achieved by Eq. (11). The [TeX:] $$R_{l}, X_{l}, B_{s} \text { and } B_{r}$$ included in the extended state variables of the transmission lines can be obtained based on PMU. Due to the high cost of PMU devices and the richness of SCADA data in power system, SCADA systems are put into use.

2.2 Online Estimation of Transmission Line Parameters Based on PMU

In practice, the active power [TeX:] $$\left(P_{s}, P_{r}\right),$$ the reactive power [TeX:] $$\left(Q_{s}, Q_{r}\right),$$ the current value[TeX:] $$\left(I_{s}, I_{r}\right)$$ and the voltage [TeX:] $$\left(U_{s}, U_{r}\right)$$ can be obtained by SCADA. The phase angle of reference node is set as phase angle [TeX:] $$\dot{I}$$ and [TeX:] $$\dot{I}=I \angle 0^{\circ},$$ the estimated state variables can be expressed as Eq. (12).

(12)
[TeX:] $$\mathbf{x}=\left[\begin{array}{llllllll} U_{s} & U_{r} & \theta_{s} & \theta_{r} & I & R_{l} & X_{l} & B_{l} \end{array}\right]^{T}$$

where [TeX:] $$\theta_{s} \text { and } \theta_{r}$$ are the voltage phase angles at both ends of the transmission lines. According to the principle of electric engineering, Eq. (13) is established [14].

(13)
[TeX:] $$\begin{aligned} &I_{s}^{2}=I^{2}+B_{l}^{2} U_{s}^{2}-2 I B_{l} U_{s} \sin \theta_{s}\\ &I_{r}^{2}=I^{2}+B_{l}^{2} U_{r}^{2}+2 I B_{l} U_{r} \sin \theta_{r} \end{aligned} \\ \begin{aligned} &0=U_{s}^{2}+U_{r}^{2}-2 U_{s} U_{r} \cos \left(\theta_{s}-\theta_{r}\right)-\left(R_{l}^{2}+X_{l}^{2}\right) I^{2}\\ &0=\frac{P_{r}^{2}+\left(Q_{r}-B_{l} U_{r}^{2}\right)^{2}}{U_{r}^{2}}-\frac{P_{s}^{2}+\left(Q_{s}-B_{l} U_{s}^{2}\right)^{2}}{U_{s}^{2}} \end{aligned} \\ \begin{aligned} &\frac{P_{s}^{2}+\left(Q_{s}-B_{l} U_{s}^{2}\right)^{2}}{U_{s}^{2}}=I^{2}\\ &P_{l a s s}=I^{2} R_{l} \end{aligned} \\ Q_{l o s s}=I^{2} X_{l}-B_{l}\left(U_{s}^{2}+U_{r}^{2}\right)$$

where [TeX:] $$P_{\text {loss }}$$ [TeX:] $$P_{\text {loss}}=P_{s}-P_{r} ; Q_{\text {loss}}$$ is the reactive power loss, and [TeX:] $$Q_{\text {loss}}=Q_{s}-Q_{r}.$$ From a mathematical point of view, Eq. (14) is shown as the measurement equation based on the least square estimation of a time section.

(14)
[TeX:] $$\mathbf{z}=\mathbf{h}(\mathbf{x})+\mathbf{v}$$

where z is the measurement vector, [TeX:] $$\mathbf{h}(\cdot)$$ is the measurement function vector, v is the disturbance vector of the random measurement error. According to each measurement data based on SCADA, such as [TeX:] $$\mathrm{P}_{\mathrm{s}},\mathrm{P}_{\mathrm{r}}, \mathrm{Q}_{\mathrm{s}}, \mathrm{Q}_{\mathrm{r}}, \mathrm{I}_{\mathrm{s}}, \mathrm{I}_{\mathrm{r}}, \mathrm{U}_{\mathrm{s}}, \mathrm{U}_{\mathrm{r}}.$$ Measurement vector based on the least square estimation of a time section is shown in Eq. (15), and the measurement function vector is shown in Eq. (16).

(15)
[TeX:] $$\mathbf{z}=\left[\begin{array}{lllllllll} U_{s} & U_{r} & I_{s}^{2} & I_{r}^{2} & 0 & P_{s}^{2} & P_{r}^{2} & P_{\text {lass}} & Q_{\text {loss}} \end{array}\right]^{\mathrm{T}}$$

(16)
[TeX:] $$\mathbf{h}(\mathbf{x})=\left[\begin{array}{l} U_{s} \\ U_{r} \\ I^{2}+B_{l}^{2} U_{s}^{2}-2 I B_{l} U_{s} \sin \theta_{s} \\ I^{2}+B_{l}^{2} U_{r}^{2}+2 I B_{l} U_{r} \sin \theta_{r} \\ U_{s}^{2}+U_{r}^{2}-2 U_{s} U_{r} \cos \left(\theta_{s}-\theta_{r}\right)-\left(R_{l}^{2}+X_{l}^{2}\right) I^{2} \\ I^{2} U_{r}^{2}-\left(Q_{r}-B_{l} U_{r}^{2}\right)^{2} \\ I^{2} U_{s}^{2}-\left(Q_{s}-B_{l} U_{s}^{2}\right)^{2} \\ I^{2} R_{l} \\ I^{2} X_{l}-B_{l}\left(U_{s}^{2}+U_{r}^{2}\right) \end{array}\right]$$

According to Eq. (14), the modified equation of least square estimation is shown in Eq. (17) and its iterative equation is shown in Eq. (18).

(17)
[TeX:] $$\Delta \hat{\mathbf{x}}^{(k)}=\left[\mathbf{H}^{\mathrm{T}}\left(\hat{\mathbf{x}}^{(k)}\right) \mathbf{R}^{-1} \mathbf{H}\left(\hat{\mathbf{x}}^{(k)}\right)\right]^{-1} \mathbf{H}^{\mathrm{T}}\left(\hat{\mathbf{x}}^{(k)}\right) \mathbf{R}^{-1}\left[\mathbf{z}-\mathbf{h}\left(\hat{\mathbf{x}}^{(k)}\right)\right]$$

(18)
[TeX:] $$\hat{\mathbf{x}}^{(k+1)}=\hat{\mathbf{x}}^{(k)}+\Delta \hat{\mathbf{x}}^{(k)}$$

The superscript (k) represents the k-th iteration. The corresponding Jacobi matrix is shown in Eq. (19).

(19)
[TeX:] $$\mathbf{H}\left(\hat{\mathbf{x}}^{(k)}\right)=\left.\frac{\partial \mathbf{z}}{\partial \mathbf{x}}\right|_{\mathbf{x}=\hat{\mathbf{x}}^{(k)}}$$

According to Eq. (14) and Eq. (18), online estimation of the parameters such as [TeX:] $$\mathrm{R}_{1}, \mathrm{X}_{1} \text { and } \mathrm{B}_{1}$$ included in the extended state variables of transmission lines can be obtained, which is known as online estimation of full parameters based on SCADA. Therefore, online estimation of the resistance, reactance and susceptance can be achieved. Due to the universality and the richness of SCADA data in power system, full parameter estimation based on SCADA system can be obtained without extra hardware facilities and provide a reliable basis for power grid analysis.

2.3 Calculation of Electromagnetic Environment

The electromagnetic environment of alternating current transmission lines includes power frequency electric field, power frequency magnetic field, radio interference and audible noise, of which power frequency magnetic field is the main research direction. The calculation of electric field intensity at point is shown in Eqs. (20) and (21).

(20)
[TeX:] $$\left[Q_{c}\right]=[C]^{-1}[U]$$

(21)
[TeX:] $$\begin{aligned} &E_{x}=\frac{1}{2 \pi \varepsilon} \sum_{k=1}^{n} Q_{c k}\left(\frac{x-x_{k}}{d_{k}^{2}}-\frac{x-x_{k}}{d_{k}^{\prime 2}}\right)\\ &E_{y}=\frac{1}{2 \pi \varepsilon} \sum_{k=1}^{n} Q_{c k}\left(\frac{y-y_{k}}{d_{k}^{2}}-\frac{y-y_{k}}{d_{k}^{\prime 2}}\right) \end{aligned}$$

where [c], [C] and [U] are the charge, potential coefficient matrix and voltage matrix, respectively. [TeX:] $$E_{x}$$ and [TeX:] $$E_{y}$$ are the abscissa and ordinate components of electric field intensity, respectively. [TeX:] $$\varepsilon$$ is the air dielectric constant. n is the sum of the number of lines and mirror lines. k is the line number; [TeX:] $$x_{k}$$ and [TeX:] $$y_{k}$$ are the abscissa and ordinate of the line . The distance from the point A(x, y) to the line k and its mirror line are demoted by [TeX:] $$d_{k} \text { and } d'_{k},$$ respectively. The potential coefficients are calculated by Eq. (22).

(22)
[TeX:] $$\begin{aligned} C_{m, m} &=\frac{1}{2 \pi \varepsilon} \ln \frac{2 H}{R} \\ C_{m, n} &=\frac{1}{2 \pi \varepsilon} \ln \frac{D_{m, n}}{D_{m, n}} \\ C_{m, n} &=C_{n, m} \end{aligned}$$

where R is the line radius. [TeX:] $$D_{m n} \text { and } D'_{m, n}$$ are the distance from the line m to the line m and its mirror line, respectively. H is the vertical distance from the line to ground. According to [TeX:] $$E=\sqrt{E_{x}^{2}+E_{y}^{2}},$$ electric field intensity can be obtained. It can be seen that the variation of the line voltage will directly affect the power frequency electric field.

According to the ampere circuit theorem, the magnetic field intensity at point A(x, y) can be calculated by Eq. (23):

(23)
[TeX:] $$\begin{array}{l} M_{x}=\sum_{k=1}^{n} \frac{\mu I_{k}}{2 \pi}\left(\frac{x_{k}-x}{d_{k}^{2}}\right) \\ M_{y}=\sum_{k=1}^{n} \frac{\mu I_{k}}{2 \pi}\left(\frac{y_{k}-y}{d_{k}^{2}}\right) \end{array}$$

where [TeX:] $$M_{x} \text { and } M_{y}$$ are the abscissa and ordinate components of magnetic field intensity, respectively. [TeX:] $$\mu$$ is the magnetic permeability. [TeX:] $$I_{k}$$ is the current of the line k. According to [TeX:] $$M=\sqrt{M_{x}^{2}+M_{y}^{2}}$$ , magnetic field intensity can be calculated. The power frequency magnetic field will be directly affected by the variation of the line current.

3. Case Study

3.1 Analysis of 500kV Transmission Line Parameters Based on PMU

In this paper, a case study of a 500kV transmission line equipped with PMU is conducted. The PMU sampling width is 35s and the PMU sampling interval is 5s. The theoretical value of parameters is [TeX:] $$\mathrm{R}=6.0951 \Omega, \mathrm{X}=57.4038 \Omega, \mathrm{B}_{1}=0.000478 \mathrm{S}.$$ According to the [TeX:] $$\Pi -type$$ equivalent line based on PMU, online estimation of R, X and Bl at sampling interval points can be calculated and the results are shown in Table 1 and Fig. 2.

According to the estimation results of transmission line parameters in Table 1 and Fig. 2, there is a gap between the actual value and the theoretical value. During the actual operation of power grid, it may lead to the inaccurate state estimation. This further verifies the necessity and the importance of online estimation of transmission line parameters.

Table 1.
Online estimation results of line parameters based on PMU
Fig. 2.
Online estimation results of line parameters based on PMU: (a) resistance, (b) reactance, and (c) susceptance.
3.2 220 kV Transmission Line Parameters Based on SCADA

A case study of a 220kV transmission line named Yuxin Line from Gangyu Station to Xingang Station in Yantai is conducted. The parameter information of the transmission line is shown in Table 2. The theoretical value of parameters is [TeX:] $$\mathrm{R}=3.95 \Omega, \mathrm{X}=20.8 \Omega, \mathrm{B}_{1}=68.325 \mu \mathrm{S},$$ which is set as the base case. The test time is from 21:00 to 22:00. The sampling interval is 5 minutes and there are 12 groups of samples.

Table 2.
The information of line parameters
Table 3.
Comparisons of online estimation results based on SCADA

Firstly, online estimation of full parameters is achieved by the proposed method in this paper. The result of full-parameter estimation is shown in the left part of Table 3, and there are 12 sampling points. The maximum difference between the actual value and the theoretical value of the resistance is 0.5055 at the 11th sampling point, and the maximum difference between the actual value and the theoretical value of the reactance is 2.8242 at the 8th sampling point.

Then, the resistance is estimated independently online while the reactance and the susceptance are set to be the theoretical value. The results of independent resistance estimation based on SCADA are shown in the right part of Table 3. Fig. 3 further describes online estimation results of the resistance under the two methods. It can also be seen from Table 3 and Fig. 3 that the maximum difference of resistance between full-parameter estimation and independent resistance estimation is 0.206 at the 2nd sampling point, which accounts for 4.8% of the total resistance. The difference does exist and it is necessary to be considered. This further illustrates the necessity of the full-parameter online estimation. The full parameter estimation of the resistance will be more appropriate with the actual operation of transmission lines.

Fig. 3.
Online estimation results of the resistance under two methods.
3.3 Analysis of Electromagnetic Environment

A 5-bus power system is taken as the example. The transmission line type is LGJ-400/50, and the voltage level is 220 kV. The conductor diameter is 27.63 mm. The phase sequence is horizontal arrangement, and the spacing is 8 m. The height of the line to ground is 10 m. The lengths of branch 1– 2, 2–3, and 1–3 are all 50 km. The electromagnetic environment indicator A is the midpoint of branch 2– 3. The network structure is shown in Fig. 4.

Based on the online estimation results of transmission line parameters in Section 3.2, the electromagnetic environment of alternating current transmission lines is obtained. The theoretical value of parameters is set as the base case. The result of full-parameter estimation at the 11th sampling point is set as the case 1, and the result of independent resistance estimation at the 6th sampling point is set as the case 2. The line parameters of three cases are shown in Table 4.

Fig. 4.
The network structure diagram.
Table 4.
Line parameters of three cases

The results of voltage and current are obtained, and then the electromagnetic environment surrounding the line is analyzed. The voltage magnitudes at point A under the three cases are all 1.04 p.u., so that the electric field intensity of three cases are close. For the currents at point A under the three cases, base case is 43.94 A, case 1 is 44.58 A, and case 2 is 44.55 A. The magnetic field intensity of three cases are different. The height of calculation point to ground is 1.5 m. The horizontal range is between -20 m and +20 m. The magnetic field intensity is calculated every 0.01 m horizontal distance. The distribution curves of the magnetic field intensity distribution at point A is shown in Fig. 5.

Fig. 5.
Power frequency magnetic field intensity at point A.

The curves are distributed unimodally, presenting lateral attenuation trend with horizontal distance increasing. The maximum magnetic intensity [TeX:] $$\mathrm{M}_{\mathrm{m}}$$ occurs at the center of the tower. Overall, for the magnetic field intensity, case 1 and case 2 are both higher than base case. The closer to the tower, the more obvious the phenomenon. The magnetic field intensity at the [TeX:] $$\mathrm{M}_{\mathrm{m}}$$ of three cases are shown in Table 5. Case 1 is [TeX:] $$1.7200 \mu \mathrm{T},$$ which is 1.46% higher than base case. Case 2 is [TeX:] $$1.7188 \mu \mathrm{T},$$ which is 1.39% higher than base case. It can be seen that online estimation of parameters have a great influence on the power frequency magnetic field.

Table 5.
The maximum magnetic field intensity of three cases

4. Conclusion

In this paper, online estimation of transmission line parameters is achieved based on PMU. On the basis of SCADA, the parameters of any transmission line in power grid can be estimated without adding extra hardware by the least square method. The electromagnetic environment surrounding the line is analyzed and compared based on the variation of line voltage and current. The actual results show that line parameters have a great influence on the electromagnetic environment. The higher accuracy of online estimation of transmission line parameters, the higher utilization and reliability of power grid. The proposed method can be adapted to any transmission line.

Acknowledgement

This paper is supported by the Science and Technology Development Project of Weihai City (No. 2019DXGJ23).

Biography

Zidan Sun
https://orcid.org/0000-0002-1502-5513

He was born in Jiangsu province, China, in 1994. He received his B.E. degree in School of Software from Shandong University in 2017. He is currently pursuing his M.E. degree in electronics and communication engineering at Shandong University, China. His main research interests include power system operation and control.

Biography

Xiaofeng Zhou
https://orcid.org/0000-0002-3545-8738

He was born in Jiangsu province, China, in 1994. He received his B.E. degree in School of Software from Shandong University in 2017. He is currently pursuing his M.E. degree in electronics and communication engineering at Shandong University, China. His main research interests include power system operation and control. He obtained Bachelor’s degree in Agricultural Machinery Engineering from Shandong University of Technology. He received his master’s degree from Shandong University of Technology and now teaches at the Department of Electrical and Mechanical Engineering in Weihai Vocational College, China. The current research fields are numerical control technology and mechatronics technology.

Biography

Likai Liang
https://orcid.org/0000-0003-1271-6906

She received Ph.D. degree in School of Electrical Engineering from Shandong University in 2013. She teaches as an associate professor in School of Mechanical Electrical and Information Engineering of Shandong University. Her current research interests include power system operation and control.

Biography

Yang Mo
https://orcid.org/0000-0001-6734-8691

He was born in Shandong Province, in China, on November 21, 1993. He obtained Bachelor in Electrical Engineering and Automation from Qingdao University in 2012. Now, he is a postgraduate in School of Mechanical Electrical and Information Engi-neering at Shandong University, China. He is major in electronics and communication engineering. His main research interests include power system operation and control.

References

  • 1 A. Costa, D. Georgiadis, T. S. Ng, M. Sim, "An optimization model for power grid fortification to maximize attack immunity," International Journal of Electrical Power & Energy Systems, vol. 99, pp. 594-602, 2018.custom:[[[-]]]
  • 2 I. S. Okrainskaya, A. I. Sidorov, S. P. Gladyshev, "Electromagnetic environment under over head power transmission lines 110–500 kV," in Proceedings of International Symposium on Power Electronics Power Electronics, Electrical Drives, Automation and Motion, Sorrento, Italy, 2012;pp. 796-801. custom:[[[-]]]
  • 3 E. Al-Bassam, A. Elumalai, A. Khan, L. Al-Awadi, "Assessment of electromagnetic field levels from surrounding high-tension overhead power lines for proposed land use," Environmental Monitoring and Assessment, vol. 188, no. 316, 2016.custom:[[[-]]]
  • 4 J. He, S. Chen, J. Guo, R. Zeng, J. Lee, S. Chang, B. Zhang, J. Zou, Z. Guan, "Electromagnetic environment analysis of a software park near transmission lines," IEEE Transactions on Industry Applications, vol. 40, no. 4, pp. 995-1002, 2004.custom:[[[-]]]
  • 5 M. Repacholi, "Concern that "EMF" magnetic fields from power lines cause cancer," Science of the Total Environment, vol. 426, pp. 454-458, 2012.custom:[[[-]]]
  • 6 I. N. Ztoupis, I. F. Gonos, I. A. Stathopulos, "Uncertainty evaluation in the measurement of power frequency electric and magnetic fields from AC overhead power lines," Radiation Protection Dosimetry, vol. 157, no. 1, pp. 11-21, 2013.custom:[[[-]]]
  • 7 Y. Mo, Y. Wang, F. Song, Z. Xu, Q. Zhang, Z. Niu, "Investigating the impacts of meteorological parameters on electromagnetic environment of overhead transmission line," Progress In Electromagnetics Research, vol. 70, pp. 177-185, 2018.custom:[[[-]]]
  • 8 Y. Wang, H. Wang, H. Xue, C. Yang, T. Yan, "Research on the electromagnetic environment of 110kV six-circuit transmission line on the same tower," in Proceedings of IEEE PES Innovative Smart Grid Technologies, Tianjin, China, 2012;pp. 1-5. custom:[[[-]]]
  • 9 M. Sakashita, K. Nishi, S. Ito, T. Mifune, T. Matsuo, "Postcorrection of current/voltage and electromagnetic force for efficient hysteretic magnetic field analysis," IEEE Transactions on Magnetics, vol. 53, no. 6, pp. 1-4, 2017.custom:[[[-]]]
  • 10 L. Zhao, J. Lu, G. Wu, "Measurement and analysis on electromagnetic environment of 1000kV UHV AC transmission line," in Proceedings of 2012 Asia-Pacific Power and Energy Engineering Conference, Shanghai, China, 2012;pp. 1-4. custom:[[[-]]]
  • 11 M. Sibanda, R. R. Van Zyl, N. Parus, "Overview of the electromagnetic environment in the vicinity of HVDC transmission lines," in Proceedings of 2013 Proceedings of the 10th Industrial and Commercial Use of Energy Conference, Cape Town, South Africa, 2013;pp. 1-7. custom:[[[-]]]
  • 12 Z. Siroma, T. Ioroi, "Expected electrochemical impedance responses of porous electrodes based on theoretical solutions of transmission-line models," Electrochemistry, vol. 83, no. 6, pp. 425-433, 2015.custom:[[[-]]]
  • 13 F. Chen, X. Han, M. Li, M. Wang, M. Yang, "Tracking estimation of transmission line temperature based on PMU measurement," Automation of Electric Power Systems, vol. 33, no. 19, pp. 25-29, 2009.custom:[[[-]]]
  • 14 F. Chen, X. Han, K. Kang, H. Li, "Tracking of dynamic thermal rating of transmission line based on SCADA," Automation of Electric Power Systems, vol. 34, no. 5, pp. 81-85, 2010.custom:[[[-]]]
  • 15 J. C. Ding, Z. X. Cai, K. Y. Wang, "An overview of state estimation based on wide-area measurement system," Automation of Electric Power Systems, vol. 30, no. 7, pp. 98-103, 2006.custom:[[[-]]]
  • 16 Y. Xue, W. Xu, Z. Dong, Q. Wan, "A review of wide area measurement system and wide area control system," Automation of Electric Power Systems, vol. 31, no. 15, pp. 1-5, 2007.custom:[[[-]]]
  • 17 Y. Tohidi, L. Olmos, M. Rivier, M. R. Hesamzadeh, "Coordination of generation and transmission development through generation transmission charges: a game theoretical approach," IEEE Transactions on Power Systems, vol. 32, no. 2, pp. 1103-1114, 2016.custom:[[[-]]]
  • 18 C. Fang, "Studies on theory and expansion of power network state estimation," PhD dissertationShandong University, Jinan, China, 2010.custom:[[[-]]]
  • 19 K. Kopsidas, A. Kapetanaki, V. Levi, "Optimal demand response scheduling with real-time thermal ratings of overhead lines for improved network reliability," IEEE Transactions on Smart Grid, vol. 8, no. 6, pp. 2813-2825, 2016.doi:[[[10.1109/TSG.2016.2542922]]]
  • 20 S. Ge, J. Li, T. Li, H. Liu, R. Li, "Integrated analysis on reliability of power distribution network and urban road network," in Proceedings of the CSEE, 2016;vol. 36, no. 6, pp. 1568-1577. custom:[[[-]]]

Table 1.

Online estimation results of line parameters based on PMU
Time (s) Resistance [TeX:] $$(\Omega)$$ Reactance[TeX:] $$(\Omega)$$ Susceptance[TeX:] $$\left(\times 10^{-3} \mathrm{S}\right)$$
0 6.0937 49.0941 0.4512
5 6.0874 49.0613 0.4509
10 6.1058 49.0173 0.4510
15 6.1163 48.9674 0.4508
20 6.1400 49.0084 0.4510
25 6.1338 48.9493 0.4510
30 6.1490 49.0507 0.4514
35 6.1061 49.1011 0.4513

Table 2.

The information of line parameters
Line parameters Information
Line name Yuxin Line
The first station Gangyu Station
The last station Xingang Station
Voltage level (kV) 220
Line type LGJ-400
Line length (km) 50
Unit resistance [TeX:] $$(\Omega)$$ 0.079
Unit reactance [TeX:] $$(\Omega)$$ 0.416
Unit susceptance [TeX:] $$(\mu \mathrm{S})$$ 2.733
Unit specific heat [TeX:] $$\left(\mathrm{J} /\left(\mathrm{km} \cdot^{\circ} \mathrm{C}\right)\right)$$ 1127
Safety current (A) 845

Table 3.

Comparisons of online estimation results based on SCADA
Sampling point Full-parameter online estimation Independent resistance online estimation
Resistance [TeX:] $$(\Omega)$$ Reactance[TeX:] $$(\Omega)$$ Susceptance[TeX:] $$(\mu \mathrm{S})$$ Resistance[TeX:] $$(\Omega)$$ Reactance[TeX:] $$(\Omega)$$ Susceptance[TeX:] $$(\mu \mathrm{S})$$
1 4.3729 19.3273 68.431 4.2890 20.8000 68.325
2 4.2335 18.8677 68.376 4.4395 20.8000 68.325
3 4.4231 19.3496 68.302 4.3859 20.8000 68.325
4 4.3804 18.5730 68.412 4.4086 20.8000 68.325
5 4.2799 19.9671 68.297 4.3275 20.8000 68.325
6 4.4238 18.7734 68.288 4.5064 20.8000 68.325
7 4.2803 18.5423 68.305 4.2076 20.8000 68.325
8 4.2640 17.9758 68.349 4.1098 20.8000 68.325
9 4.3346 17.9857 68.286 4.4683 20.8000 68.325
10 4.3709 19.3020 68.376 4.2279 20.8000 68.325
11 4.4555 19.4563 68.478 4.4203 20.8000 68.325
12 4.3644 19.8603 68.455 4.3068 20.8000 68.325

Table 4.

Line parameters of three cases
Resistance () Reactance () Susceptance [TeX:] $$(\mu S)$$ Voltage [TeX:] $$(k V)$$ Current (A)
Base case 3.9500 20.8000 68.325 228.8 43.94
Case 1 4.4555 19.4563 68.478 228.8 44.58
Case 2 4.5064 20.8000 68.325 228.8 44.55

Table 5.

The maximum magnetic field intensity of three cases
[TeX:] $$\mathbf{M}_{\mathrm{m}}(\mu \mathrm{T})$$ Difference (%)
Base case 1.6953 0
Case 1 1.7200 1.46
Case 2 1.7188 1.39
The [TeX:] $$\Pi- type$$ equivalent line.
Online estimation results of line parameters based on PMU: (a) resistance, (b) reactance, and (c) susceptance.
Online estimation results of the resistance under two methods.
The network structure diagram.
Power frequency magnetic field intensity at point A.