Research on calculation and optimization method of refrigeration and cooling pipe network in high temperature mine
-
摘要:
高温矿井采用井下集中式制冷降温系统受采掘范围的影响,冷冻水输送管网复杂,导致末端冷量供给不足,严重影响系统降温效果。为提高矿井降温系统冷量利用率,以赵楼煤矿井下制冷降温管网为例,基于图论原理建立井下降温管网拓扑模型,采用水力基本方程计算管段流量、节点阻力和水力损失,得到降温管网水力特性;通过管网节点温升的计算,确定冷冻水输送过程的冷量损失;结合管网水力、热力特性,对降温管网进行优化。结果表明:赵楼煤矿井下制冷降温管网末端工作面冷冻水流量最小为0.001 m3/s,管网水力损失大的位置为一集轨道下山、7302运输巷和中部辅运大巷。降温系统分别给五采区和七采区共4个工作面供冷,冷量损失为1.22×106 J/s,其中七采区降温系统冷量损失占88.15%,管网摩擦和传热冷量损失分别为5.39×105、6.805×105 J/s,末端冷冻水最高输水温度为13.9 ℃。提出管道−泵阀联调优化方法,采用动态平衡阀对南部1号辅助运输大巷和二集辅助巷冷冻水流量恒定在0.022~0.04 m3/s,实现末端空冷器流量稳定;采用静态平衡阀调节管网支路阻力,将工作面冷冻水流量增大至0.005 m3/s;南部1号辅助运输大巷、南部2号辅助运输大巷和二集辅助巷的管径增加至0.325 m,管网总水力损失从30.93 m减少到20.44 m,减少了35%;将离心泵扬程调整为183~195 m、流量为0.085~0.112 m3/s,离心泵高效运行。
Abstract:Adoption of centralized underground refrigeration cooling system in high-temperature mine is affected by the range of mining. The complex chilled water transmission pipeline network leads to inadequate supply of terminal cooling capacity, which significantly affects the cooling effectiveness. In order to improve the utilization of cooling capacity in mine cooling system, taken the underground refrigeration cooling pipe network of Zhaolou Coal Mine as an example. The topological model of cooling pipe network was established based on the principle of graph theory. The hydraulic basic equation was used to calculate the flow, nodal resistance and hydraulic losses of the pipe section to obtain the hydraulic characteristics of cooling pipe network. The cooling loss of chilled water transportation was determined by calculating the temperature rise of pipe network nodes. The cooling pipe network was optimized by combining the hydraulic and thermal characteristics of pipe network. The results shown that, the minimum chilled water flow was 0.001 m3/s in the working face at the end of underground cooling pipe network in Zhaolou Coal Mine. The locations with significant hydraulic losses in the pipe network were A set of Track Downhill, 7302 Transport Roadway and Central Auxiliary Transport Roadway. The system respectively supplied cooling to four working faces in the Fifth and Seventh mining areas. The total cooling loss was 1.219×106 J/s, of which the cooling system in Seventh mining area accounted for 88.15%, the friction and heat transfer losses of pipe network is 5.39×105 J/s and 6.805×105 J/s, respectively, and the maximum temperature of chilled water at the end was 13.9 °C. The optimization method of pipeline-pump-valve inter-conditioning was proposed, the dynamic balance valve was used to maintain the chilled water flow in the range of 0.022−0.04 m3/s in the South 1 # Auxiliary Transportation Roadway and the Second Set Auxiliary Roadway, to achieve the flow stability of the end air cooler. The static balance valve was used to adjust the branch resistance of the pipe network, and the chilled water flow was increased to 0.005 m3/s at the working face. The pipe diameter of the South No.1 Auxiliary Transportation Roadway, the South No.2 Auxiliary Transportation Roadway and the Second Set Auxiliary Roadway was increased to 0.325 m, and total hydraulic loss of the pipe network was reduced from 30.93 m to 20.44 m, which was a reduction of 35%. The centrifugal pump head was adjusted to 183−195 m with a flow rate of 0.085−0.112 m3/s, which ensured efficient operation of the centrifugal pump.
-
Keywords:
- mine cooling /
- refrigeration cooling pipe network /
- hydraulic loss /
- cold loss
-
0. 引 言
随着煤矿开采深度的增加[1],矿井热害问题逐渐严重[2]。矿井常采用井下集中制冷降温系统,将制冷机组的冷冻水通过管网输送至降温地点空冷器进行降温[3]。制冷降温系统中冷冻水在输送过程中受管道布置、管径和长度等因素的影响,存在末端冷量供给不足[4],制冷效率低[5]等问题。
国内外学者分别对管网水力特性、压力损失和系统优化进行研究。魏立杰等[6]、郑进福等[7]基于节点分析法分析各管段阻力分布。RAHIMI等[8]引入基图和余数图,通过计算管道流量和水头损耗来确定管道直径。窦明等[9]基于图论原理定量研究水量平衡关系,得到不同工况下最佳水系连通方案。SABZKOUHI等[10]分析了管网节点需求、管道摩擦系数、泵扬程与管网水力特性的关系。CHICHERIN等[11]研究了管道中水的压力、密度、速度等对供热管网的阻力影响,得出改变管道粗糙度可显著减小水力损失。CHANG等[12]利用海曾威廉系数和水力损失值来寻找敏感度高的管段。当管网流量供给不平衡会造成区域冷热不均,离心泵与管网不匹配会造成泵可靠性降低和故障多发等问题。由世俊等[13]分析了阀门开度和水泵转速对管网动态响应时间的影响。HONG等[14]提出将单个泵特性曲线与管网水力特性曲线联合优化管网水头和流量。高铁军等[15]建立了阀门−管网拓扑模型,评估和识别关键阀门,提高了阀门可靠性。周璇等[16]动态优化了供热系统热水循环泵的运行频率和各支路电动调节阀的开度降低系统运行能耗。
降温系统冷冻水管路随着采掘范围增加距离变长,管网复杂造成末端需冷区域冷量供给不足。笔者以赵楼煤矿井下制冷降温管网为例,基于图论原理建立管网拓扑模型,计算节点流量和阻力,得到管网水力和热力特性,提出降温管网优化方法。
1. 制冷降温管网解算
1.1 水力计算方法
冷冻水通过制冷管网由制冷机组输送到各末端需冷地点,管网输送阻力损失主要包括沿程阻力损失和局部阻力损失2部分[17]。
沿程阻力损失按式(1)计算:
$$ {h_{\mathrm{f}}} = \lambda \frac{L}{D}\frac{{{v^2}}}{{2g}} $$ (1) 式中:$ {h_{\mathrm{f}}} $为沿程阻力损失,m;$ \lambda $为沿程阻力系数;$ L $为管长,m;D为管道内径,m;$ v $为流速,m/s;$ g $为重力加速度,取9.8 m/s2。
局部阻力损失主要由管道中的阀门、弯管等组成。产生局部阻力的原因是流体与接触边界相互碰撞产生漩涡。局部损失按式(2)计算:
$$ {h_{\mathrm{m}}} = \xi \frac{{{v^2}}}{{2g}} $$ (2) 式中:$ {h_{\mathrm{m}}} $为局部阻力损失,m;$ \xi $为局部阻力系数。
则总阻力:
$$ h = S{Q^2} $$ (3) 式中:$ h $为管路总阻力,$ h = {h_{\mathrm{f}}} + {h_{\mathrm{m}}} $,m;$ S $为管路阻抗,$ S = 8\left(\lambda \dfrac{L}{D} + \xi \right){({{\text{π}} ^2}{D^4}g)^{ - 1}} $;$ Q $为流量,m3/s。
管路中总阻力、管路阻抗和流量是管网水力计算中的重要参数,也是水力平衡性的重要指标。管路阻抗大小取决于管道长度、管径、局部阻力系数等[18],但在管路设计完成后,管路阻抗是固定值,管路总阻力随流量变化而变化。
1.2 热力损失计算方法
管网中冷冻水流动时冷量损失主要包括传热损失和水阻摩擦损失2部分[19]。传热损失指冷冻水与周围环境进行对流换热损失的冷量,取决于管道热阻以及管道内冷冻水温度与周围环境温度的温差。管网包括供水管道与回水管道,供水温度低于回水温度,而供水管的冷量损失决定了末端冷量供给。水阻摩擦损失指冷冻水输送过程中摩擦生热导致冷冻水温度升高,最终以热的形式耗散在冷冻水中引起温度升高。
管网传热温升:
$$ \Delta {t_{\mathrm{f}}} = AL $$ (4) 式中:$ \Delta {t_{\mathrm{f}}} $为传热温升,℃;$ A $为单位管道温升,赵楼煤矿制冷管道为聚氨酯发泡保温无缝钢管,经计算,单位管道温升取值为0.001 ℃/m。
水阻摩擦温升由焦耳−汤姆逊效应[20]计算:
$$ \Delta {t_{\mathrm{w}}} = \mu L\Delta p $$ (5) 式中:$ \Delta {t_{\mathrm{w}}} $为水阻摩擦温升,℃;$ \Delta p $为单位长度管道的阻力损失;$ \mu $为摩擦升温系数,取2.4×10−7 ℃/Pa。
对于管道长度$ L $的管段单位面积总冷量损失:
$$ q = \frac{D}{4}v\rho {C_{\mathrm{P}}}\frac{{\Delta t}}{L} $$ (6) 式中:$ q $为传热冷量损失,kJ/(s·m2);$ \rho $为水密度,取
1000 kg/m3;$ \Delta t $为管道总温升,℃;$ {C_{\mathrm{P}}} $为水的定压比热容,取4.28 kJ/(kg·℃)。总冷量损失为
$$ {Q_{\mathrm{r}}} = q\sum\limits_{i = 1}^n {\pi {D_i}{L_i}} $$ (7) 式中:$ {Q_{\mathrm{r}}} $为总冷量损失,kJ/s。
由式中可看出,决定管网冷量损失的参数包括管长、管径、流速和管段温升。
2. 高温矿井制冷管网优化分析
2.1 赵楼井下制冷管网概况
赵楼煤矿采用井下集中式制冷降温系统,冷冻水管网布置如图1所示。井下有东部和中部2个制冷硐室,各有1台KM3000型制冷机运行,制冷功率为6.6 MW,单台制冷机可提供流量为0.04 m3/s的3 ℃冷冻水。东部制冷硐室为七采区供冷,中部制冷硐室为五采区供冷。矿井目前安装制冷管路约40 000 m,采用聚氨酯发泡保温无缝钢管,管路连接处使用保温海棉进行处理。
井下集中制冷降温系统原理如图2所示。制冷机组产生的冷冻水经输冷管路送至末端空冷器实现对风流降温,冷冻水经过空冷器后沿回水管返回至制冷机组,形成冷冻水循环。
降温系统建成初期,工作面降温幅度在6 ℃左右,降温效果良好。随着矿井开采范围扩大,需冷区域增大,造成制冷机组与热交换器之间的距离增加,输冷管道较长,沿程阻力增大,管路中含有各种产生局部损失的管路元件,局部阻力增加;水头压力变化和冷冻水沿程冷损增加,制冷效果差,需要优化管网系统。
2.2 管网拓扑模型
根据赵楼煤矿井下制冷降温管网布置图,建立了与之对应的拓扑模型,管网模型包括节点与管段。各末端支路依次编号,节点之间通过管段连接,节点为管段产生流动损失的位置,如图3所示。
为方便管网水力和热力特性辨识,对管网系统模型进行适当简化:①冷水机组及水泵等制冷机房设备合并作为管网的1个节点;②对于相邻且中间无节点的管道,均视为同一管道简化处理;③冷冻水系统是闭式循环,回水管和输水管对偶,管网节点既代表输水管节点也表示回水管节点。
七采区降温管网有1个冷冻水源,即东部制冷机组(E1),冷冻水流量0.04 m3/s。输水管网共有15个节点,14根管段。主要输水管为沿南部1号辅助运输大巷和二集辅助巷布置的管道。管网基本信息及设定见表1。
表 1 七采区降温管网信息Table 1. Information of cooling pipe network in Seventh mining area管段 管段位置 管长
L/m管径
D/m标高/m 流量
Q/(m3·s−1)Si Si+1 E1-2 南部1号辅助
运输巷178.8 0.377 −858.3 −857.1 0.04 E2-3 南部1号辅助
运输巷810.7 0.377 −857.1 −855.6 0.023 E3-4 南部1号辅助
运输巷956.22 0.325 −855.6 −853.5 0.023 E4-5 南部1号辅助
运输巷1134.5 0.219 −853.5 −923.9 0.023 E5-6 南部1号辅助
运输巷341.02 0.219 −923.9 −853.4 0.012 E5-7 二集辅助巷 168.52 0.219 −923.9 −923.3 0.012 E7-8 二集辅助巷 1481.3 0.273 −923.3 −927.1 0.012 E2-9 一集轨道
下山898.21 0.273 −857.1 −910.0 0.017 E8-10 7302运输巷 1596.5 0.159 −927.1 −792.0 0.006 E10-11 7302工作面 300 0.159 −792 −861.9 0.006 E8-12 二集辅助巷 325 0.159 −927.1 −924.9 0.006 E12-13 7302轨道
巷道1250 0.159 −924.9 −829.7 0.003 E12-14 二集辅助巷 370 0.159 −924.9 − 1015.0 0.003 E9-15 一集轨道
下山560.67 0.159 −910.0 −945.0 0.017 五采区降温管网有1个冷冻水源,即中部制冷机组(M1),冷冻水流量0.04 m3/s。输水管网共有15个节点,14根管段。管网基本信息设定见表2。
表 2 五采区降温管网信息Table 2. Information of cooling pipe network in Fifth mining area管段 管段位置 管长
L/m管径
D/m标高/m 流量
Q/(m3·s−1)Si Si+1 M1-2 南部2号辅助
运输大巷613 0.325 −848.9 −849.2 0.015 M2-3 南部2号辅助
运输大巷419 0.159 −849.2 −822.5 0.007 M3-4 5309运输巷 675 0.159 −822.5 −892.2 0.004 M3-5 南部2号辅助
运输大巷274 0.159 −822.5 −782.7 0.004 M5-6 5309轨道巷 693 0.159 −782.7 −865.1 0.004 M2-7 2号辅助6号
联络巷250 0.159 −849.2 −840 0.007 M7-8 5304轨道巷 993 0.159 −840 −808.9 0.003 M7-9 2号辅助6号
联络巷218 0.273 −840 −803.3 0.004 M9-10 5304运输巷 980 0.159 −803.3 −772.2 0.002 M9-11 2号辅助6号
联络巷207 0.273 −803.3 −778.4 0.002 M11-12 5305运输巷 1110 0.159 −778.4 −732.8 0.001 M11-13 2号辅助6号
联络巷700 0.159 −778.4 −761.2 0.001 M1-14 中部辅运
大巷855 0.159 −848.9 −900 0.010 M1-15 南部2号辅助
运输大巷1236 0.159 −848.9 −853.3 0.015 2.3 水力特性分析
2.3.1 流量分布
由图3a知,当东部制冷硐室(E1)的节点流量为0.04 m3/s时,七采区降温系统处于稳态水力变化中,冷冻水从东部制冷硐室流出沿南部1号辅助运输大巷在节点E2分别向管段E2-3-4-5和E2-9供水,流量为0.023和0.017 m3/s。在节点E5分别流向E5-6和E5-7-8,其中E5-7-8是布置在二集辅助巷流向工作面的支管,流量为0.012 m3/s。在节点E8处冷冻水流向7302工作面,其中管段E8-10、E8-12流量为0.006 m3/s,管段E12-13、E12-14流量是0.003 m3/s。流入工作面冷冻水流量小于0.008 m3/s,主要原因是主管E2-3-4-5管径逐渐缩小,同时在支路分流的作用下,远离制冷硐室的节点E10、E11、E13流量偏小;而E2-9管径较大,导致大量冷冻水流向一集轨道下山巷道。
五采区管网对5304、5305、5309工作面进行降温,如图3b所示,主管是布置在南部2号辅助运输大巷中的M1-2-3-5,供水管路支管较多。当中部制冷硐室(M1)节点流量为0.04 m3/s时,流入5305运输顺槽冷冻水流量为0.001 m3/s,3个工作面输水流量均小于0.005 m3/s,主要是管网支管多,串联布置时,靠近制冷硐室的支管流入的冷冻水流量较大。
2.3.2 管网阻力损失
七采区管网总长度达
20700 m,管段管径为0.159~0.377 m不等,如图4所示,沿程总损失为7 m,E2-3-4-5管段流量不变,管长分别为810.7、956.2、1134 m,管径由0.377 m缩小至0.219 m,冷冻水流速增加,沿程损失分别为0.10、0.24、1.96 m,千米管长损失为0.12、0.25、1.74 m。当冷冻水流量不变时,管径变化导致管道水流速变化时,千米管长损失也随之变化;而冷冻水在管道中流经的长度越长时,沿程损失越大。七采区局部损失共0.093 m,主要由弯管、三通、变径管产生。节点E2、E5、E8、E12是三通分流处,E3、E4为管段变径处,E15为弯管。局部损失最大节点为E15,损失0.036 m,原因是管径偏小,当冷冻水流量较大时,造成局部损失大。五采区管网总长度约
18400 m,沿程总损失22.7 m,如图5所示。沿程损失集中在M1-14、M1-15,分别为4.29、12.9 m;冷冻水在节点M1分别向管段M1-2、M1-14、M1-15供水,千米管长损失为1.97、5.02、10.43 m;其中M1-14、M1-15冷冻水流速较大,分别为0.18、0.26 m/s,说明当冷冻水流速越大、管径过小时,单位管长损失越大。五采区中产生局部损失的管道附件多,M1、M2、M3、M5、M7、M9、M11均为三通分流处,局部损失均小于0.001 m,原因是冷冻水流量小,流速较小,局部损失少。井下制冷降温管网中压降较大的管段有E5-6、E8-10、E10-11、E12-13、E12-14,分别为1.3×106、1.4×106、2×106、9.2×105、9.3×105 Pa,位于二集辅助巷和7302工作面。此管段两端节点高差分别为70.5、135.1、69.9、90.2、90.5 m,管路之间的高差造成冷冻水流经时所需能量大,管段压降增大。
2.4 热力特性分析
2.4.1 管网温升变化
赵楼煤矿降温系统每千米温升在2~4 ℃,其中七采区温度变化大,主要原因是七采区需冷区域与东部制冷硐室距离远、沿程损失大,末端温度较高,如图6、图7所示。
当制冷硐室冷冻水输/回水温度3 ℃/16 ℃时,七采区降温管网各节点温度变化较大,节点E11在7302工作面上,距离制冷硐室最远,输水长度
6700 m,冷冻水温度高达13.9 ℃;节点E13、E14温度大于13 ℃,布置在7302轨道巷、七采区辅助巷附近,输水长度大于5000 m,温升0.75 ℃/100 m;管道虽有保温处理,但沿程传热损失大,末端换热器输水温度高、降温效果不佳。五采区降温系统各节点温度变化较小,平均温升3 ℃;温度较高的节点M11、M8、M4、M10、M6、M13、M12在树状管网末梢,输水温度大于6 ℃,位于5304、5305、5309工作面,输水距离超过2000 m,沿程热损失大,温升0.12 ℃/100 m。2.4.2 冷量损失分析
管网中冷量损失大的管道集中在七采区,多个管道冷量损失大于105 J/s,其中管道E1-2-3-4-5为主管道,水流速度快,管道温升大,冷量损失大;管道E7-8、E2-9输送距离大于
1000 m,冷冻水传热损失较大,造成冷量在输送过程中散失。见表3,七采区管网冷量损失1.08×106 J/s,五采区管网冷量损失1.4×105 J/s,总冷量损失1.22×106 J/s。表 3 系统冷冻水管道冷量损失Table 3. Cooling loss of system chilled water pipeline七采区管道 冷量损失
Qr/(103 J·s−1)五采区管道 冷量损失
Qr/(103 J·s−1)E1-2 122.47 M1-2 24.23 E2-3 133.61 M2-3 9.22 E3-4 147.93 M3-4 7.32 E4-5 182.24 M3-5 3.65 E5-6 34.26 M5-6 7.59 E5-7 35.07 M2-7 7.83 E7-8 100.80 M7-8 6.88 E2-9 113.66 M7-9 3.74 E8-10 43.90 M9-10 3.77 E10-11 25.52 M9-11 2.13 E8-12 21.31 M11-12 2.25 E12-13 19.02 M11-13 2.14 E12-14 14.31 M1-14 21.62 E9-15 85.47 M1-15 37.58 总计 1079.59 总计 139.96 七采区管网冷量损失占总冷量损失的88.51%,这是由于管网输水范围大距离远,造成管路温升大。五采区管网各工作面比较集中,冷冻水能有效供给到需冷地点,降温效果良好。摩擦冷量损失和传热冷量损失分别为5.39×105、6.805×105 J/s。
2.5 实测数据对比
结合赵楼煤矿降温管网2022年7月末端空冷器入口冷冻水流量、温度实测数据,见表4。管网冷冻水实际流量低于计算值,主要原因是冷冻水在管道中流动时会有一定的损耗和流失,同时管道使用时间较长,管内存在沉淀物等导致水流量下降,实际测量数据偏低。管网冷冻水入口温度较计算值高2~3 ℃,这是由于冷冻水管网在管路连接处采用保温海绵,其冷量损失较大;同时管网冷量损失受多种因素影响,如流量、管长、管径等,实际冷冻水流量偏低也是造成管网冷冻水温度升高的原因。
表 4 空冷器入口水流量、温度实测数据Table 4. Measured data of inlet water flow and temperature of air cooler测点位置 流量Q/(m3·s−1) 温度t/℃ 7302运输巷 0.005 14.2 7302轨道巷 0.0024 13.4 7302工作面 0.0043 16.7 5304轨道巷 0.002 7.7 5304运输巷 0.002 8.6 5305运输巷 0.0008 9.5 5309运输巷 0.0035 7.6 5309轨道巷 0.0032 7.9 南部2号辅助运输大巷 0.001 9.2 3. 管网优化方法
3.1 阀门调节
在制冷降温管网中,由于支路较多、各个末端空冷器串联布置,当某些支路或空冷器流量调节时会引起系统冷冻水流量改变,其他设备的流量也随之发生改变,偏离设计要求,导致管网系统处于动态水力失衡状态。造成冷冻水流量不受控的现象,对制冷效果造成了巨大影响。
在管网中,各末端空冷器的流量取决于管网阻力特性系数。系统任一管段阻力特性发生变化,位于该管段之后的区段流量成一致性变化[21]。采用安装水力平衡阀来解决管网水力失调问题,将动态平衡阀设在供水总管路上、静态平衡阀设在供水支路上(图8)。阀门的直接调节对象为冷冻水流量,如果流量比目标流量小,则加大阀门开度使出流量增大,反之则减小阀门开度。
在降温管路中使用静态平衡阀对支路阻力进行调节,测量阀门前后压降变化,根据阀门开度计算调节后支路流量,使空冷器的流量达到设计值,调节后不随意改变阀门开度,实现系统的静态水力平衡。当管路总流量发生变化,由于动态平衡阀的恒定流量作用,末端设备的流量保持不变[22],即空冷器的实际流量不受其干扰,实现动态水力平衡。
阀门调节优化管网运行后,七采区降温系统调整了管网流量(表5),增加了由东向西的输水量,流入7302工作面冷冻水流量为0.009 m3/s。五采区降温系统经过水力调节阀重新分配输水量后,工作面最小流量0.005 m3/s,解决了降温系统水力失调问题。
表 5 系统管网管径改造后流量变化Table 5. Fow change after the pipe diameter transformation of the system pipe network管段 流量Q/(m3·s−1) 管径D/m 管段 流量Q/(m3·s−1) 管径D/m E1-2 0.040 0.377 M1-2 0.028 0.325 E2-3 0.032 0.377 M2-3 0.009 0.325 E3-4 0.032 0.325 M3-4 0.005 0.159 E4-5 0.032 0.325 M3-5 0.005 0.325 E5-6 0.010 0.219 M5-6 0.005 0.159 E5-7 0.022 0.325 M2-7 0.018 0.325 E7-8 0.022 0.325 M7-8 0.005 0.159 E2-9 0.011 0.159 M7-9 0.014 0.325 E8-10 0.009 0.159 M9-10 0.005 0.159 E10-11 0.009 0.159 M9-11 0.009 0.325 E8-12 0.012 0.273 M11-12 0.005 0.159 E12-13 0.009 0.159 M11-13 0.005 0.159 E12-14 0.003 0.159 M1-14 0.006 0.159 E9-15 0.008 0.159 M1-15 0.006 0.159 3.2 管径优化
结合水力损失情况,对管网中不利管段改造优化。南部1号辅助运输大巷、二集辅助巷和南部2号辅助运输大巷的管径增加至0.325 m,减小了各管段比摩阻[23],管网总水力损失从30.93 m减少到20.44 m,减少了35%。管内流体摩擦生热量与管段水力损失成正比,冷冻水温升减小,工作面降温明显。具体管径改造方案见表5。
3.3 离心泵功率
在流量和管径调节后,管网阻力发生变化,管路压降变化。最不利管路选择离制冷机组最远、埋深最深、节点水压最低和需水量最大的管路。七采区系统最不利管路:南部1号辅助运输大巷—二集辅助巷—七采区辅运巷下段,东部泵站与最不利末端相对高差90.02 m;五采区系统最不利管路:南部2号辅助运输大巷—2号辅运6号联络巷—5309运输巷,中部泵站与最不利末端相对高差81.84 m;产生局部损失管件有三通、变径管等。
七采区管路特性方程:
$$ H = 90.02 + 8.9 \times {10^3}{Q_{\mathrm{v}}^2} $$ (8) 五采区管路特性方程:
$$ H = 81.84 + 1.26 \times {10^4}{Q_{\mathrm{v}}^2} $$ (9) 式中:$ H $为管路的扬程,m;$ {Q_{\mathrm{v}}} $为管路系统输送流体的流量,m3/s。
离心泵工况点由离心泵特性曲线和管路特性曲线共同决定。将MD280-43*5型离心泵的特性曲线与管路的特性曲线绘于同一坐标图上,两线的交点成为泵在该管路上的最佳工作点。如图9所示,当流量在0.102 m3/s和0.094 m3/s为最佳工作点,这一点泵在管路中效率最高。离心泵流量调整为0.085~0.112 m3/s,扬程为183~195 m,在离心泵的效率大于90%高效区间。
3.4 调节效果对比分析
井下制冷硐室优化后冷冻水泵最佳流量是0.085~0.112 m3/s,在此范围内,降温管网最不利管路冷冻水供应满足设计流量,同时进入工作面冷冻水流量大于0.005 m3/s。当扬程在183~195 m范围内,能够克服管网因阻力损失造成的压降,使降温系统能够平稳运行。
4. 结 论
1) 赵楼煤矿井下制冷降温管网中冷冻水流量分配不均,末端最大流量为0.017 m3/s,最小流量为0.001 m3/s;管网中水头损失大的位置为一集轨道下山巷道、7302运输巷和中部辅运大巷。
2) 矿井降温系统分别给五采区和七采区共4个工作面供冷。七采区降温系统冷量损失为1.08×106 J/s,温升0.75 ℃/100 m;五采区降温系统冷量损失为1.4×105 J/s,温升0.12 ℃/100 m;摩擦和传热冷量损失分别是5.39×105 J/s和6.805×105 J/s;冷冻水在输送中末端温度最高为13.9 ℃。
3) 提出管道−泵阀联调的管网优化方法,采用动态平衡阀对南部1号辅助运输大巷和二集辅助巷冷冻水流量恒定在0.022~0.04 m3/s范围,增加了七采区自东向西输水量,实现了末端空冷器流量稳定;采用静态平衡阀调节管网支路阻力,将工作面冷冻水流量增大至0.005 m3/s;结合水力损失情况,将南部1号辅助运输大巷、南部2号辅助运输大巷和二集辅助巷的管径增加至0.325 m,管网总水力损失从30.94 m减少到20.44 m,减少了35%;离心泵流量调整为0.085~0.112 m3/s,扬程调整为183~195 m,在离心泵高效运行区间。
-
表 1 七采区降温管网信息
Table 1 Information of cooling pipe network in Seventh mining area
管段 管段位置 管长
L/m管径
D/m标高/m 流量
Q/(m3·s−1)Si Si+1 E1-2 南部1号辅助
运输巷178.8 0.377 −858.3 −857.1 0.04 E2-3 南部1号辅助
运输巷810.7 0.377 −857.1 −855.6 0.023 E3-4 南部1号辅助
运输巷956.22 0.325 −855.6 −853.5 0.023 E4-5 南部1号辅助
运输巷1134.5 0.219 −853.5 −923.9 0.023 E5-6 南部1号辅助
运输巷341.02 0.219 −923.9 −853.4 0.012 E5-7 二集辅助巷 168.52 0.219 −923.9 −923.3 0.012 E7-8 二集辅助巷 1481.3 0.273 −923.3 −927.1 0.012 E2-9 一集轨道
下山898.21 0.273 −857.1 −910.0 0.017 E8-10 7302运输巷 1596.5 0.159 −927.1 −792.0 0.006 E10-11 7302工作面 300 0.159 −792 −861.9 0.006 E8-12 二集辅助巷 325 0.159 −927.1 −924.9 0.006 E12-13 7302轨道
巷道1250 0.159 −924.9 −829.7 0.003 E12-14 二集辅助巷 370 0.159 −924.9 − 1015.0 0.003 E9-15 一集轨道
下山560.67 0.159 −910.0 −945.0 0.017 表 2 五采区降温管网信息
Table 2 Information of cooling pipe network in Fifth mining area
管段 管段位置 管长
L/m管径
D/m标高/m 流量
Q/(m3·s−1)Si Si+1 M1-2 南部2号辅助
运输大巷613 0.325 −848.9 −849.2 0.015 M2-3 南部2号辅助
运输大巷419 0.159 −849.2 −822.5 0.007 M3-4 5309运输巷 675 0.159 −822.5 −892.2 0.004 M3-5 南部2号辅助
运输大巷274 0.159 −822.5 −782.7 0.004 M5-6 5309轨道巷 693 0.159 −782.7 −865.1 0.004 M2-7 2号辅助6号
联络巷250 0.159 −849.2 −840 0.007 M7-8 5304轨道巷 993 0.159 −840 −808.9 0.003 M7-9 2号辅助6号
联络巷218 0.273 −840 −803.3 0.004 M9-10 5304运输巷 980 0.159 −803.3 −772.2 0.002 M9-11 2号辅助6号
联络巷207 0.273 −803.3 −778.4 0.002 M11-12 5305运输巷 1110 0.159 −778.4 −732.8 0.001 M11-13 2号辅助6号
联络巷700 0.159 −778.4 −761.2 0.001 M1-14 中部辅运
大巷855 0.159 −848.9 −900 0.010 M1-15 南部2号辅助
运输大巷1236 0.159 −848.9 −853.3 0.015 表 3 系统冷冻水管道冷量损失
Table 3 Cooling loss of system chilled water pipeline
七采区管道 冷量损失
Qr/(103 J·s−1)五采区管道 冷量损失
Qr/(103 J·s−1)E1-2 122.47 M1-2 24.23 E2-3 133.61 M2-3 9.22 E3-4 147.93 M3-4 7.32 E4-5 182.24 M3-5 3.65 E5-6 34.26 M5-6 7.59 E5-7 35.07 M2-7 7.83 E7-8 100.80 M7-8 6.88 E2-9 113.66 M7-9 3.74 E8-10 43.90 M9-10 3.77 E10-11 25.52 M9-11 2.13 E8-12 21.31 M11-12 2.25 E12-13 19.02 M11-13 2.14 E12-14 14.31 M1-14 21.62 E9-15 85.47 M1-15 37.58 总计 1079.59 总计 139.96 表 4 空冷器入口水流量、温度实测数据
Table 4 Measured data of inlet water flow and temperature of air cooler
测点位置 流量Q/(m3·s−1) 温度t/℃ 7302运输巷 0.005 14.2 7302轨道巷 0.0024 13.4 7302工作面 0.0043 16.7 5304轨道巷 0.002 7.7 5304运输巷 0.002 8.6 5305运输巷 0.0008 9.5 5309运输巷 0.0035 7.6 5309轨道巷 0.0032 7.9 南部2号辅助运输大巷 0.001 9.2 表 5 系统管网管径改造后流量变化
Table 5 Fow change after the pipe diameter transformation of the system pipe network
管段 流量Q/(m3·s−1) 管径D/m 管段 流量Q/(m3·s−1) 管径D/m E1-2 0.040 0.377 M1-2 0.028 0.325 E2-3 0.032 0.377 M2-3 0.009 0.325 E3-4 0.032 0.325 M3-4 0.005 0.159 E4-5 0.032 0.325 M3-5 0.005 0.325 E5-6 0.010 0.219 M5-6 0.005 0.159 E5-7 0.022 0.325 M2-7 0.018 0.325 E7-8 0.022 0.325 M7-8 0.005 0.159 E2-9 0.011 0.159 M7-9 0.014 0.325 E8-10 0.009 0.159 M9-10 0.005 0.159 E10-11 0.009 0.159 M9-11 0.009 0.325 E8-12 0.012 0.273 M11-12 0.005 0.159 E12-13 0.009 0.159 M11-13 0.005 0.159 E12-14 0.003 0.159 M1-14 0.006 0.159 E9-15 0.008 0.159 M1-15 0.006 0.159 -
[1] 何满潮,徐 敏. HEMS深井降温系统研发及热害控制对策[J]. 岩石力学与工程学报,2008,27(7):1353−1361. doi: 10.3321/j.issn:1000-6915.2008.07.007 HE Manchao,XU Min. Research and development of HEMS cooling system and heat-harm control in deep mine[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(7):1353−1361. doi: 10.3321/j.issn:1000-6915.2008.07.007
[2] 蓝 航,陈东科,毛德兵. 我国煤矿深部开采现状及灾害防治分析[J]. 煤炭科学技术,2016,44(1):39−46. LAN Hang,CHEN Dongke,MAO Debing. Current status of deep mining and disaster prevention in China[J]. Coal Science and Technology,2016,44(1):39−46.
[3] 李 伟. 深部煤炭资源智能化开采技术现状与发展方向[J]. 煤炭科学技术,2021,49(1):137−145. LI Wei. Current status and development direction of intelligent mining technology for deep coal resources[J]. Coal Science and Technology,2021,49(1):137−145.
[4] 隋秀华,曾现伟,黄云前,等. 高温采煤工作面集中供冷和分段供冷方式对比分析[J]. 煤炭科学技术,2016,44 (12):91−95. SUI Xiuhua,ZENG Xianwei,HUANG Yunqian,et al. Comparison analysis on centralized cooling and sectional cooling methods of high temperature coal mining face[J]. Coal Science and Technology,2016,44 (12):91−95.
[5] 岳丰田,刘存玉,魏京胜,等. 建井期间矿井降温系统能耗分析与优化[J]. 煤炭科学技术,2014,42(4):57−60. YUE Fengtian,LIU Cunyu,WEI Jingsheng, et al. Energy consumption analysis and optimization of mine cooling system during mine construction[J]. Coal Science and Technology. 2014,42(4):57−60.
[6] 魏杰立,张净玉,叶大海,等. 基于图论的燃油管路稳态流动计算方法[J]. 航空动力学报,2017,32(10):2344−2354. WEI Jieli,ZHANG Jingyu,YE Dahai,et al. Stable Hydraulic computation method for fuel pipe system based on the graph theory[J]. Journal of Aerospace Power,2017,32(10):2344−2354.
[7] 郑进福,周志刚,赵加宁,等. 供热系统动态热特性研究[J]. 建筑科学,2019,35(6):122−127. doi: 10.13614/j.cnki.11-1962/tu.2019.06.18 ZHENG Jinfu,ZHOU Zhigang,ZHAO Jianing,et al. Research on dynamic thermal characteristics of district heating system[J]. Building Science,2019,35(6):122−127. doi: 10.13614/j.cnki.11-1962/tu.2019.06.18
[8] RAHIMI M,HAGHIGHI A. A graph portioning approach for hydraulic analysis-design of looped pipe networks[J]. Water Resources Management,2015,29:5339−5352. doi: 10.1007/s11269-015-1121-9
[9] 窦 明,石亚欣,于 璐,等. 基于图论的城市河网水系连通方案优选:以清潩河许昌段为例[J]. 水利学报,2020,51(6):664−674. DOU Ming,SHI Yaxin,YU Lu,et al. Optimization of connecting schemes for urban river networks based on graph theory:a case study of Xuchang section of Qingying River[J]. Journal of Hydraulic Engineering,2020,51(6):664−674.
[10] SABZKOUHI A M,HAGHIGHI A. Uncertainty analysis of pipe-network hydraulics using a many-objective particle swarm optimization[J]. Journal of Hydraulic Engineering,2016,142(9):04016030. doi: 10.1061/(ASCE)HY.1943-7900.0001148
[11] CHICHERIN S,ZHUIKOV A,JUNUSSOVA L. The new method for hydraulic calculations of a district heating (DH) network[J]. Energy,2022,260:125071. doi: 10.1016/j.energy.2022.125071
[12] CHANG C,ZHANG X. Optimal local pressure control in water distribution network based on valve location[C]//2018 7th International Conference on Energy,Environment and Sustainable Development (CEESD2018). Atlantis Press,2018:8−13.
[13] 由世俊,米雷洋,王雅然,等. 集中供热管网的非稳态水力建模与动态响应分析[J]. 天津大学学报(自然科学与工程技术版),2019,52(8):849−856. YOU Shijun,MI Leiyang,WANG Yaran,et al. Unsteady hydraulic modeling and dynamic response analysis of a district heating network[J]. Journal of Tianjin University (Science and Technology),2019,52(8):849−856.
[14] HONG S P,KIM T,LEE S. A precision pump schedule optimization for the water supply networks with small buffers[J]. Omega,2019,82:24−37. doi: 10.1016/j.omega.2017.12.001
[15] 高铁军,赵 明,毛亚纯. 供水管网全网关阀预案与关键阀门分析[J]. 哈尔滨工业大学学报,2018,50(2):94−99. doi: 10.11918/j.issn.0367-6234.201608028 GAO Tiejun,ZHAO Ming,MAO Yachun. Analysis of valve closing schemes and key valves in water distribution networks[J]. Journal of Harbin Institute of Technology,2018,50(2):94−99. doi: 10.11918/j.issn.0367-6234.201608028
[16] 周 璇,刘国强,王晓佩,等. 基于动态水力平衡的集中供热系统二次管网节能优化调节方法研究[J]. 中南大学学报(自然科学版),2019,50(1):220−233. ZHOU Xuan,LIU Guoqiang,WANG Xiaopei,et al. Study on optimal adjustment method of second pipe network in central heating system based on dynamic hydraulic balance[J]. Journal of Central South University (Science and Technology),2019,50(1):220−233.
[17] 刘 洪,赵晨晓,葛少云,等. 基于精细化热网模型的电热综合能源系统时序潮流计算[J]. 电力系统自动化,2021,45(4):63−72. LIU Hong,ZHAO Chenxiao,GE Shaoyun,et al. Sequential power flow calculation of power-heat integrated energy system based on refined heat network model[J]. Automation of Electric Power Systems,2021,45(4):63−72.
[18] 张 鑫,刘绍忠. 基于阻抗理论的管网水力平衡调节新思路[C]//中国建筑学会建筑热能动力分会第二十一届学术交流大会论文集. 安徽,2019:19−22. ZHANG Xin,LIU Shaozhong. A new idea of hydraulic balance regulation of pipe network based on impedance theory[C]//Proceedings of the 21st Academic Exchange Conference of Building Thermal Power Branch of Architectural Society of China. Anhui,2019:19−22.
[19] 旷金国,王朝晖,罗曙光. 区域供冷系统外管网冷量损失分析[J]. 暖通空调,2021,51(12):16−21. KUANG Jinguo,WANG Chaohui,LUO Shuguang. Analysis on energy loss of external pipe network in district cooling systems[J]. Heating Ventilating & Air Conditioning,2021,51(12):16−21.
[20] 胡春胜. 矿井空调系统中输冷管道的热力计算方法[J]. 煤矿设计,1990(8):41−44. HU Chunsheng. Thermal calculation method of cooling pipeline in mine air conditioning system[J]. Coal Engineering,1990(8):41−44.
[21] 王 鹏,张 丽. 浅谈阀门调节对热网系统水力工况的影响[J]. 暖通空调,2022,52(S2):179−183. WANG Peng,ZHANG Li. Discuss the influence of valve adjustment on the hydraulic working conditions of heat network system[J]. Heat Ventilating & Air Conditioning,2022,52(S2):179−183.
[22] YU J,LIU Q,ZHAO A,et al. A distributed optimization algorithm for the dynamic hydraulic balance of chilled water pipe network in air-conditioning system[J]. Energy,2021,223:120059. doi: 10.1016/j.energy.2021.120059
[23] 范钟引,熊帝战,杨 玲,等. 区域供冷供热系统输配管网管径优化设计研究[J]. 暖通空调,2021,51(8):50−55. FAN Zhongyin,XIONG Dizhan,Yangling,et al. Optimal design of pipe diameters for distribution pipe networks of district cooling and heating systems[J]. Heating Ventilating & Air Conditioning,2021,51(8):50−55.