Overhead Line Lightning Strike Severity and Probability

October 15, 2009

Lightning Phenomenon:

Lightning Strike is the discharge of electric charge accumulated in the clouds to the ground. Clouds accumulate typically the negative charge, and in response the ground produces the counter charge: the positive charge. Electric discharge happens between the two, i.e. clouds and the ground, under certain electric discharge conducive conditions, effectively creating the shortest electrical path.  These conditions include accumulated charge density, humidity in the air (enabling faster dielectric break down), ground elevation (buildings, mountains, tall living things) etc.

Physics behind this phenomenon is much elaborated in several books and articles found in both libraries and internet, and their perusal is recommended.

Lightning Strike Frequency:

Although the frequency of lightning strikes might vary year to year, long term (usually many years to decades) maintained records give a statistical approximate number of expected strikes each year. This is usually referred as Isokeraunic Level [1]. Reference [2] and [3] give typical Isokeraunic Levels around the world, and gives a geography dependent empirical formulae to calculate Ground-Flash Density (or – Strike Density). The empirical relationship between Isokeraunic Level and Ground-Flash Density for UK is given as:

GFD = aT_{b}

Where, GFD is the Ground-Flash Density in Flashes/km^2/Yr, a is a factor that varies between (2.6\pm 0.2)10^{-3}, b is a factor that varies between (1.9\pm 0.1), and T is the Isokeraunic Level in Flashes/km^2/Yr. Isokeraunic Level for Ireland, north UK, mid-west UK  and southern-west UK is typically between 5-10, while for middle-east UK and southern-east UK is between 10-20.

Overhead Hit Frequency:

Strike radius is the scope of the ground elevated structures or ground itself that will be exposed to lightning. This, usually calculated in meters, is very much dependent of lightning strike peak current magnitude.

Lightning strike distance in meters [3], R_{s}=8I_{s}^{0.65}

Where, I_{s} is the Lightning Strike peak magnitude in kA.

Now, the scope of Lightning Strike area that will hit the OHL will depend on the tower structure. For a un-shielded wood pole structure (assuming all conductors are same height from the ground level), this is given as following:

Strike Area, A_{s}=(R_{s}+D_{c}+R_{s})\cdot L

Where, (R_{s}+D_{c}+R_{s}) is the sum of strike radius of left most conductor, distance between the farthest conductors, and strike radius of the right most conductor respectively. $latex L$ is the length of the exposed overhead line.

Total estimated OHL Flashes per Year if every flash was as chosen Strike Current, I_{s}, is given as, N=GFD\cdot A_{s}

For a given strike current, Cumulative Lightning Strike Probability is given as [3],

P_{c}=\frac{1}{[1+(\frac{I_{s}}{31})^{2.6}]}

For a given strike current, Estimated Annual OHL Lightning Strikes in (Flashes/Yr) = N\cdot P_{c}

Now plot a bar chart between various strike current versus expected annual lightning strike frequency, and you see typical decline of Flashes/Yr with increasing strike current magnitude. And if you invert the Flashes/Yr you get expected number of years before you get see a lightning strike for a given strike current peak.

References:

  1. Lightningtech Website. [online], Available: http://www.lightningtech.com/d~ta/faq2.html, Accessed: Oct 2009.
  2. Chowddhuri, P., Electromagnetic Transients in Power Systems, Exeter: Research Strudies Press Ltd., 1996.
  3. Grigsby, L. L., Power Systems, Boca Raton: CRC Press, 2007.

Per Unit Formulae and Conversions

September 19, 2009

Even the best of the power systems Geeks sometime stumble at the basics. Here are some per unit formulae to brush up if you are one of them or among those willing to learn like me.

A per unit value is a scaled factor of the measured system parameter to the chosen system base value. In other words, the chosen value becomes by which all other values for consider system parameter will be compared. It’s a ratio and is a unit less quantity.

Most common approach, for most studies, is to choose a system apparent power base (typically 100MVA) and system local network voltage (line to line voltage) base values and base values for all other parameters can be estimated there after.

Base Value Calculation:

I_{base}=\frac{S_{base}}{\sqrt{3}V_{base}}          (1)

Z_{base}=\frac{V_{base}^2}{S_{base}}=\frac{V_{base}}{\sqrt{3}I_{base}}         (2)

Y_{base}=\frac{1}{Z_{base}}         (3)

Per Unit Calculation:

I_{pu}=\frac{I}{I_{base}}         (4)

Z_{pu}=\frac{Z}{Z_{base}}         (5)

B_{pu}=\frac{2\pi fC}{Y_{base}}          (6)

Per Unit Conversion:

Z_{new}=Z_{old}\left (\frac{V_{old}}{V_{new}}\right )^2\left (\frac{S_{new}}{S_{old}}\right )         (7)

References:

  1. Wikipedia, Per Unit System, Wikipedia. [online]. Available: http://en.wikipedia.org/wiki/Per-unit_system. [Accessed: Sep. 19, 2009].
  2. W. D. Stevenson, Jr., Elements of Power System Analysis, 3rd ed. New York: McGraw-Hill, 1994.

Power Factor Capacitor Discharge Resistor Design

August 27, 2009

When a Power Factor Capacitor (PFC) step, a single unit in a bank of several capacitors, is disconnected or switched off, the discharging resistor connected across the capacitor will discharge it to designed retained voltage V value in discharge time t seconds. Typical discharge resistor ratings for a given power factor capacitor C  in μF or MVAR include: maximum normal operation system voltage V_0  in kV and required retained discharge voltage V  in Volts (around 50V) at discharge time t  in seconds (usually <60s). The discharge resistor R  in kΩ  is given as,

R=\left | \frac{-t}{ln(\frac{V}{V_0})C} \right |          (1)

Depending whether the PFC bank step is star or delta connected (usually star connected), for given system frequency f  and capacitor step’s rated reactive power output Q in MVAR, the capacitor step’s capacitance C  in μF is given as:

  • If star connected, C=\frac{Q}{2 \pi f V_0^{2}}
  • If delta connected, C=\frac{3Q}{2 \pi f V_0^{2}}

Explanation:

Capacitor voltage decay across the resistor is given as, V=V_{0}e^{\frac{-t}{RC}}. Rearrange this for Discharge Resistor ‘R’ and you get (1). Simple!


Power Factor Capacitor Detuning

August 17, 2009

Power Factor Capacitor (PFC) banks are commonly detuned to specific frequencies to reduce or remove harmonic resonances. This is done by introducing a detuning reactor in series with the capacitor cluster at each three phase leg.

At the detuning frequency, detuning reactor will have the same impedance as the PFC capacitance per phase leg, i.e. X_{L}=X_{C}.

Detune reactor inductance is then given as: L=\frac{1}{(2\Pi f_{d})^{2}C}Henry

Where, f_{d} is the PFC detune frequency and C is the PFC capacitance per phase leg.


R and X Values from given Grid Parameters

August 8, 2009

Wondering how to calculate R and X values from given X/R ratio, system voltage and fault level (kA or MVA) information at a power system node or bus? Here’s how.

Lets define, X/R ratio as: \frac{X}{R}=n.

Then the node or bus thevenin equivalent resistance and reactance are given as:

R=\frac{V^{2}}{S\sqrt{n^{2}+1}}, X=R\times n

Where, ‘R’ is the thevenin equivalent resistance in Ω, ‘X’ is the thevenin equivalent reactance in Ω, ‘V’ is the system voltage is kV, ‘S’ RMS break fault power in MVA, ‘I’ RMS break fault current in kA.

RMS break fault level power and current relation is given as: S=\sqrt{3}\times V \times I

Download spreadsheet to do the Math.

Explanation:

Thevenin’s equivalent impedance is given as,

\left | Z \right |=\sqrt{R^{2}+X^{2}}         →          (1)

Substituting X=R \times n in (1) gives,

R\sqrt{n^{2}+1}           →          (2)

Also, \left | Z \right |=\frac{V^{2}}{S}        →          (3)

From (1) and (3), R\sqrt{n^{2}+1}=\frac{V^{2}}{S}

=>R=\frac{V^{2}}{S\sqrt{n^{2}+1}}         →          (4)

And X=R\times n          →          (5)


Follow

Get every new post delivered to your Inbox.