一、引言：氯化钾生产工艺优化与模型精准化的技术需求
氯化钾作为保障国家粮食安全的战略性资源，其生产工艺的持续优化对农业可持续发展具有重要意义。光卤石冷分解结晶法是生产氯化钾的主流技术路线，近年来面临原料品位下降、粒度波动加剧等挑战，导致光卤石分解率降低、产品纯度下降、粒度分布恶化等问题。在光卤石冷分解结晶过程中，氯化钾晶体的沉降速度直接决定结晶器内的颗粒分级效果与产品排出效率，是导流筒-挡板结晶器结构设计与操作参数优化的核心依据。然而，现有固-液曳力模型在宽粒径、高固含率条件下的预测精度不足，难以满足工业级结晶器的精准设计需求。本文基于25°C条件下氯化钾沉降实验数据，系统对比分析现有曳力模型的适用范围，提出同时考虑晶体粒径与固相浓度的晶浆表观黏度修正模型，显著提升了氯化钾沉降速度的预测精度，为2026年氯化钾行业的结晶器优化设计与工艺控制提供理论支撑。
二、氯化钾沉降实验设计与测试方法
(一)氯化钾实验体系构建
《2026-2031年中国氯化钾行业重点企业发展分析及投资前景可行性评估报告》实验针对光卤石连续冷分解工艺中光卤石充分溶解、氯化钾晶体稳定析出的连续排料过程，首先配制25°C条件下钠-钾-镁-氯-水四元体系饱和母液。采用化学滴定与离子色谱双重校核分解液组成，确保实验体系的准确性。利用质量-体积法与旋转黏度计分别测定母液密度与黏度，随后将母液转移至双层控温量筒中，通过循环控温设备精确维持25°C恒温环境。
(二)氯化钾颗粒分级与沉降测试
采用电动振筛机对氯化钾颗粒进行分级，获取25-54μm、100-250μm、250-500μm三个特征粒径区间的颗粒样品。鉴于氯化钾颗粒呈现规则立方体形貌特征，采用上下限筛网孔径的算术平均值计算颗粒的中值粒径，分别得到39.5μm、175μm、375μm三个代表性粒径。
将定量颗粒倒入双层控温量筒内搅拌均匀并开始计时，观测浆体清浑水交界面，记录不同时刻下的澄清层高度，以交界面停止下降为沉降终点，通过Talmage-Fitch法计算颗粒沉降速度。每个工况重复测试8-10次，取平均值并确定误差范围。实验系统考察固相质量分数(5%、10%、15%、25%)与平均粒径(39.5μm、175μm、375μm)对氯化钾沉降速度的影响规律。
三、氯化钾沉降模型对比分析与偏差识别
(一)现有曳力模型的适用性评估
固-液两相流中颗粒沉降速度的计算主要依赖曳力模型的选择。典型模型包括Stokes定律、Syamlal-O'Brien模型、Wen-Yu模型及Song模型等。Stokes定律最初描述小球在黏性流体中的低速运动，后经拓展适用于不同雷诺数范围，但仅适用于固含率极低、颗粒间相互作用可忽略不计的体系。Song模型考虑了颗粒球形度对沉降速度的影响，但未考虑固含率的影响。Syamlal-O'Brien模型与Wen-Yu模型引入了与固含率相关的修正因子，其中Syamlal-O'Brien模型在高固含率体系适用性更强，Wen-Yu模型能够准确描述固含率在20%以内的颗粒沉降速度。
(二)氯化钾沉降模型的预测偏差分析
将各文献曳力模型的计算结果与实验数据进行对比。结果表明，在饱和母液中，当颗粒平均粒径从39.5μm增至375μm时，各固含率体系的平均沉降速度由7.70×10⁻⁵m/s增加至5.33×10⁻³m/s，粒径增大显著促进沉降速度提升。当固相质量分数从5%提高至25%时，各粒径颗粒的平均沉降速度由2.94×10⁻³m/s下降至9.11×10⁻⁴m/s，高固含率条件下颗粒间相互作用增强、晶浆表观黏度增加，导致动量传递过程中的动能耗损加剧，沉降速度显著降低。
更重要的是，当固相质量分数超过10%时，各模型的预测准确度均显著下降。传统Wen-Yu模型在固相质量分数为5%、10%、15%、25%条件下，对氯化钾沉降速度预测的平均偏差分别达到2.02×10⁻³m/s、2.28×10⁻³m/s、2.68×10⁻³m/s、3.23×10⁻³m/s，最大偏差超过2×10⁻³m/s，难以满足结晶器精准设计的需求。
四、氯化钾晶浆表观黏度模型的修正与验证
(一)氯化钾晶浆表观黏度的粒径效应
传统固-液体系表观黏度计算公式仅考虑固含率的影响，而实验表明相同固含率下颗粒粒径对表观黏度同样产生显著影响。细颗粒体系具有更大的比表面积，颗粒间相互作用更强，导致表观黏度显著高于粗颗粒体系。
基于实验数据，提出同时考虑液相黏度、固含率与颗粒粒径的晶浆表观黏度修正模型。模型引入粒径相关修正因子，通过最小二乘法回归确定模型参数。修正后的模型计算值与实验值的平均偏差由修正前的5.14×10⁻³Pa·s下降至1.25×10⁻³Pa·s，显著提升了晶浆表观黏度的预测精度。
(二)氯化钾表观黏度修正模型的数学表达
修正后的晶浆表观黏度模型将固相体积分数与颗粒粒径纳入统一框架，模型参数通过实验数据非线性回归确定。结果表明，在相同液相黏度和固含率条件下，颗粒直径越小，悬浮液的表观黏度越大，39.5μm细颗粒体系的表观黏度显著高于175μm与375μm粗颗粒体系。修正模型能够准确捕捉这一粒径效应，为后续曳力系数修正奠定基础。
五、氯化钾沉降模型的曳力系数修正与精度提升
(一)氯化钾沉降模型的Wen-Yu-μs修正
将修正后的表观黏度模型代入Wen-Yu曳力系数模型，用晶浆表观黏度替换雷诺数计算中的液相黏度，得到等效雷诺数，进而修正曳力系数表达式。通过GRG非线性规划求解法对模型参数进行数值回归求解，确定修正后的Wen-Yu-μs模型参数。
修正后的曳力系数模型表达式为：
CD' = (24/αlRes') × [1 + 3.02(αlRes')^0.81]
其中Res'为基于表观黏度计算的等效雷诺数，αl为液相体积分数。该模型同时考虑了固含率与粒径对曳力系数的耦合影响。
(二)氯化钾沉降预测精度的显著提升
将Wen-Yu-μs模型应用于不同固含率与粒径条件下氯化钾沉降速度的预测。结果表明，特别是在粒径大于300μm且固相质量分数大于10%的工况下，修正模型的预测能力比原模型显著提升。
具体而言，Wen-Yu-μs模型对固相质量分数为5%、10%、15%、25%晶浆体系的颗粒沉降速度预测平均偏差，分别降至8.00×10⁻⁴m/s、4.66×10⁻⁴m/s、2.15×10⁻⁴m/s、6.20×10⁻⁴m/s。对于39.5μm、175μm、375μm颗粒沉降速度预测的平均偏差，分别降低至3.65×10⁻⁵m/s、7.44×10⁻⁴m/s、7.95×10⁻⁴m/s。整体预测平均偏差从原模型的2.55×10⁻³m/s下降至5.25×10⁻⁴m/s，精度提升约79%。
六、氯化钾结晶器设计的理论支撑与应用展望
(一)氯化钾结晶器颗粒分级机制的模型支撑
导流筒-挡板结晶器是光卤石冷分解生产氯化钾的典型设备，其工作原理依赖颗粒沉降速度差异实现分级：大粒径晶体沉降至结晶器底部从产品出口排出，细小晶体沿环隙上升至溢流区进入细晶消除系统，中等粒径颗粒在导流筒内循环实现持续生长。准确的氯化钾沉降速度预测是合理设计设备结构、优化操作条件、实现流场稳定与颗粒分级的前提。
修正后的Wen-Yu-μs模型能够精确预测宽粒径分布、高固含率条件下氯化钾晶体的沉降速度，为结晶器导流筒直径、挡板位置、循环流量等关键结构参数的优化设计提供可靠的理论依据。模型预测精度的提升有助于缩小设计裕度，降低设备投资与运行成本，提高氯化钾产品质量与生产稳定性。
(二)氯化钾生产工艺的智能化控制基础
精准的沉降速度模型为氯化钾冷分解结晶过程的智能化控制奠定基础。通过实时监测晶浆固含率与粒度分布，结合修正模型预测不同粒径颗粒的沉降行为，可实现结晶器排料时机、细晶消除强度、循环流量等操作参数的动态优化。这对于应对光卤石原矿粒度波动、稳定氯化钾产品质量具有重要意义。
七、总结
2026年，氯化钾行业面临资源品位下降与品质要求提升的双重挑战，光卤石冷分解结晶工艺的精细化控制成为技术升级的关键。本文针对光卤石冷结晶体系中氯化钾沉降速度预测难题，基于25°C条件下系统实验数据，识别出现有曳力模型在固相质量分数超过10%时预测偏差超过2×10⁻³m/s的技术瓶颈。通过引入晶体粒径与固相浓度耦合影响的晶浆表观黏度修正模型，对Wen-Yu模型曳力系数进行修正，开发出Wen-Yu-μs新模型。验证结果表明，修正模型将氯化钾沉降速度预测的整体平均偏差从2.55×10⁻³m/s降至5.25×10⁻⁴m/s，精度提升近80%，特别在粒径大于300μm、固含率大于10%的工业关键工况下优势显著。研究成果为光卤石冷分解结晶器的结构设计优化与操作参数调控提供了精准的理论工具，有助于提升氯化钾生产效率、产品质量与工艺稳定性，对推动钾肥生产关键工艺技术进步、保障国家粮食安全战略资源供应具有重要实践价值。未来研究将针对固含率超过30%的更高浓度体系，结合颗粒间作用力模型进行更深入探究，进一步完善氯化钾沉降预测的理论体系。
