两个人的电影免费视频_国产精品久久久久久久久成人_97视频在线观看播放_久久这里只有精品777_亚洲熟女少妇二三区_4438x8成人网亚洲av_内谢国产内射夫妻免费视频_人妻精品久久久久中国字幕

考慮多風(fēng)場預(yù)測誤差時(shí)空相關(guān)性的儲(chǔ)能容量優(yōu)化配置方法與流程

文檔序號(hào):11233652閱讀:660來源:國知局
考慮多風(fēng)場預(yù)測誤差時(shí)空相關(guān)性的儲(chǔ)能容量優(yōu)化配置方法與流程

本發(fā)明涉及一種考慮多風(fēng)場預(yù)測誤差時(shí)空相關(guān)性的儲(chǔ)能容量優(yōu)化配置方法。



背景技術(shù):

由于資源與環(huán)境問題日益嚴(yán)峻,以風(fēng)能為首的新能源開發(fā)和利用在全球范圍內(nèi)受到廣泛關(guān)注,并持續(xù)保持快速發(fā)展水平。根據(jù)全球風(fēng)能理事會(huì)(gwec)發(fā)布的全球風(fēng)能數(shù)據(jù),截止到2015年,全球風(fēng)電裝機(jī)容量達(dá)到433gw且2016年又將增加60gw。我國風(fēng)電裝機(jī)容量截止2015年達(dá)到145.4gw,比2014年增長了21.2%。然而,由于風(fēng)電機(jī)組出力的不確定性和波動(dòng)性,風(fēng)電大規(guī)模并網(wǎng)必將對(duì)電力系統(tǒng)安全、可靠、穩(wěn)定和經(jīng)濟(jì)運(yùn)行帶來巨大挑戰(zhàn)。

隨著儲(chǔ)能技術(shù)的發(fā)展,安裝儲(chǔ)能系統(tǒng)(ess)被看作減小風(fēng)電機(jī)組隨機(jī)性和不確定性影響的有效方法。然而,由于目前多數(shù)類型的儲(chǔ)能價(jià)格依然相對(duì)較為昂貴,優(yōu)化儲(chǔ)能容量是當(dāng)前儲(chǔ)能應(yīng)用中的一個(gè)關(guān)鍵性問題。目前,關(guān)于風(fēng)電場儲(chǔ)能的優(yōu)化配置問題,學(xué)者們進(jìn)行了大量的有益的研究工作。

但是,現(xiàn)有的研究僅針對(duì)單一風(fēng)場,并未考慮多風(fēng)場之間的相關(guān)關(guān)系。而距離較近的風(fēng)電場之間的風(fēng)電功率預(yù)測誤差具有非常強(qiáng)的相關(guān)性,這一特性將直接影響風(fēng)電功率的總體波動(dòng)特性、儲(chǔ)能配置結(jié)果的經(jīng)濟(jì)性與有效性。

預(yù)測誤差的時(shí)空相關(guān)性分析時(shí)風(fēng)儲(chǔ)系統(tǒng)研究的重要基礎(chǔ),準(zhǔn)確分析風(fēng)場間時(shí)空相關(guān)性對(duì)于確定合理的儲(chǔ)能容量十分重要,因此有必要對(duì)考慮預(yù)測誤差特性和時(shí)空相關(guān)性的儲(chǔ)能優(yōu)化配置進(jìn)行進(jìn)一步研究,尤其是現(xiàn)有的考慮誤差分布的儲(chǔ)能容量優(yōu)化研究存在以下問題:

(a)并未同時(shí)考慮風(fēng)場之間的時(shí)空相關(guān)性;

(b)利用線性方法擬合預(yù)測誤差無法保持?jǐn)?shù)據(jù)的原始分布特性,用于規(guī)劃儲(chǔ)能會(huì)錯(cuò)誤估計(jì)所需容量;

(c)均未提到預(yù)測誤差時(shí)間相關(guān)性和空間相關(guān)性如何影響風(fēng)儲(chǔ)系統(tǒng)容量配置和成本效益。



技術(shù)實(shí)現(xiàn)要素:

本發(fā)明為了解決上述問題,提出了一種考慮多風(fēng)場預(yù)測誤差時(shí)空相關(guān)性的儲(chǔ)能容量優(yōu)化配置方法,本發(fā)明基于多場景分析方法,以儲(chǔ)能系統(tǒng)綜合成本效益最優(yōu)為目標(biāo),建立考慮多風(fēng)場間風(fēng)電功率預(yù)測誤差時(shí)空相關(guān)性的儲(chǔ)能容量優(yōu)化模型。

為了實(shí)現(xiàn)上述目的,本發(fā)明采用如下技術(shù)方案:

一種考慮多風(fēng)場預(yù)測誤差時(shí)空相關(guān)性的儲(chǔ)能容量優(yōu)化配置方法,通過對(duì)比分析參數(shù)和非參數(shù)兩種方法對(duì)預(yù)測誤差的擬合效果,選取擬合精度最高的kde方法進(jìn)行誤差擬合;采用copula理論的多風(fēng)電場預(yù)測誤差時(shí)空相關(guān)性建模方法,根據(jù)多風(fēng)場間實(shí)際空間相關(guān)關(guān)系,得到多風(fēng)場預(yù)測誤差聯(lián)合分布函數(shù),利用kde方法擬合預(yù)測誤差的邊緣分布函數(shù);基于多場景分析方法,以儲(chǔ)能系統(tǒng)投資成本和運(yùn)行成本最小為目標(biāo),建立考慮多風(fēng)場間預(yù)測誤差時(shí)空相關(guān)性的儲(chǔ)能容量優(yōu)化模型。

進(jìn)一步的,所述預(yù)測誤差為功率實(shí)際值減去預(yù)測值,且位于[-1,1]之間的數(shù)據(jù)。

進(jìn)一步的,分別利用正態(tài)分布、t分布和kde方法擬合預(yù)測誤差來進(jìn)行風(fēng)功率預(yù)測誤差分析,確定最優(yōu)擬合方法。

進(jìn)一步的,進(jìn)行多風(fēng)電場預(yù)測誤差時(shí)空相關(guān)性分析,風(fēng)電場之間的距離越近,其預(yù)測誤差線性相關(guān)系數(shù)越大,時(shí)空相關(guān)性越強(qiáng)。

進(jìn)一步的,采用copula理論的多風(fēng)電場預(yù)測誤差時(shí)空相關(guān)性建模的具體過程包括:

(1)生成n風(fēng)場功率預(yù)測誤差矩陣,假設(shè)每風(fēng)場有t個(gè)時(shí)間段,每個(gè)時(shí)間段內(nèi)有t個(gè)觀測值,構(gòu)建預(yù)測誤差矩陣;

(2)利用kde擬合方法計(jì)算預(yù)測誤差邊緣分布;

(3)利用最大似然估計(jì)法對(duì)normal/tcopula函數(shù)參數(shù)進(jìn)行估計(jì),其中normalcopula函數(shù)參數(shù)為t×t秩相關(guān)矩陣ρ,tcopula函數(shù)參數(shù)為t×t秩相關(guān)矩陣和自由度;

(4)利用normal/tcopula方程計(jì)算得到多風(fēng)場預(yù)測誤差聯(lián)合分布函數(shù)。

進(jìn)一步的,根據(jù)構(gòu)建的normal/tcopula生成的聯(lián)合分布函數(shù)進(jìn)行反變換,生成預(yù)測誤差場景。

進(jìn)一步的,利用向后縮減技術(shù)對(duì)原始場景進(jìn)行縮減,尋找最接近原始系統(tǒng)的場景,計(jì)算確定的每個(gè)場景的棄風(fēng)能量和損失能量。

進(jìn)一步的,蓄電池儲(chǔ)能容量優(yōu)化的目標(biāo)為補(bǔ)償風(fēng)電功率預(yù)測誤差和使儲(chǔ)能系統(tǒng)投資成本與風(fēng)儲(chǔ)系統(tǒng)運(yùn)行成本最小以實(shí)現(xiàn)綜合效益最優(yōu)。

進(jìn)一步的,利用折中方法考慮投資成本和運(yùn)行成本間的對(duì)立關(guān)系以實(shí)現(xiàn)風(fēng)儲(chǔ)系統(tǒng)的優(yōu)化運(yùn)行,運(yùn)行成本包括棄風(fēng)成本和損失成本,兩者隨額定容量和充放電功率的變化而變化。

進(jìn)一步的,蓄電池儲(chǔ)能容量優(yōu)化的約束條件包括蓄電池約束和輸出功率約束,蓄電池約束又包括充電功率約束和放電功率約束。

與現(xiàn)有技術(shù)相比,本發(fā)明的有益效果為:

(1)本發(fā)明利用kde擬合各時(shí)間段預(yù)測誤差的概率分布函數(shù)和多元聯(lián)合分布函數(shù)的邊緣分布函數(shù),擬合效果最好;

(2)本發(fā)明建立基于copula函數(shù)的預(yù)測誤差時(shí)空相關(guān)性模型,得到多風(fēng)場預(yù)測誤差的聯(lián)合分布函數(shù),該聯(lián)合分布包含了各風(fēng)場的時(shí)間相關(guān)性和風(fēng)場間的空間相關(guān)性;

(3)本發(fā)明構(gòu)建了考慮風(fēng)電場時(shí)空相關(guān)性的儲(chǔ)能容量規(guī)劃模型。該模型也考慮了投資成本與運(yùn)行成本間的對(duì)立關(guān)系,模型目標(biāo)是實(shí)現(xiàn)風(fēng)儲(chǔ)系統(tǒng)的最優(yōu)綜合效益;

(4)算例計(jì)算結(jié)果證明了對(duì)bess容量進(jìn)行優(yōu)化可以有效補(bǔ)償風(fēng)電功率預(yù)測誤差以及減小預(yù)測方法局限的影響;通過對(duì)計(jì)算結(jié)果進(jìn)行對(duì)比分析,預(yù)測誤差時(shí)空相關(guān)性對(duì)于bess配置有顯著影響,風(fēng)電場間的空間相關(guān)性越強(qiáng),風(fēng)儲(chǔ)系統(tǒng)綜合成本效益越好,對(duì)風(fēng)電功率預(yù)測誤差的補(bǔ)償效果越好。

附圖說明

構(gòu)成本申請(qǐng)的一部分的說明書附圖用來提供對(duì)本申請(qǐng)的進(jìn)一步理解,本申請(qǐng)的示意性實(shí)施例及其說明用于解釋本申請(qǐng),并不構(gòu)成對(duì)本申請(qǐng)的不當(dāng)限定。

圖1是本發(fā)明實(shí)施例中四個(gè)風(fēng)電場地理位置分布圖;

圖2(a)-圖2(d)是本發(fā)明a風(fēng)場24小時(shí)預(yù)測誤差均值、方差、偏度與峰度計(jì)算結(jié)果圖;

圖3(a)-圖3(d)是本發(fā)明a風(fēng)場分別提前1小時(shí)、提前7小時(shí)、提前8小時(shí)和提前24小時(shí)預(yù)測誤差概率密度分布直方圖和利用normal分布,t分布以及kde方法得到的擬合概率密度函數(shù);

圖4(a)-圖4(d)是本發(fā)明a風(fēng)場分別提前1小時(shí)、提前7小時(shí)、提前8小時(shí)和提前24小時(shí)預(yù)測誤差擬合分布函數(shù)與經(jīng)驗(yàn)分布函數(shù)qq圖;

圖5是本發(fā)明預(yù)測誤差的線性相關(guān)結(jié)構(gòu)圖;

圖6是本發(fā)明場景生成步驟圖;

圖7是本發(fā)明a風(fēng)場第一天功率觀測值、預(yù)測值和置信區(qū)間示意圖;

圖8(a)-圖8(d)是本發(fā)明a風(fēng)場第一天預(yù)測誤差場景縮減結(jié)果,其中,(a)不考慮時(shí)空相關(guān)性;(b)時(shí)空相關(guān)性最強(qiáng);(c)時(shí)間相關(guān)性最強(qiáng)空間相關(guān)性較強(qiáng),(d)時(shí)間相關(guān)性最強(qiáng)空間相關(guān)性最弱;

圖9是本發(fā)明a風(fēng)場原始預(yù)測誤差曲線和經(jīng)過不同時(shí)空相關(guān)性的儲(chǔ)能補(bǔ)償后的預(yù)測誤差曲線圖。

具體實(shí)施方式:

下面結(jié)合附圖與實(shí)施例對(duì)本發(fā)明作進(jìn)一步說明。

應(yīng)該指出,以下詳細(xì)說明都是例示性的,旨在對(duì)本申請(qǐng)?zhí)峁┻M(jìn)一步的說明。除非另有指明,本發(fā)明使用的所有技術(shù)和科學(xué)術(shù)語具有與本申請(qǐng)所屬技術(shù)領(lǐng)域的普通技術(shù)人員通常理解的相同含義。

需要注意的是,這里所使用的術(shù)語僅是為了描述具體實(shí)施方式,而非意圖限制根據(jù)本申請(qǐng)的示例性實(shí)施方式。如在這里所使用的,除非上下文另外明確指出,否則單數(shù)形式也意圖包括復(fù)數(shù)形式,此外,還應(yīng)當(dāng)理解的是,當(dāng)在本說明書中使用術(shù)語“包含”和/或“包括”時(shí),其指明存在特征、步驟、操作、器件、組件和/或它們的組合。

正如背景技術(shù)所介紹的,現(xiàn)有技術(shù)中存在并未同時(shí)考慮風(fēng)場之間的時(shí)空相關(guān)性以及利用線性方法擬合預(yù)測誤差無法保持?jǐn)?shù)據(jù)的原始分布特性,用于規(guī)劃儲(chǔ)能會(huì)錯(cuò)誤估計(jì)所需容量的不足,為了解決如上的技術(shù)問題,本申請(qǐng)?zhí)岢隽艘环N考慮多風(fēng)場預(yù)測誤差時(shí)空相關(guān)性的儲(chǔ)能容量優(yōu)化配置方法。

主要貢獻(xiàn)如下:

1)提出了基于copula函數(shù)的預(yù)測誤差高維建模方法,建立多風(fēng)場預(yù)測誤差的聯(lián)合分布函數(shù),該聯(lián)合分布包含了各風(fēng)場的時(shí)間相關(guān)性和風(fēng)場間的空間相關(guān)性。

2)基于預(yù)測誤差的聯(lián)合分布,本發(fā)明考慮儲(chǔ)能系統(tǒng)的投資成本與運(yùn)行成本兩者間對(duì)立關(guān)系,基于多場景技術(shù),對(duì)儲(chǔ)能容量進(jìn)行優(yōu)化配置,使風(fēng)儲(chǔ)系統(tǒng)的綜合效益達(dá)到最優(yōu)。

3)通過對(duì)比分析預(yù)測誤差時(shí)空相關(guān)性對(duì)風(fēng)電場儲(chǔ)能配置有顯著影響,風(fēng)電場間的空間相關(guān)性越強(qiáng),風(fēng)儲(chǔ)系統(tǒng)綜合成本效益越好,對(duì)風(fēng)電功率預(yù)測誤差的補(bǔ)償效果越好。

本申請(qǐng)的一種典型的實(shí)施方式中,所用數(shù)據(jù)來自冀北地區(qū)的四個(gè)風(fēng)電場,其地理位置分布如圖1所示,四個(gè)數(shù)據(jù)集合分別為a、b、c和d。每個(gè)數(shù)據(jù)集包含2013年全年數(shù)據(jù),并根據(jù)各自裝機(jī)容量進(jìn)行歸一化處理。本發(fā)明采用文獻(xiàn)liangz,liangj,wangc,etal.short-termwindpowercombinedforecastingbasedonerrorforecastcorrection[j].energyconversionandmanagement,2016,119:215-226提出的風(fēng)電功率方法預(yù)測風(fēng)功率得到8688個(gè)預(yù)測功率數(shù)據(jù)和相應(yīng)的預(yù)測誤差。預(yù)測誤差等于功率實(shí)際值減去預(yù)測值,且位于[-1,1]之間。

圖2為a風(fēng)場24小時(shí)預(yù)測誤差的均值、方差、偏度與峰度計(jì)算結(jié)果,bcd三風(fēng)場的計(jì)算結(jié)果圖如附錄圖1—圖3所示。由圖2(a)-圖2(d)可知,圖2(a)各時(shí)段預(yù)測誤差的均值均不為零;圖2(b)隨著預(yù)測區(qū)間增大,預(yù)測誤差的方差逐漸減?。粓D2(c)預(yù)測誤差大于2h后的偏度系數(shù)均小于零,表示對(duì)應(yīng)時(shí)間段的預(yù)測誤差分布為左偏,且隨著預(yù)測區(qū)間的增大,偏度系數(shù)絕對(duì)值逐漸增大,表示拖尾越長;圖2(d)峰度隨著預(yù)測區(qū)間的增大而減小,表示誤差分布形態(tài)逐漸變緩,且峰度系數(shù)均大于3,表示誤差分布有明顯的厚尾特征。因此,a風(fēng)場24小時(shí)預(yù)測誤差形態(tài)呈現(xiàn)尖峰、厚尾、拖尾的特征,為了更加清晰地分析預(yù)測誤差的特征,本發(fā)明利用不同擬合方法計(jì)算并繪制預(yù)測誤差的概率分布,如圖3(a)-圖3(d)所示,為a多風(fēng)場提前1小時(shí)、7小時(shí)、8小時(shí)和24小時(shí)的預(yù)測誤差概率密度分布直方圖和利用正態(tài)分布、t分布和kde方法擬合的概率密度函數(shù)對(duì)比圖,可知,隨著預(yù)測區(qū)間的增大,(a)預(yù)測誤差概率密度分布曲線逐漸由凸變平;(b)基于正態(tài)分布的概率密度函數(shù)與預(yù)測誤差的經(jīng)驗(yàn)概率密度函數(shù)之間的區(qū)別逐漸增大,而基于kde的概率密度分布與經(jīng)驗(yàn)概率密度函數(shù)之間無明顯差別;(e)當(dāng)預(yù)測區(qū)間大于7小時(shí)時(shí),t分布由于最大似然估計(jì)發(fā)散而無法擬合預(yù)測誤差;(f)預(yù)測誤差的概率密度分布的尾部逐漸變厚。

qq圖是用來擬合分布函數(shù)與經(jīng)驗(yàn)分布函數(shù)之間的距離,若兩者都重合或者接近于45°直線,說明擬合度較高。圖4(a)-圖4(d)為a風(fēng)場提前1小時(shí)、7小時(shí)、8小時(shí)和24小時(shí)預(yù)測誤差擬合分布函數(shù)與經(jīng)驗(yàn)分布函數(shù)qq圖,可以看出kde擬合分布函數(shù)最接近45°直線,擬合效果最好?;谝陨戏治觯景l(fā)明選用kde方法擬合copula函數(shù)的邊緣分布函數(shù)。

預(yù)測誤差線性相關(guān)結(jié)構(gòu)如圖5所示。圖5表示相鄰時(shí)間區(qū)間的預(yù)測誤差具有很強(qiáng)的時(shí)間相關(guān)性。結(jié)合圖1的四風(fēng)場地理位置分布,ab兩風(fēng)場距離最小,其最大線性相關(guān)系數(shù)為0.6425,因此ab兩風(fēng)場預(yù)測誤差時(shí)空相關(guān)性最強(qiáng)。ac兩風(fēng)場的最大線性相關(guān)系數(shù)為0.4338,其距離大于ab間距離。ad兩風(fēng)場相距最遠(yuǎn),其最大線性相關(guān)系數(shù)為0.06536,幾乎沒有相關(guān)性。因此,風(fēng)電場距離越近,其預(yù)測誤差線性相關(guān)系數(shù)越大,時(shí)空相關(guān)性越強(qiáng)。

copula函數(shù)提供了一種構(gòu)造不同獨(dú)立變量的聯(lián)合分布的方法,該分布包含了隨機(jī)變量的相關(guān)關(guān)系并且能夠捕捉其非線性相關(guān)關(guān)系。

數(shù)學(xué)上,sklar定理強(qiáng)調(diào)任何多元聯(lián)合分布可以寫成單變量的邊緣分布和一個(gè)描述變量間相關(guān)結(jié)構(gòu)的copula函數(shù)。假設(shè)x1,x2,…xn為隨機(jī)變量,其邊緣分布函數(shù)為f1(x1),f2(x2),…,ft(xt),且均遵循以下均勻分布:

f1(x1)~u(0,1),f2(x2)~u(0,1),…,fn(xn)~u(0,1)(2)

根據(jù)sklar定理,f1(x1),f2(x2),…,ft(xt)的聯(lián)合分布為:

f(x1,x2,…,xt)=c(f1(x1),f2(x2),…,ft(xt))(3)

其中c即為copula函數(shù)。

由于normal/tcopula可以用于多元建模,允許任何邊緣分布和正定矩陣,因此本發(fā)明利用normal/tcopula函數(shù)建立日前預(yù)測誤差多元分布函數(shù)。normal/tcopula函數(shù)如式(4)(5)所示:

f(x1,x2,…,xt;ρ)=c(f1(x1),f2(x2),…,ft(xt);ρ)

=φρ(φ-1(f1(x1)),φ-1(f2(x2)),…,φ-1(ft(xt)))(4)

f(x1,x2,…,xt;ρ,k)=c(f1(x1),f2(x2),…,ft(xt);ρ,k)

=tρ,k(tk-1(f1(x1)),tk-1(f2(x2)),…,tk-1(ft(xt)))(5)

其中,

φ-1標(biāo)準(zhǔn)正態(tài)分布的分布函數(shù)的反函數(shù)

φρ標(biāo)準(zhǔn)正態(tài)分布的分布函數(shù)

ρ對(duì)稱正定矩陣

t分布的分布函數(shù)的反函數(shù)

tρ,kt分布的分布函數(shù)

k自由度

利用copula理論的多風(fēng)場風(fēng)電功率預(yù)測誤差建模步驟

基于以上分析,n風(fēng)場預(yù)測誤差的高維建模步驟如下。

1)生成n風(fēng)場功率預(yù)測誤差矩陣。假設(shè)每風(fēng)場有t個(gè)時(shí)間段,每個(gè)時(shí)間段內(nèi)有t個(gè)觀測值。則預(yù)測誤差矩陣可由式(6)表示。

其中為第n個(gè)風(fēng)電場t時(shí)間段第t個(gè)預(yù)測誤差值。

2)利用kde擬合方法計(jì)算預(yù)測誤差邊緣分布。擬合結(jié)果形式如式(7)所示。

3)利用最大似然估計(jì)法對(duì)normal/tcopula函數(shù)參數(shù)進(jìn)行估計(jì)。normalcopula函數(shù)參數(shù)為t×t秩相關(guān)矩陣ρ,tcopula函數(shù)參數(shù)為t×t秩相關(guān)矩陣ρ和自由度k。

4)利用normal/tcopula方程計(jì)算得到多風(fēng)場預(yù)測誤差聯(lián)合分布函數(shù)(jcdf)。

在由normal/tcopula生成的聯(lián)合分布函數(shù)的基礎(chǔ)上,可生成預(yù)測誤差場景。圖6為多場景的生成步驟,其中ppre為風(fēng)電功率預(yù)測值。根據(jù)圖6可生成m-t維風(fēng)電功率預(yù)測誤差場景。

為減小模型規(guī)模提高運(yùn)算效率,有必要對(duì)中生成的場景進(jìn)行縮減。本發(fā)明利用向后縮減技術(shù)對(duì)原始場景進(jìn)行縮減。以尋找最接近原始系統(tǒng)的場景。

蓄電池儲(chǔ)能容量優(yōu)化的目的為:(a)補(bǔ)償風(fēng)電功率預(yù)測誤差;(b)使儲(chǔ)能系統(tǒng)投資成本與風(fēng)儲(chǔ)系統(tǒng)運(yùn)行成本最小以實(shí)現(xiàn)綜合效益最優(yōu)。本發(fā)明利用折中方法考慮投資成本和運(yùn)行成本間的對(duì)立關(guān)系以實(shí)現(xiàn)風(fēng)儲(chǔ)系統(tǒng)的優(yōu)化運(yùn)行。運(yùn)行成本包括棄風(fēng)成本和損失成本,兩者隨額定容量和充放電功率的變化而變化。為更加有效的分析風(fēng)電功率預(yù)測誤差的時(shí)空相關(guān)性對(duì)蓄電池儲(chǔ)能容量優(yōu)化配置的影響,本發(fā)明分析當(dāng)?shù)氐娘L(fēng)況并選取特征最為明顯的月份作為研究對(duì)象。

縮減場景數(shù)之后,每個(gè)場景的棄風(fēng)能量和損失能量如式(8)(9)所示。

式(8)與(9)中

ci,bessn=max(ci,bess(t))-min(ci,bess(t))(15)

pi,bessn=max(pi,besscen,|pi,bessden|)(16)

對(duì)于第i個(gè)場景,其中

si,su1(t)布爾量,用于描述由充電功率限制產(chǎn)生的棄風(fēng)能量

si,su2(t)布爾量,用于描述由蓄電池額定容量限制產(chǎn)生的棄風(fēng)能量

si,sh1(t)布爾量,用于描述由放電功率限制產(chǎn)生的損失能量

si,sh1(t)布爾量,用于描述由蓄電池最小容量限制產(chǎn)生的損失能量

t時(shí)間尺度

pi(t)t時(shí)刻的風(fēng)電功率

參考輸出功率區(qū)間

pi,besscen額定充電功率

pi,ce(t)t時(shí)刻的充電功率

pi,bessden額定放電功率

pi,de(t)t時(shí)刻的放電功率

ci,bess(t)t時(shí)刻的荷電狀態(tài)

ci,bessn額定容量

dod放電深度

xi,c(t)0-1變量,表示t時(shí)刻的充電功率

xi,d(t)0-1變量,表示t時(shí)刻的放電功率

本發(fā)明釩電池為例進(jìn)行具體分析說明,其相關(guān)參數(shù)如表1所示。

表1釩電池參數(shù)

儲(chǔ)能配置的目標(biāo)函數(shù)是風(fēng)儲(chǔ)系統(tǒng)的成本最小。每個(gè)場景下的目標(biāo)函數(shù)為:

minf=ksuρ1wi,su+kshρ2wi,sh

+kin(ρcci.bessn+ρppi.bessn+rs)(17)

對(duì)于第i個(gè)場景,其中

ρ1棄風(fēng)能量對(duì)應(yīng)單價(jià)

ρ2損失能量對(duì)應(yīng)單價(jià)

ρc額定容量對(duì)應(yīng)單價(jià)

ρp額定功率對(duì)應(yīng)單價(jià)

rs儲(chǔ)能系統(tǒng)安裝成本

ksu棄風(fēng)能量成本系數(shù)

ksh損失能量成本系數(shù)

kin投資成本系數(shù)

約束條件包括蓄電池約束和輸出功率約束:

充電功率約束:

放電功率約束:

輸出功率約束:

分別計(jì)算m個(gè)場景的優(yōu)化結(jié)果,其額定容量和充放電功率如下:

