劉連光，王開讓，魏愷，錢晨，崔曉丹

(1．華北電力大學電氣與電子工程學院，北京市102206;2．國電南瑞科技股份有限公司，南京市211106)

0 引言

與中低緯度地區(qū)相比，靠近南北極的地磁場相對較強，地磁暴的地磁擾動(geomagnetic disturbance，GMD)也相對較強，此外，高緯地區(qū)電離層的極光電集流具有東西向流動特征［1］，因此很多人認為，電網(wǎng)GMD災害是高緯度地區(qū)的問題，且東西走向輸電線路會產(chǎn)生較大的地磁感應電流(geomagnetically induced current，GIC)，該觀點忽略了電網(wǎng)條件及輸電線路參數(shù)的影響。隨著研究的深入，在世界中低緯度國家中，我國首先獲得了地磁暴侵害廣東和江蘇電網(wǎng)的GIC數(shù)據(jù)［2］，發(fā)現(xiàn)了中低緯度地區(qū)存在南北向地電場比東西向地電場大的現(xiàn)象［3-4］，以及電網(wǎng)的直流電阻小是GIC大的原因［5］。這些研究，揭示了中低緯度地區(qū)電網(wǎng)GMD災害的機理，分析了GMD對我國電網(wǎng)安全的危害。中國工程院2013年2月啟動的《我國應對復雜電磁脈沖環(huán)境威脅戰(zhàn)略研究》科技咨詢項目，在對電網(wǎng)和油氣管道地磁暴危害調(diào)研、咨詢和研討基礎上，于2014年12月向國家呈報了《關于加強油氣管網(wǎng)和電網(wǎng)地磁暴災害防御的建議》報告，建議開展電網(wǎng)災害防御和管道影響測評研究［6］。

防御電網(wǎng)GMD災害，可考慮從電網(wǎng)規(guī)劃設計、變壓器設計制造和電網(wǎng)調(diào)度運行等方面采取措施［7-9］。其中，在電網(wǎng)規(guī)劃設計上，主要是合理規(guī)劃和科學安排接入變電站的饋電線路的方向及數(shù)量［5］，使流入GMD災害風險高的變電站變壓器的GIC最小。因此，準確計算電網(wǎng)輸電線路以及流入變電站變壓器的GIC是防御電網(wǎng)GMD災害的前提。目前，在電網(wǎng)GIC計算上，主要是根據(jù)可獲得的地磁臺的GMD觀測數(shù)據(jù)計算電網(wǎng)的GIC。由于電網(wǎng)的GIC與GMD強度及變化率、深層大地的電性構造、輸電線路的結構及參數(shù)以及變電站饋電線路的接入方向及數(shù)量有關［2-5，10］，因此，根據(jù)地磁臺 GMD 數(shù)據(jù)計算電網(wǎng)GIC，首先要建立GMD的感應地電場模型，而感應地電場模型及計算涉及深層大地的電性構造，由于缺乏深層大地的電導率精細數(shù)據(jù)，基于地磁數(shù)據(jù)計算電網(wǎng)的GIC困難和誤差較大。

東半球空間環(huán)境地基綜合監(jiān)測子午鏈(簡稱子午工程)是我國空間科學領域的第1個國家重大科技基礎設施工程［11-12］。子午工程于2008年1月5日開工建設，2012年10月23日通過國家驗收。子午工程可提供地磁暴的地電場實測數(shù)據(jù)，本文利用廣東肇慶臺站的地電場數(shù)據(jù)，直接計算廣東500 kV電網(wǎng)的GIC，與電網(wǎng)GIC實測數(shù)據(jù)以及基于地磁數(shù)據(jù)的GIC計算模型及算法進行比較，研究電網(wǎng)GIC算法。

1 電網(wǎng)GIC全節(jié)點模型及算法

用地電場數(shù)據(jù)計算電網(wǎng)GIC，可不考慮空間條件和大地條件的影響，直接將地電場實測數(shù)據(jù)與電網(wǎng)條件相結合建立電網(wǎng)GIC模型，利用掌握的方法對線路和變壓器等元件等效處理。

1.1 廣東電網(wǎng)全節(jié)點GIC模型構建

2012年，美國電力科學研究院(EPRI)牽頭，聯(lián)合芬蘭氣象科學研究院等單位，共同開發(fā)了計算電網(wǎng)GIC 的“GIC-Benchmark”標準算例［13］，文獻［14］針對“GIC-Benchmark”算例，提出了計算多電壓等級電網(wǎng)GIC的全節(jié)點模型，以及輸電線路和變壓器等元件模型及參數(shù)的處理方法，據(jù)此可建立廣東電網(wǎng)500 kV及220 kV電網(wǎng)的全節(jié)點GIC模型。由于本文主要通過與廣東國安500 kV變電站GIC實測數(shù)據(jù)進行比較，來驗證直接計算電網(wǎng)GIC的精度和可行性，因此廣東電網(wǎng)GIC的計算采用全節(jié)點模型，結果取500 kV變電站 GIC計算數(shù)據(jù)，并與實測 GIC進行比較。

在建立電網(wǎng)全節(jié)點GIC模型時，網(wǎng)內(nèi)各變電站母線和中性點均作為獨立節(jié)點，輸電線路和變壓器各級繞組作為獨立支路。輸電線路支路等效模型為電源和電阻的串聯(lián)電路，其中，電阻為輸電線路直流電阻，電源為地電場在輸電線路上的等效電壓源，輸電線路的等效電路如圖1所示。變壓器繞組的支路等效模型為繞組直流電阻，因此本文中的變壓器等效電路如圖2所示，其中，公共繞組和串聯(lián)繞組分別作為一條獨立支路。根據(jù)廣東電網(wǎng)的構成與結構，對所有線路和變壓器等元件進行等效處理，可建立電網(wǎng)的全節(jié)點GIC模型。

圖1 輸電線路等效電路Fig.1 Equivalent circuit of transmission line

圖2 自耦變壓器等效電路Fig.2 Equivalent circuit of autotransformer

1.2 基于節(jié)點導納矩陣的GIC算法

本文采用在電網(wǎng)工程中得到廣泛應用的節(jié)點導納矩陣法計算電網(wǎng)GIC，并忽略了變壓器繞組中電感的影響。在節(jié)點導納矩陣法中，支路模型如圖3所示，其中i、k表示網(wǎng)絡中任意2個節(jié)點。與節(jié)點i相連的支路分為2種，與全節(jié)點模型中的支路模型相對應:(1)輸電線路，如i-m支路;(2)變壓器支路，如i-p支路，其中輸電線路上的電流源和電阻并聯(lián)支路是通過全節(jié)點模型中的電壓源和電阻串聯(lián)支路經(jīng)戴維南變換得到的。yi=0時表示節(jié)點i不接地。

圖3 節(jié)點導納計算電路Fig.3 Node admittance calculation circuit

對任意節(jié)點i，根據(jù)基爾霍夫電流定律有:

式中:Ii為由節(jié)點i注入大地的電流;N為節(jié)點數(shù);Iki為由節(jié)點k流向節(jié)點i的電流。顯然，對于不接地的節(jié)點，Ii=0。

對于支路電流Iki有:

式中:Uk和Ui分別為節(jié)點k和節(jié)點i的節(jié)點電壓;Jki和yik分別為節(jié)點i和節(jié)點k之間的支路電流源和線路直流電阻。

令Ji=，并將式(2)代入式(1)，則對于節(jié)點i有:

在式(3)中，Ji為與節(jié)點i相連的支路電流源之和，式中未知量為節(jié)點電壓和入地電流，對于節(jié)點i，其入地電流和節(jié)點電壓滿足:

將式(4)代入式(3)，式中的未知量僅剩節(jié)點電壓，即

寫成矩陣形式:

式中:Y為節(jié)點導納矩陣;J為電流源列向量;U為節(jié)點電壓列向量。

因此，電網(wǎng)全節(jié)點模型中的節(jié)點電壓由節(jié)點導納矩陣的逆矩陣和電流源列向量相乘求得:

將式(7)所得節(jié)點電壓U中i節(jié)點的電壓代入公式(4)，即可求得入地電流，即GIC。流向大地為正值，流入變壓器中性點為負值。

2 基于地電場數(shù)據(jù)的電網(wǎng)GIC計算

由于不需用地磁數(shù)據(jù)計算感應地電場，在建立了電網(wǎng)GIC模型的基礎上，電網(wǎng)GIC的計算也變得相對簡單。除規(guī)劃設計防災外，考慮到調(diào)度運行適時監(jiān)測防災需要，將利用地電場監(jiān)測數(shù)據(jù)直接計算電網(wǎng)GIC的方法分為2部分，首先確定地磁暴GMD的強度和起止時間，進而收集響應電網(wǎng)覆蓋區(qū)域內(nèi)的地電場監(jiān)測數(shù)據(jù);然后基于電網(wǎng)構建的電網(wǎng)全節(jié)點GIC模型，綜合電網(wǎng)構成條件和地電場因素對電網(wǎng)GIC進行計算，計算流程圖如圖4所示，計算方法描述如下。

(1)確定強磁暴GMD發(fā)生的起止時間，并獲取該時段的地電場監(jiān)測數(shù)據(jù)，本文取廣東肇慶監(jiān)測點該時段的時變地電場數(shù)據(jù)。

(2)根據(jù)構建的電網(wǎng)GIC模型，輸入步驟(1)中獲取的t時刻地電場數(shù)據(jù)，計算各節(jié)點間的等效電壓源，本文算例為廣東電網(wǎng)模型。

(3)結合步驟(2)中計算所得的各節(jié)點間的等效電壓源，計算t時刻電網(wǎng)中各個節(jié)點的GIC，并將數(shù)據(jù)保存。

圖4 基于全節(jié)點模型的電網(wǎng)GIC計算流程Fig.4 Power grid GIC calculation process based on full-node GIC model

(4)檢查t是否大于時間T，如小于T則讀取下一時刻的地電場數(shù)據(jù)循環(huán)計算GIC，并保存;如果大于T，則調(diào)用GIC數(shù)據(jù)保存模塊。

(5)將計算得到的所有時刻的GIC數(shù)據(jù)按時序進行整理，輸出時變GIC計算結果。

3 算例分析

3.1 計算數(shù)據(jù)

為驗證中、小地磁暴的GMD電網(wǎng)效應，本文算例取北京時間2014年9月12日的地磁暴事件來計算廣東電網(wǎng)的GIC。該次事件受9日和11日的日冕物質(zhì)拋射(CME)事件的影響，分別在12日出現(xiàn)Kp=5的小地磁暴，13日出現(xiàn)Kp=7的大地磁暴，Dst指數(shù)最小值為－75 nT。依據(jù)此次的地磁暴事件的起止時間，收集了廣東電網(wǎng)國安(N22E113)500 kV變電站監(jiān)測的電網(wǎng)GIC數(shù)據(jù)以及子午工程廣東肇慶監(jiān)測站(N23E112)的地電場的監(jiān)測數(shù)據(jù)，以驗證所建電網(wǎng)模型及GIC算法的有效性，監(jiān)測數(shù)據(jù)的持續(xù)時間共計33 h。該次地磁暴Dst指數(shù)隨時間的擾動情況如圖5所示。

圖5 Dst指數(shù)擾動情況Fig.5 Disturbance of Dstindex

根據(jù)掌握的廣東電網(wǎng)數(shù)據(jù)和資料［3-4］，收集最新增加的電網(wǎng)數(shù)據(jù)和資料，廣東500 kV電網(wǎng)接線如圖6所示，該電網(wǎng)共54個變電站，62條線路，由本文第一節(jié)中各支路等效方法，得到廣東電網(wǎng)全節(jié)點模型中線路和變壓器支路的等效參數(shù)，其中國安站的變電站等效電阻為1.6 Ω，變電站等效電阻包含變壓器支路電阻和接地電阻。與國安站相連的珠海站和桂山站，線路回路數(shù)均為2回，線路等效電阻分別是0.22 Ω和0.52 Ω，廣東電網(wǎng)各站點和線路的具體參數(shù)見文獻［15］，進而采用導納矩陣法對廣東電網(wǎng)中可能存在的GIC進行計算。

3.2 結果分析

