機電－電磁暫態混合仿真復合非對稱故障計算方法

楊 洋，肖湘寧，甄曉晨（.新能源電力系統國家重點實驗室（華北電力大學），北京006；.國網石家莊供電公司，河北石家莊050000）

機電－電磁暫態混合仿真復合非對稱故障計算方法

楊 洋1，肖湘寧1，甄曉晨2
（1.新能源電力系統國家重點實驗室（華北電力大學），北京102206；2.國網石家莊供電公司，河北石家莊050000）

為了提高機電－電磁暫態混合仿真對機電側非對稱故障的處理能力，進一步提升混合仿真對各類工況的適應性，本文提出了一種同時考慮機電側非對稱故障和電磁側非對稱故障的計算方法。首先提出了一種機電側發生非對稱故障后的機電側故障仿真方法，該方法根據非對稱故障的類型獲取故障電路導納矩陣，繼而通過修改與故障節點相關的導納元素，將故障電路導納矩陣的疊加到原導納矩陣中獲得故障后的導納矩陣，并通過求解網絡方程獲得故障后的網絡解。然后根據線性電路疊加定理，提出了接口節點戴維南等值電勢的計算方法。為減少電磁側的建模量，采用接口正序等值導納進行戴維南電勢向諾頓等值電流的轉化。最后結合含有一回直流線路的IEEE39節點系統的基于PSCAD＋C架構的機電－電磁暫態混合仿真平臺，驗證了所提方法的精度和有效性。

機電－電磁；混合仿真；非對稱故障；計算方法；精度

0 引言

近年來隨著我國直流工程陸續投運，多直流、多落點、強耦合的交直流輸電格局已經形成［1－4］。傳統LCC?HVDC直流系統具有輸送功率大、造價相對較低的特點，但是由于換流母線電壓降低、電壓不對稱等因素造成的換相失敗也為直流安全穩定運行埋下隱患［5，6］。非對稱故障占電力系統故障的90%以上［7］，在引起直流換相失敗的原因中，逆變側交流系統非對稱故障是其中的主要因素，近年來由其引起的直流換相失敗和短時間內多次換相失敗時有發生［8，9］。如何恰當地對這類交流系統故障引起的直流換相失敗問題進行仿真，是電力系統混合仿真技術的客觀需要。

機電－電磁暫態混合仿真綜合了機電暫態仿真規模化建模以及電磁暫態仿真精細建模的傳統優勢，又避免了各自的準穩態模型不精確以及仿真資源需求過多的內在不足，較好地處理了仿真規模和仿真精度的矛盾，尤其適合大規模交直流電網的仿真分析［10］。傳統的機電－電磁暫態混合仿真的研究往往重點關注于電磁側的仿真，涉及等值形式［11，12］、接口電路寬頻等值［13，14］、相量提取算法及精度［15］等，這類研究的關注點在于電磁側故障之后直流系統或柔性交流輸電系統（FACTS）裝置的動態行為，而對機電側故障以及機電側和電磁側同時發生的復合故障研究鮮見報道。文獻［16］研究了混合仿真機電側仿真方法，提出當機電側故障節點與接口節點的電氣距離小于1.6時，可以認為不用切換頻率相關網絡等值（frequency dependent network equivalent，FDNE）電路，但是沒有給出機電側故障仿真的實現形式，而且所研究的故障也不是非對稱故障。

本文分析了機電側非對稱故障的計算方法，提出了機電側非對稱故障后的接口節點諾頓等值電流的計算流程，并通過含有一回直流的IEEE39節點系統機電－電磁暫態混合仿真驗證了所提方法在處理機電側和電磁側復合非對稱故障方面的精度和合理性。

1 機電側非對稱故障計算

設機電側為一個含有N節點的系統，則存在下列節點導納方程：

其中：Y為（3×N）×（3×N）階方程；Yx為N×N維矩陣；Ux和Ix都為N×1維矩陣。式（1）中，電磁側系統在機電側等值為三序分立注入的電流源，各序電流分別考慮進I0，I1和I2中，接口節點的電壓作為狀態變量參與機電側的計算。

設某時刻在p點發生非對稱故障，則Y矩陣的變化如圖1所示：

當機電側系統沒有發生非對稱故障時，圓點位置的元素均為零，此時各序解耦可以分立求解，由于非對稱故障發生，圓點位置出現了非零元素，各序之間不再解耦，Y矩陣元素的修改量為：

圖1 Y矩陣的變化Fig.1 The change of Y matrix

YF為非對稱故障電路對應的序分量形式表示的故障導納矩陣，不同的非對稱故障對應不同的故障導納形式，三類典型非對稱故障的故障電路如圖2所示。

圖2 三種典型非對稱故障的故障電路Fig.2 Diagrams of fault circuits of three typical asymmetrical faults

圖2 中單相接地故障、兩相接地故障和兩相短路故障的故障導納矩陣的序分量形式分別為：

單相斷線故障和兩相斷線故障分別與ZF＝0時的兩相金屬性接地短路故障和單相金屬性接地故障的故障導納矩陣一致。

經過對故障前導納矩陣的修改，可以求得故障后的網絡節點電壓為：

可知，機電側進行正常情況和故障期間計算的主要差別在于所采用的導納矩陣以及相應的阻抗矩陣不同。正常情況下，三序導納矩陣以及等值阻抗矩陣解耦，因此各序的電壓只與本序對應的電流激勵相關，三序導納陣的求逆可以轉化為3個獨立的序導納矩陣的求逆。當系統發生非對稱故障尤其考慮接地的非對稱故障之后，故障導納矩陣中三序導納陣相互之間不再解耦，經過矩陣求逆之后的等值阻抗陣也不再是3個獨立的矩陣，而一般是3N ×3N（N是節點數目，這里不考慮零序導納陣的規模）的滿陣。此時，某一序的節點電壓不僅與本序的電流激勵相關，還可能與其他兩序的電流激勵相關。在電力系統中，式（6）右側電流相量的非零元素主要包括發電機注入電流、電動機注入電流、負荷的電流和功率部分等，后者主要以正序電流為主，因此故障后某個節點某一序的節點電壓與對應的該節點和上述三類節點之間該序下的互阻抗以及注入的電流相關。從而，式（6）可以相應簡化，當僅考慮發電機注入的正序電流時，式（6）可以表示為：

其中：Zxsf表示故障情況下某序下節點與發電機節點之間的互阻抗矩陣，其維數為N×Ns維，N是系統總的節點數目，Ns是系統內發電機節點數目。

2 接口諾頓等值電流形成和計算

對于接口節點戴維南等值電壓，有計算公式如式（8）所示：

其中：ZⅡ為接口三序等值阻抗矩陣，在機電側沒有發生接地故障時，式（1）的三序獨立，如果電磁側發生非對稱故障而向機電側注入負序或零序電流，可知對于機電側而言，接口的負序和零序電流是唯一的電源激勵，因此UI0＝ZⅡ0II0，UI2＝ZⅡ2II2，Udwn0和Udwn2均為零，戴維南等值電勢轉化為Udwn1＝UI1－ZⅡ1II1，為正序等值。

當機電側發生非對稱故障后，式（8）中的ZⅡ也發生變化，理論上戴維南等值電勢的計算方法為：

其中：ZⅡf為機電側發生故障后接口的序等值阻抗矩陣。按此式計算戴維南等值電勢，則在電磁側的等值電路也要相應進行切換。

由式（7），并考慮接口注入電流，接口電壓的表達式為：

其中：ZISf是接口節點與發電機節點之間的三序互阻抗矩陣，帶入式（9）可得：

也就是，故障后接口戴維南電勢的大小只與發電機注入電流和發電機與接口節點之間的互阻抗矩陣相關，而與接口節點注入電流無關。此時將戴維南電勢轉化為諾頓等值電流可得：

定義為端口－發電機故障情況下的等效互導納陣，其值為端口互阻抗矩陣的逆矩陣與端口－發電機阻抗矩陣的積。

當接口數目較多時，電磁側接口電路的切換將十分復雜。為減少電磁側建模復雜度，無論機電側是否發生非對稱故障，式（9）均表示為：也就是無論機電側是否發生節點故障或者其他接線方式的變化，戴維南電勢的求取始終采用系統正常運行時的各序等值阻抗矩陣。將式（10）帶入式（13）可得：

此時，諾頓等值電流為：

比較式（11）和式（14）可知，當考慮電磁側電路切換時，接口節點的戴維南電勢的表達式僅與發電機注入電流和端口阻抗矩陣以及端口－發電機互阻抗矩陣相關，此時僅僅需要向機電側傳遞接口正序電流即可以完成交互步長的計算和戴維南電勢的求取。當不考慮電磁側電路切換時，戴維南電勢同時取決于發電機注入電流、接口阻抗矩陣、接口－發電機互阻抗矩陣以及接口各序注入電流。因此如果采用式（11）進行戴維南等值電勢的求取，可以不用考慮接口負序和零序注入電流；但如果采用式（14）進行戴維南電勢求取，則需考慮接口注入負序和零序電流。前者傳遞的信息較少，但是電路建模復雜；后者傳遞的信息稍多，但接口電路建模稍簡單，減少了電磁側的建模復雜度，增加仿真方式的靈活性。另一方面，從第二部分仿真結果的比較可知，這種方式帶來的誤差較小，能夠滿足對精度的需求。

進一步將諾頓電流轉化為ABC三相的形式有：

其中：I為N×N維單位矩陣，a＝ej2π／3。

在式（6）的機電暫態計算過程中同時考慮由于機電側故障帶來的節點側節點導納變化以及由于可能的電磁側故障引起的接口注入電流變化，因此這種計算方式可以同時處理機電側和電磁側的各類故障，從而具備了對復合故障的處理能力。考慮機電側故障后的機電側計算流程如圖3所示。

圖3 考慮機電側故障的機電側計算流程Fig.3 Electromechanical side calculation procedure considering faults occurring in the electromechanical side

為了提高仿真效率，在正常和故障計算中用到的阻抗矩陣都在仿真初始化期間計算完成，仿真過程中根據需要直接調用即可。

3 仿真驗證

在PSCAD／EMTDC中建立如圖4所示基于PSCAD＋C架構的機電－電磁暫態仿真模型，將新英格蘭39節點系統中的第39節點的發電機替換為直流輸電線路，直流的逆變側接入39節點，并選擇39母線為接口母線。直流模型參數及控制方式參見CIGRE標準直流模型［17］，發電機采用六階模型，負荷采用恒阻抗模型，直流系統在電磁暫態側仿真，交流系統在機電側系統仿真。機電側仿真步長設定為10 ms，電磁側仿真步長設定為50 μs，交互周期設定為10 ms。對本文提出的處理復合故障的計算方法進行驗證。

3.1處理電磁側非對稱故障能力的驗證

電磁側發生非對稱故障后，以零序、正序和負序等值電流的形式傳遞給機電側進行機電暫態求解。電磁側三序電流的作用對象主要針對式（1）右側的電流列向量，由于此時機電側沒有發生故障，因此機電側采用正常運行時的節點阻抗矩陣進行計算，各序的計算仍然解耦。設t＝5 s時在直流逆變側發生A相金屬性短路故障。直流逆變側功率和36號發電機轉速如圖5所示。

圖5 只考慮電磁側非對稱故障的仿真結果Fig.5 Simulation results only considering electromagnetic side asymmetric fault

圖4 含有一回直流的IEEE39節點系統Fig.4 IEEE 39 system with a DC link

如圖5（a）所示，在故障期間和故障清除后混合仿真得到的直流功率變化與全電磁暫態仿真基本一致；5（b）中發電機轉速在最大值與全電磁暫態仿真基本相同，故障清除后發電機轉速的變化曲線的相位與全電磁有所差異，但轉速偏差量的衰減速度與全電磁仿真基本相似，最終均達到穩態，整體上當電磁側發生非對稱故障后，混合仿真與全電磁暫態仿真均有較為相似的仿真結果，說明當電磁側故障后，所提計算方法具有較好的仿真精度。

比較全電磁仿真與混合仿真計算時間。啟動到穩態、穩態運行和故障的仿真時間分別為1 s，4 s，0.1 s，分別采用全電磁仿真和混合仿真所用的計算時間如表1所示（計算機配置為Intel（R）Core（TM）i3－2310M CPU＠2.10 GHz，內存2 GB（RAM））。

表1 全電磁與混合仿真計算時間對比Table 1 Comparison of calculation time between full electromagnetic and hybrid transient simulation

從表1可見，混合仿真比全電磁仿真從啟動開始更快地進入穩態，在穩態期間具有較高的仿真速度，能更快地處理電磁側的暫態故障，混合仿真相對于全電磁暫態仿真在計算效率方面具有顯著的優勢。

3.2處理機電側非對稱故障能力的驗證

驗證處理機電側非對稱故障的能力。設t＝5 s時在第8節點發生AB相經0.1 Ω電阻相間短路故障，直流逆變側功率和36號發電機轉速如圖6所示。

圖6 只考慮機電側故障的仿真結果Fig.6 Simulation results only considering electromechanical side asymmetric fault

從圖6（a）可見，故障期間直流輸送有功功率接近零，交流故障引起了直流換相失敗，說明混合仿真對由機電側故障引起的直流換相失敗現象的仿真是有效的；從圖6（b）可見，故障期間發電機轉速的最大值與全電磁一致，盡管故障清除后發電機轉速變化曲線擺動的相位存在一定差異，但故障后發電機轉速與額定值差異量的衰減速度基本一致，并且最后均達到穩定。通過在機電側設置AB相相間短路故障，驗證了混合仿真對機電側故障的處理能力。當故障發生在機電側時，故障期間混合仿真的加速比為TM／TH＝5.04／0.38＝13.26，相比故障設置在電磁側的情況（TM／TH＝3.49），混合仿真可以獲得更大的加速比。

3.3同時處理機電側和電磁側同時發生非對稱故障能力的驗證

進一步驗證混合仿真同時處理機電側和電磁側非對稱故障的能力。設t＝5 s時電磁側的直流逆變側交流母線發生金屬性三相接地故障，與此同時機電側發生經0.1 Ω電阻AB相間短路故障，直流逆變側功率和36號發電機轉速如圖7所示。為了進行比較，同時顯示了考慮電磁側等值電路切換的仿真效果。等值電路切換仿真方法見文獻［16］。

圖7 同時考慮機電側和電磁側非對稱故障的仿真結果Fig.7 Simulation results only both considering electromagnetic and electromechanical side asymmetric faults

從圖7（a）可見，故障期間，直流功率接近為零，混合仿真與全電磁仿真具有較好的契合度。故障期間混合仿真對應的發電機轉速略高于全電磁仿真，故障清除后發電機轉速變化曲線的相位與全電磁有所差異，但總體趨勢與全電磁一致，最終二者均趨于穩定。考慮電路切換的仿真結果與未考慮電路切換的仿真結果基本一致，說明在進行機電側故障仿真時可以不考慮機電側故障給電磁側等值電路帶來的變化，而僅通過諾頓等值電流補償電路變化帶來的影響。機電側和電磁側故障的仿真驗證了混合仿真對機電側和電磁側同時故障的復合故障的處理能力。故障同時發生在機電側和電磁時，故障期間混合仿真的加速比為TM／TH＝7.72／1.58＝4.89，介于故障設置在電磁側和機電側之間。

需要指出，當接口數目不唯一，接口電路相對復雜時，采用式（13）和式（14）計算出的諾頓等值電流仍然可以保持相對較高的仿真精度。仍以圖4所示的系統為例，接口位置選擇在2－8節點，包含直流、1、9節點的系統在電磁側仿真，剩余部分在機電側仿真，此時接口數目為2個。設定在t＝5 s時在節點8位置發生AB相經0.1 Ω電阻短路故障，直流功率和36號發電機轉速如圖8所示。

圖8 采用兩個端口且不考慮非對稱故障期間接口電路切換時的仿真結果Fig.8 Simulation results using two ports without interface circuit switching during asymmetric fault

比較圖8和圖4在故障期間的仿真結果可知，2個端口的情況下，故障位置與接口存在一定電氣距離，如果不考慮故障期間的電路切換，相比于接口位置單一且近故障點的情況，其仿真精度略為下降，但整體上與全電磁暫態仿真仍然較為接近，從而進一步說明在故障期間采用式（13）、（14）計算諾頓等值電流的可行性。

4 結語

