大面积充填采矿时的流固耦合数值分析及优化.pdfVIP

第24卷第4期 矿业T程研究 V01．24No．4 2009年12月 Mineral Research Dec．2009 Engineering 大面积充填采矿时的流固耦合数值分析及优化 余伟健他，高谦3 (1．湖南科技大学煤矿安全开采技术湖南省重点实验室．湖南湘潭411201； 2．湖南科技大学能源与安全工程学院，湖南湘潭411201；3北京科技大学土木与环境工程学院，北京100083) 摘要：针对大面积充填采矿涌水型矿井的水文地质条件。进行了流固耦合的数值分析．首先，阐述了有限差分法的流固耦合的 计算基本原理；然后。根据司家营铁矿南区矿体的特点，合理地选取了各计算参数，确立了4因素3水平共9个计算方案，并分别对这 些方案进行了采矿与充填的数值模拟分析；最后。针对采矿过程中的渗流场和采场涌水量进行分析，并得到了不同方案对地表沉降的 影响程度，从而最终优选出最佳方案．该研究成果为大水型矿山的大面积开采设计和安全评价提供一定的参考价值． 关键词：大面积充填采矿；流固耦合；数值分析；优化 中图分类号：TD853．34 文献标识码：A 目前，在进行大面积地下采矿时，采场安全和围 对孔隙的影响可以通过流动本构定律加以反映，反之 岩稳定性成为研究的热点，其影响因素错综复杂，其 孔隙压力的变化可引起固体产生变形．在FLAC3D的 中之一就是水的影响．在地下水和降水较多的地区，进 有限差分计算中，将质量平衡方程和流体矢量与孑L隙 行大面积采矿时，不但会诱发岩层的移动，而且涌水 压力之间的关系运用到均质、各向同性的恒流单元 多也会直接影响矿山生产效益．唐山钢铁司家营铁矿 体．对于均质、各向同性固体和恒流液体，达西定律可 南区的开采就面临着这个问题，其地下水非常充足， 表示为岫 采矿坑内涌水量较大．水和岩石的相互作用势必会影 g{=一蠡p—pf菇{g,l， (1) 响大面积采场围岩及顶板的稳定性和安全性，有必要 式中，吼：流体的矢量，m／s；p：孔隙压力，MPa；xi：坐标矢 进行专门的流固耦合分析，并达到优化充填采矿．目前 有关地下工程的流固耦合问题大多采用数值计算方 重力矢量的3个分量，m／s2． 法进行【l-31，因此，本文根据流固耦合作用机制，采用 对于小变形问题，流体的平衡可由下式表达： nAC30进行分析与优化． 等}=_gl+g。 (2) 1 流体一固体相互作用数值表达 式中，g，：流体源头流动强度，l，s；f：流量的变化，即每 单位孔隙介质中由于流体扩散造成的流体体积变化． 在地下工程中，除开采诱发的动态应力外，渗透 渗流的本构方程可表示为 压力也会对岩石造成变形损伤，不能忽略导水裂隙生 (3) 成过程中水压力的影响．这里涉及2种耦合：(1)应力一 am．C=寺等+a等书誓， 渗流耦合问题．这种问题是将孑L隙介质中流体压力与 应力建构在统一的本构模型中，平行计算应力场与渗 热系数，l，℃． 流场；(2)固一液耦合．这种耦合的实质是材料体积应变 将(3)式代入(2)式，可以导出： 收稿日期：2009—03—17 基金项目：国家科技支撑计划课题资助项11I(2008BAB32801、 作者简介：余伟健(1978-)，男，