CN105896596B - 一种考虑需求侧响应的风电功率分层平滑系统及其方法 - Google Patents
一种考虑需求侧响应的风电功率分层平滑系统及其方法 Download PDFInfo
- Publication number
- CN105896596B CN105896596B CN201610241244.8A CN201610241244A CN105896596B CN 105896596 B CN105896596 B CN 105896596B CN 201610241244 A CN201610241244 A CN 201610241244A CN 105896596 B CN105896596 B CN 105896596B
- Authority
- CN
- China
- Prior art keywords
- load
- energy storage
- storage device
- translatable
- 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.)
- Expired - Fee Related
Links
- 238000009499 grossing Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000004146 energy storage Methods 0.000 claims abstract description 129
- 238000005457 optimization Methods 0.000 claims abstract description 59
- 230000005611 electricity Effects 0.000 claims abstract description 54
- 238000007599 discharging Methods 0.000 claims abstract description 39
- 238000013178 mathematical model Methods 0.000 claims abstract description 26
- 240000002853 Nelumbo nucifera Species 0.000 claims abstract description 11
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims abstract description 11
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims abstract description 11
- 239000000203 mixture Substances 0.000 claims description 7
- 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 claims description 6
- 238000007493 shaping process Methods 0.000 claims description 6
- 239000002245 particle Substances 0.000 abstract description 14
- 230000002452 interceptive effect Effects 0.000 abstract description 2
- 230000003993 interaction Effects 0.000 description 5
- 238000013519 translation Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- CWYNVVGOOAEACU-UHFFFAOYSA-N Fe2+ Chemical compound [Fe+2] CWYNVVGOOAEACU-UHFFFAOYSA-N 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000004378 air conditioning Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009415 formwork Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- VWDWKYIASSYTQR-UHFFFAOYSA-N sodium nitrate Chemical group [Na+].[O-][N+]([O-])=O VWDWKYIASSYTQR-UHFFFAOYSA-N 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- H02J3/386—
-
- 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/28—Arrangements for balancing of the load in a network by storage of energy
- H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
-
- 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
- Y02E70/00—Other energy conversion or management systems reducing GHG emissions
- Y02E70/30—Systems combining energy storage with energy generation of non-fossil origin
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑需求侧响应的风电功率分层平滑系统及其方法，是按如下步骤进行：1建立源‑荷‑储互动分层平滑模型，包括风电分时电价数学模型、储能装置的充放电状态模型、负荷分层优化风电功率平滑模型、储能装置充放电状态优化模型；2基于PSO粒子群优化算法对模型进行求解，求得可平移负荷运行调度以及储能装置充放电的最佳方案。本发明能够使负荷有效“跟踪”风电机组出力，降低联络线上风电功率输出的波动性，使风电功率输出更加平滑，并提高系统内用户购电的经济性。
Description
技术领域
本发明涉及需求侧响应技术领域，更具体地说就是用户侧负荷进行优化调度，减小大规模可再生能源接入电网的波动性，维持电网稳定运行的优化系统和方法。
背景技术
随着世界经济的不断发展，能源与环境问题日益突出，可再生能源已成为国家能源战略的重要组成部分。大规模的分布式电源的并网，其间歇性和波动性对电网的电能质量以及稳定性造成负面影响。研究表明用户负荷与储能作为需求侧资源的重要组成部分，可以有效地参与电网系统的运行优化，对维持电网安全、稳定的运行起着至关重要的作用。随着需求侧响应技术的不断发展，人们对电力的供求关系提出了更高的要求，用户根据电力需求和电力系统满足其需求的能力来调整其消费有助于平衡电力供求关系。然而，单一的储能技术并不能有效解决日益增长的电力需求与有限的电力资源之间的矛盾。为了充分利用电网中各种资源，实现良性互动，让用户侧更多地参与电网的优化运行，应从需求侧响应的角度，深入研究源-荷-储互动的风电功率平滑方法。
目前，平滑风电功率的方法一般是：利用超级电容储能与蓄电池储能组成混合储能系统，通过对运行控制方式的设计，使得该储能系统能够与风电系统进行功率交换；联合调度风电机组和常规机组，建立含风电机组的系统动态调度模型，使机组之间的出力互相协调。但单纯的用储能或者联合调度机组平滑风电功率，都只是从发电侧的角度优化控制，目前发电侧备用资源十分有限，这些方法不能有效的解决增长的电力需求与有限的电力资源之间的矛盾。
发明内容
本发明为克服上述现有技术中存在的不足之处，提供一种考虑需求侧响应的风电功率分层平滑系统及其方法，以期能够使负荷有效“跟踪”风电机组出力，降低供电侧联络线上风电功率输出的波动性，使风电功率平滑输出，从而维持电力系统的稳定运行，提高对用电侧系统内用户的经济性。
是把负荷作为需求侧资源的重要组成部分，通过分层优化多时间尺度可平移负荷，以实现风电功率的平滑输出。
为了达到上述目的，本发明所采用的技术方案为：
本发明一种考虑需求侧响应的风电功率分层平滑系统的特点是包括：风电机组、储能装置、参与需求侧响应的负荷群、风电预测模块、负荷调度模块、负荷优化模块；所述负荷群包括：不可控负荷和可控负荷中的可平移负荷；
所述风电预测模块获取所述风电机组的历史出力数据并预测未来出力数据，从而建立分时电价数学模型；
所述负荷调度模块获取所述可平移负荷的负荷参数信息，并按时间尺度的降序顺序对可平移负荷进行分层处理，获得N层可平移负荷曲线；记为{L1,L2,…,Ln,…,LN}；Ln表示第n层可平移负荷曲线；1≤n≤N；同时，获取所述不可控负荷的负荷参数信息并得到不可控负荷曲线，记为L*；
所述负荷调度模块从所述风电预测模块获取所述未来风电机组出力数据和分时电价，并建立所述储能装置的充放电状态数学模型；
所述负荷调度模块根据所述未来风电机组出力数据和不可控负荷曲线L*获得N层可平移负荷调度的目标曲线；
所述负荷优化模块根据所述不可控负荷曲线L*、分时电价数学模型和N层可平移负荷调度的目标曲线，建立可平移负荷的N层经济性目标函数；
所述负荷优化模块根据N层可平移负荷曲线、分时电价数学模型和N层可平移负荷调度的目标曲线，建立可平移负荷的N层波动性目标函数；
所述负荷优化模块根据N层可平移负荷曲线和N层可平移负荷调度的目标曲线，建立可平移负荷的N层平滑性目标函数；
所述负荷优化模块建立由所述N层经济性目标函数、N层波动性目标函数和N层平滑性目标函数，以及所设定的可平移负荷的时间范围约束组成的负荷分层优化风电功率平滑模型；
所述负荷优化模块求取所述负荷分层优化风电功率平滑模型的适应度函数，并利用PSO算法对所述负荷分层优化风电功率平滑模型进行求解，获得优化后的N层可平移负荷曲线；从而使得所述负荷群得到优化；
所述负荷调度模块根据所述风电机组未来出力数据、不可控负荷曲线L*和优化后的N层可平移负荷曲线获得储能装置的目标曲线；
所述负荷调度模块根据所述储能装置充放电状态数学模型、储能装置的目标曲线和分时电价数学模型获得储能装置的波动性目标函数；
所述负荷调度模块建立由所述储能装置的目标曲线、储能装置的波动性目标函数以及所设定的充放电约束条件组成的所述储能装置的充放电状态优化模型；
所述负荷优化模块利用PSO算法对所述储能装置的充放电状态优化模型进行求解，获得优化后的储能装置的充放电状态，从而使得所述储能装置得到优化；
所述风电功率分层平滑系统通过联络线实现所述风电机组的输出功率、优化后的负荷群以及优化后的储能装置的充放电功率之间的功率流交换，并通过联络线将风电功率输送至电网，从而实现风电功率的平滑输出。
本发明一种考虑需求侧响应的风电功率分层平滑方法，是应用于由风电机组、储能装置和参与需求侧响应的负荷群通过联络线输出至电网而组成的风力发电系统中；所述负荷群包括：不可控负荷和可控负荷中的可平移负荷；其特点是，所述风电功率分层平滑方法是按如下步骤进行：
步骤1、根据所述风电机组的历史出力数据预测未来出力数据，从而利用式(1)-式(3)建立分时电价数学模型：
Pmax(t)＝Ps(t)(1+K) (1)
Pmin(t)＝Ps(t)(1-K) (2)
式(1)中，Ps(t)为t时刻联络线输出至电网的目标输出功率，且为定值；K为基准变化率；Pmax(t)为t时刻所述风电机组的输出功率的截止上限；Pmin(t)为t时刻所述风电机组的输出功率的截止下限；Pwind(t)为t时刻所述风电机组的预测输出功率，Creal(t)为t时刻的分时电价；CHigh、CNomal、CLow分别为分时电价中的高中低电价；
步骤2、获取所述可平移负荷的负荷参数信息并对可平移负荷按时间尺度的降序顺序进行分层处理，获得N层可平移负荷曲线；记为{L1,L2,…,Ln,…,LN}；Ln表示第n层可平移负荷曲线；1≤n≤N；同时，获取所述不可控负荷的负荷参数信息并得到不可控负荷曲线，记为L*；
步骤3、根据分时电价和所述未来风电机组出力数据，建立如式(4)和式(5)所示的储能装置的充放电状态数学模型：
式(4)中，S(t)为t时刻储能装置的充放电状态；S(t)＝1表示储能装置放电；S(t)＝0表示储能装置不工作；S(t)＝-1表示储能装置充电；
式(5)中，SOC(t)为t时刻储能装置的电荷状态，PBESS(t)为t时刻储能装置的充放电功率，PBESS(t)为正表示放电，PBESS(t)为负表示充电；SBESS为储能装置的总容量；Δt为充放电时间；
步骤4、利用式(6)获得N层可平移负荷调度的目标曲线；
式(6)中，Pn(t)为第n层可平移负荷调度的目标曲线，Li′(t)为优化后的n-1层可平移负荷曲线中的第i层可平移负荷曲线；
步骤5、建立负荷分层优化风电功率平滑模型：
步骤5.1、利用式(7)建立可平移负荷的N层经济性目标函数：
式(7)中，f1-n为第n层可平移负荷的经济性目标函数；T为时间间隔；
步骤5.2、利用式(8)建立可平移负荷的N层波动性目标函数：
式(8)中，f2-n为第n层可平移负荷的波动性目标函数；
步骤5.3、利用式(9)建立可平移负荷的N层平滑性目标函数：
步骤5.4、利用式(10)和式(11)设定可平移负荷优化时间范围约束：
α≤ton≤α+x-η (10)
toff＝ton+η (11)
式(10)和式(11)中，ton为可平移负荷运行的起始时间；toff为可平移负荷运行的结束时间；η为可平移负荷运行的时间；
步骤6、利用式(12)求取所述负荷分层优化风电功率平滑模型的适应度函数：
式(12)中，Fn为第n层可平移负荷的适应度函数；Ns为目标函数的个数；A为第n层可平移负荷的经济性目标函数f1-n的权重系数，B为第n层可平移负荷的波动性目标函数f2-n的权重系数，C为第n层可平移负荷的平滑性目标函数f3-n的权重系数；
步骤7、利用PSO算法对所述负荷分层优化风电功率平滑模型进行求解，获得优化后的N层可平移负荷曲线；记为{L′1,L′2,…,L′n,…,L′N}；L′n表示优化后的第n层可平移负荷曲线；
步骤8、建立所述储能装置的充放电状态优化模型：
步骤8.1、利用式(13)获得储能装置t时刻的目标曲线Paim-BESS(t)：
步骤8.2、利用式(14)获得储能装置的波动性目标函数fBESS：
步骤8.3、储能装置的充放电约束条件：
利用式(15)设定储能装置的荷电状态SOC(t)的约束条件：
SOCmin≤SOC(t)≤SOCmax (15)
式(15)中，SOCmin和SOCmax分别为SOC(t)的上下限；
利用式(16)设定储能装置的充放电功率约束条件：
PBESS(t)≤PBESS-max (16)
式(16)中，PBESS-max表示储能装置的最大充放电功率；
利用式(17)设定储能装置的始末能量约束条件：
式(17)中，SOC0为储能装置的初始SOC状态；SOC′0为储能装置的结束SOC状态；
步骤9、利用PSO算法对所述储能装置的充放电状态优化模型进行求解，获得优化后的储能装置的充放电状态；从而将所述风电机组的输出功率、优化后的负荷群、优化后的储能装置的充放电功率通过联络线输出至电网，以实现风电功率的平滑输出。
与已有技术相比，本发明的有益效果体现在：
1、本发明提出的考虑需求侧响应的风电功率分层平滑系统，由风电机组、储能装置、参与需求侧响应的负荷群、风电预测模块、负荷调度模块、负荷优化模块组成，充分利用了电网中各种资源，实现了各种资源直接的良性互动，降低联络线上风电功率输出的波动性，使风电功率输出更加平滑，维持了电力系统的稳定运行，并缓解增长的电力需求与有限的电力资源之间的矛盾，平衡了电力供求关系。
2、本发明建立了风电分时电价数学模型，将电价作为负荷分层优化风电功率平滑模型和储能装置充放电状态优化模型的重要指标，充分考虑了用户的利益，在提高用户侧经济性的同时，平滑了风电输出功率，降低了联络线上风电功率输出的波动性。
3、本发明采用了多时间尺度可平移负荷分层优化方法，建立负荷分层优化风电功率平滑模型；根据各层可平移负荷的特性，侧重不同的目标函数，设置不同的适应度函数，利用PSO粒子群优化算法进行求解，从而更加有效、全面的反映可平移负荷的实际情况，达到了最佳的优化结果。利用负荷平滑风电功率，可以减少储能装置的容量配置，延长了储能装置中蓄电池的使用寿命，从而降低了风电平滑系统中储能装置的成本。
4、本发明建立了储能装置充放电状态优化模型，在负荷分层优化后，对储能装置进行了优化，进一步降低供电侧联络线上风电功率输出的波动性，并解决了可平移负荷优化调度中负荷随机性带来的系统内可平移负荷不足的问题，从而使得风电功率分层平滑系统达到最佳的优化效果。
附图说明
图1为本发明所涉及的风电功率平滑系统示意图；
图2为本发明所涉及的多时间尺度可平移负荷示意图。
具体实施方式
本实施例中，一种考虑需求侧响应的风电功率分层平滑系统，包括：风电机组、储能装置、参与需求侧响应的负荷群、风电预测模块、负荷调度模块、负荷优化模块；负荷群包括：不可控负荷和可控负荷中的可平移负荷；
风电预测模块获取风电机组的历史出力数据并预测未来出力数据，从而建立分时电价数学模型；
负荷调度模块获取可平移负荷的负荷参数信息，并按时间尺度的降序顺序对可平移负荷进行分层处理，获得N层可平移负荷曲线；记为{L1,L2,…,Ln,…,LN}；Ln表示第n层可平移负荷曲线；1≤n≤N；同时，获取不可控负荷的负荷参数信息并得到不可控负荷曲线，记为L*；
负荷调度模块从风电预测模块获取未来风电机组出力数据和分时电价，并建立储能装置的充放电状态数学模型；
负荷调度模块根据未来风电机组出力数据和不可控负荷曲线L*获得N层可平移负荷调度的目标曲线；
负荷优化模块根据不可控负荷曲线L*、分时电价数学模型和N层可平移负荷调度的目标曲线，建立可平移负荷的N层经济性目标函数；
负荷优化模块根据N层可平移负荷曲线、分时电价数学模型和N层可平移负荷调度的目标曲线，建立可平移负荷的N层波动性目标函数；
负荷优化模块根据N层可平移负荷曲线和N层可平移负荷调度的目标曲线，建立可平移负荷的N层平滑性目标函数；
负荷优化模块建立由N层经济性目标函数、N层波动性目标函数和N层平滑性目标函数，以及所设定的可平移负荷的时间范围约束组成的负荷分层优化风电功率平滑模型；
负荷优化模块求取负荷分层优化风电功率平滑模型的适应度函数，并利用PSO算法对负荷分层优化风电功率平滑模型进行求解，获得优化后的N层可平移负荷曲线；从而使得负荷群得到优化；
负荷调度模块根据风电机组未来出力数据、不可控负荷曲线L*和优化后的N层可平移负荷曲线获得储能装置的目标曲线；
负荷调度模块根据储能装置充放电状态数学模型、储能装置的目标曲线和分时电价数学模型获得储能装置的波动性目标函数；
负荷调度模块建立由储能装置的目标曲线、储能装置的波动性目标函数以及所设定的充放电约束条件组成的储能装置的充放电状态优化模型；
负荷优化模块利用PSO算法对储能装置的充放电状态优化模型进行求解，获得优化后的储能装置的充放电状态，从而使得储能装置得到优化；
风电功率分层平滑系统通过联络线实现风电机组的输出功率、优化后的负荷群以及优化后的储能装置的充放电功率之间的功率流交换，并通过联络线将风电功率输送至电网，从而实现了源-荷-储的互动使负荷有效“跟踪”风电机组出力，使得风电功率能更平滑的输出。
本实施例中，一种考虑需求侧响应的风电功率分层平滑方法，是应用于由风电机组、储能装置和参与需求侧响应的负荷群通过联络线输出至电网而组成的风力发电系统中；负荷群包括：不可控负荷和可控负荷中的可平移负荷；如图1所示，风电机组通过联络线向电网中不参与需求响应的负荷群(即用电侧)送电，同时将风电数据发送给风电预测模块；实际情况下进行优化调度的风电机组出力数据应该为实际值，当预测值与实际值很接近的情况下，可将预测值替换实际值进行优化调度，风电预测模块接收风电机组出力的历史数据进行出力预测，将预测的风电数据发送到负荷调度模块；负荷优化模块接收负荷模块的优化数据，对负荷以及储能装置进行优化；负荷调度模块接收参与需求响应的负荷信息以及储能状态信息，按照负荷优化模块优化后的信息对负荷以及储能装置出力统一进行调度；具体的说，风电功率分层平滑方法是按如下步骤进行：
步骤1、根据风电机组的历史出力数据预测未来出力数据，从而利用式(1)-式(3)建立分时电价数学模型：
需求侧响应常通过价格信号引导用户改变自身的用电方式，减少或者推移某时段的用电负荷，促进供需两侧优化平衡。当风电机组预测输出功率波动较大，且超过截止上限功率值时，其电价较低，鼓励用电侧用电，当低于截止下限功率值时，其电价较高，则不鼓励用电侧在这个时刻用电。其分时电价模型可表示为：
Pmax(t)＝Ps(t)(1+K) (1)
Pmin(t)＝Ps(t)(1-K) (2)
式(1)中，Ps(t)为t时刻联络线输出至电网的目标输出功率，且为定值，可由风电机组出力的平均功率与负荷群的平均功率的差值求出；K为基准变化率，通常把风电功率的波动控制在20％以内；Pmax(t)为t时刻风电机组的输出功率的截止上限；Pmin(t)为t时刻风电机组的输出功率的截止下限；Pwind(t)为t时刻风电机组的预测输出功率，Creal(t)为t时刻的分时电价；CHigh、CNomal、CLow分别为分时电价中的高中低电价；
步骤2、获取可平移负荷的负荷参数信息并对可平移负荷按时间尺度的降序顺序进行分层处理，获得N层可平移负荷曲线；记为{L1,L2,…,Ln,…,LN}；Ln表示第n层可平移负荷曲线；1≤n≤N；同时，获取不可控负荷的负荷参数信息并得到不可控负荷曲线，记为L*；
可平移负荷表示为Lshift-xh(α,β,η)。其中x表示可平移负荷的时间尺度；α，β分别表示可平移负荷运行的起始时间与终止时间；η表示可平移负荷运行的时间。根据可平移负荷的可平移时间尺度的长短特性，将负荷定义为多时间尺度的可平移负荷，如图2所示，本实例将可平移的时间尺度分为24h，9h以及3h，即Lshift-24h(α,β,η)，Lshift-9h(α,β,η)，Lshift-3h(α,β,η)。可平移时间尺度为24h的负荷一般为工业负荷，如：钢铁，工地施工，有色金属加工。这些负荷，可以调动的时间跨度大，负荷的量大，对电的价格十分敏感。可平移时间尺度为9h的负荷如：电动汽车，洗衣机。这类负荷可以调度的时间跨度也比较大，数量也较多，对电价也比较敏感，可以很好的平滑一段时间内风电功率的波动。可平移时间尺度为3h的负荷如：热水器，空调。这类负荷由于平移调度的时间跨度较小，数量较少，所以很难用价格引导来大尺度的调度，但其精确性很好，能很好的弥补短期内风电功率的波动。
步骤3、根据分时电价和未来风电机组出力数据，建立如式(4)和式(5)所示的储能装置的充放电状态数学模型，实现了源-储互动：
储能频繁切换充放电状态会大大缩短蓄电池的使用寿命，提高储能系统的成本，因此有必要确定储能系统的充放电状态，防止其频繁的充放电。根据启发式规则的思想，当电价低于平均电价时储能充电，当电价高于平均电价时储能放电，可以将储能装置的充放电状态表示为：
式(4)中，S(t)为t时刻储能装置的充放电状态；S(t)＝1表示储能装置放电；S(t)＝0表示储能装置不工作；S(t)＝-1表示储能装置充电；
式(5)中，SOC(t)为t时刻储能装置的电荷状态，PBESS(t)为t时刻储能装置的充放电功率，PBESS(t)为正表示放电，PBESS(t)为负表示充电；SBESS为储能装置的总容量；Δt为充放电时间；
步骤4、利用式(6)获得N层可平移负荷调度的目标曲线；
式(6)中，Pn(t)为第n层可平移负荷调度的目标曲线，Li′(t)为优化后的n-1层可平移负荷曲线中的第i层可平移负荷曲线；
步骤5、建立负荷分层优化风电功率平滑模型：
步骤5.1、利用式(7)建立可平移负荷的N层经济性目标函数：
可平移负荷的优化调度与用户的利益息息相关，所以应从用户的角度考虑，提高用户参与需求侧响应的经济性。由式(3)可知，分时电价由风电机组出力的波动性决定，因此经济性目标函数从一定程度上反应风电机组出力的波动性，实现了源-荷的互动。通过分时电价来调度可平移负荷负荷的经济性目标函数如下：
式(7)中，f1-n为第n层可平移负荷的经济性目标函数；T为时间间隔；
步骤5.2、利用式(8)建立可平移负荷的N层波动性目标函数，波动性目标函数可以使可平移负荷曲线贴近目标曲线：
式(8)中，f2-n为第n层可平移负荷的波动性目标函数；
步骤5.3、利用式(9)建立可平移负荷的N层平滑性目标函数：
考虑到可平移负荷的优化调度后会改变原有负荷的曲线，使风电输出功率曲线变得不平滑，降低电能的质量，所以有必要对输出功率的平滑性进行优化。
通过相邻两个时间段的输出功率的起伏，可以反应整个输出功率曲线的平滑性。可平移负荷的N层平滑性目标函数可以表示为：
步骤5.4、利用式(10)和式(11)设定可平移负荷优化时间范围约束：
α≤ton≤α+x-η (10)
toff＝ton+η (11)
式(10)和式(11)中，ton为可平移负荷运行的起始时间；toff为可平移负荷运行的结束时间；η为可平移负荷运行的时间；
步骤6、利用式(12)求取负荷分层优化风电功率平滑模型的适应度函数：
考虑到本实例为多目标求解问题，考虑到多时间尺度可平移负荷每一层负荷的特性不同，所对应的各层目标函数的权重是不同的，因此采用设置权重系数的方法，将多目标加权求和得到适应度值F。第1层可平移负荷平移跨度大，对电价十分敏感，因此该层的可平移负荷的经济性目标函数的权重系数应较大；此后每层平移跨度逐渐减小，电价对其约束逐渐减小，因此经济性目标函数的权重系数逐渐减小，波动性目标函数的权重系数逐渐增大；各层可平移负荷曲线都需要进行平滑性的优化，因此各层平滑性目标函数的权重系数不变。
式(12)中，Fn为第n层可平移负荷的适应度函数；Ns为目标函数的个数；A为第n层可平移负荷的经济性目标函数f1-n的权重系数，B为第n层可平移负荷的波动性目标函数f2-n的权重系数，C为第n层可平移负荷的平滑性目标函数f3-n的权重系数；本实例中，A＝3，B＝4，C＝1。
步骤7、利用PSO算法对负荷分层优化风电功率平滑模型进行求解，获得优化后的N层可平移负荷曲线；记为{L′1,L′2,…,L′n,…,L′N}；L′n表示优化后的第n层可平移负荷曲线；
步骤8、建立储能装置的充放电状态优化模型，储能装置可以弥补优化调度中负荷随机性带来的系统内可平移负荷不足的问题：
步骤8.1、利用式(13)获得储能装置t时刻的目标曲线Paim-BESS(t)：
步骤8.2、利用式(14)获得储能装置的波动性目标函数fBESS：
步骤8.3、储能装置的充放电约束条件：
利用式(15)设定储能装置的荷电状态SOC(t)的约束条件：
SOCmin≤SOC(t)≤SOCmax (15)
式(15)中，SOCmin和SOCmax分别为SOC(t)的上下限；
利用式(16)设定储能装置的充放电功率约束条件：
PBESS(t)≤PBESS-max (16)
式(16)中，PBESS-max表示储能装置的最大充放电功率；
利用式(17)设定储能装置的始末能量约束条件：
式(17)中，SOC0为储能装置的初始SOC状态；SOC′0为储能装置的结束SOC状态；
步骤9、利用PSO算法对储能装置的充放电状态优化模型进行求解，获得优化后的储能装置的充放电状态；从而将风电机组的输出功率、优化后的负荷群、优化后的储能装置的充放电功率通过联络线输出至电网，以实现风电功率的平滑输出。
步骤10、本实例按如下方法求解负荷分层优化风电功率平滑模型以及储能装置的充放电状态优化模型：
步骤10.1、负荷调度模块从风电预测模块获取预测风电机组出力信息，读取用户设定的负荷参数信息。
步骤10.2、负荷调度模块设定分时电价；将可平移负荷按时间尺度分为N层，令n＝1；根据式(4)设定储能装置的充放电状态。
步骤10.3、根据式(6)求出第n层可平移负荷的目标曲线Pn(t)。
步骤10.4、初始化粒子种群：将第n层负荷作为一个粒子，设置m个粒子，粒子的个体极值pbest和全局极值gbest均设为粒子的初始位置。
步骤10.5、根据式(18)和式(19)更新权重因子w和学习因子c1、c2：
式(18)中，wmax、wmin为惯性权重因子的最大和最小值，取wmax＝1，wmin＝0.5；令迭代次数为Niter。
式(19)中，c1f和c1i是c1的最终值和初始值，c2f和c2i是c2的最终值和初始值。
步骤10.6、对目标函数Fn进行计算，得到当前粒子的适应度。更新局部最优适应度，更新体局部最优向量；更新全局最优适应度，更新全局最优向量。
步骤10.7、在考虑约束条件的情况下，利用式(20)更新每个粒子的位置和速度。
式(20)中，i＝1，2，…，m，m为种群的规模，表示粒子在第k次迭代中第i个粒子的d维的速度；表示粒子在第k次迭代中第i个粒子的d维的位置；ω表示惯性权重；c1、c2表示学习因子，表示粒子i在第k次迭代中d维的个体极值；表示整个粒子群在第k次迭代中d维的全局极值；为(0,1)区间分布的随机数。
步骤10.8、若迭代次数Niter到达最大值Niter-max，则停止搜索，输出第n层可平移负荷平移优化后的负荷L′n，否则返回⑤继续迭代计算。若n等于N则进入⑨，否则n＝n+1，返回③继续下一层优化。
步骤10.9、根据式(13)求得储能装置的目标曲线，以式(14)为适应度函数，按步骤④到⑧的方法求解储能装置的最优运行方式，输出PBESS(t)储能系统的充放电功率。
步骤10.10、负荷调度模块根据优化后的可平移负荷{L′1,L′2,…,L′n,…,L′N}以及优化后的储能装置的充放电状态，调度参与需求响应的负荷群以及控制储能运行方式。
Claims (2)
1.一种考虑需求侧响应的风电功率分层平滑系统，其特征是包括：风电机组、储能装置、参与需求侧响应的负荷群、风电预测模块、负荷调度模块、负荷优化模块；所述负荷群包括：不可控负荷和可控负荷中的可平移负荷；
所述风电预测模块获取所述风电机组的历史出力数据并预测未来出力数据，从而建立分时电价数学模型；
所述负荷调度模块获取所述可平移负荷的负荷参数信息，并按时间尺度的降序顺序对可平移负荷进行分层处理，获得N层可平移负荷曲线；记为{L1,L2,…,Ln,…,LN}；Ln表示第n层可平移负荷曲线；1≤n≤N；同时，获取所述不可控负荷的负荷参数信息并得到不可控负荷曲线，记为L*；
所述负荷调度模块从所述风电预测模块获取所述风电机组未来出力数据和分时电价，并建立所述储能装置的充放电状态数学模型；
所述负荷调度模块根据所述风电机组未来出力数据和不可控负荷曲线L*获得N层可平移负荷调度的目标曲线；
所述负荷优化模块根据所述不可控负荷曲线L*、分时电价数学模型和N层可平移负荷调度的目标曲线，建立可平移负荷的N层经济性目标函数；
所述负荷优化模块根据N层可平移负荷曲线、分时电价数学模型和N层可平移负荷调度的目标曲线，建立可平移负荷的N层波动性目标函数；
所述负荷优化模块根据N层可平移负荷曲线和N层可平移负荷调度的目标曲线，建立可平移负荷的N层平滑性目标函数；
所述负荷优化模块建立由所述N层经济性目标函数、N层波动性目标函数和N层平滑性目标函数，以及所设定的可平移负荷的时间范围约束组成的负荷分层优化风电功率平滑模型；
所述负荷优化模块求取所述负荷分层优化风电功率平滑模型的适应度函数，并利用PSO算法对所述负荷分层优化风电功率平滑模型进行求解，获得优化后的N层可平移负荷曲线；从而使得所述负荷群得到优化；
所述负荷调度模块根据所述风电机组未来出力数据、不可控负荷曲线L*和优化后的N层可平移负荷曲线获得储能装置的目标曲线；
所述负荷调度模块根据所述储能装置充放电状态数学模型、储能装置的目标曲线和分时电价数学模型获得储能装置的波动性目标函数；
所述负荷调度模块建立由所述储能装置的目标曲线、储能装置的波动性目标函数以及所设定的充放电约束条件组成的所述储能装置的充放电状态优化模型；
所述负荷优化模块利用PSO算法对所述储能装置的充放电状态优化模型进行求解，获得优化后的储能装置的充放电状态，从而使得所述储能装置得到优化；
所述风电功率分层平滑系统通过联络线实现所述风电机组的输出功率、优化后的负荷群以及优化后的储能装置的充放电功率之间的功率流交换，并通过联络线将风电功率输送至电网，从而实现风电功率的平滑输出。
2.一种考虑需求侧响应的风电功率分层平滑方法，是应用于由风电机组、储能装置和参与需求侧响应的负荷群通过联络线输出至电网而组成的风力发电系统中；所述负荷群包括：不可控负荷和可控负荷中的可平移负荷；其特征是，所述风电功率分层平滑方法是按如下步骤进行：
步骤1、根据所述风电机组的历史出力数据预测未来出力数据，从而利用式(1)-式(3)建立分时电价数学模型：
Pmax(t)＝Ps(t)(1+K) (1)
Pmin(t)＝Ps(t)(1-K) (2)
式(1)中，Ps(t)为t时刻联络线输出至电网的目标输出功率，且为定值；K为基准变化率；Pmax(t)为t时刻所述风电机组的输出功率的截止上限；Pmin(t)为t时刻所述风电机组的输出功率的截止下限；Pwind(t)为t时刻所述风电机组的预测输出功率，Creal(t)为t时刻的分时电价；CHigh、CNomal、CLow分别为分时电价中的高中低电价；
步骤2、获取所述可平移负荷的负荷参数信息并对可平移负荷按时间尺度的降序顺序进行分层处理，获得N层可平移负荷曲线；记为{L1,L2,…,Ln,…,LN}；Ln表示第n层可平移负荷曲线；1≤n≤N；同时，获取所述不可控负荷的负荷参数信息并得到不可控负荷曲线，记为L*；
步骤3、根据分时电价和所述风电机组未来出力数据，建立如式(4)和式(5)所示的储能装置的充放电状态数学模型：
式(4)中，S(t)为t时刻储能装置的充放电状态；S(t)＝1表示储能装置放电；S(t)＝0表示储能装置不工作；S(t)＝-1表示储能装置充电；
式(5)中，SOC(t)为t时刻储能装置的电荷状态，PBESS(t)为t时刻储能装置的充放电功率，PBESS(t)为正表示放电，PBESS(t)为负表示充电；SBESS为储能装置的总容量；Δt为充放电时间；
步骤4、利用式(6)获得N层可平移负荷调度的目标曲线；
式(6)中，Pn(t)为第n层可平移负荷调度的目标曲线，L′i(t)为优化后的n-1层可平移负荷曲线中的第i层可平移负荷曲线；
步骤5、建立负荷分层优化风电功率平滑模型：
步骤5.1、利用式(7)建立可平移负荷的N层经济性目标函数：
式(7)中，f1-n为第n层可平移负荷的经济性目标函数；T为时间间隔；
步骤5.2、利用式(8)建立可平移负荷的N层波动性目标函数：
式(8)中，f2-n为第n层可平移负荷的波动性目标函数；
步骤5.3、利用式(9)建立可平移负荷的N层平滑性目标函数：
式(9)中，f3-n为第n层可平移负荷的平滑性目标函数；
步骤5.4、利用式(10)和式(11)设定可平移负荷优化时间范围约束：
α≤ton≤α+x-η (10)
toff＝ton+η (11)
式(10)和式(11)中，ton为可平移负荷运行的起始时间；toff为可平移负荷运行的结束时间；η为可平移负荷运行的时间；
步骤6、利用式(12)求取所述负荷分层优化风电功率平滑模型的适应度函数：
式(12)中，Fn为第n层可平移负荷的适应度函数；Ns为目标函数的个数；A为第n层可平移负荷的经济性目标函数f1-n的权重系数，B为第n层可平移负荷的波动性目标函数f2-n的权重系数，C为第n层可平移负荷的平滑性目标函数f3-n的权重系数；
步骤7、利用PSO算法对所述负荷分层优化风电功率平滑模型进行求解，获得优化后的N层可平移负荷曲线；记为{L′1,L′2,…,L′n,…,L′N}；L′n表示优化后的第n层可平移负荷曲线；
步骤8、建立所述储能装置的充放电状态优化模型：
步骤8.1、利用式(13)获得储能装置t时刻的目标曲线Paim-BESS(t)：
步骤8.2、利用式(14)获得储能装置的波动性目标函数fBESS：
步骤8.3、储能装置的充放电约束条件：
利用式(15)设定储能装置的荷电状态SOC(t)的约束条件：
SOCmin≤SOC(t)≤SOCmax (15)
式(15)中，SOCmin和SOCmax分别为SOC(t)的上下限；
利用式(16)设定储能装置的充放电功率约束条件：
PBESS(t)≤PBESS-max (16)
式(16)中，PBESS-max表示储能装置的最大充放电功率；
利用式(17)设定储能装置的始末能量约束条件：
式(17)中，SOC0为储能装置的初始SOC状态；SOC′0为储能装置的结束SOC状态；
步骤9、利用PSO算法对所述储能装置的充放电状态优化模型进行求解，获得优化后的储能装置的充放电状态；从而将所述风电机组的输出功率、优化后的负荷群、优化后的储能装置的充放电功率通过联络线输出至电网，以实现风电功率的平滑输出。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610241244.8A CN105896596B (zh) | 2016-04-15 | 2016-04-15 | 一种考虑需求侧响应的风电功率分层平滑系统及其方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610241244.8A CN105896596B (zh) | 2016-04-15 | 2016-04-15 | 一种考虑需求侧响应的风电功率分层平滑系统及其方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105896596A CN105896596A (zh) | 2016-08-24 |
| CN105896596B true CN105896596B (zh) | 2018-08-28 |
Family
ID=56705063
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610241244.8A Expired - Fee Related CN105896596B (zh) | 2016-04-15 | 2016-04-15 | 一种考虑需求侧响应的风电功率分层平滑系统及其方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105896596B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107665378B (zh) * | 2017-09-25 | 2021-05-28 | 燕山大学 | 基于需求侧竞价交易机制的电力系统多目标优化调度方法 |
| CN109687479B (zh) * | 2017-10-19 | 2020-12-25 | 中国南方电网有限责任公司 | 功率波动平抑方法、系统、存储介质和计算机设备 |
| CN108227488B (zh) * | 2017-12-22 | 2020-02-04 | 上海交通大学 | 基于滑模预测控制的超超临界火电机组协调控制方法 |
| CN110601260B (zh) * | 2019-09-11 | 2021-05-14 | 电子科技大学 | 一种限定联络线上功率波动的光-蓄系统容量优化方法 |
| CN110729768A (zh) * | 2019-10-09 | 2020-01-24 | 南方电网能源发展研究院有限责任公司 | 一种分布式电源出力特性的增量配电网分时配电方法 |
| CN110943477B (zh) * | 2019-11-19 | 2021-12-28 | 国网江苏省电力有限公司经济技术研究院 | 提升电动汽车优化充电对分布式电源消纳的方法和装置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102411729A (zh) * | 2011-11-04 | 2012-04-11 | 国电南京自动化股份有限公司 | 基于自适应线性逻辑网络的风电功率预测方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009065778A (ja) * | 2007-09-06 | 2009-03-26 | Tokyo Electric Power Co Inc:The | 電力系統監視制御システム |
-
2016
- 2016-04-15 CN CN201610241244.8A patent/CN105896596B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102411729A (zh) * | 2011-11-04 | 2012-04-11 | 国电南京自动化股份有限公司 | 基于自适应线性逻辑网络的风电功率预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105896596A (zh) | 2016-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105896596B (zh) | 一种考虑需求侧响应的风电功率分层平滑系统及其方法 | |
| CN106410861B (zh) | 一种基于可调度能力的微电网优化运行实时控制方法 | |
| CN105160451B (zh) | 一种含电动汽车的微电网多目标优化调度方法 | |
| Wu et al. | Dynamic economic dispatch of a microgrid: Mathematical models and solution algorithm | |
| CN106230020B (zh) | 一种微电网下考虑分布式电源消纳的电动汽车互动响应控制方法 | |
| CN106058855A (zh) | 协调储能与柔性负荷的主动配电网多目标优化调度方法 | |
| CN107633333A (zh) | 柔性变电站区域电网的源‑荷‑储调度优化方法及系统 | |
| CN111769600B (zh) | 一种基于灵活性裕度的电力系统源荷储协调滚动调度方法 | |
| CN109146201A (zh) | 基于合作博弈的充换储一体化电站微电网优化调度方法 | |
| CN105071389B (zh) | 计及源网荷互动的交直流混合微电网优化运行方法及装置 | |
| CN107069776B (zh) | 一种平滑微网联络线功率的储能前瞻分布式控制方法 | |
| CN109816171A (zh) | 一种基于模糊pid实时电价的电动汽车区域微网群双层分布式优化调度方法 | |
| CN107453381B (zh) | 基于两阶段交叉控制的电动汽车集群功率调节方法及系统 | |
| CN106451552B (zh) | 一种基于势博弈的微电网能量管理系统分布式优化算法 | |
| CN105515110B (zh) | 一种电动汽车有序充电实时控制系统 | |
| CN108875992A (zh) | 一种考虑需求响应的虚拟电厂日前优化调度方法 | |
| CN109657993A (zh) | 一种基于非合作博弈的能源局域网储能系统自动需求响应方法 | |
| CN112800658A (zh) | 一种考虑源储荷互动的主动配电网调度方法 | |
| CN110138006A (zh) | 考虑含有新能源电动汽车的多微电网协调优化调度方法 | |
| CN106953316A (zh) | 微电网变时间尺度优化调度方法 | |
| CN108233430A (zh) | 一种计及系统能源波动性的交直流混合微网优化方法 | |
| CN108197766A (zh) | 一种包含微电网群的主动配电网优化调度模型 | |
| CN110350512A (zh) | 一种智能园区新能源发电站调度优化方法及系统 | |
| CN109409615A (zh) | 基于非合作博弈的充换储一体化电站微电网优化调度方法 | |
| CN114155103A (zh) | 基于区块链合作博弈的能量共享联盟灵活性交易方法 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180828 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |