CN104037756B - 一种含复杂电力设备模型的电力***稳定评估方法 - Google Patents
一种含复杂电力设备模型的电力***稳定评估方法 Download PDFInfo
- Publication number
- CN104037756B CN104037756B CN201410204163.1A CN201410204163A CN104037756B CN 104037756 B CN104037756 B CN 104037756B CN 201410204163 A CN201410204163 A CN 201410204163A CN 104037756 B CN104037756 B CN 104037756B
- Authority
- CN
- China
- Prior art keywords
- equation
- equipment
- power system
- power
- node
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000013097 stability assessment Methods 0.000 title claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 47
- 238000004458 analytical method Methods 0.000 claims abstract description 23
- 238000005516 engineering process Methods 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000002347 injection Methods 0.000 claims description 24
- 239000007924 injection Substances 0.000 claims description 24
- 238000004088 simulation Methods 0.000 claims description 12
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 8
- 238000013461 design Methods 0.000 abstract description 2
- 230000007812 deficiency Effects 0.000 abstract 1
- 238000000205 computational method Methods 0.000 description 8
- 230000008901 benefit Effects 0.000 description 6
- 230000006872 improvement Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种含复杂电力设备模型的电力***稳定评估方法,本方法通过自动微分技术高效计算复杂电力设备对稳定评估计算中雅可比矩阵的影响,通过设计部分自动微分和解耦自动微分策略,克服了传统自动微分技术计算效率差的不足,在代码可维护性和代码效率之间获得了有效的平衡。本方法能够满足大规模电力***稳定评估中对复杂电力设备模型处理的需要,可用于现代电力***的分析、运行与控制等领域。
Description
技术领域
本发明属于电力***的稳定分析与控制技术领域,尤其涉及一种含复杂电力设备模型的电力***稳定评估方法。
背景技术
电力***安全是建设坚强智能电网的核心。动态安全评估是电力***安全分析的重要内容,通过完整的时域仿真评估电力***受扰后的***动态性能;小干扰稳定分析同样是电力***安全分析的另一个有力工具,它能够分析某运行稳态点附近的线性***稳定性,从而分析和预防可能的***振荡。在上述分析本质上均围绕电力***动态模型对应的微分代数方程组开展分析和研究,在电力***时域仿真和小干扰稳定分析中,求解该微分代数方程组的雅可比矩阵尤其重要,是整个分析和计算的基础。然后随着电力电子技术和可再生能源发电技术的发展,越来越多的复杂电力设备接入***,同时这些复杂电力设备多存在非线性强、模型阶数高等特点,为求解上述微分代数方程组的雅可比矩阵带来了困难。
针对上述技术困难,研究者们希望借用先进的计算机软件技术解决复杂电力设备的处理问题。文献《TheRoleofAutomaticDifferentiationinPowerSystemAnalysis》提出了使用自动微分技术计算电力***计算与分析中各类应用的一阶雅可比矩阵和二阶海森矩阵,但其所提出的方法将自动微分软件包视为黑盒,对应的计算效率受到了严重的降低,故无法应用于大规模电力***;发明专利《高效处理复杂电力控制设备的电力***潮流优化方法》提出了改进的自动微分技术并将其应用于电力***稳态分析与优化计算中去,通过应用矩阵叠加原理集成不同种类的复杂电力设备。尽管已有诸多研究成果,但在电力***稳定分析领域中的复杂电力设备模型处理方法上尚无高效的计算方法。本发明基于上述已有研究工作,将自动微分技术应用于电力***稳定分析中雅可比矩阵计算,并针对电力***动态模型微分代数方程组的结构特点,提出了其专有的改进策略,包括部分和解耦的自动微分策略。
发明内容
本发明的目的在于提出一种含复杂电力设备模型的电力***稳定评估方法,能够高效、通用的计算含复杂电力设备模型的电力***动态模型微分代数方程组的雅可比矩阵,可用于电力***时域仿真和小干扰稳定等稳定分析计算应用。
含复杂电力设备模型的电力***稳定评估方法包括如下步骤:
第一步:利用电力***数据采集与监控***采集***,采集电力***各节点的电压和功率数据,通过对***运行状态进行估计,确定***当前稳态运行点,同时获取复杂电力设备列表,对其中每个设备建立其动态方程和电流注入方程。记第k个设备的动态方程为hk,电流注入方程为cx,k和cy,k,其中角标x和y分别代表交流电流的实部与虚部。
第二步:根据***运行数据和每个电力设备的电流注入方程,建立节点电压方程f。将该节点电压方程与每个设备的动态方程联立,得到电力***稳定分析所使用的微分代数方程组F,对应的状态变量为z;
第三步:计算舍去节点注入电流的节点电压方程f对应的雅可比矩阵Jf;计算每个设备电流注入方程cx,k、cy,k对应的雅可比矩阵Jc,k;计算每个设备动态方程hk对应的雅可比矩阵Jh,k。其中Jf为手动计算,Jc,k和Jh,k使用自动微分技术进行计算。
第四步:将上述雅可比矩阵分量(即Jf、Jc,k和Jh,k)进行累加,得到微分代数方程组F对应的雅可比矩阵J。
第五步:利用计算得到的雅可比矩阵进行电力***时域仿真,从而实现对含复杂电力设备模型的电力***稳定评估,对***运行的失稳情况进行预警和校正。
所述的第一步中的第k个设备的动态方程hk和电流注入方程cx,k、cy,k具有如下形式:
其中,zk为第k个设备的内部动态变量,Ux,k和Uy,k为该设备所在节点的电压实部与虚部,t为仿真时间,gx,k和gy,k为该设备注入电流的计算函数的实部与虚部;
所述的第二步中的节点电压方程具有下述形式:
其中Ux和Uy分别为***节点电压向量的实部和虚部;Ix和Iy分别为***节点注入电流向量的实部和虚部;G和B分别为***节点导纳矩阵的实部和虚部。***节点注入电流向量的第i个分量,即节点i的注入电流Ix,i、Iy,i,是对连接在该节点上所有电力设备注入电流的累加,若记其中第k个设备的注入电流为cx,k、cy,k,则有:
通过联立得到的微分代数方程组F及其状态变量z:
F=[fh1…hk…hn]
z=[UxUyz1…zk…zn]
其中,n为***中复杂电力设备的数量。
所述的第四步中的各项雅可比矩阵累加具有下述形式:
J=Jf+∑(Jc,k+Jh,k)
本发明提出了一种含复杂电力设备模型的电力***稳定评估方法。该方法基于改进的自动微分技术,避免了复杂的模型导数计算,使得高效处理复杂电力设备模型成为可能,从而拓展了已有电力***稳定评估的应用范围。与已有的技术相比,本发明提出的方法主要有以下改进:
1、提出了自动微分技术应用于电力***稳定分析雅可比矩阵的计算方法,能够有效的处理大量电力设备模型带来的复杂性;
2、提出了一种部分-解耦自动微分策略,能够极大的提升自动微分的计算效率,同时保证了计算的精确度;
3、所计算得到的***雅可比矩阵可直接用于电力***稳定评估,包括小干扰稳定分析和时域仿真,提升了***运行的安全性。
附图说明
图1是含复杂电力设备模型的电力***稳定评估方法流程图;
图2是稳定分析雅可比矩阵计算中部分自动微分示意图;
图3是稳定分析雅可比矩阵计算中解耦自动微分示意图。
具体实施方式
含复杂电力设备模型的电力***稳定评估方法包括如下步骤:
第一步:利用电力***数据采集与监控***采集***,采集电力***各节点的电压和功率数据,通过对***运行状态进行估计,确定***当前稳态运行点,同时获取复杂电力设备列表,对其中每个设备建立其动态方程和电流注入方程。记第k个设备的动态方程为hk,电流注入方程为cx,k和cy,k,其中角标x和y分别代表交流电流的实部与虚部。
第二步:根据***运行数据和每个电力设备的电流注入方程,建立节点电压方程f。将该节点电压方程与每个设备的动态方程联立,得到电力***稳定分析所使用的微分代数方程组F,对应的状态变量为z;
第三步:计算舍去节点注入电流的节点电压方程f对应的雅可比矩阵Jf;计算每个设备电流注入方程cx,k、cy,k对应的雅可比矩阵Jc,k;计算每个设备动态方程hk对应的雅可比矩阵Jh,k。其中Jf为手动计算,Jc,k和Jh,k使用自动微分技术进行计算。
第四步:将上述雅可比矩阵分量(即Jf、Jc,k和Jh,k)进行累加,得到微分代数方程组F对应的雅可比矩阵J。
第五步:利用计算得到的雅可比矩阵进行电力***时域仿真,从而实现对含复杂电力设备模型的电力***稳定评估,对***运行的失稳情况进行预警和校正。
所述的第一步中的第k个设备的动态方程hk和电流注入方程cx,k、cy,k具有如下形式:
其中,zk为第k个设备的内部动态变量,Ux,k和Uy,k为该设备所在节点的电压实部与虚部,t为仿真时间,gx,k和gy,k为该设备注入电流的计算函数的实部与虚部;
所述的第二步中的节点电压方程具有下述形式:
其中Ux和Uy分别为***节点电压向量的实部和虚部;Ix和Iy分别为***节点注入电流向量的实部和虚部;G和B分别为***节点导纳矩阵的实部和虚部。***节点注入电流向量的第i个分量,即节点i的注入电流Ix,i、Iy,i,是对连接在该节点上所有电力设备注入电流的累加,若记其中第k个设备的注入电流为cx,k、cy,k,则有:
通过联立得到的微分代数方程组F及其状态变量z:
F=[fh1…hk…hn]
z=[UxUyz1…zk…zn]
其中,n为***中复杂电力设备的数量。
所述的第四步中的各项雅可比矩阵累加具有下述形式:
J=Jf+∑(Jc,k+Jh,k)
以下结合附图,对本发明的实施例作详细说明,该发明的流程图如图1所示。
实施例:
为了验证本发明所提出的处理复杂电力设备模型方法的计算效率,发明人使用C++编程语言开发实现了该雅可比矩阵计算程序,并使用一台装配有AMDA8-38502.90GHzCPU和8GB内存的PC机完成了本实施例的测试和验证。在众多电力***稳定分析计算应用中,电力***时域仿真及其灵敏度分析被作为代表选作本实施例的验证应用,因为其中雅可比矩阵计算占用了大量的计算时间,可用于验证本方法的有效性。第三方软件包IDA和ADC04被用于实现时域仿真和自动微分的功能。表1所示的一系列电力***算例被用于本实施例中的测试与验证。
表1:测试算例的***参数
测试*** | 节点 | 发电机 | 线路 | 微分代数方程组维数 |
CASE162 | 162 | 25 | 284 | 892 |
CASE300 | 300 | 69 | 411 | 1422 |
CASE678 | 678 | 170 | 919 | 3194 |
CASE2052 | 2052 | 211 | 2533 | 9170 |
CASE2746 | 2746 | 355 | 3514 | 12520 |
本发明所提出的计算方法主要利用了电力***雅可比矩阵的下述性质:
性质一:电力***雅可比矩阵可以被认为是一个线性时不变部分(网络导纳矩阵)和非线性时变部分(节点注入电流和设备动态方程)的叠加。前部分易于直接矩阵处理,而后一部分不易计算,但由于该部分维数低,可利用自动微分技术予以处理。该计算方法在本发明中称为部分自动微分策略,如图2所示;
性质二:针对性质一中的后一非线性部分,每个设备在其中的分量都是独立的。即每个设备注入电流和动态方程均只与其内部状态变量和所在节点电压有关,因此自动微分可以直接应用于每个设备的注入电流和动态方程的导数求解,从而避免了整体计算,从而提高自动微分的计算效率。该计算方法在本发明中被称为解耦自动微分策略,如图3所示。
除本方法所基于的自动微分技术之外,有限差分方法也被广泛应用于电力***计算中的导数求解和雅可比矩阵计算中。在以下的计算中,我们将展示本方法与未改进的自动微分方法、有限差分方法相比的性能优势。
表2:有限差分、未改进的自动微分与本发明改进的自动微分计算效率对比
如表2所示,本发明所提出的改进的自动微分策略能够有效的提升雅可比矩阵的计算效率,与有限差分和未改进的自动微分方法相比,具有明显的计算效率优势。
与此同时,为了进一步对本发明提出的方法的算法优势进行分析,我们针对CASE678算例开展了详细分析,如表3所示。
表3:改进自动微分计算方法的算法优势对比
根据表3中的计算数据分析,本发明所提出的改进自动微分的计算性能优势主要有下述因素形成:
1)自动微分技术的使用能够充分利用电力***稳定分析中的微分代数方程组中的稀疏结构,而传统的有限差分方法无法利用其稀疏结构,从而只能使用性能较差的稠密解法器求解;
2)与有限差分方法相比,基于自动微分技术的方法能够更精确的求解雅可比矩阵,不存在截断误差。因此更精确的雅可比矩阵导致了更少的牛顿迭代次数,从整体上降低了算法所需的计算量;
3)改进的自动微分策略与未改进的实现相比所需内存更少。“自动微分计算链”是自动微分计算中的最小计算单元,其内存量的需求反映了整个计算应用的内存用量。更少的内存消耗进一步提升了其计算效率。
综上所述,本发明提出的基于改进自动微分策略的复杂电力设备模型计算方法能够有效处理电力***稳定分析中的雅可比矩阵计算问题,与已有的计算方法相比,具有内存需求少、计算效率高等技术优势。
Claims (4)
1.一种含复杂电力设备模型的电力***稳定评估方法,其特征在于,包括如下步骤:
第一步:利用电力***数据采集与监控***采集***,采集电力***各节点的电压和功率数据,通过对***运行状态进行估计,确定***当前稳态运行点,同时获取复杂电力设备列表,对其中每个设备建立其动态方程和电流注入方程;记第k个设备的动态方程为hk,电流注入方程为cx,k和cy,k,其中角标x和y分别代表交流电流的实部与虚部;
第二步:根据***运行数据和每个电力设备的电流注入方程,建立节点电压方程f;将该节点电压方程与每个设备的动态方程联立,得到电力***稳定分析所使用的微分代数方程组F,对应的状态变量为z;
第三步:计算舍去节点注入电流的节点电压方程f对应的雅可比矩阵Jf;计算每个设备电流注入方程cx,k、cy,k对应的雅可比矩阵Jc,k;计算每个设备动态方程hk对应的雅可比矩阵Jh,k;其中Jf为手动计算,Jc,k和Jh,k使用自动微分技术进行计算;
第四步:将上述雅可比矩阵分量(即Jf、Jc,k和Jh,k)进行累加,得到微分代数方程组F对应的雅可比矩阵J;
第五步:利用计算得到的雅可比矩阵进行电力***时域仿真,从而实现对含复杂电力设备模型的电力***稳定评估,对***运行的失稳情况进行预警和校正。
2.根据权利要求1所述的含复杂电力设备模型的电力***稳定评估方法,其特征在于:所述的第一步中的第k个设备的动态方程hk和电流注入方程cx,k、cy,k具有如下形式:
其中,zk为第k个设备的内部动态变量,Ux,k和Uy,k为该设备所在节点的电压实部与虚部,t为仿真时间,gx,k和gy,k为该设备注入电流的计算函数的实部与虚部。
3.根据权利要求1所述的含复杂电力设备模型的电力***稳定评估方法,其特征在于:所述的第二步中的节点电压方程具有下述形式:
其中Ux和Uy分别为***节点电压向量的实部和虚部;Ix和Iy分别为***节点注入电流向量的实部和虚部;G和B分别为***节点导纳矩阵的实部和虚部;***节点注入电流向量的第i个分量,即节点i的注入电流Ix,i、Iy,i,是对连接在该节点上所有电力设备注入电流的累加,若记其中第k个设备的注入电流为cx,k、cy,k,则有:
通过联立得到的微分代数方程组F及其状态变量z:
F=[fh1…hk…hn]
z=[UxUyz1…zk…zn]
其中,n为***中复杂电力设备的数量。
4.根据权利要求1所述的含复杂电力设备模型的电力***稳定评估方法,其特征在于:所述的第四步中的各项雅可比矩阵累加具有下述形式:
J=Jf+Σ(Jc,k+Jh,k)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410204163.1A CN104037756B (zh) | 2014-05-14 | 2014-05-14 | 一种含复杂电力设备模型的电力***稳定评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410204163.1A CN104037756B (zh) | 2014-05-14 | 2014-05-14 | 一种含复杂电力设备模型的电力***稳定评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104037756A CN104037756A (zh) | 2014-09-10 |
CN104037756B true CN104037756B (zh) | 2015-11-18 |
Family
ID=51468420
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410204163.1A Active CN104037756B (zh) | 2014-05-14 | 2014-05-14 | 一种含复杂电力设备模型的电力***稳定评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104037756B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110619162A (zh) * | 2019-09-03 | 2019-12-27 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种联合仿真平台稳定性分析方法及设备 |
CN115758784B (zh) * | 2022-11-30 | 2023-12-12 | 南方电网数字电网研究院有限公司 | 支撑电力***时域仿真的大型雅各比矩阵低耗时迭代方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101478160A (zh) * | 2009-01-12 | 2009-07-08 | 浙江大学 | 高效处理复杂电力控制设备的电力***潮流优化方法 |
CN101707373B (zh) * | 2009-11-20 | 2011-09-14 | 河海大学 | 基于自动微分的电力***状态估计方法 |
-
2014
- 2014-05-14 CN CN201410204163.1A patent/CN104037756B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101478160A (zh) * | 2009-01-12 | 2009-07-08 | 浙江大学 | 高效处理复杂电力控制设备的电力***潮流优化方法 |
CN101707373B (zh) * | 2009-11-20 | 2011-09-14 | 河海大学 | 基于自动微分的电力***状态估计方法 |
Non-Patent Citations (2)
Title |
---|
THE RESEARCH OF PID NEURAL NET WORK DECOUPLING CONTROLLER AND ITS APPLICATION IN UNIT LOAD SYSTEM;HONG-JUN LIU;<Proceedings of the Third International Conference on Machine Learning and Cybemetics>;20040829;505-509 * |
The Role of Automatic Differentiation in Power System Analysis;V.Ajjarapu;<IEEE Transactions on Power System>;19970531;第12卷(第2期);592-597 * |
Also Published As
Publication number | Publication date |
---|---|
CN104037756A (zh) | 2014-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | A data-driven approach to linearize power flow equations considering measurement noise | |
Xu et al. | Robust dispatch of high wind power-penetrated power systems against transient instability | |
CN104269867B (zh) | 一种节点扰动功率转移分布均衡度分析方法 | |
CN101661530A (zh) | 基于相关分析求取风电场稳态等效风速与发电功率的方法 | |
CN105512808A (zh) | 一种面向大数据的电力***暂态稳定评估方法 | |
CN104113061B (zh) | 一种含分布式电源的配电网三相潮流计算方法 | |
CN107681685A (zh) | 一种考虑光伏非线性相关性的电力***概率潮流计算方法 | |
Kaplan et al. | A novel method based on Weibull distribution for short-term wind speed prediction | |
CN104810826A (zh) | 结合拉丁超立方抽样的双向迭代并行概率潮流计算方法 | |
CN103956735B (zh) | 一种分布式发电***的谐波潮流分析方法 | |
CN104156542A (zh) | 一种基于隐式投影的有源配电***稳定性仿真方法 | |
Yang et al. | Monitoring data factorization of high renewable energy penetrated grids for probabilistic static voltage stability assessment | |
CN109412161A (zh) | 一种电力***概率潮流计算方法及*** | |
CN104682387A (zh) | 一种基于多区域交互迭代的概率潮流计算方法 | |
Wang et al. | Applying probabilistic collocation method to power flow analysis in networks with wind farms | |
Ivanov et al. | State estimation for power systems with multilayer perceptron neural networks | |
CN104037756B (zh) | 一种含复杂电力设备模型的电力***稳定评估方法 | |
CN104036118B (zh) | 一种电力***并行化轨迹灵敏度获取方法 | |
Zhang et al. | Low frequency oscillation mode estimation using synchrophasor data | |
Yang et al. | The point estimate method using third-order polynomial normal transformation technique to solve probabilistic power flow with correlated wind source and load | |
Qin et al. | A modified data-driven regression model for power flow analysis | |
CN104158190B (zh) | 电网潮流安全预测装置及方法 | |
CN116316629A (zh) | 一种考虑拓扑特征学习的数据驱动最优潮流计算方法 | |
CN115048857A (zh) | 一种基于cnn的暂稳极限传输功率的确定方法 | |
CN111179110B (zh) | 虚拟电厂可变阶聚合等效鲁棒动态模型建模方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |