基于Matlab的柴油发电机组功率
交互振荡现象的仿真
闫利伟 杜海生 宋东升 刘冀
（海军工程大学天津校区，天津塘沽300450）
摘 要：交流供电使用的柴油发电机组，并联运行过程中经常出现功率交互振荡。本文先通过分析产生功率振荡的原因，得出影响功率振荡的参数，然后利用Matlab对柴油发电机组单机、并联运行进行仿真，验证分析结果。
关键词：柴油发电机组 并联运行 功率交互振荡
中图分类号：TM314 文献标识码：A 文章编号：1003-4862 (2010) 01-0021-06
The Simulation Of Interactive Power Oscillation For Diesel
Generator Sets Based On Matlab
Yan Liwei, Du Haiheng, Song Dongsheng, Liu Ji
（Naval University of Engineering Tianjin Campus, Tanggu 300450, Tianjin, China）
Abstract: Interactive power oscillation often appears in diesel generator sets for AC power supply operating in parallel. This paper analyzes the reason for power oscillation, and its the parameter. Besides that, it makes simulation for the stand-alone or paralleling connection of diesel generator by using Matlab, and verifies the result.
Key words: diesel generator sets, operation in parallel, interactive power oscillation
干扰力矩产生的原因如下: 1 引言
⑴ 各缸转矩不均匀所引起的干扰力矩
交流供电使用的柴油发电机组广泛地应用于
各缸转矩不均匀原因有：
船舶、通信等国民经济的众多领域。交流柴油发
a. 高压燃油泵和喷油嘴调节不均匀，各缸的
电机组在互相并联运行及并入电网运行时，往往
喷油大小量不等；
出现功率的周期性振荡。功率振荡时，并联运行
b. 喷油泵凸轮传动扭转变形；
的各机组的总负载不变。因此，一台机组的功率
c. 配气相位、喷油提前角、压缩比等参数调
增加和另一台机组的功率减少同时进行。通常把
整不一致；
这种现象称为功率交互振荡，俗称“游车”。无论
d. 各缸制造误差不一致。
在机组台架试验中，还是在实际运行的舰船或陆
⑵ 高压燃油管剩余压力波动引起的干扰力
上电站中都曾发现功率交互振荡的现象。
矩
X.甘切夫提出把喷油泵齿条位置不变下喷油2 柴油发电机组功率交互振荡的分析[1]
后输油管内的残余压力pocm2与喷油前输油管内
产生功率交互振荡的原因是由于并联运行的
残余压力pocm1的关系作为评价供油装置工作稳
机组中有干扰力矩的存在。
定性的标准[2]。 可以近似认为，齿条位置一定时，每循环实收稿日期: 2009-10-15
际供油量的变化（相对单位）为： 作者简介：闫利伟（1978－），男，天津校区教员，研
究方向为电力系统运行、控制与分析。
∆gu=−m0(∆pocm2−∆pocm1) (1)
21
船电技术|电机 Vol.30 No.1 2010.1
式中：∆gu—每循环实际供油量的变化（相对单 位）；m0—表征高压输油管容积大小的比例系数。
在供油装置工作正常情况下，pocm2＝pocm1，即∆gu＝ 0。实际上影响pocm2有很多因素，在供油装置工作不正常情况下，就更容易造成∆gu≠ 0。尽管喷油泵齿条处于同一位置上，进入气缸的供油量却在周期性的变化，就会产生干扰力矩。影响pocm2波动的主要因素有：
a. 燃油的可压缩性；
b. 高压油管的容积变化。
柴油发电机组在出厂前都要做供油装置调整，使各项要求都达到国家标准。但是，往往在使用一定时期后，由于种种原因，产生供油不稳定。上面几种原因只要发生一种，气缸供给的燃油量就不相同，就会造成各个气缸燃烧状况不一致，产生干扰转矩[2-5]。产生干扰力矩的原因都是因为柴油机供油出现周期变化。本文将以上各种因素对供油量产生的干扰近似的用一正弦函数代替即：
△gr＝k sin(ωt+α) (2) 式中：△gr—为干扰油量；
k—为干扰油量的幅值（相对单位）； ω—为干扰油量的角频率； α—为干扰油量的初相角。
单机运行时，供油干扰虽然使柴油机输出的机械功率变化，电磁功率相应也变化，但由于负
载不变，稳定运行点是不变的，电磁功率围绕这个点上下波动，但偏差不是很大，所以供油干扰对单机运行影响不大。在多机并联运行时，并联机组总的负载是不变的，但加在每台机组上的负荷是变化的，这样就造成电磁功率稳定运行点不断变化，又因为所加干扰为周期性干扰，每台机组调速器在某一段时间起到抑制干扰的作用，而在另一段时间起到加强干扰的作用，所以使得机组电磁功率也在周期性变化。机组电磁功率变化的程度取决于干扰油量中三个参数k、ω、α。
3 柴油发电机组并联运行功率交互振荡
仿真
柴油发电机组单机运行模型如图1所示，双机并联运行模型如图2所示，两台机组为同型号机组，参数都一样。励磁模块[6]如图3所示，调速和柴油机模块[7,8]如图4所示，加入干扰的调速原动机模块如图5所示。
柴油发电机组型号为500GF，发电机型号为TFW－500－6，柴油机型号为P28V190。发电机参数如下：Pn＝500 kW，Un＝383 V，fn＝50 Hz，
X
\"'＝2.12，Xq＝1.35，Xd＝0.163，Xd＝0.13，
\"
＝0.11，Xl＝0.0572，Td＝0.121，Td＝0.01，Xq
'\"
Tq\"＝0.01，Rs＝0.003，H(s)＝2.1。
输电线路比较短，所以忽略不计。 仿真结果如下。
22 Vol.30 No.1 2010.1 船电技术|电机 1 图3 励磁系统模块内部模型图 图4 调速器模块内部模型图 图5 加入供油干扰后的调速器模块内部模型图 3.1 单机运行仿真 负荷参数：把所有负荷等效成一个电阻和电感并联的负载，Vn＝380 V，fn＝50 Hz，PL＝3e5 W，QL＝1000 var。 说明：仿真除被考察的功率值采用有名值外，其他量都采用标幺值。所考察的功率值指6 s后功率值，机组启动过程功率值不在考察范围内。 （1） 机组供油正常时。 （2）如图5，加入干扰油量△gr＝ksin(ωt+α)，其中k＝0.02（相对单位），ω＝9π，α＝0。 由图6、图7可知，柴油发电机组单机运行时，负载上消耗的功率近似于机组发出的功率，供油正常时机组发出的功率稳定，正常运行，当加入干扰油量时，机组发出的功率出现微小等幅周期振荡，振荡幅值为0.024×105 W。 3.2 双机并联运行仿真 负荷参数：把所有负荷等效成一个电阻和电感并联的负载，Vn＝380 V，fn＝50 Hz，PL＝6 23 船电技术|电机 Vol.30 No.1 2010.1 ×105 W，QL＝2000 var。 图6 供油正常时负载上消耗的功率 图7 供油有干扰时负载上消耗的功率 图8 功率变化曲线 （1）一、二号机组供油正常时。 由图8可以看出一、二号机组发出的功率一致、稳定，负荷消耗功率稳定，表明并联机组运行正常。 （2）一号机组加入干扰油量△gr＝ksin(ωt+α)，其中k1＝0.02（相对单位），ω1＝9π，α1＝0；二号机组供油正常。 24 图9 功率变化曲线 由图9可以看出，一、二号机组发出的功率发生等幅周期振荡。一号机组振荡幅值为0.121×105 W，二号机组振荡幅值为0.096×105 W。并且一台机组发出功率增加的同时，另一台机组发出功率在减小。负载上消耗的功率出现微小的等幅周期振荡，振荡幅值为0.04×105 W，并联机组发生功率交互振荡。 （3）一号机组加入干扰油量△gr＝ksin(ωt+α)，其中k1＝0.03（相对单位），ω1＝9π，α1＝0；二号机组供油正常。 图10 功率变化曲线 由图10可以看出，一、二号机组发出的功率发生等幅周期振荡。一号机组振荡幅值为0.181×105 W，二号机组振荡幅值为0.144×105 W。并且一台机组发出功率增加的同时，另一台机组发出功率在减小。负载上消耗的功率出现微小的等幅周期振荡，振荡幅值为0.038×105 W，并联机组发生功率交互振荡。与图9相比，一、二号机组功率振荡幅值增大。 （4）一、二号机组同时加入干扰油量△gr＝ksin(ωt+α)，其中k1＝k2＝0.02（相对单位），ω1＝ω2＝9π，α1＝α2＝0。 Vol.30 No.1 2010.1 船电技术|电机 图11 功率变化曲线 由图11可以看出，一、二号机组功率变化曲线一致，虽然有微小等幅周期振荡，但没有发生功率交互振荡。 （5）一、二号机组同时加入干扰油量△gr＝ksin(ωt+α)，其中k1＝0.02（相对单位），ω1＝9π，α1＝0，k2＝0.02（相对单位），ω2＝9π，α2＝π/2。 图12 功率变化曲线 由图12可以看出，一、二号机组发出的功率 发生等幅周期振荡。一号机组振荡幅值为0.155×105 W，二号机组振荡幅值为0.155×105 W。并且一台机组发出功率增加的同时，另一台机组发出功率在减小。负载上消耗的功率出现微小的等幅周期振荡，振荡幅值为0.035×105 W，并联机组发生功率交互振荡。 （6）一、二号机组同时加入干扰油量△gr＝ksin(ωt+α)，其中k1＝0.02（相对单位），ω1＝9π，α1＝0，k2＝0.02（相对单位），ω2＝9π，α2＝π。 由图13可以看出，一、二号机组发出的功率发生等幅周期振荡。一号机组振荡幅值为0.216×105W，二号机组振荡幅值为0.216×105W。 并且一台机组发出功率增加的同时，另一台机组发出功率在减小。负载上消耗的功率出现微小的等幅周期振荡，振荡幅值为0.001×105 W，并联机组发生功率交互振荡。与图12相比，一、二号机组功率振荡幅值增大。 图13 功率变化曲线 （7）一、二号机组同时加入干扰油量△gr＝ ksin(ωt+α)，其中k1＝0.02（相对单位），ω1＝9π，α1＝0，k2＝0.02（相对单位），ω2＝7π，α2＝0。 图14 功率变化曲线 由图14可以看出，一、二号机组发出的功率发生振荡。一号机组幅值变化范围为2.793×105～3.174×105 W，二号机组幅值变化范围为2.862×105～3.229×105 W。并且一台机组发出功率增加的同时，另一台机组发出功率在减小。负载上消耗的功率出现微小的振荡，幅值变化范围为5.961×105～6.046×105 W，并联机组发生功率交互振荡。与前几种振荡情况相比，ω不同产生的振荡是非周期变化的，幅值也不是定值。 （8）一、二号机组同时加入干扰油量△gr＝ksin(ωt+α)，其中k1＝0.02（相对单位），ω1＝9π， 25 船电技术|电机 Vol.30 No.1 2010.1 α1＝0，k2＝0.02（相对单位），ω2＝6π，α2＝0。 响功率交互振荡的三个关键参数，因此通过对幅 值k、角频率ω、初相角α三个量进行调节，可以有效消除柴油发电机组在并联运行时产生的功率交互振荡，提高并联运行的可靠性。 参考文献： [1] 闫利伟 宋东升, 常鲜戎. 柴油发电机组功率交 互振荡现象的分析与仿真研究. 电气应用, 2009,(8). [2] [苏]В.И.托尔申著. 柴油发电机组并联运行的稳 定性（田治喜, 吴霞芳译）. 北京: 国防工业出 图15 功率变化曲线 版社, 1977. [3] 苏石川, 元广杰, 杨宗明编著.现代柴油发电机 组的应用与管理. 北京: 化学工业出版社，2005.8. [4] 高国权编著.电站用柴油机调速系统.人民交通出 版社, 1983.12. [5] 项国波编著.柴油机交流发电机组并联运行及其 稳定性. 北京: 国防工业出版社，1979.5 [6] 任伟.船舶电站柴油发电机组自整定PID控制与 虚拟仿真：[硕士学士论文].上海：上海海事大学，2004. [7] 韦韩英. 船舶柴油发电机系统仿真：[硕士学士论 文].上海：上海海运学院，2003 [8] 刘雨. 船舶柴油发电机组的建模与智能控制的 仿真研究：[硕士学士论文].大连：大连海事大学，2002. 由图15可以看出，一、二号机组发出的功率 发生振荡。一号机组幅值变化范围为2.804×105～3.225×105 W，二号机组幅值变化范围为2.810×105～3.220×105 W。并且一台机组发出功率增加的同时，另一台机组发出功率在减小，负载上消耗的功率出现微小的振荡，幅值变化范围为5.956×105～6.053×105 W，并联机组发生功率交互振荡，相比于图14振荡程度严重。 4 结论 本文通过对柴油发电机组功率交互振荡现象的分析、仿真，我们得出干扰油量△gr是功率交互振荡产生的关键量，并联运行的机组相互间干扰油量差别越大，功率振荡的程度就越严重，因此干扰油量的幅值k、角频率ω、初相角α是影 ======================================= [海装动态] 俄称若印度同意加价12亿美元将尽快 交付航母 据美国《空中打击》网站2009年9月3日报道 俄罗斯国防技术公司(Rostekhnologii)总裁在莫斯科的一场记者招待会上向外界透露，俄罗斯和印度双方于2009年10月中旬签署一份关于“戈尔什科夫海军元帅”号航空母舰改装所需资金的新协议。 根据俄罗斯国防产品出口公司(Rosoboronexport)和印度海军于2004年签署的一份价值15亿美元的合同，包括米格－29K舰载机26 在内，“戈尔什科夫海军元帅”号航母的改装工作 应于2008年完成。但是俄罗斯后来表示，由于低估了该航母改装的规模和所需费用因而要求再加价12亿美元，印度方面则认为俄方此要求难以接受。经过一再延期和长时间谈判，印度曾于2008年2月表示可以接受6亿美元的加价。但俄罗斯方面对此价格不满意致使该问题一直拖延至今，目前相关谈判仍在进行。 俄罗斯方面已经保证，如果自己提出的加价12亿美元的要求得到落实，将会尽快完成“戈尔什科夫海军元帅”号航母的改装并于2012年交付给印度。 因篇幅问题不能全部显示，请点此查看更多更全内容