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

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

    2023-10-07 18:36:19張上要羅軍剛石國棟景鑫連亞妮左崗崗
    人民黃河 2023年10期
    關(guān)鍵詞:渭河流域

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

    2 實(shí)例分析

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

    本試驗(yàn)采用python3.7 環(huán)境,VMD 的附加分解使用python 的vmdpy 庫進(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。渭河干流全長818 km,是黃河最大的支流,多年平均徑流量高達(dá)75.7 億m3,多年平均降水量550 mm。咸陽水文站控制面積為4.7 萬km2,多年平均徑流量為37.298 億m3,是渭河中游的控制站。華縣水文站控制面積為10.6 萬km2,多年平均徑流量為63.768 億m3,是渭河下游主要控制站。本研究選用咸陽水文站和華縣水文站1953 年1 月至2018 年12 月共66 a 的月徑流量數(shù)據(jù),按7 ∶ 3 的比例劃分訓(xùn)練集和驗(yàn)證集。

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

    2.3.1 月徑流量序列VMD 分解

    采用附加分解方式對月徑流量序列進(jìn)行VMD 分解,分別將咸陽站和華縣站的訓(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 為咸陽站訓(xùn)練數(shù)據(jù)分解之后的IMF1 與IMF8 圖,IMF2 ~IMF7 圖略,華縣站同理。

    2.3.2 特征序列滯后長度確定

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

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

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

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

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

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

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

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

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

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

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

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

    3 結(jié)論

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

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

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

    參考文獻(xiàn):

    [1] 梁浩,黃生志,孟二浩,等.基于多種混合模型的徑流預(yù)測研究[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ù)測[J].水力發(fā)電學(xué)報(bào),2020,39(3):34-44.

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

    [5] 張金萍,許敏,張鑫,等.基于CEEMDAN-ARMA 模型的年徑流量預(yù)測研究[J].人民黃河,2021,43(1):35-39.

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

    [7] 包苑村,解建倉,羅軍剛.基于VMD-CNN-LSTM 模型的渭河流域月徑流預(yù)測[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] 蘇麗娜,張亞東,王梓丞,等.基于變分模態(tài)分解和能量譜的軌道電路分路不良故障監(jiān)測[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ù)測[J].全球能源互聯(lián)網(wǎng),2022,5(1):11-22.

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

    [15] 王文川,高暢,徐雷.基于TVF-EMD 與LSTM 神經(jīng)網(wǎng)絡(luò)耦合的月徑流預(yù)測研究[J].中國農(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é)任編輯 張 帥】

    猜你喜歡
    渭河流域
    渭河流域(陜西段)面源污染總磷時(shí)空特征評估
    渭河水岸帶土壤理化特征及其與重金屬含量的關(guān)系
    基于NDVI的渭河流域時(shí)空演變分析
    甘肅省渭河流域韭菜產(chǎn)業(yè)現(xiàn)狀及發(fā)展建議
    基于DPCA方法對渭河流域降水—徑流—蒸發(fā)序列的長程互相關(guān)性研究
    世界家苑(2018年1期)2018-04-27 11:42:06
    基于RS/GIS 渭河流域植被覆蓋時(shí)空變化特征研究
    燕太子回國
    環(huán)境流域治理中的法律對策
    法制博覽(2015年10期)2015-11-06 12:16:49
    渭河流域香菜夏秋無公害栽培技術(shù)
    生態(tài)系統(tǒng)服務(wù)支付意愿及其影響因素分析
    軟科學(xué)(2015年6期)2015-07-10 02:38:40
    久久综合国产亚洲精品| 男女边吃奶边做爰视频| 麻豆久久精品国产亚洲av| 网址你懂的国产日韩在线| 丝袜喷水一区| 亚洲人成网站高清观看| 国产成人a区在线观看| 久久精品国产亚洲av涩爱| 最近中文字幕2019免费版| 日韩人妻高清精品专区| 九九久久精品国产亚洲av麻豆| 亚洲欧美日韩无卡精品| 国产精品偷伦视频观看了| 欧美亚洲 丝袜 人妻 在线| 五月伊人婷婷丁香| 亚洲国产欧美在线一区| 一级二级三级毛片免费看| 少妇猛男粗大的猛烈进出视频 | 男人添女人高潮全过程视频| 在线观看人妻少妇| 青春草视频在线免费观看| 免费av观看视频| 亚洲在久久综合| 99re6热这里在线精品视频| 99re6热这里在线精品视频| 夫妻性生交免费视频一级片| 我的老师免费观看完整版| 国产高潮美女av| 欧美日韩在线观看h| 一级毛片 在线播放| 一边亲一边摸免费视频| 又粗又硬又长又爽又黄的视频| 精品一区在线观看国产| 九色成人免费人妻av| 免费大片18禁| 欧美日韩国产mv在线观看视频 | 亚洲av成人精品一区久久| 99精国产麻豆久久婷婷| 97超视频在线观看视频| 亚洲成人久久爱视频| 黄色欧美视频在线观看| 国产乱人视频| 久久影院123| 国产乱来视频区| 亚洲精华国产精华液的使用体验| 国产精品成人在线| 一级二级三级毛片免费看| 国产精品国产三级国产专区5o| 看黄色毛片网站| 久久精品人妻少妇| 亚洲精品国产色婷婷电影| 日韩 亚洲 欧美在线| 成人欧美大片| 欧美精品一区二区大全| 国产亚洲午夜精品一区二区久久 | 久久国产乱子免费精品| 性插视频无遮挡在线免费观看| 熟妇人妻不卡中文字幕| 黄色日韩在线| 777米奇影视久久| 日本欧美国产在线视频| 国产午夜精品一二区理论片| 午夜福利在线在线| 国产黄片视频在线免费观看| 久久久久久久久久久丰满| 亚洲av二区三区四区| 亚洲欧美成人综合另类久久久| 亚洲国产精品成人综合色| 成人亚洲精品av一区二区| 欧美激情久久久久久爽电影| 久久久久久久久大av| 精品久久久精品久久久| 欧美成人午夜免费资源| 欧美最新免费一区二区三区| 五月玫瑰六月丁香| 熟女人妻精品中文字幕| 国产精品99久久久久久久久| 搡女人真爽免费视频火全软件| 欧美zozozo另类| 女人被狂操c到高潮| 99re6热这里在线精品视频| av国产免费在线观看| 亚洲欧美日韩卡通动漫| 国产日韩欧美在线精品| 2022亚洲国产成人精品| 国产精品久久久久久av不卡| 欧美丝袜亚洲另类| 如何舔出高潮| 极品少妇高潮喷水抽搐| 精品亚洲乱码少妇综合久久| 亚洲在线观看片| 另类亚洲欧美激情| 99re6热这里在线精品视频| 精品熟女少妇av免费看| 22中文网久久字幕| 久久精品国产自在天天线| 亚洲精品乱久久久久久| 欧美另类一区| 国精品久久久久久国模美| 99久久精品国产国产毛片| 欧美丝袜亚洲另类| 97在线人人人人妻| 91狼人影院| 国产一区二区三区av在线| 91在线精品国自产拍蜜月| 国产久久久一区二区三区| 国产老妇女一区| 免费不卡的大黄色大毛片视频在线观看| 好男人在线观看高清免费视频| 国产伦精品一区二区三区四那| 91午夜精品亚洲一区二区三区| 久久ye,这里只有精品| 成年免费大片在线观看| 色网站视频免费| 国产欧美另类精品又又久久亚洲欧美| 少妇 在线观看| 人人妻人人看人人澡| 亚洲av日韩在线播放| freevideosex欧美| 国产免费一区二区三区四区乱码| 久久久精品94久久精品| 久久精品综合一区二区三区| 毛片女人毛片| 极品教师在线视频| 美女被艹到高潮喷水动态| 国产成人freesex在线| 国产日韩欧美亚洲二区| 免费av不卡在线播放| 成人毛片a级毛片在线播放| 精品人妻视频免费看| 丰满乱子伦码专区| 亚洲精品影视一区二区三区av| 大香蕉97超碰在线| 国产v大片淫在线免费观看| 国产成人a区在线观看| 亚洲精华国产精华液的使用体验| 成人亚洲精品一区在线观看 | 女的被弄到高潮叫床怎么办| 亚洲人成网站高清观看| 91狼人影院| 国产精品一区二区性色av| 精品国产三级普通话版| 女人十人毛片免费观看3o分钟| 高清毛片免费看| 制服丝袜香蕉在线| 亚洲色图综合在线观看| 成人一区二区视频在线观看| 日韩精品有码人妻一区| 国产美女午夜福利| .国产精品久久| 777米奇影视久久| 国产视频内射| 欧美xxⅹ黑人| 久久精品久久精品一区二区三区| av黄色大香蕉| 久久99蜜桃精品久久| 国产免费视频播放在线视频| 看黄色毛片网站| 又大又黄又爽视频免费| 狂野欧美白嫩少妇大欣赏| 亚洲最大成人中文| 国产成人福利小说| 亚洲图色成人| 1000部很黄的大片| 男男h啪啪无遮挡| 欧美97在线视频| 国产精品偷伦视频观看了| 好男人视频免费观看在线| 国产淫语在线视频| 亚洲怡红院男人天堂| 男人爽女人下面视频在线观看| 日韩成人伦理影院| 菩萨蛮人人尽说江南好唐韦庄| 国产淫片久久久久久久久| 欧美老熟妇乱子伦牲交| 中国三级夫妇交换| 婷婷色av中文字幕| 亚洲av日韩在线播放| 一级毛片久久久久久久久女| 国产成人精品一,二区| 神马国产精品三级电影在线观看| 亚洲国产最新在线播放| 欧美激情国产日韩精品一区| 各种免费的搞黄视频| 韩国av在线不卡| 亚洲经典国产精华液单| 舔av片在线| 国产久久久一区二区三区| 中文字幕av成人在线电影| 亚洲,一卡二卡三卡| 一级爰片在线观看| 在线观看一区二区三区| 国产精品久久久久久精品电影小说 | 国产爽快片一区二区三区| 身体一侧抽搐| 性插视频无遮挡在线免费观看| 可以在线观看毛片的网站| 久久精品国产自在天天线| av女优亚洲男人天堂| 日本熟妇午夜| 久久ye,这里只有精品| 亚洲国产精品国产精品| 国产午夜精品久久久久久一区二区三区| 亚洲综合精品二区| 亚洲欧美日韩东京热| 亚洲av福利一区| 性色av一级| 国产有黄有色有爽视频| 免费观看无遮挡的男女| 日韩欧美一区视频在线观看 | 日本与韩国留学比较| 国产免费一区二区三区四区乱码| 内射极品少妇av片p| 五月玫瑰六月丁香| 亚洲精品国产av蜜桃| 国产乱人视频| 在现免费观看毛片| 亚洲成色77777| 久久ye,这里只有精品| 国产欧美日韩一区二区三区在线 | 亚洲av免费高清在线观看| 男人爽女人下面视频在线观看| 日本三级黄在线观看| 国产白丝娇喘喷水9色精品| 蜜臀久久99精品久久宅男| 亚洲国产欧美在线一区| 欧美成人一区二区免费高清观看| 内地一区二区视频在线| av在线亚洲专区| 亚洲精品视频女| 色哟哟·www| 免费观看的影片在线观看| 亚洲图色成人| h日本视频在线播放| 久久精品久久久久久噜噜老黄| 午夜免费鲁丝| 观看免费一级毛片| 97人妻精品一区二区三区麻豆| 久久久色成人| 插逼视频在线观看| 在线观看免费高清a一片| 波多野结衣巨乳人妻| 日韩制服骚丝袜av| 久久久久久久大尺度免费视频| 欧美精品一区二区大全| 婷婷色综合www| 内射极品少妇av片p| 国产男女内射视频| 丝袜喷水一区| 亚洲国产精品专区欧美| 亚洲国产成人一精品久久久| 中文欧美无线码| 麻豆成人午夜福利视频| 免费观看在线日韩| 一区二区三区免费毛片| 人妻制服诱惑在线中文字幕| 一级毛片久久久久久久久女| 97超碰精品成人国产| 最近最新中文字幕免费大全7| 久久精品久久久久久久性| 天堂中文最新版在线下载 | 日韩欧美一区视频在线观看 | 欧美bdsm另类| 国产精品久久久久久精品古装| 麻豆久久精品国产亚洲av| av卡一久久| 日日摸夜夜添夜夜添av毛片| 国产精品三级大全| 99久久精品一区二区三区| 九九在线视频观看精品| 少妇人妻久久综合中文| 好男人视频免费观看在线| 日本三级黄在线观看| 欧美日韩综合久久久久久| 国产成人免费观看mmmm| 春色校园在线视频观看| 大又大粗又爽又黄少妇毛片口| 成人毛片60女人毛片免费| 777米奇影视久久| 欧美xxxx性猛交bbbb| 91精品一卡2卡3卡4卡| 只有这里有精品99| 超碰av人人做人人爽久久| 欧美 日韩 精品 国产| 亚洲欧美日韩东京热| 草草在线视频免费看| 插阴视频在线观看视频| 联通29元200g的流量卡| 只有这里有精品99| 丝袜美腿在线中文| 久久久久久久久久久丰满| 我的老师免费观看完整版| 国产亚洲最大av| 建设人人有责人人尽责人人享有的 | 少妇猛男粗大的猛烈进出视频 | 日本wwww免费看| 日韩欧美精品免费久久| 国产色婷婷99| 国产高清不卡午夜福利| 中国三级夫妇交换| 男插女下体视频免费在线播放| 国产精品国产av在线观看| 成人黄色视频免费在线看| 性插视频无遮挡在线免费观看| 青春草视频在线免费观看| 国产精品国产三级专区第一集| 丝袜美腿在线中文| 在线观看一区二区三区激情| 亚洲一区二区三区欧美精品 | 十八禁网站网址无遮挡 | 黄色日韩在线| 国产欧美另类精品又又久久亚洲欧美| 熟女人妻精品中文字幕| 久久鲁丝午夜福利片| 精品国产三级普通话版| 亚洲高清免费不卡视频| av在线app专区| 中文天堂在线官网| 综合色丁香网| 尤物成人国产欧美一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 又粗又硬又长又爽又黄的视频| 国产69精品久久久久777片| 亚洲av二区三区四区| 国产亚洲91精品色在线| 成年免费大片在线观看| 色网站视频免费| 亚洲精品国产成人久久av| 国产爽快片一区二区三区| 亚洲av日韩在线播放| 久久ye,这里只有精品| 啦啦啦啦在线视频资源| 看十八女毛片水多多多| 国产欧美另类精品又又久久亚洲欧美| 国产中年淑女户外野战色| 看免费成人av毛片| 777米奇影视久久| 免费看日本二区| 亚洲电影在线观看av| 精品久久久久久久久av| 日韩一区二区三区影片| 成人免费观看视频高清| 亚洲自偷自拍三级| 狂野欧美白嫩少妇大欣赏| 国产毛片a区久久久久| 高清av免费在线| 午夜精品国产一区二区电影 | 亚洲国产欧美人成| 超碰97精品在线观看| 欧美日韩综合久久久久久| 亚洲欧洲日产国产| 成年av动漫网址| 日本-黄色视频高清免费观看| av一本久久久久| 尤物成人国产欧美一区二区三区| 欧美精品人与动牲交sv欧美| 国产黄片视频在线免费观看| 青春草视频在线免费观看| 欧美日韩视频精品一区| 99久国产av精品国产电影| av在线亚洲专区| 成人国产麻豆网| 一本一本综合久久| av在线亚洲专区| 国产精品秋霞免费鲁丝片| 国产精品一区www在线观看| 不卡视频在线观看欧美| 丰满人妻一区二区三区视频av| 一级黄片播放器| 看免费成人av毛片| 国产成人福利小说| 晚上一个人看的免费电影| av.在线天堂| 啦啦啦在线观看免费高清www| 欧美丝袜亚洲另类| 国产欧美亚洲国产| 黑人高潮一二区| 特大巨黑吊av在线直播| 国产精品福利在线免费观看| 极品教师在线视频| 99久久精品热视频| 又黄又爽又刺激的免费视频.| 爱豆传媒免费全集在线观看| 日韩视频在线欧美| 一级a做视频免费观看| 80岁老熟妇乱子伦牲交| 中国国产av一级| 中国美白少妇内射xxxbb| 美女xxoo啪啪120秒动态图| 亚洲国产精品999| 成人漫画全彩无遮挡| 亚洲,欧美,日韩| 91在线精品国自产拍蜜月| 国产成人精品福利久久| 一级毛片我不卡| 亚洲电影在线观看av| 99久国产av精品国产电影| 午夜福利网站1000一区二区三区| 亚洲欧美日韩另类电影网站 | 日韩制服骚丝袜av| 在线观看一区二区三区激情| 男女边吃奶边做爰视频| 人妻一区二区av| 午夜免费鲁丝| 波野结衣二区三区在线| 日本欧美国产在线视频| 久久精品人妻少妇| 建设人人有责人人尽责人人享有的 | 最近中文字幕2019免费版| 噜噜噜噜噜久久久久久91| 丝袜喷水一区| 中文欧美无线码| 99久国产av精品国产电影| 网址你懂的国产日韩在线| 国产国拍精品亚洲av在线观看| www.av在线官网国产| xxx大片免费视频| 精品亚洲乱码少妇综合久久| 国产乱人视频| 禁无遮挡网站| 国产成人freesex在线| 啦啦啦在线观看免费高清www| 久久久久久久大尺度免费视频| 国产精品福利在线免费观看| 国产 精品1| 中文天堂在线官网| 国产精品熟女久久久久浪| 亚洲精品一二三| 草草在线视频免费看| 超碰97精品在线观看| 天天躁夜夜躁狠狠久久av| videossex国产| 内射极品少妇av片p| 男女边吃奶边做爰视频| 丝袜喷水一区| 99久久中文字幕三级久久日本| 老师上课跳d突然被开到最大视频| 国产成人91sexporn| 久久久久久久久久久丰满| 日本wwww免费看| 尤物成人国产欧美一区二区三区| 亚洲国产av新网站| 久久久久久久大尺度免费视频| 日韩av不卡免费在线播放| 日韩免费高清中文字幕av| 九草在线视频观看| 青青草视频在线视频观看| 成人欧美大片| 国产av不卡久久| 97在线人人人人妻| av免费观看日本| 麻豆乱淫一区二区| 国内少妇人妻偷人精品xxx网站| 男插女下体视频免费在线播放| 成人毛片a级毛片在线播放| 一级黄片播放器| 精华霜和精华液先用哪个| 欧美三级亚洲精品| 大片免费播放器 马上看| 国产精品.久久久| 亚洲欧美日韩卡通动漫| 欧美人与善性xxx| 午夜免费鲁丝| 日韩,欧美,国产一区二区三区| 嫩草影院精品99| 欧美高清成人免费视频www| 精品久久国产蜜桃| 国产色婷婷99| 亚洲av免费在线观看| av在线app专区| 国产一区二区在线观看日韩| 欧美3d第一页| 国产精品偷伦视频观看了| 在线免费观看不下载黄p国产| 国产免费一级a男人的天堂| 最近手机中文字幕大全| 在线a可以看的网站| 欧美成人a在线观看| 亚洲色图av天堂| 91久久精品国产一区二区成人| 国产欧美亚洲国产| 校园人妻丝袜中文字幕| 人人妻人人爽人人添夜夜欢视频 | 久久精品国产a三级三级三级| 视频区图区小说| 亚洲不卡免费看| 菩萨蛮人人尽说江南好唐韦庄| 99久久人妻综合| 免费av毛片视频| 日韩欧美一区视频在线观看 | 舔av片在线| 国产精品.久久久| 日韩亚洲欧美综合| 国产精品无大码| 中文字幕人妻熟人妻熟丝袜美| 国产黄a三级三级三级人| 国产精品99久久久久久久久| 亚洲欧美清纯卡通| 2021天堂中文幕一二区在线观| 成年av动漫网址| 高清av免费在线| 日韩国内少妇激情av| 春色校园在线视频观看| 国产精品福利在线免费观看| 中文在线观看免费www的网站| 男女国产视频网站| 网址你懂的国产日韩在线| 国产精品久久久久久精品电影小说 | 精品一区二区免费观看| 久久久精品免费免费高清| 免费观看无遮挡的男女| 交换朋友夫妻互换小说| 国产欧美日韩一区二区三区在线 | 亚洲国产日韩一区二区| 色婷婷久久久亚洲欧美| 哪个播放器可以免费观看大片| 男女那种视频在线观看| 国产综合精华液| 肉色欧美久久久久久久蜜桃 | 高清日韩中文字幕在线| 91精品一卡2卡3卡4卡| 一区二区三区免费毛片| 2022亚洲国产成人精品| 国产日韩欧美亚洲二区| 中文精品一卡2卡3卡4更新| 国产色婷婷99| 麻豆国产97在线/欧美| 久久久精品免费免费高清| 久久影院123| 亚洲在久久综合| 天美传媒精品一区二区| 特级一级黄色大片| 青青草视频在线视频观看| 亚洲在久久综合| 日日摸夜夜添夜夜爱| 精品酒店卫生间| 欧美成人一区二区免费高清观看| 国产av码专区亚洲av| 免费观看的影片在线观看| 大话2 男鬼变身卡| 日韩 亚洲 欧美在线| 日产精品乱码卡一卡2卡三| 美女高潮的动态| 精华霜和精华液先用哪个| 九九爱精品视频在线观看| 免费av毛片视频| 亚洲人成网站高清观看| 日本猛色少妇xxxxx猛交久久| 午夜免费男女啪啪视频观看| 国产精品女同一区二区软件| 免费大片黄手机在线观看| 久久午夜福利片| 国产午夜福利久久久久久| 黄色视频在线播放观看不卡| 国产免费福利视频在线观看| 一级av片app| 久久人人爽人人爽人人片va| 欧美激情久久久久久爽电影| 亚洲欧美一区二区三区国产| 精品一区二区免费观看| 亚洲国产精品999| 99久久精品国产国产毛片| 免费观看a级毛片全部| 国产色婷婷99| 久久韩国三级中文字幕| freevideosex欧美| 真实男女啪啪啪动态图| 久久99热这里只频精品6学生| 两个人的视频大全免费| 欧美高清性xxxxhd video| 国产视频首页在线观看| 国产老妇女一区| 干丝袜人妻中文字幕| 国产欧美另类精品又又久久亚洲欧美| 亚洲最大成人中文| 国产成人aa在线观看| 免费av毛片视频| 久久久久久久久久成人| 亚洲国产精品成人久久小说| 一区二区av电影网| 日日摸夜夜添夜夜添av毛片| 免费观看的影片在线观看| 日韩中字成人| 人妻 亚洲 视频| 一区二区三区乱码不卡18| 亚洲熟女精品中文字幕| 秋霞伦理黄片| 亚洲精品影视一区二区三区av| 久久久午夜欧美精品| 国产亚洲av片在线观看秒播厂| 精品少妇久久久久久888优播| av网站免费在线观看视频| 日日啪夜夜爽| 天美传媒精品一区二区| 制服丝袜香蕉在线| 麻豆精品久久久久久蜜桃| 又大又黄又爽视频免费| 日韩欧美一区视频在线观看 | 亚洲av成人精品一二三区| 成人欧美大片| 18禁裸乳无遮挡免费网站照片| 免费观看a级毛片全部| 午夜视频国产福利| 美女cb高潮喷水在线观看| 少妇丰满av| 尾随美女入室| 简卡轻食公司| 免费电影在线观看免费观看| 看免费成人av毛片| 99热6这里只有精品| 精品午夜福利在线看| 久久女婷五月综合色啪小说 | 街头女战士在线观看网站| 国产淫语在线视频| 日韩欧美精品免费久久|