• <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
    探秘“大徑流”
    攻克“大徑流”
    少妇熟女aⅴ在线视频| 亚洲黑人精品在线| 亚洲第一电影网av| 黄网站色视频无遮挡免费观看| 国产精品精品国产色婷婷| 九色国产91popny在线| 色综合亚洲欧美另类图片| 国产单亲对白刺激| 好男人电影高清在线观看| 一边摸一边抽搐一进一小说| 国产免费男女视频| 色综合站精品国产| 免费在线观看日本一区| 99热只有精品国产| 正在播放国产对白刺激| 亚洲欧美日韩无卡精品| 波多野结衣一区麻豆| 一进一出好大好爽视频| 久久人妻av系列| 久久精品国产清高在天天线| 大码成人一级视频| 国产成人系列免费观看| 操美女的视频在线观看| 女同久久另类99精品国产91| 日本免费a在线| 黄色毛片三级朝国网站| 亚洲熟妇中文字幕五十中出| 亚洲精品国产色婷婷电影| 日本 av在线| √禁漫天堂资源中文www| 亚洲精品久久成人aⅴ小说| 欧美成人免费av一区二区三区| 乱人伦中国视频| aaaaa片日本免费| 99久久久亚洲精品蜜臀av| 欧美人与性动交α欧美精品济南到| 国产成人欧美| 国产国语露脸激情在线看| 一进一出抽搐gif免费好疼| 热99re8久久精品国产| 悠悠久久av| 国产精品美女特级片免费视频播放器 | 激情在线观看视频在线高清| 成在线人永久免费视频| 中文字幕av电影在线播放| 老司机福利观看| 亚洲免费av在线视频| 中国美女看黄片| 两性夫妻黄色片| bbb黄色大片| 久久久国产成人免费| 99热只有精品国产| 国产男靠女视频免费网站| 激情在线观看视频在线高清| 99国产精品99久久久久| 亚洲情色 制服丝袜| 亚洲国产精品久久男人天堂| 亚洲欧美激情综合另类| 亚洲人成电影免费在线| 国产欧美日韩一区二区精品| 国产av又大| 别揉我奶头~嗯~啊~动态视频| 日韩视频一区二区在线观看| 露出奶头的视频| 日本免费a在线| 亚洲精品中文字幕一二三四区| 亚洲七黄色美女视频| 国产高清有码在线观看视频 | 日韩欧美一区二区三区在线观看| 国产免费av片在线观看野外av| 国产精品国产高清国产av| 一边摸一边做爽爽视频免费| 欧美日本亚洲视频在线播放| 露出奶头的视频| 两性夫妻黄色片| 中文字幕色久视频| 精品卡一卡二卡四卡免费| 亚洲人成网站在线播放欧美日韩| av在线播放免费不卡| 免费高清在线观看日韩| 精品欧美国产一区二区三| 欧美成人免费av一区二区三区| 韩国av一区二区三区四区| 999久久久精品免费观看国产| 精品一品国产午夜福利视频| 亚洲一码二码三码区别大吗| 一级片免费观看大全| 日本vs欧美在线观看视频| 黄色a级毛片大全视频| 手机成人av网站| 最好的美女福利视频网| 美女高潮到喷水免费观看| 黄片大片在线免费观看| 正在播放国产对白刺激| 侵犯人妻中文字幕一二三四区| 欧美一区二区精品小视频在线| 欧美av亚洲av综合av国产av| 国产高清有码在线观看视频 | 国产精品亚洲一级av第二区| 日本 av在线| 男女之事视频高清在线观看| 成熟少妇高潮喷水视频| 91九色精品人成在线观看| 91字幕亚洲| 99国产精品免费福利视频| 色播亚洲综合网| 精品国产亚洲在线| 一进一出好大好爽视频| 老汉色av国产亚洲站长工具| 波多野结衣一区麻豆| 日韩欧美国产在线观看| 国产99白浆流出| 日本免费一区二区三区高清不卡 | 大香蕉久久成人网| 亚洲精品国产区一区二| 婷婷精品国产亚洲av在线| 亚洲片人在线观看| 成人av一区二区三区在线看| 一区二区三区高清视频在线| 日韩成人在线观看一区二区三区| 精品国产美女av久久久久小说| 一进一出好大好爽视频| 色综合站精品国产| 啪啪无遮挡十八禁网站| 搡老岳熟女国产| 又黄又爽又免费观看的视频| 99精品在免费线老司机午夜| 动漫黄色视频在线观看| 国产一卡二卡三卡精品| 婷婷六月久久综合丁香| 国产精品亚洲一级av第二区| 欧美中文日本在线观看视频| 嫁个100分男人电影在线观看| 国产成人一区二区三区免费视频网站| 久久久久久亚洲精品国产蜜桃av| 久久久国产成人精品二区| 成熟少妇高潮喷水视频| 久久青草综合色| 一区二区三区高清视频在线| 亚洲九九香蕉| 亚洲精华国产精华精| 人成视频在线观看免费观看| 亚洲成国产人片在线观看| 丝袜在线中文字幕| 亚洲久久久国产精品| 国产成人精品久久二区二区91| 欧美一级毛片孕妇| 欧美绝顶高潮抽搐喷水| 国产精品香港三级国产av潘金莲| 欧美午夜高清在线| 人成视频在线观看免费观看| 亚洲国产精品久久男人天堂| 在线观看www视频免费| 亚洲精品久久成人aⅴ小说| av视频在线观看入口| 欧美一级毛片孕妇| 在线十欧美十亚洲十日本专区| 免费少妇av软件| 欧美日韩亚洲综合一区二区三区_| 三级毛片av免费| 99国产精品99久久久久| 搡老岳熟女国产| 亚洲av第一区精品v没综合| 后天国语完整版免费观看| 国产av一区在线观看免费| 国产亚洲精品av在线| 12—13女人毛片做爰片一| 久久精品影院6| 制服诱惑二区| 怎么达到女性高潮| 免费在线观看影片大全网站| 天堂动漫精品| 极品人妻少妇av视频| 咕卡用的链子| 村上凉子中文字幕在线| 国产亚洲精品久久久久5区| 人成视频在线观看免费观看| 亚洲欧美激情综合另类| 在线播放国产精品三级| 午夜日韩欧美国产| 女人高潮潮喷娇喘18禁视频| 男女下面进入的视频免费午夜 | 脱女人内裤的视频| 每晚都被弄得嗷嗷叫到高潮| 丝袜美腿诱惑在线| 亚洲一区高清亚洲精品| 亚洲最大成人中文| 精品国产美女av久久久久小说| 丝袜人妻中文字幕| 啦啦啦韩国在线观看视频| 99久久久亚洲精品蜜臀av| 久久婷婷人人爽人人干人人爱 | 麻豆国产av国片精品| 亚洲av电影不卡..在线观看| 激情视频va一区二区三区| 国产一区二区三区综合在线观看| 国产欧美日韩综合在线一区二区| 精品卡一卡二卡四卡免费| 亚洲成人精品中文字幕电影| 欧美激情高清一区二区三区| 亚洲国产精品久久男人天堂| 日本黄色视频三级网站网址| 国产精品综合久久久久久久免费 | 美女免费视频网站| 视频在线观看一区二区三区| 欧美另类亚洲清纯唯美| 激情视频va一区二区三区| 国产欧美日韩精品亚洲av| 可以在线观看的亚洲视频| 午夜日韩欧美国产| 亚洲国产精品999在线| 久久久久久久午夜电影| 精品一区二区三区四区五区乱码| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人精品中文字幕电影| a级毛片在线看网站| 最近最新中文字幕大全电影3 | 欧美黑人欧美精品刺激| av福利片在线| 亚洲一区二区三区不卡视频| 日韩欧美一区视频在线观看| a在线观看视频网站| 成人18禁高潮啪啪吃奶动态图| 精品一区二区三区视频在线观看免费| 日本精品一区二区三区蜜桃| 制服诱惑二区| 精品久久久精品久久久| 男女下面进入的视频免费午夜 | 久热爱精品视频在线9| 久久精品国产清高在天天线| 很黄的视频免费| 亚洲一区二区三区不卡视频| 午夜影院日韩av| 99热只有精品国产| 97超级碰碰碰精品色视频在线观看| 亚洲专区字幕在线| 日韩 欧美 亚洲 中文字幕| 十八禁人妻一区二区| 精品国内亚洲2022精品成人| 亚洲欧美日韩无卡精品| av在线天堂中文字幕| 亚洲成a人片在线一区二区| 国产精品乱码一区二三区的特点 | 国产麻豆69| 久久精品国产综合久久久| 欧美在线黄色| 无限看片的www在线观看| 99国产极品粉嫩在线观看| 一进一出好大好爽视频| 亚洲成av人片免费观看| 久久香蕉激情| 久久久久久久久免费视频了| 人人妻人人澡欧美一区二区 | 欧美色欧美亚洲另类二区 | 美女高潮喷水抽搐中文字幕| 国产亚洲欧美98| 久久 成人 亚洲| 亚洲欧美日韩无卡精品| 两性夫妻黄色片| 韩国av一区二区三区四区| 久久国产乱子伦精品免费另类| а√天堂www在线а√下载| 法律面前人人平等表现在哪些方面| 亚洲中文字幕日韩| 在线天堂中文资源库| 女人精品久久久久毛片| 少妇 在线观看| 国产黄a三级三级三级人| 99久久精品国产亚洲精品| 色老头精品视频在线观看| 国产免费av片在线观看野外av| 午夜精品在线福利| 亚洲国产日韩欧美精品在线观看 | 亚洲美女黄片视频| 久久精品国产亚洲av高清一级| 亚洲天堂国产精品一区在线| www.自偷自拍.com| 99riav亚洲国产免费| 一级毛片女人18水好多| 在线观看一区二区三区| 国产精品美女特级片免费视频播放器 | 极品人妻少妇av视频| 9热在线视频观看99| 日本精品一区二区三区蜜桃| 91麻豆av在线| 久热爱精品视频在线9| 午夜免费观看网址| 午夜老司机福利片| 在线观看66精品国产| 女性生殖器流出的白浆| 自拍欧美九色日韩亚洲蝌蚪91| 国产午夜福利久久久久久| 黄色视频,在线免费观看| 在线av久久热| 久久精品91蜜桃| 国产亚洲欧美在线一区二区| 欧美中文综合在线视频| 国产精品久久久久久亚洲av鲁大| 久久天躁狠狠躁夜夜2o2o| 亚洲精品av麻豆狂野| 久久狼人影院| 一级作爱视频免费观看| 又黄又粗又硬又大视频| 久久精品国产清高在天天线| 久久香蕉激情| 变态另类成人亚洲欧美熟女 | 午夜精品久久久久久毛片777| 看黄色毛片网站| 无遮挡黄片免费观看| 国产精华一区二区三区| 日韩精品中文字幕看吧| 亚洲av日韩精品久久久久久密| 国产精品亚洲av一区麻豆| 神马国产精品三级电影在线观看 | 色综合婷婷激情| 色综合亚洲欧美另类图片| 久久久精品国产亚洲av高清涩受| 一边摸一边做爽爽视频免费| 嫩草影院精品99| 久久热在线av| av在线播放免费不卡| 两性夫妻黄色片| 精品国产乱子伦一区二区三区| 中文字幕人成人乱码亚洲影| 可以在线观看的亚洲视频| 看片在线看免费视频| 成人三级黄色视频| 一区二区三区精品91| 热re99久久国产66热| 搡老岳熟女国产| 色精品久久人妻99蜜桃| 91麻豆av在线| 一本综合久久免费| 亚洲片人在线观看| 午夜免费观看网址| 日本一区二区免费在线视频| 亚洲午夜精品一区,二区,三区| 国产精品,欧美在线| 一进一出抽搐动态| 国产亚洲av高清不卡| 久久久久久大精品| 亚洲国产中文字幕在线视频| 精品电影一区二区在线| 久久精品91无色码中文字幕| 成人精品一区二区免费| 欧美成人午夜精品| 露出奶头的视频| 叶爱在线成人免费视频播放| 一边摸一边抽搐一进一小说| 亚洲精品av麻豆狂野| 日韩免费av在线播放| 俄罗斯特黄特色一大片| 欧美日韩亚洲国产一区二区在线观看| 999久久久精品免费观看国产| 久久中文字幕人妻熟女| 伦理电影免费视频| 久久香蕉精品热| 男女做爰动态图高潮gif福利片 | 日韩一卡2卡3卡4卡2021年| 无限看片的www在线观看| 18禁国产床啪视频网站| 久久草成人影院| 亚洲 欧美 日韩 在线 免费| а√天堂www在线а√下载| 婷婷丁香在线五月| 午夜视频精品福利| 国产精品 国内视频| 国产精品一区二区免费欧美| 伊人久久大香线蕉亚洲五| 男女下面插进去视频免费观看| 久久影院123| 丁香六月欧美| 在线视频色国产色| 久久午夜亚洲精品久久| 99国产精品一区二区三区| 黄色视频不卡| 色尼玛亚洲综合影院| 少妇 在线观看| 久久久久久免费高清国产稀缺| 久久久久久久久久久久大奶| 欧美最黄视频在线播放免费| 亚洲成人国产一区在线观看| 18禁黄网站禁片午夜丰满| 成人18禁高潮啪啪吃奶动态图| 久久久久亚洲av毛片大全| 日韩大码丰满熟妇| 日韩精品中文字幕看吧| 国内毛片毛片毛片毛片毛片| 亚洲色图综合在线观看| 在线观看免费视频网站a站| 亚洲一码二码三码区别大吗| 免费观看精品视频网站| 少妇被粗大的猛进出69影院| 免费在线观看亚洲国产| 看黄色毛片网站| 美女午夜性视频免费| 99精品久久久久人妻精品| 少妇 在线观看| 亚洲精品美女久久久久99蜜臀| 免费人成视频x8x8入口观看| 老司机福利观看| 国产99白浆流出| 亚洲七黄色美女视频| 精品卡一卡二卡四卡免费| 亚洲久久久国产精品| 日本在线视频免费播放| 国产不卡一卡二| 国产精品久久视频播放| 欧美乱色亚洲激情| 成人亚洲精品一区在线观看| 啦啦啦观看免费观看视频高清 | 久久午夜亚洲精品久久| av免费在线观看网站| 久久中文字幕人妻熟女| 久久精品国产综合久久久| 人人妻人人爽人人添夜夜欢视频| 久久久国产成人免费| 怎么达到女性高潮| 国产不卡一卡二| 国内精品久久久久久久电影| 日日摸夜夜添夜夜添小说| 久久中文看片网| 午夜成年电影在线免费观看| 天堂√8在线中文| 如日韩欧美国产精品一区二区三区| 超碰成人久久| 久久精品人人爽人人爽视色| 亚洲天堂国产精品一区在线| 国产精品国产高清国产av| 亚洲男人天堂网一区| 国内久久婷婷六月综合欲色啪| 国产单亲对白刺激| 99久久综合精品五月天人人| 亚洲国产精品久久男人天堂| 亚洲国产看品久久| 久久草成人影院| 亚洲熟女毛片儿| a级毛片在线看网站| 日韩精品中文字幕看吧| 极品教师在线免费播放| 淫妇啪啪啪对白视频| 宅男免费午夜| 国产精品永久免费网站| 国产又爽黄色视频| 搡老妇女老女人老熟妇| 一本综合久久免费| 国产精品久久久久久精品电影 | 免费观看人在逋| 成人18禁在线播放| 一本大道久久a久久精品| 久久中文看片网| 国产在线精品亚洲第一网站| 99国产精品免费福利视频| 精品久久久久久久久久免费视频| 午夜福利,免费看| 九色国产91popny在线| 久久天堂一区二区三区四区| 亚洲天堂国产精品一区在线| 欧美黄色片欧美黄色片| 久久午夜综合久久蜜桃| 午夜福利,免费看| 成人欧美大片| 露出奶头的视频| 在线国产一区二区在线| 精品久久久久久久毛片微露脸| 国产成人精品久久二区二区免费| 成人亚洲精品一区在线观看| 精品卡一卡二卡四卡免费| 淫秽高清视频在线观看| 亚洲av成人一区二区三| 亚洲av日韩精品久久久久久密| 无限看片的www在线观看| 少妇裸体淫交视频免费看高清 | 久久草成人影院| 国产又爽黄色视频| 欧美另类亚洲清纯唯美| 美女大奶头视频| 性欧美人与动物交配| 国产野战对白在线观看| 欧美丝袜亚洲另类 | 国产av一区在线观看免费| 国产欧美日韩一区二区三| 中文字幕最新亚洲高清| 在线观看免费视频网站a站| 欧美中文综合在线视频| 他把我摸到了高潮在线观看| 两个人看的免费小视频| 韩国精品一区二区三区| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线观看二区| 午夜精品在线福利| 12—13女人毛片做爰片一| 国产成人影院久久av| 亚洲成人免费电影在线观看| 青草久久国产| 老司机靠b影院| 99精品久久久久人妻精品| 欧美亚洲日本最大视频资源| 99香蕉大伊视频| 国产精品99久久99久久久不卡| 1024视频免费在线观看| 亚洲精品在线美女| 制服诱惑二区| 国产一级毛片七仙女欲春2 | 免费少妇av软件| 欧美另类亚洲清纯唯美| 男女做爰动态图高潮gif福利片 | 久久国产精品影院| 国产色视频综合| 亚洲国产高清在线一区二区三 | 波多野结衣一区麻豆| 亚洲电影在线观看av| 精品国产一区二区久久| 精品电影一区二区在线| 两人在一起打扑克的视频| 伊人久久大香线蕉亚洲五| 很黄的视频免费| 亚洲欧美日韩无卡精品| av福利片在线| 一a级毛片在线观看| 国产片内射在线| 亚洲成人国产一区在线观看| 亚洲国产日韩欧美精品在线观看 | 精品国产乱码久久久久久男人| 无限看片的www在线观看| 免费观看人在逋| 亚洲国产高清在线一区二区三 | 老汉色av国产亚洲站长工具| 熟妇人妻久久中文字幕3abv| 免费搜索国产男女视频| 亚洲 欧美 日韩 在线 免费| 久久久国产成人免费| 悠悠久久av| 免费不卡黄色视频| 国产一区在线观看成人免费| 国产精品 欧美亚洲| 在线免费观看的www视频| 亚洲欧美日韩高清在线视频| 久久久国产精品麻豆| 中文字幕精品免费在线观看视频| 一级片免费观看大全| 午夜两性在线视频| 禁无遮挡网站| 亚洲专区中文字幕在线| 桃色一区二区三区在线观看| 男人舔女人下体高潮全视频| 91大片在线观看| 国产一区二区三区综合在线观看| 中文字幕高清在线视频| 亚洲欧美一区二区三区黑人| 亚洲精品一卡2卡三卡4卡5卡| 1024视频免费在线观看| 色哟哟哟哟哟哟| 国产成人av教育| 成人三级黄色视频| 丰满人妻熟妇乱又伦精品不卡| 精品乱码久久久久久99久播| 亚洲精品av麻豆狂野| 久久中文字幕人妻熟女| 亚洲av五月六月丁香网| 亚洲欧美一区二区三区黑人| 身体一侧抽搐| 夜夜夜夜夜久久久久| 啦啦啦免费观看视频1| 级片在线观看| 天堂√8在线中文| 亚洲中文日韩欧美视频| 极品教师在线免费播放| 亚洲一区中文字幕在线| 日韩高清综合在线| 精品国产亚洲在线| 人人妻人人澡人人看| 在线视频色国产色| 国内精品久久久久精免费| 国产一区二区激情短视频| 亚洲专区字幕在线| 一级,二级,三级黄色视频| 亚洲精品一区av在线观看| 精品久久久久久成人av| 黄片小视频在线播放| 脱女人内裤的视频| 久久久精品欧美日韩精品| 亚洲最大成人中文| 欧美日韩亚洲综合一区二区三区_| 国产亚洲精品综合一区在线观看 | 色在线成人网| 熟女少妇亚洲综合色aaa.| 国产精品免费一区二区三区在线| 日日爽夜夜爽网站| 国产伦人伦偷精品视频| 精品少妇一区二区三区视频日本电影| 国产精品亚洲av一区麻豆| 一进一出好大好爽视频| 真人做人爱边吃奶动态| av在线天堂中文字幕| 欧美亚洲日本最大视频资源| av片东京热男人的天堂| 精品一品国产午夜福利视频| 亚洲欧美日韩高清在线视频| 国产精品av久久久久免费| 在线观看66精品国产| 99香蕉大伊视频| 亚洲伊人色综图| 欧美另类亚洲清纯唯美| 十八禁网站免费在线| 亚洲一区二区三区不卡视频| 午夜福利在线观看吧| 伊人久久大香线蕉亚洲五| 校园春色视频在线观看| 岛国视频午夜一区免费看| 女人爽到高潮嗷嗷叫在线视频| cao死你这个sao货| 女人被躁到高潮嗷嗷叫费观| 久久久久国产一级毛片高清牌| 国产精品,欧美在线|