一种基于通用分布的含风电电力系统实时动态经济调度方法
【技术领域】
[0001] 本发明属于电力系统运行与控制领域，设及一种基于通用分布的含风电电力系统 实时动态经济调度方法。
【背景技术】
[0002] 随着风力发电的大规模开发，风电对电网安全运行和调度控制等带来巨大挑战， 一种将风电作为调度量的实时动态经济调度方法对电力安全生产至关重要。
[0003] 经济调度问题是电力系统的基本问题，在风电大规模接入之前，由于火电机组良 好的可控性及电力系统负荷成熟的预测方法，W往的经济调度通常是日前调度和实时手动 调节的组合。风电的接入，风电的自然属性(随机性和波动性)成为电力系统的不确定电源， 在风电渗透率不高的情况下和负荷一样作为电力系统扰动源。随着风电渗透率的不断增 加，风电的随机针对电力系统的影响越来越大。由于日前风功率预测误差较大，而风电的超 短期预测（化)精度较高，故国内外的研究大都是通过日前机组组合、日前经济调度和实时 经济调度配合，多尺度进行调度。在含风电调度方面，国内外学者进行了大量基础研究，可 分为四类。
[0004] (1)基于动态规划法的动态经济调度。此类算法解析度通常较高，算法逻辑清晰， 计算速度快，例如等微增率算法。不足之处为此类算法通常将风电作为净负荷处理，如前所 述，运种处理方式在风电渗透率不高的风电发展初期较为适用。在风电渗透率达到一定规 模后，此类算法的弊端凸显。
[0005] (2)基于启发式方法的动态经济调度。此类算法的数学模型相对复杂，传统的算法 难W求解，故通常适用蚁群、粒子群等启发式算法求解。启发式算法虽然智能性高，可W保 证求出解，但每次优化结果是不确定的，有事无法保证得到的解是全局最优解，部分启发式 算法不适合用于大规模系统求解。
[0006] (3)基于鲁棒优化方法的动态经济调度。鲁棒优化的核屯、思想是考虑最极端的边 界情况，达到优化方法的鲁棒性。即牺牲一定经济性的条件下，保证模型对极端情况的有效 处理。该方法可W有效求解发电计划，但由于鲁棒优化固有的鲁棒特性，在经济性等方面会 有一定的降低。
[0007] (4)基于随机优化方法的动态经济调度。随机优化算法为将风电功率作为随机变 量处理，场景法也是随机优化方法的一种常用方法。随机优化是一种处理含不确定性优化 问题的有效方法，目前已广泛应用于含不确定性的电力系统经济调度问题中。如何准确描 述风电功率的不确定性、有效求解对应的优化模型是含风电电力系统随机经济调度的关键 问题。此类算法物理概念较为明确，但在计算速度方面目前还有很大的提升空间，尤其是对 于对计算速度有较高要求的实时调度层面。
[000引动态规划法由于其计算速度快，在实时调度中具有很大的潜力，但动态规划法大 都将风电视为净负荷处理。在随机优化中，如何合理描述风电功率的不确定性并将其运用 于电力系统的经济调度与优化运行中具有重要意义。描述风电功率的分布常用正态分布。 对风电功率分布描述的准确性与对应随机动态经济调度模型的解的准确性密切相关。但目 前尚未出现相关具有实用价值的技术方案。
【发明内容】
[0009] 本发明针对现有技术的缺陷，提供基于通用分布的电力系统实时动态经济调度技 术方案。
[0010] 本发明提供的技术方案是：
[0011] 1、一种基于通用分布的含风电电力系统实时动态经济调度方法，其特征在于，包 括W下步骤：
[0012] 步骤1:输入日前系统负荷预测数据和日前风电功率预测数据、火电机组运行参 数、历史统计数据，所述历史统计数据包括不同风电功率预测水平下实际风电功率的通用 分布参数α,β, 丫；
[OOK]步骤2:定义pi,t为第i台火电机组t时刻的出力，火电机组的总数为I，i = l，2,...， I，wj, t为第j个风电场t时刻的计划出力，风电场的总数为J，j = 1，2，. . .，J，T为时刻的总数， t = l,2,. . . ,Τ；
[0014] 计算从当前时刻也就是当前计算周期到后面1小时的调度时间域每个调度周期火 电机组出力可能达到的范围，称为火电机组在t调度周期的可行域，该1小时的调度时间域 为12个调度周期：
[0015]
式一
[0016] 式中：Pi康示t时刻第i台火电机组的出力;單t和焉是表示t时刻第i台火电机组的 出力下限和上限；
[0017] 使用通用分布表征风功率，使用等微增率准则进行12个调度周期的风电火电静态 经济调度，得到风电和火电的经济调度曲线，基于W下的目标函数W及约束条件：
[001引 目标函数；
[0019] 火电机组出力功率约束条件：A 化户,·式Ξ
[0020] 火电机组爬坡约束条件:-Si如i广Pi(t+i)<Si式四
[0021] 风电场出力功率约束条件：0 <式五
[0022] 系统功率平衡约束条件
式六
[0023] 式中：ft fg,t fw,t分别表示t时刻系统总发电成本、火电发电成本和风电成本;wjt 表示*时刻第^'个风电场的出力；1表示火电机组数量，1表示风电场数量;？*^表示*时刻系统 总的负荷预测值；兵和￡康示第i台火电机组的出力上下限;S康示第i台火电机组的爬坡能 力；跨表示第j个风电场的风电场装机容量；
[0024] 其中，所述火电燃料成本Cg,i,t采用W下方式计算，
[0025]
式屯
[0026] 其中，Pit为第i台火电机组t时刻的出力，ai，bi，ci为第i台火电机组的燃料成本系 数；
[0027] 所述风电成本Cwj,t采用W下方式计算，
[002引
式八
[0029] 其中，为第j个风电场t时刻的计划出力，山为第j个风电场的运行成本系数；
[0030] 所述平均低估成本Cuw,t采用风电场的平均弃风成本，采用W下方式计算，
[0031]
式九
[0032] 其中，kun,j为第j个风电场对应的低估成本系数，Wav,j,t为第j个风电场t时刻的实 际可能出力，为第j个风电场在对应风电预测水平下实际可能出力的概率密度函 数，表达形式为通用分布的概率密度函数，根据t时刻对应通用分布参数α，β，丫确定；
[0033] 所述平均高估成本Cdv,It采用系统的平均备用成本，采用W下方式计算，
[0034]
式十
[0035] 其中，kDv,j为第j个风电场对应的高估成本系数，Wav,j,t为第j个风电场t时刻的实 际可能出力，为第j个风电场在对应风电预测水平下实际可能出力的概率密度函 数，表达形式为通用分布的概率密度函数，根据t时刻对应通用分布参数α，β，丫确定；
[0036] 步骤3:利用步骤2得到的12个调度周期各自的风电场出力，修正火电机组各个调 度周期的可行域;在新的火电机组可行域内，使用等微增率准则对火电机组的出力进行修 正，保证火电机组爬坡限制并保证火电机组发电计划的稳定性；
[0037] 步骤4:利用步骤2得到的风电场出力和步骤4得到的火电机组出力，输出火电机组 和风电场的计划出力；
[0038] 而且，所述火电燃料成本Cg,i,t采用W下方式计算，
[0039]
[0040] 其中，Pit为第i台火电机组t时刻的出力，ai，bi，ci为第i台火电机组的燃料成本系 数；
[0041] 而且，所述风电成本采用W下方式计算，
[0042] Cw,j,t(wj,t)=djWjt+Cun,j,t(wj,t)+C〇v,j,t(wj,t)式十二
[0043] 其中，为第j个风电场t时刻的计划出力，山为第j个风电场的运行成本系数；
[0044] 而且，所述平均低估成本Cun, j，t采用风电场的平均弃风成本，采用W下方式计 算，
[0045]
[0046] 其中，kun,j为第j个风电场对应的低估成本系数，Wav,j,t为第j个风电场t时刻的实 际可能出力，为第j个风电场在对应风电预测水平下实际可能出力的概率密度函 数，表达形式为通用分布的概率密度函数，根据t时刻对应通用分布参数α，β，丫确定；
[0047] 而且，所述平均高估成本Cdv,It采用系统的平均备用成本，采用W下方式计算，
[004引
[0049] 其中，kDv,j为第j个风电场对应的高估成本系数，Wav,j,t为第j个风电场t时刻的实 际可能出力，为第j个风电场在对应风电预测水平下实际可能出力的概率密度函 数，表达形式为通用分布的概率密度函数，根据t时刻对应通用分布参数α，β，丫确定。
[0050] 本发明利用通用分布刻画了风电功率的不确定性，W此为基础，提出了考虑风电 低、高估成本的实时动态经济调度技术方案，包括W风电场的历史风电功率数据为基础，利 用通用分布模型拟合不同风电功率预测水平下实际风电功率的分布，计及风电预测误差带 来的惩罚成本，建立了基于通用分布的考虑风电低、高估成本的实时动态经济调度的随机 优化模型；通过计算火电机组在各个调度周期的可行域，将对应的动态经济调度模型转化 为12个静态经济调度模型，得到12个调度周期的火电机组和风电场出力，并将随机优化模 型转化为确定性优化模型；使用动态规划算法修正火电机组的出力，求解动态经济调度模 型，得到火电机组和风电
场的实时计划出力。经验证，可知本发明技术方案的有效性，具有 良好的推广价值和应用前景。
【附图说明】
[0051] 图1是本发明实施例的通用分布与正态分布对于拟合风功率分布的拟合效果比 较。
[0052] 图2是本发明实施例的基于通用分布模型应用的含风电电力系统调度结构示意 图。
[0053] 图3是本发明实施例的风火静态经济调度等微增率算法。
[0054] 图4是本发明实施例的火电机组在调度时间域的各个调度周期的运行可行域。
[0055] 图5是本发明实施例的火电机组基本发电计划。
[0056] 图6是本发明实施例的火电机组发电计划的修正过程。
[0057] 图7是本发明实施例的火电机组发电计划的修正等微增率算法。
[005引图8是本发明实施例的ΙΕ邸10机39节点系统图。
[0059] 图9是本发明实施例的日负荷曲线。
[0060] 图10是本发明实施例的丹麦海上风电预测与实测值。
[0061] 图11是本发明实施例的丹麦海上风电通用参数表。
[0062] 图12是本发明实施例的风电调度值和实测值。
[0063] 图13是本发明实施例的两种调度模式的系统成本比较。
[0064] 图14是本发明实施例的不同发电机数目下的算法计算时间。
[0065] 图15是本发明实施例的采用不同分布表征风功率的计算时间。
[0066] 图16是本发明实施例的调度模型对于保证火电机组发电计划稳定性的结果比较 分析。
[0067] 图17是本发明实施例的采用不同方式表征风功率的随机性对系统频率的影响比 较。
【具体实施方式】
[0068] 为了使本发明实施例的目的、技术方案、优点更加清晰，下面将结合本发明实施例 和附图来介绍本发明的技术方案。
[0069] 本发明提供的技术方案是一种基于通用分布的考虑风电低、高估成本的实时动态 经济调度方法，原理如下：
[0070] 将并网风电场历史风电功率的预测和实测数据标么化，根据风电功率的预测值的 不同，对历史风电功率数据进行分箱，在不同风电功率预测水平下，利用通用分布函数拟合 不同预测箱下实测风电功率的分布，得到对应的通用分布参数；
[0071] 考虑系统有功功率平衡约束、火电机组和风电场的出力限制、火电机组爬坡速率 约束，建立基于通用分布的计及风电低、高估成本的含风电电力系统的实时动态经济调度 模型；
[0072] 基于预测风电功率的随机经济调度模型，使用通用分布表征风功率，并计算各个 调度周期火电机组的出力可行域，使用等微增率准则进行12个调度周期的风电火电静态经 济调度，得到考虑风电低、高估成本的实时随机动态经济调度模型的风电和火电的经济调 度曲线。
[0073] 得到各个调度周期风电场出力计划后，基于通用分布的随机动态经济调度问题转 化成确定性经济调度模型。利用步骤二得到的12个调度周期各自的风电场出力，修正火电 机组各个调度周期的可行域。在新的火电机组可行域内，使用等微增率准则对火电机组的 出力进行修正，保证火电机组爬坡限制并保证火电机组发电计划的稳定性。最终输出日前 火电机组和风电场的计划出力曲线。
[0074] 首先，为便于实施参考起见，介绍基于通用分布的随机优化模型：
[0075] -、通用分布模型和形状参数表的生成。
[0076] 1.1、通用分布模型。
[0077] 若连续型随机变量X服从一个形状参数为α、β和丫的通用分布，则记为：
[007引 X~ν(α，β，丫） (11)
[0079] 其中，形状参数α、β和丫满足；
[0080] a>〇，e>〇，-w< 丫 <+w (12)
[0081 ] 通用分布的概率密度函数(Probability Density F^mctionJOF)定义为：
[00 剧
（13)
[0083] 通用分布的累积分布函数(Cumulative Distribution Function,CDF)定义为：
[0084] F(x) = (i+e_a(x-丫）)-e (14)
[00化]1.2、风电概率密度直方图。
[00化]1.2.1、绘制风电概率密度直方图。
[0087] 统计风电场预测和实测历史数据，得到风电功率实测时间序列和预测时间序列W 后，W风电场最大出力为基值，将数据转化为标么值，供后文研究使用。
[0088] 本文所设及的几个关键名词的定义如下：
[0089] (1)风电预测功率等级:将风电预测功率的取值范围0~Ip.u.等分为Μ个小区间， 各区间的宽度为1/Μ，取值范围分别为[0，1朋，[1/^，2/^]^。，[(1-1)/^，1]，每个区间称作 一个"风电预测功率等级"，第m级(m=l，2，···，Μ)对应的预测功率取值范围是[(m-l)/M，m/ Μ] ο
[0090] (2)第m个预测功率等级对应的风电实测功率原始数据集合:在原始数据集合中的 每一个时段，风电场有一个预测功率值和一个实测功率值，二者组成数对。对于每个属于第 m级预测功率等级的风电预测功率值，都有一个与之成对的实测功率值，所有运些实测功率 值构成的集合称作"第m个预测功率等级对应的风电实测功率原始数据集合"。
[0091] 值得一提的是，Μ为用户自定义的正整数，可用的历史数据越多，Μ取值越大，反之 越小。另一方面，如果历史数据总量固定，Μ取值越大，每个预测功率等级对应的风电实测功 率原始数据集合中样本数量越少，反之越多。对于每个预测功率等级对应的风电实测功率 原始数据集合，如果该集合的样本数量过少，将降低统计结果的可信度，运就失去了统计的 意义。因此，为保障后文中统计分析的效果，应使每个预测功率等级对应的风电实测功率原 始数据集合中样本数量足够多。经试验，对于本文采用的风电场数据，取Μ = 25可W达到较 好的统计效果。
[0092] 使用直方图表征风电功率概率分布，概率密度直方图的绘制步骤如下：
[0093] Step 1:收集数据。对于不同的风电预测功率等级，实际风电功率的概率分布不 同。对于第m个预测功率等级，选取其对应的风电实测功率原始数据集合作为样本集合。
[0094] Step 2:确定直方图的总宽度。传统做法是：由样本集合中的最大值与最小值之差 求得极差，直方图的总宽度（即组距与组数的乘积)取略大于极差的数。由于风电出力在0~ 1P. U.范围内，为简化问题，本文假设直方图的跨度也在0~Ip. U.范围内，于是总宽度就等 于Ip.u.。
[00M] step 3:确定直方图的组数、组距和各组界限。假设直方图的组数为N;每组的宽度 为1/N，即组距;第η组(η = 1，2，···，Ν)的下界为(n-l)/N，上界为n/N。注意，组数的确定要适 当，组数过少可能会使数据的分布规律失真，组数过多可能会影响数据分布规律的明显性。 经试验，对于本文采用的风电场数据，取N=50可W达到较好的统计效果。
[0096] Step 4:确定各组频数(Frequency Counts)、概率和概率密度。把样本集合中的数 据按照大小划分到直方图各组，统计各组中数据的个数，即频数;各组频数除W总样本数， 就得到各组概率，用probn(n = l，2，-'，N)表示；各组概率除W组距1/N，就得到各组的概率 密度，即N · probn。
[0097] Step 5:绘制概率密度直方图。横坐标为风电实际功率，纵坐标为概率密度，直方 图每组对应一个矩形，矩形的宽度为组距，高度为各组的概率密度。
[0098] 由于直方图是由真实数据统计而来，没有采用任何近似方法，故本文将使用直方 图表征的风电功率概率分布称为"实际分布"。实际分布本质上是离散型分布，其PDF定义如 式(5)所示：
[0099]
[0100] 对直方图各组概率离散求和就可W得到实际分布的CDF，如式(6)所示：
[0101]
[0102] 1.2.2、曲线拟合得到通用分布形状参数表。
[0103] 正态和贝塔分布近似模型的参数均是通过所统计的样本集合的平均值μ和标准差 曰直接或间接推导求得。运使得正态和贝塔分布的PDF/CDF曲线形状与μ和σ强相关，而与实 际分布的PDF/CDF曲线形状弱相关，导致正态和贝塔分布近似模型模拟实际分布的效果较 差。
[0104] 为了能够更精确地模拟实际分布，本文提出：当使用通用分布近似表征实际分布 时，不通过μ和σ求取通用分布的形状参数α、β和γ，而是使用通用分布PDF/CDF曲线直接对 实际分布PDF/CDF曲线进行最小二乘拟合得到α、β和丫，请见图1，通用分布拟合效果相比于 其他分布能够更加准确的拟合风功率实际分布。由于实际分布是W风电预测功率为均值的 条件概率分布，故对于不同的风电预测功率等级，通用分布的形状参数不同。对于某一预测 时间尺度下的风电场，可W通过拟合各风电预测功率等级下的实际分布PDF或CDF，生成一 个通用分布的形状参数表化ook-Up Table)。
[0105] 需要强调的是，PDF拟合法和CDF拟合法得到的通用分布形状参数大小可能有一定 差异。那么，到底应该使用哪种拟合方法获取通用分布的形状参数呢？在经济调度问题中，
一些经典算法通常需要计算低/高估风电代价的偏导数并将约束条件线性化，低/高估风电 代价的偏导数和机会约束是关于风电功率概率分布CDF或其逆函数的表达式。所W，CDF拟 合法得到的通用分布形状参数更实用、更直接和更精确。当求解实际电力系统经济调度问 题时，应尽量采用CD巧ft合法求取通用分布的形状参数。
[0106] 另外需要注意的是，最小二乘拟合解(通用分布的Ξ个形状参数α、β和丫）并不唯 一，主要取决于初值如何选取。如果随机选取每个风电预测功率等级条件下的通用分布形 状参数的初值，可能使某个预测功率等级对应的拟合解与其它预测功率等级对应的拟合解 在正负号或数量级方面产生较大差异。当然，虽然有些拟合解的正负号或数量级"特殊"，但 只要相应的通用分布曲线能够较好地拟合实际分布曲线，运些看似"特殊"的拟合解就具有 实用性，所选取的初值就是可W接受的。
[0107] 二、基于通用分布的风电火电静态经济调度。
[0108] 2.1、基于通用分布的风功率随机性表征。
[0109] 基于通用分布模型应用的含风电电力系统调度结构示意图如图2所示，在调度计 划制定之前，风电场提供实测功率历史数据，风电预测系统提供预测功率历史数据，采用上 节所述的方法，通过使用通用分布CDF对实际分布CDF进行最小二乘拟合，得到通用分布形 状参数表。在调度周期内，风电功率预测系统给定当前风电功率预测值，该预测值属于某一 预测功率等级，从通用分布形状参数表中可查询当前预测功率等级条件下的形状参数，进 而得到当前预测功率等级条件下的通用分布模型，然后将其应用到调度模型和算法中。最 终得到调度系统向风电场和常规机组发送的调度计划指令。
[0110] 如果风电场的地理布局、装机容量和风机类型等条件过去没有太大的变化，则为 生成通用分布的形状参数表，图2中的统计分析环节应采用全部的风电预测和实测功率历 史数据。为提高统计分析的精度，应定期（如每天或每星期)将最新的数据补充到历史数据 库中。生成通用分布形状参数表的过程并不会影响调度系统计算和决策的效率，原因如下： 第一，生成通用分布形状参数表的耗时很短，例如，在一台Core-i72.70G化笔记本电脑中的 Matlab 2012a环境下耗时不到1S;第二，形状参数表无需实时更新，每天或每星期更新一次 即可。
[0111] 2.2、风电火电静态经济调度。
[0112] 对于实时经济调度，其目标函数和约束条件如下：
[0113] 目标函数：
[0114]
(2)
[0115] 火电机组出力功率约束：
[0116]
(3)
[0117]火电机组爬坡约束：
[011 引
[0119] 风电场出力功率约束：
[0120]
(5)
[0121] 系统功率平衡约束：
[0122]
(6)
[0123] 式中：ft fg,t fw,t分别表示t时刻系统总发电成本、火电发电成本和风电成本;Pit 表示t时刻第i台火电机组的出力;ww表示t时刻第j个风电场的出力；I表示火电机组数量J 表示风电场数量;PtL表示t时刻系统总的负荷预测值；兵和￡i表示第i台火电机组的出力上 下限;δι表示第i台火电机组的爬坡能力；朽表示第j个风电场的风电场装机容量。
[0124] 目标函数中风电成本展开为：
[0125]
(1巧
[0126] 式中：山是风电场j的运行费用;kuw是风电场j的低估代价费用;kDv,j是风电场j的 高估代价费用;Wav,W是风电场j的在t时刻的实际可能出力。
[0127] 风电成本等微增率：
[012 引
（18)
[0129] 根据(144)式的通用分布函数CDF，（18)式可化为：
[0130]
(19)
[0131] 根据火电成本的二次曲线，火电成本等微增率可按下式计算：
[0132]
(20)
[0133] 式中：ai和bi表示火电机组i的二次成本函数系数。
[0134] 根据火电机组初始运行点，可通过式(3)(4)(6)计算出调度域所有调度周期火电 机组的功率范围，即可行域，如图4。按照可行域的计算方法，在计算可行域上限的时候，风 电功率取最小值，即0;在计算可行域下限的时候，风电功率取最大值，即风电场装机容量。 暂不考虑可行域中的爬坡因素，火电机组可行域上下限可由下式表示：
[0135] a"<p"<c\ (1)
[0136] 在火电机组可行域（1)和风电场装机容量(5)内，可通过经典的等微增率计算调度 域(化）内12个调度周期(5min)各自的风电场和火电机组出力，如图3。
[0137] 3、基于动态规划法的火电机组出力修正方法。
[013引 3.1、火电机组基本发电计划。
[0139] 火电机组基本发电计划的引入是降低不同调度周期对于调度域同一调度周期发 电计划的更改，保证发电计划的稳定性。故某一调度周期的基本发电计划由上一调度周期 计算得到但未使用的11个调度周期的火电出力和此调度周期计算的最远调度周期的火电 计划组成，如图5所示。
[0140] 3.2、火电机组运行可行域和系统功率平衡修正。
[0141] 风电火电静态经济调度得到调度域各调度时间风电发电计划后，火电机组运行可 行域可W使用风电发电计划进行修正。同时，系统功率平衡也可修正为：
[0142]
(21)
[0143] 式中:ΡΛ表示t调度周期下的净负荷(系统总负荷减去风电出力）。
[0144] 3.3、火电机组出力修正。
[0145] 火电机组发电计划可W通过在每调度周期基本发电计划加上增量表示：
[0146]
(22)
[0147] 注意到姑是基于上一计算周期的净负荷P严'求得的，基于本计算周期净负荷P严 得到Pit,故A Pit,…Δ pii即表示两个计算周期对于同一调度周期的净负荷变化，体现在如 下的功率平衡中：
[0150] 为计算Δ Pit,定义爬坡率变化值：Δ化t
[0151] Δ Dit= Δ Pit-Δ pi(t+：L) (24)
[0152] ADii，Δ〇ι2. . . ΔΟιτ可W通过从较远的调度周期向较近的调度周期通过等微增率 准则进行计算。在计算过程中，当计算到t = s阶段，即计算Δ Dis，假设对于t〉s Δ化t = 〇。即 假设：^0如+1片〇
[0153] 运样，在修正过程中，只有Δ化sa = l...I)为变量。在修正到每一个调度周期时， A Dis的值决定所有的Δ pit(t = 1. . . s)值：
[0154]
(25)
[0155] 当ADis变化时，t = s...l调度周期的火电机组发电计划根据式(25)平行移动，如 图6。由于等微增率准则的计算速度非常快，故修正过程耗时很短。在每一个调度周期的修 正过程中，等微增率准则的各项参数如图7。
[0156] 其中；
[0159] 当s = 12,所有调度周期的Δ化t计算完毕，修正后的火电机组出力即：
[0160]
(28)
[0161] 四、算例及仿真。
[0162] 本发明WIEEE10机39节点系统标准算例为例进行了仿真计算。系统10台发电机 中，G1为调频机组，不进行调度，G7、G9只相应日前调度出力曲线，不参与实时动态经济调 度，如图8所示。风电使用丹麦海上风电实际数据，在8、15和39接入，设置风电运行成本24$/ MWh，风电低估成本60$/MWh，风电高估成本120$/MWh。实时调度算法计算一天的机组调度出 力，整个系统使用RTDS进行一天的仿真计算，7台化0机组（即实时调度机组）的出力由实时 调度计算结果决定，忽略系统负荷预测误差，系统一天负荷值见图9,系统一天风电预测和 实测值见图10,通用分布参数见图11。比较使用通用分布将风电作为调度量和传统将风电 作为净负荷两种调度模式的调度结果。
[0163] 4.1、风电调度出力。
[0164] 如图12,系统前Ξ小时风电调度出力与预测出力的关系。根据分析，若风电低高估 成本极大且相等，PDF曲线不偏轴时，则风电的调度量应在使低高估期望相等的点，即CDF = 0.5点，也即PDF最高点，故对于风电分布偏轴的PDF曲线，风电调度量应在尖峰偏0.5p. U. 处。而对于将风电作为净负荷的算法，其PDF曲线相当于只在其调度值的箱内，故将风电作 为调度量的结果应为风电预测值略靠近0.5P.U.。而实际系统运行过程中，低估成本要高于 高估成本，从此因素考虑，风电调度值应略低于风电预测值。
[0165] 综合如上两个因素，可W得到风电调度值与期望值的一个大致的规律：当风电低 于0.5p.U.时，风电调度值在预测值上下；当风电高于0.5p.U.时，风电调度值低于风电预测 值。图12结果证明了运个分析。5-105min时，风电调度值在预测值上下，105min-ISOmin， 风电调度值低于风电预测值。
[0166] 4.2、系统运行成本。
[0167] 如图13,系统前Ξ小时风电作为调度量和风电作为净负荷两种调度模式的成本比 较，可W看出，风电作为调度量相比
于风电作为净负荷，可W大大减少由风电随机性带来的 系统对风电低高估成本，减少总的系统发电成本。
[016引 4.3、调度模型计算时间。
[0169] 如图14,火电机组数量与算法计算时间的关系，可W看出，调度算法在火电机组数 量为7台时，计算时间只有0.097s，且计算时间和火电机组数量呈线性关系（区别于二次规 划法，计算时间和火电机组数量的平方呈线性关系）。算法计算时间满足实时调度的要求。
[0170] 如图15,通用分布和正态分布对于表征风功率，并将风功率作为调度量进行调度 算法计算时间。正态分布也可W用来描述风电功率的随机性，但正态分布相比于通用分布 的拟合效果较差。除此之外，虽然正态分布通过查表也可计算风电成本等微增率，但由于每 次迭代计算中通用分布的解析表达式计算时间比正态分布的查表计算时间快几十倍，且通 用分布由于有解析表达式，可W进行线性化减少迭代计算次数。故由于运两个方面，使用基 于通用分布的含风电经济调度计算时间远远小于基于正态分布的含风电经济调度计算时 间。
[0171] 4.4、调度模型对于保证火电机组发电计划稳定性的分析。
[0172] 在实时调度中，每个计算周期(5min)计算下后面一小时的发电计划，即后面12个 调度周期（12X5min)的发电计划，但只使用下一个调度周期的发电计划。故每个调度周期 的发电计划被重复计算了 12次才被使用，假设现在是Omin，需要计算5、10…60min运12个调 度周期的发电计划。第60min的发电计划，在0、5···55π?η的计算周期均被计算，但只使用 55min计算周期的发电计划结果。在运12个计算周期计算的第60min的计划，必须有一定的 稳定性，尽量减少不必要的发电计划更改，否则会使发电机无所适从。
[0173] 见图16,比较将风电作为净负荷进行调度的本文方法和经典的二次规划算法(将 风电作为净负荷处理），对于火电机组调度出力稳定性的比较。可W看出，9:35制定的9:40 到l〇:35的发电计划，本文方法和二次规划法的对于相同的一台发电机其处理曲线相同。在 9:40制定9:45到10:40的发电计划，本文方法对于9:45到10:35的发电计划和9:35制定的是 相同的，即没有出现发电计划变动；而使用二次规划算法制定的发电计划在10:00到10:35 发生了变动。在9:45,二次规划法制定的发电计划也出现了变动，而本文方法在Ξ个计算时 间对于相同调度周期的出力结果是相同的，即本文算法相比于传统算法能有效保证发电计 划的稳定性。
[0174] 4.5、采用不同方式表征风功率的随机性对系统频率的影响比较。
[0175] 分别计算Ξ种调度模式：1)将风电作为净负荷处理（即不进行调度），2)使用正态 分布表征风电功率并进行调度，3)使用通用分布表征风电功率并进行调度。Ξ种调度模式 均在相同的负荷值(见图9)、相同的风功率预测值和实测值(见图10)下进行一天的滚动计 算，记录系统频率偏差超过某一设定值的次数，用W表征系统有功平衡情况。
[0176] 见图17结果：1)将风电作为净负荷系统频率偏差最多，2)使用正态分布表征风电 功率次之，3)使用通用分布表征风电功率最少。经过之前分析，由于风电预测的特性，风电 作为调度量的结果应为风电预测值略靠近0.5P.U.，风电随机性对于系统有功平衡的影响 最小。1)将风电作为净负荷，即调度量结果就是预测值，故最不利于系统有功平衡;2)使用 正态分布表征风电功率有使风电预测值略靠近0.5p.U.的效果，但由于正态分布为对称分 布，不具有通用分布的偏轴特性，所W其对系统有功平衡的保证能力强于将风电作为净负 荷，但不如3)使用通用分布表征风电功率。
【主权项】
1. 一种基于通用分布的含风电电力系统实时动态经济调度方法，其特征在于，包括以 下步骤： 步骤1:输入日前系统负荷预测数据和日前风电功率预测数据、火电机组运行参数、历 史统计数据，所述历史统计数据包括不同风电功率预测水平下实际风电功率的通用分布参 数α，β， γ ; 步骤2:定义pi, t为第i台火电机组t时刻的出力，火电机组的总数为I，i = I，2，…，I，wj, t 为第j个风电场t时刻的计划出力，风电场的总数为J，j = I，2，...，J，T为时刻的总数，t = I， 2,...,T； 计算从当前时刻也就是当前计算周期到后面1小时的调度时间域每个调度周期火电机 组出力可能达到的范围，称为火电机组在t调度周期的可行域，该1小时的调度时间域为12 个调度周期：式中：P it表示t时刻第i台火电机组的出力;三it和5*是表示t时刻第i台火电机组的出力 下限和上限； 使用通用分布表征风功率，使用等微增率准则进行12个调度周期的风电火电静态经济 调度，得到风电和火电的经济调度曲线，基于以下的目标函数以及约束条件：火电机组出力功率约束条件：己< Λ., </飞式三 火电机组爬坡约束条件:U Pit-Pi(t+1) < Si式四 风电场出力功率约束条件：系统功率平衡约束条件：式中：6，匕^，仏*分别表示七时刻系统总发电成本、火电发电成本和风电成本;￥^表示七 时刻第j个风电场的出力；I表示火电机组数量，J表示风电场数量；P/*表示t时刻系统总的 负荷预测值；再和￡1表示第i台火电机组的出力上下限J 1表示第i台火电机组的爬坡能力； ^7表示第j个风电场的风电场装机容量； 其中，所述火电燃料成本Cg, 1; t采用以下方式计算，其中，Plt为第i台火电机组t时刻的出力，为第i台火电机组的燃料成本系数； 所述风电成本Cw, j,t采用以下方式计算， Cw, j, t(wj, t) = djWjt+Cun, j, t(wj, t)+C〇v, j, t(wj, t)式八 其中，Wjt为第j个风电场t时刻的计划出力，dj为第j个风电场的运行成本系数； 所述平均低估成本Cun^t采用风电场的平均弃风成本，采用以下方式计算，其中，kun,j为第j个风电场对应的低估成本系数，Wav^t为第j个风电场t时刻的实际可能 出力，A(Wav^t)为第j个风电场在对应风电预测水平下实际可能出力的概率密度函数，表 达形式为通用分布的概率密度函数，根据t时刻对应通用分布参数α，β，γ确定； 所述平均高估成本Cciv^t采用系统的平均备用成本，采用以下方式计算，其中，IWd为第j个风电场对应的高估成本系数，Wav,it为第j个风电场t时刻的实际可能 出力，A(Wav^t)为第j个风电场在对应风电预测水平下实际可能出力的概率密度函数，表 达形式为通用分布的概率密度函数，根据t时刻对应通用分布参数α，β，γ确定； 步骤3:利用步骤2得到的12个调度周期各自的风电场出力，修正火电机组各个调度周 期的可行域;在新的火电机组可行域内，使用等微增率准则对火电机组的出力进行修正，保 证火电机组爬坡限制并保证火电机组发电计划的稳定性； 步骤4:利用步骤2得到的风电场出力和步骤4得到的火电机组出力，输出火电机组和风 电场的计划出力； 而且，所述火电燃料成本Cg, 1; t采用以下方式计算，其中，Plt为第i台火电机组t时刻的出力，为第i台火电机组的燃料成本系数； 而且，所述风电成本t采用以下方式计算， Cw, j, t(wj, t) = djWjt+Cun, j, t(wj, t)+C〇v, j, t(wj, t)式十-· 其中，Wjt为第j个风电场t时刻的计划出力，dj为第j个风电场的运行成本系数； 而且，所述平均低估成本Cun，j，t采用风电场的平均弃风成本，采用以下方式计算，其中，kUy为第j个风电场对应的低估成本系数，Wav,it为第j个风电场t时刻的实际可能 出力，A(Wav^t)为第j个风电场在对应风电预测水平下实际可能出力的概率密度函数，表 达形式为通用分布的概率密度函数，根据t时刻对应通用分布参数α，β，γ确定； 而且，所述平均高估成本Cciv^t采用系统的平均备用成本，采用以下方式计算，其中，IWd为第j个风电场对应的高估成本系数，Wav,it为第j个风电场t时刻的实际可能 出力，A(Wav^t)为第j个风电场在对应风电预测水平下实际可能出力的概率密度函数，表 达形式为通用分布的概率密度函数，根据t时刻对应通用分布参数α，β，γ确定。
【专利摘要】本发明公开了一种基于通用分布的含风电电力系统实时动态经济调度方法，通过统计历史数据，离线生成通用分布的形状参数表，表征在不同风电功率预测值下的风电功率概率分布；输入12个调度周期系统负荷预测数据和风电功率预测数据等；计算火电机组在实时调度12个调度周期内的火电出力可行域；使用通用分布表征风功率，使用等微增率准则进行12个时间周期的风电火电静态经济调度，得到风电和火电的静态经济调度曲线，将随机经济调度转化为确定性经济调度模型；使用等微增率准则对火电机组的出力进行修正，在保证火电机组爬坡限制并保证火电机组发电计划的稳定性的基础上，使火电出力成本最优。以某实际风电功率为基础的仿真算例验证了所提方法的可行性和合理性。本方法具有良好的推广价值和应用前景。
【IPC分类】G06Q10/04, H02J3/46, G06Q50/06, H02J3/38
【公开号】CN105490310
【申请号】CN201610049466
【发明人】徐箭, 唐程辉, 孙元章, 江海燕, 王豹, 洪敏
【申请人】武汉大学
【公开日】2016年4月13日
【申请日】2016年1月25日
