风蚀是在风力作用下地表物质被释放搬运的过程，是我国西北干旱区最主要的地表侵蚀类型^[1-2]。风蚀一方面使地表物质被搬运，破坏土壤理化性质，使土壤养分流失，造成土地荒漠化；另一方面风蚀释放的地表颗粒物进入大气，造成环境污染及沙尘暴的频发^[3-4]。研究表明^[5-6]，风蚀主要与土壤自身物理特性、下垫面及环境状况有关，包括土壤含水率、土壤机械组成、土壤团聚体、地表覆盖度、地表紧实度、环境风速、环境湿度等。目前，国内外学者主要进行农田、退化草地、固定半固定沙漠等土壤的风蚀研究，关注于风蚀率及颗粒分布，风蚀过程及动力机制，风蚀影响因子及模型，抗风蚀措施及效果等方面，取得了较多研究成果^[7-15]。但对于大量存在的矿业废弃尾矿风蚀及其颗粒释放研究较少。尾矿是矿业生产过程中的产物，是重要的沙尘源和环境污染源^[16]。尾矿无论是干式排放、湿式排放还是膏体式排放，所形成的尾矿库都会对周边土壤及地下水造成不同程度的污染^[17]，同时尾矿库干涸后形成的干滩在风力作用下会产生尾矿风蚀。尾矿风蚀颗粒表面极为粗糙，含有原子价不饱和的离子，人若过量吸入，会导致病患，尾矿风蚀也是造成环境空气质量恶化的局部污染源之一^[18-19]。因此，尾矿风蚀对环境空气质量和人体健康的影响绝对不能低估^[3]，其风蚀率及颗粒释放特征等研究是尾矿风蚀及其相关防治研究的基础。研究结果有助于尾矿采取有效的防护措施，控制尾矿风蚀污染。

为了解尾矿风蚀特征及其对环境的影响，合理采取防治措施，本文通过风洞试验对西北干旱区尾矿风蚀特征进行研究，探讨不同风速下含水率和结皮量对尾矿风蚀的影响，研究结果对于揭示我国西北干旱区尾矿风蚀及其颗粒释放规律，拓展我国风蚀研究内容，探索我国西北干旱区尾矿风蚀防护措施具有重要意义，并可为科学制定尾矿库风蚀防护措施等提供理论依据。

1 材料与方法 1.1 研究区概况

研究所用尾矿采自于喀拉通克铜镍矿加乌尔尾矿库，该尾矿库位于新疆富蕴县西南25 km(89°41′E，46°48′N，海拔931 m)。尾矿库所在区域属于典型的干旱大陆性气候，夏季炎热干旱，冬季严寒少雪，根据当地气象资料显示：该区年平均气温2.86 ℃，负温期长达5个月(11月至次年3月)，年最高温集中在7—8月；年平均相对湿度为64%；年平均降水量195.8 mm，多集中在5—10月，有效降水量为41.6 mm，年蒸发量1 960.6 mm；春末夏初多风，风向以西北、东北为主，年平均风速2.7 m/s。

尾矿库附近地形平缓，植被较少且以骆驼刺为主。库内尾矿表层在干涸过程中会逐渐固化，形成一层结构松散但盖度较高的物理结皮。库内尾矿化学成分主要以SiO₂、Al₂O₃、MgO、Fe为主，含量约为76.58%。库内尾矿表层物质机械组成(表 1)不仅含少量黏土(1.63%)等细颗粒，还含有大量中沙(11.01%)及以上粗颗粒，但主要以极细沙和细沙为主，含量分别为30.92% 和30.50%，其次为粗粉沙，含量15.30%。平均粒径3.349φ(0.098 mm)，属于极细沙，分选性较差，属于正偏窄峰态。

1.2 试验设计 1.2.1 样品采集与处理

于2019年6月在尾矿干滩上随机采集6处50 cm×50 cm，表层20 cm的尾矿样品，分别装入样品箱，带回室内进行风洞试验及粒度分析。

粒度分析在中国科学院新疆生态与地理研究所荒漠环境与生态修复实验室完成，采用BT-2001型激光粒度分布仪(湿法)，分析结果以体积百分含量表示，根据Wentworth^[20]对粒级的划分标准统计了尾矿样品的颗粒级配，并根据Folk-Ward图解法^[21]计算粒度参数(平均粒径M_Z、分选系数σ_I、峰度K_G、偏度Sk_I)。

风洞试验前将采集的6处尾矿样品混合均匀，晾干，过2 mm筛，去除杂质备用。由于所研究尾矿为湿式排放，结皮为物理结皮，因此本试验采用混合法制备尾矿含水，喷水法制备尾矿结皮。由混合法和喷水法制备的样品含水和结皮情况与尾矿库排放的湿式尾矿及其干涸后形成的结皮较为相似，还原度较高，具有良好的代表性。

