CN113224788B - 含有风电接入的交直流混联电网静态电压稳定评估方法 - Google Patents
含有风电接入的交直流混联电网静态电压稳定评估方法 Download PDFInfo
- Publication number
- CN113224788B CN113224788B CN202110468364.2A CN202110468364A CN113224788B CN 113224788 B CN113224788 B CN 113224788B CN 202110468364 A CN202110468364 A CN 202110468364A CN 113224788 B CN113224788 B CN 113224788B
- Authority
- CN
- China
- Prior art keywords
- node
- direct current
- convex
- svsm
- formula
- 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
- 230000003068 static effect Effects 0.000 title claims abstract description 22
- 238000011156 evaluation Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 54
- 238000005457 optimization Methods 0.000 claims description 63
- 239000010410 layer Substances 0.000 claims description 54
- 230000009977 dual effect Effects 0.000 claims description 30
- 230000005540 biological transmission Effects 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 21
- 238000005516 engineering process Methods 0.000 claims description 8
- 239000002356 single layer Substances 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 230000000295 complement effect Effects 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 150000001875 compounds Chemical class 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- CGFFKDRVHZIQHL-UHFFFAOYSA-N 1-but-3-en-2-yl-3-(methylcarbamothioylamino)thiourea Chemical compound CNC(=S)NNC(=S)NC(C)C=C CGFFKDRVHZIQHL-UHFFFAOYSA-N 0.000 claims description 2
- 239000003990 capacitor Substances 0.000 claims description 2
- 239000011248 coating agent Substances 0.000 claims description 2
- 238000000576 coating method Methods 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000002360 preparation method Methods 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 5
- 238000005315 distribution function Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 3
- 239000002355 dual-layer Substances 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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/02—Circuit arrangements for ac mains or ac distribution networks using a single network for simultaneous distribution of power at different frequencies; using a single network for simultaneous distribution of ac power and of dc power
-
- 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/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/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (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
本发明公开了一种含有风电接入的交直流混联电网静态电压稳定评估方法,所述方法包括:获取数据,所述数据包括电网当前运行状态数据和网架数据,以及风电场的出力范围上/下界数据;将所获取到的数据输入交直流混联电网SVSM区间计算的双层最优潮流模型中,计算获得SVSM区间的上/下界及对应的风电场出力值;将获得的SVSM区间的上/下界及对应的风电场出力值应用至电网运行调度。本方法仅需根据风电场出力不确定波动范围的上/下边界信息,不仅能够直接计算得到交直流混联电网的SVSM区间,还可以得到在风电场出力波动区间范围内SVSM区间取下界和上界对应的风电场出力值,为交直流混联电网的安全稳定运行提供更完整的决策信息,以保证电网的安全稳定运行。
Description
技术领域
本发明涉及电力技术领域,具体涉及一种含有风电接入的交直流混联电网静态电压稳定评估方法。
背景技术
随着经济和社会的不断发展,我国电网的负荷快速增长,从2009年开始我国电力负荷一直稳居世界第一。特别是经济较好的东部地区,电力负荷普遍较重,需要超/特高压大容量交直流输电线路从西部能源基地输入大量电力,因而目前南方电网、华东电网等区域电网都已形成了大型交直流混联电网格局。在交直流混联电网运行中,由于直流馈入线路不具备常规电源的动态无功支撑能力,在负荷快速增长条件下,受端电网的电压稳定问题日益突出。另外,为了应对能源危机和降低环境污染,新能源发电技术在世界各国得到大力发展。其中,风电作为应用最广泛的新能源发电技术之一,也得到了迅速发展。风电场的出力受到风速等自然因素的影响,具有较大的不确定波动特性,这使得风电场的大规模并网给交直流混联电网的安全运行带来了极大的挑战。其中,风电场的大规模接入会引起交直流混联电网的静态电压稳定水平产生不确定波动,给系统的安全稳定运行带来很大冲击。因此,考虑风电场出力不确定波动的交直流混联电网静态电压稳定评估问题,是一个亟待解决的关键技术问题。
静态电压稳定裕度(static voltage stability margin,SVSM)是反映电力系统静态电压稳定性的一个最常用的指标,其表示着系统当前状态能够承担的最大负荷增长量。
目前,计算电网SVSM的方法通常有连续潮流法和最优潮流法,不过这两种方法都只适用于确定性系统,对于含有风电接入的交直流混联电网,他们只能得到风电场出力保持不变时系统的SVSM静态电压稳定裕度值,不能完整地反映风电场出力不确定波动对系统SVSM的影响。如果根据确定性计算得到的SVSM值进行运行决策,当由于天气变化导致风电场出力出现大范围波动时,将会给系统的安全运行带来很大的风险。
发明内容
本发明的目的在于克服上述背景技术所存在的至少一种技术问题,提供一种含有风电接入的交直流混联电网静态电压稳定评估方法。
为实现上述目的,本发明的技术方案是:
含有风电接入的交直流混联电网静态电压稳定评估方法,所述方法包括:
获取数据,所述数据包括电网当前运行状态数据和网架数据,以及风电场的出力范围上/下界数据;
将所获取到的数据输入交直流混联电网SVSM区间计算的双层最优潮流模型中,计算获得SVSM区间的上/下界及对应的风电场出力值;
将获得的SVSM区间的上/下界及对应的风电场出力值应用至电网运行调度。
进一步地,所述交直流混联电网SVSM区间计算的双层最优潮流模型包括:
计算SVSM区间上界的min-min模型和计算区间下界的max-min模型两个双层优化模型。
进一步地,通过如下方式来计算获得SVSM区间的上/下界及对应的风电场出力值:
通过凸松弛技术将内层非凸的非线性优化模型转化为凸规划模型,并基于对偶优化理论得到内层凸优化模型的对偶优化模型,将计算SVSM区间下界的max-min模型转化为max-max模型;将min-min双层优化模型和max-max双层优化模型转化为单层优化模型求解。
本发明与现有技术相比,其有益效果在于:
本方法建立了考虑风电场出力不确定波动的交直流混联电网SVSM区间计算的双层最优潮流模型,包括计算SVSM区间上界的min-min模型和区间下界的max-min模型两个双层优化模型。通过凸松弛技术将内层非凸的非线性优化模型转化为凸规划模型,并基于对偶优化理论得到内层凸优化模型的对偶优化模型,从而将计算SVSM区间下界的max-min模型转化为max-max模型。min-min双层优化模型和max-max双层优化模型都可方便转化为单层优化模型求解。本方法仅需根据风电场出力不确定波动范围的上/下边界信息,不仅能够直接计算得到交直流混联电网的SVSM区间,还可以得到在风电场出力波动区间范围内SVSM区间取下界和上界对应的风电场出力值,为交直流混联电网的安全稳定运行提供更完整的决策信息,保证电网的安全稳定运行。
与蒙特卡洛法比较,本方法计算得到的SVSM区间具有较高的精度,并且计算时间明显减少,降低了对计算机的性能要求。
附图说明
图1为本发明实施例提供的含有风电接入的交直流混联电网静态电压稳定评估方法的流程图;
图2为39节点交直流系统的组成示意图。
具体实施方式
实施例:
下面结合附图和实施例对本发明的技术方案做进一步的说明。
参阅附图1所示,本实施例提供的含有风电接入的交直流混联电网静态电压稳定评估方法具体包括如下步骤:
101、获取数据,所述数据包括电网当前运行状态数据和网架数据,以及风电场的出力范围上/下界数据;
102、将所获取到的数据输入交直流混联电网SVSM区间计算的双层最优潮流模型中,计算获得SVSM区间的上/下界及对应的风电场出力值;
103、将获得的SVSM区间的上/下界及对应的风电场出力值应用至电网运行调度。
由此可见,本方法仅需根据风电场出力不确定波动范围的上/下边界信息,不仅能够直接计算得到交直流混联电网的SVSM区间,还可以得到在风电场出力波动区间范围内SVSM区间取下界和上界对应的风电场出力值,为交直流混联电网的安全稳定运行提供更完整的决策信息,以保证电网的安全稳定运行。
具体地,上述交直流混联电网SVSM区间计算的双层最优潮流模型包括:
计算SVSM区间上界的min-min模型和计算区间下界的max-min模型两个双层优化模型。
通过如下方式来计算获得SVSM区间的上/下界及对应的风电场出力值:
通过凸松弛技术将内层非凸的非线性优化模型转化为凸规划模型,并基于对偶优化理论得到内层凸优化模型的对偶优化模型,将计算SVSM区间下界的max-min模型转化为max-max模型;将min-min双层优化模型和max-max双层优化模型转化为单层优化模型求解。
目前已有的考虑风电场出力不确定性的静态电压稳定分析方法主要是蒙特卡络法和概率分析法。蒙特卡络法需要对先对风电场出力不确定变量进行随机抽样得到一系列样本,并对每个风电场出力样本逐一进行确定性的连续潮流或最优潮流计算得到对应SVSM值,进而统计得到SVSM的概率分布特性;但是,为了获得SVSM的准确的概率分布特性,蒙特卡络法需要抽样出大量的风电场出力样本来逐一进行计算,计算量非常大。概率分析法根据已知的风电场出力的概率分布函数,通过以半不变量法为代表的解析法或以点估计法为代表的近似法得到SVSM的概率分布函数;不过,概率分析法需要事先得到风电场出力的准确概率分布函数,而要得到随机变量准确的概率分布函数在实际运行中往往已需要大量统计工作,较难获得。
而本方法是基于SVSM计算的最优潮流法和区间优化理论,采用区间数描述风电场出力的不确定波动特性,建立了考虑风电场出力不确定波动的交直流混联电网SVSM区间计算的双层最优潮流模型,包括计算SVSM区间上界的min-min模型和区间下界的max-min模型两个双层优化模型。通过凸松弛技术将内层非凸的非线性优化模型转化为凸规划模型,并基于对偶优化理论得到内层凸优化模型的对偶优化模型,从而将计算SVSM区间下界的max-min模型转化为max-max模型。min-min双层优化模型和max-max双层优化模型都可方便转化为单层优化模型求解。本方法仅需根据风电场出力不确定波动范围的上/下边界信息,不仅能够直接计算得到交直流混联电网的SVSM区间,还可以得到在风电场出力波动区间范围内SVSM区间取下界和上界对应的风电场出力值。与蒙特卡洛法比较,本方法计算得到的SVSM区间具有较高的精度,并且计算时间明显减少,降低了对计算机的性能要求。
交直流混联电网SVSM区间计算的最优潮流模型的内层模型的目标函数如下:
min-λ (1)
式中,λ为负荷增长系数,SVSM计算问题实质上就是在满足电力系统各种运行约束的情况下,求取负荷增长系数λ的最大值。
约束条件如下:
(1)各节点的功率平衡约束
普通交流节点的功率平衡方程如下:
式中,PGi和QGi为节点i的发电机有功和无功出力;PWi和QWi为节点i的风电场有功和无功出力;PLi0和QLi0为节点i的初始有功和无功负荷;和为节点i负荷有功和无功的增长方式,一般可取Gij和Bij为节点i和j之间的互电导和互电纳;Vi和Vj分别为节点i和节点j的电压幅值,θij为节点i和节点j的电压相角差;NLCC和NVSC分别表示LCC和VSC直流换流站的交流节点集合;n为系统节点数。
连接LCC直流换流站的交流节点的功率平衡方程如下:
式中,Kpi为节点i连接直流输电系统的极对数;Udi为节点i连接直流换流站的直流侧电压;Idi为节点i连接直流换流站的直流电流;为节点i连接直流换流站的功率因数角。式(4)中的正负号对于整流站为负,逆变站为正。
连接VSC直流换流站的交流节点的功率平衡方程如下:
式中,Pvi和Qvi为节点i连接VSC直流换流站从交流系统吸收的有功和无功。
(2)直流输电线路的运行特性约束
对于LCC直流输电线路,考虑换流变压器和换相电抗的影响,则运行特性约束如下:
UdRi=UdIi+IdiRdc (8c)
式中,Kdi为节点i连接直流换流站的换流变压器变比;θi为节点i连接直流换流站的换流器控制角;Xci为节点i连接直流换流站的换相电抗;Rdc为直流线路电阻;UdRi为整流侧直流电压;UdIi为逆变侧直流电压。
在求解优化模型的过程中,给定LCC换流站的控制方式。对于LCC直流输电线路采用控制直流电流和整流站直流电压的方式,结合式(8c),可以得到逆变站直流电压。即在给定控制方式后,LCC直流电路的直流电压和直流电流都可认为是已知量,则LCC直流输电线路的运行特性约束为式(8a)和式(8b)。
对于VSC直流输电线路,忽略换流站等效电阻的影响,则运行特性约束如下:
式中,μi为直流电压利用率,在PWM调制策略确定时为定值,一般为0.866;Mi为调制比;Vi为节点i连接直流换流站的交流侧电压;Yi为节点i连接直流换流站的等效导纳;δi为节点i连接的交流节点电压与换流器输入电压的相角差;gil,d为消去联络节点后直流网络节点电导矩阵中的元素;m为连接VSC直流换流站的交流节点数。
对于VSC直流输电系统,采用整流站控制直流电压和无功功率,逆变站控制有功功率和无功功率的方式。结合式(9d),可以得到整流站直流电流、逆变站直流电流和直流电压,则VSC直流输电线路的运行特性约束为式(9a)~(9c)。
(3)发电机有功出力增长方式约束:
假定除平衡节点发电机外的其他发电机按照其有功出力裕度占总有功出力裕度的比例来分配负荷有功的总增长量,如式(10b)。对于有功出力已经到达上界的发电机,不再承担负荷有功增长量,因此其有功增长系数为0。
(4)发电机机端电压和无功出力的互补约束:
(QGi-QGimin)Vai=0 (11a)
(QGi-QGimax)Vbi=0 (11b)
VGi=VGiref+Vai-Vbi (11c)
式中,QGimax和QGimin分别为第i台发电机无功出力的上界和下界;VGi和VGiref分别为第i台发电机的机端电压以及其给定参考值;Vai和Vbi是辅助变量,用于修正发电机无功出力到达下界或者上界时的机端电压。
(5)变量的上下限约束:包括发电机节点有功和无功出力的上下限,交流节点电压及辅助变量的上下限,LCC直流换流站的控制角和功率因数角的上下限,及VSC直流换流站的调制比和相角差的上下限,即:
xmin≤x≤xmax (12)
当采用区间数描述风电场出力的不确定波动,相应的SVSM也将会是一个区间数。基于区间优化理论,可将计算SVSM区间的最小值与最大值分别描述为下面两个双层优化模型,如式(13)和(14):
式中,PW=(PW1,PW2,…,PWn)T;x为系统中除风电场出力外的其他变量组成的向量。
问题(13)称为乐观模型,用于求解SVSM区间的上界;问题(14)称为悲观模型,用于求解SVSM区间的下界。由于乐观模型式(13)的两层规划问题都是计算最小值的min问题,因此可以通过直接把两层规划合并为一层的方法进行求解。合并后的单层规划问题如下:
此问题是连续非线性规划问题,可以直接调用GAMS软件中的CONOPT求解器进行求解。
在悲观模型式(14)中,两层规划问题分别是max问题和min问题,不能直接合并为一层优化问题求解。为了解决这个问题,引入凸规划的对偶优化理论,先将内层min模型转化为其对偶优化max模型。对于如式(16)的非线性规划问题,若f(x)和a(x)都是连续可微的凸函数,则称优化问题式(16)为凸规划问题。并且,式(17)的非线性规划问题为原优化问题式(16)的对偶优化问题。
式中,v=(v1,v2,…,vr)T,vi为对偶变量,r为原优化问题(16)的约束条件的数目;
原问题与对偶问题之间存在着以下关系:如果x0是原问题的最优解,那么存在一组对偶变量v0,使得(x0,v0)是对偶问题的一个最优解,并且,原问题最优解x0所对应的目标函数值f(x0)与对偶问题最优解(x0,v0)所对应的目标函数值g(x0,v0)相同,即:
由此可见,只要把悲观模型式(14)中的内层优化模型转化成其对偶优化模型,那么悲观模型将会转化为max-max问题,此时便可以直接把两层规划问题合并为一层规划问题进行求解。
但是,悲观模型式(14)的约束条件中,各节点的功率平衡约束,直流输电线路的运行特性约束,发电机端电压约束和无功出力的互补约束都含有非线性项,如三角函数和二次项,并不符合式(16)中的凸规划模型的条件。因此,在转化为对偶优化问题之前,还需要对内层规划问题的约束条件进行凸松弛处理。
凸松弛过程如下:
(1)交流线路运行约束的凸松弛
Rij=Rji (20b)
Iij=-Iji (20c)
由于新增的约束(20a)是非凸的,需要通过二阶锥松弛将其转化为如下的凸不等式约束:
||2Rij;2Iij;Ui-Uj||2≤Ui+Uj (21)
为了使凸化前后的系统等价,需要补充约束条件来保证每个独立回路中各支路两端节点电压相角差之和为0。并将此约束转化为中间变量的线性等式如式(23)。
而由于节点电压幅值变量采用其平方的形式表示,发电机端电压和无功的互补方程相应地转化如下:
(QGi-QGimin)Uai=0 (25a)
(QGimax-QGi)Ubi=0 (25b)
而此时发电机无功互补方程的式(25)中存在两个连续变量相乘的项,是非凸的,令qai=QGiUai,qbi=QGiUbi,则利用凸包络法可将式(25)可转换为线性约束如下:
式中,Uaimax和Uaimin表示Uai的上限和下限,Ubimax和Ubimin表示Ubi的上限和下限。
(2)LCC直流线路运行约束的凸松弛
在SVSM计算过程中,当确定了直流输电系统的控制方式后,直流电压和直流电流都是已知量,因此连接LCC直流换流站的交流侧节点的功率平衡方程中的直流功率项是已知常数;同时,假定换流变变比不会随着负荷增长而变化,即Kdi为定值。需要对式(8a)和(8b)进行凸松弛。
式(8a)存在非线性项Vicosθi,将cosθi整体看成一个连续变量,则该项为两个连续变量相乘,令wi=Vicosθi,可利用凸包络法将其转化为如下线性约束:
式中,Vimax和Vimin表示Vi的上限和下限。
经过上述松弛后,式(8b)可转化为如下线性等式约束:
(3)VSC直流线路运行约束的凸松弛
||2Kvi;2Lvi;Oi-Ui||2≤Oi+Ui (32)
综上,对优化模型中的非线性部分进行凸松弛处理后,交直流混联电网SVSM区间计算的最优潮流模型转化为目标函数为线性、约束条件为线性或二阶锥的凸规划模型,如下:
通过将内层优化模型进行凸松弛处理后,可基于凸规划对偶理论获得内层min模型的对偶优化max模型,进而可将求SVSM区间下界的max-min双层优化模型转化为max-max双层优化模型,并进一步转化为式(34)的单层优化模型。
此问题是连续非线性规划问题,可以直接调用GAMS软件中的CONOPT求解器进行求解。
下面结合一个仿真试验来对本方法进行进一步的验证说明:
以修改的39节点交直流混联系统为例进行仿真计算。该系统是在IEEE 39节点系统基础上将交流线路14-4改为常规LCC直流输电线路、将11-6和11-12两个交流线路改为三端柔性直流输电系统,并在节点9接入一个风电场,如附图2所示。常规LCC直流输电线路电阻Rdc=0.02,换流站的换相电抗为Xc=0.1;整流侧换流变变比为1.0,逆变侧换流变变比为1.05;换流器控制角θi取值的上/下限为15°和8°;换流站功率因数取值的上/下限为0.7和0.4。三端VSC直流输电系统中的3个VSC直流线路电阻均为Rdc=0.01,换流站等效电纳为15;调制Mi取值的上/下限为1和0.8;相角差δi取值的上/下限为5°和0°。交流节点电压的上/下限为为1.1和0.8。对常规LCC直流输电线路采用定直流电流和整流站直流电压的控制方式。对三端VSC直流输电系统,采用整流站控制直流电压和无功功率,逆变站控制有功功率和无功功率的方式。结合直流线路运行约束,可推算出直流电路的直流量,如表1所示。
表1.直流量取值
给定不同的风电场有功出力波动区间,风电场的功率因数取1.0。分别采用本发明的双层最优潮流法和传统的蒙特卡罗法来计算SVSM区间。其中,蒙特卡罗法采样次数为1000次,在给定风电出力区间中随机抽样并对每个样本进行计算得到其对应的SVSM值,进而综合得到SVSM区间。得到结果比较如表2。可以看到,所提出的双层最优潮流法得到的SVSM区间与蒙特卡罗法得到的SVSM区间结果很接近,SVSM区间上界和下界的绝对误差都小于0.0023,很接近零,表明所提出的双层最优潮流法计算得到的SVSM区间具有很高的精度。通过对比两种方法的计算时间可知,所提出的双层最优潮流法的计算耗时比蒙特卡洛法大大减少,只需耗时1s左右,而蒙特卡洛法的耗时超过140s。因此,本发明的双层最优潮流法更适合应用于实际大型交直流混联电网的SVSM区间计算,也更具有实际应用价值。
表2.考虑风电场出力不确定波动的SVSM区间计算结果比较
综上,本发明与现有技术相比具有如下技术优势:
1)相比于确定性SVSM计算方法,该方法不仅能够直接计算得到风电场出力波动区间范围对应的交直流混联电网SVSM区间,还可以得到在风电场出力波动区间范围内SVSM区间取下界和上界对应的风电场出力值,为交直流混联电网的安全稳定运行提供更完整的决策信息;
2)相比于蒙特卡洛法,该方法计算得到的交直流混联电网SVSM区间具有较高的精度,并且计算时间大大减少。
上述实施例只是为了说明本发明的技术构思及特点,其目的是在于让本领域内的普通技术人员能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡是根据本发明内容的实质所做出的等效的变化或修饰,都应涵盖在本发明的保护范围内。
Claims (8)
1.含有风电接入的交直流混联电网静态电压稳定评估方法,其特征在于,所述方法包括:
获取数据,所述数据包括电网当前运行状态数据和网架数据,以及风电场的出力范围上/下界数据;
将所获取到的数据输入交直流混联电网SVSM区间计算的双层最优潮流模型中,计算获得SVSM区间的上/下界及对应的风电场出力值;
将获得的SVSM区间的上/下界及对应的风电场出力值应用至电网运行调度;
交直流混联电网SVSM区间计算的双层最优潮流模型的内层模型的目标函数如下:
min-λ (1)
式中,λ为负荷增长系数;
所述目标函数的约束条件包括:
(1)各节点的功率平衡约束
普通交流节点的功率平衡方程如下:
式中,PGi和QGi为节点i的发电机有功和无功出力;PWi和QWi为节点i的风电场有功和无功出力;PLi0和QLi0为节点i的初始有功和无功负荷;和为节点i负荷有功和无功的增长方式,取Gij和Bij为节点i和j之间的互电导和互电纳;Vi和Vj分别为节点i和节点j的电压幅值,θij为节点i和节点j的电压相角差;NLCC和NVSC分别表示LCC和VSC直流换流站的交流节点集合;n为系统节点数;
连接LCC直流换流站的交流节点的功率平衡方程如下:
式中,Kpi为节点i连接直流输电系统的极对数;Udi为节点i连接直流换流站的直流侧电压;Idi为节点i连接直流换流站的直流电流;为节点i连接直流换流站的功率因数角;式(4)中的正负号对于整流站为负,逆变站为正;
连接VSC直流换流站的交流节点的功率平衡方程如下:
式中,Pvi和Qvi为节点i连接VSC直流换流站从交流系统吸收的有功和无功;
(2)直流输电线路的运行特性约束
对于LCC直流输电线路,考虑换流变压器和换相电抗的影响,则运行特性约束如下:
UdRi=UdIi+IdiRdc(8c)
式中,Kdi为节点i连接直流换流站的换流变压器变比;θi为节点i连接直流换流站的换流器控制角;Xci为节点i连接直流换流站的换相电抗;Rdc为直流线路电阻;UdRi为整流侧直流电压;UdIi为逆变侧直流电压;
对于VSC直流输电线路,忽略换流站等效电阻的影响,则运行特性约束如下:
式中,μi为节点i连接直流换流站的直流电压利用率,在PWM调制策略确定时为定值;Mi为节点i连接直流换流站的调制比;Vi为节点i连接直流换流站的交流侧电压;Yi为节点i连接直流换流站的等效导纳;δi为节点i连接的交流节点电压与换流器输入电压的相角差;gil,d为消去联络节点后直流网络节点电导矩阵中的元素;m为连接VSC直流换流站的交流节点数。
2.如权利要求1所述的含有风电接入的交直流混联电网静态电压稳定评估方法,其特征在于,所述交直流混联电网SVSM区间计算的双层最优潮流模型包括:
计算SVSM区间上界的min-min模型和计算区间下界的max-min模型两个双层优化模型。
3.如权利要求2所述的含有风电接入的交直流混联电网静态电压稳定评估方法,其特征在于,通过如下方式来计算获得SVSM区间的上/下界及对应的风电场出力值:
通过凸松弛技术将内层非凸的非线性优化模型转化为凸规划模型,并基于对偶优化理论得到内层凸优化模型的对偶优化模型,将计算SVSM区间下界的max-min模型转化为max-max模型;将min-min双层优化模型和max-max双层优化模型转化为单层优化模型求解。
4.如权利要求1所述的含有风电接入的交直流混联电网静态电压稳定评估方法,其特征在于,所述目标函数的约束条件还包括:
(3)发电机有功出力增长方式约束:
假定除平衡节点发电机外的其他发电机按照其有功出力裕度占总有功出力裕度的比例来分配负荷有功的总增长量,如式(10b);对于有功出力已经到达上界的发电机,不再承担负荷有功增长量,因此其有功增长系数为0;
(4)发电机机端电压和无功出力的互补约束:
(QGi-QGimin)Vai=0 (11a)
(QGi-QGimax)Vbi=0 (11b)
VGi=VGiref+Vai-Vbi (11c)
式中,QGimax和QGimin分别为第i台发电机无功出力的上界和下界;VGi和VGiref分别为第i台发电机的机端电压以及其给定参考值;Vai和Vbi是辅助变量,用于修正发电机无功出力到达下界或者上界时的机端电压;
(5)变量的上下限约束:包括发电机节点有功和无功出力的上下限,交流节点电压及辅助变量的上下限,LCC直流换流站的控制角和功率因数角的上下限,及VSC直流换流站的调制比和相角差的上下限,即:
xmin≤x≤xmax (12)
7.如权利要求6所述的含有风电接入的交直流混联电网静态电压稳定评估方法,其特征在于,在悲观模型式(14)中,先将内层min模型转化为其对偶优化max模型,对于如式(16)的非线性规划问题,若f(x)和a(x)都是连续可微的凸函数,则称优化问题式(16)为凸规划问题;并且,式(17)的非线性规划问题为原优化问题式(16)的对偶优化问题:
式中,v=(v1,v2,…,vr)T,vi为对偶变量,r为原优化问题(16)的约束条件的数目;
原问题与对偶问题之间存在着以下关系:如果x0是原问题的最优解,那么存在一组对偶变量v0,使得(x0,v0)是对偶问题的一个最优解,并且,原问题最优解x0所对应的目标函数值f(x0)与对偶问题最优解(x0,v0)所对应的目标函数值g(x0,v0)相同,即:
在将悲观模型式(14)中的内层优化模型转化成其对偶优化模型之前,需要对内层规划问题的约束条件进行凸松弛处理。
8.如权利要求7所述的含有风电接入的交直流混联电网静态电压稳定评估方法,其特征在于,所述凸松弛处理过程如下:
(1)交流线路运行约束的凸松弛
引入中间变量,令Rij=ViVjcosθij,Iij=ViVjsinθij,Ui=Vi 2,此时交流线路的功率平衡方程转化为线性形式如式(19),同时,中间变量应满足约束如式(20),
Rij=Rji (20b)
Iij=-Iji (20c)
由于新增的约束(20a)是非凸的,需要通过二阶锥松弛将其转化为如下的凸不等式约束:
||2Rij;2Iij;Ui-Uj||2≤Ui+Uj (21)
为了使凸化前后的系统等价,需要补充约束条件来保证每个独立回路中各支路两端节点电压相角差之和为0,将此约束转化为中间变量的线性等式如式(23):
发电机端电压和无功的互补方程相应地转化如下:
(QGi-QGimin)Uai=0 (25a)
(QGimax-QGi)Ubi=0 (25b)
令qai=QGiUai,qbi=QGiUbi,则利用凸包络法可将式(25)转换为线性约束如下:
式中,Uaimax和Uaimin表示Uai的上限和下限,Ubimax和Ubimin表示Ubi的上限和下限;
(2)LCC直流线路运行约束的凸松弛
在SVSM计算过程中,当确定了直流输电系统的控制方式后,直流电压和直流电流都是已知量;同时,假定换流变变比不会随着负荷增长而变化,即Kdi为定值,需要对式(8a)和(8b)进行凸松弛:
式(8a)存在非线性项Vicosθi,将cosθi整体看成一个连续变量,则该项为两个连续变量相乘,令wi=Vicosθi,利用凸包络法将其转化为如下线性约束:
式中,Vimax和Vimin表示Vi的上限和下限;
经过上述松弛后,式(8b)可转化为如下线性等式约束:
(3)VSC直流线路运行约束的凸松弛
||2Kvi; 2Lvi; Oi-Ui||2≤Oi+Ui (32)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110468364.2A CN113224788B (zh) | 2021-04-28 | 2021-04-28 | 含有风电接入的交直流混联电网静态电压稳定评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110468364.2A CN113224788B (zh) | 2021-04-28 | 2021-04-28 | 含有风电接入的交直流混联电网静态电压稳定评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113224788A CN113224788A (zh) | 2021-08-06 |
CN113224788B true CN113224788B (zh) | 2023-03-31 |
Family
ID=77089889
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110468364.2A Active CN113224788B (zh) | 2021-04-28 | 2021-04-28 | 含有风电接入的交直流混联电网静态电压稳定评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113224788B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117578551B (zh) * | 2024-01-15 | 2024-03-19 | 四川大学 | 一种连续精确的三端混合直流输电系统的无功调节方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018095123A1 (zh) * | 2016-11-24 | 2018-05-31 | 中国电力科学研究院有限公司 | 确定分布式电源接纳能力的方法和装置、存储介质 |
CN110535121A (zh) * | 2019-08-27 | 2019-12-03 | 东北大学 | 一种交直流混联电网的两阶段多目标动态优化调度方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102983573B (zh) * | 2012-11-09 | 2014-10-15 | 天津大学 | 基于安全域的安全约束经济调度方法 |
CN111541280B (zh) * | 2020-05-14 | 2022-08-09 | 东南大学 | 一种考虑静态电压稳定约束的电网风电最大渗透率评估方法 |
-
2021
- 2021-04-28 CN CN202110468364.2A patent/CN113224788B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018095123A1 (zh) * | 2016-11-24 | 2018-05-31 | 中国电力科学研究院有限公司 | 确定分布式电源接纳能力的方法和装置、存储介质 |
CN110535121A (zh) * | 2019-08-27 | 2019-12-03 | 东北大学 | 一种交直流混联电网的两阶段多目标动态优化调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113224788A (zh) | 2021-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Dynamic modeling and small signal stability analysis of distributed photovoltaic grid-connected system with large scale of panel level DC optimizers | |
Colombino et al. | Towards robustness guarantees for feedback-based optimization | |
Purba et al. | Dynamic aggregation of grid-tied three-phase inverters | |
Shah et al. | Fuzzy logic based FOGI-FLL algorithm for optimal operation of single-stage three-phase grid interfaced multifunctional SECS | |
CN113158573B (zh) | 一种基于深度学习的小水电群区域配网最优潮流获取方法 | |
Cifuentes et al. | Black-box impedance-based stability assessment of dynamic interactions between converters and grid | |
CN110797874A (zh) | 含电力电子变压器的交直流混合配电网状态估计方法 | |
CN115622053B (zh) | 一种用于考虑分布式电源的自动负荷建模方法及装置 | |
Li et al. | A three-phase unbalanced linear power flow solution with PV bus and ZIP load | |
CN113224788B (zh) | 含有风电接入的交直流混联电网静态电压稳定评估方法 | |
CN115021256A (zh) | 一种大规模交直流输电系统电磁暂态模型自动生成方法 | |
Lu et al. | An external network modeling approach for online security analysis | |
CN113612252B (zh) | 含风电场的交直流混联电网静态电压稳定优化控制方法 | |
Girigoudar et al. | Integration of centralized and distributed methods to mitigate voltage unbalance using solar inverters | |
Li et al. | Dispatchable region for active distribution networks using approximate second-order cone relaxation | |
Jo et al. | Kalman-filter-based multilevel analysis to estimate electric load composition | |
CN107230982A (zh) | 一种微电网线性化潮流计算方法 | |
CN113507128B (zh) | 一种特高压直流换流站近场无功优化配置方法 | |
Guo et al. | Online model-free estimation of the dynamic system model for a power system with renewables in ambient conditions | |
CN113258576B (zh) | 交直流互联电网pq节点静态电压稳定性评估方法及系统 | |
Kuri et al. | Dynamic model reduction based on coherency and genetic optimization methodology | |
CN113991652A (zh) | 数据驱动的含iidg配电网短路电流多输出计算方法 | |
Yang et al. | Dynamic Equivalent Research Based on Fuzzy C-means | |
Lu et al. | Impact of series compensation on operation performance of large‐scale PV plants | |
CN110707762A (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 |