CN106021768B - 含分布式电源接入的配电网简化建模方法 - Google Patents
含分布式电源接入的配电网简化建模方法 Download PDFInfo
- Publication number
- CN106021768B CN106021768B CN201610371147.0A CN201610371147A CN106021768B CN 106021768 B CN106021768 B CN 106021768B CN 201610371147 A CN201610371147 A CN 201610371147A CN 106021768 B CN106021768 B CN 106021768B
- Authority
- CN
- China
- Prior art keywords
- equation
- distribution network
- state
- current
- generation resource
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提出了含分布式电源接入的配电网简化建模方法。大规模分布式电源接入配电网,不仅分布式电源的非线性特征使得仿真计算速度受到很大限制,而且当配电网的节点、支路数较多时,网络矩阵维数的过高会进一步导致仿真计算速度下降。本发明将含分布式电源接入的配电网划分为内部线性区域和外部非线性区域,将配电网网架部分与分布式电源分别进行简化处理,得到含分布式电源接入的配电网新建模方法,可准确地模拟系统的动态过程,大大提升系统的仿真速度。该建模方法保证了仿真的准确性和快速性,为大规模分布式电源接入配电网的快速仿真提供了新技术。
Description
技术领域
本发明属于电力系统领域,具体地说是含分布式电源接入的配电网简化建模方法。
背景技术
随着经济社会的不断发展,能源消耗不断增加,不可再生资源逐渐枯竭,并且由此引发的雾霾等环境污染问题愈加严重,这直接促使发展新能源成为各国今后能源发展的战略需求。在此背景下,分布式电源(Distributed Generation,DG)近年来发展迅猛。然而大规模分布式电源接入传统配电网,会对配电网的电力电量平衡、电气设备经济运行、线路潮流、电压降落、电能损耗、电能质量、开关设备保护等情况产生重要影响,因此必须对含分布式电源接入的配电网进行建模仿真研究。
含分布式电源接入的配电网可以分为线性的配电网网架部分以及非线性的分布式电源部分。对于线性的配电网网架部分,当配电网节点、支路数较多时,会由于其网络矩阵维数较高,导致直接进行分析相对困难或数值模拟耗时过长,模型的求解时间相当可观,因此对于大型或复杂的配电网系统,有必要寻求对模型的规模或阶数进行有效降阶的方法,减少数据运算量,达到减少仿真耗时的目的。而非线性的分布式电源,由于含有大量非线性特征的电力电子器件及其控制器,仿真计算速度受到很大限制,当分布式电源大规模接入配电网时,现有的仿真技术手段会使仿真耗时过长,甚至无法进行,因此对于分布式电源,也有必要探索在不牺牲仿真精度情况下加快仿真速度、减少计算量的方法,为此可对分布式电源进行简化建模。当线性的配电网网架部分以及非线性的分布式电源模型建立后,两种仿真技术的合成也需要研究。
发明内容
本发明所要解决的技术问题是克服上述现有技术存在的缺陷,提供一种含分布式电源接入的配电网简化建模方法,该方法将配电网线性网架部分与含非线性特征的分布式电源分别进行简化处理,得到含分布式电源接入的配电网新仿真模型,可准确模拟系统的动态过程,大大提升系统的仿真速度。
为此,本发明采用如下的技术方案:将含分布式电源的配电网分为两个区域:1)内部线性区域,即配电网网架部分;2)外部非线性区域,即分布式电源部分。内部区域主要是由系统元件组成,包括线路、变压器,将负荷视作恒阻抗负荷。外部非线性分布式电源假设采用恒功率控制方式,以最大可能输出有功功率。先对内部线性的配电网网架区域推导以基本回路电流和电容电压为状态变量的配电网状态空间方程,并采用降阶方法实现线性系统的简化处理,在不改变系统的输入输出个数的情况下减少状态方程中的状态量数量,得到降阶后的配电网状态空间方程;对外部分布式电源区域建立微分代数方程,得到分布式电源的状态空间方程;并且用动态相量法提取内部线性区域以及外部非线性区域状态空间方程的主要频率特征量;最后以网络接口的形式对接,对整个系统进行互联,建立含分布式电源接入的配电网新仿真模型。
本发明的有益效果是:本发明充分考虑了含分布式电源配电网的特点,根据配电网中元件的线性和非线性特征,将配电网划分为两个研究区域分别进行处理,提供了一种配电网简化建模方法,降低了模型复杂度,加快了仿真速度。该简化建模方法保证了仿真的准确性和快速性,为大规模分布式电源接入配电网的快速仿真奠定了基础。
本发明采用以下具体步骤:
步骤1),将含分布式电源接入的配电网划分成内部、外部研究区域;
将含分布式电源接入的配电网划分为两个区域,一个是以线性的配电网网架部分为对象的内部区域和以非线性的分布式电源为对象的外部区域。
步骤2),对内部线性区域建立配电网网架部分的状态空间方程;
在对内部线性区域的建模中,由系统初始的节点关联矩阵推导配电网络的基本回路矩阵,再利用基本回路矩阵得到配电网络的基尔霍夫电流方程(KCL)、基尔霍夫电压方程(KVL)。采用ASMG方法中提出的一般支路模型,将系统中所有支路的伏安关系写成相量形式。定义相应的关联矩阵,使输入系统的电压量和电流量与一般支路模型对应,得到输入量与支路量的对应关系。对于配电网,分布式电源的输入量为电流量,将模型输入量定义为两种类型:电压输入端子和电流输入端子,并定义电压输入端子电流与电流输入端子电压为模型输出量,建立系统状态方程与输出方程,进而得到以基本回路电流和电容电压为状态变量的配电网状态空间方程,并对照标准状态方程式,得到对应的系数矩阵表达式。
步骤3),对配电网网架部分的状态空间方程降阶,求取其降阶后的状态空间方程;
在建立完配电网网架部分的状态空间方程后,对于线性系统,采用模型降阶处理方法,对模型进行降阶。即,对线性非时变的动态系统,根据其状态空间方程,计算该系统的拉普拉斯域下的传递函数,得到Krylov子空间以及子空间的正交特征向量。再利用SPRIM算法获取特征矩阵,即将子空间的正交特征向量拆分为子向量,进而由以基本回路电流和电容电压为状态变量的配电网状态空间方程得到降阶后的状态空间方程。
步骤4),对外部非线性区域建立分布式电源的状态空间方程;
对于外部非线性的分布式电源,假设逆变器型分布式电源采用PQ控制策略。在功率控制中,计算分布式电源的瞬时有功和无功功率,功率控制利用PI调节,分别列写功率PI调节的状态方程及代数方程,进而得到以电流环电流指令值为输出的功率PI调节的状态空间方程;电流控制也利用PI调节,分别列写电流PI调节的状态方程及代数方程,进而得到以逆变器出口电压指令值为输出的电流PI调节的状态空间方程;对于网侧滤波器,列写以电网侧电流为输出的状态方程,进而得到其状态空间方程。之后,将上述功率控制、电流控制、网侧滤波器的状态空间方程整理合并,得到整个分布式电源的状态空间方程。
步骤5),对内部、外部区域的状态空间方程进行动态相量建模;
运用基于信号调制理论的动态相量模型,即对时域内周期为T的信号选择宽度为T的时间窗内的主要傅立叶系数,将配电网内部区域的配电网状态空间方程以及外部区域的分布式电源状态空间方程改写成动态相量的表达形式,对原始波形作近似的估计,进一步降低模型的复杂度。由于上述分布式电源的控制均是基于dq0坐标系,进一步将分布式电源的动态相量模型转换到abc坐标系下以实现内部区域与外部区域的互联。
步骤6),内部、外部区域的动态相量模型互联;
将动态相量表达形式的内部区域的配电网状态空间方程和外部区域的分布式电源状态空间方程以网络接口的形式对接,即内部区域状态空间方程的输入量为外部区域状态空间方程的输出量和外部电源的电压量,外部区域状态空间方程的输入量为内部区域状态空间方程的输出量。
本发明提出了一种含分布式电源接入配电网的新简化建模方法,可加快仿真速度,减少仿真耗时。利用本发明提出的含分布式电源接入配电网的简化建模方法,可推广应用于分布式电源高渗透的大型配电网系统的技术验证和分析。
附图说明
图1为含分布式电源接入配电网的内外区域划分图。
图2为配电网一般支路模型。
图3为分布式电源的控制框图。
图4为IEEE34节点配电网系统接线图。
图5a为功率指令为0.5MW时大电源接入处的电流Is以及DG接入处的电流Idg波形图。
图5b为功率指令波动时大电源接入处的电流Is,DG接入处的电流Idg波形图。
图5c为节点828发生三相接地短路故障时大电源接入处的电流Is以及DG接入处的电流Idg波形图。
具体实施方式
以下结合附图对本发明作进一步说明,本发明包括以下步骤:
根据配电网中元件的线性和非线性特征,将含分布式电源的配电网分为两个区域:1)内部线性区域,即配电网网架部分;2)外部非线性区域,即分布式电源部分,如图1所示,其中vvs为外部电压源,ics为外部电流源。内部区域主要是由系统元件组成,包括线路、变压器。在暂态建模过程中,将负荷视作恒阻抗负荷。同时,系统中分布式电源假设采用恒功率控制方式,以最大可能输出有功功率。因此由于分布式电源的控制作用,使得系统的外部区域呈现非线性特性。本发明先对内部线性的配电网网架部分推导以基本回路电流和电容电压为状态变量的配电网状态空间方程,并采用降阶方法实现线性系统的简化处理,在不改变系统输入输出个数的情况下减少状态方程中的状态变量数量,得到降阶后的配电网状态空间方程;对外部分布式电源区域建立微分代数方程,得到分布式电源的状态空间方程;并且用动态相量法提取其中的主要频率特征量;最后以网络接口的形式对接,对整个系统进行互联,建立含分布式电源接入的配电网新仿真模型。具体实施步骤如下:
步骤1),将含分布式电源接入的配电网划分成内部、外部研究区域,
将含分布式电源接入的配电网划分为两个区域,一个是以线性的配电网网架部分为对象的内部区域和以非线性的分布式电源为对象的外部区域。
步骤2),对内部线性区域建立配电网网架部分的状态空间方程,
利用基本回路矩阵Bb表述配电网络的基尔霍夫电流方程(KCL)和基尔霍夫电压方程(KVL),如下式所示:
Bbvbr=0
其中ibr为支路电流向量;vbr为支路电压向量;ix为基本回路电流向量。基本回路矩阵Bb由系统初始的节点关联矩阵推导得到。
采用ASMG方法中提出的一般支路模型,如图2所示。图中ik和jk分别代表第k条支路的支路电流和电流源电流。对于系统中的所有支路将伏安关系写成向量形式:
vbr=Rbribr+Lbrpibr+Pbrqbr+ebr
其中,Pbr主对角线元素为相应支路电容的倒数,若该支路不存在电容则对应位置为0;p为微分算子;jbr为支路电流源电流向量;Lbr是支路电感矩阵,其主对角线元素表示相应支路的自感,非主对角线元素表示对应支路的互感;Rbr是支路电阻矩阵,ebr为支路电压源电压向量。
为了让输入系统的电压量和电流量与一般支路模型对应,定义相应的关联矩阵W1∈Rn_vs×n_br,W2∈Rn_cs×n_br,其中n_vs、n_cs分别表示电压、电流输入量的个数,n_br表示网络中支路个数。当第r个电压输入量位于支路k上时,则W1(r,k)=1,否则W1(r,k)=0。W2采用类似方式定义,即当第r个电流输入量位于支路k上时,则W2(r,k)=1,否则W2(r,k)=0。可得输入量与支路量的对应关系为:
其中ebr为支路电压源电压向量,jbr为支路电流源电流向量,vvs为输入电压向量,ics为输入电流向量。
对于配电网来说,分布式电源的输入量为电流量。因此为了不失一般性,将模型输入量定义为2种类型,即电压输入端子和电流输入端子,如图1所示。结合实际情况,定义电压输入端子电流ics与电流输入端子电压vvs为模型输出量,可以得到以基本回路电流ix和电容电压vc为状态变量的配电网状态空间方程如下式:
其中ix为基本回路电流向量,vc为支路电容电压向量,Cbr为支路电容矩阵,支路电容关联矩阵M∈Rn_br×n_c,n_c表示网络中电容个数,n_br表示网络中支路个数,取值方式与W1和W2相同。对照标准状态方程式:
y=LTx
可得标准状态方程对应的系数矩阵表达式如下:
步骤3),对配电网网架部分的状态空间方程降阶,求取其降阶后的状态空间方程,
在建立完内部配电网架部分的状态空间方程后,对模型进行降阶处理。以下为具体降阶方法:
假设某个线性时不变的动态系统:
y(t)=BTx(t)
计算得到该系统的拉普拉斯域下的传递函数如下:
H(s)=BT(s0E-A+(s-s0)E)-1B
=BT(I+(s-s0)M)-1R
其中
M:=(s0E-A)-1E,R:=(s0E-A)-1B
n阶Krylov子空间可以通过下面得到:
κn(M,R)=colspan[R MR M2R … MnR]
进而得到空间的正交特征向量Vn:
colspanVn=κn(M,R)
再利用SPRIM算法,获取特征矩阵
将Vn根据状态变量中的电压、电流量拆分为子向量和其中 分别代表状态变量中的电压和电流量,利用和可以构造得到特征矩阵
则以基本回路电流和电容电压为状态变量的配电网状态空间方程可以变换为下面降阶后的状态空间方程:
其中,上式方程加下标n表示降阶后的状态空间方程,该方程矩阵的阶数小于降阶前的系统。上述状态方程对应的系数矩阵表达式如下:
步骤4),对外部非线性区域建立分布式电源的状态空间方程,
对于外部非线性区域的分布式电源,它由三个子模块组成:功率控制、电流控制以及网侧滤波器,分别建立其状态空间方程。图3为逆变器型分布式电源的控制框图。
对于功率控制,假设分布式电源采用恒功率控制策略(PQ控制),瞬时有功功率p和无功功率q由下式计算得到:
p=vodiod+voqioq
q=vodioq-voqiod
其中vod、iod、voq、ioq分别为电网电压、电流在两相同步旋转坐标系dq下的d轴电压、电流和q轴电压、电流。
功率控制利用PI调节,则对应的状态方程如下:
以及代数方程为:
其中,Pref和Qref分别是有功、无功功率指令,和是状态变量,和为电流控制的输入参考值,Kp和Ki分别是功率控制中的比例调节增益和积分调节增益。
则功率控制的状态空间方程可列写如下,方程的输入量为参考功率和反馈功率:
其中:
[PQref]=[Pref Qref]T,[pq]=[pq]T,
对于电流控制,对应状态方程为:
以及代数方程为:
其中,φd和φq是状态变量,和电流控制的输出参考值,ω为电网角频率,Lc为逆变器出口侧滤波电感的大小,Kpv和Kiv分别是电流控制中的比例调节增益和积分调节增益。
则电流控制的状态空间方程可列写如下,方程的输入量为参考电流和反馈电流:
其中:
[iodq]=[iod ioq]T,[vodq]=[vod voq]T
对于网侧滤波器,在假设逆变器实际提供输出电压vid、viq等于指令电压的情况下,滤波器模型可列写如下:
则网侧滤波器的状态空间方程为:
其中:
[iodq]=[iodioq]T,[vidq]=[vidviq]T,[vodq]=[vodvoq]T
将上述功率控制、电流控制、网侧滤波器的状态空间方程整理合并,得到整个分布式电源的状态空间方程:
[idq]=Ct[xDG]
其中:
步骤5),对内部、外部区域的状态空间方程进行动态相量建模,
运用基于信号调制理论的动态相量模型,即对时域内周期为T的信号选择宽度为T的时间窗内的主要傅立叶系数,对配电网内部区域的以基本回路电流和电容电压为状态变量的配电网状态空间方程以及外部区域的分布式电源状态空间方程改写成动态相量的表达形式,分别计算如下:
<idq>k=Ct<xDG>k
其中,<·>k代表相关变量的k阶动态相量。
由于上述分布式电源的控制均是基于dq0坐标系,则将分布式电源的动态相量模型转换为abc坐标系以实现内部区域与外部区域的互联,转换关系满足以下方程:
其中,f代表可能的电量参数,例如电压或电流,DQ代表同步旋转坐标系,p、n分别代表正序和负序。
步骤6),内部、外部区域的动态相量模型互联,
将动态相量表达形式的内部区域配电网状态空间方程和外部区域分布式电源状态空间方程以网络接口的形式对接,即内部区域状态空间方程的输入量为外部区域状态空间方程的输出量和外接电压源的电压量,外部区域状态空间方程的输入量为内部区域状态空间方程的输出量。
应用例
为验证本发明所述含分布式电源配电网简化建模方法的有效性和准确性,以图4所示的典型的IEEE 34节点配电网系统为研究对象,其中包含单相、两相、三相的线路与负荷,在编号850节点上接入分布式电源,编号800节点接入电压源。应用本发明所述的建模方法对含分布式电源接入的配电网进行仿真分析,并与电磁暂态仿真结果进行对比。
以大电源接入处的电流Is以及DG接入处的电流Idg为观察对象,图5a为DG功率指令为0.5MW时相应的仿真结果。从仿真结果可以看出,简化处理的模型与最原始模型的电磁暂态仿真结果完全吻合。说明提出的仿真模型能够准确地模拟分布式电源接入下的配电网系统。
为充分验证本发明所述简化模型的准确性和适用性,下面考虑几种变化情况:
(1)分布式电源功率的波动,假设分布式电源的功率在0.3s时从0.2MW上升为0.8MW,可以得到图5b的仿真结果。从仿真结果可以看出,本发明所述的模型可以很好地模拟系统的动态过程。图中所存在的偏差是因为采用的动态相量模型只提取了其中的基波频率量,其他的谐波频率量并不考虑,因此会造成一点偏差,但是很快两者就保持一致,并不影响系统的动态行为分析。
(2)故障情况下,假设在0.3s时节点828处发生三相接地短路故障。仿真结果如图5c所示,从结果可以看出,在故障开始发生时,本发明所述的动态相量模型能够很好地跟踪系统的动态行为,最后与原系统的稳定状态相同。
由以上仿真结果可以看出,对于多种仿真场景,本发明所述模型的仿真结果均较为准确,稳态误差很小且过渡过程基本吻合,能很好地反映系统的动态过程和主要特征。
仿真计算时间比较如表1所示,可以看出本发明所述的模型简化方法显著节省了仿真用时,在仿真速度上相对电磁暂态模型具有突出的优势。
表1模型所用仿真时间对比
可见,本发明所述的含分布式电源接入的配电网简化建模方法精度高,仿真速度快,为大规模分布式电源接入配电网的快速仿真奠定了基础。
Claims (1)
1.含分布式电源接入的配电网简化建模方法,其特征在于:
将含分布式电源的配电网分为两个区域:内部线性区域,即配电网网架部分;外部非线性区域,即分布式电源部分;
内部区域主要是由系统元件组成,包括线路、变压器,将负荷视作恒阻抗负荷;外部非线性分布式电源假设采用恒功率控制方式,以最大可能输出有功功率;先对内部线性的配电网网架区域推导以基本回路电流和电容电压为状态变量的配电网状态空间方程,并采用降阶方法实现线性系统的简化处理,在不改变系统的输入输出个数的情况下减少状态方程中的状态量数量,得到降阶后的配电网状态空间方程;对外部分布式电源区域建立微分代数方程,得到分布式电源的状态空间方程;并且用动态相量法提取内部线性区域以及外部非线性区域状态空间方程的主要频率特征量;最后以网络接口的形式对接,对整个系统进行互联,建立含分布式电源接入的配电网新仿真模型;
具体步骤如下:
步骤1),将含分布式电源接入的配电网划分成内部、外部区域;
将含分布式电源接入的配电网划分为两个区域,一个是以线性的配电网网架部分为对象的内部区域和以非线性的分布式电源为对象的外部区域;
步骤2),对内部线性区域建立配电网网架部分的状态空间方程;
利用基本回路矩阵Bb表述配电网络的基尔霍夫电流方程和基尔霍夫电压方程,如下式所示:
Bbvbr=0
其中ibr为支路电流向量,vbr为支路电压向量,ix为基本回路电流向量,基本回路矩阵Bb由系统初始的节点关联矩阵推导得到;
采用ASMG方法中提出的一般支路模型,ik和jk分别代表第k条支路的支路电流和电流源电流;对于系统中的所有支路将伏安关系写成向量形式:
vbr=Rbribr+Lbrpibr+Pbrqbr+ebr
其中,Pbr主对角线元素为相应支路电容的倒数,若该支路不存在电容则对应位置为0;p为微分算子;jbr为支路电流源电流向量;Lbr是支路电感矩阵,其主对角线元素表示相应支路的自感,非主对角线元素表示对应支路的互感;Rbr是支路电阻矩阵,ebr为支路电压源电压向量;
为了让输入系统的电压量和电流量与一般支路模型对应,定义相应的关联矩阵W1∈Rn _vs×n_br,W2∈Rn_cs×n_br,其中n_vs、n_cs分别表示电压、电流输入量的个数,n_br表示网络中支路个数;当第r个电压输入量位于支路k上时,则W1(r,k)=1,否则W1(r,k)=0;W2采用类似方式定义,即当第r个电流输入量位于支路k上时,则W2(r,k)=1,否则W2(r,k)=0;可得输入量与支路量的对应关系为:
其中ebr为支路电压源电压向量,jbr为支路电流源电流向量,vvs为输入电压向量,ics为输入电流向量;
对于配电网来说,分布式电源的输入量为电流量;将模型输入量定义为两种类型,即电压输入端子和电流输入端子;结合实际情况,定义电压输入端子电流ics与电流输入端子电压vvs为模型输出量,可以得到以基本回路电流ix和电容电压vc为状态变量的配电网状态空间方程如下式:
其中ix为基本回路电流向量,vc为支路电容电压向量,Cbr为支路电容矩阵,支路电容关联矩阵M∈Rn_br×n_c,n_c表示网络中电容个数,n_br表示网络中支路个数,取值方式与W1和W2相同;对照标准状态方程式:
y=LTx
可得标准状态方程对应的系数矩阵表达式如下:
步骤3),对配电网网架部分的状态空间方程降阶,求取其降阶后的状态空间方程,
在建立完内部配电网架部分的状态空间方程后,对模型进行降阶处理;以下为具体降阶方法:
假设某个线性时不变的动态系统:
y(t)=BTx(t)
计算得到该系统的拉普拉斯域下的传递函数如下:
H(s)=BT(s0E-A+(s-s0)E)-1B
=BT(I+(s-s0)M)-1R
其中
M:=(s0E-A)-1E,R:=(s0E-A)-1B
n阶Krylov子空间可以通过下面得到:
kn(M,R)=colspan[R MR M2R … MnR]
进而得到空间的正交特征向量Vn:
colspanVn=kn(M,R)
再利用SPRIM算法,获取特征矩阵
将Vn根据状态变量中的电压、电流量拆分为子向量和其中分别代表状态变量中的电压和电流量,利用和可以构造得到特征矩阵
则以基本回路电流和电容电压为状态变量的配电网状态空间方程可以变换为下面降阶后的状态空间方程:
其中,上式方程加下标n表示降阶后的状态空间方程,该方程矩阵的阶数小于降阶前的系统;上述状态方程对应的系数矩阵表达式如下:
步骤4),对外部非线性区域建立分布式电源的状态空间方程;
对于外部非线性区域的分布式电源,它由三个子模块组成:功率控制、电流控制以及网侧滤波器,分别建立其状态空间方程;
对于功率控制,假设分布式电源采用恒功率控制策略,即PQ控制,瞬时有功功率p和无功功率q由下式计算得到:
p=vodiod+voqioq
q=vodioq-voqiod
其中vod、iod、voq、ioq分别为电网电压、电流在两相同步旋转坐标系dq下的d轴电压、电流和q轴电压、电流;
功率控制利用PI调节,则对应的状态方程如下:
以及代数方程为:
其中,Pref和Qref分别是有功、无功功率指令,和是状态变量,和为电流控制的输入参考值,Kp和Ki分别是功率控制中的比例调节增益和积分调节增益;
则功率控制的状态空间方程可列写如下,方程的输入量为参考功率和反馈功率:
其中:
对于电流控制,对应状态方程为:
以及代数方程为:
其中,φd和φq是状态变量,和电流控制的输出参考值,ω为电网角频率,Lc为逆变器出口侧滤波电感的大小,Kpv和Kiv分别是电流控制中的比例调节增益和积分调节增益;
则电流控制的状态空间方程可列写如下,方程的输入量为参考电流和反馈电流:
其中:
对于网侧滤波器,在假设逆变器实际提供输出电压vid、viq等于指令电压的情况下,滤波器模型可列写如下:
则网侧滤波器的状态空间方程为:
其中:
[iodq]=[iod ioq]T,[vidq]=[vid viq]T,[vodq]=[vod voq]T
将上述功率控制、电流控制、网侧滤波器的状态空间方程整理合并,得到整个分布式电源的状态空间方程:
[idq]=Ct[xDG]
其中:
步骤5),对内部、外部区域的状态空间方程进行动态相量建模,
运用基于信号调制理论的动态相量模型,即对时域内周期为T的信号选择宽度为T的时间窗内的主要傅立叶系数,对配电网内部区域的以基本回路电流和电容电压为状态变量的配电网状态空间方程以及外部区域的分布式电源状态空间方程改写成动态相量的表达形式,分别计算如下:
<idq>k=Ct<xDG>k
其中,<·>k代表相关变量的k阶动态相量;
由于上述分布式电源的控制均是基于dq0坐标系,则将分布式电源的动态相量模型转换为abc坐标系以实现内部区域与外部区域的互联,转换关系满足以下方程:
其中,f代表可能的电量参数,例如电压或电流,DQ代表同步旋转坐标系,p、n分别代表正序和负序;
步骤6),内部、外部区域的动态相量模型互联,
将动态相量表达形式的内部区域配电网状态空间方程和外部区域分布式电源状态空间方程以网络接口的形式对接,即内部区域状态空间方程的输入量为外部区域状态空间方程的输出量和外接电压源的电压量,外部区域状态空间方程的输入量为内部区域状态空间方程的输出量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610371147.0A CN106021768B (zh) | 2016-05-30 | 2016-05-30 | 含分布式电源接入的配电网简化建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610371147.0A CN106021768B (zh) | 2016-05-30 | 2016-05-30 | 含分布式电源接入的配电网简化建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106021768A CN106021768A (zh) | 2016-10-12 |
CN106021768B true CN106021768B (zh) | 2019-01-11 |
Family
ID=57091551
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610371147.0A Expired - Fee Related CN106021768B (zh) | 2016-05-30 | 2016-05-30 | 含分布式电源接入的配电网简化建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106021768B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108347060A (zh) * | 2017-01-22 | 2018-07-31 | 中国电力科学研究院 | 一种电力电子化接口电源降阶建模方法及系统 |
CN107204614A (zh) * | 2017-05-31 | 2017-09-26 | 中南大学 | 一种包含多并联dc‑dc转换器的直流微电网系统的稳定方法 |
CN107944179A (zh) * | 2017-12-07 | 2018-04-20 | 中国航空工业集团公司西安航空计算技术研究所 | 一种机载计算机电源系统稳定性设计方法 |
CN108448631B (zh) * | 2018-04-25 | 2020-10-27 | 浙江大学 | 基于模型降阶的含分布式电源接入配电网的动态相量建模方法 |
CN108490379B (zh) * | 2018-05-18 | 2020-04-10 | 云南电网有限责任公司电力科学研究院 | 一种基于自激振荡波的变压器绕组波过程校验方法 |
CN110492531B (zh) * | 2019-08-06 | 2020-09-11 | 清华大学 | 考虑同步旋转惯量水平的电力系统调度运行方法及系统 |
CN110994666B (zh) * | 2019-10-21 | 2022-09-27 | 合肥工业大学 | 在谐振频率处将阻抗模型降阶为rlc电路的降阶方法 |
CN112395721A (zh) * | 2020-12-09 | 2021-02-23 | 南京四方亿能电力自动化有限公司 | 一种适用于配电网的快速绘图建模方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012012897A1 (fr) * | 2010-07-28 | 2012-02-02 | HYDRO-QUéBEC | Localisation d'un défaut sur une section de ligne électrique hors tension |
CN103049617A (zh) * | 2012-12-27 | 2013-04-17 | 天津大学 | 保留无源性的大规模配电网络电磁暂态仿真模型化简方法 |
CN103870649A (zh) * | 2014-03-19 | 2014-06-18 | 国家电网公司 | 一种基于分布式智能计算的主动配电网自治化仿真方法 |
CN104283217A (zh) * | 2014-11-03 | 2015-01-14 | 广东电网有限责任公司佛山供电局 | 含分布式电源与微电网的配电网运行区间判定方法 |
-
2016
- 2016-05-30 CN CN201610371147.0A patent/CN106021768B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012012897A1 (fr) * | 2010-07-28 | 2012-02-02 | HYDRO-QUéBEC | Localisation d'un défaut sur une section de ligne électrique hors tension |
CN103049617A (zh) * | 2012-12-27 | 2013-04-17 | 天津大学 | 保留无源性的大规模配电网络电磁暂态仿真模型化简方法 |
CN103870649A (zh) * | 2014-03-19 | 2014-06-18 | 国家电网公司 | 一种基于分布式智能计算的主动配电网自治化仿真方法 |
CN104283217A (zh) * | 2014-11-03 | 2015-01-14 | 广东电网有限责任公司佛山供电局 | 含分布式电源与微电网的配电网运行区间判定方法 |
Non-Patent Citations (1)
Title |
---|
基于状态变量分析的有源配电网电磁暂态仿真自动建模方法;于浩 等;《电网技术》;20150630;第39卷(第6期);1518-1524 |
Also Published As
Publication number | Publication date |
---|---|
CN106021768A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106021768B (zh) | 含分布式电源接入的配电网简化建模方法 | |
Ma et al. | Modeling and transient stability analysis for type-3 wind turbines using singular perturbation and Lyapunov methods | |
Ajjarapu | Computational techniques for voltage stability assessment and control | |
Groß et al. | On the steady-state behavior of a nonlinear power system model | |
Zhou et al. | Small-signal stability assessment of heterogeneous grid-following converter power systems based on grid strength analysis | |
Wang et al. | Reduced-order state space model for dynamic phasors in active distribution networks | |
Yuan et al. | Assessing maximal capacity of grid-following converters with grid strength constraints | |
Davari et al. | A fault-tolerant, passivity-based controller enhanced by the equilibrium-to-equilibrium maneuver capability for the DC-voltage power port VSC in multi-infeed AC/DC modernized grids | |
Zhou et al. | DQ impedance‐decoupled network model‐based stability analysis of offshore wind power plant under weak grid conditions | |
Gao et al. | Shifted frequency‐based electromagnetic transient simulation for AC power systems in symmetrical component domain | |
Meng et al. | Stability analysis and control optimisation based on particle swarm algorithm of modular multilevel matrix converter in fractional frequency transmission system | |
CN105305392A (zh) | 适用于含电压控型iidg配电网的短路计算对称分量法 | |
Dubyago | Mathematical description of interaction of the elements in the power network system | |
Guo et al. | Single‐Phase Reactive Power Compensation Control for STATCOMs via Unknown System Dynamics Estimation | |
Nicolini et al. | PLL parameters tuning guidelines to increase stability margins in multiple three‐phase converters connected to weak grids | |
CN106505612A (zh) | 基于mmc的大功率永磁同步风力发电并网系统控制方法 | |
Chen et al. | Error estimation method of reduced-order small-signal model for multiterminal DC distribution network | |
CN115905781A (zh) | 一种双馈风电场等值方法、系统及设备 | |
Ebrahimi et al. | Numerically Efficient Average-Value Model for Voltage-Source Converters in Nodal-Based Programs | |
Levron et al. | Observable canonical forms of multi-machine power systems using dq0 signals | |
Liu et al. | Dynamic aggregation modeling for droop control inverter based on slow coherency algorithm | |
Campos et al. | Novel Gramian-based Structure-preserving Model Order Reduction for Power Systems with High Penetration of Power Converters | |
Davari et al. | Analysing dynamics and synthesising a robust vector control for the dc‐voltage power port based on the modular multilevel converter in multi‐infeed AC/DC smart grids | |
Wang et al. | Krylov subspace based model reduction method for transient simulation of active distribution grid | |
Rathore et al. | ARI and ARID control of virtual synchronous generator for frequency response improvement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190111 |
|
CF01 | Termination of patent right due to non-payment of annual fee |