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

    新安江模型與水箱模型在馬家溝水庫洪水預(yù)報(bào)的適用性研究

    2023-08-02 14:52:54馬靖航咸永財(cái)何學(xué)平韓牧原杜佰林阮炳南許鎏佳
    人民珠江 2023年7期
    關(guān)鍵詞:模型

    馬靖航,咸永財(cái),何學(xué)平,劉 明,韓牧原,杜佰林,阮炳南,許鎏佳,吳 磊*

    (1.國能陜西水電有限公司,陜西 西安 710000;2.西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)試驗(yàn)室,陜西 楊凌 712100;3.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院,陜西 楊凌 712100)

    洪水預(yù)報(bào)是直接為國民經(jīng)濟(jì)建設(shè)服務(wù)的重要基礎(chǔ)性工作[1-3]。它是根據(jù)洪水形成和運(yùn)動(dòng)規(guī)律,利用水文、氣象信息,預(yù)測洪水的發(fā)生與變化過程的應(yīng)用基礎(chǔ)性技術(shù)科學(xué)[4]。準(zhǔn)確的洪水預(yù)報(bào)是防汛搶險(xiǎn)和防洪系統(tǒng)調(diào)度運(yùn)用的決策依據(jù)[5-6]。同時(shí),也為減少洪災(zāi)損失爭取時(shí)間、有助于更好地控制和利用洪水資源,是重要的防洪減災(zāi)非工程措施,對于優(yōu)化水庫調(diào)度、確保防洪安全具有重要意義[7-8]。

    通常,把預(yù)見期在2 d以內(nèi)的洪水預(yù)報(bào)稱短期預(yù)報(bào),洪峰水位、洪峰流量、洪現(xiàn)時(shí)間和一次洪水總量等是主要的預(yù)報(bào)項(xiàng)目。短期洪水預(yù)報(bào)多采用基于物理成因、有一定理論基礎(chǔ)的經(jīng)驗(yàn)性方法。這類方法的優(yōu)點(diǎn)是結(jié)構(gòu)簡單、使用方便,但缺乏物理意義,模擬精度不高[9]。隨著計(jì)算機(jī)水平的快速發(fā)展,流域水文模型已逐漸成為現(xiàn)代實(shí)時(shí)洪水預(yù)報(bào)系統(tǒng)的核心模塊[10-12]。國內(nèi)外相繼提出了Sacramento、TANK、新安江、HBV、TUWMODEL、HEC-HMS、TOPMODEL、SHE、VIC、流溪河模型等一系列集總式和分布式水文模型,以及多種模型參數(shù)優(yōu)化方法,洪水預(yù)報(bào)也隨之發(fā)展[13-17]。

    概念性水文模型相比黑箱模型以及分布式水文模型,由于其數(shù)據(jù)需求量相對較小而常用于歷史資料欠缺的地區(qū)[18]。新安江模型與水箱模型是概念性水文模型的典型代表,它們以物理成因?yàn)榛A(chǔ),采用概化和假設(shè)的方法對水文現(xiàn)象進(jìn)行數(shù)學(xué)模擬,模型的結(jié)構(gòu)和參數(shù)具有一定物理意義。新安江模型由河海大學(xué)趙人俊教授團(tuán)隊(duì)提出,模型最初為二水源新安江模型(地表徑流和地下徑流),后引入薩克拉門托模型與水箱模型中線性水庫函數(shù)劃分水源的概念,提出了三水源新安江模型(地表徑流、壤中流、地下徑流),在中國濕潤和半濕潤地區(qū)應(yīng)用廣泛,取得了較好的效果,其主要特點(diǎn)是蓄滿產(chǎn)流結(jié)構(gòu)[19-20]。水箱模型(Tank model)是1961年由日本的菅原正已博士構(gòu)建,經(jīng)后續(xù)研究者的不斷改進(jìn)和更新?lián)Q代逐漸成為一種被廣泛應(yīng)用的水文模型[21]。水箱模型能以水箱中蓄水量出流關(guān)系來模擬復(fù)雜的降雨徑流轉(zhuǎn)換過程,結(jié)構(gòu)靈活、適用性強(qiáng),在洪水預(yù)報(bào)、水沙研究等方面應(yīng)用良好[22-23]。

    湑水河干流自上而下分觀音峽水電站、金龍水電站、黑峽子水電站、八仙園一級水電站、八仙園二級水電站、馬家溝水電站等共 11 級開發(fā)。觀音峽水電站為第一級,控制流域面積431 km2,裝機(jī)容量2.6萬 kW,為引水式水電站,年均入庫徑流量19 731 萬m3。馬家溝水電站位于湑水河中游,為第六級,壩址位于石槽河匯合口上游約3.0 km處,控制流域面積1 311 km2,具有防洪、發(fā)電等功能。電站廠房位于石槽河匯合口上游處,流域面積1 323 km2,屬壩-引水混合式水電站。

    本文以馬家溝水庫為對象,分別選用新安江模型和水箱模型模擬馬家溝水庫流域近3 a的日徑流過程與次洪過程,采用遺傳算法(GA)率定模型參數(shù),從模擬效果、指標(biāo)誤差等方面對比分析新安江和TANK模型在水文模擬過程中的適應(yīng)性,綜合評估2種洪水預(yù)報(bào)方法的可靠性。研究結(jié)果對及時(shí)發(fā)出洪水災(zāi)害預(yù)警、制定防洪措施、減少洪災(zāi)損失具有重要意義。

    1 研究方法

    1.1 流域概況

    馬家溝水庫位于陜西省城固縣,地處漢江一級支流湑水河中游,多年平均徑流量6.43億m3,水庫正常蓄水位794 m,正常蓄水位以下庫容2 706萬m3。流域?qū)俦眮啛釒駶櫦撅L(fēng)氣候,雨量年際變化大,年內(nèi)分配不均;地形以秦嶺山地丘陵為主,地勢高、起伏變化大、河谷縱橫交錯(cuò),河系支流呈羽型分布。流域全年溫暖濕潤,符合新安江模型、水箱模型的應(yīng)用條件。

    圖1 馬家溝水庫流域水系及水庫、雨量站位置分布

    1.2 新安江模型與水箱模型

    三水源新安江(XAJ)模型采用蓄滿產(chǎn)流概念,對降雨空間分布不均勻采用分單元產(chǎn)流計(jì)算,流域調(diào)蓄對坡面不同水源的差異采用分水源匯流計(jì)算,坡面匯流與河網(wǎng)匯流的差異采用分階段匯流計(jì)算,并引入了土壤蓄水量和自由水容量分布曲線[24-25]。模型整體包含蒸散發(fā)模塊、產(chǎn)流模塊、水源劃分模塊和匯流模塊[26-28]。

    水箱(TANK)模型是一種概念性徑流模型,能以簡單的形式模擬徑流的形成,并將流域產(chǎn)匯流中各個(gè)水文過程,利用多個(gè)彼此聯(lián)系的垂直串聯(lián)水箱進(jìn)行模擬[23]。除底層水箱沒有底孔外,一般而言每個(gè)水箱都由底孔和邊孔組成。其中邊孔出流表示不同的徑流產(chǎn)流,底孔出流代表土壤下滲,同時(shí)也是下層水箱的入流,因此一個(gè)簡單水箱通常有3個(gè)主要參數(shù),分別為邊孔高度,水箱邊孔出流系數(shù)和水箱底孔下滲系數(shù)[25]。

    1.3 數(shù)據(jù)來源與處理

    研究所使用的氣象、流量數(shù)據(jù)均來源于國能城固馬家溝水電有限公司。在考慮峰高量大、降雨徑流關(guān)系顯著原則的基礎(chǔ)上,根據(jù)搜集的2019—2021年水庫逐日流量、水位和泄洪資料,整理得到9場歷史洪水?dāng)?shù)據(jù),每年各3場。域內(nèi)量級和強(qiáng)度較大的暴雨一般出現(xiàn)在7—9月,洪水由暴雨形成,出現(xiàn)時(shí)間與暴雨基本一致,受暴雨特性的影響,洪水多為單式洪峰,受秋淋連陰雨影響也易出現(xiàn)復(fù)式洪峰。采用泰森多邊形方法[29],計(jì)算黃柏塬等12個(gè)雨量站的面平均雨量,見式(1):

    (1)

    式中P——面降水量;Pi——站點(diǎn)降水量;wi——站點(diǎn)面積權(quán)重。

    經(jīng)降水量、降水歷時(shí)、降水強(qiáng)度、洪量、洪峰流量、洪水歷時(shí)、徑流深、徑流系數(shù)等指標(biāo)的聚類分析,發(fā)現(xiàn)9場洪水具有一定的代表性?;诖?設(shè)定2019—2020年為日模率定期, 2021年為日模檢驗(yàn)期;選用2019—2020年的6場代表性洪水用于次洪率定, 2021年的3場洪水用于次洪檢驗(yàn)。

    1.4 遺傳算法與實(shí)時(shí)校正

    遺傳算法是Holland等于1975年模擬生物進(jìn)化規(guī)律提出的一種隨機(jī)、高度并行、自適應(yīng)的搜索方法。與傳統(tǒng)優(yōu)化算法不同,遺傳算法可直接對結(jié)構(gòu)對象進(jìn)行操作并自適應(yīng)地調(diào)整搜索方向,不需要明確的規(guī)則。本文以MATLAB R2018b為開發(fā)工具編寫遺傳算法程序,優(yōu)選模型參數(shù),設(shè)定種群規(guī)模為1 000,染色體編碼長度為50,交叉概率為0.75,變異概率為0.05,最大迭代次數(shù)為1 000,經(jīng)多次試算求解以實(shí)現(xiàn)模型參數(shù)的優(yōu)化。

    為減少洪水預(yù)報(bào)的誤差來源,在場次洪水過程模擬中引入了具有3 h預(yù)報(bào)時(shí)段的一階自回歸(AR(1))實(shí)時(shí)校正方法[30]。根據(jù)誤差相依的特性,首先計(jì)算觀測數(shù)據(jù)與模擬數(shù)據(jù)的誤差序列,然后建立自回歸模型并依據(jù)最小二乘法進(jìn)行參數(shù)估計(jì),最后計(jì)算校正后的預(yù)報(bào)值,獲取預(yù)報(bào)結(jié)果。

    1.5 預(yù)報(bào)精度評價(jià)指標(biāo)與評定依據(jù)

    洪水預(yù)報(bào)精度評定包括洪峰流量(水位)、洪峰出現(xiàn)時(shí)間、洪量(徑流量)和洪水過程等,依據(jù)GB/T 22482—2008《水文情報(bào)預(yù)報(bào)規(guī)范》[31],選用洪量預(yù)報(bào)許可誤差FR,洪峰預(yù)報(bào)許可誤差FP,峰現(xiàn)時(shí)間預(yù)報(bào)許可誤差FT和Nash-Sutcliffe (NS)效率系數(shù)作為評價(jià)指標(biāo)對模擬結(jié)果進(jìn)行精度評定。

    1990年春,見廉小花對兒子特別好,徐河開始向廉小花套近乎。廉小花覺得徐河肯干,心地也不壞,自己和小云天又天生有緣,于是廉小花不顧家人的反對,于1990年冬天,執(zhí)意嫁給了徐河。

    2 結(jié)果與分析

    2.1 日徑流過程模擬分析

    2.1.1日模參數(shù)率定

    采用遺傳算法對新安江模型的16個(gè)參數(shù)進(jìn)行多次反復(fù)調(diào)試率定。首先隨機(jī)生成初始群體,計(jì)算群體中每個(gè)個(gè)體的適應(yīng)度,依次進(jìn)行選擇、交叉、變異運(yùn)算,當(dāng)進(jìn)化代數(shù)等于最大進(jìn)化代數(shù)時(shí),將最大適應(yīng)度個(gè)體作為最優(yōu)解輸出并停止優(yōu)化。采用試錯(cuò)法和遺傳算法對水箱模型日模參數(shù)進(jìn)行設(shè)計(jì)/率定,首先采用4個(gè)邊孔的單水箱模型對流域日徑流過程進(jìn)行模擬,確定出水箱模型的10個(gè)參數(shù)初值范圍,然后分別調(diào)整水箱中上孔系數(shù)和下孔出流系數(shù),并進(jìn)行自動(dòng)參數(shù)調(diào)試率定。統(tǒng)計(jì)新安江模型和水箱模型日模參數(shù)的率定結(jié)果,見表1。

    表1 新安江模型和水箱模型的日模優(yōu)化參數(shù)

    2.1.2日模結(jié)果分析

    新安江模型和水箱模型率定期、驗(yàn)證期日模精度評定結(jié)果,見表2。

    表2 新安江模型和水箱模型日值模擬精度評定結(jié)果

    通過表2可以看出,新安江模型和水箱模型在馬家溝水庫流域均具有較好的適應(yīng)性,對于旱期及汛期均具有較好的模擬效果,二者在率定期與驗(yàn)證期的各項(xiàng)洪水指標(biāo)基本能滿足洪水預(yù)報(bào)精度要求。為更直觀地對比新安江模型與水箱模型在馬家溝水庫流域的適用性,分別將二者的日模擬結(jié)果與觀測值進(jìn)行對比,見圖2。

    圖2 2019—2021年馬家溝水庫日徑流模擬結(jié)果對比

    新安江模型在率定期對洪水模擬效果較好,僅在2020年出現(xiàn)模擬洪峰過低現(xiàn)象(圖2),這可能是因?yàn)榱饔蚓植看蟊┯甓饔蛎娼涤昶?dǎo)致與實(shí)測徑流值偏差較大。在驗(yàn)證期新安江模型對模擬洪水存在高估現(xiàn)象,部分模擬洪峰在實(shí)際觀測中峰量較低或不存在,分析可能是馬家溝上游部分水庫攔蓄洪水所致。水箱模型對于較大洪水能進(jìn)行準(zhǔn)確預(yù)報(bào),其模擬偏差主要源于200 m3/s以下的小洪水,此類模擬誤差可能是因?yàn)榱饔蛏嫌蔚挠^音峽等水電站串聯(lián)調(diào)蓄對小型洪水產(chǎn)生了削峰作用。整體來看,2種模型在汛期的洪峰及洪量基本滿足洪水預(yù)報(bào)精度要求,二者均可用于馬家溝水庫流域的水文過程預(yù)報(bào)。

    2.2 次洪模擬過程比較

    2.2.1次洪參數(shù)率定

    新安江模型在次洪模擬的參數(shù)調(diào)試中,蒸散發(fā)參數(shù)K、WUM、WUM、WLM、C和產(chǎn)流參數(shù)IMP、B的取值參考了日值模型模擬的結(jié)果,水源劃分參數(shù)SM、KG、KSS、EX和匯流參數(shù)KKG、KKSS、CS、CR、L則采用遺傳算法率定優(yōu)選得出。水箱模型在次洪模擬的參數(shù)調(diào)試中,采用了具有3、2、1個(gè)邊孔的三水箱模型,考慮到洪水開始前均處于基流階段,將模型初始蓄水量設(shè)置為0;同時(shí)忽略汛期暴雨時(shí)段的蒸發(fā)和下滲作用,以消除蒸散發(fā)觀測數(shù)據(jù)缺乏的影響。同時(shí),引入實(shí)時(shí)校正方法,依據(jù)《水文預(yù)報(bào)情報(bào)規(guī)范》得到新安江模型和水箱模型模擬馬家溝水庫流域次洪過程的最優(yōu)參數(shù),見表3。

    表3 新安江模型和水箱模型次洪模擬的優(yōu)化參數(shù)

    2.2.2次洪模擬結(jié)果分析

    引入實(shí)時(shí)校正后,新安江模型和水箱模型率定期、驗(yàn)證期次洪過程模擬精度評定結(jié)果,見表4。

    表4 新安江模型和水箱模型次洪率定期與驗(yàn)證期模擬結(jié)果

    新安江模型和水箱模型均能較好地模擬次洪過程,洪量和洪峰均能基本滿足水文預(yù)報(bào)要求(表4)。其中,率定期間,新安江模型洪量相對誤差最大值為20%,平均相對誤差為13%;洪峰相對誤差最大值為24%,平均相對誤差為13%;NS效率系數(shù)的最大、最小值分別為0.87、0.43,平均值為0.55。水箱模型洪量相對誤差最大值為10%,平均相對誤差為1%;洪峰相對誤差最大值為16.9%,平均相對誤差為13%;NS效率系數(shù)的最大、最小值分別為0.86、0.40,平均值為0.44。驗(yàn)證期間,新安江模型洪量相對誤差最大值為14%,平均相對誤差為12%;洪峰相對誤差最大值為16%,平均相對誤差為11%;NS效率系數(shù)的最大、最小值分別為0.82、0.51,平均值為0.62。水箱模型洪量相對誤差最大值為7%,平均相對誤差為3%;洪峰相對誤差最大值為25%,平均相對誤差為10%;NS效率系數(shù)的最大、最小值分別為0.81、0.46,平均值為0.59。模型峰現(xiàn)時(shí)差均在5 h以內(nèi),這與模型設(shè)置時(shí)引入3 h的實(shí)時(shí)校正有關(guān),整體來看,9場洪水均能基本滿足預(yù)報(bào)合格標(biāo)準(zhǔn)。

    為進(jìn)一步描述2種模型在馬家溝水庫流域次洪模擬中的適用性,選擇4場洪水表征其洪水過程,洪號分別為20190722、20200803、20210818、20210918,其中20190722、20210818是單式洪峰,過程陡漲陡落;20200803、20210918是受秋淋連陰雨影響的復(fù)式洪峰,模擬洪水過程見圖 3。

    新安江模型和水箱模型均能較好地模擬場次洪水過程(圖3)。新安江模型模擬的次洪過程相較水箱模型更平滑,但水箱模型更貼近實(shí)測流量過程。水箱模型模擬四場洪水的整體效果更好,洪量、洪峰甚至峰現(xiàn)時(shí)間的表現(xiàn)均優(yōu)于新安江模型,其中模擬20190722、20210818(單式洪峰)的洪量、洪峰相對誤差以及峰現(xiàn)時(shí)間分別為0.3%、22%、4 h和3%、2%、3 h;模擬20200803、20210918(復(fù)式洪峰)的洪量、洪峰相對誤差以及峰現(xiàn)時(shí)間分別為4%、10%、2h和7%、3%、4 h;總體來看,單峰洪水的模擬效果優(yōu)于多峰洪水。

    a)20190722

    3 結(jié)論

    研究新安江模型和水箱模型在馬家溝水庫洪水預(yù)報(bào)中的適應(yīng)性,不僅對洪水預(yù)報(bào)模型選擇有一定的參考作用,而且對流域抗洪措施的制定、減少洪災(zāi)損失危害具有重要意義。主要結(jié)論如下。

    a)日模方面,水箱模型的模擬效果更優(yōu),洪量相對誤差,洪峰相對誤差、峰現(xiàn)時(shí)差和NS效率系數(shù)等4項(xiàng)指標(biāo)均達(dá)到了水文情報(bào)預(yù)報(bào)規(guī)范的評定精度要求;次模方面,新安江模型的模擬效果更平滑,水箱模型的洪量、洪峰以及NS效率系數(shù)的評定精度更高,更適合馬家溝水庫流域。

    b)集總概念式模型將流域作為一個(gè)整體單元進(jìn)行模擬,未充分考慮模型輸入數(shù)據(jù)的空間分散性和不均勻性,對流域徑流的預(yù)報(bào)會(huì)造成一定影響。盡管泰森多邊形方法可充分考慮降雨站點(diǎn)中心對雨量空間分布的影響,但采用泰森多邊形計(jì)算面降雨會(huì)坦化局部雨量站降雨峰值的影響,導(dǎo)致洪水預(yù)報(bào)結(jié)果出現(xiàn)一定偏差。

    c)為提高預(yù)報(bào)結(jié)果的準(zhǔn)確性和有效性,后續(xù)將強(qiáng)化以測補(bǔ)報(bào),根據(jù)前期汛情災(zāi)情檢視雨水情監(jiān)測站點(diǎn)布設(shè)是否滿足防御要求,滾動(dòng)修訂模型,并考慮引入分布式水文模型,不斷提高智能預(yù)報(bào)水平。

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    69av精品久久久久久| 亚洲色图 男人天堂 中文字幕| 午夜老司机福利片| 看免费av毛片| 欧美成人午夜精品| 久久精品aⅴ一区二区三区四区| 国产伦人伦偷精品视频| 国产97色在线日韩免费| 又黄又粗又硬又大视频| 亚洲美女黄片视频| 久久久国产精品麻豆| 男人操女人黄网站| 黄色a级毛片大全视频| 精品福利观看| 日韩欧美一区二区三区在线观看| av超薄肉色丝袜交足视频| 欧美成人免费av一区二区三区| 久久精品国产综合久久久| 亚洲国产欧美日韩在线播放| 美女高潮到喷水免费观看| 后天国语完整版免费观看| 一夜夜www| 老鸭窝网址在线观看| 久久草成人影院| 亚洲国产欧美网| 一区二区三区高清视频在线| 亚洲成av片中文字幕在线观看| 18禁裸乳无遮挡免费网站照片 | 午夜福利视频1000在线观看 | 国产黄a三级三级三级人| 午夜福利视频1000在线观看 | videosex国产| 亚洲成a人片在线一区二区| 亚洲中文日韩欧美视频| 国产野战对白在线观看| 久久中文字幕一级| 亚洲一区高清亚洲精品| 亚洲一区中文字幕在线| 久久狼人影院| 欧美成人一区二区免费高清观看 | АⅤ资源中文在线天堂| 久久香蕉激情| 精品国产国语对白av| 韩国精品一区二区三区| 制服丝袜大香蕉在线| bbb黄色大片| 亚洲第一电影网av| 91九色精品人成在线观看| 精品国产亚洲在线| 亚洲国产看品久久| 精品久久久久久久毛片微露脸| 国产欧美日韩一区二区精品| 午夜成年电影在线免费观看| 又黄又粗又硬又大视频| 此物有八面人人有两片| 日韩欧美三级三区| 夜夜爽天天搞| 亚洲人成77777在线视频| 免费看十八禁软件| 亚洲人成电影观看| 制服诱惑二区| 免费在线观看亚洲国产| 国产精品一区二区在线不卡| 一级毛片高清免费大全| 亚洲熟女毛片儿| 婷婷精品国产亚洲av在线| 久久人人精品亚洲av| 欧美老熟妇乱子伦牲交| 啪啪无遮挡十八禁网站| 中文字幕色久视频| 看免费av毛片| 麻豆一二三区av精品| 一级a爱视频在线免费观看| 精品一区二区三区视频在线观看免费| 女人高潮潮喷娇喘18禁视频| 男女做爰动态图高潮gif福利片 | 欧美成人一区二区免费高清观看 | 91老司机精品| 国产亚洲精品av在线| 免费无遮挡裸体视频| 少妇粗大呻吟视频| 日韩成人在线观看一区二区三区| 一进一出好大好爽视频| 久久精品人人爽人人爽视色| 一级片免费观看大全| 人成视频在线观看免费观看| 日韩欧美国产在线观看| 一区二区三区激情视频| 琪琪午夜伦伦电影理论片6080| 亚洲人成电影免费在线| 无限看片的www在线观看| 亚洲国产精品999在线| 久久精品国产亚洲av高清一级| 在线av久久热| 中文字幕精品免费在线观看视频| 一区福利在线观看| 久久久国产成人免费| 91精品国产国语对白视频| 欧美+亚洲+日韩+国产| or卡值多少钱| 欧美中文综合在线视频| 国产99白浆流出| 在线天堂中文资源库| 国产精品秋霞免费鲁丝片| 久久精品国产综合久久久| 成人三级黄色视频| 国产亚洲av嫩草精品影院| 夜夜躁狠狠躁天天躁| 纯流量卡能插随身wifi吗| or卡值多少钱| 国产片内射在线| 久久久久久免费高清国产稀缺| 久久精品影院6| 日本一区二区免费在线视频| 国产免费男女视频| 一进一出好大好爽视频| 一进一出抽搐动态| 9热在线视频观看99| 色综合亚洲欧美另类图片| 午夜福利视频1000在线观看 | 精品免费久久久久久久清纯| 欧美大码av| 每晚都被弄得嗷嗷叫到高潮| 国产高清有码在线观看视频 | 嫁个100分男人电影在线观看| 岛国在线观看网站| 国产真人三级小视频在线观看| 国产精品 国内视频| 午夜久久久久精精品| 久久婷婷成人综合色麻豆| 久久久久久大精品| 亚洲男人的天堂狠狠| 亚洲三区欧美一区| 免费看美女性在线毛片视频| 欧美大码av| 久久久久久人人人人人| 两个人视频免费观看高清| 亚洲欧洲精品一区二区精品久久久| 久久久久久免费高清国产稀缺| 色综合站精品国产| www国产在线视频色| 久热爱精品视频在线9| 91精品三级在线观看| 精品国产一区二区久久| av电影中文网址| 男女下面进入的视频免费午夜 | 国产麻豆69| 日韩国内少妇激情av| 久热这里只有精品99| 后天国语完整版免费观看| 久久精品国产综合久久久| 中文字幕人成人乱码亚洲影| 好男人在线观看高清免费视频 | 色在线成人网| 亚洲五月婷婷丁香| 国产免费男女视频| 久久精品国产清高在天天线| 最新在线观看一区二区三区| 夜夜看夜夜爽夜夜摸| 久久香蕉国产精品| 欧美成人免费av一区二区三区| 久久香蕉精品热| 黄色视频不卡| √禁漫天堂资源中文www| 黄色 视频免费看| 99国产综合亚洲精品| 欧美中文综合在线视频| 久久欧美精品欧美久久欧美| 在线观看午夜福利视频| 精品久久久精品久久久| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三区四区五区乱码| 女同久久另类99精品国产91| 国内毛片毛片毛片毛片毛片| 久久人人爽av亚洲精品天堂| 精品国产一区二区三区四区第35| 久99久视频精品免费| 亚洲自偷自拍图片 自拍| 久久中文看片网| 99久久国产精品久久久| 夜夜看夜夜爽夜夜摸| 91成年电影在线观看| 亚洲avbb在线观看| 12—13女人毛片做爰片一| 国产亚洲精品久久久久5区| 伦理电影免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 真人一进一出gif抽搐免费| 国产精品亚洲美女久久久| 亚洲视频免费观看视频| 精品久久蜜臀av无| 久久久久久人人人人人| 亚洲精品久久国产高清桃花| 日日夜夜操网爽| 好男人电影高清在线观看| 1024香蕉在线观看| 一夜夜www| 国产色视频综合| 9色porny在线观看| 日本黄色视频三级网站网址| 欧美av亚洲av综合av国产av| 午夜福利免费观看在线| 中文字幕色久视频| 在线观看66精品国产| 亚洲激情在线av| 欧美乱妇无乱码| 国产精品电影一区二区三区| 涩涩av久久男人的天堂| 亚洲五月天丁香| 国内久久婷婷六月综合欲色啪| 老熟妇乱子伦视频在线观看| 搡老岳熟女国产| 国产亚洲精品久久久久5区| 国产区一区二久久| 久久影院123| 欧美亚洲日本最大视频资源| 国产蜜桃级精品一区二区三区| 嫩草影院精品99| av福利片在线| 亚洲色图综合在线观看| 12—13女人毛片做爰片一| 国产成年人精品一区二区| 一个人观看的视频www高清免费观看 | 亚洲熟女毛片儿| 一区二区日韩欧美中文字幕| 久久精品国产亚洲av高清一级| 国产精品秋霞免费鲁丝片| 精品国产国语对白av| 欧美午夜高清在线| 成熟少妇高潮喷水视频| 老熟妇仑乱视频hdxx| 夜夜爽天天搞| 99国产精品一区二区三区| 亚洲熟妇中文字幕五十中出| 色尼玛亚洲综合影院| 神马国产精品三级电影在线观看 | 制服丝袜大香蕉在线| 两个人看的免费小视频| av欧美777| 国产成人一区二区三区免费视频网站| 身体一侧抽搐| a在线观看视频网站| 亚洲电影在线观看av| АⅤ资源中文在线天堂| 校园春色视频在线观看| 国产三级在线视频| 无人区码免费观看不卡| 国产亚洲av高清不卡| 欧美日韩乱码在线| 黑人操中国人逼视频| 欧美黄色片欧美黄色片| svipshipincom国产片| 99国产极品粉嫩在线观看| 国产亚洲av嫩草精品影院| 少妇熟女aⅴ在线视频| 99在线人妻在线中文字幕| 法律面前人人平等表现在哪些方面| 好男人电影高清在线观看| 亚洲精品久久国产高清桃花| 精品福利观看| 亚洲一码二码三码区别大吗| 亚洲一区二区三区不卡视频| 一个人免费在线观看的高清视频| 亚洲 欧美一区二区三区| 婷婷六月久久综合丁香| 老司机福利观看| 大型黄色视频在线免费观看| 国产av又大| 国产精品98久久久久久宅男小说| 国产高清有码在线观看视频 | 国产成人一区二区三区免费视频网站| 免费在线观看影片大全网站| 男人的好看免费观看在线视频 | 亚洲伊人色综图| 亚洲国产精品成人综合色| 久久精品国产综合久久久| 91精品三级在线观看| svipshipincom国产片| 久久久国产欧美日韩av| 非洲黑人性xxxx精品又粗又长| 亚洲欧美激情在线| 电影成人av| 成人永久免费在线观看视频| 嫩草影院精品99| 啦啦啦观看免费观看视频高清 | 国产精品野战在线观看| 99在线人妻在线中文字幕| 亚洲精品国产色婷婷电影| 亚洲伊人色综图| 99国产精品一区二区蜜桃av| 亚洲欧洲精品一区二区精品久久久| 999久久久精品免费观看国产| 国产成人精品在线电影| 可以在线观看毛片的网站| 国产高清激情床上av| 国产成人欧美在线观看| 无遮挡黄片免费观看| 国产极品粉嫩免费观看在线| 成人av一区二区三区在线看| 后天国语完整版免费观看| 在线观看免费视频日本深夜| 色精品久久人妻99蜜桃| 欧美日韩精品网址| cao死你这个sao货| 久久久久国产精品人妻aⅴ院| 黑丝袜美女国产一区| 看免费av毛片| 久久久水蜜桃国产精品网| 亚洲第一电影网av| 少妇被粗大的猛进出69影院| 法律面前人人平等表现在哪些方面| 亚洲国产精品999在线| 免费在线观看亚洲国产| 久久欧美精品欧美久久欧美| 精品无人区乱码1区二区| 麻豆久久精品国产亚洲av| 十八禁人妻一区二区| 精品福利观看| 亚洲精华国产精华精| 久久久久国产精品人妻aⅴ院| 叶爱在线成人免费视频播放| 一本综合久久免费| 香蕉丝袜av| 亚洲av片天天在线观看| 久久久久久大精品| 黑人操中国人逼视频| 国产精品日韩av在线免费观看 | 少妇裸体淫交视频免费看高清 | 性少妇av在线| 涩涩av久久男人的天堂| 亚洲一区高清亚洲精品| 国产成人系列免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 人人妻,人人澡人人爽秒播| 精品熟女少妇八av免费久了| 国产欧美日韩一区二区三区在线| 国产极品粉嫩免费观看在线| www.自偷自拍.com| 人人妻人人澡人人看| 12—13女人毛片做爰片一| 亚洲成av人片免费观看| 免费在线观看影片大全网站| 在线观看免费午夜福利视频| 久久人人97超碰香蕉20202| 久9热在线精品视频| 在线观看www视频免费| a在线观看视频网站| 麻豆成人av在线观看| 亚洲五月天丁香| 欧美日韩精品网址| 免费看十八禁软件| 久久精品91蜜桃| 十八禁人妻一区二区| 国产视频一区二区在线看| 欧美成人一区二区免费高清观看 | 别揉我奶头~嗯~啊~动态视频| 精品国产超薄肉色丝袜足j| 国产成人影院久久av| 老司机在亚洲福利影院| 中国美女看黄片| 变态另类成人亚洲欧美熟女 | 久久久久久人人人人人| 国产成人av教育| 91字幕亚洲| 日本一区二区免费在线视频| 国产精品久久久久久亚洲av鲁大| 免费不卡黄色视频| 搡老岳熟女国产| 婷婷丁香在线五月| 久久精品国产亚洲av香蕉五月| 美女扒开内裤让男人捅视频| 国产精品久久久久久精品电影 | 日韩中文字幕欧美一区二区| 久久午夜亚洲精品久久| 老熟妇仑乱视频hdxx| 日韩有码中文字幕| 日韩免费av在线播放| 国产精品美女特级片免费视频播放器 | 欧美日本中文国产一区发布| 窝窝影院91人妻| av有码第一页| 久久久国产成人免费| 999久久久精品免费观看国产| 亚洲成av人片免费观看| 午夜福利高清视频| 热99re8久久精品国产| 此物有八面人人有两片| tocl精华| 老熟妇乱子伦视频在线观看| 亚洲中文av在线| 午夜久久久在线观看| 99久久久亚洲精品蜜臀av| 国产精品av久久久久免费| 国产高清videossex| 亚洲少妇的诱惑av| 亚洲人成电影观看| 国产熟女xx| 国产一卡二卡三卡精品| 亚洲欧美激情综合另类| 亚洲 欧美 日韩 在线 免费| 国产亚洲精品久久久久5区| 啦啦啦观看免费观看视频高清 | 成年女人毛片免费观看观看9| 国产高清视频在线播放一区| 亚洲一码二码三码区别大吗| 日日摸夜夜添夜夜添小说| 久久婷婷人人爽人人干人人爱 | 村上凉子中文字幕在线| 久久午夜亚洲精品久久| 欧美精品亚洲一区二区| 叶爱在线成人免费视频播放| 亚洲 欧美一区二区三区| 91麻豆精品激情在线观看国产| 国产精品一区二区在线不卡| 99在线人妻在线中文字幕| 真人做人爱边吃奶动态| 香蕉丝袜av| 国产精品国产高清国产av| 国产精品久久视频播放| 国产成人欧美| x7x7x7水蜜桃| 日本欧美视频一区| 国产高清videossex| 亚洲成av片中文字幕在线观看| 亚洲 欧美一区二区三区| 午夜福利免费观看在线| 757午夜福利合集在线观看| a级毛片在线看网站| 国产亚洲精品久久久久久毛片| 成人手机av| 一进一出抽搐gif免费好疼| 国产麻豆成人av免费视频| 日本五十路高清| 国产精品亚洲av一区麻豆| 中文字幕人成人乱码亚洲影| 成人18禁高潮啪啪吃奶动态图| 国产精品二区激情视频| 在线国产一区二区在线| 国产成人精品久久二区二区91| 少妇 在线观看| 久久午夜亚洲精品久久| av中文乱码字幕在线| 婷婷丁香在线五月| 国产精品免费一区二区三区在线| 一级毛片精品| 男女之事视频高清在线观看| 久热爱精品视频在线9| tocl精华| 国产亚洲精品第一综合不卡| 亚洲国产精品久久男人天堂| 天堂√8在线中文| 黑人操中国人逼视频| 满18在线观看网站| 精品国产亚洲在线| aaaaa片日本免费| 男人舔女人下体高潮全视频| 亚洲情色 制服丝袜| 变态另类丝袜制服| 又黄又粗又硬又大视频| 这个男人来自地球电影免费观看| 中出人妻视频一区二区| 免费看a级黄色片| 国产亚洲欧美精品永久| 免费在线观看完整版高清| 日日爽夜夜爽网站| 免费女性裸体啪啪无遮挡网站| 欧美日韩精品网址| 欧美成狂野欧美在线观看| 嫩草影视91久久| 老司机午夜福利在线观看视频| 国产精品久久久久久精品电影 | 老熟妇仑乱视频hdxx| 成人国语在线视频| 丰满的人妻完整版| 日日爽夜夜爽网站| 18禁美女被吸乳视频| 又紧又爽又黄一区二区| 99香蕉大伊视频| 久久久国产欧美日韩av| av天堂久久9| 女生性感内裤真人,穿戴方法视频| 一本综合久久免费| 日韩大码丰满熟妇| 人人妻,人人澡人人爽秒播| 欧美另类亚洲清纯唯美| 国产精品二区激情视频| 一本综合久久免费| 日韩大码丰满熟妇| 亚洲伊人色综图| 国产激情久久老熟女| 亚洲久久久国产精品| 国产亚洲精品av在线| 国产精品 国内视频| 悠悠久久av| 成人国语在线视频| 欧美激情 高清一区二区三区| 看黄色毛片网站| 两个人视频免费观看高清| 久久精品91蜜桃| 日本vs欧美在线观看视频| 久久精品aⅴ一区二区三区四区| 怎么达到女性高潮| 久久人妻熟女aⅴ| 亚洲欧美激情在线| 亚洲成a人片在线一区二区| 两个人免费观看高清视频| 看免费av毛片| 国产高清有码在线观看视频 | 中文字幕久久专区| 亚洲欧美一区二区三区黑人| 在线观看www视频免费| 操美女的视频在线观看| 亚洲激情在线av| 午夜日韩欧美国产| 一级,二级,三级黄色视频| 天天躁夜夜躁狠狠躁躁| 午夜久久久在线观看| 亚洲精品粉嫩美女一区| 色婷婷久久久亚洲欧美| 午夜精品在线福利| 亚洲国产毛片av蜜桃av| 50天的宝宝边吃奶边哭怎么回事| 国产伦人伦偷精品视频| 国产一区二区三区视频了| 久久久国产成人精品二区| 18禁国产床啪视频网站| 女人爽到高潮嗷嗷叫在线视频| 这个男人来自地球电影免费观看| 午夜免费成人在线视频| 久久午夜综合久久蜜桃| 久久中文字幕人妻熟女| 亚洲电影在线观看av| 在线av久久热| 天堂动漫精品| 不卡av一区二区三区| 制服诱惑二区| 欧美在线黄色| 成人欧美大片| 啪啪无遮挡十八禁网站| 自线自在国产av| 成人18禁在线播放| 日韩视频一区二区在线观看| 国产亚洲av高清不卡| 日韩一卡2卡3卡4卡2021年| 午夜免费观看网址| 欧美乱码精品一区二区三区| 久久久久国内视频| 国产亚洲精品第一综合不卡| 久久久久久久久免费视频了| 18禁美女被吸乳视频| 亚洲国产欧美日韩在线播放| 成人亚洲精品av一区二区| 国产午夜精品久久久久久| 亚洲午夜精品一区,二区,三区| 久99久视频精品免费| 国产高清激情床上av| 丝袜美足系列| 午夜免费激情av| 波多野结衣av一区二区av| 国产成人av激情在线播放| 两个人视频免费观看高清| ponron亚洲| 国产精品永久免费网站| or卡值多少钱| 亚洲 欧美 日韩 在线 免费| 中出人妻视频一区二区| 国产一区二区在线av高清观看| 成年版毛片免费区| 国产麻豆69| 精品乱码久久久久久99久播| 欧美绝顶高潮抽搐喷水| 香蕉国产在线看| 日韩 欧美 亚洲 中文字幕| 国产成人免费无遮挡视频| 电影成人av| 亚洲精品一区av在线观看| xxx96com| 国产成人精品无人区| 大型黄色视频在线免费观看| 国产人伦9x9x在线观看| 日本免费a在线| 一进一出抽搐gif免费好疼| 欧美国产日韩亚洲一区| 久久中文字幕一级| 亚洲国产欧美日韩在线播放| √禁漫天堂资源中文www| 成在线人永久免费视频| 不卡一级毛片| 亚洲精品国产一区二区精华液| 90打野战视频偷拍视频| 亚洲三区欧美一区| 中出人妻视频一区二区| 视频区欧美日本亚洲| 亚洲三区欧美一区| 欧美成人免费av一区二区三区| 亚洲av日韩精品久久久久久密| 中文亚洲av片在线观看爽| 精品久久蜜臀av无| 色综合站精品国产| 国产精品免费视频内射| 18禁裸乳无遮挡免费网站照片 | 亚洲视频免费观看视频| 欧美绝顶高潮抽搐喷水| 国产精品美女特级片免费视频播放器 | 高清毛片免费观看视频网站| 97人妻精品一区二区三区麻豆 | 久久精品亚洲熟妇少妇任你| 18禁黄网站禁片午夜丰满| 黄片小视频在线播放| 少妇熟女aⅴ在线视频| 国产麻豆69| 亚洲熟妇中文字幕五十中出| 黄色成人免费大全| 狂野欧美激情性xxxx| 女性生殖器流出的白浆|