CN113949104A - 一种多类型电力电子设备馈入系统动态稳定判断方法 - Google Patents

一种多类型电力电子设备馈入系统动态稳定判断方法 Download PDF

Info

Publication number
CN113949104A
CN113949104A CN202111054374.8A CN202111054374A CN113949104A CN 113949104 A CN113949104 A CN 113949104A CN 202111054374 A CN202111054374 A CN 202111054374A CN 113949104 A CN113949104 A CN 113949104A
Authority
CN
China
Prior art keywords
subsystem
feed
model
power electronic
electronic equipment
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
Application number
CN202111054374.8A
Other languages
English (en)
Other versions
CN113949104B (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.)
Shanghai Electric Power University
Original Assignee
Shanghai Electric Power 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 Shanghai Electric Power University filed Critical Shanghai Electric Power University
Priority to CN202111054374.8A priority Critical patent/CN113949104B/zh
Publication of CN113949104A publication Critical patent/CN113949104A/zh
Application granted granted Critical
Publication of CN113949104B publication Critical patent/CN113949104B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48—Controlling the sharing of the in-phase component
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/50—Controlling the sharing of the out-of-phase component
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/22—The renewable source being solar energy
    • H02J2300/24—The renewable source being solar energy of photovoltaic origin
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • 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
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/50—Photovoltaic [PV] energy
    • Y02E10/56—Power conversion systems, e.g. maximum power point trackers

Abstract

本发明涉及一种多类型电力电子设备馈入系统动态稳定判断方法。根据改进动态短路比IMDSCR模型对于多类型电力电子设备馈入系统建立的子系统进行判断，IMDSCR模型在传统单馈入短路比的基础上增加了模型系数，不同子系统模型系数不同，但判断方法相同。IMDSCR模型的建立方法为根据子系统建立离散时间状态模型，并建立李雅普诺夫函数，最后基于压缩映射定理，建立IMDSCR模型。与现有技术相比，本发明可同时判断所有子系统的稳定性从而得到总系统的稳定性，具有可避免多类型电力电子设备之间动态交互影响，参数具有实时性等优点。

Description

一种多类型电力电子设备馈入系统动态稳定判断方法
技术领域
本发明涉及电力系统控制领域，尤其是涉及一种多类型电力电子设备馈入系统动态稳定判断方法。
背景技术
随着科技的发展，以光伏和风电为代表的新能源将会被大量开发，进而使得电力系统电力电子化程度快速变高。随着馈入受端电网的风场和光伏场的规模增大，受端电网的支撑强度逐渐下降，使得风场或光伏场与受端电网在并网点的交互作用加强，受端电网的灵活多变的运行方式势必会影响到风场或光伏场并网的稳定性。另外，馈入的风场或光伏场与距离较近本地柔性高压直流存在复杂的动态交互作用，会进一步恶化系统动态稳定性。因此如何准确评估多类型电力电子设备多馈入系统动态稳定性显得至关重要。
在理论研究和工程应用中，通常使用单馈入短路比(short circuit ratio，SCR)来评估受端电网对于外接设备的支撑强度。在实际运行过程中，受端电网的灵活多变的运行方式(如切负荷、断线等)，受端电网的等值电感不再是恒定值，而是时变量，使得传统单馈入短路比很难应用于在线评估系统稳定性。在弱电网环境下，受端电网的变化，会导致风场、光伏电场、柔性高压直流以及受端电网之间的交互作用更加复杂。基于传统短路比所计算的实际短路比的临界值会受到上述因素的影响而发生变化，难以准确判断系统临界稳定。另外，现有的静态解耦系统模型难以考虑到受端电网出现大扰动情况，使得等效模型不准确，进而使得指标计算不准确。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种多类型电力电子设备馈入系统动态稳定判断方法，以解决现有技术无法考虑到手段电网运行方式的变化以及多类型电力电子设备间动态交互的影响。
本发明的目的可以通过以下技术方案来实现：
一种多类型电力电子设备馈入系统动态稳定判断方法，将多类型电力电子设备馈入系统中每个子系统的实时信息，分别输入至改进动态短路比模型IMDSCR并进行判断，判断每个子系统的改进动态短路比模型IMDSCR的输出数值：
若所有子系统的改进动态短路比模型IMDSCR输出数值均大于1，则多类型电力电子设备馈入系统动态稳定；
若存在子系统的改进动态短路比模型IMDSCR输出数值小于1，则判断该子系统是否满足第一条件：若是，则多类型电力电子设备馈入系统动态稳定；若否，则多类型电力电子设备馈入系统动态不稳定；
所述第一条件具体如下：在M时刻，存在一个正数λ，使得以下表达式成立：
Figure BDA0003253994070000021
Figure BDA0003253994070000022
式中，
Figure BDA0003253994070000023
分别是增益函数值和通用解值，
Figure BDA0003253994070000024
表示增益系数；
所述子系统包括单VSC-HVDC馈入子系统、单光伏场馈入子系统和单直驱风场馈入子系统，所述单光伏场馈入子系统和所述单直驱风场馈入子系统的系统参数相同，所述改进动态短路比模型IMDSCR表达式为：
Figure BDA0003253994070000025
式中，Kk,r表示模型系数，在所有下标为k，r的参数中，k表示单VSC-HVDC馈入子系统的编号，r表示单光伏场馈入子系统或单直驱风场馈入子系统的编号，逗号均表示或者，M表示每个子系统迭代次数的最小公倍数，
Figure BDA0003253994070000026
均表示增益系数，SCRk,r表示子系统的单馈入短路比。
进一步地，通用解值
Figure BDA0003253994070000027
表达式为：
Figure BDA0003253994070000031
Figure BDA0003253994070000032
Figure BDA0003253994070000033
Figure BDA0003253994070000034
uk＝[veqαk veqβk vgαk vgβk Pvk Qvk Prefk Qrefk]T
ur＝[veqαr veqβr iDCRr iDCIr vgαr vgβr VDCr Ig2r]T
式中，
Figure BDA0003253994070000035
均为计算系数矩阵，其中最后的下标数字第一位数字代表矩阵的行，第二位数字代表矩阵的列，M表示每个子系统迭代次数的最小公倍数，w表示采样时刻，uk,r(w)表示输入系数；veqαk、veqβk、vgαk和vgβk分别表示第k个单VSC-HVDC馈入子系统的动态等效电源瞬时电压和并网点瞬时电压在两相静止坐标系下的α、β轴分量；Pref、Qref分别表示第k个单VSC-HVDC馈入子系统定功率外环控制的有功参考值、无功参考值；Pv、Qv分别表示第k个单VSC-HVDC馈入子系统电力电子设备输出的有功、无功；veqαr、veqβr、vgαr和vgβr分别表示第r个单光伏场馈入子系统或单直驱风场馈入子系统的动态等效电源瞬时电压和并网点瞬时电压在两相静止坐标系下的α、β轴分量；VDCr、Ig2r分别表示第r个单光伏场馈入子系统或单直驱风场馈入子系统直流外环控制输入参考值，电流内环q轴控制输入参考值；iDCR、iDCI分别表示第r个单光伏场馈入子系统或单直驱风场馈入子系统整流侧注入直流侧电容节点电流、从直流侧电容节点流出的电流。
进一步地，所述增益系数
Figure BDA0003253994070000041
的表达式为：
Figure BDA0003253994070000042
Figure BDA0003253994070000049
Figure BDA0003253994070000043
式中，
Figure BDA0003253994070000044
均为计算系数矩阵，其中最后的下标数字第一位数字代表矩阵的行，第二位数字代表矩阵的列，M表示每个子系统迭代次数的最小公倍数，w表示采样时刻。
进一步地，所述计算系数矩阵
Figure BDA0003253994070000045
的表达式为：
Figure BDA0003253994070000046
Figure BDA0003253994070000047
Figure BDA0003253994070000048
Figure BDA0003253994070000051
Figure BDA0003253994070000052
式中，UgN、LeqN、PEN分别表示并网点额定电压，理想变压器额定电感、电力电子设备馈入的额定有功；下标k表示第k个单VSC-HVDC馈入子系统，下标r表示第r个单光伏场馈入子系统或单直驱风场馈入子系统，
Figure BDA0003253994070000053
Figure BDA0003253994070000054
均表示子系统系数。
进一步地，所述子系统系数
Figure BDA0003253994070000055
Figure BDA0003253994070000056
的表达式为：
Figure BDA0003253994070000057
Figure BDA0003253994070000058
Figure BDA0003253994070000059
Figure BDA00032539940700000510
Figure BDA0003253994070000061
Figure BDA0003253994070000062
式中，
Figure BDA0003253994070000063
分别表示外环功率控制比例积分环节的积分系数、比例系数，内环电流控制比例积分环节的积分系数、比例系数，锁相环形成角速度环节比例积分环节的积分系数、比例系数，Lc、Leq分别表示滤波电感、理想变压器电感，CDC表示直流侧电容；下标k代表第k个单VSC-HVDC馈入子系统，下标r表示第r个单光伏场馈入子系统或单直驱风场馈入子系统，Ts表示采样周期，ω0是基准角频率。
进一步地，所述改进动态短路比模型IMDSCR的建立方法如下：
S1、建立子系统动态模型；
S2、根据子系统动态模型，建立不含控制器的连续域主电路动态方程，离散化后得到主电路离散状态模型；
S3、建立控制器的离散状态方程；
S4、根据步骤S2和S3的两个离散状态模型建立子系统的离散时间状态模型；
S5、根据所述离散时间状态模型建立李雅普诺夫函数；
S5、基于压缩映射定理，根据所述李雅普诺夫函数建立IMDSCR模型。
进一步地，所述单VSC-HVDC馈入子系统在模型是基于并网点WAMS设备实时数据建立的。
进一步地，所述李雅普诺夫函数的建立步骤如下：
SA1、通过数学归纳法，求解所述离散时间状态模型在M时刻的解；
SA2、根据步骤SA1中的解，求得在M时刻的通用解值；
SA3、对于步骤SA2中得到的通用解值，选取1-范数作为测度，结合增益函数，得到李雅普诺夫函数。
进一步地，所述增益函数表达式为：
Figure BDA0003253994070000071
式中，
Figure BDA0003253994070000072
表示增益系数。
进一步地，所述离散时间状态模型表达式为：
Figure BDA0003253994070000073
Figure BDA0003253994070000074
下标k代表第k个单VSC-HVDC馈入子系统，下标r表示第r个单光伏场馈入子系统或单直驱风场馈入子系统，
Figure BDA0003253994070000075
Figure BDA0003253994070000076
均表示子系统系数，ukr(w)表示输入系数，w表示采样时刻，
Figure BDA0003253994070000077
表示不同子系统的状态变量。
与现有技术相比，本发明具有以下优点：
1、对于单VSC-HVDC馈入子系统和第r个单光伏场馈入子系统或单直驱风场馈入子系统统一建立改进动态短路比模型IMDSCR并根据模型参数的不同，分别进行判断，以判断子系统的稳定性来判断多类型电力电子设备馈入系统整体的稳定性，考虑了多类型电力电子设备之间动态交互影响，提高了模型判断的准确性。
2、本发明中的改进动态短路比模型使用的参数均为子系统的实时信息，对于系统的规划设计保证了实时性。
附图说明
图1为本发明的方法流程示意图。
图2为本发明子系统等效过程示意图。
图3为本发明第k个单VSC-HVDC馈入子系统。
图4为本发明第r个单光伏场或者直驱风场馈入子系统。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
本实施例提供了一种多类型电力电子设备馈入系统动态稳定判断方法，将多类型电力电子设备馈入系统中每个子系统的实时信息，分别输入至改进动态短路比模型IMDSCR并进行判断，判断每个子系统的改进动态短路比模型IMDSCR的输出数值：
若所有子系统的改进动态短路比模型IMDSCR输出数值均大于1，则多类型电力电子设备馈入系统动态稳定；
若存在子系统的改进动态短路比模型IMDSCR输出数值小于1，则判断该子系统是否满足第一条件：若是，则多类型电力电子设备馈入系统动态稳定；若否，则多类型电力电子设备馈入系统动态不稳定。
该模型的建立方法如下，对于不同的子系统建立不同的离散时间状态模型，子系统包括单VSC-HVDC馈入子系统、单光伏场馈入子系统和单直驱风场馈入子系统。
首先，以第k个单VSC-HVDC馈入子系统为例，其并网点母线装有WAMS设备，等效过程示意图如图2所示。在仿真环境中，通过用一个理想可调变压器和一个理想发电机来代替从第k个母线向电网看进去的剩余系统。理想可调变压器要满足其电阻为0，电感较小的条件。理想发电机要满足其内部电感为0，电阻接近于0，惯量较大，以保证初始电压幅值为1，相角为0。
基于WAMS实测数据，基于公式(1)，通过不断地改变理想变压的变比和移相角，使得在仿真环境下所得的并网点电压幅值、相角和基于WAMS所测得的电压幅值、相角相等，且馈入到并网点的有功、无功与实际测量的有功无功相同。倘若满足上述条件，则理想变压串联理想发电机这一等效电路可以近似模拟从并网点看进去的剩余系统的动态行为，可得到如下表达式：
Figure BDA0003253994070000091
其中，ηk,r，
Figure BDA0003253994070000092
veq(k,r)，δeq(k,r)，
Figure BDA0003253994070000093
和
Figure BDA0003253994070000094
分别是每一时刻下的理想变压器变比、移相角，理想发电机电压幅值、相角，并网点WAMS设备记录的电压幅值、相角。“k,r”，分别代表VSC-HVDC、新能源场站的编号k∈{1,2...,m},r∈{m+1,...,m+r1,m+r1+1,...,m+r1+r2}，光伏场站编号对应着{m+1,...,m+r1}风场编号对应着{m+r1,...,m+r1+r2}，逗号表示或者的关系。
然后建立第k个单VSC-HVDC馈入等效系统离散状态模型，第k个单VSC-HVDC馈入等效系统模型如图3所示。
其中，所有变量下标含k，代表第k个，下面不再赘述；Lc、Leq分别是滤波电感、理想变压器电感；ig代表流过滤波电感、理想变压器电感的瞬时电流；vc、vg、veq分别代表电力电子设备输出瞬时电压，并网点瞬时电压，动态等效电源瞬时电压；vcd c、vgd c、Vcd、vcq c、vgq c、Vcq分别是电力电子设备输出瞬时电压、并网点瞬时电压、控制器输出的调制电压经过park变换(即在锁相环提供的d-q旋转坐标系下，以及转子角度提供的d-q旋转坐标系)后的d轴、q轴量；igd c、igq c分别是流过滤波电感、理想变压器电感和理想变压器电阻的瞬时电流经过park变换(即在锁相环提供的d-q旋转坐标系下，以及转子角度提供的d-q旋转坐标系)后的d轴、q轴量；ωpll、θpll是锁相环环节产生的角频率、相角；Pref、Qref分别是定功率外环控制的有功参考值、无功参考值；Pv、Qv分别是电力电子设备输出的有功、无功；PI(z)离散比例积分环节和离散比例微分环节，其下标i、O、pll分别代表内环电流控制、外环功率控制、锁相环。
根据图3，可以得到静止坐标系下的不含控制器的连续域主电路动态方程，如下式所示：
Figure BDA0003253994070000101
其中，igα、vcα、veqα、igβ、vcβ、veqβ分别是流过滤波电感、理想变压器电感和理想变压器电阻的瞬时电流、电力电子设备输出瞬时电压、动态等效电源瞬时电压在两相静止坐标系下的α、β轴分量；Leq(t)是理想变压器的时变电感；
基于精确离散化方法，可将上式离散化，可得主电路离散状态模型，如下式所示：
Figure BDA0003253994070000102
其中，
Figure BDA0003253994070000103
uk，
Figure BDA0003253994070000104
的表达式如下式所示：
Figure BDA0003253994070000105
Figure BDA0003253994070000106
其中，w是采样时刻，Ts＝1/fs，fs为采样频率，其具体数值由实际系统中的采样设备确定。
基于前向欧拉法，可将VSC-HVDC的控制器进行离散化，可获得控制器的离散状态方程如下式所示：
Figure BDA0003253994070000107
其中，
Figure BDA0003253994070000112
如下式所示
Figure BDA0003253994070000113
Figure BDA0003253994070000114
Figure BDA0003253994070000115
其中，
Figure BDA0003253994070000116
分别是外环功率控制比例积分环节的积分系数、比例系数，内环电流控制比例积分环节的积分系数、比例系数，锁相环形成角速度环节比例积分环节的积分系数、比例系数；下标k代表第k个单VSC-HVDC馈入子系统。则第k个单VSC-HVDC馈入等效系统离散状态模型如下式所示：
Figure BDA0003253994070000117
上式可改写为：
Figure BDA0003253994070000118
其中，xk(w)，uk(w)，Θk(w)，Ψk(w)，
Figure BDA0003253994070000119
的表达式如下式所示：
Figure BDA0003253994070000121
其中，UgN、LeqN、PEN分别是并网点额定电压，理想变压器额定电感、电力电子设备馈入的额定有功；下标k代表第k个单VSC-HVDC馈入字系统。
然后建立第r个单光伏场或者直驱风场馈入等效系统离散状态模型，第r个单光伏场或者直驱风场馈入等效系统模型如图4所示。
其中，上图中所有变量下标含r，代表第r个；VDC、Ig2分别式直流外环控制输入参考值，电流内环q轴控制输入参考值；其他变量解释在之前已作了解释。
首先根据图4，得到静止坐标系下的不含控制器的连续域主电路动态方程，如下式所示：
Figure BDA0003253994070000122
其中，vDC是直流侧电容电压，iDCR、iDCI分别是整流侧注入直流侧电容节点电流、从直流侧电容节点流出的电流。
基于精确离散化方法，可将上式离散化，可得主电路离散状态模型，如下式所示：
Figure BDA0003253994070000123
其中，
Figure BDA0003253994070000124
uk的表达式如下式所示：
Figure BDA0003253994070000131
Figure BDA0003253994070000132
由于光伏场站和直驱风场采用相同的控制结构，因此，基于前向欧拉法，可将其控制器进行离散化，可获得控制器的离散状态方程如下式所示：
Figure BDA0003253994070000133
其中，
Figure BDA0003253994070000134
如下式所示：
Figure BDA0003253994070000141
可得到第r个单光伏场或者直驱风场馈入等效系统离散状态模型如下：
Figure BDA0003253994070000142
可上式改写为
Figure BDA0003253994070000143
其中，xr(w)，ur(w)，Θr(w)，Ψr(w)，
Figure BDA0003253994070000144
的表达式如下式所示
Figure BDA0003253994070000145
其中，PEN代表光伏电场或者直驱风场馈送到受端电网的有功。与前文关于该变量的解释本质是一样的，都是电力电子设备馈送到受端电网的有功。
由于公式(11)和公式(20)在物理意义上是相同的，都是对于单个电力电子设备馈入动态系统的离散状态模型，整体的数学表达式也类似。为了给出通式，因此，在此之后，将公式公式(11)和公式(20)一般化，如下式所示：
Figure BDA0003253994070000151
上式中出现的下标“k,r”代表VSC-HVDC、光伏电场或直驱风场的编号，逗号代表或者的意思。对整个系统而言，M是每个单馈入系统有限步迭代次数Mk,r的最小公倍数(k∈{1,2...,m},r∈{m+1,...,m+r1,m+r1+1,...,m+r1+r2})。上式在第M个时刻的通解为下式：
Figure BDA0003253994070000152
通过选取合适的基底矩阵
Figure BDA0003253994070000153
和
Figure BDA0003253994070000154
可将上式中的时变矩阵Θk,r(w)和Ψk,r(w)表达成下式
Figure BDA0003253994070000155
将上式代入公式(23)可得下式
Figure BDA0003253994070000156
Figure BDA0003253994070000157
Figure BDA0003253994070000158
进一步可从上式得到在M时刻的通用解值如下式所示
Figure BDA0003253994070000161
Figure BDA0003253994070000162
选取1-范数作为测度，且根据范数相关性质，可得最大形式的李雅普诺夫函数，如下式所示
Figure BDA0003253994070000163
其中增益函数
Figure BDA0003253994070000164
的表达式如下式所示
Figure BDA0003253994070000165
其中，
Figure BDA0003253994070000166
的表达式如下式所示
Figure BDA0003253994070000171
为了保证系统满足李雅普诺夫意义稳定，需要满足以下两个条件：
条件1：
Figure BDA0003253994070000172
和
Figure BDA0003253994070000173
在M时刻分别小于1的，则可以推导出
Figure BDA0003253994070000174
条件2：保证
Figure BDA0003253994070000175
和
Figure BDA0003253994070000176
的乘积小于1，则可以推导出
Figure BDA0003253994070000177
因此，可得到改进动态短路比IMDSCR的表达式为
Figure BDA0003253994070000178
第一条件：在M时刻，存在一个正数λk，使得下式成立
Figure BDA0003253994070000179
其中，
Figure BDA0003253994070000181
分别是增益函数和子系统解，i＝1,2。
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。

Claims (10)

1.一种多类型电力电子设备馈入系统动态稳定判断方法，其特征在于，将多类型电力电子设备馈入系统中每个子系统的实时信息，分别输入至改进动态短路比模型IMDSCR并进行判断，判断每个子系统的改进动态短路比模型IMDSCR的输出数值：
若所有子系统的改进动态短路比模型IMDSCR输出数值均大于1，则多类型电力电子设备馈入系统动态稳定；
若存在子系统的改进动态短路比模型IMDSCR输出数值小于1，则判断该子系统是否满足第一条件：若是，则多类型电力电子设备馈入系统动态稳定；若否，则多类型电力电子设备馈入系统动态不稳定；
所述第一条件具体如下：在M时刻，存在一个正数λ，使得以下表达式成立：
Figure FDA0003253994060000011
Figure FDA0003253994060000012
式中，
Figure FDA0003253994060000013
分别是增益函数值和通用解值，
Figure FDA0003253994060000014
表示增益系数；
所述子系统包括单VSC-HVDC馈入子系统、单光伏场馈入子系统和单直驱风场馈入子系统，所述单光伏场馈入子系统和所述单直驱风场馈入子系统的系统参数相同，所述改进动态短路比模型IMDSCR表达式为：
Figure FDA0003253994060000015
式中，Kk,r表示模型系数，在所有下标为k，r的参数中，k表示单VSC-HVDC馈入子系统的编号，r表示单光伏场馈入子系统或单直驱风场馈入子系统的编号，逗号均表示或者，M表示每个子系统迭代次数的最小公倍数，
Figure FDA0003253994060000016
均表示增益系数，SCRk,r表示子系统的单馈入短路比。
2.根据权利要求1所述的一种多类型电力电子设备馈入系统动态稳定判断方法，其特征在于，通用解值
Figure FDA0003253994060000017
表达式为：
Figure FDA0003253994060000021
Figure FDA0003253994060000022
Figure FDA0003253994060000023
Figure FDA0003253994060000024
uk＝[veqαk veqβk vgαk vgβk Pvk Qvk Prefk Qrefk]T
ur＝[veqαr veqβr iDCRr iDCIr vgαr vgβr VDCr Ig2r]T
式中，
Figure FDA0003253994060000025
均为计算系数矩阵，其中最后的下标数字第一位数字代表矩阵的行，第二位数字代表矩阵的列，M表示每个子系统迭代次数的最小公倍数，w表示采样时刻，uk,r(w)表示输入系数；veqαk、veqβk、vgαk和vgβk分别表示第k个单VSC-HVDC馈入子系统的动态等效电源瞬时电压和并网点瞬时电压在两相静止坐标系下的α、β轴分量；Pref、Qref分别表示第k个单VSC-HVDC馈入子系统定功率外环控制的有功参考值、无功参考值；Pv、Qv分别表示第k个单VSC-HVDC馈入子系统电力电子设备输出的有功、无功；veqαr、veqβr、vgαr和vgβr分别表示第r个单光伏场馈入子系统或单直驱风场馈入子系统的动态等效电源瞬时电压和并网点瞬时电压在两相静止坐标系下的α、β轴分量；VDCr、Ig2r分别表示第r个单光伏场馈入子系统或单直驱风场馈入子系统直流外环控制输入参考值，电流内环q轴控制输入参考值；iDCR、iDCI分别表示第r个单光伏场馈入子系统或单直驱风场馈入子系统整流侧注入直流侧电容节点电流、从直流侧电容节点流出的电流。
3.根据权利要求1所述的一种多类型电力电子设备馈入系统动态稳定判断方法，其特征在于，所述增益系数
Figure FDA0003253994060000031
的表达式为：
Figure FDA0003253994060000032
Figure FDA0003253994060000033
Figure FDA0003253994060000034
式中，
Figure FDA0003253994060000035
均为计算系数矩阵，其中最后的下标数字第一位数字代表矩阵的行，第二位数字代表矩阵的列，M表示每个子系统迭代次数的最小公倍数，w表示采样时刻。
4.根据权利要求3所述的一种多类型电力电子设备馈入系统动态稳定判断方法，其特征在于，所述计算系数矩阵
Figure FDA0003253994060000036
的表达式为：
Figure FDA0003253994060000037
Figure FDA0003253994060000038
Figure FDA0003253994060000041
Figure FDA0003253994060000042
Figure FDA0003253994060000043
式中，UgN、LeqN、PEN分别表示并网点额定电压，理想变压器额定电感、电力电子设备馈入的额定有功；下标k表示第k个单VSC-HVDC馈入子系统，下标r表示第r个单光伏场馈入子系统或单直驱风场馈入子系统，
Figure FDA0003253994060000044
Figure FDA0003253994060000045
均表示子系统系数。
5.根据权利要求4所述的一种多类型电力电子设备馈入系统动态稳定判断方法，其特征在于，所述子系统系数
Figure FDA0003253994060000046
Figure FDA0003253994060000047
的表达式为：
Figure FDA0003253994060000048
Figure FDA0003253994060000049
Figure FDA00032539940600000410
Figure FDA0003253994060000051
Figure FDA0003253994060000052
Figure FDA0003253994060000053
式中，
Figure FDA0003253994060000054
分别表示外环功率控制比例积分环节的积分系数、比例系数，内环电流控制比例积分环节的积分系数、比例系数，锁相环形成角速度环节比例积分环节的积分系数、比例系数，Lc、Leq分别表示滤波电感、理想变压器电感，CDC表示直流侧电容；下标k代表第k个单VSC-HVDC馈入子系统，下标r表示第r个单光伏场馈入子系统或单直驱风场馈入子系统，Ts表示采样周期，ω0是基准角频率。
6.根据权利要求1所述的一种多类型电力电子设备馈入系统动态稳定判断方法，其特征在于，所述改进动态短路比模型IMDSCR的建立方法如下：
S1、建立子系统动态模型；
S2、根据子系统动态模型，建立不含控制器的连续域主电路动态方程，离散化后得到主电路离散状态模型；
S3、建立控制器的离散状态方程；
S4、根据步骤S2和S3的两个离散状态模型建立子系统的离散时间状态模型；
S5、根据所述离散时间状态模型建立李雅普诺夫函数；
S5、基于压缩映射定理，根据所述李雅普诺夫函数建立改进动态短路比模型IMDSCR。
7.根据权利要求6所述的一种多类型电力电子设备馈入系统动态稳定判断方法，其特征在于，所述单VSC-HVDC馈入子系统在模型是基于并网点WAMS设备实时数据建立的。
8.根据权利要求6所述的一种多类型电力电子设备馈入系统动态稳定判断方法，其特征在于，所述李雅普诺夫函数的建立步骤如下：
SA1、通过数学归纳法，求解所述离散时间状态模型在M时刻的解；
SA2、根据步骤SA1中的解，求得在M时刻的通用解值；
SA3、对于步骤SA2中得到的通用解值，选取1-范数作为测度，结合增益函数，得到李雅普诺夫函数。
9.根据权利要求8所述的一种多类型电力电子设备馈入系统动态稳定判断方法，其特征在于，所述增益函数表达式为：
Figure FDA0003253994060000061
式中，
Figure FDA0003253994060000062
表示增益系数。
10.根据权利要求6所述的一种多类型电力电子设备馈入系统动态稳定判断方法，其特征在于，所述离散时间状态模型表达式为：
Figure FDA0003253994060000071
下标k代表第k个单VSC-HVDC馈入子系统，下标r表示第r个单光伏场馈入子系统或单直驱风场馈入子系统，
Figure FDA0003253994060000072
Figure FDA0003253994060000073
均表示子系统系数，uk,r(w)表示输入系数，w表示采样时刻，
Figure FDA0003253994060000074
表示不同子系统的状态变量。
CN202111054374.8A 2021-09-09 2021-09-09 一种多类型电力电子设备馈入系统动态稳定判断方法 Active CN113949104B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111054374.8A CN113949104B (zh) 2021-09-09 2021-09-09 一种多类型电力电子设备馈入系统动态稳定判断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111054374.8A CN113949104B (zh) 2021-09-09 2021-09-09 一种多类型电力电子设备馈入系统动态稳定判断方法

Publications (2)

Publication Number Publication Date
CN113949104A true CN113949104A (zh) 2022-01-18
CN113949104B CN113949104B (zh) 2023-09-26

Family

ID=79328314

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111054374.8A Active CN113949104B (zh) 2021-09-09 2021-09-09 一种多类型电力电子设备馈入系统动态稳定判断方法

Country Status (1)

Country Link
CN (1) CN113949104B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060282239A1 (en) * 2005-06-08 2006-12-14 Chang Gung University Method of setting-up steady state model of VSC-based multi-terminal HVDC transmission system
CN109245146A (zh) * 2018-09-14 2019-01-18 上海电力学院 基于多馈入广义运行短路比的电力系统稳定性判断方法
CN110808603A (zh) * 2019-11-18 2020-02-18 上海电力大学 适用于多台双馈式风机接入受端电网电压稳定评估方法
CN112636324A (zh) * 2020-11-09 2021-04-09 上海电力大学 多电力电子设备馈入复杂弱电网的稳定性判断方法和系统
EP3852213A1 (de) * 2020-01-14 2021-07-21 Innogy SE Steuerungseinheit für eine energieeinspeiseeinheit, energieeinspeiseeinheit, energieeinspeisesystem und netzinfrastruktur

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060282239A1 (en) * 2005-06-08 2006-12-14 Chang Gung University Method of setting-up steady state model of VSC-based multi-terminal HVDC transmission system
CN109245146A (zh) * 2018-09-14 2019-01-18 上海电力学院 基于多馈入广义运行短路比的电力系统稳定性判断方法
CN110808603A (zh) * 2019-11-18 2020-02-18 上海电力大学 适用于多台双馈式风机接入受端电网电压稳定评估方法
EP3852213A1 (de) * 2020-01-14 2021-07-21 Innogy SE Steuerungseinheit für eine energieeinspeiseeinheit, energieeinspeiseeinheit, energieeinspeisesystem und netzinfrastruktur
CN112636324A (zh) * 2020-11-09 2021-04-09 上海电力大学 多电力电子设备馈入复杂弱电网的稳定性判断方法和系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘炜;赵成勇;郭春义;陆翌;: "混合双馈入直流系统中LCC-HVDC对VSC-HVDC稳态运行区域的影响", 中国电机工程学报, no. 13 *
辛焕海;董炜;袁小明;甘德强;王康;谢欢;: "电力电子多馈入电力系统的广义短路比", 中国电机工程学报, no. 22 *

Also Published As

Publication number Publication date
CN113949104B (zh) 2023-09-26

Similar Documents

Publication Publication Date Title
CN109217371B (zh) 考虑锁相环影响的电压源型变流器并网系统稳定性分析方法、装置及系统
CN107294527A (zh) 同步旋转坐标系锁相环及其测试方法、装置
CN105978373B (zh) 实现微电网稳定的三相逆变器逆推滑模控制方法及系统
CN113285481B (zh) 并网变流器电感参数在线估计方法、预测控制方法及系统
CN110011295A (zh) 一种分布式直流电网的稳定性分析方法及系统
CN111130135B (zh) 一种适用于高比例新能源接入下的电力系统惯量计算方法
CN111327075B (zh) 两级式光伏并网系统低电压穿越动态评估方法及系统
CN112310991A (zh) 基于扰动观测的储能变换器改进有限集模型预测控制方法
CN112907075B (zh) 一种电力系统综合负荷模型参数辨识方法
CN111737919A (zh) 一种适用于次同步振荡分析的直驱式风电场分群方法
CN113949104B (zh) 一种多类型电力电子设备馈入系统动态稳定判断方法
CN113644680B (zh) 一种交直流混联配电系统不确定扰动下的可达性分析方法
CN115021256B (zh) 一种大规模交直流输电系统电磁暂态模型自动生成方法
CN114465280A (zh) 一种新能源并网系统动态等效建模方法
CN110504693B (zh) 一种基于负载参数测量扰动下的电力弹簧优化控制方法
CN113224788A (zh) 含有风电接入的交直流混联电网静态电压稳定评估方法
Gao et al. Improved extended kalman filter based dynamic equivalent method of DFIG wind farm cluster
CN113869668B (zh) 一种多电力电子设备馈入系统动态稳定判断方法
CN113675888B (zh) 基于精确离散化的变流器级联预测控制方法及系统
CN114938039B (zh) 单相并网逆变器无功调节方法、系统和计算机存储介质
CN108680792A (zh) 一种含电力电子装置新能源并网系统动态导纳测量方法
CN113452068B (zh) 一种连接风电场的vsc整流站多步模型预测控制方法
CN111697582B (zh) 一种交流微电网系统
CN116299041A (zh) 一种双馈风电机组短路电流的获取方法及系统
CN117220334A (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