王智琦 ，楊洪耕，肖楚鵬 ，冷 月

（1.四川大學(xué) 電氣信息學(xué)院，四川 成都 610065；2.南瑞（武漢）電氣設(shè)備與工程能效測評中心，湖北 武漢 430074）

0 引言

電網(wǎng)中非線性負(fù)荷的不斷接入以及電網(wǎng)中諧波源的復(fù)雜性和多變性，使諧波溯源成為諧波研究中的一個(gè)難點(diǎn)［1］。目前主要的諧波溯源方法分為兩大類：一類是阻抗參數(shù)法，另一類是諧波功率方向法。

當(dāng)采用“干預(yù)式”方法時(shí)，利用人為方式產(chǎn)生的附加擾動(dòng)量計(jì)算系統(tǒng)側(cè)諧波阻抗，得到的結(jié)果較可靠。但這種擾動(dòng)很可能會對系統(tǒng)的正常運(yùn)行產(chǎn)生不利影響，另外，在很多場合不存在可投切的電氣元件和設(shè)備［2］。因此，工程上常采用“非干預(yù)式”方法，利用公共連接點(diǎn)（PCC）處諧波電壓電流錄波數(shù)據(jù)計(jì)算系統(tǒng)側(cè)諧波阻抗。典型的方法有波動(dòng)量法［2-3］、回歸法［4-7］。這些方法有效的前提在于找到系統(tǒng)背景諧波不變的樣本，在背景諧波波動(dòng)較大時(shí)，計(jì)算結(jié)果往往不能滿足工程精度的要求。隨機(jī)獨(dú)立矢量法［8-9］、獨(dú)立分量法［10］、支持向量機(jī)［11］等方法在一定程度上消除了背景諧波波動(dòng)的影響，但并無實(shí)質(zhì)上的改變，當(dāng)PCC處諧波主要不是用戶側(cè)諧波造成時(shí)，這類方法不能得到精確的計(jì)算結(jié)果［12］，尤其是阻抗角的精確值更難獲得。由于電網(wǎng)中阻抗參數(shù)不易精確獲得，而諧波功率方向是可直接根據(jù)測量數(shù)據(jù)得到的信息，所以工程上一直廣泛采用更簡單的諧波功率方向法。

文獻(xiàn)［13］指出使用諧波有功功率方向法判斷主諧波源位置具有局限性，該方法受PCC兩側(cè)諧波源相角差影響較大，不能保證諧波有功功率方向法始終適用。為此，文獻(xiàn)［14］提出了基于無功功率方向判斷主諧波源位置的方法，但該方法不僅與兩側(cè)等值諧波電壓源幅值有關(guān)而且還與兩側(cè)諧波阻抗性質(zhì)有關(guān)。當(dāng)以兩側(cè)諧波電壓源幅值大小為判斷指標(biāo)時(shí)，在兩側(cè)阻抗同呈感性的情況下僅能判斷用戶側(cè)為主諧波源的情況，對于其他情況存在不確定性。且國內(nèi)外尚無以諧波電壓貢獻(xiàn)作為考核指標(biāo)時(shí)無功功率方向法的可行性研究。

由于阻抗參數(shù)法和諧波功率方向法均無法在任何情況下都具有較高的計(jì)算精度，除個(gè)別大型諧波源負(fù)荷工作人員能根據(jù)運(yùn)行經(jīng)驗(yàn)結(jié)合測量數(shù)據(jù)估計(jì)出主諧波源位置外，系統(tǒng)中還有許多PCC是由不同類型和容量的用電設(shè)備組成的，在這樣的實(shí)際情況中，不能直接根據(jù)測量數(shù)據(jù)判斷主諧波源位置。

針對上述問題，本文提出一種基于區(qū)間約束的諧波溯源方法。通過仿真分析得出無功功率方向法更適用于解決PCC處諧波溯源問題。以兩側(cè)諧波電壓貢獻(xiàn)度大小作為判斷兩側(cè)諧波源作用大小的依據(jù)，根據(jù)兩側(cè)諧波電壓貢獻(xiàn)度隨諧波電壓電流夾角的變化規(guī)律，分區(qū)間推導(dǎo)得到判斷系統(tǒng)側(cè)和用戶側(cè)為主諧波源時(shí)各區(qū)間的嚴(yán)格不等式約束條件。本文以諧波無功功率方向法為基礎(chǔ)，僅需根據(jù)諧波電壓電流測量信息并結(jié)合簡單的區(qū)間約束條件即可確定主諧波源。為提高溯源結(jié)果的可靠性，利用可信區(qū)間去除較嚴(yán)重的異常監(jiān)測數(shù)據(jù)。通過多場景實(shí)際工程數(shù)據(jù)驗(yàn)證了本文方法的準(zhǔn)確性和有效性。

1 諧波溯源方法及分析

1.1 諧波無功功率方向

諧波無功功率方向法的基本原理是通過測量PCC處諧波無功功率方向來追溯主諧波源。圖1為兩側(cè)諧波源分析模型。圖中，Iu、Ic分別為系統(tǒng)側(cè)和用戶側(cè)等值諧波電流源；Zu、Zc分別為系統(tǒng)側(cè)和用戶側(cè)等值諧波阻抗；以諧波電壓Upcc作為參考相量，Upcc=Upcc∠0°，則諧波電流 Ipcc=Ipcc∠φ，φ 為諧波電壓和電流的夾角。兩側(cè)諧波電流源幅值為Iu和Ic?？汕蟪鯬CC處諧波無功功率為：

其中，Q為諧波無功功率。實(shí)際電網(wǎng)中，一般系統(tǒng)側(cè)諧波阻抗呈感性，若Q>0，則用戶側(cè)為主諧波源；若Q<0，則系統(tǒng)側(cè)為主諧波源。

圖1 諧波溯源模型Fig.1 Model of harmonic source tracing

1.2 諧波溯源分析

根據(jù)圖1的諧波源模型，來用電壓矢量疊加原理對PCC處諧波電壓進(jìn)行分解，得到系統(tǒng)和用戶兩側(cè)在PCC處引起的諧波電壓貢獻(xiàn)度，進(jìn)而判斷主諧波源位置。比較兩側(cè)諧波電壓貢獻(xiàn)度大小的實(shí)質(zhì)是比較兩側(cè)諧波電流源幅值的大小，即比較Iu和Ic大小。若用戶側(cè)諧波電流源幅值較大，則用戶側(cè)為主諧波源；反之，系統(tǒng)側(cè)為主諧波源。

系統(tǒng)和用戶兩側(cè)諧波源相位變化會影響PCC處諧波電壓和電流的夾角φ，且實(shí)際情況下兩者相位均可為任意角度，需分析夾角φ為任意值的情況?？紤]實(shí)際電網(wǎng)運(yùn)行中并以兩側(cè)阻抗同為感性的情況為例，利用表1參數(shù)［13］進(jìn)行諧波溯源分析，其結(jié)果如圖2所示。

表1 實(shí)驗(yàn)參數(shù)Table 1 Experiment parameters

圖2 諧波溯源分析Fig.2 Analysis of harmonic source tracing

根據(jù)圖2仿真分析結(jié)果可知，利用諧波無功功率方向法判斷用戶側(cè)為主諧波源的條件是夾角φ位于第一、二象限，系統(tǒng)側(cè)為主諧波源的條件是夾角φ位于第三、四象限。而利用諧波電壓貢獻(xiàn)度法判斷用戶側(cè)為主諧波源的條件是夾角φ位于兩臨界角之間，判斷區(qū)域近似位于第一、二象限；判斷系統(tǒng)側(cè)為主諧波源的條件是夾角φ位于兩臨界角之外，判斷區(qū)域近似位于第三、四象限。

雖然這2種方法的諧波溯源結(jié)果并不完全一致，但可認(rèn)為諧波無功功率方向法在大部分情況下能準(zhǔn)確判斷主諧波源位置。相比有功功率方向法，利用諧波無功功率方向來判斷主諧波源位置與諧波電壓貢獻(xiàn)度法的結(jié)果更準(zhǔn)確。根據(jù)無功功率方向僅能初步估計(jì)主諧波源位置，要得到準(zhǔn)確的諧波溯源結(jié)果，需在現(xiàn)有無功功率方向法的基礎(chǔ)上附加一定約束條件。

2 基于區(qū)間約束的諧波溯源方法

2.1 諧波溯源方法的基礎(chǔ)理論

根據(jù)諧波源等值模型，利用基爾霍夫電流定律得：

從而有：

圖3 諧波電壓貢獻(xiàn)度變化規(guī)律Fig.3 Variation rules of harmonic voltage contribution

其中，αu為系統(tǒng)側(cè)阻抗角；αc為用戶側(cè)阻抗角。式（3）、（4）等號右側(cè)的第一項(xiàng)、第二項(xiàng)構(gòu)成恒定分量，第三項(xiàng)為波動(dòng)分量。因?yàn)樗缘暮愣ǚ至啃∮诘暮愣ǚ至康牟▌?dòng)分量幅值小于的波動(dòng)分量幅值。若則系統(tǒng)側(cè)為主諧波源；反之，用戶側(cè)為主諧波源。根據(jù)式（3）、（4）可得出兩側(cè)諧波貢獻(xiàn)度與夾角φ阻抗角之間的變化關(guān)系，如圖3所示。當(dāng)阻抗角一定時(shí)的曲線隨夾角φ呈正弦規(guī)律變化，而阻抗角變化時(shí)，曲線發(fā)生相移。

大多數(shù)情況下電網(wǎng)的系統(tǒng)側(cè)和用戶側(cè)諧波阻抗均為感性；當(dāng)無功補(bǔ)償電容器選用電抗率為6%串聯(lián)電抗器時(shí)，在5次以上諧波阻抗呈感性；在35 kV或10 kV的城市配電網(wǎng)中通常使用幾千米的短電纜，僅在高次諧波（如29次以上［15］）處有可能會發(fā)生諧振，當(dāng)諧波頻率低于諧振點(diǎn)時(shí)用戶側(cè)諧波阻抗為容性。電網(wǎng)諧波阻抗具有一些固有特性，在一定條件下既可能呈感性，也可能呈容性，但多數(shù)條件下呈感性［12］。且系統(tǒng)阻抗與用戶阻抗相對大小幾乎是確定的。通常阻抗幅值的估計(jì)相比阻抗角的估計(jì)更為準(zhǔn)確。

2.2 判斷系統(tǒng)側(cè)為主諧波源

根據(jù)諧波無功功率方向法，當(dāng)系統(tǒng)側(cè)諧波阻抗為感性時(shí)，若測得PCC處諧波無功功率方向?yàn)樨?fù)，則系統(tǒng)側(cè)為主諧波源，此時(shí)夾角φ位于第三、四象限。由于難以得到使恒成立的參數(shù)條件，需要根據(jù)不同阻抗性質(zhì)及測量角度詳細(xì)討論并分類推導(dǎo)。

（1）兩側(cè)諧波阻抗同呈感性。

當(dāng)兩側(cè)諧波阻抗角未知時(shí)，不能直接判斷Iu和Ic的大小。的變化范圍如圖4所示的最大值、最小值分別為 Iumax、Iumin和 Icmax、Icmin。 對于測得的夾角φ，若任意阻抗角都使得Iumin>Icmax成立，則系統(tǒng)側(cè)為主諧波源。

圖4 兩側(cè)阻抗同呈感性時(shí)的趨勢圖Fig.4 Tendency chart when both sides are of inductive impedance

a.夾角φ位于第四象限。

因?yàn)榈暮愣ǚ至看笥诘暮愣ǚ至?，由圖4可得的恒定分量的恒定分量>Icmax，由不等式的傳遞性可以得到Iumin>Icmax恒成立，即當(dāng)夾角φ位于第四象限時(shí)，可直接判斷系統(tǒng)側(cè)為主諧波源。

b.夾角φ位于第三象限。

的最小值Iumin的表達(dá)式為：

的最大值Icmax的表達(dá)式為：

設(shè)的最小值大于的最大值，即Iumin>Icmax，可得：

其中，兩側(cè)阻抗幅值的比值由式（7）可得：

當(dāng)夾角φ位于第三象限且滿足式（8）范圍時(shí)，可判斷系統(tǒng)側(cè)為主諧波源。

（2）Zu呈感性，Zc呈容性。

a.夾角 φ 位于第四象限（315°，360°）。

圖5中，因?yàn)榈暮愣ǚ至看笥诘暮愣ǚ至?，且的波?dòng)分量幅值大于的波動(dòng)分量幅值，所以Iumin>Icmax，可判斷系統(tǒng)側(cè)為主諧波源。

圖5 系統(tǒng)側(cè)阻抗呈感性、用戶側(cè)阻抗呈容性時(shí)的趨勢圖Fig.5 Tendency chart when system-side is of inductive impedance and customer-side is of capacitive impedance

b.夾角 φ 位于第四象限（270°，315°）。

的最小值Iumin的表達(dá)式為：

的最大值Icmax的表達(dá)式為：

當(dāng)滿足式（11）時(shí)，可判斷系統(tǒng)側(cè)為主諧波源。

c.夾角φ位于第三象限。

的最小值Iumin的表達(dá)式為：

的最大值Icmax的表達(dá)式為：

設(shè)的最小值大于的最大值，即Iumin>Icmax，簡化后得：

當(dāng)滿足式（14）時(shí)，可判斷系統(tǒng)側(cè)為主諧波源。

所以，當(dāng)夾角φ位于第三、四象限時(shí)有如下結(jié)論。

a.兩側(cè)諧波阻抗同為感性時(shí)，當(dāng)夾角φ位于第四象限時(shí)可直接判斷系統(tǒng)側(cè)為主諧波源；當(dāng)夾角φ位于第三象限且滿足式（8）范圍時(shí)可判斷系統(tǒng)側(cè)為主諧波源。

b.Zu呈感性、Zc呈容性時(shí)，當(dāng)夾角φ位于第四象限（315°，360°）時(shí)可直接判斷系統(tǒng)側(cè)為主諧波源；當(dāng)夾角 φ 位于第四象限（270°，315°）且滿足式（11）時(shí)可判斷系統(tǒng)側(cè)為主諧波源；當(dāng)夾角φ位于第三象限且滿足式（14）時(shí)可判斷系統(tǒng)側(cè)為主諧波源。

2.3 判斷用戶側(cè)為主諧波源

同理，若測得PCC處諧波無功功率方向?yàn)檎瑒t用戶側(cè)為主諧波源，此時(shí)夾角φ位于第一、二象限。推導(dǎo)過程與判斷系統(tǒng)側(cè)為主諧波源類似。

a.兩側(cè)諧波阻抗同為感性。

當(dāng)夾角φ位于第二象限時(shí)，需滿足區(qū)間約束：

b.Zu呈感性，Zc呈容性。

當(dāng)夾角 φ 位于第二象限（135°，180°）時(shí)，需滿足區(qū)間約束：

而其他更復(fù)雜情況不易得出區(qū)間約束公式，需根據(jù)諧波電壓電流實(shí)測值及阻抗計(jì)算值通過式（2）判斷Iu和Ic的大小，其難點(diǎn)在于如何獲得兩側(cè)諧波阻抗幅值及阻抗角的準(zhǔn)確值，其中阻抗角的準(zhǔn)確值更難得到。相比回歸法，獨(dú)立分量法［10］、支持向量機(jī)［11］和極大似然估計(jì)法［16］的諧波阻抗計(jì)算結(jié)果較準(zhǔn)確?？衫们蟪龅淖杩构烙?jì)值，并結(jié)合諧波電壓電流測量信息進(jìn)行諧波溯源。

當(dāng)兩側(cè)諧波阻抗同呈容性時(shí)推導(dǎo)過程與兩側(cè)諧波阻抗同呈感性相似，所有區(qū)間邊界值在180°處具有鏡像對稱性，且判斷條件不變；同理，當(dāng)Zu呈容性、Zc呈感性時(shí)，與Zu呈感性、Zc呈容性的判斷關(guān)系在180°處具有鏡像對稱性。并將系統(tǒng)側(cè)諧波阻抗呈感性時(shí)實(shí)際工程應(yīng)用操作方法總結(jié)至表2中，當(dāng)系統(tǒng)側(cè)諧波阻抗呈容性、測得諧波電壓電流夾角為φ時(shí)，可將其轉(zhuǎn)換為系統(tǒng)側(cè)諧波阻抗呈感性、諧波電壓電流夾角為360°-φ的情況進(jìn)行判斷。

基于區(qū)間約束的諧波溯源方法采用嚴(yán)格不等式制約，其結(jié)果是利用諧波電壓貢獻(xiàn)度法求解的子集（充分條件），說明該方法的溯源結(jié)果一定正確。其研究范圍涵蓋了實(shí)際工程中所有情況，同時(shí)也證明了諧波無功功率方向法在區(qū)間約束下的有效性。

表2 方法歸納表Table 2 Summary table of harmonic source tracing

3 實(shí)例計(jì)算

3.1 算例1

實(shí)測數(shù)據(jù)來自法國一100 MV·A工業(yè)直流電弧爐的150 kV母線，系統(tǒng)短路容量為7500 Mvar。監(jiān)測點(diǎn)系統(tǒng)側(cè)為感性，用戶側(cè)是工業(yè)電弧爐負(fù)荷，為典型的感性負(fù)荷。使用LEM TOPAS 1000電能分析儀測量PCC處諧波電壓和諧波電流，并通過快速傅里葉變換，獲得每分鐘諧波電壓電流數(shù)據(jù)。以3次諧波測量值為例，如圖6所示，圖中為PCC處測得的某日連續(xù)10 h內(nèi)的3次諧波電壓、電流測量值，圖7為測量數(shù)據(jù)散點(diǎn)圖。由圖6、7可知，在電弧爐工作時(shí)，即使諧波電壓和電流的幅值波動(dòng)較大，但其夾角始終在第二象限某范圍內(nèi)波動(dòng)。

圖6 PCC處3次諧波波形圖Fig.6 3rd harmonic waveforms of PCC

圖7 3次諧波實(shí)測數(shù)據(jù)散點(diǎn)圖Fig.7 Scatter plots of measured 3rd harmonic data

實(shí)際電網(wǎng)中，系統(tǒng)和用戶兩側(cè)的諧波發(fā)射水平是實(shí)時(shí)變化的，準(zhǔn)確有效的諧波溯源方法應(yīng)在不引入大量異常數(shù)據(jù)的前提下將監(jiān)測數(shù)據(jù)變化情況考慮在內(nèi)。當(dāng)監(jiān)測數(shù)據(jù)混入異常數(shù)據(jù)后，會使計(jì)算結(jié)果失真。為避免上述情況，利用可信區(qū)間寬度來排除部分較嚴(yán)重的異常數(shù)據(jù)［17-18］。選取95%可信度水平，在確保正常數(shù)據(jù)保留的同時(shí)排除了對結(jié)果誤差起主要作用的異常數(shù)據(jù)，即使剩余少量異常數(shù)據(jù)，其帶來的誤差可被大量正常數(shù)據(jù)削減，不會影響計(jì)算結(jié)果。

根據(jù)可信區(qū)間計(jì)算公式（17），求得諧波電壓和諧波電流的可信區(qū)間。

其中，μ為PCC處諧波電壓和諧波電流的期望值；σ0為測量數(shù)據(jù)的標(biāo)準(zhǔn)差；α=0.05；u=1.96。去除異常邊界值后，重新求得參數(shù)的期望值，如表3所示。

表3 PCC參數(shù)Table 3 Parameters of PCC

將表3參數(shù)代入式（15），可得φ的約束區(qū)間為（103°，167°），從圖 7 可以看出，夾角 φ 主要在 120°附近波動(dòng)，且其可信區(qū)間內(nèi)的數(shù)據(jù)同樣也位于約束區(qū)間內(nèi)，由此可確定用戶側(cè)為主諧波源。

此外，利用本文方法還可以對每組實(shí)測數(shù)據(jù)及電力部門或用戶比較關(guān)注的一些被污染較嚴(yán)重的時(shí)間區(qū)段進(jìn)行諧波溯源分析，找出對PCC處諧波畸變產(chǎn)生主要作用的諧波源。

3.2 算例2

圖8為我國某城市電網(wǎng)局部示意圖，其直流落點(diǎn)及10 kV、35 kV負(fù)荷接入點(diǎn)的監(jiān)測點(diǎn)信息見表4。監(jiān)測點(diǎn)系統(tǒng)側(cè)均為感性，由于城市電網(wǎng)大量使用電纜，以及裝有大容量濾波裝置，使10 kV、35 kV監(jiān)測點(diǎn)的用戶側(cè)呈容性。直流落點(diǎn)的諧波影響，以及其諧波在電網(wǎng)中的交互影響、傳導(dǎo)和放大等問題一直是工程實(shí)際遇到的問題和研究熱點(diǎn)［12，19-20］。

圖8 具有直流落點(diǎn)的局部城市電網(wǎng)示意圖Fig.8 Schematic diagram of partial urban power network with DC coupling point

表4 監(jiān)測點(diǎn)性質(zhì)Table 4 Property of monitoring points

該城市電網(wǎng)中諧波主要來自直流落點(diǎn)和電網(wǎng)中其他非線性負(fù)荷。城市電網(wǎng)中直流落點(diǎn)作為系統(tǒng)側(cè)諧波源從500 kV輸電網(wǎng)向10 kV、35 kV配電網(wǎng)傳播，用戶側(cè)非線性負(fù)荷作為諧波源從10 kV、35 kV配電網(wǎng)向500 kV輸電網(wǎng)傳播。分析城市電網(wǎng)中該10 kV、35 kV配電網(wǎng)及500 kV直流落點(diǎn)諧波的主要來源。

取某日負(fù)荷工作時(shí)10 kV監(jiān)測點(diǎn)的連續(xù)1800 s諧波電壓和電流監(jiān)測數(shù)據(jù)，其11次諧波電壓和電流監(jiān)測數(shù)據(jù)如圖9所示。電能質(zhì)量監(jiān)測平臺1 s記錄一組諧波數(shù)據(jù)，并對該監(jiān)測數(shù)據(jù)做出散點(diǎn)圖，如圖10所示。此時(shí)，諧波電壓電流夾角位于第四象限上半?yún)^(qū)域，可直接判斷系統(tǒng)側(cè)為該站的主諧波源。與10 kV監(jiān)測數(shù)據(jù)的分析方法類似，文獻(xiàn)［1］已經(jīng)指出該城市電網(wǎng)35 kV監(jiān)測點(diǎn)所監(jiān)測的11次諧波的主諧波源同樣位于系統(tǒng)側(cè)。文獻(xiàn)［12］使用Digsilent程序?qū)φ麄€(gè)城市電網(wǎng)進(jìn)行計(jì)算，直流落點(diǎn)諧波電流傳導(dǎo)至各支路的計(jì)算值占測量值的70%以上，可以確定諧波主要是由直流落點(diǎn)諧波在城市電網(wǎng)傳導(dǎo)引起的。

