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

    土地利用變化對(duì)阜平流域的徑流影響研究*

    2014-11-24 08:35:30郝振純蘇振寬
    關(guān)鍵詞:阜平徑流量水文

    郝振純,蘇振寬,鞠 琴

    (河海大學(xué)水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210098)

    阜平流域是海河流域上游的一個(gè)重要支流,近幾十年來(lái),阜平流域流量呈現(xiàn)明顯的減小趨勢(shì)。該流域除受降雨等氣象因素影響外,人類對(duì)于自然的改造和土地利用類型的變化也對(duì)徑流產(chǎn)生了重要影響。在水文循環(huán)過(guò)程中,土地利用變化對(duì)水文過(guò)程的影響主要表現(xiàn)在對(duì)水分循環(huán)和水質(zhì)水量的改變上[1]。流域土地利用變化通過(guò)影響地表蒸散發(fā)、土壤水分狀況及地表植被的截留量等,對(duì)流域水量平衡產(chǎn)生影響。

    國(guó)內(nèi)外學(xué)者對(duì)土地利用變化的水文效應(yīng)研究主要是借助于SWAT模型。郝芳華等[2]討論了土地利用變化對(duì)產(chǎn)流產(chǎn)沙的影響。陳軍鋒[3]等研究了有植被全覆蓋與無(wú)植被覆蓋情景下徑流深與蒸發(fā)量的變化。張蕾娜[4]通過(guò)情景模擬,發(fā)現(xiàn)還草比還林更能增加徑流。索安寧等[5]對(duì)黃土高原流域的水土流失效應(yīng)做了研究。Hernandez等[6]認(rèn)為SWAT模型可以很好地反映土地覆被變化條件下的多年降水 - 徑流關(guān)系。Costa 等[7]、Fohrer等[8]和Schuol等[9]分別對(duì)托坎廷斯河流域、德國(guó)阿勒河流域和瑞典Ronnea流域土地利用變化的水文效應(yīng)做了分析。因此,本研究針對(duì)阜平流域使用SWAT模型研究土地利用變化的水文效應(yīng),為該地區(qū)的土地規(guī)劃和水資源管理提供依據(jù)。

    1 研究流域概況

    阜平流域是海河上游大清河水系的一個(gè)重要支流,包括沙河、冉莊河、青羊河、下關(guān)河獨(dú)峪河、南河溝和北流河。流域地處河北省和山西省交界處,太行山東麓,經(jīng)度 113.64-114.31°E,緯度38.78 -39.37°N,總面積 2210 km2。地勢(shì)西北高東南低,最高海拔2286 m,最低200 m。

    流域出口阜平水文站處于大清河水系沙河支流上,阜平站以上以深山區(qū)為主,山地坡度一般大于25度,溝谷多深長(zhǎng)。由于地形坡度較大,表層土壤覆蓋較少,植被較差,小片成林。流域?qū)贉貛О霛駶?rùn)半干旱地區(qū),氣候?qū)儆趤喼薮箨懶约撅L(fēng)氣候,主要特點(diǎn)為冬季干旱少雨,夏季多暴雨,降雨主要集中在6-10月份。多年平均降雨量約490 mm,年平均氣溫12.7℃,年平均相對(duì)濕度52%,年平均日照時(shí)數(shù)約2800 h。土地利用類型分為為林地、草地、居住地、水域、耕地和裸地6大類,土壤類型以棕壤和褐土為主。

    2 模型介紹

    SWAT模型是20世紀(jì)90年代由美國(guó)農(nóng)業(yè)部(USDA)的Jeff Arnold博士開(kāi)發(fā)的分布式水文模型。它具有很強(qiáng)的物理機(jī)制,能夠利用遙感和地理信息等空間信息,模擬不同土地利用、多種土地管理措施對(duì)流域水文、泥沙和化學(xué)物質(zhì)的影響[10]。模型在模擬過(guò)程中,首先將流域劃分為若干個(gè)子流域,然后根據(jù)不同土地利用類型和土壤類型,將子流域劃分為若干個(gè)水文響應(yīng)單元 (HRU),單獨(dú)計(jì)算每個(gè)HRU的徑流量,最后經(jīng)河網(wǎng)匯集得出流域的總徑流量。

    模型采用的水量平衡公式為[11]:

    式中,SWt為土壤最終含水量(mm);SW0為土壤前期含水量(mm);t為時(shí)間步長(zhǎng)(d);Rday為第i天降水量(mm);Qsurf為第i天地表徑流(mm);Ea為第i天的蒸發(fā)量(mm);Wseep為第i天存在于土壤剖面底層的滲透量和測(cè)流量(mm);Qgw為第i天地下含水量(mm)。

    模型計(jì)算的蒸散發(fā)量包括樹(shù)冠截留的水分蒸發(fā)、蒸騰和升華及土壤水的蒸發(fā)。在實(shí)際計(jì)算時(shí),首先從植被冠層截留的水分蒸發(fā)開(kāi)始計(jì)算,然后計(jì)算最大蒸騰量、最大升華量和最大土壤水分蒸發(fā)量,最后計(jì)算實(shí)際的升華量和土壤水分蒸發(fā)量[12]。

    3 模型的構(gòu)建與評(píng)價(jià)

    3.1 研究區(qū)域數(shù)據(jù)收集與處理

    3.1.1 DEM數(shù)據(jù) 地形數(shù)據(jù)數(shù)字高程模型 (DEM)來(lái)源于 SRTM 數(shù)據(jù)集 (http://strm.csi.cgiar.org/),空間分辨率為90 m。

    3.1.2 土壤數(shù)據(jù) 土壤數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院地理所數(shù)據(jù)中心提供的1∶100萬(wàn)的土壤矢量數(shù)據(jù),土壤類型為棕壤 (ZONG,22.19%)和褐土 (LING,77.81%),根據(jù)土壤類型從中國(guó)土壤數(shù)據(jù)庫(kù)查得土壤的顆粒組成以及物質(zhì)成分,使用matlab利用三次樣條插值將顆粒級(jí)配轉(zhuǎn)化為美國(guó)標(biāo)準(zhǔn)[13],并通過(guò)SPAW軟件[14]計(jì)算得出相關(guān)屬性數(shù)據(jù),建立土壤數(shù)據(jù)庫(kù)[15]。

    3.1.3 土地利用數(shù)據(jù) 土地利用類型動(dòng)態(tài)變化選取Landsat遙感衛(wèi)星圖像為數(shù)據(jù)源。利用ENVI遙感圖像處理軟件對(duì)圖像進(jìn)行地物監(jiān)督分類,在信息提取過(guò)程中參考區(qū)域?qū)n}研究結(jié)果及圖件、地形圖等相關(guān)資料,將監(jiān)督分類結(jié)果進(jìn)行人工校正,最終將土地利用類型劃分為六類:林地、草地、居住地、水域、耕地、裸地,如圖2。

    3.1.4 水文與氣象資料 水文資料采用流域出口阜平水文站的1968-2010年的逐月實(shí)測(cè)流量資料。雨量資料采用冉莊、下關(guān)、莊旺、不老臺(tái)、砂窩、龍泉關(guān)和橋南溝7個(gè)雨量站的逐日降雨資料。

    氣象資料采用繁峙氣象站的同期逐日氣象數(shù)據(jù),包括最高、最低氣溫,太陽(yáng)輻射量,降水量,相對(duì)濕度,平均氣壓,平均風(fēng)速,日照時(shí)數(shù)。太陽(yáng)輻射量由模型自帶的天氣發(fā)生器模擬得到。

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

    圖1 阜平流域土壤圖Fig.1 Soil types of the Fuping basin

    由于SWAT模型輸入?yún)?shù)較多,在校準(zhǔn)之前需對(duì)參數(shù)進(jìn)行敏感性分析。模型參數(shù)敏感性分析采用LH-OAT方法,該方法結(jié)合LH和OAT,同時(shí)又具備這兩種方法的優(yōu)點(diǎn)[16-17]。LH方法由 Mckay于1979年提出,將參數(shù)分布空間分成N個(gè),每一個(gè)范圍出現(xiàn)的概率均為1/N,然后生成參數(shù)隨機(jī)值,且每一個(gè)值域僅抽樣一次,參數(shù)隨機(jī)組合,最后對(duì)結(jié)果進(jìn)行多元線性回歸分析。OAT方法由Morris于1991年提出,模型每次運(yùn)行時(shí)只改變一個(gè)參數(shù),考慮到某些參數(shù)的靈敏度可能會(huì)受到其他參數(shù)值的選取的影響,所以模型輸入若干組參數(shù)重復(fù)運(yùn)行,最終靈敏度由其平均值決定。模型首先利用LH法進(jìn)行抽樣,然后采用OAT法進(jìn)行抽樣,確保了所有參數(shù)均被采樣,精減了需要率定的參數(shù)的個(gè)數(shù),提高了計(jì)算效率。

    圖2 土地利用變化圖Fig.2 Change of land use

    圖3 阜平流域站點(diǎn)分布圖Fig.3 Distribution of stations in the Fuping Basin

    圖4 月流量模擬值與實(shí)測(cè)值對(duì)比圖Fig.4 Comparison of simulated and measured average

    本次研究中靈敏度較高的參數(shù)為CH-K2(河道有效水力傳導(dǎo)系數(shù))、CANMX(最大冠層蓄水量)、ESCO(土壤蒸發(fā)補(bǔ)償系數(shù))、ALPHA-BF(基流ɑ系數(shù))、SLOPE(平均坡度)、SOL_K(飽和水力傳導(dǎo)系數(shù))、CN2(SCS徑流曲線系數(shù))、GWQMN(淺層地下水徑流系數(shù))、GW_REVAP(地下水再蒸發(fā)系數(shù))、SOL_AWC(土壤可利用水量)。土壤類型與地形坡度反映的是流域基礎(chǔ)地質(zhì)和地理背景,它們的變化以百萬(wàn)年為周期,在研究過(guò)程中假定其不變。為了確定土地利用類型的改變對(duì)流域徑流的影響,所以在后面的模擬過(guò)程中僅改變土地利用資料,其他參數(shù)和數(shù)據(jù)均保持不變。

    選取1968-1969年為模型緩沖期,1970-1995年為校準(zhǔn)期,1996-2010年為驗(yàn)證期。模型模擬結(jié)果采用Nash-Sutcliffe確定性系數(shù)Ens、決定系數(shù)R2來(lái)評(píng)價(jià)[18]。計(jì)算公式如下:

    式中,Qobs,i為實(shí)測(cè)數(shù)據(jù),Qsim,i為模擬數(shù)據(jù),為實(shí)測(cè)數(shù)據(jù)平均值,為模擬數(shù)據(jù)平均值,n為數(shù)據(jù)個(gè)數(shù)。

    模型校準(zhǔn)和驗(yàn)證結(jié)果見(jiàn)圖4,評(píng)價(jià)指標(biāo)見(jiàn)表1

    表1 阜平水文站月徑流模擬結(jié)果評(píng)價(jià)Table 1 Evaluation of monthly runoff simulated results of the Fuping station

    綜合這兩種指標(biāo),通常認(rèn)為Ens>0.50,R2>0.65,模擬結(jié)果令人滿意[18]。根據(jù)模型模擬結(jié)果可知SWAT模型適用于阜平流域。

    4 土地利用變化影響

    4.1 降水量與徑流量變化關(guān)系

    在研究流域土地利用變化引起的水文效應(yīng)之前,首先要確定氣象因素主要是降水量的變化對(duì)其的影響,因此需對(duì)流域多年的降水量與徑流量變化趨勢(shì)進(jìn)行分析,結(jié)果如圖5和6所示。

    圖5 流域多年降水量變化圖Fig.5 Change of precipitation for years

    根據(jù)對(duì)流域多年降水量的5年滑動(dòng)平均值變化可知,其呈現(xiàn)不規(guī)則周期性變化,21世紀(jì)之前降水量變化幅度較大,在最近的10年間,變化幅度較小,多年降水量呈微弱下降趨勢(shì),年平均變化率為-0.25 mm/a,下降趨勢(shì)不明顯。

    圖6 流域多年徑流量變化圖Fig.6 Change of runoff for years

    對(duì)流域多年徑流量趨勢(shì)分析可知,其相應(yīng)于降水量變化也呈現(xiàn)不規(guī)則周期變化,但不同的是,多年徑流量下降趨勢(shì)非常明顯,年變化率為-1.97 mm/a。

    4.2 土地利用類型變化對(duì)流域徑流的影響

    為了定量分析土地利用變化對(duì)流域徑流影響過(guò)程,以校準(zhǔn)好的參數(shù)分別對(duì)模型輸入1974、1989、1993、2000、2001、2004、2006年土地利用資料設(shè)定7種情景模擬徑流。土地利用變化如圖7。

    圖7 土地利用面積變化圖Fig.7 Change of percent of land use area

    根據(jù)圖7可知,該流域土地利用主要變化趨勢(shì)是草地減少,林地、裸地、耕地和居住地增加。土地利用類型以林地和草地為主,總面積占流域90%左右。從20世紀(jì)90年代開(kāi)始,林地面積逐漸增多,草地面積在不斷減少。土地利用類型的轉(zhuǎn)化主要是草地向林地的大面積轉(zhuǎn)化。

    根據(jù)以上7種情景,分別模擬了阜平流域1970-2010年的年徑流量與蒸散發(fā)量變化,見(jiàn)圖8、9。結(jié)果表明,土地利用變化對(duì)徑流影響比較顯著。從不同土地利用時(shí)期資料模擬出的多年平均徑流量可知,7種情景下的多年平均徑流量呈現(xiàn)明顯減小的趨勢(shì),相應(yīng)的多年平均蒸散發(fā)量呈現(xiàn)增大的趨勢(shì),也就是說(shuō),隨著草地向林地的面積轉(zhuǎn)化,流域徑流減少,蒸散發(fā)量增大。

    圖8 不同年份土地利用情景下年徑流量對(duì)比圖Fig.8 Comparison of values of runoff in various year-dependent scenarios

    圖9 多年平均蒸散發(fā)量與徑流量模擬值Fig.9 Values of average annual evapotranspiration and simulated runoff

    1974年土地利用資料中,草地面積1154.32 km2,在7種情景中最大,林地面積最小922.98 km2,模擬出的多年平均徑流量最大,達(dá)到90.42 mm,實(shí)際蒸散發(fā)量最小,為389.89 mm。主要是由于林地的截留量與蒸散發(fā)量要比草地大得多。而居住地、水域、耕地和裸地所占面積較小,總體上對(duì)徑流影響不大。從土地利用變化來(lái)看,1989年和1993年之間變化不大,模擬的多年平均徑流量和蒸散發(fā)量變化也不大。1993年和2000年相比,林地面積增加了4.7%,草地面積減少了6.4%,面積變化相對(duì)較為明顯,但由于裸地面積增加了55.7%,使得模型輸入2000年土地利用資料模擬出的多年平均徑流量與前兩期相比,變化并不是特別大。2004年與2006年土地利用變化不大,模擬的徑流結(jié)果差別也相對(duì)不大。從圖7可知,1993年和2004年的土地利用資料變化最為明顯,變化量見(jiàn)表2,模擬結(jié)果見(jiàn)表3??梢?jiàn)居住地、耕地和裸地變化幅度比較大,但是所占的總面積依然很小,對(duì)徑流影響不如林地和草地影響大。與實(shí)測(cè)數(shù)據(jù)對(duì)比可知,草地面積最小,林地面積最大的2004年土地利用資料模擬出的徑流量減少12.38%,與1993年相比減少6.46%。

    表2 1993和2004年土地利用面積變化Table 2 Change of land use area in 1993 and 2004,respectively

    表3 1993和2004年徑流模擬結(jié)果Table 3 Results of simulated runoff in 1993 and 2004

    5 結(jié)論

    利用流域出口阜平水文站43 a的水文實(shí)測(cè)數(shù)據(jù)和流域7個(gè)雨量站43 a的降水?dāng)?shù)據(jù),利用7種土地利用情景,建立SWAT模型,模擬了不同情景下徑流的響應(yīng)。

    1)通過(guò)對(duì)模型的校準(zhǔn)和驗(yàn)證,校準(zhǔn)期Ens=0.88,R2=0.93,驗(yàn)證期 Ens=0.80,R2=0.83,說(shuō)明SWAT模型適合阜平流域的徑流模擬。

    2)從20世紀(jì)90年代開(kāi)始,該流域土地利用類型發(fā)生明顯變化,不同情景下模擬的流域出口徑流量明顯減少。

    3)阜平流域土地利用類型以林地和草地為主,流域徑流量呈現(xiàn)明顯的減少趨勢(shì),降水量下降趨勢(shì)不明顯,林地和草地的變化是徑流量變化的重要因素。

    4)林地和草地的變化主要是通過(guò)影響蒸散發(fā)量改變水文循環(huán)過(guò)程,從而影響徑流量的變化。不同情景下模擬多年平均蒸發(fā)量變化趨勢(shì)與徑流量變化相反。

    由此可見(jiàn),不合理的土地利用會(huì)導(dǎo)致一定的生態(tài)平衡失調(diào),在土地利用規(guī)劃過(guò)程中,要兼顧水文效應(yīng)與經(jīng)濟(jì)效益。另一方面要充分認(rèn)識(shí)人類活動(dòng)引起的土地利用變化對(duì)徑流的影響。

    [1]史培軍,宮鵬,李曉兵,等.土地利用/覆被變化研究的方法與實(shí)踐[M].北京:科學(xué)出版社,2000.

    [2]郝芳華,陳利群,劉昌明,等.土地利用變化對(duì)產(chǎn)流和產(chǎn)沙的影響分析[J].水土保持學(xué)報(bào),2004,18(3):5 -8.

    [3]陳軍鋒,李秀彬.土地覆被變化的水文響應(yīng)模擬研究[J].應(yīng)用生態(tài)學(xué)報(bào),2004,15(5):833 -836.

    [4]張蕾娜.白河流域土地覆被變化水文效應(yīng)的分析與模擬[D].北京:中國(guó)科學(xué)院地理科學(xué)與資源研究所,2004.

    [5]索安寧,李金朝,王天明,等.黃土高原流域土地利用變化的水土流失效應(yīng)[J].水利學(xué)報(bào),2008,39(7):767-772.

    [6]HERNANDEZ M,MILLER S N,GOODRIOH D C,et al.Modeling runoff response to land cover and rainfall spatial variability in semi-arid watersheds[J].Environmental Monitoring and Assessment,2000,64:285 -298.

    [7]COSTA M H,BOTTA A,CARDILLE J A.Effects of large-scale changes in land cover on the discharge of Tocantins River,Southeastern Amazonia[J].Journal of Hydrology,2003,283(1):206 -217.

    [8]FOHRER N,M?LLER D,STEINER N.An interdisciplinary modelling approach to evaluate the effects of land use change[J].Physics and Chemistry of the Earth,Parts A/B/C,2002,27(9):655 -662.

    [9]SCHUOL J,ABBASPOUR K C,SRINIVASAN R,et al.Estimation of freshwater availability in the West African sub-continent using the SWAT hydrologic model[J].Journal of Hydrology,2008,352(1/2):30 -49.

    [10]王中根,劉昌明,黃友波.SWAT模型的原理、結(jié)構(gòu)及應(yīng)用研究[J].地理科學(xué)進(jìn)展,2003,22(1):79-86.

    [11]NEITSCH S L,AONOLD J G,KINIRY J R,et al.Soil and water assessment tool theoretical documentation,version 2005 [EB/OL ].[2011-04-12 ]. http://www.brc.tamus.edu/swat/downloads/doc/swat2005/SWAT2005 theory final.pdf,2005.

    [12]肖軍倉(cāng),周文斌,羅定貴,等.非點(diǎn)源污染模型 -SWAT用戶應(yīng)用指南[M].北京:地質(zhì)出版社,2010.

    [13]蔡永明,張科利,李雙才.不同粒徑制間土壤質(zhì)地資料的轉(zhuǎn)換問(wèn)題研究[J].土壤學(xué)報(bào),2003,40(4):511-517.

    [14]SAXTON K E,RAWLS W J.Soil water characteristic estimates by texture and organic matter for hydrologic solutions[J].Soil Scence Society of America Journal,2006,70(5):1569 -1578.

    [15]魏懷斌,張占龐,楊金鵬.SWAT模型土壤數(shù)據(jù)庫(kù)建立方法[J].水利水電技術(shù),2007,38(6):15 -18.

    [16]Van GRIENSVEN.Sensitivity auto-calibration uncertainty and model evaluation in SWAT2005[Z].Grassland,soil and water research service,Temple,TX,2007.

    [17]McKAY M D,BECKMAN R J,CONOVER W J.A comparison of three methods for selecting values of input variables in the analysis of output from a computer code[J].Technometeics,1979,21(2):239 - 245.

    [18]MORIASI D N,ARNOLD J G,Van LIEW M W,et al.Model evaluation guidelines for systematic quantification of accuracy in watershed simulations[J].Transactions of the ASABE,2007,50(3):885 -900.

    猜你喜歡
    阜平徑流量水文
    2022年《中國(guó)水文年報(bào)》發(fā)布
    我家的小貓 咪咪
    阜平輝綠巖飾面石材礦廢石利用簡(jiǎn)析
    訪河北阜平農(nóng)家
    水文
    水文水資源管理
    水文比擬法在計(jì)算河川徑流量時(shí)的修正
    水文
    SCS模型在紅壤土坡地降雨徑流量估算中的應(yīng)用
    資江流域徑流量演變規(guī)律研究
    午夜免费激情av| 男插女下体视频免费在线播放| 国产日本99.免费观看| 老熟妇乱子伦视频在线观看| 俺也久久电影网| 狂野欧美激情性xxxx在线观看| 欧美日本视频| 国产精品日韩av在线免费观看| 亚洲成人久久性| av天堂中文字幕网| 日韩大尺度精品在线看网址| 精品人妻一区二区三区麻豆 | 欧美丝袜亚洲另类 | 免费高清视频大片| 少妇人妻一区二区三区视频| 如何舔出高潮| 露出奶头的视频| 久久精品国产清高在天天线| av.在线天堂| 97碰自拍视频| 国产一区二区在线观看日韩| 国内揄拍国产精品人妻在线| 看免费成人av毛片| 很黄的视频免费| 免费看美女性在线毛片视频| 欧美精品国产亚洲| 国产av在哪里看| 波多野结衣高清作品| 欧美一区二区国产精品久久精品| 国产一区二区在线观看日韩| 免费一级毛片在线播放高清视频| x7x7x7水蜜桃| 黄色视频,在线免费观看| a级毛片a级免费在线| 色精品久久人妻99蜜桃| 久久久国产成人免费| 一个人观看的视频www高清免费观看| 草草在线视频免费看| 午夜福利在线观看免费完整高清在 | 麻豆一二三区av精品| 人妻制服诱惑在线中文字幕| 午夜久久久久精精品| 日本a在线网址| 一a级毛片在线观看| 亚洲专区国产一区二区| 亚洲欧美日韩东京热| 在线播放无遮挡| 国产精品不卡视频一区二区| 国产亚洲精品久久久com| 国内揄拍国产精品人妻在线| 亚洲 国产 在线| 亚洲av五月六月丁香网| 日韩欧美国产一区二区入口| 搡女人真爽免费视频火全软件 | 欧美精品国产亚洲| 12—13女人毛片做爰片一| 在线观看66精品国产| 久久久久久久久中文| 亚洲,欧美,日韩| 嫩草影视91久久| 久久久久久国产a免费观看| 97热精品久久久久久| 亚洲人与动物交配视频| 桃色一区二区三区在线观看| 国产精品爽爽va在线观看网站| 婷婷色综合大香蕉| 免费看a级黄色片| 日本五十路高清| 亚洲一级一片aⅴ在线观看| 两个人视频免费观看高清| 日本欧美国产在线视频| 国产人妻一区二区三区在| 精品一区二区免费观看| 色在线成人网| 欧美区成人在线视频| 精品一区二区三区视频在线观看免费| 欧美日韩中文字幕国产精品一区二区三区| 蜜桃亚洲精品一区二区三区| 成人欧美大片| 国产精品乱码一区二三区的特点| 中文字幕免费在线视频6| 狂野欧美激情性xxxx在线观看| 最新中文字幕久久久久| 小说图片视频综合网站| 熟女电影av网| 亚洲av二区三区四区| 内地一区二区视频在线| 男人和女人高潮做爰伦理| 99riav亚洲国产免费| 在线观看午夜福利视频| av黄色大香蕉| 黄色日韩在线| 99久久久亚洲精品蜜臀av| 麻豆成人av在线观看| 久久久精品大字幕| 亚洲av不卡在线观看| 韩国av在线不卡| 日韩亚洲欧美综合| 成人亚洲精品av一区二区| 十八禁国产超污无遮挡网站| 欧美高清成人免费视频www| netflix在线观看网站| 老熟妇乱子伦视频在线观看| 国产午夜精品论理片| 国产午夜福利久久久久久| 一进一出好大好爽视频| 久久久久久伊人网av| 午夜福利高清视频| 国产欧美日韩一区二区精品| 免费av观看视频| 欧美+亚洲+日韩+国产| 欧美性猛交黑人性爽| 熟妇人妻久久中文字幕3abv| 无遮挡黄片免费观看| 亚洲性夜色夜夜综合| 久久久国产成人免费| 十八禁国产超污无遮挡网站| 在线观看av片永久免费下载| 亚洲最大成人av| 欧美xxxx性猛交bbbb| 赤兔流量卡办理| 嫩草影视91久久| 乱人视频在线观看| 精品久久久久久成人av| 日本黄色视频三级网站网址| 免费高清视频大片| 欧美色欧美亚洲另类二区| 日韩欧美在线乱码| 国产伦精品一区二区三区视频9| 天天躁日日操中文字幕| 日韩欧美国产在线观看| 男人舔奶头视频| 99九九线精品视频在线观看视频| 日韩欧美在线二视频| 国语自产精品视频在线第100页| 国产蜜桃级精品一区二区三区| 99riav亚洲国产免费| 亚洲在线自拍视频| 人人妻人人澡欧美一区二区| 色精品久久人妻99蜜桃| 可以在线观看毛片的网站| 97碰自拍视频| 毛片一级片免费看久久久久 | 嫩草影视91久久| 我要搜黄色片| 乱码一卡2卡4卡精品| 999久久久精品免费观看国产| 床上黄色一级片| 麻豆av噜噜一区二区三区| 亚洲欧美日韩东京热| 精品国产三级普通话版| 国产淫片久久久久久久久| 小蜜桃在线观看免费完整版高清| 亚洲精品日韩av片在线观看| 禁无遮挡网站| 色5月婷婷丁香| eeuss影院久久| 国产成人影院久久av| 联通29元200g的流量卡| 亚洲精品成人久久久久久| 嫩草影院入口| a级一级毛片免费在线观看| 精品人妻1区二区| 制服丝袜大香蕉在线| 精品99又大又爽又粗少妇毛片 | 亚洲成人久久爱视频| 日日摸夜夜添夜夜添av毛片 | 亚洲在线自拍视频| av国产免费在线观看| 久久午夜福利片| 亚洲在线观看片| 97碰自拍视频| 嫁个100分男人电影在线观看| 国产高清三级在线| 村上凉子中文字幕在线| 免费观看的影片在线观看| 日日摸夜夜添夜夜添小说| 国产男靠女视频免费网站| 岛国在线免费视频观看| 日日夜夜操网爽| 亚洲美女黄片视频| 美女cb高潮喷水在线观看| 免费人成在线观看视频色| 最好的美女福利视频网| 男人的好看免费观看在线视频| 黄色欧美视频在线观看| 亚洲成人中文字幕在线播放| 亚洲va日本ⅴa欧美va伊人久久| 欧美bdsm另类| 能在线免费观看的黄片| 欧美最新免费一区二区三区| 美女xxoo啪啪120秒动态图| 国内久久婷婷六月综合欲色啪| 亚洲一区二区三区色噜噜| 中文字幕熟女人妻在线| 国产av一区在线观看免费| 国产亚洲精品久久久久久毛片| АⅤ资源中文在线天堂| 内地一区二区视频在线| 久久精品人妻少妇| 中亚洲国语对白在线视频| 免费观看精品视频网站| avwww免费| 精品欧美国产一区二区三| 赤兔流量卡办理| 日本 av在线| 欧美色视频一区免费| 中文资源天堂在线| 国产精品一区二区三区四区久久| 特大巨黑吊av在线直播| 精品人妻偷拍中文字幕| 一个人观看的视频www高清免费观看| 久久久精品大字幕| 最新中文字幕久久久久| x7x7x7水蜜桃| 国产极品精品免费视频能看的| 国产精品自产拍在线观看55亚洲| 成人综合一区亚洲| 国产高清激情床上av| av女优亚洲男人天堂| 内射极品少妇av片p| 99热网站在线观看| 久久久国产成人精品二区| 国产精品亚洲一级av第二区| 观看美女的网站| 真实男女啪啪啪动态图| 搡老岳熟女国产| 1024手机看黄色片| 欧美最新免费一区二区三区| 搡老妇女老女人老熟妇| 中文字幕高清在线视频| 亚洲一级一片aⅴ在线观看| 亚洲成a人片在线一区二区| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区亚洲精品在线观看| 级片在线观看| 精华霜和精华液先用哪个| 国产淫片久久久久久久久| 国产精品国产高清国产av| 国产精品日韩av在线免费观看| 亚洲久久久久久中文字幕| 两个人的视频大全免费| 一区二区三区免费毛片| 老司机深夜福利视频在线观看| 十八禁网站免费在线| 成人二区视频| 大型黄色视频在线免费观看| 欧美成人一区二区免费高清观看| 国产午夜精品论理片| 亚洲国产欧洲综合997久久,| 九九爱精品视频在线观看| ponron亚洲| videossex国产| 91狼人影院| 亚洲精品色激情综合| 老司机午夜福利在线观看视频| 动漫黄色视频在线观看| 夜夜爽天天搞| 成人永久免费在线观看视频| 色哟哟·www| 午夜久久久久精精品| 日本与韩国留学比较| 春色校园在线视频观看| 精品福利观看| 我要搜黄色片| 韩国av一区二区三区四区| 午夜福利欧美成人| 午夜免费激情av| 国产成人a区在线观看| 女同久久另类99精品国产91| 最新中文字幕久久久久| 日日干狠狠操夜夜爽| 欧美高清性xxxxhd video| 人人妻,人人澡人人爽秒播| 桃红色精品国产亚洲av| 国产高清不卡午夜福利| 狠狠狠狠99中文字幕| 日日夜夜操网爽| 香蕉av资源在线| 干丝袜人妻中文字幕| 12—13女人毛片做爰片一| 少妇猛男粗大的猛烈进出视频 | 99精品在免费线老司机午夜| 色吧在线观看| 亚洲avbb在线观看| 日韩中字成人| 免费看av在线观看网站| 日日啪夜夜撸| 亚洲成人久久爱视频| 日韩中文字幕欧美一区二区| 日本 欧美在线| 久久久久免费精品人妻一区二区| 国产精品爽爽va在线观看网站| 淫妇啪啪啪对白视频| 最新中文字幕久久久久| 欧美成人免费av一区二区三区| 神马国产精品三级电影在线观看| 能在线免费观看的黄片| 国产高清视频在线播放一区| 免费观看的影片在线观看| 亚洲色图av天堂| 国产伦人伦偷精品视频| 熟妇人妻久久中文字幕3abv| 免费大片18禁| 在线播放无遮挡| 亚洲国产欧洲综合997久久,| 99久国产av精品| 亚洲色图av天堂| 波多野结衣高清作品| 免费看日本二区| 国产乱人视频| 午夜精品久久久久久毛片777| 熟女电影av网| 中亚洲国语对白在线视频| 天天一区二区日本电影三级| 嫩草影院新地址| 18+在线观看网站| 国产精品一区二区三区四区久久| 精品久久久久久久久久免费视频| 亚洲av免费在线观看| 欧美高清性xxxxhd video| 亚洲国产精品sss在线观看| x7x7x7水蜜桃| 精品人妻1区二区| 国产美女午夜福利| 久久精品人妻少妇| 国产中年淑女户外野战色| 亚洲一区二区三区色噜噜| 午夜老司机福利剧场| 亚洲色图av天堂| 欧美成人性av电影在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国内毛片毛片毛片毛片毛片| 中国美白少妇内射xxxbb| 久久精品综合一区二区三区| 天堂av国产一区二区熟女人妻| 亚州av有码| 日本三级黄在线观看| 亚洲精品日韩av片在线观看| 如何舔出高潮| 久久精品久久久久久噜噜老黄 | 精品午夜福利在线看| 午夜久久久久精精品| 久久99热这里只有精品18| 中文字幕免费在线视频6| 97人妻精品一区二区三区麻豆| 免费在线观看成人毛片| 欧美又色又爽又黄视频| 久久99热6这里只有精品| 国产高清视频在线播放一区| 一进一出抽搐动态| 亚洲 国产 在线| 日韩 亚洲 欧美在线| 精品一区二区三区人妻视频| 国产精品嫩草影院av在线观看 | av在线观看视频网站免费| 在线免费观看不下载黄p国产 | 中文字幕免费在线视频6| 精华霜和精华液先用哪个| 成人鲁丝片一二三区免费| 国产三级在线视频| 欧美丝袜亚洲另类 | 91精品国产九色| 自拍偷自拍亚洲精品老妇| 热99re8久久精品国产| 国产高清视频在线播放一区| 国产成人av教育| 嫩草影视91久久| 国产成人av教育| 91午夜精品亚洲一区二区三区 | 一a级毛片在线观看| 国产精品久久久久久久久免| 3wmmmm亚洲av在线观看| 亚洲图色成人| 我的女老师完整版在线观看| 黄色欧美视频在线观看| 欧美一区二区国产精品久久精品| 国产精品人妻久久久久久| 亚洲最大成人av| 日韩在线高清观看一区二区三区 | 免费人成在线观看视频色| 哪里可以看免费的av片| 日日摸夜夜添夜夜添小说| 琪琪午夜伦伦电影理论片6080| 国产av麻豆久久久久久久| 天堂√8在线中文| 亚洲欧美日韩高清在线视频| 丰满乱子伦码专区| 亚洲成av人片在线播放无| a级毛片免费高清观看在线播放| 日韩一区二区视频免费看| 少妇的逼好多水| a级一级毛片免费在线观看| 欧美一区二区国产精品久久精品| 韩国av一区二区三区四区| 国内毛片毛片毛片毛片毛片| 成人国产综合亚洲| a在线观看视频网站| 动漫黄色视频在线观看| 国产精品美女特级片免费视频播放器| 亚洲真实伦在线观看| 夜夜夜夜夜久久久久| 最好的美女福利视频网| 色av中文字幕| 国产免费av片在线观看野外av| 国产精品永久免费网站| 亚洲一级一片aⅴ在线观看| 久久久精品大字幕| a级毛片免费高清观看在线播放| 国产亚洲精品综合一区在线观看| 精品免费久久久久久久清纯| 国产乱人伦免费视频| 女人十人毛片免费观看3o分钟| 又黄又爽又免费观看的视频| 国产精品久久久久久av不卡| 日本-黄色视频高清免费观看| 99在线视频只有这里精品首页| 色哟哟哟哟哟哟| 国产精品美女特级片免费视频播放器| 2021天堂中文幕一二区在线观| 一级av片app| 男人舔女人下体高潮全视频| 最近最新中文字幕大全电影3| 国产精品国产三级国产av玫瑰| 中文字幕久久专区| or卡值多少钱| 两性午夜刺激爽爽歪歪视频在线观看| 久久天躁狠狠躁夜夜2o2o| 麻豆成人av在线观看| 日本欧美国产在线视频| 久久热精品热| 能在线免费观看的黄片| 亚洲中文字幕一区二区三区有码在线看| 中亚洲国语对白在线视频| 亚洲人成网站在线播放欧美日韩| 久久6这里有精品| a级一级毛片免费在线观看| 国产精品自产拍在线观看55亚洲| 99久久九九国产精品国产免费| 日韩欧美精品免费久久| 男人的好看免费观看在线视频| 亚洲中文字幕一区二区三区有码在线看| 国产人妻一区二区三区在| av在线观看视频网站免费| 人妻久久中文字幕网| 精品国产三级普通话版| 黄色欧美视频在线观看| av在线天堂中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久国内视频| 亚洲专区中文字幕在线| 99久国产av精品| 免费av毛片视频| 国产成人a区在线观看| 丰满乱子伦码专区| 国产毛片a区久久久久| 亚洲欧美清纯卡通| 亚洲专区中文字幕在线| 欧美成人免费av一区二区三区| 99热网站在线观看| 国产在线男女| 一夜夜www| 精品不卡国产一区二区三区| 久久久国产成人免费| 亚洲精品一区av在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产91精品成人一区二区三区| 人妻丰满熟妇av一区二区三区| 亚洲狠狠婷婷综合久久图片| 欧美激情久久久久久爽电影| 亚洲欧美清纯卡通| 午夜免费成人在线视频| 别揉我奶头 嗯啊视频| 亚洲18禁久久av| 国产精品不卡视频一区二区| 亚洲精品影视一区二区三区av| av.在线天堂| 国内精品久久久久久久电影| 国产精品av视频在线免费观看| av在线老鸭窝| 淫秽高清视频在线观看| 黄色一级大片看看| 麻豆国产av国片精品| 久久人人精品亚洲av| 校园春色视频在线观看| 国产精品,欧美在线| av专区在线播放| 我要搜黄色片| 成人特级黄色片久久久久久久| 久99久视频精品免费| 亚洲真实伦在线观看| 亚洲精品乱码久久久v下载方式| 亚洲av美国av| 91在线观看av| 99国产精品一区二区蜜桃av| 国产免费一级a男人的天堂| 国产69精品久久久久777片| 午夜a级毛片| 大又大粗又爽又黄少妇毛片口| 我要搜黄色片| 动漫黄色视频在线观看| 性欧美人与动物交配| 我的女老师完整版在线观看| 91在线观看av| www.www免费av| 日韩人妻高清精品专区| 国产精品久久电影中文字幕| 国产一区二区在线av高清观看| 成人高潮视频无遮挡免费网站| a在线观看视频网站| 亚洲精品在线观看二区| 日本在线视频免费播放| 综合色av麻豆| 精品福利观看| 国产麻豆成人av免费视频| 一级a爱片免费观看的视频| 亚洲人成网站高清观看| 欧美性感艳星| 在线a可以看的网站| 精品人妻熟女av久视频| 热99在线观看视频| 成人欧美大片| 亚洲av熟女| 12—13女人毛片做爰片一| 99久国产av精品| 波多野结衣高清作品| 国产探花极品一区二区| 日本一二三区视频观看| 国产男靠女视频免费网站| 日韩欧美国产在线观看| 18禁在线播放成人免费| 久久精品91蜜桃| 麻豆精品久久久久久蜜桃| 国产高清激情床上av| 午夜视频国产福利| 啪啪无遮挡十八禁网站| 久久久久久九九精品二区国产| 亚洲人成网站在线播| 亚洲国产欧美人成| 免费av不卡在线播放| 一级黄片播放器| 成人国产麻豆网| 小蜜桃在线观看免费完整版高清| 精品久久久久久久末码| 午夜免费成人在线视频| 九色国产91popny在线| 在线观看美女被高潮喷水网站| 熟女电影av网| 男人狂女人下面高潮的视频| 久久人人爽人人爽人人片va| 国产免费男女视频| 久久久国产成人免费| 久久精品久久久久久噜噜老黄 | 欧美黑人巨大hd| 12—13女人毛片做爰片一| 亚洲avbb在线观看| 亚洲精品国产成人久久av| 变态另类成人亚洲欧美熟女| 两人在一起打扑克的视频| 亚州av有码| 亚洲最大成人中文| 熟女人妻精品中文字幕| 美女大奶头视频| 久久久精品大字幕| 日韩av在线大香蕉| 亚洲三级黄色毛片| 亚洲aⅴ乱码一区二区在线播放| 在线观看av片永久免费下载| 99热这里只有精品一区| 国产一区二区亚洲精品在线观看| av福利片在线观看| 亚洲欧美激情综合另类| 亚洲av美国av| 欧美激情久久久久久爽电影| 无人区码免费观看不卡| 精品一区二区免费观看| netflix在线观看网站| 久久国内精品自在自线图片| 精品久久久久久久久亚洲 | 日本免费a在线| 国产午夜福利久久久久久| 午夜福利成人在线免费观看| 91久久精品国产一区二区三区| 一级黄色大片毛片| xxxwww97欧美| 亚洲av成人av| 黄色配什么色好看| 国产中年淑女户外野战色| 国产精品国产三级国产av玫瑰| 亚洲最大成人中文| 无人区码免费观看不卡| 亚洲一级一片aⅴ在线观看| 天天一区二区日本电影三级| 日韩 亚洲 欧美在线| 国产私拍福利视频在线观看| 日本熟妇午夜| 热99re8久久精品国产| 精品乱码久久久久久99久播| 不卡视频在线观看欧美| 亚洲精华国产精华精| 亚洲成人精品中文字幕电影| 亚洲欧美日韩高清在线视频| 欧美成人a在线观看| 村上凉子中文字幕在线| 可以在线观看毛片的网站| 久久久久国内视频| 婷婷色综合大香蕉| 欧美高清成人免费视频www| 久久午夜亚洲精品久久| 一a级毛片在线观看| 亚洲国产日韩欧美精品在线观看| 亚洲四区av| 搡女人真爽免费视频火全软件 | 日本成人三级电影网站| 97超级碰碰碰精品色视频在线观看|