基于電流注入法的光伏電源并網(wǎng)諧波電壓畸變分析

石磊磊，韓天華，唐 超，衛(wèi) 丹，趙凡凡，賈清泉，?

（1.國(guó)網(wǎng)河北省電力公司邢臺(tái)供電分公司，河北邢臺(tái)054001；2.燕山大學(xué)電氣工程學(xué)院，河北秦皇島066004）

基于電流注入法的光伏電源并網(wǎng)諧波電壓畸變分析

石磊磊1，2，韓天華1，唐 超1，衛(wèi) 丹1，趙凡凡2，賈清泉2，?

（1.國(guó)網(wǎng)河北省電力公司邢臺(tái)供電分公司，河北邢臺(tái)054001；2.燕山大學(xué)電氣工程學(xué)院，河北秦皇島066004）

通過研究光伏并網(wǎng)對(duì)電網(wǎng)電壓畸變的影響機(jī)理，本文提出了基于電流注入法的電網(wǎng)諧波電壓畸變分析方法。建立了光伏并網(wǎng)系統(tǒng)的諧波源模型，將光伏電源等效為諧波電流注入源，推導(dǎo)出光伏諧波注入源與配電網(wǎng)各點(diǎn)諧波電壓畸變率的關(guān)系式，得出光伏電源接入位置和有功出力對(duì)電網(wǎng)諧波電壓的影響關(guān)系。對(duì)含多個(gè)光伏電源接入的配電網(wǎng)諧波分布特性和相互影響關(guān)系進(jìn)行了理論分析與仿真，得出光伏電源接入后的諧波分布規(guī)律，并驗(yàn)證了本文方法的正確性。

光伏發(fā)電；電能質(zhì)量；諧波畸變；電流注入法

0 引言

近年來，隨著國(guó)家對(duì)光伏發(fā)電重視程度和政策支持力度的加大，光伏并網(wǎng)發(fā)電的規(guī)模越來越大。2015年3月，國(guó)家能源局下發(fā)《關(guān)于下達(dá)2015年光伏發(fā)電建設(shè)實(shí)施方案的通知》，2015年全國(guó)將新增光伏電站建設(shè)規(guī)模為17.8 GW。光伏電源大多通過逆變器接入配電網(wǎng)，由于電力電子設(shè)備的非線性，接入配電網(wǎng)會(huì)產(chǎn)生諧波問題。當(dāng)多個(gè)光伏電源并網(wǎng)時(shí)，不同光伏逆變器會(huì)帶來諧波交互影響［1］，進(jìn)一步提高了光伏發(fā)電系統(tǒng)注入電網(wǎng)的諧波含量，嚴(yán)重影響電網(wǎng)的供電質(zhì)量［2?3］。因此分析光伏電源的并網(wǎng)位置和輸出功率對(duì)電網(wǎng)諧波的影響具有重要理論和實(shí)際意義。

國(guó)內(nèi)外學(xué)者針對(duì)光伏電源接入配電網(wǎng)產(chǎn)生的電能質(zhì)量問題開展了一定的研究。文獻(xiàn)［4］利用解耦潮流算法對(duì)含分布式電源的配電網(wǎng)進(jìn)行了諧波分析，表明分布式電源接入配電網(wǎng)會(huì)注入大量的諧波電流，增加系統(tǒng)的諧波畸變水平。文獻(xiàn)［5?6］分析了分布式電源不同接入技術(shù)對(duì)電網(wǎng)諧波分布的影響。文獻(xiàn)［7］利用MATLAB／Simulink對(duì)分布式電源產(chǎn)生的諧波在電網(wǎng)中的傳播機(jī)理進(jìn)行了分析，表明分布式電源接入電網(wǎng)的電氣距離越遠(yuǎn)，高次諧波電流越小。文獻(xiàn)［8?10］討論了光伏并網(wǎng)逆變器的輸出諧波模型。

本文建立了光伏并網(wǎng)系統(tǒng)的諧波源模型和接入光伏電源的配電網(wǎng)等效模型，將光伏電源等效為諧波電流注入源，通過理論推導(dǎo)得出光伏電源不同接入位置和有功出力與電網(wǎng)諧波電壓的關(guān)系式。對(duì)含光伏電源的配電網(wǎng)進(jìn)行理論和仿真分析，得出各節(jié)點(diǎn)各次諧波電壓畸變率和相互影響規(guī)律，驗(yàn)證了本文分析方法的正確性。

