CN102354332A - 一种用于简化柔性交直流输电系统中rga计算的方法 - Google Patents
一种用于简化柔性交直流输电系统中rga计算的方法 Download PDFInfo
- Publication number
- CN102354332A CN102354332A CN2011102767550A CN201110276755A CN102354332A CN 102354332 A CN102354332 A CN 102354332A CN 2011102767550 A CN2011102767550 A CN 2011102767550A CN 201110276755 A CN201110276755 A CN 201110276755A CN 102354332 A CN102354332 A CN 102354332A
- Authority
- CN
- China
- Prior art keywords
- delta
- prime
- generator
- rga
- omega
- 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
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种用于简化柔性交直流输电系统中RGA计算的方法,通过先对柔性交直流输电系统模型中的微分代数方程组进行线性化,再通过简明的矩阵运算得到传递函数,再在此基础上进行RGA计算。虽然形成的线性化系数矩阵元素较多,但避免了多次矩阵求逆导致的矩阵元素表达式过于复杂的问题,易于利用程序实现,提高了运算速度,降低了技术成本和对运行环境的要求,适用于高维数柔性交直流输电系统。
Description
技术领域
本发明涉及一种柔性交直流输电系统中RGA计算的方法,尤其涉及的是一种用于简化柔性交直流输电系统中RGA计算的方法。
背景技术
柔性交流输电系统即将现代电力电子技术和现代自动控制技术引入到交直流输电系统中,以提高原有交直流输电系统的可控性、灵活性、运行的稳定性和经济性。该现代电力电子技术和现代自动控制技术即为柔性交流输电技术,即在传统的输电系统中安装FACTS(柔性交流输电)装置。然而,研究表明FACTS装置的不同控制回路之间存在着交互影响,这种交互影响对其本身的控制效果和系统的稳定性都有着较大的负面作用。因此,需要用到相对增益矩阵(RGA)方法对FACTS多个控制回路之间的交互影响进行量化分析,从而采取相应的技术手段以减少其负面作用。现有技术通常是对该柔性交直流输电系统模型的平衡方程进行节点消去从而得到传递函数,再在此基础上对柔性交直流输电系统模型进行RGA计算。该方法容易多次矩阵求逆导致的矩阵元素表达式过于复杂的问题,用程序实现比较复杂,计算成本高,尤其不适用于高维数电网模型的RGA计算。
因此,现有技术还有待于改进和发展。
发明内容
本发明的目的在于提供一种用于简化柔性交直流输电系统中RGA计算的方法,以解决现有技术中对该柔性交直流输电系统模型的平衡方程进行 节点消去从而得到传递函数,再在此基础上对柔性交直流输电系统模型进行RGA计算,从而导致的计算成本高,实现程序复杂,不适于高维数电网模型RGA计算的问题。
本发明的技术方案如下:
一种用于简化柔性交直流输电系统RGA计算的方法,其包括以下步骤:
步骤A:构建柔性交直流输电系统数学模型:
g(x,y,u)=0
其中式 为描述系统各元件动态的微分方程,包括发电机及其励磁系统和FACTS装置的动态;式g(x,y,u)=0为表示网络各个节点电压电流关系;本步骤主要目的是将柔性交直流输电系统用数学模型表示出来,包含各发电机及其励磁系统、FACTS装置的动态特性,各个节点的电压、电流之间的关系。
步骤B:将描述发电机和励磁系统动态的微分方程写成直角坐标形式并进行线性化,得电力系统的线性化模型:
本步骤主要解决了高维数微分方程的求解问题,将其变化为线性化模型,从而可以将该方法应用于含有柔性交直流输电装置的大电网的此类问题的求解。
步骤C:对步骤B中的电力系统的线性化模型再次进行线性化可得:
其中,Yphzt为2n×4m阶矩阵,Yphds为2n×2n阶矩阵,Yphkz为2n×2阶矩阵;
本步骤主要是线性化模型做技术处理,提高求解效率。
步骤D:对用于RGA计算的传递函数进行线性化,可得:
式中,Yscsc为2×2阶矩阵,Ysczt为2×4m阶矩阵,Yscds为2×2n阶矩阵,Ysckz为2×2阶矩阵。
本步骤主要是对用于RGA计算的传递函数进行线性化,解决传递函数高维数,难以求解的问题。
步骤E:联合步骤B、步骤C和步骤D中的三个线性化模型得出传递函数:
其中,A=Ywfzt-Ywfds*(Yphds-Yphzt);B=Ywfkz-Ywfds*(Yphds-Yphkz);C=(Yscsc -1Yscds)*(Yphds -1Yphzt)-Yscsc -1Ysczt;D=(Yscsc -1Yscds)*(Yphds -1Yphkz)。
本步骤将上述三个线性化模型综合变换,得出传递函数。利用本传递函数即可求解并分析柔性交直流输电系统中电网各元件相互影响的问题。
步骤F:将步骤E得到的传递函数转化为频域传递函数:
本步骤主要是将传递函数转化为频域传递函数,将其变换为RGA求解的形式,利用RGA算法求解。
所述的用于简化柔性交直流输电系统RGA计算的方法,其中,发电机系统采用的是三阶实用模型:
其中,δ为发电机功角;w为角速度;Pm为原动机输入的机械功率; 为发电机的暂态电势; 为发电机的暂态电抗;Xd为发电机的d轴同步电抗;Xq为发电机的q轴同步电抗;Tj为发电机的惯性时间常数; 为发电机的直轴暂态开路时间常数;Efd为励磁电压;D为阻尼系数。
所述的用于简化柔性交直流输电系统RGA计算的方法,其中,励磁系统的数学模型为:
本发明通过先对柔性交直流输电系统模型中的微分代数方程组进行线性化,再通过简明的矩阵运算得到传递函数,再在此基础上进行RGA计算。虽然形成的线性化系数矩阵元素较多,但避免了多次矩阵求逆导致的矩阵元素表达式过于复杂的问题,易于利用程序实现,提高了运算速度,降低了技术成本和对运行环境的要求,适用于高维数柔性交直流输电系统。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚、明确,以下举实施例对本发明进一步详细说明。
本发明中方法的基本流程如下:
步骤A:构建柔性交直流输电系统的数学模型:
g(x,y,u)=0…………(2)
式(1)为描述系统各元件动态的微分方程,包括发电机及其励磁系统和FACTS装置的动态;式(2)为描述网络各个节点电压电流关系的代数方程。
本步骤主要目的是将柔性交直流输电系统用数学模型表示出来,包含各发电机及其励磁系统、FACTS装置的动态特性,各个节点的电压、电流之间的关系。
柔性交直流输电系统动态过程可用式(1)和(2)的模型描述,x为状态变量,即节点电压的幅值和相角;y为代数变量,一般是网络中节点电压向量;u为控制变量,包括发电机有功、无功输出功率、发电机机端电压和变压器变比等。
本发明提供的发电机系统采用的是三阶实用模型,如式(3)~(5)所表示的:
其中,δ为发电机功角,w为角速度,Pm为原动机输入的机械功率, 为发电机的暂态电势、 为发电机的暂态电抗、Xd为发电机的d轴同步电抗、Xq为发电机的q轴同步电抗、Tj为发电机的惯性时间常数、 为发电机的直轴暂态开路时间常数;Efd为励磁电压;D为阻尼系数。
其中,id和iq可由下式求出:
id、iq分别表示d轴、q轴的电流(A),通过坐标变换将电力系统从abc三相坐标变换到dq0坐标,参数ua,ub,uc,ia,ib,ic,磁链a,磁链b,磁链c这些量都变换到dq0坐标中,如果有需要可以逆变换回来。从物理意义上讲,该变换就是将ia,ib,ic电流投影等效到d轴和q轴上,将定子上的电流都等效到直轴和交轴上去。对于稳态来说,等效之后,iq,id正好就是一个常数。
构建励磁系统的数学模型,如式(7):
步骤B:假设多级系统中共有n个节点,m台发电机。对式(1)中描述发电机和励磁系统动态的微分方程,即发电机系统的三阶实用模型、以及励磁系统的数学模型写成直角坐标形式并进行线性化,可得到:
其中,Δδ[Δδ1Δδ2…Δδm]T;Δω=[Δω1Δω2…Δωm]T;ΔE′q=[ΔE′q1ΔE′q2…ΔE′qm]T;ΔEfd=[ΔEfd1ΔEfd2…ΔEfdm]T;Vi=[Vix Viy]T;Ywfzt为4m×4m阶矩阵,Ywfds为4m×2n阶矩阵,Ywfkz为4m×2阶矩阵。
将描述网络各个节点电压电流关系的代数方程式(2),改写成为:
I-(Y0+Y′)V=0 (9)
其中,I=[Ix1Iy1…IxiIyi…IxnIyn]T,为节点注入电流向量;V=[Vx1Vy1…VxiVyi…VxnVyn]T,为节点电压向量;Y0为网络节点导纳矩阵,Y′=diag(Y1′,…,Yi′,…,Yn′),是由发电机、负荷和FACTS装置并入网络的导纳所构成的矩阵。
若节点i为发电机的节点,则其导纳为:
式(10)中,Gxi、Bxi、Byi、Gyi为发电机i及其凸极效应后并入网络的导纳,其中Gxi、Bxi、Byi、Gyi通过式(11)求得:
若节点i安装STATCOM,则该节点的注入电流为:
I′=I+IST (12)
其中IST可以表示为:
若节点i安装SVC,则该节点的自导纳为:
本步骤主要解决了高维数微分方程的求解问题,将其变化为线性化模型,从而可以将该方法应用于含有柔性交直流输电装置的大电网的此类问题的求解。
步骤C:式(8)是含状态变量、代数变量和控制变量的代数方程组,对其线性化可得:
式中,Yphzt为2n×4m阶矩阵,Yphds为2n×2n阶矩阵,Yphkz为2n×2阶矩阵。
由式(8)、(15)构成了含FACTS装置的多机电力系统线性化模型。用于RGA计算的传递函数的求解方法的系统输出为:
Y=g1(x,y,u) (16)
本步骤主要是线性化模型做技术处理,提高求解效率。
步骤D:对式(16)线性化,可得:
式中,Yscsc为2×2阶矩阵,Ysczt为2×4m阶矩阵,Yscds为2×2n阶矩阵,Ysckz为2×2阶矩阵。本步骤主要是对用于RGA计算的传递函数进行线性化,解决传递函数高维数,难以求解的问题。
步骤E:联立式(8)、(15)和(17),消去代数变量Δy,得到:
其中,A=Ywfzt-Ywfds*(Yphds -1Yphzt);B=Ywfkz-Ywfds*(Yphds -1Yphkz);C=(Yscsc -1Yscds)*(Yphds -1Yphzt)-Yscsc -1Ysczt;D=(Yscsc -1Yscds)*(Yphds -1Yphkz)。
本步骤将上述三个线性化模型综合变换,得出传递函数。利用本传递函数即可求解并分析柔性交直流输电系统中电网各元件相互影响的问题。
步骤F:将式(18)转化为频域传递函数的形式:
将上述频域传递函数用于RGA计算中,以简化RGA计算。
本步骤主要是将传递函数转化为频域传递函数,将其变换为RGA求解的形式,利用RGA算法求解。
本发明提出的线性化模型,未对平衡方程进行节点消去,而是先对微分代数方程组进行线性化,再通过简明的矩阵运算得到传递函数。虽然形成的线性化系数矩阵元素较多,但避免了多次矩阵求逆导致的矩阵元素表达式过于复杂的问题,易于利用程序实现,适用于多机系统。
应当理解的是,本发明的应用不限于上述的举例,对本领域普通技术 人员来说,可以根据上述说明加以改进或变换,所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (3)
1.一种用于简化柔性交直流输电系统RGA计算的方法,其包括以下步骤:
步骤A:构建柔性交直流输电系统数学模型:
g(x,y,u)=0
其中式为描述系统各元件动态的微分方程,包括发电机及其励磁系统和FACTS装置的动态;式g(x,y,u)=0为表示网络各个节点电压电流关系;
步骤B:将描述发电机和励磁系统动态的微分方程写成直角坐标形式并进行线性化,得电力系统的线性化模型:
步骤C:对步骤B中的电力系统的线性化模型再次进行线性化可得:
其中,Yphzt为2n×4m阶矩阵,Yphds为2n×2n阶矩阵,Yphkz为2n×2阶矩阵;
步骤D:对用于RGA计算的传递函数进行线性化,可得:
式中,Yscsc为2×2阶矩阵,Ysczt为2×4m阶矩阵,Yscds为2×2n阶矩阵,Ysckz为2×2阶矩阵。
步骤E:联合步骤B、步骤C和步骤D中的三个线性化模型得出传递函数:
其中,A=Ywfzt-Ywfds*(Yphds -1Yphzt);B=Ywfkz-Ywfds*(Yphds -1Yphkz)C=(Yscsc -1Yscds)*(Yphds -1Yphzt)-Yscsc -1Ysczt;D=(Yscsc -1Yscds)*(Yphds -1Yphkz)。
步骤F:将步骤E得到的传递函数转化为频域传递函数:
3.根据权利要求1所述的用于简化柔性交直流输电系统RGA计算的方法,其特征在于,励磁系统的数学模型为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110276755.0A CN102354332B (zh) | 2011-09-16 | 2011-09-16 | 一种用于简化柔性交直流输电系统中rga计算的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110276755.0A CN102354332B (zh) | 2011-09-16 | 2011-09-16 | 一种用于简化柔性交直流输电系统中rga计算的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102354332A true CN102354332A (zh) | 2012-02-15 |
CN102354332B CN102354332B (zh) | 2014-12-24 |
Family
ID=45577897
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110276755.0A Active CN102354332B (zh) | 2011-09-16 | 2011-09-16 | 一种用于简化柔性交直流输电系统中rga计算的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102354332B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102611128A (zh) * | 2012-03-09 | 2012-07-25 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种高压直流输电系统直流功率调制方法 |
CN103018583A (zh) * | 2012-08-03 | 2013-04-03 | 中国能源建设集团广东省电力设计研究院 | 基于mmc柔性直流输电系统电平数选择验证方法 |
CN104268798A (zh) * | 2014-09-25 | 2015-01-07 | 国家电网公司 | 一种分布式电源及电动汽车对配电网影响的评估方法 |
CN109038616A (zh) * | 2018-07-27 | 2018-12-18 | 中国电力科学研究院有限公司 | 一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630273A (zh) * | 2009-08-06 | 2010-01-20 | 中国电力科学研究院 | 一种电力系统小干扰稳定仿真方法 |
CN101917003A (zh) * | 2010-07-29 | 2010-12-15 | 东南大学 | 一种电力系统小干扰振荡稳定性分解式模态分析方法 |
-
2011
- 2011-09-16 CN CN201110276755.0A patent/CN102354332B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630273A (zh) * | 2009-08-06 | 2010-01-20 | 中国电力科学研究院 | 一种电力系统小干扰稳定仿真方法 |
CN101917003A (zh) * | 2010-07-29 | 2010-12-15 | 东南大学 | 一种电力系统小干扰振荡稳定性分解式模态分析方法 |
Non-Patent Citations (2)
Title |
---|
江全元 等: "基于相对增益矩阵原理的柔性交流输电系统控制器交互影响分析", 《中国电机工程学报》 * |
董萍 等: "FACTS装置中多个控制器之间的交互影响分析", 《湖南工业大学学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102611128A (zh) * | 2012-03-09 | 2012-07-25 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种高压直流输电系统直流功率调制方法 |
CN102611128B (zh) * | 2012-03-09 | 2014-03-19 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种高压直流输电系统直流功率调制方法 |
CN103018583A (zh) * | 2012-08-03 | 2013-04-03 | 中国能源建设集团广东省电力设计研究院 | 基于mmc柔性直流输电系统电平数选择验证方法 |
CN103018583B (zh) * | 2012-08-03 | 2016-02-03 | 中国能源建设集团广东省电力设计研究院有限公司 | 基于mmc柔性直流输电系统电平数选择验证方法 |
CN104268798A (zh) * | 2014-09-25 | 2015-01-07 | 国家电网公司 | 一种分布式电源及电动汽车对配电网影响的评估方法 |
CN104268798B (zh) * | 2014-09-25 | 2018-01-16 | 国家电网公司 | 一种分布式电源及电动汽车对配电网影响的评估方法 |
CN109038616A (zh) * | 2018-07-27 | 2018-12-18 | 中国电力科学研究院有限公司 | 一种计算励磁附加调差对同步发电机动态阻尼影响的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102354332B (zh) | 2014-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104866665B (zh) | 基于接口等值与交互的含电力电子设备的混合仿真方法 | |
Xia et al. | Direct nonlinear primal–dual interior-point method for transient stability constrained optimal power flow | |
Leon et al. | Multi-machine power system stability improvement using an observer-based nonlinear controller | |
CN103995734A (zh) | 基于rtds的电力系统混合实时仿真系统及仿真方法 | |
CN105048917A (zh) | 基于eso的双馈风力发电系统积分滑模控制器的控制方法 | |
Groß et al. | On the steady-state behavior of a nonlinear power system model | |
CN103401503A (zh) | 一种在谐波平面在线辨识双三相电机参数的方法及装置 | |
CN102427229B (zh) | 基于修正牛顿法的带零注入约束的电力系统状态估计方法 | |
CN110768299A (zh) | 负荷虚拟同步机的序阻抗建模与稳定性分析方法 | |
CN104901394B (zh) | 基于soc的光储式充电站准pr下垂控制方法 | |
Belikov et al. | Comparison of time-varying phasor and dq0 dynamic models for large transmission networks | |
CN102354332A (zh) | 一种用于简化柔性交直流输电系统中rga计算的方法 | |
CN107733273A (zh) | 基于分序和虚拟同步发电机复合控制的三相四桥臂逆变器的控制方法和系统 | |
CN106533289B (zh) | 一种非线性电压控制方法及系统 | |
CN104375876B (zh) | 一种输入量突变情况下的0+误差免疫电磁暂态仿真方法 | |
CN104377697A (zh) | 模拟同步发电机动态特性的电压源换流站的控制方法 | |
CN106549402A (zh) | Tcsc抑制交直流并联运行系统次同步振荡的方法 | |
CN108573094B (zh) | 同步发电机的vbr电磁暂态仿真模型的建立方法及系统 | |
CN103217896B (zh) | 基于自由权矩阵方法的多facts抗时滞协调控制方法 | |
CN115470736A (zh) | 适应储能电站变工况运行的电力系统动态行为建模方法 | |
CN104701881A (zh) | 一种基于凸多面体理论的hvdc系统广域分散协调控制方法 | |
CN110474322B (zh) | 包含虚拟同步发电机的多机系统小干扰稳定性评估方法 | |
CN204229391U (zh) | 混合仿真接口及具有该接口的电力系统混合实时仿真设备 | |
CN103973151B (zh) | 电感不平衡条件下三相pwm并网逆变器的解耦控制方法 | |
Tang et al. | Full electromagnetic transient simulation for large power systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |