CN106372339B - 电力电子化电力系统的多速率仿真方法及装置 - Google Patents
电力电子化电力系统的多速率仿真方法及装置 Download PDFInfo
- Publication number
- CN106372339B CN106372339B CN201610803055.5A CN201610803055A CN106372339B CN 106372339 B CN106372339 B CN 106372339B CN 201610803055 A CN201610803055 A CN 201610803055A CN 106372339 B CN106372339 B CN 106372339B
- Authority
- CN
- China
- Prior art keywords
- simulation
- fast
- slow
- slow system
- varying
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 92
- 238000004088 simulation Methods 0.000 claims abstract description 135
- 230000001052 transient effect Effects 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 230000005540 biological transmission Effects 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 11
- 230000008878 coupling Effects 0.000 claims description 7
- 238000010168 coupling process Methods 0.000 claims description 7
- 238000005859 coupling reaction Methods 0.000 claims description 7
- 238000000638 solvent extraction Methods 0.000 claims description 2
- 238000005312 nonlinear dynamic Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
-
- 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
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种电力电子化电力系统的多速率仿真方法及装置,其中,方法包括:将电力电子化电力系统划分为慢系统和多个快子系统;分别建立等值模型;设定仿真参数;通过对电力系统进行潮流计算得到电力系统的稳态初始值;建立多个快子系统与慢子系统间的接口,并且建立时变戴维南等效电路和时变诺顿等效电路;根据建立的时变戴维南等效电路和时变诺顿等效电路得到接口参数,并通过并行求解每个快子系统和慢子系统的等值模型对应的节点导纳方程,直到仿真结束,进而得到电磁暂态仿真结果。本发明实施例的仿真方法不仅可以保留交流系统或者直流电网内部的非线性动态特性,同时,仿真效率相对于单一步长结果得到明显提高。
Description
技术领域
本发明涉及电力电子技术领域,特别涉及一种电力电子化电力系统的多速率仿真方法及装置。
背景技术
随着柔性直流(两电平拓扑、多电平拓扑等)、传统高压直流、新能源、电力电子固态变压器等大量电力电子设备接入系统,网络规模节点数急剧增加并且电力电子频繁变拓扑的特征导致网络节点导纳矩阵需要频繁求解,导致仿真效率急剧下降。目前,包含电力电子设备的电力系统电磁暂态仿真仅能对较小规模进行仿真计算,目前仿真的规模和效率均难以满足电网运行、规划或者行为特性的要求。随着高压直流输电和中低压配单网构成的电力电子单元大规模接入交流系统以后,该问题变得尤为突出。目前,为了提高仿真规模,多采用较粗糙的电力电子单元准稳态模型或者交流系统采用戴维南/诺顿等值,仿真精度不高,数值稳定性问题突出。
国内已投运多个多回柔性直流输电工程,如南澳三端柔性直流输电、南汇柔性直流输电、舟山五端柔性直流输电等工程投运,在建的有厦门柔性直流输电工程、云南鲁西直流背靠背工程。为了研究上述实际工程问题以及提供理论支撑,需要针对大规模交流系统和直流电网进行整体建模与仿真验证。相关技术中,多采用等值方式来提高仿真效率,例如交流电网采用戴维南/诺顿等值,直流电网采用受控源模型。然而,相关技术中均无法保留交流系统或者直流电网内部的非线性动态特性,不但仿真结果可信度较差,有时甚至会出现错误的结果。
发明内容
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
为此,本发明的一个目的在于提出一种电力电子化电力系统的多速率仿真方法,该方法可以提高仿真效率,简单易实现。
本发明的另一个目的在于提出一种电力电子化电力系统的多速率仿真装置。
为达到上述目的,本发明一方面实施例提出了一种电力电子化电力系统的多速率仿真方法,包括以下步骤:将电力电子化电力系统进行网络划分,以划分为慢系统和多个快子系统,其中交流系统位于所述慢系统中,电力电子设备位于不同的快子系统中;分别建立所述慢系统和所述多个快子系统的等值模型;设定仿真参数,其中,所述仿真参数包括仿真总时长、所述慢系统的仿真步长和所述每个快子系统的仿真步长;通过对电力系统进行潮流计算得到所述电力系统的稳态初始值;建立所述多个快子系统与所述慢系统间的接口,并且在所述每个快子系统中建立时变戴维南等效电路,并且在所述慢系统中建立时变诺顿等效电路;根据所述每个快子系统中建立的所述时变戴维南等效电路和所述慢系统中建立的所述时变诺顿等效电路得到接口参数,并通过并行求解所述每个快子系统和慢系统的等值模型对应的节点导纳方程,直到仿真结束,进而得到电力电子化电力系统的电磁暂态仿真结果。
本发明实施例的电力电子化电力系统的多速率仿真方法,可以结合含有电力电子单元的电力系统自身多尺度特性,采用多速率仿真,不仅可以保留交流系统或者直流电网内部的非线性动态特性,同时,仿真效率相对于单一步长结果得到明显提高,不仅可以满足仿真精度和数值稳定性要求,同时极大地提高了仿真效率,利于工程推广实现。
另外,根据本发明上述实施例的电力电子化电力系统的多速率仿真方法还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述设定仿真参数进一步包括:在所述慢系统中,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法将所述慢系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到所述慢系统中选择支路对应数值阻抗;针对所述每个快子系统,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法分别将所述每个快子系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到所述每个快子系统中选择支路对应数值阻抗;定义误差系数为所述慢系统与每个快系统的数值阻抗之比,以根据相对误差最小要求选择速率比。
进一步地,在本发明的一个实施例中,所述根据所述每个快子系统中建立的时变戴维南等效电路得到接口参数,进一步包括:获取所述每个快子系统中时变戴维南等效电路的等效参数;获取接口等值戴维南阻抗矩阵;考虑多端口之间的耦合影响,修正所述每个快子系统在预设时间内的戴维南等效参数;获取所述慢系统在预设时间内的诺顿等效参数。
进一步地,在本发明的一个实施例中,通过分网方法进行网络划分,其中,所述分网方法包括节点分裂法、长输电线解耦法、变压器分网法、支路切割法和戴维南诺顿等值分网法。
进一步地,在本发明的一个实施例中,所述慢系统由发电机、变压器、传输线、负荷组成,所述多个快子系统由新能源模块、柔性直流模块和高压直流模块组成。
为达到上述目的,本发明另一方面实施例提出了一种电力电子化电力系统的多速率仿真装置,包括:分网模块,用于将电力电子化电力系统进行网络划分,以划分为慢系统和多个快子系统,其中交流系统位于所述慢系统中,电力电子设备位于不同的快子系统中;第一建立模块,用于分别建立所述慢系统和所述多个快子系统的等值模型;设定模块,用于设定仿真参数,其中,所述仿真参数包括仿真总时长、所述慢系统的仿真步长和所述每个快子系统的仿真步长;计算模块,用于通过对电力系统进行潮流计算得到所述电力系统的稳态初始值;第二建立模块,用于建立所述多个快子系统与所述慢系统间的接口,并且在所述每个快子系统中建立时变戴维南等效电路,并且在慢系统中建立时变诺顿等效电路;获取模块,用于根据所述每个快子系统中建立的所述时变戴维南等效电路和所述慢系统中建立的所述时变诺顿等效电路得到接口参数,并通过并行求解所述每个快子系统和慢系统的等值模型对应的节点导纳方程,直到仿真结束,进而得到电力电子化电力系统的电磁暂态仿真结果。
本发明实施例的电力电子化电力系统的多速率仿真装置,可以结合含有电力电子单元的电力系统自身多尺度特性,采用多速率仿真,不仅可以保留交流系统或者直流电网内部的非线性动态特性,同时,仿真效率相对于单一步长结果得到明显提高,不仅可以满足仿真精度和数值稳定性要求,同时极大地提高了仿真效率,利于工程推广实现。
另外,根据本发明上述实施例的电力电子化电力系统的多速率仿真装置还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述设定模块还用于:在所述慢系统中,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法将所述慢系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到所述慢系统中选择支路对应数值阻抗;针对所述每个快子系统,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法分别将所述每个快子系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到所述每个快子系统中选择支路对应数值阻抗;定义误差系数为所述慢系统与每个快系统的数值阻抗之比,以根据相对误差最小要求选择速率比。
进一步地,在本发明的一个实施例中,所述获取模块包括:第一获取单元,用于获取所述每个快子系统中时变戴维南等效电路的等效参数;第二获取单元,用于获取接口等值戴维南阻抗矩阵;修正单元,用于考虑多端口之间的耦合影响,修正所述每个快子系统在预设时间内的戴维南等效参数;第三获取单元,用于获取所述慢系统在预设时间内的诺顿等效参数。
进一步地,在本发明的一个实施例中,所述分网模块具体用于通过分网方法进行网络划分,其中,所述分网方法包括节点分裂法、长输电线解耦法、变压器分网法、支路切割法和戴维南诺顿等值分网法。
进一步地,在本发明的一个实施例中,所述慢系统由发电机、变压器、传输线、负荷组成,所述多个快子系统由新能源模块、柔性直流模块和高压直流模块组成。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明一个实施例的电力电子化电力系统的多速率仿真方法的流程图;
图2为根据本发明一个具体实施例的电力电子化电力系统的多速率仿真方法的流程图;
图3为根据本发明一个实施例的电力电子单元与对交流系统进行网络划分示意图;
图4为根据本发明一个实施例的电力电子化电力系统的多速率仿真装置的结构示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
下面参照附图描述根据本发明实施例提出的电力电子化电力系统的多速率仿真方法及装置,首先将参照附图描述根据本发明实施例提出的电力电子化电力系统的多速率仿真方法。
图1是本发明一个实施例的电力电子化电力系统的多速率仿真方法的流程图。
如图1所示,该电力电子化电力系统的多速率仿真方法包括以下步骤:
在步骤S101中,将电力电子化电力系统进行网络划分,以划分为慢系统和多个快子系统,其中交流系统位于慢系统中,电力电子设备位于不同的快子系统中。
其中,在本发明的一个实施例中,通过分网方法进行网络划分,其中,分网方法包括节点分裂法、长输电线解耦法、变压器分网法、支路切割法和戴维南诺顿等值分网法。
可选地,在本发明的一个实施例中,慢系统由发电机、变压器、传输线、负荷组成,多个快子系统由新能源模块、柔性直流模块和高压直流模块组成。
可以理解的是,如图2所示,首先可以采用节点分裂法、长输电线解耦法、变压器分网法、支路切割法和戴维南诺顿等值分网法等分网方法,将电力电子单元与对交流系统进行网络划分,如图3所示。把整个系统划分成一个慢系统S和n个快系统Fp(p=1,2,3,…,N),即多个快子系统。慢系统可以由包括发电机、变压器、传输线、负荷等的交流网络组成,n个快系统可以分别由不同的光伏或风电等新能源模块、柔性直流模块、传统高压直流模块等组成。
在步骤S102中,分别建立慢系统和多个快子系统的等值模型。
进一步地,根据上述的网络划分,整个系统的状态空间方程如下:
其中,tp,p=1,..,N和T分别表示第p个快系统以及慢系统的时间变量,xfp,p=1,..,N和xs分别对应第p个快系统和慢系统的状态变量,ufp,int和ufp,bdry分别对应第p个快系统的内部输入变量和接口变量,us,int和us,bdry分别对应慢系统的内部输入变量和接口变量,Afp,Bfp,int,Bfp,bdry分别对应第p个快系统的参数,As,Bs,int,Bs,bdry代表慢系统的参数。该节点导纳方程采用隐式梯形积分法离散差分化以后,即为慢系统和多个快子系统的节点导纳方程。
在步骤S103中,设定仿真参数,其中,仿真参数包括仿真总时长、慢系统的仿真步长和每个快子系统的仿真步长。
其中,在本发明的一个实施例中,设定仿真参数进一步包括:在慢系统中,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法将慢系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到慢系统中选择支路对应数值阻抗;针对每个快子系统,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法分别将每个快子系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到每个快子系统中选择支路对应数值阻抗;定义误差系数为慢系统与每个快系统的数值阻抗之比,以根据相对误差最小要求选择速率比。
可以理解的是,设定与仿真时间相关的参数,包括仿真总时长Ttotal;慢系统仿真步长hs和n个快系统的仿真步长hfp(p=1,2,3,…,n),其中,hs=mphfp,mp表示第p个快系统与慢系统之间的速率比,选择速率比mp的具体过程包括以下步骤:
S1:针对第p个子系统,在子系统中,选择任意支路(包括RL支路,RC支路,传输线等),建立该支路的电压电流微分方程。
举例而言,以RL支路为例,该支路对应的电压电流微分方程为:
其中,u和i分别为该支路相对应端口的电压值和电流值,R,L为支路参数,Tp对应第p个子系统的时间变量。
S2:在慢系统中,采用隐式梯形法,将上述方程转换为离散域下的差分方程如下:
其中,nk、nk+1分别为慢系统第k、k+1个仿真步长的对应时刻,u(nk+1)、u(nk)分别对应nk、nk+1时刻的电压值,i(nk)、i(nk+1)分别为nk、nk+1时刻的电流值。
利用离散傅里叶变换,得到慢系统中该支路对应数值阻抗Zs如下:
其中, 分别对应慢系统中该元件的电压和电流的离散傅里叶系数。
S3:在每一个快系统Fp(p=1,2,3,…,n)中,将上述步骤(3-1)的微分方程转换为上述步骤S2中,该慢系统积分步长内的差分方程:
利用离散傅里叶变换,在每一个快系统Fp中,得到该支路的数值阻抗Zfp(p=1,2,3,…,n)如下:
其中, 分别对应快系统Fp中该元件的电压和电流的离散傅里叶系数。
S4:定义误差系ζ(mp,w)=Zs/Zfp,(p=1,2,3,...N),根据相对误差的取值范围,选择速率比mp:
其中,fNy为系统对应的奈奎斯特频率;
mphfp=hs≤hs,max,p=1,2,3,…N;
其中,hs,max为由数值稳定条件决定的慢系统的最大仿真步长。
在步骤S104中,通过对电力系统进行潮流计算得到电力系统的稳态初始值。
也就是说,对整个系统进行潮流计算,确定整个系统的稳态初始值。
在步骤S105中,建立多个快子系统与慢系统间的接口,并且在每个快子系统中建立时变戴维南等效电路,并且在慢系统中建立时变诺顿等效电路。
具体地,建立快系统与慢系统间的接口,在每个快系统Fp(p=1,2,3,…N)中,建立时变戴维南等效电路,表征慢系统对第p个快系统的影响,其中戴维南等效电路参数包括和(p=1,2,3,…N);在慢系统S中,建立时变诺顿等效电路,表征第p个快系统对慢系统的影响,诺顿等效电路参数为和(p=1,2,3,…N)。
在步骤S106中,根据每个快子系统中建立的时变戴维南等效电路和慢系统中建立的时变诺顿等效电路得到接口参数,并通过并行求解每个快子系统和慢系统的等值模型对应的节点导纳方程,直到仿真结束,进而得到电力电子化电力系统的电磁暂态仿真结果。
进一步地,在本发明的一个实施例中,根据每个快子系统中建立的时变戴维南等效电路得到接口参数,进一步包括:获取每个快子系统中时变戴维南等效电路的等效参数;获取接口等值戴维南阻抗矩阵;考虑多端口之间的耦合影响,修正每个快子系统在预设时间内的戴维南等效参数;获取慢系统在预设时间内的诺顿等效参数。其中,预设时间可以根据实际情况进行设置。
具体而言,假定慢系统从第k步到第k+1步仿真,即仿真时间从nk时刻步进到nk+1时刻,计算接口参数,具体过程包括以下步骤:
S1:求解每个快系统Fp(p=1,2,3,…N)中,戴维南电路的等效参数(i=1,2,...mp)和(i=1,2,...mp),具体过程包括以下步骤:
S11,由于(i=1,2,...mp)在慢系统积分步长hs内的中间值不被程序计算,所以采用三次样条插值和滑动窗口算法,预测下一个快系统积分步长hfp时刻的最佳值(i=1,2,...mp),为了表达简洁,省去括号中的时间变量,例如简写为
S12,判断慢系统的网络拓扑结构是否发生变化,若发生变化,则根据网络节点阻抗矩阵,重新求解p=1,2,3,…N。
S2:假定共有N个接口,求解接口等值戴维南阻抗矩阵,记作[Zint],该矩阵行数等于N,列数等于N。
S3:考虑多端口之间的耦合影响,修正每个快系统Fp(p=1,2,3,…,N)中,戴维南等效参数对应第p个接口的修正公式如下:
其中,表示第k个接口的开路电压值,在上述步骤S1中获得;表示第n个接口的支路电流,由该步长下对应的仿真结果获取。
S4:求解慢系统在[nk,nk+1]时间内的诺顿等效参数和(p=1,2,3,…N),具体过程包括以下步骤:
S41,的值通过在所有采样点中每隔t个点选择一个点,最后计算被选择点的平均值获得,如下:
其中,[n/t]是取整函数,t代表每隔t个点取一个值。
S42,判断任意一个快系统的网络拓扑结构是否发生变化,若发生变化,则根据网络节点导纳矩阵,重新求解(p=1,2,3,…N)。
最后,快系统与慢系统并行求解上述给出的各个子系统的节点导纳方程,重复仿真,直到仿真时间到达设定时间Ttotal,整个过程结束,最终得到含有电力电子单元电力系统的电磁暂态仿真结果。
在本发明的实施例中,本发明实施例的电力电子化电力系统的多速率仿真方法既能提高仿真规模,包含交流系统和大规模电力电子单元,并且可以保证仿真效率,满足含有大规模电力电子单元电力系统的安全稳定分析的要求,可以根据耦合最小原则,将电力电子单元与交流系统切分成不同子系统。交流系统采用大步长进行仿真,构成慢系统,电力电子单元采用不同的小步长进行仿真,构成多个快系统,两者之间通过接口戴维南/诺顿等值模型进行交互,实现全网仿真,不仅可以满足仿真精度和数值稳定性要求,同时极大地提高了仿真效率,利于工程推广实现。
根据本发明实施例提出的电力电子化电力系统的多速率仿真方法,由于根据系统多尺度特性,将整个系统划分成多个子系统,采用不同速率进行全网仿真。其中交流系统采用大步长仿真,电力电子单元采用不同的小步长进行仿真,该方法可以极大地提高仿真效率,仿真规模得到进一步提升,适用于含有大规模电力电子单元的电力系统仿真,通过对若干个子系统并行地进行电磁暂态仿真,并采用接口技术在不同子网之间传输仿真数据,可实现全网并行计算,尤其适用于含有模块化多电平模块、高压直流模块等多电力电子开关器件的电力系统电磁暂态仿真计算,具有极大的工程实用价值。
其次参照附图描述根据本发明实施例提出的电力电子化电力系统的多速率仿真装置。
图4是本发明一个实施例的电力电子化电力系统的多速率仿真装置的结构示意图。
如图4所示,该电力电子化电力系统的多速率仿真装置10包括:分网模块100、第一建立模块200、设定模块300、计算模块400、第二建立模块500和获取模块600。
其中,分网模块100用于将电力电子化电力系统进行网络划分,以划分为慢系统和多个快子系统,其中交流系统位于慢系统中,电力电子设备位于不同的快子系统中。第一建立模块200用于分别建立慢系统和多个快子系统的等值模型。设定模块300用于设定仿真参数,其中,仿真参数包括仿真总时长、慢系统的仿真步长和每个快子系统的仿真步长。计算模块400用于通过对电力系统进行潮流计算得到电力系统的稳态初始值。第二建立模块500用于建立多个快子系统与慢系统间的接口,并且在每个快子系统中建立时变戴维南等效电路,并且在慢系统中建立时变诺顿等效电路。获取模块600用于根据每个快子系统中建立的时变戴维南等效电路和慢系统中建立的时变诺顿等效电路得到接口参数,并通过并行求解每个快子系统和慢系统的等值模型对应的节点导纳方程,直到仿真结束,进而得到电力电子化电力系统的电磁暂态仿真结果。本发明实施例的电力电子化电力系统的多速率仿真装置10可以结合含有电力电子单元的电力系统自身多尺度特性,采用多速率仿真,不仅可以保留交流系统或者直流电网内部的非线性动态特性,同时,仿真效率相对于单一步长结果得到明显提高,不仅可以满足仿真精度和数值稳定性要求,同时极大地提高了仿真效率,利于工程推广实现。
进一步地,在本发明的一个实施例中,设定模块300还用于:在慢系统中,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法将慢系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到慢系统中选择支路对应数值阻抗;针对每个快子系统,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法分别将每个快子系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到慢系统中选择支路对应数值阻抗;定义误差系数为慢系统与每个快系统的数值阻抗之比,以根据相对误差最小要求选择速率比。
进一步地,在本发明的一个实施例中,获取模块600包括:第一获取单元,用于获取每个快子系统中时变戴维南等效电路的等效参数;第二获取单元,用于获取接口等值戴维南阻抗矩阵;修正单元,用于考虑多端口之间的耦合影响,修正每个快子系统在预设时间内的戴维南等效参数;第三获取单元,用于获取慢系统在预设时间内的诺顿等效参数。
进一步地,在本发明的一个实施例中,分网模块100具体用于通过分网方法进行网络划分,其中,分网方法包括节点分裂法、长输电线解耦法、变压器分网法、支路切割法和戴维南诺顿等值分网法。
进一步地,在本发明的一个实施例中,慢系统由发电机、变压器、传输线、负荷组成,多个快子系统由新能源模块、柔性直流模块和高压直流模块组成。
需要说明的是,前述对电力电子化电力系统的多速率仿真方法实施例的解释说明也适用于该实施例的电力电子化电力系统的多速率仿真装置,此处不再赘述。
根据本发明实施例提出的电力电子化电力系统的多速率仿真装置,由于根据系统多尺度特性,将整个系统划分成多个子系统,采用不同速率进行全网仿真。其中交流系统采用大步长仿真,电力电子单元采用不同的小步长进行仿真,该方法可以极大地提高仿真效率,仿真规模得到进一步提升,适用于含有大规模电力电子单元的电力系统仿真,通过对若干个子系统并行地进行电磁暂态仿真,并采用接口技术在不同子网之间传输仿真数据,可实现全网并行计算,尤其适用于含有模块化多电平模块、高压直流模块等多电力电子开关器件的电力系统电磁暂态仿真计算,具有极大的工程实用价值。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
在本发明中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度小于第二特征。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (6)
1.一种电力电子化电力系统的多速率仿真方法,其特征在于,包括以下步骤:
将电力电子化电力系统进行网络划分,以划分为慢系统和多个快子系统,其中交流系统位于所述慢系统中,电力电子设备位于不同的快子系统中;
分别建立所述慢系统和所述多个快子系统的等值模型;
设定仿真参数,其中,所述仿真参数包括仿真总时长、所述慢系统的仿真步长和所述每个快子系统的仿真步长,其中所述设定仿真参数进一步包括所述设定仿真参数还包括在所述慢系统中,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法将所述慢系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到所述慢系统中选择支路对应数值阻抗;针对所述每个快子系统,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法分别将所述每个快子系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到所述每个快子系统中选择支路对应数值阻抗;定义误差系数为所述慢系统与每个快系统的数值阻抗之比,以根据相对误差最小要求选择速率比;
通过对电力系统进行潮流计算得到所述电力系统的稳态初始值;
建立所述多个快子系统与所述慢系统间的接口,并且在所述每个快子系统中建立时变戴维南等效电路,并且在所述慢系统中建立时变诺顿等效电路;以及根据所述每个快子系统中建立的所述时变戴维南等效电路和所述慢系统中建立的所述时变诺顿等效电路得到接口参数,并通过并行求解所述每个快子系统和慢系统的等值模型对应的节点导纳方程,直到仿真结束,进而得到电力电子化电力系统的电磁暂态仿真结果,其中,所述根据所述每个快子系统中建立的所述时变戴维南等效电路和所述慢系统中建立的所述时变诺顿等效电路得到接口参数进一步包括获取所述每个快子系统中时变戴维南等效电路的等效参数;获取接口等值戴维南阻抗矩阵;考虑多端口之间的耦合影响,修正所述每个快子系统在预设时间内的戴维南等效参数;获取所述慢系统在预设时间内的诺顿等效参数。
2.根据权利要求1所述的电力电子化电力系统的多速率仿真方法,其特征在于,通过分网方法进行网络划分,其中,所述分网方法包括节点分裂法、长输电线解耦法、变压器分网法、支路切割法和戴维南诺顿等值分网法。
3.根据权利要求1所述的电力电子化电力系统的多速率仿真方法,其特征在于,所述慢系统由发电机、变压器、传输线、负荷组成,所述多个快子系统由新能源模块、柔性直流模块和高压直流模块组成。
4.一种电力电子化电力系统的多速率仿真装置,其特征在于,包括:
分网模块,用于将电力电子化电力系统进行网络划分,以划分为慢系统和多个快子系统,其中交流系统位于所述慢系统中,电力电子设备位于不同的快子系统中;
第一建立模块,用于分别建立所述慢系统和所述多个快子系统的等值模型;
设定模块,用于设定仿真参数,其中,所述仿真参数包括仿真总时长、所述慢系统的仿真步长和所述每个快子系统的仿真步长,所述设定模块还用于在所述慢系统中,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法将所述慢系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到所述慢系统中选择支路对应数值阻抗;针对所述每个快子系统,选择任意支路以建立选择支路的电压电流微分方程,并且采用隐式梯形法分别将所述每个快子系统的电压电流微分方程转换为离散域下的差分方程,并利用离散傅里叶变换,以得到所述每个快子系统中选择支路对应数值阻抗;定义误差系数为所述慢系统与每个快系统的数值阻抗之比,以根据相对误差最小要求选择速率比;
计算模块,用于通过对电力系统进行潮流计算得到所述电力系统的稳态初始值;
第二建立模块,用于建立所述多个快子系统与所述慢系统间的接口,并且在所述每个快子系统中建立时变戴维南等效电路,并且在所述慢系统中建立时变诺顿等效电路;以及获取模块,用于根据所述每个快子系统中建立的所述时变戴维南等效电路和所述慢系统中建立的所述时变诺顿等效电路得到接口参数,并通过并行求解所述每个快子系统和慢系统的等值模型对应的节点导纳方程,直到仿真结束,进而得到电力电子化电力系统的电磁暂态仿真结果,其中,所述获取模块包括第一获取单元,用于获取所述每个快子系统中时变戴维南等效电路的等效参数;第二获取单元,用于获取接口等值戴维南阻抗矩阵;修正单元,用于考虑多端口之间的耦合影响,修正所述每个快子系统在预设时间内的戴维南等效参数;第三获取单元,用于获取所述慢系统在预设时间内的诺顿等效参数。
5.根据权利要求4所述的电力电子化电力系统的多速率仿真装置,其特征在于,所述分网模块具体用于通过分网方法进行网络划分,其中,所述分网方法包括节点分裂法、长输电线解耦法、变压器分网法、支路切割法和戴维南诺顿等值分网法。
6.根据权利要求4所述的电力电子化电力系统的多速率仿真装置,其特征在于,所述慢系统由发电机、变压器、传输线、负荷组成,所述多个快子系统由新能源模块、柔性直流模块和高压直流模块组成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610803055.5A CN106372339B (zh) | 2016-09-05 | 2016-09-05 | 电力电子化电力系统的多速率仿真方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610803055.5A CN106372339B (zh) | 2016-09-05 | 2016-09-05 | 电力电子化电力系统的多速率仿真方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106372339A CN106372339A (zh) | 2017-02-01 |
CN106372339B true CN106372339B (zh) | 2019-08-09 |
Family
ID=57900020
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610803055.5A Active CN106372339B (zh) | 2016-09-05 | 2016-09-05 | 电力电子化电力系统的多速率仿真方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106372339B (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107069811B (zh) * | 2017-04-12 | 2019-07-26 | 清华大学 | 基于同步参考坐标系的阻抗网络建模与稳定性分析方法 |
CN107066732A (zh) * | 2017-04-13 | 2017-08-18 | 清华大学 | 面向交直流混合网络的分布式混合仿真分析方法及装置 |
CN107368451B (zh) * | 2017-06-14 | 2019-09-17 | 清华大学 | 电路状态空间方程计算方法、装置和系统 |
CN107832504B (zh) * | 2017-10-24 | 2021-04-16 | 南方电网科学研究院有限责任公司 | 仿真多端直流系统的最大转换电流的方法及系统 |
CN108646587B (zh) * | 2018-04-18 | 2021-08-17 | 许继电气股份有限公司 | 一种emtdc并行仿真方法及系统 |
CN109165426A (zh) * | 2018-08-06 | 2019-01-08 | 清华四川能源互联网研究院 | 电磁暂态仿真方法和电磁暂态仿真装置 |
CN109066669B (zh) * | 2018-09-12 | 2020-06-12 | 清华四川能源互联网研究院 | 电力系统仿真方法、电力系统仿真装置及电子设备 |
CN109063408B (zh) * | 2018-10-31 | 2022-04-01 | 邓宏伟 | 一种基于直算法的电磁机电暂态仿真方法 |
CN109782625B (zh) * | 2018-12-20 | 2023-11-14 | 中国电力科学研究院有限公司 | 一种电路模型的实时仿真方法及系统 |
CN110362862A (zh) * | 2019-06-11 | 2019-10-22 | 中国电力科学研究院有限公司 | 一种电网信息物理系统物理层优化仿真方法和系统 |
CN110516314A (zh) * | 2019-08-02 | 2019-11-29 | 上海交通大学 | 用于大规模风电场的多速率多模态联合仿真方法 |
CN110569558B (zh) * | 2019-08-16 | 2021-09-07 | 上海交通大学 | 适用于微电网实时仿真的混合电磁暂态仿真方法 |
US11476667B2 (en) | 2019-08-16 | 2022-10-18 | Shanghai Jiao Tong University | Hybrid electromagnetic transient simulation method for microgrid real-time simulation |
CN110851998A (zh) * | 2019-12-06 | 2020-02-28 | 华北电力科学研究院有限责任公司 | 换流站分层仿真测试系统 |
CN111628495B (zh) * | 2020-05-13 | 2022-09-02 | 上海交通大学 | 含vsc-mmc交直流电网系统的电磁暂态仿真方法 |
CN111898234A (zh) * | 2020-05-18 | 2020-11-06 | 天津大学 | 一种适用于frtds的燃气机发电系统仿真建模方法 |
CN112182815A (zh) * | 2020-09-16 | 2021-01-05 | 天津大学 | 一种基于诺顿等效的多速率实时仿真方法 |
CN113706094B (zh) * | 2021-07-29 | 2024-02-20 | 国电南瑞科技股份有限公司 | 一种基于消息总线的综合能源实时协同仿真系统及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102386627A (zh) * | 2011-11-15 | 2012-03-21 | 中国海洋石油总公司 | 一种动态无功补偿器的控制系统及其控制方法 |
CN104393602A (zh) * | 2014-12-12 | 2015-03-04 | 新奥科技发展有限公司 | 一种分布式能源网络的调节方法及装置 |
CN104578402A (zh) * | 2013-10-17 | 2015-04-29 | 通用电气公司 | 用于控制电力网络的方法和系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8907655B2 (en) * | 2011-04-29 | 2014-12-09 | Analog Devices, Inc. | System and method for detecting a fundamental frequency of an electric power system |
-
2016
- 2016-09-05 CN CN201610803055.5A patent/CN106372339B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102386627A (zh) * | 2011-11-15 | 2012-03-21 | 中国海洋石油总公司 | 一种动态无功补偿器的控制系统及其控制方法 |
CN104578402A (zh) * | 2013-10-17 | 2015-04-29 | 通用电气公司 | 用于控制电力网络的方法和系统 |
CN104393602A (zh) * | 2014-12-12 | 2015-03-04 | 新奥科技发展有限公司 | 一种分布式能源网络的调节方法及装置 |
Non-Patent Citations (1)
Title |
---|
基于变量自适应分组的电力系统多速率仿真算法;高毅 等;《电力自动化设备》;20110228;第31卷(第2期);1-6,25 * |
Also Published As
Publication number | Publication date |
---|---|
CN106372339A (zh) | 2017-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106372339B (zh) | 电力电子化电力系统的多速率仿真方法及装置 | |
CN107577870B (zh) | 基于同步相量量测的配电网电压功率灵敏度鲁棒估计方法 | |
CN109802406B (zh) | 一种分析柔性直流输电接入系统谐振稳定性的方法 | |
CN107342586B (zh) | 一种用于确定配电网静态等值的方法及装置 | |
CN103810646B (zh) | 一种基于改进投影积分算法的有源配电系统动态仿真方法 | |
CN108448631B (zh) | 基于模型降阶的含分布式电源接入配电网的动态相量建模方法 | |
CN105938578A (zh) | 一种基于聚类分析的大规模光伏电站等值建模方法 | |
CN103077268B (zh) | 面向电力系统电磁暂态仿真的状态空间自动建模方法 | |
CN108429294B (zh) | 一种含能量路由器的交直流网络潮流模型及求解方法 | |
CN108959671B (zh) | 半桥和全桥型模块化多电平换流器的实时仿真建模方法 | |
CN103336866A (zh) | 一种电磁暂态仿真中含负电阻支路的处理方法 | |
CN107066732A (zh) | 面向交直流混合网络的分布式混合仿真分析方法及装置 | |
CN104899396A (zh) | 一种修正系数矩阵的快速分解法潮流计算方法 | |
CN112505473A (zh) | 柔性直流电网双极短路故障暂态电流的解析计算方法 | |
CN109149566B (zh) | 一种大功率缺失下频率最低点预测的仿真模型的建模方法 | |
CN111339624B (zh) | 基于psasp和emtp/atp短路电流直流分量计算方法 | |
CN104135038A (zh) | 一种交直流混联系统不对称故障分析方法 | |
CN112670992B (zh) | 含能量路由器的配电网稳定性分析和失稳校正方法及系统 | |
CN112072692A (zh) | 一种新能源发电场站的阻抗等值方法及装置 | |
CN113078670A (zh) | 混合级联直流输电作用下受端电网谐振稳定性的评估方法 | |
CN106849154B (zh) | 一种光伏电站节点谐波电压与谐波源灵敏度关系分析方法 | |
CN102709911A (zh) | 一种混合仿真谐波特性的接口设计方法 | |
CN109829178B (zh) | 一种基于模型分割法的变电站实时仿真方法 | |
CN112036010A (zh) | 一种基于数据驱动的光伏系统动态过程混合等效建模方法 | |
CN105095659B (zh) | 基于云计算的省地协调分布式状态估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |