• <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涩爱| 亚洲经典国产精华液单| 99视频精品全部免费 在线| 寂寞人妻少妇视频99o| 日韩人妻高清精品专区| 国产伦精品一区二区三区视频9| 亚洲内射少妇av| 久久久精品94久久精品| 下体分泌物呈黄色| 国产爽快片一区二区三区| 免费在线观看成人毛片| 大又大粗又爽又黄少妇毛片口| 日韩,欧美,国产一区二区三区| 黑丝袜美女国产一区| 天天操日日干夜夜撸| 中文欧美无线码| 桃花免费在线播放| 久久99热这里只频精品6学生| 视频中文字幕在线观看| 免费看av在线观看网站| 久久精品熟女亚洲av麻豆精品| 十八禁网站网址无遮挡 | 亚洲欧洲日产国产| 免费久久久久久久精品成人欧美视频 | 一级a做视频免费观看| 成人国产av品久久久| 99热6这里只有精品| 欧美 日韩 精品 国产| 精品亚洲乱码少妇综合久久| 日韩成人伦理影院| 精品熟女少妇av免费看| 深夜a级毛片| 97超视频在线观看视频| 啦啦啦视频在线资源免费观看| 精品国产国语对白av| 91精品国产九色| 老女人水多毛片| 亚洲,欧美,日韩| 久久人妻熟女aⅴ| 久久久午夜欧美精品| 日本猛色少妇xxxxx猛交久久| 日日爽夜夜爽网站| 少妇 在线观看| 伊人久久国产一区二区| 欧美高清成人免费视频www| 日本免费在线观看一区| 国产伦精品一区二区三区视频9| 亚洲欧洲精品一区二区精品久久久 | 九九在线视频观看精品| 久久99蜜桃精品久久| 黑人猛操日本美女一级片| 成人无遮挡网站| a级毛色黄片| 一个人免费看片子| 精品人妻熟女毛片av久久网站| 亚洲av中文av极速乱| av免费在线看不卡| 国产成人免费无遮挡视频| 午夜激情福利司机影院| 熟妇人妻不卡中文字幕| 亚洲欧美成人精品一区二区| 中文天堂在线官网| 99久久精品热视频| 久久6这里有精品| 成人二区视频| 国产毛片在线视频| 日日啪夜夜撸| 一级毛片电影观看| 男人狂女人下面高潮的视频| 少妇被粗大的猛进出69影院 | 亚洲欧美一区二区三区国产| 在线观看一区二区三区激情| 99国产精品免费福利视频| 婷婷色麻豆天堂久久| 最近最新中文字幕免费大全7| 好男人视频免费观看在线| 韩国高清视频一区二区三区| 免费看不卡的av| 美女脱内裤让男人舔精品视频| 亚洲综合色惰| 午夜91福利影院| 另类精品久久| 日韩成人伦理影院| 永久网站在线| 搡女人真爽免费视频火全软件| 欧美少妇被猛烈插入视频| 18禁在线无遮挡免费观看视频| 国产亚洲一区二区精品| 欧美激情国产日韩精品一区| 只有这里有精品99| 国产精品一区二区在线不卡| 国产免费一区二区三区四区乱码| 国产色婷婷99| 精品国产一区二区三区久久久樱花| 亚洲精品久久久久久婷婷小说| 七月丁香在线播放| 国产精品熟女久久久久浪| 自线自在国产av| 国产毛片在线视频| 日本av免费视频播放| 色视频在线一区二区三区| 国产 精品1| tube8黄色片| 亚洲中文av在线| 久久精品久久久久久噜噜老黄| 欧美性感艳星| 秋霞在线观看毛片| 中国美白少妇内射xxxbb| 国产午夜精品久久久久久一区二区三区| 国产亚洲91精品色在线| 一本久久精品| 在线亚洲精品国产二区图片欧美 | 两个人的视频大全免费| 桃花免费在线播放| 亚洲精品国产av成人精品| 久久久国产欧美日韩av| 九九在线视频观看精品| 国产伦精品一区二区三区视频9| 天美传媒精品一区二区| 国产精品国产三级国产专区5o| 尾随美女入室| 久久国产乱子免费精品| 久久精品久久精品一区二区三区| 成人无遮挡网站| 99久久人妻综合| 久久鲁丝午夜福利片| 国产亚洲5aaaaa淫片| 七月丁香在线播放| 热99国产精品久久久久久7| 精品亚洲乱码少妇综合久久| 日本-黄色视频高清免费观看| www.色视频.com| 亚洲色图综合在线观看| 老熟女久久久| 欧美激情极品国产一区二区三区 | 黄色欧美视频在线观看| 国内少妇人妻偷人精品xxx网站| 精品久久久久久电影网| 亚洲国产最新在线播放| 久久久久久人妻| 久久精品国产亚洲av涩爱| 在线观看免费日韩欧美大片 | av不卡在线播放| videos熟女内射| 亚洲经典国产精华液单| 国产精品免费大片| 人妻系列 视频| 伊人久久精品亚洲午夜| 亚洲精品色激情综合| 丝袜脚勾引网站| 免费看不卡的av| 久久国产乱子免费精品| 老司机亚洲免费影院| 久热这里只有精品99| 免费高清在线观看视频在线观看| 一区二区三区免费毛片| 午夜免费男女啪啪视频观看| 又大又黄又爽视频免费| 大香蕉97超碰在线| 伦精品一区二区三区| 久久精品国产自在天天线| 日韩欧美精品免费久久| 久久av网站| 日韩不卡一区二区三区视频在线| 亚洲四区av| 久久毛片免费看一区二区三区| 欧美亚洲 丝袜 人妻 在线| 一级二级三级毛片免费看| a 毛片基地| 夜夜爽夜夜爽视频| 亚洲激情五月婷婷啪啪| 国产成人a∨麻豆精品| 黑人巨大精品欧美一区二区蜜桃 | 男人爽女人下面视频在线观看| 国产高清三级在线| 熟女电影av网| 免费av中文字幕在线| 99久久精品国产国产毛片| 国产精品久久久久久久久免| 国产av国产精品国产| 女人精品久久久久毛片| 视频中文字幕在线观看| 大香蕉97超碰在线| 亚洲精品国产av蜜桃| 日本91视频免费播放| 久久久欧美国产精品| 高清午夜精品一区二区三区| 国产极品天堂在线| 国内揄拍国产精品人妻在线| 男女边吃奶边做爰视频| 一级a做视频免费观看| 男男h啪啪无遮挡| 男人爽女人下面视频在线观看| 亚洲不卡免费看| 久久久精品免费免费高清| 大香蕉久久网| √禁漫天堂资源中文www| 观看美女的网站| 五月天丁香电影| 亚洲国产精品国产精品| 91精品一卡2卡3卡4卡| 国产白丝娇喘喷水9色精品| 亚洲精品久久久久久婷婷小说| 又粗又硬又长又爽又黄的视频| 日日爽夜夜爽网站| 国产老妇伦熟女老妇高清| 免费少妇av软件| 国产av一区二区精品久久| 免费看不卡的av| av福利片在线观看| av播播在线观看一区| 亚洲精品乱久久久久久| 国产精品熟女久久久久浪| 人妻一区二区av| 久久精品国产自在天天线| 国产免费一区二区三区四区乱码| 久久精品熟女亚洲av麻豆精品| 哪个播放器可以免费观看大片| 免费观看无遮挡的男女| 亚洲不卡免费看| 免费av不卡在线播放| 色哟哟·www| 人妻系列 视频| 国产成人freesex在线| 日韩大片免费观看网站| 狂野欧美激情性bbbbbb| 久久鲁丝午夜福利片| 日韩av在线免费看完整版不卡| 亚洲欧美日韩卡通动漫| 国产精品麻豆人妻色哟哟久久| 三级国产精品欧美在线观看| 亚洲高清免费不卡视频| 又爽又黄a免费视频| 中文字幕人妻丝袜制服| 日韩中文字幕视频在线看片| 国产91av在线免费观看| 久久99一区二区三区| 免费av中文字幕在线| 欧美成人午夜免费资源| 国产免费一级a男人的天堂| 狂野欧美激情性xxxx在线观看| 麻豆乱淫一区二区| 午夜福利网站1000一区二区三区| 欧美日韩视频高清一区二区三区二| 99热这里只有精品一区| 国产av码专区亚洲av| 欧美精品国产亚洲| 老司机影院成人| 亚洲久久久国产精品| 久久久久精品久久久久真实原创| 国产淫片久久久久久久久| 五月开心婷婷网| 久久精品久久久久久噜噜老黄| 欧美少妇被猛烈插入视频| 国产精品久久久久久精品电影小说| 高清在线视频一区二区三区| 十八禁网站网址无遮挡 | 69精品国产乱码久久久| 大码成人一级视频| 精品少妇内射三级| 婷婷色综合www| 亚洲av不卡在线观看| 在线观看免费视频网站a站| 大陆偷拍与自拍| 国产精品国产三级专区第一集| 制服丝袜香蕉在线| 亚洲av日韩在线播放| 制服丝袜香蕉在线| 免费av不卡在线播放| 国产在线视频一区二区| 伦理电影大哥的女人| 免费看光身美女| 久久久欧美国产精品| 男人添女人高潮全过程视频| 精品人妻一区二区三区麻豆| 久久精品国产鲁丝片午夜精品| 国产91av在线免费观看| 久久婷婷青草| 亚洲国产日韩一区二区| 亚洲国产精品专区欧美| 国产av国产精品国产| 亚洲欧美日韩另类电影网站| 日本欧美视频一区| av福利片在线观看| 国产一区二区在线观看日韩| 黑人猛操日本美女一级片| 人人澡人人妻人| 国产片特级美女逼逼视频| 亚洲av在线观看美女高潮| 欧美激情国产日韩精品一区| 欧美日韩在线观看h| 观看美女的网站| 久久久久久久久久久丰满| 日本-黄色视频高清免费观看| 欧美另类一区| 极品少妇高潮喷水抽搐| 国产精品偷伦视频观看了| 黄色毛片三级朝国网站 | 午夜福利网站1000一区二区三区| 男男h啪啪无遮挡| 国产淫语在线视频| 最新中文字幕久久久久| 免费黄色在线免费观看| 亚洲经典国产精华液单| 丰满迷人的少妇在线观看| 国产日韩欧美视频二区| 日日摸夜夜添夜夜添av毛片| 精品久久久久久电影网| 亚洲欧美一区二区三区国产| www.av在线官网国产| 爱豆传媒免费全集在线观看| 99热全是精品| .国产精品久久| 国产精品三级大全| 亚洲图色成人| 男的添女的下面高潮视频| 最近中文字幕高清免费大全6| 涩涩av久久男人的天堂| 一区二区三区精品91| 美女福利国产在线| 久久毛片免费看一区二区三区| 最新中文字幕久久久久| 97精品久久久久久久久久精品| 日韩,欧美,国产一区二区三区| 亚洲av.av天堂| 国产中年淑女户外野战色| 国产一区二区三区综合在线观看 | 欧美日韩综合久久久久久| 免费久久久久久久精品成人欧美视频 | 香蕉精品网在线| 成人毛片60女人毛片免费| 久久综合国产亚洲精品| 在线观看一区二区三区激情| 国产欧美日韩综合在线一区二区 | 精品99又大又爽又粗少妇毛片| 天堂中文最新版在线下载| 亚洲国产精品一区三区| 精品久久久噜噜| 日本av免费视频播放| 春色校园在线视频观看| 国产在线视频一区二区| 亚洲欧洲日产国产| 嘟嘟电影网在线观看| 欧美日韩一区二区视频在线观看视频在线| 天堂8中文在线网| 国产精品麻豆人妻色哟哟久久| 国产精品无大码| 2021少妇久久久久久久久久久| 国产精品99久久久久久久久| 日韩三级伦理在线观看| 丝袜在线中文字幕| 成人毛片a级毛片在线播放| 日日爽夜夜爽网站| 成年女人在线观看亚洲视频| 最新中文字幕久久久久| 秋霞在线观看毛片| 在线天堂最新版资源| 亚洲av电影在线观看一区二区三区| 中文在线观看免费www的网站| 国产黄频视频在线观看| 久久久久网色| 色94色欧美一区二区| 精品久久久久久久久av| 国产精品伦人一区二区| 国产亚洲5aaaaa淫片| 一本久久精品| 伊人久久精品亚洲午夜| 激情五月婷婷亚洲| 国产高清三级在线| 深夜a级毛片| 国产片特级美女逼逼视频| 精品酒店卫生间| av天堂久久9| 免费看av在线观看网站| 久久热精品热| 国产成人免费无遮挡视频| 午夜免费男女啪啪视频观看| 在线观看三级黄色| 高清毛片免费看| 欧美bdsm另类| 日韩一区二区三区影片| 亚洲av欧美aⅴ国产| 亚洲高清免费不卡视频| 欧美老熟妇乱子伦牲交| 精品一区在线观看国产| 一级毛片aaaaaa免费看小| 中文欧美无线码| 最近中文字幕高清免费大全6| 亚洲久久久国产精品| 高清av免费在线| 丝瓜视频免费看黄片| 亚洲av国产av综合av卡| 最近手机中文字幕大全| 看免费成人av毛片| 中国国产av一级| 午夜影院在线不卡| 99久久精品热视频| 91精品国产九色| a级毛片在线看网站| 日韩欧美精品免费久久| 国产在线免费精品| 精品国产国语对白av| 欧美激情国产日韩精品一区| 免费久久久久久久精品成人欧美视频 | 日本欧美视频一区| 国产伦精品一区二区三区视频9| 国产亚洲一区二区精品| 日韩一本色道免费dvd| 亚洲精品国产av成人精品| 国产成人精品福利久久| 一级毛片久久久久久久久女| 精品一区二区三卡| 桃花免费在线播放| 99热这里只有精品一区| 七月丁香在线播放| 18禁动态无遮挡网站| 亚洲精品日韩在线中文字幕| 国产精品久久久久成人av| 伊人久久精品亚洲午夜| 成人影院久久| 精品久久久噜噜| 99热网站在线观看| 蜜桃久久精品国产亚洲av| 亚洲欧美成人综合另类久久久| 国产亚洲一区二区精品| 日日啪夜夜撸| 下体分泌物呈黄色| 久久久久久久久久人人人人人人| 国产在视频线精品| 女人久久www免费人成看片| 中国国产av一级| 国产伦精品一区二区三区四那| 十分钟在线观看高清视频www | 欧美日韩综合久久久久久| 黄色怎么调成土黄色| 啦啦啦啦在线视频资源| 夜夜爽夜夜爽视频| 99热网站在线观看| 亚洲av在线观看美女高潮| 精华霜和精华液先用哪个| 午夜老司机福利剧场| 成人免费观看视频高清| 国产高清有码在线观看视频| 99国产精品免费福利视频| 国产精品欧美亚洲77777| tube8黄色片| 99久久精品热视频| 中文字幕免费在线视频6| 亚洲天堂av无毛| 51国产日韩欧美| 午夜福利网站1000一区二区三区| 国产精品嫩草影院av在线观看| 在线观看av片永久免费下载| 在线 av 中文字幕| 两个人的视频大全免费| 男人和女人高潮做爰伦理| 久久人妻熟女aⅴ| 性色avwww在线观看| 一个人看视频在线观看www免费| 美女福利国产在线| 五月开心婷婷网| 久久精品久久精品一区二区三区| 国产精品.久久久| av国产精品久久久久影院| 午夜福利,免费看| 曰老女人黄片| 插逼视频在线观看| av女优亚洲男人天堂| 日本-黄色视频高清免费观看| 插阴视频在线观看视频| 岛国毛片在线播放| 久久久久精品性色| 热99国产精品久久久久久7| 久久热精品热| 亚洲欧美一区二区三区国产| 欧美精品一区二区免费开放| 啦啦啦啦在线视频资源| 亚洲真实伦在线观看| 中文字幕亚洲精品专区| 国产一区二区在线观看日韩| 如何舔出高潮| 99久久精品热视频| 国产日韩欧美亚洲二区| 国产精品偷伦视频观看了| 亚洲人成网站在线播| 少妇被粗大猛烈的视频| 欧美日韩亚洲高清精品| 国产成人91sexporn| 乱人伦中国视频| 99热国产这里只有精品6| 中文在线观看免费www的网站| www.色视频.com| 日日啪夜夜爽| 国产亚洲欧美精品永久| 精品一区二区三卡| 多毛熟女@视频| 精品少妇黑人巨大在线播放| 欧美日韩视频高清一区二区三区二| 亚洲精品一二三| 国产免费视频播放在线视频| 大香蕉久久网| av天堂久久9| 免费看不卡的av| 一级a做视频免费观看| 又爽又黄a免费视频| 多毛熟女@视频| 少妇被粗大的猛进出69影院 | 精品亚洲乱码少妇综合久久| 精品国产一区二区久久| 国产深夜福利视频在线观看| 久久久久久伊人网av| 一个人免费看片子| 大码成人一级视频| 99re6热这里在线精品视频| 嫩草影院新地址| 伊人久久精品亚洲午夜| 一区二区三区四区激情视频| 国产成人精品福利久久| 国产精品久久久久久av不卡| 三级经典国产精品| 精品少妇内射三级| 欧美日韩一区二区视频在线观看视频在线| 久久人人爽人人片av| 日本av免费视频播放| 亚洲国产成人一精品久久久| 国产亚洲91精品色在线| 伊人久久精品亚洲午夜| 老司机影院成人| 少妇人妻一区二区三区视频| 午夜影院在线不卡| 日韩成人av中文字幕在线观看| 色网站视频免费| 一级黄片播放器| www.色视频.com| 久久久欧美国产精品| 六月丁香七月| a级毛片在线看网站| 久久久久久人妻| 日本猛色少妇xxxxx猛交久久| 成年人免费黄色播放视频 | 亚洲无线观看免费| 日韩免费高清中文字幕av| 欧美日韩综合久久久久久| 能在线免费看毛片的网站| 亚洲精品日韩在线中文字幕| 啦啦啦啦在线视频资源| 久久精品国产自在天天线| 亚洲美女搞黄在线观看| 久久久久视频综合| 两个人免费观看高清视频 | 精品视频人人做人人爽| 97在线视频观看| 亚洲成色77777| 2021少妇久久久久久久久久久| 国产白丝娇喘喷水9色精品| 偷拍熟女少妇极品色| 精品一区二区三卡| freevideosex欧美| 简卡轻食公司| 天天操日日干夜夜撸| 久久久久国产精品人妻一区二区| 国产精品福利在线免费观看| 久久久久久久国产电影| 国产精品国产三级国产av玫瑰| 欧美高清成人免费视频www| 91在线精品国自产拍蜜月| 久久女婷五月综合色啪小说| 亚洲综合精品二区| 精品人妻偷拍中文字幕| 国产无遮挡羞羞视频在线观看| 男的添女的下面高潮视频| 精品人妻熟女av久视频| 精品久久久久久电影网| 免费少妇av软件| a级毛片免费高清观看在线播放| 麻豆成人av视频| 热re99久久国产66热| 欧美三级亚洲精品| 晚上一个人看的免费电影| 日韩制服骚丝袜av| 中文字幕久久专区| 成人漫画全彩无遮挡| 成人免费观看视频高清| 我要看黄色一级片免费的| 亚洲av日韩在线播放| 欧美xxxx性猛交bbbb| 亚洲美女黄色视频免费看| 高清毛片免费看| 在线观看三级黄色| 久久亚洲国产成人精品v| 美女福利国产在线| 熟女电影av网| 久热久热在线精品观看| 精品一区二区免费观看| 久久国产精品男人的天堂亚洲 | 国产欧美亚洲国产| 人人澡人人妻人| 黄色配什么色好看| 少妇 在线观看| 观看av在线不卡| 精品亚洲成国产av| 一级a做视频免费观看| 99国产精品免费福利视频| 丝袜脚勾引网站| 免费观看a级毛片全部| 精品人妻一区二区三区麻豆| 日韩成人av中文字幕在线观看| 欧美日韩综合久久久久久| 国产精品久久久久成人av| 这个男人来自地球电影免费观看 | 免费av中文字幕在线| 成人亚洲精品一区在线观看| 亚洲婷婷狠狠爱综合网| 久久女婷五月综合色啪小说| 日韩欧美精品免费久久|