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

    淮河洪水概率預報方法研究及應用

    2018-05-18 05:43:07蔣曉蕾
    治淮 2018年4期
    關鍵詞:王家壩貝葉斯不確定性

    王 凱 蔣曉蕾

    (淮河水利委員會水文局(信息中心) 蚌埠 233001 河海大學 南京 210098)

    1 研究背景

    水文預報是一種重要的防洪非工程措施,直接為防汛抗旱和水資源管理服務。目前廣泛使用的水文預報模型大多是確定性的,模型以確定預報值的形式輸出給用戶。實際上,由于水文過程的影響因素復雜多變,加之人類認識水平的有限,使得水文預報過程中不可避免地存在著諸多的不確定性。目前,水文預報不確定性量化已逐漸成為了研究熱點,概率預報作為其重要的表現(xiàn)形式,已成為了水文預報的發(fā)展趨勢。

    當前國內(nèi)外關于洪水概率預報的研究主要包括兩類途徑:第一類是全要素耦合途徑,分別量化降雨—徑流過程各個環(huán)節(jié)的主要不確定性,如降雨輸入不確定性、模型結構不確定性、模型參數(shù)不確定性等,并進行耦合,實現(xiàn)概率預報。如Kavetski等采用“潛在變量”雨深乘子(stormdepthmultiplier)反映降雨輸入的不確定性,并將模型的敏感性參數(shù)隨機化,應用馬爾可夫鏈蒙特卡洛(Markov Chain Monte Carlo,MCMC)采樣方法求解流量后驗分布,提出貝葉斯總誤差分析(Bayesian total error analysis,BATEA)方法。在國內(nèi),李明亮等基于層次貝葉斯模型(Bayesian hierarchical model),構建聯(lián)合概率密度函數(shù)以考慮模型參數(shù)和降雨輸入不確定性,并采用MCMC方法進行求解。梁忠民等借用抽站法原理推求降雨量的條件概率分布,進而實現(xiàn)考慮輸入不確定性的洪水概率預報。全要素耦合途徑能夠溯源預報不確定性,但計算相對耗時,無法滿足實時預報的需要。

    第二類是總誤差分析途徑,即從確定性預報結果入手,直接對預報不確定性進行量化分析,推求預報量的分布函數(shù),實現(xiàn)概率預報。Krzysztofowicz等提出的貝葉斯預報系統(tǒng)(Bayesian forecastingsystem,BFS)最具代表性,其中水文不確定性處理器(hydrologic uncertainty processor,HUP)在正態(tài)空間中對似然函數(shù)進行了線性假定,為推求預報量后驗分布的解析表達提供了可能。王善序詳細介紹了BFS的理論,并指出它可以綜合考慮預報過程的不確定性,不限定預報模型的內(nèi)部結構。近年來,很多研究表明,不同流量量級的預報不確定性存在差異。王晶晶等在原始空間中開展預報誤差規(guī)律研究,發(fā)現(xiàn)不同流量量級的誤差服從不同的分布函數(shù),為此采用極小熵確定各分布線型,采用極大熵進行參數(shù)估計,進而降低了福建池潭水庫管理的風險。Van Steenbergen等針對不同預見期、不同流量量級預報誤差統(tǒng)計規(guī)律的差異,采用頻率學方法,構建了三維誤差矩陣,量化預報不確定性。

    本文在第二類研究途徑框架下,結合研究區(qū)域(淮河王家壩斷面)預報誤差的統(tǒng)計特征,假定誤差均值隨流量量級呈分段線性變化,在貝葉斯預報系統(tǒng)(BFS)的基礎上,對其水文不確定性處理器(HUP)進行改進,提出改進的貝葉斯概率預報系統(tǒng)(PCA-HUP),定量評估單個水文模型預報結果的可靠度,實現(xiàn)概率預報。

    2 模型方法

    考慮誤差異分布的洪水概率預報方法不涉及預報的中間環(huán)節(jié),只對預報結果進行分析,定量評估預報結果的不確定性,并在此基礎上實現(xiàn)洪水概率預報。在貝葉斯預報系統(tǒng)(BFS)的基礎上,對其水文不確定性處理器(HUP)進行改進,提出改進的貝葉斯概率預報系統(tǒng)(PCA-HUP),定量評估單個水文模型預報結果的可靠度,實現(xiàn)概率預報。

    2.1 HUP基本原理

    水文不確定性處理器(Hydrologic UncertaintyProcessor,HUP)是貝葉斯概率預報系統(tǒng)(BFS)的主要組成部分,用以分析除降雨之外的其他所有不確定性。其特點是,不需要直接處理預報模型的結構與參數(shù),而是從預報結果入手,分析其與實測水文過程的誤差,再利用貝葉斯公式估計預報變量的概率分布,從而實現(xiàn)水文模型預報結果的不確定性分析及概率預報。其工作流程如圖1所示。

    2.2 改進的概率洪水預報PCA-HUP模型

    為推求預報量后驗分布的解析解,傳統(tǒng)HUP模型結合亞高斯模型,在正態(tài)空間中對先驗分布式和似然函數(shù)式進行線性假設,并采用最小二乘法對相關參數(shù)進行估計。

    然而,由于似然函數(shù)式的自變量之間存在明確的線性關系,必然導致回歸方程的多重共線性問題。若采用傳統(tǒng)最小二乘法進行參數(shù)估計,會使得估計的回歸系數(shù)不唯一,也使得回歸方程不穩(wěn)定(原始數(shù)據(jù)的極小變化可造成參數(shù)估計值和標準差的明顯變化)。因此,項目研究結合主成分分析技術(Principal components analysis,PCA),對傳統(tǒng) HUP模型進行改進,提出PCA-HUP模型。

    主成分回歸的基本思想:對原始回歸變量進行主成分分析,將線性相關的自變量,轉化為線性無關的新的綜合變量,采用新的綜合變量建立模型回歸方程。

    (1)主成分分析

    設 X=(X1,…Xp)T是 P 維隨機向量,均值 E(X)=μ,協(xié)方差陣 D(X)=Σ。

    考慮它的線性變換:

    用矩陣表示為:

    圖1 水文不確定性處理器工作流程示意圖

    由式(2)可以將 P 個 X1,X2,…Xp轉化為 P 個新變量,Z1,Z2,…Zp若新變量 Z1,Z2,…Zp滿足下列條件:

    1)Zi和 Zj相互獨立,i≠j,i,j=1,2,…,p;

    2)Var(Z1)≥Var(Z2)≥…≥Var(Zp);

    則新變量 Z1,Z2,…Zp為 X1,X2,…Xp的 P 個主成分,且Z1,Z2,…Zp線性無關。

    (2)主成分回歸

    實際問題中不同的變量經(jīng)常具有不同的量綱,變量的量綱不同會使分析結果不合理,將變量進行標準化處理可避免這種不合理的影響。

    對數(shù)據(jù)進行標準化首先要得出樣本標準差和樣本均值,記sj為xj樣本標準差,即是 xj的樣本均值,即原始數(shù)據(jù)的標準化變換為:

    標準化后的數(shù)據(jù)矩陣為:

    標準化后,X的相關系數(shù)陣也就是X的協(xié)方差陣(半正定矩陣):

    其中:

    采用Lagrange乘子法求解,可以求得:

    其中:λ1≥λ2≥…≥λp≥0 為 R 的特征值,a1,a2,…ap是相對應的單位正交特征向量,ap=(a1pa2p,…app)T。

    主成分回歸可以得到p個主成分,這p個主成分之間互相獨立,且方差呈遞減趨勢,所包含的自變量的信息也是遞減的。即主成分對因變量的貢獻率是遞減的,第i個主成分Zi的貢獻率可以用來表示。

    在實際問題的分析時,由于主成分的貢獻率是遞減的,后面的主成分貢獻率有時會非常小,所以一般不選取P個主成分,而是根據(jù)累計貢獻率來確定主成分個數(shù),即:前m個主成分的累計貢獻率達到0.85時,選取前m個主成分進行回歸。則原始回歸問題轉化為以下回歸問題:

    其中:E(εt)=0,Var(εt)=σ2,Cov(εi,εi)=0,(i≠j)

    回歸模型的矩陣形式為:

    采用最小二乘法估計參數(shù)矩陣B,根據(jù)式(1)和式(3)可以估計因變量矩陣Y與自變量矩陣之間X的回歸系數(shù)矩陣。

    由此可見,主成分回歸模型是對普通的最小二乘估計的改進,首先選取主成分,克服自變量間的多重共線性,然后對所選的主成分進行線性回歸,進而得到主成分回歸方程。

    3 實例應用

    3.1 研究區(qū)概況

    王家壩站系淮河上游總控制站,集水面積30630km2。上游干流河長360km,河道比降0.5‰,年平均降水量800~1200mm,且降水年際變化大,時空分布不均勻。年降水量的60%集中在5~8月,以6、7兩月暴雨次數(shù)較多。產(chǎn)生暴雨的主要天氣系統(tǒng)是西南低渦、切變線、低壓槽和臺風等?;锤缮嫌渭盎茨仙絽^(qū)一般是王家壩洪水的主要來源區(qū)。本項目主要針對息縣、潢川、班臺至王家壩區(qū)間(面積為7110km2)開展研究,見圖 2。

    3.2 概率預報結果

    本研究以新安江模型為確定性預報模型,在此基礎上,采用PCA-HUP模型分別對上述模型預報的可靠度進行定量,并在淮河干流主要控制斷面王家壩進行了實際應用,采取滾動預報方式實現(xiàn)了洪水概率預報。

    將新安江模型的預報結果與實測資料輸入PCA-HUP模型中,其中20場洪水用于相關參數(shù)的率定,8場洪水用于模型驗證。模型相關參數(shù)見表1。以置信度為90%(亦可采用其他置信度值)的預報區(qū)間為例,對概率預報結果進行評估,同時,對流量分布函數(shù)的中位數(shù)Q50進行分位數(shù)評價,率定期模型的模擬精度見表2。

    由表2可知,PCA-HUP模型率定期模擬結果:預報區(qū)間(置信度為90%)覆蓋率較高,且離散度在0.2以內(nèi)。此外,將每一時刻預報量概率分布的中位數(shù)預報與實測流量進行比較,確定性系數(shù)接近于1,洪峰誤差在1%以內(nèi),說明中位數(shù)預報的精度非常高,且從不同預見期的模擬過程線中可以看出,隨著預見期的逐漸增大,區(qū)間離散度呈現(xiàn)出遞增的趨勢。

    采用驗證期8場洪水對概率預報模型進行檢驗,推求預報流量的概率分布,實現(xiàn)王家壩斷面的洪水概率預報。預報精度統(tǒng)計見表3,預報流量過程線以其中兩場為例,預報流量過程線如圖3~圖4所示。

    由表3和2場洪水概率預報過程線可知,PCA-HUP模型(以新安江模型為確定性預報模型)提供的概率預報結果:預報區(qū)間(置信度為90%)覆蓋率在87%以上,且離散度在0.2以內(nèi),說明在相對較小的區(qū)間寬度內(nèi),預報區(qū)間仍然能夠覆蓋絕大多數(shù)實測數(shù)據(jù),說明概率預報精度較高。此外,將每一時刻預報量概率分布的中位數(shù)預報與實測流量進行比較,確定性系數(shù)接近于1,洪峰誤差在1%以內(nèi),說明中位數(shù)預報的精度非常高,明顯高于新安江模型預報結果,充分體現(xiàn)了貝葉斯修正原理。從不同預見期的預報過程線中可以看出,隨著預見期的逐漸增大,區(qū)間寬度呈現(xiàn)出遞增的趨勢。

    表2 PCA-HUP模型率定期模擬精度統(tǒng)計表(王家壩,Δt=2h)

    圖2 王家壩斷面控制區(qū)域圖

    表1 PCA-HUP模型參數(shù)表(王家壩)

    圖3 a 19920505號洪水預報過程線圖(王家壩,Δt=2h)

    圖3 b 19920505號洪水預報過程線圖(王家壩,Δt=6h)

    圖3 c 19920505號洪水預報過程線圖(王家壩,Δt=12h)

    圖4 a 19950707號洪水預報過程線圖(王家壩,Δt=2h)

    圖4 b 19950707號洪水預報過程線圖(王家壩,Δt=6h)

    圖4 c 19950707號洪水預報過程線圖(王家壩,Δt=12h)

    表3 PCA-HUP模型驗證期概率預報精度統(tǒng)計表(王家壩,Δt=2h)

    4 結論

    本文在貝葉斯預報系統(tǒng)(BFS)的基礎上,對其水文不確定性處理器(HUP)進行改進,提出改進的貝葉斯概率預報系統(tǒng)(PCA-HUP),PCA-HUP模型是在正態(tài)空間中對預報誤差進行分析,構造線性似然函數(shù),推求預報流量的后驗分布,實現(xiàn)洪水概率預報。以淮河王家壩斷面的洪水預報為示例進行了應用研究,表明該方法不僅可以提供不同置信度的區(qū)間預報結果,還可以獲得精度更高的定值預報結果。

    該方法可以與任何確定性的水文預報模型相耦合。作為示例,文中對相對誤差進行了正態(tài)假設,但研究方法不局限于此,同樣適用于誤差服從其他分布的情況。限于研究區(qū)的預報誤差特征,本文僅考慮了誤差均值隨流量量級的變化規(guī)律,相同思路亦適用于誤差方差等分布參數(shù)隨流量變化的情況■

    猜你喜歡
    王家壩貝葉斯不確定性
    最美王家壩,我的家
    家教世界(2022年28期)2022-10-25 01:02:36
    法律的兩種不確定性
    法律方法(2022年2期)2022-10-20 06:41:56
    最美王家壩,我的家
    價值·機理·路徑:新時代弘揚王家壩精神探析
    英鎊或繼續(xù)面臨不確定性風險
    中國外匯(2019年7期)2019-07-13 05:45:04
    貝葉斯公式及其應用
    具有不可測動態(tài)不確定性非線性系統(tǒng)的控制
    基于貝葉斯估計的軌道占用識別方法
    一種基于貝葉斯壓縮感知的說話人識別方法
    電子器件(2015年5期)2015-12-29 08:43:15
    IIRCT下負二項分布參數(shù)多變點的貝葉斯估計
    国产成人啪精品午夜网站| 一进一出抽搐gif免费好疼| 两个人看的免费小视频| 婷婷亚洲欧美| 日本与韩国留学比较| 三级国产精品欧美在线观看 | 午夜福利在线观看免费完整高清在 | 一进一出抽搐动态| 91在线精品国自产拍蜜月 | 极品教师在线免费播放| 亚洲九九香蕉| 亚洲aⅴ乱码一区二区在线播放| a级毛片在线看网站| 性色avwww在线观看| 可以在线观看的亚洲视频| 日本黄色片子视频| 亚洲欧洲精品一区二区精品久久久| 九色国产91popny在线| 国产成人精品久久二区二区免费| 欧美日韩瑟瑟在线播放| 国产亚洲av高清不卡| 国产av在哪里看| 午夜福利免费观看在线| 啪啪无遮挡十八禁网站| 久久亚洲真实| 国产伦在线观看视频一区| 国产av麻豆久久久久久久| 亚洲av中文字字幕乱码综合| 久久久久久久久免费视频了| 欧美一级a爱片免费观看看| 又黄又粗又硬又大视频| 国产乱人伦免费视频| 日本黄大片高清| 无遮挡黄片免费观看| www日本黄色视频网| 欧美在线黄色| 一区福利在线观看| 亚洲熟妇熟女久久| 免费搜索国产男女视频| 欧美色欧美亚洲另类二区| 午夜精品一区二区三区免费看| av黄色大香蕉| 特大巨黑吊av在线直播| 欧美3d第一页| 黄片大片在线免费观看| 欧美色视频一区免费| 无人区码免费观看不卡| 久久久久国产一级毛片高清牌| 18美女黄网站色大片免费观看| 99久久无色码亚洲精品果冻| 老鸭窝网址在线观看| 国产精品久久视频播放| 淫秽高清视频在线观看| 黄色视频,在线免费观看| 十八禁网站免费在线| av在线天堂中文字幕| 丝袜人妻中文字幕| 免费在线观看成人毛片| 欧洲精品卡2卡3卡4卡5卡区| 国产探花在线观看一区二区| 国产免费男女视频| 嫩草影视91久久| 男女做爰动态图高潮gif福利片| 免费高清视频大片| 国产精品久久久久久精品电影| 99精品欧美一区二区三区四区| 日本a在线网址| 成人性生交大片免费视频hd| 一本精品99久久精品77| 不卡av一区二区三区| 精品电影一区二区在线| 一区福利在线观看| 日日干狠狠操夜夜爽| 亚洲国产色片| 特大巨黑吊av在线直播| 欧美另类亚洲清纯唯美| 国产精品久久电影中文字幕| 午夜精品一区二区三区免费看| 超碰成人久久| 欧美日韩黄片免| 男女下面进入的视频免费午夜| 国产成人福利小说| 国产精品久久久久久人妻精品电影| 天堂网av新在线| 91麻豆av在线| 小蜜桃在线观看免费完整版高清| 91字幕亚洲| 亚洲狠狠婷婷综合久久图片| 亚洲国产精品久久男人天堂| 悠悠久久av| 中国美女看黄片| 天天一区二区日本电影三级| 国产高清三级在线| 日日夜夜操网爽| 日本免费一区二区三区高清不卡| 叶爱在线成人免费视频播放| 法律面前人人平等表现在哪些方面| 成人国产综合亚洲| 在线永久观看黄色视频| 别揉我奶头~嗯~啊~动态视频| 免费一级毛片在线播放高清视频| 中文亚洲av片在线观看爽| 久久这里只有精品中国| 真人一进一出gif抽搐免费| 国产激情欧美一区二区| 国产三级在线视频| 久久天躁狠狠躁夜夜2o2o| 久久午夜综合久久蜜桃| 久久精品人妻少妇| 欧美在线一区亚洲| 久久热在线av| 亚洲真实伦在线观看| 俄罗斯特黄特色一大片| 一卡2卡三卡四卡精品乱码亚洲| 欧美日韩瑟瑟在线播放| 美女黄网站色视频| 精品欧美国产一区二区三| 黄色丝袜av网址大全| 亚洲 欧美 日韩 在线 免费| 麻豆国产97在线/欧美| 久久这里只有精品19| 国产成+人综合+亚洲专区| 97人妻精品一区二区三区麻豆| cao死你这个sao货| 最近最新中文字幕大全电影3| 国产精品综合久久久久久久免费| 淫秽高清视频在线观看| 欧美国产日韩亚洲一区| 国产乱人视频| 国产高清有码在线观看视频| www.熟女人妻精品国产| 亚洲av成人不卡在线观看播放网| 一个人观看的视频www高清免费观看 | 桃红色精品国产亚洲av| 国产熟女xx| 国产精品香港三级国产av潘金莲| 男插女下体视频免费在线播放| 免费在线观看日本一区| 亚洲成人中文字幕在线播放| 久久久色成人| 免费在线观看视频国产中文字幕亚洲| 久久婷婷人人爽人人干人人爱| 久久九九热精品免费| 在线国产一区二区在线| 丰满的人妻完整版| 日韩有码中文字幕| 岛国在线观看网站| 久久香蕉国产精品| 国产精品永久免费网站| 久久久久亚洲av毛片大全| 国产黄a三级三级三级人| 免费在线观看成人毛片| 国产精品久久久久久久电影 | 午夜激情欧美在线| 不卡av一区二区三区| 国产91精品成人一区二区三区| 亚洲成人久久爱视频| 亚洲自偷自拍图片 自拍| 久久婷婷人人爽人人干人人爱| 国产精品久久视频播放| 丝袜人妻中文字幕| 丁香欧美五月| 九九久久精品国产亚洲av麻豆 | 国产成人精品无人区| 欧美国产日韩亚洲一区| 女人高潮潮喷娇喘18禁视频| 俺也久久电影网| 真人做人爱边吃奶动态| 午夜免费观看网址| 国产97色在线日韩免费| 国产极品精品免费视频能看的| a级毛片a级免费在线| 色综合欧美亚洲国产小说| 一级黄色大片毛片| 叶爱在线成人免费视频播放| 精品国产乱子伦一区二区三区| 男女午夜视频在线观看| 成人无遮挡网站| 亚洲av中文字字幕乱码综合| 1024香蕉在线观看| 免费在线观看成人毛片| 国产一区二区三区视频了| 俄罗斯特黄特色一大片| 99热精品在线国产| 一进一出好大好爽视频| 亚洲精品456在线播放app | 首页视频小说图片口味搜索| 国产一级毛片七仙女欲春2| 国产精品精品国产色婷婷| 色尼玛亚洲综合影院| 日本免费一区二区三区高清不卡| 成人午夜高清在线视频| 亚洲午夜理论影院| 男女床上黄色一级片免费看| 亚洲一区高清亚洲精品| 午夜视频精品福利| 麻豆av在线久日| 在线国产一区二区在线| 久久精品国产亚洲av香蕉五月| 可以在线观看的亚洲视频| 老熟妇乱子伦视频在线观看| 国产精品精品国产色婷婷| 欧美大码av| 97碰自拍视频| 好男人在线观看高清免费视频| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久午夜电影| 99热这里只有精品一区 | 日韩人妻高清精品专区| 我要搜黄色片| 99国产精品99久久久久| 国产成年人精品一区二区| 狂野欧美激情性xxxx| 国产一级毛片七仙女欲春2| 少妇裸体淫交视频免费看高清| 国产又黄又爽又无遮挡在线| 村上凉子中文字幕在线| 精品熟女少妇八av免费久了| 成人国产综合亚洲| 免费观看人在逋| 99久久成人亚洲精品观看| 女人高潮潮喷娇喘18禁视频| 国产精华一区二区三区| 白带黄色成豆腐渣| 在线观看美女被高潮喷水网站 | 好男人在线观看高清免费视频| 少妇熟女aⅴ在线视频| 国产黄a三级三级三级人| 99热这里只有精品一区 | 国产精品av久久久久免费| 欧美乱妇无乱码| 一二三四社区在线视频社区8| 搡老妇女老女人老熟妇| 99精品久久久久人妻精品| 窝窝影院91人妻| 亚洲 国产 在线| 欧美性猛交╳xxx乱大交人| 日本撒尿小便嘘嘘汇集6| 久久精品国产综合久久久| 午夜福利18| 精品免费久久久久久久清纯| 国产精品一及| 国产亚洲精品久久久com| 中文字幕久久专区| 亚洲 国产 在线| 噜噜噜噜噜久久久久久91| aaaaa片日本免费| av福利片在线观看| 一夜夜www| 一级毛片精品| 亚洲18禁久久av| 精品日产1卡2卡| 久久欧美精品欧美久久欧美| 久9热在线精品视频| 757午夜福利合集在线观看| 亚洲中文av在线| 久久久久精品国产欧美久久久| 人人妻人人澡欧美一区二区| 热99re8久久精品国产| 女同久久另类99精品国产91| 免费搜索国产男女视频| 欧美日本视频| 一二三四社区在线视频社区8| 亚洲精品乱码久久久v下载方式 | 久久久久久久久久黄片| 国产精品 欧美亚洲| 亚洲乱码一区二区免费版| 男人和女人高潮做爰伦理| 欧美日韩瑟瑟在线播放| 制服丝袜大香蕉在线| 精品久久久久久久久久免费视频| 精品国产美女av久久久久小说| 久久精品国产99精品国产亚洲性色| 午夜福利高清视频| 一级作爱视频免费观看| 麻豆久久精品国产亚洲av| 色精品久久人妻99蜜桃| 激情在线观看视频在线高清| av视频在线观看入口| 91在线观看av| 嫩草影院精品99| 99精品在免费线老司机午夜| 18禁国产床啪视频网站| 丰满人妻熟妇乱又伦精品不卡| 日本与韩国留学比较| 欧美色视频一区免费| 在线观看免费视频日本深夜| 国产精品久久久久久人妻精品电影| 亚洲无线观看免费| 在线免费观看的www视频| 一本久久中文字幕| 亚洲av成人av| 又爽又黄无遮挡网站| 午夜两性在线视频| 久久精品国产清高在天天线| 99在线视频只有这里精品首页| 日日摸夜夜添夜夜添小说| 色在线成人网| 麻豆成人午夜福利视频| 白带黄色成豆腐渣| 国产精品九九99| 亚洲avbb在线观看| h日本视频在线播放| 黄片小视频在线播放| 国产麻豆成人av免费视频| 国产亚洲精品av在线| 国产野战对白在线观看| 长腿黑丝高跟| 国产三级在线视频| 嫩草影院精品99| 成人鲁丝片一二三区免费| 国产亚洲精品久久久久久毛片| 日本a在线网址| 欧美乱码精品一区二区三区| 听说在线观看完整版免费高清| 啦啦啦观看免费观看视频高清| 香蕉av资源在线| 校园春色视频在线观看| 男人舔女人下体高潮全视频| 亚洲男人的天堂狠狠| 91av网一区二区| 免费在线观看成人毛片| 男女午夜视频在线观看| 国产精品亚洲av一区麻豆| 男人舔奶头视频| 午夜精品在线福利| 黄色 视频免费看| 在线观看午夜福利视频| 国产精品久久久久久久电影 | 成人av在线播放网站| 亚洲电影在线观看av| 亚洲五月天丁香| 成年版毛片免费区| 嫩草影院入口| 久久国产精品影院| 一边摸一边抽搐一进一小说| 国产精品99久久99久久久不卡| av视频在线观看入口| 欧美绝顶高潮抽搐喷水| 日本a在线网址| 老司机深夜福利视频在线观看| 久久中文字幕人妻熟女| 国产伦在线观看视频一区| 天天添夜夜摸| 国产免费男女视频| 国产成人精品久久二区二区91| 国产精品亚洲美女久久久| 极品教师在线免费播放| 搞女人的毛片| 国产精品女同一区二区软件 | 国模一区二区三区四区视频 | 18禁黄网站禁片午夜丰满| 老汉色∧v一级毛片| 男人和女人高潮做爰伦理| 亚洲精品在线美女| 99精品久久久久人妻精品| 成人av一区二区三区在线看| 午夜精品久久久久久毛片777| 999久久久精品免费观看国产| 午夜免费激情av| 全区人妻精品视频| 午夜福利视频1000在线观看| 国产99白浆流出| 欧美成狂野欧美在线观看| 成人特级黄色片久久久久久久| 久久久久九九精品影院| 身体一侧抽搐| 狂野欧美白嫩少妇大欣赏| 午夜福利高清视频| x7x7x7水蜜桃| 亚洲国产精品999在线| 国产69精品久久久久777片 | 91老司机精品| 女生性感内裤真人,穿戴方法视频| 一级毛片精品| 国产精品精品国产色婷婷| 国产熟女xx| 观看美女的网站| 亚洲最大成人中文| 色精品久久人妻99蜜桃| 久久精品91蜜桃| 国产亚洲精品久久久com| 欧美日韩黄片免| 欧美色视频一区免费| 国产高清视频在线播放一区| aaaaa片日本免费| 久久人人精品亚洲av| 最新在线观看一区二区三区| 巨乳人妻的诱惑在线观看| 成人亚洲精品av一区二区| 国产精品香港三级国产av潘金莲| 日韩欧美国产在线观看| 97超视频在线观看视频| 美女被艹到高潮喷水动态| 成人特级av手机在线观看| 国产精品电影一区二区三区| 欧美av亚洲av综合av国产av| 国内毛片毛片毛片毛片毛片| 久久久久久国产a免费观看| 国产单亲对白刺激| 黄色成人免费大全| 亚洲精品美女久久久久99蜜臀| 免费看日本二区| 国产高清视频在线播放一区| 欧美xxxx黑人xx丫x性爽| 9191精品国产免费久久| 91麻豆精品激情在线观看国产| 桃色一区二区三区在线观看| 国产一级毛片七仙女欲春2| 91久久精品国产一区二区成人 | 国产精品99久久久久久久久| 久久久久久久久久黄片| 精品日产1卡2卡| 日韩欧美国产一区二区入口| 日本一二三区视频观看| 国产毛片a区久久久久| 国产精品日韩av在线免费观看| 中文在线观看免费www的网站| 男女下面进入的视频免费午夜| 国产成人aa在线观看| 亚洲av成人不卡在线观看播放网| 久久草成人影院| 日本 av在线| 久久午夜综合久久蜜桃| 国产精品九九99| 好看av亚洲va欧美ⅴa在| 亚洲av第一区精品v没综合| 国产乱人伦免费视频| 午夜免费激情av| 久久久久久人人人人人| 女人被狂操c到高潮| 国内揄拍国产精品人妻在线| 一本精品99久久精品77| 女生性感内裤真人,穿戴方法视频| 黑人操中国人逼视频| 搡老妇女老女人老熟妇| 精品不卡国产一区二区三区| 最近在线观看免费完整版| 大型黄色视频在线免费观看| 久久久国产成人精品二区| 国产主播在线观看一区二区| 欧美一级a爱片免费观看看| 成人特级av手机在线观看| 熟妇人妻久久中文字幕3abv| 夜夜爽天天搞| 亚洲一区二区三区色噜噜| 黄色片一级片一级黄色片| 亚洲精品一区av在线观看| 欧美一级毛片孕妇| 精品乱码久久久久久99久播| 成人av一区二区三区在线看| x7x7x7水蜜桃| 国内揄拍国产精品人妻在线| 欧美一级a爱片免费观看看| a级毛片a级免费在线| 日韩国内少妇激情av| 18美女黄网站色大片免费观看| 老汉色∧v一级毛片| 免费在线观看亚洲国产| 男女下面进入的视频免费午夜| 国产又黄又爽又无遮挡在线| 亚洲无线在线观看| 在线十欧美十亚洲十日本专区| 成人无遮挡网站| 欧美一区二区精品小视频在线| 国产97色在线日韩免费| 国产精品久久久久久精品电影| 国内精品久久久久精免费| 嫩草影院精品99| 国产精品一区二区三区四区免费观看 | a在线观看视频网站| 久久久久久人人人人人| 欧美一级毛片孕妇| 深夜精品福利| 亚洲一区二区三区不卡视频| 久久久久国内视频| 伊人久久大香线蕉亚洲五| 特大巨黑吊av在线直播| 五月玫瑰六月丁香| 91在线观看av| 久久天躁狠狠躁夜夜2o2o| 日韩欧美 国产精品| 女人被狂操c到高潮| 一进一出好大好爽视频| 看黄色毛片网站| 99国产综合亚洲精品| 亚洲国产欧美网| 午夜福利免费观看在线| 日韩欧美在线二视频| 香蕉久久夜色| 俄罗斯特黄特色一大片| 好看av亚洲va欧美ⅴa在| 国产成人精品久久二区二区91| 亚洲在线自拍视频| 亚洲精品色激情综合| 国产又黄又爽又无遮挡在线| 国产精品一及| 亚洲avbb在线观看| 他把我摸到了高潮在线观看| 色尼玛亚洲综合影院| 久久天堂一区二区三区四区| 两个人看的免费小视频| 色吧在线观看| 亚洲国产中文字幕在线视频| 久久婷婷人人爽人人干人人爱| 国产人伦9x9x在线观看| 国产精品av视频在线免费观看| 午夜久久久久精精品| 在线免费观看的www视频| 亚洲 欧美一区二区三区| 成在线人永久免费视频| 少妇熟女aⅴ在线视频| 午夜福利在线观看吧| 999精品在线视频| 可以在线观看的亚洲视频| 国产三级中文精品| 国产成人欧美在线观看| 热99在线观看视频| 精品久久久久久成人av| 草草在线视频免费看| 午夜福利成人在线免费观看| 久久亚洲精品不卡| 亚洲国产精品久久男人天堂| 黄色女人牲交| av黄色大香蕉| 免费av毛片视频| 欧美另类亚洲清纯唯美| 久久99热这里只有精品18| 久久伊人香网站| 色综合亚洲欧美另类图片| 亚洲第一电影网av| 国产精品,欧美在线| 国产成人av激情在线播放| 亚洲中文av在线| 黑人操中国人逼视频| 成人午夜高清在线视频| 一级黄色大片毛片| 一本久久中文字幕| 日本黄色视频三级网站网址| 中亚洲国语对白在线视频| 亚洲精品在线美女| 热99re8久久精品国产| 999久久久精品免费观看国产| 欧美成狂野欧美在线观看| 国产99白浆流出| 久久久国产成人免费| 亚洲成人久久性| 一进一出抽搐动态| 久久久久亚洲av毛片大全| 久久久国产成人精品二区| 五月伊人婷婷丁香| 校园春色视频在线观看| 国产伦精品一区二区三区视频9 | 日韩av在线大香蕉| 黄色成人免费大全| 国产亚洲av嫩草精品影院| 免费在线观看影片大全网站| 99久久国产精品久久久| 村上凉子中文字幕在线| 午夜福利成人在线免费观看| 免费av毛片视频| 国内精品久久久久精免费| 久久精品综合一区二区三区| 亚洲无线观看免费| 亚洲五月婷婷丁香| 亚洲人成网站在线播放欧美日韩| 午夜激情福利司机影院| 91麻豆精品激情在线观看国产| 香蕉久久夜色| 真人一进一出gif抽搐免费| 搡老岳熟女国产| 欧美一级a爱片免费观看看| 99国产精品一区二区蜜桃av| 久久精品人妻少妇| 国产aⅴ精品一区二区三区波| 欧美日韩乱码在线| 午夜视频精品福利| 在线播放国产精品三级| 国产精品女同一区二区软件 | 午夜激情欧美在线| 国产男靠女视频免费网站| 久久久国产成人免费| 国产精品 欧美亚洲| 日韩 欧美 亚洲 中文字幕| 国产久久久一区二区三区| 色吧在线观看| 成人欧美大片| 观看美女的网站| 国产高潮美女av| 99久久成人亚洲精品观看| 欧美日韩乱码在线| 国内少妇人妻偷人精品xxx网站 | 男女下面进入的视频免费午夜| 小蜜桃在线观看免费完整版高清| 午夜免费激情av| 亚洲国产看品久久| 9191精品国产免费久久| 久久久色成人| 国产蜜桃级精品一区二区三区| 久久久久免费精品人妻一区二区| 后天国语完整版免费观看| 黄色女人牲交| 99国产精品一区二区三区| 国产熟女xx| 久久午夜亚洲精品久久| avwww免费| 国产亚洲精品一区二区www| 麻豆国产av国片精品| 亚洲美女黄片视频| 国产精品野战在线观看| 18美女黄网站色大片免费观看| 91av网一区二区| 中文字幕av在线有码专区| 亚洲成av人片免费观看| 男女做爰动态图高潮gif福利片|