A fast boundary protection for an AC transmission line connected to an LCC-HVDC inverter station

For AC transmission lines connected to an LCC-HVDC inverter station, commutation failure can lead to the wrong operation of traditional protection. To solve the problem, this paper proposes a fast protection scheme using transient information from one end of the AC line. The boundary frequency characteristics of the AC line connected to LCC-HVDC inverter are analyzed first. This reveals the existence of significant attenuation on both high frequency signals and some specific frequency signals. Based on the boundary characteristics, a novel boundary protection principle utilizing a backward traveling wave is then proposed. A PSCAD/EMTDC simulation model of a ± 800 kV LCC-HVDC and 500 kV AC transmission system is established, and different fault cases are simulated. The simulation results prove that the novel protection principle is immune to commutation failure, fault resistance and fault type.


Introduction
An HVDC transmission system has many advantages over an AC transmission system such as long transmission distance, high efficiency and large transmission capacity [1][2][3][4][5][6]. However, as an LCC-HVDC inverter station is severely affected by an AC transmission line (ACTL) fault, it is important that an ACTL fault is cleared quickly and reliably. Otherwise, it could lead to commutation failure or even an HVDC system blocking [7][8][9][10][11][12]. Commutation failure and complex electrical dynamics during the recovery process result in the backside system impedance exhibiting nonlinear time-varying characteristics, whereas the traditional AC system fault analysis theory is based on linear systems. Thus, the adaptability of conventional power frequency protection (PFP) in LCC-HVDC connected AC/DC hybrid systems has problems, which have been extensively and deeply studied [13][14][15][16][17][18].
When a commutation failure occurs due to an AC system fault, the amplitude of the equivalent power frequency current decreases and the fluctuation range of the phase angle is large. This may cause current differential protection to fail [14]. In addition, the directional pilot protection principle cannot adapt to an ACTL connected to an inverter station [15]. Since the positive and negative sequence impedances are not equal for an ACTL connected to an inverter station, it may lead to mis-operation of the directional element using sequence components [17]. The DC part exhibits a large shortcircuit impedance when an AC line fault occurs, and this can cause distance protection to reduce or lose its protection range [16], while the main electromagnetic transient process caused by commutation failure influences AC line distance protection [18].
However, theoretical analysis and on-site recording indicate that the time interval from fault to first commutation failure is generally more than 3 ms [19,20]. Therefore, if protection can activate a tripping order for an ACTL fault within 3 ms of the fault, the consequent commutation failure will not affect the correct action of the protection. However, none of the power frequencybased protection can meet such a requirement, and it is only possible using transient-based protection.
Transient-based protection directly utilizing fault transient information to quickly remove a fault has been studied over the past four decades [21][22][23][24][25][26][27]. It can be divided into unit and non-unit transient-based protection from the perspective of space utilization information. Compared with unit transient protection, non-unit transient protection only uses locally measured transient fault information and can accomplish ultra-high-speed operation without communication delay. Traveling wave distance protection [23] and boundary protection are both non-unit protection. Boundary protection can distinguish internal and external faults from the fault characteristic differences caused by the change of surge impedance at the line boundary.
The line boundary in a primary system and boundary characteristics determine the performance of boundary protection. The earliest line boundary is the conventional line trap and specially designed stack tuner aimed at capturing a high-frequency voltage signal in EHV transmission systems [28,29]. However, adjusting the line trap and installing a stack tuner are difficult in practice. The busbar stray capacitance is regarded as the line boundary in [30]. This has severe attenuation on high-frequency transient current. Some researchers consider both line trap and busbar stray capacitance as line boundary [31]. Since fiber optic communication is now widely used, line traps are rarely used and the scope of their application is limited. To improve the performance of boundary protection, the line boundary and corresponding characteristics in UHVAC transmission systems have been analyzed comprehensively based on real projects. Busbar stray capacitance, the shunt reactor, series capacitance and combined line boundary are used as the line boundary to distinguish internal and external faults considering fault characteristic differences, while a more precise setting method is proposed [1].
Compared with that for the AC transmission line, the development of boundary protection for the DC transmission line has made great progress. In the primary HVDC transmission system, DC filters and smoothing reactors at each line terminal constitute the HVDC line boundary. This line boundary imposes a significant attenuating or smoothing effect on transient signals, and makes significant characteristic differences between internal and external faults [32][33][34][35][36]. However, study on composition and characteristics of an ACTL connected to an inverter station is still limited.
The extent of full usage of the fault characteristic differences caused by line boundary also determines the performance of boundary protection. Digital filters [29,30], wavelet transform, neural network and mathematical morphology are applied in boundary protection algorithms [37][38][39][40], but they still rely on the line boundary composed of busbar stray capacitance [37][38][39][40].
To develop a novel protection not affected by commutation failure, a boundary protection (BP) principle utilizing a backward traveling wave is proposed because of the boundary specificity of an ACTL connected to an inverter station. The proposed method can respond to faults quickly (in under 3 ms) and is immune to commutation failure. The boundary composition and characteristics of the AC line connected to an LCC-HVDC inverter are analyzed first and based on the boundary characteristics, the novel boundary protection principle using a backward traveling wave is then proposed. A synchronous squeeze wavelet transform-based algorithm is established to extract the energy of a narrow band signal. Finally, the effectiveness of the proposed protection is validated by simulation tests of the ±800 kV HVDC system in China.
This rest of the paper is organized as follows. Boundary characteristic analysis is investigated in Section 2 and the BP algorithm using a backward traveling wave is described in Section 3. The simulation model of the ± 800 kV LCC-HVDC and 500 kV AC transmission system is established in Section 4, and extensive simulations are carried out to verify the effectiveness of the proposed protection. Finally, conclusions are drawn in Section 5.

