CN106951619A - 计及可靠性约束的海上风电场拓扑设计方法 - Google Patents

计及可靠性约束的海上风电场拓扑设计方法 Download PDF

Info

Publication number
CN106951619A
CN106951619A CN201710144999.0A CN201710144999A CN106951619A CN 106951619 A CN106951619 A CN 106951619A CN 201710144999 A CN201710144999 A CN 201710144999A CN 106951619 A CN106951619 A CN 106951619A
Authority
CN
China
Prior art keywords
path
topology
standard
ruskal
carrying capacity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710144999.0A
Other languages
English (en)
Other versions
CN106951619B (zh
Inventor
刘俊
郝旭东
赵康宁
王旭
陈业夫
苏灿
张�杰
刘佳钰
高欣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201710144999.0A priority Critical patent/CN106951619B/zh
Publication of CN106951619A publication Critical patent/CN106951619A/zh
Application granted granted Critical
Publication of CN106951619B publication Critical patent/CN106951619B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/06Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]

Abstract

本发明公开了计及可靠性约束的海上风电场拓扑设计方法,对标准Kruskal的每步搜索结果进行再约束,对其中超过载流量限制的连接线增加辅助线路保证线路载流量的均匀性;具体思路为:对标准Kruskal搜索得到的结果进行分析,找出其中是哪一条线路对该拓扑的载流量影响最大,进而选取其它线路或增加辅助线路来代替,最终使得拓扑中的每一条线路都满足载流量的要求;标准Kruskal算法求得的风电场拓扑虽然结构上显著减少了风电场内电缆线数量,但也引入了诸如可靠性与电缆线截面等问题,其在经济指标中相对于典型排布方式并无优势可言;本发明的方法则兼顾了较短的拓扑连线与良好的可靠性,其拓扑设计技术指标和经济指标都明显优于原传统海上风电场拓扑设计结构。

Description

