首页 / 专利库 / 专利权 / 专利合作条约 / 第I章 / 国际申请 / 请求书 / 保护类型 / 专利 / 首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法

首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法

阅读:454发布:2022-10-01

专利汇可以提供首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种首次 迭代 小阻抗支路端点改变雅可比矩阵的潮流计算方法,在首次迭代时小阻抗支路端点采用由给定值Pis和Qis计算出的ai和bi值计算雅可比矩阵元素,首次迭代时正常支路端点以及后续各次迭代时所有 节点 则采用传统方法计算雅可比矩阵元素。本发明通过在首次迭代过程中小阻抗支路两端节点采用与以后各次迭代过程不同的雅可比矩阵计算方法,解决了直 角 坐标 牛 顿法潮流计算在分析含有小阻抗支路电 力 系统时的收敛性问题。采用常规直角坐标牛顿法潮流计算不收敛时,本发明能够可靠收敛,且比现有 专利 技术迭代次数少。本发明同时也能对正常电力系统进行潮流计算,没有不良影响。,下面是首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法专利的具体信息内容。

1.首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法,其特征在于:包括以下步骤:
A、输入原始数据和初始化电压
根据电系统节点的特点,潮流计算把电力系统节点分成3类:节点有功功率和无功功率已知、节点电压幅值和电压相未知的节点称为PQ节点;节点有功功率和电压幅值已知、节点无功功率和电压相角未知的节点称为PV节点;节点电压幅值和电压相角已知,节点有功功率和无功功率未知的节点称为平衡节点;
电压初始化采用平启动,即PV节点和平衡节点的电压实部取给定值,PQ节点的电压实部取1.0;所有电压的虚部都取0.0;这里单位采用标幺值;
B、根据支路电阻和电抗的大小确定两端节点所连支路类型数组T
形成节点所连支路类型数组的具体步骤如下:
B1、读入支路数据,设置小电阻阈值rmin和小电抗阈值xmin;
B2、节点所连支路类型数组T清零;
B3、令q=1;
B4、取支路q的首末节点号i和j、电阻r、电抗x;
B5、判断是否满足r≤rmin且x≤xmin的条件,如果不满足,转步骤B7;
B6、令Ti=1,Tj=1;
B7、令q=q+1;
B8、判断q是否大于支路数l,如果q不大于l转步骤B4;否则转步骤C;
C、形成节点导纳矩阵
设节点i和节点j原来的自电导与自电纳分别为Gi0、Bi0、Gj0、Bj0,在它们之间增加一条小阻抗支路后的自导纳和互导纳分别为:
式中,Yii、Yjj分别为节点i和节点j的自导纳;Yij为节点i和节点j之间的互导纳;rij、xij分别为节点i和节点j之间小阻抗支路的电阻和电抗;k为节点i和节点j之间小阻抗支路的变比,如果是输电线支路,变比k为1;
D、设置迭代计数t=0;
E、计算功率及电压偏差,求最大不平衡量ΔWmax;
PQ节点的功率偏差计算公式为:
式中,Pis、Qis分别为节点i给定的注入有功功率和无功功率;ei和fi分别为节点i的电压相量的实部和虚部;ai、bi分别为节点i的计算注入电流相量的实部和虚部,为式中,n为电力系统的节点数;
潮流计算收敛时,式(4)中ΔPi、ΔQi都趋近于0,因此ai和bi等于由给定值Pis和Qis计算出的ais和bis
PV节点的有功功率及电压偏差计算公式为:
式中,Vis为节点i给定的电压幅值;
平衡节点不参与迭代计算,不需要计算功率偏差或电压偏差;
求各节点功率或电压偏差中绝对值最大的值,称为最大不平衡量ΔWmax;
F、判断最大不平衡量绝对值|ΔWmax|是否小于收敛精度ε;如果小于收敛精度ε,执行步骤J;否则,执行步骤G;
G、形成雅可比矩阵J
除首次迭代外,雅可比矩阵计算方法仍采用传统方法;首次迭代的雅可比矩阵计算方法根据节点连接支路的类型采用不同方法;对于小阻抗支路的端点,因为采用传统方法计算雅可比矩阵会导致潮流计算发散,所以计算雅可比矩阵元素时采用式(6)计算注入电流相量的实部和虚部效果较好;对于正常阻抗支路的端点,仍然按传统方法计算雅可比矩阵元素,即雅可比矩阵计算公式中的节点i注入电流相量的实部ai和虚部bi按式(5)计算;
形成雅可比矩阵元素的具体步骤如下:
G1、计算i≠j时的雅可比矩阵元素;
当i≠j时,雅可比矩阵J的元素计算公式如下:
G2、令i=1;
G3、判断是否同时满足t=0且Ti=1的条件,如果不满足转步骤G4;如果满足,则按式(6)计算节点i的注入电流相量的实部ai和虚部bi,然后转步骤G5;
G4、按式(5)计算节点i的注入电流相量的实部ai和虚部bi;
G5、计算i=j时的雅可比矩阵元素;
当i=j时,雅可比矩阵J的元素计算公式如下:
PQ节点按式(14)-(17)计算雅可比矩阵元素;PV节点按式(14)、(15)、(18)、(19)计算雅可比矩阵元素;平衡节点不计算雅可比矩阵元素;
G6、令i=i+1;
G7、判断i是否大于节点数n,如果i不大于n转步骤G3;否则转步骤H;
H、解修正方程及修正电压实部e、虚部f
潮流计算的基本方程(4)和(7)是非线性方程组,采用逐次线性化方法迭代求解;线性化得到的方程称为修正方程,用来求电压实部和虚部的修正量;
修正方程为:
式中,J为雅可比矩阵;ΔP和ΔQ分别为有功功率和无功功率偏差列向量;ΔV2为电压幅值偏差列向量;Δe和Δf分别为电压相量的实部和虚部修正量列向量; 为有功功率偏差函数列向量对电压相量实部列向量转置的偏导矩阵,上标T为转置符号;
电压修正公式为:
式中,上标t表示第t次迭代;
I、令t=t+1,返回步骤E进行下一次迭代;
J、输出节点及支路数据。

说明书全文

首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法

技术领域

[0001] 本发明涉及含小阻抗支路电系统的潮流计算方法,特别是一种电力系统的直坐标顿法潮流计算方法。

背景技术

[0002] 电力系统潮流计算是研究电力系统稳态运行的一项基本计算,它根据电力系统给定的运行条件和网络结构确定整个电力系统的运行状态。潮流计算也是电力系统其他分析的基础,如安全分析、暂态稳定分析等都要用到潮流计算。由于具有收敛可靠、计算速度较快及内存需求适中的优点,牛顿法成为当前潮流计算的主流方法。牛顿法分为极坐标和直角坐标两种形式,两种形式的牛顿法潮流计算都在电力系统中得到了广泛的应用。
[0003] 在直角坐标牛顿法潮流计算中,节点i的电压采用直角坐标表示为:
[0004] 对正常电力网络,直角坐标牛顿法潮流计算具有良好的收敛性,但遇到含有小阻抗支路的病态网络时,直角坐标牛顿法潮流计算就可能发散。电力系统小阻抗支路可分为小阻抗线路和小阻抗变压器支路,在数学模型上线路可以看作变比为1:1的变压器,因此下面分析时仅以小阻抗变压器支路为例分析。小阻抗变压器模型见图1,变压器的非标准变比k位于节点i侧,阻抗位于标准变比侧。变压器阻抗zij=rij+jxij很小,导纳为[0005]
[0006] 式中,yij、gij、bij分别为节点i和节点j之间小阻抗支路的导纳、电导和电纳;rij、xij分别为节点i和节点j之间小阻抗支路的电阻和电抗。
[0007] 由于小阻抗支路li-j的阻抗很小,支路的电压降也很小,因此变压器两端节点的电压应满足:
[0008]
[0009] 如图2所示,现有直角坐标牛顿法潮流计算方法,主要包括以下步骤:
[0010] A、输入原始数据和初始化电压
[0011] 根据电力系统节点的特点,潮流计算把电力系统节点分成3类:节点有功功率和无功功率已知、节点电压幅值和电压相角未知的节点称为PQ节点;节点有功功率和电压幅值已知、节点无功功率和电压相角未知的节点称为PV节点;节点电压幅值和电压相角已知,节点有功功率和无功功率未知的节点称为平衡节点。
[0012] 电压初始化采用平启动,即PV节点和平衡节点的电压实部取给定值,PQ节点的电压实部取1.0;所有电压的虚部都取0.0。这里单位采用标幺值。
[0013] B、形成节点导纳矩阵
[0014] 设节点i和节点j原来的自电导与自电纳分别为Gi0、Bi0、Gj0、Bj0,在它们之间增加一条小阻抗支路后的自导纳和互导纳分别为:
[0015]
[0016]
[0017]
[0018] 式中,Yii、Yjj分别为节点i和节点j的自导纳;Yij为节点i和节点j之间的互导纳;rij、xij分别为节点i和节点j之间小阻抗支路的电阻和电抗;k为节点i和节点j之间小阻抗支路的变比(如果是输电线支路,变比为1);
[0019] C、计算功率及电压偏差
[0020] PQ节点的功率偏差计算公式为:
[0021]
[0022] 式中,Pis、Qis分别为节点i给定的注入有功功率和无功功率,Pis为电源有功功率与负荷有功功率之差,Qis为电源无功功率与负荷无功功率之差;ai、bi分别为节点i的计算注入电流相量的实部和虚部,为
[0023]
[0024] 式中,n为电力系统的节点数。
[0025] 潮流计算收敛时,式(6)中ΔPi、ΔQi都趋近于0,因此ai和bi等于由给定值Pis和Qis计算出的ais和bis
[0026]
[0027] PV节点的有功功率及电压偏差计算公式为:
[0028]
[0029] 式中,Vis为节点i给定的电压幅值。
[0030] 平衡节点不参与迭代计算,不需要计算功率偏差或电压偏差。
[0031] 求各节点功率或电压偏差中绝对值最大的值,称为最大不平衡量,如果最大不平衡量的绝对值小于给定的收敛精度,转步骤F,否则执行步骤D。
[0032] D、形成雅可比矩阵J
[0033] 雅可比矩阵J的元素(i≠j时)计算公式如下:
[0034]
[0035]
[0036]
[0037]
[0038]
[0039]
[0040] 雅可比矩阵J的元素(i=j时)计算公式如下:
[0041]
[0042]
[0043]
[0044]
[0045]
[0046]
[0047] PQ节点按式(16)-(19)计算雅可比矩阵元素;PV节点按式(16)、(17)、(20)、(21)计算雅可比矩阵元素;平衡节点不计算雅可比矩阵元素。
[0048] E、解修正方程及修正电压实部e、虚部f
[0049] 潮流计算的基本方程(6)和(9)是非线性方程组,通常采用逐次线性化方法迭代求解。线性化得到的方程称为修正方程,用来求电压实部和虚部的修正量。
[0050] 修正方程为:
[0051]
[0052] 式中,J为雅可比矩阵;ΔP和ΔQ分别为有功功率和无功功率偏差列向量;ΔV2为电压幅值偏差列向量;Δe和Δf分别为电压相量的实部和虚部修正量列向量; 为有功功率偏差函数列向量对电压相量实部列向量转置的偏导矩阵,上标T为转置符号。
[0053] 电压修正公式为:
[0054]
[0055] 式中,上标t表示第t次迭代。
[0056] F、输出节点及支路数据。
[0057] 对正常电力网络,牛顿法潮流计算具有良好的收敛性,但遇到含有小阻抗支路的病态网络时,牛顿法潮流计算就可能发散。而电力系统中小阻抗支路普遍存在,收敛性是电力系统潮流计算这类非线性问题的最重要指标,计算不收敛就无法得到方程的解。因此改善直角坐标牛顿法潮流计算针对含有小阻抗支路电力系统的收敛性具有非常重要的意义。
[0058] 中国专利ZL201410299531.5披露了一种通过修改常规直角坐标牛顿法潮流计算雅可比矩阵的方法,改善了潮流计算的收敛性。该方法计算雅可比元素时采用由给定值Pis和Qis计算的ai和bi值,有效解决了含有电阻为0的小阻抗支路电力系统潮流计算的发散问题。但当小阻抗支路的电阻不为0时,该方法迭代次数增加,收敛性变差,甚至不收敛。
[0059] 中国专利ZL201410315785.1提出了一种雅可比矩阵改变的直角坐标牛顿法潮流计算方法,该方法首次迭代和后续各次迭代采用不同的雅可比矩阵计算方法,首次迭代计算雅可比元素时采用由给定值Pis和Qis计算的ai和bi值,后续各次迭代计算雅可比元素时仍采用传统方法,有效解决了含电阻不为0的小阻抗支路电力系统潮流计算的发散问题,但该方法在电力系统包含多条电阻不为0的小阻抗支路时,迭代次数增加,收敛性变差。
[0060] 中国专利ZL201611094297.8提出了一种随迭代和节点类型改变雅可比矩阵的牛顿法潮流计算方法,该方法首次迭代时所有PQ节点与后续各次迭代采用不同的雅可比矩阵计算方法,首次迭代时所有PQ节点计算雅可比元素时采用由给定值Pis和Qis计算的ai和bi值,首次迭代时所有PV节点及后续各次迭代时所有节点计算雅可比元素时仍采用传统方法,有效解决了电力系统包含多条电阻不为0的小阻抗支路的潮流计算的发散问题,但迭代次数仍然较多,需要进一步改进。

发明内容

[0061] 为解决现有技术存在的上述问题,本发明要提出一种首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法,该方法可以提高其分析含有电阻不为0的小阻抗支路电力系统的收敛速度。
[0062] 传统的计算雅可比矩阵的方法是从牛顿法的基本原理推导出来的,正常阻抗支路的端点采用传统方法计算雅可比矩阵元素是合适的,但小阻抗支路的端点采用传统方法计算雅可比矩阵元素则会导致潮流计算发散。首次迭代时,电压为平启动方法设置的电压初值,正常阻抗支路计算的支路功率与实际值相差不大,和这些支路相连的节点的注入功率计算值与给定值接近,因此首次迭代正常阻抗支路的端点使用注入功率计算值来计算雅可比矩阵元素也比较合理。小阻抗支路由于其阻抗很小,两端电压初值与实际值的不一致带来的很小的电压差就会计算出很大的支路功率,与该支路相连的节点的注入功率计算值也很大,会导致潮流计算发散,因此首次迭代小阻抗支路的端点不应该使用注入功率计算值计算雅可比矩阵元素,而应该使用注入功率给定值或初值计算雅可比矩阵元素。
[0063] 为了实现上述目的,本发明提出了一种直角坐标牛顿法潮流计算方法来改善潮流计算收敛性。本发明的首次迭代时小阻抗支路端点采用由给定值Pis和Qis计算出的ai和bi值计算雅可比矩阵元素,首次迭代时正常支路端点以及后续各次迭代时所有节点则采用传统方法计算雅可比矩阵元素。
[0064] 本发明的技术方案如下:首次迭代小阻抗支路端点改变雅可比矩阵的潮流计算方法,包括以下步骤:
[0065] A、输入原始数据和初始化电压;
[0066] B、根据支路电阻和电抗的大小确定两端节点所连支路类型T
[0067] 形成节点所连支路类型数组的具体步骤如下:
[0068] B1、读入支路数据,设置小电阻阈值rmin和小电抗阈值xmin;
[0069] B2、节点所连支路类型数组T清零;
[0070] B3、令m=1;
[0071] B4、取支路m的首末节点号i和j、电阻r、电抗x;
[0072] B5、判断是否满足r≤rmin且x≤xmin的条件,如果不满足,转步骤B7;
[0073] B6、令Ti=1,Tj=1;
[0074] B7、令m=m+1;
[0075] B8、判断m是否大于支路数l,如果m不大于l转步骤B4;否则转步骤C;
[0076] C、形成节点导纳矩阵;
[0077] D、设置迭代计数t=0;
[0078] E、计算功率及电压偏差,求最大不平衡量ΔWmax;
[0079] F、判断最大不平衡量绝对值|ΔWmax|是否小于收敛精度ε;如果小于收敛精度ε,执行步骤J;否则,执行步骤G;
[0080] G、形成雅可比矩阵;
[0081] 除首次迭代外,雅可比矩阵计算方法仍采用传统方法。首次迭代的雅可比矩阵计算方法根据节点连接支路的类型采用不同方法。对于小阻抗支路的端点,因为采用传统方法计算雅可比矩阵会导致潮流计算发散,所以计算雅可比矩阵元素时采用式(8)计算注入电流相量的实部和虚部效果较好;对于正常阻抗支路的端点,仍然按传统方法计算雅可比矩阵元素,即雅可比矩阵计算公式中的节点i注入电流相量的实部ai和虚部bi按式(7)计算。
[0082] 形成雅可比矩阵元素的具体步骤如下:
[0083] G1、按式(10)-(15)计算i≠j时的雅可比矩阵元素;
[0084] G2、令i=1;
[0085] G3、判断是否同时满足t=0且Ti=1的条件,如果不满足此条件转步骤G4;如果满足,则按式(8)计算节点i的注入电流相量的实部ai和虚部bi,然后转步骤G5;
[0086] G4、按式(7)计算节点i的注入电流相量的实部ai和虚部bi;
[0087] G5、按式(16)-(21)计算i=j时的雅可比矩阵元素;
[0088] G6、令i=i+1;
[0089] G7、判断i是否大于节点数n,如果i不大于n转步骤G3;否则转步骤H;
[0090] H、解修正方程及修正电压实部e、虚部f;
[0091] I、令t=t+1,返回步骤E进行下一次迭代;
[0092] J、输出节点及支路数据。
[0093] 与现有技术相比,本发明具有以下有益效果:
[0094] 1、本发明通过在首次迭代过程中小阻抗支路两端节点采用与以后各次迭代过程不同的雅可比矩阵计算方法,解决了直角坐标牛顿法潮流计算在分析含有小阻抗支路电力系统时的收敛性问题。采用常规直角坐标牛顿法潮流计算不收敛时,本方法能够可靠收敛,且比现有专利技术迭代次数少。
[0095] 2、由于本发明不仅能有效解决常规直角坐标牛顿法潮流计算分析含有小阻抗支路电力系统的收敛性问题,同时也能对正常电力系统进行潮流计算,没有不良影响。附图说明
[0096] 本发明共有附图6张。其中:
[0097] 图1是电力系统小阻抗变压器模型示意图。
[0098] 图2是直角坐标牛顿法潮流计算的流程图
[0099] 图3是专利方法1直角坐标牛顿法潮流计算的流程图。
[0100] 图4是专利方法2直角坐标牛顿法潮流计算的流程图。
[0101] 图5是本发明直角坐标牛顿法潮流计算的流程图。
[0102] 图6是本发明形成节点所连支路类型数组的流程图。

