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

    基于氣象水文耦合的中小河流枯季徑流短期預(yù)報(bào)

    2018-11-29 00:47:02
    中國(guó)農(nóng)村水利水電 2018年11期
    關(guān)鍵詞:天氣預(yù)報(bào)水文站水文

    王 東 升

    (云南省水文水資源局,昆明 650106)

    0 引 言

    在水文預(yù)報(bào)中引入降水預(yù)報(bào)信息,采用氣象水文耦合方法延長(zhǎng)預(yù)報(bào)預(yù)見期,可為防汛減災(zāi)、水資源調(diào)度、水力發(fā)電、生態(tài)環(huán)境保護(hù)提供重要依據(jù),贏得更多的應(yīng)急響應(yīng)時(shí)間,是近年正在形成和發(fā)展的水文預(yù)報(bào)分支,是未來水文預(yù)報(bào)發(fā)展趨勢(shì)[1-3]。張建云院士提出,需要加強(qiáng)水文氣象的學(xué)科交叉和結(jié)合,開展定量降水預(yù)報(bào)(QPF)的研究和應(yīng)用,進(jìn)一步延長(zhǎng)洪水預(yù)報(bào)的預(yù)見期,提高洪水預(yù)報(bào)的效益和作用[4]。降水量預(yù)報(bào)的精度成為洪水預(yù)報(bào)中影響有效預(yù)見期長(zhǎng)短的關(guān)鍵因素[3],近年,水文集合預(yù)報(bào)技術(shù)、應(yīng)用水文氣象概率預(yù)報(bào)的發(fā)展增加了預(yù)報(bào)附加值,并能夠延長(zhǎng)預(yù)警提前時(shí)間[1],全球/區(qū)域數(shù)值天氣預(yù)報(bào)模式已經(jīng)可以給出相對(duì)比較準(zhǔn)確的定量降水預(yù)報(bào),但對(duì)于洪水預(yù)報(bào)定時(shí)、定點(diǎn)、定量降水預(yù)報(bào)的要求,降水預(yù)報(bào)的時(shí)段、雨帶位置與量級(jí)的微量偏差都可能引起截然不同的洪水效應(yīng)[3],中小河流由于流域面積小,地形地貌復(fù)雜,下墊面對(duì)降雨的形成影響較大,汛期局部性、單點(diǎn)性暴雨常引發(fā)洪水,數(shù)值天氣預(yù)報(bào)模式更難以提供可靠的預(yù)報(bào)成果??菁咎鞖庀到y(tǒng)相對(duì)較為單一,受地形地貌影響小,數(shù)值天氣預(yù)報(bào)模式預(yù)報(bào)成果準(zhǔn)確程度較高。

    開展枯水期徑流過程的模擬對(duì)水資源配置、生態(tài)環(huán)境需水量的研究、抗旱救災(zāi)等十分重要,但一般情況下,僅利用降雨徑流模型模擬出流域的日徑流過程,對(duì)枯水期的徑流過程很少給予關(guān)注,也很少對(duì)模擬結(jié)果的精度進(jìn)行專門地評(píng)價(jià),現(xiàn)有的研究往往將豐水期和枯水期作為一個(gè)整體進(jìn)行考慮,先模擬出流域的整體日徑流過程后,再對(duì)枯水徑流進(jìn)行分析[5],結(jié)合數(shù)值天氣預(yù)報(bào)結(jié)果開展中小河流短期徑流預(yù)報(bào)研究成果更為稀少。2006年,中國(guó)建立了T213全球中期集合預(yù)報(bào)系統(tǒng)和中尺度區(qū)域集合預(yù)報(bào)系統(tǒng),預(yù)報(bào)時(shí)效達(dá)10~16 d[1],中央氣象臺(tái)對(duì)外公布未來7 d數(shù)值天氣預(yù)報(bào)結(jié)果。本文以數(shù)值天氣預(yù)報(bào)1~7 d降水量、推算的蒸發(fā)量預(yù)報(bào)結(jié)果及實(shí)時(shí)水文觀測(cè)數(shù)據(jù)為資料,以中小河流枯季徑流為分析對(duì)象,優(yōu)選預(yù)報(bào)因子,基于二次回歸正交旋轉(zhuǎn)設(shè)計(jì),研究枯季未來3、5、7 d枯季徑流預(yù)報(bào)模型,探索氣象水文耦合在中小河流枯季徑流預(yù)報(bào)方法。

    1 研究方法

    影響徑流量的因素中最主要的因素是降雨量、初始流域蓄水量和流域蒸發(fā)量[6],水文模型是根據(jù)已知的降雨時(shí)空分布、初始狀態(tài)、模型的結(jié)構(gòu)和參數(shù)等推求流域出口斷面流量過程[7]。南方地區(qū)中小河流枯季徑流以地下水補(bǔ)給為主,徑流量大小主要取決于流域蓄水的多寡,枯季降水將增加流域蓄水,蒸發(fā)將消耗流域蓄水。據(jù)此,以數(shù)值天氣預(yù)報(bào)中短期降水量、推算的蒸發(fā)量預(yù)報(bào)結(jié)果及水文觀測(cè)數(shù)據(jù)為資料,基于二次回歸正交旋轉(zhuǎn)設(shè)計(jì),構(gòu)建短期徑流預(yù)報(bào)模型,數(shù)據(jù)流程見圖1。初始流域蓄水量通過流域前期實(shí)測(cè)河道徑流量、降水量、蒸發(fā)量等因子集反映,預(yù)見期內(nèi)降水量依據(jù)數(shù)值天氣預(yù)報(bào)結(jié)果因子集反映,預(yù)見期內(nèi)蒸發(fā)主要和天氣狀況相關(guān),同天氣狀況下變化不大,可采用數(shù)值天氣預(yù)報(bào)結(jié)果結(jié)合前期同天氣狀況下蒸發(fā)觀測(cè)數(shù)據(jù)推算構(gòu)建蒸發(fā)預(yù)報(bào)因子集。

    圖1 基于氣象水文耦合的二次回歸正交旋轉(zhuǎn)設(shè)計(jì)短期徑流預(yù)報(bào)模型數(shù)據(jù)流程Fig.1 Data flow design of short-term runoff forecasting model based on quadratic regression orthogonal rotation based on meteorological and hydrological coupling

    二次回歸正交旋轉(zhuǎn)設(shè)計(jì)通過統(tǒng)計(jì)軟件對(duì)試驗(yàn)結(jié)果進(jìn)行非線性數(shù)學(xué)模型擬合,是一種將正交、回歸、均有和較高飽和程度融為一體的試驗(yàn)設(shè)計(jì)方法[8],在正交設(shè)計(jì)的基礎(chǔ)上,利用回歸分析,在給出的因素和指標(biāo)之間找出一個(gè)明確的函數(shù)表達(dá)式,建立因果關(guān)系的數(shù)學(xué)模型,以便定量地描述各因素對(duì)指標(biāo)的作用,并用該數(shù)學(xué)模型預(yù)測(cè)[11],該方法在生物學(xué)研究中廣泛應(yīng)用[8-11],二次回歸正交旋轉(zhuǎn)設(shè)計(jì)徑流預(yù)報(bào)模型數(shù)學(xué)表達(dá)式見式(1)。

    (1)

    式中:Qt為預(yù)報(bào)徑流量,t=1,2,…,n;xi為優(yōu)選的預(yù)報(bào)因子;m為預(yù)報(bào)因子個(gè)數(shù);bi、bii、bij均為參數(shù)項(xiàng);b0為常量。

    模型擬定過程中,Qt使用實(shí)測(cè)值,xi可從流域前期實(shí)測(cè)河道徑流量、降水量、蒸發(fā)量,預(yù)見期降水量、蒸發(fā)量通過因子分析法、經(jīng)驗(yàn)法等優(yōu)選確定,為已知值。因可能有多個(gè)xi形成多種組合通過相關(guān)性及顯著性檢驗(yàn),各xi間可能存在多重線性關(guān)系,會(huì)給模型帶來不合理的解釋,致模型的穩(wěn)定性差,造成預(yù)報(bào)模型可靠性差、精度低。為了得到一個(gè)可靠的預(yù)報(bào)模型,需要有效地從眾多影響Qt的因素中挑選貢獻(xiàn)最大的變量,建立最優(yōu)的預(yù)報(bào)模型,這樣的最優(yōu)模型包含所有對(duì)Qt有顯著影響的變量。在此推薦采用逐步回歸分析法[14-16]推求式(1)最優(yōu)預(yù)報(bào)模型結(jié)構(gòu)、預(yù)報(bào)因子及參數(shù)值,該方法是從一個(gè)自變量開始,視自變量對(duì)Qt作用的顯著程度,從大到小依次逐個(gè)引入模型,當(dāng)引入的自變量由于后面變量的引入而變得不顯著時(shí),要將其剔除,引入一個(gè)自變量或從模型中剔除一個(gè)自變量,為逐步回歸的一步,每一步都要進(jìn)行F檢驗(yàn)[11],以確保每次引入新的變量有統(tǒng)計(jì)意義。

    本文使用DPS數(shù)據(jù)處理系統(tǒng)逐步回歸分析法率定模型。限于篇幅,模型率定計(jì)算方法詳見參考文獻(xiàn)《DPS數(shù)據(jù)處理系統(tǒng) 第二卷 現(xiàn)代統(tǒng)計(jì)及數(shù)據(jù)挖掘》[11]關(guān)于逐步回歸分析內(nèi)容,在此不再?gòu)?fù)述。

    2 實(shí)例應(yīng)用

    2.1 流域基本情況及資料

    南汀河流域處于云貴高原西部邊緣,橫斷山脈南延段,流域內(nèi)姑老河水文站(東經(jīng)99°14′、北緯23°41′)以上流域面積4 185 km2,海拔高差2 984 m,屬中山地區(qū)典型的亞熱帶濕潤(rùn)立體氣候類型,多年平均降水量1 600 mm,其中汛期(5-10)月份降水量約占全年降水量的90%,多年平均產(chǎn)水量56.8 億m3,多年平均徑流深701 mm。流域枯季徑流以地下水補(bǔ)給為主,徑流量大小取決于流域蓄水的多寡,同時(shí),枯季降水將增加流域蓄水,蒸發(fā)將消耗流域蓄水,區(qū)域內(nèi)水旱災(zāi)害頻繁,通常是“三年一小旱,五年一大旱”。

    本文選取南汀河姑老河水文站作為研究對(duì)象,資料包括姑老河水文站及上游6個(gè)雨量站1999-2003年、2006-2008年枯季日雨量資料、日蒸發(fā)資料及對(duì)應(yīng)的實(shí)測(cè)日徑流資料,均為云南省水文水資源局整編成果,資料來源可靠。1999-2008年10 a內(nèi),南汀河流域姑老河水文站以上無對(duì)徑流產(chǎn)生較大影響的中大型水利水電工程建設(shè),流域內(nèi)土壤植被等產(chǎn)匯流條件未發(fā)生較為明顯的變化,2002年以來,姑老河水文站降雨減少較為明顯,但徑流與降雨變化趨勢(shì)一致[17],流域產(chǎn)匯流具有較好的一致性,本文選取1999-2003年、2006-2008年作為研究對(duì)象,資料可靠。

    2.2 預(yù)報(bào)模型擬定

    采用DPS數(shù)據(jù)處理系統(tǒng)軟件提供的逐步回歸分析法、收集到的1999-2003年枯季日雨量資料、日蒸發(fā)資料及對(duì)應(yīng)的實(shí)測(cè)日徑流資料分析確定最優(yōu)預(yù)報(bào)模型結(jié)構(gòu)、預(yù)報(bào)因子及參數(shù)值。姑老河站基于數(shù)值天氣預(yù)報(bào)的未來3、5、7 d枯季徑流預(yù)報(bào)模型結(jié)構(gòu)、參數(shù)值、預(yù)報(bào)因子優(yōu)選結(jié)果及模型見式(2)、式(3)、式(4)。

    (2)

    (3)

    (4)

    2.3 應(yīng)用結(jié)果及分析

    采用南汀河姑老河水文站及上游6個(gè)雨量站點(diǎn)1999-2003年枯季雨量資料、蒸發(fā)資料及相應(yīng)實(shí)測(cè)徑流資料依據(jù)預(yù)報(bào)模型進(jìn)行擬合,2006-2008年枯季雨量資料、蒸發(fā)資料及相應(yīng)實(shí)測(cè)徑流資料對(duì)預(yù)報(bào)模型進(jìn)行檢驗(yàn),結(jié)果見圖2、圖3、圖4,由于篇幅限制,本文使用1999年11月-2000年4月、2007年11月-2008年4月未來3、5、7 d平均徑流預(yù)報(bào)模擬結(jié)果與實(shí)際觀測(cè)值比較作為示例說明擬合及檢驗(yàn)情況。

    圖2 姑老河水文站未來3 d徑流模型擬合值與觀測(cè)值比較Fig.2 Comparison of the fitted values and observations of the runoff model of the Gulaohe Hydrological station in the next 3 days

    圖3 姑老河水文站未來5 d徑流模型擬合值與觀測(cè)值比較Fig.3 Comparison of the fitted values and observations of the runoff model of the Gulaohe Hydrological station in the next 5 days

    圖4 姑老河水文站枯季未來7 d徑流模型擬合值與觀測(cè)值比較Fig.4 Comparison of the fitted values and observations of the runoff model of the Gulaohe Hydrological station in the next 7 days

    根據(jù)《水文情報(bào)預(yù)報(bào)規(guī)范(GBT22482-2008)》,使用確定性系數(shù)、絕對(duì)誤差、相對(duì)誤差、合格率對(duì)姑老河站未來3、5、7 d徑流預(yù)報(bào)結(jié)果進(jìn)行分析。一次預(yù)報(bào)的誤差小于許可誤差時(shí),為合格預(yù)報(bào),徑流預(yù)報(bào)許可誤差取徑流實(shí)測(cè)值的20%。

    確定性系數(shù):

    (5)

    絕對(duì)誤差:

    ΔQ=QC(i)-Q0(i)

    (6)

    相對(duì)誤差:

    δ=ΔQ/Q0(i)×100%

    (7)

    合格率:

    (8)

    式中:ΔQ為絕對(duì)誤差;δ為相對(duì)誤差;QR為合格率;n為合格預(yù)報(bào)次數(shù);m為預(yù)報(bào)總次數(shù)。

    表1為未來3 d徑流模擬預(yù)報(bào)精度分析結(jié)果。由表1知,6個(gè)枯季模擬預(yù)報(bào)時(shí)段DC值最大為0.992、最小為0.948、平均為0.974,平均絕對(duì)誤差最大值為2.59 m3/s、最小為1.32 m3/s、平均為1.96 m3/s,平均相對(duì)誤差最大為5.6%、最小為2.9%、平均為4.2%。預(yù)報(bào)合格率最大為99.3%、最小為97.9%、平均為98.8%,達(dá)到甲級(jí)精度等級(jí)。

    表2為未來5 d徑流模擬預(yù)報(bào)精度分析結(jié)果。由表2知,6個(gè)枯季模擬預(yù)報(bào)時(shí)段DC值最大為0.989、最小為0.909、平均為0.960,平均絕對(duì)誤差最大值為3.10 m3/s、最小為1.49 m3/s、平均為2.36 m3/s,平均相對(duì)誤差最大為7.4%、最小為3.2%、平均為5.1%。預(yù)報(bào)合格率最大為100%、最小為92.0%、平均為97.8%,達(dá)到甲級(jí)精度等級(jí)。

    表3為未來7 d徑流模擬預(yù)報(bào)精度分析結(jié)果。由表3知,6個(gè)枯季模擬預(yù)報(bào)時(shí)段DC值最大為0.988、最小為0.865、平均為0.947,平均絕對(duì)誤差最大值為3.53 m3/s、最小為1.56 m3/s、平均為2.61 m3/s,平均相對(duì)誤差最大為8.8%、最小為3.3%、平均為5.7%。預(yù)報(bào)合格率最大為100%、最小為89.7%、平均為97.5%,達(dá)到甲級(jí)精度等級(jí)。

    由表1至表3,隨著預(yù)見期延長(zhǎng),預(yù)報(bào)精度有所下降,但下降不明顯,均滿足甲級(jí)預(yù)報(bào)精度要求,可用于對(duì)外發(fā)布正式預(yù)報(bào)。

    表1 姑老河水文站枯季未來3天徑流預(yù)報(bào)結(jié)果分析Tab.1 Analysis of the runoff forecast results of the Gulaohe Hydrological station in the next 3 days

    3 結(jié) 語(yǔ)

    (1)本文以數(shù)值天氣預(yù)報(bào)及前期實(shí)測(cè)河道徑流量、降水量、蒸發(fā)量為輸入,優(yōu)選預(yù)報(bào)因子,基于二次回歸正交旋轉(zhuǎn)設(shè)計(jì),建立了中小河流枯季徑流預(yù)報(bào)模型,探索了基于水文氣象耦合中小河流枯季短期徑流預(yù)報(bào)方法,經(jīng)在南汀河姑老河水文站未來3、5、7 d枯季徑流模擬中應(yīng)用,效果理想,精度較高,雖然隨著預(yù)見期延長(zhǎng),模型預(yù)報(bào)精度有所下降,但下降不明顯。

    (2)本文使用的短期徑流預(yù)報(bào)方法選用的預(yù)報(bào)因子為預(yù)報(bào)模型提供了良好的物理基礎(chǔ),同時(shí)降低了模型結(jié)構(gòu)的復(fù)雜程度及計(jì)算量,在取得良好的預(yù)報(bào)精度同時(shí),預(yù)見期由基于實(shí)測(cè)降雨的15 h左右分別延長(zhǎng)至3、5、7 d,可滿足枯季水資源科學(xué)調(diào)度、抗旱救災(zāi)、水生態(tài)保護(hù)決策重點(diǎn)關(guān)注的不同預(yù)見期短期徑流預(yù)報(bào)的需要,且該方法操作簡(jiǎn)單,易于實(shí)現(xiàn),具有較好的推廣價(jià)值。

    表2 姑老河水文站枯季未來5天徑流預(yù)報(bào)結(jié)果分析Tab.2 Analysis of the runoff forecast results of the Gulaohe Hydrological station in the next 5 days

    表3 姑老河水文站枯季未來7天徑流預(yù)報(bào)結(jié)果分析Tab.3 Analysis of the runoff forecast results of the Gulaohe Hydrological station in the next 7 days

    (3)汛期數(shù)值天氣預(yù)報(bào)模式預(yù)報(bào)成果因無法滿足中小河流洪水預(yù)報(bào)對(duì)定時(shí)、定點(diǎn)、定量降水預(yù)報(bào)的要求,一直以來主要用于洪水參考性估報(bào),枯季由于洪水災(zāi)害較少,關(guān)注度較低,需求稍顯不足,同樣相關(guān)研究成果較少。本文考慮到枯季數(shù)值天氣預(yù)報(bào)模式預(yù)報(bào)成果精度相對(duì)較高,枯季對(duì)短期徑流預(yù)報(bào)的需求大于洪水預(yù)報(bào)的實(shí)際情況,探索了中小河流氣象水文耦合短期徑流預(yù)報(bào)的方法,通過南汀河姑老河水文站實(shí)例應(yīng)用顯示,氣象水文耦合在中小河流枯季徑流預(yù)報(bào)中,可取得良好的預(yù)報(bào)精度,應(yīng)用前景較為廣闊。

    (4)本模型適用于汛期降水集中、占比大,枯季降水稀少、天氣系統(tǒng)相對(duì)較為單一的季風(fēng)氣候區(qū)的枯季徑流預(yù)報(bào)。由于本文使用了一個(gè)代表流域作為分析對(duì)象,且用于模型建立及檢驗(yàn)的資料較少,還有待進(jìn)一步深化研究。

    猜你喜歡
    天氣預(yù)報(bào)水文站水文
    天氣預(yù)報(bào)員
    2022年《中國(guó)水文年報(bào)》發(fā)布
    SL流量計(jì)在特殊河段的應(yīng)用——以河源水文站為例
    水文
    水文水資源管理
    天氣預(yù)報(bào)的前世今生
    水文
    郭家屯水文站單斷沙關(guān)系分析
    草壩水文站兩種蒸發(fā)器對(duì)比分析
    中期天氣預(yù)報(bào)
    久久97久久精品| 青春草国产在线视频| 精品亚洲乱码少妇综合久久| 亚洲精品美女久久av网站| 边亲边吃奶的免费视频| 国产一区二区在线观看av| 欧美日韩视频高清一区二区三区二| 99久久精品国产国产毛片| 亚洲精品国产av蜜桃| 日日摸夜夜添夜夜添av毛片| 欧美激情 高清一区二区三区| 亚洲久久久国产精品| 久久女婷五月综合色啪小说| 22中文网久久字幕| 久久久久精品久久久久真实原创| 国产国拍精品亚洲av在线观看| 午夜老司机福利剧场| 国产片特级美女逼逼视频| 天堂俺去俺来也www色官网| 一边摸一边做爽爽视频免费| 国产av一区二区精品久久| 日韩,欧美,国产一区二区三区| 国产男女超爽视频在线观看| 夫妻性生交免费视频一级片| av网站免费在线观看视频| 成人手机av| 大话2 男鬼变身卡| 九九久久精品国产亚洲av麻豆| 18禁裸乳无遮挡动漫免费视频| 欧美国产精品一级二级三级| 啦啦啦中文免费视频观看日本| 亚洲欧美清纯卡通| 建设人人有责人人尽责人人享有的| 丝瓜视频免费看黄片| 老女人水多毛片| 国产成人午夜福利电影在线观看| 亚洲av二区三区四区| 满18在线观看网站| 视频区图区小说| 日韩av免费高清视频| 日本猛色少妇xxxxx猛交久久| 男女高潮啪啪啪动态图| 我要看黄色一级片免费的| 亚洲少妇的诱惑av| 国产毛片在线视频| 伊人亚洲综合成人网| 哪个播放器可以免费观看大片| 日本vs欧美在线观看视频| 黑人欧美特级aaaaaa片| 青春草国产在线视频| 久久久国产欧美日韩av| 女人久久www免费人成看片| 欧美精品一区二区大全| 久久久久久伊人网av| 国产视频首页在线观看| 丝袜美足系列| 亚洲精品第二区| 高清毛片免费看| 欧美老熟妇乱子伦牲交| 人妻夜夜爽99麻豆av| av免费观看日本| 一边亲一边摸免费视频| 一二三四中文在线观看免费高清| 最近手机中文字幕大全| 蜜桃久久精品国产亚洲av| 男人添女人高潮全过程视频| 伦理电影大哥的女人| 亚洲国产精品999| 中文字幕人妻丝袜制服| 亚洲欧美成人综合另类久久久| 精品久久蜜臀av无| 国产av国产精品国产| 久久国内精品自在自线图片| 欧美丝袜亚洲另类| 夜夜爽夜夜爽视频| 免费av不卡在线播放| 成人毛片60女人毛片免费| 免费人妻精品一区二区三区视频| 最新中文字幕久久久久| 制服人妻中文乱码| 狠狠精品人妻久久久久久综合| 免费人成在线观看视频色| 高清黄色对白视频在线免费看| 观看av在线不卡| 国产综合精华液| 亚洲国产色片| 十八禁高潮呻吟视频| 久久婷婷青草| 亚洲精品久久久久久婷婷小说| 菩萨蛮人人尽说江南好唐韦庄| 国产一区亚洲一区在线观看| 欧美日韩av久久| 久久精品国产亚洲av天美| 美女中出高潮动态图| 欧美精品一区二区大全| 亚洲国产欧美在线一区| 黄色毛片三级朝国网站| 日本91视频免费播放| 国产在线视频一区二区| 这个男人来自地球电影免费观看 | 男女边吃奶边做爰视频| 少妇猛男粗大的猛烈进出视频| 九九爱精品视频在线观看| 免费看不卡的av| 人体艺术视频欧美日本| 男女边吃奶边做爰视频| 午夜激情av网站| 久久毛片免费看一区二区三区| 夫妻性生交免费视频一级片| 观看av在线不卡| 热re99久久国产66热| 久久精品久久久久久久性| 国产成人一区二区在线| 99久久中文字幕三级久久日本| 久久97久久精品| 久久精品国产鲁丝片午夜精品| 亚洲在久久综合| 精品人妻偷拍中文字幕| 十分钟在线观看高清视频www| 美女福利国产在线| 蜜桃国产av成人99| av视频免费观看在线观看| 中国国产av一级| 黄色毛片三级朝国网站| 日本欧美国产在线视频| 国产乱来视频区| 美女中出高潮动态图| kizo精华| 国产永久视频网站| 一级a做视频免费观看| 丰满少妇做爰视频| 亚洲国产精品国产精品| 波野结衣二区三区在线| 亚洲精品久久久久久婷婷小说| 爱豆传媒免费全集在线观看| 午夜激情av网站| 亚洲欧美一区二区三区国产| 少妇精品久久久久久久| 国产欧美亚洲国产| 亚洲国产欧美在线一区| 男人爽女人下面视频在线观看| 特大巨黑吊av在线直播| 在线天堂最新版资源| 欧美日韩视频精品一区| 欧美日韩视频高清一区二区三区二| 内地一区二区视频在线| 欧美最新免费一区二区三区| 久久久久精品久久久久真实原创| 亚洲精品久久午夜乱码| 一本一本综合久久| 高清毛片免费看| 欧美+日韩+精品| 精品人妻偷拍中文字幕| 在线观看国产h片| 一本一本综合久久| 欧美老熟妇乱子伦牲交| 日韩一本色道免费dvd| 日韩大片免费观看网站| 一区二区三区免费毛片| 日本欧美视频一区| 国产有黄有色有爽视频| 九九久久精品国产亚洲av麻豆| 一区二区三区四区激情视频| 嘟嘟电影网在线观看| 国产毛片在线视频| 精品99又大又爽又粗少妇毛片| 中文字幕人妻熟人妻熟丝袜美| 精品酒店卫生间| 3wmmmm亚洲av在线观看| 91精品国产国语对白视频| 性色avwww在线观看| 国产精品久久久久久精品古装| a级片在线免费高清观看视频| 我的老师免费观看完整版| 亚洲伊人久久精品综合| 丰满少妇做爰视频| 国产精品麻豆人妻色哟哟久久| .国产精品久久| 精品国产乱码久久久久久小说| 国产成人av激情在线播放 | av天堂久久9| 国产精品久久久久久精品电影小说| 国产色爽女视频免费观看| 在线观看一区二区三区激情| 在线观看一区二区三区激情| av.在线天堂| 久久久久久伊人网av| 18禁在线无遮挡免费观看视频| 亚洲国产成人一精品久久久| 国产探花极品一区二区| av免费观看日本| 久久青草综合色| 亚洲高清免费不卡视频| 亚洲国产精品国产精品| 国产午夜精品久久久久久一区二区三区| 国产精品免费大片| 亚洲精品成人av观看孕妇| 成年女人在线观看亚洲视频| 午夜日本视频在线| 国产片内射在线| 欧美日韩av久久| 国产淫语在线视频| 边亲边吃奶的免费视频| 中文字幕最新亚洲高清| 性色av一级| 日日摸夜夜添夜夜添av毛片| 免费观看a级毛片全部| 黑人欧美特级aaaaaa片| 香蕉精品网在线| 人妻一区二区av| 男女边摸边吃奶| 免费av中文字幕在线| 久久久久网色| 99久久人妻综合| 狠狠婷婷综合久久久久久88av| 免费人妻精品一区二区三区视频| 中文精品一卡2卡3卡4更新| 国产精品嫩草影院av在线观看| 精品亚洲成国产av| 日韩人妻高清精品专区| 国国产精品蜜臀av免费| 欧美丝袜亚洲另类| 日韩熟女老妇一区二区性免费视频| 欧美日韩视频精品一区| 久久青草综合色| 国产精品久久久久久久久免| 欧美人与善性xxx| 街头女战士在线观看网站| 久久影院123| 久久久久国产精品人妻一区二区| 久热这里只有精品99| 91午夜精品亚洲一区二区三区| 一级a做视频免费观看| 亚洲精品成人av观看孕妇| 亚洲精品一二三| 这个男人来自地球电影免费观看 | 人妻人人澡人人爽人人| 久久久a久久爽久久v久久| 夜夜看夜夜爽夜夜摸| av免费在线看不卡| 午夜激情av网站| av.在线天堂| 亚洲欧洲日产国产| 日日啪夜夜爽| 欧美日韩综合久久久久久| 国产成人91sexporn| 最后的刺客免费高清国语| 极品少妇高潮喷水抽搐| 高清黄色对白视频在线免费看| 国产亚洲欧美精品永久| 久久久精品区二区三区| 成人免费观看视频高清| 18禁在线播放成人免费| 亚洲精品456在线播放app| 免费大片黄手机在线观看| 久久久久精品久久久久真实原创| 一本大道久久a久久精品| 国产日韩欧美亚洲二区| 中文天堂在线官网| 欧美少妇被猛烈插入视频| 亚洲五月色婷婷综合| 国产精品国产av在线观看| 欧美日韩成人在线一区二区| 日韩成人av中文字幕在线观看| 国产精品久久久久久精品古装| 黄色一级大片看看| 国产精品国产三级国产av玫瑰| 女性被躁到高潮视频| 成人午夜精彩视频在线观看| 九九在线视频观看精品| 啦啦啦视频在线资源免费观看| 久久精品国产鲁丝片午夜精品| 国产一区二区三区av在线| 女性被躁到高潮视频| 亚洲av二区三区四区| 免费大片18禁| 热99久久久久精品小说推荐| 午夜久久久在线观看| 国产亚洲精品第一综合不卡 | 91久久精品国产一区二区三区| 精品人妻一区二区三区麻豆| 亚洲美女视频黄频| 亚洲五月色婷婷综合| 高清视频免费观看一区二区| 街头女战士在线观看网站| 哪个播放器可以免费观看大片| 插阴视频在线观看视频| 欧美变态另类bdsm刘玥| 26uuu在线亚洲综合色| 51国产日韩欧美| 成人亚洲精品一区在线观看| 人人妻人人澡人人看| 日日爽夜夜爽网站| 国产成人一区二区在线| 亚洲一级一片aⅴ在线观看| 毛片一级片免费看久久久久| 欧美成人午夜免费资源| 欧美精品一区二区大全| 夜夜骑夜夜射夜夜干| 久久精品国产亚洲av天美| 亚洲成色77777| 久久青草综合色| 久久久久久久国产电影| 亚洲无线观看免费| 一区二区av电影网| 国产无遮挡羞羞视频在线观看| 妹子高潮喷水视频| 国产毛片在线视频| 午夜av观看不卡| 国产精品免费大片| 夜夜爽夜夜爽视频| 国产免费一级a男人的天堂| 亚洲欧美一区二区三区黑人 | 久久午夜福利片| 日韩亚洲欧美综合| 女的被弄到高潮叫床怎么办| 自拍欧美九色日韩亚洲蝌蚪91| 免费大片18禁| 丰满迷人的少妇在线观看| 免费高清在线观看日韩| www.av在线官网国产| 亚洲三级黄色毛片| 黄色怎么调成土黄色| 久久99热这里只频精品6学生| 午夜福利影视在线免费观看| 国产精品蜜桃在线观看| 国产免费视频播放在线视频| 国产av国产精品国产| 热re99久久国产66热| 97精品久久久久久久久久精品| 老熟女久久久| 大又大粗又爽又黄少妇毛片口| 欧美人与善性xxx| 在线天堂最新版资源| 特大巨黑吊av在线直播| 高清在线视频一区二区三区| 欧美人与善性xxx| 一级,二级,三级黄色视频| 少妇人妻久久综合中文| 色视频在线一区二区三区| 精品一区二区三区视频在线| 蜜桃国产av成人99| 伊人久久精品亚洲午夜| 你懂的网址亚洲精品在线观看| 26uuu在线亚洲综合色| 飞空精品影院首页| 永久免费av网站大全| 黑人猛操日本美女一级片| av专区在线播放| 少妇丰满av| 久久久久网色| 亚洲精品美女久久av网站| 亚洲欧洲精品一区二区精品久久久 | 亚洲婷婷狠狠爱综合网| 亚洲第一av免费看| 女人精品久久久久毛片| 日本黄色片子视频| 寂寞人妻少妇视频99o| 一级a做视频免费观看| 久久精品国产亚洲av天美| 男女啪啪激烈高潮av片| 99久久综合免费| 大香蕉久久网| 亚洲av电影在线观看一区二区三区| 一边摸一边做爽爽视频免费| 18在线观看网站| 美女内射精品一级片tv| 亚洲,欧美,日韩| 国产精品一区二区在线不卡| 精品视频人人做人人爽| 国产成人免费观看mmmm| 极品少妇高潮喷水抽搐| 日韩中文字幕视频在线看片| 久久精品熟女亚洲av麻豆精品| 免费观看a级毛片全部| 欧美成人精品欧美一级黄| 制服丝袜香蕉在线| 久久韩国三级中文字幕| 少妇高潮的动态图| 狂野欧美激情性xxxx在线观看| 成人二区视频| 国产熟女午夜一区二区三区 | 蜜桃在线观看..| 国产亚洲午夜精品一区二区久久| 亚洲精品一区蜜桃| 日韩欧美一区视频在线观看| 男女高潮啪啪啪动态图| 天天影视国产精品| 中文欧美无线码| 丝袜在线中文字幕| 婷婷色麻豆天堂久久| 亚洲,一卡二卡三卡| 日本黄色日本黄色录像| 久久久精品94久久精品| 免费黄频网站在线观看国产| 国产精品国产三级国产av玫瑰| 纯流量卡能插随身wifi吗| 欧美精品一区二区大全| 国语对白做爰xxxⅹ性视频网站| 中国美白少妇内射xxxbb| 成人国语在线视频| 欧美 日韩 精品 国产| 国产精品成人在线| 亚洲五月色婷婷综合| 熟女电影av网| 精品久久久精品久久久| 国产又色又爽无遮挡免| 少妇丰满av| 精品亚洲乱码少妇综合久久| 国产 一区精品| 欧美精品亚洲一区二区| 日韩人妻高清精品专区| 亚洲av免费高清在线观看| 九草在线视频观看| 卡戴珊不雅视频在线播放| 色5月婷婷丁香| 我的女老师完整版在线观看| 中文字幕久久专区| 亚洲国产av影院在线观看| 亚洲图色成人| 中文欧美无线码| 欧美xxxx性猛交bbbb| 免费高清在线观看日韩| 啦啦啦视频在线资源免费观看| 午夜福利,免费看| 欧美bdsm另类| 国产日韩欧美视频二区| 国产女主播在线喷水免费视频网站| 国语对白做爰xxxⅹ性视频网站| 男女免费视频国产| 99九九在线精品视频| 精品久久久噜噜| 啦啦啦在线观看免费高清www| 建设人人有责人人尽责人人享有的| 国产有黄有色有爽视频| 欧美少妇被猛烈插入视频| 亚洲国产精品999| 一区在线观看完整版| 午夜激情福利司机影院| 免费观看av网站的网址| 成人无遮挡网站| 成年人免费黄色播放视频| √禁漫天堂资源中文www| 国产精品人妻久久久久久| 大片电影免费在线观看免费| 欧美老熟妇乱子伦牲交| 啦啦啦在线观看免费高清www| 夜夜骑夜夜射夜夜干| 成人免费观看视频高清| 最近中文字幕高清免费大全6| 亚洲国产精品999| 午夜久久久在线观看| 亚洲国产成人一精品久久久| 国产女主播在线喷水免费视频网站| 美女xxoo啪啪120秒动态图| 国产淫语在线视频| 亚洲av电影在线观看一区二区三区| 亚洲精品久久成人aⅴ小说 | 交换朋友夫妻互换小说| 亚洲人与动物交配视频| 亚洲av欧美aⅴ国产| 美女内射精品一级片tv| 天天操日日干夜夜撸| 一本大道久久a久久精品| 日韩一本色道免费dvd| 夜夜骑夜夜射夜夜干| 国产欧美亚洲国产| 一本大道久久a久久精品| 伊人久久精品亚洲午夜| av天堂久久9| 女性生殖器流出的白浆| 视频中文字幕在线观看| 久久国产精品大桥未久av| 成人国产av品久久久| 男的添女的下面高潮视频| 最近2019中文字幕mv第一页| 边亲边吃奶的免费视频| 18禁在线播放成人免费| 精品国产一区二区三区久久久樱花| 大片免费播放器 马上看| 国产成人精品无人区| av不卡在线播放| 欧美激情极品国产一区二区三区 | 亚洲精品第二区| 国产免费一级a男人的天堂| 热re99久久国产66热| 纵有疾风起免费观看全集完整版| 99久国产av精品国产电影| 欧美激情国产日韩精品一区| 一区二区三区免费毛片| 久久国产精品大桥未久av| 午夜激情av网站| 人成视频在线观看免费观看| 能在线免费看毛片的网站| 久久久亚洲精品成人影院| 老司机影院成人| 久久99精品国语久久久| a级毛色黄片| 国产熟女午夜一区二区三区 | 天堂中文最新版在线下载| 水蜜桃什么品种好| 日本黄大片高清| videossex国产| a级毛片在线看网站| 国产 精品1| 成人亚洲欧美一区二区av| 久久国产亚洲av麻豆专区| 99视频精品全部免费 在线| 午夜福利影视在线免费观看| 全区人妻精品视频| 在线免费观看不下载黄p国产| 久久热精品热| 久久午夜福利片| 一级,二级,三级黄色视频| 九色亚洲精品在线播放| 三级国产精品欧美在线观看| 国产69精品久久久久777片| 国产黄频视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 色婷婷av一区二区三区视频| 纯流量卡能插随身wifi吗| av福利片在线| 丰满饥渴人妻一区二区三| 桃花免费在线播放| 精品人妻一区二区三区麻豆| 国产亚洲午夜精品一区二区久久| 久久久欧美国产精品| 国产国拍精品亚洲av在线观看| 国产在线视频一区二区| 99九九线精品视频在线观看视频| 欧美激情国产日韩精品一区| 国产综合精华液| 亚洲久久久国产精品| 校园人妻丝袜中文字幕| 日韩,欧美,国产一区二区三区| 国产日韩一区二区三区精品不卡 | 国产乱人偷精品视频| 日韩av免费高清视频| 日韩欧美一区视频在线观看| 亚洲精品国产色婷婷电影| 女的被弄到高潮叫床怎么办| 少妇的逼好多水| 日本欧美视频一区| 777米奇影视久久| 狠狠精品人妻久久久久久综合| 肉色欧美久久久久久久蜜桃| 亚洲高清免费不卡视频| www.av在线官网国产| 黑人猛操日本美女一级片| 国产高清三级在线| a级毛片在线看网站| 男男h啪啪无遮挡| 久久久久久久久久久久大奶| 九色亚洲精品在线播放| 午夜日本视频在线| 一边摸一边做爽爽视频免费| av网站免费在线观看视频| 亚洲美女搞黄在线观看| 国产伦精品一区二区三区视频9| 99久久人妻综合| 色婷婷av一区二区三区视频| 黑人巨大精品欧美一区二区蜜桃 | 97超碰精品成人国产| av一本久久久久| 午夜福利视频在线观看免费| 少妇人妻久久综合中文| 少妇精品久久久久久久| 18禁裸乳无遮挡动漫免费视频| 免费黄色在线免费观看| 美女大奶头黄色视频| 久久久久久人妻| 免费播放大片免费观看视频在线观看| 高清毛片免费看| 国产成人91sexporn| 老熟女久久久| 欧美性感艳星| 国产成人一区二区在线| 国产综合精华液| 日韩,欧美,国产一区二区三区| 最近中文字幕高清免费大全6| av网站免费在线观看视频| 久久精品夜色国产| 视频在线观看一区二区三区| 久久久久国产精品人妻一区二区| 精品国产露脸久久av麻豆| 久久久久精品性色| 国产精品国产三级专区第一集| 日本黄色片子视频| 99热国产这里只有精品6| 国产精品人妻久久久久久| 国产精品久久久久成人av| 看非洲黑人一级黄片| 一级毛片aaaaaa免费看小| 精品少妇久久久久久888优播| 两个人的视频大全免费| 蜜桃久久精品国产亚洲av| 青春草亚洲视频在线观看| 免费观看av网站的网址| 三上悠亚av全集在线观看| 狂野欧美激情性xxxx在线观看| 日本-黄色视频高清免费观看| 欧美成人精品欧美一级黄| 一级毛片电影观看| 人体艺术视频欧美日本| 晚上一个人看的免费电影| 国产亚洲最大av| 你懂的网址亚洲精品在线观看| 18禁动态无遮挡网站| 欧美另类一区| 精品熟女少妇av免费看| 欧美日韩亚洲高清精品| 搡老乐熟女国产| 制服丝袜香蕉在线| 国产深夜福利视频在线观看| 国产成人av激情在线播放 |