Issue 
J. Space Weather Space Clim.
Volume 4, 2014
Space Weather and Challenges for Modern Society



Article Number  A03  
Number of page(s)  7  
DOI  https://doi.org/10.1051/swsc/2013057  
Published online  09 January 2014 
Research Article
Observations and modeling of GIC in the Chinese largescale highvoltage power networks
^{1}
School of Electrical and Electronic Engineering, North China Electric Power University, Beijing, China
^{2}
Geomagnetic Laboratory, Natural Resources Canada, Ottawa, Canada
^{3}
Finnish Meteorological Institute, Helsinki, Finland
^{*} Corresponding author: LiuChunming@ncepu.edu.cn; Chunming.Liu@NRCanRNCan.gc.ca
Received:
28
March
2013
Accepted:
17
December
2013
During geomagnetic storms, the geomagnetically induced currents (GIC) cause bias fluxes in transformers, resulting in halfcycle saturation. Severely distorted exciting currents, which contain significant amounts of harmonics, threaten the safe operation of equipment and even the whole power system. In this paper, we compare GIC data measured in transformer neutrals and magnetic recordings in China, and show that the GIC amplitudes can be quite large even in midlow latitude areas. The GIC in the Chinese Northwest 750 kV Power Grid are modeled based on the plane wave assumption. The results show that GIC flowing in some transformers exceed 30 A/phase during strong geomagnetic storms. GIC are thus not only a highlatitude problem but networks in middle and low latitudes can be impacted as well, which needs careful attention.
Key words: electric circuit / geomagnetically induced currents (GIC) / modelling / engineering / space weather
© C. Liu et al., Published by EDP Sciences 2014
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
During strong space weather storms, which are caused by the activity of the Sun, the Earth’s magnetic field is intensely disturbed by the space current system in the magnetosphere and ionosphere. The electric fields induced by time variations of the geomagnetic field drive geomagnetically induced currents (GIC) in electric power transmission networks. The frequencies of GIC are in the range of 0.0001 ~ 0.1 Hz. Such quasiDC currents cause bias fluxes in transformers, which result in halfcycle saturation due to the nonlinear response of the core material (e.g., Kappenman & Albertson 1990; Molinski 2002; Kappenman 2007). The sharply increased magnetizing current with serious waveform distortion may lead to temperature rise and vibration in transformers, reactive power fluctuations, voltage sag, protection relay malfunction, and possibly even a collapse of the whole power system (e.g., Kappenman 1996; Bolduc 2002).
Large GIC are usually considered to occur at high latitudes such as North America and Scandinavia, where tripping problems and even blackouts of power systems due to GIC have been experienced (Bolduc 2002; Pulkkinen et al. 2005; Wik et al. 2009). Large currents in transformer neutrals have been monitored in the Chinese highvoltage power system many times during geomagnetic storms although China is a midlowlatitude country. At the same time, transformers have had abnormal noise and vibration. Those events have been shown to be caused by GIC based on analyses of simultaneous magnetic data and GIC recordings (Liu & Xie 2005; Liu et al. 2009a). The power grids are using higher voltages, longer transmission distances, and larger capacity with the developing economy in China. So, the risk that the power systems would suffer from GIC problems may obviously increase. The Chinese Northwest 750 kV power grid has long transmission lines with small resistances making it prone to large GIC during geomagnetic storms. Thus it is important to model GIC particularly in that network.
2. GIC observations in Chinese highvoltage power grid
We acquire GIC data through the neutral point of the transformer at the Ling’ao nuclear power plant (22.6° N, 114.6° E) in the Guangdong Province. Besides, geomagnetic field data are collected from the Zhaoqing Geomagnetic Observatory (23.1° N, 112.3° E) which is not very far from Ling’ao. Figure1 shows the neutral point current (top panel), the horizontal component of the geomagnetic field (bottom panel), and its variation rate (middle panel) during the magnetic storms on 7–8 (a) and 9–10 (b) November 2004. The occurrence times of the current peaks match with those of the geomagnetic field variation rate. It is confirmed that there is no HVDC (highvoltage direct current) monopole operation during that time. So it is reasonable to believe that the currents are really GIC induced by geomagnetic storms. The maximum value of GIC is up to 75.5 A/3 phases, which is much higher than the DC bias caused by monopole operation of HVDC.
Fig. 1. GIC data at the Ling’ao nuclear power plant on 7–8 (a) and 9–10 (b) November 2004. The horizontal component of the geomagnetic field and its variation rate are also shown based on data from the Zhaoqing Geomagnetic Observatory. 
3. Modeling GIC in power grids
The modeling of GIC in a power grid can be divided into two steps (e.g., Pirjola 2000): step 1, calculating the geoelectric field induced by a magnetic storm; step 2, calculating the GIC in the power grid. The effect of the induced geoelectric field is equivalent to voltage sources in the transmission lines, which enables converting the GIC calculation into a circuit problem in step 2.
3.1. Calculating the electric field using a layered earth model
We use the standard conventional Cartesian geomagnetic coordinate system in which the x, y and z axes point northwards, eastwards, and downwards, respectively. According to the plane wave assumption (e.g., Boteler 1999), the relation between perpendicular horizontal components of the geoelectric (E) and geomagnetic (B) fields at the earth’s surface can be expressed as(1) (2)where μ _{ 0 } is the vacuum permeability and Z is surface impedance of the earth which depends on the conductivity structure of the earth and on the angular frequency ω.
In a previous study about GIC in China, Liu et al. (2009b) used a uniform halfspace model for the earth. However, onedimensional layered earth models are more accurate descriptions for the real situations. Figure 2 shows a layered earth model which contains n layers with conductivities σ _{1}, σ _{2}, …, σ _{ n } and thicknesses h _{1}, h _{2}, …, h _{ n }→∞.
Fig. 2. Layered Earth model for calculating the induced geoelectric field. 
The thickness of the bottom layer is h _{ n }→∞, and E _{ x } = 0 and B _{ y } = 0 when z→∞. Hence the impedance at the top of the layer of the nth layer is(3)where k _{ n } is the propagation constant given by . The impedance at the top of the layer within the mth layer (m = 1, 2, …, n − 1) can be expressed as(4)where and and .
In the model, the bottom of mth layer is the top of (m + 1)th layer, so equation (4) can be seen as a recursive formula for the impedance at the top of each layer, through which we can calculate the surface impedance of the Earth Z. The geoelectric field in frequency domain can be calculated from geomagnetic data according to equations (1) and (2). Then the result has to be inverse Fourier transformed back to the time domain.
3.2. Calculating GIC
The frequencies of GIC are very low from the view point of power systems. Thus the GIC can be treated as a direct current. The effect of the geoelectric field on a power grid is equivalent to a set of voltage sources in the transmission lines between the substations. The value of the voltage is the integral of the electric field along the line, i.e.:(5)
If the geoelectric field is uniform, the integrals are independent of the paths. Therefore equation (5) can be simplified to(6)Where L _{ AB } is the direct distance between nodes A and B; θ is the “compass angles” i.e. clockwise from geographic North.
The GIC flowing from the power grid to the earth can be expressed as a column matrix I, which has the following formula (e.g., Pirjola & Lehtinen 1985)(7)where 1 is a unit (identity) matrix; Y and Z are the network admittance matrix and the earthing impedance matrix respectively. The elements of column matrix J are defined by(8)
The matrix J gives the GIC between the power grid and the earth in the case of ideal groundings, i.e. the grounding resistances are zero making Z a zero matrix.
4. Modeling GIC in Chinese Northwest 750 kV power grid
The problem of GIC should be considered more serious in the Chinese Northwest 750 kV power grid because of the high voltage implying low transmission line resistances and because of the low earth conductivity increasing geoelectric field values. The power grid (shown in Fig. 3) for which GIC calculations are made in this paper is mainly located in the Gansu Province in the Northwest of China. We ignore the lower voltage part connected to the 750 kV power grid when modeling the GIC, because the resistances of that part are much larger, and so it is considered to have little influence on GIC flowing in the 750 kV system.
Fig. 3. Chinese Northwest 750 kV power grid. Three geomagnetic observatories (GRM, LZH, and JYG) are also shown on the map. (The WMQ observatory is not located in the area of this map.) 
4.1. Geoelectric field calculation
We use data of the geomagnetic storm on 29–30 May 2005. The power grid is very large, extending more than 2 000 km in an eastwest direction and 1 500 km in a NorthSouth direction, so the geomagnetic variations cannot be considered to be the same all over the network. The magnetic data from four geomagnetic observatories, whose locations are shown in Figure 3 and in Table 1, are used to calculate the geoelectric field. The local magnetic data are interpolated by using the spherical elementary current systems (SECS) method (Amm 1997). The method uses geomagnetic field data to inverse the ionosphere equivalent current according to which the geomagnetic field data of every location can be calculated. Therefore the interpolation of magnetic data at different locations during a storm can be acquired. As examples, Figure 4a shows the measured data from JYG and the SECSderived magnetic data for Jiuquan Substation, and Figure 4b shows the measured data from LZH and the SECSderived magnetic data for Yongdeng Substation on 29–30 May 2005. It can be seen that the differences between measured magnetic data and the SECSderived data are little except for the base line values which have no effect on the induced electric fields.
Fig. 4. Measured magnetic data and the SECSderived magnetic data on 29–30 May 2005. The horizontal axis is the UT time in hours (a) magnetic data from JYG observatory and the SECSderived magnetic data for Jiuquan substation and (b) magnetic data from LZH observatory and the SECSderived magnetic data for Yongdeng substation. 
Locations of geomagnetic observatories in the area of the Chinese Northwest 750 kV power grid.
The earth conductivities are quite different across the power grid considered, so the geoelectric field values are calculated segment by segment according to the local magnetic data and the local layered earth model. In other words, we utilize the piecewise layered earth model. The earth resistivity in the region where the Chinese Northwest 750 kV power grid is located was provided by Prof. Liu GuoXing, a geologist at the Jilin University (private communication). Figure 5 shows a section of the earth resistivity in Ω·m from Xining to Yinchuan along the 750 kV power lines (see Fig. 3). The resistances of some places are given within a range such as 500–570 at Yinchuan in Figure 5. The upper limit values were used to calculate the induced electric fields because they stand for the most disadvantageous situation to the power grid.
Fig. 5. Resistivity for the section XiningYinchuan along 750 kV power transmission lines. 
As mentioned, the geoelectric fields have been calculated all over the Chinese Northwest 750 kV system based on the Piecewise layered earth models during the geomagnetic storm on 29–30 May 2005. As examples, Figure 6 shows the geoelectric field at Jiuquan and Yongdeng (whose locations are shown in Fig. 3). Our calculation results indicate that the largest E _{ x } value is 0.36 V/km and the largest E _{ y } value is 0.668 V/km in the area of the Northwest 750 kV grid during the geomagnetic storm considered. It is also shown by Figure 6 that the electric fields calculated for Yongdeng and Jiugan are quite different because the Earth conductivity at Yongdeng is much lower than that at Jiuquan.
Fig. 6. Calculated geoelectric fields at two sites (Jiuquan and Yongdeng) of the Chinese Northwest 750 kV power grid on 29–30 May 2005. The horizontal axis is the UT time in hours (a) EJiuquan and (b) EYongdeng. 
4.2. GIC calculation
The GIC through all neutral points of the transformers to the Earth and in all transmission lines of the Chinese Northwest 750 kV network have been calculated. Figure 7 shows the GIC through two typical substations: Jiuquan and Yondeng (also referred to in Fig. 6). The largest GIC at Jiuquan is 25.08 A/phase at 21:35 UT on 29 May 2005, and the largest GIC at Yongdeng is 38.63 A/phase at 22:51 UT on 29 May 2005.
Fig. 7. Calculated GIC at two sites (Jiuquan and Yongdeng) of the Chinese Northwest 750 kV power grid on 29–30 May 2005. The horizontal axis is the UT time in hours (a) calculated GIC at Jiuquan substation and (b) calculated GIC at Youngdeng substation. 
As snapshots, Figure 8 shows the GIC through every node and line at 21:35 UT (panel a) and at 22:51UT (panel b) on 29 May 2005 when the GIC through some of the nodes reach their peaks. It can be seen that the largest GIC through a neutral point is 38.63 A/phase, which is obtained at theYongdeng substation at 22:51 as already mentioned above (see also Fig. 7). The peak GIC through a transmission line is 68.84 A/phase, which occurs in the line from Yongdeng to Jinchang at 21:35 UT. It should be note that there is one singlephase transformer bank in a 750 kV substation except Guanting and Yinchuan where the number of transformer banks is two.
Fig. 8. Snapshots at 21:35 on 29 May 2005 (a) and at 22:51UT on 29 May 2005 (b) of calculated GIC at different sites of the Chinese Northwest 750 kV power grid. The solid circle represents that the GIC flow into the power network from ground, the hollow one means that GIC flow into the ground. The dashed line with an arrow represents the direction of electric field at that substation (a) calculated GIC results at 21:35UT 29 May 2005 and (b) calculated GIC results at 22:51UT 29 May 2005. 
5. Conclusions
The highvoltage power grid in China may experience large GIC during geomagnetic storms, which has been concluded from monitoring the current through the neutral point at Ling’ao nuclear power plant. The GIC in the Chinese Northwest 750 kV power grid during a specific geomagnetic storm have been modeled based on calculating the geoelectric field using the piecewise layered earth models. It can be seen from the results that some sites are sensitive to geomagnetic storms, and the magnitude of GIC can be quite large (> 30 A/phase) during strong geomagnetic storms. Our studies thus clearly demonstrate that GIC are not only a highlatitude problem but networks in middle and low latitudes can be impacted as well. Factors increasing GIC risks in China include the large size of the power network, the small resistances of the transmission lines, and the high resistivity of the earth.
Acknowledgments
This work was supported by the Nature Science Foundation of Hebei Province (E2011502085), National High Technology Research and Development of China 863 Program (2012AA121005).
References
 Amm, O., Ionospheric elementary current systems in spherical coordinates and their application, J. Geomag. Geoelectr., 49, 947–955, 1997. [Google Scholar]
 Bolduc, L., GIC observations and studies in the HydroQuébec power system, J. Atmos. Sol. Terr. Phys., 64 (16), 1793–1802, 2002. [Google Scholar]
 Boteler, D.H., Calculating the voltages induced in technological systems during a geomagnetic disturbance, IEEE Trans. Electromagn. Compat., 41 (4), 398–402, 1999. [CrossRef] [Google Scholar]
 Kappenman, J.G., Geomagnetic storms and their impact on power systems, IEEE Power Eng. Rev., 16 (5),5–8, 1996. [Google Scholar]
 Kappenman, J.G., Geomagnetic disturbances and impacts upon power system operation. In: The Electric Power Engineering Handbook, L.L., Grigsby (Ed.), 2nd Edn., Chapter 16, 161 – 1622, CRC Press/IEEE Press, 2007. [Google Scholar]
 Kappenman, J.G., and V.D. Albertson, Bracing for the geomagnetic storms, IEEE Spectrum, The Electric Power Engineering Handbook, 27 (3), 27–33, 1990. [Google Scholar]
 Liu, L.Y., and X.W. Xie, Analysis of increase of noise of 500 kV transformer (in Chinese), High Voltage Engineering, 31 (4), 85–87, 2005. [Google Scholar]
 Liu, C.M., L.G. Liu, and R. Pirjola, Geomagnetically induced currents in the highvoltage power grid in China, IEEE Trans. Power Delivery, 24 (4), 2368–2374, 2009a. [Google Scholar]
 Liu, C.M., L.G. Liu, R. Pirjola, and Z.Z. Wang, Calculation of geomagnetically induced currents in mid to lowlatitude power grids based on the plane wave method: A preliminary case study, Space Weather, 7 (4), S04005, DOI: 10.1029/2008SW000439, 9 pp., 2009b. [CrossRef] [Google Scholar]
 Molinski, T.S., Why utilities respect geomagnetically induced currents, J. Atmos. Sol. Terr. Phys., 64 (16), 1765–1778, 2002. [Google Scholar]
 Pirjola, R., Geomagnetically Induced Currents during Magnetic Storms, IEEE Trans. Plasma Sci., 28 (6), 1867–1873, 2000. [Google Scholar]
 Pirjola, R., and M. Lehtinen, Currents produced in the Finnish 400 kV power transmission grid and in the Finnish natural gas pipeline by geomagnetically induced electric fields, Ann. Geophys., 3 (4), 485–491, 1985. [Google Scholar]
 Pulkkinen, A., S. Lindahl, A. Viljanen, and R. Pirjola, Geomagnetic storm of 29–31 October 2003: Geomagnetically induced currents and their relation to problems in the Swedish highvoltage power transmission system, Space Weather, 3 (8), S08C03, DOI: 10.1029/2004SW000123, 19 pp., 2005. [Google Scholar]
 Wik, M., R. Pirjola, H. Lundstedt, A. Viljanen, P. Wintoft, and A. Pulkkinen, Space weather events in July 1982 and October 2003 and the effects of geomagnetically induced currents on Swedish technical systems, Ann. Geophys., 27 (4), 1775–1787, 2009. [Google Scholar]