1 光伏并網(wǎng)系統(tǒng)諧波源建模

光伏發(fā)電系統(tǒng)主要的諧波來源是電力電子逆變器，其輸出的諧波含量取決于自身特性和工作條件，因此可以將其看成內(nèi)阻很大的電流源。光伏并網(wǎng)逆變器的諧波特性可以用式（1）來表示。

本文采用電流注入法，即假定由某些諧波電流源注入節(jié)點(diǎn)的各次諧波電流幅值和相角，都僅僅與流過該電流源的基波電流的幅值和相角成線性關(guān)系。由此可以得出諧波電流源接入配電網(wǎng)的諧波分布規(guī)律。

光伏電源并入配電網(wǎng)如圖1所示。假設(shè)電網(wǎng)中三相負(fù)荷平衡，忽略線路集膚效應(yīng)和分布電容。光伏電源在距離B母線a處接入，設(shè)L為B母線到負(fù)荷的距離，U0為B母線端電壓，Ug為光伏電源接入點(diǎn)的電壓，U為輸電線路末端電壓。

圖1 接入光伏電源的配電網(wǎng)模型Fig.1 The distribution network with photovoltaic connected

光伏電源并網(wǎng)的諧波網(wǎng)絡(luò)模型如圖2所示。輸電線路第h次諧波的單位長(zhǎng)度阻抗為r＋jhx，光伏電源和負(fù)荷諧波源向電網(wǎng)注入的諧波用電流源gh和h表示，其幅值與基波電流幅值呈線性關(guān)系，線性系數(shù)分別為kgh和kh，如式（2）、（3）所示：

式中，Pg和cosφg分別為光伏電源輸出的有功功率和功率因數(shù)，PL和cosφL分別為負(fù)荷的有功功率和功率因數(shù)。

圖2 接入光伏電源配電網(wǎng)諧波模型Fig.2 Harmonic model of the distribution network with photovoltaic connected

2 光伏電源接入配電網(wǎng)諧波畸變率分析

2.1 光伏電源接入容量對(duì)線路諧波電壓的影響

圖1中，光伏電源未接入時(shí)，距離B母線a處 h次諧波電壓幅值為

接入光伏電源后，第h次（h≥2）諧波電壓幅值為

光伏電源接入后諧波電壓的變化量ΔUh為

將式（6）代入式（3）可得

由式（7）可知，諧波電壓變化量受光伏電源容量與接入位置的影響。在光伏接入位置不變情況下，如果光伏電源接入容量Pg增大，由式（7）可知，諧波電壓增大，且增加量ΔUh與Pg成正比關(guān)系。

2.2 光伏電源接入位置對(duì)諧波電壓的影響

1）光伏電源接入點(diǎn)上游線路的諧波電壓

光伏電源接入點(diǎn)上游線路距離B母線L（L＜a）處的h次諧波電壓為

在光伏電源接入容量Pg和負(fù)荷PL不變情況下，由式（2）和式（3）可知，諧波電流Igh和Ih也不變。根據(jù)式（8），隨著L增大，光伏電源接入點(diǎn)上游線路諧波電壓也隨之增大。

2）光伏電源接入點(diǎn)下游線路的諧波電壓

光伏電源接入點(diǎn)下游線路距離B母線a處的h次諧波電壓為

由式（9）可知，對(duì)于光伏電源接入點(diǎn)下游線路，隨著L增大，諧波電壓也隨之增大。

當(dāng)多個(gè)光伏電源共同接入配電網(wǎng)時(shí)，諧波電壓的變化量為光伏電源單獨(dú)接入時(shí)產(chǎn)生諧波電壓的疊加，第k個(gè)節(jié)點(diǎn)處h次諧波電壓為

式中，ai為第i個(gè)光伏電源接入點(diǎn)與B母線之間的距離，N為配電網(wǎng)接入光伏電源的數(shù)量。

2.3 配電網(wǎng)接入光伏電源后電壓諧波畸變率的計(jì)算方法

本文中光伏電源接入配電網(wǎng)后各節(jié)點(diǎn)電壓諧波畸變率的計(jì)算方法如下：

1）將某時(shí)刻配電網(wǎng)各節(jié)點(diǎn)光伏電源接入容量（未接入光伏的節(jié)點(diǎn)，其光伏電源容量為0）、功率因數(shù)、線路節(jié)點(diǎn)間長(zhǎng)度、單位長(zhǎng)度電阻、單位長(zhǎng)度電抗等參數(shù)代入式（10）求出節(jié)點(diǎn)k的諧波電壓。

2）根據(jù)式（11）計(jì)算各節(jié)點(diǎn)電壓諧波畸變率：

式中，Uh為第h次諧波電壓，U1為基波電壓。

3）根據(jù)步驟（2）的結(jié)果，判斷分布式光伏電源接入后電網(wǎng)各節(jié)點(diǎn)的諧波畸變率是否滿足國(guó)家標(biāo)準(zhǔn)的要求，進(jìn)行匯總可得電壓諧波畸變率越限的節(jié)點(diǎn)。

3 算例分析

算例選取某實(shí)際配電網(wǎng)的一條出線，沿饋線將每一個(gè)集中負(fù)荷視為一個(gè)節(jié)點(diǎn)，電網(wǎng)母線端編號(hào)為1，依次編號(hào)，構(gòu)成一條10節(jié)點(diǎn)的網(wǎng)絡(luò)，如圖3所示。根據(jù)2.3節(jié)所述方法計(jì)算各節(jié)點(diǎn)電壓畸變率的理論值，并將其與MATLAB／Simulink仿真得到的結(jié)果進(jìn)行比較。

3.1 光伏電源接入不同節(jié)點(diǎn)時(shí)諧波電壓畸變率分析

1）接入光伏電源前后各節(jié)點(diǎn)諧波電壓畸變率對(duì)比

在圖3中的節(jié)點(diǎn)2加入容量為30 kW的光伏電源，計(jì)算光伏電源接入前后，各節(jié)點(diǎn)的電壓畸變率如圖4所示。從圖中可知，電壓畸變率變化量ΔTHDV隨著光伏電源接入配電網(wǎng)距離的增加而增大。光伏電源接入前的諧波電壓畸變率遠(yuǎn)小于光伏電源接入后的諧波電壓畸變率，說明光伏電源接入配電網(wǎng)后，極大地增加了諧波電壓畸變率，成為配電網(wǎng)中重要的諧波污染源。圖4中各節(jié)點(diǎn)諧波電壓畸變率理論計(jì)算值與仿真值之間最大相對(duì)誤差約為4.2%，主要原因是本文方法忽略了電網(wǎng)諧波與光伏電源注入諧波之間的相互作用，因此造成理論值與仿真值之間的誤差，但可以滿足光伏電源接入配電網(wǎng)電壓諧波畸變率評(píng)估計(jì)算的需要。

2）單個(gè)光伏電源接入不同節(jié)點(diǎn)時(shí)諧波電壓畸變率

在圖3中4、10節(jié)點(diǎn)處依次加入有功功率為30 kW的單一光伏電源，仿真得到饋線各節(jié)點(diǎn)的諧波電壓畸變率THDV曲線如圖5所示。

圖3 含光伏電源的配電網(wǎng)結(jié)構(gòu)圖Fig.3 Structure of distribution network with photovoltaic power connected

圖4 加入光伏電源前后諧波電壓畸變率Fig.4 Harmonic voltage distortion rate before and after the connection of photovoltaic power

圖5 單個(gè)光伏電源接入不同節(jié)點(diǎn)時(shí)諧波電壓畸變率Fig.5 The THDVof single photovoltaic power connected with different nodes

如圖5所示，相同容量光伏電源接入電網(wǎng)的不同位置，各節(jié)點(diǎn)的諧波電壓畸變率差別較大。光伏電源接入點(diǎn)上游各節(jié)點(diǎn)，諧波電壓畸變率較??；光伏電源接入點(diǎn)下游各節(jié)點(diǎn)，諧波電壓畸變率較大，其中某些節(jié)點(diǎn)的諧波電壓畸變率可能會(huì)超過《電能質(zhì)量 公用電網(wǎng)諧波》（GB／T 14549—93）規(guī)定的限值。圖5中各節(jié)點(diǎn)諧波電壓畸變率理論值與仿真值最大相對(duì)誤差約為4.6%。

