[发明专利]一种基于随机鲁棒优化的电力系统调度方法

摘要

本发明涉及一种基于随机鲁棒优化的电力系统调度方法，属于电力系统运行控制技术领域。该方法首先建立由目标函数和约束条件构成的电力系统调度模型；对约束条件进行转化，将原模型转化为一个非线性优化问题；对转化后的模型求解，得到发电机出力，并根据求解结果进行相应调度决策实现电力系统调度。本发明考虑了电力系统中可再生能源的波动，减小了校正控制的成本，保证了电力系统安全，适合应用于电力系统的滚动调度等场景之中。

主权项

1.一种基于随机鲁棒优化的电力系统调度方法，其特征在于包括以下步骤：(1)利用时间序列分析方法对电力系统的功率进行预测，分别得到在电力系统的所有节点处风电的预测功率上限和下限w，将电力系统的第i个节点处风电的预测功率上限记为第i个节点处风电的预测功率下限记为w_i，记电力系统中风电的实际功率为w，w为一个多维不确定性变量；(2)建立基于随机鲁棒优化的电力系统调度模型，具体过程包括以下步骤：(2‑1)确定电力系统调度模型的目标函数：电力系统调度的目标为总成本最小化，总成本是发电成本、弃风惩罚和后续的校正控制成本之和，其中：发电成本函数的表达式如下：上式中，G_i表示电力系统第i个节点处发电机的发电成本函数，P_Gi表示电力系统第i个节点处发电机的有功功率，是待求解量，表示接入发电机的电力系统节点，从系统所有发电机的接入位置获取，是已知量；第i个节点处发电机的发电成本函数G_i(P_Gi)的表达式如下：G_i(P_Gi)＝a_iP_Gi²+b_iP_Gi+c_i其中a_i,b_i,c_i代表第i个节点处发电机的发电成本参数，a_i,b_i,c_i由发电厂根据发电成本给出，为已知量；弃风惩罚的表达式如下：上式中，w_i,max,w_i,min分别表示电力系统第i个节点处风电的允许出力最大值和允许出力最小值，均为待求解量，表示接入风电的电力系统节点集合，从系统所有风电的接入位置信息获取，是已知量；第i个节点处风电的弃风惩罚函数C_i(w_i,max,w_i,min)的表达式如下：其中，k_i代表第i个节点处风电的弃风惩罚系数，k_i是已知量，由电力系统控制中心的调度系统给出，w_i代表第i个节点处风电的实际功率，为一个随机变量，p(w_i)代表第i个节点处风电的实际功率w_i的概率密度函数，p(w_i)为已知量，由风电场的概率预测得到；第i个节点处发电机的校正控制成本的表达式如下：其中，UR_i,DR_i分别代表第i个节点处发电机的上调控制成本和下调控制成本；第i个节点处发电机的上调控制成本UR_i表达式如下：其中，g_i表示第i个节点处发电机的上调成本系数，由发电厂根据发电机调节能力给出，是已知量，α_i表示第i个节点处发电机的自动发电控制调节系数，由发电厂根据发电机调节能力给出，是已知量，u₀表示所有风电的有功功率设定值总和，u₀是待求解量，u_min表示所有风电的允许出力最小值总和，u_min是待求解量，u表示所有风电的实际功率总和，为一随机变量，p(u)代表所有风电的实际功率总和u的概率密度函数，p(u)是已知量，由风电场的概率预测得到；第i个节点处发电机的下调控制成本DR_i表达式如下：其中，h_i表示第i个节点处发电机的下调成本系数，由发电厂根据发电机调节能力给出，h_i是已知量，u_max表示所有风电的允许出力最大值总和，是待求解量，表示所有风电的预测功率上限总和，通过对步骤(1)中电力系统的所有节点处风电的预测功率上限求和得到；(2‑2)确定电力系统调度模型的约束条件，包括：(2‑2‑1)发电机和风电的功率约束：P_Gi,min≤P_Gi‑α_i(u_max‑u₀)P_Gi‑α_i(u_min‑u₀)≤P_Gi,maxw_i,min≤w0≤w_i,min≤w_i,0≤w_i,max其中，w_i,0表示第i个节点处风电的功率设定值，是待求解量，P_Gi,min,P_Gi,max分别表示第i个节点处发电机的功率最小值和功率最大，是已知量，由发电厂根据发电机的发电能力给出；(2‑2‑2)电力系统的功率平衡方程约束：其中，表示接入负载的电力系统节点集合，P_Di表示电力系统第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)拆分成两部分，包括f₁(x)和f₂(x)，其中f₁(x)代表步骤(2‑1)中的发电成本，f₁(x)为二次函数，f₂(x)代表步骤(2‑1)中的弃风惩罚和校正控制成本之和，f₂(x)为非线性函数；(4‑2)采用二次规划求解方法，通过下式求解得到最优解x₀：Ax≤b(4‑3)采用二次规划求解方法，通过下式求解得到最优解为x₁：Ax≤b其中，表示f₂(x)函数在x₀的导数向量；(4‑4)设定一个迭代允许误差ε，对步骤(4‑2)和步骤(4‑3)的计算结果进行判断，若abs(f(x₁)‑f(x₀))＜ε，则判定x₁为电力系统调度最优解，停止迭代，若abs(f(x₁)‑f(x₀))≥ε，则将x₁的值代替x₀，返回步骤(4‑3)，其中abs表示取绝对值，；(5)将与步骤(4)得到的电力系统调度最优解中相对应的发电机有功功率作为电力系统各发电机的发电量，制定发电机的发电计划，同时将与步骤(4)得到的电力系统调度最优解中相对应的所有风电的允许出力最小值和允许出力最大值下发给风电场，制定风电场的发电计划，完成电力系统的调度。