張孝中, 王文龍, 李建明

(1.陜西省水土保持勘測規(guī)劃研究所, 陜西 西安 710004; 2.西北農(nóng)林科技大學(xué) 水土保持研究所 黃土高原土壤

侵蝕與旱地農(nóng)業(yè)國家重點實驗室, 陜西 楊凌 712100; 3.中國科學(xué)院 水利部 水土保持研究所, 陜西 楊凌 712100)

神府東勝煤田棄土棄渣體徑流水動力學(xué)特性研究

張孝中1, 王文龍2,3, 李建明2

(1.陜西省水土保持勘測規(guī)劃研究所, 陜西 西安 710004; 2.西北農(nóng)林科技大學(xué) 水土保持研究所 黃土高原土壤

侵蝕與旱地農(nóng)業(yè)國家重點實驗室, 陜西 楊凌 712100; 3.中國科學(xué)院 水利部 水土保持研究所, 陜西 楊凌 712100)

資助項目：國家自然科學(xué)基金項目“煤礦開發(fā)建設(shè)中人為侵蝕機理及新增水土流失預(yù)報模型：以神府煤田為例”(40771127); 水利部公益性行業(yè)專項(201201047； 201201048)

第一作者：張孝中(1961—),男(漢族),陜西省澄城縣人,碩士，高級工程師，主要從事水土流失動態(tài)規(guī)律及生態(tài)恢復(fù)方面的研究。E-mail:zxz211@126.com。

摘要：[目的] 對神府煤田開發(fā)建設(shè)過程中造成的棄土棄渣體的徑流水動力學(xué)特性進行研究，并分析土壤剝蝕率與各水動力學(xué)參數(shù)之間的關(guān)系。 [方法] 采用野外模擬降雨試驗方法。[結(jié)果] (1) 棄土棄渣體的土壤剝蝕率Dspan隨雨強增大顯著線性遞增； (2) 雷諾數(shù)Rspan總體上隨著雨強增大而增大，隨降雨歷時先增大后保持穩(wěn)定；佛羅德數(shù)Fr受雨強影響不顯著，在產(chǎn)流0~12 min內(nèi)波動較大，隨后在某一常數(shù)附近上下波動；阻力系數(shù)f隨降雨歷時呈遞增趨勢，受雨強影響不顯著，次降雨平均阻力系數(shù)f大小表現(xiàn)為：沙少石多棄渣體>棄土體>沙多石少棄渣體； (3) 各下墊面土壤剝蝕率Dspan與水流剪切力 、水流功率P之間均呈顯著的線性函數(shù)關(guān)系。 [結(jié)論] 研究區(qū)不同類型下墊面次降雨徑流水動力學(xué)參數(shù)變化各具特點，土壤剝蝕率與各水動力學(xué)參數(shù)之間存在顯著相關(guān)關(guān)系。

關(guān)鍵詞：棄土棄渣體; 模擬降雨; 水動力學(xué)參數(shù); 土壤剝蝕率

神府東勝煤田的開發(fā)，確保了國家優(yōu)質(zhì)動力煤的穩(wěn)定供應(yīng)，促進國民經(jīng)濟的快速發(fā)展，但煤炭開采過程中擾動地表、修路、修建臨時措施、棄土棄渣隨意堆放等都會對礦區(qū)的環(huán)境造成嚴(yán)重破壞，產(chǎn)生一系列的環(huán)境問題[1-3]。為了更加準(zhǔn)確地預(yù)測礦區(qū)新增水土流失，許多學(xué)者針對神府東勝煤田的土壤侵蝕現(xiàn)狀展開了大量研究，張勝利[4]等在調(diào)查研究及室內(nèi)外試驗基礎(chǔ)上，分析了煤炭開采對侵蝕和產(chǎn)沙的影響，認(rèn)為煤田開發(fā)會加重水土流失；李強、羅婷等[5-6]針對神府煤田不同下墊面徑流產(chǎn)沙特征進行了詳細(xì)研究。調(diào)查發(fā)現(xiàn)，礦區(qū)中堆放的大量棄土棄渣體是礦區(qū)新增水土流失的主要策源地之一。

對礦區(qū)棄土棄渣的研究主要通過野外調(diào)查及野外放水沖刷試驗[7-8]。為了更加明確土壤侵蝕的內(nèi)在機理，學(xué)者們開展了大量的水動力學(xué)試驗研究，魏霞等[9]采用室內(nèi)放水沖刷試驗研究水動力學(xué)特性，結(jié)果表明雷諾數(shù)隨沖刷歷時增大而增大，流速沿程有先增大后減小趨勢；張科利等[10]運用水力學(xué)理論及試驗結(jié)果，提出了水動力學(xué)參數(shù)計算公式；王貞等[11]對比分析了神府煤田擾動地面及原生地面的水動力學(xué)參數(shù)，發(fā)現(xiàn)擾動地面的雷諾數(shù)和阻力系數(shù)均大于原地面，佛羅德數(shù)差異不顯著。目前，針對礦區(qū)不同下墊面侵蝕產(chǎn)沙特性的研究較多，對礦區(qū)棄土棄渣體的徑流產(chǎn)沙規(guī)律也進行了大量的試驗研究，但針對棄土棄渣體的水動力學(xué)特性試驗研究鮮見有報道。本文采用野外人工模擬降雨試驗，研究礦區(qū)棄土棄渣體徑流水動力學(xué)特性，以期為明確礦區(qū)棄土棄渣體土壤侵蝕內(nèi)在機理研究提供一定的科學(xué)理論依據(jù)。

1研究區(qū)概況與研究方法

1.1　研究區(qū)概況

神府東勝煤田地處晉陜蒙三省交界處(北緯37°20′—40°16′，東經(jīng)108°36′—110°36′)，位于黃土丘陵與鄂爾多斯臺地的過渡地帶，包括黃土丘陵、風(fēng)沙丘陵與砒沙巖地貌類型，是典型的生態(tài)脆弱區(qū)。總面積約3 481 km2，煤炭地質(zhì)儲量3.54×1010t，探明儲量占全國煤炭保有儲量的1/4。該區(qū)為典型的干旱半干旱大陸季風(fēng)性氣候，年平均氣溫6.1～9.1 ℃，年平均降水量不足400 mm，降水多集中在多集中于7—9月，占全年降雨的65%～70%，分布趨勢由東南向西北逐漸減少。年蒸發(fā)量為1 636～2 535 mm，蒸發(fā)強烈。黃土是本區(qū)的最主要地面物質(zhì)，分布面積約占63.63%以上，一般質(zhì)地較粗，抗蝕性、結(jié)構(gòu)性極差，是易侵蝕的主要物質(zhì)之一。

1.2　試驗布設(shè)與方法

試驗區(qū)位于神木縣西溝鄉(xiāng)六道溝煤礦。選取棄土體、沙多石少棄渣體、沙少石多棄渣體3種不同類型下墊面，土壤機械組成見表1。試驗小區(qū)長3 m，寬1 m，四周采用高30cm，厚2mm的鋼板控制邊界條件，鋼板露出地面10cm，以防止徑流滲漏。為降低邊界效應(yīng)，靠近板壁處土體應(yīng)盡量壓實。下墊面坡度根據(jù)研究區(qū)棄土棄渣體典型坡度選擇為35°。降雨強度按照神府東勝地區(qū)多年降雨頻率氣象資料分析，根據(jù)20年一遇頻率標(biāo)準(zhǔn)，采用1.0，1.5，2.0，2.5，3.0 mm/min 5個雨強級別。

采用水泵將水抽至小區(qū)上部放置的2 m3的儲水桶內(nèi)，通過60 m揚程的潛水泵將水抽出進行降雨，設(shè)置雙閥門和壓力表調(diào)節(jié)降雨強度。試驗露天進行，采用鋼管搭建一個高3 m的降雨架，為防止風(fēng)對降雨分布的影響，在小區(qū)四周用彩條布圍一圈，小區(qū)正上方布設(shè)下噴式降雨器。為保證降雨強度達到試驗要求，在降雨試驗開始前，需要率定降雨強度，保證率定雨強與設(shè)計雨強誤差不超過5%，雨強率定完畢，迅速揭開遮雨布并開始計時。自降雨開始至坡面有較明顯徑流形成并流出小區(qū)出口時記錄產(chǎn)流起始時間，產(chǎn)流開始后的3 min內(nèi)每1 min接一個水樣，3 min后每隔3 min接一個水樣。設(shè)置產(chǎn)流開始后降雨總時間為45 min。同時在各小區(qū)中間部位分別設(shè)置兩個1 m長的觀測斷面，以便測定流速、水深及流寬：徑流流速測定采用高錳酸鉀示蹤法，由于測得的流速是徑流表面的最大流速，雷廷武等[12]將染色劑測量得到的流速乘以系數(shù)0.75來得到較為理想的水流平均流速，本試驗亦以此法得到水流平均速度；用薄鋼尺在測量斷面處測徑流寬度；在觀測斷面內(nèi)選取三點分別用薄鋼尺測定徑流深，以平均值作為整個小區(qū)的徑流深。試驗結(jié)束后，用量筒測定各個徑流樣體積，采用烘干法測定泥沙量。采用SPSS18.0進行數(shù)據(jù)處理分析。

表1　棄土棄渣體顆粒機械組成

2結(jié)果與分析

2.1　棄土棄渣體平均土壤剝蝕率隨雨強變化

坡面水蝕過程的土壤剝蝕率為單位時間單位面積上土體在侵蝕動力的作用下被剝蝕掉的土壤顆粒的量[13]。次降雨平均土壤剝蝕率是指次降雨過程中各時段土壤剝蝕率的平均值，以Dr表示，可以反映出徑流對下墊面表層的沖刷、剝蝕能力。表2為各下墊面在不同雨強下次降雨平均土壤剝蝕率。分析表2可知，棄土棄渣體的Dr隨雨強增大顯著遞增(p<0.01)，隨降雨強度的增大，一方面，單位時間單位面積內(nèi)下墊面的承雨量增加，徑流強度增大，沖刷能力增強；另一方面，雨滴直徑增大，對地表的沖擊能力增強，濺蝕作用強烈。3種不同下墊面Dr值的大小順序為：沙多石少棄渣體>棄土體>沙少石多棄渣體，這主要是由于3種下墊面的物質(zhì)組成不同，棄土體中礫石含量較少，其土壤結(jié)構(gòu)類似于原生地面，具有一定的抗侵蝕能力；沙多石少棄渣體中含有一定量(約45%)礫石，而下墊面表層礫石存在使得土粒之間的黏結(jié)性減弱，在雨水沖刷過后很快隨徑流流失，礫石存在使得下墊面的抗蝕能力降低，且表層易侵蝕性顆粒多，各雨強條件下土壤剝蝕率最大；沙少石多棄渣體中礫石含量較大(約70%)，在下墊面的表層形成一定的固相骨架，而棄土則分布在礫石空隙之間，由于可供侵蝕的物質(zhì)源較少，導(dǎo)致總體的Dr值最小。與棄土體相比，沙多石少棄渣體的Dr為棄土體的1.85～2.76倍，而沙少石多棄渣體僅為棄土體的0.57%～21.30%?？梢姡V區(qū)嚴(yán)重的土壤侵蝕可能是由于含礫石棄渣體的隨意堆放造成的，在礦區(qū)水土流失防治過程中尤其要注重土石混合堆積體的治理。

不同的降雨強度代表著不同降雨侵蝕營力，因此，分析次降雨土壤剝蝕率與雨強的關(guān)系對于礦區(qū)進行水土流失預(yù)測具有重要的現(xiàn)實指導(dǎo)意義。分析可知，棄土體、沙多石少棄渣體、沙少石多棄渣體的Dr與雨強I之間均呈顯著地線性關(guān)系，其相關(guān)系數(shù)分別為0.929(p=0.05)，0.963(p=0.01)和0.980(p=0.01)。

表2　棄土棄渣體平均土壤剝蝕率隨雨強的變化

2.2　棄土棄渣體水動力學(xué)特性

2.2.1雷諾數(shù)Re雷諾數(shù)Re(Reynolds number)是一種可用來表征流體流動情況的無量綱數(shù)，從力學(xué)上講是慣性力與黏性力之比，可以用來判定水流的流態(tài)。當(dāng)Re<500，水流為層流，當(dāng)Re>2 000，水流為紊流，考慮到含沙量對水流運動黏滯性的影響，雷諾數(shù)計算式為[14]：

式中：V——平均流速(m/s)；R——水力半徑(m)，R=A/L；A——過水?dāng)嗝?；L——濕周；v——含沙水流的運動黏滯系數(shù)(m2/s)，其中，

式中：t——水溫(℃)。

圖1分別描繪了不同下墊面在不同雨強下的雷諾數(shù)隨降雨時間的變化過程。從圖1可以看出，棄土棄渣體的Re總體上隨著雨強增大呈遞增趨勢，隨降雨歷時呈現(xiàn)先增大后穩(wěn)定的趨勢。各下墊面Re在雨強為1.0 mm/min時維持在某一個常數(shù)上下小范圍波動，可認(rèn)為此時的水流大部分處于層流狀態(tài)。而當(dāng)雨強達到3.0 mm/min時，Re隨降雨歷時在大范圍內(nèi)起伏波動，處于過渡流與紊流不斷交換的狀態(tài)。而在1.5，2.0，2.5 mm/min雨強下，Re隨產(chǎn)流時間的持續(xù)總體上呈先波動后保持相對穩(wěn)定的狀態(tài)。由于在產(chǎn)流初始，坡面上有較多的可蝕性顆粒，水流的剝蝕能力與下墊面的抗蝕能力處于非平衡狀態(tài)，造成水流不穩(wěn)定，進而導(dǎo)致水流流速的不斷變化，使得Re不斷波動；而當(dāng)產(chǎn)流開始后6 min時下墊面達到相對穩(wěn)定，此時的水流剝蝕能力與下墊面的抗蝕能力相當(dāng)，水流處于較平穩(wěn)狀態(tài)，Re的變化也較小。經(jīng)計算，3種不同下墊面的平均雷諾數(shù)的大小順序為：沙多石少棄渣體>沙少石多棄渣體>棄土體。

2.2.2弗勞德數(shù)Fr流體力學(xué)中弗勞德數(shù)表征流體慣性力和重力相對大小的一個無量綱參數(shù)，記為Fr，它反映了過水?dāng)嗝嫔纤骶哂械膭幽芎蛣菽艿膶Ρ汝P(guān)系。按照清水明渠水流的標(biāo)準(zhǔn)，當(dāng)Fr<1時，水流為緩流；當(dāng)Fr>1時，水流為急流，其計算公式為[9]：

式中：u——水流的平均流速(m/s)； g——重力加速度(m/s)； l——斷面平均徑流深(m)。

圖2分別為棄土棄渣體在不同降雨強度條件下的Fr隨降雨歷時的變化。由圖2可知，3種不同下墊面的Fr隨雨強的變化不顯著，但隨產(chǎn)流歷時的變化總體可以分為兩個階段：產(chǎn)流0~12min內(nèi)Fr隨降雨歷時上下波動；產(chǎn)流12min后Fr在某一常數(shù)值附近上下小幅度波動，其中棄土體和沙少石多棄渣體的Fr在0.50左右，而沙多石少棄渣體的Fr在0.60左右。造成3種下墊面Fr的這種變化的主要原因是下墊面組成的不一致性，即土石比的不一致性。在產(chǎn)流初期，由于徑流對表土進行不斷地沖刷、剝離、搬運、沉積，導(dǎo)致了水流的不穩(wěn)定，在后期由于下墊面的抗蝕能力與水流的侵蝕能力達到平衡狀態(tài)，水流進而趨于穩(wěn)定。經(jīng)分析，在整個降雨過程中，3種下墊面85%以上Fr<1.0，即水流處于緩流狀態(tài)。計算可知，3種不同下墊面的Fr的大小順序為：沙多石少棄渣體>棄土體>沙少石多棄渣體。由此可知，當(dāng)下墊面中的礫石含量達到一定值時，有利于減緩坡面水流的運動速度，進而降低水土流失量。

圖1　不同下墊面(棄土棄渣體)雷諾數(shù)Re隨時間的變化

圖2　不同下墊面(棄土棄渣體)弗勞德數(shù)Fr隨時間的變化

2.2.3阻力系數(shù)fDarcy—Weisbach阻力系數(shù)是徑流向下運動過程中受到的來自水土界面的阻滯水流運動力的總稱。其表達式為[15]:

式中：V——斷面平均流速； J——水力能坡，計算式為：

圖3分別為棄土棄渣體在不同降雨強度下的阻力系數(shù)f隨降雨歷時的變化過程。由圖3可以看出，3種下墊面的阻力系數(shù)隨著降雨歷時呈遞增趨勢，而隨雨強的變化不顯著。從計算公式中可以看出，f主要受到水流流速的影響。分析圖3可知，沙多石少棄渣體下墊面的f隨降雨歷時的波動情況較棄土體和沙少石多棄渣體大，主要是由于棄土體中砂礫含量少，土粒之間的黏結(jié)性較強，流速在整個降雨過程中的變化趨勢較為平緩；而沙少石多棄渣體中由于含有大量的礫石，當(dāng)?shù)[石之間的可蝕性顆粒被沖走后，形成一個穩(wěn)定的固相骨架，有利于坡面的穩(wěn)定，而礫石存在改變了水流的流速，但水流流動速度較均勻，導(dǎo)致了流速在降雨過程中的變化不大，進而使得f在侵蝕過程中變化較平緩。經(jīng)計算，平均阻力系數(shù)f大小順序為：沙少石多棄渣體>棄土體>沙多石少棄渣體。

2.3　棄土棄渣體土壤剝蝕率與水流剪切力、水流功率關(guān)系

為了更加明確棄土棄渣體土壤侵蝕的力學(xué)性質(zhì)，對棄土棄渣體的平均土壤剝蝕率與水流剪切力及水流功率關(guān)系進行相關(guān)分析。其中，水流剪切力τ為沿著坡面梯度方向運動的水流在其運動方向上產(chǎn)生的一個作用力[13]，計算公式為[16]:

τ=ρgRJ=rRJ

式中：r——水的容重，r=ρg，水流密度ρ取平均含沙量時的密度。

圖3　不同下墊面(棄土棄渣體)阻力系數(shù)f隨時間的變化

水流功率P表征了一定高度的水體順坡流動時具有的勢能。Bagnold[17]提出了水流功率的概念，即作用于單位面積的水流所消耗的功率，表達式為：

τ=rqs=rhVs=τV

式中：q——單寬流量(m3)；h——水深(m)。

2.3.1平均土壤剝蝕率與水流剪切力關(guān)系圖4反映了3種不同下墊面次降雨平均土壤剝蝕率Dr與平均水流剪切力τ之間關(guān)系。從圖4可以看出，不同下墊面之間的平均土壤剝蝕率Dr隨平均水流剪切力 的增大整體呈遞增趨勢。計算可知，各下墊面的Dr和τ均表現(xiàn)為：沙多石少棄渣體>棄土體>沙少石多棄渣體?？梢酝ㄟ^計算水流剪切力τ進而求得次降雨的土壤剝蝕率Dr。經(jīng)分析，棄土體、沙多石少棄渣體和沙少石多棄渣體的次降雨的平均水流剪切力τ與平均土壤剝蝕率Dr之間呈顯著線性關(guān)系，相關(guān)系數(shù)分別為0.995(p=0.01)，0.914(p=0.05)和0.955(p=0.05)。

圖4　平均土壤剝蝕率Dr與水流剪切力τ的關(guān)系

2.3.2平均土壤剝蝕率與水流功率關(guān)系圖5反映了3種不同下墊面的平均土壤剝蝕率Dr隨水流功率P變化規(guī)律。從圖5可以看出，棄土體和沙少石多棄渣體的平均土壤剝蝕率Dr隨著水流功率的增大呈顯著地遞增趨勢，且變化范圍較沙多石少棄渣體都較小。主要是由于沙多石少棄渣體中含有一定含量的礫石，礫石的存在破壞了原有土壤結(jié)構(gòu)，導(dǎo)致土壤的抗蝕能力降低，在徑流的沖刷作用下產(chǎn)生較大的土壤侵蝕量，同時由于沙多石少棄渣體坡面徑流平均流速也較棄土體和沙少石多棄渣體大，進而使得沙多石少棄渣體的平均水流功率P最大。計算可知，3種不同下墊面的次降雨平均水流功率P的大小順序為：沙多石少棄渣體>棄土體>沙少石多棄渣體。3種下墊面次降雨平均土壤剝蝕率Dr均可用平均水流功率P描述。經(jīng)分析，棄土體、沙多石少棄渣體和沙少石多棄渣體的次降雨平均土壤剝蝕率Dr與水流功率P之間呈顯著地線性關(guān)系，相關(guān)系數(shù)分別為0.992(p=0.01)，0.912(p=0.05)和0.973(p=0.01)。

圖5　平均土壤剝蝕率Dr與水流功率P的關(guān)系

3結(jié) 論

(1) 棄土棄渣體的土壤剝蝕率Dr隨雨強增大顯著遞增，3種下墊面土壤剝蝕率Dr的大小順序為：沙多石少棄渣體>棄土體>沙少石多棄渣體，且各下墊面Dr值與雨強I之間呈顯著線性關(guān)系。因此，在進行礦區(qū)水土流失治理過程中，尤其要注意對含有一定礫石含量的棄渣體采用水土保持措施來減輕礦區(qū)嚴(yán)重的水土流失現(xiàn)狀。

(2) 分析棄土棄渣體的3個水動力學(xué)特性參數(shù)(Re，F(xiàn)r，f)可知，棄土棄渣體的Re總體上隨著雨強增大呈遞增趨勢，隨降雨歷時呈現(xiàn)先增大后穩(wěn)定的趨勢，3種不同下墊面的平均雷諾數(shù)的大小順序為：沙多石少棄渣體>沙少石多棄渣體>棄土體；3種不同下墊面的Fr隨雨強的變化不顯著，但隨產(chǎn)流開始后時間的變化總體可以分為兩個階段：產(chǎn)流0~12 min內(nèi)Fr隨降雨歷時上下波動；產(chǎn)流12 min后Fr在某一常數(shù)值附近上下小幅度波動；阻力系數(shù)f隨著降雨歷時呈遞增趨勢，而隨雨強的變化不顯著，平均阻力系數(shù)f大小順序為：沙少石多棄渣體>棄土體>沙多石少棄渣體。

(3) 次降雨過程中，3種不同下墊面的平均土壤剝蝕率Dr隨著平均水流剪切力 和水流功率P均呈線性遞增關(guān)系。

[參考文獻]

[1]王文龍,李占斌,張平倉.神府東勝煤田開發(fā)中誘發(fā)的環(huán)境災(zāi)害問題研究[J].生態(tài)學(xué)雜志,2004,23(1):34-38.

[2]張發(fā)旺,趙紅梅,宋亞新,等.神府東勝礦區(qū)采煤塌陷對水環(huán)境影響效應(yīng)研究[J].地球?qū)W報,2007,28(6):521-527.

[3]武強,安永會,劉文崗,等.神府東勝礦區(qū)水土環(huán)境問題及其調(diào)控技術(shù)[J].煤田地質(zhì)與勘探,2005,33(3):54-57.

[4]張勝利,時明立,張利銘.神府東勝煤田開發(fā)對侵蝕產(chǎn)沙的影響[J].水土保持學(xué)報,1992,6(2):68-75.

[5]李強,李占斌,魯克新.神府東勝礦區(qū)不同下墊面產(chǎn)流產(chǎn)沙試驗研究[J].水土保持研究,2005,15(3):1-3.

[6]羅婷,王文龍,李宏偉,等.神府東勝煤田不同下墊面侵蝕特征野外試驗[J].中國水土保持科學(xué),2012,10(6):52-57.

[7]王文龍,李占斌,李鵬.神府東勝煤田開發(fā)建設(shè)棄土棄渣沖刷試驗研究[J].水土保持學(xué)報,2004,18(5):68-71.

[8]王玲玲,劉蘭玉,左仲國.開發(fā)建設(shè)項目棄土棄渣水土流失調(diào)查分析[J].水力發(fā)電,2009,35(1):1-3.

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

[10]張科利,唐克麗.黃土坡面細(xì)溝侵蝕能力的水動力學(xué)試驗研究[J].土壤學(xué)報,2000,37(1):9-15.

[11]王貞,王文龍,金劍,等.神東煤田擾動地面與原地面產(chǎn)流產(chǎn)沙及水動力學(xué)參數(shù)對比[J].中國水土保持科學(xué),2010,8(6):69-74.

[12]雷廷武,張晴雯,閆麗娟.細(xì)溝侵蝕物理模型[M].北京:科學(xué)出版社,2009:117-129.

[13]王瑄,李占斌,尚佰曉,等.坡面土壤剝蝕率與水蝕因子關(guān)系室內(nèi)模擬試驗[J].農(nóng)業(yè)工程學(xué)報,2008,24(9):22-26.

[14]沙際德,白清俊.黏性土坡面細(xì)溝流的水力特性試驗研究[J].泥沙研究,2001(6):39-43.

[15]Foster G R, Huggins L F, Meyer L D. Laboratory study of rill hydraulics(Ⅱ): Shear stress relationships[J]. Transactions of the ASAE, 1984,27(3):797-804.

[16]王文龍,雷阿林,李占斌.土壤侵蝕鏈內(nèi)細(xì)溝淺溝切溝流動力機制研究[J].水科學(xué)進展,2003,14(4):471-475.

[17]Bagnold R. An approach to the Sediment Transport Problem from General Physics[M]. Washington: United States Government Printing Office, 1966.

Hydrodynamic Characteristics of Runoff over Abandoned Slags in Shenfu-Dongsheng Coalfield

ZHANG Xiaozhong1, WANG Wenlong2,3, LI Jianming2

(1.ShaanxiProvincialPlanningandResearchInstituteofSoilandWaterConservation,Xi’an,Shaanxi710004,China; 2.StateKeyLaboratoryofSoilErosionandDrylandFarmingontheLoessPlateau,InstituteofSoilandWaterConservation,NorthwestA&FUniversity,Yangling,Shaanxi712100,China; 3.InstituteofSoilandWaterConservation,CASandMWR,Yangling,Shaanxi712100,China)

Abstract：[Objective] Hydrodynamic characteristics of runoff over abandoned slags derived from the development and construction of the Shenfu—Dongsheng coal field were studied and the relationships between soil detachment rate and hydrodynamic parameters were analyzed. [Methods] Field simulated rainfall experiment was conducted. [Results] (1) Soil detachment rate(Dspan) for the abandoned slags significantly linearly increased with increasing rainfall intensity； (2) Reynolds number(Rspan) for the abandon slags generally increased as rainfall intensity increased, but its increase with rainfall duration showed an increasing trend at first and then maintained a stable trend. However, rainfall intensity did not affect Froude number(Fspan) significantly. The Fspanfluctuated violently 12 min after runoff generation and then fluctuated around a constant value. The drag coefficient(f) presented an increasing trend with rainfall duration, but the effect of rainfall intensity on the f was not significant. The abandoned slags studied, in order of decreasing average drag coefficient for a single rainfall, were the slag with less sand and more stone, abandoned soil body and slag with more sand and less stone； (3) The averaged soil detachment rate for each of the three underlying surfaces was a significantly linear function of shear force(τ) and stream power(P). [Conclusion] For the different types of underlying surfaces in the studied area, the variations of hydrodynamic parameters of runoff induced by single rainfall have their own characteristics. Soil detachment rate is significantly correlated with the hydrodynamic parameters.

Keywords:abandoned slag; simulated rainfall; hydrodynamic parameter; soil detachment rate

文獻標(biāo)識碼：A

文章編號：1000-288X(2015)02-0076-06

中圖分類號：S157.1

收稿日期：2014-03-19修回日期：2014-04-01