含水率制备方法：试验设置4个含水率梯度，分别为对照组、2%、3%、4%，对照组为自然晾干样品，含水率0.3%，其余每个梯度准备足够量的晾干样品，并根据样品量计算所需水量(试验用水为自来水)，将计算好的水量分别加入各样品中，混合均匀，密封保存并静置12 h。试验前将各梯度样品装入长57 cm、宽37 cm、深5 cm的沙盘中，铺满、表面刮平，并重新测量含水率，真实含水率梯度为1.74%、2.72%、3.57%。含水率测量采用烘干法，温度控制在120 ℃，恒温烘干8 h。

结皮制备方法：试验设置4个结皮量梯度，分别为对照组，25、50、100 g，对照组为自然晾干样品，结皮量0 g。试验前将自然晾干的样品装入长57 cm、宽37 cm、深5 cm的沙盘中，铺满、表面刮平，然后均匀喷洒相应水量，静置并晾干，晾干后的沙盘表面会形成一层均匀的物理结皮，但由于喷洒水量不同，形成的结皮厚度及硬度不同，然而结皮厚度及硬度不易测量且误差较大，因此，本文以喷洒的水量来量化结皮梯度。

1.2.2 风洞试验设计

风洞试验于2019年8月在中国科学院新疆生态与地理研究所环境风洞实验室(位于莫索湾沙漠研究站)完成。该风洞为直流开口吹气式，全长16.2 m，实验段长8 m，截面宽1.3 m，高1 m，采用侧壁扩散式结构，每侧壁扩散角0.2°，风速0 ~ 20 m/s可调。尾矿库年平均风速2 ~ 3 m/s，起沙风大多分布在6 ~ 16 m/s，因此试验设置风速梯度为8、10、12、14 m/s，每个梯度样品分别用8、10、12、14 m/s的风速吹蚀30、15、7、3 min。

将制备好的沙盘放入风洞实验段中段中轴线上，沙盘下风向5 cm处放置2个梯度集沙仪，用以测量风蚀颗粒分布；1个美国Sensit H14-LIN型风蚀传感器，用以测量风蚀颗粒碰撞数和碰撞能量；1个美国014A型风速传感器，用以记录实时风速^[22-23]。梯度集沙仪高度20 cm，共20层，每层开口1 cm×1 cm。风蚀传感器和风速传感器同步贴地测量，数据采集和记录频率为1 s，采集器CR1000自动记录。试验布设如图 1所示。每个沙盘试验前后用精度1 g的电子秤称重，用以计算样品风蚀率。

2 结果与分析 2.1 风速对尾矿风蚀率的影响

根据各组试验沙盘吹蚀前后质量差，计算尾矿风蚀率，并与吹蚀风速进行幂函数(y=ax^b)、指数函数(y=ae^bx)、二次多项式(y=ax²+bx+c)和对数函数(y=a-bln^(x))拟合，得到不同结皮量与含水率尾矿风蚀率与风速的关系(图 2)及最佳拟合函数结果(表 2)。结果显示，风速对风蚀率的影响十分显著。风蚀率与风速的最佳拟合函数为二次多项式，其在0.01显著性水平上的相关系数均在0.98以上。对照组风蚀率随风速增加由4.03 g/(m²·min) 增加到8 171.17 g/(m²·min)，增加2 028.60倍；含水率1.74% 组风蚀率随风速增加201.29倍，结皮量25 g组增加4 245.25倍；含水率2.72% 组增加52.49倍，结皮量50 g组增加145.67倍；含水率3.57% 组增加10.33倍，结皮量100 g组增加89.00倍。无论何种处理，风蚀率均随着风速的增加而增加，但随着含水率和结皮量的增加，风蚀率随风速的增长率在降低。目前，部分学者^{[2, 12-13]}给出的草地、沙地、农田风蚀量随风速变化遵循幂函数分布，而本研究与部分学者^{[19, 22-23]}给出的工业弃渣、煤矿土风蚀量随风速分布结果一致，说明土壤差异对风蚀量随风速变化规律影响很大，主要体现在土壤的理化性质、抗风蚀性和表层粗糙度等方面。

2.2 结皮量与含水率对尾矿风蚀率的影响

通过对各试验组风蚀率与结皮量、含水率进行拟合，得到风蚀率与结皮量、含水率的关系(图 3)，以及最佳拟合结果(表 3)。结果发现，结皮量和含水率变化对风蚀率的影响很大。风蚀率与结皮量、含水率的最佳拟合函数均为二次多项式，除8 m/s风速下结皮量的相关系数较低外，其余相关系数均在0.93以上。8 m/s风速下各处理和对照的风蚀率均较低，不超过30 g/(m²·min)，含水率和结皮量增加对尾矿风蚀率的影响不明显。10 m/s风速下含水率和结皮量增加对风蚀率的影响逐渐显现，随着含水率的提高，尾矿风蚀率降低了765.58 g/(m²·min)；随着结皮量的提高，尾矿风蚀率降低了760.43 g/(m²·min)。12 m/s风速下，随着含水率和结皮量的提高，尾矿风蚀率分别降低了3 336.98和3 317.65 g/(m²·min)。14 m/s情况下，随着含水率和结皮量的提高，尾矿风蚀率分别降低了8 116.41和8 069.71 g/(m²·min)。除8 m/s外，其他风速条件下风蚀率随结皮量和含水率的增加而降低，且风速越大，结皮量与含水率增加对风蚀率的影响越大。风蚀影响因子主要与风速、土壤抗蚀性和地表粗糙度等有关，结皮量和含水率对尾矿风蚀率产生的影响显然是由于结皮量增加了地表粗糙度，增强了表层抗蚀能力；含水率则增加了土壤颗粒粘结性，增强了土壤抗蚀性。

2.3 结皮量与含水率对尾矿风蚀颗粒分布的影响

本文利用梯度集沙仪收集各试验组样品吹蚀出的颗粒物，得到各结皮量和含水率条件下的风蚀颗粒随高度分布情况(图 4)及总风蚀颗粒量差异(图 5)。结果显示，总风蚀颗粒量随着风速的增大而增加，随着含水率和结皮量的增加而减小。各处理风蚀颗粒基本分布在10 cm内，部分处理在3 cm高度处达到最大，各处理3 cm以下颗粒物占总风蚀颗粒物的50% 以上，最大能到76.92%。风速8 m/s时，各处理颗粒分布随高度变化不明显，介于0.011 ~ 0.039 g/min。随着风速的增大，颗粒分布结构逐渐出现差异。10 m/s时，对照组，含水率1.74% 和2.72% 处理下的颗粒分布随高度先增加后减小；12 m/s时，对照组，含水率1.74%、2.72%，以及结皮量25 g处理下的颗粒分布结构先增加后减小；14 m/s时，对照组，含水率1.74%、2.72%，以及结皮量25和50 g处理下的颗粒分布结构先增加后减小。随着风速的增大、含水率的降低，结皮量的减少，这种先增加后减小的分布结构越容易出现。出现这种现象的原因是当风速增大时，表层颗粒受到的剪切力增大，颗粒间碰撞能量增大，而含水率和结皮量保护不足，导致更多颗粒发生跃移，跃移高度增加，从而增大了风蚀颗粒物释放，改变了颗粒分布结构。

2.4 结皮量与含水率对风蚀颗粒碰撞的影响 2.4.1 结皮量与含水率对颗粒碰撞数及能量的影响

各试验组在吹蚀时间内的颗粒碰撞总数及能量如图 6所示，颗粒碰撞能量与数量同步变化。风速、含水率、结皮量变化均会对颗粒碰撞产生影响。风速越大，颗粒碰撞数及能量越大；除对照组外，含水率和结皮量越高，颗粒碰撞数及能量越小；与风速、含水率和结皮量对尾矿风蚀率、风蚀颗粒分布的规律一致。风速8 m/s试验组由于风速较小，结皮量100 g及含水率3.57% 试验组由于抗风蚀能力较强，风蚀颗粒释放较少，颗粒碰撞数和能量也较小。

2.4.2 颗粒碰撞数随时间的变化

由于14 m/s风速下风蚀效果最显著，颗粒碰撞数及能量最大，且颗粒碰撞数和能量同步变化，因此本文选择14 m/s风速下不同含水率及结皮量吹蚀试验进行颗粒碰撞数随时间变化规律的分析，结果如图 7所示。不同处理下的尾矿颗粒释放随时间变化呈现相同的规律。除结皮量25 g处理外，其他处理在风速增加段(约前25 s)，颗粒碰撞数迅速增大到最大值，当风速到达14 m/s并稳定后，颗粒碰撞数随时间逐渐降低，且随着含水率和结皮量的增加，颗粒碰撞数峰值逐渐减小，每秒颗粒碰撞数也逐渐减小。含水率3.57% 的峰值颗粒碰撞数仅有2个，总颗粒碰撞数仅有65个。在风速上升过程中，不同处理首次出现颗粒碰撞的风速也有所差异，对照组首次出现颗粒碰撞的风速为8.45 m/s，随着含水率和结皮量的增加，首次出现颗粒碰撞的风速逐渐增大到10.05 m/s，说明随着含水率和结皮量的增加，尾矿颗粒释放所需风速越大。

3 讨论

