袁 林 ，楊洪耕 ，王智琦 ，肖楚鵬

（1.四川大學(xué) 電氣信息學(xué)院，四川 成都 610065；2.南瑞集團公司（國網(wǎng)電力科學(xué)研究院），江蘇 南京 211106；3.南瑞（武漢）電氣設(shè)備與工程能效測評中心，湖北 武漢 430074）

0 引言

流過線路元件（如變壓器）的諧波電流是造成其損耗增大、運行效率降低等不良影響的主要因素之一，不利于電網(wǎng)經(jīng)濟高效運行，更使得供電部門承受巨大的經(jīng)濟損失［1-6］。目前，少有文獻從機理上研究判別流過元件的諧波電流主要來源是系統(tǒng)側(cè)還是用戶側(cè)的方法，且工程中也時常關(guān)注該問題［4-5，7］?；诖?，本文從機理出發(fā)，結(jié)合基本電路定律推導(dǎo)出能直接通過簡單條件定性判別流過線路元件主諧波電流來源的方法，為進一步量化兩側(cè)貢獻的諧波電流提供理論基礎(chǔ)。

主諧波電流來源判別的核心是比較系統(tǒng)側(cè)或用戶側(cè)各自發(fā)射諧波電流在公共連接點（PCC）的貢獻度大小，即以諧波電流指標為基本準則［2，8-9］。相關(guān)的判別方法多數(shù)先估計系統(tǒng)側(cè)或用戶側(cè)諧波阻抗［10-15］，然后進行判別。主要方法有波動量法［10-12］、門限電壓法［13］、獨立隨機矢量協(xié)方差法［14］和快速獨立分量分析法［15］等。波動量法根據(jù)PCC測量的諧波電壓和電流信息變化估計諧波阻抗，根據(jù)電網(wǎng)阻抗的基本屬性區(qū)分估計結(jié)果是系統(tǒng)側(cè)還是用戶側(cè)。獨立隨機矢量協(xié)方差法利用電網(wǎng)中諧波電壓和電流近似獨立、協(xié)方差為0的基本屬性，估計得到諧波阻抗值，在一定程度上削弱背景諧波的影響。門限電壓法是利用估計的諧波阻抗值來判斷的方法?？焖侏毩⒎至糠治龇ㄊ敲ぴ捶蛛x方法的一種，此方法應(yīng)用于諧波阻抗的估計，相比前述的方法有更好的效果。主諧波電流來源判別的本質(zhì)是得到定性結(jié)論，但上述方法無論判斷系統(tǒng)側(cè)還是用戶側(cè)為主諧波電流來源均需先對歷史數(shù)據(jù)進行計算估計，無法通過測量的諧波電壓和電流信息直接得出判別結(jié)論。

諧波有功功率方向法是在諧波狀態(tài)估計的基礎(chǔ)上提出的，該方法簡單易用，曾在實際電網(wǎng)中廣泛使用［2，8，16］，但有相關(guān)文獻指出，該方法判斷結(jié)果受系統(tǒng)側(cè)和用戶側(cè)諧波源相角差影響嚴重，并且兩側(cè)諧波源相角差可在整個圓周區(qū)間任意變化且無法直接測量，因此該方法的使用受到限制［16］。對于該方法能完全正確適用的廣泛區(qū)間，尚未有文獻進行研究。

基于此，本文基于諧波有功功率方向法的思路，在電網(wǎng)物理條件限制的基礎(chǔ)上討論系統(tǒng)側(cè)和用戶側(cè)諧波阻抗性質(zhì)的各種情況，理論分析并得到了對應(yīng)情況下的區(qū)間約束判別方法。對于主諧波電流來源為用戶側(cè)的情況，僅根據(jù)諧波電壓和電流的相角差就可直接判別；而對于主諧波電流來源為系統(tǒng)側(cè)的情況，在諧波電壓和電流的相角差的基礎(chǔ)上結(jié)合相關(guān)文獻的方法大致估計用戶側(cè)阻抗范圍后進行綜合判別。最后，利用多個工程場景的實測數(shù)據(jù)對本文方法進行了驗證分析，結(jié)果表明本文方法能有效地進行判別。

1 諧波有功功率方向法及分析

1.1 諧波有功功率方向法

圖1為諾頓等效電路模型，是有功功率方向法的基礎(chǔ)。圖中，Is和Ic分別為系統(tǒng)側(cè)和用戶側(cè)諧波源電流，Is和Ic分別為兩者模值大??；Zs和Zc分別為系統(tǒng)側(cè)和用戶側(cè)諧波阻抗；Upcc和Ipcc分別為監(jiān)測點PCC處測得的諧波電壓和諧波電流，Upcc和Ipcc分別為兩者模值大小。

圖1 諾頓等效電路模型Fig.1 Norton equivalent circuit model

式（1）為諧波有功功率的計算公式，可以看出，諧波有功功率方向法的本質(zhì)是根據(jù)監(jiān)測點（一般為PCC）諧波電壓Upcc和電流Ipcc的相角差φ的余弦值正負來判斷哪一側(cè)的諧波電流貢獻大。

諧波有功功率方向法認為，諧波的有功功率由系統(tǒng)側(cè)流向用戶側(cè)，即P＞0時，系統(tǒng)側(cè)諧波電流貢獻大；當諧波的有功功率由用戶側(cè)流向系統(tǒng)側(cè)，即P＜0時，用戶側(cè)諧波電流貢獻大［2，16］。

1.2 諧波電流指標

根據(jù)電路的疊加定理，PCC處的諧波電流Ipcc是由系統(tǒng)側(cè)貢獻諧波電流Ispcc和用戶側(cè)貢獻諧波電流Icpcc兩部分疊加而成，如式（2）所示。諧波電流Ipcc相量分解如圖2所示，Isf和Icf分別為Ispcc和Icpcc在Ipcc方向的投影大小，其各自與Ipcc的比值為諧波電流指標。由余弦定理可得：

圖2 電流相量分解Fig.2 Decomposition of current phasor

由上式可以得出，諧波電流指標的本質(zhì)是判斷和的大小。

1.3 諧波有功功率方向法適用分析

文獻［16］指出，諧波有功功率方向法并不完全有效，存在判斷誤區(qū)。實際電網(wǎng)能得到PCC處測量的各次諧波電壓和電流大小及相角差，且相角差可在0°～360°區(qū)間變化。以表1實驗參數(shù)為例，以PCC處諧波電壓Upcc和電流Ipcc的相角差φ為變化量，步長取1°，得到的實驗結(jié)果如圖3所示。

從圖中結(jié)果可以看出，區(qū)間1中，且P＞0；區(qū)間2內(nèi)，且 P＜0，符合諧波有功功率方向法的判斷。即諧波有功功率方向法在一定的區(qū)間內(nèi)是能完全正確判斷哪一側(cè)諧波電流指標更大，存在一定的適用區(qū)間。

表1 實驗參數(shù)Table 1 Experimental parameters

圖3 諧波有功功率方向法判斷條件Fig.3 Judgment conditions of harmonic active power direction method

2 基于區(qū)間約束的主諧波電流來源判別

2.1 基本原理

主諧波電流來源判別的目的是比較兩側(cè)諧波電流指標的大小，實質(zhì)是比較和的大小。根據(jù)基本電路原理，從圖1可以得出：

上式兩邊同時取模值，并平方可得：

其中，θs和θc分別為系統(tǒng)側(cè)和用戶側(cè)諧波阻抗角；k為用戶側(cè)與系統(tǒng)側(cè)諧波阻抗模值的比值，一般k＞1。

實際電網(wǎng)物理限制使得諧波阻抗的實部一定為正值，也就是諧波阻抗角 θs和θc介于［-90°，90°］之間，理論上，阻抗角介于［0°，90°］為感性，介于［-90°，0°］為容性［10］。而電網(wǎng)中元件和負荷大多呈感性，對于系統(tǒng)側(cè)阻抗而言，其性質(zhì)通常也為感性且感抗分量與電阻分量的比值大于1，并且電壓等級越高這個比值越大［17］，易得出其阻抗角介于［45°，90°］區(qū)間。對于用戶側(cè)阻抗而言，文獻［18-19］給出了負荷的諧波阻抗模型及其計算方法，在不考慮無功補償裝置和城市配電網(wǎng)大量電纜的情況下，配電網(wǎng)負荷的功率因數(shù)大多數(shù)小于等于 0.85［20］，取負荷功率因數(shù)為 0.85，各次諧波的阻抗角如圖4所示?？梢钥闯?，各次諧波的阻抗角均介于［45°，90°］區(qū)間且阻抗性質(zhì)為感性，在功率因數(shù)小于0.85的其他情況也能得到相同的結(jié)論。據(jù)文獻［21-22］，在大量使用電纜和中低壓電容器等無功補償裝置的城市配電網(wǎng)中，諧振點僅可能出現(xiàn)在高次諧波，而導(dǎo)致中低次諧波的阻抗性質(zhì)為容性，但此時一般容抗分量小于電阻分量，因此諧波阻抗角介于［-45°，0°］區(qū)間。

圖4 各次諧波阻抗相角Fig.4 Phase angle of each harmonic impedance

綜上，為考慮兩側(cè)諧波阻抗特性，滿足電網(wǎng)中常見的工程場景，且使得方法在使用過程中簡單易行，在本文的分析研究中，諧波阻抗性質(zhì)為感性時，視作其阻抗角介于［45°，90°］區(qū)間；諧波阻抗性質(zhì)為容性時，視作其阻抗角介于［-45°，0°］區(qū)間。

2.2 主諧波電流來源為用戶側(cè)的約束條件

為得到主諧波電流來源為用戶側(cè)的約束條件，將式（6）與式（5）作差可得式（8），若其等號右邊部分大于0成立，則，即主諧波電流來源為用戶側(cè)。以諧波阻抗性質(zhì)相同和不同討論以下4種情況。

a.系統(tǒng)側(cè)和用戶側(cè)諧波阻抗性質(zhì)同為感性，即θs和θc均位于［45°，90°］，且相角差φ位于［180°，270°）區(qū)間，則可得 cos（φ-θs）＜0 且 cos（φ-θc）＜0；由k＞1 可得 1-1/k2＞0。因此，式（8）等號右邊部分大于0恒成立，得到主諧波電流來源為用戶側(cè)的結(jié)論。

b.系統(tǒng)側(cè)和用戶側(cè)諧波阻抗性質(zhì)均為容性，即θs和θc均位于［-45°，0°］，主諧波電流來源為用戶側(cè)的區(qū)間約束為相角差φ位于（90°，225°］區(qū)間。

c.系統(tǒng)側(cè)諧波阻抗性質(zhì)為感性和用戶側(cè)為容性，即 θs介于［45°，90°］，θc介于［-45°，0°］，相角差φ介于［180°，225°］區(qū)間，得cos（φ-θc）＜0、cos（φ-θs）＜0；由k＞1可得 1-1/k2＞0。故式（8）等號右邊部分大于0恒成立。所以，主諧波電流來源為用戶側(cè)的區(qū)間約束為相角差φ位于［180°，225°］區(qū)間。

d.系統(tǒng)側(cè)諧波阻抗性質(zhì)為容性和用戶側(cè)為感性，即 θs位于［-45°，0°］，θc位于［45°，90°］，主諧波電流來源為用戶側(cè)區(qū)間約束為相角差φ位于［180°，225°］區(qū)間。

綜上，基于圖1所示的參考方向和兩側(cè)各種阻抗性質(zhì)的各類情況，得到主諧波電流來源為用戶側(cè)的區(qū)間約束條件，如表2所示。