Boundary composition
The line boundary is the discontinuous point of the surge impedance of a transmission line, and is usually located at both ends of a protected line. Based on the parameters of the Tian-Zhong ±800 kV / 8000 MW HVDC system in China, a typical DC-AC hybrid system is shown in Fig. 1, where Cs represents the busbar stray capacitance. As shown, in order to compensate for the reactive power absorption of the LCC converter and restrain the generated AC harmonics, AC filters and a reactive power compensation device (RPCD) are installed on the busbar connected to the AC transmission lines. DC filters and a smoothing reactor (SR) are installed at both ends of the DC transmission line to retrain DC harmonics. These reactive power elements form both the AC boundary and the DC boundary of the inverter station as shown in Fig. 1.

Frequency characteristics of AC boundary of inverter station
As can be seen in Fig. 1, the AC boundary consists of the AC filters, RPCD and EHV busbar capacitance. The parameters and structure of the AC filters and RPCD can be obtained from [41]. The configuration of the inverter station AC filters and shunt capacitors is shown in Table 1.
The reactive power compensation device is a high voltage shunt capacitor (SC), which generally consists of a high voltage capacitor and a small damping reactance connected in series. The purpose of the damping reactance is to limit inrush current when the capacitor is switched in. As indicated in some studies, the inrush current generated when switching in a single group of capacitors is not severe, and thus a single shunt capacitor can be switched directly without the use of the series reactance. However, when an existing capacitor group is already connected and another group is switched in, the inrush current can be very large and the small damping reactance is required. Considering the economical aspect, there is only one set of shunt capacitors not having the small damping reactance in the whole station [42]. The small damping reactance ranges from 1 to 5 mH, with a typical value of around 2 mH. Therefore, the high voltage shunt capacitors are divided into two cases: a shunt capacitor with a small damping reactance (SC) and shunt capacitor without the damping small reactance (SC0).
The detailed filter arrangements for the 4 large filter groups at the inverter station are as follows: Groups 1 and 2: each contains two HP12/24 filters, one HP3 filter and 2 SC. Group 3: two HP12/24 filters and 3 SC. Group 4: two HP12/24 filters and 2 SC, of which one SC has no small damping reactance (SC0).
The AC filters and RPDC are all connected in parallel, which is the same as the ground stray capacitance for the AC substation [1]. The case of the line boundary is shown in Fig. 2a and the lumped parameter equivalent circuit in the Laplace domain is depicted in Fig. 2b. Because the surge impedance is discontinuous at the line boundary of a transmission line whose surge impedance is marked as Z 1 in Fig. 2, the incident traveling wave u 1b generated by a fault is reflected to form the reflected wave u 1f . It is then refracted to form the refracted wave u 2f , which enters the other transmission line whose surge impedance is Z 2 . Therefore, these devices connected in parallel have the generic boundary transfer characteristic expression: where sign // represents parallel relation; H(ω) denotes the transfer characteristic of the line boundary, i.e. the refractive coefficient, which is related to the parameters  of the line boundary and the line surge impedances Z 1 and Z 2 .
Assuming Z 1 = Z 2 = 250 Ω, then based on the Peterson principle and the system parameters, the AC boundary characteristics of the inverter station when there is an individual AC filter, the SC and filter group are considered separately can be respectively obtained, as shown in Fig. 3 and Fig. 4.
The amplitude-frequency characteristics of the AC boundary are summarized as follows: a) The reactive power compensation shunt capacitor without small damping reactance (SC0) is equivalent to directly connecting a large capacitor to the ground, which is tens to hundreds of times the stray capacitance of the busbar to ground in an AC system, and the line boundary formed has significant attenuation on high frequency signals above 10 kHz. b) The line boundary formed by the shunt capacitor with small reactance (SC) only has attenuation on a narrow band signal centered on the resonance frequency generated by the small reactance and capacitor, while it has no obvious attenuation on other frequency components. This resonance frequency and the narrow band signal are called the specific frequency and specific frequency component, respectively. c) The line boundary formed by the AC filter (HP12/ 24) has attenuation only on a narrow band signal centered on the tuned frequency, while it has no obvious attenuation on other frequency components. Similarly, this tuned frequency is the specific frequency while the narrow band signal is the specific frequency component. The line boundary formed by the low-order harmonic filter (HP3) has limited attenuation on the narrow band signal centered at the tuned frequency of 150 Hz. d) The boundary characteristics of each of the 3 largegroup filters comprehensively reflect and strengthen the boundary characteristics when considered alone and have significant attenuation on a specific frequency signal and high frequency signals.
Based on the above analysis, it can be concluded that the AC boundary of an LCC-HVDC inverter station has significant attenuation both on high frequency signals and specific frequency signals. These will cause different fault characteristics between internal and external faults.

Frequency characteristics of DC boundary of the inverter station
DC filters and the smoothing reactor form the DC boundary of the inverter station. The parameters and structure of the DC filters and SR can be obtained from [41]. The DC filters are all connected in parallel, which is the same as the AC filters and RPDC, while the smoothing reactor is connected in series. The configuration of the inverter station DC filters is shown in Table 2. The case of the line boundary is shown in Fig. 5a and the lumped parameter equivalent circuit in the Laplace domain is depicted in Fig. 5 (b). Therefore, these devices connected in series have the generic boundary transfer characteristic expression: The DC boundary characteristics of the inverter station are studied in the following ways: The amplitude-frequency characteristics of the DC boundary are summarized as follows: a) The boundary formed by the DC filters (HP12/24 and HP3/39) has no attenuation effect on high frequency signals, whereas it has significant attenuation on a narrow band signal centered on the tuned frequency (called as specific frequency). The line boundary formed by the DC filter (HP3/ 39) has limited attenuation on the narrow band signal centered at the tuned frequency of 100 Hz. b) The line boundary formed by SR has a significant attenuation effect on high frequency signals above 10 kHz. A larger SR inductance value leads to higher attenuation. c) The combination of the DC filters and SR comprehensively reflects and strengthens the boundary characteristics. This makes the attenuation of specific frequency components and high frequency components more severe.
Based on the above analysis, it can be concluded that the DC boundary of an LCC-HVDC inverter station has significant attenuation on high frequency signals and specific signals. The characteristics will make significant differences between high frequency or specific frequency components of the internal fault and the external fault. These can be used to form the boundary protection.
The AC boundary and DC boundary of the inverter station together form the boundary of the ACTL connected to the LCC-HVDC inverter station. By considering the frequency characteristics indicating that they have significant attenuation on both the specific frequency components and high frequency components, a new boundary protection principle can be formed.

Protection principle
Assuming the current reference direction at the relaying point is from the busbar to the protected line and observing the relaying point, the forward traveling wave Δu f and the backward traveling wave Δu b are: where Δu and Δi are the fault components of voltage and current measured at the relaying point, respectively. Z denotes the surge impedance of the protected line.  Assuming boundary protection is located at the line terminal M of the protected line MN as shown in Fig. 1, in the case of an internal fault at f 1 , the Bewley lattice diagram of the fault generated traveling waves is shown in Fig. 7. B 1 and F 1 are the initial backward and forward traveling waves, respectively. The subsequent backward and forward traveling waves mainly include the following: In the case of a forward external fault at f 2 , the Bewley lattice diagram of the fault generated traveling waves is depicted in Fig. 8. B 1 and F 1 are the initial backward and forward traveling waves, respectively. The subsequent backward and forward traveling waves mainly include the following:   Boundary protection detects the fault at time t 0 . Before B 3 arrives at the relaying point, i.e., within the time interval [t 0 , t 0 + 2l MN /c), where l MN denotes the length of line MN, all the backward traveling waves detected by boundary protection (B 1 , B 2 , when the line NP is short, B 4 may also be included) do pass through the line boundary N. Therefore, amplitudes of the high frequency and specific frequency components in the backward traveling waves will be attenuated by the line boundary significantly.
In the case of a backward external fault at f 3 , the Bewley lattice diagram of the fault generated traveling waves is shown in Fig. 9. At the instant t 0 , the initial forward traveling wave F 1 reaches the relaying point M, and boundary protection detects a fault occurrence at time t 0 . Before B 1 (stemming from the opposite terminal boundary N) arrives at the relaying point, i.e., within the time interval [t 0 , t 0 + 2 l MN /c), the backward traveling wave does not appear. Thus, the backward traveling wave detected by boundary protection is zero and theoretically, E f /E b becomes infinite. However, for a forward fault, the backward wave detected by boundary protection is the initial traveling wave and the subsequent backward and forward traveling waves meet the law of reflection, so within the time interval [t 0 , t 0 + 2 l MN /c), the backward and forward traveling waves can be expressed as: where k f is the reflection coefficient of the bus. Therefore, it can be concluded that within the time interval [t 0 , t 0 + 2 l MN /c), the time-domain energy ratio of the forward and the backward traveling waves E f /E b is less than 1 because |k f | ≤ 1.
From the above analysis, there exist significant differences in the backward traveling waves between internal and external faults. Thus, a novel principle of boundary protection is proposed to effectively utilize such differences.
After detecting an initial traveling wave generated by a fault, in the period of Δt, the time-domain energy ratio of the forward and backward traveling waves E f /E b , the energies of the high frequency and specific frequency components of the backward traveling waves E h and E r , are calculated respectively. If E f /E b > ε 0 , the fault is identified as a backward external fault, otherwise it is a forward fault and the product of the high frequency energy and specific frequency energy (E h E r ) is calculated. If E h E r > ε 1 , the fault is determined to be internal, otherwise it is a forward external fault.
The protection principles can be further described as follows: (1) The threshold ε 0 is used to distinguish the fault direction. In the case of backward faults, within the time interval Δt, the backward traveling wave is theoretically zero. Thus, considering sensitivity and reliability, the threshold value ε 0 is 2 in this paper.
To ensure that the fault direction characteristic is strictly established, Δt < 2l sh /c, where l sh is the length of the shortest line among all the lines connected to the same busbar.
(2) Transient current is a high frequency signal, which is easily influenced by high frequency noise. A mathematical morphology filter (MMF) with a flat structure element is applied to remove high frequency noise and eliminates the possibility of incorrect operation caused by the sharp part of the waveform.  It has no effect on the accuracy of the boundary protection and can ensure reliability. (4) The boundary elements of R 2 depend on the characteristics of the line boundary N. For ultrahigh voltage AC lines, the line boundary N consists of only the busbar ground stray capacitance which has significant attenuation on high frequency components [1]. Thus, the time-domain energy of the high frequency components can be used to identify the fault. (5) Clarke phase-modal transformation is applied in this paper to obtain the α, β and 0 modal voltages and currents. Boundary protection uses the aerialmode components due to the severe attenuation of the zero-mode component. Comparing the absolute magnitudes of the α mode and β mode, the larger is selected as the significant fault mode to cover all ten types of fault and form the backward traveling wave in the above boundary protection principle.

Synchronous squeeze wavelet transform-based algorithm
A synchronous squeeze wavelet transform (SST) transforms the time-scale (a,b) plane onto the time-frequency (ω l, b) planes to make W x (a,b) limited to the frequency range with the center frequency of ω l , by compressing the wavelet transform coefficients W x (a,b). In the frequency range where the center frequency is ω l , the timefrequency energy is concentrated near the center frequency and there is no cross term in each frequency curve, thereby effectively reducing band aliasing. SST can obtain the band energy of the signal over a certain period ([t 1 , t 2 ]) in a certain frequency interval ([f 1 , f 2 ]), as shown in (6). Compared with DWT, SST divides the frequency space more precisely, so it can be used to extract the energy of a narrow band signal, such as the specific frequency components and high frequency components.
For the shunt capacitor with small damping reactance at the inverter station in the Tianzhong AC-DC hybrid system, its resonant frequency f r is: A flowchart of the boundary protection algorithm is shown in Fig. 10.

