岩溶区煤矿工程堆积体边坡细沟发育及其水沙关系研究*

林 姿 史东梅† 娄义宝 蒋光毅 彭旭东 李叶鑫

(1 西南大学资源环境学院，重庆 400715）

（2 重庆市水土保持生态环境监测总站，重庆 401147）

（3 贵州大学林学院，贵阳 550025）

（4 辽宁工程技术大学环境科学与工程学院，辽宁阜新 123000）

工矿区水土流失是由于人为扰动地面或堆置固体废弃物而造成水资源和土地资源的破坏和损失[1］。作为我国主要能源之一的煤矿工程在建设和生产运行期间的各种活动会大量占压土地资源并产生各种废弃堆积物，对周边水资源、土地资源、植被资源、水文循环和生态环境质量造成严重影响甚至破坏。煤矿工程建设所产生的大量弃土、石、渣，在水力和重力作用下极易发生诱发性崩塌、滑坡、泥石流等地质灾害[2］。一旦坡面水流将地面切割出一定高差后，沟道两侧的物质势能得到释放机会，形成水力与重力相伴发生的复合侵蚀[3］。当坡面侵蚀形式由面状侵蚀发展为细沟侵蚀时，侵蚀量会成倍或数十倍增长[4-5］，可达到坡面总侵蚀量70%以上[6］。目前对工程堆积体坡面土壤侵蚀过程研究集中在工程堆积体陡坡土壤侵蚀过程[7］、边坡面蚀[8］、弃土场侵蚀产沙过程[9］、边坡土壤流失量的预测[10-11］、煤矿工程建设松散弃渣诱发人为泥石流发育特征[12-13］及工程堆积体水土流失防治等[14-15］，较少涉及重力作用对坡面细沟形态特征的影响。坡面上跌水的形成是侵蚀方式变化的重要分界点，标志着细沟开始形成[16］，重力作用对其形成和坡面产流产沙变化有直接影响。因此本文以煤矿工程堆积体为研究对象，采用野外实地放水冲刷法，重点研究：（1）工程堆积体坡面细沟侵蚀发生特征，分析坡面土壤侵蚀发展阶段特征；（2）重力作用对不同堆积体坡面水沙关系的影响大小及形成原因；（3）重力作用对细沟发育过程的影响。研究结果对于正确认识煤矿工程堆积体坡面产流产沙特征、影响因素及土壤侵蚀量预测模型建立具有重要科学意义，也有助于科学认识重力作用在工程堆积体人为诱发性崩塌、滑坡和边坡稳定性分析中的作用。

1 材料与方法

1.1 试验材料

工程堆积体来源于重庆市5大矿区之一的天府煤矿，该矿区处于喀斯特岩溶区；重庆市喀斯特岩溶区分布面积可达3×104km2，在该地区分布有大量锰矿、煤矿、石灰石等建设项目。当地多年平均降雨量在1 000 mm以上，工程堆积体在暴雨条件下极易发生诱发性崩塌、滑坡等严重土壤侵蚀现象。采样区成土母岩主要是石灰岩、泥灰岩，土壤主要是暗紫色水稻土、暗紫泥土、矿子黄泥土、冷沙黄泥土等，工程堆积体由土壤、不同径级的石灰岩及泥灰岩碎屑及块石所组成的岩土混合物，这种物质组成在喀斯特矿区普遍存在。调查发现，其边坡坡度多在25.5°～38°之间变化，为客观反映不同堆放坡度条件下煤矿工程堆积体边坡的侵蚀特征，本文设计了四种不同坡度25°、30°、35°、40°工程堆积体下垫面开展试验研究。将同一批煤矿工程堆积体均匀分层填入4个不同坡度（25°、30°、35°、40°）试验小区，试验小区坡长10 m，宽1 m。各小区的布置情况见表1。试验所用工程堆积体容重约为1.23 g·cm-3，总孔隙度为52.96%，毛管孔隙度为13.45%，非毛管孔隙为39.31%，黏聚力为9.82 kPa，内摩擦角为29.68°。每次放水冲刷试验前，用铁耙将小区坡面平整，同时在小区内采集工程堆积体土壤样品，测定其物理性质，试验小区基本情况如表1所示。

1.2 野外放水冲刷试验设计

试验于2013年7—8月在西南大学生产建设项目水土流失定位试验基地进行（图1）。根据当地暴雨历时以及暴雨频率在试验小区上所产生的单宽流量特征，为揭示流量对煤矿工程堆积体边坡细沟发育特征的影响，本文放水流量采用10、15、20、25、30 L·min-1五个梯度，冲刷时间设计为60 min。

小区顶端放置1个簸箕型溢流槽，以保证水流以薄层形式均匀向下流动。水流经稳压水泵及流量计后流入溢流槽，并通过溢流槽上端阀门控制流量。试验前，率定放水流量若干次，保证放水流量与设计流量误差不超过5%。记录产流时间，产流10 min内每1 min接取1次径流泥沙样，10 min后每隔3 min接取1次泥沙样，用烘干法（105℃）测定泥沙量。流速的测定是在坡面上、中、下断面中间位置，采用染色法测定坡面表层最大流速并根据径流流态进行修正（过渡流修正系数为0.7，紊流修正系数为0.8）。流宽、流深采用直尺测量，分上、中、下三个断面测定。

2 结 果

2.1 工程堆积体坡面细沟侵蚀发生特征

细沟是坡面侵蚀的一种类型，细沟形成是坡面土壤侵蚀程度加剧的主要标志，也是造成坡面产流产沙量变化的主要原因[17］；由于径流汇集面积和流速不断增大，坡面水流冲刷能力大大增加，引起工程堆积体坡面不均匀凹陷，从而在地表上逐渐形成细沟，对工程堆积体细沟发育程度的影响因不同下垫面类型而有很大差异。由不同坡度条件下工程堆积体坡面侵蚀过程的几个关键时间可见（表2），工程堆积体坡面产流时间随放水流量增大而缩短且坡度越大其产流时间越短，细沟出现时间越短；在各场冲刷试验中，产流时间最迟（2 558 s）发生在25°坡面、10 L·min-1径流冲刷条件下，而产流时间最早（76s）发生在40°坡面、25 L·min-1径流冲刷条件下，这种现象主要与前期土壤含水率大小有关。各种工程堆积体坡面细沟出现时间也与放水流量、坡度呈负相关关系，即放水流量越大、坡度越大，细沟出现时间越短；细沟最快出现为4s、最慢为97s，二者时间相差24倍，这主要因为冲刷径流量越大、则坡面泥沙流失越快，细沟出现时间越短。微地形（坡度）条件也是决定坡面细沟下切快慢和发育程度的原因之一，坡度越大、则重力作用对边坡细沟发育影响越大、细沟下切越快，出现时间越短。

表1 试验小区基本情况Table 1 Basic information of the experimental plots

图1 工程堆积体实验装置及冲刷过程示意图Fig. 1 Schematic of the field scouring experiment apparatus and its scouring process on engineering waste piles

工程堆积体坡面细沟形成后，由于细沟壁土体崩塌、脱落、下滑的随机性，使得坡面径流含沙量相应增大且呈现波动性变化趋势（图2）。不同坡度条件下工程堆积体坡面含沙量波动变化随冲刷时间呈现变化幅度减小趋势，工程堆积体坡面侵蚀过程可分为片蚀和细沟侵蚀2个阶段。

（1）砂砾化面蚀阶段，即坡面开始产流至细沟形成时的坡面侵蚀过程。该阶段发生在产流后3 min以内，此时坡面表面形态被径流侵蚀破坏程度小，未形成跌坎和细沟[18］，随着径流侵蚀力增大达到足以剥离和分散坡面土壤时便发生侵蚀，这对坡面后续跌坎及细沟形成产生直接影响。

（2）细沟状面蚀阶段，即坡面跌坎形成细沟后的坡面侵蚀过程，可概括为细沟扩张、过渡和稳定3个发展阶段。细沟扩张阶段主要发生在产流3～24 min之间，该阶段坡面流沿坡面低凹处汇集，造成堆积体坡面集中冲刷现象并开始形成细沟；同时由于在坡面细沟发育早期阶段，细沟沟头崩塌、沟壁崩塌等重力作用过程严重，造成坡面含沙量剧烈波动变化现象，表现出明显的水力—重力复合侵蚀特点。过渡阶段主要发生在产流24～30 min之间，此时径流冲刷形成的水力侵蚀作用减弱，坡面细沟长宽深基本保持稳定，细沟沟壁重力作用现象减弱，这造成径流含沙量逐渐降低。稳定阶段则发生在产流30 min后，细沟沟壁重力作用成为该阶段泥沙的主要来源，坡面径流含沙量较小。

表2 不同工程堆积体坡面细沟侵蚀过程特征Table 2 Characteristics of rill erosion process on slopes of the piles different in slope gradient

2.2 重力作用对堆积体坡面水沙关系的影响

重力作用是影响工程堆积体坡面产流产沙的重要因素。坡面细沟壁土体崩塌、滑落是造成工程堆积体坡面产流产沙过程中存在波动现象的重要原因，也是坡面不稳定的重要因素[19］。分析坡面产流产沙变异系数及重力作用对总产沙量贡献的关系（表3）可知，平均产流率与放水流量呈正相关关系（R2= 0.94），而与坡度关系不明显；坡面产流率变异系数总体上随着放水流量增大而增加，变异系数最大（374. 3%）为40°坡面、25 L·min-1径流条件，最小（7.774%）为35°坡面、10 L·min-1径流条件，二者相差48倍以上；主要原因在于当径流条件增大时，坡面细沟冲刷强烈，造成细沟沟壁崩塌，崩塌体在短时间内会堵塞水流；而当径流汇集达到一定程度后，细沟内崩塌体被迅速冲开，这使得坡面径流量在短时间内发生剧烈的波动变化。工程堆积体坡面径流含沙量总体上随放水流量增大而增加，在各次冲刷条件下径流含沙量最大为0.480 kg·L-1，最小为0.020 kg·L-1，二者相差24倍，这种现象主要是坡面水力侵蚀和细沟重力作用共同造成的。

图2 不同工程堆积体坡面含沙量变化特征Fig.2 Variation of sediment yield on the slope of the pile relative to feature of the pile

工程堆积体坡面径流含沙量变异系数在一定程度可反映重力对细沟作用的程度，变异系数越小则说明在坡面侵蚀过程中径流含沙量比较稳定，重力对细沟作用较少或没有发生；反之，重力对细沟作用严重。由图3可知，在不同坡度坡面条件下，径流含沙量变异系数均随着放水流量增大呈现先增加、后减小的变化趋势，且20 L·min-1放水条件下径流含沙量变异系数将达到最大值，这表明重力作用明显影响堆积体坡面细沟发育的临界径流条件为20 L·min-1。当放水流量由10 L·min-1增加至20 L·min-1时，重力作用对坡面细沟程度逐渐增加并达到最大值，此时坡面径流含沙量以重力作用产沙为主、水力侵蚀产沙为辅；当放水流量由20 L·min-1增加至30 L·min-1时，重力作用对坡面细沟程度减弱，此时坡面含沙量以水力侵蚀和重力对细沟沟壁综合作用为主。在不同放水流量作用下，径流含沙量变异系数随坡度增加也表现出先增加、后减小的变化特点，且均在35°坡度条件下，径流含沙量变异系数达到最大值，这表明重力作用明显影响堆积体坡面细沟发育的临界坡度条件为35°。因此，生产建设项目区工程堆积体水土保持措施布置时应关注工程堆积体边坡的坡度大小，避免在暴雨诱发条件下人为崩塌、滑坡、泥石流的发生。

表3 坡面产流产沙变异系数及重力作用对总产沙量的贡献Table 3 Variation coef fi cient of runoff and sediment yields and contribution of gravity to the total sediment yield on the slopes of the piles

图3 径流含沙量变异系数与放水流量和坡度的关系Fig. 3 Relationship of variation coefficient of the sediment concentration in runoff with discharged flow rate and slope gradient

在工程堆积体坡面细沟长、宽、深发育过程中，重力作用对细沟侵蚀发生具有很大的偶然性，很难及时观测和定量测定由水力作用单纯主导的临界侵蚀产沙量。因此本文采用细沟发育稳定的径流含沙量[21］（即每次冲刷试验中最后10 min的平均径流含沙量）作为水力侵蚀含沙量上限值和重力作用对细沟侵蚀含沙量的下限值。由表3可见，在各次冲刷试验中重力作用对坡面产沙量贡献最大为99.60%，最小为17.41%，这表明重力作用对工程堆积体坡面侵蚀是影响坡面产沙的重要因素，同时也是导致坡面径流含沙量波动变化的重要原因。

2.3 重力作用对堆积体坡面细沟形态发育的影响

相关研究表明，重力作用是造成坡面细沟由浅变深、由窄变宽的主要原因之一[19,22］，此处采用重力作用产沙贡献表明重力作用对细沟形成发育的影响。由表4可知，在不同放水流量、坡度条件下，工程堆积体坡面细沟发育过程和特征差异明显。坡面细沟平均沟宽、沟深均随放水流量增大而增加，且坡度越大则细沟宽度、深度越大；这是由于在坡面细沟发育过程中，水流侵蚀会加剧沟底下切，使细沟变深，同时重力对细沟沟壁的作用会加速沟岸扩张，使细沟变宽。在各次冲刷试验中，细沟沟宽最小7.89 cm，最大19.73 cm；细沟沟深最大6.73 cm，最小2.17 cm，细沟的沟深、沟宽极值均相差3倍左右，这说明在堆积体坡面侵蚀过程中，径流作用对沟深发展起主导作用，重力作用对沟宽发展起主导作用且两种作用程度相当，所造成的侵蚀产沙量均为坡面径流泥沙的主要来源。

表4 不同工程堆积体坡面细沟发育特征Table 4 Characteristics of rill development in the engineering accumulation slopes

在相同径流条件下，坡面细沟密度越大则相应的侵蚀产沙量也越大；重力作用坡面产沙贡献随坡度和放水流量增加呈波动趋势变化。在坡面细沟发育过程中，由于冲刷初期细沟长、宽、深和细沟密度均较小（1.35～1.90 m·m-2），水流冲刷形成的产沙量（1.74～2.39 kg）高于细沟沟壁重力坍塌形成的产沙量（0.46～1.77 kg）；在冲刷后期由于细沟长、宽、深的扩大和细沟密度变大（2.07～3.00 m·m-2），重力作用对细沟沟壁产沙量增加至149.00～162.7 kg，细沟沟壁崩塌泥沙可能堵塞水流，造成细沟水流流路发生改变使坡面细沟呈现S型发育。细沟沟壁崩塌主要与沟道下切、侧蚀以及由于土壤含水率增大所致，造成的土壤黏聚力降低及土壤膨胀现象有关。细沟宽深比反映了坡面细沟下切程度及细沟沟壁发生崩塌等重力作用可能性，工程堆积体坡面细沟宽深比最大为4.36，最小为2.12，二者相差2倍；细沟宽深比随放水流量和坡度增大而减小，说明坡度增大会加剧沟底下切、增大重力作用对细沟沟壁侵蚀发生可能性。

3 讨 论

坡面侵蚀过程及坡度对工程堆积体边坡稳定性有重要影响。相关研究表明，坡面径流侵蚀对边坡破坏形式主要包括边坡浅部冲刷、淘蚀破坏及深部入渗溶蚀破坏3种形式[23-24］。在一次径流冲刷后，工程堆积体坡面变得更加粗糙，这使得径流流路更为集中，径流动力条件不均匀程度迅速增加；因此当降雨径流再次冲刷时，重力作用对坡面细沟侵蚀和沟壁崩塌等侵蚀更剧烈，造成了松散堆积坡面不稳定程度加剧。同时坡面入渗水流会促使边坡内部土体含水率不断增大，而入渗水的溶蚀作用会对土体内部可溶性盐类和胶结物质产生化学侵蚀，这在一定程度上也增加了边坡不稳定性。同时工程堆积体中块石会直接影响堆积体坡面侵蚀类型和形式，也会导致重力作用对坡面侵蚀和水力侵蚀贡献率差异性。坡度作为影响坡面侵蚀发生与否的地形因子，可直接决定坡面土体稳定状态。坡度越大则坡面土体所受到的下滑力越大，最大静摩擦力越小，同时坡度越大径流流速越大，侵蚀潜能增加[25-27］。在不同坡度条件下，侵蚀发生形式及侵蚀量存在较大差别。根据本文研究，当放水流量为15、20、25 L·min-1时，35°坡面的径流含沙量变异系数最大、重力作用对坡面产沙贡献最大（图4），这表明小于35°的坡面以水力侵蚀产沙为主，而大于35°坡面以重力作用产沙为主；在最小（10 L·min-1）和最大放水流量（30 L·min-1）条件下，35°坡面径流含沙量变异系数较小、重力作用对坡面产沙贡献也较小，初步判定，该流量条件下径流侵蚀对坡面含沙量贡献最大的临界坡度为35°。王文龙等[28］对5°、15°、25°及35°条件下坡面侵蚀过程的研究也表明，坡沟系统在5°、15°和25°坡度条件下，坡面土壤侵蚀方式以溅蚀、片蚀和沟蚀为主，而大于或等于35°条件时，坡面土壤侵蚀方式以切沟为主，且有崩塌、滑塌等侵蚀发生。研究均表明在坡面细沟侵蚀过程中，细沟沟壁崩塌等侵蚀发生不仅取决于微地形条件，还决定于细沟被径流冲刷淘蚀的程度，因此今后应加强重力与降雨、径流、微地形（坡度和坡长）、坡面细沟发育形态之间的定量研究，不同土石比工程堆积体坡面径流连续冲刷过程及其边坡稳定性的模拟研究，为重力作用对坡面侵蚀过程模拟研究提供基础数据。

土壤抗剪强度是影响坡面侵蚀发生的重要原因之一，它决定了坡面土体内部力学稳定性。研究表明，坡面土体稳定性取决于土壤摩擦阻力和土壤内聚力[29］，而土壤含水率及密实度则是影响重力作用的主要因素[30］。根据本文研究，工程堆积体土壤黏聚力在16.18%含水率时可达最大14.64 kPa，随含水率增大呈先增加后减小变化趋势，而内摩擦角则随含水率增加而减小（图5）。重力作用明显影响堆积体坡面细沟发育的临界坡度条件为35°，其数值大于工程堆积体内摩擦角（30.79°）。相关研究也表明，土壤抗剪强度和水分含量呈反比关系[31］。当坡面土壤含水率增大时，不仅会增大土体沿坡向的重力分力，而且会减小土体黏聚力和内摩擦角，使土体或土层处于失稳状态。因此，根据上述坡面土体抗剪强度和含水率量化关系，在雨季及时采取边坡防护措施，以提高工程堆积体稳定性，防止暴雨诱发性崩塌、滑坡、泥石流灾害发生的可能性。植被在提高边坡土壤的抗侵蚀能力、增强土壤的稳定性中起着至关重要的作用，细根、中等直径的根系分泌物对土体颗粒起到黏结作用，从而提高土壤的黏聚力c，进而提高土壤的抗剪强度[32］，今后应加强植物根系、土壤含水率、土壤紧实度与土体抗剪强度相互作用关系研究，寻求减弱重力作用对坡面侵蚀有效途径。

图4 重力作用产沙贡献与坡度的关系Fig.4 Relationship between contribution of gravity and slope gradient

图5 土壤含水率与黏聚力和内摩擦角的关系Fig.5 Relationships of soil water content with cohesive strength and internal friction angle

4 结 论

工程堆积体坡面侵蚀过程可分为面蚀（产流3 min内）和细沟侵蚀阶段，其中细沟侵蚀存在扩张（3～24 min）、过渡（24～30 min）和稳定（30 min后）3个发展过程。工程堆积体坡面细沟出现时间与放水流量、坡度呈负相关关系，细沟出现时间最快为4 s，最慢为97 s，二者约相差24倍。重力作用是影响工程堆积体边坡产流产沙及波动变化的重要原因。在坡面侵蚀过程中产流率变异系数在7.77%～374.3%之间，产沙率变异系数在1.75%～1021%之间；重力作用对工程堆积体边坡总产沙量贡献在17.41%以上，最高可达99.60%；初步判定，重力作用明显影响堆积体边坡细沟发育的临界径流条件为20 L·min-1，临界坡度条件为35°。径流作用主导细沟深度发展，而重力作用主导细沟宽度发展，两种主导作用对工程堆积体边坡细沟发育影响程度相当。工程堆积体边坡细沟沟宽在7.89～19.73 cm之间，沟深介于2.17～6.73 cm，宽深比介于2.12～4.36，细沟密度介于1.35～3.00 m·m-2；边坡细沟平均沟宽、沟深均随放水流量增大而增加，且坡度越大、细沟宽度、深度越大；在相同径流条件下，侵蚀产沙量随边坡细沟密度增大而增大。