基于环管试验的粗骨料膏体管输阻力模型及优化

杨晓炳 闫泽鹏 尹升华 杨航

摘 要：為探明粗骨料膏体的管输阻力特性，开展了关于粗骨料膏体输送特性的工业级 环管输送试验，利用白金汉方程对环管数据拟合处理后获得膏体的实际流变参数，并建立了基于流动度测试的粗骨料膏体阻力方程.通过响应面分析法（RSM-BBD）对膏体输送行为进行分析，分别得到单因素和多因素耦合对阻力损失的影响规律，并优化了粗骨料膏体的管道 输送参数.结果表明：流变测试得到的流变参数普遍大于环管拟合参数，其中屈服应力相差10.4%～12.2%、塑性黏度相差 21%～32.8%.通过响应面分析得到了单因素对阻力损失影响的敏 感程度为质量分数>管径>流量，并发现因素间的交互作用影响显著，其中质量分数和管径的交互影响起到决定性作用.以金川实际充填需求为背景采用响应面优化后的输送参数为：管 径 220mm、流量142.4 m3/h、质量分数73.2%.研究成果可为粗骨料膏体充填系统运行参数的选 取提供参考.

关键词：粗骨料膏体;管输阻力;环管试验;流变特性;响应面分析

中图分类号：TD862.2 文献标志码：A

Transportation Pressure Model and Optimization of Coarse Aggregate Paste Backfill Based on Pipe Loop Test

YANG Xiaobing1，2，3，YAN Zepeng1，3?，YIN Shenghua1，3，YANG Hang2

（1.School of Civil and Resources Engineering，University of Science and Technology Beijing，Beijing100083，China;

2.State Key Laboratory of Mineral Processing，Beijing102628，China;

3.Key Laboratory of High-Efficient Mining and Safety of Metal Mines of Ministry of Education（University of Science and Technology Beijing），Beijing100083，China）

Abstract：To determine the pipe pressure drop characteristics in the coarse aggregate paste backfill（CPB）pro-cess.An industrial-scale loop test on the transportation characteristics of CPB was carried out.The actual rheological parameters of CPB are obtained by fitting the Buckingham equation to the loop test data，and the resistance equationwas established based on the fluidity test.The response surface method（RSM-BBD）was used to analyze the trans-portation behavior of the CPB，and the effects of single-factor and multi-factor coupling on the pressure loss were ob-tained，respectively.Finally，RSM-BBD was used to optimize the pipeline transportation parameters of CPB.The re-sults show that the rheological parameters obtained by the rheological test were generally larger than those obtained by the loop test，where the differences in yield stress and plastic viscosity between the two are10.4%～12.2% and 21%～32.8%，respectively.Response surface analysis showed that the sensitivity of a single factor on pressure loss was as follows concentration > pipe diameter > flow rate.And it was found that the interaction between factors had a significant effect，in which the interaction between mass fraction and pipe diameter played a decisive role.Taking the actual filling demand of Jinchuan as the background，the transport parameters optimized by response surface were as follows：pipe diameter 220mm，flow rate142.4 m3/h，and concentration 73.2%.The research results are of great significance to the selection of the CPB system operating parameters.

Key words：coarse aggregate paste;pressure drop;loop test;rheological properties;response surface analysis

矿产资源开发利用过程中产生大量的废石、尾 砂等固体废物[1-2].据统计[3]，每产出1万 t矿石，平均 要产出0.93万 t尾矿、0.34万 t废石.矿山固体废料的排放不仅占用大量的田地，还会带来一系列环境问 题，如废石场的泥石流、尾矿库的渗漏和废水对环境的污染等[4-6].经过众多专家学者的努力，粗骨料膏 体充填技术成为解决矿山固体废物问题的有效方式，并在世界各地采矿实践中越来越受重视[7-8].采用粗骨料膏体充填，通常的做法是在地面充填站将 粗骨料、尾砂、胶凝材料和水混合制成复合材料，然 后通过管道系统输送至采空区[9].膏体输送至采空区一段时间后形成具有强度的胶结充填体，这可以限制围岩变形和提高矿石回收率.

管道输送是粗骨料膏体充填的关键环节[10].由于膏体质量分数较高，管输阻力大，导致管道输送难度偏大.因此，有必要对粗骨料膏体的输送特性进行研究.流变特性是膏体输送特性研究中的一个重要的分支，主要研究方法有L管法、倾斜管法和流变仪 测试等[11].其中L管法和倾斜管法由于流动距离、循 环时间等因素的影响导致测试结果存在较大偏差，流变仪测试成为比较普遍的研究方法[12].但是，相关研究表明流变仪测试所得流变参数与实际膏体参数存在较大差距.阻力损失也是表征粗骨料膏体可输 送性能的一个重要依据[13]，国内外对阻力损失的研究方法分为三种：经验公式法[14]、数值模拟法[15]和环 管试验法[16].相比经验公式和数值模拟而言，环管试 验充分考虑了各种因素对输送的影响，所得试验结果与工程实际的吻合度最高[17].例如：王洪江等[18]采用环管试验法对全尾砂料浆的输送特性进行分析，并得到了流变参数与灰砂比和浓度相关联的管道阻力预测方程.王勇等[19]揭示了料浆管输实际值和理论计算值之间存在的关系.郑伯坤等[20]研究了改性 全尾砂料浆在环管试验条件下的输送性能.但是，前 述研究主要侧重于全尾砂浆，粗骨料膏体在物料组成与性质上和全尾砂充填料浆存在着本质区别，基于全尾砂料浆所得试验结果并不能直接应用于粗骨 料膏体.

因此，本文以金川二矿区粗骨料膏体充填系统为研究背景，基于工业级环管试验系统对全尾砂-粗 骨料膏体的输送特性进行研究，分析了流变仪测试结果与环管试验结果的区别，并建立考虑流动度的输送阻力方程.通过响应面设计法（Box Behnken De-sign）对粗骨料膏体的管道输送参数进行设计，研究了单因素和多因素耦合对阻力损失的影响，并确定 各因素对阻力损失影响程度的大小，最后对管道输 送参数进行优化.研究成果可为粗骨料充填系统的设计提供理论依据.

1试验材料和方法

1.1试验材料

进行环管试验的试验材料均取自金川二矿区，其中细骨料为全尾砂，粗骨料为废石和棒磨砂.

1）骨料.金川二矿区全尾砂密度为2.645 g/cm3，对其进行XRF 元素测试分析（结果见表1），发现尾 矿中主要化学成分为SiO2、MgO.采用激光粒度仪测试得到全尾砂粒级分布曲线见图1（a），尾砂中极细 颗粒（粒径<20μm）质量分数为54.7%，粒径<45 μm 颗粒质量分数为77.3%，粒径<74 μm 颗粒质量分数为91%;计算得到尾砂的Cu 值为7.82，Cc 值为1.61，可知全尾砂级配良好.试验所用废石密度为2.876 g/ cm3，松散堆积密实度和密实孔隙率分别为0.582 4和31.57%.采用筛分法测试得到废石的粒级分布曲线（见图1（b）），可以看出废石的主要粒径分布集中在2～15mm，粒径小于10mm的细颗粒质量分数超过80%.棒磨砂密度为2.794 g/cm3，松散堆积密实度和密实孔隙率分别为0.5576和40.23%.筛分结果表明 棒磨砂颗粒比废石颗粒小，6 mm以下的颗粒质量分数为83.1%.

2）水泥.采用金川矿用非标38.5水泥（于充填站水泥 倉处取得），测得密度为3.105 g/cm3，容 重为1.12 g/cm3.

1.2  试验装置

1.2.1  环管试验系统

充填环管试验以高可靠度的方式对充填料浆在矿山现场类似工况条件下的运行参数进行模拟，为膏体充填工艺设计提供依据.本研究采用的环管试 验系统主要包括：配料机、计量秤、螺旋输送机、皮带给料机等设备，图2给出了环管系统结构示意图.试 验物料按照配比设计，通过配料系统进行精确计量，然后进入搅拌系统，同时根据设计浓度加入清水;搅 拌系统分为二段卧式搅拌，分别采用德国 BHS 卧式双轴强力搅拌机和双螺旋卧式搅拌机.泵送系统采用德国普斯迈斯特公司（Putzmeister）的KOS1070P 型柱塞泵，额定流量60m3/h，最大工作压力8.5mPa，泵送流量通过频率来调节.管道系统由3种尺寸的无缝钢管组成.料浆在环管系统中完成测试后返回双螺旋卧式搅拌机中，实现物料的闭路循环.通过隔 膜式压力传感器（共12个，编号为P1～P12，未全部在示意图中标明）采集料浆输送过程中管段内的压力值，数据采集频率为1次/s，通过DCS 控制系统实时采集、记录.

1.2.2 流变仪

选用产自丹麦的ICAR 混凝土流变仪对环管系统搅拌槽的取样浆体进行流变测试.ICAR的转子形式为四片装叶轮（转速为0.001～0.6 r/s），转子高12.7 cm、直径为6.35 cm，可用来测量骨料尺寸范围为6～32mm的粗骨料膏体.ICAR 混凝土流变仪的优点是可以对稳定阶段的流变特性曲线进行自动拟合，从而 计算出流变参数（初始剪切应力τ0和黏度系数μ）.

流变测试流程如下：首先设置转子转速为40r/s保持100s，然后放置 200s，以消除随机因素的干扰;然后设置转子转速在100s 内从0r/s增加到150r/s，在变剪切速率的条件下获取料浆的流变性能.

1.2.3水泥胶砂流动度仪

依据《水泥胶砂流动度测定方法》（GB/T 2419— 2005），選择 NLD-3型水泥胶砂流动度仪（跳桌）对粗 骨料膏体的振动扩散度进行测试，进而表征其流动性能.

1.3 环管试验方案

由于环管试验单次试验物料消耗大，不能实现所有配比的环管试验.为此，选择具有代表性的配比进行试验，以达到预测相近配比膏体管输特性的目的，试验方案如表2所示（由于实际质量分数不易控 制，表2中仅给出质量分数的调控范围）.环管试验的具体步骤在文献[16]和文献[18]中进行了详细的叙述，不赘述.下面针对本次环管试验过程中的工况参数调节过程进行说明.

1）设备调试.系统开始使用前先做清水试验，其 目的是检测密封性、检测仪器精度可靠性以及泵对流量的调节能力，检测数据如图3所示.

2）同配比料浆质量分数调节.通过添加外部清水进行浓度调节，试验配料质量分数设定为3种，在74%～78%范围内调节，每改变一次浓度，必须经过10～20min的连续循环混合，使膏体料浆混合均匀后 进行数据读取.

3）柱塞泵流量调节.通过调节泵送控制系统流量，泵送频率范围为55%～100%，对应理论流量范围应为24～58 m3/h，每种流量有效测试时间不少于5min.

2  基于环管试验的管输阻力计算

2.1 环管试验结果

经过环管输送试验后，选取压力表P8和P9所记 录的水平管段的数据进行处理.限于篇幅，本节内容 仅给出 m 尾砂 ∶m 废石 ∶m 棒磨砂为1∶0.5 ∶0.5的数据进行分析，其余组均具有相似的变化规律，且实际质量分数以现场取样测定为准（实测质量分数分别为73.4%、75.2%和77.1%）.图4为不同浓度条件下压力数据和流量数据的变化情况，可以看出：压力值随着泵送流量的降低而降低，且随着料浆浓度的增加，压力表P8和P9所测得的压力数据的差值越来越大.

2.2  环管试验结果拟合分析

由于压力传感器的敏感度较高，当废石、棒磨砂 粗颗粒撞击传感器触头时会造成数值的波动（见图4），这将导致测试所得数据因受到“噪声”的干扰而 较为离散、数据精度低.因此，对无效数据进行剔除 后，取膏体平稳运行段的数据进行数据拟合，得到不同配比条件下膏体输送阻力和流速之间的关系式，如表3所示.同时，因膏体管输阻力与流速间呈现显 著的线性关系，膏体充填料可近似视为宾汉塑性体，适用于管流的流变方程可用白金汉方程[21-22]描述：

将式（2）与表中拟合曲线方程进行联立求解，将 管道直径代入后就可获得实际条件下膏体的屈服应力和塑性黏度，计算结果见表3.

将环管测试所得流变参数与流变测试所得结果绘制成柱状图进行对比分析，如图5所示.

