电力系统低频减载的单调控制特性_刘印.pdf
《电力系统低频减载的单调控制特性_刘印.pdf》由会员分享,可在线阅读,更多相关《电力系统低频减载的单调控制特性_刘印.pdf(8页珍藏版)》请在咨信网上搜索。
1、第 43 卷 第 7 期2023 年 7 月电 力 自 动 化 设 备Electric Power Automation EquipmentVol.43 No.7Jul.2023电力系统低频减载的单调控制特性刘印,陈民权,李京,张谦,甘德强(浙江大学 电气工程学院,浙江 杭州 310027)摘要:新能源占比的不断提高使得频率稳定问题日渐突出,因此迫切需要研究系统的频率稳定量化评估分析方法,也需对频率安全控制的最后一道防线低频减载进行更深入的分析。将单调控制系统理论运用到低频减载分析过程中,根据单调控制系统的输入-输出判定条件,推导单机和多机系统低频减载时频率分别需要满足的保序特性量化关系。此外
2、,分析时采用全状态模型,在保证电网规律性和解析性的同时获得更高的可扩展性,综合频率与功角、电压以及相关参数间的耦合关系。并结合灵敏度分析了电网内部变量之间的相关关系,以指导电网故障时的参数设置。最后,分析目前光伏参与调频的3种常见情况,分别说明单调控制理论的适用性,并通过算例进行验证。关键词:电力系统;低频减载;单调控制;保序;全状态模型;灵敏度中图分类号:TM761+.2 文献标志码:ADOI:10.16081/j.epae.2022100060 引言截至2021年底,我国风力发电和太阳能发电装机容量已分别达330 GW和310 GW,发电量同比分别增长16.6%和20.9%1。随着新能源占
3、比的不断提高,电力系统电力电子化特征愈发凸显,严重影响系统的频率稳定性。电力系统频率响应能力主要由系统的转动惯量与调频能力决定,由于新能源一般通过电力电子换流器与电网接口,在常规控制下其输送功率与电网频率解耦,不具备传统电机按照惯性时间常数分配扰动功率的能力,不能主动响应电网频率变化,从而不具备惯量及一次调频能力,使得系统频率调节能力变差23。并且新能源出力受电压影响显著,在出现扰动时进一步恶化系统频率特性,导致系统频率快速跌落或飙升,可能触发低频减载,严重时甚至会造成系统解列。为提升高占比新能源电网频率的认识水平,迫切需要深入研究系统的频率稳定分析方法。目前频率稳定分析方法主要包括时域仿真法
4、、模型分析法、机器学习法。时域仿真法是将电力系统等效为1组微分代数方程,用仿真软件进行时域仿真,得到扰动后系统的频率变化曲线,判断频率变化曲线最低点、频率变化率等关键指标是否越限4。模型分析法通过对时域仿真的详细模型进行合理的假设和简化,进而建立表征频率动态特征量的解析化表达式,求解频率特征。机器学习法指通过对数值仿真获得的样本数据进行特征提取和离线学习,将预想故障及运行数据输入已训练好的机器学习模型,以实现频率响应快速计算5。时域仿真法可利用详细的系统频率模型,计算结果精度相对较高,具有较好的可扩展性,但存在计算量大、耗时长的特点,并且由于新能源的加入,电网不确定因素增加,复杂度更高,系统更
5、难以满足频率安全在线评估需求6。机器学习法的样本数据必须由数值仿真获得,其精度本质上仍依赖于时域仿真计算。此外,上述2种方法对电网的可解释性和泛化能力较弱。系统频率动态及相关参数的解析关系以及电网规律性和解析性结论的获取主要依靠模型分析法7。为提升分析效率,模型分析法中通常采用平均系统频率(average system frequency,ASF)模型8和系统频率响应(system frequency response,SFR)模型9。这2种模型都具有降低系统规模、计算量小以及直观快速地求解频率特征的特点,但这2种模型都忽略了电压动态变化、网络拓扑和机组之间的影响,并忽略了频率与功角、电压之间
6、的耦合关系。例如当机组间频率振荡或系统电压发生较大偏移时,计算结果与真实值之间往往存在较大偏差。只有当系统采用全状态模型10,即综合考虑机组、负荷、网络以及控制器的动态特性时,才能在保证电网规律性和解析性的同时获得更高的精确度。此外,模型分析法常结合灵敏度分析探究物理量间的微分关系,即利用灵敏度刻画因变量对自变量的敏感程度。比如文献 11 采用轨迹灵敏度分析再热时间常数、惯性常数等不同参数对系统频率及其变化率的影响。文献 12 利用特征根灵敏度分析调速器参数对频率偏移峰值的影响。文献 13 根据阻尼灵敏度大小选择最佳光伏并网节点。而本文利用灵敏度分析在低频减载故障下相关变量间的关系,表征当系统
7、收稿日期:20220518;修回日期:20220818在线出版日期:20221012基金项目:国家自然科学基金资助项目(U2166601)Project supported by the National Natural Science Foundation of China(U2166601)182第 7 期刘印,等:电力系统低频减载的单调控制特性处于某一初始条件或参数发生微小变化时,利用系统中物理量之间的微分或代数关系,得到系统某一因变量对另一因变量的变化趋势及变化大小,以探究系统各变量之间的相关性。本文将单调控制理论运用到低频减载问题的分析中,论证频率响应曲线对应不同切载比率的保序关系。
8、采用全状态模型,综合考虑系统频率与功角、电压间的相互作用,尽可能真实地分析扰动作用下系统的动态行为。并通过灵敏度分析方法,分析电网内部变量之间的相关关系,以指导电网故障时的控制参数设置。1 单调控制系统理论单调动态系统理论已有悠久的历史,最早的经典结果是Muller-Kamke定理,其发现常微分方程解关于初值存在保序性质,即单调性,并将满足此性质的系统称为单调系统。但此时单调性仅指系统的解与状态量的对应关系。经过不断发展,Angeli 和Sontag将其推广到包含输入控制的系统中,即单调控制系统,并得到其输入量和输出量之间满足单调保序特性的条件和特征14。考虑输入-输出的控制系统如式(1)所示
9、。x?=f(x,u)y=h(x,u)(1)式中:状态向量x定义在X Rn上,n为状态向量的维数;输入向量u定义在U Rp上,p为控制向量的维数;输出向量y定义在Y Rm上;映射 f:X UIRn(IRn指有序欧几里得空间)在(x,u)上连续,在x上为局部李普希茨连续;映射h:XY 连续。若式(1)所示系统满足式(2)所示条件,则称其为单调控制系统。x1,x2X,u1,u2U,x1x2,u1u2 h(t,x1,u1)h(t,x2,u2),t0(2)式中:x1、x2和u1、u2分别为x、u中的元素。对于系统单调控制系统保序特性的确定可参考文献 15 中的结论:在任意时间断面上,输出向量y均可由输入
10、向量 u 和状态向量 x 经过映射 h 得到。假设在 t时刻,系统状态向量中包含 x1、x2这 2个元素,对应的输入量为u1、u2,可得单调控制系统的保序特性判定条件,如式(3)所示。yi(t,x2,u2)-yi(t,x1,u1)=01j=1nhixj(x1+r(x2-x1),u1+r(u2-u1)(xj2-xj1)dr+01j=1phiuj(x1+r(x2-x1),u1+r(u2-u1)(uj2-uj1)dr0(3)式中:yi为第i个输出量;上标 j 表示各变量的第j 个分量;r为积分变量。2 电力系统频率全状态模型2.1单机系统电力系统频率全状态模型是一个闭环控制系统,通过考虑全网机组、网
11、络结构的动态特性,从而尽可能地反映系统在扰动作用下的动态行为。由于单机系统是分析频率稳定的最简系统,并能有效反映系统有功频率控制的动态过程,以附录A图A1所示单机系统为例开展频率全状态模型研究。利用网络间代数关系消去代数方程部分,从而得到1组表征电网动态的纯微分方程。原动机调速器的动态模型为16:TSVdPSVdt=-PSV+PC-1RD(-1)(4)TCHdPMdt=-PM+PSV(5)式中:PSV、PC、PM分别为汽室输入功率、功率改变值、机械功率;TSV为调速器的时间常数,其取值通常小于 100 ms;TCH为汽室的时间常数,其取值范围为 0.2,2 s;为同步发电机转子转速;1/RD为
12、调速器增益,RD为静调压系数,其典型的取值范围为0.03,0.06。功率改变值PC作为1个控制输入项,可以是固定常数或是自动发电控制的输出。当PC为固定常数时,考虑稳态平衡点可得:PC=PMref(6)式中:PMref为参考机械功率。同步发电机采用 3 阶模型,并忽略凸极效应。当需要考虑暂态过程中转速变化对转子运动方程的影响时,同步发电机的转子运动方程为:Mddt=PM-Pe-D()-ref(7)式中:M为惯性时间常数;Pe为电磁功率;D为负荷-阻尼常数;ref为同步发电机转子额定转速。文献17 分析了同步发电机建模中转子运动方程考虑转速变化与否对仿真计算结果的影响:绝大多数情况下忽略转速变化
13、会使稳定性评估偏保守;但对于同步发电机3阶模型,其是否考虑转速变化对结果影响很小。为简化分析,本文忽略转速变化因素,当同步发电机配置比例型励磁电压调节器时,其动态模型为:|?=s()-refM?=PM-Pe-D()-refTd0E?q=Efd-Eq-()Xd-XdIdTAE?fd=KAUref-Efd-KAUG(8)式中:s为同步发电机转子标称转速;为同步发电机功角;Xd、Xd分别为同步发电机d轴电抗和d轴暂态电抗;Eq、Efd分别为同步发电机q轴暂态电动势相183电 力 自 动 化 设 备第 43 卷量Eq、励磁电压相量Efd的幅值;UG为同步发电机机端电压幅值;Id为d轴定子电流幅值;Td
14、0、TA分别为同步发电机d轴暂态时间常数和励磁调节器时间常数;KA为调节器放大系数;Uref为内部电压调节参考值,由潮流数据确定,如式(9)所示。Uref=UG0+Efd0/KA(9)式中:UG0为稳态机端电压幅值;Efd0为稳态励磁电压幅值。单机系统中同步发电机的电磁功率Pe为:Pe=EqULXsin(10)式中:UL为负荷节点电压相量UL的幅值;为Eq和UL的相角差;X为系统电抗,X=Xd+Xe,Xe为系统外接电抗。将式(10)代入式(8),消去代数方程,从而得到式(4)、(5)、(8)所示微分方程组以描述单机系统频率的全状态模型。2.2多机系统多机系统参数设置与单机系统一致。只需将式(4
15、)、(5)、(8)所示单机系统模型的微分方程扩展成矩阵形式,详细推导过程见附录A式(A1)(A7),最终可得多机系统频率的全状态模型如附录 A 式(A8)所示。3 单机系统频率单调验证3.1单机系统低频减载中的单调控制现象基于单机系统,设置同步发电机参数如附录A表A1所示,并取外接电抗(包括变压器和输电线路的电抗)Xe=0.1 p.u.,负荷为恒阻抗类型,其有功功率为80 MW,无功功率为60 Mvar;同步发电机容量等于系统基准容量SB,取SB=100 MV A。由单机系统全状态模型可知:状态向量x包含汽室输入功率、机械功率、同步发电机功角、转子转速、q轴暂态电动势Eq、励磁电压Efd;控制
16、量(即输入量)为负荷节点切载比率u;输出量y为系统频率。故单机系统数学模型可写成式(1)所示的输入-输出系统形式。首先分析代数方程,由电路方程可得考虑切载比率的负荷电压相量UL的解析表达式为:UL=RL+jXLj()Xd+Xe()1-u+()RL+jXLEqej(11)式中:RL+jXL为等值负荷阻抗。令Xu=()Xd+Xe()1-u+XL,由式(11)可得负荷电压幅值UL的简化表达式为:UL=()R2L+X2L()R2L+X2uR2L+X2uEq=K1Eq(12)将式(12)代入式(10),以q轴暂态电动势Eq表示电磁功率。对于稳态平衡点,将式(4)、(5)、(8)所示微分方程的等号左端项置
17、0,可得输出量表达式为:y=RD|PMref-K1(Eq)2Xsin +ref(13)由式(13)可知:稳态频率的变化量由系数K1、Eq和sin 决 定。不 同 切 载 比 率 下K1、Eq、sin 和K1(Eq)2sin 的变化情况见表 1,表中Eq为标幺值。由表可知:系数K1随切载比率的增加而上升,逐渐趋近于 1;Eq与sin 均随切载比率的增加而下降。结合式(13)和表1可知:切载后频率的稳态平衡点随切载比率的上升而增加,与附录A图A2所示不同切载比率频率变化特征相互印证。并由图A2可知:系统除稳态过程满足保序特性外,在过渡过程中同样严格满足保序特性。3.2基于单调控制理论的动态过程分析
18、下面基于单调控制理论,对上述现象展开模型分析。对于不同输入量对频率输出量的影响,其输出方程由转子运动方程决定,如式(14)所示。h=1M0 t|PM()-K1(Eq()2Xsin d(14)由式(14)可知:过渡过程输出量与状态量(PM、Eq)、控制量(切载比率u)相关。考虑具体的两节点系统减载场景,构造式(15)所示输入-输出系统。|TCHP?M=-PM+PSVTSVP?SV=-PSV+PC-()-ref/RDTd0E?q=Efd-Eq-()Xd-XdIdh=1M0 t|PM()-K1(Eq()2Xsin d(15)结合式(3)、(15),可得过渡过程中任意t时刻对应不同输入量u1、u2(u
19、2u1)时的判定条件表达式为:h()t,x2,u2-h()t,x1,u1=(PM2-PM1)01hPM(PM1+r(PM2-PM1)dr+(Eq2-Eq1)01hEq(Eq1+r(Eq2-Eq1)dr+(u2-u1)01hu(u1+r(u2-u1)dr(16)式中:变量下标1、2分别表示不同输入量的对应值。由于判断式(16)的正负符号特征涉及积分后求表1不同切载比率下K1、Eq、sin和K1(Eq)2sin的变化Table 1Variations of K1,Eq,sin and K1(Eq)2sinunder different load-shedding ratiosu%020406080
20、100K10.805 20.840 80.878 40.917 60.958 31.000 0Eq1.242 01.232 11.221 91.211 21.200 21.188 6sin 0.225 50.188 40.147 60.102 80.053 70K1(Eq)2sin0.280 00.240 40.193 50.138 40.074 10184第 7 期刘印,等:电力系统低频减载的单调控制特性导,且没有关于积分项的显式表达式。为简化分析,可先判断积分号内表达式在任意时刻t、不同切载比率下的输出量变化量h的正负符号特征,如式(17)所示。通过判断其大小,再从初始时刻到t时刻积分,得
21、到式(16)的符号特征,从而判断系统是否满足单调保序特性。h()t,x2,u2-h()t,x1,u1=1M(PM2()t-PM1()t+Pe1()t-Pe2()t)=PM2()tM-PM1()tM+K11(Eq1()t)2MXsin1-K21(Eq2()t)2MXsin2(17)对于式(17),t时刻输出量变化量h的单调保序性判定条件为:h()t,x2,u2-h()t,x1,u1=()PM2-PM101hPMdr+()Eq2-Eq101hEqdr+()u2-u101hudr(18)其中t时刻输出量变化量h对应各状态量的灵敏度表达式分别为:hPM=1M(19)hEq=-2K1EqMXsin(20
22、)设置在0.5 s时切除不同比例有功、无功负荷。对于h/PM,由式(19)可知其值为常数,且恒为正。对于h/Eq,在不同切载比率其仿真波形如附录A图A3所示,由图可知其值恒小于0。h/Eq与系数K1、sin 及当前时刻暂态电动势幅值相关,而系数K1和sin 在切载时发生突变,相应地h/Eq也在切载后发生突变。此外,h/PM、h/Eq的正负符号固定,说明输出量变化量h与各状态量之间相互作用明确。以h/PM为例,机械功率的增加将导致系统频率的提升,这符合传统同步发电机的运行特性。对于灵敏度h/u,因切载比率和所有状态量均相关,其关于频率的显式表达式求解难度大。下面通过差分法分别计算输出量和输出量变
23、化量对输入量的偏导,结果如附录A表A2所示。根据微积分基本公式,式(18)中机械功率时域响应差值PM2()t,u2-PM1()t,u1关于输入量的初值差存在式(21)所示关系。PM2()t,u2-PM1()t,u1=01PMu(0)()u2-u1dr(21)式中:u(0)为输入量初值。由于TCH=0.2 s,TSV=0.02 s,可近似认为TCHTSV。如果考虑再热器系统,一般再热器时间常数的取值范围为 6,12 s,此时TCH将更大。因此,在求解机械功率变化的主要特征时,可近似忽略调速器调节过程,从而获得式(22)所示的降阶模型。|TCHP?M=-PM+PSV0=-PSV+PMref-()-
24、ref/RD(22)从而可得:TCHP?M=-PM+PMref-()-ref/RD(23)参考文献 15 的推导,可得:minPMu(0)=-1RDhu(0)(24)借助仿真数据可得PM/u(0)的时域变化结果如附录 A 图 A4 所示。将式(24)所得的最终稳态值h/u与表A2数据进行比较,可验证式(24)的有效性。式(18)中暂态电动势时域响应差值可直接通过仿真数据得到。通过以上分析可将式(18)化简,由于仅取值点(取值点为积分项中偏导函数的括号内信息)为积分变量的函数,在输入量相差不大的情况下,可近似取被积函数变量中值,从而将偏导结果视为恒定进行积分估算,如式(25)所示。h()t,x2
25、,u2-h()t,x1,u1=-()u2-u1RDMhu(0)+()Eq2-Eq1(-2K1EqMXsin)+()u2-u1hu(25)设式(25)等号右侧依次为机械功率项、暂态电动势项和控制量项,以上变量在不同切载比率区间下的时域变化情况如图1所示。由图可知:机械功率项恒为负值,这说明机械功率具有维持系统频率不变的作用,可以提升频率稳定性;由暂态电动势项波形可看出电力系统频率与电压存在动态耦合关系;控制量项对频率变化量的影响恒为正,这说明切载比率和频率变化量呈正相关性,且切载比率越高,频率上升速率越快;通过纵向对比可知输出频率变化量的正负符号主要由机械功率和控制量项决定。图1机械功率项、暂态
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 电力系统 低频 单调 控制 特性 刘印
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。