[发明专利]一种谐波污染用户的概率潮流评估方法

权利要求书

1.一种谐波污染用户的概率潮流评估方法，其特征在于包括下步骤：

步骤S01：建立节点注入功率和谐波频谱特性概率模型，生成样本集P~={Pij,i=2,...,n,j=1,2,...,N},Q~={Qij,i=2,...,N,j=1,2,...,N},K~{Kih,i=1,2,...,N},]]>

其中，分别表示节点注入有功功率和无功功率的样本集，P_ij、Q_ij分别表示第i个节点的第j个样本；表示谐波频谱的样本集，表示第i个节点的第h次谐波含有率样本；N为样本数，n为节点数，记平衡节点为1；

步骤S02：从所述样本集中取一组计算数据

步骤S03：采用P-Q解耦进行谐波潮流计算：

ΔPΔQ=H00LΔθΔU---(1)]]>

式中，ΔU为节点电压幅值偏差量矩阵，Δθ为节点电压相角偏差量矩阵。

解之可得式(2)；

Δθ=H^-1ΔP(2)

ΔU=L^-1ΔQ

上式中的功率偏差量通过式(3)考虑谐波功率对基波功率的修正量；

ΔP=PiS-Pi1-Σk=2hPikΔQ=QiS-Qi1-Σk=2hQik---(3)]]>

式(3)中，表示节点i的注入有功功率和无功功率给定值；分别表示节点i的k谐波有功功率和无功功率计算值；h表示计算的最高次谐波。

在基波潮流解的基础上，谐波功率根据式(4)计算；

Pik=RiikKik[(ai1)2+(bi1)2]Qik=XiikKik[(ai1)2+(bi1)2]---(4)]]>

其中，和分别为节点i每次迭代过程中的基波电流的实部和虚部，其计算如式(5)所示；和分别为k次谐波导纳矩阵的逆矩阵的对应元素的实部和虚部；为节点i的k次谐波电流含有率；谐波功率只在有谐波注入电流的节点其值存在，其他节点的计算值为0；

ai1=Σj=1n(GijUjcosθj-BijUjsinθj)bi1=Σj=1n(GijUjsinθj+BijUjcosθj)---(5)]]>

根据式(1)到式(4)进行迭代求解直到前后两次电压差满足收敛条件即可完成当次迭代过程；

max(|ΔU_i|)≤ε    (6)

根据谐波源特性计算得到谐波注入电流，利用谐波节点电压方程(7)以及支路电流方程(8)即可求得谐波节点电压和支路电流；

Unh=(Ynh)-1Inh---(7)]]>

Ibh=YbhATUnh---(8)]]>

其中，表示h次谐波的节点电压向量、节点注入电流向量和支路电流向量；表示h次谐波节点导纳矩阵和支路导纳矩阵；A表示节点支路关联矩阵；

步骤S04：将节点功率和谐波特性的样本集重复代入式(1)到式(4)的谐波潮流迭代模型，每一组的节点功率和谐波特性都能得到一组对应的谐波电压和谐波电流计算结果，用n表示节点数，N表示样本数，N1表示支路数，表示h次谐波节点电压计算结果和支路电流计算结果，则计算结果可表示如下；

U~nh={Ui,jh,i=2,3,...,n,j=1,2,...,N}I~bh={Ii,jh,i=2,3,...,N1,j=1,2,...,N}.]]>

2.根据权利要求1所述的谐波污染用户的概率潮流评估方法，其特征在于：所述样本集获取方式为：对电网按其网络拓扑结构建模，按网络的接线情况形成用以描述谐波网络的谐波导纳矩阵；在形成各次谐波导纳矩阵时，基波分量的导纳元素与网络的基波阻抗相对应，谐波分量的导纳元素与谐波的频率相对应。