Research Article | Open Access
Calculation of Transient Potential Rise on the Wind Turbine Struck by Lightning
A circuit model is proposed in this paper for calculating the transient potential rise on the wind turbine struck by lightning. The model integrates the blade, sliding contact site, and tower and grounding system of the wind turbine into an equivalent circuit. The lightning current path from the attachment point to the ground can be fully described by the equivalent circuit. The transient potential responses are obtained in the different positions on the wind turbine by solving the circuit equations. In order to check the validity of the model, the laboratory measurement is made with a reduced-scale wind turbine. The measured potential waveform is compared with the calculated one and a better agreement is shown between them. The practical applicability of the model is also examined by a numerical example of a 2 MW Chinese-built wind turbine.
As wind power generation undergoes rapid growth, wind turbines (WTs) erected in large wind farms are regarded as a realistic alternative to conventional power plants. Due to their large height of towers, long rotating blades, and exposed locations, WTs can be easily struck by lightning. When a WT suffers a lightning stroke, a large lightning current will continue to flow from the attachment point to the grounding system and cause the transient potential rise on the WT. The high potential can sparkover the insulation distance and severely damage the WT components inside the structure. It can also give rise to the overvoltage surge that endangers the sensitive electronics and microprocessors in the WT control system. In the light of the serious problem arising from the transient potential rise, the protection design for the potential rise phenomenon has been paid more and more attention. In the protection design, the need exists for calculating the transient potential responses on the WT under lightning stroke. A few calculation methods were reported in literature [1–3]. In the previous methods, the tower, that is, the longest traveling path of lightning current on the WT, was modeled as a uniform transmission line. They are simple and easy to be adopted in the lightning transient calculation but incapable of calculating the transient potential responses in the different positions on the WT owing to their neglecting of the structural feature of the tower body. Recently, a more detailed method has been presented. This method represents the tower as a cage-like multiconductor system [4, 5] and can give the potential distribution on the WT; however, it causes a significant increase in the circuit complexity. For an improvement on the lightning transient calculation of WTs, a simplified circuit model is proposed in this paper. The proposed model divides the tower into a series of hollow cylindrical sections and represents it as a π-type circuit chain instead of the complicated cage-like multiconductor system. Based on the simplified treatment for the tower, a WT is converted into an equivalent circuit that can give a complete description of lightning current path including the blade, sliding contact site, and tower and grounding system. Then, transient potential responses are obtained in the different positions on the WT by performing the transient calculation for the equivalent circuit. In order to confirm the validity of the proposed model, a laboratory experiment has also been carried out on a reduced-scale WT.
2. Circuit Parameters
An internal down conductor is usually installed inside the blade for conducting the lightning current from the attachment point to the hub, as shown in Figure 1. To take account of the propagation phenomenon of lightning current, the down conductor is divided into a suitable number of segments. The segment length needs to fulfill the following criterion : where is the velocity of light and the maximum frequency likely to affect the system transient. An arbitrary segment in Figure 1 can be represented by a -type circuit composed of resistance, inductance, and capacitance, as shown in Figure 2. The resistance is estimated by  where and are the material permeability and conductivity of the down conductor, respectively, and is the conductor radius. The inductance is calculated by  where is the permeability of free space ( H/m) and the two parameters in the round brackets are where . According to the electromagnetic analogy [8–10], the capacitance is determined by where is the permittivity of free space [(36π)−1 × 10−9 F/m].
2.2. Sliding Contact Site
The conductive path between the blade root and tower top is a sliding contact site and mainly includes the brushes, sliding contact systems, and main shaft bearings. The brushes have been widely used in the multimegawatt WTs to divert the lightning current from the blade root to the tower top. The circuit parameters of the brushes and main shaft bearings are simply represented as the contact resistance and equivalent capacitance [3, 5], respectively. can be evaluated by the formula given in . Since the nacelle is rarely tuning, the yaw bearing is not considered here. In view of the actual shunting route of lightning current, the sliding contact site is modeled as a simple parallel circuit, as shown in Figure 3.
2.3. Tower and Grounding System
An actual tower takes the shape of the hollow circular truncated cone, as shown in Figure 4, and is the longest traveling path of lightning current. As is the case in the blade, consideration of the propagation phenomenon of the lightning current needs the tower to be divided into a certain number of sections. Each section is approximately taken as a hollow cylinder and its length is determined by (1). An arbitrary section () in Figure 4 can also be represented by a π-type circuit, as shown in Figure 5. The resistance is roughly evaluated by [12, 13] where and is the DC resistance where and are the material permeability and conductivity of the tower, respectively. The inductance is given by  where is the geometrical factor depending on the ratio of . Table 1 gives the values of in the range of . By substituting into (5), the capacitance can be obtained.
3. Circuit Model
After obtaining the circuit parameters of the blade, sliding contact site, and tower and grounding system, a complete equivalent circuit can be built for a WT, as shown in Figure 6. The lightning current source is injected to the top node of the equivalent circuit to simulate a lightning stroke to the blade tip and the impedance in parallel with is the surge impedance of the lightning channel. For the capacitance in the equivalent circuit, as shown in Figure 7(a), its circuit equation is written as Integration of (9) between gives where is the time step size. By means of the trapezoidal rule, the integration is evaluated by This can be rewritten as where The circuit model for (12) is a parallel circuit unit, as shown in Figure 7(b). For the resistance-inductance (-) branch in the equivalent circuit, as shown in Figure 8(a), its circuit equation is that is Using the trapezoidal rule to integrate both sides of (15) between yields  Rearranging the terms in (16) leads to Where In terms of (17), the parallel circuit unit is depicted for the - branch, as shown in Figure 8(b). The current sources and are known from the values for the preceding time step. After all capacitances and - branches are replaced by their respective parallel circuit units, Figure 6 is further converted into a time discretization circuit consisting only of resistances and current sources, as shown in Figure 9. The node voltage equations are set up for the time discretization circuit where is the node conductance matrix, the node voltage vector, and the node current source vector. At , all capacitances and inductances are set in zero initial conditions. The node voltage vector is obtained by solving (19) in each time step [17, 18]. Thus, the transient potential responses can be given for the different positions on the WT under lightning stroke.
4. Experimental Verification
A reduced-scale WT is built in the high voltage laboratory. Its dimensions are shown in Figure 10. The grounding resistance is 4 Ω. The fast impulse current provided by an impulse generator is injected to the blade tip. The lead wire for measuring the transient potential is stretched perpendicular to the current lead wire and grounded (connected to the steel plate) at a point 8.0 m apart from the reduced-scale WT, which can restrain the electromagnetic induction between the two wires. In the experimental measurement, the current and potential signals are recorded by a digital oscilloscope. The measured waveforms of injected current and transient potential at the tower top are shown in Figure 11, where the corresponding waveform calculated from the circuit model proposed above is also given for comparison. It can be seen from Figure 11 that the calculated waveform is close to the measured one.
5. Calculated Results
A 2MW Chinese-built WT is considered here. The dimensions of the WT are m (see Figure 1), m, m, m, and m (see Figure 4). The grounding resistance is 3 Ω and the lightning current is taken as 10/350 μs, 100 kA (see Figure 6) according to Chinese design standard [15, 19]. The peak potential distribution on the WT is plotted in Figure 12. The transient potential waveforms in three typical positions on the WT are also given in Figure 13. As seen from Figures 12 and 13, the transient potential rise on the WT is very serious and can do damage to the components and equipment inside the structure during a lightning stroke.
Calculation of the transient potential rise has been performed for the WT under lightning stroke. The lightning current path on the WT is described by a proposed circuit model that integrates the blade, sliding contact site, and tower and grounding system into a complete equivalent circuit. The model can predict the transient potential levels in the different positions on the WT. The calculated results obtained from the model indicate that the transient potential rise is serious and very harmful to the components inside the structure during a lightning stroke. The validity of the model has been confirmed by the laboratory experiment on a reduced-scale WT. The model is useful in lightning transient analysis of WTs and has the capability of providing a sound basis for the lightning protection design of WTs.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
This work was financially supported by National Natural Science Foundation of China under Award no. 509770926.
- R. B. Rodrigues, V. M. F. Mendes, and J. P. S. Catalão, “Analysis of transient phenomena due to a direct lightning strike on a wind energy system,” Energies, vol. 5, no. 7, pp. 2545–2558, 2012.
- T. Q. Nguyen1, T. Pham, and T. V. Tran, “Electromagnetic transient simulation of lightning overvoltage in a wind form,” in Proceedings of the Electrical Insulation Conference, Ottawa, Ontario, 2013.
- D. Romero, J. Montanyà, and A. Candela, “Behavior of the wind-turbines under lightning strikes including nonlinear grounding system,” in Proceedings of the International Conference on Renewable Energy and Power Quality, Barcelona, Spain, March-April 2004.
- X. Q. Zhang and C. Liu, “Lightning transient modeling of wind turbine towers,” International Review of Electrical Engineering, vol. 7, no. 1, pp. 3505–3511, 2012.
- X. Q. Zhang, Y. Z. Zhang, and C. H. Liu, “A complete model of wind turbines for lightning transient analysis,” Journal of Renewable and Sustainable Energy, vol. 6, no. 1, Article ID 013113, 2014.
- G. Celli and F. Pilo, “EMTP models for current distribution evaluation in LPS for high and low buildings,” in Proceedings of the 25th International Conference on Lightning Protection, Rhodes, Greece, September 2000.
- R. B. Standler, Protection of Electronic Circuits from Overvoltages, John Wiley & Sons, New York, NY, USA, 1989.
- C. Z. Feng, Electromagnetic Field Theory, Higher Education Press, Beijing, China, 2002.
- A. Ametani, Y. Kasai, J. Sawada, A. Mochizuki, and T. Yamada, “Frequency-dependent impedance of vertical conductors and a multiconductor tower model,” IEE Proceedings—Generation, Transmission and Distribution, vol. 141, no. 4, pp. 339–345, 1994.
- U. Y. Iosseli, A. S. Kothanof, and M. G. Stlyrski, Handbook of Capacitance Calculation, Energy Press, Moscow, Russia, 1987 (Russian).
- M. Paolone, F. Napolitano, A. Borghetti et al., “Models of wind-turbine main shaft bearings for the development of specific lightning protection systems,” in IEEE Lausanne Power Tech, pp. 1–5, Lausanne, Switzerland, July 2007.
- A. H. M. Arnold, “The alternating-current resistance of tubular conductors,” Journal of the Institution of Electrical Engineers, vol. 78, no. 473, pp. 580–596, 1936.
- V. T. Morgan, R. D. Findlay, and S. Derrah, “formula to calculate the skin effect in isolated tubular conductors at low frequencies,” IEE Proceedings—Science, Measurement and Technology, vol. 147, no. 4, pp. 169–171, 2000.
- B. L. Kalentrarov and L. A. Zeitlin, Handbook of Inductance Calculation, Energy Press, Moscow, Russia, 1992 (Russian).
- Chinese Industrial Standard of Energy Source, “Technical specification for lightning protection system of wind turbines,” Tech. Rep. NB/T 31039-2012, Standards Press of China, Beijing, China, 2012.
- H. W. Dommel, Electromagnetic Transients Program Theory Book, BPA, Portland, Ore, USA, 1995.
- L. V. D. Sluis, Transients in Power Systems, John Wiley & Sons, New York, NY, USA, 2001.
- W. Shi, Overvoltage Calculation in Power Systems, High Education Press, Beijing, China, 2009.
- Chinese National Standard; GB 50057-2010, Design Code for Lightning Protection of Structures, China Planning Press, Beijing, China, 2010.
Copyright © 2014 Zhang Xiaoqing. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.