[发明专利]基于改进连续潮流法的配电网PV曲线生成方法

权利要求书

1.一种基于改进连续潮流法的配电网PV曲线生成方法，其特征是所述方法包括：

步骤1：计算当前配电网导纳矩阵、潮流初始解U_i、K_ij和L_ij，并令负荷变化参数λ＝0；

U_i为节点i的电压相量幅值的平方，即

K_ij为节点i的和节点j的电压相量的内积，即K_ij＝V_iV_jcosθ_ij；

L_ij为节点i的和节点j的电压相量的外积，即L_ij＝V_iV_jsinθ_ij；

V_i为节点i的电压相量，V_j为节点j的电压相量；

步骤2：令k＝1，分别将潮流初始解U_i、K_ij和L_ij以及负荷变化参数λ作为潮流已知解和λ^k；

步骤3：根据潮流已知解和λ^k，计算潮流预测解和λ^(k+1)；

步骤4：对潮流预测解进行校正得到潮流实际解和λ^k+1；

步骤5：判断λ^k+1是否小于λ^k，如果λ^k+1小于λ^k，则执行步骤6；否则，令k＝k+1，将潮流实际解作为潮流已知解，返回步骤3；

步骤6：根据潮流实际解生成PV曲线。

2.根据权利要求1所述的方法，其特征是所述计算潮流预测解的过程包括：

子步骤A1：根据潮流已知解和λ^k计算潮流已知解的切相量dUikdKijkdLijkλkT,]]>其计算公式为：

HJNJPλiIMDJQλiRSTJfλiemdUikdKijkdLijkdλk=000±1;]]>

其中，H、I、R、J、M、S、N、D、T、J_Pλi、J_Qλi、J_fλi和e_m均为切相量计算矩阵参数；

H=∂Pi∂Uik=-Gii,]]>I=∂Qi∂Uik=Bii,]]>R=∂fij∂Uik=-Ujk,]]>J=∂Pi∂Kijk=-Gij,]]>M=∂Qi∂Kijk=Bij,]]>S=∂fij∂Kijk=2Kijk,]]>N=∂Pi∂Lijk=-Bij,]]>D=∂Qi∂Lijk=-Gij,]]>T=∂fij∂Lijk=2Lijk,]]>J_Pλi＝-P_iL，J_Qλi＝-Q_iL，J_fλi＝0；

和为潮流已知解，G_ii、G_ij为由节点导纳矩阵各元素的实部组成的矩阵，B_ii、B_ij为由节点导纳矩阵各元素的虚部组成的矩阵，P_iL为节点i的有功负荷，Q_iL为节点i的无功负荷，Pi=PiG-PiL(1+λk)-UikGii-Σj∈i(GijKijk+BijLijk),]]>Qi=-QiG+QiL(1+λk)-UikBii+Σj∈i(GijLijk-BijKijk),]]>fij=(Kijk)2+(Lijk)2-UikUjk,]]>为对与节点i相连的节点求和；

当k＝1时，e_m为最后一个元素为1且其余元素为0的行向量；当k≠1时，e_m为第m个元素为1且其余元素为0的行向量；其中，m为中的最大值所对应的位置；

子步骤A2：根据切相量计算潮流预测解，计算公式为：

Ui(k+1)=Uik+σdUikKij(k+1)=Kijk+σdKijkLij(k+1)=Lijk+σdLijkλ(k+1)=λk+σdλk;]]>σ为设定值。

3.根据权利要求2所述的方法，其特征是所述对潮流预测解进行校正得到潮流实际解包括如下子步骤：

子步骤B1：令t＝1；

子步骤B2：判断P(Ut,Kt,Lt,λt)Q(Ut,Kt,Lt,λt)f(Ut,Kt,Lt,λt)<ϵ]]>是否成立，如果P(Ut,Kt,Lt,λt)Q(Ut,Kt,Lt,λt)f(Ut,Kt,Lt,λt)<ϵ,]]>则执行子步骤B6；否则，执行子步骤B3；

P(Ut,Kt,Lt,λt)Q(Ut,Kt,Lt,λt)f(Ut,Kt,Lt,λt)<ϵ]]>是指列向量P(Ut,Kt,Lt,λt)Q(Ut,Kt,Lt,λt)f(Ut,Kt,Lt,λt)]]>中的所有元素的值都小于ε，ε为设定阈值；

P(Ut,Kt,Lt,λt)=PiG-PiL(1+λt)-UitGii-Σj∈i(GijKijt+BijLijt);]]>

Q(Ut,Kt,Lt,λt)=QiG-QiL(1+λt)+UitBii-Σj∈i(GijLijt-BijKijt);]]>

f(Ut,Kt,Lt,λt)=(Kijt)2+(Lijt)2-UitUjt;]]>

和λ^t为潮流预测解G_ii、G_ij为由节点导纳矩阵各元素的实部组成的矩阵，B_ii、B_ij为由节点导纳矩阵各元素的虚部组成的矩阵，P_iG为节点i的有功出力，P_iL为节点i的有功负荷，Q_iG为节点i的无功出力，Q_iL为节点i的无功负荷，为对与节点i相连的节点求和；

子步骤B3：根据公式HJNJPλiIMDJQλiRtStTtJfλitemΔUitΔKijtΔLijtΔλt=ΔPitΔQitΔfit0]]>求解修正量和Δλ^t；

ΔPit=-PiG+PiL(1+λt)+UitGii+Σj∈i(GijKijt+BijLijt);]]>

ΔQit=-QiG+QiL(1+λt)-UitBii+Σj∈i(GijLijt-BijKijt);]]>

Δfit=UitUjt-(Kijt)2-(Lijt)2;]]>

H、I、R^t、J、M、S^t、N、D、T^t、J_Pλi、J_Qλi、和e_m均为修正量计算矩阵参数；

H=∂Pi∂Uit=-Gii,]]>I=∂Qi∂Uit=Bii,]]>Rt=∂fij∂Uit=-Ujt,]]>J=∂Pi∂Kijt=-Gij,]]>M=∂Qi∂Kijt=Bij,]]>St=∂fij∂Kijt=2Kijt,]]>N=∂Pi∂Lijt=-Bij,]]>D=∂Qi∂Lijt=-Gij,]]>Tt=∂fij∂Lijt=2Lijt,]]>J_Pλi＝-P_iL，J_Qλi＝-Q_iL，Jfλit=0;]]>

Pi=PiG-PiL(1+λt)-UitGii-Σj∈i(GijKijt+BijLijt);]]>

Qi=-QiG+QiL(1+λt)-UitBii+Σj∈i(GijLijt-BijKijt);]]>

fij=(Kijt)2+(Lijt)2-UitUjt;]]>

当t＝1时，e_m为最后一个元素为1且其余元素为0的行向量；当t≠1时，e_m为第m个元素为1且其余元素为0的行向量；其中，m为中的最大值所对应的位置；

子步骤B4：根据公式Uit+1=Uit+ΔUitKijt+1=Kijt+ΔKijtLijt+1=Lijt+ΔLijtλk+1=λt+Δλt]]>计算潮流实际解；

子步骤B5：令t＝t+1，返回子步骤B2；

子步骤B6：结束。