郭雁斌，張傳偉，田 光，劉 宇，何以劍，李兆霖，4

（1.國家能源集團(tuán) 神華哈爾烏素露天煤礦，內(nèi)蒙古 鄂爾多斯 010300；2.昆明煤炭設(shè)計研究院有限公司，云南 昆明 650000；3.中國礦業(yè)大學(xué) 礦業(yè)工程學(xué)院，江蘇 徐州 221116；4.中國礦業(yè)大學(xué)煤炭資源與安全開采國家重點實驗室，江蘇 徐州 221116）

礦山生態(tài)修復(fù)作為我國生態(tài)保護(hù)修復(fù)的重要領(lǐng)域之一，在“雙碳目標(biāo)”的新背景下，對相關(guān)煤炭生產(chǎn)企業(yè)的高質(zhì)量、高標(biāo)準(zhǔn)的企業(yè)快速轉(zhuǎn)型提供了全新的理論依據(jù)和及指導(dǎo)方案[1-2]。如何實好轉(zhuǎn)型期企業(yè)發(fā)展與生態(tài)修復(fù)并行策略成為礦山可持續(xù)發(fā)展的關(guān)鍵所在。露天開采作為我國草原區(qū)域的重要煤炭開采方式，多分布于新疆、內(nèi)蒙古等生態(tài)較為脆弱區(qū)域[3]，植被覆蓋率較低，其排土場邊坡受自然環(huán)境及人為擾動影響，受侵蝕程度較大，水土流失現(xiàn)象較為嚴(yán)重[4]。

目前國內(nèi)外學(xué)者針對礦山邊坡生態(tài)修復(fù)方面已經(jīng)進(jìn)行了諸多研究，主要包括：不同影響因素下的土壤侵蝕特征[5-7]；煤巖滲流特征[8-10]；土壤理化性質(zhì)重構(gòu)技術(shù)研究[11-13]；水土流失評價模型[14]等方面。但是傳統(tǒng)階梯式排棄方案下的生態(tài)修復(fù)措施相較于礦區(qū)周邊自然環(huán)境仍存在一定差異性，難以從根本上解決邊坡水土流失情況。因此，在此基礎(chǔ)上有學(xué)者提出仿自然地貌重塑技術(shù)[15-17]，通過“引導(dǎo)型生態(tài)修復(fù)技術(shù)”，結(jié)合不同礦區(qū)實際自然邊坡分布情況布置排土場，能有效減緩邊坡的水土流失狀況，降低坡面受侵蝕程度。為此，基于反S 型近自然邊坡修復(fù)模型[18]，以我國東部草原區(qū)某大型露天礦內(nèi)排土場為例，對不同高度臺階重塑后邊坡表面的侵蝕特征變化情況、滲流路徑變化情況以及孔隙水壓力重分布情況進(jìn)行模擬分析。

1 排土場邊坡侵蝕機理

1.1 降雨侵蝕機理

我國東部草原區(qū)年平均降雨量較少，且多為短時強降雨，因此排土場邊坡坡面受水力侵蝕現(xiàn)象較為嚴(yán)重，其中水力侵蝕過程主要分為濺蝕、地表徑流及水流下滲3 個過程。通過不同影響方式造成坡面及坡體內(nèi)部土壤等物質(zhì)發(fā)生不同程度的搬運、剝蝕以及破壞[19-20]。由于露天礦排土場通常呈現(xiàn)階梯式布置，因此臺階坡面相較于水平面傾斜程度較大，受重力作用以及降雨沖刷作用影響，大量土壤沉積在臺階坡頂線與上級臺階坡底線交界處。且因為降雨強度較大，在坡面各處形成大小不同的侵蝕溝，使得坡面降水難以充分入滲，坡頂部水土大量流失，通常在坡面呈現(xiàn)出臺階底部植被生長旺盛而上部植被較為稀少的生長情況，整體復(fù)墾程度較低。

1.2 土體抗蝕性特征

邊坡抗侵蝕特征分析如圖1。

圖1 邊坡抗侵蝕特征分析Fig.1 Analysis of anti-erosion characteristics of slope

除降雨所產(chǎn)生水力侵蝕外，邊坡土體自身抗蝕性特征也很大程度上決定了礦山土體復(fù)墾效果以及最終水土流失程度。對圖1 反S 型復(fù)合重塑邊坡坡體不同位置土體分布情況進(jìn)行理論分析，選取邊坡不同高度范圍內(nèi)的3 個區(qū)域作為樣本進(jìn)行比較，其中區(qū)域Ⅰ及區(qū)域Ⅲ范圍內(nèi)的臺階坡面角β 與傳統(tǒng)單一平直坡面相比較為平緩，重力及降雨沖擊力所產(chǎn)生的下滑力分支受表層土壤阻水性影響較大；區(qū)域Ⅱ為坡頂部下降至坡底部的過渡階段，其區(qū)域內(nèi)的臺階坡面角β 較傳統(tǒng)單一平直坡面相對較陡，但因其范圍內(nèi)土壤下降高度與區(qū)域Ⅰ相比較為有限，因此僅有少部分土壤隨降雨滑落至區(qū)域Ⅲ內(nèi)以及坡底與下級臺階交界處。

其中，受邊坡土體理化性質(zhì)、降雨強度以及坡度等因素影響，坡面不同高度徑流流速的計算公式為：

式中：V 為流速，m/s；q 為單寬流量，m2/s；J 為臺階坡度；V0為徑流啟動速度，m/s；D 為土體顆粒直徑，mm；ρs為土體顆粒密度，g/cm3；a、K 為系數(shù)；m、n為指數(shù)。

隨著邊坡自身抗侵蝕能力增強，坡面流速降低，徑流通路形成時間增加，雨水在坡面上部區(qū)域滯留時間增加，表面降雨有更充足時間下滲至坡體內(nèi)部，從而有效地減緩了因入滲時間不足而導(dǎo)致的頂部區(qū)域內(nèi)的水土流失。

2 邊坡巖體參數(shù)

1）自然邊幫巖體參數(shù)。根據(jù)現(xiàn)場實際邊坡勘察研究，結(jié)合當(dāng)?shù)氐刭|(zhì)勘探報告及滑坡防治技術(shù)研究報告進(jìn)行巖體物理力學(xué)試驗，確定的邊幫巖體自然狀態(tài)下主要力學(xué)指標(biāo)見表1。其中可采煤層以5#、6#煤層為主，6#煤層底板多為泥巖、砂巖，整體強度較低。內(nèi)排土場基底無明顯弱層存在，其產(chǎn)狀特征對內(nèi)排土場穩(wěn)定性影響較小。

表1 自然邊坡巖體力學(xué)參數(shù)Table 1 Mechanical parameters of natural slope rock mass

2）內(nèi)排土場巖體參數(shù)。通過試驗發(fā)現(xiàn)，礦區(qū)排土物料整體黏聚力較大，受排棄方法影響：臺階上部5 m 范圍內(nèi)巖體受車輛動載荷，壓實效果明顯，5 m 范圍內(nèi)巖體整體密實度較高。對其進(jìn)行重新取樣測量后得到的內(nèi)排土場巖體物理力學(xué)參數(shù)見表2。

表2 內(nèi)排土場巖體力學(xué)參數(shù)Table 2 Rock mechanics parameters of inner dump

3 重塑邊坡降雨模擬分析

3.1 理論計算模型

結(jié)合該礦區(qū)實際生產(chǎn)情況，設(shè)計模擬所選理論模型參數(shù)為內(nèi)排土場邊坡高度180 m，邊坡角度約10°。分析分別選取單臺階高度為15 m 及20 m 時的內(nèi)排土場邊坡分布情況進(jìn)行數(shù)值模擬計算。其中臺階坡面角18°，平臺寬度32 m。經(jīng)計算確定單臺階高度為20 m 時，凸面曲率為7.46×10-3，凹面曲率為3.39×10-3；單臺階高度為15 m 時，凸面曲率為7.91×10-3，凹面曲率為3.72×10-3。凸面水平占比約46.8%，凸面豎直占比約38%。

3.2 重塑邊坡滲流及孔隙水壓力特征

我國東部草原區(qū)整體降雨較少，且多集中在6—8 月，因此模擬選用降雨較為集中的雨季進(jìn)行特征計算，選取降雨強度最大的30 min 進(jìn)行計算，其所產(chǎn)生的侵蝕力R 可通過式（3）計算得出。由于該地區(qū)地震烈度較低，因此根據(jù)規(guī)范，不考慮地震對邊坡穩(wěn)定性的影響。

式中：E 為單次降雨動能，J/（m2·mm）；I30為30 min 的最大降雨強度，mm。

