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

    土地利用變化對盤古河流域徑流的影響

    2017-03-27 10:13:38郝振純
    水土保持通報 2017年1期
    關(guān)鍵詞:盤古徑流量水文

    趙 霞, 郝振純

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

    土地利用變化對盤古河流域徑流的影響

    趙 霞, 郝振純

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

    [目的] 分析大興安嶺地區(qū)盤古河流域土地利用變化情況,研究流域徑流對不同土地利用情景的響應(yīng),為合理規(guī)劃土地利用提供依據(jù)。[方法] 利用盤古河流域 1987,2000和2011年 3期土地利用數(shù)據(jù)和1988—2012年的水文氣象資料,結(jié)合SWAT(soil and water assessment tool)模型,分析盤古河流域土地利用變化的徑流響應(yīng)。[結(jié)果] 1987—2011年盤古河流域主要土地利用變化為裸地向林地、草地轉(zhuǎn)化,林地面積增加,裸地面積減少;汛期徑流量占年總徑流量的比例隨著林地的增加而減少;林地的增加使得多年平均年徑流量、汛期徑流量和最大月經(jīng)流量均呈減少趨勢。[結(jié)論] SWAT模型能夠較好地模擬盤古河流域的月流量過程;流域徑流隨著林地面積的增加而減少。

    土地利用變化; SWAT模型; 徑流變化; 盤古河流域

    文獻(xiàn)參數(shù): 趙霞, 郝振純.土地利用變化對盤古河流域徑流的影響[J].水土保持通報,2017,37(1):083-087.DOI:10.13961/j.cnki.stbctb.2017.01.015; Zhao Xia, Hao Zhenchun. Land use change effects on runoff in Pangu River basin[J]. Bulletin of Soil and Water Conservation, 2017,37(1):083-087.DOI:10.13961/j.cnki.stbctb.2017.01.015

    近年來變化環(huán)境下的水文水資源問題日益成為研究重點(diǎn)[1],其中土地利用/覆蓋變化作為環(huán)境變化的重要部分,對流域水文循環(huán)有著重大影響。土地利用變化的水文響應(yīng)研究也日益成為熱點(diǎn)問題[2]。

    流域水文模型是模擬流域水文過程、研究流域水文規(guī)律的重要理論基礎(chǔ)[3]。SWAT模型作為一種具有很強(qiáng)物理機(jī)制的流域尺度、連續(xù)時段、針對過程的綜合水文模型[4-5],目前已被廣泛應(yīng)用于流域土地利用/覆被變化的水文響應(yīng)研究。然而近年來,特定流域上土地利用變化的水文響應(yīng)結(jié)論卻很難取得一致。張圣微等[6]通過對拉薩河的研究表明隨著流域林草地面積的增加,年徑流量呈減少趨勢。郝振純等[7]設(shè)置了4種土地利用情景,發(fā)現(xiàn)海河流域典型區(qū)域的徑流量隨著草地面積的減少和林地面積的增加而減少。史曉亮等[8]分析諾敏河流域土地利用變化對徑流的影響后發(fā)現(xiàn),林草地面積的減少和耕地面積的增加使得多年平均徑流量增加了17.78%。但是部分學(xué)者的研究卻的出完全相反的結(jié)論。郝芳華等[9]通過對洛河上游流域的研究表明,森林的存在使得流域徑流量增加。王根緒等[10]在馬營河流域的研究發(fā)現(xiàn)隨著林草地向耕地的大規(guī)模轉(zhuǎn)化,流域年均徑流量減少了28.12%。因此,針對流域土地利用變化的水文響應(yīng),有必要進(jìn)行大量實(shí)例研究。

    大興安嶺是中國最大的原始林區(qū),林地面積7.30×106hm2,被譽(yù)為祖國的“綠色寶庫”,是東北華北地區(qū)環(huán)境安全的天然屏障。1987年“五·六”特大森林火災(zāi)燒毀林地1.01×106hm2,使得大興安嶺的森林覆蓋率由原來的76%下降到61.5%,下墊面土地利用/覆被發(fā)生巨大變化,一定程度上影響了該地區(qū)的生態(tài)環(huán)境安全。因此,大興安嶺地區(qū)土地利用變化的水文響應(yīng)研究十分必要。本研究擬選取大興安嶺地區(qū)盤古河流域?yàn)檠芯繀^(qū),結(jié)合3期實(shí)測遙感土地利用數(shù)據(jù)分析流域土地利用變化情況,構(gòu)建SWAT模型研究不同情景下的徑流響應(yīng),從而為該地區(qū)合理規(guī)劃土地利用提供依據(jù)。

    1 研究區(qū)概況

    盤古河為黑龍江上游右岸支流,發(fā)源于大興安嶺白卡魯山東麓,位于黑龍江省塔河縣西部,河流全長165 km,寬20~45 m,水深1.2 m,流域面積3 308km2,在23站東北10 km處注入黑龍江。盤古河流域地處中高緯度,大興安嶺東側(cè),屬于寒溫帶季風(fēng)氣候,雨熱同季,冬季干燥寒冷且漫長,夏季多雨高溫。流域多年平均氣溫-2.0 ℃左右,無霜期在80 d左右。流域降雨年內(nèi)分配較不均勻,降雨多集中于6—9月,多年平均降雨量約420 mm。

    2 數(shù)據(jù)準(zhǔn)備與模型構(gòu)建

    2.1 基礎(chǔ)數(shù)據(jù)

    模型所需地形數(shù)據(jù)來源于國際科學(xué)數(shù)據(jù)服務(wù)平臺提供的90 m分辨率的SRTM(shuttle radar topography mission)數(shù)據(jù),利用ArcGIS軟件對其進(jìn)行拼接、投影、剪切等一系列預(yù)處理操作,生成模型所需的DEM。土壤數(shù)據(jù)來自于中科院地理所數(shù)據(jù)中心的1∶100萬的土壤矢量數(shù)據(jù)集。土壤屬性數(shù)據(jù)參數(shù)由中國土壤數(shù)據(jù)庫查詢得到,以及采用SPAW軟件[11]計(jì)算獲得。水文資料和雨量資料來源于水文年鑒,包括流域內(nèi)盤古水文站1988—2012年逐日流量資料以及盤古站、23站逐日降雨資料。氣象資料采用流域臨近塔河氣象站1988—2012年逐日流量資料以及盤古站、23站逐日降雨資料。氣象資料采用流域臨近塔河氣象站1988—2012年的氣象數(shù)據(jù),數(shù)據(jù)來源于中國氣象科學(xué)數(shù)據(jù)共享網(wǎng),包括逐日相對濕度、日照時數(shù)、平均風(fēng)速、最低氣溫和最高氣溫。日太陽輻射量由模型內(nèi)置的天氣發(fā)生器模擬得到。

    研究區(qū)的土地利用數(shù)據(jù)選取1987,2000和2011年3期Landsat TM影像圖,利用ENVI軟件監(jiān)督分類功能進(jìn)行人機(jī)交互目視解譯,將土地利用類型劃分為6類:林地(FRST)、草地(PAST)、耕地(AGRL)、居民區(qū)(低密度)(URLD)、水體(WATR)、裸地(BARR),流域土地利用分布見附圖1。

    2.2 模型構(gòu)建

    SWAT模型在模擬過程中,基于DEM柵格將研究流域空間上離散化成多個子流域單元,各子流域能夠體現(xiàn)出氣溫、降水、蒸發(fā)等氣候因素異性,土壤、土地利用、海拔、坡度等下墊面因子的空間變化以及管理方式的不同。模型基于最小集水面積劃分子流域,本研究最小集水面積取60 km2,將盤古河流域劃分為 31個子流域。HRU的生成采用multiple hydrologic response units方法,將土地利用類型最小面積比重設(shè)為5%,土壤類型最小面積比重設(shè)為10%,坡度最小面積比重設(shè)為10%。潛在蒸散發(fā)選擇彭曼—蒙蒂斯?jié)撛谡羯l(fā)公式計(jì)算,地表徑流模擬采用 SCS 徑流曲線數(shù)法,河道匯流演算采用Variable Storage法。

    3 研究結(jié)果

    3.1 模型驗(yàn)證及評價

    以2000年土地利用數(shù)據(jù)為輸入,利用1988—2012年盤古站月徑流數(shù)據(jù)對模型進(jìn)行校準(zhǔn)及驗(yàn)證,其中1988—2002年為校準(zhǔn)期,2003—2012年為驗(yàn)證期。由于SWAT模型輸入?yún)?shù)較多,與徑流有關(guān)的25個參數(shù)對模型結(jié)果影響程度不同,本文通過敏感性分析選取對模擬精度影響較為顯著的10個參數(shù)進(jìn)行參數(shù)率定。表1即為敏感性前10的參數(shù),以及各參數(shù)物理意義、調(diào)參方式和最優(yōu)參數(shù)值。

    本文選用決定系數(shù)R2,Nash-Sutcliffe系數(shù)(ENS)[12]來檢驗(yàn)?zāi)P偷哪M效果,具體計(jì)算公式如下:

    (1)

    模型校準(zhǔn)及驗(yàn)證結(jié)果見圖1,評價指標(biāo)見表2。

    表1 參數(shù)敏感性分析結(jié)果及參數(shù)最優(yōu)值

    表2 盤古水文站月徑流模擬結(jié)果評價

    納什系數(shù)ENS表示模擬值與實(shí)測值之間的擬合度,決定系數(shù)R2表示實(shí)測值與模擬值之間的吻合程度,通常認(rèn)為ENS越大,R2越高,模型模擬結(jié)果越好,通常認(rèn)為R2高于0.65,ENS高于0.60時,模型模擬精度令人滿意[13-14]。由表2可以看出,無論是校準(zhǔn)期還是驗(yàn)證期,ENS和R2都大于0.70,說明SWAT模型適用于盤古河流域盤古站的徑流模擬。

    圖1 校準(zhǔn)期和驗(yàn)證期月徑流模擬值與實(shí)測值對比

    3.2 1987—2011年間盤古河流域土地利用變化

    根據(jù)對解譯出的盤古河流域1987,2000和2011年3期土地利用數(shù)據(jù)的統(tǒng)計(jì)對比分析,得到各土地利用類型的面積、百分比和動態(tài)度(表3)。

    表3 1987—2011年盤古河流域土地利用類型面積變化及動態(tài)度

    長期以來,林地和草地一直是盤古河流域主導(dǎo)景觀,然而由于1987年5月6日至6月2日大興安嶺地區(qū)的特大森林火災(zāi),1987年土地利用數(shù)據(jù)中林地面積僅為1 132.85 km2,僅占流域面積的35.33%,裸地面積達(dá)1 305.48 km2,占流域總面積的40.71%。隨著多年的恢復(fù),2000年土地利用數(shù)據(jù)中林地面積為1 745.96 km2,草地面積為1 275.54 km2,相較于1987年面積分別增加了613.11和534.61 km2,年變化面積分別為47.16和41.12 km2/a,動態(tài)度分別為4.16%和5.55%。2011年土地利用數(shù)據(jù)中林地面積為2 414.52 km2,草地面積為689.41 km2,相較于2000年林地面積增加了668.56 km2,草地面積減少了586.09 km2,年變化面積分別為60.78和-53.29 km2/a,動態(tài)度分別為3.48%和-4.18%。1987—2011年期間,與林地面積增長相對應(yīng),裸地面積不斷減少,由1987年的1305.48 km2減少到2011年的90.46 km2,同時耕地面積不斷增長,居民區(qū)和水體面積總體呈下降趨勢。

    由于盤古河流域出口23站自1986年撤銷流量觀測,本研究采用徑流數(shù)據(jù)為流域內(nèi)盤古站觀測數(shù)據(jù),因而需進(jìn)一步分析盤古站集水面積上流域土地利用變化情況。盤古站集水面積上各土地利用類型面積及百分比見表4。

    表4 1987—2011年盤古站集水面積

    由表4可以看出,1987—2011年期間,林地面積由856.90 km2增長到1 056.83 km2,面積增加了近總集水面積的20%;裸地面積由83.65 km2下降為28.75 km2,面積由占集水面積的6.61%下降為2.27%;草地面積則經(jīng)歷了先增加后減少的過程,由1987年的318.64 km2增加到2000年的401.13 km2,再減少到2011年的176.35 km2;水體及居民區(qū)面積變化較小。

    3.3 盤古河流域土地利用變化的水文響應(yīng)分析

    利用盤古河流域1987,2000和2011年3期土地利用數(shù)據(jù),采用相同的氣象、土壤數(shù)據(jù),僅改變土地利用輸入,分別模擬盤古站1988—2012年的徑流過程,得到3期土地利用條件下月徑流過程(圖2),通過對比分析得到流域土地利用變化對水文過程的影響。

    圖2 3期土地利用情景下月徑流模擬值對比

    模擬得到的不同土地利用情景下的多年平均徑流深、多年平均汛期徑流深(6—9月)和多年平均最大月徑流深(7月)結(jié)果及變化情況見表5。在1987年土地利用情景下,盤古站多年平均徑流深為217.81 mm,其中多年平均汛期徑流深為155.84 mm,多年平均最大月徑流深為48.49 mm;在2000年土地利用情景下,盤古站多年平均徑流深為216.87 mm,其中多年平均汛期徑流深為153.31 mm,多年平均最大月徑流深為47.50 mm;在2011年土地利用情景下,盤古站多年平均徑流深為215.36 mm,其中多年平均汛期徑流深為151.37 mm,多年平均最大月徑流深為46.07 mm。

    表5 3期土地利用情景下模擬的徑流深

    3種情景下汛期(6—9月)流量占全年總徑流量的比例都在70%左右,分別為71.55%,70.69%和70.29%,表明隨著林草地面積的增加,占比呈下降趨勢。1987年到2000年,主要土地利用類型變化為由裸地向草地林地轉(zhuǎn)化,多年平均徑流量下降了0.43%,汛期徑流深下降了1.62%,最大月徑流量下降了2.04%;2000—2011年,主要土地利用類型變化為由草地向林地轉(zhuǎn)化,多年平均徑流量下降了0.70%,汛期徑流量下降了1.27%,最大月徑流量下降了3.01%。1987—2011年,盤古站集水面積上林地面積增加了20%,多年平均徑流量下降了1.12%,汛期徑流量下降了2.87%,最大月徑流量下降了4.99%。結(jié)果表明,隨著流域林地面積的增加,流域徑流呈下降趨勢;年徑流量、汛期徑流量和最大月徑流量對土地利用變化的響應(yīng)程度呈增加趨勢。

    4 結(jié) 論

    為研究土地利用變化對流域徑流的影響,利用1987,2000和2011年3期土地利用數(shù)據(jù)和流域1988—2012年水文氣象資料,建立SWAT模型對徑流進(jìn)行模擬。模型在校準(zhǔn)期納什系數(shù)為0.71,決定系數(shù)為0.71;驗(yàn)證期納什系數(shù)為0.73,決定系數(shù)為0.76,表明SWAT模型適用于盤古河流域的徑流模擬。

    1987—2000年期間,盤古河流域主要土地利用變化為裸地面積減少,林地、草地面積增加,年變化量分別為-88.07,47.16和41.12 km2/a。2000—2011年期間,盤古河流域主要土地利用變化為草地面積減少,林地面積增加,年變化量分別為-53.29,60.78 km2/a。多年主要土地利用變化為裸地向草地、林地轉(zhuǎn)化。

    針對1987,2000和2011年3期土地利用數(shù)據(jù)進(jìn)行的徑流模擬結(jié)果說明:當(dāng)流域林地面積增加時,盤古站年徑流量、汛期徑流量及最大月徑流量均呈下降趨勢且響應(yīng)程度增加,同時汛期流量占全年總徑流量的比例呈下降趨勢。

    [1] 王茂運(yùn),謝朝勇.變化環(huán)境下的水資源安全問題探討[J].中國水運(yùn),2013(4):210-211.

    [2] 李麗娟,姜德娟,楊俊偉,等.陜西大理河流域土地利用/覆被變化的水文效應(yīng)[J].地理研究,2010,29(7):1233-1243.

    [3] 王中根,鄭紅星,劉昌明,等.黃河典型流域分布式水文模型及應(yīng)用研究[J].中國科學(xué)(E輯):技術(shù)科學(xué),2004,34(S1):49-59.

    [4] Srinivasan R, Arnold J G, Jones C A. Hydrologic modelling of the United States with the soil and water assessment tool[J]. International Journal of Water Resources Development, 1998,14(3):315-325.

    [5] Arnold J G, Allen P M. Estimating hydrologic budgets for three Illinois watersheds[J]. Journal of Hydrology, 1996,176(1):57-77.

    [6] 張圣微,雷玉平,姚琴,等.土地覆被和氣候變化對拉薩河流域徑流量的影響[J].水資源保護(hù),2010,26(2):39-44.

    [7] 郝振純,蘇振寬.土地利用變化對海河流域典型區(qū)域的徑流影響[J].水科學(xué)進(jìn)展,2015,26(4):491-499.

    [8] 史曉亮,李穎,趙凱,等.諾敏河流域土地利用與覆被變化及其對水文過程的影響[J].水土保持通報,2013,33(1):23-28.

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

    [10] 王根緒,張鈺,劉桂民,等.馬營河流域1967—2000年土地利用變化對河流徑流的影響[J]. 中國科學(xué)(D輯):地球科學(xué),2005,35(7):671-681.

    [11] 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.

    [12] Moriasi D N, Arnold J G, Van Liew M W, et al. Model evaluation guidelines for system atic-quantification of accuracy in watershed simulations[J]. Transactions of the ASABE, 2007,50(3):885-900.

    [13] 龐靖鵬,劉昌明,徐宗學(xué).密云水庫流域土地利用變化對產(chǎn)流和產(chǎn)沙的影響[J].北京師范大學(xué)學(xué)報:自然科學(xué)版,2010,46(3):290-299.

    [14] 張利平,曾思棟,王任超,等.氣候變化對灤河流域水文循環(huán)的影響及模擬[J].資源科學(xué),2011,33(5):966-974.

    Land Use Change Effects on Runoff in Pangu River Basin

    ZHAO Xia, HAO Zhenchun

    (StateKeyLaboratoryofHydrology&WaterResourcesandHydraulicEngineering,HohaiUniversity,Nanjing,Jiangsu210098,China)

    [Objective] The objective of the study is to analyze the land use change in Pangu River basin in Daxinganling region and explore the responses of runoff to different land use scenarios and provide the basis for the rational land use planning. [Methods] Using the land use data of 1987, 2000, 2011 and the hydro meteorology data from 1988 to 2012, combining the soil and water assessment tool model, we analyzed response of runoff to land use changes. [Results] The land use changes during period from 1987 to 2011 included a transformation from barren land to forest and grassland, an increase in forest land and a decrease in barren land. The percentage of runoff in flooding season decreases with increasing forest area. The increase in forest area reduced the long-time average annual runoff, the runoff in flooding season and the maximum monthly peak discharge. [Conclusion] The SWAT model could well simulate monthly flow process of the Pangu River basin and the runoff decreased with increasing forest area.

    land use change; SWAT model; runoff changes; Pangu River basin

    2016-06-21

    2016-07-03

    國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目“多尺度水文水資源預(yù)報預(yù)測預(yù)警關(guān)鍵技術(shù)及應(yīng)用研究”(2016YFC0402704); 水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室專項(xiàng)經(jīng)費(fèi)資助項(xiàng)目(1069-514031112)

    趙霞(1992—),女(漢族),江蘇省洪澤縣人,碩士研究生,研究方向?yàn)樗奈锢硪?guī)律模擬及水文預(yù)報。E-mail:541978989@qq.com。

    A

    1000-288X(2017)01-0083-05

    TV121+.2

    猜你喜歡
    盤古徑流量水文
    2022年《中國水文年報》發(fā)布
    水文
    水文水資源管理
    盤古開天辟地
    盤古
    盤古
    盤古開天辟地
    水文比擬法在計(jì)算河川徑流量時的修正
    水文
    SCS模型在紅壤土坡地降雨徑流量估算中的應(yīng)用
    精品视频人人做人人爽| 日韩电影二区| 99九九在线精品视频| 99热国产这里只有精品6| 一区二区av电影网| 一区二区三区精品91| 国产精品偷伦视频观看了| 香蕉丝袜av| 91成年电影在线观看| 亚洲成av片中文字幕在线观看| 19禁男女啪啪无遮挡网站| 欧美激情 高清一区二区三区| 国产深夜福利视频在线观看| 韩国高清视频一区二区三区| 欧美乱码精品一区二区三区| 欧美人与性动交α欧美软件| 超碰成人久久| 久久久久久久精品精品| 午夜91福利影院| 无限看片的www在线观看| svipshipincom国产片| 久久久久视频综合| 国产亚洲av高清不卡| 夜夜骑夜夜射夜夜干| 亚洲va日本ⅴa欧美va伊人久久 | 9191精品国产免费久久| 精品少妇一区二区三区视频日本电影| 侵犯人妻中文字幕一二三四区| 亚洲精品成人av观看孕妇| 天天躁日日躁夜夜躁夜夜| www日本在线高清视频| 亚洲自偷自拍图片 自拍| 国产精品熟女久久久久浪| 色老头精品视频在线观看| 亚洲欧美精品综合一区二区三区| 老汉色∧v一级毛片| 极品少妇高潮喷水抽搐| 在线观看一区二区三区激情| 久久久久久免费高清国产稀缺| 欧美性长视频在线观看| 国产激情久久老熟女| 如日韩欧美国产精品一区二区三区| 女性生殖器流出的白浆| 国产伦理片在线播放av一区| 日韩精品免费视频一区二区三区| 亚洲成人免费电影在线观看| 丰满饥渴人妻一区二区三| 亚洲专区国产一区二区| 国产成人欧美| 丝瓜视频免费看黄片| 国产亚洲一区二区精品| 亚洲精品国产av成人精品| 久久国产亚洲av麻豆专区| 99久久精品国产亚洲精品| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利在线观看吧| 香蕉丝袜av| 久久亚洲精品不卡| 亚洲国产毛片av蜜桃av| 亚洲精品国产区一区二| 九色亚洲精品在线播放| 亚洲avbb在线观看| 日韩制服骚丝袜av| 男女边摸边吃奶| 亚洲人成77777在线视频| 亚洲,欧美精品.| 久热这里只有精品99| 在线观看免费视频网站a站| 国产淫语在线视频| 国产在线视频一区二区| 叶爱在线成人免费视频播放| 满18在线观看网站| 美女高潮到喷水免费观看| 好男人电影高清在线观看| 精品久久久精品久久久| 我要看黄色一级片免费的| 我要看黄色一级片免费的| 不卡av一区二区三区| 色婷婷久久久亚洲欧美| 久久久久网色| 婷婷成人精品国产| 亚洲国产看品久久| 一区福利在线观看| 无遮挡黄片免费观看| 国产主播在线观看一区二区| 国产在线观看jvid| 巨乳人妻的诱惑在线观看| 在线观看一区二区三区激情| 美女午夜性视频免费| 久热这里只有精品99| 成年人免费黄色播放视频| 国产精品自产拍在线观看55亚洲 | 亚洲av片天天在线观看| av在线老鸭窝| 老鸭窝网址在线观看| 亚洲第一欧美日韩一区二区三区 | 亚洲一区中文字幕在线| 人妻 亚洲 视频| 永久免费av网站大全| 色综合欧美亚洲国产小说| 欧美人与性动交α欧美软件| 国产成人系列免费观看| 欧美精品一区二区大全| 国产片内射在线| 狂野欧美激情性bbbbbb| 免费在线观看完整版高清| 大香蕉久久成人网| 国产成人免费无遮挡视频| 亚洲专区中文字幕在线| 大陆偷拍与自拍| 少妇人妻久久综合中文| 欧美精品啪啪一区二区三区 | 桃红色精品国产亚洲av| 国产av国产精品国产| 亚洲国产欧美日韩在线播放| 90打野战视频偷拍视频| 国产激情久久老熟女| 下体分泌物呈黄色| 美女高潮喷水抽搐中文字幕| 久久精品久久久久久噜噜老黄| 美女福利国产在线| 丁香六月天网| 午夜福利在线观看吧| 久久久久精品国产欧美久久久 | 青草久久国产| 国产高清国产精品国产三级| 12—13女人毛片做爰片一| 美女高潮喷水抽搐中文字幕| 极品人妻少妇av视频| av电影中文网址| 国产精品麻豆人妻色哟哟久久| 国产av精品麻豆| 极品人妻少妇av视频| 高清欧美精品videossex| 亚洲欧美色中文字幕在线| 国产伦人伦偷精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲自偷自拍图片 自拍| 成人国产一区最新在线观看| 精品久久久久久久毛片微露脸 | 秋霞在线观看毛片| 国产日韩欧美亚洲二区| 最近最新中文字幕大全免费视频| 国产一区二区 视频在线| 飞空精品影院首页| 狠狠婷婷综合久久久久久88av| 极品人妻少妇av视频| 咕卡用的链子| 久久毛片免费看一区二区三区| 男女免费视频国产| 精品少妇久久久久久888优播| 三上悠亚av全集在线观看| 精品少妇黑人巨大在线播放| 少妇人妻久久综合中文| 亚洲国产精品成人久久小说| 亚洲男人天堂网一区| 视频区图区小说| 夜夜骑夜夜射夜夜干| 精品一品国产午夜福利视频| 超色免费av| 一级毛片女人18水好多| 亚洲欧美色中文字幕在线| 最近最新免费中文字幕在线| 真人做人爱边吃奶动态| 巨乳人妻的诱惑在线观看| 欧美在线一区亚洲| 久久青草综合色| 日本五十路高清| 人成视频在线观看免费观看| 亚洲国产精品一区三区| 一个人免费看片子| 青草久久国产| 欧美日韩国产mv在线观看视频| 亚洲国产看品久久| 日韩有码中文字幕| 水蜜桃什么品种好| 国产免费现黄频在线看| 精品欧美一区二区三区在线| 久久久精品区二区三区| 香蕉丝袜av| www.熟女人妻精品国产| av线在线观看网站| 男女下面插进去视频免费观看| 国产高清国产精品国产三级| 精品亚洲成国产av| 国产黄频视频在线观看| 又黄又粗又硬又大视频| 欧美激情 高清一区二区三区| 日本一区二区免费在线视频| 亚洲专区国产一区二区| 精品国产一区二区三区久久久樱花| av片东京热男人的天堂| 999久久久精品免费观看国产| 精品久久久久久久毛片微露脸 | 精品人妻一区二区三区麻豆| 久久性视频一级片| 一边摸一边抽搐一进一出视频| 日韩电影二区| 人人妻人人澡人人看| 女性被躁到高潮视频| 99久久人妻综合| 97在线人人人人妻| 欧美精品一区二区大全| 嫁个100分男人电影在线观看| 久久久精品国产亚洲av高清涩受| 欧美乱码精品一区二区三区| 欧美在线黄色| 满18在线观看网站| 久久久久视频综合| 热re99久久国产66热| 啦啦啦啦在线视频资源| 欧美黑人精品巨大| 免费看十八禁软件| 精品一区二区三区av网在线观看 | 这个男人来自地球电影免费观看| 午夜激情久久久久久久| 国产成人a∨麻豆精品| 一二三四在线观看免费中文在| 老熟妇仑乱视频hdxx| 国产精品国产av在线观看| 国产欧美亚洲国产| 欧美精品一区二区免费开放| 日本vs欧美在线观看视频| 黄色怎么调成土黄色| netflix在线观看网站| 久久久久久久大尺度免费视频| 久久av网站| 国产av精品麻豆| 超色免费av| 国产高清videossex| 50天的宝宝边吃奶边哭怎么回事| 夜夜骑夜夜射夜夜干| 在线亚洲精品国产二区图片欧美| 91成人精品电影| 亚洲成人免费av在线播放| 久久精品国产a三级三级三级| 老熟妇仑乱视频hdxx| 精品一区二区三卡| 69av精品久久久久久 | av一本久久久久| 亚洲国产成人一精品久久久| 欧美日韩亚洲综合一区二区三区_| av一本久久久久| 伊人久久大香线蕉亚洲五| 男女之事视频高清在线观看| 黑人猛操日本美女一级片| 99热全是精品| av有码第一页| av一本久久久久| 久久久久久久国产电影| 国产精品 国内视频| 91精品国产国语对白视频| 在线观看一区二区三区激情| bbb黄色大片| 99久久99久久久精品蜜桃| 80岁老熟妇乱子伦牲交| 日韩一卡2卡3卡4卡2021年| 在线亚洲精品国产二区图片欧美| 国产亚洲欧美在线一区二区| 亚洲第一欧美日韩一区二区三区 | 欧美+亚洲+日韩+国产| 久久人妻熟女aⅴ| 国产日韩一区二区三区精品不卡| 搡老岳熟女国产| 久久精品亚洲av国产电影网| 成年人午夜在线观看视频| 一二三四社区在线视频社区8| 日本撒尿小便嘘嘘汇集6| 精品少妇一区二区三区视频日本电影| 少妇裸体淫交视频免费看高清 | 免费观看人在逋| 成人影院久久| 日韩 欧美 亚洲 中文字幕| 久9热在线精品视频| 高清欧美精品videossex| 午夜福利,免费看| 欧美一级毛片孕妇| 999精品在线视频| 亚洲国产欧美一区二区综合| 咕卡用的链子| 国产精品九九99| 男女国产视频网站| 亚洲第一av免费看| 18禁黄网站禁片午夜丰满| a 毛片基地| 一个人免费在线观看的高清视频 | 涩涩av久久男人的天堂| 一进一出抽搐动态| 少妇粗大呻吟视频| 国产精品偷伦视频观看了| 国产老妇伦熟女老妇高清| 亚洲av日韩精品久久久久久密| 亚洲欧洲精品一区二区精品久久久| av在线app专区| 久久九九热精品免费| 免费日韩欧美在线观看| 欧美日韩亚洲综合一区二区三区_| 黄色怎么调成土黄色| 国产主播在线观看一区二区| 欧美一级毛片孕妇| 一个人免费看片子| 99久久精品国产亚洲精品| 超碰97精品在线观看| 欧美精品人与动牲交sv欧美| 在线十欧美十亚洲十日本专区| av有码第一页| 亚洲色图综合在线观看| 国产伦理片在线播放av一区| 国产高清videossex| 国产一区有黄有色的免费视频| 女警被强在线播放| 成人av一区二区三区在线看 | 亚洲 国产 在线| 黄片大片在线免费观看| 秋霞在线观看毛片| 中国美女看黄片| 免费在线观看视频国产中文字幕亚洲 | 欧美黑人欧美精品刺激| 男女高潮啪啪啪动态图| 午夜福利一区二区在线看| 狂野欧美激情性xxxx| 嫩草影视91久久| 一区二区三区精品91| 国产精品免费视频内射| 亚洲av日韩在线播放| 99久久99久久久精品蜜桃| 久久久久精品人妻al黑| 高潮久久久久久久久久久不卡| 他把我摸到了高潮在线观看 | 欧美国产精品va在线观看不卡| 岛国毛片在线播放| 日日夜夜操网爽| 国产伦理片在线播放av一区| 国产精品一区二区在线不卡| 黄片小视频在线播放| 熟女少妇亚洲综合色aaa.| 精品人妻1区二区| 在线观看免费午夜福利视频| 激情视频va一区二区三区| 69av精品久久久久久 | av在线老鸭窝| 首页视频小说图片口味搜索| 国产精品久久久av美女十八| 首页视频小说图片口味搜索| 无遮挡黄片免费观看| 色婷婷av一区二区三区视频| 乱人伦中国视频| 最黄视频免费看| 国产精品久久久久久精品古装| 十八禁人妻一区二区| 国产精品99久久99久久久不卡| 欧美xxⅹ黑人| 热99国产精品久久久久久7| 日韩视频一区二区在线观看| 中文字幕高清在线视频| 亚洲中文av在线| 久久99热这里只频精品6学生| a级片在线免费高清观看视频| 亚洲精品乱久久久久久| 少妇精品久久久久久久| 亚洲国产av新网站| 女人高潮潮喷娇喘18禁视频| e午夜精品久久久久久久| 飞空精品影院首页| 最黄视频免费看| 黄色 视频免费看| 免费高清在线观看视频在线观看| 建设人人有责人人尽责人人享有的| 电影成人av| 在线 av 中文字幕| cao死你这个sao货| 一边摸一边抽搐一进一出视频| 天堂中文最新版在线下载| 中文字幕高清在线视频| videosex国产| 欧美国产精品一级二级三级| 动漫黄色视频在线观看| 亚洲av欧美aⅴ国产| 亚洲专区字幕在线| 美女脱内裤让男人舔精品视频| netflix在线观看网站| 亚洲性夜色夜夜综合| 欧美乱码精品一区二区三区| 国产av精品麻豆| 在线观看免费高清a一片| 国产精品九九99| 日日摸夜夜添夜夜添小说| 99国产精品免费福利视频| 男女高潮啪啪啪动态图| 久久毛片免费看一区二区三区| 欧美少妇被猛烈插入视频| 搡老岳熟女国产| 国产精品亚洲av一区麻豆| 丁香六月天网| 国产成人欧美在线观看 | 下体分泌物呈黄色| 国产精品久久久久久人妻精品电影 | 一本综合久久免费| 免费观看人在逋| 天天躁日日躁夜夜躁夜夜| 成人亚洲精品一区在线观看| 91精品国产国语对白视频| av有码第一页| 丝袜美足系列| 国产一区二区三区在线臀色熟女 | 狠狠婷婷综合久久久久久88av| 亚洲欧美一区二区三区久久| 俄罗斯特黄特色一大片| 中文字幕精品免费在线观看视频| 国产伦人伦偷精品视频| 亚洲精品国产av成人精品| 久久国产精品大桥未久av| 国产1区2区3区精品| 亚洲精品国产精品久久久不卡| xxxhd国产人妻xxx| videosex国产| 热re99久久国产66热| 久久精品人人爽人人爽视色| 欧美中文综合在线视频| 欧美97在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品麻豆人妻色哟哟久久| 最新在线观看一区二区三区| 精品亚洲成a人片在线观看| 久久这里只有精品19| 欧美日韩亚洲高清精品| av福利片在线| av在线app专区| 狂野欧美激情性xxxx| 午夜福利视频在线观看免费| 精品欧美一区二区三区在线| 两个人免费观看高清视频| 欧美少妇被猛烈插入视频| 两个人免费观看高清视频| 丰满人妻熟妇乱又伦精品不卡| 80岁老熟妇乱子伦牲交| av天堂久久9| 在线永久观看黄色视频| 黄色a级毛片大全视频| 免费在线观看完整版高清| 精品国产乱子伦一区二区三区 | 中文字幕精品免费在线观看视频| 免费观看av网站的网址| 99国产精品一区二区三区| 人成视频在线观看免费观看| 久久久久国产精品人妻一区二区| 久久精品成人免费网站| 热re99久久精品国产66热6| 一区在线观看完整版| 韩国高清视频一区二区三区| 黑人巨大精品欧美一区二区mp4| 日韩制服丝袜自拍偷拍| 久久精品国产a三级三级三级| 国产精品自产拍在线观看55亚洲 | 亚洲国产精品999| 久久国产精品人妻蜜桃| 亚洲精品美女久久av网站| 亚洲精品久久午夜乱码| 久久久水蜜桃国产精品网| 啦啦啦视频在线资源免费观看| 亚洲av日韩精品久久久久久密| 国产亚洲精品一区二区www | 热99久久久久精品小说推荐| 黑人欧美特级aaaaaa片| 精品人妻熟女毛片av久久网站| 国产av国产精品国产| 久久久国产精品麻豆| 又黄又粗又硬又大视频| 亚洲av欧美aⅴ国产| 亚洲三区欧美一区| 99热国产这里只有精品6| 国产日韩欧美亚洲二区| 久久久国产一区二区| 另类亚洲欧美激情| 啦啦啦在线免费观看视频4| 视频区欧美日本亚洲| 亚洲精品国产av成人精品| 亚洲av欧美aⅴ国产| 韩国高清视频一区二区三区| 热99re8久久精品国产| 国产成人系列免费观看| 欧美精品av麻豆av| 国产亚洲欧美在线一区二区| 欧美 亚洲 国产 日韩一| 精品福利永久在线观看| 亚洲精品第二区| 国产一区二区三区av在线| 国产深夜福利视频在线观看| 久久久久国产一级毛片高清牌| 精品一品国产午夜福利视频| 青青草视频在线视频观看| 成人亚洲精品一区在线观看| 久久久久网色| 欧美成狂野欧美在线观看| 国产在线一区二区三区精| 大片电影免费在线观看免费| 久久中文看片网| 精品高清国产在线一区| 男女之事视频高清在线观看| 国产亚洲av高清不卡| 久久精品久久久久久噜噜老黄| av在线老鸭窝| 中国国产av一级| 亚洲av日韩精品久久久久久密| 久久精品亚洲熟妇少妇任你| 嫁个100分男人电影在线观看| 狠狠精品人妻久久久久久综合| 亚洲免费av在线视频| 久久精品国产亚洲av香蕉五月 | 欧美另类亚洲清纯唯美| 国产精品二区激情视频| 国产亚洲精品一区二区www | 狠狠婷婷综合久久久久久88av| 亚洲欧美精品综合一区二区三区| 精品熟女少妇八av免费久了| 一个人免费看片子| 久久影院123| 国产一区二区激情短视频 | 18在线观看网站| 91字幕亚洲| 午夜激情久久久久久久| www.999成人在线观看| 91老司机精品| 国产精品一区二区在线观看99| av免费在线观看网站| bbb黄色大片| 午夜免费观看性视频| 黄色毛片三级朝国网站| 大型av网站在线播放| 女人被躁到高潮嗷嗷叫费观| av网站在线播放免费| 亚洲国产日韩一区二区| 欧美人与性动交α欧美软件| 香蕉国产在线看| 午夜影院在线不卡| 久久青草综合色| 国产福利在线免费观看视频| 亚洲国产成人一精品久久久| 人妻久久中文字幕网| 成年人免费黄色播放视频| 成年人午夜在线观看视频| 99精品久久久久人妻精品| 精品视频人人做人人爽| 午夜福利视频精品| 亚洲精品自拍成人| 男女午夜视频在线观看| 国产精品一区二区免费欧美 | 亚洲中文字幕日韩| 国产成人免费观看mmmm| 18禁国产床啪视频网站| 天天躁日日躁夜夜躁夜夜| 亚洲五月婷婷丁香| 午夜影院在线不卡| 亚洲色图综合在线观看| 亚洲国产欧美一区二区综合| 国产精品影院久久| 亚洲av电影在线观看一区二区三区| 亚洲九九香蕉| 亚洲三区欧美一区| 精品久久蜜臀av无| 精品少妇黑人巨大在线播放| 色老头精品视频在线观看| 99热国产这里只有精品6| 久久精品国产亚洲av香蕉五月 | 久久人人爽av亚洲精品天堂| 国产伦人伦偷精品视频| 成人av一区二区三区在线看 | 在线天堂中文资源库| 亚洲国产精品成人久久小说| 热99re8久久精品国产| 国产欧美日韩一区二区三 | 日本黄色日本黄色录像| 欧美日韩av久久| 菩萨蛮人人尽说江南好唐韦庄| 另类亚洲欧美激情| 亚洲欧美清纯卡通| 午夜精品久久久久久毛片777| 国产主播在线观看一区二区| netflix在线观看网站| 这个男人来自地球电影免费观看| 无遮挡黄片免费观看| 法律面前人人平等表现在哪些方面 | 午夜激情av网站| 成年人免费黄色播放视频| 1024视频免费在线观看| 97在线人人人人妻| 可以免费在线观看a视频的电影网站| 欧美精品高潮呻吟av久久| 中文字幕人妻熟女乱码| av线在线观看网站| 国产精品一区二区免费欧美 | 色精品久久人妻99蜜桃| 下体分泌物呈黄色| 一区二区av电影网| 嫩草影视91久久| 亚洲欧洲精品一区二区精品久久久| 欧美激情高清一区二区三区| 男人爽女人下面视频在线观看| 午夜老司机福利片| 天天躁狠狠躁夜夜躁狠狠躁| xxxhd国产人妻xxx| 国产一区二区激情短视频 | 国产精品久久久人人做人人爽| 黑人操中国人逼视频| 亚洲精品美女久久av网站| 少妇 在线观看| 亚洲色图 男人天堂 中文字幕| 丁香六月天网| 亚洲精品av麻豆狂野| 99精国产麻豆久久婷婷| 老汉色∧v一级毛片| 国产欧美日韩一区二区精品| 久久久精品免费免费高清| 久久久久视频综合| videosex国产|