CN106295915B - 具有最大容量准则约束的含清洁能源电网优化调度方法 - Google Patents
具有最大容量准则约束的含清洁能源电网优化调度方法 Download PDFInfo
- Publication number
- CN106295915B CN106295915B CN201610883400.0A CN201610883400A CN106295915B CN 106295915 B CN106295915 B CN 106295915B CN 201610883400 A CN201610883400 A CN 201610883400A CN 106295915 B CN106295915 B CN 106295915B
- Authority
- CN
- China
- Prior art keywords
- power
- node
- energy
- maximum
- determined
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 238000004146 energy storage Methods 0.000 claims description 42
- 230000005611 electricity Effects 0.000 claims description 20
- 230000005540 biological transmission Effects 0.000 claims description 6
- 239000012141 concentrate Substances 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 238000005457 optimization Methods 0.000 description 15
- 238000010248 power generation Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 238000004140 cleaning Methods 0.000 description 4
- 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 3
- 238000003860 storage Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 206010044565 Tremor Diseases 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000012010 growth Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Classifications
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
-
- 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
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种具有最大容量准则约束的含清洁能源电网优化调度方法,步骤包括:步骤1,根据电网模型确定系统的拉普拉斯矩阵;步骤2,确定参数;步骤3,确定电网中各个节点的最大变化容量约束;步骤4,确定目标函数和待优化参数;步骤5,基于最大容量准则确定约束条件;步骤6,采用NSGA‑II在线滚动算法对上述多目标问题进行优化。本发明的方法,给出了电网能够维持稳定的清洁能源接入容量的极限值,在实现合理调度的基础上,保证了系统的运行稳定性。
Description
技术领域
本发明属于电网优化调度技术领域,涉及一种具有最大容量准则约束的含清洁能源电网优化调度方法。
背景技术
近年来,世界各国清洁能源如太阳能、风能等装机规模都快速增长,预计到2020年,中国清洁能源发电容量占总装机容量的28.4%,然而这些清洁能源发电具有间歇性与“不可控性”,当其大规模并网运行时,严重影响现有电网的电能质量,甚至导致电压、频率崩溃。如何协调全网可调度资源,实现资源的高效配置、清洁能源发电的合理消纳、电网的稳定运行,是电力行业面临的重大挑战。
传统发电计划的编制是以常规能源出力的可靠性和负荷预测的准确性为基础。然而,清洁能源发电具有随机性和间歇性,为确保并网后电力系统的稳定运行,必须对调度计划进行调整和优化。当清洁能源发电容量在总发电总量中所占比例较小时,可将清洁能源发电出力的预测值作为“负”的负荷与原始负荷曲线叠加,形成净负荷曲线,然后进行发电计划的优化编制;然而,当清洁能源的并网容量占总发电需求的比例达到一定值时,清洁能源的可信赖电力和电量贡献不可忽视。这种情况下,现有的一些调度方法[1-6]并不能保证电力系统的频率稳定性。
发明内容
本发明的目的是提供一种具有最大容量准则约束的含清洁能源电网优化调度方法,解决了现有技术在大量清洁能源并网过程中可能导致电网频率失稳的问题,给出了保证电网稳定前提下可接入波动功率的最大变化(调度)容量。
本发明所采用的技术方案是,一种具有最大容量准则约束的含清洁能源电网优化调度方法,按照以下步骤实施:
步骤1,根据电网模型确定系统的拉普拉斯矩阵
包含风/光/储节点的电网频率分析模型:
其中,i,j分别为电网中第i,j个节点的编号,n为电网中节点总数,Mi为节点i的惯性常数,Di为节点i的阻尼常数,θi,θj分别为节点i,j的相角,Ei,Ej分别为节点i,j的电压,Yij为节点i与节点j之间的导纳,gi(t)为节点i的波动功率;根据上述模型得到系统的拉普拉斯矩阵为:
其中,αij=EiEj|Yij|为两节点之间的耦合强度;
步骤2,确定参数
将步骤1所确定的拉普拉斯矩阵L去掉零特征值及其零特征值对应的特征向量后,重构为一个新矩阵La,解方程PLa+La TP=-I求得唯一解P,其中I为单位矩阵,进而得到:
c1=λmin(P),c2=λmax(P),c3=λmin(I)=1,c4=2λmax(P), (3)
其中,λmin(P),λmax(P),λmin(I)分别表示P的最小特征值、P的最大特征值、I的最小特征值;c1,c2,c3,c4同是步骤3为确定电网中各个节点的容许波动功率的最大范围|gi(k)|所需的参数;
步骤3,确定电网中各个节点的最大变化容量约束|gi(k)|
电网中各个节点的最大变化容量约束:
将步骤2确定的四个参数c1,c2,c3,c4代入式(4),得到式(5):
其中,为节点i的物理最大传输功率,Pij max=EiEj|Yij|为边εij的物理最大传输功率,r为允许频率偏差,Di为节点i的阻尼常数,k表示第k个时刻;
步骤4,确定目标函数和待优化参数
以常规机组出力与储能单元充放电功率为优化变量,确定待优化参数数量为m,针对电网运行过程中要获得的不同目标选取合适的目标函数J1,J2,…,Jq,q为待优化目标函数个数;
步骤5,基于最大容量准则确定约束条件
5.1)有功平衡约束:
其中,min函数表示取两个变量中较小的一个,ΩG为常规机组节点的集合,ΩB为储能节点的集合,ΩV为光伏节点的集合,ΩW为风电节点的集合,ΩL为负荷节点的集合, 分别为k时刻常规机组ig、光伏节点iv、风电节点iw的功率;为k时刻储能单元ib的充放电功率,为k时刻的负荷节点il的消耗功率,分别为k时刻光伏节点iv、风电节点iw、负荷节点il的功率预测值;
5.2)常规机组的出力约束:
其中,分别为常规机组出力上、下界;
5.3)储能单元的功率约束:
其中,放电为正,充电为负;分别为储能单元ib的最大充、放电功率;
步骤6,采用NSGA-II在线滚动算法对上述多目标问题进行优化
由于优化变量的时变性,因此以传统机组每小时出力、储能单元每小时充放电功率为优化变量;采用NSGA-II在线滚动算法,经过多次优化即能够给出某一天的功率调度方案。
本发明的有益效果是,不仅考虑了清洁能源如太阳能、风能、负荷等波动功率,而且考虑各个节点波动功率的范围对于电网稳定性的影响,在电网的能量优化调度中引入最大波动功率限制,解决了在大量清洁能源并网过程中可能引起频率失稳的问题。与一般的电网优化调度方法相比,本发明方法基于最大容许容量准则约束,考虑了清洁能源波动功率对稳定性的影响,解决了现行优化调度研究中不考虑稳定性影响的问题。本发明能够在保证电网稳定的情况下,计算出各节点容许容量变化的最大值,并能够与传统约束条件相结合以便合理安排各常规机组的出力及储能单元充放电功率,为电力行业的自主、智能发展提供了理论指导。
附图说明
图1是IEEE14的节点标准测试系统连线图;
图2是一天中6:00-18:00时间段内太阳能、风能、负荷的预测波动功率;
图3是基于最大容量准则给出的风电节点容许波动功率;
图4是基于最大容量准则给出的光伏节点容许波动功率;
图5是基于最大容量准则给出的负荷节点容许波动值;
图6是具有上述最大容量准则约束的电网优化调度方法给出的三个储能节点充放电功率调度方案;
图7是具有上述最大容量准则约束的电网优化调度方法给出的常规机组出力方案;
图8是具有上述最大容量准则约束的电网优化调度方案下电网系统的频率变化图;
图9是不考虑最大容量准则约束的一般调度方法给出的三个储能节点的功率调度方案;
图10是不考虑最大容量准则约束的一般优化调度方法给出的常规机组出力方案;
图11是不考虑最大容量准则约束的一般优化调度方案下电网系统的频率变化图。
具体实施方式
下面结合附图及具体实施例对本发明进行进一步的说明。
本发明方法通过李亚普诺夫逆定理给出各个节点容许波动功率的最大值,并将该理论运用在电力系统的功率优化调度中,在合理安排各常规机组出力及储能单元充放电的前提下,能够保证电网的运行稳定性,实现电网的智能调控。能够在基于最大容量准则的约束条件下合理安排各常规机组出力及储能单元充放电功率,实现清洁能源发电功率的合理消纳。
下面对基于二阶swing方程的电网模型及节点容许波动功率界限值的计算进行简要的描述与分析:
根据方程式(1)给出的包含风/光/储节点的电网频率分析模型:
其中,h为光伏发电机输出功率的幅值,γ为描述光伏发电机输出功率波动的参量,ρ为空气密度,cp为风能利用系数,λ为叶尖速比,θp为桨距角,Ar为扫风面积,vω为风速;Pmi为传统发电机机械功率,Pdi为储能单元存储功率;对于太阳能发电机/储能,Mi=0。
将式(1)写为矩阵形式,则有:
其中, 为方程式(10)的非线性部分,Gi(t)=gi(t)/Mi。
则方程式(10)相当于式(11)的扰动系统:
假设式(11)的扰动系统存在一个指数稳定平衡点;则存在一个Lyapunov函数V满足下列不等式:
其中,参数c1,c2,c3,c4均为正实数。
继续使用为式(10)扰动系统的备选Lyapunov函数,则沿式(10)扰动系统轨线的导数满足:
将不等式(13)、式(14)一起代入式(15)得:
上述方程式第一项小于0,只要第二项与第三项之和小于0,则
即:
由等式(12)得到:在集合之外,
即:在集合之外,
则:在集合之外,
其中,r为允许频率偏差。
定义:始于该集合Ωρ的解在未来时刻不会离开该集合,因为在边界V=ρ上是负的,
由于:则有
对于任意时刻t,则得到:
将式(18)代入式(17)后得:
当则式(10)扰动系统的指数稳定。
图1为IEEE14节点标准测试系统连线图,节点2、节点7分别为风电、光伏节点,节点4、节点5、节点6为储能节点,节点1为常规发电机节点,其它节点为负荷节点。
参照图1,结合前述的实施例对象描述,本发明的方法,按照以下步骤具体实施:
步骤1,根据电网模型确定系统的拉普拉斯矩阵
根据方程式(1)给出的包含风/光/储节点的电网频率分析模型:
其中,i,j分别为电网中第i,j个节点的编号,n为电网中节点总数,Mi为节点i的惯性常数,Di为节点i的阻尼常数,θi,θj分别为节点i,j的相角,Ei,Ej分别为节点i,j的电压;Yij为节点i与节点j之间的导纳,gi(t)为节点i的波动功率。
当研究式(1)所代表的电网系统的频率同步稳定性时,惯性项只会影响整个系统达到同步的时间,不会影响同步状态的存在,因此忽略掉惯性项后得一阶模型:将式(20)线性化得:写成矩阵形式为:
其中,则称:
为系统的拉普拉斯矩阵;其中,αij=EiEj|Yij|。
结合图1的线路数据[7]最终可得系统的拉普拉斯矩阵如下:
步骤2,确定参数
由于式(2)所表示的拉普拉斯矩阵L是行和为0的Metzle矩阵,因此系统除零特征值外,其余特征值都具有严格的负实部,首先去掉矩阵L的零特征值及其零特征值对应的特征向量,然后重构一个矩阵La如下:
该矩阵的所有特征值都具有负实部,再解方程PLa+LaTP=-I,其中I为单位矩阵,可求得唯一解P如下:
则存在一个李雅普诺夫函数V(t,θ)=θTPθ满足下式:
λmin(P)||θ||2≤V(t,θ)≤λmax(P)||θ||2, (23)
对照式(12)、式(13)、式(14)得到下式(3):
c1=λmin(P),c2=λmax(P),c3=λmin(I)=1,c4=2λmax(P), (3)
然后求P的特征值,最终得到四个参数的值:
c1=λmin(P)=0.0026,c2=λmax(P)=0.1357,c3=λmin(I)=1,c4=2λmax(P)=0.2714。
步骤3,确定电网中各个节点的最大变化容量约束(容许波动功率的最大范围)|gi(k)|
式(19)是基于李亚普诺夫逆定理给出的清洁能源的允许波动功率的界限值。因此将步骤2中根据式(3)确定的参数c1,c2,c3,c4代入式(19)中,得到波动功率界限值:
为了考虑波动功率与风/光/储节点的接入位置及其接入数量的关系,定义:
其中,为节点i的物理最大传输功率,Pij max=EiEj|Yij|为边εij的物理最大传输功率,对于发电机节点,Pi(t-1)为t-1时刻发电功率,对于储能节点,Pi(t-1)为t-1时刻充放电功率,对于负载节点,Pi(t-1)为t-1时刻消耗功率。
则得到各个节点的最大功率波动范围为:
由于在上述实施例中:Gi(t)=gi(t)/Di, (29)
并且为了与后面步骤中调度时间一致,将t换作k表示离散的时间,于是得:
其中,r=0.5为允许频率偏差,Di为节点i的阻尼常数,k表示第k个时刻。结合图1的线路数据及节点数据[7]最终计算得到表1如下:
表1、第一次迭代中各节点波动功率的界限值
节点i | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
|g<sub>i</sub>(1)| | 4.2682 | 6.6809 | 2.0912 | 7.9845 | 7.2726 | 4.0980 | 4.2475 |
节点i | 8 | 9 | 10 | 11 | 12 | 13 | 14 |
|g<sub>i</sub>(1)| | 1.2994 | 3.0243 | 1.0246 | 2.0088 | 1.4516 | 2.6796 | 1.2349 |
步骤4,确定目标函数和待优化参数
以一个常规机组出力与三个储能单元充放电功率为优化变量。对于图1为IEEE14的节点标准测试系统,本实施例仅优化上午6点至下午18点的调度方案(其它时刻的优化以此类似进行),滚动优化次数12次,确定待优化参数的数量为(4*12),设待优化目标函数为2个,则有:
优化目标1:传统发电机与储能单元的运行成本最小:
其中,为传统发电机的成本函数;为储能单元的成本函数,为k时刻常规机组发电功率;为k时刻储能单元的充放电功率;aGi,bGi,cGi,aBi为成本函数的参量:aGi=0.04,bGi=2,cGi=0;aBi=0.18。
优化目标2:实时调度功率最小:
其中,为k时刻的调度功率。
该实施例中常规机组数目1个、储能单元数目3个、成本函数的参量取值为:aGi=0.04,bGi=2,cGi=0;aBi=0.18,代入式(30)、式(31)得到目标函数为:
步骤5,基于最大容量准则确定约束条件
5.1)有功平衡约束:
其中,min函数表示取两个变量中较小的一个,ΩG为常规机组节点的集合,ΩB为储能节点的集合,ΩV为光伏节点的集合,ΩW为风电节点的集合,ΩL为负荷节点的集合, 分别为k时刻常规机组ig、光伏节点iv、风电节点iw的功率;为k时刻储能单元ib的充放电功率,为k时刻的负荷节点il的负荷功率,分别为k时刻光伏节点ig、风电节点iw、负荷节点il的功率预测值。
5.2)常规机组出力约束:
其中,为常规机组出力上、下界。
5.3)储能单元的运行约束:
功率约束:
其中,放电为正,充电为负;分别为储能单元ib最大充、放电功率。
该实施例中约束条件的具体取值如下:
传统发电机出力:
储能单元充放电功率:
步骤6,采用NSGA-II在线滚动算法对上述多目标问题进行优化
由于优化变量的时变性,因此以常规机组每小时出力、储能单元每小时充放电功率为优化变量;对于多目标优化问题,可能会产生多个最优解组成Pareto最优解集,因此再根据实际情况,选出实际需求的最优解。
本实施例中优化时段为6点到18点(时间间隔1小时),共有优化变量(4*12)个。本实施例中,采用NSGA-II在线滚动算法,共经过12次优化(或根据需要进行24次优化),每次优化产生10个Pareto解组成的解集,然后在这个解集中按如下原则选择一组最优解:
6.1)出于节能、环保考虑,首先以常规机组出力最小为原则,假如解集中所有常规机组出力相同则执行下一步;
6.2)以储能单元1的调度功率最小为原则,假如解集中所有储能单元1的调度功率相同,则以储能单元2的调度功率最小为原则,以此类推,直到找出本次的最优解。
本实施例第一次优化产生10个Pareto解组成的解集如表2所示,然后按照上述步骤得到第3组解:即为本次优化的最优解;其他时刻最优解的选择也如上述步骤所示。经过12次优化得到出6点到18点这个时间段内的功率调度方案如图3、图4所示。
本发明将最大容量准则与优化调度相结合,采用NSGA-II在线滚动算法,在清洁能源大范围并网运行的情况下,实现功率的合理调度,并保证了电网运行的安全与稳定性。
表2、第一次优化产生的Pareto解集
下面对于一般优化调度方法进行简要的描述:
以1个常规机组出力与3个储能单元充放电功率为优化变量。优化6点至18点的调度方案,确定待优化参数数量为(4*12),设待优化目标函数2个,与步骤4中的相同。
约束条件:
1)有功平衡约束:
其中,ΩG为常规机组节点的组合,ΩB为储能节点的组合,ΩV为光伏节点的组合,ΩW为风电节点的组合,ΩL为负荷节点的组合,分别为k时刻常规机组ig、光伏节点iv、风电节点iv的功率;为k时刻储能单元ib的充放电功率,为k时刻的负荷节点il的负荷功率,分别为k时刻光伏节点iv、风电节点iw、负荷节点il的功率预测值。
2)常规机组出力约束:
其中,为常规机组出力上、下界。
3)储能单元的功率约束:
其中,放电为正,充电为负;分别为储能单元ib最大充、放电功率。
约束条件的具体取值如下:
采用NSGA-II在线滚动算法,共经过12次优化,每次优化产生10个最优解,然后按照步骤6中给出的原则在最优解集中选出最优解得出6点到18点这个时间段内的功率调度方案如图9、图10所示。
验证本发明方法的效果
图2是某一天6:00—18:00时间段内对风能、太阳能、负荷等三种不同类型节点波动功率的预测图。负荷预测参考文献[6]:负荷功率波动一般会在±20%以内。图3、图4、图5分别为基于最大容量准则给出在IEEE14标准测试系统的线路数据和节点数据下的风能发电机、光伏发电机、负荷可消纳的波动功率。图6、图7为具有最大容量准则约束的电网优化调度方法给出的储能节点和常规机组出力的调度方案,图8和图11分别为电网在本发明优化调度方法和一般不考虑最大功率波动约束优化调度方法下电网的频率变化图,在仿真的过程中,由于频率稳定性的时间尺度较小,因此,1小时时间的稳定性可以简单地在100s时间内体现,可以简化仿真过程。图8和图11都是每隔100s,变化一次调度功率,总共变化了12次。由图8可以看出,上述优化调度方案不仅实现了一般优化调度的功能,而且能够确保电网的频率稳定性。
图9、图10为在图2的所示清洁能源波动功率预测值下,一般不考虑最大功率波动约束优化调度方法给出的储能节点、常规机组出力的功率调度方案,图11为一般不考虑最大功率波动约束优化调度方案下电网的频率变化图。
对比图3、图4与图11可见,在图3、图4中,由于清洁能源节点注入功率在11点到13点时较大超出稳定范围,分析可知在这个时间段内电网会出现较大的频率不稳,则可预测在图11中,第500s到700s(对应图3、图4中的11点到13点)可能会出现频率不稳定情况。而在随后的时间段内,图3、图4中节点注入功率基本满足稳定性要求,因此,可预测在图11中,电网频率在随后的时间段内是同步稳定的。事实上,图11的结果与上述预测分析结果一致,因此也验证了本发明方法的正确性。对于实际电网,由于受到线路参数的限制,电网接纳清洁能源的容量是有限,超出容量限制,会导致电能质量快速下降,甚至导致不稳定。现行解决方法是通过多次仿真试验得出电网稳态运行时清洁能源电站的出力,逐渐增加清洁能源电站出力,同时观察电网电能质量是否越限,当电能质量不满足要求时,采用保守调度,往往使得清洁能源电站的出力小于稳定极限值运行。本发明方法是在一般优化调度的基础上,计算出电网能够接纳清洁能源的最大容量,当清洁能源发电量预测值在最大容量范围内,采用基于预测调度的方法,假如清洁能源发电量预测值超出最大容量范围,采用我们给出的最大容量约束的电网调度方法,因此,本发明给出了电网能够维持稳定的清洁能源接入容量的极限值,解决了现行方法通过仿真验证很难穷举、可靠性差的问题,在实现合理调度的基础上,保证了系统的运行稳定性。
Claims (2)
1.一种具有最大容量准则约束的含清洁能源电网优化调度方法,其特征在于,按照以下步骤实施:
步骤1,根据电网模型确定系统的拉普拉斯矩阵
包含风/光/储节点的电网频率分析模型:
其中,i,j分别为电网中第i,j个节点的编号,n为电网中节点总数,Mi为节点i的惯性常数,Di为节点i的阻尼常数,θi,θj分别为节点i,j的相角,Ei,Ej分别为节点i,j的电压,Yij为节点i与节点j之间的导纳,gi(t)为节点i的波动功率;根据上述模型得到系统的拉普拉斯矩阵为:
其中,αij=EiEj|Yij|为两节点之间的耦合强度;
步骤2,确定参数
将步骤1所确定的拉普拉斯矩阵L去掉零特征值及其零特征值对应的特征向量后,重构为一个新矩阵La,解方程PLa+La TP=-I求得唯一解P,其中I为单位矩阵,进而得到:
c1=λmin(P),c2=λmax(P),c3=λmin(I)=1,c4=2λmax(P), (3)
其中,λmin(P),λmax(P),λmin(I)分别表示P的最小特征值、P的最大特征值、I的最小特征值;c1,c2,c3,c4同是步骤3为确定电网中各个节点的容许波动功率的最大范围|gi(k)|所需的参数;
步骤3,确定电网中各个节点的最大变化容量约束|gi(k)|
电网中各个节点的最大变化容量约束:
将步骤2确定的四个参数c1,c2,c3,c4代入式(4),得到式(5):
其中,为节点i的物理最大传输功率,Pij max=EiEj|Yij|为边εij的物理最大传输功率,r为允许频率偏差,Di为节点i的阻尼常数,k表示第k个时刻;
步骤4,确定目标函数和待优化参数
以常规机组出力与储能单元充放电功率为优化变量,确定待优化参数数量为m,针对电网运行过程中要获得的不同目标选取合适的目标函数J1,J2,…,Jq,q为待优化目标函数个数;
步骤5,基于最大容量准则确定约束条件,
5.1)有功平衡约束:
其中,min函数表示取两个变量中较小的一个,ΩG为常规机组节点的集合,ΩB为储能节点的集合,ΩV为光伏节点的集合,ΩW为风电节点的集合,ΩL为负荷节点的集合,
分别为k时刻常规机组ig、光伏节点iv、风电节点iw的功率;为k时刻储能单元ib的充放电功率,为k时刻的负荷节点il的消耗功率,分别为k时刻光伏节点iv、风电节点iw、负荷节点il的功率预测值;
5.2)常规机组的出力约束:
其中,分别为常规机组出力上、下界;
5.3)储能单元的功率约束:
其中,放电为正,充电为负;分别为储能单元ib的最大充、放电功率;
步骤6,采用NSGA-II在线滚动算法对上述多目标问题进行优化。
2.根据权利要求1所述的具有最大容量准则约束的含清洁能源电网优化调度方法,其特征在于:所述的步骤6中,采用NSGA-II在线滚动算法,共经过12次优化,每次优化产生10个Pareto解组成的解集,然后在这个解集中按如下原则选择一组最优解:
6.1)以常规机组出力最小为原则,假如解集中所有常规机组出力相同则执行下一步;
6.2)以储能单元1的调度功率最小为原则,假如解集中所有储能单元1的调度功率相同,则以储能单元2的调度功率最小为原则,以此类推,直到找出本次的最优解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610883400.0A CN106295915B (zh) | 2016-10-10 | 2016-10-10 | 具有最大容量准则约束的含清洁能源电网优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610883400.0A CN106295915B (zh) | 2016-10-10 | 2016-10-10 | 具有最大容量准则约束的含清洁能源电网优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106295915A CN106295915A (zh) | 2017-01-04 |
CN106295915B true CN106295915B (zh) | 2019-05-28 |
Family
ID=57718101
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610883400.0A Active CN106295915B (zh) | 2016-10-10 | 2016-10-10 | 具有最大容量准则约束的含清洁能源电网优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106295915B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112418533A (zh) * | 2020-11-25 | 2021-02-26 | 江苏电力交易中心有限公司 | 一种清洁能源电量分解预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102545258A (zh) * | 2012-01-12 | 2012-07-04 | 西安交通大学 | 一种大规模风电并网的电网优化规划方法 |
CN103915839A (zh) * | 2014-04-08 | 2014-07-09 | 华北电力大学 | 一种用于分析含风电电力系统随机稳定性的方法 |
CN104852373A (zh) * | 2015-04-29 | 2015-08-19 | 南方电网科学研究院有限责任公司 | 一种含间歇性能源独立微电网的动态经济调度优化方法 |
CN105117797A (zh) * | 2015-08-13 | 2015-12-02 | 浙江工业大学 | 基于孤岛划分的微网日前经济最优调度方法 |
CN105160440A (zh) * | 2015-10-09 | 2015-12-16 | 国家电网公司 | 考虑电网经济调度的储能平抑波动功率的容量规划方法 |
-
2016
- 2016-10-10 CN CN201610883400.0A patent/CN106295915B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102545258A (zh) * | 2012-01-12 | 2012-07-04 | 西安交通大学 | 一种大规模风电并网的电网优化规划方法 |
CN103915839A (zh) * | 2014-04-08 | 2014-07-09 | 华北电力大学 | 一种用于分析含风电电力系统随机稳定性的方法 |
CN104852373A (zh) * | 2015-04-29 | 2015-08-19 | 南方电网科学研究院有限责任公司 | 一种含间歇性能源独立微电网的动态经济调度优化方法 |
CN105117797A (zh) * | 2015-08-13 | 2015-12-02 | 浙江工业大学 | 基于孤岛划分的微网日前经济最优调度方法 |
CN105160440A (zh) * | 2015-10-09 | 2015-12-16 | 国家电网公司 | 考虑电网经济调度的储能平抑波动功率的容量规划方法 |
Non-Patent Citations (1)
Title |
---|
"风光储联合发电系统的组合建模与等值";王皓怀等;《中国电机工程学报》;20111205;第31卷(第34期);第1-9页 |
Also Published As
Publication number | Publication date |
---|---|
CN106295915A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Krishan et al. | Grid-independent PV system hybridization with fuel cell-battery/supercapacitor: Optimum sizing and comparative techno-economic analysis | |
Sedghi et al. | Optimal storage planning in active distribution network considering uncertainty of wind power distributed generation | |
Ma et al. | Optimal allocation of hybrid energy storage systems for smoothing photovoltaic power fluctuations considering the active power curtailment of photovoltaic | |
Injeti et al. | Optimal planning of distributed generation for improved voltage stability and loss reduction | |
CN110581571A (zh) | 一种主动配电网动态优化调度方法 | |
CN106130004A (zh) | 一种综合考虑稳定特性的并网点新能源接纳能力的评估方法 | |
Wang et al. | An improved min-max power dispatching method for integration of variable renewable energy | |
Beus et al. | A model predictive control approach to operation optimization of an ultracapacitor bank for frequency control | |
CN107492903A (zh) | 一种基于统计学模型的混合储能系统容量优化配置方法 | |
Xu et al. | A hybrid energy storage strategy based on multivariable fuzzy coordinated control of photovoltaic grid‐connected power fluctuations | |
CN106295915B (zh) | 具有最大容量准则约束的含清洁能源电网优化调度方法 | |
Akbari-Zadeh et al. | Dstatcom allocation in the distribution system considering load uncertainty | |
CN110323779B (zh) | 一种分布式发电与储能装置的功率动态聚合的方法和系统 | |
Liu et al. | A comprehensive control strategy for photovoltaic virtual synchronous generator considering frequency regulation capability | |
Zhang et al. | Optimization method of under frequency load shedding schemes for systems with high permeability new energy | |
Di et al. | Multi-objective collaborative control scheduling optimization considering wind power grid-connected energy storage access | |
Wang et al. | Research on Coordinated Reactive Power and Voltage Control Strategy for Regional Power Grids with High Penetration of Renewable Energy | |
Ding et al. | Multi-Objective optimial configuration of distributed wind-solar generation considering energy storage | |
Kan et al. | Optimal configuration of the hybrid energy storage system for reducing the amount of discarded photovoltaic | |
Abegaz et al. | Optimal dispatch control of energy storage systems using forward-backward induction | |
CN114239252B (zh) | 电力系统运行方式生成方法、系统、计算机及存储介质 | |
CN111162565B (zh) | 基于多源数据融合的中低压网络在线拼接方法及系统 | |
Cai et al. | PV cluster partition based coordinated voltage control for active distribution network with high penetrated PVs | |
CN111130101B (zh) | 一种多端口能量路由器多场景容量配置优化方法 | |
Abd et al. | Power Grid Frequency Control Improvement Using Pumped Storage Based Se-PSO/LADRC. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |