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

    基于SWAT模型的天目湖流域徑流模擬分析

    2022-07-05 09:04:40許延昌
    地下水 2022年3期
    關(guān)鍵詞:天目湖徑流土地利用

    許延昌

    (1.江蘇海洋大學(xué)海洋技術(shù)與測繪學(xué)院,江蘇 連云港 222000;2.中國科學(xué)院南京地理與湖泊研究所,江蘇 南京 210000)

    0 引言

    隨著社會(huì)經(jīng)濟(jì)的快速發(fā)展,人們的生活水平日益提高,但出現(xiàn)的問題也日趨明顯,其中比較突出的一個(gè)就是水環(huán)境問題[1]。近幾年來,天目湖流域隨著種植業(yè)、水產(chǎn)養(yǎng)殖業(yè)等產(chǎn)業(yè)的蓬勃發(fā)展,產(chǎn)生的污染物致使流域內(nèi)很多湖泊、河流等都受到了不同程度的污染,出現(xiàn)了藻類大量繁殖、富營養(yǎng)化等現(xiàn)象,甚至有進(jìn)一步惡化的趨勢,這就亟需要一種合適的方法對該流域的污染情況等進(jìn)行分析,從而提出有效的對策。在研究水環(huán)境的眾多方法中,SWAT( Soil and Water Assessment Tool)水文模型模擬分析法應(yīng)用較為廣泛。Woo So-Young[2]等利用 SWAT 模型,選取 8 種水量水質(zhì)相關(guān)情景,分析了 IBWT(跨流域調(diào)水)對受納流域水質(zhì)的影響; M’Hamed Boufala[3]等采用水農(nóng)SWAT模型,對M’dez流域等高線耕作、梯田等最佳管理措施( BMPs )進(jìn)行評價(jià),量化了水土平衡和土壤侵蝕及其空間分布; 馬新萍[4]等基于 SWAT 模型預(yù)測了未來 2 種土地利用變化模式下的徑流量變化。故本次研究選擇SWAT模型,結(jié)合天目湖流域土地利用等相關(guān)數(shù)據(jù),對該流域的徑流進(jìn)行了模擬。

    1 研究區(qū)概況

    天目湖流域主要位于江蘇省溧陽市南部,小部分位于安徽廣德、郎溪境內(nèi),介于119°18′E~119°29′E、31°8′N~31°25′N之間,見圖1,地勢整體上南高北低。其中的天目湖是國家5A級景區(qū)、生態(tài)旅游示范區(qū),流域面積246 km2,有大溪、沙河兩座國家級水庫。流域內(nèi)有多種地貌類型,且四季分明,物種豐富,雨量充沛,有包括水體、林地、茶園、農(nóng)田等多種覆被類型[5]。

    圖1 研究區(qū)示意圖

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

    2.1 數(shù)據(jù)來源

    SWAT模型的構(gòu)建主要需要兩種數(shù)據(jù),即屬性數(shù)據(jù)和空間數(shù)據(jù)[6]。本次研究所使用的空間數(shù)據(jù)有DEM高程數(shù)據(jù)、土地利用數(shù)據(jù)、土壤數(shù)據(jù)等,數(shù)據(jù)詳細(xì)情況見表1;屬性數(shù)據(jù)主要用到的是所在研究區(qū)的降水、風(fēng)速、相對濕度、溫度以及太陽輻射等氣象數(shù)據(jù)。

    表1 數(shù)據(jù)說明

    2.2 研究方法

    SWAT模型是由美國農(nóng)業(yè)部農(nóng)業(yè)研究服務(wù)中心( USDA-ARS) 開發(fā)的一種基于流域尺度的分布式水文模型,該模型的原理是將參與水文循環(huán)模擬的各要素的變化進(jìn)行模擬并分析,模型的物理基礎(chǔ)較強(qiáng),能綜合考慮土壤、土地利用等數(shù)據(jù)來模擬多種復(fù)雜的情況[7-8]。模型的率定及驗(yàn)證使用的是SWAT-CUP 2012軟件,通過實(shí)測徑流數(shù)據(jù)對模擬數(shù)據(jù)進(jìn)行率定,通過Sensitivity analysis模塊中的Global Sensitivity功能進(jìn)行敏感性分析,挑選出對模型模擬結(jié)果影響最敏感的參數(shù),用新的參數(shù)及參數(shù)范圍重復(fù)率定,直到滿足要求為止。將率定好的參數(shù)帶回SWAT模型,重新寫入數(shù)據(jù)再次模擬,用實(shí)測徑流數(shù)據(jù)驗(yàn)證模擬結(jié)果。本次研究選擇的評價(jià)指標(biāo)是決定系數(shù)R2和Nash-Suttcliffe效率系數(shù)NS,公式如下。式中:Qm為觀測流量(m3/s),Qs為模擬流量(m3/s),Qavgm為觀測流量平均值 (m3/s) ,Qavgs為模擬流量平均值(m3/s)[9]。

    (1)

    (2)

    3 模型構(gòu)建與模擬

    3.1 流域劃分

    在地理空間數(shù)據(jù)云下載天目湖流域分辨率為30 m的DEM數(shù)據(jù),在ArcSWAT的Watershed Delineation(流域劃分)功能模塊中進(jìn)行加載。由于洼地的存在,水流方向、匯流累積量的計(jì)算以及河道的生成都會(huì)受到影響,所以要進(jìn)行填洼預(yù)處理。填洼后要確定最小流域面積,面積越小生成的子流域越多,反之生成的子流域則越少,所以選擇一個(gè)合適的最小流域面積至關(guān)重要,經(jīng)多次測試,發(fā)現(xiàn)最小面積設(shè)置在1 000 hm2是比較合適的,然后選擇流域總出口后即可開始子流域的劃分,最后計(jì)算子流域的參數(shù)。

    3.2 水文響應(yīng)單元的定義

    水文響應(yīng)單元(Hydrologic Response Unit,HRU),由土地利用數(shù)據(jù)、土壤數(shù)據(jù)以及坡度數(shù)據(jù)共同組合而成。根據(jù)二調(diào)土地利用SWAT分類對照表,將土地利用類型分為耕地 (AGRL) 、林地 (FRST) 、牧場 (PAST) 、水域 (WATR) 、居民區(qū) (URLD) 以及荒地裸地 (BARR),見圖2(受dem影響,研究區(qū)北部一小部分區(qū)域在ArcSWAT生成流域時(shí)未生成,下同)。用天目湖流域的掩膜數(shù)據(jù)提取HWSD土壤數(shù)據(jù)庫,流域內(nèi)共有4種土壤類型,其中人為堆積土(Cumulic Anthrosols)占比30.369%,簡育低活性強(qiáng)酸土(Haplic Acrisols)占比14.644%,簡育高活性淋溶土(Haplic Luvisols)占比47.653%,水體(Water bodies)占比7.334%,見圖3。坡度代表地形陡緩的程度,用坡面的垂直高度與水平距離的比來計(jì)算,參考坡度的統(tǒng)計(jì)參數(shù):平均值為14.20、標(biāo)準(zhǔn)偏差為11.58,將坡度分為三級:0~9、9~15、15~9999,見圖4。加載并設(shè)置完三種類型的數(shù)據(jù)后分別對其進(jìn)行重分類,后續(xù)計(jì)算中會(huì)根據(jù)分類代碼從數(shù)據(jù)庫中檢索相關(guān)屬性值進(jìn)行計(jì)算。

    HRU的定義有四種方式:Dominant Land Use,Soils,Slope、Dominant HRU、Target Number of HRUs和Multiple HRUs,本次研究選用使用頻率最高的第四種,該方式會(huì)通過設(shè)置閾值的方式,將小于閾值的部分按比重分配到大于該閾值的土地利用、土壤和坡度數(shù)據(jù)中[10]。通過面積比重的比較,將土地利用、土壤和坡度數(shù)據(jù)的閾值分別設(shè)置為:1%、0%、0%,最終將研究區(qū)劃分為385個(gè)HRU。

    圖2 土地利用空間分布圖 圖3 土壤空間分布圖 圖4 坡度空間分布圖

    3.3 氣象數(shù)據(jù)的定義

    CMADS(The China Meteorological Assimilation Driving Datasets for the SWAT model)數(shù)據(jù)集引入了LAPSSTMAS算法[11],采用雙線性插值、數(shù)據(jù)循環(huán)嵌套以及重采樣等多種科學(xué)方法和技術(shù)手段建立,并且針對SWAT模型進(jìn)行了格式等方面的一系列修改,使得SWAT模型的研究者可以直接使用本套中國大氣同化驅(qū)動(dòng)數(shù)據(jù)集而不需要做任何格式上的處理[12],這降低了數(shù)據(jù)集使用的難度,減少了研究人員的工作量,也使得CMADS的應(yīng)用變的更加廣泛。本次研究使用的是CMADS V1.2版本,空間分辨率為1/8°*1/8°,時(shí)間分辨率為逐日,時(shí)間尺度為2008年至2018年,主要包括日平均相對濕度 (Relative-Humidity) 、日平均風(fēng)速 (Wind) 、 日24h累計(jì)降水量 (Precipitation) 、日最高最低氣溫 (Temperature) 以及日平均太陽輻射 (Solar radiation) 等數(shù)據(jù)。

    3.4 模型模擬

    目前,CMADS公布的最新氣象數(shù)據(jù)是從2008年1月1日至2018年12月31日,故本次模擬也在該時(shí)間段內(nèi)。為了使模擬結(jié)果更加準(zhǔn)確,設(shè)置2008年為預(yù)熱期進(jìn)行預(yù)熱,然后對天目湖流域2009至2018年的徑流進(jìn)行模擬。

    3.5 率定與驗(yàn)證

    3.5.1 模型率定

    SWAT-CUP2012軟件是專為SWAT模型的率定、驗(yàn)證工作開發(fā)的計(jì)算機(jī)軟件,該軟件集成了SUFI-2、MCMC、GLUE、ParaSol和PSO程序,可用于SWAT模型的不確定性分析、敏感性分析、率定和驗(yàn)證等。為使率定的效果更好,本次研究使用的是最常用的SUFI-2算法,該算法通過拉丁超立方體抽樣法判斷參數(shù)的敏感性并通過95PPU(95%置信水平上的不確定性區(qū)間帶)來表示參數(shù)的不確定性[13]。

    SWAT-CUP里的Sensitivity analysis模塊中的Global Sensitivity功能可進(jìn)行參數(shù)的敏感性分析,由t-Stat和P-Value兩個(gè)參數(shù)來判定,t-Stat代表敏感性的程度,絕對值越大說明參數(shù)越敏感,P-Value表示了敏感性的顯著性,值越接近于零就越顯著。率定數(shù)據(jù)采用的是2010年,2011年和2012年的實(shí)測徑流數(shù)據(jù),經(jīng)率定及敏感性分析,測得ALPHA_BF(基流消退系數(shù))、GW_REVAP(地下水再蒸發(fā)系數(shù))等8個(gè)參數(shù)對天目湖流域徑流的模擬結(jié)果影響最大,見表2。本次研究選擇的評價(jià)指標(biāo)是決定系數(shù)R2和Nash-Suttcliffe效率系數(shù)NS,其中R2反映的是模擬值與實(shí)測值的吻合程度,取值范圍是[0,1],當(dāng)R2>0.5時(shí)說明模擬結(jié)果是合理的,另一個(gè)評價(jià)指標(biāo)NS反映的是觀測徑流過程與模擬徑流過程間的擬合程度,取值范圍是(-∞,1],當(dāng)NS>0.5時(shí)表明模擬結(jié)果是滿足要求的。本次率定的R2為0.68,NS為0.66,說明率定效果滿足要求,率定效果圖見圖5。

    表2 敏感參數(shù)說明

    3.5.2 模型驗(yàn)證

    率定完成后,將SWAT-CUP軟件校準(zhǔn)輸出中Best_Par里敏感參數(shù)的Fitted_Value代回SWAT模型中。ArcSWAT中Manual Calibration功能中有三種參數(shù)校正方法,分別是Multiply By(r)、Add(a)和Replace Value(v)。Multiply By運(yùn)算是將原參數(shù)值乘以(Fitted_Value + 1),Add運(yùn)算是將原參數(shù)值加Fitted_Value,Replace Value是將原參數(shù)值用Fitted_Value替換。

    參數(shù)校正完成后重新寫入輸入文件中,再次執(zhí)行模擬,用2016年和2018年兩年的實(shí)測徑流數(shù)據(jù)在SWAT-CUP軟件中驗(yàn)證模擬結(jié)果。本次驗(yàn)證決定系數(shù)R2為0.70,Nash-Suttcliffe效率系數(shù)NS為0.65,驗(yàn)證效果見圖6。

    圖5 率定結(jié)果 圖6 驗(yàn)證結(jié)果

    4 結(jié)語

    (1)本次研究通過DEM數(shù)據(jù)、土地利用數(shù)據(jù)、土壤數(shù)據(jù)和氣象數(shù)據(jù)等構(gòu)建了天目湖流域SWAT水文模型,并利用該流域的實(shí)測徑流數(shù)據(jù)對模擬結(jié)果進(jìn)行了率定及驗(yàn)證,結(jié)果表明:率定階段決定系數(shù)R2為0.68,Nash-Suttcliffe效率系數(shù)NS為0.66,驗(yàn)證階段R2為0.70,NS為0.65,4值均大于0.5,說明本次模擬結(jié)果是滿足要求的,構(gòu)建的SWAT模型能較好的反映天目湖流域的產(chǎn)匯流規(guī)律,可以用于以后的水文規(guī)律的預(yù)測等工作。另外,之所以會(huì)有較好的模擬結(jié)果,一個(gè)主要的原因便是SWAT模型采用了把流域劃分成下墊面相似的子流域進(jìn)行分塊計(jì)算的計(jì)算方式。

    (2)從參數(shù)敏感性分析結(jié)果中可以看出,敏感度最高的參數(shù)是基流消退系數(shù)ALPHA_BF,它是地下水對補(bǔ)給變化的一個(gè)直接指數(shù),其t-Stat值為6.36,遠(yuǎn)大于地下水再蒸發(fā)系數(shù)GW_REVAP等7個(gè)參數(shù),P-Value為0;而敏感度最低的參數(shù)是土壤飽和水力傳導(dǎo)系數(shù)SOL_K,其t-Stat值為0.52, P-Value為0.61。

    (3)率定期(2010-2012年)和驗(yàn)證期(2016年、2018年)的觀測值和模擬值能夠較好的吻合,且觀察折線圖能夠看出,每一年的 5-9月份的徑流量都明顯高于其他月份,一個(gè)很重要的原因是5-9月份降雨比較多,這符合現(xiàn)實(shí)情況,也從一個(gè)側(cè)面反映了本次模擬結(jié)果較好。

    (4)本次研究后續(xù)可以根據(jù)模擬結(jié)果中hru(水文響應(yīng)單元)數(shù)據(jù)表、rch(河道)數(shù)據(jù)表以及sub(子流域)數(shù)據(jù)表中模擬的總磷、總氮等數(shù)據(jù)進(jìn)行污染物輸移等相關(guān)分析,為找出污染天目湖流域水質(zhì)的因子及為相關(guān)部門制定水質(zhì)保護(hù)措施等方面提供幫助。

    猜你喜歡
    天目湖徑流土地利用
    美麗的天目湖
    美麗的天目湖
    天目湖日月潭演繹的佳話
    TЫCяЧEЛETний КУЛьTУрнЫй гороДоК
    中國(俄文)(2018年2期)2018-03-12 08:37:46
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    濱??h土地利用挖潛方向在哪里
    探秘“大徑流”
    攻克“大徑流”
    江埡水庫降雨徑流相關(guān)圖的建立
    在线观看三级黄色| 日韩亚洲欧美综合| 亚洲精品乱久久久久久| 国产精品福利在线免费观看| 永久网站在线| 国产在线男女| 熟妇人妻不卡中文字幕| 99热6这里只有精品| 久久综合国产亚洲精品| 亚洲av男天堂| 欧美精品一区二区免费开放| 欧美区成人在线视频| 在线观看一区二区三区激情| 热99国产精品久久久久久7| 女人十人毛片免费观看3o分钟| 又大又黄又爽视频免费| 波野结衣二区三区在线| 国产成人精品一,二区| 久久人人爽人人片av| 美女内射精品一级片tv| 日韩强制内射视频| 日韩大片免费观看网站| 国产黄色免费在线视频| 性高湖久久久久久久久免费观看| 男女无遮挡免费网站观看| 久久久久久久国产电影| 国产人妻一区二区三区在| 丰满少妇做爰视频| 亚洲欧美日韩卡通动漫| 国产精品三级大全| 丰满人妻一区二区三区视频av| 国产一区二区在线观看日韩| 在线观看免费日韩欧美大片 | 久久久久视频综合| 观看av在线不卡| 黄色日韩在线| 久久久欧美国产精品| 一级毛片aaaaaa免费看小| 精品一区二区三卡| 国产一区有黄有色的免费视频| 亚洲图色成人| 亚洲成人一二三区av| 国产成人aa在线观看| 2021少妇久久久久久久久久久| 久久午夜福利片| 99久久综合免费| 国产免费视频播放在线视频| 久久97久久精品| 一级毛片我不卡| 亚洲国产最新在线播放| 日韩免费高清中文字幕av| 六月丁香七月| 国产av国产精品国产| 国产精品一二三区在线看| 乱系列少妇在线播放| 久久久久网色| 亚洲四区av| 久久国内精品自在自线图片| 一区二区三区免费毛片| 毛片一级片免费看久久久久| 国产av国产精品国产| 能在线免费看毛片的网站| 夜夜骑夜夜射夜夜干| 亚洲四区av| av不卡在线播放| 色婷婷久久久亚洲欧美| 精品久久久久久久末码| 91在线精品国自产拍蜜月| av国产精品久久久久影院| 丰满迷人的少妇在线观看| 男女国产视频网站| 成人特级av手机在线观看| 男女边吃奶边做爰视频| 亚洲国产欧美人成| 啦啦啦在线观看免费高清www| 十分钟在线观看高清视频www | 久久av网站| 99国产精品免费福利视频| 91狼人影院| 亚洲无线观看免费| 伦理电影大哥的女人| 成人高潮视频无遮挡免费网站| 高清视频免费观看一区二区| 熟女av电影| 国产伦理片在线播放av一区| 亚洲国产精品999| 亚洲,欧美,日韩| 国产伦在线观看视频一区| 国产成人午夜福利电影在线观看| 日日摸夜夜添夜夜添av毛片| 欧美人与善性xxx| 三级经典国产精品| av福利片在线观看| 肉色欧美久久久久久久蜜桃| 国产黄片美女视频| 精品国产露脸久久av麻豆| 精品国产露脸久久av麻豆| 全区人妻精品视频| 国产成人精品久久久久久| 精品人妻偷拍中文字幕| 国产极品天堂在线| 亚洲欧美精品专区久久| 久久精品国产自在天天线| av又黄又爽大尺度在线免费看| 多毛熟女@视频| 国产精品精品国产色婷婷| 一区二区三区四区激情视频| 好男人视频免费观看在线| av国产久精品久网站免费入址| 欧美日韩国产mv在线观看视频 | 人妻 亚洲 视频| 香蕉精品网在线| 成人二区视频| 水蜜桃什么品种好| 日本猛色少妇xxxxx猛交久久| 亚洲美女视频黄频| 国产精品一及| 777米奇影视久久| 搡老乐熟女国产| 最近中文字幕高清免费大全6| 综合色丁香网| 久久6这里有精品| 日日摸夜夜添夜夜爱| 91在线精品国自产拍蜜月| 国产综合精华液| 校园人妻丝袜中文字幕| 久久97久久精品| 免费观看性生交大片5| 中文字幕精品免费在线观看视频 | 久久久久人妻精品一区果冻| 国产精品欧美亚洲77777| 国产视频首页在线观看| 青青草视频在线视频观看| 亚洲精品国产成人久久av| 一级毛片电影观看| 深爱激情五月婷婷| 91在线精品国自产拍蜜月| 久久韩国三级中文字幕| 春色校园在线视频观看| 我要看黄色一级片免费的| 五月天丁香电影| a 毛片基地| 在线免费观看不下载黄p国产| 麻豆精品久久久久久蜜桃| 亚洲精品456在线播放app| 一级av片app| 高清在线视频一区二区三区| 日日啪夜夜爽| tube8黄色片| 99久久精品一区二区三区| 成人午夜精彩视频在线观看| 精品一品国产午夜福利视频| 午夜免费男女啪啪视频观看| 国产成人午夜福利电影在线观看| 国产av一区二区精品久久 | 欧美日韩一区二区视频在线观看视频在线| 国产有黄有色有爽视频| 亚洲欧美精品自产自拍| 日韩av在线免费看完整版不卡| 欧美精品人与动牲交sv欧美| 亚洲精品久久久久久婷婷小说| 国产真实伦视频高清在线观看| 在线观看一区二区三区| 久久久久久伊人网av| 夫妻午夜视频| 久久影院123| 日本爱情动作片www.在线观看| 欧美成人a在线观看| 色婷婷久久久亚洲欧美| 青春草国产在线视频| 国产91av在线免费观看| 久久婷婷青草| 99视频精品全部免费 在线| 久久久久久九九精品二区国产| 亚洲精品久久久久久婷婷小说| 日本午夜av视频| 少妇人妻 视频| av免费在线看不卡| 国产精品三级大全| 激情 狠狠 欧美| 欧美性感艳星| 一级二级三级毛片免费看| av国产精品久久久久影院| 毛片女人毛片| 国产成人aa在线观看| 热re99久久精品国产66热6| 如何舔出高潮| 欧美人与善性xxx| 日本av免费视频播放| 少妇裸体淫交视频免费看高清| 精品人妻熟女av久视频| 极品少妇高潮喷水抽搐| 免费看光身美女| av播播在线观看一区| 欧美国产精品一级二级三级 | 久久6这里有精品| 国产美女午夜福利| 国产日韩欧美亚洲二区| 亚洲在久久综合| 久久久久久九九精品二区国产| 日本wwww免费看| 观看美女的网站| 男人爽女人下面视频在线观看| 特大巨黑吊av在线直播| 精品亚洲成a人片在线观看 | 日本av手机在线免费观看| 在线 av 中文字幕| 亚洲av不卡在线观看| 精品久久久噜噜| 亚洲中文av在线| 亚洲精品视频女| 亚洲自偷自拍三级| 亚洲成人手机| 久久99热6这里只有精品| 久久久久久久久久久丰满| 国产精品人妻久久久久久| 交换朋友夫妻互换小说| av又黄又爽大尺度在线免费看| 超碰av人人做人人爽久久| 亚洲欧美中文字幕日韩二区| av黄色大香蕉| 国产成人精品福利久久| 亚洲真实伦在线观看| 一级毛片电影观看| 国产在视频线精品| 久久久久久久国产电影| 亚洲久久久国产精品| 日韩av不卡免费在线播放| 国产成人精品婷婷| 亚洲婷婷狠狠爱综合网| 国产无遮挡羞羞视频在线观看| 国产成人a区在线观看| 欧美区成人在线视频| 久久精品国产自在天天线| 久久青草综合色| 欧美日韩亚洲高清精品| 国产一区二区三区综合在线观看 | 午夜日本视频在线| 国产一区亚洲一区在线观看| 99九九线精品视频在线观看视频| 国产片特级美女逼逼视频| 蜜桃亚洲精品一区二区三区| 另类亚洲欧美激情| 国产免费福利视频在线观看| 色吧在线观看| 国产黄片视频在线免费观看| 国产成人aa在线观看| 亚洲av成人精品一区久久| 亚洲中文av在线| 亚洲欧美精品专区久久| 欧美精品国产亚洲| 男男h啪啪无遮挡| 国产一级毛片在线| av专区在线播放| 国产精品久久久久久久电影| av免费观看日本| 欧美极品一区二区三区四区| 精品一品国产午夜福利视频| 欧美人与善性xxx| 汤姆久久久久久久影院中文字幕| 日韩三级伦理在线观看| 天天躁日日操中文字幕| 日韩一区二区三区影片| 亚洲av中文字字幕乱码综合| 日本免费在线观看一区| 日韩伦理黄色片| 我要看日韩黄色一级片| 最近手机中文字幕大全| 国产日韩欧美亚洲二区| 国产成人精品福利久久| 99热国产这里只有精品6| 一级毛片久久久久久久久女| 亚洲国产精品成人久久小说| 亚洲国产毛片av蜜桃av| 久久女婷五月综合色啪小说| 亚洲美女视频黄频| 亚洲aⅴ乱码一区二区在线播放| 99热这里只有是精品在线观看| 午夜视频国产福利| 午夜福利网站1000一区二区三区| 在线观看三级黄色| 91久久精品国产一区二区成人| 国产老妇伦熟女老妇高清| 国产精品久久久久久精品电影小说 | 又粗又硬又长又爽又黄的视频| 91精品国产国语对白视频| 麻豆乱淫一区二区| 国产色爽女视频免费观看| 久久久久久久久久人人人人人人| 日产精品乱码卡一卡2卡三| 97超视频在线观看视频| 91精品国产国语对白视频| 成年免费大片在线观看| 亚洲丝袜综合中文字幕| 国产成人精品久久久久久| 久久久成人免费电影| 高清不卡的av网站| 少妇熟女欧美另类| 乱码一卡2卡4卡精品| 女性生殖器流出的白浆| 搡女人真爽免费视频火全软件| 97热精品久久久久久| 乱码一卡2卡4卡精品| 亚洲欧美一区二区三区黑人 | 久久久a久久爽久久v久久| 网址你懂的国产日韩在线| 26uuu在线亚洲综合色| 久久97久久精品| 日韩一区二区三区影片| 亚洲人成网站在线观看播放| 男女边摸边吃奶| 国产精品爽爽va在线观看网站| 久久人人爽人人片av| 国产伦理片在线播放av一区| 看非洲黑人一级黄片| 在线观看免费高清a一片| 久久精品国产a三级三级三级| 亚洲中文av在线| 亚洲av.av天堂| 国产成人a∨麻豆精品| 丝瓜视频免费看黄片| 国产精品一区二区在线观看99| 黄片wwwwww| 99久久综合免费| 涩涩av久久男人的天堂| 高清午夜精品一区二区三区| 日本vs欧美在线观看视频 | 99久久综合免费| 国产爽快片一区二区三区| 亚洲av二区三区四区| 22中文网久久字幕| 欧美日韩精品成人综合77777| av免费在线看不卡| 有码 亚洲区| 伦理电影免费视频| 在线观看美女被高潮喷水网站| 免费观看性生交大片5| 六月丁香七月| 国产精品一二三区在线看| 99久久人妻综合| 久久久久久久亚洲中文字幕| 国产毛片在线视频| 国产在线免费精品| 免费观看在线日韩| 自拍偷自拍亚洲精品老妇| 日韩制服骚丝袜av| 一级a做视频免费观看| 日本爱情动作片www.在线观看| 亚洲精品第二区| 51国产日韩欧美| 国产爽快片一区二区三区| 精品亚洲乱码少妇综合久久| 乱码一卡2卡4卡精品| 亚洲无线观看免费| 亚洲经典国产精华液单| 美女脱内裤让男人舔精品视频| 制服丝袜香蕉在线| 国产精品福利在线免费观看| 日韩制服骚丝袜av| 成人无遮挡网站| 观看免费一级毛片| 永久网站在线| 免费不卡的大黄色大毛片视频在线观看| 日韩欧美精品免费久久| 精品亚洲乱码少妇综合久久| av.在线天堂| 少妇的逼好多水| 欧美xxxx性猛交bbbb| 少妇的逼水好多| 只有这里有精品99| 国产黄片视频在线免费观看| 一区在线观看完整版| 高清av免费在线| 99热全是精品| 国产成人一区二区在线| 国产真实伦视频高清在线观看| 我要看日韩黄色一级片| 国产精品伦人一区二区| 人妻制服诱惑在线中文字幕| 成人国产麻豆网| 人妻少妇偷人精品九色| 中文字幕人妻熟人妻熟丝袜美| 国产女主播在线喷水免费视频网站| 免费久久久久久久精品成人欧美视频 | 大香蕉久久网| 免费不卡的大黄色大毛片视频在线观看| 麻豆成人午夜福利视频| 国产午夜精品一二区理论片| 2022亚洲国产成人精品| 18禁在线无遮挡免费观看视频| 免费少妇av软件| av在线观看视频网站免费| 国精品久久久久久国模美| 久久久a久久爽久久v久久| 日产精品乱码卡一卡2卡三| 亚洲,一卡二卡三卡| 夜夜爽夜夜爽视频| 亚洲av成人精品一二三区| 夜夜骑夜夜射夜夜干| 欧美变态另类bdsm刘玥| 国产精品成人在线| 亚洲高清免费不卡视频| 午夜免费观看性视频| 日本色播在线视频| 亚洲综合精品二区| 久久久欧美国产精品| av国产免费在线观看| 一级毛片我不卡| av专区在线播放| 91aial.com中文字幕在线观看| av卡一久久| 成人午夜精彩视频在线观看| 欧美精品一区二区免费开放| av卡一久久| 精品视频人人做人人爽| 国产精品成人在线| 一级毛片我不卡| 18+在线观看网站| 久久精品国产自在天天线| 高清欧美精品videossex| 赤兔流量卡办理| 国产一区二区三区av在线| 能在线免费看毛片的网站| 亚洲丝袜综合中文字幕| 欧美日韩在线观看h| 精品久久久精品久久久| 久久精品国产亚洲av涩爱| 一级毛片 在线播放| 欧美xxxx性猛交bbbb| 岛国毛片在线播放| 一边亲一边摸免费视频| 欧美日韩亚洲高清精品| 国产精品一区二区性色av| 视频中文字幕在线观看| 高清在线视频一区二区三区| 寂寞人妻少妇视频99o| 在线观看国产h片| 伊人久久国产一区二区| 男人舔奶头视频| av国产精品久久久久影院| 欧美97在线视频| 观看av在线不卡| 久久久久精品性色| 久久久久国产网址| 九色成人免费人妻av| 国产无遮挡羞羞视频在线观看| 插逼视频在线观看| 国产精品精品国产色婷婷| 久久久久久人妻| 人妻夜夜爽99麻豆av| 嫩草影院新地址| 又爽又黄a免费视频| 久久99蜜桃精品久久| 丰满少妇做爰视频| 蜜桃在线观看..| 草草在线视频免费看| 少妇人妻精品综合一区二区| 亚洲中文av在线| 在线观看三级黄色| 男人和女人高潮做爰伦理| 国产伦在线观看视频一区| 久久久久久伊人网av| 干丝袜人妻中文字幕| 我的老师免费观看完整版| 欧美zozozo另类| 99久久精品一区二区三区| 免费观看a级毛片全部| 国产片特级美女逼逼视频| 高清视频免费观看一区二区| 亚洲欧洲国产日韩| 亚洲精品成人av观看孕妇| 黑丝袜美女国产一区| 免费大片18禁| 日韩电影二区| 色吧在线观看| 亚洲av不卡在线观看| 国产女主播在线喷水免费视频网站| 日日摸夜夜添夜夜爱| 国产成人一区二区在线| 国产男女内射视频| 中文在线观看免费www的网站| 日本黄色日本黄色录像| 啦啦啦视频在线资源免费观看| 国产成人精品一,二区| 久久国产乱子免费精品| 国产亚洲精品久久久com| 亚洲精品一二三| 爱豆传媒免费全集在线观看| 亚洲精品成人av观看孕妇| 成人亚洲精品一区在线观看 | 高清日韩中文字幕在线| 亚洲国产色片| 国产在线免费精品| 啦啦啦在线观看免费高清www| 狂野欧美激情性bbbbbb| 久久久久久久久久人人人人人人| 欧美少妇被猛烈插入视频| 久久久久久久国产电影| 交换朋友夫妻互换小说| 在线亚洲精品国产二区图片欧美 | 五月天丁香电影| 黄色配什么色好看| 国产精品免费大片| 免费黄频网站在线观看国产| 亚洲国产av新网站| 狠狠精品人妻久久久久久综合| 男女啪啪激烈高潮av片| 卡戴珊不雅视频在线播放| 精品一区二区三卡| 亚洲精品乱久久久久久| 亚洲不卡免费看| 久久久亚洲精品成人影院| 国产免费一区二区三区四区乱码| 日韩在线高清观看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 啦啦啦啦在线视频资源| 日本黄大片高清| 91在线精品国自产拍蜜月| 国产午夜精品一二区理论片| 在线观看三级黄色| 成人综合一区亚洲| 又粗又硬又长又爽又黄的视频| 精品99又大又爽又粗少妇毛片| 99热国产这里只有精品6| 成人二区视频| 精华霜和精华液先用哪个| 亚洲精品第二区| 边亲边吃奶的免费视频| 一本一本综合久久| 亚洲av二区三区四区| 国产精品久久久久久av不卡| 一个人看视频在线观看www免费| 91久久精品国产一区二区三区| 久久这里有精品视频免费| 啦啦啦在线观看免费高清www| 又黄又爽又刺激的免费视频.| 国产一区二区在线观看日韩| 欧美丝袜亚洲另类| 国产女主播在线喷水免费视频网站| 国语对白做爰xxxⅹ性视频网站| 久久久久久久久久人人人人人人| 午夜免费鲁丝| 国产毛片在线视频| 欧美一区二区亚洲| 亚洲性久久影院| 男女国产视频网站| 久久久欧美国产精品| a级毛片免费高清观看在线播放| 一级毛片aaaaaa免费看小| 欧美精品国产亚洲| 免费少妇av软件| 一级毛片 在线播放| 在线观看人妻少妇| 国产精品一及| 精品亚洲成国产av| 我的女老师完整版在线观看| 尤物成人国产欧美一区二区三区| 最新中文字幕久久久久| 中文字幕制服av| 哪个播放器可以免费观看大片| 99热这里只有是精品在线观看| 日韩av不卡免费在线播放| 水蜜桃什么品种好| 国产永久视频网站| 黄色视频在线播放观看不卡| 交换朋友夫妻互换小说| 老司机影院成人| 身体一侧抽搐| 午夜福利视频精品| 国产一区二区三区综合在线观看 | 好男人视频免费观看在线| 国产免费一级a男人的天堂| 麻豆成人午夜福利视频| 三级国产精品欧美在线观看| 亚洲精品成人av观看孕妇| 国产亚洲5aaaaa淫片| 99热全是精品| 少妇人妻 视频| 亚洲第一av免费看| 久久人人爽av亚洲精品天堂 | 亚洲综合精品二区| 久久久亚洲精品成人影院| 欧美成人精品欧美一级黄| 国产精品国产三级国产av玫瑰| 亚洲成人一二三区av| 亚洲综合色惰| a级一级毛片免费在线观看| 狂野欧美激情性xxxx在线观看| 亚洲国产精品999| av卡一久久| 国产熟女欧美一区二区| 一本久久精品| 夜夜看夜夜爽夜夜摸| 毛片女人毛片| 国产黄色免费在线视频| 特大巨黑吊av在线直播| 各种免费的搞黄视频| 国产av一区二区精品久久 | 精品99又大又爽又粗少妇毛片| 七月丁香在线播放| 亚洲av中文字字幕乱码综合| 久久人人爽av亚洲精品天堂 | 国模一区二区三区四区视频| 日韩一本色道免费dvd| 妹子高潮喷水视频| 一区二区av电影网| 成年女人在线观看亚洲视频| 国产综合精华液| 夫妻性生交免费视频一级片| 99久国产av精品国产电影| 蜜臀久久99精品久久宅男| 内地一区二区视频在线| 国产一区亚洲一区在线观看| 久热这里只有精品99| 日韩成人av中文字幕在线观看|