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

    模擬降雨條件下棄渣體邊坡不同防護措施的減水減沙效益*

    2017-12-11 05:32:43王文龍郭明明康宏亮李建明
    土壤學(xué)報 2017年6期
    關(guān)鍵詞:渣體魚鱗坑堆積體

    楊 波 王文龍,2,? 郭明明 康宏亮 白 蕓 李建明

    (1 西北農(nóng)林科技大學(xué)水土保持研究所黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室,陜西楊凌 712100)

    (2 中國科學(xué)院水利部水土保持研究所,陜西楊凌 712100)

    (3 榆林學(xué)院陜西省陜北礦區(qū)生態(tài)修復(fù)重點實驗室,陜西榆林 719000)

    (4 長江科學(xué)院水土保持研究所,武漢 430010)

    模擬降雨條件下棄渣體邊坡不同防護措施的減水減沙效益*

    楊 波1王文龍1,2,3?郭明明1康宏亮1白 蕓3李建明4

    (1 西北農(nóng)林科技大學(xué)水土保持研究所黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室,陜西楊凌 712100)

    (2 中國科學(xué)院水利部水土保持研究所,陜西楊凌 712100)

    (3 榆林學(xué)院陜西省陜北礦區(qū)生態(tài)修復(fù)重點實驗室,陜西榆林 719000)

    (4 長江科學(xué)院水土保持研究所,武漢 430010)

    采用野外人工模擬降雨方法,以未防護棄渣體邊坡為對照,研究了1.0、1.5 mm min-1降雨強度條件下神府礦區(qū)種草和魚鱗坑措施對偏土質(zhì)、偏石質(zhì)和煤矸石棄渣體產(chǎn)流產(chǎn)沙的調(diào)控作用。結(jié)果表明:(1)不同措施下3種棄渣體邊坡徑流率均在產(chǎn)流6~9 min后趨于穩(wěn)定,產(chǎn)流過程中棄渣體邊坡侵蝕速率均呈波動減小趨勢,且未防護坡面減小趨勢較防護更明顯。(2)植草對偏土質(zhì)、偏石質(zhì)、煤矸石棄渣體減水和減沙效益分別為42.91%~51.21%、26.28%~55.20%、10.33%和97.54%~97.95%、41.87%~42.26%、7.80%;魚鱗坑的減水和減沙效益則分別為:51.89%~72.72%、22.37%~42.92%、21.32%和98.41% ~99.30%、94.90%~91.84%、39.50%。(3)魚鱗坑措施對偏土質(zhì)棄渣體的減水和減沙效益較種草分別提高8.98%~21.51%和0.46%~1.76%;種草措施對偏石質(zhì)和煤矸石棄渣體的減水和減沙效益較魚鱗坑措施分別提高3.91%~12.28%和0.28%~3.06%及10.99%和31.70%。(4)3種未防護棄渣體的侵蝕速率和徑流率呈顯著線性關(guān)系,種草和魚鱗坑措施改變了坡面水沙關(guān)系,侵蝕速率和徑流率相關(guān)性減弱或無顯著關(guān)系。研究結(jié)果可為礦區(qū)棄渣體邊坡生態(tài)恢復(fù)措施布設(shè)提供科學(xué)指導(dǎo)。

    神府煤田;棄渣體;魚鱗坑;植被措施;徑流率;侵蝕速率

    神府礦區(qū)煤炭資源富集,露天開采是該區(qū)主要的煤礦開采方式。資料顯示[1],露天煤礦每采萬噸煤約破壞土地面積0.l hm2,約排放剝離物20 000~61 000 m3。大面積的開采和排廢不僅破壞和擾動原始地面,還產(chǎn)生大量的棄土棄渣體,引發(fā)了劇烈的新增水土流失,導(dǎo)致礦區(qū)生態(tài)環(huán)境遭到嚴(yán)重破壞。

    當(dāng)前礦區(qū)生態(tài)環(huán)境問題已經(jīng)得到了社會的廣泛關(guān)注和重視[2-4],礦區(qū)水土流失控制已成為解決當(dāng)?shù)厣鷳B(tài)環(huán)境問題的重點內(nèi)容之一。棄土棄渣等工程堆積體類型繁多、形態(tài)各異[5-6],且具有物質(zhì)組成復(fù)雜,堆積坡度大,結(jié)構(gòu)疏松,穩(wěn)定性差等特點,堆積體邊坡坡面可蝕性是破壞前土壤的10倍~100倍[7],Zhou等[8]研究表明露天煤礦復(fù)墾土地有機碳、氮含量需要20a時間才可恢復(fù)至開采前的水平。因此諸多學(xué)者針對礦區(qū)棄土棄渣體土壤侵蝕問題做了大量研究。多通過人工模擬降雨[9-14]、放水沖刷[6,15-17]以及小區(qū)定位觀測[18]的方法,集中研究了棄土棄渣堆積體產(chǎn)流產(chǎn)沙特性、水動力學(xué)參數(shù)變化以及侵蝕速率預(yù)測。郭明明等[10]認(rèn)為棄渣體坡面徑流率隨產(chǎn)流歷時呈現(xiàn)“增長-下降-穩(wěn)定”的過程。王貞等[13]研究了神東煤田擾動地面徑流特性,發(fā)現(xiàn)擾動地面坡面徑流雷諾數(shù)和阻力系數(shù)均較原始坡面高。張樂濤等[17]認(rèn)為水流功率是預(yù)測公路棄土堆積體坡面侵蝕的最佳水動力學(xué)參數(shù)。郭明明等[14]采用了中值粒徑、分形維數(shù)、降雨強度、坡度及徑流參數(shù)建立了棄土棄渣體土壤侵蝕預(yù)測模型。而在針對礦區(qū)棄土棄渣體水土流失治理措施防護方面,田野[12]研究了棄渣體邊坡措施的水土保持效益,認(rèn)為水平階護坡改變坡面物理特性,減流效果可達15.84%。劉瑞順等[18]在年限為1a的排土場邊坡通過野外徑流小區(qū)定位觀測得出,稻草簾子和沙打旺措施的減水減沙效益分別達60%和90%以上。以往的研究大都集中在礦區(qū)棄土棄渣堆積體水土流失規(guī)律及機制層面上,而在針對礦區(qū)棄渣體生態(tài)修復(fù)方面研究不夠深入,尤其是在棄土棄渣堆積初期短歷時強降雨條件下工程措施與植被措施對坡面水沙關(guān)系的研究還少有報道。

    鑒于此,本文通過模擬降雨的方法研究植被措施(種草)和工程措施(魚鱗坑)對偏土質(zhì)和偏石質(zhì)棄渣體以及煤矸石棄渣體這3種礦區(qū)典型堆積體邊坡減水減沙效益及水沙關(guān)系調(diào)控作用,探索礦區(qū)不同棄渣場最佳治理方式,為神府礦區(qū)初期的棄土棄渣堆積體能夠快速有效地進行生態(tài)恢復(fù)提供一定的參考和指導(dǎo),也為該區(qū)建立開發(fā)建設(shè)項目水土流失模型提供基礎(chǔ)數(shù)據(jù)。

    1 材料與方法

    1.1 試驗區(qū)概況

    神府東勝煤田地處晉陜蒙三省交界處(37°20'~40°16' N,108°36'~110°3' E),包括陜西的榆林市、內(nèi)蒙古的鄂爾多斯市和山西省忻州市的部分地區(qū);目前是我國已探明的最大煤田,占全國探明儲量的15%;該區(qū)位于黃土高原地區(qū)的風(fēng)蝕水蝕交錯地帶,屬于干旱半干旱氣候,年均氣溫為6.1~9.1 ℃,多年平均降雨量為325~460 mm,且多集中于7—9 月,常以暴雨形式發(fā)生,降雨量占到全年的65%~70%,年蒸發(fā)量為1636~2535 mm,自然條件惡劣,生態(tài)環(huán)境十分脆弱,該區(qū)土壤多以風(fēng)沙土和黃土為主,結(jié)構(gòu)復(fù)雜,易風(fēng)化,抗侵蝕能力差,神府東勝礦區(qū)內(nèi)的植被主要是冰草(Agropyroncristatum(Linn.)Gaertn)、沙蒿(Artemisia desertorum Spreng. Syst. Veg.)、沙柳(Salix cheilophila)、檸條(CaraganaKorshinskii Kom.)等耐旱的沙生植物,區(qū)內(nèi)植被覆蓋度很低,加之人類活動影響,使該地區(qū)生態(tài)環(huán)境十分脆弱,是我國水土流失最為嚴(yán)重的地區(qū),也是黃河泥沙的主要來源地。

    1.2 試驗因素的選取

    調(diào)查發(fā)現(xiàn),礦區(qū)最常見的堆積體是由煤礦開采過程中形成的大量棄土、棄渣、棄石及其他廢棄物的混合堆砌而成,顆粒成分以礫石、碎石、風(fēng)沙土和黃綿土為主,組成成分種類復(fù)雜,顆粒粒徑差異大,含礫石較多,疏松多孔,易發(fā)生侵蝕。試驗選取該礦區(qū)中具有代表性的偏土質(zhì)棄渣體(土石比約2∶1)、偏石質(zhì)棄渣體(土石比約1∶2)、及經(jīng)過3~5a風(fēng)化后形成的細碎石渣和土粒混合的煤矸石堆積體。前2種棄渣體顆粒組成見表1;實地調(diào)查發(fā)現(xiàn),3種堆積體坡度多集中在33°~39°,所占比例在90%以上,因此試驗小區(qū)坡度選擇為35°;根據(jù)當(dāng)?shù)囟嗄昵治g性短歷時典型暴雨特性,雨強設(shè)計為1.0 mm min-1和1.5 mm min-1,降雨歷時為24 min。

    1.3 試驗布設(shè)

    本次試驗在中國科學(xué)院水利部水土保持研究所神木侵蝕與環(huán)境試驗站進行,徑流小區(qū)為在棄渣體上安置3m×1m×1m的鋼槽,內(nèi)置棄土棄渣體后小區(qū)表面需人工平整處理,小區(qū)下端設(shè)鋼制集流槽,并在距離頂部1m和2m處設(shè)置測流斷面;棄渣體坡面防護措施為人工移植種草和布設(shè)魚鱗坑2種,移植的草為礦區(qū)周圍農(nóng)地上人工種植的兩年生冰草,挑選長勢良好的株高相近的冰草,將其植入棄渣體約20~30 cm深,從小區(qū)頂端開始每隔50 cm水平布設(shè)草帶,草帶長寬為0.9 m×0.2 m,采取移植的主要原因是為了保證在試驗時3種棄渣體坡面的植被條件一致。魚鱗坑處理是在坡面上臨時開挖魚鱗坑,坑的直徑為50 cm,深30 cm。魚鱗坑呈三角形排列,沿小區(qū)坡面等距布設(shè),坑與坑之間的縱向距離為0.5 m。移植種草小區(qū)布設(shè)完畢后,周圍設(shè)立警示標(biāo)志,以防止人為干擾破壞,每天用噴壺進行澆水養(yǎng)護并觀察長勢,2周左右,待冰草扎根較好能正常生長且每株草周圍土壤固結(jié)良好時進行試驗。開始實驗前,在小區(qū)周圍用3.0 m高的鋼管搭建降雨棚,在正上方布設(shè)下噴式模擬降雨器,降雨機型號為槽式人工模擬降雨機,噴頭高度為3.0 m,噴頭間距1.0 m,降雨雨滴終速接近天然降雨,降雨均勻系數(shù)在80%以上,并在小區(qū)上方放置2 m3儲水桶,采用60 m揚程水泵抽取儲水桶中的水,通過閥門和壓力表來控制雨強(如圖1),小區(qū)周圍用帆布圍住以減少風(fēng)對試驗的影響。本研究共設(shè)計3×1×2×3=18場降雨試驗。

    表1 土壤顆粒組成Table 1 Particle size composition of the soil

    圖1 試驗?zāi)P褪疽鈭DFig. 1 Sketch of experimental model

    1.4 試驗過程與數(shù)據(jù)采集

    試驗前用遮雨布遮蓋小區(qū),降雨器開啟后先率定雨強,率定雨強與設(shè)計值誤差控制在5%以內(nèi),均勻度在80%以上即滿足要求[19],雨強達到要求后,快速掀起小區(qū)上方的遮雨布,待坡面開始產(chǎn)流后使用秒表計時,并在集流槽出口處用徑流桶收集徑流泥沙樣,同時記錄接樣時間,前3 min內(nèi)每1 min接1次泥沙樣,3 min后每3 min接1次泥沙樣,產(chǎn)流歷時24 min,接樣同時用高錳酸鉀溶液和秒表測流速,取2個斷面的平均值作為坡面流速,將所測流速乘以0.75得到較為理想的徑流流速[20],使用精度1 mm的鋼尺測量徑流寬度和深度,同一個斷面多次測量取平均值,用普通溫度計測量水溫,采用烘干法確定泥沙重量。

    1.5 數(shù)據(jù)計算與處理

    采用Excel 2010和SPSS17.0進行數(shù)據(jù)統(tǒng)計分析,用Origin 2016 進行繪圖。

    2 結(jié)果與討論

    2.1 棄渣體產(chǎn)流特征

    2.1.1 徑流率隨時間的變化 圖2~圖4分別為1.0、1.5 mm min-1雨強下偏土質(zhì)、偏石質(zhì)和煤矸石棄渣體在未防護、種草、布設(shè)魚鱗坑措施下邊坡的徑流率隨降雨歷時的變化。由圖可知,3種棄渣體的徑流率在產(chǎn)流6~9 min后逐漸趨于穩(wěn)定。對于偏土質(zhì)棄渣體(圖2),1.0、1.5 mm min-1雨強情況下,未防護坡面徑流率隨降雨歷時呈“下降—穩(wěn)定”的變化趨勢,這與坡面結(jié)皮的產(chǎn)生和破壞過程密切相關(guān),當(dāng)一定含量的礫石鑲嵌于地表土壤時,結(jié)皮強度將大大增加,坡面入滲速率較?。?1],而當(dāng)徑流對坡面經(jīng)過一定時間的沖刷后,使得結(jié)皮破壞,入滲速率增大,并隨著土壤水分的增加逐漸趨于穩(wěn)定;種草和魚鱗坑2種防護措施下坡面徑流率隨降雨歷時則呈“增大—穩(wěn)定”的變化趨勢。2種雨強條件下,種草和魚鱗坑2種措施下坡面起始徑流率(產(chǎn)流前1 min的徑流率)較未防護坡面分別減小80.40%、90.71%和83.65%、89.65%。對于偏石質(zhì)棄渣體(圖3),2種雨強下不同措施徑流率均呈先增大后穩(wěn)定的變化趨勢。對于煤矸石堆積體(圖4),1.0 mm min-1雨強下,與防護坡面不同,在未防護坡面未觀測到產(chǎn)流現(xiàn)象,可能是由于煤矸石堆積體大孔隙較多,Wang等[22]研究表明粉煤灰和煤矸石的混合物具有與土壤相似的水力特性,尤其是具有較高的含水能力和滲透性能,加上相對雨強較小時,會導(dǎo)致降雨全部就地入滲,而在種草和布設(shè)魚鱗坑時對煤矸石棄渣體坡面表層物質(zhì)結(jié)構(gòu)進行了重塑,使得棄渣體邊坡表面大孔隙數(shù)量大大降低,從而促進了產(chǎn)流。在2種降雨強度條件下坡面徑流率隨降雨歷時均呈現(xiàn)先增加后穩(wěn)定的變化趨勢。

    雖然煤矸石在1.0 mm min-1表面未觀測到產(chǎn)流現(xiàn)象,說明徑流全部下滲形成壤中流,有可能進一步入滲至堆積體和原始地面界面容易形成滑動面,長時間的降雨也使堆積體內(nèi)部摩擦力減小,從而容易誘發(fā)大面積的滑坡和坍塌事故,尤其在長歷時強降雨條件下更易發(fā)生此類重力侵蝕事件;其次,該區(qū)露天煤礦排矸量較大,不僅壓占大量土地而且影響該區(qū)的景觀格局和生態(tài)環(huán)境質(zhì)量[23],因此,煤矸石堆積體不僅需要重點防護滑坡等重力侵蝕事件,還需對煤矸石堆積體進行必要的生態(tài)恢復(fù)。

    圖2 偏土質(zhì)棄渣體徑流率隨時間的變化Fig. 2 Temporal variation of soil runoff rate on the earthy slope

    圖3 偏石質(zhì)棄渣體徑流率隨時間的變化Fig. 3 Temporal variation of soil runoff rate on the rocky slope

    圖4 煤矸石棄渣體坡面徑流率隨時間的變化Fig. 4 Temporal variation of soil runoff rate on the slope of gangue

    2.1.2 不同防護措施的減水效益 圖5為1.0、1.5 mm min-1雨強下不同類型棄渣體在不同防護措施下的平均徑流率變化特征。偏土質(zhì)棄渣體在未防護、種草、魚鱗坑措施下的平均徑流率分別為3.39、1.94、1.63 L min-1和4.64、2.27、1.27 L min-1,2種措施的減水效益分別為42.91%、51.89%和51.21%、72.72%,可見,2種措施的減水效益十分明顯,且魚鱗坑的減水效益較種草高8.98%~21.51%,這是因為魚鱗坑具有通過改變坡面局部地形,增加地面粗糙度以削弱徑流的能力[11];1.5 mm min-1雨強時,2種措施的減水效益較1.0 mm min-1雨強增加8.3%和20.83%,這表明雨強越大,魚鱗坑的蓄水潛力較草被更大。偏石質(zhì)棄渣體在未防護、種草、魚鱗坑措施下的平均徑流率分別為1.71、1.26、1.33 L min-1和3.97、1.78、2.26 L min-1,種草和魚鱗坑的減水效率分別為26.28%、22.37%和55.20%、42.92%,種草對于偏石質(zhì)棄渣體坡面徑流的消減效果較魚鱗坑好,其減水效益較魚鱗坑措施高3.91%~12.28%。對于煤矸石棄渣體而言,1.0 mm min-1雨強下種草坡面徑流率較魚鱗坑措施降低17.88%;雨強1.5 mm min-1時,未防護和2種措施防護的坡面平均徑流率為2.31、1.82、2.07 L min-1,2種措施減水效益分別為21.32%和10.33%,草被的減水效益也優(yōu)于魚鱗坑,其減水效益較魚鱗坑措施高10.99%。與偏土質(zhì)相反,偏石質(zhì)、石質(zhì)(煤矸石)棄渣體種草措施減水效益均較魚鱗坑措施好,一方面可能是由于植被上部葉片攔截降雨,枝干又?jǐn)r截了一定的徑流,使流速降低,徑流能量減少,從而降低其挾沙能力。另一方面可能是由于土粒在沿坡面根莖連接處沉積構(gòu)成微型慮水土體,對徑流產(chǎn)生攔截和過濾作用[24];而礫石也具有改變徑流路徑削弱徑流的能力[25-26],二者的共同作用使得下滲作用更強;而魚鱗坑措施僅僅是改變了小地形,短時間內(nèi)攔截了降雨,效果略低于種草。

    圖5 不同防護措施棄渣體邊坡的平均徑流率Fig. 5 Variation of average runoff ratio on the slopes of slag

    2.2 棄渣體產(chǎn)沙特征

    2.2.1 侵蝕速率隨時間的變化 圖6~圖8分別為1.0、1.5 mm min-1雨強下偏土質(zhì)、偏石質(zhì)和煤矸石棄渣體在未防護、種草、布設(shè)魚鱗坑措施下邊坡的侵蝕速率隨產(chǎn)流歷時的變化。

    由圖可知,3種棄渣體在不同雨強條件下邊坡侵蝕速率均呈波動減小趨勢,且未防護減小趨勢較防護的明顯。對于偏土質(zhì)堆積體(圖6),未防護坡面侵蝕速率在0~6 min波動劇烈,而兩種防護措施下的侵蝕速率則波動較小,1.0 mm min-1和1.5 mm min-1雨強下侵蝕速率峰值分別達到9521 g min-1和15531 g min-1,遠大于2種防護措施的侵蝕速率,這是因為棄渣體結(jié)構(gòu)性、穩(wěn)定性差,不具有良好土壤的各項物理化學(xué)性質(zhì)[27],在產(chǎn)流初期表層粘結(jié)性較好的細顆粒容易被侵蝕,而留下穩(wěn)定性較差的大顆粒物質(zhì),因此導(dǎo)致降雨前期侵蝕速率波動較大。對于偏石質(zhì)堆積體(圖7),在1.0 mm min-1和1.5 mm min-1雨強條件下未防護坡面侵蝕速率隨時間的變化差異較大,種草和魚鱗坑措施防護坡面侵蝕速率差異較小,總體上呈現(xiàn)下降的趨勢;在1.5 mm min-1雨強時,偏石質(zhì)棄渣體同樣是在產(chǎn)流前期起伏較大,后期較小,原因可能是因為隨著坡面細顆粒逐漸被搬運,棄渣體中未被侵蝕的大顆粒之間相互支撐形成了穩(wěn)定渣床面[10],侵蝕速率逐漸穩(wěn)定下來,但相比偏土質(zhì)棄渣體則波動大,說明礫石的存在影響了降雨的侵蝕作用[24],使堆積體侵蝕差異性增大;對于煤矸石棄渣體(圖8),其侵蝕速率隨時間的變化呈多峰多谷型,這與煤矸石棄渣體的物質(zhì)結(jié)構(gòu)的復(fù)雜性關(guān)系密切。

    圖6 偏土質(zhì)棄渣體侵蝕速率隨時間的變化Fig. 6 Temporal variation of soil erosion rate on the earthy slope

    圖7 偏石質(zhì)棄渣體侵蝕速率隨時間的變化Fig. 7 Temporal variation of soil erosion rate on the rocky slope

    圖8 煤矸石棄渣體坡面侵蝕速率隨時間的變化Fig. 8 Temporal variation of soil erosion rate on the slope of gangue

    2.2.2 不同防護措施的減沙效益 圖9為1.0和1.5 mm min-1雨強下不同類型棄渣體在不同防護措施下平均侵蝕速率變化。由圖可知,偏土質(zhì)棄渣體邊坡平均侵蝕速率均遠大于偏石質(zhì)和煤矸石棄渣體。對于偏土質(zhì),2種防護措施的減沙效益十分明顯,1.0和1.5 mm min-1雨強時未防護坡面、種草和魚鱗坑坡面平均侵蝕速率分別為2411、49.31、37.75 g min-1和4355、107、27.43 g min-1,2種措施在2種雨強條件下減沙效益分別高達97.95%、98.41%和97.54%、99.30%,這表明魚鱗坑措施減沙效益較種草高出0.46%~1.76%,具有更高的防護效益[12]。雨強由1.0 mm min-1增大至1.5 mm min-1時,未防護坡面侵蝕速率增加81.64%,因此,對于礦區(qū)裸露偏土質(zhì)棄渣體的防護顯得尤為重要。與偏土質(zhì)棄渣體相比,偏石質(zhì)棄渣堆積體邊坡平均侵蝕速率要小很多,其未防護坡面在2種雨強下的平均侵蝕速率較偏土質(zhì)堆積體降低99.6%、96.2%,說明礫石的存在明顯減小了土壤可蝕性;2種雨強條件下2種措施的減沙效益分別為42.26%、41.98%和94.90%、91.84%,種草措施的減沙效益較魚鱗坑措施高0.28%~3.06%。煤矸石棄渣體由于所含成份復(fù)雜、含細小煤灰較多,黏性較大,故在1.0 mm min-1雨強下2種措施的侵蝕速率較小,分別為0.88、3.79 g min-1,雨強為1.5 mm min-1時,未防護坡面侵蝕速率為5.12 g min-1,相比其他2類棄渣體其侵蝕速率降低了很多,2種措施的減沙效益分別為39.50%和7.80%,種草措施的減沙效益較魚鱗坑措施高31.70%。

    2.3 棄渣體坡面水沙關(guān)系對不同防護措施的響應(yīng)

    圖9 不同防護措施棄渣體邊坡的平均侵蝕速率Fig. 9 Variation of average erosion rate on the slopes of slag

    表2 不同措施條件下3類棄渣體侵蝕速率與徑流率相關(guān)性Table 2 Correlation between erosion rate and runoff rate on slopes of slag relative to type of the slope and protection measure taken

    表2為3種棄渣體在未防護、種草、魚鱗坑措施下,侵蝕速率與徑流率的擬合方程。由表2可知,3種棄渣體坡面在未采取防護措施的情況下,侵蝕速率與徑流率均呈極顯著線性關(guān)系。其中,偏土質(zhì)和偏石質(zhì)棄渣體侵蝕速率隨徑流率的增大而增大,而煤矸石棄渣體坡面上侵蝕速率隨徑流率增大而減小,原因是因為粉煤灰和煤矸石的混合物具有與土壤相似的水力學(xué)特性[22],在降雨初期,水分下滲較多,導(dǎo)致坡面徑流量小,而坡面易蝕細顆粒大量存在,故侵蝕量較大;隨著降雨歷時的延長,土層含水量逐漸達到飽和,水分入滲減弱,坡面細顆粒減少,致使徑流率大而侵蝕量小。

    在種草措施下,偏土質(zhì)棄渣體侵蝕速率和徑流率呈顯著對數(shù)關(guān)系,而在偏石質(zhì)和煤矸石棄渣體坡面上二者相關(guān)性不顯著;魚鱗坑措施下,除偏石質(zhì)棄渣體的侵蝕速率和徑流率具有極顯著線性關(guān)系外,在偏土質(zhì)和煤矸石棄渣體坡面上,二者相關(guān)性未達到顯著性水平。這表明,這2種防護措施對棄土棄渣體邊坡土壤侵蝕速率與徑流率的關(guān)系產(chǎn)生了重要影響[12,18]。對于植被措施而言,(1)冰草是一種良好的水土保持和防風(fēng)固沙植物,其根為須狀,密生且具有入土較深的特點[28],因而可深入土壤深層,提高土壤透水性和保水能力,降低地表徑流和沖刷。(2)植被根系能夠改變土壤疏松程度,增加土壤孔隙度,根系對土體具有拉伸和固持作用,有利于增加降雨入滲能力[29]和土壤持水能力,同時使土壤表層出現(xiàn)結(jié)皮[30],植被促進入滲能力主要是因為地表覆蓋層和根系土壤層對水分的截留和吸收,進一步延緩徑流流速,從而減弱了降雨徑流對表層土壤的沖刷作用;(3)植被冠層對降雨的攔截作用較強,可降低雨滴的動能[31],減輕雨滴對地面的直接打擊作用,降低擊濺侵蝕,致使侵蝕速率隨之降低,最終導(dǎo)致侵蝕速率和徑流率的相關(guān)性較小。就魚鱗坑措施而言,其作為水保工程措施能夠改變局部小地形,在降雨期間能夠有效蓄積儲存降雨,攔蓄坡面降雨徑流,增強降雨的疊加入滲能力,降低土壤侵蝕[32-33]。因此,種草和魚鱗坑措施處理坡面可以打破侵蝕速率與徑流率的一般關(guān)系,改變坡面水沙關(guān)系。

    3 結(jié) 論

    采用野外模擬降雨的方法,以神府礦區(qū)不同類型棄土棄渣體為研究對象,研究1.0和1.5 mm min-1雨強條件下種草和魚鱗坑措施對棄土棄渣堆積體坡面徑流率、侵蝕速率及水沙關(guān)系的影響。不同措施下3種棄渣體邊坡徑流率均在產(chǎn)流6~9 min后趨于穩(wěn)定,產(chǎn)流過程中棄渣體邊坡侵蝕速率均呈波動減小趨勢,且未防護坡面減小趨勢較防護更明顯;種草對3種棄渣體的減水和減沙效益分別為10.33%~55.20%和7.80%~97.95%;魚鱗坑分別為21.32%~72.72%和39.50%~99.30%。魚鱗坑措施對偏土質(zhì)棄渣體的防護效益優(yōu)于種草措施,而種草措施對偏石質(zhì)棄渣體和煤矸石的防護效益優(yōu)于魚鱗坑措施。未防護條件下3種棄渣體坡面的侵蝕速率與徑流率呈現(xiàn)顯著的線性關(guān)系(p<0.05),而種草和魚鱗坑措施改變了坡面水沙關(guān)系;盡管煤矸石堆積體在1.0 mm min-1雨強下坡面無徑流產(chǎn)生,但考慮到生態(tài)安全,進行相關(guān)的生態(tài)修復(fù)是十分必要的。本試驗未進行更多坡度、更大雨強及更長歷時的試驗,且措施種類也比較單一,研究中也未考慮到壤中流的變化對土壤侵蝕的影響,今后應(yīng)加強這幾個方面的研究。

    [1] 楊選民,丁長印. 神府東勝礦區(qū)生態(tài)環(huán)境問題及對策.能源環(huán)境保護,2000,14(1):69-72 Yang X M,Ding C Y. Ecological environmental problems and countermeasures in ShenfuDongsheng mining area(In Chinese). Energy Environmental Protection,2000,14(1):69—72

    [2] 徐占軍,侯湖平,張紹良,等. 采礦活動和氣候變化對煤礦區(qū)生態(tài)環(huán)境損失的影響. 農(nóng)業(yè)工程學(xué)報,2012,28(5):232—240 Xu Z J,Hou H P,Zhang S L,et al. Effects of mining activity and climatic change on ecological losses in coal mining areas(In Chinese). Transactions of the Chinese Society of Agricultural Engineering,2012,28(5):232—240

    [3] Ahirwal J,Maiti S K,Singh A K.Changes in ecosystem carbon pool and soil CO2flux following post-mine reclamation in dry tropical environment,India. Science of the Total Environment,2017,583(1):153—162

    [4] Chalupnik S,Wysocka M,Janson E,et al. Long term changes in the concentration of radium in discharge waters of coal mines and Upper Silesian rivers. Journal of Environmental Radioactivity,2017,171(2):117—123

    [5] 趙暄,謝永生,王允怡,等. 模擬降雨條件下棄土堆置體侵蝕產(chǎn)沙試驗研究. 水土保持學(xué)報,2013,27(3):1—8,76 Zhao X,Xie Y S,Wang Y Y,et al. Experimental study on soil erosion and sediment yield of spoil bank under simulated rainfall condition(In Chinese).Journal of Soil and Water Conservation,2013,27(3):1—8,76

    [6] 馬春艷,王占禮,寇曉梅,等. 工程建設(shè)棄土棄渣水土流失過程試驗研究. 水土保持通報,2009,23(3):78—82 Ma C Y,Wang Z L,Kou X M,et al. Experimental study of soil and water loss processes on waste soil and residue in project construction(In Chinese). Bulletin of Soil and Water Conservation,2009,23(3):78—82

    [7] Riley S J.Aspects of the differences in the erodibility of the waste rock dump and natural surfaces,Ranger Uranium Mine,Northern Territory,Australia. Applied Geography,1995,15(4):309—323

    [8] Zhou W,Yang K,Bai Z,et al. The development of topsoil properties under different reclaimed land uses in the Pingshuo opencast coalmine of Loess Plateau of China. Ecological Engineering,2017,100:237—245

    [9] 康宏亮,王文龍,薛智德,等. 陜北風(fēng)沙區(qū)含礫石工程堆積體坡面產(chǎn)流產(chǎn)沙試驗. 水科學(xué)進展,2016,27(2):256—265 Kang H L,Wang W L,Xue Z D,et al. Experimental study on runoff and sediment yield from engineering deposition with gravel in the northern windy-sandy region,Shaanxi(In Chinese). Advances in Water Science,2016,27(2):256-265

    [10] 郭明明,王文龍,李建明,等. 神府礦區(qū)棄土棄渣體侵蝕特征及預(yù)測. 土壤學(xué)報,2015,52(5):1044—1057 Guo M M,Wang W L,Li J M,et al.Erosion on dunes of overburden and waste slag in Shenfu coal field and prediction(In Chinese). Acta Pedologica Sinica,2015,52(5):1044—1057

    [11] 王貞. 神東煤田不同下墊面侵蝕產(chǎn)沙規(guī)律及水動力參數(shù)特征. 陜西:西北農(nóng)林科技大學(xué),2011 Wang Z. The feature of sediment yield and hydrodynamics of different underlying surfaces in the construction of the Shenfu-Dongsheng coalfield(In Chinese).Shaanxi:Northwest Aamp;F Universtiy,2011

    [12] 田野. 堆渣邊坡兩種護坡措施的水土保持效益研究. 北京:北京林業(yè)大學(xué),2015 Tian Y. Study on benefits of soil and water conservation slag heap slope revetment two measures(In Chinese).Beijing:Beijing Forestry University,2015

    [13] 王貞,王文龍,金劍,等. 神東煤田擾動地面與原地面產(chǎn)流產(chǎn)沙及水動力學(xué)參數(shù)對比. 中國水土保持科學(xué),2010,8(6):69—74 Wang Z,Wang W L,Jin J,et al. Comparison of runoff generation,sediment yield and hydrodynamic parameters on the undisturbed and the disturbed land surface in the Shenfu-Dongsheng coalfield(In Chinese). Science of Soil and Water Conservation,2010,8(6):69—74

    [14] 郭明明,王文龍,李建明,等. 神府煤田土壤顆粒分形及降雨對徑流產(chǎn)沙的影響. 土壤學(xué)報,2014,51(5):983—992 Guo M M,Wang W L,Li J M,et al. The influence of soil particles fractal and rainfall on runoff and sediment yield in Shenfu coalfield(In Chinese). Acta Pedologica Sinica,2014,51(5):983—992

    [15] 丁文斌,史東梅,何文健,等. 放水沖刷條件下工程堆積體邊坡徑流侵蝕水動力學(xué)特性. 農(nóng)業(yè)工程學(xué)報,2016,32(18):153—161 Ding W B,Shi D M,He W J,et al. Hydrodynamic characteristics of engineering accumulation erosion under side slope runoff erosion process in field scouring experiment(In Chinese). Transactions of the Chinese Society of Agricultural Engineering,2016,32(18):153—161

    [16] 丁文斌,李葉鑫,史東梅,等. 兩種工程堆積體邊坡模擬徑流侵蝕對比研究. 土壤學(xué)報,2017,54(3):558—569 Ding W B,Li Y X,Shi D M,et al. Contrast study on simulated runoff erosion of two engineering accumulation slopes(In Chinese). Acta Pedologica Sinica,2017,54(3):558—569

    [17] 張樂濤,高照良,田紅衛(wèi). 工程堆積體陡坡坡面土壤侵蝕水動力學(xué)過程. 農(nóng)業(yè)工程學(xué)報,2013,29(24):94—102 Zhang L T,Gao Z L,Tian H W. Hydrodynamic process of soil erosion in steep slope of engineering accumulation(In Chinese). Transactions of the Chinese Society of Agricultural Engineering,2013,29(24):94—102

    [18] 劉瑞順,王文龍,廖超英,等. 露天煤礦排土場邊坡防護措施減水減沙效益分析. 西北林學(xué)院學(xué)報,2014(4):59—64 Liu R S,Wang W L,Liao C Y,et al. Benefits in runfall and sediment reductions of the protection methods for the dump slope of opencast coal mine(In Chinese). Journal of Northwest Forestry University,2014(4):59—64

    [19] 倪含斌,張麗萍,張登榮. 模擬降雨試驗研究神東礦區(qū)不同階段堆積棄土的水土流失. 環(huán)境科學(xué)學(xué)報,2006,26(12):2065—2071 Ni H B,Zhang L P,Zhang D R. Simulated rainfall experiments on soil and water erosion in different phases of resources exploitation in Shendong mine region(In Chinese). Acta Scientiae Circumstantiae,2006,26(12):2065—2071

    [20] 雷廷武,張晴雯,趙軍,等. 確定侵蝕細溝集中水流剝離速率的解析方法. 土壤學(xué)報,2002,39(6):788—793 Lei T W,Zhang Q W,Zhao J,et al. Analytic method for determination of detachment rate of concentrated flow in erosion rills(InChinese). Acta Pedologica Sinica,2002,39(6):788—793

    [21] 符素華. 土壤中礫石存在對入滲影響研究進展. 水土保持學(xué)報,2005,19(1):171—175 Fu S H. Effect of soil containing rock fragment on infiltration(In Chinese). Journal of Soil and Water Conservation,2005,19(1):171—175

    [22] Wang J,Li X,Bai Z,et al. The effects of coal gangue and fly ash on the hydraulic properties and water content distribution in reconstructed soil profiles of coalmined land with a high groundwater table. Hydrological Processes,2017,31:687—697

    [23] 勾立新. 煤矸石的環(huán)境問題及對策初探. 鐵法科技,2015,B12:208—211 Gou L X. Environment problem and discussion on the countermeasures of coal gangue(InChinese). Science and Technology of Iron Processes. 2015,B12:208—211

    [24] 楊亞川,莫永京,王芝芳,等. 土壤-草本植被根系復(fù)合體抗水蝕強度與抗剪強度的試驗研究. 中國農(nóng)業(yè)大學(xué)學(xué)報,1996,1(2):31—38 Yang Y C,Mo Y J,Wang Z F,et al. Experimental study on anti-water erosion and shear strength of soil-root composite(In Chinese). Journal of China Agricultural University,1996,1(2):31—38

    [25] 史東梅,蔣光毅,彭旭東,等. 不同土石比的工程堆積體邊坡徑流侵蝕過程. 農(nóng)業(yè)工程學(xué)報,2015,31(17):152—161 Shi D M,Jiang G Y,Peng X D,et al. Runoff erosion process on slope of engineering accumulation-with different soil-rock ratio(In Chinese). Transactions of the Chinese Society of Agricultural Engineering,2015,31(17):152—161

    [26] Wilcox B P,Wood M K. Factors influencing interrill erosion from semiarid slopes in New Mexico Journal of Range Management,1989,42(1):66—70

    [27] 李建明,王文龍,王貞,等. 神府東勝煤田棄土棄渣體徑流產(chǎn)沙過程的野外試驗. 應(yīng)用生態(tài)學(xué)報,2013,24(12):3537—3545 Li J M,Wang W L,Wang Z,et al. A field experiment of runoff and sediment yielding processes from residues in Shenfu-Dongsheng coalfield(In Chinese). Chinese Journal of Applied Ecology,2013,24(12):3537—3545

    [28] 曾亮,袁慶華,王方,等. 冰草屬植物種質(zhì)資源遺傳多樣性的ISSR分析. 草業(yè)學(xué)報,2013,22(1):260—267 Zeng L,Yuan Q H,Wang F,et al. Gentic diversity analysis of Agropyron germplasm resource by ISSR(In Chinese). Acta Prataculturae Sinica,2013,22(1):260—267

    [29] Rpc M,Mcintyre K,Vickers AW,et al. A rainfall simulation study of soil erosion rangeland in Swaziland.Soil Technology,1997,11(3):291—299

    [30] 吳發(fā)啟,范文波. 土壤結(jié)皮對降雨入滲和產(chǎn)流產(chǎn)沙的影響. 中國水土保持科學(xué),2005,3(2):97—101 Wu F Q,F(xiàn)an W B. Effects of soil encrustation on rainfall infiltration,runoff and sediment generation(In Chinese). Science of Soil and Water Conservation,2005,3(2):97—101

    [31] 楊勝利,王云鵬. 排土場穩(wěn)定性影響因素分析. 露天采礦技術(shù),2009(3):4—7 Yang S L,Wang Y P. Analysis of factors affecting dump stability(In Chinese). Opencast Mining Technology,2009(3):4—7

    [32] 李虹辰,趙西寧,高曉東,等. 魚鱗坑與覆蓋組合措施對陜北旱作棗園土壤水分的影響. 應(yīng)用生態(tài)學(xué)報,2014,25(8):2297—2303 Li H C,Zhao X N,Gao X D,et al. Effects of scalelike pit and mulching measures on soil moisture of dryland jujube orchard in North Shaanxi Province,China(In Chinese). Chinese Journal of Applied Ecology,2014,25(8):2297—2303

    [33] 張北贏,徐學(xué)選,劉文兆. 黃土丘陵溝壑區(qū)不同水保措施條件下土壤水分狀況. 農(nóng)業(yè)工程學(xué)報,2009,25(4):54—58 Zhang B Y,Xu X X,Liu W Z. Soil water condition under different measures of soil and waterconservation in loess hilly and gully region(In Chinese). Transactions of the Chinese Society of Agricultural Engineering,2009,25(4):54—58

    (責(zé)任編輯:檀滿枝)

    Study on Runoff and Sediment Reducing Effects of Erosion Control Measures on Slopes of Slag Dumps with Simulated Rainfall

    YANG Bo1WANG Wenlong1,2,3?GUO Mingming1KANG Hongliang1BAI Yun3LI Jianming4
    (1 State Key Laboratory of Soil Erosion and Dryland Farming on Loess Plateau,Institute of Soil and Water Conservation,Northwest Aamp;F Universtiy,Yangling,Shaanxi 712100,China)
    (2 Institute of Soil and Water Conservation,Chinese Academy of Sciences and Ministry of Water Resources,Yangling,Shaanxi 712100,China.)
    (3 Key Laboratory of Ecological Restoration of Mining Area in Northern Shaanxi Province;Yulin University;Yulin,Shaanxi 719000,China)
    (4 Department of Soil and Water Conservation,Yangtze River Scientific Research Institute,Wuhan 430010,China)

    【Objective】The ecological and environmental problems caused by mining,especially coal mining,have aroused worldwide concerns and attention. It is,therefore,essential to explore rapid and effective ways to rehabilitate the ecological systems in the mining areas.【Method】Outdoor simulation of rainfall events was conducted with rainfall intensity controlled at 1.0 mm min-1or 1.5 mm min-1on slopes,different in texture,1)earthy slope(more soil than rock);2)rocky slope(more rock than soil):and 3)slope of gangue,to evaluate effects of the erosion control measures,planting grass and digging fishscale pit. The two years old grass of Agropyron was formerly planted in the farmlands around the mining and transplanted on the slopes and the fish-scale pits dug temporarily. The rainfall simulator was a troughtype artificial simulator. When the nozzle is 3 m above the ground,the raindrops may reach the end point speed. Plots of the artificial slopes were 3 m×1 each in area and 35° in slope gradient. Before the start of the experiment,rainfall was calibrated in intensity until it reached 85% or higher in uniformity coefficient. During the rainfall events,runoff flow velocity was measured with the dye tracing method. During the first 3 minutes of rainfall,samples of runoff and sediment were gathered every minute,and after that once every 3 minutes.Sediment samples were oven-dried for weighing.【Result】(1)Runoff leveled off in rate in 6minutes after the initiation of runoff on all the three types of slopes,regardless of type of control measure. The erosion rate on the slope of gangue tended to fluctuate,while declining and the trend was more obvious than those on the protected slopes.(2)Planting grass reduced runoff by 42.91%~51.21%,26.28%~55.20% and 10.33% and sediment by 97.54%~97.95%,41.87%~42.26% and7.80%,respectively,on earthy slopes,rocky slopes and slopes of gangue,while digging fish-scale pits did by 51.89%~72.72%,22.37%~42.92% and 21.32%and by 98.41%~99.30%,94.90%~91.84% and 39.50%,respectively.(3)On earthy slopes,digging fishscale pits was 8.98%~21.51% and 0.46%~1.76%,respectively,higher in runoff and sediment reduction rate than planting grass,while on rocky slopesand slopes of gangue,planting grass was 3.91%~12.28% and 0.28%~3.06% higher and 10.99% and 31.70% higher,respectively.(4)On all the three types of slopes without protection measures taken,linear relationships(R2= 0.775,p<0.01;R2= 0.649,p<0.01;R2=0.450,p<0.05)were found between erosion rate and runoff rate,but they were weakened or fading out altered when the protection measures were taken.【Conclusion】All findings in this study may serve as scientific basis in guiding ecological restoration of slopes of gangue in the mining area.

    Shenfu coalfield;Mining slags;Fish-scale pit;Grass planting;Runoff rate;Erosion rate

    S157.1

    A

    10.11766/trxb201705270151

    * 國家重點研發(fā)計劃重點專項項目(YS2016YFSF030019)和國家自然科學(xué)基金項目(40771127)資助 Supported by the National Key Research and Development Program of China(No.YS2016YFSF030019)and the National Natural Science Foundation of China(No.40771127)

    ? 通訊作者 Corresponding author:wlwang@nwsuaf.edu.cn

    楊 波(1993—),男,陜西商南人,碩士研究生。主要從事開發(fā)建設(shè)項目水土保持。E-mail:584431911@qq.com

    2017-05-27;

    2017-07-10;優(yōu)先數(shù)字出版日期(www.cnki.net):2017-07-25

    猜你喜歡
    渣體魚鱗坑堆積體
    云南某冶煉廢渣堆滲濾液產(chǎn)生及其影響因素研究
    太行山區(qū)困難立地“五步造林法”
    不同規(guī)格魚鱗坑坡面侵蝕過程及特征研究
    后退式注漿技術(shù)在隧道突涌堆積體加固處治中的應(yīng)用
    棄渣體潛在失穩(wěn)滑動面探討
    魚鱗坑生態(tài)恢復(fù)措施對寧夏典型草原植物群落特征的影響
    隧道黃土堆積體施工技術(shù)
    金沙江某水電站100m級高棄渣體穩(wěn)定性研究與評價
    大型堆積體滑坡治理施工技術(shù)研究
    高密度電法在尋找泥石流堆積體中的應(yīng)用研究
    在线免费十八禁| 成人一区二区视频在线观看| 免费观看在线日韩| 日日摸夜夜添夜夜添小说| av黄色大香蕉| 如何舔出高潮| 一卡2卡三卡四卡精品乱码亚洲| 别揉我奶头~嗯~啊~动态视频| 久久久久久国产a免费观看| 久久婷婷人人爽人人干人人爱| 国产视频一区二区在线看| 一区二区三区激情视频| 又黄又爽又免费观看的视频| 色综合色国产| 国产精品国产高清国产av| 亚洲经典国产精华液单| 久久久久国产精品人妻aⅴ院| 网址你懂的国产日韩在线| av天堂中文字幕网| 特级一级黄色大片| a级一级毛片免费在线观看| 国产单亲对白刺激| 国产精品一区www在线观看 | 国产色爽女视频免费观看| 亚洲欧美精品综合久久99| 国产单亲对白刺激| 91麻豆av在线| 一本精品99久久精品77| 精品99又大又爽又粗少妇毛片 | 男女啪啪激烈高潮av片| 不卡视频在线观看欧美| 免费av不卡在线播放| 亚洲,欧美,日韩| 久久久精品大字幕| 国产精品美女特级片免费视频播放器| 亚洲aⅴ乱码一区二区在线播放| 美女黄网站色视频| 久久精品久久久久久噜噜老黄 | 亚洲精品乱码久久久v下载方式| 很黄的视频免费| 亚洲av不卡在线观看| 中文字幕精品亚洲无线码一区| 午夜福利在线观看吧| 国产精品三级大全| 国产精品人妻久久久影院| 老熟妇仑乱视频hdxx| 色综合婷婷激情| 日日撸夜夜添| 人人妻人人澡欧美一区二区| 2021天堂中文幕一二区在线观| 国产极品精品免费视频能看的| 久久久久久久久中文| 国产精品一区二区三区四区免费观看 | 久久人妻av系列| 天堂影院成人在线观看| 国产在视频线在精品| 久久久精品大字幕| 国产乱人伦免费视频| 日本免费一区二区三区高清不卡| 日韩欧美国产一区二区入口| 国产高清激情床上av| 国内毛片毛片毛片毛片毛片| 制服丝袜大香蕉在线| 欧美黑人欧美精品刺激| 欧美bdsm另类| 网址你懂的国产日韩在线| 波多野结衣巨乳人妻| 国产精品永久免费网站| 在线免费观看不下载黄p国产 | 最近在线观看免费完整版| 国产精品一区二区三区四区久久| 国产高清视频在线播放一区| 91狼人影院| 淫妇啪啪啪对白视频| 亚洲人成伊人成综合网2020| 18禁在线播放成人免费| 色综合亚洲欧美另类图片| 久久久久久九九精品二区国产| 99热这里只有是精品50| 男人狂女人下面高潮的视频| 国产成人福利小说| 亚洲电影在线观看av| 午夜视频国产福利| 久久精品国产鲁丝片午夜精品 | 一本精品99久久精品77| 精品99又大又爽又粗少妇毛片 | 亚洲无线观看免费| 国产成人影院久久av| 男女之事视频高清在线观看| 日韩强制内射视频| 国产蜜桃级精品一区二区三区| 成人国产一区最新在线观看| 一进一出抽搐gif免费好疼| 亚洲欧美日韩东京热| 看黄色毛片网站| 啦啦啦韩国在线观看视频| 欧美另类亚洲清纯唯美| 69人妻影院| 天堂√8在线中文| 香蕉av资源在线| 亚洲人成网站在线播放欧美日韩| 亚洲avbb在线观看| www.www免费av| 亚洲av美国av| 不卡视频在线观看欧美| 国产 一区 欧美 日韩| 成年人黄色毛片网站| 有码 亚洲区| 日本色播在线视频| 精品一区二区三区av网在线观看| 久久久国产成人免费| 久久6这里有精品| 欧美黑人欧美精品刺激| 免费av不卡在线播放| 在线免费观看不下载黄p国产 | 久久亚洲真实| 久久人人精品亚洲av| 国产色爽女视频免费观看| 一夜夜www| 久久精品国产亚洲av香蕉五月| 人妻夜夜爽99麻豆av| 在线播放无遮挡| 女人被狂操c到高潮| 人妻少妇偷人精品九色| 一边摸一边抽搐一进一小说| 午夜a级毛片| 欧美精品啪啪一区二区三区| 国产精品不卡视频一区二区| 亚洲久久久久久中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 99热这里只有是精品50| 在线观看免费视频日本深夜| 中文资源天堂在线| 日本撒尿小便嘘嘘汇集6| 别揉我奶头~嗯~啊~动态视频| 女人被狂操c到高潮| 99热这里只有是精品50| 联通29元200g的流量卡| 亚洲欧美日韩高清在线视频| 国产精品一区二区三区四区免费观看 | 非洲黑人性xxxx精品又粗又长| 亚洲欧美精品综合久久99| 精品久久久久久久久久免费视频| 午夜精品一区二区三区免费看| 国产精品不卡视频一区二区| 熟女电影av网| 亚洲专区中文字幕在线| 婷婷六月久久综合丁香| 一夜夜www| 久久九九热精品免费| 欧美色欧美亚洲另类二区| 别揉我奶头~嗯~啊~动态视频| 久久久国产成人精品二区| 国产精品久久电影中文字幕| 国产亚洲欧美98| 热99re8久久精品国产| 欧美xxxx性猛交bbbb| 三级毛片av免费| 日韩欧美在线二视频| 美女被艹到高潮喷水动态| 国产精品98久久久久久宅男小说| 99久久久亚洲精品蜜臀av| 又紧又爽又黄一区二区| 午夜福利在线观看免费完整高清在 | 国产精品亚洲一级av第二区| 麻豆成人av在线观看| 国产午夜福利久久久久久| 黄色配什么色好看| 亚洲va日本ⅴa欧美va伊人久久| 国产在线精品亚洲第一网站| 波多野结衣巨乳人妻| 九色国产91popny在线| 精品一区二区三区av网在线观看| 特大巨黑吊av在线直播| 亚洲av电影不卡..在线观看| 麻豆精品久久久久久蜜桃| 特大巨黑吊av在线直播| 欧美成人性av电影在线观看| 搞女人的毛片| 99热这里只有是精品50| 中亚洲国语对白在线视频| 久久人人精品亚洲av| 嫩草影院入口| avwww免费| 亚洲美女黄片视频| 在线观看免费视频日本深夜| 亚洲精品一卡2卡三卡4卡5卡| 日韩一区二区视频免费看| 一个人免费在线观看电影| 久久久久久大精品| 日韩欧美一区二区三区在线观看| 久久中文看片网| 97碰自拍视频| 亚洲成人免费电影在线观看| 2021天堂中文幕一二区在线观| 久久99热这里只有精品18| or卡值多少钱| 日韩欧美免费精品| 久久久久久久久久黄片| 亚洲性夜色夜夜综合| 男人和女人高潮做爰伦理| 91麻豆av在线| 免费一级毛片在线播放高清视频| 尾随美女入室| 国内精品美女久久久久久| 搡女人真爽免费视频火全软件 | 日本黄色视频三级网站网址| 精品一区二区免费观看| 免费不卡的大黄色大毛片视频在线观看 | 久久精品国产亚洲av香蕉五月| 九九热线精品视视频播放| 国产成年人精品一区二区| 99热网站在线观看| a级毛片a级免费在线| 99精品久久久久人妻精品| 内地一区二区视频在线| 国产免费男女视频| 精品人妻1区二区| 日日啪夜夜撸| 亚洲av美国av| 亚洲av成人av| 国产精品久久久久久av不卡| 日韩亚洲欧美综合| 日本欧美国产在线视频| 啦啦啦韩国在线观看视频| 国产精华一区二区三区| 夜夜爽天天搞| 好男人在线观看高清免费视频| 欧美3d第一页| 天堂av国产一区二区熟女人妻| 国产毛片a区久久久久| 成人三级黄色视频| 中文字幕熟女人妻在线| 午夜福利视频1000在线观看| 亚洲美女搞黄在线观看 | 中文字幕高清在线视频| 日韩中文字幕欧美一区二区| 久久久久久久久大av| 狠狠狠狠99中文字幕| 别揉我奶头 嗯啊视频| 老熟妇乱子伦视频在线观看| 国产免费av片在线观看野外av| 亚洲七黄色美女视频| 亚洲男人的天堂狠狠| 色综合婷婷激情| 少妇人妻一区二区三区视频| 久久6这里有精品| 91精品国产九色| 日韩,欧美,国产一区二区三区 | 少妇熟女aⅴ在线视频| 可以在线观看毛片的网站| 两性午夜刺激爽爽歪歪视频在线观看| 日韩中文字幕欧美一区二区| 亚洲国产高清在线一区二区三| 精华霜和精华液先用哪个| 亚洲av日韩精品久久久久久密| 乱人视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产视频内射| 国产精品不卡视频一区二区| 亚洲精品亚洲一区二区| 亚洲av免费在线观看| 亚洲欧美日韩卡通动漫| 亚洲精品久久国产高清桃花| 精品无人区乱码1区二区| 久久人妻av系列| 国产精品1区2区在线观看.| 少妇丰满av| 久久久精品欧美日韩精品| 日本欧美国产在线视频| 国产精品伦人一区二区| 国产精品电影一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 村上凉子中文字幕在线| 五月玫瑰六月丁香| 亚洲精品久久国产高清桃花| 黄色欧美视频在线观看| 色哟哟哟哟哟哟| 亚洲av成人av| av女优亚洲男人天堂| 真实男女啪啪啪动态图| 亚洲第一区二区三区不卡| 国产精品,欧美在线| 啪啪无遮挡十八禁网站| 精品午夜福利视频在线观看一区| 日韩在线高清观看一区二区三区 | 黄色欧美视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 91久久精品电影网| 亚洲av免费在线观看| 国产亚洲精品久久久com| 日日夜夜操网爽| 国内精品一区二区在线观看| 少妇高潮的动态图| 国产精品福利在线免费观看| 亚洲精品粉嫩美女一区| 人妻少妇偷人精品九色| 免费一级毛片在线播放高清视频| 12—13女人毛片做爰片一| 亚洲精品粉嫩美女一区| 国内精品久久久久久久电影| 国产午夜精品论理片| 日韩欧美三级三区| 日本撒尿小便嘘嘘汇集6| 国产精品98久久久久久宅男小说| a级毛片免费高清观看在线播放| 高清在线国产一区| 国产精品98久久久久久宅男小说| 国产成人福利小说| 51国产日韩欧美| 亚洲色图av天堂| 99在线人妻在线中文字幕| 99久国产av精品| 日本黄大片高清| 深爱激情五月婷婷| 在线看三级毛片| 白带黄色成豆腐渣| 狂野欧美激情性xxxx在线观看| 久久久国产成人精品二区| 白带黄色成豆腐渣| 亚洲av成人av| 久久久国产成人精品二区| 在线天堂最新版资源| 欧美性感艳星| 九九爱精品视频在线观看| 国产高潮美女av| 国产黄色小视频在线观看| 国产精品永久免费网站| 亚洲专区国产一区二区| 18禁黄网站禁片免费观看直播| 毛片女人毛片| 少妇高潮的动态图| 床上黄色一级片| 非洲黑人性xxxx精品又粗又长| 99久久久亚洲精品蜜臀av| 给我免费播放毛片高清在线观看| 全区人妻精品视频| 欧美黑人巨大hd| 色综合色国产| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美日韩东京热| 制服丝袜大香蕉在线| 色5月婷婷丁香| 我要看日韩黄色一级片| 国产成人一区二区在线| 久久久精品欧美日韩精品| 看片在线看免费视频| 国产成人福利小说| 波多野结衣巨乳人妻| 国产一区二区亚洲精品在线观看| 看黄色毛片网站| 欧美一区二区国产精品久久精品| 免费一级毛片在线播放高清视频| 国产一区二区在线av高清观看| 亚洲av成人精品一区久久| .国产精品久久| 91久久精品国产一区二区成人| 成人高潮视频无遮挡免费网站| 校园人妻丝袜中文字幕| 国产精品乱码一区二三区的特点| 国产黄a三级三级三级人| 十八禁国产超污无遮挡网站| 亚洲av成人av| 真人做人爱边吃奶动态| a级毛片免费高清观看在线播放| 国产激情偷乱视频一区二区| or卡值多少钱| 伦精品一区二区三区| 舔av片在线| 日韩欧美在线乱码| 我要搜黄色片| 国产精品国产高清国产av| 久久6这里有精品| 日日摸夜夜添夜夜添av毛片 | 亚洲三级黄色毛片| 久久久久久国产a免费观看| 99久久精品一区二区三区| 国产精品三级大全| netflix在线观看网站| 国产精品自产拍在线观看55亚洲| 天堂网av新在线| 久久精品影院6| 成人欧美大片| 乱人视频在线观看| 欧美精品啪啪一区二区三区| 性插视频无遮挡在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 欧美xxxx黑人xx丫x性爽| 亚洲精品一卡2卡三卡4卡5卡| 色精品久久人妻99蜜桃| 色综合色国产| 日韩欧美免费精品| 免费在线观看日本一区| 国产蜜桃级精品一区二区三区| 国产成人aa在线观看| 亚洲国产日韩欧美精品在线观看| 精品一区二区免费观看| 如何舔出高潮| ponron亚洲| 亚洲五月天丁香| 国产亚洲精品综合一区在线观看| www.色视频.com| 久久99热6这里只有精品| 在线a可以看的网站| 国产精品精品国产色婷婷| av黄色大香蕉| 少妇人妻一区二区三区视频| 国产真实乱freesex| 在线观看舔阴道视频| 88av欧美| 成人性生交大片免费视频hd| 国产精品乱码一区二三区的特点| 国产一区二区三区在线臀色熟女| 赤兔流量卡办理| 国产男人的电影天堂91| 天堂影院成人在线观看| 国产老妇女一区| 老熟妇乱子伦视频在线观看| 亚洲18禁久久av| 麻豆精品久久久久久蜜桃| 精品不卡国产一区二区三区| 我要看日韩黄色一级片| 性欧美人与动物交配| 日韩欧美国产一区二区入口| 在现免费观看毛片| 国产乱人伦免费视频| 亚洲熟妇中文字幕五十中出| 久久精品国产鲁丝片午夜精品 | 成人无遮挡网站| 黄色欧美视频在线观看| 亚洲av美国av| 99久久成人亚洲精品观看| 亚洲成人精品中文字幕电影| 久久99热6这里只有精品| 欧美日韩综合久久久久久 | 一本久久中文字幕| 91在线观看av| 色吧在线观看| 成人亚洲精品av一区二区| 22中文网久久字幕| 亚洲四区av| 日韩国内少妇激情av| 日韩人妻高清精品专区| 精品久久久久久久末码| 国产亚洲精品综合一区在线观看| 嫩草影院精品99| 天堂av国产一区二区熟女人妻| 亚州av有码| 免费一级毛片在线播放高清视频| 日本免费a在线| 成人国产综合亚洲| 高清毛片免费观看视频网站| 又紧又爽又黄一区二区| 男人和女人高潮做爰伦理| 97碰自拍视频| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产清高在天天线| 国产精品电影一区二区三区| 午夜亚洲福利在线播放| 特级一级黄色大片| 99riav亚洲国产免费| 日本黄大片高清| 久久久久性生活片| 搡女人真爽免费视频火全软件 | 黄色欧美视频在线观看| 国产精品爽爽va在线观看网站| 久久久久免费精品人妻一区二区| av在线天堂中文字幕| 成年版毛片免费区| 精品一区二区三区视频在线| 深夜a级毛片| av黄色大香蕉| 国产精品三级大全| 看免费成人av毛片| www.www免费av| 在线观看av片永久免费下载| 欧美丝袜亚洲另类 | 免费看日本二区| 午夜福利在线在线| 欧美在线一区亚洲| 国产淫片久久久久久久久| 亚洲av五月六月丁香网| 欧美性猛交黑人性爽| 可以在线观看毛片的网站| 亚州av有码| 女人被狂操c到高潮| 国产av在哪里看| 午夜福利视频1000在线观看| 日本成人三级电影网站| 久久久久久久精品吃奶| 中文字幕久久专区| 天堂动漫精品| 亚洲三级黄色毛片| 成人特级av手机在线观看| 国产色婷婷99| 久久久成人免费电影| 一卡2卡三卡四卡精品乱码亚洲| 久久国产乱子免费精品| 蜜桃亚洲精品一区二区三区| 久久人人精品亚洲av| 亚洲av不卡在线观看| 黄色女人牲交| 久久久久久久亚洲中文字幕| 男女下面进入的视频免费午夜| 成人二区视频| 亚洲自拍偷在线| 精品午夜福利视频在线观看一区| 一级a爱片免费观看的视频| 成人午夜高清在线视频| 中文字幕高清在线视频| 丰满乱子伦码专区| 联通29元200g的流量卡| 久久亚洲精品不卡| 亚洲精品456在线播放app | 午夜亚洲福利在线播放| 国产探花极品一区二区| 精品免费久久久久久久清纯| netflix在线观看网站| www日本黄色视频网| 国产黄色小视频在线观看| 亚洲欧美清纯卡通| 99久久九九国产精品国产免费| 亚洲五月天丁香| 中文字幕高清在线视频| 日韩强制内射视频| 免费看光身美女| 亚洲精华国产精华精| 日韩,欧美,国产一区二区三区 | 性插视频无遮挡在线免费观看| 老司机午夜福利在线观看视频| 久久婷婷人人爽人人干人人爱| 亚洲欧美日韩高清在线视频| 搡老岳熟女国产| 午夜久久久久精精品| 精品久久久久久久久av| 国产午夜精品论理片| 天天躁日日操中文字幕| h日本视频在线播放| 啪啪无遮挡十八禁网站| av.在线天堂| 午夜激情福利司机影院| 亚洲精品成人久久久久久| 一级av片app| 男女做爰动态图高潮gif福利片| 老师上课跳d突然被开到最大视频| 亚洲av成人av| 国产色爽女视频免费观看| 欧美一区二区国产精品久久精品| 欧美+日韩+精品| 天天躁日日操中文字幕| 欧美最新免费一区二区三区| 亚洲av熟女| 999久久久精品免费观看国产| 很黄的视频免费| 日本色播在线视频| 老女人水多毛片| 日韩强制内射视频| 国产免费男女视频| 99riav亚洲国产免费| 成人国产一区最新在线观看| 女的被弄到高潮叫床怎么办 | 欧美最黄视频在线播放免费| 女生性感内裤真人,穿戴方法视频| 亚洲性久久影院| 久久久久国内视频| 午夜福利高清视频| 女的被弄到高潮叫床怎么办 | 亚洲中文日韩欧美视频| 91在线精品国自产拍蜜月| 国产熟女欧美一区二区| 亚洲人与动物交配视频| 午夜福利欧美成人| 久久久久久久亚洲中文字幕| 噜噜噜噜噜久久久久久91| 亚洲国产精品合色在线| 毛片女人毛片| 亚洲乱码一区二区免费版| 中文字幕人妻熟人妻熟丝袜美| 久久精品综合一区二区三区| 欧美精品啪啪一区二区三区| 禁无遮挡网站| 午夜福利在线在线| 99热这里只有精品一区| 少妇的逼好多水| 国产精品99久久久久久久久| 毛片一级片免费看久久久久 | 精品久久国产蜜桃| 人妻制服诱惑在线中文字幕| 精品久久久久久久末码| 免费在线观看影片大全网站| 麻豆av噜噜一区二区三区| 国内揄拍国产精品人妻在线| 精品久久久噜噜| 一个人免费在线观看电影| 亚洲中文字幕日韩| 国产日本99.免费观看| 久久精品国产亚洲av香蕉五月| 天堂网av新在线| 俄罗斯特黄特色一大片| 欧美最新免费一区二区三区| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久,| 国产大屁股一区二区在线视频| 日韩强制内射视频| 又黄又爽又免费观看的视频| 性色avwww在线观看| 日日啪夜夜撸| 别揉我奶头 嗯啊视频| 国产国拍精品亚洲av在线观看| 韩国av一区二区三区四区| a在线观看视频网站| 亚洲av中文av极速乱 | 九九久久精品国产亚洲av麻豆| 真人做人爱边吃奶动态| 亚洲黑人精品在线|