1 Introduction

In order to reduce the impact of electricity generation on alarmingly unfolding climate change, policies have been implemented to increase the use of renewable energy resources for bulk electricity generation. Some of the best resources such as wind and tidal energy are located in offshore sites. High voltage direct current (HVDC) transmission is the only feasible mode of importing energy from large-scale offshore generation sites through long submarine cables [1]. Moreover, the HVDC technology based on voltage source converters (VSCs) can facilitate asynchronous connection between fluctuating and weak energy sources and AC load centers [1]. VSC-HVDC schemes are implemented using the efficient multilevel modular converter (MMC) technology. Multi-terminal HVDC (MT-HVDC) grids are proposed to collect energy from geographically dispersed green energy resources and transport to load centers located beyond a single political zone.

There are few MT-HVDC grids which have been commissioned [2, 3], while several others are being planned. However, protection systems of those grids are not capable of preserving the availability of the healthy part of an MT-HVDC grid during DC side faults. Unlike in current source converters (CSCs) based on thyristors, the fault currents through the MMCs based on half bridge switching cells cannot be easily interrupted. As any means of controlling fault current via reduction of DC side voltage, fault tolerant convertors such as those proposed in [4, 5] are not suitable for large MT-HVDC grids. Such fault current control strategies cause to cease the operation of the whole grid via voltage collapsing. Therefore, DC circuit Breakers (DCCB) such as [6, 7] are proposed to clear DC side faults. However, fault current interruption in HVDC transmission circuits is challenging due to rapid rate of rise of currents and the absence of zero crossings [8]. Therefore, MT-HVDC grids require ultra-fast fault detection and isolation.

Fault discrimination is a key factor in MT-HVDC grid protection. A boundary inductor, an inductor placed between the converter and the transmission line end, is proposed to lower the rate of rising fault current [9]. Those boundary inductors ease the ability to discriminate faults in adjacent lines up to certain extend [9]. Commonly, the magnitudes of fault indicators observed or derived from terminal measurements are utilized to discriminate faults in the protected line from those on adjacent lines. Fault discrimination has been carried out using the magnitude of observed rate of rise of current [10], magnitude of the rate of change of line side voltage of the boundary inductor [9, 10], the time rate of voltage across the boundary inductor [11], or the magnitude of backward travelling wave [12]. However, large discrimination margins are required when the magnitudes of the fault disturbances are used for discriminating the faults, and as a result sensitivity of protection is reduced. Although communication based fault detection schemes such as current differential fault detection [13] and energy based fault detection methods [14, 15] are superior in terms of fault discrimination, any fault detection scheme which includes communication is not a likely to be satisfactory for primary protection of MT-HVDC grids as communication delays may not be tolerated.

Identification of the direction of a fault with respected the location of the protection relay (here after referred to as fault direction identification) is helpful to enhance the security. Directional relays such as [16] are widely used to improve fault discrimination in AC grids. The most common method of detecting the direction of a fault is comparison of the magnitudes of forward and backward waves which are calculated using both voltage and current measurements [17, 18]. Current derivatives have been used to discriminate forward and backward faults in [19]. However, both of above methods of detecting direction of fault requires terminal current measurements taken through DC current transformers (CT).

Since a fault created wave attenuates when travelling through a boundary inductor [12], the side of the fault with respect to boundary inductor can be easily and reliably identified if the magnitudes of disturbance at two sides are compared. A number of advantages are apparent in this method when compared with aforementioned methods of fault direction identification. As the rate of change of voltage (ROCOV) based fault detection methods [9, 11] applied for detection of transmission line and bus faults demands both line side voltage and bus side voltages, no additional sensors are necessary. Fault directions can be readily identified without using costly high bandwidth CTs that needed to measure fast current transients. Furthermore, the rate of dropping of voltages during DC side faults is much faster than the rate of rising of fault currents: this is because the rate of rising of current is approximately equals to the time integral of voltage variation. Therefore, a voltage based directional element can detect the direction of a fault promptly when a travelling wave reaches the terminal. Rapid discharge of filter and cable capacitance are dominant during the early phase rising fault current [20]. Since the maximum/steady state fault currents for short circuit faults are likely to be much higher than the circuit breaker capabilities [21], a fault must be detected and the DC circuit breakers (DCCBs) must be opened before the fault current rises above the breaker’s maximum interruption current. This aspect has not been tested in the fault detection schemes such as the one presented in [12]. The rate of change of voltage based HVDC grid protection scheme presented in [9] can face difficulties during high resistance faults.

In this paper, a method to identify the direction of a fault using local measurements and a basic arrangement to use the identified fault directions to improve the sensitivity, speed and reliability of the ROCOV based HVDC grid protection scheme described in [9] are presented. As no communication delay is involved, the local measurement based directional protection scheme can detect faults well before the fault currents reach the maximum breakable current. As shown in this paper, discrimination of high resistance faults at the far end of a long cable or transmission line is challenging [22]. To overcome such situations, a double ended ROCOV based directional comparison scheme is also presented.

2 Test grid and ROCOV based protection

The three-terminal VSC-HVDC test grid model developed in [9] is utilized in this study. The test grid shown in Fig. 1 is a bipolar-HVDC grid having a metallic return with a single earthing point. This configuration is one of the most likely practical HVDC grid configurations [9].

Fig. 1
figure 1

Structure of the MT-HVDC test grid

A bipolar low resistance single earth scheme requires a faster protection scheme than bipolar high resistance grounding and symmetrical monopole HVDC grids [23]. The test grid includes a long underground cable (Cable-14), a short underground cable (Cable-24), as well as a long overhead line (Line-24) to illustrate contrasting protection challenges. The grid was simulated in PSCAD/EMTDC electromagnetic transient (EMT) simulation software. Parameters of the test grid are shown in Table 1.

Table 1 Test grid details

The frequency dependent phase domain models available in PSCAD were used to simulate the cables and transmission lines. The simulation time step used is 10 μs. The MMCs connected to Buses 1, 2 and 3 were modeled with decoupled controls and converter transformers. It was assumed that basic fault protection is provided using the ROCOV based protection proposed in [9]. In this scheme, each line/cable section is protected using DC circuit breakers, and di/dt limiting inductors (also referred to as terminal inductors) placed at either end [24]. These di/dt limiting terminal inductors (values given in Table 1) relax the DC circuit breaker requirement and considered as vital for DC breaker application [24]. In simulations, the hybrid DC breaker model described in [6] was used with 2 ms operating delay and a maximum breaker current of 8 kA as in [9]. The PSCAD surge arrester model with default characteristic and 243 kV rating was used in the hybrid breaker circuit.

2.1 ROCOV based protection Principle

Sudden occurrence of high dv/dt in the voltage measured between the line side of a terminal inductor and ground is an indication of a transmission line fault [9, 19]. The fault detection and discrimination algorithm used in [9] is concisely depicted in Fig. 2. The process of selecting dv/dt settings is fully described in [9]. Since the fault generated voltage surge attenuate when passing through a di/dt limiting inductor, the peak dv/dt observed at the line side of terminal inductor is significantly smaller for bus faults and the faults in the adjacent lines, compared to the peak dv/dt observed for a fault in the protected line. This feature is used to discriminate the line faults from the faults on buses and adjacent branches [9].

Fig. 2
figure 2

ROCOV based protection algorithm proposed in [9]

A bus fault can be identified by comparing the dv/dt values on both sides of the all terminal inductors connected to the bus, as a bus fault generate a higher peak dv/dt on the bus side of all terminal inductors.

2.2 Signal processing in ROCOV based protection

The performance of a ROCOV based protection scheme can be significantly affected by the sensor bandwidth and the characteristics of the signal processing scheme. Assuming that practical high voltage sensors have at least 8 kHz bandwidth, any effects from the variations in the bandwidth of different sensors and measurement noise were eliminated by low pass filtering the measured voltage signals with a 8 kHz third order Butterworth filter.

Impact of the filtering on two voltage signals, one caused by a nearby fault and the other caused by a faraway fault, is shown in Fig. 3.

Fig. 3
figure 3

Measurements at line side of L 14 and L 41 with/without filtering during a solid fault (0.1 Ω) on Cable-14 at 400 km away from L 14

According to Fig. 3, when the fault is far from the measurement point, the difference between the peak dv/dt values calculated from the filtered (V* 41) and unfiltered (V 41) voltage waveform is negligible. However, for nearby faults, peak dv/dt calculated with filtered signal (V* 14) is significantly lower than that obtained from the unfiltered signal (V 14).

Although filtering removes the dependency of peak dv/dt on the location of the fault, this does not affect the protection scheme: only the attenuation of voltage transients by the di/dt limiting inductors is important for fault discrimination. This is because it is not required to discriminate between the close up faults and the far end faults on the same line. Initial investigations showed that 8 kHz bandwidth is sufficient to capture the maximum |dv/dt| for a far end fault on the shortest cable/overhead line in the test grid considered. The band limited voltage signal is then sampled at 32 kHz (sampling interval 31.25 μs > simulation time step 10 μs) and digitized with 12-bit ADC to make sure difference between peak dv/dt before and after the sampling remains less than 5%. ROCOV or dv/dt is digitally calculated between two time steps.

3 Challenges for ROCOV based DC grid protection

A dependable DC line protection scheme should be capable of detecting a fault at any point on the protected line regardless of the fault resistance. The faults must be detected and discriminated very quickly because if the fault current rises above the maximum breaker current, fault current cannot be interrupted subsequently. Furthermore, a secure protection scheme should not unnecessarily trip a circuit breaker. The setting of the ROCOV based DC line protection scheme in [9] is a compromise between these two aspects. When the threshold value of peak dv/dt is set higher, the operation is more secure, but fails to detect high resistance faults, especially when they are far from the relay. In the following sections, it is shown through simulations of the test MT-HVDC grid in Fig. 1 that a purely ROCOV based DC line protection strategy can fail for certain line faults.

All simulations were performed in PSCAD/EMTDC with a simulation time step of 10 μs and the fault was applied when t = 2.5 s in each simulation case. In the analysis, single pole to ground faults (negative pole to ground, for illustration purpose negative pole voltage is measured by keeping the reference probe on the ground and tip on the negative pole so that a fault shows a drop in voltage) were considered, because independent protection design for each pole is expected to be common in DC grid protection [25] and the pole to ground faults are considered as significantly more frequent [26]. The impact of the signal processing circuit used for capturing peak dv/dt was also taken in to account. All voltages are measured with respect to the ground.

3.1 High resistance faults at far end of a line

In order to examine the response of the typical ROCOV based protection scheme for faults with fairly high fault resistance, a fault having 6 Ω resistance is applied on Cable-14 at a point 450 km away from breaker B41. With the suggested ROCOV setting of 800 kV/ms for the breakers B14 and B41 as in [9], the expected response was obtained using the test grid model simulated in PSCAD/EMTDC. Fig. 4a shows the variations of voltages at line side of the di/dt limiting inductors at Bus-1 and Bus-4 (labeled as V 14 and V 41 respectively).

Fig. 4
figure 4

Measurements at a 6 Ω fault on Cable-14 at 450 km away from B41

The corresponding calculated dv/dt values are shown in Fig. 4b. The peaks in dv/dt values correspond with the arrival of travelling waves. According to simulation results, the relay at breaker B14 generates a trip signal because the observed peak dv/dt is greater than its setting. As shown in Fig. 5, the breaker B14 breaks the fault current before current through the relay reaches to its maximum breakable current 8 kA.

Fig. 5
figure 5

Calculated negative-pole conductor current through B14 and B41

According to Fig. 4b, the peak dv/dt seen at breaker B41 for the same fault is 790 kV/ms, which is less than its ROCOV setting 800 kV/ms. Therefore, no trip signal is generated for the breaker B41 to clear the fault (non-detection zone). However, the fault current eventually exceeds the maximum breakable current of breaker B41, about 126 ms after the fault as can be seen in Fig. 5.

This simulation case demonstrates that simple ROCOV based protection scheme can fail to operate for high resistance faults far from the relay point. Since the peak dv/dt values observed by a relay during a fault are dependent on the location as well as the resistance of the fault, it is possible to construct a surface plot of peak |dv/dt| as in Fig. 6.

Fig. 6
figure 6

Peak |dv/dt| surface of Cable-14

The plot in Fig. 6 was generated for Cable-14 in the test grid. The setting of the ROCOV relay can be shown as a plane on the same plot. The faults with distance -resistance combinations resulting in peak |dv/dt| values below the plane of the setting will not be detected.

3.2 Impact of external Events

In addition to the faults in the protected line, there are other external events such as faults on the adjacent lines and opening of circuit breakers protecting the adjacent lines cause high values of dv/dt at the measurement point of a ROCOV relay. To avoid inadvertent operation of the relay for such external events, ROCOV threshold need to be set sufficiently above the highest value of dv/dt observed due to external events. First, the impact of external faults on the operation B43 (on overhead line) is examined by disabling the operation of breakers B14 and B41 for a 0.1 Ω fault on Cable-14 at 50 km away from B41. The dv/dt values observed at B41, B43 and Bus-4 are shown in Fig. 7. The peak dv/dt observed at B43 is 632 kV/ms as indicated on the graph. Next, the impact of opening of breaker B41 on the operation of relay at B43 is examined by repeating the simulation without disabling the breakers B14 and B41. This time a high resistance fault of 5 Ω is applied on Cable-41, 400 km away from B41.The observed dv/dt values are shown in Fig. 8 with the respective peak values. A peak dv/dt value of 643 kV/ms appears at B43 when the breaker B41 opened to clear the fault. The peak dv/dt of transient caused by the fault (between 2.5035 s and 2.504 s) is much smaller than this value.

Fig. 7
figure 7

Peak dv/dt values of V 14, V 41, V 43, and V 4 for a 0.1 Ω fault at 50 km away from B41 on Cable-14

Fig. 8
figure 8

Peak dv/dt values of V 14, V 41, V 43, and V 4 for a 5 Ω fault on Cable-14, 400 km away from Bus-4 and subsequent opening of B41

The results presented in Fig. 7 and Fig. 8 show that the peak dv/dt created when opening an adjacent breaker has the same degree of impact as a solid close up fault in an adjacent line. In this particular case, if the threshold for relay at B43 is set only considering the peak dv/dt values due to faults, it may malfunction during the opening of breaker B41. Some measures need to be taken to improve the security of protection during external events [9, 10].

4 Augmenting ROCOV based DC line protection with directional characteristics

Two potential vulnerabilities of the simple ROCOV based DC grid protection proposed in [9] and possible challenges in selecting protection settings were illustrated: ① higher ROCOV settings can cause non-detection of high resistance faults far from the relay point as described in Section 3.1, ② a lower ROCOV settings can cause false tripping during the opening of breakers protecting adjacent lines as described in Section 3.2. In order to overcome the above mentioned shortcomings, it is proposed to augment the basic ROCOV based protection scheme with additional supervision by an under voltage element and by incorporating directional properties to the ROCOV relay through local and remote measurements. Supervision with under voltage elements helps avoiding false operation of ROCOV based detection due to noise in sensing and signal processing circuits [16].

If the ROCOV relay can be supervised by a directional signal that indicates faults in the forward direction, the ROCOV settings can be lowered to cover high resistant faults without reducing the security. Such a directional relay will not be affected by the high dv/dt created during the opening of circuit breakers in the reverse direction. Furthermore, the directional ability greatly simplifies the coordination of settings.

Fortunately, there is a simple method to accurately determine the fault direction with respect to a di/dt limiting terminal inductor (relay location) using only the local signals. For practical overhead lines and cables, a fault generated voltage wave propagating along the line is mostly reflected at the terminal inductor [16]. The transmitted component is much smaller. Therefore:

  1. 1)

    the peak |dv/dt| observed at the line side of terminal inductor is much greater than that observed at the bus side when the voltage wave is initiated from a fault on the protected line.

  2. 2)

    the peak |dv/dt| observed at the line side of terminal inductor is much smaller than that observed at the bus side for faults on the bus or lines behind the inductor location.

Therefore, the fault direction with respect to a relay location can be determined by comparing the magnitudes of peak |dv/dt| values observed at either side of the terminal inductor associated with the respective breaker.

For example, consider the simulation results shown in Fig. 7 for a fault on Cable-14. The corresponding peak |dv/dt| values observed at the measurement points near Bus-4 and Bus-1 are shown in Fig. 9.

Fig. 9
figure 9

Peak |dv/dt| values for a fault 50 km away from L 41 on Cable-14

The peak |dv/dt| observed at line side of L 41 is more than four time times the peak |dv/dt| observed on Bus-4, indicating that fault is in the direction of Cable-41. The observed peak |dv/dt| at Bus-4 is about twenty eight times greater than that observed at the line side of L 42. Thus the fault is not in the direction of Cable-42. Similarly the peak |dv/dt| measurements across L 43 confirm that fault is not in the direction of Line-43. At the other end of the faulted cable, the peak |dv/dt| observed on line side of L 14 is much higher than that observed on Bus-1, indicating that the fault is in the direction of Cable-41. The decision process of the new ROCOV based protection solution can be explained considering a generic bus (labeled as Bus-Y) in a HVDC grid, as shown in Fig. 10. The notations for various voltage signals, lines, and breakers are indicated on the diagram. Any fault in the direction of the protected line from a breaker is defined as a forward fault and a fault on the bus side from a breaker is defined as a reverse fault. For example, faults on both Line-XY and Line-YP are forward faults for relay associated with breaker BXY whereas a fault on Line-XY is in the reverse direction for the relay associated with breaker BYP.

Fig. 10
figure 10

Notations and voltage measurements utilized at relays

4.1 Local measurement based directional ROCOV relay

Figure 11 shows the scheme proposed to implement the proposed directional ROCOV relay for line/cable protection. The relay protecting Line-XY only uses local measurements available at Bus-Y.

Fig. 11
figure 11

Local measurement based directional ROCOV scheme

The peak |dv/dt| at line side of the terminal inductor is compared with threshold “High-setting”, and the comparator output is supervised by the proposed directional signal. The directional element calculates the ratio between the peak |dv/dt| values observed at the line and bus sides of the terminal inductor L YX . If this ratio is greater than one, theoretically the fault is in the direction of the protected line. However, to increase the security, the direction is declared as forward only if the ratio of peak |dv/dt| values is greater than 1.5. In addition to the direction, output is also supervised by an under voltage element to enhance the security [19, 23]. The ROCOV threshold “High-setting” in this scheme can be set lower than a non-directional ROCOV relay [9]. The lower limit of the “High-setting” is dictated by the peak |dv/dt| created at the relay point during a solid fault on the remote bus in the forward direction. For example, “High-setting” of relay RYX should be higher than the peak |dv/dt| recorded for a fault on Bus-X.

4.2 Communication based high resistance fault detector

As the minimum value of “High-setting” in the local measurements based directional ROCOV relay is constrained, it may not be able to detect faults with a significantly high resistance. This happens when the peak |dv/dt| observed at the relay (on line side of terminal inductor) during high resistance fault is less than the peak |dv/dt| values caused by a solid remote bus fault. In order to deal with such significantly high resistance faults, a communication assisted ROCOV based directional comparison scheme is proposed.

In this scheme, the relay at one end of the line determines whether the fault is in the forward direction with respect to it, and communicates that information to the other end if the peak |dv/dt| observed at the line side of the terminal inductor is above a minimum threshold. A trip signal is issued when both relays indicate that the fault is in the forward direction. There is no supervision by under voltage in this scheme, because for very high resistance faults there may not be a significant low voltage. Figure 12 shows the logic scheme for generating the direction signal using relay RYX as an example.

Fig. 12
figure 12

Communication based high resistance fault detector

Direction is determined similar to the previous case by comparing the ratio of peak |dv/dt| values across the terminal inductor (L YX ) with a threshold of 1.2. The threshold was set slightly lower because the ratio between line side and bus side peak |dv/dt| values tend to be lower for high resistance faults due to absence of high frequency components in the transients generated by these faults. Directional signal is enabled only if the peak |dv/dt| observed at the line side of L YX is larger than the threshold “Low-setting”. Although communication based protection schemes cannot be used to protect HVDC transmission from low resistance faults due to high rate of rising current [10], communication delay can be tolerated in the case of high resistance faults which cause a slower rate of rise of current. For example, in the simulation results shown in Fig. 5, the fault current takes about 126 ms to exceed the limit of breaker B41. This gives sufficient time to apply a communication based protection with no additional cost as communication between convertor stations is expected to be an integral part of future HVDC grids [25]. The two line fault detection schemes operating in parallel satisfy both reliability and speed requirements for a range of fault resistances.

4.3 ROCOV Comparison Based Bus Fault Detection

As mentioned before, the voltage disturbance created by a bus fault is highly attenuated when it is passing through the di/dt limiting inductors at line ends. Therefore, the rate of change of bus voltage is much greater than the rate of change of voltages measured at the line sides of the di/dt limiting inductors. Thus, for all lines emanating from a bus, if the dv/dt values observed at line side of the inductors are smaller than the dv/dt observed on the bus, a bus fault can be declared. This bus fault detection logic is shown in Fig. 13. To improve the security, a fault is declared only if the maximum dv/dt value observed at the bus is greater than Setting_B. This fault detection logic can be supervised by an under voltage element to further improve the security as in [19, 22]. Drop of bus voltage below VB_Setting confirms the occurrence of a bus fault. Although a number of voltage measurements are taken to detect a bus fault, all measurements are local and no additional delays. The breaker BYX opens when F YX_1, F YX_2, or F Bus-Y becomes TRUE.

Fig. 13
figure 13

Bus fault detection logic

5 Protection Settings and Coordination

The di/dt limiting inductors at the line ends performs a key role in the proposed protection scheme for discrimination of faulted zones. Thus the DC grid design should be done in conjunction with its protection design. The terminal inductors should be determined considering the breaker constraints (this paper considers 8 kA maximum breaking current and 2 ms maximum operating time [6]). An initial guess for terminal inductance is found by applying a solid fault at the line end of the inductor and at the critical distance defined in [9] (distance travelled by waves along the cable or overhead line within half the breaker operating time). The simulations must be used to verify that selected terminal inductances ensure that the fault currents do not exceed the breaker limit within the maximum breaker operating time. The terminal inductors values given in Table 1 were calculated using the above procedure.

The thresholds for the under voltage elements supervising the fault detection schemes are to 85% of the nominal DC voltage as in [19]. The observations used to calculate the “High-setting” values are shown in Table 2. The aim is to protect the entire transmission line using the local measurement based directional ROCOV relay protection scheme for short circuit faults. Thus the peak |dv/dt| for a solid fault at the far end of the protected line is recorded. Then the peak |dv/dt| values observed for a solid fault on the remote bus (in forward direction) and that observed when clearing the fault are recorded. Sometimes, the peak |dv/dt| observed when clearing the fault is higher than that observed when the fault occurs. The threshold “High-setting” is found by rounding up the result of (1).

$$High\_setting = A - 0.3\left[ {A\text{ - max}\left( {B\;\text{or}\;C} \right)} \right]$$
(1)

A, B and C in (1) are column labels in Table 2.

Table 2 High-settings for line protection

In a ±320 kV DC grid, any fault with resistance higher than 40 Ω will not result in steady state current above the breaker maximum current. Any faults with higher resistance can be protected using slower differential schemes or over current magnitude relays. The “Low-setting” should be selected to make sure any fault that could exceeds the breaker maximum current is detected and cleared in time. Leaving some margin, “Low-settings” are set to cover faults with up to 50 Ω resistance in this example. Thus the “Low-setting” for the cables is set at 75% of the observed peak |dv/dt| for a 50 Ω fault at far end of the protected line. Although cable faults are usually low resistance in nature, the high resistance fault protection is provided for all cables so that any eventuality covered. The “Low-setting” of the overhead lines can be set to cover higher fault resistance due to low attenuation of travelling waves. Thus “Low-setting” for overhead lines is fixed at 75% of the peak |dv/dt| observed for a 200 Ω fault at far end of the protected line. This is beneficial as high impedance faults are more likely in overhead lines. The “Low-setting” values calculated for each relay is shown in Table 3.

Table 3 Low settings line protection and bus settings

The threshold “Bus-setting” is set below the peak |dv/dt| observed during a 50 Ω on the protected bus, and above the highest peak |dv/dt| observed during the opening of breakers in adjacent lines or during the faults on remote converter buses. Exception is the bus at which metallic return wire is grounded. If the metallic return is grounded close to a bus, that bus can be more effectively protected using overcurrent or di/dt relays, as peak |dv/dt| measured at such a bus (example bus-3 in test system) would not change significantly due to strong dc sources (VSCs) connected between the ground and poles. For example, if “Bus-setting” for Bus-3 is set at 200 kV/ms, it can only detect faults with a resistance of 4 Ω or less. If the fault current rises above the breaker limit within the operating delay of the communication based scheme (direction detection time + communication delay + breaker operating time), the value of the terminal inductors of the line need to be increased or the safety margin of “High-setting” of the local measurement based schemes need to be decreased. In this paper, the communication delay (T D in ms) is approximately modeled as:

$$T_{D} = 5 + 0.5v_{o} l_{tx}$$
(2)

where v o is the velocity of light in free space (km/ms) and l tx is length of the line (km). The constant term in (2) represents the processing delays in equipment while the second term represents the signal travel time. Speed of wave propagation in communication media is assumed half of the speed of light, which is a worst case scenario [13]. It is required to test whether the settings are satisfactory under both directions of power flows (if allowed). This is because the time taken by fault currents to rise above the maximum breakable current depends on the direction of pre-fault load currents. If any inductor value is decreased, the all settings need to be re-evaluated.

6 Validation

The proposed ROCOV based protection scheme was implemented for the simple ± 320 kV HVDC grid system presented in Section 2, with the settings in Tables 2 and 3. Simulation model in PSCAD/EMTDC included sufficient details including the signal processing blocks. In order to demonstrate the operation and need for the communication based protection scheme, a 60 Ω fault was applied on the overhead line (Line-34), 1200 km away from B34. Some illustrative results are shown in Figs. 14 and 15.

Fig. 14
figure 14

Measurements during a 60 Ω fault on OHL-34 at 1200 km away from B34

Fig. 15
figure 15

Signals at R34 and R43

The local measurement based protection scheme at B43 declares a transmission line fault as the observed peak |dv/dt|, 3223 kV/ms is higher than the “High-setting” (1550 kV/ms), and the direction is found to be forward. Also, detection of forward fault is communicated to the remote relay at B34. Since the peak |dv/dt| observed at B34 (2295 kV/ms) is lower than the “High-setting” of R34 (2650 kV/ms), the local measurement based ROCOV unit at R34 did not detect this fault. However, directional element of the communication based scheme in R34 detects a forward fault, and the observed peak |dv/dt| is above its “Low-setting” (1450 kV/ms). Thus upon reception of the directional signal from R43 at the other end, occurrence of a fault in the protected line was confirmed.

The fault detection capability, satisfaction of breaker constraints, and mal-operation of proposed protection schemes are evaluated for different faults on the test HVDC grid. Table 4 summarizes the performance of the combined protection scheme for faults on overhead line. The fault resistances are given in the column labeled as RF. The columns labeled R 43-1 and R 43-2 indicate detection by the local measurement based and communication based protection schemes at R43 respectively. A tick (√) indicates correct detection while a cross (x) indicates a failure. The columns R 34-1 and R 34-3 corresponds to the protection schemes at R34. Columns I R43 and I R34 indicate whether the fault is cleared before the fault current exceed the breaker maximum current at the respective ends. The last column of Table 4 (labeled as S) shows the security of the overall protection and a tick indicates that protection schemes of Cable-14, Cable-24 or bus protection did not incorrectly responded to the fault. The performance of new protection scheme is especially investigated at high fault resistances to evaluate the improvement. Very low resistances faults are applied to check security issues.

Table 4 Protection performance during faults on overhead line