图5表明，流变仪测试所得流变参数普遍比环 管拟合参数要大，其中屈服应力相差10.4%～12.2%、塑性黏度相差 21%～32.8%.分析认为，在进行流变仪 测试时采用的是桨式转子，该转子可有效地减少壁 面滑移效应的发生，但是壁面滑移在实际输送过程中无法避免.同时，由于实际输送的时间远大于流变测试的时间，这就导致环管中料浆的触变性对实际 流变参数的影响较大.由图5 还可以看出屈服应力和塑性黏度均随着质量分数的增加而增大，且当质量分数固定时料浆的屈服应力随着物料组分的增加而降低.这主要是因为粗骨料膏体会形成一种由全尾砂包裹粗颗粒的稳定浆体.而由粗骨料粒径分析可知，棒磨砂6 mm以下的颗粒质量分数为83.1%，相对废石而言是一种细骨料.因此，当骨料为三组分时，废石和棒磨砂的级配更为合理，全尾砂浆体对二 者的悬浮包裹效果更佳，减少了颗粒间的运动阻力，表现为粗骨料膏体的屈服应力减小.

式（1）表明传统的阻力计算通常根据测出的流变参数并结合相应的输送参数（如管径、流速等）进行计算.然而矿山并不具备相应测试条件，通常委托 试验机构进行流变参数的获取[23].但是图5结果表明实验室所得参数和实际参数存在较大的误差.为了更加有效、准确且便捷地评估粗骨料膏体的输送 性能，在传统的阻力计算公式中引入相对容易获得的流动度参数.将膏体的流动度值与屈服应力、塑性黏度进行拟合，进而得到基于流动度的阻力计算新模型，流变参数拟合结果如图6所示（R2 均大于0.9）.将图6中的拟合方程代入式（1）即可得到不同骨料条件下的输送阻力经验公式，如式（3）～式（5）.

3 管道输送参数优化分析

3.1  响应面法试验水平设计

本节内容以全尾砂-废石-棒磨砂膏体为主要研究对象，选择管道直径、输送流量和固体质量分数作为自变量，输送阻力损失作为响应值，结合矿山实际 充填需求确定本次试验中管径及输送流量的取值范围，管道直径为140mm、180mm、220mm;输送流量为110m3/h、130m3/h、150m3/h，为了便于计算膏体的质量分数，根据环管试验时的参数进行取整，膏体质量分数分别为73%、75%、77%.最后，采用Design Expert 软件对方案进行设计，进而分析多因素交互 作用对阻力损失的影响，设计后的因素水平及编码 如表4所示.试验方案如表5所示，将试验所对应参数值代入式（5）中进行计算，分别得到不同条件下的输送阻力.

Design Export对表5中数据进行多元回归拟合后得到响应面函数，如式（6）所示.对所建立的响应 面回归模型进行可靠性检验（R2 检验），得到每个影 响因子的显著程度（P 值），发现所建立响应面模型的R2为0.9973且 P 值小于0.0001，说明响应面回归模型的回归效果极好、可靠性高.

3.2 单因素对膏体输送阻力的影响

图7（a）为膏体质量分数为73%时管道直径与阻力损失的关系.由图7（a）可知，当流量和质量分数不变时，随着管道直径的增大阻力损失呈下降的趋势，但是下降幅度在逐渐降低.这是因为管径增大时，管 壁剪切作用力减小，边界层的滑移效果减弱，此时阻力损失主要受到管内颗粒间摩擦作用的影響.管径增大时各流层间的距离增大[22]，颗粒的碰撞行为减 少，导致阻力损失快速减小;但是管径达到220mm 后临界速度发生变化，料浆流动形态由层流向紊流 转变，颗粒间的碰撞行为加剧，导致阻力损失的降低 幅度变缓.同时，结合表6发现管道直径 X1的F 值较大，说明管径对阻力损失的影响较为显著.图7（b）为膏体质量分数为73%时管道流量与阻力损失的关系.由图7（b）可知，当管径和质量分数不变时，随着 流量的增大，阻力损失呈上升的趋势.同时，结合表6发现，管道流量的F 值最小，说明其对阻力损失的影响程度最低.流速的增大导致膏体物料所承受的作用力发生变化，极大地提高了颗粒间的碰撞概率，相应的阻力损失增大.但是图7（b）曲线增长缓慢和F 值均表明流量对阻力损失的影响明显不如管径及膏体质量分数.图7（c）为管径110mm时质量分数与阻力损失的关系.由图7（c）可知，膏体质量分数与阻力损失呈正相关关系，且随着质量分数的增加，阻力损失的增幅变大.分析认为，膏体中的水分起到润滑的效果，当质量分数增加时相应的膏体含水量降低，导致浆体在流经管道时所受的摩擦增大，从而提高了膏体输送的阻力损失.在满足膏体可输送性能的同时适当提高质量分数可以降低水泥的用量，有助于降低生产成本.同时，结合表6发现，质量分数的F 值最大，说明其对阻力损失的影响程度最大，最终得到管径、流量和质量分数对阻力损失的敏感程度为：质量分数>管径>流量.

