趙君,陳繼開
(1.國網(wǎng)黑龍江省電力有限公司,哈爾濱150001;2.東北電力大學,吉林吉林132012)
風電場集中區(qū)域一般都處于電網(wǎng)末端,電網(wǎng)架結(jié)構(gòu)相對薄弱[1],而風能具有隨機性和間歇性的特點[2],所以風電機組有功出力將隨風速隨機變化,繼而引起風電場區(qū)域電網(wǎng)無功分布的改變,如果對該區(qū)域電網(wǎng)的無功補償策略選擇不當,不但會影響電網(wǎng)電壓質(zhì)量和運行的經(jīng)濟性,嚴重時甚至會危及風電機組的持續(xù)聯(lián)網(wǎng)運行[3-8]。因此如何提高多風場區(qū)域電網(wǎng)電壓穩(wěn)定性和經(jīng)濟運行效率將是無功補償協(xié)調(diào)策略需要解決的關(guān)鍵問題。文獻[9]針對風電場自備無功補償設(shè)備投運率低、運行可靠性低、動態(tài)響應速度慢、額定容量不符、控制策略不當以及經(jīng)濟性差等問題,提出了提升動態(tài)無功補償裝置技術(shù)性能及運行管理水平的措施,有效的解決了上述問題。文獻[10]針對風電接入給電網(wǎng)帶來的電壓控制和經(jīng)濟調(diào)節(jié)問題,提出了相應的無功補償方案及電壓控制措施,仿真結(jié)果驗證了該方法的有效性。文獻[11]提出了綜合考慮動態(tài)無功響應能力的大規(guī)模風電場全過程無功電壓緊急控制策略,解決了大規(guī)模風電機組連鎖脫網(wǎng)事故。
根據(jù)上述分析可知,目前許多專家學者針對風電場無功補償策略的問題做了大量研究,但從電壓質(zhì)量和經(jīng)濟性兩個方面對各種無功補償策略進行有效性評價并從中選擇最優(yōu)的無功補償策略目前還沒有相關(guān)的報道。鑒于此,首先,本文構(gòu)建了電壓質(zhì)量指標和經(jīng)濟性指標并對其進行準確數(shù)學描述。而后,以三能級系統(tǒng)復雜度統(tǒng)計為例,分析了不同參數(shù)條件下生存指數(shù)熵的統(tǒng)計特性?;谏嬷笖?shù)熵對小概率事件敏感、具有靈活的統(tǒng)計特性等特點,將其引入到風電場群級無功補償策略的評價中,提出了基于生存指數(shù)熵權(quán)的多風電場無功補償策略評價方法。最后,基于實際風場群區(qū)域電網(wǎng)算例,運用本文提出的評價方法對無功補償設(shè)備、以電壓的安全性為目標的無功協(xié)調(diào)控制策略、以電網(wǎng)運行的經(jīng)濟性為目標的無功協(xié)調(diào)控制策略、改進POS的無功協(xié)調(diào)控制策略進行評估。通過理論分析和仿真結(jié)果證明了該方法的有效性。
由于利用各種控制策略對風電場集中區(qū)域進行無功補償協(xié)調(diào)控制目的在于保證系統(tǒng)電壓質(zhì)量合格并降低區(qū)域網(wǎng)損,所以本文提出評價無功補償控制策略優(yōu)劣的綜合效益指標主要包括:電壓質(zhì)量指標和經(jīng)濟效益指標[12]。
1.1.1 端口電壓合格率η
端口電壓合格率是指端口電壓合格的節(jié)點數(shù)占總點數(shù)的比例,提高節(jié)點電壓合格率是無功補償最主要的目的,所以該指標最能夠直觀反映無功補償前后各端口電壓的水平。根據(jù)相關(guān)的技術(shù)規(guī)定端口電壓允許偏差值為額電壓的-10%~+7%,如果某一端口電壓不在該范圍內(nèi),則認為該節(jié)點電壓不合格,計算公式如下:
式中a2表示總節(jié)點數(shù);a1表示電壓合格的節(jié)點數(shù)。
1.1.2 端口電壓最大值 Ui,max
式中 Ui,max(i=1,2,…,a2)表示端口電壓的最大值;Ut1,Ut2,…,Utn表示風電場端口電壓的瞬時值。
1.1.3 端口電壓最大偏移量 ΔUi,max
端口電壓最大偏移量反映端口電壓的節(jié)點與額定電壓之間偏移量大小。該指標越大說明端口電壓偏離額定電壓越大。計算公式如下:
式中 Ua,i(i=1,2,…,a2)表示端口的額定電壓;Ui,max(i=1,2,…,a2)表示端口電壓的最大值;Ui,min(i=1,2,…,a2)表示端口電壓的最小值,端口電壓最大偏移量越小說明端口電壓的最小值越靠近額定電壓。
1.2.1 無功補償裝置成本費G1
無功補償裝置成本包括初始成本以及運行過程中運行費。其中初始成本是指建設(shè)階段的設(shè)備采購成本以及施工安裝成本。則無功補償裝置成本費G1計算如下:
式中l(wèi)1是初始成本,包括兩部分,一部分是購買無功補償裝置的費用,第二部分為施工安裝及調(diào)試費用;l2是無功補償運行成本費用,主要包括電容器等無功補償裝置自身消耗電能造成的經(jīng)濟損失費用。
1.2.2 風電場群級變壓器有功功率損失費G2
當風電場向電網(wǎng)供電時,無功功率占用變壓器的容量,由于變壓器的容量有限,會導致有功功率傳送降低,進而導致經(jīng)濟上損失。定義風電場群級變壓器有功功率損失費為:
式中m為風電場群區(qū)域中變壓器的數(shù)量;Qk為風電場群區(qū)域中第k個變壓器運行時的無功功率;Ca為風電場群級變壓器有功功率損失的單位價格。
1.2.3 風電場群級無功補償網(wǎng)損費G3
在風電場群級,定義風電場群區(qū)域電網(wǎng)的網(wǎng)損和風電場群接入點的離散設(shè)備的調(diào)節(jié)費用為經(jīng)濟性指標:
式中Ploss為風電場群區(qū)域內(nèi)有功功率傳輸?shù)木W(wǎng)損;Qloss為風電場群區(qū)域內(nèi)無功功率傳輸?shù)木W(wǎng)損;Cc、Ct分別為電容器組和有載調(diào)壓變壓器的單位調(diào)節(jié)成本;na、nb、nc、nt、分別為風電場有功功率的單位價格、風電場無功功率的單位價格、電容器組和有載調(diào)壓變壓器的調(diào)節(jié)次數(shù)。
由于分布函數(shù)比密度函數(shù)更一般,并且影響測度的是分布函數(shù),所以這里生存指數(shù)熵用分布函數(shù)替代密度函數(shù),生存指數(shù)熵的連續(xù)形式為[13]:
式中x∈Rm+;X表示隨機向量;|X|(x)表示隨機變量|X|的多元生存函數(shù)。對于,定義:P(|X|>x),a為生存指數(shù)熵的統(tǒng)計參數(shù)。
基于式(7)可得到離散形式的生存指數(shù)熵權(quán),表示為:
式中F(i)為分布函數(shù)。
設(shè)三能級系統(tǒng)由3個獨立的子系統(tǒng)A、B和C組成時,根據(jù)式(8),繪制不同α指數(shù)下生存指數(shù)熵隨子系統(tǒng)存在概率分布的變化關(guān)系圖(如圖1,圖中p1、p2分別為系統(tǒng)中A、B子系統(tǒng)的存在概率)。觀察圖1可知,當0<α<1且α≠0,隨著值的增大,生存指數(shù)熵對系統(tǒng)狀態(tài)中小概率事件的統(tǒng)計范圍不斷擴大,而對小概率事件的統(tǒng)計敏感度相應降低;相反,當時,隨著值的增大,生存指數(shù)熵對系統(tǒng)狀態(tài)中小概率事件的統(tǒng)計范圍不斷縮小,而對小概率事件的統(tǒng)計敏感度相應升高。這意味著通過調(diào)整α可以改變對被測系統(tǒng)的統(tǒng)計范圍并提高辨識度。通過研究發(fā)現(xiàn)分布函數(shù)的選擇對生存指數(shù)熵的影響較小,所以本文在此不做論述。
在多種無功補償?shù)姆椒ㄖ?,最?yōu)的無功補償方法屬于小概率事件,利用生存指數(shù)熵權(quán)進行評估時,基于上述分析和已有的先驗知識,可根據(jù)無功補償策略的數(shù)量n來確定α的取值,如表1所示。
圖1 生存指數(shù)熵的三能級系統(tǒng)能量圖Fig.1 Three-level system energy diagram of survival exponential entropy
表1 無功補償方法數(shù)量n和生存指數(shù)熵權(quán)參數(shù)α的取值Tab.1 Number n ofmethods for reactive power compensation and parameterαof survival exponential entropy
基于上述分析,這里提出一種結(jié)合生存指數(shù)熵權(quán)的無功補償策略的評價方法,評判流程如圖2所示,具體步驟如圖2所示。
圖2 評判流程圖Fig.2 Flow chart of evaluation
步驟一:設(shè)共有n種無功補償策略參與評估,評估指標數(shù)量為m個,dij為第i個無功補償策略指標j的參數(shù),根據(jù)如下原則進行指標參數(shù)歸一化:
對于越大越好的指標j對應參數(shù)按照式(10)進行歸一化:
對于越小越好的指標j對應參數(shù)按照式(11)進行歸一化:
繼而組成規(guī)范指標矩陣χ=[χ1χ2Lχm?,其中相量 Xj=1,…,m={xij|x≤ i≤ n};
步驟二:利用式(12)對規(guī)范指標矩陣χ=χ1,χ2Lχm進行計算。
步驟三:根據(jù)式(13)確定綜合權(quán)重:
式中 i=1,2,…,m,z=[z(1),z(2),…,z(m)]。
步驟四:根據(jù)式(14)確定綜合評估結(jié)果:
式中D表示無功補償策略指標的客觀值矩陣;采用線性加權(quán)法,將各指標權(quán)重z與各指標參數(shù)xij進行線性加權(quán),得到?jīng)Q策值并進行排序。由于本文選取指標參數(shù)值為越小越好,所以最后得到的決策值越小,說明無功補償策略無功補償?shù)男Ч胶谩?/p>
為了驗證本文提出無功補償策略評價方法的有效性,首先根據(jù)東北某風電集中區(qū)域網(wǎng)架參數(shù),搭建一個含5風電場在內(nèi)的區(qū)域電網(wǎng)模型,如圖3所示。
圖3 多風電場區(qū)域網(wǎng)架圖Fig.3 Topology diagram of grid in multiple wind farms area
基于此模型,針對四種無功補償策略進行無功控制仿真,利用本文提出的評價方法對仿真結(jié)果進行無功補償協(xié)調(diào)控制策略的評估。
首先將不采用任何無功補償手段的控制策略、以電壓的安全性為目標的無功協(xié)調(diào)控制策略、以電網(wǎng)運行的經(jīng)濟性為目標的無功協(xié)調(diào)控制策略、改進POS的無功協(xié)調(diào)控制策略和各指標按順序進行編號,分別為1、2、3、4,并基于 RT-LAB仿真工具搭建包含多風電場的區(qū)域電網(wǎng)模型,測得4種控制策略下各指標對應的參數(shù)值,如表2~表5所示。
表2 電壓合格率η參數(shù)值Tab.2 Parameterfor the rate of qualified voltage
表3 風場端口電壓最大值Tab.3 Maximum voltage value of wind farm port
表4 風場端口電壓最大偏移量Tab.4 Maximum voltage offset ofwind farm
表5 經(jīng)濟指標參數(shù)值Tab.5 Parameter value of economic index
令表2~表5所提供的4種策略的指標數(shù)值為di,j(i=1,…,4;j=1,…,14),i為第 i種無功補償策略,j為第i種無功補償策略下的第j個指標。將各指標參數(shù)進行歸一化處理,可得到規(guī)范指標矩陣X=[X1X2… X14],其中Xj=1,…,14={xij|1≤i≤4}。再利用式(12)對規(guī)范指標矩陣X=[X1X2…X14]進行計算,可得到y(tǒng)。然后根據(jù)式(13)確定綜合權(quán)重,如表6所示。
最后根據(jù)式(14)得出運算結(jié)果:IP=[71.05,69.75,69.72,68.52],矩陣 D的數(shù)值為:
表6 各指標綜合權(quán)重Tab.6 Comprehensive weight of each index
通過分析圖4可知(橫坐標為控制策略序號,縱坐標為綜合評估結(jié)果的決策值),決策值從大到小的順序為:不采用任何無功補償手段的控制策略,以電網(wǎng)運行的經(jīng)濟性為目標的無功協(xié)調(diào)控制策略,以電壓的安全性為目標的無功協(xié)調(diào)控制策略,改進POS的無功協(xié)調(diào)控制策略。因此,改進POS的無功協(xié)調(diào)控制策略無功補償?shù)男Ч詈?,以電網(wǎng)運行經(jīng)濟性為目標的無功協(xié)調(diào)控制策略和以電壓安全性為目標的無功協(xié)調(diào)控制策略無功補償效果其次(以經(jīng)濟性為目標的控制策略的效果優(yōu)于以電壓安全性為目標的控制策略的效果),不采用任何無功補償手段的控制策略效果最差,此評估結(jié)果與實際情況吻合。
圖4 多風電場區(qū)域無功補償策略綜合評估結(jié)果Fig.4 Comprehensive evaluation of reactive power compensation strategy formultiple wind farms
本文以評價風電場集中區(qū)域無功補償策略有效性為目標,從電壓質(zhì)量和經(jīng)濟性兩個方面提出了電壓質(zhì)量評價指標和經(jīng)濟評價指標?;谏嬷笖?shù)熵對小概率事件敏感、具有靈活的統(tǒng)計特性的特點,提出一種基于生存指數(shù)熵權(quán)的多風場區(qū)域電網(wǎng)無功補償策略的評價方法,該方法能夠快速準確的完成對各種無功補償控制策略的評價和排序,繼而為風場集中區(qū)域電網(wǎng)安全高效運行提供了技術(shù)支持。