• <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)圖的建立
    国产免费又黄又爽又色| 亚洲在久久综合| 欧美一区二区亚洲| 老司机影院毛片| 免费看av在线观看网站| 日本猛色少妇xxxxx猛交久久| 亚洲中文字幕日韩| 内射极品少妇av片p| 国产亚洲精品av在线| 久久这里有精品视频免费| 亚洲久久久久久中文字幕| 亚洲av电影不卡..在线观看| 国产美女午夜福利| 免费观看的影片在线观看| 亚洲av不卡在线观看| 国产成人午夜福利电影在线观看| 中文字幕av成人在线电影| 久久99热这里只有精品18| 一夜夜www| 精品久久久久久久末码| 99热6这里只有精品| 久久久久久久亚洲中文字幕| 寂寞人妻少妇视频99o| 国产精品av视频在线免费观看| 亚洲av电影在线观看一区二区三区 | 亚洲av免费高清在线观看| 99热这里只有精品一区| 国语自产精品视频在线第100页| 精品人妻视频免费看| 欧美性猛交黑人性爽| 精品欧美国产一区二区三| 亚洲国产精品国产精品| 99视频精品全部免费 在线| 看非洲黑人一级黄片| 久久综合国产亚洲精品| 美女高潮的动态| 国产精品不卡视频一区二区| 国产午夜精品久久久久久一区二区三区| 九九热线精品视视频播放| 国产精品一区www在线观看| 国产免费又黄又爽又色| 国产爱豆传媒在线观看| 乱人视频在线观看| 少妇人妻精品综合一区二区| 欧美激情国产日韩精品一区| 黄片wwwwww| 国产黄色小视频在线观看| 麻豆av噜噜一区二区三区| 久久久久免费精品人妻一区二区| 亚洲成人中文字幕在线播放| 亚洲国产最新在线播放| 免费观看精品视频网站| 国产极品天堂在线| 嫩草影院新地址| 干丝袜人妻中文字幕| 色5月婷婷丁香| 亚洲真实伦在线观看| 亚洲精品色激情综合| 国产精华一区二区三区| 99久久精品一区二区三区| 成人毛片a级毛片在线播放| 床上黄色一级片| 99在线视频只有这里精品首页| 天堂网av新在线| 午夜福利在线观看免费完整高清在| 美女黄网站色视频| 成人av在线播放网站| 国产欧美日韩精品一区二区| 男女啪啪激烈高潮av片| 我要看日韩黄色一级片| 热99re8久久精品国产| 久久精品熟女亚洲av麻豆精品 | 青春草视频在线免费观看| 成人午夜精彩视频在线观看| 久久久精品欧美日韩精品| 九九在线视频观看精品| 毛片女人毛片| 国产精品精品国产色婷婷| 精品一区二区三区人妻视频| videos熟女内射| 欧美激情在线99| 亚洲成人中文字幕在线播放| 国产伦理片在线播放av一区| 日韩一区二区三区影片| 亚洲av电影在线观看一区二区三区 | 亚洲性久久影院| 午夜爱爱视频在线播放| 99在线视频只有这里精品首页| 久久久久九九精品影院| 99久久精品热视频| 国产精品国产高清国产av| 亚洲最大成人手机在线| 亚洲精品国产av成人精品| 2021天堂中文幕一二区在线观| 中文资源天堂在线| 午夜爱爱视频在线播放| 国产精品日韩av在线免费观看| 天堂中文最新版在线下载 | 国产片特级美女逼逼视频| 亚洲国产精品专区欧美| av在线亚洲专区| 国产色爽女视频免费观看| 日本欧美国产在线视频| 亚洲av成人av| 欧美一区二区精品小视频在线| 青春草视频在线免费观看| 亚洲国产色片| 我要看日韩黄色一级片| 国产午夜精品久久久久久一区二区三区| 日韩一区二区视频免费看| 亚洲欧美一区二区三区国产| 亚洲av电影在线观看一区二区三区 | 亚洲欧美日韩东京热| 亚洲人成网站在线播| 99九九线精品视频在线观看视频| 黄片无遮挡物在线观看| 亚洲,欧美,日韩| 老司机福利观看| 18禁在线无遮挡免费观看视频| 99久久中文字幕三级久久日本| 身体一侧抽搐| 国产成人午夜福利电影在线观看| 欧美丝袜亚洲另类| 国产精品久久久久久精品电影小说 | 久久国内精品自在自线图片| 精品久久久久久久久av| 波多野结衣高清无吗| 欧美激情国产日韩精品一区| 别揉我奶头 嗯啊视频| 国产精华一区二区三区| 丰满乱子伦码专区| 日产精品乱码卡一卡2卡三| 国产成人freesex在线| 又爽又黄a免费视频| 精品一区二区三区视频在线| 观看美女的网站| 菩萨蛮人人尽说江南好唐韦庄 | 能在线免费观看的黄片| 91在线精品国自产拍蜜月| 91av网一区二区| 日本免费一区二区三区高清不卡| 成人午夜精彩视频在线观看| 又黄又爽又刺激的免费视频.| 能在线免费观看的黄片| 久久久久九九精品影院| 综合色av麻豆| 亚洲欧洲日产国产| 亚洲自偷自拍三级| 日韩大片免费观看网站 | 久久久精品94久久精品| 亚洲av福利一区| 久久久久久久久大av| 久久精品久久久久久噜噜老黄 | 插逼视频在线观看| 精品一区二区三区人妻视频| 亚洲精品久久久久久婷婷小说 | 2021天堂中文幕一二区在线观| 91久久精品电影网| 亚洲乱码一区二区免费版| 精品人妻偷拍中文字幕| 成人毛片60女人毛片免费| 日本免费在线观看一区| 男人舔奶头视频| 中文天堂在线官网| 国产乱人视频| 两个人的视频大全免费| 亚洲中文字幕一区二区三区有码在线看| 纵有疾风起免费观看全集完整版 | 成人二区视频| 久久久久性生活片| 99热全是精品| 国产人妻一区二区三区在| 成人鲁丝片一二三区免费| 国产精品乱码一区二三区的特点| 精品少妇黑人巨大在线播放 | 亚洲一级一片aⅴ在线观看| 男人舔女人下体高潮全视频| 日本黄色片子视频| 国产精品一二三区在线看| 国产乱来视频区| 亚洲欧美一区二区三区国产| 国产成人aa在线观看| 亚洲av不卡在线观看| 如何舔出高潮| 国产在线一区二区三区精 | 日韩成人伦理影院| 在线观看66精品国产| av在线观看视频网站免费| 综合色av麻豆| 能在线免费看毛片的网站| 久久99热6这里只有精品| 国产成人午夜福利电影在线观看| 国产精品久久久久久久电影| 亚洲av男天堂| 一边亲一边摸免费视频| 国产高清有码在线观看视频| 简卡轻食公司| 黄色配什么色好看| 日韩一本色道免费dvd| 黄片wwwwww| 天美传媒精品一区二区| 最后的刺客免费高清国语| 日韩 亚洲 欧美在线| av线在线观看网站| 老司机影院毛片| 国产三级在线视频| 国产欧美另类精品又又久久亚洲欧美| 丰满少妇做爰视频| 超碰97精品在线观看| 国产免费又黄又爽又色| 伦精品一区二区三区| 免费av不卡在线播放| 日本wwww免费看| 亚洲国产欧美人成| 国模一区二区三区四区视频| 内地一区二区视频在线| 国产又色又爽无遮挡免| 亚洲电影在线观看av| 91av网一区二区| 一二三四中文在线观看免费高清| 国产精品久久久久久精品电影小说 | 欧美性猛交黑人性爽| 18禁在线无遮挡免费观看视频| 国产中年淑女户外野战色| 成人欧美大片| 色网站视频免费| 亚洲精品久久久久久婷婷小说 | 在现免费观看毛片| 欧美成人免费av一区二区三区| 一个人看视频在线观看www免费| 久久欧美精品欧美久久欧美| 久久精品影院6| 大香蕉97超碰在线| 亚洲国产成人一精品久久久| 精品一区二区三区人妻视频| 亚洲av成人精品一二三区| 亚洲精品一区蜜桃| kizo精华| 婷婷色麻豆天堂久久 | 久久久久久久久大av| 日韩成人av中文字幕在线观看| 最近的中文字幕免费完整| 国产乱人偷精品视频| 尾随美女入室| 全区人妻精品视频| 亚洲精品成人久久久久久| 狂野欧美激情性xxxx在线观看| 国产乱来视频区| 成人一区二区视频在线观看| 欧美区成人在线视频| 亚洲av熟女| 日韩一区二区三区影片| 欧美潮喷喷水| 久久欧美精品欧美久久欧美| 在线免费观看的www视频| av在线播放精品| 欧美色视频一区免费| 成人午夜精彩视频在线观看| 国产午夜福利久久久久久| 哪个播放器可以免费观看大片| 国产黄片视频在线免费观看| 亚洲人与动物交配视频| 精品久久久久久久久久久久久| 亚洲国产高清在线一区二区三| 五月伊人婷婷丁香| 国产视频首页在线观看| 毛片女人毛片| 国产精品电影一区二区三区| 人妻夜夜爽99麻豆av| 国国产精品蜜臀av免费| 亚洲精品国产av成人精品| 色尼玛亚洲综合影院| 天堂中文最新版在线下载 | 国模一区二区三区四区视频| 一级爰片在线观看| 在线观看66精品国产| 黄片wwwwww| 日本av手机在线免费观看| 蜜桃久久精品国产亚洲av| 免费一级毛片在线播放高清视频| 免费观看在线日韩| 国国产精品蜜臀av免费| 亚洲最大成人手机在线| 直男gayav资源| 久久久久久久久中文| 欧美日韩在线观看h| 淫秽高清视频在线观看| 久99久视频精品免费| 青春草视频在线免费观看| 免费观看在线日韩| 九九在线视频观看精品| ponron亚洲| 国产欧美另类精品又又久久亚洲欧美| 51国产日韩欧美| 国产激情偷乱视频一区二区| 亚洲中文字幕一区二区三区有码在线看| 国产一区二区亚洲精品在线观看| 大香蕉97超碰在线| 十八禁国产超污无遮挡网站| 人人妻人人澡人人爽人人夜夜 | 91狼人影院| 亚洲欧洲国产日韩| av视频在线观看入口| 亚洲av不卡在线观看| 国产精品国产三级国产av玫瑰| 女人十人毛片免费观看3o分钟| 干丝袜人妻中文字幕| 亚洲成人av在线免费| 午夜a级毛片| 国产免费视频播放在线视频 | 噜噜噜噜噜久久久久久91| 欧美成人免费av一区二区三区| 一个人看视频在线观看www免费| 自拍偷自拍亚洲精品老妇| 亚洲人与动物交配视频| 美女被艹到高潮喷水动态| 成人毛片60女人毛片免费| 99在线人妻在线中文字幕| 女人被狂操c到高潮| 亚洲人成网站在线观看播放| 亚洲av日韩在线播放| 最近最新中文字幕免费大全7| 中文字幕精品亚洲无线码一区| 麻豆一二三区av精品| 日韩精品青青久久久久久| 日本猛色少妇xxxxx猛交久久| 成人鲁丝片一二三区免费| 亚洲成av人片在线播放无| 国产极品天堂在线| 欧美成人a在线观看| 久久久久久久久久久丰满| 国产黄片美女视频| 国产亚洲91精品色在线| 亚洲av免费高清在线观看| 午夜精品一区二区三区免费看| 中文精品一卡2卡3卡4更新| 亚洲国产欧美人成| 日韩 亚洲 欧美在线| 联通29元200g的流量卡| 日本午夜av视频| 嫩草影院新地址| 少妇熟女欧美另类| 久久久a久久爽久久v久久| 亚洲av中文av极速乱| 国产精品日韩av在线免费观看| 国产精品久久视频播放| 秋霞在线观看毛片| 国产伦一二天堂av在线观看| 波野结衣二区三区在线| 色综合亚洲欧美另类图片| 69人妻影院| 丝袜喷水一区| 寂寞人妻少妇视频99o| 亚洲国产最新在线播放| 尤物成人国产欧美一区二区三区| 大香蕉久久网| 91在线精品国自产拍蜜月| 欧美成人精品欧美一级黄| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品人妻久久久影院| 久久久色成人| 国产在线男女| 久久这里只有精品中国| 99久久精品国产国产毛片| 黄色欧美视频在线观看| 3wmmmm亚洲av在线观看| 一二三四中文在线观看免费高清| 99热这里只有精品一区| 边亲边吃奶的免费视频| 色综合站精品国产| 少妇熟女欧美另类| 尾随美女入室| av又黄又爽大尺度在线免费看 | 建设人人有责人人尽责人人享有的 | 亚洲av.av天堂| 日本欧美国产在线视频| 亚洲国产精品合色在线| 色噜噜av男人的天堂激情| 午夜老司机福利剧场| 亚洲,欧美,日韩| 日韩精品有码人妻一区| 禁无遮挡网站| 国产探花在线观看一区二区| 欧美又色又爽又黄视频| 嫩草影院新地址| 亚洲欧美精品综合久久99| 亚洲精品aⅴ在线观看| 亚洲最大成人av| 国产一级毛片在线| 九草在线视频观看| АⅤ资源中文在线天堂| 非洲黑人性xxxx精品又粗又长| 午夜精品一区二区三区免费看| .国产精品久久| 国产伦精品一区二区三区视频9| 人妻系列 视频| 久久久久久九九精品二区国产| 亚洲av免费高清在线观看| 日韩欧美在线乱码| 男人和女人高潮做爰伦理| 国产色婷婷99| 国产亚洲5aaaaa淫片| 99国产精品一区二区蜜桃av| 成人性生交大片免费视频hd| 欧美不卡视频在线免费观看| 国产成人91sexporn| 亚洲va在线va天堂va国产| 亚洲在线观看片| 中文欧美无线码| 亚洲高清免费不卡视频| av在线老鸭窝| 亚洲精品成人久久久久久| 晚上一个人看的免费电影| 最近中文字幕2019免费版| 亚洲18禁久久av| 久久精品国产亚洲av天美| 国产精品一区二区性色av| 好男人在线观看高清免费视频| 精品久久久久久久人妻蜜臀av| videossex国产| 亚洲av中文字字幕乱码综合| 我的女老师完整版在线观看| 午夜爱爱视频在线播放| 一个人看视频在线观看www免费| 亚洲av中文av极速乱| 久久久色成人| 麻豆成人av视频| 美女黄网站色视频| 内地一区二区视频在线| 亚洲av日韩在线播放| 成人漫画全彩无遮挡| 中文字幕制服av| 村上凉子中文字幕在线| 欧美成人一区二区免费高清观看| 最近最新中文字幕大全电影3| 两个人视频免费观看高清| 日本-黄色视频高清免费观看| 91久久精品国产一区二区三区| 波野结衣二区三区在线| 色综合站精品国产| 91精品国产九色| 国产成人精品久久久久久| 嫩草影院入口| 国产淫语在线视频| 少妇被粗大猛烈的视频| 午夜免费激情av| 久久人妻av系列| 久久欧美精品欧美久久欧美| 久久久久久大精品| 内地一区二区视频在线| 九九在线视频观看精品| 欧美性猛交黑人性爽| 久久精品熟女亚洲av麻豆精品 | 热99re8久久精品国产| 成人国产麻豆网| 欧美成人a在线观看| 国产淫语在线视频| 免费不卡的大黄色大毛片视频在线观看 | 中文字幕人妻熟人妻熟丝袜美| 男的添女的下面高潮视频| 亚洲av熟女| 国产日韩欧美在线精品| 97人妻精品一区二区三区麻豆| 国产精品美女特级片免费视频播放器| 男女下面进入的视频免费午夜| 亚洲国产色片| 亚洲第一区二区三区不卡| 在线观看美女被高潮喷水网站| 又爽又黄a免费视频| 亚洲欧美精品自产自拍| 亚洲va在线va天堂va国产| 日韩成人伦理影院| 亚州av有码| 久久人妻av系列| 亚洲成人精品中文字幕电影| 亚洲四区av| 夜夜看夜夜爽夜夜摸| 久久久精品欧美日韩精品| 精品人妻一区二区三区麻豆| 婷婷色麻豆天堂久久 | 国产色爽女视频免费观看| 久久久久久久久大av| 非洲黑人性xxxx精品又粗又长| 国产黄片美女视频| 亚洲国产精品成人综合色| 国产精品一区二区性色av| 亚洲av熟女| 国产激情偷乱视频一区二区| 午夜a级毛片| 久久久久网色| 我的老师免费观看完整版| 91久久精品国产一区二区成人| 欧美最新免费一区二区三区| 亚洲av电影在线观看一区二区三区 | 国产探花极品一区二区| 亚洲精品乱久久久久久| 国产单亲对白刺激| 国产精品无大码| 我的老师免费观看完整版| 听说在线观看完整版免费高清| 天堂影院成人在线观看| 亚洲伊人久久精品综合 | 精品人妻一区二区三区麻豆| 久久99蜜桃精品久久| 韩国高清视频一区二区三区| 国产熟女欧美一区二区| 一夜夜www| 啦啦啦观看免费观看视频高清| 成人三级黄色视频| 99久久精品国产国产毛片| 1024手机看黄色片| 超碰97精品在线观看| 亚洲伊人久久精品综合 | 亚洲欧美中文字幕日韩二区| 久久热精品热| 久久久国产成人免费| 亚洲成人av在线免费| 综合色av麻豆| av黄色大香蕉| 在线免费观看的www视频| 欧美成人免费av一区二区三区| 亚洲精品自拍成人| 色播亚洲综合网| 久久久欧美国产精品| 99久久中文字幕三级久久日本| 干丝袜人妻中文字幕| av国产免费在线观看| 国产黄色视频一区二区在线观看 | 久久久久久久久大av| 观看免费一级毛片| 欧美bdsm另类| 国产成人aa在线观看| 国产av码专区亚洲av| 免费观看性生交大片5| 国产探花极品一区二区| 亚洲av福利一区| 特大巨黑吊av在线直播| 天堂中文最新版在线下载 | 久久精品国产亚洲网站| 高清午夜精品一区二区三区| 欧美一区二区亚洲| 中文字幕久久专区| 亚洲欧洲日产国产| 国产精品麻豆人妻色哟哟久久 | 亚洲精品乱码久久久久久按摩| 美女内射精品一级片tv| 久久亚洲国产成人精品v| 一个人免费在线观看电影| 国产午夜精品一二区理论片| 久久人人爽人人爽人人片va| 赤兔流量卡办理| 亚洲熟妇中文字幕五十中出| 日韩欧美国产在线观看| 国产高清视频在线观看网站| 国产久久久一区二区三区| 亚洲精品日韩av片在线观看| 如何舔出高潮| 亚洲精品自拍成人| 亚洲四区av| 国产av码专区亚洲av| 久久久国产成人免费| 日本-黄色视频高清免费观看| 两个人视频免费观看高清| 国产精品1区2区在线观看.| 日本爱情动作片www.在线观看| 日日摸夜夜添夜夜添av毛片| 一区二区三区四区激情视频| 日韩欧美在线乱码| 久久久久久国产a免费观看| 亚洲欧美成人精品一区二区| 国产成人福利小说| 美女被艹到高潮喷水动态| 99九九线精品视频在线观看视频| 日韩视频在线欧美| 亚洲精品色激情综合| 麻豆一二三区av精品| 黑人高潮一二区| 亚洲天堂国产精品一区在线| 国产男人的电影天堂91| 丝袜喷水一区| 国产精品一及| av在线老鸭窝| 亚洲av.av天堂| 日韩欧美国产在线观看| 亚洲精品日韩av片在线观看| 亚洲国产精品国产精品| 少妇人妻一区二区三区视频| 水蜜桃什么品种好| 日本黄色视频三级网站网址| 嫩草影院入口| 亚洲av日韩在线播放| 国语对白做爰xxxⅹ性视频网站| 国内精品美女久久久久久| 能在线免费看毛片的网站| 欧美激情久久久久久爽电影| 青青草视频在线视频观看| 大香蕉久久网| 国产黄片视频在线免费观看| 波多野结衣巨乳人妻| 嫩草影院新地址| 亚洲无线观看免费| 99国产精品一区二区蜜桃av| 免费看光身美女| 国产精品爽爽va在线观看网站| 搞女人的毛片| 日本av手机在线免费观看| 高清午夜精品一区二区三区| 非洲黑人性xxxx精品又粗又长| 最新中文字幕久久久久| 夜夜爽夜夜爽视频| 在线观看美女被高潮喷水网站| 精品国产露脸久久av麻豆 | 嫩草影院精品99| 亚洲经典国产精华液单| 久久久国产成人免费|