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

    基于灰色關(guān)聯(lián)分析和K均值聚類的短期負荷預(yù)測

    2021-07-09 09:26:00黃冬梅莊興科胡安鐸孫錦中時帥孫園唐振
    電力建設(shè) 2021年7期
    關(guān)鍵詞:分析法灰色關(guān)聯(lián)

    黃冬梅, 莊興科,胡安鐸,孫錦中,時帥,孫園,唐振

    (1.上海電力大學(xué)電子與信息工程學(xué)院, 上海市 201306;2.上海電力大學(xué)電氣工程學(xué)院,上海市 200090;3.上海電力大學(xué)數(shù)理學(xué)院, 上海市 201306)

    0 引 言

    電力負荷預(yù)測是電力系統(tǒng)安全、穩(wěn)定、經(jīng)濟運行的保障,也是編排調(diào)度計劃、交易計劃的重要基礎(chǔ)[1]。目前智能電網(wǎng)的建設(shè)在我國逐步展開,電力系統(tǒng)中分布式電源的配置也隨之增長,使得負荷變得更加復(fù)雜[2]。另一方面,電動汽車的大量使用,也增大了負荷消耗的隨機性,這些都給短期電力負荷預(yù)測帶來了巨大挑戰(zhàn)。

    短期負荷預(yù)測大致可以分為傳統(tǒng)方法、人工智能方法和相似日方法[3]。傳統(tǒng)方法可分為線性回歸、自回歸移動平均值、灰色模型[4-5]等。這些方法模型簡單、計算速度快,但是由于電力負荷具有非線性的特點,這些傳統(tǒng)方法并不能很好地解決非線性問題。

    近些年,人工智能方法快速發(fā)展,突破了傳統(tǒng)方法不能很好地解決非線性問題的不足。人工智能方法主要可分為人工神經(jīng)網(wǎng)絡(luò)[6]、支持向量機[7-9]、深度學(xué)習(xí)[10]等。隨著人工智能技術(shù)的快速發(fā)展,基于深度學(xué)習(xí)方法的負荷預(yù)測技術(shù)日益受到重視。深度學(xué)習(xí)中循環(huán)神經(jīng)網(wǎng)絡(luò)(recurrent neural network,RNN)[11-12]通過將上一時刻神經(jīng)元狀態(tài)輸入到當(dāng)前時刻神經(jīng)元中,實現(xiàn)了對時序序列的處理。但是它存在梯度爆炸或者梯度消失問題,不能有效地處理長時間序列。長短期記憶(long short-term memory,LSTM)神經(jīng)網(wǎng)絡(luò)在RNN的基礎(chǔ)上做了改進,有效解決了RNN存在的問題,實現(xiàn)了對長時間序列的預(yù)測。文獻[13-15]均采用LSTM神經(jīng)網(wǎng)絡(luò)實現(xiàn)對負荷的預(yù)測,充分發(fā)揮了LSTM神經(jīng)網(wǎng)絡(luò)能夠處理長時間非線性時序序列的特點,取得了較好效果。

    采用神經(jīng)網(wǎng)絡(luò)的負荷預(yù)測精度會受到樣本質(zhì)量的影響。采用相似日的方法可以篩選得到相關(guān)性高的訓(xùn)練樣本,并且能夠有效降低數(shù)據(jù)樣本維度,從而利用較少的訓(xùn)練樣本就能達到較高的預(yù)測精度。但是訓(xùn)練樣本的選取是相似日方法的關(guān)鍵,其相似程度直接決定了預(yù)測算法的精度。文獻[16-17]均采用灰色關(guān)聯(lián)分析(grey relation analysis,GRA)法選取相似日,并結(jié)合各自的神經(jīng)網(wǎng)絡(luò)模型實現(xiàn)負荷預(yù)測。但是灰色關(guān)聯(lián)分析法選取的相似日的負荷曲線與待預(yù)測日的負荷曲線存在較大差異。文獻[18]在灰色關(guān)聯(lián)法基礎(chǔ)上引入加權(quán)和投影概念,但是該改進方法所選取相似日的負荷曲線與待預(yù)測日負荷曲線形狀相似性不佳。文獻[19]結(jié)合灰色關(guān)聯(lián)分析法,提出一種綜合相似度的相似日選取方法并結(jié)合梯度提升決策樹實現(xiàn)負荷預(yù)測,但是該方法未考慮不同因素所衡量相似度的重要性差異,會影響所選取相似日集合的大小。文獻[20]采用聚類的方法選取相似日,并結(jié)合預(yù)測模型完成負荷預(yù)測,但是對聚類數(shù)目的選擇存在主觀性,可能導(dǎo)致所選取的相似日與待預(yù)測日的相似程度不高。

    綜上所述,現(xiàn)有相似日方法中相似日的選取主要是采用灰色關(guān)聯(lián)分析或聚類,但所選取的相似日與待預(yù)測日的相似度有待進一步提升。本文針對該問題,提出一種基于灰色關(guān)聯(lián)分析和K均值聚類(GRA-K-means)選取相似日的短期負荷預(yù)測方法。該方法采用灰色關(guān)聯(lián)分析法選取相似日的粗集,然后建立綜合聚類指標(biāo),確定聚類數(shù)目,并對相似日粗集進行K-means聚類,得到最終的相似日集合,從而提升所選取的相似日與待預(yù)測日的相似度,并且降低樣本維度。最后,采用LSTM神經(jīng)網(wǎng)絡(luò)進行負荷預(yù)測。算例表明,本文提出的采用GRA-K-means選取相似日的短期負荷預(yù)測方法能夠?qū)崿F(xiàn)較高的負荷預(yù)測精度。

    1 相似日選取

    1.1 影響負荷的因素分析

    短期負荷預(yù)測模型的建立需要大量的電力負荷數(shù)據(jù)和其他外部影響因素的歷史數(shù)據(jù)。外部因素中包含了氣象因素和日期類型因素等,負荷會隨外部因素的變化而改變。常見的外部因素有溫度、節(jié)假日等。

    圖1展示了某地區(qū)半年內(nèi)日最高溫度與日總負荷的變化趨勢。

    圖1 某地區(qū)2014年上半年最高溫度與日總負荷變化Fig.1 Variations of the maximum temperature and total daily load in a certain area in the first half of 2014

    由圖1可知,溫度的變化和負荷的變化具有明顯的相關(guān)性,兩者的變化趨勢大致相同,溫度信息可以用來提升負荷預(yù)測的性能??紤]到負荷還受濕度的影響,本文采用溫濕度指數(shù)來描述溫度和濕度之間的耦合關(guān)系,計算公式為[21]:

    ITH=1.8T+0.55(1-H)×(1.8T-26)

    (1)

    式中:ITH為溫濕度指數(shù);T為溫度;H為相對濕度。

    日期類型也是短期負荷預(yù)測中需要考慮的因素,周一到周日標(biāo)記為1到7,若這一天為節(jié)假日時,則特征編碼為1,否則為0。同理可對是否為工作日進行編碼。具體如表1所示。

    表1 日期類型編碼Table 1 Date type encoding

    文獻[3]考慮前日負荷峰值,所選取的相似日比未考慮前日負荷峰值更接近待預(yù)測日,并且所得預(yù)測效果較好,故本文也加入這一因素。本文選取了10個外部因素作為影響因素,基于某地區(qū)2014年1月到6月的數(shù)據(jù)進行Pearson相關(guān)系數(shù)分析,結(jié)果如表2所示。

    表2 相似日影響因素的相關(guān)系數(shù)Table 2 Correlation coefficients of influencing factors on similar days

    從表2中可以看到,前日負荷峰值的相關(guān)性最大,相關(guān)系數(shù)的絕對值為0.898 6,日降雨量的相關(guān)性最低,相關(guān)系數(shù)的絕對值為0.103 1。本文選擇相關(guān)系數(shù)絕對值大于0.3的8個因素作為相似日的影響因素。

    1.2 灰色關(guān)聯(lián)分析法

    灰色關(guān)聯(lián)分析法是分析各因素之間相似程度的一種統(tǒng)計方法。該方法通過一系列的數(shù)學(xué)公式來挖掘數(shù)據(jù)之間的關(guān)系從而分析數(shù)據(jù)之間的相似程度[22]。該方法通過以下步驟,計算各個樣本的灰色關(guān)聯(lián)值。

    1) 選取影響電力負荷的若干因素,則第i天樣本的特征向量可以表示為:

    (2)

    式中:i=1,2,…,n,n為歷史樣本總數(shù);yim為第i個樣本的第m個因素的特征值。待預(yù)測的特征向量表示為:

    (3)

    式中:y0m為待預(yù)測日特征向量的第m個影響因素的特征值。

    2) 構(gòu)建灰色關(guān)聯(lián)判斷矩陣F,如式(4)所示:

    (4)

    各元素按式(5)計算:

    (5)

    式中:Fik(i=1,2,…,n;k=1,2,…,m)為第i個樣本的第k個因素對應(yīng)的關(guān)聯(lián)系數(shù);yik(i=1,2,…,n;k=1,2,…,m)為第i個樣本的第k個因素的特征值;y0k(k=1,2,…,m)為待預(yù)測日特征向量的第k個影響因素的特征值;ρ為分辨系數(shù),通常取ρ=0.5。采用相關(guān)系數(shù)法確定各影響因素的權(quán)重:

    (6)

    (7)

    式中:pk為第k個影響因素的Pearson相關(guān)系數(shù)絕對值;ωk第k個影響因素的權(quán)重,k=1,2,…,m。

    3) 利用上述權(quán)重對灰色關(guān)聯(lián)判斷矩陣進行加權(quán),得到加權(quán)灰色關(guān)聯(lián)決策陣F′:

    (8)

    4) 得到每個樣本的灰色關(guān)聯(lián)值:

    (9)

    式中:Di為第i個樣本的灰色關(guān)聯(lián)度。

    5) 將每個歷史日特征向量的灰色關(guān)聯(lián)度按照數(shù)值從大到小進行排序,設(shè)置閾值,選擇滿足閾值的樣本作為相似日的粗集。本文將閾值設(shè)為0.7。

    1.3 K均值聚類算法

    K均值聚類算法是根據(jù)距離相近的數(shù)據(jù)作為一類的原則,將數(shù)據(jù)樣本進行分類的統(tǒng)計分析方法,優(yōu)點是實現(xiàn)簡單且高效[23]。

    將灰色關(guān)聯(lián)分析法找到的相似日粗集的8個外部影響因素通過K均值聚類,計算待預(yù)測日與聚類中心的歐式距離,計算公式如(10)所示。當(dāng)dj數(shù)值最小時,此時聚類中心所對應(yīng)的一類即為最終相似日集合。

    (10)

    式中:xk(k=1,2,…,m)為待預(yù)測日特征向量的第k個因素的特征值;pjk(j=1,2,…,l;k=1,2,…,m)為第j組聚類中心的第k個因素所對應(yīng)的值。

    2 基于LSTM網(wǎng)絡(luò)的預(yù)測模型

    2.1 LSTM網(wǎng)絡(luò)

    LSTM是一種在循環(huán)神經(jīng)網(wǎng)絡(luò)的基礎(chǔ)上進行改進的神經(jīng)網(wǎng)絡(luò)。它解決了RNN模型存在的問題,并且可以學(xué)習(xí)長時間序列中隱藏信息,充分利用歷史數(shù)據(jù)中的有用信息,具有更強的適應(yīng)性。LSTM網(wǎng)絡(luò)基本單元如圖2所示。

    圖2 LSTM神經(jīng)網(wǎng)絡(luò)基本單元Fig.2 Basic unit of LSTM neural network

    遺忘門、輸入門和輸出門3個部分共同組成了LSTM神經(jīng)網(wǎng)絡(luò)的基本單元[24]。輸入xt和中間輸出ht-1通過sigmoid函數(shù)生成變量ft,并且與上一個狀態(tài)記憶單元St-1決定遺忘部分。輸入xt和中間輸出ht-1在輸入門中經(jīng)過sigmoid函數(shù)和tanh函數(shù)后,再與遺忘門生成的遺忘部分一同決定需要保留的部分。變量St與ot決定了中間輸出ht,計算公式如下所示:

    ft=σ(Wfxxt+Wfhht-1+bf)

    (11)

    it=σ(Wixxt+Wihht-1+bi)

    (12)

    gt=φ(Wgxxt+Wghht-1+bg)

    (13)

    ot=σ(Woxxt+Wohht-1+bo)

    (14)

    St=gt?it+St-1?ft

    (15)

    ht=φ(St)?ot

    (16)

    式中:ft、it、gt、ot、St分別為遺忘門、輸入門、輸入節(jié)點、輸出門和狀態(tài)單元的狀態(tài);Wfx、Wfh、Wfh、Wix、Wih、Wgx、Wgh、Wox和Woh為相應(yīng)門與變量之間的權(quán)重;bf、bi、bg、bo為偏置項;?表示元素按位相乘;σ(·)和φ(·)分別表示sigmoid函數(shù)和tanh函數(shù)。

    2.2 預(yù)測模型

    本文在采用灰色關(guān)聯(lián)分析和K均值聚類選取了最終的相似日集合的基礎(chǔ)上,訓(xùn)練LSTM網(wǎng)絡(luò),實現(xiàn)短期負荷預(yù)測。這一方法的模型如圖3所示,它包括以下幾個步驟:

    圖3 模型框架圖Fig.3 Diagram of the model frame

    步驟1:輸入歷史負荷數(shù)據(jù)和各種外部因素,并對數(shù)據(jù)進行預(yù)處理,選取主要相似日影響因素;

    步驟2:利用灰色關(guān)聯(lián)分析法,找到關(guān)聯(lián)度大于等于閾值的相似日,將這些相似日作為相似日的粗集;

    步驟3:使用K均值聚類對找到的相似日粗集的外部因素進行聚類,找到待預(yù)測日與聚類中心距離最小的一類作為最終相似日集合;

    步驟4:根據(jù)最終相似日集合的日期,將對應(yīng)的歷史負荷數(shù)據(jù)與外部影響因素組合為輸入特征,如表3所示,模型輸入為t時刻的13個特征,輸出為t時刻的負荷值,并將輸入特征和輸出負荷值歸一化之后依次帶入模型訓(xùn)練,最后得到預(yù)測結(jié)果。

    表3 模型的輸入特征Table 3 Input features of the model

    2.3 預(yù)測模型的評價指標(biāo)

    采用平均絕對百分比誤差MAPE和均方根誤差RMSE來評價模型的好壞,公式如下:

    (17)

    (18)

    式中:y′i為預(yù)測值;yi為真實值。

    3 算例分析

    為驗證本文所提方法的科學(xué)性和可靠性,本文實驗選取使用2016年電工數(shù)學(xué)建模競賽中2014年1月1日至2014年6月29日的負荷數(shù)據(jù),數(shù)據(jù)的采樣周期為15 min,共180天,每天共有96組數(shù)據(jù),共17 280條負荷數(shù)據(jù),以2014年1月1日至2014年6月23日的數(shù)據(jù)為歷史負荷數(shù)據(jù),預(yù)測2014年6月24日到6月29日的負荷數(shù)據(jù)。LSTM模型采用基于Python編寫的Keras深度學(xué)習(xí)框架[25],K均值采用常用的sklearn機器學(xué)習(xí)框架。

    3.1 相似日選取有效性分析

    以待預(yù)測日6月24日為例,進行分析。首先使用灰色關(guān)聯(lián)分析法,計算待預(yù)測日與歷史樣本的灰色關(guān)聯(lián)度,將灰色關(guān)聯(lián)度大于等于0.7的歷史樣本作為相似日的粗集,可得到94天相似日。在相似日粗集的基礎(chǔ)上,將每天作為一個數(shù)據(jù)對象,每天特征維度為8,將其輸入到K均值聚類模型中。聚類個數(shù)的選擇是聚類效果重要的一環(huán),傳統(tǒng)聚類個數(shù)確定的方法都是根據(jù)輪廓系數(shù)(Silhouette coefficient, SIL)[26]和戴維森堡丁指數(shù)(Davies-Bouldin Index, DBI)[27]2個指標(biāo)進行決定,需要人為決定聚類個數(shù)。當(dāng)SIL數(shù)值越大,DBI數(shù)值越小時,說明此時聚類的效果最好。圖4為2個聚類指標(biāo)的折線圖,可以看到此時的最佳聚類數(shù)為2,但是SIL和DBI數(shù)值都隨聚類個數(shù)的變化而波動。

    圖4 聚類評價指標(biāo)Fig.4 evaluation index of clustering

    本文將2個聚類指標(biāo)進行統(tǒng)一,提出了綜合聚類指標(biāo),當(dāng)綜合聚類指標(biāo)最小時所對應(yīng)的聚類個數(shù)就是最佳聚類個數(shù),計算公式如下所示:

    (19)

    式中:E為綜合聚類指標(biāo);CSIL和IDB分別為SIL和DBI指標(biāo)的數(shù)值。

    圖5為綜合聚類指標(biāo)隨聚類個數(shù)的變化情況,聚類個數(shù)可根據(jù)綜合聚類指標(biāo)直接決定,不需要進行人為分析決定。故本文將待預(yù)測日(6月24日)的聚類個數(shù)選取為2。計算待預(yù)測日與各聚類中心的歐式距離,結(jié)果如表4所示,故選距離最小一類為最終相似日,共72天。

    圖5 綜合聚類評價指標(biāo)Fig.5 Comprehensive evaluation index of clustering

    表4 聚類中心與待預(yù)測日的歐氏距離Table 4 Euclidean distance between cluster center and the day to be predicted

    為比較本文方法與灰色關(guān)聯(lián)方法所選取的相似日的有效性,本文將形狀系數(shù)[25]作為比較兩者的標(biāo)準。形狀系數(shù)越小,說明兩者曲線在形態(tài)方面越接近。計算公式如下所示:

    xijk=xik-xjk

    (20)

    (21)

    (22)

    式中:xi、xj分別為長度為l的等長序列,分別表示為xi={xi1,xi2,…,xil},xj={xj1,xj2,…,xjl},xik表示第i序列的第k個元素,xjk表示第j序列的第k個元素;Sij為序列i和j的形狀系數(shù)。

    計算兩組相似日負荷與待預(yù)測日之間的形狀系數(shù),其測試集的形狀系數(shù)平均值的計算結(jié)果如表5所示。

    表5 形狀系數(shù)平均值Table 5 Average of shape coefficient

    從表5中可知,本文方法在相似日的選取上,更接近相似日,故本文方法在相似日的選取上比灰色關(guān)聯(lián)分析法效果更好。

    3.2 不同相似日選取方法對比分析

    為驗證本文相似日選取方法是否能提升精度,采用LSTM模型,對不選取相似日(LSTM)、灰色關(guān)聯(lián)選取相似日(GRA-LSTM)和基于灰色關(guān)聯(lián)分析與K均值聚類選取相似日(GRA-K-means-LSTM)3種情況進行對比分析,共同預(yù)測6月24日的負荷。LSTM模型使用2層隱藏層,神經(jīng)元個數(shù)為32 100,最后使用激活函數(shù)為sigmoid的全連接層輸出預(yù)測值,神經(jīng)網(wǎng)絡(luò)采用Adam算法優(yōu)化,學(xué)習(xí)率為0.001。3種情況的誤差如表6所示。

    表6 不同相似日選取方法誤差對比Table 6 Error comparison of different similar-day selection methods

    從表6中可知,采用相似日的方法在MAPE這一項指標(biāo)上比不選取相似日的方法均有降低,傳統(tǒng)的灰色關(guān)聯(lián)分析法選取相似日,在MAPE上相比未選取相似日降低了0.634%,本文的方法在MAPE上降低了1.198%。傳統(tǒng)的灰色關(guān)聯(lián)分析法在RMSE上效果不佳,但是本文方法卻降低了81.124 MW。由此可見,本文提出的采用GRA-K-means選取相似日方法,在誤差精度上均比傳統(tǒng)的灰色關(guān)聯(lián)分析法要好,也證明了方法的有效性。圖6為負荷預(yù)測曲線,從圖6中可知本文方法比未采用相似日方法和采用灰色關(guān)聯(lián)分析方法的LSTM模型更加貼合實際負荷曲線,在峰谷處本文方法比其他2種方法更接近真實值,擬合度更高,預(yù)測誤差更小。

    圖6 不同相似日選取方法的預(yù)測負荷曲線Fig.6 Load forecasting curves of different similar-day search methods

    3.3 不同模型的對比分析

    為驗證本文方法模型的可行性與合理性,采用不選取相似日的LSTM模型和本文所提模型(GRA-K-means-LSTM)、XGBoost模型(GRA-K-means-XGBoost)、LightGBM模型(GRA-K-means-LightGBM)相結(jié)合的模型進行對比分析,預(yù)測6月24日到6月29日共6天的負荷數(shù)據(jù),并對結(jié)果取平均值,結(jié)果如表7所示。XGBoost[28]與LightGBM[29]是各種競賽中常用的算法且取得了比較好的成績,將這兩者當(dāng)作本文的對照模型,來驗證本文方法的有效性。這2種模型均采用sklearn框架。

    表7 不同模型誤差對比Table 7 Error comparison of different models

    如表7所示,本文所提出的GRA-K-means-LSTM方法與其他4種方法相比,在MAPE指標(biāo)上,本文方法比LSTM、GRA-LSTM、GRA-Kmeans-XGBoost、GRA-Kmeans-LightGBM分別降低了0.911%、0.637%、0.503%、0.389%,在RMSE上本文方法比LSTM、GRA-LSTM、GRA-K-means-XGBoost分別降低了55.739、54.579和12.264 MW,與GRA-K-means-LightGBM結(jié)果相差很小。

    在幾種對比模型中,本文所提出的GRA-K-means-LSTM是最優(yōu)的,其原因在于采用灰色關(guān)聯(lián)分析法和K均值聚類的方法選取相似日,刪除了相似度不高的數(shù)據(jù),提高了數(shù)據(jù)的相關(guān)性,并結(jié)合LSTM神經(jīng)網(wǎng)絡(luò)模型,實現(xiàn)對負荷的短期預(yù)測,充分發(fā)揮了LSTM神經(jīng)網(wǎng)絡(luò)能夠處理時間序列和非線性數(shù)據(jù)的特性,提高了預(yù)測精度。

    4 結(jié) 論

    1) 本文選取相似日時加入了前日負荷峰值和溫濕度指數(shù),并使用皮爾遜相關(guān)系數(shù)法證明其與日總負荷值具有明顯的相關(guān)性,可以用于相似日的選取。

    2) 在傳統(tǒng)的灰色關(guān)聯(lián)分析法的基礎(chǔ)上加入了K均值聚類,并通過綜合聚類指標(biāo)選取聚類個數(shù),解決了人工選取聚類中心的缺點。通過改進的方法得到了與待預(yù)測日相似程度更高的相似日,減少了訓(xùn)練樣本的數(shù)量和差異程度,有助于提升負荷預(yù)測精度。

    3) 通過將本文相似日選取方法與LSTM模型相結(jié)合,有效地提高了短期負荷預(yù)測的精度,且LSTM模型能夠更好地擬合負荷數(shù)據(jù)。

    未來可以考慮更多導(dǎo)致電網(wǎng)氣象災(zāi)害的因素,進一步提升短期負荷預(yù)測精度,應(yīng)用于電網(wǎng)氣象災(zāi)害預(yù)警研究。

    猜你喜歡
    分析法灰色關(guān)聯(lián)
    異步機傳統(tǒng)分析法之困難及其克服
    防爆電機(2022年4期)2022-08-17 05:59:06
    淺灰色的小豬
    “一帶一路”遞進,關(guān)聯(lián)民生更緊
    奇趣搭配
    基于時間重疊分析法的同車倒卡逃費探析
    灰色時代
    Coco薇(2017年2期)2017-04-25 17:59:38
    她、它的灰色時髦觀
    Coco薇(2017年2期)2017-04-25 17:57:49
    智趣
    讀者(2017年5期)2017-02-15 18:04:18
    層次分析法在SWOT分析法中的應(yīng)用
    感覺
    日本a在线网址| 性色av乱码一区二区三区2| 在线a可以看的网站| 国产毛片a区久久久久| 国产精品乱码一区二三区的特点| 老司机深夜福利视频在线观看| 在线观看舔阴道视频| 真人一进一出gif抽搐免费| 一个人免费在线观看电影| 69人妻影院| 精品免费久久久久久久清纯| 欧美一区二区国产精品久久精品| 亚洲成人精品中文字幕电影| 亚洲 国产 在线| 国产真实伦视频高清在线观看 | 成人国产综合亚洲| 日本在线视频免费播放| 首页视频小说图片口味搜索| 别揉我奶头~嗯~啊~动态视频| 一区二区三区四区激情视频 | 精品人妻1区二区| 一进一出好大好爽视频| 亚洲最大成人av| 好男人电影高清在线观看| 国产熟女xx| 人妻久久中文字幕网| 日日干狠狠操夜夜爽| 欧美精品国产亚洲| 日韩欧美国产一区二区入口| 黄色一级大片看看| 国内精品久久久久久久电影| 乱人视频在线观看| 丁香欧美五月| 久久精品91蜜桃| 91在线观看av| 99久久成人亚洲精品观看| 中文字幕精品亚洲无线码一区| 日韩av在线大香蕉| 久久这里只有精品中国| av在线蜜桃| 老熟妇乱子伦视频在线观看| 午夜免费男女啪啪视频观看 | 可以在线观看的亚洲视频| 免费大片18禁| 黄片小视频在线播放| av在线观看视频网站免费| 久久国产精品影院| 日日摸夜夜添夜夜添av毛片 | 少妇丰满av| 国产白丝娇喘喷水9色精品| 午夜免费男女啪啪视频观看 | 久久草成人影院| 久久精品综合一区二区三区| 宅男免费午夜| 亚洲精品成人久久久久久| 亚洲中文字幕一区二区三区有码在线看| 熟妇人妻久久中文字幕3abv| 国产伦精品一区二区三区四那| 51国产日韩欧美| 欧美成人一区二区免费高清观看| 欧美xxxx黑人xx丫x性爽| 亚洲第一电影网av| 免费在线观看成人毛片| 免费在线观看亚洲国产| 最后的刺客免费高清国语| 岛国在线免费视频观看| 国产精品1区2区在线观看.| 成年免费大片在线观看| 国产精华一区二区三区| 国内毛片毛片毛片毛片毛片| 宅男免费午夜| 国产免费av片在线观看野外av| 国产男靠女视频免费网站| 久久久久亚洲av毛片大全| 国产精品综合久久久久久久免费| 成人特级黄色片久久久久久久| 国产探花极品一区二区| 久久午夜亚洲精品久久| 一级毛片久久久久久久久女| 欧美黄色片欧美黄色片| 男女做爰动态图高潮gif福利片| 日韩成人在线观看一区二区三区| 我的老师免费观看完整版| 欧美区成人在线视频| 国产高清视频在线观看网站| 精品久久久久久,| 最好的美女福利视频网| 午夜福利在线在线| 最好的美女福利视频网| 中文字幕人成人乱码亚洲影| 男女做爰动态图高潮gif福利片| 国产免费av片在线观看野外av| 麻豆成人av在线观看| 亚洲国产精品久久男人天堂| 最新在线观看一区二区三区| 又爽又黄无遮挡网站| 国内精品久久久久久久电影| 老司机午夜十八禁免费视频| 亚洲人成伊人成综合网2020| 首页视频小说图片口味搜索| 亚洲国产精品久久男人天堂| 十八禁国产超污无遮挡网站| 天堂av国产一区二区熟女人妻| 亚洲av五月六月丁香网| 少妇高潮的动态图| 757午夜福利合集在线观看| 亚洲欧美日韩东京热| 级片在线观看| 亚洲天堂国产精品一区在线| 欧美3d第一页| 高清在线国产一区| 国产激情偷乱视频一区二区| 一本综合久久免费| 亚洲欧美清纯卡通| or卡值多少钱| 国产不卡一卡二| 国产91精品成人一区二区三区| 美女黄网站色视频| 亚洲av电影在线进入| 国产三级中文精品| 欧美又色又爽又黄视频| 美女大奶头视频| 精品国产三级普通话版| 精华霜和精华液先用哪个| 欧美黄色片欧美黄色片| 天堂动漫精品| a级毛片免费高清观看在线播放| 在线免费观看的www视频| 久久人人爽人人爽人人片va | 亚洲激情在线av| www.999成人在线观看| or卡值多少钱| 亚洲aⅴ乱码一区二区在线播放| 99国产精品一区二区蜜桃av| 69av精品久久久久久| 青草久久国产| 在线免费观看不下载黄p国产 | 亚洲avbb在线观看| av黄色大香蕉| 精品一区二区三区av网在线观看| 好看av亚洲va欧美ⅴa在| 又紧又爽又黄一区二区| 91字幕亚洲| 亚洲精品色激情综合| 国产精品自产拍在线观看55亚洲| 国产精品久久视频播放| 亚洲人成电影免费在线| 亚洲一区高清亚洲精品| 日本撒尿小便嘘嘘汇集6| 亚洲自偷自拍三级| 综合色av麻豆| 成人毛片a级毛片在线播放| 久久午夜亚洲精品久久| 男女视频在线观看网站免费| av在线观看视频网站免费| 中文字幕av成人在线电影| 欧美区成人在线视频| 欧美乱色亚洲激情| 国产黄片美女视频| 此物有八面人人有两片| www.色视频.com| 国产精品自产拍在线观看55亚洲| 亚洲五月天丁香| 午夜激情福利司机影院| 久久国产精品人妻蜜桃| 窝窝影院91人妻| 亚洲精品一卡2卡三卡4卡5卡| 日本与韩国留学比较| 亚洲激情在线av| 免费观看人在逋| 1024手机看黄色片| 午夜福利成人在线免费观看| 啦啦啦观看免费观看视频高清| 亚洲欧美日韩东京热| 麻豆久久精品国产亚洲av| 久久精品国产自在天天线| 亚洲av日韩精品久久久久久密| 午夜免费激情av| 国产亚洲欧美98| 免费在线观看日本一区| 网址你懂的国产日韩在线| 此物有八面人人有两片| 免费电影在线观看免费观看| 久久久精品欧美日韩精品| 亚洲av免费高清在线观看| 欧美日韩亚洲国产一区二区在线观看| 99精品在免费线老司机午夜| 亚洲精品一卡2卡三卡4卡5卡| 久久人人精品亚洲av| 国产欧美日韩精品亚洲av| 国模一区二区三区四区视频| 久久久久性生活片| 国产精品一区二区免费欧美| 国产成+人综合+亚洲专区| 婷婷精品国产亚洲av在线| 国产免费av片在线观看野外av| 一个人免费在线观看电影| 真人一进一出gif抽搐免费| 国产成人啪精品午夜网站| 国产在线男女| 99久久九九国产精品国产免费| 国产精品久久久久久亚洲av鲁大| 男女做爰动态图高潮gif福利片| 99国产综合亚洲精品| 真人一进一出gif抽搐免费| 校园春色视频在线观看| 久久久久久久久中文| 一级黄色大片毛片| 亚洲精品成人久久久久久| 欧美一区二区国产精品久久精品| 中亚洲国语对白在线视频| 国产欧美日韩一区二区精品| 亚洲avbb在线观看| 高清在线国产一区| 91九色精品人成在线观看| 国内精品久久久久精免费| 伊人久久精品亚洲午夜| ponron亚洲| 热99在线观看视频| 我要看日韩黄色一级片| 一边摸一边抽搐一进一小说| 97超级碰碰碰精品色视频在线观看| 十八禁国产超污无遮挡网站| 欧美精品啪啪一区二区三区| 久久欧美精品欧美久久欧美| 亚洲精品色激情综合| 性色av乱码一区二区三区2| 久久午夜福利片| 免费在线观看亚洲国产| 国产色爽女视频免费观看| 亚洲欧美日韩东京热| 91午夜精品亚洲一区二区三区 | 高潮久久久久久久久久久不卡| 99精品在免费线老司机午夜| 女同久久另类99精品国产91| 少妇人妻一区二区三区视频| 亚洲熟妇熟女久久| 色精品久久人妻99蜜桃| 男人狂女人下面高潮的视频| 免费观看精品视频网站| 成年版毛片免费区| 日本五十路高清| 97超视频在线观看视频| 亚洲精品影视一区二区三区av| 精品日产1卡2卡| 久久久久久久久中文| 色哟哟哟哟哟哟| 九九久久精品国产亚洲av麻豆| av视频在线观看入口| 直男gayav资源| 91在线观看av| 亚洲精品一区av在线观看| 国产私拍福利视频在线观看| 国产精品自产拍在线观看55亚洲| 欧美成狂野欧美在线观看| 精品久久久久久久久久久久久| 一本综合久久免费| 久久精品人妻少妇| 高清日韩中文字幕在线| 99久久无色码亚洲精品果冻| 亚洲精品456在线播放app | 日本三级黄在线观看| 啪啪无遮挡十八禁网站| 欧美3d第一页| 日韩欧美国产在线观看| 日本a在线网址| 精品人妻熟女av久视频| 亚洲av成人精品一区久久| 日本成人三级电影网站| 色综合亚洲欧美另类图片| 色播亚洲综合网| 亚洲精品456在线播放app | 女人被狂操c到高潮| 国产亚洲精品综合一区在线观看| 久久精品91蜜桃| 男女做爰动态图高潮gif福利片| 亚洲aⅴ乱码一区二区在线播放| x7x7x7水蜜桃| 波多野结衣巨乳人妻| 麻豆久久精品国产亚洲av| 久久精品夜夜夜夜夜久久蜜豆| 午夜亚洲福利在线播放| 热99re8久久精品国产| 亚洲在线观看片| 国产精品一区二区三区四区免费观看 | 亚洲精品乱码久久久v下载方式| 国产精品国产高清国产av| 亚洲内射少妇av| 两个人的视频大全免费| 国产蜜桃级精品一区二区三区| 啪啪无遮挡十八禁网站| 欧美日韩亚洲国产一区二区在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产午夜精品久久久久久一区二区三区 | 亚洲精品色激情综合| 97热精品久久久久久| 亚洲人成网站在线播| 亚洲经典国产精华液单 | 中出人妻视频一区二区| 成人精品一区二区免费| 欧美一区二区精品小视频在线| 无遮挡黄片免费观看| 两个人的视频大全免费| 超碰av人人做人人爽久久| 久久这里只有精品中国| 亚洲av免费在线观看| ponron亚洲| 久久精品国产亚洲av涩爱 | 久久久精品大字幕| 在线观看舔阴道视频| 老熟妇仑乱视频hdxx| 少妇的逼好多水| 老司机福利观看| 男人舔奶头视频| 久久久久久久久大av| 青草久久国产| 国产精品亚洲一级av第二区| 99久久无色码亚洲精品果冻| 赤兔流量卡办理| 国产欧美日韩一区二区三| 在线看三级毛片| 日本黄色视频三级网站网址| 亚洲av熟女| 亚洲中文日韩欧美视频| 国产黄a三级三级三级人| 亚洲av电影不卡..在线观看| 少妇人妻一区二区三区视频| 成年免费大片在线观看| 日本免费一区二区三区高清不卡| 美女cb高潮喷水在线观看| 免费一级毛片在线播放高清视频| 亚洲,欧美精品.| 99久国产av精品| 国产高潮美女av| 18美女黄网站色大片免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 俄罗斯特黄特色一大片| 亚洲中文日韩欧美视频| 午夜免费成人在线视频| 日韩欧美国产一区二区入口| 99精品在免费线老司机午夜| 欧美高清成人免费视频www| 天堂av国产一区二区熟女人妻| 日日夜夜操网爽| 国产亚洲精品av在线| netflix在线观看网站| 欧美乱色亚洲激情| 国产精品亚洲一级av第二区| 身体一侧抽搐| 人妻夜夜爽99麻豆av| 一区福利在线观看| 国产精品av视频在线免费观看| 久久精品久久久久久噜噜老黄 | 国产美女午夜福利| 中文亚洲av片在线观看爽| 亚洲国产精品成人综合色| 成人永久免费在线观看视频| av在线天堂中文字幕| 国产欧美日韩精品一区二区| 国产国拍精品亚洲av在线观看| 精品久久久久久久久久久久久| 欧美在线黄色| 特大巨黑吊av在线直播| 18禁黄网站禁片免费观看直播| 成年免费大片在线观看| 日韩欧美精品v在线| 亚洲欧美激情综合另类| 国产男靠女视频免费网站| 国产美女午夜福利| 午夜福利在线在线| 在线看三级毛片| 神马国产精品三级电影在线观看| 三级毛片av免费| 亚洲成人久久性| 成年人黄色毛片网站| 在线观看免费视频日本深夜| 性色avwww在线观看| 国产精品1区2区在线观看.| 国产高清有码在线观看视频| 看十八女毛片水多多多| 性插视频无遮挡在线免费观看| 757午夜福利合集在线观看| 在线观看美女被高潮喷水网站 | 一二三四社区在线视频社区8| 国产综合懂色| 91麻豆精品激情在线观看国产| 亚洲精品乱码久久久v下载方式| 一级a爱片免费观看的视频| 首页视频小说图片口味搜索| 日韩成人在线观看一区二区三区| 国产精品伦人一区二区| 婷婷色综合大香蕉| 亚洲人成网站在线播| 精品久久久久久久久亚洲 | 搞女人的毛片| 嫩草影院新地址| 内射极品少妇av片p| 搡老妇女老女人老熟妇| 欧美日韩瑟瑟在线播放| 99久久无色码亚洲精品果冻| 国产精品亚洲一级av第二区| 男人和女人高潮做爰伦理| 日本黄色视频三级网站网址| 国产真实乱freesex| 国产乱人视频| aaaaa片日本免费| 人人妻,人人澡人人爽秒播| 非洲黑人性xxxx精品又粗又长| 在线观看av片永久免费下载| 中文字幕精品亚洲无线码一区| 日韩欧美国产一区二区入口| 国产淫片久久久久久久久 | 午夜福利在线观看吧| 亚洲精品一卡2卡三卡4卡5卡| 日本 欧美在线| 91狼人影院| 亚洲av一区综合| 内射极品少妇av片p| 欧美成人性av电影在线观看| 欧美绝顶高潮抽搐喷水| 亚洲av熟女| 国产黄片美女视频| 十八禁人妻一区二区| 亚洲国产高清在线一区二区三| 欧美成人性av电影在线观看| 欧美日韩综合久久久久久 | 午夜福利高清视频| 十八禁网站免费在线| 日韩欧美免费精品| 亚洲欧美激情综合另类| 久久精品国产亚洲av涩爱 | 99riav亚洲国产免费| 日本三级黄在线观看| 国产老妇女一区| 成人无遮挡网站| 欧美国产日韩亚洲一区| 桃色一区二区三区在线观看| 亚洲av免费高清在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲最大成人中文| 午夜福利高清视频| 欧美性猛交黑人性爽| 婷婷丁香在线五月| 精品久久久久久久人妻蜜臀av| 欧美激情国产日韩精品一区| 国产午夜精品久久久久久一区二区三区 | 欧美zozozo另类| 狂野欧美白嫩少妇大欣赏| 可以在线观看毛片的网站| 亚洲精品在线美女| 亚洲av免费在线观看| 国产亚洲精品久久久com| 午夜视频国产福利| 最好的美女福利视频网| 岛国在线免费视频观看| 99久国产av精品| 99视频精品全部免费 在线| 欧美黄色片欧美黄色片| 国产成人aa在线观看| 国产午夜精品论理片| 久久精品久久久久久噜噜老黄 | 亚洲国产精品sss在线观看| 亚洲无线在线观看| 午夜免费激情av| 国产v大片淫在线免费观看| 国产熟女xx| 男女床上黄色一级片免费看| 一级av片app| 日韩欧美 国产精品| 国产精品久久视频播放| 一个人看视频在线观看www免费| 亚洲内射少妇av| 九色国产91popny在线| 欧美bdsm另类| 老司机午夜十八禁免费视频| 成人特级黄色片久久久久久久| 亚洲aⅴ乱码一区二区在线播放| 久久草成人影院| 看十八女毛片水多多多| 亚洲av成人av| 久久性视频一级片| 99久久99久久久精品蜜桃| 深夜a级毛片| 长腿黑丝高跟| 最近最新免费中文字幕在线| 一卡2卡三卡四卡精品乱码亚洲| 日本熟妇午夜| 国产免费av片在线观看野外av| 禁无遮挡网站| 亚洲色图av天堂| а√天堂www在线а√下载| 在线观看舔阴道视频| 嫩草影院精品99| 成人精品一区二区免费| 久久人人爽人人爽人人片va | 国产一区二区三区视频了| 国产亚洲精品久久久com| 午夜福利在线观看免费完整高清在 | 一个人免费在线观看电影| 免费av观看视频| 日本a在线网址| 久久香蕉精品热| 18禁黄网站禁片午夜丰满| 午夜福利欧美成人| 国产 一区 欧美 日韩| 欧美一区二区国产精品久久精品| 国产精品乱码一区二三区的特点| 成年人黄色毛片网站| 午夜福利在线在线| 精品久久久久久久久久久久久| 久久精品综合一区二区三区| 757午夜福利合集在线观看| 久久午夜亚洲精品久久| 亚洲欧美日韩无卡精品| 亚洲人成网站高清观看| 制服丝袜大香蕉在线| 麻豆成人av在线观看| 国产美女午夜福利| 草草在线视频免费看| 免费一级毛片在线播放高清视频| 欧美在线一区亚洲| 亚洲专区中文字幕在线| 亚洲美女搞黄在线观看 | 日韩欧美三级三区| 中文资源天堂在线| 18美女黄网站色大片免费观看| 91午夜精品亚洲一区二区三区 | 亚洲经典国产精华液单 | 国内揄拍国产精品人妻在线| 丁香六月欧美| 亚洲aⅴ乱码一区二区在线播放| 悠悠久久av| 免费在线观看成人毛片| 自拍偷自拍亚洲精品老妇| 一本久久中文字幕| 一级作爱视频免费观看| 1000部很黄的大片| 亚洲五月婷婷丁香| 亚洲精品在线美女| 成人性生交大片免费视频hd| 波多野结衣高清无吗| 在线免费观看不下载黄p国产 | 久久精品人妻少妇| 简卡轻食公司| 91午夜精品亚洲一区二区三区 | 久久伊人香网站| 日韩欧美三级三区| 在线观看av片永久免费下载| 在线免费观看不下载黄p国产 | 成人三级黄色视频| 国产精品亚洲av一区麻豆| 2021天堂中文幕一二区在线观| 麻豆av噜噜一区二区三区| 日韩成人在线观看一区二区三区| 精品午夜福利视频在线观看一区| 一区福利在线观看| 日韩欧美国产在线观看| 国产精品一区二区三区四区久久| 中文字幕久久专区| 久久久久久久午夜电影| 欧美+日韩+精品| 男女视频在线观看网站免费| 国产黄片美女视频| 久久性视频一级片| 日韩人妻高清精品专区| 日本黄色片子视频| 亚洲精品日韩av片在线观看| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美一区二区三区在线观看| 又爽又黄a免费视频| 免费电影在线观看免费观看| 亚洲人成网站高清观看| 国产三级黄色录像| 国产午夜精品论理片| 日日干狠狠操夜夜爽| 99视频精品全部免费 在线| 首页视频小说图片口味搜索| 哪里可以看免费的av片| 国产精品98久久久久久宅男小说| 非洲黑人性xxxx精品又粗又长| 我的女老师完整版在线观看| 欧美日韩国产亚洲二区| 精品一区二区三区视频在线| 亚洲成人中文字幕在线播放| 真人做人爱边吃奶动态| 精品久久久久久久末码| 村上凉子中文字幕在线| 国产在线精品亚洲第一网站| 嫩草影院新地址| 亚洲专区中文字幕在线| 中文在线观看免费www的网站| 一夜夜www| 国产高潮美女av| 黄色丝袜av网址大全| 观看美女的网站| av在线蜜桃| 一卡2卡三卡四卡精品乱码亚洲| 露出奶头的视频| 999久久久精品免费观看国产| 精品99又大又爽又粗少妇毛片 | 夜夜夜夜夜久久久久| 久久人人爽人人爽人人片va | 亚洲七黄色美女视频| 97碰自拍视频| 亚洲av成人精品一区久久| 日韩av在线大香蕉| 国产男靠女视频免费网站| 亚洲三级黄色毛片| 亚洲国产精品成人综合色| 免费看a级黄色片| 国产精品电影一区二区三区| 亚洲精品影视一区二区三区av| 国产午夜精品久久久久久一区二区三区 | 亚洲av成人av| 变态另类成人亚洲欧美熟女|