基于潮流、可靠性线性化分析的海上风电场集电系统集中式拓扑优化.pdf
《基于潮流、可靠性线性化分析的海上风电场集电系统集中式拓扑优化.pdf》由会员分享,可在线阅读,更多相关《基于潮流、可靠性线性化分析的海上风电场集电系统集中式拓扑优化.pdf(9页珍藏版)》请在咨信网上搜索。
1、第 43 卷 第 8 期2023 年 8 月电 力 自 动 化 设 备Electric Power Automation EquipmentVol.43 No.8Aug.2023基于潮流、可靠性线性化分析的海上风电场集电系统集中式拓扑优化王邦彦1,王秀丽1,帅轩越1,王一飞1,王帅2,辛超山3(1.西安交通大学 电气工程学院(陕西省智能电网重点实验室),陕西 西安 710049;2.国网经济技术研究院有限公司,北京 102209;3.国网新疆电力有限公司经济技术研究院,新疆 乌鲁木齐 830047)摘要:海上风电是新型电力系统的重要支撑,其集电系统需要大量昂贵海缆,占总成本比重较高,优化规划较
2、为困难。为此,提出适宜集电系统的含网损线性潮流模型,基于决策变量构建潮流约束。考虑开关的优化配置,将可靠性综合成本引入模型,并对场景进行解耦,基于网络流模型构建故障功率流变量实现具有可靠性的数学约束化建模分析。综合多方面成本建立集电系统规划模型,将风机聚类后,进行集中式优化得到最优的拓扑设计。同时,提出标准成本的概念以综合评判设计方案与求解模型的优劣。最后,通过该方法与现有方法的算例对比,说明该方法得到的方案具有更低的标准成本,该模型具有更好的计算精度。综上,所提出的模型及方法更适宜于集电系统的规划分析问题。关键词:海上风电场;集电系统;线性潮流;拓扑优化;可靠性分析中图分类号:TM614;T
3、M71 文献标志码:ADOI:10.16081/j.epae.2023020100 引言构建新型电力系统与2030碳达峰、2060碳中和两项任务将会促进可再生能源持续高速发展。风能作为可再生能源的重要组成部分,已经得到了大规模的开发。2021年年底我国并网风电容量达到了3.3108 kW,较2020年增长了16.6%1。由于海上具有更加充沛的未开发风能,海上风电是更具前景的发展方向。2021年全球海上风电并网装机容量增加了21.1 GW,我国的装机增量占据了全球的80%,连续4年居世界第一位2。随着沿海大省纷纷将海上风电发展写入“十四五”规划,海上风电将在未来持续蓬勃发展。海上风电场一般具有2
4、025 a的项目周期,建设投资巨大,需要合理规划风电场的建设方案以达到更好的经济性3。集电系统承担了汇集风机发出的电能的任务,将电能通过大量千米级海底电缆汇流至中压母线,随后由送出系统进行传输并网。集电系统的建设成本可达到风电场总成本的 15%30%4。除投资成本外,故障损失、网损等电量损失成本也应考虑在集电系统的全寿命周期成本之中5。由于集电系统拓扑复杂多变,且拓扑和可靠性问题、潮流问题存在紧密的耦合,综合考虑上述成本的规划问题十分复杂。实际工程中常采用链形拓扑进行简单设计,但是这样的设计方案的经济性必然不能达到最优。寻找最佳的拓扑设计方案为一个全局优化问题。集电系统拓扑规划中的可行解空间极
5、大,所以部分研究采用智能算法进行求解,如遗传算法6。但是智能算法的使用不能保证解的精度。此外,有大量研究采用聚类算法7进行风机聚类分组,削减搜索空间,提高求解效率。大部分研究应用了 CPLEX、Gurobi等商业求解器进行求解8。混合整数线性规划(mixed integer linear programming,MILP)模型是求解此问题的基础模型912,部分研究构建了混合整数二次规划(mixed integer quadratic programming,MIQP)模型以更加精确地计算网损13。由于潮流约束的非线性,采用完整的交流潮流模型将使整体模型变为一个混合整数非线性规划模型,求解较为困
6、难14。除此之外,许多研究提出采用迭代的方法求解拓扑规划问题1112。在可靠性方面,海上风电场不同于陆上风电场,运行工况恶劣,受大风大浪天气的影响,可靠性问题尤为突出1516,风电机组的实时数据采集与可靠性监测较为困难17。电缆故障导致的电量损失在运行成本中占据了一定比例。配置断路器以隔离故障是提高海上风电场的可靠性,减少故障损失的重要方法18。有研究指出部分开关配置方式相较于传统开关配置方式和完全开关配置方式具有更高的经济性19。更进一步地,若希望系统整体经济性达到最优,则应建立开关优化配置模型进行分析20。若在规划模型的目标函数中引入可靠性相关的成本,则可以避免反复校核调整,直接得出经济性
7、最优的方案21。值得注意的是,目前针对集电系统优化规划已有大量的研究,但仍存在以下不足之处:已有研收稿日期:20220814;修回日期:20221106在线出版日期:20230209基金项目:国家电网公司科技项目(5100202056458A0000)Project supported by the Science and Technology Project of SGCC(5100202056458A0000)第 8 期王邦彦,等:基于潮流、可靠性线性化分析的海上风电场集电系统集中式拓扑优化究建立的潮流模型较为粗略,以仅考虑有功的无网损线性潮流(即传输模型8)为主;尚未有研究将开关配置作为
8、变量引入集电系统拓扑优化问题,忽略了此方面的成本优化空间;尚未有研究评估潮流模型、开关优化对优化模型计算精度、寻优能力的影响。鉴于此,本文开展了相关针对性研究,主要贡献及创新点如下:基于提出的含网损线性潮流模型构建了潮流约束,提升了计算精度;提出故障功率流模型,通过数学约束计算集电系统可靠性;提出了标准成本的概念以分析各模型的优劣。1 含网损线性潮流建模及分析海上风电场集电系统大多呈放射状结构,无论是拓扑还是电压等级都与配电网非常相似。假设集电系统仅具有一个汇流站且拓扑为放射状,则可利用配电网的特性进行类比分析。由于海上风电场集电系统具有确定的最大运行方式,工程设计中电缆往往采取最小化截面的选
9、型。电缆电阻引起的电量损失成本在20余年的项目周期中叠加,占据总成本的部分比例,因而不能忽视。为此,应当合理选择所有电缆的截面,使得系统投资成本与网损成本达到综合最低。而对于网损成本的分析需进行所有运行方式下的潮流计算。为了构成集中式优化架构,需对潮流模型进行线性化。本文根据集电系统的特点,基于线性潮流模型构建了面向优化问题的考虑拓扑不确定的潮流约束。1.1复数量的近似文献 22 提出了一种计算电压的线性潮流模型,将配电网中节点电压表示为节点注入复功率的线性函数:U=U0(E+ZS*/U20)(1)式中:U为不包括平衡节点的(n-1)1维电压向量,n为配电网节点总数;U0为平衡节点电压幅值;E
10、为元素全是1的向量;Z为由部分导纳阵取逆得到的(n-1)(n-1)维阻抗矩阵;S为(n-1)1维节点注入复功率向量;上标“*”表示取共轭。考虑集电系统的特点,该关系可进行简化和重新推导。首先,集电系统中仅有风机作为电源,且风机发出的功率几乎全为有功功率。由于电抗产生的无功功率相对风机发出的有功功率极小,所以整个系统中流动的无功功率相较有功功率可以近似忽略。据此,忽略无功功率、电压相角、电抗得到实数方程,并将式(1)全部取标幺值得到:U=E+RP(2)式中:U为电压幅值向量;R为电阻矩阵;P为节点注入有功功率向量。若无特别说明,则后文中所有电压、电流、电阻、功率都默认为标幺值。值得注意的是,式(
11、1)实际上等同于配电网潮流的前推回代法第一次迭代的计算结果,存在计算误差,误差来源是注入功率和注入电流有如下关系(对于某一节点):I=S*(3)式中:I为某节点的注入电流(其为复数);S为某节点的注入复功率。更精确的关系应考虑节点的电压不等于额定电压:I=S*/U*(4)式中:U为某节点的电压(其为复数)。类比式(4)推导式(2),则某节点注入电流、电压的表达式分别为:I=P/U(5)U=E+RI(6)式中:I为某节点的注入电流幅值;P为某节点的注入有功功率;U为某节点的电压幅值;I为节点注入电流幅值向量。1.2基于泰勒展开的近似式(5)、(6)所示的约束使潮流计算除复数量近似外,不引入系统误
12、差,但I、U都为变量,其引入了分母为变量的分式,不为线性约束。考虑集电系统内节点电压都在11.05 p.u.范围内,在这样小的变化范围内可以对分式进行线性近似。在U=1 p.u.处将分式做泰勒展开23得:I=P(2-U)(7)利用该近似关系继续构建潮流约束。首先,风电场中的风速及风向不定,共同构成了多个场景,计算风机出力P在每个场景下的值并引入模型,最终形成多场景的随机规划。由于尾流效应对风机出力的影响是非常显著的24,为保证整体规划的精度,这一因素必须考虑在模型中。利用式(7)构建场景 s下节点i的注入电流表达式为:Ii,s=Pwti,s(2-Ui,s)iNW,sS(8)式中:变量Ii,s、
13、系数Pwti,s和变量Ui,s分别为场景s下节点i的注入电流、节点i处风机的出力和节点i的电压幅值;NW为风机节点集合;S为场景集合。1.3潮流约束构建及网损计算集电系统规划模型中,网架结构不确定,式(6)中的矩阵R、I都是由变量构成的矩阵,不能直接引入模型。同样地,如果构建关联矩阵描述注入电流与支路电流间的关系,则也会形成变量相乘的复杂情况。因此,利用连续非负变量Ii,j,s描述场景s下节点i向节点 j 通过节点间的支路传输的电流,则根据基尔霍夫电流定律有:jNIi,j,s=Pwti,s(2-Ui,s)+aNWIa,i,s iNW,sS(9)式中:N为风机节点集合与汇流站节点集合的并集。设y
14、i,j,t表示描述节点i、j 之间是否存在型号为t的电缆的 0-1 变量,yi,j,t=1 表示存在型号为 t 的电缆,电 力 自 动 化 设 备第 43 卷yi,j,t=0表示不存在型号为t的电缆。根据式(6)采用大M法建立考虑拓扑不确定的电压、电流关系如下:Uj,s-M(1-yi,j,t)Ui,s-Ii,j,tdi,jRtUj,s+M(1-yi,j,t)(i,j)B,tT,sS(10)式中:M为一个极大的数;di,j为节点i、j 间的距离;Rt表示型号为t的电缆的单位长度电阻值;T为电缆型号集合;集合B的定义如式(11)所示。B=(i,j)|i,jN(11)对于式(10),若两节点间不存在
15、型号为 t的电缆,则约束被松弛;若存在,则变为一个等式约束。潮流问题总会存在一个平衡节点,在集电系统中,潮流的平衡节点应为可以提供充足无功功率支撑的汇流站节点。同时,为保证风机的安全运行,且各节点电压应在11.05 p.u.范围内。用U1,s(下标1表示汇流站节点,后同)表示汇流站节点电压,增加约束条件如下:U1,s=1 p.u.sS(12)1 p.u.Ui,s1.05 p.u.iNW,sS(13)随后构建其他类型的支路约束。若节点i、j 间不存在支路,则不应存在支路电流。采用大M法构建如下约束:Ii,j,sMtTyi,j,t (i,j)B,sS(14)另外,每条电缆的运行电流都存在固定的上限
16、值,构建电缆运行电流约束条件如下:Ii,j,stTItyi,j,t (i,j)B,sS(15)式中:It表示型号为t的电缆的载流量。需要注意的是,式(15)实际上包含了式(14)所示的约束。至此,考虑拓扑不确定的新型潮流约束已经建立完毕。需要指出的是,由于该潮流模型没有忽略网损,实际计算网损时可以将汇集功率与所有风机输出功率作差以计算网损。2 集电系统可靠性数学约束化建模分析海底电缆故障导致的系统停运电量损失也是总成本中的重要组成部分。配置更多的断路器或减少风机串的电缆长度都可实现故障损失的下降。为此,本节考虑拓扑变量及开关配置变量建立了以数学约束为具体形式的可靠性模型,基于故障功率流构建了可
17、靠性综合成本的计算表达式,并将可靠性因素反映到集中优化模型中。该方法是可靠性理论及评估逻辑的数学约束形式化,以优化模型的形式完成较复杂的集电系统可靠性评估。2.1基于分块及网络流模型的可靠性建模集电系统的故障损失计算需要对集电系统进行可靠性评估。模拟法对于优化问题是一个黑箱模型,且耗时较长,所以可采用的可靠性分析方法只能为解析法。电缆的故障概率极小,对于项目全寿命周期内网损成本的影响极小;且电缆故障时损失功率远远大于网损功率。因此可采用故障场景与风速风向场景(后文简称潮流场景)解耦的方式进行分析。若考虑故障场景与潮流场景的叠加耦合,则将大幅增加计算量和分析难度。因此,本节提出故障场景与潮流场景
18、解耦的分析方法,即稳态潮流场景不考虑电缆故障,电缆故障场景不考虑网损及风机出力波动。由于故障场景的时间占比不超过5%,且可靠性计算中线路网损不超过 0.4%,场景解耦对于整体计算精度的影响可以忽略。更重要的是,该方法减少了潮流场景及故障场景的数量,大幅减少了模型中的变量、约束的数量。断路器可以隔断故障,减小故障损失,但是断路器的配置同样需要一定的成本。为了实现系统总成本的最优,本文令开关配置情况必须满足传统开关配置要求并且可继续灵活新增,从而建立开关配置的优化模型。在配电网中,断路器应配置在潮流上游,即线路更靠近母线的一侧(简称末端)。在集电系统中由于潮流方向和配电网相反,断路器应起到保护风机
19、正常并网的功能,同样应配置在线路更靠近母线的一侧。此时断路器恰好在潮流下游。由于风机出口存在断路器,所以线路两端并不需要完全配置2台断路器。故障的功率损失期望可以基于电缆流过的平均功率进行描述。利用网络流模型建立用于可靠性分析的潮流约束如下:jNPi,j=Pwti+aNWPa,i iNW(16)Pi,j Mt Tyi,j,t (i,j)B(17)式中:Pi,j为无故障时节点i、j 间传输的平均功率;Pwti为节点i处风机的平均出力。基于配电网可靠性评估理论,集电系统可仅考虑一阶事件,电缆的故障频率可以用故障率近似描述。放射状拓扑中功率的流通是单向的,根据此性质,每个节点汇集功率后仅通过一条支路
20、输出。设bi为描述节点i的出线末端是否存在断路器的0-1变量,bi=1表示存在断路器,bi=0表示不存在断路器。集电系统分块示例如图1所示,图中WT表示风机。根据断路器将集电系统分块25,计算每块停运事件的频率,再乘以该块停运导致的总功率损失(即分块图1集电系统分块示例Fig.1Example of collector system sectioning第 8 期王邦彦,等:基于潮流、可靠性线性化分析的海上风电场集电系统集中式拓扑优化出口处流出的功率)及故障时间,可得到集电系统故障总电量损失,这是最直观的方法。类比式(16)的网络流模型,结合大M法,可构建模型如下:j Ni,j=0j Nt T
21、yi,j,tdi,j+a NWa,i i NW(18)i,j-Mbii,ji,j+Mbi (i,j)A(19)i,jM(1-bi)(i,j)A(20)i,j,i,j,Pi,jMtTyi,j,t (i,j)B(21)Pi,j-M(1-bi)Pi,jPi,j+M(1-bi)(i,j)A(22)Pi,jMbi (i,j)A(23)tTyi,1,tbi iNW(24)W=TMTTRiNW()jNPi,jjNi,j(25)式中:A=(i,j)|iNW,jN;i,j为节点i、j间的故障流,表示故障率在分块中的汇集效应;0为电缆单位长度故障率,其为常数;i,j为节点i、j 间的传导故障流,若电缆无断路器,则
22、其值等于故障流,若电缆具备断路器,则其值为0;Pi,j为节点i、j 间的虚拟功率流,仅在分块出口处有取值,表示出口处流出的功率;TMTTR为电缆平均故障修复时间;W为年故障损失电量期望。式(24)确保系统满足最基本开关配置要求。故障流、传导故障流与虚拟功率流都是网络流变量。这些变量及相关的约束实现了可靠性分析中的逻辑推理及计算过程,使得本文提出的集中式拓扑及开关优化成为可能。故障流、传导故障流的单位都为次a。传导故障流描述了断路器隔离故障的逻辑,即故障不会影响到断路器处的潮流下游。该变量使得故障率顺着潮流方向汇集叠加得到故障流,且遇到断路器会停止汇集。因此,整个分块的故障率之和即为分块出口处电
23、缆的故障流。虚拟功率流描述了分块故障导致的停运功率,即分块出口处电缆流过的平均功率,单位为MW。由于故障损失仅需计算分块出口处的功率与故障流之积,为了区分出口处线路和普通线路,可令普通线路的虚拟功率流为0。该变量实现了对于线路是否处于分块出口处的逻辑判断。上述分析建立了故障流模型和虚拟功率流模型,对可靠性问题及开关优化问题进行了详细的评估,但是最终功率损失表达式仍是变量相乘的形式。如果将该项引入目标函数,则构成二次规划模型0.5xTQx+cTx(x为模型决策变量,c为二次规划模型中含决策变量一次项的系数矩阵),其系数矩阵Q不定,固定或松弛所有整数变量后下层问题仍非凸,无法求解。因此,需要对上述
24、方法进行进一步改进。2.2基于故障功率流的线性化方法上述建模过程说明了不能采用故障流与虚拟功率流相乘的方法,需要采用另外的建模方式,将原模型故障流变量与虚拟功率流相乘的方式改为线性表达式,并只能利用一个变量表示故障所带来的功率损失。虚拟功率只能利用变量表示,不能修改,但是分块出口处的故障流可以表示为分块内各电缆故障率之和。可使虚拟功率流变量修改为与分块出口处的功率等值,并将分块内每段电缆的故障率与虚拟功率流相乘,得到一个新的流变量,利用此流变量表示功率损失在馈线中的流通和叠加,因而不妨将此流变量定义为故障功率流。修改后的虚拟功率流在分块内具有相同的值,需要根据节点间相连情况和断路器情况构建虚拟
25、功率流的等式约束。使用大M法构建所有流变量的网络流约束,最后得到故障功率流模型如下:jNPi,j-M(1-bi)PijNPi,j+M(1-bi)iNW(26)Pj-M(1-tTyi,j,t+bi)PiPj+M(1-tTyi,j,t+bi)(i,j)A(27)Pi0di,j-M()1-tTyi,j,taNPwi,a-fNWPwf,iPi0di,j+M()1-tTyi,j,t (i,j)A(28)Pwi,jMtTyi,j,t (i,j)B(29)jNPi,j=Pwti+aNWPa,i iNWPi,jMtTyi,j,t (i,j)BtTyi,1,tbi iNW(30)式中:Pi为修改后的节点i出线的
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 潮流 可靠性 线性化 分析 海上 电场 系统 集中 拓扑 优化
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。