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

    基于HEC-HMS模型的溫榆河流域水文模擬

    2018-02-01 15:15:10廖如婷胡珊珊杜龍剛黃振芳
    南水北調(diào)與水利科技 2018年6期

    廖如婷 胡珊珊 杜龍剛 黃振芳

    摘要:水文模擬對于流域防洪減災(zāi)、水資源規(guī)劃管理有重要的意義。將分布式水文模型HEC-HMS應(yīng)用到北京溫榆河流域,模擬該區(qū)域1980-2005年間6場暴雨洪水過程和1983-1995年日降雨徑流過程,評估HEC-HMS模型在模擬場次洪水和日徑流過程兩種時(shí)間尺度上的適用性。結(jié)果表明,場次洪水模擬的洪峰流量和洪量相對誤差均在20%以內(nèi),峰現(xiàn)時(shí)差均不超過2 h,平均Nash效率系數(shù)為0.82,平均相關(guān)系數(shù)r為0.92;日徑流模擬率定期和驗(yàn)證期的平均Nash效率系數(shù)為0.6,平均相關(guān)系數(shù)r為0.78,平均相對誤差為3.95 %。分析結(jié)果表明,HEC-HMS模型在溫榆河流域的適用性較好,能夠有效模擬北方地區(qū)短時(shí)和長時(shí)降雨徑流過程,可用于該區(qū)域的洪水預(yù)報(bào)和水資源評價(jià)與管理。

    關(guān)鍵詞:HEC-HMS模型;溫榆河流域;次洪模擬;日徑流模擬;參數(shù)率定

    中圖分類號:TV121文獻(xiàn)標(biāo)志碼:A開放科學(xué)(資源服務(wù))標(biāo)識(shí)碼(OSID):[TP寥如婷.TIF]

    Hydrological simulation of Wenyu River basin based on HEC-HMS model

    LIAO Ruting1,2,HU Shanshan1,2,DU Longgang.3,HUANG Zhenfang.3

    (1.College of Resource Environment and Tourism,Capital Normal University,Beijing 100048,China;2.Beijing Laboratory of Water Resources Safety,Capital Normal University,Beijing 100048,China;3.Beijing Hydrological Station,Beijing 100089,China)

    Abstract:Hydrological simulation is important for flood hazard mitigation and water resources management.In this study,we adopted the distributed hydrological model HEC-HMS to simulate the flood events and daily runoff of the Wenyu River basin in Beijing,China.The model was calibrated and validated with the data of 6 storm events during 1980-2005 and daily stream flow during 1983-1995.The results showed that for all of the 6 floods,the relative errors of peak flow and runoff amount were less than 20%,and the peak time difference was within 2 h.The average Nash efficiency coefficient and average correlation coefficient r were 0.82 and 0.92,respectively.For the daily runoff simulation,the average Nash efficiency coefficient was 0.6,the average correlation coefficient r was 0.78,and the average relative error was 3.95%.Preliminary analysis showed that the HEC-HMS model applied well to the Wenyu River basin.It can effectively simulate the short-term and long-term rainfall and runoff processes in northern China,and can be used for flood forecast,water resources evaluation and management in this area.

    Key words:HEC-HMS model;Wenyu River basin;flood simulation;daily runoff simulation;calibration

    20世紀(jì)以來,隨著地[HJ2.15mm]理信息系統(tǒng)和遙感技術(shù)的迅速發(fā)展,分布式水文模型成為水文水資源研究的有效手段。HEC-HMS是由美國陸軍兵團(tuán)水文工程中心 (Hydrologic Engineering Center) 開發(fā)的半分布式水文模型,因其物理機(jī)制明確、結(jié)構(gòu)簡單、能概化模擬流域的復(fù)雜下墊面特性,已經(jīng)被廣泛的應(yīng)用于不同國家和地區(qū)。國外學(xué)者們應(yīng)用HEC-HMS模型進(jìn)行了一系列流域水文相關(guān)的研究,在沙特阿拉伯Makkah[1-2]、印度Chenab[3]、伊朗Jafar-Abad[4]等流域模擬效果良好。此外,HEC-HMS模型評估城市化對洪水的影響研究在印度Thirusoolam[5]、斯里蘭卡Kelani[6]、沙特阿拉伯Jazan[7]、孟買Oshiwara[8]等流域取得了很大進(jìn)展。在國內(nèi),HEC-HMS水文模型先后應(yīng)用于米塞、萬安流域[9]、晉江流域[10]、江河流域[11]、紫荊關(guān)流域[12]等,效果都較為理想。近些年,我國學(xué)者開始將HEC-HMS模型應(yīng)用于流域土地覆被變化的洪水響應(yīng)研究[13-14]、山洪預(yù)報(bào)[15-17]和長期徑流過程模擬[18-19]等方面,都取得了較好的研究結(jié)果。從國內(nèi)外研究現(xiàn)狀來看,HEC-HMS模型目前主要應(yīng)用于流域洪水預(yù)報(bào)及下墊面變化影響分析,在流域長時(shí)間尺度的水資源模擬與管理中應(yīng)用較少。本文將HEC-HMS模型應(yīng)用到北京溫榆河流域,分別進(jìn)行場次洪水模擬和日徑流模擬,分析評價(jià)HEC-HMS模型在兩個(gè)時(shí)間尺度上的模擬效果,探討并總結(jié)HEC-HMS模型在北方城市化影響區(qū)域的適用性和局限性,為HEC-HMS模型在類似區(qū)域的應(yīng)用研究提供參考。

    1研究區(qū)概況

    溫榆河流域位于北京市北部,發(fā)源于燕山南麓,屬于海河流域的北運(yùn)河水系,全長約47.5 km,流域面積2 478 km.2[20]。本文選取溫榆河沙河閘水文站控制流域?yàn)檠芯繀^(qū),匯水面積約1 099 km.2(圖1),其中山區(qū)(高程大于100 m)面積約708.3 km.2,占研究區(qū)面積的64.4 %;平原區(qū)面積約390.7 km.2,占35.6 %,跨北京市昌平、延慶和海淀三個(gè)區(qū)縣,20世紀(jì)80年代以來隨著北京衛(wèi)星城鎮(zhèn)的擴(kuò)展,流域內(nèi)的土地利用發(fā)生了明顯改變,建成區(qū)面積大大增加,目前已占到整個(gè)流域面積的19.5 %。流域?qū)儆跍貛Т箨懶约撅L(fēng)氣候,多年平均降雨量為642.5 mm,降雨多集中在6月-9月,約占全年的84 %,多年平均蒸發(fā)量為1 815.5 mm (20 cm口徑蒸發(fā)皿),其中4月-6月占全年的42 %[21]。

    2溫榆河流域HEC-HMS模型構(gòu)建

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

    研究區(qū)DEM(數(shù)字高程模型)數(shù)據(jù)和Landsat遙感影像數(shù)據(jù)來源于中國科學(xué)院地理空間數(shù)據(jù)云平臺(tái)[22],空間分辨率均為30 m。DEM數(shù)據(jù)用于子流域的劃分和流域特征的提取,遙感影像用于解譯得到研究區(qū)土地利用類型,本文將研究區(qū)土地利用類型分為草地、耕地、林地、水域、建設(shè)用地和未利用地6大類。土壤數(shù)據(jù)是由南京土壤所提供的第二次全國土地調(diào)查1∶400萬土壤數(shù)據(jù),主要土壤類型為褐土和黃壚土。通過查找土壤中黏粒、砂粒、有機(jī)質(zhì)等的含量,建立土壤類型屬性數(shù)據(jù)表。

    研究區(qū)水文氣象資料主要來源于北京市水文總站,部分摘自歷年《海河流域水文資料(第3卷)》。水文氣象資料包括溫榆河流域上游1980-2005年間8個(gè)雨量站及流域出口沙河閘水文站的實(shí)時(shí)、日降雨徑流數(shù)據(jù),其中流量數(shù)據(jù)根據(jù)沙河水庫水位庫容曲線反推入庫流量得到[23-24]。蒸發(fā)數(shù)據(jù)采用蘇莊站1980-2005年的20 cm蒸發(fā)皿日蒸發(fā)數(shù)據(jù)。

    2.2模型構(gòu)建

    本文在HEC-HMS模型的擴(kuò)展模塊HEC-GeoHMS的支持下,依托ArcGIS軟件將研究區(qū)劃分為6個(gè)子流域,并提取子流域及河道的特征參數(shù)。HEC-HMS模型主要由流域模塊(Basin Model)、氣象模塊(Meteorological Model)、控制運(yùn)行模塊(Control Specifications)以及時(shí)間序列數(shù)據(jù)模塊(Time-series Data)組成。其中,流域模塊是對實(shí)際流域上降雨徑流過程的模型概化。本文根據(jù)降雨徑流的形成過程,將其劃分為降雨損失、直接徑流、基流和河道匯流四個(gè)部分,分別采用SCS-CN曲線法、SCS單位線法、指數(shù)退水法和馬斯京根法進(jìn)行計(jì)算。氣象模塊采用權(quán)重因子法對降雨數(shù)據(jù)分析處理,建立雨量站與各個(gè)子流域的地理關(guān)系??刂七\(yùn)行模塊設(shè)定模擬計(jì)算的時(shí)間步長為1 h。時(shí)間序列數(shù)據(jù)模塊輸入研究區(qū)中相關(guān)站點(diǎn)的坐標(biāo)、雨量、流量以及蒸發(fā)等實(shí)測資料,為整個(gè)模型的模擬計(jì)算提供數(shù)據(jù)基礎(chǔ)。

    本文建立的模型中主要需要確定的參數(shù)為CN、流域滯時(shí)以及馬斯京根法中的蓄量常數(shù)K和流量比重X。其中,CN值根據(jù)研究區(qū)土地利用和土壤數(shù)據(jù)查表[25]計(jì)算得出。流域滯時(shí)、蓄量常數(shù)K和流量比重X則需要在已有流域、河道特征參數(shù)的基礎(chǔ)上計(jì)算,然后進(jìn)行參數(shù)率定確定。

    2.3模型率定與驗(yàn)證

    為了使建立的HEC-HMS模型模擬結(jié)果達(dá)到最優(yōu),必須對模型參數(shù)進(jìn)行優(yōu)化。HEC-HMS模型提供了單變量梯度法(Univariate Gradient)與Nelder Mead法2種方法及七種目標(biāo)函數(shù)供用戶選擇,目標(biāo)函數(shù)包括峰值加權(quán)均方根誤差、洪峰流量百分比誤差和洪量百分比誤差等[26]。本文采用單變量梯度參數(shù)優(yōu)選法,并以峰值加權(quán)均方根誤差為目標(biāo)函數(shù)來進(jìn)行模型參數(shù)的優(yōu)化。

    在驗(yàn)證過程中,選用洪峰流量相對誤差(Dv[WTBX])、洪量相對誤差(Dp)、峰現(xiàn)時(shí)差(ΔT)、Nash效率系數(shù)(NSE)、相關(guān)系數(shù)(r)等5個(gè)指標(biāo)進(jìn)行模擬精度評價(jià)。跟據(jù)SL250-2000《水文情報(bào)預(yù)報(bào)規(guī)范》[27],洪峰流量相對誤差Dv和洪量相對誤差Dp在實(shí)測值的 20%以內(nèi)為合格;峰現(xiàn)時(shí)差ΔT的范圍為 3 h,其中Dv、Dp、ΔT的絕對值越小,表明模型的模擬精準(zhǔn)度越高;Nash效率系數(shù)NSE和相關(guān)系數(shù)r的范圍是0到1,值越大,表示模擬效果越好[28]。

    3結(jié)果與分析

    根據(jù)上述方法構(gòu)建溫榆河流域HEC-HMS模型,其中場次洪水模擬采用1980-2005年間的4場暴雨洪水?dāng)?shù)據(jù)對模型的參數(shù)進(jìn)行率定,2場暴雨洪水?dāng)?shù)據(jù)進(jìn)行驗(yàn)證;日降雨徑流模擬采用1983-1988年的日徑流數(shù)據(jù)對模型的參數(shù)進(jìn)行率定,1990-1995年的日徑流數(shù)據(jù)對模型進(jìn)行驗(yàn)證。

    3.1場次洪水模擬結(jié)果分析

    從模擬結(jié)果可以看出(表1),在率定期和驗(yàn)證期這6場洪水中,洪峰流量相對誤差Dv和洪量相對誤差Dp基本上都控制在20%以內(nèi),峰現(xiàn)時(shí)差ΔT均不超過2 h;相關(guān)系數(shù)r均大于0.85;率定期和驗(yàn)證期的平均Nash系數(shù)NSE為0.82。根據(jù)《水文情報(bào)預(yù)報(bào)規(guī)范》(SL 250-2000)[27],本文構(gòu)建的HEC-HMS模型次洪模擬結(jié)果比較理想。

    結(jié)合實(shí)測與模擬流量過程對比圖(圖2)可以看出,模型對于峰高量大的單峰洪水模擬效果較好,洪號19880802、19940711和19980705這三場洪水的相關(guān)系數(shù)和Nash效率系數(shù)都在0.9左右,峰現(xiàn)時(shí)間上多表現(xiàn)為雨停峰現(xiàn)。洪號19860626和19940717模擬洪水過程與實(shí)測洪水過程誤差較大,分析其原因可能是受實(shí)測流量精度和時(shí)間步長影響,由于實(shí)測流量數(shù)據(jù)為經(jīng)過反推的入庫流量,與真實(shí)值之間可能存在一定誤差;且實(shí)測流量的時(shí)間步長不一致,模型中為線性插值結(jié)果,坦化了實(shí)際的流量過程線,導(dǎo)致局部模擬過程與實(shí)測過程相差較大。

    3.2日徑流過程模擬結(jié)果分析

    從表2中可以看出,率定期日徑流相關(guān)系數(shù)為0.74、Nash效率系數(shù)為0.55,相對誤差為-5.66 %,率定期日流量模擬效果較好,總體模擬值略低于實(shí)測值;驗(yàn)證期日徑流相關(guān)系數(shù)為0.81、Nash效率系數(shù)為0.65,相對誤差為2.23 %,表明驗(yàn)證期總體模擬值略高于實(shí)測值,且模擬效果較優(yōu)于率定期。

    圖3和圖4分別為率定期和驗(yàn)證期的模擬與實(shí)測逐日流量過程線的比較,從圖中可以看出,日徑流模擬值與實(shí)測日流量過程線擬合程度較好,模擬效果整體令人滿意。部分日流量峰值模擬偏大,可能是由于流域內(nèi)水保措施的修建、人類活動(dòng)頻繁等多方面因素造成了徑流量減少。

    4結(jié)論與討論

    (1)本文利用研究區(qū)氣象、水文、土壤、土地利用數(shù)據(jù),構(gòu)建了溫榆河流域HEC-HMS模型。通過對模型參數(shù)的不斷優(yōu)化與率定,場次洪水模擬的平均Nash效率系數(shù)為0.82,平均相關(guān)系數(shù)r為0.92;日徑流模擬的平均Nash效率系數(shù)為0.6,平均相關(guān)系數(shù)r為0.78。表明HEC-HMS模型能夠較為準(zhǔn)確的反映研究區(qū)的水文過程,具有較好的適用性。

    (2)HEC-HMS模型對于場次洪水模擬和日降雨徑流模擬都具有較好的適用性,可以用于流域洪水預(yù)報(bào)和水資源評價(jià)與管理。

    (3)流域內(nèi)高強(qiáng)度的人類活動(dòng)和下墊面改變導(dǎo)致的產(chǎn)匯流變化是影響模型模擬精度的重要原因,HEC-HMS模型在城市化地區(qū)的改進(jìn)是進(jìn)一步研究的方向。

    參考文獻(xiàn)(References):

    [1]DAWOD G M,MIRZA M N,ALGHAMDI K A.GIS-Based spatial mapping of flash flood hazard in Makkah City,Saudi Arabia[J].Journal of Geographic Information System,2011,3(3):225-231.

    [2]AL-GHAMDI K A,ELZAHRANY R A,MIRZA M N,et al.Impacts of urban growth on flood hazards in Makkah City,Saudi Arabia[J].International Journal of Water Resources & Environmental Engineering,2012,4:23-34.

    [3]SHAHID M A,BOCCARDO P,USMAN M,et al.Predicting peak flows in real time through event based hydrologic modeling for a trans-boundary river catchment[J].Water Resources Management,2017,31(3):793-810.DOI:10.1007/s11269-016-1435-2.

    [4]MOSTAFAZADEH R,SADODDIN A,BAHREMAND A,et al.Scenario analysis of flood control structures using a multi-criteria decision-making technique in Northeast Iran[J].Natural Hazards,2017,87(3):1827-1846.DOI:10.1007/s11069-017-2851-1.

    [5]SURIYA S,MUDGAL B V.Impact of urbanization on flooding:The thirusoolam sub watershed-A case study[J].Journal of Hydrology,2012,412(1):210-219.

    [6]SILVA M M G T D,WEERAKOON S B,HERATH S.Modeling of event and continuous flow hydrographs with HEC-HMS:Case study in the Kelani River Basin,Sri Lanka[J].Journal of Hydrologic Engineering,2014,19(4):800-806.DOI:10.1061 / (ASCE) HE.1943-5584.0000846.

    [7]ALFY M E.Assessing the impact of arid area urbanization on flash floods using GIS,remote sensing,and HEC-HMS rainfall-runoff modeling[J].Hydrology Research,2016,47(6).DOI:10.2166/nh.2016.133.

    [8]ZOPE P E,ELDHO T I,JOTHIPRAKASH V.Impacts of land use-land cover change and urbanization on flooding:A case study of Oshiwara River Basin in Mumbai,India[J].Catena,2016,145:142-154.DOI:10.1016/j.catena.2016.06.009.

    [9]OLEYIBLO J O.Application of HEC-HMS for flood forecasting in Misai and Wan′an catchments in China[J].Water Science and Engineering,2010,3(1):14-22.

    [10][ZK(#]陳芬,林峰,陳興偉.采用分布式HEC-HMS水文模型的晉江流域暴雨次洪模擬[J].華僑大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,33(3):325-329.(CHEN F,LIN F,CHEN X W.Simulation of storm in Jinjiang Basin using HEC-HMS hydrological model[J].Journal of Huaqiao University (Natural Science),2012,33(3):325-329.(in Chinese))

    [11]JIN H,LIANG R,WANG Y,et al.Flood-Runoff in Semi-Arid and Sub-Humid Regions,a case study&58;A simulation of Jianghe Watershed in Northern China[J].Water,2015,7(9):5155-5172.DOI:10.3390/w7095155.

    [12]趙永超,王加虎,梁菊平,等.HEC-HMS模型在紫荊關(guān)流域水文模擬中的應(yīng)用[J].水電能源科學(xué),2017(12):10-13.(ZHAO Y C,WANG J H,LIANG J P,et al.Application of HEC-HMS model in hydrological simulation of the Zijinguan Basin[J].Hydroelectric Energy Science,2017(12):10-13.(in Chinese)) DOI:1000-7709 (2017) 12-0010-04.

    [13]林木生.東南沿海晉江西溪流域洪水徑流演變與模擬分析[D].福州:福建師范大學(xué),2015.(LIN M S.Flood Runoff Evolution and Simulation Analysis of Xixi River Basin in Southeast Coast of Jinjiang[D].Fuzhou:Fujian Normal University,2015.(in Chinese))

    [14]REN J,ZHENG X,CHEN P,et al.An Investigation into Sub-Basin Rainfall Losses in different underlying surface conditions using HEC-HMS:A case study of a loess hilly region in Gedong Basin in the Western Shanxi Province of China[J].Water,2017,2(1):1-27.DOI:10.3390/w9110870.

    [15]ZHANG H L,WANG Y J,WANG Y Q,et al.Quantitative comparison of semi-?and fully-distributed hydrologic models in simulating flood hydrographs on a mountain watershed in southwest China[J].Journal of Hydrodynamics,2013,25(6):877-885.

    [16]趙睿.山丘區(qū)小流域暴雨洪水分析計(jì)算方法應(yīng)用研究[D].濟(jì)南:山東大學(xué),2015.(ZHAO R.Application of analysis and calculation method for storm flood in small watershed at hill area[D].Jinan:Shandong University,2015.(in Chinese))

    [17]胡國華,陳肖,于澤興,等.基于HEC-HMS的郴江流域山洪預(yù)報(bào)研究[J].自然災(zāi)害學(xué)報(bào),2017(3):147-155.(HU G H,CHEN X,YU Z X,et al.Research on mountain flood forecasting in Lancang River Basin based on HEC-HMS[J].Journal of Natural Disasters,2017(3):147-155.(in Chinese)) DOI:10.13577/j.jnd.2017.0318.

    [18]吳建紅,朱積軍.HEC-HMS模型及其應(yīng)用比較研究[J].科技創(chuàng)新導(dǎo)報(bào),2010(4):6-7.(WU J H,ZHU J J.Application and comparison of HEC-HMS model[J].Science and Technology Innovation Herald,2010(4):6-7.(in Chinese))

    [19]DU J K,LI Q,HAN Y R,et al.Assessing the effects of urbanization on annual runoff and flood events using an integrated hydrological modeling system for Qinhuai River basin,China[J].Journal of Hydrology,2012,464-465(5):127-139.

    [20]郁達(dá)偉,于淼,魏源送,等.1980~2010年溫榆河的水環(huán)境質(zhì)量時(shí)空演變特征[J].環(huán)境科學(xué)學(xué)報(bào),2012,32(11):2803-2813.(YU D W,YU M,WEI Y S,et al.Temporal and spatial evolution characteristics of water environmental quality in Wenyu River from 1980 to 2010[J].Journal of Environmental Science,2012,32(11):2803-2813.(in Chinese))

    [21]唐莉華,張思聰,呂賢弼.溫榆河上游流域的洪水模擬及不確定性分析[C]// 第二屆全國水問題研究學(xué)術(shù)研討會(huì)論文集,2009:168-174.(TANG L H,ZHANG S C,LV X B.Simulation and uncertainty analysis of flood in upper reaches of Wenyu River[C]// National Water Research Symposium.2009:168-174.(in Chinese))

    [22]地理空間數(shù)據(jù)云[EB/OL]http://www.gscloud.cn.(Geospatial data Cloud[EB/OL]http://www.gscloud.cn.(in Chinese))

    [23]張帆,張君友.淺談反推入庫洪水方法[J].治淮,2010(2):24-25.(ZHANG F,ZHANG J Y.Discussion on the method of ?backstepping?reservoir-inflow[J].Governance Huai River,2010(2):24-25.(in Chinese))

    [24]陳昌杰,姜和平.應(yīng)用Excel軟件擬合庫容曲線[J].小水電,2002(6):19-20.(CHEN C J,JIANG H P.Application of Excel software to fit reservoir capacity curves[J].Small Hydro Power,2002(6):19-20.(in Chinese))

    [25]李向新.HEC-HMS水文建模系統(tǒng)原理·方法·應(yīng)用[M].中國水利水電出版社,2015.(LI X X.Principles,methods and applications of HEC-HMS hydrological modeling system[M].China Water Power Press,2015.(in Chinese))

    [26]鄧霞,董曉華,薄會(huì)娟.目標(biāo)函數(shù)對HEC-HMS模型參數(shù)率定的影響研究[J].水電能源科學(xué),2010,28(8):23-25.(DENG X,DONG X H,BO H J.Influence of objective function on parameter calibration of HEC-HMS model[J].Hydroelectric Energy,2010,28(8):23-25.(in Chinese))

    [27]SL 250-2000,水文情報(bào)預(yù)報(bào)規(guī)范[S].(SL 250-2000,Specification for hydrological information forecast[S].(in Chinese))

    [28]梁睿.HEC-HMS水文模型在北張店流域的應(yīng)用研究[D].太原:太原理工大學(xué),2012.(LIANG R.Application of HEC-HMS hydrological model in the Beizhangdian watershed[D].Taiyuan:Taiyuan University of Technology,2012.(in Chinese))

    亚洲人成网站在线观看播放| 国产黄片视频在线免费观看| 亚洲精品自拍成人| 乱人伦中国视频| 综合色丁香网| 亚洲av国产av综合av卡| 亚洲中文av在线| 欧美精品高潮呻吟av久久| 美女cb高潮喷水在线观看| 亚洲av福利一区| 亚洲精品自拍成人| 一个人看视频在线观看www免费| 亚洲成人一二三区av| 欧美日韩国产mv在线观看视频| xxx大片免费视频| 69精品国产乱码久久久| 夜夜看夜夜爽夜夜摸| 国产成人午夜福利电影在线观看| 中文字幕制服av| 亚洲精品国产av成人精品| 久久人人爽人人片av| 亚洲精品乱久久久久久| 桃花免费在线播放| 欧美日韩成人在线一区二区| 日本与韩国留学比较| 亚洲成人av在线免费| 丝袜脚勾引网站| 亚洲精品国产av成人精品| 婷婷色综合www| 卡戴珊不雅视频在线播放| 精品少妇久久久久久888优播| 欧美精品人与动牲交sv欧美| 精品一区二区三卡| 一个人免费看片子| 免费人成在线观看视频色| 成年女人在线观看亚洲视频| 午夜视频国产福利| 欧美老熟妇乱子伦牲交| av线在线观看网站| 嘟嘟电影网在线观看| 亚洲人成网站在线播| 精品亚洲乱码少妇综合久久| 国产视频内射| 国产日韩欧美视频二区| 日本欧美国产在线视频| 春色校园在线视频观看| 成年人午夜在线观看视频| 亚洲内射少妇av| 26uuu在线亚洲综合色| 成年美女黄网站色视频大全免费 | 精品视频人人做人人爽| 在线观看www视频免费| 18禁在线无遮挡免费观看视频| 欧美精品国产亚洲| av福利片在线| 午夜日本视频在线| 美女脱内裤让男人舔精品视频| 久久久欧美国产精品| 日日啪夜夜爽| 建设人人有责人人尽责人人享有的| 美女内射精品一级片tv| 五月玫瑰六月丁香| 婷婷色综合大香蕉| 日韩av在线免费看完整版不卡| 精品一品国产午夜福利视频| 丝袜喷水一区| 乱码一卡2卡4卡精品| 极品少妇高潮喷水抽搐| tube8黄色片| 黄色欧美视频在线观看| 80岁老熟妇乱子伦牲交| 中国美白少妇内射xxxbb| 亚洲人成77777在线视频| 色婷婷久久久亚洲欧美| 好男人视频免费观看在线| 夜夜骑夜夜射夜夜干| 亚洲欧美中文字幕日韩二区| 日日爽夜夜爽网站| 国产一区二区在线观看av| 在线观看免费日韩欧美大片 | 最近的中文字幕免费完整| 色婷婷久久久亚洲欧美| 久久影院123| 亚洲成人手机| 亚洲精品av麻豆狂野| 精品少妇久久久久久888优播| 女人精品久久久久毛片| av又黄又爽大尺度在线免费看| 十分钟在线观看高清视频www| 美女大奶头黄色视频| 免费人成在线观看视频色| 少妇人妻精品综合一区二区| 我要看黄色一级片免费的| 亚洲精品av麻豆狂野| 精品一区在线观看国产| 十分钟在线观看高清视频www| av福利片在线| 大香蕉久久成人网| 精品人妻偷拍中文字幕| 美女脱内裤让男人舔精品视频| 日韩,欧美,国产一区二区三区| 亚洲精品久久成人aⅴ小说 | 国产日韩欧美在线精品| 免费黄网站久久成人精品| 99热全是精品| 亚洲精品国产av成人精品| 亚洲婷婷狠狠爱综合网| 亚洲欧洲日产国产| 亚洲综合色惰| 国产精品免费大片| 久久97久久精品| 99久久精品国产国产毛片| 久久久精品区二区三区| 中文字幕制服av| 大又大粗又爽又黄少妇毛片口| 国产探花极品一区二区| 色网站视频免费| 成人毛片a级毛片在线播放| 狠狠婷婷综合久久久久久88av| 在线天堂最新版资源| 欧美成人午夜免费资源| 亚洲经典国产精华液单| 亚洲欧美成人综合另类久久久| 男女边吃奶边做爰视频| 美女内射精品一级片tv| 美女福利国产在线| 三级国产精品欧美在线观看| 成人亚洲精品一区在线观看| 久久久久久久久久人人人人人人| 欧美日韩一区二区视频在线观看视频在线| 欧美人与性动交α欧美精品济南到 | 丝瓜视频免费看黄片| 日韩一区二区视频免费看| 欧美日本中文国产一区发布| av视频免费观看在线观看| 国产成人精品在线电影| av网站免费在线观看视频| 又大又黄又爽视频免费| 成人影院久久| 亚洲国产最新在线播放| 啦啦啦啦在线视频资源| 天天影视国产精品| 久久久久久久久久成人| 成人毛片60女人毛片免费| 国产精品偷伦视频观看了| 26uuu在线亚洲综合色| 久久久精品94久久精品| videossex国产| 简卡轻食公司| 国产精品一区www在线观看| 久久久久久久大尺度免费视频| 精品视频人人做人人爽| 久久精品人人爽人人爽视色| 久久99精品国语久久久| 七月丁香在线播放| 激情五月婷婷亚洲| 天堂8中文在线网| 国产精品三级大全| 成人影院久久| 久久精品久久精品一区二区三区| 视频中文字幕在线观看| av国产精品久久久久影院| 亚洲综合精品二区| 美女cb高潮喷水在线观看| 水蜜桃什么品种好| 婷婷色av中文字幕| 精品人妻一区二区三区麻豆| 欧美日本中文国产一区发布| 边亲边吃奶的免费视频| 极品人妻少妇av视频| 午夜免费鲁丝| 丝袜脚勾引网站| 免费黄频网站在线观看国产| 伦精品一区二区三区| 中文天堂在线官网| av福利片在线| 久久久久久久久久人人人人人人| 国产又色又爽无遮挡免| 岛国毛片在线播放| 日产精品乱码卡一卡2卡三| 精品午夜福利在线看| 欧美日韩国产mv在线观看视频| 久久精品久久久久久久性| 熟妇人妻不卡中文字幕| av专区在线播放| 99国产综合亚洲精品| 国产一区二区三区av在线| 最新中文字幕久久久久| √禁漫天堂资源中文www| 精品国产国语对白av| 中文字幕免费在线视频6| 日日摸夜夜添夜夜添av毛片| 国产免费福利视频在线观看| 国产伦精品一区二区三区视频9| 日本猛色少妇xxxxx猛交久久| 在线观看国产h片| 久久影院123| 亚洲,一卡二卡三卡| 热re99久久国产66热| 一本久久精品| 国产片特级美女逼逼视频| 熟女人妻精品中文字幕| 秋霞在线观看毛片| 国产精品蜜桃在线观看| 精品99又大又爽又粗少妇毛片| av线在线观看网站| 国产乱人偷精品视频| 97在线人人人人妻| 国产日韩欧美在线精品| 赤兔流量卡办理| 久久久久国产网址| 成人亚洲欧美一区二区av| 女性生殖器流出的白浆| 欧美日韩一区二区视频在线观看视频在线| 久久免费观看电影| 久久久久久久久大av| 亚洲综合色网址| 美女xxoo啪啪120秒动态图| 亚洲人与动物交配视频| 亚洲三级黄色毛片| 欧美三级亚洲精品| 色吧在线观看| 91精品国产国语对白视频| 三级国产精品片| 免费看不卡的av| 国产精品国产三级专区第一集| 婷婷色麻豆天堂久久| 欧美少妇被猛烈插入视频| 九色亚洲精品在线播放| 蜜桃在线观看..| 这个男人来自地球电影免费观看 | 中文字幕制服av| 久久婷婷青草| 日本-黄色视频高清免费观看| 极品人妻少妇av视频| 在线观看免费高清a一片| 欧美成人精品欧美一级黄| 亚洲性久久影院| av在线观看视频网站免费| 69精品国产乱码久久久| 黑人欧美特级aaaaaa片| 国产 精品1| 美女中出高潮动态图| 十八禁网站网址无遮挡| 精品人妻偷拍中文字幕| 色吧在线观看| 中文乱码字字幕精品一区二区三区| 蜜桃久久精品国产亚洲av| 欧美老熟妇乱子伦牲交| 亚洲精品久久久久久婷婷小说| 日韩三级伦理在线观看| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产成人一精品久久久| 久久国产精品男人的天堂亚洲 | 青春草视频在线免费观看| 不卡视频在线观看欧美| 69精品国产乱码久久久| 黑人欧美特级aaaaaa片| 国产黄色视频一区二区在线观看| 国产淫语在线视频| 精品熟女少妇av免费看| 国产精品国产三级国产av玫瑰| 亚洲成人手机| 伊人久久国产一区二区| 国产精品久久久久久久电影| 99热网站在线观看| 亚洲成人手机| 欧美3d第一页| 香蕉精品网在线| 欧美激情 高清一区二区三区| 亚洲精品,欧美精品| 久久人人爽av亚洲精品天堂| 妹子高潮喷水视频| 精品少妇黑人巨大在线播放| 国产 一区精品| 91久久精品国产一区二区成人| 精品一区二区三区视频在线| 天天影视国产精品| 我的老师免费观看完整版| 国产男人的电影天堂91| 久久国产亚洲av麻豆专区| 老司机亚洲免费影院| 久久人人爽人人片av| 日本欧美视频一区| 你懂的网址亚洲精品在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧美一区二区三区国产| 丰满乱子伦码专区| 大话2 男鬼变身卡| 99热国产这里只有精品6| 久久久久精品国产欧美久久久| 精品国产国语对白av| 天堂俺去俺来也www色官网| 在线观看人妻少妇| 1024香蕉在线观看| 国产免费福利视频在线观看| 国产欧美日韩综合在线一区二区| av一本久久久久| 又黄又粗又硬又大视频| 国产高清视频在线播放一区| 91精品三级在线观看| 51午夜福利影视在线观看| 91麻豆精品激情在线观看国产 | 热99国产精品久久久久久7| avwww免费| 国产av精品麻豆| 18禁黄网站禁片午夜丰满| 精品国产乱子伦一区二区三区| 亚洲中文字幕日韩| 亚洲视频免费观看视频| 中文字幕人妻熟女乱码| 国产成人精品久久二区二区免费| 亚洲精品在线观看二区| 两个人免费观看高清视频| 十八禁网站网址无遮挡| 大香蕉久久成人网| √禁漫天堂资源中文www| 国产精品99久久99久久久不卡| 久久av网站| 精品国产乱码久久久久久小说| 咕卡用的链子| avwww免费| 超色免费av| 久久热在线av| 午夜激情av网站| 亚洲,欧美精品.| av福利片在线| 一本综合久久免费| 亚洲欧美日韩另类电影网站| 少妇被粗大的猛进出69影院| 狠狠狠狠99中文字幕| 亚洲情色 制服丝袜| 老汉色∧v一级毛片| 精品国产一区二区三区久久久樱花| 亚洲成a人片在线一区二区| 熟女少妇亚洲综合色aaa.| 热99国产精品久久久久久7| 人人妻,人人澡人人爽秒播| 岛国毛片在线播放| 日本wwww免费看| 两个人看的免费小视频| 免费看十八禁软件| 亚洲国产欧美一区二区综合| 成人国产av品久久久| 欧美日韩黄片免| 另类亚洲欧美激情| 国产伦人伦偷精品视频| 亚洲国产av新网站| 午夜福利视频精品| 黄色 视频免费看| 精品乱码久久久久久99久播| 中国美女看黄片| 美女扒开内裤让男人捅视频| 中文字幕精品免费在线观看视频| 韩国精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 欧美日韩视频精品一区| 天堂8中文在线网| 飞空精品影院首页| 极品教师在线免费播放| 大片电影免费在线观看免费| 黄网站色视频无遮挡免费观看| 国产成人精品在线电影| 欧美中文综合在线视频| 9191精品国产免费久久| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲九九香蕉| 人妻久久中文字幕网| 激情视频va一区二区三区| 啦啦啦 在线观看视频| 国产精品免费大片| 国产亚洲精品一区二区www | 少妇 在线观看| av线在线观看网站| 纵有疾风起免费观看全集完整版| 色视频在线一区二区三区| av天堂在线播放| 日韩大码丰满熟妇| 91精品国产国语对白视频| 精品亚洲成国产av| 国产伦人伦偷精品视频| 淫妇啪啪啪对白视频| 丁香六月天网| 这个男人来自地球电影免费观看| 少妇猛男粗大的猛烈进出视频| 老熟女久久久| 一夜夜www| www.精华液| 午夜精品国产一区二区电影| 亚洲黑人精品在线| 黄色丝袜av网址大全| 水蜜桃什么品种好| 久久久国产精品麻豆| 熟女少妇亚洲综合色aaa.| 国产亚洲一区二区精品| 91麻豆精品激情在线观看国产 | 亚洲午夜精品一区,二区,三区| avwww免费| 亚洲av国产av综合av卡| 久久人妻av系列| 悠悠久久av| 国产精品一区二区在线观看99| 亚洲伊人久久精品综合| 一本大道久久a久久精品| 国产又爽黄色视频| a级毛片在线看网站| 国产高清videossex| 欧美 日韩 精品 国产| 啦啦啦 在线观看视频| 视频区图区小说| 侵犯人妻中文字幕一二三四区| 天天添夜夜摸| 一本久久精品| 日韩欧美三级三区| 少妇的丰满在线观看| 女人爽到高潮嗷嗷叫在线视频| 精品一区二区三卡| 制服诱惑二区| 亚洲专区字幕在线| 美女高潮喷水抽搐中文字幕| 国产单亲对白刺激| 纯流量卡能插随身wifi吗| 久久精品国产亚洲av高清一级| 91老司机精品| 亚洲伊人色综图| 久久亚洲真实| 久久久久久久久免费视频了| 丝瓜视频免费看黄片| 久久精品亚洲av国产电影网| 日韩免费高清中文字幕av| 国产91精品成人一区二区三区 | 每晚都被弄得嗷嗷叫到高潮| 日韩欧美免费精品| 正在播放国产对白刺激| 女人高潮潮喷娇喘18禁视频| 国产黄频视频在线观看| www日本在线高清视频| 亚洲熟妇熟女久久| 一区二区av电影网| 亚洲欧美精品综合一区二区三区| 黄色视频不卡| 国产高清videossex| 精品亚洲成a人片在线观看| 国产黄频视频在线观看| 曰老女人黄片| 90打野战视频偷拍视频| 亚洲情色 制服丝袜| 国产精品香港三级国产av潘金莲| 亚洲国产欧美一区二区综合| 怎么达到女性高潮| 精品亚洲成国产av| 1024视频免费在线观看| 亚洲精品av麻豆狂野| 丁香六月天网| 亚洲av国产av综合av卡| 久久 成人 亚洲| 日韩视频在线欧美| 久久九九热精品免费| 高潮久久久久久久久久久不卡| 极品少妇高潮喷水抽搐| 人成视频在线观看免费观看| 人人妻人人添人人爽欧美一区卜| 黑人巨大精品欧美一区二区蜜桃| 国产成人精品久久二区二区91| 天堂8中文在线网| 久久久国产欧美日韩av| 亚洲黑人精品在线| 啦啦啦视频在线资源免费观看| 在线观看一区二区三区激情| 在线观看免费日韩欧美大片| 欧美日韩国产mv在线观看视频| 在线观看免费视频日本深夜| 免费在线观看日本一区| kizo精华| 日韩中文字幕视频在线看片| 在线亚洲精品国产二区图片欧美| 少妇裸体淫交视频免费看高清 | 免费在线观看黄色视频的| 九色亚洲精品在线播放| 亚洲男人天堂网一区| 中文字幕最新亚洲高清| 在线观看免费视频网站a站| 精品久久蜜臀av无| 国产深夜福利视频在线观看| 国产av又大| 一二三四在线观看免费中文在| 欧美黑人欧美精品刺激| 多毛熟女@视频| 国产成人欧美| 久久99一区二区三区| 人人妻人人添人人爽欧美一区卜| 精品人妻1区二区| 在线观看免费日韩欧美大片| 久久精品亚洲精品国产色婷小说| 宅男免费午夜| 日韩一卡2卡3卡4卡2021年| 欧美成人午夜精品| 91av网站免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久网色| 国产三级黄色录像| 国产片内射在线| 无遮挡黄片免费观看| 亚洲天堂av无毛| 成年人午夜在线观看视频| 成人免费观看视频高清| 久久久精品94久久精品| 国产亚洲精品第一综合不卡| 后天国语完整版免费观看| 一边摸一边抽搐一进一出视频| 日韩熟女老妇一区二区性免费视频| 伦理电影免费视频| 国产一区有黄有色的免费视频| 欧美日韩一级在线毛片| 国产精品99久久99久久久不卡| 黄色视频,在线免费观看| 成人18禁高潮啪啪吃奶动态图| 国产在线一区二区三区精| 精品国产乱码久久久久久小说| 成人特级黄色片久久久久久久 | 菩萨蛮人人尽说江南好唐韦庄| 亚洲专区中文字幕在线| 欧美日韩一级在线毛片| 美女视频免费永久观看网站| 国产色视频综合| 日本黄色视频三级网站网址 | 一边摸一边抽搐一进一小说 | 变态另类成人亚洲欧美熟女 | 日本撒尿小便嘘嘘汇集6| 欧美日韩黄片免| 久久精品成人免费网站| 在线十欧美十亚洲十日本专区| 成人国产av品久久久| 精品久久久久久久毛片微露脸| 一级毛片精品| 亚洲熟妇熟女久久| 黑人欧美特级aaaaaa片| 国产不卡av网站在线观看| 国产精品欧美亚洲77777| 美女午夜性视频免费| 狂野欧美激情性xxxx| 国产日韩欧美亚洲二区| 国产成人一区二区三区免费视频网站| 久久久久久免费高清国产稀缺| 亚洲精品自拍成人| 国产精品一区二区精品视频观看| 老鸭窝网址在线观看| 黄色成人免费大全| 又大又爽又粗| 国产亚洲一区二区精品| 欧美精品一区二区免费开放| 999精品在线视频| 69精品国产乱码久久久| 狠狠精品人妻久久久久久综合| 操出白浆在线播放| 999久久久精品免费观看国产| 在线看a的网站| 99九九在线精品视频| 国产一区二区激情短视频| 国产一区二区 视频在线| 人妻 亚洲 视频| 午夜激情久久久久久久| 国产视频一区二区在线看| 国产精品.久久久| 91麻豆av在线| 另类精品久久| 欧美激情高清一区二区三区| 国产片内射在线| 成人国产一区最新在线观看| 国产成人精品久久二区二区91| 黑人欧美特级aaaaaa片| 乱人伦中国视频| 久久天躁狠狠躁夜夜2o2o| 高潮久久久久久久久久久不卡| 天天躁夜夜躁狠狠躁躁| av网站免费在线观看视频| 精品少妇内射三级| 午夜福利,免费看| 久久久久国产一级毛片高清牌| 在线亚洲精品国产二区图片欧美| 一级片免费观看大全| 精品国产国语对白av| 亚洲,欧美精品.| 大片电影免费在线观看免费| 伊人久久大香线蕉亚洲五| 亚洲一区二区三区欧美精品| 午夜福利影视在线免费观看| 高清视频免费观看一区二区| www.精华液| 男女午夜视频在线观看| 别揉我奶头~嗯~啊~动态视频| 在线十欧美十亚洲十日本专区| 欧美黄色片欧美黄色片| 一区二区av电影网| 一二三四在线观看免费中文在| 欧美日韩亚洲综合一区二区三区_| 久久影院123| av在线播放免费不卡| 18禁黄网站禁片午夜丰满| 丰满饥渴人妻一区二区三| 高清黄色对白视频在线免费看| 亚洲专区字幕在线| 久久久久久久精品吃奶| 又紧又爽又黄一区二区| 超色免费av| 999精品在线视频| 一区在线观看完整版| 丝瓜视频免费看黄片| 亚洲第一欧美日韩一区二区三区 | 久久久欧美国产精品| 国产高清国产精品国产三级| 精品少妇黑人巨大在线播放| 亚洲精品av麻豆狂野| 免费高清在线观看日韩| 91九色精品人成在线观看|