3. 膏体充填料浆的流变模型

膏体充填料浆不是清水，也不是稀砂浆。它由尾砂、水泥、粉煤灰等细颗粒与水混合而成，浓度高、颗粒间作用强，在外力作用下表现出明显的非牛顿特征。要计算管输阻力，第一步就是选对描述“剪切应力—剪切速率”关系的流变模型。模型选错，屈服应力和粘度参数就失去物理意义，后续沿程阻力计算也会系统性偏离。

流变测试通常用旋转流变仪或管式流变装置，得到一系列剪切速率 γ̇（单位 s⁻¹）与对应剪切应力 τ（单位 Pa）的数据点。把这些点画在 τ–γ̇ 坐标上，观察曲线形状，再决定用哪个模型拟合。下面这张图对比了四种典型模型的曲线形态及其参数含义。

牛顿流体：理想参照，不是膏体的真实描述

牛顿流体的本构关系是 τ = μ·γ̇，剪切应力与剪切速率成正比，比例系数 μ 为动力粘度，曲线是过原点的直线。水、低浓度砂浆在工程近似上可视为牛顿流体。膏体料浆由于固相体积分数高、颗粒间存在絮网结构，即使在极低剪切速率下也需要一定应力才能启动流动，因此 τ–γ̇ 曲线不过原点，牛顿模型只能作为对照基准，不能用于膏体阻力计算。

宾汉姆模型：最常用的工程近似

宾汉姆模型写作 τ = τ₀ + μₚ·γ̇。它假设料浆存在一个屈服应力 τ₀，只有当外加剪切应力超过 τ₀ 时料浆才开始流动；流动后应力与剪切速率呈线性关系，斜率即塑性粘度 μₚ。对膏体充填料浆，宾汉姆模型在中等剪切速率区间（约 1～60 s⁻¹，接近实际管道剪切范围）拟合效果较好，参数少、物理意义直观，是矿山充填管输阻力计算中应用最广的模型。

参数物理意义要分清：屈服应力 τ₀ 反映料浆形成结构网络的强度，与浓度、细颗粒含量、水泥水化程度、静置时间密切相关。τ₀ 越大，料浆越“稠”、启动越困难，管道停输后再启动的峰值压力也越高。塑性粘度 μₚ 反映流动后内部剪切阻力的大小，主要受颗粒浓度和颗粒级配影响，浓度升高 μₚ 显著增大。工程上常用坍落度筒和扩展度试验间接判断 τ₀ 与 μₚ 的变化趋势：坍落度小、扩展度低，通常对应高 τ₀ 和高 μₚ。

赫巴模型：兼顾屈服与剪切变稀

赫巴模型写作 τ = τ₀ + K·γ̇ⁿ，在宾汉姆模型基础上把线性项换成幂律项。它同时保留屈服应力 τ₀、稠度系数 K 和流变指数 n 三个参数，能描述“有屈服应力且流动后仍呈非线性”的料浆。膏体料浆在较宽剪切速率范围内常表现出剪切变稀（n<1），此时赫巴模型拟合精度优于宾汉姆模型，尤其适合流变曲线弯曲明显、剪切速率跨度大的测试数据。代价是参数增多，拟合和物理解释更复杂，工程应用不如宾汉姆模型普遍。

幂律模型：无屈服应力的非线性描述

幂律模型写作 τ = K·γ̇ⁿ，曲线过原点，不含屈服应力项。K 为稠度系数，数值上等于 γ̇ = 1 s⁻¹ 时的剪切应力，反映料浆整体稠稀程度；n 为流变指数，n<1 为剪切变稀（假塑性），n=1 退化为牛顿流体，n>1 为剪切增稠（膨胀性）。幂律模型适合描述高分子溶液、部分细颗粒悬浮液，但对具有明显屈服行为的膏体料浆，低剪切速率段会低估应力，因此多用于近似分析或特定浓度区间的拟合，不宜作为充填管输阻力计算的主模型。

模型选择与参数获取的实际做法

实际工作中，先做流变测试得到 τ–γ̇ 散点，再用最小二乘法分别拟合四种模型，比较决定系数 R² 和残差分布。若数据在低剪切段明显不过原点、中高剪切段近似直线，优先选宾汉姆模型；若整条曲线弯曲显著，选赫巴模型；幂律模型仅在确认无明显屈服行为时使用。拟合得到的 τ₀、μₚ（或 K、n）代入管输阻力公式前，要核对测试剪切速率范围是否覆盖管道实际剪切速率范围，避免外推。

归纳起来，四种模型构成从简到繁的描述谱系：牛顿模型是基准，宾汉姆模型是膏体工程计算的主力，赫巴模型用于提高宽剪切范围的拟合精度，幂律模型适用于无屈服行为的近似分析。理解每个参数的物理意义，比记住公式本身更重要——它决定了你能否判断拟合结果是否合理、能否把实验室参数正确迁移到管道阻力计算中。