3）多個(gè)光伏電源接入時(shí)諧波電壓畸變率

在圖3中節(jié)點(diǎn)4和10同時(shí)接入兩個(gè)有功出力均為30 kW的光伏電源，與單獨(dú)接入時(shí)得到各次諧波的電壓畸變率進(jìn)行對(duì)比。圖6和圖7分別為節(jié)點(diǎn)2各次諧波畸變率的理論值和仿真值，圖8和圖9分別為節(jié)點(diǎn)4各次諧波畸變率的理論值和仿真值。由于實(shí)際光伏電源輸出的偶次諧波和19次以上的奇次諧波含量很少，因此這些諧波忽略不計(jì)。

圖6 節(jié)點(diǎn)4各次諧波THDV理論值Fig.6 THDVtheoretical value of node 4

圖7 節(jié)點(diǎn)4各次諧波THDV仿真值Fig.7 THDVsimulation value of node 4

圖8 節(jié)點(diǎn)10各次諧波THDV理論值Fig.8 THDVtheoretical value of node 10

圖9 節(jié)點(diǎn)10各次諧波THDV仿真值Fig.9 THDVsimulation value of node 10

光伏電源并網(wǎng)后產(chǎn)生的諧波在傳播過程中以5次、9次和11次諧波為主。當(dāng)兩個(gè)光伏電源同時(shí)接入配電網(wǎng)時(shí)，二者產(chǎn)生的各次諧波會(huì)互相作用，5次、9次和11次等低次諧波變化較明顯。通過對(duì)比可以發(fā)現(xiàn)，各次諧波電壓畸變率理論計(jì)算值與仿真值最大相對(duì)誤差約為4.5%。

3.2 光伏電源輸出不同有功時(shí)諧波電壓畸變率分析

配電網(wǎng)接入負(fù)荷總量為80 kW，在節(jié)點(diǎn)3處分別接入有功功率為50 kW和90 kW的光伏電源，仿真得到諧波電壓畸變率THDV曲線如圖10所示。從圖中可以看出，當(dāng)光伏電源接入電網(wǎng)的位置不變時(shí)，光伏電源輸出有功功率為90 kW時(shí)（超過了配電網(wǎng)接入負(fù)荷總量），諧波電壓畸變率THDV較高，有些節(jié)點(diǎn)（節(jié)點(diǎn)3～10）會(huì)超出國(guó)家標(biāo)準(zhǔn)規(guī)定的限值。通過比較可知，理論計(jì)算值與仿真值最大相對(duì)誤差約為4.8%。

圖10 光伏電源輸出不同有功時(shí)諧波電壓畸變率Fig.10 THDVof different photovoltaic active power output

4 結(jié)論

本文運(yùn)用電流注入法建立了光伏電源接入配電網(wǎng)的諧波電壓畸變率分析方法，得出光伏電源接入位置和有功出力對(duì)電網(wǎng)諧波電壓影響關(guān)系的計(jì)算式。對(duì)某配電網(wǎng)接入光伏電源前后配電網(wǎng)中諧波電壓畸變率變化情況進(jìn)行了理論分析和仿真。理論計(jì)算和仿真結(jié)果之間誤差較小，驗(yàn)證了所提出方法的有效性。本文的主要結(jié)論如下：

1）有功出力相同的光伏電源接入配電網(wǎng)的位置越接近配電網(wǎng)末端，接入點(diǎn)下游饋線各節(jié)點(diǎn)的電壓畸變?cè)絿?yán)重。反之，越接近配電網(wǎng)首端，造成的諧波影響越小。因此，光伏電源不適合接入配電網(wǎng)饋線末端。

2）光伏電源接入配電網(wǎng)位置固定時(shí)，饋線各負(fù)荷節(jié)點(diǎn)上的諧波電壓畸變率由光伏電源的有功出力決定，有功出力越大，饋線各負(fù)荷節(jié)點(diǎn)電壓諧波畸變率就越大，畸變率增加量與光伏電源的有功出力成正比，當(dāng)光伏電源有功出力超過負(fù)荷時(shí)會(huì)超過國(guó)家標(biāo)準(zhǔn)允許的限值。

