一种计及新能源发电预测误差的改进随机潮流算法
【技术领域】
[0001] 本发明设及电力系统技术领域，尤其是一种计及新能源发电预测误差的改进随机 潮流算法。
【背景技术】
[0002] 随机潮流(PLF)能反映电力系统中各种因素的随机变化对系统运行的影响。它可 W综合考虑风电出力的随机波动、负荷的变化、发电机的强迫停运及线路的故障等各种不 确定情况，得到系统节点电压和支路潮流的概率统计特性。相比于常规潮流，它大大减少了 计算量。对计算结果做适当处理，可推出线路过负荷概率、电压越限概率等系统静态安全的 评价指标，W便及时发现电网中的薄弱环节，为电力系统规划和决策人员提供有价值的信 息。
[0003] 蒙特卡罗仿真法(Monte Carlo Simulation,MCS)是求解随机潮流问题最直观、最 基础的算法。它是一种通过随机变量的变化规律或者物理现象本身的统计规律求取数学、 物理、工程技术问题近似解的数值方法。具有适用性广、计算精度高且原理简单的优势。然 而计算规模大、耗时长限制了其实际应用。非线性蒙特卡罗仿真法一般都是作为基准方法 与其他算法进行比较，W验证算法的准确性。
[0004] 与非线性蒙特卡罗法相比，线性蒙特卡罗法主要在其基础上增加了对非线性潮流 方程的线性化处理，避免了多次迭代计算，从而可W显著减少计算量，提高运算速度。然而 当节点注入量变化范围较大时，该方法可能存在较大的误差。
【发明内容】
[0005] 发明目的：针对上述现有技术存在的缺陷，本发明旨在提供一种计及新能源发电 预测误差的改进随机潮流算法。
[0006] 技术方案:一种计及新能源发电预测误差的改进随机潮流算法，包括如下步骤：
[0007] S1:输入电网原始数据及风电场相关数据；
[000引S2:根据S1中的输入数据，建立风速的We化ull分布模型、负荷的正态分布模型W 及常规发电机组的0-1分布模型，建立新能源预测误差随机模型；
[0009] S3:采用蒙特卡罗法进行抽样;根据风速的Weibul 1分布模型、负荷的正态分布模 型及常规发电机组的0-1分布模型，通过反函数变换法分别产生风速、负荷和发电机组出力 的随机样本值；
[0010] S4:根据尾流效应建立风电场综合模型，得到不同位置处风机的风速并由风电场 功率输出特性求得单机输出功率;最后将所有风机功率叠加，得到整个风电场的输出功率；
[0011] S5:用牛顿拉夫逊法进行确定性潮流计算;求出节点电压和支路潮流的期望值W 及灵敏度矩阵So和To;
[0012] S6:求出状态变量和支路潮流的随机扰动量ΔΧ和ΔΖ，再得到最终的电压状态向 量X和支路潮流的功率向量Z;
[0013] S7:判断抽样次数是否达到5000;若小于5000,返回步骤S6,否则执行步骤S8;
[0014] S8:统计得到节点电压及支路潮流的概率分布，输出结果。
[0015] 进一步的，步骤S1中所述电网原始数据及风电场相关数据包括常规潮流计算数 据、负荷分布情况、常规发电机的出力、强迫停运率、风电场历史风速和风电场功率输出特 性。
[0016] 进一步的，步骤S2中所述建立新能源预测误差随机模型包括如下子步骤：
[0017] S2.1:用正态分布分别近似描述有功和无功功率的概率密度函数：
[001 引
W
[0019] 式中，μρ、σρ和yQ、〇Q分别为有功功率的期望值和方差和无功功率的期望值和方差；
[0020] S2.2:用双参数Weibul 1函数来拟合风速，其概率密度函数表示如下：
[0021]
(2)
[0022] 式中，V为风速;k和C分别为Weibull分布的两个参数;将式（2)产生的风速序列代 入式(3)的风机输出特性表达式，求出风力发电机有功出力的概率分布：
[0023]
(3)
[0024] 式中，Pr为风机的额定功率;VR、Vin和VDut分别为额定风速、切入风速和切出风速；
[0025] S2.3:假设太阳福射强度和气溫均服从正态分布；W福射强度GsTc = 1000W/m2和电 池板表面溫度Tr=化°C作为标准状态，则任意状态下PV系统的输出功率如下：
[0026]
(4)
[0027] 式中，PsTG为标准状态PV系统的额定功率;Ging、Gstg为实际和标况下的福射强度;k 为最大功率溫度系数;T。和Tr分别为实际和参考溫度。
[0028] 进一步的，步骤S4中所述根据尾流效应建立风电场综合模型具体为:采用适用于 不同地理环境的尾流模型，综合运用Jensen模型及Lissaman模型，得到大型风电场综合模 型;通过该模型求得位于不同位置的风电机组的输出功率，通过整个风电场中风机输出功 率的叠加，得到整个风电场的功率输出。
[0029] 进一步的，所述位于不同位置的风电机组的输出功率包括位于平坦地形风电机组 的出力和位于复杂地形风电机组的出力，所述位于平坦地形风电机组的出力的计算方法如 下：
[0030]
I其中X是两个风机之间的距离，R和Rw分别为叶轮半径和尾流半径，ν〇、ντ、 VX分别对应自然风速、通过叶轮的风速和受尾流效应影响的风速;通过动量理论计算得到 如下表达式：
[0031]
(5)
[0032] 式中：Ct为风机的推力系数;k为尾流衰减系数。
[0033] 进一步的，所述位于复杂地形风电机组的出力计算方法如下：
[0034] 两台相同型号的风机分别位于两个场地的边缘，它们的所处的坐标位置都是x = 0，风速都为VO;当x = 0处没有风机时，平坦地形X处风速仍为VO，而海拔高度Η处的风速为：
[0035]
(6)
[0036] 式中:h为风机塔筒高度;αι为风速随高度的变化系数，一般％ = ^ ;海拔越高，风速 越大；
[0037] 当在X = 0处安装风机后，受尾流效应影响，X=Χ处风速分别为：
[00：3 引
（7)
[0039] 式中：Vix为地形平坦即风机安装处高度相同时的风机风速，V2X为地形复杂即风机 安装处高度不同时的风机风速;山和d2是相应的风速下降系数；由Jensen模型求得由Jensen 模型求得：
[0040]
(8)
[0041] 假设两种地形的尾流损耗相同，尾流中的压力相同，则根据无损耗伯努利方程得 到：
[0042]
(9)。
[0043] 进一步的，步骤S5中所述用牛顿拉夫逊法进行确定性潮流计算方法如下：
[0044] 用极坐标表示的电力系统交流潮流的节点功率和支路潮流方程分别如式（10)和 (11)所示：
[0047] 式中：Pi和化分别为节点i的有功功率注入量和无功功率注入量;Gu和Βυ分别为节 点导纳矩阵相应元素的实部和虚部;Pij和Qij分别为线路ij中由节点i流向节点j的有功功 率和无功功率;tu为支路变压器变比的标么值;biw为支路容纳的一半;Ui、Uj分别为节点i、 j的电压，θι功节点i、j间的相角差；
[004引可将上式概括为：
[0049] W=f(X) (12)
[0050] Z = g(X) (13)
[0051 ]式中：W为节点功率注入向量，包括有功功率和无功功率;X为节点的状态向量，包 括节点电压的幅值和相角;Z为支路潮流的功率向量;在随机潮流模型中，节点功率注入量 和状态变量都是随机变量，所W有：
[0052] W=Wo+AW (14)
[0053] Χ=Χ〇+ΔΧ (15)
[0054] 其中，Wo、Xo分别为节点功率注入量W和状态变量X的期望值，AW、Δ X为随机扰动， 可W看作