CN108767898A - 一种考虑风电与直流综合作用的电力系统暂态稳定分析方法 - Google Patents
一种考虑风电与直流综合作用的电力系统暂态稳定分析方法 Download PDFInfo
- Publication number
- CN108767898A CN108767898A CN201810609425.0A CN201810609425A CN108767898A CN 108767898 A CN108767898 A CN 108767898A CN 201810609425 A CN201810609425 A CN 201810609425A CN 108767898 A CN108767898 A CN 108767898A
- Authority
- CN
- China
- Prior art keywords
- wind
- power
- planes
- mechanical output
- equivalent
- 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
- 230000005611 electricity Effects 0.000 title claims abstract description 73
- 230000001052 transient effect Effects 0.000 title claims abstract description 53
- 238000004458 analytical method Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000010354 integration Effects 0.000 claims abstract description 25
- 230000000694 effects Effects 0.000 claims abstract description 23
- 230000001133 acceleration Effects 0.000 claims abstract description 15
- 230000008859 change Effects 0.000 claims abstract description 9
- 238000004445 quantitative analysis Methods 0.000 claims abstract description 8
- 238000004451 qualitative analysis Methods 0.000 claims abstract description 6
- 230000001360 synchronised effect Effects 0.000 claims description 23
- 238000002347 injection Methods 0.000 claims description 19
- 239000007924 injection Substances 0.000 claims description 19
- 230000009471 action Effects 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 13
- 239000013598 vector Substances 0.000 claims description 9
- 230000007423 decrease Effects 0.000 claims description 8
- 230000007246 mechanism Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 238000011084 recovery Methods 0.000 claims description 5
- 230000005684 electric field Effects 0.000 claims description 4
- 238000013016 damping Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- XOOUIPVCVHRTMJ-UHFFFAOYSA-L zinc stearate Chemical compound [Zn+2].CCCCCCCCCCCCCCCCCC([O-])=O.CCCCCCCCCCCCCCCCCC([O-])=O XOOUIPVCVHRTMJ-UHFFFAOYSA-L 0.000 claims description 3
- 230000008569 process Effects 0.000 description 4
- 230000004888 barrier function Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000013256 coordination polymer Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Classifications
-
- H02J3/386—
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种考虑风电与直流综合作用的电力系统暂态稳定分析方法,通过建立风电接入交直流非并联送出的系统暂态稳分析模型,提出大规模风电接入对系统暂态稳定影响的量化指标,提出将风电与直流综合作用等效为系统等值机械功率增量的定量分析方法,基于等面积定则提出了考虑风电与直流综合作用对加速面积与减速面积的改变的定性分析方法。本发明实现方便,分析方法更明确有效;适用性高,分析简便;适用范围广,分析效果好;应用场景广,利于解决实际工程问题。
Description
技术领域
本发明涉及电力系统稳定性分析技术领域,具体为一种考虑风电与直流综合作用的电力系统暂态稳定分析方法。
背景技术
双馈式感应发电机(doubly-fed induction generator,DFIG)具有机电解耦等特性在风力发电中得到了广泛应用,风火打捆接入送端可抑制有功大幅度波动,同时大容量、远距离的高压直流输电技术可满足大规模新能源送出的要求,但能源多样性与复杂的混联输电结构也给电网的稳定问题带来了新的挑战。风电场本身没有功角稳定问题,但不同规模和渗透率的风电会通过影响邻近火电机组从而影响系统的稳定运行。
目前综合考虑风电与直流输电对系统稳定性影响的研究较少,且大多是通过时域仿真得到,而当直流系统异步运行时,仅考虑风电的作用是不够的,直流与风电综合作用下电网的安全稳定问题亟待进一步研究。
如分析方法1:提出在故障期间将风电作用等效为阻抗的暂态稳定分析方法;可参考文献《双馈风电机组接入对电力系统功角稳定性的影响》;分析方法2:提出将双馈风机采用直流潮流法判断风电接入后对系统暂态稳定影响的判据;可参考文献《大规模风电集中接入对电力系统暂态功角稳定性的影响(一):理论基础》;分析方法3:提出在传输功率不同时将直流系统等效为恒定负荷的增量的分析方法,可参考文献《大规模直流异步互联系统受端故障引发送端稳定破坏的机理分析》;分析方法4:提出一种基于扩展等面积定则的含风电场的电力系统暂态稳定性分析方法;可参考文献《基于扩展等面积定则的含大规模风电场电力系统暂态稳定性分析》。
但现有分析方法均单独从定性或者定量分析风电接入对系统暂态稳定性的影响,分析层面单一;虽然考虑的影响因素较为全面,但实际在在某些情况下提出的判据仍然不能确定与风电接入前相比系统暂态稳定水平提高或降低,适用性较低;对交直流并联运行的情况,在分析风电的作用时忽略了直流系统的影响,目前暂无用于风电与直流综合作用对电网暂态稳定影响机理的分析方法。
发明内容
针对上述问题,本发明的目的在于提供一种既定量分析了风电与直流两者综合作用对系统暂态稳定的影响,同时又能基于等面积定则对两者综合作用对系统暂态稳定影响机理进行定性分析,分析过程明确且适用性高,更便于实际工程利用的考虑风电与直流综合作用的电力系统暂态稳定分析方法。技术方案如下:
一种考虑风电与直流综合作用的电力系统暂态稳定分析方法,包括以下步骤:
步骤1:根据扩展等面积法,得到单机无穷大系统的运动方程:
将受扰严重的同步机组定义为临界S机群,剩下的同步机组称为余下A机群,对全系统进行双机等值得到转子运动方程:
式中:Ms和Ma分别为S机群和A机群的惯性时间常数;δs和δa分别为S机群和A机群的转子角;Pms和Pma分别为S机群和A机群的等值机械功率;Pes和Pea分别为S机群和A机群的等值电磁功率,其表达式为:
其中,Mi、Pmi、Pei分别为S机群中第i个机组的惯性时间常数、等值机械功率和等值电磁功率;Mj、Pmj和Pej分别为A机群中第j个机组的惯性时间常数、等值机械功率和等值电磁功率;
进一步作单机无穷大等值并定义转子角为:
δ=δs-δa (6)
从而单机无穷大系统的运动方程为:
其中,
式中:M、Pm、Pe分别为单机无穷大系统的等值惯量、等值机械功率和等值电磁功率;G为网络节点导纳矩阵中的元素;γ为熄弧角;Ei、Ek、Ej、El分别为S机群与A机群中各台发电机的等值电势;Gij为发电机i与发电机j的互导纳;MT为整个系统的惯量,PC、C、D分别为等值电磁功率表达式中的各分项;
步骤2:建立含大规模风电的交直流混联三区域系统模型:
将大规模双馈风力机组成的风电场等效为一台双馈风机,送端系统为区域1,由双馈风机与常规火电厂的同步机按比例组合构成;区域1与区域2同步弱互联,区域2为交流受端区域,区域3为直流受端区域,其中区域2与区域3非并联运行;
步骤3:确定风电接入对系统暂态稳定影响分析指标:
假定系统失稳模式已知,将区域1内的同步机组称为受扰严重的临界S机群,区域2内的同步机组称为余下A机群;
当风电未接入时,系统的网络方程为
式中:分别为S机群和A机群节点注入电流列向量;分别为S机群和A机群节点注入电压列向量;Yss、Ysa、Yas、Yaa是已包含负荷阻抗的导纳矩阵的元素;
DFIG接入后,即等效电流源从风电并网点注入,则系统的网络方程变为:
式中:分别为风机注入电流和节点电压列向量;Ysw、Yww、Yws、Ywa、Yaw为风机接入后节点导纳矩阵中新增的元素,且Ywa、Yaw均为零;
为分析风电场接入对原系统同步机的影响,将DFIG等值纳入S机群,则对以上网络方程进行变换,得到:
式中:Y′ss为风电接入后区域1的自阻抗矩阵;△Iw为注入电流;
设风电机组接入比例为k,风电输出有功功率为Pw=kPs0,风电并网母线电压为Uw,则注入电流为:
△Iw=kPs0/Uw (11)
对公式(7)进行修正,其中与S机群自导纳相关的量只有PC,而C、D、Pem、γ均不变,则
式中,P′C为对式(7)中PC进行修正之后的值;△PC为修正量;
此外,风电接入后S机群的机械功率变为:
式中:Mw、Pmw分别为风电场等值惯量与机械功率;△Pm为风电作用下对机械功率的修正量;令则风机接入后的系统运动方程为:
对比式(7)与式(14),风电接入后对系统稳定性的影响等效表现为对S机群机械功率或者电磁功率的增量△Pmw-△Pew,变换为标准形式:
式中:P″m=P′m-△PC=Pm+△Pmw-△Pew;
步骤4:记及风电与直流综合作用的三区域系统暂态稳定的定量分析:
进一步考虑区域2直流系统的影响,得到三区域系统转子运动方程:
式中,PD13为直流输送功率;
将直流系统的影响等效为对S机群机械功率或者电磁功率的增量PD13,忽略发电机调速器动作以及负荷功率变化,式(16)简化为
其中,
式中:P″′m为考虑风电与直流综合作用后的三区域系统等值机械功率;P′em为考虑风电与直流综合作用后的三区域系统等值电磁功率最大值;
将风电与直流的综合作用等效为对系统等值机械功率的增量△Pmw-△Pew-PD13,基于该模型定量分析风电与直流对系统暂态功角稳定的影响。
进一步的,还包括基于扩展等面积定则的暂态稳定机理定性分析:
根据扩展等面积定则,基于所述式(17)与(18)绘制风电与直流综合作用下系统功角特性曲线:故障发生时刻为δ0;故障切除时刻为δc;系统功角最大时刻为δm;将三区域系统正常运行时的电磁功率曲线记为PI,系统等值机械功率为P′m,稳态运行点为a点;故障发生时的系统功角曲线记为PII,稳态运行为b点;故障切除时的时的系统功角曲线记为记为PIII;c点为故障切除瞬间系统运行点;d点为故障切除后电磁功率恢复至PII时系统运行点;g、f点为故障后未考虑风电与直流的制动作用时系统在δ0与δc时刻的等值机械功率;e点为故障后考虑风电与直流制动作用后系统在δc时刻的等值机械功率;h点为系统在δm时刻的运行点;i、k为故障切除后未考虑风电与直流的制动作用时系统在δm与δc时刻的等值机械功率;j为故障切除后考虑风电与直流制动作用后系统在δc时刻的等值机械功率;
故障发生时,同步发电机的机械功率不变,系统等值机械功率的变化分为两部分:由于区域1母线电压的下降,风电场的电磁功率下降△Pew,风机由于转速控制与变桨距控制的作用,风机机械功率在故障期间近似线性变化,因此△Pmw-△Pew使得机械功率不再是平行的直线,而是按一定斜率下降,该斜率由DFIG的桨距角控制特性以及转速-转矩特性等共同决定;其次,考虑严重故障下直流功率瞬降为零,机械功率增大PD13,总的等值机械功率增大△Pmw-△Pew-PD13变为P″′m;故风电场与直流系统综合作用下的系统等值机械功率P″′m在δ0~δc故障期间的沿着直线ge变化,且大于电磁功率,转子沿着PII加速至c点;在δc时刻,故障切除,系统电磁功率曲线恢复为PIII,同步电机运行于d点,加速面积为Sabcdeg;可知风电与直流系统的作用增大了加速面积,而风电场的制动作用减小了加速面积,减小量为Sefg,但两者在故障期间的综合作用增大了总加速面积,增加量为Sakdeg,不利于系统的暂态稳定;
故障切除后,等值机械功率的恢复也分为两部分:直流功率能够快速恢复到稳态值,机械功率P″′m减小PD13;风电场的无功控制能力使得并网点母线电压恢复,则△Pew逐渐恢复,并且受风电场控制特性的影响,△Pmw-△Pew作用下P″′m以一定斜率线性上升,故风电场与直流系统综合作用下的系统等值机械功率P″′m在δc~δm的故障清除期间的沿着直线ji变化,且系统等值机械功率小于电磁功率;转子减速功角δ不断增大,运行点由d点沿着曲线PIII运动到h点,减速面积为Sdhij;直流系统的阻尼作用和风电场的制动作用增大了减速面积,增加量为Sijk,有利于系统的暂态稳定。
本发明的有益效果是:本发明实现方便,分析方法更明确有效,方法1所采用的电阻在故障过程中不是单一变量,对动态分析带来困难,本发明采用的功率分析模型数据更加容易得到,得到的量化指标能定量分析风电接入对系统暂态稳定性影响;相比方法2,适用性高,分析简便,本方法所提量化指标适用性更高,且考虑了直流系统的作用;相比分析方法3:适用范围广,分析效果好。方法3仅考虑了直流系统的作用,本发明综合考虑了风电机组与直流系统的作用,适用的分析应用场景范围更宽广;相比分析方法4:适用性高,应用场景广。本发明还考虑了直流系统的作用,针对直流系统的换相失败预直流闭锁故障均有效,因此相比于仅考虑风电作用本分析方法更利于解决实际工程问题。
附图说明
图1为双馈异步风力发电机结构。
图2为含大规模风电的交直流混联三区域系统模型。
图3为风电与直流综合作用下系统功角特性曲线图。
图4为送端故障时系统的暂态响应曲线图;(a)发电机输出功率;(b)直流输送功率。
图5为三区域系统实际功角曲线。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细说明。如图1所示,由双馈异步风力发电机结构可知风力机捕获风能并将其转化为机械功率Pmw:
式中:ρ为空气密度;A为叶轮扫过面积;vw为风速;CP为功率转换系数,由叶尖速比λ和桨距角β共同确定;PB为功率基准值。当风速确定时,机械功率由CP决定。
以双馈异步风力发电机的定子电压矢量为q轴,在同步旋转dq0坐标系下建立的双馈异步电机模型为:
式中:uds、uqs、ψds、ψqs分别为定子电压与磁链在d轴和q轴的分量;udr、uqr、ψdr、ψqr为转子电压与磁链在d轴和q轴的分量;Rs、Rr分别为定、转子电阻,ids、iqs、idr、iqr分别为定、转子电流在d轴和q轴的分量,ωs为角频率。
高压直流外送系统采用准稳态模型如下:
式中:Ud0为空载直流电压,EL为交流母线电压;
Ud0r、Ud0i分别为整流侧和逆变侧空载直流电压;Id为直流电流;Rcr、Rci分别代表整流侧和逆变侧换相电抗;Pac为直流功率;α为触发延迟角;β为触发超前角;γ为熄弧角;μ为换相角。
根据扩展等面积定则(EEAC),系统受到大扰动后的失稳以两机群功角相对摆开为主导。假定系统失稳模式已知,将受扰严重的同步机组定义为临界S机群,剩下的同步机组称为余下A机群,对全系统进行双机等值得到转子运动方程:
式中:Ms和Ma分别为S机群和A机群的惯性时间常数;δs和δa分别为S机群和A机群的转子角;Pms和Pma分别为S机群和A机群的等值机械功率;Pes和Pea分别为S机群和A机群的等值电磁功率,其表达式为:
其中,Mi、Pmi、Pei分别为S机群中第i个机组的惯性时间常数、等值机械功率和等值电磁功率;Mj、Pmj和Pej分别为A机群中第j个机组的惯性时间常数、等值机械功率和等值电磁功率;
进一步作单机无穷大等值并定义转子角为:
δ=δs-δa (6)
从而单机无穷大系统的运动方程为
其中,
式中:M、Pm、Pe分别为单机无穷大系统的等值惯量、等值机械功率和等值电磁功率;G为网络节点导纳矩阵中的元素;γ为熄弧角;Ei、Ek、Ej、El分别为S机群与A机群中各台发电机的等值电势;Gij为发电机i与发电机j的互导纳;MT为整个系统的惯量,PC、C、D分别为等值电磁功率表达式中的各分项,无实际的物理意义。
基于上述模型,本实施例提出了大规模风电接入对系统暂态稳定影响的量化指标,建立了风电接入交直流非并联送出的系统暂态稳分析模型,提出一种将风电与直流综合作用等效为系统等值机械功率增量的定量分析方法,基于等面积定则提出了考虑风电与直流综合作用对加速面积与减速面积的改变的定性分析方法。具体步骤如下:
1)建立含大规模风电的交直流混联三区域系统模型:
大规模风电集中接入经交直流外送的系统一般有相同受端或不同受端,本发明主要研究交直流系统非并联运行不同受端的情况,建立含大规模风电接入的交直流混联三区域系统的模型。
如图2所示,将大规模双馈风力机组成的风电场等效为一台双馈风机,送端系统由双馈风机与常规火电厂的同步机按一定比例组合构成。区域1与区域2同步弱互联,其中E1∠δ1、E2∠δ2、E3∠δ3分别为等值发电机G1、G2、G3的电压及转子角,Uw为风机汇集点的电压,x12为区域1与区域2的互阻抗,x1、x2、x3、xs分别为G1、G2、G3和风机端节点到各母线的等值电抗,PL1、PL2、PL3为负荷,PA12、PD13分别为交流与直流输电系统输送的功率。
2)风电接入对系统暂态稳定影响分析指标
鉴于DFIG的快速恢复特性,相比机电暂态过程可忽略不计,其故障前和故障清除后的功率基本恒定,且假定DFIG满足:(1)风电场及并网点配置有无功补偿装置,能维持故障时电压跌落不低于0.2p.u;(2)风机具备良好的低电压穿越能力,并配置Crowbar保护装置;(3)DFIG的转速能在桨距角控制和转速-转矩特性的作用下维持在允许转速范围内,可将风机看成注入区域1的电流源,基于式(7),深入研究DFIG与直流对系统转子运动方程中具体参数与暂态稳定的影响。
当风电未接入时,系统的网络方程为:
式中:分别为S机群和A机群节点注入电流列向量;分别为S机群和A机群节点注入电压列向量;Yss、Ysa、Yas、Yaa是已包含负荷阻抗的导纳矩阵的元素。
DFIG接入后,即等效电流源从风电并网点注入,则系统的网络方程变为:
式中:分别为风机注入电流和节点电压列向;Ysw、Yww、Yws、Ywa、Yaw为风机接入后节点导纳矩阵中新增的元素,且Ywa、Yaw均为零。
为分析风电场接入对原系统同步机的影响,将DFIG等值纳入S机群,则可以对以上网络方程进行变换,得到:
式中:Y′ss为风电接入后区域1的自阻抗矩阵,△Iw为注入电流。
比较式(8)与式(10),由于DFIG与A机群并未直接相连,接入后改变了系统网络方程中区域1的注入电流和自导纳矩阵Yss。考虑风电等容量接入的情况,即保持送端总的输出功率Ps0不变,火电机组出力减少,设风电机组接入比例为k,风电输出有功功率为Pw=kPs0,风电并网母线电压为Uw,则注入电流为:
△Iw=kPs0/Uw (11)
对公式(7)进行修正,其中与S机群自导纳相关的量只有PC,而C、D、Pem、γ均不变,则
式中,P′C为对式(7)中PC进行修正之后的值;△PC为修正量;
此外,风电接入后S机群的机械功率变为
式中:Mw、Pmw分别为风电场等值惯量与机械功率,△Pm为风电作用下对机械功率的修正量。
令则风机接入后的系统运动方程为
对比式(7)与式(14),风电接入后对系统稳定性的影响可等效表现为对S机群机械功率或者电磁功率的增量(△Pmw-△Pew),变换为标准形式:
式中:P″′m=P′m-△PC=Pm+△Pmw-△Pew。
3)记及风电与直流综合作用的三区域系统暂态稳定的定量分析
进一步考虑区域2直流系统的影响,得到三区域系统转子运动方程
式中,PD13为直流输送功率。
类比于风电对系统暂态稳定的影响,直流系统的影响也可等效为对对S机群机械功率或者电磁功率的增量PD13,为简化推导,忽略发电机调速器动作以及负荷功率变化,式(16)可以简化为:
其中,
式中:P″′m为考虑风电与直流综合作用后的三区域系统等值机械功率;P′em为考虑风电与直流综合作用后的三区域系统等值电磁功率最大值。
式(17)与式(18)构成了图2所示含大规模风电接入的三区域系统暂态稳定分析模型,将风电与直流的综合作用等效为对系统等值机械功率的增量(△Pmw-△Pew-PD13),基于该模型分析可定量分析风电与直流对系统暂态功角稳定的影响。
4)基于扩展等面积定则的暂态稳定机理定性分析
根据扩展等面积定则,基于式(17)与(18),对大规模风电与直流系统综合作用下的送端系统暂态功角稳定进行分析。
如图3所示,故障发生时刻为δ0;故障切除时刻为δc;系统功角最大时刻为δm;将三区域系统正常运行时的电磁功率曲线记为PI,系统等值机械功率为P′m,稳态运行点为a点;故障发生时的系统功角曲线记为PII,稳态运行为b点;故障切除时的时的系统功角曲线记为记为PIII;c点为故障切除瞬间系统运行点;d点为故障切除后电磁功率恢复至PII时系统运行点;g、f点为故障后未考虑风电与直流的制动作用时系统在δ0与δc时刻的等值机械功率;e点为故障后考虑风电与直流制动作用后系统在δc时刻的等值机械功率;h点为系统在δm时刻的运行点;i、k为故障切除后未考虑风电与直流的制动作用时系统在δm与δc时刻的等值机械功率;j为故障切除后考虑风电与直流制动作用后系统在δc时刻的等值机械功率。
故障发生时,系统功角曲线变为PII,运行于b点,此时,同步发电机的机械功率不变,系统等值机械功率的变化分为两部分:由于区域1母线电压的下降,风电场的电磁功率下降△Pew,风机由于转速控制与变桨距控制的作用,风机机械功率在故障期间近似线性变化,因此△Pmw-△Pew使得机械功率不再是平行的直线,而是按一定斜率下降,该斜率由DFIG的桨距角控制特性以及转速-转矩特性等共同决定;其次,考虑严重故障下直流功率瞬降为零,机械功率增大PD13,总的等值机械功率增大△Pmw-△Pew-PD13变为P″′m。故风电场与直流系统综合作用下的系统等值机械功率P″′m在δ0~δc故障期间的沿着直线ge变化,且大于电磁功率,转子沿着PII加速至c点。在δc时刻,故障切除,系统电磁功率曲线恢复为PIII,同步电机运行于d点,加速面积为Sabcdeg。可知风电与直流系统的作用增大了加速面积,而风电场的制动作用减小了加速面积,减小量为Sefg,但两者在故障期间的综合作用增大了总加速面积,增加量为Sakdeg,不利于系统的暂态稳定。
故障切除后,等值机械功率的恢复也分为两部分:直流功率能够快速恢复到稳态值,机械功率P″′m减小PD13;风电场的无功控制能力使得并网点母线电压恢复,则△Pew逐渐恢复,并且受风电场控制特性的影响,△Pmw-△Pew作用下P″′m一定斜率线性上升,故风电场与直流系统综合作用下的系统等值机械功率P″′m在δc~δm的故障清除期间的沿着直线ji变化,且系统等值机械功率小于电磁功率。转子减速但是功角δ不断增大,运行点由d点沿着曲线PIII运动到h点,减速面积为Sdhij。直流系统的阻尼作用和风电场的制动作用增大了减速面积,增加量为Sijk,有利于系统的暂态稳定。
综上,对于含大规模风电接入的交直流混联系统,风电与直流系统的综合作用使得故障期间的加速面积增大,不利于系统暂态稳定。故障切除后两者的快速恢复特性与风机的制动作用也增加了减速面积,有利于系统的暂态稳定,但加速面积增量大于减速面积增量,因此风电与异步互联直流综合作用下系统的暂态稳定性降低。
5)控制效果验证
故障期间风机与直流系统的输出功率如图4所示,(△Pmw-△Pew-PD13)作用下三区域系统暂态功角曲线如图5所示。区域1母线发生故障,风电未接入时,三区域系统初始功角为16°,首摆最大功角达到19.4°,非周期振荡约6s后恢复稳定。风电接入后,等值的三区域系统初始功角增大至19.2°,首摆最大功角增大为32.5°,5s后恢复至稳定状态。因此,根据以上仿真结果可知,对于交直流非并联运行的三区域系统,大规模风电接入后与直流的相互作用降低了系统的暂态功角稳定性,这也验证了机理分析的正确性。
6)分析方法分析过程
a.基于故障期间风机测量得到的实际电磁功率与机械功率,得到量化指标△Pmw-△Pew;
b.根据故障期间直流系统的直流功率瞬降计算△Pmw-△Pew-PD13;
c.根据△Pmw-△Pew-PD13,基于等面积定则,分析加速面积与减速面积,判断对系统暂态稳定性的影响程度。
Claims (2)
1.一种考虑风电与直流综合作用的电力系统暂态稳定分析方法,其特征在于,包括以下步骤:
步骤1:根据扩展等面积法,得到单机无穷大系统的运动方程:
将受扰严重的同步机组定义为临界S机群,剩下的同步机组称为余下A机群,对全系统进行双机等值得到转子运动方程:
式中:Ms和Ma分别为S机群和A机群的惯性时间常数;δs和δa分别为S机群和A机群的转子角;Pms和Pma分别为S机群和A机群的等值机械功率;Pes和Pea分别为S机群和A机群的等值电磁功率,其表达式为:
其中,Mi、Pmi、Pei分别为S机群中第i个机组的惯性时间常数、等值机械功率和等值电磁功率;Mj、Pmj和Pej分别为A机群中第j个机组的惯性时间常数、等值机械功率和等值电磁功率;
进一步作单机无穷大等值并定义转子角为:
δ=δs-δa (6)
从而单机无穷大系统的运动方程为:
其中,
式中:M、Pm、Pe分别为单机无穷大系统的等值惯量、等值机械功率和等值电磁功率;G为网络节点导纳矩阵中的元素;γ为熄弧角;Ei、Ek、Ej、El分别为S机群与A机群中各台发电机的等值电势;Gij为发电机i与发电机j的互导纳;MT为整个系统的惯量,PC、C、D分别为等值电磁功率表达式中的各分项;
步骤2:建立含大规模风电的交直流混联三区域系统模型:
将大规模双馈风力机组成的风电场等效为一台双馈风机,送端系统为区域1,由双馈风机与常规火电厂的同步机按比例组合构成;区域1与区域2同步弱互联,区域2为交流受端区域,区域3为直流受端区域,其中区域2与区域3非并联运行;
步骤3:确定风电接入对系统暂态稳定影响分析指标:
假定系统失稳模式已知,将区域1内的同步机组称为受扰严重的临界S机群,区域2内的同步机组称为余下A机群;
当风电未接入时,系统的网络方程为
式中:分别为S机群和A机群节点注入电流列向量;分别为S机群和A机群节点注入电压列向量;Yss、Ysa、Yas、Yaa是已包含负荷阻抗的导纳矩阵的元素;DFIG接入后,即等效电流源从风电并网点注入,则系统的网络方程变为:
式中:分别为风机注入电流和节点电压列向量;Ysw、Yww、Yws、Ywa、Yaw为风机接入后节点导纳矩阵中新增的元素,且Ywa、Yaw均为零;
为分析风电场接入对原系统同步机的影响,将DFIG等值纳入S机群,则对以上网络方程进行变换,得到:
式中:Ys's为风电接入后区域1的自阻抗矩阵;△Iw为注入电流;
设风电机组接入比例为k,风电输出有功功率为Pw=kPs0,风电并网母线电压为Uw,则注入电流为:
△Iw=kPs0/Uw (11)
对公式(7)进行修正,其中与S机群自导纳相关的量只有PC,而C、D、Pem、γ均不变,则
式中,P′C为对式(7)中PC进行修正之后的值;△PC为修正量;
此外,风电接入后S机群的机械功率变为:
式中:Mw、Pmw分别为风电场等值惯量与机械功率;△Pm为风电作用下对机械功率的修正量;
令则风机接入后的系统运动方程为:
对比式(7)与式(14),风电接入后对系统稳定性的影响等效表现为对S机群机械功率或者电磁功率的增量△Pmw-△Pew,变换为标准形式:
式中:P″m=P′m-△PC=Pm+△Pmw-△Pew;
步骤4:记及风电与直流综合作用的三区域系统暂态稳定的定量分析:
进一步考虑区域2直流系统的影响,得到三区域系统转子运动方程:
式中,PD13为直流输送功率;
将直流系统的影响等效为对S机群机械功率或者电磁功率的增量PD13,忽略发电机调速器动作以及负荷功率变化,式(16)简化为
其中,
式中:P″′m为考虑风电与直流综合作用后的三区域系统等值机械功率;P′em为考虑风电与直流综合作用后的三区域系统等值电磁功率最大值;
将风电与直流的综合作用等效为对系统等值机械功率的增量△Pmw-ΔPew-PD13,基于该模型定量分析风电与直流对系统暂态功角稳定的影响。
2.根据权利要求1所述的考虑风电与直流综合作用的电力系统暂态稳定分析方法,其特征在于,还包括基于扩展等面积定则的暂态稳定机理定性分析:
根据扩展等面积定则,基于所述式(17)与(18)绘制风电与直流综合作用下系统功角特性曲线:故障发生时刻为δ0;故障切除时刻为δc;系统功角最大时刻为δm;将三区域系统正常运行时的电磁功率曲线记为PI,系统等值机械功率为P′m,稳态运行点为a点;故障发生时的系统功角曲线记为PII,稳态运行为b点;故障切除时的时的系统功角曲线记为记为PIII;c点为故障切除瞬间系统运行点;d点为故障切除后电磁功率恢复至PII时系统运行点;g、f点为故障后未考虑风电与直流的制动作用时系统在δ0与δc时刻的等值机械功率;e点为故障后考虑风电与直流制动作用后系统在δc时刻的等值机械功率;h点为系统在δm时刻的运行点;i、k为故障切除后未考虑风电与直流的制动作用时系统在δm与δc时刻的等值机械功率;j为故障切除后考虑风电与直流制动作用后系统在δc时刻的等值机械功率;
故障发生时,同步发电机的机械功率不变,系统等值机械功率的变化分为两部分:由于区域1母线电压的下降,风电场的电磁功率下降△Pew,风机由于转速控制与变桨距控制的作用,风机机械功率在故障期间近似线性变化,因此△Pmw-△Pew使得机械功率不再是平行的直线,而是按一定斜率下降,该斜率由DFIG的桨距角控制特性以及转速-转矩特性等共同决定;其次,考虑严重故障下直流功率瞬降为零,机械功率增大PD13,总的等值机械功率增大△Pmw-△Pew-PD13变为P″′m;故风电场与直流系统综合作用下的系统等值机械功率P″′m在δ0~δc故障期间的沿着直线ge变化,且大于电磁功率,转子沿着PII加速至c点;在δc时刻,故障切除,系统电磁功率曲线恢复为PIII,同步电机运行于d点,加速面积为Sabcdeg;可知风电与直流系统的作用增大了加速面积,而风电场的制动作用减小了加速面积,减小量为Sefg,但两者在故障期间的综合作用增大了总加速面积,增加量为Sakdeg,不利于系统的暂态稳定;
故障切除后,等值机械功率的恢复也分为两部分:直流功率能够快速恢复到稳态值,机械功率P″′m减小PD13;风电场的无功控制能力使得并网点母线电压恢复,则△Pew逐渐恢复,并且受风电场控制特性的影响,△Pmw-△Pew作用下P″′m以一定斜率线性上升,故风电场与直流系统综合作用下的系统等值机械功率P″′m在δc~δm的故障清除期间的沿着直线ji变化,且系统等值机械功率小于电磁功率;转子减速功角δ不断增大,运行点由d点沿着曲线PIII运动到h点,减速面积为Sdhij;直流系统的阻尼作用和风电场的制动作用增大了减速面积,增加量为Sijk,有利于系统的暂态稳定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810609425.0A CN108767898B (zh) | 2018-06-13 | 2018-06-13 | 一种考虑风电与直流综合作用的电力系统暂态稳定分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810609425.0A CN108767898B (zh) | 2018-06-13 | 2018-06-13 | 一种考虑风电与直流综合作用的电力系统暂态稳定分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108767898A true CN108767898A (zh) | 2018-11-06 |
CN108767898B CN108767898B (zh) | 2021-07-06 |
Family
ID=64021708
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810609425.0A Active CN108767898B (zh) | 2018-06-13 | 2018-06-13 | 一种考虑风电与直流综合作用的电力系统暂态稳定分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108767898B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110890751A (zh) * | 2019-10-18 | 2020-03-17 | 东南大学 | 一种适用含高比例风电电力系统的暂态稳定特征保全方法 |
CN112054550A (zh) * | 2020-09-11 | 2020-12-08 | 武汉大学 | 一种多回特高压直流馈入的交流系统的直流功率互济方法 |
CN114221339A (zh) * | 2021-12-21 | 2022-03-22 | 上海电力大学 | 一种高风电占比系统的暂态稳定性判断方法 |
CN114268100A (zh) * | 2022-03-01 | 2022-04-01 | 国网江西省电力有限公司电力科学研究院 | 一种交直流混联电力系统暂态稳定性能分析方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070114796A1 (en) * | 2005-11-18 | 2007-05-24 | Garces Luis J | System and method for integrating wind and hydroelectric generation and pumped hydro energy storage systems |
CN102403721A (zh) * | 2011-11-23 | 2012-04-04 | 昆明理工大学 | 一种基于暂态功角稳定裕度的超实时重合时序整定方法 |
CN103094921A (zh) * | 2013-01-06 | 2013-05-08 | 宁夏电力公司电力科学研究院 | 直驱式风力发电机组机电暂态建模方法及机电暂态模型 |
CN103746404A (zh) * | 2014-01-24 | 2014-04-23 | 国家电网公司 | 一种风光火打捆直流外送系统稳定性评价方法 |
CN103904678A (zh) * | 2014-03-31 | 2014-07-02 | 中国南方电网有限责任公司电网技术研究中心 | 高压直流输电分段变速率低压限流单元的控制方法 |
CN106356885A (zh) * | 2016-09-19 | 2017-01-25 | 国家电网公司 | 一种风电并网系统暂态功角前二摆稳定性判断方法 |
CN107317354A (zh) * | 2017-06-30 | 2017-11-03 | 天津大学 | 一种含风电场的多机系统暂态功角稳定性分析方法 |
-
2018
- 2018-06-13 CN CN201810609425.0A patent/CN108767898B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070114796A1 (en) * | 2005-11-18 | 2007-05-24 | Garces Luis J | System and method for integrating wind and hydroelectric generation and pumped hydro energy storage systems |
CN102403721A (zh) * | 2011-11-23 | 2012-04-04 | 昆明理工大学 | 一种基于暂态功角稳定裕度的超实时重合时序整定方法 |
CN103094921A (zh) * | 2013-01-06 | 2013-05-08 | 宁夏电力公司电力科学研究院 | 直驱式风力发电机组机电暂态建模方法及机电暂态模型 |
CN103746404A (zh) * | 2014-01-24 | 2014-04-23 | 国家电网公司 | 一种风光火打捆直流外送系统稳定性评价方法 |
CN103904678A (zh) * | 2014-03-31 | 2014-07-02 | 中国南方电网有限责任公司电网技术研究中心 | 高压直流输电分段变速率低压限流单元的控制方法 |
CN106356885A (zh) * | 2016-09-19 | 2017-01-25 | 国家电网公司 | 一种风电并网系统暂态功角前二摆稳定性判断方法 |
CN107317354A (zh) * | 2017-06-30 | 2017-11-03 | 天津大学 | 一种含风电场的多机系统暂态功角稳定性分析方法 |
Non-Patent Citations (3)
Title |
---|
屠竞哲 等: "多送出直流系统送端故障引发稳定破坏机理分析", 《电力系统自动化》 * |
李佳 等: "感应电动机模型和机械转矩参数对暂态电压稳定评估的影响", 《电力系统保护与控制》 * |
林俐 等: "基于扩展等面积定则的含大规模风电场电力系统暂态稳定性分析", 《电力系统保护与控制》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110890751A (zh) * | 2019-10-18 | 2020-03-17 | 东南大学 | 一种适用含高比例风电电力系统的暂态稳定特征保全方法 |
CN110890751B (zh) * | 2019-10-18 | 2022-07-05 | 东南大学 | 一种适用含高比例风电电力系统的暂态稳定特征保全方法 |
CN112054550A (zh) * | 2020-09-11 | 2020-12-08 | 武汉大学 | 一种多回特高压直流馈入的交流系统的直流功率互济方法 |
CN112054550B (zh) * | 2020-09-11 | 2022-04-15 | 武汉大学 | 一种多回特高压直流馈入的交流系统的直流功率互济方法 |
CN114221339A (zh) * | 2021-12-21 | 2022-03-22 | 上海电力大学 | 一种高风电占比系统的暂态稳定性判断方法 |
CN114221339B (zh) * | 2021-12-21 | 2024-03-26 | 上海电力大学 | 一种高风电占比系统的暂态稳定性判断方法 |
CN114268100A (zh) * | 2022-03-01 | 2022-04-01 | 国网江西省电力有限公司电力科学研究院 | 一种交直流混联电力系统暂态稳定性能分析方法及装置 |
CN114268100B (zh) * | 2022-03-01 | 2022-07-12 | 国网江西省电力有限公司电力科学研究院 | 一种交直流混联电力系统暂态稳定性能分析方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108767898B (zh) | 2021-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108767898A (zh) | 一种考虑风电与直流综合作用的电力系统暂态稳定分析方法 | |
CN106410870B (zh) | 一种海上风电接入多端柔性直流输电系统的频率控制方法 | |
CN101937483B (zh) | 基于桨距角动作情况的双馈机组风电场动态等值建模方法 | |
Zhang et al. | Impacts of PLL on the DFIG-based WTG's electromechanical response under transient conditions: analysis and modeling | |
CN103715696B (zh) | 一种基于多种风电机组共济的风电场无功电压集群控制方法 | |
Zhang et al. | Control of DFIG-based wind farms for power network frequency support | |
Zhai et al. | Reactive power control strategy of DFIG wind farms for regulating voltage of power grid | |
CN106374524A (zh) | 一种upfc对异步风力发电机系统次同步谐振的抑制方法 | |
Siswanto et al. | Steady state stability limit assessment when wind turbine penetrated to the systems using REI approach | |
Shi et al. | Modelling analysis of transient stability simulation with high penetration of grid‐connected wind farms of DFIG type | |
CN109713661A (zh) | 风电场接入对多机系统故障极限切除时间影响的分析方法 | |
Krpan et al. | Towards the new low-order system frequency response model of power systems with high penetration of variable-speed wind turbine generators | |
Sravanthi et al. | Critical clearing time and transient stability analysis of SCIG based wind farm with STATCOM | |
Lin et al. | Transient stability mechanism of DFIG wind farm and grid-connected power system | |
Wang et al. | Wind farm model with DFIG for small signal stability study | |
Chandra et al. | Small signal stability of power system with SCIG, DFIG wind turbines | |
CN104135207B (zh) | 一种大型双馈风力发电机并网联合优化控制方法 | |
Li | Power flow modeling considering detailed constraints of the DFIGs and collector networks based on 3-layer BFS and convergence improvement | |
Elkomy et al. | Enhancement of wind energy conversion systems active and reactive power control via flywheel energy storage systems integration | |
Mensou et al. | A robust speed control of a doubly fed induction generator using in WECS by the nonlinear Backstepping controller | |
Ali | Step towards enriching frequency support from wind-driven permanent-magnet synchronous generator for power system stability | |
Kawale et al. | Comparative study of converter topologies used for PMSG based wind power generation | |
Ebrahimi et al. | A robust fuzzy-based control technique for wind farm transient voltage stability using SVC and STATCOM: comparison study | |
Lin | Transient stability study of power system with large-scale wind farm integration | |
Shao et al. | Digital Mirroring System for DFIG Wind Turbine Generator with Parallel Converters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |