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

    含風(fēng)光出力隨機(jī)性的獨(dú)立微電網(wǎng)二次頻率控制

    2024-08-21 00:00:00鐘誠姜志富張翔宇陳繼開李揚(yáng)
    太陽能學(xué)報(bào) 2024年1期
    關(guān)鍵詞:模型預(yù)測控制

    摘要:針對風(fēng)、光發(fā)電隨機(jī)性的問題,提出一種考慮新能源減載參與二次調(diào)頻的模型預(yù)測控制方法。建立包含隨機(jī)功率擾動的擴(kuò)展?fàn)顟B(tài)矩陣,采用卡爾曼濾波估算隨機(jī)未知擾動;依據(jù)風(fēng)、光最大可用功率,建立實(shí)時(shí)變約束,避免機(jī)組功率越限;設(shè)置合理的權(quán)重系數(shù),優(yōu)先風(fēng)、光發(fā)電出力參與二次調(diào)頻;通過求解變約束二次規(guī)劃問題,獲得各個(gè)機(jī)組的優(yōu)化調(diào)頻功率。最后,建立含多個(gè)光伏、風(fēng)電的微電網(wǎng)模型,在不同場景下與常規(guī)二次調(diào)頻方法進(jìn)行對比仿真。仿真結(jié)果表明,所提方法能提高系統(tǒng)頻率恢復(fù)速度,減小系統(tǒng)頻率波動,尤其在風(fēng)、光發(fā)電劇烈波動場景下。

    關(guān)鍵詞:二次調(diào)頻;模型預(yù)測控制;獨(dú)立微電網(wǎng);減載控制;隨機(jī)擾動觀測

    中圖分類號:TM727 文獻(xiàn)標(biāo)志碼:A

    0引言

    隨著國家“雙碳”戰(zhàn)略目標(biāo)制定,新能源發(fā)電得到進(jìn)一步重視,將成為未來主要發(fā)電形式1。微電網(wǎng)能有效整合各種分布式發(fā)電和儲能,實(shí)現(xiàn)負(fù)荷就近供電,為提高新能源接入水平提供一種途徑,是未來電網(wǎng)的重要組成部分2]。

    微電網(wǎng)可工作在并網(wǎng)和孤島兩種模式[3]。孤島模式下,微電網(wǎng)頻率穩(wěn)定需依靠其包含的分布式發(fā)電機(jī)組調(diào)節(jié)[4]。但這些電源慣性小,出力隨機(jī)性大(例如風(fēng)力、光伏發(fā)電)。因此,微電網(wǎng)需求更加先進(jìn)靈活控制方法,來提高頻率控制能力5。類似于常規(guī)大電網(wǎng),微電網(wǎng)頻率調(diào)節(jié)可分成一次調(diào)頻、二次調(diào)頻,甚至三次調(diào)頻L?-7]。一次調(diào)頻為本地有差調(diào)節(jié),無法實(shí)現(xiàn)頻率恢復(fù)。引入二次調(diào)頻控制,可進(jìn)一步調(diào)整分布式發(fā)電出力,提高頻率品質(zhì)。

    常規(guī)二次調(diào)頻通過采用PI控制器幫助系統(tǒng)頻率恢復(fù)。但微電網(wǎng)中電源種類差異較大,集中PI控制不能靈活整合各種調(diào)頻資源。文獻(xiàn)[8]采用PI控制器來實(shí)現(xiàn)二次調(diào)頻控制,但其參數(shù)整定受限于非機(jī)理模型的準(zhǔn)確性;文獻(xiàn)[9]提出一種基于自適應(yīng)動態(tài)規(guī)劃(adaptive dynamic programming,ADP)的控制策略,該方法考慮了分布式電源出力的不確定性以及負(fù)荷的隨機(jī)性,調(diào)整柴油機(jī)和儲能系統(tǒng)的出力來抑制頻率波動;文獻(xiàn)[10]提出能將一、二次調(diào)頻進(jìn)行切換的虛擬同步發(fā)電機(jī)(virtual synchronous generator,VSG)控制,但未考慮變流器響應(yīng)速度不同對頻率調(diào)整的影響;文獻(xiàn)[11]提出一種多微電網(wǎng)系統(tǒng)魯棒模型預(yù)測控制方法,利用系統(tǒng)之間的耦合來減小微電網(wǎng)系統(tǒng)的頻率波動,從而提高系統(tǒng)的穩(wěn)定性;文獻(xiàn)[12]提出一種自適應(yīng)模型預(yù)測控制方法,用于確保有新能源發(fā)電裝置的兩區(qū)域互聯(lián)的頻率控制。

    近年來,模型預(yù)測控制(model predictive control,MPC)得到迅速發(fā)展,其滾動優(yōu)化特點(diǎn)可較好處理風(fēng)光實(shí)時(shí)變化的出力所引起的頻率波動。文獻(xiàn)[13]提出一種模糊控制結(jié)合MPC來控制VSG的控制方式,模糊MPC控制器通過修正虛擬慣量和阻尼系數(shù),從而提高孤島微電網(wǎng)的頻率穩(wěn)定性,但模糊控制主要依賴于經(jīng)驗(yàn)和湊試,且控制規(guī)則一旦確定,不能實(shí)時(shí)調(diào)整,不太適合處理隨機(jī)性較強(qiáng)的環(huán)境。文獻(xiàn)[14]提出一種可優(yōu)化微電網(wǎng)中儲能電池之間的功率潮流模型預(yù)測控制策略,同時(shí)使用求解器進(jìn)行快速處理儲能電池的非線性變化,可提高實(shí)時(shí)MPC的運(yùn)算速度,保證控制器更加快速的調(diào)節(jié)頻率,但僅考慮了儲能系統(tǒng),若系統(tǒng)新能源滲透率較高、波動較大時(shí),只有儲能參與調(diào)頻可能會出現(xiàn)頻率越限問題。文獻(xiàn)[15]針對于VSG二次調(diào)頻能力不足的問題提出了MPC的控制方式,來減少頻率恢復(fù)時(shí)間,但也僅僅是結(jié)合VSG和MPC的控制方式,針對于微電網(wǎng)的新能源隨機(jī)性引起的頻率問題未深入的討論,上述研究多集中于傳統(tǒng)儲能裝置,未涉及到新能源發(fā)電機(jī)組參與調(diào)頻情況。近年來,風(fēng)、光發(fā)電有功備用參與系統(tǒng)頻率得到關(guān)注[16-18]。風(fēng)、光發(fā)電采用減載模式預(yù)留部分有功功率,參與系統(tǒng)頻率調(diào)節(jié)改善頻率響應(yīng)。但風(fēng)、光發(fā)電具有隨機(jī)性,使其參與二次調(diào)頻增加難度。

    綜上,提出一種考慮新能源減載參與二次調(diào)頻控制方法。該方法采用模型預(yù)測控制(MPC)構(gòu)架,實(shí)時(shí)優(yōu)化各個(gè)分布式單元出力。為緩解風(fēng)、光發(fā)電隨機(jī)的影響,建立包含隨機(jī)擾動的擴(kuò)展?fàn)顟B(tài)矩陣,采用卡爾曼濾波估算隨機(jī)擾動;依據(jù)實(shí)時(shí)估算風(fēng)、光最大功率估算值,建立實(shí)時(shí)變約束。另外,通過合理權(quán)重設(shè)計(jì),優(yōu)先風(fēng)、光發(fā)電出力。最后,建立微電網(wǎng)仿真模型,在不同場景下與傳統(tǒng)控制策略進(jìn)行對比驗(yàn)證。

    1微電網(wǎng)模型

    1.1微電網(wǎng)結(jié)構(gòu)

    獨(dú)立微電網(wǎng)結(jié)構(gòu)如圖1所示。主要包含兩個(gè)分布式光伏發(fā)電機(jī)組(PV1,PV2)、兩個(gè)分布式風(fēng)力發(fā)電機(jī)組(WT1,WT2)、儲能單元(battery energy storage system,BESS)和柴油機(jī)組(diesel unit,DU),及其對應(yīng)的負(fù)荷。各發(fā)電單元的容量配置如表1所示。

    微電網(wǎng)獨(dú)立運(yùn)行時(shí),需依靠自身發(fā)電單元調(diào)節(jié),維持電壓和頻率穩(wěn)定。微電網(wǎng)頻率控制包含一次調(diào)頻和二次調(diào)頻兩層。一次調(diào)頻為本地控制,發(fā)電單元通過下垂控制,依據(jù)頻率波動改變機(jī)組有功出力,減少功率擾動下的頻率偏差。但一次調(diào)頻為有差調(diào)節(jié),為維持微電網(wǎng)頻率恒定,需增加二次調(diào)頻控制。二次調(diào)頻控制協(xié)調(diào)多個(gè)電源出力,進(jìn)一步消除頻率偏差。常規(guī)二次調(diào)頻采用低帶寬通信,速度較慢集中PI控制器實(shí)現(xiàn)。

    通常,獨(dú)立微電網(wǎng)頻率調(diào)節(jié)由儲能和柴油機(jī)承擔(dān)。而光伏、風(fēng)電等單元受發(fā)電隨機(jī)性影響,常采用最大功率控制,不參與微電網(wǎng)頻率調(diào)節(jié)。

    近年來,風(fēng)、光主動參與系統(tǒng)頻率調(diào)節(jié)的研究逐漸興起。風(fēng)、光發(fā)電采用減載控制,預(yù)留部分有功功率。當(dāng)系統(tǒng)頻率變化時(shí),可靈活調(diào)節(jié)該部分預(yù)留功率,改善系統(tǒng)頻率響應(yīng)。但是,目前風(fēng)、光減載調(diào)頻策略主要集中在一次調(diào)頻策略。受風(fēng)、光發(fā)電隨機(jī)性的影響,其備用功率處在實(shí)時(shí)變化中,給風(fēng)、光參與二次調(diào)頻增加難度。

    1.2微電網(wǎng)狀態(tài)空間模型

    主要聚焦考慮風(fēng)、光減載下的微電網(wǎng)二次調(diào)頻控制。結(jié)合圖1,微電網(wǎng)的負(fù)載頻率控制(load frqeuncycontrol,LFC)模型如圖2所示。圖2中,Tpv?為光伏陣列時(shí)間常數(shù),Tpv?為光伏逆變器控制時(shí)間常數(shù),△Upv為光伏控制輸入?yún)⒖脊β?;Tw為風(fēng)力機(jī)時(shí)間常數(shù),△Uw為風(fēng)力機(jī)控制輸入;R和Tou?分別為柴油機(jī)下垂系數(shù)和調(diào)速時(shí)間常數(shù),Tpu?為柴油機(jī)的時(shí)間常數(shù)。Tgess為儲能系統(tǒng)的時(shí)間常數(shù),R為柴油機(jī)下垂系數(shù)。上述參數(shù)具體值如表2所示。

    微電網(wǎng)中,分布式光伏、風(fēng)電和儲能單元采用并網(wǎng)逆變器接口接入電網(wǎng)。通過并網(wǎng)逆變器控制,DG單元可快速跟蹤給定參考功率。在LFC模型中,為了簡化,采用一階模型表示并網(wǎng)逆變器模型。關(guān)于更多LFC模型的細(xì)節(jié),可參考文獻(xiàn)[19]。

    進(jìn)一步,圖2微電網(wǎng)LFC模型整理為狀態(tài)空間形式,如式(1)所示。

    式中:x——狀態(tài)變量;u——控制變量;d——風(fēng)、光出力波動以及負(fù)荷擾動;y——輸出變量;A。、B、C、D?!B續(xù)狀態(tài)方程的狀態(tài)常數(shù)矩陣、控制常數(shù)矩陣,擾動常數(shù)矩陣和輸出常數(shù)矩陣。式(1)具體展開式見附錄。

    2隨機(jī)性的觀測器設(shè)計(jì)

    風(fēng)光的不確定性以及負(fù)荷擾動對于MPC控制器來說是未知擾動d(k),會直接影響MPC的控制效果,所以首先要對未知擾動進(jìn)行觀測處理,將系統(tǒng)的狀態(tài)空間表達(dá)式進(jìn)行解耦,分解為兩部分,一部分為系統(tǒng)中的已知狀態(tài)變量x和未知狀態(tài)變量x?,其中x?為受未知擾動影響的頻率偏差△f(k)。

    將式(1)進(jìn)行離散化處理,可得:

    式中:A、B、C、D——離散狀態(tài)方程的狀態(tài)矩陣、控制矩陣、擾動矩陣和輸出矩陣。上述系數(shù)矩陣通過連續(xù)狀態(tài)方程(式(1))中的系數(shù)矩陣離散化獲得。采用DU Hamel方法進(jìn)行離散化[20]:

    式中:T,——系統(tǒng)采樣時(shí)間,T,=0.2s。

    對于解耦的狀態(tài)變量,定義其中N為任意矩陣,使得業(yè)=[N D]是非奇異矩陣,將業(yè)與離散化之后的式(2)相乘,得:

    式中:

    由此可得到對應(yīng)于已知狀態(tài)變量x?的表達(dá)式為:

    狀態(tài)變量x?可從測量輸出y(k)中獲得,式(7)可表示為一個(gè)線性表達(dá)式。

    在轉(zhuǎn)移矩陣U=[CD F]中,CD是一個(gè)全列秩矩陣,F(xiàn)是一個(gè)任意矩陣,因此U是一個(gè)非奇異矩陣。因此,U-1=[U?U?],將式(6)中的測量方程乘以U-1,可得:

    將式(8)代入式(7)并與式(9)結(jié)合,可得:

    式中:?=A?-A??U?CN——修正狀態(tài)矩陣;C=U?CN——修正測量矩陣;y(k)=U?y(k)——修正測量向量;E=A??U?。

    如果?、C是可觀測的,則可設(shè)計(jì)卡爾曼濾波器。系統(tǒng)可觀測性的存在條件見文獻(xiàn)[26]??柭鼮V波器如式(11)所示。

    其中:L=LU?+E——卡爾曼增益矩陣

    利用式(8)和式(11)即可估計(jì)出系統(tǒng)中所有的狀態(tài)變量,即:

    式中:我(alk)=U?y(l)-U?CN成,(alk)—△f可從測量值和剩余的估計(jì)狀態(tài)x?(k)和測量值y(k)中估計(jì)出來。

    為了識別未知擾動d(k),擾動包括風(fēng)力機(jī)出力、光伏出力和負(fù)荷擾動。將式(12)代入式(6)中的第一個(gè)表達(dá)式,可得:

    3考慮隨機(jī)性的模型預(yù)測控制器設(shè)計(jì)

    3.1 MPC控制器結(jié)構(gòu)

    微電網(wǎng)二次調(diào)頻控制為多輸入單輸出模型,系統(tǒng)存在強(qiáng)耦合。且風(fēng)力機(jī)、光伏出力及負(fù)荷的不確定性波動會引起系統(tǒng)頻率波動。常規(guī)PI控制并不適合該類系統(tǒng)的控制。

    模型預(yù)測控制(MPC)在處理多輸入多輸出耦合系統(tǒng)時(shí),可獲得令人滿意的控制效果。采用MPC控制來實(shí)現(xiàn)微電網(wǎng)二次調(diào)頻控制,其控制結(jié)構(gòu)如圖3所示。將系統(tǒng)輸出和控制器控制量反饋到觀測器,通過式(10)進(jìn)行不確定性觀測。首先將估計(jì)出來的d(k)輸入到模型預(yù)測環(huán)節(jié)中,然后在考慮到約束條件和頻率波動參考值△f,設(shè)置優(yōu)化目標(biāo)進(jìn)行求解,最后對優(yōu)化目標(biāo)進(jìn)行滾動式優(yōu)化,進(jìn)一步輸出控制量。

    圖3中,MPC控制器頻率波動參考值為0,輸出為各個(gè)分布式電源的調(diào)節(jié)功率輸入?yún)⒖贾怠鱑。采用卡爾曼濾波器獲得系統(tǒng)當(dāng)前狀態(tài)量x(klk)。模型預(yù)測控制依據(jù)當(dāng)前狀態(tài)量和系統(tǒng)離散模型,預(yù)測有限時(shí)域內(nèi)系統(tǒng)狀態(tài)量;采用滾動最優(yōu)控制的思想,求解滿足約束下最優(yōu)控制量,輸出△U調(diào)整電源單元出力,改善系統(tǒng)頻率響應(yīng)。

    3.2 MPC控制器設(shè)計(jì)

    3.2.1模型預(yù)測

    考慮到MPC的準(zhǔn)確性以及計(jì)算的復(fù)雜性,設(shè)定預(yù)測時(shí)域p=10,控制時(shí)域m=3。結(jié)合系統(tǒng)離散模型式(2),可得到預(yù)測時(shí)域p內(nèi)的頻率偏差為:

    式中:Y,(k+1k)——預(yù)測時(shí)域內(nèi)第k時(shí)刻頻率偏差,k=1,2,3,…,P。

    3.2.2約束條件

    對于柴油機(jī)和儲能,其調(diào)頻可用功率受其功率容量的約束。

    式中:Ppu,、Ppu,mx、PBEss,和PBEss,mx——柴油機(jī)或儲能單元的功率上、下限值;Pu、Pgss——當(dāng)前時(shí)刻上層調(diào)度功率值。

    風(fēng)、光發(fā)電可用功率受外部環(huán)境和減載水平(d)影響。選取風(fēng)電和光伏減載水平d=10%。

    風(fēng)電機(jī)組發(fā)出功率可簡化表示為[21]:

    式中:Pwr——雙饋風(fēng)力機(jī)功率輸出;p——空氣密度;R——雙饋葉片半徑;v ——風(fēng)速;Cp——雙饋槳葉功率系數(shù),其值與槳距角β和葉尖速比λ有關(guān),當(dāng)β不變時(shí),可通過改變λ來調(diào)整Cp的大小,來獲得風(fēng)能利用系數(shù)最優(yōu)值C-。C隨著槳距角β減小而增大,通過貝茲理論可知道,風(fēng)能利用系數(shù)的極限值為16/27;C?~C?——風(fēng)力機(jī)特性的相關(guān)參數(shù);w——葉片角速度,其值與轉(zhuǎn)子轉(zhuǎn)速w,成正比。

    其風(fēng)力機(jī)的發(fā)出功率Pw可表示為:

    光伏與風(fēng)力機(jī)的控制方式類似,也采用最大功率跟蹤控制策略,其最大功率估算表達(dá)式為21:

    式中:PMpr,pv——光伏陣列最大輸出功率;n,、n?!夥M件的串聯(lián)和并聯(lián)數(shù);P——光伏組件的額定峰值功率;G.——有效太陽輻照度;G——標(biāo)準(zhǔn)太陽輻照度,取值為1000 W/m2;γ——功率溫度系數(shù);T。——光伏組件的工作溫度;T——標(biāo)準(zhǔn)測量溫度,取值為25℃;T?——環(huán)境溫度;Cr——溫度輻射系數(shù);eNocr——光伏組件的額定溫度。

    光伏機(jī)組的發(fā)出功率Ppv可表示為:

    由式(22)和式(21)可知,獲得風(fēng)、光發(fā)電機(jī)組的調(diào)頻可用功率風(fēng)電關(guān)鍵是獲得機(jī)組最大可用功率PMppr。

    具體的風(fēng)力、光伏發(fā)電機(jī)組減載運(yùn)行控制,已有較多文獻(xiàn)討論[2-25],不再贅述。據(jù)此,風(fēng)、光發(fā)電機(jī)組可用調(diào)頻功率約束設(shè)置為:

    3.2.3優(yōu)化目標(biāo)

    當(dāng)頻率出現(xiàn)波動時(shí),MPC需對系統(tǒng)頻率偏差Yp和輸出控制量△U進(jìn)行權(quán)衡,即存在頻率誤差時(shí),控制器發(fā)出控制指令△U使得頻率恢復(fù)到額定值,具體表達(dá)式為:

    式中:α和β——對應(yīng)的懲罰因子;F,——用來懲罰系統(tǒng)頻率偏差;F.——用來懲罰控制輸出,即發(fā)電單元輸出功率。α=1.6596;βw=βpv=0.2894 ,風(fēng)、光取相同權(quán)重。βou=βBss=0.3762,柴、儲單元比風(fēng)、光發(fā)電單元權(quán)重值大,這樣處理可優(yōu)先風(fēng)、光機(jī)組參與二次調(diào)頻。

    3.2.4含約束最優(yōu)化問題求解

    由于約束條件的存在,不能直接得到目標(biāo)函數(shù)的最優(yōu)解。因此,需將變約束的MPC優(yōu)化問題轉(zhuǎn)換成二次規(guī)劃問題。將預(yù)測方程式(14)代入目標(biāo)函數(shù)中,由于MPC只有控制變量U有關(guān),所以對于優(yōu)化問題而言,目標(biāo)函數(shù)簡化為:

    將約束條件轉(zhuǎn)換為不等式形式:

    為了求解二次規(guī)劃形式的最優(yōu)問題,定義新的變量,即:

    優(yōu)化問題就變成了minp'p,由p1p=(Az-b)'(Az-b)的極值條件,求導(dǎo)可得極值解為:

    可知式(31)取得最小值的解,即k時(shí)刻控制序列最優(yōu)解為:

    只取最優(yōu)控制序列中第一個(gè)元素作為輸出。

    3.2.5 MPC控制器流程

    綜上,MPC控制器流程圖如圖4所示。

    4算例分析

    為驗(yàn)證所提出控制策略的正確性,依照圖1,建立孤島型風(fēng)光柴儲微電網(wǎng)仿真模型??刂品椒ǚ謩e采用PI控制(只有柴、儲參與)、PI控制(風(fēng)、光、柴、儲共同參與)和所提MPC控制3種方式。

    4.1風(fēng)速、光照不變,負(fù)荷突變

    假定風(fēng)速,光照恒定,負(fù)荷階躍突變仿真結(jié)果見圖5,其中實(shí)線代表實(shí)際負(fù)荷波動,虛線代表對于擾動的觀測辨識。

    圖5b~圖5e中,兩種PI控制的二次調(diào)頻出力按機(jī)組容量分配,各個(gè)發(fā)電單元出力波形形狀一致。相比于PI二次調(diào)頻控制,減載+PI控制中,新能源機(jī)組參與調(diào)頻,改善了頻率效果。進(jìn)一步,所提MPC控制方法并結(jié)合約束和權(quán)重,通過滾動優(yōu)化協(xié)調(diào)各個(gè)發(fā)電單元出力。由于光伏和風(fēng)機(jī)出力懲罰因子小于柴油機(jī)和儲能。因此光伏和風(fēng)機(jī)調(diào)頻出力大于減載+PI控制方法,柴油機(jī)和儲能調(diào)頻出力相對應(yīng)小,優(yōu)先新能源機(jī)組調(diào)頻。

    圖5a中,虛線代表預(yù)測的負(fù)荷波動,當(dāng)負(fù)荷變化時(shí),由于預(yù)測需系統(tǒng)輸出值,所以導(dǎo)致預(yù)測會有一定的時(shí)滯性,但是最終會與實(shí)際值一致。圖5f中,相對于兩種PI控制方法,所提出的MPC控制器的頻率波動更小,收斂速度更快。尤其是在90s時(shí)刻,系統(tǒng)出現(xiàn)大擾動,3種控制器的最大頻率偏差分別是1.926×10-3、1.42×10-3和2.453×10-?pu,標(biāo)準(zhǔn)差為5.194×10-3、3.595×10-3、3.788×10?pu。

    4.2風(fēng)速、光照和負(fù)荷波動場景

    實(shí)際工作中,風(fēng)速、光照和系統(tǒng)負(fù)荷處在波動狀態(tài)。本場景將3 min實(shí)際風(fēng)速、光照數(shù)據(jù)和負(fù)荷波動數(shù)據(jù)導(dǎo)入仿真系統(tǒng)。風(fēng)、光伏波動和負(fù)荷波動如圖6a所示。該情景下的仿真結(jié)果見圖6b~圖6f,其中實(shí)線代表實(shí)際功率波動,虛線代表對于擾動的觀測辨識;星形曲線表示負(fù)荷擾動,圓型曲線表示光伏功率波動,角形曲線表示風(fēng)力機(jī)功率波動。

    受風(fēng)、光伏波動性的影響,其可用調(diào)頻功率限值也出現(xiàn)波動,見圖6b~圖6c。兩種PI控制方法中,各個(gè)發(fā)電單元的出力形狀一致,見圖6b~圖6e。而MPC控制,由于風(fēng)、光懲罰因子更小,優(yōu)先風(fēng)、光發(fā)電單元參與調(diào)頻。圖6a中,當(dāng)考慮風(fēng)、光波動較平穩(wěn)時(shí),觀測器所預(yù)測的值與實(shí)際值基本一致。圖6b~圖6c的部分時(shí)間段,風(fēng)、光調(diào)頻出力達(dá)到限值。相應(yīng)的,柴油機(jī)、儲能調(diào)頻出力得到減少,提高微電網(wǎng)經(jīng)濟(jì)運(yùn)行。

    3種控制方法中,MPC二次調(diào)頻控制方法最優(yōu),如圖8f所示。3種方法最大頻率偏差分別是2.885×10-3、1.541×10-3和7.732×10-?pu,頻率波動標(biāo)準(zhǔn)差分別是2.477×10-3、7.392×10-3和1.385×10-2pu。

    4.3風(fēng)速、光照快速變化場景

    某些惡劣天氣下,風(fēng)速、光照可能會處于快速變化中。例如云層的快速遮蔽、風(fēng)速的突然變化等。負(fù)荷波動與場景2保持一致,但風(fēng)速,光照選取實(shí)測數(shù)據(jù)中快速變化的部分。本場景的風(fēng)速、光照和負(fù)荷波動曲線如圖7a所示,其中實(shí)線代表實(shí)際功率波動,虛線代表對于擾動的觀測辨識;星形曲線表示負(fù)荷擾動,圓型曲線表示光伏功率波動,角形曲線表示風(fēng)力機(jī)功率波動。

    受光照、風(fēng)速快速變化的影響。風(fēng)電,光伏單元的可用調(diào)頻功率存在劇烈波動動作。常規(guī)PI控制只依照機(jī)組額定容量,分配調(diào)頻功率。PI+減載控制,部分區(qū)域的風(fēng)、光機(jī)組的功率達(dá)到限值。由于飽和限值消減PI控制器部分輸出,導(dǎo)致這些區(qū)域段的頻率波動增大。

    相比較,提出的MPC控制方法,由于優(yōu)先光伏、風(fēng)電出力,部分區(qū)域風(fēng)、光發(fā)電單元調(diào)頻出力也達(dá)到上限值。但是,MPC控制是在考慮功率約束下的最優(yōu)控制輸出。因此,這些區(qū)域段,柴油、儲能輸出功率會適當(dāng)提高,來保證頻率調(diào)節(jié)性能。因此,相比于場景2,場景3中MPC二次調(diào)頻方法的頻率響應(yīng)得到更明顯程提升。3種控制方法的最大頻率偏差分別為4.563×10-3、2.970×10-3和1.584×10-3pu,標(biāo)準(zhǔn)差分別為2.6×10-2、1.818×10-2和5.934×10-3pu。

    為了更加直觀的體現(xiàn)所提出控制器在二次調(diào)頻控制響應(yīng)性能,表3描述了3種控制器的最大頻率偏差以及標(biāo)準(zhǔn)差。

    5結(jié)論

    微電網(wǎng)二次調(diào)頻控制是個(gè)多輸入單輸出且輸入存在隨機(jī)性的系統(tǒng)控制問題。針對隨機(jī)性的問題,引入不確定輸入擾動觀測,對未知負(fù)荷和發(fā)電擾動進(jìn)行卡爾曼濾波觀測;增加風(fēng)、光伏發(fā)電實(shí)時(shí)約束,避免機(jī)組調(diào)頻超出功率限額;通過權(quán)重系數(shù)設(shè)置,優(yōu)先風(fēng)、光發(fā)電參與二次調(diào)頻。所提MPC控制器在線優(yōu)化調(diào)整各個(gè)分布式電源的出力,提供微電網(wǎng)的頻率控制能力。建立了含多臺風(fēng)、光發(fā)電的獨(dú)立微電網(wǎng)二次調(diào)頻模型。風(fēng)、光伏恒定,典型風(fēng)、光伏波動和劇烈風(fēng)、光波動場景下,對所提方法進(jìn)行了對比驗(yàn)證。結(jié)果顯示所提方法能加快系統(tǒng)的頻率恢復(fù)速度,減小頻率的波動峰值,尤其在劇烈風(fēng)、光波動場景下,頻率調(diào)節(jié)提升效果更為明顯。

    [參考文獻(xiàn)]

    [1]陳國平,李明節(jié),許濤,等.我國電網(wǎng)支撐可再生能源發(fā)展的實(shí)踐與挑戰(zhàn)[J].電網(wǎng)技術(shù),2017,41(10):3095-3103.

    CHEN GP,LIM J,XU T,etal.Practice and challengeof renewable energy development based on interconnectedpower grids[J].Power system technology,2017,41(10):3095-3103.

    [2]丁圣康,劉璐裔.微電網(wǎng)研究綜述[J].大眾標(biāo)準(zhǔn)化,2021(22):225-227.

    DING S K,LIU L Y.Review of microgrid research [J].Popular standardization,2021(22):225-227.

    [3]畢永健,徐丙垠,趙艷雷,等.同步定頻微電網(wǎng)的并網(wǎng)/孤島無縫切換控制策略[J].電網(wǎng)技術(shù),2022,46(3):923-933.

    BI Y J,XU B Y,ZHAO YL,etal.Seamlesstransfercontrol strategy between grid-connected and islandingoperation for synchronous fixed-frequency microgrid [J].Power system technology,2022,46(3):923-933.

    [4]蘭征,刁偉業(yè),涂春鳴,等.含儲能和氫燃料電池的孤島微電網(wǎng)混合運(yùn)行模式與功率協(xié)調(diào)策略研究[J].電網(wǎng)技術(shù),2022,46(1):156-164.

    LAN Z,DIAO WY,TU C M,etal.Research on hybridoperation mode and power coordination strategy of islandmicrogrid with energy storage and hydrogen fuel cell[J].Power system technology,2022,46(1):156-164.

    [5]許志榮,楊蘋,鄭成立,等.孤島型風(fēng)柴儲微電網(wǎng)運(yùn)行情況分析[J].電網(wǎng)技術(shù),2016,40(7):1978-1984.

    XUZ R,YANG P,ZHENG C L,etal.Operationanalysisof isolated microgrid including wind turbine,dieselgenerator and battery storage[J].Power system technology,2016,40(7):1978-1984.

    [6]胡家欣,胥國毅,畢天姝,等.減載風(fēng)電機(jī)組變速變槳協(xié)調(diào)頻率控制方法[J].電網(wǎng)技術(shù),2019,43(10):3656-3663.

    HU JX,XU G Y,BI T S,etal.A strategy of frequencycontrol for deloaded wind turbine generator based oncoordination between rotor speed and pitch angle [J].Power system technology,2019,43(10):3656-3663.

    [7]劉漂慶,吳寧.張煥早.等微申網(wǎng)經(jīng)濟(jì)型一次頻率和電壓控制的多目標(biāo)優(yōu)化模型及仿真驗(yàn)證[J].電網(wǎng)技術(shù),2019,43(2):521-530.

    LIU L Q,WU N,ZHANG H H,etal.Multi-objectiveoptimization model and its simulation verification foreconomical secondary frequency and voltage controls ofmicrogrids[J].Power system technology,2019,43(2):521-530

    [8]施永,徐冬,于鴻儒,等.基于系統(tǒng)辨識建模的微網(wǎng)二次電壓頻率控制器參數(shù)設(shè)計(jì)方法[J].電力系統(tǒng)自動化,2020,44(13):89-97.

    SHIY,XUD,YUHR,etal.Parameter design method ofsecondary voltage and frequency regulation controller inmicrogrid based on system identification modeling [J].Automation ofelectric power systems,2020,44(13):89-97.

    [9]MU CX,ZHANG Y,JIA HJ,etal.Energy-storage-basedintelligent frequency control of microgrid with stochasticmodel uncertainties[J].IEEE transactions on smart grid,2020,11(2):1748-1758.

    [10]李斌,周林,余希瑞,等.基于改進(jìn)虛擬同步發(fā)電機(jī)算法的微網(wǎng)逆變器二次調(diào)頻方案[J].電網(wǎng)技術(shù),2017,41(8):2680-2687.

    LI B,ZHOU L,YUX R,etal.Secondaryfrequencyregulation for microgrid inverters based on improvingvirtual synchronous generator[J].Power system technology,2017,41(8):2680-2687.

    [11]MA L,ZHANG X D,RAO YQ,etal.A robust modelpredictive control for multi-microgrids frequency control[C]//2020 IEEE 4th Conference on Energy Internet andEnergy System Integration(EI2).Wuhan,China,2021:1637-1641.

    [12]GBADEGA PA,SAHA A K.Load frequency control of atwo-area power system with a stand-alone microgrid basedon adaptive model predictive control[J].IEEE journal ofemerging and selected topics in power electronics,2021,9(6):7253-7263.

    [13]LONG B,LIAO Y,CHONG KT,etal.Enhancement of frequency regulation in AC microgrid:a fuzzy-MPCcontrolled virtual synchronous generator[J].IEEEtransactions on smart grid,2021,12(4):3138-3149.

    [14]MORSTYN T,HREDZAK B,AGUILERA R P,etal.Model predictive control for distributed microgrid batteryenergy storage systems [J].IEEE transactions on controlsystems technology,2018,26(3):1107-1114.

    [15]李帥虎,向麗珍,向振宇,等.用于改善VSG頻率響應(yīng)的模型預(yù)測控制方法[J].高電壓技術(shù),2021,47(8):2856-2864.

    LI S H,XIANG LZ,XIANG Z Y,etal.MPCcontrolmethod for improving VSG frequency response [J].Highvoltage engineering,2021,47(8):2856-2864.

    [16]WANG Y,DELILLE G,BAYEM H,etal.Highwindpower penetration in isolated power systems-assessment ofwind inertial and primary frequency responses [J].IEEEtransactions on power systems,2013,28(3):2412-2420.

    [17]吳巖,王瑋,曾國宏,等.四開關(guān)Buck-Boost變換器的多模式模型預(yù)測控制策略[J].電工技術(shù)學(xué)報(bào),2022,37(10):2572-2583.

    WU Y,WANG W,ZENG G H,etal.Multi-mode modelpredictive control strategy for the four-switch buck-boostconverter[J].Transactions of China ElectrotechnicalSociety,2022,37(10):2572-2583.

    [18]KROPOSKI B,JOHNSON B,ZHANG Y C,etal.Achieving a 100%renewable grid:operatingelectricpower systems with extremely high levels of variablerenewable energy [J].IEEE power and energy magazine,2017,15(2):61-73.

    [19]楊海濤,江晶晶,趙敏,等.基于模型預(yù)測控制的區(qū)域綜合能源系統(tǒng)運(yùn)行優(yōu)化方法[J].電氣技術(shù),2022,23(4):7-13.

    YANG HT,JIANG J J,ZHAO M,etal.Operationaloptimization method of regional integrated energy systembased on model predictive control[J].Electrical engineering,2022,23(4):7-13.

    [20]CHEN XT,WU WM,GAO N,etal.Finite control setmodel predictive control for LCL-filtered grid-tied inverter with minimum sensors [J].IEEE transactions on industrialelectronics,2020,67(12):9980-9990.

    [21]ZHANG Q,LI Y,DING Z W,etal.Self-adaptivesecondary frequency regulation strategy of micro-grid withmultiple virtual synchronous generators[J].IEEEtransactions on industry applications,2020,56(5):6007-6018.

    [22]張偉亮,張輝,支娜,等.基于節(jié)點(diǎn)源荷電流差分的直流微電網(wǎng)儲能變換器控制策略[J].電工技術(shù)學(xué)報(bào),2022,37(9):2199-2210.

    ZHANG WL,ZHANG H,ZHI N,etal.Controlstrategyof DC microgrid energy storage converter based on nodedifferential current[J].Transactions of ChinaElectrotechnical Society,2022,37(9):2199-2210.

    [23]顏湘武,徐韻,李若瑾,等.基于模型預(yù)測控制含可再生分布式電源參與調(diào)控的配電網(wǎng)多時(shí)間尺度無功動態(tài)優(yōu)化[J].電工技術(shù)學(xué)報(bào),2019,34(10):2022-2037.

    YANXW,XU Y,LI R J,etal.Multi-time scale reactivepower optimization of distribution grid based on modelpredictive control and including RDG regulation[J].Transactions of China Electrotechnical Society,2019,34(10):2022-2037.

    [24]馬偉棟,高丙朋,楊武幫,等.引入權(quán)值修正預(yù)測控制的風(fēng)電葉片自適應(yīng)組合抑振策略研究[J].太陽能學(xué)報(bào),2022,43(8):382-390.

    MA W D,GAO B P,YANG W B,etal.Researchonadaptive combined vibration suppression strategy of windpower blade based on weight modified predictive control[J].Acta energiaesolaris sinica,2022,43(8):382-390.

    [25]蔡鵬程,文傳博.基于比率一致性算法的孤島微電網(wǎng)分布式二次頻率控制[J].太陽能學(xué)報(bào),2020,41(10):74-81.

    CAI P C,WEN C B.Secondary frequency control ofislanded microgrids based on ratio consensue algorithm[J].Acta energiaesolaris sinica,2020,41(10):74-81.

    [26]KHALGHANI M R,SOLANKI J,KHUSHALANISOLANKI S,etal.Stochastic secondary frequency controlof islanded microgrid under uncertainties[J].IEEE systemsjournal,2021,15(1):1056-1065.

    SECONDARY FREQUENCY CONTROL OF ISLANDED MICROGRIDCONSIDERING WIND AND SOLAR STOCHASTICS

    Zhong Cheng1,2,Jiang Zhifu2,Zhang Xiangyu2,Chen Jikai1,2,Li Yang2

    (1.Key Laboratory of Modern Power System Simulation and Control amp;Renewable Energy Technology,Ministry of Education,Northeast Electric Power University,Jilin 132012,China;

    2.School of Electrical Engineering,Northeast Electric Power University,Jilin 132012,China)

    Abstract:This paper proposed a model predictive control(MPC)secondary frequency control method considering wind and solarpower generation stochastics.The extended state-space matrix including unknown stochastic power disturbance is established,andaKalmanflter is used to observe the unknown disturbance.The maximum available power of wind and solar DGs is estimated forestablishing real-time variable constraints that prevent DGs output power from exceding the limits.Throughseting proper weightcoeffcients,wind and photovoltaic DGs are given priority to participate in secondary frequency control.The distributed restorativepower of each DG is obtained by solving the quadratic programming(QP)optimal problem with variable constraints.Finally,amicrogridsimulation model including multiple PV and wind DGs is built and performed in various scenarios compared to the traditional secondaryfrequency control method.The simulation results validated that the proposed method can enhance the frequency recovery speed andreduce the frequency deviation,especially in severe photovoltaic and wind fluctuations scenarios.

    Keywords:secondary frequency eontrol;model predictive control;islandedmicrogrid;deloadingcontrol;stochastic input observer

    附錄

    式中:△f——頻率偏差;△x?——光伏陣列輸出功率;△Py——光伏逆變器輸出功率;△Pw——風(fēng)力機(jī)輸出功率;△x?——調(diào)速器增量;△Pou——柴油機(jī)輸出功率;△PBss——儲能系統(tǒng)輸出功率;H,——等效慣性。

    猜你喜歡
    模型預(yù)測控制
    基于預(yù)測控制綜合算法的控制器設(shè)計(jì)與仿真
    科技視界(2017年35期)2018-03-09 20:37:34
    室溫預(yù)測控制的仿真研究
    基于虛擬矢量的FCS—MPC在單相PWM整流器中的應(yīng)用
    科技視界(2017年23期)2017-12-09 07:56:32
    核電站蒸汽發(fā)生器水位模型預(yù)測控制方法研究
    預(yù)測控制在電力電子技術(shù)中的應(yīng)用
    三相電壓型PWM整流器直接功率控制策略研究
    三相電壓型PWM整流器直接功率控制策略研究
    AECS—2000在機(jī)組AGC控制中的應(yīng)用
    非完整移動機(jī)器人的模型預(yù)測控制技術(shù)的研究
    科技資訊(2015年4期)2015-07-02 17:55:56
    飽和離散隨機(jī)非線性系統(tǒng)的H模型預(yù)測控制
    国产免费又黄又爽又色| 啦啦啦在线观看免费高清www| 国产白丝娇喘喷水9色精品| a级毛色黄片| www.熟女人妻精品国产 | 国产日韩一区二区三区精品不卡| 欧美日韩综合久久久久久| 亚洲少妇的诱惑av| 2021少妇久久久久久久久久久| 欧美精品国产亚洲| 亚洲av欧美aⅴ国产| 80岁老熟妇乱子伦牲交| 久久精品久久久久久久性| 人人妻人人添人人爽欧美一区卜| 男人操女人黄网站| av有码第一页| 精品人妻一区二区三区麻豆| 国产精品久久久久成人av| 一级片'在线观看视频| 天天操日日干夜夜撸| 久久精品aⅴ一区二区三区四区 | 日韩大片免费观看网站| 国产熟女欧美一区二区| 精品人妻一区二区三区麻豆| 老司机影院成人| www.熟女人妻精品国产 | 少妇人妻精品综合一区二区| www.熟女人妻精品国产 | 男人爽女人下面视频在线观看| 爱豆传媒免费全集在线观看| 国产福利在线免费观看视频| 久久久精品免费免费高清| 国产69精品久久久久777片| 色婷婷久久久亚洲欧美| 欧美国产精品一级二级三级| 只有这里有精品99| 少妇高潮的动态图| 丁香六月天网| 日韩 亚洲 欧美在线| 国产综合精华液| 搡女人真爽免费视频火全软件| 日本-黄色视频高清免费观看| 午夜免费观看性视频| 免费看av在线观看网站| 久久韩国三级中文字幕| 中国三级夫妇交换| 亚洲精品乱码久久久久久按摩| 人妻一区二区av| 丝瓜视频免费看黄片| 亚洲av免费高清在线观看| 一级毛片电影观看| 国产一区二区三区av在线| 大陆偷拍与自拍| 久久精品熟女亚洲av麻豆精品| 国产精品一区www在线观看| 欧美性感艳星| av福利片在线| 永久网站在线| 亚洲国产精品999| 国产精品久久久久久精品古装| 卡戴珊不雅视频在线播放| 一级片'在线观看视频| 中文欧美无线码| 亚洲精品乱码久久久久久按摩| 伦精品一区二区三区| 熟妇人妻不卡中文字幕| www日本在线高清视频| 亚洲av电影在线观看一区二区三区| 免费观看a级毛片全部| 少妇的丰满在线观看| 久久免费观看电影| 在线精品无人区一区二区三| 久久久久久久亚洲中文字幕| 少妇 在线观看| 卡戴珊不雅视频在线播放| 午夜福利影视在线免费观看| 国国产精品蜜臀av免费| 国产精品一国产av| 国产有黄有色有爽视频| 宅男免费午夜| 久久精品国产a三级三级三级| 男人爽女人下面视频在线观看| 校园人妻丝袜中文字幕| 桃花免费在线播放| 国产黄色免费在线视频| 欧美日韩综合久久久久久| 欧美老熟妇乱子伦牲交| 国产有黄有色有爽视频| 伦理电影大哥的女人| 国产精品人妻久久久久久| 久久久久久人人人人人| av国产久精品久网站免费入址| 免费久久久久久久精品成人欧美视频 | 亚洲国产欧美日韩在线播放| 一本大道久久a久久精品| 午夜激情久久久久久久| 香蕉国产在线看| 久久人妻熟女aⅴ| 国产不卡av网站在线观看| 久久久国产欧美日韩av| 亚洲国产精品999| 久久久久久久亚洲中文字幕| 另类精品久久| 久久国内精品自在自线图片| av天堂久久9| 女人精品久久久久毛片| 日产精品乱码卡一卡2卡三| 夫妻性生交免费视频一级片| 午夜免费男女啪啪视频观看| av黄色大香蕉| 妹子高潮喷水视频| 亚洲婷婷狠狠爱综合网| 国产视频首页在线观看| 久久久久视频综合| 国产极品粉嫩免费观看在线| 欧美精品人与动牲交sv欧美| 国产免费现黄频在线看| 久久这里有精品视频免费| 免费黄网站久久成人精品| 亚洲欧美中文字幕日韩二区| 好男人视频免费观看在线| 欧美日韩视频高清一区二区三区二| 精品久久国产蜜桃| 美女国产视频在线观看| 少妇被粗大猛烈的视频| 国产极品天堂在线| 满18在线观看网站| 日韩熟女老妇一区二区性免费视频| 欧美老熟妇乱子伦牲交| 在线观看免费视频网站a站| 精品亚洲成国产av| 99久久人妻综合| 欧美亚洲日本最大视频资源| 亚洲精品第二区| 国产精品熟女久久久久浪| videossex国产| 久久国产亚洲av麻豆专区| 亚洲国产精品一区三区| 男女边吃奶边做爰视频| 丰满乱子伦码专区| 国产国拍精品亚洲av在线观看| 久久这里只有精品19| 一区在线观看完整版| 免费观看av网站的网址| 自拍欧美九色日韩亚洲蝌蚪91| 人妻少妇偷人精品九色| 亚洲伊人久久精品综合| 国产欧美亚洲国产| 国产成人一区二区在线| 国产精品一区www在线观看| 99热这里只有是精品在线观看| 在线观看人妻少妇| 成人二区视频| 欧美精品亚洲一区二区| 黄色怎么调成土黄色| 成人国产麻豆网| 午夜视频国产福利| 大香蕉97超碰在线| 午夜激情久久久久久久| 亚洲成国产人片在线观看| 国产麻豆69| 免费高清在线观看视频在线观看| 十八禁网站网址无遮挡| 美女中出高潮动态图| 丁香六月天网| 欧美丝袜亚洲另类| 免费日韩欧美在线观看| 91久久精品国产一区二区三区| 国产极品粉嫩免费观看在线| 亚洲欧美精品自产自拍| 熟女av电影| 欧美3d第一页| av不卡在线播放| 十分钟在线观看高清视频www| av在线观看视频网站免费| 久久精品人人爽人人爽视色| 国产又色又爽无遮挡免| 9热在线视频观看99| av卡一久久| 久久ye,这里只有精品| 国产亚洲最大av| 我的女老师完整版在线观看| 在线观看免费日韩欧美大片| 中文字幕另类日韩欧美亚洲嫩草| 深夜精品福利| 91成人精品电影| 欧美xxxx性猛交bbbb| 9热在线视频观看99| 9191精品国产免费久久| 欧美精品av麻豆av| 亚洲精品乱久久久久久| 中文欧美无线码| 满18在线观看网站| 国产精品.久久久| 午夜免费观看性视频| 亚洲av在线观看美女高潮| 午夜福利乱码中文字幕| 亚洲少妇的诱惑av| 国产日韩欧美视频二区| av不卡在线播放| 国产精品久久久久久av不卡| av一本久久久久| 观看av在线不卡| 91午夜精品亚洲一区二区三区| 成人亚洲欧美一区二区av| 免费在线观看完整版高清| 亚洲三级黄色毛片| 久久精品国产鲁丝片午夜精品| 欧美xxxx性猛交bbbb| 丰满少妇做爰视频| 国产精品久久久久成人av| av在线播放精品| 久久这里只有精品19| 亚洲精品,欧美精品| 欧美日韩精品成人综合77777| 黄色一级大片看看| 国产亚洲av片在线观看秒播厂| 性色av一级| 免费久久久久久久精品成人欧美视频 | 99久久精品国产国产毛片| 日韩免费高清中文字幕av| 亚洲高清免费不卡视频| 国产激情久久老熟女| 一级毛片 在线播放| 一区二区av电影网| 男人操女人黄网站| 日韩av在线免费看完整版不卡| 日本-黄色视频高清免费观看| 大陆偷拍与自拍| 国产精品.久久久| 亚洲精品美女久久av网站| 中国美白少妇内射xxxbb| 少妇被粗大的猛进出69影院 | 超碰97精品在线观看| 国产乱来视频区| 水蜜桃什么品种好| 久久精品国产综合久久久 | www日本在线高清视频| 成人影院久久| 精品午夜福利在线看| 亚洲av男天堂| 人人妻人人澡人人看| 观看av在线不卡| 美女脱内裤让男人舔精品视频| 黄色配什么色好看| 十八禁网站网址无遮挡| 久久亚洲国产成人精品v| 亚洲第一av免费看| 老熟女久久久| 下体分泌物呈黄色| 青青草视频在线视频观看| 一级片'在线观看视频| 一级毛片 在线播放| 久久久久久伊人网av| 国产黄频视频在线观看| 视频在线观看一区二区三区| 大话2 男鬼变身卡| 精品熟女少妇av免费看| 免费高清在线观看日韩| 美女脱内裤让男人舔精品视频| 国产在线免费精品| 免费看光身美女| 免费看不卡的av| 老熟女久久久| 丝瓜视频免费看黄片| 久久久久久久国产电影| 久久99热6这里只有精品| 超碰97精品在线观看| 丰满饥渴人妻一区二区三| 国产精品嫩草影院av在线观看| 精品少妇内射三级| 亚洲av电影在线观看一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 一区二区三区乱码不卡18| 久久午夜福利片| 最近中文字幕高清免费大全6| 久久精品久久精品一区二区三区| 欧美精品av麻豆av| 久久影院123| 又粗又硬又长又爽又黄的视频| 色吧在线观看| 男的添女的下面高潮视频| 午夜福利在线观看免费完整高清在| 国产淫语在线视频| 国产av一区二区精品久久| 久久狼人影院| av黄色大香蕉| 日韩成人伦理影院| 尾随美女入室| 亚洲在久久综合| 伦理电影大哥的女人| 国产av精品麻豆| 亚洲国产精品一区二区三区在线| 交换朋友夫妻互换小说| 天天躁夜夜躁狠狠久久av| 亚洲欧洲国产日韩| 日本vs欧美在线观看视频| 五月玫瑰六月丁香| 有码 亚洲区| av一本久久久久| 美女xxoo啪啪120秒动态图| 天美传媒精品一区二区| 2018国产大陆天天弄谢| 国产精品免费大片| 最新中文字幕久久久久| 日韩不卡一区二区三区视频在线| 日本色播在线视频| 天堂中文最新版在线下载| 亚洲精品aⅴ在线观看| 人成视频在线观看免费观看| 亚洲激情五月婷婷啪啪| 亚洲一区二区三区欧美精品| 免费高清在线观看视频在线观看| 如日韩欧美国产精品一区二区三区| 日韩av不卡免费在线播放| 欧美日韩成人在线一区二区| av福利片在线| 亚洲欧美成人综合另类久久久| 搡女人真爽免费视频火全软件| 日韩三级伦理在线观看| 久久久精品94久久精品| 成人18禁高潮啪啪吃奶动态图| 制服人妻中文乱码| 日韩成人av中文字幕在线观看| 黑人欧美特级aaaaaa片| 久久人人爽av亚洲精品天堂| 日本爱情动作片www.在线观看| 黄色毛片三级朝国网站| 王馨瑶露胸无遮挡在线观看| 亚洲欧美清纯卡通| 精品一区二区三区视频在线| 国产成人免费无遮挡视频| 性色avwww在线观看| 97人妻天天添夜夜摸| 成人国产麻豆网| 最近最新中文字幕大全免费视频 | 一级毛片电影观看| 色视频在线一区二区三区| 国产成人91sexporn| 亚洲一级一片aⅴ在线观看| 少妇的逼好多水| 日本av手机在线免费观看| 国精品久久久久久国模美| 91精品伊人久久大香线蕉| 丰满少妇做爰视频| 国产精品嫩草影院av在线观看| 国产日韩欧美视频二区| 丝袜在线中文字幕| 久久精品国产a三级三级三级| 日韩视频在线欧美| 欧美变态另类bdsm刘玥| 亚洲成av片中文字幕在线观看 | 肉色欧美久久久久久久蜜桃| 久久免费观看电影| 哪个播放器可以免费观看大片| 久久久久久伊人网av| 高清av免费在线| 中文精品一卡2卡3卡4更新| 免费av中文字幕在线| 日韩成人av中文字幕在线观看| 丰满少妇做爰视频| 侵犯人妻中文字幕一二三四区| 国产爽快片一区二区三区| 亚洲欧美色中文字幕在线| 亚洲久久久国产精品| 久久国产精品男人的天堂亚洲 | 日韩人妻精品一区2区三区| 亚洲精品,欧美精品| av又黄又爽大尺度在线免费看| 色哟哟·www| 看免费av毛片| 亚洲欧美清纯卡通| 建设人人有责人人尽责人人享有的| 如何舔出高潮| 欧美成人午夜免费资源| 韩国高清视频一区二区三区| 飞空精品影院首页| 又黄又爽又刺激的免费视频.| 两个人免费观看高清视频| 蜜桃在线观看..| 久久久久久久久久人人人人人人| 男的添女的下面高潮视频| 亚洲天堂av无毛| 亚洲成国产人片在线观看| av国产精品久久久久影院| 欧美精品亚洲一区二区| 日韩制服丝袜自拍偷拍| 日韩制服丝袜自拍偷拍| 男女无遮挡免费网站观看| 久久久久精品人妻al黑| 九九在线视频观看精品| av网站免费在线观看视频| 最新的欧美精品一区二区| 99香蕉大伊视频| 国产亚洲午夜精品一区二区久久| 夫妻午夜视频| 亚洲成人手机| 国产精品久久久久久久久免| 国产淫语在线视频| 日本av免费视频播放| 亚洲伊人久久精品综合| 国产视频首页在线观看| 亚洲伊人久久精品综合| 美女国产视频在线观看| 精品第一国产精品| 女性被躁到高潮视频| 极品少妇高潮喷水抽搐| 99re6热这里在线精品视频| 免费观看a级毛片全部| 色哟哟·www| 中文字幕精品免费在线观看视频 | 亚洲久久久国产精品| 亚洲一码二码三码区别大吗| 一区二区av电影网| 青青草视频在线视频观看| 国产精品一区www在线观看| 国产乱人偷精品视频| 国产精品.久久久| 91aial.com中文字幕在线观看| 国产欧美另类精品又又久久亚洲欧美| xxx大片免费视频| 热re99久久国产66热| 99热全是精品| 欧美97在线视频| 国产有黄有色有爽视频| 国产极品粉嫩免费观看在线| 在线观看国产h片| 免费日韩欧美在线观看| 香蕉精品网在线| 在线精品无人区一区二区三| 超色免费av| 国产精品一区二区在线不卡| 韩国av在线不卡| 18+在线观看网站| 22中文网久久字幕| 一级a做视频免费观看| 下体分泌物呈黄色| 看非洲黑人一级黄片| 中国国产av一级| 色婷婷久久久亚洲欧美| tube8黄色片| 免费av不卡在线播放| 欧美bdsm另类| 亚洲 欧美一区二区三区| 免费看av在线观看网站| 国产午夜精品一二区理论片| 美女大奶头黄色视频| 一边亲一边摸免费视频| 日韩免费高清中文字幕av| 乱人伦中国视频| 看免费成人av毛片| 国产熟女午夜一区二区三区| 欧美日韩视频高清一区二区三区二| 99热这里只有是精品在线观看| 91国产中文字幕| 亚洲色图综合在线观看| 伊人久久国产一区二区| 国产 一区精品| 亚洲精品日本国产第一区| 久久av网站| 18禁裸乳无遮挡动漫免费视频| 国产69精品久久久久777片| 这个男人来自地球电影免费观看 | 国产一区亚洲一区在线观看| 人妻系列 视频| 九草在线视频观看| 国产精品免费大片| 人妻系列 视频| 欧美最新免费一区二区三区| 成人国产av品久久久| av国产精品久久久久影院| 乱码一卡2卡4卡精品| 22中文网久久字幕| 天堂8中文在线网| 日本91视频免费播放| 一区二区三区精品91| 高清不卡的av网站| 伊人久久国产一区二区| 人成视频在线观看免费观看| 制服人妻中文乱码| 日韩一区二区三区影片| 亚洲av在线观看美女高潮| 国产极品天堂在线| 免费黄色在线免费观看| 亚洲精品美女久久久久99蜜臀 | 免费久久久久久久精品成人欧美视频 | 成年女人在线观看亚洲视频| 国产成人精品福利久久| 成年美女黄网站色视频大全免费| 精品亚洲成a人片在线观看| 巨乳人妻的诱惑在线观看| 丝袜喷水一区| 一区二区日韩欧美中文字幕 | 18禁裸乳无遮挡动漫免费视频| 视频在线观看一区二区三区| 桃花免费在线播放| 久久午夜综合久久蜜桃| 国产精品嫩草影院av在线观看| 九色亚洲精品在线播放| 高清毛片免费看| 十八禁网站网址无遮挡| 成人黄色视频免费在线看| 黑人欧美特级aaaaaa片| 国产亚洲精品久久久com| 一个人免费看片子| 看免费成人av毛片| 22中文网久久字幕| 2018国产大陆天天弄谢| 欧美丝袜亚洲另类| 国产熟女欧美一区二区| 中国国产av一级| 天天躁夜夜躁狠狠躁躁| 精品国产露脸久久av麻豆| 在线观看美女被高潮喷水网站| 国产精品久久久av美女十八| 你懂的网址亚洲精品在线观看| 国产欧美另类精品又又久久亚洲欧美| 久久久久网色| 国产不卡av网站在线观看| 国产精品国产三级国产av玫瑰| 曰老女人黄片| 免费av不卡在线播放| xxxhd国产人妻xxx| 日韩欧美一区视频在线观看| 18禁裸乳无遮挡动漫免费视频| 三上悠亚av全集在线观看| 超色免费av| 熟女人妻精品中文字幕| 日本免费在线观看一区| 亚洲高清免费不卡视频| 久久久久精品人妻al黑| 国产亚洲最大av| 男人添女人高潮全过程视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产欧美日韩在线播放| 久久免费观看电影| 欧美国产精品va在线观看不卡| 国产精品一区二区在线观看99| 18禁裸乳无遮挡动漫免费视频| 久久人人爽人人片av| 国产成人一区二区在线| 最黄视频免费看| 色婷婷久久久亚洲欧美| 18禁动态无遮挡网站| 亚洲国产毛片av蜜桃av| 男女边摸边吃奶| 久久韩国三级中文字幕| 亚洲国产看品久久| 国产免费一级a男人的天堂| 黄片播放在线免费| 一个人免费看片子| 成人手机av| 欧美 亚洲 国产 日韩一| 国产毛片在线视频| av片东京热男人的天堂| 少妇 在线观看| 夫妻性生交免费视频一级片| 色哟哟·www| 久久久久久久大尺度免费视频| 欧美xxⅹ黑人| 天美传媒精品一区二区| 伊人久久国产一区二区| 亚洲精品,欧美精品| 99热全是精品| 亚洲av电影在线进入| 日本爱情动作片www.在线观看| 欧美97在线视频| 久久精品久久久久久噜噜老黄| 亚洲精品日本国产第一区| 精品酒店卫生间| www.熟女人妻精品国产 | 久热这里只有精品99| 国产欧美日韩综合在线一区二区| 中文乱码字字幕精品一区二区三区| 免费看不卡的av| 少妇人妻久久综合中文| 国产日韩欧美视频二区| 中文字幕人妻丝袜制服| 日本av手机在线免费观看| 免费人成在线观看视频色| 视频区图区小说| 麻豆精品久久久久久蜜桃| 一级毛片黄色毛片免费观看视频| 在线观看国产h片| 香蕉精品网在线| 午夜福利在线观看免费完整高清在| 黄色怎么调成土黄色| 亚洲丝袜综合中文字幕| 精品久久久久久电影网| 热re99久久精品国产66热6| 满18在线观看网站| 国产免费视频播放在线视频| 久久久精品免费免费高清| 精品一区二区三区四区五区乱码 | 国产精品久久久久久久久免| 中文天堂在线官网| 午夜老司机福利剧场| 在线天堂中文资源库| 久久久精品94久久精品| av黄色大香蕉| 99热国产这里只有精品6| 少妇人妻精品综合一区二区| 中文字幕人妻丝袜制服| 9色porny在线观看| 国产xxxxx性猛交| 日韩伦理黄色片| 在现免费观看毛片| 亚洲在久久综合| 美女国产视频在线观看| 啦啦啦啦在线视频资源| 日本av免费视频播放| 黑人欧美特级aaaaaa片| 精品国产一区二区三区四区第35| 成人毛片a级毛片在线播放| 女性被躁到高潮视频|