Next Article in Journal
Innovative Membrane Electrode Assembly (MEA) Fabrication for Proton Exchange Membrane Water Electrolysis
Previous Article in Journal
Towards Overcoming the Curse of Dimensionality: The Third-Order Adjoint Method for Sensitivity Analysis of Response-Coupled Linear Forward/Adjoint Systems, with Applications to Uncertainty Quantification and Predictive Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Overvoltage and Insulation Coordination in Mixed HVDC Transmission Lines Exposed to Lightning Strikes

Department of Electronic Engineering, Hanyang University, Hanyangdaehak-ro 55, Ansan 15588, Korea
*
Author to whom correspondence should be addressed.
Energies 2019, 12(21), 4217; https://doi.org/10.3390/en12214217
Submission received: 2 October 2019 / Revised: 31 October 2019 / Accepted: 3 November 2019 / Published: 5 November 2019
(This article belongs to the Section F: Electrical Engineering)

Abstract

:
Many geographical constraints and aesthetic concerns necessitate the partial use of cable sections in the High Voltage DC (HVDC) transmission line, resulting in a mixed transmission line. The overhead sections of mixed lines are exposed to lightning strikes. The lightning strikes can not only result in flashover of overhead line (OHL) insulators but can enter the cable and permanently damage its insulation if adequate insulation coordination measures are not taken. In this work, we have analyzed the factors that affect the level of overvoltage inside the cable by simulating a fast front model in PSCAD. It has been determined that surge arresters must be provided at cable terminals when the length of cable sections is less than 16 km to limit the core-ground overvoltage within the lightning impulse protective level (LIPL). The level of sheath-ground overvoltage is independent of the length of cable; however, it can be limited within LIPL by lowering the sheath grounding impedance to 1.2 Ω. Insulation coordination measures do not impact the likelihood of OHL insulators’ flashover. The flashover performance of OHL can be improved by lowering the footing impedance of the second tower closest to the cable terminals, which is otherwise most likely to flashover.

1. Introduction

HVDC transmission found its first commercial application over five decades ago. The practicality of HVDC projects at that time was limited due to high cost and technical challenges. HVDC transmission projects were deemed suitable for either long distance terrestrial power transmission where they would offer economy over the HVAC transmission, or submarine energy transmission where HVAC transmission was excessively inefficient [1,2]. With the advancement in converter technologies HVDC systems are becoming more and more cost effective and technically viable. Therefore, the application of HVDC transmission has experienced a tremendous growth over the past two decades.
The use of cable section between overhead line sections is sometimes necessitated due to geographical constraints or concern of visually polluting scenic landscapes. Such transmission lines are called mixed HVDC transmission lines. Nearly 25% of the existing HVDC transmission projects fall into this category. The lightning can strike the overhead sections of the mixed line. The resulting overvoltage surge can not only cause the flashover of OHL insulators but can also travel towards the cable. If the overvoltage inside the cable exceeds the lightning impulse withstand level (LIWL) of insulation layers of the cable a permanent failure of cable is possible [3,4,5]. Therefore, the possibility of such events should be thoroughly analyzed, and overvoltage protection strategies should be considered.
Many researchers have studied the impact of lightning strikes on mixed transmission lines since the 1990s and several resulting articles have attempted to address various aspects of lightning transients in mixed HVDC as well as HVAC lines [6,7,8,9,10,11,12,13,14,15,16]. Many of these studies focused on specific design aspects such as the shunt reactor on cable terminal, presence of surge arrester/sheath voltage limiter, OHL geometries and tower footing impedance. Others focused on the modelling details such as corona, frequency dependent sheath grounding impedance as well as tower footing impedance. In our previous study [15], we presented a thorough theoretical analysis supported by fast front transient simulation regarding the transient behavior of mixed HVDC transmission lines. The insulators of OHL most vulnerable to flashover were identified and the reason for their high vulnerability was explained.
In this paper we have analyzed the maximum overvoltage stresses experienced by the cables’ insulation layer by simulating a high frequency model of mixed HVDC transmission line in time domain simulation software PSCAD (X4, v4.6.2). We have found that the core-ground overvoltage in a cable is dependent on the length of the cable. Due to the positive reflection coefficient at the receiving end boundary of the cable, the voltage surges reflected from the boundary have the same polarity as the incident surge and they superimpose resulting in amplification of voltage inside the cable. Sheath-ground overvoltage on the other hand are independent of cable length due to the negative reflection coefficient at their receiving end boundary. The surge arrester must be provided at the terminals of the small cables to provide an adequate protection margin to the cable’s main insulation layer, whereas in longer cables the use of a surge arrester is not necessary. The sheath-ground overvoltage can be limited within the lighting impulse protection level (LIPL) of the jacket by lowering the sheath grounding impedance.
It has been found that the insulation coordination measures do not alter the transient behavior of the mixed transmission line studied in our previous work [15]. The second tower closest to the cable boundary remains most vulnerable to flashover despite the installation of a surge arrester at cable terminals and lowering the sheath grounding impedance. However, the possibility of flashover of the positive pole insulator can be eliminated by lowering the footing impedance of the aforementioned tower.
In the following sections we have explained the system and details of fast front transient models, results and discussion and finally the conclusion.

2. System Description and Fast Front Transient Models

In this work we have considered a ±500 kV bipolar mixed HVDC transmission system shown in Figure 1. A cable with finite length up to 100 km is used for analysis. The sheath of the cable is grounded at the terminals. Five towers closest to both of the cable terminals are represented in detail. The overhead line is terminated with matching impedance at both ends to emulate an infinite transmission line.
Fast front models of various components used to represent the mixed HVDC transmission system are described below.

2.1. Tower

To emulate the overvoltage stress experienced by the insulators accurately, the vertical section of the tower is split into two sections connected at cross arms as shown in Figure 2b.
Due to the relatively small structure of HVDC towers, the multistorey tower models [17,18,19] yield inaccurate results. Therefore, we have represented the HVDC tower by a frequency independent distributed parameter line model composed of single conductor. The characteristic impedances of upper and lower tower sections are considered to be same [20,21].
Cross arms and braces retard the wave propagation on the tower. Therefore, the wave velocity of the vertical tower section is 85% of the speed of light [22]. The characteristic impedance of the waisted tower structure, shown in Figure 2, is calculated using Equation (1) [23].
Z T = 60   l n { cot ( 1 2 t a n 1 T ) }
where,
Tintermediate value →{r1h2 + r2(h1+ h2) + r3h1}/(h1+ h2)2
h1height from base to waist,
h2height from waist to top.
r1radius of tower at top,
r2radius of tower at waist,
r3radius of tower at base.

2.2. Footing Impedance of Tower

The level of overvoltage on the top of the tower and overvoltage stress experienced by OHL insulators is strongly dependent on the tower footing impedance. Therefore, prudently selecting a model of footing impedance is vital to obtaining reliable results.
The impedance of high resistivity soil has strong frequency dependence [24,25]. However, frequency dependent models have shown optimistic results regarding lightning overvoltage [14]. Moreover, the extent of frequency dependence of soil depends on resistivity which varies considerably according to the time of year. Therefore, to conservatively estimate overvoltage, we have represented the footing impedance by a fixed resistance rf [15,17,18].

2.3. Insulators

Suspension insulators shown in Figure 2 are comprised of twenty-nine disks of serially connected cap and pin type ceramic insulators (80 pF/disk). The air gap between arc horns is 5 m [26].
The flashover of an insulator is predicted using a leader propagation model (LPM). The time taken to breakdown is calculated by taking into account the physics of the discharge process and the gap configuration [27]. The total breakdown time, tc is given by Equation (2):
t c = t i + t s + t l
where,
tiCorona inception time, µs
tlLeader propagation time, µs
tsStreamer propagation time, µs
The corona inception voltage is very small compared to the breakdown voltage. The corona inception time ti is negligible, owing to the swift increase in voltage, and is assumed to be zero. If the voltage keeps increasing after the corona inception, a streamer begins to grow and cross the gap after time ts, in Equation (3):
1 t s = 1.25 ( E E 50 ) 0.95
where,
E50Average gradient at the critical flashover (CFO) voltage, kV/m
EMaximum gradient in the gap before breakdown, kV/m
After ts a leader starts to develop, with an exponentially increasing velocity given by (4),
d L d t = K V ( t ) [ V ( t ) g L E 0 ]
where,
LLeader length, m
gGap length, m
V(t)Voltage across gap, kV
KLeader coefficient, m2/kV2.s
E0Critical leader inception gradient, kV/m
K and E0 for cap and pin type suspension insulators under negative polarity stress are given in Table 1. Leader crosses the gap after time tl and the breakdown occurs.
The breakdown of the insulator is emulated by a switch. Equivalent capacitance cins, of the insulator as shown in Figure 2b is connected in parallel to the switch. The inductive characteristics of the arc path are represented by an inductance Larc connected in series with the switch [28].

2.4. Transmission Line

The cables and overhead lines are represented using the frequency dependent (phase) model commonly known as the universal line model (ULM). The ULM can adequately represent the transmission line for a wide range of frequencies [29,30].
The cable and overhead conductor are chosen to have a power rating of 660 MW/pole of the HVDC line. Midspan conductor sag for Overhead Power Conductors (OHPCs) as well as Overhead Ground Wires (OHGWs) is chosen as 10 m. OHGWs are modelled using ½ high strength steel, whereas Parrot ACSR 54/19 is used to model the OHPCs.
The 500 kV underground cables with a conductor cross section of 2100 mm2 are based on [31]. The two cables of the bipolar system are placed horizontally as shown in Figure 3. The main insulation layer of the cable has a lightning impulse withstand level (LIWL) of 2.7 p.u., whereas the jacket has a LIWL of 0.2 p.u. To ensure reliability and longevity, a protection margin of 20% should be provided [5].
Wave velocity and characteristic impedance of cable and overhead line conductors are calculated using the line constant program (LCP) of PSCAD and are listed in Table 1.

2.5. Lightning Strike

The Lightning strike is represented by a current source connected in parallel with the resistance. The parallel resistance emulates the impedance of the lightning channel [32].
A lightning current waveform with an upwardly concave front, as shown in Figure 4, resembles closely the records of lightning current wave-shapes given in [27,33]. The mathematical representation of this waveshape can be seen in Equation (5):
i ( t ) = { A t + B t n , t < t n I 1 e ( t t n ) t 1 I 2 e ( t t n ) t 2 , t > t n
constants of the wave shape listed in Table 2, have been derived according to [27].
The total probability of a lightning current exceeding 200 kA is 0.77% [33]. Therefore, the peak magnitude of the lightning strike applied to OHGWs or the tower top is chosen to be 200 kA. The IEEE modelling guideline also recommends the same for insulation coordination studies [26].
I M S F = [ γ ( h m + h p ) / 2 A ( 1 γ   s i n α ) ] 1 / B
The, maximum shielding failure current (IMSF) is calculated by using an electro-geometric model (EGM) proposed by IEEE Std [34]. By solving Equation (6) based on the geometric model given in Figure 5, the IMSF is calculated to be 20 kA. Where, hm: Height of OHGWs, hp: Height of OHPCs, α, β and γ are as recommended by IEEE Standard are as follows. α: Shielding angle, β = 0.36 + 0.17 ln (43 − hp) ⇒ hp ≤ 40 m, γ: 1/β: the parameters of the lightning impulse wave shape are listed in Table 2.

2.6. Surge Arrester

The internal overvoltages in the line commutated converter (LCC) HVDC system are not expected to exceed 1.5 p.u. for the maximum duration of 100 ms. The insulation of the HVDC cable is capable of withstanding this stress and therefore protection does not need to be provided against internal overvoltages [5].
The lightning strikes on the overhead sections of the line would expose the cable to overvoltage stress, which can lead to permanent damage of the cable if the said overvoltage exceeds the LIWL of the insulation layer of cable. The level of overvoltage in the cable resulting from lightning strikes depends on several design parameters of the mixed HVDC transmission line.
Since there is no need to provide protection against internal overvoltage, a surge arrester with a maximum continuous operating voltage (MCOV) of 750 kV, i.e., 1.5 p.u. of the system voltage, and the following protective characteristics will be employed if the voltage inside the unprotected cable exceeds 2.25 p.u. [6].
  • Lightning impulse protective levels (LIPL), V10: 1130 kV (10 kA, 8 × 20 µs)
  • Switching impulse protective levels (SIPL), Vss: 960 kV (3 kA, 30 × 60 µs)
The lightning impulse protection margin for this surge arrester can be calculated using Equation (7).
P M = ( L I W L L I P L 1 ) 100
This provides a desired protection margin (PM) of 20% which is in accordance with the recommended minimum level of 20% for HVDC cables [5].
Metal oxide surge arresters (MOSA) are used universally for providing protection against surge voltages in a range of power system applications. However, the residual voltage across the arrester is a function of both the rate of rise and the magnitude of the current conducted by the arrester. The waveshape of the voltage impressed on the surge arrester can vary considerably depending on the location of incidence of the lightning strike and on the characteristics of the lightning strike and of mixed transmission line.
In order to emulate the characteristics of a real MOSA, a frequency dependent model proposed by IEEE WG 3.4.11 shown in Figure 6 will be used in this paper [35].
The nonlinear V-I characteristic of an arrester is represented using two nonlinear resistors A0 and A1, shown in the Figure 7. The nonlinear resistors are separated by an RL filter, which offers very low impedance to slow front surges, resulting in parallel conduction of the two resistors. However, for fast front surges the impedance of the RL section become quite significant. Therefore, more current is conducted by the nonlinear resistance ‘A0’. This results in a higher voltage across arrester terminals because the nonlinear resistor ‘A0’ has a higher voltage for a given current than the resistor ‘A1’, as can be seen in Figure 7. These V-I characteristics are obtained by modifying the characteristics proposed in [35] to match the V10 and VSS required in this study.
The inductance ‘L0’ in the model represents the magnetic field in the vicinity of the surge arrester, whereas the capacitor ‘C’ represents the terminal-to-terminal capacitance of the arrester.
The parameters of the surge arrester model can be calculated using Equations (8)–(12).
L 1 = 15 d n   µ H
R 1 = 65 d n  
L 0 = 0.2 d n   µ H
R 0 = 100 d n  
C = 100 n d   pF
where,
dheight of the arrester [m]
nnumber of parallel columns of metal oxide in arrester
The parameters of the surge arrester model are listed in Table 3.

3. Results and Discussion

Lightning can strike the overhead power conductors (OHPC) or the overhead ground wire (OHGW) and tower top as discussed previously. In this section we will present and analyze in detail the overvoltage stresses experienced by the insulation layers of the cable resulting from lightning strikes. We will further analyze whether these overvoltages exceed the lightning impulse withstand level (LIWL) of the cable and whether insulation coordination measures are necessary. Various insulation coordination measures will be analyzed. Finally, the impact of insulation coordination on transient behavior of a mixed HVDC transmission line will be presented and analyzed in depth.

3.1. Overvoltage Analysis

The lightning strikes closest to the cable, i.e., on the transition tower will cause the maximum overvoltage in the cable. The strike on the OHPC of the transition tower will enter the cable and cause core-ground (CG) as well as sheath-ground (SG) overvoltage, whereas a strike on the tower top or OHGW will result in sheath-ground (SG) overvoltage.

3.1.1. Lightning Strike on an Overhead Power Conductor

Core-Ground Overvoltage

Negative lightning strikes are most likely to occur [27]. A negative lightning strike on a negative pole OHPC at twr1S will result in the constructive superposition of voltage, hence subjecting the cable to the highest overvoltages. The equivalent circuit is shown in Figure 8. Only two towers are represented in detail in this equivalent circuit. The OHPC can be seen to be isolated from the tower by an open switch. Therefore, a negative lightning strike on a negative OHPC can initiate a voltage surge ef1 travelling on the OHPC in the forward direction towards the cable and another voltage surge er0 traveling towards the DC source in the reverse direction as shown on the equivalent circuit in Figure 8. A portion of the wave ef2 is refracted inside the cable, whereas the rest, er1, is reflected backward at boundary a. The wave ef2 travelling forward inside the cable will be again reflected and refracted at the receiving end boundary ‘b’ of the cable. The magnitude of reflected and refracted waves can be determined by a traveling wave theory as per Equations (13)–(15).
e f 2 = e f 1 2 Z c o r e Z O H P C + Z c o r e
e r 1 = e f 1 Z c o r e Z O H P C Z O H P C + Z c o r e
e r 2 = e f 2 Z O H P C Z c o r e Z O H P C + Z c o r e
The simulation results for the resulting overvoltage at the boundary ‘a’, i.e., the entrance of a cable according to cable length are shown in Figure 9.
It can be seen that, as the length of cable decreases, the magnitude of overvoltage increases. The characteristic impedance of the OHPC is much higher than that of the cable core conductor, making the reflection coefficient positive. Therefore, as the voltage surge travels from the cable to the OHPC, the same polarity surge is reflected from their boundary ‘b’. In a shorter cable section, the reflected wave arrives at the sending end very quickly and without sufficient attenuation. This same polarity reflected voltage surge superimposes with the incident surge resulting in an increase in the magnitude of the overvoltage. It can be seen that in the cable sections shorter than 12.5 km, multiple reflections superimpose on the incident surge even before the tail of the incident surge reaches its 50% value.
On the other hand, in a 100 km long cable section, the reflected wave arrives just before 8.1 ms. The incident wave reduces to less than 10% before the arrival of the reflected wave. Moreover, due to the longer distance, the reflected wave is highly attenuated. Therefore, it can be seen that in longer cable sections the maximum overvoltage depends only on the magnitude of the lightning strike and characteristics impedances of the OHL and the cable.
However, the maximum overvoltage does not always occur at the entrance of the cable. Therefore, to determine the maximum overvoltage, the cables with the lengths shown in Figure 9 are divided into ten segments. The voltage at the boundaries of each section is calculated. The results are presented in Figure 10.
To easily visualize the results of Figure 10, the plot of the maximum overvoltage along the length of the cable is shown in Figure 11. The location and magnitude of the maximum overvoltage in each case is listed in Table 4.
It can be seen that in longer cable sections, i.e., 100 km, 50 km and 25 km, the maximum overvoltage occurs at the cable entrance. The magnitude of the maximum overvoltage decreases towards the center of the cable before increasing again near the receiving end of the cable. A close comparison of 100 km and 90 km plots in Figure 10a shows that in the case of the former, the reflected wave instantly superimposes on the incident wave upon its arrival, whereas in the case of the latter there is a displacement between the incident and the reflected wave equal to 20% of the travel time of the cable. This leads to the decrease in amplitude of the incident wave before the arrival of the reflected wave. Therefore, the magnitude of the maximum voltage continues to decrease as we move away from the receiving end.
In the shorter cable sections, i.e., 12.5 km, 6.25 km and 3.125 km, the maximum overvoltage occurs at 10% from the receiving end. This can be attributed to the shorter travel times of the reflected wave, lower attenuation of the reflected wave and insufficient damping of the incident wave.

Sheath-Ground Overvoltage

The impulse current flowing through the core of the cable induces a voltage in the sheath of the cable. The voltage in the sheath of the cable stresses the jacket of the cable. If the voltage exceeds the LIWL of the jacket (0.2 p.u.), it may rupture. Therefore, in this subsection we have presented a detailed analysis of sheath-ground (SG) overvoltage considering the various factors that influence it.
The voltage at the sending end terminal of the sheath resulting from a direct lightning strike on the sending end OHPC is shown in Figure 12. It can be seen that the maximum overvoltage is the same regardless of the length of the cable. This behavior is in contradiction with the length dependent behavior of the CG overvoltage.
The SG overvoltage along the length of the cable at equally spaced intervals is shown in Figure 13. Compared to the CG overvoltage results shown in Figure 10, the incident wave takes much longer to arrive at the receiving end. This is due to the lower wave velocity in sheath. Moreover, due to the higher attenuation constant of the sheath, the surge voltage gets significantly attenuated away from the cable entrance.
In the smaller cable sections (50 km, 25 km, 12.5 km, 6.25 km and 3.25 km) the opposite polarity reflected surge arrives at the sending end within the time frame shown in Figure 13.
The sheath is grounded at both ends. Since the sheath grounding impedance is smaller than the characteristics impedance of the sheath, the reflection coefficient is negative, resulting in the positive polarity reflected waves. Therefore, the portion of the voltage surge reflected from the receiving end terminal cannot cause the amplification of the voltage even in very small cable sections. The maximum sheath voltage will always appear at the cable entrance regardless of its length.
The sheath-ground voltage along the length of the cable in cables with different total lengths is shown in Figure 14. Contrary to the core-ground voltages seen in Figure 11, the SG voltage is maximum at the cable entrance regardless of the cable length and reduces significantly at the receiving ends of the cable.

3.1.2. Lightning Strikes on an Overhead Ground Wire/Tower Top

A lightning strike on an OHGW/tower top results in voltage surge travelling along the OHGW in the reverse and forward directions, along with part of the wave travelling along the tower in the downward direction, as shown in the equivalent circuit of Figure 15.
When lightning strikes the OHGW/tower top far away from the cable, the major portion of the surge travels through the towers and is absorbed by the ground. Therefore, in this section we will only consider the OHGW strikes at the transition tower. A part of the voltage surge travels towards the cable through the riser section resulting in sheath-ground overvoltage. The simulation results of sheath-ground overvoltage at the entrance of the cable according to cable lengths are shown in Figure 16.
It can be seen that SG overvoltage at the sending end terminal is the same in all cases regardless of the cable length. Moreover, it can be noticed that the magnitude of the voltage in this case is about four times higher than that resulting from the OHPC strike. The SG voltage at equally spaced intervals along the length of cable are shown in Figure 17.
It can be noted again that the SG voltage decreases exponentially along the length of cable due to a high attenuation constant. Moreover, even in the smallest cable sections, the reflected waves cannot cause amplification of the voltage inside the cable due to the negative reflection coefficient at boundary ‘b’ as shown in Figure 15. The comparison of voltage magnitude along the length of the cable in cables with different lengths is shown in Figure 18.

3.2. Insuation Coordination Requirement

The stress on the insulation layers of the cable must be limited to acceptable levels to prevent the insulation failure under the influence of lightning strikes. If the CG overvoltage exceeds the desired LIPL the surge arrester should be installed at the terminals of the cable between the core and the ground, whereas if the SG voltage exceeds the desired LIPL of the jacket, the sheath grounding impedance must be lowered. The schematic of the surge arrester installed at the cable terminals is shown in Figure 19.

3.2.1. Core-Ground Overvoltage

The LIWL of the HVDC cable in this work is 2.7 p.u. Surge arresters as shown in Figure 19 will be applied at cable terminals if the CG voltage inside the cable exceeds the LIPL, i.e., 2.25 p.u., ensuring a protection margin of 20%.
The direct lightning strike has been applied at the negative OHPC of the sending end transition tower. Cable lengths in the range of 2.5 km to 100 km (with an increment of 2.5 km) are considered. The transient voltage in each case is measured at ten equally spaced intervals along the length of the cable. The maximum overvoltage in each case is recorded and plotted according to the length of the cable as shown in Figure 20.
It can be seen that the overvoltage exceeds 2.25 p.u. only if the cable length is shorter than 16 km. For cable lengths greater than 22.5 km, the cable overvoltage does not vary by changing the length of the cable. This implies that the reflection of the voltage surge from the receiving end of the cable takes sufficiently long to arrive at the sending end that amplification of the voltage is not possible. In such cases the maximum overvoltage always occurs at the sending end of the cable as can also be seen in Figure 11. Therefore, for this particular case of mixed HVDC transmission lines, the surge arrester should always be installed at the terminal of cables shorter than 16 km.

3.2.2. Sheath-Ground Overvoltage

Overhead Power Conductor Strike

It can be seen that in the worst case OHPC strike, i.e., at the sending end transition tower, SG overvoltage is about 0.22 p.u., significantly higher than the desired LIPL of 0.1667 p.u. However, the magnitude of the overvoltage is independent of the length of the cable.
The SG overvoltage at the terminal of the sheath according to the sheath grounding impedance is shown in Figure 21. It can be seen that the maximum overvoltage can be limited to the desired protective level if the sheath grounding impedance is limited up to 6.5 Ω.

Overhead Ground Wire/Tower Top Strike

In the worst case OHGW strike, i.e., at the sending end transition tower, overvoltage is more than five times the desired protective level. The maximum overvoltage in this case is also independent of the cable length and occurs at the sending end of the cable. The SG overvoltage at the terminal of the cable according to resistance of local ground is shown in Figure 22. It can be seen that the maximum overvoltage can be limited to the desired protective level if sheath grounding resistance is limited up to 1.2 Ω.

3.3. Insulation Coordination: Impact and Analysis

3.3.1. Core-Ground Overvoltage

The surge arresters are installed at both terminals of the cable as shown in Figure 19. The direct lightning strike is applied at the OHPC of the sending end transition tower. The maximum overvoltage along the length of cable is shown in Figure 23.
It can be seen that the overvoltage exceeds the desired protective level of the surge arrester only in the case of the smallest cable. However, it is still within the LIWL of the insulation. Therefore, damage to the cable is unlikely.

3.3.2. Sheath-Ground Overvoltage

The SG overvoltage along the length of the cable resulting from a direct lightning strike on the negative OHPC at the transition tower is shown in Figure 24. Despite the presence of surge arresters, the SG overvoltage exceeds the desired protective level of 0.16 p.u. in all cases.
The energy dissipation in surge arresters installed at both ends of the line is shown in Figure 25. It can be seen that energy dissipation in the receiving end arrester is larger than the sending end regardless of the length of the cable. Moreover, energy dissipation becomes higher with the decrease in the length of the cable.

3.3.3. Flashover Performance of Overhead Line Insulators

In our previous research [15] we determined that the negative pole insulator of twr2S is most vulnerable to flashover in cases of direct lightning strike on the negative OHPC on twr2S, whereas the positive pole insulator is most vulnerable in cases of indirect lightning strike on the OHGW or tower top. In this section, we have analyzed the impact of insulation coordination measures on the flashover performance of twr2S insulators.

Strike on Overhead Power Conductors

A lightning strike is applied to a negative OHPC at the subsequent tower, i.e., twr2S. The flashover performance with or without a surge arrester at the cable terminals is shown in Figure 26.
In our previous study [15] negative pole insulators of twr2S and twr3S were found to be most vulnerable. It can be seen in Figure 26 that the presence of a surge arrester at the cable terminal has virtually no impact on the insulator voltage. The insulator flashover occurs at both towers regardless of the presence of a surge arrester.

Strike on Overhead Ground Wire/Tower Top

In our previous study [15] we determined that in cases of lightning strike on the top of various towers, the plus pole insulators were most vulnerable to flashover. In the case of the lightning strike on the top of twr2S, the plus pole insulator of twr2S always experienced a flashover. The sheath grounding impedance ‘rsg’ was assumed to be constant, i.e., 10 ohms in that study. We have however, concluded that a low value of ‘rsg’ is necessary in order to achieve the desired protective level for the sheath of HVDC cables.
The flashover performance of the twr2S insulator considering varying sheath grounding impedances is shown in Figure 27.
It can be seen that lowering the sheath grounding impedance has virtually no effect on the flashover performance of a plus pole insulator. Due to the much higher characteristic impedance of OHGW compared to the sheath grounding impedance, any change in its value has a marginal effect on the reflection coefficient at the boundary of the cable. The impact of footing impedance of twr2S on flashover performance of its plus pole insulator is presented in Figure 28.
It can be seen that by lowering the tower footing impedance to 6 ohms the flashover of the insulator can be prevented. Lowering the footing impedance increases the reflection coefficient at the boundary c in Figure 8, resulting in an increased damping of voltage on the top of the tower.

4. Conclusions

Lightning strikes on the overhead section of the line subject the components of mixed HVDC transmission line to tremendous overvoltage stresses. These overvoltage stresses can cause flashover of OHL insulators and can cause cable insulation to fail permanently if they exceed the LIWL of the insulation layers.
The following outcomes have been drawn from this paper.
  • The maximum core-ground overvoltage increases with a decrease in the length of the cable. This behavior is experienced due to a positive reflection coefficient at the receiving end boundary of the cable. For cables, shorter than 16 km, the magnitude of overvoltage is expected to exceed the LIPL of the cable. Therefore, surge arresters should be installed at the terminals of such cables.
  • The maximum sheath-ground overvoltage is independent of the length of the cable due to the negative reflection coefficient at the receiving end boundary of the sheath.
  • The sheath-ground overvoltage can exceed the LIWL of the jacket significantly for high sheath grounding impedances. The SG overvoltage can be limited within the LIPL of the jacket by limiting the sheath grounding impedance up to 1.2 ohms.
  • As concluded in our previous work [15] the sending end subsequent tower “twr2S” is most vulnerable to flashover in the case if OHPC as well as OHGW strikes. The flashover due to OHGW strikes can be prevented by limiting the footing impedance of twr2S to 6 ohms.
The presence of surge arresters on cable terminals does not impact the flashover performance of OHL insulators.

Author Contributions

Conceptualization, M.A. and B.-W.L.; methodology, M.A.; software, M.A. and A.S.; validation, M.A.; resources, K.-H.P. and H.-Y.L.; writing—original draft preparation, M.A.; writing—review and editing, B.-W.L.; project administration, H.-Y.L.; funding acquisition, B.-W.L.; supervision, B.-W.L.

Funding

This work was supported by the Korea Electric Power Corporation, Republic of Korea (Grant number. R17XA05-3). The APC was funded by Korea Electric Power Corporation.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Jovcic, D.; Ahmed, K. Introduction to line commutated HVDC. In High Voltage Direct-Current Transmission Converters, 2nd ed.; Wiley: West Sussex, UK, 2019; Volume 1, pp. 1–3. [Google Scholar]
  2. Siepmann, W.; Kirby, N.M.; Xu, L.; Luckett, M. HVDC transmission for large offshore wind farms. Power Eng. J. 2002, 16, 135–141. [Google Scholar]
  3. CIGRÉ Brochure 496. Recommendations for Testing DC Extruded Cable Systems for Power Transmission at a Rated Voltage up to 500kV; Working Group B1.32; CIGRÉ: Paris, France, 2012. [Google Scholar]
  4. IEC 62067. Power Cables with Extruded Insulation and Their Accessories for Rated Voltages above 150 kV (170 kV) up to 500 kV (550 kV)—Test Methods and Requirements; Report; IEC: Geneva, Switzerland, 2011. [Google Scholar]
  5. CIGRÉ Brochure 86. Overvoltages on HVDC Cables Final Report; Joint Working Group 33/21/14.16; CIGRÉ: Paris, France, 1994. [Google Scholar]
  6. Bui-Van, Q.; Beaulieu, G.; Huynh, H.; Rosenqvist, R. Overvoltage studies for the St. Lawrence river 500-kV DC cable crossing. IEEE Trans. Power Deliv. 1991, 6, 1205–1215. [Google Scholar] [CrossRef]
  7. Goertz, M.; Wenig, S.; Gorges, S.; Kahl, M.; Beckler, S.; Christian, J.; Suriyah, M.; Leibfried, T. Lightning overvoltages in a HVDC Transmission System comprising Mixed Overhead-Cable Lines. In Proceedings of the International Conference on Power System Transients, Seoul, Korea, 26–29 June 2017. [Google Scholar]
  8. Goertz, M.; Wenig, S.; Hirsching, C.; Kahl, M.; Suriyah, M.; Leibfried, T. Analysis of Extruded HVDC Cable Systems Exposed to Lightning Strokes. IEEE Trans. Power Deliv. 2018, 33, 3009–3018. [Google Scholar] [CrossRef]
  9. Silva, F.F.; Pedersen, K.S.; Bak, C.L. Lightning in hybrid cable overhead lines and consequent transient overvoltages. In Proceedings of the International Conference on Power System Transients, Seoul, Korea, 26–29 June 2017. [Google Scholar]
  10. Henriksen, T. Calculation of lightning overvoltages using EMTP. In Proceedings of the International Conference on Power System Transients, Lisbon, Portugal, 3–7 September 1995. [Google Scholar]
  11. Colla, L.; Gatta, F.M.; Geri, A.; Lauria, S. Lightning overvoltages in hv-ehv “mixed” overhead-cable lines. In Proceedings of the International Conference on Power System Transients, Lyon, France, 4–7 June 2007. [Google Scholar]
  12. Massaro, F.; Morana, G.; Musca, R. Transient behaviour of a “mixed” overhead-cable EHV line under lightning events. In Proceedings of the International Universities Power Engineering Conference, Glasgow, UK, 1–4 September 2009. [Google Scholar]
  13. Hoogendorp, G.; Popov, M.; van der Sluis, L. Lightning Induced Overvoltages in Mixed 380 kV OHL-Cable-OHL connections. In Proceedings of the International Conference on Power System Transients, Vancouver, BC, Canada, 18–20 July 2013. [Google Scholar]
  14. Gomes, T.V.; Schroeder, M.A.O.; Alipio, R.; de Lima, A.C.S.; Piantini, A. Investigation of Overvoltages in HV Underground Sections Caused by Direct Strokes Considering the Frequency-Dependent Characteristics of Grounding. IEEE Trans. Electromagn. Compat. 2018, 60, 2002–2010. [Google Scholar] [CrossRef]
  15. Asif, M.; Lee, H.-Y.; Khan, U.A.; Park, K.-H.; Lee, B.W. Analysis of Transient Behavior of Mixed High Voltage DC Transmission Line Under Lightning Strikes. IEEE Access 2019, 7, 7194–7205. [Google Scholar] [CrossRef]
  16. Lennerhag, O.; Lundquist, J.; Engelbrecht, C.; Karmokar, T.; Bollen, M.H.J. An Improved Statistical Method for Calculating Lightning Overvoltages in HVDC Overhead Line/Cable Systems. Energies 2019, 12, 3121. [Google Scholar] [CrossRef]
  17. Ametani, A.; Kawamura, T. A Method of a Lightning Surge Analysis Recommended in Japan Using EMTP. IEEE Trans. Power Deliv. 2005, 20, 867–875. [Google Scholar] [CrossRef]
  18. Ametani, A.; Nagaoka, N.; Ohno, T.; Baba, Y. Transient on Overhead Lines. In Power System Transients: Theory and Applications, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2017; Volume 1, pp. 176–193. [Google Scholar]
  19. Ishii, M.; Kawamura, T.; Kouno, T.; Ohsaki, E.; Shiokawa, K.; Murotani, K.; Higuchi, T. Multistory transmission tower model for lightning surge analysis. IEEE Trans. Power Deliv. 1991, 6, 1327–1335. [Google Scholar] [CrossRef]
  20. Datsiso, Z.G.; Mikropoulous, P.N. Effect of tower modelling on the minimum backflashover current of overhead transmission lines. In Proceedings of the 19th International Symposium on High Voltage Engineering, Pilsem, Czech Republic, 23–28 August 2015. [Google Scholar]
  21. Martinez, J.A.; Castro-Aranda, F. Tower modeling for lightning analysis of overhead transmission lines. In Proceedings of the IEEE Power Engineering Society General Meeting, 22–25 May 2005; IEEE: San Francisco, CA, USA, 2005; pp. 1345–1350. [Google Scholar]
  22. A Simplified Method for Estimating Lightning Performance of Transmission Lines. IEEE Trans. Power Appar. Syst. 1985, PAS-104, 918–932. [CrossRef]
  23. Chisholm, W.; Chow, Y.; Srivastava, K. Travel Time of Transmission Towers. IEEE Trans. Power Appar. Syst. 1985, PAS-104, 2922–2928. [Google Scholar] [CrossRef]
  24. Alipio, R.; Visacro, S. Frequency Dependence of Soil Parameters: Effect on the Lightning Response of Grounding Electrodes. IEEE Trans. Electromagn. Compat. 2013, 55, 132–139. [Google Scholar] [CrossRef]
  25. Cavka, D.; Mora, N.; Rachidi, F. A Comparison of Frequency-Dependent Soil Models: Application to the Analysis of Grounding Systems. IEEE Trans. Electromagn. Compat. 2014, 56, 177–187. [Google Scholar] [CrossRef]
  26. Imece, A.F.; Durbak, D.W.; Elahi, H.; Kolluri, S.; Lux, A.; Mader, D.; McDemott, T.E.; Morched, A.; Mousa, A.M.; Natarajan, R.; et al. Modeling guidelines for fast front transients. IEEE Trans. Power Deliv. 1996, 11, 493–506. [Google Scholar]
  27. CIGRÉ Technical Brochure 63. Guide to Procedures for Estimating the Lightning Performance of Transmission Lines; Working Group 33.01; CIGRÉ: Paris, France, 1991. [Google Scholar]
  28. IEC Technical Report 60071-4. Insulation Co-Ordination—Part 4: Computational Guide to Insulation Co-Ordination and Modelling of Electrical Networks, 1st ed.; Report; IEC: Geneva, Switzerland, 2004. [Google Scholar]
  29. Morched, A.; Gustavsen, B.; Tartibi, M. A universal model for accurate calculation of electromagnetic transients on overhead lines and underground cables. IEEE Trans. Power Deliv. 1999, 14, 1032–1038. [Google Scholar] [CrossRef]
  30. Gustavsen, B.; Irwin, G.; Mangelrød, R.; Brandt, D.; Kent, K. Transmission Line Models for the Simulation of Interaction Phenomena Between Parallel AC and DC Overhead Lines. In Proceedings of the International Conference on Power System Transients, Budapest, Hungary, 20–24 June 1999. [Google Scholar]
  31. Stern, E.; Nash, C.; Shhoeniger, C.; Bartzsch, C.; Acquaotta, G.; Bacchini, M.; Orini, A. The Neptune regional transmission system 500 kV HVDC project. In Proceedings of the Conseil International des Grands Réseaux Electriques, Paris, France, 24–29 August 2008. [Google Scholar]
  32. Bewley, L.V. Travelling Waves on Transmission System, 1st ed.; Dover: New York, NY, USA, 1963; Volume 1. [Google Scholar]
  33. Chowdhuri, P.; Anderson, J.G.; Chisholm, W.A.; Field, T.E.; Ishii, M.; Martinez, J.A.; Marz, M.B.; McDaniel, J.; McDermott, T.E.; Mousa, A.M.; et al. Parameters of lightning strokes: A review. IEEE Trans. Power Deliv. 2005, 20, 346–358. [Google Scholar] [CrossRef]
  34. Mikropoulos, P.N.; Tsovilis, T.E. Lightning attachment models and maximum shielding failure current of overhead transmission lines: Implications in insulation coordination of substations. IET Gener. Transm. Distrib. 2010, 4, 1299–1313. [Google Scholar] [CrossRef]
  35. Modeling of metal oxide surge arresters. IEEE Trans. Power Deliver. 1992, 7, 302–309. [CrossRef]
Figure 1. Fast front model of mixed High Voltage DC (HVDC) transmission line for analysis of lightning transients. A finite cable section is present between two overhead line sections terminating in matching impedance to avoid reflection of voltage surge from remote ends of the transmission line (description of symbols is provided in Table 1).
Figure 1. Fast front model of mixed High Voltage DC (HVDC) transmission line for analysis of lightning transients. A finite cable section is present between two overhead line sections terminating in matching impedance to avoid reflection of voltage surge from remote ends of the transmission line (description of symbols is provided in Table 1).
Energies 12 04217 g001
Figure 2. Overhead section of mixed HVDC transmission line. (a) Structure. (b) Fast front model for electromagnetic transient (EMT) simulation (description of symbols is provided in Table 1 and Equation (1)) [15].
Figure 2. Overhead section of mixed HVDC transmission line. (a) Structure. (b) Fast front model for electromagnetic transient (EMT) simulation (description of symbols is provided in Table 1 and Equation (1)) [15].
Energies 12 04217 g002
Figure 3. Underground section of mixed HVDC transmission line (a) layout, (b) dimensions and electrical properties.
Figure 3. Underground section of mixed HVDC transmission line (a) layout, (b) dimensions and electrical properties.
Energies 12 04217 g003
Figure 4. Upwardly concave front lightning current wave shape [27].
Figure 4. Upwardly concave front lightning current wave shape [27].
Energies 12 04217 g004
Figure 5. Geometric model of tower for application of an electro-geometric model (EGM) to calculate the shielding failure current (IMSF).
Figure 5. Geometric model of tower for application of an electro-geometric model (EGM) to calculate the shielding failure current (IMSF).
Energies 12 04217 g005
Figure 6. Frequency dependent surge arrester model.
Figure 6. Frequency dependent surge arrester model.
Energies 12 04217 g006
Figure 7. Residual voltage across non-linear resistors as a function of the impulse current.
Figure 7. Residual voltage across non-linear resistors as a function of the impulse current.
Energies 12 04217 g007
Figure 8. An equivalent circuit of a mixed HVDC transmission line for analysis of direct lightning strikes on overhead power conductors (OHPCs) at a negative pole (only the two towers closest to the sending end of the cable are shown) (description of symbols is provided in Table 1). [15].
Figure 8. An equivalent circuit of a mixed HVDC transmission line for analysis of direct lightning strikes on overhead power conductors (OHPCs) at a negative pole (only the two towers closest to the sending end of the cable are shown) (description of symbols is provided in Table 1). [15].
Energies 12 04217 g008
Figure 9. Core-ground overvoltage at the entrance of a cable according to the length of the cable.
Figure 9. Core-ground overvoltage at the entrance of a cable according to the length of the cable.
Energies 12 04217 g009
Figure 10. Core-ground overvoltage at ten equally spaced intervals along the length of the cable (3.125–100 km).
Figure 10. Core-ground overvoltage at ten equally spaced intervals along the length of the cable (3.125–100 km).
Energies 12 04217 g010
Figure 11. Maximum core-ground voltage along the length in cables with different lengths.
Figure 11. Maximum core-ground voltage along the length in cables with different lengths.
Energies 12 04217 g011
Figure 12. Sheath-ground overvoltage at the sending end terminal of the cables considering different cable lengths.
Figure 12. Sheath-ground overvoltage at the sending end terminal of the cables considering different cable lengths.
Energies 12 04217 g012
Figure 13. Sheath-ground overvoltage at ten equally spaced intervals along the length of a cable (3.125–100 km).
Figure 13. Sheath-ground overvoltage at ten equally spaced intervals along the length of a cable (3.125–100 km).
Energies 12 04217 g013
Figure 14. Maximum sheath-ground voltage along the length in cables with different lengths.
Figure 14. Maximum sheath-ground voltage along the length in cables with different lengths.
Energies 12 04217 g014
Figure 15. An equivalent circuit of a mixed HVDC transmission line for analysis of a lightning strike on the OHGW/tower top (only the two towers closest to the sending end of the cable are shown) (description of symbols is provided in Table 1).
Figure 15. An equivalent circuit of a mixed HVDC transmission line for analysis of a lightning strike on the OHGW/tower top (only the two towers closest to the sending end of the cable are shown) (description of symbols is provided in Table 1).
Energies 12 04217 g015
Figure 16. Sheath-ground overvoltage at the sending end terminal of the cables considering different cable lengths.
Figure 16. Sheath-ground overvoltage at the sending end terminal of the cables considering different cable lengths.
Energies 12 04217 g016
Figure 17. Sheath-ground overvoltage at ten equally spaced intervals along the length of the cable (3.125–100 km).
Figure 17. Sheath-ground overvoltage at ten equally spaced intervals along the length of the cable (3.125–100 km).
Energies 12 04217 g017
Figure 18. Maximum sheath-ground voltage along the length in cables with different lengths.
Figure 18. Maximum sheath-ground voltage along the length in cables with different lengths.
Energies 12 04217 g018
Figure 19. Location of surge arresters.
Figure 19. Location of surge arresters.
Energies 12 04217 g019
Figure 20. Maximum core-ground overvoltage according to the cable length.
Figure 20. Maximum core-ground overvoltage according to the cable length.
Energies 12 04217 g020
Figure 21. Maximum sheath-ground overvoltage at the entrance of the cable resulting from a lightning strike on the negative OHPC at the transition tower.
Figure 21. Maximum sheath-ground overvoltage at the entrance of the cable resulting from a lightning strike on the negative OHPC at the transition tower.
Energies 12 04217 g021
Figure 22. Maximum sheath-ground overvoltage at the entrance of the cable resulting from a lightning strike on OHGWs at the transition tower.
Figure 22. Maximum sheath-ground overvoltage at the entrance of the cable resulting from a lightning strike on OHGWs at the transition tower.
Energies 12 04217 g022
Figure 23. Core-ground overvoltage along the length of the cable with a surge arrester installed at the terminals.
Figure 23. Core-ground overvoltage along the length of the cable with a surge arrester installed at the terminals.
Energies 12 04217 g023
Figure 24. Sheath-ground (SG) overvoltage along the length of cable with a surge arrester installed at the terminals.
Figure 24. Sheath-ground (SG) overvoltage along the length of cable with a surge arrester installed at the terminals.
Energies 12 04217 g024
Figure 25. Surge arrester energy according to the length of the cable.
Figure 25. Surge arrester energy according to the length of the cable.
Energies 12 04217 g025
Figure 26. Lightning strike on the negative pole OHPC at twr2S. (a) Voltage across negative pole insulators of twr2S and twr3S. (b) Leader propagation between insulator arc horns of twr2S and twr3S.
Figure 26. Lightning strike on the negative pole OHPC at twr2S. (a) Voltage across negative pole insulators of twr2S and twr3S. (b) Leader propagation between insulator arc horns of twr2S and twr3S.
Energies 12 04217 g026
Figure 27. Lightning strike on the top of twr2S. (a) Voltage across positive pole insulators of twr2S according to sheath grounding impedance. (b) Leader growth between insulator gap of twr2S according to sheath grounding impedance.
Figure 27. Lightning strike on the top of twr2S. (a) Voltage across positive pole insulators of twr2S according to sheath grounding impedance. (b) Leader growth between insulator gap of twr2S according to sheath grounding impedance.
Energies 12 04217 g027
Figure 28. Lightning strike on the top of twr2S. (a) Voltage across insulators of the positive pole of twr2S according to footing impedance. (b) Leader propagation between the insulator gap of twr2S according to footing impedance.
Figure 28. Lightning strike on the top of twr2S. (a) Voltage across insulators of the positive pole of twr2S according to footing impedance. (b) Leader propagation between the insulator gap of twr2S according to footing impedance.
Energies 12 04217 g028
Table 1. Values of electromagnetic transient (EMT) simulation model parameters.
Table 1. Values of electromagnetic transient (EMT) simulation model parameters.
SymbolQuantityValue
hHeight of tower47 m
lt1Length of upper tower segment4 m
lt2Length of lower tower segment43 m
ZTCharacteristic impedance of tower segments139 Ω
rfFooting impedance10 Ω
rsgSheath-grounding impedance10 Ω
cinsCapacitance of insulator2.75 pF
LarcInsulator arc inductance5 µH
ltwr_spanLength of tower span350 m
vsheathWave velocity on cable sheath2.7 × 107 m/s
vcoreWave velocity on cable core6.4 × 107 m/s
vOHPCWave velocity on OHPC2.8 × 108 m/s
vOHGWWave velocity on OHGW2.8 × 108 m/s
ZsheathCharacteristic impedance of cable sheath27.1 Ω
ZcoreCharacteristic impedance of cable core72.7 Ω
ZOHGWCharacteristic impedance of OHGW541.8 Ω
ZOHPCCharacteristic impedance of OHPC465.3 Ω
vtwrWave velocity on tower segments2.55 × 108 m/s
lriserLength of riser section50 m
KLeader coefficient1.3 m2/kV2.s
E0Critical leader inception gradient600 kV/m
yTime of lightning strike0.007 s
Table 2. Lightning impulse current parameters.
Table 2. Lightning impulse current parameters.
Parameters of StrikeOHPC Strike OHGW/Tower Top Strike
Ippeak current20 kA200 kA
t30front (30%–90%)3 µs8 µs
tfeffective front time5 µs13.33 µs
tttail time77.5 µs200 µs
Smmaximum steepness20 kA/µs48 kA/µs
nwave-shape constant18.611.171
A//1.061 kA/µs4.325 kA/µs
B//3.257 × 10−17 kA/µs9.247 × 10−14 kA/µs
tn//8.653 µs21.858 µs
t1//99.324 µs257.004 µs
t2//0.1 µs0.416 µs
I1//20.020 kA200.324 kA
I2//2.0201 kA20.324 kA
Table 3. Parameters of the surge arrester model
Table 3. Parameters of the surge arrester model
Components of Surge Arrester Modeln = 1, d = 6 m
L190 µH
R1390 Ω
L01.2 µH
R0600 Ω
C16.66667 pF
Table 4. Maximum core-ground overvoltage in cables with different lengths.
Table 4. Maximum core-ground overvoltage in cables with different lengths.
Cable Length (km)Maximum Overvoltage (p.u.)Distance from Sending End (%)
100−2.1410
50−2.1410
25−2.1410
12.5−2.34510
6.25−3.03210
3.125−3.94810

Share and Cite

MDPI and ACS Style

Asif, M.; Lee, H.-Y.; Park, K.-H.; Shakeel, A.; Lee, B.-W. Assessment of Overvoltage and Insulation Coordination in Mixed HVDC Transmission Lines Exposed to Lightning Strikes. Energies 2019, 12, 4217. https://doi.org/10.3390/en12214217

AMA Style

Asif M, Lee H-Y, Park K-H, Shakeel A, Lee B-W. Assessment of Overvoltage and Insulation Coordination in Mixed HVDC Transmission Lines Exposed to Lightning Strikes. Energies. 2019; 12(21):4217. https://doi.org/10.3390/en12214217

Chicago/Turabian Style

Asif, Mansoor, Ho-Yun Lee, Kyu-Hoon Park, Ayesha Shakeel, and Bang-Wook Lee. 2019. "Assessment of Overvoltage and Insulation Coordination in Mixed HVDC Transmission Lines Exposed to Lightning Strikes" Energies 12, no. 21: 4217. https://doi.org/10.3390/en12214217

APA Style

Asif, M., Lee, H. -Y., Park, K. -H., Shakeel, A., & Lee, B. -W. (2019). Assessment of Overvoltage and Insulation Coordination in Mixed HVDC Transmission Lines Exposed to Lightning Strikes. Energies, 12(21), 4217. https://doi.org/10.3390/en12214217

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop