高淑萍，李元澤，宋國兵，呂宇星，左俊杰，沈渠旺

（1. 西安科技大學電氣與控制工程學院，西安 710600；2. 西安交通大學電氣工程學院，西安 710049）

0 引言

目前輸電網(wǎng)絡(luò)正快速發(fā)展，高壓直流輸電技術(shù)也隨之不斷進步。傳統(tǒng)高壓直流輸電是基于電網(wǎng)換向換流器高壓直流輸電（line-commutated converter based high volatage direct current，LCC-HVDC）技術(shù)，這種直流輸電技術(shù)成熟，建設(shè)運行成本低。隨著高壓直流輸電對換流器的要求越來越高，便誕生了一種基于模塊化多電平換流器高壓直流輸電（modular multilevel converter based high voltage direct current，MMC-HVDC）技術(shù)。這種直流輸電功率控制靈活，不存在換相失敗的問題。為了將這兩種直流輸電的優(yōu)勢結(jié)合起來，近年來出現(xiàn)了混合直流輸電系統(tǒng)。但在目前世界上的混合直流輸電系統(tǒng)實際工程中，除了廣泛應用的混合雙端和混合多端直流輸電系統(tǒng)，僅在挪威和丹麥之間存在一條混合四極直流輸電系統(tǒng)，還沒有混合雙極直流輸電系統(tǒng)投入運營，針對混合雙極直流輸電系統(tǒng)的研究較少?；旌纤臉O直流輸電系統(tǒng)占地面積大，經(jīng)濟性較差，混合雙極直流輸電系統(tǒng)相較于混合四級直流輸電系統(tǒng)控制方式簡單，造價低，占地面積小，具有很好的應用前景［1-2］。LCC-MMC 型混合雙極直流輸電系統(tǒng)正極由電網(wǎng)換相換流器（LCC）構(gòu)成，負極由模塊化多電平換流器（MMC）構(gòu)成。由于其兼有兩者的優(yōu)點，所以混合直流輸電比傳統(tǒng)直流輸電和柔性直流輸電具有更明顯的優(yōu)勢，它的誕生將改善未來直流輸電工程的傳輸特性，其應用領(lǐng)域會越來越大［3-4］。

反向行波通常被用來作為行波保護的判斷依據(jù)，當線路發(fā)生故障時會產(chǎn)生以故障點為中心向線路兩側(cè)傳播的行波，通過檢測反向行波即可判定直流線路是否發(fā)生區(qū)內(nèi)故障［5-7］。目前大多數(shù)輸電線路將行波保護和微分欠壓保護作為高壓直流輸電線路保護的主保護，電流差動保護和低電壓保護作為后備保護［8-11］。SIMENS 和ABB 制定的行波保護作為主流保護方案被廣泛使用，但兩種保護在實際高壓直流輸電工程的運行經(jīng)驗表明其在耐過渡電阻的能力上均存在一定欠缺［12-15］。

文獻［16］在分析行波傳輸特性的基礎(chǔ)上基于高頻阻抗在正反方向故障時的不同提出了一種混合直流輸電縱聯(lián)保護原理。文獻［17］根據(jù)直流線路區(qū)內(nèi)外故障時線路兩端電流有無線性關(guān)系基于相關(guān)系數(shù)提出了一種混合直流輸電線路縱聯(lián)保護方法。文獻［18］在傳統(tǒng)行波保護基礎(chǔ)上加以小步長采樣提出了一種可提升保護靈敏性的新型直流線路保護方案。文獻［19］通過比較區(qū)內(nèi)外故障時線路兩端測量波阻抗與測點背側(cè)波阻抗和線路波阻抗的差異提出一種基于波阻抗的縱聯(lián)保護原理。文獻［20］根據(jù)區(qū)內(nèi)外故障時故障前、反行波到達保護點時刻的不同構(gòu)成了VSC-HVDC 縱聯(lián)方向保護判據(jù)。文獻［21］利用直流線路區(qū)內(nèi)外故障時的電壓特征信號的暫態(tài)高頻分量差異，通過小波包變換得到的低頻與高頻能量和的比值構(gòu)成保護判據(jù)。文獻［22］利用小波包變換，基于保護安裝處獲得的暫態(tài)能量差異，提出了一種單端電氣量的行波保護新方法。文獻［23］通過對故障發(fā)生后的行波波前信息以及波前形狀的分析和過波前相關(guān)系數(shù)的計算提出了一種直流輸電線路超高速保護原理。

混合雙極直流輸電系統(tǒng)與普通直流輸電系統(tǒng)相比在輸電線路上的波過程一致，在行波方向的判別上沒有差異。其差異主要體現(xiàn)在：由于LCC 和MMC換流站的拓撲結(jié)構(gòu)和換流站元件參數(shù)不同，導致其邊界對波過程的響應特性存在差異，即行波波頭到達換流站后的折反射是不一樣的，所以電壓行波和電流行波幅值的整定應該是與邊界特性相關(guān)的。

本文在分析直流輸電系統(tǒng)中行波傳播特性的基礎(chǔ)上提出了一種新的利用暫態(tài)電流行波突變的縱聯(lián)方向保護原理，該原理包括故障方向判別算法和故障極點識別算法。對所提方法在PSCAD/MATLAB中進行仿真實驗，大量的仿真結(jié)果表明此保護原理不受故障類型以及故障發(fā)生距離的影響，并且具有很好的耐過渡電阻能力。

1 LCC-MMC 混合雙極直流輸電系統(tǒng)結(jié)構(gòu)分析

圖1 所示為一個±500 kV 的LCC-MMC 混合雙極直流輸電系統(tǒng)結(jié)構(gòu)框圖。不同于LCC-MMC 混合雙端直流輸電系統(tǒng)的整流側(cè)為LCC型換流器，逆變側(cè)為MMC 型換流器，LCC-MMC 混合雙極直流系統(tǒng)的正極采用LCC 型換流器，負極采用MMC 型換流器。圖1 中a、b、c、d 分別為區(qū)內(nèi)線路保護測量安裝位置，L為平波電抗器，以平波電抗器為邊界，將直流輸電系統(tǒng)劃分為區(qū)內(nèi)部分和區(qū)外部分，其中直流線路為區(qū)內(nèi)部分，其余的為區(qū)外部分。f1—f7分別為不同的故障，其中f1、f3、f4、f7為區(qū)外故障，f2、f5為區(qū)內(nèi)單極故障，f6為區(qū)內(nèi)雙極故障。輸電系統(tǒng)各部分具體參數(shù)如表1所示。

圖1 LCC-MMC混合雙極直流輸電系統(tǒng)結(jié)構(gòu)框圖Fig. 1 Structure block diagram of LCC-MMC hybrid bipolar DC transmission system

表1 LCC-MMC混合雙極直流輸電系統(tǒng)各部分參數(shù)Tab. 1 Parameters for each part of LCC-MMC hybrid bipolar DC transmission system

2 行波傳輸原理分析

2.1 行波傳輸過程分析

行波保護原理是通過故障時故障點產(chǎn)生的行波來檢測故障。當線路上發(fā)生故障時故障點處將產(chǎn)生行波，該行波會以故障點為中心向線路兩側(cè)傳播。一般情況下故障信息反映在電壓、電流或兩者的組合中，因此在極端情況下這些參數(shù)可以用來檢測故障［24］。

式中：u2、i2分別為行波傳播線路末端電壓和電流；v為行波波速；Z為輸電線路波阻抗，其值為其中為L0、C0分別為感抗和容抗，根據(jù)表1參數(shù)計算可得線路波阻抗Z=319.556 Ω；γ為線路傳播衰減常數(shù)，其值由線路的阻抗導納決定；x為行波傳輸距離。

2.2 行波的折射與反射

圖2 為行波在阻抗不連續(xù)線路上的傳輸示意圖，A為阻抗不連續(xù)點，l1、l2分別為兩條阻抗不同的線路，其阻抗值分別為Z1和Z2，ua、ia分別為電壓和電流入射波，ub、ib分別為入射波通過阻抗不連續(xù)點后的電壓和電流折射波，uc、ic分別為電壓和電流反射波。

圖2 行波傳輸示意圖Fig.2 Schematic diagram of traveling wave transmission

故障發(fā)生后故障點將產(chǎn)生故障行波，該行波會以故障點為中心沿著線路向兩側(cè)傳輸，隨后在行波到達阻抗不連續(xù)點（如故障點或電容電感處）時發(fā)生行波折射和反射現(xiàn)象。

行波傳輸反射系數(shù)和折射系數(shù)分別為：

式中λ、β分別為行波傳輸?shù)姆瓷湎禂?shù)與折射系數(shù)， -1≤λ≤1，β≥0。

3 故障特征分析

