第5次课矿井通风网络优化

1、2. 3应用计算机解算通风网络M的：已知风网各分支风阻和主通风机的特性，求算主要通风机的工况点，各分 支的风量和风向，以便验算各用风地点的风量和风整速是否符合规程要求。原P 依据风量平衡定律、风压平衡定律、阻力定律方法：回路法假设风网中每一回路内各分支的风向和风量开始，逐渐修正风量，使之满足风压平衡定律3节点法假设风网中每一回路内各分支节点压力值开始，逐渐修正压力分布， 使之满足风量平衡定律。一、改进的斯考德-恒斯雷试算法一一回路法把风流在风网中的流动看成是在一些互不重复的独立的闭合回路中各有一定的风量在循环，这种风量称为回路风量。如图：回路:ABDEF （风BCDB（只属于一个回路的分支。反

2、之，为非独立分支。且满足：独立分支伽）=分支总数一节点数（J）+ 1如：BCv CEvEFAB-独立分支，BD、DE、CD-非独立分支华北科技学院2.3应用计算机解算通风网络基本思；：初拟风网中各回路风量9iD（如q1q2。3），使其满足风网中节 点风量风量平衡定律，然后利用 风压平衡定律对其逐一进行修正，从而将各分支假设风量，经把迭 代计算修正，各回路风压逐渐趋 于平衡，这样各分支风量逐渐接真实值.回路风量修正值（）: Aq回路中各分支阻力代数和，当分支流向与回路流向一致时，取 “ + ”，反之，取。当回路中有和总时：想二落、2郭2 卜_ _ _故分支风量为：Q： = Qi 土 AQ华北科技

3、学院（二）通风网络解算数学模型由（5-5-1）式,设F（2HQyl diagCTQy-CHf-HN=Q <5-5-4）由于上式为非线性方程组，只能采用迭代法进行求解。设第k次迭代后的风量近似值为；将（5-5-4）式用泰勒级数展开并忽略二次 以上微分项得F（e；+，） = F（Q：） + -1。AQ： = o（5-5-5）g北科技学应故&=一（和）。：出尸（。；）（5-5-6）式中，&a为独立回路风量的修正值列向量；z第k+1次迭代后的风量可由下式计算：aM =a + AQ：5-5-9）上述式（556）和式（559） Bp为牛顿法解算通风网 的迭代计算模型。若雅可比矩阵具行

4、对角线优势，即则可略去雅可比矩阵的非对角元素，（556）式可简化为：式（5511）即为克劳斯（H, cross）近似试算法的计算公式, 也有的称其为斯考德-恒斯雷（Scott-Hinsley）法。（1）绘制通风网路图，标定风流方自。（2）输入网络结构及数据。（3）选择独立回路，确定独立结构的分支构成。（4）拟定初始风量。通常，先给余树边赋值，再计算各树支初 始风量。拟定的初始风量应尽量接近真实风量，以加快计算速度。（5）计算修正风量。使用牛顿法时由（5-56）式计算出各独立回路的风量修 正值，然后用（559）式对独立分支的风量进行修正，再用（5 -5- 2）式计算所有分支的新风量。使用克劳斯法

5、时，可用(5-5-11)分别计算各回路是风量修正值。当计算出一个回路的风量修正值后，立即对该回路所有分支的风量进行 修正，以加快收敛速度。(6)计算修正风量每计算完一次所有分支的修正风量，称为迭代一次。每次迭代后应 判断计算结果是否已满足给定的精度要求，即检查是否满足。max I Aqyi <s,<i<N(5-512)若上式巳满足，迭代计算终止。否则，转第5步继续迭代。精 度，根据需要，人为给定，华北科技学院克劳斯法的回路风量修正值可逐个回路独立计算， 因而可手工进行计算。而牛顿的通风量修正值的计算必须 先求雅可比矩阵的逆矩阵，实质上是线形方程的求解，一 般不能手工计算。由于

6、克劳撕法实质上是牛顿法在满足 (5-5-10)式条件下的一种简化，因此在选择独立回路的 分支结构时有特殊的要求。上述迭代计算过程容易编写成计算机程序，目前已多 种成熟的通风网络解算软件。例：某简单通风网络如图5-5-1所示。网络分 支数n=6,节点数m=4,分支风阻向量为区=（0.01, 0.38, 0.5, 0.65, 0.20, 0.085）。设在分支1中有风机，风机风压特性方程：hf= -0.3433 +16.2027% 61.2432图 5-5-1风阻、风量和风压的单位均采用国际单位制。不考虑 自然风压，式计算通风网络中各分支的风量。解：“克劳撕”法。由（5-5-11）式可得风量的修正

7、 计算式为：华北科技学院将已知风阻和初始风量代人上式得:=2&MI +2&|密一9+2.一如I + 0. 686691 - 16- 2027& W夕2 + &3一船I如一姨一凤© 一九I 一 3zETJil + 2凡 1他一。$I + 2凡 <?1 <h将已知风阻和初始风量代人上式得:当 ql<n 二舄际 一 q"（外 一 qg）二处出 _ qj'z二城2&I或I + 2七皿-如1 + 2&也一 9sI用上述修正风量值对应余树分支风量进行修正后得到:重复上述迭代过程，经10次迭代后可使修正风.量的最大

8、 值小于0. 001I113/S；最后得到与上述“牛顿法”完全一致的 计算结果。（三）风量调节计算数学模型L调节点的数量与位置华北科技学旄由前面节点风量平衡方程（552）式可见，只要网 络中的条独立分支的风量已确定，则其他分支的风量也就被 确定。假设在每'一条独立分支中安设一个调节设施，就可以 实现对风网风量的完全控制。因此通风网络中的最多有效调 节点数量为N。若网络仅靠自然分风就能满足全部风量要求，则网络 内部不需节。但每个矿井罕少应有一台主通风机提供通风动 力，而通风机可视为一个调节设施，因此通风网络的最少调 节点数为1。在实际通风网络中，常有一些不需调节的自然分风子 网络，设自然

9、分风子网窄包含k个独立回路，则网络中应具 有的独立调节点数目为：S = N k = m k + 1（5513）虽然在通风网络中仅需S个调节点，且这S个调节点必须包含在网络的同一棵余树中。但由于通风网 络的生成树数目非常之多，因此实际上网络中任何 一条分支都有可能成为调节点。如何在通风网络的全部',条分支中选出所 需的S个调节分支，就是确定风量调节方 案的重要内容之一。选择调节点位.皆应遵守的总的原则是：安全可靠、技 术可行和经济合理。为保证所选择的调节点安全可靠和技 术可行，一般应注意以下问题：(1)在矿井总进风巷道和各类运输频繁的井巷中，应避免 安装调节设施，如调打风窗、带风墙的辅助

10、通风机等；(2)在矿井总回风巷道、在倾斜和垂直井巷中，不宜安装 增阻调节设施；(3)在采区内部服务年限短的巷道中，不宜采用工程量很 大的减阻调节措施；(4)在封闭火区附近采用调节措施时，要考虑有利于火区 的均压灭火。华北科技学院“回路法”是事先人为地确定需进行调节的分支，将令:Hr = Rdias IQ dias Q = R飒 1 CTQy I diagCTQy显然这种方法求出的调节方案的优劣完全取决于调节分 支的选择。适用于调节点位置固定后计算其调节量的大小。3 .固定风量法这是一种在自然分风程序中加固定风量分支，以求解按 需分风问题的方法。其实质就是“回路法”。其算法原理是: 求解时将固定

11、风量分支选取为独立分支（余树枝），初拟风量 时先将固定风量加到其所在回路的所有分支中；此后在修正。华北科技学院各回路风量的迭代计算时把它从网络中分裂开，不进 行迭代计算；在自然分风计算结果后，将固定分量分支所 在回路的风压闭合差作为该分支的阻力调节量。前一部分 的自然分风计算，实际玉就是计算自然分风弋M络66风不 门然分配；后一部分计算调节量时，是假设需进行调节的 分支就是需风分支，因而根据回路风压平衡方程即可计算 出调节量。自固定风量法的算法原理可发现这种方法的两个主要问题：（1）为保证除固定风量分支以外的自然分风子网络是一 个连通的闭合网络，就不允许穿过用风区的割集全部由固 定风区的割集全

12、部由固定风量分支组成，而必须至少有一 条自然分风分支。但实际矿井在正常生产时用风区都必须 按需分风，因而计算时不得不反复试算，以满足放松为自 然分风分支的需风分支的需风量。（2）只能在有需风要求的分支本身进行调节，很难保证 调节方案的可行性。华北科技学院4 .通路法“通路法”是通过构造由网络总进风节点（源点）至总回风节点（汇点） 的所有独立风路（通路、道路），由各通路的风压差来计算调方量。若 对各风机子系统分别选取最大风压的通路作为计算标准，则在其他通路 上必须采取增阻调节的方法，因而有时也称之为“最大阻力路线法”。 当然也可以根据风机能力选择合适风压的通路为基准通路，风压大于基 准通路者应减

13、阻，小于基准通路者应增阻。由于每一条独立通路可以加上一条连结挥点与源点的地面伪分支， 构成个独立问路，因此调节分支的选择与问路相同。通路法也必须事 先选择好调节分支，然后将调节分支选为独立分支，包含在同一棵余树 中。每一条通路中仅包含一条独立分支，各通路风压与基准通路的风压 差就是通路独立分支的调节量。通路法虽然可保证风机压最小，但也不能保证调节点位置的合理性。 目前有些通路法程序仅能得到一种调节方案；有西通路法程序则能自动 得到全部可能的调节方案；还有的通路法程序在选择调节分支时，尽量 使一个调节能调平较多的通路，从而可得到调节设施数最少的调节方案。五、通风网络解算程序的使用说明(1)输入仃

14、关参数通常需输入分支的始末节点、风阻、风机风压和固定风量； 若网络风阻值未知，则还须输入管道断面几何形状特征尺寸、摩 擦阻力系数、局部阻力系数、管道长度等参数，以便计算风阻值。 (2)各分支风阻按降序排列风网各分支参数的输入先后次序是任意的，程序按其风阻大 小以降序排列。(3)确定余树弦。华北科技学院余树弦是做为选择回路的基底分支，除固定风量分支和风机分 支做为法定的余数弦外，还要从一般分支中选择风阻较大者作为余 树弦，使余树弦的总数为(n-m+1)个.这类分支构成回路的非公共分 支，每个回路只有一个，并且其风阻值在回路的各分支中最大，这 类分支数与独立回路数相对应.选余树弦的方法是在网络中形

15、成最 小树，而未包含在最小生成树中的分支即为余树弦.本程序采用 “避图法”，寻找最小树。即把各分支的风阻值按升序排列，然后 从最小风阻分支开始，按风阻升序逐个取出各分支添加到原图的位 置上，检查该分支，如果与前面添加的分支构成回路，则将此从图 中移出，成为余树弦；否则将此填入树中，成为树枝.如此继续下 去，直到找出(m-l)条分支所构成的连通因为止，此图即为所要求的 一棵最小树.同时也找出了(n-m+1)个余树弦.华北科技学院(4)选回路选回路是以余树弦为基础，从一般分支中风阻最大的分支开始，逐一查找可与余树弦出风节点相联的分支，反复操作直到构成回路为止°程序中规定余树弦的方向为回路

16、的方向，回路中与余树弦方向一致的分支为正向分支，否则为反向分支，并逐一记录。程序中只选择一般分支与余树弦相联，从而保证了 一个回路只含一个余树弦。华北科技学班图2-3-1通风网络解算程序框图华北科技学班通风机特性曲线的拟合根据通风工程解题精度的要求，通常采用二次或三次多 项式拟合风机曲线。采用二次多项式拟合风机风压曲线，应 在风机特性曲线上选三个点，即风机稳定工作的上限点、高 效点和下限点.三个点的效率均应大于70%,根据三点的风 量，风压座标值(，)可求出拟合系数a, b和c按式(2-3-5), 只要确定出风机风量即可算出风机的风压值。(6)风量赋初值在迭代计算之前应将风网的余树弦赋风量初值

17、，然后根 据风量平衡方程对每一回路中的各分值赋风量初值，作为迭 代计算的开始。程序规定：固定风量分支以固定风量值为初 值，风机分支以风机特性曲线第二点的风量为初值；其它余 树弦均以0. 5为初值。然后按回路对各个分支赋初值。华北科技学院(7 )迭代计算迭代计算是在风网各分支风阻和风机性能曲线已知的情况下，对 风网中的每个回路的余树弦风量进行迭代修正计算，为保证各用风地点 的固定风量不变，程序规定有固定风量的不参与迭代计算，由于这类回 路中除固定风量分支外，其余分支都与其他回路公用的分支，因而实际 上都参与了迭代计算。迭代计算以回路为单位，直到所有回路的风量修 正值都达到预定精度为止。如在给定的

18、最大迭代次数内未达预定的精度， 为加快收敛，程序可将各分支的(RQ)值赋给R,然后返回，重新排序 确定余树弦和选回路，再做迭代计算。(8)计算各分支阻力并打印计算结果按已算出的风量和输入的风阻值计算各风道的阻力，同时按风量 和风道断面计算风速，最后以表格形式打印各分支号、始末节点号、风 阻、风量、阻力、断面和风速。 (9)计算固定风量分支的阻力值由于固定风量分支不参与迭代计算，固定风量分支的 阻力是按回路风压平衡式(2-3-15)计算，然后用阻力定律 算出固定风量分支的风阻值。该风阻值是保证固定风量值 不变，所必须具有的风阻值。当该值小于实际的风阻值时, 需进行减阻调节，反之需进行增阻调节。斗

19、(10)处理风机参数根据风机曲线拟合系数和风机风量，计算风机风压， 然后打印风机风量、风压。如果风机工作点落在合理工作 范围之外，则打印超限警告，必要时可改变风机转速或更 换风机，重新进行通风网络调节计算。华北科技学医(11)上机操作注意事项：将上述数据文件愉入计算机并存入硬盘，如果编写和输入的数据文 件有误，计算机将提示出错的类型或得出错误的结果。因此，上机前一 定要反复核对，如果运算中出现错误也要按上述要求检查核对。一般上 机过程中出现的错误有以下几种：通风网络分支数、固定风量个数和风机个数与实际输入的不符；分支节点号出错，将导致通风网络结构发生变化。这种错误较为隐 蔽，不易被发现，只有当

20、回路选不出，程序中断运行时，方可发现。因 此，必须认真检验分支节点号的正确性，避免此类错误的发生。风机分支号出错，这种错误从计算机输出的数据栏目中较易发现。总之，无论出现何种错误，只要计算机不能正常运算或出现非预想的 结果，都要重新审查数据，确认输入数据无误后方可进行通风网络解算。华北科技学庞2. 4 通风管网系统优化设计随着工业的发展，通风除尘、通风空调系统具有大型化、集中化的 趋势。大型化的复杂工业通风除尘和通风空涧系统存在两个突出问题， 即能耗问题和风量不平衡问题。70年代发展起来的工业通风和通风空调 系统优化设计法，对降低系统能耗，减少系统总费用方面起到了积极的 作用.传统的设计方法，

21、一般把求解费用最/氏问题和解决风平衡问题作 为两个独立的问题进行处理，即首先根据目标函数求解最优参数，然后 进行风量平衡调节。这样的优化设计方法，在减少费用方面取得了一定 的效果，但存在风量平衡精度低的缺点，尤其对日益复杂的通风管网系 统的设计，风量不平衡问题更为突出。因此：把费用最低目标函数和系 统风量平衡方程组成一个非线性的数学规划问题，使设计优化参数能同 时满足目标函数和系统风量平衡条件，设计中不需进行风量平衡调节的 计算。华北科技学院一、通风管网系统优化设计的数学模型L目标函数目前，一般以系统使用年限(即系统寿命周期)内的总费用最低作 为目标，如果不考虑初期投资的贷款利率和货币的时间价

22、值，总费 用G可表示为：式中：G =GaSGQ (2-4-1)Ga系统的初期投资，元；Gh系统的年运行与维护费用，元/年；Gc一系统的年废料回收费用，元/年；t系统的使用年限，年。系统使用年限的取值：排含腐蚀性气体的系统t=5年；一般除 尘系统t = 10年；一般通风系统t=20年。系统初期投资Ga乞括管道费用Ga】、风机费用Ga2和除尘器费 用勾3；系统运行与维修费用Gb包括动力消耗费用Ga和维修费用 Gb2根据其是否随“变量”变化，可分为可变费用和不变费用， 随“变量”变化而变化的费用称为可变费用，反之则为不变费 用。对于通风除尘系统的优化设计，系统费用的变量是风管直 径。显然，Gal.

23、Ga2, Gb均为可变变费用；Ga3和Gc为不变费用。 所以，系统寿命周期总费用中的可变费用Gv为：华北科技学院G、，= 0+Ga2 , Gh（2-4-2）式中各项按以下各式确定：G0=Z4 ,勿。L（243）Gb = V a-b-1（245）PI<式中 D.第i分支管道的管径，m；4 第i分支管道的长度，m；A 园管单位面积的价格，元/N.第i分支管道分担的风机功率，kW,风机 总功率0 出风机的价格系数，由风机产品说明书统 计求得，见表2-3-1；a 电价，元/（左卬/）b-风机的年运行小时，ho华北科技学院表2-3-1 4-72 I型风机价格系数传动方式 ABC40015004400606895每分支管道分担的风机功率N,为:式中：Qi第i分支管道的风量，加3s风机的总效率，一般取 =0.60.7； 第i分支管道的通风阻力，尸a华北科技学院按阻力计算公式有：区=(2今 + 0)哼I(2-4-7)式中：2风管的摩擦阻力