CN110224434A - 一种基于随机鲁棒优化的电力系统调度方法 - Google Patents
一种基于随机鲁棒优化的电力系统调度方法 Download PDFInfo
- Publication number
- CN110224434A CN110224434A CN201910519809.8A CN201910519809A CN110224434A CN 110224434 A CN110224434 A CN 110224434A CN 201910519809 A CN201910519809 A CN 201910519809A CN 110224434 A CN110224434 A CN 110224434A
- Authority
- CN
- China
- Prior art keywords
- power
- wind
- node
- powered electricity
- generator
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005457 optimization Methods 0.000 title claims abstract description 15
- 230000005611 electricity Effects 0.000 claims description 84
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000003827 upregulation Effects 0.000 claims description 9
- 230000001105 regulatory effect Effects 0.000 claims description 8
- 230000003828 downregulation Effects 0.000 claims description 6
- 238000005553 drilling Methods 0.000 claims description 6
- 238000010248 power generation Methods 0.000 claims description 6
- 108010074506 Transfer Factor Proteins 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims description 3
- 238000012731 temporal analysis Methods 0.000 claims description 3
- 238000000700 time series analysis Methods 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims 1
- 230000035945 sensitivity Effects 0.000 claims 1
- 238000005096 rolling process Methods 0.000 abstract description 2
- 238000006243 chemical reaction Methods 0.000 abstract 1
- 230000005684 electric field Effects 0.000 description 2
- 239000013641 positive control Substances 0.000 description 2
- 241000219000 Populus Species 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000005259 measurement 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/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
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
-
- 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/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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- 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
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Power Engineering (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明涉及一种基于随机鲁棒优化的电力系统调度方法,属于电力系统运行控制技术领域。该方法首先建立由目标函数和约束条件构成的电力系统调度模型;对约束条件进行转化,将原模型转化为一个非线性优化问题;对转化后的模型求解,得到发电机出力,并根据求解结果进行相应调度决策实现电力系统调度。本发明考虑了电力系统中可再生能源的波动,减小了校正控制的成本,保证了电力系统安全,适合应用于电力系统的滚动调度等场景之中。
Description
技术领域
本发明涉及一种基于随机鲁棒优化的电力系统调度方法,属于电力系统运行控制技术领域。
背景技术
电力系统的调度对电力系统运行具有重要的指导作用,负责制定发电机组的出力计划,在满足用电负荷和保证电力系统安全的前提下实现最优运行。近年来越来越多的风电接入电力系统,而风电的出力具有随机性,因此需要在电力系统的调度中引入考虑不确定性的优化算法。
目前常用的调度方法主要是基于直流潮流的有功调度,通过求解一个确定性的二次规划问题得到调度计划,其不足之处在于无法考虑新能源接入电力系统带来的不确定性下的系统安全,且没有考虑发电机的功率调节成本。例如,由陈庆攀,刘洋,许立雄,周博,杨得洋等分别的“含风电的电力系统鲁棒优化调度[J/OL].电测与仪表”,参见1-7[2019-06-14].http://kns.cnki.net/kcms/detail/23.1202.TH.20190318.1205.004.html。其中,涉及风电且以最低发电成本为目标的鲁棒调度方法,考虑机组的自动发电控制响应来应对风电出力波动,维持系统功率平衡,其不足之处在于未能使用随机调度模型考虑风电的概率特性。
发明内容
本发明的目的是提出一种基于随机鲁棒优化的电力系统调度方法,以克服已有技术的不足之处,利用风电的概率预测信息描述可再生能源的不确定性,通过求解非线性优化问题得到电力系统的调度方案。
本发明提出的基于随机鲁棒优化的电力系统调度方法,包括以下步骤:
(1)利用时间序列分析方法对电力系统的功率进行预测,分别得到在电力系统的所有节点处风电的预测功率上限和下限w,将电力系统的第i个节点处风电的预测功率上限记为第i个节点处风电的预测功率下限记为wi ,记电力系统中风电的实际功率为w,w为一个多维不确定性变量;
(2)建立基于随机鲁棒优化的电力系统调度模型,具体过程包括以下步骤:
(2-1)确定电力系统调度模型的目标函数:
电力系统调度的目标为总成本最小化,总成本是发电成本、弃风惩罚和后续的校正控制成本之和,其中:
发电成本函数的表达式如下:
上式中,Gi表示电力系统第i个节点处发电机的发电成本函数,PGi表示电力系统第i个节点处发电机的有功功率,是待求解量,表示接入发电机的电力系统节点,从系统所有发电机的接入位置获取,是已知量;
第i个节点处发电机的发电成本函数Gi(PGi)的表达式如下:
Gi(PGi)=aiPGi 2+biPGi+ci
其中ai,bi,ci代表第i个节点处发电机的发电成本参数,ai,bi,ci由发电厂根据发电成本给出,为已知量;
弃风惩罚的表达式如下:
上式中,wi,max,wi,min分别表示电力系统第i个节点处风电的允许出力最大值和允许出力最小值,均为待求解量,表示接入风电的电力系统节点集合,从系统所有风电的接入位置信息获取,是已知量;
第i个节点处风电的弃风惩罚函数Ci(wi,max,wi,min)的表达式如下:
其中,ki代表第i个节点处风电的弃风惩罚系数,ki是已知量,由电力系统控制中心的调度系统给出,取值范围为大于零的正数,wi代表第i个节点处风电的实际功率,为一个随机变量,p(wi)代表第i个节点处风电的实际功率wi的概率密度函数,p(wi)为已知量,由风电场的概率预测得到;
第i个节点处发电机的校正控制成本的表达式如下:
其中,URi,DRi分别代表第i个节点处发电机的上调控制成本和下调控制成本;
第i个节点处发电机的上调控制成本URi表达式如下:
其中,gi表示第i个节点处发电机的上调成本系数,由发电厂根据发电机调节能力给出,是已知量,αi表示第i个节点处发电机的自动发电控制调节系数,由发电厂根据发电机调节能力给出,是已知量,u0表示所有风电的有功功率设定值总和,u0是待求解量,umin表示所有风电的允许出力最小值总和,umin是待求解量,u表示所有风电的实际功率总和,为一随机变量,p(u)代表所有风电的实际功率总和u的概率密度函数,p(u)是已知量,由风电场的概率预测得到;
第i个节点处发电机的下调控制成本DRi表达式如下:
其中,hi表示第i个节点处发电机的下调成本系数,由发电厂根据发电机调节能力给出,hi是已知量,umax表示所有风电的允许出力最大值总和,是待求解量,表示所有风电的预测功率上限总和,通过对步骤(1)中电力系统的所有节点处风电的预测功率上限求和得到;
(2-2)确定电力系统调度模型的约束条件,包括:
(2-2-1)发电机和风电的功率约束:
PGi,min≤PGi-αi(umax-u0)
PGi-αi(umin-u0)≤PGi,max
wi,min≤w
0≤wi,min≤wi,0≤wi,max
其中,wi,0表示第i个节点处风电的功率设定值,是待求解量,PGi,min,PGi,max分别表示第i个节点处发电机的功率最小值和功率最大,是已知量,由发电厂根据发电机的发电能力给出;
(2-2-2)电力系统的功率平衡方程约束:
其中,表示接入负载的电力系统节点集合,PDi表示电力系统第i个节点处负载的有功功率,是已知量,由电力系统控制中心根据负荷预测给出;
(2-2-3)正向支路功率容量约束和反向支路功率容量约束:
其中,分别表示电力系统中编号为L的支路相对于电力系统节点i和节点j的功率分布转移因子,为电力系统中编号为L的支路相对于节点i的有功功率的灵敏度:
其中,表示编号为L的支路的最大允许功率,是已知量,由电力系统控制中心根据支路的设计参数给出;
(3)将上述步骤(2)的目标函数和约束条件构成的随机鲁棒优化的电力系统调度模型改写成如下的数学形式:
Ax≤b
其中,代表对x取最小值,x表示由电力系统中所有发电机的有功功率、所有风电的功率设定值、所有风电的允许出力最小值和所有风电的允许出力最大值组成的向量,f(x)表示电力系统调度目标函数,将上述步骤(2)中的约束条件写成矩阵不等式形式,矩阵不等式的左侧为与x对应的矩阵A,矩阵不等式的右侧为常数矩阵b;
(4)利用序列线性化方法,求解上述步骤(3)的矩阵不等式,得到最优解,包括如下步骤:
(4-1)将f(x)拆分成两部分,包括f1(x)和f2(x),其中f1(x)代表步骤(2-1)中的发电成本,f1(x)为二次函数,f2(x)代表步骤(2-1)中的弃风惩罚和校正控制成本之和,f2(x)为非线性函数;
(4-2)采用二次规划求解方法,通过下式求解得到最优解x0:
Ax≤b
(4-3)采用二次规划求解方法,通过下式求解得到最优解为x1:
Ax≤b
其中,表示f2(x)函数在x0的导数向量;
(4-4)设定一个迭代允许误差ε,对步骤(4-2)和步骤(4-3)的计算结果进行判断,若abs(f(x1)-f(x0))<ε,则判定x1为电力系统调度最优解,停止迭代,若abs(f(x1)-f(x0))≥ε,则将x1的值代替x0,返回步骤(4-3),其中abs表示取绝对值;
(5)将与步骤(4)得到的电力系统调度最优解中相对应的发电机有功功率作为电力系统各发电机的发电量,制定发电机的发电计划,同时将与步骤(4)得到的电力系统调度最优解中相对应的所有风电的允许出力最小值和允许出力最大值下发给风电场,制定风电场的发电计划,完成电力系统的调度。
本发明提出的基于随机鲁棒优化的电力系统调度方法,其优点是:
本发明的基于随机鲁棒优化的电力系统调度方法,其目标函数采用随机模型考虑风电的概率特性,约束条件采用鲁棒模型,利用风电的概率预测信息描述可再生能源的不确定性,通过求解非线性优化问题得到调度方案,通过鲁棒优化保证了电力系统在风电波动下的安全性。本发明方法考虑了电力系统中可再生能源的波动,减小了校正控制的成本,保证了电力系统安全,适合应用于电力系统的滚动调度等场景之中。
具体实施方式
本发明提出的基于随机鲁棒优化的电力系统调度方法,包括以下步骤:
(1)利用时间序列分析方法对电力系统的功率进行预测,分别得到在电力系统的所有节点处风电的预测功率上限和下限w,将电力系统的第i个节点处风电的预测功率上限记为第i个节点处风电的预测功率下限记为wi ,记电力系统中风电的实际功率为w,w为一个多维不确定性变量;
(2)建立基于随机鲁棒优化的电力系统调度模型,具体过程包括以下步骤:
(2-1)确定电力系统调度模型的目标函数:
电力系统调度的目标为总成本最小化,总成本是发电成本、弃风惩罚和后续的校正控制成本之和,其中:
发电成本函数的表达式如下:
上式中,Gi表示电力系统第i个节点处发电机的发电成本函数,PGi表示电力系统第i个节点处发电机的有功功率,是待求解量,表示接入发电机的电力系统节点,从系统所有发电机的接入位置获取,是已知量;
第i个节点处发电机的发电成本函数Gi(PGi)的表达式如下:
Gi(PGi)=aiPGi 2+biPGi+ci
其中ai,bi,ci代表第i个节点处发电机的发电成本参数,ai,bi,ci由发电厂根据发电成本给出,为已知量;
弃风惩罚的表达式如下:
上式中,wi,max,wi,min分别表示电力系统第i个节点处风电的允许出力(即风电的允许有功功率)最大值和允许出力最小值,均为待求解量,表示接入风电的电力系统节点集合,从系统所有风电的接入位置信息获取,是已知量;
第i个节点处风电的弃风惩罚函数Ci(wi,max,wi,min)的表达式如下:
其中,ki代表第i个节点处风电的弃风惩罚系数,ki是已知量,由电力系统控制中心的调度系统给出,取值范围为大于零的正数,wi代表第i个节点处风电的实际功率,为一个随机变量,p(wi)代表第i个节点处风电的实际功率wi的概率密度函数,p(wi)为已知量,由风电场的概率预测得到;
第i个节点处发电机的校正控制成本的表达式如下:
其中,URi,DRi分别代表第i个节点处发电机的上调控制成本和下调控制成本;
第i个节点处发电机的上调控制成本URi表达式如下:
其中,gi表示第i个节点处发电机的上调成本系数,由发电厂根据发电机调节能力给出,是已知量,αi表示第i个节点处发电机的自动发电控制调节系数,由发电厂根据发电机调节能力给出,是已知量,u0表示所有风电的有功功率设定值总和,u0是待求解量,umin表示所有风电的允许出力最小值总和,umin是待求解量,u表示所有风电的实际功率总和,为一随机变量,p(u)代表所有风电的实际功率总和u的概率密度函数,p(u)是已知量,由风电场的概率预测得到;
第i个节点处发电机的下调控制成本DRi表达式如下:
其中,hi表示第i个节点处发电机的下调成本系数,由发电厂根据发电机调节能力给出,hi是已知量,umax表示所有风电的允许出力最大值总和,是待求解量,表示所有风电的预测功率上限总和,通过对步骤(1)中电力系统的所有节点处风电的预测功率上限求和得到;
(2-2)确定电力系统调度模型的约束条件,包括:
(2-2-1)发电机和风电的功率约束:
PGi,min≤PGi-αi(umax-u0)
PGi-αi(umin-u0)≤PGi,max
wi,min≤w
0≤wi,min≤wi,0≤wi,max
其中,wi,0表示第i个节点处风电的功率设定值,是待求解量,PGi,min,PGi,max分别表示第i个节点处发电机的功率最小值和功率最大,是已知量,由发电厂根据发电机的发电能力给出;
(2-2-2)电力系统的功率平衡方程约束:
其中,表示接入负载的电力系统节点集合,PDi表示电力系统第i个节点处负载的有功功率,是已知量,由电力系统控制中心根据负荷预测给出;
(2-2-3)正向支路功率容量约束和反向支路功率容量约束:
其中,分别表示电力系统中编号为L的支路相对于电力系统节点i和节点j的功率分布转移因子,为电力系统中编号为L的支路相对于节点i的有功功率的灵敏度:
其中,表示编号为L的支路的最大允许功率,是已知量,由电力系统控制中心根据支路的设计参数给出;
(3)将上述步骤(2)的目标函数和约束条件构成的随机鲁棒优化的电力系统调度模型改写成如下的数学形式:
Ax≤b
其中,代表对x取最小值,x表示由电力系统中所有发电机的有功功率、所有风电的功率设定值、所有风电的允许出力最小值和所有风电的允许出力最大值组成的向量,f(x)表示电力系统调度目标函数,将上述步骤(2)中的约束条件写成矩阵不等式形式,矩阵不等式的左侧为与x对应的矩阵A,矩阵不等式的右侧为常数矩阵b;
(4)利用序列线性化方法,求解上述步骤(3)的矩阵不等式,得到最优解,包括如下步骤:
(4-1)将f(x)拆分成两部分,包括f1(x)和f2(x),其中f1(x)代表步骤(2-1)中的发电成本,f1(x)为二次函数,f2(x)代表步骤(2-1)中的弃风惩罚和校正控制成本之和,f2(x)为非线性函数;
(4-2)采用二次规划求解方法,通过下式求解得到最优解x0:
Ax≤b
(4-3)采用二次规划求解方法,通过下式求解得到最优解为x1:
Ax≤b
其中,表示f2(x)函数在x0的导数向量;
(4-4)设定一个迭代允许误差ε,对步骤(4-2)和步骤(4-3)的计算结果进行判断,若abs(f(x1)-f(x0))<ε,则判定x1为电力系统调度最优解,停止迭代,若abs(f(x1)-f(x0))≥ε,则将x1的值代替x0,返回步骤(4-3),其中abs表示取绝对值,本发明的一个实施例中,ε的取值为0.001;
(5)将与步骤(4)得到的电力系统调度最优解中相对应的发电机有功功率作为电力系统各发电机的发电量,制定发电机的发电计划,同时将与步骤(4)得到的电力系统调度最优解中相对应的所有风电的允许出力最小值和允许出力最大值下发给风电场,制定风电场的发电计划,完成电力系统的调度。
Claims (1)
1.一种基于随机鲁棒优化的电力系统调度方法,其特征在于包括以下步骤:
(1)利用时间序列分析方法对电力系统的功率进行预测,分别得到在电力系统的所有节点处风电的预测功率上限和下限w,将电力系统的第i个节点处风电的预测功率上限记为第i个节点处风电的预测功率下限记为w i,记电力系统中风电的实际功率为w,w为一个多维不确定性变量;
(2)建立基于随机鲁棒优化的电力系统调度模型,具体过程包括以下步骤:
(2-1)确定电力系统调度模型的目标函数:
电力系统调度的目标为总成本最小化,总成本是发电成本、弃风惩罚和后续的校正控制成本之和,其中:
发电成本函数的表达式如下:
上式中,Gi表示电力系统第i个节点处发电机的发电成本函数,PGi表示电力系统第i个节点处发电机的有功功率,是待求解量,表示接入发电机的电力系统节点,从系统所有发电机的接入位置获取,是已知量;
第i个节点处发电机的发电成本函数Gi(PGi)的表达式如下:
Gi(PGi)=aiPGi 2+biPGi+ci
其中ai,bi,ci代表第i个节点处发电机的发电成本参数,ai,bi,ci由发电厂根据发电成本给出,为已知量;
弃风惩罚的表达式如下:
上式中,wi,max,wi,min分别表示电力系统第i个节点处风电的允许出力最大值和允许出力最小值,均为待求解量,表示接入风电的电力系统节点集合,从系统所有风电的接入位置信息获取,是已知量;
第i个节点处风电的弃风惩罚函数Ci(wi,max,wi,min)的表达式如下:
其中,ki代表第i个节点处风电的弃风惩罚系数,ki是已知量,由电力系统控制中心的调度系统给出,wi代表第i个节点处风电的实际功率,为一个随机变量,p(wi)代表第i个节点处风电的实际功率wi的概率密度函数,p(wi)为已知量,由风电场的概率预测得到;
第i个节点处发电机的校正控制成本的表达式如下:
其中,URi,DRi分别代表第i个节点处发电机的上调控制成本和下调控制成本;
第i个节点处发电机的上调控制成本URi表达式如下:
其中,gi表示第i个节点处发电机的上调成本系数,由发电厂根据发电机调节能力给出,是已知量,αi表示第i个节点处发电机的自动发电控制调节系数,由发电厂根据发电机调节能力给出,是已知量,u0表示所有风电的有功功率设定值总和,u0是待求解量,umin表示所有风电的允许出力最小值总和,umin是待求解量,u表示所有风电的实际功率总和,为一随机变量,p(u)代表所有风电的实际功率总和u的概率密度函数,p(u)是已知量,由风电场的概率预测得到;
第i个节点处发电机的下调控制成本DRi表达式如下:
其中,hi表示第i个节点处发电机的下调成本系数,由发电厂根据发电机调节能力给出,hi是已知量,umax表示所有风电的允许出力最大值总和,是待求解量,表示所有风电的预测功率上限总和,通过对步骤(1)中电力系统的所有节点处风电的预测功率上限求和得到;
(2-2)确定电力系统调度模型的约束条件,包括:
(2-2-1)发电机和风电的功率约束:
PGi,min≤PGi-αi(umax-u0)
PGi-αi(umin-u0)≤PGi,max
wi,min≤w
0≤wi,min≤wi,0≤wi,max
其中,wi,0表示第i个节点处风电的功率设定值,是待求解量,PGi,min,PGi,max分别表示第i个节点处发电机的功率最小值和功率最大,是已知量,由发电厂根据发电机的发电能力给出;
(2-2-2)电力系统的功率平衡方程约束:
其中,表示接入负载的电力系统节点集合,PDi表示电力系统第i个节点处负载的有功功率,是已知量,由电力系统控制中心根据负荷预测给出;
(2-2-3)正向支路功率容量约束和反向支路功率容量约束:
其中,分别表示电力系统中编号为L的支路相对于电力系统节点i和节点j的功率分布转移因子,为电力系统中编号为L的支路相对于节点i的有功功率的灵敏度:
其中,表示编号为L的支路的最大允许功率,是已知量,由电力系统控制中心根据支路的设计参数给出;
(3)将上述步骤(2)的目标函数和约束条件构成的随机鲁棒优化的电力系统调度模型改写成如下的数学形式:
Ax≤b
其中,代表对x取最小值,x表示由电力系统中所有发电机的有功功率、所有风电的功率设定值、所有风电的允许出力最小值和所有风电的允许出力最大值组成的向量,f(x)表示电力系统调度目标函数,将上述步骤(2)中的约束条件写成矩阵不等式形式,矩阵不等式的左侧为与x对应的矩阵A,矩阵不等式的右侧为常数矩阵b;
(4)利用序列线性化方法,求解上述步骤(3)的矩阵不等式,得到最优解,包括如下步骤:
(4-1)将f(x)拆分成两部分,包括f1(x)和f2(x),其中f1(x)代表步骤(2-1)中的发电成本,f1(x)为二次函数,f2(x)代表步骤(2-1)中的弃风惩罚和校正控制成本之和,f2(x)为非线性函数;
(4-2)采用二次规划求解方法,通过下式求解得到最优解x0:
Ax≤b
(4-3)采用二次规划求解方法,通过下式求解得到最优解为x1:
Ax≤b
其中,表示f2(x)函数在x0的导数向量;
(4-4)设定一个迭代允许误差ε,对步骤(4-2)和步骤(4-3)的计算结果进行判断,若abs(f(x1)-f(x0))<ε,则判定x1为电力系统调度最优解,停止迭代,若abs(f(x1)-f(x0))≥ε,则将x1的值代替x0,返回步骤(4-3),其中abs表示取绝对值,;
(5)将与步骤(4)得到的电力系统调度最优解中相对应的发电机有功功率作为电力系统各发电机的发电量,制定发电机的发电计划,同时将与步骤(4)得到的电力系统调度最优解中相对应的所有风电的允许出力最小值和允许出力最大值下发给风电场,制定风电场的发电计划,完成电力系统的调度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910519809.8A CN110224434B (zh) | 2019-06-17 | 2019-06-17 | 一种基于随机鲁棒优化的电力系统调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910519809.8A CN110224434B (zh) | 2019-06-17 | 2019-06-17 | 一种基于随机鲁棒优化的电力系统调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110224434A true CN110224434A (zh) | 2019-09-10 |
CN110224434B CN110224434B (zh) | 2020-09-01 |
Family
ID=67817360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910519809.8A Active CN110224434B (zh) | 2019-06-17 | 2019-06-17 | 一种基于随机鲁棒优化的电力系统调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110224434B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110985294A (zh) * | 2019-12-10 | 2020-04-10 | 华北电力大学 | 一种鲁棒概率tube联合使用的随机模型预测控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130066791A1 (en) * | 2011-09-09 | 2013-03-14 | Kabushiki Kaisha Toshiba | Device and method for determining storage battery rental capacity |
CN103793612A (zh) * | 2014-02-18 | 2014-05-14 | 广西大学 | 一种适合于考虑风电随机特性的电力系统的电网规划方法 |
CN104242356A (zh) * | 2014-09-26 | 2014-12-24 | 国家电网公司 | 考虑风电场集电线故障的鲁棒区间风电调度方法及装置 |
CN108667010A (zh) * | 2018-05-04 | 2018-10-16 | 清华大学 | 一种基于分布鲁棒优化的配电网经济调度方法 |
-
2019
- 2019-06-17 CN CN201910519809.8A patent/CN110224434B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130066791A1 (en) * | 2011-09-09 | 2013-03-14 | Kabushiki Kaisha Toshiba | Device and method for determining storage battery rental capacity |
CN103793612A (zh) * | 2014-02-18 | 2014-05-14 | 广西大学 | 一种适合于考虑风电随机特性的电力系统的电网规划方法 |
CN104242356A (zh) * | 2014-09-26 | 2014-12-24 | 国家电网公司 | 考虑风电场集电线故障的鲁棒区间风电调度方法及装置 |
CN108667010A (zh) * | 2018-05-04 | 2018-10-16 | 清华大学 | 一种基于分布鲁棒优化的配电网经济调度方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110985294A (zh) * | 2019-12-10 | 2020-04-10 | 华北电力大学 | 一种鲁棒概率tube联合使用的随机模型预测控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110224434B (zh) | 2020-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ma et al. | Optimal allocation of hybrid energy storage systems for smoothing photovoltaic power fluctuations considering the active power curtailment of photovoltaic | |
Khatod et al. | Evolutionary programming based optimal placement of renewable distributed generators | |
Suchitra et al. | Optimal design of hybrid power generation system and its integration in the distribution network | |
CN106532778B (zh) | 一种计算分布式光伏并网最大准入容量的方法 | |
CN110970912A (zh) | 一种含储能的新能源电力系统运行模拟方法 | |
CN107910863A (zh) | 综合考虑光伏出力和负荷需求预测区间的配电网调度方法 | |
CN111092429A (zh) | 一种柔性互联配电网的优化调度方法、存储介质及处理器 | |
CN114707767B (zh) | 一种新能源电力系统低谷时段可调峰功率预测方法 | |
CN109840636A (zh) | 一种基于牛顿法的电力系统随机滚动调度方法 | |
Maheshwari et al. | Flow direction algorithm-based optimal power flow analysis in the presence of stochastic renewable energy sources | |
Hong et al. | Interactive multi-objective active power scheduling considering uncertain renewable energies using adaptive chaos clonal evolutionary programming | |
CN103904664B (zh) | 一种基于有效静态安全域的agc机组实时调度方法 | |
Maulik | Probabilistic power management of a grid-connected microgrid considering electric vehicles, demand response, smart transformers, and soft open points | |
CN112821470A (zh) | 基于小生境混沌粒子群算法的微网群优化调度策略 | |
CN109755959A (zh) | 基于风/光出力柯西分布的火电机组动态实时调度方法 | |
Arief et al. | ANALYTICAL METHOD FOR REACTIVE POWER COMPENSATORS ALLOCATION. | |
CN116388278A (zh) | 一种微电网群协同控制方法、装置、设备及介质 | |
CN108667010B (zh) | 一种基于分布鲁棒优化的配电网经济调度方法 | |
CN110224434A (zh) | 一种基于随机鲁棒优化的电力系统调度方法 | |
CN116526511B (zh) | 一种多源协同参与系统负荷频率控制的方法 | |
Haddi et al. | Improved optimal power flow for a power system incorporating wind power generation by using Grey Wolf Optimizer algorithm | |
Thorat et al. | Optimization of fuel cost incorporating with wind, solar PV and Electric vehicle energy sources using improved artificial bee colony algorithm | |
CN111753438B (zh) | 一种基于时序特性的配电网分布式电源及储能双层优化规划方法 | |
CN112861376B (zh) | 一种基于单元调度模型的评估方法及装置 | |
CN115764900A (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 |