• <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
    国产精品乱码一区二三区的特点| 亚洲国产精品专区欧美| 国产欧美日韩精品一区二区| 蜜桃亚洲精品一区二区三区| 免费观看的影片在线观看| 久久久久久国产a免费观看| 26uuu在线亚洲综合色| 白带黄色成豆腐渣| 午夜精品一区二区三区免费看| videos熟女内射| 久久久久久久久久久丰满| 午夜精品一区二区三区免费看| 女人久久www免费人成看片 | 三级国产精品片| 别揉我奶头 嗯啊视频| 中文资源天堂在线| 国产亚洲5aaaaa淫片| 在线免费观看的www视频| 在线观看美女被高潮喷水网站| 国产精品蜜桃在线观看| 国产伦理片在线播放av一区| 亚洲精品自拍成人| 乱人视频在线观看| 国产亚洲一区二区精品| 国产伦理片在线播放av一区| 欧美日韩在线观看h| 国产一区二区在线av高清观看| 亚洲综合精品二区| 久久久a久久爽久久v久久| 欧美激情久久久久久爽电影| 免费播放大片免费观看视频在线观看 | 精品久久久久久久久亚洲| 丰满乱子伦码专区| 午夜激情福利司机影院| 欧美日本视频| 亚洲真实伦在线观看| 一边亲一边摸免费视频| 熟女人妻精品中文字幕| 亚洲av中文字字幕乱码综合| 欧美又色又爽又黄视频| 欧美性猛交╳xxx乱大交人| 国产淫片久久久久久久久| 最近最新中文字幕大全电影3| 成年版毛片免费区| 97人妻精品一区二区三区麻豆| 色综合站精品国产| 非洲黑人性xxxx精品又粗又长| 天堂网av新在线| 99视频精品全部免费 在线| 男女下面进入的视频免费午夜| 国产精品精品国产色婷婷| 成人av在线播放网站| 亚洲av成人精品一区久久| 在线a可以看的网站| 26uuu在线亚洲综合色| 视频中文字幕在线观看| 亚洲av电影不卡..在线观看| 国语自产精品视频在线第100页| 中文欧美无线码| 乱人视频在线观看| 女的被弄到高潮叫床怎么办| 国产白丝娇喘喷水9色精品| 日韩欧美精品v在线| 久久久精品欧美日韩精品| 精品人妻熟女av久视频| kizo精华| 亚洲精品日韩av片在线观看| 亚洲国产精品成人综合色| 欧美成人精品欧美一级黄| 国产午夜精品论理片| 亚洲激情五月婷婷啪啪| 国产精品日韩av在线免费观看| 1024手机看黄色片| a级一级毛片免费在线观看| 欧美精品一区二区大全| 久久久精品欧美日韩精品| 久久精品夜色国产| 亚洲国产日韩欧美精品在线观看| 91aial.com中文字幕在线观看| or卡值多少钱| 国产一级毛片在线| 99久久九九国产精品国产免费| 国产 一区精品| 少妇裸体淫交视频免费看高清| 久久精品久久久久久久性| 亚洲精华国产精华液的使用体验| 日本与韩国留学比较| 久久综合国产亚洲精品| 成年av动漫网址| 国产白丝娇喘喷水9色精品| 国内揄拍国产精品人妻在线| 午夜老司机福利剧场| 亚洲精品一区蜜桃| 免费观看在线日韩| 色哟哟·www| 国产人妻一区二区三区在| 大又大粗又爽又黄少妇毛片口| 国产成人精品久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 久久午夜福利片| 99久久成人亚洲精品观看| 亚洲成色77777| 26uuu在线亚洲综合色| 亚洲国产精品成人综合色| 亚洲不卡免费看| 一级黄色大片毛片| 高清av免费在线| 简卡轻食公司| 亚洲国产日韩欧美精品在线观看| 午夜a级毛片| 日韩av在线大香蕉| 亚洲av成人精品一区久久| 永久免费av网站大全| 国产成人精品婷婷| 最近2019中文字幕mv第一页| 伦精品一区二区三区| 亚洲自拍偷在线| 午夜亚洲福利在线播放| 18禁动态无遮挡网站| 搞女人的毛片| 男女啪啪激烈高潮av片| 亚洲av福利一区| 成人一区二区视频在线观看| 日本三级黄在线观看| 美女高潮的动态| 成人漫画全彩无遮挡| 黄片wwwwww| 22中文网久久字幕| 午夜福利高清视频| 国产精品国产高清国产av| 中文字幕av在线有码专区| 欧美成人免费av一区二区三区| 六月丁香七月| 国产成人免费观看mmmm| 一夜夜www| 亚洲成人中文字幕在线播放| 中国美白少妇内射xxxbb| videos熟女内射| av天堂中文字幕网| 亚洲精华国产精华液的使用体验| 亚洲在久久综合| 成人亚洲欧美一区二区av| 黄片wwwwww| 亚洲国产欧美人成| 卡戴珊不雅视频在线播放| 嫩草影院入口| videossex国产| 国产免费一级a男人的天堂| av国产久精品久网站免费入址| 久久精品久久精品一区二区三区| 国产成人一区二区在线| 亚洲欧美成人精品一区二区| 亚洲成人精品中文字幕电影| 亚洲av成人精品一二三区| 菩萨蛮人人尽说江南好唐韦庄 | 91狼人影院| 99在线人妻在线中文字幕| 日韩av在线大香蕉| 最近2019中文字幕mv第一页| 亚洲av不卡在线观看| 国产在线一区二区三区精 | 久久精品国产99精品国产亚洲性色| 人妻系列 视频| 五月玫瑰六月丁香| 亚洲熟妇中文字幕五十中出| 日本免费a在线| 亚洲国产精品成人久久小说| 久久久久性生活片| 亚洲美女搞黄在线观看| 国产精品人妻久久久久久| 国产精品乱码一区二三区的特点| 婷婷色麻豆天堂久久 | 亚洲最大成人av| 91久久精品国产一区二区成人| 久久精品久久久久久久性| 99久久成人亚洲精品观看| 国产单亲对白刺激| 观看免费一级毛片| 午夜视频国产福利| 色吧在线观看| 国产精品无大码| 丰满人妻一区二区三区视频av| 久久99蜜桃精品久久| 国产麻豆成人av免费视频| 日韩亚洲欧美综合| 国语对白做爰xxxⅹ性视频网站| 中文字幕av成人在线电影| 99视频精品全部免费 在线| 亚洲精品日韩在线中文字幕| 久99久视频精品免费| 欧美97在线视频| 小说图片视频综合网站| 亚洲精品成人久久久久久| 91久久精品国产一区二区三区| 国模一区二区三区四区视频| 干丝袜人妻中文字幕| 国产又色又爽无遮挡免| 精品国产露脸久久av麻豆 | 国产免费又黄又爽又色| 精品少妇黑人巨大在线播放 | 国产高清国产精品国产三级 | 少妇熟女欧美另类| 亚洲av不卡在线观看| 啦啦啦观看免费观看视频高清| 黄色配什么色好看| 日韩一本色道免费dvd| 免费看a级黄色片| 国产欧美另类精品又又久久亚洲欧美| 一级二级三级毛片免费看| 一本久久精品| 欧美高清性xxxxhd video| 97超碰精品成人国产| 亚洲精品日韩在线中文字幕| 美女大奶头视频| 亚洲无线观看免费| 免费无遮挡裸体视频| 久久国产乱子免费精品| 免费看美女性在线毛片视频| 国产麻豆成人av免费视频| 99热精品在线国产| 中文天堂在线官网| 日本猛色少妇xxxxx猛交久久| 嫩草影院新地址| 国内精品宾馆在线| av线在线观看网站| 天堂av国产一区二区熟女人妻| 天堂影院成人在线观看| kizo精华| 日本免费一区二区三区高清不卡| 久久久久久久久大av| 久久综合国产亚洲精品| 最新中文字幕久久久久| 精品久久久久久久久av| 日韩亚洲欧美综合| 精品久久久久久成人av| 日本一二三区视频观看| 亚洲五月天丁香| 成人鲁丝片一二三区免费| 免费黄网站久久成人精品| 两性午夜刺激爽爽歪歪视频在线观看| 小说图片视频综合网站| 26uuu在线亚洲综合色| 日日干狠狠操夜夜爽| 久久99热6这里只有精品| 亚洲av电影不卡..在线观看| av在线亚洲专区| 精品国产露脸久久av麻豆 | 黄色欧美视频在线观看| 黄色配什么色好看| 精品人妻视频免费看| 欧美性猛交黑人性爽| 国产精品美女特级片免费视频播放器| 99久久精品一区二区三区| 91aial.com中文字幕在线观看| 亚洲av免费高清在线观看| 久久久久国产网址| 嫩草影院精品99| 国产综合懂色| 国产伦精品一区二区三区视频9| 国产黄色视频一区二区在线观看 | a级毛片免费高清观看在线播放| 麻豆成人午夜福利视频| 国产成人午夜福利电影在线观看| 3wmmmm亚洲av在线观看| 成人午夜高清在线视频| av黄色大香蕉| 一区二区三区高清视频在线| 最近视频中文字幕2019在线8| 青春草国产在线视频| 欧美日本视频| 最新中文字幕久久久久| 大香蕉97超碰在线| 日韩欧美在线乱码| 五月伊人婷婷丁香| 毛片女人毛片| 中国美白少妇内射xxxbb| 日本与韩国留学比较| 69人妻影院| 中文字幕免费在线视频6| av天堂中文字幕网| 久久久久久久午夜电影| 精品人妻熟女av久视频| 国产在视频线精品| 精品人妻视频免费看| 国内揄拍国产精品人妻在线| 欧美性猛交黑人性爽| 建设人人有责人人尽责人人享有的 | 波多野结衣高清无吗| 久久久成人免费电影| av在线播放精品| 午夜老司机福利剧场| 亚洲图色成人| 白带黄色成豆腐渣| 精品国产三级普通话版| 亚洲高清免费不卡视频| 免费黄色在线免费观看| 美女被艹到高潮喷水动态| av福利片在线观看| 亚洲怡红院男人天堂| 精品久久久久久久末码| 中文字幕免费在线视频6| 一级黄色大片毛片| 久久99热这里只有精品18| 一本久久精品| 午夜久久久久精精品| 欧美三级亚洲精品| 国产男人的电影天堂91| 中文字幕精品亚洲无线码一区| 日韩在线高清观看一区二区三区| 看非洲黑人一级黄片| 春色校园在线视频观看| 一级黄色大片毛片| 国模一区二区三区四区视频| 小说图片视频综合网站| 能在线免费看毛片的网站| 国产精品国产高清国产av| 在线观看美女被高潮喷水网站| 人人妻人人看人人澡| 亚洲18禁久久av| 在线播放国产精品三级| 美女国产视频在线观看| av线在线观看网站| 少妇丰满av| 91精品国产九色| 亚洲欧美精品专区久久| 亚洲国产色片| 日韩欧美 国产精品| 一边摸一边抽搐一进一小说| 人体艺术视频欧美日本| 午夜免费激情av| 国产黄片美女视频| 精品久久久久久久末码| 国产老妇女一区| 国产精品久久久久久av不卡| .国产精品久久| 国产免费视频播放在线视频 | 深夜a级毛片| 麻豆精品久久久久久蜜桃| 亚洲激情五月婷婷啪啪| 午夜激情欧美在线| 久久久亚洲精品成人影院| 大香蕉97超碰在线| 麻豆久久精品国产亚洲av| 18禁裸乳无遮挡免费网站照片| 白带黄色成豆腐渣| 一区二区三区免费毛片| 高清午夜精品一区二区三区| 亚洲四区av| 国产熟女欧美一区二区| 99久国产av精品国产电影| 在线a可以看的网站| 六月丁香七月| 插逼视频在线观看| 男人的好看免费观看在线视频| 三级男女做爰猛烈吃奶摸视频| 综合色av麻豆| 亚洲国产精品成人久久小说| 青春草亚洲视频在线观看| videos熟女内射| 日韩精品有码人妻一区| 青春草国产在线视频| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品成人综合色| 97超碰精品成人国产| 日韩三级伦理在线观看| 欧美成人午夜免费资源| 欧美日本亚洲视频在线播放| 亚洲国产色片| 国产午夜福利久久久久久| 欧美xxxx黑人xx丫x性爽| 伊人久久精品亚洲午夜| 中文字幕熟女人妻在线| 免费观看精品视频网站| 日韩大片免费观看网站 | 精品酒店卫生间| 亚洲av一区综合| 欧美变态另类bdsm刘玥| 免费电影在线观看免费观看| 亚洲国产色片| 青春草视频在线免费观看| 欧美日本亚洲视频在线播放| 在线观看美女被高潮喷水网站| 国产淫片久久久久久久久| 中文字幕久久专区| 波多野结衣巨乳人妻| 黄色日韩在线| 日韩大片免费观看网站 | 麻豆av噜噜一区二区三区| 亚洲第一区二区三区不卡| 一个人免费在线观看电影| 高清日韩中文字幕在线| 午夜福利视频1000在线观看| 天堂影院成人在线观看| 黄片无遮挡物在线观看| av女优亚洲男人天堂| 欧美不卡视频在线免费观看| 欧美一区二区国产精品久久精品| 精品久久久久久久久久久久久| 日产精品乱码卡一卡2卡三| 国产精品久久久久久av不卡| 亚洲欧美成人精品一区二区| 国产精品人妻久久久影院| 国产精品永久免费网站| 村上凉子中文字幕在线| 最近手机中文字幕大全| 伊人久久精品亚洲午夜| 亚洲伊人久久精品综合 | 老女人水多毛片| 热99在线观看视频| 亚洲一区高清亚洲精品| 日韩av不卡免费在线播放| 99在线视频只有这里精品首页| 精品不卡国产一区二区三区| 精品久久国产蜜桃| 日本黄色片子视频| 国内精品美女久久久久久| videossex国产| 国产成人一区二区在线| 热99在线观看视频| 秋霞在线观看毛片| 久久精品91蜜桃| 国产精品女同一区二区软件| 亚洲在线自拍视频| 日本午夜av视频| 精品欧美国产一区二区三| 少妇高潮的动态图| 狂野欧美激情性xxxx在线观看| av黄色大香蕉| av免费观看日本| 2021少妇久久久久久久久久久| 亚洲最大成人中文| 久久精品综合一区二区三区| 久久人妻av系列| eeuss影院久久| 秋霞在线观看毛片| 亚洲经典国产精华液单| 听说在线观看完整版免费高清| 成人美女网站在线观看视频| 久久精品久久久久久噜噜老黄 | 午夜福利在线在线| 在线免费十八禁| 国内精品一区二区在线观看| 我要看日韩黄色一级片| 日日干狠狠操夜夜爽| 国产精品久久电影中文字幕| 亚洲色图av天堂| 精品久久国产蜜桃| 国产在线一区二区三区精 | 久久99热这里只频精品6学生 | 亚洲av免费在线观看| av黄色大香蕉| 丰满人妻一区二区三区视频av| 国产成人精品婷婷| 免费看av在线观看网站| 日本免费一区二区三区高清不卡| 日本欧美国产在线视频| .国产精品久久| 国产成人a∨麻豆精品| 国产精品一区二区三区四区免费观看| 日韩强制内射视频| 国产在线男女| 精品免费久久久久久久清纯| 全区人妻精品视频| 最近2019中文字幕mv第一页| 黄片wwwwww| 国产 一区 欧美 日韩| 成年女人永久免费观看视频| 国产精品99久久久久久久久| 一本久久精品| 亚洲精品日韩在线中文字幕| 亚洲人成网站在线播| 日韩大片免费观看网站 | 99久久中文字幕三级久久日本| 国产高清三级在线| 噜噜噜噜噜久久久久久91| 伊人久久精品亚洲午夜| 美女cb高潮喷水在线观看| 国产爱豆传媒在线观看| 晚上一个人看的免费电影| 免费电影在线观看免费观看| 午夜免费男女啪啪视频观看| 一个人观看的视频www高清免费观看| 永久网站在线| 成人鲁丝片一二三区免费| 色综合亚洲欧美另类图片| 青春草视频在线免费观看| 99九九线精品视频在线观看视频| 夫妻性生交免费视频一级片| 国产成人aa在线观看| 亚洲成人中文字幕在线播放| 一级爰片在线观看| 最近视频中文字幕2019在线8| 网址你懂的国产日韩在线| 欧美色视频一区免费| 村上凉子中文字幕在线| 国产综合懂色| 中国国产av一级| 亚洲经典国产精华液单| 精品人妻熟女av久视频| 26uuu在线亚洲综合色| 国产日韩欧美在线精品| 免费大片18禁| 欧美区成人在线视频| 少妇人妻精品综合一区二区| 久久精品国产自在天天线| 精品久久久噜噜| 成人漫画全彩无遮挡| 网址你懂的国产日韩在线| 简卡轻食公司| 男女下面进入的视频免费午夜| 日韩视频在线欧美| 国产黄片视频在线免费观看| 国产免费福利视频在线观看| 日本黄色片子视频| 国产精品国产三级国产av玫瑰| 日韩在线高清观看一区二区三区| videossex国产| 日韩av不卡免费在线播放| 久久久久久久国产电影| 麻豆国产97在线/欧美| 日韩,欧美,国产一区二区三区 | av在线观看视频网站免费| 国产在线一区二区三区精 | 一级爰片在线观看| 观看美女的网站| 天堂中文最新版在线下载 | 男女国产视频网站| .国产精品久久| 三级毛片av免费| 日本wwww免费看| 亚洲av中文字字幕乱码综合| 嫩草影院入口| 日日啪夜夜撸| videos熟女内射| 成人漫画全彩无遮挡| 国产av码专区亚洲av| 免费在线观看成人毛片| 国产精品一及| 91狼人影院| 毛片女人毛片| 精品免费久久久久久久清纯| 麻豆久久精品国产亚洲av| 国产v大片淫在线免费观看| 亚洲av成人av| 人妻制服诱惑在线中文字幕| 国产真实乱freesex| 欧美bdsm另类| 欧美潮喷喷水| 天天躁日日操中文字幕| 亚洲精品乱码久久久久久按摩| 久久精品国产亚洲av天美| a级毛片免费高清观看在线播放| 国产精品一区二区在线观看99 | 久久精品国产99精品国产亚洲性色| 亚洲av免费在线观看| 国产成人freesex在线| 一本一本综合久久| 欧美日韩精品成人综合77777| 99热精品在线国产| 色播亚洲综合网| 国产三级中文精品| 国产大屁股一区二区在线视频| 欧美一区二区亚洲| 小说图片视频综合网站| 日本黄色片子视频| 国产男人的电影天堂91| 精品酒店卫生间| 午夜激情欧美在线| 丰满人妻一区二区三区视频av| 亚洲在久久综合| 在线免费十八禁| 日本wwww免费看| 亚洲国产日韩欧美精品在线观看| 久久国内精品自在自线图片| 精品久久久久久久久亚洲| 天堂av国产一区二区熟女人妻| 欧美xxxx黑人xx丫x性爽| 亚洲真实伦在线观看| 亚洲av成人精品一区久久| 亚洲av熟女| 日韩av在线免费看完整版不卡| 中文字幕人妻熟人妻熟丝袜美| 国产精品精品国产色婷婷| 午夜福利在线观看免费完整高清在| 特大巨黑吊av在线直播| 色综合站精品国产| 在线天堂最新版资源| 婷婷色av中文字幕| 日本猛色少妇xxxxx猛交久久| 夫妻性生交免费视频一级片| 亚洲一区高清亚洲精品| 变态另类丝袜制服| 波野结衣二区三区在线| 啦啦啦啦在线视频资源| 高清av免费在线| 夫妻性生交免费视频一级片| 欧美丝袜亚洲另类| 97超视频在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 乱系列少妇在线播放| 久久99热这里只有精品18| 国产伦一二天堂av在线观看| 国产真实乱freesex| 国产精品久久久久久精品电影| 高清视频免费观看一区二区 | 啦啦啦韩国在线观看视频| 卡戴珊不雅视频在线播放| 久久99精品国语久久久| 大又大粗又爽又黄少妇毛片口| 亚洲婷婷狠狠爱综合网| 91aial.com中文字幕在线观看| av专区在线播放| 成人无遮挡网站| 哪个播放器可以免费观看大片| 小蜜桃在线观看免费完整版高清|