基于原对偶内点法的非线性规划潮流计算方法技术

本发明专利技术提供了一种基于原对偶内点法的非线性规划潮流计算方法，包括如下步骤：(1)构造并求解无调整量的非线性规划潮流模型，得到目标函数值；(2)判断目标函数值是否为0，若为0则直接输出潮流结果；否则执行步骤3；(3)对可调整的调整量进行调整；(4)构造并求解有调整量的非线性规划潮流模型，得到目标函数值；(5)判断目标函数值是否为0，若为0，则输出调整量的调整信息和潮流计算结果，程序结束；否则跳转到步骤3。本发明专利技术提供的基于原对偶内点法的非线性规划潮流计算方法，在传统算法不收敛时，判断潮流解是否存在，准确得到潮流有解还是无解的信息；而且在潮流解不存在的情况下，给出潮流调整措施，并能自动进行潮流数据调整，得到正常的潮流解。

【技术实现步骤摘要】
基于原对偶内点法的非线性规划潮流调整方法
本专利技术属于电力系统仿真分析领域，具体涉及一种基于原对偶内点法的非线性规划潮流调整方法。
技术介绍
潮流计算是电力系统最基本的分析计算。牛顿-拉夫逊法(简称牛顿法)是求解非线性方程组的经典方法，特别是稀疏存储技术和节点优化编号技术的引入后，因其收敛性好、计算速度快、内存需求少的优点，广泛应用于潮流计算领域，至今仍然是求解潮流计算问题的主流算法。但牛顿法难以处理潮流计算不收敛情况。潮流计算不收敛的原因从原理上可分为两种情况：一种是潮流解(实数解)不存在，这种情况不收敛是正常的，需要计算人员调整潮流数据本身；一种是潮流解存在，但是由于计算方法本身的问题没有计算出来，这种情况是不正常的，需要更换或改进潮流算法以得到潮流解。牛顿法无法进行不收敛原因分析，如果出现不收敛情况无法判断是潮流解不存在还是算法本身的问题。根据非线性规划的原理提出了最优乘子法，将最优乘子与常规的牛顿法相结合，在牛顿法状态变量的方向上乘以一个最优乘子，通过最优乘子信息可判断潮流解是否存在。如果潮流迭代计算过程中，最优乘子趋近于0，则说明是潮流解不存在；如果最优乘子趋近于1，则说明是算法本身问题。现有经典的潮流计算方法--牛顿法最大的缺点是无法进行不收敛原因分析，如果出现不收敛情况无法判断是否是算法本身的问题。尤其是近年来随着我国电网规模的不断扩大，在实际电网分析中常常出现采用牛顿法及其改进方法都不能很好收敛的情况，这种情况下牛顿法本身无法做出任何处理。最优乘子法理论上能给出潮流解存在与否的信息，但仅仅得到潮流解不存在的信息是不够的，潮流解不存在并不是潮流计算人员希望的最终结果，潮流计算人员往往只能自行调整潮流数据以期望得到正常的潮流解。在潮流解不存在的情况下，最优乘子法具有不能给出进一步的信息来指导用户进行潮流数据调整以得到潮流解的问题。
技术实现思路
为克服上述缺陷，本专利技术提供了一种基于原对偶内点法的非线性规划潮流计算方法，在传统算法不收敛时，判断潮流解是否存在，准确得到潮流有解还是无解的信息；而且在潮流解不存在的情况下，给出潮流调整措施，并能自动进行潮流数据调整，得到正常的潮流解。为实现上述目的，本专利技术提供一种基于原对偶内点法的非线性规划潮流计算方法，其改进之处在于，所述方法包括如下步骤：(1).构造并求解无调整量的非线性规划潮流模型，得到目标函数值；(2).判断目标函数值是否为0，若为0则直接输出潮流结果；否则执行步骤3；(3).对可调整的调整量进行调整；(4).构造并求解有调整量的非线性规划潮流模型，得到目标函数值；(5).判断目标函数值是否为0，若为0，则输出调整量的调整信息和潮流计算结果，程序结束；否则跳转到步骤3。本专利技术提供的优选技术方案中，在所述步骤1中，无调整量的情况下只求解潮流方程；电力系统的潮流方程可表述为：其中，Pi、Qi分别为第i个节点的注入有功功率、注入无功功率，Gij、Bij为第i个节点和第j个节点之间的节点导纳的电导分量和电纳分量，为已知量；Vi、θi分别为第i个节点的电压幅值和电压相角，是待求量；NB为节点总数；将待求量统一用X表示，潮流方程即可统一用g(X)＝0表示；引入潮流方程误差量，构造方程误差量平方和的目标函数，可构造非线性规划潮流模型：其中，ε为每个潮流方程的误差量，当潮流方程有解时，每个潮流方程的ε均为0，即目标函数为0，求解的待求量X为正常的潮流结果；若目标函数不为0，则说明原始潮流方程无解，求解的待求量X为潮流方程误差最小的最小二乘解，不是正常的潮流结果，潮流解不存在。本专利技术提供的第二优选技术方案中，在所述步骤1中，采用原对偶内点法求解潮流模型。本专利技术提供的第三优选技术方案中，原对偶内点法是在非线性规划问题中的不等式约束加入松弛变量，如公式4所示：由式(4)可得拉格朗日函数：其中，λ为拉格朗日乘子，Sl、Su为上、下限松弛变量，πl、πu为对偶变量，μ为障碍因子。根据K-T条件求解式(5)可得：其中，H阵被称为海森阵：求解以上线性方程组，更新原、对偶变量：其中，σ为常数，一般0.99≤σ≤0.9995；αP，αD为原、对偶变量的步长因子，其设置是为了保证所有变量可行性，一般由下式确定：障碍因子μ的选取直接影响算法的收敛性，无论原变量还是对偶变量，其转移方向都受障碍因子影响；障碍因子采用互补间隙法确定：ρ(K+1)＝-Sl(K+1)Tπl(K+1)+Su(K+1)Tπu(K+1)(12)μ(K+1)＝ρ(K+1)/2nβ(13)其中，ρ为互补间隙，β为加速因子，n为同一类型变量的数目。本专利技术提供的第四优选技术方案中，在所述步骤3中，调整量包括：变压器变比、发电机有功无功功率，负荷有功无功功率和无功补偿容量。本专利技术提供的第五优选技术方案中，对变压器的变比的调整包括：取变比的倒数(1/Tk)为调整量；目标函数、节点平衡方程约束与变比调整量的关系为一次函数关系，避免变比为分母的强非线性问题；设(i，k)支路为变压器，变比为1:Tk，令Trec＝1/Tk，节点i注入的有功功率为：则潮流等式方程对Tk、Trec的二阶导数为：可见，采用1/Tk为调整量避免了Tk为调整量的强非线性。本专利技术提供的第六优选技术方案中，发电机有功功率Pg和无功功率Qg作为调整量，根据发电机功率的上下限约束条件对发电机有功无功功率进行调整：Pgmin≤Pg≤Pgmax(17)Qgmin≤Qg≤Qgmax(18)本专利技术提供的第六优选技术方案中，负荷有功功率Pl和无功功率Ql作为调整量，根据负荷功率的上下限约束条件对负荷有功、无功功率进行调整：Plmin≤Pl≤Plmax(19)Qlmin≤Ql≤Qlmax(20)本专利技术提供的第七优选技术方案中，无功补偿包括：并联电容电抗器的投切和未知的负荷侧无功补偿。本专利技术提供的第八优选技术方案中，对并联电容电抗器的投切的调整通过控制支路的导纳值实现；设i支路为可控支路，阻抗为R+jX，导纳值G+jB，则节点i注入的功率为：其中，则潮流有功等式方程对X、B的一阶导数为：取导纳B为调整量避免了X为调整量的强非线性。本专利技术提供的第九优选技术方案中，利用负荷的无功功率Ql，根据补偿设备总的容量限制，调整未知的负荷侧的无功补偿。本专利技术提供的第十优选技术方案中，在所述步骤4中，通过优化各个调整量使潮流调整到有解的状态，其数学模型如下：其中，X为系统待求量，u为调整量；g(X,u)为潮流方程；h(u)为不等式约束方程，即调整量的调节范围约束，hl、hu为调整量的下限和上限。与现有技术比，本专利技术提供的一种基于原对偶内点法的非线性规划潮流计算方法，解决判断潮流解是否存在的问题以及潮流解不存在的潮流调整问题。与传统的牛顿法相比，可以准确得到潮流有解还是无解的信息；与最优乘子法相比，可以在潮流解不存在的情况下，可以给出潮流调整措施，并能自动进行潮流数据调整，得到正常的潮流解。附图说明图1为一般矩阵排列方式。图2为改进后矩阵排列方式。图3为基于原对偶内点法的非线性规划潮流计算方法的流程图。具体实施方式如图3所示，一种基于原对偶内点法的非线性规划潮流计算方法，包括如下步骤：(1).构造并求解无调整量的非线性规划潮流模型，得到目标函数值；(2).判断目标函数值是否为0，若为本文档来自技高网...

【技术保护点】
一种基于原对偶内点法的非线性规划潮流计算方法，其特征在于，所述方法包括如下步骤：?(1).构造并求解无调整量的非线性规划潮流模型，得到目标函数值；?(2).判断目标函数值是否为0，若为0则直接输出潮流结果；否则执行步骤3；?(3).对可调整的调整量进行调整；?(4).构造并求解有调整量的非线性规划潮流模型，得到目标函数值；?(5).判断目标函数值是否为0，若为0，则输出调整量的调整信息和潮流计算结果，程序结束；否则跳转到步骤3。

【技术特征摘要】
1.一种基于原对偶内点法的非线性规划潮流调整方法，其特征在于，所述方法包括如下步骤：(1).构造并求解无调整量的非线性规划潮流模型，得到目标函数值；(2).判断目标函数值是否为0，若为0则直接输出潮流结果；否则执行步骤(3)；(3).对可调整的调整量进行调整；所述调整量包括：变压器变比，对所述变压器变比的调整包括：取变比的倒数1/Tk为调整量；目标函数、节点平衡方程约束与变比调整量的关系为一次函数关系，避免变比为分母的强非线性问题；设(i，k)支路为变压器，令Trec＝1/Tk，节点i注入的功率为：Vi、Vj分别为第i个节点、第j个节点的电压幅值；NB为节点总数；Gij为第i个节点和第j个节点之间的导纳的电导分量，Bij为第i个节点和第j个节点之间的导纳的电纳分量；θij为第i个节点和第j个节点之间的电压相角差，Vk为第k个节点的电压幅值，Gik和Bik为第i个节点和第k个节点之间的导纳的电导分量和电纳分量；则潮流等式方程对Tk、Trec的二阶导数为：其中，θik为第i个节点和第k个节点的电压相角差，可见，采用1/Tk为调整量避免了Tk为调整量的强非线性；(4).构造并求解有调整量的非线性规划潮流模型，得到目标函数值；(5).判断目标函数值是否为0，若为0，则输出调整量的调整信息和潮流计算结果，程序结束；否则跳转到步骤(3)。2.根据权利要求1所述的方法，其特征在于，在所述步骤(1)中，无调整量的情况下求解潮流方程；电力系统的潮流方程可表述为：其中，Pi、Qi分别为第i个节点的注入有功功率、注入无功功率，Gij、Bij为第i个节点和第j个节点之间的导纳的电导分量和电纳分量，为已知量；Vi、θi分别为第i个节点的电压幅值和电压相角，是待求量,Vj为第j个节点的电压幅值；NB为节点总数；θj为第j个节点的电压相角；将待求量统一用X表示，潮流方程即可统一用g(X)＝0表示；引入潮流方程误差量，构造方程误差量平方和的目标函数，可构造非线性规划潮流模型：其中，ε为每个潮流方程的误差量，当潮流方程有解时，每个潮流方程的ε均为0，即目标函数为0，求解的待求量X为正常的潮流结果；若目标函数不为0，则说明原始潮流方程无解，求解的待求量X为潮流方程误差最小的最小二乘解，不是正常的潮流结果，潮流解不存在。3.根据权利要求1所述的方法，其特征在于，采用原对偶内点法求解所述步骤(1)的潮流模型。4.根据权利要求3所述的方法，其特征在于，原对偶内点法是在非线性规划问题中的不等式约束加入松弛变量，如公式(7)所示：由式(7)可得拉格朗日函数：其中，X为系统待求量，hl和hu为不等式约束的下限和上限，λ为拉格朗日乘子，Sl、Su为上、下限松弛变量，πl、πu为对偶变量，μ为障碍因子；根据K-T条件求解式(8)可得：其中，H阵被称为海森阵：

【专利技术属性】
技术研发人员：陈兴雷，李芳，
申请(专利权)人：国家电网公司，中国电力科学研究院，
类型：发明
国别省市：