1.一種雙向交直流變換器模型預(yù)測電流控制方法,其特征在于:步驟如下,
步驟S1,定義雙向交直流變換器的開關(guān)狀態(tài)Si;
其中,i為交流電網(wǎng)的相,且i∈(a,b,c);
S2,獲取αβ兩相靜止坐標(biāo)下雙向交直流變換器的輸出電壓矢量Uj與開關(guān)狀態(tài)Si的表達(dá)式,具體如下:
其中,uα為輸出電壓的α分量;uβ為輸出電壓的β分量;Udc為直流母線電壓,Sa為a相的開關(guān)狀態(tài)值;Sb為b相的開關(guān)狀態(tài)值;Sc為c相的開關(guān)狀態(tài)值;
S3,構(gòu)造雙向交直流變換器與輸出電壓矢量Uj有關(guān)的二次功率預(yù)測模型和輸出電流二次預(yù)測模型;
輸出電流二次預(yù)測模型為,
式中,iα(k+1)為tk+1時(shí)刻輸出電流預(yù)測值的α分量;iβ(k+1)為tk+1時(shí)刻輸出電流預(yù)測值的β分量;iα(k+2)為tk+2時(shí)刻輸出電流的α分量;iβ(k+2)為tk+2時(shí)刻輸出電流的β分量;eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;uα(k+1)為tk+1時(shí)刻輸出電壓的α分量;uβ(k+1)為tk+1時(shí)刻輸出電壓的β分量;L為電感;R為電阻;Ts為采樣頻率;
二次功率預(yù)測模型,
其中,P(k+1)為tk+1時(shí)刻有功功率預(yù)測值;Q(k+1)為tk+1時(shí)刻無功功率預(yù)測值;P(k+2)為tk+2時(shí)刻有功功率預(yù)測值;Q(k+2)為tk+2時(shí)刻無功功率預(yù)測值;iα(k+1)為tk+1時(shí)刻輸出電流預(yù)測值的α分量;iβ(k+1)為tk+1時(shí)刻輸出電流預(yù)測值的β分量;eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;uα(k+1)為tk+1時(shí)刻輸出電壓的α分量;uβ(k+1)為tk+1時(shí)刻輸出電壓的β分量;L為電感;Ts為采樣頻率;
S4,計(jì)算αβ靜止坐標(biāo)下的輸出電流參考值iαref和iβref;
與有功功率有關(guān)的輸出電流參考值iαref和iβref的計(jì)算公式如下;
式中,iαref為輸出電流參考值的α分量;iβref為輸出電流參考值的β分量;ex為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;eα′為電網(wǎng)電壓α分量的90°延遲信號;eβ′為電網(wǎng)電壓β分量的90°延遲信號;pref為有功功率給定值;
與無功功率有關(guān)的輸出電流參考值iαref和iβref的計(jì)算公式如下;
式中,iαref為輸出電流參考值的α分量;iβref為輸出電流參考值的β分量;eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;eα′為電網(wǎng)電壓α分量的90°延遲信號;eβ′為電網(wǎng)電壓β分量的90°延遲信號;qref為有功功率給定值;
S5,構(gòu)造價(jià)值函數(shù)g;
其中,iαref為輸出電流參考值的α分量;iβref為輸出電流參考值的β分量;λ為權(quán)重系數(shù);iα(k+2)為tk+2時(shí)刻輸出電流的α分量;iβ(k+2)為tk+2時(shí)刻輸出電流的β分量;P(k+2)為tk+2時(shí)刻有功功率預(yù)測值;Q(k+2)為tk+2時(shí)刻無功功率預(yù)測值;
S6,初始化,給定價(jià)值函數(shù)g的比較變量m,并給比較變量m和開關(guān)狀態(tài)Si賦初值;
S7,采集電網(wǎng)電壓ea、eb、ec,進(jìn)行Clark變換得到電網(wǎng)電壓的α分量eα和β分量eβ,并對電網(wǎng)電壓的α分量eα、電網(wǎng)電壓的β分量eβ分別進(jìn)行90°延遲,得到電網(wǎng)電壓α分量的90°延遲信號和電網(wǎng)電壓β分量的90°延遲信號;采集雙向交直流變換器的輸出電流ia、ib、ic并進(jìn)行Clark變換得到雙向交直流變換器輸出電流的α分量iα和β分量iβ,并對輸出電流的α分量iα和β分量iβ分別進(jìn)行90°延遲,得到輸出電流α分量的90°延遲信號和輸出電流β分量的90°延遲信號;
S8,結(jié)合步驟S2和S7計(jì)算當(dāng)前開關(guān)狀態(tài)下的雙向交直流變換器的輸出電壓Uj;
S9,結(jié)合步驟S3和步驟S8計(jì)算雙向交直流變換器的第一次輸出電流預(yù)測值和第一次功率預(yù)測值;
S10,結(jié)合步驟S3、步驟S8和步驟S9計(jì)算雙向交直流變換器的第二次輸出電流預(yù)測值和第二次功率預(yù)測值;
S11,結(jié)合步驟S4和步驟S7計(jì)算αβ靜止坐標(biāo)下的輸出電流參考值iαref和iβref;
S12,結(jié)合步驟S5、步驟S9和步驟S10計(jì)算價(jià)值函數(shù)g;
S13,比較價(jià)值函數(shù)g與比較變量m的大小,并將最小值賦值給比較變量m;
S14,判斷循環(huán)次數(shù)是否達(dá)到設(shè)定值,當(dāng)循環(huán)次數(shù)小于設(shè)定值時(shí),改變開關(guān)狀態(tài)值,重復(fù)步驟S7-S13;當(dāng)循環(huán)次數(shù)等于設(shè)定值時(shí),輸出最小價(jià)值函數(shù)g所對應(yīng)的輸出電壓矢量Uj;輸出電壓矢量Uj所對應(yīng)的開關(guān)狀態(tài)應(yīng)用于下一時(shí)刻,實(shí)現(xiàn)直接功率控制。
2.根據(jù)權(quán)利要求1所述的雙向交直流變換器模型預(yù)測電流控制方法,其特征在于:在步驟S2中,具體步驟為,S2.1,在abc三相靜止坐標(biāo)系下,獲取雙向交直流變換器的輸出電壓與開關(guān)狀態(tài)Si的計(jì)算公式,具體如下:
其中,Udc為直流母線電壓,uan為雙向交直流變換器的a相輸出電壓;ubn為雙向交直流變換器的b相輸出電壓;ucn為雙向交直流變換器的c相輸出電壓;Sa為a相的開關(guān)狀態(tài)值;Sb為b相的開關(guān)狀態(tài)值;Sc為c相的開關(guān)狀態(tài)值;
S2.2,對步驟S2.1中的公式2進(jìn)行Clark變換,得到αβ兩相靜止坐標(biāo)下雙向交直流變換器輸出電壓Uj與開關(guān)狀態(tài)Si的表達(dá)式,具體如下:
其中,uα為輸出電壓的α分量;uβ為輸出電壓的β分量;Udc為直流母線電壓,Sa為a相的開關(guān)狀態(tài)值;Sb為b相的開關(guān)狀態(tài)值;Sc為c相的開關(guān)狀態(tài)值。
3.根據(jù)權(quán)利要求1所述的雙向交直流變換器模型預(yù)測電流控制方法,其特征在于:在步驟S3中,具體步驟為,S3.1,根據(jù)基爾霍夫定律,得到并網(wǎng)逆變器在三相靜止坐標(biāo)系下的狀態(tài)方程;
其中,uan為雙向交直流變換器的a相輸出電壓;ubn為雙向交直流變換器的b相輸出電壓;ucn為雙向交直流變換器的c相輸出電壓;ia為雙向交直流變換器的a相輸出電流;ib為雙向交直流變換器的b相輸出電流;ic為雙向交直流變換器的c相輸出電流;ea為電網(wǎng)a相電壓;eb為電網(wǎng)b相電壓;ec為電網(wǎng)c相電壓;L為電感;R為電阻;
S3.2,對步驟S3.1中的公式4進(jìn)行Clark變換,得到αβ兩相靜止坐標(biāo)下的狀態(tài)方程;具體如下:
其中,L為電感;R為電阻;eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;iα為雙向交直流變換器的輸出電流的α分量;iβ為雙向交直流變換器的輸出電流的β分量;uα為輸出電壓的α分量;uβ為輸出電壓的β分量;
S3.3,對步驟S3.2中的公式5進(jìn)行離散化并整理,得到雙向交直流變換器在tk+1時(shí)刻預(yù)測電流:
式中,iα(k+1)為tk+1時(shí)刻輸出電流預(yù)測值的α分量;iβ(k+1)為tk+1時(shí)刻輸出電流預(yù)測值的β分量;iα(k)為tk時(shí)刻輸出電流的α分量;iβ(k)為tk時(shí)刻輸出電流的β分量;eα(k)為tk時(shí)刻電網(wǎng)電壓的α分量;eβ(k)為tk時(shí)刻電網(wǎng)電壓的β分量;uα(k)為tk時(shí)刻輸出電壓的α分量;uβ(k)為tk時(shí)刻輸出電壓的β分量;L為電感;R為電阻;Ts為采樣頻率;
S3.4,根據(jù)步驟S3.3中的公式6,得到tk+2時(shí)刻雙向交直流變換器在tk+2時(shí)刻預(yù)測電流;
式中,iα(k+1)為tk+1時(shí)刻輸出電流預(yù)測值的α分量;iβ(k+1)為tk+1時(shí)刻輸出電流預(yù)測值的β分量;iα(k+2)為tk+2時(shí)刻輸出電流的α分量;iβ(k+2)為tk+2時(shí)刻輸出電流的β分量;eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;uα(k+1)為tk+1時(shí)刻輸出電壓的α分量;uβ(k+1)為tk+1時(shí)刻輸出電壓的β分量;L為電感;R為電阻;Ts為采樣頻率;
S3.5,根據(jù)瞬時(shí)功率理論,得到電網(wǎng)側(cè)的有功功率p和無功功率q的計(jì)算公式,具體為:
式中:eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;iα為輸出電流的α分量;iβ為輸出電流的β分量;p為有功功率,q為無功功率;
S3.6,對于三相平衡電網(wǎng),當(dāng)采樣頻率Ts較高時(shí),有:
S3.7,將步驟S3.6中的公式9代入步驟S3.5的公式8中,得到tk+1時(shí)刻雙向交直流變換器的功率預(yù)測模型:
式中,P(k+1)為tk+1時(shí)刻有功功率預(yù)測值;Q(k+1)為tk+1時(shí)刻無功功率預(yù)測值;P(k)為tk時(shí)刻有功功率預(yù)測值;Q(k)為tk時(shí)刻無功功率預(yù)測值;iα(k)為tk時(shí)刻輸出電流預(yù)測值的α分量;iβ(k)為tk時(shí)刻輸出電流預(yù)測值的β分量;eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;uα(k)為tk時(shí)刻輸出電壓的α分量;uβ(k)為tk時(shí)刻輸出電壓的β分量;L為電感;Ts為采樣頻率;
S3.8,根據(jù)步驟S3.7中的公式10得到tk+2時(shí)刻并網(wǎng)逆變器與輸出電壓Uj有關(guān)的功率預(yù)測模型;具體為:
其中,P(k+1)為tk+1時(shí)刻有功功率預(yù)測值;Q(k+1)為tk+1時(shí)刻無功功率預(yù)測值;P(k+2)為tk+2時(shí)刻有功功率預(yù)測值;Q(k+2)為tk+2時(shí)刻無功功率預(yù)測值;iα(k+1)為tk+1時(shí)刻輸出電流預(yù)測值的α分量;iβ(k+1)為tk+1時(shí)刻輸出電流預(yù)測值的β分量;eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;uα(k+1)為tk+1時(shí)刻輸出電壓的α分量;uβ(k+1)為tk+1時(shí)刻輸出電壓的β分量;L為電感;Ts為采樣頻率。
4.根據(jù)權(quán)利要求1所述的雙向交直流變換器模型預(yù)測電流控制方法,其特征在于:在步驟S4中,具體步驟為,S4.1,在不平衡電網(wǎng)下,分別計(jì)算電網(wǎng)電壓e、輸出電流i的正序分量和負(fù)序分量;
式中:ω為dq坐標(biāo)系旋轉(zhuǎn)角速度,為電網(wǎng)電壓在αβ坐標(biāo)系的正序分量;為電網(wǎng)電壓在αβ坐標(biāo)系的負(fù)序分量;為電網(wǎng)電壓在dq坐標(biāo)系的正序分量;為電網(wǎng)電壓在dq坐標(biāo)系的負(fù)序分量;為輸出電流在αβ坐標(biāo)系的正序分量;為輸出電流在αβ坐標(biāo)系的負(fù)序分量;為輸出電流在dq坐標(biāo)系的正序分量;為輸出電流在dq坐標(biāo)系的負(fù)序分量;ed+為電網(wǎng)電壓在dq坐標(biāo)系的d軸正序分量數(shù)值;eq+為電網(wǎng)電壓在dq坐標(biāo)系的q軸正序分量數(shù)值;ed-為電網(wǎng)電壓在dq坐標(biāo)系的d軸負(fù)序分量數(shù)值;eq-為電網(wǎng)電壓在dq坐標(biāo)系的q軸負(fù)序分量數(shù)值;id+為輸出電流在dq坐標(biāo)系的d軸正序分量數(shù)值;iq+為輸出電流在dq坐標(biāo)系的q軸正序分量數(shù)值;id-為輸出電流在dq坐標(biāo)系的d軸負(fù)序分量數(shù)值;iq-為輸出電流在dq坐標(biāo)系的q軸負(fù)序分量數(shù)值;
S4.2,獲得dq坐標(biāo)下的有功功率和無功功率與正負(fù)序分量之間的關(guān)系式;
具體步驟為:S4.2.1,根據(jù)瞬時(shí)功率理論,電網(wǎng)側(cè)功率表示如下:
S=ei*=p+jq (14);
式中:
其中,p為有功功率,q為無功功率;p0為有功功率的基準(zhǔn)值;pc2為有功功率的余弦脈動(dòng)分量;ps2為有功功率的正弦脈動(dòng)分量;q0為無功功率的基準(zhǔn)值;qc2為無功功率的余弦脈動(dòng)分量;qs2為無功功率的正弦脈動(dòng)分量;
S4.2.2,將步驟S4.1中公式12和公式13代入步驟S4.2.1中的公式15,計(jì)算整理,得到dq坐標(biāo)下的有功功率和無功功率與正負(fù)序分量之間的關(guān)系式:
式中:p0為有功功率的基準(zhǔn)值;pc2為有功功率的余弦脈動(dòng)分量;ps2為有功功率的正弦脈動(dòng)分量;q0為無功功率的基準(zhǔn)值;qc2為無功功率的余弦脈動(dòng)分量;qs2為無功功率的正弦脈動(dòng)分量;ed+為電網(wǎng)電壓在dq坐標(biāo)系的d軸正序分量數(shù)值;eq+為電網(wǎng)電壓在dq坐標(biāo)系的q軸正序分量數(shù)值;ed-為電網(wǎng)電壓在dq坐標(biāo)系的d軸負(fù)序分量數(shù)值;eq-為電網(wǎng)電壓在dq坐標(biāo)系的q軸負(fù)序分量數(shù)值;id+為輸出電流在dq坐標(biāo)系的d軸正序分量數(shù)值;iq+為輸出電流在dq坐標(biāo)系的q軸正序分量數(shù)值;id-為輸出電流在dq坐標(biāo)系的d軸負(fù)序分量數(shù)值;iq-為輸出電流在dq坐標(biāo)系的q軸負(fù)序分量數(shù)值;
S4.3,得到在αβ靜止坐標(biāo)系下,有功功率p、無功功率q與電網(wǎng)電壓、輸出電流及電網(wǎng)電壓的90°延遲信號、輸出電流的90°延遲信號的關(guān)系式。
具體步驟如下:S4.3.1,在αβ靜止坐標(biāo)系下,計(jì)算90°延遲信號與正負(fù)序分量之間的關(guān)系:
假設(shè)αβ靜止坐標(biāo)系下的變量為x,則其90°延遲信號表示為x′,延遲信號與正負(fù)序分量之間的關(guān)系為:
x′=xαβ+′+xαβ-′=-jxαβ++jxαβ- (17);
則x、x′與正負(fù)序分量的關(guān)系表示為:
S4.3.2,對步驟S4.3.1中的公式18求逆可得:
整理后,得到dq旋轉(zhuǎn)坐標(biāo)系和αβ靜止坐標(biāo)系的正負(fù)序分量之間的關(guān)系為:
S4.3.3,結(jié)合步驟S4.3.2中的公式19和公式20,得到dq坐標(biāo)系下正負(fù)序分量與αβ坐標(biāo)下變量及延遲信號之間的表達(dá)式:
S4.3.4,將步驟S4.3.3中的公式21代入步驟S4.2中的公式16,得到dq坐標(biāo)下的有功功率和無功功率與正負(fù)序分量之間的關(guān)系式:
其中:
式中:iα為輸出電流的α分量;iβ為輸出電流的β分量;iα′為輸出電流α分量的90°延遲信號;iβ′為輸出電流β分量的90°延遲信號;eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;eα′為電網(wǎng)電壓α分量的90°延遲信號;eβ′為電網(wǎng)電壓β分量的90°延遲信號。
S4.4,為消除雙向交直流變換器功率的二次脈動(dòng),分為有功功率的二次脈動(dòng)和無功功率的二次脈動(dòng);
為消除有功功率的二次脈動(dòng),實(shí)現(xiàn)雙向交直流變換器有功功率的穩(wěn)定輸出,令:
根據(jù)步驟S4.3中的公式22和公式23求解公式24,得到輸出電流參考值與有功功率、電網(wǎng)電壓的α分量、β分量及延遲信號之間的表達(dá)式:
式中,iaref為輸出電流參考值的α分量;iβref為輸出電流參考值的β分量;eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;eα′為電網(wǎng)電壓α分量的90°延遲信號;eβ′為電網(wǎng)電壓β分量的90°延遲信號;pref為有功功率給定值;
為消除無功功率的二次脈動(dòng),實(shí)現(xiàn)雙向交直流變換器無功功率的穩(wěn)定輸出,令:
根據(jù)步驟S4.3中的公式22和公式23求解公式26,輸出電流參考值與無功功率、電網(wǎng)電壓的α分量、β分量及延遲信號之間的表達(dá)式:
式中,iαref為輸出電流參考值的α分量;iβref為輸出電流參考值的β分量;eα為電網(wǎng)電壓的α分量;eβ為電網(wǎng)電壓的β分量;eα′為電網(wǎng)電壓α分量的90°延遲信號;eβ′為電網(wǎng)電壓β分量的90°延遲信號;qref為有功功率給定值。