CN108021999A - 一种快速逼近最大负荷功率点的方法及装置 - Google Patents
一种快速逼近最大负荷功率点的方法及装置 Download PDFInfo
- Publication number
- CN108021999A CN108021999A CN201610974527.3A CN201610974527A CN108021999A CN 108021999 A CN108021999 A CN 108021999A CN 201610974527 A CN201610974527 A CN 201610974527A CN 108021999 A CN108021999 A CN 108021999A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- load
- delta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000004364 calculation method Methods 0.000 claims abstract description 58
- 238000005516 engineering process Methods 0.000 claims abstract description 9
- 239000013598 vector Substances 0.000 claims description 25
- 238000012937 correction Methods 0.000 claims description 14
- 238000013459 approach Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000000205 computational method Methods 0.000 abstract 1
- 238000009966 trimming Methods 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011437 continuous method Methods 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供了一种快速逼近最大负荷功率点的方法及装置，其方法按高阶泰勒级数展开法，利用预测‑校正过程追踪潮流方程平衡解曲线计算预测点、确定泰勒展开系数、计算步长逐步逼近最大负荷功率点；其装置包括预测模块、潮流计算模块和解列模块。本发明提供的技术方案基于泰勒级数展开技术的严格变步长计算方法，涉及的计算量相当于计算量不多的潮流计算，计算量小，计算效率高。
Description
技术领域
本发明涉及大电网安全稳定分析领域，具体讲涉及一种快速逼近最大负荷功率点的方法及装置。
背景技术
电压稳定性是指系统维持电压的能力.当负荷导纳增大时,负荷功率亦随之增大,并且功率和电压都是可控的.电压崩溃是指由于电压不稳定导致系统内大面积、大幅度的电压下降的过程。电压稳定性分析则是对这一过程进行理论分析，使得这个过程变得可控。
随着现代电力系统规模逐渐增大，安全稳定问题也日益突出，作为电压稳定分析重要内容的最大负荷功率点(即PV曲线的鼻尖点)的计算，是后续计算负荷裕度及灵敏度等计算的基础。
目前电力系统一般采用的固定步长的预测-校正方法计算传统连续潮流，其计算效率低、速度慢，在大电网分析和在线应用中难以胜任。一般变步长策略采用启发式方法，如简单地根据校正过程迭代次数修正步长，其数学基础存在不足。
发明内容
为获取最大负荷功率点精确值，克服现有技术计算效率低的缺点，本发明提供了一种基于泰勒级数展开技术快速逼近最大负荷功率点的一种严格的变步长计算方法，其计算量仅仅相当于几次潮流计算。
实现本发明的目的技术方案如下：
本发明提供的快速逼近最大负荷功率点的方法，其改进之处在于，所述方法包括：
(1)用高阶泰勒展开预测-校正过程的预测步并计算预测点；
(2)由潮流方程和以预测点为参数的连续潮流扩张方程确定泰勒展开系数；
(3)简化泰勒展开系数，逼近最大负荷功率点。
进一步的，所述步骤(1)中的预测-校正过程如下式(1)扩张方程G(u,λ,Δs)所示：
其中，u和λ：分别为节点电压和负荷因子；Δs为预测步长；f(u,λ):潮流方程；E(u,λ,Δs):表示一维扩张方程。
进一步的，所述步骤(1)中按下式进行所述预测-校正过程的高级泰勒展开预测：
其中，j为当前计算步；(uj,λj)为当前点；为预测点；及分别为高阶泰勒级数展开式系数。
进一步的，所述步骤(2)中泰勒展开系数的确定包括：
(2-1)在直角坐标系下将潮流方程分解为下(3)所示的形式：
f(u,λ):L(u)+q(u,u)+λF＝0 (3)
其中，L：线性部分；q：双线性函数部分；F：负荷变化功率方向；
(2-2)局部参数化的连续潮流扩张方程如下式(4)所示：
其中，ei∈Rn+1是第i个元素为1，其它元素为0的单位坐标向量；[]T：坐标向量的转置；为预测点；u和λ：分别为节点电压和负荷因子。
进一步的，各阶泰勒展开系数如下：
1)当P＝1时，
其中，是向量v的扩展，且v满足 为雅可比矩阵；σ＝±1：
2)当P≥2时，
其中，是向量vP的扩展，vP满足且 为右端向量，为双线性函数； 为单位坐标向量的转置向量。
进一步的，所述步骤(3)中用潮流计算法逼近最大功率负荷点，包括：选择单位坐标向量ei＝en+1＝(0,0,…,1)T，泰勒展开系数如下式所示：
3-1)当P＝1时，泰勒展开系数：
3-2)当P≥2时，泰勒展开系数：
进一步的，按下式计算最大负荷功率点：
其中，K：泰勒展开截断阶数；δ：控制精度；为K阶和K-1阶的预测解。
进一步的，简化计算步长ΔSmax如下式所示：
其中，K为泰勒展开截断阶数；δ为控制精度。
根据最大负荷功率点与当前运行点的向量距离计算负荷裕度并判断负荷裕度是否充足；当判断当前系统的负荷裕度不足时，通过切除指定母线的部分负荷，保持系统的负荷裕度大于门槛值。
本发明提供一种快速逼近最大负荷功率点并根据负荷裕度进行切负荷的装置，所述装置包括：预测模块、潮流计算模块和解列模块
预测模块用高阶泰勒级数展开法展开预测步并计算预测点；
潮流计算模块采用自适应变步长技术快速逼近最大功率负荷点。
解列模块用于切除指定母线的负荷，保持系统具备足够的电压负荷裕度。当判断当前系统的负荷裕度不足时，通过切除母线部分负荷，保持系统的负荷裕度大于门槛值。
与最接近的现有技术比，本发明提供的技术方案具有以下优异效果：
1、本发明提供的技术方案基于泰勒级数展开技术，给出一种严格的变步长计算的方法，可以很快逼近最大负荷功率点，计算量相当于不多几次的潮流计算，减小计算量，提高运算效率，简化运算过程。
2、本发明提供的技术方案中的计算步长是根据泰勒展开系数计算，采用高阶泰勒展开预测方法，减小计算步一般预测解与真实解误差，提高计算精确度。
3、本发明提供的技术方案在发电机无功越限发生PV/PQ节点类型转换时，只需不多几次的迭代修正，且技术方案用自适应变步长技术，计算灵活性和可扩展性强，进一步提高计算的准确性和计算速率。
附图说明
图1为本发明提供的快速逼近最大负荷功率点的方法流程图；
图2为本发明提供的技术方案与固定步长传统方法的对比图；
图3为本发明提供的技术方案计算步长曲线。
具体实施方式
以下将结合说明书附图，以具体实施例的方式详细介绍本发明提供的技术方案。
传统连续潮流计算最大负荷功率点采用固定步长效率低，一般变步长策略采用启发式方法，如简单地根据校正过程迭代次数修正步长，数学基础并不严格。本发明提供的技术方案基于泰勒级数展开技术的严格变步长计算方法，可以很快逼近最大负荷功率点，计算量相当于不多几次的潮流计算，计算量小，计算效率高。
如图1所示，本发明提供的方法具体包括：
a)基于高阶泰勒展开的连续方法
连续方法采用预测-校正(predictor-corrector)过程追踪潮流方程的平衡解曲线，得到PV曲线上的一系列点。预测-校正过程的扩张系统如下所示
其中，u、λ分别为节点电压和负荷因子；Δs为预测步长；f(u,λ):潮流方程；E(u,λ,Δs):表示一维扩张方程。
在预测模块中用高阶泰勒展开预测步，具体方法如下式所示：
其中，j为当前计算步；(uj,λj)为当前点；为预测点；及分别为高阶泰勒级数展开式系数。
以为初始点，采用牛顿方法求解(1)得到(j+1)步计算点。由于采用高阶预测，预测点与真实解比切线预测更接近，校正过程求解线性方程组迭代次数减少，相对于切线预测方法，初始点距离真实解更接近，牛顿迭代次数减小，因此计算量很小。
b)泰勒展开系数的确定
在直角坐标系下，潮流方程分解为如下形式：
f(u,λ):L(u)+q(u,u)+λF＝0 (3)
其中，L为线性部分，q为双线性函数部分，F为负荷变化功率方向。
对于局部参数化的连续潮流，扩张系统式(1)的扩张方程如下所示：
其中，ei∈Rn+1是第i个元素为1，其它元素为0的单位坐标向量；[]T：坐标向量的转置；为预测点；u和λ：分别为节点电压和负荷因子。
式(2)中各阶泰勒展开系数如下式所示：
1)当p＝1时,
其中，是向量v的扩展，且v满足 为雅可比矩阵；σ＝±1：
能保证解曲线路径走向的一致性。
2)当p≥2时，
其中，是向量vP的扩展，vP满足且 为右端向量，为双线性函数； 为单位坐标向量的转置向量。
c)用潮流计算方法逼近最大负荷功率点；
选择ei＝en+1＝(0,0,…,1)T，则泰勒展开系数变为
c-1)当P＝1时，泰勒展开系数变为：
c-2)当P≥2时，泰勒展开系数为：
逼近最大负荷功率点的过程实质上就是不断改变λ而执行常规潮流计算的过程。
计算步长满足下式要求：
其中，K为泰勒展开截断阶数；δ为控制精度；为采用不同阶数(K,K-1)的预测解。
忽略(9)分母中高阶项，可得计算步长为：
本发明提供的技术方案与常规固定步长连续潮流计算不同的是本发明中的计算步长是根据泰勒展开系数计算。由于采用高阶泰勒展开预测方法，除初始运行点执行常规潮流计算需多次迭代求解线性修正方程外，其它计算步一般预测解与真实解差别很小，实际上不需进行潮流计算。即使考虑发电机无功越限发生PV/PQ节点类型转换，方程形式发生变化，也只需不多的几次迭代修正，因此本方法相比其它方法计算量很小。
在逼近最大负荷功率点时，计算步长越来越小，在判断Δsmax小于预先设定的门槛值η时，计算结束，得到最大负荷功率点。根据最大负荷功率点与当前运行点的向量距离计算负荷裕度并判断负荷裕度是否充足；当判断当前系统的负荷裕度不足时，通过切除指定母线的部分负荷，保持系统的负荷裕度大于门槛值。
本发明提供一种快速逼近最大负荷功率点并根据负荷裕度进行切负荷的装置，所述装置包括：预测模块、潮流计算模块和解列模块
预测模块用高阶泰勒级数展开法展开预测步并计算预测点；
潮流计算模块采用自适应变步长技术快速逼近最大功率负荷点。
解列模块用于切除指定母线的负荷，保持系统具备足够的电压负荷裕度。当判断当前系统的负荷裕度不足时，通过切除母线部分负荷，保持系统的负荷裕度大于门槛值。
图2为采用自适应步长变化系列潮流计算与固定步长连续潮流方法逼近最大负荷功率点的对比，计算对象为一实际大电网。
将潮流雅可比矩阵三角分解次数作为计算量的比较依据。传统方法采用Δs＝3的固定步长，需计算18步才能经过最大负荷功率点，潮流雅可比矩阵需分解163次，计算时间为3.5s；而本发明提供的方法需计算14步逼近最大负荷功率点，潮流雅可比矩阵需分解33次，计算时间仅为1.2s。改变传统方法的设定步长虽获得相同的计算结果，但计算时间都不如所提出的自适应变步长技术。
如图3的整个计算过程的计算步长曲线所示，随着计算步的增加，越接近最大负荷功率点，步长越小。在判断步长ΔSmax＜δ时，计算可终止，得到最大负荷功率点。
以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的权利要求保护范围之内。
Claims (10)
1.一种快速逼近最大负荷功率点的方法，其特征在于，所述方法包括：
(1)用高阶泰勒展开预测-校正过程的预测步并计算预测点；
(2)由潮流方程和以预测点为参数的连续潮流扩张方程确定泰勒展开系数；
(3)简化泰勒展开系数，逼近最大负荷功率点。
2.如权利要求1所述的方法，其特征在于，所述步骤(1)中的预测-校正过程如下式(1)扩张方程G(u,λ,Δs)所示：
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>,</mo>
<mi>&lambda;</mi>
<mo>,</mo>
<mi>&Delta;</mi>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>,</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>,</mo>
<mi>&lambda;</mi>
<mo>,</mo>
<mi>&Delta;</mi>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，u和λ：分别为节点电压和负荷因子；Δs为预测步长；f(u,λ):潮流方程；E(u,λ,Δs):表示一维扩张方程。
3.如权利要求1所述的方法，其特征在于，所述步骤(1)中按下式进行所述预测-校正过程的高级泰勒展开预测：
<mrow>
<msub>
<mover>
<mi>u</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
<mi>&Delta;</mi>
<mi>s</mi>
<mo>+</mo>
<mfrac>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;&CenterDot;&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
<mn>6</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>+</mo>
<mo>...</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&lambda;</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
<mi>&Delta;</mi>
<mi>s</mi>
<mo>+</mo>
<mfrac>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>&CenterDot;&CenterDot;&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
<mn>6</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>+</mo>
<mn>...</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，j为当前计算步；(uj,λj)为当前点；为预测点；及分别为高阶泰勒级数展开式系数。
4.如权利要求1所述的方法，其特征在于，所述步骤(2)中泰勒展开系数的确定包括：
(2-1)在直角坐标系下将潮流方程分解为下(3)所示的形式：
f(u,λ):L(u)+q(u,u)+λF＝0 (3)
其中，L：线性部分；q：双线性函数部分；F：负荷变化功率方向；
(2-2)局部参数化的连续潮流扩张方程如下式(4)所示：
<mrow>
<msubsup>
<mi>e</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>-</mo>
<msub>
<mover>
<mi>u</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<mi>&lambda;</mi>
<mo>-</mo>
<msub>
<mover>
<mi>&lambda;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，ei∈Rn+1是第i个元素为1，其它元素为0的单位坐标向量；[]T：坐标向量的转置；为预测点；u和λ：分别为节点电压和负荷因子。
5.如权利要求4所述的方法，其特征在于，各阶泰勒展开系数如下：
1)当P＝1时，
其中，是向量v的扩展，且v满足 为雅可比矩阵；σ＝±1：
2)当P≥2时，
其中，是向量vP的扩展，vP满足且 为右端向量，为双线性函数； 为单位坐标向量的转置向量。
6.如权利要求1所述的方法，其特征在于，所述步骤(3)中用潮流计算法逼近最大功率负荷点，包括：选择单位坐标向量ei＝en+1＝(0,0,…,1)T，泰勒展开系数如下式所示：
3-1)当P＝1时，泰勒展开系数：
3-2)当P≥2时，泰勒展开系数：
7.如权利要求6所述的方法，其特征在于，按下式计算最大负荷功率点：
<mrow>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>u</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>K</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>u</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>u</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>K</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msubsup>
<mi>u</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msup>
<mi>&Delta;s</mi>
<mi>K</mi>
</msup>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
<mi>&Delta;</mi>
<mi>s</mi>
<mo>+</mo>
<mo>...</mo>
<mo>+</mo>
<msubsup>
<mi>u</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msup>
<mi>&Delta;s</mi>
<mi>K</mi>
</msup>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo><</mo>
<mi>&delta;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，K：泰勒展开截断阶数；δ：控制精度；为K阶和K-1阶的预测解。
8.如权利要求7所述的方法，其特征在于，简化计算步长ΔSmax如下式所示：
<mrow>
<msub>
<mi>&Delta;S</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>&delta;</mi>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>j</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<mo>/</mo>
<mo>|</mo>
<mo>|</mo>
<msubsup>
<mi>u</mi>
<mi>j</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>|</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，K为泰勒展开截断阶数；δ为控制精度。
9.如权利要求7所述的方法，其特征在于，根据最大负荷功率点与当前运行点的向量距离计算负荷裕度并判断负荷裕度是否充足；
当判断当前系统的负荷裕度不足时，通过切除指定母线的部分负荷，保持系统的负荷裕度大于门槛值。
10.一种应用权利要求1-9任一所述方法的装置，其特征在于，所述装置包括：预测模块、潮流计算模块和解列模块；
预测模块用高阶泰勒级数展开法展开预测步并计算预测点；
潮流计算模块采用自适应变步长技术快速逼近最大功率负荷点；
解列模块用于切除指定母线的负荷，保持系统具备足够的电压负荷裕度；当判断当前系统的负荷裕度不足时，通过切除母线部分负荷，保持系统的负荷裕度大于门槛值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610974527.3A CN108021999B (zh) | 2016-11-04 | 2016-11-04 | 一种快速逼近最大负荷功率点的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610974527.3A CN108021999B (zh) | 2016-11-04 | 2016-11-04 | 一种快速逼近最大负荷功率点的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108021999A true CN108021999A (zh) | 2018-05-11 |
| CN108021999B CN108021999B (zh) | 2022-02-22 |
Family
ID=62083642
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610974527.3A Active CN108021999B (zh) | 2016-11-04 | 2016-11-04 | 一种快速逼近最大负荷功率点的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108021999B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101582589A (zh) * | 2009-06-18 | 2009-11-18 | 华东电网有限公司 | 基于负荷裕度最大化的有功出力方式优化方法 |
| CN102545207A (zh) * | 2011-12-22 | 2012-07-04 | 河海大学 | 基于预测-校正内点法的vsc-hvdc交直流最优潮流方法 |
| CN103646152A (zh) * | 2013-12-23 | 2014-03-19 | 南方电网科学研究院有限责任公司 | 一种基于矩阵指数的电力系统电磁暂态仿真方法 |
| CN104217090A (zh) * | 2013-06-04 | 2014-12-17 | 刘光晔 | 一种用泰勒级数解析计算电力系统电压稳定临界点的方法 |
| CN104392285A (zh) * | 2014-11-28 | 2015-03-04 | 河海大学 | 一种含混合直流输电的电力系统最优潮流获取方法 |
| CN105762826A (zh) * | 2016-04-15 | 2016-07-13 | 华北电力大学 | 一种含vsc-hvdc交直流系统最大输电能力计算方法 |
-
2016
- 2016-11-04 CN CN201610974527.3A patent/CN108021999B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101582589A (zh) * | 2009-06-18 | 2009-11-18 | 华东电网有限公司 | 基于负荷裕度最大化的有功出力方式优化方法 |
| CN102545207A (zh) * | 2011-12-22 | 2012-07-04 | 河海大学 | 基于预测-校正内点法的vsc-hvdc交直流最优潮流方法 |
| CN104217090A (zh) * | 2013-06-04 | 2014-12-17 | 刘光晔 | 一种用泰勒级数解析计算电力系统电压稳定临界点的方法 |
| CN103646152A (zh) * | 2013-12-23 | 2014-03-19 | 南方电网科学研究院有限责任公司 | 一种基于矩阵指数的电力系统电磁暂态仿真方法 |
| CN104392285A (zh) * | 2014-11-28 | 2015-03-04 | 河海大学 | 一种含混合直流输电的电力系统最优潮流获取方法 |
| CN105762826A (zh) * | 2016-04-15 | 2016-07-13 | 华北电力大学 | 一种含vsc-hvdc交直流系统最大输电能力计算方法 |
Non-Patent Citations (3)
| Title |
|---|
| 容文光: "基于高阶泰勒级数的ATC计算", 《硕士学位论文》 * |
| 杨小煜 等: "利用混合块消去算法的电压稳定负荷裕度灵敏度计算", 《中国电机工程学报》 * |
| 钟浩: "电力系统静态电压稳定性快速评估及其预防控制_钟浩", 《博士学位论文》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108021999B (zh) | 2022-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Alves et al. | Continuation fast decoupled power flow with secant predictor | |
| US20120242152A1 (en) | Method for searching global maximum power point | |
| CN108023356A (zh) | 光伏电站电力系统的电压稳定分析方法及装置 | |
| CN102170131B (zh) | 基于并行计算模式的关联输电断面极限功率计算方法 | |
| CN106249599A (zh) | 一种基于神经网络预测的网络化控制系统故障检测方法 | |
| CN101222140B (zh) | 预想事故电力系统电压稳定快速在线分析与预防控制方法 | |
| CN110401184A (zh) | 多直流馈入受端电网紧急控制优化方法及系统 | |
| CN104156609A (zh) | 基于连续潮流的电网潮流可解性辨识及初值生成方法 | |
| CN111900718B (zh) | 基于多级优化追赶变分迭代法的有源配电网动态仿真方法 | |
| CN104022507A (zh) | 一种直角坐标牛顿法潮流计算方法 | |
| Yang et al. | Asymptotic numerical method for continuation power flow | |
| JP2018133951A (ja) | 太陽光発電量推定装置および方法 | |
| CN103401234B (zh) | 一种基于发电机节点类型的潮流计算方法 | |
| CN105552960A (zh) | 风电场电力系统的电压稳定分析方法及装置 | |
| WO2016067438A1 (ja) | 系統安定化制御装置及び方法 | |
| CN104679968A (zh) | 一种含分布式电源的配电网潮流计算器 | |
| CN108021999B (zh) | 一种快速逼近最大负荷功率点的方法及装置 | |
| CN103326367B (zh) | 一种基于在线电压稳定指标的集中式低压减载方法 | |
| CN110445130B (zh) | 考虑最优无功支撑的静态电压稳定裕度计算装置 | |
| CN103219719B (zh) | 一种利用线路出串运行限制短路电流的分析方法 | |
| CN102593876A (zh) | 一种大规模交直流互联电力系统的连续潮流算法 | |
| CN111740425A (zh) | 一种基于改进型连续潮流法的电力系统静态电压稳定分析方法及系统 | |
| Mori et al. | Continuation Newton-GMRES power flow with linear and nonlinear predictors | |
| CN105356472B (zh) | 基于源荷特性的在线频率综合控制方法 | |
| CN114400662A (zh) | 一种基于连续型潮流计算求解p-v曲线的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |