• <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ī)律研究
    国产男靠女视频免费网站| 精品日产1卡2卡| 国产免费av片在线观看野外av| 一区二区三区国产精品乱码| av片东京热男人的天堂| 少妇人妻一区二区三区视频| 国产精品永久免费网站| 十八禁网站免费在线| 老司机在亚洲福利影院| 看免费av毛片| 久9热在线精品视频| 人人妻,人人澡人人爽秒播| 国产精品一区二区免费欧美| 欧美在线一区亚洲| 女人高潮潮喷娇喘18禁视频| 九九久久精品国产亚洲av麻豆 | 全区人妻精品视频| 成年女人毛片免费观看观看9| 视频区欧美日本亚洲| 国产99白浆流出| 亚洲第一电影网av| 国产成人精品久久二区二区91| 亚洲真实伦在线观看| 亚洲av电影不卡..在线观看| 精品久久久久久久毛片微露脸| 老汉色∧v一级毛片| 久久午夜综合久久蜜桃| 91久久精品国产一区二区成人 | 国产亚洲av高清不卡| 亚洲天堂国产精品一区在线| 别揉我奶头~嗯~啊~动态视频| 国产精品精品国产色婷婷| 99热只有精品国产| 亚洲国产欧美网| 久久久久九九精品影院| 久久久久久国产a免费观看| 国产真人三级小视频在线观看| 国产精品电影一区二区三区| 美女 人体艺术 gogo| 日韩欧美 国产精品| 欧美丝袜亚洲另类 | 婷婷精品国产亚洲av| 亚洲avbb在线观看| 免费大片18禁| 五月伊人婷婷丁香| 黑人巨大精品欧美一区二区mp4| www国产在线视频色| 精品国产乱码久久久久久男人| 女人高潮潮喷娇喘18禁视频| 天天一区二区日本电影三级| 97超级碰碰碰精品色视频在线观看| 国产一区二区三区在线臀色熟女| 亚洲欧美日韩无卡精品| 色老头精品视频在线观看| 中文字幕最新亚洲高清| 亚洲性夜色夜夜综合| 国产三级中文精品| 白带黄色成豆腐渣| 天堂av国产一区二区熟女人妻| 成人无遮挡网站| 琪琪午夜伦伦电影理论片6080| 99久久久亚洲精品蜜臀av| 真人做人爱边吃奶动态| 久久精品国产99精品国产亚洲性色| 免费无遮挡裸体视频| 91麻豆av在线| 天天躁日日操中文字幕| 久久这里只有精品中国| 全区人妻精品视频| 国产成人av激情在线播放| 日本 欧美在线| 久久精品国产亚洲av香蕉五月| 成年版毛片免费区| 精品久久久久久成人av| а√天堂www在线а√下载| 五月伊人婷婷丁香| 欧美不卡视频在线免费观看| 欧美一级毛片孕妇| 99国产综合亚洲精品| 2021天堂中文幕一二区在线观| 小蜜桃在线观看免费完整版高清| 99热6这里只有精品| 国产一区二区三区在线臀色熟女| 亚洲精品一区av在线观看| 国产午夜精品论理片| 国产精品美女特级片免费视频播放器 | 亚洲国产精品久久男人天堂| or卡值多少钱| 色吧在线观看| 午夜福利欧美成人| 757午夜福利合集在线观看| 亚洲人成伊人成综合网2020| 亚洲精品久久国产高清桃花| 成在线人永久免费视频| 亚洲美女黄片视频| 亚洲成人免费电影在线观看| 精品国产乱子伦一区二区三区| 日本撒尿小便嘘嘘汇集6| 欧美xxxx黑人xx丫x性爽| 亚洲性夜色夜夜综合| 999久久久国产精品视频| 欧美精品啪啪一区二区三区| 亚洲第一欧美日韩一区二区三区| 最近最新中文字幕大全免费视频| 国产免费av片在线观看野外av| 亚洲中文字幕一区二区三区有码在线看 | 中亚洲国语对白在线视频| 欧美午夜高清在线| 日韩欧美 国产精品| 国产精品亚洲美女久久久| 精品日产1卡2卡| 午夜福利视频1000在线观看| 麻豆国产97在线/欧美| 亚洲自拍偷在线| 小说图片视频综合网站| 精品久久久久久成人av| 一级作爱视频免费观看| 一本久久中文字幕| 亚洲av第一区精品v没综合| 亚洲av电影在线进入| 88av欧美| 亚洲专区中文字幕在线| 人人妻人人看人人澡| 九色国产91popny在线| 亚洲国产中文字幕在线视频| 99在线视频只有这里精品首页| 免费观看人在逋| 人人妻,人人澡人人爽秒播| 国产视频一区二区在线看| av黄色大香蕉| 精品久久蜜臀av无| 欧美日韩一级在线毛片| a级毛片在线看网站| 欧美日本亚洲视频在线播放| 免费无遮挡裸体视频| 淫妇啪啪啪对白视频| 亚洲人成伊人成综合网2020| 日本 欧美在线| 亚洲自拍偷在线| 国产精品乱码一区二三区的特点| 国产成人av激情在线播放| 亚洲国产精品久久男人天堂| 欧美3d第一页| 别揉我奶头~嗯~啊~动态视频| 亚洲,欧美精品.| 成人国产综合亚洲| 亚洲午夜理论影院| 日日夜夜操网爽| 99热这里只有精品一区 | 美女 人体艺术 gogo| 桃红色精品国产亚洲av| 午夜福利视频1000在线观看| 最近最新免费中文字幕在线| 九色成人免费人妻av| 午夜福利18| 国产亚洲精品一区二区www| 校园春色视频在线观看| 天堂动漫精品| 好男人电影高清在线观看| 亚洲 欧美一区二区三区| 麻豆成人午夜福利视频| 免费大片18禁| 国内毛片毛片毛片毛片毛片| 亚洲av日韩精品久久久久久密| 99久久精品热视频| 无遮挡黄片免费观看| 成年女人看的毛片在线观看| 在线观看舔阴道视频| 国产伦精品一区二区三区四那| 午夜久久久久精精品| 国产成人欧美在线观看| 十八禁网站免费在线| 欧美日本亚洲视频在线播放| 国产午夜福利久久久久久| 欧美日韩黄片免| 亚洲欧美精品综合一区二区三区| 99久久精品一区二区三区| 少妇的丰满在线观看| 欧美成狂野欧美在线观看| 国产成人精品无人区| 美女高潮喷水抽搐中文字幕| 在线视频色国产色| 悠悠久久av| 成人国产综合亚洲| 91老司机精品| 精品欧美国产一区二区三| 国产一区二区三区在线臀色熟女| 亚洲av美国av| 久久精品91蜜桃| 成年人黄色毛片网站| 日韩中文字幕欧美一区二区| 一进一出抽搐动态| 免费看十八禁软件| 999久久久国产精品视频| 欧美最黄视频在线播放免费| aaaaa片日本免费| 国产男靠女视频免费网站| 国内精品一区二区在线观看| 亚洲 欧美一区二区三区| 国产激情久久老熟女| av女优亚洲男人天堂 | 精品一区二区三区av网在线观看| 女人被狂操c到高潮| 精品一区二区三区四区五区乱码| 少妇裸体淫交视频免费看高清| 一进一出好大好爽视频| 他把我摸到了高潮在线观看| 精品国产超薄肉色丝袜足j| 岛国在线观看网站| 国产av一区在线观看免费| 欧美中文日本在线观看视频| 久久天堂一区二区三区四区| 午夜福利在线观看免费完整高清在 | 国产精品一及| 久久人人精品亚洲av| 国产久久久一区二区三区| 香蕉国产在线看| 亚洲性夜色夜夜综合| 极品教师在线免费播放| 免费看美女性在线毛片视频| av在线天堂中文字幕| 视频区欧美日本亚洲| 亚洲精品粉嫩美女一区| 国产欧美日韩精品亚洲av| 俄罗斯特黄特色一大片| 91久久精品国产一区二区成人 | 日韩中文字幕欧美一区二区| 免费看日本二区| 精品一区二区三区视频在线 | 日韩欧美精品v在线| 亚洲avbb在线观看| 成人性生交大片免费视频hd| 色播亚洲综合网| 叶爱在线成人免费视频播放| 国产精品av视频在线免费观看| 日本五十路高清| 国产精品98久久久久久宅男小说| 人妻丰满熟妇av一区二区三区| 级片在线观看| 熟女少妇亚洲综合色aaa.| 久久九九热精品免费| 国产伦精品一区二区三区视频9 | 黄色女人牲交| 久久婷婷人人爽人人干人人爱| www日本黄色视频网| 国产伦人伦偷精品视频| 一本一本综合久久| 一a级毛片在线观看| 视频区欧美日本亚洲| 真人一进一出gif抽搐免费| av片东京热男人的天堂| 午夜久久久久精精品| 国产精品99久久久久久久久| 人人妻人人澡欧美一区二区| 91麻豆精品激情在线观看国产| 精品免费久久久久久久清纯| xxx96com| 91字幕亚洲| 性色av乱码一区二区三区2| 日本在线视频免费播放| 99国产极品粉嫩在线观看| 日韩高清综合在线| 欧美乱码精品一区二区三区| 国内揄拍国产精品人妻在线| 亚洲欧美激情综合另类| 国产亚洲av高清不卡| 中文字幕人妻丝袜一区二区| 又爽又黄无遮挡网站| 一级黄色大片毛片| 亚洲精品美女久久久久99蜜臀| 色播亚洲综合网| 一卡2卡三卡四卡精品乱码亚洲| 久久国产精品影院| 哪里可以看免费的av片| 欧美中文综合在线视频| ponron亚洲| 97超视频在线观看视频| 99久久综合精品五月天人人| 三级国产精品欧美在线观看 | 伦理电影免费视频| 成人国产一区最新在线观看| 国产亚洲精品久久久com| 亚洲av成人一区二区三| 免费高清视频大片| 色av中文字幕| 日韩中文字幕欧美一区二区| 99国产精品一区二区蜜桃av| 脱女人内裤的视频| 国产精品香港三级国产av潘金莲| 亚洲成人久久爱视频| 岛国在线观看网站| 欧美性猛交╳xxx乱大交人| 美女扒开内裤让男人捅视频| 国产伦在线观看视频一区| 成年女人永久免费观看视频| www.熟女人妻精品国产| av欧美777| 成年人黄色毛片网站| 看黄色毛片网站| 国产精品久久久人人做人人爽| 最近在线观看免费完整版| 亚洲在线观看片| 嫩草影视91久久| 69av精品久久久久久| 一级毛片女人18水好多| 人人妻,人人澡人人爽秒播| 国内毛片毛片毛片毛片毛片| 看黄色毛片网站| 久久久成人免费电影| 国产精品一区二区三区四区久久| 国产av一区在线观看免费| 伊人久久大香线蕉亚洲五| 欧美日韩中文字幕国产精品一区二区三区| 国产高清videossex| 国产精品综合久久久久久久免费| 男插女下体视频免费在线播放| 国产亚洲精品av在线| 欧美av亚洲av综合av国产av| 亚洲欧美激情综合另类| 日本成人三级电影网站| 天天躁日日操中文字幕| 91九色精品人成在线观看| 巨乳人妻的诱惑在线观看| 国产精品美女特级片免费视频播放器 | av在线蜜桃| 91在线精品国自产拍蜜月 | 亚洲av成人精品一区久久| 黄色 视频免费看| 男插女下体视频免费在线播放| 亚洲av中文字字幕乱码综合| 国产av不卡久久| 亚洲精品456在线播放app | 国产欧美日韩精品亚洲av| 99久久成人亚洲精品观看| 手机成人av网站| 国产精品久久久av美女十八| 无限看片的www在线观看| 999精品在线视频| 亚洲精品色激情综合| 一边摸一边抽搐一进一小说| 狂野欧美激情性xxxx| 久久精品国产综合久久久| 久久人人精品亚洲av| 国产精品自产拍在线观看55亚洲| 午夜激情欧美在线| 一本精品99久久精品77| 久久久久久人人人人人| 亚洲熟妇熟女久久| 欧美成狂野欧美在线观看| 男女那种视频在线观看| 国产99白浆流出| 露出奶头的视频| 国内精品久久久久久久电影| 亚洲av中文字字幕乱码综合| 神马国产精品三级电影在线观看| 91av网站免费观看| e午夜精品久久久久久久| 搞女人的毛片| 一a级毛片在线观看| 国产精品美女特级片免费视频播放器 | 麻豆一二三区av精品| 久久精品影院6| 国产精品美女特级片免费视频播放器 | 久久欧美精品欧美久久欧美| x7x7x7水蜜桃| 黄色 视频免费看| 欧美色视频一区免费| 老熟妇乱子伦视频在线观看| 国产高清有码在线观看视频| 国产伦在线观看视频一区| 色噜噜av男人的天堂激情| 老司机午夜福利在线观看视频| 成年女人看的毛片在线观看| 国产精品日韩av在线免费观看| 国产av一区在线观看免费| 亚洲自拍偷在线| 国产单亲对白刺激| 1000部很黄的大片| 久久久久九九精品影院| 亚洲av五月六月丁香网| 香蕉丝袜av| 国产伦精品一区二区三区四那| 国产精品免费一区二区三区在线| 757午夜福利合集在线观看| 97超视频在线观看视频| 欧美一级a爱片免费观看看| www.自偷自拍.com| 好男人在线观看高清免费视频| 白带黄色成豆腐渣| 国产精品香港三级国产av潘金莲| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品av在线| 香蕉丝袜av| 国产1区2区3区精品| 99久久精品热视频| 国产精品亚洲av一区麻豆| 久久精品人妻少妇| 欧美乱妇无乱码| 超碰成人久久| 日韩av在线大香蕉| 欧美乱色亚洲激情| 老汉色av国产亚洲站长工具| 丁香六月欧美| 欧美黑人欧美精品刺激| 一本精品99久久精品77| 亚洲18禁久久av| 蜜桃久久精品国产亚洲av| 午夜福利在线在线| 精品午夜福利视频在线观看一区| 亚洲乱码一区二区免费版| 麻豆成人午夜福利视频| 国产av在哪里看| 国产真实乱freesex| 曰老女人黄片| 亚洲国产欧洲综合997久久,| 国内久久婷婷六月综合欲色啪| 国产av在哪里看| 国产成人影院久久av| 99久久精品一区二区三区| 国产一区二区三区在线臀色熟女| 脱女人内裤的视频| 国产高清视频在线观看网站| 免费一级毛片在线播放高清视频| 99精品欧美一区二区三区四区| 特大巨黑吊av在线直播| 久久精品人妻少妇| 成人鲁丝片一二三区免费| 欧美一区二区精品小视频在线| 这个男人来自地球电影免费观看| 日韩欧美 国产精品| 男插女下体视频免费在线播放| 757午夜福利合集在线观看| 国产不卡一卡二| 欧美在线一区亚洲| 9191精品国产免费久久| 国产三级黄色录像| 中文在线观看免费www的网站| 午夜a级毛片| 国产又黄又爽又无遮挡在线| 最新中文字幕久久久久 | 亚洲av片天天在线观看| 九色成人免费人妻av| 久久精品国产综合久久久| 长腿黑丝高跟| 亚洲天堂国产精品一区在线| 琪琪午夜伦伦电影理论片6080| 99热这里只有是精品50| 宅男免费午夜| www国产在线视频色| 久久精品aⅴ一区二区三区四区| 亚洲一区高清亚洲精品| 国产av在哪里看| 欧美极品一区二区三区四区| 最好的美女福利视频网| xxxwww97欧美| 小说图片视频综合网站| 夜夜爽天天搞| 又黄又爽又免费观看的视频| 淫秽高清视频在线观看| 在线十欧美十亚洲十日本专区| 美女高潮喷水抽搐中文字幕| 午夜成年电影在线免费观看| 少妇人妻一区二区三区视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲乱码一区二区免费版| 亚洲人与动物交配视频| 国内精品一区二区在线观看| 久久九九热精品免费| 国产精品久久久久久久电影 | 国产亚洲精品久久久com| 曰老女人黄片| 日韩三级视频一区二区三区| 国产激情欧美一区二区| 欧美黑人欧美精品刺激| 最新美女视频免费是黄的| 琪琪午夜伦伦电影理论片6080| 亚洲在线自拍视频| 偷拍熟女少妇极品色| 看黄色毛片网站| 亚洲欧美日韩卡通动漫| 黑人操中国人逼视频| 我要搜黄色片| 国产熟女xx| 欧美在线一区亚洲| www日本在线高清视频| 国模一区二区三区四区视频 | 中文亚洲av片在线观看爽| 午夜免费激情av| 亚洲熟女毛片儿| 久久这里只有精品19| 日韩精品青青久久久久久| 精品久久蜜臀av无| 午夜福利欧美成人| 国产精品一区二区免费欧美| 亚洲熟妇中文字幕五十中出| 久久精品亚洲精品国产色婷小说| 午夜亚洲福利在线播放| 午夜两性在线视频| 两性夫妻黄色片| 国产精品久久久久久亚洲av鲁大| bbb黄色大片| 国产高潮美女av| 三级毛片av免费| 国产伦一二天堂av在线观看| 久久久久九九精品影院| 国产一级毛片七仙女欲春2| 亚洲黑人精品在线| 国产精品国产高清国产av| 在线观看美女被高潮喷水网站 | 男人舔女人的私密视频| 操出白浆在线播放| 一a级毛片在线观看| 成人三级做爰电影| 午夜a级毛片| 99精品在免费线老司机午夜| www.999成人在线观看| 最新美女视频免费是黄的| 长腿黑丝高跟| 一区二区三区国产精品乱码| 成年女人永久免费观看视频| 国产激情偷乱视频一区二区| 免费搜索国产男女视频| 成年版毛片免费区| 狠狠狠狠99中文字幕| 99视频精品全部免费 在线 | 白带黄色成豆腐渣| 51午夜福利影视在线观看| 国语自产精品视频在线第100页| 夜夜爽天天搞| 亚洲午夜精品一区,二区,三区| 中国美女看黄片| 夜夜看夜夜爽夜夜摸| www.www免费av| 一级毛片精品| 欧美成人性av电影在线观看| 成年版毛片免费区| 日本 欧美在线| 日韩免费av在线播放| 欧美日韩一级在线毛片| av天堂在线播放| 国产精品永久免费网站| 99久久精品热视频| 国产高清videossex| 天天添夜夜摸| 特级一级黄色大片| 成在线人永久免费视频| 不卡av一区二区三区| 亚洲七黄色美女视频| 欧美日韩黄片免| 香蕉国产在线看| 波多野结衣高清无吗| 神马国产精品三级电影在线观看| 最新在线观看一区二区三区| 97人妻精品一区二区三区麻豆| 99热这里只有精品一区 | 亚洲av第一区精品v没综合| 亚洲精品在线美女| 99在线视频只有这里精品首页| 一区二区三区激情视频| 国产97色在线日韩免费| 搡老岳熟女国产| 免费看日本二区| 亚洲熟妇熟女久久| 国产成人精品久久二区二区免费| 免费看美女性在线毛片视频| 国产精品一区二区精品视频观看| 给我免费播放毛片高清在线观看| 黄色成人免费大全| 性色avwww在线观看| 午夜福利在线观看免费完整高清在 | 99视频精品全部免费 在线 | 神马国产精品三级电影在线观看| 国产97色在线日韩免费| 久久精品影院6| 一级毛片高清免费大全| 亚洲av电影在线进入| 黑人巨大精品欧美一区二区mp4| 禁无遮挡网站| 在线观看免费视频日本深夜| 老鸭窝网址在线观看| 伊人久久大香线蕉亚洲五| 又爽又黄无遮挡网站| 亚洲美女黄片视频| 精品国内亚洲2022精品成人| 午夜福利免费观看在线| 一区二区三区激情视频| 精品日产1卡2卡| 国产主播在线观看一区二区| 一级毛片女人18水好多| 日本 欧美在线| 黄色日韩在线| 在线看三级毛片| 国产精品一区二区三区四区久久| 日韩欧美在线乱码| 中国美女看黄片| 麻豆久久精品国产亚洲av| 成人亚洲精品av一区二区| 他把我摸到了高潮在线观看| 搞女人的毛片| 亚洲国产欧美一区二区综合| 伊人久久大香线蕉亚洲五| tocl精华| 国产av在哪里看| 在线播放国产精品三级| 18禁裸乳无遮挡免费网站照片| 国产高清激情床上av| 亚洲成av人片免费观看| 亚洲成av人片在线播放无| 美女午夜性视频免费| 好男人在线观看高清免费视频| 一个人看视频在线观看www免费 | 日韩大尺度精品在线看网址| 五月伊人婷婷丁香| 美女高潮的动态| 午夜福利高清视频|