Simulation and results analysis
Based on the Tian-Zhong (in China) ± 800 kV / 8000 MW HVDC system, a typical AC-DC hybrid simulation model is established as shown in Fig. 1. Twenty-fourpulse converter valves are installed in the converter station. The lengths of the HVDC transmission line and protected AC transmission line, both modelled using frequency dependent models, are 2190 km and 62.5 km, respectively. The other AC transmission lines are 21.27 km and 16.87 km long. The transmission line structures are shown in Fig. 11 (a) and (b) for AC and DC, respectively. As for control strategies, the rectifier station adopts the conventional constant current control strategy, and the inverter station adopts conventional constant current control and constant extinction angle control strategies.
The sampling rate is 400 kHz and the total data window Δt = 1.28 ms. The specific frequencies take the two center frequencies closest to f r , i.e., f 0 = 2040.  is arranged at the N end. For MMF, a flat structure element is employed and the size is 3. The data window of the direction element Δt dir = 0.15 ms (N dir = 60) and the threshold ε 0 = 2. The setting method for boundary protection of an adaptive fault type is adopted, and the threshold ε 1 is respectively determined in terms of fault types: ε ϕg = 4.56e13, ε ϕϕg = 5.65e15, ε ϕϕ = 7.32e15, and ε ABC = 1.25e16.
For an internal metallic phase A grounding fault at the middle of line MN, which is 32.5 km from busbar N with the inception angle of 90°, as shown in Fig. 12, the fault occurs at 0.32 ms and the total fault data is 1.28 ms (0.32-1.6 ms). The fault direction is identified as a forward fault because of the energy ratio of the forward and backward traveling wave being E f /E b = 0.401034. At 1.6 ms, it calculates E h E r = 6.05e19, which is larger than the threshold ε ϕg . Thus, at only 1.28 ms after the fault occurrence, boundary protection operates.
For a forward external metallic phase A grounding fault at the middle of another line MQ at 0.32 ms, which is 15 km from busbar M with phase A inception angle of    90°, the results are shown in Fig. 13. The fault direction is identified as a forward fault because of the energy ratio of the backward and forward traveling waves E f /E b = 0.534551. At 1.6 ms, it calculates E h E r = 4.12e12, which is less than the threshold ε ϕg . Thus, boundary protection does not operate. For a backward external phase A grounding fault at busbar N at 0.32 ms, with the fault resistance of 5 Ω and the phase A inception angle of 90°, the results are shown in Fig. 14. As can be seen, the data window time for calculating E f /E b is 0.32-0.47 ms, and only 0.32 ms after the fault occurrence. It calculates E f /E b = 1.6e6, which is larger than the threshold ε 0 up to six orders of magnitude. Thus, the fault is reliably determined to be backward, and boundary protection does not operate.
The boundary protection is further tested under the following different fault conditions (f indicates the fault point) and Tables 3, 4 Tables 3, 4 and 5. The smaller the initial fault angle is, the smaller the E h E r value is, and the sensitivity of boundary element discrimination will decrease. When the fault with high resistance and small inception angle occurs at the end of the line, the boundary element will refuse to operate.
The extensive simulations in Tables 3, 4 and 5 show that most forward internal and external faults can be distinguished correctly within 1.28 ms after a fault occurrence. Only for single-phase grounding faults at the end of the protected line with high fault resistance or small fault inception angle, the sensitivity of boundary protection is reduced and boundary protection may fail to operate (such as the phase A grounding fault with l Nf = 62.5 km and the inception angle of 5°in Table 5). It also indicates that it is difficult to reliably protect the whole line by only utilizing the non-unit transient fault information. The simulation results prove that the novel protection principle is immune to commutation failure, fault resistance and fault type.

Conclusions
This paper comprehensively analyzes the inherent frequency characteristics of an inverter station boundary. Both the AC boundary and DC boundary have significant attenuation on high frequency components and specific frequency components. The AC filters, shunt capacitors with small damping reactance and DC filters have attenuation on specific frequency components dependent on the structure of the filters, whereas shunt capacitors without small damping reactance and SR have significant attenuation on high frequency components.
Based on the boundary characteristics, a novel boundary protection principle using a backward traveling wave is proposed for an ACTL connected to an inverter station. As this protection scheme does not need a  communication channel to receive information from the remote end, it can identify a fault in under 3 ms. Therefore, the consequent commutation failure will not affect its correct detection of the fault. Simulation results show that the boundary protection operates stably under various fault conditions. However, it is difficult to reliably protect the whole line under some heavy fault conditions and further studies and explorations are needed to improve this scheme.