• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    考慮風(fēng)能隨機(jī)性的V S C - H V D C風(fēng)電并網(wǎng)優(yōu)化方法研究

    2016-03-30 05:45:40于永生馮延暉江紅鑫邱穎寧南京理工大學(xué)能源與動力工程學(xué)院江蘇南京210094
    電網(wǎng)與清潔能源 2016年1期
    關(guān)鍵詞:單純形法風(fēng)電場遺傳算法

    于永生,馮延暉,江紅鑫,邱穎寧(南京理工大學(xué)能源與動力工程學(xué)院,江蘇南京 210094)

    ?

    考慮風(fēng)能隨機(jī)性的V S C - H V D C風(fēng)電并網(wǎng)優(yōu)化方法研究

    于永生,馮延暉,江紅鑫,邱穎寧
    (南京理工大學(xué)能源與動力工程學(xué)院,江蘇南京210094)

    KEY W0RDS:VSC-HVDC;wind farm;genetic a1gorithm;simP1ex method;random wind

    摘要:隨著能源危機(jī)的日益臨近和電力電子器件技術(shù)的大幅進(jìn)步與應(yīng)用,風(fēng)電場經(jīng)VSC-HVDC并網(wǎng)已經(jīng)變成了一種可能。相對于其他諸如火力發(fā)電和水利發(fā)電等傳統(tǒng)發(fā)電方式,風(fēng)力發(fā)電不能夠控制能量的注入與供給,所以風(fēng)力發(fā)電的輸出穩(wěn)定與否主要取決于風(fēng)速的變化。建立了風(fēng)電場等效聚合模型和柔性直流輸電(VSC-HVDC)的詳細(xì)模型。為了平抑風(fēng)電并網(wǎng)時風(fēng)電場的功率波動,分別采用了單純形法和遺傳算法對柔性直流輸電的控制器進(jìn)行了優(yōu)化,優(yōu)化結(jié)果顯示,兩種優(yōu)化方法均能平抑風(fēng)電場的輸出功率波動。

    關(guān)鍵詞:柔性直流輸電;風(fēng)電場;遺傳算法;單純形法;隨機(jī)風(fēng)

    以風(fēng)力發(fā)電為代表的新能源發(fā)電技術(shù)與傳統(tǒng)的發(fā)電型式如火力發(fā)電的區(qū)別之一就是能量來源的不可控性,這就造成了風(fēng)電場出力的不穩(wěn)定性,表現(xiàn)在風(fēng)電場輸出功率穩(wěn)定性的好壞取決于風(fēng)能資源的穩(wěn)定與否。基于IGBT的新型輸電方式——柔性直流輸電,以其高度可控性得到了風(fēng)電并網(wǎng)領(lǐng)域內(nèi)眾多專家學(xué)者和工程師廣泛關(guān)注[1]。

    文獻(xiàn)[2]對于采用雙饋機(jī)組的風(fēng)電場經(jīng)柔性直流輸電并網(wǎng)設(shè)計了一種新型控制方式,并通過風(fēng)電場波動和負(fù)荷波動進(jìn)行仿真驗證該控制方式的正確性。文獻(xiàn)[3]對柔性直流輸電的風(fēng)電場側(cè)變流器設(shè)計了一種新的交流電壓—功角控制方式,能夠?qū)崿F(xiàn)對波動風(fēng)速的同步傳送。文獻(xiàn)[4]對于兩端VSC-HVDC采用單純形算法進(jìn)行了優(yōu)化,仿真結(jié)果顯示單純形算法對于VSC-HVDC優(yōu)化有重要的作用,但該文所做研究是基于間接電流控制。間接電流控制動態(tài)響應(yīng)慢,工程實(shí)際大多采用響應(yīng)速度快、魯棒性好的直接電流控制。文獻(xiàn)[5]采用遺傳算法對柔性直流輸電控制器參數(shù)進(jìn)行優(yōu)化,優(yōu)化后的系統(tǒng)響應(yīng)特性較優(yōu)化前有較明顯改善。文獻(xiàn)[6]采用單純性法對風(fēng)電場經(jīng)多端柔性直流輸電進(jìn)行了優(yōu)化,但文獻(xiàn)[6]僅考慮到了穩(wěn)態(tài)情況下的響應(yīng)情況,并沒有對隨機(jī)風(fēng)況下的柔性直流輸電模型進(jìn)行分析討論。

    鑒于大多數(shù)文獻(xiàn)都沒有涉及到風(fēng)電場和柔性直流輸電系統(tǒng)之間的相互作用,考慮到隨機(jī)風(fēng)速,即風(fēng)速隨機(jī)變化情況下風(fēng)電場的輸出特性,本文對于風(fēng)電場在隨機(jī)風(fēng)條件下輸出特性進(jìn)行了討論,并采用遺傳算法和單純形法兩種啟發(fā)式優(yōu)化算法分別對柔性直流輸電控制器參數(shù)進(jìn)行了優(yōu)化,最后將2種優(yōu)化方法的優(yōu)化結(jié)果進(jìn)行了比較分析。

    1 風(fēng)電場模型

    由于技術(shù)的限制和風(fēng)力發(fā)電的特性,與其他傳統(tǒng)的大型發(fā)電機(jī)組不同,風(fēng)機(jī)的單機(jī)容量目前最大只有幾個MW,一個風(fēng)電場中風(fēng)力發(fā)電機(jī)組數(shù)量眾多,從幾十臺到上百臺不等,如果所有風(fēng)力發(fā)電機(jī)均采用發(fā)電機(jī)的全暫態(tài)模型來仿真,這樣的仿真模型是一種很復(fù)雜的高階模型[10]。仿真分析時,計算時間長,內(nèi)存需求大,對計算機(jī)的配置要求較高。事實(shí)上,研究某個風(fēng)電場對電力系統(tǒng)的影響,可以將整個風(fēng)電場看作一個整體,并不需要考慮每一臺風(fēng)力發(fā)電機(jī)對電網(wǎng)的影響。因此在大規(guī)模風(fēng)電場并網(wǎng)仿真計算中,不適合也沒有必要采用詳細(xì)模型來研究整個風(fēng)電場對電網(wǎng)的動態(tài)性能的影響。研究電網(wǎng)和風(fēng)電場通過輸電系統(tǒng)之間的交互作用,可以將相同類型的風(fēng)力發(fā)電機(jī)等效為一臺機(jī)組[10-11]。對于有上百臺機(jī)組的風(fēng)電場,并網(wǎng)分析中采用這種聚合模型是很有必要的。在對于含風(fēng)電場的電力系統(tǒng)暫態(tài)分析時,可以通過這種化簡等值,降低仿真模型階數(shù),節(jié)約計算時間[13]。

    將風(fēng)電場等效為一臺發(fā)電機(jī)的仿真方法的精度雖然略低于將風(fēng)電場所有發(fā)電機(jī)的數(shù)學(xué)模型全部代入仿真模型中,但這種等效的方法仍然可以正確反映整個風(fēng)電場在并網(wǎng)點(diǎn)的動態(tài)特性,顯著減少了計算時間,具有很高的實(shí)用性[12]。本文將風(fēng)電場等效為同一個聚合模型,風(fēng)電場的聚合模型包含了風(fēng)速模型,風(fēng)力機(jī)的機(jī)械部分,等效的機(jī)組、母線、變壓器等電氣元件和等效風(fēng)速[12]。引入等效風(fēng)速是為了研究風(fēng)能資源的特性對于風(fēng)電場并網(wǎng)功率的影響。

    對于等效風(fēng)速的計算方法,文獻(xiàn)[12]將整個風(fēng)場機(jī)組功率求取平均值后,通過等效風(fēng)功率曲線逆函數(shù)求取等效風(fēng)速。

    風(fēng)力機(jī)模型

    式中:Pw為風(fēng)輪從風(fēng)中捕獲的風(fēng)能;p為空氣密度,kg/m3;CP為風(fēng)力機(jī)的風(fēng)能轉(zhuǎn)換效率系數(shù);AR為風(fēng)輪掃過的面積;vw為風(fēng)速,m/s;λ為葉尖速比;β為葉片槳距角,°。

    目前風(fēng)力發(fā)電的研究主要采用風(fēng)速的四分量模型[21]。為了準(zhǔn)確地模擬風(fēng)電場風(fēng)速對風(fēng)電場出力的影響,以歐洲某個風(fēng)電場的風(fēng)速統(tǒng)計為基準(zhǔn),按四分量模型進(jìn)行擬合,得到一個風(fēng)速的隨機(jī)數(shù)序列,作為風(fēng)電場風(fēng)速輸入。

    考慮到風(fēng)電場需要能夠反映風(fēng)電場風(fēng)機(jī)出力變化對風(fēng)電場輸出功率的影響,本文用一臺大容量風(fēng)力發(fā)電機(jī)來代表整個風(fēng)電場的出力狀況[10]。該等值機(jī)組在等效風(fēng)速作用下所產(chǎn)生的風(fēng)電功率與所有機(jī)組輸出功率相等。等效風(fēng)力發(fā)電機(jī)采用直驅(qū)型永磁同步發(fā)電機(jī)(PMSG)。由于實(shí)際中風(fēng)力機(jī)轉(zhuǎn)速很慢,需要安裝齒輪箱等變速環(huán)節(jié)。本文主要研究風(fēng)電場風(fēng)速變化對風(fēng)電場輸出的影響,為了簡化,本文省略了齒輪箱環(huán)節(jié),即風(fēng)力機(jī)中齒輪變比為1。風(fēng)力發(fā)電機(jī)的極對數(shù)設(shè)置為100,此時在額定恒風(fēng)速下運(yùn)行風(fēng)力發(fā)電機(jī)輸出電壓和電流頻率為50 Hz。變流器和風(fēng)電場以及變流器和電網(wǎng)之間的聯(lián)結(jié)處含有等效變壓器以及聯(lián)結(jié)電抗和等效電阻。

    風(fēng)電場模型如圖1所示,四分量風(fēng)速模型輸出為風(fēng)速。風(fēng)力機(jī)模型輸入為風(fēng)速和永磁同步發(fā)電機(jī)的轉(zhuǎn)速,輸出為轉(zhuǎn)矩。風(fēng)力機(jī)槳距角控制模型用來調(diào)節(jié)風(fēng)力機(jī)的槳距角。

    2 柔性直流輸電模型

    2.1換流器數(shù)學(xué)模型及其控制模型

    目前工程上投入應(yīng)用較多的柔性直流輸電(VSC-HVDC)主要采用三相兩電平結(jié)構(gòu)[8]。

    由Kirchhoff定律,可得如下暫態(tài)方程

    圖1 PSCAD中的風(fēng)電場模型Flg. 1 Wlnd farm model ln PSCAD

    圖2 柔性直流輸電換流站模型Flg. 2 VSC-HVDC system dlagram

    式中:ia,ib,ic為三相電流;Usa,Usb,Usc為風(fēng)電場變壓器一次側(cè)三相電壓;Uca,Ucb,Ucc為變流器交流側(cè)基波三相電壓。式(2)經(jīng)過dq變換后為:

    直流側(cè)動態(tài)方程:

    式中:ir為電容和變流器之間直流線路的電流;iL為電容右側(cè)直流線路的電流;Ur為直流電壓,該公式描述直流電容流入流出電流與電容電壓的關(guān)系。

    根據(jù)瞬時無功功率理論[8],abc三相靜止坐標(biāo)系下?lián)Q流站與交流系統(tǒng)交換的有功功率ps和無功功率qs為:

    經(jīng)過dq變換,在旋轉(zhuǎn)坐標(biāo)系中,風(fēng)電場輸出的有功功率ps和無功功率qs為:

    當(dāng)d軸以電網(wǎng)電壓向量定位時,即Usq=0,則式(6)變?yōu)?/p>

    在式(7)中,分別控制id和iq可以調(diào)節(jié)換流站與交流系統(tǒng)交換的有功和無功功率。

    對于直接電流控制,令:

    又有

    當(dāng)標(biāo)量函數(shù)V.L<0時,系統(tǒng)處在穩(wěn)定平衡點(diǎn),ur能夠跟蹤urref,因此可以得出控制公式:

    即為:

    在大多數(shù)情況下,u.rref=0,即可得到直流電壓PI控制器。同理可以推導(dǎo)出其他控制量的控制器。

    采用前饋解耦控制,可以得到電網(wǎng)側(cè)換流站(G-VSC)的電流內(nèi)環(huán)控制和風(fēng)電場側(cè)換流站(WFVSC)的電流內(nèi)環(huán)控制。為了減小風(fēng)電場側(cè)換流站輸出波動,風(fēng)電場側(cè)換流站采用改進(jìn)的電流內(nèi)環(huán)控制[17]。

    圖3 風(fēng)電場側(cè)電流內(nèi)環(huán)控制Flg. 3 Inner current control of the WF-VSC

    圖4 電網(wǎng)側(cè)電流內(nèi)環(huán)控制Flg. 4 Inner current control of G-VSC

    內(nèi)環(huán)電流控制的傳遞函數(shù)[19]為:

    式中:(kP+kis)為比例積分環(huán)節(jié);1/(R+sXL/ω)為等效電感和電阻的傳遞函數(shù);KPWM/(1+(1/2)Tss)為變流器傳遞函數(shù);KPWM為PWM的等效增益,KPWM=Udc/(2urms),urms為交流線電壓有效值。1/(1+Tss)為采樣濾波環(huán)節(jié)。

    本文所采用模型通過外環(huán)控制來區(qū)別不同換流站的控制職能。在連接有風(fēng)電場的多端系統(tǒng)中,由于風(fēng)電場無法提供一個穩(wěn)定的電壓,所以離線仿真時需要大電網(wǎng)先并入網(wǎng)絡(luò)給柔性直流輸電系統(tǒng)提供一個穩(wěn)定的直流電壓,然后風(fēng)電場才能并入網(wǎng)絡(luò),否則容易引起風(fēng)電場輸出的波動[7]。

    圖5 外環(huán)電壓控制器Flg. 5 Voltage outer control

    外環(huán)電壓控制的傳遞函數(shù)[19]為:

    式中:(kP+kis)為比例積分環(huán)節(jié);1/(1+Tss)為采樣濾波環(huán)節(jié);usq為交流側(cè)三相電壓經(jīng)dq變換得到的q軸電壓;Udc為直流電壓;C為直流側(cè)電容;1/(1+3Tss)為電流內(nèi)環(huán)的等效傳遞函數(shù)。

    圖6 外環(huán)有功控制器Flg. 6 Actlve power outer control

    圖7 外環(huán)無功控制器Flg. 7 Reactlve power outer control

    外環(huán)有功和無功的傳遞函數(shù)[19]為:

    式中:(kP+kis)為比例積分環(huán)節(jié);1/(1+Tss)為采樣濾波環(huán)節(jié);usq為交流側(cè)三相電壓經(jīng)dq變換得到的q軸電壓;Udc為直流電壓,1/(1+3Tss)為電流內(nèi)環(huán)的等效傳遞函數(shù)。

    2.2電纜線路模型

    電纜線路暫態(tài)模型[20]如圖8所示。

    圖8 電纜線路的分布模型Flg. 8 Dlstrlbutlon model of the cable llne

    式(22)為電纜線路的微分方程組。由于該方程組在電纜線路中很難找到解析解。在含有多導(dǎo)體的同軸電纜中,這個模型會更加復(fù)雜。因此,大多數(shù)情況下主要采用行波模型。行波模型如式(23)所示。其中ω為輸電線傳輸正弦波角頻率,χ和t分別為傳播距離和時間。

    省略對時間的導(dǎo)數(shù),可以得到僅關(guān)于距離χ的微分方程組:

    令Z0=R0+jωL0,Y0=G0+jωC0。對于含多導(dǎo)體的同軸電纜來說:

    該方程的通解即為行波模型:

    由于直流線路不同于交流線路,直流線路沒有基波頻率ω,所以為了降低模型的復(fù)雜度,電纜線路模型采用Bergeron模型,Bergeron模型基于行波理論,Bergeron是一連串π形等效電路,該模型能夠準(zhǔn)確反映電纜線路的電感和電容的分布特性。Bergeron模型的電阻參數(shù)基于半集總模型,電阻值有50%集中在電纜中心點(diǎn),另外兩端各占25%的阻值。

    表1 電纜線路參數(shù)Tab. 1 Cable parameters

    3 單純形法

    假設(shè)搜索空間中有m個點(diǎn),每個點(diǎn)包括d個分量,記X={χ1,χ2,…,χd},其中m=d+l。將m個頂點(diǎn)按照其對應(yīng)的目標(biāo)函數(shù)值從小到大排序(優(yōu)劣排序),第i個頂點(diǎn)記作Xm,i,m個頂點(diǎn)的單純形Sm表示為:S={Xm,1,Xm,2,…,Xm,m}。單純形Sm的最差頂點(diǎn)記為Xm,w,Xm,w=Xm,m;最好頂點(diǎn)記作Xm,b,Xm,b=Xm,1;除最差頂點(diǎn)Xm,w之外的其他所有頂點(diǎn)的形心記作Xcm,w;除第i個頂點(diǎn)Xm,i之外的其他所有頂點(diǎn)的形心記作Xcm,i。

    單純形法的搜索步驟為:

    1)給定初始頂點(diǎn)=X'm,0=[χ0,1,χ0,2,…,χ0,m-1]

    則單純形其他初始頂點(diǎn)為:

    2)a)檢查收斂準(zhǔn)則,若滿足則終止優(yōu)化。b)計算各頂點(diǎn)的目標(biāo)函數(shù),根據(jù)目標(biāo)函數(shù)值從小到大順序排序(優(yōu)劣順序)。排序后各頂點(diǎn)依次為Xm,1,Xm,2,…,Xm,m;確定最好頂點(diǎn)Xm,b、最差頂點(diǎn)Xm,w,計算形心設(shè)置一次成功搜索標(biāo)志變量F1ag=F。

    3)單純法尋優(yōu)過程:

    a)X*=Xm,w;F1ag=F;

    b)計算反射點(diǎn)

    及其對應(yīng)的目標(biāo)函數(shù)之f(Xrm,w);

    c)如果f(Xrm,w)<f(Xm,b),則計算擴(kuò)張點(diǎn)

    及其目標(biāo)函數(shù)值f(Xem,w),

    如果f(Xem,w)<f(Xm,b),則

    否則X*=Xrm,w;

    d)如果f(Xrm,w)<f(Xm,m-1),則

    e)如果

    則計算內(nèi)壓縮點(diǎn)

    及其目標(biāo)函數(shù)f(Xicm,w),

    如果f(Xicm,w)<f(Xm,m),則

    否則X*=Xrm,w;

    f)如果f(Xrm,w)≥f(Xm,m),

    則計算外部壓縮點(diǎn)

    及其目標(biāo)函數(shù)值f(Xocm,w);如果

    則X*=Xocm,w,F(xiàn)1ag=T;

    g)返回計算結(jié)果F1ag和X*;

    4)if(F1ag=T){Xm,w=X*;Goto(2)}

    e1se{

    a)單純性各頂點(diǎn)以最好頂點(diǎn)Xm,b為中心進(jìn)行整體壓縮,即

    圖9 單純形法優(yōu)化過程流程圖Flg. 9 Flowchart of the slmplex method

    b)Goto(2)

    由于目標(biāo)函數(shù)是個具有凸性的超平面,故它的一個最優(yōu)值(若存在的話)必在單純形的某個頂點(diǎn)上得到[14]。

    4 遺傳算法

    遺傳算法作為一種啟發(fā)式算法,模擬“物競天擇,適者生存”的生物進(jìn)化理論,通過維護(hù)一個包含潛在解的群體執(zhí)行多方向的搜索。遺傳算法的步驟如下[16]:

    1)對優(yōu)化問題的解進(jìn)行編碼。按生物學(xué)術(shù)語,我們稱一個解的編碼為一個染色體,組成編碼的元素成為基因。

    2)適應(yīng)度函數(shù)的構(gòu)造和應(yīng)用,適應(yīng)度函數(shù)一般是依據(jù)優(yōu)化問題的目標(biāo)函數(shù)而定。適度函數(shù)對應(yīng)于生物學(xué)中的環(huán)境,適度函數(shù)是用來淘汰掉劣等基因,即不滿足適度函數(shù)的染色體將被淘汰掉。

    3)(隨機(jī))產(chǎn)生初始種群并計算各染色體的適應(yīng)度。

    4)父母雙親的選擇。衡量標(biāo)準(zhǔn)就是染色體的適應(yīng)度。從生物學(xué)的角度講,后代具有優(yōu)良特性的先決條件是選出的雙親具有相對優(yōu)良的特性。不滿足給定條件的個體不能作為下一代的雙親,即被淘汰。

    5)繁衍后代。應(yīng)遵循兩條自然法則:父母的優(yōu)良特性得以遺傳、后代能夠產(chǎn)生新的優(yōu)良特性。這可通過交叉和變異的協(xié)同作用來實(shí)現(xiàn)。

    6)群體的換代。重復(fù)第4,5步直到新的群體全部產(chǎn)生,并計算各染色體的適應(yīng)度。

    7)進(jìn)化過程終止判別。依據(jù)某種收斂準(zhǔn)則進(jìn)行判斷,該收斂準(zhǔn)則可以以目標(biāo)函數(shù)的誤差值,也可以用固定代數(shù)作為判別條件。如果滿足則終止進(jìn)化,否則重復(fù)第4~7步進(jìn)行下一輪迭代。

    本文所采用的遺傳的參數(shù)為:群體規(guī)模為100,交叉率為0.8,變異率為0.1,選擇方式為隨機(jī)選擇,以40代作為收斂依據(jù)。

    遺傳算法并不能保證獲得問題的最優(yōu)解,但是使用遺傳算法的最大優(yōu)點(diǎn)在于你不必去了解和操心如何去尋找最優(yōu)解。而只要簡單地淘汰掉一些表現(xiàn)不好的個體就行了。這樣,經(jīng)過若干次淘汰,剩下的個體即為滿足條件的優(yōu)良個體。

    5 仿真結(jié)果與分析

    由ITAE最小化準(zhǔn)則[15],本文設(shè)置優(yōu)化系統(tǒng)的目標(biāo)函數(shù)為:

    式中:eP_rec為海上換流站(整流側(cè))有功功率給定值與測量值之間的差值;eP_inv為陸上換流站(逆變側(cè))有功功率給定值與測量值之間的差值;eUdc為直流側(cè)電壓給定值與測量值之間的差值。

    由圖10和圖11我們可以看出柔性直流輸電輸送功率在趨勢上跟隨風(fēng)速變化。由于風(fēng)速的大幅波動導(dǎo)致能量來源的大幅波動,不可能使得風(fēng)電場的輸出波動為零。本文所做工作即在風(fēng)速隨機(jī)變化的情況下使得柔性直流輸電的功率波動盡量減小。

    圖10 風(fēng)速Flg. 10 Wlnd speed

    圖11 風(fēng)電場側(cè)有功功率Flg. 11 The actlve power ln WF-VSC

    由圖11可以看出,在系統(tǒng)啟動階段,采用遺傳算法優(yōu)化的控制器能夠使得整個系統(tǒng)的響應(yīng)速度優(yōu)于單純形法的優(yōu)化結(jié)果。系統(tǒng)啟動階段的最大超調(diào)量98.254小于系統(tǒng)初始值運(yùn)行狀態(tài)下的最大超調(diào)量139.463。而單純形算法優(yōu)化后的結(jié)果沒有超調(diào)量,屬于過阻尼系統(tǒng),響應(yīng)速度太慢。同時,采用遺傳算法優(yōu)化后的控制器輸出功率在0.562 s后即接近風(fēng)電場的額定功率75 MW,而采用單純形法優(yōu)化后的控制器則需要在1.373 s后才接近風(fēng)電場的額定輸出功率。

    采用初始參數(shù)的控制器,柔性直流輸電整流側(cè)輸入功率1~3 s之間有功功率標(biāo)準(zhǔn)差為13.316 7。整流側(cè)控制器采用單純形優(yōu)化后,整流側(cè)輸入功率1~3 s之間有功功率標(biāo)準(zhǔn)差為4.451 7。整流側(cè)控制器采用遺傳算法優(yōu)化后,整流側(cè)輸入功率1~3 s之間有功功率標(biāo)準(zhǔn)差為9.357。

    由于柔性直流輸電逆變側(cè)控制柔性直流輸電直流側(cè)電壓,因此在系統(tǒng)啟動階段,交流側(cè)會向直流側(cè)輸送大量的有功功率,如圖12所示。此時的初始控制器參數(shù)下的反向功率超調(diào)量達(dá)到367.278 MW。該反向超調(diào)量可以通過增設(shè)柔性直流輸電啟動電阻來消除,該問題不在本文討論范圍。

    圖12 逆變側(cè)有功功率Flg. 12 The actlve power ln G-VSC

    由圖12可以看出,在系統(tǒng)啟動階段,采用遺傳算法優(yōu)化的系統(tǒng)的響應(yīng)與單純形類似,但系統(tǒng)的穩(wěn)態(tài)誤差要比單純形小,而在超調(diào)量指標(biāo)遺傳算法要小于控制系統(tǒng)初始參數(shù)的動態(tài)性能。

    采用初始參數(shù)的控制器,柔性直流輸電逆變側(cè)輸出功率1~3 s之間有功功率的標(biāo)準(zhǔn)差為9.673 3。逆變側(cè)控制器采用單純形優(yōu)化后,逆變側(cè)輸出功率1~3 s之間有功功率的標(biāo)準(zhǔn)差為3.714。逆變側(cè)控制器采用遺傳算法優(yōu)化后,逆變側(cè)輸出功率1~3 s之間有功功率的標(biāo)準(zhǔn)差為7.956 5。

    圖13 直流電壓Flg. 13 DC Voltage

    由于柔性直流輸電系統(tǒng)的直流電壓是由逆變側(cè)來來控制,而逆變側(cè)系統(tǒng)連接的是無窮大電網(wǎng),所以柔性直流輸電直流側(cè)電壓不會隨風(fēng)電場的風(fēng)速波動而波動。由于優(yōu)化的目標(biāo)函數(shù)中包含直流電壓一項,所以直流電壓的最大超調(diào)量由優(yōu)化前的233.452降低至176.425(遺傳算法)和151.63(單純形法)。而優(yōu)化后的系統(tǒng)中直流電壓的震蕩次數(shù)也明顯減小。

    由表2可以看出,無論在何種控制器參數(shù)下,逆變側(cè)有功功率的波動都小于整流側(cè)有功功率的波動,即柔性直流輸電本身的特性即能夠平抑風(fēng)電場的輸出波動,起到削峰填谷的效果,具有魯棒性,適合風(fēng)力發(fā)電并網(wǎng)。為了便于比較,本文中所做仿真采用同一個隨機(jī)風(fēng)速,優(yōu)化后的參數(shù)能夠使得風(fēng)力發(fā)電機(jī)的輸出波動減小,對于波動方差這一指標(biāo)來說,單純形法的優(yōu)化結(jié)果要優(yōu)于遺傳算法。

    表2 不同優(yōu)化結(jié)果的輸出方差Tab. 2 0utput varlance of dlfferent optlmlzatlon results

    所做仿真中穩(wěn)定額定風(fēng)速情況下風(fēng)電場有功功率額定值為75 MW。由表3可以看出,單純形法優(yōu)化后的結(jié)果均值偏離額定值較多,即單純形優(yōu)化后的控制系統(tǒng)使得風(fēng)力發(fā)電機(jī)經(jīng)柔性直流輸電并入電網(wǎng)的能量相對減少,而遺傳算法優(yōu)化后的結(jié)果偏離額定值較小。由于積分時間常數(shù)是用來減小系統(tǒng)的穩(wěn)態(tài)誤差,這與表4中優(yōu)化后的積分時間常數(shù)相一致。同時,表3中整流側(cè)有功功率和逆變側(cè)有功功率之間的差值是由換流站,輸電線路以及變壓器的損耗造成的。

    表3 不同優(yōu)化結(jié)果的輸出均值Tab. 3 0utput mean values of dlfferent optlmlzatlon results

    表4 不同優(yōu)化結(jié)果的參數(shù)值Tab. 4 Parameters of dlfferent optlmlzatlon results

    表4中,P1和I1是柔性直流輸電風(fēng)電場側(cè)有功功率控制的PI控制器參數(shù),P3和I3是柔性直流輸電電網(wǎng)側(cè)直流電壓PI控制器的控制參數(shù),P2,I2和P4,I4分別是整流和逆變側(cè)的電流內(nèi)環(huán)控制器參數(shù)。

    為了測試該系統(tǒng)在故障情況下的動態(tài)響應(yīng),本文在系統(tǒng)仿真中的3 s時,在系統(tǒng)逆變側(cè)交流電網(wǎng)引入三相接地短路故障,故障持續(xù)時間0.05 s。

    如圖14所示,在逆變側(cè)發(fā)生故障后,整流側(cè)有功功率輸出迅速下跌,當(dāng)柔性直流輸電系統(tǒng)控制器采用初始值參數(shù)時,系統(tǒng)發(fā)生震蕩,故障切除后最大超調(diào)量達(dá)到160 MW,系統(tǒng)失去穩(wěn)定(恒定風(fēng)速則不會失去穩(wěn)定),而參數(shù)優(yōu)化后,系統(tǒng)能夠在故障切除后跟隨風(fēng)速變化,與圖11穩(wěn)態(tài)情況相一致。采用遺傳算法優(yōu)化后的控制器在發(fā)生故障后波動較為頻繁,而采用單純形法優(yōu)化后的控制器能夠恢復(fù)到穩(wěn)定運(yùn)行狀態(tài),與圖11相一致。

    圖14 引入故障情況下整流側(cè)有功功率Flg. 14 WF-VSC actlve power after fault ls lntroduced

    如圖15所示,系統(tǒng)受到故障擾動后,采用初始值參數(shù)的控制器使得系統(tǒng)在故障切除后失去穩(wěn)定,發(fā)生震蕩,最大超調(diào)量達(dá)到220 MW,采用單純形優(yōu)化后的控制器在故障切除后的0.01 s內(nèi)恢復(fù)到穩(wěn)定運(yùn)行能夠狀態(tài),且系統(tǒng)幾乎無超調(diào)量,但在3.5 s后有功功率的波動又增加。采用遺傳算法優(yōu)化后的控制器需要在故障切除后的1s后恢復(fù)到穩(wěn)定運(yùn)行狀態(tài),遺傳算法的有功功率的高頻波動很少。

    圖15 引入故障情況下逆變側(cè)有功功率Flg. 15 G-VSC actlve power after fault ls lntroduced

    如圖16所示,在原始參數(shù)情況下,故障切除后,直流電壓發(fā)生振蕩,最大超調(diào)量達(dá)到215 kV。采用單純形改進(jìn)的控制器超調(diào)量很小,能夠在故障切除后的0.01 s內(nèi)恢復(fù)到穩(wěn)定狀態(tài)。采用遺傳算法優(yōu)化后的控制器需要在故障切除后0.8 s后回到穩(wěn)定運(yùn)行狀態(tài),超調(diào)量為174.25 kV。

    圖16 引入故障情況下直流側(cè)電壓Flg. 16 DC Voltage after fault ls lntroduced

    6 結(jié)論

    本文建立了風(fēng)電場的聚合模型以及柔性直流輸電模型,風(fēng)電場經(jīng)柔性直流輸電并入交流電網(wǎng)。采用單純形和遺傳算法兩種啟發(fā)式優(yōu)化算法分別對該系統(tǒng)進(jìn)行了參數(shù)優(yōu)化,通過仿真結(jié)果可以得出如下結(jié)論:

    1)在風(fēng)電并網(wǎng)中,柔性直流輸電逆變側(cè)輸出功率的波動要小于整流側(cè)風(fēng)電場輸入功率的波動,仿真結(jié)果顯示其波動方差減小了27.6%,柔性直流輸電較好地平抑了風(fēng)電場的功率波動。這說明柔性直流輸電能夠減小風(fēng)電場的輸出波動對電網(wǎng)的影響。

    2)采用啟發(fā)式優(yōu)化算法優(yōu)化后的控制器較好地減小了風(fēng)電場輸出功率波動對并網(wǎng)功率的影響,仿真結(jié)果顯示采用單純形法和遺傳算法優(yōu)化控制器參數(shù)后,并網(wǎng)有功功率方差比優(yōu)化前減小了61.5% 和17.7%。使得風(fēng)力發(fā)電并網(wǎng)功率變得更加平滑,風(fēng)電場輸出功率能夠更好地適應(yīng)電網(wǎng)所需要的并網(wǎng)準(zhǔn)則。

    3)2種優(yōu)化算法各有優(yōu)點(diǎn):在啟動階段,遺傳算法優(yōu)化后的結(jié)果相比單純形法優(yōu)化的結(jié)果有更好的響應(yīng)速度,控制器動態(tài)響應(yīng)性能增強(qiáng);啟動完成后,單純形法的優(yōu)化結(jié)果對于平抑功率的波動好于遺傳算法,控制器魯棒性增強(qiáng)。

    4)在故障響應(yīng)情況下,故障切除后,優(yōu)化后的控制器能夠很快恢復(fù)穩(wěn)定狀態(tài),而未優(yōu)化參數(shù)的系統(tǒng)會容易失去穩(wěn)定。單純形法在故障切除后能夠更快的恢復(fù),且有較小的超調(diào)量。同時,直流電壓的平衡決定著風(fēng)電場并網(wǎng)的穩(wěn)定與否。風(fēng)速波動情況下,在遇到電網(wǎng)側(cè)故障擾動時會使得風(fēng)電場并入系統(tǒng)的功率發(fā)生振蕩,通過優(yōu)化參數(shù)能夠改善這一情況。

    綜上所述,采用啟發(fā)式優(yōu)化算法能夠改善風(fēng)電場經(jīng)柔性直流輸電系統(tǒng)的動態(tài)響應(yīng)能力。在風(fēng)速大幅變化的情況中,優(yōu)化后的柔性直流輸電系統(tǒng)更能平抑風(fēng)電場輸出功率波動。2種遺傳算法各有優(yōu)劣,對于采用何種優(yōu)化需要根據(jù)工程實(shí)際進(jìn)行選擇。

    參考文獻(xiàn)

    [1]趙成勇.柔性直流輸電建模和仿真技術(shù)[M].北京:中國電力出版社,2014.

    [2]廖勇,王國棟.雙饋風(fēng)電場的柔性高壓直流輸電系統(tǒng)控制[J].中國電機(jī)工程學(xué)報,2012,32(28):7-15. LIAO Yong,WANG Guodong. VSC-HVDC system contro1 for grid-connection of DIFG wind farms[J]. Proceedings of the CSEE. 2012,32(28):7-15(in Chinese).

    [3]范心明,管霖,夏成軍,等.風(fēng)電場交直流混合輸電并網(wǎng)中VSC-HVDC的控制[J].中國電機(jī)工程學(xué)報,2014,34(28):4781-4790. FAN Xinming,GUAN Lin,XIA Chengjun,et a1. Contro1 of VSC-HVDC in AC/DC hybrid transmission with wind farms integrated[J]. Proceedings of the CSEE,2014,34 (28):4781-4790(in Chinese).

    [4]郭春義,趙成勇,李廣凱,等.基于SimP1ex算法的VSCHVDC控制參數(shù)優(yōu)化[J].電力自動化設(shè)備,2010,30(9):13-17. GUO Chunyi,ZHAO Chengyong,LI Guangkai,et a1. Contro1 Parameter oPtimization based on SimP1ex a1gorithm for VSC-HVDC[J]. E1ectric Power Au-tomation EquiPment,2010,30(9):13-17(in Chinese).

    [5]陳蔓,陸繼明,毛承雄,等.基于遺傳算法的優(yōu)化控制在VSC-HVDC中的應(yīng)用[J].電力系統(tǒng)及其自動化學(xué)報,2006,18(4):19-23. CHEN Man,LU Jiming,MAO Chengxiong,et a1. APP1ication of genetic a1gorithm based oPtima1 contro1 in VSC-HVDC[J]. Proceedings of the CSU-EPSA. 2006,18 (4):19-23(in Chinese).

    [6]于永生,馮延暉,江紅鑫,等.海上風(fēng)電經(jīng)VSC-MTDC并網(wǎng)研究[J].高壓電器,2015,51(10):26-33. YU Yongsheng,F(xiàn)ENG Yanhui,JIANG Hongxin,et a1. Research on offshore wind farm integration via VSCMTDC[J]. High Vo1tage APParatus,2015,51(10):26-33 (in Chinese).

    [7]陳樹勇,徐林巖,孫栩,等.基于多端柔性直流輸電的風(fēng)電并網(wǎng)控制研究[J].中國電機(jī)工程學(xué)報,2014,34:32-38. CHEN Shuyong,XU Linyan,SUN Xu,et a1. The contro1 of wind Power integration based on mu1ti -termina1 high vo1tage DC transmission with vo1tage source converter[J]. Proceedings of the CSEE. 2014,34:32-38(in Chinese).

    [8]湯廣福.基于電壓源換流器的高壓直流輸電技術(shù)[M].北京:中國電力出版社,2010:132-144.

    [9]張興,張崇巍. PWM整流器及其控制[M].北京:機(jī)械工業(yè)出版社,2012:88-106.

    [10]李環(huán)平,楊金明.基于PSCAD_EMTDC的大型并網(wǎng)風(fēng)電場建模與仿真[J].電力系統(tǒng)保護(hù)與控制,2009,37(21):62-66. LI HuanPing,YANG Jinming. Mode1 and simu1ation of 1arge grid connected wind farm based on PSCAD_EMTDC[J]. Power System Protection and Contro1,2009,37(21):62-66.

    [11] CONROY J,WATSON R. Aggregate mode11ing of wind farms containing fu11-converter wind turbine generators with Permanent magnet synchronous machines -transient stabi1ity studies[J]. IET Renewab1e Power Generation,2009,3(1):39-52.

    [12]夏玥,李征,蔡旭,等.基于直驅(qū)式永磁同步發(fā)電機(jī)組的風(fēng)電場動態(tài)建模[J].電網(wǎng)技術(shù),2014,38(6):1439-1445. XIA Yue,LI Zheng,CAI Xu,et a1. Dynamic mode1ing of wind farm comPosed of direct—driven Permanent magnet synchronous generators[J]. Power System Techno1ogy,2014,38(6):1439-1445(in Chinese).

    [13]楊勇波,江波,楊華云,等.風(fēng)電場動態(tài)等值及其在暫態(tài)穩(wěn)定分析中的應(yīng)用[J].中國電力,2011,44(7):66-70. YANG Yongbo,JIANG Bo,YANC Huayun,et a1. Dynamic equiva1ence of wind farms and its aPP1ication to Power systemtransient stabi1ity ana1ysis[J]. E1ectric Power,2011,44(7):66-70(in Chinese).

    [14]肖宏峰.基于單純形多向搜索的大規(guī)模進(jìn)化優(yōu)化算法[D].長沙:中南大學(xué),2009.

    [15]項國波. ITAE最佳控制[M].北京:機(jī)械工業(yè)出版社,1986.

    [16]張勇軍.電力系統(tǒng)無功優(yōu)化的災(zāi)變遺傳算法及MAS模型研究[D].廣州:華南理工大學(xué),2004:17-30.

    [17]于永生,馮延暉,江紅鑫,等.一種風(fēng)電場經(jīng)柔性直流輸電并網(wǎng)控制方法改進(jìn)探討[J].電網(wǎng)與清潔能源,2015,31(7):81-87. YU Yongsheng,F(xiàn)ENG Yanhui,JIANG Hongxin,et a1. An imProvement on 1arge-sca1e wind Power grid integration through VSC-HVDC[J]. Power System and C1ean Energy,2015,31(7):81-87(in Chinese).

    [18] ZHAO Xiaodong,LI Kang. Contro1 of VSC-HVDC for wind farm integration based on adaPtive backstePPing method[C]. IEEE Internationa1 WorkshoP on Inte11igent Energy System(IWIES),2013:64-69.

    [19] WANG Liying,ERTUGRUL N. Se1ection of PI comPensator Parameters for VSC-HVDC system using decouP1ed contro1 strategy[C]. 20th Austra1asian Universities Power Engineering Conference(AUPEC),2010:1-7.

    [20] YANG Jin,O'REILLY J,F(xiàn)LETCHER J E. An overview of DC cab1e mode1ing for fau1t ana1ysis of VSC -HVDC transmission systems[C]. 20th Austra1asian Universities Power Engineering Conference(AUPEC),2010:1-7.

    [21]陜?nèi)A平,肖登明,薛愛東.大型風(fēng)電場的風(fēng)資源評估[J].華東電力,2006,34(2):16-17. SHAN HuaPing,XIAO Dengming,XUE Aidong. Wind resource assessment for 1arge wind Power fie1ds[J]. East China E1ectric Power,2006,34(2):16-17(in Chinese).

    于永生(1989—),男,碩士研究生,主要從事柔性直流輸電,風(fēng)電并網(wǎng)等方面的研究工作;

    馮延暉(1977—),男,博士,副教授,主要研究方向為可再生能源發(fā)電,電力電子與電力傳動,新能源并網(wǎng)控制。

    (編輯徐花榮)

    Project SuPPorted by Jiangsu ToP Six Ta1ent Summit Fund(ZBZZ-045);Natura1 Science Fund of Jiangsu Province(BK 20131350);the Fundamenta1 Research Funds for the Centra1 Universities(309150 11324);Returned Overseas Students Preferred Funding.

    0ptlmlzatlon of Wlnd Farm Grld Integratlon vla VSC-HVDC Conslderlng Randomness of Wlnd Energy

    YU Yongsheng,F(xiàn)ENG Yanhui,JIANG Hongxin,QIU Yingning
    (Schoo1 of Energy and Power Engineering,Nanjing University of Science and Techno1ogy,Nanjing 210094,Jiangsu,China)

    ABSTRACT:As energy crisis has draws more and more attention a11 around the wor1d,wind energy is getting more and more imPortant as a kind of renewab1e and c1ean energy. The outPut of the wind energy,however,is subjected to changes of wind sPeed,thus,the wind farm outPut often f1uctuates. This PaPer bui1ds the aggregate wind farm mode1 and the VSC -HVDC mode1. In order to stabi1ize the f1uctuation of the wind farm outPut,the PaPer adoPts the SimP1ex Method and Genetic A1gorithm resPective1y to oPtimize the contro11er of VSC-HVDC. The resu1t shows that both of these two different ways he1P to constrain the f1uctuation of the outPut of the wind farm and have strong engineering aPP1ication va1ue.

    作者簡介:

    收稿日期:2015-06-17。

    基金項目:江蘇省六大人才高峰項目(ZBZZ-045);江蘇省自然科學(xué)基金面上項目(BK20131350);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項資金資助(30915011324);留學(xué)回國人員擇優(yōu)資助項目。

    文章編號:1674-3814(2016)01-0075-10

    中圖分類號:TM743

    文獻(xiàn)標(biāo)志碼:A

    猜你喜歡
    單純形法風(fēng)電場遺傳算法
    基于單純形法的TLE軌道確定
    基于單純形法的簡單問題的研究與應(yīng)用
    青年生活(2019年35期)2019-09-10 00:13:32
    基于PSS/E的風(fēng)電場建模與動態(tài)分析
    電子制作(2018年17期)2018-09-28 01:56:44
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
    線性規(guī)劃最優(yōu)解研究
    基于遺傳算法和LS-SVM的財務(wù)危機(jī)預(yù)測
    基于改進(jìn)單純形法的冗余證券的判別
    含風(fēng)電場電力系統(tǒng)的潮流計算
    基于改進(jìn)的遺傳算法的模糊聚類算法
    国产精品久久久人人做人人爽| 亚洲精品久久国产高清桃花| 亚洲成人久久性| 国产亚洲精品一区二区www| 99热这里只有精品一区| 精品福利观看| 中文字幕久久专区| 一本一本综合久久| 在线观看66精品国产| 人妻丰满熟妇av一区二区三区| 免费看十八禁软件| svipshipincom国产片| 国产av在哪里看| 丰满的人妻完整版| 久久久精品欧美日韩精品| 天美传媒精品一区二区| 午夜福利免费观看在线| 村上凉子中文字幕在线| 欧美区成人在线视频| 人人妻人人澡欧美一区二区| 国产精品亚洲美女久久久| 真人做人爱边吃奶动态| 欧美一级a爱片免费观看看| 色综合婷婷激情| 免费av不卡在线播放| 久久久久亚洲av毛片大全| 免费电影在线观看免费观看| 亚洲av一区综合| 真人一进一出gif抽搐免费| 成人鲁丝片一二三区免费| 最新美女视频免费是黄的| 午夜免费男女啪啪视频观看 | 亚洲专区国产一区二区| 男人舔奶头视频| 一进一出好大好爽视频| 两个人视频免费观看高清| 亚洲成人久久爱视频| 校园春色视频在线观看| 国产精品野战在线观看| 国产老妇女一区| 久久欧美精品欧美久久欧美| 免费av毛片视频| 国产精品精品国产色婷婷| 禁无遮挡网站| 欧美另类亚洲清纯唯美| 波多野结衣巨乳人妻| 熟女电影av网| 色在线成人网| www.999成人在线观看| 91久久精品国产一区二区成人 | 久久久久免费精品人妻一区二区| 欧美三级亚洲精品| 18禁黄网站禁片免费观看直播| 在线免费观看的www视频| 日本a在线网址| 女警被强在线播放| 最新美女视频免费是黄的| 久久精品国产综合久久久| 国语自产精品视频在线第100页| 一区二区三区高清视频在线| av天堂在线播放| 九九久久精品国产亚洲av麻豆| 免费高清视频大片| 欧美成人性av电影在线观看| 国产亚洲av嫩草精品影院| 国产日本99.免费观看| 国产精品久久电影中文字幕| 国产 一区 欧美 日韩| 国产伦在线观看视频一区| 国产视频一区二区在线看| 嫩草影视91久久| 亚洲美女视频黄频| 天堂√8在线中文| 免费大片18禁| 青草久久国产| 日本撒尿小便嘘嘘汇集6| 一卡2卡三卡四卡精品乱码亚洲| 国产色婷婷99| 色老头精品视频在线观看| 久久精品91蜜桃| 波多野结衣高清无吗| 五月伊人婷婷丁香| 久久精品国产亚洲av香蕉五月| 99久久精品热视频| 亚洲精品美女久久久久99蜜臀| 国产亚洲欧美在线一区二区| 国内少妇人妻偷人精品xxx网站| 国产美女午夜福利| 久久午夜亚洲精品久久| 亚洲国产日韩欧美精品在线观看 | 国产探花极品一区二区| 天天一区二区日本电影三级| 国产单亲对白刺激| 国产三级中文精品| 欧美三级亚洲精品| h日本视频在线播放| 成年免费大片在线观看| 国产成人啪精品午夜网站| 18禁黄网站禁片午夜丰满| 亚洲av免费高清在线观看| 美女高潮的动态| 成人18禁在线播放| 国产中年淑女户外野战色| 99在线视频只有这里精品首页| 最近最新免费中文字幕在线| 高清毛片免费观看视频网站| 久久久久精品国产欧美久久久| 亚洲中文日韩欧美视频| 观看免费一级毛片| 国产伦在线观看视频一区| 亚洲人成网站在线播放欧美日韩| 人人妻,人人澡人人爽秒播| 日本 av在线| 久久久久久久久大av| 亚洲精品粉嫩美女一区| 97超视频在线观看视频| 噜噜噜噜噜久久久久久91| 性色av乱码一区二区三区2| 欧美在线一区亚洲| 成人无遮挡网站| 久久久久性生活片| 久久精品国产99精品国产亚洲性色| 欧美国产日韩亚洲一区| 亚洲久久久久久中文字幕| 亚洲精品乱码久久久v下载方式 | 久久精品91蜜桃| 日韩欧美 国产精品| 亚洲电影在线观看av| 身体一侧抽搐| av中文乱码字幕在线| 九九久久精品国产亚洲av麻豆| 亚洲av免费在线观看| 中文字幕熟女人妻在线| 欧美激情在线99| 久久久成人免费电影| 日韩精品青青久久久久久| 亚洲欧美日韩无卡精品| 日韩成人在线观看一区二区三区| 啦啦啦观看免费观看视频高清| 99热这里只有是精品50| 欧美国产日韩亚洲一区| 久久久国产成人精品二区| 亚洲欧美精品综合久久99| 亚洲国产中文字幕在线视频| 国产激情欧美一区二区| 欧美激情久久久久久爽电影| 国产黄片美女视频| 老熟妇仑乱视频hdxx| 亚洲国产精品sss在线观看| 免费观看人在逋| 国产成人欧美在线观看| 免费无遮挡裸体视频| 午夜福利欧美成人| 亚洲五月天丁香| 久久久久精品国产欧美久久久| 免费av不卡在线播放| 男女床上黄色一级片免费看| 一级作爱视频免费观看| 搡老熟女国产l中国老女人| 亚洲av成人不卡在线观看播放网| 日韩av在线大香蕉| 国产成人啪精品午夜网站| 成人特级黄色片久久久久久久| 欧美黑人欧美精品刺激| 99精品久久久久人妻精品| 一级毛片高清免费大全| 一级a爱片免费观看的视频| 日日干狠狠操夜夜爽| 法律面前人人平等表现在哪些方面| 日韩大尺度精品在线看网址| 亚洲va日本ⅴa欧美va伊人久久| 91字幕亚洲| 久久久久久国产a免费观看| 非洲黑人性xxxx精品又粗又长| 精品久久久久久久人妻蜜臀av| 女人高潮潮喷娇喘18禁视频| 中亚洲国语对白在线视频| 内地一区二区视频在线| 久久精品91蜜桃| 有码 亚洲区| 日本撒尿小便嘘嘘汇集6| 天天一区二区日本电影三级| 色吧在线观看| e午夜精品久久久久久久| 亚洲va日本ⅴa欧美va伊人久久| av国产免费在线观看| 夜夜爽天天搞| 成人一区二区视频在线观看| 精品久久久久久,| 一个人观看的视频www高清免费观看| 高清在线国产一区| 成人av一区二区三区在线看| 观看美女的网站| 99久久精品国产亚洲精品| 色尼玛亚洲综合影院| 在线观看66精品国产| 午夜影院日韩av| 成人国产一区最新在线观看| 欧美成人一区二区免费高清观看| 国产一区二区在线观看日韩 | 国产真人三级小视频在线观看| 亚洲熟妇中文字幕五十中出| 欧美性感艳星| 精品欧美国产一区二区三| 精品国内亚洲2022精品成人| 日韩欧美精品v在线| 国产黄色小视频在线观看| 精品午夜福利视频在线观看一区| 成年女人看的毛片在线观看| 18禁在线播放成人免费| 99热这里只有精品一区| 国产午夜福利久久久久久| 免费高清视频大片| 午夜激情福利司机影院| 日韩国内少妇激情av| 久久精品国产自在天天线| 色精品久久人妻99蜜桃| 丁香六月欧美| 国产精品香港三级国产av潘金莲| 麻豆国产av国片精品| 91九色精品人成在线观看| 一a级毛片在线观看| 午夜免费激情av| 啦啦啦韩国在线观看视频| 大型黄色视频在线免费观看| 日日干狠狠操夜夜爽| 亚洲国产精品久久男人天堂| av天堂在线播放| 两个人视频免费观看高清| 国产在线精品亚洲第一网站| 一进一出抽搐gif免费好疼| 欧美一区二区亚洲| 国产久久久一区二区三区| 国产一区二区三区视频了| 淫秽高清视频在线观看| 内地一区二区视频在线| 国产精品综合久久久久久久免费| 国产高清有码在线观看视频| 免费在线观看日本一区| 婷婷精品国产亚洲av在线| 岛国在线免费视频观看| 亚洲在线观看片| 此物有八面人人有两片| 母亲3免费完整高清在线观看| 国产av一区在线观看免费| 日本免费一区二区三区高清不卡| 69av精品久久久久久| 两个人看的免费小视频| 国产午夜精品久久久久久一区二区三区 | 搡老岳熟女国产| 18禁国产床啪视频网站| 国产精品女同一区二区软件 | 俄罗斯特黄特色一大片| 国产精品 欧美亚洲| 每晚都被弄得嗷嗷叫到高潮| 久久久久久久精品吃奶| 美女被艹到高潮喷水动态| 女人高潮潮喷娇喘18禁视频| 青草久久国产| 欧美乱色亚洲激情| 老汉色av国产亚洲站长工具| 亚洲欧美精品综合久久99| 久久精品91蜜桃| 午夜福利免费观看在线| 精品人妻一区二区三区麻豆 | 久久精品国产清高在天天线| 国产视频内射| 久久精品国产亚洲av涩爱 | 国产精品久久电影中文字幕| 婷婷亚洲欧美| 精品国产三级普通话版| 99久久无色码亚洲精品果冻| 欧洲精品卡2卡3卡4卡5卡区| 黄色日韩在线| 亚洲国产精品sss在线观看| 深夜精品福利| 日韩亚洲欧美综合| 啦啦啦观看免费观看视频高清| 国产精品1区2区在线观看.| 男女下面进入的视频免费午夜| 51午夜福利影视在线观看| 99久久精品国产亚洲精品| 国产伦精品一区二区三区四那| 日本成人三级电影网站| av欧美777| 最近在线观看免费完整版| 国产精品国产高清国产av| 国产av不卡久久| 老司机午夜十八禁免费视频| 国产视频一区二区在线看| 俄罗斯特黄特色一大片| 精品久久久久久久久久免费视频| 欧美在线黄色| 黄色片一级片一级黄色片| 欧美激情久久久久久爽电影| 成年女人永久免费观看视频| 国产又黄又爽又无遮挡在线| 日本熟妇午夜| 中文字幕人成人乱码亚洲影| 黄色片一级片一级黄色片| 9191精品国产免费久久| 欧美日韩瑟瑟在线播放| 日本黄大片高清| 欧美日韩国产亚洲二区| 色哟哟哟哟哟哟| 亚洲乱码一区二区免费版| 丁香欧美五月| av中文乱码字幕在线| 国产成人福利小说| 麻豆国产av国片精品| 国产又黄又爽又无遮挡在线| 亚洲专区国产一区二区| 亚洲中文字幕日韩| 亚洲18禁久久av| 国产一区在线观看成人免费| 亚洲国产精品成人综合色| 桃红色精品国产亚洲av| 国产成人av教育| netflix在线观看网站| 国产伦在线观看视频一区| 国产三级中文精品| 国产高清视频在线观看网站| 亚洲精品久久国产高清桃花| 成熟少妇高潮喷水视频| av福利片在线观看| 搡女人真爽免费视频火全软件 | 极品教师在线免费播放| 黄色片一级片一级黄色片| 91在线观看av| 国产乱人视频| 少妇裸体淫交视频免费看高清| 成熟少妇高潮喷水视频| 亚洲午夜理论影院| 亚洲av电影不卡..在线观看| 亚洲欧美日韩东京热| 欧美av亚洲av综合av国产av| 欧美性感艳星| 久久久久久久精品吃奶| 亚洲av电影不卡..在线观看| 天堂av国产一区二区熟女人妻| 日本五十路高清| 亚洲av第一区精品v没综合| 亚洲avbb在线观看| 国产乱人伦免费视频| 欧美又色又爽又黄视频| 深爱激情五月婷婷| 女警被强在线播放| 日韩av在线大香蕉| 无人区码免费观看不卡| 成人精品一区二区免费| 精品人妻1区二区| 亚洲av第一区精品v没综合| 日本在线视频免费播放| 国产一区二区三区视频了| 国产精品女同一区二区软件 | 国产熟女xx| 精品一区二区三区人妻视频| 香蕉av资源在线| 中文字幕精品亚洲无线码一区| 最新美女视频免费是黄的| 国产一区二区三区视频了| 日本精品一区二区三区蜜桃| 麻豆国产av国片精品| 国产探花在线观看一区二区| 国产亚洲精品久久久久久毛片| 亚洲自拍偷在线| 国产成人福利小说| 欧美色视频一区免费| 高清在线国产一区| 国产精品野战在线观看| 欧美最新免费一区二区三区 | 俺也久久电影网| 久久伊人香网站| 国产精品乱码一区二三区的特点| 91在线观看av| 少妇高潮的动态图| 日本五十路高清| 国产黄a三级三级三级人| 欧美乱妇无乱码| 级片在线观看| 国语自产精品视频在线第100页| 久久久久性生活片| 亚洲av成人av| 午夜精品在线福利| 美女被艹到高潮喷水动态| www.熟女人妻精品国产| 国产激情欧美一区二区| 色av中文字幕| 免费高清视频大片| 日韩亚洲欧美综合| 中文字幕久久专区| 久久精品91蜜桃| 日韩 欧美 亚洲 中文字幕| 国产av不卡久久| 男女视频在线观看网站免费| 老司机福利观看| 在线播放无遮挡| 成人av一区二区三区在线看| 精品乱码久久久久久99久播| 99在线视频只有这里精品首页| 老司机福利观看| 少妇高潮的动态图| 深夜精品福利| 欧美极品一区二区三区四区| 成人欧美大片| 国产亚洲精品一区二区www| 免费看光身美女| 亚洲美女黄片视频| 精华霜和精华液先用哪个| av视频在线观看入口| 免费在线观看日本一区| 动漫黄色视频在线观看| 一进一出好大好爽视频| 欧美黑人巨大hd| 国产 一区 欧美 日韩| 婷婷精品国产亚洲av| 神马国产精品三级电影在线观看| 国产精品女同一区二区软件 | 国产久久久一区二区三区| 久久草成人影院| 国内精品一区二区在线观看| 亚洲成人精品中文字幕电影| 精品久久久久久久毛片微露脸| 婷婷精品国产亚洲av在线| 亚洲乱码一区二区免费版| 制服丝袜大香蕉在线| 成人18禁在线播放| 天天一区二区日本电影三级| 身体一侧抽搐| 欧美黑人欧美精品刺激| 一a级毛片在线观看| 日本精品一区二区三区蜜桃| a在线观看视频网站| 国产v大片淫在线免费观看| 免费看光身美女| 日韩中文字幕欧美一区二区| 亚洲自拍偷在线| a级毛片a级免费在线| 99riav亚洲国产免费| 日本一本二区三区精品| 成年女人看的毛片在线观看| 国产精品日韩av在线免费观看| 成年免费大片在线观看| 精品久久久久久久人妻蜜臀av| 国产一区二区三区视频了| 精品国产超薄肉色丝袜足j| 国产精品亚洲美女久久久| 小说图片视频综合网站| 亚洲激情在线av| 亚洲成av人片免费观看| 小蜜桃在线观看免费完整版高清| 在线观看一区二区三区| 国产精品亚洲av一区麻豆| 国产精品1区2区在线观看.| 美女高潮喷水抽搐中文字幕| 亚洲欧美精品综合久久99| 国产精品电影一区二区三区| 最好的美女福利视频网| 国产亚洲精品久久久久久毛片| 精品久久久久久久毛片微露脸| xxxwww97欧美| 国产午夜精品论理片| 国产精品香港三级国产av潘金莲| 黄色视频,在线免费观看| 一级黄片播放器| 在线十欧美十亚洲十日本专区| 国产一区二区在线观看日韩 | 亚洲精品日韩av片在线观看 | 国产真实乱freesex| 成人欧美大片| 国产一级毛片七仙女欲春2| 18禁在线播放成人免费| 麻豆国产av国片精品| 午夜久久久久精精品| 亚洲美女黄片视频| 国产精品久久视频播放| 久久久久久久亚洲中文字幕 | 欧美高清成人免费视频www| 亚洲精品在线观看二区| 日本免费a在线| 老司机午夜十八禁免费视频| 亚洲成人久久性| 18禁在线播放成人免费| 熟女少妇亚洲综合色aaa.| 国产精品久久久久久久电影 | 日本与韩国留学比较| 在线免费观看的www视频| 99久久成人亚洲精品观看| 国内少妇人妻偷人精品xxx网站| 免费大片18禁| 少妇的丰满在线观看| 国产午夜精品久久久久久一区二区三区 | 欧美av亚洲av综合av国产av| 日本 欧美在线| 国产伦精品一区二区三区四那| 国产欧美日韩一区二区三| 三级毛片av免费| 啦啦啦韩国在线观看视频| 精品久久久久久,| 免费看十八禁软件| 亚洲成av人片在线播放无| 久久久成人免费电影| 色视频www国产| 黄色日韩在线| 婷婷精品国产亚洲av| 12—13女人毛片做爰片一| 夜夜看夜夜爽夜夜摸| 少妇丰满av| 91麻豆精品激情在线观看国产| 国产欧美日韩一区二区精品| 国产av麻豆久久久久久久| 色在线成人网| 一本一本综合久久| 中文亚洲av片在线观看爽| 欧美日韩乱码在线| 精品99又大又爽又粗少妇毛片 | 久久久久性生活片| 久久久久免费精品人妻一区二区| ponron亚洲| 制服丝袜大香蕉在线| 免费大片18禁| 成年女人永久免费观看视频| 国内精品久久久久久久电影| 免费av观看视频| 搡女人真爽免费视频火全软件 | tocl精华| 在线观看66精品国产| 综合色av麻豆| 内地一区二区视频在线| 国产一级毛片七仙女欲春2| 国产精品精品国产色婷婷| 在线播放无遮挡| 日韩欧美国产一区二区入口| 国产成人福利小说| 免费搜索国产男女视频| 亚洲av成人av| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久末码| 久久伊人香网站| 中文字幕人成人乱码亚洲影| 男人的好看免费观看在线视频| 天堂影院成人在线观看| 欧美不卡视频在线免费观看| 色综合亚洲欧美另类图片| 亚洲人成伊人成综合网2020| 人人妻人人看人人澡| 亚洲国产欧洲综合997久久,| 久久草成人影院| 在线天堂最新版资源| 亚洲熟妇熟女久久| 亚洲欧美日韩高清在线视频| 亚洲色图av天堂| 麻豆一二三区av精品| 欧美乱码精品一区二区三区| 亚洲国产精品999在线| 国产一区二区在线av高清观看| 国内久久婷婷六月综合欲色啪| 美女大奶头视频| 校园春色视频在线观看| 一区二区三区高清视频在线| 国产一区二区在线av高清观看| 99在线人妻在线中文字幕| 国产亚洲欧美98| 精品国产亚洲在线| 99热只有精品国产| 乱人视频在线观看| 非洲黑人性xxxx精品又粗又长| 97碰自拍视频| 中国美女看黄片| 国产亚洲精品一区二区www| 99热只有精品国产| 亚洲在线自拍视频| 成人三级黄色视频| 国产精品 国内视频| h日本视频在线播放| 中文字幕熟女人妻在线| 亚洲不卡免费看| 我的老师免费观看完整版| 国产伦人伦偷精品视频| 午夜日韩欧美国产| 中文字幕熟女人妻在线| 国产 一区 欧美 日韩| 国产精品香港三级国产av潘金莲| 国产精华一区二区三区| 欧美乱色亚洲激情| 国产精品日韩av在线免费观看| 在线观看免费午夜福利视频| 亚洲专区中文字幕在线| a在线观看视频网站| 黄色成人免费大全| 亚洲五月婷婷丁香| 亚洲av二区三区四区| 久久6这里有精品| 精品人妻一区二区三区麻豆 | 国产精品久久久久久亚洲av鲁大| 亚洲狠狠婷婷综合久久图片| 人人妻人人看人人澡| 久久99热这里只有精品18| 精品欧美国产一区二区三| 岛国在线观看网站| 欧美在线一区亚洲| 久久精品国产综合久久久| 久久精品影院6| 亚洲中文日韩欧美视频| 黄色女人牲交| 97超级碰碰碰精品色视频在线观看| 12—13女人毛片做爰片一| 国内精品久久久久久久电影| 日本黄色视频三级网站网址| 乱人视频在线观看| 在线观看舔阴道视频| 亚洲精品成人久久久久久| 麻豆国产av国片精品| 又黄又粗又硬又大视频| 亚洲国产欧洲综合997久久,|