[发明专利]基于风/光出力柯西分布的火电机组动态实时调度方法

摘要

本发明涉及一种基于风/光出力柯西分布的火电机组动态实时调度方法，属于电力系统的运行技术领域。本通过分析风电出力的历史数据，利用统计或拟合软件进行联合柯西分布拟合。针对确定的电力系统参数，建立机会约束的随机动态实时调度模型；然后，利用柯西分布的数学性质，把原问题转化为容易求解的线性约束凸优化问题；最后求解调度模型，得到调度策略。本发明充分利用了柯西分布在风电/光伏电站出力的短期预测方面的优越性和优良的数学特性，有效提高了模型的求解效率，同时风险水平可调的机会约束模型消除了传统鲁棒经济调度的保守性，为决策者提供了更合理的调度依据。本发明可应用于包含大规模风电并网的电力系统有功实时经济调度中。

主权项

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维列向量，那么随机变量服从位置参数和尺度参数为(a^Tμ_t,a^TΣ_ta)的一维柯西分布，其概率密度函数、累积分布函数以及累积分布函数的逆函数可以表达为以下形式：其中，概率密度函数为：累积分布函数为：累积分布函数的逆函数为：其中，tan为正切函数，arctan为反正切函数，F为分位数；(2)建立一个基于风电/光伏电站出力柯西分布的随机动态实时调度模型，该随机动态实时调度模型由目标函数和约束条件组成，具体步骤如下：(2‑1)建立基于风电/光伏电站出力柯西分布的随机动态实时调度模型的目标函数f：目标函数为使运行成本最小化，表达式如下：其中，T、N和J分别表示调度时段t的数量、火电机组数量以及发电量自动控制机组的数量，t、i和j分别为调度时段、火电机组编号和发电量自动控制机组的编号，上标s表示该变量描述火电机组，上标“+”表示该变量描述正旋转备用，上标“‑”表示该变量描述负旋转备用，表示第i台火电机组在t调度时段的计划出力，表示第j个发电量自动控制机组在t调度时段内计划出力，表示t调度时段所有风电/光伏电站实际出力的总和，和分别表示火电机组和发电量自动控制机组的燃料成本：其中，a_i,t,b_i,t,c_i,t分别为t时段火电机组i的燃料成本的二次项系数、一次项系数和常数项，a_j,t,b_j,t,c_j,t分别为t时段发电量自动控制机组j的燃料成本的二次项系数、一次项系数和常数项，E(·)表示随机变量的期望值，表示在t调度时段内，由于风电/光伏电站发电实际出力低于计划出力引起的正旋转备用的需求成本，即高估风电/光伏电站出力的惩罚，当风电/光伏电站的实际出力小于计划值时会调度发电量自动控制机组的正旋转备用以维持功率平衡，具体表达式如下：其中，为正旋转备用的成本系数，α_j为第j个发电量自动控制机组的功率分配系数，根据发电量自动控制机组的额定容量占发电量自动控制机组的总容量的比例确定，w_t表示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表示相邻两个调度时段之间的调度间隔，β为调度员设定的可接受风险水平，表示w_t,t‑1的累积分布函数的逆函数，的表达式如下：(2‑2‑4)发电量自动控制机组的旋转备用约束，具体表达式如下：对其中，和分别表示t调度时段第j台发电量自动控制机组提供的正负旋转备用的数量，和分别表示t调度时段电网所需要的最小正、负旋转备用的数量，ε为调度员设定的可接受的风险水平；(2‑2‑5)电网线路潮流约束，表达式如下：对其中，G_i,l为电网中第l条线路对第i台火力发电机组有功出力的转移分布因子，G_j,l为第l条线路对第j个发电量自动控制机组有功出力的转移分布因子，G_k,l为第l条线路对第k个风电/光伏电站有功出力的转移分布因子，G_d,l为第l条线路对第d个节点负荷功率的转移分布因子，各转移分布因子分别从电网调度中心获取，L_l,t为t调度时段第l条线路上的有功功率上限，η为电网线路上的有功功率超过线路有功功率上限的风险水平，由调度员设定，为的累积分布函数的逆函数，满足：α_l为K维向量，其第k个元素为(3)采用内点法，求解上述步骤(2)中的目标函数和约束条件组成的随机动态实时调度模型，得到和将其中的作为t调度时段第i个火电机组的计划出力，第j个发电量自动控制机组的计划出力，作为第k个风电/光伏电站的参考出力，实现基于风/光出力柯西分布的火电机组动态实时调度。