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

    淮河上游流域SCS_CN模型初損取值與CN值確定方法的研究

    2017-03-21 02:04:16黃兆歡張銀雪張友靜
    中國農(nóng)村水利水電 2017年12期
    關(guān)鍵詞:漸近線損率徑流

    黃兆歡,劉 陽,張銀雪,曾 天,王 欣,張友靜

    (河海大學(xué)地球科學(xué)與工程學(xué)院,南京 211100)

    目前,由降雨產(chǎn)流導(dǎo)致的土壤侵蝕和水土流失已經(jīng)成為國際社會關(guān)注的熱點(diǎn)問題,如何建立準(zhǔn)確的流域水文模型來確定下墊面徑流量,對水土流失的治理、預(yù)測具有重要的意義[1,2]。美國農(nóng)業(yè)部水土保持局提出的SCS-CN經(jīng)驗(yàn)?zāi)P涂梢杂脕矸治霾煌捌谕寥罎穸取⒉煌寥李愋?、不同土地利用方式及不同坡度等因素對降雨徑流的影響情況。由于此模型結(jié)構(gòu)簡單,使用方便,現(xiàn)已被國際上廣泛運(yùn)用于不同流域地區(qū)的徑流量計(jì)算中[3]。該模型基于兩個(gè)參數(shù)計(jì)算,其中一個(gè)是徑流曲線數(shù)(CN),Ponce等[4]指明,CN值對SCS-CN模型影響重大,CN值變化±10%可引起-45%~55%的徑流變化。雖然美國農(nóng)業(yè)部提供了CN值查找表[5],但由于各地區(qū)氣候、土壤、地形等因素的差異,許多研究表明直接利用該表誤差較大。目前可利用次降雨徑流資料來反推流域的CN值,國內(nèi)外采用的相關(guān)方法主要有平均值法、漸近線法、S對數(shù)頻率分布法等。張鈺嫻等[6]利用62場降雨徑流數(shù)據(jù)對黃土丘陵區(qū)不同坡度下的CN值進(jìn)行了研究。符素華、王向亮等[3]人在北京密云石匣CN方法的研究中,認(rèn)為算數(shù)平均值法和漸近線法計(jì)算降雨徑流效果最好。

    初損率λ是SCS-CN模型的另一個(gè)重要參數(shù),它是降雨徑流過程初損量與最大潛在蓄水量的比值。在模型研究過程中,美國農(nóng)業(yè)部水土保持局提出λ值大致分布在0.095~0.38之間的概率占50%,因此取λ=0.2作為模型初損率參數(shù)的標(biāo)準(zhǔn)值被公認(rèn)推廣[7]。而大量研究表明初損率標(biāo)準(zhǔn)參數(shù)λ=0.2并不適用于其他流域下墊面的研究分析,其取值的確定應(yīng)具有地域性。穆宏強(qiáng)、張秀英等[8-9]研究表明,λ取值主要在0.00~0.30之間變化,并且不同區(qū)域變化明顯。目前通過許多相關(guān)文獻(xiàn)調(diào)研得出λ=0.05更加適合預(yù)測降雨徑流[10],特別是在降雨較少或者CN值較小的條件下[11]。

    本文以淮河上游流的大坡嶺、長臺關(guān)、息縣和淮濱四個(gè)小流域,選取2000-2010年的次降雨觀測數(shù)據(jù),采用平均值法、中值法、算數(shù)平均值法、S對數(shù)頻率分布法、最小二乘法及漸近線法六種方法推求各流域在初損率λ=0.2及λ=0.05時(shí)相應(yīng)的CN值,并結(jié)合R2、NS、RMSE、PBIAS參數(shù)評價(jià)淮河上游流域的最適λ參數(shù)取值和CN值測定方法,為SCS-CN模型在本流域的進(jìn)一步應(yīng)用提供有效參考。

    1 研究區(qū)概況

    選擇淮河上游淮濱水文站以上流域作為研究區(qū)域(圖1),流域范圍為東經(jīng)113°15′~115°27′,北緯31°28′~32°44′。流域內(nèi)地形以低平的平原和緩丘為主,小部分為平原洼地,呈現(xiàn)西高東低趨勢,由西南逐漸向東北傾斜。流域地處北亞熱帶與暖溫帶氣候過渡地帶,季風(fēng)氣候明顯,光熱資源豐實(shí)。降雨在時(shí)間和空間上分布不均勻,主要集中在夏季與研究區(qū)的西部。選取流域內(nèi)4個(gè)水文站點(diǎn)(大坡嶺、長臺關(guān)、息縣和淮濱),各對應(yīng)的流域面積為1 640、3 090、10 190、16 005 km2。

    圖1 研究區(qū)概況Fig.1 Situations of study region

    2 研究數(shù)據(jù)和方法

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

    數(shù)據(jù)來源于淮河流域水文年鑒,根據(jù)4個(gè)流域2000-2010年逐日降雨和徑流數(shù)據(jù),選取了單個(gè)次降雨總量大于25 mm的次降雨事件[10],具體統(tǒng)計(jì)結(jié)果為表1。其中需要將實(shí)測的日徑流數(shù)據(jù)通過基流分割,減去地下徑流,為SCS模型計(jì)算出的真實(shí)徑流量[12]。從上游大坡嶺至下游淮濱水文站,平均徑流系數(shù)呈遞減趨勢。

    表1 各流域次降雨徑流數(shù)據(jù)統(tǒng)計(jì)Tab.1 Statistics of rainfall runoff data in each watershed

    2.2 SCS-CN模型

    本研究采用的是20世紀(jì)50年代美國農(nóng)業(yè)部水土保持局建立SCS-CN模型,該模型的建立基于一個(gè)水平衡方程和兩個(gè)假設(shè)條件。水平衡方程為:

    P=Ia+F+Q

    (1)

    一個(gè)基本假設(shè)是實(shí)際地表徑流量Q與流域可能最大徑流量(P-Ia)的比值等于實(shí)際入滲量F與潛在蓄水能力S之比,表示公式為:

    (2)

    另一個(gè)基本假設(shè)是初損Ia和潛在蓄水能力S之間存在一定的關(guān)系,公式表示為:

    Ia=λS

    (3)

    由式(1)~(3)可得:

    (4)

    式中:Q,P,Ia,λ,S和CN分別代表徑流量,降雨量,初損,初損系數(shù),潛在蓄水能力和流域曲線參數(shù)。

    在式(4)中以λ=0.2為例,代入后經(jīng)過簡單的計(jì)算,S的值能通過公式表示出來結(jié)果為:

    S0.20=5[P+2Q-(4Q2+5PQ)0.5]

    (5)

    當(dāng)初損λ=0.05時(shí),對應(yīng)的S0.20和S0.05有如下轉(zhuǎn)換關(guān)系:

    S0.05=0.818 7S1.150.20

    (6)

    根據(jù)公式(7),得出CN值。

    (7)

    2.3 CN值確定方法

    (1) 平均值法。利用篩選后的降雨徑流數(shù)據(jù),計(jì)算潛在蓄水能力S值,如果λ=0.2,利用公式(5)來決定S值;如果λ=0.05,利用公式(5)和(6)確定S值。然后算出S的平均值,將S值的平均值代入(7)中得到CN值。

    (2)中位數(shù)法。首先利用上述方法得出潛在蓄水能力S值,將S值重新排序,選出其中的中位數(shù)值作為最終的S值,利用公式(7)得到CN值。

    (3)算術(shù)平均值法。計(jì)算出潛在蓄水能力S值,對于每一個(gè)降雨徑流數(shù)據(jù)對應(yīng)的S值首先利用式(7)計(jì)算CN,然后取CN均值為最后結(jié)果。

    (4)S對數(shù)頻率分布法。計(jì)算出每對降雨徑流數(shù)據(jù)對應(yīng)的潛在蓄水能力S值,分別對S值求其對數(shù),并得出對數(shù)頻率經(jīng)驗(yàn)分布曲線中50%對應(yīng)的S值,并算出相應(yīng)的CN。

    (5)最小二乘法。最小二乘法使用迭代的非線性最小二乘擬合方法與NRCS CN方程[10]。Q測、Q計(jì)分別為直接測量所得降雨量和經(jīng)方法計(jì)算后所得降雨量,該方法是通過代入λ后降低Q測和Q計(jì)之間的方差來求出所求CN的值,所采用公式如下:

    (8)

    (6)漸近線法。漸近線法具體可參考Richard H. Hawkins[13]。首先將所選地區(qū)的P,Q值按照大小排序,排序結(jié)果按順序依次對應(yīng)。然后利用式(7)得到每對CN值,并對所求得CN值與對應(yīng)的P值配對繪制出P-CN的散點(diǎn)圖 ,以其漸近線值作為最后的CN值。

    2.4 模型評價(jià)參數(shù)

    研究中分別用了6種方法計(jì)算CN值,考慮到多種因素對結(jié)果的影響,采用四種模型評估參數(shù)R2、NS、RMSE和PBIAS,通過綜合對比不同評估參數(shù)獲得最佳結(jié)果。

    (1)R2模型:

    (2)NS模型:

    NS的取值范圍是-∞~1.0,值越大代表模型適用性越好。模型的適用性可參考兩個(gè)標(biāo)準(zhǔn),分別是2007年Moriasi[14]提出的NS>0.5時(shí)模型可接受和2013年Ritter和Munoz-Carpena[15]提出的水文模型標(biāo)準(zhǔn)NS>0.65時(shí)模型可接受。

    (3)RMSE模型:

    RMSE是用于描述測量值與估算值之間的誤差的評價(jià)參數(shù)模型,其計(jì)算值越小代表相應(yīng)的模型適用性越好,最佳取值為RMSE=0,這是理想情況下的取值。

    (4)PBIAS模型:

    PBIAS是用來量化實(shí)測值與估算值之間誤差百分比和顯示過高或過低估算趨勢的水文模型。理想的模型取值為PBIAS=0,一般不會取到該值,因此PBIAS取值越趨近于0,偏離程度越小,模型效果越好。

    3 結(jié)果與分析

    表2給出了在不同初損率(λ=0.2和λ=0.05)和6種CN值推算方法下得出的4個(gè)流域的CN值。初損率與降雨產(chǎn)生的徑流呈反比,CN值與徑流呈正比。因此,在λ=0.05條件下得出的CN值小于λ=0.2時(shí)得出的CN值。六種CN值推算方法得出的結(jié)果存在一定差異,其中S對數(shù)頻率分布法對應(yīng)的CN值最大。λ=0.2時(shí),最小二乘法對應(yīng)的CN值最小,λ=0.05時(shí),漸近線法對應(yīng)的CN值最小。表3給出了不同參數(shù)下計(jì)算的次降雨徑流均值。平均值法、中值法、算數(shù)平均值法和S對數(shù)頻率分布法在λ為0.2時(shí)模擬的平均徑流都大于λ為0.05的模擬結(jié)果,最小二乘法和漸進(jìn)線法則相反。平均值法最接近實(shí)測值,漸近線法次之。λ為0.05時(shí)模擬的平均徑流相比λ為0.2時(shí)模擬的平均徑流更靠近實(shí)測值。

    表2 不同初損和不同方法下的CN值Tab.2 CN values under different λ and different methods

    表3 不同初損不同方法下的徑流均值Tab.3 the Mean of runoff under different λ and different methods

    采用4種評價(jià)指標(biāo)(R2,NSE,RMSE和PBIAS)對初損率為0.05時(shí),4個(gè)流域在6種方法下次降雨徑流模擬結(jié)果評價(jià)(表4)。R2對6種方法評價(jià)的差異性較小,而其他3種評價(jià)指標(biāo)可反映出不同方法的估算效果。雖然6種方法在4個(gè)流域的評價(jià)效果略有差異,但總體呈現(xiàn)為:在NSE和RMSE評價(jià)指標(biāo)中,漸近線法和最小二乘法對4個(gè)流域的徑流值計(jì)算更具適用性;在PBIAS評價(jià)指標(biāo)中,平均值法和漸近線法模擬徑流的效果最佳。因此綜合上述4種評價(jià)指標(biāo),可得出:漸近線法最好,平均值法次之。

    表5給出了初損率為0.2時(shí),不同方法的次降雨徑流模擬結(jié)果評價(jià)指標(biāo)。6種方法在4個(gè)流域數(shù)據(jù)中的總體估算效果,基本呈現(xiàn):漸近線法最好,平均值法次之。這與初損率為0.05時(shí)的趨勢一致。但在對比這兩種不同的初損率的條件下,初損率為0.05在NSE和RMSE指標(biāo)下,其次降雨徑流模擬結(jié)果明顯優(yōu)于初損率為0.2的結(jié)果。因此,將原SCS模型中的初損率設(shè)為0.05,可進(jìn)一步提高徑流模擬效果。Ajmal M[16]也曾在南韓的10個(gè)山地流域研究中指出利用初損λ=0.05和漸近線法預(yù)測降雨徑流的精度會更好。

    4 結(jié) 語

    水文模型參數(shù)的準(zhǔn)確識別獲取對提高水文模型應(yīng)用效果具有意義,本文根據(jù)淮河上游4個(gè)流域的次降雨數(shù)據(jù),針對SCS

    表4 初損率為0.05下不同方法的模型評價(jià)參數(shù)的比較Tab.4 Comparison of model evaluation parametersunder different methods with λ=0.05

    表5 初損率為0.2下不同方法的模型評價(jià)參數(shù)的比較Tab.5 Comparison of model evaluation parametersunder different methods with λ=0.2

    模型中的主要參數(shù)CN值和初損率λ進(jìn)行分析,得到的結(jié)論如下:

    (1)將λ=0.2改為λ=0.05,可提高降雨徑流模擬精度。

    (2)不同評價(jià)指標(biāo)(R2、NSE、RMSE和PBIAS)評價(jià)結(jié)果略有差異,需綜合考慮。

    (3)6種CN值推求方法(平均值法、中值法、算數(shù)平均值法、S對數(shù)頻率分布法、最小二乘法和漸近線法)中,漸近線法最好,平均值法次之。

    [1] 李智廣,曹 煒,劉秉正,等. 我國水土流失狀況與發(fā)展趨勢研究[J]. 中國水土保持科學(xué),2008,(1):57-62.

    [2] Nearing M A, Liu B Y, Risse L M, et al. Curve numbers and Green-Ampt effective hydraulic conductivities[J]. Jawra Journal of the American Water Resources Association, 1996,32(1):125-136.

    [3] 符素華,王向亮. SCS-CN徑流模型中CN值確定方法研究[J]. 干旱區(qū)地理,2012,(3):415-421.

    [4] Victor M. Ponce. Runoff Curve Number: Has It Reached Maturity?[J]. Journal of Hydrologic Engineering,1996,1(1):11-19.

    [5] 魏文秋,謝淑琴. 遙感資料在SCS模型產(chǎn)流計(jì)算中的應(yīng)用[J]. 環(huán)境遙感,1992,(4):243-250.

    [6] Allen T. Hjelmfelt. Investigation of Curve Number Procedure[J]. Journal of Hydraulic Engineering,1991,117(6):725-737.

    [7] 張鈺嫻,穆興民,王飛. 徑流曲線數(shù)模型(SCS-CN)參數(shù)λ在黃土丘陵區(qū)的率定[J]. 干旱地區(qū)農(nóng)業(yè)研究,2008,(5):124-128.

    [8] 穆宏強(qiáng). SCS模型在石橋鋪流域的應(yīng)用研究[J]. 水利學(xué)報(bào),1992,(10):79-83,89.

    [9] 張秀英,孟飛,丁寧. SCS模型在干旱半干旱區(qū)小流域徑流估算中的應(yīng)用[J]. 水土保持研究,2003,(4):172-174,249.

    [10] Hawkins R H, Ward T, Woodward D E, et al. Curve Number Hydrology: State of the Practice[J]. Journal of Hydrologic Engineering,2009,14(9):978-0-7844-7257-6.

    [11] Samuel J. Lamont,Robert N. Eli,Jerald J. Fletcher. Continuous Hydrologic Models and Curve Numbers: A Path Forward[J]. Journal of Hydrologic Engineering,2008,13(7):621-635.

    [12] 雷曉輝. 基于二元演化模式的流域水文模型初探[D]. 北京:中國水利水電科學(xué)研究院,2001.

    [13] Richard H. Hawkins. Asymptotic Determination of Runoff Curve Numbers from Data[J]. Journal of Irrigation and Drainage Engineering,1993,119(2):334-345.

    [14] DN Moriasi,JL Fouss,RL Bengtson.Modeling the effects of deep chiseling with DRAINMOD for alluvial soils[J]. Transactions of the Asabe, 2007,50 (2):543-556.

    [15] Munoz-Carpena R,Regalado C M, Ritter A.TDR estimation of electrical conductivity and saline solute concentration in a volcanic soil[J].Geoderma, 2005,124(3/4):399-413.

    [16] Ajmal M, Kim T W, Ahn J H. Stability assessment of the curve number methodology used to estimate excess rainfall in forest-dominated watersheds[J]. Arabian Journal of Geosciences, 2016,9(5):1-14.

    猜你喜歡
    漸近線損率徑流
    我國水庫淤損情勢分析
    關(guān)于Pα漸近線
    無功補(bǔ)償極限線損率分析及降損措施探究
    電子制作(2017年2期)2017-05-17 03:55:17
    供電企業(yè)月度實(shí)際線損率定量計(jì)算方法
    電子制作(2016年1期)2016-11-07 08:42:53
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    漸近線,你值得擁有
    探秘“大徑流”
    攻克“大徑流”
    漸近線問題研究
    江埡水庫降雨徑流相關(guān)圖的建立
    亚洲av中文av极速乱| 我的女老师完整版在线观看| 精品福利观看| 18+在线观看网站| 欧美一级a爱片免费观看看| 欧美成人免费av一区二区三区| 国产成人影院久久av| 久久久精品94久久精品| 国产 一区精品| 中文字幕av在线有码专区| 丝袜喷水一区| 日韩亚洲欧美综合| 综合色丁香网| 天堂√8在线中文| 国产成人91sexporn| 亚洲乱码一区二区免费版| 日本黄大片高清| 欧美日韩国产亚洲二区| 特大巨黑吊av在线直播| 在线免费观看不下载黄p国产| 99久国产av精品国产电影| av在线播放精品| 国产精品人妻久久久影院| 18禁黄网站禁片免费观看直播| 亚洲专区国产一区二区| 搡老熟女国产l中国老女人| 联通29元200g的流量卡| 99久国产av精品| 成年免费大片在线观看| 在现免费观看毛片| 日韩成人av中文字幕在线观看 | 一个人免费在线观看电影| 亚洲国产欧洲综合997久久,| 日韩精品有码人妻一区| 两个人的视频大全免费| 亚洲专区国产一区二区| 午夜视频国产福利| 一本精品99久久精品77| 国产又黄又爽又无遮挡在线| 国产精品伦人一区二区| 国产亚洲欧美98| 亚洲人与动物交配视频| 18禁在线无遮挡免费观看视频 | 国产色婷婷99| 直男gayav资源| 国产美女午夜福利| 免费观看的影片在线观看| 嫩草影院入口| 丝袜美腿在线中文| 大又大粗又爽又黄少妇毛片口| 国国产精品蜜臀av免费| 国内揄拍国产精品人妻在线| 蜜臀久久99精品久久宅男| 亚洲精品日韩在线中文字幕 | 国产精品不卡视频一区二区| 九九爱精品视频在线观看| 成人二区视频| 欧美成人精品欧美一级黄| 亚洲av电影不卡..在线观看| 精品无人区乱码1区二区| 啦啦啦啦在线视频资源| 热99re8久久精品国产| 麻豆国产av国片精品| 最近在线观看免费完整版| 男人的好看免费观看在线视频| 亚洲七黄色美女视频| 国产视频一区二区在线看| 亚洲综合色惰| 国产成人a区在线观看| 久久久久久久久久久丰满| 最近视频中文字幕2019在线8| 搡女人真爽免费视频火全软件 | 久久久久久久午夜电影| 成人特级av手机在线观看| a级毛片a级免费在线| 九九热线精品视视频播放| 黑人高潮一二区| 国产成年人精品一区二区| 亚洲av成人av| 精品久久国产蜜桃| 久久久久久国产a免费观看| 床上黄色一级片| 亚洲av成人av| 亚洲内射少妇av| 精品久久国产蜜桃| 一级毛片电影观看 | 一级黄片播放器| 午夜激情欧美在线| 亚洲人与动物交配视频| 热99re8久久精品国产| 最新在线观看一区二区三区| 美女高潮的动态| 亚洲美女搞黄在线观看 | 精品人妻视频免费看| 美女免费视频网站| 一a级毛片在线观看| 亚洲美女视频黄频| 国产精品,欧美在线| 波多野结衣高清作品| 床上黄色一级片| 国产精品爽爽va在线观看网站| av在线播放精品| 日韩av在线大香蕉| 99久久无色码亚洲精品果冻| 丝袜美腿在线中文| 国产精品久久视频播放| 日本欧美国产在线视频| 99在线视频只有这里精品首页| 麻豆国产av国片精品| 在线免费观看不下载黄p国产| 一进一出抽搐动态| 久久精品夜色国产| 在线看三级毛片| 1024手机看黄色片| 午夜激情福利司机影院| 99热这里只有是精品在线观看| 美女高潮的动态| 国产精品久久久久久久久免| 亚洲av不卡在线观看| 国产av一区在线观看免费| 真实男女啪啪啪动态图| 联通29元200g的流量卡| 成人欧美大片| 日日摸夜夜添夜夜添小说| 亚洲中文日韩欧美视频| 真人做人爱边吃奶动态| 久久亚洲国产成人精品v| 久久久欧美国产精品| 国产午夜精品论理片| 搞女人的毛片| 欧美日韩国产亚洲二区| 看十八女毛片水多多多| 熟女电影av网| 1000部很黄的大片| 国产免费一级a男人的天堂| 欧美一区二区精品小视频在线| 亚洲激情五月婷婷啪啪| 日韩精品青青久久久久久| 自拍偷自拍亚洲精品老妇| 亚洲乱码一区二区免费版| 18禁黄网站禁片免费观看直播| 亚洲成人久久性| av.在线天堂| av天堂在线播放| 日韩三级伦理在线观看| 哪里可以看免费的av片| 日韩欧美精品v在线| 国产一区二区三区av在线 | 我的老师免费观看完整版| 国产一区二区在线观看日韩| 精品国产三级普通话版| 我的女老师完整版在线观看| 一级毛片电影观看 | 国产精品一二三区在线看| 欧美另类亚洲清纯唯美| 欧美一级a爱片免费观看看| 国产午夜精品论理片| 国产精品1区2区在线观看.| 日本-黄色视频高清免费观看| 神马国产精品三级电影在线观看| 草草在线视频免费看| 精品午夜福利在线看| 亚洲最大成人av| 欧美zozozo另类| 国产片特级美女逼逼视频| 久久精品国产99精品国产亚洲性色| av黄色大香蕉| 亚洲成av人片在线播放无| 国产毛片a区久久久久| 在线观看免费视频日本深夜| 国产午夜福利久久久久久| 婷婷精品国产亚洲av在线| 久久久久久伊人网av| 欧美一区二区国产精品久久精品| 超碰av人人做人人爽久久| 少妇猛男粗大的猛烈进出视频 | 国产三级中文精品| 美女免费视频网站| 波多野结衣高清作品| 一区二区三区四区激情视频 | 一个人观看的视频www高清免费观看| 日韩av在线大香蕉| 一本一本综合久久| 韩国av在线不卡| 精品人妻一区二区三区麻豆 | 三级国产精品欧美在线观看| 亚洲成人中文字幕在线播放| 长腿黑丝高跟| 亚洲av熟女| 国产美女午夜福利| 亚洲美女搞黄在线观看 | 五月玫瑰六月丁香| 成人午夜高清在线视频| av女优亚洲男人天堂| 精品午夜福利视频在线观看一区| 国产精品野战在线观看| 国产成人freesex在线 | 欧美成人a在线观看| 欧美色视频一区免费| 在线看三级毛片| 国产亚洲欧美98| 国产精品电影一区二区三区| 黄色日韩在线| 亚洲av中文av极速乱| 亚洲精品色激情综合| 免费人成视频x8x8入口观看| 国产人妻一区二区三区在| 成熟少妇高潮喷水视频| 国产精品一区二区免费欧美| 国产av麻豆久久久久久久| 三级经典国产精品| ponron亚洲| 亚洲精品色激情综合| 偷拍熟女少妇极品色| 成人av在线播放网站| 十八禁国产超污无遮挡网站| 国产精品国产高清国产av| 此物有八面人人有两片| 免费不卡的大黄色大毛片视频在线观看 | 色5月婷婷丁香| 欧美最新免费一区二区三区| 日本与韩国留学比较| 天堂动漫精品| 一级毛片久久久久久久久女| 久久九九热精品免费| 又爽又黄无遮挡网站| 欧美潮喷喷水| 日本熟妇午夜| ponron亚洲| 久久精品影院6| 一本久久中文字幕| 欧美色欧美亚洲另类二区| 中文在线观看免费www的网站| 伊人久久精品亚洲午夜| 久久精品人妻少妇| 亚洲av.av天堂| 女人十人毛片免费观看3o分钟| 麻豆精品久久久久久蜜桃| 蜜桃久久精品国产亚洲av| 亚洲欧美日韩东京热| 亚洲精品日韩在线中文字幕 | 中文字幕av成人在线电影| 草草在线视频免费看| 狠狠狠狠99中文字幕| 中文字幕精品亚洲无线码一区| 国产高清三级在线| 亚洲国产精品合色在线| 午夜福利在线在线| 乱人视频在线观看| 男人和女人高潮做爰伦理| 别揉我奶头 嗯啊视频| 婷婷亚洲欧美| www日本黄色视频网| 国产一区二区在线av高清观看| 高清毛片免费观看视频网站| 色视频www国产| 日日摸夜夜添夜夜添av毛片| 悠悠久久av| 99热网站在线观看| 一区二区三区四区激情视频 | 色视频www国产| 欧美一级a爱片免费观看看| 国产午夜精品论理片| 嫩草影视91久久| 女生性感内裤真人,穿戴方法视频| 午夜免费激情av| 欧美3d第一页| 亚洲av中文字字幕乱码综合| 悠悠久久av| 在线免费观看的www视频| 69av精品久久久久久| 桃色一区二区三区在线观看| 精品久久久久久久久久免费视频| 成人永久免费在线观看视频| 久久热精品热| 日韩欧美三级三区| 97热精品久久久久久| 99久国产av精品| 国产精品乱码一区二三区的特点| 国产av在哪里看| 精品久久国产蜜桃| 99在线视频只有这里精品首页| 精品久久久久久久久久免费视频| 精品不卡国产一区二区三区| 国内揄拍国产精品人妻在线| 深爱激情五月婷婷| 一级黄片播放器| 国产伦精品一区二区三区四那| 在线天堂最新版资源| 五月伊人婷婷丁香| 校园人妻丝袜中文字幕| 免费在线观看成人毛片| 日日摸夜夜添夜夜添小说| 12—13女人毛片做爰片一| 久久人人爽人人爽人人片va| 午夜福利视频1000在线观看| 国内久久婷婷六月综合欲色啪| 伦精品一区二区三区| 中文字幕av成人在线电影| 有码 亚洲区| 国内精品宾馆在线| 久久久久久久久久久丰满| 欧美不卡视频在线免费观看| 无遮挡黄片免费观看| 露出奶头的视频| 最近的中文字幕免费完整| 成人美女网站在线观看视频| 波野结衣二区三区在线| 久久精品久久久久久噜噜老黄 | 精品久久久噜噜| 黄色一级大片看看| 美女大奶头视频| 亚洲精品国产成人久久av| 国产真实乱freesex| 日本一本二区三区精品| 三级男女做爰猛烈吃奶摸视频| 18禁黄网站禁片免费观看直播| 偷拍熟女少妇极品色| 在线看三级毛片| 赤兔流量卡办理| 国产女主播在线喷水免费视频网站 | 天天躁日日操中文字幕| 高清毛片免费观看视频网站| 久久久久久伊人网av| 日日摸夜夜添夜夜爱| 99久国产av精品| 久久鲁丝午夜福利片| 少妇丰满av| 99国产精品一区二区蜜桃av| 97超级碰碰碰精品色视频在线观看| 亚洲精华国产精华液的使用体验 | 在线天堂最新版资源| 97超视频在线观看视频| 午夜精品在线福利| 国产久久久一区二区三区| 国内久久婷婷六月综合欲色啪| 日韩国内少妇激情av| 3wmmmm亚洲av在线观看| 日本爱情动作片www.在线观看 | 免费看a级黄色片| 欧美日韩乱码在线| 又爽又黄无遮挡网站| 不卡视频在线观看欧美| 亚洲精品成人久久久久久| 丝袜喷水一区| h日本视频在线播放| 亚洲精品色激情综合| 美女高潮的动态| 又黄又爽又免费观看的视频| 精品一区二区免费观看| 最近在线观看免费完整版| 亚洲精华国产精华液的使用体验 | 热99在线观看视频| 禁无遮挡网站| 尾随美女入室| 久久久久久久久大av| 国产美女午夜福利| 成人高潮视频无遮挡免费网站| 午夜激情福利司机影院| 国产黄片美女视频| 亚洲自偷自拍三级| 国产69精品久久久久777片| 黄色日韩在线| 国产精品一及| 久久精品国产亚洲网站| 国产精品一区二区性色av| 日本爱情动作片www.在线观看 | 国产精品久久久久久av不卡| 国产精品一区二区三区四区免费观看 | 狂野欧美白嫩少妇大欣赏| 国产午夜福利久久久久久| 男人的好看免费观看在线视频| 亚洲成人av在线免费| 亚洲高清免费不卡视频| 欧美区成人在线视频| 成人无遮挡网站| 亚洲专区国产一区二区| av卡一久久| 黄色一级大片看看| 久久天躁狠狠躁夜夜2o2o| 网址你懂的国产日韩在线| 日日干狠狠操夜夜爽| 日韩成人伦理影院| 欧美一区二区亚洲| 欧美绝顶高潮抽搐喷水| 校园人妻丝袜中文字幕| 国产精品乱码一区二三区的特点| 免费搜索国产男女视频| 成人亚洲欧美一区二区av| 欧美成人a在线观看| 日韩在线高清观看一区二区三区| 欧美一区二区亚洲| 国产精品国产高清国产av| 嫩草影院新地址| 成人永久免费在线观看视频| 成人午夜高清在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 偷拍熟女少妇极品色| 亚洲国产日韩欧美精品在线观看| а√天堂www在线а√下载| 国产成人freesex在线 | 寂寞人妻少妇视频99o| 99久久精品国产国产毛片| 亚洲中文字幕日韩| 久久午夜亚洲精品久久| 97碰自拍视频| 亚洲人与动物交配视频| 人人妻,人人澡人人爽秒播| 国产精品一二三区在线看| 午夜影院日韩av| 国产亚洲欧美98| 久久人人精品亚洲av| 欧美中文日本在线观看视频| 婷婷六月久久综合丁香| 国产在视频线在精品| 精品一区二区免费观看| 男女那种视频在线观看| 色综合站精品国产| 午夜日韩欧美国产| 国产精品女同一区二区软件| 人妻制服诱惑在线中文字幕| 99久久无色码亚洲精品果冻| 精品福利观看| 日本一二三区视频观看| 国产高清视频在线播放一区| 日本黄色片子视频| 18禁黄网站禁片免费观看直播| 精品午夜福利视频在线观看一区| 91麻豆精品激情在线观看国产| 级片在线观看| 亚洲自偷自拍三级| 麻豆一二三区av精品| 最近视频中文字幕2019在线8| 欧美日本视频| 精品免费久久久久久久清纯| 国产白丝娇喘喷水9色精品| 一边摸一边抽搐一进一小说| 亚洲中文字幕日韩| 亚洲国产色片| 亚洲自拍偷在线| 极品教师在线视频| 男女做爰动态图高潮gif福利片| av黄色大香蕉| 亚洲av成人av| 精品久久久久久久人妻蜜臀av| 免费看美女性在线毛片视频| 人妻制服诱惑在线中文字幕| 欧美成人a在线观看| 亚洲人成网站在线播放欧美日韩| 可以在线观看的亚洲视频| 国产精品久久久久久亚洲av鲁大| 女人被狂操c到高潮| avwww免费| 亚洲国产精品国产精品| av免费在线看不卡| 男女那种视频在线观看| 亚洲最大成人av| 日韩国内少妇激情av| 免费人成视频x8x8入口观看| 欧美三级亚洲精品| 午夜激情欧美在线| 国产三级在线视频| 亚洲av电影不卡..在线观看| 五月玫瑰六月丁香| a级一级毛片免费在线观看| 舔av片在线| 伦理电影大哥的女人| 变态另类丝袜制服| 亚洲三级黄色毛片| 国产成年人精品一区二区| av专区在线播放| 亚洲最大成人av| 伊人久久精品亚洲午夜| 美女cb高潮喷水在线观看| 国产高潮美女av| 简卡轻食公司| 一级a爱片免费观看的视频| 亚洲美女搞黄在线观看 | 大又大粗又爽又黄少妇毛片口| 综合色av麻豆| 搞女人的毛片| 插阴视频在线观看视频| 国产精品久久久久久亚洲av鲁大| 免费看日本二区| 色综合亚洲欧美另类图片| 国产又黄又爽又无遮挡在线| 变态另类丝袜制服| 成人永久免费在线观看视频| 高清毛片免费观看视频网站| 久久久久久九九精品二区国产| 99热这里只有是精品在线观看| 色噜噜av男人的天堂激情| 亚洲在线自拍视频| 麻豆国产97在线/欧美| 国产精品野战在线观看| 精品免费久久久久久久清纯| 精品一区二区三区视频在线观看免费| 三级毛片av免费| 国产一区二区三区av在线 | 国产极品精品免费视频能看的| 别揉我奶头~嗯~啊~动态视频| 日本-黄色视频高清免费观看| 亚洲最大成人手机在线| 欧美+亚洲+日韩+国产| 中文字幕免费在线视频6| 在线免费观看的www视频| 日日摸夜夜添夜夜添小说| 亚洲av第一区精品v没综合| 日本成人三级电影网站| 日韩欧美精品免费久久| 男女啪啪激烈高潮av片| 亚洲不卡免费看| 日韩欧美一区二区三区在线观看| 国产色婷婷99| 亚洲国产精品合色在线| 日本欧美国产在线视频| 久久人人精品亚洲av| 别揉我奶头 嗯啊视频| 亚洲熟妇熟女久久| 国产毛片a区久久久久| 久久国内精品自在自线图片| 国产三级在线视频| 十八禁国产超污无遮挡网站| 亚洲最大成人av| 亚洲av五月六月丁香网| 简卡轻食公司| 18禁黄网站禁片免费观看直播| 精品99又大又爽又粗少妇毛片| 国产精品一区二区三区四区久久| 搡老熟女国产l中国老女人| av卡一久久| 欧美一区二区精品小视频在线| 亚洲欧美精品综合久久99| 精品乱码久久久久久99久播| 国产精品人妻久久久久久| 午夜福利在线观看吧| 午夜福利高清视频| 亚洲av不卡在线观看| 波多野结衣高清作品| 国产午夜精品论理片| 日韩高清综合在线| 婷婷精品国产亚洲av在线| 99久久成人亚洲精品观看| 国产一区二区亚洲精品在线观看| 日韩,欧美,国产一区二区三区 | 一边摸一边抽搐一进一小说| 中文字幕av成人在线电影| 国产v大片淫在线免费观看| 长腿黑丝高跟| 亚洲最大成人中文| 亚洲av免费高清在线观看| 18禁黄网站禁片免费观看直播| 亚洲在线观看片| 超碰av人人做人人爽久久| 永久网站在线| 女人被狂操c到高潮| 亚州av有码| 少妇的逼水好多| 国产成人a∨麻豆精品| 永久网站在线| 校园春色视频在线观看| 国产黄a三级三级三级人| 日韩,欧美,国产一区二区三区 | 最近手机中文字幕大全| 国产男靠女视频免费网站| 国产aⅴ精品一区二区三区波| 我要搜黄色片| 熟女电影av网| 亚洲成a人片在线一区二区| 在线免费观看的www视频| 久久久久免费精品人妻一区二区| 在线观看午夜福利视频| 午夜影院日韩av| 国产 一区 欧美 日韩| 最近视频中文字幕2019在线8| av.在线天堂| 啦啦啦啦在线视频资源| 亚洲18禁久久av| 欧美不卡视频在线免费观看| 久99久视频精品免费| 草草在线视频免费看| 中国美白少妇内射xxxbb| 亚洲va在线va天堂va国产| 日本五十路高清| 欧美一级a爱片免费观看看| 精品免费久久久久久久清纯| 免费在线观看影片大全网站| 久久精品国产亚洲av天美| 91久久精品国产一区二区三区| 国产高清不卡午夜福利| 成人二区视频| 国产亚洲91精品色在线| 婷婷亚洲欧美| 美女 人体艺术 gogo| 日韩欧美一区二区三区在线观看| 亚洲,欧美,日韩| 久久国内精品自在自线图片| 亚洲精品日韩av片在线观看| 麻豆一二三区av精品| 最后的刺客免费高清国语| 亚洲成人精品中文字幕电影| 欧美潮喷喷水| 99riav亚洲国产免费| 欧美日韩精品成人综合77777| 成人精品一区二区免费| 99久久精品国产国产毛片| 午夜a级毛片| 日日摸夜夜添夜夜添av毛片| 日本撒尿小便嘘嘘汇集6| 国产精品野战在线观看| 久久精品国产清高在天天线| 老司机福利观看| 久久久午夜欧美精品| 久久精品夜色国产| av中文乱码字幕在线|