[发明专利]一种计及变量相关性的概率潮流计算方法

摘要

本发明提供一种计及变量相关性的概率潮流计算方法，所述方法包括(1)获取常规潮流计算数据，确定性潮流计算；(2)求取PV节点和PQ节点功率注入的变化对节点电压影响的计算模型；(3)采用正态分布的多项式表示非正态分布的随机变量；(4)采用Cholesky分解将相关的随机变量转化为不相关的随机变量；(5)对变量进行概率潮流计算。本发明在保证计算精度的同时提高了计算速度，能够适用于大规模间歇性能源并网后的复杂电网系统分析、安全评估等方面，有助于提升电力系统的新能源接纳能力。

主权项

1.一种计及变量相关性的概率潮流计算方法，其特征在于，所述方法步骤如下：(1)获取常规潮流计算数据；(2)潮流计算；(3)求取PV节点功率注入的变化对节点电压影响的计算模型；所述潮流计算模型可表示如下：式中：ΔP_G为PV节点有功注入变化向量，ΔV_Gr为PV节点实部电压变化量，ΔV_Lr为PQ节点实部电压变化量，ΔV_Lm为PQ节点虚部电压变化量，计算模型的系数矩阵[H^/]由和四部分组成；ΔPG的求取方法如下:式中，ΔPk为k号节点有功注入变化量，Vmk为k号节点的虚部电压；的求取方法如下：式中，h_ii为矩阵的i行i列元素，h_ij为矩阵的i行j列元素,b_ii为矩阵的i行i列元素的虚部，g_ij为矩阵的i行j列元素的实部，b_ij为矩阵的i行j列的虚部；的求取方法如下：式中，g_kn为矩阵的k行n列元素的实部，b_kn为矩阵的k行n列元素的虚部，V_rk为k号节点的实部电压，V_mk为k号节点的虚部电压；的求取方法如下：式中，g_kn为矩阵的k行n列元素的实部，b_kn为矩阵的k行n列元素的虚部，V_rk为k号节点的实部电压，V_mk为k号节点的虚部电压；的求取方法如下：当i＝j时：式中，g_ii为矩阵的i行i列元素的实部，b_ii为矩阵的i行i列元素的虚部，V_ri为i号节点的实部电压，V_mi为i号节点的虚部电压，I_ri为i号节点的实部注入电流，I_mi为i号节点的虚部注入电流，A_ii、为矩阵的系数；当i≠j时：(4)求取PQ节点功率注入的变化对节点电压影响的计算模型；所述PQ节点功率注入的变化对节点电压影响的计算模型，其定义与求取方式与常规潮流计算完全一致；(5)采用正态分布的多项式表示非正态分布的随机变量；(6)将原随机变量的相关性系数矩阵转化为正态分布随机变量间的系数矩阵；(7)对相关性系数矩阵进行Cholesky分解；(8)修正对PV节点注入功率变化对节点电压影响的潮流计算模型；(9)修正对PQ节点注入功率变化对节点电压影响的潮流计算模型；(10)求取灵敏度矩阵；(11)求取节点电压的半不变量；(12)求取节点电压的各阶中心距；(13)求取节点电压的概率分布；设系统模型中节点1,2,...,k号节点为PV节点，k+1,k+2,...,n号节点为PQ节点；所述步骤(12)由节点电压实部和虚部的各阶半不变量求出各阶中心距。