CN111310972A - 一种考虑尾流效应的海上风电机组维护路径随机规划方法 - Google Patents

一种考虑尾流效应的海上风电机组维护路径随机规划方法 Download PDF

Info

Publication number
CN111310972A
CN111310972A CN202010052324.5A CN202010052324A CN111310972A CN 111310972 A CN111310972 A CN 111310972A CN 202010052324 A CN202010052324 A CN 202010052324A CN 111310972 A CN111310972 A CN 111310972A
Authority
CN
China
Prior art keywords
maintenance
node
ship
unit
under
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
CN202010052324.5A
Other languages
English (en)
Other versions
CN111310972B (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.)
Shanghai University of Electric Power
Original Assignee
Shanghai University of Electric Power
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 Shanghai University of Electric Power filed Critical Shanghai University of Electric Power
Priority to CN202010052324.5A priority Critical patent/CN111310972B/zh
Publication of CN111310972A publication Critical patent/CN111310972A/zh
Application granted granted Critical
Publication of CN111310972B publication Critical patent/CN111310972B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/047Optimisation of routes or paths, e.g. travelling salesman problem
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种考虑尾流效应的海上风电机组维护路径随机规划方法,包括以下步骤:1)描述海洋天气中风、浪、涌的相关性和不确定性;2)考虑人员、船只、惩罚和停机损失的各项成本构成,以维护成本最小作为目标函数,并基于风、浪、涌对船只航行时间以及出海可及窗口的影响构建对应的约束条件;3)由于维护状态会改变机组间尾流分布,结合尾流模型与维护状态描述各机组输入风速的变化,并精细化目标函数中的机组停机损失;4)对目标函数进行求解,得到最优的维护路径。与现有技术相比,本发明具有适应性强、经济性高、求解速度快等优点。

Description

一种考虑尾流效应的海上风电机组维护路径随机规划方法
技术领域
本发明涉及海上风电场检修调度领域,尤其是涉及一种考虑尾流效应的海上风电机组维护路径随机规划方法。
背景技术
海上风电相比于陆上风电,具有平均风速高,发电利用小时数高等显著优势,近年来发展迅速。然而受到海洋环境的制约,海上运行维护难度和成本均较大。据资料统计表明,海上风电机组的运维成本约占全生命周期成本的40%,是陆上风电运维成本的2-3倍。并且随着海上风电场的规模化发展和远离海岸的建设,运维船在风电场内的航行时间增加,进出风电场的路径变长,这使得海上运维成本进一步增加。所以合理地优化维护路径,降低运维成本是目前海上风电运维管理的重要要求。
在开展海上风电维护工作时,风、浪、涌等因素会影响运维船的出海可及窗口和航行时间,是制约维护作业开展的重要因素。一些文献在制定维护策略时,为了简化计算仅计及风速对海上风电场可及性的影响,对天气状态进行划分,以此对维护的天气等待时间进行计算。此外,也有一些文献考虑风速和浪高这两种主要天气因素,采用马尔科夫法分别描述两种因素的时序变化关系,并构建了可及度指标用于评估海上风电机组维修的可及性。这些研究中的风、浪等天气因素大多按照预测值或转化为可及因子,构建了确定性的维护模型。但是,考虑到受海洋环境的影响,风、浪、涌具有较强的相关性和难以预测性,确定性的维护策略无法适应所有可能的情况。为此,综合考虑风、浪、涌三者相关性和不确定性的维护模型仍需进一步研究。
海上风电机组的维护路径规划属于短期维护策略,是指在满足各类运维约束下,构建运维船的最优维护路径以实现成本最低。现有停机损失模型均忽略了风速变化对停机损失的影响,在短期的维护路径规划中,这种模型较为粗糙。同时由于海洋占地资源的有限性,风电机组之间具有较强的尾流效应。因此,为了准确评估停机损失对路径决策制定的影响,需进一步研究在风速变化下考虑尾流效应的停机损失模型。
此外,由于海上风电机组维护时需要考虑船只、备件、人员、天气的不确定性等诸多因素,且运维船维护路径的组合会随着维护任务的完成而发生改变,处于维护状态的机组将使得机组间尾流分布发生变化,下风向机组的输入风速受到影响,这使得短期维护下路径的求取变得非常复杂,这是一个包含连续变量和离散变量的多约束非线性优化问题。
因此,急需一种考虑尾流效应的海上风电机组维护路径随机规划方法,能够针对复杂的海洋环境,描述风、浪、涌相关性和随机性,同时能考虑到维护状态与机组间尾流分布的关联性。建立相应的维护路径随机规划模型,来为风电场维护获取更优的经济效益并且为决策者提供更多决策支持。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种考虑尾流效应的海上风电机组维护路径随机规划方法。
本发明的目的可以通过以下技术方案来实现:
一种考虑尾流效应的海上风电机组维护路径随机规划方法,包括以下步骤:
1)描述海洋天气中风、浪、涌的相关性和不确定性;
2)考虑人员、船只、惩罚和停机损失的各项成本构成,以维护成本最小作为目标函数,并基于风、浪、涌对船只航行时间以及出海可及窗口的影响构建对应的约束条件;
3)由于维护状态会改变机组间尾流分布,结合尾流模型与维护状态描述各机组输入风速的变化,并精细化目标函数中的机组停机损失;
4)对目标函数进行求解,得到最优的维护路径。
所述的步骤1)中,采用Copula函数描述风、浪、涌的相关性,并且采用多场景模型描述各时刻下风、浪、涌的不确定性,则有:
E(x1,x2,x3)=C[F1(x1),F2(x2),F3(x3)]
其中,E(x1,x2,x3)为风、浪、涌随机变量x1、x2、x3的联合分布函数,C[·]为Copula函数,F1(x1),F2(x2),F3(x3)分别为风、浪、涌随机变量x1、x2、x3的边缘分布。
所述的步骤2)中,目标函数C的表达式为:
minC=cq+cb+cp+cs
Figure BDA0002371621820000031
Figure BDA0002371621820000032
Figure BDA0002371621820000033
Figure BDA0002371621820000034
其中,cq为维护人员成本,cb为船只成本,cp为惩罚成本,cs为停机损失,Ω为场景总数,ω为场景编号,D为总天数,d为天数编号,L为运维船只数,l为运维船编号,H为技术人员类型总数,h为技术人员类型编号,πω为场景ω的概率,Ql,d,h为运维船l在第d天离开节点0时船上类型为h的技术人员数量,
Figure BDA0002371621820000035
为场景ω下类型为h的单个技术人员每天所需费用,μi,h为维护节点i所需要类型为h的技术人员数量,
Figure BDA0002371621820000036
为场景ω下单个维护人员单位时间的海上等待成本,
Figure BDA0002371621820000037
Figure BDA0002371621820000038
分别为场景ω下运维船l在第d天到达放下的节点i和接取的节点n+i的时间段,τ为人员和备件从运维船转移到节点机组所在处所需时间段,
Figure BDA00023716218200000312
为维护节点i所需时间段,Z-={1,2,3,...,n}为放下节点的集合,Z+={n+1,n+2,...,2n}为接取节点的集合,则Z=Z-∪Z+,并将一天内维护作业开始和结束的节点分别记为节点0和节点2n+1,即港口;
N为所有待维修机组所在节点的集合,其包括港口且i、j∈N,
Figure BDA0002371621820000039
为场景ω下运维船l在第d天海上航行的单位时间成本,Tω,l,i,j,d为场景ω下运维船l在第d天从节点i航行到节点j所需时间,xl,i,j,d为路径决策变量,即当船l在第d天从节点i驶向节点j的路径被选取时,其值为1,否则为0,ξj表示是否需要船只驻留,若节点j维护时不需要船只驻留,则ξj为1,否则为0,
Figure BDA00023716218200000310
为场景ω下船l在第d天在节点j维护时所需的驻留成本;
yi为机组所在节点i超过推荐的维护时间后延迟维护的天数,
Figure BDA00023716218200000311
为场景ω下机组所在节点i延迟维护后每天的惩罚成本;
t为时段编号,λi为机组所在节点i待维护部件故障停机的概率,cf为海上风电的上网电价,Pω,i,t为场景ω下机组所在节点i在时段t下正常运行时的输出功率,Δt为时段t的时长,
Figure BDA0002371621820000041
为场景ω下由于机组停机造成的损失电量。
所述的场景ω下运维船l在第d天从节点i航行到节点j所需时间Tω,l,i,j,d根据考虑风浪涌海洋天气建模,具体为:
Figure BDA0002371621820000042
Figure BDA0002371621820000043
Figure BDA0002371621820000044
其中,
Figure BDA0002371621820000045
为自然运行条件下运维船l从节点i到节点j的航行时间,aω,l,i,j,d为场景ω下船l在第d天从节点i到节点j适宜航行的概率,γ为风浪涌下的航行时间综合影响系数,
Figure BDA0002371621820000046
为场景ω下船l在第d天从节点i到节点j航行所需各时段下的海上平均风速,
Figure BDA0002371621820000047
为场景ω下时段t的海上风速。
所述的步骤2)中,约束条件具体包括:
维护任务约束:
在维护调度周期内,待维护机组都要得到维护,并确保每台待维护的机组只能放下和接取运维人员一次,且每艘运维船每天仅离开和返回港口一次,则有:
Figure BDA0002371621820000048
Figure BDA0002371621820000049
Figure BDA00023716218200000410
节点交通流守恒约束:
Figure BDA00023716218200000411
Figure BDA0002371621820000051
船只驻留约束:
Figure BDA0002371621820000052
其中,Zl为执行维护作业时需要运维船存在的机组所在节点的集合;
接取时间约束:
Figure BDA0002371621820000053
航行时间兼容性约束:
Figure BDA0002371621820000054
Figure BDA0002371621820000055
备件约束:
Figure BDA0002371621820000056
Figure BDA0002371621820000057
Figure BDA0002371621820000058
其中,Rj为维护机组j所需的备件重量,
Figure BDA0002371621820000059
为运维船l的负载容量,
Figure BDA00023716218200000510
为运维船l在第d天离开节点i后运维船上的备件重量;
人员约束:
Figure BDA00023716218200000511
Figure BDA00023716218200000512
(ql,i,d,hj,h-ql,j,d,h)xl,i,j,d=0,
Figure BDA00023716218200000513
(ql,i,d,hj,h-ql,j,d,h)xl,i,j,d=0,
Figure BDA00023716218200000514
其中,ql,i,d,h、ql,j,d,h为分别运维船l在第d天离开节点i、节点j后船上类型为h的技术人员数量,ql,0,d,h为运维船l在第d天离开港口节点0后船上类型为h的技术人员数量,
Figure BDA00023716218200000515
为运维船l的人员容量,φd,h为第d天可用的类型为h的技术人员数量,μj,h为节点j所需要类型为h的技术人员数量;
操作窗口约束:
Figure BDA0002371621820000061
其中,
Figure BDA0002371621820000062
为场景ω下运维船l在第d天的最大可工作时间段;
延迟维护约束:
Figure BDA0002371621820000063
Figure BDA0002371621820000064
式中:
Figure BDA0002371621820000065
为推荐维护机组j的最迟天数,yj为机组j超过推荐的维护时间后延迟维护的天数。
所述的步骤3)中,精细化后的机组停机损失中的损失电量Wω s表示为:
Figure BDA0002371621820000066
其中,NT为总时段数,
Figure BDA0002371621820000067
为在风电场无维护任务时场景ω下在时段t机组所在节点r的输出功率,Iω,r,t为场景ω下在时段t机组所在节点r的维护状态,取值为1时表示维护停机,取值为0时表示正常运行,S为所有的机组所在节点总数。
所述的步骤4)具体包括以下步骤:
对各场景下各船只每日最大可维护机组数ηl,d进行预评估,在完成所有场景下ηl,d的预评估后,选取各ηl,d中的最小值代入每日最大可维护机组数与维护路径决策变量之间的关系式,对运维船每日的维护路径产生约束,并完成模型求解。
各场景下各船只每日最大可维护机组数ηl,d的预评估具体为:
Figure BDA0002371621820000068
其中,
Figure BDA0002371621820000069
为考虑航行时间和机组维护时间并在计及风、浪、涌的出海可及窗口下的最大可维护机组数;
Figure BDA00023716218200000610
为满足船只容量下的最大可维护机组数;
Figure BDA00023716218200000611
为满足备件和载重量下的最大可维护机组数。
每日最大可维护机组数与维护路径决策变量之间的关系式具体为:
Figure BDA00023716218200000612
与现有技术相比,本发明具有以下优点:
一、适应性强:在维护计划周期内考虑海上天气环境风、浪、涌的相关性和不确定性制定的维护路径决策能更好地应对随机的海洋天气环境,提高了模型路径决策结果的适应性,能为决策者提供更加切实可行的方案。
二、经济性高:在制定维护路径决策时,考虑机组间的尾流效应能更准确地刻画机组停机损失,从而更确切地描述其对路径决策制定的影响,以此得到的维护路径结果更具经济性。
三、求解速度快:所提基于每日可维护机组数目预评估的求解方法,能够有效地约束维护路径优化的搜索空间,在保证求解精度的同时降低了求解规模,加快了求解速度。
附图说明
图1为海上风浪涌散点图。
图2为考虑风浪涌相关性下抽样的时间窗口描述。
图3为海上风电机组布局图。
图4为风电机组风速-输出功率拟合曲线。
图5为模型求解的整体思路。
图6为风、浪、涌预测期望值。
图7为不同策略的维护路径决策图,其中,图(7a)为策略1的维护路径决策图,图(7b)为策略2的维护路径决策图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
本发明提出了一种考虑尾流效应的海上风电机组维护路径随机规划方法,该模型考虑风、浪、涌的相关性,同时综合考虑各成本以及维护相关约束,海洋天气环境对船只航行时间,计及维护状态对尾流分布的影响,构建考虑尾流效应的维护路径随机规划模型。具体的建模步骤如下:
步骤1:如图1所示,风、浪、涌的相关性可以用Copula函数予以描述。设风、浪、涌随机变量x1、x2、x3的联合分布函数为E,边缘分布为F1(x1),F2(x2),F3(x3),那么存在一个Copula函数C,使得
E(x1,x2,x3)=C[F1(x1),F2(x2),F3(x3)]
考虑到各时刻下风、浪、涌的不确定性,运用多场景模型予以描述。假设各时段下风、浪、涌的预测值服从正态分布,运用基于Copula相关性的拉丁超立方抽样(latinhypercube sampling,LHS)进行大量采样,然后利用场景缩减法在抽样样本的基础上做缩减运算,得到各场景下的风、浪、涌数据。由于风、浪、涌任一值较大时均不可出海维护。若不考虑三者相关性下进行抽样,所得场景中可及的操作窗口将较为分散。如图2所示。
步骤2:对海上风电机组维护安排的总成本由维护人员成本cq,船只成本cb,惩罚成本cp和停机损失cs四部分构成,则目标函数C表示如下:
minC=cq+cb+cp+cs
步骤3:针对调度周期内需要维护的风电机组,引入检修机组索引变量i、j∈{1,2,…,n},那么标记机组所在节点i为运维船放下维护人员的节点,相应地接取维护人员的节点记为(n+i)。同时在一天的维护作业中开始和结束的节点分别记为0和2n+1,均指代港口。由此得放下节点的集合为Z-={1,2,3,...,n},接取节点的集合为Z+={n+1,n+2,...,2n},并且令放下节点与接取节点的集合Z=Z-∪Z+。各成本函数表达如下:
1)维护人员成本cq
Figure BDA0002371621820000081
2)船只成本cb
船只成本包括船只航行成本和驻留成本。
Figure BDA0002371621820000082
其中,船只航行时间Tω,l,i,j,d需要综合考虑风浪涌等海洋天气予以建模:
Figure BDA0002371621820000083
Figure BDA0002371621820000084
Figure BDA0002371621820000091
3)惩罚成本cp
当机组迟于其推荐时期维护时,会产生惩罚成本。
Figure BDA0002371621820000092
4)停机损失cs
停机损失包括机组维护带来的停机损失和潜在停机损失。
Figure BDA0002371621820000093
步骤4:构建考虑尾流效应的机组停机损失模型。
风电场布局如图3所示,为便于计算各机组的输出功率,引入机组索引变量r、u∈{1,2,…,s}对所有机组依次编号。
步骤5:建立风电机组输出功率与输入风速的三次函数模型,为将其线性化,如图4所示,采用分段函数予以拟合,那么输出功率可被表示为:
Figure BDA0002371621820000094
Figure BDA0002371621820000095
Figure BDA0002371621820000096
Figure BDA0002371621820000097
步骤6:考虑尾流效应的输入风速可表达为:
Figure BDA0002371621820000098
Figure BDA0002371621820000099
Iω,r,t=Nr,iIω,i,t
步骤7:采用MISOCP和松弛法将考虑尾流效应的输入风速模型转化为MILP问题,可线性化为:
Figure BDA0002371621820000101
Figure BDA0002371621820000102
Figure BDA0002371621820000103
Figure BDA0002371621820000104
Figure BDA0002371621820000105
Figure BDA0002371621820000109
步骤8:机组停机损失式中风电场的损失电量可表达为:
Figure BDA0002371621820000106
步骤9:在对模型求解时,引入每日可维护机组数目预评估策略,为运维船的路径规划确定有效的搜索方向,减少计算维度。
运维船l在第d天最大可维护的机组数ηl,d,可以根据各场景下的天气情况,出海窗口,运维船的参数等限制条件进行预评估,即:
Figure BDA0002371621820000107
步骤10:由于运维船需要到达放下节点和接取节点各一次,与维护路径决策变量之间的关系如下所示:
Figure BDA0002371621820000108
步骤11:对各场景下各船只每日最大可维护机组数ηl,d进行预评估,完成所有场景ηl,d的预评估后,由于本文制定的维护路径决策需要对所有场景均能适用,故选取各ηl,d的最小值代入ηl,d与维护路径决策变量之间的关系式中,以对运维船每日的维护路径产生约束,基于此完成模型求解。
本发明首先对复杂的海洋环境,运用场景法描述了风、浪、涌相关性和随机性。在此基础上,建立了维护成本最小化的目标函数,同时考虑到海洋天气对船只航行时间、出海可及窗口的影响构建了相关的维护约束。此外,考虑到维护状态与机组间尾流分布的关联性,构建了将尾流模型、维护状态与各机组输入风速的耦合函数,进而精细地量化了停机损失。针对所建立的考虑尾流效应的非线性模型,运用二阶锥规划结合松弛法转化为MILP模型。为了提高求解效率,在求解过程中,引入各船只每日最大可维护的机组数目预评估策略,确定了维护路径优化的有效搜索空间,降低了计算的规模。对某海上风电场机组的维护路径规划表明本发明所提方法的可行性和有效性。
具体应用场景1:对如图3所示布局的海上风电场进行算例仿真,共有8行8列共64台风电机组。机组的额定功率为3MW,叶轮半径63m,机组行间距、列间距均为1km,且风电机组的切入风速为3m/s,额定风速为13m/s,切出风速为25m/s。运维基地有两艘规格不同的运维船只,参数如表1所示,基地可用的每种类型技术人员数量均为6。依据状态监测系统,海上风电场需要维护的风电机组有12台,维护优化计划的时间范围为3天,且以15分钟为一个时间段,由于不建议夜晚出海维护,以7a.m.-7p.m.定义为白天,每日风、浪、涌预测的期望值如图6所示,开始时间为7a.m.,预测误差不超过10%。另外,维护风电机组所需的备件充足,人员和备件从运维船到风电机组的转移时间为15分钟,且不同运维船每天的最大可工作时间和海上环境天气有关。
表1运维船参数
运维船 载重量(kg) 人员容量(人)
船1 2000 12
船2 15000 12
结合上述基本数据,分别对如下3类情景进行建模分析:情景1:以风、浪、涌各期望值建立场景进行分析;情景2:对风、浪、涌不考虑Copula相关性进行抽样,产生多场景集进行分析;情景3:对风、浪、涌基于Copula相关性抽样,产生多场景集进行分析。通过本发明所提方法进行仿真计算,情景1和情景3下可以得到维护路径的优化方法,而情景2无法得到一种可以适应所有场景的维护路径,这是由于情景2下不考虑风、浪、涌相关性抽样,使得产生的某些场景出海可及窗口较为分散,且不同场景出海可及窗口分布的差异较大,求解时虽能各场景分别可以得到可行解,但却无法获得能满足所有场景约束下通用的可行路线,这将无法适应天气的不确定性带来的影响,故需要基于三者相关性进行抽样。
对情景1和情景3下所得的维护策略结果对比如表2所示,维护路径结果如表3所示,表3中维护路径下的数字代表检修机组索引下的机组编号,0指代港口,一天中的维护总是从港口出发,最后返回港口结束,维护周期为3天。由表3可以看出,情景1和情景3下所得维护路径结果不同,情景3由于考虑风、浪、涌不确定性产生的多场景,出海可及窗口会变化,为满足各场景,需要船2维护更多机组以充分利用时间窗。结合表2可知,情景3下的维护总成本要高于情景1,但将情景1下所得的各运维船维护路径决策代入情景3中是不可行的。它无法满足情景3下各场景的约束,而情景3所得的维护决策在情景1中仍然适用。因此,考虑风、浪、涌相关性及不确定性所安排的维护策略可以更好地应对不确定性的海洋天气环境,为决策者提供更加切实可行的方案。
表2不同情景下结果对比
Figure BDA0002371621820000121
表3运维船的维护路径结果
Figure BDA0002371621820000122
具体应用场景2:利用本发明所提模型对海上风电场中10台机组进行仿真,由于风电场中尾流效应的存在,下风向机组的输入风速会受到影响,机组输出功率发生变化,在考虑和不考虑尾流时机组停机损失的计算会有差异。为分析这种差异下尾流效应对路径决策制定的影响,以一天中对机组的维护为例进行研究。由于第一日预测的风、浪、涌较小,由上述可得,第一日维护的机组数较多,更适合对路径决策进行分析,故选取上述情景3下第一日并仅以船1维护4台机组(按检修机组索引下编号为4,7,9,10)为例,分别对以下两种策略进行分析:策略1,制定路径决策时不考虑机组间的尾流效应;策略2,制定路径决策时考虑机组间的尾流效应。两种策略下获得的运维船的路径决策结果如图7和表4所示。图7对待维护机组的地理位置及其交通路径特征进行抽象,忽略不同风机位置及航运通道的几何特征。图7中实线为运维船放下技术人员的路径,虚线为运维船的接取路径。对于策略1,当天整体的维护次序从港口出发到机组10,再从路径A依次到G以完成当天的维护和人员接取,最后返回港口,策略2同理,其从港口出发先到达机组4,然后沿着维护路径完成当天的维护任务。由于只研究4台机组当天的维护,故表4中惩罚成本为0,且其中的船只成本主要指代船只航行成本。另外,策略1下的总成本要高于策略2,两者的差异主要是在停机损失上,这是由于策略1不考虑机组间的尾流效应,处于下风向的机组停机损失也会较高。
表4两种策略下的各成本
Figure BDA0002371621820000131
由图7可得,在策略2下机组10的维护次序靠后,而在策略1下维护次序靠前。这是由于在考虑尾流效应下,处于下风向的机组输出功率低,潜在停机损失小,尽管机组10维护所需时间较长,但维护造成的停机损失会较小。在综合考虑航行路线和人员海上等待成本后,将机组10的维护次序排后能得到最佳的维护路径。而不考虑尾流效应时,各机组的潜在停机损失相差不大,此时综合考虑各机组维护成本后,将维护时间较长的机组排前有利于停机损失的下降,并且策略1的决策会尽快接取已经完成维护作业的一组人员以减小停机损失和人员海上等待成本,但这样航行路线对路径决策的影响被弱化,使得船只航行成本会更大。
因此,两种策略下的维护路径决策是不同的,策略1更侧重于机组会造成较大的停机损失制定维护决策,但这不符合海上风电场中下风向后排机组的实际运行情况,侧重于此作维护路径并不合理,故在路径决策时需要考虑机组间的尾流效应以使得运维船的路径安排更为经济可行。

Claims (9)

1.一种考虑尾流效应的海上风电机组维护路径随机规划方法,其特征在于,包括以下步骤:
1)描述海洋天气中风、浪、涌的相关性和不确定性;
2)考虑人员、船只、惩罚和停机损失的各项成本构成,以维护成本最小作为目标函数,并基于风、浪、涌对船只航行时间以及出海可及窗口的影响构建对应的约束条件;
3)由于维护状态会改变机组间尾流分布,结合尾流模型与维护状态描述各机组输入风速的变化,并精细化目标函数中的机组停机损失;
4)对目标函数进行求解,得到最优的维护路径。
2.根据权利要求1所述的一种考虑尾流效应的海上风电机组维护路径随机规划方法,其特征在于,所述的步骤1)中,采用Copula函数描述风、浪、涌的相关性,并且采用多场景模型描述各时刻下风、浪、涌的不确定性,则有:
E(x1,x2,x3)=C[F1(x1),F2(x2),F3(x3)]
其中,E(x1,x2,x3)为风、浪、涌随机变量x1、x2、x3的联合分布函数,C[·]为Copula函数,F1(x1),F2(x2),F3(x3)分别为风、浪、涌随机变量x1、x2、x3的边缘分布。
3.根据权利要求1所述的一种考虑尾流效应的海上风电机组维护路径随机规划方法,其特征在于,所述的步骤2)中,目标函数C的表达式为:
minC=cq+cb+cp+cs
Figure FDA0002371621810000011
Figure FDA0002371621810000012
Figure FDA0002371621810000013
Figure FDA0002371621810000014
其中,cq为维护人员成本,cb为船只成本,cp为惩罚成本,cs为停机损失,Ω为场景总数,ω为场景编号,D为总天数,d为天数编号,L为运维船只数,l为运维船编号,H为技术人员类型总数,h为技术人员类型编号,πω为场景ω的概率,Ql,d,h为运维船l在第d天离开节点0时船上类型为h的技术人员数量,
Figure FDA0002371621810000021
为场景ω下类型为h的单个技术人员每天所需费用,μi,h为维护节点i所需要类型为h的技术人员数量,
Figure FDA0002371621810000022
为场景ω下单个维护人员单位时间的海上等待成本,
Figure FDA0002371621810000023
Figure FDA0002371621810000024
分别为场景ω下运维船l在第d天到达放下的节点i和接取的节点n+i的时间段,τ为人员和备件从运维船转移到节点机组所在处所需时间段,Ti m为维护节点i所需时间段,Z-={1,2,3,...,n}为放下节点的集合,Z+={n+1,n+2,...,2n}为接取节点的集合,则Z=Z-∪Z+,并将一天内维护作业开始和结束的节点分别记为节点0和节点2n+1,即港口;
N为所有待维修机组所在节点的集合,其包括港口且i、j∈N,
Figure FDA0002371621810000025
为场景ω下运维船l在第d天海上航行的单位时间成本,Tω,l,i,j,d为场景ω下运维船l在第d天从节点i航行到节点j所需时间,xl,i,j,d为路径决策变量,即当船l在第d天从节点i驶向节点j的路径被选取时,其值为1,否则为0,ξj表示是否需要船只驻留,若节点j维护时不需要船只驻留,则ξj为1,否则为0,
Figure FDA0002371621810000026
为场景ω下船l在第d天在节点j维护时所需的驻留成本;
yi为机组所在节点i超过推荐的维护时间后延迟维护的天数,
Figure FDA0002371621810000027
为场景ω下机组所在节点i延迟维护后每天的惩罚成本;
t为时段编号,λi为机组所在节点i待维护部件故障停机的概率,cf为海上风电的上网电价,Pω,i,t为场景ω下机组所在节点i在时段t下正常运行时的输出功率,Δt为时段t的时长,
Figure FDA0002371621810000028
为场景ω下由于机组停机造成的损失电量。
4.根据权利要求3所述的一种考虑尾流效应的海上风电机组维护路径随机规划方法,其特征在于,所述的场景ω下运维船l在第d天从节点i航行到节点j所需时间Tω,l,i,j,d根据考虑风浪涌海洋天气建模,具体为:
Figure FDA0002371621810000029
Figure FDA0002371621810000031
Figure FDA0002371621810000032
其中,
Figure FDA0002371621810000033
为自然运行条件下运维船l从节点i到节点j的航行时间,aω,l,i,j,d为场景ω下船l在第d天从节点i到节点j适宜航行的概率,γ为风浪涌下的航行时间综合影响系数,
Figure FDA0002371621810000034
为场景ω下船l在第d天从节点i到节点j航行所需各时段下的海上平均风速,
Figure FDA0002371621810000035
为场景ω下时段t的海上风速。
5.根据权利要求3所述的一种考虑尾流效应的海上风电机组维护路径随机规划方法,其特征在于,所述的步骤2)中,约束条件具体包括:
维护任务约束:
在维护调度周期内,待维护机组都要得到维护,并确保每台待维护的机组只能放下和接取运维人员一次,且每艘运维船每天仅离开和返回港口一次,则有:
Figure FDA0002371621810000036
Figure FDA0002371621810000037
Figure FDA0002371621810000038
节点交通流守恒约束:
Figure FDA0002371621810000039
Figure FDA00023716218100000310
船只驻留约束:
Figure FDA00023716218100000311
其中,Zl为执行维护作业时需要运维船存在的机组所在节点的集合;
接取时间约束:
Figure FDA0002371621810000041
航行时间兼容性约束:
Figure FDA0002371621810000042
备件约束:
Figure FDA0002371621810000043
Figure FDA0002371621810000044
其中,Rj为维护机组j所需的备件重量,
Figure FDA0002371621810000045
为运维船l的负载容量,
Figure FDA0002371621810000046
为运维船l在第d天离开节点i后运维船上的备件重量;
人员约束:
Figure FDA0002371621810000047
Figure FDA0002371621810000048
Figure FDA0002371621810000049
Figure FDA00023716218100000410
其中,ql,i,d,h、ql,j,d,h为分别运维船l在第d天离开节点i、节点j后船上类型为h的技术人员数量,ql,0,d,h为运维船l在第d天离开港口节点0后船上类型为h的技术人员数量,
Figure FDA00023716218100000411
为运维船l的人员容量,φd,h为第d天可用的类型为h的技术人员数量,μj,h为节点j所需要类型为h的技术人员数量;
操作窗口约束:
Figure FDA00023716218100000412
其中,
Figure FDA00023716218100000413
为场景ω下运维船l在第d天的最大可工作时间段;
延迟维护约束:
Figure FDA00023716218100000414
Figure FDA0002371621810000051
式中:
Figure FDA0002371621810000052
为推荐维护机组j的最迟天数,yj为机组j超过推荐的维护时间后延迟维护的天数。
6.根据权利要求3所述的一种考虑尾流效应的海上风电机组维护路径随机规划方法,其特征在于,所述的步骤3)中,精细化后的机组停机损失中的损失电量
Figure FDA0002371621810000053
表示为:
Figure FDA0002371621810000054
其中,NT为总时段数,
Figure FDA0002371621810000055
为在风电场无维护任务时场景ω下在时段t机组所在节点r的输出功率,Iω,r,t为场景ω下在时段t机组所在节点r的维护状态,取值为1时表示维护停机,取值为0时表示正常运行,S为所有的机组所在节点总数。
7.根据权利要求3所述的一种考虑尾流效应的海上风电机组维护路径随机规划方法,其特征在于,所述的步骤4)具体包括以下步骤:
对各场景下各船只每日最大可维护机组数ηl,d进行预评估,在完成所有场景下ηl,d的预评估后,选取各ηl,d中的最小值代入每日最大可维护机组数与维护路径决策变量之间的关系式,对运维船每日的维护路径产生约束,并完成模型求解。
8.根据权利要求7所述的一种考虑尾流效应的海上风电机组维护路径随机规划方法,其特征在于,各场景下各船只每日最大可维护机组数ηl,d的预评估具体为:
Figure FDA0002371621810000056
其中,
Figure FDA0002371621810000057
为考虑航行时间和机组维护时间并在计及风、浪、涌的出海可及窗口下的最大可维护机组数;
Figure FDA0002371621810000058
为满足船只容量下的最大可维护机组数;
Figure FDA0002371621810000059
为满足备件和载重量下的最大可维护机组数。
9.根据权利要求7所述的一种考虑尾流效应的海上风电机组维护路径随机规划方法,其特征在于,每日最大可维护机组数与维护路径决策变量之间的关系式具体为:
Figure FDA00023716218100000510
CN202010052324.5A 2020-01-17 2020-01-17 一种考虑尾流效应的海上风电机组维护路径随机规划方法 Active CN111310972B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010052324.5A CN111310972B (zh) 2020-01-17 2020-01-17 一种考虑尾流效应的海上风电机组维护路径随机规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010052324.5A CN111310972B (zh) 2020-01-17 2020-01-17 一种考虑尾流效应的海上风电机组维护路径随机规划方法

