南水北调工程采空区三维宾汉姆流体紊流模拟 Three-dimensional turbulent numerical simulation of Bingham fluid in the goaf grouting of the South-to-North Water Transfer Project.pdfVIP

水 利 学 报 2013年11月 SHUILI XUEBAO 第44卷 第11期 文章编号：0559-9350（2013）11-1295-08 南水北调工程采空区三维宾汉姆流体紊流模拟 王晓玲，王青松，周正印，敖雪菲 （天津大学 水利工程仿真与安全国家重点实验室，天津 300072） 摘要：南水北调中线工程穿越多个煤矿采空区，对工程的安全施工与运行造成了严重威胁，因此对煤矿采空区进 行注浆加固处理已成为工程迫切需要解决的问题。本文根据注浆浆液的宾汉姆流体特性，建立宾汉姆流体本构模 型方程。基于采空区三维几何地质信息模型，提出三维宾汉姆流体k-ε紊流数学模型。针对注浆工程中岩层、断 裂带、冒落区、采空区、煤层底板等地质层采用不同的地质参数、孔隙率和阻力系数，对宾汉姆流体浆液在采空 区中的紊流扩散进行模拟，实现了宾汉姆流体浆液的采空区多孔分序注浆模拟，并分析了浆液的扩散规律。结果 表明：在一定时间范围内，浆液扩散半径和采空区的充填率随时间的增加而增加，但其变化率呈逐渐减少的趋 势，并最终趋于零，数值模拟的注浆量与实际注浆量相比误差为11.8%。研究可为采空区注浆的设计与优化提供 理论依据与技术支持，对确保南水北调渠道长期安全稳定运行具有重要的理论意义和实践价值。 关键词：注浆；宾汉姆浆液；k-ε紊流模型；采空区；南水北调 中图分类号：O373 文献标识码：A 1 研究背景 南水北调中线一期工程总干渠途经多个煤矿采空区，严重威胁着工程的安全施工和运行，因此 对采空区进行注浆加固处理已成为工程迫切需要解决的问题。注浆技术已被广泛地应用于地基工 程、桩基工程、边坡工程、隧道工程和煤矿开采，由于注浆具有掩蔽性和岩石地层的复杂性，使得 越来越多的学者对注浆技术进行了研究［1］。数值模拟技术比传统试验方法具有更大的灵活性、经济 性，能综合考虑更多的影响因素，随着计算机技术的发展，已经越来越多的被应用到注浆工程中， 显示出巨大的优势。 目前，国内外在注浆数值模拟方面主要是以裂隙或者单个的注浆孔为研究对象，而针对施工面 积大、注浆孔众多、注浆工序繁杂、复杂地质条件的实际注浆工程进行注浆模拟的较少。在数学模 型上，多数研究认为注浆液流体特性符合非牛顿宾汉姆流体，且流态近似做层流处理，采用达西定 ［2］ ［3］ 律，并在此基础上推导出了一维、二维或者三维的渗流控制模型。在国外，Houlsby 、Yang 、Fi⁃ ［4］ ［5］ delibus 等研究了裂隙网络中浆液流动过程的模拟方法。Bouchelaghem 等选用渗流模型进行了以细 沙为多孔介质的注浆实验，并建立了二维和三维数学模型分别模拟浆液的渗透规律。Wang［6］等进行 了以土壤为介质的注浆试验，并将介质等效为连续介质，采用二维有限元模型对注浆试验模型进行 ［7］ ［8］ ［9］ 了模拟。在国内，杨米加 、湛铠瑜 、罗平平 等在裂隙结构模拟的基础上，认为浆液流态为层 ［10］ ［11］ 流，基于质量守恒和达西定律，建立了裂隙网络非牛顿流体渗流模型。王万顺 、孙斌堂 、程鹏 ［12］ ［13］ 达 、闫常赫 等基于达西定律和质量守恒定律，建立了注浆渗流有限元模型。上述国内外的研究 收稿日期：2012-12-12 基金项目：国家自然科学基金创新研究群体科学基金；国家自然科学基金项目；“十二五”国家科技支撑 计划重点项目（2011BAB10B06） 作者简介：王晓玲（1968-），女，重庆人，博士，教授，主要从事环境工程模拟与优化研究。E-mail：wangxl@ 通讯作者：周正印（1986-），男，山东人，博士生，主要从事水环境工程模拟与优化研究。E-mail：zhouzhengyin@126.c