其中pi為第i個(gè)場景的概率。

為驗(yàn)證提出的儲(chǔ)能配置模型,在求出式(21)—式(23)額定容量和充放電功率的基礎(chǔ)上,將其作為b中模型的已知量,利用風(fēng)電功率實(shí)測數(shù)據(jù)作為模型的輸入量(即上述模型中的pi(t)),計(jì)算棄風(fēng)成本、損失成本和最優(yōu)綜合成本。

為分析預(yù)測誤差時(shí)空相關(guān)性對(duì)風(fēng)儲(chǔ)系統(tǒng)綜合成本效益的影響以及證明上述方法的有效性和可行性,本發(fā)明對(duì)風(fēng)電場30天數(shù)據(jù),設(shè)計(jì)以下4個(gè)算例為a風(fēng)電場配置儲(chǔ)能:(a)不考慮預(yù)測誤差的時(shí)空相關(guān)性;(b)考慮ab兩風(fēng)場間預(yù)測誤差時(shí)空相關(guān)性最強(qiáng)的情況;(c)考慮ac兩風(fēng)場間預(yù)測誤差較強(qiáng)空間相關(guān)性以及風(fēng)場自身的時(shí)間相關(guān)性的情況;(d)考慮ad兩風(fēng)場間預(yù)測誤差最弱空間相關(guān)性和風(fēng)場自身的時(shí)間相關(guān)性的情況。根據(jù)第四部分提出的蓄電池儲(chǔ)能優(yōu)化模型,需要首先確定輸出功率參考值。假設(shè)輸出功率參考值為一個(gè)區(qū)間,該區(qū)間由相應(yīng)的預(yù)測誤差分布。圖7為a風(fēng)場第一天的參考功率,由圖7可以看出輸出功率的范圍隨著預(yù)測區(qū)間增加而增大。

根據(jù)基于normal/tcopula得到的多元聯(lián)合分布函數(shù),可分別得到每個(gè)算例中各個(gè)時(shí)間段的1000個(gè)場景,并利用向后縮減法將1000個(gè)場景縮減為10個(gè)。

圖8(a)-圖8(d)為場景縮減后的10個(gè)場景:圖8(a)不考慮時(shí)空相關(guān)性;圖8(b)時(shí)空相關(guān)性最強(qiáng);圖8(c)時(shí)間相關(guān)性最強(qiáng)空間相關(guān)性較強(qiáng),圖8(d)時(shí)間相關(guān)性最強(qiáng)空間相關(guān)性最弱。可以看出,所有場景均隨預(yù)測區(qū)間的增大而呈現(xiàn)發(fā)散趨勢。對(duì)比四個(gè)算例的場景縮減圖,不考慮時(shí)空相關(guān)性時(shí)的場景波動(dòng)最為劇烈,隨著空間相關(guān)性的增強(qiáng),場景波動(dòng)隨之減弱。

本發(fā)明利用gams平臺(tái)對(duì)第三部分提出的優(yōu)化模型進(jìn)行計(jì)算。充電功率、放電功率和額定容量計(jì)算結(jié)果如表3所示。由表3可以看出:(a)與其他三種算例相比,不考慮預(yù)測誤差時(shí)空相關(guān)性時(shí)額定充電功率和額定放電功率絕對(duì)值最大;(b)考慮時(shí)空相關(guān)性時(shí),額定容量比不考慮時(shí)空相關(guān)性的算例1大得多;(c)比較算例2、算例3和算例4,額定充電功率和額定放電功率的絕對(duì)值隨著空間相關(guān)性的減弱而增大;(d)考慮時(shí)空相關(guān)性時(shí),額定容量隨著空間相關(guān)性的減弱而減小;(e)基于tcopula計(jì)算得到的額定充電功率和額定放電功率的絕對(duì)值大于由normalcopula求得的結(jié)果;(f)由tcopula計(jì)算得到的額定容量比基于normalcopula得到的要小。

表3儲(chǔ)能優(yōu)化模型額定充放電功率和額定容量優(yōu)化計(jì)算結(jié)果

表4為運(yùn)行成本和最優(yōu)綜合成本。由表4可以看出,(a)算例1基于kde計(jì)算得到的運(yùn)行成本和綜合成本最大;(b)考慮預(yù)測誤差時(shí)空相關(guān)性的算例中,運(yùn)行成本與綜合成本均隨著空間相關(guān)性的減弱而增大;(c)對(duì)于算例2、算例3和算例4,基于tcopula的計(jì)算結(jié)果均小于由normalcopula計(jì)算得到的結(jié)果。

表4運(yùn)行成本和綜合成本計(jì)算結(jié)果

圖9為原始預(yù)測誤差曲線和經(jīng)過儲(chǔ)能補(bǔ)償?shù)念A(yù)測誤差曲線。由上圖可明顯看出為風(fēng)電場配置儲(chǔ)能可以有效補(bǔ)償預(yù)測誤差,且風(fēng)電場間預(yù)測誤差的空間相關(guān)性越強(qiáng),補(bǔ)償效果越好。

總之,若不考慮風(fēng)電功率短期預(yù)測誤差的時(shí)空相關(guān)性,則會(huì)導(dǎo)致額定容量被高估、充放電功率被低估。若考慮預(yù)測誤差時(shí)空相關(guān)性,不同空間相關(guān)性的綜合成本效益差異明顯,風(fēng)電場間的空間相關(guān)性越強(qiáng),風(fēng)儲(chǔ)系統(tǒng)綜合成本效益越好,對(duì)風(fēng)電功率預(yù)測誤差的補(bǔ)償效果越好。

以上所述僅為本申請(qǐng)的優(yōu)選實(shí)施例而已,并不用于限制本申請(qǐng),對(duì)于本領(lǐng)域的技術(shù)人員來說,本申請(qǐng)可以有各種更改和變化。凡在本申請(qǐng)的精神和原則之內(nèi),所作的任何修改、等同替換、改進(jìn)等,均應(yīng)包含在本申請(qǐng)的保護(hù)范圍之內(nèi)。

上述雖然結(jié)合附圖對(duì)本發(fā)明的具體實(shí)施方式進(jìn)行了描述,但并非對(duì)本發(fā)明保護(hù)范圍的限制,所屬領(lǐng)域技術(shù)人員應(yīng)該明白,在本發(fā)明的技術(shù)方案的基礎(chǔ)上,本領(lǐng)域技術(shù)人員不需要付出創(chuàng)造性勞動(dòng)即可做出的各種修改或變形仍在本發(fā)明的保護(hù)范圍以內(nèi)。

當(dāng)前第1頁1 2 
網(wǎng)友詢問留言 已有0條留言
  • 還沒有人留言評(píng)論。精彩留言會(huì)獲得點(diǎn)贊!
1
婷婷色综合www| 亚洲精品456在线播放app| 婷婷色综合www| 韩国高清视频一区二区三区| 18+在线观看网站| 亚洲欧美一区二区三区黑人 | 成人影院久久| 欧美xxxx性猛交bbbb| 精品久久久噜噜| 网址你懂的国产日韩在线| 亚洲欧美日韩无卡精品| 国产精品女同一区二区软件| 亚洲成人一二三区av| 国产精品人妻久久久久久| 男的添女的下面高潮视频| 一级毛片 在线播放| 美女福利国产在线 | 国产乱人视频| 噜噜噜噜噜久久久久久91| 女性被躁到高潮视频| 国产一区有黄有色的免费视频| 夫妻性生交免费视频一级片| 免费不卡的大黄色大毛片视频在线观看| av卡一久久| 欧美成人精品欧美一级黄| 国产伦在线观看视频一区| 这个男人来自地球电影免费观看 | 国产v大片淫在线免费观看| 欧美高清成人免费视频www| 91久久精品电影网| 久久久久久久久久久免费av| 日产精品乱码卡一卡2卡三| 中国三级夫妇交换| 啦啦啦在线观看免费高清www| 亚洲av福利一区| 日韩精品有码人妻一区| 日本wwww免费看| 十分钟在线观看高清视频www | 身体一侧抽搐| 国产女主播在线喷水免费视频网站| 国产成人精品福利久久| 赤兔流量卡办理| 久久久亚洲精品成人影院| 日本vs欧美在线观看视频 | 久久精品国产a三级三级三级| 性高湖久久久久久久久免费观看| 七月丁香在线播放| 毛片女人毛片| 日韩欧美一区视频在线观看 | 国国产精品蜜臀av免费| a级毛色黄片| 亚洲精品一二三| 久久久午夜欧美精品| 精品一区在线观看国产| 亚洲精品456在线播放app| 亚洲无线观看免费| 国产真实伦视频高清在线观看| 国产亚洲最大av| 国产精品久久久久成人av| 国产在线一区二区三区精| 亚洲色图av天堂| 少妇丰满av| 啦啦啦中文免费视频观看日本| 国产精品国产三级国产专区5o| 国产在线男女| 国产精品成人在线| 日本爱情动作片www.在线观看| 国产69精品久久久久777片| 久久毛片免费看一区二区三区| 边亲边吃奶的免费视频| 国产高清国产精品国产三级 | 亚洲真实伦在线观看| 亚洲精品国产成人久久av| 日韩 亚洲 欧美在线| 秋霞伦理黄片| 婷婷色综合大香蕉| 18禁裸乳无遮挡免费网站照片| 观看免费一级毛片| 天天躁夜夜躁狠狠久久av| 欧美少妇被猛烈插入视频| 久久精品夜色国产| 蜜桃在线观看..| 国产精品偷伦视频观看了| 久久人人爽av亚洲精品天堂 | 国产精品国产三级国产av玫瑰| 高清欧美精品videossex| 久久久久久久亚洲中文字幕| 婷婷色麻豆天堂久久| 丰满迷人的少妇在线观看| 国产无遮挡羞羞视频在线观看| 欧美+日韩+精品| 婷婷色综合大香蕉| 国语对白做爰xxxⅹ性视频网站| 丝袜喷水一区| 男女边摸边吃奶| 国产精品一及| 国产亚洲91精品色在线| 美女xxoo啪啪120秒动态图| 亚洲va在线va天堂va国产| 亚洲精品第二区| 免费黄频网站在线观看国产| av视频免费观看在线观看| 国产精品一区二区三区四区免费观看| 成年美女黄网站色视频大全免费 | 一区在线观看完整版| 国精品久久久久久国模美| 欧美成人午夜免费资源| 边亲边吃奶的免费视频| 91在线精品国自产拍蜜月| 亚洲人成网站在线播| 亚洲成色77777| 国产精品一及| 国产视频内射| 岛国毛片在线播放| 精品一区二区免费观看| 高清av免费在线| 午夜免费鲁丝| av黄色大香蕉| 女性生殖器流出的白浆| 蜜臀久久99精品久久宅男| a级毛片免费高清观看在线播放| 麻豆精品久久久久久蜜桃| 亚洲国产日韩一区二区| 亚洲av男天堂| 韩国av在线不卡| 日本黄色日本黄色录像| xxx大片免费视频| 欧美xxxx性猛交bbbb| 亚洲久久久国产精品| 亚洲av成人精品一二三区| 一级毛片我不卡| 中文字幕久久专区| 国产极品天堂在线| 啦啦啦视频在线资源免费观看| 尤物成人国产欧美一区二区三区| 亚洲精品自拍成人| 国产精品精品国产色婷婷| 国产91av在线免费观看| 日韩欧美 国产精品| 久久ye,这里只有精品| 国产精品蜜桃在线观看| 插逼视频在线观看| 日韩av在线免费看完整版不卡| 一边亲一边摸免费视频| 三级国产精品片| 女人久久www免费人成看片| 最近最新中文字幕免费大全7| 亚洲欧美中文字幕日韩二区| 国产精品一区二区三区四区免费观看| 日本猛色少妇xxxxx猛交久久| 国产视频内射| 亚洲第一区二区三区不卡| 有码 亚洲区| 美女高潮的动态| 亚洲中文av在线| 一本色道久久久久久精品综合| 日韩 亚洲 欧美在线| 不卡视频在线观看欧美| 国产精品免费大片| 精品99又大又爽又粗少妇毛片| 免费观看在线日韩| 校园人妻丝袜中文字幕| 一个人免费看片子| 黄片无遮挡物在线观看| 高清欧美精品videossex| 久久久欧美国产精品| 最黄视频免费看| 免费观看的影片在线观看| 国产午夜精品一二区理论片| 亚洲伊人久久精品综合| 国产淫片久久久久久久久| 亚洲,欧美,日韩| 亚洲国产精品专区欧美| 韩国高清视频一区二区三区| 人妻夜夜爽99麻豆av| 天堂8中文在线网| 久久久久久久久久久丰满| 插逼视频在线观看| 免费观看性生交大片5| 最近最新中文字幕大全电影3| 日韩av不卡免费在线播放| 久久青草综合色| 久久久a久久爽久久v久久| 亚洲精品日韩在线中文字幕| 亚洲精品国产av成人精品| 午夜免费男女啪啪视频观看| 国产精品秋霞免费鲁丝片| 国产精品一区www在线观看| 国产探花极品一区二区| 高清日韩中文字幕在线| 亚洲精品456在线播放app| 免费观看在线日韩| 校园人妻丝袜中文字幕| 午夜视频国产福利| 最近最新中文字幕免费大全7| 国产中年淑女户外野战色| 黄色一级大片看看| 国产精品爽爽va在线观看网站| 熟妇人妻不卡中文字幕| 国产精品成人在线| 水蜜桃什么品种好| 五月开心婷婷网| 精品久久国产蜜桃| 精品99又大又爽又粗少妇毛片| av在线蜜桃| 日本av手机在线免费观看| 一个人免费看片子| av免费观看日本| 女性被躁到高潮视频| 国产一区二区三区综合在线观看 | 91精品一卡2卡3卡4卡| 一区二区av电影网| 超碰97精品在线观看| 久久人妻熟女aⅴ| 青青草视频在线视频观看| 老师上课跳d突然被开到最大视频| 这个男人来自地球电影免费观看 | 高清不卡的av网站| 国产精品一区二区三区四区免费观看| 成人黄色视频免费在线看| 黄色怎么调成土黄色| 亚洲欧美日韩另类电影网站 | 精品亚洲成a人片在线观看 | 久久ye,这里只有精品| 日韩亚洲欧美综合| 一个人免费看片子| 在线天堂最新版资源| 九九久久精品国产亚洲av麻豆| 在线观看免费高清a一片| 狠狠精品人妻久久久久久综合| 亚洲性久久影院| 亚洲激情五月婷婷啪啪| 精品国产三级普通话版| 午夜免费鲁丝| 秋霞在线观看毛片| 日韩在线高清观看一区二区三区| 一区二区三区精品91| 免费观看无遮挡的男女| 亚洲av成人精品一区久久| 边亲边吃奶的免费视频| 亚洲精品国产色婷婷电影| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品亚洲一区二区| 99久国产av精品国产电影| 国产成人freesex在线| 国产男女内射视频| 国产一区有黄有色的免费视频| 婷婷色综合大香蕉| 日韩 亚洲 欧美在线| 亚洲成人av在线免费| 中国国产av一级| tube8黄色片| 日本色播在线视频| 国产乱来视频区| 国产成人精品久久久久久| 高清不卡的av网站| 精品少妇黑人巨大在线播放| 日韩强制内射视频| h视频一区二区三区| 国产成人91sexporn| 国产无遮挡羞羞视频在线观看| 永久免费av网站大全| 毛片一级片免费看久久久久| 亚洲国产日韩一区二区| 欧美日韩视频高清一区二区三区二| 国产色婷婷99| 十分钟在线观看高清视频www | 少妇精品久久久久久久| 亚洲天堂av无毛| 国产精品99久久99久久久不卡 | 国产精品99久久久久久久久| 狂野欧美激情性xxxx在线观看| 丝袜喷水一区| 亚洲欧洲国产日韩| 婷婷色av中文字幕| 亚洲av中文av极速乱| 亚洲精品国产av成人精品| www.av在线官网国产| 亚洲精品国产色婷婷电影| 久久久久久久亚洲中文字幕| 我的女老师完整版在线观看| 欧美极品一区二区三区四区| 国产一区二区三区av在线| 国产人妻一区二区三区在| 美女脱内裤让男人舔精品视频| 啦啦啦中文免费视频观看日本| 国产又色又爽无遮挡免| 亚州av有码| 国国产精品蜜臀av免费| 免费人成在线观看视频色| 菩萨蛮人人尽说江南好唐韦庄| 免费看日本二区| 91久久精品电影网| 亚洲精品日本国产第一区| 亚洲成人av在线免费| 少妇裸体淫交视频免费看高清| 亚洲av中文字字幕乱码综合| 久久久国产一区二区| 国产日韩欧美亚洲二区| 国产伦精品一区二区三区视频9| 免费看av在线观看网站| 亚洲精品久久久久久婷婷小说| 大香蕉97超碰在线| 老司机影院毛片| 亚洲国产精品一区三区| 久久亚洲国产成人精品v| 2018国产大陆天天弄谢| av黄色大香蕉| 久久久久久久久大av| 亚洲国产成人一精品久久久| 看非洲黑人一级黄片| 青春草国产在线视频| 久久久久久九九精品二区国产| 婷婷色综合www| 国产在线免费精品| 99九九线精品视频在线观看视频| av黄色大香蕉| 在线天堂最新版资源| kizo精华| 久久精品国产亚洲av涩爱| 汤姆久久久久久久影院中文字幕| 插阴视频在线观看视频| 国产成人精品一,二区| 久久99精品国语久久久| 王馨瑶露胸无遮挡在线观看| 国产色婷婷99| 三级国产精品欧美在线观看| 欧美性感艳星| av在线app专区| 国产成人精品久久久久久| 天堂8中文在线网| 色视频www国产| 亚洲中文av在线| 国产精品免费大片| 欧美成人精品欧美一级黄| av.在线天堂| 一级二级三级毛片免费看| 国产在线一区二区三区精| 日韩一本色道免费dvd| 女人十人毛片免费观看3o分钟| 国产一区有黄有色的免费视频| 亚洲精品久久久久久婷婷小说| 一区二区三区四区激情视频| 亚洲av成人精品一区久久| 久久99蜜桃精品久久| 最后的刺客免费高清国语| 97超碰精品成人国产| 啦啦啦在线观看免费高清www| 久久久a久久爽久久v久久| 日本vs欧美在线观看视频 | 中文字幕制服av| 国产真实伦视频高清在线观看| 性色av一级| av在线观看视频网站免费| 亚洲精品成人av观看孕妇| 国产永久视频网站| 久久精品国产亚洲网站| 亚洲欧美日韩无卡精品| 国产欧美日韩一区二区三区在线 | av在线观看视频网站免费| 韩国av在线不卡| 国产黄色免费在线视频| 美女高潮的动态| 中文欧美无线码| 美女高潮的动态| 国产精品成人在线| 亚洲av不卡在线观看| 一区在线观看完整版| 香蕉精品网在线| av国产免费在线观看| 女的被弄到高潮叫床怎么办| 日韩精品有码人妻一区| 新久久久久国产一级毛片| 99热网站在线观看| 亚洲人成网站高清观看| 免费久久久久久久精品成人欧美视频 | 国产成人午夜福利电影在线观看| 青春草亚洲视频在线观看| 中文资源天堂在线| 亚洲精品久久午夜乱码| 国产一区二区三区综合在线观看 | 我的老师免费观看完整版| 老司机影院成人| 日本黄大片高清| 国产片特级美女逼逼视频| 中文天堂在线官网| 日韩一区二区三区影片| 日本av免费视频播放| 亚洲成人一二三区av| 婷婷色综合大香蕉| 伊人久久国产一区二区| 日韩中文字幕视频在线看片 | 久久人人爽av亚洲精品天堂 | 国产精品久久久久成人av| 寂寞人妻少妇视频99o| 在线 av 中文字幕| 亚洲精品亚洲一区二区| 亚洲精品中文字幕在线视频 | 亚洲精品久久久久久婷婷小说| 男人爽女人下面视频在线观看| 成人免费观看视频高清| 超碰av人人做人人爽久久| 黄色欧美视频在线观看| 精品久久国产蜜桃| 91狼人影院| 国产黄片视频在线免费观看| 亚洲经典国产精华液单| 2021少妇久久久久久久久久久| 人妻夜夜爽99麻豆av| 亚洲欧美一区二区三区国产| 国产精品偷伦视频观看了| 美女视频免费永久观看网站| 精品久久久久久电影网| 久久久久久久亚洲中文字幕| 欧美少妇被猛烈插入视频| 欧美日本视频| 亚洲色图综合在线观看| 黑丝袜美女国产一区| 午夜福利影视在线免费观看| 国内少妇人妻偷人精品xxx网站| 秋霞在线观看毛片| 国产一区二区在线观看日韩| 国产大屁股一区二区在线视频| 自拍欧美九色日韩亚洲蝌蚪91 | 2022亚洲国产成人精品| 国产一级毛片在线| 精品一区二区免费观看| 联通29元200g的流量卡| 美女xxoo啪啪120秒动态图| 人妻制服诱惑在线中文字幕| 插阴视频在线观看视频| 日韩人妻高清精品专区| 国产无遮挡羞羞视频在线观看| 婷婷色综合大香蕉| 老师上课跳d突然被开到最大视频| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲真实伦在线观看| 久热久热在线精品观看| 七月丁香在线播放| 美女中出高潮动态图| 成人影院久久| 五月开心婷婷网| 亚洲精品国产av成人精品| 亚洲av欧美aⅴ国产| 成年免费大片在线观看| 亚洲欧美成人综合另类久久久| 国产色爽女视频免费观看| 一级a做视频免费观看| 男女免费视频国产| 国产亚洲精品久久久com| 身体一侧抽搐| 成年免费大片在线观看| 日韩制服骚丝袜av| 我要看日韩黄色一级片| 亚洲国产日韩一区二区| 日韩免费高清中文字幕av| 亚洲精品aⅴ在线观看| 久久人人爽人人片av| 欧美日韩精品成人综合77777| 又爽又黄a免费视频| 日韩av免费高清视频| 97超碰精品成人国产| 高清不卡的av网站| 制服丝袜香蕉在线| 久久久久性生活片| 搡女人真爽免费视频火全软件| 又黄又爽又刺激的免费视频.| 一级毛片电影观看| 人妻制服诱惑在线中文字幕| 我要看黄色一级片免费的| 久久久亚洲精品成人影院| 国产一级毛片在线| 国产精品国产三级国产av玫瑰| 成人午夜精彩视频在线观看| 国产av国产精品国产| 亚洲电影在线观看av| 亚洲精品第二区| 久久精品国产自在天天线| av在线老鸭窝| 成人18禁高潮啪啪吃奶动态图 | 极品少妇高潮喷水抽搐| 亚洲欧洲日产国产| 国产欧美亚洲国产| 国产精品嫩草影院av在线观看| 亚洲av中文av极速乱| 亚洲精品国产av蜜桃| 18禁裸乳无遮挡动漫免费视频| 亚洲人成网站在线观看播放| 亚洲欧美日韩卡通动漫| 亚洲欧美一区二区三区黑人 | 久久亚洲国产成人精品v| 国产成人a∨麻豆精品| 国产免费一级a男人的天堂| 国产视频内射| 日日撸夜夜添| 亚洲精品一二三| 嫩草影院入口| 久久精品人妻少妇| 亚洲人成网站在线播| 午夜激情久久久久久久| 国产精品99久久久久久久久| 日本一二三区视频观看| 日本与韩国留学比较| 国产在线一区二区三区精| 男女边吃奶边做爰视频| 久久鲁丝午夜福利片| 春色校园在线视频观看| 久久99精品国语久久久| 又爽又黄a免费视频| av一本久久久久| 伦理电影大哥的女人| 一本一本综合久久| 成人黄色视频免费在线看| 少妇精品久久久久久久| 青青草视频在线视频观看| 插逼视频在线观看| 国产精品一二三区在线看| 亚洲伊人久久精品综合| kizo精华| 一区二区三区乱码不卡18| 国产黄色免费在线视频| 97在线人人人人妻| 成人影院久久| 精品国产露脸久久av麻豆| 狂野欧美白嫩少妇大欣赏| 国内揄拍国产精品人妻在线| 国产欧美日韩精品一区二区| 91精品国产国语对白视频| 欧美精品人与动牲交sv欧美| 久久久久久久国产电影| 男女边吃奶边做爰视频| 成人亚洲精品一区在线观看 | 制服丝袜香蕉在线| 91精品国产九色| 777米奇影视久久| 亚洲在久久综合| 男人和女人高潮做爰伦理| 成人黄色视频免费在线看| 精品久久久噜噜| 免费观看性生交大片5| 三级国产精品欧美在线观看| 欧美日韩亚洲高清精品| 大片电影免费在线观看免费| 免费大片黄手机在线观看| 在线免费十八禁| 日韩精品有码人妻一区| 国产午夜精品一二区理论片| 丰满少妇做爰视频| 亚洲欧美精品专区久久| tube8黄色片| 校园人妻丝袜中文字幕| 亚洲精品aⅴ在线观看| 精品99又大又爽又粗少妇毛片| 国产亚洲午夜精品一区二区久久| 国产精品一区二区三区四区免费观看| 日韩 亚洲 欧美在线| 亚洲国产av新网站| 日韩在线高清观看一区二区三区| 日韩免费高清中文字幕av| 国产精品无大码| 夜夜爽夜夜爽视频| 91精品国产国语对白视频| 亚洲精品自拍成人| 国产大屁股一区二区在线视频| 超碰av人人做人人爽久久| 身体一侧抽搐| 久久久久久人妻| 国产亚洲最大av| 精品国产一区二区三区久久久樱花 | 人人妻人人看人人澡| videossex国产| 久久99蜜桃精品久久| 国产色婷婷99| 少妇猛男粗大的猛烈进出视频| 午夜日本视频在线| 丝袜喷水一区| 亚洲,一卡二卡三卡| 天天躁夜夜躁狠狠久久av| 高清在线视频一区二区三区| 欧美性感艳星| 1000部很黄的大片| av线在线观看网站| 一级毛片久久久久久久久女| 在线看a的网站| 成年av动漫网址| 久久影院123| 亚洲第一区二区三区不卡| 中文字幕制服av| 国产视频首页在线观看| 夜夜看夜夜爽夜夜摸| 久热久热在线精品观看| 欧美另类一区| 国产精品一区二区三区四区免费观看| 一级毛片我不卡| 日韩欧美 国产精品| 亚洲精品中文字幕在线视频 | 九色成人免费人妻av| av天堂中文字幕网| 成人二区视频| 欧美国产精品一级二级三级 | av又黄又爽大尺度在线免费看| 少妇高潮的动态图| 午夜激情福利司机影院| 国产欧美另类精品又又久久亚洲欧美| 精品酒店卫生间| 三级国产精品片| 一级片'在线观看视频| 免费av不卡在线播放| 成人二区视频| 一级毛片电影观看| 亚洲电影在线观看av| 熟妇人妻不卡中文字幕| 偷拍熟女少妇极品色| 国产爱豆传媒在线观看| 国产亚洲av片在线观看秒播厂|