数值方法在尾矿库调洪演算中的应用及研究初探.docVIP

数值方法在尾矿库调洪演算中的应用及研究初探 　　摘要：在某些矿山中，尾矿库对矿山安全具有重要作用。通过数学化、根据已知条件最终确定尾矿库调洪库容和设计坝高。提出具有较强实用性的调洪演算方法，并可根据条件不同及时变化。 　　关键词：尾矿库、排洪、调洪演算 　　1.前言 　　尾矿库调洪演算的目的是根据既定的排洪系统确定所需的调洪库容及泄洪流量，对一定的来水过程线，排水构筑物愈小，所需的调洪库容就愈大，坝也就愈高。尾矿库设计中常常需要通过几种不同尺寸排水系统的调洪演算，来合理地确定坝高及排水构筑物的尺寸，从而使尾矿库在满足使用功能的条件下达到整个工程造价最小的目的。 　　2.调洪演算的数学化 　　尾矿库调洪演算的基本思路是：在已知尾矿库的标高-库容关系曲线和洪水流量过程线等前提下，针对预先制定的几种调洪方案（包括泄流条件、调洪方式、排水构筑物参数等），由该尾矿库的调洪初始水位开始依次计算，直到洪水过程结束，分别给出各时段末排洪系统的泄流量、尾矿库蓄洪量等，由这些参数最终确定尾矿库调洪库容和设计坝高。 　　3.计算方法及其描述 　　从尾矿库洪水入泄关系和水量平衡方程出发，将数值计算方法应用于尾矿库调洪演算，采用拉格朗日插值法将尾矿库库容曲线及洪水入泄曲线解析化，并采用精度较高的四阶龙格-库塔法求解水量平衡微分方程与入泄关系，实现尾矿库库容的调洪演算。 　　3.1尾矿库入泄关系和水量平衡计算 　　3.1.1尾矿库排洪建筑物的泄洪能力 　　目前，通常选用排水井-涵管（或隧洞）或排水斜槽-涵管（或隧洞）作为尾矿库的排洪系统。排洪系统的泄流能力随尾矿库内水位的变化而变化：当水位较低时，排水井（或排水斜槽）的进水能力小于排水涵管无压满流时的最大泄流能力；反之亦然。前者以排水井（或排水斜槽）的进水能力作为排水系统的泄流量，后者以排水涵管的泄流能力作为排水系统的泄流量。 　　（1）排水涵管泄流能力 　　（1） 　　式中：为排水涵管的泄流量，m3；为过水断面，m2；为水力半径，m；为排水涵管的敷设坡度；为谢才系数。 　　a）排水井的泄流能力 　　当时， （2） 　　当时， （3） 　　式中：为井顶溢流的泄流量，m3/s；为井顶上的水头压力，m；为排水井直径，m；为重力加速度，9.8m/s2；为井顶溢流排水井的流量系数，采用0.51； 为流量系数，随的值而定。 　　b）窗口式排水井 　　（4） 　　式中：为水面以下进水窗口的面积和，m2；为水面至被淹没窗口中心的距离，m；为流量系数，矩形和圆形窗口取0.6。 　　（2）排水斜槽的泄流能力 　　（5） 　　式中：为斜槽断面面积，m2；为水面至斜槽中心的距离，m；为流量系数，按水流非全部收缩计算。 　　3.1.2尾矿库水量平衡方程 　　有限差形式的水量平衡方程，即在某一时段内（），入库流量与尾矿库下泄流量之差应等于该时刻内尾矿库库内水量的变化值，水量平衡方程的有限差分形式可表示为： 　　（6） 　　式中：、分别为时段初的入、出库流量；、分别为、时段初的尾矿库库内水量；为时刻长。 　　3.2尾矿库调洪演算的数值方法 　　3.2.1数值法演算的基本原理 　　在采用数值方法进行尾矿库调洪演算时，通常采用尾矿库水量平衡方程的微分形式。与式（6）所表示的尾矿库水量平衡方程的有限差分形式相比，尾矿库水量平衡微分方程主要从尾矿库水位与库面积或库容的动态变化关系出发，来描述某一时段内尾矿库积水量的变化情况。尾矿库水量平衡微分方程如式（7）所示。 　　（7） 　　式中：为尾矿库库容，是尾矿库水位Z的函数；t为时间；为入库流量；为出库流量，是尾矿库库容的函数；、和分别为尾矿库的起调水位、起调库容和起调时刻。 　　3.2.2尾矿库水量平衡微分方程入泄关系的解析法 　　基于尾矿库的泄洪能力，依照尾矿库入泄关系曲线和尾矿库水量平衡微分方程，可求解尾矿库各时段的泄洪量过程线及库容、库水位的变化过程线。尾矿库水量平衡微分方程与入泄关系的数值解析法求解就是对式（7）的一阶常微分方程初值求解，可采用四阶龙格-库塔（Runge-Kutta）法。四阶龙格-库塔法是4级计算方法，求解精度较高、简便，相应的尾矿库库容Vn，各级系数K1，K2，K3和K4如式（8）所示。 　　公式（8） 　　式中， 为插值计算时的时间步长；为四阶龙格-库塔法中每完成一步需四次计算的函数值； 、分别为时刻的入库流量和出库流量，其他符号意义同前。 　　3.3尾矿库特性曲线的处理方法 　　由尾矿库调洪演算可知，尾矿库的水位-库容曲线和库容-泄流曲线等尾矿库特性曲线通常都是单一的曲线，尾矿库库容-泄流曲线、水位-库容曲线等特性曲线，一般是以离散的列表函数形式给出。目前，从给定的尾矿库水位-库