Prediction of radio interference from HVDC transmission lines based on corona discharge characteristics

High VoltageVolume 5, Issue 6 p. 679-687 Research ArticleOpen Access Prediction of radio interference from HVDC transmission lines based on corona discharge characteristics Caowei Huang, State Key Laboratory of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing, People's Republic of ChinaSearch for more papers by this authorHan Yin, State Key Laboratory of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing, People's Republic of ChinaSearch for more papers by this authorPengfei Xu, State Key Laboratory of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing, People's Republic of ChinaSearch for more papers by this authorBo Zhang, Corresponding Author shizbcn@tsinghua.edu.cn State Key Laboratory of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing, People's Republic of ChinaSearch for more papers by this authorJinliang He, State Key Laboratory of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing, People's Republic of ChinaSearch for more papers by this authorJianben Liu, State Key Laboratory of Power Grid Environmental Protection, China Electric Power Research Institute, Wuhan, People's Republic of ChinaSearch for more papers by this author Caowei Huang, State Key Laboratory of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing, People's Republic of ChinaSearch for more papers by this authorHan Yin, State Key Laboratory of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing, People's Republic of ChinaSearch for more papers by this authorPengfei Xu, State Key Laboratory of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing, People's Republic of ChinaSearch for more papers by this authorBo Zhang, Corresponding Author shizbcn@tsinghua.edu.cn State Key Laboratory of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing, People's Republic of ChinaSearch for more papers by this authorJinliang He, State Key Laboratory of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing, People's Republic of ChinaSearch for more papers by this authorJianben Liu, State Key Laboratory of Power Grid Environmental Protection, China Electric Power Research Institute, Wuhan, People's Republic of ChinaSearch for more papers by this author First published: 01 December 2020 https://doi.org/10.1049/hve.2019.0050Citations: 2AboutSectionsPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShare Give accessShare full text accessShare full-text accessPlease review our Terms and Conditions of Use and check box below to share full-text version of article.I have read and accept the Wiley Online Library Terms and Conditions of UseShareable LinkUse the link below to share a full-text version of this article with your friends and colleagues. Learn more.Copy URL Share a linkShare onEmailFacebookTwitterLinked InRedditWechat Abstract Many research works have proved that it is viable to predict the radio interference (RI) of high-voltage alternating current transmission lines by corona cage test. However, it is uncertain whether the same method can be adopted for RI prediction of high-voltage direct current (HVDC) transmission lines, since the corona performance of HVDC lines is obviously affected by the widely distributed space charge. In this study, the detailed corona current characteristics, such as the current pulse amplitude, rise time, half-wave time, and repetition frequency are systematically studied by the reduced-scale experiment and the relationship between the average corona current and the RI excitation current is confirmed. Furthermore, based on the experimental study, a method for predicting the RI level of HVDC transmission lines by the corona cage test is proposed. The RI level of a ±800 kV bipolar HVDC transmission line is predicted and compared with the measured result, which verifies the validity of the proposed method. 1 Introduction CORONA performances such as radio interference (RI) and audible noise (AN) are important restrictions for constructing new transmission lines [1]. Corona cages are usually used to study and predict the corona performance of high-voltage (HV) transmission lines. As there is no widely distributed space charge existing around the HV alternating current (HVAC) conductors, an excitation function, which is independent of conductor configurations, is defined to evaluate the RI generation quantities of HVAC transmission lines [2–4]. Therefore, the cage result is available for RI prediction of actual HVAC lines [5, 6]. However, since the space charge greatly affect the corona discharge of direct current (DC) conductors, it is hard to define a configuration-independent excitation function for HVDC transmission lines [7], which means that the cage result cannot be directly used. Thus some scholars have tried to make the RI prediction by monitoring overhead transmission lines [8–12]. However, using a corona cage is the most convenient and economical way to study corona performance. So lots of scholars, till now, have conducted researches on the corona cage test of HVDC conductors. Djeumen et al. [13] measured the space charge distribution in a corona cage under the influence of different atmospheric conditions. Liao et al. [14] also measured space charge in negative corona on a laboratory-scale corona cage. Nakano and Sunaga [7] proposed that corona cage tests were available for RI prediction based on the assumption that RI quantities were determined by the maximum bundle gradient in the presence of space charge. Otto and Reader [15] found that a small corona cage could not be used to accurately predict the RI performance of actual HVDC transmission lines. These studies in some extent proved that it was difficult to convert the corona cage result to actual line result under DC voltage. The characteristic of corona current pulse is the foundation in studying the relationship between corona cage corona and actual line corona. This characteristic has been studied by many scholars. Trichel [16] first discovered the periodically regular corona pulse waveform which is called Trichel pulse. Loeb [17] further developed the theory of Trichel pulse and categorised the negative corona mechanism and positive corona mechanism into three phases and four phases, respectively. Denholm [18] studied the positive corona current features, the peak value and pulse repetition frequency. However, these studies were limited to qualitative descriptions of the corona current characteristics. Recently, more and more specific corona currents were studied. Li et al.[19] studied the statistical characteristic in time-domain of DC corona and AN from the conductor in corona cage. Bian et al. [20, 21] studied the corona current characteristics under different air pressures and humidities. Yin et al. [22] made a research on the detailed characteristics of positive corona current for coaxial geometry electrodes, such as the current pulse amplitude, rise time, half-wave time and repetition frequency. Liu et al.[23] found the characteristics of intermittent positive corona current pulse and provided a reasonable explanation for this intermittent characteristics. Liu et al. [24] measured a large amount of corona current waveforms based on the largest outdoor HVDC corona cage, and found related statistical regularities between a corona current spectrum and an AN spectrum. Liang et al. [25] measured corona discharge characteristics for HVDC conductors with different height to the ground, but he mainly focused on the relationship between corona inception voltage, conductor structure and height. Mo et al. [26] studied the corona characteristics between the stranded and the shaped wires of HVDC transmission lines. Pfeiffer and Frank [27], Pfeiffer et al. [28] studied the relationship between the corona and the rain drop in raining condition. Although a large amount of valuable work has been reported on the corona current characteristics, few scholars investigated the corona current characteristic relationship between corona cage and actual line, which may be the key to investigating the RI of HVDC transmission line based on the corona cage test. In this paper, a reduce-scale experimental platform is set up and the detailed corona current features such as current pulse amplitude, rise time, half-wave time and repetition frequency are systematically studied in both corona cages and actual lines. Furthermore, the correlations among the average corona current, the RI excitation current and pulse repetition frequency are discovered and confirmed based on the experimental results. Finally, based on the corona cage test results, a method for predicting RI from HVDC lines is proposed and validated by the measurement on Xiangjiaba-Shanghai ±800 kV bipolar HVDC transmission line. 2 Traditional method As for HVAC conductors, the equivalence between the RI quantity measured in the cage and those on the line is based on the excitation function theory. According to the Shockley–Ramo theorem, the current induced in the electrode due to a particular charge in the inter-electrode space can be given as [29, 30] i e = q E q 0 ⋅ v (1)where ie is the induced current; q is the charge; Eq0 is the space-charge-free electric field at the charge located point due to unit potential applied to the electrode; v is the charge velocity. In the cylindrical corona cage, as shown in Fig. 1a, the induced current ic can be calculated as i c = q r q ln ( R / r ) v ⋅ E q 0 E q 0 = C c 2 π ε 0 q r q v ⋅ E q 0 E q 0 = C c 2 π ε 0 Γ c (2)where rq is the radial distance of the charge; r and R are the radii of the conductor and the cage, respectively; Cc is the capacitance per unit length of the conductor in the cage; ɛ0 is the permittivity of free space; Гc is the defined excitation function for the conductor in the cage, given as Γ c = q r q v ⋅ E q 0 E q 0 (3)As for the conductor above ground, as shown in Fig. 1b, the induced current il can be calculated as i l = 2 h q r q 1 r q 2 ln ( 2 h / r ) v E q 0 E q 0 = C l 2 π ε 0 2 h q r q 1 r q 2 v ⋅ E q 0 E q 0 = C l 2 π ε 0 Γ l (4)where h is the conductor height; rq1 and rq2 are the distances of the charge to the conductor and the image, respectively; Cl is the capacitance per unit length of the conductor above ground; Гl is the corresponding excitation function, given as Γ l = 2 h q r q 1 r q 2 v ⋅ E q 0 E q 0 (5)If the charge is in the vicinity of the conductor, as for HVAC lines, which means r q 2 ≃ 2 h, the excitation functions Гl and Гc can be simplified as the same. Therefore, the current il on the line can be predicted as i l = C l C c i c (6)However, as for HVDC lines, the existence of space charge region will significantly affect the electric field around the conductor. The space charge distribution around conductors in the corona cage is significantly different from that above the ground [31]. Therefore, the excitation functions Гl and Гc cannot be simplified as the same. Equation (6) is not valid for HVDC conductors, more studies should be carried out. Fig. 1Open in figure viewer Currents induced due to the space charge (a) Conductor in the cage, (b) Conductor above the ground 3 Reduced-scale test 3.1 Experimental platform The experimental platform consists of two systems as shown in Fig. 2a, one is the scaled conductor-plane test system and the other is the scaled cage test system. Both systems contain one corona onset section and two conjunct sections. The corona onset section is used to produce corona in both cage and line-to-plane experiments; the two conjunct sections are the power supply section providing energy for coronas and current measuring section, respectively. The measuring section can measure as well as record the corona pulse characteristics in both time domain and frequency domain. The photo of the real platform is illustrated in Fig. 2b. Fig. 2Open in figure viewer Experimental platform (a) Schematic representation of the platform, (b) Photograph of the platform The corona onset section for corona cage tests include a series of one-meter long corona cages, whose radii are 10, 20 and 30 cm, respectively. The corona onset section for corona line tests includes two insulating brackets designed to fix the reduced-scale lines so that the line-to-plane height can be adjusted from 20 cm to 1 m. The conductors used in both cage and line tests are the same and the radii are 0.4, 0.6 and 0.8 mm, respectively. The power supply section consists of a rectified DC supply which can provide high DC voltages and an inductive trap providing high DC voltages and blocking the interference from the power source. The current measuring section is the most complicated section in the experiment facilities, comprising a coupling circuit to decouple the radio frequency components of corona currents, an electromagnetic interference (EMI) receiver containing the CISPR standard and the peak (P), quasi-peak (QP) as well as root-mean-square (RMS) detectors, a Faraday's cage to protect a high-speed data acquisition card and a non-inductive resistor in its interior, a photoelectric conversion device with an optical fibre to transfer as well as deliver the current signals, and a computer to record the corona current characteristics. The high-speed data acquisition card has two channels and a 200 MHz sampling frequency which is high enough to analyse the radio frequency components of corona currents. The non-inductive resistor, whose resistance is 400 Ω, is connected in series with the line conductor to measure the corona currents. In order to verify the validity of the measurement system, measured waveforms of a typical positive corona current pulse from the non-inductive resistor and a current sensor are compared in Fig. 3. In the figure, A represents the peak value of the current amplitude, tr is the rise time, th is the half-wave time. The similar waveforms show that the measurement system has a good performance in measuring corona current. Fig. 3Open in figure viewer Typical positive corona current pulse measured by the non-inductive resistor and current sensor As demonstrated by Maruvada et al. [31], the RI is mostly induced by positive corona on bipolar HVDC transmission lines, therefore the focus point in this paper is the positive corona discharge. 3.2 Time-domain results The typically measured corona pulses when the conductor radius is 0.4 mm are shown in Figs. 4 and 5. In order to provide a common basis for any comparison, instead of using applied voltage to the conductor as the independent variable, the corresponding space-charge free surface electric field is considered as the independent variable. When the static electric field is 100 kV/cm, the corona current pulses change obviously with the size of the electrode. As demonstrated in Figs. 4a–c, with the increases of line-to-plane height h, the repetition frequency of the current pulses reduces evidently. The explanation for this observation might be that the space charge movement distance for larger h is longer than that in smaller h, and the larger one accumulates more charges resulting in the decrease of the total electric field. This decrease in total electric field weakens the discharge on the conductor surface. From another perspective, the average charge velocity is smaller due to the decrease of the total electric field, the region between line and plane can accumulate more charges, leading to further decrease of total electric field and the attenuation of discharges on the conductor surface. These two factors, together, result in less current pulses in a larger h model. For the same reasons, with the increases of corona cage radius R, the repetition frequency of the current pulses reduces apparently which has been demonstrated in Figs. 5a–c. Moreover, the current pulse repetition frequency in the corona cage is usually much higher than that on the actual line. Fig. 4Open in figure viewer Typical measured corona current pulses (a) Results on the reduced-scale line of 20 cm height, (b)Results on the reduced-scale line of 30 cm height, (c) Results on the reduced-scale line of 40 cm height Fig. 5Open in figure viewer Typical measured corona current pulses (a) Results in the cage of 10 cm radius, (b) Results in the cage of 20 cm radius, (c) Results in the cage of 30 cm Although the pulse repetition frequency decreases apparently with the increase in the line-to-plane height or the corona cage radius, the average pulse amplitude almost keeps unchanged at a constant value of a little >20 mA in both actual lines and corona cages, as illustrated in Figs. 4 and 5. Experiments are repeated more than ten times and the statistical deviations are illustrated in Fig. 6. Fig. 6Open in figure viewer Statistics of corona current pulses (a) Pulse amplitude, (b) Rise time, (c) Half-wave time, (d) Repetition frequency As shown in Fig. 6a, the maximum amplitude of corona pulses increases with the electric field but the average amplitude remains almost the same. As either the line-to-plane height or the cage radius increases, the maximum amplitude reduces to some extent, however, the average amplitude changes on a negligible scale. As the electrode spacing changes, the average rise time and half-wave time of corona currents change slightly, as shown in Figs. 6b and c. The average rising time is about 50 ns and the average half-wave time is around 160 ns. As demonstrated in Fig. 6d, when the electric field increases, the repetition frequency increases obviously. Additionally, the electrode size greatly affects the repetition frequency of corona pulses. When either the line-to-plane height or the cage radius increases, the repetition frequency reduces evidently under the same electric field condition. The repetition frequency of corona currents in the cages is much larger than that on the reduced-scale lines. However, this disparity grows smaller as the cage radius increases. The effect of the conductor radius on the corona characteristics is also studied. Similarly, the average amplitude, rising time and half-wave time of corona pulses change slightly with the variation of the conductor radius. However, the repetition frequency changes obviously with the conductor radius, as shown in Fig. 7. The reason for this observation can be explained in two aspects. On one hand, the conductor with a larger radius may have more discharge points on the conductor surface when the conductor surface roughness and conductor surface electric field are the same. On the other hand, the electric field is stronger in the surrounding area of larger radius conductors when the conductor surface electric field is identical, leading to a higher velocity of the ambient electrons and thus the ionisation process is easy to engender. Similar corona cage test results are also mentioned in [32]. Fig. 7Open in figure viewer Repetition frequency of current pulses in reduced-scale line corona of different conductor radii In conclusion, the electrode configuration greatly affects the repetition frequency of corona pulses. Therefore, the corona currents measured in the cage cannot be directly equivalent to those on the actual line. In the meanwhile, the statistics of a single corona pulse including the average amplitude, rising time and half-wave time are slightly affected by the electric field and the electrode configuration. 3.3 Frequency-domain result The radiofrequency components measured by the EMI receiver with the P, QP and RMS detectors are shown in Fig. 8. The correction factors for the circuit attenuation and resistive network are considered. As shown in the figure, the RMS responses are much smaller than the P and QP responses. For the purpose of comparison, the corresponding FFT values are also given. The RI results measured with a 200 Hz bandwidth, which are much lower than that with a 9 kHz bandwidth, are closer to the FFT value. Fig. 8Open in figure viewer Radiofrequency components measured with different EMI detectors According to the CISPR standard [32], QP values are recommended to evaluate the RI level. However, the excitation function, which is defined to describe the power spectral density of corona currents on the transmission line, should be expressed only in terms of RMS values. In the past two decades, RMS detectors have been widely adopted in the RI measurement [33]. In this paper, evaluating the RI performance of the transmission line by RMS values is recommended. The 0.5 MHz RMS components of corona pulse current which is defined as the RI excitation current, measured in corona cages as well as actual lines, are shown in Fig. 9. When the electric field is large, the RMS values in dBμA increase almost linearly with the electric field. When the electric field is close to the onset value, the RMS values decrease rapidly. As for the cage test, the RI level is higher when the cage radius is smaller. As for the reduced-scale line test, the RI level is higher when the conductor height is lower. Additionally, the RI level in the cage is much larger than that on the reduced-scale line. Fig. 9Open in figure viewer 0.5 MHz RMS components of corona currents measured in the corona cage and on the reduced-scale line 4 Correlations of corona currents After a long-time observation, a large-sized corona current pulse sequence can be acquired. The average corona current can be calculated as I 0 = ∫ 0 T all i ( t ) d t T all (7)where I 0 is the average corona current, T all is the duration time of the pulse sequence, i ( t ) is the waveform of corona current. According to the previous conclusion, the waveform of a single corona pulse remains unchanged. It is assumed there are n pulses in the sequence, then (7) can be modified as I 0 = ∫ 0 T all ⁡ i ( t ) d t T all = n ∫ 0 t 0 ⁡ i s ( t ) d t n T av = ∫ 0 t 0 ⁡ i s ( t ) d t T av = A f (8)where n is the number of pulses in the sequence, t 0 is a single corona pulse duration time; i s ( t ) is the single corona pulse waveform function; T av is the average time interval between two consecutive corona pulses; f is the pulse repetition frequency and A is the scale factor. The factor A represents the average charge provided by a single corona pulse. From (8), it can be seen that a single corona pulse waveform can be used to calculate the I 0 of the entire pulse sequence. The relationship between the average corona current I 0 and the pulse repetition frequency f is acquired. The RI excitation current, an important characterisation of the corona discharge intensity, is generally measured by EMI receivers. To find out the relationship between the RI excitation current and the current repetition frequency, a simulation for the receivers' measurement process is established. Substantially, this measurement process is a filtering process and furthermore, the receiver can be simplified to a Hanning band pass filter [34] whose characteristic function is given as H f k = 0.5 − 0.5 cos 2 π ( f k − f L + ( f B / 2 ) ) 2 f B , f L ≤ f k ≤ f U 0 , f k f L or f k f U (9)where fIF = 0.5 MHz is the filter centre frequency; fB = 9 kHz is the filter bandwidth; fL = 0.4955 MHz is the lower cut-off frequency and fU = 0.5045 MHz is the upper cut-off frequency. The filtered current if is defined as i f = Re F − 1 ( H ( f k ) F ( i t ) ) (10)where F is the fast Fourier transform, F − 1 is the inverse fast Fourier transform. Similarly as (8), the RMS value of i f can be expressed as I RMS = ∫ 0 T all ⁡ i f 2 ( t ) d t T all = n ∫ 0 T f ⁡ i f 2 ( t ) d t n T av = ∫ 0 T f i f 2 ( t ) d t f = B f (11)where I RMS is the RMS value of RI excitation current; T all is the duration time of the pulse sequence, n is the number of pulses in the sequence, T f is the duration time of filtered current i f , f is the repetition frequency of i f , which is also the repetition frequency of corona current i ( t ); B is the scale factor and B2 represents the RI component of average power provided by a single corona pulse. To confirm (8) and (11), the relationships between measured average corona current I 0 , measured RMS value I RMS and the current pulse repetition frequency f are shown in Fig. 10. It is obvious that the average corona current I0 varies linearly with the increase of frequency and meanwhile the RMS value of filtered corona current I RMS presents a linear dependence with f .This relationship is universal for all electrode configurations. Fig. 10Open in figure viewer Variation of measured corona current components with the repetition frequency Once the (8) and (11) are confirmed, the correlation between I0 and IRMS can be rewritten as I RMS = B A I 0 (12)For corona cages and reduced-scale lines, respectively, (12) can be expressed as I RMS , l = B l A l I 0 , l (13) I RMS , c = B c A c I 0 , c (14)where Bl, Al and Bc, Ac are the scale factors for reduced-scale lines and corona cages respectively; I0,l, IRMS,l and I0,c, IRMS,c are the average corona current and the RMS value of RI excitation current for reduced-scale lines and corona cages, respectively. As already mentioned in (8) and (11), the factor A represents the average charge provided by a single corona pulse and the factor B2 represents RI component of average power provided by a single corona pulse. Therefore, the parameters A and B depend only on the waveform of a single corona current pulse. Moreover, as shown in Fig. 6, whether the electric field or the electrode configuration changes, the average amplitude, rise time and half-wave time of corona pulses practically remain unchanged. Therefore, if the surface condition of the conductor and ambient conditions are same, the waveform of a single current pulse is almost the same whether in line coronas or in cage coronas. This result can be expressed as A l ≈ A c , B l ≈ B c (15)According to (13), (14) and (15), the relationship between IRMS,l and IRMS,c can be represented as I RMS , l = I 0 , l I 0 , c I RMS , c (16)Based on the measured experimental results of I0,l, IRMS,l and I0,c, IRMS,c, the further comparisons between the ratio of IRMS,l/IRMS,c and the square root ratio of I0,l /I0,c is demonstrated in Fig. 11. The measured values with subscript c are acquired from corona cage results for different values of radii from 20, 30 and 40 cm. The measured values with subscript l are acquired based on the test results from the reduced-scale line of height 20 cm. It is obvious that among all electrode configurations, the square root ratio of average corona current matches with the corresponding ratio of RI excitation current well. Besides, the result of Cl/Cc which presents the corona current relationship between the cage coronas and line coronas as shown in (6), is also provided as a contrast. As shown in Fig. 11, it has been clearly demonstrated that because the neglecting of the space charge, the traditional excitation function cannot be used directly in DC corona cage test, otherwise a misguiding of much larger line excitation current would mislead a much larger predicted corona level in transmission lines. This observation is consistent with the research of Maruvada et al. [33]. Fig. 11Open in figure viewer Relationship between average corona current and RI excitation current in corona cages as well as on reduced-scale lines Based on the (16), steps to confirm the relationship of excitation currents from corona cage and the actual line can be concluded as follows: Step 1: Acquire RI excitation currents in different space-charge-free surface electric field values based on the corona cage test, the conductor should be the same as used in the actual line. Step 2: Among the corona cage test results, the ones with the same space-charge-free surface electric field values as the actual lines are selected as the IRMS,c. Step 3: I0,c can be obtained by corona cage test or numerical method. Step 4: I0,l can be obtained by numerical method. Step 5: Based on (16), the actual line RMS value of RI excitation current IRMS,l can be calculated as I dB , l = 20 lo g 10 I RMS , l = 20 lo g 10 I 0 , l I 0 , c I RMS , c = 20 lo g 10 I RMS , c + 10 lo g 10 I 0 , l I 0 , c (17)where IdB,l is the decibel of IRMS,l. In comparison, for the traditional method, the actual line RI excitation function can be calculated as I dB , l = 20 lo g 10 I RMS , c + 20 lo g 10 C l C c (18) 5 Calculation of average corona current Equation (17) provides a method of converting the RI excitation current from corona cage to transmission lines. In order to obtain the average corona current I0,c and I0,l in step 3 and step 4, an ef