不同直流拓扑和敷设环境下35 kV交流XLPE电缆直流载流量分析

张添胤; 王启隆; 王新; 金泱; 杨敏; 陈向荣

doi:10.3969/j.issn.0258-2724.20220704

不同直流拓扑和敷设环境下35 kV交流XLPE电缆直流载流量分析

doi: 10.3969/j.issn.0258-2724.20220704

张添胤^{1, 2,},
王启隆^{1, 2},
王新³,
金泱³,
杨敏³,
陈向荣^{1, 2, 4, 5, ,}

1.
浙江大学电气工程学院，浙江杭州 310027
2.
浙江大学杭州国际科创中心，浙江杭州 311200
3.
浙江浙能技术研究院有限公司浙江省火力发电高效节能与污染物控制技术研究重点实验室，浙江杭州 311121
4.
浙江大学浙江省宽禁带功率半导体材料与器件重点实验室，浙江杭州 311200
5.
浙江大学先进电气国际研究中心，浙江海宁 314400

详细信息

作者简介:
张添胤（2000—），男，博士研究生，研究方向为高压直流输电和电力电缆技术，E-mail：zhangtianyin_zju@163.com

通讯作者:
陈向荣（1982—），男，研究员，研究方向为先进电气材料与高压绝缘测试技术、先进电力装备与新型电力系统、高电压新技术，E-mail：chenxiangrongxh@zju.edu.cn

中图分类号: TM247
计量
- 文章访问数: 186
- HTML全文浏览量: 57
- PDF下载量: 42
- 被引次数: 0
出版历程
- 收稿日期: 2022-10-19
- 修回日期: 2023-04-17
- 网络出版日期: 2024-08-31
- 刊出日期: 2023-05-08

Direct Current Ampacity of 35 kV Alternating Current Cross-Linked Polyethylene Cables Under Various Direct Current Topologies and Laying Environments

ZHANG Tianyin^{1, 2
,},
WANG Qilong^{1, 2},
WANG Xin³,
JIN Yang³,
YANG Min³,
CHEN Xiangrong^{1, 2, 4, 5
, ,}

1.
College of Electrical Engineering, Zhejiang University, Hangzhou 310027, China
2.
ZJU-Hangzhou Global Scientific and Technological Innovation Center, Hangzhou 311200, China
3.
Zhejiang Provincial Key Laboratory of High-Efficiency Energy-Saving and Pollutant Control Technology for Thermal Power Generation, Zhejiang Energy Research Institute Co., Ltd., Hangzhou 311121, China
4.
Zhejiang Provincial Key Laboratory of Power Semiconductor Materials and Devices, Zhejiang University, Hangzhou 311200, China
5.
International Research Center for Advanced Electrical Engineering, Zhejiang University, Haining 314400, China

摘要

摘要:
交流电缆改为直流运行对于新能源发电并网以及提高供电容量有着重要意义. 以35 kV交流交联聚乙烯（cross-linked polyethylene，XLPE）电缆为研究对象，采用有限元方法对三线双极式、单极式和双极式3种直流拓扑结构下的三芯交流电缆进行温度场仿真；同时考虑配电网中常见的影响因素，如靠近热水管道、电缆集群敷设、电流不平衡等，研究其对交流电缆改为直流运行下的电缆载流量的影响. 研究发现：在相同工作温度下，35 kV交流XLPE电缆在单极式直流拓扑结构下运行时的载流量最小，在三线双极式直流拓扑下运行时的载流量最大；当交流电缆与城市供水管道水平间距为2.0 m，垂直间距为0.5 m时，在上述3种直流拓扑下运行时的载流量下降了约3.0%；在电缆集群敷设时，加载非同步峰值能使载流量最大提升90 A；电缆载流量随电流不平衡度先增大后减小，在不平衡度为0时，电缆载流量达到最大. 研究结果为35 kV交流XLPE电缆的直流改造提供了参考依据.
- XLPE电缆 /
- 直流拓扑 /
- 敷设环境 /
- 直流载流量
Abstract:
Converting alternating current (AC) cables into direct current (DC) operation is of great significance to achieve new energy generation connected to the grid and increase power supply capacity. A 35 kV AC cross-linked polyethylene (XLPE) cable was taken as the research object, and temperature field simulation of the three-core AC cable was carried out through the finite element method under three kinds of DC topologies: three-wire bipole, monopole, and bipole. At the same time, common influencing factors in the distribution network were considered, such as distance from hot water pipes, cable cluster laying, and current imbalance, so as to explore the influence of these factors on the ampacity of the cable when the AC cable was changed to the DC operation. The results show that under the same operating temperature, the 35 kV AC XLPE cable has the smallest ampacity when operating under the monopolar DC topology and the largest ampacity under the three-wire bipolar DC topology. When the horizontal distance between AC cable and urban water supply pipes is 2.0 m, and the vertical distance is 0.5 m, the ampacity of cables operating under the above kinds of DC topologies is reduced by about 3.0%. The maximum ampacity can be increased by 90 A as the asynchronous peak value is loaded when the cable cluster is laid. The cable ampacity increases first and then decreases with increasing current imbalance, and it reaches the maximum value when the current imbalance is 0. The research results can provide some reference for converting 35 kV AC XLPE cables into DC operation.
- XLPE cables /
- DC topologies /
- laying environment /
- DC ampacity

HTML全文

随着我国西部山区大跨桥梁蓬勃发展^[1-2]，复杂山区桥梁桥址区风场环境模拟日益受到重视^[3-4]. 而山区地形高低起伏，高度差异较大^[5]，在地形风场试验中不可避免会截断地形，截断地形的高度突变对来流特性将产生明显影响^[6-8]. 此时，在地形模型边界添加过渡段统一高低不平的地形边线是一种有效方法，所以，研究地形边界过渡段对复杂地形区桥梁风场研究具有重要意义^[9-11].

目前，国内外学者对地形模型边界过渡段已有一定研究：胡朋等^[12]推导出一类基于势流理论下圆柱绕流流线的地形模型过渡段曲线线型，并和斜坡过渡段对比，验证了其作为过渡段的适用性和优越性；靖洪淼等^[13]对基于双曲正切函数的过渡段流动特性进行了研究；Uchida等^[14]基于大涡模拟的方法对复杂地形上风场减速效应和流动分离进行了模拟研究.

然而，目前过渡段的设置原则及其合理形式仍未达成统一认识，绕兰金物体不可压缩无旋流动流场线及平方正弦曲线作为过渡段线型的可行性有待明确. 本文通过数值模拟的方式对比了上述2类线型与目前常用的双曲正切函数及平板风场线型的流动特性，并提出了地形模型边界过渡段线型构造的新原则，可以为地形桥址区的CFD （计算流体力学）模拟提供一定参考价值.

1. 过渡段设置原则

在桥址区风特性研究中，需先截取桥位处附近的一部分地形，后将其生成为几何固定壁面供后续的风洞试验或数值模拟使用. 截取地形常存在较大的起伏，致使地形模型与风洞地面或计算域底面存在高差，造成了不真实的台阶，如图1所示. 图中：x为水平方向位置，y为高程（L，H）为地形与过渡段曲线f(x)与地形平台的接触点坐标. 地形模型边缘台阶的存在使来流接触台阶，而后产生流动分离、局部空化现象和湍流转捩等与实际风场不符的现象，导致流场和实际来流产生明显不同. 该现象对地形范围较小、地形起伏大的桥址区风场结果影响尤其显著. 为此，设置过渡段来衔接不同高度的模型边界和风洞地板或数值模拟计算域底面，以实现来流从入口到地形区域的风场贴近实际情况.

图 1 地形模型边界

Figure 1. Terrain model boundary

下载: 全尺寸图片幻灯片

在桥址区地形风洞试验中常以风洞初始风场为参照，因此，试验中要求风洞来流由地面过渡到截取地形边缘时，其风场应尽量与参考风场保持相同的特性^[12-13]，避免发生来流与壁面分离.

存在逆压梯度和黏性力是来流接触地形边缘台阶后产生分离的原因. 以斜坡线型为例，流体动量方程如式（1）.

$\begin{array}{c} \dfrac{{\rm{d}}{U}_{\rm{i}}}{{\rm{d}}t}={f}_{\mathrm{i}}-\dfrac{1}{\rho }\dfrac{\partial p}{\partial {x}_{\mathrm{i}}} + \dfrac{\mu }{\rho }\dfrac{\partial }{\partial {y}_{{\rm{i}}}}\left(\dfrac{\partial {U}_{\mathrm{i}}}{\partial {y}_{{\rm{i}}}}\right)\text{，}\end{array}$

(1)

式中：U_i为流体速度，f_i为流体所受体积力，p为当前位置压强，ρ为流体密度，μ为黏性力系数， $x_{\rm{i}}$ 、 $y_{\rm{i}}$ 为 $x$ 、 $y$ 方向的2个位移分量，t为时间.

速度的突然变化将造成压力的突然变化，由此可以得出2种过渡段构造思路：

1）为使来流速度流经过渡段时平稳连续，可认为流线的过渡段线型应保持一阶导数连续. 类似地，考虑到黏性力的影响，流体加速度也应避免突变，过渡段线型应保持二阶导数连续. 由此，过渡段曲线f(x)需满足二阶导数连续和式（2）边界条件.

$\begin{array}{c}\left\{\begin{array}{c}f\left(0\right)=0,\;f\left(L\right)=H,\\ \dot{f}\left(0\right)=0,\;\dot{f}\left(L\right)=0,\\ \ddot{f}\left(0\right)=0,\;\ddot{f}\left(L\right)=0. \end{array}\right.\end{array}$

(2)

2）在势流场假定中各流线上的点满足流动不分离，各流线可视为彼此独立互不影响的气流，可认为气流在各流线上经过时性质不受影响，因此，可将流线线型作为过渡段线型. 图2所示为圆柱绕流流线. 图中：a为圆球半径. 通常选取已有解析表达式的流线，先考虑过渡段的边界条件，再进行截取，构造过渡段曲线.

图 2 理想流体圆柱绕流流线

Figure 2. Streamlines around an ideal fluid cylinder

下载: 全尺寸图片幻灯片

由上述2种构造过渡段思路出发，提出构造过渡段线型的2个原则：

1）过渡后的气流特性（如风速、风攻角和湍流强度等）应尽量与未受扰动的参考风场特性保持一致. 不能出现流动分离，因为气流分离会使边界层增厚、出现气流不稳定等现象.

2）气流从过渡段开始到达到稳定时的过渡总长度应尽可能短，较长的气流过渡总长度会进一步减小地形模型的有效区域，对于数值模拟方法，也会相应地增大计算区域尺寸、增加计算工作量.

2. 过渡段线型选取

目前，在地形模型风洞试验和CFD模拟中，主要的过渡段线型有：斜坡、理想流体圆柱绕流场流线、维多辛斯基曲线和双曲正切函数. 相较于其他主要过渡段线型，平方正弦曲线斜率变化均匀、一阶和二阶连续可导，对来流风场影响较小，符合上文提出的过渡段线型构造原则，是比较理想的过渡段线型. 为探究不同流线形式作为过渡段线型的适用性，本节利用平方正弦曲线分别使用第1节的思路1、2绘制典型过渡段线型并进行对比，研究其在作为过渡段时展现的气动特性. 结合上文提出的2个原则综合判断上述流线及其对应解析函数在地形模型过渡段问题上的适用性，本节选用双曲正切函数作为思路1的典型过渡段线型. 常见的具有解析式的流线线型：绕金兰物体的不可压缩无旋流动、均匀来流绕圆球流动等，选用绕兰金物体的不可压缩无旋流动函数作为思路2的过渡段线型.

2.1 平方正弦函数及双曲正切函数的过渡段曲线

用于过渡段曲线的平方正弦函数、双曲正切函数分别如式（3）、（4）.

$f_1\left(x\right)={h\mathrm{s}\mathrm{i}\mathrm{n}}^{2}\dfrac{{\text{π}} x}{2l}\text{，}$

(3)

$f_2\left(x\right)=\dfrac{h}{2}\left[\mathrm{t}\mathrm{a}\mathrm{n}\mathrm{h}\left(\dfrac{2cx}{l}-c\right) + 1\right]\text{，}$

(4)

式中：h为台阶高度差；l为过渡段长度；c为在tanh x中截取的单边长度，本节取4（一般取 $c\geqslant 3$ ）.

平方正弦函数、双曲正切函数及其一、二阶导数如图3所示. 由图可见，平方正弦函数及其一阶导数在边界处连续，满足式（2）边界条件的前2项；双曲正切函数满足式（2）中的全部条件.

图 3 用于过渡段的2种函数及其导函数

Figure 3. Two functions and their derivatives used for transition segments

下载: 全尺寸图片幻灯片

2.2 绕兰金物体过渡段曲线

绕兰金物体不可压缩无旋流动流函数为

$\begin{split} &\varPsi =-\frac{1}{2}{U}_{\mathrm{\infty }}{y}^{2} + \frac{1}{4{\text{π}} }\frac{Q\left(x + d\right)}{\sqrt{{\left(x + d\right)}^{2} + {y}^{2}}}-\\ &\quad \frac{1}{4{\text{π}} }\frac{Q\left(x-d\right)}{\sqrt{{\left(x-d\right)}^{2} + {y}^{2}}}, \end{split}$

(5)

式中： ${U}_{\infty }$ 为给定的来流速度，Q为源强，d为点源到点汇距离的一半.

令 ${b}^{2}=\dfrac{Q}{2{\text{π}} {U}_{\infty }}$ ，则有

$-\frac{2\varPsi }{{U}_{\mathrm{\infty }}}={y}^{2}-\frac{{b}^{2}\left(x + d\right)}{\sqrt{{\left(x + d\right)}^{2} + {y}^{2}}} + \frac{{b}^{2}\left(x-d\right)}{\sqrt{{\left(x-d\right)}^{2} + {y}^{2}}}.$

(6)

对流函数求导，可得其与速度场之间的关系为

$\begin{split} &{v}_{x}=\frac{1}{y}\frac{\partial {\varPsi }}{\partial y}=-\frac{{U}_{\mathrm{\infty }}}{2y}\times\\ &\quad \left(2y + \frac{{b}^{2}y\left(d-x\right)}{{\left({y}^{2} + {\left(d-x\right)}^{2}\right)}^{\tfrac{3}{2}}} + \frac{{b}^{2}y\left(d + x\right)}{{\left({\left(d + x\right)}^{2} + {y}^{2}\right)}^{\tfrac{3}{2}}}\right)\text{，} \end{split}$

(7)

$\begin{split} &{v}_{y}=-\frac{1}{x}\frac{\partial \varPsi }{\partial x}=\left(\frac{{b}^{2}}{\sqrt{{y}^{2} + {\left(d-x\right)}^{2}}}-\frac{{b}^{2}}{\sqrt{{y}^{2} + {\left(d + x\right)}^{2}}}-\right.\\ &\quad \left.\frac{{b}^{2}\left(d-x\right)\left(2b-2x\right)}{2{\left({y}^{2} + {\left(d-x\right)}^{2}\right)}^{\tfrac{3}{2}}} + \frac{{b}^{2}\left(d + x\right)\left(2b + 2x\right)}{2{\left({y}^{2} + {\left(d + x\right)}^{2}\right)}^{\tfrac{3}{2}}} \right) \frac{{U}_{\mathrm{\infty }}}{2y}, \end{split}$

(8)

式中：v_x、v_y分别为横向速度场和垂直速度场.

兰金物体半长轴长度R、半短轴长度S与b、d的关系为

$\left\{\begin{array}{l} \left({R}^{2}-{d}^{2}\right)=2{b}^{2}Rd,\\ {S}^{2}=\dfrac{2{b}^{2}d}{\sqrt{{S}^{2} + {d}^{2}}}. \end{array}\right.$

(9)

式（6）引入参数m可得流线方程通式为

$\begin{split} & E\left(x,y\right)={y}^{2}-\frac{{b}^{2}\left(x + d\right)}{\sqrt{{\left(x + d\right)}^{2} + {y}^{2}}} +\\ &\quad \frac{{b}^{2}\left(x-d\right)}{\sqrt{{\left(x-d\right)}^{2} + {y}^{2}}}-m=0. \end{split}$

(10)

由于过渡段曲线在其起点和终点处斜率为0，可得 ${x}_{1}=\pm \mathrm{\infty }$ ， ${x}_{2}=0$ ，起点不可能取到无穷远处，故取斜率较小的点作为起点. 当 $x=\pm \mathrm{\infty }$ 时，式（10）变为 ${y}^{2}-m=0$ ，可推得 $y=\pm \sqrt{m}$ 为流场线的渐近线.

在构造过渡段曲线时，通常已知来流边界和地形模型之间的高差即台阶高度差h和过渡段等效斜率l. 因此，不妨建立绕流兰金体的S、R和h、l之间的经验关系式，经式（9），求得b和d的值，再结合过渡段曲线应满足的边界条件确定出可作为过渡段的流线方程. 本节中取h=0.6 m，l=1/3，S=4/(3H)=0.8 m，R=h/(2l) + h/10=0.96 m，由式（9）联立方程组，解得b=0.745，d=0.565，将b、d代入式（10），绘制流场线，如图4所示.

图 4 绕兰金物体不可压缩无旋流动流场线

Figure 4. Streamline of incompressible irrotational flow around Rankin body

下载: 全尺寸图片幻灯片

求解方程 $y\left(0\right)-y\left(-\infty \right)=h$ ，可得m=0.04，绘制该流场线及其一阶、二阶导数线型如图5所示. 其中，流场线可作为地形边界过渡段曲线使用.

图 5 用于过渡段的流场线及其导函数

Figure 5. Flow field lines and derivatives for transition section

下载: 全尺寸图片幻灯片

由可见：虽然该线型是基于第1节中第2种思路为指导绘出，但在 $\left|x\right|\geqslant 2$ 时，该线型也满足第1种思路中的二阶连续可导. 因此，将式（10）作为地形边界过渡段曲线是比较合理的.

3. 数值模型及工况设置

本研究采用CFD数值模拟的方法对比分析气流经过平方正弦函数、绕兰金物体不可压缩无旋流动流线、双曲正切函数流线过渡段前后的流场分布、变化过程，计算工况如表1所示.

表 1 计算工况

Table 1. Working conditions for calculation

计算工况	A	B	C	D
过渡段形式	平板参考风场	兰金物体绕流流线	双曲正切函数	平方正弦函数

下载: 导出CSV

| 显示表格

为加快运算速度，设置二维理想台阶模型为地形，以未设置台阶处平板风场为参考风场. 取h=0.6 m，平台长度为9 m，过渡段等效斜率l=1/3，过渡段长度为1.8 m，计算域长度为50 m，高度为18 m，阻塞率设为3.3%以满足规范的要求（小于5%）.

流体经过渡段流经区域，选用结构化四边形网格进行离散. 为准确计算黏性子层的流动，要求壁面首层网格 ${y}^ + \leqslant 1$ . 基于光滑平板湍流边界层理论，估算边界层第1层网格高度为4 × 10⁻⁴ m，边界层内网格高度增长因子取为1.05，边界层以上网格高度增长因子取为1.10. 网格水平方向大小在过渡段处适当加密，各工况计算模型的网格数量见表2.

表 2 网格数量

Table 2. Number of grids 万个

计算工况	A	B	C	D
网格数	244.000	390.230	100.750	281.212

下载: 导出CSV

| 显示表格

设置模型计算域的上边界为对称边界，设置底边为无滑移固定壁面，左边界设置为速度入口，计算域右边界面设置为零压力出口，来流设置为1 m/s的均匀流，方向垂直于入口边界，湍流强度根据一般情况下风洞试验均匀流要求设置为0.5%. 计算区域及边界条件设置如图6所示.

图 6 计算区域及边界条件（单位：m）

Figure 6. Calculation area and boundary conditions (unit: m)

下载: 全尺寸图片幻灯片

数值模拟使用FLUENT软件，为获得更好的近壁面流动解析度，湍流模型使用SST k-ω模型（k为湍流动能密度，ω 为湍流比耗散率），流体设置为不可压缩定常流计算. 采用压力基耦合求解器进行计算，动量方程、湍动能方程及湍流耗散率方程均采用二阶离散格式. 计算收敛精度设置为 $1\times {10}^{-6}$ . 首层网格 $y^+$ 值如所示. 可见：对工况A、B、D，满足采用SST k-ω模型精确计算黏性底层流动时对 $y^+$ 值的要求（<1.0）；对工况C， $y^+$ 值大于1.0的壁面范围较小，仅出现在局部个别单元；所有模型的 $y^+$ 均值均小于0.6，网格划分情况符合计算要求.

表 3 不同工况首层网格

$y^+$ 值

Table 3.

$y^+$ values of the first layer grid under different working conditions

计算工况	$y^+$ 最大值	$y^+$ 平均值
A	0.7615	0.5674
B	0.9946	0.5334
C	1.3259	0.5251
D	0.9808	0.5338

下载: 导出CSV

| 显示表格

为分析流体流经过渡段后不同工况下的风剖面、风攻角等沿程变化情况，即不同过渡段线型的过渡效果，在高度为0.602～1.100的高度区间内选取了26个监控点：在y=0.602～0.610以0.002为间距设置5个监控点；在y=0.610～0.650以0.005为间距设置8个监控点；在y=0.650～1.100以0.010为间距设置13个监控点. 每组监控点沿壁面法向在x=−4.5，−4.0，−3.5，−2.5，−1.5，0处分别设置了6组监控点，分别设置在过渡段的终点到模型中点之间.

4. 不同过渡段曲线流动特性分析

4.1 不同线型过渡段的气流分离特性

是否能够有效抑制来流产生流动分离是衡量地形边界过渡段线型是否适用的重要标准，因为来流流动分离后将可能产生诸如转捩成湍流和局部空化等现象，从而对风场模拟的精度和计算效率产生很大影响.

流体流经几何壁面时存在足够大的逆压梯度是产生来流流动分离的主要原因，此时流体流速将沿流向逐渐降低且有可能降低至负数（式（1）），从而出现回流现象，导致来流产生分离. 据此，可以通过气体与壁面的接触情况判断是否发生分离，当几何壁面表面剪应力 ${\tau }_{{\rm{w}}}\leqslant 0$ 时，代表流体已经产生分离. 因此，确定壁面剪应力的分布情况就能确定来流与壁面的分离情况，图7为基于FLUENT计算结果绘出的壁面剪应力分布情况.

图 7 各工况在过渡段前后水平剪应力分布

Figure 7. Distribution of horizontal shear stress before and after transition section in each working condition

下载: 全尺寸图片幻灯片

由可知：1）在平板风场中， ${\tau }_{{\rm{w}}}$ 沿程基本保持不变. 2） 3类被过渡段线型工况中， ${\tau }_{{\rm{w}}}$ 曲线在距离设置的过渡段起点较远的距离就开始远离参考风场的剪应力曲线；剪应力在到达过渡段前沿程分布基本重合，到达过渡段后剪应力曲线开始产生分化，数值在一段区间内继续减小后迅速增大，超过参考风场的剪应力，达到最大点后又快速降低，回到参考风场工况附近. 3）总体上看，3类过渡段工况的沿程剪应力变化曲线变化趋势相近，而流经过渡段时， ${\tau }_{\rm{w}}$ 波动的陡峭程度不同，出现向上增长、向下回归的先后次序也不同. 4）采用兰金物体绕流流线作为过渡段的工况B和采用平方正弦函数作为过渡段的工况D， ${\tau }_{\rm{w}}$ 沿程未出现负值，即没有出现流动分离，2个工况 ${\tau }_{\rm{w}}$ 延程波动幅度都较小，而工况D的 ${\tau }_{\rm{w}}$ 波动程度相较工况B更小一些. 采用双曲正切函数作为过渡段的工况C， ${\tau }_{\rm{w}}$ 在过渡段开始的一段距离出现负值，发生了流动分离，且 ${\tau }_{\rm{w}}$ 沿程波动剧烈，波动幅度较大.

过渡段线型的不同是造成上述差异的原因：1）对于双曲正切函数，其在连接边界严格满足式（2），过渡段曲线在开头处斜率较小，而在中间段距离内陡然增加，斜率较大；过渡段结尾接头处斜率又较小，曲线平缓；相应的双曲正切函数工况B中的 ${\tau }_{\rm{w}}$ 在初始过渡段平缓增长，而后流体在较短距离内较快抬升，对应部分的 ${\tau }_{\rm{w}}$ 快速增大，后快速回落，而后又经过较长较平缓的阶段. 2）对于平方正弦函数，其连接边界条件仅满足式（2）中的前2项，但其斜率变化平缓、分布均匀；相应的工况D中平方正弦的 ${\tau }_{\rm{w}}$ 在过渡段各阶段波动幅度较小；对于兰金物体绕流流线，其在连接边界基本满足式（2），各段斜率变化较双曲正切函数均匀一些，但不如平方正弦函数；相应的工况B中兰金物体绕流流线的 ${\tau }_{\rm{w}}$ 在过渡段波动程度较工况C小，波动幅度和工况D基本一致.

4.2 气流经过过渡段后平均风速的沿程分布

本节主要考察不同过渡段线型中沿程平均风特性变化情况，并与平板参考风场对照后，进一步评判各类过渡段曲线的优劣.

参考风场工况在典型监测位置处的风剖面如图8所示. 图中，v为风速. 由图可见：沿程不同位置处的风剖面已基本重合，这说明平板边界层风场已经充分发展、基本稳定，可作为参考风场. 不同过渡段曲线下，典型监测位置处的风剖面如图9所示.

图 8 平板参考风场沿程不同位置处风剖面

Figure 8. Wind profiles at different positions of plane reference wind field

下载: 全尺寸图片幻灯片

图 9 不同位置处各工况平均风剖面对比

Figure 9. Comparison of average wind profiles at different locations under each working condition

下载: 全尺寸图片幻灯片

由图9可见：1）随着流动沿程发展，各过渡段工况的风剖面和参考风场风剖面的差距逐渐缩小，但直到模型中心位置（x=0）处，各过渡段工况与平板参考风场仍有一定的差距. 2）在不同位置处的风剖面均有平方正弦函数（工况D）与参考风场的差距最小，工况B、C与参考风场的差距较大，且两工况基本重合. 3）工况B、C、D的梯度风速基本一致，风剖面在距离壁面较远处趋于重合. 4）在距离壁面较近处工况D的风剖面和参考风场重合，在壁面附近风速较小，随着流体远离过渡段，重合的部分增多；工况B、C的风剖面在壁面附近处风速较大，随着远离过渡段，风速有一定的减小，显然和参考风场差距较大.

采用平方正弦函数作为过渡段曲线时，在平均风剖面的保持方面性能较采用兰金物体绕流流线和双曲正切函数效果更优. 其原因在于：平方正弦函数的导函数即斜率沿程分布较均匀，曲线变化较平缓，对来流的扰动较小，没有发生流动分离，实现了气流的平稳过渡.

4.3 气流经过过渡段后风攻角的沿程分布

各工况在典型监测位置处的风攻角 $\alpha$ 剖面如图10所示.

图 10 不同位置处各工况风攻角剖面对比

Figure 10. Comparison of wind attack angle profiles at different locations under each working condition

下载: 全尺寸图片幻灯片

由图10可知：1）参考风场的风攻角基本为0；2）随着来流沿过渡段流动，3类过渡段工况的风攻角逐渐减小；3）随着距离壁面高度的增加，风攻角逐渐增大，增大到一定程度后，保持不变；4）在接近过渡段终点的位置处，工况C的风攻角剖面较工况B、D更为接近参考风场，这显然是由于工况C的过渡段在接头处平缓段较长的原因；5）对于工况D，在距离壁面较近的一段，风攻角始终保持为0，这说明来流经过平方正弦曲线过渡后，边界层内始终保持为层流.

4.4 气流经过过渡段后湍流动能的沿程分布

本节通过考察不同过渡段线型中沿程湍流动能 $\kappa$ 的变化规律（图11所示），以来流流经过渡段后产生流动分离时流体的稳定程度为标准进一步评判各类过渡段曲线的优劣.

图 11 不同位置处各工况湍流动能对比

Figure 11. Comparison of turbulent kinetic energy at different locations under each working condition

下载: 全尺寸图片幻灯片

由图11可见：1）随着距离壁面高度的增加，各工况湍流动能基本重合且趋近于0. 2）平板参考风场（工况A）距离壁面相同高度处湍流动能沿程变化较小，比较平稳，沿高度变化为先增后减，最终趋于平稳. 3）在贴近壁面的区域，除工况B在流经过渡段后湍流动能在平台初段有所增加外，各过渡段工况沿程湍动能度逐渐减少，来流流动分离后转捩成湍流的部分逐渐降低. 这是来流流经过渡段后接触平直壁面时沿程流速逐渐下降，同距离壁面较高处风速差有所提高而引起的. 4）对于工况B、C，湍流动能沿高度变化趋势为逐渐降低，与平板参考风场变化趋势不同. 5）对于工况D，湍流动能沿程变化较大，但来流流经过渡段后变化趋势与参考风场工况比较贴合.

5. 结　论

1）以桥位为中心截取部分地形进行风场流动模拟，对于造成的台阶问题，可采取添加地形边界过渡段的形式加以解决.

2）构造地形模型边界过渡段线型宜遵循以下2个原则：其一是过渡段曲线本身满足二阶导函数连续，在连接边界处满足一阶导函数连续；其二是过渡段曲线一阶导函数即斜率应沿程分布较均匀，避免出现较大的窄带尖峰.

3）在可以作为过渡段曲线线型的兰金物体绕流流线、双曲正切函数和平方正弦函数中，平方正弦函数曲线绕流性能最优，可作为合理的过渡段曲线线型被选用进行桥址区风场模拟.

本文尚存在一些不足之处，研究工作均基于数值模拟开展，后续将通过风洞试验加以验证.

图 1 三相电缆的双极式直流运行方案

Figure 1. Bipolar DC operation scheme of three-phase cables

下载: 全尺寸图片幻灯片

图 2 三相电缆的单极式直流运行方案

Figure 2. Monopolar DC operation scheme of three-phase cables

下载: 全尺寸图片幻灯片

图 3 三相电缆的TWBS-HVDC直流运行方案

Figure 3. TWBS-HVDC operation scheme of three-phase cables

下载: 全尺寸图片幻灯片

图 4 电缆的敷设环境

Figure 4. Laying environment of cable

下载: 全尺寸图片幻灯片

图 5 网格剖分

Figure 5. Mesh settings

下载: 全尺寸图片幻灯片

图 6 不同直流拓扑和直流负荷的缆芯温度

Figure 6. Cable core temperature under different DC topologies and loads

下载: 全尺寸图片幻灯片

图 7 通过210.0 A电流时单极式电缆的温度场分布

Figure 7. Temperature field distribution of monopolar cable at current of 210.0 A

下载: 全尺寸图片幻灯片

图 8 通过240.0 A电流时双极式电缆的温度场分布

Figure 8. Temperature field distribution of bipolar cable at current of 240.0 A

下载: 全尺寸图片幻灯片

图 9 不同直流拓扑和工作温度的电缆载流量

Figure 9. Cable ampacity under different DC topologies and operating temperatures

下载: 全尺寸图片幻灯片

图 10 热水管与电缆的相对位置

Figure 10. Distance from hot water pipe to cable

下载: 全尺寸图片幻灯片

图 11 电缆载流量随热水管道和电缆距离的变化

Figure 11. Variation of cable ampacity with distance from hot water pipe to cable

下载: 全尺寸图片幻灯片

图 12 集群电缆加载相同直流负荷时的温度场分布

Figure 12. Temperature field distribution of cluster cables under the same DC load

下载: 全尺寸图片幻灯片

图 13 单根电缆运行时的温度场分布

Figure 13. Temperature field distribution of a single cable in operation

下载: 全尺寸图片幻灯片

图 14 集群电缆加载不同负荷时的温度场分布

Figure 14. Temperature field distribution of cluster cables under different loads

下载: 全尺寸图片幻灯片

图 15 电缆温度场分布

Figure 15. Temperature field distribution of cable

下载: 全尺寸图片幻灯片

图 16 电流不平衡度和三芯电流和之间的关系

Figure 16. Relationship between current unbalance and three-core current sum

下载: 全尺寸图片幻灯片

表 1 电缆材料与尺寸参数

Table 1. Cable materials and size parameters

结构名称	材料	厚度/mm	直径/mm
导体	铜		13.0
导体屏蔽层	半导体混合物	0.8	14.6
XLPE 屏蔽层	XLPE	10.5	35.6
绝缘屏蔽层	半导体混合物	1.0	37.6
铜屏蔽层	铜	0.2	38.0
包带	聚丙烯	0.4	83.2
内护套	聚氯乙烯	2.0	87.2
铠装层	钢	1.6	90.4
外护套	聚氯乙烯	5.0	100.4

下载: 导出CSV

表 2 电缆材料的物理参数

Table 2. Physical parameters of cable materials

材料	密度/ （kg·m⁻³）	比热容/ （J·（kg·K）⁻¹）	热导率/ （W·（m·K）⁻¹）
铜	8900	380	385.00
XLPE	1200	1000	0.29
聚氯乙烯	1380	1000	0.16
半导电材料	1200	1100	0.28
填充材料	550	1900	0.25
钢	7850	450	45.00

下载: 导出CSV

表 3 热水管道的材料和结构参数

Table 3. Hot water pipe materials and structural parameters

结构名称	材料	厚度/mm	半径/mm
工作区	水		516
保温层	硬质聚氨酯	65	581
外护管	高密度聚乙烯	15	596
螺旋焊缝钢管	钢	14	610

下载: 导出CSV

参考文献(21)

[1]	郑欢,江道灼,杜翼. 交流配电网与直流配电网的经济性比较[J]. 电网技术,2013,37(12): 3368-3374. ZHENG Huan, JIANG Daozhuo, DU Yi. Economic comparison of AC and DC distribution system[J]. Power System Technology, 2013, 37(12): 3368-3374.
[2]	ALASSI A, BAÑALES S, ELLABBAN O, et al. HVDC transmission: technology review, market trends and future outlook[J]. Renewable and Sustainable Energy Reviews, 2019, 112: 530-554. doi: 10.1016/j.rser.2019.04.062
[3]	WANG M, AN T, ERGUN H K, et al. Review and outlook of HVDC grids as backbone of transmission system[J]. CSEE Journal of Power and Energy Systems, 2021, 7(4): 797-810.
[4]	徐政,许烽. 输电线路交改直的关键技术研究[J]. 高电压技术,2016,42(1): 1-10. XU Zheng, XU Feng. Research on key technologies of AC-to-DC transmission lines conversion[J]. High Voltage Engineering, 2016, 42(1): 1-10.
[5]	YU J, SMITH K, URIZARBARRENA M, et al. Initial designs for the ANGLE DC project; converting existing AC cable and overhead line into DC operation[C]//13th IET International Conference on AC and DC Power Transmission (ACDC 2017). Manchester: IET, 2017: 1-6.
[6]	刘英,曹晓珑,何子兰,等. 现役交流XLPE电缆配电线路改为直流运行的技术方案及实例分析[J]. 中国电机工程学报,2016,36(1): 96-103. LIU Ying, CAO Xiaolong, HE Zilan, et al. Technical scheme and case study of the uprating renovation of existing XLPE cables from AC distribution system to DC operation[J]. Proceedings of the CSEE, 2016, 36(1): 96-103.
[7]	LIU Y, CAO X L, FU M L. The upgrading renovation of an existing XLPE cable circuit by conversion of AC line to DC operation[J]. IEEE Transactions on Power Delivery, 2017, 32(3): 1321-1328. doi: 10.1109/TPWRD.2015.2496178
[8]	郑雁翎. 电缆线路的电缆群载流量优化数值计算模型研究[J]. 高电压技术,2015,41(11): 3760-3765. ZHENG Yanling. Research on optimal numerical caculation model of cable groups ampacity in cable route[J]. High Voltage Engineering, 2015, 41(11): 3760-3765.
[9]	郑雁翎,许志亮,张冠军,等. 采用MATLAB仿真的变电站高压进线温度场和载流量数值计算[J]. 高电压技术,2012,38(3): 566-572. ZHENG Yanling, XU Zhiliang, ZHANG Guanjun, et al. Numerical calculation of temperature field and ampacity of high-voltage lines using MATLAB[J]. High Voltage Engineering, 2012, 38(3): 566-572.
[10]	王启隆,王国海,陈向荣,等. 10 kV交流XLPE电缆改为直流运行的热电耦合仿真[J]. 西南交通大学学报,2022,57(1): 46-54. doi: 10.3969/j.issn.0258-2724.20200111 WANG Qilong, WANG Guohai, CHEN Xiangrong, et al. Thermo-electric coupling simulation for 10 kV AC XLPE cable in DC operation[J]. Journal of Southwest Jiaotong University, 2022, 57(1): 46-54. doi: 10.3969/j.issn.0258-2724.20200111
[11]	刘士利,李宁,蔡国伟,等. 66kV交流交联聚乙烯电缆线路改为直流运行的直流载流量[J]. 高电压技术,2017,43(5): 1664-1669. LIU Shili, LI Ning, CAI Guowei, et al. DC ampacity of 66 kV AC XLPE cable line transformed into DC operation[J]. High Voltage Engineering, 2017, 43(5): 1664-1669.
[12]	LI J X, YANG F, FANG Z K, et al. Assessment on capacity expansion of AC XLPE cable under DC operation with the constraint of temperature and electric field intensity[J]. IET Generation, Transmission & Distribution, 2022, 16(13): 2612-2622.
[13]	许烽,徐政. 一种适用于交流线路改造成直流的扩展式双极直流输电结构[J]. 中国电机工程学报,2014,34(33): 5827-5835. XU Feng, XU Zheng. An extended bipole HVDC structure with three wires for conversion of AC lines to HVDC[J]. Proceedings of the CSEE, 2014, 34(33): 5827-5835.
[14]	LIU Y, ZHANG S D, CAO X L, et al. Simulation of electric field distribution in the XLPE insulation of a 320 kV DC cable under steady and time-varying states[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2018, 25(3): 954-964. doi: 10.1109/TDEI.2018.006973
[15]	陈铮铮,赵健康,欧阳本红,等. 直流电缆料工作温度和击穿特性的纳米改性研究[J]. 高电压技术,2015,41(4): 1214-1227. CHEN Zhengzheng, ZHAO Jiankang, OUYANG Benhong, et al. Study on the nanoparticle-modification of working temperature and breakdown characteristics for insulating materials in DC cables[J]. High Voltage Engineering, 2015, 41(4): 1214-1227.
[16]	杨黎明,朱智恩,杨荣凯,等. 柔性直流电缆绝缘料及电缆结构设计[J]. 电力系统自动化,2013,37(15): 117-124. doi: 10.7500/AEPS20130607005 YANG Liming, ZHU Zhien, YANG Rongkai, et al. Insulation material and structure design of HVDC flexible cables[J]. Automation of Electric Power Systems, 2013, 37(15): 117-124. doi: 10.7500/AEPS20130607005
[17]	Internation Electrotechnical Commission. Electric cables—calculation of the current rating—part 1-3: current rating equations (100% load factor) and calculation of losses—current sharing between parallel single-core cables and calculation of circulating current losses: IEC 60287-1-3[S]. Geneva: IEC National Committees, 2014.
[18]	中华人民共和国住房和城乡建设部. 城镇供热直埋热水管道技术规程:CJJ/T 81—2013 [S]. 北京:中国建筑工业出版社,2013.
[19]	陈向荣,王启隆,于竞哲,等. 10 kV交流XLPE电缆在不同直流拓扑结构和敷设方式下的直流载流量仿真研究[J]. 高电压技术,2021,47(11): 4044-4054. CHEN Xiangrong, WANG Qilong, YU Jingzhe, et al. Simulation research on DC ampacity of 10 kV AC XLPE cable under different DC operation topologies and laying modes[J]. High Voltage Engineering, 2021, 47(11): 4044-4054.
[20]	曾嵘,赵宇明,赵彪,等. 直流配用电关键技术研究与应用展望[J]. 中国电机工程学报,2018,38(23): 6790-6801. ZENG Rong, ZHAO Yuming, ZHAO Biao, et al. A prospective look on research and application of DC power distribution technology[J]. Proceedings of the CSEE, 2018, 38(23): 6790-6801.
[21]	阮羚,赵艾萱,邓丹,等. 三芯电缆不平衡电流对温度分布的影响[J]. 高电压技术,2018,44(8): 2704-2709. RUAN Ling, ZHAO Aixuan, DENG Dan, et al. Influence of unbalanced current in 3-core power cable on temperature distribution[J]. High Voltage Engineering, 2018, 44(8): 2704-2709.