考慮配電網(wǎng)電能質(zhì)量改善的分布式光伏優(yōu)化調(diào)度方法

張敏，祗會強，張世鋒，樊瑞，李冉，王金浩

(國網(wǎng)山西省電力公司電力科學研究院，太原市 030001)

0 引 言

為構(gòu)建清潔低碳的能源體系，應對能源資源緊張、氣候變化等挑戰(zhàn)，可再生能源發(fā)電受到了各國政府的重視并得到迅速發(fā)展。分布式光伏發(fā)電作為利用可再生能源的一種重要形式，已在中低壓配電網(wǎng)中得到了推廣應用。然而，隨著分布式光伏容量的增加以及大量單相分布式光伏的接入，配電網(wǎng)也將面臨電壓越限[1]、三相不平衡[2]等電能質(zhì)量問題。另一方面，不斷增加且分布廣泛的非線性負荷以及單相負荷也導致配電網(wǎng)的諧波及不平衡問題日趨嚴重[3]，影響其正常運行。為解決配電網(wǎng)電壓越限、三相不平衡和諧波等問題，通常需要安裝各種電能質(zhì)量補償裝置，但這樣會增加配電網(wǎng)的投資運行成本。由于分布式光伏的并網(wǎng)逆變器與有源濾波器、靜止同步補償器等裝置具有一致的拓撲結(jié)構(gòu)，其在適當?shù)目刂葡乱部砂l(fā)揮電能質(zhì)量治理功能[4-5]，因此利用分布式光伏自身改善配電網(wǎng)電能質(zhì)量的方案因其經(jīng)濟性已經(jīng)在近年來受到了較多關(guān)注。

為充分發(fā)揮分布式光伏的電能質(zhì)量補償能力，國內(nèi)外學者們針對其優(yōu)化調(diào)度問題開展了廣泛研究。文獻[6]針對高比例分布式光伏接入導致的電壓越限問題，提出一種光伏逆變器的無功調(diào)度方法。文獻[7]通過對配電網(wǎng)中的可調(diào)串聯(lián)電抗器的電抗值及光伏逆變器的輸出無功進行協(xié)調(diào)優(yōu)化，在解決電壓越限問題的同時還減小了線路有功損耗。文獻[8-10]提出了用于配電網(wǎng)電壓控制的分布式光伏有功無功優(yōu)化調(diào)度方法。然而，上述研究均在配電網(wǎng)三相平衡的假設(shè)下進行，并采用單相配電網(wǎng)模型進行分析，因此建立的優(yōu)化模型不適用于不平衡配電網(wǎng)。文獻[11]針對單相分布式光伏接入配電網(wǎng)引起的不平衡問題，提出了一種光伏選相接入及優(yōu)化其有功無功輸出的方法。而文獻[12-13]以改善網(wǎng)損、不平衡等為目標，對中低壓配電網(wǎng)中分布式光伏的有功無功輸出進行了協(xié)調(diào)優(yōu)化。但上述文獻并未考慮利用分布式光伏進行配電網(wǎng)諧波問題的優(yōu)化治理。由于目前配電網(wǎng)中的非線性負荷逐漸呈現(xiàn)出數(shù)量多、分布廣的特點，傳統(tǒng)的單點諧波治理方法不再適用，因此從優(yōu)化的角度利用分布式光伏改善配電網(wǎng)諧波問題同樣值得研究。

分布式光伏優(yōu)化調(diào)度本質(zhì)上是建立并求解配電網(wǎng)最優(yōu)潮流(optimal power flow, OPF)模型。在同時考慮改善不平衡以及諧波問題的情況下，OPF模型將包括配電網(wǎng)三相基波潮流和諧波潮流約束，以及對應的三相電壓電流變量，是一個復雜的非凸非線性約束優(yōu)化模型，因此有必要研究此OPF模型的求解算法。由于OPF模型的非凸性質(zhì)主要源于潮流方程，文獻[14]利用一階泰勒近似方法將不平衡配電網(wǎng)的三相潮流方程線性化。文獻[15-16]采用支路潮流等式作為配電網(wǎng)的運行約束條件，并采用變量替換的方式將其松弛為線性不等式約束。文獻[17]在三相潮流約束中加入輔助變量，將非凸非線性O(shè)PF模型轉(zhuǎn)換為半定規(guī)劃模型求解。然而，上述近似方法的精度受配電網(wǎng)運行狀態(tài)的影響較大，根據(jù)近似模型優(yōu)化結(jié)果可能不能滿足配電網(wǎng)某些運行約束[18]。此外，遺傳算法、粒子群算法等智能算法也在配電網(wǎng)OPF模型的求解中得到了一些應用[19-20]，但當變量及約束條件數(shù)量較多時，智能算法往往存在計算效率低、解的最優(yōu)性差等問題，導致其在配電網(wǎng)優(yōu)化調(diào)度中的應用受到限制。

綜上所述，本文從利用分布式光伏逆變器輸出有功功率后的剩余容量進行配電網(wǎng)電能質(zhì)量改善的角度出發(fā)，研究以綜合改善配電網(wǎng)電壓越限、三相不平衡、諧波等電能質(zhì)量問題的分布式光伏優(yōu)化調(diào)度方法，主要貢獻如下：

1)以可調(diào)度分布式光伏(schedulable distributed photovoltaic, SDPV)三相有功、無功以及諧波電流輸出為控制變量，并在考慮配電網(wǎng)三相基波潮流及諧波潮流約束、節(jié)點電壓約束、光伏逆變器運行約束以及滿足光伏三相總有功功率等于直流側(cè)可輸出最大有功功率的基礎(chǔ)上，建立最小化網(wǎng)絡(luò)損耗、諧波電壓、不平衡電壓的配電網(wǎng)OPF模型。該模型采用統(tǒng)一的形式描述相線、中性線以及連接在不同相間的負荷，可適用于不同結(jié)構(gòu)配電網(wǎng)。

2)為穩(wěn)定高效地求解所建立的配電網(wǎng)OPF模型，提出一種潮流計算與優(yōu)化計算的兩階段迭代優(yōu)化求解算法，通過潮流計算對OPF模型中電流變量進行削減，并提出節(jié)點電壓幅值約束的幾何近似方法，將原始非凸非線性O(shè)PF模型轉(zhuǎn)換為凸二次約束二次規(guī)劃模型，使其可以被原對偶內(nèi)點法快速求解。最后通過潮流計算與優(yōu)化計算的交替迭代減小模型凸化誤差，從而在降低OPF模型求解難度的同時保證解的精度。

3)建立一個162節(jié)點三相配電網(wǎng)測試模型并進行算例分析，對采用所提分布式光伏優(yōu)化調(diào)度方法前后的配電網(wǎng)損耗、節(jié)點電壓越限情況、網(wǎng)絡(luò)最大負序不平衡度以及最大諧波畸變率進行對比，然后對兩階段優(yōu)化求解算法的求解誤差以及計算效率進行分析，驗證所提方法的有效性。

1 多功能分布式光伏發(fā)電系統(tǒng)

首先從簡單的單相光伏發(fā)電系統(tǒng)開始分析，其常見結(jié)構(gòu)如圖1所示，主要包括光伏板、直流變流器、逆變器以及輸出濾波器。其中，直流側(cè)電容用于實現(xiàn)直流變流器與逆變器的解耦運行。而直流變流器一般采用最大功率點跟蹤(maximum power point tracking，MPPT)控制，保證光伏輸出最大有功功率。逆變器則采用電壓外環(huán)電流內(nèi)環(huán)的雙環(huán)控制策略，以維持直流側(cè)電壓的穩(wěn)定以及輸入輸出功率的平衡。由于光伏逆變器通常不會持續(xù)工作在額定狀態(tài)，因此在其具有剩余容量時，即可在適當?shù)目刂葡螺敵鲆欢康臒o功功率及諧波電流，用于改善配電網(wǎng)電能質(zhì)量?？蓪崿F(xiàn)上述功能的控制方案示意如圖2、3所示。

圖1 單相分布式光伏發(fā)電系統(tǒng)結(jié)構(gòu)

圖2 直流電壓外環(huán)示意圖

在圖2中，Qref、Ihref及udcref分別為光伏逆變器無功功率輸出、諧波電流輸出以及直流側(cè)電壓參考值；uac、U1及θU分別為并網(wǎng)點電壓及其基波分量的幅值和相位；I1x、I1y及Ihx、Ihy分別為光伏輸出基波和h次諧波電流相量參考信號的實部與虛部，通過乘以正弦信號，可最終生成光伏輸出電流的正弦參考信號iref，然后由圖3所示的電流內(nèi)環(huán)對iref進行無差跟蹤。圖2中，Qref、Ihref可提前設(shè)定，也可由調(diào)度中心生成并下發(fā)給配電網(wǎng)中的可調(diào)度分布式光伏。另外，調(diào)度中心通過收集配電網(wǎng)運行數(shù)據(jù)，建立并求解配電網(wǎng)最優(yōu)潮流模型，還可計算得到各SDPV輸出參考值的最優(yōu)組合，實現(xiàn)對SDPV的優(yōu)化調(diào)度。

圖3 電流內(nèi)環(huán)示意圖

將圖1中的逆變器替換為三相逆變器，即可得到三相分布式光伏發(fā)電系統(tǒng)。三相逆變器具有三橋臂、四橋臂以及3個單相逆變器并聯(lián)而成的組式逆變器等不同結(jié)構(gòu)，其中四橋臂及組式逆變器具有獨立調(diào)節(jié)各相輸出的能力，但前者控制相對復雜。為改善配電網(wǎng)三相不平衡問題，本文考慮SDPV的并網(wǎng)接口為組式逆變器，如圖4所示，并且其中2臺單相逆變器的外環(huán)采用圖5所示的控制方案，其中，Pref為逆變器有功功率輸出參考。由于光伏三相逆變器其中的兩相根據(jù)Pref的值輸出固定有功功率，因此剩余一相逆變器為了維持直流側(cè)電壓穩(wěn)定，其輸出有功功率會自動調(diào)節(jié)從而滿足三相逆變器輸出有功功率之和等于直流側(cè)輸出功率，這使得光伏各相輸出功率可在總有功功率輸出約束下進行獨立調(diào)節(jié)。在上述條件下，本文后續(xù)將對SDPV的有功功率、無功功率及諧波電流的輸出進行優(yōu)化，從而實現(xiàn)配電網(wǎng)的運行優(yōu)化及諧波畸變、不平衡等電能質(zhì)量問題的綜合優(yōu)化治理。

圖4 三相分布式光伏發(fā)電系統(tǒng)結(jié)構(gòu)

圖5 功率外環(huán)示意圖

2 改善電能質(zhì)量的最優(yōu)潮流模型

由于配電網(wǎng)具有單相、三相三線、三相四線等不同網(wǎng)絡(luò)結(jié)構(gòu)，因此為建立具有通用性的最優(yōu)潮流模型，本文以三相四線制配電網(wǎng)為例進行建模，對建立的模型進行簡化即可得到其他結(jié)構(gòu)配電網(wǎng)的OPF模型。

2.1 目標函數(shù)

1)從配電網(wǎng)運行優(yōu)化的角度出發(fā)，要求線路有功功率損耗floss最小化，該項目標函數(shù)表達式為：

(1)

2)當配電網(wǎng)中存在非線性負荷時，從改善諧波問題的角度出發(fā)，應最小化網(wǎng)絡(luò)整體的諧波電壓水平。

(2)

3)當配電網(wǎng)中三相功率不平衡或網(wǎng)絡(luò)參數(shù)本身不對稱時，從改善不平衡問題的角度出發(fā)，應抑制網(wǎng)絡(luò)各節(jié)點負序電壓大小。

(3)

由于多目標優(yōu)化問題的直接求解較為困難，因此本文通過線性加權(quán)[21-23]的方法將上述多個目標函數(shù)進行組合得到單目標函數(shù)f：

(4)

式中：Floss、Fh、Fu為不同目標函數(shù)的歸一化系數(shù)；ωloss、ωh、ωu為對應目標函數(shù)的權(quán)重系數(shù)，其值越大表示對應的目標函數(shù)項優(yōu)先級越高。在實際應用中，可以根據(jù)對電能質(zhì)量問題的關(guān)注度或者嚴重程度主觀設(shè)置各權(quán)重系數(shù)，也可以通過層次分析法[23]等權(quán)重決策方法得到更為客觀的權(quán)重系數(shù)。

2.2 約束條件

配電網(wǎng)OPF模型的約束條件包括等式約束及不等式約束，其中，等式約束主要為網(wǎng)絡(luò)潮流方程，不等式約束則包括節(jié)點電壓幅值約束、SDPV運行約束等。

1)基波潮流約束如式(5)—(7)所示，其中，式(5)、(6)為非線性等式約束。

(5)

(6)

(7)

對于三相SDPV，其三相逆變器輸出有功功率應等于直流側(cè)輸出有功功率：

(8)

2)諧波潮流約束如式(9)、(10)所示。本文主要考慮低壓臺區(qū)負荷產(chǎn)生的諧波電流，并采用諧波恒流源模型作為低壓臺區(qū)負荷的諧波模型。

(9)

(10)

3)節(jié)點電壓約束。為保證配電網(wǎng)的安全運行，通常要求節(jié)點電壓幅值不能超出一定范圍，即滿足：

(11)

4)線路電流約束。一般地，需要保證線路上流過的總電流不能過大，因此需要滿足：

(12)

5)電能質(zhì)量約束。為保證良好的電能質(zhì)量，要求節(jié)點諧波電壓及不平衡電壓不能過大，即滿足：

(13)

(14)

4)SDPV并網(wǎng)逆變器運行約束。為了不影響SDPV的正常工作，應保證其逆變器不過載，因此需要對其輸出功率及電流進行限制。

(15)

(16)

結(jié)合上述目標函數(shù)及約束條件可得，以降低配電網(wǎng)網(wǎng)損，治理諧波及不平衡等電能質(zhì)量問題為目標，對SDPV進行優(yōu)化調(diào)度的配電網(wǎng)OPF模型為：

(17)

需要注意的是，在OPF模型式(17)中，不等式約束式(11)—(14)和式(15)—(16)之間存在沖突，這是因為在光伏逆變器輸出有功功率后剩余容量較少時，其無功功率和諧波電流補償量也會相應減少，在這種情況下可能就無法保證配電網(wǎng)節(jié)點電壓或者線路電流等不超標，即出現(xiàn)模型無解的情況，此時可適當放寬或者取消部分電能質(zhì)量約束，從而在不影響光伏逆變器有功輸出的前提下盡力改善配電網(wǎng)電能質(zhì)量。

3 模型求解流程

3.1 模型凸化過程

求解式(17)即可得到配電網(wǎng)中各SDPV的優(yōu)化調(diào)度方案。然而，式(17)是一個典型的非凸非線性約束優(yōu)化模型，直接對其進行求解存在計算效率低、收斂難等問題，因此本節(jié)中提出一種凸化近似方法，將式(17)轉(zhuǎn)化為更易于求解的凸優(yōu)化模型。

為了消除非線性等式約束式(5)、(6)、(9)，首先根據(jù)當前配電網(wǎng)中低壓臺區(qū)以及SDPV的運行數(shù)據(jù)進行一次基波潮流與諧波潮流計算，得到各低壓臺區(qū)的基波電流和諧波電流初值，然后將電流初值代入式(7)和(10)中，得到直角坐標系下的線性化潮流方程，如式(18)所示。

(18)

(19)

對于SDPV逆變器容量約束式(15)及(16)，也可通過類似的方式進行凸化，得到新的不等式約束：

(20)

(21)

經(jīng)過上述近似簡化后，只剩節(jié)點電壓約束式(11)為非凸約束。直角坐標系下的節(jié)點電壓幅值約束如圖6所示。

圖6 直角坐標下的電壓幅值約束

由圖6可見，節(jié)點電壓幅值約束表現(xiàn)為電壓上限與下限之間的圓環(huán)區(qū)域。由于實際配電網(wǎng)各節(jié)點電壓相角變化范圍通常不大，因此可采用某一小段區(qū)域近似表示電壓約束，例如圖6中的藍色矩形區(qū)域。圖中的每個矩形區(qū)域均可由一組線性不等式描述，以a相為例：

(22)

式中：δ為節(jié)點電壓相角的變化范圍，為保證電壓約束近似的合理性，其值應大于節(jié)點電壓相量相角與參考相角間的最大相角差。但δ的值選得過大，近似精度會下降，在實際應用中，可以通過在優(yōu)化前對配電網(wǎng)進行潮流計算，得到網(wǎng)絡(luò)中節(jié)點電壓相角的變化范圍，并據(jù)此合理設(shè)置δ的值。在本文算例分析中，設(shè)置δ=15°。另外，選擇合適的參考相角也有助于提高近似精度。

經(jīng)過上述凸化過程，可得到OPF模型式(17)的凸化近似模型為：

(23)

由于式(23)的目標函數(shù)為凸二次函數(shù)，等式約束均為線性約束，而不等式約束由線性約束和凸二次約束組成，因此其為凸二次約束二次規(guī)劃(quadratic constrained quadratic program，QCQP)模型，凸QCQP模型的最優(yōu)解可通過非線性優(yōu)化算法，如原對偶內(nèi)點法[24]，可靠求得。

3.2 兩階段迭代求解

對3.1節(jié)中的OPF模型凸化過程進行分析可知，在進行一次潮流計算后，模型中低壓臺區(qū)的電流將不再變化，相當于完全將其視為恒流源，這顯然忽略了節(jié)點電壓對臺區(qū)輸出電流的影響，因此可能導致凸QCQP模型式(23)的優(yōu)化結(jié)果與原OPF模型式(17)之間存在較大差異。為了減小凸化帶來的誤差，本節(jié)提出一種潮流計算與優(yōu)化計算之間的兩階段迭代優(yōu)化求解方法，即令潮流計算與優(yōu)化計算的輸出結(jié)果互為對方的輸入，不斷循環(huán)迭代，直至計算結(jié)果不再發(fā)生明顯變化。具體步驟如下：

步驟1：根據(jù)配電網(wǎng)中低壓臺區(qū)及SDPV當前的運行數(shù)據(jù)，進行基波潮流及諧波潮流計算，得到低壓臺區(qū)輸出基波電流及諧波電流初值，以及SDPV并網(wǎng)點基波電壓及諧波電壓初值。其中，潮流計算依靠OpenDSS配電網(wǎng)潮流計算軟件進行。

步驟2：根據(jù)3.1節(jié)中的過程，將非凸非線性O(shè)PF模型轉(zhuǎn)化為凸QCQP模型。

步驟3：使用原對偶內(nèi)點算法求解凸QCQP模型。

步驟4：計算步驟3前后模型中各變量值的相對誤差。如果誤差小于某一閾值，例如0.1%，則將當前計算結(jié)果作為最終優(yōu)化結(jié)果輸出。否則，將當前優(yōu)化結(jié)果中SDPV的基波功率及諧波電流輸出作為輸入，返回步驟1進行新一輪計算，直至結(jié)果收斂。

