Multiterminal Hybrid DC Line Protection Based on Intrinsic Mode Energy Entropy

In view of the UHV multiterminal hybrid DC transmission system, the DC line protection with universal applicability, fast response speed, stability, and reliability is particularly important for its safe and economic operation. In this article, the boundary frequency characteristics of the UHV multiterminal hybrid DC transmission system are analyzed by considering the DC boundary of the rectifier side and inverter side, line frequency variation, and other factors. According to the fact that the boundary of the fault line has a strong attenuation effect on the high-frequency component of the transient current, a protection method based on the intrinsic mode energy entropy is proposed to distinguish the faults inside and outside. The criterion to initiate the protection is constructed by using the amplitude of the transient power; the fault direction criterion is constructed by using the positive and negative characteristics of the transient energy of the fault power detected by power direction elements on both sides of the T-zone; the fault pole selection criterion is constructed by using the ratio of low-frequency transient power changes of positive and negative poles. Finally, the Kun–Liu–Long line UHV three-terminal hybrid DC transmission system model is built on the PSCAD/EMTDC simulation platform, and a large number of simulation examples verify that the protection can operate reliably under different fault poles, fault positions, and transition resistances.


INTRODUCTION
In recent years, large-scale new energy grid connection, regional power grid interconnection, and central load power supply have become hot issues at home and abroad, and it is urgent to build a DC network with higher power supply reliability on the basis of existing practical experience and theoretical research on DC transmission projects (Tang et al., 2013;He et al., 2017). Most of the traditional direct current transmission systems are two-terminal systems, which can only realize point-to-point power transmission whose remarkable feature is to realize AC collection-DC transmission-AC dispersion of electric energy (Wang et al., 2014). Compared with the traditional two-terminal DC transmission system, the UHV hybrid DC transmission system contains multiple sending terminals or receiving terminals, which can realize multipower supply and multi-drop power reception, and the operation mode is more flexible (Xu et al., 2013). The sending terminal is usually connected with a line-commutated converter (LCC), while the receiver is usually connected in parallel with a modular multilevel converter (MMC) with fault self-clearing capability (He et al., 2020). This structure makes comprehensive use of the advantages of low loss and good economy of conventional DC power transmission as well as flexible and fast control of flexible DC power transmission without commutation failure, such as Kun-Liu-Long line UHV multiterminal hybrid DC transmission demonstration project planned and constructed in China (Li et al., 2019).
The long distance of DC transmission lines and the complex environment of transmission corridors lead to a high fault rate Liu et al., 2020). When the DC fault occurs, it develops very quickly. The superposition of fault output of multiple converter stations in the UHV multiterminal hybrid DC transmission system makes the fault hazards such as line overcurrent more prominent (Muniappan, 2021). At present, the main protection widely used for the DC transmission lines is traveling wave protection (Shen and Raksincharoensak, 2021a). However, in actual operation, traveling wave protection is easily affected by factors such as inaccurate extraction of wave head information, selection of setting value, and high resistance ground fault, while its protection reliability is poor (Zheng et al., 2018). Xia et al. (2018) use the longitudinal impedance at both ends of the transmission line to propose a multiterminal line differential current protection, but there is a problem of too many measured values, and the communication time required increases when the line is long. Zhou et al. (2017) put forward a multiterminal flexible DC line protection scheme based on the voltage amplitude of a current-limiting reactor by using the characteristics of DC reactors as line boundary elements. However, there is a problem of poor transient resistance, and the multiterminal DC lines are interconnected by DC buses, and the T-zone structure formed by interconnection can be used directly with relatively few line boundary elements. Li et al. (2019) study the unique structure of T-zone in the parallel multiterminal hybrid DC transmission system and propose a fault direction discrimination principle based on the comparison of wavelet energy of fault current on both sides of T-zone. When applied to the three-terminal hybrid DC system, there is no need for interstation communication, but the selection of its setting value depends on the result of wavelet decomposition. In reference , the transient current fault components at the three ends of the T-junction bus are measured, and the correlation coefficient is calculated to serve as the basis for fault zone identification of the hybrid multiterminal DC transmission system. No interstation communication is required, but the noise tolerance is poor. Therefore, it is necessary to further study the DC line protection which is suitable for the UHV multiterminal hybrid DC transmission system and has fast response, stability, and reliability. Yang and Liu (2018) take the current mutation of each phase in the AC distribution network as the characteristic quantity, calculate its inherent modal energy entropy, and use the maximum modal energy maximum value for fault line selection. It is less affected by grounding the resistance and system structure and is more suitable for systems with complex structures. This principle provides a certain idea for the study of multiterminal hybrid HVDC transmission line protection.
In this article, the Kun-Liu-Long line UHV three-terminal hybrid DC transmission system is selected as the research object, and the boundary frequency characteristics of the rectifier side and the inverter side are analyzed theoretically. On this basis, a protection method is proposed to distinguish faults, inside and outside, based on the intrinsic mode energy entropy. The amplitude of transient power is used as the criterion to initiate the protection; the positive and negative characteristics of transient energy of fault power detected by power direction elements on both sides of the T-zone are used to identify the fault direction; the ratio of the variable quantity of low-frequency transient power of positive and negative poles is used to select the fault poles. The model of the Kun-Liu-Long line UHV threeterminal hybrid DC transmission system is built on the PSCAD/ EMTDC simulation platform, and the effectiveness of the protection scheme is verified by a large number of simulation examples. Figure 1 shows the simulation model of the Kun-Liu-Long line UHV three-terminal hybrid DC transmission system. The system adopts the true bipolar connection mode, with independent converter equipment, transmission lines, and control systems for both positive and negative poles. The LCC is used in the Kunbei Converter Station on the rectifier side, while the bipolar MMC with fault self-clearing capability is used in the Liubei Converter Station and Longmen Converter Station on the inverter side. The sending terminal Kunbei Converter Station and the receiving terminal Longmen Converter Station are connected to the bus bar in the receiving terminal Liubei Converter Station by line L 1 and line L 2 , respectively (Yu et al., 2020). The T-zone is formed by the DC line L 1 and L 2 connected in parallel with the bus bar, and the protection devices M1 and M2 are installed on the left and right sides of the T-zone, respectively.

The Analysis of Boundary Frequency Characteristics on the Kunbei Side Line
Because the rectifier side of the Kun-Liu-Long line UHV threeterminal hybrid DC transmission system adopts LCC, during the normal operation, due to the nonlinear characteristics of the converter, a large number of harmonic waves will be generated in the DC transmission line (Shen and Raksincharoensak, 2021b). Therefore, it is necessary to configure two smoothing reactors and two groups of triple-tuned DC filters at the head of the DC line L 1 to suppress the harmonic current generated by the converter and injected into the DC line (Shen et al., 2017, Shen et al., 2020a. The smoothing reactors and DC filters together form the physical boundary on the rectifier side, as shown in Figure 2. The specific parameters of each component are shown in Table 1. In Figure 2, the symbol (I 1 ) represents the transient current outside the protection zone of the DC line L 1 , while the symbol (I 2 ) represents the transient current flowing into the DC line L 1 after attenuation by the line boundary on the Kunbei side line. According to the parameters in Table 1, the expressions of equivalent impedance of the smoothing reactor and DC filter can be obtained, respectively, as follows: (1) In the formula, the symbol (Z 2 (jω)) represents the equivalent impedance of the smoothing reactor, while the symbol (Z 2 (jω)) represents the equivalent impedance of the DC filter. The expression defining the line boundary transfer function on the Kunbei side line is as follows : According to Eqs 1-3, the expression of the line boundary transfer function on the north side of Kunming can be further obtained as: In order to quantitatively analyze the characteristics of G 1 (jω), the amplitude frequency characteristics of the boundary transfer function on the Kunbei side line can be obtained by substituting the parameters of each element in Table 1, as shown in Figure 3.
It can be seen from Figure 3 that when 0 Hz < f ≤ 32 Hz |G 1 (jω)| = 1, and it can be considered that the transient current has not decayed. When f > 1472 H, |G 1 (jω)| ≪ 1. It can be seen that the line boundary on the Kunbei side line has a strong attenuation effect on the high-frequency component of the transient current.

The Analysis of Boundary Frequency Characteristics on the Longmen Side
The MMC is used in the inverter side of the Kun-Liu-Long line UHV three-terminal hybrid DC transmission system. Because MMC uses the sub-module cascade to output multilevel stepped waves and the output wave form is of high quality, then only one current-limiting reactor is required at the end of the DC line L 2 (Shen et al., 2020b. Because there is a ground capacitance between overhead lines and the ground, the current-limiting reactor and the ground capacitance of an overhead line at the end of the DC line L 2 constitute the realistic boundary on the Longmen side. The boundary equivalent circuit on the Longmen side is shown in Figure 4. In Figure 4, the symbol (I 3 ) represents the transient current outside the protection zone of the DC line L 2 , while the symbol (I 4 ) represents the transient current flowing into the DC line L 2 after being attenuated by the line boundary. The expression for defining the line boundary transfer function on the Longmen side is as follows: In the formula, Z 3 (jω) jωL gd and Z 4 (jω) 1/jωC L . By substituting the element parameters L gd = 150 mH and C L = 0.00124 µF, the amplitude frequency characteristics of the line boundary transfer function on the Longmen side can be obtained, as shown in Figure 5.
It can be seen from Figure 5 that when 0Hz＜f ≤ 434 Hz |G 2 (jω)| = 1, it can be considered that the transient current has not decayed. When f > 4,248 Hz, |G 2 (jω)| ≪ 1. Therefore, the line boundary on the Longmen side has a strong attenuation effect on the high-frequency component of the transient current.
From the aforementioned analysis of amplitude frequency characteristics of line boundary transfer functions on rectifier and inverter sides, it can be seen that high-frequency components of the transient current passing through boundary elements have obvious attenuation effect. When the fault occurs outside the DC line zone, the high-frequency energy of the fault transient current detected by the protection devices M1 and M2 is small due to the blocking effect of the line boundary. When the fault occurs in the DC line zone, the high-frequency signal is transmitted to the protection installation without obstruction, and the highfrequency energy of the fault transient current detected by the protection devices M1 and M2 is large. Therefore, the attenuation effect of high-frequency components of the transient current can be fully utilized, and the faults inside and outside the zone can be effectively distinguished.

The Complementary Ensemble Empirical Mode Decomposition
Because a certain wavelet base needs to be selected in the wavelet transform, the selection of the wavelet base has a great influence on the whole wavelet analysis result , Yang et al., 2019a. Once the wavelet base is determined, it cannot be replaced in the whole analysis process. Even though the wavelet base may be the best in the whole situation, it may  have poor effect in some parts (Yang et al., 2019b. In this article, complementary ensemble empirical mode decomposition (CEEMD) is selected as the signal processing method. CEEMD decomposes the line mode component of the fault transient current into several intrinsic mode functions (IMFs) and a residual margin according to the fluctuation or trend of different scales, and the intrinsic mode function reflects the fluctuation characteristics of the fault transient current in different time scales (Yang et al., 2021c, Yang et al., 2021b. The residual margin reflects the long-term trend characteristics of the fault transient current (Yeh et al., 2010). CEEMD does not need to predetermine or force a given basis function in advance but depends on the characteristics of the signal itself to decompose adaptively. The IMF components obtained generally have an obvious physical meaning, which is very suitable for the feature extraction of transient signals (Zhao, 2004).
Therefore, CEEMD is used to decompose the line mode component of the fault transient current inside and outside the DC line zone, and the decomposition result is shown in Eq. 6 (Ye and Liu, 2011): In the formula, the symbol x(t) represents the original signal, c j (t) represents the intrinsic mode function, and r n (t) represents the residual term. c 1 (t), c 2 (t), . . ., and c n (t) are obtained after multiple decomposition of x(t), which have only a single frequency component at any time. They respectively represent n intrinsic mode function components which decrease from high frequency to low frequency, and the changing trend of r n (t) is similar to that of the original signal.

The Intrinsic Mode Energy Entropy
The intrinsic mode energy entropy is a measure that reflects the degree of energy distribution disorder of signals in various frequency bands (Zhu et al., 2020). It can be used to extract fault characteristics of the transient current. For the fault transient current, after CEEMD decomposition, the intrinsic mode energy of each frequency band in the decomposition scale on the i layer can be obtained as E c E c1 , E c2 , . . . , E ci , and the energy of the i layer is defined as: By dividing the intrinsic mode function of each layer obtained by decomposition into n parts on the time axis and combining it with information entropy, we can define the intrinsic mode energy entropy H as (Ning et al., 2017): In the formula, p i E ci /(E rn + n i 1 E ci ) (i = 1, 2, . . . ,n); p n+1 E rn /(E rn + n i 1 E ci ).
The intrinsic mode energy entropy effectively combines complementary ensemble empirical mode decomposition, signal energy, and information entropy and effectively quantifies the discrete degree of energy distribution of the fault transient current with frequency (Guan and Zhang, 2011). It can be seen from Eq. 8 that the more uneven the distribution of energy of the fault transient current with frequency is, the smaller the intrinsic mode energy entropy is, and the greater the energy entropy is.

The Criterion to Initiate the Protection
When the DC line L 1 or the DC line L 2 fails, the protection devices on both sides of the T-zone will detect larger transient current (Δi) and voltage (Δu); so, the amplitude of the transient power can be used to initiate the criterion of line protection. Amid the actual operation of the system, the current oscillation generally is allowed to be less than 10% of its rated value, and the voltage oscillation, less than 20% of its rated value. In this context, a certain margin should be reserved for the starting value of the protection device. In order to ensure the anti-interference performance of the protection device, the protection device operates after continuously detecting that the starting values of three data points are greater than the setting values since there are faults. The fault protection device at any position in the zone can be ensured to initiate. For that to happen, the reliability coefficient of the system, according to the analysis of the fault simulation data, is equal to 1.2, and the setting value is equal to 0.0028. By doing so, the starting is shown in the formula as follows (Zhang et al., 2016):

The Criterion to Identify the Direction of the Fault
As there is the T-zone in the Kun-Liu-Long UHV three-terminal hybrid DC transmission system, protection devices are installed on the left and right sides of the zone; thus, directional components need to be added to identify where the faults are in an effective way. In this article, the power directional element as well as the positive and negative characteristics of the fault transient power is applied to identify the direction of faults. We select the time window of 2 ms after the protection is initiated, and we sum the transient power in the time window to form the transient energy. The formula is shown as follows: The The transient energy detected by the power directional element at the protection device (M1) is negative, while the one detected by the power directional element at the protection device (M2) is positive. When the grounding fault emerges from the line on the right side of the T-zone, the transient energy detected at the protection device (M1) is positive, but the one detected at the protection device (M2) is negative. Therefore, the specific discrimination results of the fault direction are as follows: In the formula mentioned previously, the symbol (ΔE M1 ) represents the transient energy detected by the power directional element at the protection device (M1), while the symbol (ΔE M2 ) refers to the transient energy detected by the power direction element at the protection device (M2). As a result, the line mode component of the fault transient current is decomposed by CEEMD (the Complementary Ensemble Empirical Mode Decomposition) before the intrinsic mode energy entropy of each frequency band is calculated, respectively. The energy entropy can be obtained through Eq. 8. When a fault occurs outside the DC line L 1 , the intrinsic mode energy entropy of the fault transient current is smaller. However, when a fault occurs in the DC line L 2 , the intrinsic mode energy entropy of the fault transient current is larger. In order to ensure the protection of the DC line L 1 , the minimum intrinsic mode energy entropy of faults in the zone is supposed to be greater than the maximum one of faults outside the zone. Thus, the criterion to identify the faults, inside and outside, on the left side of the T-zone is as follows:

The Criterion to Identify Faults Inside and Outside the Left Side of the T-Zone
In the formula mentioned previously, the letter N represents the number of sampling points, collected within the 2-ms time window after the protection is initiated, the letter K stands for the number of sampling points at the current moment, the symbol K rel refers to the reliable coefficient, and the symbol H set1 is the threshold value of the intrinsic mode energy entropy of the protection device (M1). Given the interference of the signal noise and sampling errors, the reliability coefficient is equivalent to 1.3. In order to distinguish the faults inside and outside of the zone in a reliable manner, the value of H set1 is selected with the most serious situation taken into account. The setting principle is that the high-resistance grounding fault in the DC line avoids the metallic grounding fault outside the line. The analysis of fault simulation data shows that, taking a certain margin into consideration, the threshold value of the intrinsic mode energy entropy is reasonably equivalent to 0.0008.

The Criterion to Identify Faults Inside and Outside the Right Side of the T-Zone
According to the frequency characteristic analysis of the boundary of the Longmen side line in Section 2.2, the boundary has a strong attenuation effect on the highfrequency component of the transient current. In the same way mentioned previously, the fault transient current line mode component on the right side of the T-zone is decomposed by CEEMD, thus calculating its intrinsic mode energy entropy. As there are faults emerging in and out of the DC line L 2 , the intrinsic mode energy entropy detected by the protection device (M2) is used to distinguish the faults in and out of the DC-zone. The criterion to identify faults, inside and outside, in the right area of the T-zone is expressed in the following: In the formula mentioned previously, the letter N represents the number of sampling points, collected within 2-ms time window after the protection is initiated, the letter K stands for the number of sampling points at the current moment, the symbol K rel refers to the reliable coefficient, and the sign H set2 means the threshold value of the intrinsic mode energy entropy of the protection device (M2). Taking the interference of the signal noise and sampling errors into consideration, the reliable coefficient is equivalent to 1.3. In order to distinguish the faults inside and outside of the zone in a reliable manner, the value of H set2 is selected with the most serious situation taken into account. The setting principle is that the high-resistance grounding fault in the DC line L 2 avoids the metallic grounding fault outside the line. The analysis of fault simulation data shows that, reflecting a certain margin, the threshold value of the natural modal energy entropy is reasonably equivalent to 0.001.

The Criterion of Fault Pole Selection
Normally, the positive and negative poles of the DC power transmission system operate symmetrically. When a fault occurs in the line, a protection device is needed to quickly and accurately identify the fault line as a way of ensuring that the non-fault pole line can normally transmit power without interference from it. Due to the electromagnetic coupling between the positive and negative poles of the line, the non-fault pole will produce an induced current when the grounding fault of the single pole occurs in the DC line. Amid the increase of the high-frequency component of the fault transient, the coupling effect between line poles is enhanced, narrowing the divide of high-frequency components between non-fault poles and fault ones. Therefore, the low-frequency component of the fault transient must be the option for the criterion of fault pole selection (Chu et al., 2017).
In this article, the fault transient power signal with a frequency below 1,000 Hz is extracted in a time window of 2 ms after the protection is operated. Also, the variation of the positive and negative transient power is viewed as the criterion of pole selection, and it is expressed in the following: In the formula mentioned previously, the letter N represents the number of sampling points, collected within the time window of 2 ms after the protection is initiated (the letter i is equal to 1, 2, and so on.). The symbols ΔP 1 (i) and ΔP 2 (i) stand for the amplitude variation of the fault transient power with positive and negative frequency below 1,000 Hz, respectively.
The analysis of fault simulation data suggests that, in order to make the protection devices against faults at any position in the zone, they can accurately select the poles; so, the setting value (K set1 ) of the fault positive pole is equivalent to 150, while the value of the fault negative pole is equivalent to 60. Comparing the ratio (K) of the positive transient power variation to the negative one with both setting values K set1 and K set2 , the result of the criterion is as follows: Positive pole fauls K < K set2 Negative pole fauls K set2 ≤ K ≤ K set1 Faults between poles . (15)

The Protection Scheme
According to the aforementioned analysis, the transient scheme to protect the UHV multiterminal hybrid DC transmission line based on intrinsic mode energy entropy is introduced, by virtue of the criterion to initiate the protection, identify the direction of faults, to protect the criterion of fault pole selection inside and outside the zone. The scheme is shown in Figure 6. First, the protection device detects the transient current (Δi) and voltage (Δu) of the DC lines on both sides of the T-zone, which are substituted into Eq. 9 to calculate the transient power. As the result, when (ΔP star ) is greater than its setting value (ΔP set ), the protection device works. By doing so, the power directional element is used to further identify the decomposed by CEEMD as a way of calculating its intrinsic mode energy entropy. When the minimum intrinsic mode energy entropy of potential faults in DC line L 2 is greater than the threshold value (H set2 ), the right side of the T-zone is dysfunctional. That makes the protection device M2 to operate. On the contrary, when the maximum intrinsic mode energy entropy of the fault outside the DC line L 2 is less than the threshold value (H set2 ), the fault occurs outside the right side of the T-zones. In this context, the protection device M2 does not operate. In the final analysis, the low-frequency component of the fault transient power is extracted to identify in which pole the malfunction appears through the ratio (K) of the positive transient power variation to the negative one.

The Simulation Model
In light of the engineering parameters of the Kun-Liu-Long line UHV three-terminal hybrid DC transmission system, the simulation model of the system as shown in Figure 1 is established in the simulation platform PSCAD/EMTDC (Power Systems Computer-Aided Design/Electromagnetic Transients including DC). The main parameters of the simulation model are shown in Table 2.
In the Kunbei Converter Station, the main connection mode, featured by two unipolar 12-pulse converter units in series, is adopted. Each converter unit bears 400 kV voltage, and the unipolar series voltage is distributed according to (400 plus 400) kV. The constant current control is handled. The Liubei and Longmen converter stations are connected in series with two MMC units in a single pole. The two units are connected in series to form high-and low-valve groups, and the single-pole series voltage is distributed according to (400 plus 400) kV. Each converter chain unit inside the MMC applies the main connection mode that is mixing and cascading half-bridge and full-bridge in proportion. The upper and lower bridge arms are connected to 200 sub-modules at any time, with some left behind. The output voltage of a single converter valve is maintained at 400 kV.
In Figure 1, the total length of the DC line L 1 is 932 km, with an average soil resistivity of 1,750 along the line, whereas the total length of the DC line L 2 is 554 km, with an average soil resistivity of 2,500 along the line. Those where both adopt the frequency-dependent (Phase) model options f 1 , f 3 , f 5 , and f 7 stand for the location of faults on the positive and negative DC line L 1 and the DC line L 2 , the location of faults falls outside the positive and negative DC lines. The protection range of the protection device M1 includes the whole line L 1 , whereas the protection device M2 responses to the whole line L 2 .

The Simulation Verification
Under the amplitude frequency characteristics of the boundary transfer function gained from Sections 2.1, 2.2, in this article, the sampling frequency of 50 kHz is adopted, and the time window of 2 ms is taken after the protection is started. The fault appears at 1.1 s and lasts for 0.1 s.
The faults, inside and outside of the zone, can be distinguished by the value of intrinsic mode energy entropy in an effective way. For that to happen, the grounding fault at the point f 1 , 932 km away from the protection device on the left side of the T-zone, and another one at the point f 2 , outside the DC line L 1 on the left side of the T-zone, are viewed as typical examples in this article. The intrinsic mode energy entropy that is calculated is made to slide in the sampling window of 2 ms. By doing so, the curve, showed in Figure 7, in which the intrinsic mode energy entropy changes with time, is accessible. From Figure 7, the intrinsic mode energy entropy responses to the overall amplitude of the signal in a less sensitive manner. If faults occur inside and outside the DC line L 1 , the intrinsic mode energy entropy at the position where the protection devices are installed will suddenly change. Due to the attenuation effect of the line boundary on the high-frequency component of the fault transient, the intrinsic mode energy entropy in the DC line is obviously larger than that of outside the line.
With regard to the transient protection scheme of the UHV multiterminal hybrid DC transmission line based on intrinsic mode energy entropy, the simulation verification of the protective principle is conducted, taking different positions of the fault anode and cathode (the points f 2 and f 6 outside the DC line L 1 , the points f 1 and f 5 at 932, 466, and 0 km away from the protection device on the left of the T-zone, the points f 3 and f 7 at 0, 277, and 554 km away from the protection device on the right side of the T-zone, and the points f 4 and f 8 outside the DC line L 2 ) as well as various transition resistances (0.01 Ω, 300 Ω, and 500 Ω) into account. The setting values of the protection are shown as follows: K set1 is equivalent to 150, K set2 , 60, H set1 , 0.0008, and H set2 , 0.001.
From the simulation results in Tables 3, 4, the location of faults and the transition resistance are expanding, and the energy entropy of the natural mode is decreasing. The entropy values, however, of faults inside and outside the zone are still quite different. That said, the protection is not affected and resistant to the transition resistance to some extent. In this article, the transient protection scheme of the UHV multiterminal hybrid DC transmission line based on the intrinsic mode energy entropy boasts a good performance on protection and high reliability. That means the protection devices (M1) and (M2) can operate in a reliable way as they are in face of different fault poles, positions, and transition resistances.

CONCLUSION
Regarding to the analysis and summary of boundary frequency characteristics of the Kun-Liu-Long UHV three-terminal hybrid DC transmission line, and on the basis of the complementary ensemble empirical mode decomposition (CEEMD), the intrinsic mode energy entropy, and the transient power polarity, a transient protection scheme, suitable for the UHV multiterminal hybrid DC transmission system, is put forward in this article. The protection scheme handles the positive and negative characteristics of transient energy of the fault power detected by power directional elements on both sides of the T-zone to identify the fault direction. As the boundary of the fault line is capable of reducing the high-frequency component of the transient current in a strong manner, the value of the intrinsic mode energy entropy is applied to further distinguish the faults inside and outside the zone. After the analysis mentioned previously, the following conclusions are attainable: 1) The transient current high-frequency component has obviously attenuated as it passes through the line boundary element, but the potential of the attenuation is supposed to be fully tapped to effectively distinguish the faults in and out of the zone.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.