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

    基于長短期記憶神經(jīng)網(wǎng)絡(luò)的風(fēng)力發(fā)電功率預(yù)測方法

    2019-11-06 12:33:18李相俊許格健
    發(fā)電技術(shù) 2019年5期
    關(guān)鍵詞:方法模型

    李相俊,許格健

    基于長短期記憶神經(jīng)網(wǎng)絡(luò)的風(fēng)力發(fā)電功率預(yù)測方法

    李相俊,許格健

    (新能源與儲能運行控制國家重點實驗室(中國電力科學(xué)研究院有限公司),北京市 海淀區(qū) 100192)

    風(fēng)力發(fā)電過程具有較強(qiáng)的隨機(jī)性,導(dǎo)致風(fēng)力發(fā)電功率的預(yù)測準(zhǔn)確度不高。針對上述問題,提出了一種融合深度學(xué)習(xí)算法的風(fēng)力發(fā)電功率預(yù)測方法。以歷史風(fēng)力發(fā)電功率數(shù)據(jù)作為輸入,建立風(fēng)力發(fā)電功率預(yù)測模型,實現(xiàn)對未來一個時間刻度的風(fēng)力發(fā)電功率預(yù)測。算例結(jié)果表明,與傳統(tǒng)時序預(yù)測方法相比,基于長短期記憶神經(jīng)網(wǎng)絡(luò)的風(fēng)力發(fā)電功率預(yù)測結(jié)果在各項指標(biāo)中誤差更小,驗證了上述方法在風(fēng)力發(fā)電功率預(yù)測中的可行性和有效性,提升了風(fēng)力發(fā)電功率預(yù)測的準(zhǔn)確性。

    深度學(xué)習(xí);時序預(yù)測;風(fēng)力發(fā)電;長短期記憶神經(jīng)網(wǎng)絡(luò)

    0 引言

    隨著能源與環(huán)境問題的日益突出[1],風(fēng)能作為一種無污染、可持續(xù)的清潔能源成為了能源發(fā)展的關(guān)注點。國內(nèi)風(fēng)力發(fā)電水平也保持著穩(wěn)步增長,2018年風(fēng)力發(fā)電量達(dá)到了3 660億kW·h,占全部發(fā)電量的5.2%[2]。在風(fēng)力發(fā)電規(guī)模穩(wěn)步增長的同時,其本身存在的問題也被逐漸放大。由于受多種環(huán)境因素的影響,風(fēng)力發(fā)電的波動性及隨機(jī)性[3-5]成為了制約其發(fā)展的重要因素。為了解決這一問題,需進(jìn)行風(fēng)力發(fā)電功率預(yù)測研究,以判斷其變化趨勢并提前采取相應(yīng)措施。

    物理模型預(yù)測方法與數(shù)據(jù)模型預(yù)測方法是目前針對風(fēng)力發(fā)電功率預(yù)測的兩大主要方法。物理模型預(yù)測方法主要考慮數(shù)值天氣預(yù)報(numerical weather prediction,NWP) 及其他多種因素的影響?;谠摲椒ǖ娘L(fēng)力發(fā)電功率預(yù)測需要較多的參數(shù),如歷史數(shù)據(jù)、地形地貌、氣象數(shù)據(jù)等[6-13],導(dǎo)致其在短期預(yù)測能力上表現(xiàn)較為一般,其結(jié)果更適合作為長期預(yù)測的參考標(biāo)準(zhǔn)。

    數(shù)學(xué)模型預(yù)測方法是目前研究最廣泛的風(fēng)力發(fā)電功率預(yù)測方法,在數(shù)學(xué)預(yù)測模型中,針對風(fēng)力發(fā)電功率的預(yù)測分為2種形式。一種預(yù)測方法是考慮多種因素對風(fēng)力發(fā)電功率的影響,通過建立數(shù)學(xué)模型得到相應(yīng)的非線性回歸曲線,并依此進(jìn)行風(fēng)力發(fā)電功率預(yù)測。傳統(tǒng)的數(shù)學(xué)模型預(yù)測方法,如小波算法(wavelet)[14-15]、支持向量機(jī)(support vector machine,SVM)[16]等,都是基于風(fēng)力發(fā)電功率的相關(guān)特征進(jìn)行短時預(yù)測,通過對影響風(fēng)速以及風(fēng)電力發(fā)電功率的特征進(jìn)行分析,預(yù)測未來一段時間內(nèi)的風(fēng)力發(fā)電功率。另一種預(yù)測方法是不考慮其余環(huán)境條件的影響,只考慮風(fēng)力發(fā)電功率本身隨時間變化的特性趨勢,并依此進(jìn)行風(fēng)力發(fā)電功率時序預(yù)測分析。自回歸移動平均(autoregressive moving average,ARMA)模型[17]是時序預(yù)測的主要應(yīng)用算法,通過量化當(dāng)前數(shù) 據(jù)與之前數(shù)據(jù)的關(guān)系,同時解決隨機(jī)變化項的 求解問題,進(jìn)一步求得需求預(yù)測的下一時刻 數(shù)據(jù)。

    除傳統(tǒng)預(yù)測方法外,人工神經(jīng)網(wǎng)絡(luò)也常應(yīng)用于風(fēng)力發(fā)電功率的短時預(yù)測中。由于神經(jīng)網(wǎng)絡(luò)形式多樣,其對應(yīng)的預(yù)測方法也不同。如李彬等[18]提出了基于人工神經(jīng)網(wǎng)絡(luò)的預(yù)測方法,以多種特征因素作為神經(jīng)網(wǎng)絡(luò)輸入,預(yù)測得到其當(dāng)月的風(fēng)力發(fā)電功率概率分布曲線。江岳春等[19]將傳統(tǒng)BP神經(jīng)網(wǎng)絡(luò)與多種算法相結(jié)合,對風(fēng)力發(fā)電功率的時間序列進(jìn)行短時時序預(yù)測,取得了良好的預(yù)測效果。由于神經(jīng)網(wǎng)絡(luò)出色的非線性數(shù)據(jù)擬合能力,在許多針對風(fēng)力發(fā)電功率的短時及超短時預(yù)測方法中,都會應(yīng)用神經(jīng)網(wǎng)絡(luò)算法對歷史數(shù)據(jù)進(jìn)行分析,再采取其余預(yù)測方法進(jìn)行風(fēng)力發(fā)電功率及其分布概率的預(yù)測[20-24]。

    隨著人工智能及深度學(xué)習(xí)的不斷推廣應(yīng) 用[25-26],各種針對不同數(shù)據(jù)的深度學(xué)習(xí)算法逐漸成為各領(lǐng)域?qū)<谊P(guān)注的焦點,也在新能源領(lǐng)域得到了廣泛的應(yīng)用[27-28]。其中循環(huán)神經(jīng)網(wǎng)絡(luò)(recurrent neural networks,RNN)作為深度學(xué)習(xí)中一種典型的序列分析模型,在語言翻譯、流量預(yù)測等方面都得到了廣泛的應(yīng)用。但在應(yīng)用過程中,RNN算法的序列數(shù)據(jù)處理能力在面對大規(guī)模數(shù)據(jù)進(jìn)行跨度較大的序列預(yù)測時,易出現(xiàn)梯度消失和梯度爆炸的問題。在RNN的發(fā)展過程中,一種特殊的循環(huán)神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu),即長短期記憶(long short-term memory,LSTM)神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu),很大程度上解決了這一問題,更加適用于數(shù)據(jù)規(guī)模較大的序列數(shù)據(jù)[26],并且在預(yù)測精度上有著良好的表現(xiàn)。但是該方法在風(fēng)力發(fā)電功率預(yù)測方面的應(yīng)用國內(nèi)外鮮見報道。

    本文提出一種基于 LSTM神經(jīng)網(wǎng)絡(luò)的風(fēng)力發(fā)電功率時序預(yù)測方法,采用某風(fēng)電場的發(fā)電數(shù)據(jù)進(jìn)行時序預(yù)測分析,并與傳統(tǒng)時序預(yù)測方法進(jìn)行對比分析。結(jié)果表明,該方法可有效提高風(fēng)力發(fā)電功率的預(yù)測精度,并且針對不同規(guī)模數(shù)據(jù)均有良好的預(yù)測效果。

    1 時序預(yù)測方法

    1.1 差分整合移動平均自回歸模型

    在傳統(tǒng)時序算法中,差分整合移動平均自回歸(autoregressive integrated moving average,ARIMA)算法是一種常用的時序預(yù)測算法[29-31],該算法通過將ARMA模型與差分模型相結(jié)合,提升了原有算法的時序預(yù)測能力。ARIMA算法主要分為3部分:1)通過自回歸(autoregressive,AR)模型描述當(dāng)前值與歷史值之間的關(guān)系;2)用變量自身的歷史時間數(shù)據(jù)對自身進(jìn)行預(yù)測;3)通過移動平均法(moving average,MA)消除預(yù)測中的隨機(jī)波動,并通過適當(dāng)?shù)牟罘执螖?shù)使數(shù)據(jù)序列更加平穩(wěn)。

    在ARIMA模型中,AR模型可表示為:

    式中:y表示當(dāng)前值;與r分別為常數(shù)項與常系數(shù);e為誤差;為階數(shù),表示以之前多少歷史值為參考。在AR模型中,主要關(guān)注的是當(dāng)前值與歷史值之間的關(guān)系。

    在ARIMA模型中,MA模型可表示為:

    式中:l為常系數(shù)。在MA模型中,主要關(guān)注的是計算當(dāng)前值的過程中誤差的累加。

    將AR模型與MA模型相結(jié)合所得即為ARMA模型,通過2種模型的結(jié)合,可觀測到在一組時間序列中穩(wěn)定因素的趨勢與作為波動因素的誤差趨勢,將兩者相結(jié)合即可預(yù)測未來時刻會出現(xiàn)的數(shù)值,即

    在ARIMA模型中,為了保證序列數(shù)據(jù)的平穩(wěn),添加了差分模型對原有數(shù)據(jù)進(jìn)行差分運算,整體運算過程仍滿足ARMA模型的運算方法。由于其本身以歷史數(shù)據(jù)預(yù)測自身的特性,因此常用于能源領(lǐng)域的時序功率預(yù)測。ARIMA在計算過程中要經(jīng)過復(fù)雜的數(shù)據(jù)處理過程,而且由于本身對輸入數(shù)據(jù)的要求,其對個別輸入數(shù)據(jù)無法識別,導(dǎo)致其應(yīng)用受到了一定程度的限制。

    1.2 LSTM神經(jīng)網(wǎng)絡(luò)模型

    1.2.1 LSTM神經(jīng)網(wǎng)絡(luò)單元

    深度學(xué)習(xí)中的LSTM神經(jīng)網(wǎng)絡(luò)本質(zhì)上是循環(huán)神經(jīng)網(wǎng)絡(luò)的一種,與普通循環(huán)神經(jīng)網(wǎng)絡(luò)所不同的是,LSTM神經(jīng)網(wǎng)絡(luò)是一種時間尺度較長的短時記憶網(wǎng)絡(luò)。其主要運算方法仍是基于循環(huán)神經(jīng)網(wǎng)絡(luò)的“記憶”基礎(chǔ),其網(wǎng)絡(luò)結(jié)構(gòu)雖然與循環(huán)神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)相似,但其內(nèi)部構(gòu)成及運算更加復(fù)雜。

    LSTM神經(jīng)網(wǎng)絡(luò)單元結(jié)構(gòu)如圖1所示,其網(wǎng)絡(luò)單元相較于傳統(tǒng)循環(huán)神經(jīng)網(wǎng)絡(luò)運算方法的優(yōu)勢在于其特殊的“遺忘”過程,該過程確保LSTM神經(jīng)網(wǎng)絡(luò)可以進(jìn)行時間跨度更長的預(yù)測。

    圖1 LSTM神經(jīng)網(wǎng)絡(luò)單元結(jié)構(gòu)

    LSTM神經(jīng)網(wǎng)絡(luò)的遺忘率是由前一時刻的輸出與該時刻網(wǎng)絡(luò)的輸入共同計算得出的,計算方法如下:

    分別為權(quán)重及截距,將-1時刻的神經(jīng)網(wǎng)絡(luò)輸出-1與時刻的輸入數(shù)據(jù)相結(jié)合進(jìn)行線性變換。表示sigmoid激活函數(shù),通過sigmoid激活函數(shù)將線性變換結(jié)果映射為0~1的衰減函數(shù),得到f。

    在計算遺忘率后,需進(jìn)行同一時刻輸入的記憶運算。所需的衰減系數(shù)記為,當(dāng)前的記憶狀態(tài)記為¢,其對應(yīng)的求取方法為:

    式中:W、W為權(quán)重;b、b為截距。

    1.2.2 LSTM神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

    LSTM神經(jīng)網(wǎng)絡(luò)進(jìn)行預(yù)測通常有2種方法:1)與傳統(tǒng)深度神經(jīng)網(wǎng)絡(luò)相似的反向傳播誤差算法,這種方法考慮多因素輸入,在不同時刻將LSTM展開為隨時間推移的深度神經(jīng)網(wǎng)絡(luò);2)只考慮單一時序序列的實時遞歸學(xué)習(xí)算法。在應(yīng)用上,2種方法各有優(yōu)勢。本文的預(yù)測過程是基于時序算法的風(fēng)力發(fā)電功率預(yù)測,由于只考慮單一變量的時序變化,因此主要應(yīng)用實時遞歸學(xué)習(xí)算法。

    在實時遞歸學(xué)習(xí)算法中,以1.2.1中所述的神經(jīng)網(wǎng)絡(luò)單元結(jié)構(gòu)為LSTM基本網(wǎng)絡(luò)單元,生成的神經(jīng)網(wǎng)絡(luò)單元的連接形式與輸入形式如圖2所示。

    圖2 LSTM神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

    圖2中,神經(jīng)網(wǎng)絡(luò)在每個時刻都進(jìn)行數(shù)據(jù)輸入,經(jīng)過“遺忘”與“記憶”的處理后,保留-1時刻數(shù)據(jù)特征C-1,并結(jié)合時刻的輸入數(shù)據(jù)x計算本時刻的數(shù)據(jù)特征,記作C,其對應(yīng)的計算公式為:

    求取時刻的數(shù)據(jù)特征后,進(jìn)行LSTM神經(jīng)網(wǎng)絡(luò)本身記憶階段的算法實現(xiàn)。時刻的記憶輸出記為h,與之前所求取遺忘部分的衰減系數(shù)相似,計算得到輸出門的衰減系數(shù):

    2 LSTM網(wǎng)絡(luò)風(fēng)力發(fā)電功率預(yù)測

    由于LSTM網(wǎng)絡(luò)的結(jié)構(gòu)需求,在進(jìn)行預(yù)測分析之前,需先進(jìn)行輸入數(shù)據(jù)的預(yù)處理以及LSTM神經(jīng)網(wǎng)絡(luò)超參數(shù)設(shè)定。

    2.1 數(shù)據(jù)預(yù)處理

    使用風(fēng)力發(fā)電功率數(shù)據(jù)進(jìn)行時序預(yù)測時,考慮到神經(jīng)網(wǎng)絡(luò)預(yù)測模型中使用的非線性激活函數(shù)的輸入、輸出范圍,為避免神經(jīng)元飽和,需要對風(fēng)力發(fā)電功率數(shù)據(jù)進(jìn)行歸一化處理。本文所采用的歸一化標(biāo)準(zhǔn)為最小、最大值標(biāo)準(zhǔn)化方法,其對應(yīng)的公式如下:

    式中:為當(dāng)前時刻的風(fēng)力發(fā)電功率數(shù)據(jù)值;min、max分別為所有數(shù)據(jù)中風(fēng)力發(fā)電功率的最小值、最大值;¢為求得的風(fēng)力發(fā)電功率歸一化結(jié)果。

    在數(shù)據(jù)輸入階段完成歸一化后,對實際輸出的預(yù)測數(shù)據(jù)需應(yīng)用歸一化還原函數(shù)將數(shù)據(jù)還原至其原有大小,最小、最大值標(biāo)準(zhǔn)化方法對應(yīng)的還原函數(shù)為:

    式中:¢pre為未經(jīng)還原的LSTM神經(jīng)網(wǎng)絡(luò)預(yù)測輸出值;pre為還原后的風(fēng)力發(fā)電功率預(yù)測值。

    數(shù)據(jù)輸入階段除了對數(shù)據(jù)進(jìn)行歸一化處理之外,還要對輸入數(shù)據(jù)進(jìn)行數(shù)據(jù)分類處理。風(fēng)力發(fā)電功率數(shù)據(jù)分為訓(xùn)練集和測試集2個集合。通過訓(xùn)練集數(shù)據(jù)對LSTM神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練,在訓(xùn)練的過程不斷優(yōu)化神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)及內(nèi)部參數(shù),再用測試集數(shù)據(jù)測試訓(xùn)練后的LSTM神經(jīng)網(wǎng)絡(luò)的預(yù)測效果。

    2.2 確定LSTM神經(jīng)網(wǎng)絡(luò)超參數(shù)

    在建立 LSTM 神經(jīng)網(wǎng)絡(luò)時序風(fēng)電預(yù)測模型的過程中,需進(jìn)行LSTM神經(jīng)網(wǎng)絡(luò)超參數(shù)的設(shè)定。其中網(wǎng)絡(luò)步長(epoch)以及隱層個數(shù)為神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)參數(shù),另外,為了防止網(wǎng)絡(luò)出現(xiàn)過擬合,還需設(shè)定神經(jīng)網(wǎng)絡(luò)中神經(jīng)元失活系數(shù)(dropout)與批尺寸(batch_size)。

    在LSTM神經(jīng)網(wǎng)絡(luò)中,網(wǎng)絡(luò)步長表示在神經(jīng)網(wǎng)絡(luò)運算過程中從輸入到輸出再返回輸入的一次過程。隨著步長的增加,神經(jīng)網(wǎng)絡(luò)的運算精度也得到提高,同時運算時間也會增加。在LSTM神經(jīng)網(wǎng)絡(luò)的數(shù)據(jù)輸入階段,輸入維度表示每次運算輸入數(shù)據(jù)的個數(shù)。神經(jīng)元失活系數(shù)表示在進(jìn)行神經(jīng)網(wǎng)絡(luò)運算的過程中,棄用其中一部分神經(jīng)元進(jìn)行計算的概率。通過應(yīng)用該參數(shù),取代了原本在神經(jīng)網(wǎng)絡(luò)中需要的正則化方法,也避免了過擬合的出現(xiàn)。在神經(jīng)網(wǎng)絡(luò)運算的過程中,通過設(shè)定適當(dāng)?shù)呐叽?,可以?jié)省運算時間,避免全數(shù)據(jù)的同時計算。

    2.3 風(fēng)力發(fā)電功率預(yù)測

    將經(jīng)過歸一化處理的風(fēng)力發(fā)電功率數(shù)據(jù)作為LSTM神經(jīng)網(wǎng)絡(luò)的輸入,以LSTM單元結(jié)構(gòu)以及網(wǎng)絡(luò)連接結(jié)構(gòu)為基礎(chǔ),并完成LSTM神經(jīng)網(wǎng)絡(luò)超參數(shù)的設(shè)定,得到LSTM神經(jīng)網(wǎng)絡(luò)風(fēng)力發(fā)電功率預(yù)測流程,如圖3所示。其運算公式如下:

    式中:pre表示LSTM神經(jīng)網(wǎng)絡(luò)預(yù)測模型;epoch,batch_size,dropout根據(jù)網(wǎng)絡(luò)超參數(shù)確定。

    圖3 LSTM神經(jīng)網(wǎng)絡(luò)風(fēng)力發(fā)電功率預(yù)測流程

    ¢由式(10)計算得出,為風(fēng)力發(fā)電功率經(jīng)數(shù)據(jù)預(yù)處理后所得結(jié)果,以¢作為輸入求取對應(yīng)的預(yù)測值¢pre。在進(jìn)行預(yù)測的過程中,LSTM神經(jīng)網(wǎng)絡(luò)會根據(jù)設(shè)定的網(wǎng)絡(luò)步長進(jìn)行循環(huán)計算,不斷更新LSTM神經(jīng)網(wǎng)絡(luò)內(nèi)部參數(shù)在循環(huán)結(jié)束后,根據(jù)式(11)將輸出結(jié)果還原為原有風(fēng)力發(fā)電功率大小的pre,即為所求取的LSTM神經(jīng)網(wǎng)絡(luò)風(fēng)力發(fā)電功率預(yù)測數(shù)據(jù)。

    2.4 誤差計算

    本文選用2種不同的誤差標(biāo)準(zhǔn)進(jìn)行預(yù)測效果比較,其中一種誤差標(biāo)準(zhǔn)為均方根誤差(root mean squared error,RMSE),該誤差標(biāo)準(zhǔn)可以表明預(yù)測值與真實值之間的偏差,其計算公式如下:

    除了上述誤差標(biāo)準(zhǔn),本文同樣選擇了平均絕對誤差(mean absolute error,MAE)作為另一種誤差標(biāo)準(zhǔn),對預(yù)測結(jié)果進(jìn)行誤差計算。與RMSE不同, MAE主要表現(xiàn)的是預(yù)測值誤差的實際情況,其計算公式如下:

    3 算例分析

    3.1 實驗設(shè)計

    本文采用國內(nèi)某風(fēng)電場不同季節(jié)典型風(fēng)力發(fā)電功率作為數(shù)據(jù),數(shù)據(jù)采樣時間間隔為1 min。數(shù)據(jù)總體為2類:1)一天24 h的風(fēng)力發(fā)電功率數(shù)據(jù);2)一周168 h的風(fēng)力發(fā)電功率數(shù)據(jù)。通過2組不同規(guī)模的數(shù)據(jù),與ARIMA算法分別進(jìn)行預(yù)測精度對比,以驗證LSTM神經(jīng)網(wǎng)絡(luò)在不同時間尺度以及不同數(shù)據(jù)規(guī)模下的預(yù)測效果。

    為了優(yōu)化神經(jīng)網(wǎng)絡(luò)運算效率及精度,在所應(yīng)用的LSTM神經(jīng)網(wǎng)絡(luò)中,激活函數(shù)選擇線性整流函數(shù)(rectified linear unit,ReLU),優(yōu)化函數(shù)選擇自適應(yīng)矩估計(adaptive moment estimation,ADAM),網(wǎng)絡(luò)訓(xùn)練過程中的誤差函數(shù)選擇均方誤差(mean squared error,MSE)。

    3.2 預(yù)測性能分析

    分別采用ARIMA算法與LSTM神經(jīng)網(wǎng)絡(luò)進(jìn)行春季日風(fēng)力發(fā)電功率時序預(yù)測分析,并對所得的真實值與預(yù)測值進(jìn)行對比。針對一天1 440個風(fēng)力發(fā)電功率數(shù)據(jù),本文將0:00—16:00的864個數(shù)據(jù)作為訓(xùn)練集,16:00—24:00的576個數(shù)據(jù)作為測試集,開展了2種方法的預(yù)測效果分析。圖4為基于測試集數(shù)據(jù)得到的日風(fēng)力發(fā)電功率ARIMA、LSTM神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果,及應(yīng)用2種方法的預(yù)測殘差對比??梢钥闯?,LSTM神經(jīng)網(wǎng)絡(luò)的時序預(yù)測結(jié)果在整體趨勢上更接近風(fēng)力發(fā)電功率真實值,且不會產(chǎn)生較大偏差,相較于傳統(tǒng)ARIMA算法有更好的預(yù)測效果。

    為了驗證LSTM的預(yù)測能力,應(yīng)用時間尺度更長的春季周風(fēng)力發(fā)電功率曲線進(jìn)行時序預(yù)測分析。針對一周10080個風(fēng)力發(fā)電功率數(shù)據(jù),本文將初始約100h的6048個數(shù)據(jù)作為訓(xùn)練集,后續(xù)約68h的4032個數(shù)據(jù)作為測試集,開展了預(yù)測效果分析。圖5為基于測試集數(shù)據(jù)得到的周風(fēng)力發(fā)電功率ARIMA、LSTM神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果,及應(yīng)用2種方法的預(yù)測殘差對比??梢钥闯觯捎谟?xùn)練數(shù)據(jù)增加,LSTM神經(jīng)網(wǎng)絡(luò)與ARIMA算法的預(yù)測效果均有明顯提升;但LSTM神經(jīng)網(wǎng)絡(luò)的預(yù)測結(jié)果更符合原時序序列的趨勢,且不會出現(xiàn)較大波動,在長時間序列的周數(shù)據(jù)預(yù)測中計算效果更好。應(yīng)用式(14)、(15)分別對LSTM神經(jīng)網(wǎng)絡(luò)與ARIMA算法的預(yù)測結(jié)果進(jìn)行了對比,如表1所示。在春季日風(fēng)力發(fā)電功率曲線中,LSTM神經(jīng)網(wǎng)絡(luò)的2種預(yù)測誤差RMSE、MAE比ARIMA算法分別降低了7.7%、7%。在數(shù)據(jù)規(guī)模更大的春季周風(fēng)力發(fā)電功率曲線中,LSTM神經(jīng)網(wǎng)絡(luò)的2種預(yù)測誤差RMSE、MAE比ARIMA算法分別降低了4.3%、3.1%。隨著數(shù)據(jù)規(guī)模的增大,LSTM神經(jīng)網(wǎng)絡(luò)本身的預(yù)測性能也得到了提高。與春季日風(fēng)力發(fā)電功率LSTM神經(jīng)網(wǎng)絡(luò)的預(yù)測誤差相比,春季周風(fēng)力發(fā)電功率預(yù)測誤差RMSE降低了7.1%,MAE降低了12.5%。由此可知,在針對不同數(shù)據(jù)規(guī)模的風(fēng)力發(fā)電功率預(yù)測中,深度學(xué)習(xí)算法均有良好表現(xiàn),且隨著數(shù)據(jù)規(guī)模的提升,可獲得更好預(yù)測效果。

    表1 基于LSTM神經(jīng)網(wǎng)絡(luò)與ARIMA算法的春季數(shù)據(jù)預(yù)測結(jié)果對比

    3.3 四季風(fēng)力發(fā)電功率預(yù)測分析結(jié)果

    為了體現(xiàn)LSTM神經(jīng)網(wǎng)絡(luò)更優(yōu)的預(yù)測效果,針對同一風(fēng)電場的風(fēng)電數(shù)據(jù),選取另外3個季節(jié)的日風(fēng)力發(fā)電功率數(shù)據(jù)與周風(fēng)力發(fā)電功率數(shù)據(jù)進(jìn)行了對比分析,如表2所示??梢钥闯觯琇STM神經(jīng)網(wǎng)絡(luò)在各季節(jié)風(fēng)力發(fā)電功率的時序預(yù)測結(jié)果中,相較于ARIMA算法的預(yù)測精度都有明顯提升;同時隨著數(shù)據(jù)規(guī)模變大,LSTM神經(jīng)網(wǎng)絡(luò)可以獲得更好的預(yù)測結(jié)果,在短時數(shù)據(jù)基礎(chǔ)以及長時數(shù)據(jù)基礎(chǔ)中均有良好表現(xiàn)。

    表2 基于LSTM神經(jīng)網(wǎng)絡(luò)與ARIMA算法的風(fēng)力發(fā)電功率數(shù)據(jù)預(yù)測結(jié)果對比

    4 結(jié)論

    通過對LSTM神經(jīng)網(wǎng)絡(luò)進(jìn)行分析,提出了一種基于深度學(xué)習(xí)算法的風(fēng)力發(fā)電功率預(yù)測方法,以多種誤差精度標(biāo)準(zhǔn)作為參考,通過與傳統(tǒng)時序預(yù)測算法ARIMA進(jìn)行對比分析得到了以下結(jié)論:

    1)基于深度學(xué)習(xí)算法的風(fēng)力發(fā)電功率時序預(yù)測方法可較好擬合原風(fēng)力發(fā)電功率曲線,驗證了該方法的有效性;

    2)通過與傳統(tǒng)ARIMA時序預(yù)測算法進(jìn)行對比分析可知,LSTM神經(jīng)網(wǎng)絡(luò)的預(yù)測誤差相對更小,整體效果優(yōu)于ARIMA方法;

    3)通過對典型四季日風(fēng)力發(fā)電功率與周風(fēng)力發(fā)電功率的預(yù)測分析可知,在數(shù)據(jù)規(guī)模明顯提升的情況下,LSTM神經(jīng)網(wǎng)絡(luò)的預(yù)測效果優(yōu)于ARIMA算法。

    然而,若風(fēng)力發(fā)電功率發(fā)生連續(xù)突變時,深度學(xué)習(xí)算法的跟蹤變化能力有待進(jìn)一步完善。下一步的研究中將融合多種優(yōu)化算法,以提升深度學(xué)習(xí)算法的預(yù)測性能。

    [1] 陳國初,公維祥,馮兆紅.熵判別人工蜂群算法優(yōu)化的風(fēng)電功率[J].電力系統(tǒng)及其自動化學(xué)報,2017,29(7):41-45,67.

    [2] 國家能源局.2018年風(fēng)電并網(wǎng)運行情況[EB/ OL].[2019-01-28].http://www.nea.gov.cn/2019-01/28/c_137780779.htm.

    [3] 周云,嚴(yán)正,李乃湖,等.BESS平抑風(fēng)功率波動算法中多重SOC反饋策略設(shè)計[J].太陽能學(xué)報,2017,38(6):1459-1467.

    [4] 薛禹勝,郁琛,趙俊華,等.關(guān)于短期及超短期風(fēng)電功率預(yù)測的評述[J].電力系統(tǒng)自動化,2015,39(6):141-151.

    [5] 閆鶴鳴,李相俊,麻秀范,等.基于超短期風(fēng)電預(yù)測功率的儲能系統(tǒng)跟蹤風(fēng)電計劃出力控制方法[J].電網(wǎng)技術(shù),2015,39(2):433-439.

    [6] Liu C Y,Chao Q,Wei L L.Wind-storage coupling based on actual data and fuzzy control in multiple time scales for real-time rolling smoothing of fluctuation [J].Electric Power Automation Equipment,2015,35(2):35-41.

    [7] Fang S,Chiang H D.A high-accuracy wind power forecasting model[J].IEEE Transactions on Power Systems,2017,32(2):1589-1590.

    [8] Qian Z,Pei Y,Cao L X,et al.Review of wind power forecasting method[J].High Voltage Engineering,2016, 42(4):1047-1060.

    [9] Hu M Y,Hu Z J,Zhang M L,et al.Research on wind power forecasting method based on improved AdaBoost.RT and KELM algorithm[J].Power System Technology,2017,37(2):536-542.

    [10] 張俊,賀旭,陸春良,等.基于數(shù)值天氣預(yù)報的光伏功率短期預(yù)測分類組合算法[J].廣東電力,2019,32(6):55-60.

    [11] 楊正瓴,劉仍祥,李真真.基于季風(fēng)和大氣壓分布的我國風(fēng)電功率預(yù)測研究[J].分布式能源,2018,3(2):29-38.

    [12] 吳宇杭,溫步瀛,江岳文,等.基于變分模式分解和深度門控循環(huán)網(wǎng)絡(luò)的風(fēng)速短期預(yù)測模型[J].電網(wǎng)與清潔能源,2018,34(12):59-64.

    [13] 雷爽,王鵬卉,張亞剛.基于集合經(jīng)驗?zāi)B(tài)分解的組合風(fēng)速預(yù)測方法[J].電力科學(xué)與工程,2018,34(1):18-23.

    [14] 楊明莉,劉三明,王致杰,等.卡爾曼小波神經(jīng)網(wǎng)絡(luò)風(fēng)速預(yù)測[J].電力系統(tǒng)及其自動化學(xué)報,2015(12):42-46.

    [15] 方必武,劉滌塵,王波,等.基于小波變換和改進(jìn)螢火蟲算法優(yōu)化 LSSVM 的短期風(fēng)速預(yù)測[J].電力系統(tǒng)保護(hù)與控制,2016(8):37-43.

    [16] 李國慶,張鈺,張明江,等.基于 MRMR 的集合經(jīng)驗?zāi)B(tài)分解和支持向量機(jī)的風(fēng)電功率實時預(yù)測[J].東北電力大學(xué)學(xué)報,2017,37(2):39-44.

    [17] 曹俊波,周任軍,鄧學(xué)華,等.考慮優(yōu)化 ARIMA 模型差分次數(shù)的風(fēng)功率預(yù)測[J].電力系統(tǒng)及其自動化學(xué)報,2019,31(1):105-111.

    [18] 李彬,彭曙蓉,彭君哲,等.基于深度學(xué)習(xí)分位數(shù)回歸模型的風(fēng)電功率概率密度預(yù)測[J].電力自動化設(shè)備,2018,38(9):15-20.

    [19] 江岳春,張丙江,邢方方,等.基于混沌時間序列 GA-VNN 模型的超短期風(fēng)功率多步預(yù)測[J].電網(wǎng)技術(shù),2015,39(8):2160-2166.

    [20] 陳寧,薛禹勝,丁杰,等.利用空間相關(guān)性的超短期風(fēng)速預(yù)測[J].電力系統(tǒng)自動化,2017(12):124-130.

    [21] 陳杰,沈艷霞,陸欣,等.一種風(fēng)電功率概率區(qū)間多目標(biāo)智能優(yōu)化預(yù)測方法[J].電網(wǎng)技術(shù),2016,40(8):2281-2287.

    [22] 吳曉剛,孫榮富,喬穎,等.基于風(fēng)電場功率特性的風(fēng)電預(yù)測誤差分布估計[J].電網(wǎng)技術(shù),2017,41(6):1802-1807.

    [23] Li D,Ren Z Y,Yan W,et al.Month-ahead wind power curve probabilistic prediction based on factor analysis and quantile regression neural network[J].Proceedings of the CSEE,2017,37(18):5238-5248.

    [24] 董偉,楊強(qiáng),葛磊蛟,等.基于混合智能模型的分布式風(fēng)力發(fā)電預(yù)測方法[J].供用電,2019(1):34-39.

    [25] Li Y,Li X J,Jia X C,et al.Monitoring and control for hundreds megawatt scale battery energy storage station based on multi-agent:methodology and system design[C]//2018 IEEE International Conference of Safety Produce Informatization.Chongqing,China:IEEE,2018.

    [26] Li X J,Wang S X.A review on energy management,operation control and application methods for grid battery energy storage systems[J].CSEE Journal of Power and Energy Systems,DOI:10.17775/CSEEJPES,2019.00160.

    [27] 馮澤磊,吳美鳳.基于大數(shù)據(jù)平臺的火電機(jī)組實時運行狀態(tài)監(jiān)測[J].華電技術(shù),2019,41(4):29-32,37.

    [28] 李冰,張妍,劉石.基于LSTM的短期風(fēng)速預(yù)測研究[J].計算機(jī)仿真,2018,35(11):456-461.

    [29] 梁榮,王洪濤,吳奎華,等.基于神經(jīng)網(wǎng)絡(luò)和 ARIMA 模型的冷熱電短期負(fù)荷預(yù)測[J].電力系統(tǒng)及其自動化學(xué)報,DOI:10.19635/j.cnki.csu-epsa.000260.

    [30] 赫衛(wèi)國,郝向軍,郭雅娟,等.基于ARIMA和SVR的光伏電站超短期功率預(yù)測[J].廣東電力,2017,30(8):32-37.

    [31] 王旭,陳瀟一.基于AHP、ARIMA算法的電力信息系統(tǒng)負(fù)載預(yù)測研究與應(yīng)用[J].電網(wǎng)與清潔能源,2017,33(8):20-25.

    Wind Power Prediction Method Based on Long Short-term Memory Neural Network

    LI Xiangjun, XU Gejian

    (State Key Laboratory of Control and Operation of Renewable Energy and Storage Systems (China Electric Power Research Institute), Haidian District, Beijing 100192, China)

    Wind power generation process has strong randomness, which leads to low accuracy of wind power prediction. In view of the above phenomenon,a wind power generation power prediction method based on deep learning algorithm was proposed. Taking the historical wind power data as input, a wind power prediction model was established to realize the wind power prediction on a time scale in the future. The results of the example show that compared with the traditional time-series prediction method, the average absolute error of the wind power prediction results based on long short-term memory neural network is smaller in each index, which verifies the feasibility and effectiveness of the above method in wind power generation prediction, and improves the accuracy of wind power generation prediction.

    deep learning; time-series prediction; wind

    10.12096/j.2096-4528.pgt.19108

    國家電網(wǎng)公司科技項目(DG71-19-015)。

    Project Supported by Science and Technology Project of State Grid Corporation of China (DG71-19-015).power generation; long short-term memory (LSTM) neural network

    2019-07-16。

    李相俊(1979),男,博士,教授級高級工程師,研究方向為大規(guī)模儲能技術(shù)、新能源與分布式發(fā)電、電力系統(tǒng)運行與控制,li_xiangjun@126.com。

    李相俊

    (責(zé)任編輯 尚彩娟)

    猜你喜歡
    方法模型
    一半模型
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    學(xué)習(xí)方法
    可能是方法不對
    3D打印中的模型分割與打包
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    亚洲国产精品合色在线| 日韩亚洲欧美综合| 日本五十路高清| 啦啦啦韩国在线观看视频| 欧美成人一区二区免费高清观看| 97在线视频观看| 日产精品乱码卡一卡2卡三| 三级男女做爰猛烈吃奶摸视频| 精品久久久久久久久久久久久| 欧美xxxx性猛交bbbb| 国产精品久久视频播放| 久久精品影院6| 男人舔女人下体高潮全视频| 三级男女做爰猛烈吃奶摸视频| 国产在线男女| 国产色爽女视频免费观看| 亚洲人与动物交配视频| 亚洲av熟女| 亚洲伊人久久精品综合 | 久久久精品大字幕| 欧美xxxx黑人xx丫x性爽| 爱豆传媒免费全集在线观看| 又粗又爽又猛毛片免费看| www.av在线官网国产| 尤物成人国产欧美一区二区三区| 97热精品久久久久久| 色综合站精品国产| 国产精品一及| 只有这里有精品99| 欧美区成人在线视频| 久久久久免费精品人妻一区二区| 男人狂女人下面高潮的视频| 久久久色成人| 欧美zozozo另类| 偷拍熟女少妇极品色| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲欧美日韩东京热| 国产综合懂色| 乱系列少妇在线播放| 日本欧美国产在线视频| 丰满少妇做爰视频| 岛国在线免费视频观看| 亚洲国产欧美在线一区| 亚洲自拍偷在线| 永久网站在线| 黄色配什么色好看| 精品久久久噜噜| 日本三级黄在线观看| 一级黄片播放器| 99久久无色码亚洲精品果冻| 乱人视频在线观看| 高清日韩中文字幕在线| 自拍偷自拍亚洲精品老妇| 天堂中文最新版在线下载 | 国产成人freesex在线| 免费人成在线观看视频色| 啦啦啦啦在线视频资源| 国产人妻一区二区三区在| 99久国产av精品国产电影| 国产精品爽爽va在线观看网站| 亚洲综合色惰| 精品久久久久久久久亚洲| 亚洲精品日韩在线中文字幕| 精品国产一区二区三区久久久樱花 | 国产一级毛片七仙女欲春2| 国产伦在线观看视频一区| 国产av一区在线观看免费| 三级国产精品片| 久久国内精品自在自线图片| 中国国产av一级| 性插视频无遮挡在线免费观看| 国产午夜精品久久久久久一区二区三区| 老司机影院毛片| 天天一区二区日本电影三级| 国产午夜精品一二区理论片| 观看免费一级毛片| 久久精品国产亚洲av天美| 最近2019中文字幕mv第一页| 久久久久九九精品影院| 18禁在线播放成人免费| 九色成人免费人妻av| 亚洲美女视频黄频| 久久久成人免费电影| 日韩制服骚丝袜av| 美女国产视频在线观看| 国产 一区精品| 午夜福利在线在线| 国产av码专区亚洲av| 国产午夜精品论理片| 91精品一卡2卡3卡4卡| av专区在线播放| 亚洲精品一区蜜桃| 天堂影院成人在线观看| 女的被弄到高潮叫床怎么办| 99热这里只有是精品在线观看| 欧美高清性xxxxhd video| 99热这里只有是精品在线观看| 国产精品美女特级片免费视频播放器| 午夜老司机福利剧场| 午夜亚洲福利在线播放| 男女国产视频网站| 精品久久国产蜜桃| 亚洲av福利一区| 国产成人91sexporn| 免费播放大片免费观看视频在线观看 | 欧美另类亚洲清纯唯美| 三级国产精品片| av女优亚洲男人天堂| 在线免费十八禁| 日本欧美国产在线视频| 亚洲高清免费不卡视频| 一二三四中文在线观看免费高清| 国产黄片视频在线免费观看| 超碰av人人做人人爽久久| 国产精品一区二区三区四区久久| 精品一区二区免费观看| 国产伦精品一区二区三区视频9| 精品久久久久久久久亚洲| 99久久精品热视频| 男女下面进入的视频免费午夜| 国产爱豆传媒在线观看| 国产极品精品免费视频能看的| 夫妻性生交免费视频一级片| 精品少妇黑人巨大在线播放 | 91av网一区二区| 天堂√8在线中文| 国产免费视频播放在线视频 | 亚州av有码| 亚洲美女视频黄频| 久久精品夜夜夜夜夜久久蜜豆| 成年女人永久免费观看视频| 亚洲五月天丁香| 国产亚洲精品久久久com| 18禁裸乳无遮挡免费网站照片| 最近中文字幕高清免费大全6| 校园人妻丝袜中文字幕| 久久久色成人| 在线播放国产精品三级| 日韩一区二区三区影片| 国产淫片久久久久久久久| 国产精华一区二区三区| 性色avwww在线观看| 嫩草影院入口| 美女内射精品一级片tv| 亚洲aⅴ乱码一区二区在线播放| 免费大片18禁| 少妇裸体淫交视频免费看高清| 听说在线观看完整版免费高清| 亚洲精品国产成人久久av| 听说在线观看完整版免费高清| 我要看日韩黄色一级片| 最近视频中文字幕2019在线8| 激情 狠狠 欧美| 99九九线精品视频在线观看视频| 深爱激情五月婷婷| 偷拍熟女少妇极品色| 少妇裸体淫交视频免费看高清| 午夜亚洲福利在线播放| 69人妻影院| av在线播放精品| 99久久无色码亚洲精品果冻| 性插视频无遮挡在线免费观看| 91精品伊人久久大香线蕉| 亚洲精品影视一区二区三区av| 性插视频无遮挡在线免费观看| 三级毛片av免费| 国产亚洲av嫩草精品影院| 禁无遮挡网站| 99在线视频只有这里精品首页| 人妻少妇偷人精品九色| 99热6这里只有精品| 建设人人有责人人尽责人人享有的 | 亚洲成人av在线免费| 亚洲中文字幕日韩| 国产爱豆传媒在线观看| 91午夜精品亚洲一区二区三区| av在线观看视频网站免费| 最新中文字幕久久久久| 最近视频中文字幕2019在线8| 国产片特级美女逼逼视频| av线在线观看网站| 搞女人的毛片| 亚洲av中文av极速乱| 精品久久久久久久久av| 亚洲欧美清纯卡通| 中文字幕亚洲精品专区| 亚洲av成人精品一区久久| 久久久久久久久久久免费av| 国产三级在线视频| 国产精品精品国产色婷婷| 熟女人妻精品中文字幕| 久久精品熟女亚洲av麻豆精品 | 欧美xxxx性猛交bbbb| 精品久久久久久久久久久久久| 亚洲欧美日韩东京热| 中文欧美无线码| 舔av片在线| 国产精品一二三区在线看| 成人综合一区亚洲| 三级男女做爰猛烈吃奶摸视频| 亚洲成人精品中文字幕电影| 男女国产视频网站| 99在线视频只有这里精品首页| 又爽又黄a免费视频| 日韩在线高清观看一区二区三区| 亚洲欧美精品自产自拍| 国产成人91sexporn| 91精品伊人久久大香线蕉| 国产精品久久久久久久久免| 亚洲国产日韩欧美精品在线观看| 久久欧美精品欧美久久欧美| 国产精品三级大全| 午夜福利视频1000在线观看| 亚洲国产成人一精品久久久| 国产伦一二天堂av在线观看| 日本五十路高清| 亚洲国产欧美人成| 亚洲成人久久爱视频| 啦啦啦韩国在线观看视频| 久久精品91蜜桃| av免费在线看不卡| 丰满少妇做爰视频| 99久久人妻综合| 国语自产精品视频在线第100页| 2022亚洲国产成人精品| 国产精品爽爽va在线观看网站| 午夜激情福利司机影院| 天堂网av新在线| 国产精品久久久久久久电影| 日日撸夜夜添| 色播亚洲综合网| 欧美成人免费av一区二区三区| 亚洲丝袜综合中文字幕| 国产欧美日韩精品一区二区| 高清视频免费观看一区二区 | av福利片在线观看| 男女下面进入的视频免费午夜| 久久午夜福利片| 国产 一区 欧美 日韩| 国产老妇女一区| 亚洲欧美成人综合另类久久久 | 国产伦理片在线播放av一区| 午夜福利网站1000一区二区三区| 国产黄a三级三级三级人| 精品一区二区三区视频在线| 国产综合懂色| 精品熟女少妇av免费看| 日韩大片免费观看网站 | 直男gayav资源| 插逼视频在线观看| www.色视频.com| 国模一区二区三区四区视频| 久久久精品大字幕| 成人毛片a级毛片在线播放| 亚洲精品国产成人久久av| 青春草国产在线视频| 国产伦在线观看视频一区| 亚洲欧美日韩无卡精品| 观看美女的网站| 国产一级毛片在线| 免费av观看视频| 亚洲精品久久久久久婷婷小说 | 日本欧美国产在线视频| 真实男女啪啪啪动态图| 欧美日韩一区二区视频在线观看视频在线 | 国产激情偷乱视频一区二区| 高清视频免费观看一区二区 | 国产毛片a区久久久久| 永久网站在线| 免费观看a级毛片全部| 色吧在线观看| 青青草视频在线视频观看| 日本欧美国产在线视频| 国产精品蜜桃在线观看| 亚洲欧美中文字幕日韩二区| 人妻少妇偷人精品九色| 国产v大片淫在线免费观看| 亚洲三级黄色毛片| 国产精品一区二区三区四区免费观看| 三级国产精品欧美在线观看| 亚洲伊人久久精品综合 | 国产精品一二三区在线看| 国产亚洲精品av在线| 国产成人精品久久久久久| 久久国产乱子免费精品| 高清av免费在线| 欧美xxxx黑人xx丫x性爽| 国产精品一二三区在线看| 亚洲性久久影院| 哪个播放器可以免费观看大片| 久久久精品欧美日韩精品| 亚洲四区av| 少妇人妻精品综合一区二区| 精品人妻一区二区三区麻豆| kizo精华| 精品一区二区三区人妻视频| 成人国产麻豆网| 乱人视频在线观看| 免费观看人在逋| 亚洲欧美日韩高清专用| 国产av一区在线观看免费| 国产精品不卡视频一区二区| 久久精品91蜜桃| 26uuu在线亚洲综合色| 六月丁香七月| 久久久久久久亚洲中文字幕| 国产黄片美女视频| 国产极品天堂在线| 丰满人妻一区二区三区视频av| 免费观看人在逋| 国产美女午夜福利| 久久久久久久亚洲中文字幕| 国产精品一区二区三区四区免费观看| 午夜福利高清视频| 亚洲精品乱久久久久久| 日本免费a在线| 国产麻豆成人av免费视频| 精品人妻一区二区三区麻豆| 少妇熟女aⅴ在线视频| 精品人妻视频免费看| 一级二级三级毛片免费看| 国产人妻一区二区三区在| 国产 一区 欧美 日韩| 99久久九九国产精品国产免费| 只有这里有精品99| 国产黄a三级三级三级人| 久久久久网色| 青春草视频在线免费观看| av播播在线观看一区| 国产黄色小视频在线观看| 麻豆乱淫一区二区| 欧美高清成人免费视频www| 中文字幕av成人在线电影| 中文精品一卡2卡3卡4更新| 六月丁香七月| 直男gayav资源| 国产精品一区二区在线观看99 | 人妻系列 视频| 欧美高清成人免费视频www| 中文字幕av成人在线电影| 亚洲内射少妇av| 2021少妇久久久久久久久久久| 亚洲国产精品国产精品| 国内精品一区二区在线观看| 日韩中字成人| 69av精品久久久久久| 一级av片app| 国产午夜精品论理片| 亚洲乱码一区二区免费版| 国产高清国产精品国产三级 | 久久韩国三级中文字幕| 亚洲av不卡在线观看| 亚洲精华国产精华液的使用体验| 欧美精品一区二区大全| 亚洲av成人av| 国产精品一区www在线观看| 午夜精品一区二区三区免费看| 精品熟女少妇av免费看| 国产精品,欧美在线| 久久国内精品自在自线图片| 乱码一卡2卡4卡精品| 国产熟女欧美一区二区| 免费不卡的大黄色大毛片视频在线观看 | 99久久九九国产精品国产免费| 亚洲乱码一区二区免费版| 超碰av人人做人人爽久久| 桃色一区二区三区在线观看| 久久久久久久久久成人| 神马国产精品三级电影在线观看| 亚洲最大成人av| 国产亚洲一区二区精品| av在线亚洲专区| 国产真实伦视频高清在线观看| 久久综合国产亚洲精品| 免费看av在线观看网站| 国产免费男女视频| 亚洲精品成人久久久久久| 国产成人福利小说| 老司机影院成人| 波多野结衣巨乳人妻| www日本黄色视频网| 插阴视频在线观看视频| 中文资源天堂在线| 免费av毛片视频| 丰满少妇做爰视频| 啦啦啦观看免费观看视频高清| 男的添女的下面高潮视频| 精品久久久久久久久av| 亚洲最大成人手机在线| 少妇的逼水好多| 免费观看a级毛片全部| av在线播放精品| 亚洲av不卡在线观看| 伦理电影大哥的女人| 国产成人91sexporn| 亚洲av免费在线观看| 亚洲五月天丁香| 成人欧美大片| 97人妻精品一区二区三区麻豆| 久久午夜福利片| 欧美一级a爱片免费观看看| 2022亚洲国产成人精品| 欧美成人a在线观看| 国产一区有黄有色的免费视频 | 久热久热在线精品观看| 国产探花极品一区二区| 亚洲欧洲国产日韩| 大香蕉97超碰在线| 欧美激情国产日韩精品一区| 精品久久久久久电影网 | 日韩,欧美,国产一区二区三区 | 91精品国产九色| 免费av不卡在线播放| 国产成人免费观看mmmm| 狂野欧美白嫩少妇大欣赏| 搞女人的毛片| 久久99热这里只有精品18| 老司机影院成人| 日韩在线高清观看一区二区三区| 国产不卡一卡二| 性插视频无遮挡在线免费观看| 十八禁国产超污无遮挡网站| 国产爱豆传媒在线观看| 18禁在线播放成人免费| 精品久久久久久成人av| 成人亚洲欧美一区二区av| 精品午夜福利在线看| 日韩一区二区三区影片| www.av在线官网国产| 欧美性猛交黑人性爽| 18禁动态无遮挡网站| 亚洲色图av天堂| 久久这里只有精品中国| 亚洲成人久久爱视频| 欧美激情久久久久久爽电影| 国产精品一区二区三区四区久久| 久久久国产成人精品二区| 亚洲中文字幕日韩| 天美传媒精品一区二区| 老师上课跳d突然被开到最大视频| 午夜福利在线在线| 免费一级毛片在线播放高清视频| 成人欧美大片| 91aial.com中文字幕在线观看| 成人无遮挡网站| 国产午夜精品论理片| 色哟哟·www| 欧美高清性xxxxhd video| 身体一侧抽搐| 国产黄色小视频在线观看| 级片在线观看| 男人舔奶头视频| 亚洲精品影视一区二区三区av| 免费观看精品视频网站| 免费看a级黄色片| 99国产精品一区二区蜜桃av| 精品国内亚洲2022精品成人| 日本五十路高清| av国产久精品久网站免费入址| 成人欧美大片| 九九久久精品国产亚洲av麻豆| 国产麻豆成人av免费视频| 亚洲成色77777| 中文字幕熟女人妻在线| 国产又黄又爽又无遮挡在线| 成人三级黄色视频| 国产精品美女特级片免费视频播放器| 中文亚洲av片在线观看爽| 卡戴珊不雅视频在线播放| 高清视频免费观看一区二区 | 免费av毛片视频| 一卡2卡三卡四卡精品乱码亚洲| 国产精品福利在线免费观看| 亚洲欧美中文字幕日韩二区| 你懂的网址亚洲精品在线观看 | 亚洲色图av天堂| 欧美性感艳星| 十八禁国产超污无遮挡网站| 久久精品国产自在天天线| 国产精品一二三区在线看| 午夜福利网站1000一区二区三区| 免费搜索国产男女视频| 久久亚洲国产成人精品v| 国产高清视频在线观看网站| 别揉我奶头 嗯啊视频| 国产三级在线视频| 成人av在线播放网站| 精品国产露脸久久av麻豆 | www日本黄色视频网| 久99久视频精品免费| .国产精品久久| 少妇猛男粗大的猛烈进出视频 | 久久6这里有精品| kizo精华| 又粗又硬又长又爽又黄的视频| 亚洲国产精品国产精品| 黄片无遮挡物在线观看| 水蜜桃什么品种好| 国产午夜福利久久久久久| 美女内射精品一级片tv| 国产成人freesex在线| 久久久精品94久久精品| 欧美高清性xxxxhd video| 国产亚洲5aaaaa淫片| av女优亚洲男人天堂| 国产精品一区二区性色av| 日本欧美国产在线视频| 亚洲中文字幕一区二区三区有码在线看| 91久久精品国产一区二区成人| 亚洲欧美成人精品一区二区| 高清日韩中文字幕在线| 别揉我奶头 嗯啊视频| 亚洲国产精品sss在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲中文字幕一区二区三区有码在线看| 欧美成人a在线观看| 黄色日韩在线| 一二三四中文在线观看免费高清| 波多野结衣巨乳人妻| 中文精品一卡2卡3卡4更新| 日韩欧美国产在线观看| 在线观看一区二区三区| 韩国高清视频一区二区三区| 日韩欧美 国产精品| 春色校园在线视频观看| 国产精品久久久久久久电影| 在线天堂最新版资源| 国产免费视频播放在线视频 | 99久久九九国产精品国产免费| 女人十人毛片免费观看3o分钟| 黄色配什么色好看| 日日摸夜夜添夜夜爱| 成年版毛片免费区| 美女内射精品一级片tv| 精品不卡国产一区二区三区| 亚洲欧美清纯卡通| 亚洲欧美精品自产自拍| 精品一区二区三区视频在线| 久久热精品热| 秋霞在线观看毛片| 黄片无遮挡物在线观看| 亚洲欧美精品综合久久99| 国内揄拍国产精品人妻在线| 久久精品国产鲁丝片午夜精品| 日韩大片免费观看网站 | 麻豆久久精品国产亚洲av| 久99久视频精品免费| 欧美又色又爽又黄视频| 欧美丝袜亚洲另类| 精品人妻视频免费看| 特大巨黑吊av在线直播| 日韩制服骚丝袜av| 十八禁国产超污无遮挡网站| 变态另类丝袜制服| 国产亚洲精品av在线| 两个人视频免费观看高清| 国语自产精品视频在线第100页| 91久久精品国产一区二区成人| 人妻夜夜爽99麻豆av| 久久久久久久国产电影| a级毛片免费高清观看在线播放| 国产精品1区2区在线观看.| 少妇熟女欧美另类| 日韩欧美三级三区| 亚洲最大成人中文| 亚洲成人久久爱视频| 日韩制服骚丝袜av| 久久久久免费精品人妻一区二区| 毛片女人毛片| 色噜噜av男人的天堂激情| 麻豆一二三区av精品| 成人无遮挡网站| 村上凉子中文字幕在线| 青春草视频在线免费观看| 免费搜索国产男女视频| 国产极品天堂在线| 又黄又爽又刺激的免费视频.| 精品免费久久久久久久清纯| 91久久精品电影网| 直男gayav资源| 一级二级三级毛片免费看| 国产一级毛片七仙女欲春2| 亚洲一区高清亚洲精品| 草草在线视频免费看| 欧美变态另类bdsm刘玥| 男人的好看免费观看在线视频| 美女大奶头视频| 日韩一区二区视频免费看| 一边摸一边抽搐一进一小说| 深夜a级毛片| 久久这里有精品视频免费| 欧美成人一区二区免费高清观看| 91aial.com中文字幕在线观看| 联通29元200g的流量卡| 91精品伊人久久大香线蕉| 亚洲真实伦在线观看| 夜夜爽夜夜爽视频| 国产免费视频播放在线视频 | 男女国产视频网站| 日本黄色视频三级网站网址| 亚洲国产成人一精品久久久| 淫秽高清视频在线观看| 日本-黄色视频高清免费观看| 中文字幕制服av| 中文字幕熟女人妻在线| 日本黄大片高清| 别揉我奶头 嗯啊视频| 在线观看美女被高潮喷水网站| 成人毛片a级毛片在线播放| 在线播放国产精品三级| 一夜夜www| 男人舔奶头视频| 欧美一级a爱片免费观看看| 亚洲高清免费不卡视频| 精品人妻视频免费看|