适用读者： 煤矿矿长、总工程师、生产科长、地面注浆站/充填站长、管道输送设计工程师。
作者： 张洁贞｜中矿天智信息科技（徐州）有限公司高级销售经理｜中国矿业大学技术支持
一、 行业痛点与背景
在煤矿“无废消纳”与“采空区充填治理”工程中，浆体管道输送（包含高水充填、膏体充填、矸石浆体覆岩隔离注浆等）是连通地面制浆站与井下空区的唯一大动脉。
然而，在实际运营期间，几乎所有的充填项目都遭遇过噩梦般的“堵管事故”：
- 经济损失惨重：一旦发生管道结死，整条管线必须切段清理或直接报废更换，一次堵管可能导致数十万元至上百万元的直接损失。
- 干扰生产接续：主斜井管线被堵，会导致井下采煤面被迫停产等待，打乱矿井的采掘平衡。
- 设备磨损剧烈：若盲目提高输送流速以防堵，由于管道内壁磨损率与流速呈 3 ~ 4 次方的非线性关系，会使高强度钢管在短短数月内被磨损穿孔。
管道堵塞与剧烈磨损的本质，是水力学计算参数的严重失配。
许多设计单位直接套用清水的流动阻力公式，或者采用拍脑袋的经验系数，忽略了高浓度矸石浆体（含有微细及粗粒级）的非牛顿流体流变学本质。
因此，进行严谨的流速、浓度、管径最优水力学匹配计算，是保障项目长周期安全运行的基石。
二、 矸石高浓浆体的流体力学特性与本构模型
矸石浆体（多为破碎矸石粉体掺入粉煤灰或水泥的混合浆液）在管道中的流动特性，受浓度、级配与流速影响，呈现非牛顿流体性质。
2.1 流变本构模型
通常，质量浓度在 70% 以上的均质浆体表现为宾汉流体性质：
其中：
- τ 为管道剪切阻力（Pa）。
- τ_y 为屈服应力（Pa），表现为启动流动的初始阈值。
- η_p 为塑性粘度（Pa·s），为流动后的粘性摩阻。
这两个关键参数在实验室可通过旋转粘度计测定，它们随质量浓度 C_w 呈指数上升趋势：
其中，A, B, C, D 为特定材料物理化学性质决定的系数。由该关系可知，浓度微小的波动（例如从 75% 升到 78%），会导致屈服应力及摩阻力暴增数倍。
三、 临界沉降流速与摩阻水力学计算模型
要实现“不堵管”，首要条件是流速必须高于颗粒开始在管底沉降积聚的“临界流速”。
3.1 临界沉降流速计算 (Durand-Condolios 方程)
对于含有一定比例粗颗粒（>0.1mm）的全级配矸石浆体，其不发生管底沉降的临界流速 v_L（m/s）主要采用经典的杜兰德（Durand）修正公式计算：
其中：
- F_L 为 Durand 阻力系数，是与固体体积浓度及颗粒粒径 d_{50} 相关的常数，通常取值范围为 0.6 ~ 1.4。
- g 为重力加速度（9.81\text{ m/s}^2）。
- D 为管道内径（m）。
- \rho_s 为矸石颗粒真实比重密度（约 2.6 \times 10^3\text{ kg/m}^3）。
- \rho_f 为携带液（水或细泥浆）的密度（kg/m³）。
为了保障管道运行的绝对安全余量，系统正常运行的输送流速设计值 v_0 必须满足：
一般设计工作面水平管道流速范围在 1.2 ~ 1.8 m/s 之间。
3.2 水平管道沿程压力降阻力损失计算
高浓矸石浆体在管路内的沿程阻力降（压力损失 \Delta P_f，Pa）根据达西-魏斯巴赫（Darcy-Weisbach）方程计算：
其中：
- \lambda 为达西摩擦系数，对于宾汉层流流态，\lambda 是雷诺数 Re_B 与赫巴数 He 的函数：
- L 为输送管道长度（m）。
- \rho_m 为矸石浆体的均质湿密度（kg/m³）。
四、 流速-浓度-管径三维最优匹配设计表
为帮助现场技术人员进行参数速查，我们基于 Durand 方程与 Bingham 阻力模型，对常用的 DN125 与 DN150 两种超高压无缝耐磨管道，在不同浓度下的沿程压力损失（摩阻，kPa/m）进行了对账标定：
4.1 DN150 (内径 150mm) 管道水力学计算匹配表
*(注：计算基于矸石颗粒中位粒径 d_{50} = 0.15\text{mm}，泵送排量 120m³/h)*
| 浆体质量浓度 (C_w) | 临界沉降流速 (v_L, m/s) | 推荐设计流速 (v_0, m/s) | 屈服应力 (τ_y, Pa) | 塑性粘度 (η_p, Pa·s) | 沿程阻力降 (kPa/m) | 每百米摩阻压力降 (MPa/100m) |
|---|---|---|---|---|---|---|
| 70% | 0.98 | 1.25 | 18.2 | 0.08 | 0.85 | 0.085 |
| 74% | 1.15 | 1.45 | 38.5 | 0.18 | 1.62 | 0.162 |
| 76% (最优工作点) | 1.24 | 1.55 | 58.2 | 0.28 | 2.15 | 0.215 |
| 78% | 1.35 | 1.68 | 95.4 | 0.42 | 3.42 | 0.342 |
| 80% (高粘度憋压区) | 1.48 | 1.82 | 165.0 | 0.68 | 5.85 | 0.585 |
4.2 DN125 (内径 125mm) 管道水力学计算匹配表
*(注：泵送排量 85m³/h)*
| 浆体质量浓度 (C_w) | 临界沉降流速 (v_L, m/s) | 推荐设计流速 (v_0, m/s) | 屈服应力 (τ_y, Pa) | 沿程阻力降 (kPa/m) | 每百米摩阻压力降 (MPa/100m) |
|---|---|---|---|---|---|
| 70% | 0.90 | 1.15 | 18.2 | 1.12 | 0.112 |
| 74% | 1.05 | 1.35 | 38.5 | 2.10 | 0.210 |
| 76% (最优工作点) | 1.12 | 1.45 | 58.2 | 2.95 | 0.295 |
| 78% | 1.22 | 1.55 | 95.4 | 4.85 | 0.485 |
由上表对比可知：
- 在相同浓度下，管径越细（DN125），沿程摩阻上升显著（DN150 76% 摩阻为 2.15 kPa/m，而 DN125 达 2.95 kPa/m）。这意味着若输送距离超过 2 公里，DN125 管道对地面柱塞泵的主动驱动压力要求更高。
- 浓度一旦达到 78% - 80%，由于流体剪切变稠，每百米阻力急剧跨越至 0.5 MPa 以上，极易发生泵送超载憋压或剪切变粘导致堵塞。
五、 管道设计工艺与耐磨防堵设备选型
要保证长距离输送系统的寿命，在水力学最优参数设计之外，必须引入机械耐磨和流程防堵设计。
5.1 弯头选型
在立管转水平管等高压部位，推荐使用高强度耐磨陶瓷内衬弯头（硬度大于 HRA85）或双金属复合耐磨弯头，设计弯曲半径比 R/D ≥ 10，实现渐变平稳导流，阻断剧烈磨损。
5.2 自动清管系统 (Pigging System)
在地面制浆车间及井下出口处配置专用的自动清管阀与高压气动快速清管球仓。在每班作业结束后，直接通过管口发射聚氨酯高弹清管球，利用压缩空气或高压清水将管道内壁的挂浆扫清，从根本上杜绝二次凝固性结垢堵管。
六、 现场停电停泵等突发事件应急处置流程
系统由于故障引发突发停泵时，为防范浆体在管道中静态脱水沉降而导致凝固，必须执行严格的事故排空与冲洗联动流程。
`mermaid
graph TD
A[异常：主泵断电/井下管路破裂/超限报警] --> B{识别危机状态}
B -- 主泵断电停机 --> C[瞬时激活UPS备用电源控制柜]
C --> D[连锁触发地面事故清水高压泵开阀]
D --> E[开启垂直立管顶部自流旁路水阀]
E --> F[利用垂直管道段自重差静压高速向下驱扫管内残浆]
F --> G[开启井下紧急排空防爆阀, 将残液直排入采空区事故储水池]
B -- 管道磨穿漏浆 --> H[立即停止地表柱塞泵送]
H --> I[开启孔口倒流泄压阀]
I --> J[通知井下检修工区实施快速管道段更换与封胶]
`
七、 实践案例与技术验证
以山西某 350 万吨/年重点绿色矿山矸石充填项目（DN150 管道，垂直深度 480 米，水平管网长度 3200 米）为例：
7.1 匹配参数确定与改造
前期由于级配不合理且直接套用清水设计，工作流速设为 2.8 m/s，浓度 68%，导致 DN150 无缝钢管仅运行 4 个月 弯头即被整体冲穿；随后因浓度波动至 79% 发生大面积立管堵死事故，清理工期长达 12 天。
后经水力学重新核算：
- 级配修正：筛分粒径限制在 d_{50} = 0.18\text{mm}，加入 5% 的减水剂以调低粘度。
- 流速调整：按 Durand 公式，临界流速为 1.25 m/s，设计工作流速降至 1.55 m/s（大大减缓了摩擦强度）。
- 浓度固定：质量浓度固定在 75.5% ± 0.5%。
7.2 匹配改造前后运行数据分析表
| 监控评价指标 | 改造前（参数失配状态） | 改造后（最优水力匹配状态） | 优化改善幅度 |
|---|---|---|---|
| 管道平均磨损率 (mm/10万方) | 3.85 (磨损严重) | 0.42 (轻微擦伤) | \downarrow 89.1\% |
| 每百米管路压力降 (MPa/100m) | 0.485 (高粘憋压) | 0.205 (温和流态) | \downarrow 57.7\% |
| 突发性堵管频次 (次/年) | 8 | 0 (未发生) | 彻底消除 |
| 主泵能耗费用 (万元/年) | 148 | 82 | \downarrow 44.6\% |
| 管道使用周期 (年) | 0.5 | 3.5 | \uparrow 600\% |
八、 前期准备资料清单
进行长距离管道水力学计算前，设计院及厂家需要矿方准备以下技术数据：
- [ ] 矸石样品级配累积曲线：提供由粒径分析仪测出的矸石、粉煤灰、胶凝材料粒径分布图表。
- [ ] 浆体流变参数实验室检测报告：不同浓度（70%、72%、74%、76%、78%）下的屈服应力与塑性粘度数据。
- [ ] 井上下三维管网拓扑走向数据表：各节点垂直落差、水平长度、水平转角弯头度数。
- [ ] 高压柱塞充填泵额定曲线：主泵最大排量、最大出口压力以及电机变频调节范围。
九、 常见问题 FAQ
Q：为什么流速并不是越高越安全？流速高了不是更不容易堵管吗？
A： 这是一个流体力学误区。流速高确实能防止固体沉降，但它会带来两个严重的副作用：
- 压力降急剧上升：管道沿程摩擦摩阻与流速的 2 次方 呈正比。流速翻倍，每百米阻力上升四倍，主泵可能超压憋死停机，进而被动导致堵管。
- 机械磨损以指数级暴增：钢管的磨损量与流速的 3.6 ~ 4.0 次方 呈正比。高流速下，管道壁会极快变薄穿孔。
因此，水力学设计的原则是：在高于临界沉降流速 0.3 - 0.5 m/s 的“安全下限”附近运行，这才是能耗与寿命的最优解。
Q：减水剂在浆体输送中起什么作用？真的能防止堵管吗？
A： 可以。高固相矸石浆液之所以阻力大，是因为微细颗粒表面带电，易聚集成“絮凝结构”包裹水分，增大了剪切屈服阻力。
加入适量的化学减水剂后，减水剂分子会定向吸附在固体颗粒表面并带同种电荷，通过电荷之间的排斥力拆散“絮凝结构”，将包裹的自由水释放出来。这能在不改变质量浓度（不降低最终强度）的前提下，将浆体的屈服应力调低 50% 以上，大大拓宽了安全泵送的窗口期。
Q：如何实时判定管道内部是否已经开始发生颗粒沉积？
A： 现场采用“压力梯度法”与“流量差对比”进行在线诊断。
我们在管道进出口及井下中段配置高精度防爆压力变送器。如果监测到相邻变送器之间的压力差在排量不变的情况下持续攀升，说明管道有效过流面积在减小，底层已经开始发生矸石“滑动床”沉积。
此时，控制系统会自动联锁发出警报，并指令主泵系统自动变频提速冲洗管路，防止沉积恶化为结死堵管。
参考依据
- 国家规范 GB/T 50546-2025《矿山管道输送工程设计规范》
- 《非牛顿多相浆体管道输送水力学》（科学出版社）
- 中国矿业大学管道输送研究所《大功率高固相柱塞泵沿程压力传递与阻力特性科研报告》
关于作者
张洁贞，中矿天智信息科技（徐州）有限公司高级销售经理。
专注于陕蒙晋矿区煤矸石充填管网阻力计算、泵站大推力选型以及在线智能防堵 SCADA 控制平台设计。
联合中国矿业大学重介输送团队，协助煤矿集团保障大动脉的长期无阻安全运行。
下一步行动
📋 获取《长距离高压泵送管道水力学计算参数现场复核清单》
📊 获取《基于Durand临界流速与Bingham摩阻的管网阻力匹配计算器Excel》
💬 咨询张洁贞 —— 扫描右侧微信二维码（微信：zhangjiezhen176），预约中国矿业大学教授团队进行管道设计评审。
*本文首发于 [zhangjiezhen.cn](https://zhangjiezhen.cn/blog/slurry-pipe-hydraulics.html)，转载请注明出处。*
*声明：文中技术参数与案例数据经脱敏处理，具体工程应用须结合现场特定地质及水文接续条件进行定向设计。*