為比較地電場監(jiān)測數(shù)據(jù)、電網(wǎng)GIC數(shù)據(jù)和計算數(shù)據(jù)的形態(tài)及特征，首先繪制地電場監(jiān)測數(shù)據(jù)的北向(x)和東向(y)分量波形，如圖7所示。

圖6 廣東500 kV電網(wǎng)接線示意圖Fig.6 Guangdong 500 kV power grid

圖7 廣東肇慶實測地電場的x和y分量Fig.7Measured x and y of geoelectric field in Guangdong Zhaoqing

根據(jù)圖7數(shù)據(jù)及本文模型和算法，可對廣東電網(wǎng)GIC進行仿真計算。在計算中，由于廣東電網(wǎng)覆蓋區(qū)域的地理跨度不大，假設地電場在電網(wǎng)覆蓋區(qū)域內(nèi)的擾動情況是一致的，區(qū)域內(nèi)均采用肇慶監(jiān)測站的地電場數(shù)據(jù)，據(jù)此計算圖6電網(wǎng)中54個變電站節(jié)點和62條線路的 GIC。其中，500 kV國安站GIC計算值(I'GIC)與實測值(IGIC)的對比結果如圖8所示。

此次地磁暴的強度相對較小，實測的國安站節(jié)點GIC最大值為2.35 A，最小值為－5.56 A，計算所得的GIC最大值為2.26 A，最小值為－3.70 A。分析認為，計算誤差為肇慶監(jiān)測點與國安站實際地電場不同所致。但從圖8可看出，計算值與實測值在GIC形態(tài)、突變時間上，與圖7地電場形態(tài)、突變時間相吻合，并且比用地磁數(shù)據(jù)算法的誤差?。?-4］，表明地電場算法更能準確反應GIC的形態(tài)和突變時間特征。

圖8 廣東國安變電站實測與仿真GIC比較Fig.8 Comparison of measured(IGIC)and calculated(I＇GIC)GIC at Guangdong Guo'an substation

4 結論

(1)利用子午工程提供的地電場數(shù)據(jù)計算電網(wǎng)GIC的方法，比利用地磁場數(shù)據(jù)計算電網(wǎng)GIC的方法簡單，且更能準確反應GIC的形態(tài)和突變時間等特征。由于方法簡單以及識別GIC時變特征的精度高，利用地電場數(shù)據(jù)計算GIC的方法更適合在運行調(diào)度中應用，通過地電場監(jiān)測數(shù)據(jù)或GIC計算數(shù)據(jù)識別GIC引發(fā)的故障，并據(jù)此制定防災策略。

(2)由于肇慶地電場監(jiān)測站與國安變電站的位置不同，以及大地深層電阻率不同等原因，導致國安變電站GIC計算值存在一定的誤差，證明了由于大地深層電性構造的不同，不同位置的GMD的感應地電場不相同，準確計算GIC，需要考慮大地電性構造以及電阻率橫向變化的影響。

(3)子午工程提供的地電場監(jiān)測數(shù)據(jù)，為計算電網(wǎng)的GIC提供了一種簡單方法，但子午工程地電場監(jiān)測臺站的數(shù)量有限，在臺站的布局上，還不能適應我國大規(guī)模電網(wǎng)GMD安全分析及防御地磁暴災害的需要，建議加強對地磁暴感應地電場的監(jiān)測力度，為防御地磁暴災害提供數(shù)據(jù)及服務。

(4)利用地電場監(jiān)測數(shù)據(jù)計算電網(wǎng)GIC是地磁暴電網(wǎng)效應研究的新進展，加強地電場監(jiān)測數(shù)據(jù)深化應用和直接計算電網(wǎng)GIC的方法，不僅對電網(wǎng)規(guī)劃設計和運行調(diào)度具有重要意義，直接計算GIC的方法還可用于變壓器等電力設備的GIC衍生效應的研究，這是下一步的研究任務。

致 謝

本項成果使用國家重大科技基礎設施子午工程科學數(shù)據(jù)，感謝中國南方電網(wǎng)公司對電網(wǎng)GIC監(jiān)測工作的支持。

［1］ Bolduc L，Langlois P，Boteler D，et al．A study of geomagnetic disturbance in Quebec 1．genral results［J］．IEEE Transactions on Power Delivery，1998，13(4):1251–1256．

［2］劉連光，劉春明，張冰，等．中國廣東電網(wǎng)的幾次強磁暴影響事件［J］．地球物理學報，2008，51(4):976-981．Liu lianguang，Liu Chunming，Zhang Bing，et al．Strong magneticstorm sinfluenceon Chinas Guangdong power grid［J］．Chinese Journal of Geophysics，2008，51(4):976-981．

［3］Liu C M，Liu L G，Pirjola R，et al．Calculation of geomagnetically induced currents in mid-to low-latitude power grids based on the plane wave method:A preliminary case study［J］．Space Weather，2009，7(4):1291-1298

［4］ Liu C M，Liu L G，Pirjola R，et al．Geomagnetically induced currents in the high-voltage power grid in China［J］． IEEE Transactions on Power Delivery，2009，24(4):2368-2374．

［5］劉連光，郭世曉，魏愷，等．基于全節(jié)點GIC模型的三華電網(wǎng)地磁感應電流計算［J］．電網(wǎng)技術，2014，38(7):1946-1952．Liu Lianguang， Guo Shixiao， Wei Kai， et al．Calculation of geomagnetically induced currents in interconnected North China Central China-East China power grid based on full-node GIC model［J］．Power System Technology，2014，38(7):1946-1952．

［6］中國工程院．關于加強油氣管網(wǎng)和電網(wǎng)地磁暴災害防御的建議［R］．北京:中國工程院，2014．

［7］ Zhang B，Liu L G，Liu Y L，et al．Effect of geomagnetically induced current on the loss of transformer tank［J］．IET Electric Power Applications，2010，4(5):373-379．

［8］張冰，劉連光，肖湘寧．地磁感應電流對變壓器振動噪聲的影響［J］．高電壓技術，2009，35(4):900-904．Zhang Bing， Liu Lianguang， Xiao Xiangning． Effects of geomagnetically induced current on the vibration and noise of transformer［J］．High Voltage Engineering，2009，35(4):900-904

［9］劉連光．災害空間天氣對我國電網(wǎng)安全的影響及風險［J］．中國工程科學，2010，12(9):29-33．Liu Lianguang．Influence and hazard of disastrous space weather on power grid in China［J］．Engineering Sciences，2010，12(9):29-33．

［10］董博，王澤忠，劉連光，等．大地電導率橫向突變處磁暴感應地電場的鄰近效應［J］．地球物理學報，2015，58(1):238-246．Dong Bo，Wang Zezhong，Liu Lianguang，et al．The proximity effect on the induced geoelectric field at the interface of different conductivity structures with lateral variations during geomagnetic storms［J］．Chinese Journal of Geophysics，2015，58(1):238-246．

［11］Wang C．New chains of space weather monitoring stations in China［J］．Space Weather，2010，8(8):5．

［12］王赤，馮學尚，萬衛(wèi)星，等．東半球空間環(huán)境地基綜合監(jiān)測子午鏈簡介［J］．國際地震動態(tài)，2009(6):32-38．Wang Chi，F(xiàn)eng Xueshang，Wan Weixing，et al．The meridian chain of comprehensive gound-based space environment monitors in the eastern hemisphere［J］． Recent Developments in World Seismology，2009(6):32-38．

［13］ Horton R，Boteler D，Overbye T J，et al．A test case for the calculation ofgeomagnetically induced currents［J］．IEEE Transactions on Power Delivery，2012，27(4):2368-2373．

［14］鄭寬，劉連光，Boteler D H，等．多電壓等級電網(wǎng)的 GICBenchmark建模方法［J］．中國電機工程學報，2013，33(16):179-186．Zheng Kuan，Liu Lianguang，Boteler D H，et al．Modelling geomagnetically induced currents in multiple voltage levels of a power system illustrated usingthe GIC-benchmark case［J］．Proceedings of the CSEE，2013，33(16):179-186．

［15］李云龍．中低緯電網(wǎng)地磁感應電流影響因素研究［D］．北京:華北電力大學，2014．Li Yunlong．Research on the effect factors of geomagnetic induced current in mid-low latitude power grid［D］．Beijing:North China Electric Power University，2014．