表2 主諧波電流來源為用戶側(cè)的約束條件Table 2 Restrictions of harmonic current mainly origins from user side

2.3 主諧波電流來源為系統(tǒng)側(cè)的約束條件

為得到主諧波電流來源為系統(tǒng)側(cè)的約束條件，將式（5）與式（6）作差可得：

若式（9）等號右邊部分大于 0成立，則，即主諧波電流來源為系統(tǒng)側(cè)。以諧波阻抗性質(zhì)相同和不同討論以下4種情況。

a.系統(tǒng)側(cè)和用戶側(cè)諧波阻抗性質(zhì)同為感性，即θs和θc均位于［45°，90°］，相角差 φ 位于［0°，90°）區(qū)間，易得 cos（φ-θs）＞0 且 cos（φ-θc）＞0 但 1 /k2-1＜0，因此無法直接得出結(jié)論。

根據(jù)余弦相關(guān)知識易得θs在式（10）條件下，取最小值；同理，θc在式（11）條件下，取最大值。只要滿足條件就可得到結(jié)論。

令式（9）大于 0，結(jié)合式（10）、（11），可得到化簡后附加約束條件：

b.系統(tǒng)側(cè)和用戶側(cè)諧波阻抗性質(zhì)同為容性，即θs和θc均在［-45°，0°］，主諧波電流來源為系統(tǒng)側(cè)的約束是相角差φ位于（-90°，-22.5°］區(qū)間且滿足式（13）所示附加約束條件。

c.系統(tǒng)側(cè)諧波阻抗性質(zhì)為感性和用戶側(cè)為容性，即 θs位于［45°，90°］，θc位于［-45°，0°］，且相角差 φ位于［0°，45°］區(qū)間。易得 θs為 90°時，??；θc為 -45°時，取。令式（9）大于 0可得到化簡后附加約束條件如式（14）所示。

d.系統(tǒng)側(cè)諧波阻抗性質(zhì)為容性和用戶側(cè)為感性，即 θs在［-45°，0°］，θc在［45°，90°］，主諧波電流來源為系統(tǒng)側(cè)的約束是相角差φ介于［0°，45°］區(qū)間且滿足式（15）所示附加約束條件。

通過上文的分析可知，在主諧波電流來源為系統(tǒng)側(cè)的約束的附加條件中含有用戶側(cè)諧波阻抗模值，僅在用戶側(cè)諧波源波動較小時，可以得到較好的阻抗估計結(jié)果［10-11，13］。因此，本文僅利用相關(guān)方法估計用戶側(cè)諧波阻抗模值范圍進行輔助判斷。

實際中低壓配電網(wǎng)的電壓波形通常為平頂波，因此可以用諧波諾頓電路來等效用戶側(cè)模型。文獻［14-15，23］是近年來新提出的諧波阻抗計算的方法，在測得的諧波電壓和電流的相角差φ滿足區(qū)間約束條件，且用戶側(cè)諧波波動較小時，可以較精確地計算用戶側(cè)諧波阻抗模值。同時為盡量使本文方法在應(yīng)用中簡單易行，本文采用文獻［15］方法對用戶側(cè)諧波阻抗的模值進行計算。

綜上所述，根據(jù)圖1所示的參考方向和兩側(cè)各種阻抗性質(zhì)分類，得到主諧波電流來源為系統(tǒng)側(cè)的區(qū)間約束條件，如表3所示。

表3 主諧波電流來源為系統(tǒng)側(cè)的約束條件Table 3 Restrictions of harmonic current mainly origins from system side

3 實例驗證

3.1 電弧爐負荷算例

圖5為來自某工業(yè)電弧爐廠150 kV母線（PCC）的3次諧波電壓和電流實測數(shù)據(jù)，這一段時間內(nèi)有約1 h的停機時段，且已知PCC兩側(cè)的該次諧波阻抗性質(zhì)均呈感性。

圖5 3次諧波電壓和電流實測數(shù)據(jù)Fig.5 Measured data of 3rd harmonic voltage and current

