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

    采用分形插值的典型日負(fù)荷曲線改進(jìn)預(yù)測(cè)方法

    2015-07-18 11:06:00李萌程浩忠楊宗麟韓新陽(yáng)楊鏡非上海交通大學(xué)電力傳輸與功率變換控制教育部重點(diǎn)實(shí)驗(yàn)室上海0040華東電網(wǎng)有限公司上海0040北京國(guó)網(wǎng)經(jīng)濟(jì)技術(shù)研究院北京0005
    關(guān)鍵詞:分形插值典型

    李萌,程浩忠,楊宗麟,韓新陽(yáng),楊鏡非(.上海交通大學(xué)電力傳輸與功率變換控制教育部重點(diǎn)實(shí)驗(yàn)室,上海0040;.華東電網(wǎng)有限公司,上海0040;.北京國(guó)網(wǎng)經(jīng)濟(jì)技術(shù)研究院,北京0005)

    采用分形插值的典型日負(fù)荷曲線改進(jìn)預(yù)測(cè)方法

    李萌1,程浩忠1,楊宗麟2,韓新陽(yáng)3,楊鏡非1
    (1.上海交通大學(xué)電力傳輸與功率變換控制教育部重點(diǎn)實(shí)驗(yàn)室,上海200240;2.華東電網(wǎng)有限公司,上海200240;3.北京國(guó)網(wǎng)經(jīng)濟(jì)技術(shù)研究院,北京100052)

    提出了一種結(jié)合粒子群算法的改進(jìn)分形預(yù)測(cè)方法。針對(duì)各年典型日負(fù)荷曲線形態(tài)相近且具有上移趨勢(shì)的特點(diǎn),采用調(diào)整向量來(lái)描述該趨勢(shì),在生成迭代函數(shù)系吸引子的過(guò)程中利用粒子群算法對(duì)調(diào)整向量進(jìn)行優(yōu)化。針對(duì)傳統(tǒng)分形預(yù)測(cè)中迭代初始點(diǎn)經(jīng)驗(yàn)性選取的問(wèn)題,提出了利用“時(shí)序平移”的思想來(lái)計(jì)算迭代初始點(diǎn)的方法。結(jié)合調(diào)整向量?jī)?yōu)化和時(shí)序平移思想,建立改進(jìn)的分形預(yù)測(cè)模型。最后,通過(guò)實(shí)例計(jì)算說(shuō)明了該方法的有效性。

    負(fù)荷曲線預(yù)測(cè);典型日負(fù)荷;分形插值;迭代函數(shù)系;粒子群算法

    中長(zhǎng)期的負(fù)荷預(yù)測(cè)對(duì)于電力設(shè)備的投資、輸變電系統(tǒng)的規(guī)劃[1]等具有前瞻性的指導(dǎo)作用,是電力系統(tǒng)規(guī)劃的基礎(chǔ),而典型日負(fù)荷曲線的預(yù)測(cè)則是中長(zhǎng)期負(fù)荷預(yù)測(cè)的重要組成部分,其預(yù)測(cè)的準(zhǔn)確性對(duì)于把握未來(lái)負(fù)荷特性、模擬系統(tǒng)生產(chǎn)以及在市場(chǎng)環(huán)境下輔助售電公司選擇優(yōu)質(zhì)用戶[2]等具有重要意義。目前典型日負(fù)荷曲線的預(yù)測(cè)思路大體可以歸納為兩類:第1類是直接預(yù)測(cè)負(fù)荷曲線形狀,并利用典型日電量、最大負(fù)荷以及負(fù)荷率等特征參數(shù)的預(yù)測(cè)值對(duì)曲線進(jìn)行修正[3-4];第2類則是通過(guò)分別獨(dú)立預(yù)測(cè)每一時(shí)刻的負(fù)荷值來(lái)預(yù)測(cè)曲線[5]。第1類方法對(duì)特征參數(shù)的預(yù)測(cè)精度要求較高,易造成預(yù)測(cè)誤差的積累;第2類方法雖然可以避免誤差累積的問(wèn)題,但由其對(duì)不同時(shí)刻負(fù)荷值獨(dú)立預(yù)測(cè),需要分別設(shè)定模型參數(shù),操作復(fù)雜,也無(wú)法考慮到負(fù)荷的短期相關(guān)性[6]。

    分形理論是非線性科學(xué)的前沿分支,近幾年被逐漸應(yīng)用到了電力負(fù)荷預(yù)測(cè)之中,具有不依賴任何負(fù)荷特征參數(shù)、操作簡(jiǎn)便的優(yōu)點(diǎn)。該理論基于部分與整體的自相似性,從部分出發(fā)來(lái)確定整體的性質(zhì),可以直接預(yù)測(cè)曲線的整體形狀。分形插值則是構(gòu)建分形曲線的方法,相比于Lagrange等插值方法更能反映出插值點(diǎn)之間的突變情況,具有較好的擬合強(qiáng)波動(dòng)性的曲線的能力[7]。文獻(xiàn)[8]利用重標(biāo)極差分析方法闡明了電力負(fù)荷的持久性,驗(yàn)證了分形插值法在電力負(fù)荷預(yù)測(cè)中的可行性;文獻(xiàn)[9]指出迭代函數(shù)系IFS(iterative function system)的迭代收斂過(guò)程與初始點(diǎn)的選取有關(guān),提出了基于均方偏差最小化的分形外推差值算法;文獻(xiàn)[10]利用分形維數(shù)計(jì)算垂直比例因子的方法,避免了計(jì)算受到人為影響。這些文獻(xiàn)的預(yù)測(cè)對(duì)象均為短期內(nèi)日負(fù)荷曲線。相比于短期負(fù)荷預(yù)測(cè),典型日的負(fù)荷曲線的預(yù)測(cè)時(shí)間跨度較長(zhǎng)。典型日負(fù)荷不但受當(dāng)日氣象特征等短期隨機(jī)因素的影響,還會(huì)受到經(jīng)濟(jì)發(fā)展、電力需求增加等長(zhǎng)期趨勢(shì)因素的影響,不同年份的典型日負(fù)荷曲線總體上呈現(xiàn)出增長(zhǎng)趨勢(shì),即后年的典型日負(fù)荷曲線一般要高于前年,造成各年典型日負(fù)荷曲線的IFS差距較大。傳統(tǒng)分形預(yù)測(cè)方法(適用于短期預(yù)測(cè))直接利用歷史曲線的迭代函數(shù)系的加權(quán)平均值來(lái)近似未來(lái)曲線的IFS,該方法應(yīng)用于典型日負(fù)荷預(yù)測(cè)時(shí)將產(chǎn)生較大誤差。

    針對(duì)傳統(tǒng)分形預(yù)測(cè)的這種局限,本文提出了利用調(diào)整向量來(lái)自適應(yīng)處理歷史負(fù)荷曲線的方法,并通過(guò)粒子群算法來(lái)優(yōu)化調(diào)整向量以減小預(yù)測(cè)誤差。利用“時(shí)序平移”[11]思想來(lái)求取最優(yōu)的迭代初始點(diǎn),進(jìn)而預(yù)測(cè)典型日負(fù)荷曲線。最后結(jié)合某省級(jí)電網(wǎng)的負(fù)荷數(shù)據(jù)進(jìn)行方法的驗(yàn)證。

    1 分形插值預(yù)測(cè)

    1.1 分形插值原理

    分形拼貼定理[7]指出,對(duì)于完備度量空間中的非空緊集S,存在一個(gè)迭代函數(shù)系IFS{R2∶wi,i= 1,2,…,N},使得該IFS的吸引子A近似于該集合,其中wi為空間R2→R2的壓縮映射,其具有的仿射變換[12]為

    即A經(jīng)過(guò)壓縮映射后的象仍為其本身。

    分形插值方法[13]即為構(gòu)造上述IFS的過(guò)程。對(duì)于選定的插值點(diǎn)集合{(xi,yi)∶i=0,1,…,N},xi表示第i個(gè)插值點(diǎn)所對(duì)應(yīng)的時(shí)刻,xi∈{1,2,…,24},yi表示xi時(shí)刻對(duì)應(yīng)的負(fù)荷值;壓縮映射wi的各參數(shù)可由插值點(diǎn)集確定,即

    式中:N為插值點(diǎn)個(gè)數(shù);di為自由參數(shù),又稱為垂直尺度因子,一般可憑經(jīng)驗(yàn)選取。為了避免di選取的隨意性,文獻(xiàn)[10]依據(jù)分形維數(shù)定理提出了一種計(jì)算方法,假設(shè)各垂直比例因子di大小相等,則有

    式中,H為Hurst指數(shù),是檢驗(yàn)電力負(fù)荷數(shù)據(jù)自相似性的重要指數(shù)。Hurst指數(shù)的計(jì)算方法在文獻(xiàn)[14-15]中有詳細(xì)的介紹,本文不再贅述。

    1.2 分形插值的預(yù)測(cè)流程

    一個(gè)系統(tǒng)的部分與整體之間某種形式的相似稱為自相似性[7]。若將負(fù)荷曲線構(gòu)成的時(shí)間序列視為一個(gè)系統(tǒng),以日為時(shí)間尺度時(shí),則該系統(tǒng)具有自相似性[9]。電力負(fù)荷具有自相似性是分形理論應(yīng)用于負(fù)荷預(yù)測(cè)的前提。分形拼貼定理說(shuō)明,如果以未來(lái)負(fù)荷曲線為集合S,那么可以通過(guò)歷史負(fù)荷曲線來(lái)構(gòu)造一個(gè)統(tǒng)計(jì)意義上的IFS[13],并且使得這個(gè)IFS的吸引子近似于未來(lái)負(fù)荷曲線。利用分形插值法進(jìn)行負(fù)荷曲線預(yù)測(cè)的過(guò)程如下。

    步驟1選取各歷史典型日負(fù)荷曲線的插值點(diǎn)對(duì)。通常可選擇曲線中的峰谷點(diǎn)、極值點(diǎn)、拐點(diǎn)等信息較為豐富的點(diǎn)[8];

    步驟2依據(jù)式(1)~式(4)計(jì)算各歷史曲線的IFS;

    步驟3對(duì)各歷史曲線的迭代函數(shù)系加權(quán)平均,形成統(tǒng)計(jì)意義上的IFS,用該IFS來(lái)近似表示未來(lái)的典型日負(fù)荷曲線的IFS;

    步驟4選取統(tǒng)計(jì)意義上的IFS的迭代初始點(diǎn),利用確定型迭代算法或隨機(jī)算法[7]計(jì)算該函數(shù)系的吸引子,所得的吸引子即為未來(lái)負(fù)荷曲線的近似值。

    文獻(xiàn)[8~10]均采用了以上的思路。在進(jìn)行短期日負(fù)荷曲線預(yù)測(cè)時(shí),由于未來(lái)曲線與歷史曲線在縱坐標(biāo)方向的位置相近,步驟3中的近似才得以成立。由于受趨勢(shì)性因素的影響,不同年份的典型日負(fù)荷曲線在縱坐標(biāo)上的位置呈現(xiàn)出上移的趨勢(shì),如果直接根據(jù)步驟1~步驟4進(jìn)行預(yù)測(cè)將會(huì)產(chǎn)生較大的誤差。因此,有必要對(duì)歷史曲線進(jìn)行適當(dāng)?shù)恼{(diào)整(如對(duì)曲線進(jìn)行上下平移),使得調(diào)整后由歷史曲線所生成的統(tǒng)計(jì)意義上的IFS可以近似表示未來(lái)曲線的IFS。本文引入粒子群優(yōu)化PSO(particle swarm optimization)算法來(lái)實(shí)現(xiàn)歷史曲線的自動(dòng)調(diào)整。

    2 調(diào)整向量的計(jì)算

    本文提出的“調(diào)整向量”的作用為:反映各歷史曲線間趨勢(shì)變動(dòng)信息,對(duì)歷史負(fù)荷曲線進(jìn)行調(diào)整,并利用調(diào)整后的曲線實(shí)現(xiàn)對(duì)未來(lái)年負(fù)荷曲線的預(yù)測(cè)。要實(shí)現(xiàn)以上預(yù)測(cè)思路必須解決2個(gè)問(wèn)題:①如何計(jì)算調(diào)整向量;②如何判斷計(jì)算所得的調(diào)整向量是合理的。

    本文采用粒子群算法來(lái)求解調(diào)整向量,以預(yù)測(cè)誤差是否達(dá)最小來(lái)判斷調(diào)整向量的合理性。

    PSO是一種進(jìn)化計(jì)算技術(shù),該算法通過(guò)構(gòu)造“粒子群”進(jìn)行系統(tǒng)尋優(yōu)。以各調(diào)整向量構(gòu)成的矩陣(簡(jiǎn)稱調(diào)整矩陣,下同)作為個(gè)體粒子,并設(shè)歷史負(fù)荷曲線的集合為{L1,L2,…,Lm},其中Li=(Li,1,Li,2,…,Li,24)為第i年典型日的24個(gè)整點(diǎn)負(fù)荷值所構(gòu)成的向量。利用該集合利用PSO求解調(diào)整向量并實(shí)現(xiàn)預(yù)測(cè)的步驟如下。

    (1)初始化參數(shù),令迭代次數(shù)k=1。

    (2)對(duì)于第k次迭代,設(shè)第p個(gè)粒子的位置為

    式中,xkpi為第i條歷史曲線的調(diào)整向量,則調(diào)整后的歷史曲線集合為=+,i=1,2,…,m}。利用該集合依據(jù)第1.2節(jié)中的步驟1~步驟4生成一條預(yù)測(cè)曲線,并計(jì)算預(yù)測(cè)誤差。

    (3)若迭代次數(shù)達(dá)設(shè)定的上限或滿足條件,

    則尋優(yōu)結(jié)束,輸出預(yù)測(cè)曲線;否則轉(zhuǎn)步驟(4)。

    (4)更新粒子位置與運(yùn)動(dòng)速度為

    3 典型日負(fù)荷曲線預(yù)測(cè)的實(shí)現(xiàn)

    本文提出預(yù)測(cè)方法主要分為2個(gè)部分:求取IFS的迭代初始點(diǎn)和預(yù)測(cè)未來(lái)年負(fù)荷曲線。

    3.1 迭代初始點(diǎn)的求取

    文獻(xiàn)[9]指出,迭代初始點(diǎn)的選取將影響IFS的迭代收斂速度。當(dāng)所選取的初始點(diǎn)落在分析的內(nèi)區(qū)間時(shí),由已知的IFS迭代產(chǎn)生的所有點(diǎn)將全部落在吸引子上。因此如果初始點(diǎn)能足夠接近待預(yù)測(cè)負(fù)荷曲線中的某點(diǎn)時(shí),迭代將加快收斂。

    在預(yù)測(cè)過(guò)程中,待預(yù)測(cè)曲線的吸引子是未知的,無(wú)法直接選擇初始點(diǎn),但如果能得到或預(yù)測(cè)得到待預(yù)測(cè)曲線上的某一點(diǎn)(該點(diǎn)必位于或者接近于吸引子),則可將該點(diǎn)作為迭代初始點(diǎn)。為了獲得迭代初始點(diǎn),本文采用了“時(shí)序平移”的思想?!皶r(shí)序平移”是指利用歷史負(fù)荷曲線構(gòu)造一個(gè)新集合

    式中,li,j為第i條歷史曲線的第j個(gè)負(fù)荷值。經(jīng)過(guò)“時(shí)序平移”,第i條歷史曲線的第1個(gè)時(shí)刻的負(fù)荷值變成了新集合中第i-1條曲線的最后一個(gè)時(shí)刻的負(fù)荷值。

    以日為時(shí)間尺度下,歷史負(fù)荷曲線集合具有自相似性,且自相似性與曲線的起點(diǎn)無(wú)關(guān),因此所生成的新曲線集合仍然具有自相似性。對(duì)新曲線集合利用PSO算法優(yōu)化其調(diào)整向量,設(shè)在一次迭代尋優(yōu)過(guò)程中生成的曲線為=(,…),則PSO算法的優(yōu)化目標(biāo)為的前23個(gè)點(diǎn)與的后23個(gè)點(diǎn)的相對(duì)誤差的平均值達(dá)最小,即

    當(dāng)PSO算法尋優(yōu)結(jié)束后,由于電力負(fù)荷具有自相似性與持久性,的最后一個(gè)點(diǎn)將近似于待預(yù)測(cè)負(fù)荷曲線的第1個(gè)時(shí)刻的負(fù)荷值,于是可將點(diǎn)(1,)作為迭代初始點(diǎn)。歷史曲線集合L、新曲線集合Lnew和迭代初始點(diǎn)在時(shí)間軸上的位置關(guān)系如圖1所示。

    圖1 L、Lnew以及初始點(diǎn)在時(shí)間軸上的位置Fig.1 L、Lnewand startpointon time axes

    圖2 本文預(yù)測(cè)算法流程Fig.2 Flow chartofproposed forecasting calculation

    3.2 負(fù)荷曲線的預(yù)測(cè)

    4 算例分析

    以某省2008—2011年夏季典型日負(fù)荷曲線為歷史數(shù)據(jù),對(duì)2012年夏季典型日負(fù)荷進(jìn)行預(yù)測(cè)。算法流程如圖2所示。本例中,選取形態(tài)與夏季平均日負(fù)荷曲線最為接近,且無(wú)畸變的典型日負(fù)荷曲線。

    本文在計(jì)算中發(fā)現(xiàn),分形插值點(diǎn)的選取對(duì)預(yù)測(cè)的精度和迭代收斂速度有著較大影響。插值點(diǎn)選取過(guò)少,將無(wú)法描述曲線的細(xì)節(jié)變化,將產(chǎn)生較大的預(yù)測(cè)誤差,插值點(diǎn)選取過(guò)多則將影響迭代速度。為了兼顧預(yù)測(cè)精度與迭代速度,本文選定兩組插值點(diǎn)分別進(jìn)行預(yù)測(cè):第1次插值的點(diǎn)為1、5、7、12、13、20、22共7個(gè)整點(diǎn)負(fù)荷,第2次插值的點(diǎn)為1,7,12,13,18,21,24共7個(gè)整點(diǎn)負(fù)荷。

    預(yù)測(cè)的第1步為計(jì)算第2步所需要的迭代初始點(diǎn)。利用2008-2010年每年24個(gè)時(shí)刻的負(fù)荷值以及2011年第1個(gè)時(shí)刻負(fù)荷值,采用本文的方法預(yù)測(cè)2011年的后23個(gè)時(shí)刻和2012年第1個(gè)時(shí)刻的負(fù)荷值。預(yù)測(cè)中插值點(diǎn)對(duì)按上文選取,迭代初始點(diǎn)取為2011年第2個(gè)時(shí)刻負(fù)荷值。預(yù)測(cè)得到2012年夏季典型日負(fù)荷曲線的第1時(shí)刻的負(fù)荷值為15 362 MW。

    預(yù)測(cè)的第2步為2012年負(fù)荷曲線的預(yù)測(cè)。利用第1步的預(yù)測(cè)結(jié)果,設(shè)定此步中的迭代初始點(diǎn)對(duì)為(1,15 362),插值點(diǎn)對(duì)與上一步相同。2012年夏季典型日的負(fù)荷曲線預(yù)測(cè)結(jié)果見(jiàn)表1。

    為了說(shuō)明本文提出的分形預(yù)測(cè)改進(jìn)方法的有效性,本文還利用不同方法對(duì)2012年的典型日負(fù)荷曲線進(jìn)行預(yù)測(cè),由于篇幅有限僅展示各類方法的預(yù)測(cè)曲線如圖3所示,其中方法1為傳統(tǒng)分形方法,方法2在方法1的預(yù)測(cè)結(jié)果之上進(jìn)行了電量修正,方法3為文獻(xiàn)[4]介紹的方法。不同方法的預(yù)測(cè)誤差見(jiàn)表2,具體結(jié)果說(shuō)明如下。

    (1)本文提出的相對(duì)誤差的平均值為1.671%,低于其他3種方法。相比于傳統(tǒng)分形方法而言,其預(yù)測(cè)誤差大幅度降低。傳統(tǒng)分形方法因無(wú)法考慮各年間負(fù)荷曲線的增長(zhǎng)信息,當(dāng)僅通過(guò)歷史負(fù)荷曲線生成的統(tǒng)計(jì)意義上的IFS進(jìn)行預(yù)測(cè)時(shí)將產(chǎn)生較大的誤差。

    (2)方法2在方法1的基礎(chǔ)上進(jìn)行了電量修正,本例中利用2012年夏季典型日用電量的實(shí)際值對(duì)方法1中獲得的曲線進(jìn)行比例縮放,使修正后曲線對(duì)應(yīng)的日用電量與實(shí)際日用電量相等。進(jìn)行電量修正后,預(yù)測(cè)精度得到提高。但由于該方法需提前預(yù)測(cè)日用電量,當(dāng)日用電量的預(yù)測(cè)存在誤差時(shí),該方法的誤差將增大。

    (3)方法3在個(gè)別時(shí)刻負(fù)荷值的預(yù)測(cè)精度高于本文提出的方法,但其并不能較好地反映實(shí)際曲線的形態(tài),總體上預(yù)測(cè)精度較低。

    表1 某省電網(wǎng)2012年夏季典型日負(fù)荷預(yù)測(cè)值Tab.1 Load forecasting ofsummer typicaldaily load of one certain province in 2012

    圖3 不同預(yù)測(cè)方法的結(jié)果對(duì)比Fig.3 Result comparison ofdifferentforecasting methods

    表2 不同預(yù)測(cè)方法的誤差對(duì)比Tab.2 Error comparison of different forecasting methods %

    5 結(jié)論

    (1)電力負(fù)荷曲線的自相似性與時(shí)間起點(diǎn)的選取無(wú)關(guān)是本文利用“時(shí)序平移”來(lái)求取迭代初始點(diǎn)的理論基礎(chǔ)。

    (2)傳統(tǒng)分形插值方法無(wú)法直接應(yīng)用于典型日負(fù)荷曲線預(yù)測(cè),利用調(diào)整向量對(duì)歷史曲線進(jìn)行處理,可以克服這種局限。算例證明了本文提出的改進(jìn)分形預(yù)測(cè)方法具有較高的精度,可以滿足工程實(shí)踐的要求。

    (3)與神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)等方法相比,本文所提算法對(duì)樣本容量要求低,可同時(shí)預(yù)測(cè)1 d 24點(diǎn)的負(fù)荷值。此外,本文的方法并不需要提前對(duì)電量、負(fù)荷率等參數(shù)進(jìn)行預(yù)測(cè),具有較強(qiáng)的操作性。

    [1]Hyndman R J,Shu Fan.Density forecasting for long-term peak electricity demand[J].IEEE Trans on Power Systems,2010,25(2):1142-1153.

    [2]楊科,譚倫農(nóng)(Yang Ke,Tan Lunnong).電力市場(chǎng)環(huán)境下中長(zhǎng)期負(fù)荷預(yù)測(cè)的應(yīng)用分析(Analysis oflong-term load forecasting application in power market environment)[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)(Proceedings ofthe CSU-EPSA),2011,23(3):54-57.

    [3]趙暉(Zhao Hui).用樣條插值法模擬典型日負(fù)荷曲線(Simulation oftypicaldaily load curve with spline interpolation)[J].電網(wǎng)技術(shù)(Power System Technology),1998,22(5):39-41,45.

    [4]何光宇,郭家春,鮑毅,等(He Guangyu,Guo Jiachun,Bao Yi,etal).一種中長(zhǎng)期日負(fù)荷曲線預(yù)測(cè)的新方法-雙向夾逼法(A novel method of middle and long term forecasting for daily load curve-approaching algorithm in tow directions)[J].電網(wǎng)技術(shù)(Power System Technology),2004,28(7):27-29.

    [5]蔡國(guó)偉,杜毅,李春山,等(Cai Guowei,Du Yi,Li Chunshan,etal).基于支持向量機(jī)的中長(zhǎng)期日負(fù)荷曲線預(yù)測(cè)(Middle and long-term daily load curve forecasting based on support vector machine)[J].電網(wǎng)技術(shù)(Power System Technology),2006,30(23):56-60.

    [6]Al-Hamadi H M,Soliman S A.Long-term/mid-term electric load forecasting based on short-term correlation and annualgrowth[J].Electric Power Systems Research,2005,74(3):353-361.

    [7]孫洪泉.分形與分形插值[M].北京:科學(xué)出版社,2011.

    [8]李翔,關(guān)勇,喬艷芬(Li Xiang,Guan Yong,Qiao Yanfen).基于分形理論的電力負(fù)荷持久性分析及預(yù)測(cè)(Fractal theory based durability analysis and forecasting of power load)[J].電網(wǎng)技術(shù)(Power System Technology),2006,30(16):84-88.

    [9]薛萬(wàn)磊,于繼來(lái)(Xue Wanlei,Yu Jilai).分形外推插值算法在電力負(fù)荷預(yù)測(cè)中的應(yīng)用(Application of fractal extrapolation algorithm in lfad Forecasting)[J].電網(wǎng)技術(shù)(Power System Technology),2006,30(13):49-54.

    [10]劉映,翟明岳,劉永笑,等(Liu Ying,Zhai Mingyue,Liu Yongxiao,etal).一種改進(jìn)的基于分形插值和小波分析的電力負(fù)荷短期預(yù)測(cè)方法(An improved short-term load forecasting based on fractal interpretation and wavelet analysis)[J].電網(wǎng)技術(shù)(Power System Technology),2013,37(1):276-280.

    [11]Xie Kaigui,Billiton Rov.Considering wind speed correlation ofWECS in reliability evaluation using the time-shifting technique[J].Electric Power Systems Research,2009,79(4):687-693.

    [12]張宏偉,陸仁強(qiáng),牛志廣(Zhang Hongwei,Lu Renqiang,Niu Zhiguang).基于分形理論的城市日用水量預(yù)測(cè)方法(Prediction method of urban daily water consumption based on fractional theory)[J].天津大學(xué)學(xué)報(bào)(Journalof Tianjin University),2009,42(1):56-59.

    [13]樊福梅,梁平(Fan Fumei,Liang Ping).基于分形的社會(huì)總用電量及其構(gòu)成預(yù)測(cè)(Forecasting about national electric consumption and its constitution based on fractal)[J].中國(guó)電機(jī)工程學(xué)報(bào)(Proceedings ofthe CSEE),2004,24(11):91-95.

    [14]Hurst H E.The long-term scorage capacity ofreservoirs[J].Transactions of the American Society of Civil Engineers,1951,(116):770-799.

    [15]史永東(Shi Yongdong).上海證券市場(chǎng)的分形結(jié)構(gòu)(Structure of stock returns in Shanghai stock market)[J].預(yù)測(cè)(Forecasting),2000,19(5):78-80,50.

    [16]劉玲,嚴(yán)登俊,龔燈才,等(Liu Ling,Yan Dengjun,Gong Dengcai,etal).基于粒子群模糊神經(jīng)網(wǎng)絡(luò)的短期電力負(fù)荷預(yù)測(cè)(New method for short term load forecasting based on particle swarm optimization and fuzzy neuralnetwork)[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào)(Proceedings of the CSU-EPSA),2006,18(3):47-50.

    Improved Forecasting Method of TypicalDaily Load Curve Based on FractalInterpolation

    LIMeng1,CHENG Haozhong1,YANG Zonglin2,HAN Xinyang3,YANG Jingfei1
    (1.Key Laboratory ofControlofPower Transmission and Transformation ofMinistry ofeducation,ShanghaiJiao Tong University,Shanghai200240,China;2.EastChina Grid Company Limited,Shanghai200240,China;3.State Power Economic Research Institute,Beijing 100052,China)

    This paper proposes an improved method of forecasting typical daily load curve combined with particle swarm optimization.In point of the typicaldaily load curve characteristics of similar shape and upward trend year by year,this paper adopts adjustment vector to express the trend,and optimizes the adjustment vector using particle swarm algorithm in the progress ofgenerating attractors ofiterative function system.Aiming atthe experientialselection problems ofthe initialiteration pointin the traditionalfractalprediction,a method based on time-shifting technique is proposed to calculate the initialpointfor the fractalinterpretation iteration.Combined with the adjustmentvector optimization and time-shifting thought,an improved fractal prediction model is then established.Case study results show thatthe proposed algorithm possesses advantages ofmore accurate forecasting results.

    load curve forecasting;typical day load;fractal interpretation;iterative function system;particle swarm optimization(PSO)

    TM715

    A

    1003-8930(2015)03-0036-06

    10.3969/j.issn.1003-8930.2015.03.07

    李萌(1988—),女,碩士研究生,研究方向?yàn)殡娋W(wǎng)規(guī)劃與優(yōu)化、電力系統(tǒng)負(fù)荷預(yù)測(cè)。Email:lmfighting@sina.com

    2013-08-12;

    2013-11-15

    2012年度上海市優(yōu)秀學(xué)術(shù)帶頭人計(jì)劃項(xiàng)目(12XD1402900);上海領(lǐng)軍人才資助項(xiàng)目(2012020)

    程浩忠(1962—),男,博士,教授,研究方向?yàn)殡娋W(wǎng)規(guī)劃、電壓穩(wěn)定性、電能質(zhì)量。Email:hzcheng@sjtu.edu.cn

    楊宗麟(1956—),男,本科,高級(jí)工程師,研究方向?yàn)殡娏κ袌?chǎng)分析預(yù)測(cè)、電力規(guī)劃。Email:wmy0817@163.com

    猜你喜歡
    分形插值典型
    用最典型的事寫(xiě)最有特點(diǎn)的人
    多項(xiàng)式求值題的典型解法
    典型胰島素瘤1例報(bào)道
    感受分形
    分形之美
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    分形空間上廣義凸函數(shù)的新Simpson型不等式及應(yīng)用
    一種改進(jìn)FFT多譜線插值諧波分析方法
    基于四項(xiàng)最低旁瓣Nuttall窗的插值FFT諧波分析
    Blackman-Harris窗的插值FFT諧波分析與應(yīng)用
    成人18禁高潮啪啪吃奶动态图| 久久婷婷人人爽人人干人人爱 | 校园春色视频在线观看| 中文字幕精品免费在线观看视频| 国产精品影院久久| 欧美日韩中文字幕国产精品一区二区三区 | 日本vs欧美在线观看视频| 91老司机精品| 一区二区三区激情视频| 他把我摸到了高潮在线观看| 一边摸一边抽搐一进一小说| 色av中文字幕| 亚洲av五月六月丁香网| 久久精品国产亚洲av香蕉五月| 性色av乱码一区二区三区2| 女警被强在线播放| 啦啦啦观看免费观看视频高清 | 操出白浆在线播放| 热re99久久国产66热| 久久久久久免费高清国产稀缺| 国产免费av片在线观看野外av| 亚洲欧美一区二区三区黑人| 校园春色视频在线观看| 一级黄色大片毛片| 精品人妻在线不人妻| 18禁美女被吸乳视频| 97超级碰碰碰精品色视频在线观看| 欧美最黄视频在线播放免费| 亚洲全国av大片| www.999成人在线观看| 色播亚洲综合网| 狂野欧美激情性xxxx| 免费一级毛片在线播放高清视频 | 欧美色欧美亚洲另类二区 | 亚洲精品中文字幕一二三四区| 欧美色欧美亚洲另类二区 | 日韩av在线大香蕉| 色尼玛亚洲综合影院| 亚洲色图av天堂| 日韩 欧美 亚洲 中文字幕| 制服丝袜大香蕉在线| 日韩有码中文字幕| 一级毛片女人18水好多| 日韩精品青青久久久久久| 国产一区二区在线av高清观看| 给我免费播放毛片高清在线观看| 高潮久久久久久久久久久不卡| 午夜两性在线视频| 国产黄a三级三级三级人| 国产精品久久久av美女十八| 久久久久久久精品吃奶| 色综合欧美亚洲国产小说| 午夜福利影视在线免费观看| 美女大奶头视频| 丝袜美腿诱惑在线| 欧美成人免费av一区二区三区| 国产熟女xx| 亚洲伊人色综图| 淫秽高清视频在线观看| 国产单亲对白刺激| 欧美不卡视频在线免费观看 | 黄色 视频免费看| 18禁国产床啪视频网站| 成年人黄色毛片网站| 精品久久久久久成人av| 欧美成人午夜精品| 国产欧美日韩一区二区三| 久久久久国内视频| av天堂久久9| 又紧又爽又黄一区二区| 国产精品二区激情视频| 精品熟女少妇八av免费久了| netflix在线观看网站| 亚洲伊人色综图| 国产精品自产拍在线观看55亚洲| 国产97色在线日韩免费| 美国免费a级毛片| 国产三级在线视频| 久久亚洲真实| 久久午夜综合久久蜜桃| 午夜精品国产一区二区电影| 日韩欧美免费精品| 欧美国产日韩亚洲一区| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品美女久久久久99蜜臀| 亚洲成av人片免费观看| 久久国产亚洲av麻豆专区| 999久久久国产精品视频| 欧美激情极品国产一区二区三区| 久久这里只有精品19| 精品高清国产在线一区| 久久人人精品亚洲av| 香蕉丝袜av| 少妇的丰满在线观看| 香蕉久久夜色| 国产免费av片在线观看野外av| 又大又爽又粗| 人人澡人人妻人| 悠悠久久av| 亚洲中文av在线| 国产精品秋霞免费鲁丝片| 亚洲午夜理论影院| 男女下面进入的视频免费午夜 | 久久久久国内视频| 我的亚洲天堂| 午夜激情av网站| 91成人精品电影| 变态另类丝袜制服| 中亚洲国语对白在线视频| 午夜成年电影在线免费观看| 91九色精品人成在线观看| 天天躁夜夜躁狠狠躁躁| 黑丝袜美女国产一区| 久久草成人影院| 亚洲av五月六月丁香网| 国产黄a三级三级三级人| 精品国产一区二区久久| 亚洲va日本ⅴa欧美va伊人久久| 精品国产一区二区三区四区第35| 欧美成狂野欧美在线观看| 精品国内亚洲2022精品成人| 亚洲欧美日韩无卡精品| 久久久久久久久久久久大奶| 久久人人爽av亚洲精品天堂| 亚洲少妇的诱惑av| 日本a在线网址| 欧美中文综合在线视频| 少妇粗大呻吟视频| 久久久久久大精品| www国产在线视频色| 97人妻天天添夜夜摸| 国产精品,欧美在线| 麻豆国产av国片精品| 黄色毛片三级朝国网站| 男女下面进入的视频免费午夜 | 日韩有码中文字幕| 激情在线观看视频在线高清| 国产午夜精品久久久久久| 国产精品秋霞免费鲁丝片| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产综合久久久| 亚洲国产毛片av蜜桃av| 一边摸一边抽搐一进一小说| 色综合婷婷激情| 欧美在线黄色| 18禁裸乳无遮挡免费网站照片 | 女人爽到高潮嗷嗷叫在线视频| 免费少妇av软件| 亚洲精华国产精华精| 黑人欧美特级aaaaaa片| 国产欧美日韩一区二区三区在线| 啪啪无遮挡十八禁网站| xxx96com| 美女大奶头视频| 欧美成人一区二区免费高清观看 | 视频区欧美日本亚洲| 中文字幕色久视频| 激情视频va一区二区三区| 久久香蕉国产精品| 亚洲男人的天堂狠狠| 亚洲熟女毛片儿| 老汉色av国产亚洲站长工具| 中出人妻视频一区二区| 亚洲成人免费电影在线观看| 免费观看人在逋| 日韩欧美在线二视频| 国产午夜福利久久久久久| 久久久国产精品麻豆| 色av中文字幕| 成人av一区二区三区在线看| 亚洲九九香蕉| 9色porny在线观看| 国产在线观看jvid| www.www免费av| 怎么达到女性高潮| av欧美777| 欧美成狂野欧美在线观看| 国产亚洲精品av在线| 日韩精品免费视频一区二区三区| 国产亚洲av嫩草精品影院| 亚洲av成人一区二区三| 久久天堂一区二区三区四区| 女人精品久久久久毛片| 夜夜看夜夜爽夜夜摸| 国产欧美日韩精品亚洲av| 性少妇av在线| 日韩中文字幕欧美一区二区| 丝袜美腿诱惑在线| 欧美大码av| 少妇粗大呻吟视频| 极品教师在线免费播放| 久久久久久久久免费视频了| 91成人精品电影| 丝袜人妻中文字幕| 熟妇人妻久久中文字幕3abv| 黄色视频不卡| 欧美国产日韩亚洲一区| 欧美+亚洲+日韩+国产| 亚洲黑人精品在线| 丁香六月欧美| 1024香蕉在线观看| 欧美中文日本在线观看视频| 久久久久久免费高清国产稀缺| 99国产综合亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 好男人电影高清在线观看| 久久精品国产亚洲av高清一级| 久久热在线av| 一a级毛片在线观看| 午夜精品久久久久久毛片777| 搞女人的毛片| 亚洲精品美女久久av网站| 欧美成人免费av一区二区三区| 亚洲精品在线美女| 一个人观看的视频www高清免费观看 | 夜夜看夜夜爽夜夜摸| 久久中文字幕人妻熟女| 一级毛片女人18水好多| 免费av毛片视频| 欧美色欧美亚洲另类二区 | 无人区码免费观看不卡| 日本黄色视频三级网站网址| 国产精品久久久久久人妻精品电影| 国产在线精品亚洲第一网站| 国产不卡一卡二| 午夜精品久久久久久毛片777| www日本在线高清视频| 国产精品二区激情视频| 亚洲avbb在线观看| 嫩草影视91久久| av超薄肉色丝袜交足视频| 免费一级毛片在线播放高清视频 | 国产99久久九九免费精品| 日本欧美视频一区| 精品熟女少妇八av免费久了| 一区二区三区激情视频| 电影成人av| 真人做人爱边吃奶动态| 国产精品久久久久久精品电影 | 黑人操中国人逼视频| 宅男免费午夜| 国产激情久久老熟女| 午夜精品国产一区二区电影| 成人手机av| 欧美绝顶高潮抽搐喷水| 免费高清在线观看日韩| 亚洲精品av麻豆狂野| 久久亚洲精品不卡| 亚洲成国产人片在线观看| 夜夜看夜夜爽夜夜摸| 免费女性裸体啪啪无遮挡网站| 淫妇啪啪啪对白视频| 亚洲国产精品sss在线观看| 久久午夜亚洲精品久久| 又紧又爽又黄一区二区| 69av精品久久久久久| 麻豆一二三区av精品| 欧美午夜高清在线| 久久草成人影院| 精品电影一区二区在线| 久久精品91无色码中文字幕| 搞女人的毛片| 久久中文字幕人妻熟女| 色综合婷婷激情| 岛国视频午夜一区免费看| 欧美成人免费av一区二区三区| 大型黄色视频在线免费观看| 欧美色视频一区免费| 国产乱人伦免费视频| 天天添夜夜摸| 热re99久久国产66热| 99热只有精品国产| svipshipincom国产片| 最好的美女福利视频网| bbb黄色大片| 亚洲欧美激情在线| 欧美成人性av电影在线观看| 在线国产一区二区在线| 老汉色av国产亚洲站长工具| xxx96com| 国产91精品成人一区二区三区| 可以免费在线观看a视频的电影网站| 亚洲国产欧美日韩在线播放| 在线观看午夜福利视频| 精品久久久精品久久久| 久久精品影院6| 日韩欧美国产在线观看| 国产av一区二区精品久久| 国产97色在线日韩免费| 99国产极品粉嫩在线观看| 视频在线观看一区二区三区| 亚洲精品av麻豆狂野| 美女午夜性视频免费| 在线观看日韩欧美| 免费在线观看视频国产中文字幕亚洲| 亚洲熟妇中文字幕五十中出| 91老司机精品| 美女高潮喷水抽搐中文字幕| www.熟女人妻精品国产| 成人三级黄色视频| 夜夜夜夜夜久久久久| 精品电影一区二区在线| 久99久视频精品免费| 亚洲电影在线观看av| 69av精品久久久久久| 一本综合久久免费| 欧美日韩亚洲综合一区二区三区_| 亚洲色图综合在线观看| 99国产精品免费福利视频| 国产精品 欧美亚洲| 亚洲国产精品999在线| 非洲黑人性xxxx精品又粗又长| 在线观看www视频免费| 久久久久久久久久久久大奶| 成年版毛片免费区| 日韩av在线大香蕉| 老司机深夜福利视频在线观看| 国产成人精品在线电影| 久久国产亚洲av麻豆专区| 国产精品一区二区三区四区久久 | 日韩大尺度精品在线看网址 | 一区二区三区国产精品乱码| 国产三级黄色录像| 99国产精品99久久久久| 成人亚洲精品一区在线观看| 久久 成人 亚洲| 日本精品一区二区三区蜜桃| 我的亚洲天堂| 老汉色av国产亚洲站长工具| 脱女人内裤的视频| АⅤ资源中文在线天堂| 怎么达到女性高潮| 在线天堂中文资源库| xxx96com| 国产私拍福利视频在线观看| 天堂影院成人在线观看| 午夜免费观看网址| 午夜久久久久精精品| 国产主播在线观看一区二区| 日本在线视频免费播放| 国产免费av片在线观看野外av| 国产亚洲欧美精品永久| 一a级毛片在线观看| 老司机午夜福利在线观看视频| 亚洲伊人色综图| 国产成人av教育| 久久久久久人人人人人| bbb黄色大片| 一进一出抽搐gif免费好疼| 大香蕉久久成人网| 国产99白浆流出| 男人的好看免费观看在线视频 | 精品不卡国产一区二区三区| 亚洲精品久久国产高清桃花| 一边摸一边做爽爽视频免费| 18禁裸乳无遮挡免费网站照片 | 久久久久久人人人人人| bbb黄色大片| ponron亚洲| 久久香蕉激情| 亚洲专区中文字幕在线| 亚洲人成伊人成综合网2020| 99国产极品粉嫩在线观看| 久久久精品国产亚洲av高清涩受| av中文乱码字幕在线| 无人区码免费观看不卡| 亚洲少妇的诱惑av| 中文字幕色久视频| 男人舔女人下体高潮全视频| 亚洲男人的天堂狠狠| 男女下面进入的视频免费午夜 | xxx96com| 老汉色av国产亚洲站长工具| 丁香六月欧美| 国产一区二区三区在线臀色熟女| 亚洲一区二区三区不卡视频| 国产亚洲精品久久久久5区| 亚洲专区字幕在线| 狠狠狠狠99中文字幕| 女人爽到高潮嗷嗷叫在线视频| 嫩草影院精品99| 国产亚洲精品av在线| 大陆偷拍与自拍| 好看av亚洲va欧美ⅴa在| 久久中文字幕人妻熟女| 日韩国内少妇激情av| 狠狠狠狠99中文字幕| 日韩精品中文字幕看吧| 啪啪无遮挡十八禁网站| 在线观看免费视频日本深夜| 两性夫妻黄色片| 不卡一级毛片| 无人区码免费观看不卡| 村上凉子中文字幕在线| 91麻豆av在线| 黄色女人牲交| 亚洲在线自拍视频| 在线天堂中文资源库| 黄色毛片三级朝国网站| av在线天堂中文字幕| 精品久久久精品久久久| 亚洲av熟女| 免费无遮挡裸体视频| 国产精品 国内视频| 在线观看免费视频日本深夜| 真人一进一出gif抽搐免费| 欧美日韩黄片免| 欧美色视频一区免费| 国产欧美日韩一区二区三区在线| av片东京热男人的天堂| 亚洲国产看品久久| 亚洲欧洲精品一区二区精品久久久| 亚洲激情在线av| 制服人妻中文乱码| 亚洲精品粉嫩美女一区| 久久久久久久久中文| 国产精品av久久久久免费| 超碰成人久久| 麻豆av在线久日| 一区二区三区激情视频| 国产欧美日韩综合在线一区二区| 99国产极品粉嫩在线观看| 女人被狂操c到高潮| 99re在线观看精品视频| 午夜久久久久精精品| 91精品三级在线观看| 国产精品野战在线观看| 少妇裸体淫交视频免费看高清 | 十分钟在线观看高清视频www| 天天躁狠狠躁夜夜躁狠狠躁| 欧美另类亚洲清纯唯美| 午夜福利,免费看| 国产精品电影一区二区三区| or卡值多少钱| 色婷婷久久久亚洲欧美| 日韩高清综合在线| 午夜精品国产一区二区电影| 国产精品亚洲av一区麻豆| 国产xxxxx性猛交| √禁漫天堂资源中文www| 久久精品国产清高在天天线| 亚洲欧美日韩高清在线视频| 久久久精品欧美日韩精品| netflix在线观看网站| 搡老妇女老女人老熟妇| 国产亚洲av高清不卡| 99在线视频只有这里精品首页| 韩国精品一区二区三区| 亚洲欧美精品综合久久99| 黄色女人牲交| 亚洲 欧美一区二区三区| 久久精品人人爽人人爽视色| www.999成人在线观看| 日韩大尺度精品在线看网址 | 国产亚洲精品一区二区www| 18禁美女被吸乳视频| 久久国产精品影院| 欧美日韩福利视频一区二区| 久久久国产欧美日韩av| 在线播放国产精品三级| 久久久水蜜桃国产精品网| 中文字幕av电影在线播放| 亚洲九九香蕉| 亚洲专区国产一区二区| 在线观看免费午夜福利视频| 精品久久久久久,| 国产精品久久电影中文字幕| www.999成人在线观看| 亚洲国产欧美网| 满18在线观看网站| 亚洲成人免费电影在线观看| 国产精品秋霞免费鲁丝片| 久久国产精品男人的天堂亚洲| 亚洲第一电影网av| 国产午夜精品久久久久久| 无限看片的www在线观看| 中亚洲国语对白在线视频| 欧美亚洲日本最大视频资源| 村上凉子中文字幕在线| 十分钟在线观看高清视频www| 日本 av在线| 黄片播放在线免费| 俄罗斯特黄特色一大片| 一二三四在线观看免费中文在| 男人操女人黄网站| 精品久久久久久,| 亚洲精品国产色婷婷电影| 黄色视频不卡| 亚洲最大成人中文| 久久久久九九精品影院| 亚洲伊人色综图| 亚洲男人的天堂狠狠| 91精品三级在线观看| 好看av亚洲va欧美ⅴa在| 色精品久久人妻99蜜桃| 久久久久久人人人人人| 精品日产1卡2卡| 国产精品,欧美在线| 丝袜美足系列| 国产精品 国内视频| 女人被躁到高潮嗷嗷叫费观| 欧美一级毛片孕妇| 99精品欧美一区二区三区四区| 色综合婷婷激情| 啪啪无遮挡十八禁网站| 老熟妇仑乱视频hdxx| 欧美一级毛片孕妇| 99国产精品免费福利视频| 女人爽到高潮嗷嗷叫在线视频| 一级,二级,三级黄色视频| 黄色 视频免费看| 在线av久久热| 欧美激情 高清一区二区三区| 成人欧美大片| 亚洲av熟女| 欧美亚洲日本最大视频资源| 老司机午夜福利在线观看视频| 免费在线观看完整版高清| 黑人操中国人逼视频| 一级,二级,三级黄色视频| 999久久久国产精品视频| 中出人妻视频一区二区| www.999成人在线观看| 一本大道久久a久久精品| 99riav亚洲国产免费| 亚洲伊人色综图| 国产精品永久免费网站| 国产成人精品在线电影| 在线观看66精品国产| 热99re8久久精品国产| 精品一品国产午夜福利视频| 欧美最黄视频在线播放免费| 国产成人精品无人区| 身体一侧抽搐| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲久久久国产精品| 搞女人的毛片| 男人操女人黄网站| 国产一区在线观看成人免费| 中文字幕久久专区| 99国产精品一区二区蜜桃av| 99国产极品粉嫩在线观看| √禁漫天堂资源中文www| 亚洲中文av在线| 777久久人妻少妇嫩草av网站| 久99久视频精品免费| 一区福利在线观看| 中文字幕高清在线视频| 18美女黄网站色大片免费观看| 国产精品久久久久久亚洲av鲁大| 亚洲欧洲精品一区二区精品久久久| 久久久精品欧美日韩精品| 国产亚洲av高清不卡| 两人在一起打扑克的视频| 欧美黄色片欧美黄色片| 又黄又粗又硬又大视频| 午夜精品在线福利| 日本精品一区二区三区蜜桃| 国产精品乱码一区二三区的特点 | 国产精品影院久久| 欧美大码av| 国内毛片毛片毛片毛片毛片| 在线观看www视频免费| 老司机午夜福利在线观看视频| 老汉色av国产亚洲站长工具| 亚洲av五月六月丁香网| 国产成人一区二区三区免费视频网站| 国产欧美日韩一区二区三| 麻豆一二三区av精品| 午夜福利在线观看吧| 女人高潮潮喷娇喘18禁视频| 九色国产91popny在线| 性色av乱码一区二区三区2| 丝袜在线中文字幕| 亚洲男人的天堂狠狠| 淫秽高清视频在线观看| 中文亚洲av片在线观看爽| 国产成人欧美在线观看| 久久精品影院6| 国产乱人伦免费视频| 日本vs欧美在线观看视频| 亚洲男人天堂网一区| www日本在线高清视频| 免费在线观看黄色视频的| 久久精品国产亚洲av香蕉五月| 欧美日本视频| 精品欧美国产一区二区三| av欧美777| 精品国产国语对白av| 国产亚洲av嫩草精品影院| 午夜福利欧美成人| 一区二区三区国产精品乱码| 国产成人影院久久av| 97人妻精品一区二区三区麻豆 | 国产熟女xx| 男人操女人黄网站| 日韩欧美一区视频在线观看| 国产熟女xx| 99在线视频只有这里精品首页| 一区二区日韩欧美中文字幕| 精品久久久精品久久久| 一区在线观看完整版| 国产av一区二区精品久久| 国产伦一二天堂av在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 美女午夜性视频免费| 在线国产一区二区在线| 久久九九热精品免费| 在线观看一区二区三区| 国产亚洲精品综合一区在线观看 | cao死你这个sao货| av天堂在线播放| 日韩大尺度精品在线看网址 | 男人舔女人的私密视频| 亚洲精品国产色婷婷电影|