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

    張峰水庫(kù)入庫(kù)流量模擬及產(chǎn)流分析

    2019-01-26 10:25:04侯芳李曉暉胡文慧李光永
    科技創(chuàng)新與應(yīng)用 2019年1期
    關(guān)鍵詞:土地利用

    侯芳 李曉暉 胡文慧 李光永

    摘 要:本研究依托GIS,構(gòu)建張峰水庫(kù)上游流域SWAT(Soil and Water Assessment Tool)模型。基于構(gòu)建的模型,創(chuàng)新性采用從上游到下游依次對(duì)孔家坡、飛嶺、張峰3水文站1993~2008年的月徑流進(jìn)行多站點(diǎn)校準(zhǔn)和驗(yàn)證,以1993~2002年為校準(zhǔn)期,2003~2008年為驗(yàn)證期。3站點(diǎn)校準(zhǔn)期和驗(yàn)證期Nash-Sutcliffe效率系數(shù)(Ens)均大于0.70,相關(guān)性系數(shù)(R2)達(dá)到0.80以上,表明模型能夠適用于集水區(qū)徑流過(guò)程?;谀M結(jié)果表明,耕地單位面積產(chǎn)流比林地和草地大,是主要的產(chǎn)流土地利用方式;此外土壤、坡度對(duì)產(chǎn)流影響也非常大,因此土地利用方式、土壤、坡度綜合,導(dǎo)致飛嶺以上區(qū)域產(chǎn)流大于飛嶺以下。

    關(guān)鍵詞:SWAT模型;張峰水庫(kù);徑流模擬;土地利用

    中圖分類號(hào):P332 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):2095-2945(2019)01-0064-03

    Abstract: This study relies on GIS to build the SWAT (Soil and Water Assessment Tool) model of the upstream basin of Zhangfeng Reservoir. Based on the model, the monthly runoff of Kongjiapo, Feiling and Zhangfeng hydrometric stations from 1993 to 2008 is calibrated and verified from upstream to downstream. The calibration period is 1993-2002, and the validation period is 2003-2008. The Nash-Sutcliffe efficiency coefficients (Ens) of the calibration period and validation period of three stations are respectively more than 0.70 and the correlation coefficient(R2) is higher than 0.80, which indicates that the model can be applied to runoff process in catchment area. Based on the simulation results, the runoff yield per unit area of cultivated land is larger than that of woodland and grassland, which is the main way of land use. In addition, the influence of soil and slope on runoff production is also very large, so the land use pattern, soil and slope are integrated. The runoff in the area above Feiling is greater than that below it.

    Keywords: SWAT model; Zhangfeng Reservoir; runoff simulation; land use

    1 概述

    張峰水庫(kù)位于山西省晉城市,是黃河一級(jí)支流沁河上游一座重要的水利樞紐工程,主要供給以農(nóng)村人畜飲水、城市生活用水和工業(yè)用水,兼顧發(fā)電、防洪。由于氣候變化,城市化建設(shè)造成下墊面的改變,洪水日益頻發(fā)。有歷史文獻(xiàn)記載的洪水139次,其中張峰測(cè)站1892年洪峰流量達(dá)到4150~4810m3/s,1937年為2300~3140m3/s。

    為指導(dǎo)水庫(kù)調(diào)度,僅通過(guò)單點(diǎn)降雨徑流觀測(cè)數(shù)據(jù)分析已經(jīng)不能滿足要求。近年來(lái)采用分布式水文模型對(duì)沁河流域水文規(guī)律的研究較少,亟待基于研究流域內(nèi)土地利用、土壤、坡度等空間分布規(guī)律對(duì)產(chǎn)流的影響,從而分析入庫(kù)徑流的影響因素,進(jìn)而指導(dǎo)流域綜合治理,防洪抗災(zāi)。采用分布式水文模型,可提高評(píng)價(jià)水資源的精度[1],為水土保持提供重要理論依據(jù)。

    分布式水文模型Soil and Water Assessment Tool(SWA

    T)模型在國(guó)內(nèi)外均得到廣泛的應(yīng)用,模擬適用性較好[2~4]。本文針對(duì)張峰水庫(kù)上游流域,以GIS為平臺(tái),基于SWAT模型構(gòu)建分布式水文模型。創(chuàng)新性的采用孔家坡、飛嶺、張峰3個(gè)水文站月徑流資料,進(jìn)行多站點(diǎn)校準(zhǔn)及驗(yàn)證,既對(duì)模型的適用性上進(jìn)行實(shí)踐,更進(jìn)一步研究相關(guān)參數(shù)的取值,探討產(chǎn)流關(guān)鍵區(qū)及其影響因素,以期為當(dāng)?shù)厮帘3旨八畮?kù)調(diào)度提供依據(jù)。

    2 研究區(qū)概況

    張峰水庫(kù)總庫(kù)容5.05億m3,壩址以上干流長(zhǎng)224 km,控制面積4990km2。研究區(qū)北高南低,北部沁河河源至飛嶺以上以土石丘陵和山石林區(qū)為主,南部飛嶺至張峰為土石山區(qū),東西群山林立,溝壑縱橫,中部較低,沁河干流于中部縱貫?zāi)媳?。流域年均氣?.4~17.5°C。多年平均降水量587.3mm,多集中在6~9月。多年平均蒸發(fā)量為1000 mm,由南向北遞減,為降水量的2.3~3.7倍。研究區(qū)的徑流比降水滯后一個(gè)月,年內(nèi)分配極不均勻,7~10月占全年徑流量59%~69%。據(jù)1958~1980年徑流資料顯示,沁河來(lái)水量自上游到下游逐漸增加,每平方公里年產(chǎn)水量,上游孔家坡9.2m3,飛嶺到張峰站增加到11.18~15.04m3[5]。

    3 模型構(gòu)建

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

    模型需空間數(shù)據(jù)包括DEM圖、土地利用圖、土壤圖以及河網(wǎng)水系圖。屬性數(shù)據(jù)包括土壤、氣象、水文等數(shù)據(jù)。本文投影坐標(biāo)采用WGS 1984 UTM Zone 49N。

    本文選用流域內(nèi)及周邊的介休、陽(yáng)城、安澤、沁源4個(gè)氣象站,王村、衛(wèi)村等15個(gè)雨量站,孔家坡、飛嶺、張峰3個(gè)水文站日觀測(cè)值(見(jiàn)圖1)。數(shù)據(jù)系列為1991~2008年連續(xù)18年逐日觀測(cè)數(shù)據(jù)。

    本文采用中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心提供的2000年,1:10萬(wàn)的土地利用數(shù)據(jù)。采用中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心提供的2000年,1:100萬(wàn)的土壤數(shù)據(jù)。

    3.2 子流域和水文響應(yīng)單元?jiǎng)澐?/p>

    本文采用90m×90m分辨率的數(shù)字高程圖(DEM)。以張峰水庫(kù)為出口,將研究區(qū)劃分成67個(gè)子流域和785個(gè)水文響應(yīng)單元(HRUs)。

    3.3 參數(shù)敏感性分析及率定

    本文采用綜合LH(Latin Hypercube)抽樣方法和OAT(One factor At a Time)靈敏度分析方法的拉丁超立方隨機(jī)取樣(LH-OAT)法進(jìn)行敏感性分析。結(jié)合模型敏感性分析結(jié)果和前人研究成果[6],選取對(duì)集水區(qū)徑流影響較大的5個(gè)參數(shù)進(jìn)行率定(見(jiàn)表1)。

    3.4 模型校準(zhǔn)及驗(yàn)證

    本文采用SCS曲線數(shù)法計(jì)算地表徑流,Penman-Monteith計(jì)算潛在蒸騰蒸發(fā)量,變動(dòng)存儲(chǔ)系數(shù)法演算河道。采用孔家坡、飛嶺、張峰3個(gè)水文站進(jìn)行多站點(diǎn)法校準(zhǔn)和驗(yàn)證模型。

    利用構(gòu)建好的分布式水文模型對(duì)張峰水庫(kù)上游流域月徑流進(jìn)行模擬。1993~2002年為校準(zhǔn)期,2003~2008為驗(yàn)證期。采用Nash-Sutcliffe效率系數(shù)Ens、相關(guān)性系數(shù)R2、相對(duì)平均誤差RME評(píng)價(jià)模型的適用性。

    根據(jù)模擬結(jié)果3站點(diǎn)校準(zhǔn)期和驗(yàn)證期Nash-Sutcliffe效率系數(shù)均大于0.70。3站月實(shí)測(cè)與模擬徑流相關(guān)系數(shù)R2均達(dá)到0.80以上。相對(duì)平均誤差RME均在20%以內(nèi)(表2)。

    孔家坡、飛嶺、張峰站的月徑流模擬結(jié)果表明,在校準(zhǔn)期和驗(yàn)證期,各水文站徑流模擬值和實(shí)測(cè)值趨勢(shì)基本一致,構(gòu)建的模型可以較好的模擬張峰水庫(kù)上游月徑流過(guò)程。

    4 研究區(qū)產(chǎn)流規(guī)律分析

    4.1 不同土地利用產(chǎn)流貢獻(xiàn)

    本研究區(qū)內(nèi)城鄉(xiāng)、工礦、居民用地、水域四種土地利用占比為0.33%,占比較少。本文主要探討耕地、林地、草地3種土地利用方式的產(chǎn)流貢獻(xiàn)率(見(jiàn)表3)。

    由表3可知林地產(chǎn)流量最小,其原因是植被的存在使得林地地表有效糙率提高,同時(shí)植被的根系改良土壤結(jié)構(gòu)提高起入滲強(qiáng)度和蓄水性能,進(jìn)而有效地?cái)r蓄徑流[7]。

    4.2 不同水文年產(chǎn)流空間分布

    本文將產(chǎn)流空間分為“飛嶺以上”和“飛嶺-張峰”兩部分。在豐水年、平水年、枯水年,“飛嶺以上”產(chǎn)流均比“飛嶺-張峰”產(chǎn)流量大(見(jiàn)圖2)。

    產(chǎn)流與降雨、土地利用、土壤、坡度等因素有關(guān)。由于流域產(chǎn)流主要是由降雨引起。由表4可知,各水文年兩段流域范圍內(nèi)降雨量?jī)H相差5%以內(nèi),而產(chǎn)流量相差高達(dá)800%以上,由此可得,降雨對(duì)產(chǎn)流關(guān)鍵區(qū)影響不大,可看成降雨量在流域內(nèi)均等分布。

    從土地利用類型分析,研究區(qū)以林地、耕地、草地為主?!帮w嶺以上”林地大于“飛嶺-張峰”,耕地小于“飛嶺-張峰”,僅從土地利用方式分析,產(chǎn)流應(yīng)該增加,而實(shí)際情況卻相反。因此產(chǎn)流量不僅取決于土地利用方式。

    從坡度角度分析,“飛嶺以上”坡度在25%以上的面積占比為30.48%,大于“飛嶺-張峰”的27.80%。對(duì)于低海拔區(qū)降水主要以蒸發(fā)等形式耗散,較高海拔區(qū)更難形成徑流。隨海拔降低,產(chǎn)流能力逐漸減弱,低海拔區(qū)產(chǎn)流能力要明顯低于高海拔山區(qū)。因此導(dǎo)致“飛嶺以上”產(chǎn)流較“飛嶺-張峰”多。

    考慮土壤類型對(duì)產(chǎn)流分布的影響“飛嶺以上”主要土壤類型褐土;“飛嶺-張峰”主要土壤類型褐土占比較上游減少,且粗骨土、山地草甸土占比提升。根據(jù)土壤性質(zhì),褐土滲透性能更差,而粗骨土、草甸土滲透性能更好,因此進(jìn)一步導(dǎo)致“飛嶺以上”產(chǎn)流較“飛嶺-張峰”大。

    綜上所述,“飛嶺以上”產(chǎn)流大于“飛嶺-張峰”,主要原因:(1)“飛嶺以上”坡度大,有利于產(chǎn)流(2)“飛嶺以上”土壤滲透性能小,有利于產(chǎn)流。通過(guò)探究產(chǎn)流空間分布,本研究區(qū)的地形坡度和土壤綜合作用較土地利用方式對(duì)產(chǎn)流影響更大。

    5 結(jié)論

    本文主要結(jié)論如下:(1)本文基于ArcGIS構(gòu)建了張峰水庫(kù)上游流域SWAT模型。(2)采用孔家坡、飛嶺、張峰站的月徑流觀測(cè)值,對(duì)模型進(jìn)行了多站點(diǎn)校準(zhǔn)及驗(yàn)證。模擬值與實(shí)測(cè)值在校準(zhǔn)期和驗(yàn)證期Ens均大于0.70,R2均達(dá)到0.80以上,RME在20以內(nèi)。表明對(duì)于月徑流模擬,模型在張峰水庫(kù)上游流域具有較好的適用性。(3)通過(guò)不同土地利用產(chǎn)流貢獻(xiàn)的分析。表明林地產(chǎn)流貢獻(xiàn)率最小,耕地產(chǎn)流貢獻(xiàn)率最高。(4)研究區(qū)主要產(chǎn)流區(qū)域在飛嶺以上。這跟單位面積山區(qū)產(chǎn)流量高于平原地區(qū)有關(guān),隨海拔降低,產(chǎn)流能力逐漸減弱。此外產(chǎn)流量也受土壤、坡度等下墊面因素的影響。

    參考文獻(xiàn):

    [1]彭文英,張科利.不同土地利用產(chǎn)流產(chǎn)沙與降雨特征的關(guān)系[J].水土保持通報(bào),2001(04):25-29.

    [2]賴格英,吳敦銀,鐘業(yè)喜,等.SWAT模型的開(kāi)發(fā)與應(yīng)用進(jìn)展[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2012(03):243-251.

    [3]Bekiaris I G., Panagopoulos I N., Mimikou M. A. Application of the SWAT (Soil and Water Assessment Tool) model in the Ronnea catchment of Sweden[M]//2005:A160-A165.

    [4]張雪松,郝芳華,楊志峰,等.基于SWAT模型的中尺度流域產(chǎn)流產(chǎn)沙模擬研究[J].水土保持研究,2003(04):38-42.

    [5]沁河水資源利用規(guī)劃報(bào)告[M].水利部黃河水利委員會(huì)勘測(cè)規(guī)劃設(shè)計(jì)院,1991.

    [6]林炳青,陳瑩,陳興偉.SWAT模型水文過(guò)程參數(shù)區(qū)域差異研究[J].自然資源學(xué)報(bào),2013(11):1988-1999.

    [7]肖金強(qiáng),張志強(qiáng),武軍.坡面尺度林地植被對(duì)地表徑流與土壤水分的影響初步研究[J].水土保持研究,2006(05):227-231.

    猜你喜歡
    土地利用
    略談鄉(xiāng)鎮(zhèn)土地利用效益
    土地利用規(guī)劃的環(huán)境影響評(píng)價(jià)
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    土地利用規(guī)劃在房地產(chǎn)宏觀調(diào)控中的作用探討
    濱??h土地利用挖潛方向在哪里
    近20年廈門(mén)市土地利用動(dòng)態(tài)變化及預(yù)測(cè)
    基于ArcGIS 10的土地利用總體規(guī)劃圖斑自動(dòng)化綜合
    重慶市村級(jí)土地利用規(guī)劃編制實(shí)踐探討
    龍泉市MCDA土地利用功能分區(qū)
    基于PSR概念模型的土地利用系統(tǒng)健康評(píng)價(jià)
    婷婷色麻豆天堂久久| 久久午夜福利片| 少妇被粗大猛烈的视频| a级毛片在线看网站| 99国产精品免费福利视频| 亚洲av成人精品一区久久| 中文字幕人妻熟人妻熟丝袜美| 丝袜美足系列| 一本—道久久a久久精品蜜桃钙片| 欧美精品一区二区免费开放| 2022亚洲国产成人精品| 午夜av观看不卡| 99热国产这里只有精品6| 欧美日韩视频精品一区| 国产成人免费观看mmmm| 国产男人的电影天堂91| 777米奇影视久久| 色视频在线一区二区三区| 少妇猛男粗大的猛烈进出视频| 黑人欧美特级aaaaaa片| 啦啦啦中文免费视频观看日本| 91久久精品国产一区二区成人| 久久 成人 亚洲| 亚洲高清免费不卡视频| 美女cb高潮喷水在线观看| 高清欧美精品videossex| 久久人人爽人人片av| 老司机亚洲免费影院| 丰满少妇做爰视频| 一本大道久久a久久精品| 成人影院久久| 精品国产乱码久久久久久小说| 制服诱惑二区| 国产精品无大码| 国产一区二区三区av在线| 视频在线观看一区二区三区| 亚洲成色77777| 国产av精品麻豆| 婷婷色麻豆天堂久久| 亚洲精品日韩av片在线观看| 亚洲一级一片aⅴ在线观看| 精品一区在线观看国产| 多毛熟女@视频| 又粗又硬又长又爽又黄的视频| 少妇人妻久久综合中文| 久久久久久伊人网av| 免费不卡的大黄色大毛片视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 九色成人免费人妻av| 色吧在线观看| a级毛片免费高清观看在线播放| 亚洲av电影在线观看一区二区三区| 国产成人免费无遮挡视频| 亚洲国产色片| 中文字幕免费在线视频6| 91成人精品电影| 亚洲国产精品999| 99国产综合亚洲精品| 91aial.com中文字幕在线观看| 成人国语在线视频| 国产日韩一区二区三区精品不卡 | 菩萨蛮人人尽说江南好唐韦庄| 国产伦精品一区二区三区视频9| 激情五月婷婷亚洲| 久久国产亚洲av麻豆专区| 99久久精品国产国产毛片| 亚洲五月色婷婷综合| 在线观看免费日韩欧美大片 | 性色avwww在线观看| 亚洲精品成人av观看孕妇| 乱码一卡2卡4卡精品| 国产不卡av网站在线观看| 国产免费现黄频在线看| 日本黄大片高清| 成人国产麻豆网| av又黄又爽大尺度在线免费看| 伦理电影免费视频| 91在线精品国自产拍蜜月| 国产老妇伦熟女老妇高清| 亚洲精品456在线播放app| 免费看av在线观看网站| 久久久a久久爽久久v久久| 高清在线视频一区二区三区| 久久久精品区二区三区| 久久99热6这里只有精品| 亚洲人与动物交配视频| 免费观看性生交大片5| 91在线精品国自产拍蜜月| 卡戴珊不雅视频在线播放| 亚洲国产精品专区欧美| 少妇高潮的动态图| 两个人免费观看高清视频| av在线app专区| 纯流量卡能插随身wifi吗| 久久99热6这里只有精品| 国产精品久久久久久久久免| 亚洲不卡免费看| 老熟女久久久| 国产成人精品一,二区| 亚洲精品一区蜜桃| av在线老鸭窝| 国产日韩欧美亚洲二区| 国产高清三级在线| 晚上一个人看的免费电影| 亚洲成人一二三区av| 99久国产av精品国产电影| 少妇的逼好多水| 国产国拍精品亚洲av在线观看| 国产一区亚洲一区在线观看| 爱豆传媒免费全集在线观看| 久久久久视频综合| 自拍欧美九色日韩亚洲蝌蚪91| 黑人高潮一二区| av又黄又爽大尺度在线免费看| 国产永久视频网站| 久久久久久久久久久免费av| 一个人看视频在线观看www免费| 国产伦理片在线播放av一区| 人妻一区二区av| 大香蕉久久网| 免费观看av网站的网址| 亚洲人成77777在线视频| 在线观看人妻少妇| av视频免费观看在线观看| 精品久久久噜噜| 五月开心婷婷网| freevideosex欧美| 免费高清在线观看日韩| 人妻一区二区av| 卡戴珊不雅视频在线播放| 91精品伊人久久大香线蕉| 一本一本综合久久| 18禁动态无遮挡网站| 在线观看人妻少妇| 热99久久久久精品小说推荐| 大又大粗又爽又黄少妇毛片口| 久久久精品免费免费高清| 亚洲精品一二三| 日韩欧美一区视频在线观看| 亚洲第一av免费看| 晚上一个人看的免费电影| 99久国产av精品国产电影| 丝袜美足系列| 久久人人爽人人爽人人片va| 免费不卡的大黄色大毛片视频在线观看| 日韩中文字幕视频在线看片| 中文精品一卡2卡3卡4更新| 日日撸夜夜添| 亚洲伊人久久精品综合| 国产伦精品一区二区三区视频9| 欧美日韩精品成人综合77777| 国产精品成人在线| 又大又黄又爽视频免费| 观看av在线不卡| 久久精品久久久久久噜噜老黄| 一区在线观看完整版| 国产极品粉嫩免费观看在线 | 国产午夜精品久久久久久一区二区三区| 欧美亚洲 丝袜 人妻 在线| 久久精品国产亚洲av涩爱| 国产精品不卡视频一区二区| 精品酒店卫生间| 亚洲在久久综合| 精品国产一区二区三区久久久樱花| 欧美三级亚洲精品| 色5月婷婷丁香| 国产高清不卡午夜福利| 91国产中文字幕| 欧美bdsm另类| 七月丁香在线播放| 久久久国产一区二区| 亚洲综合色惰| 各种免费的搞黄视频| 久久久久人妻精品一区果冻| 久久97久久精品| 亚洲精品日本国产第一区| 老司机影院成人| 黄色视频在线播放观看不卡| 欧美精品亚洲一区二区| 飞空精品影院首页| 亚洲三级黄色毛片| 午夜免费男女啪啪视频观看| 色网站视频免费| 母亲3免费完整高清在线观看 | 日本午夜av视频| 日产精品乱码卡一卡2卡三| 一区二区三区精品91| 人人妻人人澡人人看| 国产亚洲最大av| av卡一久久| 2021少妇久久久久久久久久久| 天堂俺去俺来也www色官网| 久久99热6这里只有精品| 999精品在线视频| 久久精品国产亚洲网站| 18禁动态无遮挡网站| 成人综合一区亚洲| 97超碰精品成人国产| 少妇高潮的动态图| 交换朋友夫妻互换小说| 日本色播在线视频| 亚洲av福利一区| 日韩强制内射视频| 精品卡一卡二卡四卡免费| 丝袜脚勾引网站| 成年av动漫网址| 人妻一区二区av| 亚州av有码| 3wmmmm亚洲av在线观看| 18禁裸乳无遮挡动漫免费视频| 激情五月婷婷亚洲| 亚洲欧美清纯卡通| 性色av一级| 午夜福利视频在线观看免费| 日日摸夜夜添夜夜添av毛片| 亚洲欧美一区二区三区黑人 | 亚洲精品av麻豆狂野| 五月天丁香电影| 精品一区二区三区视频在线| 国产一级毛片在线| 国产亚洲av片在线观看秒播厂| 伊人久久国产一区二区| 欧美变态另类bdsm刘玥| 大香蕉久久网| 国产成人精品无人区| 婷婷色综合大香蕉| 久久久精品94久久精品| 好男人视频免费观看在线| 夜夜骑夜夜射夜夜干| 晚上一个人看的免费电影| h视频一区二区三区| 夫妻性生交免费视频一级片| 黄片无遮挡物在线观看| 欧美国产精品一级二级三级| 高清不卡的av网站| 成人漫画全彩无遮挡| 国产综合精华液| 久久97久久精品| 亚洲精品日本国产第一区| 欧美成人午夜免费资源| 十八禁高潮呻吟视频| 久久av网站| 欧美三级亚洲精品| 久久久亚洲精品成人影院| 各种免费的搞黄视频| videossex国产| 亚洲av综合色区一区| 国产精品久久久久久久久免| 午夜精品国产一区二区电影| 中文字幕免费在线视频6| 精品人妻一区二区三区麻豆| 十八禁网站网址无遮挡| 交换朋友夫妻互换小说| 最新中文字幕久久久久| 欧美国产精品一级二级三级| av不卡在线播放| 在线亚洲精品国产二区图片欧美 | 只有这里有精品99| 在线观看国产h片| 久久国产亚洲av麻豆专区| 婷婷色av中文字幕| 精品亚洲成a人片在线观看| 欧美人与善性xxx| 如日韩欧美国产精品一区二区三区 | 欧美97在线视频| 日本91视频免费播放| 国产在线视频一区二区| 人妻夜夜爽99麻豆av| 少妇人妻 视频| 久久久精品94久久精品| 亚洲国产最新在线播放| 特大巨黑吊av在线直播| 国产精品不卡视频一区二区| 国产黄色免费在线视频| 中文字幕最新亚洲高清| 色婷婷久久久亚洲欧美| 久久精品国产a三级三级三级| 欧美激情 高清一区二区三区| 日韩一本色道免费dvd| 成人国产av品久久久| 国产日韩欧美在线精品| 黄色配什么色好看| 国产 一区精品| 大码成人一级视频| 99久久综合免费| 成人二区视频| 黄色一级大片看看| 亚洲精品,欧美精品| 欧美少妇被猛烈插入视频| 久久久久网色| 啦啦啦在线观看免费高清www| 国产综合精华液| 国产乱来视频区| 国产精品不卡视频一区二区| 母亲3免费完整高清在线观看 | 亚洲国产最新在线播放| 亚洲不卡免费看| 久久女婷五月综合色啪小说| 日韩精品有码人妻一区| 视频中文字幕在线观看| 国产精品嫩草影院av在线观看| 精品亚洲乱码少妇综合久久| 99re6热这里在线精品视频| 三级国产精品欧美在线观看| 人人妻人人澡人人看| 久久久久久久精品精品| 国产成人aa在线观看| 亚洲国产毛片av蜜桃av| 黄色配什么色好看| 国产精品无大码| 校园人妻丝袜中文字幕| 日韩 亚洲 欧美在线| 秋霞伦理黄片| 日韩一本色道免费dvd| 免费久久久久久久精品成人欧美视频 | av视频免费观看在线观看| 寂寞人妻少妇视频99o| 人人妻人人澡人人看| 国产成人a∨麻豆精品| 欧美日韩亚洲高清精品| 亚洲成色77777| 成人国产av品久久久| 成年人免费黄色播放视频| 日韩人妻高清精品专区| 亚洲综合色网址| 高清不卡的av网站| 国产精品久久久久久久久免| 国产成人freesex在线| 最黄视频免费看| 精品少妇久久久久久888优播| 热re99久久精品国产66热6| 一个人看视频在线观看www免费| 大香蕉久久成人网| 丰满迷人的少妇在线观看| 国产无遮挡羞羞视频在线观看| 国产精品99久久久久久久久| 观看av在线不卡| 一级a做视频免费观看| av不卡在线播放| 又黄又爽又刺激的免费视频.| 天堂中文最新版在线下载| 国产男人的电影天堂91| 国产精品熟女久久久久浪| 97在线人人人人妻| 欧美激情国产日韩精品一区| 国产黄色免费在线视频| √禁漫天堂资源中文www| 亚洲精品第二区| 色5月婷婷丁香| 制服诱惑二区| 亚洲精品av麻豆狂野| 亚洲不卡免费看| 美女脱内裤让男人舔精品视频| 国产成人免费观看mmmm| 亚洲av免费高清在线观看| 蜜桃在线观看..| 天天躁夜夜躁狠狠久久av| 99九九线精品视频在线观看视频| 久久久久久久久久成人| 免费观看a级毛片全部| 久久97久久精品| 如何舔出高潮| 大香蕉97超碰在线| 国产一区二区在线观看日韩| 一级毛片aaaaaa免费看小| 欧美+日韩+精品| 母亲3免费完整高清在线观看 | 成人影院久久| 一级毛片 在线播放| 国产老妇伦熟女老妇高清| 免费人成在线观看视频色| 91久久精品电影网| 桃花免费在线播放| 久久精品久久精品一区二区三区| 麻豆精品久久久久久蜜桃| 亚洲在久久综合| 我的女老师完整版在线观看| 久久久久国产精品人妻一区二区| 久久亚洲国产成人精品v| 欧美日韩视频精品一区| 最近手机中文字幕大全| 午夜福利视频在线观看免费| 日本黄色日本黄色录像| 亚洲熟女精品中文字幕| 亚洲婷婷狠狠爱综合网| 国产亚洲av片在线观看秒播厂| 纵有疾风起免费观看全集完整版| 九色成人免费人妻av| 丝袜脚勾引网站| 日日摸夜夜添夜夜添av毛片| 亚洲欧洲国产日韩| 国产在视频线精品| 夜夜爽夜夜爽视频| 亚洲国产av新网站| 国产亚洲精品第一综合不卡 | 高清不卡的av网站| 免费高清在线观看日韩| 一区二区三区乱码不卡18| 免费高清在线观看日韩| 久久国内精品自在自线图片| 亚洲欧美精品自产自拍| 国产片特级美女逼逼视频| 在线看a的网站| 国产精品嫩草影院av在线观看| 亚洲欧美精品自产自拍| 51国产日韩欧美| 男女边摸边吃奶| 热99国产精品久久久久久7| 欧美人与善性xxx| 国产精品国产三级国产av玫瑰| 欧美日韩av久久| 夫妻性生交免费视频一级片| 一区二区三区免费毛片| 中文字幕人妻丝袜制服| a级毛片黄视频| 91久久精品电影网| 国产视频内射| 免费观看a级毛片全部| 飞空精品影院首页| 久久久久久久久大av| 国产亚洲精品第一综合不卡 | 一本一本综合久久| 永久免费av网站大全| 51国产日韩欧美| 在线观看免费日韩欧美大片 | 十分钟在线观看高清视频www| 国产亚洲av片在线观看秒播厂| av又黄又爽大尺度在线免费看| 欧美一级a爱片免费观看看| 午夜激情久久久久久久| 赤兔流量卡办理| 国产精品一区二区在线不卡| 久久精品国产亚洲av涩爱| 国产高清有码在线观看视频| 亚洲经典国产精华液单| 春色校园在线视频观看| 成人毛片a级毛片在线播放| 水蜜桃什么品种好| 国产黄频视频在线观看| 中文字幕亚洲精品专区| 美女国产视频在线观看| 日韩视频在线欧美| 国产视频首页在线观看| 一个人看视频在线观看www免费| av国产精品久久久久影院| 女性生殖器流出的白浆| 国产在视频线精品| 一级,二级,三级黄色视频| 人妻系列 视频| 3wmmmm亚洲av在线观看| 精品人妻一区二区三区麻豆| 一级a做视频免费观看| 久久毛片免费看一区二区三区| 国产片内射在线| 国产在线视频一区二区| 国模一区二区三区四区视频| 国产毛片在线视频| 午夜激情久久久久久久| 精品亚洲成国产av| 亚洲av中文av极速乱| 久久国产亚洲av麻豆专区| 国产深夜福利视频在线观看| 欧美日韩成人在线一区二区| 久久久精品区二区三区| 18+在线观看网站| 天美传媒精品一区二区| 波野结衣二区三区在线| 五月天丁香电影| 亚洲激情五月婷婷啪啪| 国产一级毛片在线| 街头女战士在线观看网站| 春色校园在线视频观看| 日韩熟女老妇一区二区性免费视频| 中文字幕最新亚洲高清| 亚洲人成网站在线播| 亚洲精品aⅴ在线观看| 国产免费一区二区三区四区乱码| 赤兔流量卡办理| 各种免费的搞黄视频| 欧美亚洲日本最大视频资源| 永久网站在线| 成人影院久久| 久久国产精品男人的天堂亚洲 | 精品人妻偷拍中文字幕| 久久97久久精品| 在线观看人妻少妇| 免费观看a级毛片全部| 国产69精品久久久久777片| 久久精品国产自在天天线| 涩涩av久久男人的天堂| 伦理电影免费视频| 国产成人精品福利久久| 国产老妇伦熟女老妇高清| 成人手机av| 十八禁高潮呻吟视频| 大片电影免费在线观看免费| 伊人久久精品亚洲午夜| 亚洲四区av| 欧美精品一区二区大全| 成人二区视频| 我的女老师完整版在线观看| 97超碰精品成人国产| 国产成人免费无遮挡视频| 亚洲国产精品一区二区三区在线| a级毛色黄片| 九九久久精品国产亚洲av麻豆| 久久婷婷青草| 亚洲精品久久午夜乱码| 国产有黄有色有爽视频| 秋霞在线观看毛片| 搡女人真爽免费视频火全软件| 亚洲国产毛片av蜜桃av| 亚洲av电影在线观看一区二区三区| 欧美最新免费一区二区三区| 黄色视频在线播放观看不卡| 日本午夜av视频| 天天影视国产精品| 国产视频首页在线观看| 午夜视频国产福利| 国产欧美另类精品又又久久亚洲欧美| videos熟女内射| 亚洲熟女精品中文字幕| 国产精品国产三级国产专区5o| 日韩av不卡免费在线播放| 夜夜骑夜夜射夜夜干| 国产成人精品福利久久| 国产av精品麻豆| 中国国产av一级| 国产色婷婷99| 桃花免费在线播放| 中文字幕人妻熟人妻熟丝袜美| 哪个播放器可以免费观看大片| 99热这里只有是精品在线观看| 久久久午夜欧美精品| 少妇 在线观看| 寂寞人妻少妇视频99o| 黄色视频在线播放观看不卡| 日产精品乱码卡一卡2卡三| 日韩 亚洲 欧美在线| 亚洲av国产av综合av卡| 国产精品一二三区在线看| av播播在线观看一区| 蜜桃久久精品国产亚洲av| 欧美bdsm另类| 久热久热在线精品观看| 免费观看a级毛片全部| 老熟女久久久| 亚洲av男天堂| 国产精品麻豆人妻色哟哟久久| 有码 亚洲区| 赤兔流量卡办理| 一级二级三级毛片免费看| 国产一区二区三区综合在线观看 | 久久ye,这里只有精品| 日韩在线高清观看一区二区三区| 在线观看三级黄色| 在线观看免费视频网站a站| 国产精品三级大全| 国产精品99久久久久久久久| 国产黄片视频在线免费观看| 少妇高潮的动态图| 午夜老司机福利剧场| 三级国产精品欧美在线观看| 黄色视频在线播放观看不卡| 18在线观看网站| 婷婷成人精品国产| 国产熟女午夜一区二区三区 | 五月伊人婷婷丁香| 蜜桃久久精品国产亚洲av| 国产精品国产三级国产av玫瑰| 午夜免费观看性视频| 天堂8中文在线网| 性高湖久久久久久久久免费观看| av国产精品久久久久影院| 99国产综合亚洲精品| 免费日韩欧美在线观看| 精品人妻偷拍中文字幕| 在线免费观看不下载黄p国产| 亚洲精品日韩在线中文字幕| 免费看av在线观看网站| 最新中文字幕久久久久| 中文字幕免费在线视频6| 成人国产麻豆网| 老女人水多毛片| 在线观看三级黄色| 国产午夜精品一二区理论片| 国产高清国产精品国产三级| 人妻制服诱惑在线中文字幕| 18禁在线无遮挡免费观看视频| 亚洲成人av在线免费| 美女国产高潮福利片在线看| 男女啪啪激烈高潮av片| 嘟嘟电影网在线观看| 女人精品久久久久毛片| 日韩视频在线欧美| 一级二级三级毛片免费看| 亚洲中文av在线| 国语对白做爰xxxⅹ性视频网站| 边亲边吃奶的免费视频| 女性被躁到高潮视频| 2022亚洲国产成人精品| 性色av一级| 国产 一区精品| 日本91视频免费播放| 精品久久蜜臀av无| 免费看不卡的av| 一区二区三区四区激情视频| 亚洲成人av在线免费| 人妻少妇偷人精品九色| 久久99蜜桃精品久久| 免费黄网站久久成人精品| 我的女老师完整版在线观看| 人妻 亚洲 视频| 国产高清有码在线观看视频| 国产精品久久久久久久电影|