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

張孝中, 王文龍, 李建明

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

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

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

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

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

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

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

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

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

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

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

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

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

1.1　研究區(qū)概況

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

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

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

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

表1　棄土棄渣體顆粒機(jī)械組成

2結(jié)果與分析

2.1　棄土棄渣體平均土壤剝蝕率隨雨強(qiáng)變化

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

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

表2　棄土棄渣體平均土壤剝蝕率隨雨強(qiáng)的變化

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

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

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

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

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

2.2.2弗勞德數(shù)Fr流體力學(xué)中弗勞德數(shù)表征流體慣性力和重力相對(duì)大小的一個(gè)無(wú)量綱參數(shù)，記為Fr，它反映了過(guò)水?dāng)嗝嫔纤骶哂械膭?dòng)能和勢(shì)能的對(duì)比關(guān)系。按照清水明渠水流的標(biāo)準(zhǔn)，當(dāng)Fr<1時(shí)，水流為緩流；當(dāng)Fr>1時(shí)，水流為急流，其計(jì)算公式為[9]：

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

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

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

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

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

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

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

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

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

τ=ρgRJ=rRJ

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

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

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

τ=rqs=rhVs=τV

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

2.3.1平均土壤剝蝕率與水流剪切力關(guān)系圖4反映了3種不同下墊面次降雨平均土壤剝蝕率Dr與平均水流剪切力τ之間關(guān)系。從圖4可以看出，不同下墊面之間的平均土壤剝蝕率Dr隨平均水流剪切力 的增大整體呈遞增趨勢(shì)。計(jì)算可知，各下墊面的Dr和τ均表現(xiàn)為：沙多石少棄渣體>棄土體>沙少石多棄渣體。可以通過(guò)計(jì)算水流剪切力τ進(jì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隨著水流功率的增大呈顯著地遞增趨勢(shì)，且變化范圍較沙多石少棄渣體都較小。主要是由于沙多石少棄渣體中含有一定含量的礫石，礫石的存在破壞了原有土壤結(jié)構(gòu)，導(dǎo)致土壤的抗蝕能力降低，在徑流的沖刷作用下產(chǎn)生較大的土壤侵蝕量，同時(shí)由于沙多石少棄渣體坡面徑流平均流速也較棄土體和沙少石多棄渣體大，進(jìn)而使得沙多石少棄渣體的平均水流功率P最大。計(jì)算可知，3種不同下墊面的次降雨平均水流功率P的大小順序?yàn)椋荷扯嗍贄壴w>棄土體>沙少石多棄渣體。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隨雨強(qiáng)增大顯著遞增，3種下墊面土壤剝蝕率Dr的大小順序?yàn)椋荷扯嗍贄壴w>棄土體>沙少石多棄渣體，且各下墊面Dr值與雨強(qiáng)I之間呈顯著線性關(guān)系。因此，在進(jìn)行礦區(qū)水土流失治理過(guò)程中，尤其要注意對(duì)含有一定礫石含量的棄渣體采用水土保持措施來(lái)減輕礦區(qū)嚴(yán)重的水土流失現(xiàn)狀。

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

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

[參考文獻(xiàn)]

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

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

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

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

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

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

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

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

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

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

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

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

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

[14]沙際德,白清俊.黏性土坡面細(xì)溝流的水力特性試驗(yàn)研究[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ì)溝淺溝切溝流動(dòng)力機(jī)制研究[J].水科學(xué)進(jìn)展,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

文獻(xiàn)標(biāo)識(shí)碼：A

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

中圖分類號(hào)：S157.1

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