Preprint
Article

Hybrid Multi-Infeed Receiver Line Longitudinal Protection Scheme Based on Voltage Waveform Comprehensive Distance Similarity

Altmetrics

Downloads

67

Views

12

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

09 January 2024

Posted:

09 January 2024

You are already at the latest version

Alerts
Abstract
The coupling of AC and DC power will impact the protective actions on the AC side and pose a threat to the stable operation of the interconnection system. Therefore, a new longitudinal protection method is proposed based on the comprehensive distance similarity of voltage waveforms. Initially, the measured voltage and current data are extracted to calculate the reference voltage, and the voltage waveform fitting is optimized. Subsequently, the Euclidean dynamic time warp (DTW) distance and entropy weight method are utilized to process the voltage waveform, enabling the calculation of its comprehensive distance similarity. This similarity is adopted to determine fault location. A hybrid DC multi-feed AC/DC interconnection system, incorporating line commutated converter-voltage source converter (LCC-VSC) and line commutated converter-modular multilevel converter (LCC-MMC), was established in PSCAD, and fault data was simulated and output. The effectiveness of the protection scheme was validated using MATLAB. Simulation results demonstrate that the proposed method can accurately distinguish between faults inside and outside a region. When compared to existing protection methods, it demonstrates superior performance in resisting transition resistance and noise interference, while also mitigating the impact of data asynchronicity. The speed and reliability of the method are further enhanced.
Keywords: 
Subject: Engineering  -   Electrical and Electronic Engineering

1. Introduction

With the continuous development of China's socio-economic landscape, the electrical demand is rapidly increasing. In order to balance the rising demand for electricity and resource allocation, China is accelerating the development of the strategies of "west to east power transmission " and " north to south power transmission ". High voltage direct current (HVDC) transmission technology is being utilized in grid development to meet the inevitable requirements of long distances, high capacity, and high voltage levels. [1,2].
Line commutated converter (LCC) is suitable for long-distance and high-capacity scenarios. However, due to its semiconductor control, it suffers from the drawback of commutation failure and is generally used at the rectifier end. Conversely, voltage source converter (VSC), employing fully controlled devices, can resist commutation failure on the inverter side and flexibly control the active and reactive power of the system. Its application in multi-feed systems can enhance operational flexibility at the receiving end and expand transmission capacity [3,4]. However, the application of VSC as a fully controlled device means a higher economic investment Therefore, the combination of LCC and VSC in multi-infeed DC transmission harnesses the advantages of both to compensate for their respective shortcomings. This becomes especially advantageous in aspects such as multi-AC grid interconnection, multiple power supply, and multiple load connection, fully leveraging the economic and flexible nature of DC transmission. However, the introduction of multi-DC lines for feed-in complicates the system and makes protection prone to false operations [5]. Hence, research is necessary for the relay protection of the receiving end AC system.
Currently, the protection methods for AC transmission lines primarily rely on traditional relay protection principles, including directional protection, distance protection, and differential protection [6]. In comparison to traditional AC transmission, the fault characteristics for AC lines in hybrid multi-infeed DC systems change, prompting scholars to research the applicability of traditional relay protection [7,8,9,10,11,12,13,14]. For instance, Hao et al. [7] discovered that due to the fact that the protection unit near the inverter side can only identify continuous commutation failures occurring after a fault, it is ineffective at identifying faults where there is no commutation failure after the fault or when there is only one commutation failure. Therefore, when a single transmission corridor structure is used on the AC output side of the inverter side, traditional directional protection is not suitable for AC/DC hybrid interconnected networks. Additionally, continuous commutation failures might lead to severe faults on the AC side of multi-infeed DC systems, causing a sharp increase in the reactive power demand for system fault recovery [8,9,10,11]. This can lead to transient power reversals in neighboring multiple AC lines and the DC system [12,13], while the reverse power flow can lead to potential misoperation of the longitudinal directional protection [14]. Meanwhile, Huang et al. [15] proposed improvements to the protection setting principles of the second zone of distance protection based on fault characteristics, but were unable to address errors caused by high-order harmonics. Additionally, Huang et al. [16] presented a rapid fault detection method based on instantaneous active power and fault transient components to alleviate the effects of commutation failures and reactive power compensation. Ma et al. [17] used time-domain differential equations to calculate fault point location, enhancing the protective performance of fault distance measurement, but with poor anti-interference capabilities that are unable to protect the entire length of the line. Moreover, Zheng et al. [18] analyzed the impact of non-periodic components and harmonic components caused by DC feeding on the accuracy of algorithms, resulting in decreased protective operation performance. In summary, the multi-DC system feeding complicates the research for AC line protection at the receiving end of hybrid AC/DC systems, currently highlighting the following limitations: (1) Traditional directional and distance protection strategies are challenging to ensure the reliability of AC protection, and existing protection schemes still have shortcomings. (2) Higher harmonics on the DC side can also affect the reliability of AC protection operations. (3) Currently, there is limited research and analysis on the protection at the AC receiving end, with most of it being qualitative analysis.
This article introduces a new protection principle based on comprehensive voltage waveform analysis for multi-infeed AC/DC hybrid interconnected systems. The key contributions are as follows:
(1)
A novel protection method is proposed that relies on the similarity of voltage waveforms between measuring and reference points. It addresses the issue of unreliable operation of AC protection in AC/DC hybrid systems.
(2)
The Euclidean distance and dynamic time warping (DTW) distance are investigated to assess waveform similarity, thus avoiding the impact of data asynchrony or protection misoperation caused by a single distance metric.
(3)
The least squares method (LSM) is employed to fit the voltage waveform, effectively mitigating interference from high-order harmonics.
The article is structured as follows: Section 2 introduces the measurement method for voltage waveform distance similarity. Section 3 proposes the topology of the hybrid AC/DC multi-infeed interconnected system and analyzes the system's fault characteristics, then describes the protection scheme steps and processes. In Section 4, a model is established in PSCAD, and MATLAB is used to simulate and verify the feasibility of the protection scheme. The conclusion is presented in Section 5.

2. Comprehensive Voltage Waveform Similarity Calculation

2.1. Waveform Processing

In the event of an AC fault in the AC/DC hybrid network, the fault phase voltage amplitude rapidly decreases. Additionally, due to the introduction of high-order harmonics from DC injections, waveform preprocessing is necessary before comparing the similarity of the measured and reference voltage waveforms, in order to more intuitively reflect the trend of change. As for the high-order harmonic interference, common engineering solutions currently include wavelet transform, which suffers from drawbacks such as high computational complexity, poor speed, and low accuracy [19]. Filtering methods like amplitude limiting and median filtering perform poorly in handling periodic disturbances in AC/DC power systems and are prone to eliminating crucial waveform information [20]. Therefore, waveform signal processing must ensure both rapid response and the retention of waveform information that reflects trends. In this study, the use of the LSM to fit voltage waveforms meets both of these requirements.
LSM is a mathematical fitting technique that connects discrete points on a plane through a curve. It generally assumes that the curve can be represented by:
f ( x ) = a 0 + a 1 x + a 2 x 2 + a 3 x 3 + + a N x b
he unknown coefficients are denoted as a0, a1, a2aN, and b in xb represents the power to which we want to fit.
Assuming the voltage waveform is ( x i , y i ) , i = ( 1 , 2 , 3 , , m ) , where m denotes the number of voltage values, the fitted curve is then given by:
f ( x ) = min i m [ f ( x i ) y i ] 2
By fitting the measured voltage and reference voltage waveforms, high-order harmonic interference can be eliminated while retaining the variation trend.

2.2. Euclidean Distance Calculation

Euclidean distance, also known as Euclidean metric, is commonly used to measure the true distance between two points in a space, reflecting the overall distribution characteristics of waveform similarity [21]. The instantaneous values of measured voltages at the protection installation points M and N are extracted. Corresponding reference voltages are calculated using Equation (12), defining the sample time sequences of the measured voltages and the reference voltages as dM, dN, drM, and drN, respectively. Here, d M = ( x 1 , x 2 , x 3 , , x n ) d N = ( y 1 , y 2 , y 3 , , y n ) d r M = ( z 1 , z 2 , z 3 , , z n ) d r N = ( q 1 , q 2 , q 3 , , q n ) , n represents the number of sample points, chosen based on the precision requirements of the protection method. The Euclidean distance of the voltage waveforms between the measured voltages and the corresponding reference voltages at the protection installation points on both sides is given by:
D W M ( d M , d r M ) = i = 1 n ( x i z i ) 2
D W N ( d N , d r N ) = i = 1 n ( y i q i ) 2
where DWM and DWN represent the Euclidean distance between the measured voltage and reference voltage waveforms at the installation points on sides M and N, respectively.

2.3. DTW Distance

The DTW method employs the concept of dynamic programming, adjusting the different sample points of sequences through compression and expansion, thus enabling a point in one sequence to match multiple points in another, optimizing the warped path. Therefore, when the DTW algorithm computes distance, it doesn't require a one-to-one correspondence of values based on the same sample points, unlike the calculation method for Euclidean distance. The DTW algorithm can address issues of sequence asynchrony, better capturing the overall dynamic characteristics of curves, and measuring the overall similarity of shapes between sequences [22]. For the sampled sequences dM and drM, an n×m distance matrix DM is constructed, where the elements are DM(i, j) as per Equation (5). Similarly, for the sampled sequences dN and drN, an n×m distance matrix DN is built, with elements being DN(i, j) as per Equation (4).
D M ( i , j ) = ( x i z j ) 2 1 i n , 1 j m
D N ( i , j ) = ( y i q j ) 2 1 i n , 1 j m
where Equation (5) represents the Euclidean distance between the sequences (dM) and (drM) at the sequence points (xi) and (zj), while Equation (6) represents the Euclidean distance between the sequences (dN) and (drN) at the sequence points (yi) and (qj).
The set composed of each pair of adjacent elements from DM(i, j) and DN(i, j) is defined as the bent path of the adopted waveform, denoted by L. Taking path L as an example, the element ls refers to the coordinates of the s-th point on the path, i.e., ls = (i, j). Path L must satisfy the following two-part constraints:
(I)The selected path must include all sampling points.
(II)Each sample point needs to be matched with the adjacent sample point, such that if ls=(i, j), ls+1=(a, b)satisfies { 0 a i 1 0 b j 1 .
This yields the DTW distance between the sequences dM and drM, and the DTW distance between the sequences dN and drN:
D M d t w ( d M , d r M ) = min L s = 1 k D ( L s )
D N d t w ( d N , d r N ) = min R s = 1 k D ( R s )
where D(Ls) and D(Rs) represent the cumulative distance of the bent path L and R, while DMdtw and DNdtw respectively denote the waveform DTW distance between the measured voltage and the reference voltage at the installation locations of the protections on both sides.

2.4. Local Dynamic Mapping

In order to better describe the local dynamic characteristics of the voltage waveform, reflecting its trend information at each sampling interval and fully considering the accuracy of the algorithm, the discrete sampling sequences d M d N d r M and d r N are transformed into shape sequences of length d M d N d r M d r N 。Taking the element xʹi as an example, as shown in Equation (9), where Δt represents the time interval:
x i = x i + 1 x i Δ t i = 1 , 2 , , n 1

2.5. Comprehensive Voltage Waveform Similarity Description

This chapter combines the Euclidean distance with the DTWdistance to propose a measurement method for the comprehensive voltage waveform similarity that encompasses overall, local, and dynamic characteristics.
D M φ = α D W M ( d M , d r M ) + β D M d t w ( d M , d r M ) + γ D M d t w ( d M , d r M )
D N φ = α D W N ( d N , d r N ) + β D N d t w ( d N , d r N ) + γ D N d t w ( d N , d r N )
where DM𝜙 and DN𝜙 represent the comprehensive distance similarity of the measured voltage and reference voltage waveforms at the protection installation points on M and N sides, where φ = A, B, C represents the phase sequence index. The smaller the values of D and D, the higher the similarity of the voltage waveforms. DWM(dM, drM) and DWN(dN, drN) reflect the overall distribution characteristics of the voltage waveforms, while DMdtw(dM, drM) and DNdtw(dN, drN) reflect the overall dynamic characteristics of the voltage waveforms. Additionally, DMdtw(dʹM, dʹrM) and DNdtw(dʹN, dʹrN) reflect the local dynamic characteristics of the voltage waveforms. α, β, and γ are the weights assigned to each distance measure, and based on the actual sampling accuracy requirements and the principle of entropy weighting, their respective values are set to 0.4, 0.3, and 0.3 [22].

3. System Fault Characteristics and Protection Process

3.1. Hybrid DC Multi-Infeed AC/DC Interconnected System Topology Structure

The constructed hybrid AC/DC interconnected system with multiple DC inputs, as shown in Fig. 1, features the transmitting end's DC section consisting of two distinct ±500kV hybrid DC transmission lines: a bipolar double-end line commutated converter-modular multilevel converter (LCC-MMC) transmission system and a pseudo-bipolar double-end LCC-VSC transmission system. The receiving end's AC section comprises a 230kV dual-source system, with M and N denoting the positions of the protection units on both sides of the line. rM represents the calculation reference point at end M, and rN represents the calculation reference point at end N, defining the protection scope from the protection installation point to the reference point on that side. The overlap of the protection scopes at both ends indicates the range for determining internal faults, denoted by the green area in Fig. 1. This study focuses primarily on the central section, MN, where f1 and f5 represent faults occurring outside the protection scope, f2 represents a fault occurring solely within the protection scope at end N but outside the reverse area, f3 represents a fault occurring within the protection scope at both ends M and N, and f4 represents a fault occurring solely within the protection scope at end M but outside the forward area.
In an LCC-VSC system, the rectifier side comprises two sets of twelve-pulse LCC converters connected in series. The rectifier side direct current is input into the constant current control circuit, which controls the α angle to regulate the rectifier-side converters. On the inverter side, low-cost and highly stable two-level VSC modules are employed. Through vector control, the modulated waves are input into the pulse width modulation (PWM) controller, which then produces the control signals.
In an LCC-MMC system, the rectifier side is similar to that of the LCC-VSC system, whereas the inverter side consists of two MMC modules connected in series. Input signals related to inverter-side power and current are sent to the inner and outer loop control circuits, and the output signals are directed to the bridge arms to control the inverter-side converters.
Figure 1. Topology diagram of the hybrid DC multi-feed system.
Figure 1. Topology diagram of the hybrid DC multi-feed system.
Preprints 95869 g001

3.2. Reference Voltage Calculation and System Fault Characteristics

The similarity of the voltage waveforms between the measuring point and the reference point is calculated to determine if a fault has occurred, and the position of the reference point needs to be determined. The reference point is usually located outside the protection line to ensure sufficient sensitivity. Additionally, to reduce errors, the reference point should not be too far away; scholars both domestically and abroad typically position it within the range of 1.2 to 2.0 times the length of the protected line [23]. After extensive simulation experiments, it has been found that placing the reference point at 1.3 times the length of the line can meet the experimental requirements outlined in this paper, while allowing for flexible adjustments within the range based on specific protection needs and interference considerations.
Compared to traditional purely AC lines, in a multi-infeed interconnected grid, significant non-periodic components are introduced into the DC system after a fault occurs. This results in a more complex variation of electrical quantities in the AC lines. Obtaining accurate voltage phasors through traditional Fourier algorithms no longer guarantees accuracy. Instead, this paper proposes using instantaneous measurement methods to calculate the reference voltage, taking the M-side as an example.
In Fig. 2 (a), ZSM and ZSN represent the equivalent system impedances on both sides of the AC line, Im and Um represent the instantaneous measured currents and voltages, and Ir and Ur represent the instantaneous reference currents and voltages. Due to the rectifying sides of the hybrid multi-infeed DC system using constant current control, and the inverting side voltage being determined by the receiving-end AC bus, both can be considered equivalent to a voltage-controlled current source controlled by the receiving end. ΔIdc.MMC and ΔIdc.VSC represent the equivalent fault component currents for the LCC-MMC and LCC-VSC hybrid DC systems, respectively.
Based on the situation of the completed and ongoing multi-infeed AC/DC hybrid power grid projects, this paper sets a total line length of 160 km, using an equivalent π-type line model considering the electric shunt, as depicted in Fig. 2(b). After measuring the instantaneous current value Im and the instantaneous voltage Um, the reference voltage Ur can be obtained through the differential equation in Equation (12).
U r = U m [ R ( I m C 2 d U m d t ) + L d ( I m C 2 d U m d t ) d t ]
Figure 2. Schematic diagram of system and line model (a) Simplify the system model; (b) Schematic diagram of system and line model.
Figure 2. Schematic diagram of system and line model (a) Simplify the system model; (b) Schematic diagram of system and line model.
Preprints 95869 g002
Under normal operation, the voltage waveforms at the measurement point and reference point are similar. To analyze the fault characteristics of the system at different locations, we take the occurrence of a ground fault as an example and analyze the measured voltage Um and the reference voltage Ur. Assuming the fault occurs at location k1 as illustrated in Fig. 2(a), which represents an external fault in the opposite direction of MrM, as shown in Fig. 3:
Figure 3. Schematic diagram of out-of-reverse fault system.
Figure 3. Schematic diagram of out-of-reverse fault system.
Preprints 95869 g003
At this point in time, the actual current direction is opposite to the measured current direction during normal operation, hence Equation (12) can be rewritten as:
U m = [ R f ( I m C f 2 d U m d t ) + L f d ( I m C f 2 d U m d t ) d t ]
U r = [ I m ( R f R r ) d U m d t ( R f C f 2 R r C r 2 ) + L f d ( I m C f 2 d U m d t ) d t L r d ( I m C r 2 d U m d t ) d t ]
where Rʹf, Lʹf, and Cʹfrepresent the resistance, inductance, and capacitance values from measurement point M to fault point k1, respectively. Rʹr, Lʹr, and Cʹr represent the resistance, inductance, and capacitance values from measurement point M to reference point rM, respectively. By comparing the two equations, it can be concluded that although there are differences in the amplitude of the measured voltage Uʹm and the reference voltage Uʹr, their variation trends are consistent. The computational relationship deduced from Equation (14) is identical to that of Equation (12), indicating that when an external fault occurs in the reverse direction area, the measured voltage and the reference voltage waveforms are similar.
Assuming that the fault occurs at position k3 as shown in Fig. 2 (a), i.e., an external fault in the positive direction of MrM, as illustrated in Fig. 4:
Figure 4. Schematic diagram of out-of-zone fault system.
Figure 4. Schematic diagram of out-of-zone fault system.
Preprints 95869 g004
At this time, the actual current direction is the same as the measured current direction under normal operation, rewriting Equation (1) as follows:
U m = R f ( I m C f 2 d U m d t ) + L f d ( I m C f 2 d U m d t ) d t
U r = I m ( R f R r ) d U m d t ( R f C f 2 R r C r 2 ) + L f d ( I m C f 2 d U m d t ) d t L r d ( I m C r 2 d U m d t ) d t
where Rʺf, Lʺf, and Cʺfrepresent the resistance, inductance, and capacitance values from measurement point M to fault point k3, respectively. Rʺr, Lʺr, and Cʺr respectively represent the resistance, inductance, and capacitance values from measurement point M to reference point rM. By comparing Equations (15) and (16), it can be observed that although there is a difference in the magnitude of the measured voltage Uʺm and the reference voltage Uʺr, their variation follows a consistent pattern. The computational relationship derived from Equation (16) is the same as Equation (12), indicating that when an external fault occurs, the measured voltage and the reference voltage waveforms are similar.
Suppose that the fault has occurred at k2 in the area of MrM, as shown in Fig. 2(a), depicting an internal fault in the MrM area, as illustrated in Fig. 5:
Figure 5. Schematic diagram of fault system in the area.
Figure 5. Schematic diagram of fault system in the area.
Preprints 95869 g005
At this time, the fault branch between Mr and M is increased, and the line structure is disrupted. Based on the characteristic that the distance from measuring point M to reference point rM is greater than to the fault point during an internal fault, Equation (1) is modified as follows:
U m = R f ( I m C f 2 d U m d t ) + L f d ( I m C f 2 d U m d t ) d t
U r = [ R r ( I m C r 2 d U m d t ) + L d ( I m C r 2 d U m d t ) d t ] U m = d U m d t ( R f C f 2 R r C r 2 ) I m ( R f R r ) + L f d ( C f 2 d U m d t I m ) d t L r d ( C f 2 d U m d t I m ) d t
where Rʺʹf, Lʺʹf, and Cʺʹfrespectively represent the resistance value, inductance value, and capacitance value from measuring point M to fault point k2. Rʺʹr, Lʺʹr, and Cʺʹr represent the resistance value, inductance value, and capacitance value from measuring point M to reference point rM. By comparing Equations (17) and (18), it can be inferred that although there are differences in the amplitude of the measured voltage Uʺʹm and the reference voltage Uʺʹr, their variation trends remain consistent. Equation (18) deduces an operational relationship opposite to that of Equation (12), indicating that when an internal fault occurs, the measured voltage and the reference voltage waveforms are dissimilar.
In summary, it can be concluded that during an internal fault, the measured voltage and the reference voltage waveforms are dissimilar, whereas during normal operation and external fault occurrence, the voltage waveforms are similar. This provides a theoretical basis for the subsequent use of voltage comprehensive waveform similarity.

3.3. Protection Action Criteria

By employing the calculation method in Section 2.5, the comprehensive waveform similarity at both ends is obtained. From the computed results for the three phases, the maximum values are denoted as MAX(D) and MAX(D), and a threshold value Dset is established. When the voltage comprehensive waveform similarity exceeds the threshold value, it indicates dissimilarity between the two voltage waveforms. Considering the topological structure and operational conditions of the AC/DC hybrid system, along with practical engineering requirements and a large volume of simulation data, the threshold value Dset is chosen as 1.2. By comparing the relationship between MAX(D), MAX(D), and the protection threshold value Dset, fault locations can be further determined.
When MAX(D) < Dset and MAX(D) > Dset, the fault occurs between protection device M and reference point rM, that is, the position f2 in Fig. 1 represents an external fault.
When MAX(D) > Dset and MAX(D) < Dset, the fault occurs between protection device N and the reference point rN, indicating that the fault location is external, as illustrated in position f4 in Fig. 1.
When MAX(DMφ) < Dset and MAX(D) < Dset, the fault occurs in areas outside protection devices M and N, i.e., positions f1 and f5 in Fig. 1 indicate an external fault zone.
When MAX(D) > Dset and MAX(D) > Dset, the fault is determined to have occurred between the two protective devices, i.e., the location f3 in Fig. 1 indicates an internal fault, triggering the protective action.
When a fault occurs, the fault component of the current will increase or decrease. Based on the magnitude of this sudden change, protection initiation criteria can be established, and the fault component of the current can be extracted using the superposition theorem.
Δ I s ( t ) = I s ( t ) I s ( t T )
where t denotes the sampling instant. Is(t) represents the actual current value at location s at instant t. ΔIs(t) stands for the fault component of the current extracted at location s at instant t for protection measurement. T represents the selected pre-fault sampling instant. Considering the system's operational stability, typically the two sampling cycles before the fault are chosen. Based on the variation in the fault component of the current obtained from Equation (19), the initiation criterion can be obtained as follows:.
Δ I s M φ > 0.01 k I n
Δ I s N φ > 0.01 k I n
where ΔIsMφ and ΔIsNφ respectively denote the change in current fault components at the protection installation points on sides M and N. Here, φ=A, B, C represents the phase sequence index, k is the reliability factor, and In stands for the rated line current. As long as any of the judgment conditions in Equations (20) and (21) are satisfied, the calculation of the comprehensive voltage waveform distance similarity can be initiated to further determine the fault area.
Based on the above fault analysis process and the introduction of the protection scheme, the overall protection process can be illustrated as shown in Fig. 6, with the specific steps outlined below:
(1) Run the system model, debug the protection measurement devices, collect electrical signal data, and extract the sampled values of voltage and current for the M and N ends of the AC transmission line using preset sensors.
(2) Extract transient changes in fault component from the sampled current values at ends M and N, and determine if a fault has occurred using a protection initiation criterion. When the criterion is met, it indicates a fault in the system.
(3) Calculate the reference voltages at both ends based on the sampled voltage and current values, and use the method of least squares to fit the waveform trend of the sampled and calculated values, thereby filtering out high-order harmonic interference.
(4) Derive the comprehensive waveform distance similarity metrics D and D between the measured voltage and the reference voltage at both ends of the line. Compare the maximum values of D and D with the protection set value Dset to determine the fault occurrence region.
When D > Dset and D > Dset, the fault is determined to occur between the two protective devices, indicating an internal fault, located at position f3 in Fig. 1; otherwise, it is determined to be an external fault.

4. Experimental Verification

4.1. Introduction to the Experimental Platform

To verify the theoretical analysis and the feasibility of the protection scheme mentioned above, the PSCAD/EMTDC simulation platform is utilized to construct the hybrid DC multi-infeed system, as depicted in Fig. 1, for power system simulation. The electrical parameters are adjusted, and various fault occurrences are simulated under different conditions by controlling the fault area, line positions, and other interfering factors. Subsequently, the validation of the protection algorithm is carried out in MATLAB.
Figure 6. Schematic diagram of fault system in the area.
Figure 6. Schematic diagram of fault system in the area.
Preprints 95869 g006
The sampling window for calculating the comprehensive voltage distance similarity in this paper is 1ms, and the sampling frequency for line protection is 10 kHz. Therefore, fault detection requires only 10 sampling points. In practical engineering applications, the synchronization error of the commonly used global positioning system (GPS) synchronization method is less than 2 μs. The protection algorithm investigated in this paper calculates the similarity value within the sampling interval and compares to find the maximum value, making it possible to accurately identify the fault position within a 1ms delay, without strict bi-directional timing and communication synchronization.
In conclusion, the simulated duration in this paper is set as 3s, and the fault occurs at 2.25s and lasts for 0.05s. Due to the rapid change in the current fault component, the protection can be activated in about 0.5ms, while the simulation data collection and the calculation of comprehensive voltage distance similarity can be completed within 2ms. Considering a certain margin, the proposed method can operate quickly and reliably within 3ms.

4.2. In-Zone Fault Simulation Result

In the positive direction of the AC line, a phase A-to-ground fault is set at a distance of 50 km. To stabilize the feed-in of the hybrid DC system, the fault is initiated at 2.25s and lasts for 0.05s. Taking the data at the M-end as an example, the simulation compares the waveform differences between the normal and fault states of both a pure AC system and a hybrid multi-feed-in system. The simulation window length is 0.025s, and the simulation results are shown in Fig. 4.
Figure 7. Comparison of measured voltage waveforms between pure AC system and hybrid multi-feed system (a) Comparison of measured voltage waveforms under normal operation; (b) Comparison of measured voltage waveforms under faults in the area.
Figure 7. Comparison of measured voltage waveforms between pure AC system and hybrid multi-feed system (a) Comparison of measured voltage waveforms under normal operation; (b) Comparison of measured voltage waveforms under faults in the area.
Preprints 95869 g007
In Fig. 7, it is evident that after the introduction of multi-feed-in, both the voltage amplitude and phase undergo changes. Additionally, the system, besides experiencing voltage reduction after the fault, is also influenced by higher harmonics, leading to irregular oscillations in the fault waveform. By extracting data from the N-end and computing the reference voltage at both ends, Figure 8 compares the measured voltage and reference voltage waveforms.
Figure 8. Comparison of measured voltage waveforms between pure AC system and hybrid multi-feed system (a) Comparison of measured voltage waveforms under normal operation; (b) Comparison of measured voltage waveforms under faults in area.
Figure 8. Comparison of measured voltage waveforms between pure AC system and hybrid multi-feed system (a) Comparison of measured voltage waveforms under normal operation; (b) Comparison of measured voltage waveforms under faults in area.
Preprints 95869 g008
In Figs 8(a) and 8(b), the measured voltage and reference voltage waveforms at terminals M and N both exhibit an upward trend before the fault occurs, with high waveform similarity. After the fault occurs, the measured voltage continues to rise, while the reference voltage decreases. The magnitudes of voltage fluctuations in both waveforms are quite similar, exhibiting an almost symmetrical relationship around 0kV, albeit with opposite directions. This aligns with the conclusion drawn in Section 1.2.
For irregular variations in fault waveforms in mixed multi-feed systems, directly applying a similarity algorithm focusing on detailed levels can lead to significant errors. However, after applying the waveform processing method in Section 2.1, waveform optimization can be achieved while retaining trend information, as illustrated in Fig. 9.
Figure 9. Comparison of voltage waveforms after treatment (a) Comparison of M-terminal voltage waveforms after optimization; (b) Comparison of N-terminal voltage waveforms after optimization.
Figure 9. Comparison of voltage waveforms after treatment (a) Comparison of M-terminal voltage waveforms after optimization; (b) Comparison of N-terminal voltage waveforms after optimization.
Preprints 95869 g009
As shown in Fig. 9, the processed waveform is smoother and largely retains the trend information of the waveform changes, reducing errors in the subsequent calculation of the comprehensive distance similarity of voltage waveforms.
To study the changes in the comprehensive distance similarity of internal fault waveforms under various operating conditions, faults were set at 50km in the positive direction of lines M and N, with controlled transition resistance size. This involved simulating single-phase-to-ground faults (A-G), two-phase-to-ground faults (BC-G), three-phase-to-ground faults (ABC-G), and phase-to-phase short-circuit faults (BC). The fault initiation time was set at 2.25s, with a duration of 0.05s. Simulation data were tabulated in Table 1.
In Table 1, bolded data indicates values greater than the set value, and data marked in red, such as DMA=4.4298 and DNA=4.6376, represent the maximum values in this dataset. As Table 1 demonstrates, during the occurrence of faults, the comprehensive distance similarity of the faulted phase consistently exceeds the set value Dset, enabling correct action and non-operation for non-faulted phases. All of these are in line with the simulation results and theoretical analysis, meeting the protection requirements.

4.3. Out-of-Zone Fault Simulation Result

Referring to Fig. 1, a single-phase bolted fault (f1) is established 40 km outside the M and N line in the opposite direction. Voltage and current data from both ends of M and N are extracted, the reference voltage is computed, and waveforms are processed. A comparison of the measured voltages at both ends and the reference voltage waveform is shown in Fig. 10.
In Fig. 10, the trends in the measured voltage and the reference voltage waveform are consistent, aligning with the conclusion derived from theoretical analysis concerning the similarity of the two voltage waveforms during an external fault. To further investigate the comprehensive distance similarity characteristics of the voltage waveform during external faults under various operating conditions, with reference to Fig. 1, f1 and f2 are faults set respectively 30 km and 10 km outside the opposite direction of the M, N line, and f4 and f5 are faults set respectively 10 km and 30 km outside the forward direction of the M, N line.
These faults simulate single-phase earth faults, two-phase earth faults, and three-phase earth faults. The fault initiation time is set at 2.25s, with a duration of 0.05s. The simulated data is summarized in Table 2.
Figure 10. Comparison of measured voltage and reference voltage waveform during out-of-zone fault (a) Comparison of M-terminal voltage waveforms; (b) Comparison of N-terminal voltage waveforms.
Figure 10. Comparison of measured voltage and reference voltage waveform during out-of-zone fault (a) Comparison of M-terminal voltage waveforms; (b) Comparison of N-terminal voltage waveforms.
Preprints 95869 g010
Based on theoretical analysis, when an external fault occurs, the measured voltage and the reference voltage waveform are similar. In Table 2, bold data indicates values greater than the set value, while red data represents the maximum value within the dataset. For external faults occurring at f1 and f5, the data in Table 2 are all less than the set value (Dset), thus neither end of the protection devices act falsely. In the case of an f2 external fault, as MAX(D) remains consistently less than the set value Dset, both ends of the protection devices remain unaffected. Similarly, when an f4 external fault occurs, as MAX(D) consistently remains less than the set value Dset, both ends of the protection devices remain unaffected. In conclusion, the simulation results for external faults meet the protection requirements.

4.4. Protection Action Status

Currently, in practical engineering of AC transmission lines, the highest proportion of faults is single-phase-to-ground faults. Taking the occurrence of single-phase high-impedance ground faults as an example, simulated protection data is illustrated in Fig. 11.
Figure 11. Simulation results for internal and external faults identification.
Figure 11. Simulation results for internal and external faults identification.
Preprints 95869 g011
From Fig. 11, it can be visually observed that the operating condition of D > Dset and D > Dset is satisfied only when an internal fault f3 occurs. Selecting the setting value Dset = 1.2 provides a certain margin for the calculation results when the waveforms are similar. The proposed method can identify internal and external faults reliably and can meet the protection requirements.

4.5. Anti-Interference Analysis

In actual engineering scenarios, there exist numerous interferences. To test the method's resistance to interference, this paper introduced a 10dB noise ratio interference under a 100Ω grounding fault and simulated different positions within the zone. These positions were at 10%, 50%, and 90% of the total line length. The simulation data is summarized in Table 3:
In Table 3, bold data indicates values greater than the setting value, and red data represents the maximum value within the dataset. When comparing the interference-added comprehensive distance similarity calculation results at the 50% total line length with the data in Table 2, it is observed that there are varying degrees of increase or decrease in results after the addition of noise interference. For the 10% and 90% total length positions, the comprehensive distance similarity of the fault phase is greater than the setting value Dset, indicating non-rejection of misoperation. Conversely, the comprehensive distance similarity of the non-fault phase is lower than the setting value Dset, demonstrating non-operation. This shows that the proposed method can correctly operate even in the presence of noise interference, meeting the requirements for protection.

4.6. Comparison and Discussion

To demonstrate the superiority of the proposed method, it is compared with [24] and the longitudinal differential protection scheme. [24] uses the conclusion that reference voltage and measured voltage waveforms have different characteristics during faults inside and outside the area to roughly determine the fault location through the Pearson correlation coefficient. Then, adopting the longitudinal protection logic, it provides a more accurate judgment result. However, [24] lacks an analysis of the protection scheme's ability to withstand transition resistance. The longitudinal current differential protection scheme utilizes the assigned values and phases of both sides of the protected line to determine whether a fault has occurred. The criterion is shown in Equation (22).
I d ˙ > K I r ˙
where the differential current I d ˙ = I ˙ M + I ˙ N and the braking current I r ˙ = I ˙ M I ˙ N are defined. I ˙ M I ˙ N represent the measured currents at both ends of M and N, respectively. K is the braking coefficient, and the operating threshold is set as the maximum unbalanced current.
To investigate the operational capabilities of various protection schemes for high-impedance grounding faults, three methods are applied to a self-built experimental model. A phase high-impedance grounding fault is set at 50km away from M end. The fault occurs at 2.25s and lasts for 0.05s. The transition resistors are set to 100Ω, 200Ω, and 1000Ω, [25] respectively. The comparison results of the three methods are shown in Table 5.
Since [24]'s sampling frequency is 4kHz, it requires 25 sampling points with a data window of 6.25ms, while the proposed method uses 10 sampling points with a data window of 1ms. Therefore, [24] has a larger computational burden and slower speed. The current-based differential protection scheme acts at 4.6ms under a 100Ω transition resistor, while the proposed method acts at 3ms. Therefore, the proposed method has faster speed. By simulating and recording the extreme nodes where the three methods show significant errors under different transition resistors, it is found that the proposed method only experiences computational errors at 1000Ω, indicating stronger tolerance to transition resistors. In conclusion, the proposed method has faster speed and stronger tolerance to transition resistors compared to the other two methods.

5. Conclusions

Based on the analysis of the voltage fault characteristics of mixed AC-DC multi-feed-in receiving-end AC lines, this paper proposes a series-connected protection scheme based on voltage waveform similarity. This proposal has practical significance for the reliable operation of hybrid AC/DC grids. Through simulation verification, it is evident that the proposed protection scheme can operate correctly under different fault distances, fault types, transition resistances, and noise interferences. Compared to other protection schemes, proposed method operates more rapidly within 3ms, exhibiting improved reliability, speed, and adaptability, while remaining unaffected by data asynchrony, thus meeting protective requirements. The proposed method is based on the analysis of multi-fed AC-DC system fault characteristics, theoretically applicable to fault protection in all types of mixed AC-DC interconnected systems under various fault conditions. Future research may further investigate the applicability of the protection scheme.

Author Contributions

Conceptualization, S.G.; methodology, S.G., and X.L.; software, X.L; validation, H.Z., and Y.D.; formal analysis, X.L.; investigation, S.G.; resources, X.L.; data curation, X.L.; writing—original draft preparation, X.L; writing—review and editing, X.L.; visualization, H.Z.; supervision, G.S.; project administration, G.S.; funding acquisition, G.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (51777166), National Natural Science Foundation of China International Cooperation and Exchange Program (52061635105), State Grid Corporation of China headquarters science and technology Program(52094020006U).

Data Availability Statement

Not applicable.

References

  1. A, Alassi; S, Bañales; O, Ellabban; G, Adam and C, MacIver. HVDC transmission: technology review, market trends and future outlook. Renewable and Sustainable Energy Reviews 2019, 112, 530–554. [CrossRef]
  2. Haleem, NM; Rajapakse, AD; Gole, AM; Fernando, IT. Investigation of fault ride-through capability of hybrid VSC-LCC multi-terminal HVDC transmission systems. IEEE Transactions on Power Delivery 2019, 34, 241–250. [CrossRef]
  3. Xiao, H; Zhang, Y; Xue, Yand Yao, W. Correcting the Calculation Method of Commutation Failure Immunity Index for LCC-HVDC Inverters. IEEE Transactions on Power Delivery 2022, 37, 5469-5472. [CrossRef]
  4. Zheng, H; Yuan, X; Cai, J; Sun, P and Zhou, L. Large-Signal Stability Analysis of DC Side of VSC-HVDC System Based on Estimation of Domain of Attraction. IEEE Transactions on Power Systems 2022, 37, 3630-3641. [CrossRef]
  5. Guo, C; Liu, W; Zhao, J; et al. Impact of control system on small-signal stability of hybrid multi-infeed HVDC system. IET Generation, Transmission & Distribution 2018, 12, 4233–4239.
  6. S. V, Unde; P. Gawande and S. Dambhare. New Algorithm for Protection of Double Circuit Transmission Lines Using Modal Currents. IEEE Transactions on Power Delivery 2019, 34, 1967-1977. [CrossRef]
  7. Hao, X; Li, Y; He, X et al. A rapid assessment method of commutation failure immunity level for hybrid multi-infeed HVDC transmission systems. Proceedings of the CSEE 2017, 37, 4986-4998+5218. [CrossRef]
  8. Zhou, H; Yao, W; Ai, X; Li, D; Wen, J; Li, C. Comprehensive review of commutation failure in HVDC transmission systems. Electric Power Systems Research 2022, 205, 107768. [CrossRef]
  9. Jing, L; Wang, B; Dong, X. Review of consecutive commutation failure research for HVDC transmission system. Electric Power Automation Equipment 2019, 39, 116-123. [CrossRef]
  10. Hao, X; Li, Y; He, X; Yi, X and Duan, X. A Rapid Assessment Method of Commutation Failure Immunity Level for Hybrid Multi-infeed HVDC Transmission Systems. Proceedings of the Chinese Society of Electrical Engineering 2017, 37. 4986-4998. [CrossRef]
  11. Zheng, Z; Song, D; Du, K; Xiao, X; Ren, J; and Xie, Q. A Continuous Fault Ride-through Scheme for DFIGs Under Commutation Failures in LCC-HVDC Transmission Systems. Journal of Modern Power Systems and Clean Energy 2023, 11, 1126-135. [CrossRef]
  12. Li, X; Cai, Z; Huang, M et al. Transient Power Converse in an AC/DC Interconnected Power Grid. Automation of Electric Power Systems 2012, 36, 61-66+85.
  13. J. Ma, W. Ma, Y. Qiu, and J. S. Thorp. An adaptive distance protection scheme based on the voltage drop equation,” IEEE Transactions on Power Delivery 2015, 30, 1931–1940. [CrossRef]
  14. Huang, K; Jiao, B; Yu, Y; Liu, H; Hao, Q and Liu, Y. Performance Improvement for Pilot Protection of Double Circuit Lines on the Same Tower Based on Power Converse Accurate Identification, 2023 IEEE International Conference on Advanced Power System Automation and Protection (APAP) China, 2023, 260-264. [CrossRef]
  15. Huang, N; Fan, C; Jiang, S. Research of Protection for AC Transmission Lines Based on Characteristic of Grid-side Inverter. Journal of Engineering 2019, 16, 3240-3245. [CrossRef]
  16. Huang, J; Gao, H; Zhao, L and Feng, Y. Instantaneous Active Power Integral Differential Protection for Hybrid AC/DC Transmission Systems Based on Fault Variation Component, IEEE Transactions on Power Delivery 2020, 35, 2791-2799. [CrossRef]
  17. J. Ma, Y. Wu, and A. G. Phadke, A novel AC line distance protection scheme for AC/DC hybrid system based on fault likeness factor. IET Generation, Transmission & Distribution 2020, 15, 912–926. [CrossRef]
  18. Zheng, J; Wen, M; Chen, Y; Shao, X. A novel differential protection scheme for HVDC transmission lines. International Journal of Electrical Power & Energy Systems 2018, 94, 171-178. [CrossRef]
  19. Psaras, Vasileios; Dimitrios, Tzelepis; Dimitrios, Vozikis; Grain Philip Adam; and Graeme, Burt. Non-Unit Protection for HVDC Grids: An Analytical Approach for Wavelet Transform-Based Schemes. IEEE Transactions on Power Delivery 2021, 36.5, 2634-645. [CrossRef]
  20. Tian, C; Li, K. A method of ambiguity resolution based on matrix ordering. Journal of Navigation and Positioning 2021, 9, 65-70.
  21. Song, J; Cui, Y; Li, X; et al. Load Curve Clustering Method Based on Euclidean Dynamic Time Warping Distance and Entropy Weight. Automation of Electric Power System 2020, 44, 87-94.
  22. Cai, S; Lu, Z; Chen, B; Guo, Lin; Qing, Z; and Yao, L. Dynamic Gesture Recognition of A-Mode Ultrasonic Based on the DTW Algorithm. IEEE Sensors Journal 2022, 18, 17924-17931. [CrossRef]
  23. Ma, J; Kang, J; Zhou, Y and Cheng, P. A Novel Pilot Protection Scheme for AC Line Connected to LCC-HVDC Inverter Based on Composite-Mode Model. IEEE Transactions on Power Delivery 2022, 37, 1627-1639. [CrossRef]
  24. Zhang, H; Cong, W; Kong, H; Chen, M and Wei, Z. Longitudinal Protection Method based on Voltage Waveform Comparison for AC / DC Hybrid System. 2020 IEEE/IAS Industrial and Commercial Power System Asia (I&CPS Asia), 2020, 937-942. [CrossRef]
  25. Gao, W; He, W; Wai, R; Zeng, X; Guo, M. High-impedance arc fault modeling for distribution networks based on dynamic geometry dimension. Electric Power Systems Research 2024, 229, 110109. [CrossRef]
Table 1. Model parameters.
Table 1. Model parameters.
Terminal voltage EM/EN(/kV) 230
Short-circuit ratio 5
MN line length (/km) 100
Two-terminal out-of-zone line length (/km) 30
Total length at receiving end (/km) 160
Positive-sequence parameters of line R1(Ω/km) 0.01782
Positive-sequence parameters of line ωL1 (Ω/km) 0.3139
Positive-sequence parameters of line ωC1 (S/km) 3.626×10-6
Zero-sequence parameters of line R0(Ω/km) 0.2952
Zero-sequence parameters of line ωL0 (Ω/km) 1.039
Zero-sequence parameters of line ωC0 (S/km) 2.414×10-6
Table 2. Comprehensive distance similarity of fault waveforms in area DMA.
Table 2. Comprehensive distance similarity of fault waveforms in area DMA.
Fault Type Transition resistance (Ω) DMA DMB DMC DNA DNB DNC
Single-phase ground fault 0.01 4.4298 0.3175 0.8093 4.6376 0.2579 0.2792
50 3.1682 0.3248 0.7994 3.2402 0.2459 0.2690
100 2.7807 0.3268 0.7998 2.6619 0.2279 0.2612
Double-phase ground fault 0.01 0.1387 3.3960 2.5062 0.4000 7.5495 4.4422
50 0.4916 6.4778 3.7474 0.1258 7.6393 4.1362
100 0.5539 6.5028 4.8863 0.1434 7.7454 4.1295
Three-phase ground fault 0.01 5.2445 8.6719 3.0781 6.3565 5.2792 6.3929
50 4.2174 8.2139 3.0262 4.4252 4.1751 6.8504
100 3.6431 3.5486 2.8565 3.5486 7.6422 6.8159
Phase-to-phase short circuit fault 0.4916 7.4778 4.8882 0.1258 6.6393 4.1362
Table 3. Simulation results of anti-jamming ability.
Table 3. Simulation results of anti-jamming ability.
Fault location Fault type DMA DMB DMC DNA DNB DNC
f1 (A-G) 0.4191 0.3818 0.8793 0.2577 0.2016 0.3326
(BC-G) 0.3387 0.7960 0.4305 0.3298 0.6940 0.4132
(ABC-G) 0.5454 0.5267 0.3940 0.6781 0.8451 0.6219
f2 (A-G) 0.4237 0.3780 0.8514 8.2461 0.2011 0.3291
(BC-G) 0.3398 0.7862 0.4297 0.3127 4.6844 5.4057
(ABC-G) 0.5401 0.5227 0.3899 4.6774 5.8057 4.6211
f4 (A-G) 3.4668 0.3669 0.7471 0.1971 0.2014 0.3279
(BC-G) 0.3598 5.7681 4.4179 0.3177 0.6873 0.4098
(ABC-G) 4.5129 4.5115 4.4187 0.6228 0.8047 0.6047
f5 (A-G) 0.4795 0.3215 0.7258 0.0904 0.1995 0.2525
(BC-G) 0.3671 0.7543 0.3938 0.1097 0.6574 0.3561
(ABC-G) 0.4911 0.4981 0.3477 0.6129 0.7987 0.5917
Table 4. Comprehensive distance similarity of out-of-zone fault waveform.
Table 4. Comprehensive distance similarity of out-of-zone fault waveform.
Fault location Fault type DMA DMB DMC DNA DNB DNC
10% (A-G) 6.7142 0.2533 0.0574 5.7713 0.2314 0.2732
(BC-G) 0.4351 4.4981 4.3163 0.1473 6.8017 4.1781
(ABC-G) 7.0384 6.9961 6.6941 7.4620 8.0194 6.9814
50% (A-G) 4.5695 0.2654 0.7600 5.1755 0.2522 0.2696
(BC-G) 0.4101 3.4312 3.6343 0.1347 7.2423 3.8155
(ABC-G) 4.8054 4.5203 3.3509 5.4671 6.2985 6.5022
90% (A-G) 1.9995 0.2467 0.7250 7.2727 0.1836 0.2470
(BC-G) 0.2174 3.3058 3.1673 0.0571 3.1982 2.0070
(ABC-G) 4.3491 4.1847 2.7319 5.4894 6.7844 6.0876
Table 5. Execution result of different models.
Table 5. Execution result of different models.
Model [24] Longitudinal differential protection Ours
Data window 6.25ms 1ms 1ms
Protection operation time 8ms 4.6ms 3ms
Error occurrence conditions 100Ω 200Ω 1000Ω
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated