桃山矿胶结充填管道输送系统充填试验

聂文波 李凤义 陈 雷

(1.黑龙江科技大学安全工程学院，黑龙江 哈尔滨 150022,2.黑龙江科技大学矿业研究院，黑龙江 哈尔滨 150022)

桃山矿胶结充填管道输送系统充填试验

聂文波1李凤义2陈 雷1

(1.黑龙江科技大学安全工程学院，黑龙江 哈尔滨 150022,2.黑龙江科技大学矿业研究院，黑龙江 哈尔滨 150022)

为了解决桃山矿三井工作面充填开采过程中出现的堵管以及管道的高磨损率等问题，运用流体力学理论，结合该工作面的实际情况，采用自流+泵送的管道输送技术进行了工业试验，并利用ANSYS模拟分析了在速度2.0 m/s时，不同料浆浓度、管径对全程管压的影响。模拟结果表明：既定输送速度下，在合理的可控范围内且易于输送时的充填料浆浓度、管径取值范围分别为52%、140～180 mm时，输送管压分别为134.561 kPa、130.053～138.759 kPa。根据模拟结果进行的工业试验表明，堵管和高磨损率等问题得到了解决，充填开采取得良好的效果，有效地控制了顶底板的变形，降低了顶板的下沉量，减少对环境的破坏，验证了该管道输送系统是一种“安全、经济、环保、可靠”的充填采煤系统，为该矿以及类似条件下的矿山充填开采提供了借鉴。

充填开采 流体力学 管道输送 全程管压

充填开采作为绿色开采的一部分，因其在解决环境灾害、“三下”压煤问题上的突出效果[1],越来越受到各方面的重视。但在充填开采过程中管道输送出现的问题经常会导致采煤工作面开采能力下降，甚至会使煤矿开采无法进行。因此，充填管路输送问题也成为了影响煤矿开采与充填匹配的主要因素之一。

管道输送充填料浆技术以其具有的建造速度快、劳动强度低、运输速度快、连续性强且建造费用低、充填效果好等优点，成为充填开采的首选。常用的管道输送充填料浆的方法主要有自流输送和泵压输送2种[2-3]，料浆在管道中自流依靠的仅仅是料浆本身的重力，不需要外界动力源，降低了输送成本，被多数矿山采纳推广应用[4]。近半个世纪以来，国内外专家虽然对充填料浆管道流阻测试的研究[6]取得了很大进展，但由于充填材料的各异性、充填管路的特殊性、料浆流动的复杂性[6]，在实际充填项目操作的过程中，仍会出现管道高磨损、堵塞爆裂等问题。为提高充填料浆管道输送的安全可靠度，实现充填管道输送料浆的安全高效性，利用数值模拟软件ANSYS动态仿真模拟管道输送系统，查出原方案所存在的问题，并对参数优化，再利用模拟优化措施和实际经验对管道输送系统参数进行设计，最后将其应用到工业试验，成为矿山充填开采管道输送系统的主要发展方向和动态。

ANSYS 作为广泛应用的有限元软件充分综合了CAD等图像处理工具的优点[7-8]，是建立复杂计算模型有效而又方便快捷的平台，其含有的FLOTRAN流体分析模块[9-10]，通过模拟分析，对管道输送参数进一步优化，从而降低充填系统的事故率，节约生产成本。

1 充填工艺

1.1 充填材料选取

在前期新安矿巷旁充填[11]的基础上，黑龙江矿业研究院制定的桃山矿技术指标为：选取当地的粉煤灰，充填体能很快凝固膨胀实现主动接顶，凝结体的强度在6～8 h达到0.8～1 MPa，24 h后达到2～3 MPa，72 h后达到3.5～4 MPa，且强度还能够缓慢增长。即使充填体遭受顶板来压后内部出现裂纹，在集中应力转移以后，充填材料可以胶结自动愈合，且胶结后的抗压强度高于破坏前，以至于能良好有效地预防顶板覆岩移动。另外，还要求材料中的外加剂全部无毒，不能污染地下水源，真正实现绿色充填。

在满足充填体强度及料浆良好流动性的前提下，采用正交试验的方法和料浆扩散度实验方法对充填材料进行配比试验，最后获得了料浆浓度的合理范围为48%～58%。确定的固体物料合理配比(质量比)为：粗粉煤灰∶水泥∶石灰∶石膏∶早强剂∶J85∶硅酸钠∶元明粉=4 000∶500∶125∶75∶15∶10∶4∶20。

1.2 充填工艺

胶结充填能够在短时间内实现接顶，让充填体承压，这对于防止顶板冒落、采空区有害气体逸出等起到了关键作用。其充填工艺如图1。

图1 桃山矿三井胶结充填的工艺流程

2 数值模拟及结果分析

根据桃山矿实际情况及减压池减压机理，决定运用“自流加泵送”的管道输送方式。即井上活化站到井下制浆站采用自流输送系统，井上料浆在24 h内不凝固，避免由于长距离自流输送堵管的可能性。井下成浆制备到工作面充填采用泵送系统，能实现循环连续制浆，减少充填时间，提高矿山的充填效率[12]，更方便于管道输送充填系统在矿山上的推广应用，因此，本研究通过计算机软件模拟自流加泵送的管道输送方式，借以确定该方式下合理压力范围内适宜的料浆浓度、管路半径。

2.1 浆体参数

本次模拟的充填料浆主要有粗粉煤灰、水、水泥、活化剂和速凝剂等组成，水泥与粉煤灰的质量比为1∶8，选择浆体质量浓度Cw=52%。

(1)浆体密度ρJ。根据水泥与粉煤灰的质量比为1∶8，质量浓度为52%，各种活化剂和速凝剂的质量总和为0.038 66 t，则：

(2)浆体体积浓度CV。充填1 m3体积料浆需要水0.68 t，水的体积为VS=0.68 m3，可 得：

CV=32%.

(3)料浆黏度μm。以托马斯方程为基础，通过大量实验室试验，粗粉煤灰料浆固体物料的特性系数k、B分别取93、0.002 73。根据以上数据计算出质量浓度为52%的料浆与清水的黏度比：

0.322+93e0.002 73×0.32=95.72.

在实际测定下，该充填区气温均为20 ℃，清水的黏度1.005 mPa·s，粗粉煤灰料浆黏度μm=96.2 mPa·s[5]。

采用上面同样的方法，对质量浓度为48%、55%、58%料浆的密度、体积浓度和黏度进行计算，获得的参数如表1所示。

2.2 建立模型

利用分析软件ANSYS/FLOTRAN[13]进行数值模拟研究。桃山矿三井管道输送模型的建立以桃山矿充填系统管路铺设实际情况为原型，原型总长度为894 m，垂直高度差为300 m。基于模拟软件及计算机的局限性和矿山料浆管道输送系统的复杂性，考虑模型误差的影响，决定对管路系统进行简化。在假设浆体是宾汉体条件下进行料浆管道输送模拟，认定料浆黏性是恒定不变的，不考虑热交换，不考虑地压波等振动带来的影响。

表1 不同质量浓度下浆体的物理参数

2.3 数值模拟

在本次模拟过程中，选定浆体流速为2.0 m/s，从浆体浓度、管径2方面对输送系统进行模拟优化。

2.3.1 浆体浓度的模拟

模拟4组不同浆液质量浓度48%、52%、55%和58%的流动状态，模拟结果如图2(以质量浓度58%为例)和表2。

图2 浓度58%的管路全程压力和速度矢量

表2 不同质量浓度模拟的结果

2.3.2 管径的模拟

本次模拟选用100、140、160和180 mm 这4种直径进行模拟，浆液浓度定为52%，浆体流速为2.0 m/s，模拟结果如图3(以直径180 mm为例)和表3。

图3 180 mm管径的沿程压力分布和弯管处速度矢量

表3 不同管径模拟的结果数据

由图2和图3可知，管道最大压力随着料浆质量浓度和输送速度的增加而增大，但是料浆浓度对管压的增幅较小，管道最大压力仅增大4 kPa，确立最适宜的浓度为52%。管径越大，管压就越小，管径从100 mm升高到180 mm时，管道所受压力降了近28.7 kPa。而管径140～180 mm时，最大压力变化不大，故选择140～180 mm管径均可。

3 试验结果分析

充填体强度是衡量充填效果的重要指标，因此，充填体的强度是充填系统计算中要考虑的重要参数[16]。通过与以前充填项目中充填体质量的对比，本次充填项目中充填体的质量非常稳定、强度也增加了，对顶板及上覆岩层的控制取得了良好的效果。充填效果对比如图4所示。

图4 桃山矿工业试验充填前后效果对比

4 结 论

(1)采用矸石电厂粉煤灰，降低了充填成本，有效地解决了当地粉煤灰的污染问题，符合生态矿山的要求，同时得到的胶结充填材料的最佳质量配比，即粗粉煤灰(4 000)∶水(4 490)∶水泥(500)∶石灰(125)∶石膏(75)∶早强剂(15)∶ J85(10)∶硅酸钠(4)∶元明粉(20)。

(2) 管道最大压力在一定阶段内随管径增加而减小，随料浆浓度增加而增加，同时管道所受最大的压力与管径成非线性关系，但是在140～180 mm的管径对管道最大压力影响程度波动最小。

(3)充填体的良好密实性及高强性，不仅有效地阻止采空区有害气体的逸出，还在短时间内稳定了上覆岩层，为矿山安全开采提供了新的技术方向。

[1] 钱鸣高，许家林,缪协兴.煤矿绿色开采技术[J].中国矿业大学学报，2003,32(4)：343-348. Qian Minggao，Xu Jialin，Miao Xiexing.Green technique in coal mining[J].Journal of China University of Mining and Technology,2003,32(4)：343-348.

[2] 王新民,肖卫国,张钦礼.深井矿山充填理论与技术[M].长沙：中南大学出版社,2005:16-21. Wang Xinmin，Xiao Weiguo，Zhang Qinli.Deep Well Mining Filling Theory and Technology[M].Changsha:Central South University Press,2005:16-21.

[3] 郑晶晶,张钦礼,王新民.充填管道系统失效模式与影响分析(FMEA) 及失效影响模糊评估[J].中国安全科学学报,2009,19(6):166-171. Zheng Jingjing，Zhang Qinli，Wang Xinmin.FMEA analysis of backfilling pipeline system and fuzzy evaluation of failure effects[J].China Safety Science Journal,2009,19(19)：166-171.

[4] 李少辉,赵 澜，包先成，等.粉煤灰的特性及其资源化综合利用[J].混凝土.2010，25(4)：45-47. Li Shaohui，Zhao Lan，Bao Xiancheng，et al.Characteristics and comprehensive resources utilization of fly ash[J].Concrete,2010，25(4)：45-47.

[5] 金子桥,吕宪俊,胡术刚,等.膏体充填料浆水力坡度研究进展[J].金属矿山，2009(11):35-37. Jin Ziqiao，Lu Xianjun，Hu Shugang，et al.Development of research on hydraulic gradient of paste backfill slurry[J].Metal Mine，2010，25(4):35-37.

[6] 李永业，孙西欢,阎庆绂,等.不同导流条安放角条件下筒装料管道水力输送试验研究[J].水动力学研究与进展:A 辑,2008,23(1):86-89. Li Yongye，Sun Xihuan，Yan Qingfu，et al.Experimental research on the piped hydraulic transportation with the different diversion angle[J].Journal of Hydrodynamics:Ser.A,2008,23(1):86-89.

[7] 杨家冕,刘人恩,王 星.数值模拟在分层充填法采场参数选择中的应用[J].金属矿山，2013(3):29-35. Yang Jiamian，Liu Ren'en，Wang Xing.Application of numerical simulation to stope parameters selection of the slicing and back-filling mining method[J].Metal Mine，2013(3):29-35.

[8] 李 勇，刘志友，安亦然.介绍计算流体力学通用软件-FLUENT[J].水动力学研究与进展，2001，16(2):254-258. Li Yong，Liu Zhiyou，An Yiran．A brief introduction to FLUENT:a general purpose CFD code[J].Journal of Hydrodynamics,2001,16(2):254-258.

[9] 胡红军.ANSYS10.0 材料工程有限元分析实例教程[M].北京：电子工业出版社,2008:45-50. Hu Hongjun.Tutorial ANSYS10.0 materials engineering finite element analysis of an example[M].Beijing:Publishing House of Electronics Industry,2008:45-50.

[10] 谢龙汉.ANSYS FLOTRAN流体及热分析[M].北京:电予工业出版社，2012:89-92. Xie Longhan.ANSYS FLOTRAN Fluid and Thermal Analysis[M].Beijing:Publishing House of Electronics Industry,2012:89-92.

[11] 李凤义,王伟渊.新安煤矿井下沿空留巷巷旁充填实验[J].黑龙江科技大学学报,2013,23(5):409-412. Li Fengyi，Wang Weiyuan.Filling test filled in gob-side entry retaining in Xin'an Coal Mine[J].Journal of Heilongjiang Institute of Science & Technology,2013,23(5):409-412.

[12] 陈维新,关显华.粉煤灰基胶结充填采煤技术及应用[J].煤矿现代化,2014,120(3):15-18. Chen Weixin，Guan Xianhua.Fly ash based cemented filling mining technique and application[J].Coal Mine Modernization,2014,120(3):15-18.

[13] 袁 博，张召千，张百胜.破碎围岩注浆加固数值模拟分析与工程应用[J].金属矿山,2013(7):45-48. Yuan Bo，Zhang Zhaoqian，Zhang Baisheng，et al.Numerical simulation of grouting reinforcement of cracked surrounding rocks and its application in engineering[J].Metal Mine，2013(7):45-48.

(责任编辑 徐志宏)

Back-filling Experiment of Cemented Filling Slurry Pipeline System for Taoshan Mine

Nie Wenbo1Li Fengyi2Chen Lei2

(1.SafetyEngineeringInstitute，HeilongjiangUniversityofScienceandTechnology，Harbin150022，China;2.MiningResearchInstitute，HeilongjiangUniversityofScienceandTechnology，Harbin150022，China)

In order to solve the problems of pipe blocking and high wear rate of pipeline during filling mining process at Sanjing working face of Taoshan Mine，industrial tests of gravity and pumping transport are carried out by adopting the fluid mechanics theory，in combination with the practical situation of the working face.The ANSYS simulation was made to analyze the effect of different slurry concentration and pipeline diameters at the given speed 2.0 m/s on the entire pipe pressure.Simulation results showed that under the given delivery speed，when the filling slurry concentration values at 52% and pipe diameter ranges from 140 mm to 180 mm，which are within the reasonable control and makes slurry transport easy，the pipeline pressure are 134.561 kPa，and ranges from 130.053 to 138.759 kPa respectively.Industrial tests together with the simulation indicated that the problems of pipe blocking and high rate of pipeline of slurry transportation are solved and backfill mining obtained good effect.The pilpline system can effectively control the deformation of roof and floor，reduce the roof subsidence and lower the environment.It is approved that it is a kind of “safe，economical，environment-friendly and reliable” filling mining system.It provides reference for filling mining of the mine and similar mines.

Filling mining，Fluid dynamics，Pipeline，The entire pressure pipeline

2014-05-04

黑龙江省发展高新技术产业(非信息产业)专项(编号:FW12A018)。

聂文波(1987—)，男，助教，硕士研究生。

TD823.7

A

1001-1250(2014)-09-035-04