可再生能源

截面积	传输容量/MW	价格/万元	挂载风机数
70	24	77.8	2
120	33	91.16	3
185	41	106.76	4
240	47	117.23	5
400	60	155.29	6
800	150	550
1000	200	734
1 200	300	887
1400	350	1200
1600	400	1500

截面积	传输容量/MW	价格/万元	挂载风机数
70	24	77.8	2
120	33	91.16	3
185	41	106.76	4
240	47	117.23	5
400	60	155.29	6
800	150	550
1000	200	734
1 200	300	887
1400	350	1200
1600	400	1500

	截面积mm	拓扑 1	拓扑 2
优化前	优化后	节约	优化前	优化后	节约
不同截面积对应线路长度/km	70	70	66.2	3.8	64.6	64.6	0
120	65.8	25	40.8	32.3	32.3	0
185	31.2	25	6.2	32.8	64	-31.2
240	77	20	57	36.5	61.9	-25.4
400	0	71.8	-71.8	104.5	63.6	40.9
总长度 $/\mathrm{{km}}$	244	208	36	270.7	286.4	-15.7

	截面积mm	拓扑 1	拓扑 2
优化前	优化后	节约	优化前	优化后	节约
不同截面积对应线路长度/km	70	70	66.2	3.8	64.6	64.6	0
120	65.8	25	40.8	32.3	32.3	0
185	31.2	25	6.2	32.8	64	-31.2
240	77	20	57	36.5	61.9	-25.4
400	0	71.8	-71.8	104.5	63.6	40.9
总长度 $/\mathrm{{km}}$	244	208	36	270.7	286.4	-15.7

	截面积mm	拓扑 1	拓扑 2
优化前	优化后	节约	优化前	优化后	节约
不同截	70	5 446.00	5 150.36	295.64	5025.88	5 025.88	0.00
面积对	120	5 998.33	2 279.00	3 719.33	2944.47	2 944.47	0.00
应线路	185	3 330.91	2 669.00	661.91	3501.73	6 832.64	-3330.91
成本/	240	9026.71	2 344.60	6 682.11	4 278.90	7 256.54	-2977.64
万元	400	0.00	11 149.82	11 149.82	16 227.81	9 876.44	6 351.36
成本/亿元	2.38	2.36	0.0209	3.20	3.19	0.004 3

	截面积mm	拓扑 1	拓扑 2
优化前	优化后	节约	优化前	优化后	节约
不同截	70	5 446.00	5 150.36	295.64	5025.88	5 025.88	0.00
面积对	120	5 998.33	2 279.00	3 719.33	2944.47	2 944.47	0.00
应线路	185	3 330.91	2 669.00	661.91	3501.73	6 832.64	-3330.91
成本/	240	9026.71	2 344.60	6 682.11	4 278.90	7 256.54	-2977.64
万元	400	0.00	11 149.82	11 149.82	16 227.81	9 876.44	6 351.36
成本/亿元	2.38	2.36	0.0209	3.20	3.19	0.004 3

编号			装机容量/MW
1	-53.09	51.11	10
2	-48.16	51.39	20
3	-33.27	-24.20	30
4	-17.11	-22.55	20
5	-13.16	51.29	30
6	-19.80	-45.94	10
7	100.49	-31.29	20

编号			装机容量/MW
1	-53.09	51.11	10
2	-48.16	51.39	20
3	-33.27	-24.20	30
4	-17.11	-22.55	20
5	-13.16	51.29	30
6	-19.80	-45.94	10
7	100.49	-31.29	20

截面积mm	优化前	优化后	节约
长度/km	成本/万元	长度/km	成本/万元	长度/km	成本/万元
240	123.70	14501.35	148.90	17 455.55	-25.20	-2954.20
400	203.90	31 663.63	50.90	7904.26	153.00	23 759.37
800	94.10	51 755.00	97.20	53 460.00	-3.10	-1705.00
1600	141.40	212 100.00	119.30	178 950.00	22.10	33 150.00
总计	563.10	310 019.98	416.30	257 769.81	146.80	52 250.17
收益/亿元		10.35		15.50		-5.16

截面积mm	优化前	优化后	节约
长度/km	成本/万元	长度/km	成本/万元	长度/km	成本/万元
240	123.70	14501.35	148.90	17 455.55	-25.20	-2954.20
400	203.90	31 663.63	50.90	7904.26	153.00	23 759.37
800	94.10	51 755.00	97.20	53 460.00	-3.10	-1705.00
1600	141.40	212 100.00	119.30	178 950.00	22.10	33 150.00
总计	563.10	310 019.98	416.30	257 769.81	146.80	52 250.17
收益/亿元		10.35		15.50		-5.16

考虑相关性与平滑效应的海上风电场集群汇集系统拓扑优化

PDF下载

谢珍建 ¹ , 许偲轩 ¹ , 李辰 ¹ , 李海波 ² , 江坷滕 ² , 鲁宗相 ²

可再生能源 | 2024,42(1): 79-86

收起

可再生能源 | 2024, 42(1): 79-86

考虑相关性与平滑效应的海上风电场集群汇集系统拓扑优化

全屏

谢珍建¹, 许偲轩¹, 李辰¹, 李海波², 江坷滕², 鲁宗相²

作者信息

¹ 国网江苏省电力有限公司经济技术研究院江苏南京 610095

² 清华四川能源互联网研究院四川成都 610000

谢珍建(1980-),男,博士研究生,研究员级高级工程师,主要从事新型电力系统规划与决策仿真分析相关方面研究。E-mail: 13584095327@163.com。

Topology optimization of offshore wind farm cluster collection system considering correlation and smoothing effects

Zhenjian Xie¹, Sixuan Xu¹, Chen Li¹, Haibo Li², Keteng Jiang², Zongxiang Lu²

Affiliations

¹ Institute of Economics and Technology State Grid Jiangsu Electric Power Co., Ltd. Nanjing 610095 China

² Tsinghua Sichuan Energy Internet Research Institute Chengdu 610000 China

出版时间: 2024-01-20

文章导航

摘要

收起

海上风力发电在时间尺度上具有间歇性和不确定性,空间尺度上具有互补性。传统规划方法线路容量利用不足,风力发电成本高。文章结合大规模风力发电出力特性,考虑风电平滑效应、输电工程建设成本、风力弃电成本、输送风力发电收益,建立了海上风电场汇集系统的经济性评估模型,采用改进的遗传算法对模型进行求解。基于动态变权最小生成树对遗传算法进行改进。最后,以江苏某海上风电场集群为例,对风机拓扑结构和风电场汇集连接方式进行了优化求解。结果表明,所提算法具有较好的寻优性和收敛性,所建立的海上风电场及内部风机拓扑优化模型能有效减少工程投资。

关键词

平滑效应 / 海上风电场汇集系统 / 遗传算法 / 拓扑优化

Abstract

收起

Due to the intermittency and uncertainty of wind power on the time scale and complementarity on the spatial scale, the line capacity is underutilized and the cost of wind power is high in traditional planning. Combined with the output characteristics of largescale wind power, considering the smoothing effect, the construction cost of transmission projects, the cost of wind power curtailment, and the income, the economic evaluation model of offshore wind farm aggregation system is established, and the model is solved by using an improved genetic algorithm. The genetic algorithm is improved based on the minimum spanning tree of dynamic weight variation, and the coding, selection, crossing, and mutation links in the genetic algorithm are improved to promote the efficiency of the algorithm. Taking a wind farm cluster in Jiangsu Province as an example, the topology of wind turbines and the convergence and connection mode of wind farms are optimized, and the results show that the proposed algorithm has good optimization and convergence. And the established topology optimization model of wind farm and internal wind turbines can effectively reduce the engineering investment.

Key words

smoothing effect / collection system of offshore wind generation / genetic algorithm / topology optimization

引用本文

谢珍建, 许偲轩, 李辰, 李海波, 江坷滕, 鲁宗相. 考虑相关性与平滑效应的海上风电场集群汇集系统拓扑优化. 可再生能源, 2024 , 42 (1) : 79 -86 .

Zhenjian Xie, Sixuan Xu, Chen Li, Haibo Li, Keteng Jiang, Zongxiang Lu. Topology optimization of offshore wind farm cluster collection system considering correlation and smoothing effects[J]. Renewable Energy Resources, 2024 , 42 (1) : 79 -86 .

正文

收起

0 引言

收起

与内陆风力资源相比, 海上风力资源更加稳定, 风能质量更好, 发电效率更高, 对人类生活的直接影响较小 ^{[ 1 - 3 ]} 。风电汇集系统承担着汇集风机电能并向内陆电网输送的任务。在海上风电汇集系统规划过程中需保证风电发电效率和可靠性需求, 因此海上风电场从环境和技术的角度相较于其他发电方式具有更昂贵的建造、运维费用, 其电网接入成本也是重要的经济成本之一。

海上风电场集电系统拓扑优化是对风机机组、汇集站和集群中心站的空间位置和连接方式进行优化 ^{[ 4 , 5 ]} 。一些学者首先从不改变现有拓扑的角度出发,针对“弃风”问题设计新设备来提高系统可靠性。文献[ 6 ]提出了直流集电器的概念, 设计了一个集汇集升压、控制保护、故障隔离的综合性平台取代传统的分布式升压站。文献[ 7 ]设计了一种分流电路来保证风机在最大功率追踪条件下的输出电压稳定, 降低了风机间的耦合性和脱机率。另一种优化方向为集电系统的拓扑结构多层设计优化。已有的优化方案有变风机选址和变集电平台选址的双层优化 ^{[ 8 - 11 ]} 和子区域风电汇集站选址和中心汇集站选址的双层优化 ^{[ 4 ]} 等。针对拓扑结构的优化, 文献[ 8 - 14 ]提出了改进遗传算法、模糊化算法、混合生态共生算法、拓扑连线优化算法和基于改进最小生成树的经典 Prim 算法。这些算法均以实现综合经济效益为目标, 给出了最终的拓扑结构。文献[ 8 - 10 ]将风机划分为集群子区域,多个区域之间采用星型结构连接。文献[ 11 - 14 ]尽管给出了可靠性和可行性分析, 但由于没有考虑集群海上风电场出力之间的互补性与相关性, 其规划结果有进一步优化的空间。

现有的研究更侧重从拓扑结构的经济性、可靠性等角度开展优化研究, 但是并未考虑海上风电出力特性与平滑效应 ^{[ 15 , 16 ]} 。根据实际风电场数据可知,海上风电场出力利用小时数较高,且具有较强的相关性和空间平滑效应。相关性大导致风电场集群出力仍有较强的波动性, 空间平滑效应则削弱了波动性, 使集群出力表现出更加稳定的特征。实际工程中, 若不考虑海上风电场集群出力特性对汇集和输送系统容量的影响, 则会导致规划容量过大, 经济性较差。因此在汇集系统优化中,应结合风电场出力时序数据,通过容量优化合理平衡弃风电量与经济性。

本文提出了一种考虑相关性与平滑效应的海上风电场集群汇集系统拓扑优化方法。首先基于我国某沿海省份风电场实际历史出力数据, 开展相关性与平滑效应分析, 进一步考虑系统投资成本、运行成本与弃风成本等,以净收益最大为目标,建立了风电场集群汇集系统优化规划模型。采用改进遗传算法对风电场内部风机拓扑结构和风电场汇集拓扑进行了优化分析, 并以沿海某风电场集群作为对象, 验证了本文所提出方法的有效性与正确性。

1 海上风电汇集系统

收起

本文简要介绍海上风电汇集系统, 并结合风电场实际历史出力数据分析其相关性与空间平滑效应。为简化分析, 本文做出如下合理假设: ①海上风电场的装机容量及其年出力特性已知; ②陆上电网足以输送并消纳通过海上风电场群汇集系统注入的功率；③暂不考虑多年滚动规划以及电源、电网规划的时序不匹配问题。

1.1 集电系统拓扑优化模型概述

海上风电汇集系统首先通过输电网将风机发出的电能通过海上升压站统一汇集, 再通过线缆传输到陆上电网。江苏某海上风电场群采用树形拓扑,模型如图 1 所示。

由图 1 可知,风电从各个风电场出发,经过多个风电场汇集后, 通过高压输电网络传输至海上升压站。在风电场和海上升压站选址确定的条件下, 以风电场全生命周期效益最大为目标, 优化风电场分组汇集结构, 汇集过程中线路和变电站的技术经济参数是优化结果的重要影响因素。

1.2 海上风电场集群特性分析

在 “风电机组-风电场-海上升压站-输电系统一陆上电网”的汇集传递过程中,海上风电的有功输出在空间上具有互补性, 通过汇集不同区域的海上风电有利于减少出力波动性, 使得风电机群出力特性曲线更加平滑, 也称作风电出力的平滑效应或汇聚效应。对江苏某地区的数个海上风电场组成的集群进行历史数据统计, 得到风电场集群出力曲线,如图 2 所示。

由图 2 可以看出, 多个风电场出力具有明显的相关性, 同时由于空间平滑效应具有一定的互补性, 相比单个风电场波动性有所减少, 最大出力标幺值降低约为 0.85 p.u.,验证了风电场集群的平滑效应。在风电场汇集系统容量优化中,应充分利用风电场集群出力这一特性, 对汇集线路的容量进行优化, 在满足风电消纳目标的前提下最大化系统经济效益。

本文提出了一种考虑平滑效应海上风电场汇集系统拓扑优化方法, 根据海上风电场集群出力特性, 建立系统接入协调的规划模型, 以经济性最优为目标, 依据未送出及送出风电电量的损失及效益、建设成本等因素进行优化分析, 对风电场内部的风机拓扑结构进行优化, 并对风电场汇集拓扑进行优化, 从而完成整个海上风电场汇集系统的拓扑优化。

2 海上风电汇集系统优化建模

收起

2.1 目标函数

对于风电场或风电场集群的送出线路, 如按照装机容量进行配置, 则会经常导致输电线路轻载,资产利用率低下；如果送出线路容量过低,又会造成大量弃风。因此需要建立经济性综合优化模型, 协调输电线路的建设容量与弃风损失。输电线路综合效益的成本及收益为

(1)${I}_{\mathrm{{NPV}}}= \mathop{\sum }\limits_{{t = 0}}^{N}{\left({B}_{\text{total }}- {C}_{\mathrm{L}}\right)}_{t}{\left( 1 + r\right)}^{-t}$

式中: ${B}_{\text{total }}$ 为年度收益; ${C}_{\mathrm{L}}$ 为年度成本; $r$ 为利率; $N$ 为运行年限。

风电场送出电能为线路的收益来源, 因此传输容量限制造成的弃风,即未送出电量,也应作为罚款计入线路的收益。忽略风电场送出电能不同时段的价格差异,风电场年度收益 ${B}_{\text{total }}$ 由风电收益 ${B}^{\mathrm{{TE}}}$ 和弃电补偿成本 ${C}^{\text{cur }}$ 组成,其中, ${B}^{\mathrm{{TE}}}$ 主要由风电站年发电量决定, ${C}^{\mathrm{{cur}}}$ 主要由输电容量限制导致的风电年弃电量决定,计算方法为

(2)${C}^{\mathrm{{cur}}}= \left\lbrack {{\left( 1 + r\right)}^{N}- 1}\right\rbrack /\left\lbrack {r{\left( 1 + r\right)}^{N}}\right\rbrack {p}_{\mathrm{b}}{E}^{\mathrm{{cur}}}$

(3)${B}^{\mathrm{{TE}}}= \left\lbrack {{\left( 1 + r\right)}^{N}- 1}\right\rbrack /\left\lbrack {r{\left( 1 + r\right)}^{N}}\right\rbrack {p}_{\mathrm{o}}E $

式中: ${B}^{\mathrm{{TE}}}$ 为风电电量收益; ${C}^{\text{cur }}$ 为弃电补偿成本; ${p}_{\mathrm{o}}$ , ${p}_{\mathrm{b}}$ 为价格参数； ${E}^{\text{cur }}$ 为风电场的年弃电量； $E$ 为风电场的年发电量。

图 3 为某海上风电场的年持续出力曲线,即 ${P}^{\mathrm{{dur}}}\left( t\right)$ 。

当风电场理论出力大于其出口线路有功功率限制 ${P}^{\mathrm{{LL}}}$ 时,超过 ${P}^{\mathrm{{LL}}}$ 的部分将被弃电。根据上述定义,风电场的年弃电量 ${E}^{\text{cur }}$ 和年发电量 $E$ 分别为

(4)${E}^{\text{cur }}= \left\{\begin{array}{l} 0,{P}^{\mathrm{{IL}}}\geq {P}^{\mathrm{{dur}},\max }\\{\int }_{0}^{{T}^{\mathrm{{IL}}}}{P}^{\mathrm{{dur}}}\left( t\right)\mathrm{d}t -{P}^{\mathrm{{IL}}}{T}^{\mathrm{{IL}}},0 \leq {P}^{\mathrm{{IL}}}< {P}^{\mathrm{{dur}},\max }\end{array}\right.$

(5)$ E =\left\{\begin{array}{l}{\int }_{0}^{8760}{P}^{\mathrm{{dur}}}\left( t\right)\mathrm{d}t,{P}^{\mathrm{{LL}}}\geq {P}^{\mathrm{{dur}},\max }\\{P}^{\mathrm{{LL}}}{T}^{\mathrm{{LL}}}+ {\int }_{{T}^{\mathrm{{LL}}}}^{8760}{P}^{\mathrm{{dur}}}\left( t\right)\mathrm{d}t,0 \leq {P}^{\mathrm{{LL}}}< {P}^{\mathrm{{dur}},\max }\end{array}\right.$

式中: ${T}^{\mathrm{{LL}}}$ 为出力曲线大于 ${P}^{\mathrm{{LL}}}$ 的时间; ${P}^{\mathrm{{dur}},\max }$ 为风电场的年最大出力; ${E}^{0}$ 为风电场年理论可发电量。

根据持续出力曲线,可得:

(6)$ E ={E}^{0}- {E}^{\text{cur }}= {\int }_{0}^{8760}{P}^{\text{dur }}\left( t\right)\mathrm{d}t -{E}^{\text{cur }}$

2.2 约束条件

在对海上风电汇集系统进行优化建模时, 应分别对海底线缆的最大负载电流、线缆的最大传输功率、线缆能够挂载的最大风机数设置约束条件,并且考虑海底线缆不允许交叉铺设的问题。

①海缆传输容量约束

海缆所传输的电能不超过海缆的最大传输容量,即:

(7)$\mathop{\sum }\limits_{{j = 1}}^{{n}_{i}}{S}_{ij}\leq {S}_{i\max }\left({i = 1,2\cdots , N}\right)$

式中: ${S}_{ij}$ 为与海缆 $i$ 相连的风机 $j$ 的电能; ${n}_{i}$ 为与海缆 $i$ 相连的风机数量； ${S}_{i\max }$ 为海缆 $i$ 的最大传输容量; $N$ 为海缆数量。

②海缆载流量约束

海缆载流量约束包括单风机海缆和汇集海缆两部分, 其表达式为

(8)$\left\{{\begin{array}{l}{I}_{ij}\geq {P}_{\mathrm{w}j}/\left({\sqrt{3}{U}_{ij}\cos {\varphi }_{j}}\right)\\{I}_{i}\geq \mathop{\sum }\limits_{{j = 1}}^{{n}_{i}}{P}_{\mathrm{w}{ij}}/\left({\sqrt{3}{U}_{i}\cos {\varphi }_{j}}\right)\end{array}\left({i = 1,2\cdots , N}\right)}\right.$

式中: ${I}_{ij}$ 为连接风机 $j$ 的电缆 $i$ 的载流量; ${P}_{\mathrm{w}j}$ 为风机 $j$ 的额定出力; ${U}_{ij}$ 为海缆 ${ij}$ 的额定电压; $\cos {\varphi }_{j}$ 为风机 $j$ 的功率因数; ${I}_{i}$ 为连接了 ${n}_{i}$ 台风机的电缆 $i$ 的载流量; ${U}_{i}$ 为汇集海缆 $i$ 的额定电压。

③海缆承受风机数量约束

海缆载流量有限,结合式 (8) 可知,每条海缆能汇集的风机数量 ${n}_{i}$ 有限,其表达式为

(9)${n}_{i}\leq \sqrt{3}{U}_{i}{I}_{\max i}\cos {\varphi }_{i}/{P}_{\mathrm{w}}$

式中: ${I}_{\max i}$ 为汇集电缆 $i$ 的最大载流量; $\cos {\varphi }_{i}$ 为汇集海缆 $i$ 的功率因数; ${P}_{\mathrm{w}}$ 为与汇集海缆 $i$ 相连的风机的平均额定出力。

④海缆交叉规避约束

在实际工程中, 海缆之间不允许交叉。将风机位置的经纬度转化为坐标, 根据风机坐标信息进行海缆交叉判断, 其表达式为

(10)$\left({{Q}_{1}- {P}_{1}}\right)\otimes \left({{P}_{2}- {P}_{1}}\right)\cdot \left({{Q}_{2}- {P}_{1}}\right)\otimes \left({{P}_{2}- {P}_{1}}\right)> 0 $

式中: ${P}_{1},{P}_{2}$ 和 ${Q}_{1},{Q}_{2}$ 为 4 个风机的坐标; 和 $\otimes$ 分别表示点积和叉积计算。

⑤汇集变电站容量约束

海上风电场群的并网运行需要先汇集到海上升压站再进行输送, 海上升压站的容量不小于各汇集电缆容量, 其表达式为

(11)$\mathop{\sum }\limits_{{i = 1}}^{{N}_{i}}{S}_{i}\leq {S}_{\text{asc }}$

式中: ${S}_{i}$ 为电缆容量; ${N}_{i}$ 为汇集到升压站的电缆数量; ${S}_{\text{asc }}$ 为海上升压站的容量。

3 海上风电汇集系统模型求解算法

收起

风机数量多, 主要通过线缆和汇集电站汇集到海上升压站, 其物理接线形式类似一个树形结构,将海上风机和电缆分别抽象成顶点和边,由于海上升压站是所有风机电能的汇集点, 因此对于以最小生成树算法为基础的遗传算法优化模型而言,设置海上升压站为生成树的根节点,有利于简化求解过程获得最优解。对于汇集系统而言, 拓扑结构主要由海缆成本所决定, 海缆成本与其型号相关,海缆型号与其传输的功率相关,传输功率由拓扑结构确定,这就形成了海缆成本和拓扑结构相互耦合的过程, 可以通过动态变权最小生成树 (Dynamic Variable Weight Minimum Spanning Tree, DVW_MST)算法对遗传算法进行改进, 求解得到海上风电场优化拓扑结构。

3.1 算法原理

使用动态变权最小生成树算法生成初始种群, DVW_MST 算法的流程如图 4 所示。

交叉判定的具体步骤: 每次使用 DVW_MST 进行选边时, 按照净收益变化量大小对全部可选边进行排序,之后从净收益变化量最大的边开始, 利用式(10)所给的交叉规避约束对模型的树形拓扑进行交叉识别,如果判定算法产生交叉线缆,则按照净收益大小排序依次寻找不产生交叉的次最优边；如果计算过程中未出现交叉线缆, 则将此边添加进树形结构中。

动态边权是 DVW_MST 算法最重要的环节之一, 建立动态边权最重要方法在于充分考虑海缆的建设成本和拓扑结构的耦合关系。算法过程首先需产生 3 个分别与树形拓扑相关的集合:集合 $U$ 表示还没有纳入树形结构的风机、节点集合 $S$ 中包括所有已经接入到树形结构中的风机、已经加入树形结构的边则存入节点集合 $E$ 。具体 DVW_MST 算法大致可分解为以下几个步骤: ① 输入海上输电网络完全图 $G$ 中的风机坐标作为顶点集合 $V$ 以及输入输电线缆作为边的集合 $E$ , 并初始化树形拓扑节点集合 $S$ ,此时 $S$ 只存入了以海上升压站为根节点的元素, 即第一个元素为 1 的集合,用数学关系表达为 $U =\{ i \mid i \in V, i \neq 1\}$ , $H \neq \Phi$ ; 同时初始化邻接矩阵 ${\mathbf{A}}_{\mathrm{T}}$ 为 0 ; ②在节点集合 $U$ 中把一个未加入树形拓扑的第 $i$ 个风机与树形拓扑集合 $S$ 中的第 $j$ 个风机关联,通过使用集电拓扑成本计算公式, 建立连接后的成本变化量 $\Delta {C}_{ij}$ ; ③通过比较挑选出 $\Delta {C}_{ij}$ 最小的顶点和相应连接的边,分别加入树形拓扑集合 $S$ 和集合 $H$ 中,生成新树的邻接矩阵 ${\mathbf{A}}_{\mathrm{T}}$ 。重复以上步骤,不断迭代更新集合 $U, S, E$ ,直到遍历完所有完全图中的节点,所有集合 $U$ 中的点清空并且风机坐标全都加入到了树形结构中,判断标准为点集 $S$ 和点集 $V$ 相等；④生成最小树,得到包括净收益和邻接矩阵 $A$ 的输出结果,结束程序。

本文设计了一种新的链表编码方式, 通过生成最小树形结构求解问题, 反映了海上风电的优化问题, 编写过程兼具了完整性和精确性。矩阵 ${\mathbf{A}}_{m \times n}$ 用来表示一个含有 $m$ 行元素、 $n$ 个节点的生成树。在本文的优化模型中 $m = 3$ ,第一行元素表示节点的连接情况, 第二行表示节点负载的风机个数,第三行包含了两节点之间电缆的长度。在生成树矩阵中,每个“个体”,即每个风机点都有与之对应的信息。例如,生成树矩阵第一行第 $i$ 个元素如果等于 $j$ ,则表示第 $i$ 个海缆节点与第 $j$ 个海缆节点直接相连, 并且电能的流通是单向的, 通过节点 $i$ 的电能必须流过节点 $j$ 才能到达海上升压站；第二行第 $i$ 个元素等于 $k$ ,表示节点之间的负载有 $k$ 个风机; 同理第三行第 $i$ 个元素等于 $d$ ,表示节点间距为 $d$ 。

为了判断遗传算法中种群个体的表现, 需定义适应度函数以判断个体表现的优劣性。在本文中, 成本较低的个体定义为优秀, 成本较高的个体则为低劣。因此适应度函数可表示为成本的倒数, 即:

(12)$ f\left({x}_{i}\right)= 1/{E}_{\text{cost }}$

本文使用精英选择、轮盘选择两种算子。精英选择保留最好种群中适应度函数最高的个体,直接进入下一次迭代, 这样遗传算法的趋优性可以得到保障,并且加快遗传算法的收敛速度。

交叉算子是将两个个体进行叠加得到可行解空间, 再利用 DVW_MST 算法在可行解空间中生成子树, 从而得到子代。交叉算子的求解过程是模拟两个父代个体的染色体信息遗传给子代。在本文中父代的信息可以由邻接矩阵表示, 任意两个父代生成树链表信息可表达为邻接矩阵。个体 ${\mathbf{A}}_{1}\left({3, n}\right)$ 和 ${\mathbf{A}}_{2}\left({3, n}\right)$ 的邻接矩阵 ${\mathbf{B}}_{1}\left({n, n}\right)$ 和 ${\mathbf{B}}_{2}(n$ , $n)$ 交叉重组得到结果 $\mathbf{B}\left({n, n}\right)$ ,再利用 DVW_MST 生成子代种群。

变异算子是在个体邻接矩阵中随机选择 $k$ 条边加入邻接矩阵中, 再利用动态变权最小生成树算法重新生成子代。父代个体 $\mathbf{A}\left({n,3}\right)$ 的遗传信息同样可以转化为邻接矩阵 $\mathbf{B}\left({n, n}\right)$ 。为了实现个体遗传信息的突变过程, 在生成子代之前, 在此个体邻接矩阵中加入任意条边 $p$ 数得到。

3.2 计算流程

结合海上风电场的经济成本模型, 采用上述的基于动态变权最小生成树的改进遗传算法, 对海上风电场集群的汇集系统进行规划求解, 分别进行初始种群生成、编码、适应度计算、选择、交叉、变异等操作,详细求解过程如图 5 所示。首先输入海上升压站和风电场坐标生成初始种群, 其次分别计算各个风电场的净收益, 通过适应度函数的处理, 判断是否满足终止准则, 如果不满足, 则一步选取两个个体进行交叉变异,利用 DVW_MST 算法生成子树得到子代, 重新迭代求解；如果满足, 则输出种群作为最优个体, 纳入最优解并进一步产生效果更好的下一代种群。

4 案例分析

收起

4.1 海上风机拓扑优化

为了验证本文所建模型和所提算法的有效性, 选取江苏某地区正在规划的 3 种风电场拓扑结构为例,其集电海缆经济参数如表 1 所示 ^{[ 16 ]} 。拓扑结构如图 6 所示。

针对上述 3 种拓扑结构进行优化, 优化后的拓扑见图 7 ,优化结果如表 2 ,3 所示。

通过优化结果可知, 拓扑 3 的结构未发生变化, 所需成本均为 3.6369 亿元, 验证了本文所提智能算法与实际工程经验结果一致。拓扑 1,2 的优化结果与现有规划方案不同, 在相同功率传输的约束下, 本文所提优化算法成本更低, 资产利用率更高。通过本文所提优化算法, 拓扑 1 增加了 ${71.8}\mathrm{\;{km}}$ 截面积为 ${400}{\mathrm{\;{mm}}}^{2}$ 的大功率传输电缆,但是减少了 ${107.8}\mathrm{\;{km}}$ 的小功率传输电缆,提升了现有传输线路的资产利用率,并且降低了工程总投资,线缆成本减少了 209 万元,提升了系统经济性。拓扑 2 通过优化风机分组,提升了汇集电缆传输功率,虽然增加了 ${56.6}\mathrm{\;{km}}$ 的较小截面积电缆,但是减少了 ${40.9}\mathrm{\;{km}}$ 截面积为 ${400}{\mathrm{\;{mm}}}^{2}$ 的大功率传输电缆,增加了线缆利用率,线缆成本减少了 43 万元。说明本文所提出的方法能够对风电场风机拓扑进行有效优化, 提高设备利用率, 降低成本。

4.2 风电集群拓扑优化

本文以江苏某地区 7 个海上风电场为对象, 相关参数如表 1 。对 7 个风电场的经纬度进行换算, 得到风电场的地理位置信息, 陆上升压站设在 (100,100)处,具体如表 4 所示。

风电场的出力信息利用风电场历史出力数据进行截取, 作为其场景年的年出力数据, 同时也反映了风电场集群内部的相关性信息。同时给定输电价格及投资信息:输电价格为 0.08 元/( $\mathrm{{kW}}\cdot \mathrm{h}$ ), 因输电有功功率限制而导致风电场弃电的补偿电价为 0.32 元 $/\left({\mathrm{{kW}}\cdot \mathrm{h}}\right)$ ,集电系统使用寿命为 ${20}\mathrm{a}$ , 基准折现率为 0.08。

使用本文所提出的基于 DVW_MST 的改进遗传算法对优化模型进行求解, 海上升压站位置经过优化得到坐标为(0,35),拓扑优化结果如图 8 所示。图中,风车为风电场,黑色方块分别为海上升压站和陆上升压站。

将各个风电场直接接入海上升压站默认坐标 (0,0),再通过海上升压站与陆上升压站的连接方式与本文优化后的结果进行对比,对比结果如表 5 所示。

由表 5 可知, 使用本文的改进遗传算法后, 集电系统线缆长度由 ${563.1}\mathrm{\;{km}}$ 减少到 ${416.3}\mathrm{\;{km}}$ ,通过适当减少截面积较小的 ${240}{\mathrm{\;{mm}}}^{2}$ 和 ${800}{\mathrm{\;{mm}}}^{2}$ 的线缆,改用截面积更大的 ${400}{\mathrm{\;{mm}}}^{2}$ 和 ${1600}{\mathrm{\;{mm}}}^{2}$ 的截面线缆, 能够对各风电场出力进行汇集输送, 避免为了满足单个风电场传输需求而采用多条小功率电缆, 从而提高线路利用率。改进风电场的汇集方式后,系统成本从 31 亿元降低至 25.78 亿元,成本降低了 16.84%,集电系统的收益从 10.35 亿元提升至 15.5 亿元,总体收益提高 5.15 亿元。

最后,为了验证本文所提改进遗传算法的寻优寻优效果, 采用传统遗传算法对上述优化模型进行求解, 统计不同迭代次数中得到最优解的次数如图 9 所示。

由图 9 可知,与传统的遗传算法相比,本文所提的改进遗传算法在相同迭代次数中, 获得最优解的次数明显更多。

5 结论

收起

本文考虑了海上风电出力在空间上的互补特性, 并综合考虑多种约束条件, 提出了一种海上风电场汇集系统拓扑优化方法。算例结果表明:①多个风电场的平滑效应使风电场集群的最大出力降低, 持续出力曲线更加平缓; ②风电场或风电场集群的送出线路的优化容量通常低于其装机容量; ③本文优化后得到的拓扑结构相比传统直接连接所得拓扑结构成本明显降低,净收益显著提高。验证了改进遗传算法的可行性, 体现了本文所提出模型的有效性。

基金

收起

国网江苏省电力有限公司咨询研究项目(B710K02135JB)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

刘万琨, 张志英, 李银凤, 等. 风能与风力发电技术[M]. 北京: 化学工业出版社, 2006.

[2]

王恰. 2020-2060年中国风电装机规模及其 CO₂ 减排预测[J]. 生态经济, 2021, 37(7): 13-21.

[3]

Yang

, Bai

. Development and prospects of offshore wind power[A]. 2010 World Non-Grid-Connected Wind Power and Energy Conference[C].Nanjing:WNWEC, 2010.1-6.

[4]

穆永铮, 鲁宗相, 申洪. 汇集站和集群中心站分级协调的风场群接入系统优化规划[J]. 中国电机工程学报, 2014, 34(S1): 1-8.

[5]

孙瑞娟, 梁军, 王克文, 等. 海上风电集电系统研究综述[J]. 电力建设, 2021, 42(6): 105-115.

[6]

赵彪, 安峰, 屈鲁, 等. 多功能直流集电器概念及其全直流海上风电系统[J]. 中国电机工程学报, 2021, 41(18): 39-51.

[7]

荣飞, 徐业事, 黄守道, 等. 新型直流串联型海上风电场拓扑结构设计[J]. 电力系统及其自动化学报, 2018, 30(2): 36-42.

[8]

黄永生, 黄玲玲. 大型海上风电场风机微观选址与集电系统联合优化[EB/OL]. http://kns.cnki.net/kcms/detail/23.1202.TH.20210419.0905.002.html2021-08-20. http://kns.cnki.net/kcms/detail/23.1202.TH.20210419.0905.002.html

[9]

汪惟源, 乔颖, 窦飞, 等. 基于改进遗传算法的海上风电场集电系统拓扑优化[J]. 中国电力, 2019, 52(1): 63-68.

[10]

Huang

, Chen

, Zhang

, et al. Optimization of large-scale offshore wind farm electrical collection systems based on improved FCM[A]. International Conference on Sustainable Power Generation and Supply[C].Hangzhou: IEEE, 2012.1-5.

[11]

吴伊雯. 基于混合生态共生算法的海上风电场集电线路优化方法研究[D]. 北京: 华北电力大学, 2020.

[12]

樊潇, 卢永魁, 黄玲玲, 等. 大型海上风电场集电系统网络拓扑优化设计[J]. 电力系统及其自动化学报, 2016, 28(7): 51-56.

[13]

谭任深. 海上风电场集电系统的优化设计[D]. 广州: 华南理工大学, 2013.

[14]

车文学. 海上风电场微观选址及集电系统优化研究[D]. 北京: 华北电力大学, 2017.

[15]

刘燕华, 田茹, 张东英, 等. 风电出力平滑效应的分析与应用[J]. 电网技术, 2013, 37(4): 987-991.

[16]

黄明煌, 王秀丽, 刘沈全, 等. 分频输电应用于深远海风电并网的技术经济性分析[J]. 电力系统自动化, 2019, 43(5): 167-174.

2024年第42卷第1期

PDF下载

281

146

引用本文

BibTeX

文章信息

接收时间：2021-10-20
首发时间：2025-07-22
出版时间：2024-01-20

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2021-10-20

基金

国网江苏省电力有限公司咨询研究项目(B710K02135JB)

作者信息

¹ 国网江苏省电力有限公司经济技术研究院江苏南京 610095

² 清华四川能源互联网研究院四川成都 610000

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/kzsny/CN/1154432827316756515

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

截面积	传输容量/MW	价格/万元	挂载风机数
70	24	77.8	2
120	33	91.16	3
185	41	106.76	4
240	47	117.23	5
400	60	155.29	6
800	150	550
1000	200	734
1 200	300	887
1400	350	1200
1600	400	1500

截面积

传输容量/MW

价格/万元

挂载风机数

77.8

120

91.16

185

106.76

240

117.23

400

155.29

800

150

550

1000

200

734

1 200

300

887

1400

350

1200

1600

400

1500

	截面积mm	拓扑 1	拓扑 2
优化前	优化后	节约	优化前	优化后	节约
不同截面积对应线路长度/km	70	70	66.2	3.8	64.6	64.6	0
120	65.8	25	40.8	32.3	32.3	0
185	31.2	25	6.2	32.8	64	-31.2
240	77	20	57	36.5	61.9	-25.4
400	0	71.8	-71.8	104.5	63.6	40.9
总长度 $/\mathrm{{km}}$	244	208	36	270.7	286.4	-15.7

截面积mm

拓扑 1

拓扑 2

优化前

优化后

节约

优化前

优化后

节约

不同截面积对应线路长度/km

66.2

3.8

64.6

120

65.8

40.8

32.3

185

31.2

6.2

32.8

-31.2

240

36.5

61.9

-25.4

400

71.8

-71.8

104.5

63.6

40.9

总长度 $/\mathrm{{km}}$

244

208

270.7

286.4

-15.7

	截面积mm	拓扑 1	拓扑 2
优化前	优化后	节约	优化前	优化后	节约
不同截	70	5 446.00	5 150.36	295.64	5025.88	5 025.88	0.00
面积对	120	5 998.33	2 279.00	3 719.33	2944.47	2 944.47	0.00
应线路	185	3 330.91	2 669.00	661.91	3501.73	6 832.64	-3330.91
成本/	240	9026.71	2 344.60	6 682.11	4 278.90	7 256.54	-2977.64
万元	400	0.00	11 149.82	11 149.82	16 227.81	9 876.44	6 351.36
成本/亿元	2.38	2.36	0.0209	3.20	3.19	0.004 3

截面积mm

拓扑 1

拓扑 2

优化前

优化后

节约

优化前

优化后

节约

不同截

5 446.00

5 150.36

295.64

5025.88

5 025.88

0.00

面积对

120

5 998.33

2 279.00

3 719.33

2944.47

2 944.47

0.00

应线路

185

3 330.91

2 669.00

661.91

3501.73

6 832.64

-3330.91

成本/

240

9026.71

2 344.60

6 682.11

4 278.90

7 256.54

-2977.64

万元

400

0.00

11 149.82

16 227.81

9 876.44

6 351.36

成本/亿元

2.38

2.36

0.0209

3.20

3.19

0.004 3

编号			装机容量/MW
1	-53.09	51.11	10
2	-48.16	51.39	20
3	-33.27	-24.20	30
4	-17.11	-22.55	20
5	-13.16	51.29	30
6	-19.80	-45.94	10
7	100.49	-31.29	20

编号

装机容量/MW

-53.09

51.11

-48.16

51.39

-33.27

-24.20

-17.11

-22.55

-13.16

51.29

-19.80

-45.94

100.49

-31.29

截面积mm	优化前	优化后	节约
长度/km	成本/万元	长度/km	成本/万元	长度/km	成本/万元
240	123.70	14501.35	148.90	17 455.55	-25.20	-2954.20
400	203.90	31 663.63	50.90	7904.26	153.00	23 759.37
800	94.10	51 755.00	97.20	53 460.00	-3.10	-1705.00
1600	141.40	212 100.00	119.30	178 950.00	22.10	33 150.00
总计	563.10	310 019.98	416.30	257 769.81	146.80	52 250.17
收益/亿元		10.35		15.50		-5.16

截面积mm

优化前

优化后

节约

长度/km

成本/万元

长度/km

成本/万元

长度/km

成本/万元

240

123.70

14501.35

148.90

17 455.55

-25.20

-2954.20

400

203.90

31 663.63

50.90

7904.26

153.00

23 759.37

800

94.10

51 755.00

97.20

53 460.00

-3.10

-1705.00

1600

141.40

212 100.00

119.30

178 950.00

22.10

33 150.00

总计

563.10

310 019.98

416.30

257 769.81

146.80

52 250.17

收益/亿元

10.35

15.50

-5.16