選用GEO 軟件中的SLOPE/W 模塊以及SEEP/W 模塊進(jìn)行擬合分析，結(jié)合實際情況對15 m 及20 m 2 種不同形況下的排土場邊坡臺階進(jìn)行重塑，通過對比分析重塑前后內(nèi)排土場邊坡內(nèi)部的滲流路徑變化以及孔隙水壓力重分布情況。滲流路徑分布圖如圖2，孔隙水壓力分布圖如圖3。

圖2 滲流路徑分布圖Fig.2 Diagrams of seepage path distribution

圖3 孔隙水壓力分布圖Fig.3 Pore water pressure distribution

由圖2 可知：在相同降雨條件下，圖2（a）中滲入坡體內(nèi)部的降雨量整體較小，且下滲的降雨量均有向排土場坡面方向的運動趨勢，不同高度臺階坡面均出現(xiàn)不同程度的水流滲出現(xiàn)象，滲入坡體內(nèi)部的降雨量較少。受我國東部草原區(qū)域氣候較為干旱，蒸發(fā)量較高等因素影響，留在邊坡內(nèi)部的降雨量十分有限，因此復(fù)墾效果通常會受到影響；與原內(nèi)排土場邊坡相比，圖2（b）、圖2（c）重塑后的邊坡滲入坡體內(nèi)部的降雨量明顯增加，其中20 m 臺階重塑后的滲流密集程度更大，且整個內(nèi)排土場坡面沒有出現(xiàn)明顯的水流滲出情況；這主要是因為經(jīng)過重塑后的邊坡，坡面抗蝕性提高，降雨難以快速形成徑流通路，從而無法沿通路快速運移至邊坡底部，降雨在坡面停留時間增加，整體入滲程度較高。

由圖3 可知：圖3（b）、圖3（c）中的排土場邊坡在重塑后靠近邊坡上部區(qū)域均出現(xiàn)孔隙水壓力重分布現(xiàn)象，負(fù)孔隙水壓力分布范圍擴大，局部土壤基質(zhì)吸力增加；過渡區(qū)域面積較為均勻，上部邊坡未出現(xiàn)明顯的因巖層下滲速度不同而導(dǎo)致的積水現(xiàn)象，從而減緩了坡面滲水情況的出現(xiàn)；靠近邊坡底部以及基底等深度較深區(qū)域受下滲作用影響較小，其孔隙水壓力分布與重塑前的內(nèi)排土場相比沒有表現(xiàn)出明顯的差異性。

4 結(jié) 語

1）降雨過程中，受水力侵蝕作用影響，原內(nèi)排土場邊坡易在坡面處形成大小不一的侵蝕溝道，進(jìn)而形成徑流通路，導(dǎo)致大量坡面土壤隨水流快速運移至邊坡底部。進(jìn)行重塑后的坡面降低了單級臺階上部及下部區(qū)域臺階坡面角，增加坡面阻水性能，降低了形成徑流通路所需時間，有效地減緩了邊坡的水土流失。

2）通過對相同降雨強度下的內(nèi)排土場重塑前后邊坡進(jìn)行數(shù)值模擬發(fā)現(xiàn)，原內(nèi)排土場邊坡內(nèi)部下滲水流多集中在排土場巖層區(qū)域，且在不同臺階處均表現(xiàn)出不同程度的滲水現(xiàn)象。當(dāng)重塑高度為15 m時，坡體內(nèi)部滲流密集程度增加，其中水平方向滲流趨勢下降，豎直方向滲流趨勢明顯增加，滲流過程多集中在原始內(nèi)部，內(nèi)排土場坡面未出現(xiàn)明顯水流滲出現(xiàn)象；當(dāng)重塑高度為20 m 時，滲流密集程度進(jìn)一步增加，水平方向滲流趨勢進(jìn)一步下降。

3）經(jīng)過重塑后的邊坡孔隙水壓力實現(xiàn)重分布，邊坡上部負(fù)孔隙水壓力分布范圍增加，受雨水下滲影響靠近邊坡上部局部土壤處于過飽和狀態(tài)，整體土壤飽和度較高。因此，出于邊坡穩(wěn)定性角度考慮，建議使用方案Ⅰ作為實際生產(chǎn)過程中的實施方案對內(nèi)排土場邊坡進(jìn)行重塑，從而能夠在保證安全的前提下，提高邊坡抗侵蝕能力以及保水特性，有效緩解因干旱以及強降雨氣候所造成的礦山邊坡水土流失。