CN110676841B - 基于直接法的电力电子化电力系统暂态稳定性分析方法 - Google Patents

基于直接法的电力电子化电力系统暂态稳定性分析方法 Download PDF

Info

Publication number
CN110676841B
CN110676841B CN201910862142.1A CN201910862142A CN110676841B CN 110676841 B CN110676841 B CN 110676841B CN 201910862142 A CN201910862142 A CN 201910862142A CN 110676841 B CN110676841 B CN 110676841B
Authority
CN
China
Prior art keywords
power
function
lyapunov
energy function
stability
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.)
Expired - Fee Related
Application number
CN201910862142.1A
Other languages
English (en)
Other versions
CN110676841A (zh
Inventor
车延博
贾晶晶
徐健梅
李明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201910862142.1A priority Critical patent/CN110676841B/zh
Publication of CN110676841A publication Critical patent/CN110676841A/zh
Application granted granted Critical
Publication of CN110676841B publication Critical patent/CN110676841B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及一种基于直接法的电力电子化电力系统暂态稳定性分析方法,包括下列步骤:步骤1:确定系统模型;步骤2:构建李雅普诺夫能量函数:首先根据合理简化确定的系统模型建立其混合势函数P(i,v),其次结合混合势函数理论中的稳定性判别定理,构建系统的李雅普诺夫型能量函数并求其微分值;最后通过能量函数和其微分值的特征判断系统的暂态稳定性;步骤3:估计系统的临界能量。

Description

基于直接法的电力电子化电力系统暂态稳定性分析方法
技术领域
本发明涉及电力电子化电力系统暂态稳定性分析领域,具体来说是基于直接法提出一套完整的暂态稳定性分析流程,对电力电子化电力系统的暂态稳定性进行分析的方法。
背景技术
当前,结合新能源并网发电、柔性直流输电以及微电网等技术的发展,电力电子变流器在电力系统中的应用越来越广泛。在现代电力系统中,发电、输电、变电、配电、用电和储能等各方面基本上都会涉及电力电子设备,电力系统呈现出明显的电力电子化趋势。
电力电子装置在电力系统中大范围的应用,极大地增强了系统的可控性,给系统带来了较高的传输性能。与传统交流电力系统相比,电力电子化电力系统的特点是系统的拓扑随着电力电子器件的开关动作发生变化,整个系统是时变(非自治)的。此外,电力电子变流器的强非线性、动态性以及彼此之间的级联、并联等复杂的连接方式也给电力电子化电力系统的稳定性带来了严重的影响。因此,电力电子化电力系统稳定性分析已成为重要的研究课题。对于电力电子化电力系统的小信号稳定性分析,目前国内外学者已经进行了大量的研究工作。但小信号稳定性仅可以判断系统在稳态工作点附近的稳定情况,无法判断系统稳定域的边界以及平衡点的稳定裕度,难以全面地分析系统的稳定性问题,因此有必要对电力电子化电力系统的暂态稳定性进行研究。
应用于传统交流电力系统暂态稳定性分析的方法主要有时域仿真法(也称逐步积分法)、人工智能法和直接法等。时域仿真法优点是无论待分析的系统多么复杂,系统组成的元件模型多么详细,都可以通过该方法进行系统给定设计方案的检验与性能分析。但时域仿真法无法揭示系统参数与系统性能、稳定性之间的关系,无法得到系统的稳定裕度等定量信息。人工智能法在电力系统动态安全评估方面已有丰富的研究成果,计算速度快,常应用于数据的预处理和后处理。但在分析过程中由于实际数据与预设数据不一致时,应用人工智能法会造成分析结果与实际的稳定指标之间存在偏差。直接法通过构造类似于反映系统能量的标量函数,研究函数随时间变化的趋势,进而判断非线性系统的稳定性。该方法的主要优点是判别速度快,并可以给出系统失稳的模式和程度,但目前仍存在详细模型下能量函数构造困难及判别结果偏保守的局限。电力电子化电力系统和传统的交流电力系统存在很大的差别,但如何运用和改进应用于传统交流电力系统暂态稳定性分析的方法对电力电子化电力系统的暂态稳定性分析具有非常重要的借鉴意义。
发明内容
本发明以直接法为基础,引入混合势函数理论构建系统全局的李雅普诺夫能量函数,提出了电力电子化电力系统暂态稳定性分析方法。本发明的技术方案如下:
一种基于直接法的电力电子化电力系统暂态稳定性分析方法,包括下列步骤:
步骤1:确定系统模型
基于目前对数学模型简化方法的研究,对装置层次的电力电子化电力系统暂态稳定性分析,采用电力电子变流器模型,电网简化为带阻抗的理想电源;对子系统层次和全局电力系统层次的分析,则需要对变流器进行等效建模;
步骤2:构建李雅普诺夫能量函数
首先根据合理简化确定的系统模型建立其混合势函数P(i,v),其次结合混合势函数理论中的稳定性判别定理,构建系统的李雅普诺夫型能量函数并求其微分值;最后通过能量函数和其微分值的特征判断系统的暂态稳定性;具体步骤描述如下:
1)引入混合势函数
首先根据系统中的元件和拓扑关系建立混合势函数P(i,v):
P(i,v)=-A(i)+B(v)+D(i,v) (1)
其中,i,v分别为电路中的电感电流和电容电压;A(i)为电路的电流势函数,B(v)为电路的电压势函数;D(i,v)=iT·γ·v为电路中电容的能量和部分非储能元件的能量,γ为与电路拓扑有关的常系数矩阵;
2)建立李雅普诺夫函数
结合混合势函数P(i,v),令
Figure BDA0002200108150000021
Figure BDA0002200108150000022
则构建系统的李雅普诺夫型能量函数为:
Figure BDA0002200108150000023
其中,L为电路中电感元件的对角矩阵;C为电路中电容元件的对角矩阵;u1为矩阵L-1/2·Aii·L-1/2的最小特征值;u2为矩阵C-1/2·Bvv·C-1/2的最小特征值;
根据式(3)求出其微分值:
Figure BDA0002200108150000031
如果P*(i,v)正定,而
Figure BDA0002200108150000032
是负定的,则这个系统是渐进稳定的,函数P*(i,v)可以作为系统的全局李雅普诺夫能量函数;但该条件仅为判断系统稳定性的充分条件,如果不满足P*(i,v)正定,
Figure BDA0002200108150000033
负定也不能说系统不稳定,通过进一步合理设计参数重新建立满足条件的李雅普诺夫函数;
步骤3:估计系统的临界能量
首先根据电路结构求出系统的稳态平衡点,得出系统处于稳态工作的前提条件;
其次,求解系统的暂态稳定工作条件,结合公式(3),若对电路中所有属于某区域的i,v均有
u1+u2>0 (5)
且当|i|+|v|→∞时,满足
P*(i,v)→∞ (6)
则当时间t→∞时,被研究系统的所有解都会趋于稳态平衡工作点,系统最终能够稳定运行,根据公式(5)就可以求出系统的暂态稳定条件;
综合系统稳态工作的前提条件及其暂态稳定条件,求出系统整体的暂态稳定性判据;
最后估计临界能量,系统处于临界稳定时u1+u2=0,得出临界电压vmin,将vmin带入李雅普诺夫能量函数(3),估计出被研究系统的临界能量。
本发明以直接法为基础,引入混合势函数理论构建电力电子化电力系统全局的李雅普诺夫能量函数,提出了较为具体的电力电子化电力系统暂态稳定性分析流程。同时以电压源型换流器构成的直流输电系统(VSC-HVDC)为实例,分析其暂态稳定性,验证本发明方法的可行性。通过构造一般性的理论和方法对电力电子化电力系统的暂态稳定性分析,使设计过程得以简化,为具体的工程设计提供了一种实用的判定准则。
附图说明
图1:本发明具体流程图
图2:两端VSC-HVDC直流输电系统结构图
图3:两端VSC-HVDC系统简化电路图
图4:两端VSC-HVDC系统电源-负载特性曲线图
具体实施方式
基于直接法的暂态稳定性分析具体流程如图1所示。
本发明以直接法为基础,引入混合势函数理论构建系统全局的李雅普诺夫能量函数,提出了较为具体的电力电子化电力系统暂态稳定性分析方法。本发明的技术方案如下:
步骤1:确定系统模型
基于目前对数学模型简化方法的研究,对装置层次的电力电子化电力系统暂态稳定性分析,采用电力电子变流器模型,电网简化为带阻抗的理想电源;对子系统层次和全局电力系统层次的分析,则需要对变流器进行等效建模。
步骤2:构建李雅普诺夫能量函数
首先根据合理简化确定的系统模型建立其混合势函数P(i,v),其次结合混合势函数理论中的稳定性判别定理,构建系统的李雅普诺夫型能量函数并求其微分值。最后通过能量函数和其微分值的特征判断系统的暂态稳定性。具体步骤描述如下:
1)引入混合势函数
首先根据系统中的元件和拓扑关系建立混合势函数P(i,v):
P(i,v)=-A(i)+B(v)+D(i,v) (1)
其中,i,v分别为电路中的电感电流和电容电压;A(i)为电路的电流势函数,B(v)为电路的电压势函数;D(i,v)=iT·γ·v为电路中电容的能量和部分非储能元件的能量,γ为与电路拓扑有关的常系数矩阵。
其次,通过混合势函数P(i,v)与电路的状态方程之间是否满足公式(2)检验所建立函数的正确性。
Figure BDA0002200108150000041
2)建立李雅普诺夫函数
结合混合势函数P(i,v),令
Figure BDA0002200108150000042
Figure BDA0002200108150000043
则可以构建系统的李雅普诺夫型能量函数为:
Figure BDA0002200108150000051
其中,L为电路中电感元件的对角矩阵;C为电路中电容元件的对角矩阵;u1为矩阵L-1/2·Aii·L-1/2的最小特征值;u2为矩阵C-1/2·Bvv·C-1/2的最小特征值。
根据式(3)求出其微分值:
Figure BDA0002200108150000052
如果P*(i,v)正定,而
Figure BDA0002200108150000053
是负定的,则这个系统是渐进稳定的,该函数P*(i,v)可以作为系统的全局李雅普诺夫能量函数。但该条件仅为判断系统稳定性的充分条件,如果不满足P*(i,v)正定,
Figure BDA0002200108150000054
负定也不能说系统不稳定,此时可以通过进一步合理设计参数重新建立满足条件的李雅普诺夫函数。
步骤3:估计系统的临界能量
首先根据电路结构求出系统的稳态平衡点,得出系统处于稳态工作的前提条件。
其次,求解系统的暂态稳定工作条件,结合公式(3),若对电路中所有属于某区域的i,v均有
u1+u2>0 (5)
且当|i|+|v|→∞时,满足
P*(i,v)→∞ (6)
则当时间t→∞时,被研究系统的所有解都会趋于稳态平衡工作点,系统最终能够稳定运行,根据公式(5)就可以求出系统的暂态稳定条件。
综合系统稳态工作的前提条件及其暂态稳定条件,可以求出系统整体的暂态稳定性判据。
最后估计临界能量,系统处于临界稳定时u1+u2=0,得出临界电压vmin,将vmin带入李雅普诺夫能量函数(3),便可以估计出被研究系统的临界能量如公式(7)所示。
P*(i,v)=minP*(i,vmin) (7)
考虑到电压源型换流器构成的高压直流输电系统(VSC-HVDC)是在电压源型换流器(VSC)技术和门极可关断晶闸管(GTO)及绝缘栅双极晶体管(IGBT)等全控性功率器件基础上发展起来的新型直流输电技术,其中包含了大量的电力电子变换装置,可以实现发、输、配、变、用环节中的特定功能,属于典型的电力电子化电力系统。
故现结合VSC-HVDC系统为实施例、附图对本发明作进一步描述:
步骤1:确定系统模型
如图2所示为两端VSC-HVDC系统结构图。
假定运行方式为左侧的整流器控制直流电压恒定,右侧的逆变器控制直流功率恒定。若不考虑换流器开关特性的影响,并且假定控制器带宽无限高,换流器有足够的响应速度,即整流端可以始终维持直流电压为恒定值,且逆变端的直流功率也可以保持恒定,则可以简化系统结构为如图3所示的理想电路。
Veq为直流输出电压;R、L表示直流线路的电阻和电感;电容C用于模拟逆变器直流侧电容和直流线路充电电容;iL、v分别为流过电感的电流和电容C两端的电容;受控电流源表示负荷恒定为Pcpl的逆变器。
步骤2:构建李雅普诺夫能量函数
依据步骤1确定的系统等效电路模型,得出系统的状态方程为
Figure BDA0002200108150000061
得到整个系统的混合势函数为
Figure BDA0002200108150000062
进一步求得系统的全局李雅普诺夫型能量函数及其微分值分别为
Figure BDA0002200108150000063
Figure BDA0002200108150000064
根据公式(10)、(11)对系统的暂态稳定性进行判定,如果系统的李雅普诺夫能量函数(10)包含平衡点附近大范围内是正定的,且同时其微分值(11)为负定,则判断该系统此时是渐进稳定的,公式(10)即为所求的两端VSC-HVDC系统的李雅普诺夫能量函数。
步骤3:估计系统的临界能量
首先求出系统的稳态平衡点。系统电源特性曲线和负载特性曲线分别如图4中曲线1和曲线2所示,当电压略微升高或略微降低时的电源特性曲线分别为曲线3和曲线4,电源电流等于负载电流时可以得到系统的平衡工作点,即图中交点A,B。如果运行在平衡工作点的系统在微小扰动下能够恢复稳态,则该平衡工作点就是稳定的。显然,图4中A点不是稳态平衡工作点,B点是稳态平衡工作点。
因此可以求得系统的稳态平衡工作点判据为
Figure BDA0002200108150000071
该判据描述了等效电源电压、等效电阻和恒定负荷功率之间的关系,是系统稳定工作的前提条件。
同时,根据所建立的系统混合势函数及其稳定性定理,可以得出系统暂态稳定工作的充分条件是
Figure BDA0002200108150000072
综合平衡点判定及混合势函数稳定性定理,可以得到系统的稳定性判据为
Figure BDA0002200108150000073
同时求得系统的临界能量为
P*(iL,v)=minP*(iL,vmin) (15)
其中,vmin为最小临界电压。

Claims (1)

1.一种基于直接法的电力电子化电力系统暂态稳定性分析方法,包括下列步骤:
步骤1:确定系统模型
对装置层次的电力电子化电力系统暂态稳定性分析,采用电力电子变流器模型,电网简化为带阻抗的理想电源;对子系统层次和全局电力系统层次的分析,则需要对变流器进行等效建模;
步骤2:构建李雅普诺夫能量函数
根据所确定的系统模型建立两端VSC-HVDC系统的混合势函数P(i,v),其次结合混合势函数理论中的稳定性判别定理,构建系统的李雅普诺夫型能量函数并求其微分值;通过能量函数和其微分值的特征判断系统的暂态稳定性;具体步骤描述如下:
1)引入混合势函数
根据两端VSC-HVDC系统中的元件和拓扑关系,建立整个系统的混合势函数:
Figure FDA0003528126680000011
P(iL,v)为混合势函数,Veq为直流输出电压;R、L表示直流线路的电阻和电感;电容C用于模拟逆变器直流侧电容和直流线路充电电容;iL、v分别为流过电感的电流和电容C两端的电压;Pcpl为逆变器恒定负荷;
2)建立李雅普诺夫函数
求得系统的全局李雅普诺夫型能量函数及其微分值分别为
Figure FDA0003528126680000012
Figure FDA0003528126680000013
C是模拟逆变器直流侧电容和直流线路充电电容;vmin为流过电容C两端的电压最小值;
根据公式(2)、(3)对系统的暂态稳定性进行判定,如果系统的李雅普诺夫能量函数(2)包含平衡点附近大范围内是正定的,且同时其微分值(3)为负定,则判断该系统此时是渐进稳定的,公式(2)即用来求得两端VSC-HVDC系统的李雅普诺夫能量函数;
步骤3:估计系统的临界能量
求出系统的稳态平衡点,电源电流等于负载电流时可得到系统的平衡工作点,如果运行在平衡工作点的系统在微小扰动下能够恢复稳态,则该平衡工作点就是稳定的,系统的稳态平衡工作点判据为
Figure FDA0003528126680000021
其中,Rmax为线路阻抗最大值;
该判据描述了等效电源电压、等效电阻和恒定负荷功率之间的关系,是系统稳定工作的前提条件;
根据所建立的系统混合势函数及其稳定性定理,得出系统暂态稳定工作的充分条件是
Figure FDA0003528126680000022
其中,Rmin为线路阻抗最小值;
从而得到系统的稳定性判据为
Figure FDA0003528126680000023
其中,P0为逆变器额定功率;
求得系统的临界能量为
P*(iL,v)=min P*(iL,vmin) (7)
其中,vmin为最小临界电压。
CN201910862142.1A 2019-09-12 2019-09-12 基于直接法的电力电子化电力系统暂态稳定性分析方法 Expired - Fee Related CN110676841B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910862142.1A CN110676841B (zh) 2019-09-12 2019-09-12 基于直接法的电力电子化电力系统暂态稳定性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910862142.1A CN110676841B (zh) 2019-09-12 2019-09-12 基于直接法的电力电子化电力系统暂态稳定性分析方法

Publications (2)

Publication Number Publication Date
CN110676841A CN110676841A (zh) 2020-01-10
CN110676841B true CN110676841B (zh) 2022-06-14

Family

ID=69077846

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910862142.1A Expired - Fee Related CN110676841B (zh) 2019-09-12 2019-09-12 基于直接法的电力电子化电力系统暂态稳定性分析方法

Country Status (1)

Country Link
CN (1) CN110676841B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111914486B (zh) * 2020-08-07 2023-09-15 中国南方电网有限责任公司 一种基于图注意力网络的电力系统暂态稳定评估方法
CN112510688B (zh) * 2020-11-20 2023-01-06 东北电力大学 基于相位校正李雅普诺夫指数的暂态电压稳定监测方法
CN112559960B (zh) * 2020-12-10 2022-06-21 清华大学 微电网的小干扰安全域构建方法及系统
CN113642143A (zh) * 2021-06-16 2021-11-12 南方电网能源发展研究院有限责任公司 电力系统控制方法、装置、计算机设备和存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106294993A (zh) * 2016-08-09 2017-01-04 南京理工大学 一种考虑逆变器电流饱和的暂态能量函数分析方法
CN106451419A (zh) * 2016-09-18 2017-02-22 天津大学 一种飞机电力系统平均值模型及大扰动稳定域构建方法
CN107908829A (zh) * 2017-10-25 2018-04-13 天津大学 基于统一大信号模型的机载电力系统稳定性分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106294993A (zh) * 2016-08-09 2017-01-04 南京理工大学 一种考虑逆变器电流饱和的暂态能量函数分析方法
CN106451419A (zh) * 2016-09-18 2017-02-22 天津大学 一种飞机电力系统平均值模型及大扰动稳定域构建方法
CN107908829A (zh) * 2017-10-25 2018-04-13 天津大学 基于统一大信号模型的机载电力系统稳定性分析方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Stability Analysis of Aircraft Power Systems Based on a Unified Large Signal Model;Yanbo Che;《energies》;20171030;第1-15页 *
电力电子化机载电力系统稳定性分析;车延博;《电力自动化设备》;20180604;第152-156页 *

Also Published As

Publication number Publication date
CN110676841A (zh) 2020-01-10

Similar Documents

Publication Publication Date Title
CN110676841B (zh) 基于直接法的电力电子化电力系统暂态稳定性分析方法
Jiang et al. A conservatism-free large signal stability analysis method for DC microgrid based on mixed potential theory
Al-Saedi et al. Voltage and frequency regulation based DG unit in an autonomous microgrid operation using Particle Swarm Optimization
Patnaik et al. Real-time performance analysis and comparison of various control schemes for particle swarm optimization-based shunt active power filters
Cao et al. Low-frequency stability analysis of inverter-based islanded multiple-bus AC microgrids based on terminal characteristics
CN103916039A (zh) 一种基于反步滑模控制的光伏并网逆变器的控制方法
CN103296905B (zh) 三相电压型功率因数校正变换器的自适应控制方法
Dang et al. Cost function‐based modulation scheme of model predictive control for VIENNA rectifier
CN106300426A (zh) 基于mmc的并网逆变系统的自适应无源性pi控制方法
CN103746583A (zh) 一种mmc的桥臂等效仿真方法
CN105006839A (zh) 并网发电系统的弱电网网源荷阻抗模型分析方法
CN106208782A (zh) 基于模型预测控制的级联h桥光伏逆变器漏电流抑制方法
Canciello et al. Multi-objective adaptive sliding manifold control for More Electric Aircraft
CN111262466B (zh) 基于奇异摄动的模块化多电平变流器的建模控制方法
Molaee et al. A unified power flow approach using VSC-efficiency for AC-DC distribution systems operating at grid connected and islanded modes
Xu et al. Dynamic average-value modeling of three-level T-type grid-connected converter system
Yu et al. Bifurcation analysis of cascaded H-bridge converter controlled by proportional resonant
Amini et al. Interface converter control of distributed generation in microgrids using fractional proportional—Resonant controller
Khan et al. Intelligent transition control approach for different operating modes of photovoltaic inverter
Wang et al. Hierarchical voltage imbalance control for single‐/three‐phase hybrid multimicrogrid
CN106374452A (zh) 一种直流微电网变流器的反馈无源化控制方法
Brito et al. A new three-phase voltage sourced converter laplace model
Gao et al. Comparative Study of Four Droop Control Strategies in DC Microgrid
CN112003271B (zh) 基于分散式阻抗判据的换流器接入交流微电网稳定性分析方法
CN112199914B (zh) 一种电力电子开关恒导纳模型建立方法及系统

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220614