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

    基于注意力機制的CNN-BIGRU 短期電價預測

    2024-04-17 09:16:36楊超冉啟武羅德虎豆旺
    關(guān)鍵詞:特征提取機制模型

    楊超,冉啟武,羅德虎,豆旺

    (陜西理工大學電氣工程學院,漢中 723001)

    隨著我國電力市場體系穩(wěn)步有序推進,多元化的競爭格局已初具雛形,市場化發(fā)售電量已大幅度增加,在資源優(yōu)化配置中市場占據(jù)的作用日益增大。電價作為電力市場中不可或缺的重要因素,對電價的準確預測不僅可以對市場參與者提供指導,還可以了解電力市場的運轉(zhuǎn)狀況,為電力監(jiān)管部門提供客觀依據(jù)[1-2]。因此,電力系統(tǒng)中電價預測對電力系統(tǒng)的正常運行起到重要的作用。

    目前,國內(nèi)外對電力市場中電價預測的方法主要有統(tǒng)計學方法和機器學習方法。統(tǒng)計學方法主要包括差分整合移動平均自回歸ARIMA(autoregressive integrated moving average)模型[3]、廣義自回歸條件異方差GARCH(generalized auto regressive conditional heteroscedasticity)[4]等。但是傳統(tǒng)的統(tǒng)計學方法在處理非線性問題或波動較大的時間序列時性能表現(xiàn)較差。隨著人工智能和大數(shù)據(jù)技術(shù)的發(fā)展,機器學習方法被大量應(yīng)用在電價預測的研究中。機器學習方法主要包括支持向量機SVM(support vector machine)[5]、隨機森林RF(random forest)[6]、循環(huán)神 經(jīng)網(wǎng)絡(luò)RNN(recurrent neural network)[7]、長短期記憶LSTM(long short-term memory)網(wǎng)絡(luò)[8]、門控循環(huán)單元GRU(gated recurrent unit)網(wǎng)絡(luò)[9]等。但是單一的機器學習模型往往由于模型的局限性,難以反映數(shù)據(jù)的變化趨勢,無法對電價進行更加準確的預測,所以國內(nèi)外學者往往通過將多種模型進行組合來對電價進行預測。組合多種模型可以充分發(fā)揮單一模型優(yōu)勢,故組合預測模型預測準確性往往高于單一預測模型[10]。

    卷積神經(jīng)網(wǎng)絡(luò)CNN(convolutional neural network)被廣泛應(yīng)用于自然語言處理和計算機視覺,但在處理非線性電價序列時表現(xiàn)欠佳,故常用CNN模型和RNN 模型結(jié)合來對時間序列進行預測。但是鑒于數(shù)據(jù)量龐大時,使用RNN 會存在梯度消失或梯度爆炸等問題,故Hochreiter等[11]提出了LSTM以解決上述問題。文獻[12]通過引入風荷載比特征,利用LSTM對丹麥電力市場進行預測,雖然引入了風荷載比特征可以提高模型的精確度,但是在處理影響因素過多、特征量較大的電價序列時,單一的LSTM 模型就會產(chǎn)生過擬合現(xiàn)象而影響精確度。為提高預測電價模型的準確性,文獻[13]將CNN 和LSTM結(jié)合,提出了一種CNN-LSTM預測模型,先利用CNN 提取特征向量,再將特征輸入進LSTM 進行訓練并預測,得到最終預測結(jié)果,有效地避免了上述過擬合的問題。但是因為LSTM 參數(shù)較多、訓練時間較長,故Bahdanan[14]在2014 年提出了GRU。GRU 是LSTM 的一種變形,GRU 結(jié)構(gòu)相比LSTM 做了簡化,訓練速度相對更快。文獻[15]利用GRU 模型并引入注意力機制以突出關(guān)鍵性的輸入特征,在計負荷和不計負荷兩種情況下對多種算法模型分別進行分析,驗證了該模型在預測精度和訓練時間上都優(yōu)于LSTM,但是單向GRU 捕獲序列中的信息時能力弱于雙向門控循環(huán)單元BIGRU(bidirectional gated recurrent unit)。文獻[16-17]都采用了BIGRU,先利用CNN進行特征提取,然后將特征向量構(gòu)建成時序并輸入到BIGRU或BIGRU-全連接神經(jīng)網(wǎng)絡(luò)模型中訓練,得到最終預測結(jié)果,通過對比實驗驗證了BIGRU比GRU更具優(yōu)勢。

    綜上所述,本文提出了一種基于注意力機制的CNN-BIGRU模型,對西班牙電力市場中電價進行預測。該模型通過對電價影響因素的引入,例如上周期的電價、負荷、能源價格、天氣、溫度等,訓練大量歷史數(shù)據(jù),提取和映射輸入輸出間的內(nèi)在關(guān)系[18],最終得到精確的預測結(jié)果。在本文預測電價實驗中,首先對原始電價數(shù)據(jù)進行數(shù)據(jù)預處理,利用由一維卷積層和池化層搭建的CNN進行特征提??;其次,將CNN提取的特征向量構(gòu)建成時間序列并輸入到BIGRU,利用BIGRU 在處理時間序列方面的優(yōu)勢進行數(shù)據(jù)訓練;然后,尋找輸入數(shù)據(jù)中與預測量之間有強相關(guān)的變量并賦予合理的權(quán)值,采用注意力機制對BIGRU每個時間步的輸出進行加權(quán)求和;最后在全連接層利用激活函數(shù)計算輸出最終預測值。將本文所提模型與其他多種模型進行比較,實驗結(jié)果驗證了本文模型具有更好的精確度和穩(wěn)定性。

    1 深度學習模型介紹

    1.1 CNN 模型

    CNN借助卷積運算的特性,可以對數(shù)據(jù)進行特征提取[19-20],并且模型采用權(quán)值共享能大幅度減少訓練時間,降低錯誤率。CNN 主要由卷積層、池化層和全連接層組成,其結(jié)構(gòu)如圖1所示。

    圖1 CNN 結(jié)構(gòu)Fig.1 Structure of CNN

    卷積層是CNN的核心,主要是對數(shù)據(jù)進行特征提取,以實現(xiàn)局部連接和權(quán)值共享,其核心是進行卷積運算。池化層的主要作用是減小卷積核的尺寸,減少參數(shù)以防止過擬合,對運算進行降維處理,主要有平均池化、最大池化和隨機池化。全連接層是將經(jīng)過卷積和池化后的數(shù)據(jù)映射在CNN 的隱層特征空間,以進行維度轉(zhuǎn)換,最終輸出一維向量。全連接層主要采用sigmoid激活函數(shù)。

    1.2 BIGRU 模型

    LSTM 解決了RNN 無法處理長期依賴問題,通過引入自循環(huán),解決了梯度消失的問題。相比于RNN,LSTM主要在其隱含層做出了改進,通過門機制改變隱含層中的神經(jīng)元。GRU 作為LSTM 的一種變形,在結(jié)構(gòu)上進行了簡化,將LSTM的三個門機制優(yōu)化成GRU 的兩個門機制,減少了訓練時間,同時保證了預測精度[21]。

    GRU 模型包括更新門和重置門,其結(jié)構(gòu)如圖2所示。更新門用于控制前一時刻的狀態(tài)信息保留到當前狀態(tài)的程度,重置門決定是否要結(jié)合當前狀態(tài)與先前的信息[22]。更新門和重置門的輸入都是由當前時間步xt和上一個隱藏狀態(tài)時間步ht-1共同決定,輸出由全連接層通過sigmoid 函數(shù)計算得到。重置門有助于捕獲時間序列中的短期依賴關(guān)系,更新門有助于捕獲時間序列中的長期依賴關(guān)系。GRU內(nèi)部的數(shù)學關(guān)系可表示為

    圖2 GRU 結(jié)構(gòu)Fig.2 Structure of GRU

    式中:xt為輸入;zt為更新門;rt為重置門;ht-1為之前的隱藏層狀態(tài);ht為此層的輸出;為當前候選集狀態(tài);σ 為sigmoid函數(shù);tanh為雙曲正切函數(shù);W(z)、W(r)、W、U(z)、U(r)、U為訓練參數(shù)矩陣。

    BIGRU由一個正向GRU和一個反向GRU共同組成,使信息可以雙向傳播。BIGRU的輸出受前向狀態(tài)和后向狀態(tài)共同影響,相比單向GRU更能反映過去和未來時刻對當前的影響[23]。BIGRU 結(jié)構(gòu)如圖3所示,其中xt-1、xt、xt+1為輸入,ht-1、ht、ht+1為由正向GRU和反向GRU共同作用的隱藏狀態(tài)。

    圖3 BIGRU 結(jié)構(gòu)Fig.3 Structure of BIGRU

    1.3 注意力機制

    注意力機制是通過模擬人腦對外界信息注意力的分配而提出的一種資源分配機制。人腦的注意力會根據(jù)對某些事物的關(guān)注程度而產(chǎn)生不同的注意程度,在需要被注意力高度集中的信息分配更多注意力,在不重要的信息上產(chǎn)生更少的注意力,以此讓人腦獲取更多有用的信息,減少無用信息,實現(xiàn)資源合理分配利用[14]。

    注意力機制對重要強相關(guān)的數(shù)據(jù)賦予較高的權(quán)重值,以突出強相關(guān)因素對實驗結(jié)果的影響;對弱相關(guān)或者無關(guān)的數(shù)據(jù)賦予較小的權(quán)重值甚至忽略,以減少無關(guān)信息對模型的影響,提高模型的準確性。注意力機制是由傳統(tǒng)編碼-解碼框架衍生而來。但傳統(tǒng)的編碼-解碼模型具有局限性,即在編碼和解碼之間的唯一聯(lián)系是固定長度的語義向量。注意力機制結(jié)構(gòu)如圖4所示。

    圖4 注意力機制結(jié)構(gòu)Fig.4 Structure of Attention mechanism

    圖4中,x1、x2、x3、xn為輸入值,為經(jīng)過BIGRU 運算后的輸出值,at,1、at,2、at,3、at,n為注意力機制分配概率分布值,yt、yn為最終的輸出值。

    2 基于注意力機制的CNN-BIGRU 預測模型

    2.1 預測模型敘述

    實驗通過收集西班牙電力市場中歷史電價數(shù)據(jù),對歷史電價序列進行特征提取來學習其潛在的變化規(guī)律并預測出未來電價的變化趨勢。對于龐大的數(shù)據(jù)集,人工提取歷史電價中特征的工程量過大,并且人工特征提取會破壞歷史電價的時序性和潛在規(guī)律,無法準確提取重要信息作為特征值。因此,本文借助CNN的優(yōu)勢,充分利用CNN的一維卷積層優(yōu)勢對原始數(shù)據(jù)進行特征提取,不僅節(jié)省了時間,還提高了模型預測的準確性。然后將提取的特征向量構(gòu)建成時間序列輸入至BIGRU,利用GRU在處理時間序列方向的優(yōu)勢建模訓練,發(fā)掘歷史電價序列中的變化規(guī)律。由于影響電價變化的因素太多,故引入注意力機制,通過對特征重要程度的大小來賦予不同的權(quán)重值,保留或突出重要信息,減少或者丟棄無關(guān)信息,以便提高預測模型的精確度。采用注意力機制對BIGRU 每個時間步的輸出進行加權(quán)求和。最后在全連接層通過sigmoid 激活函數(shù)計算輸出最終預測值。

    2.2 預測模型結(jié)構(gòu)

    基于注意力機制的CNN-BIGRU預測電價的模型結(jié)構(gòu)如圖5所示。

    圖5 基于注意力機制的CNN-BIGRU 模型結(jié)構(gòu)Fig.5 Structure of CNN-BIGRU model based on Attention mechanism

    具體算法步驟如下。

    步驟1輸入層對原始數(shù)據(jù)進行預處理,對數(shù)據(jù)中的異常值和空缺值處理后,對其余數(shù)據(jù)進行歸一化,將處理后的數(shù)據(jù)作為預測模型的輸入。設(shè)序列長度等于n的輸入數(shù)據(jù)為X=[x1,x2,x3,…,xn]T。

    步驟2CNN 層對電價序列進行特征提取,首先利用一維卷積層(Conv1D)對歷史電價序列的特征進行提?。蝗缓罄米畲蟪鼗姆绞竭^濾信息,防止過擬合風險;最后通過sigmoid 激活函數(shù)得到CNN 層的輸出。設(shè)CNN 層的輸出長度為i,則其輸出Hc=[Hc1,…,Hci-1,Hci]T可表示為

    卷積層輸出C可表示為

    其中

    式中:?表示卷積運算;P為池化層的輸出;W1、W2為權(quán)重值;b、b1為偏置向量。

    步驟3BIGRU 層對經(jīng)過CNN 提取的特征進行訓練。GRU的具體計算方法見第1.2節(jié)。BIGRU層在t時刻的輸出ht可表示為

    式中,BIGRU()表示經(jīng)過BIGRU層運算。

    步驟4注意力機制層的輸入是經(jīng)過BIGRU層計算得到的輸出向量ht,通過注意力機制可以得到不同輸入特征向量的權(quán)重。注意力機制的權(quán)重計算方法如下。

    相似度計算公式為

    權(quán)值歸一化計算公式為

    加權(quán)求和計算公式為

    式中:V、W3、U1為權(quán)重參數(shù);et為t時刻注意力機制的概率分布值;at為t時刻BIGRU隱藏層輸出值對應(yīng)注意力機制的權(quán)重;st為t時刻注意力機制輸出值。

    步驟5全連接層通過sigmoid 激活函數(shù)計算出預測值。設(shè)t時刻的預測值為yt,可表示為

    式中:yt為t時刻預測值;w0為權(quán)重矩陣;b0為偏置向量。

    步驟6輸出層是將本文模型訓練所得的數(shù)據(jù)進行反歸一化并最終輸出預測值Ypredict。設(shè)預測步長為m,則Ypredict=[y1,y2,…,ym]T。

    2.3 損失函數(shù)

    損失函數(shù)用來判斷模型的預測值和真實值之間的差距。模型的性能通??梢酝ㄟ^損失函數(shù)來評價,損失函數(shù)值越小則模型性能越好。本文基于注意力機制的CNN-BIGRU預測模型的損失函數(shù)選用均方誤差MSE(mean squared error),其計算公式為

    式中:Loss為損失函數(shù);n為樣本個數(shù);yi為第i個采樣點的實際值;為第i個采樣點的預測值。

    3 算例分析

    本文實驗的硬件平臺使用Intel(R)Xeon(R)Platinum 8358P 處理器 和NVIDIA RTX 3080顯卡。本文所提出基于注意力機制的CNN-BIGRU模型由Tensorflow2.9、Keras、Scikit-learn 等框架搭建,使用python3.8編寫。

    3.1 數(shù)據(jù)集與數(shù)據(jù)分析

    數(shù)據(jù)集選取西班牙電力市場2015 年1 月1 日—2018年12月31日共4 a的歷史電價數(shù)據(jù)[24],每日采樣24 次數(shù)據(jù),時間間隔為1 h,將訓練集、測試集和驗證集按照8∶1∶1進行劃分。

    電價的波動是多因素共同影響形成的結(jié)果,例如歷史電價、負荷、能源價格、假期、氣溫、氣壓、降水濕度等因素,本文通過皮爾遜相關(guān)系數(shù)計算并給出相關(guān)性大于0.25 的影響因素。電價相關(guān)性分析如表1 所示。利用皮爾遜相關(guān)系數(shù)計算出各外部影響因素對電價的影響程度。皮爾遜相關(guān)系數(shù)是一種衡量向量相似度的指標,其取值范圍為[-1,1],其中0 代表無相關(guān)性,正值代表正相關(guān),負值代表負相關(guān)。通常情況下,(0.5,1.0]代表強相關(guān)性,(0,0.5)代表弱相關(guān)性。皮爾遜相關(guān)系數(shù)的數(shù)學表達式為

    表1 電價相關(guān)性分析Tab.1 Correlation analysis of electricity price

    式中:R為皮爾遜相關(guān)系數(shù);xi和yi為兩個變量;和為兩個變量的均值。

    由于同一地區(qū)電價變化具有周期性,故本文給出了西班牙2016年5月—2017年5月的電價曲線,如圖6 所示。由圖6 可知,西班牙電力市場的電價序列具有非線性特點,且電價大約在第240 天(即2017年1月)達到最高。

    圖6 西班牙市場年電價曲線Fig.6 Curve of annual electricity price on Spanish market

    3.2 數(shù)據(jù)預處理

    通常用直接刪除、平均值補全、K 最近鄰法KNN(K-nearest neighbor)[25]對原始數(shù)據(jù)集中的個別異常值、缺失值進行處理,以更好學習歷史電價序列中的潛在規(guī)律和提高預測的準確性。為了將數(shù)據(jù)映射在[-1,1],消除異常樣本對預測結(jié)果產(chǎn)生的不良影響,降低訓練結(jié)果的誤差,需要對數(shù)據(jù)進行歸一化。歸一化常用方法有min-max 歸一化和Zscore 歸一化。由于min-max 歸一化在性能上優(yōu)于Z-score 歸一化,故本文采用min-max 歸一化,其數(shù)學表達式為

    式中:x為歸一化后的數(shù)據(jù);x1為原始數(shù)據(jù);xmin和xmax分別為樣本數(shù)據(jù)的最小值和最大值。

    3.3 評價指標

    本文選用均方誤差根RMSE(root mean squared error)、平均絕對百分比誤差MAPE(mean absolute percentage error)和平均絕對誤差MAE(mean absolute error)作為算法的評價指標,模型的RMSE、MAPE、MAE 越小代表預測結(jié)果越準確。3 種評價指標的數(shù)學表達式分別為

    式中:N為樣本數(shù)量;yi為歷史電價實際值;為電價預測值。

    3.4 滑動窗口不同寬度對比分析

    本文實驗采用滑動窗口處理,滑動窗口寬度對預測結(jié)果的精確度存在一定影響?;瑒哟翱趯挾冗^小可能會導致模型無法捕捉到長期依賴關(guān)系,從而影響預測精度;而滑動窗口寬度過大可能會使模型過于依賴歷史信息,導致對未來的預測結(jié)果不準確。本文將時間步設(shè)置為4,滑動窗口寬度分別設(shè)置為24的倍數(shù)進行對比分析。由對滑動窗口進行尋優(yōu)實驗可知,當時間步為4 時,將窗口寬度設(shè)置為24,本文模型的評價指標最小、預測結(jié)果最優(yōu)。首先利用滑動窗口處理截取前24 h 的歷史電價序列,從而對下一個時刻的電價進行預測;然后滑動到下一個新窗口,不斷重復此過程進行迭代;最后通過不斷迭代得到未來72 h 的電價預測值。需要說明的是,原始電價序列限定在指定的長度。

    表2 為不同滑動窗口寬度情況下本文所提模型在測試集上的評價指標對比。

    3.5 模型是否引入注意力機制對比分析

    為驗證本文所提出模型引入注意力機制后的優(yōu)越性,將去除注意力機制的本文模型記作CNNBIGRU,對是否引入注意力機制進行對比,實驗結(jié)果如圖7和表3所示。

    表3 是否引入注意力機制的評價指標對比Tab.3 Comparison of evaluation indicators with and without the introduction of Attention mechanism

    圖7 是否引入注意力機制對比Fig.7 Comparison of electricity price with and without the introduction of Attention mechanism

    由圖7 可知,本文所提模型的預測電價曲線與真實電價曲線更接近,而CNN-BIGRU 模型的預測精度明顯低于本文所提模型。表3 為是否引入注意力機制的評價指標對比。對比分析可知,本文模型比CNN-BIGRU 的RMSE、MAPE、MAE 分別減少了14.155%、7.962%、6.058%。結(jié)果表明,注意力機制對電價的影響因素的權(quán)重進行自適應(yīng)調(diào)整能有效提高模型的預測精度。

    3.6 模型是否采用CNN 對比分析

    為驗證本文模型采用CNN 進行特征提取對預測精度的影響,將去除CNN 的本文所提模型記作BIGRU-Attention,兩者進行對比,實驗結(jié)果如圖8和表4所示。

    表4 是否采用CNN 的評價指標對比Tab.4 Comparison of evaluation indicators with and without the use of CNN

    圖8 是否采用CNN 對比Fig.8 Comparison of electricity price with and without of the use of CNN

    由圖8 可知,本文模型的預測電價曲線與真實電價曲線更接近,而BIGRU-Attention的預測精度低于本文所提模型。表4 為是否采用CNN 的評價指標對比。對比分析可知,本文模型比BIGRU-Attention 的RMSE、MAPE、MAE 分別減少了11.385%、14.793%、16.910%,驗證了CNN 特征提取功能的重要性。

    3.7 模型是否采用雙向GRU 對比分析

    為驗證本文模型采用雙向GRU 的優(yōu)越性,將本文模型中的雙向GRU 改為單向GRU,記作CNNGRU-Attention,兩者進行對比,實驗結(jié)果如圖9 和表5所示。

    表5 是否采用雙向GRU 的評價指標對比Tab.5 Comparison of evaluation indicators with and without the use of BIGRU

    圖9 是否采用雙向GRU 對比Fig.9 Comparison of electricity price with and without the use of BIGRU

    由圖9 可知,本文模型的預測電價曲線與真實電價曲線更接近,而CNN-GRU-Attention 的預測精度略低于本文模型。表5 為是否采用雙向GRU 的評價指標對比。對比分析可知,本文模型比CNNGRU-Attention 的RMSE、MAPE、MAE 分別減少了0.202%、1.370%、1.322%。結(jié)果表明,相比于單向GRU,雙向GRU 更能挖掘和學習數(shù)據(jù)序列中的規(guī)律,但是其學習能力有限。

    3.8 不同模型對比分析

    為進一步驗證本文模型的優(yōu)越性,將本文模型與ARIMA、CNN-LSTM、支持向量回歸SVR(support vector regression)、多層感知機MLP(Multilayer Perceptron)模型進行對比,實驗結(jié)果如圖10 和表6所示。

    表6 不同模型的評價指標對比Tab.6 Comparison of evaluation indicators among different models

    圖10 不同模型對比Fig.10 Comparison among different models

    由圖10 可知,相比于其他模型,本文模型的預測電價曲線與真實電價曲線更接近,對歷史電價有更強的學習能力。表6 為不同模型的評價指標對比,可以看出,本文模型的RMSE、MAPE、MAE均低于其他模型,并且本文模型比ARIMA 模型的RMSE、MAPE、MAE分別減少了31.875%、16.280%、19.379%。上述實驗結(jié)果進一步驗證了本文模型的精確性。

    4 結(jié)論

    為了提高電價預測的精確度,本文提出了一種基于注意力機制的CNN-BIGRU 預測電價模型,將本文模型與多種模型進行對比分析,驗證了本文模型的精確性。主要結(jié)論如下。

    (1)針對龐大的數(shù)據(jù)集,利用CNN模型對特征進行提取,充分挖掘電價序列的特征,提高預測模型的精確度,避免了因人工提取特征導致的失誤。

    (2)采用BIGRU模型進行訓練,使得信息可以雙向傳遞,更好地學習電價序列的變化規(guī)律,以提高預測模型的精確度。

    (3)引入注意力機制對影響電價變化的因素的權(quán)重進行合理設(shè)置,從而進一步提高了預測模型的精確度。

    后續(xù)研究可以從考慮電力市場中的政策變化、用戶行為等因素來進一步提高預測的精確度。

    猜你喜歡
    特征提取機制模型
    一半模型
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    基于Gazebo仿真環(huán)境的ORB特征提取與比對的研究
    電子制作(2019年15期)2019-08-27 01:12:00
    自制力是一種很好的篩選機制
    文苑(2018年21期)2018-11-09 01:23:06
    一種基于LBP 特征提取和稀疏表示的肝病識別算法
    3D打印中的模型分割與打包
    破除舊機制要分步推進
    基于MED和循環(huán)域解調(diào)的多故障特征提取
    注重機制的相互配合
    亚洲欧美日韩卡通动漫| 亚洲av熟女| 日本与韩国留学比较| 人妻制服诱惑在线中文字幕| 国产成人a区在线观看| 成年女人毛片免费观看观看9| 久久午夜福利片| 男女床上黄色一级片免费看| 国产精品久久久久久亚洲av鲁大| 一本综合久久免费| 欧美高清性xxxxhd video| 一级作爱视频免费观看| 丰满乱子伦码专区| eeuss影院久久| 1000部很黄的大片| 在线播放无遮挡| 欧美黑人欧美精品刺激| 久久精品人妻少妇| 看十八女毛片水多多多| 看免费av毛片| 又爽又黄无遮挡网站| 午夜老司机福利剧场| 欧美一区二区国产精品久久精品| 91av网一区二区| 永久网站在线| 黄色配什么色好看| 中文字幕久久专区| 如何舔出高潮| 久99久视频精品免费| 天天躁日日操中文字幕| 变态另类丝袜制服| 精品久久久久久久久亚洲 | bbb黄色大片| 99国产精品一区二区蜜桃av| 国产精品一区二区三区四区久久| 99国产极品粉嫩在线观看| 在线a可以看的网站| 国产色婷婷99| 亚洲国产精品久久男人天堂| 99在线视频只有这里精品首页| 欧美一区二区亚洲| 亚洲av免费在线观看| 18禁黄网站禁片免费观看直播| 国内精品美女久久久久久| 免费看光身美女| 国产综合懂色| 99热这里只有是精品在线观看 | 国产伦一二天堂av在线观看| 午夜久久久久精精品| 桃色一区二区三区在线观看| 十八禁人妻一区二区| 国产一区二区三区在线臀色熟女| 国产精品免费一区二区三区在线| 国产午夜精品久久久久久一区二区三区 | 18禁裸乳无遮挡免费网站照片| 亚洲av电影在线进入| 亚洲自偷自拍三级| 美女免费视频网站| 好男人电影高清在线观看| 国产午夜福利久久久久久| 99国产精品一区二区三区| 看免费av毛片| 高清在线国产一区| 国产精品永久免费网站| 亚洲av免费高清在线观看| 国产黄色小视频在线观看| 国产精品三级大全| 国产主播在线观看一区二区| 99国产极品粉嫩在线观看| 国内精品美女久久久久久| 国产免费一级a男人的天堂| .国产精品久久| 亚洲无线观看免费| 亚洲三级黄色毛片| 日本在线视频免费播放| 成年免费大片在线观看| 亚洲成av人片免费观看| 色播亚洲综合网| 露出奶头的视频| 亚洲自偷自拍三级| 久久精品国产亚洲av香蕉五月| 欧美又色又爽又黄视频| 国产极品精品免费视频能看的| 日韩精品中文字幕看吧| 在线播放无遮挡| 人妻丰满熟妇av一区二区三区| 人妻久久中文字幕网| 国产成人aa在线观看| 亚洲五月婷婷丁香| 丝袜美腿在线中文| 亚洲最大成人av| 久久人人爽人人爽人人片va | 成人特级av手机在线观看| 宅男免费午夜| av欧美777| 久久这里只有精品中国| 丰满的人妻完整版| 久久国产乱子免费精品| 99国产极品粉嫩在线观看| 国产乱人视频| 久久久久亚洲av毛片大全| 日韩欧美精品v在线| 国产精品久久久久久人妻精品电影| 国产精品一区二区三区四区久久| 丰满的人妻完整版| 国产私拍福利视频在线观看| 美女大奶头视频| 中文字幕人成人乱码亚洲影| 精品久久久久久久久亚洲 | 熟女电影av网| 国产毛片a区久久久久| 在线播放无遮挡| 国产av麻豆久久久久久久| 成人欧美大片| 99久久精品一区二区三区| 日本一本二区三区精品| 欧美最黄视频在线播放免费| 又黄又爽又刺激的免费视频.| 99热这里只有精品一区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品456在线播放app | 天天躁日日操中文字幕| 90打野战视频偷拍视频| 久久久久久国产a免费观看| 一进一出好大好爽视频| av欧美777| 黄色女人牲交| 亚洲国产精品999在线| 精品人妻偷拍中文字幕| 国产精品一区二区性色av| 欧美不卡视频在线免费观看| 欧美精品啪啪一区二区三区| 99国产精品一区二区三区| 欧美+日韩+精品| 欧美高清成人免费视频www| 在线播放无遮挡| 久久久久久久久久成人| 亚洲一区二区三区不卡视频| 丁香六月欧美| 99精品在免费线老司机午夜| 国产色婷婷99| 亚洲在线观看片| 看免费av毛片| 国产精品一区二区免费欧美| 国产综合懂色| 国产色婷婷99| 国产激情偷乱视频一区二区| 啦啦啦韩国在线观看视频| 国产又黄又爽又无遮挡在线| 国产在线精品亚洲第一网站| 成人美女网站在线观看视频| 国产精品电影一区二区三区| 国产探花极品一区二区| 国内精品美女久久久久久| 国内精品美女久久久久久| 日韩人妻高清精品专区| 自拍偷自拍亚洲精品老妇| 丁香欧美五月| 日韩成人在线观看一区二区三区| eeuss影院久久| 国产精品一及| av在线天堂中文字幕| 亚洲国产精品成人综合色| 90打野战视频偷拍视频| 色哟哟·www| 露出奶头的视频| 国产精品自产拍在线观看55亚洲| 国产精品爽爽va在线观看网站| 欧美性感艳星| 狠狠狠狠99中文字幕| 久久亚洲真实| 嫩草影院入口| 久久天躁狠狠躁夜夜2o2o| 久久午夜亚洲精品久久| 精品熟女少妇八av免费久了| 日本撒尿小便嘘嘘汇集6| 亚洲第一欧美日韩一区二区三区| 高清在线国产一区| 窝窝影院91人妻| 精品福利观看| 在线观看av片永久免费下载| 亚洲成人久久性| 在线观看一区二区三区| 亚洲激情在线av| 久久久久久国产a免费观看| 狠狠狠狠99中文字幕| 色精品久久人妻99蜜桃| 午夜精品久久久久久毛片777| 99精品久久久久人妻精品| 欧美色欧美亚洲另类二区| 特大巨黑吊av在线直播| 久久久久精品国产欧美久久久| 成人一区二区视频在线观看| 免费一级毛片在线播放高清视频| 9191精品国产免费久久| 亚洲国产精品合色在线| 999久久久精品免费观看国产| 亚洲美女视频黄频| 十八禁网站免费在线| 亚洲七黄色美女视频| 国产精品久久视频播放| 国产精品美女特级片免费视频播放器| 亚洲成a人片在线一区二区| 亚洲人成电影免费在线| 99精品久久久久人妻精品| 一级黄片播放器| 欧美3d第一页| 成人欧美大片| 九九在线视频观看精品| av在线观看视频网站免费| 亚洲av免费在线观看| 制服丝袜大香蕉在线| 91在线精品国自产拍蜜月| 亚洲精华国产精华精| 少妇人妻精品综合一区二区 | www.999成人在线观看| 亚洲精品日韩av片在线观看| a级毛片a级免费在线| 亚洲精品久久国产高清桃花| 国产中年淑女户外野战色| 老司机午夜福利在线观看视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久久久久久久中文| 婷婷亚洲欧美| 一区二区三区激情视频| 亚洲真实伦在线观看| 色吧在线观看| 国产三级黄色录像| av视频在线观看入口| 成年人黄色毛片网站| 国产日本99.免费观看| 精品久久久久久久久久久久久| 欧美潮喷喷水| 99久久九九国产精品国产免费| 亚洲一区二区三区色噜噜| 亚洲美女黄片视频| 久久久久久久久大av| 十八禁网站免费在线| 日韩欧美免费精品| 国产精华一区二区三区| 很黄的视频免费| 国产精品久久电影中文字幕| 欧美潮喷喷水| 波野结衣二区三区在线| 亚洲精品亚洲一区二区| 久久精品夜夜夜夜夜久久蜜豆| 色哟哟·www| 欧美黑人巨大hd| 欧美xxxx性猛交bbbb| 国产成人啪精品午夜网站| 亚洲av中文字字幕乱码综合| 国内精品一区二区在线观看| 国产一级毛片七仙女欲春2| 国产真实伦视频高清在线观看 | 自拍偷自拍亚洲精品老妇| 高清日韩中文字幕在线| 乱人视频在线观看| a在线观看视频网站| 日本撒尿小便嘘嘘汇集6| 国产 一区 欧美 日韩| 99久久99久久久精品蜜桃| 麻豆成人av在线观看| 大型黄色视频在线免费观看| 国产高清激情床上av| 成年女人永久免费观看视频| 亚洲精品影视一区二区三区av| 国产视频内射| 国产激情偷乱视频一区二区| 午夜福利在线观看免费完整高清在 | 色av中文字幕| 黄色丝袜av网址大全| 99久久99久久久精品蜜桃| av福利片在线观看| 日本五十路高清| 宅男免费午夜| 精品免费久久久久久久清纯| 成人性生交大片免费视频hd| 久久久久精品国产欧美久久久| 欧美日本视频| 动漫黄色视频在线观看| av天堂在线播放| 欧美黄色片欧美黄色片| 91在线精品国自产拍蜜月| 99在线视频只有这里精品首页| 亚洲久久久久久中文字幕| 亚洲成av人片在线播放无| 亚洲中文字幕日韩| 中文字幕人成人乱码亚洲影| 亚洲av电影在线进入| 亚洲欧美清纯卡通| 亚洲av不卡在线观看| 国产激情偷乱视频一区二区| 国产精品久久久久久久电影| 国产一区二区在线观看日韩| 久久精品综合一区二区三区| 日韩欧美精品免费久久 | 日韩欧美 国产精品| 日韩欧美一区二区三区在线观看| 九九热线精品视视频播放| 免费观看的影片在线观看| 在线看三级毛片| 亚洲在线自拍视频| 国内精品一区二区在线观看| 此物有八面人人有两片| 99热只有精品国产| 亚洲黑人精品在线| 99热这里只有是精品在线观看 | 草草在线视频免费看| 午夜亚洲福利在线播放| 欧美最新免费一区二区三区 | 国产欧美日韩一区二区三| 深爱激情五月婷婷| 热99re8久久精品国产| 人妻丰满熟妇av一区二区三区| 三级男女做爰猛烈吃奶摸视频| 亚洲国产精品999在线| 国产精品日韩av在线免费观看| 性色avwww在线观看| 色av中文字幕| 搞女人的毛片| 神马国产精品三级电影在线观看| 少妇的逼好多水| 男人和女人高潮做爰伦理| 国产在视频线在精品| 国产乱人伦免费视频| 91av网一区二区| 免费看美女性在线毛片视频| bbb黄色大片| 免费一级毛片在线播放高清视频| 国产精品久久久久久久电影| 我的女老师完整版在线观看| 国产精品久久久久久人妻精品电影| 亚洲一区高清亚洲精品| 成人精品一区二区免费| 国产一级毛片七仙女欲春2| 国产91精品成人一区二区三区| 久久6这里有精品| 性欧美人与动物交配| 欧美高清成人免费视频www| 国内少妇人妻偷人精品xxx网站| www.999成人在线观看| 午夜免费激情av| 国产aⅴ精品一区二区三区波| 一级黄色大片毛片| avwww免费| 免费人成视频x8x8入口观看| 蜜桃亚洲精品一区二区三区| 老司机深夜福利视频在线观看| av在线观看视频网站免费| 一个人免费在线观看的高清视频| 精品人妻熟女av久视频| 免费在线观看日本一区| 午夜福利在线观看免费完整高清在 | 国产激情偷乱视频一区二区| 国产午夜福利久久久久久| 宅男免费午夜| 性色avwww在线观看| 又爽又黄a免费视频| 女同久久另类99精品国产91| 国产私拍福利视频在线观看| 桃红色精品国产亚洲av| 九色国产91popny在线| 欧美区成人在线视频| 国产精品免费一区二区三区在线| 国模一区二区三区四区视频| 久久精品国产亚洲av天美| 婷婷亚洲欧美| 搡女人真爽免费视频火全软件 | 中出人妻视频一区二区| 熟妇人妻久久中文字幕3abv| 精品久久久久久,| 国产精品精品国产色婷婷| 女人被狂操c到高潮| 99国产精品一区二区三区| 久久精品久久久久久噜噜老黄 | 免费无遮挡裸体视频| 99热只有精品国产| 欧美成人a在线观看| 别揉我奶头 嗯啊视频| 国产黄片美女视频| 久久这里只有精品中国| 中国美女看黄片| 久久人妻av系列| 日本黄色视频三级网站网址| 亚州av有码| 好男人在线观看高清免费视频| 日本一本二区三区精品| 日韩人妻高清精品专区| or卡值多少钱| 免费看日本二区| 老熟妇乱子伦视频在线观看| 男女下面进入的视频免费午夜| 日本 欧美在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品久久久久久久电影| 我的女老师完整版在线观看| 噜噜噜噜噜久久久久久91| 国产亚洲精品久久久com| 99久久精品国产亚洲精品| 淫妇啪啪啪对白视频| 夜夜看夜夜爽夜夜摸| 久久久久久久久久黄片| 给我免费播放毛片高清在线观看| 国产色爽女视频免费观看| 精品国产亚洲在线| 麻豆国产97在线/欧美| 久久久国产成人免费| 日日摸夜夜添夜夜添小说| 最好的美女福利视频网| 亚洲七黄色美女视频| 久久伊人香网站| 亚洲国产高清在线一区二区三| 麻豆国产97在线/欧美| 久久99热6这里只有精品| 久久久久国产精品人妻aⅴ院| 他把我摸到了高潮在线观看| 精品欧美国产一区二区三| 天堂av国产一区二区熟女人妻| 国产三级在线视频| 最近中文字幕高清免费大全6 | 久久久精品欧美日韩精品| 老司机午夜福利在线观看视频| 免费av不卡在线播放| 长腿黑丝高跟| aaaaa片日本免费| 欧美高清性xxxxhd video| 亚洲av五月六月丁香网| 美女黄网站色视频| 简卡轻食公司| 欧美绝顶高潮抽搐喷水| 丰满的人妻完整版| 国产淫片久久久久久久久 | 亚洲国产高清在线一区二区三| 欧美一区二区国产精品久久精品| 国产单亲对白刺激| .国产精品久久| 欧美精品国产亚洲| 国产视频内射| 又紧又爽又黄一区二区| 在线播放国产精品三级| 国产亚洲精品av在线| 亚洲片人在线观看| 18禁裸乳无遮挡免费网站照片| 国产精品电影一区二区三区| av在线观看视频网站免费| 91久久精品国产一区二区成人| 国产 一区 欧美 日韩| 久久热精品热| 日本a在线网址| 中文在线观看免费www的网站| 99久久久亚洲精品蜜臀av| 最近视频中文字幕2019在线8| 亚洲欧美精品综合久久99| 熟女电影av网| 九九热线精品视视频播放| 国产单亲对白刺激| 久久香蕉精品热| 天堂动漫精品| 性欧美人与动物交配| av福利片在线观看| 一二三四社区在线视频社区8| 在线免费观看不下载黄p国产 | 日日夜夜操网爽| 2021天堂中文幕一二区在线观| 欧美黑人巨大hd| 美女被艹到高潮喷水动态| 亚洲精品成人久久久久久| 久久九九热精品免费| 18+在线观看网站| 搡女人真爽免费视频火全软件 | 国产在视频线在精品| 神马国产精品三级电影在线观看| 在线观看午夜福利视频| 亚洲在线自拍视频| 深夜a级毛片| 亚洲欧美清纯卡通| 9191精品国产免费久久| www.色视频.com| 可以在线观看的亚洲视频| 人人妻人人看人人澡| 国产伦精品一区二区三区视频9| 能在线免费观看的黄片| 两性午夜刺激爽爽歪歪视频在线观看| 黄片小视频在线播放| 国产亚洲欧美98| 日韩中文字幕欧美一区二区| 成人特级av手机在线观看| 久久精品国产99精品国产亚洲性色| 亚洲中文字幕日韩| av在线天堂中文字幕| 精品久久久久久久久久久久久| 久久久久亚洲av毛片大全| 久久久久九九精品影院| 国产白丝娇喘喷水9色精品| 国产精品精品国产色婷婷| 在线播放国产精品三级| 99国产极品粉嫩在线观看| 久久天躁狠狠躁夜夜2o2o| 中文资源天堂在线| 美女黄网站色视频| 亚洲一区高清亚洲精品| 亚洲中文字幕日韩| 欧美日韩乱码在线| 99热这里只有是精品50| 少妇人妻一区二区三区视频| 国产精华一区二区三区| 亚洲最大成人中文| 国产亚洲欧美在线一区二区| 国产成人啪精品午夜网站| 波多野结衣高清无吗| 毛片女人毛片| 看十八女毛片水多多多| 国产精品自产拍在线观看55亚洲| 中文在线观看免费www的网站| 亚洲美女黄片视频| 成人av在线播放网站| 少妇高潮的动态图| 日韩欧美免费精品| 99久久精品一区二区三区| 国产私拍福利视频在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲一区二区三区色噜噜| av视频在线观看入口| 小蜜桃在线观看免费完整版高清| 在线观看一区二区三区| 国产成人aa在线观看| 国产老妇女一区| 丝袜美腿在线中文| 一个人看视频在线观看www免费| 国产精品三级大全| 国产色爽女视频免费观看| 午夜福利欧美成人| 日本黄色片子视频| 欧美日韩黄片免| 我要搜黄色片| 高潮久久久久久久久久久不卡| 草草在线视频免费看| 成年女人毛片免费观看观看9| 久久精品人妻少妇| 亚洲五月婷婷丁香| 国产白丝娇喘喷水9色精品| 超碰av人人做人人爽久久| a级一级毛片免费在线观看| 久久精品国产清高在天天线| 欧美区成人在线视频| 成年女人永久免费观看视频| 久久精品国产亚洲av天美| 日韩成人在线观看一区二区三区| 成人性生交大片免费视频hd| 在线a可以看的网站| 日韩欧美免费精品| 99久久精品一区二区三区| 深夜精品福利| 久久国产精品人妻蜜桃| 久久亚洲精品不卡| АⅤ资源中文在线天堂| 欧美国产日韩亚洲一区| 中文字幕人成人乱码亚洲影| 99热6这里只有精品| 日本免费一区二区三区高清不卡| 午夜精品在线福利| 亚洲精品一卡2卡三卡4卡5卡| 女人十人毛片免费观看3o分钟| 国产精品亚洲美女久久久| 亚洲av电影在线进入| 亚洲av中文字字幕乱码综合| 天堂√8在线中文| 脱女人内裤的视频| 十八禁网站免费在线| 久久久久久久亚洲中文字幕 | 国产成人aa在线观看| 亚洲一区二区三区色噜噜| 亚洲av熟女| 少妇熟女aⅴ在线视频| 黄色视频,在线免费观看| 99久久精品一区二区三区| 国产av不卡久久| 精品一区二区三区av网在线观看| 亚洲激情在线av| 首页视频小说图片口味搜索| 亚洲美女黄片视频| 最新在线观看一区二区三区| 亚洲自拍偷在线| 国产精华一区二区三区| 在线天堂最新版资源| 免费看美女性在线毛片视频| 久久久久免费精品人妻一区二区| 午夜福利免费观看在线| 小蜜桃在线观看免费完整版高清| 久久久久久久午夜电影| 9191精品国产免费久久| 我要看日韩黄色一级片| 亚洲一区二区三区不卡视频| 制服丝袜大香蕉在线| 日韩国内少妇激情av| 久久久久久久精品吃奶| 看十八女毛片水多多多| 亚洲成av人片免费观看| 久久久久亚洲av毛片大全| 国产精品久久久久久久久免 | 一边摸一边抽搐一进一小说| 波野结衣二区三区在线| 久久久久国内视频| 精品国产亚洲在线| 亚洲欧美日韩高清专用| 国产高清视频在线播放一区| 亚洲成人久久爱视频| 成人国产一区最新在线观看| 欧美色欧美亚洲另类二区| 村上凉子中文字幕在线| 老司机午夜福利在线观看视频| 午夜福利欧美成人| 别揉我奶头 嗯啊视频| 宅男免费午夜| 啦啦啦观看免费观看视频高清|