试验结果表明，风速、含水率和结皮量对尾矿风蚀和颗粒释放均有较大影响。风速增加对风蚀和颗粒释放有促进作用，含水率和结皮增加对风蚀和颗粒释放有抑制作用，二者耦合作用下，含水率和结皮量增加能有效抵消风速增加对风蚀和颗粒释放的影响，且风速越大效果越明显。但这仅仅从整体上说明了风速与含水率、结皮量对尾矿风蚀和颗粒释放的影响，没有体现含水率和结皮量之间的差异。图 8比较了含水率和结皮量在不同风速下的抗风蚀效率，高含水率和结皮量(3.57% 和100 g)在任何风速下防护效果均能达到98% 以上，差异不大；中含水率和结皮量(2.72% 和50 g)在任何风速下，结皮的防护效果均比含水率好，且呈现出随着风速的增大，含水率防护效果增强，结皮效果略微降低的趋势；低含水率和结皮量(1.74% 和25 g)，这种趋势更加明显，低、中风速(10、12 m/s)结皮的防护效果比含水率好，但高风速(14 m/s)含水率的防护效果比结皮好。造成同一结皮量防护效果随风速增大而降低这一现象是由于结皮量与风速耦合作用有临界值，当结皮量小于或风速大于临界值，结皮会被破坏，一旦结皮被吹破，破口处防护效力基本为零，因此，25 g和50 g结皮量下，防护效力随风速增大而减小，25 g结皮量能有效防护的临界风速为12 m/s，50 g以上结皮量能有效防护的风速大于14 m/s。同时，造成相同含水率防护效果随风速增加而增加这种现象可能跟吹蚀时间有关，10 m/s吹蚀15 min，12 m/s吹蚀7 min，14 m/s吹蚀3 min，新疆夏季高温干燥，样品暴露于空气时间越长，水分散失越多，抗风蚀效率越低，且含水率越低对吹蚀时长越敏感，当含水率增大后，吹蚀时长造成的水分散失对抗风蚀效率的影响逐渐减弱。

综上，在结皮量足够的情况下，结皮的防护效果、适应性和稳定性比含水率高。一般来说，为防止尾矿风蚀，可采用物理法、植物法和化学法进行防护。物理法主要是采用物理方法覆盖尾矿表面，比如覆盖砾石、树皮、稻草，向尾矿喷水形成固化层等；植物法主要是在尾矿上栽种永久性植物；化学法主要是利用可与尾矿发生反应的化学试剂，比如石灰、硅酸钠等，在尾矿表面形成坚固硬壳^[24-26]。相比较而言，植物法持续效果好、环境友好度高，但植物筛选难、持续投入多、管护难；化学法和物理法类似，但化学法对环境污染大、成本高、化学试剂选择难；物理法是较为适合西北干旱区的一种低成本防护方式，且本文验证了物理结皮法的可行性及防护效果。对于本研究的尾矿，200 mm的全年降水量，相当于在试验沙盘上喷洒了42 g的水，能承受14 m/s的风吹蚀，已可以发挥较好的防护效果，但这是多次降雨的总降水量，单次降水量以及两次降水之间结皮的老化问题未考虑，因此，结皮野外实际的防护效果还需更多试验加以研究。

4 结论

风速、含水率和结皮量对尾矿风蚀影响显著。风蚀率随风速的增加而增加，随结皮和含水率的增加而减小，均符合二次多项式函数拟合，二者耦合作用下，含水率和结皮量增加能有效抵消风速增加对风蚀的影响，且风速越大效果越明显。

风速、结皮量和含水率对尾矿风蚀颗粒释放也会产生一定影响。尾矿风蚀颗粒主要分布在3 cm以下，占总风蚀颗粒的50% 以上。颗粒分布结构会随着风速的增加、含水率和结皮量的降低，逐渐出现先增加后减小的现象。总风蚀颗粒物、风蚀颗粒物碰撞数量和能量随着风速的增加而增加，随着含水率和结皮量的增加而减小。首次出现颗粒碰撞的风速随着含水率和结皮量的增加而增加。颗粒碰撞数随时间变化呈现一定的规律，在风速增加段，颗粒碰撞数迅速增大到最大值，当到达最大风速并稳定后，颗粒碰撞数随时间逐渐降低，且随着含水率和结皮量的增加，颗粒碰撞数峰值逐渐减小，每秒颗粒碰撞数也逐渐减小。

在结皮不被吹破的情况下，结皮的抗风蚀效果要优于含水率，因此，结皮抗风蚀效果与风速和结皮量有关，10 m/s以下风速任何结皮量均能有效抗风蚀，12 m/s需25 g以上结皮量，14 m/s需要50 g以上结皮量。而含水率在任何试验风速下需达到2.72% 以上才有较好的抗风蚀效果。