4 算例分析

4.1 算例信息

為了驗證所提考慮電能質(zhì)量改善的分布式光伏優(yōu)化調(diào)度算法的有效性，以162節(jié)點三相四線制配電網(wǎng)為測試算例進行分析。算例中配電網(wǎng)根據(jù)北美某電壓等級為13.8 kV的實際配電網(wǎng)改編得到，其拓撲如圖7所示，詳細的網(wǎng)絡(luò)參數(shù)信息可見文獻[25]。在本算例中，設(shè)置各目標函數(shù)的權(quán)重ωloss、ωh、ωu分別為0.70、0.15、0.15。

圖7 測試配電網(wǎng)結(jié)構(gòu)

算例中低壓臺區(qū)負荷總功率及接入低壓臺區(qū)分布式光伏的總發(fā)電功率如圖8所示，其中負荷用電高峰主要出現(xiàn)在傍晚，而光伏發(fā)電高峰出現(xiàn)在中午。算例的其他基本信息如表1所示。另外，由于中性線的存在，該網(wǎng)絡(luò)中具有零序諧波電流通路，因此針對諧波問題，本算例將考慮低壓臺區(qū)負荷注入的3、5、7、9次諧波電流的影響，其諧波電流頻譜按照低壓臺區(qū)典型諧波頻譜進行設(shè)置[26]。

圖8 低壓側(cè)總負荷與總光伏發(fā)電功率

表1 算例基本信息

4.2 結(jié)果分析

本算例中，每15 min進行一次優(yōu)化計算，得到各SDPV的功率及諧波電流優(yōu)化輸出方案，并且認為優(yōu)化前SDPV僅輸出有功功率，且三相輸出功率相同。采用所提優(yōu)化調(diào)度方法，得到的12:00以及20:00 SDPV 1—5的各相有功、無功及諧波電流輸出分別如圖9及圖10所示。

圖9 可調(diào)度光伏在12:00的輸出

圖10 可調(diào)度光伏在20:00的輸出

在12:00，低壓臺區(qū)接入的分布式光伏發(fā)電功率三相不平衡問題較為突出，其中b相的光伏發(fā)電功率明顯小于其他兩相，這將導致配電網(wǎng)不平衡問題加劇。理想情況下，即SDPV逆變器容量充足時，可以依靠在b相輸出更多有功功率的方法來改善三相不平衡問題，但是由于此時SDPV自身有功功率輸出的需求以及各相逆變器容量的限制，上述方案并不可行。而為了在這種情況下改善網(wǎng)絡(luò)不平衡問題，在實際優(yōu)化結(jié)果中，可見各SDPV減小了c相有功輸出并將該相逆變器容量用于輸出更多無功，使得網(wǎng)絡(luò)各節(jié)點c相電壓相角偏移減小，三相電壓之間的相角差更接近，從而達到減小負序電壓的目的。以節(jié)點107為例，優(yōu)化前后電壓變化情況如表2所示。另外，可以發(fā)現(xiàn)，由于與SDPV5距離更近的低壓臺區(qū)節(jié)點數(shù)量更多，因此其輸出的諧波電流相比其他SDPV也更多一些，從而抵消流向上游線路的負荷諧波電流，改善諧波電流流經(jīng)線路阻抗造成的電壓諧波畸變問題。

表2 優(yōu)化前后節(jié)點107電壓變化

在20：00，各SDPV三相總有功為0，此時逆變器容量較為充足，因此可以在不同相之間進行有功功率的轉(zhuǎn)移來改善負荷三相功率不平衡的問題，同時也能輸出充足的無功功率及諧波電流來改善網(wǎng)絡(luò)中的電壓越限以及諧波畸變等問題。