本文提出了一種處理機電－電磁暫態混合仿真機電側和電磁側非對稱復合故障的計算流程和方法。首先分析了機電側計算非對稱故障的計算過程，然后結合接口節點的電壓解和等值阻抗矩陣提出各序戴維南等值電勢的計算方法，為了減小電磁側接口電路的復雜程度并提高混合仿真對機電側故障的適應性，提出采用穩態等值阻抗計算戴維南電勢和諾頓等值電流的方法。結合含有一回直流的IEEE39節點系統，通過分別在電磁側、機電側以及機電側和電磁側同時設置非對稱故障，驗證了本文所提方法的精度和合理性。

對比了啟動到穩態，穩態運行和故障期間全電磁仿真和混合仿真的計算時間。結果表明，混合仿真在計算效率上比全電磁仿真具有顯著優勢，并且故障設置在機電側比故障設置在電磁側在故障期間具有更大的加速比。

［1］覃 琴，郭 強，周勤勇，等.國網“十三·五”規劃電網面臨的安全穩定問題及對策［J］.中國電力，2015，48（1）：25－32.

QIN Qin，GUO Qiang，ZHOU Qinyong，et al.The security and stability of power grids in 13 th Five?Year planning and counter?measures［J］.Electric Power，2015，48（1）：25－32.

［2］國家電網公司.國家電網公司“十二·五”電網發展規劃［R］.北京：國家電網公司，2012. State Grid Corporation.State grid corporation“Twelf?Five”power grid development plan［R］.Beijign：State Grid Corporation，2012.

［3］劉振亞，張啟平.國家電網發展模式研究［J］.中國電機工程學報，2013，33（7）：1－10.

LIU Zhenya，ZHANG Qiping.Study on the development mode of national power grid of China［J］.Proceedings of the CSEE，2013，33（7）：1－10.

［4］劉振亞.特高壓電網［M］.北京：中國經濟出版社，2005：159－166.

LIU Zhenya.Ultra high voltage grid［M］.Beijing：China Eco?nomic Press，2005：159－166.

［5］朱鑫要，汪惟源，汪成根.江蘇特高壓交直流輸電交互影響研究［J］.江蘇電機工程，2016，35（2）：39－42.

ZHU Xinyao，WANG Weiyuan，WANG Chenggen.Analysis of interaction between UHVAC and UHVDC for Jiangsu Power Grid［J］.Jiangsu Electrical Engineering，2016，35（2）：39－42.

［6］歐開健，任 震，荊 勇.直流輸電系統換相失敗的研究（一）——換相失敗的影響因素分析［J］.電力自動化設備，2003，23（5）：5－8. OU Kaijian，REN Zhen，JING Yong.Research on commutation failure in HVDC transmission system Part 1：Commutation failure factors analysis［J］.Electric Power Automation Equip?ment，2003，23（5）：5－8.

［7］VAAHEDI E，LI W，CHIA T，et al.Large scale probabilistic transient stability assessment using BC Hydro s on?line tool［J］. IEEE Transactions on Power Systems，2000，15（2）：661－667.

［8］李新年，易 俊，李柏青，等.直流輸電系統換相失敗仿真分析及運行情況統計［J］.電網技術，2012，36（6）：266－271.

LI Xinnian，YI Jun，LI Baiqing，et al.Simulation analysis and operation statistics of commutation failure in HVDC transmission system［J］.Power System Technology，2012，36（6）：266－271.

［9］王 晶，梁志峰，江 木，等.多饋入直流同時換相失敗案例分析及仿真計算［J］.電力系統自動化，2015，39（4）：141－146.

WANG Jing，LIANG Zhifeng，JIANG Mu，et al.Case analysis and simulation of commutation failure in multi?infeed HVDC transmission systems［J］.AutomationofElectricPower Systems，2015，39（4）：141－146.

［10］SU Hongtian.Real?time electromagnetic electromechanical hy?brid transient simulation for large power systems［D］.Hong Kong：Hong Kong Polytechnic University，2005.

［11］楊 洋，肖湘寧，陶 順，等.混合仿真電磁側功率源等效誤差原理分析及改進［J］.電力系統自動化，2015，39（24）：104－109.

YANG Yang，XIAO Xiangning，TAO Shun，et al.Electromag?netic side power source equivalent error principle analysis and its improvement for hybrid simulation［J］.Automation of Elec?tric Power Systems，2015，39（24）：104－109.

