[發(fā)明專利]大規(guī)模風(fēng)電接入電網(wǎng)后的網(wǎng)損預(yù)測(cè)方法

無(wú)效

• 申請(qǐng)?zhí)枺?/td>: 201210184250.6
• 申請(qǐng)日: 2012-06-05
• 公開(kāi)（公告）號(hào): CN102709908A
• 公開(kāi)（公告）日: 2012-10-03
• 發(fā)明（設(shè)計(jì)）人: 劉文穎;秦睿;文晶;姚旭;葛潤(rùn)東;智勇;但揚(yáng)清;梁琛;吳曉丹;邢延?xùn)|;楊楠;崔力心;劉茜;邢晶;趙子蘭;謝昶;徐鵬;李波;金娜;杜珣;李亞龍;盧甜甜;梁才;周海洋;門德月;王久成;李揚(yáng);劉景延;王佳明;楊斌;曹俊龍;梁紀(jì)峰;杜波;吳耀昊;鐘佳辰
• 申請(qǐng)（專利權(quán)）人: 華北電力大學(xué);甘肅省電力公司電力科學(xué)研究院
• 主分類號(hào): H02J3/00
• 分類號(hào): H02J3/00;G06F19/00
• 代理公司: 北京眾合誠(chéng)成知識(shí)產(chǎn)權(quán)代理有限公司 11246
• 代理人: 朱琨

1.一種大規(guī)模風(fēng)電接入電網(wǎng)后的網(wǎng)損預(yù)測(cè)方法，其特征是所述方法包括：

步驟1：輸入系統(tǒng)數(shù)據(jù)和風(fēng)電場(chǎng)數(shù)據(jù)，包括風(fēng)力發(fā)電機(jī)的額定功率、風(fēng)力發(fā)電機(jī)的額定風(fēng)速、風(fēng)速、切入風(fēng)速、切出風(fēng)速、風(fēng)速分布參數(shù)、發(fā)電機(jī)組的有功出力、發(fā)電機(jī)組正常運(yùn)行的概率、發(fā)電機(jī)組正常運(yùn)行時(shí)的有功出力值、負(fù)荷有功功率的期望值、負(fù)荷有功功率的方差和負(fù)荷有功功率值；所述風(fēng)速分布參數(shù)包括形狀參數(shù)、尺度參數(shù)和位置參數(shù)；

步驟2：確定風(fēng)電場(chǎng)出力的概率分布曲線、發(fā)電機(jī)出力的概率分布曲線和負(fù)荷的概率分布曲線；

所述風(fēng)電場(chǎng)出力的概率分布曲線為概率密度函數(shù)為f(PW)=exp[-(PW-k1v0-k2k1c)k]kk1c(PW-k1v0-k2k1c)k-1;]]>其中，P_W為風(fēng)力發(fā)電機(jī)輸出功率且PW=0,v≤vcik1v+k2,vci<v≤vrPr,vr<v≤vco0,v>vco,]]>v為風(fēng)速，P_r為風(fēng)力發(fā)電機(jī)的額定功率，v_r為風(fēng)力發(fā)電機(jī)的額定風(fēng)速，v_ci為切入風(fēng)速，v_co為切出風(fēng)速，k₁=P_r/(v_r-v_ci)，k₂=-k₁v_ci，f(v)為風(fēng)速的概率密度函數(shù)且k為形狀參數(shù)，c為尺度參數(shù)，v₀為位置參數(shù)；

所述發(fā)電機(jī)出力的概率分布曲線為P(PG=xi)=pp,xi=Cp1-ppxi=0;]]>其中，P_G為發(fā)電機(jī)組的有功出力，p_p為發(fā)電機(jī)組正常運(yùn)行的概率，C_p是發(fā)電機(jī)組正常運(yùn)行時(shí)的有功出力值；

所述負(fù)荷的概率分布曲線為其概率密度函數(shù)為其中，μ_P為負(fù)荷有功功率的期望值，δ_P為負(fù)荷有功功率的方差，P_L為負(fù)荷有功功率值；

步驟3：根據(jù)風(fēng)電場(chǎng)出力的概率分布曲線、發(fā)電機(jī)出力的概率分布曲線和負(fù)荷的概率分布曲線，計(jì)算得到各節(jié)點(diǎn)注入功率的概率分布曲線，對(duì)各節(jié)點(diǎn)的概率分布曲線在[0,T]時(shí)段內(nèi)積分并且平均化，得到各節(jié)點(diǎn)在[0,T]時(shí)段注入功率的均值，也即各節(jié)點(diǎn)注入功率的期望值;其中，T為設(shè)定值；

步驟4：根據(jù)已經(jīng)求得的各節(jié)點(diǎn)注入功率的期望值，利用公式求得系數(shù)矩陣J₀，即將網(wǎng)損功率P_LOSS對(duì)各節(jié)點(diǎn)的注入功率P分別求偏導(dǎo)，之后將各節(jié)點(diǎn)注入功率的期望值代入求偏導(dǎo)后的式子得到系數(shù)矩陣J₀；其中，PLOSS=PLOSS0+J0ΔP,PLOSS0]]>為潮流條件下的網(wǎng)損功率，ΔP=[P1-P10,P2-P20,...,Pn-Pn0],]]>P=[P₁,P₂,…,P_n]，P_i(i=1,2,...,n)為各節(jié)點(diǎn)注入功率，為各節(jié)點(diǎn)注入功率的期望值，n為節(jié)點(diǎn)數(shù)量；

步驟5：根據(jù)風(fēng)電場(chǎng)出力的概率分布曲線、發(fā)電機(jī)出力的概率分布曲線和負(fù)荷的概率分布曲線，分別計(jì)算各節(jié)點(diǎn)的風(fēng)電場(chǎng)注入功率、發(fā)電機(jī)注入功率和負(fù)荷功率的各階半不變量；

步驟6：根據(jù)公式計(jì)算各節(jié)點(diǎn)的注入功率ΔP的各階半不變量；其中，ΔP^(k)為各節(jié)點(diǎn)的注入功率ΔP的第k階半不變量，為各節(jié)點(diǎn)的風(fēng)電場(chǎng)注入功率的第k階半不變量，為各節(jié)點(diǎn)的發(fā)電機(jī)注入功率的第k階半不變量，為各節(jié)點(diǎn)的負(fù)荷功率的各階半不變量；k為半不變量的階數(shù)；

步驟7：利用公式計(jì)算各節(jié)點(diǎn)網(wǎng)損功率增量ΔP_LOSS的各階半不變量，并利用公式計(jì)算各節(jié)點(diǎn)的網(wǎng)損功率P_LOSS；其中，為系數(shù)矩陣J₀中元素的k次冪所構(gòu)成的矩陣，為潮流條件下的網(wǎng)損功率；

步驟8：利用Gram-Charlier級(jí)數(shù)展開(kāi)求得網(wǎng)損功率P_LOSS的概率密度函數(shù)和累積分布函數(shù)；

步驟9：利用公式求出[0,T]時(shí)段內(nèi)的網(wǎng)損電量W_LOSS；P_LOSS(p)為網(wǎng)損功率P_LOSS的累積分布函數(shù)。