圖3—5 分別為直流輸電系統(tǒng)在區(qū)內(nèi)故障和區(qū)外故障下的行波傳播過程。M、N 分別為整流站、逆變站。M1、M2、M3、M4 為整流側(cè)不同的行波突變點。N1、N2、N3、N4為逆變側(cè)不同的行波突變點。

圖3 直流輸電系統(tǒng)區(qū)內(nèi)故障時行波傳播過程Fig.3 Traveling wave propagation process of DC transmission system when an internal fault occurs

表2—4 分別給出了直流輸電系統(tǒng)在區(qū)內(nèi)故障和區(qū)外故障時整流側(cè)和逆變側(cè)的電流行波。λM、λf和λN分別為電流行波在不同位置（整流側(cè)、故障側(cè)和逆變側(cè)）的反射系數(shù)，βM、βf、βN分別為電流行波在不同位置（整流側(cè)、故障側(cè)和逆變側(cè)）的折射系數(shù)，is為初始電流行波。

表2 區(qū)內(nèi)故障時兩端電流行波Tab. 2 Current traveling wave at both ends when an internal fault occurs

由表1 參數(shù)并結(jié)合式（6）計算可得，正極線路整流側(cè)和逆變側(cè)反射系數(shù)λM+=λN+=-0.961，正極線路整流側(cè)和逆變側(cè)折射系數(shù)βM+=βN+=1.961，負極線路整流側(cè)反射系數(shù)λM-=-0.491，負極線路逆變側(cè)反射系數(shù)λN-=-0.846，負極線路整流側(cè)折射系數(shù)βM-=1.491，負極線路逆變側(cè)折射系數(shù)βN-=1.846。設(shè)線路參考方向為母線指向線路，由圖3—5 可知，is<0。

為滿足繼電保護速動性要求，取區(qū)內(nèi)故障時暫態(tài)電流行波第一次到達保護安裝處發(fā)生折反射為標準（即圖3 和表2 中M1、N1 時）。當發(fā)生區(qū)內(nèi)正極故障時，此時整流側(cè)電流iM=-1.961is，逆變側(cè)電流iN=-1.961is，因為is<0，故iM>0，iN>0。同理，當發(fā)生區(qū)內(nèi)負極故障時，iM=-1.491is>0，iN=-1.846is>0。

同理，取區(qū)外整流側(cè)故障時暫態(tài)電流行波第一次到達保護安裝處發(fā)生折反射為標準（即圖4 和表3中M1、N1 時），此時iM=1.961is，iN=-3.846is。因為is<0，故iM<0，iN>0。同理，當發(fā)生區(qū)外整流側(cè)負極故障時亦有此結(jié)論，在此不再贅述。

圖4 直流輸電系統(tǒng)區(qū)外整流側(cè)故障時行波傳播過程Fig.4 Traveling wave propagation process of DC transmission system at rectifier side when an external fault occurs

表3 區(qū)外整流側(cè)故障時兩端電流行波Tab. 3 Current traveling waves at both ends when an external fault occurs at rectifier side

同理，取區(qū)外逆變側(cè)故障時暫態(tài)電流行波第一次到達保護安裝處發(fā)生折反射為標準（即圖5 和表4中M1、N1 時），此時iM=-3.846is，iN=1.961is。因為is<0，故iM>0，iN<0。同理，當發(fā)生區(qū)外逆變側(cè)負極故障時亦有此結(jié)論，在此不再贅述。

圖5 直流輸電系統(tǒng)區(qū)外逆變側(cè)故障時行波傳播過程Fig.5 Traveling wave propagation process of DC transmission system at inverter side when an external fault occurs

表4 區(qū)外逆變側(cè)故障時兩端電流行波Tab. 4 Current traveling wave at both ends at inverter side when an external fault occurs

4 基于暫態(tài)電流行波突變的縱聯(lián)方向保護原理

4.1 保護啟動判據(jù)

啟動功能的主要要求是可靠性，因為保護算法必須針對每一個可能的故障啟動。通過使用一種可以區(qū)分正常運行和故障的電流判據(jù)可以快速檢測直流線路故障，該判據(jù)可以表示為：

式中：di（n）和di（n+1）分別為采樣點編號n、n+1的差分電流行波；M為累計次數(shù)（經(jīng)驗值為3）；Δset為基于保護標準的門檻值。

保護啟動判據(jù)門檻值設(shè)定一般需要保證在最不靈敏情況下的故障時保護也可以可靠啟動，此處設(shè)定為2 A。當式（4）中的條件滿足時保護被激活。

