![]()
Transmission and Distribution Engineering and Technology
输配电工程与技术, 2012, 1, 7-12
http://dx.doi.org/10.12677/tdet.2012.12002
Published Online December 2012 (//www.abtbus.com/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
1
State Key Laboratory for Alternate Electrical Power System with Renewable Energy Sources, Key Laboratory of HV and EMC Beijing,
North China Electric Power University, Beijing
2
Zhongtian Technology Group Company, Ltd., Shanghai
Email: zhiqiang721@163.com, yswang@ncepu.edu.cn
Received: Oct. 15
th
, 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 LN
2
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 (LN
2
), 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 LN
2
flow in corrugated pipes is systematically studied, and therefore value
of pumping loss, which is caused by LN
2
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
P
D
(1)
其中,
λ为通道内液氮对应的摩擦系数,l为通道长度
(m)
,ρ为液氮密度(Kg/m3),V为液氮在通道内的流速
(m/s)
,Dh为液氮通道对应的水利学直径(m),其一般
定义为
由于波纹
管具有比光滑管更好的热补偿性能,所
以波纹管在高温超导电缆中得到更广泛应用
[4]。高温
超导电缆用波纹管的换热特性决定了高温超导电缆
外保护
铜屏蔽
PPLP
绝缘
低温容器外壁
低温容器内壁
热绝缘
保护层
HTS
导体层
骨架
液氮
液氮
4
h
w
A
D
P
(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
LN
2
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 Pr
Nu (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)
外径
d
0
(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)的关系
Cop
yright © 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 LN
2 (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 × 10
3 0.15 150 2 800
0
50010001500 2000 25003000
0
0.5
1
1.5
2
2.5
3
3.5
4
x 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/(m
2
·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/
(m
2
·K)。
综上,
220 kV/3 kA冷绝缘高温超导电缆同向液
氮流动产生的泵浦损耗为:
pumppump1pump2
0.0295 QQQ Wm
4.
低温容器漏热与电缆介质损耗
4.1.
低温容器漏热计算
根据经验公式
[13]
0
2
π
ln
eff i
co ci
kTT
Q
DD
(15)
其中,
keff为低温绝热材料对应的有效热导率
(W/
(m·K))
,Ti、T∞分别为低温容器内、外壁温度(K),Dio
、
D
co
分别为低温容器内、外壁直径(m)。低温容器绝热
为高真空多层绝热,绝热材料采用
10
层GS-80 绝热
材料。
Cop
yright © 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.
电缆介质损耗
绝缘材料在交变电磁场中产生介质损耗,一般采
用下式计算:
2
0
ln
r
o
iiio
QU tg
rr
(16)
其中
Uo为高温超导电缆对地电压,
为圆频率,
rii 和
r
io
为绝缘层内外半径,
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
d
2
d
d
2
d
p
p
T
CMqx Q
x
T
CMqx QQ
x
(17)
边界条件如下:
12
00
TT 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)
Q
t
= 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
2
2
2
π
112π1
12
π12π
2
2
e
2
π
112π1
12
π12π
p
p
K
x
CM
p p
K
x
CM
p
Q QQQQQQQ
QQQ Q
Tx xTex
CM KKCM
KK
Q QQQQQQQ
QQ
Tx Tx
KKCM
KK
0
其中,
,
elec pump1
2
QQ Q0cryo pump2 pump1
QQ QQ
各项损耗值已在上述介绍中给出,可以得到 λ =
![]()
220 kV/3 kA
冷绝缘高温超导电缆的轴向温度分布
0.05
与
λ
= 0.2
时的高温超导电缆温度分布曲线如图7
所示。
由温度分布图可以看出,相对于液氮反向流动的
制冷方式,
3000 m
高温超导电缆两端温差可以维持在
较小的
(10 K左右)范围内,这也验证了高温超导电缆
采用液氮同向流动制冷方式的可行性;对于绝缘材料
PPLP
的热导率取值较大时对应的高温超导电缆两端
温差较小,这也说明了高温超导电缆两端温差随介质
热导率的增大而减小的规律。需要指出的是,在满足
绝缘强度的前提下,对于长度为大于
3000 m、采用不
同热导率绝缘材料的高温超导电缆温度分布也可以
进行计算。
0
500 10001500 2000 2500 3000
64
66
68
70
72
74
76
m
K
1
2
(a)
0
50010001500 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.
Cop
yright © 2012 Hanspub
12
|