Publications (2)

Publication Number Publication Date
CN111310972A true CN111310972A (zh) 2020-06-19
CN111310972B CN111310972B (zh) 2022-06-03

Family

ID=71161475

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010052324.5A Active CN111310972B (zh) 2020-01-17 2020-01-17 一种考虑尾流效应的海上风电机组维护路径随机规划方法

Country Status (1)

Country Link
CN (1) CN111310972B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112101897A (zh) * 2020-08-25 2020-12-18 华能国际电力股份有限公司浙江清洁能源分公司 海上风电智慧管理方法、装置及系统
CN112347694A (zh) * 2020-10-26 2021-02-09 重庆大学 含洋流-海上风电-潮汐流发电的海岛微网电源规划方法
CN113469379A (zh) * 2021-06-24 2021-10-01 中国能源建设集团广东省电力设计研究院有限公司 一种基于大数据中心的海上风电场运维管理方法及装置
CN113642937A (zh) * 2021-10-13 2021-11-12 深圳市信润富联数字科技有限公司 风机集群的运维排程方法、装置、电子设备及存储介质
CN114048930A (zh) * 2022-01-13 2022-02-15 广东电网有限责任公司揭阳供电局 一种超短期风电功率预测方法和装置
CN114157026A (zh) * 2021-11-19 2022-03-08 广东精铟海洋工程股份有限公司 一种海上风电运维的调度方法
CN114997644A (zh) * 2022-05-31 2022-09-02 中国能源建设集团广东省电力设计研究院有限公司 海上风电场运维调度方法、系统、计算机设备和存储介质
CN115271684A (zh) * 2022-09-26 2022-11-01 北京千尧新能源科技开发有限公司 一种海上风电场智慧运维管理方法及系统
CN117726060A (zh) * 2024-02-18 2024-03-19 华北电力大学 一种基于Vincenty公式的海上风机运维船路径优化方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150302332A1 (en) * 2013-05-03 2015-10-22 Gansu Electric Power Corporation Wind Power Technology Center Maintenance schedule optimization method for electric power system including large-scale wind power
CN105354632A (zh) * 2015-10-26 2016-02-24 江苏省电力公司电力经济技术研究院 一种考虑尾流效应的风电场功率优化分配策略
CN110175684A (zh) * 2019-05-30 2019-08-27 上海电力学院 一种海上风电场运维方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150302332A1 (en) * 2013-05-03 2015-10-22 Gansu Electric Power Corporation Wind Power Technology Center Maintenance schedule optimization method for electric power system including large-scale wind power
CN105354632A (zh) * 2015-10-26 2016-02-24 江苏省电力公司电力经济技术研究院 一种考虑尾流效应的风电场功率优化分配策略
CN110175684A (zh) * 2019-05-30 2019-08-27 上海电力学院 一种海上风电场运维方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨洪明 等: "考虑多风电场出力相关性的电力系统随机优化调度", 《电力自动化设备》 *
符杨 等: "考虑天气因素的海上风电机组预防性机会维护策略优化方法", 《中国电机工程学报》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112101897B (zh) * 2020-08-25 2022-03-22 华能(浙江)能源开发有限公司清洁能源分公司 海上风电智慧管理方法、装置及系统
CN112101897A (zh) * 2020-08-25 2020-12-18 华能国际电力股份有限公司浙江清洁能源分公司 海上风电智慧管理方法、装置及系统
CN112347694B (zh) * 2020-10-26 2023-05-12 重庆大学 含洋流-海上风电-潮汐流发电的海岛微网电源规划方法
CN112347694A (zh) * 2020-10-26 2021-02-09 重庆大学 含洋流-海上风电-潮汐流发电的海岛微网电源规划方法
CN113469379A (zh) * 2021-06-24 2021-10-01 中国能源建设集团广东省电力设计研究院有限公司 一种基于大数据中心的海上风电场运维管理方法及装置
CN113642937A (zh) * 2021-10-13 2021-11-12 深圳市信润富联数字科技有限公司 风机集群的运维排程方法、装置、电子设备及存储介质
CN114157026A (zh) * 2021-11-19 2022-03-08 广东精铟海洋工程股份有限公司 一种海上风电运维的调度方法
CN114157026B (zh) * 2021-11-19 2024-05-28 广东精铟海洋工程股份有限公司 一种海上风电运维的调度方法
CN114048930A (zh) * 2022-01-13 2022-02-15 广东电网有限责任公司揭阳供电局 一种超短期风电功率预测方法和装置
CN114997644A (zh) * 2022-05-31 2022-09-02 中国能源建设集团广东省电力设计研究院有限公司 海上风电场运维调度方法、系统、计算机设备和存储介质
CN115271684A (zh) * 2022-09-26 2022-11-01 北京千尧新能源科技开发有限公司 一种海上风电场智慧运维管理方法及系统
CN117726060A (zh) * 2024-02-18 2024-03-19 华北电力大学 一种基于Vincenty公式的海上风机运维船路径优化方法
CN117726060B (zh) * 2024-02-18 2024-05-24 华北电力大学 一种基于Vincenty公式的海上风机运维船路径优化方法

Also Published As

Publication number Publication date
CN111310972B (zh) 2022-06-03

Similar Documents

Publication Publication Date Title
CN111310972B (zh) 一种考虑尾流效应的海上风电机组维护路径随机规划方法
Deveci et al. A study on offshore wind farm siting criteria using a novel interval-valued fuzzy-rough based Delphi method
Bórawski et al. Development of wind energy market in the European Union
Latif et al. Comparative performance evaluation of WCA‐optimised non‐integer controller employed with WPG–DSPG–PHEV based isolated two‐area interconnected microgrid system
Scheu et al. Maintenance strategies for large offshore wind farms
Scholz-Reiter et al. A MILP for installation scheduling of offshore wind farms
CN110245795A (zh) 一种用于海上风电运维决策仿真模拟的方法及系统
CN111401664A (zh) 一种综合能源系统鲁棒优化调度方法及装置
CN107785935A (zh) 计及风险及机组组合的电网风电消纳能力评估方法
CN104037776A (zh) 随机惯性因子粒子群优化算法的电网无功容量配置方法
Pambudi et al. Wind turbine site selection in Indonesia, based on a hierarchical dual data envelopment analysis model
Zhou et al. Research on offshore wind power system based on Internet of Things technology
CN114301089A (zh) 一种风光联合发电系统储能容量配置优化方法
Leng et al. A novel evaluation method of renewable energy development based on improved rough set theory and grey cloud model
Moness et al. A critical review of research trends for wind energy in Egypt: recent progress and promising opportunities
Anna et al. Community smart grid utilizing dynamic demand response and tidal power for grid stabilization
CN116611628A (zh) 计及动态性能的综合能源生产单元两阶段容量配置方法
Anis et al. Challenges facing the implementation of wind energy in Algeria: a review paper
Lu et al. Annual Benefit Analysis of Integrating the Seasonal Hydrogen Storage into the Renewable Power Grids
Fuchs et al. Assessment of Offshore Wind Energy Opportunities and Challenges in the US Gulf of Mexico
Jäger-Roschko et al. Modelling of spare parts storage strategies for offshore wind
Cao et al. A Year-ahead Prediction Method of Monthly Wind Power Generation
Alkinani et al. Wind Energy in Saudi Arabia Opportunities and Challenges
Zhuang et al. Robust Optimization-Driven Strategy for Dynamic Planning of Power Transmission Networks
Adefisan OPTIMAL SITING AND CONNECTION OF WIND ENERGY FARMS IN VALENCIA COMMUNITY.

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