考虑壁面滑移效应的充填料浆管道输送阻力研究

考虑壁面滑移效应的充填料浆管道输送阻力研究

近年来，浅部矿产源逐渐减少，矿山开采深度不断增加，“三公”（高温、高压、深井）问题日益严重。影响管道输送阻力的因素众多，由于受到测量手段和试验装置限制，常规的环管试验受到多种条件制约，粗略的测量值很难作为论证依据由于高浓度料浆属于非牛顿体范畴，因而其在管道输送过程中必然存在表观滑移现象1数学模型1.1模型假设管道内部主流区与滑移层接触部分的流体速度分布是连续的，壁面滑移速度是剪切流动区和滑移层连接处的速度与管道内壁的速度差1.2流场基本流动方程(1）滑移流动的基本方程。根据管流中剪切速率的定义γ=4vR及管流基本公式，得出管径为R的管道中料浆产生壁滑移时的基本流动方程式中，R为管道半径，m;v为浆体平均流速，m/s；τ为壁面剪切应力，Pa；τ(2）黏性耗散。管道输送过程存在的黏性耗散量是一个与速度梯度相关的物理量，可表示为式中，ϕ为黏性耗散量；η为流体黏度，Pa/s;u、v、w分别为X、Y、Z轴方向的速度分量，m/s;T为料浆温度，℃；1.3流动量表的控制方程对于不可压缩的流体，Comsol软件中控制流体流动的基本方程主要有能量守恒方程、质量方程和状态方程，控制流动的变量有流体压力P、密度ρ、温度T和流速u。Navier-Stokes动量守恒方程（忽略外部体积力）可表示为式中，u质量守恒方程可表示为式中，ρ为料浆密度，kg/m状态方程可表示为式中，ρ为料浆密度，kg/m2模拟试验设计2.1试验材料本研究全尾砂样品取自河北某矿选矿厂，从底堆处取料，经烘干后备用测定其密度为2850kg/m2.2该项目的内部设计根据国内矿山自流充填系统的发展情况根据雷诺数Re计算公式2.3网格剖分单元本研究基于Comsol软件三维空间维度对充填管路进行建模，弯管曲率半径为0.5m，由于充填倍线为5，设置竖直管长度为1.5m，水平管长度为9.5m。Comsol软件自带的网格剖分工具有很强的编辑能力，用户可根据需求调整网格的形状、大小等。考虑到滑移层极薄，为使计算结果更精确，设置模型的边界层数为5，边界层拉伸因子为2，最终剖分单元数目为350623。建模效果如图1所示。3基于数值模型的可靠性验证和结果分析3.1数值模拟结果与公式计算结果对比吴爱祥等式中，ΔP为管道压力损失，Pa;D和L分别为管道直径和管道长度，m；μ文献式中，i为验证Comsol数值模型在充填料浆管道输送阻力损失研究方面的可靠性，将上述试验条件代入式（7）和式（8）进行计算，对比分析数值模拟结果和公式计算结果，如表2所示。由表2可知：式（8）计算结果与数值模型的计算结果相对误差为48.1%～386.1%，误差极大。原因是该式在多次试验的基础上所得，试验条件对公式的适用性影响较大，且该式未考虑壁面滑移对管内料浆流动的影响；式（7）计算结果与数值模拟计算结果的相对误差在20%以内，说明本研究试验所建立的考虑壁面滑移效应的数值模型相对可靠。3.2壁面滑移效应试验结果分析将数值模拟数据导入DPS数据处理软件中，采用极差分析法研究料浆浓度、灰砂比、管径对阻力损失的影响程度，并建立了阻力损失关于三者的回归方程。试验结果如表3所示。由表3可知：考虑壁面滑移效应时，各因素对管道阻力的影响程度依次为管径>质量浓度>灰砂比。通过回归拟合建立的阻力损失与浓度、灰砂比、管径的函数关系式为式中，Y为阻力损失，Pa/m;x4填充料浆管道输送阻力损失的变化特征4.1管道输送阻力管径与管道阻力损失之间的关系如图2所示，可以看出管径与阻力损失呈负相关。图2(a）表明管径增大时，灰砂比对阻力损失影响的显著程度减小，3条曲线逐渐发散；图2(b）表明，灰砂比不变，浓度增大时，料浆流动的“黏性”增大，输送阻力越大。高浓度料浆中尾砂颗粒迁移主要受到壁面剪切作用力影响，管径增大，壁面剪切作用力减小，颗粒的迁移运动变缓，滑移效应减弱，对摩擦阻力的抵消作用减小。管径较小时，滑移层的减阻效果明显，可以看出，管径为100～150mm时，管道输送阻力快速降低；当管径增大到一定值时，料浆中的尾砂颗粒所受到的剪切作用力较弱，致其无法进行迁移运动，不能形成覆盖管道内壁的滑移层，此时滑移作用效果较弱，随着管径增大，管道输送阻力缓慢减小。4.2浓度对输送阻力的影响料浆浓度与管道阻力损失之间的关系曲线如图3所示。由图3可知：料浆浓度与管道阻力损失基本呈正相关。图3(a）表明，管径不变料浆浓度增加时，管道输送阻力逐渐增大，同时可以看出灰砂比对输送阻力影响较大，水泥含量越高，阻力损失值越大。图3(b）反映灰砂比一定时浓度对阻力损失的影响规律，管径越大，管道阻力损失越小，同时随着浓度增大，图中3条曲线开始聚拢，说明浓度变化导致输送阻力的差异越来越小。浓度较低时，料浆管道输送过程中滑移层的形成主因为静态壁面损耗效应4.3管道阻力下降砂灰比与管道阻力损失之间的关系如图4所示，可以看出砂灰比与管道阻力损失呈负相关。图4(a）表明，管径一定时，管道阻力损失随着砂灰比增大而呈线性减小；图4(b）表明料浆浓度一定时，料浆中的水泥含量降低，管道阻力损失快速下降，当水泥含量低于一定值后，管道阻力损失减小的速率快速减小。水泥遇水发生水化反应的过程会大量消耗料浆中的自由水，使得高浓度料浆浓度增大，表观黏度及屈服应力增加，导致料浆的流动性变差，料浆的黏性耗散量增大。同时由于滑移层的组成成分主要为自由水，当大量的自由水被消耗时，不能形成可以填充整个管道粗糙内壁的滑移层，此时管道输送阻力较大。当水泥含量下降（砂灰比增大）时，水化反应需水量下降，料浆流动性较好，可以形成有效的滑移层，此时管道输送阻力快速下降；当水泥含量较少时，其水化反应消耗的水对滑移层结构的影响较小，管道输送阻力降低速率减小。由图4(b）可知：料浆在管道输送过程中，灰砂比与阻力损失的关系曲线存在明显转折，当灰砂比超过1∶8时，阻力损失的降低速率明显变小，曲线趋于水平，灰砂比低于1∶10时会影响充填体强度的形成，因此本研究推荐的最佳灰砂比范围为1∶8～1∶10。5充填料浆配比参数对充填管道输送阻力的影响冀东地区某矿为地下矿山，由于生产需要建有1座地面充填站，标高+110m。该矿山采用全尾砂胶结充填技术进行充填作业。充填选用的尾砂来自选厂浓缩池，经浓缩后在搅拌桶中与水泥混合制备成浓度为70%～73%、灰砂比为1∶4～1∶6的充填料浆。充填料浆经砂浆泵加压后通过充填钻孔（钻孔长度110m）进入地下采场。其中，充填管路内径180mm，砂浆流量为100～140m根据矿山充填料浆实际配比参数选择灰砂比为1∶6，浓度分别为72%、73%的料浆进行模拟试验。同时根据前述分析所得的参数最优取值范围，选择灰砂比为1∶8，浓度为66%、68%、70%的料浆进行对比验证，不同料浆所对应的流变参数如表4所示。选择浆体初始速度为1.5m/s进行模拟试验，结果如图6所示。由图6可知：当采用矿山实际充填参数时，管道输送阻力较大。矿山充填管道的输送阻力随着料浆浓度增加而快速增大，当料浆浓度小于70%时，可大大减小料浆管道输送过程中的能量损耗，同时也会降低堵管、爆管的风险。考虑到矿山充填结构体的强度需求，为有效提高矿山生产效率和减小渗水量，充填料浆浓度不宜低于66%，因此，得到本试验条件下料浆的最佳输送浓度为66%和68%。灰砂比由1∶6降低为1∶8时充填管道的输送阻力减小，但是水泥含量过低时会影响矿山充填体早期强度形成。当灰砂比大于1∶10时，水泥含量变化对充填体强度的影响较小，同时基于管道输送阻力的变化情况得出该试验条件下的最佳灰砂比为1∶8。6数值模型的建立(1）通过Comsol软件建立了考虑壁面滑移效应的管道输送模型，将数值模拟结果与长沙矿冶研究院的经验公式和考虑滑移效应的推导公式进行比较，相对误差分别为48.1%～386.1%和4.8%～16.7%，证明本研究所建立的数值模型用于计算考虑壁面滑移效应的阻力损失具有可靠性。(2）基于壁面滑移效应分析管道阻力损失的变化特征，认为阻力损失与灰砂比、浓度呈正