999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于VMD-TCN 模型的渭河流域月徑流量預(yù)測(cè)研究

2023-10-07 18:36:19張上要羅軍剛石國(guó)棟景鑫連亞妮左崗崗
人民黃河 2023年10期

張上要 羅軍剛 石國(guó)棟 景鑫 連亞妮左崗崗

摘 要:為了有效減小徑流序列非線性、非平穩(wěn)性對(duì)徑流預(yù)測(cè)精度的影響,提出一種將變分模態(tài)分解(VMD)與時(shí)間卷積網(wǎng)絡(luò)(TCN)相耦合的VMD-TCN 模型,并用于渭河流域咸陽(yáng)和華縣水文站的月徑流量預(yù)測(cè)。采用多種模型性能評(píng)價(jià)指標(biāo)對(duì)VMD-TCN 模型與其他模型(TCN、EEMD-TCN、ARIMA)的預(yù)測(cè)性能進(jìn)行比較,并分析了VMD-TCN 模型在不同預(yù)見(jiàn)期下的預(yù)測(cè)表現(xiàn)。結(jié)果顯示VMDTCN模型較其他模型具有預(yù)測(cè)誤差更小、預(yù)測(cè)精度更高以及峰谷值擬合更優(yōu)的特點(diǎn),且隨著預(yù)見(jiàn)期的增大,模型的預(yù)測(cè)性能會(huì)逐漸降低但預(yù)測(cè)效果保持良好。

關(guān)鍵詞:徑流預(yù)測(cè);變分模態(tài)分解;時(shí)間卷積網(wǎng)絡(luò);渭河流域

中圖分類號(hào):P338 文獻(xiàn)標(biāo)志碼:A doi:10.3969/ j.issn.1000-1379.2023.10.005

引用格式:張上要,羅軍剛,石國(guó)棟,等.基于VMD-TCN 模型的渭河流域月徑流量預(yù)測(cè)研究[J].人民黃河,2023,45(10):25-29.

在變化環(huán)境影響下,徑流變化逐漸呈現(xiàn)出非線性、非平穩(wěn)性特征,傳統(tǒng)的線性時(shí)間序列模型預(yù)測(cè)精度不高且模型方法的改進(jìn)空間有限[1-2] 。隨著人工智能的發(fā)展,機(jī)器學(xué)習(xí)領(lǐng)域的非線性模型開(kāi)始應(yīng)用于徑流預(yù)測(cè),并且取得不錯(cuò)的效果。同時(shí),信號(hào)分解技術(shù)能夠?qū)⑺臅r(shí)間序列分解成若干相對(duì)穩(wěn)定的固有模態(tài)分量(IMFs),將其與徑流預(yù)測(cè)模型相耦合能夠顯著提升徑流預(yù)測(cè)的精度[3] 。為克服徑流序列非平穩(wěn)性導(dǎo)致預(yù)測(cè)精度差的問(wèn)題,周婷等[4] 將小波分解(WD)與支持向量機(jī)(SVR)進(jìn)行耦合并應(yīng)用于徑流預(yù)測(cè),研究表明WD 能夠有效提升徑流預(yù)測(cè)的精度。張金萍等[5] 將自回歸滑動(dòng)平均模型(ARMA)與完全集合經(jīng)驗(yàn)?zāi)B(tài)分解(CEEMDAN) 結(jié)合,發(fā)現(xiàn)CEEMDAN 能夠顯著提高ARMA 的擬合優(yōu)度。桑宇婷等[6] 將CEEMDAN 與BP神經(jīng)網(wǎng)絡(luò)進(jìn)行耦合,并用于汾河流域的月徑流量預(yù)測(cè),預(yù)測(cè)結(jié)果顯示CEEMDAN-BP 模型在模擬期和驗(yàn)證期的徑流預(yù)測(cè)精度都達(dá)到甲級(jí)。包苑村等[7] 將變分模態(tài)分解(VMD)與卷積-長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)(CNNLSTM)結(jié)合開(kāi)展月徑流量預(yù)測(cè),結(jié)果表明VMD-CNNLSTM模型對(duì)徑流序列的峰值和谷值擬合更優(yōu)。He等[8] 將VMD 與深度神經(jīng)網(wǎng)絡(luò)(DNN)進(jìn)行耦合,用于渭河流域張家山水文站的日徑流量預(yù)測(cè),其預(yù)測(cè)結(jié)果表現(xiàn)良好。Zuo 等[9] 采用一種兩階段分解集成水文預(yù)測(cè)方法,解決信號(hào)分解過(guò)程中未來(lái)信息的引入導(dǎo)致模型預(yù)測(cè)結(jié)果失真的問(wèn)題。

鑒于上述研究背景,筆者將變分模態(tài)分解(VMD)與時(shí)間卷積網(wǎng)絡(luò)(TCN)相結(jié)合構(gòu)建了VMD-TCN 耦合預(yù)測(cè)模型,將其應(yīng)用于渭河流域咸陽(yáng)和華縣水文站的月徑流量預(yù)測(cè)。通過(guò)與其他3 種模型(TCN、EEMD-TCN、ARIMA)對(duì)比,驗(yàn)證VMD-TCN 的預(yù)測(cè)性能,同時(shí)分析VMD-TCN 在不同預(yù)見(jiàn)期的預(yù)測(cè)表現(xiàn)。

1.2 時(shí)間卷積網(wǎng)絡(luò)(TCN)

時(shí)間卷積網(wǎng)絡(luò)(TCN)是一種新型的時(shí)間序列預(yù)測(cè)模型,它是在卷積神經(jīng)網(wǎng)絡(luò)(CNN)的基礎(chǔ)上進(jìn)行改進(jìn)的,其模型架構(gòu)包括因果卷積(causal convolution)、空洞卷積(dilated convolution)和殘差模塊(residual block)。

TCN 的梯度穩(wěn)定使其能夠有效避免模型訓(xùn)練過(guò)程中發(fā)生梯度消失或梯度爆炸導(dǎo)致模型訓(xùn)練失敗的問(wèn)題,同時(shí)具有靈活的感受野[13] 。

TCN 模型卷積結(jié)構(gòu)的層與層之間都是相互依賴的,后一層的輸出依賴上一層的輸入,每一層都間隔神經(jīng)元對(duì)上一層進(jìn)行信息提取,通過(guò)間隔采樣的形式,逐層膨脹系數(shù)以2 為指數(shù)增長(zhǎng),以便使用更少的層數(shù)獲得更大的感受野。為保證每一層的信息不丟失,需要對(duì)每一層的邊緣進(jìn)行填充。TCN 模型梯度穩(wěn)定的原因是引入殘差模塊,殘差模塊的引入可以保證在進(jìn)行特征提取時(shí)信息相對(duì)完整,保證模型預(yù)測(cè)準(zhǔn)確[14] 。

1.3 VMD-TCN 預(yù)測(cè)模型構(gòu)建

本研究將VMD 分解與TCN 模型相結(jié)合構(gòu)建VMD-TCN 模型用于月徑流量預(yù)測(cè),使用貝葉斯優(yōu)化算法進(jìn)行超參數(shù)調(diào)優(yōu)。開(kāi)展徑流預(yù)測(cè)主要分3 個(gè)階段,分別是數(shù)據(jù)集劃分階段、預(yù)測(cè)樣本生成階段和徑流預(yù)測(cè)階段,具體流程見(jiàn)圖1。VMD-TCN 模型在進(jìn)行序列分解時(shí)將驗(yàn)證集數(shù)據(jù)逐一添加至訓(xùn)練集中能夠有效避免模型訓(xùn)練引入未來(lái)信息。同時(shí)采用直接法生成預(yù)測(cè)樣本能夠提高模型的預(yù)測(cè)效率,節(jié)省計(jì)算資源[9] 。

1.4 模型性能評(píng)價(jià)指標(biāo)

本文使用4 種模型性能評(píng)價(jià)指標(biāo)用于評(píng)估模型的預(yù)測(cè)能力,分別是納什效率系數(shù)(NSE)、標(biāo)準(zhǔn)均方根誤差(NRMSE)、平均絕對(duì)百分比誤差(MAPE)和峰值預(yù)測(cè)性能指標(biāo)(PPTS)。

NSE 反映預(yù)測(cè)結(jié)果所含信息量占預(yù)測(cè)目標(biāo)所含信息量的比重,其值越接近1 表示模型的預(yù)測(cè)效果越好。NRMSE 是對(duì)預(yù)測(cè)結(jié)果均方根誤差的歸一化,適用于多模型對(duì)比, 其值越接近0 則預(yù)測(cè)效果越好。MAPE 直接反映模型的預(yù)測(cè)誤差,越接近0 表示模型的預(yù)測(cè)性能越好。PPTS 能夠評(píng)價(jià)不同比例峰值的預(yù)測(cè)效果,越接近0 表示峰值預(yù)測(cè)性能越好。按照《水文情報(bào)預(yù)報(bào)規(guī)范》(GB/ T 22482—2008)規(guī)定,NSE≥0.9,預(yù)測(cè)精度為甲級(jí);0.7≤NSE<0.9,預(yù)測(cè)精度為乙級(jí);0.5≤NSE<0.7,預(yù)測(cè)精度為丙級(jí)。

2 實(shí)例分析

2.1 試驗(yàn)環(huán)境搭建

本試驗(yàn)采用python3.7 環(huán)境,VMD 的附加分解使用python 的vmdpy 庫(kù)進(jìn)行,TCN 模型采用keras2.8.0搭建,超參數(shù)優(yōu)化框架采用keras-tuner1.1.0。

2.2 研究區(qū)域概況及數(shù)據(jù)資料

渭河流域面積為134 766 km2,其中陜西省境內(nèi)流域面積為67 108 km2。渭河干流全長(zhǎng)818 km,是黃河最大的支流,多年平均徑流量高達(dá)75.7 億m3,多年平均降水量550 mm。咸陽(yáng)水文站控制面積為4.7 萬(wàn)km2,多年平均徑流量為37.298 億m3,是渭河中游的控制站。華縣水文站控制面積為10.6 萬(wàn)km2,多年平均徑流量為63.768 億m3,是渭河下游主要控制站。本研究選用咸陽(yáng)水文站和華縣水文站1953 年1 月至2018 年12 月共66 a 的月徑流量數(shù)據(jù),按7 ∶ 3 的比例劃分訓(xùn)練集和驗(yàn)證集。

2.3 月徑流量數(shù)據(jù)預(yù)處理

2.3.1 月徑流量序列VMD 分解

采用附加分解方式對(duì)月徑流量序列進(jìn)行VMD 分解,分別將咸陽(yáng)站和華縣站的訓(xùn)練數(shù)據(jù)進(jìn)行分解。將兩個(gè)水文站的訓(xùn)練數(shù)據(jù)分別進(jìn)行VMD 分解預(yù)試驗(yàn),兩站在k =8 時(shí)分解效果最好,且恰好不出現(xiàn)模態(tài)混疊現(xiàn)象,圖2 為咸陽(yáng)站訓(xùn)練數(shù)據(jù)分解之后的IMF1 與IMF8 圖,IMF2 ~IMF7 圖略,華縣站同理。

2.3.2 特征序列滯后長(zhǎng)度確定

采用直接法生成徑流預(yù)測(cè)樣本前,需要確定每一個(gè)特征序列輸入的滯后長(zhǎng)度,以確定t 時(shí)刻的流量與其前期流量的相關(guān)性[15] 。考慮到完整原始徑流序列存在以12 個(gè)月為周期的規(guī)律性,將未分解的TCN 模型特征序列輸入的滯后長(zhǎng)度設(shè)為12。采用自相關(guān)函數(shù)法(ACF)和偏自相關(guān)函數(shù)法(PACF)確定VMD 分解后各個(gè)模態(tài)分量的滯后長(zhǎng)度,具體方法見(jiàn)文獻(xiàn)[16-17]。

2.4 VMD-TCN 預(yù)測(cè)結(jié)果分析

2.4.1 不同模型預(yù)測(cè)結(jié)果對(duì)比分析

徑流序列分解并確定特征序列滯后長(zhǎng)度后可構(gòu)建徑流預(yù)測(cè)樣本,將預(yù)測(cè)樣本輸入TCN 模型可開(kāi)展徑流預(yù)測(cè)。為充分驗(yàn)證模型的預(yù)測(cè)性能,將TCN、EEMDTCN、ARIMA 與VMD-TCN 進(jìn)行對(duì)比,其中EEMDTCN模型同樣使用附加分解以避免未來(lái)信息引入。圖3、圖4 分別為咸陽(yáng)站預(yù)見(jiàn)期1 個(gè)月的多模型實(shí)測(cè)與預(yù)測(cè)徑流量折線圖和散點(diǎn)分布圖,華縣站同理。表1 為兩站多模型預(yù)測(cè)性能評(píng)價(jià)指標(biāo)對(duì)比。

以咸陽(yáng)站為例,由圖3、圖4 可知:1)相對(duì)TCN、EEMD-TCN 和ARIMA,VMD-TCN 對(duì)實(shí)測(cè)序列的擬合效果最好,特別是峰值的擬合,散點(diǎn)分布更加集中,理想擬合線與線性擬合線的夾角最小,說(shuō)明VMD-TCN 模型的預(yù)測(cè)值與實(shí)測(cè)值具有高度的一致性。2)EEMD-TCN模型的擬合效果次之,但較VMD-TCN 對(duì)峰值的擬合不佳,散點(diǎn)更為分散。3)單一的TCN 和ARIMA 模型擬合效果不佳,TCN 與ARIMA 模型對(duì)峰值的擬合有明顯的滯后且難以預(yù)測(cè)高峰值,說(shuō)明單一TCN 和ARIMA 難以預(yù)測(cè)峰值的到達(dá)時(shí)間及量值。

結(jié)合表1 咸陽(yáng)站的預(yù)測(cè)性能評(píng)價(jià)結(jié)果可知:1)VMDTCN模型的預(yù)測(cè)效果最好,NSE 在0.9 以上,預(yù)測(cè)精度達(dá)到了甲級(jí)。VMD-TCN 模型與單一的TCN 模型相比,NSE 提升了318%,NRMSE、MAPE 和PPTS 分別下降了67%、40%、85%。采用VMD 分解之后的TCN 模型預(yù)測(cè)性能有了極大的提升,主要原因是VMD 分解能夠得到平穩(wěn)的信號(hào)分量,從而有效降低徑流序列中噪聲對(duì)預(yù)測(cè)的影響,進(jìn)而提升預(yù)測(cè)的精度。2)EEMDTCN模型相比TCN,NSE 提高了205%,NRMSE、MAPE和PPTS 分別下降了35%、26%、55%,說(shuō)明EEMD 分解能對(duì)非平穩(wěn)時(shí)間序列起到一定的降噪作用并提升預(yù)測(cè)精度。采用附加分解的EEMD-TCN 的NSE 只有0.67,預(yù)測(cè)精度只有丙級(jí),主要原因是EEMD 在分解過(guò)程中避免了未來(lái)信息的引入,預(yù)測(cè)精度有所降低但更可靠,同時(shí)EEMD 分解會(huì)出現(xiàn)模態(tài)混疊現(xiàn)象,較VMD分解效果更差。3)使用分解之后的TCN 預(yù)測(cè)效果遠(yuǎn)高于TCN 和ARIMA,由PPTS 指標(biāo)可以看出使用信號(hào)分解的模型峰值的預(yù)測(cè)效果更好。分析華縣站預(yù)測(cè)結(jié)果可得到與上述相同的結(jié)論。

