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

    黃土區(qū)工程堆積體陡坡坡面徑流調(diào)控工程措施的減沙效應

    2019-09-24 11:34:38張樂濤董俊武朱雅琴黃增玉李婉璐王慧慧唐林昊田紅衛(wèi)高照良
    農(nóng)業(yè)工程學報 2019年15期
    關鍵詞:調(diào)沙沙量水沙

    張樂濤,董俊武,袁 琳,朱雅琴,黃增玉,李婉璐,王慧慧,唐林昊,田紅衛(wèi),高照良

    黃土區(qū)工程堆積體陡坡坡面徑流調(diào)控工程措施的減沙效應

    張樂濤1,董俊武1,袁 琳1,朱雅琴1,黃增玉1,李婉璐1,王慧慧1,唐林昊1,田紅衛(wèi)2,高照良3※

    (1. 河南大學環(huán)境與規(guī)劃學院環(huán)境與規(guī)劃國家級實驗教學示范中心,開封 475004; 2. 長江勘測規(guī)劃設計研究有限責任公司,武漢 430010; 3. 西北農(nóng)林科技大學水土保持研究所,楊凌 712100)

    定量分析減少徑流和改變水沙關系在泥沙調(diào)控中的不同作用,對于深刻理解徑流調(diào)控措施的水土保持效益具有重要意義。以黃土區(qū)工程堆積體陡坡坡面(36°)為例,探討了野外模擬徑流沖刷試驗條件下,不同工程措施及其組合調(diào)控坡面徑流的水沙效應及其作用效率。結果表明:1)不同工程措施均能較好地調(diào)控坡面徑流侵蝕過程,不同情形下的產(chǎn)流時間控制比為2~20,徑流量控制比為0.45~0.78,產(chǎn)沙量控制比為0.20~0.59;平均含沙量控制比為0.38~0.79;2)減流控沙作用是工程措施調(diào)控坡面侵蝕產(chǎn)沙的主要原因,水沙關系調(diào)沙作用則受減流控沙作用的制約;3)水平階類措施的水沙關系調(diào)沙量與減流控沙量呈線性正相關,當減流控沙量超過一定臨界值時,水沙關系才開始發(fā)揮調(diào)沙作用;水平溝類措施的水沙關系調(diào)沙量與減流控沙量呈二次函數(shù)關系,水沙關系調(diào)沙量存在極大值;4)水平溝類措施調(diào)控泥沙的作用效率高于水平階類措施,水平溝+魚鱗坑的組合可很好地發(fā)揮減流控沙和水沙關系調(diào)沙作用的潛力(55%),使二者在較高的水平上維持相對平衡;因此,不同工程措施與組合的實際應用應以具體的水土保持效益和防治目標為布設依據(jù)。研究可為堆積體陡坡治理的工程措施優(yōu)化提供理論參考。

    侵蝕;泥沙;徑流;水沙關系;泥沙調(diào)控能力;效率系數(shù);棄土堆積體

    0 引 言

    工程堆積體是由生產(chǎn)建設活動所產(chǎn)生的棄(廢)土、棄(廢)渣等廢棄材料積聚而成的人工技術結構體,是人類地貌活動的產(chǎn)物和典型的人工地貌景觀。從地貌過程的規(guī)模和維度看,工程堆積體是最為壯觀的人造積聚形態(tài)之一[1]。從景觀功能的角度看,工程堆積體的形成是人工選擇作用主導下的低洼地、溝谷甚至農(nóng)田改造為人造山丘的景觀逆轉(zhuǎn)過程。土體、地形等景觀要素的人為重組喪失了自然生態(tài)系統(tǒng)原有景觀功能的整體性和持久性。與自然坡面相比,工程堆積體地邊坡失調(diào)、地表組成物質(zhì)松散、穩(wěn)定性結構體缺失、土壤組分的質(zhì)地及肥力均遭到破壞,植被及有機質(zhì)匱乏、水文地質(zhì)條件惡化,其對降雨、風速等自然營力的抵抗能力下降[2-4]。因此,與自然過程相比,工程堆積體的地表過程發(fā)生速度與強度發(fā)生了明顯改變[5-6],是生產(chǎn)建設區(qū)新增水土流失的重要策源地[7-8]。

    黃土高原地區(qū)生態(tài)環(huán)境脆弱,水分是制約該區(qū)生態(tài)修復的限制性因子。因此,“28字方略”仍是黃土區(qū)工程堆積體土地資源整治的指導性策略,以促進地表徑流攔蓄、增加土壤入滲,為植被的迅速恢復創(chuàng)造水分條件。在此背景下,實施以工程措施為基礎的陡坡治理進而恢復土壤水文功能成為促進該區(qū)工程堆積體生態(tài)修復的關鍵和前提。然而,筆者在野外考察中發(fā)現(xiàn),相當數(shù)量的工程堆積體高陡邊坡在治理、恢復過程中嚴重忽視植被分布的地帶性規(guī)律,片面強調(diào)植被的生態(tài)治理功能,一味盲目地追求提高植被覆蓋度,導致植物成活率低和局部生態(tài)環(huán)境惡化的惡性循環(huán)。

    為合理有效優(yōu)化工程堆積體的管理維護、維持并提升工程堆積體的生態(tài)景觀功能,工程堆積體的微地形整治、生態(tài)系統(tǒng)的恢復與重建日益成為關注焦點與研究熱點[9-15]。水平階、水平溝、魚鱗坑等傳統(tǒng)陡坡治理工程廣泛應用于工程堆積體微地形整治中,其具有縮短徑流流線、降低徑流流速、攔截坡面上方來水,促進局地降水徑流的富集疊加等徑流調(diào)節(jié)功能,進而在減輕土壤沖刷的同時有效改善土壤水文條件[16-22],提高植被存活幾率、恢復生態(tài)環(huán)境。目前,在工程堆積體陡坡土地整治背景下,研究者多聚焦于工程措施減水減沙的綜合效益分析[10,23-25],未能密切關注陡坡條件下工程措施的水沙調(diào)控效率,即徑流與水沙關系變化對坡面侵蝕泥沙輸出的不同調(diào)控作用,因而不能進一步精細地量化表達工程措施的水土保持效益,以為促進坡面水資源富集、提高水資源利用效率、提升堆積體陡坡治理水平奠定理論基礎。

    為進一步探討工程措施對堆積體陡坡坡面侵蝕泥沙的調(diào)節(jié)作用,本文在對“神府高速公路”沿線棄土堆積體全面調(diào)查的基礎上,以黃土區(qū)典型工程堆積體陡坡坡面為例,設計水平階、水平溝、魚鱗坑等不同的工程措施及其組合形式,采用野外放水沖刷的試驗方法,對工程措施及其組合影響下的徑流減沙(因徑流量改變引起的減沙)和水沙關系調(diào)沙(因水沙關系改變引起的減沙)過程進行了研究,以期為堆積體陡坡工程治理過程中的措施優(yōu)化與管理維護提供參考依據(jù)。

    1 材料與方法

    1.1 試驗區(qū)概況

    試驗區(qū)位于黃河的一級支流—石馬川小流域,流域面積為243 km2。該區(qū)原生地貌主要為黃土丘陵溝壑,屬中溫帶半干旱大陸性季風氣候,處于440 mm等降水量線附近,降水量不足,降水主要集中在6-9月,時間分布不均,夏季平均降水284.6 mm,占全年62.77%,極易形成山洪。本區(qū)氣候干燥,日照充足,年平均氣溫8.0~8.5 ℃,無霜期平均173 d。石馬川小流域水量貧乏,年徑流總量約為0.2144億m3,年輸沙量約為597萬t。原生植主要包括檸條(Kom)、沙棘(Linn.)、酸棗(Mill. var.(Bunge) Hu ex H. F. Chow)等落葉闊葉灌叢,原生土壤類型主要為黃綿土和風沙土,該流域土壤侵蝕形式以水蝕為主,多年平均年土壤侵蝕模數(shù)為3.04×104t/(km2·a),水土流失嚴重[26]。

    神木至府谷高速公路在建設過程中產(chǎn)生頂部平臺面積大于250 m2的工程堆積體120座,坡度變化范圍為30°~50°,陡坡斜長變化范圍為5~80 m,其中,35°~40°和10~25 m的堆積體數(shù)量最多,占比分別為67.5%(81座)和62.5%(75座)。

    布設徑流試驗小區(qū)的工程堆積體位于神木—府谷高速公路沿線的一處棄土場(38°56′15.5″N,110°53′10.0″E),由依溝傾倒棄土形成,距石馬川收費站約0.8 km,海拔高度966 m。堆積體為“頂部平臺+高陡邊坡”結構,于2012年達到堆積的穩(wěn)定狀態(tài)。平臺面積約為2.8 hm2,邊坡平均坡度為36°,陡坡平均斜長為78 m,坡面未布設水土保持措施。自然恢復的植被主要以沙蓬((Linn.)Moq)、沙柳()等沙地植被為主,植被覆蓋率低。坡面地表組成物質(zhì)以黃綿土為主,土石比超過9:1,粒徑多在1 mm以下(表1)。

    表1 土壤顆粒組成

    1.2 試驗設計與觀測

    徑流小區(qū)的頂端與頂部平臺的距離為63 m,沿坡長方向布設,長12 m,寬2.5 m,坡度為36°。工程措施布設方式及規(guī)格如下:

    1)為紊亂坡面徑流流線、改變徑流運動方向,緩解徑流強度,沿坡寬方向修筑的長2.5 m、寬30 cm的水平階,布設方式為:①均勻型,4個水平階,沿坡長方向的間距為3 m;②增加型(間距),4個水平階,沿坡長方向的間距分別為1、2、3、4 m;

    2)為存蓄坡面徑流、攔截侵蝕泥沙,沿坡寬方向修建長2.5 m、頂寬30 cm、深10 cm的水平溝,斷面為梯形,布設方式為:①均勻型,4條水平溝,沿坡長方向的間距為3 m;②增加型(間距),4條水平溝,沿坡長方向的間距分別為1、2、3、4 m;

    3)工程措施組合方式:①水平階+水平溝組合(規(guī)格與前述統(tǒng)一),按“水平階+水平溝+水平階+水平溝”的格局布設,沿坡長方向的間距分別為1、2、3、4 m;②水平階+魚鱗坑組合,水平階的規(guī)格與前述統(tǒng)一,魚鱗坑的半徑為25 cm,深為10 cm,魚鱗坑呈“品”字形布設,總體沿坡長方向按“水平階+魚鱗坑+魚鱗坑+水平溝”的格局布設,沿坡長方向的間距分別為1、2、3、4 m;③水平溝+魚鱗坑組合(水平溝、魚鱗坑的規(guī)格與前述統(tǒng)一),沿坡長方向按“水平溝+魚鱗坑+魚鱗坑+水平溝”的格局布設,沿坡長方向的間距分別為1、2、3、4 m。為對比試驗結果,設置裸露坡面的對照試驗,不布設任何措施。試驗過程中,共有8種試驗情形產(chǎn)生有效觀測,各場次試驗的放水時間視情況持續(xù)時間約為30 min,各情形設置重復一次。

    模擬放水試驗的放水流量按如下方式計算:放水流量=(模擬寬度/實際溝間距)×降水強度×匯水面積×產(chǎn)流系數(shù),徑流系數(shù)取值0.65;最大1 h降水強度取100 mm/h,溝蝕發(fā)生的臨界匯水面積變化取300 m2,溝間距30 m,試驗中的模擬寬度為2.5 m,利用上述算法計算得出放水流量為27 L/min。因此,根據(jù)前期試驗結果[27]和該區(qū)暴雨產(chǎn)流在試驗小區(qū)產(chǎn)生的單寬流量,設計供水流量為25 L/min。試驗開始前,利用流量計和用標有刻度的徑流桶率定放水流量2~3次,沿坡長方向設置4個等間距的觀測斷面(圖2)。試驗開始后,待小區(qū)出口產(chǎn)生徑流時,記錄出流時間,并在小區(qū)出口處接取徑流泥沙樣品。徑流泥沙樣在產(chǎn)流后3 min內(nèi)每1 min取1次,3 min后每隔2 min取1次。徑流桶的半徑為0.49 m,高1 m,不設分流設施,徑流桶接收全部泥沙和徑流。徑流桶內(nèi)設置一把量程為150 mm、精度為1 mm的直尺,根據(jù)徑流桶內(nèi)的水位變化過程推算小區(qū)徑流、泥沙過程。泥沙樣品烘干后用天平稱其質(zhì)量(精度為0.01 g),徑流均換算為清水體積,含沙量為渾水樣品中泥沙質(zhì)量與清水體積的比值。試驗結束后,將小區(qū)全部進行翻新、填土,重新布設工程措施,盡可能使不同場次試驗的地表狀況一致,表面進行灑水(翻新期間如遇晚間降雨,則不灑水),晾置12 h后,重新進行試驗。試驗期間,土壤容重變化范圍為1.32~1.41 g/cm3,前期質(zhì)量含水率變化范圍為13.6%~14.6%。

    試驗系統(tǒng)、試驗過程和觀測參照參考文獻[2,7],不同工程措施斷面的設計規(guī)格及坡面措施布局如圖1和圖2所示。

    圖1 工程措施的橫斷面設計

    圖2 坡面工程措施布設

    1.3 數(shù)據(jù)分析方法

    1)控制比()

    為表征工程措施對坡面徑流侵蝕的實際控制作用,采用控制比反映工程措施實施前后堆積體坡面徑流、侵蝕泥沙的變化情況。

    式中表示產(chǎn)流時間,min;徑流量,m3;產(chǎn)沙量,kg;平均徑流含沙量,kg/m3等觀測指標,1表示裸露坡面的觀測結果,2表示不同工程措施情形下的觀測結果。

    2)泥沙調(diào)控能力()

    徑流驅(qū)動的侵蝕泥沙輸出過程受徑流量和水沙關系的制約[20],水土保持措施調(diào)節(jié)下的侵蝕泥沙變化可以分解為兩種來源:徑流量變化引起的侵蝕泥沙變化和水沙關系改變引起的侵蝕泥沙變化[28-30]。

    工程措施實施前后的產(chǎn)沙量變化可用下式表示

    式中Δ為工程措施發(fā)揮作用后的總減沙量,kg;1為工程措施實施前的總產(chǎn)沙量,kg;2為工程措施實施后的總產(chǎn)沙量,kg;Δ為工程措施發(fā)揮作用后的產(chǎn)流變化總量,m3;1為工程措施實施前的總產(chǎn)流量,m3;2為工程措施實施后的總產(chǎn)流量,m3;Δ為工程措施發(fā)揮作用后徑流含沙量的變化量,kg/m3;1為工程措施實施前的徑流平均含沙量,kg/m3;2為工程措施實施后的徑流平均含沙量,kg/m3。

    式(2)中,1Δ項可被認為是徑流量改變所引起的泥沙量的變化,定義為減流控沙量(SR);而Δ2項則可被認為是徑流含沙量也即水沙關系改變引起的泥沙量的變化,定義為水沙關系調(diào)沙量(SR)。

    式(1)左右兩邊同時除以Δ,可得

    式(3)中,定義SRC、1和2Δ/Δ分別為不同措施條件下單位徑流的泥沙調(diào)控能力、減流控沙能力(SRC)和水沙關系調(diào)沙能力(SRC),kg/m3。

    式(3)表明,不同措施條件下單位徑流對泥沙的調(diào)控作用表現(xiàn)為兩個方面:一為對原始坡面含沙水流的阻滯,包括攔蓄泥沙、增加“沙匯”和減少徑流量,表現(xiàn)為減流控沙能力;一為降低徑流含沙量,改善水沙傳遞關系,表現(xiàn)為水沙關系調(diào)沙能力。

    3)泥沙調(diào)控效率系數(shù)(、)

    令Δ=1,即為徑流量改變所能引起的減沙潛力,定義為減流控沙潛力

    SRP

    R=

    S

    1

    Q

    1

    (4)

    令Δ=1,此時即為徑流含沙量即水沙關系改變所能引起的泥沙調(diào)節(jié)潛力,定義為水沙關系調(diào)沙潛力

    SRP

    S

    =

    S

    1

    Q

    2

    (5)

    式中SRP、SRP分別表示徑流的控沙潛力和水沙關系的調(diào)沙潛力,kg。

    因此,式(2)可變?yōu)?/p>

    本文定義、分別為減流控沙和水沙關系調(diào)沙作用的效率系數(shù)。

    數(shù)據(jù)分析使用SPSS 18.0,繪圖使用SigmaPlot 13.0。

    2 結果與分析

    2.1 不同工程措施的徑流泥沙調(diào)控效益

    受工程措施特點和布設方式的影響,不同工程措施配置對徑流、侵蝕的控制比差距較大(表2)。水平溝+魚鱗坑組合對延遲坡面出流歷時(水流從坡頂(進水端)匯流至小區(qū)出口(出水端)所需要的時間)、控制產(chǎn)水量的作用最佳(控制比分別為20.0、0.45),而均勻型水平溝攔蓄坡面泥沙、降低徑流含沙量的效果最佳(控制比分為0.20、0.38);相反地,均勻型水平階對延長坡面出流歷時、減小坡面徑流量的作用最差(控制比分別為2.0、0.78),而增加型水平階控制坡面產(chǎn)沙量、減小徑流含沙量的效果最差(控制比分別為0.59、0.79)。上述結果表明,與水平階類措施相比,水平溝類措施的坡面徑流泥沙調(diào)控效益較佳。

    簡單回歸分析發(fā)現(xiàn),在工程措施的調(diào)控下,隨著CR的提高,CR、CR、CR呈下降趨勢,表現(xiàn)出顯著的負相關關系。出流時間、徑流量、產(chǎn)沙量、平均含沙量間的控制比存在如下相關關系

    CR

    W

    =1.23CR

    Q

    –0.37

    R

    2

    =0.89 (10)

    式中CR、CR、CR、CR分別表示工程措施調(diào)控下出流時間、徑流量、產(chǎn)沙量和平均含沙量的控制比,擬合的數(shù)據(jù)樣本量為7。

    表2 不同工程措施影響下的坡面徑流侵蝕響應

    2.2 減流控沙和水沙關系調(diào)沙動態(tài)過程

    根據(jù)式(2)分別計算不同工程措施配置影響下不同觀測時段內(nèi)的減流控沙量和水沙關系調(diào)沙量,得到減流控沙量和水沙關系調(diào)沙量隨時間變化的變化趨勢(圖3)。

    不同工程措施布局情形影響下的減流控沙量和水沙關系調(diào)沙量呈現(xiàn)出不同的時間變化特征:1)減流控沙量隨時間呈現(xiàn)先迅速增大后迅速減小的單峰“拋物線”變化趨勢,控沙量峰值均出現(xiàn)在坡面產(chǎn)流后的9~12 min(以對照組產(chǎn)流時間計),且減流控沙量均為正值;2)水沙關系調(diào)沙量隨時間呈現(xiàn)非單峰的波動變化趨勢,不同措施組合的水沙關系調(diào)沙量峰值出現(xiàn)在坡面產(chǎn)流后的不同時段,水平階類措施的水沙關系調(diào)沙量峰值出現(xiàn)在坡面產(chǎn)流后的6~12 min,水平溝類措施的水沙關系調(diào)沙量峰值出現(xiàn)在坡面產(chǎn)流后的15~18 min(以對照組產(chǎn)流時間計),而到試驗后期,水平階類措施的水沙關系調(diào)沙量出現(xiàn)負值;3)針對不同的措施組合,減流控沙量始終明顯高于同時段的水沙關系調(diào)沙量。上述現(xiàn)象表明,不同措施組合通過控制徑流減少侵蝕泥沙輸出的作用過程類似,對減流控沙量的影響具有一致性(如峰現(xiàn)時間、總體趨勢等);而通過調(diào)節(jié)水流攜沙力減少泥沙輸移數(shù)量的作用過程則因不同措施組合而有不同,對水沙關系調(diào)沙量的影響具有不確定性(如峰現(xiàn)時間、正負效應等);通過控制徑流減少泥沙,也即發(fā)揮減流控沙能力是不同工程措施組合發(fā)揮泥沙調(diào)控作用的主要方面。

    圖3 不同工程措施下減流控沙量和水沙關系調(diào)沙量的時間變化特征

    通過回歸分析可得到不同工程措施布局情形下水沙關系調(diào)沙量與減流控沙量之間的數(shù)量關系(表3)。

    表3 水沙關系調(diào)沙量與減流控沙量間的回歸關系

    注:回歸結果均在=0.01水平上具有顯著性。

    Note: All the regression results are statistically significant at 0.01 levels.

    表3結果顯示,水平階類措施影響下的水沙關系調(diào)沙量與減流控沙量的變化同頻(正相關),由受水平階單一影響的線性關系轉(zhuǎn)變?yōu)槭芩诫A+水平溝、水平階+魚鱗坑復合影響的對數(shù)關系,表明此類措施及其組合情形下徑流量的調(diào)控對水沙關系的改善具有重要影響?;貧w關系表明,當減流控沙量達到某一臨界時,水沙關系的調(diào)節(jié)才能發(fā)揮調(diào)控泥沙的作用。試驗條件下,均勻型水平階、增加型水平階、水平階+水平溝和水平階+魚鱗坑四種情形下這一減流控沙量的臨界值分別為7.5、16.8、14.6、25 kg。

    水平溝類措施影響下的水沙關系調(diào)沙量與減流控沙量的變化則出現(xiàn)異頻現(xiàn)象(呈現(xiàn)先增大后減小的單峰拋物線趨勢),表明此類措施及其組合情形下當徑流的調(diào)節(jié)功能達到某一臨界時,水沙關系的調(diào)沙功能達到極限,當減流控沙量超越這一臨界時,水沙關系調(diào)沙量開始下降,此時徑流量的調(diào)控已不能顯著改善坡面水沙關系。試驗條件下,均勻型水平溝、增加型水平溝和水平溝+魚鱗坑3種情形下的臨界減流控沙量值和最大水沙關系調(diào)沙量值分別為:40.5、24.4 kg;47.7、22.5 kg;44.5、31.2 kg。

    2.3 不同工程措施的泥沙調(diào)控效率

    不同措施影響下減流控沙和水沙關系調(diào)沙作用的效率系數(shù)變化及其相互關系列于表4。

    綜合泥沙調(diào)控過程可以發(fā)現(xiàn),水平階類措施影響下減流控沙作用的效率系數(shù)波動范圍與變異性整體低于水平溝類措施,水沙關系調(diào)沙作用的效率系數(shù)波動范圍與變異性則呈現(xiàn)相反的變化(表4),表明水平階類措施能夠維持較為穩(wěn)定的減流控沙功能,而水平溝類措施則能夠維持較為穩(wěn)定的水沙關系調(diào)沙功能。水平階類措施條件下水沙關系調(diào)沙作用的效率系數(shù)與減流控沙作用的效率系數(shù)呈對數(shù)函數(shù)關系;水平溝類措施條件下則呈線性函數(shù)關系。就綜合值而言,水平階類措施條件下的減流控沙和水沙關系調(diào)沙作用均處于較低水平(效率系數(shù)分別為0.22~0.47、0.21~0.28);與水平階類措施相比,水平溝類措施能明顯提高減流控沙作用和水沙關系調(diào)沙作用的作用效率(效率系數(shù)分別為0.48~0.55、0.55~0.62),表明水平溝類措施更容易發(fā)揮減流控沙作用和水沙關系調(diào)沙作用。從、綜合值的對比關系來看,不同措施類型及組合條件下減流控沙與水沙關系調(diào)沙作用的對比關系不同:增加型水平階、水平階+魚鱗坑條件下減流控沙作用的調(diào)控效率高于水沙關系調(diào)沙作用(:>1),均勻型水平階、水平階+水平溝、均勻型水平溝、增加型水平溝條件下減流控沙作用的調(diào)控效率低于水沙關系調(diào)沙作用(:<1),水平溝+魚鱗坑條件下二者作用效率相當(:=1)。

    不同工程措施情形下的減流控沙與水沙關系調(diào)沙占比及其減沙能力列于表5。表5顯示水平溝+魚鱗坑組合情形下減流控沙占比最大,水沙關系調(diào)沙占比最小;均勻型水平階與均勻型水平溝情形下減流控沙占比最小,水沙關系調(diào)沙占比最大。均勻型水平階與水平階+水平溝情形下的單位徑流調(diào)控能力與水沙關系調(diào)沙能力最大,而水平階+魚鱗坑情形下的單位徑流調(diào)控能力與水沙關系調(diào)沙能力最小。

    注:回歸結果均在=0.01水平上具有顯著性。

    Note: All the regression results are statistically significant at 0.01 levels.

    表5 不同工程措施的減沙能力

    3 討 論

    地表徑流是誘發(fā)坡面侵蝕產(chǎn)沙的根本原因和直接動力,坡面徑流調(diào)控是防治土壤侵蝕的關鍵。減流控沙作用和水沙關系調(diào)沙作用是坡面徑流調(diào)控水沙效應的兩個重要方面。各類工程措施及其組合均能通過明顯削減洪峰流量、有效控制徑流總量,降低徑流含沙量、減少侵蝕泥沙輸出。從泥沙調(diào)控的時間變化特征中可以看出,不同工程措施類型影響下的減流控沙量及其在總減沙量中的占比明顯高于水沙關系調(diào)沙量(圖3),當水沙關系調(diào)沙量為負值時,減流控沙量甚至能達到總減沙量的400倍以上;而從綜合作用結果看,不同工程措施的減流控沙量在總減沙量中的占比均超過50%,最低為均勻型水平階的52%,最高為水平溝+魚鱗坑的77%。因此,工程措施影響下的坡面徑流調(diào)控主要是通過發(fā)揮減流控沙作用來減少坡面侵蝕產(chǎn)沙的數(shù)量。

    從作用效率及其對比關系上看,減流控沙作用與水沙關系調(diào)沙作用相互制約、相互影響。除水平溝+魚鱗坑的組合類型外,水平階、水平溝、魚鱗坑(水毀無效)等單一的工程措施均不能在較高的效率水平上維持減流控沙和水沙關系調(diào)沙作用的相對平衡。在此情形下,坡面徑流調(diào)控的減流控沙潛力和水沙關系調(diào)沙潛力均得到了相當程度的發(fā)揮(55%),徑流量和產(chǎn)沙量的控制比也達到最低,是最優(yōu)措施布設方案(調(diào)控效率和調(diào)控效果而言)。

    從不同措施徑流調(diào)控的作用機制上看,水平階類措施影響下的徑流調(diào)控作用主要體現(xiàn)在改變坡面水流流線(如改變運動方向、延長等路徑長度)而非徑流攔蓄,因而消減徑流量的作用并不突出,盡管提高了單位徑流的泥沙調(diào)控能力和水沙關系調(diào)沙能力(水平階+魚鱗坑較低除外),由于減流控沙的作用效率處于較低水平,水沙關系調(diào)沙的作用水平也較低(減流控沙是水沙關系調(diào)沙的前提)。水平溝類措施則與之相反,其徑流調(diào)控作用主要為徑流攔蓄,徑流量的控制比較低,減流控沙作用效率較高。然而,與水平階類措施持續(xù)穩(wěn)定的發(fā)揮調(diào)節(jié)作用不同,水平溝類措施有限的攔蓄庫容決定了水沙關系調(diào)節(jié)作用的上限,措施被蓄滿、淤滿后,水沙關系調(diào)沙作用不再變化,甚至在措施發(fā)生水毀后加速泥沙的釋放。因此,盡管對泥沙的控制比較低、作用效率較高,單位徑流的泥沙調(diào)控能力和水沙關系調(diào)沙能力卻并不突出。因此,如以保有坡面水資源、提高減沙效率為目標,應考慮以布設水平溝類措施、降低水文連通度為主;如以增加坡面產(chǎn)水量、提高單位徑流泥沙調(diào)控能力并達到一定的減沙效果為目標,則應考慮以布設水平階類措施、提升坡面徑流的緩沖效果為主。

    4 結 論

    本文以黃土區(qū)工程堆積體陡坡坡面為例,設計不同的水土保持工程措施,通過野外模擬徑流沖刷試驗的方法,分析了工程措施徑流調(diào)控的水沙效應及其作用效率。主要結論為:1)減少徑流、發(fā)揮減流控沙作用是工程措施抑制坡面侵蝕、調(diào)控泥沙的主要原因,水沙關系的調(diào)沙過程受減流控沙過程的制約;2)水平階類措施條件下水沙關系調(diào)沙量與減流控沙量之間呈線性函數(shù)關系,水平溝類措施下則呈二次函數(shù)關系;3)水平溝+魚鱗坑的組合措施能在較高的效率水平上維持減流控沙和水沙關系調(diào)沙作用的相對平衡;4)水平階類措施下水沙關系調(diào)沙作用的效率系數(shù)與減流控沙作用的效率系數(shù)之間呈對數(shù)函數(shù)關系,水平溝類措施下則呈線性函數(shù)關系;5)水平階類措施的優(yōu)點為單位徑流泥沙調(diào)控能力和水沙關系調(diào)沙能力較大,不足為減流控沙和水沙關系調(diào)沙的作用效率較低;水平溝類措施的優(yōu)點為減流控沙量大、徑流調(diào)控作用效率高,缺點為徑流的消耗數(shù)量大;在生產(chǎn)實踐中,不同工程措施的選擇及組合措施的布設應以具體目標為依據(jù)。

    [1] Anghel T, Balazsi K. Modelling processes on spoil heaps[J]. Revista de geomorfologie, 2005, 7: 133-140.

    [2] Zhang Letao, Gao Zhaoliang, Yang Shiwei, et al. Dynamic processes of soil erosion by runoff on engineered landforms derived from expressway construction: A case study of typical steep spoil heap[J]. Catena, 2015, 128: 108-121.

    [3] Zhang Letao, Gao Zhaoiang, Li Zhanbin, et al. Downslope runoff and erosion response of typical engineered landform to variable inflow rate patterns from upslope[J]. Natural Hazards, 2016, 80(2): 775-796.

    [4] Hancock G R, Crawter D, Fityus S G, et al. The measurement and modelling of rill erosion at angle of repose slopes in mine spoil[J]. Earth Surface Processes and Landforms, 2008, 33: 1006-1020.

    [5] Lu Huizhong, Cao Longxi, Liang Yin, et al. Mineral-leaching chemical transport with runoff and sediment from severely eroded rare-earth tailings in southern China[J]. Solid Earth, 2017, 8(4): 845-855.

    [6] 史東梅,蔣光毅,彭旭東,等. 不同土石比的工程堆積體邊坡徑流侵蝕過程[J]. 農(nóng)業(yè)工程學報,2015,31(17):152-161.

    Shi Dongmei, Jiang Guangyi, Peng Xudong, et al. Runoff erosion process on slope of engineering accumulation with different soil-rock ratio[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(17): 152-161. (in Chinese with English abstract)

    [7] 張樂濤,高照良,田紅衛(wèi). 工程堆積體陡坡坡面土壤侵蝕水動力學過程[J]. 農(nóng)業(yè)工程學報,2013,29(24):94-102.

    Zhang Letao, Gao Zhaoliang, Tian Hongwei. Hydrodynamic process of soil erosion in steep slope of engineering accumulation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(24): 94-102. (in Chinese with English abstract)

    [8] 牛耀彬,高照良,李永紅,等. 工程堆積體坡面細溝形態(tài)發(fā)育及其與產(chǎn)流產(chǎn)沙量的關系[J]. 農(nóng)業(yè)工程學報,2016,32(19):154-161.

    Niu Yaobin, Gao Zhaoliang, Li Yonghong, et al. Rill morphology development of engineering accumulation and its relationship with runoff and sediment[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(19): 154-161. (in Chinese with English abstract)

    [9] 楊婭雙,王金滿,萬德鵬. 人工堆墊地貌微地形改造及其水土保持效果研究進展[J]. 生態(tài)學雜志,2018,37(2):569-579.

    Yang Yashuang, Wang Jinman, Wan Depeng. Mcro- topography modification and its effects on the conservation of soil and water in artificially piled landform area: A review[J]. Chinese Journal of Ecology, 2018, 37(2): 569-579. (in Chinese with English abstract)

    [10] 郭建英,何京麗,李錦榮,等. 典型草原大型露天煤礦排土場邊坡水蝕控制效果[J]. 農(nóng)業(yè)工程學報,2015,31(3):296-303.

    Guo Jianying, He Jingli, Li Jinrong, et al. Effects of different measures on water erosion control of dump slope at opencast coal mine in typical steppe[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(3): 296-303. (in Chinese with English abstract)

    [11] Feng Tianjiao, Wei Wei, Chen Liding, et al. Effects of land preparation and artificial vegetation on soil moisture in a hilly loess catchment in China[J]. Land Degradation & Development, 2018, 29: 1427-1441.

    [12] Novianti V, Marrs R H, Choesin D N, et al. Natural regeneration on land degraded by coal mining in a tropical climate: Lessons for ecological restoration from Indonesia[J]. Land Degradation & Development, 2018, 29: 4050-4060.

    [13] Brofas O, Varelides C. Hydro-seeding and mulching for establishing vegetation on mining spoils in Greece[J]. Land Degradation & Development, 2000, 11: 375-382.

    [14] Tizado E J, Nú?ez-Pérez E. Terrestrial arthropods in the initial restoration stages of anthracite coal mine spoil heaps in northwestern Spain: potential usefulness of higher taxa as restoration indicators[J]. Land Degradation & Development, 2016, 27: 1131-1140.

    [15] Harantová L, Mudrák O, Kohout P, et al. Development of microbial community during primary succession in areas degraded by mining activities[J]. Land Degradation & Development, 2017, 28: 2574-2584.

    [16] 李萍,朱清科,趙磊磊,等. 黃土丘陵溝壑區(qū)魚鱗坑雨季土壤水分狀況[J]. 農(nóng)業(yè)工程學報,2011,27(7):76-81.

    Li Ping, Zhu Qingke, Zhao Leilei, et al. Soil moisture of fish-scale pit during rainy season in Loess hilly and gully region[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(7): 76-81. (in Chinese with English abstract)

    [17] 周瑤,馬紅彬,賈希洋,等. 不同生態(tài)恢復措施下寧夏黃土丘陵典型草原土壤質(zhì)量評價[J]. 農(nóng)業(yè)工程學報,2017,33(18):102-110.

    Zhou Yao, Ma Hongbin, Jia Xiyang, et al. Soil quality assessment under different ecological restoration measures in typical steppe in loess hilly area in Ningxia[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(18): 102-110. (in Chinese with English abstract)

    [18] 張海,張立新,柏延芳,等. 黃土峁狀丘陵區(qū)坡地治理模式對土壤水分環(huán)境及植被恢復效應[J]. 農(nóng)業(yè)工程學報,2007,23(11):108-113.

    Zhang Hai, Zhang Lixin, Bai Yanfang, et al.Effects of management models in sloping fields on soil moisture and vegetation restoration in the hilly and gully regions of the Loess Plateau[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2007, 23(11): 108-113. (in Chinese with English abstract)

    [19] 張北贏,徐學選,劉文兆. 黃土丘陵溝壑區(qū)不同水保措施條件下土壤水分狀況[J]. 農(nóng)業(yè)工程學報,2009,25(4):54-58.

    Zhang Beiying, Xu Xuexuan, Liu Wenzhao. Soil water condition under different measures of soil and water conservation in loess hilly and gully region[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(4): 54-58. (in Chinese with English abstract)

    [20] 任柯蒙,衛(wèi)偉,趙西寧,等. 基于水蝕預報模型的黃土高原水平階減流阻蝕效應模擬[J]. 生態(tài)學報,2018,38(14):5067-5077.

    Ren Kemeng, Wei Wei, Zhao Xining, et al. Simulation of the effects of level bench terraces on water erosion reduction in the Loess Plateau based on the WEPP model[J]. Acta Ecologica, 2018, 38(14): 5067-5077. (in Chinese with English abstract)

    [21] 張巖,楊松,李鎮(zhèn),等. 陜北黃土區(qū)水平條帶整地措施對切溝發(fā)育的影響[J]. 農(nóng)業(yè)工程學報,2015,31(7):125-130.

    Zhang Yan, Yang Song, Li Zhen, et al. Effect of narrow terrace on gully erosion in Northern Shaanxi Loess area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(7): 125-130. (in Chinese with English abstract)

    [22] Li Lusheng, Gao Xiaodong, Wu Pute, et al. Soil water content and root patterns in a rain-fed jujube plantation across stand ages on the loess plateau of China[J]. Land Degradation & Development, 2017, 28: 207-216.

    [23] 楊波,王文龍,郭明明,等. 模擬降雨條件下棄渣體邊坡不同防護措施的減水減沙效益[J]. 土壤學報,2017,54(6):1357-1368.

    Yang Bo, Wang Wenlong, Guo Mingming, et al. Study on runoff and sediment reducing effects of erosion control measures on slopes of slag dumps with simulated rainfall[J]. Acta Pedologica Sinica, 2017, 54(6): 1357-1368. (in Chinese with English abstract)

    [24] Gilley J E, Gee G W, Bauer A, et al. Runoff and erosion characteristics of surface-mined sites in western North Dakota[J]. Transactions of the ASAE, 1977, 20(4): 697-700, 704.

    [25] Martín-Moreno C, Martín Duque J F M, Nicolau Ibarra J M, et al. Effects of topography and surface soil cover on erosion for mining reclamation: The experimental spoil heap at EI Machorro mine (Central Spain)[J]. Land Degradation & Development, 2016, 27: 145-159.

    [26] 府谷縣志編纂委員會. 府谷縣志[M]. 西安:陜西人民出版社,1994.

    [27] 張樂濤,高照良,李永紅,等. 模擬徑流條件下工程堆積體陡坡土壤侵蝕過程[J]. 農(nóng)業(yè)工程學報,2013,29(8):145-153.

    Zhang Letao, Gao Zhaoliang, Li Yonghong, et al. Soil erosion process of engineering accumulation in steep slope under simulated runoff conditions[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(8): 145-153. (in Chinese with English abstract)

    [28] Zhang Letao, Li Zhanbin, Wang He, et al. Influence of intra-event-based flood regime on sediment flow behavior from a typical agro-catchment of the Chinese Loess Plateau[J]. Journal of Hydrology, 2016, 538: 71-81.

    [29] Zhang Jianjun, Zhang Xiaoping, Li Rui, et al. Did streamflow or suspended sediment concentration changes reduce sediment load in the middle reaches of the Yellow River?[J]. Journal of Hydrology, 2017, 546: 357-369.

    [30] Zhang Jianjun, Gao Guangyao, Fu Bojie, et al. Explanation of climate and human impacts on sediment discharge change in Darwinian hydrology: Derivation of a differential equation[J]. Journal of Hydrology, 2018, 559: 827-834.

    Sediment-reducing benefits by runoff regulation under engineering measures in steep slope of abandoned soil deposits in Chinese loessial region

    Zhang Letao1, Dong Junwu1, Yuan Lin1, Zhu Yaqin1, Huang Zengyu1, Li Wanlu1, Wang Huihui1, Tang Linhao1, Tian Hongwei2, Gao Zhaoliang3※

    (1.475004,;2.,,430010,; 3.,,712100,)

    It is of significant sense to quantitatively distinguish the different role of decreasing runoff and changing flow-sediment relationship in sediment control, which will be helpful to profoundly understand the soil and water conservation benefits of runoff regulation measures. In order to explore the water-sediment effects and associated control efficiency of different engineering measures and their combinations in regulating runoff yield and sediment export from the deposit slope, a series of simulated runoff scouring experiments under the control of different engineering measures were conducted on a steep deposit slope (36°) along the Shenmu-Fugu Expressway in the loessial region of China. Several types of traditional engineering measures were designed in the study, among which individual measures included level terraces, level trenches, and fish-scale pits. In order to optimize the layouts of engineering measures on the slope, combined measures of level terraces + level trenches, level trenches + fish-scale pits, and level trenches+ fish-scale pits were also designed. Closed runoff plots with a size of 2.5×12 m2were established on the deposit slope, clean water was added at the top of the plot with an inflow rate of 25 L/min to simulate surface runoff processes. Runoff and sediment was funneled through a PVC tube down to a bucket with a radius of 0.49 m and a height of 1 m. Runoff and sediment samples from the plots were collected by were taken using1000 mL plastic bottles every 1 min in initial 3 min of each run and then every 2 min for the rest of the duration, for determining sediment concentration. Therefore, the runoff and sediment process was calculated under different scenarios. The results showed that slope runoff and erosion were well controlled by engineering measures, and the control ratio of runoff commencement time, runoff, sediment export, and mean sediment concentration varied from 2 to 20, 0.45 to 0.78, 0.20 to 0.59, and 0.38 to 0.79 under different situations, respectively. Overall, sediment reduction induced by decreasing runoff was the main cause for sediment reduction benefits by engineering measures. However, sediment reduction induced by flow-sediment relationship changes was highly dependent on the behavior of decreasing runoff in reducing sediment output. For level terrace-based measures, sediment reduction induced by flow-sediment relationship changes was positively and linearly correlated to sediment reduction induced by decreasing runoff. The sediment regulation effects of altering flow-sediment relationship initiated, when sediment reduction induced by decreasing runoff exceeded a certain threshold. For level trench-based measures, there was a quadratic function relationship between sediment reduction induced by flow-sediment relationship changes and sediment reduction induced by decreasing runoff. Furthermore, a climax existed for sediment reduction induced by flow-sediment changes. In terms of sediment control, the level trench-based measures were more efficient than level terrace-based measures. In the case of the combination of level trench and fish-scale pit, the potential of decreasing runoff and altering flow-sediment relationship in reducing sediment output was well exerted with the efficiency coefficient being 55%. It indicated that relative equilibrium between the two causes for sediment reduction was kept at a high level under this combination. Therefore, the application of different engineering measures and their combinations should be considered based on specific soil and water conservation benefits and erosion control goals. The study provides theoretical references for the optimization of engineering measure in soil and water conservation management of steep deposit slope.

    erosion; sediment; runoff; flow-sediment relationship; sediment-reducing capacity; efficiency coefficient; abandoned soil deposits

    10.11975/j.issn.1002-6819.2019.15.014

    S157.1

    A

    1002-6819(2019)-15-0101-09

    2019-01-26

    2019-07-28

    國家自然科學基金項目(41807066,41671283,51309093);國家重點研發(fā)計劃項目(2016YFC0501706-02);河南大學一流學科創(chuàng)新團隊培育項目(2018YLTD16)

    張樂濤,講師,博士,研究方向為土壤侵蝕與水土保持。Email:letao_20062006@126.com

    高照良,研究員,博士生導師,研究方向為土壤侵蝕與水土保持。Email:gzl@ms.iswc.ac.cn

    張樂濤,董俊武,袁琳,朱雅琴,黃增玉,李婉璐,王慧慧,唐林昊,田紅衛(wèi),高照良. 黃土區(qū)工程堆積體陡坡坡面徑流調(diào)控工程措施的減沙效應[J]. 農(nóng)業(yè)工程學報,2019,35(15):101-109. doi:10.11975/j.issn.1002-6819.2019.15.014 http://www.tcsae.org

    Zhang Letao, Dong Junwu, Yuan Lin, Zhu Yaqin, Huang Zengyu, Li Wanlu, Wang Huihui, Tang Linhao, Tian Hongwei, Gao Zhaoliang. Sediment-reducing benefits by runoff regulation under engineering measures in steep slope of abandoned soil deposits in Chinese loessial region[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(15): 101-109. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.15.014 http://www.tcsae.org

    猜你喜歡
    調(diào)沙沙量水沙
    鄱陽湖入江水道輸沙能力與出湖沙量關系及沖淤調(diào)整規(guī)律
    水利學報(2023年8期)2023-09-19 02:49:08
    迎戰(zhàn)調(diào)水調(diào)沙
    1996-2020年黃河口清水溝流路海岸線動態(tài)演變及其與水沙量的關系研究
    海洋通報(2022年4期)2022-10-10 07:41:22
    2011—2020年調(diào)水調(diào)沙前后黃河口海域浮游動物年間變化及影響因子
    海洋科學(2022年12期)2022-03-01 07:25:40
    黃河寧夏段主要支流水沙量變化分析
    大型水利樞紐下游水沙變異特征
    山區(qū)河流上下雙丁壩回流區(qū)水沙特性淺探
    江西建材(2018年1期)2018-04-04 05:26:28
    潼關水文站近十年水沙情勢分析
    長江上中下游河道水沙特征和水沙關系
    黃河調(diào)水調(diào)沙對濮陽河段的影響分析
    av免费观看日本| 中文字幕免费在线视频6| 夜夜爽夜夜爽视频| 99视频精品全部免费 在线| av国产免费在线观看| 美女xxoo啪啪120秒动态图| av女优亚洲男人天堂| 国产色爽女视频免费观看| 少妇高潮的动态图| 亚洲国产av新网站| 久久午夜福利片| 狠狠精品人妻久久久久久综合| 熟女av电影| 高清黄色对白视频在线免费看 | 久久精品久久久久久噜噜老黄| 夫妻午夜视频| a级毛色黄片| av专区在线播放| 三级国产精品欧美在线观看| 少妇的逼好多水| 日韩亚洲欧美综合| 国产淫语在线视频| 日韩av不卡免费在线播放| 国产成人精品婷婷| 99久国产av精品国产电影| 国产日韩欧美亚洲二区| 国产免费福利视频在线观看| 国产男人的电影天堂91| 在线免费十八禁| 久久久a久久爽久久v久久| 黑丝袜美女国产一区| 26uuu在线亚洲综合色| 国产黄片视频在线免费观看| 中文字幕免费在线视频6| 免费大片黄手机在线观看| 欧美精品一区二区免费开放| 校园人妻丝袜中文字幕| 亚洲人成网站在线播| 亚洲国产日韩一区二区| 日韩不卡一区二区三区视频在线| 久久精品久久精品一区二区三区| 美女视频免费永久观看网站| 免费看光身美女| 日本av免费视频播放| 一级二级三级毛片免费看| 国产乱来视频区| 欧美老熟妇乱子伦牲交| 国产精品成人在线| 日日摸夜夜添夜夜爱| 亚洲av不卡在线观看| 国产在线视频一区二区| 男女下面进入的视频免费午夜| 久久久久久久精品精品| 丝瓜视频免费看黄片| 成人免费观看视频高清| 十分钟在线观看高清视频www | 亚洲人与动物交配视频| 久久久精品免费免费高清| 91久久精品电影网| 亚洲欧洲日产国产| 干丝袜人妻中文字幕| 色5月婷婷丁香| 国产黄频视频在线观看| 91精品国产国语对白视频| 国产精品蜜桃在线观看| 老司机影院毛片| 国产成人91sexporn| 国模一区二区三区四区视频| 日韩电影二区| 国产黄色免费在线视频| 青春草亚洲视频在线观看| 国产中年淑女户外野战色| 少妇的逼水好多| 91午夜精品亚洲一区二区三区| 欧美少妇被猛烈插入视频| 精品人妻熟女av久视频| 91狼人影院| 国产爱豆传媒在线观看| 成人黄色视频免费在线看| 又大又黄又爽视频免费| 一本色道久久久久久精品综合| www.色视频.com| xxx大片免费视频| 国产一区亚洲一区在线观看| 草草在线视频免费看| 国产免费又黄又爽又色| 国产成人午夜福利电影在线观看| 麻豆成人av视频| 亚洲经典国产精华液单| 亚洲精品乱码久久久久久按摩| 亚洲av成人精品一二三区| 国国产精品蜜臀av免费| 女性被躁到高潮视频| 亚洲精品国产av蜜桃| 国产在线一区二区三区精| 久久久精品免费免费高清| 欧美+日韩+精品| 亚洲成人中文字幕在线播放| 好男人视频免费观看在线| av在线老鸭窝| 国产高潮美女av| 多毛熟女@视频| 久久久久性生活片| 国产白丝娇喘喷水9色精品| 国产伦精品一区二区三区视频9| 18禁裸乳无遮挡免费网站照片| 日本欧美视频一区| 欧美日韩在线观看h| 一二三四中文在线观看免费高清| 国产在线一区二区三区精| 免费观看av网站的网址| 青春草亚洲视频在线观看| 久热久热在线精品观看| 日本免费在线观看一区| 成人18禁高潮啪啪吃奶动态图 | 亚洲av.av天堂| 亚洲四区av| 黑人猛操日本美女一级片| 韩国av在线不卡| 中国美白少妇内射xxxbb| 成人黄色视频免费在线看| 国产成人a∨麻豆精品| 免费黄频网站在线观看国产| 日产精品乱码卡一卡2卡三| 黄色怎么调成土黄色| 久久久久久人妻| 插逼视频在线观看| 狂野欧美激情性xxxx在线观看| 免费看av在线观看网站| 亚洲电影在线观看av| 高清av免费在线| 青春草国产在线视频| 久久久久国产网址| 在线观看三级黄色| 精品人妻熟女av久视频| 成人亚洲精品一区在线观看 | 国产欧美日韩一区二区三区在线 | 看免费成人av毛片| 成人国产av品久久久| 狂野欧美激情性xxxx在线观看| 国产一级毛片在线| 人妻夜夜爽99麻豆av| 狂野欧美激情性xxxx在线观看| 欧美成人精品欧美一级黄| 成人亚洲欧美一区二区av| 亚洲国产欧美在线一区| 特大巨黑吊av在线直播| 国产免费一级a男人的天堂| 成人影院久久| 老熟女久久久| 这个男人来自地球电影免费观看 | 久久ye,这里只有精品| 美女视频免费永久观看网站| 色网站视频免费| 一级a做视频免费观看| 干丝袜人妻中文字幕| 久久国产精品大桥未久av | 国产美女午夜福利| 成人亚洲精品一区在线观看 | av在线app专区| 91久久精品电影网| 啦啦啦视频在线资源免费观看| 老师上课跳d突然被开到最大视频| 91久久精品电影网| 啦啦啦中文免费视频观看日本| 精品亚洲成a人片在线观看 | 中文在线观看免费www的网站| 交换朋友夫妻互换小说| 亚洲aⅴ乱码一区二区在线播放| 国产日韩欧美在线精品| 日本色播在线视频| 丝袜脚勾引网站| 欧美+日韩+精品| 国产精品一及| 少妇人妻精品综合一区二区| 高清午夜精品一区二区三区| 最近中文字幕高清免费大全6| 在线观看人妻少妇| 精品久久久久久久久亚洲| 欧美日韩国产mv在线观看视频 | 啦啦啦视频在线资源免费观看| 在线观看一区二区三区激情| 日韩伦理黄色片| 色婷婷久久久亚洲欧美| 国产在线免费精品| 国产精品爽爽va在线观看网站| 精品99又大又爽又粗少妇毛片| 男女免费视频国产| 啦啦啦视频在线资源免费观看| 精品一区二区三卡| 精华霜和精华液先用哪个| 亚洲四区av| 久久久久久久久久成人| 日本猛色少妇xxxxx猛交久久| 亚洲精品456在线播放app| 国产熟女欧美一区二区| 成人18禁高潮啪啪吃奶动态图 | 毛片女人毛片| 久久久久久久久久久免费av| 欧美性感艳星| av专区在线播放| 又大又黄又爽视频免费| 日本黄色片子视频| 国产精品一区二区性色av| 日本爱情动作片www.在线观看| 国产大屁股一区二区在线视频| 免费高清在线观看视频在线观看| 日本欧美国产在线视频| 网址你懂的国产日韩在线| 精品久久久久久久久av| 极品少妇高潮喷水抽搐| 汤姆久久久久久久影院中文字幕| 中文乱码字字幕精品一区二区三区| 99re6热这里在线精品视频| xxx大片免费视频| 一级av片app| 欧美xxⅹ黑人| 国产黄片视频在线免费观看| 精品视频人人做人人爽| 视频区图区小说| 狂野欧美激情性bbbbbb| 日本黄大片高清| 91在线精品国自产拍蜜月| 成人一区二区视频在线观看| 日本vs欧美在线观看视频 | 伦精品一区二区三区| 日本wwww免费看| 久久ye,这里只有精品| 大话2 男鬼变身卡| 亚洲欧美日韩无卡精品| 国产视频首页在线观看| 日本欧美国产在线视频| videos熟女内射| .国产精品久久| 日本vs欧美在线观看视频 | 激情 狠狠 欧美| 亚洲av综合色区一区| 国产探花极品一区二区| 3wmmmm亚洲av在线观看| 久久久久久久精品精品| 欧美 日韩 精品 国产| 国产 一区 欧美 日韩| 亚洲精品一二三| 久久热精品热| 国产精品一区二区在线不卡| 免费看日本二区| www.av在线官网国产| 欧美bdsm另类| 男人和女人高潮做爰伦理| 国产一级毛片在线| 最新中文字幕久久久久| 亚洲国产欧美在线一区| 18禁在线播放成人免费| 一区二区三区免费毛片| 少妇 在线观看| 高清在线视频一区二区三区| 久久精品人妻少妇| 亚洲av电影在线观看一区二区三区| 91在线精品国自产拍蜜月| 久久ye,这里只有精品| 又粗又硬又长又爽又黄的视频| 日韩伦理黄色片| 街头女战士在线观看网站| 国产乱人视频| 日本免费在线观看一区| 精品酒店卫生间| 国产伦理片在线播放av一区| 亚洲三级黄色毛片| 老熟女久久久| 妹子高潮喷水视频| 在线免费十八禁| 免费观看的影片在线观看| 久久 成人 亚洲| 免费大片18禁| 欧美国产精品一级二级三级 | 久久热精品热| 国产免费福利视频在线观看| 国产精品一区二区在线观看99| 亚洲性久久影院| 日韩欧美一区视频在线观看 | 国产人妻一区二区三区在| 夜夜骑夜夜射夜夜干| 亚洲,一卡二卡三卡| 在线观看美女被高潮喷水网站| 麻豆国产97在线/欧美| 一级毛片 在线播放| 国产精品秋霞免费鲁丝片| 精华霜和精华液先用哪个| 国产亚洲5aaaaa淫片| av网站免费在线观看视频| 亚洲电影在线观看av| 大陆偷拍与自拍| 国产男女超爽视频在线观看| 免费大片黄手机在线观看| 91狼人影院| 亚洲国产精品999| 蜜桃在线观看..| 国产在线男女| 成人一区二区视频在线观看| 国产精品麻豆人妻色哟哟久久| 狠狠精品人妻久久久久久综合| 婷婷色av中文字幕| 干丝袜人妻中文字幕| 天堂俺去俺来也www色官网| 一级黄片播放器| 狂野欧美白嫩少妇大欣赏| av女优亚洲男人天堂| 成年人午夜在线观看视频| 高清欧美精品videossex| 日本一二三区视频观看| 最后的刺客免费高清国语| av国产久精品久网站免费入址| 在线观看免费高清a一片| 国产成人免费无遮挡视频| 国产成人免费无遮挡视频| av国产精品久久久久影院| 国产高清有码在线观看视频| 免费观看性生交大片5| 一级毛片久久久久久久久女| 欧美一区二区亚洲| 青春草亚洲视频在线观看| 成人漫画全彩无遮挡| 天堂8中文在线网| 观看美女的网站| 一区二区三区四区激情视频| 国产av国产精品国产| 日日摸夜夜添夜夜爱| 天天躁日日操中文字幕| 男男h啪啪无遮挡| 99久久综合免费| 成人无遮挡网站| 日韩免费高清中文字幕av| 国产亚洲精品久久久com| 大香蕉97超碰在线| 少妇人妻精品综合一区二区| 人人妻人人澡人人爽人人夜夜| 欧美另类一区| 久久久久久久久久久丰满| 日韩一区二区视频免费看| 午夜视频国产福利| 中国美白少妇内射xxxbb| 国产探花极品一区二区| 国产精品.久久久| av天堂中文字幕网| 亚洲第一区二区三区不卡| 日韩欧美 国产精品| 人妻夜夜爽99麻豆av| 亚洲精品乱码久久久久久按摩| 只有这里有精品99| 精品久久久久久久久亚洲| 午夜日本视频在线| 久久99精品国语久久久| 欧美日韩精品成人综合77777| 国产伦精品一区二区三区视频9| 高清日韩中文字幕在线| 啦啦啦在线观看免费高清www| a级毛片免费高清观看在线播放| 亚洲av电影在线观看一区二区三区| 国产黄片视频在线免费观看| 好男人视频免费观看在线| 在现免费观看毛片| 亚洲人与动物交配视频| 免费看日本二区| 中文字幕人妻熟人妻熟丝袜美| 婷婷色麻豆天堂久久| av国产久精品久网站免费入址| 赤兔流量卡办理| 少妇人妻 视频| 成人漫画全彩无遮挡| 卡戴珊不雅视频在线播放| 99久国产av精品国产电影| 免费在线观看成人毛片| 内射极品少妇av片p| 日韩三级伦理在线观看| 国产精品久久久久成人av| 最近最新中文字幕大全电影3| 你懂的网址亚洲精品在线观看| 七月丁香在线播放| 亚洲综合精品二区| 在线观看一区二区三区激情| 亚洲美女黄色视频免费看| 国产精品人妻久久久影院| 久久99热这里只频精品6学生| av不卡在线播放| 国产无遮挡羞羞视频在线观看| 乱码一卡2卡4卡精品| 欧美老熟妇乱子伦牲交| 国产精品一区二区在线观看99| 中文字幕av成人在线电影| 少妇被粗大猛烈的视频| 国产亚洲91精品色在线| 久久青草综合色| 久久久久网色| 亚洲精品国产av蜜桃| 午夜免费观看性视频| h日本视频在线播放| 国产精品一区二区在线不卡| 久久久a久久爽久久v久久| 黄色日韩在线| 日韩av免费高清视频| 欧美3d第一页| 欧美人与善性xxx| 观看av在线不卡| 久久精品久久精品一区二区三区| 亚洲成人一二三区av| 国产精品一区二区三区四区免费观看| 日韩精品有码人妻一区| 成人无遮挡网站| 国产久久久一区二区三区| 日韩成人伦理影院| 日本-黄色视频高清免费观看| 五月伊人婷婷丁香| av一本久久久久| 最近手机中文字幕大全| 蜜桃亚洲精品一区二区三区| 成人一区二区视频在线观看| 久久女婷五月综合色啪小说| 国产精品偷伦视频观看了| 大片电影免费在线观看免费| 欧美zozozo另类| 少妇熟女欧美另类| 在线观看一区二区三区激情| 黄色一级大片看看| 亚洲av中文av极速乱| 精品午夜福利在线看| 91精品伊人久久大香线蕉| 日韩一区二区三区影片| tube8黄色片| 波野结衣二区三区在线| 91aial.com中文字幕在线观看| 国产男女内射视频| 中文字幕精品免费在线观看视频 | 精品少妇久久久久久888优播| 国产淫片久久久久久久久| 国产日韩欧美亚洲二区| www.av在线官网国产| 在线观看免费日韩欧美大片 | 在线免费观看不下载黄p国产| 亚洲人成网站高清观看| 人人妻人人澡人人爽人人夜夜| 91aial.com中文字幕在线观看| 欧美国产精品一级二级三级 | 亚洲av电影在线观看一区二区三区| 老司机影院成人| 国产大屁股一区二区在线视频| 亚洲中文av在线| 久久国产乱子免费精品| 青春草亚洲视频在线观看| 熟女人妻精品中文字幕| 亚洲av成人精品一区久久| 国产精品一区二区在线观看99| 免费观看a级毛片全部| 偷拍熟女少妇极品色| 在线免费观看不下载黄p国产| 日韩电影二区| 国产免费一级a男人的天堂| 一本色道久久久久久精品综合| 卡戴珊不雅视频在线播放| 欧美成人精品欧美一级黄| 亚洲,欧美,日韩| 亚洲精品,欧美精品| 日韩强制内射视频| 久久久久国产精品人妻一区二区| 亚洲丝袜综合中文字幕| 欧美高清成人免费视频www| 三级国产精品欧美在线观看| 免费不卡的大黄色大毛片视频在线观看| 精品久久久久久久末码| 国产久久久一区二区三区| 少妇高潮的动态图| 欧美bdsm另类| 午夜福利视频精品| 国产爽快片一区二区三区| 成人高潮视频无遮挡免费网站| 亚洲熟女精品中文字幕| 日本av手机在线免费观看| 国产一级毛片在线| 亚洲欧美日韩东京热| 国产午夜精品久久久久久一区二区三区| 激情 狠狠 欧美| 建设人人有责人人尽责人人享有的 | 国产精品熟女久久久久浪| 一级爰片在线观看| 男女免费视频国产| 国产亚洲欧美精品永久| 99久久精品国产国产毛片| 少妇精品久久久久久久| 欧美日韩一区二区视频在线观看视频在线| 国产成人免费观看mmmm| 最近手机中文字幕大全| 国产成人精品一,二区| 国产永久视频网站| 久久久精品94久久精品| 麻豆乱淫一区二区| 久久亚洲国产成人精品v| 国产精品福利在线免费观看| 观看美女的网站| 狠狠精品人妻久久久久久综合| 国产精品一区二区在线观看99| 色吧在线观看| 亚洲精品日本国产第一区| 99久久精品一区二区三区| 汤姆久久久久久久影院中文字幕| 超碰av人人做人人爽久久| 女的被弄到高潮叫床怎么办| h日本视频在线播放| 一本色道久久久久久精品综合| 亚洲精品日本国产第一区| a级一级毛片免费在线观看| 汤姆久久久久久久影院中文字幕| 日本av手机在线免费观看| 高清毛片免费看| 亚洲国产高清在线一区二区三| 免费久久久久久久精品成人欧美视频 | 精品99又大又爽又粗少妇毛片| 国产av精品麻豆| 免费久久久久久久精品成人欧美视频 | 乱码一卡2卡4卡精品| 纵有疾风起免费观看全集完整版| 成人影院久久| 成年av动漫网址| 久久热精品热| 亚洲国产色片| 国产爱豆传媒在线观看| 婷婷色综合www| 国产精品.久久久| 天堂8中文在线网| 80岁老熟妇乱子伦牲交| 日韩大片免费观看网站| 欧美成人精品欧美一级黄| 日韩一本色道免费dvd| 中文乱码字字幕精品一区二区三区| 91精品一卡2卡3卡4卡| 亚洲国产av新网站| 51国产日韩欧美| 性高湖久久久久久久久免费观看| 国产亚洲91精品色在线| 亚洲av综合色区一区| av.在线天堂| 免费大片18禁| 一个人免费看片子| 免费高清在线观看视频在线观看| 亚洲成人手机| 免费少妇av软件| 国产亚洲5aaaaa淫片| 18禁在线无遮挡免费观看视频| 秋霞伦理黄片| 日韩不卡一区二区三区视频在线| 国产淫片久久久久久久久| 日韩,欧美,国产一区二区三区| 亚洲一级一片aⅴ在线观看| 国产亚洲一区二区精品| 国产有黄有色有爽视频| 日韩免费高清中文字幕av| 永久网站在线| 在现免费观看毛片| av.在线天堂| 亚洲va在线va天堂va国产| 人人妻人人澡人人爽人人夜夜| 国产av国产精品国产| 在线精品无人区一区二区三 | 精品视频人人做人人爽| 国产欧美另类精品又又久久亚洲欧美| 精品人妻一区二区三区麻豆| 久久99热这里只有精品18| 成年av动漫网址| 汤姆久久久久久久影院中文字幕| 亚洲婷婷狠狠爱综合网| 尾随美女入室| 一级片'在线观看视频| 亚洲图色成人| 中文欧美无线码| 色婷婷久久久亚洲欧美| 色综合色国产| 在线观看国产h片| 又黄又爽又刺激的免费视频.| 观看av在线不卡| 99热6这里只有精品| 亚洲国产高清在线一区二区三| 国产黄片视频在线免费观看| 国产伦精品一区二区三区四那| 国产成人91sexporn| 美女高潮的动态| 国内少妇人妻偷人精品xxx网站| 亚洲婷婷狠狠爱综合网| 欧美人与善性xxx| 欧美国产精品一级二级三级 | 国产黄色视频一区二区在线观看| 欧美区成人在线视频| 女人十人毛片免费观看3o分钟| 国产午夜精品一二区理论片| 亚洲国产欧美在线一区| 亚洲国产精品国产精品| 观看av在线不卡| 韩国高清视频一区二区三区| 国产免费一级a男人的天堂| videossex国产| 晚上一个人看的免费电影| 汤姆久久久久久久影院中文字幕| 日韩亚洲欧美综合| 久久国产精品大桥未久av | 成人国产av品久久久| 最近的中文字幕免费完整| 成年美女黄网站色视频大全免费 | 99国产精品免费福利视频| 丰满少妇做爰视频| 久久热精品热| 欧美日韩亚洲高清精品| 高清日韩中文字幕在线| 亚洲av在线观看美女高潮| 26uuu在线亚洲综合色| 免费看av在线观看网站| 日日啪夜夜撸| 久久久久国产精品人妻一区二区| 亚洲av日韩在线播放|