圖11展示了優(yōu)化前后網(wǎng)絡(luò)節(jié)點電壓幅值最小值的日變化趨勢，圖中虛線為設(shè)定的電壓幅值約束0.95 pu以及1.05 pu。可見，優(yōu)化前網(wǎng)絡(luò)最小電壓已低于0.92 pu。而按照所提方法對SDPV的功率輸出進行優(yōu)化后，在SDPV逆變器輸出無功的作用下，網(wǎng)絡(luò)中最小電壓始終高于0.95 pu，保證了各節(jié)點電壓均滿足運行約束條件。另外，在SDPV輸出無功作用下，可見優(yōu)化后節(jié)點電壓幅值最大值也有所升高，但由于設(shè)置了電壓上限約束，因此這種電壓幅值的升高將被控制在合理的范圍之內(nèi)，不會造成新的電壓過高問題。

圖11 節(jié)點電壓幅值日變化情況

網(wǎng)絡(luò)各時刻有功功率損耗情況如圖12所示?？梢妰?yōu)化后各時刻有功功率損耗有不同程度的降低，其中，負荷功率較大且SDPV逆變器容量充足時，例如20：00，由于可以較多地補償負荷無功，因此有功損耗降幅相對明顯，而在SDPV逆變器容量不足時，例如12:00，其對網(wǎng)絡(luò)損耗的改善作用也會受到一定限制。另外，將各時刻有功損耗乘以計算時間間隔，可得到優(yōu)化前配電網(wǎng)總能量損耗約為3.53 MW·h，而優(yōu)化后，損耗降至2.94 MW·h。

圖12 配電網(wǎng)有功損耗日變化情況

圖13為網(wǎng)絡(luò)中負序電壓不平衡度最大值的日變化情況?？梢钥吹剑鶕?jù)所提優(yōu)化調(diào)度方法的計算結(jié)果調(diào)整可調(diào)度分布式光伏各相有功及無功功率輸出后，網(wǎng)絡(luò)中的負序電壓不平衡度最大值從約1.5%降低至0.6%左右。

圖13 負序電壓不平衡度最大值日變化情況

圖14展示了網(wǎng)絡(luò)電壓總諧波畸變率最大值的變化情況。經(jīng)過優(yōu)化后，電壓總諧波畸變率最大值從3.57%降至0.48%。另外可以觀察到當SDPV并網(wǎng)逆變器的剩余容量充足時，例如16:00—20:00時，配電網(wǎng)節(jié)點電壓總諧波畸變率在SDPV進行諧波電流補償后顯著降低。然而，當SDPV逆變器容量不足時，例如12:00，此時SDPV需要優(yōu)先滿足其自身有功功率輸出需求，因此對諧波畸變問題的改善效果有所減弱。

圖14 電壓總諧波畸變率最大值日變化情況

為了觀察目標函數(shù)權(quán)重系數(shù)對優(yōu)化結(jié)果的影響，考慮權(quán)重系數(shù)ωloss、ωh、ωu分別為(0.70, 0.15, 0.15)、(0.15, 0.70, 0.15)、(0.15, 0.15, 0.70)的3種場景，分別為場景1、2、3，得到20:00的優(yōu)化結(jié)果如表3所示。從表3中可見，權(quán)重系數(shù)越大，其對應的單項目標函數(shù)值越小，表明了通過增大權(quán)重系數(shù)可以有效增加目標的優(yōu)先級。另外，通過比較不同場景中目標函數(shù)值的差異，可以發(fā)現(xiàn)目標函數(shù)floss對其權(quán)重變化敏感程度較低，這意味著在增大其余目標函數(shù)項的權(quán)重有助于獲得更好的綜合優(yōu)化結(jié)果。

表3 不同權(quán)重系數(shù)下優(yōu)化結(jié)果

4.3 算法性能分析