圖9 11次諧波波形圖Fig.9 11th harmonic waveforms

圖10 負(fù)荷工作時(shí)監(jiān)測數(shù)據(jù)散點(diǎn)圖Fig.10 Scatter plots of monitored data when load is working

對500 kV直流落點(diǎn)處監(jiān)測數(shù)據(jù)進(jìn)行分析，如圖11所示?？芍?00 kV處的主諧波源也為系統(tǒng)側(cè)。進(jìn)一步確認(rèn)了該城市電網(wǎng)中各支路諧波主要是由直流落點(diǎn)側(cè)自上而下傳遞。

圖11 諧波監(jiān)測數(shù)據(jù)Fig.11 Monitored harmonic data

4 結(jié)論

a.分析了兩側(cè)諧波阻抗角為任意值時(shí)，諧波電壓貢獻(xiàn)度隨諧波電壓電流夾角的變化規(guī)律，得到基于區(qū)間約束的諧波溯源方法，其結(jié)果是諧波電壓貢獻(xiàn)度直接求解的充分條件，求得的溯源結(jié)果準(zhǔn)確可靠。

b.區(qū)間約束法根據(jù)不同區(qū)間范圍結(jié)合嚴(yán)格不等式制約，多數(shù)區(qū)間得出了不需要阻抗角參數(shù)的諧波溯源方法，部分區(qū)間得出了不需要阻抗幅值和阻抗角僅需諧波電壓電流測量夾角的判斷方法。

c.區(qū)間約束法研究范圍涵蓋了實(shí)際工程中的所有情況。利用多場景實(shí)際工程監(jiān)測數(shù)據(jù)驗(yàn)證了諧波無功功率方向法在區(qū)間約束條件制約下的有效性，并給出實(shí)際工程應(yīng)用的操作方法。

d.基于區(qū)間約束的諧波溯源方法為實(shí)現(xiàn)諧波責(zé)任劃分及諧波經(jīng)濟(jì)管理奠定了基礎(chǔ)，為解決供用電雙方矛盾提供了理論依據(jù)。

參考文獻(xiàn)：

［1］王智琦，楊洪耕，肖楚鵬，等.基于阻抗約束的諧波溯源方法［J］.電力系統(tǒng)自動(dòng)化，2016，40（21）：182-188.WANG Zhiqi，YANG Honggeng，XIAO Chupeng，et al.Harmonic traceability method based on impedance constraints［J］.Automation of Electric Power Systems，2016，40（21）：182-188.

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

［3］HUI J，F(xiàn)REITAS W，VIEIRA J C M，et al.Utility harmonic impedance measurement based on data selection［J］.IEEE Transactions on Power Delivery，2012，27（4）：2193-2202.

［4］張巍，楊洪耕.基于二元線性回歸的諧波發(fā)射水平估計(jì)方法［J］.中國電機(jī)工程學(xué)報(bào)，2004，24（6）：50-53.ZHANG Wei，YANG Honggeng.A method for assessing harmonic emission level based on binary linear regression［J］.Proceedings of the CSEE，2004，24（6）：50-53.

［5］黃舜，徐永海.基于偏最小二乘回歸的系統(tǒng)諧波阻抗與諧波發(fā)射水平的評估方法［J］. 中國電機(jī)工程學(xué)報(bào)，2007，27（1）：93-97.HUANG Shun，XU Yonghai.Assessing harmonic impedance and the harmonic emission levelbased on partialleast-squares regression method［J］.Proceedings of the CSEE，2007，27（1）：93-97.

［6］解紹鋒.基于秩次回歸的系統(tǒng)諧波阻抗分析與諧波評估方法［J］.電力自動(dòng)化設(shè)備，2010，30（11）：29-33.XIE Shaofeng.Harmonic impedance analysis and harmonic assessment based on rank regression［J］.Electric Power Automation Equipment，2010，30（11）：29-33.

