Zidan Sun* , Xiaofeng Zhou** , Likai Liang* and Yang Mo*Electromagnetic Environment of Transmission Line Based on Full Parameter Online EstimationAbstract: 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. IntroductionWith 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 Model2.1 Online Estimation of Transmission Line Parameters Based on PMUOnline 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}$$Matrix form of Eq. (1) is shown in Eq. (2).
(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.
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}}$$
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 PMUIn 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.
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]$$
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 EnvironmentThe 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).
(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 Study3.1 Analysis of 500kV Transmission Line Parameters Based on PMUIn 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.
3.2 220 kV Transmission Line Parameters Based on SCADAA 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.
Table 3.
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. 3.3 Analysis of Electromagnetic EnvironmentA 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. Table 4.
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. 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. 4. ConclusionIn 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. BiographyZidan Sunhttps://orcid.org/0000-0002-1502-5513He 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. BiographyXiaofeng Zhouhttps://orcid.org/0000-0002-3545-8738He 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. BiographyYang Mohttps://orcid.org/0000-0001-6734-8691He 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
|