總體上,4 種月徑流量預(yù)測(cè)模型預(yù)測(cè)性能的高低排序如下:VMD-TCN>EEMD-TCN>TCN≈ARIMA。耦合信號(hào)分解算法的組合模型預(yù)測(cè)效果高于單一模型,上述分析結(jié)果說(shuō)明了VMD-TCN 對(duì)復(fù)雜非線性的徑流序列的處理和預(yù)測(cè)是可行的。

2.4.2 不同預(yù)見(jiàn)期預(yù)測(cè)結(jié)果分析

為進(jìn)一步驗(yàn)證VMD-TCN 模型在不同預(yù)見(jiàn)期下的預(yù)測(cè)性能,以華縣站為例,針對(duì)預(yù)見(jiàn)期分別為1、3、5、7個(gè)月開(kāi)展徑流預(yù)測(cè)研究。圖5 為華縣站VMD-TCN 模型在預(yù)見(jiàn)期為1、5 個(gè)月的預(yù)測(cè)與實(shí)測(cè)徑流量折線圖及散點(diǎn)分布圖,表2 為不同預(yù)見(jiàn)期預(yù)測(cè)性能評(píng)價(jià)指標(biāo)對(duì)比。

不同預(yù)見(jiàn)期下,模型的預(yù)測(cè)能力會(huì)有所變化,由圖5 可知:1)預(yù)見(jiàn)期為1 個(gè)月時(shí),實(shí)測(cè)值折線與預(yù)測(cè)值折線能夠很好地貼合,對(duì)預(yù)測(cè)序列的峰值和谷值的擬合效果較好。從散點(diǎn)圖也可以看出,線性擬合線與理想擬合線的夾角較小,說(shuō)明預(yù)見(jiàn)期為1 個(gè)月時(shí)實(shí)測(cè)值與預(yù)測(cè)值的一致性較高。2)隨著預(yù)見(jiàn)期增大,實(shí)測(cè)值與預(yù)測(cè)值之間的擬合誤差增大,線性擬合線和理想擬合線的角度增大,說(shuō)明VMD-TCN 模型的預(yù)測(cè)能力會(huì)隨著預(yù)見(jiàn)期的增大而降低。

結(jié)合表2 不同預(yù)見(jiàn)期的預(yù)測(cè)性能評(píng)價(jià)指標(biāo)可知:1)預(yù)見(jiàn)期為1 個(gè)月時(shí),NSE 在0.9 以上,預(yù)見(jiàn)期不斷增大NSE 值逐漸減小,當(dāng)預(yù)見(jiàn)期為7 個(gè)月時(shí)NSE 接近0.8,說(shuō)明隨著預(yù)見(jiàn)期增大,VMD-TCN 模型的預(yù)測(cè)精度會(huì)下降,但依舊保持著較好的預(yù)測(cè)能力。2)隨著預(yù)見(jiàn)期增大,NSE 逐漸減小而NRMSE、MAPE 和PPTS 總體上逐漸增大,說(shuō)明預(yù)見(jiàn)期增大而預(yù)測(cè)精度降低,部分原因是模型的峰值預(yù)測(cè)能力下降。預(yù)見(jiàn)期增大而模型預(yù)測(cè)能力降低的原因是隨著預(yù)見(jiàn)期的增大,預(yù)測(cè)值與目標(biāo)值的相關(guān)性逐漸降低,模型無(wú)法獲得更為準(zhǔn)確的信息進(jìn)行精確預(yù)測(cè)。

綜上所述,VMD-TCN 模型的預(yù)測(cè)性能會(huì)隨著預(yù)見(jiàn)期的增大而降低,特別是在峰值的預(yù)測(cè)上。預(yù)見(jiàn)期為1 個(gè)月時(shí),預(yù)測(cè)效果最好,預(yù)測(cè)精度達(dá)到甲級(jí);預(yù)見(jiàn)期為7 個(gè)月時(shí)效果較差,但預(yù)測(cè)精度仍達(dá)到了乙級(jí)。可見(jiàn),VMD-TCN 在不同預(yù)見(jiàn)期的徑流預(yù)測(cè)也是可行且有效的。

3 結(jié)論

1)VMD 對(duì)徑流序列進(jìn)行預(yù)處理并耦合TCN 模型能夠顯著提升月徑流量預(yù)測(cè)的精度,同時(shí)VMD 采用附加分解方式能夠避免模型預(yù)測(cè)過(guò)程中未來(lái)信息的引入,預(yù)測(cè)結(jié)果更為可靠。

2)VMD-TCN 與EEMD-TCN 較之單一的TCN 和ARIMA 模型的預(yù)測(cè)表現(xiàn)更佳,且EEMD 分解會(huì)出現(xiàn)影響預(yù)測(cè)精度的模態(tài)混疊問(wèn)題,VMD 分解能夠通過(guò)預(yù)試驗(yàn)確定模態(tài)分量的個(gè)數(shù)從而有效避免模態(tài)混疊問(wèn)題,VMD-TCN 預(yù)測(cè)精度高于EEMD-TCN。

3)隨著預(yù)見(jiàn)期的增大,構(gòu)建的徑流預(yù)測(cè)樣本所包含的信息會(huì)有所損失,預(yù)測(cè)因子與預(yù)測(cè)目標(biāo)的相關(guān)性逐漸降低,VMD-TCN 模型的預(yù)測(cè)性能會(huì)逐漸下降。

參考文獻(xiàn):

[1] 梁浩,黃生志,孟二浩,等.基于多種混合模型的徑流預(yù)測(cè)研究[J].水利學(xué)報(bào),2020,51(1):112-125.

[2] 雷曉輝,王浩,廖衛(wèi)紅,等.變化環(huán)境下氣象水文預(yù)報(bào)研究進(jìn)展[J].水利學(xué)報(bào),2018,49(1):9-18.

[3] 李文武,石強(qiáng),王凱,等.基于變分模態(tài)分解和深度門控網(wǎng)絡(luò)的徑流預(yù)測(cè)[J].水力發(fā)電學(xué)報(bào),2020,39(3):34-44.

[4] 周婷,金菊良,李榮波,等.基于小波支持向量機(jī)的徑流預(yù)測(cè)性能優(yōu)化分析[J].水力發(fā)電學(xué)報(bào),2017,36(10):45-55.

[5] 張金萍,許敏,張?chǎng)危龋冢茫牛牛停模粒危粒遥停?模型的年徑流量預(yù)測(cè)研究[J].人民黃河,2021,43(1):35-39.

[6] 桑宇婷,趙雪花,祝雪萍,等.基于CEEMD-BP 模型的汾河上游月徑流預(yù)測(cè)[J].人民黃河,2019,41(8):1-5.

[7] 包苑村,解建倉(cāng),羅軍剛.基于VMD-CNN-LSTM 模型的渭河流域月徑流預(yù)測(cè)[J].西安理工大學(xué)學(xué)報(bào),2021,37(1):1-8.

[8] HE X X,LUO J G,ZUO G G,et al.Daily Runoff ForecastingUsing a Hybrid Model Based on Variational Mode Decompo?sition and Deep Neural Networks[J].Water Resources Man?agement,2019,33(4):1571-1590.

[9] ZUO G G,LUO J G,WANG N,et al.Two?Stage VariationalMode Decomposition and Support Vector Regression for Stre?amflow Forecasting[J].Hydrology and Earth System Sciences,2020,24(11):5491-5518.

[10] DRAGOMIRETSKIY K,ZOSSO D.Variational Mode Decom?position[J].IEEE Transactions on Signal Processing,2013,62(3):531-544.

[11] 蘇麗娜,張亞?wèn)|,王梓丞,等.基于變分模態(tài)分解和能量譜的軌道電路分路不良故障監(jiān)測(cè)[J].鐵道標(biāo)準(zhǔn)設(shè)計(jì),2019,63(2):155-160.

[12] 方江雄,溫志平,顧華奇,等.基于變分模態(tài)分解的地震隨機(jī)噪聲壓制方法[J].石油地球物理勘探,2019,54(4):757-767,722.