［12］肖湘寧，楊 洋，陶 順，等.混合仿真電流源等值誤差機理分析及改進［J］.電力建設，2016，37（6）：38－42.

XIAO Xiangning，YANG Yang，TAO Shun，et al.Current sourceequivalenterrormechanismanalysisandits improvement for hybrid simulation［J］.Electric Power Con?struction，2016，37（6）：38－42.

［13］張 怡，吳文傳，張伯明，等.基于頻率相關網絡等值的電磁－機電暫態解耦混合仿真［J］.中國電機工程學報，2012，32（16）：107－114.

ZHANG Yi，WUWenchuan，ZHANGBoming，etal. Frequency dependent network equivalent based electromagnetic and electromechanical decoupled hybrid simulation［J］.Pro?ceedings of the CSEE，2012，32（16）：107－114.

［14］LIN Xi，GOLE A M，YU Ming.A wide?band multi?port system equivalent for real?time digital power system simulators［J］. IEEE Transactions on Power Systems，2009，24（1）：237－249

［15］楊 洋，肖湘寧，陶 順，等.考慮衰減直流分量的dq?120改進算法及其在混合仿真中的應用［J］.電力建設，2016，37（6）：43－48.

YANG Yang，XIAO Xiangning，TAO Shun，et al.An improved dq?120 algorithm considering decaying DC component and Its application in hybrid simulation［J］.Electric Power Construc?tion，2016，37（6）：43－48.

［16］張 怡，吳文傳，張伯明，等.電磁－機電暫態混合仿真中機電側故障的仿真方法［J］.中國電機工程學報，2012，32（19）：81－88.

ZHANG Yi，WU Wenchuan，ZHANG Boming，et al.Simula?tion method of faults on electromechanical side in electromag?netic and electromechanical hybrid simulation［J］.Proceedings of the CSEE，2012，32（19）：81－88.

［17］鄧廣靜，周 威.CIGRE直流輸電標準模型的建模及控制策略仿真研究［J］.江蘇電機工程，2009，28（1）：31－34.

DENG Guangjing，ZHOU Wei.Study on modeling and control strategy of CIGRE HVDC benchmark system［J］.Jiangsu Elec?trical Engineering，2009，28（1）：31－34.

Calculation Method of Complex Asymmetrical Fault
in Electromechanical and Electromagnetic Transient Hybrid Simulation

YANG Yang1，XIAO Xiangning1，ZHEN Xiaochen2
（1.State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources North China Electric Power University（NCEPU），Beijing 102206，China；2.State Grid Shijiazhuang Power Supply Company，Shijiazhuang 050000，China）

To improve the processing ability of asymmetrical eletromechanical side fault as well as the flexibility of electromechanical and electromagnetic transient hybrid simulation，a calculation method of complex asymmetrical fault in the electromechanical and electromagnetic side at the same time is proposed.This method obtains the fault circuit admittance matrix according to the type of asymmetric fault.By modifying the admittance element associated with the faulty node，the fault circuit admittance matrix is superimposed on the normal admittance matrix to obtain the admittance matrix after fault.The network solution is finally achieved by solving the network equation after fault.Based on the linear circuit superposition theorem，the Thevenin equivalent voltage of the interface node is calculated.To reduce the modeling scale of the electromagnetic side，positive sequence equivalent conductance of the interface is used for the conversion of the Thevenin equivalent voltage into Norton equivalent current.At last，through the electromechanical and electromagnetic transient hybrid simulation platform of a IEEE 39 system containing a HVDC line based on PSCAD＋C frame，the simulation accuracy and effectiveness of the method is validated.

Electromechanical and electromagnetic；hybrid simulation；asymmetrical Fault；accuracy

TM73

：A

：2096－3203（2017）02－0014－07

楊 洋

楊 洋（1989—），男，河北邢臺人，博士生，主要研究方向為電力系統仿真及分析；

肖湘寧（1953—），男，湖南澧縣人，教授，博士生導師。主要研究方向為新能源電網中的電力電子技術、電力系統電能質量等；

甄曉晨（1989—），男，河北石家莊人，工程師，從事電力調度控制及電能質量經濟性方面研究工作。

（編輯 劉曉燕）

2016－12－10；

2017－01－28

中央高校基本科研業務費專項資金資助項目（2015XS22）