• <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)
    亚洲少妇的诱惑av| www.精华液| 亚洲一区二区三区欧美精品| 国产欧美日韩精品亚洲av| 国产深夜福利视频在线观看| 欧美人与性动交α欧美软件| 国产单亲对白刺激| 亚洲人成电影观看| 一区二区日韩欧美中文字幕| 国产高清视频在线播放一区| 99国产精品99久久久久| 一级毛片电影观看| 中文字幕色久视频| 别揉我奶头~嗯~啊~动态视频| 国产亚洲精品久久久久5区| 久久毛片免费看一区二区三区| 另类精品久久| 9191精品国产免费久久| 国产高清videossex| 操美女的视频在线观看| 大型av网站在线播放| 人人妻人人澡人人爽人人夜夜| 人妻 亚洲 视频| www.自偷自拍.com| 天堂俺去俺来也www色官网| 国产免费视频播放在线视频| 欧美另类亚洲清纯唯美| 成人黄色视频免费在线看| 一级,二级,三级黄色视频| 一二三四社区在线视频社区8| 91成人精品电影| 久久久久久免费高清国产稀缺| 天天躁日日躁夜夜躁夜夜| www.精华液| av网站在线播放免费| 黄色片一级片一级黄色片| 纵有疾风起免费观看全集完整版| 国产欧美亚洲国产| 99re在线观看精品视频| 黄色 视频免费看| 欧美乱码精品一区二区三区| 久久久久久久久免费视频了| 狂野欧美激情性xxxx| 国产精品偷伦视频观看了| 亚洲av第一区精品v没综合| 午夜福利在线免费观看网站| 久久久欧美国产精品| 久久久精品免费免费高清| 国产真人三级小视频在线观看| 美女视频免费永久观看网站| 亚洲专区字幕在线| 超碰成人久久| 精品久久久精品久久久| 亚洲精品av麻豆狂野| 人人妻人人添人人爽欧美一区卜| 亚洲avbb在线观看| 中文字幕人妻熟女乱码| 大香蕉久久成人网| 精品福利永久在线观看| 成年人免费黄色播放视频| 久久久精品免费免费高清| 999久久久国产精品视频| 一级a爱视频在线免费观看| 黄色成人免费大全| 热99国产精品久久久久久7| 黑人操中国人逼视频| 国产人伦9x9x在线观看| 国产单亲对白刺激| 国产极品粉嫩免费观看在线| 汤姆久久久久久久影院中文字幕| 中亚洲国语对白在线视频| a级毛片在线看网站| 免费女性裸体啪啪无遮挡网站| 免费高清在线观看日韩| 久久久水蜜桃国产精品网| 久久亚洲精品不卡| 日韩视频在线欧美| 十八禁高潮呻吟视频| 自线自在国产av| 亚洲全国av大片| 男人舔女人的私密视频| 精品人妻在线不人妻| 亚洲自偷自拍图片 自拍| 黄片播放在线免费| 免费一级毛片在线播放高清视频 | 十八禁网站网址无遮挡| 日韩欧美一区二区三区在线观看 | 中文字幕最新亚洲高清| 18禁裸乳无遮挡动漫免费视频| 免费观看a级毛片全部| 一级片免费观看大全| 免费女性裸体啪啪无遮挡网站| 国产深夜福利视频在线观看| av在线播放免费不卡| 国产97色在线日韩免费| 老司机影院毛片| 久久人人97超碰香蕉20202| 亚洲午夜精品一区,二区,三区| 国产在视频线精品| 天天躁日日躁夜夜躁夜夜| 又黄又粗又硬又大视频| 女人爽到高潮嗷嗷叫在线视频| 大片免费播放器 马上看| 亚洲自偷自拍图片 自拍| 国产精品久久久久久精品电影小说| 亚洲午夜精品一区,二区,三区| 久久国产精品影院| 精品国产亚洲在线| 欧美日韩黄片免| 一级黄色大片毛片| 侵犯人妻中文字幕一二三四区| 欧美日韩av久久| 午夜激情av网站| 成年人午夜在线观看视频| 老熟妇乱子伦视频在线观看| 不卡av一区二区三区| 欧美精品啪啪一区二区三区| 国产精品熟女久久久久浪| 又黄又粗又硬又大视频| 丰满少妇做爰视频| 天天躁日日躁夜夜躁夜夜| 国产成人精品久久二区二区免费| 日日摸夜夜添夜夜添小说| 久久久久久久久免费视频了| 亚洲欧美日韩高清在线视频 | 国产亚洲精品第一综合不卡| 丁香欧美五月| 国产真人三级小视频在线观看| av在线播放免费不卡| 欧美老熟妇乱子伦牲交| 99在线人妻在线中文字幕 | 亚洲视频免费观看视频| 日韩欧美三级三区| 亚洲中文av在线| videosex国产| 国产有黄有色有爽视频| 成在线人永久免费视频| 日韩欧美免费精品| 免费观看av网站的网址| 久久久久久久精品吃奶| av网站免费在线观看视频| 天天操日日干夜夜撸| 亚洲国产欧美一区二区综合| 色精品久久人妻99蜜桃| 成年人黄色毛片网站| tocl精华| 色婷婷av一区二区三区视频| 一区二区av电影网| 国产精品免费大片| 正在播放国产对白刺激| 精品第一国产精品| 国产aⅴ精品一区二区三区波| 亚洲精品久久午夜乱码| 久久 成人 亚洲| 我要看黄色一级片免费的| 又紧又爽又黄一区二区| 久久久久久免费高清国产稀缺| 12—13女人毛片做爰片一| 嫩草影视91久久| 亚洲伊人色综图| 国产又色又爽无遮挡免费看| 一本一本久久a久久精品综合妖精| 久久精品人人爽人人爽视色| 国产免费福利视频在线观看| 男女边摸边吃奶| tocl精华| 午夜激情久久久久久久| 精品一区二区三区av网在线观看 | 日本黄色视频三级网站网址 | 丝袜在线中文字幕| 男女午夜视频在线观看| 99国产精品免费福利视频| www.999成人在线观看| 久久久国产精品麻豆| av超薄肉色丝袜交足视频| 日韩有码中文字幕| 精品一区二区三区四区五区乱码| 首页视频小说图片口味搜索| 成人18禁在线播放| 亚洲精品av麻豆狂野| 在线十欧美十亚洲十日本专区| 久久国产精品影院| 美女视频免费永久观看网站| 精品亚洲成a人片在线观看| 80岁老熟妇乱子伦牲交| 女人久久www免费人成看片| 久久中文字幕人妻熟女| 国产一区二区三区综合在线观看| 男女午夜视频在线观看| 精品国产超薄肉色丝袜足j| 一本综合久久免费| 高清在线国产一区| 中文字幕色久视频| 中国美女看黄片| 久久精品熟女亚洲av麻豆精品| 999精品在线视频| 99九九在线精品视频| 亚洲国产中文字幕在线视频| 免费人妻精品一区二区三区视频| 欧美激情极品国产一区二区三区| 亚洲美女黄片视频| 欧美国产精品va在线观看不卡| 国产人伦9x9x在线观看| 中文欧美无线码| 日本a在线网址| 99re6热这里在线精品视频| 少妇精品久久久久久久| 久久精品亚洲熟妇少妇任你| 国产aⅴ精品一区二区三区波| 精品少妇久久久久久888优播| 天天躁日日躁夜夜躁夜夜| 变态另类成人亚洲欧美熟女 | 亚洲精品一卡2卡三卡4卡5卡| 超碰97精品在线观看| 午夜福利视频在线观看免费| 欧美激情 高清一区二区三区| 91av网站免费观看| 国产色视频综合| 国产成人av教育| 成人18禁在线播放| 男男h啪啪无遮挡| 少妇 在线观看| 97在线人人人人妻| 久久毛片免费看一区二区三区| 国产精品自产拍在线观看55亚洲 | 母亲3免费完整高清在线观看| 欧美精品一区二区免费开放| 日韩中文字幕视频在线看片| 欧美大码av| 欧美人与性动交α欧美软件| 欧美 日韩 精品 国产| 一区福利在线观看| 大香蕉久久网| 啪啪无遮挡十八禁网站| 久久这里只有精品19| 中文字幕最新亚洲高清| 一本久久精品| 国产一区二区三区视频了| 少妇精品久久久久久久| 亚洲男人天堂网一区| 91av网站免费观看| 国产熟女午夜一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 欧美乱码精品一区二区三区| 在线观看www视频免费| 巨乳人妻的诱惑在线观看| 精品少妇黑人巨大在线播放| 亚洲av第一区精品v没综合| 欧美激情极品国产一区二区三区| 久久av网站| 久久中文看片网| 久久毛片免费看一区二区三区| av视频免费观看在线观看| avwww免费| 99re在线观看精品视频| tube8黄色片| 嫩草影视91久久| 精品国产亚洲在线| 久久久精品区二区三区| 天堂俺去俺来也www色官网| 在线观看免费视频日本深夜| 老司机在亚洲福利影院| 国产精品国产高清国产av | 激情在线观看视频在线高清 | 丝瓜视频免费看黄片| kizo精华| 久久亚洲精品不卡| 90打野战视频偷拍视频| 日本撒尿小便嘘嘘汇集6| 久久精品亚洲精品国产色婷小说| 国产日韩欧美在线精品| 亚洲七黄色美女视频| 老司机午夜十八禁免费视频| 国产1区2区3区精品| 少妇 在线观看| 成年动漫av网址| 亚洲国产欧美一区二区综合| svipshipincom国产片| 女性被躁到高潮视频| av欧美777| 国产福利在线免费观看视频| 精品午夜福利视频在线观看一区 | 亚洲视频免费观看视频| 欧美精品一区二区免费开放| 91精品国产国语对白视频| 91成人精品电影| 亚洲黑人精品在线| 高潮久久久久久久久久久不卡| 日韩精品免费视频一区二区三区| 亚洲免费av在线视频| 在线亚洲精品国产二区图片欧美| 十八禁人妻一区二区| 色婷婷久久久亚洲欧美| 天堂中文最新版在线下载| 在线观看免费高清a一片| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| av在线播放免费不卡| 国产97色在线日韩免费| 国产日韩欧美亚洲二区| 在线观看免费午夜福利视频| 午夜精品国产一区二区电影| 别揉我奶头~嗯~啊~动态视频| 一边摸一边做爽爽视频免费| 亚洲一码二码三码区别大吗| 色婷婷久久久亚洲欧美| 狠狠狠狠99中文字幕| 大片免费播放器 马上看| 国产一区二区三区综合在线观看| 一级毛片女人18水好多| 国产欧美日韩精品亚洲av| 国产精品一区二区在线不卡| 国产精品国产av在线观看| 少妇的丰满在线观看| 久久精品国产a三级三级三级| 久久香蕉激情| 中文字幕制服av| 成人手机av| 水蜜桃什么品种好| 亚洲一码二码三码区别大吗| a级毛片黄视频| 久久人妻熟女aⅴ| 在线观看舔阴道视频| 亚洲精品国产精品久久久不卡| 最新美女视频免费是黄的| 男女无遮挡免费网站观看| 老司机午夜福利在线观看视频 | 波多野结衣一区麻豆| 女同久久另类99精品国产91| 桃红色精品国产亚洲av| 日本黄色视频三级网站网址 | 欧美久久黑人一区二区| 高清视频免费观看一区二区| 午夜免费鲁丝| 在线天堂中文资源库| 人人妻人人添人人爽欧美一区卜| 精品人妻熟女毛片av久久网站| 中文字幕最新亚洲高清| 免费在线观看完整版高清| 精品欧美一区二区三区在线| 日韩一区二区三区影片| 国产欧美日韩一区二区三| 下体分泌物呈黄色| 露出奶头的视频| 精品国产一区二区三区久久久樱花| 亚洲伊人久久精品综合| 精品国产超薄肉色丝袜足j| 男女下面插进去视频免费观看| 9色porny在线观看| 国产日韩欧美在线精品| 日本一区二区免费在线视频| 九色亚洲精品在线播放| 免费女性裸体啪啪无遮挡网站| 成人18禁高潮啪啪吃奶动态图| 久久人人97超碰香蕉20202| 久久天躁狠狠躁夜夜2o2o| 欧美日韩国产mv在线观看视频| 考比视频在线观看| 精品久久久久久电影网| 欧美激情久久久久久爽电影 | 视频区欧美日本亚洲| 久久久国产一区二区| 一边摸一边做爽爽视频免费| 免费高清在线观看日韩| 高清在线国产一区| 国产精品影院久久| 日韩大码丰满熟妇| 香蕉久久夜色| www.999成人在线观看| 色综合欧美亚洲国产小说| 国产成人欧美| 电影成人av| 热re99久久国产66热| 19禁男女啪啪无遮挡网站| 国精品久久久久久国模美| 黄色视频不卡| 99久久99久久久精品蜜桃| 欧美日韩中文字幕国产精品一区二区三区 | 精品国产超薄肉色丝袜足j| 91国产中文字幕| 超碰97精品在线观看| 国产欧美日韩一区二区精品| 亚洲专区中文字幕在线| av片东京热男人的天堂| 国产精品麻豆人妻色哟哟久久| 成人国产av品久久久| 久久久精品国产亚洲av高清涩受| 丁香六月天网| 国产男女内射视频| 亚洲视频免费观看视频| 精品国产乱子伦一区二区三区| 99在线人妻在线中文字幕 | 欧美+亚洲+日韩+国产| 99久久国产精品久久久| 波多野结衣av一区二区av| 国产成人欧美在线观看 | av网站免费在线观看视频| av免费在线观看网站| 日韩大码丰满熟妇| 亚洲国产欧美日韩在线播放| 久久这里只有精品19| 精品亚洲成a人片在线观看| 国产单亲对白刺激| 国产高清激情床上av| svipshipincom国产片| 亚洲视频免费观看视频| 黄色a级毛片大全视频| 国产欧美日韩精品亚洲av| 丁香六月欧美| 久久中文字幕人妻熟女| 曰老女人黄片| 国产精品成人在线| 伊人久久大香线蕉亚洲五| 免费不卡黄色视频| 精品国产一区二区三区四区第35| 高清在线国产一区| 国产一区二区三区综合在线观看| 亚洲av片天天在线观看| 丝袜美腿诱惑在线| 精品福利观看| 制服诱惑二区| 人人妻,人人澡人人爽秒播| 操美女的视频在线观看| 国产精品香港三级国产av潘金莲| 久久人妻av系列| 精品人妻熟女毛片av久久网站| 国产免费视频播放在线视频| 狠狠婷婷综合久久久久久88av| 国产老妇伦熟女老妇高清| 国产精品 国内视频| 一级黄色大片毛片| 美女主播在线视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成伊人成综合网2020| 亚洲国产看品久久| 成人影院久久| 国产午夜精品久久久久久| netflix在线观看网站| 热99re8久久精品国产| 亚洲 国产 在线| 亚洲成a人片在线一区二区| 无遮挡黄片免费观看| 欧美一级毛片孕妇| 一进一出好大好爽视频| 亚洲精品av麻豆狂野| 成在线人永久免费视频| 免费少妇av软件| 精品一区二区三区四区五区乱码| 午夜福利欧美成人| 最新美女视频免费是黄的| 欧美激情极品国产一区二区三区| 夜夜骑夜夜射夜夜干| 国产精品久久久久久人妻精品电影 | 国产免费现黄频在线看| 一区二区三区国产精品乱码| 在线看a的网站| 在线观看66精品国产| 后天国语完整版免费观看| 久久人人97超碰香蕉20202| 亚洲精品乱久久久久久| 91麻豆精品激情在线观看国产 | 搡老乐熟女国产| 啦啦啦 在线观看视频| av电影中文网址| 一级片免费观看大全| 亚洲精品国产区一区二| 亚洲免费av在线视频| 久热这里只有精品99| 亚洲情色 制服丝袜| 色94色欧美一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利乱码中文字幕| 精品国产乱子伦一区二区三区| 国产一区二区三区在线臀色熟女 | 欧美日韩国产mv在线观看视频| 男女下面插进去视频免费观看| 一区二区三区激情视频| 日本一区二区免费在线视频| 午夜福利欧美成人| 国产成人av激情在线播放| 国产精品.久久久| 国产精品熟女久久久久浪| a在线观看视频网站| 18禁观看日本| 美女扒开内裤让男人捅视频| 十八禁网站免费在线| 99久久99久久久精品蜜桃| 交换朋友夫妻互换小说| 久久久久久免费高清国产稀缺| aaaaa片日本免费| 一区二区日韩欧美中文字幕| 亚洲精品av麻豆狂野| 男女之事视频高清在线观看| 黄网站色视频无遮挡免费观看| 黄色片一级片一级黄色片| 国产日韩欧美在线精品| 精品一区二区三区av网在线观看 | 中文字幕av电影在线播放| 亚洲av第一区精品v没综合| 天堂俺去俺来也www色官网| 亚洲国产av新网站| 国产av一区二区精品久久| 日日摸夜夜添夜夜添小说| 狠狠狠狠99中文字幕| 在线观看www视频免费| 免费不卡黄色视频| 蜜桃在线观看..| 黄色毛片三级朝国网站| 男女下面插进去视频免费观看| 亚洲精品中文字幕在线视频| 18禁裸乳无遮挡动漫免费视频| 宅男免费午夜| 欧美日韩国产mv在线观看视频| 日韩人妻精品一区2区三区| 蜜桃国产av成人99| 日韩中文字幕视频在线看片| 亚洲久久久国产精品| 我要看黄色一级片免费的| 亚洲五月婷婷丁香| 国产亚洲精品第一综合不卡| 丰满迷人的少妇在线观看| 黑丝袜美女国产一区| 又紧又爽又黄一区二区| 精品久久久久久电影网| 一进一出好大好爽视频| 欧美中文综合在线视频| 老汉色av国产亚洲站长工具| 久久午夜亚洲精品久久| 日本a在线网址| 亚洲欧洲精品一区二区精品久久久| 国产精品久久久久成人av| 精品人妻熟女毛片av久久网站| 国产亚洲午夜精品一区二区久久| 亚洲成人免费av在线播放| 成人国产av品久久久| 国产免费福利视频在线观看| 亚洲午夜理论影院| 深夜精品福利| www.精华液| 亚洲成人免费av在线播放| 极品人妻少妇av视频| 亚洲午夜精品一区,二区,三区| 精品人妻熟女毛片av久久网站| 精品卡一卡二卡四卡免费| 久久国产精品影院| 露出奶头的视频| 伦理电影免费视频| 大型黄色视频在线免费观看| 国产一区二区三区在线臀色熟女 | 69精品国产乱码久久久| 丝袜喷水一区| 青青草视频在线视频观看| 丝袜美腿诱惑在线| 国产亚洲av高清不卡| 亚洲国产av影院在线观看| 国产免费av片在线观看野外av| 精品高清国产在线一区| 母亲3免费完整高清在线观看| 亚洲国产成人一精品久久久| 熟女少妇亚洲综合色aaa.| 国产欧美亚洲国产| 手机成人av网站| 最黄视频免费看| 久久精品熟女亚洲av麻豆精品| 啦啦啦在线免费观看视频4| 久久久久久久精品吃奶| 午夜精品国产一区二区电影| 国产男女内射视频| 久久久国产精品麻豆| 97人妻天天添夜夜摸| 国产极品粉嫩免费观看在线| 两个人免费观看高清视频| 真人做人爱边吃奶动态| 久久亚洲真实| 丝瓜视频免费看黄片| av一本久久久久| 黄色视频,在线免费观看| 欧美成人午夜精品| 日韩一卡2卡3卡4卡2021年| 久久天堂一区二区三区四区| 香蕉久久夜色| 中文欧美无线码| 丁香欧美五月| 一级毛片精品| 一级,二级,三级黄色视频| 极品少妇高潮喷水抽搐| 久久人人爽av亚洲精品天堂| 亚洲av日韩在线播放| 国产欧美日韩一区二区三区在线| 成人18禁在线播放| 国产av精品麻豆| 日本黄色视频三级网站网址 | 亚洲va日本ⅴa欧美va伊人久久| 欧美精品一区二区免费开放| 国产一区二区激情短视频| 变态另类成人亚洲欧美熟女 | 咕卡用的链子| 久久人妻熟女aⅴ| 水蜜桃什么品种好| 欧美乱妇无乱码| 黑人巨大精品欧美一区二区蜜桃| 国产成人欧美在线观看 | 99精品欧美一区二区三区四区| 一级毛片女人18水好多| 一个人免费在线观看的高清视频| 高清欧美精品videossex| 纵有疾风起免费观看全集完整版| 一边摸一边做爽爽视频免费| 国产国语露脸激情在线看| 国产亚洲精品第一综合不卡| 国产成人啪精品午夜网站| 97在线人人人人妻| bbb黄色大片| 国产深夜福利视频在线观看| 久久久精品国产亚洲av高清涩受|