计及可靠性约束的海上风电场拓扑设计方法
技术领域
本发明属于电源规划领域,具体涉及计及可靠性约束的海上风电场拓扑设计方法。
背景技术
风力发电凭借其清洁、低成本等特点,已经成为当前可再生能源的重要利用方式之一,越来越受到各个国家的重视。而相比于陆地风电,海上风电具有无需占用土地资源、能量效益高、运行稳定、还可建设于沿海经济发达地区以利于就近消化等诸多优点,利用前景非常广阔,已经成为未来风能利用的必然趋势。但同时,海上风电在利用上技术难度大、同时投资与维护成本高,因而对海上风电的技术研究也就显得尤为重要。在当前世界范围内,已经投入运营或在建的海上风电场一般容量不大、电气连接方式并不复杂,因此可以由设计人员结合海上风电场工程项目的实际需要根据经验设计其拓扑结构。而在海上风电场建设规模日益趋向大型化的未来,该种设计方式对设计人员的水平依赖将越来越强,不同的设计可能导致的误差也将越来越大。
根据海上风电场典型排布研究与电源规划理论,可以采用标准Kruskal算法来重新设计海上风电场的拓扑结构。针对风电场内部结构的特点,把风电场拓扑的问题转化为最小生成树的问题,利用标准Kruskal算法获得风机之间连接线最短的拓扑路径,并通过算例证明了风电场在采取Kruskal算法之后对原来的潮流分布影响不大,而降低连接线长度的效果十分显著。
标准Kruskal算法中出现的电缆线载流量分布极不均匀的问题,仍是海上风电场拓扑结果设计中亟待解决的难题。
发明内容
为了解决上述现有技术存在的问题,进一步提高风电场拓扑的综合价值,需在尽量保留Kruskal算法在连接线路长度上的优势基础上充分考虑如何均匀化线路的载流量,即需要对Kruskal算法进行改进,本发明提供计及可靠性约束的海上风电场拓扑设计方法,得到满足系统可靠性约束的拓扑结构。
为达到以上目的,本发明采用如下技术方案:
计及可靠性约束的海上风电场拓扑设计方法,对已经算得的标准Kruskal结果进行再处理,形成一种计及电缆线路载流量约束的改进Kruskal算法;依据此思路,能够对标准Kruskal算法求得的结果进行分析,找出其中对该拓扑的载流量影响最大的那条线路,进而选取其他线路或增加辅助线路代替,最终使得拓扑中的所有线路都满足载流量的要求;
标准Kruskal算法忽视了原本典型排布中各个接线可能的载流量限制与可靠性问题,由此引出对标准Kruskal算法的改进方法,提出了计及可靠性约束的改进Kruskal算法对海上风电场拓扑进行优化设计;所述计及可靠性约束的海上风电场拓扑设计方法,包括如下步骤:
步骤一:录入初始数据、载流量限制Cmax,进行标准Kruskal计算;标准Kruskal算法为:初始化计算数据,主要是录入所有节点的坐标与节点间的路径权重,形成初始路径权重集,并建立已选节点集与已选路径集供储存计算结果备用;根据上述初始化的结果,从路径权重集中挑选权重最小的路径加入已选路径集,将该路径连接起来的点加入已选节点集;对每一个已经由多个节点和路径组成起来的树,将除了已选路径以外其他节点之间的路径从路径权重集中删除,以防止形成环网;重复上述选取最小路径和更新路径权重集的步骤,直到所有节点都加入到已选节点集中;这时已选路径集中的路径即为该加权连通图最小生成树的所有路径;
步骤二:统计当前拓扑中每个线路载流量,找出一条载流量Ci≥Cmax的线路,标记该线路功率注入点为O;
步骤三:从初始路径权重集中分别选出能够通过O点及与O点相邻点将功率送出去的最短路径,将添加后功率更平均的线路加入拓扑中;
步骤四:若产生了多重环网,则去除环网的内部多余拓扑;
步骤五:判断对所有线路编号i,均有Ci<Cmax,如果不满足,返回步骤二;
步骤六:输出当前拓扑结果。
本发明方法并未对Kruskal算法的使用本身做出任何更改,而是对标准Kruskal算法计算出的拓扑根据原定线路载流量的约束做了少许必要更改使其满足条件,而又不至于像典型排布那样具有许多冗余线路。这种对标准Kruskal的拓扑更改的方法最主要的特点是保留了原来拓扑远端或相距非常近的点之间的树,使其在满足条件内的部分以最小生成树的形式继续留在拓扑中;而对于靠近汇流点的电缆线,基于标准Kruskal算法的拓扑使得远端的风机功率在送至汇流点的过程中会越来越集中,而该再约束算法会从靠近不满足线路最大载流量的线路节点开始对拓扑结构进行更改,筛选出其他可供输送功率的多条支路,选出其中最小的一条最为辅助路径加入到拓扑中,再去除辅助路径加入后可能出现的不必要的多重环网结构,以此逐步对每一条超过载流量的线路进行优化,最终得到满足系统可靠性约束的拓扑结构。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为标准Kruskal算法的流程框图。
图2为计及可靠性约束的Kruskal优化方式流程图。
图3为算例风机排布情况。
图4为典型拓扑。
图5为标准Kruskal算法拓扑。
图6为采用计及可靠性约束的Kruskal算法拓扑。
具体实施方式
以下结合附图及详细模型构建原理,对本发明作进一步的详细描述。
标准Kruskal算法的基本步骤为:初始化计算数据,主要是录入所有节点的坐标与节点间的路径权重,形成路径权重集,并建立已选节点集与已选路径集供之后储存计算结果备用;根据上述初始化的结果,从路径权重集中挑选权重最小的路径加入已选路径集,将该路径连接起来的点加入已选节点集;对每一个已经由多个节点和路径组成起来的树,把除了已选路径以外其他节点之间的路径从路径权重集中删除,以防止形成环网;重复上述选取最小路径和更新路径权重集的步骤,直到所有节点都加入到已选节点集中。这时已选路径集中的路径即为该加权连通图最小生成树的所有路径。具体流程图如图1。
算例风机排布情况如图3,该风电场典型拓扑如图4,经过标准Kruskal算法计算后拓扑如图5,图中数字为风机编号,横纵坐标为风机所在地理位置(适用于图3‐图6)。该算例经过标准Kruskal算法计算后风电场拓扑所需要的连接线路长度降低得十分明显,对成本降低十分有帮助。同时潮流计算的结果中,经过标准Kruskal算法计算后的风电场拓扑由于电缆线长度的降低,虽然线路载流量上升了,但整体电压升高量与典型排布相差无几,网络中各个节点依然易于控制在允许范围内,而由于电压相角对整个风电场网络运行影响不大,其微小偏移对电网稳定性几乎没有影响。
但与此同时,由于Kruskal算法的结果中所有风机只经过一条电缆线汇流到汇流母线进行升压,不同的电缆线路载流量差距巨大,潮流计算前提出的均一网假设很可能无法同时匹配两种拓扑方式对电缆线的要求。同时,标准Kruskal算法求得的拓扑可靠性相比于放射链式更差,其中关键的某一条或几条电缆故障后造成的损失可能比放射链式故障时更大。除此之外,采用标准Kruskal算法计算的拓扑结构中,由于线路载流量高,总功率损耗也有少量升高,抵消了一部分减少电缆的铺设线路所带来的成本降低,具体效果需要采用完善的经济比较模型进行进一步对比。因此总的来说,使用标准Kruskal算法是一个较为直观且简单的优化方法,效果一目了然,但在可行性方面有待商榷。
本发明设计计及可靠性约束的改进Kruskal算法进行优化,具体流程图如图2。对已经算得的标准Kruskal结果进行再处理。标准Kruskal算法中亟待解决的核心问题是线路载流量极为不均匀。依据此思路,可以对标准Kruskal算法求得的结果进行分析,找出其中是哪一条线路对该拓扑的载流量影响最大,进而选取其他线路或增加辅助线路代替,最终使得拓扑中的每一条线路都满足载流量的要求。
该方法并未对Kruskal算法的使用本身做出任何更改,而是对标准Kruskal算法计算出的拓扑根据原定线路载流量的约束做了少许必要更改使其满足条件,而又不至于像典型排布那样具有许多冗余线路。这种对标准Kruskal的拓扑更改的方法最主要的特点是保留了原来拓扑远端或相距非常近的点之间的树,使其在满足条件内的部分以最小生成树的形式继续留在拓扑中;而对于靠近汇流点的电缆线,基于标准Kruskal算法的拓扑使得远端的风机功率在送至汇流点的过程中会越来越集中,而该再约束算法会从靠近不满足线路最大载流量的线路节点开始对拓扑结构进行更改,筛选出其他可供输送功率的多条支路,选出其中最小的一条最为辅助路径加入到拓扑中,再去除辅助路径加入后可能出现的不必要的多重环网结构,以此逐步对每一条超过载流量的线路进行优化,最终得到满意的拓扑结构。采用计及可靠性约束的Kruskal算法对上述算例进行优化后结果如图6。
该方法在对Kruskal结果进行优化后,由于电路中形成了环网,每条线路的载流量无法再简单地用该线路所连接的风机数量表示,这里近似认为功率顺着环网的不同两条汇流路线平均分配。在汇流的过程中,每到一个载流量将要超过允许值的线路时,按照算法线路一般会在近处分为两条不同的传输路径并尽可能地平均带载,以此来尽量发挥每条电路的载流能力,可以想象,在节点数较多时,风电场拓扑将趋于网状。
本发明采用电力系统成本评估中常用的等年值法来对不同拓扑结构的组网方案进行经济比较。等年值法是选取等年值费用最小的方案为最佳方案,该方案将电网设备在使用寿命内的所有费用均换算成每年等值费用。
根据计算,可算得本发明中给出的几种拓扑结构固定投资成本如下表1。
表1不同拓扑结构固定投资成本表
这里根据线路最大载流量而选取了不同的电缆规格,可以看出,虽然标准Kruskal算法的电缆线长度很短,但由于其对电缆线横截面规格要求很高,实际投资成本并没有显著降低,而优化算法在投资成本上明显降低。
这几种拓扑结构的总年损耗费用如下表2。
表2不同拓扑结构年损耗费用表
根据计算,这几种拓扑结构年故障损失功率与总平均功率成本如下表3。
表3不同拓扑结构年平均功率成本表
根据以上计算结果可以看出,标准Kruskal在计算经济指标后,其平均功率成本反而不如典型放射链式排布,这是由于在光考虑其电缆线投资建设的同时忽略了可能影响经济指标的其他方面因素,如电缆线截面取值、由于拓扑结构不可靠带来的故障功率损失等。而对于经过优化后的拓扑结构而言,其平均功率成本较典型排布明显降低。
和现有技术相比较,本发明具备如下优点:
(1)计及可靠性约束的改进Kruskal算法是从标准Kruskal算法的结果入手,从单独连接的风机向里找出超过载流量限制的线路,在其附近增加辅助电缆线,分担功率输送的同时在风电场拓扑中形成环网,加强了整个风电场拓扑的稳定性。
(2)本发明解决了以下问题:标准Kruskal算法求得的风电场拓扑虽然结构上显著减少了风电场内电缆线数量,但由于带来的诸如可靠性与电缆线截面等问题,其在经济指标中相对于典型排布方式并无优势可言。而将标准Kruskal算法改进后的优化拓扑结构则兼顾了拓扑连线短与良好的可靠性,无论从拓扑直观结果还是从经济指标的对比来看都优于原传统排布结构,值得在未来的工程中作为参考。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案,本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都属于本发明的保护范围。