首先，對電壓幅值約束式(22)的近似合理性進行分析。通過潮流計算與優(yōu)化計算，得到優(yōu)化前后網(wǎng)絡(luò)各節(jié)點電壓相角的變化范圍如圖15所示，其中，實線和虛線各自圍成的區(qū)域分別是優(yōu)化前和優(yōu)化后節(jié)點電壓相角的變化范圍?？梢?，在本算例中，選取δ=15°可以完全覆蓋各節(jié)點電壓相角的變化范圍，意味著此時δ值的選取是合理的。另外，從圖中可以發(fā)現(xiàn)，僅針對某一時刻而言，配電網(wǎng)中各節(jié)點電壓相角的變化范圍會更小，因此也可以在不同時刻選取不同的參考相角及δ值，從而進一步提高對原始電壓約束式(11)的近似精度。

圖15 各相節(jié)點電壓相角變化范圍

另外，由于在優(yōu)化求解過程中對原始OPF模型還采用了其他的凸化近似手段，因此也有必要計算求解過程中凸QCQP模型解與原始OPF模型解之間的誤差，誤差計算過程如圖16所示。其中，將凸QCQP模型解中的控制變量值作為輸入代入到基波及諧波潮流計算中，以得到滿足原始OPF模型潮流約束條件的一組解。定義節(jié)點電壓相對誤差eU以及約束違反度eh如下：

圖16 所提優(yōu)化算法誤差分析流程

(24)

不同迭代次數(shù)下的eU、eh以及目標函數(shù)值如表4所示。可見在迭代次數(shù)較少時，凸QCQP模型解與原始OPF模型解之間將存在較大誤差，而隨著潮流計算與優(yōu)化計算交替迭代的次數(shù)增加，模型誤差逐漸減小。在第8次迭代結(jié)束后，凸QCQP模型的解已經(jīng)基本滿足原始OPF模型的所有約束，且目標函數(shù)值也不再變化。這說明凸QCQP模型的解在原OPF模型中的最優(yōu)性以及可行性能夠隨交替迭代次數(shù)的增加而得到保證。

表4 不同迭代次數(shù)下的模型誤差

表4所示迭代過程的計算總用時情況如圖17所示，圖中也展示了每次迭代過程中使用原對偶內(nèi)點算法求解凸QCQP模型時的誤差下降過程?？梢钥闯?，每當原對偶內(nèi)點算法的內(nèi)部誤差下降到10-6時，即認為本次迭代中凸QCQP模型的最優(yōu)解已被求得，隨后新一輪迭代開始。經(jīng)過8輪潮流計算與優(yōu)化計算迭代后，算法最終停止，總用時約為150 s，能夠滿足通常為15 min進行一次配電網(wǎng)調(diào)度的時間要求。

圖17 所提兩階段優(yōu)化算法計算總用時

5 結(jié) 論

本文針對如何充分發(fā)揮可調(diào)度分布式光伏在改善配電網(wǎng)網(wǎng)損、電壓越限、三相不平衡以及諧波畸變等問題的能力展開了研究。首先，分析了利用可調(diào)度分布式光伏進行配電網(wǎng)電能質(zhì)量改善的可行性。然后，在配電網(wǎng)三相潮流約束以及其余運行約束條件基礎(chǔ)上，以可調(diào)度分布式光伏三相有功功率、無功功率以及諧波電流輸出為控制變量，建立了一個考慮電能質(zhì)量改善的配電網(wǎng)最優(yōu)潮流模型，并提出了一種潮流計算與優(yōu)化計算交替迭代的兩階段迭代優(yōu)化算法用于求解最優(yōu)潮流模型解。算例分析結(jié)果表明，本文所提分布式光伏優(yōu)化調(diào)度方法可有效改善配電網(wǎng)電壓越限、三相不平衡以及諧波畸變等電能質(zhì)量問題。下一步將進一步研究多種類型分布式電源以及電能質(zhì)量治理裝置共同參與的配電網(wǎng)優(yōu)化運行策略。