圖5中，由于電弧爐在點弧和熄弧工作狀態(tài)頻繁變化，0～600 min工作時間段諧波電壓和電流幅值波動十分明顯，但諧波電壓和電流的相角差幾乎全集中在［210°，270°］的 60°區(qū)間。根據(jù)表2 的約束條件，可得用戶側(cè)在PCC貢獻的諧波電流更大，主諧波電流來源為用戶側(cè)的結(jié)論。眾所周知，采用晶閘管整流的電弧爐負荷在工作時會發(fā)射大量的諧波電流，特別是中低次諧波［6］。同時在600～660 min的停機時段，PCC處的諧波電壓和電流幅值急劇減小且兩者相角差集中在［330°，360°］區(qū)間，不在表2 的約束條件內(nèi)，也可得到用戶側(cè)對PCC處該次諧波電流指標貢獻更大；并且文獻［11］的評估結(jié)果顯示用戶側(cè)發(fā)射的諧波水平更大。所以，以上的分析與結(jié)果驗證了本文方法的正確性與有效性。

3.2 城市配電網(wǎng)算例

圖6為某城市電網(wǎng)臨近直流落點的66 kV母線（PCC）11次諧波電壓和電流某天的5 h實測數(shù)據(jù)，整個數(shù)據(jù)區(qū)間有約1h的低谷用電時段。已知PCC等效的系統(tǒng)側(cè)諧波阻抗性質(zhì)為感性，由于該城市電網(wǎng)大量使用地下電纜，同時為保證用戶電能質(zhì)量裝有大量的電容器等裝置，而導(dǎo)致用戶側(cè)該次諧波阻抗性質(zhì)為容性［17］。

圖6 11次諧波電壓和電流實測數(shù)據(jù)Fig.6 Measured data of 11th harmonic voltage and current

從圖6中可以看出，0～4 h時間段的正常用電時段，諧波電壓和電流幅值變化很小，這也是用戶側(cè)諧波阻抗估計的基本條件，同時可知諧波電壓和電流的相角差集中在［0°，30°］區(qū)間。利用表3的約束條件，可知當前情況滿足主諧波電流來源為系統(tǒng)側(cè)的相角差φ的范圍，還需要滿足式（14）所示附加條件。采用文獻［15］的方法對用戶側(cè)諧波阻抗模值進行計算，得到該時間段內(nèi)95%概率大值的PCC相關(guān)參數(shù)如表4所示，代入式（14），易知滿足此附加約束條件，因此可以判斷當前情況下主諧波電流來源為系統(tǒng)側(cè)。

表4 PCC相關(guān)參數(shù)Table 4 Relative parameters of PCC

從4～5 h的低谷時段可以看出，其相比于正常時段PCC處該次諧波電壓和電流幅值均沒有明顯的減小，諧波電壓與電流的相角差仍集中在［0°，30°］區(qū)間。文獻［17］采用Digsilent對該城市電網(wǎng)進行仿真計算，得到臨近直流落點的各支路母線11次諧波電流測量值的約70%由500 kV直流落點貢獻，因此可確定系統(tǒng)側(cè)對PCC處11次諧波電流貢獻更大，主諧波電流來源為系統(tǒng)側(cè)。所以，上述結(jié)果驗證了本文方法的正確性與有效性。

4 結(jié)論

對諧波有功功率方向法進行分析，并從原理上推導(dǎo)得到基于區(qū)間約束的主諧波電流來源判別方法，得出以下結(jié)論。

a.分析并得出諧波有功功率方向法在一定區(qū)域內(nèi)可以完全準確地得到結(jié)論。結(jié)合電路基本原理和嚴格不等式約束推證得到的區(qū)間約束條件，可以有效地解決實際工程中難以判斷流過線路元件諧波電流主要來源的問題。

b.在電網(wǎng)基本特性下，廣泛研究了兩側(cè)諧波阻抗性質(zhì)的各類情況，能滿足工程使用的大多數(shù)場景。最后，使用多個工程場景實測數(shù)據(jù)分析結(jié)果驗證了方法的正確性與有效性。