Claims (1)

1.计及可靠性约束的海上风电场拓扑设计方法,其特征在于:对已经算得的标准Kruskal结果进行再处理,形成一种计及电缆线路载流量约束的改进Kruskal算法;依据此思路,能够对标准Kruskal算法求得的结果进行分析,找出其中对该拓扑的载流量影响最大的那条线路,进而选取其他线路或增加辅助线路代替,最终使得拓扑中的所有线路都满足载流量的要求;
所述计及可靠性约束的海上风电场拓扑设计方法,包括如下步骤:
步骤一:录入初始数据、载流量限制Cmax,进行标准Kruskal计算;标准Kruskal算法为:初始化计算数据,主要是录入所有节点的坐标与节点间的路径权重,形成初始路径权重集,并建立已选节点集与已选路径集供储存计算结果备用;根据上述初始化的结果,从路径权重集中挑选权重最小的路径加入已选路径集,将该路径连接起来的点加入已选节点集;对每一个已经由多个节点和路径组成起来的树,将除了已选路径以外其他节点之间的路径从路径权重集中删除,以防止形成环网;重复上述选取最小路径和更新路径权重集的步骤,直到所有节点都加入到已选节点集中;这时已选路径集中的路径即为该加权连通图最小生成树的所有路径;
步骤二:统计当前拓扑中每个线路载流量,找出一条载流量Ci≥Cmax的线路,标记该线路功率注入点为O;
步骤三:从初始路径权重集中分别选出能够通过O点及与O点相邻点将功率送出去的最短路径,将添加后功率更平均的线路加入拓扑中;
步骤四:若产生了多重环网,则去除环网的内部多余拓扑;
步骤五:判断对所有线路编号i,均有Ci<Cmax,如果不满足,返回步骤二;
步骤六:输出当前拓扑结果。
CN201710144999.0A 2017-03-13 2017-03-13 计及可靠性约束的海上风电场拓扑设计方法 Expired - Fee Related CN106951619B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710144999.0A CN106951619B (zh) 2017-03-13 2017-03-13 计及可靠性约束的海上风电场拓扑设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710144999.0A CN106951619B (zh) 2017-03-13 2017-03-13 计及可靠性约束的海上风电场拓扑设计方法

Publications (2)

Publication Number Publication Date
CN106951619A true CN106951619A (zh) 2017-07-14
CN106951619B CN106951619B (zh) 2020-01-17

Family

ID=59467302

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710144999.0A Expired - Fee Related CN106951619B (zh) 2017-03-13 2017-03-13 计及可靠性约束的海上风电场拓扑设计方法

Country Status (1)

Country Link
CN (1) CN106951619B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109104328A (zh) * 2018-06-14 2018-12-28 北京航空航天大学 一种基于复杂网络的集群系统可靠性测评方法
CN110059885A (zh) * 2019-04-23 2019-07-26 南京林业大学 一种光纤铺设路径的优化方法
CN112131770A (zh) * 2020-09-15 2020-12-25 北京化工大学 一种考虑可靠性的功能梯度连续体结构轻量化设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010029548A1 (en) * 2008-09-11 2010-03-18 Yissum Research Development Company Of The Hebrew University Of Jerusalem, Ltd. Method for identifying genetic loci invovled in hybrid vigor
CN102760195A (zh) * 2012-08-01 2012-10-31 中国能源建设集团广东省电力设计研究院 海上风电场集电系统拓扑设计方法及装置
CN103219725A (zh) * 2013-01-18 2013-07-24 云南电网公司大理供电局 一种基于实时运行数据的风电场等值建模方法
CN104281737A (zh) * 2014-08-06 2015-01-14 华南理工大学 一种海上风电场集电系统拓扑结构选型方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010029548A1 (en) * 2008-09-11 2010-03-18 Yissum Research Development Company Of The Hebrew University Of Jerusalem, Ltd. Method for identifying genetic loci invovled in hybrid vigor
CN102760195A (zh) * 2012-08-01 2012-10-31 中国能源建设集团广东省电力设计研究院 海上风电场集电系统拓扑设计方法及装置
CN103219725A (zh) * 2013-01-18 2013-07-24 云南电网公司大理供电局 一种基于实时运行数据的风电场等值建模方法
CN104281737A (zh) * 2014-08-06 2015-01-14 华南理工大学 一种海上风电场集电系统拓扑结构选型方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109104328A (zh) * 2018-06-14 2018-12-28 北京航空航天大学 一种基于复杂网络的集群系统可靠性测评方法
CN109104328B (zh) * 2018-06-14 2020-09-18 北京航空航天大学 一种基于复杂网络的集群系统可靠性测评方法
CN110059885A (zh) * 2019-04-23 2019-07-26 南京林业大学 一种光纤铺设路径的优化方法
CN110059885B (zh) * 2019-04-23 2022-01-14 南京林业大学 一种光纤铺设路径的优化方法
CN112131770A (zh) * 2020-09-15 2020-12-25 北京化工大学 一种考虑可靠性的功能梯度连续体结构轻量化设计方法
CN112131770B (zh) * 2020-09-15 2023-12-15 北京化工大学 一种考虑可靠性的功能梯度连续体结构轻量化设计方法

Also Published As

Publication number Publication date
CN106951619B (zh) 2020-01-17

Similar Documents

Publication Publication Date Title
CN103514570B (zh) 一种含分布式电源的配电网扩展规划综合优化方法
CN106503839B (zh) 一种海上风电场环形集电网络分层规划方法
WO2016141683A1 (zh) 适用于消纳大规模分布式电源的能源互联网的实现方法
CN106951619A (zh) 计及可靠性约束的海上风电场拓扑设计方法
CN107706893B (zh) 一种dg接入配网最佳容量的计算方法
CN102760195A (zh) 海上风电场集电系统拓扑设计方法及装置
CN105846422B (zh) 一种大容量电网短路电流控制措施优化配置方法
CN104200263A (zh) 一种基于禁忌差分进化和gis的配电网络线路规划方法
CN105528653A (zh) 考虑经济性和安全性的光伏接入优化方法
Hopewell et al. Optimising the design of offshore wind farm collection networks
Zuo et al. Collector system topology design for offshore wind farm's repowering and expansion
CN105656076A (zh) 海上风电场海缆布线的获取方法
CN109598377A (zh) 一种基于故障约束的交直流混合配电网鲁棒规划方法
CN109740975A (zh) 考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法
CN105740966A (zh) 一种含分布式电源的配电网扩展规划方法
CN114662799A (zh) 一种输电线路检修计划优化方法、装置、设备及存储介质
CN106410842B (zh) 一种基于全寿命周期的海上风电场集电系统成本计算方法
CN111799845B (zh) 电源集群送出的电力系统安控策略自适应建模方法及系统
CN106849068B (zh) 基于分区组合Kruskal算法的海上风电场拓扑设计方法
CN102426665A (zh) 一种配电网规划建设标准的选择方法
CN104809340A (zh) 一种多端柔性直流输电系统运行点的快速计算方法
CN110197302B (zh) 一种考虑接线模式拓扑约束的配电网规划方法
CN103258301A (zh) 一种适用于县域电网的供电可靠性评估方法
CN104821590A (zh) 智能电网节能优化综合测评方法
CN112398129B (zh) 电网事故后风险控制辅助决策方法、装置和存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200117

CF01 Termination of patent right due to non-payment of annual fee