CN113515827B - 一种电力系统全过程动态仿真方法及系统 - Google Patents
一种电力系统全过程动态仿真方法及系统 Download PDFInfo
- Publication number
- CN113515827B CN113515827B CN202110434096.2A CN202110434096A CN113515827B CN 113515827 B CN113515827 B CN 113515827B CN 202110434096 A CN202110434096 A CN 202110434096A CN 113515827 B CN113515827 B CN 113515827B
- Authority
- CN
- China
- Prior art keywords
- simulation
- electromechanical
- electromagnetic
- model
- transient
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Networks & Wireless Communication (AREA)
- Operations Research (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种电力系统全过程动态仿真方法及系统,所述方法包括:利用电力系统的元件参数及网络拓扑结构构建电力仿真混合系统,利用电磁暂态仿真模型和机电暂态仿真模型同步进行电力系统全过程的混合仿真,在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,控制所述电磁暂态仿真模型退出仿真,所述机电暂态仿真模型继续进行电力系统全过程的仿真,在所述电力仿真混合系统进入中期动态过程后,对状态变量满足预设的机电稳态条件的一部分机电仿真子模型进行断开操作,对状态变量不满足预设的机电稳态条件的另一部分所述机电仿真子模型进行向前仿真。本发明能够通用于电力系统不同暂态过程的仿真分析,实现仿真效率、规模和速度的统一。
Description
技术领域
本发明涉及电力系统仿真领域,尤其涉及一种交直流混联电力系统全过程动态仿真方法及系统。
背景技术
根据电力系统受到扰动后的动态过程作用时间的长短,暂态过程通常可以分为三部分:电磁暂态过程、机电暂态过程和中长期过程。电力系统受扰动后的电磁、机电和中长期动态过程是相互联系和相互影响的,并没有明确的界线。因此,电力系统从电磁暂态、机电暂态到中长期动态过程是一个连续的过程,电磁过程、机电暂态过程、中长期过程对电力系统不同时间尺度下的刻画,相互之间相辅相成,往往是密不可分的。
电磁暂态仿真受限于仿真规模和仿真速度,仅能对局部的区域进行精细化的模拟,难以对完整的交直流混联系统进行仿真,且仿真速度过慢,不利于进行多故障扫描,而中长期动态仿真、机电暂态仿真可以对大规模电网进行较快的仿真计算,但均无法精确模拟现代电力系统中越来越多的电力电子设备的动态行为,其难以反映直流的换相失败过程,对电网的运行和控制造成较大的误差。上述各种仿真分析技术都各自存在缺陷,若能进一步融合各种分析技术,对于不同的仿真局部采用不同的仿真分析技术,集成各种分析技术的优势,实现全过程动态的一体化仿真,具有重大的现实意义。
发明内容
本发明提供一种电力系统全过程动态仿真方法及系统,以解决现有技术不具备适用于不同暂态过程的融合性仿真技术,本发明能够通用于电力系统不同暂态过程的仿真分析,实现仿真效率、规模和速度的统一。
本发明实施例提供一种电力系统全过程动态仿真方法,包括如下步骤:
利用电力系统的元件参数及网络拓扑结构构建电力仿真混合系统,获得所述电力系统的初始运行状态参数;其中,所述电力仿真混合系统包括电磁暂态仿真模型和机电暂态仿真模型,所述机电暂态仿真模型由电磁仿真子模型和多个机电仿真子模型组成,所述初始运行状态参数包括所述电力系统各个节点的初始电压和初始注入功率;
根据所述初始运行状态参数,利用所述电磁暂态仿真模型和所述机电暂态仿真模型同步进行电力系统全过程的混合仿真,并输出所述电磁暂态仿真模型的电磁暂态仿真结果、所述机电暂态仿真模型的机电暂态仿真结果;
在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,控制所述电磁暂态仿真模型退出仿真,所述机电暂态仿真模型继续进行电力系统全过程的仿真,所述电磁仿真子模型输出电磁暂态仿真子结果,所述机电仿真子模型输出机电暂态仿真子结果;
在所述电力仿真混合系统进入中期动态过程后,对状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型进行断开操作,对状态变量不满足预设的机电稳态条件的另一部分所述机电仿真子模型进行向前仿真直至仿真结束,并输出所述电力仿真混合系统的仿真结果;其中,所述仿真结果包括所述电力仿真混合系统各个节点的电压、功率、发电机功角曲线和发电机相对摇摆角曲线。
优选地,在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,所述电磁暂态仿真模型向所述机电暂态仿真模型发送稳态进入信号;
所述机电暂态仿真模型响应于所述稳态进入信号,向所述电磁暂态仿真模型反馈退出信号,以使所述电磁暂态仿真模型退出仿真。
优选地,所述预设的电磁稳态进入条件如公式(1)所示:
|Is,n+1-Is,n|<ε1,n=i,i+1,...,i+m1 (1)
式中,Is,n+1、Is,n分别为第n+1步和第n步的所述电磁暂态仿真模型注入到所述机电暂态仿真模型的注入电流,ε1为第一预设偏差判定阈值,m1为满足所述预设的电磁稳态进入条件的连续步数的个数,需要连续m1个步长注入电流值均满足所述公式(1),电磁暂态计算达到稳态。
优选地,所述状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型包括励磁机模型;
对应所述励磁机模型的预设的机电暂态稳态条件如公式(2)所示:
|EFD,n+1-EFD,n|<ε2,n=i,i+1,...,i+m2 (2)
式中,EFD,n+1、EFD,n分别为第n+1步和第n步的励磁电压,ε2为第二预设偏差判定阈值,m2为满足所述预设的机电稳态条件的连续步数的个数,需要连续m2个步长励磁电压值均满足所述公式(2),励磁电压暂态计算达到稳态。
优选地,所述机电暂态仿真结果和所述电磁暂态仿真结果通过数据交互接口进行等值数据交互通信。
相应地,本发明实施例还提供了一种电力系统全过程动态仿真系统,包括:
构建单元,用于利用电力系统的元件参数及网络拓扑结构构建电力仿真混合系统,获得所述电力系统的初始运行状态参数;其中,所述电力仿真混合系统包括电磁暂态仿真模型和机电暂态仿真模型,所述机电暂态仿真模型由电磁仿真子模型和多个机电仿真子模型组成,所述初始运行状态参数包括所述电力系统各个节点的初始电压和初始注入功率;
混合仿真单元,用于根据所述初始运行状态参数,利用所述电磁暂态仿真模型和所述机电暂态仿真模型同步进行电力系统全过程的混合仿真,并输出所述电磁暂态仿真模型的电磁暂态仿真结果、所述机电暂态仿真模型的机电暂态仿真结果;
第一机电单元,用于在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,控制所述电磁暂态仿真模型退出仿真,所述机电暂态仿真模型继续进行电力系统全过程的仿真,所述电磁仿真子模型输出电磁暂态仿真子结果,所述机电仿真子模型输出机电暂态仿真子结果;
第二机电单元,用于在所述电力仿真混合系统进入中期动态过程后,对状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型进行断开操作,对状态变量不满足预设的机电稳态条件的另一部分所述机电仿真子模型进行向前仿真直至仿真结束,并输出所述电力仿真混合系统的仿真结果;其中,所述仿真结果包括所述电力仿真混合系统各个节点的电压、功率、发电机功角曲线和发电机相对摇摆角曲线。
优选地,在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,所述电磁暂态仿真模型向所述机电暂态仿真模型发送稳态进入信号;
用于所述机电暂态仿真模型响应于所述稳态进入信号,向所述电磁暂态仿真模型反馈退出信号,以使所述电磁暂态仿真模型退出仿真。
优选地,所述预设的电磁稳态进入条件如公式(1)所示:
|Is,n+1-Is,n|<ε1,n=i,i+1,...,i+m1 (1)
式中,Is,n+1、Is,n分别为第n+1步和第n步的所述电磁暂态仿真模型注入到所述机电暂态仿真模型的注入电流,ε1为第一预设偏差判定阈值,m1为满足所述预设的电磁稳态进入条件的连续步数的个数,需要连续m1个步长注入电流值均满足所述公式(1),电磁暂态计算达到稳态。
优选地,所述状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型包括励磁机模型;
对应所述励磁机模型的预设的机电暂态稳态条件如公式(2)所示:
|EFD,n+1-EFD,n|<ε2,n=i,i+1,...,i+m2 (2)
式中,EFD,n+1、EFD,n分别为第n+1步和第n步的励磁电压,ε2为第二预设偏差判定阈值,m2为满足所述预设的机电稳态条件的连续步数的个数,需要连续m2个步长励磁电压值均满足所述公式(2),励磁电压暂态计算达到稳态。
优选地,所述机电暂态仿真结果和所述电磁暂态仿真结果通过数据交互接口进行等值数据交互通信。
本发明实施例提供的电力系统全过程动态仿真方法及系统,通过构建电磁暂态仿真模型和机电暂态仿真模型,对不同的仿真局部采用不同的仿真分析技术,通过融合电磁暂态仿真、机电暂态仿真和中长期动态仿真,实现全过程动态的一体化仿真。本发明实施例能够通用于电力系统不同暂态过程的仿真分析,实现仿真效率、规模和速度的统一。
附图说明
图1是本发明实施例提供的一种仿真平台架构说明示意图;
图2是本发明实施例提供的一种电力系统全过程动态仿真方法流程示意图;
图3是本发明实施例提供的一种电力系统全过程动态仿真系统结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参见图1,是本发明实施例提供的一种仿真平台架构说明示意图。
本发明实施例提供的仿真平台架构,包括:电磁暂态仿真程序、机电暂态仿真程序和电磁-机电暂态混合仿真接口。具体的,电磁-机电暂态混合仿真接口包括机电和电磁侧的数据交互接口和等值数据交互通信设计两个部分,机电暂态仿真和电磁暂态仿真通过数据交互接口进行等值数据交互通信。以下结合实施例对本发明提供的电力系统全过程动态仿真方法进行具体说明。
参见图2,是本发明实施例提供的一种电力系统全过程动态仿真方法流程示意图,所述方法包括以下步骤:
S1,利用电力系统的元件参数及网络拓扑结构构建电力仿真混合系统,获得所述电力系统的初始运行状态参数;其中,所述电力仿真混合系统包括电磁暂态仿真模型和机电暂态仿真模型,所述机电暂态仿真模型由电磁仿真子模型和多个机电仿真子模型组成,所述初始运行状态参数包括所述电力系统各个节点的初始电压和初始注入功率;
需要说明的是,在本发明实施例中,将电力仿真混合划分为电磁暂态仿真模型和机电暂态仿真模型,电磁暂态仿真模型用于对应刻画电力系统中局部动态响应的特性,机电暂态仿真模型用于对应刻画电力系统中其他如暂态稳定或静态稳定的稳定特性。
S2,根据所述初始运行状态参数,利用所述电磁暂态仿真模型和所述机电暂态仿真模型同步进行电力系统全过程的混合仿真,并输出所述电磁暂态仿真模型的电磁暂态仿真结果、所述机电暂态仿真模型的机电暂态仿真结果;
需要说明的是,利用所述电磁暂态仿真模型和所述机电暂态仿真模型同步进行电力系统全过程的混合仿真是以机电-电磁暂态混合仿真为出发点,开展初始化、故障期间及故障恢复过程中的动态仿真。在机电-电磁暂态混合仿真中进行机电侧仿真时,对电磁侧采用注入电流源进行等效。如果电磁侧的注入电流进入稳态,电磁暂态仿真将退出,仿真全部由机电暂态仿真完成。在电磁暂态仿真达到稳态退出后电磁侧的部分将变为机电暂态仿真,因此该部分的模型需要在初始时即在机电暂态仿真中建立相应的模型并同步开展仿真,其依据电磁暂态侧计算结果不断进行状态变量的调整,只是该部分机电暂态仿真的结果不与系统中的其他部分进行交互;在电磁暂态仿真退出后,该部分的机电暂态模型开展起作用,取代电磁暂态仿真与系统中的其他部分进行交互(注入电流形式)。
由于在电磁暂态仿真中,导纳阵是积分步长的函数,电磁暂态仿真方法只能采用固定步长。因此在机电暂态和电磁暂态的混合仿真阶段,机电暂态仿真方法不采用变步长积分算法,采用定步长隐式梯形积分算法。优选的,在混合仿真阶段,电磁仿真方法采用固定步长积分算法,机电暂态仿真采用定步长隐式梯形积分算法。
S3,在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,控制所述电磁暂态仿真模型退出仿真,所述机电暂态仿真模型继续进行电力系统全过程的仿真,所述电磁仿真子模型输出电磁暂态仿真子结果,所述机电仿真子模型输出机电暂态仿真子结果;
需要说明的是,在电磁暂态仿真结束,采用机电暂态仿真进行全系统的仿真的阶段是系统剧烈变化阶段结束后的相对变快阶段,由于系统变量变化仍比较快,采用变步长积分算法的情况下,误差控制和计算速度的矛盾比较突出。如果误差控制严格,步长容易过小,例如小于1ms,会使此暂态阶段仿真速度很慢;反之,如果误差控制条件放宽,虽然在此阶段的步长控制和计算速度可能比较合适,但在随后的慢变(中长期动态)过程中,由于累积误差会很大,校正迭代中往往不收敛,需要强制减小步长,计算量反而增加很多。因此,在此阶段采用定步长加速牛顿算法进行机电暂态仿真。此时系统模型需考虑常规机电暂态仿真中包含的模型,包括发电机及其各类控制器、各类型直流输电系统、各类FACTS(FlexibleAlternative Current Transmission Systems,柔性交流输电系统)装置、负荷等;还包括中长期动态仿真中需考虑的元件模型,包括继电保护和自动装置模型、发电厂机组的过励和低励限制等限制器模型、锅炉汽机协调控制模型、水力系统模型、核反应堆模型、变电站有载调压模型、AGC模型、AVC模型等。
其中,采用定步长加速牛顿算法的机电暂态仿真进行全系统的仿真,包括:
利用电力系统的元件参数和网络拓扑结构,形成描述所述电力系统机电暂态过程的第一微分-代数方程组;其中,所述第一微分-代数方程组如下式所示:
式中,微分方程表示电力系统元件的动态特性,是电力系统的状态方程;代数方程表示电力系统元件的静态特性,主要是系统的网络方程;x为n个状态向量,V为m个代数向量,微分代数方程组的阶数为n和m之和;Y是电力系统的导纳矩阵。
利用所述电力系统初始化时的运行参数,获得初始仿真结果,其中,所述初始仿真结果包括状态变量的初始值和网络变量的初始值;
基于牛顿积分算法,根据所述初始仿真结果计算所述第一微分-代数方程组的雅可比矩阵,并在相邻两次迭代的雅可比矩阵解的残差小于第一预设阈值时,不更新雅可比矩阵,获得所述雅可比矩阵、所述第一微分-代数方程组的第一状态变量和所述第一微分-代数方程组的第一网络变量;
利用所有设备的注入电流,形成描述所述电力系统机电暂态过程的第二微分-代数方程组;其中,所述电力系统机电暂态过程的第二微分-代数方程组如下式所示:
式中,AG、BG、CG和YG是雅克比矩阵,雅可比矩阵AG、BG、CG和YG是(x,V)的函数;YGn是电力系统所有设备对应于非凸极效应的常数矩阵。
基于牛顿积分算法,根据所述雅可比矩阵、所述第一状态变量和所述第一网络变量对所述第二微分-代数方程组进行迭代求解,在相邻两次迭代的状态变量的残差小于第二预设阈值,且相邻两次迭代的网络变量的残差小于第三预设阈值时,计算收敛,进行下一步机电仿真计算,直至完成此阶段的机电暂态仿真。
S4,在所述电力仿真混合系统进入中期动态过程后,对状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型进行断开操作,对状态变量不满足预设的机电稳态条件的另一部分所述机电仿真子模型进行向前仿真直至仿真结束,并输出所述电力仿真混合系统的仿真结果;其中,所述仿真结果包括所述电力仿真混合系统各个节点的电压、功率、发电机功角曲线和发电机相对摇摆角曲线。
需要说明的是,在系统仿真进入中期后,系统相对慢变阶段,快速变化的状态变量例如励磁机、PSS(Power System Stabilize,电力系统稳定器)、电力电子设备等都进入相对稳定阶段,在此阶段将电力电子设备作为恒定注入电流,并忽略励磁机、PSS等装置的影响,以提高系统的计算速度,缩短仿真计算时间。因此,可以判断励磁机模型和PSS模型的状态变量是否满足预设的机电稳态条件。
在本发明实施例中,优选地,在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,所述电磁暂态仿真模型向所述机电暂态仿真模型发送稳态进入信号;
所述机电暂态仿真模型响应于所述稳态进入信号,向所述电磁暂态仿真模型反馈退出信号,以使所述电磁暂态仿真模型退出仿真。在一个具体实施例中,机电暂态仿真程序监测电磁暂态仿真是否满足预设的电磁稳态进入条件。若不满足,机电暂态仿真和电磁暂态仿真都进行到下一时刻;若满足,所述电磁暂态仿真模型向所述机电暂态仿真模型发送稳态进入信号,所述机电暂态仿真模型响应于所述稳态进入信号,向所述电磁暂态仿真模型反馈电压幅值为9999.0的退出信号,这样当电磁仿真程序接收到这个9999.0的电压幅值信号后,所述电磁暂态仿真模型退出仿真。
在本发明实施例中,优选地,所述预设的电磁稳态进入条件如公式(1)所示:
|Is,n+1-Is,n|<ε1,n=i,i+1,...,i+m1 (1)
式中,Is,n+1、Is,n分别为第n+1步和第n步的所述电磁暂态仿真模型注入到所述机电暂态仿真模型的注入电流,ε1为第一预设偏差判定阈值,m1为满足所述预设的电磁稳态进入条件的连续步数的个数,需要连续m1个步长注入电流值均满足所述公式(1),电磁暂态计算达到稳态。
需要说明的是,在机电-电磁暂态混合仿真中进行机电侧仿真时,所述电磁暂态仿真模型满足预设的电磁稳态进入条件即当电磁侧的电流注入满足预设的稳态进入条件,电磁暂态仿真退出,仿真全部由机电暂态仿真完成。
在本发明实施例中,优选地,所述状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型包括励磁机模型;
对应所述励磁机模型的预设的机电暂态稳态条件如公式(2)所示:
|EFD,n+1-EFD,n|<ε2,n=i,i+1,...,i+m2 (2)
式中,EFD,n+1、EFD,n分别为第n+1步和第n步的励磁电压,ε2为第二预设偏差判定阈值,m2为满足所述预设的机电稳态条件的连续步数的个数,需要连续m2个步长励磁电压值均满足所述公式(2),励磁电压暂态计算达到稳态。需要说明的是,因为PSS的输出是加在励磁机的输入上,因此不必对PSS做单独判断。在一个具体实施例中,机电暂态仿真自动向前仿真并监测励磁和PSS是否满足预设的机电稳态条件。若不满足,机电暂态仿真进行到下一时刻;若满足,励磁和PSS已进入稳态,将励磁机模型和PSS模型断开,用这一时刻的励磁电压并保持其恒定,机电暂态仿真自动向前仿真直到仿真结束。
在本发明实施例中,优选地,所述机电暂态仿真结果和所述电磁暂态仿真结果通过数据交互接口进行等值数据交互通信。
参见图3,是本发明实施例提供的一种电力系统全过程动态仿真系统结构示意图,所述系统包括:
构建单元11,用于利用电力系统的元件参数及网络拓扑结构构建电力仿真混合系统,获得所述电力系统的初始运行状态参数;其中,所述电力仿真混合系统包括电磁暂态仿真模型和机电暂态仿真模型,所述机电暂态仿真模型由电磁仿真子模型和多个机电仿真子模型组成,所述初始运行状态参数包括所述电力系统各个节点的初始电压和初始注入功率;
混合仿真单元12,用于根据所述初始运行状态参数,利用所述电磁暂态仿真模型和所述机电暂态仿真模型同步进行电力系统全过程的混合仿真,并输出所述电磁暂态仿真模型的电磁暂态仿真结果、所述机电暂态仿真模型的机电暂态仿真结果;
第一机电单元13,用于在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,控制所述电磁暂态仿真模型退出仿真,所述机电暂态仿真模型继续进行电力系统全过程的仿真,所述电磁仿真子模型输出电磁暂态仿真子结果,所述机电仿真子模型输出机电暂态仿真子结果;
第二机电单元14,用于在所述电力仿真混合系统进入中期动态过程后,对状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型进行断开操作,对状态变量不满足预设的机电稳态条件的另一部分所述机电仿真子模型进行向前仿真直至仿真结束,并输出所述电力仿真混合系统的仿真结果;其中,所述仿真结果包括所述电力仿真混合系统各个节点的电压、功率、发电机功角曲线和发电机相对摇摆角曲线。
优选地,在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,所述电磁暂态仿真模型向所述机电暂态仿真模型发送稳态进入信号;
用于所述机电暂态仿真模型响应于所述稳态进入信号,向所述电磁暂态仿真模型反馈退出信号,以使所述电磁暂态仿真模型退出仿真。
优选地,所述预设的电磁稳态进入条件如公式(1)所示:
|Is,n+1-Is,n|<ε1,n=i,i+1,...,i+m1 (1)
式中,Is,n+1、Is,n分别为第n+1步和第n步的所述电磁暂态仿真模型注入到所述机电暂态仿真模型的注入电流,ε1为第一预设偏差判定阈值,m1为满足所述预设的电磁稳态进入条件的连续步数的个数,需要连续m1个步长注入电流值均满足所述公式(1),电磁暂态计算达到稳态。
优选地,所述状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型包括励磁机模型;
对应所述励磁机模型的预设的机电暂态稳态条件如公式(2)所示:
|EFD,n+1-EFD,n|<ε2,n=i,i+1,...,i+m2 (2)
式中,EFD,n+1、EFD,n分别为第n+1步和第n步的励磁电压,ε2为第二预设偏差判定阈值,m2为满足所述预设的机电稳态条件的连续步数的个数,需要连续m2个步长励磁电压值均满足所述公式(2),励磁电压暂态计算达到稳态。
优选地,所述机电暂态仿真结果和所述电磁暂态仿真结果通过数据交互接口进行等值数据交互通信。
可以理解的是,上述的所述装置的实施例的相关内容还可以参考上述的对应的所述方法的实施例的内容,在此不做赘述。
本发明实施例提供的电力系统全过程动态仿真方法及系统,通过构建电磁暂态仿真模型和机电暂态仿真模型,对不同的仿真局部采用不同的仿真分析技术,通过融合电磁暂态仿真、机电暂态仿真和中长期动态仿真,实现全过程动态的一体化仿真。本发明实施例能够通用于电力系统不同暂态过程的仿真分析,实现仿真效率、规模和速度的统一。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。
Claims (10)
1.一种电力系统全过程动态仿真方法,其特征在于,包括如下步骤:
利用电力系统的元件参数及网络拓扑结构构建电力仿真混合系统,获得所述电力系统的初始运行状态参数;其中,所述电力仿真混合系统包括电磁暂态仿真模型和机电暂态仿真模型,所述机电暂态仿真模型由电磁仿真子模型和多个机电仿真子模型组成,所述初始运行状态参数包括所述电力系统各个节点的初始电压和初始注入功率;
根据所述初始运行状态参数,利用所述电磁暂态仿真模型和所述机电暂态仿真模型同步进行电力系统全过程的混合仿真,并输出所述电磁暂态仿真模型的电磁暂态仿真结果、所述机电暂态仿真模型的机电暂态仿真结果;
在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,控制所述电磁暂态仿真模型退出仿真,所述机电暂态仿真模型继续进行电力系统全过程的仿真,所述电磁仿真子模型输出电磁暂态仿真子结果,所述机电仿真子模型输出机电暂态仿真子结果;
在所述电力仿真混合系统进入中期动态过程后,对状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型进行断开操作,对状态变量不满足预设的机电稳态条件的另一部分所述机电仿真子模型进行向前仿真直至仿真结束,并输出所述电力仿真混合系统的仿真结果;其中,所述仿真结果包括所述电力仿真混合系统各个节点的电压、功率、发电机功角曲线和发电机相对摇摆角曲线。
2.如权利要求1所述的电力系统全过程动态仿真方法,其特征在于,在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,所述电磁暂态仿真模型向所述机电暂态仿真模型发送稳态进入信号;
所述机电暂态仿真模型响应于所述稳态进入信号,向所述电磁暂态仿真模型反馈退出信号,以使所述电磁暂态仿真模型退出仿真。
3.如权利要求1或2任一项所述的电力系统全过程动态仿真方法,其特征在于,所述预设的电磁稳态进入条件如公式(1)所示:
|Is,n+1-Is,n|<ε1,n=i,i+1,...,i+m1 (1)
式中,Is,n+1、Is,n分别为第n+1步和第n步的所述电磁暂态仿真模型注入到所述机电暂态仿真模型的注入电流,ε1为第一预设偏差判定阈值,m1为满足所述预设的电磁稳态进入条件的连续步数的个数,需要连续m1个步长注入电流值均满足所述公式(1),电磁暂态计算达到稳态。
4.如权利要求1所述的电力系统全过程动态仿真方法,其特征在于,所述状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型包括励磁机模型;
对应所述励磁机模型的预设的机电暂态稳态条件如公式(2)所示:
|EFD,n+1-EFD,n|<ε2,n=i,i+1,...,i+m2 (2)
式中,EFD,n+1、EFD,n分别为第n+1步和第n步的励磁电压,ε2为第二预设偏差判定阈值,m2为满足所述预设的机电稳态条件的连续步数的个数,需要连续m2个步长励磁电压值均满足所述公式(2),励磁电压暂态计算达到稳态。
5.如权利要求1所述的电力系统全过程动态仿真方法,其特征在于,所述机电暂态仿真结果和所述电磁暂态仿真结果通过数据交互接口进行等值数据交互通信。
6.一种电力系统全过程动态仿真系统,其特征在于,包括:
构建单元,用于利用电力系统的元件参数及网络拓扑结构构建电力仿真混合系统,获得所述电力系统的初始运行状态参数;其中,所述电力仿真混合系统包括电磁暂态仿真模型和机电暂态仿真模型,所述机电暂态仿真模型由电磁仿真子模型和多个机电仿真子模型组成,所述初始运行状态参数包括所述电力系统各个节点的初始电压和初始注入功率;
混合仿真单元,用于根据所述初始运行状态参数,利用所述电磁暂态仿真模型和所述机电暂态仿真模型同步进行电力系统全过程的混合仿真,并输出所述电磁暂态仿真模型的电磁暂态仿真结果、所述机电暂态仿真模型的机电暂态仿真结果;
第一机电单元,用于在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,控制所述电磁暂态仿真模型退出仿真,所述机电暂态仿真模型继续进行电力系统全过程的仿真,所述电磁仿真子模型输出电磁暂态仿真子结果,所述机电仿真子模型输出机电暂态仿真子结果;
第二机电单元,用于在所述电力仿真混合系统进入中期动态过程后,对状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型进行断开操作,对状态变量不满足预设的机电稳态条件的另一部分所述机电仿真子模型进行向前仿真直至仿真结束,并输出所述电力仿真混合系统的仿真结果;其中,所述仿真结果包括所述电力仿真混合系统各个节点的电压、功率、发电机功角曲线和发电机相对摇摆角曲线。
7.如权利要求6所述的电力系统全过程动态仿真系统,其特征在于,在所述电磁暂态仿真模型满足预设的电磁稳态进入条件之时,所述电磁暂态仿真模型向所述机电暂态仿真模型发送稳态进入信号;
用于所述机电暂态仿真模型响应于所述稳态进入信号,向所述电磁暂态仿真模型反馈退出信号,以使所述电磁暂态仿真模型退出仿真。
8.如权利要求6或7任一项所述的电力系统全过程动态仿真系统,其特征在于,所述预设的电磁稳态进入条件如公式(1)所示:
|Is,n+1-Is,n|<ε1,n=i,i+1,...,i+m1 (1)
式中,Is,n+1、Is,n分别为第n+1步和第n步的所述电磁暂态仿真模型注入到所述机电暂态仿真模型的注入电流,ε1为第一预设偏差判定阈值,m1为满足所述预设的电磁稳态进入条件的连续步数的个数,需要连续m1个步长注入电流值均满足所述公式(1),电磁暂态计算达到稳态。
9.如权利要求6所述的电力系统全过程动态仿真系统,其特征在于,所述状态变量满足预设的机电稳态条件的一部分所述机电仿真子模型包括励磁机模型;
对应所述励磁机模型的预设的机电暂态稳态条件如公式(2)所示:
|EFD,n+1-EFD,n|<ε2,n=i,i+1,...,i+m2 (2)
式中,EFD,n+1、EFD,n分别为第n+1步和第n步的励磁电压,ε2为第二预设偏差判定阈值,m2为满足所述预设的机电稳态条件的连续步数的个数,需要连续m2个步长励磁电压值均满足所述公式(2),励磁电压暂态计算达到稳态。
10.如权利要求6所述的电力系统全过程动态仿真系统,其特征在于,所述机电暂态仿真结果和所述电磁暂态仿真结果通过数据交互接口进行等值数据交互通信。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110434096.2A CN113515827B (zh) | 2021-04-21 | 2021-04-21 | 一种电力系统全过程动态仿真方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110434096.2A CN113515827B (zh) | 2021-04-21 | 2021-04-21 | 一种电力系统全过程动态仿真方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113515827A CN113515827A (zh) | 2021-10-19 |
CN113515827B true CN113515827B (zh) | 2022-12-06 |
Family
ID=78062665
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110434096.2A Active CN113515827B (zh) | 2021-04-21 | 2021-04-21 | 一种电力系统全过程动态仿真方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113515827B (zh) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103995734A (zh) * | 2014-06-18 | 2014-08-20 | 中国南方电网有限责任公司电网技术研究中心 | 基于rtds的电力系统混合实时仿真系统及仿真方法 |
CN104732033A (zh) * | 2015-03-31 | 2015-06-24 | 河南行知专利服务有限公司 | 一种电力系统机电暂态过程的仿真分析方法 |
CN112001057B (zh) * | 2020-06-30 | 2024-01-23 | 中国电力科学研究院有限公司 | 基于序分量提高电力系统混合仿真稳定性的方法及装置 |
-
2021
- 2021-04-21 CN CN202110434096.2A patent/CN113515827B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113515827A (zh) | 2021-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mohagheghi et al. | Optimal wide area controller and state predictor for a power system | |
Azmy et al. | Artificial neural network-based dynamic equivalents for distribution systems containing active sources | |
CN103810646B (zh) | 一种基于改进投影积分算法的有源配电系统动态仿真方法 | |
US10971931B2 (en) | Decentralized hardware-in-the-loop scheme | |
CN108873733B (zh) | 电力信息物理系统中信息预想事故影响的分析方法 | |
US20200091764A1 (en) | System and method for a fast power network simulator | |
Shafiee-Rad et al. | Robust decentralized voltage control for uncertain DC microgrids | |
Shi et al. | Data-driven model-free adaptive damping control with unknown control direction for wind farms | |
Chen et al. | Security assessment for intentional island operation in modern power system | |
CN113515827B (zh) | 一种电力系统全过程动态仿真方法及系统 | |
Suvire et al. | Wind farm: Dynamic model and impact on a weak power system | |
Agarwal et al. | Continuously differentiable analytical models for implicit control within power flow | |
CN114188945B (zh) | 一种含光伏电源的配电网短路电流计算方法及装置 | |
JP2008040034A (ja) | 電力系統のシミュレーション方法 | |
KR101848993B1 (ko) | 급전원 훈련 시뮬레이터 | |
van der Meer et al. | Computationally efficient transient stability modeling of multi-terminal VSC-HVDC | |
Tarasov et al. | The use of a topological list for calculating the parameters of electric-power systems on the basis of the node-voltage method | |
CN103729574B (zh) | 海上油田群电网预想操作的在线仿真计算方法 | |
Brandl et al. | HIL methods supporting the development process from simulations to real environment testing | |
Lino | Robust recurrent neural network-based dynamic equivalencing in power system | |
Sowa et al. | Dynamic equivalents for calculation of power system restoration | |
Lv et al. | Power System Load Restoration Based on Data-Driven Constraints Linearization Method | |
Sacasqui et al. | Simulation Study of a GPC Controller for a Hydroelectric Plant during Electric Faults | |
CN117458622A (zh) | 一种基于CatBoost的暂态电压安全域在线修正方法及平台 | |
Lee et al. | Reinforcement Learning Approach to Velocity and Position Control of Metro Trains |
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 |