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

    基于CEEMD-BP模型的汾河上游月徑流預(yù)測

    2019-08-24 08:53:18桑宇婷趙雪花祝雪萍席東潔
    人民黃河 2019年8期
    關(guān)鍵詞:汾河水文站徑流

    桑宇婷,趙雪花,祝雪萍,席東潔

    (太原理工大學(xué)水利科學(xué)與工程學(xué)院,山西太原030024)

    隨著經(jīng)濟(jì)社會(huì)的發(fā)展,溫室氣體大量排放,使得全球平均氣溫升高,降水分布不均勻度加劇,極端天氣愈發(fā)頻繁,洪水、干旱等自然災(zāi)害增多[1-2]。準(zhǔn)確、合理地對(duì)徑流進(jìn)行分析預(yù)測,對(duì)水資源的合理開發(fā)利用及生態(tài)環(huán)境保護(hù)具有重要意義。受氣候、地理環(huán)境、社會(huì)發(fā)展及人類活動(dòng)的影響,徑流的變化規(guī)律錯(cuò)綜復(fù)雜,具體表現(xiàn)為突變性、非線性、隨機(jī)性等特性[3]。為克服徑流序列非平穩(wěn)性的影響,Huang等于1998年提出了經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,簡稱EMD)法,該方法是一種自適應(yīng)的局部化分析方法,可將非平穩(wěn)的徑流序列分解為有限個(gè)具有一定周期及平穩(wěn)性的固有模態(tài)函數(shù)(IMF),但EMD會(huì)出現(xiàn)模態(tài)混疊現(xiàn)象[4-5],影響后續(xù)徑流預(yù)測的精度。為解決EMD的模態(tài)混疊問題,Huang等于2009年提出了集合經(jīng)驗(yàn)?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition,簡稱EEMD)法,將具有均勻分布特征的白噪聲加入到待分解的信號(hào)中[6],但其分解效率較低。因此,Yeh等于2010年提出了互補(bǔ)集合經(jīng)驗(yàn)?zāi)B(tài)分解(Complete Ensemble Empirical Mode Decomposition,簡稱 CEEMD)法[7],主要是通過向待分解信號(hào)中加入兩個(gè)互補(bǔ)的白噪聲信號(hào)并進(jìn)行分解,此方法不僅能得到與EEMD法相同的分解效果,而且減小了由白噪聲引起的重構(gòu)誤差,提高了分解效率,使徑流序列得到了有效的平穩(wěn)化處理。

    在目前已有的多種徑流預(yù)測方法中,人工神經(jīng)網(wǎng)絡(luò)是一種模擬動(dòng)物神經(jīng)網(wǎng)絡(luò)進(jìn)行信息處理的算法數(shù)學(xué)模型,具有自適應(yīng)、自學(xué)習(xí)及容錯(cuò)性高等特點(diǎn),可有效克服徑流序列非線性的影響。BP神經(jīng)網(wǎng)絡(luò)是應(yīng)用較為普遍的人工神經(jīng)網(wǎng)絡(luò)之一[8],其結(jié)構(gòu)簡單,具有較強(qiáng)的泛化能力和良好的非線性映射能力[9],能夠在復(fù)雜的徑流序列中找出一定的變化規(guī)律。

    隨著預(yù)測模型的不斷改進(jìn)和預(yù)測機(jī)制認(rèn)識(shí)的不斷深化,為了降低預(yù)測誤差,人們更趨向于利用多種模型進(jìn)行組合預(yù)測。本文通過構(gòu)建CEEMD-BP模型,對(duì)汾河上游月徑流進(jìn)行了預(yù)測,并與單一BP模型及EMDBP模型的預(yù)測結(jié)果進(jìn)行對(duì)比,以驗(yàn)證該模型預(yù)測非平穩(wěn)、非線性徑流序列的可行性與有效性。

    1 方法理論簡介

    1.1 CEEMD法的基本理論

    CEEMD法通過對(duì)月徑流序列進(jìn)行除噪分解,得到多個(gè)具有一定物理意義的IMF分量和趨勢項(xiàng)(Res),進(jìn)而確定月徑流序列的周期性及變化趨勢。CEEMD法可依據(jù)各層分解的頻譜提取有效信號(hào),消除高頻隨機(jī)噪聲,使徑流序列平穩(wěn)化。其主要分解步驟如下[10-11]。

    (1)設(shè) x(t)為原序列,則 IMF1的計(jì)算公式為

    (2)求一階殘差 r1(t),計(jì)算公式為

    (3)計(jì)算 IMF2,公式為

    (4)求 k 階殘差 rk(t),計(jì)算公式為

    (5)計(jì)算原序列的 IMFk+1(t),公式為

    (6)重復(fù)步驟(4)~(5),直至殘差不能再分解,求得最終殘差 R(t):

    則原序列x(t)的表達(dá)式為

    式中:I為集總次數(shù),即每次求取分量時(shí)進(jìn)行I次分解;ωi(t)為單位方差的零均值高斯白噪聲;εk為控制所添加噪聲和原始信號(hào)信噪比的噪聲系數(shù);Ek(·)為定義好的算子;K為IMF(t)的總數(shù)。

    1.2 BP 神經(jīng)網(wǎng)絡(luò)

    1986年,Rumelhart D E提出了BP神經(jīng)網(wǎng)絡(luò),它是一種按誤差逆?zhèn)鞑ニ惴ㄓ?xùn)練的多層前饋網(wǎng)絡(luò),普遍應(yīng)用于徑流預(yù)測。BP神經(jīng)網(wǎng)絡(luò)包括輸入層、隱含層和輸出層,網(wǎng)絡(luò)結(jié)構(gòu)見圖1,通過各層之間的信息正向傳播和誤差反向傳播來確定權(quán)值與閾值。具體訓(xùn)練過程見文獻(xiàn)[12-15]。

    1.3 CEEMD-BP 模型

    CEEMD-BP模型進(jìn)行徑流預(yù)測的基本原理:月徑流序列通過CEEMD法分解,提取月徑流序列中具有一定物理意義的信號(hào),并降低噪音干擾的影響,得到不同頻率的子序列,再通過BP神經(jīng)網(wǎng)絡(luò)對(duì)各子序列進(jìn)行預(yù)測,最后對(duì)各子序列預(yù)測結(jié)果進(jìn)行重構(gòu),得到月徑流預(yù)測值。模型預(yù)測過程見圖2。

    2 實(shí)例分析

    2.1 汾河上游月徑流資料分析

    汾河位于黃河中游,全長716 km,東西寬188 km,南北長412 km,從寧武縣管涔山雷鳴寺到太原市蘭村為汾河上游,全長216.9 km,流域面積為7 705 km2,屬溫帶大陸性氣候區(qū)。汾河上游水文站中,上靜游站位于汾河支流嵐河上,此站上游集水面積為1 140 km2;汾河水庫站位于嵐河與汾河干流交匯處,集水面積約32 km2,水庫總庫容為7億m3;寨上站坐落于古交市寨上村,位于汾河干流,集水面積約6 800 km2,與汾河水庫站之間的流域面積約1 500 km2;蘭村水文站位于太原市西北22.5 km的上蘭村,地處汾河中上游,集水面積約7 700 km2。4個(gè)水文站1958—2000年月徑流數(shù)據(jù)資料來自山西省水文局,已被運(yùn)用于多個(gè)科研項(xiàng)目及工程,因此可以保證可靠性,其變差系數(shù)、偏度、峰度見表1。

    表1 各站月徑流特征值

    由表1可知:①上靜游站位于汾河支流,其月徑流量均值在4個(gè)水文站中最?。虎诜诤由嫌胃髡镜淖儾钕禂?shù)為1.35~1.61,其月徑流量變化劇烈,徑流年內(nèi)分配不均,上靜游站變差系數(shù)最大,這主要與汾河水庫對(duì)下游的調(diào)節(jié)作用有關(guān);③4個(gè)水文站的偏度均大于0,分布為右偏態(tài);④4個(gè)水文站的峰度均大于40,遠(yuǎn)遠(yuǎn)大于正態(tài)分布及均勻分布的峰度,說明月徑流量數(shù)值與均值相差較大,分布有更多的極端值,具有尖峰厚尾的特征。綜上所述,汾河上游4個(gè)水文站的月徑流序列是具有較多極端值、存在右偏現(xiàn)象且變化劇烈的非正態(tài)序列,即為非平穩(wěn)、非線性序列。

    2.2 月徑流資料的CEEMD法分解

    利用CEEMD法對(duì)4個(gè)水文站的月徑流資料進(jìn)行分解,由分解結(jié)果可知4個(gè)水文站從IMF1到趨勢項(xiàng)的子序列均呈現(xiàn)頻率降低、波長變長、振幅變小的特征。各水文站的分解結(jié)果中,由于IMF1~I(xiàn)MF3中包含了原徑流序列中的主要成分,因此其仍表現(xiàn)出較高的頻率、較短的波長及較大的振幅;上靜游站、汾河水庫站、寨上站的IMF4~I(xiàn)MF7及蘭村站的IMF4~I(xiàn)MF6已明顯表現(xiàn)出一定的周期性及比較規(guī)律的波動(dòng)性,其非線性、非平穩(wěn)性均有所下降;各站的Res已基本表現(xiàn)出徑流序列的長期變化趨勢,即上靜游站及汾河水庫站為下降趨勢、寨上站為先下降后上升趨勢、蘭村站總體為先下降后小幅上升再下降趨勢。由此可見,經(jīng)過CEEMD法分解可得到比原徑流序列更加平穩(wěn)的子序列。

    為進(jìn)一步說明CEEMD法分解效果,將CEEMD法分解結(jié)果與EMD分解結(jié)果進(jìn)行了對(duì)比,可知EMD法分解結(jié)果的變化趨勢與CEEMD法的相似。但是,上靜游站、汾河水庫站、寨上站、蘭村站月徑流CEEMD法分解得到的IMF分別為8、8、8、7個(gè),而EMD法得到的IMF分別為8、9、9、8個(gè)。除上靜游站外,其他3個(gè)水文站月徑流CEEMD法分解出的IMF比EMD法的少,且EMD法分解的各分量頻率曲線中頻率變化比CEEMD法的多,即CEEMD法分解減少了EMD法各分量中的某些頻率信號(hào),可以推測CEEMD法能解決EMD法的模態(tài)混疊問題。通過CEEMD法分解的各分量比前一個(gè)分量頻率減小程度比EMD法的更加明顯,且CEEMD法分解的趨勢項(xiàng)也更加趨于線性平穩(wěn),由此驗(yàn)證了CEEMD法在保留了EMD法處理非平穩(wěn)問題的優(yōu)勢的基礎(chǔ)上,能夠有效解決EMD法的模態(tài)混疊問題,這為后續(xù)BP神經(jīng)網(wǎng)絡(luò)預(yù)測精度的提高奠定了基礎(chǔ)。

    2.3 預(yù)測結(jié)果分析

    將上靜游站、汾河水庫站、寨上站及蘭村站1958—2000年516個(gè)月的徑流資料分為兩部分,1958年1月—1999年6月498個(gè)月的徑流數(shù)據(jù)用于模型建模,1999年7月—2000年12月18個(gè)月的徑流數(shù)據(jù)用于模型驗(yàn)證。

    本文使用的BP神經(jīng)網(wǎng)絡(luò)輸入層節(jié)點(diǎn)數(shù)為31,輸出層節(jié)點(diǎn)數(shù)為18,根據(jù)前文分解出的IMF的不同頻率選擇合適的隱含層節(jié)點(diǎn)數(shù),高頻項(xiàng)、低頻項(xiàng)、趨勢項(xiàng)對(duì)應(yīng)的隱含層節(jié)點(diǎn)數(shù)分別為18、15、17。確定好BP神經(jīng)網(wǎng)絡(luò)參數(shù)后,各IMF同樣被分為兩部分,前498個(gè)數(shù)據(jù)進(jìn)行模擬,對(duì)后18個(gè)數(shù)據(jù)進(jìn)行預(yù)測,再對(duì)各IMF的后18個(gè)預(yù)測數(shù)據(jù)進(jìn)行重構(gòu),與原徑流數(shù)據(jù)進(jìn)行對(duì)比,以驗(yàn)證模型的預(yù)測效果,結(jié)果見圖3。

    由圖3可知:①使用單一BP模型預(yù)測時(shí),上靜游站2000年2—7月的預(yù)測值明顯大于實(shí)測值,其中2月偏差達(dá)到50%,BP模型預(yù)測結(jié)果是3個(gè)模型中偏差最大的;汾河水庫站除了2000年8月外,其余月份預(yù)測值均與實(shí)測值有較大偏差;寨上站出現(xiàn)兩處偏差大于50%的值,分別為1999年7月、2000年4月;蘭村站在2000年6月同樣出現(xiàn)偏差大于50%的值,但其預(yù)測值與實(shí)測值的變化趨勢基本相同,因此雖然BP模型預(yù)測精度較另外兩種模型低,但其也適用于非線性徑流序列的預(yù)測。②EMD-BP模型的預(yù)測效果整體較單一BP模型更好,均更加接近實(shí)測值,變化趨勢也基本相同,說明徑流經(jīng)過平穩(wěn)化處理后再進(jìn)行預(yù)測,精度會(huì)有所提高,但汾河水庫站2000年9—12月通過EMD-BP模型預(yù)測出的結(jié)果是3種模型中偏差最大的,由此推測EMD法存在一定的模態(tài)混疊問題,對(duì)徑流預(yù)測有一定影響;改進(jìn)EMD法分解過程的CEEMD-BP模型模擬各時(shí)段徑流的擬合程度均比EMD-BP模型更高。③CEEMD-BP模型的預(yù)測結(jié)果與實(shí)測徑流曲線擬合度最好,變化趨勢均與原序列相同,且無極端值出現(xiàn)。

    綜上,CEEMD法在保留平穩(wěn)化處理的優(yōu)勢下,可有效克服EMD模態(tài)混疊問題對(duì)預(yù)測效果的影響,且BP模型適用于非線性徑流序列的預(yù)測,這說明利用CEEMD-BP模型對(duì)非平穩(wěn)、非線性徑流序列進(jìn)行預(yù)測是可行、有效的。

    為了更加清楚地反映模型預(yù)測結(jié)果的誤差及預(yù)測精度,采用平均絕對(duì)誤差(MAE)、均方根誤差(RMSE)及確定性系數(shù)(NS)這3個(gè)指標(biāo)進(jìn)行分析,3個(gè)指標(biāo)的計(jì)算公式分別為

    式中:σi為i時(shí)刻的預(yù)測值;σ0為i時(shí)刻的實(shí)測值;σ為實(shí)測值的均值。

    根據(jù)《水文情報(bào)預(yù)報(bào)規(guī)范》(SL 250—2000)的規(guī)定,當(dāng) NS≥0.9 時(shí),預(yù)測精度為甲級(jí);0.7≤NS<0.9 時(shí),預(yù)測精度為乙級(jí);0.5≤NS<0.7時(shí),預(yù)測精度為丙級(jí);NS<0.5時(shí),預(yù)測結(jié)果不可信。汾河上游4個(gè)水文站的誤差結(jié)果見表2。

    由表2可知:①模擬期CEEMD-BP模型的誤差與BP模型、EMD-BP模型相比大幅減小,CEEMD-BP模型的MAE、RMSE與單一BP模型相比分別減小75%~82%、57%~70%;與 EMD-BP模型相比分別減小62%~77%、40%~64%。 BP 模型的 NS為 0.83~0.87,預(yù)測精度為乙級(jí);EMD-BP 模型的 NS 為 0.87~0.91,其中上靜游站、汾河水庫站的模擬精度為乙級(jí),寨上站及蘭村站的模擬精度為甲級(jí);CEEMD-BP模型的NS在0.95以上,模擬精度均為甲級(jí)。由此可知,CEEMDBP模型的預(yù)測精度最高;3種模型的預(yù)測精度都達(dá)到了乙級(jí)以上,因此模型預(yù)測效果良好。②驗(yàn)證期預(yù)測結(jié)果與模擬期相似,CEEMD-BP模型的誤差與BP模型、EMD-BP模型相比同樣大幅減小,CEEMD-BP模型的MAE、RMSE與單一BP模型相比分別減小53%~62%、48%~65%,與EMD-BP模型相比分別減小34%~46%、30%~43%。由于數(shù)據(jù)較多,因此BP神經(jīng)網(wǎng)絡(luò)出現(xiàn)了一定的過擬合現(xiàn)象,導(dǎo)致驗(yàn)證期各模型的誤差均大于對(duì)應(yīng)的模擬期誤差,但驗(yàn)證期4個(gè)水文站通過BP 單一模型預(yù)測得出的 NS 為 0.81~0.84,EMD-BP模型的 NS 為0.84~0.86,CEEMD-BP 模型的 NS 均在0.90以上,3種模型的預(yù)測精度均達(dá)到乙級(jí)以上,因此這3種模型仍然適用于非平穩(wěn)、非線性的徑流序列預(yù)測,其中CEEMD-BP模型的預(yù)測精度均達(dá)到甲級(jí),是3種模型中預(yù)測誤差最小的。

    表2 汾河上游4個(gè)水文站的預(yù)測結(jié)果

    綜上所述,無論是模擬期還是驗(yàn)證期,CEEMDBP模型的模擬預(yù)測效果都比單一BP模型及EMDBP模型的要好,且預(yù)測精度均達(dá)到甲級(jí),這進(jìn)一步說明其用于非平穩(wěn)、非線性的徑流序列預(yù)測是可行、有效的。

    3 結(jié) 論

    (1)在對(duì)非平穩(wěn)、非線性徑流序列進(jìn)行預(yù)測時(shí),單一BP模型的預(yù)測精度雖然達(dá)到乙級(jí),但是預(yù)測值與實(shí)測值易有偏差較大的情況出現(xiàn),因此需要先對(duì)徑流序列進(jìn)行平穩(wěn)化處理。

    (2)CEEMD法在保留EMD法對(duì)徑流序列平穩(wěn)化的優(yōu)勢下,解決了EMD法的模態(tài)混疊問題,從而提高了徑流處理的平穩(wěn)化程度。

    (3)CEEMD-BP模型先通過CEEMD法將原徑流序列分解成若干子序列,再通過BP神經(jīng)網(wǎng)絡(luò)對(duì)各子序列進(jìn)行預(yù)測,其預(yù)測結(jié)果與單一BP模型及EMDBP模型相比,誤差最小,預(yù)測精度最高,等級(jí)均達(dá)到甲級(jí)。因此,該方法用于復(fù)雜的非平穩(wěn)、非線性徑流序列預(yù)測是可行、有效的。

    猜你喜歡
    汾河水文站徑流
    一圖讀懂《山西省汾河保護(hù)條例》
    山西水利(2022年1期)2022-06-07 11:19:24
    SL流量計(jì)在特殊河段的應(yīng)用——以河源水文站為例
    汾河水韻
    郭家屯水文站單斷沙關(guān)系分析
    草壩水文站兩種蒸發(fā)器對(duì)比分析
    韓府灣水文站報(bào)汛方案
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    依托百里汾河工程構(gòu)建體育健身綠道的思考
    體育科技(2016年2期)2016-02-28 17:06:05
    探秘“大徑流”
    攻克“大徑流”
    日韩av在线免费看完整版不卡| 插逼视频在线观看| 国产成人aa在线观看| 婷婷色麻豆天堂久久| 国产精品一区二区三区四区免费观看| 国产伦精品一区二区三区四那| 伊人久久国产一区二区| 久久久久久久国产电影| 亚洲美女搞黄在线观看| 内地一区二区视频在线| 精品亚洲乱码少妇综合久久| 午夜精品国产一区二区电影| 国产91av在线免费观看| 高清在线视频一区二区三区| 观看免费一级毛片| 国产高清三级在线| 人人妻人人看人人澡| 99国产精品免费福利视频| 欧美+日韩+精品| 精品国产一区二区三区久久久樱花| 少妇人妻一区二区三区视频| 国产精品女同一区二区软件| 国产精品一二三区在线看| 亚洲av成人精品一区久久| 国产精品一区二区三区四区免费观看| 自拍偷自拍亚洲精品老妇| 国产成人精品婷婷| 国产免费一区二区三区四区乱码| 十八禁高潮呻吟视频 | 亚洲欧美成人综合另类久久久| 黑人猛操日本美女一级片| 最近手机中文字幕大全| 熟妇人妻不卡中文字幕| 岛国毛片在线播放| 欧美日韩视频精品一区| 国产综合精华液| 狂野欧美白嫩少妇大欣赏| av在线app专区| 十八禁网站网址无遮挡 | 99热这里只有是精品在线观看| 青春草视频在线免费观看| 免费av不卡在线播放| 国产极品粉嫩免费观看在线 | 99久久综合免费| 嘟嘟电影网在线观看| 久久狼人影院| 亚洲av欧美aⅴ国产| 国产精品国产三级国产av玫瑰| 卡戴珊不雅视频在线播放| 国产无遮挡羞羞视频在线观看| 精品一品国产午夜福利视频| freevideosex欧美| 狂野欧美激情性xxxx在线观看| 曰老女人黄片| 免费少妇av软件| 黄色日韩在线| 亚洲欧洲国产日韩| av国产精品久久久久影院| 人妻系列 视频| 十八禁高潮呻吟视频 | 亚洲国产精品成人久久小说| 久久婷婷青草| 日日啪夜夜撸| 天堂8中文在线网| 一级黄片播放器| 大又大粗又爽又黄少妇毛片口| 国产av码专区亚洲av| 精品国产一区二区久久| 一本久久精品| 亚洲精品一区蜜桃| 亚洲av二区三区四区| 国产黄频视频在线观看| 久久午夜综合久久蜜桃| 另类精品久久| 久久久久国产网址| 亚洲怡红院男人天堂| 国产av国产精品国产| 晚上一个人看的免费电影| 99re6热这里在线精品视频| 最近中文字幕2019免费版| 精品人妻一区二区三区麻豆| 亚洲不卡免费看| 在线观看人妻少妇| 日韩精品免费视频一区二区三区 | 一级二级三级毛片免费看| 91精品国产国语对白视频| 欧美日韩国产mv在线观看视频| av女优亚洲男人天堂| 久热久热在线精品观看| 国产乱人偷精品视频| a级一级毛片免费在线观看| 少妇被粗大的猛进出69影院 | 日韩一本色道免费dvd| 中文字幕人妻丝袜制服| 中文字幕人妻丝袜制服| av在线观看视频网站免费| 女人精品久久久久毛片| 成人漫画全彩无遮挡| 免费观看的影片在线观看| 韩国高清视频一区二区三区| av在线老鸭窝| 尾随美女入室| av.在线天堂| 亚洲一区二区三区欧美精品| 国产午夜精品久久久久久一区二区三区| 蜜臀久久99精品久久宅男| 国产精品福利在线免费观看| 亚洲天堂av无毛| 免费不卡的大黄色大毛片视频在线观看| 国产精品一区二区在线不卡| 国产免费视频播放在线视频| 下体分泌物呈黄色| 夜夜看夜夜爽夜夜摸| 亚洲av日韩在线播放| 免费观看a级毛片全部| 99久久中文字幕三级久久日本| 国产亚洲欧美精品永久| 有码 亚洲区| 天堂8中文在线网| 大香蕉97超碰在线| 超碰97精品在线观看| 丝袜喷水一区| 不卡视频在线观看欧美| 在线看a的网站| 久久精品久久久久久噜噜老黄| 夫妻午夜视频| 一区二区三区免费毛片| 国产精品一二三区在线看| 黄色欧美视频在线观看| 五月开心婷婷网| 日产精品乱码卡一卡2卡三| 成人影院久久| 国产欧美日韩一区二区三区在线 | 精华霜和精华液先用哪个| 一区二区三区四区激情视频| 免费看av在线观看网站| 国产高清不卡午夜福利| 男女边吃奶边做爰视频| 成人二区视频| 久久久国产精品麻豆| 26uuu在线亚洲综合色| 91精品伊人久久大香线蕉| 自线自在国产av| 亚洲成色77777| 国产精品国产三级专区第一集| 一区二区av电影网| 女人精品久久久久毛片| 久久久国产一区二区| 黑人猛操日本美女一级片| 男男h啪啪无遮挡| a级毛色黄片| 日韩电影二区| 波野结衣二区三区在线| 亚洲精品视频女| 亚洲av.av天堂| 校园人妻丝袜中文字幕| 国产 精品1| 亚洲av在线观看美女高潮| 久久狼人影院| 国产有黄有色有爽视频| 国产极品粉嫩免费观看在线 | 亚洲欧洲国产日韩| 国产老妇伦熟女老妇高清| h视频一区二区三区| 亚洲精品日韩在线中文字幕| 亚洲,一卡二卡三卡| 亚洲人成网站在线播| 免费av不卡在线播放| 爱豆传媒免费全集在线观看| 国产 精品1| 精品人妻熟女毛片av久久网站| av福利片在线| 亚洲精品日本国产第一区| 熟女av电影| 久久精品国产自在天天线| 国模一区二区三区四区视频| a级一级毛片免费在线观看| 婷婷色麻豆天堂久久| 黑人猛操日本美女一级片| 欧美3d第一页| 久久国产精品大桥未久av | 18+在线观看网站| 国产黄频视频在线观看| 一本大道久久a久久精品| 日韩欧美一区视频在线观看 | 久久久久久久国产电影| 久久久久久久久久久久大奶| 国产成人a∨麻豆精品| 免费少妇av软件| 午夜日本视频在线| 亚洲婷婷狠狠爱综合网| 91精品一卡2卡3卡4卡| 哪个播放器可以免费观看大片| 国产免费一区二区三区四区乱码| 尾随美女入室| 日韩中文字幕视频在线看片| 18禁动态无遮挡网站| 大陆偷拍与自拍| av免费观看日本| 国产一区二区在线观看av| 一级毛片 在线播放| 9色porny在线观看| 色婷婷久久久亚洲欧美| 99热这里只有是精品50| 夜夜看夜夜爽夜夜摸| 中国三级夫妇交换| 性色avwww在线观看| 亚洲av在线观看美女高潮| 黄色日韩在线| 成人黄色视频免费在线看| 人人妻人人添人人爽欧美一区卜| 秋霞在线观看毛片| 日韩成人伦理影院| 国产成人精品无人区| 九九久久精品国产亚洲av麻豆| 交换朋友夫妻互换小说| 婷婷色麻豆天堂久久| 久久久精品94久久精品| 日本av手机在线免费观看| 久久久国产欧美日韩av| 亚洲一区二区三区欧美精品| av国产精品久久久久影院| 十八禁高潮呻吟视频 | 国产极品粉嫩免费观看在线 | av视频免费观看在线观看| 久久久久精品久久久久真实原创| 中文字幕精品免费在线观看视频 | 中文天堂在线官网| 亚洲精品久久午夜乱码| 亚洲图色成人| 最新的欧美精品一区二区| av国产久精品久网站免费入址| 少妇人妻精品综合一区二区| 两个人的视频大全免费| 免费观看性生交大片5| 国产高清不卡午夜福利| 精品久久久久久电影网| 国产一级毛片在线| 亚洲欧美精品自产自拍| 久久99热这里只频精品6学生| videos熟女内射| 欧美xxⅹ黑人| 亚洲精品久久午夜乱码| 亚洲综合精品二区| 在线观看www视频免费| 欧美成人午夜免费资源| 国产男人的电影天堂91| 亚洲av.av天堂| 亚洲国产精品专区欧美| 国产男人的电影天堂91| 天堂8中文在线网| 嫩草影院新地址| 国产精品99久久99久久久不卡 | 国产精品99久久久久久久久| 99视频精品全部免费 在线| 最黄视频免费看| 一二三四中文在线观看免费高清| 亚洲欧美中文字幕日韩二区| 亚洲av在线观看美女高潮| 在线精品无人区一区二区三| 中文天堂在线官网| 18禁在线播放成人免费| 丰满饥渴人妻一区二区三| 免费人成在线观看视频色| 在线天堂最新版资源| 在线观看国产h片| 女性被躁到高潮视频| 建设人人有责人人尽责人人享有的| 高清毛片免费看| 久久婷婷青草| 国产精品国产av在线观看| 成年美女黄网站色视频大全免费 | 亚洲欧美成人精品一区二区| 全区人妻精品视频| 亚洲精品日韩av片在线观看| 亚洲欧洲国产日韩| 亚洲自偷自拍三级| 黑人高潮一二区| 国产黄片视频在线免费观看| 男人和女人高潮做爰伦理| 久久久久久久国产电影| 不卡视频在线观看欧美| 国产精品秋霞免费鲁丝片| 成人特级av手机在线观看| 亚洲欧美中文字幕日韩二区| 少妇的逼水好多| 国产男女超爽视频在线观看| 国产精品偷伦视频观看了| 91久久精品国产一区二区三区| 久久国内精品自在自线图片| 精品视频人人做人人爽| 高清黄色对白视频在线免费看 | 亚洲色图综合在线观看| 九九爱精品视频在线观看| 国产熟女欧美一区二区| 青春草亚洲视频在线观看| 日本欧美视频一区| 极品人妻少妇av视频| 女的被弄到高潮叫床怎么办| 纵有疾风起免费观看全集完整版| 国内揄拍国产精品人妻在线| 国产日韩欧美视频二区| 少妇高潮的动态图| 久久99热这里只频精品6学生| 亚洲国产精品国产精品| 一区在线观看完整版| av国产久精品久网站免费入址| 观看免费一级毛片| 免费高清在线观看视频在线观看| xxx大片免费视频| 在线观看免费高清a一片| 一区二区三区乱码不卡18| 2021少妇久久久久久久久久久| 精品人妻熟女av久视频| 国产色婷婷99| 女人久久www免费人成看片| 日本色播在线视频| 最近中文字幕高清免费大全6| av在线app专区| 亚洲中文av在线| 日韩制服骚丝袜av| 国产av国产精品国产| 在线 av 中文字幕| 中文字幕制服av| 性色avwww在线观看| 热re99久久精品国产66热6| 亚洲国产日韩一区二区| 亚洲精品乱久久久久久| 99久久精品国产国产毛片| 国产淫语在线视频| 国产一区二区三区av在线| 青春草视频在线免费观看| 亚洲高清免费不卡视频| 91在线精品国自产拍蜜月| 麻豆精品久久久久久蜜桃| 亚洲第一av免费看| 国产有黄有色有爽视频| 国产精品三级大全| 最近2019中文字幕mv第一页| 80岁老熟妇乱子伦牲交| 我的老师免费观看完整版| 午夜免费鲁丝| 香蕉精品网在线| 久久鲁丝午夜福利片| 看非洲黑人一级黄片| 一本一本综合久久| 黄色欧美视频在线观看| 国产永久视频网站| 精品久久国产蜜桃| 久久久久久久国产电影| 黄色日韩在线| 大码成人一级视频| 成年人午夜在线观看视频| 国产黄片美女视频| 国产视频首页在线观看| 美女中出高潮动态图| 亚洲精品色激情综合| 日本与韩国留学比较| 色婷婷久久久亚洲欧美| 欧美精品人与动牲交sv欧美| 男人爽女人下面视频在线观看| 三级经典国产精品| 亚洲色图综合在线观看| 国产伦理片在线播放av一区| 狂野欧美激情性xxxx在线观看| 免费少妇av软件| 午夜av观看不卡| 日韩一区二区三区影片| 亚洲精品自拍成人| 啦啦啦啦在线视频资源| av.在线天堂| 免费观看的影片在线观看| 岛国毛片在线播放| 卡戴珊不雅视频在线播放| 亚洲精品一二三| 成人影院久久| 国产在线男女| 国产亚洲最大av| 久久久a久久爽久久v久久| 欧美bdsm另类| 亚洲美女搞黄在线观看| 亚洲精品,欧美精品| 男人爽女人下面视频在线观看| 人妻夜夜爽99麻豆av| 亚洲成人手机| 中文资源天堂在线| 在线观看美女被高潮喷水网站| 男人添女人高潮全过程视频| 国产成人免费观看mmmm| 国产精品麻豆人妻色哟哟久久| 亚洲av成人精品一二三区| 日韩一区二区三区影片| 在线看a的网站| 国产免费视频播放在线视频| 国产精品国产av在线观看| 在线观看国产h片| 国产色爽女视频免费观看| 哪个播放器可以免费观看大片| 日韩av免费高清视频| 欧美日韩视频高清一区二区三区二| 一区在线观看完整版| 丰满人妻一区二区三区视频av| 国产高清有码在线观看视频| 最近的中文字幕免费完整| 久久久a久久爽久久v久久| 女人久久www免费人成看片| 美女国产视频在线观看| 一本一本综合久久| 秋霞伦理黄片| 国产精品成人在线| 国产免费视频播放在线视频| 中文在线观看免费www的网站| 国产成人精品无人区| 2018国产大陆天天弄谢| 久久久久久久国产电影| 成人亚洲精品一区在线观看| 99热网站在线观看| 国产精品熟女久久久久浪| 亚洲人成网站在线观看播放| 欧美精品一区二区免费开放| 国产综合精华液| 亚洲人成网站在线播| 交换朋友夫妻互换小说| 精品久久久精品久久久| 欧美日韩一区二区视频在线观看视频在线| 啦啦啦视频在线资源免费观看| 国产精品人妻久久久久久| 亚洲av电影在线观看一区二区三区| 精品视频人人做人人爽| 久久久国产一区二区| 久久久精品免费免费高清| 亚洲av福利一区| 亚洲精品乱码久久久v下载方式| 桃花免费在线播放| 免费看av在线观看网站| 只有这里有精品99| 99热这里只有是精品在线观看| 日本av免费视频播放| 日韩 亚洲 欧美在线| 大话2 男鬼变身卡| 曰老女人黄片| 日本-黄色视频高清免费观看| 性色avwww在线观看| 国产成人精品久久久久久| 欧美日韩国产mv在线观看视频| 99re6热这里在线精品视频| 美女福利国产在线| 亚洲经典国产精华液单| 2021少妇久久久久久久久久久| 国产黄片美女视频| 精品一区二区三区视频在线| 欧美xxⅹ黑人| a级毛色黄片| 久久99蜜桃精品久久| videossex国产| 久久精品国产鲁丝片午夜精品| 国产乱来视频区| 久久这里有精品视频免费| 亚洲av中文av极速乱| 久久久久人妻精品一区果冻| 亚洲国产精品国产精品| 天天躁夜夜躁狠狠久久av| 亚洲欧洲精品一区二区精品久久久 | 亚州av有码| 人妻系列 视频| 免费观看在线日韩| 国产免费视频播放在线视频| 亚洲国产毛片av蜜桃av| 69精品国产乱码久久久| 久久99热6这里只有精品| 欧美最新免费一区二区三区| 一二三四中文在线观看免费高清| 观看免费一级毛片| 一级,二级,三级黄色视频| 我要看黄色一级片免费的| 国产又色又爽无遮挡免| 视频区图区小说| 国产无遮挡羞羞视频在线观看| 老司机影院毛片| 王馨瑶露胸无遮挡在线观看| 高清在线视频一区二区三区| 国产av码专区亚洲av| 免费观看无遮挡的男女| 丝瓜视频免费看黄片| 国产精品一区二区性色av| 国产成人a∨麻豆精品| 国产一区有黄有色的免费视频| 天天躁夜夜躁狠狠久久av| 中文字幕av电影在线播放| 在线观看一区二区三区激情| 国产精品国产三级国产av玫瑰| 免费看av在线观看网站| 国产探花极品一区二区| 成人特级av手机在线观看| 少妇熟女欧美另类| 久久国产精品大桥未久av | 亚洲成人手机| 国产成人精品一,二区| 午夜日本视频在线| 成人18禁高潮啪啪吃奶动态图 | 如日韩欧美国产精品一区二区三区 | 欧美少妇被猛烈插入视频| 简卡轻食公司| 免费看不卡的av| 国产视频首页在线观看| 久久精品国产鲁丝片午夜精品| 久久国产乱子免费精品| 国产在视频线精品| 亚洲无线观看免费| 日韩成人av中文字幕在线观看| 免费看av在线观看网站| 女人久久www免费人成看片| 日日摸夜夜添夜夜爱| 精品人妻一区二区三区麻豆| 99热国产这里只有精品6| 视频中文字幕在线观看| 韩国高清视频一区二区三区| 美女大奶头黄色视频| 亚洲av电影在线观看一区二区三区| 日本黄色片子视频| 边亲边吃奶的免费视频| 免费观看性生交大片5| 少妇的逼水好多| 亚洲av中文av极速乱| 精品少妇久久久久久888优播| videos熟女内射| 亚洲欧美成人综合另类久久久| 精品久久久精品久久久| 乱码一卡2卡4卡精品| 91久久精品电影网| 久久青草综合色| 在线精品无人区一区二区三| 亚洲熟女精品中文字幕| 亚洲色图综合在线观看| 亚洲av不卡在线观看| 自拍偷自拍亚洲精品老妇| 啦啦啦中文免费视频观看日本| 80岁老熟妇乱子伦牲交| √禁漫天堂资源中文www| 涩涩av久久男人的天堂| 国产亚洲91精品色在线| 伊人久久国产一区二区| 欧美精品国产亚洲| 国产乱人偷精品视频| 国产亚洲91精品色在线| 国产在线一区二区三区精| 一本色道久久久久久精品综合| 18禁动态无遮挡网站| 国产精品福利在线免费观看| 亚洲人成网站在线播| 日日啪夜夜撸| 国产精品久久久久久久久免| 亚洲精品乱码久久久久久按摩| 亚洲av综合色区一区| 欧美日韩在线观看h| 国产爽快片一区二区三区| 街头女战士在线观看网站| 99热这里只有精品一区| 99热国产这里只有精品6| 国产色爽女视频免费观看| 最黄视频免费看| 亚洲av在线观看美女高潮| 国产极品天堂在线| 建设人人有责人人尽责人人享有的| 色网站视频免费| 久久影院123| 亚洲在久久综合| 狂野欧美激情性xxxx在线观看| 欧美高清成人免费视频www| 一级毛片 在线播放| 国产日韩欧美亚洲二区| 中文欧美无线码| 亚洲欧美中文字幕日韩二区| 街头女战士在线观看网站| 好男人视频免费观看在线| 肉色欧美久久久久久久蜜桃| 久久97久久精品| 欧美区成人在线视频| 视频中文字幕在线观看| 成人亚洲精品一区在线观看| 最近中文字幕高清免费大全6| 精品国产一区二区三区久久久樱花| 在线观看美女被高潮喷水网站| 美女中出高潮动态图| av免费在线看不卡| 午夜免费观看性视频| 九色成人免费人妻av| 国产毛片在线视频| 精品久久国产蜜桃| 一个人免费看片子| 伦理电影免费视频| 欧美一级a爱片免费观看看| 91久久精品国产一区二区三区| 久久久精品免费免费高清| 精品少妇久久久久久888优播| 国产淫语在线视频| 人人澡人人妻人| 80岁老熟妇乱子伦牲交| 精品少妇久久久久久888优播| 18+在线观看网站| 亚洲色图综合在线观看| 国产精品国产av在线观看| 国产综合精华液| 熟女av电影| 国产黄片视频在线免费观看| 91成人精品电影| 少妇丰满av| 三上悠亚av全集在线观看 | 少妇猛男粗大的猛烈进出视频| 久久av网站| 亚洲国产最新在线播放| 免费观看在线日韩| 人妻系列 视频| 成人毛片a级毛片在线播放| 免费看av在线观看网站|