Cite this article as: Liu C, Li Y & Pirjola R: Observations and modeling of GIC in the Chinese largescale highvoltage power networks. J. Space Weather Space Clim., 2014, 4, A03.
All Tables
Locations of geomagnetic observatories in the area of the Chinese Northwest 750 kV power grid.
All Figures
Fig. 1. GIC data at the Ling’ao nuclear power plant on 7–8 (a) and 9–10 (b) November 2004. The horizontal component of the geomagnetic field and its variation rate are also shown based on data from the Zhaoqing Geomagnetic Observatory. 

In the text 
Fig. 2. Layered Earth model for calculating the induced geoelectric field. 

In the text 
Fig. 3. Chinese Northwest 750 kV power grid. Three geomagnetic observatories (GRM, LZH, and JYG) are also shown on the map. (The WMQ observatory is not located in the area of this map.) 

In the text 
Fig. 4. Measured magnetic data and the SECSderived magnetic data on 29–30 May 2005. The horizontal axis is the UT time in hours (a) magnetic data from JYG observatory and the SECSderived magnetic data for Jiuquan substation and (b) magnetic data from LZH observatory and the SECSderived magnetic data for Yongdeng substation. 

In the text 
Fig. 5. Resistivity for the section XiningYinchuan along 750 kV power transmission lines. 

In the text 
Fig. 6. Calculated geoelectric fields at two sites (Jiuquan and Yongdeng) of the Chinese Northwest 750 kV power grid on 29–30 May 2005. The horizontal axis is the UT time in hours (a) EJiuquan and (b) EYongdeng. 

In the text 
Fig. 7. Calculated GIC at two sites (Jiuquan and Yongdeng) of the Chinese Northwest 750 kV power grid on 29–30 May 2005. The horizontal axis is the UT time in hours (a) calculated GIC at Jiuquan substation and (b) calculated GIC at Youngdeng substation. 

In the text 
Fig. 8. Snapshots at 21:35 on 29 May 2005 (a) and at 22:51UT on 29 May 2005 (b) of calculated GIC at different sites of the Chinese Northwest 750 kV power grid. The solid circle represents that the GIC flow into the power network from ground, the hollow one means that GIC flow into the ground. The dashed line with an arrow represents the direction of electric field at that substation (a) calculated GIC results at 21:35UT 29 May 2005 and (b) calculated GIC results at 22:51UT 29 May 2005. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.