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]提出了巷道围岩协同锚固的研究思路,靖洪文等[8]以物理模拟试验探索深部巷道围岩锚固结构从开始承载至整体失稳全过程,王襄禹等[9]通过引入弹黏塑性本构模型,理论分析了锚注结构的承载特性以及巷道流变变形特性;在围岩改性方面,康红普等[10]介绍了水泥、化学注浆加固的应用条件、注浆加固设计及注浆效果评价,张庆松等[11]研究浆液扩散路径的多孔介质渗透注浆机理,设计一套多孔介质渗透注浆扩散模拟实验装置,张振峰等[12]开展了千米深井巷道高压劈裂注浆改性技术研究与实践,测试评价了注浆改性后的效果;在围岩卸压方面,康红普等[13]介绍了巷道卸压法分类及卸压机理、影响因素和适用条件,吴拥政等[14]提出了深部冲击地压巷道“卸压−支护−防护”协同防控原理及技术,徐佑林等[15]提出受强动压影响的“三软”煤层巷道“卸−转−固”围岩综合控制理论及技术。
上述研究基本都是围绕巷道围岩与支护结构的作用机理开展的,形成了围岩松动载荷[16]、控制围岩变形[17]、承载结构[18]、围岩改性[19]及应力控制理论[20]等众多围岩控制理论,并根据巷道围岩条件和控制理论提出了一系列的支护技术,并且在现场工程实践取得了较好的验证。然而,近年来,随着材料、施工机械设备的发展,深部复杂地质条件巷道围岩的支护方法却没有较大的变革与发展,现有的巷道支护理论与技术均有各自限制与适用条件,如强构造影响、“三软”煤层、破碎围岩及强动压影响的大变形巷道,围岩可锚性较差,无法形成有效的承载结构并发挥围岩自身承载能力,难以保证支护效果。
综上所述,现有的支护手段不能完全满足深部巷道围岩的支护要求,亟需通过新材料、新工艺优化现有支护技术与工艺。笔者提出高强度再造承载结构巷道支护技术,并以具体煤矿为例,通过对两帮松软煤体进行置换,让顶板、两帮充填墙和底板构成高强度承载结构,采用智能算法优化了承载结构支护参数与工艺并进行井下工业性试验,取得了良好的工业试验效果,同时为相似条件巷道支护提供了参考。
1. 典型难支护巷道类型
在煤矿巷道支护领域,众多学者和工程技术人员开展了大量的研究,采用支护法、加固法、应力控制法和联合支护法等方法已经解决了大量的煤矿巷道支护难题。团队经过长时间的现场调研及分析,发现由于地质构造的复杂程度不同、开采深度和强度不同、开采布局不合理等多种原因,当前存在以下一些典型难支护巷道,亟需进行支护技术攻关。现将调研的部分煤矿支护难题列举如下。
1.1 受强构造影响的难支护巷道
贵州水城县格木底向斜边缘分布了大树脚、吉源、阿嘎、小牛等20余个煤矿,受该向斜的影响,这些煤矿基本上都是急倾斜煤层,煤层最大倾角达到86°,煤层埋深大,基本在600~800 m,开采难度极大,巷道支护极为困难。主要原因是受构造影响的地区,特别是靠近向斜轴部等构造线的部位,构造运动形成的应力场在这些地区特别是深部岩层中个点的应力(包括大小和方向)将是自重应力场和构造应力场在该点的综合叠加。国外开采深度较大的矿井进行了井下应力实测,测量结果表明水平方向最大主应力可比垂直应力高出3~5倍,最高可达8~10倍。同时,受构造运动影响的原始应力场是一个应力分布很不均匀且在一定程度上受时间影响的应力场,在这样的原始应力场中进行采掘工作,必须考虑时间、地点和构造条件[20-24]。根据传统矿压理论推导,垂直应力为15~20 MPa的构造水平应力最高可达200 MPa,巷道围岩难以抵抗。
对该地区煤矿(吉源煤矿)进行调研发现基本都存在巷道修复困难的问题。以水城县吉源煤矿为例,吉源煤矿位于扬子准地台黔北台隆六盘水断陷或威宁北西向构造变形区的格目底向斜南翼东端仰起处,二采区轨道下山在1400轨道石门开口往里702 m处以318°方位,−26°坡度至标高+1 333 m落平。掘处以318°方位,−26°坡度至标高+1 333 m落平,掘进全断面为直墙圆弧拱形,采用U型棚+锚网索支护。目前巷道所采用的复合支护方式,效果不理想,矿压显现剧烈,局部严重地段维修周期在1~2个月,浪费大量的人力和财力,严重影响矿井的正常生产。经过调研分析可知,主要有以下方面:矿井构造应力大,加上二采区埋深较深(660 m),巷道处于高应力状态。同时,由于岩层破碎,巷道压力大,导致反复维修,扩帮、拉底维修量大,巷道两帮最大变形量达1.5 m,底鼓达1.0 m,局部顶板下沉达1.0 m,如图1所示。围岩主要为泥岩,强度低、遇风、水容易风化破碎,同时局部受到采区的采动影响,巷道维护难度大。
1.2 近距离煤层群开采强动压影响巷道
课题组对汪家寨煤矿、土城矿、山角树煤矿、湾田煤矿、金河煤矿、打牛场煤矿、仲恒煤矿等30余个煤矿进行调研发现,受井田范围限制,多数巷道布置在中部,采用两翼开采,加之近距离煤层群赋存条件,不断跳采,同时部分老矿井由于开采年限长,系统极其复杂,形成复合应力状态,在采动应力作用下,巷道围岩受多次动压影响,已经完全破碎,围岩对应力变得极为敏感,变形极为严重。出现支护构件严重破坏,顶板、两帮变形严重(图2)。
汪家寨煤矿位于六盘水煤田大河边向斜中段,开采上二叠系宣威组煤层,可采及局部可采煤层11层,属缓倾斜近距离煤层群,煤系地层主要由细砂岩、粉砂岩、黏土岩、泥岩等组成,软岩约占65%。井田地质条件复杂,属煤与瓦斯突出矿井。开拓、准备巷道一般采用锚网喷或U型棚支护,采区巷道采用联合布置方式,各煤层采用区段石门联络,井筒、大巷、采区上、下山均布置于煤组底板,绝大部分处于软岩中。
1)由于巷道围岩受采动的原因,巷道矿压显现极其明显,造成巷道断面严重变形,出现顶板破碎、垮落、两帮及顶底板位移量大的现象,部分巷道断面收缩率达到60%以上甚至闭口,给通风、行人及安全造成重大影响。
2)平硐井区巷道失修严重,主干巷道尤为明显,巷道破坏表现具有显著的高应力软岩特性,围岩松散破碎,而且受历史条件及支护手段限制,造成支护方式与围岩条件不相适应。
3)采用锚网喷、U型钢可缩性支架等支护形式对巷道进行多次翻修,均未达到理想效果。
现以2105轨道石门为例进行简要介绍,该巷位于+1 547 m水平,埋深153 m,长度约为520 m,岩性为泥质细粉砂岩,倾斜上方分别有210706采空区、211106回采工作面、2105回风石门、2105运输石门,受采动影响,矿压显现明显,巷道底鼓严重,由于水沟不畅局部地点巷道积水,受裂隙影响,2105回风石门积水沿着围岩裂隙进入该巷,造成局部顶板淋水,如图3所示。该巷目前采用断面积13.9 m2(29 kg/m)U型棚进行支护,顶帮用钢筋网、铁丝网及铁背板穿严接实,棚距700 mm,施工方法采用风镐配合手镐进行施工,排矸方式为人工装矸上矿车进行出矸,每月进度20 m,施工缓慢。
1.3 “三软”煤层等软岩巷道
软岩支护是当前巷道支护的普遍性难题,“三软”煤层巷道支护尤其突出。据统计,在我国“三软”煤层开采中,每年因顶板事故死亡的人数占煤矿事故的45%左右,其中工作面顶板事故占70%~80%。课题组对火铺、金佳等该类条件的32个煤矿相关巷道进行调研发现,各个矿区均对“三软”煤层综合机械化开采进行积极的尝试与探索,如盘江精煤股份有限公司火铺矿对其“三软”煤层-17号煤层工作面的ZY3800/15/33支架进行改进,力图解决“三软”煤层液压支架埋头的问题;六盘水恒鼎实业有限公司积极总结洪兴煤矿111703综采工作面的开采经验,改进“三软”煤层生产工艺,使得111703工作面的月推进度显著提高。但是由于生产矿井长期以来缺乏对“三软”煤层围岩力学性状、矿压显现规律以及覆岩移动特征的探索,面临的主要技术难点在于:
1)该类巷道底板黏土矿物质量分数59%~85.69%,其中伊/蒙间层矿物(I/S)质量分数65%~92%,顶底板岩性为泥质粉砂岩和粉砂岩,部分为煤,遇水易软化、膨胀和泥化,围岩可锚性、稳定性差,常规锚杆(索)等主动支护方式在极软弱和破碎围岩条件下其可锚性较差,无法提供有效的约束力,如图4所示。
2)采用架棚、液压支柱等被动支护措施,而由于“三软”煤层底板软弱,U型棚以及液压支柱在支承压力作用下被压入软弱底板,不仅无法提供有效的支撑力,而且严重影响工作面的推进速度。因此,以目前的支护手段仍旧无法对“三软”煤层采动巷道围岩进行有效控制,亟需综合采用新的支护手段、工艺与材料方能以解决“三软”煤层采动巷道围岩控制困难问题。
3)“三软”煤层覆岩移动规律尚未完全掌握,其覆岩移动规律与非“三软”煤层覆岩移动规律存在一定差异;破碎围岩条件下支架−围岩关系不明确,“三软”煤层采场围岩力学性状较为特殊。
综合分析可知,采用架棚、单体液压支柱等被动支护措施,而由于“三软”煤层底板软弱,U型棚以及液压支柱在支承压力作用下被压入软弱底板,不仅无法提供有效的支撑力,而且导致巷道两帮和顶板进一步变形破坏,易出现底鼓、顶底和两帮大变形,支架扭曲变形、“U”卡崩坏和出现尖顶等破坏变形情况(图5)。
2. 现有巷道围岩控制技术及原理
巷道在开挖后,原有的应力平衡状态被打破,围岩应力重新分布,受采动应力的影响,巷道围岩出现应力集中现象,超过围岩的承载能力,引起巷道围岩的变形破坏,形成了破碎区、塑性区及弹性区3个较为明显的分区破坏[24-25],同时巷道围岩变形具有流变与扩容特征,持续弱化围岩承载能力降低,围岩体变形加剧,严重影响深部巷道围岩的长期稳定。康红普院士按照围岩控制部位和原理将支护类型分为5种,分别是巷道围岩表面支护型、巷道围岩锚固型、巷道围岩改性型、巷道围岩卸压型及联合控制型[26]。
从形式上,我国现有的支护主要划分为被动支护型和主动支护型2种,被动支护型主要采用支架、喷射混凝土、砌碹等巷道表面支护方式,都是靠支护本身的强度、刚度来支撑围岩体;主动支护型主要采用锚固技术、围岩注浆改性及卸压技术等当中的一种或多种相互组合等巷道内部支护形式,通过对围岩体强度的改变,使围岩体自身对围岩进行支撑,强化围岩应力状态,提高围岩承载能力。从围岩控制原理划分为加固和应力控制2大类,加固主要是给巷道表面提供约束,以及锚网、注浆加固等提高围岩强度、刚度、完整性和稳定性的支护加固法。应力控制主要是通过卸压等手段降低围岩应力、转移围岩高应力、均匀化围岩应力避免出现过高集中应力、降低围岩偏应力和应力梯度等。现分别简单总结如下:
2.1 巷道表面支护(被动支护)
对巷道表面施加约束力来控制围岩变形,主要有支架、砌碹及喷射混凝土等支护方式,如图6所示。主要是依靠支护本身的强度、刚度来支撑围岩体,在开采阶段初期可以很好地解决支护难题,但没有最大程度地调动巷道围岩的自承载能力,随着时间的推移,破坏后难以修复进行二次利用,且支护成本较高、工人劳动强度较大,且架棚支护强度有限,无法完全控制巷道围岩变形。
2.2 巷道围岩内部支护(主动支护)
注浆加固主要是通过注浆管将浆液扩散到加固区域,经过渗透→充填胶结→压密充实→凝结固化的过程,从而恢复围岩完整性,使锚杆锚索施加的预应力能有效传导和扩散,形成有效的支护应力场,起到支护作用。锚杆锚索支护主要是通过打设锚杆(索)压紧锚固区围岩,改善巷道围岩的黏聚力、内摩擦角等物理力学参数,在浅部围岩形成一个“加固拱”,发挥围岩自身承载能力,保证巷道围岩稳定,如图7所示。但若巷道覆岩中含有大量易膨胀的黏土矿物的松、散、软、弱岩层,锚杆、锚索的可锚性大大降低,无法有效调动巷道围岩深部自身强度,从而限制了锚杆、锚索等主动支护手段在松软破碎围岩巷道的应用。
2.3 卸压支护(应力控制)
主要采用卸压(爆破、切顶、钻孔、水力压裂等)来改善围岩应力状态,降低或转移应力集中区域,如图8a所示。通过卸压技术,巷道围岩应力重新分布,使应力集中的弹性区转移到围岩更深处,降低浅部围岩应力集中区域,起到了较好的控制作用,如图8b所示。
2.4 联合支护(主动+被动)
针对难支护巷道,目前大部分矿井巷道治理基本上采取联合控制方式,主要是主动支护方式结合被动支护方式,如U型钢棚+喷浆、锚网索等联合支护方式,锚杆、锚网索及注浆不适应现场地质条件的情况下,采取U型钢棚支护,弥补了锚杆、锚网索及注浆支护的缺点,一定程度上保证现场的安全生产需要,延长了巷道变形破坏周期,但在后期维修时影响较大。
3. 再造高强度承载结构支护技术思路
3.1 研究思路的提出
根据上述分析,部分难支护巷道采用现有的支护理论及技术不能彻底解决巷道围岩大变形的问题。现有支护技术原理不论是主动支护还是被动支护,其目的都是为了给开挖造成的应力环境变化提供一个抵抗变形的约束力,挤压锚固区内围岩,提高其物理力学性能,在围岩表面或者浅部形成一个具有自身承载能力的“加固拱”结构,让这个结构具备抵抗变形的能力,从而保证巷道围岩稳定。但U型钢等支护形式存在的缺点是单点受力,没有整体承载,虽然其材料强度足够高,但是受力条件不好,容易变形;而锚网索支护等经过多年的研究和工程实践证明具有较好的支护作用,但锚网索等主动支护施工质量是关键,由于锚网索支护是“隐蔽工程”,加上工人素质参差不齐,无法确保其施工质量,因而无法保证支护强度和支护效果。笔者正是基于这个观点和现实工程问题,结合主被动支护特征,提出系列再造承载结构快速支护成巷技术和工艺,重新构建的“加固拱”结构,依靠再造承载结构实现整体承载,整体抵抗变形,从根本上实现对围岩“彻底改性”,达到再造承载结构与围岩体强度、刚度的可适性,提高巷道整体稳定性。主要从以下3个方面进行研究。
1)基础理论研究:研究煤层开采后采动空间围岩应力分布、承载结构、岩体损伤的突变特征;建立巷道高强度再造承载结构拱力学模型,揭示巷道围岩−支护结构的耦合效应,形成巷道再造承载体围岩结构失效判别的新方法,确定受采动影响下高强度承载结构厚度、强度及结构。
2)工艺和技术:针对煤矿巷道掘进和返修过程效率低,安全性不高等问题,研究煤矿巷道快速成巷(掘进和支护)关键技术和工艺。
3)材料和装备:研发轻质高强新材料,高效快速的施工机具和装备。通过上述研究,形成集理论、技术、工艺、材料和自动化装备于一体的成套技术。
通过再造承载结构将整条巷道构建成一个具有高强度承载能力的结构体,承载结构整体受力,发挥再造岩体自身承载能力和围岩整体抗变形能力,保证巷道正常使用,为煤矿巷道支护提供参考。
3.2 巷道高强度承载结构再造方法
1)模板台车浇筑成型(图9)。借鉴隧道支护方法,针对煤矿巷道支护特点和空间特征,研究煤矿巷道模板台车,研发高韧性高强度的支护材料进行壁后充填,在围岩表面形成高强度可让压的再造承载结构,根据矿压大小考虑是否需要结合锚杆(索)支护,实现安全高效支护。
2)装配式预制件快速支护。在研发高强度高韧性支护材料的基础上,借鉴建筑上的装配式建筑思路,将支护材料做成预制结构(图10),在井下进行安装,利用传统的榫卯结构和自锁原理,形成具有高强承载能力的承载结构,结合锚杆(索)支护,实现快速支护。
3)堆喷快速支护。研发高强度、高韧性和高黏结能力的新型材料,采用喷射混凝土的方式进行堆喷,如图11所示,在巷道围岩表面形成一层高性能的承载结构,具备高承载能力,再根据矿压大小和巷道变形破坏因素,结合锚杆(索)支护和注浆加固等实现快速、安全和有效支护。
4. 井下工业性试验及效果分析
4.1 工程背景
龙宝煤矿11205运输下山布置在C12煤层中,平均厚度2.18 m,巷道顶板为粉砂岩,底板为粉砂质泥岩,由于巷道布置在煤层中,在水平应力作用下,两帮巷道变形大,每个月的巷道两帮移近量达到600 mm。
1)巷道原支护方案与状况。 11205运输下山巷原设计断面为矩形,巷道断面净宽4 200 mm,中高2 700 mm。巷道原支护采用锚杆(索)+铁丝网+钢带联合支护,支护设计如图12所示。锚杆间排距800 mm×800 mm,顶锚杆使用ø20 mm×2 500 mm的螺纹钢筋锚杆,两帮采用ø20 mm×2 000 mm的树脂锚杆,锚杆托盘规格为:150 mm×150 mm×10 mm(厚);锚索间排距1 600 mm×1 600 mm,每向前掘进1.6 m打2根锚索加强支护,锚索规格为ø17.8 mm×7 200 mm,锚索托盘采用300 mm×300 mm×15 mm(厚)的托盘。采用上述支护后,巷道仍出现剧烈变形,主要表现为巷帮大变形,两帮每月收缩400~600 mm。
2)巷道变形破坏原因分析。①巷道沿煤层布置,煤坚固性系数为0.6,稳定性差,巷道水平应力占主导,受到强烈水平应力挤压而产生变形破坏。②巷道顶板较好,但是巷道布置在煤层中,煤层软弱,没有较好的承载能力,在水平应力和垂直应力作用下,巷道两帮破碎区扩大锚杆锚索难以发挥有效支护作用。③巷道围岩锚固力低及注浆加固效果有限。由于围岩松软破碎,造成钻孔塌孔严重,导致锚杆预紧力不能得到有效扩散、锚杆锚固力较低,支护强度不够,严重影响巷道支护效果。
4.2 围岩控制方法确定
针对11205运输下山巷道变形严重问题,采用现有型钢和锚网索等支护手段难以取得较好的支护效果。采用井下实测、理论计算和数值计算等研究方法,结合加固拱、卸压支护等理论和技术,研究巷道整体结构受力特征及破坏规律,研发出充填新材料,提出墙体置换再造承载结构快速支护技术,即对两帮软弱煤体进行一定厚度的置换,如图13所示,在巷道两帮置换掉原来的软弱岩体,再造一定厚度的高强度充填墙,使顶板、两帮和底板形成一个闭合的完整结构,充分发挥承载结构的自身承载能力,抵抗围岩压力,同时,根据现场地质情况,结合理论分析在再造结构与围岩之间预留一定宽度作为卸压区,来达到减缓水平应力作用。
4.3 巷旁充填墙承载力学模型建立
为便于实际工程应用中计算巷旁充填墙承载能力,开展理论分析,如图14a所示,x为充填墙到围岩的距离,$ H $为矩形巷道埋深,$ \gamma $为上覆岩层容重,在巷道两侧建立巷旁充填墙,$ l_{1} $为墙体厚度,$ h $为高度,$ l_{2} $为巷旁充填墙与巷帮距离。
如图14a所示,该模型为对称模型,选取图14a中模型左边一半作为研究对象,进行受力分析,将该模型简化为梁模型,以巷道左上角为梁的支点,梁受垂直向下的力,为上覆岩层重力$ \gamma H $,而施工的巷旁充填墙提供向上的支撑力,如图14b所示。
根据力矩平衡可得,上覆岩层产生的力矩与巷旁充填墙产生的力矩平衡,从而使整体保证力学平衡。根据图14b,由上覆岩层产生向下的力矩MM可表示为
$$ M_{{\mathrm{M}}}=\int_0^{l_2+l_1+W / 2} x \gamma H {\mathrm{d}} x $$ (1) 根据图14c,墙体的支撑力q所形成向上力矩MZ可表示为
$$ M_{{\mathrm{Z}}}=\int_{l_2}^{l_2+l_1} x q {\mathrm{d}} x $$ (2) 根据力矩平衡,式(1)与式(2)相等,可得如下关系式:
$$ \int_0^{l_2+l_1+W / 2} x \gamma H {\mathrm{d}} x=\int_{l_2}^{l_2+l_1} x q {\mathrm{d}} x $$ (3) 式(3)进一步简化为:
$$ \gamma H\left(l_2+l_1+W / 2\right)^2=q\left(2 l_2 l_1+l_1^2\right) $$ (4) 可得巷旁充填墙所须提供的支撑力q为
$$ q=\frac{\gamma H\left(l_2+l_1+W / 2\right)^2}{\left(2 l_2 l_1+l_1^2\right)} $$ (5) 根据巷旁充填墙、巷道的几何参数及上覆岩层容重,结合式(5)可快速得到巷旁充填墙的最小强度q。根据煤和岩石物理力学性质测定方法(GB/T 23561.2—2009),在现场取芯、切割、打磨至光滑,然后在实验室进行单轴压缩试验得到相关力学参数,见表1,由上述计算公式,得到q值为39.78 MPa。因此,根据理论分析结果当巷旁充填墙体支撑能力达到39.78 MPa,能够控制巷道围岩变形,并保证巷旁充填墙体不被破坏。
表 1 岩土体和巷旁充填墙的力学参数Table 1. Mechanical parameters of rock and soil mass and roadside backfill wall岩性 密度/(kg·m−3) 体积模量/GPa 剪切模量/GPa 黏聚力/MPa 内摩擦角/(°) 抗拉强度/MPa 岩土体 2500 11.18 11.90 0.71 20 0.55 巷旁充填墙 2800 22.00 22.00 60.00 50 40.00 4.4 数值模拟分析
4.4.1 数值模拟计算模型建立
选取龙宝煤矿11205运输下山作为研究对象,11205运输下山上覆岩层厚度H为418 m,巷道高度h为2.8 m。边宽l2、巷旁充填墙厚度l1及巷道一半W/2的总和为3.7 m。工程应用中需要确定合理巷旁充填墙的厚度l1,墙体越厚,则其提供的支护强度越大,而经济上不合理;如果巷旁充填墙厚度l1偏小,则巷道有可能会坍塌(数值模拟表现为计算结果不收敛)。因此确定满足工程稳定的巷旁充填墙的最小厚度是本次数值模拟计算的目标。根据工程实际要求,该数值模拟计算优化目标函数及约束条件可表示为
$$ \begin{aligned} & \min f=l_1 \\ & \text { sub. }\left\{\begin{array}{l} 0 \leqslant l_2 \leqslant 0.5 \\ 0 < l_1 < 2 \\ l_1+l_2 \leqslant 2 \end{array}\right. \end{aligned} $$ (6) 为确定合理的巷旁充填墙、边宽几何参数,本文采用FLAC3D数值模拟软件开展数值模拟计算,根据工程实践,建立FLAC3D数值模拟计算模型,图15给出了FLAC3D数值模拟计算模型建立过程,选取模型一半(图15a)开展数值模拟计算,数值模拟计算边界条件如图15b所示。
数值模拟计算采用的模型为Mohr-Coloumb模型,岩土体和巷旁充填墙的力学参数见表1。
为实现巷道稳定前提条件下(FLAC3D中认为数值模拟计算模型的最大不平衡比率小于10−5时,则认为该模型稳定),满足约束条件(式(1)),求得巷旁充填墙最小厚度l1,采用Python编程语言结合优化算法与FLAC3D数组模拟计算软件,根据数值模拟计算模型收敛状态及目标函数,基于优化算法不断调整巷旁充填墙几何参数。数值模拟计算过程中计算模型的最大不平衡比率小于10−5时,则巷旁充填墙对巷道进行了有效支护,巷道将处于稳定状态,否则认为巷旁充填墙支护失败,巷道垮塌。数值模拟计算过程中运行步数为70 000步,确保数值模拟充分计算。采用飞蛾火焰优化算法确定合理巷旁充填墙厚度,该算法具有超参数少,计算收敛快的特点。
4.4.2 飞蛾火焰优化算法
飞蛾火焰优化算法(moth-flame optimization algorithm, MFO)[27-28]是由Seyedali Mirjalili在2015年提出,是一种新颖的群体智能算法,该算法由飞蛾和火焰2部分构成,通过定位导航机制来解决探索和寻优之间的平衡问题。飞蛾火焰优化算法具有参数简单、容易实现、鲁棒性好等优点,因此该算法自提出以来,便收到国内外学者的广泛关注,在各领域中得到广泛的应用。
在飞蛾火焰算法中,每一个飞蛾有一组最优的参数,而飞蛾的位置是由最优参数确定的,同时飞蛾的位置根据其相应的适应值进行排序,最优的位置则求解问题的最优解。
在飞蛾火焰优化算法中,飞蛾的位置是由优化参数确定,求解问题的维数就是变量的个数,其中飞蛾的位置矩阵可表示为
$$ L_M=\left[\begin{array}{cccc} m_{1,1} & m_{1,2} & \ldots & m_{1, d} \\ m_{2,1} & m_{2,2} & \ldots & m_{2, d} \\ \vdots & \vdots & & \vdots \\ m_{n, 1} & m_{n, 2} & \ldots & m_{n, d} \end{array}\right] $$ (7) 式中:n为飞蛾的数量(种群大小);d为变量的数量。飞蛾对应位置的适应度可根据目标函数计算得到,飞蛾位置适应度矩阵可表示为:
$$ F_M=\left[\begin{array}{c} F_{M, 1} \\ F_{M, 2} \\ \vdots \\ F_{M, n} \end{array}\right] $$ (8) 在飞蛾火焰算法中,每一个飞蛾对应一个火焰位置,飞蛾最新的位置会根据火焰的位置更新,一旦最新最合适的火焰位置发现时,此时最新的火焰位置可表示为
$$ F_f=\left[\begin{array}{cccc} f_{1,1} & f_{1,2} & \ldots & f_{1, d} \\ f_{2,1} & f_{2,2} & \ldots & f_{2, d} \\ \vdots & \vdots & & \vdots \\ f_{n, 1} & f_{n, 2} & \ldots & f_{n, d} \end{array}\right] $$ (9) 此时火焰位置的适应度函数可表示为
$$ F_f=\left[\begin{array}{c} F_{f, 1} \\ F_{f, 2} \\ \vdots \\ F_{f, n} \end{array}\right] $$ (10) 飞蛾的位置根据火焰的位置进行相应的更新,可表示为
$$ L_{{\mathrm{M}}, P}=S\left(L_{{\mathrm{M}}, P}, L_{{\mathrm{f}}, q}\right) $$ (11) 式中:$ L_{{\mathrm{M}}, P} $为第P个飞蛾位置;$ L_{{\mathrm{f}}, q} $为第q个火焰位置;$ S $为螺旋函数;该螺旋函数可表示为
$$ S\left(L_{{\mathrm{M}}, P}, L_{{\mathrm{f}}, q}\right)=D_p {\mathrm{e}}^{b t} \cos (2 \text{π} t)+L_{{\mathrm{f}}, q} $$ (12) 式中:$ D_{p} $为第p个飞蛾与第q个火焰的距离;b为螺旋系数;t为−1~1的随机数。
飞蛾火焰算法的具体流程可描述为
步骤1:初始化飞蛾种群个数及最大迭代次数,本文中飞蛾的数量为50,最大迭代次数为100。
步骤2:确定变量的上限及下限,基于变量上限及下限,产生初始解。
步骤3:计算飞蛾对应位置的适应度,作为初始火焰的适应值。
步骤4:根据式(22)更新飞蛾的位置,计算飞蛾位置的适应度,当找到更好的位置时,将该位置作为火焰的位置。
步骤5:根据下式减少火焰的数量,有助于提高飞蛾火焰算法局部搜索能力:
$$ f_n=\operatorname{round}\left(N-I \frac{N-1}{T}\right) $$ (13) 式中:N为火焰位置数量的最大值;I为当前的迭代次数;T为最大的迭代次数。
步骤6:当终止条件达到时,则停止计算,进入步骤7;否则进入步骤4。
步骤7:结束。
飞蛾优化算法具有超参数少、收敛速度较快的特点,因此本文中采用飞蛾火焰优化算法确定巷道巷旁充填墙的最优几何参数。
为便于实现飞蛾火焰优化算法实现巷旁充填墙几何参数优化,给出了具体的实现过程,如图16所示。
根据上述数值模拟计算流程,编写Python程序,实现FLAC3D数值模拟计算机飞蛾火焰优化算法结合。
4.4.3 数值模拟计算结果分析
通过飞蛾火焰优化算法100次迭代后,得到边宽l2为0.49 m,巷旁充填墙厚度l1为1.45 m。得到的最优巷旁充填墙几何参数满足约束条件,同时图17给出了巷旁充填墙垂直方向的应力分布。
分析图18可知,在巷旁充填墙右上端出现较为明显的应力集中现象,且巷旁充填墙的承载应力为30~36.098 MPa。为进一步确定巷旁充填墙最大压应力,编写FISH语言,获取巷旁充填墙最大压应力为36.098 MPa,数值模拟计算结果(36.098 MPa)与理论计算结果(39.78 MPa)吻合较好,且两者之间的误差为9.25%,该误差在可接受的范围之内。
为了进一步监测巷旁充填墙支护后,巷道顶底板变形量,在巷道布置了4个监测点,如图19a所示。其中p1、p2点与巷道边界距离为0.245 m(巷旁充填墙左边界与巷道边界中间),p3、p4则布置在巷道中间位置。
通过计算得到监测点p1、p2、p3及p4位移,如图19所示,其中p1和p3点变形较为接近,均为25 mm左右,其中p3点位移稍大。而底板监测点p2和p4变形量较小,其中又以p4点位移量较大,最大位移为10 mm左右,而p2点最大位移为5 mm左右。通过数值模拟试验可得,采用该巷旁充填墙几何参数,巷道顶底板变形较小,且变形范围在工程实践中是可以接受的。因此,说明当巷旁充填墙承载能力达到36.098 MPa时,通过采用巷旁充填墙体置换破碎围岩巷帮,能够有效控制破碎围岩巷道的变形。
4.5 充填材料研发
巷道处于高应力下且围岩较为松软,易造成巷道墙体大变形且围岩破碎严重,故团队研发了一种低成本、凝结时间短、塑性较好且承载强度较高的新材料,在破碎围岩巷道巷旁构建再造高强承载结构为充填墙。
材料主要以硫铝酸盐水泥、石膏及石灰为主,主要有A料、B料2部分,分别加入1.5∶1的水料比进行混合搅拌,如图20所示,分别对A料、B料进行单浆搅拌不凝固,两者混合后,A料和B料迅速发生反应并凝固,分别进行1 d、3 d、7 d及15 d的标准养护,进行单轴抗压强度测试,发现充填墙材料在标准养护7 d后,强度达到12 MPa左右,呈现出明显的塑性变形特征,可以维持较高的残余强度,可以从根本上满足巷道巷旁充填墙承载结构支护的要求。
4.6 巷旁充填墙支护工艺
根据数值模拟和理论计算分析,得到充填后方卸压区预留宽度490 mm,充填体宽度1 450 mm,再造承载结构支护如图21所示。
1)施工顶板锚杆,确保巷道扩刷安全。顶板采用锚杆支护,锚杆使用ø20 mm×2 500 mm螺纹钢锚杆,托盘规格为150 mm×150 mm×10 mm,间排距800 mm×800 mm。锚索规格为ø17.8 mm×7 200 mm,锚索托盘规格300 mm×300 mm×15 mm,间排距
1600 mm×1600 mm。2)将两帮松软围岩扩刷到设计范围,预留出充填和卸压位置。
3)立模进行充填。将充填袋吊挂在设计充填位置,以钢筋网连接并穿上对拉锚杆,打牢单体液压支柱进行约束。
4)充填材料按设计比例加水进行搅拌,利用管道输送到充填袋,浆液反应完后形成墙体,与巷道顶板和底板重新形成一个完整的承载结构。
4.7 现场控制效果分析
11205运输下山充填支护承载结构工程于2021年1月实施,至今已经将近3年时间(图22)。为验证再造承载结构设计的合理性,对该再造承载结构工程进行了相关的监测工作,采用了常规的监测方式对围岩的收敛进行定期测量。根据现场情况,一共进行了4个测点(p1,p2,p3,p4)的监测,求其平均值得到最终的监测数据,如图23所示。
从图23中可以看到,巷旁充填墙再造承载结构支护后的巷道经过120 d的监测后,总体变形经历了加速期、减速期及稳定期,其监测断面的收敛速率基本都小于0.2 mm/d,表面围岩变形进入了减速阶段,且各收敛速率具有随着监测时间的增长而减小的趋势并趋于稳定,有利于巷道的较长时期的稳定。对比数值模拟实验发现,监测点p1、p3出均出现较大的变形,且数值模拟试验结果(25 mm)与工程实践监测(20 mm)结果较为吻合,其次p2点位移较大,而p4点位移最小。数值模拟实验结果与工程实践监测变形在同一数量级,且最大变形量值较为接近,变化趋势两者均吻合,表明数值模拟实验的合理性及有效性。
5. 结 论
1)在复杂地质条件下,传统支护理论与技术难以根本克服巷道围岩大变形问题。针对巷道多样变形破坏机制,深入分析现有技术原理,创新提出系列再造承载结构拱支护技术与工艺体系,结合主被动支护原理,利用再造承载结构拱技术,构建具有高强承载能力的结构体,从整体上调动围岩体的自身承载能力与围岩整体抗变形能力,发挥整体承载和整体抵抗变形的能力,保证巷道正常使用,为复杂条件下煤矿巷道支护提供参考。
2)以龙宝煤矿为研究对象,通过调研分析了巷道围岩大变形、支护构件失效的原因,采用Python脚本编程语言实现对FLAC3D数值模拟计算软件的控制,结合飞蛾火焰优化算法、巷旁充填墙优化目标函数及FLAC3D数值模拟计算,确定巷旁充填墙几何最优参数(边宽0.49 m,巷旁充填墙厚度1.45 m),在确保巷道稳定的前提下,获取最优的巷旁充填墙位置及巷旁充填墙厚度。
3)推导了巷旁充填墙承重强度公式,结合龙宝煤矿11205运输下山施工,计算得到巷旁充填墙承载强度q值为39.78 MPa,这与数值模拟计算结果(36.098 MPa)吻合较好,且两者之间的误差为9.25%,该误差在可接受的范围之内。
4)基于研究成果确定对龙宝煤矿软弱围岩通过再造高强度墙体进行置换,井下工业实验表明,充填墙再造承载结构支护后的巷道经 过120 d的监测发现,总体变形经历了加速期、减速期及稳定期,其监测断面的收敛速率基本都小于0.2 mm/d,无明显变形,验证了支护的有效性及合理性。
6. 展 望
在工程实际中,应根据不同巷道围岩破坏特征,结合高强度高韧性充填新材料开发不同形式的再造承载结构快速成巷技术还需继续深入研究和验证。其中充填材料及设备还有提高的空间,再造承载结构与围岩体强度、刚度的可适性及有效性还需进一步优化,同时需不断补充高精度、多源化的矿压监测,从而形成巷道再造承载体围岩结构失效判别的新方法,为形成一套集技术、工艺、材料到装备制造的全产业链技术体系进一步丰富和完善,则有望解决深部煤矿巷道松软岩层的稳定控制技术难题,提高煤矿巷道整体稳定性,主要包括以下方面:
1)多次冲击荷载作用下巷道围岩变形破坏机理研究。基于围岩基本物理力学性质基础上,研究多次冲击荷载的特点,包括冲击荷载的大小、频率、冲击形式等因素,研究多次冲击荷载下巷道围岩的变形规律以及不同冲击荷载作用下围岩的变形特征。
2)多次冲击荷载下巷道稳定性评价及预测方法研究。通过基本物理力学试验、现场实测以及数值模拟的分析,在探究围岩变形规律的基础上,进一步评价多次冲击荷载下巷道的稳定性。主要是对巷道变形和破坏的临界条件进行研究,考虑不同冲击荷载条件下的巷道稳定性,研究巷道围岩、支护结构等因素对巷道围岩稳定性的影响,通过稳定性评价,提出预测模型及方法,与现场实测进行对比分析,验证预测模型的准确性。
3)巷道再造承载结构与围岩、注浆加固、锚网索支护等组合条件下的耦合机理和协同支护效果还需进一步研究。进一步研究再造承载结构力学模型与搭建实验室模拟试验平台,开展等比例再造承载结构相似模拟试验,分析巷道围岩−支护结构的耦合效应,形成巷道再造承载结构围岩结构失效判别的新方法;研究再造承载拱强度和刚度随时间变化规律、再造承载拱与围岩体刚度的可适性及有效性。
4)进一步研发低成本、更加高强轻质的新型支护材料。同时考虑到井下巷道空间狭小、机械设备较多的现实情况,研发喷射等更加简单的施工工艺,研制更加实用的施工装备。解决再造结构在边界不变化的规则几何截面下的快速、高效率的施工,提高施工经济性和安全性,实现各个工序之间的平行问题,实现快速成巷。本文只提出了一些思路,所做工作也还需进一步加强,因此难免有考虑不全和不妥之处,还请大家多多包含并给予批评指正为谢!
-
表 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.