李國棟，李洪安，葛磊蛟，呂金炳，劉亞麗

（1.國網(wǎng)天津市電力公司電力科學(xué)研究院，天津 300384；2.天津大學(xué)電氣與自動化工程學(xué)院，天津 300072）

阻抗矩陣的元素基本不變，因此有

這樣，可以消去最后一項，即得

主動配電網(wǎng)的諧波源定位方法及諧波責(zé)任劃分

李國棟1，李洪安2，葛磊蛟2，呂金炳1，劉亞麗1

（1.國網(wǎng)天津市電力公司電力科學(xué)研究院，天津 300384；2.天津大學(xué)電氣與自動化工程學(xué)院，天津 300072）

針對主動配電網(wǎng)非線性負(fù)荷的諧波特性，提出了一種應(yīng)用于主動配電網(wǎng)的諧波源定位方法及責(zé)任劃分規(guī)則。首先，建立負(fù)荷時域等值模型，利用公共耦合點上電壓、電流的物理關(guān)系，采用參數(shù)辨識法定位主動配電網(wǎng)中的諧波源；其次，通過分析主動配電網(wǎng)中所有諧波源共同作用的母線諧波電壓，并計算該母線與諧波源所在節(jié)點間的諧波互阻抗（自阻抗）而確認(rèn)各個諧波源的貢獻，從而利用諧波貢獻法定量劃分各個諧波源的責(zé)任。最后，通過18節(jié)點主動配電網(wǎng)模型案例，驗證該方法的有效性和合理性。

主動配電網(wǎng)；諧波源；公共耦合點；非線性度；諧波貢獻

隨著配電網(wǎng)中大量分布式電源DG（distributed generation）[1]、微電網(wǎng)的接入，主動配電網(wǎng)ADN（ac?tive distribution network）技術(shù)也應(yīng)運而生。風(fēng)能、太陽能等可再生能源間歇性和隨機性的特性，不僅其發(fā)電量存在不確定性和波動性，而且通過電力電子控制裝置接入電網(wǎng)的可再生能源，將會產(chǎn)生大量諧波；隨著配電網(wǎng)中DG的增多，將導(dǎo)致配電網(wǎng)電能質(zhì)量的進一步惡化[2-3]。這種情況下，“獎懲性措施”的實施更加困難。為了給諧波治理和實施“獎懲性措施”提供可靠依據(jù)[18-20]，ADN中諧波源定位和諧波責(zé)任定量劃分成為電網(wǎng)公司一項重要的工作內(nèi)容，是一個值得深入研究的課題。

目前諧波定位的方法主要從兩個層面出發(fā)：基于功率潮流的檢測方法和基于諧波阻抗的檢測方法?；诠β食绷鞯臋z測方法主要有有功功率方向法[4-5]、無功功率方向法[6]以及無功功率變化法[7]等，但是，由于某些諧波源能吸收諧波功率，這些方法的準(zhǔn)確性一般難以保證；基于諧波阻抗的檢測方法是從測量用戶和系統(tǒng)側(cè)的諧波阻抗入手，非線性元件導(dǎo)致電壓和電流波形畸變，從而出現(xiàn)高次諧波的原理，基于此提出了參數(shù)辨識法[8]、最小二乘系統(tǒng)辨識法[9]等方法。由于ADN中含有大量DG和非線性負(fù)荷，功率具有雙向流動的特點，基于諧波功率的辨識法可能失效，因此本文選用不考慮背景諧波的參數(shù)辨識法。

關(guān)于諧波責(zé)任劃分方面，國內(nèi)外學(xué)者的研究都立足于用戶和供電方的共同責(zé)任，旨在分清公共耦合點PCC（point of common coupling）兩側(cè)的諧波責(zé)任[10]，但是大量事實證明，配電網(wǎng)中存在不止一個非線性負(fù)荷，因此每個節(jié)點的電壓波形畸變并不僅僅由一個非線性負(fù)荷引起，而是由配電網(wǎng)中所有非線性負(fù)荷共同作用的結(jié)果[11]。當(dāng)前，國內(nèi)外的學(xué)者在多個諧波源的諧波責(zé)任定量劃分方面一般以諧波源注入的諧波電流各自在PCC點所引起的諧波電壓降作為劃分依據(jù)。

本文首先建立了負(fù)荷時域等值模型，并利用公共耦合點上電壓、電流的物理關(guān)系，采用參數(shù)辨識法定位主動配電網(wǎng)中的諧波源；然后，通過分析主動配電網(wǎng)中所有諧波源共同作用的母線諧波電壓，并計算該母線與諧波源所在節(jié)點間的諧波互阻抗（自阻抗），確認(rèn)各個諧波源的貢獻，從而利用諧波貢獻法定量劃分各個諧波源的責(zé)任；最后，通過18節(jié)點主動配電網(wǎng)模型案例，驗證該方法的有效性和合理性。

1 諧波源定位原理

1.1 諧波產(chǎn)生原理

非線性元件在工頻電壓下，除了會產(chǎn)生工頻電流外，還會產(chǎn)生工頻整倍次數(shù)的電流，即諧波污染[12]，如圖1所示。圖1（a）為一個簡單的配電網(wǎng)示意，Z1為線性負(fù)荷，Z2為非線性負(fù)荷，i1（t）、i2（t）分別為其支路上電流，u（t）為PCC處電壓。由圖（b）、（c）、（d）可見，工頻電壓加在線性負(fù)荷與非線性負(fù)荷所產(chǎn)生電流的不同。由此可見，非線性負(fù)荷是產(chǎn)生諧波的主要原因。

1.2 參數(shù)辨識法原理

ADN中含有大量的非線性元件，諧波問題非常突出，為了有效解決這些問題，有效實施“獎懲性措施”，進行諧波源定位非常重要。本文將諧波源的定位問題轉(zhuǎn)換為負(fù)荷參數(shù)辨識的問題。以阻感負(fù)荷為例，其等值模型如圖2所示。

圖1 線性負(fù)荷與非線性負(fù)荷對比Fig.1 Comparison between linear load and nonlinear load

圖2 負(fù)荷參數(shù)模型Fig.2 Load parameter model

圖中，Z為阻感負(fù)載，u（t）和i（t）為PCC的測量電壓和測量電流。若阻抗參數(shù)呈線性關(guān)系，則有

式中：R(t)、L(t)分別為t時刻感性負(fù)荷的電阻和電感部分等值參數(shù)；uL(t)為t時刻感性負(fù)荷的電壓。

當(dāng)負(fù)荷為非線性負(fù)荷時，按照微元法的思想，采用分段線性化法計算，則在整個采樣周期內(nèi)其阻抗參數(shù)呈非線性關(guān)系。

采樣頻率遠大于負(fù)荷參數(shù)的變化頻率，在一個微元時間Δt內(nèi)，可認(rèn)為相鄰的4個采樣點上的R和L不變。根據(jù)梯形等效法則，設(shè)Δt=t3-t1，則有

矩陣形式為

通過對公共耦合點的電壓與電流的測量，利用式（5）可計算出相應(yīng)的R（t）和L（t）。若兩者不呈線性，則直接可以判定為諧波源。

對于阻容負(fù)載，可進行類似分析。需要指出的一點是，非線性計算出的R、L和C，并不完全代表負(fù)荷的實際參數(shù)，而是體現(xiàn)出負(fù)荷電壓和電流間的數(shù)學(xué)關(guān)系。由于實際負(fù)荷大多呈感性，一般采用式（5）計算負(fù)荷的等值阻抗參數(shù)。

ADN中，可以通過測量負(fù)荷或DG的電壓電流，利用式（5），獲得多組負(fù)荷或DG的電阻電感值，繪制電阻電感的時域變化曲線，根據(jù)曲線形狀辨識諧波源。

2 諧波貢獻法

上述利用參數(shù)辨識法獲得了ADN中的諧波源，為了有效進行諧波責(zé)任劃分，利用諧波貢獻法進行定量分析，其原理如下。

ADN的節(jié)點阻抗矩陣體現(xiàn)諧波源對各母線諧波電壓的影響，有

式中：Zh為h次諧波對應(yīng)的阻抗矩陣，下文簡化為Z。

若關(guān)注母線m，則節(jié)點阻抗矩陣的第m行的h次諧波電壓為ADN中各個諧波源的貢獻之和，即

式中：Zmj為節(jié)點m與節(jié)點j之間的諧波互阻抗（m≠j）或自阻抗（m=j）；Umj為節(jié)點j諧波源對母線m的諧波貢獻電壓。用Um的相角作為參考相角，可以區(qū)分各個諧波源對母線m的h次諧波電壓起助增還是抑制作用。

電力負(fù)荷的負(fù)荷曲線一般包含緩慢變化分量和快速變化分量[13-14]兩部分，其中，快速變化分量之間是相互獨立的。事實證明，諧波源的諧波電流發(fā)射量也是相互獨立的[12]。因此可以把阻抗矩陣方程線性地劃分為緩慢變化分量方程和快速變化分量方程兩部分，即

根據(jù)快速變化分量的相互獨立特性，提取節(jié)點i對應(yīng)的諧波源在母線m上的諧波貢獻。式（9）的第m行可以表示為

根據(jù)前文可知，關(guān)鍵是把ZmiIfasti也就是Umj分離出來，

即去除式（11）的最右邊項，即對式（11）左右同時求期望值，可得

由于ADN的阻抗矩陣元素和各個諧波源注入諧波電流是相互獨立的，因此式（12）可以改寫為

用式（11）減去式（13），可得

阻抗矩陣的元素基本不變，因此有

這樣，可以消去最后一項，即得

由此，可以獲得母線m與節(jié)點i間的諧波互阻抗，即

因此，節(jié)點i的諧波源對母線m的諧波貢獻百分比為

綜上所述，諧波貢獻法是通過對關(guān)注母線諧波電壓和諧波源節(jié)點的諧波電流進行采樣分析，利用式（19）獲得兩者之間的諧波互（自）阻抗，進而計算各個諧波源對關(guān)注母線的諧波貢獻百分比，定量劃分諧波責(zé)任。

3 算例

為了驗證上述方法的有效性，本文搭建了18節(jié)點的10 kV ADN模型[16]，如圖3所示。在節(jié)點3、6和10處分別接入DG，采用DG模型為Detail模型[17]，該模型主體為直流逆變器，直接體現(xiàn)以電力電子元件為主體的逆變器并網(wǎng)DG，忽略DG自身功率的波動性和隨機性及其對直流側(cè)電壓的影響，認(rèn)為直流側(cè)電壓已控制恒定，故在直流側(cè)選用恒定電壓源代替并網(wǎng)逆變器模型。

圖3 18節(jié)點ADN模型Fig.3 18-node ADN model

首先利用參數(shù)辨識法定位諧波源，提取各個節(jié)點在一個工頻周期內(nèi)的電壓和電流，采樣時間為2×10-4s。根據(jù)式（5），可以獲得一個周期內(nèi)的電阻電感變化情況，并據(jù)此繪制電阻電感的時域變化曲線，根據(jù)這些曲線可以確認(rèn)諧波源所在節(jié)點為節(jié)點3、6和10，即DG所在節(jié)點。節(jié)點3、6和10的電阻電感的時域變化曲線如圖4～圖6所示。

從圖4～圖6可知，DG產(chǎn)生諧波的原理與非線性負(fù)荷相似，跟DG本身阻抗在運行中會發(fā)生波動有關(guān)。但DG產(chǎn)生的波動除去DG自身的間歇性和不穩(wěn)定性外，其接入電網(wǎng)所需的電力電子調(diào)控器件也是DG非線性度增加的主要原因。

利用參數(shù)辨識法定位諧波源后，再利用諧波貢獻法確認(rèn)各個諧波源對關(guān)注母線的諧波責(zé)任。選取關(guān)注母線為節(jié)點6所在母線，以5次諧波為例。模型中每0.02 s取一個采樣點，采集100個樣本。然后根據(jù)式（19），可以獲得關(guān)注母線與3個諧波源節(jié)點間的諧波互阻抗，如表1所示。

根據(jù)所獲得的諧波互阻抗、各個諧波源的5次諧波電流以及式（20），可以計算獲得各個諧波源對關(guān)注母線5次諧波的貢獻百分比，如圖7所示。由圖可見，節(jié)點3處諧波源對關(guān)注母線的諧波貢獻百分比為73.26%，節(jié)點6處為66.78%，節(jié)點10處為80.34%。

圖4 節(jié)點3的電阻電感Fig.4 Resistance and inductance at node 3

圖5 節(jié)點6的電阻電感Fig.5 Resistance and inductance at node 6

由圖7和表1可知，雖然節(jié)點3處諧波源與關(guān)注母線間諧波互阻抗大于節(jié)點10，但由于節(jié)點10處諧波電流較大，其諧波貢獻大于節(jié)點3。由這一比較結(jié)果可知通過改變系統(tǒng)結(jié)構(gòu)，可以降低某些諧波源的諧波貢獻，從而降低ADN中特定母線上的諧波含量。

圖6 節(jié)點10的電阻電感Fig.6 Resistance and inductance at node 10

表1 關(guān)注母線與諧波源間的諧波互阻抗Tab.1 Harmonic mutual impedance between harmonic source buses and concerned buses

圖7 各諧波源諧波貢獻Fig.7 Contribution of all harmonic sources

4 結(jié)論

（1）本文利用測量負(fù)荷瞬時電流電壓值的參數(shù)辨識方法，能有效辨識ADN中的諧波源，具有原理簡單、不受背景諧波影響等優(yōu)點。

（2）一般單純從電壓和電流的波形畸變率，很難準(zhǔn)確判斷諧波源，必須從分析負(fù)荷特性入手，且等值阻抗參數(shù)R、L、C組成的非線性負(fù)荷是造成公共連接點電壓波形畸變的根本原因。

（3）本文考慮到諧波源發(fā)射諧波的隨機性，利用獨立隨機矢量的統(tǒng)計特性，提取了各個諧波源對關(guān)注母線的諧波貢獻，獲得各個諧波源對系統(tǒng)的助增或抑制作用大小，能夠在不干擾ADN正常運行的前提下定量劃分各個諧波源的諧波責(zé)任。

[1]Djapic P，Ramsay C，Pudjianto D，et al.Taking an active ap?proach[J].IEEE Power and Energy Magazine，2007，5（4）：68-77.

[2]張勇（Zhang Yong）.分布式發(fā)電對電網(wǎng)繼電保護的影響綜述（Review of impact of distributed generation on distribu?tion system protection）[J].電力系統(tǒng)及其自動化學(xué)報（Pro?ceedings of the CSU-EPSA），2010，22（2）：145-151.

[3]于建成，遲福建，徐科，等（Yu Jiancheng，Chi Fujian，Xu Ke，et al）.分布式電源接入對電網(wǎng)的影響分析（Analysis of the impact of distributed generation on power grid）[J].電力系統(tǒng)及其自動化學(xué)報（Proceedings of the CSU-EPSA），2012，24（1）：138-141.

[4]Yang H，Porotte P，Robert A，et al.Assessing the harmonic emission level from one particular customer[C]//3rd Inter?national Conference on Power Quality.Amsterdam，Nether?lands，1994：160-166.

[5]艾永樂，鄭建云.（Ai Yongle，Zheng Jianyun）.基于諧波有功功率貢獻量的主諧波源定位（Localization of the main harmonic source based on the harmonic active power contri?butions）[J].電力系統(tǒng)保護與控制（Power System Protec?tion and Control），2015，43（7）：16-22.

[6]Li Chun，Xu Wilsun，Tayjasanant T.A“critical impedance”-based method for identifying harmonic sources[J].IEEE Trans on Power Delivery，2004，19（2）：671-678.

[7]Hamzah N，Mohamed A，Hussain A.Methods for determin?ing utility and customer harmonic contributions at the point of common coupling[C]//National Power and Energy Confer?ence.Bangi，Malaysia，2003：167-171.

[8]蔡明，唐昆明，張?zhí)冢–ai Ming，Tang Kunming，Zhang Tai?qin）.基于參數(shù)辨識方法的諧波源定位（Harmonic source localization based on parameter identification）[J].電網(wǎng)技術(shù)（Power System Technology），2011，35（6）：134-138.

[9]趙勇，沈紅，李建華，等（Zhao Yong，Shen Hong，Li Jianhua，et al）.諧波源的識別及其與非諧波源的分離方法（Ap?proach of identification and separation of harmonic sources）[J].中國電機工程學(xué)報（Proceedings of the CSEE），2002，22（5）：84-87.

[10]楊洪耕，王磊（Yang Honggeng，Wang Lei）.基于拉蓋爾多項式的非線性負(fù)荷諧波發(fā)射水平估計（Assessing harmon?ic emission levels from non-loads using Laguerre polynomi?als）[J].中國電機工程學(xué)報（Proceedings of the CSEE），2005，25（7）：81-85．

[11]馬歷，劉開培，雷肖（Ma Li，Liu Kaipei，Lei Xiao）.配電網(wǎng)諧波源定位的支持向量機估計算法（Support vector ma?chine estimating algorithm of harmonic source localizationin power distribution network）[J].中國電機工程學(xué)報（Pro?ceedings of the CSEE），2008，28（10）：111-116．

[12]林海雪，范明天，薛惠，等.電力系統(tǒng)諧波[M].北京：中國電力出版社，2008．

[13]Gursoy E，Niebur D.Harmonic load identification using complex independent component analysis[J].IEEE Trans on Power Delivery，2009，24（1）：285-292．

[14]Liao Huaiwei，Niebur D.Load profile estimation in electric transmission networks using independent component analy?sis[J].IEEE Trans on Power Systems，2003，18（2）：707-715．

[15]Hui J，Yang H，Lin S，et al.Accessing utility harmonic im?pedance based on the covariance characteristic of random vectors[J].IEEE Trans on Power Delivery，2010，25（3）：1778-1786．

[16]Lee San-Yi，Wu Chi-Jui.Combined compensation structure of an SVC and an active filter for unbalanced three phase distribution feeders with harmonic distortion[C]//Fourth In?ternational Conference on Advances in Power System Con?trol，Operation and Management.Hong Kong，China，1997：543-548.

[17]周文珊（Zhou Wenshan）.主動配電網(wǎng)建模及諧波特性研究（Modeling and Harmonic Characteristics Research of Ac?tive Distribution Network）[D].北京：北京交通大學(xué)電氣工程學(xué)院（Beijing：School of Electrical Engineering，Beijing Jiaotong University），2012.

[18]葛少云，董佳霖，劉福潮，等（Ge Shaoyun，Dong Jialin，Liu Fuchao，et al）.配電網(wǎng)諧波治理降損效果實證方法（Em?pirical method of harmonic control loss reduction effect in distribution network）[J].電力系統(tǒng)及其自動化學(xué)報（Pro?ceedings of the CSU-EPSA），2015，27（3）：1-4.

[19]喬和，劉陽，單錦寧，等（Qiao He，Liu Yang，Shan Jinning，et al）.改進的變步長自適應(yīng)諧波檢測算法（Adaptive har?monic detecting algorithm based on improved variable stepsize）[J].電力系統(tǒng)及其自動化學(xué)報（Proceedings of the CSU-EPSA），2015，27（5）：23-27.

[20]李國華，張影，孫強，等（Li Guohua，Zhang Ying，Sun Qiang，et al）.有源電力濾波器選擇性諧波補償方法（Se?lective harmonic compensation method for active power fil?ter）[J].電力系統(tǒng)及其自動化學(xué)報（Proceedings of the CSU-EPSA），2015，27（6）：32-36.

Harmonic Source Location Method and Division of Harmonic Responsibility in Active Distribution Network

LI Guodong1，LI Hong’an2，GE Leijiao2，LYU Jinbing1，LIU Yali1
（1.Tianjin Electrical Power Research Institute，State Grid Tianjin Electric Power Company，Tianjin 300384，China；2.School of Electrical Engineering and Automation，Tianjin University，Tianjin 300072，China）

Due to the feature of nonlinear load in active distribution network，a harmonic source location method and di?vision of harmonic responsibility in active distribution network are put forward in this paper.Firstly，an equivalent load model in time domain is built.Using the physical relationship between voltage and current at the point of common cou?pling（PCC），the harmonic source in active distribution network is located using parameter identification method.Sec?ondly，it analyzes the harmonic voltages on the same bus in active distribution network，calculates accurate harmonic mutual impedance between harmonic source buses and the concerned buses，obtains the harmonic contribution for all harmonic sources，and determines the responsibility for all harmonic sources using quantitative analysis method.Final?ly，an 18-node active distribution network is used to verify that the proposed method is reasonable and effective.

active distribution network（ADN）；harmonic source；point common coupling；nonlinearity degree；har?monic contribution

TM72

A

1003-8930（2016）11-0095-05

10.3969/j.issn.1003-8930.2016.11.016

2015-04-16；

2016-03-17

國網(wǎng)天津市電力公司科技資助項目（KJ15-1-10）

李國棟（1978—），男，本科，高級工程師，研究方向為電能質(zhì)量和新能源。Email：tjlgd@163.com

李洪安（1990—），男，碩士研究生，研究方向為主動配電網(wǎng)諧波源定位與諧波責(zé)任定量劃分。Email：lihongantj@163.com

葛磊蛟（1984—），男，通信作者，博士，講師，研究方向為智能配用電與不確定性仿真。Email：legendglj99@tju.edu.cn