一种计算变质量破碎岩体渗透性参量的方法-应用力学学报.PDF

第 31 卷 第 6 期 应 用 力 学 学 报 Vol.31 No.6 2014 年 12 月 CHINESE JOURNAL OF APPLIED MECHANICS Dec. 2014 文章编号：1000- 4939(2014) 06-0927-06 一种计算变质量破碎岩体渗透性参量的方法 陈占清 1,2 王路珍 1,2,3 孔海陵 3 倪晓燕 1,2 姚邦华 4 (1 中国矿业大学 深部岩土力学与地下工程国家重点实验室 221008 徐州；2 中国矿业大学力学与建筑工程学院 221008 徐州； 3 盐城工学院土木工程学院 224051 盐城；4 河南理工大学安全科学与工程学院 454000 焦作) 摘要：建立了一种基于压力梯度和渗流速度时间序列提取变质量破碎岩体渗透性参量（渗透率、 非 Darcy 流β 因子、加速度系数）的方法。首先，根据破碎岩体渗透性参量之间的幂指数关系和 Forchheimer 关系，建立了采样时刻渗透率的代数方程，并利用Newton 切线法求得此代数方程的 根；其次，对渗透性参量的参考值及幂指数进行优化。通过算例分析了该方法用于计算变质量破 碎岩体渗透性参量的可行性及准确性。研究结果表明：随着细小颗粒的迁移流失，破碎岩体在恒 定的压力梯度下，渗流速度和渗透率增加，非 Darcy 流β 因子和加速度系数降低。 关键词：变质量；破碎岩体；渗流；渗透性参量；方法 中图分类号：TU452 文献标识码：A DOI ：10.11776/cjam.31.06.A026 析突水机理，并提出相应的防治措施。因此，研究 1 引 言 质量迁移和流失过程中破碎岩体渗透性参量的变化 规律，对于研究破碎岩体渗流的非线性行为及渗流 在破碎岩体的渗透过程中，由于水的溶蚀、磨 失稳诱发的灾害具有重要的工程意义。 蚀、冲刷作用，渗透质（破碎岩体）中细小颗粒随 在早期的岩石渗透试验中，人们考虑了变形引 渗透剂（水）一起迁移。这种质量迁移必然造成破 起的孔隙度变化，但未考虑岩样内部细小颗粒的迁 碎岩体孔隙度、渗透性参量的时变[1] 。根据微分动 移和流失引起的孔隙度的变化，故渗透性参量可以 力系统理论，渗透性参量决定破碎岩体渗流系统的 通过线性回归而获得。在标准圆柱形岩样破坏过程 拓扑结构，当渗透性参量发生变化时，渗流系统因 （应力-应变全程）的瞬态渗透试验中，文献[11]利 为结构失稳而产生分岔和突变等行为[2-4] 。破碎岩体 用液体性质建立岩样中质量流量与压力（梯度（变 中质量迁移和流失，不仅决定渗流场的时变，而且 化率））之间的关系式（流体的质量守恒方程），基 还会引起流动模式的转换，即由渗流转换为管流[5] 。 于流体的质量守恒方程和动量守恒方程，通过线性 渗流失稳造成的动力灾害在采矿工程中已经 回归可提取出岩样的渗透率、非Darcy 流β 因子、 得到重视，相关学者从断层活化[6-7] 、含水层和/或 加速度系数；文献[12]利用孔隙压力差单一时间序 其邻近岩层的变形和破坏[8-9]、陷落柱质量流失和载 列提取非 Darcy 流渗透性参数。在质量恒定的破碎 荷变化[1,10]等方面解释采动岩体渗流失稳原因，分 岩体稳态渗透试验中，岩样的渗透率和非 Darcy 流 基金项目：国家重点基础研究发展计划项目(973 计划)(2013CB227900) ；国家自然科学基金委员会与神华集团有限责任公司联合资助项目(U1261201) ； 国家自然科学基金(5097