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

    VIC模型在西江流域徑流模擬中的應(yīng)用

    2016-03-23 05:30:33山紅翠孫飛飛
    中國(guó)農(nóng)村水利水電 2016年4期
    關(guān)鍵詞:西江流域匯流水文站

    山紅翠,袁 飛, 盛 東,孫飛飛

    (1.湖南省水利水電科學(xué)研究所,長(zhǎng)沙 410007;2.河海大學(xué)水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,南京 210098;3.寧波弘泰水利信息科技有限公司,浙江 寧波 315100)

    水文模型是水文科學(xué)與當(dāng)今計(jì)算機(jī)科學(xué)有機(jī)結(jié)合產(chǎn)生的應(yīng)用產(chǎn)物,是信息革命創(chuàng)造的水文學(xué)的一個(gè)全新的研究領(lǐng)域[1]。根據(jù)對(duì)復(fù)雜水文過(guò)程的概化方式、概化程度及概化假設(shè)的不同,可將流域水文模型分為集總式水文模型和分布式水文模型。分布式水文模型因其可以同時(shí)考慮降水的空間分布不均勻性和流域下墊面條件的不均勻性,對(duì)徑流預(yù)報(bào)具有很大的優(yōu)越性[2,3]。為了盡可能提高對(duì)流域產(chǎn)匯流預(yù)報(bào)精度,分布式水文模型近年來(lái)得到了廣泛的應(yīng)用[4]。

    本文利用DEM、植被覆蓋/土地利用、土壤等下墊面資料,基于0.25°×0.25°網(wǎng)格分辨率,選取西江流域?yàn)槔捎肰IC模型與匯流模型耦合,進(jìn)行參數(shù)率定與流量過(guò)程的模擬,探討模型在西江流域的徑流模擬效果與適用性。鑒于氣候變化對(duì)水資源影響研究中常選取1961-1990年作為氣候基準(zhǔn)期[5-7],故本文的研究時(shí)期為1961-1990年,以期為氣候變化對(duì)流域水文過(guò)程的影響研究提供科學(xué)參考。

    1 研究區(qū)概況

    西江全長(zhǎng)2 214 km,流域總面積為34.6萬(wàn)km2,占珠江全流域的75.6%[8]。西江流域?qū)贌釒Ъ皝啛釒夂颍?-8月降水總量占全年降水量的62.5%,各月降水量均超過(guò)200 mm,其中6月是一年中降水量最大的月份,達(dá)到261.5 mm[9]。西江流域徑流以降水作為補(bǔ)給來(lái)源,降水的空間分布決定了徑流的空間分布[8]。西江流域暴雨多由地面冷鋒、靜止鋒、高空切變線、低渦和熱帶氣旋等天氣系統(tǒng)形成,具有強(qiáng)度大、次數(shù)多、歷時(shí)長(zhǎng)等特點(diǎn)。暴雨、大暴雨或特大暴雨多出現(xiàn)在每年的4-10月,約占全年暴雨總次數(shù)的58%。西江流域洪水多由暴雨形成,且出現(xiàn)時(shí)間與暴雨一致。洪水過(guò)程以多峰型為主,有峰高、量大、歷時(shí)長(zhǎng)的特點(diǎn),一次較大的洪水過(guò)程約歷時(shí)30~40 d。

    本文選取西江流域武宣水文站以上集水區(qū)域?yàn)檠芯繀^(qū),研究區(qū)集水面積196 255 km2。研究區(qū)內(nèi)共有氣象站點(diǎn)19個(gè),主要水文控制站12個(gè),如圖1所示。

    圖1 西江流域武宣站以上區(qū)域站網(wǎng)分布圖Fig.1 The stations network distribution in Xijiang River basin

    2 模型簡(jiǎn)介

    可變下滲能力水文模型VIC[10](Variable Infiltration Capacity)是由美國(guó)華盛頓大學(xué)、加利福尼亞大學(xué)伯克利分校、普里斯頓大學(xué)共同發(fā)展的基于物理機(jī)制的、基于網(wǎng)格的大尺度分布式水文模型。該模型考慮了大氣-植被-土壤之間的相互作用過(guò)程,反映土壤、植被、大氣中水熱狀態(tài)變化和水熱傳輸[4],并通過(guò)網(wǎng)格化,分別考慮每個(gè)計(jì)算網(wǎng)格內(nèi)多種植被覆蓋類型,以及土壤特性和降雨的空間變異性對(duì)徑流的影響。

    VIC模型由一層土壤的VIC-1L模型和兩層土壤的VIC-2L模型,發(fā)展為具有三層土壤的VIC-3L模型。VIC-3L模型在VIC-2L的土壤頂層分出一個(gè)頂薄層[11](常取0.1 m),它允許土壤層與層之間土壤水的擴(kuò)散。VIC模型是一個(gè)基于空間分布網(wǎng)格化的水文模型,其格網(wǎng)結(jié)構(gòu)便于與氣候模型嵌套以評(píng)價(jià)氣候變化對(duì)水資源的影響[12]。VIC模型結(jié)構(gòu)如圖2所示,VIC模型與水文相關(guān)的變量有土壤蒸發(fā)E,植物散發(fā)Et、地表截流蒸發(fā)Ec、潛熱通量L、感熱通量S、長(zhǎng)波輻射RL、短波輻射RS、地表熱通量G、下滲i、滲透Q、地表徑流R、基流B。

    圖2 VIC模型結(jié)構(gòu)圖[13]Fig.2 VIC model structure chart

    實(shí)際應(yīng)用中,一般先采用VIC模型進(jìn)行水量平衡計(jì)算,輸出研究區(qū)各網(wǎng)格的徑流和蒸發(fā),再與匯流模型耦合,將網(wǎng)格上的產(chǎn)流轉(zhuǎn)化為流域出口斷面的流量過(guò)程。

    3 數(shù)據(jù)與方法

    3.1 數(shù)據(jù)準(zhǔn)備

    本文將研究區(qū)劃分為351個(gè)0.25°×0.25°的網(wǎng)格,VIC模型在每個(gè)網(wǎng)格內(nèi)計(jì)算產(chǎn)流。驅(qū)動(dòng)VIC模型需要準(zhǔn)備的輸入文件包括:植被參數(shù)庫(kù)文件、植被參數(shù)文件、土壤參數(shù)文件,氣象強(qiáng)迫數(shù)據(jù)文件和全局控制文件[5]。

    VIC模型將陸地表面用N+1種陸面覆蓋類型來(lái)描述,n=1,2,…,N表示N種植被覆蓋類型,n=N+1代表裸土。VIC模型中植被參數(shù)通過(guò)植被參數(shù)庫(kù)文件和植被參數(shù)文件描述。本文中植被參數(shù)庫(kù)中各參數(shù)主要根據(jù)陸面數(shù)據(jù)同化系統(tǒng)(Land Data Assimilation Systems,LADS)確定。植被參數(shù)文件描述各個(gè)網(wǎng)格內(nèi)各種植被類型的面積比例、葉面積指數(shù)等信息。在模型中,網(wǎng)格內(nèi)總的蒸散發(fā)通過(guò)對(duì)各種地表覆蓋類型上的蒸散發(fā)進(jìn)行面積加權(quán)平均來(lái)計(jì)算。

    VIC模型的土壤參數(shù)文件用來(lái)描述土壤的空間差異性。土壤分類根據(jù)NOAA(the National Oceanic and Atmospheric Administration)辦公室提供的全球5′土壤質(zhì)地分類描述。土壤參數(shù)文件中包含每個(gè)網(wǎng)格的主要土壤類型參數(shù)、平均高程、平均年降水量等信息。本文中土壤參數(shù)主要根據(jù)Cosby[14]的成果確定。

    本文中氣象強(qiáng)迫數(shù)據(jù)通過(guò)將19個(gè)氣象站點(diǎn)1961-1990年的實(shí)測(cè)日最高、日最低氣溫和日降水?dāng)?shù)據(jù)按反距離平方插值方法插值得到。

    全局控制文件描述模型運(yùn)行的時(shí)間步長(zhǎng)、模擬起始時(shí)間、各參數(shù)文件的路徑等信息,用以引導(dǎo)模型運(yùn)行。

    3.2 模擬方法

    上述文件準(zhǔn)備完畢后則可驅(qū)動(dòng)VIC模型產(chǎn)流部分的運(yùn)行。VIC模型雖具有一定的物理機(jī)制,但由于自然界能量及水文過(guò)程的復(fù)雜性,模型往往采用簡(jiǎn)化的物理方程或經(jīng)驗(yàn)方程來(lái)描述水分及能量的變化過(guò)程,因此有部分參數(shù)需要水文資料進(jìn)行率定[15]。需要率定的產(chǎn)流參數(shù)主要有模型頂薄層、上層和下層土壤深度、可變下滲曲線指數(shù)binf、最大基流速度Dsmax、產(chǎn)生非線性基流與Dsmax的比值Ds等。底層土壤采用Arno[16]模型的土壤水模塊計(jì)算基流。

    本文的模型匯流部分,采用自由水蓄水庫(kù)將總徑流劃分為地表徑流、壤中流和地下徑流,采用馬斯京根分段連續(xù)演算法進(jìn)行河道匯流演算。匯流過(guò)程需率定參數(shù)包括自由水蓄水容量SM、自由水蓄水容量曲線指數(shù)EX、地下水和壤中流的日出流系數(shù)KG和KI、地面徑流消退系數(shù)CS、壤中流消退系數(shù)CI、地下水消退系數(shù)CG、馬斯京根法參數(shù)Kmus和Xmus。

    水文站點(diǎn)之間按照匯流次序演算,上游站點(diǎn)的實(shí)測(cè)徑流值采用河道演算方法演算至下游站點(diǎn)。按照所收集到12水文站點(diǎn)各自的徑流資料序列,將其劃分為參數(shù)率定期和驗(yàn)證期。本文采用確定性系數(shù)(Nash-Sutcliffe效率系數(shù),NSE)最大為目標(biāo)函數(shù)率定VIC模型產(chǎn)、匯流參數(shù),采用確定性系數(shù)NSE與相對(duì)誤差Bias判別其模擬效果。確定性系數(shù)與相對(duì)誤差計(jì)算方法分別如下式。

    (2)

    4 結(jié)果分析

    4.1 多年徑流模擬

    按照上述方法進(jìn)行產(chǎn)匯流計(jì)算,得到12個(gè)水文站率定期與驗(yàn)證期日徑流模擬結(jié)果,如表1所示。

    表1 VIC模型模擬西江流域12個(gè)水文站參數(shù)率定期和驗(yàn)證期日徑流過(guò)程Tab.1 The result of day runoff processes during parameters calibrationperiod and verification period simulated by VIC model

    由表1可知VIC模型模擬研究區(qū)12個(gè)水文站的多年徑流過(guò)程結(jié)果如下:①率定期7個(gè)站點(diǎn)徑流模擬的確定性系數(shù)(NSE)在0.9以上,3個(gè)站點(diǎn)在0.7~0.9,2個(gè)站點(diǎn)在0.6~0.7;率定期12個(gè)水文站點(diǎn)模擬值與實(shí)測(cè)值的相對(duì)誤差(Bias)均在±5%之間。②驗(yàn)證期8個(gè)站點(diǎn)的確定性系數(shù)(NSE)高于0.75;除小龍?zhí)墩竞桶嗣镜腂ias偏高, 其余10個(gè)站點(diǎn)的Bias均在±10%。根據(jù)表1,本研究采用的VIC模型能夠較準(zhǔn)確地模擬西江流域多年日流量過(guò)程。

    小龍?zhí)?、八茂、這洞等站模擬值與實(shí)測(cè)徑流值相比,NSE偏低(或相對(duì)誤差偏大),分析其原因,可能是這些站點(diǎn)區(qū)間內(nèi)雨量站數(shù)目較少,降水輸入誤差會(huì)較大程度地影響徑流模擬的精度。

    4.2 年內(nèi)徑流模擬

    以研究區(qū)中、下游區(qū)域、資料系列較長(zhǎng)且模擬效果相對(duì)較好的天峨、武宣兩個(gè)水文站點(diǎn)為例,分析其年內(nèi)(以模擬相對(duì)穩(wěn)定期1975年為例)日徑流模擬效果,結(jié)果如圖3和圖4所示。

    圖3 VIC模型模擬天峨站1975年日流量過(guò)程Fig.3 Observed and simulated hydrograph by VIC model at Tian′e Station in 1975

    圖4 VIC模型模擬武宣站1975年日流量過(guò)程Fig.4 Observed and simulated hydrograph by VIC model at Wuxuan Station in 1975

    由圖3和圖4顯示,VIC模型能較好地模擬西江流域1975年的日徑流過(guò)程:天峨站徑流過(guò)程模擬NSE為0.92,相對(duì)誤差為6.78%;武宣站徑流過(guò)程模擬NSE為0.95,相對(duì)誤差為2.71%。但由圖3和圖4可以看出VIC模型對(duì)兩個(gè)站點(diǎn)的洪峰模擬系統(tǒng)偏小。主要原因是這些峰值高的洪水過(guò)程多是由暴雨形成,而研究區(qū)內(nèi)數(shù)量有限的氣象站點(diǎn)可能不能夠準(zhǔn)確記錄暴雨中心的雨量,因此輸入的降水值可能偏小,很可能引起洪峰模擬系統(tǒng)偏小。對(duì)枯季徑流,VIC模型模擬值與實(shí)測(cè)值之間的相對(duì)偏差較為明顯。主要原因可能為,本文采用確定性系數(shù)最大為目標(biāo)函數(shù),通常能有效地率定高流量過(guò)程,所率定的參數(shù)可能不能有效表征枯季徑流的水文特征。

    總體而言,VIC模型能較好地再現(xiàn)西江流域的水文過(guò)程,模擬效果較好。

    5 結(jié) 語(yǔ)

    西江流域水資源呈現(xiàn)年際間豐枯交替、年內(nèi)水資源空間分布不均等特點(diǎn),并導(dǎo)致流域水旱災(zāi)害的威脅與水資源短缺、水污染嚴(yán)重困擾并存[17]。此外,由于氣候變化或氣候異常的影響,水資源的供需矛盾日益突出。因此,氣候變化對(duì)流域水資源的影響研究已成為目前水文水資源研究的熱點(diǎn)。而評(píng)價(jià)氣候變化對(duì)水文水資源的影響時(shí),水文模型是最為常用且有效的工具,故研究不同區(qū)域水文模型的適用性可為氣候變化影響研究奠定科學(xué)基礎(chǔ)。

    本文采用VIC模型模擬西江流域基準(zhǔn)期(1961-1990年)多年日徑流過(guò)程,確定性系數(shù)較高,相對(duì)誤差較 ??;模型對(duì)兩個(gè)水文站點(diǎn)的洪峰估計(jì)均偏低,但整體誤差較小,模擬精度較高。因此,可采用VIC模型模擬西江流域的徑流過(guò)程,并做水文預(yù)報(bào)分析及氣候變化的水文響應(yīng)研究。

    [1] 芮孝芳, 蔣成煜, 張金存.流域水文模型的發(fā)展[J].水文,2006,26(3):22-26.

    [2] 李恒鵬,王旭強(qiáng),楊桂山,等.基于單元格網(wǎng)的STREAM分布式水文模型及其應(yīng)用----以太湖上游西苕溪流域?yàn)槔齕J].長(zhǎng)江流域資源與環(huán)境,2007,16(6):715-720.

    [3] 吳 波,張萬(wàn)昌.基于柵格產(chǎn)匯流的TOPMODEL[J].長(zhǎng)江流域資源與環(huán)境,2007,16(5) :661-666.

    [4] 張利平,陳小鳳,張曉琳,等.VIC模型與SWAT模型在中小流域徑流模擬中的對(duì)比研究[J].長(zhǎng)江流域資源與環(huán)境,2009,18(8):745-752.

    [5] 劉 謙.VIC大尺度陸面水文模型在中國(guó)區(qū)域的應(yīng)用[D].長(zhǎng)沙:湖南大學(xué),2004.

    [6] 趙芳芳,徐宗學(xué). 黃河源區(qū)未來(lái)氣候變化的水文響應(yīng)[J]. 資源環(huán)境,2009,31(5):722-730.

    [7] 張利平,秦琳琳,胡志芳,等.南水北調(diào)中線工程水源區(qū)水文循環(huán)過(guò)程對(duì)氣候變化的響應(yīng)[J].水利學(xué)報(bào),2010,41(11):1 261-1 271.

    [8] 張 峰.西江流域分布式水文模型及其應(yīng)用研究[D].上海:東華大學(xué),2012.

    [9] 何 慧,陸 虹,歐 藝.1959-2008年廣西西江流域洪澇氣候特征[J].氣候變化研究進(jìn)展,2009,5(3):134-138.

    [10] Liang X,Wood E F,Lettenmaier D P,et al.The project for intercomparison of land-surface parameterization schemes Phase2(c) Red-arkansas River basin experiment: spatial and temporal analysis of energy fluxes[J].Global and Planet Change,1998,19:137-159.

    [11] Liang X,Lettenmaier D P,Wood E F.Surface soil moisture parameterization of the VIC-2L model:evalutation and modification[J].Global Planet Change, 1996,13:195-206.

    [12] 胡彩虹,郭生練,彭定志,等.VIC模型在流域徑流模擬中的應(yīng)用[J].人民黃河, 2005,27(10): 22-28.

    [13] Variable Inflitration Capacity(VIC) Macroscale Hydrological Model[EB/OL].http:∥www.hydro. washington.edu/Lettenmaier/Models/VIC/.

    [14] Cosby B J,Hornberger G M,Clapp R B,et al.A statistical exploration of the relationships of soil moisture characteristics to the physical properties of soils[J].Water Resources Research,1984,20:682-690.

    [15] 張磊磊,郝振純,童 凱,等.VIC模型在三江源地區(qū)產(chǎn)匯流模擬中的應(yīng)用[J].水電能源科學(xué)[J].2013,31(1):18-20.

    [16] Todini E.The ARNO rain-runoff model[J].Journal of Hydrology,1996,175:339-382.

    [17] 張繼軍.VIC模型在中國(guó)濕潤(rùn)地區(qū)的應(yīng)用研究[D].南京:河海大學(xué),2006.

    猜你喜歡
    西江流域匯流水文站
    廣西西江流域旅游經(jīng)濟(jì)與生態(tài)環(huán)境 耦合協(xié)調(diào)發(fā)展時(shí)空分異研究
    SL流量計(jì)在特殊河段的應(yīng)用——以河源水文站為例
    體驗(yàn)經(jīng)濟(jì)視角下西江流域非物質(zhì)文化遺產(chǎn)旅游開(kāi)發(fā)研究
    廣西西江流域生態(tài)補(bǔ)償問(wèn)題與對(duì)策研究
    流域匯流計(jì)算方法在風(fēng)峪沙河段的應(yīng)用探討
    郭家屯水文站單斷沙關(guān)系分析
    草壩水文站兩種蒸發(fā)器對(duì)比分析
    韓府灣水文站報(bào)汛方案
    基于ARM CortexM0的智能光伏匯流采集裝置設(shè)計(jì)及應(yīng)用
    一種球載雷達(dá)匯流環(huán)設(shè)計(jì)
    男男h啪啪无遮挡| 人妻 亚洲 视频| 97在线人人人人妻| 如何舔出高潮| 国产免费视频播放在线视频| 青春草亚洲视频在线观看| 欧美国产精品一级二级三级 | 青春草亚洲视频在线观看| 国产黄a三级三级三级人| 亚洲美女视频黄频| 午夜福利高清视频| 国产精品久久久久久久电影| 天堂中文最新版在线下载 | 欧美激情国产日韩精品一区| 久久久a久久爽久久v久久| 国产在线一区二区三区精| 国产精品国产三级国产专区5o| 国产精品熟女久久久久浪| 热99国产精品久久久久久7| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品国产成人久久av| 99热网站在线观看| 草草在线视频免费看| 国产真实伦视频高清在线观看| 亚洲欧美精品自产自拍| 女人久久www免费人成看片| 看十八女毛片水多多多| 在线精品无人区一区二区三 | 成人美女网站在线观看视频| 特大巨黑吊av在线直播| 人妻系列 视频| 国产精品99久久99久久久不卡 | 在线观看av片永久免费下载| 免费av不卡在线播放| 国产探花在线观看一区二区| 晚上一个人看的免费电影| 久久久色成人| 美女xxoo啪啪120秒动态图| 国产亚洲91精品色在线| 草草在线视频免费看| 日韩中字成人| 久久久午夜欧美精品| 你懂的网址亚洲精品在线观看| 熟女av电影| 亚洲图色成人| 免费黄网站久久成人精品| 成年免费大片在线观看| 国产精品一区www在线观看| 国产av国产精品国产| 国产成人精品婷婷| 国产精品秋霞免费鲁丝片| 成人黄色视频免费在线看| 亚洲精品乱码久久久久久按摩| 大香蕉久久网| 一级毛片aaaaaa免费看小| 久久久a久久爽久久v久久| 伊人久久国产一区二区| 成年av动漫网址| 亚洲精品中文字幕在线视频 | 国产人妻一区二区三区在| 女人十人毛片免费观看3o分钟| 听说在线观看完整版免费高清| 搡女人真爽免费视频火全软件| 久久久亚洲精品成人影院| 日韩欧美 国产精品| 亚洲国产精品成人综合色| 久久精品综合一区二区三区| 国产成人福利小说| 亚洲精品乱久久久久久| 亚洲精品色激情综合| 草草在线视频免费看| 草草在线视频免费看| 精品国产乱码久久久久久小说| 久久久午夜欧美精品| 欧美亚洲 丝袜 人妻 在线| 自拍偷自拍亚洲精品老妇| 天堂中文最新版在线下载 | 精华霜和精华液先用哪个| 日本午夜av视频| 色吧在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 2018国产大陆天天弄谢| 日韩av免费高清视频| 男女国产视频网站| 夫妻午夜视频| 一级毛片久久久久久久久女| 久久精品久久精品一区二区三区| 一级二级三级毛片免费看| 国产免费视频播放在线视频| av国产久精品久网站免费入址| 午夜福利视频1000在线观看| 精品久久久久久久末码| 精品久久久久久久末码| 99久久精品国产国产毛片| 国产淫片久久久久久久久| 国产精品一区www在线观看| 少妇 在线观看| 超碰97精品在线观看| 狠狠精品人妻久久久久久综合| 亚洲人成网站高清观看| 免费在线观看成人毛片| 麻豆精品久久久久久蜜桃| 免费看av在线观看网站| 国产日韩欧美亚洲二区| 国产大屁股一区二区在线视频| tube8黄色片| 国产日韩欧美在线精品| 国产日韩欧美在线精品| 国产伦精品一区二区三区四那| 久热这里只有精品99| 人人妻人人看人人澡| 亚洲,欧美,日韩| 少妇高潮的动态图| 人妻夜夜爽99麻豆av| 久久久午夜欧美精品| 亚洲精品成人久久久久久| 丰满人妻一区二区三区视频av| 色哟哟·www| 亚洲色图综合在线观看| 最后的刺客免费高清国语| 色婷婷久久久亚洲欧美| 韩国av在线不卡| 如何舔出高潮| 男女那种视频在线观看| 国产日韩欧美在线精品| 另类亚洲欧美激情| 久久久久国产网址| 人妻夜夜爽99麻豆av| 日韩一本色道免费dvd| 亚洲精品日韩av片在线观看| 最新中文字幕久久久久| 啦啦啦在线观看免费高清www| 九色成人免费人妻av| 91久久精品国产一区二区三区| 欧美激情在线99| 亚洲,一卡二卡三卡| 好男人视频免费观看在线| 色婷婷久久久亚洲欧美| 久久99热6这里只有精品| 国产欧美日韩一区二区三区在线 | 国产精品一区www在线观看| 成年免费大片在线观看| 在线观看一区二区三区激情| 久久精品国产亚洲av涩爱| 国产老妇女一区| 国产成人精品久久久久久| 只有这里有精品99| 久久精品久久久久久噜噜老黄| 久久久久久久久大av| 97在线视频观看| 国产精品女同一区二区软件| 欧美xxxx黑人xx丫x性爽| 免费黄色在线免费观看| 欧美日韩亚洲高清精品| 夜夜爽夜夜爽视频| 黄色配什么色好看| 寂寞人妻少妇视频99o| 久久99精品国语久久久| av免费在线看不卡| 又粗又硬又长又爽又黄的视频| 久久人人爽av亚洲精品天堂 | 99九九线精品视频在线观看视频| 日本爱情动作片www.在线观看| 亚洲aⅴ乱码一区二区在线播放| 热99国产精品久久久久久7| 久久精品国产a三级三级三级| 青春草国产在线视频| 人妻系列 视频| 国产日韩欧美亚洲二区| 一级爰片在线观看| 精品久久久久久久人妻蜜臀av| 久久久久久久亚洲中文字幕| 午夜免费男女啪啪视频观看| .国产精品久久| 亚洲精品亚洲一区二区| 亚洲欧美成人综合另类久久久| 国产精品一及| 噜噜噜噜噜久久久久久91| 婷婷色综合www| 美女视频免费永久观看网站| 亚洲在线观看片| 白带黄色成豆腐渣| 神马国产精品三级电影在线观看| 亚洲av成人精品一二三区| 日韩av免费高清视频| 色播亚洲综合网| 亚洲欧美一区二区三区国产| 香蕉精品网在线| 亚洲三级黄色毛片| 久久久久久久久久久免费av| 边亲边吃奶的免费视频| 男女下面进入的视频免费午夜| 激情 狠狠 欧美| 午夜福利视频1000在线观看| 蜜桃久久精品国产亚洲av| 亚洲av中文字字幕乱码综合| 在线a可以看的网站| 国产极品天堂在线| 日日撸夜夜添| 超碰av人人做人人爽久久| 久久久久国产精品人妻一区二区| av黄色大香蕉| 日韩,欧美,国产一区二区三区| 日韩av不卡免费在线播放| 在线观看三级黄色| 色播亚洲综合网| 免费在线观看成人毛片| 精品人妻一区二区三区麻豆| 成年女人在线观看亚洲视频 | 亚洲欧美日韩卡通动漫| 51国产日韩欧美| 真实男女啪啪啪动态图| 一个人看视频在线观看www免费| 男的添女的下面高潮视频| 色吧在线观看| 免费观看无遮挡的男女| 免费黄色在线免费观看| 欧美97在线视频| 亚洲经典国产精华液单| 国产黄频视频在线观看| 午夜激情久久久久久久| 97热精品久久久久久| 91精品国产九色| 在线播放无遮挡| 久久久亚洲精品成人影院| 亚洲精品乱久久久久久| 男男h啪啪无遮挡| 一级av片app| 亚洲av中文字字幕乱码综合| 亚洲不卡免费看| 亚洲av一区综合| 看十八女毛片水多多多| 色婷婷久久久亚洲欧美| 久久久午夜欧美精品| 黄色欧美视频在线观看| 精品视频人人做人人爽| 欧美激情久久久久久爽电影| 亚洲精品aⅴ在线观看| 国产精品久久久久久精品电影| 午夜免费鲁丝| 国产亚洲一区二区精品| 成人午夜精彩视频在线观看| 午夜福利网站1000一区二区三区| 五月玫瑰六月丁香| 一二三四中文在线观看免费高清| 久久精品夜色国产| a级一级毛片免费在线观看| 1000部很黄的大片| 七月丁香在线播放| 嫩草影院入口| 又爽又黄无遮挡网站| 小蜜桃在线观看免费完整版高清| 成人毛片a级毛片在线播放| 亚洲av免费在线观看| 亚洲婷婷狠狠爱综合网| 精品午夜福利在线看| 免费看光身美女| 中文字幕av成人在线电影| 久久久久久久午夜电影| 国产老妇伦熟女老妇高清| 精品久久久噜噜| 插逼视频在线观看| 国产精品久久久久久精品电影| 欧美日韩国产mv在线观看视频 | 国产欧美日韩一区二区三区在线 | 日本熟妇午夜| 精华霜和精华液先用哪个| 69人妻影院| 香蕉精品网在线| 国产成人精品福利久久| 男女下面进入的视频免费午夜| av在线播放精品| 国产v大片淫在线免费观看| 伊人久久精品亚洲午夜| 黄片wwwwww| 欧美丝袜亚洲另类| 狠狠精品人妻久久久久久综合| 成人免费观看视频高清| 亚洲精品日韩av片在线观看| 丝袜美腿在线中文| 97在线人人人人妻| 午夜福利在线观看免费完整高清在| 国产欧美另类精品又又久久亚洲欧美| 亚洲最大成人手机在线| 国产在线男女| 国产精品无大码| 91精品伊人久久大香线蕉| 熟女av电影| 九色成人免费人妻av| videos熟女内射| 日韩中字成人| 国产精品一二三区在线看| 国产精品一区二区性色av| 亚洲精品影视一区二区三区av| 国产毛片在线视频| 欧美另类一区| 黄色欧美视频在线观看| 久久久久久久久大av| 欧美日韩亚洲高清精品| 自拍偷自拍亚洲精品老妇| 水蜜桃什么品种好| 制服丝袜香蕉在线| 精品久久久精品久久久| 网址你懂的国产日韩在线| 午夜免费鲁丝| 日韩成人伦理影院| 欧美国产精品一级二级三级 | 美女视频免费永久观看网站| 日日啪夜夜撸| 久久久色成人| 高清午夜精品一区二区三区| 亚洲精品影视一区二区三区av| 内射极品少妇av片p| 激情 狠狠 欧美| 成人欧美大片| 中文精品一卡2卡3卡4更新| 涩涩av久久男人的天堂| 国产一区二区在线观看日韩| 干丝袜人妻中文字幕| 国产乱人视频| 男人爽女人下面视频在线观看| 国产精品久久久久久久电影| 久久午夜福利片| 精品酒店卫生间| 校园人妻丝袜中文字幕| a级一级毛片免费在线观看| 一区二区av电影网| 中文字幕亚洲精品专区| 精品人妻视频免费看| 韩国av在线不卡| 中文字幕av成人在线电影| 中文乱码字字幕精品一区二区三区| 全区人妻精品视频| 成人无遮挡网站| 白带黄色成豆腐渣| 69av精品久久久久久| 一个人看视频在线观看www免费| 久久影院123| 另类亚洲欧美激情| 国产一区二区亚洲精品在线观看| 男人添女人高潮全过程视频| 五月玫瑰六月丁香| 日本一本二区三区精品| 一本色道久久久久久精品综合| av黄色大香蕉| 美女国产视频在线观看| 伊人久久国产一区二区| 久久99热这里只有精品18| 欧美高清性xxxxhd video| 亚洲精品日韩av片在线观看| 亚洲精品中文字幕在线视频 | 国产成人freesex在线| 青青草视频在线视频观看| 日韩成人伦理影院| 欧美性猛交╳xxx乱大交人| 国产精品人妻久久久久久| a级毛片免费高清观看在线播放| 青春草视频在线免费观看| 各种免费的搞黄视频| 国产大屁股一区二区在线视频| www.色视频.com| 久久精品国产亚洲av天美| 日韩精品有码人妻一区| 中文字幕亚洲精品专区| 欧美bdsm另类| 免费黄网站久久成人精品| 美女主播在线视频| 18禁动态无遮挡网站| 黄片wwwwww| 在线亚洲精品国产二区图片欧美 | 大又大粗又爽又黄少妇毛片口| 最近最新中文字幕免费大全7| 97在线人人人人妻| 高清日韩中文字幕在线| 观看免费一级毛片| 亚洲精品色激情综合| 国产成人a区在线观看| 亚洲一级一片aⅴ在线观看| 高清在线视频一区二区三区| av免费在线看不卡| 国产色爽女视频免费观看| 日韩三级伦理在线观看| 夜夜爽夜夜爽视频| 国产男女内射视频| 小蜜桃在线观看免费完整版高清| 亚洲最大成人av| 久久亚洲国产成人精品v| 久久ye,这里只有精品| 国产精品不卡视频一区二区| 久久综合国产亚洲精品| 国产成人freesex在线| 国产高清有码在线观看视频| 国产精品一区二区三区四区免费观看| 91精品伊人久久大香线蕉| 波野结衣二区三区在线| 精品久久久精品久久久| 亚洲欧美精品自产自拍| 国产精品久久久久久精品电影小说 | 尾随美女入室| 午夜免费鲁丝| 国产成人精品福利久久| av卡一久久| 人体艺术视频欧美日本| 免费观看的影片在线观看| 免费观看av网站的网址| 亚洲国产精品成人久久小说| 精品人妻偷拍中文字幕| 三级国产精品片| 国产真实伦视频高清在线观看| 超碰97精品在线观看| 少妇被粗大猛烈的视频| 亚洲久久久久久中文字幕| 日韩不卡一区二区三区视频在线| 中文字幕人妻熟人妻熟丝袜美| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产午夜精品久久久久久一区二区三区| 国产精品久久久久久精品古装| 看非洲黑人一级黄片| 国产久久久一区二区三区| 国产成人免费观看mmmm| 国产一区有黄有色的免费视频| 国产成人精品婷婷| 人人妻人人看人人澡| 超碰av人人做人人爽久久| 国产欧美亚洲国产| 97热精品久久久久久| 国产午夜福利久久久久久| 99久国产av精品国产电影| 22中文网久久字幕| 久久久久性生活片| 3wmmmm亚洲av在线观看| 97在线人人人人妻| 又粗又硬又长又爽又黄的视频| 亚洲av中文字字幕乱码综合| 国产综合精华液| 五月开心婷婷网| 国语对白做爰xxxⅹ性视频网站| 亚洲av日韩在线播放| 最近的中文字幕免费完整| 91aial.com中文字幕在线观看| 午夜精品国产一区二区电影 | 美女视频免费永久观看网站| 久久精品国产a三级三级三级| 男人舔奶头视频| 国产精品久久久久久精品古装| 少妇人妻 视频| 又黄又爽又刺激的免费视频.| 国产成人精品婷婷| 国产精品人妻久久久久久| 91在线精品国自产拍蜜月| 成人亚洲精品一区在线观看 | 99热这里只有是精品在线观看| av国产免费在线观看| 极品少妇高潮喷水抽搐| 一边亲一边摸免费视频| 久久久久久久久久成人| 午夜老司机福利剧场| 大话2 男鬼变身卡| 日本黄色片子视频| 少妇人妻久久综合中文| 中文字幕久久专区| 欧美日韩视频高清一区二区三区二| 国产白丝娇喘喷水9色精品| 久久久久九九精品影院| 亚洲精品乱码久久久久久按摩| 嘟嘟电影网在线观看| 精品一区二区免费观看| 中文字幕久久专区| 观看美女的网站| 成人亚洲欧美一区二区av| 在线免费十八禁| 黄片无遮挡物在线观看| 直男gayav资源| 日韩三级伦理在线观看| 亚洲欧美成人精品一区二区| 性插视频无遮挡在线免费观看| 大话2 男鬼变身卡| 免费看av在线观看网站| 我的女老师完整版在线观看| 另类亚洲欧美激情| 国产成人a∨麻豆精品| 天堂中文最新版在线下载 | 午夜福利网站1000一区二区三区| 成人鲁丝片一二三区免费| 日本wwww免费看| 一级毛片aaaaaa免费看小| 久久人人爽av亚洲精品天堂 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲综合色惰| 久久久精品免费免费高清| 啦啦啦中文免费视频观看日本| 高清午夜精品一区二区三区| 日韩一区二区三区影片| 少妇人妻久久综合中文| 另类亚洲欧美激情| 嫩草影院入口| 久久精品人妻少妇| 亚洲精品国产av蜜桃| 亚州av有码| 国产精品不卡视频一区二区| 欧美xxxx黑人xx丫x性爽| 午夜福利高清视频| 精品99又大又爽又粗少妇毛片| 中文精品一卡2卡3卡4更新| 麻豆成人av视频| 一区二区三区免费毛片| 亚洲av国产av综合av卡| 99热这里只有精品一区| 麻豆国产97在线/欧美| 看免费成人av毛片| 日韩亚洲欧美综合| av免费观看日本| 免费少妇av软件| 综合色丁香网| 亚洲欧美精品自产自拍| 亚洲精品456在线播放app| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩大片免费观看网站| 免费看不卡的av| 91精品伊人久久大香线蕉| av网站免费在线观看视频| 18禁动态无遮挡网站| 欧美精品国产亚洲| 免费大片黄手机在线观看| 久久久色成人| 日韩中字成人| 可以在线观看毛片的网站| 嫩草影院新地址| 亚洲av日韩在线播放| 美女视频免费永久观看网站| 中文精品一卡2卡3卡4更新| 特级一级黄色大片| 久久久色成人| 黄片wwwwww| 亚洲丝袜综合中文字幕| 少妇猛男粗大的猛烈进出视频 | 国产精品蜜桃在线观看| 少妇的逼好多水| 岛国毛片在线播放| 在线免费十八禁| 国内揄拍国产精品人妻在线| 少妇的逼水好多| videossex国产| 国产极品天堂在线| 国产日韩欧美亚洲二区| 久久人人爽av亚洲精品天堂 | 欧美精品国产亚洲| 九草在线视频观看| 国产精品99久久99久久久不卡 | 国产午夜精品一二区理论片| 一区二区三区免费毛片| 免费看a级黄色片| 丝袜喷水一区| 日本爱情动作片www.在线观看| 人妻夜夜爽99麻豆av| 成人亚洲精品av一区二区| 色哟哟·www| 一个人看的www免费观看视频| 免费黄色在线免费观看| 我的女老师完整版在线观看| 亚洲,欧美,日韩| 高清欧美精品videossex| 亚洲天堂av无毛| 日韩电影二区| 性插视频无遮挡在线免费观看| 亚洲精品日韩在线中文字幕| 精品久久久久久久久亚洲| 国产淫语在线视频| 性插视频无遮挡在线免费观看| 久久久午夜欧美精品| 最近最新中文字幕大全电影3| 色视频在线一区二区三区| 国产免费一级a男人的天堂| 国产精品蜜桃在线观看| 国国产精品蜜臀av免费| 最近中文字幕2019免费版| 欧美丝袜亚洲另类| 在线观看三级黄色| 成人免费观看视频高清| 韩国av在线不卡| 亚州av有码| www.av在线官网国产| 久久久久国产精品人妻一区二区| av免费在线看不卡| 亚洲最大成人中文| 国产精品一区www在线观看| 99久久精品热视频| 欧美人与善性xxx| 国产精品熟女久久久久浪| 日韩精品有码人妻一区| 欧美老熟妇乱子伦牲交| 在现免费观看毛片| 亚洲欧洲国产日韩| 欧美高清成人免费视频www| 春色校园在线视频观看| 日韩中字成人| 午夜免费男女啪啪视频观看| 99久久精品国产国产毛片| 日韩伦理黄色片| 97超碰精品成人国产| 热99国产精品久久久久久7| 国内揄拍国产精品人妻在线| 熟女av电影| 精品一区二区三区视频在线| 菩萨蛮人人尽说江南好唐韦庄| 国产毛片在线视频| 国产人妻一区二区三区在| 久久久久久久久大av| 亚洲最大成人中文| 亚洲精品第二区| xxx大片免费视频| 一个人看的www免费观看视频| 成年免费大片在线观看| 一边亲一边摸免费视频| 欧美区成人在线视频|