從機理推證的區(qū)間約束條件判別主諧波電流來源的方法，為諧波有功功率方向法的使用和定量劃分電力設(shè)備因諧波電流造成損失提供了理論基礎(chǔ)，有較大的理論意義和工程實用價值。

參考文獻：

［1］WILSUN X，YILU L.A method for determining customer and utility harmonic contributions at the point of common coupling［J］.IEEE Transactions on Power Delivery，2000，15（2）：804-811.

［2］WILSUN X，XIAN L，YILU L.An investigation on the validity of power-direction method for harmonic source determination［J］.IEEE Transactions on Power Delivery，2003，18（1）：214-219.

［3］LI C ，XU W，TAYJASANANT T.A critical impedance based method for identifying harmonic sources［J］.IEEE Transactions on Power Delivery，2004，19（2）：671-678.

［4］劉華，劉銳，付洋洋，等.基于阻抗邊界法的換流變壓器箱殼諧波損耗計算［J］.高電壓技術(shù)，2015，41（9）：3171-3176.LIU Hua，LIU Rui，F(xiàn)U Yangyang，et al.Calculation of harmonic loss in converter transformer tank shells based on the impedance boundary condition［J］.High Voltage Engineering，2015，41（9）：3171-3176.

［5］葉建新，黃健豪，張祿亮，等.東莞配電網(wǎng)諧波對配電變壓器危害分析［J］.電力自動化設(shè)備，2011，31（7）：146-149.YE Jianxin，HUANG Jiaohao，ZHANG Luliang，et al.Analysis of harmonic harm to distribution transformers［J］.Electric Power Automation Equipment，2011，31（7）：146-149.

［6］郎維川.供電系統(tǒng)諧波的產(chǎn)生、危害及其防護對策［J］.高電壓技術(shù)，2002，28（6）：30-31.LANG Weichuan.Generation and harm of harmonic in power supply system and the countermeasure［J］.High Voltage Engineering，2002，28（6）：30-31.

［7］張良縣，陳模生，彭宗仁，等.非正弦負載電流下特高壓換流變壓器繞組的諧波損耗分析［J］.中國電機工程學(xué)報，2014，34（15）：2452-2459.ZHANG Liangxian，CHEN Mosheng，PENG Zongren，et al.Study on the harmonic losses of UHV converter transformer windings subject to non-sinusoidal load current［J］.Proceedings of the CSEE，2014，34（15）：2452-2459.

［8］周林，張鳳，栗秋華，等.配電網(wǎng)中諧波源定位方法綜述［J］.高電壓技術(shù)，2007，33（5）：103-108.ZHOU Lin，ZHANG Feng，LI Qiuhua，et al.Methods for localizing the harmonic source in power distribution network ［J］.High Voltage Engineering，2007，33（5）：103-108.

［9］華回春，賈秀芳，曹東升，等.電能質(zhì)量數(shù)據(jù)交換格式下的諧波責(zé)任估計［J］.電網(wǎng)技術(shù)，2013，37（11）：3110-3117.HUA Huichun，JIA Xiufang，CAO Dongsheng，etal.Harmonic contribution estimation underpowerqualitydatainterchange format［J］.Power System Technology，2013，37（11）：3110-3117.

［10］YANG H，PIROTTE P，ROBERT A.Harmonic emission levels of industrial loads statistical assessment［J］.CIGRE，1996，96：1936-1996.

［11］惠錦，楊洪耕，葉茂清.基于阻抗歸一化趨勢判別的諧波發(fā)射水平估計［J］.中國電機工程學(xué)報，2011，31（10）：73-80.HUIJin，YANG Honggeng，YE Maoqing.Assessingharmonic emission level based on the impedance gathering trend discrimination［J］.Proceedings of the CSEE，2011，31（10）：73-80.

