数值模拟技术在风氧化岩裂隙水充水矿区水文地质条件定量计算与分析应用研究.docVIP

数值模拟技术在风氧化岩裂隙水充水矿区水文地质条件定量计算与分析应用研究 　　摘要：陕北矿区井下主要水害类型为直罗组风氧化岩裂隙水，其补给条件复杂，给矿井实际生产过程中带来巨大威胁。为了准确掌握该种含水层的水文地质条件，利用刻画地下水系统空间结构和水力特征的数学模型作为工具，以数字模拟方法为手段来定量分析、评价、预测地下水系统的水文地质条件、参数结构、行为规律及其在扰动条件下的变化与响应。分析结果可以更加科学精确的指导煤矿企业在实际生产过程中对类似水害采取有效措施。 　　关键词：数值模拟 矿井涌水量 水文地质 陕北矿区 　　中图分类号：P641.4 　　一、 前 言 　　陕北矿区煤层赋存稳定，地质结构简单，充水水源包括地表水、松散含水层水、基岩含水层水和采空区积水等。陕北矿区大多矿井的主要水害威胁来自基岩风氧化裂隙水，其补给条件较为复杂??补给源丰富、宽泛，大多无法采取常规的堵、截等措施，只能采取强行疏排的方式。为了有效指导矿井疏排工作，确保煤矿安全生产，杜绝由于排水能力不够引发水害事故的发生，就需要对主要含水层的水文地质条件进行进一步分析，尽量准备的计算矿机开采过程中的矿井涌水量，从而使矿井排水系统的设计依据更加充分、科学。 　　二、 数学模型与数值方法 　　采用国际上先进的三维有限差分软件Visual MODFLOW来模拟陕北某煤矿S1210工作面放水试验过程（放水I、II、III、 IV及水位恢复阶段）地下水系统的动态变化，通过放水试验成果进一步校正模型，进而以校正模型为基础预测地下水流场发展趋势与涌水量计算预测，分析矿井水文地质条件和水害问题。 　　1、三维流数学模型 　　根据地下水流动的水文地质概念模型，区内放水试验过程中含水系统地下水流动的三维非稳定流数学模型可描述如下： 　　边界条件： 　　式中， x，y，z为笛卡尔坐标轴；t为时间；H0为初始统测水头，H为已知水头；Kxx为坐标轴方向的主渗透系数；？s为比弹性给水度；？d为重力给水度；W为单位体积井流量，抽水时取负号；Γ1为第一类边界；Γ2-1为潜水面边界；Γ2-2为零流量边界；ε’为降雨入渗补给量。 　　2、 数值方法 　　有限差分方法（FDM）是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，用网格节点上的函数值的差商代替控制方程中的导数进行离散，从而建立以网格节点上的值为未知数的代数方程组。基于质量守恒（水流连续性方程）与渗流基本定律（达西定律）的有限差分与有限元技术的数学微分方程（数学模型）表达式相同。 　　有限差分方法的基本思想是：用渗流区内的有限个离散点的集合代替连续的渗流区，在离散点上用差商近似代替微商，将微分方程及其定解条件化为未知函数在离散点上的近似值为未知量的代数方程，然后求解差分方程，进而得到微分方程的解在离散点上的近似值。该方法是根据地下水流动的连续性方程进行的，按照连续性方程，流入和流出某个计算单元的水量之差等于该单元储水量的变化。 　　地下水连续性方程可表示为： 　　式中：――单位时间内流进或流出该计算单元的水量； 　　SS――含水层的贮水率； 　　――计算单元的体积； 　　――单位时间内水头的变化。 　　将渗流区进行剖分离散后，可以确定计算单元，根据连续性方程及达西公式： 　　可以得到在行方向上由计算单元（i，j-1，k）流入单元（i，j，k）的流量，表达为： 　　式中：――水头在格点（i，j，k）处的值； 　　――水头在格点（i，j-1，k）的值； 　　――通过格点（i，j，k）和格点（i，j-1，k）之间界面的流量（L3 T-1）； 　　―― 格点（i，j，k）和格点（i，j-1，k）之间的渗透系数（LT-1）； 　　―― 横断面面积（L2）； 　　――格点（i，j，k）和格点（i，j-1，k）之间的距离（L）。 　　通过其他五个界面的地下水流量，均可类推，例如，沿行方向的格点（i，j+1，k）至格点（i，j，k）的地下水流量为 　　同理，可依次得出： 　　以上公式表示了通过计算单元（i，j，k）的六个界面的地下水流量，将格点间距和渗透系数合并为一个变量，成为水力传导系数： 　　将水力传导系数（12）应用于公式（6）～（11）得到： 　　这些公式用来计算单元（i，j，k）的六个边界面流入的地下水流量，此外，还需要考虑其他各种外部源和汇对计算单元的影响，例如河流，沟渠，生产井，注水井，蒸发蒸腾等，这些源汇流入单元的量可用一个通式表达： 　　式中： ―― 第n项外部源对计算单元（i，j，k）的补给量（L3 T-1）； 　　，――常数。单位为（L2 T-1），（L3