渗流—应力耦合作用下尾矿库回采的稳定性分析研究

Stability Analysis of Tailings Ponds Reclamation Stoping Under the Coupled Stress and Seepage Fields

摘要

尾矿库虽然是矿山重大危险源与潜在环境污染源，但具有巨大开发潜力。推进尾矿资源安全回采及二次利用不仅能够解决尾矿库的安全与环境问题，而且对矿山可持续发展具有重大意义。基于渗流—应力耦合理论，运用COMSOL多物理场仿真软件建立尾矿库坝体二维模型。通过理论分析与数值模拟相结合的方法，系统研究了尾矿库回采过程中坝体的渗流规律、应力变化、变形特征及稳定性演化规律。结果表明：在渗流—应力耦合作用下，随着尾矿库回采的推进，尾矿库的干滩长度减小，饱和区域逐渐扩大，孔隙水压力、von Mises应力和最大剪应力均逐渐增大，导致尾矿坝的整体位移和等效塑性应变逐渐增加，安全系数逐步减小，尾矿库的稳定性呈下降趋势。尽管如此，尾矿库回采过程中的安全系数仍符合相关规范要求，在安全生产管理条件下，坝体能够安全稳定运行。该研究成果为尾矿库的安全回采提供了重要理论依据和工程技术参考。

Abstract

While tailings ponds are notable sources of environmental pollution in mining operations, they also offer considerable potential for development. Advancing the safe mining and secondary utilization of tailings resources not only mitigates the safety and environmental challenges associated with tailings ponds but also significantly contributes to the sustainable development of mining activities. To examine the impact of seepage-stress coupling on the stability of tailings ponds during reclamation, the interaction between the seepage field and the stress field was analyzed, leading to the construction of a seepage-stress coupling analysis model for tailings ponds. Utilizing COMSOL multiphysics simulation software, a two-dimensional model of the tailings dam was developed. By integrating the seepage-stress coupling theory with the strength reduction method, the study systematically investigated the seepage patterns, stress variations, deformation characteristics, and stability evolution of the tailings dam throughout the stoping process. Through numerical simulations, variations in pore water pressure, von Mises stress, shear stress, overall displacement, equivalent plastic strain, and safety factor of the tailings dam were evaluated across six stoping scenarios: 0, 50, 100, 150, 200, 230 m. The findings indicate that under the influence of seepage-stress coupling, the reclamation of the tailings pond results in a reduction in the length of the dry beach, while pore water pressure, von Mises stress, and maximum shear stress increase. Consequently, there is a gradual rise in the overall displacement and equivalent plastic strain of the tailings dam, accompanied by a decrease in the safety factor and a declining trend in the stability of the tailings pond. Nonetheless, the safety factor throughout the reclamation process remains within regulatory standards, ensuring that, under conditions of safe production management, the dam can function normally. These findings provide a significant theoretical foundation and engineering technical reference for the safe reclamation of tailings ponds.

Graphical abstract

关键词

尾矿资源 / 尾矿库回采 / 渗流—应力耦合 / 尾矿坝 / 稳定性分析 / 数值模拟

Key words

tailings resources / tailings ponds reclamation / seepage-stress coupling / tailings dams / stability analysis / numerical simulation

尾矿库是用于堆存尾砂或其他工业废渣的场所，具有安全风险高、环境污染风险大和管理难度大等特点，是高势能人造泥石流重大危险源与潜在环境污染源（王昆等，2024）。我国现有各类大小尾矿库1万多座，一旦发生尾矿库溃坝和尾矿泄漏等事故，不仅会严重威胁人民群众的生命财产安全，还会对生态环境造成极大的破坏，甚至可能会影响社会安定（张家荣等，2019）。因此，尾矿坝稳定性分析一直是尾矿库安全领域的重要课题。与此同时，推进尾矿资源的安全回采与二次综合利用，不仅能从根源上化解尾矿库的安全隐患与环境问题，还能实现资源循环利用，具有巨大的发展潜力。

尾矿坝稳定性是矿山企业十分关注的问题，国内外学者从静力稳定性、动力稳定性及渗流稳定性等多个方面对其开展研究（杨春和等，2021；许云美等，2023；朱远乐等，2024）。静力稳定性主要研究在尾矿坝自重、堆积荷载等静力条件下，尾矿材料的物理力学性质、坝体结构特征等因素对尾矿坝稳定性的影响（杨春和等，2022）；动力稳定性则研究尾矿坝在面临地震和爆破振动等动力荷载冲击时的稳定性特征（王文松等，2018；李生伟等，2023）。其中，尾矿坝渗流稳定性一直是研究的一个难点，学者们通过室内渗流试验、现场监测以及数值模拟等方法对尾矿坝内渗流场的分布进行分析，研究渗流对坝体稳定性的作用机制，为制定合理的排水减压措施提供科学依据（吴高杰等，2022；Du et al，2024；梁冰等，2025）。在尾矿坝渗流稳定性研究基础上，许多学者开始考虑多场耦合作用对尾矿坝稳定性的影响，并利用数值模拟软件建立尾矿坝多场耦合数值模型，研究尾矿坝在渗流—应力耦合（邓红卫等，2016；Popescu et al，2024）、流—固—化耦合（樊嘉庆，2020）以及热—流—固耦合（李向阳等，2022）作用下的稳定性，分析库水位、化学腐蚀、降雨强度和温度等因素对尾矿坝稳定性的影响规律，为尾矿坝的稳定性分析及后期管理提供理论支持。

近年来，尾矿资源利用受到国内外学者的广泛关注。以往研究表明，尾矿资源潜力巨大，通过推进尾矿综合利用，不仅能够消除尾矿库这一重大危险源，还能够有效促进矿区生态环境保护，对于实现资源循环利用和矿山可持续发展具有重要的战略意义（Araya et al，2021；李鑫等，2024；伍天华等，2024）。随着选矿技术装备的不断进步，国内外学者对尾矿资源安全回采及二次利用进行了研究（Alcalde et al，2018；任珩等，2024）。例如：通过对尾矿库回采方法、回采顺序及回采工艺的系统研究，对比分析了不同回采方案的优缺点，并提出了一种改进的尾矿库回采方案（王新岩，2023；付宪吕等，2024；郭大平等，2023）。在此基础上，通过室内试验（周罕等，2023）和数值模拟（段晓阳等，2023；Pan et al，2023；肖容等，2023）等方法，研究了在尾矿库回采的不同阶段尾矿坝稳定性特征的演化规律，确保尾矿坝在整个回采过程中保持稳定运行。

综上所述，国外学者关于尾矿库资源回采的研究主要集中在回采技术和工艺方法等方面，国内学者则更多关注尾矿库回采过程中的静力稳定性分析，以及不同回采方法对尾矿坝稳定性的影响。然而，在综合考虑渗流—应力耦合作用对尾矿库回采过程中稳定性影响的研究方面，目前国内外学者开展的研究工作相对较少，缺乏系统性成果。在尾矿库回采过程中，干滩长度与浸润线的变化会引发渗流场的改变，从而破坏尾矿库原有的应力平衡状态，导致尾矿坝产生变形破坏。因此，开展尾矿库回采过程中渗流—应力耦合作用下的稳定性研究，对于揭示坝体变形破坏机制和保障尾矿库安全回采具有重要的理论意义和工程参考价值。

1 渗流—应力耦合分析模型研究

在尾矿库回采过程中，尾砂介质内部水头分布发生变化，地下水在水头差驱动下沿着孔隙发生渗流运动，产生渗流压力。渗流压力以体积力的形式作用于尾砂介质，改变其应力场，进而引起位移场的变化。位移场的改变反过来又影响介质的渗透系数和渗流场分布，形成渗流—应力耦合作用。

1.1 渗流场对应力场的作用

渗流场对应力场的作用主要表现为流体在孔隙介质中运动时对固体骨架施加的孔隙水压力，包括静水压力和超静水压力（渗流压力）。

一方面，尾矿库原有地下水及回采工作面基坑中的积水会给尾矿库施加静水压力，它以压力荷载的形式作用于接触面上，其函数表达式为

P={\gamma }_{\mathrm{w}}(H-z)（1）

式中：P为静水压力；{\gamma }_{\mathrm{w}}为水的容重；H为位置水头；z为水头高度。

另一方面，尾砂介质中的流体在水头差的作用下发生渗流运动，产生渗流压力，以渗流体积力的形式作用于尾砂介质中，从而影响尾矿库内原有的应力场和位移场，改变其原有的应力平衡状态。根据水力学理论，渗流体积力与水力梯度成正比（李雪梅，2017），其函数表达式为

\left\{\begin{array}{c}{f}_x={\gamma }_{\mathrm{w}}{J}_x=-{\gamma }_{\mathrm{w}}\frac{\partial H}{\partial x}\\ f_y={\gamma }_{\mathrm{w}}{J}_y=-{\gamma }_{\mathrm{w}}\frac{\partial H}{\partial y}\end{array}\right.（2）

f=\sqrt[]{{f_x}^{\mathrm{2}}+{f_y}^{\mathrm{2}}}（3）

式中：f为渗流体积力；f_x、f_y分别为渗流体积力沿x、y方向的分力；J_x、J_y分别为水力坡降沿x、y方向的分量。

1.2 应力场对渗流场的作用

随着尾矿库表面尾砂的回采，其内部应力将重新分布，导致某些部分卸荷或荷载减少，形成卸荷效应（如应力释放和体积膨胀），从而使尾矿库原有的应力场和位移场发生变化。应力场和位移场的变化会改变尾砂介质的孔隙结构和分布特征，多孔介质骨架的体积应变会导致孔隙率等孔隙特征发生改变，进而影响尾砂介质的渗透系数，最终导致渗流场的变化。因此，应力场对渗流场的影响主要是通过改变尾砂介质的孔隙特征和渗透系数来实现的。多孔介质的渗透系数可用应力状态函数表示为

K=K({\sigma }_{ij})（4）

此时，受应力场影响的二维渗流微分方程可表示为

\frac{\partial }{\partial x}[K({\sigma }_{ij})]\frac{\partial H}{\partial x}+\frac{\partial }{\partial y}[K({\sigma }_{ij})]\frac{\partial H}{\partial y}=\mathrm{0}（5）

式中：{\sigma }_{ij}为应力张量；K为渗透系数。

1.3 渗流—应力耦合模型在COMSOL中的实现

COMSOL多物理场仿真软件可通过添加多个物理场接口实现多物理场耦合分析，具有易操作和可视化强的特点，被广泛应用于渗流—应力耦合分析模型研究中。本文研究中假设地下水服从达西定律，尾矿库尾砂为各向同性等效连续介质，因此渗流场选择达西定律接口，应力场选择固体力学接口，以此搭建尾矿库渗流—应力耦合分析模型。

在COMSOL软件进行渗流—应力耦合计算中，考虑孔隙水压力作用的达西定律表达式为

V=-\frac{k}{\mu }(\nabla p+\rho \mathrm{g})（6）

式中：V为达西速度；k为渗透率；\mu为流体的黏度；p为孔隙水压力；\rho为流体密度；\mathrm{g}为重力加速度。

固体力学和达西定律接口进行耦合计算时，通过达西定律接口计算出模型中的渗流场分布情况，并将其作为初始条件导入到固体力学接口中，以此实现渗流场与应力场的耦合分析。由于尾矿库浸润线以上部分尾砂处于非饱和状态，因此用理查兹方程来描述流体在尾砂介质中的非饱和渗流过程，考虑有效应力的控制方程为

\rho (\frac{C_{\mathrm{m}}}{{\rho }_g}+S_{\mathrm{e}}S)\frac{\partial p}{\partial t}+\nabla \bullet \rho [-\frac{k_{\mathrm{s}}}{\mu }{k}_{\mathrm{r}}(\nabla p+{\rho }_{\mathrm{g}}\nabla D)]=Q_{\mathrm{m}}-\rho \alpha \frac{\partial \epsilon }{\partial t}（7）

式中：C_{\mathrm{m}}为含水量；S_{\mathrm{e}}为有效饱和度；S为储水系数；k_{\mathrm{s}}为介质饱和渗透率；k_{\mathrm{r}}为相对渗透率；Q_{\mathrm{m}}为渗流源汇项；\alpha为Biot系数；\epsilon为体积应变量。

2 尾矿库数值模型的构建

2.1 工程概况

河南省某钼矿尾矿库总体上为山谷型尾矿库，局部属傍山型。该尾矿库目前堆积坝顶标高约为1 420.76 m，现状总坝高为60.50 m，全库容约为85×10⁴ m³，有效库容约为76.5×10⁴ m³。初期坝为碾压式土石坝，坝高为14.76 m，坝顶长度为69.6 m，顶宽为3.17 m，外坡比为1∶2.48，内坡比为1∶1.75。堆积坝采用上游式筑坝，堆积坝高约为45.74 m，水平长度为260.62 m，堆积坝整体外坡比约为1∶5.62，共堆有18级子坝。子坝坝面覆土植草，植被生长茂密，护坡良好。该尾矿库现状等别为四等库，根据《尾矿库安全规程》（GB39496-2020），正常运行时尾矿坝稳定性计算的荷载包括运行期正常库水位时的稳定渗透压力、坝体自重以及坝体和坝基中的孔隙水压力，坝体正常运行的最小安全系数为1.25。

本次设计回采的尾砂按照成分可划分为2个部分进行利用，通过取样化验，具备再选价值的尾砂输送至选矿厂，通过二次造浆后进行再选；不具备再选价值的尾砂则用作建筑材料，用于建设或修整矿区内道路。

2.2 尾矿库回采方式

尾矿库采用坝前干式与库内湿式的联合回采方式，按照先内后外、先库后坝、先上后下和分层开采的顺序进行回采，如图1所示。库内湿式回采区域选择水枪与漂浮式砂泵船联合工艺进行尾砂回采，由库内向坝前方向推进，设计回采分层高度为2 m，回采坑坡比按1∶3放坡以防止坡面垮塌，回采采坑水面以上的边坡高度应小于等于3 m。坝前干滩部分选择干式回采，采用推土机将干滩区堆积的尾砂推送至湿采区再回采，在回采过程中需保持滩面坡度向库内倾斜大于等于1%，干滩长度大于等于50 m。在回采初始阶段，库内湿式回采区域先开采1~2层，形成必要滩面形态，然后按库内湿采区与坝前干采区有序循环进行回采，回采过程中干采区域始终高于湿采区域2 m以上。

2.3 模型建立及参数设定

结合矿山企业提供的工程地质剖面图及尾砂沉积规律，建立尾矿库坝体稳定性二维分析模型，如图2所示。构建模型时，固体力学模块选择德鲁克—普拉格本构模型匹配莫尔—库伦准则进行计算，对模型底部施加固定约束，左右两侧施加辊支撑约束，同时对整个模型施加y方向的重力。渗流场则选择达西定律接口开展计算研究，材料选择非饱和多孔介质，持水模型选择V-G模型，根据工程地质剖面图中浸润线的位置确定压力水头边界，模型底部为无流动边界。尾矿库坝体材料物理力学参数根据室内试验成果和现有的工勘报告等相关资料综合确定，计算所需材料的物理力学参数见表1。V-G模型计算需要用到模型参数α、n、l和残余含水率θr等参数，相关参数根据室内试验并参考经验数值选取，如表2所示。在渗流—应力耦合计算时，首先通过达西定律模块计算出尾矿库内部孔隙水压力分布，然后将渗流场计算得到的孔隙水压力作为荷载施加到应力场。通过固体力学和达西定律接口进行耦合计算，基于有效应力原理和强度折减原理分析尾矿库在回采过程中的渗流场、应力场及稳定性变化情况。

3 尾矿库回采过程渗流与应力耦合分析

根据尾矿库回采方案，尾矿库回采过程中干滩长度应大于等于50 m，因此取回采过程推进0，50，100，150，200，230 m共6种工况，利用COMSOL软件分析回采过程中尾矿库坝体的渗流规律、应力变化、变形特征和坝体稳定性。在尾矿库回采过程中，不同推进距离工况下尾矿库坝体的最大孔隙水压力、最大Mises等效应力、最大剪应力、最大位移、最大等效塑性应变和安全系数见表3。由于尾矿库回采尚处于可行性研究阶段，目前矿山企业已根据相关研究成果开展尾矿库回采先期准备工作，后续将采用尾矿库回采过程中的现场监测数据与模拟结果进行对比，验证数值模拟的可靠性。

3.1 尾矿库回采过程渗流规律分析

在尾矿库回采过程中，随着回采区域由库内向坝前推进，基坑中的水将向尾砂介质中渗流，导致尾矿库的渗流场发生变化。图3所示为尾矿库回采过程中尾矿库坝体的饱和度分布情况。图中红色区域为饱和区域，蓝色区域为非饱和区域。由图3可以看出，随着回采区域向前推进，尾矿库原本非饱和区域变为饱和区域，饱和区域范围逐渐增大，且有向堆积坝表层蔓延的趋势；非饱和区域范围逐渐减小，主要集中在堆积坝表层和初期坝区域。随着尾矿库从0 m回采至230 m，尾矿库干滩长度由280 m减小至50 m，饱和区域面积由10 811 m²增大至12 455 m²（不考虑基岩），增幅为15.21%。图4所示为尾矿库回采推进0 m和230 m这2种工况下尾矿库坝体的孔隙水压力分布云图。为了查看尾矿库回采过程中坝体应力情况，对模型中的基岩部分进行隐藏，该操作不影响相关数值计算结果。由图4可知，尾矿库坝体的孔隙水压力随着地层深度的增加而增大，最大孔隙水压力位于尾粉土层底部与基岩结合处。

这是由于尾粉土的渗透性较低，阻碍了地下水的快速排出，导致孔隙水压力在该区域积聚。图5中灰色曲线为尾矿库回采过程中尾矿库坝体的最大孔隙水压力变化情况。随着尾矿库从0 m回采至230 m，尾矿库坝体的最大孔隙水压力逐渐增大，压力值从225.402 kPa增大至277.866 kPa，增幅为23.28%。孔隙水压力增大的原因主要是随着回采工作的推进，基坑中的水增多并向深部渗流，导致下部尾砂中的水无法及时排出，从而形成超静水压力。孔隙水压力的增大会降低尾砂介质的有效应力，导致土体的抗剪强度降低，甚至引发液化现象，进而增加尾矿库变形和破坏的风险，严重影响尾矿坝的整体稳定性。

3.2 尾矿库回采过程应力变化分析

von Mises应力是一种基于剪切应变能的等效应力，可用于快速确定尾矿坝的最危险区域。图6所示为尾矿库回采推进0 m和230 m这2种工况下坝体的von Mises应力分布云图。由图6可以看出，von Mises应力随土层深度的增加而增大，其最大值出现在初期坝内侧及堆积坝最大埋深处。这表明该区域为应力集中区，是尾矿坝的最危险区域。由于初期坝的强度远高于堆积坝，因此堆积坝最大埋深处更容易发生变形和破坏。图5中红色曲线为尾矿坝的最大von Mises应力变化曲线。由图5可以看出，尾矿坝的最大von Mises应力随回采推进距离的增加而增大。其中，0~50 m和200~230 m回采段的斜率较大，主要原因是：回采0~50 m区域时，回采工作面中的积水开始渗流到尾砂介质中，尾矿库原本非饱和区的饱和度提高，孔隙水压力增大，导致von Mises应力显著增大；回采200~230 m区域时，孔隙水压力的进一步增大导致尾砂介质的有效应力降低，从而加速了von Mises应力的增大。von Mises应力的增大会降低尾矿坝的抗剪强度，增大变形和破坏的风险，进而影响其整体稳定性。

Tresca剪应力则主要用来表征尾矿坝的最大剪应力，从而判断其发生剪切破坏的最危险区域。图7所示为尾矿库回采推进0 m和230 m这2种工况下坝体的Tresca剪应力分布云图。与von Mises应力分布特征类似，Tresca剪应力也随着土层深度的增加而增大，最大值出现在初期坝内侧及堆积坝最大埋深处。图5中蓝色曲线为尾矿坝的最大Tresca剪应力变化曲线。由图5可知，尾矿坝的最大Tresca剪应力随着回采推进距离的增加而增大，其变化趋势与von Mises应力基本一致。其中，0~50 m和200~230 m回采段的斜率较大，主要原因是：回采0~50 m时，回采工作面中的积水导致孔隙水压力增大，土体的饱和度提高，剪应力显著增大；回采200~230 m区域时，浸润线的上升进一步加剧了孔隙水压力的增大，导致剪应力加速上升。Tresca剪应力的增大会增加尾矿坝发生剪切破坏的风险，特别是在堆积坝最大埋深处，由于土体强度较低，更容易发生破坏。

3.3 尾矿库回采过程变形特征分析

在实际工程中，尾矿坝的位移是判断其稳定性的重要指标。图8所示为尾矿库回采推进0 m和230 m这2种工况下坝体的位移分布云图。由图8可以看出，尾矿坝最大位移出现在堆积坝底部与初期坝交界处。这是由于该区域的孔隙水压力、von Mises应力和Tresca剪应力均较大，对土体的变形和破坏产生了显著影响。随着回采区域向前推进，尾矿坝最大位移所在区域从堆积坝底部向堆积坝顶部扩展。图9中灰色曲线展示了回采过程中尾矿坝的最大位移变化情况。由图9可以看出，尾矿坝的最大位移随着回采推进距离的增加而增大，当回采区域从0 m推进至230 m时，尾矿坝的最大位移从123.960 mm增大至153.616 mm，增幅为23.92%。位移的显著增大会增加尾矿坝发生失稳破坏的风险。在尾矿库回采过程中，由于渗流场与应力场的耦合作用，坝体的潜在滑动面表现为典型的圆弧形滑动面，并有向坝体顶部蔓延的趋势。这种滑动面的形成主要是由于坝体内部发生剪切破坏，特别是在堆积坝底部与初期坝交界处，由于应力集中和孔隙水压力增大，土体的抗剪强度显著降低，因此更容易发生滑动破坏。

3.4 尾矿库回采过程坝体稳定性分析

图10所示为尾矿库回采0 m和230 m这2种工况下坝体的等效塑性应变云图。由图10可以看出，塑性区主要位于尾矿库堆积坝与基岩及初期坝结合处，塑性应变最大值出现在尾矿库堆积坝底部与初期坝顶部交界处。这是由于该区域的应力集中和孔隙水压力增大，导致土体的抗剪强度显著降低，从而容易发生塑性变形。图9中红色曲线展示了尾矿库回采过程中坝体的最大等效塑性应变变化情况。由图9可以看出，尾矿库坝体的最大等效塑性应变随着回采推进距离的增加而增大。当回采区域从0 m推进至230 m时，最大等效塑性应变由0.0407增大至0.0500，增幅为22.85%。最大等效塑性应变的变化趋势与最大位移的变化趋势一致，表明塑性变形是导致坝体位移增大的主要原因。

运用强度折减理论，借助COMSOL软件的参数辅助扫描功能，计算得到尾矿库回采过程中尾矿库坝体的安全系数，如表3所示。图9中蓝色曲线为尾矿库回采过程中尾矿库坝体的安全系数变化情况。由表3和图9可知，随着回采区域由0 m推进至230 m，尾矿库坝体的安全系数从1.881下降至1.743，降幅为7.34%。安全系数的下降主要是由于回采过程中干滩长度减小，导致孔隙水压力增大和潜在滑动面的摩擦系数降低。此外，安全系数的下降趋势在回采后期呈现加速现象。这是由于随着回采区域的推进，尾矿库湿采区域的水渗流到尾砂介质中，导致坝体内部侵蚀加剧和浸润线抬高，有效应力进一步降低，进而引起安全系数加速下降。根据数值模拟结果，尾矿库回采过程中的安全系数均超过《尾矿库安全规程》（国家市场监督管理总局，2020）规定的最小值（1.25），且具有一定富余。这说明尾矿坝在回采过程中处于稳定状态，研究结果可为尾矿库回采工程提供重要参考。

4 结论

以河南省某钼矿尾矿库作为研究对象，基于渗流—应力耦合理论，运用COMSOL软件建立了尾矿库的渗流—应力耦合分析模型，系统研究了尾矿库回采过程中坝体的渗流规律、应力变化、变形特征及稳定性演化规律，得出如下结论：

（1）在渗流—应力耦合作用下，随着尾矿库回采向前推进，尾矿库的干滩长度减小，饱和区域逐渐扩大，孔隙水压力、von Mises应力和最大剪应力均逐渐增大，导致尾矿坝整体位移和等效塑性应变逐渐增大，尾矿坝的安全系数逐渐减小，尾矿库的稳定性下降。

（2）在尾矿库回采过程中，尾矿坝的安全系数由1.881下降至1.743，但仍远大于《尾矿库安全规程》（国家市场监督管理总局，2020）规定的四等坝安全系数（1.25），说明尾矿库回采过程中的稳定性满足规范要求，在满足安全生产管理条件下，坝体可正常运行。本文研究成果能够为尾矿库回采稳定性及可行性分析提供理论支持。

参考文献

[1] Alcalde J， Kelm U， Vergara D，2018.Historical assessment of metal recovery potential from old mine tailings：a study case for porphyry copper tailings，Chile［J］. Minerals Engineering，127：334-338.

[2] Araya N， Mamani Quiñonez O， Cisternas L A，et al，2021.Sustainable development goals in mine tailings management：targets and indicators［J］. Materials Proceedings，5（1）：82.

[3] Du C B， Tao H， Yi F，2024.Seepage stability analysis of geogrid reinforced tailings dam［J］. Scientific Reports，14：1814.

[4] Pan Y W， Chen J P， Zuo X H，et al，2023.The stability of dams with different stoping elevations in the Tongling valley-type tailings impoundment：a case study in Yunnan China［J］. Minerals，13（11）：1365.

[5] Popescu F D， Radu S M， Andras A，et al，2024.Stability assessment of the dam of a tailings pond using computer modeling-case study：Coroiesti，Romania［J］. Applied Sciences，14（1）：268.

[6] 邓红卫，李爽，邓畯仁，2016.渗流—应力耦合作用下尾矿库稳定性的三维数值分析［J］.安全与环境学报，16（4）：133-138.

[7] Deng Hongwei， Li Shuang， Deng Junren，2016.3D-numerical simulation on the stability of tailings dam under the coupled stress and seepage fields［J］. Journal of Safety and Environment，16（4）：133-138.

[8] 段晓阳，李学锋，吴高善，等，2023.某尾矿库回采技术及坝体稳定性分析［J］.矿业研究与开发，43（1）：48-54.

[9] Duan Xiaoyang， Li Xuefeng， Wu Gaoshan，et al，2023. Mining technology and dam stability analysis of a tailings pond［J］.Mining Research and Development，43（1）：48-54.

[10] 樊嘉庆，2020.基于流—固—化多场耦合分析的尾矿坝长期稳定性研究［D］.衡阳：南华大学. Fan Jiaqing，2020.Study on long-term stability of tailings dam based on fluid-solid-chemical multi-field coupling analysis［D］.Hengyang：University of South China.

[11] 付宪吕，程立家，傅灿，2024.某尾矿回采工程实践［J］.有色金属设计，51（4）：5-9.

[12] Fu Xianlü， Cheng Lijia， Fu Can，2024.Practice of a tailings stoping project［J］. Nonferrous Metals Design，51（4）：5-9.

[13] 郭大平，胡良才，李玉雷，等，2023.一种干式尾矿库安全高效回采方法［J］.铀矿冶，42（1）：15-22，54.

[14] Guo Daping， Hu Liangcai， Li Yulei，et al，2023.A safe and efficient method of tailings recovery for dry tailings pond［J］.Uranium Mining and Metallurgy，42（1）：15-22，54.

[15] 国家市场监督管理总局，2020. 尾矿库安全规程：GB39496-2020 ［S］.北京：中国质检出版社.State Administration for Market Regulation，2020.Safety regulation for tailings pond：GB39496-2020 ［S］.Beijing：China Planning Press.

[16] 李生伟，陈文婷，朱海丽，等，2023.达尔乌拉铁矿尾矿库三维有限元静动力分析［J］.安全与环境工程，30（2）：136-141.

[17] Li Shengwei， Chen Wenting， Zhu Haili，et al，2023.Static and dynamic analysis of tailing reservoir in Daerwula iron mine based on three dimensional finite element［J］.Safety and Environmental Engineering，30（2）：136-141.

[18] 李向阳，陈艳，洪昌寿，等，2022.热—流—固耦合作用下铀尾矿库坝体稳定性数值模拟［J］.安全与环境学报，22（6）：3135-3142.

[19] Li Xiangyang， Chen Yan， Hong Changshou，et al，2022.Numerical simulation of dam stability of a uranium tailings pond under thermal-fluid-solid coupling［J］.Journal of Safety and Environment，22（6）：3135-3142.

[20] 李鑫，王锦荣，邓荣东，等，2024.我国铅锌尾矿综合利用研究现状及其开发利用意义［J］.地质与勘探，60（4）：724-734.

[21] Li Xin， Wang Jinrong， Deng Rongdong，et al，2024.Research status and significance of comprehensive utilization of lead-zinc tailings in China［J］.Geology and Exploration，60（4）：724-734.

[22] 李雪梅，2017.基于流固耦合作用的尾矿坝稳定性研究［D］.昆明：昆明理工大学. Li Xuemei，2017.Study on stability of tailings dam based on fluid-solid coupling［D］.Kunming：Kunming University of Science and Technology.

[23] 梁冰，曾勇，易富，等，2025.尾矿坝渗流稳定性研究进展［J］.长江科学院院报，42（3）：107-117.

[24] Liang Bing， Zeng Yong， Yi Fu，et al，2025.Advances in research of seepage stability of tailings dams［J］.Journal of Changjiang River Scientific Research Institute，42（3）：107-117.

[25] 任珩，李沛霖，李金潞，2024.基于黄金尾矿处理与综合利用研究的文献计量分析［J］.黄金科学技术，32（5）：939-948.

[26] Ren Heng， Li Peilin， Li Jinlu，2024. Bibliometric analysis of gold tailings treatment and comprehensive utilization research［J］.Gold Science and Technology，32（5）：939-948.

[27] 王昆，张峥， Hudson-Edwards K A，等，2024.尾矿的特性、灾害及其资源潜力［J］.金属矿山，53（8）：216-227.

[28] Wang Kun， Zhang Zheng， Hudson-Edwards K A，et al，2024. Pro-perties，disasters，and resource potential of mine waste tailings［J］. Metal Mine，53（8）：216-227.

[29] 王文松，尹光志，魏作安，等，2018.基于时程分析法的尾矿坝动力稳定性研究［J］.中国矿业大学学报，47（2）：271-279.

[30] Wang Wensong， Yin Guangzhi， Wei Zuo’an，et al，2018.Study of dynamic stability of tailings dam based on time-history analysis method［J］.Journal of China University of Mining and Technology，47（2）：271-279.

[31] 王新岩，2023.尾矿库回采方案研究：以北京某山谷尾矿库为例［J］.科技创新与应用，13（11）：9-12.

[32] Wang Xinyan，2023.Research on mining plan of tailings pond：a case study of a valley tailings pond in Beijing［J］. Technology Innovation and Application，13（11）：9-12.

[33] 吴高杰，任伟中，占美杰，等，2022.武阳高速某尾矿库坝体渗流稳定性及溃坝演进过程分析［J］.安全与环境工程，29（4）：74-84.

[34] Wu Gaojie， Ren Weizhong， Zhan Meijie，et al，2022. Analysis of seepage stability and dam break evolution treatment of tailings pond dam on Wuyang expressway［J］. Safety and Environmental Engineering，29（4）：74-84.

[35] 伍天华，高永涛，宋振，等，2024.稀土尾矿资源综合利用研究现状及展望［J］.中国稀土学报，42（4）：635-665.

[36] Wu Tianhua， Gao Yongtao， Song Zhen，et al，2024.Comprehensive utilization of rare earth tailings resources：a review［J］. Journal of the Chinese Society of Rare Earths，42（4）：635-665.

[37] 肖容，韩追，袁利伟，等，2023.尾矿库资源回采过程中坝体的稳定性分析［J］.有色金属（矿山部分），75（5）：11-16.

[38] Xiao Rong， Han Zhui， Yuan Liwei，et al，2023.Stability an-alysis of dam body during resource extraction of tailings pond［J］.Nonferrous Metals（Mining Section），75（5）：11-16.

[39] 许云美，袁利伟，龙皓楠，2023.干堆尾矿库稳定性影响因素的敏感性分析［J］.黄金科学技术，31（6）：1014-1022.

[40] Xu Yunmei， Yuan Liwei， Long Haonan，2023. Sensitivity an-alysis of stability influencing factors of dry heap tailings reservoir［J］.Gold Science and Technology，31（6）：1014-1022.

[41] 杨春和，张超，李全明，等，2021.大型高尾矿坝灾变机制与防控方法［J］.岩土力学，42（1）：1-17.

[42] Yang Chunhe， Zhang Chao， Li Quanming，et al，2021.Disaster mechanism and prevention methods of large-scale high tailing dam［J］. Rock and Soil Mechanics，42（1）：1-17.

[43] 杨春和，张超，马昌坤，等，2022.高应力条件下尾矿破碎特性及坝体稳定性研究［J］.中国安全生产科学技术，18（2）：20-26，2.

[44] Yang Chunhe， Zhang Chao， Ma Changkun，et al，2022.Study on tailings breakage characteristics and dam stability under high stress conditions［J］.Journal of Safety Science and Technology，18（2）：20-26，2.

[45] 张家荣，刘建林，2019.中国尾矿库溃坝与泄露事故统计及成因分析［J］.中国钼业，43（4）：10-14.

[46] Zhang Jiarong， Liu Jianlin，2019.The statistics and causes of dam break and leakage in Chinese tailings pond［J］.China Molybdenum Industry，43（4）：10-14.

[47] 周罕，付俊，陈永贵，等，2023.上游式尾矿库回采利用坝坡动态稳定性研究［J］.金属矿山，52（4）：199-202.

[48] Zhou Han， Fu Jun， Chen Yonggui，et al，2023.Research on variational stability of dam and slope for upstream tailings pond utilization［J］. Metal Mine，52（4）：199-202.

[49] 朱远乐，王淇萱，刘文岗，等，2024.尾矿库灾害风险防控研究现状及发展态势［J］.武汉大学学报（工学版），57（8）：1054-1064.

[50] Zhu Yuanle， Wang Qixuan， Liu Wengang，et al，2024.Research status and development trend of disaster risk prevention and control of tailings pond［J］.Engineering Journal of Wuhan University，57（8）：1054-1064.

基金资助

国家自然科学基金面上项目“寒区岩质散体冻胀裂解孕育排土场灾变机理及干预机制研究”(51874352)