[13] 孫蓉,李強(qiáng),羅海峰,等.基于自適應(yīng)移動(dòng)平滑與時(shí)間卷積網(wǎng)絡(luò)誤差修正的風(fēng)電功率預(yù)測(cè)[J].全球能源互聯(lián)網(wǎng),2022,5(1):11-22.

[14] 楊海晶,孫運(yùn)全,朱偉,等.基于CEEMD-TCN 模型的變壓器油中溶解氣體濃度預(yù)測(cè)方法[J].電子器件,2021,44(4):887-892.

[15] 王文川,高暢,徐雷.基于TVF-EMD 與LSTM 神經(jīng)網(wǎng)絡(luò)耦合的月徑流預(yù)測(cè)研究[J].中國(guó)農(nóng)村水利水電,2022(2):76-81,89.

[16] HE Z B,WEN X H,LIU H,et al.A Comparative Study ofArtificial Neural Network,Adaptive Neuro Fuzzy InferenceSystem and Support Vector Machine for Forecasting RiverFlow in the Semiarid Mountain Region[J].Journal of Hy? drology,2014,509:379-386.

[17] YU X,ZHANG X Q,QIN H.A Data?Driven Model Based onFourier Transform and Support Vector Regression forMonthly Reservoir Inflow Forecasting[J].Journal of Hydro?Environment Research,2018,18:12-24.

【責(zé)任編輯 張 帥】

主站蜘蛛池模板: 久青草免费在线视频| 婷婷午夜影院| 日韩高清中文字幕| 国产一二三区视频| 中文字幕无线码一区| 青草视频网站在线观看| 日本人真淫视频一区二区三区| 喷潮白浆直流在线播放| 亚洲综合经典在线一区二区| 精品国产自在现线看久久| 日本成人一区| 欧美精品亚洲精品日韩专区| 久久99精品国产麻豆宅宅| 成人午夜久久| 国产精品成人观看视频国产| 久久精品日日躁夜夜躁欧美| 1024你懂的国产精品| 亚洲av无码专区久久蜜芽| 91精品国产综合久久香蕉922 | 欧美亚洲网| 婷婷六月在线| 亚洲第一色视频| 亚洲国产亚综合在线区| a在线观看免费| www.99在线观看| 无码人妻免费| 欧美一级高清免费a| 成年人国产网站| 久久成人免费| 中文字幕亚洲乱码熟女1区2区| 国产麻豆永久视频| 亚洲天堂777| 老色鬼久久亚洲AV综合| 99这里精品| 国产自在线播放| 18禁不卡免费网站| 无码中文字幕精品推荐| 国产毛片片精品天天看视频| 91视频99| 在线亚洲小视频| 韩日免费小视频| 日本在线国产| 99在线视频精品| 538国产视频| 亚洲视频黄| 国产特级毛片aaaaaa| 成人国产精品网站在线看| 岛国精品一区免费视频在线观看| 欧美午夜小视频| 日韩毛片视频| 亚洲人成成无码网WWW| 国产视频一二三区| 亚洲国产av无码综合原创国产| 色悠久久久| 大乳丰满人妻中文字幕日本| 人妻中文久热无码丝袜| www中文字幕在线观看| 日韩欧美视频第一区在线观看| 免费在线a视频| 亚洲国产AV无码综合原创| 日本午夜网站| 国外欧美一区另类中文字幕| 欧美久久网| 日本91视频| 国产第一色| 国产对白刺激真实精品91| 国产91高清视频| 激情午夜婷婷| 欧美日韩精品一区二区视频| 国产精品久久久免费视频| 久久成人国产精品免费软件| 久久这里只有精品国产99| 最新日韩AV网址在线观看| 国产一区二区免费播放| 国产微拍一区二区三区四区| yjizz国产在线视频网| 91偷拍一区| 国产乱子伦手机在线| 四虎精品国产AV二区| 日韩无码黄色网站| 亚洲欧美另类视频| 欧美一级夜夜爽|