3.3 因素交互作用对输送阻力的影响

结合表6中模型参数的方差分析结果发现，X1X3的F 值最大，说明管径和质量分数的交互作用对膏 体阻力损失的影响最为显著.图8（a）为膏体质量分数为73%时，管径与流量间的交互作用对膏体管道 输送阻力损失的影响.可以看出，当流量较低（110m3/h）时，管径由140mm增大至220mm，阻力损失降 低了54.9%，这主要是由于各流层间的距离增大所致;当管径为180mm，流量从110m3/h增大至150m3/h时，阻力损失提升了9.35%，可以看出随着管径的增大，阻力损失对流量的敏感度逐渐降低，这主要是由于流量提升导致颗粒碰撞的效果小于流层间距增大导致的减阻效果.综上，在一定范围内同时增加管径和流量有助于矿山充填在提升产能的同时降低 损耗.图8（b）为流量110m3/h时，管径和质量分数的交互作用对膏体管道输送阻力损失的影响.图8（b）表明当管径为140mm，膏体质量分数由 73%增加至77%时，阻力损失提高了163.2%;而当质量分数为73%时，管径由140mm增大至220mm，阻力损失下降了55.1%.同时可以看出，阻力损失对管径的敏感度随着质量分数的增加而降低.这主要是因为膏体的流动性能主要依赖于其内部含水量，质量分数的增大导致含水量的减少，抵消了管径增大所带来的减阻效果.图8（c）为管径180mm时，流量和质量分数的交互作用对膏体管道输送阻力损失的影响.图8（c）表明当流量为110m3/h，膏体质量分数由 73%增加至77%时，阻力损失上涨了177.9%;而当质量分数为73%时，流量由110m3/h增大至150m3/h，阻力损失增加了9.89%.由于流量增大时会导致膏体内部 颗粒的紊动进而提高输送阻力，而质量分数的提高导致膏体内部固体颗粒增多，因此阻力损失对流量的敏感度随着质量分数的增加而变得显著.

3.4 响应面模型对输送参数的优化分析

金川全尾砂-废石膏体充填站设计充填能力200万 m3/a（两套充填管路），其中单套充填系统承担能力100万 m3/a，经换算后每小时的充填能力不应小于115m3，将其代入响应面模型进行优化设计，响应面 优化结果如图9和图10所示.可以看出，在管道参数为管径 220mm、流量142.4 m3/h、质量分数73.2%时，管道阻力损失最低为1.939kPa/m，此结果的可信度为0.996（最大值为1）.

4结论

1）利用白金汉方程对环管试验数据拟合后得到膏体的实际流变参数，发现流变仪测试所得流变参数普遍比环管拟合参数要大，其中屈服应力相差10.4%～12.2%、塑性黏度相差 21%～32.8%，这是由实 际输送过程中壁面滑移和触变性的影响所致.

2）由于流变测试所得参数和实际参数存在较大误差，在传统阻力计算公式中引入相对容易获得的流动度参数，建立了基于流动度的阻力计算新模型.

3）基于RSM-BBD开展全尾砂-粗骨料膏体输送阻力行为分析，研究了质量分数、管径和流量及其交 互作用条件下管道输送阻力的变化规律.得到了单 因素对阻力损失影响的敏感度大小为：质量分数>管 径>流量，并发现因素间的交互作用影响显著，其中质量分数和管径的交互影响起决定性的作用.

4）以金川二矿区粗骨料膏体充填系统为背景，采用响应面对输送参数进行优化，最终得到最佳的输送参数为：管径 220mm、流量142.4 m3/h、质量分数73.2%.该成果对于粗骨料膏体充填系统运行参数的选取具有一定的指导意义.

参考文献

[1]李夕兵，周健，王少锋，等.深部固体资源开采评述与探索[J].中国有色金属学报，2017，27（6）：1236-1262.

LI X B，ZHOU J，WANG S F，et al.Review and practice of deep mining for solid mineral resources [J].The Chinese Journal of Nonferrous Metals，2017，27（6）：1236-1262.（In Chinese）

[2]YANG L，YILMAZ E，LI J W，et al.Effect of superplasticizer

type and dosage on fluidity and strength behavior of cemented tail-ings backfill with different solid contents [J].Construction and Building Materials，2018，187：290-298.

[3]黄志伟，古德生.我国矿山无废开采的现状[J].矿业研究与开发，2002，22（4）：9-10.

HUANG Z W，GU D S.The state-of arts of wasteless mining in China[J].Mining Research and Development，2002，22（4）：9-10.（In Chinese）

[4]尹升华，刘家明，陈威，等.不同粗骨料对膏体凝结性能的影响及配比优化[J].工程科学学报，2020，42（7）：829-837.

YIN S H，LIU J M，CHEN W，et al.Optimization of the effect and formulation of different coarse aggregates on performance of the paste backfill condensation[J].Chinese Journal of Engineering，2020，42（7）：829-837.（In Chinese）

[5]BENZAAZOUA M，BUSSI?RE B，DEMERS I，et al.Integrated

mine tailings management by combining environmental desulphu-rization and cemented paste backfill：application to mine Doyon，Quebec，Canada [J]. Minerals Engineering，2008，21（4） ：330-340.

[6]杨晓炳，肖柏林，高谦.非牛顿流体充填料浆的管输摩擦阻力预测[J].湖南大学学报（自然科学版），2020，47（6）：125-131.YANG X B，XIAO B L，GAO Q.Prediction of friction pressure for non-Newtonian backfilling pipe-flow slurry [J].Journal of Hunan University（Natural Sciences），2020，47（6）：125-131.（In Chinese）

[7]杨天雨，乔登攀，王俊，等.废石-风砂高浓度料浆管道输送数值模拟及管输阻力新模型[J].中国有色金属学报，2021，31（1）：234-244.

YANG T Y，QIAO D P，WANG J，et al.Numerical simulation and new model of pipeline transportation resistance of waste rock-aeolian sand high concentration slurry[J].The Chinese Journal of Nonferrous Metals，2021，31（1）：234-244.（In Chinese）

[8]甘德清，闫泽鹏，薛振林，等.考虑壁面滑移效应的充填料浆管道输送阻力研究[J].金属矿山，2020（9）：26-32.

GAN D Q，YAN Z P，XUE Z L，et al.Study on the resistance of cemented paste backfill slurry pipeline transport considering wall slip effect[J].Metal Mine，2020（9）：26-32.（In Chinese）

[9]WU D，YANG B G，LIU Y C.Transportability and pressure drop

of fresh cemented coal gangue-fly ash backfill（CGFB）slurry in pipe loop[J].Powder Technology，2015，284：218-224.

[10]颜丙恒，李翠平，吴爱祥，等.膏体料浆管道输送中粗骨料颗粒运动规律分析[J].中南大学学报（自然科学版），2019，50（1）：172-179.

YAN B H，LI C P，WU A X，et al.Analysis of law of movement of coarse aggregate particles in pipeline transportation of paste[J].Journal of Central South University（Science and Technology），2019，50（1）：172-179.（In Chinese）

[11]胡华，孙恒虎，黄玉诚，等.似膏体粘弹塑性流變模型与流变方程研究[J].中国矿业大学学报，2003，32（2）：119-122.

HU H，SUN H H，HUANG Y C，et al.Rheological model and equation of viscoelastic-plasticity of paste-like backfill pulp[J].Journal of China University of Mining & Technology，2003，32（2）：119-122.（In Chinese）

[12]苗苗，雪凯旺，苗芳，等.石灰石粉对水泥浆体水化特性及流变性能的影响[J].湖南大学学报（自然科学版），2018，45（12）：90-96.

MIAO M，XUE K W，MIAO F，et al.Influence of limestone pow-der on hydration characteristics and rheological properties of ce-ment paste[J].Journal of Hunan University（Natural Sciences），2018，45（12）：90-96.（In Chinese）

[13]张修香，乔登攀 .粗骨料高浓度充填料浆的管道输送模拟及试验[J].中国有色金属学报，2015，25（1）：258-266.

ZHANG X X，QIAO D P.Simulation and experiment of pipeline transportation of high density filling slurry with coarse aggregates[J].The Chinese Journal of Nonferrous Metals，2015，25（1）： 258-266.（In Chinese）

[14]陈秋松，张钦礼，王新民，等.全尾砂似膏体管输水力坡度计算模型研究[J].中国矿业大学学报，2016，45（5）：901-906.   CHEN Q S，ZHANG Q L，WANG X M，et al.Pipeline hydraulic gradient model of paste-like unclassified tailings backfill slurry[J].Journal of China University of Mining & Technology，2016，45（5）：901-906.（In Chinese）

[15]WU A X，RUAN Z E，WANG Y M，et al.Simulation of long-distance pipeline transportation properties of whole-tailings paste with high sliming[J].Journal of Central South University，2018，25（1）：141-150.

[16]王少勇，吴爱祥，阮竹恩，等.基于环管实验的膏体流变特性及影响因素[J].中南大学学报（自然科学版），2018，49（10）： 2519-2525.

WANG S Y，WU A X，RUAN Z E，et al.Rheological properties of paste slurry and influence factors based on pipe loop test [J].Journal of Central South University（Science and Technology），2018，49（10）：2519-2525.（In Chinese）

[17]YANG X B，XIAO B L，GAO Q，et al.Determining the pressure

drop of cemented Gobi sand and tailings paste backfill in a pipe flow [J]. Construction and Building Materials，2020，255：119371.

[18]王洪江，王小林，寇云鹏，等.全尾砂高浓度胶结充填的环管试验[J].工程科学学报，2021，43（2）：215-222.

WANG H J，WANG X L，KOU Y P，et al.Loop test study on the high-concentration cemented filling of full tailings [J].Chinese Journal of Engineering，2021，43（2）：215-222.（In Chinese）

[19] 王勇，那庆，吴爱祥，等.粉煤灰-炉渣浆体环管实验和理论计算阻力对比研究）[J/OL].中南大学学报（自然科学版）：1-11（2021-07-26）[2021-10-12].

https：//kns.cnki.net/kcms/detail/ 43.1426.N.20210723.1649.002.html.

WANG Y，NA Q，WU A X，et al.Comparison of resistance be-tween pipe loop test and theoretical calculations in fly ash?slag slurry [J/OL].

Journal of Central South University（Science and Technology）：1-11（2021-07-26）[2021-10-12].

https：//kns.cnki.net/kcms/detail/43.1426.N.20210723.1649.002.html.（In Chinese）

[20]郑伯坤，姚维，黄腾龙，等.基于环管试验的改性全尾砂充填料浆输送性能[J].中国有色金属学报，2021，31（2）：520-529. ZHENG B K，YAO W，HUANG T L，et al.Transportable perfor-mance of modified unclassed tailings filling slurry based on loop test[J].The Chinese Journal of Nonferrous Metals，2021，31（2）： 520-529.（In Chinese）

[21]吴爱祥，刘晓辉，王洪江，等.考虑时变性的全尾膏体管输阻力计算[J].中国矿业大学学报，2013，42（5）：736-740.

WU A X，LIU X H，WANG H J，et al.Calculation of resistance in total tailings paste piping transportation based on time-varying be-havior[J].Journal of China University of Mining & Technology，2013，42（5）：736-740.（In Chinese）

[22]CHEN Q S，ZHANG Q L，WANG X M，et al.A hydraulic gradi-ent model of paste-like crude tailings backfill slurry transported by a pipeline system[J].Environmental Earth Sciences，2016，75（14）：1-9.

[23]侯永強，尹升华，戴超群，等.尾矿膏体流变特性和管输阻力计算模型[J].中国有色金属学报，2021，31（2）：510-519.

HOU Y Q，YIN S H，DAI C Q，et al.Rheological properties and pipeline resistance calculation model in tailings paste [J].The Chinese Journal of Nonferrous Metals，2021，31（2）：510-519.（In Chinese）