4.2 保護方向識別判據(jù)

由圖3 可知，當發(fā)生區(qū)內(nèi)故障時整流側(cè)與逆變側(cè)均為正方向故障。由表2 分析可知，整流側(cè)和逆變側(cè)的暫態(tài)電流行波發(fā)生正向突變。

由圖4 可知，當故障發(fā)生為整流側(cè)區(qū)外故障時，此時對于整流側(cè)為反方向故障，對于逆變側(cè)為正方向故障。由表3 分析可知，整流側(cè)的暫態(tài)電流行波發(fā)生負向突變，逆變側(cè)的暫態(tài)電流行波發(fā)生正向突變。

由圖5 可知，當逆變側(cè)發(fā)生區(qū)外故障時對于逆變側(cè)為反方向故障，對于整流側(cè)則為正方向故障。由表4 分析可知，逆變側(cè)的暫態(tài)電流行波發(fā)生負向突變，整流側(cè)的暫態(tài)電流行波發(fā)生正向突變。即當電流行波i>0時為正向故障，當電流行波i<0時為反向故障。

由此可構(gòu)成縱聯(lián)方向保護判據(jù)，當故障發(fā)生后，整流側(cè)與逆變側(cè)檢測到的暫態(tài)電流行波突變方向相同時則判斷為區(qū)內(nèi)故障，保護動作；整流側(cè)與逆變側(cè)檢測到的暫態(tài)電流行波突變方向相反時則判斷為區(qū)外故障，保護不動作。

4.3 故障極識別算法

當發(fā)生區(qū)內(nèi)正極故障時正極作為故障極其電流行波突變幅值大于負極電流突變幅值，其值大于1；當發(fā)生區(qū)內(nèi)負極故障時，負極作為故障極其電流行波突變幅值大于負極電流突變幅值，其值遠小于1；當發(fā)生區(qū)內(nèi)雙極故障時，正負極電流行波均有一定突變。因此，可得到故障極識別公式如式（5）所示。

式中：β為故障極識別參數(shù)；|iPmax|和|iNmax|分別為正極和負極暫態(tài)電流行波的最大值。

選極判據(jù)的門檻值設(shè)定一般需要考慮故障距離以及所設(shè)過渡電阻最大情況下的影響，并留有一定的裕度，故本文β的門檻值取0.6和1。

則故障極點識別算法設(shè)置為：

4.4 影響因素分析

4.4.1 過渡電阻影響

過渡電阻的存在會削弱初始行波的幅值，不同過渡電阻對初始行波幅值的削弱程度不一樣，且過渡電阻值越高則初始行波幅值越低，所以當初始行波到達邊界發(fā)生折反射后保護安裝處的電壓電流行波幅值也會降低。

4.4.2 故障距離影響

故障距離對行波的影響體現(xiàn)在行波波頭的上升沿陡度值，故障距離不同則上升沿陡度不同。故障距離保護安裝處越遠，則初始行波在線路上傳播的距離越遠，保護安裝處檢測的電壓電流行波幅值的上升沿陡度越平緩。

綜上所述，過渡電阻和故障距離僅僅對行波幅值和上升沿陡度產(chǎn)生影響，而對行波的突變方向不會產(chǎn)生影響。

4.4.3 保護流程圖

圖6 給出了基于暫態(tài)電流行波突變的保護流程框圖。

圖6 保護流程圖Fig.6 Flow chart of the protection

5 仿真驗證

如圖1 所示，通過PSCAD/EMTDC 構(gòu)建一個500 kV 的LCC-MMC 混合直流輸電系統(tǒng)。此系統(tǒng)采用頻域相關(guān)模型，輸電線路采用架空線，線路長度為1 500 km，采樣頻率設(shè)置為100 kHz，采樣周期為10 μs，仿真步長為10 μs，設(shè)置1 s 時刻發(fā)生故障。

5.1 不同典型故障仿真分析

圖7 給出了區(qū)內(nèi)正極故障（距整流側(cè)750 km，如圖1中f2）的暫態(tài)電流行波波形圖。

圖7 區(qū)內(nèi)正極故障時暫態(tài)電流行波波形圖Fig.7 Waveform diagram of transient current traveling wave when an internal fault occurs in positive pole

由圖7 可知整流側(cè)正極暫態(tài)電流行波與逆變測正極暫態(tài)電流行波均發(fā)生正向突變，即故障為區(qū)內(nèi)故障。

