白 蕓24，康宏亮，王文龍26，速 歡，李建明248，馬春艷

排土場(chǎng)平臺(tái)-邊坡系統(tǒng)溝蝕形態(tài)演變與產(chǎn)沙特征

白 蕓1,2,3,4，康宏亮5，王文龍1,2,6※，速 歡7，李建明1,2,4,8，馬春艷3

（1. 中國(guó)科學(xué)院教育部水土保持與生態(tài)環(huán)境研究中心，楊凌 712100；2. 中國(guó)科學(xué)院水利部水土保持研究所，楊凌 712100；3. 榆林學(xué)院陜西省陜北礦區(qū)生態(tài)修復(fù)重點(diǎn)實(shí)驗(yàn)室，榆林 719000；4. 中國(guó)科學(xué)院大學(xué)，北京 100049；5. 長(zhǎng)安大學(xué)土地工程學(xué)院，西安 710054；6. 西北農(nóng)林科技大學(xué)黃土高原土壤侵蝕與旱地農(nóng)業(yè)國(guó)家重點(diǎn)實(shí)驗(yàn)室，楊凌 712100；7. 南京市水利規(guī)劃設(shè)計(jì)院股份有限公司，南京 210000；8. 長(zhǎng)江科學(xué)院水土保持研究所，武漢 430010）

排土場(chǎng)是露天采礦區(qū)主要泥沙來(lái)源的人造地貌之一，控制排土場(chǎng)土壤侵蝕對(duì)礦區(qū)高質(zhì)量發(fā)展具有重要意義。該研究采用野外放水沖刷試驗(yàn)，研究排土場(chǎng)平臺(tái)-邊坡系統(tǒng)溝蝕演變及產(chǎn)沙特征。結(jié)果表明：1）溝蝕演變中的主導(dǎo)侵蝕方式存在階段性轉(zhuǎn)變，侵蝕溝形態(tài)發(fā)育特征也呈階段性差異，平臺(tái)侵蝕溝分為溝頭形成階段，溯源-拓寬階段和穩(wěn)定階段3個(gè)發(fā)育演變階段；邊坡侵蝕溝依次經(jīng)歷覆土層下切階段，覆土層拓寬階段，紅土層下切階段和侵蝕減緩階段。2）邊坡是平臺(tái)-邊坡系統(tǒng)主要沙源，其累積產(chǎn)沙量占平臺(tái)-邊坡系統(tǒng)的88.15%～90.16%；覆土層下切階段和紅土層下切階段是邊坡的主要產(chǎn)沙階段，其累積產(chǎn)沙量分別占邊坡的29.72%～53.36%和19.06%～48.88%。3）平臺(tái)和邊坡侵蝕速率均與徑流功率的響應(yīng)關(guān)系較優(yōu)，在平臺(tái)溝蝕的溯源-拓寬階段和穩(wěn)定階段為線(xiàn)性響應(yīng)；在邊坡溝蝕的覆土層下切階段和紅土層下切階段為指數(shù)響應(yīng)，拓寬階段和侵蝕減緩階段為線(xiàn)性響應(yīng)。模型建立中需進(jìn)一步考慮此種響應(yīng)規(guī)律隨溝蝕發(fā)育演變產(chǎn)生的變化，研究結(jié)果可為排土場(chǎng)水土保持措施布設(shè)和科學(xué)認(rèn)識(shí)溝蝕過(guò)程與特征提供參考。

土壤；沖刷試驗(yàn)；溝蝕；徑流功率；侵蝕演變；排土場(chǎng)

0 引 言

露天開(kāi)采是晉陜蒙能源基地一種重要的煤炭開(kāi)采形式，形成了一種大型人工重塑地貌——排土場(chǎng)，其以“平臺(tái)-邊坡”為基本地貌單元，單座排土場(chǎng)呈一級(jí)至多級(jí)地貌單元疊置的形態(tài)。排土場(chǎng)平臺(tái)是平臺(tái)-邊坡系統(tǒng)主要的匯水區(qū)域[1]，在平臺(tái)匯水作用下，平臺(tái)-邊坡系統(tǒng)往往遭受強(qiáng)烈的溝蝕[2-3]，使得排土場(chǎng)成為露天礦區(qū)主要的泥沙策源地。有研究顯示排土場(chǎng)平臺(tái)-邊坡系統(tǒng)土壤侵蝕是水力和重力共同作用的結(jié)果[4]，在暴雨條件下直接威脅邊坡穩(wěn)定和安全[5]，嚴(yán)重影響礦區(qū)生態(tài)安全和地區(qū)經(jīng)濟(jì)的綠色發(fā)展。

侵蝕溝是在集中股流侵蝕形成的線(xiàn)狀延伸凹地，其發(fā)生發(fā)展具有明顯的階段性特征。朱顯謨[6]根據(jù)溝蝕的發(fā)展階段、演變時(shí)期和侵蝕強(qiáng)度將其分為細(xì)溝侵蝕、淺溝侵蝕和切溝侵蝕，其中淺溝的發(fā)育經(jīng)歷不連續(xù)溝頭形成、不連續(xù)跌坑形成、連續(xù)跌坑、溝道連通等發(fā)育階段[7]；淺溝發(fā)育的不同階段產(chǎn)沙特征也不相同，郭軍權(quán)等[8]的研究結(jié)果表明淺溝系統(tǒng)產(chǎn)沙率、含沙量和總產(chǎn)沙量隨時(shí)間均出現(xiàn)先增大—波動(dòng)減小—穩(wěn)定趨勢(shì)。坡面在經(jīng)歷片蝕向細(xì)溝演變階段、細(xì)溝發(fā)育階段、細(xì)溝向切溝演變階段和切溝發(fā)育階段后，發(fā)育形成切溝；坡度的增加和上方含沙水流的匯入增加侵蝕溝內(nèi)徑流剪切力、水流功率、單位水流功率、過(guò)水?dāng)嗝鎲挝荒芰?，加快坡面侵蝕方式演變的速度與進(jìn)程[9]；溝頭形成、溝頭溯源、溝底下切和溝岸擴(kuò)張各個(gè)子過(guò)程相繼發(fā)生[10]，坡面侵蝕產(chǎn)沙強(qiáng)度變化的實(shí)質(zhì)是侵蝕方式的演變。溝蝕的發(fā)育演化過(guò)程和特征研究已獲得堅(jiān)實(shí)的基礎(chǔ)，為進(jìn)一步認(rèn)識(shí)溝蝕機(jī)制與布設(shè)溝蝕防治戰(zhàn)略提供強(qiáng)有力的理論支撐。排土場(chǎng)平臺(tái)的巨大匯水進(jìn)入邊坡后[11]，勢(shì)必加快邊坡侵蝕溝的發(fā)育演變。速歡等[12]認(rèn)為排土場(chǎng)平臺(tái)-邊坡系統(tǒng)經(jīng)歷片蝕、細(xì)溝侵蝕和切溝侵蝕3個(gè)侵蝕階段；Bai等[13]按照侵蝕溝形態(tài)特征和主導(dǎo)侵蝕方式進(jìn)一步將平臺(tái)和邊坡溝蝕階段進(jìn)行劃分，并建立了不同侵蝕階段溝深和溝寬與水動(dòng)力學(xué)參數(shù)間關(guān)系。但排土場(chǎng)平臺(tái)-邊坡系統(tǒng)侵蝕演變的不同階段內(nèi)侵蝕產(chǎn)沙過(guò)程及特征的變化尚不清晰，仍需進(jìn)一步深入探索。

坡面侵蝕產(chǎn)沙與水力條件關(guān)系緊密。已有研究表明，輸沙速率與徑流率、雷諾數(shù)、弗汝德數(shù)、單位徑流功率、過(guò)水?dāng)嗝鎲挝荒芰康葏?shù)均存在顯著相關(guān)關(guān)系[14-17]，其中針對(duì)基于WEPP（Water Erosion Prediction Project）模型的徑流剪切力和基于GUEST（Griffith University Erosion System Template）模型的徑流功率在侵蝕產(chǎn)沙預(yù)測(cè)的研究較為深入。徑流剪切力可以較好地預(yù)測(cè)不同尺度侵蝕溝的泥沙輸移特征，在小區(qū)或坡面的細(xì)溝尺度上，輸沙速率與徑流剪切力間關(guān)系多為線(xiàn)性關(guān)系和冪函數(shù)關(guān)系，在不同研究區(qū)域間均取得較為一致的研究結(jié)果[17-19]。部分學(xué)者針對(duì)淺溝侵蝕的研究結(jié)果也發(fā)現(xiàn)，徑流剪切力與淺溝系統(tǒng)的輸沙量間存在最佳的相關(guān)性[20-21]。Istanbulluoglu等[22]利用基于徑流剪切力的模型，驗(yàn)證了模型在切溝尺度上輸沙率的空間變化，結(jié)果表明模型可以反映83%的輸沙率空間變異性；Vanwalleghem等[23]建立了侵蝕產(chǎn)沙與徑流剪切力間關(guān)系，野外觀(guān)測(cè)數(shù)據(jù)的驗(yàn)證結(jié)果顯示關(guān)系式在從細(xì)溝到切溝尺度的線(xiàn)性侵蝕面上均有較好的預(yù)測(cè)結(jié)果。徑流功率在細(xì)溝和淺溝侵蝕的侵蝕輸沙量預(yù)測(cè)方面也得到廣泛應(yīng)用[24-25]。不同尺度上針對(duì)徑流剪切力和徑流功率的研究結(jié)果對(duì)深度認(rèn)識(shí)土壤侵蝕的水動(dòng)力機(jī)制和實(shí)踐侵蝕產(chǎn)沙預(yù)測(cè)具有重要意義。排土場(chǎng)平臺(tái)-邊坡系統(tǒng)作為露天礦開(kāi)采中的主要泥沙來(lái)源地，其侵蝕水動(dòng)力學(xué)機(jī)制尚不清晰，侵蝕產(chǎn)沙對(duì)水動(dòng)力學(xué)參數(shù)變化的響應(yīng)在溝蝕演變過(guò)程中的轉(zhuǎn)變特征仍需進(jìn)一步深入研究。

綜上，排土場(chǎng)平臺(tái)-邊坡系統(tǒng)的徑流條件和產(chǎn)沙特征對(duì)溝蝕演變存在怎樣的響應(yīng)，產(chǎn)沙水動(dòng)力過(guò)程是否對(duì)侵蝕溝發(fā)育演變存在某種響應(yīng)特征是關(guān)鍵科學(xué)問(wèn)題。為此，本文以排土場(chǎng)平臺(tái)-邊坡系統(tǒng)為研究對(duì)象，基于野外調(diào)查數(shù)據(jù)，構(gòu)建排土場(chǎng)平臺(tái)-邊坡系統(tǒng)徑流小區(qū)，采用人工放水沖刷試驗(yàn)，研究平臺(tái)匯水條件下平臺(tái)-邊坡系統(tǒng)溝蝕的發(fā)生發(fā)展過(guò)程，厘清平臺(tái)-邊坡系統(tǒng)侵蝕溝發(fā)育演變的過(guò)程性特征，分析侵蝕溝發(fā)育演變過(guò)程中的徑流特征和產(chǎn)沙響應(yīng)，探索侵蝕溝發(fā)育演變過(guò)程中的產(chǎn)沙水動(dòng)力學(xué)規(guī)律，以期為能源區(qū)排土場(chǎng)土壤侵蝕防治提供研究支撐。

1 材料與方法

1.1 研究區(qū)概況

試驗(yàn)在陜西省榆林市花園溝水土流失監(jiān)測(cè)站進(jìn)行，該試驗(yàn)站位于中國(guó)陜西省榆林市。野外放水沖刷試驗(yàn)前，共調(diào)查神府-東勝煤田范圍內(nèi)27個(gè)露天煤礦的80座排土場(chǎng)，調(diào)查區(qū)位于晉陜蒙三省交界區(qū)，地處毛烏素沙地與黃土高原過(guò)渡帶，為干旱半干旱大陸性季風(fēng)氣候，年平均降水量約為348～400 mm，主要發(fā)生在7－9月，多為短歷時(shí)高強(qiáng)度暴雨。排土場(chǎng)由“平臺(tái)+邊坡”地形單元組成，經(jīng)調(diào)查邊坡坡度23.5°～45.0°（平均35°），邊坡坡長(zhǎng)10.0～40.0 m（平均22.5 m），大面積平臺(tái)匯水傾瀉進(jìn)入邊坡后，平臺(tái)與邊坡轉(zhuǎn)折處向下3～8 m的邊坡段溝蝕最為嚴(yán)重，因此將此坡段作為邊坡段研究對(duì)象。排土場(chǎng)土層分為表層覆土層和下伏土層，覆土層厚20～30 cm，下伏土層主要為第三紀(jì)紅土（又名三趾馬紅土或保德紅土），該下伏土層類(lèi)型占被調(diào)查排土場(chǎng)的62.5%，其主要礦物為石英和粘土礦物，伴有少量長(zhǎng)石、方解石和重礦物，主要植被有藜（L.）和紫花苜蓿（L.）。

1.2 試驗(yàn)設(shè)計(jì)與裝置

試驗(yàn)流量設(shè)計(jì)依據(jù)張漢雄等[26]對(duì)1 710個(gè)典型強(qiáng)降水事件的統(tǒng)計(jì)分析，得到黃土高原平均降雨強(qiáng)度與降雨持續(xù)時(shí)間的方程為依據(jù)，以5年一遇暴雨在120、145、170 m2平臺(tái)匯水面積上匯集的徑流量，徑流系數(shù)取0.7，經(jīng)計(jì)算后數(shù)值取整，對(duì)應(yīng)的放水流量分別為60、72和84 L/min。依據(jù)前期調(diào)查中排土場(chǎng)“平臺(tái)+邊坡”的地形和侵蝕特征，平臺(tái)-邊坡系統(tǒng)徑流小區(qū)由平臺(tái)和邊坡兩部分組成，邊坡坡度和坡長(zhǎng)為35°和6 m，平臺(tái)坡度和坡長(zhǎng)為2°和2 m，寬1.5 m（圖1）。

圖1 徑流小區(qū)及試驗(yàn)設(shè)備示意圖

徑流小區(qū)用高60 cm和厚2 mm的鋼板插入地下50 cm作為邊界，以防止水分的側(cè)向交換。試驗(yàn)中所用覆土層和第三紀(jì)紅土層土樣均采自西灣露天礦排土場(chǎng)（38°37′N(xiāo)，109°57′E），使用Mastersizer 3000 激光粒度分析儀（Malvern, UK, ±1%）測(cè)定的覆土層土樣的砂粒、粉粒和黏粒的含量分別為67.81%、17.55%和14.64%，根據(jù)國(guó)際土壤分類(lèi)屬于砂壤土，第三紀(jì)紅土土樣的砂粒、粉粒和黏粒含量分別為42.72%、33.09%和24.19%，為黏壤土。共分8層裝填試驗(yàn)用土，每層10 cm，底部6層裝填黏壤土，上部2層裝填砂壤土。覆土層和紅土層的平臺(tái)容重控制在1.40～1.45和1.50～1.60 g/cm3，邊坡容重控制在1.25～1.35和1.40～1.50 g/cm3，填土?xí)r覆土層和紅土層的土壤含水量分別控制在11.00%～14.00%和16.00%～19.00%；裝填過(guò)程中用工具扒松兩層土間接觸面以保證土層間土壤顆粒的自然接觸，裝填完成后自然沉降2周。試驗(yàn)前將表層土修整為微內(nèi)凹，以保證侵蝕溝從小區(qū)中間開(kāi)始發(fā)育。

試驗(yàn)裝置由供水設(shè)備、溢流槽、集流槽、攝像系統(tǒng)等組成（圖1）。蓄水池（長(zhǎng)×寬×高為3 m×3 m×1 m）和潛水泵位于小區(qū)上方的平臺(tái)，通過(guò)水閥和電磁流量計(jì)調(diào)節(jié)流量。溢流槽（長(zhǎng)×寬×高為1.5 m×0.5 m×0.3 m）安裝在小區(qū)平臺(tái)段頂部，通過(guò)一塊與小區(qū)同寬，長(zhǎng)15 cm的防水布與小區(qū)連接。小區(qū)邊坡段出口處安裝集流槽，采集徑流和泥沙樣品。在垂直于徑流小區(qū)平臺(tái)和邊坡正上方安裝兩個(gè)攝像機(jī)，記錄平臺(tái)和邊坡侵蝕溝的實(shí)時(shí)形態(tài)變化，平臺(tái)處架設(shè)榮耀V20，視頻錄制像素為4 K（3 840× 2 160像素，30幀/s），邊坡處架設(shè)羅技C920 PRO攝像頭，視頻錄制像素為720 P（1 280×720像素，30幀/s）。

試驗(yàn)前，校準(zhǔn)電磁流量計(jì)，調(diào)節(jié)水閥將流量調(diào)節(jié)至設(shè)計(jì)流量，試驗(yàn)流量與設(shè)計(jì)流量的相對(duì)誤差低于5%。試驗(yàn)開(kāi)始后，采用自制接樣裝置在平臺(tái)和邊坡轉(zhuǎn)折處接取平臺(tái)段徑流泥沙樣，用于測(cè)定平臺(tái)徑流含沙量，同時(shí)在集流槽出口接取徑流泥沙樣，前5 min每1 min接一次樣品，5 min后每2 min接一次。接取徑流泥沙樣的同時(shí)測(cè)定徑流寬和徑流流速，邊坡上每隔2 m設(shè)置一個(gè)觀(guān)測(cè)斷面，共3個(gè)斷面，平臺(tái)上在距離溢流口1～2 m處設(shè)置一個(gè)觀(guān)測(cè)斷面，流速采用顏色蹤法（KMnO4）測(cè)定，重復(fù)2次；流寬在各個(gè)觀(guān)測(cè)斷面的三等分點(diǎn)上定點(diǎn)觀(guān)測(cè)，重復(fù)2次。3 min后每2 min測(cè)量一次侵蝕溝深度，平臺(tái)每間隔0.25 m設(shè)置一個(gè)溝深觀(guān)測(cè)斷面，邊坡每隔0.5m設(shè)置一個(gè)，共計(jì)19個(gè)斷面，試驗(yàn)時(shí)長(zhǎng)45 min。試驗(yàn)結(jié)束后，截取錄制視頻中溝深測(cè)量時(shí)刻的圖片，將截圖導(dǎo)入Digimizer 2.0，利用放置在小區(qū)內(nèi)邊界的彩色軟尺確定1.0 m標(biāo)尺，量取各溝深觀(guān)測(cè)斷面上對(duì)應(yīng)的溝寬，將連通溝道中部的平滑曲線(xiàn)長(zhǎng)度作為溝長(zhǎng)。

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

平臺(tái)剝蝕速率為平臺(tái)段單位時(shí)間單位面積徑流剝蝕土壤的質(zhì)量，計(jì)算式為

式中M為平臺(tái)侵蝕速率，g/(m2·s)；S為平臺(tái)產(chǎn)沙率，g/s，如式（2）；D為平臺(tái)溝長(zhǎng)，m；W為平臺(tái)侵蝕溝平均寬度，m。

式中為試驗(yàn)中實(shí)際平均流量（每5～10 min記錄一次），L/min；m為平臺(tái)和邊坡轉(zhuǎn)折處徑流泥沙樣的泥沙量，g；V為平臺(tái)和邊坡轉(zhuǎn)折處徑流泥沙樣的渾水體積，mL。

邊坡侵蝕速率為邊坡段單位時(shí)間單位面積徑流剝蝕土壤的質(zhì)量，計(jì)算式為

式中M為邊坡侵蝕速率，g/(m2·s)；為集流槽出口處接樣時(shí)間內(nèi)的泥沙量，為徑流泥沙量烘干后稱(chēng)量數(shù)值，g；為接樣時(shí)間，s；D為邊坡各斷面的侵蝕溝長(zhǎng)度，m；W為邊坡各斷面侵蝕溝平均寬度，m。

徑流剪切力，用于表征水流作用于土壤顆粒的剪切應(yīng)力，計(jì)算式為式（4）～（7）[27]：

其中

式中為徑流剪切力，N/m2；γ為渾水密度，kg/m3；為水力半徑，m；為水力坡度，約為坡度的正弦值；為徑流含沙量，kg/m3；為清水密度，kg/m3；γ為泥沙顆粒密度，kg/m3；為水流寬度，m；為水流深度，m；為接樣時(shí)間內(nèi)的徑流體積，m3；為徑流流速，m/s。

徑流功率是單位面積水體勢(shì)能隨時(shí)間的變化率，用于表征作用于單位面積的水流所消耗的功率，N/(m·s)，計(jì)算式為[28]

2 結(jié)果與分析

2.1 平臺(tái)-邊坡侵蝕溝發(fā)育演變過(guò)程

圖2為平臺(tái)和邊坡侵蝕溝發(fā)育演變的形態(tài)變化過(guò)程。0～2 min平臺(tái)-邊坡系統(tǒng)的侵蝕經(jīng)歷短暫的層狀面蝕-不連續(xù)跌坎-連續(xù)跌坎的發(fā)育過(guò)程，在3 min時(shí)形成較為清晰的溝沿形態(tài)，滿(mǎn)足通過(guò)垂直于小區(qū)上方的攝像資料進(jìn)行后期溝寬和溝長(zhǎng)參數(shù)的提取條件，因此侵蝕溝形態(tài)變化的過(guò)程分析從3 min開(kāi)始。60、72、84 L/min流量下的3～45 min平臺(tái)溝長(zhǎng)為0.13～3.30、0.15～3.18和0.12～3.26 m，隨試驗(yàn)歷時(shí)呈現(xiàn)緩慢增大、快速增大和趨于穩(wěn)定的變化過(guò)程；平臺(tái)溝寬為0.11～0.26、0.13～0.25和0.15～0.32 m，隨試驗(yàn)歷時(shí)均呈現(xiàn)波動(dòng)增大、快速增大和趨于穩(wěn)定的變化；平臺(tái)溝深為0.13～0.20、0.07～0.20和0.16～0.20 m，隨試驗(yàn)歷時(shí)在3～5 min內(nèi)迅速增大到0.2 m，之后保持穩(wěn)定。根據(jù)平臺(tái)侵蝕溝在長(zhǎng)、寬、深變化過(guò)程，將平臺(tái)侵蝕劃分為3個(gè)階段：第一階段（3～5 min）溝深和溝寬迅速增大，溝長(zhǎng)緩慢增大，此階段內(nèi)平臺(tái)末端（平臺(tái)和邊坡轉(zhuǎn)折處）上層砂壤質(zhì)覆土層被不斷下切，形成溝頭，溝頭跌水高度與覆土層厚度相同，稱(chēng)為溝頭形成階段（圖3a）；第二階段在溝頭形成后，平臺(tái)溝長(zhǎng)經(jīng)歷2～4 min的緩慢增長(zhǎng)后進(jìn)入快速增長(zhǎng)時(shí)段，溝寬呈波動(dòng)增大趨勢(shì)，該階段內(nèi)溝長(zhǎng)和溝寬增量占總增長(zhǎng)值的91.84%～95.92%和68.31%～83.21%，侵蝕溝同時(shí)向長(zhǎng)和向?qū)挵l(fā)展，溝深不再變化，稱(chēng)為溯源-拓寬階段；第三階段溝長(zhǎng)、溝寬和溝深均不在發(fā)生變化，侵蝕溝溝頭到達(dá)平臺(tái)頂端與溢流槽相接處，平臺(tái)侵蝕溝形態(tài)發(fā)育趨于穩(wěn)定，稱(chēng)為穩(wěn)定階段。放水流量對(duì)平臺(tái)侵蝕溝的發(fā)育演變影響較小，不同流量下各個(gè)發(fā)育階段的持續(xù)時(shí)長(zhǎng)無(wú)顯著性差異。

注：平臺(tái)系列中P1、P2、P3分別代表溝頭形成階段、溯源-拓寬階段和穩(wěn)定階段；邊坡系列圖中S1、S2、S3、S4分別代表覆土層下切階段、覆土層拓寬階段、紅土層下切階段和侵蝕減緩階段，下同。

圖3 排土場(chǎng)平臺(tái)和邊坡侵蝕溝發(fā)育階段圖

60、72、84 L/min流量下3～45 min的邊坡溝寬為0.09～0.25、0.10～0.31和0.09～0.27 m，隨試驗(yàn)歷時(shí)均呈現(xiàn)快速增大、階梯式增大和趨于穩(wěn)定的變化；3～45 min的邊坡溝深為0.09～0.37、0.11～0.39和0.15～0.44 m，隨試驗(yàn)歷時(shí)呈快速增大、緩慢增大、二次快速增大、趨于穩(wěn)定的趨勢(shì)（圖2）。在試驗(yàn)歷時(shí)內(nèi)，邊坡溝寬和溝深的發(fā)展不同步，侵蝕溝演變的主導(dǎo)發(fā)展方向呈階段性變化：不同流量下第一階段的溝深為0.20～0.23 m，溝寬為0.14～0.15 m，侵蝕溝寬深比小于1，侵蝕溝以向深發(fā)展為主要特點(diǎn)，且主要發(fā)生在上層砂壤質(zhì)覆土，稱(chēng)為覆土層下切階段（圖3b）；第二階段侵蝕溝深度為0.22～0.25 m，溝寬為0.20～0.22 m，與覆土層下切階段相比溝深增大9.00%～24.47%，溝寬增大38.47%～46.90%，侵蝕溝以向?qū)挵l(fā)展為主，且溝壁拓寬主要發(fā)生在覆土層，稱(chēng)為覆土層拓寬階段；第三階段溝深再次快速增大，溝寬呈小幅階梯式增大并趨于穩(wěn)定，較覆土層拓寬階段的溝深增大47.05%～75.10%，溝寬增大10.51%～40.26%，侵蝕溝以向深發(fā)展為主，且主要發(fā)生在第三紀(jì)紅土層，稱(chēng)為紅土層下切階段；第四階段溝深和溝寬均緩慢增大，較紅土層下切階段溝深和溝寬分別增大2.75%～12.42%和2.60%～9.13%，溝蝕速度減緩，稱(chēng)為侵蝕減緩階段。侵蝕溝演變速率受流量影響較大，隨流量增大，侵蝕階段的轉(zhuǎn)變時(shí)間提前（圖2），侵蝕溝在更短的時(shí)間內(nèi)發(fā)育形成較大的溝深和溝寬。

2.2 平臺(tái)-邊坡侵蝕產(chǎn)沙的階段性特征

圖4為平臺(tái)和邊坡侵蝕速率隨試驗(yàn)歷時(shí)的變化過(guò)程。

圖4 不同流量下平臺(tái)和邊坡侵蝕速率隨試驗(yàn)歷時(shí)的變化

在平臺(tái)侵蝕溝3個(gè)不同發(fā)育階段，各試驗(yàn)流量的平臺(tái)侵蝕速率分別呈穩(wěn)定、上升和下降的變化過(guò)程，各個(gè)發(fā)育階段內(nèi)的變異系數(shù)分別為0.07～0.12、0.17～0.51和0.29～0.48，階段平均侵蝕速率隨侵蝕溝發(fā)育演化先增大后減小（表1）。3個(gè)階段的累積產(chǎn)沙貢獻(xiàn)比為6.95%～8.21%、39.98%～74.38%和18.67%～52.15%。溯源-拓寬階段是主要的侵蝕產(chǎn)沙階段，侵蝕速率和累計(jì)產(chǎn)沙量貢獻(xiàn)比均最大，分別為104.14～175.69 g/(m2·s)和39.98%～74.38%。

在邊坡4個(gè)侵蝕溝發(fā)育演變階段，邊坡侵蝕速率分別呈現(xiàn)急劇上升后出現(xiàn)拐點(diǎn)、急劇下降、波動(dòng)持平、基本持平的變化特征（圖4）。不同流量下，階段平均侵蝕速率隨侵蝕溝發(fā)育演化存在下降趨勢(shì)（表1），以第一階段為基準(zhǔn)，后3個(gè)發(fā)育階段的平均侵蝕速率分別減小42.25%～68.32%、39.98%～74.61%和79.91%～83.46%。4個(gè)發(fā)育階段的累積產(chǎn)沙量逐漸減小，累積產(chǎn)沙量貢獻(xiàn)比分別為29.72%～53.36%、11.41%～21.60%、19.06%～48.88%和2.64%～20.04%（表1），覆土層下切階段和紅土層下切階段是邊坡主要的產(chǎn)沙階段。邊坡侵蝕速率分別是平臺(tái)的0.21～67.20倍，試驗(yàn)結(jié)束時(shí)邊坡累積產(chǎn)沙量占平臺(tái)-邊坡系統(tǒng)的88.15%～90.16%，是平臺(tái)的7.44～9.16倍，單位坡長(zhǎng)邊坡累積產(chǎn)沙量是平臺(tái)的2.48～3.05倍，邊坡是平臺(tái)-邊坡系統(tǒng)的主要產(chǎn)沙來(lái)源。

表1 平臺(tái)和邊坡不同侵蝕溝發(fā)育階段平均侵蝕速率和累積產(chǎn)沙量

2.3 平臺(tái)-邊坡徑流的階段性特征

表2為不同侵蝕溝發(fā)育階段平臺(tái)徑流指標(biāo)值。分析中將前2 min的各徑流指標(biāo)的平均值作為初始值，以更清晰的說(shuō)明侵蝕溝發(fā)育演化中徑流的階段性變化特征。平臺(tái)不同發(fā)育階段平均徑流流速隨侵蝕溝發(fā)育演化呈先減小后增大的變化過(guò)程，最小值較初始值減小?0.51%～21.12%，穩(wěn)定階段的平均流速較初始值增大?15.75%～58.31%。徑流剪切力和徑流功率隨侵蝕溝發(fā)育演化變逐漸增大，與初始值相比，3個(gè)發(fā)育階段的徑流剪切力分別增大0.11～0.31倍、1.10～1.94倍和1.86～3.21倍，徑流功率分別增大?0.04～0.14倍、1.15～3.25倍和1.47～6.15倍，表明在平臺(tái)溝蝕發(fā)育的后兩個(gè)階段，徑流對(duì)土壤顆粒的分離和搬運(yùn)能力增加。平臺(tái)徑流一旦進(jìn)入邊坡，流速、徑流剪切力和水流功率均增大，不同流量下邊坡徑流剪切力較平臺(tái)分別增大7.11～96.93倍、7.55～120.86倍和12.21～83.03倍，邊坡水流功率較平臺(tái)增大14.69～116.70倍、7.59～239.59倍和21.30～130.91倍，徑流對(duì)土壤顆粒的剝離和輸移能力急劇增大。

表2 不同侵蝕溝發(fā)育階段平臺(tái)徑流指標(biāo)值

注：表中P0代表平臺(tái)徑流指標(biāo)的初始值。

Note: P0represents the initial values of runoff index on the platform.

表3為不同侵蝕溝發(fā)育階段邊坡徑流指標(biāo)值。不同流量下，階段平均邊坡徑流流速隨侵蝕溝發(fā)育演化呈下降趨勢(shì)，四個(gè)發(fā)育階段的流速分別較初始值減小7.92%～17.50%、26.66%～32.18%、31.03%～62.38%和43.68%～60.00%，四個(gè)階段的流速下降幅度占比分別為14.28%～31.58%、15.28%～42.11%、0～60.71%和0～29.17%。整體上，徑流剪切力隨侵蝕溝發(fā)育演化呈波動(dòng)增大的變化，徑流功率呈先增大后減小的變化特點(diǎn)。在覆土層下切階段，徑流剪切力和徑流功率較初始值的增大顯著，不同流量下徑流剪切力分別較初始值增大4.57倍、1.05倍和1.54倍，徑流功率分別增大3.12倍、1.19倍和0.86倍，增大幅度隨流量增大而減??；覆土層下切階段之后，徑流剪切力呈波動(dòng)式增大，穩(wěn)定階段的值較覆土層下切階段增大3.22%～38.38%。徑流功率則逐漸減小，穩(wěn)定階段的徑流功率較覆土層下切階段減小19.07%～48.84%；說(shuō)明在溝蝕發(fā)育后期徑流分離土壤顆粒能力并無(wú)顯著減小，但單位面積徑流的能量消耗明顯減小。分析顯示不同階段間徑流剪切力和水流功率均值無(wú)顯著性差異。隨放水流量增加，徑流剪切力和徑流功率隨之增大（<0.05），徑流流速在前兩發(fā)育階段也隨之增大，但在后兩發(fā)育階段無(wú)增大趨勢(shì)。

表3 不同侵蝕溝發(fā)育階段邊坡徑流指標(biāo)值

注：表中S0代表邊坡徑流指標(biāo)的初始值。

Note: S0represents the initial values of runoff index on the steep slope.

2.4 平臺(tái)和邊坡侵蝕速率與水動(dòng)力參數(shù)關(guān)系

圖5為平臺(tái)侵蝕速率與徑流剪切力和徑流功率的關(guān)系分析圖。平臺(tái)侵蝕速率與徑流剪切力和徑流功率在溯源-拓寬階段（P2）和穩(wěn)定階段（P3）存在顯著線(xiàn)性正相關(guān)（<0.01），而在溝頭形成階段相關(guān)性差。由回歸關(guān)系獲得的可蝕性參數(shù)在溯源-拓寬階段到穩(wěn)定階段呈增大趨勢(shì)。溯源拓寬階段和穩(wěn)定階段的臨界徑流剪切力為-1.43和1.31 N/m2，臨界徑流功率為-0.74和0.78 N/(m·s)。

圖5 平臺(tái)侵蝕速率與徑流剪切力和徑流功率關(guān)系

圖6為邊坡侵蝕速率與徑流剪切力和徑流功率的關(guān)系分析圖。覆土層下切階段（S1）和紅土層下切階段（S3）的侵蝕速率與徑流剪切力和徑流功率間指數(shù)關(guān)系顯著（<0.01），覆土層拓寬階段和侵蝕減緩階段的侵蝕速率與徑流剪切力和徑流功率存在顯著線(xiàn)性關(guān)系，且侵蝕速率與徑流功率間相關(guān)性?xún)?yōu)于徑流剪切力。與覆土層下切階段相比，紅土層下切階段的指數(shù)函數(shù)前因子和指數(shù)前常數(shù)均減小，表明覆土層下切階段的侵蝕速率對(duì)徑流剪切力和徑流功率的響應(yīng)更敏感，單位徑流剪切力和徑流功率的變化導(dǎo)致更大幅度的侵蝕速率變化。與覆土層拓寬階段相比，侵蝕減緩階段的可蝕性參數(shù)減小，臨界徑流剪切力和臨界徑流功率增大，同樣表明覆土層拓寬階段的侵蝕速率對(duì)徑流剪切力和徑流功率的響應(yīng)更敏感。

圖6 邊坡侵蝕速率與徑流剪切力和徑流功率關(guān)系

3 討 論

3.1 排土場(chǎng)平臺(tái)-邊坡系統(tǒng)侵蝕溝發(fā)育的階段性差異

平臺(tái)侵蝕溝發(fā)育分為溝頭形成階段，溯源-拓寬階段和穩(wěn)定階段3個(gè)演變階段，邊坡溝蝕演變過(guò)程依次為覆土層下切階段，覆土層拓寬階段，紅土層下切階段和侵蝕減緩階段，由于平臺(tái)段和邊坡段地形、土壤抗蝕性等因素的差異，溝蝕發(fā)育演變特征不同步。與平臺(tái)相比，邊坡土壤抗蝕性較小，坡度較陡，徑流在緊臨地形轉(zhuǎn)折點(diǎn)的邊坡段迅速形成侵蝕凹地。一方面，侵蝕凹地向平臺(tái)段發(fā)展，下切形成溝頭，當(dāng)下切侵蝕至紅土層時(shí)土壤抗蝕性增大，限制溝頭繼續(xù)下切[29]，一旦形成溝頭，則在溝頭流的作用下溯源侵蝕[7]，并與側(cè)向侵蝕耦合[30]，侵蝕溝進(jìn)入溯源-拓寬階段，當(dāng)溝頭溯源至溢流槽與平臺(tái)段連接處時(shí)，侵蝕溝進(jìn)入穩(wěn)定階段（圖3a）。另一方面，侵蝕凹地加速了邊坡徑流的匯聚，加之邊坡徑流沿程加速，在覆土層上形成以下切侵蝕為主導(dǎo)的侵蝕階段，即覆土層下切階段；當(dāng)下切侵蝕至紅土層時(shí)，紅土層的抗蝕性較覆土層增大，侵蝕溝深度增長(zhǎng)較慢，延長(zhǎng)了徑流對(duì)覆土層底部的淘蝕時(shí)間，形成利于溝壁崩塌的臨空面進(jìn)而促進(jìn)溝壁的拓寬[31]，侵蝕溝呈現(xiàn)以側(cè)向侵蝕為主導(dǎo)侵蝕方式的發(fā)育演變特征（圖3b），即為覆土層拓寬階段；當(dāng)徑流完全集中于紅土層床面后，在次級(jí)溝頭的不斷溯源侵蝕和跌水沖淘作用下[32]，表現(xiàn)出溝深重新迅速增大，在紅土層形成以下切侵蝕為主導(dǎo)侵蝕方式的侵蝕階段，為紅土層下切階段；最后，侵蝕溝形態(tài)和徑流能量間通過(guò)互饋?zhàn)饔眠_(dá)到均衡，溝蝕發(fā)育速度減緩甚至達(dá)到穩(wěn)定階段。平臺(tái)-邊坡系統(tǒng)侵蝕溝是在溯源侵蝕、下切侵蝕和側(cè)向侵蝕的共同作用下完成的，但主導(dǎo)侵蝕方式的階段性轉(zhuǎn)變致使侵蝕溝發(fā)育演變呈現(xiàn)階段性差異，其中平臺(tái)和邊坡轉(zhuǎn)折處的地形條件、土壤條件顯著不同，轉(zhuǎn)折處侵蝕凹地成為平臺(tái)段和邊坡段溝蝕發(fā)育的分隔點(diǎn)，呂剛等[4]基于平臺(tái)與邊坡轉(zhuǎn)折處裂隙優(yōu)先流的研究也表明平臺(tái)-邊坡系統(tǒng)地形的特殊性。試驗(yàn)開(kāi)始階段在轉(zhuǎn)折處形成的侵蝕凹地逐漸發(fā)育形成跌水，跌水的沖淘和侵蝕是邊坡侵蝕溝發(fā)育演變的重要特征之一，在本試驗(yàn)放水條件下促進(jìn)了邊坡溝底的下切侵蝕，最終形成坡面切溝[33]，與工程堆積體的細(xì)溝侵蝕[27,34]相比，造成的侵蝕強(qiáng)度更大。因此做好轉(zhuǎn)折處水土保持措施十分重要，在不均勻沉降和機(jī)械施工產(chǎn)生的地勢(shì)低凹處和具有多級(jí)平臺(tái)-邊坡地貌單元的排土場(chǎng)的道路連接處應(yīng)重點(diǎn)防護(hù)，通過(guò)配備合理的工程措施和林草措施，以減緩和防止排土場(chǎng)平臺(tái)-邊坡系統(tǒng)的侵蝕溝發(fā)育。

3.2 排土場(chǎng)平臺(tái)-邊坡系統(tǒng)侵蝕產(chǎn)沙的階段性演變

在平臺(tái)和邊坡侵蝕溝發(fā)育演變的不同階段，侵蝕速率具有階段性變化特征，邊坡作為平臺(tái)-邊坡系統(tǒng)主要的泥沙來(lái)源，其侵蝕速率隨侵蝕溝發(fā)育演變呈減小趨勢(shì)，這與前人針對(duì)淺溝侵蝕產(chǎn)沙特征[8]和工程堆積體陡坡的侵蝕產(chǎn)沙特征的研究結(jié)果一致[34]。侵蝕速率隨溝蝕發(fā)育演變的階段性變化與土壤性質(zhì)、主導(dǎo)侵蝕方式轉(zhuǎn)變、徑流條件等因素密切相關(guān)[35]。在覆土層下切階段，侵蝕溝深度和寬度變化最為迅速，此階段主要發(fā)生在砂壤質(zhì)的表層覆土，土壤抗侵蝕能力差；同時(shí)，此階段內(nèi)徑流剪切力和水流功率較初始值迅速增大，徑流迅速集中成股流，對(duì)土壤顆粒的分離和輸移能力加劇，導(dǎo)致此階段內(nèi)邊坡侵蝕速率呈增大趨勢(shì)（圖4），該階段的侵蝕速率均值在4個(gè)發(fā)育階段中最大（表1）。在覆土層拓寬階段，隨著溝底床面的粗糙化，徑流流速較前一階段減小，但徑流剪切力和徑流功率較覆土層下切階段無(wú)顯著減?。ū?），徑流對(duì)土壤顆粒的分離和輸移能力無(wú)顯著減弱，侵蝕速率逐漸減小主要?dú)w因于主導(dǎo)侵蝕方式的轉(zhuǎn)變和土壤抗蝕性的差異，此階段的主導(dǎo)侵蝕方式是側(cè)向侵蝕，伴有紅土層的下切侵蝕，溝底股流對(duì)溝壁的側(cè)向侵蝕作用隨著下切侵蝕中徑流的進(jìn)一步集中逐漸弱減，此外受黏壤質(zhì)紅土層較大的抗蝕性和容重的影響，土壤侵蝕速率逐漸減小。紅土層下切階段的泥沙主要來(lái)源于紅土層溝床的下切侵蝕，由于黏壤質(zhì)土層較強(qiáng)的抗侵蝕能力導(dǎo)致下切速率減緩，同時(shí)溝底二級(jí)溝頭的發(fā)育[36]，在溝底形成跌水-深潭結(jié)構(gòu)（圖3b），跌水坑能耗消耗增大，徑流流速和徑流功率大幅減小（表2），土壤侵蝕速率持續(xù)減小。侵蝕減緩階段的侵蝕速率維持在較低水平，此時(shí)侵蝕溝發(fā)育十分緩慢，侵蝕溝形態(tài)無(wú)明顯變化，但床面形態(tài)更趨復(fù)雜，徑流能量持續(xù)消耗，此階段流速和徑流功率是試驗(yàn)過(guò)程中的最低值，徑流剪切力并無(wú)顯著減小，與紅土層下切階段相比，徑流剝蝕能力無(wú)顯著下降，但徑流攜帶的泥沙顆粒在跌水潭中受射流沖淘和能量消散后，部分泥沙顆粒沉積在跌水潭，輸移的泥沙顆粒被沉積過(guò)程抵消后，剩余泥沙顆粒方可輸移出小區(qū)，因此泥沙輸移量急劇減小。從平臺(tái)-邊坡系統(tǒng)溝蝕發(fā)育演變的全過(guò)程來(lái)看，其發(fā)育演變速度極快，在初始階段即表現(xiàn)為高產(chǎn)沙、高徑流能量的特點(diǎn)，相較于邊坡，平臺(tái)在大流量條件下仍表現(xiàn)為低侵蝕速率，但平臺(tái)徑流一旦進(jìn)入邊坡，徑流剪切力和徑流功率均急劇增大，邊坡溝蝕的防控就較為困難，做好平臺(tái)徑流的分控是控制排土場(chǎng)平臺(tái)-邊坡系統(tǒng)溝蝕的重要方向。

3.3 排土場(chǎng)平臺(tái)-邊坡系統(tǒng)侵蝕產(chǎn)沙水動(dòng)力學(xué)特征

在平臺(tái)和邊坡的不同侵蝕溝發(fā)育演變階段，侵蝕速率和徑流剪切力與徑流功率均存在顯著相關(guān)關(guān)系，其中徑流功率是預(yù)測(cè)各階段侵蝕速率的較優(yōu)參數(shù)，Kang等[7,20]針對(duì)淺溝的研究結(jié)果顯示最優(yōu)預(yù)測(cè)參數(shù)為徑流剪切力，且侵蝕速率與水動(dòng)力學(xué)參數(shù)間呈線(xiàn)性正相關(guān)，這可能與文獻(xiàn)中采用單場(chǎng)次試驗(yàn)的侵蝕速率與徑流指標(biāo)的平均值作為分析基礎(chǔ)有關(guān)[25]，也可能與試驗(yàn)不同坡度有關(guān)[37]。試驗(yàn)所得的邊坡覆土層拓寬階段的臨界徑流剪切力和臨界徑流功率與郭明明等[21]獲得的未耕作淺溝的臨界剪切力（17.58 N/m2）和臨界徑流功率（5.03 N/(m·s)）較為接近，但大于Xu等[20]在淺溝中的獲得的值，遠(yuǎn)大于細(xì)溝和層狀面蝕階段的臨界值，這可能是因?yàn)閯?chuàng)建淺溝水流通道所需的徑流剪切力、徑流功率和能量比細(xì)溝或?qū)訝顝搅鞯母骩20]。

邊坡侵蝕速率對(duì)徑流剪切力和徑流功率的響應(yīng)隨侵蝕溝發(fā)育演變呈現(xiàn)階段性差異，覆土層下切階段和紅土層下切階段的邊坡侵蝕速率與徑流剪切力和徑流功率間呈指數(shù)正相關(guān)，覆土層拓寬階段和侵蝕減緩階段的邊坡侵蝕速率與徑流剪切力和徑流功率間呈線(xiàn)性正相關(guān)。侵蝕產(chǎn)沙對(duì)水動(dòng)力學(xué)參數(shù)的響應(yīng)差異在由層狀面蝕向細(xì)溝侵蝕和由細(xì)溝侵蝕向淺溝侵蝕轉(zhuǎn)變中也存在，有研究表明層狀面蝕階段侵蝕速率與徑流動(dòng)能間呈線(xiàn)性相關(guān)，而在細(xì)溝侵蝕階段則存在對(duì)數(shù)相關(guān)[27]，不同侵蝕階段侵蝕速率與水動(dòng)力學(xué)參數(shù)的相關(guān)關(guān)系一致但關(guān)鍵參數(shù)也存在不同，細(xì)溝侵蝕階段的臨界徑流剪切力和臨界徑流功率大于片蝕階段[38]，淺溝侵蝕的臨界徑流剪切力和臨界徑流功率高于層狀面蝕和細(xì)溝侵蝕[39]。這種響應(yīng)特征差異與不同主導(dǎo)侵蝕方式作用下的下墊面條件變化有關(guān)。在覆土層下切階段和紅土層下切階段，侵蝕泥沙主要來(lái)源于侵蝕溝向深發(fā)展，泥沙來(lái)源充足（邊坡的2個(gè)主要產(chǎn)沙階段），侵蝕速率隨徑流侵蝕動(dòng)力增大呈指數(shù)函數(shù)增大，說(shuō)明當(dāng)徑流剪切力和徑流功率超過(guò)某一閾值后，相較于線(xiàn)性相關(guān)，侵蝕產(chǎn)沙的增幅更急劇，側(cè)面反映出在以下切侵蝕為主導(dǎo)的侵蝕階段溝蝕發(fā)育更加劇烈；紅土層下切階段的侵蝕速率與徑流剪切力和徑流功率的相關(guān)性較覆土層下切階段減小，這與跌水-深潭結(jié)構(gòu)對(duì)剪切應(yīng)力的分配和徑流動(dòng)能消耗有關(guān)，跌水-深潭結(jié)構(gòu)會(huì)使深潭顆粒粗化，粒徑較大不移動(dòng)的顆粒承擔(dān)了總剪切應(yīng)力的很大一部分，并減少了可用于運(yùn)輸較細(xì)顆粒的應(yīng)力[40]。在覆土層拓寬階段，溝壁崩塌較集中（圖3b），崩積體墜落對(duì)溝底含沙水流的擾動(dòng)和擁堵會(huì)形成阻力，且徑流沖刷床面已抵達(dá)紅土層，土壤顆粒間黏結(jié)力大[41]，致使土壤可蝕性減小，坡面泥沙供應(yīng)量減??；而在侵蝕減緩階段由于溝底跌水-深潭結(jié)構(gòu)更加發(fā)育，有效剪切應(yīng)力減小，徑流流速和徑流動(dòng)能減小，徑流的剝蝕和輸移能力減弱，同時(shí)侵蝕溝形態(tài)（溝長(zhǎng)、寬和深）趨于穩(wěn)定，供沙能力進(jìn)一步減小，因此侵蝕速率對(duì)徑流剪切力和徑流功率呈較平緩的線(xiàn)性響應(yīng)，且侵蝕減緩階段線(xiàn)性擬合式的斜率遠(yuǎn)小于覆土層拓寬階段（圖6）。由此可見(jiàn)，溝蝕輸沙模型的建立需考慮侵蝕產(chǎn)沙與水動(dòng)力學(xué)參數(shù)間響應(yīng)關(guān)系隨主導(dǎo)侵蝕方式和侵蝕演變階段的轉(zhuǎn)變，以進(jìn)一步提高模型的預(yù)測(cè)能力和適用性。

4 結(jié) 論

采用野外放水沖刷試驗(yàn)，以排土場(chǎng)平臺(tái)-邊坡系統(tǒng)為研究對(duì)象，分析平臺(tái)匯流條件下系統(tǒng)的侵蝕溝發(fā)育演變過(guò)程及過(guò)程中徑流和產(chǎn)沙特征，揭示了侵蝕演變過(guò)程中侵蝕速率對(duì)水動(dòng)力學(xué)參數(shù)間響應(yīng)規(guī)律，得出以下結(jié)論：

1）平臺(tái)-邊坡系統(tǒng)溝蝕發(fā)育演變是在溯源侵蝕、下切侵蝕、側(cè)向侵蝕三種侵蝕方式共同作用下完成的，但主導(dǎo)侵蝕方式存在階段性轉(zhuǎn)變，侵蝕溝形態(tài)發(fā)育特征也呈階段性差異。平臺(tái)侵蝕溝發(fā)育分為溝頭形成階段，溯源-拓寬階段和穩(wěn)定階段等3個(gè)發(fā)育演變階段；邊坡侵蝕溝依次經(jīng)歷覆土層下切階段，覆土層拓寬階段，紅土層下切階段和侵蝕減緩階段。

2）平臺(tái)和邊坡的產(chǎn)沙特征在不同溝蝕發(fā)育演變階段呈現(xiàn)差異，溯源-拓寬階段是平臺(tái)主要的產(chǎn)沙階段，占總產(chǎn)沙量的39.98%～74.38%；覆土層下切階段和紅土層下切階段是邊坡主要產(chǎn)沙階段，其累積產(chǎn)沙量貢獻(xiàn)比分別為29.72%～53.36%和19.06%～48.88%。邊坡是平臺(tái)-邊坡系統(tǒng)主要沙源，其累積產(chǎn)沙量占平臺(tái)-邊坡系統(tǒng)的88.15%～90.16%。

3）邊坡侵蝕速率對(duì)徑流剪切力和徑流功率的響應(yīng)在不同溝蝕演變階段間呈現(xiàn)差異，在覆土層下切階段和紅土層下切階段為指數(shù)響應(yīng)關(guān)系，在覆土層拓寬階段和侵蝕減緩階段為線(xiàn)性關(guān)系，模型建立中應(yīng)進(jìn)一步考慮此種響應(yīng)規(guī)律隨溝蝕發(fā)育演變產(chǎn)生的變化。

[1] Zhang L, Wang J M, Bai Z K, et al. Effects of vegetation on runoff and soil erosion on reclaimed land in an opencast coal-mine dump in a loess area[J]Catena, 2015, 128: 44-53.

[2] Gong C, Lei S, Bian Z F, et al. Analysis of the development of an erosion gully in an open-pit coal mine dump during a winter freeze-thaw cycle by using low-cost UAVs[J]Remote Sensing, 2019, 11(11): 1356.

[3] Lowry J B C, Narayan M, Hancock G R, et al. Understanding post-mining landforms: utilising pre-mine geomorphology to improve rehabilitation outcomes [J]Geomorphology, 2019, 328: 93-107.

[4] 呂剛，李葉鑫，寧寶寬，等. 暴雨作用下排土場(chǎng)平臺(tái)-邊坡系統(tǒng)土壤侵蝕過(guò)程模擬研究[J]煤炭學(xué)報(bào)，2021，46(5)：1463-1476.

Lv Gang, Li Yexin, Ning Baokuan, et al. Soil erosion process of platform-slope system of dump under heavy rain[J]. Journal of China Coal Society, 2021, 46(5): 1463-1476. (in Chinese with English abstract)

[5] 李葉鑫，史東梅，呂剛，等. 降雨強(qiáng)度對(duì)紫色土工程堆積體邊坡穩(wěn)定性的影響[J]土壤通報(bào)，2017，48(6)：1475-1480.

Li Yexin, Shi Dongmei, Lv Gang, et al. Effects of rainfall intensity on slope stability of purple engineering spoil[J]. Chinese Journal of Soil Science, 2017, 48(6): 1475- 1480. (in Chinese with English abstract)

[6] 朱顯謨. 黃土區(qū)土壤侵蝕的分類(lèi)[J]土壤學(xué)報(bào)，1956，4(2)：99-115.

Zhu Xianmo. Classification on the soil erosion in the loess region[J]Acta Pedologica Sinica, 1956, 4(2): 99-115. (in Chinese with English abstract)

[7] Kang H, Guo M, Wang W. Ephemeral gully erosion in concentrated flow channels induced by rainfall and upslope inflow on steep loessial slopes[J]Land Degradation & Development, 2021, 32(17): 5037-5051.

[8] 郭軍權(quán)，王文龍. 坡度對(duì)淺溝侵蝕產(chǎn)沙的野外放水沖刷試驗(yàn)影響[J]水土保持學(xué)報(bào)，2019，33(4)：87-92，212.

Guo Junquan, Wang Wenlong. Experimental study on the effects of field scouring slope on ephemeral gully erosion and sediment yield[J]Journal of Soil and Water Conservation, 2019, 33(4): 87-92, 212. (in Chinese with English abstract)

[9] 張新和. 黃土坡面片蝕-細(xì)溝侵蝕-切溝侵蝕演變與侵蝕產(chǎn)沙過(guò)程研究[D]. 楊凌：西北農(nóng)林科技大學(xué)，2007.

Zhang Xinhe. Experimental Study on Evolution Process of Sheet Erosion-Rill Erosion-Gully Erosion and Sediment Yield Process on Loess Hillslope[D]. Yangling: Northwest A & F University, 2007. (in Chinese with English abstract)

[10] Thomas J T, Iverson N R, Burkart M R, et al. Long-term growth of a valley‐bottom gully, western Lowa[J]Earth Surface Processes Landforms, 2010, 29(8): 995-1009.

[11] 胡興定，白中科，張靈，等. 黃土區(qū)大型露天礦排土場(chǎng)水力侵蝕計(jì)算與防治[J]. 水土保持研究，2017，24(5)：21-26.

Hu Xingding, Bai Zhongke, Zhang Ling, et al. Hydraulic erosion calculation and control in dumping site of lager opencast coal mine in loess area[J]. Research of Soil and Water Conservation, 2017, 24(5): 21-26. (in Chinese with English abstract)

[12] 速歡，王文龍，康宏亮，等. 露天礦排土場(chǎng)平臺(tái)-邊坡系統(tǒng)侵蝕形態(tài)及徑流產(chǎn)沙特征[J]應(yīng)用生態(tài)學(xué)報(bào)，2020，31(9)：3194-3206.

Su Huan, Wang Wenlong, Kang Hongliang, et al. Erosion morphology and the characteristics of runoff and sediment yielding in platform-slope system of opencast coal mine[J]Chinese Journal of Applied Ecology, 2020, 31(9): 3194-3206. (in Chinese with English abstract)

[13] Bai Y, Guo M, Kang H, et al. Morphodynamics of gully development on the platform–slope system of spoil dumps under platform concentrated flow[J]Land, 2021, 10(11): 1270.

[14] Aliakbar N S, Chen Q W, Shahram K, et al. Assessment of land use impact on hydraulic threshold conditions for gully head cut initiation[J]Hydrology Earth System Sciences, 2016, 20(7): 3005-3012.

[15] Zhang Q W, Dong Y, Li F, et al. Quantifying detachment rate of eroding rill or ephemeral gully for WEPP with flume experiments[J]Journal of Hydrology, 2014, 519: 2012-2019.

[16] 肖培青，鄭粉莉，姚文藝. 坡溝系統(tǒng)侵蝕產(chǎn)沙及其耦合關(guān)系研究[J]泥沙研究，2007(2)：30-35.

Xiao Peiqing, Zheng Fenli, Yao Wenyi. Study on the sediment yield and coupling mechanism of slope-gully system[J]. Journal of Sediment Research, 2007(2): 30-35. (in Chinese with English abstract)

[17] 黃鈺涵，楊夢(mèng)格，雷廷武，等. 飽和狀態(tài)下黃綿土坡面細(xì)溝侵蝕可蝕性和臨界剪切應(yīng)力特征[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2022，38(17)：83-90.

Huang Yuhan, Yang Mengge, Lei Tingwu, et al. Characteristics of soil erodibility parameter and critical shear stress on saturated loess soil slope[J]Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(17): 83-90. (in Chinese with English abstract)

[18] Léonard J, Richard G. Estimation of runoff critical shear stress for soil erosion from soil shear strength[J]Catena, 2004, 57(3): 233-249.

[19] Govers G, Giménez R, Oost K V. Rill erosion: exploring the relationship between experiments, modelling and field observations[J]Earth-Science Reviews, 2007, 84(3): 87-102.

[20] Xu X M, Zheng F L,Wilson G V. Flow hydraulics in an ephemeral gully system under different slope gradients, rainfall intensities and inflow conditions[J]Catena, 2021, 203: 105359.

[21] 郭明明，王文龍，李建明，等. 黃土區(qū)坡耕地耕作對(duì)淺溝徑流產(chǎn)沙及其形態(tài)發(fā)育特征的影響[J]農(nóng)業(yè)工程學(xué)報(bào)，2015，31(15)：114-123.

Guo Mingming, Wang Wenlong, Li Jianming, et al. Effect of tillage on runoff and sediment yields and morphology development characteristic of ephemeral gully in loessial region[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(15): 114-123. (in Chinese with English abstract)

[22] Istanbulluoglu E, Tarboton D G, Pack R T, et al. A sediment transport model for incision of gullies on steep topography[J]Water Resources Research, 2003, 39(4): 1103.

[23] Vanwalleghem T, Giráldez J V, Jiménezhornero F J, et al. Evaluating a general sediment transport model for linear incisions under field conditions[J]Earth Surface Processes Landforms, 2010, 34(14): 1852-1857.

[24] 王清源，張樂(lè)濤，魏儀媛，等. 基于能量的坡面侵蝕性徑流及其水沙傳遞關(guān)系[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2022，38(1)：151-157.

Wang Qingyuan, Zhang Letao, Wei Yiyuan, et al. Relationship between erosive runoff and flow-sediment on loess slope: An erosive energy perspective [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(1): 151-157. (in Chinese with English abstract)

[25] 肖海，劉剛，劉普靈. 集中流作用下黃土坡面剝蝕率對(duì)侵蝕動(dòng)力學(xué)參數(shù)的響應(yīng)[J]農(nóng)業(yè)工程學(xué)報(bào)，2016，32(17)：106-111.

Xiao Hai, Liu Gang, Liu Puling. Response of detachment rate of loess slope to hydrodynamic characteristics under concentrate flow condition[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(17): 106-111. (in Chinese with English abstract)

[26] 張漢雄，王萬(wàn)忠. 黃土高原的暴雨特性及分布規(guī)律[J]水土保持通報(bào)，1982(1)：35-43.

Zhang Hanxiong, Wang Wanzhong. Characteristics and distribution of rainstorms on the Loess Plateau[J]. Bulletin of Soil and Water Conservation, 1982(1): 35-43. (in Chinese with English abstract)

[27] Zhang L T, Gao Z L, Yang S W, 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.

[28] Bagnold R A. An approach to the sediment transport problem from general physics[C]//The physics of sediment transport by wind and water. Washington, DC: US Gov. Print. Off, 1966, 422, I.

[29] Wells R R, Bennett S J, Alonso C V. Effect of soil texture, tailwater height, and pore-water pressure on the morphodynamics of migrating headcuts in upland concentrated flows[J]Earth Surface Processes & Landforms, 2010, 34(14): 1867-1877.

[30] 鄭粉莉. 坡面細(xì)溝侵蝕過(guò)程量化研究[M]北京：科學(xué)出版社，2020.

[31] Stefanovic J R, Bryan R B. Experimental study of rill bank collapse[J]Earth Surface Processes Landforms, 2007, 32: 180-196.

[32] Bryan R B, Poesen J. Laboratory experiments on the influence of slope length on runoff, percolation and rill development[J]Earth Surface Processes & Landforms, 2010, 14(3): 211-231.

[33] 李鎮(zhèn)，張巖，楊松，等. QuickBird影像目視解譯法提取切溝形態(tài)參數(shù)的精度分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2014，30(20)：179-186.

Li Zhen, Zhang Yan, Yang Song, et al. Error assessment of extracting morphological parameters of bank gullies by manual visual interpretation based on QuickBird imagery[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(20): 179-186. (in Chinese with English abstract)

[34] 牛耀彬，吳旭，高照良，等. 降雨和上方來(lái)水條件下工程堆積體坡面土壤侵蝕特征[J]農(nóng)業(yè)工程學(xué)報(bào)，2020，36(8)：69-77.

Niu Yaobin, Wu Xu, Gao Zhaoliang, et al. Characteristics of soil erosion on engineering accumulation slope under the rainfall and inflow conditions[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(8): 69-77. (in Chinese with English abstract)

[35] Torri D, Borselli L. Equation for high-rate gully erosion[J]Catena, 2003, 50(2/3/4): 449-467.

[36] Slattery M C, Bryan R B. Hydraulic conditions for rill incision under simulated rainfall: A laboratory experiment[J]Earth Surface Processes & Landforms, 1992, 17(2): 127-146.

[37] 王秋霞，丁樹(shù)文，鄧羽松，等. 花崗巖崩崗區(qū)不同土層的侵蝕水動(dòng)力學(xué)特征[J]土壤學(xué)報(bào)，2017，54(3)：570-580.

Wang Qiuxia, Ding Shuwen, Deng Yusong, et al. Hydrodynamic characteristics of erosion in different soil layers in granite collapse region[J]. Acta Pedologica Sinica, 2017, 54(3): 570-580. (in Chinese with English abstract)

[38] 丁文斌，史東梅，何文健，等. 放水沖刷條件下工程堆積體邊坡徑流侵蝕水動(dòng)力學(xué)特性[J]農(nóng)業(yè)工程學(xué)報(bào)，2016，32(18)：153-161.

Ding Wenbin, Shi Dongmei, He Wenjian, et al. Hydrodynamic characteristics of engineering accumulation erosion under side slope runoff erosion process in field scouring experiment[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(18): 153-161. (in Chinese with English abstract)

[39] Xu X M, Zheng F L, Qin C, et al. Impact of cornstalk buffer strip on hillslope soil erosion and its hydrodynamic understanding[J]Catena, 2017, 149: 417-425.

[40] Yager E M, Kirchner J W, Dietrich W E. Calculating bed load transport in steep boulder bed channels[J]Water Resources Research, 2007, 43(7): 256-260.

[41] Wang J G, Feng S Y, Ni S M, et al. Soil detachment by overland flow on hillslopes with permanent gullies in the Granite area of southeast China[J]Catena, 2019, 183: 104235.

Characteristics of gully topography and sediment on the platform-steep slope system of spoil dump

Bai Yun1,2,3,4, Kang Hongliang5, Wang Wenlong1,2,6※, Su Huan7, Li Jianming1,2,4,8, Ma Chunyan3

(1.712100,; 2.,,712100,; 3.,,719000,; 4.,100049,; 5.,,710054,; 6.,,712100,; 7.210000; 8.,,430010,)

The spoil dumps are the main sediment sources in the open-pit mining area. The erosion control is of great significance to the high-quality development of the energy zone. This article aims to analyze the gully development and sediment process of spoil dumps. The plot construction of platform-steep slope system and field scouring experiments were also carried out, where the flow rate was ranged from 60 to 84 L/min, and the duration was 45 mins each run. The results showed that: 1) The gully was developed by the combination of different erosion progresses, including the headcut migration, bed incision, and lateral erosion on the platform-steep slope system. Notably, there were the stage differences in the topographic characteristics of gully, due to the variation of dominant process. A three-stage development was experienced for the gully on the platform, including the headcut formation, migration-expansion, and stable stage. In the gully on the steep slope, four developmental processes were divided into the incision of sandy loam layer, expansion of sandy loam layer, incision of clay loam layer, and deceleration of erosion. No significant differences were found between the various flow rates of gully development speed on the platform. Whereas, the transformed time increased with the flow rate between adjacent gully development stages on the steep slope. 2) In terms of sediment process, the migration - expansion stage was the main erosion period on the platform. The incision of sandy and clay loam layer were the main erosion stages on the steep slope, particularly with the accumulated sediment yields accounting for 29.72%-53.36% and 19.06%-48.88%, respectively. Spatially, the runoff shear force and stream power increased by 7.11-120.86 times, and 7.59-239.59 times after the runoff flowed from the platform into the steep slope. Correspondingly, the steep slope was the main sediment source of the platform and steep slope system. The cumulative sediment yield was accounted for 88.15%-90.16% of the total amount of platform and steep slope system. Hence, the separation and control of runoff on the platform was a vital way to control the gully erosion in the platform-steep slope system of spoil dumps. 3) In terms of the flow hydraulics, the platform velocity decreased first and then increased with the gully development. Conversely, the slope velocity decreased gradually. Meanwhile, the runoff shear stress and the runoff stream power increased gradually with the gully development on the platform, but increased first, finally decreased on the steep slope. 4) Comparatively, the erosion rate presented a more sensitive response to the stream power than the runoff shear force, where the response regularity was different among the gully development stages. The erosion rate was linearly responded to the stream power in the migration-expansion stage and stable stage of platform. On the steep slope, the responses were still linear in the expansion of sandy loam layer and deceleration of erosion, but the exponential in the incision of sandy and clay loam layer. Therefore, it is necessary to consider the difference of responses in the various gully development stages, in order to improve the applicability of sediment transport model for the gully erosion. As such, the results revealed the staged development of gully and sediment process on the platform-steep slope system of spoil dumps. The finding can also offer the strong reference for the gully process, further to implement the soil and water conservation measures in the spoil dump.

soils; scouring test; gully erosion; stream power; erosion development; spoil dump

10.11975/j.issn.1002-6819.2022.23.009

S157.1

A

1002-6819(2022)-23-0081-10

白蕓，康宏亮，王文龍，等. 排土場(chǎng)平臺(tái)-邊坡系統(tǒng)溝蝕形態(tài)演變與產(chǎn)沙特征[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2022，38(23)：81-90.doi：10.11975/j.issn.1002-6819.2022.23.009 http://www.tcsae.org

Bai Yun, Kang Hongliang, Wang Wenlong, et al. Characteristics of gully topography and sediment on the platform-steep slope system of spoil dump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(23): 81-90. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2022.23.009 http://www.tcsae.org

2022-07-25

2022-11-22

國(guó)家自然科學(xué)基金資助項(xiàng)目（41761062，41867015，41701316）

白蕓，博士生，講師，研究方向?yàn)榈V區(qū)土壤侵蝕與水土保持。Email：baiyun21@mails.ucas.ac.cn

王文龍，博士，研究員，研究方向?yàn)橥寥狼治g與水土保持。Email：wlwang@nwsuaf.edu.cn