具体实施方式

[0103] 下面结合附图对本发明进行进一步地说明。根据图1所示的小阻抗变压器模型,采用图5-6所示的直角坐标牛顿法潮流计算的流程图,对一个实际大型电网进行了潮流计算。该实际大型电网有445个节点,含有大量的小阻抗支路。其中,x≤0.001的小阻抗支路有49条,x≤0.0001的小阻抗支路有41条,x≤0.00001的小阻抗支路有22条。其中阻抗值最小的是节点118和节点125之间的小阻抗支路l118-125为x=0.00000001,变比k=0.9565,k位于节点118侧。潮流计算的收敛精度为0.00001。为了验证本发明计算含电阻不为0的小阻抗支路电力系统的收敛性,把小阻抗支路l118-125、l60-122及l287-310的电阻改为r=0.0001。
[0104] 作为对比,同时采用以下3种对比方法对该实际大型电网进行了潮流计算:
[0105] 常规方法:常规的直角坐标牛顿法潮流方法;
[0106] 专利方法1:专利号为ZL201410315785.1的专利方法;
[0107] 专利方法2:申请号为ZL201611094297.8的专利方法。
[0108] 迭代次数结果见表1。
[0109] 表1不同潮流方法的迭代结果
[0110]方法 常规方法 专利方法1 专利方法2 本发明方法
迭代结果 不收敛 7次收敛 6次收敛 5次收敛
[0111] 由表1可见,对于修改后的445节点实际电力系统算例,常规直角坐标牛顿法潮流方法不收敛,本发明方法和现有专利方法都能够收敛,但本发明方法的迭代次数比现有专利方法1少2次,比现有专利方法2少1次。
[0112] 不同潮流计算方法各次迭代最大不平衡量见表2。单位为标幺值。
[0113] 表2不同潮流方法各次迭代最大不平衡量
[0114]迭代序号 常规方法 专利方法1 专利方法2 本发明方法
0 -4754.570367135 -4754.570367135 -4754.570367135 -4754.570367135
1 -3451593.823720038 -11.138394991 -3.264368583 23.913925681
2 -886651.468310079 -6.163450054 -0.715148045 3.050019341
3 -222023.112200678 -1.441071252 -0.076847277 0.102604201
4 -55754.415245002 -0.106199006 -0.002294590 -0.000454516
5 -13972.568194423 -0.006353455 -0.000017499 -0.000000012
6 -6386.835620506 -0.000141863 -0.000000001  
7 -6585.38761914 -0.000000062    
8 -378994.776907351      
9 -98508.025841226      
10 -37917.863557986      
[0115] 由表2可知,4种方法首次迭代前最大不平衡量相同且很大。首次迭代后,现有专利方法和本专利方法最大不平衡量明显减少,现有专利方法1迭代7次收敛;现有专利方法2迭代6次收敛;本专利方法最大不平衡量减少速度更快,迭代5次收敛;而常规方法的最大不平衡量则变大,最终发散。
[0116] 小阻抗支路的端点为PV节点的各节点电源无功功率输入值与潮流计算收敛时的计算值及初始计算值见表3。单位为标幺值。
[0117] 表3PV节点的电源无功功率输入值与收敛时计算值及初始计算值
[0118]节点 输入值 收敛计算值 初始计算值
22 1.80000 1.36829 0.15000
400 0.10000 0.69586 204.87496
439 0.80000 0.48861 525.04000
440 0.80000 0.48861 525.04000
[0119] 由表3可见,对于小阻抗支路端点是PV节点的各节点,电源无功功率输入值与潮流计算收敛时的计算值有较大的差值,但初始计算值与潮流计算收敛时的计算值的差值更大。因此首次迭代时,小阻抗支路的端点用输入值计算注入电流相量的实部和虚部要更合理。
[0120] 本发明可以采用任何一种编程语言和编程环境实现,如C语言、C++、FORTRAN、Delphi等。开发环境可以采用Visual C++、Borland C++Builder、Visual FORTRAN等。
[0121] 本发明不局限于本实施例,任何在本发明披露的技术范围内的等同构思或者改变,均列为本发明的保护范围。
相关专利内容
标题 发布/更新时间 阅读量
专利检索系统 2020-05-11 634
一种专利申请平台 2020-05-11 704
一种专利评估系统 2020-05-12 668
专利信息自动采集方法 2020-05-12 495
专利数据显示方法 2020-05-12 626
专利检索方法 2020-05-11 719
专利技术分布分析方法 2020-05-13 183
专利交易系统 2020-05-11 607
核心专利挖掘方法 2020-05-11 827
一种专利申请审批系统 2020-05-11 954
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