一种交流潮流线性灵敏度分析方法技术

本发明专利技术涉及电力分析领域，更具体地，涉及一种交流潮流线性灵敏度分析方法。本发明专利技术基于线性化潮流计算和节点功率‑节点电压相关性矩阵，通过对节点有功功率、节点无功功率求导，得到节点功率‑潮流灵敏度因子和节点功率‑电压灵敏度因子，把灵敏度因子转化为矩阵形式，得到线性灵敏度矩阵。本发明专利技术的线性灵敏度分析方法在交流潮流方程的基础上做了线性化简化，该简化保留了所有的原始潮流变量，在最大限度上保证了线性化方程的准确性，能够准确快速求解电力系统中的潮流。本发明专利技术的线性灵敏度分析方法含有多个灵敏度因子，可用于基于线性规划的电力系统问题，实现更加全面准确的电力系统优化决策。

A linear sensitivity analysis method of AC power flow

The invention relates to the field of electric power analysis, in particular to an AC power flow linear sensitivity analysis method. Based on the linearized power flow calculation and the node power \u2011 node voltage correlation matrix, the invention obtains the node power \u2011 current sensitivity factor and the node power \u2011 voltage sensitivity factor by deriving the node active power and the node reactive power, converts the sensitivity factor into the matrix form, and obtains the linear sensitivity matrix. The linear sensitivity analysis method of the invention makes the linearization simplification on the basis of the AC power flow equation, which retains all the original power flow variables, ensures the accuracy of the linearization equation to the maximum extent, and can accurately and quickly solve the power flow in the power system. The linear sensitivity analysis method of the invention contains a plurality of sensitivity factors, which can be used for power system problems based on linear programming to realize more comprehensive and accurate power system optimization decisions.

【技术实现步骤摘要】
一种交流潮流线性灵敏度分析方法
本专利技术涉及电力分析领域，更具体地，涉及一种交流潮流线性灵敏度分析方法。
技术介绍
随着电力系统的蓬勃发展和逐步完善，国内外研究学者对电力系统运行进行了详尽的分析。在电力系统灵敏度分析领域，已有不少研究理论。在现有的研究成果中，功率传输分配因子包括SF(shiftfactor转移因子)，DF(deliveryfactor，传输因子)，LF(lossfactor，损耗因子)等。在转移因子的分析中，过去的研究只关注节点注入的有功功率对有功潮流的影响，没有分析节点注入的有功功率对无功潮流的影响、节点注入的无功功率对有功潮流、无功潮流的影响，也没有分析节点注入有功功率和无功功率对网络节点的电压幅值和相位的影响。缺乏这些分析，难以对电力系统的运行状态有全面的把握。
技术实现思路
为了解决现有技术中对电力系统灵敏度分析只关注节点注入的有功功率对有功潮流的影响的不足，本专利技术提供了一种交流潮流线性灵敏度分析方法。为解决上述技术问题，本专利技术的技术方案如下：一种交流潮流线性灵敏度分析方法，包括以下步骤：步骤S1：对电力系统进行潮流线性化处理，并对潮流线性化公式进行化简；步骤S2：根据潮流线性化处理的结果构建节点功率-节点电压相关性矩阵；步骤S3：基于节点功率-节点电压相关性矩阵求解节点功率-潮流灵敏度因子以及节点功率-电压灵敏度因子；步骤S4：将节点功率-潮流灵敏度因子以及节点功率-电压灵敏度因子转化为矩阵，得到线性灵敏度矩阵的集合；步骤S5：基于线性灵敏度的集合分析电力系统的运行状态。优选的，在步骤S3中，对节点功率-节点电压相关性矩阵中的节点有功功率以及节点无功功率求导得到求解节点功率-潮流灵敏度因子以及节点功率-电压灵敏度因子。优选的，所述的节点功率-节点电压相关性矩阵表征了节点有功功率、节点无功功率与节点电压幅值、节点电压相角的线性关联。优选的，步骤S1的具体步骤如下：原始节点电压方程为：其中，Pi，Qi分别表示母线i注入的有功功率和无功功率；Yij＝Gij+jBij是节点导纳矩阵Y中第i行第j列的元素；Y′ij＝G′ij+jB′ij是忽略了并联元件的节点导纳矩阵Y′ij中第i行第j列的元素；gij+jbij是线路(i,j).的导纳；假设Gij＝G′ij，与电纳相比，并联电导忽略不计；Vi,Vj分别为节点i，j的电压幅值，θi,θj分别为节点i，j的电压相位，n是节点的索引，Gij，Bij分别是节点电导矩阵、节点电纳矩阵中第i行第j列的元素；在电力系统中，母线电压幅值标幺值为1，电压的相角绝对值不超过30°，因此作以下假设，对潮流进行线性化处理：θi≈0°,θj≈0°(3)sin(θi-θj)≈(θi-θj)(4)cos(θi-θj)≈1(5)假设Vi,Vj的幅值接近标幺值1，即Vi＝1+ΔVi，ΔVi≈0，Vj＝1+ΔVj，ΔVj≈0对式(1)进行化简：其中，B′ij是忽略了并联元素的节点电纳矩阵中第i行第j列的元素。同理，对式(2)进行化简，得到：其中，G′ij是忽略了并联元素的节点电纳矩阵中第i行第j列的元素。线性化后的潮流方程如下所示：pfij＝gij(Vi-Vj)-bij(θi-θj)(11)qfij＝-bij(Vi-Vj)-gij(θi-θj)(12)pfij，qfij分别是线路(i,j)的有功潮流和无功潮流。优选的，步骤S2构建相关性矩阵的具体步骤如下：将式(9)与(10)转换为矩阵形式，对母线i进行分析：其中，Gi,B′i和Bi分别是矩阵G＝[Gij]n×n，B'＝[B′ij]n×n和B＝[Bij]n×n的第i行元素；V＝[Vi]n×1和θ＝[θi]n×1分别是母线电压的幅值和相角，[*]i,:表示矩阵*的第i行的所有元素；将式(14)转化为适用于所有母线的形式，同时删除参考节点所在的行与列，得到：对式(15)的常数矩阵求逆，得到：其中，K(n-1)×(n-1),L(n-1)×(n-1),M(n-1)×(n-1)和N(n-1)×(n-1)是逆阵的分块矩阵。优选的，在步骤S4中，线性灵敏度矩阵的集合包括节点注入有功功率对节点电压幅值的转移因子SFVp，节点注入无功功率对节点电压幅值的转移因子SFVq，节点注入有功功率对节点电压相角的转移因子SFθp，节点注入无功功率对节点电压相角的转移因子SFθq，节点注入有功功率对线路有功潮流的转移因子SFPp，节点注入有功功率对线路无功潮流的转移因子SFPq，节点注入无功功率对线路有功潮流的转移因子SFQp，节点注入无功功率对线路无功潮流的转移因子SPQq；其中，K、L、M、N是常数矩阵，P、Q分别是节点注入有功功率、节点注入无功功率，Pinj，Qinj分别是节点注入有功功率、节点注入无功功率；SFVp和SFVq是(n-1)×(n-1)的灵敏度矩阵，表征节点注入的单位有功功率和无功功率引起的电压幅值变化，SFθp和SFθq是(n-1)×(n-1)的灵敏度矩阵，表征节点注入的单位有功功率和无功功率引起的电压相角变化；其中，K、L、M、N是常数矩阵，P、Q分别是节点注入有功功率、节点注入无功功率，Pinj，Qinj分别是节点注入有功功率、节点注入无功功率；SFPp和SFPq是(n-1)×l的灵敏度矩阵，表征节点注入的单位有功功率和无功功率引起的支路有功潮流变化，[*]l表示矩阵*的l行，i和j分别是线路l的首端和末端；其中，SFQp和SFQq是(n-1)×l的灵敏度矩阵，表征节点注入的单位有功功率和无功功率引起的支路无功潮流变化。与现有技术相比，本专利技术技术方案的有益效果是：本专利技术的线性灵敏度分析方法在交流潮流方程的基础上做了线性化简化，该简化保留了所有的原始潮流变量，在最大限度上保证了线性化方程的准确性，能够准确快速求解电力系统中的潮流。本专利技术所提出的灵敏度矩阵能够基于网络拓扑和运行信息，全面分析电力系统的运行状态。本专利技术的线性灵敏度分析方法含有多个灵敏度因子，可用于基于线性规划的电力系统问题，实现更加全面准确的电力系统优化决策。附图说明图1为本专利技术的方法流程图。具体实施方式附图仅用于示例性说明，不能理解为对本专利的限制；为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。下面结合附图和实施例对本专利技术的技术方案做进一步的说明。实施例1如图1所示，一种交流潮流线性灵敏度分本文档来自技高网...

【技术保护点】
1.一种交流潮流线性灵敏度分析方法，其特征在于，包括以下步骤：/n步骤S1：对电力系统进行潮流线性化处理，并对潮流线性化公式进行化简；/n步骤S2：根据潮流线性化处理的结果构建节点功率-节点电压相关性矩阵；/n步骤S3：基于节点功率-节点电压相关性矩阵求解节点功率-潮流灵敏度因子以及节点功率-电压灵敏度因子；/n步骤S4：将节点功率-潮流灵敏度因子以及节点功率-电压灵敏度因子转化为矩阵，得到线性灵敏度矩阵的集合；/n步骤S5：基于线性灵敏度的集合分析电力系统的运行状态。/n

【技术特征摘要】
1.一种交流潮流线性灵敏度分析方法，其特征在于，包括以下步骤：
步骤S1：对电力系统进行潮流线性化处理，并对潮流线性化公式进行化简；
步骤S2：根据潮流线性化处理的结果构建节点功率-节点电压相关性矩阵；
步骤S3：基于节点功率-节点电压相关性矩阵求解节点功率-潮流灵敏度因子以及节点功率-电压灵敏度因子；
步骤S4：将节点功率-潮流灵敏度因子以及节点功率-电压灵敏度因子转化为矩阵，得到线性灵敏度矩阵的集合；
步骤S5：基于线性灵敏度的集合分析电力系统的运行状态。


2.根据权利要求1所述的一种交流潮流线性灵敏度分析方法，其特征在于，在步骤S3中，对节点功率-节点电压相关性矩阵中的节点有功功率以及节点无功功率求导得到求解节点功率-潮流灵敏度因子以及节点功率-电压灵敏度因子。


3.根据权利要求2所述的一种交流潮流线性灵敏度分析方法，其特征在于，所述的节点功率-节点电压相关性矩阵表征了节点有功功率、节点无功功率与节点电压幅值、节点电压相角的线性关联。


4.根据权利要求3所述的一种交流潮流线性灵敏度分析方法，其特征在于，步骤S1的具体步骤如下：
原始节点电压方程为：






其中，Pi，Qi分别表示母线i注入的有功功率和无功功率；Yij＝Gij+jBij是节点导纳矩阵Y中第i行第j列的元素；Y′ij＝G′ij+jB′ij是忽略了并联元件的节点导纳矩阵Y′ij中第i行第j列的元素；gij+jbij是线路(i,j).的导纳；假设Gij＝G′ij，与电纳相比，并联电导忽略不计；Vi,Vj分别为节点i，j的电压幅值，θi,θj分别为节点i，j的电压相位，n是节点的索引，Gij，Bij分别是节点电导矩阵、节点电纳矩阵中第i行第j列的元素；
在电力系统中，母线电压幅值标幺值为1，电压的相角绝对值不超过30°，因此作以下假设，对潮流进行线性化处理：
θi≈0°,θj≈0°(3)
sin(θi-θj)≈(θi-θj)(4)
cos(θi-θj)≈1(5)
假设Vi,Vj的幅值接近标幺值1，即Vi＝1+ΔVi，ΔVi≈0，Vj＝1+ΔVj，ΔVj≈0



对式(1)进行化简：



其中，B′ij是忽略了并联元素的节点电纳矩阵中第i行第j列的元素。
同理，对式(2)进行化简，得到：






其中，G′ij是忽略了并联元素的节点电纳矩阵中第i行第j列的元素。
线性化后的潮流方程如下所示：






pfij＝gij(Vi-Vj)-bij(θi-θj)(11)
qfij...

【专利技术属性】
技术研发人员：李镇豪，王逸飞，赖来利，
申请(专利权)人：广东工业大学，
类型：发明
国别省市：广东;44