FDEM-CFD coupling analysis of spatiotemporal evolution of mining-induced overburden fracture-water inflow in shallow and thick coal seam under river
YANG Zebin;LI Hao;MA Liqiang;MAI Longquan
河下浅埋厚煤层采动覆岩裂隙分布与涌水量是工作面安全生产的决定性因素之一。数值模拟是二者重要预测方法,其合理性的关键在于建立岩体破坏−裂隙流体耦合理论及相应的模拟方法。以路家村矿15404工作面为研究背景,构建拉/剪应力下非贯通裂隙开裂、贯通裂隙的法向与切向本构关系,并基于二相流质量守恒、动量守恒与状态方程,结合增强的浸没边界算法识别流−固界面,通过流体体积法实现裂隙内流体自由面的追踪和重构。在此基础上形成FDEM-CFD河下采煤覆岩裂隙与涌水量预测数值模型耦合程序,通过相邻工作面地面钻孔冲洗液消耗量观测法验证导水裂隙带发育高度,以及采用大井法理论对涌水量结果进行对比。结果表明:采动岩体破坏−裂隙流体耦合理论及相应的FDEM-CFD程序,可数值实现河下浅埋厚煤层采动覆岩裂隙形成,以及裂隙内流体运移过程。当工作面推进至80~120 m时覆岩内形成贯通地表的导水裂隙。招山河水的主要下泄路径为工作面后方8~20 m位置处的导水裂隙,其斜向采空区、倾角65°~72°。模拟所得采空区涌水量为18.78 m3/h,与大井法计算结果接近。上述成果在路家村矿区得到初步应用,为进一步开展河下浅埋厚煤层防治水工程提供理论支撑。
The distribution of overlying rock fractures and water inflow during mining in shallow and thick coal seams under the river is one of the decisive factors for the safe production of working faces. Numerical simulation is an important prediction method for both, and the key to its rationality lies in the establishment of a rock mass failure fracture fluid coupling theory and corresponding simulation methods. Taking the 15404 working face of Lujiacun Mine as the research background, the normal and tangential constitutive relationships of non through crack cracking and through crack under tensile/shear stress are constructed. Based on the conservation of mass, momentum, and state equations of two-phase flow, the enhanced immersion boundary algorithm is combined to identify the fluid solid interface. The fluid volume method is used to track and reconstruct the fluid free surface inside the crack. On this basis, a coupling program of FDEM-CFD numerical model for predicting overlying rock fractures and water inflow in coal mining under the river is formed. The development height of the water conducting fracture zone is verified through the observation of the consumption of flushing fluid in adjacent working faces, and the results of water inflow are compared using the large well method theory. The results indicate that the coupling theory of mining rock mass failure and fracture fluid, as well as the corresponding FDEM-CFD program, can numerically achieve the formation of overlying rock fractures in shallow and thick coal seams under rivers during mining, as well as the fluid transport process within the fractures. When the working face advances to 80-120 m, a water conducting crack that runs through the surface is formed within the overlying rock. The main discharge path of the Zhaoshan River is the water diversion fissure located about 8-20 m behind the working face, which is inclined towards the goaf and has an inclination angle of about 65°-72°. The simulated water inflow in the goaf is 18.78 m3/h, which is close to the calculation results of the large well method. The above achievements have been preliminarily applied in the Lujiacun mining area, providing theoretical support for further carrying out water prevention and control projects in shallow and thick coal seams under the river.
shallowly buried coal seam beneath a river;overlying strata fractures;water inflow;constitutive model;FDEM-CFD coupling
主办单位:煤炭科学研究总院有限公司 中国煤炭学会学术期刊工作委员会