CN113688567B - 一种考虑冲击负荷的虚拟电厂两阶段优化调度方法 - Google Patents
一种考虑冲击负荷的虚拟电厂两阶段优化调度方法 Download PDFInfo
- Publication number
- CN113688567B CN113688567B CN202110922966.0A CN202110922966A CN113688567B CN 113688567 B CN113688567 B CN 113688567B CN 202110922966 A CN202110922966 A CN 202110922966A CN 113688567 B CN113688567 B CN 113688567B
- Authority
- CN
- China
- Prior art keywords
- power
- gas turbine
- energy
- load
- daily
- 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 42
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000004146 energy storage Methods 0.000 claims abstract description 40
- 238000003860 storage Methods 0.000 claims abstract description 31
- 230000008859 change Effects 0.000 claims abstract description 28
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 28
- 238000005096 rolling process Methods 0.000 claims abstract description 13
- 239000002245 particle Substances 0.000 claims abstract description 8
- 230000010485 coping Effects 0.000 claims abstract description 4
- 238000012937 correction Methods 0.000 claims abstract description 3
- 230000005611 electricity Effects 0.000 claims description 21
- 230000001105 regulatory effect Effects 0.000 claims description 14
- 238000012423 maintenance Methods 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 12
- 238000011835 investigation Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000010248 power generation Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000008901 benefit Effects 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 5
- 230000014509 gene expression Effects 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 description 5
- 230000000295 complement effect Effects 0.000 description 4
- 230000009194 climbing Effects 0.000 description 3
- 230000001276 controlling effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000007599 discharging Methods 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based 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
- 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)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Evolutionary Computation (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Medical Informatics (AREA)
- Educational Administration (AREA)
- Geometry (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑冲击负荷的虚拟电厂两阶段优化调度方法。本发明的技术方案。步骤包括：首先，建立包括风电、燃气轮机、蓄电池储能的能源供给侧模型，设计冲击负荷下虚拟电厂优化调度的目标函数及约束条件，基于粒子群算法形成日前优化调度计划；然后，日内调度以功率平衡为目的，基于调控时段风电功率与系统净负荷变化的一致性指标，确定日内可控资源的调控量；进一步基于集合经验模态分解将日内调控量分解为高低频分量，高频分量作为蓄电池的功率指令，低频分量作为燃气轮机的功率指令；最后，基于模型预测控制的日内滚动优化校正策略，应对风电及负荷预测误差带来的联络线功率波动，同时确保储能满足日运行能量平衡约束。本发明可通过日前日内两阶段优化调度，实现冲击负荷下的虚拟电厂经济运行的同时，有效的抑制了联络线的功率波动。
Description
技术领域
本发明涉及电力系统领域，特别是一种考虑冲击负荷的虚拟电厂两阶段优化调度方法
技术背景
大力发展和使用可再生分布式电源，提高综合能源利用效率成为当今的热点研究话题。新能源发电具有天然的间歇性和波动性，其大规模并网，必将加剧电源供给和负荷需求的不平衡。虚拟电厂(virtual power plant，VPP)技术为解决此类问题应运而生。VPP利用先进的通信技术，可将不同地域的新能源以及其它分布式电源联系在一起，通过协调优化调度策略，使得出力不稳定的新能源在电力市场竞价中保有一定的竞争力
发明内容
本发明提供一种考虑冲击负荷的虚拟电厂两阶段优化调度方法，其有利于通过虚拟电厂优化调度应对冲击负荷，并提高虚拟电厂运行经济性。
本发明采用技术方案：一种考虑冲击负荷的虚拟电厂两阶段优化调度方法，其包括步骤：
(1)建立包括风电、燃气轮机、蓄电池储能的能源供给侧模型，设计冲击负荷下虚拟电厂优化调度的目标函数及约束条件，基于粒子群形成日前优化调度计划；
(2)日内调度以功率平衡为目的，基于调控时段风电功率与系统净负荷变化的一致性指标，确定日内可控资源的调控量；
(3)基于集合经验模态分解将日内调控量分解为高低频分量，高频分量作为蓄电池的功率指令，低频分量作为燃气轮机的功率指令；
(4)基于模型预测控制的日内滚动优化校正策略，应对风电及负荷预测误差带来的联络线功率波动，同时确保储能满足日运行能量平衡约束。
将多类型分布式资源集成为虚拟电厂，能够发挥不同类型资源的优势，实现资源的互补共济。本发明中并网型虚拟电厂资源包括：风电、燃气轮机、蓄电池储能。
虚拟电厂优化调度以最大化利润为目标，目标函数为：
式中：f1为虚拟电厂利润；IL,t、IE,t、IES,t分别为负荷、电力交易、蓄电池储能收益；CGT,t、CWP,t为燃气轮机、风机运维成本；△t为一个调度时段，日前调度为1h；αL为负荷电价，PL,t为负荷功率；λe为0-1变量，1表示售电，0表示购电，αse、αbe为与大电网售电、购电价格，Pse,t、Pbe,t为与大电网售电、购电功率；λES为0-1变量，1表示放电，0表示充电，αc、αd为充放电电价，PES,t为充放电功率，放电为正，充电为负，μES为蓄电池折旧补偿系数；αGT为燃气轮机运维成本系数，α′GT为燃气轮机应对冲击负荷增发功率运维成本系数，PGT,t、△PGT,t为燃气轮机发电功率、应对冲击负荷增发功率，λGT为0-1变量，1表示起机，0表示停机，Cq、Ct为启停机费用；αWP为风机运维成本系数，PWP,t风机功率。
虚拟电厂运行主要包括以下约束条件
(1)虚拟电厂功率平衡约束
式中：PGT,t、PWP,t、PL,t分别为燃气轮机、风电、负荷功率；为蓄电池充放电功率；为虚拟电厂从大电网购入、出售功率。
(2)燃气轮机约束
其中：为燃气轮机出力上下限；为燃气轮机功率上升、下降的爬坡率。
(3)蓄电池约束
式(5)为储能设备运行功率约束，其中为蓄电池最大充放电功率；为0-1变量，表征储能设备的充放电状态。
式(6)为蓄电池剩余容量约束条件，其中，EES,t、EES,t-△t分别为t时刻、t-△t时刻蓄电池剩余容量；EES,min、EES,max为剩余容量上下限；为充放电功率；ηES为充放电效率；△t为计算步长。
虚拟电厂运行以经济性最优为目标，考虑约束条件，采用粒子群算法，完成日前优化调度，得到燃气轮机，风电，蓄电池，联络线功率出力计划。
模型预测控制是一种基于模型的闭环优化控制方法，其算法的核心思想是滚动时域策略。该策略主要包括以下几步：
(1)在当前时刻k和当前状态x(k)，基于一定的预测模型，预测系统未来的状态，同时考虑当前和未来的约束条件，通过解决优化问题，得到未来k+1,k+2,L,k+M时刻的控制指令序列；
(2)将控制指令序列的第一个值应用于控制系统；
(3)在k+1时刻，更新状态为x(k+1)，重复上述操作步骤。
可见，经济模型预测控制就是不断滚动的局部优化和不断滚动实施控制作用的交替过程。
具体到虚拟电厂日内优化调度问题，根据虚拟电厂每时段的功率平衡方程及储能剩余电量迭代方程，选取燃气轮机出力、储能充放电功率、储能剩余电量，以及虚拟电厂与外电网的联络线交换功率构成的向量x(k)＝[PMT(k),PES(k),EES(k),PE(k)]为状态变量；以燃气轮机和蓄电池储能的出力增量构成的向量u(k)＝[△PMT(k),△PES(k)]为控制变量；以负荷、风电机组超短期预测功率增量构成的向量r(k)＝[△PL(k),△PWP(k)]为扰动输入；以燃气轮机出力及联络线交换功率构成的向量y(k)＝[PMT(k),PE(k)]为输出变量，则可建立如下的多输入、多输出状态空间模型：
其中
基于风电和负荷的超短期功率预测数据，通过对该状态空间预测模型反复迭代，直到向前预测m步，便可得到燃气轮机和联络线功率在预测时长m△t内的预估输出值构成的向量Y，表述如下式所示：
Y＝[PMT(k+△t),PE(k+△t),L,PMT(k+m△t),PE(k+m△t)] (8)
取当前时刻向前m△t时段内，燃气轮机出力和联络线功率日前计划值构成的向量G为跟踪控制目标，其中G可描述如下：
接着，以燃气轮机和联络线功率预估输出值与日前计划值之间的误差最小为目标，同时保证各机组控制调节增量尽量小，这样便可将对应的日内滚动优化调度转化为如下所示的二次规划问题：
J＝(G-Y)TW(G-Y)+uTQu (10)
其中为燃气轮机和联络线功率跟踪误差的权重系数矩阵；Q为控制量的权重系数矩阵。
虚拟电厂优化日内优化调度调控策略制定的基本前提条件是：在判断有关考察指标是否超出阈值及是否进行后续优化调控的初期阶段，调控系统根据日前优化调度计划可获得未来时段风电功率整体预测值、电网负荷预测值、燃气轮机发电计划，以作为优化调控建模和调控量计算的基础信息；系统内有一定容量的燃气轮机旋转备用、储能设备资源，作为调控量再分配的具体调控措施或对象。
定义虚拟电厂调控一致性指标，某功率信号i的变化速率为：
分别为信号i在t时刻、t-1时刻的采样值；△t为计算功率变化速率周期；n为考察时窗包含的总采样点数，时窗跨度为n△t。
功率信号i变化速率的相对值为：
式中：Pib为信号i功率的最大值，对待研究的风电场功率信号，该值取装机容量，对负荷功率信号，该值可取日最大负荷或更长时间级的最大负荷值。
对功率信号i、j定义集合
一致性指标定义为：
n为考察时窗包含的总采样点数，IB越接近于0，说明2个功率信号在考察时窗内的变化特性越一致；反之，两信号变化的不一致性越高。
基于功率平衡的虚拟电厂调控策略的目的是为了冲击负荷下使系统发电功率与负荷需求尽可能地保持平衡。
调控初始阶段，利用风电功率、负荷的预测信息及燃气轮机发电计划，计算下一调控时段中第t时刻系统功率不平衡量，第t-1时刻系统功率为平衡状态，即：
式中：分别为负荷、联络线功率、风电出力、燃气轮机出力的变化量；为时段t相对于时段t-1的待调控量，没有调控时，其值为0。
功率变化化为相对值：
式中：PWb为日最大负荷值
考虑系统功率平衡需求时，风电功率与负荷变化一致性指标阈值为
εpmin为功率不平衡量限制，根据系统频率误差及待考察区域的系统功–频静特性系数整定。
该调控时段风电功率与系统净负荷变化的一致性指标为
n为考察时窗包含的总采样点数。当IAL≤IALcr时，认为该调控时段风电功率的变化情况满足系统功率平衡要求。
基于EEMD的调控功率分配，大功率脉动负荷有短时的峰值，功率分布不均，使得待调控功率存在较大的功率波动。燃气轮机出力平稳，能量密度高，能很好满足电力高峰时的需求，但燃气轮机受爬坡率限制，不能实现快速的升降负荷；蓄电池储能在轻负荷时将能量存储起来，在冲击负荷到来时释放能量，蓄电池放电迅速，短时间内可以放出较大电量，可应对频繁波动的冲击负荷。
采用基于集合经验模态分解Ensemble Empirical Mode Decomposition,EEMD方法，实现调控功率的自适应分解，得到燃气轮机出力变化和储能系统充、放电功率指令，对原始调控功率进行EEMD分解步骤如下：
(1)在原始功率信号中加入高斯白噪声信号δt，得到一个常用信号△Pt即：
(2)对信号△Pt进行EMD分解，得到固有模态函数IMF分量
式中，为残余信号；为第i(i＝1,2,L,n)层固有模态函数IMF，并按照频率从高到低的顺序进行分布。
(3)高斯白噪声信号第j次加入原始信号重复进行步骤(1)、步骤(2)
式中，为第j次加入高斯白噪声后分解得到的第i个IMF分量。
(4)根据不相关随机序列统计均值为零的原理，对进行整体平均，使得多次加入的白噪声对IMF的影响相互抵消，则最终功率信号可分解为i个IMF分量。
在进行EEMD分解时，加入的高斯白噪声信号数量满足统计规律。
式中，ε为高斯白噪声信号的幅值；N为加入高斯白噪声信号的整体平均次数；εn为原始信号与经EEMD分解后重构信号之间的误差。从式(19)中可知：若高斯白噪声信号幅值ε选取一定时，整体平均次数N越大，εn越小，即高斯白噪声信号对原始风功率信号的影响越小，最终分解的IMF经过重构后越接近原始信号；若N不变，ε选取越小，会使信噪比变高，削弱噪声对原始信号高频部分产生的互补作用。
(5)经EEMD分解后重构信号可表示为：
固有模态函数IMF本身是具有能量的，从能量熵的角度出发，根据EEMD的算法，能量越大的IMF分量所占权重比越大。因此在分量和余项表达式求解之后，引入能量熵的概念，将IMF能量差值最大的作为一个分界点。假设每一阶IMF能量熵分别为Ei，总的能量为E，在公式计算中不加入余项的能量。
式中N为采样点数。将IMF能量做归一化处理得到:
相邻两个IMF的能量熵的差值越大，表明其能量占比也就越大，用于区别高频段和低频段。通过相邻IMF能量熵的变换值作为一个标准阶次k，将熵值变化最大的阶次作为此时高、低频率对应的标准阶次k，将k之前的高频分量作为蓄电池的功率指令，k之后的低频分量作为燃气轮机的功率指令。
式中：为蓄电池功率指令；为燃气轮机功率指令；n为IMF的总数。
本发明提供的技术方案具有的有益效果如下：通过引入虚拟电厂技术，将区域风电、燃气轮机与储能设备集成为虚拟电厂，设计冲击负荷下虚拟电厂日前优化调度，日内调度考虑风电功率与系统净负荷变化的一致性指标，利用集合经验模态分解分配调控量功率指令，运用模型预测制定日内调度计划，实现虚拟电厂经济运行的同时，有效的抑制了联络线的功率波动。
附图说明
下面结合附图对本发明作进一步说明：
图1为本发明的流程图；
图2为集合经验模态分解算法步骤流程图；
图3为风电负荷预测曲线图；
图4为虚拟电厂日前优化调度曲线图；
图5为燃气轮机日前日内出力曲线图；
图6为联络线功率日前日内曲线图；
具体实施方案
为了更好地理解本发明的目的、技术方案及技术效果，以下结合附图对本发明进行进一步的讲解说明。
将多类型分布式资源集成为虚拟电厂，能够发挥不同类型资源的优势，实现资源的互补共济。本发明中并网型虚拟电厂资源包括，风电、燃气轮机、蓄电池储能。
本发明提出了一种考虑冲击负荷的虚拟电厂两阶段优化调度方法，如图1所示，其实施流程包括如下详细步骤。
步骤1建立包括风电、燃气轮机、蓄电池储能的能源供给侧模型，设计冲击负荷下虚拟电厂优化调度的目标函数及约束条件，基于粒子群形成日前优化调度计划；
虚拟电厂优化调度以最大化利润为目标，目标函数为：
式中：f1为虚拟电厂利润；、IE,t、IES,t分别为负荷、电力交易、蓄电池储能收益；CGT,t、CWP,t为燃气轮机、风机运维成本；△t为一个调度时段，日前调度为1h；αL为负荷电价，PL,t为负荷功率；λe为0-1变量，1表示售电，0表示购电，αse、αbe为与大电网售电、购电价格，Pse,t、Pbe,t为与大电网售电、购电功率；λES为0-1变量，1表示放电，0表示充电，αc、αd为充放电电价，PES,t为充放电功率，放电为正，充电为负，μES为蓄电池折旧补偿系数；αGT为燃气轮机运维成本系数，α′GT为燃气轮机应对冲击负荷增发功率运维成本系数，PGT,t、△PGT,t为燃气轮机发电功率、应对冲击负荷增发功率，λGT为0-1变量，1表示起机，0表示停机，Cq、Ct为启停机费用；αWP为风机运维成本系数，PWP,t风机功率。
虚拟电厂运行主要包括以下约束条件
(1)虚拟电厂功率平衡约束
式中：PGT,t、PWP,t、PL,t分别为燃气轮机、风电、负荷功率；为蓄电池充放电功率；为虚拟电厂从大电网购入、出售功率。
(2)燃气轮机约束
其中：为燃气轮机出力上下限；为燃气轮机功率上升、下降的爬坡率。
(3)蓄电池约束
式(5)为储能设备运行功率约束，其中为蓄电池最大充放电功率；βt c、βtd为0-1变量，表征储能设备的充放电状态。
式(6)为蓄电池剩余容量约束条件，其中，EES,t、EES,t-△t分别为t时刻、t-△t时刻蓄电池剩余容量；EES,min、EES,max为剩余容量上下限；为充放电功率；ηES为充放电效率；△t为计算步长。
虚拟电厂运行以经济性最优为目标，考虑约束条件，采用粒子群算法，完成日前优化调度，得到燃气轮机，风电，蓄电池，联络线功率出力计划。
模型预测控制是一种基于模型的闭环优化控制方法，其算法的核心思想是滚动时域策略。该策略主要包括以下几步：
(1)在当前时刻k和当前状态x(k)，基于一定的预测模型，预测系统未来的状态，同时考虑当前和未来的约束条件，通过解决优化问题，得到未来k+1,k+2,L,k+M时刻的控制指令序列；
(2)将控制指令序列的第一个值应用于控制系统；
(3)在k+1时刻，更新状态为x(k+1)，重复上述操作步骤。
可见，经济模型预测控制就是不断滚动的局部优化和不断滚动实施控制作用的交替过程。
具体到虚拟电厂日内优化调度问题，根据虚拟电厂每时段的功率平衡方程及储能剩余电量迭代方程，选取燃气轮机出力、储能充放电功率、储能剩余电量，以及虚拟电厂与外电网的联络线交换功率构成的向量x(k)＝[PMT(k),PES(k),EES(k),PE(k)]为状态变量；以燃气轮机和蓄电池储能的出力增量构成的向量u(k)＝[△PMT(k),△PES(k)]为控制变量；以负荷、风电机组超短期预测功率增量构成的向量r(k)＝[△PL(k),△PWP(k)]为扰动输入；以燃气轮机出力及联络线交换功率构成的向量y(k)＝[PMT(k),PE(k)]为输出变量，则可建立如下的多输入、多输出状态空间模型：
其中
基于风电和负荷的超短期功率预测数据，通过对该状态空间预测模型反复迭代，直到向前预测m步，便可得到燃气轮机和联络线功率在预测时长m△t内的预估输出值构成的向量Y，表述如下式所示：
Y＝[PMT(k+△t),PE(k+△t),L,PMT(k+m△t),PE(k+m△t)](8)
取当前时刻向前m△t时段内，燃气轮机出力和联络线功率日前计划值构成的向量G为跟踪控制目标，其中G可描述如下：
接着，以燃气轮机和联络线功率预估输出值与日前计划值之间的误差最小为目标，同时保证各机组控制调节增量尽量小，这样便可将对应的日内滚动优化调度转化为如下所示的二次规划问题：
J＝(G-Y)TW(G-Y)+uTQu(10)
其中为燃气轮机和联络线功率跟踪误差的权重系数矩阵；Q为控制量的权重系数矩阵。
步骤2日内调度以功率平衡为目的，基于调控时段风电功率与系统净负荷变化的一致性指标，确定日内可控资源的调控量；
虚拟电厂优化日内优化调度调控策略制定的基本前提条件是：在判断有关考察指标是否超出阈值及是否进行后续优化调控的初期阶段，调控系统根据日前优化调度计划可获得未来时段风电功率整体预测值、电网负荷预测值、燃气轮机发电计划，以作为优化调控建模和调控量计算的基础信息；系统内有一定容量的燃气轮机旋转备用、储能设备资源，作为调控量再分配的具体调控措施或对象。
定义虚拟电厂调控一致性指标，某功率信号i的变化速率为：
分别为信号i在t时刻、t-1时刻的采样值；△t为计算功率变化速率周期；n为考察时窗包含的总采样点数，时窗跨度为n△t。
功率信号i变化速率的相对值为：
式中：Pib为信号i功率的最大值，对待研究的风电场功率信号，该值取装机容量，对负荷功率信号，该值可取日最大负荷或更长时间级的最大负荷值。
对功率信号i、j定义集合
一致性指标定义为：
n为考察时窗包含的总采样点数，IB越接近于0，说明2个功率信号在考察时窗内的变化特性越一致；反之，两信号变化的不一致性越高。
基于功率平衡的虚拟电厂调控策略的目的是为了冲击负荷下使系统发电功率与负荷需求尽可能地保持平衡。
调控初始阶段，利用风电功率、负荷的预测信息及燃气轮机发电计划，计算下一调控时段中第t时刻系统功率不平衡量，第t-1时刻系统功率为平衡状态，即：
式中：分别为负荷、联络线功率、风电出力、燃气轮机出力的变化量；为时段t相对于时段t-1的待调控量，没有调控时，其值为0。
功率变化化为相对值：
式中：PWb为日最大负荷值
考虑系统功率平衡需求时，风电功率与负荷变化一致性指标阈值为
εpmin为功率不平衡量限制，根据系统频率误差及待考察区域的系统功–频静特性系数整定。
该调控时段风电功率与系统净负荷变化的一致性指标为
n为考察时窗包含的总采样点数。当IAL≤IALcr时，认为该调控时段风电功率的变化情况满足系统功率平衡要求。
步骤3基于集合经验模态分解将日内调控量分解为高低频分量，高频分量作为蓄电池的功率指令，低频分量作为燃气轮机的功率指令；
基于集合经验模态分解Ensemble Empirical Mode Decomposition,EEMD的调控功率分解流程图如图2，大功率脉动负荷有短时的峰值，功率分布不均，使得待调控功率存在较大的功率波动。燃气轮机出力平稳，能量密度高，能很好满足电力高峰时的需求，但燃气轮机受爬坡率限制，不能实现快速的升降负荷；蓄电池储能在轻负荷时将能量存储起来，在冲击负荷到来时释放能量，蓄电池放电迅速，短时间内可以放出较大电量，可应对频繁波动的冲击负荷。
采用基于EEMD方法，实现调控功率的自适应分解，得到燃气轮机出力变化和储能系统充、放电功率指令，对原始调控功率进行EEMD分解步骤如下：
(1)在原始功率信号中加入高斯白噪声信号δt，得到一个常用信号△Pt即：
(2)对信号△Pt进行EMD分解，得到固有模态函数IMF分量
式中，为残余信号；为第i(i＝1,2,L,n)层固有模态函数IMF，并按照频率从高到低的顺序进行分布。
(3)高斯白噪声信号第j次加入原始信号重复进行步骤(1)、步骤(2)
式中，为第j次加入高斯白噪声后分解得到的第i个IMF分量。
(4)根据不相关随机序列统计均值为零的原理，对进行整体平均，使得多次加入的白噪声对IMF的影响相互抵消，则最终功率信号可分解为i个IMF分量。
在进行EEMD分解时，加入的高斯白噪声信号数量满足统计规律。
式中，ε为高斯白噪声信号的幅值；N为加入高斯白噪声信号的整体平均次数；εn为原始信号与经EEMD分解后重构信号之间的误差。从式(19)中可知：若高斯白噪声信号幅值ε选取一定时，整体平均次数N越大，εn越小，即高斯白噪声信号对原始风功率信号的影响越小，最终分解的IMF经过重构后越接近原始信号；若N不变，ε选取越小，会使信噪比变高，削弱噪声对原始信号高频部分产生的互补作用。
(5)经EEMD分解后重构信号可表示为：
固有模态函数IMF本身是具有能量的，从能量熵的角度出发，根据EEMD的算法，能量越大的IMF分量所占权重比越大。因此在分量和余项表达式求解之后，引入能量熵的概念，将IMF能量差值最大的作为一个分界点。假设每一阶IMF能量熵分别为Ei，总的能量为E，在公式计算中不加入余项的能量。
式中N为采样点数。将IMF能量做归一化处理得到:
相邻两个IMF的能量熵的差值越大，表明其能量占比也就越大，用于区别高频段和低频段。通过相邻IMF能量熵的变换值作为一个标准阶次k，将熵值变化最大的阶次作为此时高、低频率对应的标准阶次k，将k之前的高频分量作为蓄电池的功率指令，k之后的低频分量作为燃气轮机的功率指令。
式中：为蓄电池功率指令；为燃气轮机功率指令；n为IMF的总数。
为了进一步理解本发明，以某开发区能源系统作为算例。
系统内包括30MW风电装机，蓄电池储能为20MW·h，燃气轮机50MW；将负荷端用电情况分为峰平谷3个时段，总时间段为一天24小时。其中峰时段为10至15和18至21时段、平时段为7至10时段、15至18时段和21至24时段、谷时段为0至7时段，峰平谷3个时段电价分别为1.055、0.633、0.291元/kW·h。虚拟电厂单元其它参数如表1。
图3为虚拟电厂负荷及风电预测曲线，根据其预测结果，运用粒子群算法得到的日前优化调度结果如图4描述。
根据图4可知，储能装置主要在低电价和低负荷时段进行充电，在高电价或高负荷时段则进行放电，既保证了储能日运行SOC的平衡，又有效发挥了储能的“削峰填谷”作用，提高了系统运行的经济性；微型燃气轮机的运行成本较低，故微型燃气轮机在大部分时段保持较高的出力水平；同时通过与大电网进行交易，合理利用储能进行“套利”，提高系统运行经济性。
日内模型预测控制取预测时长为1h，控制时长为30min，滚动优化调度执行周期为5min一次，在日内一共滚动优化288次，滚动优化调整后得到的燃气轮机、联络线优化出力计划日前日内对比如图5、图6所示。
由于日前计划以小时为时间尺度，调度粗放，系统不能及时响应可再生能源及负荷波动。施加MPC滚动优化调度后，联络线功率与日前计划值基本吻合，率波动较小，实现微电网接入配电网的平稳、可控调度。由图6燃气轮机日前日内功率曲线可知，受燃气轮机爬坡率约束，日内调度计划难以较好跟踪日前计划，存在一定的功率偏差，而储能蓄电池能够很好的弥补这一不足，各能源转换设备能更快速响应系统功率波动，系统协调能力加强，更容易达到平衡。且对每次滚动优化的计算时长进行统计发现，基于MPC的单次优化计算时长仅为0.2s左右，完全可满足在线应用的需求，也充分说明了本文协调优化调度方案的有效性。
Claims (3)
1.一种考虑冲击负荷的虚拟电厂两阶段优化调度方法，其特征在于，包括步骤：
(1)建立包括风电、燃气轮机、蓄电池储能的能源供给侧模型，设计冲击负荷下虚拟电厂优化调度的目标函数及约束条件，基于粒子群算法形成日前优化调度计划；
(2)日内调度以功率平衡为目的，基于调控时段风电功率与系统净负荷变化的一致性指标，确定日内可控资源的调控量；
(3)基于集合经验模态分解将日内调控量分解为高低频分量，高频分量作为蓄电池的功率指令，低频分量作为燃气轮机的功率指令；
(4)基于模型预测控制的日内滚动优化校正策略，应对风电及负荷预测误差带来的联络线功率波动，同时确保储能满足日运行能量平衡约束；
其中步骤(1)中虚拟电厂日前调度以利润最大化为目标，目标函数为：
式中：f1为虚拟电厂利润，IL,t、IE,t、IES,t分别为负荷、电力交易、蓄电池储能收益，CGT,t、CWP,t为燃气轮机、风机运维成本，△t为一个调度时段，日前调度为1h，αL为负荷电价，PL,t为负荷功率，λe为0-1变量，1表示售电，0表示购电，αse、αbe为与大电网售电、购电价格，Pse,t、Pbe,t为与大电网售电、购电功率，λES为0-1变量，1表示放电，0表示充电，αc、αd为充放电电价，PES,t为充放电功率，放电为正，充电为负，μES为蓄电池折旧补偿系数，αGT为燃气轮机运维成本系数，α′GT为燃气轮机应对冲击负荷增发功率运维成本系数，PGT,t、△PGT,t为燃气轮机发电功率、应对冲击负荷增发功率，λGT为0-1变量，1表示起机，0表示停机，Cq、Ct为启停机费用，αWP为风机运维成本系数，PWP,t风机功率；
虚拟电厂运行主要包括以下约束条件：
(1)虚拟电厂功率平衡约束
式中：PGT,t、PWP,t、PL,t分别为燃气轮机、风电、负荷功率，为蓄电池充放电功率，为虚拟电厂从大电网购入、出售功率；
(2)燃气轮机约束
其中：为燃气轮机出力上下限，为燃气轮机功率上升、下降的爬坡率；
(3)蓄电池约束
式(5)为储能设备运行功率约束，其中为蓄电池最大充放电功率，为0-1变量，表征储能设备的充放电状态；
式(6)为蓄电池剩余容量约束条件，其中，EES,t、EES,t-△t分别为t时刻、t-△t时刻蓄电池剩余容量，EES,min、EES,max为剩余容量上下限，为充放电功率，ηES为充放电效率，△t为计算步长；
虚拟电厂运行以经济性最优为目标，考虑约束条件，采用粒子群算法，完成日前优化调度，得到燃气轮机，风电，蓄电池，联络线功率出力计划；
所述步骤(4)中模型预测控制方法如下：
模型预测控制具体到虚拟电厂日内优化调度问题如下：根据虚拟电厂每时段的功率平衡方程及储能剩余电量迭代方程，选取燃气轮机出力、储能充放电功率、储能剩余电量，以及虚拟电厂与外电网的联络线交换功率构成的向量x(k)＝[PMT(k),PES(k),EES(k),PE(k)]为状态变量，以燃气轮机和蓄电池储能的出力增量构成的向量u(k)＝[△PMT(k),△PES(k)]为控制变量，以负荷、风电机组超短期预测功率增量构成的向量r(k)＝[△PL(k),△PWP(k)]为扰动输入，以燃气轮机出力及联络线交换功率构成的向量y(k)＝[PMT(k),PE(k)]为输出变量，则可建立如下的多输入、多输出状态空间模型：
其中
基于风电和负荷的超短期功率预测数据，通过对该状态空间预测模型反复迭代，直到向前预测m步，便可得到燃气轮机和联络线功率在预测时长m△t内的预估输出值构成的向量Y，表述如下式所示：
Y＝[PMT(k+△t),PE(k+△t),L,PMT(k+m△t),PE(k+m△t)] (8)
取当前时刻向前m△t时段内，燃气轮机出力和联络线功率日前计划值构成的向量G为跟踪控制目标，其中G可描述如下：
接着，以燃气轮机和联络线功率预估输出值与日前计划值之间的误差最小为目标，同时保证各机组控制调节增量尽量小，这样便可将对应的日内滚动优化调度转化为如下所示的二次规划问题：
J＝(G-Y)TW(G-Y)+uTQu(10)
其中为燃气轮机和联络线功率跟踪误差的权重系数矩阵，Q为控制量的权重系数矩阵。
2.根据权利要求1所述的一种考虑冲击负荷的虚拟电厂两阶段优化调度方法，其特征在于，所述步骤(2)中，根据一致性指标数学模型，确定日内可控资源的调控量：
调控初始阶段，利用风电功率、负荷的预测信息及燃气轮机发电计划，计算下一调控时段中第t时刻系统功率不平衡量，第t-1时刻系统功率为平衡状态，即：
式中：分别为负荷、联络线功率、风电出力、燃气轮机出力的变化量，为时段t相对于时段t-1的待调控量；
功率变化为相对值：
式中：PWb为日最大负荷值
考虑系统功率平衡需求时，风电功率与负荷变化一致性指标阈值为：
εpmin为功率不平衡量限制，根据系统频率误差及待考察区域的系统功–频静特性系数整定；
该调控时段风电功率与系统净负荷变化的一致性指标为
n为考察时窗包含的总采样点数，当IAL≤IALcr时，认为该调控时段风电功率的变化情况满足系统功率平衡要求。
3.根据权利要求1所述的一种考虑冲击负荷的虚拟电厂两阶段优化调度方法，其特征在于，所述步骤(3)中基于集合经验模态分解的调控功率分配：
对原始调控功率进行集合经验模态分解Ensemble Empirical Mode Decomposition,EEMD步骤如下：
(1)在原始功率信号中加入高斯白噪声信号δt，得到一个常用信号△Pt即：
(2)对信号△Pt进行分解，得到固有模态函数IMF分量
式中，为残余信号，为第i(i＝1,2,L,n)层固有模态函数IMF，并按照频率从高到低的顺序进行分布；
(3)高斯白噪声信号第j次加入原始信号重复进行步骤(1)、步骤(2)
式中，为第j次加入高斯白噪声后分解得到的第i个IMF分量；
(4)根据不相关随机序列统计均值为零的原理，对进行整体平均，使得多次加入的白噪声对IMF的影响相互抵消，则最终功率信号可分解为i个IMF分量；
在进行EEMD分解时，加入的高斯白噪声信号数量满足统计规律；
式中，ε为高斯白噪声信号的幅值，N为加入高斯白噪声信号的整体平均次数，εn为原始信号与经EEMD分解后重构信号之间的误差；
(5)经EEMD分解后重构信号可表示为：
引入能量熵的概念，将IMF能量差值最大的作为一个分界点，假设每一阶IMF能量熵分别为Ei，总的能量为E，在公式计算中不加入余项的能量；
式中N为采样点数，将IMF能量做归一化处理得到:
相邻两个IMF的能量熵的差值越大，表明其能量占比也就越大，用于区别高频段和低频段，通过相邻IMF能量熵的变换值作为一个标准阶次k，将熵值变化最大的阶次作为此时高、低频率对应的标准阶次k，将k之前的高频分量作为蓄电池的功率指令，k之后的低频分量作为燃气轮机的功率指令；
式中：为蓄电池功率指令，为燃气轮机功率指令，n为IMF的总数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110922966.0A CN113688567B (zh) | 2021-08-10 | 2021-08-10 | 一种考虑冲击负荷的虚拟电厂两阶段优化调度方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110922966.0A CN113688567B (zh) | 2021-08-10 | 2021-08-10 | 一种考虑冲击负荷的虚拟电厂两阶段优化调度方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113688567A CN113688567A (zh) | 2021-11-23 |
| CN113688567B true CN113688567B (zh) | 2023-08-11 |
Family
ID=78579534
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110922966.0A Active CN113688567B (zh) | 2021-08-10 | 2021-08-10 | 一种考虑冲击负荷的虚拟电厂两阶段优化调度方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113688567B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114037191A (zh) * | 2022-01-07 | 2022-02-11 | 深圳江行联加智能科技有限公司 | 基于大数据的虚拟电厂优化调度方法、装置、设备及介质 |
| CN115021327B (zh) * | 2022-05-19 | 2024-03-01 | 合肥工业大学 | 基于虚拟储能的互联微电网系统的协同日前优化调度方法 |
| CN116231765B (zh) * | 2023-05-09 | 2023-09-19 | 上海融和元储能源有限公司 | 一种虚拟电厂出力控制方法 |
| CN116599023B (zh) * | 2023-05-17 | 2024-02-27 | 中国人民解放军海军工程大学 | 一种船用燃气轮机电网混合储能功率分配方法及设备 |
| CN116796146B (zh) * | 2023-05-24 | 2024-06-07 | 国家电网有限公司华东分部 | 抽蓄电厂的电量平衡方法、装置、设备及存储介质 |
| CN116488223A (zh) * | 2023-06-26 | 2023-07-25 | 湖南大学 | 家庭光-蓄-柔双层多时间尺度控制方法、装置及介质 |
| CN117200184B (zh) * | 2023-08-10 | 2024-04-09 | 国网浙江省电力有限公司金华供电公司 | 一种虚拟电厂负荷侧资源多时段调控潜力评估预测方法 |
| CN117411083A (zh) * | 2023-10-31 | 2024-01-16 | 东北电力大学 | 一种考虑电弧炉调控与风电模态分解的两阶段调度方法 |
| CN117748501B (zh) * | 2024-02-19 | 2024-05-07 | 西安热工研究院有限公司 | 一种储能辅助黑启动的风功率预测方法和系统 |
| CN117791742B (zh) * | 2024-02-27 | 2024-05-14 | 广东电网有限责任公司中山供电局 | 基于物联网的风电机组的发电功率的控制方法和装置 |
| CN118014315B (zh) * | 2024-04-08 | 2024-07-05 | 太湖能谷(杭州)科技有限公司 | 一种用户侧共享储能系统参与辅助服务的调度方法和系统 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007028158A2 (en) * | 2005-09-02 | 2007-03-08 | Lightridge Resources Llc | Energy and chemical species utility management system |
| CN103617455A (zh) * | 2013-11-29 | 2014-03-05 | 广东电网公司电力科学研究院 | 基于虚拟机组子群的网厂两级负荷优化调度方法 |
| CN103824134A (zh) * | 2014-03-06 | 2014-05-28 | 河海大学 | 一种用于虚拟电厂的两阶段优化调度方法 |
| CN108388973A (zh) * | 2018-01-11 | 2018-08-10 | 河海大学 | 一种虚拟电厂自适应鲁棒调度优化方法 |
| CN110620380A (zh) * | 2019-09-05 | 2019-12-27 | 深圳供电局有限公司 | 一种虚拟电厂控制系统 |
| CN111342500A (zh) * | 2020-03-19 | 2020-06-26 | 国网福建省电力有限公司宁德供电公司 | 小水电虚拟电厂的多时间尺度优化调度方法 |
| CN111784055A (zh) * | 2020-07-02 | 2020-10-16 | 长沙理工大学 | 一种考虑电力需求侧响应特性的虚拟电厂优化调度模型 |
| CN113177323A (zh) * | 2021-05-14 | 2021-07-27 | 华北电力大学 | 基于矩不确定分布式鲁棒的电热综合系统优化调度方法 |
| CN115603317A (zh) * | 2022-11-01 | 2023-01-13 | 国网湖北省电力有限公司襄阳供电公司（Cn） | 一种基于两阶段风险约束的虚拟电厂最优调度方法 |
| CN115622020A (zh) * | 2022-10-21 | 2023-01-17 | 国网电力科学研究院有限公司 | 一种面向直流负荷的虚拟电厂系统及运行方法 |
-
2021
- 2021-08-10 CN CN202110922966.0A patent/CN113688567B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007028158A2 (en) * | 2005-09-02 | 2007-03-08 | Lightridge Resources Llc | Energy and chemical species utility management system |
| CN103617455A (zh) * | 2013-11-29 | 2014-03-05 | 广东电网公司电力科学研究院 | 基于虚拟机组子群的网厂两级负荷优化调度方法 |
| CN103824134A (zh) * | 2014-03-06 | 2014-05-28 | 河海大学 | 一种用于虚拟电厂的两阶段优化调度方法 |
| CN108388973A (zh) * | 2018-01-11 | 2018-08-10 | 河海大学 | 一种虚拟电厂自适应鲁棒调度优化方法 |
| CN110620380A (zh) * | 2019-09-05 | 2019-12-27 | 深圳供电局有限公司 | 一种虚拟电厂控制系统 |
| CN111342500A (zh) * | 2020-03-19 | 2020-06-26 | 国网福建省电力有限公司宁德供电公司 | 小水电虚拟电厂的多时间尺度优化调度方法 |
| CN111784055A (zh) * | 2020-07-02 | 2020-10-16 | 长沙理工大学 | 一种考虑电力需求侧响应特性的虚拟电厂优化调度模型 |
| CN113177323A (zh) * | 2021-05-14 | 2021-07-27 | 华北电力大学 | 基于矩不确定分布式鲁棒的电热综合系统优化调度方法 |
| CN115622020A (zh) * | 2022-10-21 | 2023-01-17 | 国网电力科学研究院有限公司 | 一种面向直流负荷的虚拟电厂系统及运行方法 |
| CN115603317A (zh) * | 2022-11-01 | 2023-01-13 | 国网湖北省电力有限公司襄阳供电公司（Cn） | 一种基于两阶段风险约束的虚拟电厂最优调度方法 |
Non-Patent Citations (1)
| Title |
|---|
| 李星雨 等.考虑虚拟电厂和分时电价的风光火储系统两阶段优化调度策略.《南方电网技术》.2017,第11卷(第6期),第70-77页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113688567A (zh) | 2021-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113688567B (zh) | 一种考虑冲击负荷的虚拟电厂两阶段优化调度方法 | |
| CN110188950B (zh) | 基于多代理技术的虚拟电厂供电侧和需求侧优化调度建模方法 | |
| Yang et al. | Optimal sizing and placement of energy storage system in power grids: A state-of-the-art one-stop handbook | |
| Athari et al. | Operational performance of energy storage as function of electricity prices for on-grid hybrid renewable energy system by optimized fuzzy logic controller | |
| CN111738497B (zh) | 计及需求侧响应的虚拟电厂双层优化调度方法 | |
| CN112465181A (zh) | 一种支持源-网-荷-储多元泛在协调的两阶段优化调度方法 | |
| CN111934360B (zh) | 基于模型预测控制的虚拟电厂-储能系统能量协同优化调控方法 | |
| CN104376385A (zh) | 一种微电网电价优化方法 | |
| CN112800658A (zh) | 一种考虑源储荷互动的主动配电网调度方法 | |
| Maulik | Probabilistic power management of a grid-connected microgrid considering electric vehicles, demand response, smart transformers, and soft open points | |
| CN111697578A (zh) | 多目标含储能区域电网运行控制方法 | |
| CN112508325B (zh) | 一种家庭微电网多时间尺度电能调度方法 | |
| CN112821470A (zh) | 基于小生境混沌粒子群算法的微网群优化调度策略 | |
| Zhang et al. | Energy optimal management of microgrid with high photovoltaic penetration | |
| CN114926254A (zh) | 储能电站参与调频辅助服务市场的竞价方法 | |
| CN116109076A (zh) | 能量和调峰市场下考虑需求响应的虚拟电厂优化调度方法 | |
| CN112598232A (zh) | 一种改进量子粒子群的vpp中储能系统优化配置方法 | |
| CN116307505A (zh) | 一种企业微电网能量经济优化调度方法 | |
| CN115498624A (zh) | 基于两阶段鲁棒优化的含灵活互动资源配电网低碳经济调度方法 | |
| CN117578537A (zh) | 基于碳交易和需求响应的微电网优化调度方法 | |
| CN105930919A (zh) | 基于两阶段随机规划的虚拟电厂风险规避优化运行方法 | |
| Mohamed et al. | Residential battery energy storage sizing and profitability in the presence of PV and EV | |
| Ren et al. | Multi-objective optimal dispatching of virtual power plants considering source-load uncertainty in V2G mode | |
| Dai et al. | An equilibrium model of the electricity market considering the participation of virtual power plants | |
| CN114493143A (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 |