CN114928057B - 一种适用于新能源多端柔直***的多目标潮流优化方法 - Google Patents
一种适用于新能源多端柔直***的多目标潮流优化方法 Download PDFInfo
- Publication number
- CN114928057B CN114928057B CN202210651904.5A CN202210651904A CN114928057B CN 114928057 B CN114928057 B CN 114928057B CN 202210651904 A CN202210651904 A CN 202210651904A CN 114928057 B CN114928057 B CN 114928057B
- Authority
- CN
- China
- Prior art keywords
- direct current
- power
- new energy
- converter station
- power flow
- 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.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims abstract description 56
- 230000003068 static effect Effects 0.000 claims abstract description 8
- 230000005540 biological transmission Effects 0.000 claims description 22
- 238000011217 control strategy Methods 0.000 claims description 13
- 239000002699 waste material Substances 0.000 claims description 12
- 230000002068 genetic effect Effects 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 4
- 238000005286 illumination Methods 0.000 claims description 3
- 238000007665 sagging Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种适用于新能源多端柔直汇集***的多目标潮流优化方法,属于电力技术领域。当新能源经MMC‑MTDC***并网时,容易引起MTDC***换流站的潮流变化,且汇集***结构主要为环网,而此拓扑会削弱MMC换流站站级控制对潮流的控制能力,很难实现并网***的最优运行。该方法以直流电压偏差最小、弃风、弃光量最小及静态电压稳定裕度最大作为优化目标,搭建新能源多端柔性直流汇集***多目标潮流优化模型,并通过内点法与NSGA‑III算法结合的混合交替算法对多目标优化模型进行求解。实现风、光资源的优化调度,迭代次数少且收敛性好,能够保证对多目标有效优化。同时最大程度提高新能源多端柔性直流汇集***的稳态运行能力。
Description
技术领域
本发明涉及电力技术领域,具体涉及一种适用于新能源多端柔直***的多目标潮流优化方法。
背景技术
发展可再生能源已成为许多国家积极应对气候变化的核心内容和推进能源转型的重要战略途径。
柔性直流输电技术(Voltage-Sourced Converter Based HVDC,VSC-HVDC)可控性强,功率潮流翻转迅速,不受距离限制,无需交流侧提供换相电流。在远距离传输时可实现送端与受端的解耦,减小新能源并网对输电***的影响。因此,柔性直流输电技术成为新能源接入电网的理想传输方式,其中基于MMC的柔性直流输电技术应用更加广泛。
目前,新能源的利用和开发受到了世界各国的重视,尤其是光伏发电、风能发电等在近年来的能源供给中所占比例不断增加。然而,风电、光伏等可再生能源都具有间歇性、波动性、随机性和不可控性等特点,将其用于新能源发电属于间歇式电源。近年来,我国逐步将各种大规模可再生能源接入电网,可再生能源的高比例快速渗透,导致需要调整电力***规划方案,需要更好的预测方法以及潮流优化方法。再加上各种分布式能源采集存储装置、“即插即用”设备以及大型城市负荷需求的不断增加,传统的电网结构、运行技术和电力装备等难以满足我国大规模电力需求,因此需要采用新的电网结构、先进技术以及相应辅助装备来满足未来能源格局的深刻变化。
发明内容
本发明的目的是提供一种适用于新能源多端柔性直流汇集***的多目标潮流优化方法。
本发明是一种适用于新能源多端柔直***的多目标潮流优化方法,共有4个换流站,其中采用下垂控制的换流站有2个,其余2个换流站采用定有功功率控制,并在***中加装直流潮流控制器;其步骤为:
步骤(1):***处于正常运行时,确定采用下垂控制以及定有功功率控制的换流站i的节点电压直流电压Ui(i=1,2,3,4),确定各直流支路电阻R12、R13、R24、R34;
步骤(2):***处于正常运行时,计算各换流站支路电流I12、I13、I24、I34;
步骤(3):根据步骤(2)求解得到各支路传输有功功率P12、P13、P24、P34以及各换流站的输入有功功率Pi(i=1,2,3,4);
步骤(4):针对新能源出力不确定性与相关性进行风电场与光伏电站不确定性建模;
步骤(5):针对新能源多端柔性直流汇集***进行协调控制策略建模;
步骤(6):建立包含直流潮流控制器的新能源多端柔性直流汇集***多目标潮流优化模型:
步骤(6-1):所建立数学模型中包含三个目标函数,分别为直流电压偏差最小、弃风、弃光量以及静态电压稳定裕度最大;
步骤(6-2):确定模型中的等式约束条件,分别包括步骤(4)及步骤(5)中的新能源并网约束和MMC-MTDC***协调控制约束,并且确定***直流方程;
步骤(6-3):确定模型中的不等式约束条件,分别包括S3中所提的换流站传输功率约束以及***运行约束,同时确定***交流方程;
步骤(7):利用内点法与遗传算法对多目标潮流优化模型进行求解,对整个新能源多端柔性直流汇集***进行多目标优化。
本发明的有益效果为:
(1)本发明在传统的多目标优化方法中引入了新能源出力不确定性建模以及MMC-MTDC***协调控制策略建模,综合考虑了新能源的出力和换流站的协同控制在潮流优化模型中的约束条件,可以实现风、光资源的优化调度,迭代次数少且收敛性好,能够保证对多目标有效优化。同时所搭建模型能够根据***不同的运行目标给出新能源基地出力的协调运行方案,提高了***运行的稳定性与可靠性。
(2)采用内点法与NSGA-III遗传算法相结合的混合算法进行多目标潮流优化模型的求解,合理解决了优化模型在求解中的离散部分与连续部分的求解不当问题,能够精确处理模型中的优化问题,使得模型的迭代次数减小且容易收敛,能够实现鲁棒运行。
附图说明
图1是本发明一个实施例的包含直流潮流控制器的四端新能源柔性直流汇集***图,图 2是本发明一个实施例的内点法与NSGA-III遗传算法结合的混合算法流程图,图3是本发明一个实施例的NSGA-III遗传算法流程图,图4是本发明一个实施例的多目标潮流优化结果图。
具体实施方式
本发明是一种适用于新能源多端柔直***的多目标潮流优化方法,共有4个换流站,其中采用下垂控制的换流站有2个,其余2个换流站采用定有功功率控制,并在***中加装直流潮流控制器;其步骤为:
步骤(1):***处于正常运行时,确定采用下垂控制以及定有功功率控制的换流站i的节点电压直流电压Ui(i=1,2,3,4),确定各直流支路电阻R12、R13、R24、R34;
步骤(2):***处于正常运行时,计算各换流站支路电流I12、I13、I24、I34;
步骤(3):根据步骤(2)求解得到各支路传输有功功率P12、P13、P24、P34以及各换流站的输入有功功率Pi(i=1,2,3,4);
步骤(4):针对新能源出力不确定性与相关性进行风电场与光伏电站不确定性建模;
步骤(5):针对新能源多端柔性直流汇集***进行协调控制策略建模;
步骤(6):建立包含直流潮流控制器的新能源多端柔性直流汇集***多目标潮流优化模型:
步骤(6-1):所建立数学模型中包含三个目标函数,分别为直流电压偏差最小、弃风、弃光量以及静态电压稳定裕度最大;
步骤(6-2):确定模型中的等式约束条件,分别包括步骤(4)及步骤(5)中的新能源并网约束和MMC-MTDC***协调控制约束,并且确定***直流方程;
步骤(6-3):确定模型中的不等式约束条件,分别包括S3中所提的换流站传输功率约束以及***运行约束,同时确定***交流方程;
步骤(7):利用内点法与遗传算法对多目标潮流优化模型进行求解,对整个新能源多端柔性直流汇集***进行多目标优化。
以上所述的方法,在所述步骤(2)中,***稳态运行时根据基尔霍夫环路定律可得换流站支路电流I12和I34,由I12同理可得支路电流I13以及I24;
其中,R12、R34为支路电阻;Ui(i=1,2,3,4)为换流站i的节点电压;c1为直流潮流控制器串入直流线路中的等效变压器变比。
以上所述的方法,在所述步骤(3)中,换流站1与换流站2支路传输有功功率P12为:
由式(2)同理得剩余支路有功功率P13、P24、P34;
进一步得到含有DCPFC的四端MMC-MTDC的换流站1的输入功率P1为:
由式(3)同理得到其余换流站的输入功率P2、P3、P4。
以上所述的方法,在所述步骤(4)中,首先通过风电出力的随机模糊建模可得风电出力模型:
其中,Pdcn、PWP是风电机组的额定功率与输出功率;vN、vin、vout分别为风机的额定风速、切入风速及切出风速;
然后进行光伏出力的随机建模同样可得光伏出力模型:
其中,PPV,N、PPV分别为光伏额定功率和输出功率;rN为额定光照强度。
以上所述的方法,在所述步骤(5)中,进行MMC-MTDC***协调控制策略建模,将换流站电压调制比Zi设置为优化变量;给定约束范围为:
同时设置下垂系数约束范围为:
换流站协调控制策略给出了由式(7)唯一决定,因此给定以下协调控制策略约束条件:
其中,为电网有功功率期望值;/>与/>分别为直流电压参考值最小与最大值。
以上所述的方法,在所述步骤(6)中,建立包括目标函数、等式约束条件以及不等式约束条件的多目标潮流优化模型;
所设定的目标函数为:
(1)直流电压偏差最小
其中,ΔU为直流电压偏移量;为直流节点i的期望直流电压值;/>和/>分别为直流节点i直流电压最小与最大值;ndc为汇集***直流节点集合;
(2)弃风、弃光量最小
其中,nN表示***的送端节点集合;ni为***连接节点i的个数;Href,i(s)为第s个***的出力预测值;Hwp为直流节点i的弃风弃光量;Hx,i与为并网实际值与预测值;
(3)静态电压稳定裕度最大
maxVSM=ξmin (11)
其中,ξmin为收敛潮流的雅克比矩阵的最小奇异值;
所设定的等式约束条件:
(1)MMC-MTDC协调控制约束,包括式(5)-式(7);
(2)新能源并网***出力约束,包括式(10);
(3)含直流潮流控制器的新能源多端柔性直流汇集***功率平衡方程为:
其中,Pl1,i、Pl2,i、Ll3,i分别为***的直流损耗、MMC运行损耗及线路的交流损耗;PDCPFC,i为直流潮流控制器的损耗;Ux为前述DCPFC的换流站节点电压;ndc为前述汇集***直流节点集合;Gij、Gi分别为直流与交流线路的电导;Ic,i为换流站的电流;Us,i∠θs,i、Uc,i分别为并网***输入的交流电压相量与换流站直流侧电压;ni(n=a,b,c)为换流站损耗的多项式系数;
所设定的不等式约束条件:
(1)汇集***运行约束:
式(14)表示直流电压约束、有功功率约束、相角差约束及直流线路电流约束;
(2)换流站i传输功率约束:
其中,Sdi为换流站i的传输功率;表示直流支路传输功率的上限值;
综合上述目标函数及约束条件,提出新能源直流汇集***多目标潮流优化模型可表示为:
其中,f(x,u)为目标函数;h(x,u)为等式约束条件;g(x,u)为不等式约束条件;x主要由直流控制量及状态变量组成,其核心变量包括[Us,Ud,δ,Ic,Idij];u为控制变量构成的集合,主要包括[Z,k,Pd,Qs,Qc,Px]。
以上所述的方法,在所述步骤(7)中,针对本发明所构建的含直流潮流控制器的新能源直流汇集***多目标潮流优化问题,拟采用快速非支配排序遗传算法与内点法相结合的方法交替求解。
下面结合附图和具体实施例对本发明的具体实施方式做进一步说明:
如图1所示,MTDC***为四端新能源MTDC***,共有4个换流站,其中采用下垂控制的换流站有2个,其余2个换流站采用定有功功率控制。***构成主要包括:风电场(WP1、WP2、WP3)、光伏电站(PV)、受端负荷中心、MMC换流站以及直流潮流控制器(m1、m2)。
如图2所示,本发明的控制原理为:采用快速非支配排序遗传算法与内点法相结合的方法交替求解。首先利用加权法将多目标转换为单一目标,然后将该目标用数学公式统一表示,由此多目标潮流的优化转变成对单目标的潮流优化。同时将该单目标潮流中离散变量视为连续变量。最后利用混合交替算法进行求解,具体先采用内点法对连续部分的***模型进行求解,再利用NSGA-III算法对离散部分进行求解。
具体实施步骤如下:
本发明的适用于新能源多端柔性直流汇集***的多目标潮流优化方法,共有4个换流站,其中采用下垂控制的换流站有2个,其余2个换流站采用定有功功率控制,且为改善并网***的潮流控制自由度不足及调整环网各支路潮流,在***中加装直流潮流控制器(DC Power Flow Controller,DCPFC)来实现潮流调控;
步骤(1):在MMC-MTDC***正常运行的情况下,根据基尔霍夫环路定律可得支路电流I12和I34,由I12同理可得支路电流I13以及I24;
其中,R12、R34为各支路电阻;Ui(i=1,2,3,4)为换流站i的节点电压;c1为直流潮流控制器串入直流线路中的等效变压器变比;
步骤(2):根据式(1)可得换流站1与换流站2支路传输有功功率P12为:
由式(2)同理可得剩余支路有功功率P13、P24、P34;
进一步得到含有DCPFC的四端MMC-MTDC***换流站1的输入功率P1为:
由式(3)同理可得到其余换流站的输入功率P2、P3、P4;
步骤(3):对于新能源出力不确定性建模,首先进行风电出力的随机模糊建模,将风速当作随机模糊变量,可得风电出力模型:
其中,Pdcn、PWP是风电机组的额定功率与输出功率;vN、vin、vout分别为风机的额定风速、切入风速及切出风速;
然后进行光伏出力的随机建模结合太阳能电池方阵总面积与总光电转换率可得到光伏出力模型:
其中,PPV,N、PPV分别为光伏额定功率和输出功率;rN为额定光照强度;
步骤(4):针对MMC-MTDC***进行协调控制策略建模,将换流站电压调制比Zi设置为优化变量。Zi为MMC换流站i的电压调制比,可以体现为对***的运行调节能力,并且可以根据运行需求改变Zi的值,因此给定约束范围为:
下垂系数决定了汇集***中的不平衡有功及直流电压分配,且其大小由MMC换流站有功功率及直流电压偏差决定,使得换流站的传输容量得到了充分利用。下垂系数由于MMC- MTDC***设定的***运行界限而设定很小的约束范围,很难实现MMC换流站的灵活调节,从而引起弃风、弃光现象。但约束范围亦不应设定过大,以防止由于***功率扰动导致的直流电压越限等运行工况。因此设置下垂系数约束范围为:
换流站协调控制策略给出了换流站直流电压、有功功率参考值及下垂系数等三个控制参数,但上述三个参数不能由式(7)唯一决定,基于此给定以下协调控制策略约束条件,将有功功率参考值设置为定值,而直流电压参考值设定为优化变量:
其中,为电网有功功率期望值;/>与/>分别为直流电压参考值最小与最大值;
步骤(5):新能源汇集***可简化为交直流混合***,包含多个交直流节点。其多目标优化旨在综合考虑新能源的优化调度与换流站之间的协调控制,提高新能源消纳水平,进而提升***的运行稳定性与灵活性。所搭建含直流潮流控制器的新能源多端柔性直流汇集***多目标潮流优化数学模型,主要包括:三个目标函数、等式约束条件以及不等式约束条件;
本发明所设定的目标函数为:
(1)直流电压偏差最小
针对新能源直流汇集***中的所有节点,为改善各MMC换流站直流节点的电能质量,将直流电压偏差最小作为优化目标,设置为:
其中,ΔU为直流电压偏移量;为直流节点i的期望直流电压值;/>和/>分别为直流节点i直流电压最小与最大值;ndc为汇集***直流节点集合;
(2)弃风、弃光量最小
汇集***为提高风、光资源消纳水平的优化目标设置为:
其中,nN表示新能源并网***的送端节点集合;ni为新能源并网***连接节点i的个数; Href,i(s)为第s个并网***的出力预测值;Hwp为直流节点i的弃风弃光量;Hx,i与为新能源并网实际值与预测值;
(3)静态电压稳定裕度最大
当雅克比(Jacobi)矩阵的某一奇异值接近零时,任意小的扰动会引起***状态变量的大幅度变化,当***接近稳定运行极限,即直流电压接近运行极限值,会造成***的不稳定。因此***是否稳定将由雅克比矩阵的最小奇异值(Minimum Singular Value,MSV)决定,同时雅克比矩阵的MSV的大小亦表明***稳定运行的电压相对临近程度。综合考虑上述因素,本发明将MSV作为静态电压稳定指标:
maxVSM=ξmin (11)
其中,ξmin为收敛潮流的雅克比(Jacobi)矩阵的最小奇异值;
本发明所设定的等式约束条件:
(1)MMC-MTDC协调控制约束,包括式(5)-式(7);同时电网侧换流站采用定有功功率控制,其功率设定为:
Px,i=Pl i∈nl (12)
其中,nl为负荷中心的节点;Pl为电网发出的调度值;
(2)新能源并网***出力约束,包括式(10);
(3)含直流潮流控制器的新能源多端柔性直流汇集***功率平衡方程为:
其中,Pl1,i、Pl2,i、Ll3,i分别为汇集***的直流损耗、MMC换流站运行损耗及MMC线路的交流损耗;PDCPFC,i为直流潮流控制器的损耗;Ux为前述DCPFC的换流站节点电压;ndc为前述汇集***直流节点集合;Gij、Gi分别为直流线路与交流线路的电导;Ic,i为换流站的电流;Us,i∠θs,i、Uc,i分别为新能源并网***输入的交流电压相量与换流站直流侧电压;ni(n=a,b,c) 为换流站损耗的多项式系数;
本发明所设定的不等式约束条件:
(1)汇集***运行约束:
式(15)表示***的直流电压约束、有功功率约束、相角差约束及直流线路电流约束;
(2)换流站i传输功率约束:
其中,Sdi为换流站i的传输功率;表示直流支路传输功率的上限值;
新能源多端柔性直流汇集***的稳定运行有其范围界限,需要***控制变量以及状态变量的约束限制。综合上述的目标函数及约束条件,提出新能源直流汇集***多目标潮流优化模型可表示为:
其中,f(x,u)为目标函数;h(x,u)为等式约束条件;g(x,u)为不等式约束条件;x主要由直流控制量及状态变量组成,其核心变量包括[Us,Ud,δ,Ic,Idij];u为控制变量构成的集合,主要包括[Z,k,Pd,Qs,Qc,Px]。潮流优化模型的优化变量主要包括换流站节点电压、换流站下垂系数、换流站传输功率等控制参数;
步骤(6):利用内点法与遗传算法对多目标潮流优化模型进行求解,对整个新能源多端柔性直流汇集***进行多目标优化;
(1)内点法求解
内点法鲁棒性强且容易收敛,但很难处理离散变量的优化问题。因此利用内点法与加权法相结合来求解连续部分的优化问题,通过对目标函数进行加权处理,设置相应的权重系数:
minφ1ΔU+φ2Hwp-φ3ξmin (18)
综合考虑权重系数对结果的影响,在第一次进行加权计算时,令Φ1=Φ2=Φ3=1/3,在后续的迭代计算中依据各个权重所表示的目标函数对其所对应的权重系数进行不断更新;
为保证离散变量不变,引入松弛离散约束,将最优潮流模型简化为一般非线性优化模型:
对不等式约束条件gmin与gmax分别引入非负松弛变量l和w,对于等式约束h(x)引入拉格朗日乘子为消除引入松弛变量的非负性约束,构造对数壁垒函数μ,则构成的拉格朗日函数为:
对上式中控制变量及状态变量利用偏导解变量求解,可得:
联立上式来求解此非线性规划问题,即利用内点法求解该连续变量的多目标潮流优化;
(2)NSGA-III遗传算法求解
现代最优潮流算法用于解决寻找全局最优解的问题,能够精确处理多目标优化问题中的离散变量。而NSGA-III算法作为使用较为广泛的潮流优化算法是在NSGA-Ⅱ算法基础上进行改进,主要是在选择操作上需要同时满足参考点与其欧氏距离,即需要求取每个种群个体在参考平面上的欧氏距离并进行比较。作为多目标优化问题的求解方法,NSGA-III算法在求解 3个或3个以上目标的优化问题方面得到了有效的利用,且其有效性已经得到验证;
因此,对前述多目标优化模型中的离散变量本文采用NSGA-III算法进行求解。在模型连续变量不发生改变的前提下,基本步骤如下且具体的算法流程图如图3所示:
1)初始化新能源多端柔性直流汇集***相关参数,设置离散变量优化目标及状态变量约束条件;
2)设定NSGA-III算法的初始参数,假设种群数量为N并且将其初始化,同时设置迭代次数t为0。利用种群数量N生成H个参考点,设置Im为最大迭代次数;
3)利用决策变量生成初始化种群Pt;
4)进行交叉、变异等操作进行子代更新,从父代种群Pt产生子代种群Qt;
5)合并两个种群得到新的种群Rt,同时进行潮流计算得出种群Rt中每个个体的适应度;
6)对种群Rt进行快速非支配排序,得到非支配层F1,F2,F3……;
7)从F1开始构造种群St,令St=St∪Fi(i=1,2,…),直到St≥N。如果St=N,则Pt+1=St并转到步骤(9);否则,继续以下步骤;
8)自适应归一化St中的目标值,将每一解与参考点相互关联维持多样性。从St中选择 K=|St|-N个个体,然后,令Pt+1=St;
9)增加迭代次数t,即t=t+1。若t=Im,则输出多目标潮流优化的解;否则重复步骤(4)~(9);
10)求得最优折衷解。
由图4可得,潮流优化之后,***的弃风弃光量明显减小,直流电压偏差减小,静态电压稳定裕度升高。因此,***目标函数得到优化,保证了***在安全稳定的范围内运行,且具有一定的经济性。
综上所述,新能源多端柔性直流汇集***多目标潮流优化模型可以实现风、光资源的优化调度,迭代次数少且收敛性好,能够保证对多目标有效优化。同时所搭建模型能够根据***不同的运行目标给出新能源基地出力的协调运行方案,提高了***运行的稳定性与可靠性。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (5)
1.一种适用于新能源多端柔直***的多目标潮流优化方法,其特征在于,设共有4个换流站,其中MMC1、MMC2换流站采用下垂控制,其余2个换流站采用定有功功率控制,同时,在MMC1与MMC2之间加装直流潮流控制器m1,在MMC1与MMC3之间加装直流潮流控制器m2,MMC4与MMC1分别同时连接MMC2和MMC3;
包括如下步骤:
步骤(1)、***处于正常运行时,确定采用下垂控制以及定有功功率控制的换流站i的节点电压直流电压Ui,i=1,2,3,4,确定各直流支路电阻R12、R13、R24、R34;
步骤(2)、***处于正常运行时,计算各换流站支路电流I12、I13、I24、I34;
步骤(3)、根据步骤(2)求解得到各支路传输有功功率P12、P13、P24、P34以及各换流站的输入有功功率Pi,i=1,2,3,4;
步骤(4)、针对新能源出力不确定性与相关性进行风电场与光伏电站不确定性建模;
步骤(5)、针对新能源多端柔直汇集***进行协调控制策略建模;
步骤(6)、建立包含直流潮流控制器的新能源多端柔性直流汇集***多目标潮流优化模型:
所述步骤(6)中,建立包括目标函数、等式约束条件以及不等式约束条件的多目标潮流优化模型;
所设定的目标函数为:
(1)直流电压偏差最小
其中,ΔU为直流电压偏移量;为直流节点i的期望直流电压值;/>和/>分别为直流节点i直流电压最小与最大值;ndc为汇集***直流节点集合;
(2)弃风、弃光量最小
其中,nN表示***的送端节点集合;ni为***连接节点i的个数;Href,i(s)为第s个***的出力预测值;Hwp为直流节点i的弃风弃光量;Hx,i与为并网实际值与预测值;
(3)静态电压稳定裕度最大
maxVSM=ξmin (3)
其中,ξmin为收敛潮流的雅克比(Jacobi)矩阵的最小奇异值;
所设定的等式约束条件:
(1)MMC-MTDC协调控制约束:
进行MMC-MTDC***协调控制策略建模,将换流站电压调制比Zi设置为优化变量,给定约束范围为:
同时设置下垂系数约束范围为:
换流站协调控制策略由式(5)唯一决定,因此给定以下协调控制策略约束条件:
其中,为电网有功功率期望值,Pref.i为电网有功功率参考值;/>与/>分别为直流电压参考值最小与最大值;
(2)新能源并网***出力约束,包括式(2);
(3)含直流潮流控制器的新能源多端柔性直流汇集***功率平衡方程为:
其中,Pl1,i、Pl2,i、Ll3,i分别为***的直流损耗、MMC运行损耗及线路的交流损耗;PDCPFC,i为直流潮流控制器的损耗;Ux为前述DCPFC的换流站节点电压;ndc为前述汇集***直流节点集合;Gij、Gi分别为直流与交流线路的电导;Ic,i为换流站的电流;Us,i∠θs,i、Uc,i分别为并网***输入的交流电压相量与换流站直流侧电压;ai、bi、ci为换流站损耗的多项式系数;
所设定的不等式约束条件:
(1)汇集***运行约束:
式(9)表示直流电压约束、有功功率约束、相角差约束及直流线路电流约束;
(2)换流站i传输功率约束:
其中,Sdi为换流站i的传输功率;表示直流支路传输功率的上限值;
综合上述目标函数及约束条件,提出新能源多端柔直汇集***多目标潮流优化模型可表示为:
其中,f(x,u)为目标函数;h(x,u)为等式约束条件;g(x,u)为不等式约束条件;x由直流控制量及状态变量组成,其核心变量包括[Us,Ud,δ,Idc,Id,ij];u为控制变量构成的集合,包括
[Z,k’,Pd,Qs,Qc,Px]。
2.根据权利要求1所述的适用于新能源多端柔直***的多目标潮流优化方法,其特征在于,在所述步骤(2)中,***稳态运行时根据基尔霍夫环路定律可得换流站支路电流I12和I34,由I12同理可得支路电流I13以及I24;
其中,R12、R34为支路电阻;Ui为换流站i的节点电压,i=1,2,3,4;c1为直流潮流控制器m1串入直流线路中的等效变压器变比。
3.根据权利要求2所述的适用于新能源多端柔直***的多目标潮流优化方法,其特征在于,在所述步骤(3)中,换流站1与换流站2支路传输有功功率P12为:
由式(13)同理可得剩余支路有功功率P13、P24、P34;
进一步得到含有DCPFC的四端MMC-MTDC的换流站1的输入功率P1为:
由式(14)同理可得到其余换流站的输入功率P2、P3、P4。
4.根据权利要求1所述的适用于新能源多端柔直***的多目标潮流优化方法,其特征在于,在所述步骤(4)中,首先通过风电出力的随机模糊建模可得风电出力模型:
其中,Pdcn、PWP是风电机组的额定功率与输出功率;vN、vin、vout分别为风机的额定风速、切入风速及切出风速;
然后进行光伏出力的随机建模同样可得光伏出力模型:
其中,PPV,N、PPV分别为光伏额定功率和输出功率;rN为额定光照强度。
5.根据权利要求1所述的适用于新能源多端柔直***的多目标潮流优化方法,其特征在于,在步骤(6)中,针对本发明所构建的含直流潮流控制器的新能源直流汇集***多目标潮流优化问题,采用快速非支配排序遗传算法与内点法相结合的方法交替求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210651904.5A CN114928057B (zh) | 2022-06-08 | 2022-06-08 | 一种适用于新能源多端柔直***的多目标潮流优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210651904.5A CN114928057B (zh) | 2022-06-08 | 2022-06-08 | 一种适用于新能源多端柔直***的多目标潮流优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114928057A CN114928057A (zh) | 2022-08-19 |
CN114928057B true CN114928057B (zh) | 2024-06-04 |
Family
ID=82812142
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210651904.5A Active CN114928057B (zh) | 2022-06-08 | 2022-06-08 | 一种适用于新能源多端柔直***的多目标潮流优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114928057B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115276019B (zh) * | 2022-09-22 | 2022-12-27 | 东南大学溧阳研究院 | 一种基于自适应下垂控制的潮流优化方法 |
CN115438522B (zh) * | 2022-11-08 | 2023-03-24 | 清华大学 | 城市轨道交通混合式直流牵引供电***最优潮流建模方法 |
CN117879051A (zh) * | 2023-12-04 | 2024-04-12 | 国家电网有限公司华东分部 | 一种受端电网受电能力评估方法、装置、介质及电子设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107546734A (zh) * | 2017-10-23 | 2018-01-05 | 国网天津市电力公司 | 一种电压型直流潮流控制器变比优化计算方法 |
CN109510231A (zh) * | 2019-01-07 | 2019-03-22 | 兰州理工大学 | 一种新型的直流电压自适应下垂控制策略 |
CN110912177A (zh) * | 2019-12-15 | 2020-03-24 | 兰州交通大学 | 一种多端柔性直流输电***多目标优化设计方法 |
CN111884195A (zh) * | 2020-07-27 | 2020-11-03 | 山东大学 | 考虑概率潮流的直流电网下垂控制裕度优化方法及*** |
-
2022
- 2022-06-08 CN CN202210651904.5A patent/CN114928057B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107546734A (zh) * | 2017-10-23 | 2018-01-05 | 国网天津市电力公司 | 一种电压型直流潮流控制器变比优化计算方法 |
CN109510231A (zh) * | 2019-01-07 | 2019-03-22 | 兰州理工大学 | 一种新型的直流电压自适应下垂控制策略 |
CN110912177A (zh) * | 2019-12-15 | 2020-03-24 | 兰州交通大学 | 一种多端柔性直流输电***多目标优化设计方法 |
CN111884195A (zh) * | 2020-07-27 | 2020-11-03 | 山东大学 | 考虑概率潮流的直流电网下垂控制裕度优化方法及*** |
Non-Patent Citations (1)
Title |
---|
基于改进电压下垂控制的MMC-MTDC ***协调控制策略;宋汶秦等;《电气应用》;20211231;第40卷(第10期);第70-77页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114928057A (zh) | 2022-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114928057B (zh) | 一种适用于新能源多端柔直***的多目标潮流优化方法 | |
CN111953013B (zh) | 真双极柔性直流输电***故障下的自适应优化调控方法 | |
CN109950907B (zh) | 含电力电子变压器的交直流混合配电网的调度方法及*** | |
CN108134401B (zh) | 交直流混合***多目标潮流优化及控制方法 | |
CN110265991B (zh) | 一种直流微电网的分布式协调控制方法 | |
CN108304972B (zh) | 一种基于供需互动和dg运行特性的主动配电网网架规划方法 | |
CN108808715B (zh) | 考虑直流网络故障功率的多端柔直***静态安全分析方法 | |
CN105809265A (zh) | 含分布式可再生能源配电网柔性互联装置的容量配置方法 | |
Arani et al. | A novel control method based on droop for cooperation of flywheel and battery energy storage systems in islanded microgrids | |
CN110912177A (zh) | 一种多端柔性直流输电***多目标优化设计方法 | |
Sahri et al. | Performance improvement of Hybrid System based DFIG-Wind/PV/Batteries connected to DC and AC grid by applying Intelligent Control | |
Kaysal et al. | Hierarchical energy management system with multiple operation modes for hybrid DC microgrid | |
CN112467748A (zh) | 三相不平衡主动配电网双时标分布式电压控制方法及*** | |
CN113644693A (zh) | 含可再生能源/氢能交直流混联***分散式运行管控方法 | |
CN116667325A (zh) | 一种基于改进布谷鸟算法的微电网并网运行优化调度方法 | |
CN116388153A (zh) | 一种有源配电网中柔性互联装备优化配置方法 | |
CN106786766B (zh) | 一种基于p2g技术的提高风电最大并网容量的方法 | |
CN115276019B (zh) | 一种基于自适应下垂控制的潮流优化方法 | |
CN109193634B (zh) | 基于多端柔性直流的海岛电网运行优化方法及*** | |
CN116707023A (zh) | 基于源荷相关性聚类的主动配电网分层分区综合优化方法 | |
CN114996908B (zh) | 考虑智能软开关接入的主动配电网扩展规划方法及*** | |
CN116031889A (zh) | 基于鲁棒深度强化学习的多智能体分布式电压控制方法 | |
CN113595135B (zh) | 一种并弱电网型电氢耦合直流微电网双层同步控制方法 | |
CN115358479A (zh) | 一种虚拟电厂协同电网规划方法及*** | |
CN112909982A (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 |