［1］宣文華，牛益國(guó)，謝小英，等.基于頻域分析和低通濾波的光伏并網(wǎng)逆變器諧振抑制研究［J］.燕山大學(xué)學(xué)報(bào)，2015，39（4）：347?351.

［2］丁明，王偉勝，王秀麗，等.大規(guī)模光伏發(fā)電對(duì)電力系統(tǒng)影響綜述［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2014，34（1）：1?14.

［3］趙爭(zhēng)鳴，雷一，賀凡波，等.大容量并網(wǎng)光伏電站技術(shù)綜述［J］.電力系統(tǒng)自動(dòng)化，2011，35（12）：101?107.

［4］Ladjavardi M，Masoum M A S，Islam S.Impact of a SG nonliner model on the harmonic distortion of a distribution generation system［C］／／2008 Australasian Universities Power Engineering Conference，Sydney，Australia，2008：226?231.

［5］鐘清，高新華，余南華，等.諧波約束下的主動(dòng)配電網(wǎng)分布式電源準(zhǔn)入容量與接入方式［J］.電力系統(tǒng)自動(dòng)化，2014，38（24）：108?113.

［6］Fortes M Z，F(xiàn)erreira V H，Machado I S，et al.Harmonic analysis of distributed generation in smart city Búzios project［C］／／2015 IEEE Workshop on Power Electronics and Power Quality Applica?tions（PEPQA），Bogotá，Colombia，2015：1?5.

［7］謝秉鑫，王宗，范士林.分布式電源對(duì)于配電網(wǎng)諧波特性的影響［J］.電力系統(tǒng)保護(hù)與控制，2012，40（2）：116?119.

［8］許得志，汪飛，毛華龍，等.多并網(wǎng)逆變器與電網(wǎng)的諧波交互建模與分析［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2013，33（12）：64?71.

［9］謝寧，羅安，陳燕生，等.大型光伏電站動(dòng)態(tài)建模及諧波特性分析［J］.中國(guó)電機(jī)工程學(xué)報(bào)，2013，33（36）：10?17.

［10］胡偉，孫建軍，馬謙，等.多個(gè)并網(wǎng)逆變器間的交互影響分析［J］.電網(wǎng)技術(shù)，2014，38（9）：2511?2518.

Harmonic distortion analysis of grid?connected photovoltaic source based on current injection method

SHI Lei?lei1，2，HAN Tian?hua1，TANG Chao1，WEI Dan1，ZHAO Fan?fan2，JIA Qing?quan2
（1.Xingtai Power Supply Company，State Grid Hebei Electric Power Company，Xingtai，Hebei 054001，China；2.School of Electrical Engineering，Yanshan University，Qinhuangdao，Hebei 066004，China）

The voltage distortion influence mechanism of photovoltaic grid?connected system is studied in this paper.The harmonic voltage distortion analysis method is presented based on the current injection method.The harmonic source model of the photovoltaic grid?connected system which is equivalent to the harmonic current injection source is established.The relation formula of the photo?voltaic harmonic injection sources and voltage distortion rate of each node of the distribution network is obtained.The relationship between the interconnected location of photovoltaic source and the active power output to the harmonic voltage of the power network is also obtained.The harmonic distribution characteristics and interaction relationship are studied when multiple photovoltaic sources are connected to distribution network，and the harmonic distribution law is also obtained.The correctness of the method presented in this paper is verified by comparing theoretical analysis with simulation results.

photovoltaic generation；power quality；harmonic distortion；current injection method

TM46

A

10．3969／j．issn．1007?791X．2015．05．010

1007?791X（2015）05?0443?05

2015?08?10 基金項(xiàng)目：國(guó)家自然科學(xué)基金資助項(xiàng)目（51477147）

石磊磊（1985?），男，河北邯鄲人，博士研究生，工程師，主要研究方向?yàn)榉植际焦夥⒕W(wǎng)及其電能質(zhì)量問題；?通信作者：賈清泉（1970?），男，吉林白山人，博士，教授，博士生導(dǎo)師，主要研究方向?yàn)殡娔苜|(zhì)量、新能源電力系統(tǒng)，Email：jiaqingquan＠sina.com。