• <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
    探秘“大徑流”
    攻克“大徑流”
    热99国产精品久久久久久7| 99久久综合免费| 青春草亚洲视频在线观看| 最近中文字幕高清免费大全6| 春色校园在线视频观看| 精品少妇一区二区三区视频日本电影 | 午夜福利在线免费观看网站| 国产高清不卡午夜福利| 视频在线观看一区二区三区| 99久久中文字幕三级久久日本| av天堂久久9| 国产1区2区3区精品| 在现免费观看毛片| 久久久久精品性色| 亚洲精品日本国产第一区| 亚洲四区av| 最近的中文字幕免费完整| 国产精品国产三级国产专区5o| 麻豆乱淫一区二区| 一二三四在线观看免费中文在| 亚洲五月色婷婷综合| 国产一级毛片在线| 亚洲少妇的诱惑av| 日韩一区二区视频免费看| 亚洲成av片中文字幕在线观看 | 国产国语露脸激情在线看| 精品久久蜜臀av无| 国产日韩一区二区三区精品不卡| 啦啦啦在线免费观看视频4| 性色avwww在线观看| 熟妇人妻不卡中文字幕| 熟女电影av网| 国产男人的电影天堂91| 婷婷色综合www| av在线app专区| 亚洲人成77777在线视频| 久久精品久久久久久噜噜老黄| xxxhd国产人妻xxx| 欧美成人午夜免费资源| 国产一区二区激情短视频 | 日本欧美视频一区| 午夜福利网站1000一区二区三区| 亚洲精品国产色婷婷电影| 狂野欧美激情性bbbbbb| 中文字幕人妻丝袜一区二区 | 亚洲精品久久午夜乱码| 婷婷成人精品国产| 国产精品一区二区在线观看99| 亚洲欧洲国产日韩| 久久国产亚洲av麻豆专区| 如何舔出高潮| 国产高清不卡午夜福利| 大片电影免费在线观看免费| 国产成人aa在线观看| 午夜免费观看性视频| 天天操日日干夜夜撸| av片东京热男人的天堂| 欧美人与性动交α欧美精品济南到 | 丰满少妇做爰视频| 国产在视频线精品| 人妻系列 视频| 伊人久久国产一区二区| 久久久精品免费免费高清| 老汉色av国产亚洲站长工具| 天堂中文最新版在线下载| 黄片无遮挡物在线观看| 亚洲伊人久久精品综合| 中国国产av一级| 王馨瑶露胸无遮挡在线观看| 女性被躁到高潮视频| 欧美97在线视频| 人成视频在线观看免费观看| 熟女少妇亚洲综合色aaa.| 黄色视频在线播放观看不卡| 成年动漫av网址| 男女啪啪激烈高潮av片| 精品第一国产精品| 熟妇人妻不卡中文字幕| 国产成人精品无人区| 日韩av免费高清视频| 欧美日韩av久久| 赤兔流量卡办理| 亚洲人成电影观看| av在线播放精品| 18禁观看日本| 中文字幕人妻丝袜制服| 亚洲综合色网址| 男人爽女人下面视频在线观看| 美女午夜性视频免费| 电影成人av| 久久午夜综合久久蜜桃| 老熟女久久久| 国产成人免费观看mmmm| 99久久精品国产国产毛片| 日韩视频在线欧美| 中文欧美无线码| 观看美女的网站| 久久久久久人妻| 99精国产麻豆久久婷婷| 国产亚洲最大av| 激情视频va一区二区三区| 日韩欧美精品免费久久| 老女人水多毛片| 国产女主播在线喷水免费视频网站| 各种免费的搞黄视频| 国产日韩一区二区三区精品不卡| 十八禁高潮呻吟视频| 宅男免费午夜| 91aial.com中文字幕在线观看| av国产精品久久久久影院| www.自偷自拍.com| 国产有黄有色有爽视频| 国产日韩一区二区三区精品不卡| 免费黄网站久久成人精品| 国产熟女欧美一区二区| av天堂久久9| 晚上一个人看的免费电影| 90打野战视频偷拍视频| 国产精品 欧美亚洲| 国产精品一二三区在线看| 国产探花极品一区二区| 一本—道久久a久久精品蜜桃钙片| 2021少妇久久久久久久久久久| 久久精品熟女亚洲av麻豆精品| 丝袜人妻中文字幕| 久久午夜福利片| 国产一区有黄有色的免费视频| 午夜福利网站1000一区二区三区| 一区二区三区四区激情视频| 999精品在线视频| 人妻系列 视频| 欧美日韩精品成人综合77777| 亚洲欧美成人综合另类久久久| 又大又黄又爽视频免费| 丰满乱子伦码专区| 日日爽夜夜爽网站| 亚洲综合精品二区| 丰满饥渴人妻一区二区三| 熟妇人妻不卡中文字幕| 在线天堂最新版资源| 欧美bdsm另类| 久久久久久久精品精品| 好男人视频免费观看在线| 亚洲国产成人一精品久久久| 1024香蕉在线观看| 你懂的网址亚洲精品在线观看| 一区在线观看完整版| 国产视频首页在线观看| 精品久久久精品久久久| 麻豆乱淫一区二区| 国产伦理片在线播放av一区| 成人毛片a级毛片在线播放| 免费播放大片免费观看视频在线观看| 91精品伊人久久大香线蕉| 国产精品熟女久久久久浪| 中文字幕亚洲精品专区| 在线精品无人区一区二区三| 精品福利永久在线观看| 国产国语露脸激情在线看| 一区二区三区四区激情视频| 在线免费观看不下载黄p国产| 美女午夜性视频免费| 永久网站在线| 日韩制服骚丝袜av| 免费观看无遮挡的男女| 欧美亚洲日本最大视频资源| 尾随美女入室| 国产黄频视频在线观看| 亚洲综合色网址| 中文字幕精品免费在线观看视频| 男女午夜视频在线观看| 看十八女毛片水多多多| 久久精品国产自在天天线| 亚洲av电影在线观看一区二区三区| 国产免费一区二区三区四区乱码| 一级黄片播放器| 国产一区有黄有色的免费视频| 性高湖久久久久久久久免费观看| 麻豆av在线久日| 日韩大片免费观看网站| av在线播放精品| 黄网站色视频无遮挡免费观看| 人妻 亚洲 视频| 国产免费福利视频在线观看| 亚洲人成电影观看| 黄色怎么调成土黄色| 一区二区三区激情视频| 蜜桃在线观看..| 一级,二级,三级黄色视频| 亚洲欧洲日产国产| 亚洲一区中文字幕在线| 波多野结衣av一区二区av| 亚洲av电影在线进入| 国产精品秋霞免费鲁丝片| 啦啦啦中文免费视频观看日本| 九色亚洲精品在线播放| 男人添女人高潮全过程视频| 亚洲国产精品一区二区三区在线| 黄色毛片三级朝国网站| 热99久久久久精品小说推荐| 日韩一卡2卡3卡4卡2021年| 欧美精品亚洲一区二区| 丝袜美腿诱惑在线| 看免费av毛片| 国产不卡av网站在线观看| 亚洲欧美成人综合另类久久久| a 毛片基地| 精品酒店卫生间| av网站在线播放免费| 欧美 日韩 精品 国产| 日本-黄色视频高清免费观看| 9色porny在线观看| 男女无遮挡免费网站观看| 中文欧美无线码| 欧美人与性动交α欧美软件| 女人久久www免费人成看片| 日本猛色少妇xxxxx猛交久久| 久久久久久久国产电影| 免费人妻精品一区二区三区视频| 中文乱码字字幕精品一区二区三区| 国产精品二区激情视频| 黑人猛操日本美女一级片| 99久久人妻综合| 这个男人来自地球电影免费观看 | 少妇人妻精品综合一区二区| 欧美xxⅹ黑人| 久久精品国产综合久久久| 亚洲欧美成人精品一区二区| 日韩三级伦理在线观看| www.自偷自拍.com| 成人黄色视频免费在线看| 在现免费观看毛片| 99国产综合亚洲精品| 丰满迷人的少妇在线观看| 亚洲少妇的诱惑av| 亚洲熟女精品中文字幕| a 毛片基地| 色哟哟·www| 激情视频va一区二区三区| 美女脱内裤让男人舔精品视频| 久久午夜综合久久蜜桃| 国产日韩欧美在线精品| 日韩一卡2卡3卡4卡2021年| 丰满饥渴人妻一区二区三| 国产成人精品久久久久久| 成人漫画全彩无遮挡| 丝袜喷水一区| 日本wwww免费看| 9热在线视频观看99| 成人亚洲精品一区在线观看| 亚洲色图综合在线观看| 国产毛片在线视频| 欧美另类一区| 男人添女人高潮全过程视频| 少妇的丰满在线观看| 欧美成人精品欧美一级黄| 久久久久久久久久久久大奶| 99久国产av精品国产电影| 美女高潮到喷水免费观看| 亚洲欧洲日产国产| 精品国产一区二区三区四区第35| 18禁裸乳无遮挡动漫免费视频| 蜜桃国产av成人99| 大香蕉久久成人网| 久久人人爽人人片av| 欧美xxⅹ黑人| 亚洲欧美清纯卡通| 精品国产一区二区三区久久久樱花| 校园人妻丝袜中文字幕| 国产成人午夜福利电影在线观看| 午夜久久久在线观看| 新久久久久国产一级毛片| 国产精品无大码| 美女国产高潮福利片在线看| 寂寞人妻少妇视频99o| 中文字幕av电影在线播放| 午夜免费男女啪啪视频观看| 色哟哟·www| 日本91视频免费播放| 一区在线观看完整版| 侵犯人妻中文字幕一二三四区| 欧美激情高清一区二区三区 | 熟妇人妻不卡中文字幕| 青草久久国产| 亚洲精品成人av观看孕妇| 久久精品国产亚洲av高清一级| 国产成人免费无遮挡视频| 尾随美女入室| 成人漫画全彩无遮挡| 日韩电影二区| 亚洲av福利一区| 日日啪夜夜爽| 国产成人91sexporn| 久久久久久人人人人人| 99香蕉大伊视频| 一边摸一边做爽爽视频免费| 亚洲一区二区三区欧美精品| 最近的中文字幕免费完整| 亚洲国产精品一区二区三区在线| 亚洲四区av| 成人亚洲欧美一区二区av| 中文字幕av电影在线播放| 久久av网站| 国产精品三级大全| 午夜免费观看性视频| 国产精品人妻久久久影院| 老熟女久久久| 日韩一卡2卡3卡4卡2021年| 女人被躁到高潮嗷嗷叫费观| 欧美日韩精品成人综合77777| 丰满饥渴人妻一区二区三| 国产精品 国内视频| 免费日韩欧美在线观看| 大话2 男鬼变身卡| 中文乱码字字幕精品一区二区三区| 国产熟女欧美一区二区| 一本大道久久a久久精品| 国产精品一区二区在线观看99| 看免费av毛片| www.av在线官网国产| 美女xxoo啪啪120秒动态图| 国产av一区二区精品久久| 欧美精品国产亚洲| 久久精品熟女亚洲av麻豆精品| 日韩三级伦理在线观看| 欧美在线黄色| 午夜久久久在线观看| 久久国产精品男人的天堂亚洲| 99热网站在线观看| 国产野战对白在线观看| 免费不卡的大黄色大毛片视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产av新网站| 久久久久久久大尺度免费视频| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久精品电影小说| 天堂中文最新版在线下载| 亚洲综合精品二区| 午夜福利在线观看免费完整高清在| 一级毛片我不卡| 国产男女超爽视频在线观看| 咕卡用的链子| 久久人妻熟女aⅴ| 亚洲欧美清纯卡通| 亚洲国产精品一区二区三区在线| 亚洲一区中文字幕在线| 亚洲精品aⅴ在线观看| 欧美 亚洲 国产 日韩一| 你懂的网址亚洲精品在线观看| av国产久精品久网站免费入址| 精品酒店卫生间| 久久97久久精品| 老汉色av国产亚洲站长工具| 久久 成人 亚洲| 日韩,欧美,国产一区二区三区| 国产成人av激情在线播放| 1024香蕉在线观看| 亚洲情色 制服丝袜| 欧美老熟妇乱子伦牲交| 成人漫画全彩无遮挡| 亚洲第一青青草原| 天美传媒精品一区二区| 日韩,欧美,国产一区二区三区| 黄色怎么调成土黄色| 精品少妇黑人巨大在线播放| xxxhd国产人妻xxx| 久久久久久久久久久免费av| 午夜福利视频精品| 人体艺术视频欧美日本| 婷婷色综合www| 黄频高清免费视频| 欧美日韩精品网址| a级毛片在线看网站| 一级a爱视频在线免费观看| 亚洲激情五月婷婷啪啪| 亚洲精品在线美女| 中文字幕人妻熟女乱码| 老熟女久久久| 成年女人毛片免费观看观看9 | 人妻少妇偷人精品九色| 18在线观看网站| 亚洲精品第二区| 汤姆久久久久久久影院中文字幕| 韩国高清视频一区二区三区| 午夜av观看不卡| 美女视频免费永久观看网站| 又大又黄又爽视频免费| 国产野战对白在线观看| 久久久国产一区二区| 美女午夜性视频免费| 亚洲精品国产一区二区精华液| 国产淫语在线视频| 中国国产av一级| 欧美人与善性xxx| xxxhd国产人妻xxx| 夫妻性生交免费视频一级片| av国产久精品久网站免费入址| 精品国产露脸久久av麻豆| 一本大道久久a久久精品| 国产免费福利视频在线观看| 狂野欧美激情性bbbbbb| 国产一区二区三区av在线| 黄色 视频免费看| 最黄视频免费看| 日韩制服丝袜自拍偷拍| 成年女人毛片免费观看观看9 | 人妻系列 视频| 91aial.com中文字幕在线观看| 亚洲国产欧美日韩在线播放| 极品人妻少妇av视频| 波多野结衣av一区二区av| 在线观看三级黄色| 亚洲视频免费观看视频| 日韩成人av中文字幕在线观看| 国产在线免费精品| 夜夜骑夜夜射夜夜干| 丝袜美腿诱惑在线| 黄色视频在线播放观看不卡| 欧美精品一区二区免费开放| www.自偷自拍.com| xxx大片免费视频| 巨乳人妻的诱惑在线观看| 2018国产大陆天天弄谢| 一区二区三区四区激情视频| 色网站视频免费| 在现免费观看毛片| 美女大奶头黄色视频| 欧美精品亚洲一区二区| 亚洲熟女精品中文字幕| 女人高潮潮喷娇喘18禁视频| 国产亚洲欧美精品永久| kizo精华| 国产1区2区3区精品| 熟女电影av网| 国产精品免费视频内射| 如日韩欧美国产精品一区二区三区| 欧美日韩亚洲高清精品| 亚洲国产av影院在线观看| 国产亚洲av片在线观看秒播厂| 精品久久蜜臀av无| 1024香蕉在线观看| 日韩成人av中文字幕在线观看| 伦理电影大哥的女人| 亚洲欧美清纯卡通| 午夜日本视频在线| 久热久热在线精品观看| 最近最新中文字幕免费大全7| 男人添女人高潮全过程视频| 国产精品女同一区二区软件| 夫妻性生交免费视频一级片| 国产精品久久久久久精品电影小说| freevideosex欧美| 久久久久久人人人人人| 中文字幕色久视频| 熟女av电影| 国产免费视频播放在线视频| 久久鲁丝午夜福利片| 亚洲,一卡二卡三卡| 国产精品av久久久久免费| 不卡av一区二区三区| 久久免费观看电影| 黄色 视频免费看| 亚洲av男天堂| 男人添女人高潮全过程视频| 我的亚洲天堂| 大陆偷拍与自拍| 国产 精品1| 成人黄色视频免费在线看| 十八禁高潮呻吟视频| 免费久久久久久久精品成人欧美视频| 国产毛片在线视频| 一级毛片我不卡| 日韩人妻精品一区2区三区| 亚洲av综合色区一区| 国产精品秋霞免费鲁丝片| 精品一区二区三卡| 国产综合精华液| 午夜福利在线免费观看网站| 午夜激情av网站| 中文欧美无线码| 伊人久久国产一区二区| 欧美少妇被猛烈插入视频| 老汉色av国产亚洲站长工具| 2021少妇久久久久久久久久久| 婷婷色麻豆天堂久久| 狠狠精品人妻久久久久久综合| 免费不卡的大黄色大毛片视频在线观看| av视频免费观看在线观看| 精品少妇久久久久久888优播| 叶爱在线成人免费视频播放| 亚洲精品第二区| 成人国产麻豆网| 美女中出高潮动态图| 我的亚洲天堂| 久久人妻熟女aⅴ| 国产精品一二三区在线看| 国产一区二区在线观看av| 午夜福利乱码中文字幕| 黄色毛片三级朝国网站| 永久免费av网站大全| 性色av一级| 精品一区二区免费观看| 国产av一区二区精品久久| 国产1区2区3区精品| 成人毛片60女人毛片免费| 欧美亚洲日本最大视频资源| 国产精品人妻久久久影院| 午夜福利,免费看| 国语对白做爰xxxⅹ性视频网站| 1024视频免费在线观看| 国产在线一区二区三区精| 免费av中文字幕在线| 在线精品无人区一区二区三| 99久国产av精品国产电影| 成人手机av| 十八禁高潮呻吟视频| 91久久精品国产一区二区三区| av有码第一页| 亚洲欧美精品综合一区二区三区 | 妹子高潮喷水视频| 日韩三级伦理在线观看| 亚洲精品美女久久久久99蜜臀 | 精品卡一卡二卡四卡免费| 一级,二级,三级黄色视频| 日韩制服丝袜自拍偷拍| 国产片内射在线| 免费高清在线观看日韩| 一级黄片播放器| 日本vs欧美在线观看视频| 午夜福利一区二区在线看| 亚洲成人手机| 最新的欧美精品一区二区| 女人高潮潮喷娇喘18禁视频| 国产日韩一区二区三区精品不卡| 在线天堂最新版资源| 久久精品久久久久久久性| 亚洲精品乱久久久久久| 亚洲精品美女久久久久99蜜臀 | 一级毛片我不卡| 一区二区三区精品91| 久久久a久久爽久久v久久| 黄片播放在线免费| 亚洲欧美一区二区三区黑人 | 国产一区亚洲一区在线观看| 91精品三级在线观看| www日本在线高清视频| 波多野结衣一区麻豆| 日本欧美视频一区| 日日爽夜夜爽网站| 男女啪啪激烈高潮av片| 亚洲少妇的诱惑av| 黄片播放在线免费| 高清黄色对白视频在线免费看| 亚洲国产成人一精品久久久| 国产精品欧美亚洲77777| 久久久久久久久久久免费av| 一级黄片播放器| 18禁裸乳无遮挡动漫免费视频| 99国产精品免费福利视频| 老司机影院成人| 欧美日韩一级在线毛片| 成人午夜精彩视频在线观看| 丝袜美足系列| 大香蕉久久成人网| 国产xxxxx性猛交| 日日摸夜夜添夜夜爱| 老鸭窝网址在线观看| 老熟女久久久| 亚洲情色 制服丝袜| 在线看a的网站| 久久人人97超碰香蕉20202| 亚洲欧洲精品一区二区精品久久久 | 亚洲视频免费观看视频| 日韩免费高清中文字幕av| 两性夫妻黄色片| 美女中出高潮动态图| 精品少妇内射三级| 超色免费av| 日韩视频在线欧美| 少妇猛男粗大的猛烈进出视频| 丝袜人妻中文字幕| 亚洲欧美色中文字幕在线| 啦啦啦视频在线资源免费观看| 亚洲一区二区三区欧美精品| 波多野结衣av一区二区av| 久久97久久精品| 亚洲精品自拍成人| 天天躁狠狠躁夜夜躁狠狠躁| 日日爽夜夜爽网站| 最新的欧美精品一区二区| 亚洲欧美色中文字幕在线| 国产1区2区3区精品| 欧美亚洲 丝袜 人妻 在线| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕av电影在线播放| 18在线观看网站| 国产免费又黄又爽又色| 99久久综合免费| 国产精品无大码| 久热这里只有精品99| 男男h啪啪无遮挡| 亚洲精品,欧美精品| 美女中出高潮动态图| 午夜91福利影院| 免费观看a级毛片全部| 日韩欧美一区视频在线观看| 久久午夜福利片| 一区在线观看完整版| 热99久久久久精品小说推荐| 国产精品熟女久久久久浪| 精品国产一区二区久久| 国产精品二区激情视频| 亚洲欧美一区二区三区国产| 成人国语在线视频| 女性生殖器流出的白浆|