Transmission and Distribution Engineering and Technology 输配电工程与技术, 2012, 1, 7-12 http://dx.doi.org/10.12677/tdet.2012.12002 Published Online December 2012 (http://www.hanspub.org/journal/tdet.html) Longitudinal Temperature Distribution for 220 kV/3 kA CD HTS Cable* Zhiqiang Zheng 1, Y inshun Wa ng1, Xiao Li1, Yuan Cheng1, Changhui Dai1, Chi Xue2, Lianqi Zhao2 1State Key Laboratory for Alternate Electrical Power System with Renewable Energy Sources, Key Laboratory of HV and EMC Beijing, North China Electric Power University, Beijing 2Zhongtian Technology Group Company, Ltd., Shanghai Email: zhiqiang721@163.com, yswang@ncepu.edu.cn Received: Oct. 15th, 2012; revised: Nov. 2nd, 2012; accepted: Nov. 5th, 2012 Abstract: For long distance transmission HTS cables, co-flow cooling of LN2 should be adopted, since that counter- flow cooling of LN2 in cable former and cryostat may lead to large temperature difference between cable terminals. Due to their extraordinary property of mechanical strength and thermal compensation, corrugated pipes can be widely used in flowing channels of cooling liquid nitrogen (LN2), the coolant of superconductors. Hydraulic characteristics of LN2 in corrugated pipes are needed to be studied in design of high temperature superconducting (HTS) cable. Variation law of friction factors with Reynolds number for LN2 flow in corrugated pipes is systematically studied, and therefore value of pumping loss, which is caused by LN2 flowing in the corrugated pipe channels, for 220 kV/3 kA HTS cables can be obtained. Based on all items of losses of HTS cable calculated, equations and curve of temperature distribution for HTS cable can be finally obtained. Keywords: High Temperature Superconducting (HTS) Cable; Cold Dielectric (CD); Pumping Loss; Heat Transfer Characteristics; Temperature Distribution 220 kV/3 kA冷绝缘高温超导电缆的 轴向温度分布* 郑志强 1,王银顺 1,李 肖1,程 远1,代常会 1,薛 驰2,赵连岐 2 1华北电力大学高电压技术与电磁兼容北京市重点实验室,新能源电力系统国家重点实验室,北京 2中天科技集团股份有限公司,上海 Email: zhiqiang721@163.com, yswang@ncepu.edu.cn 收稿日期:2012 年10 月15日;修回日期:2012年11月2日;录用日期:2012 年11 月5日 摘 要:对远距离输电高温超导电缆而言,液氮在电缆导体骨架与低温容器内如采用反向循环流动制冷方式会 导致电缆两端温差过大,因此高温超导电缆应采用液氮同向流动的制冷方式。由于波纹管具有优良的机械强度 与热补偿特性,使其在高温超导电缆中得到了广泛应用,它可以形成冷却超导体用的制冷液氮的流通通道。在 高温超导电缆设计之前,有必要对液氮在波纹管中的流动特性进行研究。文章对波纹管通道内液氮流动的摩擦 系数随雷诺数的变化规律进行了研究分析,得到了220 kV/3 kA高温超导电缆液氮流动产生的泵浦损耗值。在 得到了高温超导电缆的各项损耗值之后,对高温超导电缆的温度分布进行了研究,并得到高温超导电缆温度分 布方程以及温度分布规律曲线图。 关键词:高温超导电缆;冷绝缘;泵浦损耗;换热特性;温度分布 *基金项目:国家自然科学基金项目(51077051)。 Copyright © 2012 Hanspub 7 220 kV/3 kA冷绝缘高温超导电缆的轴向温度分布 Copyright © 2012 Hanspub 8 1. 引言 的结构强度、长度、循环压力等。目前,对液氮在波 纹管内的流动特性的研究还很少,使得该研究具有重 要的意义。采用液氮双向流动的制冷方式会在高温超 导电缆两端产生较大温差[5],因此本文采用液氮同向 流动的冷却方式进行分析计算。考虑到运输条件、实 现技术等因素的影响,目前高温超导电缆总长度一般 都在 3000 m以内,所以本文中高温超导电缆的长度 取3000 m。论文通过对波纹管内液氮的流动特性进行 了系统地研究分析,得到了3000 m,220 kV/3 kA高 温超导电缆内液氮流动特性,并最终完成了电缆本体 的温度分布研究。 高温超导电缆因其损耗低、载流能力强和电流密 度大等优点,是未来实现大容量、低损耗的主输电通 道之一[1]。按其绝缘介质工作温区的不同,可以将高 温超导电缆分为室温绝缘(Warm Dielectric)电缆和冷 绝缘(Cold Dielectric)电缆两种类型。图1分别为两种 高温超导电缆结构示意图。冷绝缘高温超导电缆的运 行成本要远低于室温介质绝缘超导电缆,它可以更好 地输送大功率,输送功率为常规电缆五倍时损耗只有 常规电缆的2/3[2]。自 20 世纪90 年代以来,美国、日 本、中国和韩国等都相继开展高温超导电缆的研究[3]。 美国超导公司(AMSC)于2008 年研制成功并投入运行 了600 m,138 kV/2.4 kA 的高温超导电缆项目。这也 是迄今为止成功实现的运行容量最大的高温超导电 缆项目。本文基于华北电力大学与中天科技集团共同 合作进行的 220kV/3kA 冷绝缘高温超导电缆项目开 展研究。 2. 不同类型管道内液氮的流动特性 液氮通道两端压强差可由下式计算[6]: 2 2 h lV PD (1) 其中,λ为通道内液氮对应的摩擦系数,l为通道长度 (m),ρ为液氮密度(Kg/m3),V为液氮在通道内的流速 (m/s),Dh为液氮通道对应的水利学直径(m),其一般 定义为 由于波纹管具有比光滑管更好的热补偿性能,所 以波纹管在高温超导电缆中得到更广泛应用[4]。高温 超导电缆用波纹管的换热特性决定了高温超导电缆 外保护 铜屏蔽 PPLP绝缘 低温容器外壁 低温容器内壁 热绝缘 保护层 HTS导体层 骨架 液氮 液氮 4 h w A DP (2) 其中,A为液氮通道的横截面积(m2),Pw为液氮流通 部分的通道直径(m)。 液氮在通道内的流速: M V A (3) 其中,M为液氮的质量流速(Kg/s)。 (a) 通道内液氮的摩擦系数λ对于不同管道类型有不 同的计算方法,下面以光滑管(如图 2所示)与波纹管 (如图 3所示)为例进行分析: 骨架 HTS 导体层 PPLP绝缘层 HTS屏蔽层 保护层 液氮 低温容器内壁 低温容器外壁 外保护 铜屏蔽/稳定层 液氮 (a) 光滑管的 Koo 摩擦系数按下式计算[7]: 0.32 0.00140.125 Re sm , (4) (b) Blasius 等效摩擦系数按下式计算[8]: 0.2 Blasius 0.096 Re (5) (b) (c) 波纹管 Hex 摩擦系数按下式计算[9]: Figure 1. Schematic view of HTS cable (a) WD (b) CD 0.2 Hex 0.451 Re (6) 图1. 高温超导电缆结构 (a)室温绝缘 (b)冷绝缘 220 kV/3 kA冷绝缘高温超导电缆的轴向温度分布 Figure 2. Schematic view of smooth pipe 图2. 光滑管示意图 S D LN 2 d LN2 Figure 3. Schematic view of corrugated pipe 图3. 波纹管示意图 (d) 波纹管摩擦系数在雷诺数(Re)为10 4左右时 与只与液氮通道几何参数有关,计算公式为[10]: 2 2 1 40.438 cor dd sd s (7) (e) 波纹管的摩擦系数与雷诺数以及管道几何结 构的关系通式如下[11]: 0.179 0.164 0.239 1.15Re corh h tD sD (8) 其中,t,s分别为波纹管波纹高度与波纹间距(m),d 为波纹管的内径(m),Re为雷诺数,其计算公式如下: Re Vd (9) 其中,μ为液氮的动力粘度(Pa·s)。 由文献[12]知道,波纹管内摩擦系数为光滑管内 摩擦系数的4倍左右。 在得到了波纹管摩擦系数的条件下,结合高温超 导电缆结构参数,即可对液氮流动产生的泵浦损耗 (Pumping loss)进行计算。 对流换热系数 h是表征流体与壁面对流换热的强 弱程度的一个重要参数,其计算方法为: h hNukD (10) 其中,k为流体热导率(W/(m·K)),Nu为努塞特(Nusselt) 数,其计算方法为: 0.8 0.4 0.023Re PrNu (11) 其中,Pr为普朗特(Prandtl)数,其计算方法为: Pr p Ck (12) 其中,Cp为流体比热容(J/(Kg·K)) 。 3. 220 kV/3 kA高温超导电缆泵浦损耗 220 kV/3 kA冷绝缘高温超导电缆内液氮流通通 道参数如下表1所示。 结合 220 kV/3 kA高温超导电缆的结构参数,波 纹管内液氮流动的摩擦系数与雷诺数的关系如下图 4 所示。其中,摩擦系数下标1、2分别表示的是骨架 通道与低温容器和电缆本体形成的“环形”通道。 在雷诺数为 104附近变化时,液氮流动特性逐渐 由层流(Laminar flow)特性过渡为湍流(Turbulent flow) 特性[8]。由波纹管内液氮流动的摩擦系数与雷诺数的 关系曲线可以看出,在雷诺数介于3 × 103~105之间 时,λHex 总比相同雷诺数下的其他波纹管摩擦系数计 算值要大,因此液氮在高温超导电缆波纹管内流动的 摩擦系数可以按公式(6)进行计算。 Table 1. Parameters of LN2 flowing channels of 220 kV/3 kA HTS cable 表1. 220 kV/3 kA高温超导电缆内液氮流通通道参数 液氮 通道 内径 d (mm) 外径 d0 (mm) 波纹间距 s (mm) 波纹高度 t (mm) 通道长 度l (m) 骨架 通道 20 - 3 1 3000 “环形” 通道 80 90 3 1 3000 Figure 4. Relation of friction factor λ with reynolds number Re 图4. 摩擦系数 λ与雷诺数(Re)的关系 Copyright © 2012 Hanspub 9 220 kV/3 kA冷绝缘高温超导电缆的轴向温度分布 对于液氮同向流动的高温超导电缆,其骨架通道 与低温容器“环形”通道两端压差应大小相等,即满 足 。结合式(1)、(3)、(9)可以得到: 1 PP 2 23 2 22 111 h h d MA MAd (13) 计算得到,β = 2.677。 液氮在标准大气压下的热力学参数如表 2所示。 选定不同质量流速下,电缆两侧压强差与电缆(冷 却通道)长度的关系曲线如图 5所示。 选用的与两缸斯特林(Stirling) 制冷机配套使 用的 液氮泵可以提供的液氮流量为900 L/h,即通入高温超 导电缆的液氮总质量流速为M = 0.26 Kg/s。 3.1. 骨架通道内液氮的泵浦损耗 在电缆铜骨架所形成通道内,结合式(13)可得液 氮流速 V1 = 0.2813 m/s。 对应的液氮流动雷诺数Re1 = 3.0007 × 104,对应 的波纹管骨架通道内液氮流动的摩擦系数 λ1 = 0.0574。 高温超导电缆波纹管骨架通道两端压强差:ΔP1 = Table 2. Thermo-physical properties of LN2 (at 77 K, 0.1 MPa) 表2. 液氮在 77 K,0.1 MPa下的热力学参数 比热容, Cp (J/(Kg·K)) 热导率, k (W/m/K) 动力粘度, µ (µPa) 普朗特数, Pr 密度, ρ (Kg/m3) 2 × 103 0.15 150 2 800 050010001500 2000 25003000 0 0.5 1 1.5 2 2.5 3 3.5 4x 10 5 m Pa M=0.20 M=0.22 M=0.24 M=0.26 M=0.28 M=0.30 Figure 5. Relation of pressure difference with length of HTS cable 图5. 高温超导电缆两侧压强差与电缆长度的关系 0.27245 MPa。 高温超导电缆两端的压强差与液氮流动产生的 泵浦损耗之间的关系如下: pump QVAP (14) 因此,单位长度泵浦损耗为Qpump1 = 0.0080 W/m。 由式(11)可得,努塞特数Nu1 = 115.8573。 因此,波纹管骨架通道内对应的对流换热系数h1 = 868.929W/(m2·K)。 3.2. 低温容器内壁与电缆本体形成的“环形” 通道内液氮的泵浦损耗 低温容器内壁与电缆绝缘层之间形成的“环形” 通道横截面积A2 = 0.0013 m2,流速V2 = 0.1772 m/s。 液氮流动雷诺数Re2 = 9.4518 × 103,对应的“环 形”通道内液氮流动的摩擦系数λ2 = 0.0723。 高温超导电缆低温容器内壁与电缆本体形成的 “环形”通道两端压强差 ΔP2 = 0.27245 MPa,计算结 果与 ΔP1相等,即两通道两端压强差相等。 单位长度“环形”通道内的液氮泵浦损耗Qpump2 = 0.0215 W/m。 由式(11)可得,努塞特数:Nu2 = 45.9780。 因此电缆低温容器内壁与电缆本体形成的“环 形”通道对应的对流换热系数 h2 = 689.6699 W/ (m2·K)。 综上,220 kV/3 kA冷绝缘高温超导电缆同向液 氮流动产生的泵浦损耗为: pumppump1pump2 0.0295 QQQ Wm 4. 低温容器漏热与电缆介质损耗 4.1. 低温容器漏热计算 根据经验公式[13] 0 2π ln eff i co ci kTT QDD (15) 其中,keff为低温绝热材料对应的有效热导率(W/ (m·K)),Ti、T∞分别为低温容器内、外壁温度(K),Dio 、 Dco 分别为低温容器内、外壁直径(m)。低温容器绝热 为高真空多层绝热,绝热材料采用10 层GS-80 绝热 材料。 Copyright © 2012 Hanspub 10 220 kV/3 kA冷绝缘高温超导电缆的轴向温度分布 Copyright © 2012 Hanspub 11 Table 3. Electrical dissipation parameter of 220 kV/3 kA HTS cable 计算可得高温超导电缆低温容器漏热: 表3. 220 kV/3 kA超导电缆电气损耗参数 cryo 0.1732 QWm 磁滞损耗 (W/m) 带材涡流 损耗(W/m) 骨架涡流 损耗(W/m) 介质(PPLP) 损耗(W/m) 电气总损 耗(W/m) 1.442 0.00616 0.0214 0.0364 1.5055 4.2. 电缆介质损耗 绝缘材料在交变电磁场中产生介质损耗,一般采 用下式计算: 20 ln r o iiio QU tg rr (16) 其中 Uo为高温超导电缆对地电压, 为圆频率,rii 和 rio 为绝缘层内外半径, r为绝缘层相对介电常数,tg 为绝缘介质损耗角。对于PPLP 绝缘材料,液氮温度 下相对介电常数为2.21,介质损耗正切为8 10–4。 计算可得,高温超导电缆介质损耗 Qδ = 0.0364 W/m。 Figure 6. Heat balance of HTS cable 5. 高温超导电缆总体损耗 图6. 高温超导电缆热平衡示意图 高温超导电缆电气损耗包括电缆磁滞损耗、带材 涡流损耗、铜骨架涡流损耗以及介质损耗等。对于液 氮泵本身产生的损耗,因其在高温超导电缆末端,其 温度变化对电缆无影响,因此可忽略不计。结合之前 完成的 220 kV/3 kA冷绝缘高温超导电缆电气损耗计 算结果(计算成果已整理成文并已投其他期刊,尚未发 表),额定运行情况下高温超导电缆各种电气损耗参数 如下表 3所示。 1 1 2 10 d2 d d2 d p p T CMqx Q x T CMqx QQ x (17) 边界条件如下: 12 00TT 0 T 其中 电缆总电气损耗Qelec = 1.5055 W/m。 21 2πqxKTxT x (18) 综上,220 kV/3 kA高温超导电缆总损耗主要包 括电缆电气损耗、液压损耗以及低温容器漏热等,即: 21 112 2 ln 11 1 rr K hr hr (19) Qt = Qpump + Qcryo + Qelec = 1.7082 W/m,折合成每 kA 电流下电缆损耗为0.5694 W/(m·kA)。 式中,r为液氮通道的半径,h为热传导系数,λ为介 质热导率,本文中即 PPLP 的热导率,由文献[14,15] 可知,λ可分别取值为 0.2 W/m/K 与0.05 W/m/K。 6. 高温超导电缆轴向温度分布 由热平衡方程可以得到骨架通道与低温容器和 缆本体形成的“环形”通道的温度分布方程: 高温超导电缆处于稳态工作时的热平衡示意图 如图 6所示,结合稳态下的热平衡微分方程组: 电 1 1 2π12 00 0 00 10 2 2 1 1 2π1 00 00 0 20 2 2 1 22 2π112π1 12π12π 22e 2π112π1 12π12π p p Kx CM p p Kx CM p Q QQQQQQQ QQQ Q Tx xTex CM KKCM KK Q QQQQQQQ QQ Tx Tx KKCM KK 0 其中, , elec pump1 2QQ Q0cryo pump2 pump1 QQ QQ 各项损耗值已在上述介绍中给出,可以得到 λ = 220 kV/3 kA冷绝缘高温超导电缆的轴向温度分布 0.05与λ = 0.2时的高温超导电缆温度分布曲线如图7 所示。 由温度分布图可以看出,相对于液氮反向流动的 制冷方式, 3000 m 高温超导电缆两端温差可以维持在 较小的(10 K左右)范围内,这也验证了高温超导电缆 采用液氮同向流动制冷方式的可行性;对于绝缘材料 PPLP的热导率取值较大时对应的高温超导电缆两端 温差较小,这也说明了高温超导电缆两端温差随介质 热导率的增大而减小的规律。需要指出的是,在满足 绝缘强度的前提下,对于长度为大于 3000 m、采用不 同热导率绝缘材料的高温超导电缆温度分布也可以 进行计算。 0500 10001500 2000 2500 3000 64 66 68 70 72 74 76 m K 1 2 (a) 050010001500 2000 25003000 65 66 67 68 69 70 71 72 73 74 75 m K 1 2 (b) Figure 7. Longitudinal temperature distribution in HTS cable (a) λ = 0.05 (b) λ = 0.2 图7. 高温超导电缆轴向温度分布 7. 结论 综合分析了波纹管摩擦系数随雷诺数变化的规 律,得到了220 kV/3kA高温超导电缆波纹管内液氮 流动的摩擦系数计算方法,并最终计算得到了液氮同 向流动情形下高温超导电缆通道内液氮的泵浦损耗 值0.0295 W/m。然后结合高温超导电缆的电气损耗与 低温容器漏热计算结果,建立了高温超导电缆温度分 布方程,并得到温度分布曲线。从方程以及曲线中可 以看出,随热导率增大电缆两端温差有减小的趋势; 同时也验证了液氮同向流动方式下电缆两端温差较 小的结论,对实用高温超导电缆温度分布研究有重要 的实际参考价值。 参考文献 (References) [1] 王银顺. 超导电力技术基础[M]. 北京: 科学出版社, 2011. [2] 林良真. 我国超导技术研究现状及展望[J]. 电工技术学报, 2005, 20(1): 1-7. [3] 林良真, 肖立业. 高温超导输电电缆的现状与发 展[J]. 电力 设备, 2007, 8(1): 1-4. [4] 孙凤玉, 徐烈等. 高温超导电缆用波纹管内液氮流动特性的 数值研究[J]. 低温技术, 2007, 35(5): 376-380. [5] S. Fuchino, et al. Longitudinal temperature distribution in super- conducting power cables with counter-flow cooling. IEEE Transaction on Applied Superconductivity, 2002, 12(1): 1339- 1342. [6] 张英主. 工程流体力学[M]. 北京: 中国水利水电出版社, 2002: 69-77. [7] T. S. Drew, E. C. Koo and W. H. McAdams. The friction factor in clean, round pipes. Transactions of the American Institute of Chemical Engineers, 1932, 28: 56. [8] S. Fuchino, et al. Hydraulic characteristics in superconducting power transmission cables. Physica C, 2001, 354: 125-128. [9] JSME. JSME data book: Heat transfer. 4th Edition, Tokyo: Ma- ruzen, 1991: 255. [10] F. D. Yeaple, et al. Hydraulic and pneumatic power and control. McGraw-Hill, New York, 1966. [11] S. Pethkool, et al. Turbulent heat transfer enhancement in a heat exchanger using helically corrugated tube. International Com- munications in Heat and Mass Transfer, 2011, 38: 340-347. [12] J. G. Weisend and S.W. Van Sciver. Pressure drop from flow of cryogens in corrugated bellows. Cryogenics, 1990, 30: 935-941. [13] J. A. Demko, et al. Practical AC loss and thermal considerations for HTS power transmission cable systems. IEEE Transactions on Applied Superconductivity, 2001, 11(1): 1789-1792. [14] S. Akita, K. Ishihara, J. Ogata and M. Nemoto. Development of external cooling cryoresistive cable system. Komae Research Laboratory Report No. T86081, Nihon Kikai Gakkai, 1987. [15] T. Masuda et al. Study of a 66kV 3-core high-T superconducting cable. IEEE Japanese Proceedings of 11th Annual Conference on Power & Energy Society, 2000: 592-593. Copyright © 2012 Hanspub 12 |