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

    模擬降雨下植被蓋度對坡面流水動力學(xué)特性的影響

    2015-03-13 03:28:45孫佳美余新曉樊登星梁洪儒李瀚之
    生態(tài)學(xué)報 2015年8期
    關(guān)鍵詞:流態(tài)蓋度黑麥草

    孫佳美, 余新曉, 樊登星, 梁洪儒, 常 玉, 李瀚之

    北京林業(yè)大學(xué), 水土保持與荒漠化防治教育部重點實驗室, 北京 100083

    模擬降雨下植被蓋度對坡面流水動力學(xué)特性的影響

    孫佳美, 余新曉*, 樊登星, 梁洪儒, 常 玉, 李瀚之

    北京林業(yè)大學(xué), 水土保持與荒漠化防治教育部重點實驗室, 北京 100083

    通過模擬降雨實驗的方法,分析研究了坡度10°和20°,降雨強度30mm/h和60mm/h條件下不同蓋度黑麥草對坡面產(chǎn)流產(chǎn)沙的調(diào)控過程,并從雷諾數(shù)、弗勞德數(shù)和阻力系數(shù)三個方面對水流運動過程和黑麥草調(diào)控坡面流的水力學(xué)特性進行了剖析。研究結(jié)果表明:雷諾數(shù)隨坡度增加而相對增大,隨降雨強度增大有明顯增大趨勢。黑麥草覆蓋能夠明顯減小坡面徑流雷諾數(shù),在各降雨強度和坡度條件下,雷諾數(shù)隨黑麥草蓋度增加而減小,雷諾數(shù)大小一般呈現(xiàn):裸坡>20%>40%>60%>80%。黑麥草蓋度對坡面流弗勞德數(shù)有顯著影響,隨著黑麥草蓋度增加弗勞德數(shù)減小,并且弗勞德數(shù)隨蓋度變化為:裸地>20%>40%>60%>80%,坡面阻力系數(shù)與坡面產(chǎn)沙率有良好的擬合關(guān)系,隨坡面阻力系數(shù)的增大,坡面產(chǎn)沙率呈對數(shù)減小,并且阻力系數(shù)在0—1時減小速率很大,阻力系數(shù)大于1以后減小曲線較為平緩。

    黑麥草; 坡面流; 模擬降雨; 水力學(xué); 流態(tài)

    坡面徑流是造成水土流失的主導(dǎo)因子,在坡面布設(shè)合理的水保措施對防治土壤侵蝕,保持水土,改善生態(tài)環(huán)境十分重要。目前對于坡面流流態(tài)的研究很多,包括對水流參數(shù)的計算和坡面產(chǎn)流產(chǎn)沙機制的研究[1]。國內(nèi)外眾多學(xué)者致力于坡面流特性研究,國外部分學(xué)者從土壤入滲、坡面層流特征[2]、斜坡坡面流水深和流速測量和預(yù)報研究[3]、坡面流細溝侵蝕的水力要素[4]、土壤顆粒分離過程[5-6]、泥沙輸移啟動及運移過程、坡面流動能和功率等揭示坡面流流態(tài)機理的研究做了大量工作[7]。一般用雷諾數(shù)和弗勞德數(shù)來表征坡面流態(tài)的特性,坡面流流態(tài)主要有混合流、擾動流、過度流、擾動層流、偽層流等多種流態(tài)[8],Horton認為坡面流流態(tài)是在湍流中間偶有層流的混合流態(tài)[9];Emmett研究表明坡面流流態(tài)為擾動流,及同時具有紊流和層流的特性[10]。

    國內(nèi)學(xué)者利用室內(nèi)人工模擬降雨實驗[11]和變坡水槽放水實驗[12-13]對坡面流流態(tài)及其水動力學(xué)參數(shù)做了很多研究[14-15],對坡面流水力學(xué)特性進行了侵蝕機理研究并取得了很多研究成果。吳普特等認為降雨造成的坡面流流態(tài)是擾動層流,這是由于雨滴的擾動引起的[16];姚文藝等經(jīng)過研究認為降雨時的坡面流不是簡單的層流,而是同時存在著層流、紊流和過度流3種流態(tài)水流[17];張光輝等認為坡面流流態(tài)主要為過度流和紊流[3]。眾多研究均表明坡面薄層水流水力學(xué)特性是一個及其復(fù)雜的過程,當(dāng)?shù)乇碛胁荼局参锔采w時問題會更加復(fù)雜[18],草本植物對于坡面流有明顯的阻礙作用[19],可以減小坡面流流速[20]。

    學(xué)者們對于坡面土壤侵蝕過程的研究主要集中在黃土高原[12],但是褐土作為北京市主要土壤類型,2005年調(diào)查數(shù)據(jù)顯示,北京市水土流失面積有54621 hm2,分布在昌平,延慶,房山等區(qū)縣,因此對于褐土坡面土壤侵蝕的研究十分重要[19]。本文通過模擬降雨實驗研究褐土坡面不同降雨強度,不同坡度和不同蓋度黑麥草(LoliumperenneL.)草地坡面流形態(tài)變化規(guī)律及特征,對揭示坡面侵蝕機理非常有意義。

    1 試驗設(shè)計與方法

    1.1 試驗設(shè)計

    模擬降雨實驗在北京林業(yè)大學(xué)降雨大廳進行,降雨大廳位于北京市鷲峰實驗林場,主要由降水系統(tǒng)、坡面系統(tǒng)、控制系統(tǒng)構(gòu)成,共4個8m×8m的降雨區(qū),實際降雨控制總面積256m2。劃分為4個獨立降雨區(qū),可以分別降雨,也可以同時降雨,計算機可進行不同降雨狀態(tài)的自動調(diào)控,有效降雨高度12m,降雨過程計算機全自動控制,支持手動操作臺、計算機監(jiān)控、計算機全自動模擬降雨控制。采用旋轉(zhuǎn)下噴式噴頭,疊加式噴頭模擬自然降雨。降雨分布均勻度達到85%以上,坡面系統(tǒng)為可自由調(diào)節(jié)坡度的土槽,土槽規(guī)格為2m×0.5m×0.3m,槽內(nèi)盛放褐土,取自北京市昌平區(qū)。土槽尾部設(shè)置V形收集口,用來收集形成的徑流和泥沙,底部打孔,使土壤水可以滲透出土槽外。其他3邊額外加10cm高度,防止雨滴使侵蝕的泥沙濺出土槽之外,造成泥沙量的損失。實驗坡度為10°和20°,降雨強度選定30mm/h和60mm/h。

    黑麥草種植控制蓋度分別為0, 20%, 40%, 60%, 80%和100%,采用條帶狀種植的方法達到控制效果,首先通過種業(yè)公司調(diào)查黑麥草種植的最優(yōu)密度,將按照最優(yōu)密度種植作為100%的蓋度,然后按照控制蓋度比例確定單位面積需要的黑麥草種子質(zhì)量,將這些種子條帶狀均勻種植,每天保證其日曬時間確??梢哉9夂献饔蒙L,根據(jù)其實際生長狀況并咨詢草業(yè)專家認為一般種植2個月后即可用于實驗,在進行降水實驗時要進一步確定其真實的蓋度,對草進行修剪使其最大化符合實驗控制的蓋度。

    1.2 試驗方法

    按照設(shè)計容重分層回填土,設(shè)計實驗容重為1.34g/cm3,每10cm為一層,下面兩層分層填天然細沙,然后鋪一層紗布,上面裝實驗用土。裝土?xí)r要注意土面保持平整,每層裝完要刻意制造一定的地表糙度,防止土的滑移;為保證實驗初始含水量相同,對坡面進行預(yù)降雨,待土壤達到飽和剛開始產(chǎn)流停止降雨,用塑料布蓋上靜置24h用于降雨實驗[21];降雨過程中,記錄產(chǎn)流時間,產(chǎn)流開始后前10min每2min接1組徑流,以后每5min 1組徑流,稱量徑流體積進行分析。降雨過程中同時采用染色劑法測定坡面流速,由于坡面流在剛產(chǎn)生初期不易穩(wěn)定,因此流速測量斷面選擇在坡下部,具體位置從坡底部向上30—130cm,并沿坡面流水平方向均勻布設(shè)5個點,沿坡面流流動方向測定5個流速,求其平均值,并乘以修正系數(shù)0.75,作為水流斷面平均流速。

    1.3 分析方法

    坡面流極薄,實測方法測定的水流深誤差較大,并不能代表真實狀況,因此采用反推的方法求得坡面薄層水流水深。因此在假定坡面薄層水流是均勻分布的條件下,采用以下公式計算坡面薄層水流水深[7]:

    h=q/U=Q/(uBt)

    (1)

    式中,h為薄層水流水深(m);q為單寬流量(m3min-1m-1);Q為t時間內(nèi)徑流量(m3);t徑流取樣間隔時間(min);U為斷面平均流速(m/min);B為過水?dāng)嗝鎸挾?m)。

    雷諾數(shù)Re在流體力學(xué)中是流體慣性力與黏性力比值的量度,它是一個無量綱量,通過雷諾數(shù)可以判定水流流態(tài)。根據(jù)明渠均勻流的水力學(xué)理論,當(dāng)雷諾數(shù)>500時,坡面流流態(tài)為紊流;當(dāng)雷諾數(shù)<500時,坡面流流態(tài)為層流;當(dāng)雷諾數(shù)在500左右時,流態(tài)為過度流。采用以下公式計算雷諾數(shù)Re[22]:

    Re=Uh/v

    (2)

    式中,v為水運動黏性系數(shù),取值為7.0×10-5m2/s。

    弗勞德數(shù)Fr是流體力學(xué)中表征流體慣性力和重力相對大小的一個無量綱參數(shù),它表示慣性力和重力量級的比。當(dāng)Fr>1時,慣性力對水流起主導(dǎo)作用,水流為急流;Fr< 1時,重力起主導(dǎo)作用,水流為緩流;Fr=1時,重力、慣性力作用相等,水流為臨界流。采用以下公式計算弗勞德數(shù)Fr[7]:

    (3)

    式中,g為重力加速度(m/s2)。

    坡面流阻力是分析坡面流水力學(xué)特性的重要因子,主要來自四個方面:即顆粒阻力fg、形態(tài)阻力ff、波阻力fw和降雨阻力fr,這4個阻力的加和即為坡面流阻力f,計算公式為[7]:

    f=fg+ff+fw+fr

    (4)

    式中,fg=3.19Re-0.45;ff=4.8∑Ai/Ab,Ai為局部地形改變后的面積(m2),Ab為過水?dāng)嗝娴目偯娣e(m2);fw=2.8C,C為地表覆被蓋度(%)。

    2 結(jié)果與分析

    2.1 黑麥草覆蓋調(diào)控坡面流雷諾數(shù)

    各蓋度在不同坡度和降雨強度條件下雷諾數(shù)隨降雨歷時變化過程如圖1所示。從圖中可以看出,在降雨強度為30mm/h,坡度為10°和20°,各個蓋度覆被條件下的雷諾數(shù)均呈現(xiàn)在前30min穩(wěn)定增長,然后趨于穩(wěn)定狀態(tài)階段,雷諾數(shù)始終小于500,流態(tài)始終為層流。而當(dāng)降雨強度為60mm/h,坡度為10°和20°,各個蓋度覆被條件下的雷諾數(shù)隨降雨歷時增大,但是變化不明顯,整體處于穩(wěn)定狀態(tài),流態(tài)也始終為層流,說明降雨強度相對與坡度和蓋度對坡面流流態(tài)影響更大,并且在降雨剛開始時,坡面較為平整,坡面流流速不穩(wěn)定,雷諾數(shù)較小,隨著降雨的進行,坡面阻力增大,雷諾數(shù)也增大到穩(wěn)定狀態(tài),坡面流更加不穩(wěn)定。當(dāng)降雨強度較大時,坡面流會很快趨于穩(wěn)定,整個降雨過程雷諾數(shù)變化不大,組間雷諾數(shù)差異性顯著(P=1.01×10-7<0.05)。

    雷諾數(shù)隨坡度增加而相對增大,隨降雨強度增大有明顯增大趨勢。在各降雨強度和坡度條件下,雷諾數(shù)隨黑麥草蓋度增加而減小,雷諾數(shù)大小一般呈現(xiàn):裸坡>20%>40%>60%>80%。并且在降雨強度為30mm/h時,雷諾數(shù)隨蓋度變化不明顯,在降雨強度為60mm/h時,蓋度從0增大到40%雷諾數(shù)變化不大,當(dāng)黑麥草蓋度為60%和80%雷諾數(shù)減小非常明顯,說明當(dāng)蓋度在60%以上時對坡面流流態(tài)有較大干擾,此時黑麥草影響坡面產(chǎn)流過程和產(chǎn)沙過程較明顯,能夠顯著減少徑流和泥沙的產(chǎn)生,減少土壤侵蝕。

    圖1 不同黑麥草蓋度下徑流雷諾數(shù)變化過程Fig.1 Changeable process of Reynolds number in different ryegrass coverage

    圖2 不同草本蓋度下徑流弗勞德數(shù)變化過程Fig.2 Changeable process of Froude number in different herb cover

    圖3 坡面流阻力與坡面產(chǎn)沙率關(guān)系Fig.3 Relationship of flow resistance and sediment yield rate

    2.2 黑麥草覆蓋調(diào)控坡面流弗勞德數(shù)

    各蓋度在不同坡度和降雨強度條件下弗勞德數(shù)隨降雨歷時變化過程如圖2所示。由圖2可以看出,降雨強度、坡度和黑麥草的蓋度均對坡面流弗勞德數(shù)有較大影響。與雷諾數(shù)變化相似,降雨強度為30mm/h時,弗勞德數(shù)隨著降雨進行表現(xiàn)為開始明顯增大,然后逐漸達到穩(wěn)定狀態(tài),當(dāng)降雨強度為60mm/h,弗勞德數(shù)隨降雨歷時變化不顯著,并且弗勞德數(shù)均小于1,流態(tài)為緩流。黑麥草蓋度對坡面流弗勞德數(shù)有顯著影響,隨著黑麥草蓋度增加弗勞德數(shù)減小,在坡度20°時,降雨強度30mm/h和60mm/h均呈現(xiàn)弗勞德數(shù)隨蓋度變化為:裸地>20%>40%>60%>80%,當(dāng)坡面有黑麥草覆蓋時增加了坡面流阻力,延緩了流速,造成弗勞德數(shù)減小。但是在坡度為10°時這種隨黑麥草蓋度增加弗勞德數(shù)減小的差異不明顯,這個現(xiàn)象主要是由于坡度較小時,坡面流速相對較小,組間弗勞德數(shù)差異性顯著(P=3.59×10-9<0.05)。

    2.3 黑麥草覆蓋調(diào)控坡面流阻力

    坡面流阻力能夠影響產(chǎn)流匯流的方式,影響土壤侵蝕過程,因此坡面流阻力的研究對于表征坡面產(chǎn)流和產(chǎn)沙過程有重要意義,影響坡面流阻力的因子有很多,包括降雨強度、地表特征、坡度等。本文采用阻力疊加的方法計算坡面流阻力,一般坡面流阻力由4部分組成,顆粒阻力、形態(tài)阻力、波阻力和降雨阻力,根據(jù)力學(xué)基本原理將4個阻力疊加即為坡面流的阻力。綜合4個阻力計算得到次降雨坡面流阻力大小,并和坡面產(chǎn)沙率進行擬合見下圖3。計算結(jié)果顯示坡面阻力受波阻力影響最大,波阻力是坡面流受地表障礙物影響產(chǎn)生的阻力,對于黑麥草覆蓋實驗,波阻力的大小與黑麥草的蓋度成正比,因此坡面流阻力受蓋度影響很大。地表為裸坡時,阻力系數(shù)在0.22-0.48之間,黑麥草蓋度為20%時坡面流阻力系數(shù)在0.78—1.18之間,蓋度為40%時坡面流阻力系數(shù)在1.34—1.68之間,蓋度為60%時坡面流阻力系數(shù)在1.94—2.21之間,蓋度為80%時坡面流阻力系數(shù)在2.50—2.97之間。在次降雨中坡面流阻力直接影響坡面產(chǎn)沙量,坡面流阻力越大,坡面產(chǎn)沙越少,且阻力系數(shù)在0—1時產(chǎn)沙率隨阻力系數(shù)變化劇烈,阻力系數(shù)大于1以后變化幅度減小。經(jīng)過對次降雨阻力系數(shù)和產(chǎn)沙率進行擬合發(fā)現(xiàn)存在良好的擬合關(guān)系,擬合方程如下:

    S= -2.502ln (f)+3.5672R2=0.8532

    (5)

    式中,S為坡面產(chǎn)沙率,f為坡面阻力系數(shù)。

    3 討論

    對坡面流的水力學(xué)特性研究是研究土壤侵蝕機理的重要內(nèi)容,經(jīng)過實驗研究證明降雨強度、坡度和黑麥草的蓋度對坡面流雷諾數(shù)、弗勞德數(shù)和阻力系數(shù)有顯著影響。本研究表明,褐土坡面隨著蓋度的增加Re和Fr有顯著的減小,與吳淑芳等[7]和肖培青等[20]進行的苜蓿草地對黃土坡面流水力學(xué)特性實驗結(jié)論相同[4,13],Re是坡面流流體慣性力與黏性力的比值,Re減小則流體黏性力占主導(dǎo)地位,坡面流擾動因素會因黏性力的阻礙而減弱,流動穩(wěn)定,對地表土壤顆粒的擾動也比較弱,帶來的土壤侵蝕相應(yīng)的比較輕微;Fr是慣性力和重力的比值,隨蓋度增加Fr減小,則重力其主導(dǎo)作用,流速相對穩(wěn)定,對地表產(chǎn)生擾動較少,土壤侵蝕程度降低。

    文中采用雷諾數(shù)和弗勞德數(shù)來表征坡面流流態(tài)特征,但是坡面流水力學(xué)特征極為復(fù)雜,在實際研究當(dāng)中只用明渠水力學(xué)的方法有很多局限性,尤其在分析過程中有難度,敬向鋒等利用二元流雷諾數(shù)和擾流雷諾數(shù)以及其與阻力系數(shù)的關(guān)系判斷坡面流流態(tài)更為準確[23]。同時模擬實驗受多種因素的影響,填土的壓實均勻性、草被措施分布的狀況以及測量誤差都會造成實驗結(jié)果的不準確,而各水動力學(xué)參數(shù)受沖刷過程中坡面條件以及各個外界條件之間的相互關(guān)聯(lián)和影響都會造成實驗結(jié)果的誤差,因此模擬實驗不僅要研究單一因素對坡面流的影響,要著手因素間互相干擾效應(yīng)。在實際生產(chǎn)當(dāng)中,種植草被是常見的水土保持措施之一,實驗研究表明,草本蓋度越大坡面流速越小,流態(tài)越趨于層流對土壤表明侵蝕力度越小,但是當(dāng)蓋度增大到一定程度以后其減流效益不再增加,因此也就沒有繼續(xù)增加蓋度的意義,本實驗中蓋度60%以上其水土保持效果達到穩(wěn)定。

    4 結(jié)論

    (1)雷諾數(shù)和弗勞德數(shù)隨降雨歷時呈現(xiàn)先急劇增長,然后隨著坡面入滲和產(chǎn)流的穩(wěn)定達到穩(wěn)定的雷諾數(shù)和弗勞德數(shù),降雨強度較大時增長的時間較短,坡面流更快達到穩(wěn)定狀態(tài)。

    (2)雷諾數(shù)隨坡度增加而相對增大,隨降雨強度增大有明顯增大趨勢。黑麥草覆蓋能夠明顯減小坡面雷諾數(shù),在各降雨強度和坡度條件下,雷諾數(shù)隨黑麥草蓋度增加而減小,雷諾數(shù)大小一般呈現(xiàn):裸坡>20%>40%>60%>80%,當(dāng)蓋度大于60%以后減弱幅度很小。

    (3)黑麥草蓋度對坡面流弗勞德數(shù)有顯著影響,隨著黑麥草蓋度增加弗勞德數(shù)減小,并且弗勞德數(shù)隨蓋度變化為:裸地>20%>40%>60%>80%,

    (4)黑麥草覆蓋條件下,隨坡面阻力系數(shù)的增大,坡面產(chǎn)沙率減小,并且阻力系數(shù)在0—1時減小速率很大,阻力系數(shù)大于1以后減小曲線較為平緩。坡面阻力系數(shù)與坡面產(chǎn)沙率有良好的擬合關(guān)系,隨阻力系數(shù)增大產(chǎn)沙率呈對數(shù)減小,如公式(5)所示。

    [1] 張光輝, 衛(wèi)海燕, 劉寶元. 坡面流水動力學(xué)特性研究. 水土保持學(xué)報, 2001, 15(1): 58-61.

    [2] 田凱, 李小青, 申震洲, 蔣云鐘, 魯帆. 不同床沙下的坡面流水力學(xué)特性試驗研究. 水土保持學(xué)報, 2010, 24(2): 70-73.

    [3] 張光輝. 坡面薄層流水動力學(xué)特性的實驗研究. 水科學(xué)進展, 2002, 13(2): 159-165.

    [4] 肖培青, 鄭粉莉, 姚文藝. 坡溝系統(tǒng)坡面徑流流態(tài)及水力學(xué)參數(shù)特征研究. 水科學(xué)進展, 2009, 20(2): 236-240.

    [5] 張光輝, 劉寶元, 張科利. 坡面徑流分離土壤的水動力學(xué)實驗研究. 土壤學(xué)報, 2002, 39(6): 882-886.

    [6] 張光輝. 國外坡面徑流分離土壤過程水動力學(xué)研究進展. 水土保持學(xué)報, 2000, 14(3): 112-115.

    [7] 吳淑芳, 吳普特, 原立峰. 坡面徑流調(diào)控薄層水流水力學(xué)特性試驗. 農(nóng)業(yè)工程學(xué)報, 2010, 26(3): 14-19.

    [8] 張永東, 吳淑芳, 馮浩, 原立峰. 土壤侵蝕過程中坡面流水力學(xué)特性及侵蝕動力研究評述. 土壤, 2013, 45(1): 26-33.

    [9] Horton R E, Leach H R, Vliet V R. Laminar sheet-flow. Transactions of the American Geophysical Union, 1934, 15: 393-404.

    [10] Emmett W W. The hydraulics of overland flow on hillslopes. U S Geological Professional Paper, 1970, 662-A, A-1-A-68.

    [11] 徐向舟, 劉大慶, 張紅武, 董占地, 朱明東. 室內(nèi)人工模擬降雨試驗研究. 北京林業(yè)大學(xué)學(xué)報, 2006, 28(5): 52-58.

    [12] 潘成忠, 上官周平. 不同坡度草地含沙水流水力學(xué)特性及其攔沙機理. 水科學(xué)進展, 2007, 18(4): 490-495.

    [13] 潘成忠, 上官周平. 降雨和坡度對坡面流水動力學(xué)參數(shù)的影響. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報, 2009, 17(6): 843-851.

    [14] 汪邦穩(wěn), 肖勝生, 張光輝, 楊潔, 張利超. 南方紅壤區(qū)不同利用土地產(chǎn)流產(chǎn)沙特征試驗研究. 農(nóng)業(yè)工程學(xué)報, 2012, 28(2): 239-243.

    [15] 魏霞, 李勛貴, 李占斌, 沈冰. 黃土高原坡溝系統(tǒng)徑流水動力學(xué)特性試驗. 農(nóng)業(yè)工程學(xué)報, 2009, 25(10): 19-24.

    [16] 吳淑芳, 吳普特, 宋維秀, 卜崇峰. 坡面調(diào)控措施下的水沙輸出過程及減流減沙效應(yīng)研究. 水利學(xué)報, 2010, 41(7): 870-875.

    [17] 李勉, 姚文藝, 楊劍鋒, 陳江南, 丁文峰, 李莉, 楊春霞. 草被覆蓋對坡面流流態(tài)影響的人工模擬試驗研究. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報, 2009, 17(4): 513-523.

    [18] 曹穎, 張光輝, 唐科明, 羅榕婷. 地表模擬覆蓋率對坡面流阻力的影響. 水土保持學(xué)報, 2010, 24(4): 86-89.

    [19] 李毅, 邵明安. 草地覆蓋坡面流水動力參數(shù)的室內(nèi)降雨試驗. 農(nóng)業(yè)工程學(xué)報, 2008, 24(10): 1-5.

    [20] 肖培青, 姚文藝, 申震洲, 楊春霞. 草被減流減沙效應(yīng)及其力學(xué)機制分析. 中國水土保持科學(xué), 2010, 8(2): 15-19.

    [21] 劉振波, 史學(xué)正, 于東升, 王洪杰, 張向炎. 模擬降雨下土壤前期含水量對土壤可蝕性的影響. 生態(tài)環(huán)境, 2008, 17(1): 397-402.

    [22] 張寬地, 王光謙, 王占禮, 劉俊娥, 呂宏興. 人工加糙床面薄層滾波流水力學(xué)特性試驗. 農(nóng)業(yè)工程學(xué)報, 2011, 27(4): 28-34.

    [23] 敬向鋒, 呂宏興, 潘成忠, 潘志寶, 雒天峰, 吉麗娜. 坡面薄層水流流態(tài)判定方法的初步探討. 農(nóng)業(yè)工程學(xué)報, 2007, 23(5): 56-61.

    Impact of vegetation cover on surface runoff hydraulic characteristics with simulated rainfall

    SUN Jiamei, YU Xinxiao*, FAN Dengxing, LIANG Hongru, CHANG Yu, LI Hanzhi

    KeyLaboratoryofSoil&WaterConservationandDesertificationCombatingofMinistryofEducation,BeijingForestryUniversity,Beijing100083,China

    Soil and water loss has caused severe land resources depletion and eco-environment degradation in many areas of China, which leads to reduction of local agricultural and industrial productivity. Furthermore, as a severe environmental problem, soil and water loss will affect many respects, such as nature, economy, society and science. Many factors are influencing the degree of soil erosion, Plant vegetation coverage, rainfall intensity and duration, soil type and texture, etc. In this study we chose slope, rainfall intensity and vegetation coverage as impact factors to the degree of soil erosion. Plant ryegrass is one of the feasible measures to avoid soil and water loss. In order to understand the process of runoff movement under the complex surface conditions and the effect of flow and different surface conditions on hydraulics parameters, we utilized the artificial rainfall approach to study the process of runoff and sediment yield and its hydraulic characteristics to explore how ryegrass play a role in the work of soil and water conservation. The experiment was designed with 30 mm/h and 60 mm/h rainfall intensity and slopes of 10° and 20° with different coverage of ryegrass (0, 20%, 40%, 60%, 80% and 100%) in Jiufeng experimental forest farm of Beijing Forestry University. During the simulated rainfall process, the runoff samples were collected at fixed time: in the first ten minutes, the runoff was sampled every two minutes due to the fast increase of runoff amount, and then the runoff was sampled every five minutes in the last fifty minutes as the change of runoff amount went slow. At the end of the experiment, the volume of runoff was measured, and then standing the runoff samples were placed and dried for 24 hours, respectively, to obtain the amount of sediment. During the simulated rainfall process, the velocity of overland flow was measured every ten minutes. Then the surface runoff and hydraulic parameters were calculated with runoff, sediment and surface runoff velocity. By choosing Reynolds number, Froude number and drag coefficient as factors to characterize the water advance and water recession process, the results showed that: the Reynolds number increases with the slope degree and rainfall intensity. The increase of ryegrass could play important roles in both the decrease of the Reynolds number and the Froude number, and the Reynolds number and the Froude number under different coverage are: bare land>20%>40%>60%>80% and bare land>20%>40%>60%>80%, respectively. The drag coefficient and sediment yield rate were well-fitting, and the sediment yield rate decreased as the drag coefficient increased. When the drag coefficient was between zero and one, the growth rate of sediment yield was larger than the drag coefficient. By this study, we know that ryegrass has a great significance to soil and water conservation. Moreover, the increase of ryegrass coverage leads to the enhancement of the soil and water conservation function. The results of the ryegrass regulating runoff and preventing soil erosion in this study have theoretical significance to understanding the process of soil erosion.

    ryegrass; overland flow; simulated rainfall; hydraulics; flow pattern

    國家科技支撐項目(“三北”地區(qū)防護林體系結(jié)構(gòu)定向調(diào)控技術(shù)研究與示范)

    2013-10-30;

    日期:2014-07-25

    10.5846/stxb201310302620

    *通訊作者Corresponding author.E-mail: yuxinxiao11@126.com

    孫佳美, 余新曉, 樊登星, 梁洪儒, 常 玉, 李瀚之.模擬降雨下植被蓋度對坡面流水動力學(xué)特性的影響.生態(tài)學(xué)報,2015,35(8):2574-2580.

    Sun J M, Yu X X, Fan D X, Liang H R, Chang Y, Li H Z.Impact of vegetation cover on surface runoff hydraulic characteristics with simulated rainfall.Acta Ecologica Sinica,2015,35(8):2574-2580.

    猜你喜歡
    流態(tài)蓋度黑麥草
    認識草,也是一門學(xué)問
    側(cè)邊機組故障對泵站前池流態(tài)的影響
    不同行內(nèi)生草影響葡萄果實品質(zhì)
    黃山市近30 a植被蓋度時空變化遙感分析
    黑麥草的種植、青貯及應(yīng)用
    黃土高原地區(qū)植被蓋度對產(chǎn)流產(chǎn)沙的影響
    種苗根熒光反應(yīng)在多年生黑麥草與一年生黑麥草兩個種測定中可靠性驗證
    改進邊界條件的非恒定流模型在城市河流橡膠壩流態(tài)模擬中的應(yīng)用
    坡位與植被蓋度對楓香造林成效的影響
    防護林科技(2015年5期)2015-06-10 11:42:52
    基于數(shù)學(xué)模型研究三亞市區(qū)綠地生態(tài)破壞后蓋度恢復(fù)
    国产精品女同一区二区软件| 亚洲情色 制服丝袜| 国产又色又爽无遮挡免| 成人无遮挡网站| 哪个播放器可以免费观看大片| 日日摸夜夜添夜夜添av毛片| av在线播放精品| 蜜臀久久99精品久久宅男| 欧美 日韩 精品 国产| 亚洲一级一片aⅴ在线观看| 肉色欧美久久久久久久蜜桃| 中文字幕人妻熟人妻熟丝袜美| 一二三四中文在线观看免费高清| 制服诱惑二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美bdsm另类| 女的被弄到高潮叫床怎么办| 国产国拍精品亚洲av在线观看| 免费观看在线日韩| 成年av动漫网址| 日韩三级伦理在线观看| 久久久久国产网址| 欧美精品亚洲一区二区| 人人妻人人澡人人看| 麻豆成人av视频| 久久久久国产精品人妻一区二区| 日韩av不卡免费在线播放| 国产一区有黄有色的免费视频| 高清不卡的av网站| 国产精品偷伦视频观看了| 久久毛片免费看一区二区三区| 一本一本综合久久| 成人国语在线视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久国产亚洲av麻豆专区| 亚洲熟女精品中文字幕| 国产精品.久久久| 欧美精品高潮呻吟av久久| 丝袜脚勾引网站| 亚洲精品久久午夜乱码| 精品国产一区二区久久| 亚洲欧美精品自产自拍| 美女中出高潮动态图| 日韩一区二区三区影片| 亚洲,一卡二卡三卡| 3wmmmm亚洲av在线观看| 精品国产露脸久久av麻豆| 欧美一级a爱片免费观看看| 亚洲av免费高清在线观看| 欧美另类一区| 国产精品一区二区在线不卡| 国产精品秋霞免费鲁丝片| 波野结衣二区三区在线| 欧美一级a爱片免费观看看| 狂野欧美激情性bbbbbb| 在线观看国产h片| 国产精品一国产av| 99精国产麻豆久久婷婷| av不卡在线播放| 一级黄片播放器| 国产男女内射视频| 老司机亚洲免费影院| 日本猛色少妇xxxxx猛交久久| 精品久久久久久久久亚洲| 大香蕉久久成人网| 男女啪啪激烈高潮av片| 一级毛片aaaaaa免费看小| av国产久精品久网站免费入址| 男女国产视频网站| 我的女老师完整版在线观看| 国产极品天堂在线| 你懂的网址亚洲精品在线观看| 免费看光身美女| 日韩欧美一区视频在线观看| 91精品一卡2卡3卡4卡| 青春草视频在线免费观看| 国产亚洲一区二区精品| 热99国产精品久久久久久7| 精品少妇内射三级| 久久97久久精品| 国产亚洲最大av| 亚洲一级一片aⅴ在线观看| 欧美bdsm另类| 久久精品国产自在天天线| 亚洲精品久久久久久婷婷小说| 国产高清不卡午夜福利| 亚洲国产日韩一区二区| 在线播放无遮挡| 久久婷婷青草| 日日啪夜夜爽| 久久精品久久久久久久性| 丝袜脚勾引网站| 中文字幕精品免费在线观看视频 | 日韩人妻高清精品专区| 日韩熟女老妇一区二区性免费视频| 亚洲精品一二三| 亚洲人成网站在线观看播放| 欧美日韩一区二区视频在线观看视频在线| 午夜福利在线观看免费完整高清在| 午夜福利,免费看| 国产又色又爽无遮挡免| 少妇精品久久久久久久| 成人综合一区亚洲| 精品久久久噜噜| 精品99又大又爽又粗少妇毛片| 亚洲,欧美,日韩| 欧美一级a爱片免费观看看| 欧美成人午夜免费资源| 18禁动态无遮挡网站| 999精品在线视频| 国产亚洲一区二区精品| 丰满饥渴人妻一区二区三| 美女国产高潮福利片在线看| 成人无遮挡网站| 久久精品久久精品一区二区三区| 男女边吃奶边做爰视频| www.色视频.com| 国产精品一国产av| 亚洲精品亚洲一区二区| 日韩一本色道免费dvd| 熟女电影av网| 女性生殖器流出的白浆| 欧美激情 高清一区二区三区| 国产日韩欧美亚洲二区| 性色av一级| 亚洲欧美中文字幕日韩二区| 国产在视频线精品| 一区二区三区四区激情视频| 在线观看美女被高潮喷水网站| 国产精品三级大全| 高清av免费在线| 美女xxoo啪啪120秒动态图| 丰满乱子伦码专区| 美女主播在线视频| 日韩一本色道免费dvd| 看十八女毛片水多多多| 日韩成人伦理影院| 精品亚洲成国产av| 免费看av在线观看网站| 成年人午夜在线观看视频| 青青草视频在线视频观看| 大码成人一级视频| 菩萨蛮人人尽说江南好唐韦庄| 久久99一区二区三区| 99国产综合亚洲精品| 热99久久久久精品小说推荐| 精品人妻一区二区三区麻豆| 飞空精品影院首页| 久久人人爽人人片av| 色婷婷久久久亚洲欧美| 亚洲国产最新在线播放| 啦啦啦视频在线资源免费观看| 男女高潮啪啪啪动态图| 啦啦啦在线观看免费高清www| 色婷婷av一区二区三区视频| 观看av在线不卡| 国产 一区精品| 国产 一区精品| 日韩欧美精品免费久久| 成人无遮挡网站| 简卡轻食公司| 国产乱来视频区| 国产亚洲一区二区精品| 超碰97精品在线观看| 国产精品国产三级国产专区5o| 午夜激情福利司机影院| 国产 精品1| 九色亚洲精品在线播放| 在线观看www视频免费| 国产精品久久久久久精品古装| 久久人人爽人人片av| 丰满饥渴人妻一区二区三| 日本黄色片子视频| 欧美成人精品欧美一级黄| 18禁动态无遮挡网站| 熟妇人妻不卡中文字幕| 精品亚洲成国产av| 一级爰片在线观看| 另类亚洲欧美激情| 伦精品一区二区三区| 岛国毛片在线播放| 少妇高潮的动态图| 插逼视频在线观看| a级毛片在线看网站| 亚洲av中文av极速乱| 精品久久久久久久久亚洲| 一本色道久久久久久精品综合| 亚洲欧洲日产国产| 18禁在线播放成人免费| 国产成人免费无遮挡视频| 精品久久久精品久久久| 亚洲五月色婷婷综合| 亚洲精品中文字幕在线视频| 女的被弄到高潮叫床怎么办| 日韩,欧美,国产一区二区三区| 99视频精品全部免费 在线| 天天躁夜夜躁狠狠久久av| 2022亚洲国产成人精品| 亚洲欧洲国产日韩| 国产精品免费大片| 日本爱情动作片www.在线观看| 日韩免费高清中文字幕av| 在线免费观看不下载黄p国产| 亚洲av不卡在线观看| 精品久久久噜噜| 亚洲国产精品999| 欧美最新免费一区二区三区| 永久免费av网站大全| 一区二区三区精品91| 亚洲国产精品专区欧美| 一级毛片电影观看| 99久久综合免费| 国产成人aa在线观看| 亚洲av不卡在线观看| 少妇精品久久久久久久| 欧美人与善性xxx| 国产精品久久久久久久电影| 最新的欧美精品一区二区| 日本色播在线视频| 国产精品国产三级专区第一集| 久久精品久久久久久噜噜老黄| 色婷婷av一区二区三区视频| 有码 亚洲区| 9色porny在线观看| 美女内射精品一级片tv| 女性被躁到高潮视频| 在线观看www视频免费| 久久精品国产a三级三级三级| 久久久a久久爽久久v久久| 国产av一区二区精品久久| 考比视频在线观看| 亚洲国产日韩一区二区| 国产在线免费精品| 人妻夜夜爽99麻豆av| 只有这里有精品99| 国产极品粉嫩免费观看在线 | 午夜久久久在线观看| 国产精品一区二区在线不卡| 丝袜美足系列| 国产成人aa在线观看| 18禁观看日本| 免费少妇av软件| 尾随美女入室| 亚洲成色77777| 天天操日日干夜夜撸| 大香蕉久久网| 亚洲精品乱码久久久久久按摩| 成人国产麻豆网| 国产成人av激情在线播放 | 美女国产高潮福利片在线看| 啦啦啦视频在线资源免费观看| 男男h啪啪无遮挡| 亚洲欧美日韩卡通动漫| 午夜福利视频精品| 视频中文字幕在线观看| 一级片'在线观看视频| 午夜久久久在线观看| 久久人人爽人人爽人人片va| 2022亚洲国产成人精品| 亚洲av日韩在线播放| 男人操女人黄网站| 18禁裸乳无遮挡动漫免费视频| 2018国产大陆天天弄谢| 视频中文字幕在线观看| 最近手机中文字幕大全| 99热这里只有是精品在线观看| 国产欧美日韩一区二区三区在线 | 免费av不卡在线播放| 成年美女黄网站色视频大全免费 | 精品国产国语对白av| 99热全是精品| 精品人妻一区二区三区麻豆| 久久亚洲国产成人精品v| 在线精品无人区一区二区三| 精品国产露脸久久av麻豆| 国产精品成人在线| 制服人妻中文乱码| 亚洲四区av| 另类精品久久| 一区二区日韩欧美中文字幕 | 曰老女人黄片| 亚洲不卡免费看| 日韩精品免费视频一区二区三区 | 少妇熟女欧美另类| 亚洲av成人精品一区久久| 在线精品无人区一区二区三| 亚洲精品aⅴ在线观看| 亚洲熟女精品中文字幕| 国产亚洲午夜精品一区二区久久| 能在线免费看毛片的网站| 国产日韩欧美视频二区| 男女无遮挡免费网站观看| 少妇高潮的动态图| 欧美人与性动交α欧美精品济南到 | 国产欧美日韩一区二区三区在线 | 国产一区二区三区av在线| 18禁裸乳无遮挡动漫免费视频| 国产综合精华液| 国产精品久久久久久精品电影小说| 亚洲国产成人一精品久久久| 大又大粗又爽又黄少妇毛片口| 国产成人精品在线电影| 国产精品 国内视频| 桃花免费在线播放| 国产成人av激情在线播放 | 久久久久久久大尺度免费视频| 亚洲欧美成人精品一区二区| 久久久久久久久久久久大奶| 桃花免费在线播放| 高清不卡的av网站| 建设人人有责人人尽责人人享有的| 免费少妇av软件| 国产有黄有色有爽视频| 日本vs欧美在线观看视频| 最近中文字幕高清免费大全6| 天堂俺去俺来也www色官网| 91成人精品电影| 波野结衣二区三区在线| 亚洲精品日韩av片在线观看| 一级a做视频免费观看| 亚洲精品av麻豆狂野| 国产乱来视频区| 国产色婷婷99| 色94色欧美一区二区| 国产黄色免费在线视频| 午夜激情久久久久久久| 王馨瑶露胸无遮挡在线观看| 亚洲精品第二区| 免费少妇av软件| 亚洲第一区二区三区不卡| 成人黄色视频免费在线看| 伦理电影免费视频| 国产片特级美女逼逼视频| 丝袜喷水一区| 免费av不卡在线播放| 欧美人与性动交α欧美精品济南到 | 91午夜精品亚洲一区二区三区| 亚洲av男天堂| 国产高清有码在线观看视频| 秋霞在线观看毛片| 欧美97在线视频| .国产精品久久| 国产 一区精品| 中国国产av一级| 久久久久国产精品人妻一区二区| 国产一区有黄有色的免费视频| 人妻制服诱惑在线中文字幕| 日产精品乱码卡一卡2卡三| 大码成人一级视频| 少妇猛男粗大的猛烈进出视频| 少妇人妻精品综合一区二区| 日本91视频免费播放| 国产高清不卡午夜福利| 国产精品偷伦视频观看了| 色哟哟·www| 一本大道久久a久久精品| 欧美成人午夜免费资源| 人妻制服诱惑在线中文字幕| 日产精品乱码卡一卡2卡三| 国国产精品蜜臀av免费| 国产片特级美女逼逼视频| 日本欧美视频一区| 制服丝袜香蕉在线| 美女主播在线视频| 午夜免费鲁丝| 在线观看www视频免费| av福利片在线| 丰满乱子伦码专区| 2021少妇久久久久久久久久久| 在线观看人妻少妇| 欧美日韩一区二区视频在线观看视频在线| 日日撸夜夜添| 性色avwww在线观看| 老女人水多毛片| 制服丝袜香蕉在线| 欧美另类一区| 一级黄片播放器| 欧美激情国产日韩精品一区| 日本午夜av视频| 卡戴珊不雅视频在线播放| 亚洲不卡免费看| 中文字幕人妻丝袜制服| 久久人人爽人人爽人人片va| 国产成人aa在线观看| 色5月婷婷丁香| 制服诱惑二区| 久久韩国三级中文字幕| 亚洲少妇的诱惑av| 高清黄色对白视频在线免费看| 啦啦啦在线观看免费高清www| 美女国产高潮福利片在线看| 两个人免费观看高清视频| 91午夜精品亚洲一区二区三区| 亚洲情色 制服丝袜| 狠狠精品人妻久久久久久综合| 欧美日韩综合久久久久久| 亚洲av福利一区| 在线免费观看不下载黄p国产| 人妻一区二区av| 欧美bdsm另类| 插阴视频在线观看视频| 国产淫语在线视频| 汤姆久久久久久久影院中文字幕| 久久精品久久精品一区二区三区| 九九爱精品视频在线观看| 免费观看在线日韩| 国产精品三级大全| 亚洲av综合色区一区| 亚洲经典国产精华液单| 亚洲一级一片aⅴ在线观看| 建设人人有责人人尽责人人享有的| 一个人看视频在线观看www免费| 少妇人妻 视频| 老司机影院成人| 国产不卡av网站在线观看| 日本午夜av视频| 亚洲五月色婷婷综合| 免费av不卡在线播放| 一本—道久久a久久精品蜜桃钙片| 亚洲四区av| kizo精华| 久久精品国产自在天天线| 18禁在线无遮挡免费观看视频| 黄色欧美视频在线观看| 午夜日本视频在线| av在线播放精品| 少妇精品久久久久久久| 亚洲精品自拍成人| 中文字幕最新亚洲高清| 美女视频免费永久观看网站| 大香蕉久久网| 女人精品久久久久毛片| 狂野欧美激情性bbbbbb| 熟女av电影| 国产极品粉嫩免费观看在线 | 亚洲国产色片| 在线播放无遮挡| 精品久久久久久久久亚洲| 中文字幕精品免费在线观看视频 | 国产深夜福利视频在线观看| 欧美最新免费一区二区三区| 亚洲一级一片aⅴ在线观看| 美女国产高潮福利片在线看| 国产伦理片在线播放av一区| 国产欧美日韩综合在线一区二区| 亚洲四区av| 亚洲精品美女久久av网站| 亚洲精品色激情综合| 久久精品久久久久久久性| 大香蕉久久网| 九色成人免费人妻av| 久久久精品区二区三区| 国产成人精品福利久久| 九色成人免费人妻av| 午夜视频国产福利| 久久久久久久久久久免费av| 大码成人一级视频| 韩国高清视频一区二区三区| 在线天堂最新版资源| 97在线视频观看| 亚洲av不卡在线观看| 最近中文字幕2019免费版| 2018国产大陆天天弄谢| 久久99热6这里只有精品| 黄色一级大片看看| 十分钟在线观看高清视频www| 国产成人aa在线观看| 国产精品99久久久久久久久| 精品久久久久久电影网| 国产精品熟女久久久久浪| kizo精华| 日韩制服骚丝袜av| 美女国产高潮福利片在线看| 亚洲,欧美,日韩| 亚洲综合色惰| 啦啦啦啦在线视频资源| 亚洲av二区三区四区| 国产精品人妻久久久影院| 99热6这里只有精品| 精品少妇黑人巨大在线播放| 色吧在线观看| 麻豆精品久久久久久蜜桃| 国产白丝娇喘喷水9色精品| 大码成人一级视频| 妹子高潮喷水视频| 精品国产露脸久久av麻豆| 99九九线精品视频在线观看视频| 精品99又大又爽又粗少妇毛片| 日本免费在线观看一区| 男女高潮啪啪啪动态图| 国产一级毛片在线| 国产精品久久久久久久电影| 男女国产视频网站| 中文字幕免费在线视频6| 亚洲欧美一区二区三区黑人 | 亚洲婷婷狠狠爱综合网| 亚洲精品久久成人aⅴ小说 | 亚洲av国产av综合av卡| 国产免费视频播放在线视频| 久久久国产欧美日韩av| 日韩av不卡免费在线播放| 亚洲精品久久午夜乱码| 97在线人人人人妻| 这个男人来自地球电影免费观看 | 大陆偷拍与自拍| a级毛片免费高清观看在线播放| 欧美丝袜亚洲另类| 天堂中文最新版在线下载| 国产综合精华液| 久久这里有精品视频免费| 国产成人免费观看mmmm| 中文字幕亚洲精品专区| 蜜臀久久99精品久久宅男| 国产精品嫩草影院av在线观看| 久久鲁丝午夜福利片| 夫妻午夜视频| 亚洲五月色婷婷综合| 超碰97精品在线观看| 亚洲精华国产精华液的使用体验| 国产成人精品无人区| 视频在线观看一区二区三区| 少妇人妻 视频| 中国国产av一级| 亚洲国产精品成人久久小说| 国产熟女欧美一区二区| 久久人人爽av亚洲精品天堂| 蜜桃在线观看..| 亚洲国产欧美在线一区| 在线亚洲精品国产二区图片欧美 | 国产成人a∨麻豆精品| 伊人久久精品亚洲午夜| 91久久精品国产一区二区三区| 国产成人精品一,二区| 又大又黄又爽视频免费| 蜜桃久久精品国产亚洲av| 人人妻人人澡人人看| 嘟嘟电影网在线观看| 精品人妻偷拍中文字幕| 国内精品宾馆在线| 国产熟女欧美一区二区| 人人妻人人澡人人看| 国产亚洲一区二区精品| 国产女主播在线喷水免费视频网站| av国产久精品久网站免费入址| 欧美激情极品国产一区二区三区 | 欧美少妇被猛烈插入视频| 久久久国产欧美日韩av| 国产一区二区在线观看av| 18禁动态无遮挡网站| 国产精品无大码| 国产精品一二三区在线看| 亚洲精品乱码久久久久久按摩| 在线 av 中文字幕| 男女免费视频国产| 大香蕉97超碰在线| 欧美日韩一区二区视频在线观看视频在线| 丝袜美足系列| 亚洲av二区三区四区| 一二三四中文在线观看免费高清| 久久久久国产网址| 国精品久久久久久国模美| 国产伦理片在线播放av一区| 国产片内射在线| 热99久久久久精品小说推荐| 伦理电影免费视频| av视频免费观看在线观看| 男女国产视频网站| 久久毛片免费看一区二区三区| 一级a做视频免费观看| 久久狼人影院| 亚洲精品色激情综合| 涩涩av久久男人的天堂| 99久久中文字幕三级久久日本| 久久久久久久久久久丰满| 亚洲av.av天堂| 丁香六月天网| 中国三级夫妇交换| 一级爰片在线观看| 免费久久久久久久精品成人欧美视频 | 91aial.com中文字幕在线观看| 最近2019中文字幕mv第一页| 交换朋友夫妻互换小说| 美女国产高潮福利片在线看| 美女大奶头黄色视频| 美女国产高潮福利片在线看| av网站免费在线观看视频| a 毛片基地| 亚洲精品乱码久久久久久按摩| 九色亚洲精品在线播放| 日韩av不卡免费在线播放| 久久精品夜色国产| 视频在线观看一区二区三区| 丰满饥渴人妻一区二区三| 国产一区二区三区av在线| 91精品国产九色| 国产国语露脸激情在线看| 精品99又大又爽又粗少妇毛片| 欧美bdsm另类| 国产深夜福利视频在线观看| 国产成人91sexporn| 免费黄频网站在线观看国产| 中国三级夫妇交换| 免费少妇av软件| 精品久久国产蜜桃| 少妇被粗大的猛进出69影院 | 成人毛片a级毛片在线播放| 亚洲人与动物交配视频| 亚洲av国产av综合av卡| 亚洲精品456在线播放app| 最新中文字幕久久久久| 久久精品国产a三级三级三级| www.色视频.com| 久久久久久久亚洲中文字幕| 一级a做视频免费观看| 国产精品人妻久久久久久|