• <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ù)測控制
    他把我摸到了高潮在线观看| 亚洲全国av大片| 人妻久久中文字幕网| 男女之事视频高清在线观看| 此物有八面人人有两片| 一本综合久久免费| e午夜精品久久久久久久| 久久中文字幕一级| 神马国产精品三级电影在线观看 | 天天一区二区日本电影三级 | 亚洲五月色婷婷综合| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久久久久久大奶| 精品国产美女av久久久久小说| 精品不卡国产一区二区三区| 9色porny在线观看| 久9热在线精品视频| 成人亚洲精品一区在线观看| 中文字幕人成人乱码亚洲影| 麻豆成人av在线观看| xxx96com| 精品一区二区三区四区五区乱码| 丝袜美腿诱惑在线| 日本欧美视频一区| 午夜两性在线视频| 啦啦啦韩国在线观看视频| 久久久精品国产亚洲av高清涩受| 国产精品av久久久久免费| 99re在线观看精品视频| 日本黄色视频三级网站网址| 啦啦啦韩国在线观看视频| 欧美黄色淫秽网站| 亚洲精华国产精华精| 亚洲成人国产一区在线观看| 国产一级毛片七仙女欲春2 | 少妇的丰满在线观看| 久久天躁狠狠躁夜夜2o2o| 免费高清视频大片| 亚洲成人精品中文字幕电影| av免费在线观看网站| 我的亚洲天堂| 亚洲第一欧美日韩一区二区三区| 久久亚洲真实| 变态另类丝袜制服| 人人妻人人澡欧美一区二区 | 美女高潮到喷水免费观看| 午夜福利18| 老司机在亚洲福利影院| 国产精品一区二区免费欧美| 在线观看日韩欧美| 日本免费a在线| 精品人妻在线不人妻| 欧美日韩黄片免| 韩国av一区二区三区四区| 正在播放国产对白刺激| 国产精品秋霞免费鲁丝片| 99国产精品一区二区蜜桃av| 男女做爰动态图高潮gif福利片 | 巨乳人妻的诱惑在线观看| 国产一区二区三区视频了| 99国产精品一区二区三区| 人人澡人人妻人| 桃色一区二区三区在线观看| 色精品久久人妻99蜜桃| 黑丝袜美女国产一区| 亚洲片人在线观看| 纯流量卡能插随身wifi吗| 亚洲自偷自拍图片 自拍| 国产aⅴ精品一区二区三区波| 国产麻豆69| a级毛片在线看网站| 国产成人精品在线电影| 国产私拍福利视频在线观看| 久久人人精品亚洲av| 激情视频va一区二区三区| 久久精品国产综合久久久| 久久青草综合色| 欧美国产日韩亚洲一区| 操美女的视频在线观看| 国产成人一区二区三区免费视频网站| 麻豆av在线久日| 高清在线国产一区| 欧美日韩一级在线毛片| 一夜夜www| 免费高清视频大片| 久久天堂一区二区三区四区| 日韩免费av在线播放| 国产精品1区2区在线观看.| 波多野结衣巨乳人妻| 久久午夜综合久久蜜桃| 操美女的视频在线观看| 国产欧美日韩综合在线一区二区| 又大又爽又粗| 男女之事视频高清在线观看| 免费高清在线观看日韩| 国产精品二区激情视频| 99re在线观看精品视频| 国产高清视频在线播放一区| 女同久久另类99精品国产91| 国产熟女xx| 99精品在免费线老司机午夜| 久久久久久久久免费视频了| 十分钟在线观看高清视频www| 中文字幕人妻熟女乱码| 国产成人影院久久av| 亚洲五月天丁香| 久久人妻福利社区极品人妻图片| 国产三级黄色录像| 亚洲精品美女久久av网站| 99热只有精品国产| avwww免费| www.自偷自拍.com| 国产高清视频在线播放一区| 少妇粗大呻吟视频| 欧美成人午夜精品| av免费在线观看网站| 亚洲国产日韩欧美精品在线观看 | 精品少妇一区二区三区视频日本电影| 中文字幕精品免费在线观看视频| 国产成人免费无遮挡视频| 怎么达到女性高潮| 精品国产乱子伦一区二区三区| xxx96com| 午夜久久久久精精品| 国产av一区在线观看免费| 婷婷丁香在线五月| 亚洲免费av在线视频| x7x7x7水蜜桃| 免费不卡黄色视频| 国产麻豆成人av免费视频| 黄色片一级片一级黄色片| 亚洲成av人片免费观看| 成年人黄色毛片网站| 一级片免费观看大全| 成人免费观看视频高清| 国产一卡二卡三卡精品| 一级a爱片免费观看的视频| 人人妻人人澡欧美一区二区 | 国产精品 国内视频| 亚洲欧美一区二区三区黑人| av在线天堂中文字幕| 国产一区二区在线av高清观看| 国产私拍福利视频在线观看| 久久国产乱子伦精品免费另类| 又紧又爽又黄一区二区| 亚洲精品一卡2卡三卡4卡5卡| 精品卡一卡二卡四卡免费| 我的亚洲天堂| 午夜免费鲁丝| 脱女人内裤的视频| 啦啦啦韩国在线观看视频| 久久青草综合色| 欧美最黄视频在线播放免费| 久久国产精品人妻蜜桃| 黑人欧美特级aaaaaa片| 成人18禁高潮啪啪吃奶动态图| 国产av一区在线观看免费| svipshipincom国产片| www.精华液| 国产高清videossex| 国产激情久久老熟女| 大码成人一级视频| 中文字幕色久视频| 首页视频小说图片口味搜索| 日本vs欧美在线观看视频| 国产欧美日韩综合在线一区二区| 久久香蕉激情| 国产色视频综合| 亚洲五月婷婷丁香| 青草久久国产| 日韩高清综合在线| 99精品久久久久人妻精品| 亚洲人成伊人成综合网2020| 国产一卡二卡三卡精品| 成人特级黄色片久久久久久久| 欧美日韩瑟瑟在线播放| 97人妻天天添夜夜摸| 精品福利观看| 激情在线观看视频在线高清| 美女国产高潮福利片在线看| 国产精品乱码一区二三区的特点 | 在线视频色国产色| 欧美av亚洲av综合av国产av| 精品久久久久久,| 亚洲国产精品sss在线观看| 国产男靠女视频免费网站| 国产在线观看jvid| 国产蜜桃级精品一区二区三区| 1024香蕉在线观看| av视频免费观看在线观看| 精品乱码久久久久久99久播| 亚洲,欧美精品.| 国产精品一区二区三区四区久久 | 午夜日韩欧美国产| 窝窝影院91人妻| 久久国产亚洲av麻豆专区| 韩国av一区二区三区四区| 亚洲国产欧美网| 免费在线观看日本一区| 中文字幕av电影在线播放| 亚洲色图综合在线观看| 中文字幕av电影在线播放| 高潮久久久久久久久久久不卡| 丁香六月欧美| 亚洲av片天天在线观看| 亚洲av成人一区二区三| 丝袜美足系列| 黑人操中国人逼视频| 无限看片的www在线观看| 亚洲最大成人中文| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品一卡2卡三卡4卡5卡| 一级黄色大片毛片| 色婷婷久久久亚洲欧美| 又黄又爽又免费观看的视频| 欧美日韩瑟瑟在线播放| 一级黄色大片毛片| 精品国产亚洲在线| 亚洲三区欧美一区| 999久久久精品免费观看国产| 免费在线观看影片大全网站| 每晚都被弄得嗷嗷叫到高潮| 日本免费一区二区三区高清不卡 | 高清黄色对白视频在线免费看| 久久亚洲真实| 一进一出抽搐动态| 国内久久婷婷六月综合欲色啪| 亚洲精品av麻豆狂野| 免费无遮挡裸体视频| 不卡一级毛片| 欧美日本亚洲视频在线播放| 一级a爱片免费观看的视频| 黄片小视频在线播放| 成在线人永久免费视频| 一边摸一边抽搐一进一小说| 精品无人区乱码1区二区| 午夜福利成人在线免费观看| 日日爽夜夜爽网站| 欧美日韩精品网址| 国产精品综合久久久久久久免费 | 两个人免费观看高清视频| 变态另类丝袜制服| cao死你这个sao货| 操出白浆在线播放| 国内久久婷婷六月综合欲色啪| 男男h啪啪无遮挡| 曰老女人黄片| 纯流量卡能插随身wifi吗| 中文字幕人成人乱码亚洲影| 婷婷六月久久综合丁香| 美女高潮到喷水免费观看| 久久亚洲真实| 操美女的视频在线观看| 法律面前人人平等表现在哪些方面| 波多野结衣巨乳人妻| 国产精品免费视频内射| 此物有八面人人有两片| 人人妻人人爽人人添夜夜欢视频| 色综合站精品国产| 久久久久久久精品吃奶| 制服诱惑二区| 成年女人毛片免费观看观看9| 亚洲在线自拍视频| 色播在线永久视频| 色老头精品视频在线观看| 丝袜在线中文字幕| 美女大奶头视频| 岛国视频午夜一区免费看| 午夜亚洲福利在线播放| 久9热在线精品视频| 美女扒开内裤让男人捅视频| 日日爽夜夜爽网站| 手机成人av网站| 天天躁夜夜躁狠狠躁躁| 亚洲,欧美精品.| cao死你这个sao货| 久久久久国内视频| 如日韩欧美国产精品一区二区三区| 两人在一起打扑克的视频| 亚洲在线自拍视频| 老司机深夜福利视频在线观看| 亚洲全国av大片| av欧美777| 亚洲av电影在线进入| 久久久久九九精品影院| 日本撒尿小便嘘嘘汇集6| 91大片在线观看| 搡老熟女国产l中国老女人| 九色亚洲精品在线播放| 激情在线观看视频在线高清| 色哟哟哟哟哟哟| avwww免费| 黑丝袜美女国产一区| 精品欧美国产一区二区三| 成人亚洲精品av一区二区| 国产亚洲精品一区二区www| 99久久久亚洲精品蜜臀av| 国产一区在线观看成人免费| 国产黄a三级三级三级人| 18美女黄网站色大片免费观看| 首页视频小说图片口味搜索| 一级毛片女人18水好多| 亚洲人成网站在线播放欧美日韩| 亚洲精品一区av在线观看| 亚洲精品美女久久久久99蜜臀| 久久人妻福利社区极品人妻图片| 在线观看免费午夜福利视频| 日本在线视频免费播放| 国产免费男女视频| 18禁观看日本| 亚洲,欧美精品.| 男女下面插进去视频免费观看| 亚洲人成电影免费在线| 在线播放国产精品三级| 中文字幕精品免费在线观看视频| 日韩免费av在线播放| 国产成年人精品一区二区| 色精品久久人妻99蜜桃| 亚洲第一欧美日韩一区二区三区| 亚洲av片天天在线观看| 精品电影一区二区在线| 亚洲av美国av| 久久天躁狠狠躁夜夜2o2o| 老司机午夜十八禁免费视频| 最近最新中文字幕大全电影3 | 精品国产乱码久久久久久男人| 男人的好看免费观看在线视频 | 不卡一级毛片| 久久欧美精品欧美久久欧美| 久久久久久人人人人人| 在线永久观看黄色视频| 国产高清有码在线观看视频 | 亚洲成av片中文字幕在线观看| 视频区欧美日本亚洲| 久久久久九九精品影院| 成人18禁高潮啪啪吃奶动态图| 女性被躁到高潮视频| 日本五十路高清| 亚洲熟妇中文字幕五十中出| 久久天躁狠狠躁夜夜2o2o| 国产97色在线日韩免费| 高潮久久久久久久久久久不卡| 午夜两性在线视频| avwww免费| 久久精品人人爽人人爽视色| 久久九九热精品免费| 欧美另类亚洲清纯唯美| 91在线观看av| 亚洲 国产 在线| АⅤ资源中文在线天堂| 一本久久中文字幕| 欧美丝袜亚洲另类 | 亚洲成人免费电影在线观看| 国产av在哪里看| 十八禁网站免费在线| 一本久久中文字幕| 久久草成人影院| 黑人操中国人逼视频| 久久久久久免费高清国产稀缺| 亚洲片人在线观看| av网站免费在线观看视频| 欧美一级毛片孕妇| 久久精品亚洲精品国产色婷小说| 亚洲av成人不卡在线观看播放网| 黄色视频不卡| 午夜精品在线福利| 中文字幕最新亚洲高清| 午夜两性在线视频| 91麻豆精品激情在线观看国产| av网站免费在线观看视频| 久久久久久亚洲精品国产蜜桃av| 亚洲国产日韩欧美精品在线观看 | 亚洲av美国av| 一夜夜www| 十分钟在线观看高清视频www| 一级毛片高清免费大全| 日本免费一区二区三区高清不卡 | 国产成人欧美在线观看| 在线国产一区二区在线| 亚洲全国av大片| 成人亚洲精品一区在线观看| 少妇的丰满在线观看| 国产精品国产高清国产av| 搡老妇女老女人老熟妇| 两人在一起打扑克的视频| 亚洲一区中文字幕在线| 亚洲国产看品久久| 免费观看人在逋| 精品欧美一区二区三区在线| 搡老熟女国产l中国老女人| 一级a爱片免费观看的视频| 精品久久久久久久毛片微露脸| 日韩有码中文字幕| 真人做人爱边吃奶动态| 自线自在国产av| 极品教师在线免费播放| aaaaa片日本免费| 亚洲av成人一区二区三| 国产男靠女视频免费网站| 免费观看人在逋| 亚洲第一欧美日韩一区二区三区| 色婷婷久久久亚洲欧美| 日本 av在线| 中亚洲国语对白在线视频| 久久国产精品人妻蜜桃| 欧美一区二区精品小视频在线| 变态另类丝袜制服| 亚洲欧美激情综合另类| 亚洲专区字幕在线| 国产成人系列免费观看| 国产伦人伦偷精品视频| 91麻豆av在线| 亚洲国产精品合色在线| 国产免费av片在线观看野外av| www.www免费av| 夜夜躁狠狠躁天天躁| 91字幕亚洲| 性色av乱码一区二区三区2| 熟女少妇亚洲综合色aaa.| 国产精品一区二区精品视频观看| 国产91精品成人一区二区三区| 精品熟女少妇八av免费久了| 中文字幕av电影在线播放| ponron亚洲| 国产一区二区激情短视频| 黄色片一级片一级黄色片| 欧美 亚洲 国产 日韩一| 午夜a级毛片| 99riav亚洲国产免费| xxx96com| 精品电影一区二区在线| 黑丝袜美女国产一区| 人妻久久中文字幕网| 女同久久另类99精品国产91| 日本精品一区二区三区蜜桃| 国产欧美日韩综合在线一区二区| 精品久久久久久久人妻蜜臀av | 亚洲精品中文字幕在线视频| 久久久国产欧美日韩av| 两个人视频免费观看高清| 欧美中文综合在线视频| 91麻豆av在线| 一区二区三区国产精品乱码| 女人高潮潮喷娇喘18禁视频| 51午夜福利影视在线观看| 国产日韩一区二区三区精品不卡| 丝袜在线中文字幕| 天天添夜夜摸| 精品久久久久久成人av| 女同久久另类99精品国产91| 黄色视频不卡| 中出人妻视频一区二区| 一二三四在线观看免费中文在| 亚洲精品美女久久久久99蜜臀| 亚洲免费av在线视频| 国产精品99久久99久久久不卡| 深夜精品福利| av片东京热男人的天堂| 欧美日韩瑟瑟在线播放| av视频免费观看在线观看| 亚洲成av人片免费观看| 午夜福利影视在线免费观看| 黑人巨大精品欧美一区二区mp4| 亚洲精品久久成人aⅴ小说| 亚洲人成电影观看| 亚洲国产欧美网| 国产av一区二区精品久久| 国产极品粉嫩免费观看在线| 搞女人的毛片| 亚洲av美国av| 女人被狂操c到高潮| 女同久久另类99精品国产91| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看一区二区三区| 好男人在线观看高清免费视频 | 亚洲 欧美一区二区三区| 波多野结衣av一区二区av| 免费看十八禁软件| 怎么达到女性高潮| 亚洲国产毛片av蜜桃av| 嫁个100分男人电影在线观看| 91精品国产国语对白视频| 国产不卡一卡二| 久9热在线精品视频| 亚洲专区字幕在线| 欧美老熟妇乱子伦牲交| 老司机午夜十八禁免费视频| 亚洲av五月六月丁香网| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲人成电影观看| 久久精品国产综合久久久| 丰满的人妻完整版| x7x7x7水蜜桃| 一级黄色大片毛片| 日韩视频一区二区在线观看| 久久国产亚洲av麻豆专区| 91精品三级在线观看| 国内久久婷婷六月综合欲色啪| 亚洲国产精品sss在线观看| 极品人妻少妇av视频| 一边摸一边做爽爽视频免费| av在线天堂中文字幕| 欧美久久黑人一区二区| 婷婷六月久久综合丁香| 国产国语露脸激情在线看| 村上凉子中文字幕在线| 久久九九热精品免费| 国产精品98久久久久久宅男小说| 极品教师在线免费播放| 免费在线观看日本一区| 男女做爰动态图高潮gif福利片 | 非洲黑人性xxxx精品又粗又长| 又黄又爽又免费观看的视频| 国产精品九九99| 国产精品久久电影中文字幕| 自线自在国产av| 99riav亚洲国产免费| 欧美激情久久久久久爽电影 | 色尼玛亚洲综合影院| 黑人巨大精品欧美一区二区蜜桃| a在线观看视频网站| 亚洲国产日韩欧美精品在线观看 | 久久这里只有精品19| 亚洲成人免费电影在线观看| 久久精品人人爽人人爽视色| 亚洲精品国产色婷婷电影| 人人妻,人人澡人人爽秒播| 男人的好看免费观看在线视频 | 国产精品九九99| 精品第一国产精品| 中文字幕人成人乱码亚洲影| 国产精品美女特级片免费视频播放器 | 禁无遮挡网站| 无限看片的www在线观看| 欧美黑人欧美精品刺激| 久久草成人影院| 麻豆久久精品国产亚洲av| 免费人成视频x8x8入口观看| 国产熟女xx| 亚洲精品美女久久av网站| 丰满人妻熟妇乱又伦精品不卡| 在线观看66精品国产| 精品国产乱码久久久久久男人| 麻豆一二三区av精品| 亚洲av美国av| 欧洲精品卡2卡3卡4卡5卡区| 久久久久亚洲av毛片大全| 日韩有码中文字幕| 一边摸一边抽搐一进一出视频| 久久久国产精品麻豆| av天堂久久9| av视频免费观看在线观看| 婷婷丁香在线五月| 国产在线观看jvid| av片东京热男人的天堂| 欧美亚洲日本最大视频资源| 国产成人精品久久二区二区91| 中文字幕av电影在线播放| 1024香蕉在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 午夜影院日韩av| 免费不卡黄色视频| 亚洲av成人不卡在线观看播放网| 精品无人区乱码1区二区| 国产亚洲精品第一综合不卡| 美女 人体艺术 gogo| 日韩欧美国产一区二区入口| 国产亚洲欧美98| 中文字幕人成人乱码亚洲影| 老司机午夜福利在线观看视频| 欧美av亚洲av综合av国产av| 亚洲免费av在线视频| 波多野结衣av一区二区av| 桃色一区二区三区在线观看| 午夜影院日韩av| 国产日韩一区二区三区精品不卡| 国产不卡一卡二| 欧美黄色淫秽网站| 欧美+亚洲+日韩+国产| 中文字幕人妻熟女乱码| 久久香蕉精品热| 啪啪无遮挡十八禁网站| 色婷婷久久久亚洲欧美| 中文字幕另类日韩欧美亚洲嫩草| 麻豆一二三区av精品| 欧美激情极品国产一区二区三区| 国产精品一区二区免费欧美| 男女做爰动态图高潮gif福利片 | www.精华液| 欧美黑人精品巨大| 成熟少妇高潮喷水视频| 午夜影院日韩av| 精品无人区乱码1区二区| 久久人人97超碰香蕉20202| 99在线视频只有这里精品首页| 亚洲国产高清在线一区二区三 | 国产精品久久久久久亚洲av鲁大| 久久午夜亚洲精品久久| 亚洲熟女毛片儿| 久久久久九九精品影院| 国产精品98久久久久久宅男小说| 日日干狠狠操夜夜爽| 757午夜福利合集在线观看| 叶爱在线成人免费视频播放| 亚洲中文日韩欧美视频| 国产精品九九99| 麻豆一二三区av精品| 国产精品九九99| 99精品在免费线老司机午夜| 在线天堂中文资源库| 日韩大尺度精品在线看网址 | 亚洲欧美精品综合久久99| 涩涩av久久男人的天堂| av片东京热男人的天堂| 亚洲成人精品中文字幕电影|