城市地下综合管廊通风传热模型
Ventilation and heat transfer modeling in urban utility tunnel
收稿日期: 2020-10-22
基金资助: |
|
Received: 2020-10-22
Fund supported: | 国家自然科学基金资助项目(51678518) |
作者简介 About authors
何国青(1975—),男,副教授,博士,从事建筑环境和节能研究.orcid.org/0000-0002-7667-2335.E-mail:
为了优化地下管廊的通风设计,降低运行费用,建立苏州某地下管廊(UUT)一典型段的通风传热的解析模型. 考虑舱间传热、空气和廊壁的耦合传热以及土壤边界厚度的影响. 通过与计算流体动力学(CFD)模型的比较,验证解析模型的可靠性. 结果表明,廊内热源、通风换热和向土壤传热是影响舱内温度的主要因素,舱间传热的影响相对较小. 考虑降温需求的各舱室最低通风量受热源强度、进风温度(大气温度)和土壤温度影响较大,实际管廊最低通风量的计算须结合当地气象条件和管廊热源强度. 对于苏州地区地下综合管廊试验段,满负荷运行的电力舱和热力舱的最低通风换气次数分别为10、5次/h,而无热源的水信舱和燃气舱均可以按现行规范要求的最低换气次数设计通风量. 对于其他布置形式的管廊,可以采用类似的解析模型进行分析.
关键词:
An analytical model for ventilation and heat transfer in one typical urban utility tunnel (UUT) section in Suzhou was developed in order to aid the design and optimization of the ventilation system of the UUT for reduction in ventilation energy. The model includes the heat transfer among neighboring cabins, the conjugate heat transfer at the tunnel walls, and the thickness of the overlaying soil layer. The analytical model was validated by comparing it with the computational fluid dynamics (CFD) model. Results show that the internal heat source, ventilation, and soil temperature are the most important factors on cabin temperature, with the inter-cabin heat transfer a minor factor. The minimum ventilation rate to prevent over-heating of cabins is strongly dependent on the heat source strength, the inlet temperature (atmospheric temperature), and the soil temperature. Therefore, the minimum ventilation rate of UUT in other areas should be calculated according to the local meteorological conditions and the actual heat source intensity of the tunnel. The minimum ventilation rates, in air change rate, of the fully loaded cable cabin and the steam pipe cabin were 10 h−1 and 5 h−1, respectively, in the test section of UUT in Suzhou area. The minimum ventilation rate specified in the current code for UUT can meet the temperature control requirement for cabins without any heat sources. A similar analytical modelling approach can be used to analyze the ventilation requirement for UUT with other layouts of cabins.
Keywords:
本文引用格式
何国青, 赵文杰, 王亮.
HE Guo-qing, ZHAO Wen-jie, WANG Liang.
为了维持廊内空气品质,须保证足够的通风换气量[13]. 通风系统使用了大量的风机,能耗较大,有较大的节能空间,因此通风系统的正确选择关系到管廊运营的安全性和经济性[14]. 我国城市综合管廊工程技术规范[15]规定,正常通风换气次数不应小于2次/h,事故通风换气次数不应小于6次/h. 天然气管道舱正常通风换气次数不应小于6次/h,事故通风换气次数不应小于12次/h[16]. 规范并未针对管廊温度的控制规定通风换气次数,但要求廊内温度不高于40 ℃. 因此,对于有较强稳定内热源的舱,比如投入运行的电力舱和热力舱,还须考虑维持廊内热环境所需的通风量. 若热源强度大,2次/h的换气次数未必能满足廊内温度低于40 ℃的要求. 如果换气量过大,通风能耗势必会显著增加,因此,正确估算通风量对管廊节能运行的意义重大.
近年来,有关管廊通风的研究逐渐增多. 唐宏辉等[17]对广州某管廊的通风阻力进行测量和讨论,王雪梅等[18]模拟研究热力舱通风下的温度分布,闵绚等[19-20]模拟研究通风对电力舱降温的效果,王振芹等[21-22]模拟研究电力舱的温度场,刘承东等[23]用数值模拟的方法研究用缩尺模型代替原尺寸模型的可行性. 也有研究使用解析模型来为管廊通风提供设计依据的,如文献[24]、[25]. 计算流体动力学(computational fluid dynamics,CFD)模拟技术是成熟的空气流动和传热计算技术,可以研究管廊通风效果的细节,了解温度分布,但其需要迭代计算,比较耗时,难以满足管廊通风优化设计的快速计算需求. 相比较而言,解析模型可以快速计算,适合设计计算. 目前,在地下管廊的通风传热研究中,对管廊壁边界的处理并没有统一的方法,有的文献[18, 22]采用绝热边界,有的文献[21, 24]采用定温边界,尤其对土壤层的影响缺乏考虑. 文献[25]中的解析模型也采用了单舱模型和简化计算方法,既没有考虑空气对流和固体传热的耦合,也没有考虑隔壁非发热舱状态的影响. 针对上述问题,本研究建立了苏州综合管廊的传热模型,不仅考虑了发热舱和非发热舱之间的传热,而且考虑了空气对流传热和廊壁面传热的耦合,以及土壤边界层的影响,通过和CFD模拟结果的对比,证实解析模型的可靠性.
1. 解析模型
取苏州某综合地下管廊的一个标准管廊段作为研究对象,管廊截面如图1所示. 图中,从左到右分别为电力舱(宽2.4 m)、水信舱(宽2.8 m)、燃气舱(右上,宽2.1 m)和热力舱(右下,宽2.1 m). 廊内高度为4.05 m. 规范要求除燃气舱外各舱通风次数至少为2次/h,通风维持舱室出风温度在40 ℃以内. 假设散热所需的通风量为
管廊内表面和空气的传热系数分解为两部分:沿长度方向的强迫对流换热系数h1和自然对流换热系数h2. 本研究按定温度平板流动计算自然对流换热系数h2. 以下所采取的传热系数计算公式均来自文献[26].
图 1
图 1 苏州某地下管廊一标准段的传热模型
Fig.1 Heat transfer modeling of an underground utility tunnel in Suzhou
1.1. 管廊内强迫流动传热系数
沿管廊长度方向的气流流动一般已经达到湍流,管廊内平均流速为
式中:L为管廊长度,取200 m.
沿管长方向的流动,按定温平板流动计算强迫对流换热系数h1:
式中:
1.2. 圆管表面自然对流换热系数
对于热力舱,蒸汽圆管表面自然对流换热系数h2表达式如下:
式中:
如图2所示为蒸汽管道表面换热示意图,蒸汽钢管外径D0=0.377 m,保温材料为厚度为0.140 m的玻璃棉.
图 2
1.3. 各管廊壁面的自然对流换热系数
按定温度平板流动计算管廊壁面和空气的自然对流换热系数:
式中:
1.4. 蒸汽管热平衡方程
蒸汽管热平衡方程为
式中:tz为蒸汽管内的温度,根据设计温度取230 ℃;t为蒸汽管外壁面温度;Rq为沿着保温壁的径向传热热阻;
1.5. 舱室热平衡方程
对于舱室i,通风带走的热量和向各周边土壤或舱室传递的热量是平衡的:
式中:CP为比热容;V为舱室体积;ρ为舱内空气密度;tin为管廊入口温度;tout为管廊出口温度;Qc为除湿输入的冷量(正值)或热量(负值);ti为各舱室平均温度,i=d、r、q、s;j表示与舱室i传热的相邻舱室或土壤,tj为相邻舱室的平均温度或土壤温度,Aj为对应的传热面积,Rj为对应的传热热阻.
式中:Lk为围护结构中各层材料的厚度,λk为相应材料的导热系数,hk为从i舱到相邻舱室(或土壤)传热路径中包含的传热面上的对流换热系数,ti,out为舱室i的出口温度. 苏州管廊标准段的侧面有厚0.04 m的膨胀聚苯板护板,外侧壁混凝土厚度为0.35 m,水信舱和电力舱的隔墙厚度为0.20 m,水信舱和热力舱及燃气舱的隔墙厚度为0.25 m,地面混凝土层厚度为0.40 m,顶面混凝土厚度为0.35 m,热力舱和燃气舱隔板厚度为0.25 m.
1.6. 求解方法
式(8)、(9)包含了5个方程,将其重组得到如下函数:
方程组求解是对上述函数优化,使得f=0. 通过自编Fortran方程,采用穷举法搜索求解. 对各舱室温度(td、ts、tq、tr、t)设定搜索范围,对每一组温度计算f,数值最接近0的那一组温度为解. 可以根据精度要求,缩短搜索步长. 求解过程如下:1)依次选取一组舱室温度(td、ts、tq、tr、t). 2)对每一个舱的每一个传热面,迭代求解壁面温度:(a)假设壁面温度tw;(b)计算Nu=Nu(tw);(c)计算各表面对流换热系数h;(d)对每一个传热面,根据h计算传热量q1;(e)对同一个传热面,根据围护结构两侧的壁面温度计算传热量q2;(f)调整tw,重复步骤(c)~(e),直至∆q= q1−q2≈0;(g)获得所有传热面的壁面温度. 3)计算每个舱的净传热量,比如舱i:
2. CFD模拟
使用Phoenics软件对该管廊段进行CFD模拟仿真,计算管廊和土壤及空气的传热. 按如图1所示的截面构建模型,顶面为0.350 m厚的水泥,地面为0.550 m厚的水泥,两侧为0.500 m厚的水泥和0.040 m 厚的XPS保护板,蒸汽管保温材料厚度为0.14 m. 各舱室进口和出口都置于顶面,与舱同宽,长度为1.0 m,如图3(a)所示. 按平板表面流动,当换气次数n=2次/h时,距离管廊入口2 m处的雷诺数为104,因此,视廊内流动为湍流. 采用标准k-e湍流模型,自然对流采用常用的Boussinesq近似模型. 模型采用正交网格,网格数为1018400,在蒸汽管截面附近进行局部加密,以获得自然对流换热的计算精度,如图3(b)所示.
图 3
图 3 地下管廊的CFD模型以及网格示意图
Fig.3 CFD modeling of underground utility tunnel and mesh schematic diagram
在CFD模型中,在进口边界给定与解析模型计算相同的质量流量,即
3. 计算案例和结果
以苏州为例,如图4所示为某段管廊中水信舱一年的廊内温度、湿度(Rh)监测数据,廊内空气温度(Pt100,±0.5 ℃)在8月中旬达到最高,为26.6 ℃. 对管廊温度和管廊壁面温度(Pt1000,±0.1 ℃)的测量表明,空气温度非常接近廊壁温度,温差∆t中间值为1.0 ℃. 因此,为了计算方便,取土壤温度tsoil=25.6 ℃. 管廊夏季进风温度取当地空调设计室外计算温度ta=34.4 ℃. 各舱室通风量按各舱室的换气次数(燃气舱6次/h,其他舱室2次/h)以及进口空气参数计算. 蒸汽管保温材料(玻璃棉)导热系数取0.047 W/(m∙K),水泥导热系数取1.630 W/(m∙K),土壤导热系数取1.460 W/(m∙K),XPS板材导热系数取0.035 W/(m∙K),空气导热系数取0.027 W/(m∙K).
图 4
图 4 水信舱温度、湿度监测数据
Fig.4 Measured temperature and humidity in cabin for water pipes
3.1. 解析模型计算结果
管廊内平均流速按式(2)计算. 计算结果如图5所示. 图中,Qs为各舱和土壤间的换热量;Qa为通风带走的热量,负数表示管廊从空气获得热量. 当管廊仅维持基本的通风换气次数(燃气舱6次/h,其他舱2次/h)时,在所示燃气舱和电力舱100%投入使用后,电力舱的平均廊内温度达到47.2 ℃,出口温度超过60 ℃,无法满足规范中低于40 ℃的要求. 除了电力舱,热力舱的平均温度也达到42 ℃. 其他舱的廊内温度均在40 ℃以内. 注意到蒸汽管释放的热量约为18.4 kW,外壁面温度达到56 ℃,对流换热系数为3.34 W/(m2∙K),而各垂直壁面的自然对流换热系数约为0.80~1.94 W/(m2∙K),此换热系数远高于文献采用的0.27 W/(m2∙K)的土壤侧传热系数[24]. 这说明各壁面的实际对流换热系数差异较大,考虑到该值对管廊向土壤传热的影响较大,不宜用一个统一的值作为输入参数.
图 5
图 5 苏州某地下管廊段热量和温度计算结果
Fig.5 Calculated heat fluxes and temperatures in an underground utility tunnel in Suzhou
3.2. CFD模拟结果
如图6所示为管廊通风下的出口截面和电力舱中心剖面的温度t0分布. 电力舱温度最高,热力舱次之,两舱的平均温度均超过40 ℃. 值得注意的是电力舱廊内最高温度并不在出口处,而是在距出口约1/4管廊长度处的管廊顶部.
图 6
图 6 管廊通风下的出口截面和电力舱中心剖面温度分布图
Fig.6 Temperature contour at exit of ventilated tunnel and at center section of cable cabin
表 1 计算案例解析模型和CFD模型计算的热量和温度的比较
Tab.1
模型 | Qa/kW | Qs/kW | Qr/kW | td/℃ | tr/℃ | ts/℃ | tq/℃ |
解析模型 | 36.84 | 37.48 | 18.40 | 60.0 | 49.4 | 33.3 | 33.7 |
CFD模型 | 37.04 | 35.96 | 17.98 | 58.1 | 47.6 | 35.8 | 34.1 |
| −1% | 4% | 2% | 3% | 4% | −7% | −1% |
3.3. 苏州管廊通风分析
3.3.1. 影响廊内温度的因素
以计算案例为基准,研究热源强度、土壤温度、气温对舱室温度的影响. 如图7所示为电力舱最低通风量(满足出口温度≤40 ℃)和电缆发热量q的关系. 热源强度对舱内温度有决定性的影响. 当实际电缆发热量只有计算值的一半时,2次/h的通风量就可以将电力舱的温度控制在40 ℃以内.
图 7
图 7 电力舱最低通风换气次数和电缆发热量的关系
Fig.7 Relation between minimum ventilation rate and heat source strength in cable cabin
如图8所示,电力舱最低通风量(满足出口温度≤40 ℃)随气温和土壤温度的降低而降低.
图 8
图 8 电力舱最低通风量和进气温度、土壤温度的关系
Fig.8 Effects of inlet and soil temperatures on minimum ventilation rate for cable cabin
3.3.2. 计算案例下廊内温度
根据所给定计算案例,利用解析模型对苏州某标准管廊段进行通风需求分析. 为了维持舱内温度,有热源的舱室必须时刻通风,而非发热舱则不一定需要通风(比如未投入运行时考虑节能而不开启风机). 为了分析廊内温度和通风换气次数的关系,分析2种极端情况:非发热舱均按要求通风和非发热舱均不通风. 计算结果如图9所示. 可以看出,为了使电力舱和热力舱的温度控制在40 ℃以内,在非发热舱(水信舱和燃气舱)均不通风的情况下,热力舱的换气次数至少为5次/h,电力舱换气次数至少为10次/h. 当非发热舱按规范最低要求通风时,电力舱和热力舱的温度略有提高,但变化不大,为了使电力舱和热力舱的温度在40 ℃以内,热力舱和电力舱的换气次数分别至少为5、10次/h. 可以看出,非发热舱通风与否对热力舱和电力舱2个发热舱的温度影响可以忽略. 从另一个角度看,向土壤传热和通风仍是影响各舱温度的主要因素.
图 9
图 9 非发热舱不通风和通风状态下舱内温度和换气次数的关系
Fig.9 Effect of ventilation rates on cabin temperature with non-heated cabin unventilated and ventilated
4. 结 语
建立苏州标准地下管廊通风传热的解析模型,可以用于分析管廊通风设计. 通过与CFD模型比较,发现解析模型的解是可靠的. 结果显示,廊内热源、通风换热和土壤传热是影响舱内温度的主要因素,舱间传热的影响相对较小. 考虑降温需求的各舱室最低通风量受热源强度、进风温度、和土壤温度影响较大,因此实际管廊最低通风量的计算须结合当地的气象条件和管廊热源强度. 对于苏州地区地下综合管廊试验段,满负荷投入运行的电力舱和热力舱的最低通风换气次数分别为10、5次/h. 对于其他布置形式的管廊,可以采用类似文中的解析模型进行分析.
本研究假设1 m外的土壤温度不受管廊影响,这点仍有待验证. 下一步可以通过实测数据来确定这一土壤厚度,提高模型的可靠性.
参考文献
Effects of fire location on the capacity of smoke exhaust from natural ventilation shafts in urban tunnels
[J].DOI:10.1002/fam.2683 [本文引用: 1]
Low carbon effects of urban underground space
[J].
Feasibility study on the geothermal utility tunnel system
[J].DOI:10.1016/j.scs.2019.101445 [本文引用: 1]
哈尔滨地下空间开发利用现状与对策分析
[J].DOI:10.3321/j.issn:0367-6234.2007.10.028 [本文引用: 1]
The current situation and strategy for development and utility of underground space in Harbin
[J].DOI:10.3321/j.issn:0367-6234.2007.10.028 [本文引用: 1]
Indoor airborne fungal levels in selected comprehensive compartments of the urban utility tunnel in Nanjing, southeast China
[J].DOI:10.1016/j.scs.2019.101723 [本文引用: 1]
Indoor atmosphere hazard identification in person entry urban utility tunnels
[J].
Resolution of sick building syndrome in a high-security facility
[J].DOI:10.1080/10473220050075644 [本文引用: 1]
Positive associations between respiratory outcomes and fungal index in rural inhabitants of a representative sample of French dwellings
[J].DOI:10.1016/j.ijheh.2012.02.011 [本文引用: 1]
Study on the potential relationships between indoor culturable fungi, particle load and children respiratory health in Xi’an, China
[J].DOI:10.1016/j.buildenv.2014.05.029
Impact of a water-damaged indoor environment on kindergarten student absences due to upper respiratory infection
[J].
城市综合管廊防结露问题探讨
[J].
Discussion on problems of anti condensation in urban utility tunnel
[J].
城市综合管廊通风系统设计分析研究
[J].DOI:10.3969/j.issn.1003-0344.2017.05.022 [本文引用: 1]
Discussion and study of ventilation system design for urban utility tunnel
[J].DOI:10.3969/j.issn.1003-0344.2017.05.022 [本文引用: 1]
Ventilation and environmental control of underground spaces: a short review
[J].DOI:10.1051/e3sconf/201911101039 [本文引用: 1]
城市地下综合管廊通风设计探讨
[J].
Ventilation design of urban underground utility tunnel
[J].
某综合管廊通风系统实地测试及分析
[J].
Field test and analysis of ventilation system of a utility tunnel
[J].
综合管廊热力舱在机械通风模式下温度场的模拟分析
[J].
Simulation analysis of temperature field for thermal compartment in the utility tunnel
[J].
管线敷设与风机室布置对综合管廊通风阻力影响研究
[J].
Influence of conduit deployment and fan room location on ventilation flow resistance in a utility tunnel
[J].
综合管廊电缆舱通风数值模拟研究
[J].DOI:10.3969/j.issn.1003-0344.2016.11.007 [本文引用: 1]
Study on numerical simulation of utility cable tunnel ventilation system
[J].DOI:10.3969/j.issn.1003-0344.2016.11.007 [本文引用: 1]
GIl管廊通风温度场研究
[J].
Study on ventilation temperature field of Gil tunnel
[J].
综合管廊电缆舱室通风形式数值模拟研究
[J].DOI:10.3969/j.issn.1003-0344.2020.02.008 [本文引用: 2]
Simulation study on ventilation form of utility tunnel cable cabin
[J].DOI:10.3969/j.issn.1003-0344.2020.02.008 [本文引用: 2]
城市地下综合管廊通风缩尺模型的相似性分析
[J].
Similarity analysis on the ventilation reduced-scale model of urban underground utility tunnel
[J].
苏州城北路综合管廊通风系统设计
[J].DOI:10.3969/j.issn.1000-4416.2016.11.003 [本文引用: 4]
Design of ventilation system for Chengbei Road utility tunnel in Suzhou city
[J].DOI:10.3969/j.issn.1000-4416.2016.11.003 [本文引用: 4]
城市综合管廊通风系统设计
[J].
Ventilation system design of urban utility tunnel
[J].
/
〈 |
|
〉 |
