CN109755959A - 基于风/光出力柯西分布的火电机组动态实时调度方法 - Google Patents
基于风/光出力柯西分布的火电机组动态实时调度方法 Download PDFInfo
- Publication number
- CN109755959A CN109755959A CN201811509293.0A CN201811509293A CN109755959A CN 109755959 A CN109755959 A CN 109755959A CN 201811509293 A CN201811509293 A CN 201811509293A CN 109755959 A CN109755959 A CN 109755959A
- Authority
- CN
- China
- Prior art keywords
- power
- output
- wind power
- scheduling
- photovoltaic power
- 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 37
- 230000014509 gene expression Effects 0.000 claims description 26
- 238000005315 distribution function Methods 0.000 claims description 22
- 230000001186 cumulative effect Effects 0.000 claims description 17
- 239000013598 vector Substances 0.000 claims description 15
- 230000009194 climbing Effects 0.000 claims description 10
- 239000000446 fuel Substances 0.000 claims description 9
- 238000010248 power generation Methods 0.000 claims description 4
- 101100272279 Beauveria bassiana Beas gene Proteins 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims 1
- 238000005457 optimization Methods 0.000 abstract description 5
- 230000010354 integration Effects 0.000 abstract description 2
- 230000005611 electricity Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
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
-
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/10—Photovoltaic [PV]
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种基于风/光出力柯西分布的火电机组动态实时调度方法,属于电力系统的运行技术领域。本通过分析风电出力的历史数据,利用统计或拟合软件进行联合柯西分布拟合。针对确定的电力系统参数,建立机会约束的随机动态实时调度模型;然后,利用柯西分布的数学性质,把原问题转化为容易求解的线性约束凸优化问题;最后求解调度模型,得到调度策略。本发明充分利用了柯西分布在风电/光伏电站出力的短期预测方面的优越性和优良的数学特性,有效提高了模型的求解效率,同时风险水平可调的机会约束模型消除了传统鲁棒经济调度的保守性,为决策者提供了更合理的调度依据。本发明可应用于包含大规模风电并网的电力系统有功实时经济调度中。
Description
技术领域
本发明涉及一种基于风/光出力柯西分布的火电机组动态实时调度方法,尤其涉及一种基于风电/光伏出力柯西分布的火电机组动态实时调度方法,属于电力系统的运行技术领域。
背景技术
开发利用风电资源、实现能源的可持续发展是我国能源发展战略的重大举措。随着风电大规模接入电网,其波动性和随机性给电力系统的有功调度带来了两个方面的难题。
一方面,精确、灵活的风电出力预测是实现安全、经济有功调度的基础,传统的预测方法包括给定出力上下限的区间描述法和简单的高斯概率密度函数描述法,虽然像贝塔分布、通用分布、混合高斯分布等模型也被用在了风电预测出力的拟合中,但是它们或者无法精确拟合风电预测出力,或者给有功调度模型的求解带来了极大的困难,因此一种准确、灵活的预测模型亟需得到应用。
另一方面,风电的波动性和随机性使得传统的确定性调度方法难以适用。鲁棒经济调度通常是一种可行的方案,然而由于鲁棒优化具有保守性,会给调度带来不必要的成本;机会约束的随机经济调度是兼顾系统运行风险和减少成本的有效建模策略,该方法把风险发生的概率限制在预先给定的置信水平下,通过目标函数值的最小化得到成本最低的调度策略。然而约束和目标函数中存在的随机变量使得机会约束优化问题的求解变得非常困难,现有的求解方法普遍存在计算量大的缺点,然而松弛的方法又使得求解结果不够精确,无法实现经济调度的高效性。
综上所述,计及风电出力随机性的动态经济调度的建模以及快速求解仍然是影响风电资源利用率的一大难题。
发明内容
本发明的目的是提出一种基于风/光出力柯西分布的火电机组动态实时调度方法,基于联合柯西分布,对风电短期出力进行精确拟合,以充分利用机会约束随机经济调度的优点,有效降低系统的风险,节约电网调度的成本。
本发明提出的基于风/光出力柯西分布的火电机组动态实时调度方法,包括以下步骤:
(1)建立电网中风电/光伏电站短期预测出力的多随机变量联合柯西分布模型,该柯西分布模型包括以下内容:
a、一个多风电/光伏电站短期预测出力的联合柯西分布的概率密度函数:
其中,PDF(·)表示随机变量的概率密度函数,概率密度函数通过对风电/光伏电站历史出力的拟合得到,t为调度时段,K为风电/光伏电站的数目,为调度时段t所有风电/光伏电站实际出力的列向量,上标T为矩阵转置,上标w表示该变量描述风电/光伏电站,表示t调度时段第k个风电/光伏电站的实际出力,μt=(μ1,t,μ2,t,...,μk,t)T表示t调度时段联合柯西分布概率密度函数的位置参数,μk,t表示t调度时段第k个风电/光伏电站出力的边缘柯西分布的位置参数,Σt表示t调度时段联合柯西分布概率密度函数的尺度参数;
b、形式为各风电/光伏电站出力线性组合的随机变量的概率描述,包括概率密度函数、累积分布函数以及累积分布函数的逆函数:
设a为k维列向量,那么随机变量服从位置参数和尺度参数为(aTμt,aTΣta)的一维柯西分布,其概率密度函数、累积分布函数以及累积分布函数的逆函数可以表达为以下形式:
其中,概率密度函数为:
累积分布函数为:
累积分布函数的逆函数为:
其中,tan为正切函数,arctan为反正切函数,F为分位数;
(2)建立一个基于风电/光伏电站出力柯西分布的随机动态实时调度模型,该随机动态实时调度模型由目标函数和约束条件组成,具体步骤如下:
(2-1)建立基于风电/光伏电站出力柯西分布的随机动态实时调度模型的目标函数f:
目标函数为使运行成本最小化,表达式如下:
其中,T、N和J分别表示调度时段t的数量、火电机组数量以及发电量自动控制机组的数量,t、i和j分别为调度时段、火电机组编号和发电量自动控制机组的编号,上标s表示该变量描述火电机组,上标“+”表示该变量描述正旋转备用,上标“-”表示该变量描述负旋转备用,表示第i台火电机组在t调度时段的计划出力,表示第j个发电量自动控制机组在t调度时段内计划出力,表示t调度时段所有风电/光伏电站实际出力的总和,和分别表示火电机组和发电量自动控制机组的燃料成本:
其中,ai,t,bi,t,ci,t分别为t时段火电机组i的燃料成本的二次项系数、一次项系数和常数项,aj,t,bj,t,cj,t分别为t时段发电量自动控制机组j的燃料成本的二次项系数、一次项系数和常数项,
E(·)表示随机变量的期望值,表示在t调度时段内,由于风电/光伏电站发电实际出力低于计划出力引起的正旋转备用的需求成本,即高估风电/光伏电站出力的惩罚,当风电/光伏电站的实际出力小于计划值时会调度发电量自动控制机组的正旋转备用以维持功率平衡,具体表达式如下:
其中,为正旋转备用的成本系数,αj为第j个发电量自动控制机组的功率分配系数,根据发电量自动控制机组的额定容量占发电量自动控制机组的总容量的比例确定,wt表示t调度时段所有风电/光伏电站的计划出力总和,表示t时段所有风电/光伏电站的实际出力总和,满足如下关系:
其中,为第j个发电量自动控制机组的实际出力,表示t调度时段第k个风电/光伏电站的计划出力,K表示风电/光伏电站的数量;
表示在t时段内由于风电/光伏电站实际出力超过计划出力而引起的负旋转备用的需求成本,即低估风电/光伏电站出力的惩罚成本,当风电/光伏电站的实际出力大于计划值时会调度发电量自动控制机组的负旋转备用以维持功率平衡,具体表达式如下:
其中,为负旋转备用的成本系数,为实际风电出力总和的概率密度函数,为所有风电/光伏电站在t调度时段出力的上界之和, 为t时段第k个风电/光伏电站的出力的上界;
根据上述步骤(1)的风电/光伏电站短期预测出力的多随机变量联合柯西分布模型,得到目标函数f的表达式中的后两项为:
其中, 为所有元素全为1的K维列向量,K为风电/光伏电站的数目,A,B,C均为常数具体表达式如下:
(2-2)上述基于风电/光伏电站出力柯西分布的随机动态实时调度模型的约束条件,包括:
(2-2-1)电网功率平衡约束,表达式如下:
其中,为t调度时段、风电/光伏电站和火电机组所在电网的第d个节点的负荷量,D表示负荷的总数,同时表示电网中节点的个数;
(2-2-2)电网机组出力的上、下限约束,包括:
其中,t=1,...,T,i=1,...,N,j=1,...,J,k=1,...,K,和分别为第j个发电量自动控制机组在t调度时段出力的上下界,和分别为第i个火电机组在t时段出力的上下界,δ为调度员设定的可接受的风险水平,为的累积分布函数的逆函数,
表示t时段所有风电/光伏电站的实际出力总和,
为所有元素全为1的K维列向量,K表示风电/光伏电站的数量;
(2-2-3)电网中火电机组和发电量自动控制机组的爬坡约束,表达式如下:
对i=1,2,...,N;j=1,2,...,J:
其中,和分别为t调度时段第i台火电机组向上、向下爬坡率,和分别表示t时段第j台发电量自动控制机组向上、向下爬坡率,ΔT表示相邻两个调度时段之间的调度间隔,β为调度员设定的可接受风险水平,表示wt,t-1的累积分布函数的逆函数, 的表达式如下:
(2-2-4)发电量自动控制机组的旋转备用约束,具体表达式如下:
对i=1,2,...,N,j=1,2,...,J,
其中,和分别表示t调度时段第j台发电量自动控制机组提供的正负旋转备用的数量,和分别表示t调度时段电网所需要的最小正、负旋转备用的数量,ε为调度员设定的可接受的风险水平;
(2-2-5)电网线路潮流约束,表达式如下:
对l=1,2,...,L:
其中,Gi,l为电网中第l条线路对第i台火力发电机组有功出力的转移分布因子,Gj,l为第l条线路对第j个发电量自动控制机组有功出力的转移分布因子,Gk,l为第l条线路对第k个风电/光伏电站有功出力的转移分布因子,Gd,l为第l条线路对第d个节点负荷功率的转移分布因子,各转移分布因子分别从电网调度中心获取,Ll,t为t调度时段第l条线路上的有功功率上限,η为电网线路上的有功功率超过线路有功功率上限的风险水平,由调度员设定,为的累积分布函数的逆函数,满足:
αl为K维向量,其第k个元素为
(3)采用内点法,求解上述步骤(2)中的目标函数和约束条件组成的随机动态实时调度模型,得到和将其中的作为t调度时段第i个火电机组的计划出力,第j个发电量自动控制机组的计划出力,作为第k个风电/光伏电站的参考出力,实现基于风/光出力柯西分布的火电机组动态实时调度。
本发明提出的基于风/光出力柯西分布的火电机组动态实时调度方法,其优点是:
本发明方法首先通过多随机变量的联合柯西分布精确刻画了风电/光伏短期预测的出力特性和相关性,以该分布为基础,本方法建立了考虑确定性约束和机会约束下的成本期望值最小化的随机动态实时调度模型,机会约束把调度过程中由于风电/光伏电站和AGC机组出力的随机性带来的安全风险限制在一定的置信水平内。同时,由于柯西分布优良的数学性质使得实时有功调度模型被解析地表达为线性约束的凸优化模型,模型优化的结果是在控制调度风险和减少调度成本下的传统火电机组、可参与调频的AGC机组以及风电/光伏电站出力的最优实时调度决策。本发明充分利用了柯西分布在风电/光伏电站出力的短期预测方面的优越性和优良的数学特性,有效提高了模型的求解效率,同时风险水平可调的机会约束模型消除了传统鲁棒经济调度的保守性,为决策者提供了更合理的调度依据。本发明方法可应用于包含大规模风电并网的电力系统有功实时经济调度中。
具体实施方式
本发明提出的基于风/光出力柯西分布的火电机组动态实时调度方法,包括以下步骤:
(1)建立电网中风电/光伏电站短期预测出力的多随机变量联合柯西分布模型,该柯西分布模型包括以下内容:
(1)建立电网中风电/光伏电站短期预测出力的多随机变量联合柯西分布模型,该柯西分布模型包括以下内容:
a、一个多风电/光伏电站短期预测出力的联合柯西分布的概率密度函数:
其中,Γ为一类基本的数学函数,具体表达式为式中s为任意大于0的自变量,PDF(·)表示随机变量的概率密度函数,概率密度函数通过对风电/光伏电站历史出力的拟合得到,t为调度时段,K为风电/光伏电站的数目,为调度时段t所有风电/光伏电站实际出力的列向量,上标T为矩阵转置,上标w表示该变量描述风电/光伏电站,表示t调度时段第k个风电/光伏电站的实际出力,μt=(μ1,t,μ2,t,...,μk,t)T表示t调度时段联合柯西分布概率密度函数的位置参数,μk,t表示t调度时段第k个风电/光伏电站出力的边缘柯西分布的位置参数,Σt表示t调度时段联合柯西分布概率密度函数的尺度参数;
b、形式为各风电/光伏电站出力线性组合的随机变量的概率描述,包括概率密度函数、累积分布函数以及累积分布函数的逆函数:
设a为k维列向量,那么随机变量服从位置参数和尺度参数为(aTμt,aTΣta)的一维柯西分布,其概率密度函数、累积分布函数以及累积分布函数的逆函数可以表达为以下形式:
其中,概率密度函数为:
累积分布函数为:
累积分布函数的逆函数为:
其中,tan为正切函数,arctan为反正切函数,F为分位数;
(2)建立一个基于风电/光伏电站出力柯西分布的随机动态实时调度模型,该随机动态实时调度模型由目标函数和约束条件组成,具体步骤如下:
(2-1)建立基于风电/光伏电站出力柯西分布的随机动态实时调度模型的目标函数f:
目标函数为使运行成本最小化,表达式如下:
其中,T、N和J分别表示调度时段t的数量、火电机组数量以及发电量自动控制机组(以下简称AGC机组)的数量,t、i和j分别为调度时段、火电机组编号和发电量自动控制机组的编号,上标s表示该变量描述火电机组,上标“+”表示该变量描述正旋转备用,上标“-”表示该变量描述负旋转备用,表示第i台火电机组在t调度时段的计划出力,表示第j个AGC机组在t调度时段内计划出力,表示t调度时段所有风电/光伏电站实际出力的总和,和分别表示火电机组和AGC机组的燃料成本:
其中,ai,t,bi,t,ci,t分别为t时段火电机组i的燃料成本的二次项系数、一次项系数和常数项,aj,t,bj,t,cj,t分别为t时段AGC机组j的燃料成本的二次项系数、一次项系数和常数项,
E(·)表示随机变量的期望值,表示在t调度时段内,由于风电/光伏电站发电实际出力低于计划出力引起的正旋转备用的需求成本,即高估风电/光伏电站出力的惩罚,当风电/光伏电站的实际出力小于计划值时会调度AGC机组的正旋转备用以维持功率平衡,具体表达式如下:
其中,为正旋转备用的成本系数,根据AGC机组的成本特性确定,本发明的一个实施例中,取值为1000,单位为“元/MW”,αj为第j个AGC机组的功率分配系数,根据AGC机组的额定容量占发电量自动控制机组的总容量的比例确定,wt表示t调度时段所有风电/光伏电站的计划出力总和,表示t时段所有风电/光伏电站的实际出力总和,满足如下关系:
其中,为第j个AGC机组的实际出力,表示t调度时段第k个风电/光伏电站的计划出力,K表示风电/光伏电站的数量;
表示在t时段内由于风电/光伏电站实际出力超过计划出力而引起的负旋转备用的需求成本,即低估风电/光伏电站出力的惩罚成本,当风电/光伏电站的实际出力大于计划值时会调度AGC机组的负旋转备用以维持功率平衡,具体表达式如下:
其中,为负旋转备用的成本系数,根据AGC机组的成本特性确定,本发明的一个实施例中,取值为1000,单位为“元/MW”,为实际风电出力总和的概率密度函数,为所有风电/光伏电站在t调度时段出力的上界之和, 为t时段第k个风电/光伏电站的出力的上界;
根据上述步骤(1)的风电/光伏电站短期预测出力的多随机变量联合柯西分布模型,结合方程(2),得到目标函数f的表达式中的后两项为:
其中, 为所有元素全为1的K维列向量,K为风电/光伏电站的数目,A,B,C均为常数具体表达式如下:
基于风电/光伏电站出力柯西分布的随机动态实时调度模型由方程(5),(6),(11)和(12)确定。
2-2)上述基于风电/光伏电站出力柯西分布的随机动态实时调度模型的约束条件,包括:
(2-2-1)电网功率平衡约束,表达式如下:
对
其中,为t调度时段、风电/光伏电站和火电机组所在电网的第d个节点的负荷量,电网中每个节点都有一个负荷,D表示负荷的总数,同时表示电网中节点的个数;
由于由风电/光伏电站实际出力偏离计划出力引起的功率不平衡最终被AGC机组所平衡,因此表达式(13)可写为方程(14):
(2-2-2)电网机组出力的上、下限约束,包括:
对i=1,...,N;j=1,...,J;k=1,...,K:
其中,和分别为第j个AGC机组在t调度时段出力的上下界,和分别为第i个火电机组在t时段出力的上下界。
同时,AGC机组t时段的实际出力以一定的风险水平不得超出其出力的上下界,具体表达式如下:
其中δ为调度员设定的可接受的风险水平。
结合方程(4),机会约束(16)可以转化为确定性的线性约束:
其中,为随机变量的累积分布函数的逆函数,表示t时段所有风电/光伏电站的实际出力总和,函数中各参数为 为所有元素全为1的k维列向量。
(2-2-3)电网中火电机组和发电量自动控制机组的爬坡约束,表达式如下:
对i=1,2,...,N;j=1,2,...,J:
其中,和分别为t调度时段第i台火电机组向上、向下爬坡率;和分别表示t调度时段第j台AGC机组向上、向下爬坡率;ΔT表示相邻两个调度时段之间的调度间隔;β为调度员设定的可接受风险水平;方程(18)和方程(20)分别表示火电机组和AGC机组在相邻时段的爬坡限制,方程(19)表示AGC机组在t调度时段发挥调频作用时的爬坡限制。
结合方程(4),方程(19)和(20)可以转化为如下形式:
其中, 表示随机变量wt,t-1的累积分布函数的逆函数,表达式如下:
(2-2-4)发电量自动控制机组的旋转备用约束,为了平衡由于各种不确定性因素引起的功率波动,系统需要留有一定数量的正、负旋转备用容量,该备用容量受AGC机组的爬坡率和出力上、下界限制,具体表达式如下:
对i=1,2,...,N,j=1,2,...,J,
其中,和分别表示t调度时段第j台AGC机组提供的正、负旋转备用的数量,和分别表示t调度时段电网所需要的最小正、负旋转备用的数量,ε为调度员设定的可接受的风险水平。
结合方程(4),约束(25)可以转化为线性约束(26)
(2-2-5)电网线路潮流约束,表达式如下:
对l=1,2,...,L:
其中,Gi,l为电网中第l条线路对第i台火力发电机组有功出力的转移分布因子,Gj,l为第l条线路对第j个发电量自动控制机组有功出力的转移分布因子,Gk,l为第l条线路对第k个风电/光伏电站有功出力的转移分布因子,Gd,l为第l条线路对第d个节点负荷功率的转移分布因子,各转移分布因子分别从电网调度中心获取,Ll,t为t调度时段第l条线路上的有功功率上限,η为电网线路上的有功功率超过线路有功功率上限的风险水平,由调度员设定。结合方程(4),约束(27)可以转化为线性约束(28),
其中,为随机变量的累积分布函数的逆函数,满足:
αl为k维向量,其第k个元素为
(3)采用内点法,求解由式(5),(6),(8),(9),(11),(12),(14)~(29)确定的随机动态实时调度模型,得到和将其中的作为t调度时段第i个火电机组的计划出力,作为第j个发电量自动控制机组的计划出力,作为第k个风电/光伏电站的参考出力,实现基于风/光出力柯西分布的火电机组动态实时调度。
Claims (1)
1.一种基于风/光出力柯西分布的火电机组动态实时调度方法,其特征在于该方法包括以下步骤:
(1)建立电网中风电/光伏电站短期预测出力的多随机变量联合柯西分布模型,该柯西分布模型包括以下内容:
a、一个多风电/光伏电站短期预测出力的联合柯西分布的概率密度函数:
其中,PDF(·)表示随机变量的概率密度函数,概率密度函数通过对风电/光伏电站历史出力的拟合得到,t为调度时段,K为风电/光伏电站的数目,为调度时段t所有风电/光伏电站实际出力的列向量,上标T为矩阵转置,上标w表示该变量描述风电/光伏电站,表示t调度时段第k个风电/光伏电站的实际出力,μt=(μ1,t,μ2,t,...,μk,t)T表示t调度时段联合柯西分布概率密度函数的位置参数,μk,t表示t调度时段第k个风电/光伏电站出力的边缘柯西分布的位置参数,Σt表示t调度时段联合柯西分布概率密度函数的尺度参数;
b、形式为各风电/光伏电站出力线性组合的随机变量的概率描述,包括概率密度函数、累积分布函数以及累积分布函数的逆函数:
设a为k维列向量,那么随机变量服从位置参数和尺度参数为(aTμt,aTΣta)的一维柯西分布,其概率密度函数、累积分布函数以及累积分布函数的逆函数可以表达为以下形式:
其中,概率密度函数为:
累积分布函数为:
累积分布函数的逆函数为:
其中,tan为正切函数,arctan为反正切函数,F为分位数;
(2)建立一个基于风电/光伏电站出力柯西分布的随机动态实时调度模型,该随机动态实时调度模型由目标函数和约束条件组成,具体步骤如下:
(2-1)建立基于风电/光伏电站出力柯西分布的随机动态实时调度模型的目标函数f:
目标函数为使运行成本最小化,表达式如下:
其中,T、N和J分别表示调度时段t的数量、火电机组数量以及发电量自动控制机组的数量,t、i和j分别为调度时段、火电机组编号和发电量自动控制机组的编号,上标s表示该变量描述火电机组,上标“+”表示该变量描述正旋转备用,上标“-”表示该变量描述负旋转备用,表示第i台火电机组在t调度时段的计划出力,表示第j个发电量自动控制机组在t调度时段内计划出力,表示t调度时段所有风电/光伏电站实际出力的总和,和分别表示火电机组和发电量自动控制机组的燃料成本:
其中,ai,t,bi,t,ci,t分别为t时段火电机组i的燃料成本的二次项系数、一次项系数和常数项,aj,t,bj,t,cj,t分别为t时段发电量自动控制机组j的燃料成本的二次项系数、一次项系数和常数项,
E(·)表示随机变量的期望值,表示在t调度时段内,由于风电/光伏电站发电实际出力低于计划出力引起的正旋转备用的需求成本,即高估风电/光伏电站出力的惩罚,当风电/光伏电站的实际出力小于计划值时会调度发电量自动控制机组的正旋转备用以维持功率平衡,具体表达式如下:
其中,为正旋转备用的成本系数,αj为第j个发电量自动控制机组的功率分配系数,根据发电量自动控制机组的额定容量占发电量自动控制机组的总容量的比例确定,wt表示t调度时段所有风电/光伏电站的计划出力总和,表示t时段所有风电/光伏电站的实际出力总和,满足如下关系:
其中,为第j个发电量自动控制机组的实际出力,表示t调度时段第k个风电/光伏电站的计划出力,K表示风电/光伏电站的数量;
表示在t时段内由于风电/光伏电站实际出力超过计划出力而引起的负旋转备用的需求成本,即低估风电/光伏电站出力的惩罚成本,当风电/光伏电站的实际出力大于计划值时会调度发电量自动控制机组的负旋转备用以维持功率平衡,具体表达式如下:
其中,为负旋转备用的成本系数,为实际风电出力总和的概率密度函数,为所有风电/光伏电站在t调度时段出力的上界之和, 为t时段第k个风电/光伏电站的出力的上界;
根据上述步骤(1)的风电/光伏电站短期预测出力的多随机变量联合柯西分布模型,得到目标函数f的表达式中的后两项为:
其中, 为所有元素全为1的K维列向量,K为风电/光伏电站的数目,A,B,C均为常数具体表达式如下:
(2-2)上述基于风电/光伏电站出力柯西分布的随机动态实时调度模型的约束条件,包括:
(2-2-1)电网功率平衡约束,表达式如下:
其中,为t调度时段、风电/光伏电站和火电机组所在电网的第d个节点的负荷量,D表示负荷的总数,同时表示电网中节点的个数;
(2-2-2)电网机组出力的上、下限约束,包括:
其中,t=1,...,T,i=1,...,N,j=1,...,J,k=1,...,K,和分别为第j个发电量自动控制机组在t调度时段出力的上下界,和分别为第i个火电机组在t时段出力的上下界,δ为调度员设定的可接受的风险水平,为的累积分布函数的逆函数,表示t时段所有风电/光伏电站的实际出力总和,
为所有元素全为1的K维列向量,K表示风电/光伏电站的数量;
(2-2-3)电网中火电机组和发电量自动控制机组的爬坡约束,表达式如下:
对
其中,和分别为t调度时段第i台火电机组向上、向下爬坡率,和分别表示t时段第j台发电量自动控制机组向上、向下爬坡率,ΔT表示相邻两个调度时段之间的调度间隔,β为调度员设定的可接受风险水平,表示wt,t-1的累积分布函数的逆函数, 的表达式如下:
(2-2-4)发电量自动控制机组的旋转备用约束,具体表达式如下:
对
其中,和分别表示t调度时段第j台发电量自动控制机组提供的正负旋转备用的数量,和分别表示t调度时段电网所需要的最小正、负旋转备用的数量,ε为调度员设定的可接受的风险水平;
(2-2-5)电网线路潮流约束,表达式如下:
对
其中,Gi,l为电网中第l条线路对第i台火力发电机组有功出力的转移分布因子,Gj,l为第l条线路对第j个发电量自动控制机组有功出力的转移分布因子,Gk,l为第l条线路对第k个风电/光伏电站有功出力的转移分布因子,Gd,l为第l条线路对第d个节点负荷功率的转移分布因子,各转移分布因子分别从电网调度中心获取,Ll,t为t调度时段第l条线路上的有功功率上限,η为电网线路上的有功功率超过线路有功功率上限的风险水平,由调度员设定,为的累积分布函数的逆函数,满足:
αl为K维向量,其第k个元素为
(3)采用内点法,求解上述步骤(2)中的目标函数和约束条件组成的随机动态实时调度模型,得到和将其中的作为t调度时段第i个火电机组的计划出力,第j个发电量自动控制机组的计划出力,作为第k个风电/光伏电站的参考出力,实现基于风/光出力柯西分布的火电机组动态实时调度。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811509293.0A CN109755959B (zh) | 2018-12-11 | 2018-12-11 | 基于风/光出力柯西分布的火电机组动态实时调度方法 |
PCT/CN2019/100600 WO2020119159A1 (zh) | 2018-12-11 | 2019-08-14 | 基于风/光出力柯西分布的火电机组动态实时调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811509293.0A CN109755959B (zh) | 2018-12-11 | 2018-12-11 | 基于风/光出力柯西分布的火电机组动态实时调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109755959A true CN109755959A (zh) | 2019-05-14 |
CN109755959B CN109755959B (zh) | 2020-07-10 |
Family
ID=66403544
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811509293.0A Active CN109755959B (zh) | 2018-12-11 | 2018-12-11 | 基于风/光出力柯西分布的火电机组动态实时调度方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN109755959B (zh) |
WO (1) | WO2020119159A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020119159A1 (zh) * | 2018-12-11 | 2020-06-18 | 清华大学 | 基于风/光出力柯西分布的火电机组动态实时调度方法 |
CN111416396A (zh) * | 2020-03-31 | 2020-07-14 | 清华大学 | 辅助服务市场中考虑电采暖的火电与风电联合调度方法 |
CN113721476A (zh) * | 2021-06-29 | 2021-11-30 | 武汉大学 | 100mw级可变速海水抽水蓄能机组与可再生能源联合运行系统硬件在环仿真平台及方法 |
CN114069621A (zh) * | 2021-11-16 | 2022-02-18 | 南京邮电大学 | 计及多能源系统稳定性的多目标协同优化安全调度方法 |
CN114188942A (zh) * | 2021-12-09 | 2022-03-15 | 国网甘肃省电力公司电力科学研究院 | 一种包含大规模新能源基地的电网调度方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015031331A1 (en) * | 2013-08-26 | 2015-03-05 | Robert Bosch Gmbh | Dispatch controller for an energy system |
CN105591407A (zh) * | 2016-03-03 | 2016-05-18 | 国家电网公司 | 可再生能源电站有功出力预测误差间相关性的研究方法 |
CN106327091A (zh) * | 2016-08-26 | 2017-01-11 | 清华大学 | 基于鲁棒联络线计划的多区域异步协调动态经济调度方法 |
CN106485362A (zh) * | 2016-10-18 | 2017-03-08 | 江苏省电力试验研究院有限公司 | 一种基于高维风电预测误差模型及降维技术的发电调度方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103050998B (zh) * | 2012-11-26 | 2014-11-05 | 西安理工大学 | 一种风电并网的火电系统动态调度方法 |
CN104242354B (zh) * | 2014-06-18 | 2018-01-05 | 国家电网公司 | 计及风光出力相关性的新能源集中外送运行特性评估方法 |
US20170371306A1 (en) * | 2016-06-27 | 2017-12-28 | Ecole Polytechnique Federale De Lausanne (Epfl) | System and Method for Dispatching an Operation of a Distribution Feeder with Heterogeneous Prosumers |
CN107330546A (zh) * | 2017-06-14 | 2017-11-07 | 武汉大学 | 一种考虑风电出力和需求响应不确定性的优化调度方法 |
CN108879787B (zh) * | 2018-08-17 | 2021-04-13 | 合肥工业大学 | 一种含风电的电网随机调度优化模型与方法 |
CN109755959B (zh) * | 2018-12-11 | 2020-07-10 | 清华大学 | 基于风/光出力柯西分布的火电机组动态实时调度方法 |
-
2018
- 2018-12-11 CN CN201811509293.0A patent/CN109755959B/zh active Active
-
2019
- 2019-08-14 WO PCT/CN2019/100600 patent/WO2020119159A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015031331A1 (en) * | 2013-08-26 | 2015-03-05 | Robert Bosch Gmbh | Dispatch controller for an energy system |
CN105591407A (zh) * | 2016-03-03 | 2016-05-18 | 国家电网公司 | 可再生能源电站有功出力预测误差间相关性的研究方法 |
CN106327091A (zh) * | 2016-08-26 | 2017-01-11 | 清华大学 | 基于鲁棒联络线计划的多区域异步协调动态经济调度方法 |
CN106485362A (zh) * | 2016-10-18 | 2017-03-08 | 江苏省电力试验研究院有限公司 | 一种基于高维风电预测误差模型及降维技术的发电调度方法 |
Non-Patent Citations (5)
Title |
---|
GUOQIANG ZHANG, BOMING ZHANG: "Ultra-Short Term Probabilistic Transmission Congestion Forecasting Considering Wind Power Integration", 《8TH INTERNATIONAL CONFERENCE ON ADVANCES IN POWER SYSTEM CONTROL, OPERATION AND MANAGEMENT (APSCOM 2009)》 * |
LEI XU, LAN ZHANG, JIEDANG: "Quantum-behaved particle swarm optimization with Cauchy disturbance for power economic dispatch", 《2012 4TH INTERNATIONAL CONFERENCE ON INTELLIGENT HUMAN-MACHINE SYSTEMS AND CYBERNETICS》 * |
李志刚: "消纳大规模风电的鲁棒区间经济调度", 《电力系统自动化》 * |
汲国强: "一种适用于可靠性评估的电网设备时变停运模型", 《中国电机工程学报》 * |
黄麒元: "基于智能算法的大型风电场功率爬坡预测方法研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020119159A1 (zh) * | 2018-12-11 | 2020-06-18 | 清华大学 | 基于风/光出力柯西分布的火电机组动态实时调度方法 |
CN111416396A (zh) * | 2020-03-31 | 2020-07-14 | 清华大学 | 辅助服务市场中考虑电采暖的火电与风电联合调度方法 |
CN113721476A (zh) * | 2021-06-29 | 2021-11-30 | 武汉大学 | 100mw级可变速海水抽水蓄能机组与可再生能源联合运行系统硬件在环仿真平台及方法 |
CN114069621A (zh) * | 2021-11-16 | 2022-02-18 | 南京邮电大学 | 计及多能源系统稳定性的多目标协同优化安全调度方法 |
CN114069621B (zh) * | 2021-11-16 | 2023-08-22 | 南京邮电大学 | 计及多能源系统稳定性的多目标协同优化安全调度方法 |
CN114188942A (zh) * | 2021-12-09 | 2022-03-15 | 国网甘肃省电力公司电力科学研究院 | 一种包含大规模新能源基地的电网调度方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2020119159A1 (zh) | 2020-06-18 |
CN109755959B (zh) | 2020-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109755959B (zh) | 基于风/光出力柯西分布的火电机组动态实时调度方法 | |
CN107276127B (zh) | 考虑联络线电力交易计划的多区域互联电力系统的风电消纳优化方法 | |
CN104242356B (zh) | 考虑风电场集电线故障的鲁棒区间风电调度方法及装置 | |
CN108092324B (zh) | 一种风电参与调峰调频的agc控制系统和控制方法 | |
CN109840636B (zh) | 一种基于牛顿法的电力系统随机滚动调度方法 | |
CN108695857B (zh) | 风电场自动电压控制方法、装置及系统 | |
CN113644670B (zh) | 一种储能容量优化配置的方法及系统 | |
CN106549420B (zh) | 考虑风险和风力发电的电力系统运行备用优化方法 | |
US11689024B2 (en) | Bilateral stochastic power grid dispatching method | |
CN109713713B (zh) | 一种基于机会约束凸松弛的机组启停随机优化方法 | |
CN104037805A (zh) | 一种计及电网安全约束的光伏电站可发电裕度分配方法 | |
CN111541244A (zh) | 考虑储能装置度电成本的电网侧储能装置容量计算方法 | |
CN110826773A (zh) | 一种考虑新能源接入的火电机组月度发电计划优化方法 | |
CN111641233A (zh) | 一种考虑新能源及负荷不确定性的电力系统日内灵活调峰方法 | |
CN111626470A (zh) | 一种电热综合协调优化调度方法及系统 | |
CN110867907B (zh) | 一种基于多类型发电资源同质化的电力系统调度方法 | |
CN109299862A (zh) | 一种风电最大消纳能力的凸松弛评估方法 | |
CN109657898B (zh) | 一种基于凸松弛的可再生能源随机动态经济调度方法 | |
CN115360706A (zh) | 考虑dr和灵活性供需平衡的源荷储联合调度方法及系统 | |
CN115114854A (zh) | 一种虚拟电厂分布式资源两级自组织优化聚合方法及系统 | |
CN112994087B (zh) | 基于条件风险约束的多源电力系统中期优化调度方法 | |
CN115473284B (zh) | 区域电力交换约束下配电系统运行的鲁棒优化方法、系统及计算机设备 | |
CN116454944A (zh) | 一种基于随机生产模拟的储能装置优化配置方法及系统 | |
CN115800276A (zh) | 一种考虑机组爬坡的电力系统应急调度方法 | |
CN116341819A (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 |