[发明专利]一种基于最小输入误差准则的含风电系统的概率潮流计算方法

摘要

本发明公开了一种基于最小输入误差准则的含风电系统概率潮流计算方法，包括如下步骤步骤1获取原始数据；步骤2对不确定因素进行概率建模，获得参考分布函数；步骤3基于最小输入误差准则，判断最优计算阶数kopt；步骤4，采用kopt阶半不变量与Gram‑Charlier级数展开计算概率潮流。本发明能快速确定使输入误差最小的级数阶数，全面考虑各不确定因素对分布函数拟合效果的影响，从而在保证计算精度的前提下，提高含风电系统概率潮流的计算效率。

主权项

一种基于最小输入误差准则的含风电系统的概率潮流计算方法，所述含风电系统包括火电机组、风电场和负荷；将所述火电机组记为G＝{G1,G2,...,Gi,...,GNG}，Gi表示第i台火电机组，所述第i台火电机组Gi的有功功率与无功功率分别记为PGi和QGi，i＝1,2,…,NG；将所述风电场记为W＝{W1,W2,...,Wj,...,WNW}；Wj表示第j个风电场；所述第j个风电场Wj的有功功率与无功功率分别记为PWj和QWj，j＝1,2,…,NW；将所述负荷记为L＝{L1,L2,...,Lm,...,LNL}，Lm表示第m个负荷，所述第m个负荷Lm的有功功率与无功功率分别记为PLm和QLm，m＝1,2,…,NL；其特征在于，所述概率潮流计算方法是按如下步骤进行：步骤1：分别建立所述第i台火电机组Gi的有功功率PGi的参考分布函数CDFPGi、所述第i台火电机组Gi的无功功率QGi的参考分布函数CDFQGi、所述第j个风电场Wj的有功功率PWj的参考分布函数CDFPWj、所述第j个风电场Wj的无功功率QWj的参考分布函数CDFQWj、所述第m个负荷Lm的有功功率PLm的参考分布函数CDFPLm、以及所述第m个负荷Lm的无功功率QLm的参考分布函数CDFQLm；步骤2：分别计算所述第i台火电机组Gi的有功功率PGi的1‑kmax阶原点矩与1‑kmax阶半不变量所述第i台火电机组Gi的无功功率QGi的1‑kmax阶原点矩与1‑kmax阶半不变量所述第j个风电场Wj的有功功率PWj的1‑kmax阶原点矩与1‑kmax阶半不变量所述第j个风电场Wj的无功功率QWj的1‑kmax阶原点矩与1‑kmax阶半不变量所述第m个负荷Lm的有功功率PLm的1‑kmax阶原点矩与1‑kmax阶半不变量以及所述第m个负荷Lm的无功功率QLm的1‑kmax阶原点矩与1‑kmax阶半不变量步骤3：采用最小输入误差准则，寻找使拟合误差最小的Gram‑Charlier级数最优阶数kopt；步骤3‑1：定义测试阶数为ktest；1≤ktest≤kmax；并初始化ktest＝3；步骤3‑2：分别提取所述第i台火电机组Gi的有功功率PGi的1‑ktest阶半不变量所述第i台火电机组Gi的无功功率QGi的1‑ktest阶半不变量所述第j个风电场Wj的有功功率PWj的1‑ktest阶半不变量所述第j个风电场Wj的无功功率QWj的1‑ktest阶半不变量所述第m个负荷Lm的有功功率PLm的1‑ktest阶半不变量以及所述第m个负荷Lm的有功功率PLm的1‑ktest阶半不变量并分别运用ktest阶A型Gram‑Charlier级数生成所述第i台火电机组Gi的有功功率PGi的ktest阶拟合分布函数所述第i台火电机组Gi的无功功率QGi的ktest阶拟合分布函数所述第j个风电场Wj的有功功率PWj的ktest阶拟合分布函数所述第j个风电场Wj的无功功率QWj的ktest阶拟合分布函数所述第m个负荷Lm的有功功率PLm的ktest阶拟合分布函数所述第m个负荷Lm的无功功率QLm的ktest阶拟合分布函数步骤3‑3：分别计算所述第i台火电机组Gi的有功功率PGi的ktest阶拟合分布函数的拟合误差指标所述第i台火电机组Gi的无功功率QGi的ktest阶拟合分布函数的拟合误差指标所述第j个风电场Wj的有功功率PWj的ktest阶拟合分布函数的拟合误差指标所述第j个风电场Wj的无功功率QWj的ktest阶拟合分布函数的拟合误差指标所述第m个负荷Lm的有功功率PLm的ktest阶拟合分布函数的拟合误差指标以及所述第m个负荷Lm的无功功率QLm的ktest阶拟合分布函数的拟合误差指标步骤3‑4：计算ktest阶拟合分布函数的总体拟合误差指标并将ktest+1赋值给ktest，判断ktest≤kmax是否成立，若成立，则返回步骤3‑2顺序执行；否则，表示获得总体分布函数拟合误差指标进入步骤3‑5；步骤3‑5：从所述总体分布函数拟合误差指标中选出最小值所对应的阶数作为最优阶数，记为kopt；当存在多个相同的最小值，则选择多个最小值所对应的阶数中最低阶数作为最优阶数，记为kopt；步骤4：以所述第i台火电机组Gi的有功功率PGi的1‑kopt阶半不变量所述第i台火电机组Gi的无功功率QGi的1‑kopt阶半不变量所述第j个风电场Wj的有功功率PWj的1‑kopt阶半不变量所述第j个风电场Wj的无功功率QWj的1‑kopt阶半不变量所述第m个负荷Lm的有功功率PLm的1‑kopt阶半不变量以及所述第m个负荷Lm的无功功率QLm的1‑kopt阶半不变量作为概率潮流计算模型的输入变量并进行概率潮流计算，分别获得节点电压的1‑kopt阶半不变量与支路功率的1‑kopt阶半不变量；再利用kopt阶A型Gram‑Charlier级数拟合所述节点电压与支路功率的分布函数。