复杂矿井通风网络的分支风流计算方法.doc

复杂矿井通风网络的分支风流计算方法 刘泽功 摘要： 目前广泛应用的解算矿井复杂通风网路方法是以风网的各分支风阻值和扇风机工作特性参数作为解算基础数据，而获得风网各分支的风阻力值是费时费力的。为此，本文提出了通过测算风网风量求分支风阻的方法解算矿井复杂通风网路。从理论上分析了以分支风阻为未知数的线性方程组建立过程及方程组解存在的唯一性。用BASIC语言编写了测风求阻法解算复杂通风网路的计算机程序，并在386微机上通过。通过具体实例说明程序的实用性。该解算方法与其它解算方法不同的是本法以风网的风量和扇风机工况点作为解算网路的基本依据，对生产矿井的通风网路解算更能显示其优越性。 关键词：通风阻力，测量，计算，电子计算机 概述 目前，比较广泛使用的解算复杂通风网路方法的关键问题之一是对通风网路图G=（V, E）（其中V是节点的集合、V=m, E是边的集合、E=n）中的n条分支的风阻值R（固定风量分支除外）都必须确知，才能获得正确的解算结果。对于大型的矿井通风网络，分支少则几十条、多则几百几千条，本文提出了一种用分支风量和扇风机工况点作为基础数据解算矿井复杂通风网路的新方法—测风求阻法。并编写了用计算机解算网路的解算程序。 基本定律与基本原理 基本定律 （1）通风阻力定律 根据通风学原理，通风阻力定律为： （1） 该式中，h代表的是巷道通风阻力，Pa；R代表的是分支风阻，kg/m7；Q代表的是风量，m3/s。 （2）节点风量平衡定律 下面的公式对于通风网络的任意节点均是适用的： （2） 该式中，i代表的是网络中的节点编号（i=1，2，3，……，m）；j代表的是网络中分支的编号（j=1，2，3，……，n）；Qj代表的是第j条分支的风量，m3/s；aij代表的是与第i号节点相连接的第j条分支的风向函数： 1 第j条分支风流流入节点i aij= 0 第j条分支不与节点i相连 -1 第j条分支风流流出节点i （3）网孔风压平衡定律 （3） 该式中，i代表回路的编号（i=1，2，3，……，n-m+1）；j代表的是网孔的分支编号（j=1，2，3，……，n）；HFi代表的是第i回路中风机工作风压，Pa；HNi代表的是第i回路的自然风压，Pa；Rj代表的是分支j的风阻值，kg/m7；Qij代表的是通过分支j的风量，m3/s；bij代表的是表示第i条回路中第j条分支的风向函数： 1 回路i中第j条分支风向与回路方向一致 bij= 0 第j条分支不属于回路i -1 回路i中第j条分支风向与回路方向相反 1.2 基本原理 以通风网路的有向连通图G0=（V0，E0）（│V0│=m，│E0│=n）作为基础网路图，把各分支的已知风量Qi（0）（i=1，2，3，……，n）作为基础数据，圈划独立回路，然后根据（3）式建立以分支风阻为未知数的线性方程组[2]。由图论理论可知，对于有向连通图G0的独立回路个数为r0(r0=:n-m+1)[1,2]，因此，仅由图G0不能建立起n个独立回路方程。但是，生产矿井经常对通风系统进行调整，有时为了解算通风网路，在不影响生产的情况下，也可以而且必需人为地对通风系统进行调整。以G0作为基本网路图，通风系统进行第一次调整时，则G0变为G1=（V1，E1）（│V1│=m，│E1│=n）此时各分支风量变为Qi(1)（i=1，2，3，……，n1），独立回路方程个数r1（r1=n1-m1+1）；第二次调整后，网络变为G2=（V2，E2）（│V2│=m，│E2│=n），分支风量为Qi(2) （i=1，2，3，……，n2），独立回路方程个数r2（r2=n2-m2+1），……，进行k次调整时，网络变为Gk=（Vk，Ek）（│Vk│=m，│Ek│=n），Qi(k)（i=1，2，3，……，nk），独立回路方程个数rk（rk=nk-mk+1）。对于基本网络图G0，r1，r2，r3，……rk回路方程中，对应有r1’，r2’，r3’，……，rk’共k个回路方程含有被调阻分支，则在建立分支风阻回路方程组时，应首先将r’（r’=r1’+r2’+r3’+……rk’）个方程从r（r=r1+r2+r3+……rk）中剔除。对k的取值原则是保证r-r’n，一般情况下：只要对网路图G0进行1~3次调风，建立起来的独立回路方程数就会超过未知数个数n。 2 数学模型的建立 2.