一种电力系统静态电压稳定临界点的计算方法
专利汇可以提供一种电力系统静态电压稳定临界点的计算方法专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种电 力 系统静态 电压 稳定 临界点 的计算方法，解决了雅克比矩阵奇异导致静态电压稳定临界点处潮流不可解的问题。主要步骤包括：计算 基础 方式潮流、增加 节点 负荷至经典 牛 拉法不收敛、采用AQ节点法计算系统静态电压稳定临界点等。本发明通过在静态电压 稳定性 分析中引入一种新的节点类型——AQ节点，改善了传统连续潮流分析方法需要通过反复预测校正以求取静态电压稳定临界点，计算过程复杂、耗时长且在校正过程中仍可能出现不收敛现象的 缺陷 ，该方法可引起矩阵奇异点在静态电压稳定临界点处的“偏移”，完全消除雅克比矩阵在静态电压稳定临界点附近的奇异性，直接对系统静态电压稳定临界点进行快速、准确求取。，下面是一种电力系统静态电压稳定临界点的计算方法专利的具体信息内容。
1.一种电力系统静态电压稳定临界点的计算方法，包括以下步骤：
第一步、给定电力系统的基础运行方式，计算初始潮流P0、初始电压U0、功率因数cosφi，确定负荷和发电机增加的步长λLi、λG；
第二步、逐步增加节点负荷，采用经典牛顿-拉夫逊潮流方法计算至结果不收敛，选择负荷增量最大的节点为系统中的AQ节点，记录最后一次收敛结果中AQ节点和平衡节点的相角差θ；
第三步、采用AQ节点法，增加AQ节点和平衡节点的相角差θ至系统达到最大负荷点P0max；
第四步、列举系统中可能出现的故障运行方式集合C(C1,C2，…，Cn)，重复步骤3中的方法，求取各故障运行方式下的系统静态电压稳定临界点，若故障方式下采用AQ节点法得不到系统的最大负荷点Pimax或潮流无法收敛，则P0的负荷水平不是故障后合理的运行状态，需减小各节点的负荷计算新的合理初始潮流 重新计算最大负荷；
第五步、根据步骤4中计算得到的不同运行方式下的系统静态电压稳定临界点，确定系统运行的静态电压稳定裕度Pm＝min{Pim},i＝1,...,n。
2.根据权利要求1所述的一种电力系统静态电压稳定临界点的计算方法，其特征在于，所述第二步中，定义了电力系统中一种新的节点类型：AQ节点。AQ节点定义如下：AQ节点是电力系统中一类特殊的负荷节点，其节点相位角θ和无功功率Q均为已知量，电压幅值U和有功功率P为未知量，可以通过求解系统潮流得到该节点的有功功率。
3.根据权利要求1所述的一种电力系统静态电压稳定临界点的计算方法，其特征在于，所述第三步中，提出了一种含AQ节点电力系统的潮流计算方法和其中所包含的一种负荷和发电机的增长策略模型，该计算方法的步骤以及模型如下：
步骤3.1、考虑一个含有NG个发电机节点和NL个负荷节点电力系统，节点总数N＝NG+NL。
假设节点1为平衡节点，节点2，3，…，NG为除平衡节点外的PV类型发电机节点，节点NG+1，NG+
2，…，N为PQ类型负荷节点。得到待求解的有功和无功潮流方程：
ΔPi＝Pi-fPi(θ,U)＝0,i＝2,...,N
ΔQi＝Qi-fQi(θ,U)＝0,i＝NG+1,...,N
其中，Pi和Qi是给定的节点有功和无功功率，fPi(U,θ)＝0和fQi(U,θ)＝0是计算得到的节点有功和无功注入功率，ΔPi和ΔQi是功率不平衡量，U和θ分别为节点电压的幅值和相角；
步骤3.2、根据步骤3.1中列出的系统有功和无功潮流方程，得到系统牛顿-拉夫逊计算方法的潮流迭代方程：
其中，J是维度为(2N-NG+1)的系统雅克比矩阵；
ΔU和Δθ为U和θ的修正量，θ＝[θ2 … θN]T，
步骤3.3、假设选择节点N为AQ节点，其无功功率为QN，相位角为 将牛顿-拉夫逊法的潮流迭代方程修改为：
其中，
步骤3.4、建立符合电力系统实际情况的负荷和发电机的增长策略模型：
模型一、电力系统负荷及发电机增长模型
本发明中考虑多节负荷同时增长，同时由多台发电机增长出力以平衡负荷。假设NP，…，N为负荷增长的节点；1，…，Nq为增加出力的发电机节点。以AQ节点衡量其它节点的负荷增量，以平衡节点衡量其它发电机节点的出力增长，功率因数为cosφ：
Pm-Pm0＝αm(PN-PN0),m＝NP,...,N-1
Pn-Pn0＝βn(P1-P10),n＝2,...,q
其中，αm和βn分别为负荷和发电机的相对变化因子，将有功和无功潮流方程中的函数项修正为：
fPm(θ,U)＝αmfPN(θ,U),m＝Np,...,N-1
fQm(θ,U)＝αmfPN(θ,U)tanφm,m＝Np,...,N-1
fPn(θ,U)＝βnfP1(θ,U),n＝2,...,q
受发电机无功出力的限制，发电机无功功率达到最大值时，将由PV节点转换为PQ节点；
步骤3.5、求取fP1和fPN对θ*和U的偏导数，形成新的含AQ节点的简化雅克比矩阵：
其中，Jk是雅克比矩阵的第k列；同时，将雅克比矩阵J*部分列修改为： i＝Np-1,...,N-2；m＝Np,...,N-1 m＝
Np,...,N-1 i＝1,...,q-1；n＝2,...,q
其中， 为矩阵J*的维度，J*矩阵的其余列保持不变，得到含AQ节点
的电力系统潮流计算方法最终的牛拉法迭代方程：
在上述的一种电力系统静态电压稳定临界点的计算方法，所述步骤3中，提出了另一种衍生的含AQ节点电力系统的潮流计算方法，该方法的思路为：计算过程中直接用J*计算静态电压稳定的临界点，避免对潮流计算中的雅克比矩阵进行修改。计算的基本思路为：每次增加AQ节点和平衡节点的相角差θ时，保持其它负荷节点以及除平衡机外的发电机节点功率不变，计算得到潮流解P*，按计划的模型增加其它负荷和发电机节点，并以P*为初始功率，重新求解潮流，逐步增大θ按此方法迭代，直至AQ节点法不收敛。
一种电力系统静态电压稳定临界点的计算方法
技术领域
背景技术
发明内容
第二步、逐步增加节点负荷，采用经典牛顿-拉夫逊潮流方法计算至结果不收敛，选择负荷增量最大的节点为系统中的AQ节点，记录最后一次收敛结果中AQ节点和平衡节点的相角差θ；
第三步、采用AQ节点法，增加AQ节点和平衡节点的相角差θ至系统达到最大负荷点P0max；
第四步、列举系统中可能出现的故障运行方式集合C(C1,C2，…，Cn)，重复步骤3中的方法，求取各故障运行方式下的系统静态电压稳定临界点，若故障方式下采用AQ节点法得不到系统的最大负荷点Pimax或潮流无法收敛，则P0的负荷水平不是故障后合理的运行状态，需减小各节点的负荷计算新的合理初始潮流P0/，重新计算最大负荷；第五步、根据步骤4中计算得到的不同运行方式下的系统静态电压稳定临界点，确定系统运行的静态电压稳定裕度Pm＝min{Pim},i＝1,...,n。
假设节点1为平衡节点，节点2，3，…，NG为除平衡节点外的PV类型发电机节点，节点NG+1，NG+
2，…，N为PQ类型负荷节点。得到待求解的有功和无功潮流方程：
ΔPi＝Pi-fPi(θ,U)＝0,i＝2,...,N
ΔQi＝Qi-fQi(θ,U)＝0,i＝NG+1,...,N
其中，Pi和Qi是给定的节点有功和无功功率，fPi(U,θ)＝0和fQi(U,θ)＝0是计算得到的节点有功和无功注入功率，ΔPi和ΔQi是功率不平衡量，U和θ分别为节点电压的幅值和相角；
步骤3.2、根据步骤3.1中列出的系统有功和无功潮流方程，得到系统牛顿-拉夫逊计算方法的潮流迭代方程：
其中，J是维度为(2N-NG+1)的系统雅克比矩阵；
ΔU和Δθ为U和θ的修正量，θ＝[θ2···θN]T，
步骤3.3、假设选择节点N为AQ节点，其无功功率为QN，相位角为 将牛顿-拉夫逊法的潮流迭代方程修改为：
其中，
步骤3.4、建立符合电力系统实际情况的负荷和发电机的增长策略模型：
模型一、电力系统负荷及发电机增长模型
本发明中考虑多节负荷同时增长，同时由多台发电机增长出力以平衡负荷。假设NP，…，N为负荷增长的节点；1，…，Nq为增加出力的发电机节点。以AQ节点衡量其它节点的负荷增量，以平衡节点衡量其它发电机节点的出力增长，功率因数为cosφ：
Pm-Pm0＝αm(PN-PN0),m＝NP,...,N-1
Pn-Pn0＝βn(P1-P10),n＝2,...,q
其中，αm和βn分别为负荷和发电机的相对变化因子，将有功和无功潮流方程中的函数项修正为：
fPm(θ,U)＝αmfPN(θ,U),m＝Np,...,N-1
fQm(θ,U)＝αmfPN(θ,U)tanφm,m＝Np,...,N-1
fPn(θ,U)＝βnfP1(θ,U),n＝2,...,q
受发电机无功出力的限制，发电机无功功率达到最大值时，将由PV节点转换为PQ节点；
步骤3.5、求取fP1和fPN对θ*和U的偏导数，形成新的含AQ节点的简化雅克比矩阵：
其中，Jk是雅克比矩阵的第k列；同时，将雅克比矩阵J*部分列修改为：
其中， 为矩阵J*的维度，J*矩阵的其余列保持不变，得到含AQ节点
的电力系统潮流计算方法最终的牛拉法迭代方程：
在上述的一种电力系统静态电压稳定临界点的计算方法，所述步骤3中，提出了另一种衍生的含AQ节点电力系统的潮流计算方法，该方法的思路为：计算过程中直接用J*计算静态电压稳定的临界点，避免对潮流计算中的雅克比矩阵进行修改。计算的基本思路为：每次增加AQ节点和平衡节点的相角差θ时，保持其它负荷节点以及除平衡机外的发电机节点功率不变，计算得到潮流解P*，按计划的模型增加其它负荷和发电机节点，并以P*为初始功率，重新求解潮流，逐步增大θ按此方法迭代，直至AQ节点法不收敛。
图3为AQ节点法计算流程图；
图4为系统机组无功功率和相角差的Q-θ曲线图；
图5为正常运行方式及不同故障运行方式下的P-V曲线；
图6为矩阵J和JAQ在临界点处奇异值的分布；
图7为根据静态电压稳定临界点计算得到的负荷裕度；
图8为本发明计算性能和传统计算方法的效果比较。
具体实施方式
一种电力系统静态电压稳定临界点的计算方法，其分析步骤，包括
1.本实施例以太原220kV实际等值电网算例进行分析：
仿真软件为matlab，系统接线图如图1所示，节点2为3×200MW机组，节点3为4×300MW机组，节点15为600MW机组，系统平均总有功、无功负荷分别为1383MW和453MVar；
在系统的基础运行方式下，计算初始潮流P0、初始电压U0、功率因数cosφi等初始量；
2.逐步增加系统负荷，采用经典牛拉法计算至系统潮流不收敛在系统平均负荷水平下，增加节点1、节点4、节点11和节点14的负荷，同时增加发电机组2、3、15的出力以平衡负荷的增长，其它节点负荷水平不变。选择节点3为系统的平衡节点，节点14为系统的AQ节点，负荷节点和发电机节点的线性增长系数αm、βn如图2所示，节点2无功补偿容量上限为
240MVar，节点15无功补偿容量上限为200MVar，平衡机无功出力不受限。节点1、4、11、14的功率因数cosφi分别为0.81(滞后)、0.99(滞后)、0.97(滞后)、0.97(滞后)；
基础方式下，节点14功率达到6.156pu时，经典牛拉法出现不收敛的现象，此时AQ节点相角差θ14-3＝89.2°；
3.跳转至AQ节点法计算系统静态电压稳定极限点
在步骤2经典牛拉法潮流不收敛的基础上，转至AQ节点法逐步增加节点3和节点14之间的相角差，计算系统的潮流解，详细计算流程如图3所示。当节点14功率达到6.714pu时，θ14-3＝94.1°，达到静态电压稳定临界点，此时临界失稳电压为0.718pu；进一步增加相角差θ14-3至101.7°，AQ节点法接近修正矩阵JAQ的奇异点值，潮流不收敛；同时追踪求解过程中发电机组无功功率和相角差的Q-θ曲线如图4所示，同时绘制AQ节点P-V曲线如图5黑色曲线所示；由图4可知，当节点15和节点2机组无功功率输出分别达到2.0pu和2.4pu时，Q保持恒定，发生了PV节点向PQ节点的转换，负荷进一步增长时的无功均由平衡机组提供；根据图5的结果可知，采用AQ节点法可以直接、准确地计算系统静态电压稳定临界点，而不需要处理一直以来，静态电压稳定性分析中必须面对的牛拉法潮流不收敛的问题。其根本原因是避免了潮流计算中雅克比矩阵在系统静态电压稳定临界点附近的“病态”难题；为更清晰地展示该问题，对矩阵进行奇异值分解，研究矩阵J和JAQ在临界点处奇异值的分布，如图6所示，计算结果按升序排列；由图6的计算结果可见，系统静态电压稳定临界点处，矩阵J的“病态”程度高于JAQ；因此，AQ节点法在临界点处潮流可解，而经典牛拉法潮流无法收敛；
4.求取系统在不同故障方式下的静态电压稳定临界点
为求取系统在不同故障方式下的静态电压稳定临界点，本实施例中研究系统中重要联络线的N-2掉闸故障，采用AQ节点法计算不同方式下系统的P-V曲线，并分析系统的静态电压稳定性，仿真结果如图5所示。根据静态电压稳定临界点得到的负荷裕度如图7所示；图5和图7的计算结果表明，线路6-7由于距离AQ节点较远，其故障情况下负荷裕度与基础方式下的差异最小。线路5-12、线路12-10故障时，由于故障点距离AQ节点距离更近，节点负荷裕度较基础方式明显降低。线路10-8故障时，静态电压稳定裕度最小，由于故障导致平衡机对AQ节点有功传输降低，系统支撑负荷增长的能力减弱。因此综合考虑，系统重要联络线的N-
2掉闸故障下AQ节点的裕度约为266MW；本实施例中，所提出的AQ节点方法对雅克比矩阵J*作了进一步修正，将平衡节点的增量作为其它发电机节增量的基准，同时将AQ节点的负荷增量作为其它节点负荷增量的基准，得到新的雅克比矩阵JAQ，使得计算的迭代次数大大减少。计算过程中也可以直接用J*计算静态电压稳定的临界点，避免对潮流计算中的雅克比矩阵进行修改。计算的基本思路为：每次增加AQ节点和平衡节点的相角差θ时，保持其它负荷节点以及除平衡机外的发电机节点功率不变，计算得到潮流解P*。按计划的模型增加其它负荷和发电机节点，并以P*为初始功率，重新求解潮流，逐步增大θ按此方法迭代，直至AQ节点法不收敛；图8以基础方式为例，对比分析了本实施例中的两种方法和传统CPF方法的性能，其中CPF步长λ＝0.05。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 舵机临界点锁位方法和装置 | 2020-05-11 | 733 |
| 一种金针菇菌种的试管无基与有基临界点育种、提纯方法 | 2020-05-13 | 731 |
| 一种用于不锈钢临界点蚀温度测试的控温实验装置 | 2020-05-12 | 959 |
| 非破坏检查系统及临界点检测系统 | 2020-05-11 | 879 |
| 基于单样本sKLD指标检测复杂生物系统相变临界点的方法 | 2020-05-13 | 513 |
| 一种使用同步磁感线最大临界点的力量制造的设备 | 2020-05-13 | 234 |
| 液晶显示器临界点不良测试方法 | 2020-05-11 | 531 |
| 包含波形临界点控制器和电压传感器的电路断流器 | 2020-05-12 | 923 |
| 舵机临界点锁位方法和装置 | 2020-05-11 | 926 |
| 一种风机沉降到临界点的报警装置 | 2020-05-13 | 80 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