［12］車權(quán)，楊洪耕.基于穩(wěn)健回歸的諧波發(fā)射水平估計方法［J］.中國電機工程學(xué)報，2004，24（4）：39-42.CHE Quan，YANG Honggeng.Assessing the harmonic emission level based on robust regression method［J］.Proceedings of the CSEE，2004，24（4）：39-42.

［13］孫媛媛，李培鑫，尹志明.基于門限電壓法的諧波源定位［J］.電力系統(tǒng)自動化，2015，39（23）：145-151.SUN Yuanyuan，LI Peixin，YIN Zhiming.Harmonic source identification based on threshold voltage method［J］.Automation of Electric Power Systems，2015，39（23）：145-151.

［14］惠錦，楊洪耕，林順富，等.基于獨立隨機矢量協(xié)方差特性的諧波發(fā)射水平評估方法［J］.電力系統(tǒng)自動化，2009，33（7）：27-31.HUI Jin，YANG Honggeng，LIN Shunfu，et al.Assessment method ofharmonic emission level based on covariance characteristic ofrandom vectors［J］.Automation ofElectric Power Systems，2009，33（7）：27-31.

［15］趙熙，楊洪耕.基于快速獨立分量分析的系統(tǒng)側(cè)諧波阻抗計算方法［J］.電力系統(tǒng)自動化，2015，39（23）：139-144.ZHAO Xi，YANG Honggeng.Method of calculating system-side harmonic impedance based on FastICA［J］.Automation of Electric Power Systems，2015，39（23）：139-144.

［16］XU W.Power direction method cannot be used for harmonic source detection［C］∥Power Engineering Society Summer Meeting.［S.l.］：IEEE，2000：873-876.

［17］李勇，徐友平，肖華，等.華中電網(wǎng)穩(wěn)定計算用負荷模型參數(shù)仿真研究［J］.電網(wǎng)技術(shù)，2007，31（5）：17-20.LI Yong，XU Youping，XIAO Hua，et al.Simulation study on load model parameters for stability analysis of central china power grid［J］.Power System Technology，2007，31（5）：17-20.

［18］OBERT A，DEFLANDRE T.Guide for assessing the network harmonic impedance［J］.Electra，1996（167）：96-131.

［19］WU X，SADULLAH S，MATTHEWS B，et al.Nodal harmonic impedance derivation of AC network in PSS /E［C］∥9th IET International Conference on AC and DC Power Transmission.［S.l.］：IET，2010：1-5.

［20］李建超.配電網(wǎng)諧波調(diào)查及諧波網(wǎng)損計算［D］.濟南：山東大學(xué)，2008.LIJianchao.Investigation ofharmonics and calculation of harmonic losses in distribution systems［D］.Ji’nan：Shandong University，2008.

［21］陳甜甜，胡蓉，金祖洋，等.鄰近直流輸電落點的城市電網(wǎng)諧波傳導(dǎo)特性分析［J］.電網(wǎng)技術(shù)，2015，39（10）：3000-3004.CHEN Tiantian，HU Rong，JIN Zuyang，et al.Analysis research on harmonic transmission characteristics of urban power network adjacent to UHVDC terminal location［J］.Power System Technology，2015，39（10）：3000-3004.

［22］王建勛，劉會金，劉春陽.基于復(fù)合判據(jù)的諧波/間諧波源識別方法［J］.電力自動化設(shè)備，2013，33（7）：63-69.WANG Jianxun，LIU Huijin，LIU Chunyang.Identification of harmonic and interharmonic sources based on composite criterion［J］.Electric Power Automation Equipment，2013，33（7）：63-69.

［23］邱思語，楊洪耕.改進的加權(quán)支持向量機回歸的諧波發(fā)射水平估計方法［J］.電工技術(shù)學(xué)報，2016，31（5）：85-90.QIU Siyu，YANG Honggeng.Assessment method of harmonic emission level based on the improved weighted support vector machine regression［J］.Transactions of China Electrotechnical Society，2016，31（5）：85-90.