李 楠，吳超君，郝 喆，滕 達(dá)，王曉明

(1.遼寧大學(xué) 環(huán)境學(xué)院，遼寧 沈陽(yáng) 110036；2.遼寧有色勘察研究院有限責(zé)任公司，遼寧 沈陽(yáng) 110013)

0 引言

隨著我國(guó)經(jīng)濟(jì)的發(fā)展，每年全國(guó)礦山所排棄的尾礦越來(lái)越多，這些尾礦僅小部分被利用，剩下的大部分都存放在尾礦庫(kù)中，其中很多尾礦庫(kù)處于無(wú)人看管的廢棄狀況。結(jié)構(gòu)極其不穩(wěn)定的廢棄尾礦庫(kù)水土流失嚴(yán)重，致使礦區(qū)及其附近環(huán)境受到污染。明確廢棄尾礦壩坡面的產(chǎn)流產(chǎn)沙特征以及初始產(chǎn)流時(shí)間，可為提出合理水土保持措施提供理論依據(jù)。

國(guó)內(nèi)外諸多學(xué)者通過(guò)人工模擬降雨對(duì)產(chǎn)流產(chǎn)沙規(guī)律進(jìn)行了探究。Yuli Zhao[1]等通過(guò)模擬降雨對(duì)不同雨強(qiáng)下的植被產(chǎn)流產(chǎn)沙特征進(jìn)行了研究；Xiang Li[2]等利用人工降雨裝置研究了森林植被對(duì)徑流產(chǎn)沙的影響，并表明植被在減少產(chǎn)流量和延緩產(chǎn)流時(shí)間方面有顯著影響；Changjia Li[3]等在模擬降雨條件下對(duì)坡地地表徑流侵蝕進(jìn)行了動(dòng)態(tài)研究。把握模擬降雨中的初始產(chǎn)流時(shí)間是計(jì)算坡面產(chǎn)流產(chǎn)沙量的重要環(huán)節(jié)，而初始產(chǎn)流時(shí)間受到降雨強(qiáng)度、坡度、植被等多種因素的影響[4]。蔣秋玲[5]等對(duì)坡面初始產(chǎn)流時(shí)間進(jìn)行了探究；車明軒[6]等在對(duì)紫色土坡面產(chǎn)流影響的研究中表明了不同雨強(qiáng)和坡度下產(chǎn)流時(shí)間的變化趨勢(shì)；鄧景成[7]開展的灌草植被蓋度對(duì)產(chǎn)流產(chǎn)沙影響的模擬降雨試驗(yàn)中表明增加植被覆蓋率可有效延緩初始產(chǎn)流時(shí)間。

目前，對(duì)產(chǎn)流產(chǎn)沙的研究多為自然形成的土壤基質(zhì)，而針對(duì)廢棄多年無(wú)人管理的尾礦庫(kù)開展的研究較少。廢棄尾礦庫(kù)的表層土質(zhì)疏松孔隙率高，保水、保肥的能力差，且有毒有害物質(zhì)常存在于廢棄尾礦中，植物難以生長(zhǎng)，再加上結(jié)構(gòu)穩(wěn)定性較差，使得廢棄尾礦壩坡面極易發(fā)生徑流及侵蝕產(chǎn)沙。本文采用人工模擬降雨裝置，將廢棄尾礦土作為研究對(duì)象，選擇灰色關(guān)聯(lián)度的方法，研究植被覆蓋率、雨強(qiáng)、坡度對(duì)產(chǎn)流率和產(chǎn)沙率的影響，以及對(duì)初始產(chǎn)流時(shí)間的影響分析。

1 材料與方法

1.1 試驗(yàn)土壤與填土方案

試驗(yàn)?zāi)Ｐ筒捎玫奈驳V土壤取自遼陽(yáng)縣寒嶺鎮(zhèn)一廢棄尾礦庫(kù)，見圖1。

圖1 遼陽(yáng)縣寒嶺鎮(zhèn)一廢棄尾礦庫(kù)

進(jìn)行模型填土?xí)r，將厚度為10 cm的碎石放置在試驗(yàn)槽底部，旨在讓底部試驗(yàn)土壤同現(xiàn)場(chǎng)盡可能相仿，以減少誤差。填土?xí)r，方法為分層填筑，填充和壓實(shí)兩道工序同步進(jìn)行，使土層從上到下土壤的緊實(shí)度增加。在試驗(yàn)槽邊壁進(jìn)行適當(dāng)壓實(shí)，以降低邊壁效應(yīng)使水流正常流動(dòng)。單位容積的土體重量控制在1.8～1.9 t。

1.2 模擬降雨設(shè)計(jì)

人工降雨試驗(yàn)?zāi)Ｐ偷闹黧w規(guī)格為：7.5 m×1.2 m×1.5 m(長(zhǎng)×寬×高)。利用設(shè)置在裝置頂部的降雨架模擬降雨，降雨過(guò)程中的有效高度為6 m，保持恒定水壓，雨滴大小及到達(dá)坡面的速度與天然降雨?duì)顟B(tài)相近。裝置底部設(shè)防滲系統(tǒng)及排水管，末端出水處設(shè)計(jì)有V型集流槽，徑流泥沙樣通過(guò)放置在集流口處的空桶來(lái)收集。導(dǎo)水管上設(shè)置瞬時(shí)電子流量計(jì)和控制閥，降雨裝置噴頭并聯(lián)，通過(guò)打開不同的閥門，控制雨強(qiáng)。

為使每次試驗(yàn)的初始含水率基本維持在10%左右并消除坡體表面的差異性。試驗(yàn)前，在坡面進(jìn)行預(yù)降雨，當(dāng)徑流出現(xiàn)渾濁時(shí)停止降雨。試驗(yàn)中開始記錄降雨的時(shí)間為雨滴落在坡面上的時(shí)刻，記錄產(chǎn)流的時(shí)間以產(chǎn)流出現(xiàn)渾濁的現(xiàn)象為起始，共降雨30 min。開始記錄產(chǎn)流時(shí)即進(jìn)行徑流樣的收集，每次試驗(yàn)共采集15個(gè)樣品。待樣品靜置24 h徑流泥沙樣完全澄清后倒掉上層清液，之后將泥沙樣自然風(fēng)干后，取一部分土壤樣品烘干稱重得到產(chǎn)沙重量并記錄。試驗(yàn)實(shí)物如下圖2所示。

圖2 降雨模型實(shí)物圖

1.3 坡度及植被覆蓋設(shè)計(jì)

遼寧省雨強(qiáng)空間的分布極不均勻，范圍在21～212 mm/h之間，其南部地區(qū)中心雨強(qiáng)值達(dá)到100 mm/h，北部地區(qū)雨強(qiáng)值在21～67 mm/h之間[8]，故本試驗(yàn)設(shè)置20、60、100、220 mm/h四種雨強(qiáng)。坡度因素對(duì)產(chǎn)流產(chǎn)沙的作用最主要[9]。15°是尾礦壩的常見坡度，且楊開寶[10]等在試驗(yàn)中發(fā)現(xiàn)坡度在10°～24°時(shí)侵蝕量隨坡度增加而增加，但當(dāng)坡度大于24°以后侵蝕量同坡度的關(guān)系呈負(fù)相關(guān)。根據(jù)尾礦壩的常見坡度以及試驗(yàn)地的條件，設(shè)置了兩種坡度(15°和30°)，并將坡長(zhǎng)控制在2 m。植被可有效緩解水土流失現(xiàn)象[11]。本試驗(yàn)通過(guò)鋪設(shè)草皮模擬植被覆蓋。大部分相關(guān)研究認(rèn)為超過(guò)70%的覆蓋率，可很大程度上控制水土流失。覆蓋率在60%～80%之間產(chǎn)生的防治效果較在20%～40%之間的穩(wěn)定[12]，根據(jù)當(dāng)?shù)厍闆r，采用0、40%、80%、100%四種植被覆蓋率。

2 結(jié)果與分析

2.1 產(chǎn)沙率與影響因素之間的關(guān)系

當(dāng)坡度和植被覆蓋率一定時(shí)，產(chǎn)沙率原始數(shù)據(jù)可分為8組。每組數(shù)據(jù)中隨著一定的時(shí)間間隔一種雨強(qiáng)對(duì)應(yīng)15個(gè)產(chǎn)沙率數(shù)據(jù)，對(duì)8組數(shù)據(jù)中每一種雨強(qiáng)對(duì)應(yīng)的15個(gè)產(chǎn)沙率數(shù)據(jù)取均值后，共得到32個(gè)數(shù)據(jù)作為灰色分析的原始數(shù)據(jù)。設(shè)置影響因素坡度(X1)、植被覆蓋率(X2)、雨強(qiáng)(X3)為比較數(shù)列，產(chǎn)沙率(X0)為參考數(shù)列。計(jì)算灰色關(guān)聯(lián)度具體步驟如下[13]：

無(wú)量綱化(均值化法)處理原始數(shù)據(jù)即：

(1)

逐一計(jì)算參考數(shù)列(均值化后的第一列)與比較數(shù)列的絕對(duì)差值即：

|x0(k)-xi(k)|

(2)

計(jì)算關(guān)聯(lián)系數(shù)：

ζi(k)=

(3)

計(jì)算關(guān)聯(lián)度：

(4)

據(jù)此，求得影響因素對(duì)于產(chǎn)沙率的關(guān)聯(lián)度，詳見表1。

表1 影響因素與產(chǎn)沙率的關(guān)聯(lián)度

灰色關(guān)聯(lián)度的數(shù)值大小可表明影響因素對(duì)參考數(shù)列的影響程度，數(shù)值越大則影響程度越大。由表1可知，各因素對(duì)產(chǎn)沙率均有較大的影響，對(duì)產(chǎn)沙率的關(guān)聯(lián)度影響最大的因素為雨強(qiáng)，其數(shù)值是0.872 48，雨強(qiáng)的增大致使坡面徑流變大，土壤的滲透性由于徑流沖刷力的增強(qiáng)而下降，致使侵蝕產(chǎn)沙現(xiàn)象加重。坡度對(duì)于產(chǎn)沙率的關(guān)聯(lián)度，為0.862 22，坡度增大會(huì)增加土層被徑流沖走的幾率。植被覆蓋率較另兩種因素的數(shù)值小，為0.845 46，但對(duì)產(chǎn)沙率同樣有較大的影響，地表在植被的阻擋下受到雨水的擊濺侵蝕作用減小，另一方面植被可增加土壤的抗蝕能力，兩者共同作用使產(chǎn)沙率降低。

2.2 產(chǎn)流率與影響因素之間的關(guān)系

把參考數(shù)列Y0設(shè)為產(chǎn)流率，比較數(shù)列不變，利用公式(1)～(3)對(duì)原始數(shù)據(jù)進(jìn)行計(jì)算，得出關(guān)聯(lián)系數(shù)后，利用公式(4)得到影響因素對(duì)于產(chǎn)流率的灰色關(guān)聯(lián)度，計(jì)算結(jié)果見表2。

表2 影響因素與產(chǎn)流率的關(guān)聯(lián)度

對(duì)灰色關(guān)聯(lián)度進(jìn)行排序：雨強(qiáng)>坡度>植被覆蓋率,即對(duì)坡面產(chǎn)流率影響程度相對(duì)較大的首先為雨強(qiáng)，次之是坡度，影響程度相對(duì)最弱的為植被覆蓋率。雨強(qiáng)對(duì)于產(chǎn)流率的關(guān)聯(lián)度為0.742 54，雨強(qiáng)愈大，動(dòng)能愈大，落在坡面的雨滴速度愈大，尾礦壩坡面接受單位時(shí)間內(nèi)的雨量愈大，使得坡面產(chǎn)流率增加，該因素對(duì)產(chǎn)流率的影響最大。坡度對(duì)于產(chǎn)流率的關(guān)聯(lián)度為0.722 97，產(chǎn)流率受其影響較大，坡面徑流沿重力方向移動(dòng)的分力同坡度呈正相關(guān)，所以提高坡度會(huì)加快水流匯集速度，產(chǎn)流率隨之增大。與坡度因素和雨強(qiáng)因素相比，植被覆蓋率對(duì)于產(chǎn)流率的影響較小，關(guān)聯(lián)度為0.670 63。

2.3 尾礦壩坡面初始產(chǎn)流時(shí)間分析

不同試驗(yàn)條件下尾礦壩坡面初始產(chǎn)流時(shí)間的結(jié)果見圖3。由圖3可見，控制植被覆蓋率的條件下，坡度和雨強(qiáng)增加，可縮短尾礦壩坡面的初始產(chǎn)流時(shí)間?？刂破露群陀陱?qiáng)兩種變量，初始產(chǎn)流時(shí)間在植被覆蓋率增加的條件下被推遲。說(shuō)明同初始產(chǎn)流時(shí)間呈負(fù)相關(guān)的為雨強(qiáng)及坡度，呈正相關(guān)的為植被覆蓋率。其中裸地初始產(chǎn)流時(shí)間的最小值為130 s(220 mm/h，15°)，最大值為600 s(20 mm/h，15°)；植被覆蓋率為100%時(shí)初始產(chǎn)流時(shí)間的最小值為165 s(220 mm/h，30°)，最大值為1 080 s(20 mm/h，15°)。雨滴連續(xù)不斷的擊打坡面，使得坡面出現(xiàn)不均勻的雨穴，雨強(qiáng)增大使得雨穴的數(shù)量增加，增大了坡體表面的粗糙程度，局部徑流流速因此延緩，從而提前初始產(chǎn)流時(shí)間。坡度提高，雨水滲到土壤的幾率下降，雨滴加速在坡面上匯集，徑流流速增大，從而使初始產(chǎn)流時(shí)間提前。

圖3 不同試驗(yàn)條件下的初始產(chǎn)流時(shí)間

經(jīng)回歸分析，得出試驗(yàn)條件下初始產(chǎn)流時(shí)間和三種因素的關(guān)系式為：

T=799.981 2+3.079 237D-2.857 14I-9.183 33C

(5)

式中：T—初始產(chǎn)流時(shí)間，s；D—植被覆蓋率，%；I—雨強(qiáng)，mm/h；C—坡度，°。

回歸方程得到樣本決定系數(shù)R2=0.846，觀測(cè)值n=32，達(dá)到極顯著(p<0.01)，說(shuō)明該回歸方程能夠用于預(yù)測(cè)試驗(yàn)條件下的尾礦壩坡面初始產(chǎn)流時(shí)間。

對(duì)初始產(chǎn)流時(shí)間與植被覆蓋率、雨強(qiáng)、坡度做偏相關(guān)性分析。得到三者與初始產(chǎn)流時(shí)間的標(biāo)準(zhǔn)偏回歸系數(shù)分別為0.757、-0.902、-0.559，且植被覆蓋率、雨強(qiáng)、坡度與初始產(chǎn)流時(shí)間的相關(guān)性均達(dá)到極顯著(P<0.01)，由此得出試驗(yàn)條件下的三個(gè)因素都對(duì)尾礦壩坡面初始產(chǎn)流時(shí)間有較大影響且影響極顯著。可確定三個(gè)因素中對(duì)初始產(chǎn)流時(shí)間影響最大的因素為雨強(qiáng)，其次為植被覆蓋率，最后為坡度。

與裸地相比較，植被的存在可推遲初始產(chǎn)流時(shí)間。三種植被覆蓋率對(duì)尾礦壩坡面初始產(chǎn)流時(shí)間的延緩情況如圖4所示，由圖4可知，初始產(chǎn)流時(shí)間能被三種植被覆蓋率有效延緩，其中40%的植被覆蓋率延緩產(chǎn)流時(shí)長(zhǎng)在5～140 s范圍內(nèi)，雨強(qiáng)為100 mm/h、坡度15°時(shí)延緩產(chǎn)流時(shí)間達(dá)140 s，其他條件下延緩產(chǎn)流時(shí)間在5～110 s內(nèi)。80%的覆蓋率延緩產(chǎn)流時(shí)長(zhǎng)在15～300 s之間，較40%覆蓋率的延長(zhǎng)效果好，雨強(qiáng)為20 mm/h、坡度15°時(shí)延緩產(chǎn)流時(shí)間達(dá)300 s。100%的覆蓋率延緩產(chǎn)流時(shí)長(zhǎng)在30～612 s之間，雨強(qiáng)為100 mm/h、坡度15°時(shí)延緩產(chǎn)流時(shí)間達(dá)612 s，其延緩產(chǎn)流時(shí)間較80%、40%覆蓋率效果都要好。這是由于植被的覆蓋減小了雨滴對(duì)于坡面的沖擊，同時(shí)植被也起到阻擋作用，使徑流的速度變慢，覆蓋在地表的植物承受了原來(lái)作用在地表的力。因此，延緩產(chǎn)流時(shí)間隨著植被覆蓋率的增加相對(duì)變長(zhǎng)。

圖4 不同試驗(yàn)條件下的延緩產(chǎn)流時(shí)間

3 結(jié)論

利用人工模型實(shí)驗(yàn)裝置，研究不同植被覆蓋率、雨強(qiáng)和坡度下的廢棄尾礦坡面初始產(chǎn)流時(shí)間以及與產(chǎn)流產(chǎn)沙的關(guān)聯(lián)度，所得結(jié)論如下：

(1)植被覆蓋率、雨強(qiáng)和坡度對(duì)于產(chǎn)沙率的關(guān)聯(lián)順序?yàn)椋河陱?qiáng)>坡度>植被覆蓋率，對(duì)產(chǎn)流率的關(guān)聯(lián)順序?yàn)椋河陱?qiáng)>坡度>植被覆蓋率。當(dāng)雨強(qiáng)、坡度、植被共同作用時(shí)，廢棄尾礦坡面產(chǎn)流產(chǎn)沙率受雨強(qiáng)影響最大，坡度次之，植被覆蓋的影響相對(duì)其他兩個(gè)影響因素較弱。

(2)初始產(chǎn)流時(shí)間與植被覆蓋率、雨強(qiáng)及坡度呈一定的線性關(guān)系，與植被覆蓋率之間表現(xiàn)為正相關(guān)關(guān)系，與雨強(qiáng)及坡度之間表現(xiàn)為負(fù)相關(guān)關(guān)系。裸地、雨強(qiáng)為220 mm/h、坡度為15°時(shí)初始產(chǎn)流時(shí)間最早；覆蓋率為100%、雨強(qiáng)為20 mm/h、坡度為15°時(shí)初始產(chǎn)流時(shí)間最晚。

(3)對(duì)初始產(chǎn)流時(shí)間產(chǎn)生最主要影響的是雨強(qiáng)和植被覆蓋，坡度在試驗(yàn)條件下對(duì)初始產(chǎn)流時(shí)間的影響較另外兩種因素都小，三個(gè)因素對(duì)初始產(chǎn)流時(shí)間的作用效果都較大且均達(dá)到極顯著。試驗(yàn)條件下設(shè)置的三種植被覆蓋率均可有效延緩初始產(chǎn)流時(shí)間，植被覆蓋率的提高對(duì)產(chǎn)流時(shí)間延緩有積極作用。