圖8給出了區(qū)外整流側(cè)正極故障（如圖1中f1）的暫態(tài)電流行波波形圖。

圖8 區(qū)外整流側(cè)正極故障時暫態(tài)電流行波波形圖Fig.8 Waveform diagram of transient current traveling wave when an external fault occurs at rectifier side in positive pole

由圖8 可知，當發(fā)生區(qū)外整流側(cè)正極故障時整流側(cè)正極暫態(tài)電流行波發(fā)生負向突變，逆變側(cè)正極暫態(tài)電流行波發(fā)生正向突變。

圖9給出了區(qū)外逆變側(cè)正極故障（如圖1中f3）的暫態(tài)電流行波波形圖。

圖9 區(qū)外逆變側(cè)正極故障時暫態(tài)電流行波波形圖Fig.9 Waveform diagram of transient current traveling wave when an external fault occurs at inverter side in positive pole

由圖9 可知，當發(fā)生區(qū)外逆變側(cè)正極故障時逆變側(cè)正極暫態(tài)電流行波發(fā)生負向突變，整流側(cè)正極暫態(tài)電流行波發(fā)生正向突變。

5.2 直流線路不同過渡電阻故障選極仿真分析

4.4.1 節(jié)分析了過渡電阻對所提出的保護的影響。以下對不同故障類型分別設(shè)置不同過渡電阻進行仿真驗證，如表5和圖10所示。

圖10 不同故障類型不同過渡電阻下故障選極仿真結(jié)果圖Fig.10 Simulation results of fault pole selection under different fault types and different transition resistors

表5 不同過渡電阻作用下故障選極仿真結(jié)果分析（750 km）Tab. 5 Analysis of simulation results of fault pole selection under different transition resistors（750 km）

表5中iM+和iM-分別為整流側(cè)正極和負極暫態(tài)電流行波突變的最大值。iN+和iN-分別為逆變側(cè)正極和負極暫態(tài)電流行波突變的最大值。βM和βN分別為整流側(cè)和逆變側(cè)故障極識別參數(shù)。

本文以故障發(fā)生在距離整流側(cè)750 km 為例，考慮了不同故障類型（如圖1 中f2、f5、f6）分別在不同故障電阻（0 Ω、100 Ω、200 Ω、300 Ω、500 Ω）情況下的仿真結(jié)果。

由表5和圖10可知，系統(tǒng)在不同故障類型時在不同過渡電阻下β值均處于設(shè)定的選極算法門檻值范圍區(qū)間中，故所提出的保護對于不同的故障類型、不同的故障電阻下的故障均具有正確識別故障方向和故障極點的能力。

5.3 直流線路不同故障距離位置選極仿真分析

4.4.2 節(jié)分析了故障距離位置對所提出的保護的影響。分別對不同類型故障以及不同距離位置故障進行仿真驗證，如表6和圖11所示。

圖11 不同故障距離位置下故障選極仿真結(jié)果圖Fig.11 Simulation results of fault pole selection under different fault distance positions

表6 不同故障距離位置下故障選極仿真結(jié)果分析Tab. 6 Analysis of simulation results of fault pole selection at different fault distances

本節(jié)考慮了不同故障類型（如圖1 中f2、f5、f6）分別在不同故障距離位置（距離整流側(cè)150 km、375 km、1 125 km、1 350 km）下的仿真結(jié)果。（距離整流側(cè)750 km 情況下的故障已在上節(jié)中詳細敘述，在此不再贅述。）

由表6和圖11可知，系統(tǒng)在不同故障類型以及不同故障距離下的β值均處于設(shè)定的選極算法門檻值范圍區(qū)間中，故所提出的保護對于不同類型故障下、不同距離故障均具有正確識別故障方向和故障極點的能力。

6 結(jié)語

通過對行波在直流輸電系統(tǒng)中的傳播特性進行分析，本文提出了一種基于暫態(tài)電流行波突變的LCC-MMC 混合雙極直流輸電系統(tǒng)縱聯(lián)方向保護原理，利用整流側(cè)和逆變側(cè)在不同故障下暫態(tài)電流行波變化方向的不同來判別區(qū)內(nèi)外故障，并通過暫態(tài)電流行波突變最大值的比值來對故障極進行判別。與其他傳統(tǒng)直流線路保護原理相比，該保護原理簡單，故障后保護動作迅速，且門檻值易整定，保護受故障類型、過渡電阻以及故障發(fā)生距離位置的影響較小，具有一定的實用性。