［7］孫媛媛，尹志明.基于M估計(jì)穩(wěn)健回歸的多諧波源責(zé)任區(qū)分［J］.中國電機(jī)工程學(xué)報(bào)，2012，32（31）：166-173.SUN Yuanyuan，YIN Zhiming.Quantifyingharmonicresponsibilitiesofmultiple harmonic sourcesbased on M-estimation robust regression［J］.Proceedings of the CSEE，2012，32（31）：166-173.

［8］HUI J，YANG H，LIN S，et al.Assessing utility harmonic impedance based on the covariance characteristic ofrandom vectors［J］.IEEE Transactions on Power Delivery，2010，25（3）：1778-1786.

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

［10］趙熙，楊洪耕.基于快速獨(dú)立分量分析的系統(tǒng)側(cè)諧波阻抗計(jì)算方法［J］. 電力系統(tǒng)自動(dò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.

［11］邱思語，楊洪耕.改進(jìn)的加權(quán)支持向量機(jī)回歸的諧波發(fā)射水平估計(jì)方法［J］. 電工技術(shù)學(xué)報(bào)，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.

［12］陳甜甜，胡蓉，金祖洋，等.鄰近直流輸電落點(diǎn)的城市電網(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.

［13］XU W，LIU X，LIU Y.An investigation on the validity of power-direction method for harmonic source determination ［J］.IEEE Transactions on Power Delivery，2003，18（1）：214-219.

［14］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.

［15］陳海需，楊文煥，李榮高.電纜分布電容引起的港口供配電系統(tǒng)諧波諧振分析［J］. 電力科學(xué)與工程，2014，30（1）：38-41.CHEN Haixu，YANG Wenhuan，LI Ronggao.The analysis of harmonic resonance on distributed capacitance of power cable in harbor power distribution system［J］.Electric Power Science and Engineering，2014，30（1）：38-41.

［16］華回春，賈秀芳，曹東升，等.系統(tǒng)諧波阻抗估計(jì)的極大似然估計(jì)方法［J］. 中國電機(jī)工程學(xué)報(bào)，2014，34（10）：1692-1699.HUA Huichun，JIA Xiufang，CAO Dongsheng，et al.A maximum likelihood method for harmonic impedance estimation［J］.Proceedings of the CSEE，2014，34（10）：1692-1699.

［17］劉興杰，謝春雨.基于貝塔分布的風(fēng)電功率波動(dòng)區(qū)間估計(jì)［J］.

（continued on page 165）（continued from page 160）電力自動(dòng)化設(shè)備，2014，34（12）：26-30.LIU Xingjie，XIE Chunyu.Wind power fluctuation interval estimation based on beta distribution［J］.Electric Power Automation Equipment，2014，34（12）：26-30.

［18］呂曉祿，梁軍，贠志皓，等.風(fēng)電場出力的縱向時(shí)刻概率分布特性［J］. 電力自動(dòng)化設(shè)備，2014，34（5）：40-45.Lü Xiaolu，LIANG Jun，YUN Zhihao，et al.Longitudinal instant probability distribution of wind farm output power［J］.Electric Power Automation Equipment，2014，34（5）：40-45.

［19］楊光亮，邰能靈，鄭曉冬，等.多饋入直流輸電系統(tǒng)諧波交互影響分析［J］. 電力自動(dòng)化設(shè)備，2016，36（1）：105-110.YANG Guangliang，TAI Nengling，ZHENG Xiaodong，et al.Harmonic interaction analysis formulti-infeed HVDC system［J］.Electric Power Automation Equipment，2016，36（1）：105-110.

［20］王燦，羅隆福，陳躍輝，等.一種改進(jìn)型感應(yīng)濾波高壓直流輸電系統(tǒng)及其諧波傳遞特性分析［J］. 電力自動(dòng)化設(shè)備，2015，35（10）：127-132.WANG Can，LUO Longfu，CHEN Yuehui，et al.Harmonic interaction analysis for multi-infeed HVDC system［J］.Electric Power Automation Equipment，2015，35（10）：127-132.