CN107069811A - 基于同步参考坐标系的阻抗网络建模与稳定性分析方法 - Google Patents
基于同步参考坐标系的阻抗网络建模与稳定性分析方法 Download PDFInfo
- Publication number
- CN107069811A CN107069811A CN201710237346.7A CN201710237346A CN107069811A CN 107069811 A CN107069811 A CN 107069811A CN 201710237346 A CN201710237346 A CN 201710237346A CN 107069811 A CN107069811 A CN 107069811A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- impedance
- mtd
- mtr
- 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 56
- 239000011159 matrix material Substances 0.000 claims abstract description 113
- 238000004458 analytical method Methods 0.000 claims abstract description 27
- 230000003534 oscillatory effect Effects 0.000 claims abstract description 22
- 238000006243 chemical reaction Methods 0.000 claims abstract description 17
- 230000008901 benefit Effects 0.000 claims abstract description 11
- 238000013461 design Methods 0.000 claims abstract description 8
- 230000002776 aggregation Effects 0.000 claims abstract description 7
- 238000004220 aggregation Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 230000005540 biological transmission Effects 0.000 claims description 30
- 238000011217 control strategy Methods 0.000 claims description 26
- 238000010248 power generation Methods 0.000 claims description 19
- 230000010355 oscillation Effects 0.000 claims description 17
- 230000005611 electricity Effects 0.000 claims description 16
- 238000004146 energy storage Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 11
- 230000001360 synchronised effect Effects 0.000 claims description 9
- 238000005286 illumination Methods 0.000 claims description 6
- 238000006116 polymerization reaction Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 4
- 230000021615 conjugation Effects 0.000 claims description 3
- 210000004027 cell Anatomy 0.000 description 6
- 230000007246 mechanism Effects 0.000 description 5
- 230000014509 gene expression Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000035699 permeability Effects 0.000 description 3
- 238000005094 computer simulation Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 206010008190 Cerebrovascular accident Diseases 0.000 description 1
- 206010068052 Mosaicism Diseases 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000003765 sex chromosome Anatomy 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- H02J3/383—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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/18—Arrangements for adjusting, eliminating or compensating reactive power in networks
-
- 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/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- H02J3/386—
-
- 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/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
- H02J2003/365—Reducing harmonics or oscillations in HVDC
-
- 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]
-
- 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
-
- 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/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
-
- 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
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Abstract
本发明提出一种基于同步参考坐标系的阻抗网络建模与稳定性分析方法,包括以下步骤:S1:采集目标系统的机网参数;S2:根据机网参数进行潮流计算,以得到用于构建系统中各设备阻抗矩阵模型的状态变量初值;S3:根据状态变量初值建立系统中各设备的阻抗矩阵模型;S4:将系统中各设备的阻抗矩阵模型转换到全网统一的坐标系中;S5:构建整体系统的阻抗网络模型;S6:将阻抗网络模型聚合为组合阻抗矩阵;S7:根据组合阻抗矩阵进行振荡稳定性分析。本发明能够实现多变流器‑大电网的阻抗网络建模,具有适用范围广,物理含义明确,简洁高效的优点。
Description
技术领域
本发明涉及电力系统稳定分析技术领域,特别涉及一种基于同步参考坐标系的阻抗网络建模与稳定性分析方法。
背景技术
近年来,电力电子变流器在现代电力系统中得到了广泛应用,典型的应用领域包括可调速的电机拖动、高压直流输电、风电/光伏等新能源发电系统等。作为一种具有良好控制性能的电力电子接口,变流器使得电气设备与电网之间实现“柔性连接”。然而,电力电子变流器渗透率的迅速升高导致新能源发电基地出现了一些新型的稳定性问题,例如,分布式发电和微电网系统中出现的高频谐波振荡;变流器型恒功率负载因其负电阻特性导致的低频功率振荡;风电场和光伏电站发生的新型次同步振荡(SubsynchronousOscillation,SSO)问题。这种新兴的振荡问题已在国内外的多个系统中发生,比如美国的Ajo风力发电系统、中国河北沽源系统、中国新疆地区的新能源发电基地等等。
学术界和工业界已经对新能源发电基地中出现的新型振荡问题开展了大量的研究。为分析方便,建立了单台大容量变流器-交流并网系统的简化模型,并采用特征值分析法和时域仿真法研究了系统中出现的振荡稳定性问题。近年来,有学者建立了单台大容量变流器-交流并网系统的阻抗模型,并采用Nyquist判据来评估系统稳定性。但所有这些研究都是将多个电力电子变流器建模为一台大容量聚合变流器,同时将实际复杂的电网部分简化为一条等效输电线路模型。很明显,这种简化的单变流器-交流并网系统模型跟实际复杂多变流器系统相比,有较大差距,无法考虑变流器地理位置和参数多样化的影响。采用特征值分析法或时域仿真法开展实际系统稳定性分析时将遇到一些困难,例如,建立实际复杂系统的小信号分析模型相对困难,且复杂系统中往往包含多个模式,需要辨识出关心的振荡模式;建立复杂系统的电磁暂态模型工作量大,变流器的多时间尺度动态需要很小的仿真步长,因此开展仿真分析十分耗时。
阻抗分析法具有物理概念清晰,可扩展性强等优点,可应用于复杂多变流器系统的稳定性分析。传统的阻抗分析法在abc坐标系下建立了系统的正负序阻抗模型,但在建模过程中忽略了交-直流侧的动态耦合、交直轴控制器的结构或参数不对称等问题,难以解决实际变流器dq轴阻抗不对称的难题。在稳定性判据方面,传统的Nyquist判据和范数判据虽然能够定性的判断系统的稳定性,但不能量化的分析振荡的频率和阻尼,在实际应用中受到限制。已有的基于dq旋转坐标系的阻抗建模方法也仅适用于单机无穷大系统的分析,无法扩展应用于包含数万个节点和数千台变流器的实际复杂系统。
发明内容
本发明旨在至少解决上述技术问题之一。
为此,本发明的目的在于提出一种基于同步参考坐标系的阻抗网络建模与稳定性分析方法,该方法能够实现多变流器-大电网的阻抗网络建模,具有适用范围广,物理含义明确,简洁高效的优点。
为了实现上述目的,本发明的实施例提出了一种基于同步参考坐标系的阻抗网络建模与稳定性分析方法,包括以下步骤:S1:采集目标系统的机网参数;S2:根据所述机网参数进行潮流计算,以得到用于构建系统中各设备阻抗矩阵模型的状态变量初值;S3:根据所述状态变量初值建立系统中各设备的阻抗矩阵模型;S4:将系统中各设备的阻抗矩阵模型转换到全网统一的坐标系中;S5:构建整体系统的阻抗网络模型;S6:将所述阻抗网络模型聚合为组合阻抗矩阵;S7:根据所述组合阻抗矩阵进行振荡稳定性分析。
另外,根据本发明上述实施例的基于同步参考坐标系的阻抗网络建模与稳定性分析方法还可以具有如下附加的技术特征:
在一些示例中,所述目标系统包括:发电子系统、输电子系统和配用电子系统,所述机网参数至少包括:发电子系统中:火/水电基地中水电机组或者火电机组的台数,水/汽轮机参数,发电机、励磁系统、调速器、轴系系统和机端变压器的结构和参数,储能基地的场站网络拓扑结构,储能基地变流器的控制策略和参数,可再生能源发电基地中风机和太阳能光伏电池板的类型和台数、变流器的控制策略和控制参数、风速和光照强度;输电子系统中:传统/柔性直流输电系统换流站的控制策略和控制器参数,无功补偿装置的结构和容量,直流输电线路的型号和长度,交流输电系统的拓扑结构、各条线路的型号和参数;配用电子系统中:交直流配电系统的拓扑结构、变流器的控制策略和参数,分布式发电系统中风机和太阳能光伏电池板的类型和台数、变流器的控制策略和控制参数、风速和光照强度,变流器型恒功率负载中变流器的控制策略和参数。
在一些示例中,所述S2,进一步包括:根据系统拓扑和相应的元件参数,构建系统的节点导纳矩阵;选取系统一个预设的运行工况,结合各元件的已知运行条件,求取系统未知的运行状态,得到全网的稳态运行状态,并将所述全网的稳态运行状态作为后续构建系统中各设备阻抗矩阵模型的状态变量初始值,其中,所述全网的稳态运行状态至少包括:火电机组、风电机组、光伏单元、储能单元、柔直、FACTS装置的状态变量初值。
在一些示例中,所述S3,进一步包括:
S31:建立电力设备的阻抗矩阵模型,包括:
基于第i个电力设备的dq旋转坐标系,根据其电路结构、详细的控制策略和锁相环,建立该设备的小信号状态空间模型,具体如下式:
式(1)中:x,y,u分别表示系统的状态变量、输出变量和输入变量,A,B,C,D表示具有相应维度的系数矩阵,Δ表示增量计算;
选取设备出口处的电压udqi和电流idqi作为接口变量,即
式(2)经过拉普拉斯变换可整理为s域内的代数方程的形式:
式(3)中:s表示拉普拉斯算子;
推导得到电力设备出口电压和电流之间的关系,即阻抗矩阵模型:
Δudqi=Zi(s)Δidqi=[C(sI-A)-1B+D]Δidqi (4)
式(4)中:Zi表示第i个电力设备的阻抗矩阵模型,阻抗矩阵的阶数是2力设阶;
阻抗矩阵模型Zi可表示为下式:
式(5)中:阻抗矩阵的每个元素Zkmi(s)=akmi(s)/bkmi(s),k=d,q,m=d,q,且akmi(s)和bkmi(s)都是关于拉普拉斯算子s的多项式;
S32:建立交流电网的阻抗矩阵模型,包括:
将系统中交流线路的阻抗转换到同步旋转坐标系下,并进行推导后得到线路的阻抗矩阵模型可表示为:
式(6)中,ZL是一个2一个阶矩阵,ω0=2阵f,f为系统频率,r∑表示线路总电阻,L∑表示线路总电抗。
在一些示例中,所述S4,进一步包括:
将系统中所有设备的阻抗模型转换到全网统一的坐标系的坐标变换公式为:
ZSRFi(s)=T(s)·Zdqi(s)·T-1(s) (7)
式(7)中,ZSRFi表示在同步旋转坐标系下第i个设备的阻抗矩阵模型,T(s)表示2×2阶的坐标转换矩阵。
在一些示例中,所述S5,进一步包括:根据所述系统中各设备的阻抗矩阵模型,在已知的系统拓扑结构基础上,将各设备采用其阻抗矩阵模型来代替,以将各阻抗矩阵模型联立起来,得到所述整体系统的阻抗网络模型。
在一些示例中,在所述S6中,通过分析实际系统的录波数据以确定振荡路径,根据预设的化简方法,沿着所述振荡路径将所述整体系统的阻抗网络模型进行相应的化简,以实现阻抗网络模型的聚合。
在一些示例中,所述预设的化简方法包括:串并联结变换、星三角网络变换或者阻抗矩阵间的变换。
在一些示例中,所述S6,进一步包括:
S61:通过分析现场的录波数据,确定系统中的振荡路径;
S62:将阻抗网络模型中的储能基地阻抗与其相连的线路阻抗串联、光伏基地阻抗与其相连的线路阻抗串联、风电基地阻抗与其相连的线路阻抗串联、传统直流系统阻抗与带串补的交流线路阻抗与交流线路阻抗相并联、柔性直流输电系统阻抗与第一预设交流系统等效阻抗串联,分布式发电阻抗与其相连的线路阻抗串联、串联后的分布式发电阻抗与变流器型恒功率负载阻抗与交直流配电网等效阻抗并联,即:
式(8)中:“+”表示两个阻抗矩阵的串联,“//”表示两个阻抗矩阵的并联,其中,
当Z1=[Zdd-1 Zdq-1;Zqd-1 Zqq-1],Z2=[Zdd-2 Zdq-2;Zqd-2 Zqq-2]时,具体计算过程为:
其中:表示Z1矩阵的逆运算;
S63:当ZLine4(s)、ZLine5(s)与ZLine6(s)采用三角形接法时,将这三个阻抗进行三角-星变换,具体为:
S64:将阻抗Z∑2(s)与ZY2(s)串联,将阻抗ZTG(s)与ZY3(s)串联,进一步将ZTG(s)与ZY3(s)并联后与阻抗ZY1(s)串联,即:
Z∑6(s)={[Z∑2(s)+ZY2(s)]//[ZTG(s)+ZY3(s)]}+ZY1(s) (12)
将阻抗ZS2(s)与ZLine8(s)串联,并将串联后得到的阻抗与阻抗Z∑4(s)并联,即:
Z∑7(s)=[ZS2(s)+ZLine8(s)]//Z∑4(s) (13)
S65:将阻抗Z∑6(s)、Z∑3(s)、Z∑7(s)串联可得到组合阻抗矩阵,即:
Z∑(s)=Z∑6(s)+Z∑3(s)+Z∑7(s) (14)
其中:Z∑(s)=[Zdd-∑(s)Zdq-∑(s);Zqd-∑(s)Zqq-∑(s)]。
在一些示例中,所述S7,进一步包括:
S71:组合阻抗矩阵行列式的求解,包括:
求取得到的系统组合阻抗矩阵Z∑(s)的行列式,可表示为如下式:
D(s)=Zdd-∑(s)Zqq-∑(s)-Zdq-∑(s)Zqd-∑(s) (15)
S72:基于组合阻抗矩阵行列式的稳定性判别,包括:
选取系统中存在的一对共轭零点λ1,2=αo±jωo,且|αo|<<|ωo|,当ω位于λ1,2的微小邻域内时,D(s)可表示为:
D(jω)=(jω-λ1)(jω-λ2)G(jω) (16)
式(16)中,G(jω)=a+jb,a、b为常数;
将式(16)的实部和虚部相分离,可得如下式:
通过求解Im[D(jω)]=0,辨识出系统的过零点频率ωr,通过分析过零点频率ωr处Re[D(jω)]的正负可判断系统的稳定性,其稳定判据如下:
a:当D(s)虚部曲线从负向正穿越过零点时,如果Re{D(jωr)}>0,振荡模态稳定,反之,振荡模态不稳定;
b:当D(s)虚部曲线从正向负穿越过零点时,如果Re{D(jωr)}>0,振荡模态不稳定,反之,振荡模态稳定。
根据本发明实施例的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,具有如下优点:
1、适用于包含大量电力电子装备的现代电力系统的阻抗网络建模,尤其适用于风、光发电渗透率高的电力系统,在建模过程中克服了变流器dq轴阻抗耦合的影响,能够实现多变流器-大电网的阻抗网络建模,即适用范围广;
2、各设备的阻抗矩阵模型是一个2各设阶的阻抗矩阵,其表征的是电力设备的一种外特性模型,阻抗模型不仅能够通过详细的机理建模得到,也可以通过现场实验辨识方法来确定模型的阶数和参数,即通过获取设备端口的电压和电流信息就能简单辨识;
3、大型交直流电网的聚合可通过阻抗的串并联和星-三角变换实现,克服了传统小信号状态方程建模面临的维数灾难题;
4、将整体系统的阻抗网络模型聚合成组合阻抗矩阵后,根据提出的稳定判据,通过分析组合阻抗矩阵行列式的阻抗-频率特性评估振荡模态的稳定性,具有物理含义明确,简洁高效的优点。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是根据本发明实施例的基于同步参考坐标系的阻抗网络建模与稳定性分析方法的流程图;
图2是根据本发明一个实施例的目标系统的结构示意图;
图3是根据本发明一个实施例的目标系统的阻抗网络模型示意图;
图4是根据本发明一个实施例的目标系统的阻抗网络模型的等效聚合过程示意图;
图5是根据本发明一个实施例的目标系统组合阻抗矩阵行列式的阻抗(实部、虚部)-频率特性曲线示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
以下结合附图描述根据本发明实施例的基于同步参考坐标系的阻抗网络建模与稳定性分析方法。
图1是根据本发明一个实施例的基于同步参考坐标系的阻抗网络建模与稳定性分析方法的流程图。如图1所示,该方法包括以下步骤:
步骤S1:采集目标系统的机网参数。
在本发明的一个实施例中,图2展示了目标系统(即电力系统)的结构示意图。针对该系统,首先应仔细分析系统结构并收集相关系统参数。如图2所示,目标系统包括:发电子系统、输电子系统和配用电子系统。基于此,机网参数(即各子系统应收集的具体参数)至少包括:
发电子系统中:火/水电基地中水电机组或者火电机组的台数,水/汽轮机参数,发电机、励磁系统、调速器、轴系系统和机端变压器的结构和参数,储能基地的场站网络拓扑结构,储能基地变流器的控制策略和参数,可再生能源发电基地中风机和太阳能光伏电池板的类型和台数、变流器的控制策略和控制参数、风速和光照强度。换言之,由于发电系统由大规模火/水电基地、储能基地和可再生能源发电(例如风电、光伏)基地构成。因此对于火/水电基地,应收集的参数包括:水电或者火电机组的台数,水/汽轮机参数,发电机、励磁系统、调速器、轴系系统和机端变压器的结构和参数等;对于储能基地,应收集的参数包括:储能基地的场站网络拓扑结构、变流器的控制策略和参数;对于可再生能源发电基地,应收集的参数包括:风机和太阳能光伏电池板的类型和台数(块数),变流器的控制策略和控制参数,风速和光照强度等。
输电子系统中:传统/柔性直流输电系统换流站的控制策略和控制器参数,无功补偿装置的结构和容量,直流输电线路的型号和长度,交流输电系统的拓扑结构、各条线路的型号和参数。换言之,由于输电系统由传统直流、柔性直流和交流输电系统构成。对于传统/柔性直流输电系统,应该收集的参数包括:换流站的控制策略和控制器参数,无功补偿装置的结构和容量,直流输电线路的型号和长度等;对于交流输电系统,应收集的参数包括:系统拓扑结构、各条线路(包含串补输电线路)的型号和参数等。
配用电子系统中:交直流配电系统的拓扑结构、变流器的控制策略和参数,分布式发电系统中风机和太阳能光伏电池板的类型和台数、变流器的控制策略和控制参数,风速和光照强度,变流器型恒功率负载中变流器的控制策略和参数。换言之,由于配用电系统由交直流配电系统、分布式发电系统和变流器型恒功率负载构成。对于交直流配电系统,需要收集的参数包括:系统拓扑结构、变流器的控制策略和参数等;对于分布式发电系统,需要收集的参数与可再生能源发电基地相似;对于变流器型恒功率负载,需要收集的参数包括变流器的控制策略和参数等。
步骤S2:根据机网参数进行潮流计算,以得到用于构建系统中各设备阻抗矩阵模型的状态变量初值。
具体地,在步骤S2中,针对如图2所示的电力系统,根据上述已收集的目标系统参数,进行潮流计算,其具体过程如下:根据系统拓扑和相应的元件参数,构建系统的节点导纳矩阵;选取系统一个预设的运行工况,结合各元件的已知运行条件,求取系统未知的运行状态,得到全网的稳态运行状态,并将全网的稳态运行状态作为后续构建系统中各设备阻抗矩阵模型的状态变量初始值,其中,全网的稳态运行状态至少包括:火电机组、风电机组、光伏单元、储能单元、柔直、FACTS装置等的状态变量初值。
步骤S3:根据状态变量初值建立系统中各设备的阻抗矩阵模型。
具体地,以如图2所示的电力系统为例,需要对该系统中各设备建立阻抗矩阵模型。常用的阻抗矩阵建模有小信号模型分析法、模型辨识法等,实现方法多样。此处以小信号模型分析法为例,说明系统中各设备阻抗矩阵建模过程。也即,步骤S3进一步包括:
S31:建立电力设备的阻抗矩阵模型。电力设备的阻抗矩阵建模方法包括机理建模法和模型实测辨识法,其中,前者适用于部件为“白箱”的情况,后者适用于“黑箱”或者“灰箱”情形。
以下以机理建模法为例描述电力设备阻抗矩阵模型的建模过程,具体步骤包括:
基于第i个电力设备的dq旋转坐标系,根据其电路结构、详细的控制策略和锁相环等因素,建立该设备的小信号状态空间模型,具体如下式:
式(1)中:x,y,u分别表示系统的状态变量,输出变量和输入变量,A,B,C,D表示具有相应维度的系数矩阵,Δ表示增量计算;
为了得到该电气设备的阻抗模型,可以选取设备出口处的电压udqi和电流idqi作为接口变量,即
式(2)经过拉普拉斯变换可整理为s域内的代数方程的形式:
式(3)中:s表示拉普拉斯算子;
进一步推导得到电力设备出口电压和电流之间的关系,即阻抗矩阵模型:
Δudqi=Zi(s)Δidqi=[C(sI-A)-1B+D]Δidqi (4)
式(4)中:Zi表示第i个电力设备的阻抗矩阵模型,阻抗矩阵的阶数是2×2阶;
阻抗矩阵模型Zi可表示为下式:
式(5)中:阻抗矩阵的每个元素Zkmi(s)=akmi(s)/bkmi(s),k=d,q,m=d,q,且akmi(s)和bkmi(s)都是关于拉普拉斯算子s的多项式。
S32:建立交流电网的阻抗矩阵模型,包括:
将系统中交流线路的阻抗转换到同步旋转坐标系下,并进行推导后得到线路的阻抗矩阵模型可表示为:
式(6)中,ZL是一个2×2阶矩阵,ω0=2πf,f为系统频率,r∑表示线路总电阻,L∑表示线路总电抗。
除上述机理建模法外,电力设备的阻抗矩阵模型Zi(s)和交流电网的阻抗矩阵模型ZL(s)也可以通过实验辨识的方法得到,此处不再赘述。
需要说明的是,系统中的各种设备均可通过以上类似的方法得到,具体包括各种类型的风机、光伏发电模块、储能单元、传统/柔性直流输电系统、汽轮机组、各类型输电线路(包括串补输电线路)和变流器型负载等。但需要注意的是,系统中各设备的阻抗模型的建立方法和表示形式均多样,不仅限于上述一种。
步骤S4:将系统中各设备的阻抗矩阵模型转换到全网统一的坐标系中。
具体地,系统中各个设备的阻抗矩阵模型均是基于各自的dq旋转坐标系建立的,为了实现电力设备阻抗矩阵模型的互联,需要将所有设备的阻抗矩阵模型转换到全网统一的同步旋转坐标系(Synchronous Reference Frame,SRF)中。
基于此,在本发明的一个实施例中,将系统中所有设备的阻抗模型转换到全网统一的坐标系的坐标变换公式为:
ZSRFi(s)=T(s)·Zdqi(s)·T-1(s) (7)
式(7)中,ZSRFi表示在同步旋转坐标系下第i个设备的阻抗矩阵模型,T(s)表示2×2阶的坐标转换矩阵。
步骤S5:构建整体系统的阻抗网络模型。
具体地说,将系统中所有设备的阻抗矩阵模型转换到全网统一的同步旋转坐标系后,需要对整体系统构建其阻抗网络模型,具体实现方法如下:在步骤S5中,以图2所示的大型交直流电力系统为例,根据上述建立的系统中各设备的阻抗矩阵模型,在已知的系统拓扑结构基础上,将各设备采用其阻抗矩阵模型来代替,以将各阻抗矩阵模型联立起来,得到整体系统的阻抗网络模型,例如图3所示。
步骤S6:将阻抗网络模型聚合为组合阻抗矩阵。
具体地,在步骤S6中,通过分析实际系统的录波数据可以确定振荡路径,根据预设的化简方法,沿着振荡路径可将整体系统的阻抗网络模型进行相应的化简,以实现阻抗网络模型的聚合。其中,预设的化简方法包括:串并联结变换、星三角网络变换或者阻抗矩阵间的变换。进而,经过一系列的化简计算,最后得到整个网络的组合阻抗矩阵Z(s),一般可由传递函数形式表示,其复杂程度与系统的复杂程度有关。
以下结合具体示例对步骤S6的具体过进行详细说明。针对图3的大型交直流电力系统阻抗网络模型,经过图4所示化简过程后得到整个系统的组合阻抗矩阵Z(s),其具体聚合过程如下:
S61:振荡路径的辨识。通过分析现场的录波数据,确定系统中的振荡路径。此例中假设振荡路径如图4所示。
S62:将阻抗网络模型中的储能基地阻抗与其相连的线路阻抗串联、光伏基地阻抗与其相连的线路阻抗串联、风电基地阻抗与其相连的线路阻抗串联、传统直流系统阻抗与带串补的交流线路阻抗与交流线路阻抗相并联、柔性直流输电系统阻抗与第一预设交流系统(即图2中所示的交流系统1)等效阻抗串联,分布式发电阻抗与其相连的线路阻抗串联、串联后的分布式发电阻抗与变流器型恒功率负载阻抗与交直流配电网等效阻抗并联,即:
式(8)中:“+”表示两个阻抗矩阵的串联,“//”表示两个阻抗矩阵的并联,其中,在具体示例中,
当Z1=[Zdd-1 Zdq-1;Zqd-1 Zqq-1],Z2=[Zdd-2 Zdq-2;Zqd-2 Zqq-2]时,具体计算过程为:
其中:表示Z1矩阵的逆运算。
S63:当ZLine4(s)、ZLine5(s)与ZLine6(s)采用三角形接法时,将这三个阻抗进行三角-星变换,具体为:
S64:将阻抗Z∑2(s)与ZY2(s)串联,将阻抗ZTG(s)与ZY3(s)串联,进一步将ZTG(s)与ZY3(s)并联后与阻抗ZY1(s)串联,即:
Z∑6(s)={[Z∑2(s)+ZY2(s)]//[ZTG(s)+ZY3(s)]}+ZY1(s) (12)
将阻抗ZS2(s)与ZLine8(s)串联,并将串联后得到的阻抗与阻抗Z∑4(s)并联,即:
Z∑7(s)=[ZS2(s)+ZLine8(s)]//Z∑4(s) (13)。
S65:将阻抗Z∑6(s)、Z∑3(s)、Z∑7(s)串联可得到组合阻抗矩阵,即:
Z∑(s)=Z∑6(s)+Z∑3(s)+Z∑7(s) (14)
其中:Z∑(s)=[Zdd-∑(s)Zdq-∑(s);Zqd-∑(s)Zqq-∑(s)]。
步骤S7:根据组合阻抗矩阵进行振荡稳定性分析。
具体地说,在得到上述整个阻抗网络模型的组合阻抗矩阵后,可通过分析组合阻抗矩阵的频率特性实现系统的振荡稳定性分析。本发明的实施例采用系统组合阻抗矩阵行列式的阻抗-频率特性评估系统的稳定性。基于此,具体地,步骤S7进一步包括:
S71:组合阻抗矩阵行列式的求解,包括:
求取上述得到的系统组合阻抗矩阵Z∑(s)的行列式,可表示为如下式:
D(s)=Zdd-∑(s)Zqq-∑(s)-Zdq-∑(s)Zqd-∑(s) (15)。
S72:基于组合阻抗矩阵行列式的稳定性判别。研究表明Z∑(s)行列式零点等于系统特征值,可通过分析行列式的零点评估振荡模态的稳定性,具体包括:
选取系统中存在的一对共轭零点λ1,2=,o±j,o,且|αo|<<|ωo|,当ω位于λ1,2的微小邻域内时,D(s)可表示为:
D(jω)=(jω-λ1)(jω-λ2)G(jω) (16)
式(16)中,G(jω)=a+jb,a、b为常数;
将式(16)的实部和虚部相分离,可得如下式:
通过求解Im[D(jω)]=0,辨识出系统的过零点频率ωr,通过分析过零点频率ωr处Re[D(jω)]的正负可判断系统的稳定性,其稳定判据如下:
a:当D(s)虚部曲线从负向正穿越过零点时,如果Re{D(jωr)}>0,振荡模态稳定,反之,振荡模态不稳定;
b:当D(s)虚部曲线从正向负穿越过零点时,如果Re{D(jωr)}>0,振荡模态不稳定,反之,振荡模态稳定。
如图5所示为某一典型工况下,整体系统组合阻抗矩阵模型行列式的阻抗-频率特性曲线,由图5可以看出,D(s)虚部为0时的过零点频率为ωr=30.96Hz,且在ωr附近D(s)虚部曲线由负向正穿越过零点。由于ωr处的D(s)实部小于0,因此,在该工况下振荡模态不稳定。
综上,根据本发明实施例的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,分别建立电力系统中各个元件的基于同步参考坐标系的阻抗矩阵模型,根据系统拓扑连接成阻抗网络模型,进而将整个阻抗网络模型聚合为一个组合阻抗矩阵,根据所提出的稳定判据评估系统的振荡稳定性。因此,该方法具有如下优点:
1、适用于包含大量电力电子装备的现代电力系统的阻抗网络建模,尤其适用于风、光发电渗透率高的电力系统,在建模过程中克服了变流器dq轴阻抗耦合的影响,能够实现多变流器-大电网的阻抗网络建模,即适用范围广;
2、各设备的阻抗矩阵模型是一个2各设阶的阻抗矩阵,其表征的是电力设备的一种外特性模型,阻抗模型不仅能够通过详细的机理建模得到,也可以通过现场实验辨识方法来确定模型的阶数和参数,即通过获取设备端口的电压和电流信息就能简单辨识;
3、大型交直流电网的聚合可通过阻抗的串并联和星-三角变换实现,克服了传统小信号状态方程建模面临的维数灾难题;
4、将整体系统的阻抗网络模型聚合成组合阻抗矩阵后,根据提出的稳定判据,通过分析组合阻抗矩阵行列式的阻抗-频率特性评估振荡模态的稳定性,具有物理含义明确,简洁高效的优点。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同限定。
Claims (10)
1.一种基于同步参考坐标系的阻抗网络建模与稳定性分析方法,其特征在于,包括以下步骤:
S1:采集目标系统的机网参数;
S2:根据所述机网参数进行潮流计算,以得到用于构建系统中各设备阻抗矩阵模型的状态变量初值;
S3:根据所述状态变量初值建立系统中各设备的阻抗矩阵模型;
S4:将系统中各设备的阻抗矩阵模型转换到全网统一的坐标系中;
S5:构建整体系统的阻抗网络模型;
S6:将所述阻抗网络模型聚合为组合阻抗矩阵;
S7:根据所述组合阻抗矩阵进行振荡稳定性分析。
2.根据权利要求1所述的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,其特征在于,所述目标系统包括:发电子系统、输电子系统和配用电子系统,所述机网参数至少包括:
发电子系统中:火/水电基地中水电机组或者火电机组的台数,水/汽轮机参数,发电机、励磁系统、调速器、轴系系统和机端变压器的结构和参数,储能基地的场站网络拓扑结构,储能基地变流器的控制策略和参数,可再生能源发电基地中风机和太阳能光伏电池板的类型和台数、变流器的控制策略和控制参数、风速和光照强度;
输电子系统中:传统/柔性直流输电系统换流站的控制策略和控制器参数,无功补偿装置的结构和容量,直流输电线路的型号和长度,交流输电系统的拓扑结构、各条线路的型号和参数;
配用电子系统中:交直流配电系统的拓扑结构、变流器的控制策略和参数,分布式发电系统中风机和太阳能光伏电池板的类型和台数、变流器的控制策略和控制参数、风速和光照强度,变流器型恒功率负载中变流器的控制策略和参数。
3.根据权利要求1所述的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,其特征在于,所述S2,进一步包括:
根据系统拓扑和相应的元件参数,构建系统的节点导纳矩阵;选取系统一个预设的运行工况,结合各元件的已知运行条件,求取系统未知的运行状态,得到全网的稳态运行状态,并将所述全网的稳态运行状态作为后续构建系统中各设备阻抗矩阵模型的状态变量初始值,其中,所述全网的稳态运行状态至少包括:火电机组、风电机组、光伏单元、储能单元、柔直、FACTS装置的状态变量初值。
4.根据权利要求1所述的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,其特征在于,所述S3,进一步包括:
S31:建立电力设备的阻抗矩阵模型,包括:
基于第i个电力设备的dq旋转坐标系,根据其电路结构、详细的控制策略和锁相环,建立该设备的小信号状态空间模型,具体如下式:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>+</mo>
<mi>B</mi>
<mi>&Delta;</mi>
<mi>u</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>y</mi>
<mo>=</mo>
<mi>C</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>+</mo>
<mi>D</mi>
<mi>&Delta;</mi>
<mi>u</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)中:x,y,u分别表示系统的状态变量、输出变量和输入变量,A,B,C,D表示具有相应维度的系数矩阵,Δ表示增量计算;
选取设备出口处的电压udqi和电流idqi作为接口变量,即
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>B&Delta;i</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;u</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>C</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>D&Delta;i</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式(2)经过拉普拉斯变换可整理为s域内的代数方程的形式:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>=</mo>
<mi>A</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>B&Delta;i</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;u</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>C</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>D&Delta;i</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式(3)中:s表示拉普拉斯算子;
推导得到电力设备出口电压和电流之间的关系,即阻抗矩阵模型:
Δudqi=Zi(s)Δidqi=[C(sI-A)-1B+D]Δidqi (4)
式(4)中:Zi表示第i个电力设备的阻抗矩阵模型,阻抗矩阵的阶数是2力设阶;
阻抗矩阵模型Zi可表示为下式:
<mrow>
<msub>
<mi>Z</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>d</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>q</mi>
<mi>d</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>q</mi>
<mi>q</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式(5)中:阻抗矩阵的每个元素Zkmi(s)=akmi(s)/bkmi(s),k=d,q,m=d,q,且akmi(s)和bkmi(s)都是关于拉普拉斯算子s的多项式;
S32:建立交流电网的阻抗矩阵模型,包括:
将系统中交流线路的阻抗转换到同步旋转坐标系下,并进行推导后得到线路的阻抗矩阵模型可表示为:
<mrow>
<msub>
<mi>Z</mi>
<mi>L</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mi>&Sigma;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>sL</mi>
<mi>&Sigma;</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>L</mi>
<mi>&Sigma;</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<msub>
<mi>L</mi>
<mi>&Sigma;</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>r</mi>
<mi>&Sigma;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>sL</mi>
<mi>&Sigma;</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式(6)中,ZL是一个2一个阶矩阵,ω0=2阵f,f为系统频率,r∑表示线路总电阻,L∑表示线路总电抗。
5.根据权利要求1所述的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,其特征在于,所述S4,进一步包括:
将系统中所有设备的阻抗模型转换到全网统一的坐标系的坐标变换公式为:
ZSRFi(s)=T(s)·Zdqi(s)·T-1(s) (7)
式(7)中,ZSRFi表示在同步旋转坐标系下第i个设备的阻抗矩阵模型,T(s)表示2示的阶的坐标转换矩阵。
6.根据权利要求1所述的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,其特征在于,所述S5,进一步包括:
根据所述系统中各设备的阻抗矩阵模型,在已知的系统拓扑结构基础上,将各设备采用其阻抗矩阵模型来代替,以将各阻抗矩阵模型联立起来,得到所述整体系统的阻抗网络模型。
7.根据权利要求1所述的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,其特征在于,在所述S6中,通过分析实际系统的录波数据以确定振荡路径,根据预设的化简方法,沿着所述振荡路径将所述整体系统的阻抗网络模型进行相应的化简,以实现阻抗网络模型的聚合。
8.根据权利要求7所述的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,其特征在于,所述预设的化简方法包括:串并联结变换、星三角网络变换或者阻抗矩阵间的变换。
9.根据权利要求7所述的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,其特征在于,所述S6,进一步包括:
S61:通过分析现场的录波数据,确定系统中的振荡路径;
S62:将阻抗网络模型中的储能基地阻抗与其相连的线路阻抗串联、光伏基地阻抗与其相连的线路阻抗串联、风电基地阻抗与其相连的线路阻抗串联、传统直流系统阻抗与带串补的交流线路阻抗与交流线路阻抗相并联、柔性直流输电系统阻抗与第一预设交流系统等效阻抗串联,分布式发电阻抗与其相连的线路阻抗串联、串联后的分布式发电阻抗与变流器型恒功率负载阻抗与交直流配电网等效阻抗并联,即:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>S</mi>
<mi>t</mi>
<mi>o</mi>
<mi>r</mi>
<mi>a</mi>
<mi>g</mi>
<mi>e</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>P</mi>
<mi>V</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>/</mo>
<mo>&lsqb;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>W</mi>
<mi>T</mi>
<mi>G</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>C</mi>
<mi>C</mi>
<mo>-</mo>
<mi>H</mi>
<mi>V</mi>
<mi>D</mi>
<mi>C</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mo>/</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mo>-</mo>
<mi>S</mi>
<mi>C</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mo>/</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>7</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>4</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>9</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>/</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>C</mi>
<mi>P</mi>
<mi>L</mi>
</mrow>
</msub>
<mo>/</mo>
<mo>/</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>M</mi>
<mi>G</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>5</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>V</mi>
<mi>S</mi>
<mi>C</mi>
<mo>-</mo>
<mi>H</mi>
<mi>V</mi>
<mi>D</mi>
<mi>C</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>S</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式(8)中:“+”表示两个阻抗矩阵的串联,“//”表示两个阻抗矩阵的并联,其中,
当Z1=[Zdd-1 Zdq-1;Zqd-1 Zqq-1],Z2=[Zdd-2 Zdq-2;Zqd-2 Zqq-2]时,具体计算过程为:
<mrow>
<msub>
<mi>Z</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>d</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>q</mi>
<mi>d</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>q</mi>
<mi>q</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>d</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>q</mi>
<mi>d</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>Z</mi>
<mrow>
<mi>q</mi>
<mi>q</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>d</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>d</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>q</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>q</mi>
<mi>d</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>q</mi>
<mi>d</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>q</mi>
<mi>q</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>q</mi>
<mi>q</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>Z</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<mo>/</mo>
<msub>
<mi>Z</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>Z</mi>
<mn>1</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>Z</mi>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:表示Z1矩阵的逆运算;
S63:当ZLine4(s)、ZLine5(s)与ZLine6(s)采用三角形接法时,将这三个阻抗进行三角-星变换,具体为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>Y</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>4</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>5</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>6</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>4</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>6</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>Y</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>4</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>5</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>6</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>4</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>5</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>Y</mi>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>4</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>5</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>l</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>6</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>5</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mi>i</mi>
<mi>n</mi>
<mi>e</mi>
<mn>6</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
S64:将阻抗Z∑2(s)与ZY2(s)串联,将阻抗ZTG(s)与ZY3(s)串联,进一步将ZTG(s)与ZY3(s)并联后与阻抗ZY1(s)串联,即:
Z∑6(s)={[Z∑2(s)+ZY2(s)]//[ZTG(s)+ZY3(s)]}+ZY1(s) (12)
将阻抗ZS2(s)与ZLine8(s)串联,并将串联后得到的阻抗与阻抗Z∑4(s)并联,即:
Z∑7(s)=[ZS2(s)+ZLine8(s)]//Z∑4(s) (13)
S65:将阻抗Z∑6(s)、Z∑3(s)、Z∑7(s)串联可得到组合阻抗矩阵,即:
Z∑(s)=Z∑6(s)+Z∑3(s)+Z∑7(s) (14)
其中:Z∑(s)=[Zdd-∑(s) Zdq-∑(s);Zqd-∑(s) Zqq-∑(s)]。
10.根据权利要求9所述的基于同步参考坐标系的阻抗网络建模与稳定性分析方法,其特征在于,所述S7,进一步包括:
S71:组合阻抗矩阵行列式的求解,包括:
求取得到的系统组合阻抗矩阵Z∑(s)的行列式,可表示为如下式:
D(s)=Zdd-∑(s)Zqq-∑(s)-Zdq-∑(s)Zqd-∑(s) (15)
S72:基于组合阻抗矩阵行列式的稳定性判别,包括:
选取系统中存在的一对共轭零点λ1,2=αo±jωo,且|αo|<<|ωo|,当ω位于λ1,2的微小邻域内时,D(s)可表示为:
D(jω)=(jω-λ1)(jω-λ2)G(jω) (16)
式(16)中,G(jω)=a+jb,a、b为常数;
将式(16)的实部和虚部相分离,可得如下式:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>Re</mi>
<mo>&lsqb;</mo>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>o</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&omega;</mi>
<mi>o</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>b&alpha;</mi>
<mi>o</mi>
</msub>
<mi>&omega;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Im</mi>
<mo>&lsqb;</mo>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>=</mo>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>o</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>&omega;</mi>
<mi>o</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>a&alpha;</mi>
<mi>o</mi>
</msub>
<mi>&omega;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
通过求解Im[D(jω)]=0,辨识出系统的过零点频率ωr,通过分析过零点频率ωr处Re[D(jω)]的正负可判断系统的稳定性,其稳定判据如下:
a:当D(s)虚部曲线从负向正穿越过零点时,如果Re{D(jωr)}>0,振荡模态稳定,反之,振荡模态不稳定;
b:当D(s)虚部曲线从正向负穿越过零点时,如果Re{D(jωr)}>0,振荡模态不稳定,反之,振荡模态稳定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710237346.7A CN107069811B (zh) | 2017-04-12 | 2017-04-12 | 基于同步参考坐标系的阻抗网络建模与稳定性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710237346.7A CN107069811B (zh) | 2017-04-12 | 2017-04-12 | 基于同步参考坐标系的阻抗网络建模与稳定性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107069811A true CN107069811A (zh) | 2017-08-18 |
CN107069811B CN107069811B (zh) | 2019-07-26 |
Family
ID=59602276
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710237346.7A Active CN107069811B (zh) | 2017-04-12 | 2017-04-12 | 基于同步参考坐标系的阻抗网络建模与稳定性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107069811B (zh) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107688722A (zh) * | 2017-09-27 | 2018-02-13 | 清华大学 | 电压源型变流器的导纳模型和阻抗模型获取方法及装置 |
CN108154315A (zh) * | 2018-01-19 | 2018-06-12 | 国电南瑞科技股份有限公司 | 一种考虑锁相环影响的并网变流器次同步振荡风险分析方法 |
CN108281986A (zh) * | 2018-03-26 | 2018-07-13 | 湖南大学 | 电压控制型虚拟同步发电机的阻抗建模与稳定性分析方法 |
CN108418253A (zh) * | 2018-03-26 | 2018-08-17 | 湖南大学 | 电流控制型虚拟同步发电机的阻抗建模与稳定性分析方法 |
CN108462192A (zh) * | 2018-04-28 | 2018-08-28 | 清华大学 | 电力系统宽频带振荡的模式分析方法 |
CN108667048A (zh) * | 2018-05-31 | 2018-10-16 | 清华大学 | 新能源并网系统振荡稳定性的频域判稳方法及装置 |
CN108667008A (zh) * | 2018-04-20 | 2018-10-16 | 上海交通大学 | 一种用于电力系统次同步振荡分析的通用小信号建模方法 |
CN108808699A (zh) * | 2018-07-10 | 2018-11-13 | 华北电力大学(保定) | 一种适用于双向储能设备的双象限频率特性分析方法 |
CN109494806A (zh) * | 2018-12-28 | 2019-03-19 | 四川科陆新能电气有限公司 | 一种与储能变流器并联均流的柴油发电机系统阻抗检测方法 |
CN109672217A (zh) * | 2018-12-13 | 2019-04-23 | 华北电力大学 | 风电机组并网系统次同步振荡稳定性量化分析方法 |
CN109742748A (zh) * | 2019-02-28 | 2019-05-10 | 东南大学 | 一种适用于多电压等级直流配电系统的稳定判据方法 |
WO2019091130A1 (zh) * | 2017-11-09 | 2019-05-16 | 清华大学 | 次/超同步谐振/振荡的系统级保护系统及方法 |
CN110263377A (zh) * | 2019-05-21 | 2019-09-20 | 上海交通大学 | 一种基于频域映射的风电场单机等值聚合建模方法 |
CN110429611A (zh) * | 2019-07-19 | 2019-11-08 | 国网新疆电力有限公司电力科学研究院 | 一种静止无功补偿器序阻抗建模及控制参数调整方法 |
CN111082441A (zh) * | 2020-01-20 | 2020-04-28 | 清华大学 | 一种考虑限幅非线性的变流器大信号阻抗计算方法 |
CN111146804A (zh) * | 2020-01-02 | 2020-05-12 | 南方电网科学研究院有限责任公司 | 风电-柔性直流输电系统振荡稳定性判别方法和装置 |
CN111327061A (zh) * | 2020-02-12 | 2020-06-23 | 南方电网科学研究院有限责任公司 | 基于视在阻抗的振荡稳定性判别方法和装置 |
CN112014642A (zh) * | 2020-08-03 | 2020-12-01 | 清华大学 | 静止坐标系下电网频率耦合阻抗模型聚合计算方法和装置 |
CN112214887A (zh) * | 2020-09-28 | 2021-01-12 | 国网河南省电力公司电力科学研究院 | 一种蓄电池模型的改进方法 |
CN112217191A (zh) * | 2019-07-11 | 2021-01-12 | 华北电力大学(保定) | 一种基于节点阻抗矩阵的直流配电网稳定性分析方法 |
CN112380797A (zh) * | 2020-11-05 | 2021-02-19 | 中国第一汽车股份有限公司 | 电机建模方法、装置、设备和介质 |
CN113224782A (zh) * | 2021-06-03 | 2021-08-06 | 南方电网科学研究院有限责任公司 | 一种高压直流输电系统电磁暂态解析方法 |
CN113725910A (zh) * | 2021-10-08 | 2021-11-30 | 南通大学 | 一种风电场并网系统的稳定性分析与量化评估方法 |
CN114362144A (zh) * | 2021-12-07 | 2022-04-15 | 南方电网科学研究院有限责任公司 | 一种新能源多馈入系统稳定性的评估方法、装置及设备 |
WO2022127172A1 (zh) * | 2020-12-18 | 2022-06-23 | 国网江苏省电力有限公司经济技术研究院 | 一种评估海上风电场柔性直流输电系统谐振稳定性的方法 |
CN115688344A (zh) * | 2022-12-29 | 2023-02-03 | 国网江西省电力有限公司电力科学研究院 | 多变流器并网系统的宽频带振荡分析方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021682A (zh) * | 2016-05-13 | 2016-10-12 | 清华大学 | 基于阻抗网络模型的次同步振荡分析方法和装置 |
CN106372339A (zh) * | 2016-09-05 | 2017-02-01 | 清华大学 | 电力电子化电力系统的多速率仿真方法及装置 |
-
2017
- 2017-04-12 CN CN201710237346.7A patent/CN107069811B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106021682A (zh) * | 2016-05-13 | 2016-10-12 | 清华大学 | 基于阻抗网络模型的次同步振荡分析方法和装置 |
CN106372339A (zh) * | 2016-09-05 | 2017-02-01 | 清华大学 | 电力电子化电力系统的多速率仿真方法及装置 |
Cited By (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107688722B (zh) * | 2017-09-27 | 2020-03-10 | 清华大学 | 电压源型变流器的导纳模型和阻抗模型获取方法及装置 |
CN107688722A (zh) * | 2017-09-27 | 2018-02-13 | 清华大学 | 电压源型变流器的导纳模型和阻抗模型获取方法及装置 |
WO2019091130A1 (zh) * | 2017-11-09 | 2019-05-16 | 清华大学 | 次/超同步谐振/振荡的系统级保护系统及方法 |
CN108154315B (zh) * | 2018-01-19 | 2021-11-05 | 国电南瑞科技股份有限公司 | 一种考虑锁相环影响的并网变流器次同步振荡风险分析方法 |
CN108154315A (zh) * | 2018-01-19 | 2018-06-12 | 国电南瑞科技股份有限公司 | 一种考虑锁相环影响的并网变流器次同步振荡风险分析方法 |
CN108418253A (zh) * | 2018-03-26 | 2018-08-17 | 湖南大学 | 电流控制型虚拟同步发电机的阻抗建模与稳定性分析方法 |
CN108418253B (zh) * | 2018-03-26 | 2021-06-04 | 湖南大学 | 电流控制型虚拟同步发电机的阻抗建模与稳定性分析方法 |
CN108281986B (zh) * | 2018-03-26 | 2021-03-02 | 湖南大学 | 电压控制型虚拟同步发电机的阻抗建模与稳定性分析方法 |
CN108281986A (zh) * | 2018-03-26 | 2018-07-13 | 湖南大学 | 电压控制型虚拟同步发电机的阻抗建模与稳定性分析方法 |
CN108667008A (zh) * | 2018-04-20 | 2018-10-16 | 上海交通大学 | 一种用于电力系统次同步振荡分析的通用小信号建模方法 |
CN108667008B (zh) * | 2018-04-20 | 2020-07-14 | 上海交通大学 | 一种用于电力系统次同步振荡分析的通用小信号建模方法 |
CN108462192A (zh) * | 2018-04-28 | 2018-08-28 | 清华大学 | 电力系统宽频带振荡的模式分析方法 |
CN108667048A (zh) * | 2018-05-31 | 2018-10-16 | 清华大学 | 新能源并网系统振荡稳定性的频域判稳方法及装置 |
CN108667048B (zh) * | 2018-05-31 | 2021-04-27 | 清华大学 | 新能源并网系统振荡稳定性的频域判稳方法及装置 |
CN108808699A (zh) * | 2018-07-10 | 2018-11-13 | 华北电力大学(保定) | 一种适用于双向储能设备的双象限频率特性分析方法 |
CN109672217A (zh) * | 2018-12-13 | 2019-04-23 | 华北电力大学 | 风电机组并网系统次同步振荡稳定性量化分析方法 |
CN109672217B (zh) * | 2018-12-13 | 2022-03-01 | 华北电力大学 | 风电机组并网系统次同步振荡稳定性量化分析方法 |
CN109494806B (zh) * | 2018-12-28 | 2022-03-18 | 四川科陆新能电气有限公司 | 一种与储能变流器并联均流的柴油发电系统阻抗检测方法 |
CN109494806A (zh) * | 2018-12-28 | 2019-03-19 | 四川科陆新能电气有限公司 | 一种与储能变流器并联均流的柴油发电机系统阻抗检测方法 |
CN109742748A (zh) * | 2019-02-28 | 2019-05-10 | 东南大学 | 一种适用于多电压等级直流配电系统的稳定判据方法 |
CN110263377A (zh) * | 2019-05-21 | 2019-09-20 | 上海交通大学 | 一种基于频域映射的风电场单机等值聚合建模方法 |
CN112217191A (zh) * | 2019-07-11 | 2021-01-12 | 华北电力大学(保定) | 一种基于节点阻抗矩阵的直流配电网稳定性分析方法 |
CN112217191B (zh) * | 2019-07-11 | 2023-07-28 | 华北电力大学(保定) | 一种基于节点阻抗矩阵的直流配电网稳定性分析方法 |
CN110429611A (zh) * | 2019-07-19 | 2019-11-08 | 国网新疆电力有限公司电力科学研究院 | 一种静止无功补偿器序阻抗建模及控制参数调整方法 |
CN111146804A (zh) * | 2020-01-02 | 2020-05-12 | 南方电网科学研究院有限责任公司 | 风电-柔性直流输电系统振荡稳定性判别方法和装置 |
CN111082441A (zh) * | 2020-01-20 | 2020-04-28 | 清华大学 | 一种考虑限幅非线性的变流器大信号阻抗计算方法 |
CN111327061B (zh) * | 2020-02-12 | 2021-08-03 | 南方电网科学研究院有限责任公司 | 基于视在阻抗的振荡稳定性判别方法和装置 |
CN111327061A (zh) * | 2020-02-12 | 2020-06-23 | 南方电网科学研究院有限责任公司 | 基于视在阻抗的振荡稳定性判别方法和装置 |
CN112014642A (zh) * | 2020-08-03 | 2020-12-01 | 清华大学 | 静止坐标系下电网频率耦合阻抗模型聚合计算方法和装置 |
CN112214887A (zh) * | 2020-09-28 | 2021-01-12 | 国网河南省电力公司电力科学研究院 | 一种蓄电池模型的改进方法 |
CN112380797A (zh) * | 2020-11-05 | 2021-02-19 | 中国第一汽车股份有限公司 | 电机建模方法、装置、设备和介质 |
CN112380797B (zh) * | 2020-11-05 | 2023-03-14 | 中国第一汽车股份有限公司 | 电机建模方法、装置、设备和介质 |
WO2022127172A1 (zh) * | 2020-12-18 | 2022-06-23 | 国网江苏省电力有限公司经济技术研究院 | 一种评估海上风电场柔性直流输电系统谐振稳定性的方法 |
CN113224782A (zh) * | 2021-06-03 | 2021-08-06 | 南方电网科学研究院有限责任公司 | 一种高压直流输电系统电磁暂态解析方法 |
CN113725910A (zh) * | 2021-10-08 | 2021-11-30 | 南通大学 | 一种风电场并网系统的稳定性分析与量化评估方法 |
CN114362144A (zh) * | 2021-12-07 | 2022-04-15 | 南方电网科学研究院有限责任公司 | 一种新能源多馈入系统稳定性的评估方法、装置及设备 |
CN114362144B (zh) * | 2021-12-07 | 2023-09-26 | 南方电网科学研究院有限责任公司 | 一种新能源多馈入系统稳定性的评估方法、装置及设备 |
CN115688344A (zh) * | 2022-12-29 | 2023-02-03 | 国网江西省电力有限公司电力科学研究院 | 多变流器并网系统的宽频带振荡分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107069811B (zh) | 2019-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107069811B (zh) | 基于同步参考坐标系的阻抗网络建模与稳定性分析方法 | |
Kunjumuhammed et al. | Electrical oscillations in wind farm systems: Analysis and insight based on detailed modeling | |
Kalcon et al. | Small-signal stability analysis of multi-terminal VSC-based DC transmission systems | |
CN103197559B (zh) | 一种改善双馈机组接入后系统小干扰稳定性的方法 | |
Li et al. | Aggregated models and transient performances of a mixed wind farm with different wind turbine generator systems | |
Zeddini et al. | PSO-based MPPT control of wind-driven Self-Excited Induction Generator for pumping system | |
CN105680483A (zh) | 一种虚拟同步机阻尼的配置方法及装置 | |
CN112671006A (zh) | 一种评估海上风电场柔性直流输电系统谐振稳定性的方法 | |
CN107666155A (zh) | 基于Markov模型的多能互补系统随机稳定性分析方法 | |
CN109921421B (zh) | 双馈风电机组输出谐波电流模型建立方法 | |
CN105162099A (zh) | 一种确定分布式发电接入电网不对称短路电流的运算曲面法 | |
CN103809650B (zh) | 一种光伏发电系统的等效建模方法 | |
CN107403047A (zh) | 一种含风机和光伏发电的直流微电网数字物理仿真平台及仿真方法 | |
CN105720573A (zh) | 基于实测数据的风光储电站有功及无功控制系统建模方法 | |
Chao et al. | A comprehensive review on dynamic equivalent modeling of large photovoltaic power plants | |
CN112072692A (zh) | 一种新能源发电场站的阻抗等值方法及装置 | |
Li et al. | Dynamic aggregation of doubly-fed induction generators (DFIGs) for stability analysis of wind power systems | |
Altimania | Modeling of doubly-fed induction generators connected to distribution system based on eMEGASim® real-time digital simulator | |
CN111817359B (zh) | 一种基于等效电路的微电网可求解边界分析的方法 | |
Shen et al. | Modeling and Parameter Identification of the Photovoltaic Inverter based on VSG | |
CN111293713B (zh) | 一种风电场风机变流器仿真模型的建立方法及仿真方法 | |
Li et al. | Structure Preserving Aggregation Method for Doubly-Fed Induction Generators in Wind Power Conversion | |
Ye et al. | Simulation and Analysis of PMSG-based Wind Energy Conversion System using Different Coverter Models | |
Efe et al. | A novel approach to power flow analysis for grid connected micro grid | |
Zhu et al. | Modeling of Large-scale Offshore Wind Farm Based on Vector Modeling Technique |
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 |