According to the results shown in Fig. 6 the peak |dv/dt| for a fault with a given fault resistance does not smoothly decrease with the distance from the measuring bus. This is prominent in the case of low resistance faults, and the highest peak |dv/dt| is not observed at zero distance, but at a critical distance from the bus. This is due to the voltage doubling effect due to reflected waves and the phenomenon is described in [9]. The same effect can be observed for cables in Tables 8 and 9. Also, when the fault is very close to the remote end, the remote breaker operates immediately by its local protection and the peak |dv/dt| due to breaker operation may become higher than the peak |dv/dt| due to fault, especially when the fault resistance is higher. This is evident from the last two rows of Table 4: higher peak |dv/dt| is observed by R34 for a 200 Ω fault compared to a 110 Ω fault.

By applying faults on Bus-1, security of the line protection scheme at R14 and fault detection performances of the bus-protection scheme is evaluated. The results are summarized in Table 5. The calculated peak |dv/dt| value at line side of the terminal inductor on B14 is much very smaller than bus side. Column F Bus-1 indicates the output of the bus protection scheme for Bus-1.

Table 5 Protection performance for Bus-1 during bus faults

Table 6 shows performances bus protection scheme in the other converter bus, Bus-2. Table 7 summarizes the fault detection performances of the bus-protection scheme applied to the non-converter bus, Bus-4. The calculated peak |dv/dt| values at line side of the terminal inductors on cables (at B42 and B41) are not included in Table 7 as they were much smaller (< 150 kV/ms).

Table 6 Protection performance for Bus-2 during bus faults
Table 7 Protection performance for Bus-4 during bus faults

Although the peak |dv/dt| values observed at the line side of L43 during bus faults is higher than the “High-setting” of R43, the directional properties of the new local measurement based fault detection scheme prevents issuing of trip signals. The performance of the protection schemes for Cable-24 is shown in Table 8.

Table 8 Protection performance during faults on Cable-24

The communication based protection is required at both ends of Cable-24. When compared with overhead lines, the peak |dv/dt| values observed in cables drop rapidly with the increasing fault resistance. Therefore, the maximum detectable fault resistance in cables is lower than the overhead lines (especially in long cables). However, this is not an issue as high resistance faults are not generally expected in cables. However, high attenuation of wave travelling through cables creates a challenge to protect long cables using ROCOV based protection. Table 9 shows the observed performances of the protection scheme for faults on long cable, Cable-14.

Table 9 Protection performance during faults on Cable-41

According to Tables 8 and 9, at a relay that protecting a cable, a small value of peak |dv/dt| is observed. Therefore, communication assisted protection scheme is very important to increase the sensitivity of the relay without compromising the security.

The proposed algorithms responds when the fault generated travelling waves arrive at its measurement points. Only delays involved are the measurement and filtering delays. The detection times are less than 200 μs even for the longest transmission line in the case of solid faults. High impedance faults result in longer detection times: for example, it takes about 200 μs to detect a 50 Ω fault on Cable 14, at 120 km from B14. If this fault was a solid fault, it is detectable within 50 μs.

7 Conclusion

The paper investigated ROCOV based protection strategies for a simple HVDC grid, and showed that simple ROCOV based protection can fail to detect high resistance faults at locations far from the relay. The side of the fault with respect to the location of the relay can be determined by comparing the peak |dv/dt| values observed at two sides of the terminal inductors. Directional element can be used to create a local measurement based fast and secure line protection relay. The communication assisted ROCOV based directional comparison scheme proposed in this paper can be designed to clear high resistance faults before fault currents exceeds the DC breaker limit. In this paper, it is demonstrated that entire spectrum of fault resistances can be covered by using the proposed local measurement based directional ROCOV relay in parallel with a communication assisted two-end ROCOV relay. Importance of communication assisted scheme for cable protection is demonstrated. For a given breaker design and typical MT-HVDC grid, the proposed protection is evaluated and found that protection is capable of clearing faults in any location of the transmission system before breaker current exceeds its rating. The proposed protection scheme can be implemented with practical hardware: a voltage sensor having a bandwidth of 8 kHz, 32 kHz sampler, and a 12 bit analog to digital converter.