发明内容
本发明是为避免上述现有技术所存在的不足之处,提供一种综合考虑负荷与风电时序特性的随机生产模拟方法,以期能快速、合理地研究风电时序特性和发电机组处于不同运行状态下各种运行因素对随机生产模拟的影响,从而能建立起更加全面的数学模型,提高计算效率。
本发明为达到上述发明目的采用如下技术方案:
本发明综合考虑负荷与风电时序特性的随机生产模拟方法,是应用于电力***中,其特点是按如下步骤进行:
步骤1、获取研究周期T内电力***的发电机组数据、风电资源数据和时序负荷数据,由所述时序负荷数据获得原始时序负荷曲线,并记为Load1;
步骤2、建立四种机组模型,包括:基荷机组模型、启停调峰机组模型、压负荷调峰机组模型和两分段启停调峰机组模型;
步骤3、由所述发电机组数据和时序负荷数据,初步修正所述原始时序负荷曲线Load1,获得初步修正时序负荷曲线,记为Load2;
步骤4、设定风电修正负荷次数为n,并初始化n=1;
步骤5、由所述发电机组数据获得供热机组强迫出力,根据所述风电资源数据和供热机组强迫出力,修正所述初步修正时序负荷曲线Load2,获得第n次风电接纳量Wn和第n次时序净负荷曲线Load3n;
步骤6、将所述第n次时序净负荷曲线Load3n转化为第n次等效持续负荷曲线,并求取所述第n次等效持续负荷曲线上,负荷值为X处的初始累积概率P0,n(X)、初始累积频率F0,n(X)、初始电量不足期望值EENS0,n(X)和初始最大增长速率
步骤7、根据所述电力***的运行需求确定所述电力***的初始旋转备用容量和初始爬坡容量值
步骤8、由所述发电机组数据获得总投运次数I;设定投运机组次数为i,并初始化i=1;
步骤9、根据电力***的约束条件,以最低发电成本为目标,判断是否需要重新修正风电,若需要,则将n+1赋值给n,并返回步骤5执行,若不需要,则确定第i次投运的发电机组号;
步骤10、根据所述第i次投运的发电机组类型及带基荷的状况,确定第i次投运的发电机组所属的机组模型;
当i=1时,由所述初始累积概率P0,n(X)和所述初始累积频率F0,n(X)获得第i次投运发电机组时所述电力***需要机组投运的转移率和不需要机组投运的转移率
当i≠1时,由第i-1次投运发电机组后的累积概率Pi-1,n(X)和第i-1次累计频率Fi-1,n(X)获得第i次投运发电机组时所述电力***需要机组投运的转移率和不需要机组投运的转移率
步骤11、根据所述第i次投运的发电机组号,修正获得第i次投运发电机组后所述电力***的旋转备用容量和爬坡容量
步骤12、将所述第i次投运的发电机组所属的机组模型转化为等效两状态机组模型或等效三状态机组模型;
步骤13、将所述等效两状态机组模型或等效三状态机组模型进行随机生产模拟,获得第i次投运发电机组后的累积概率Pi,n(X)、累积频率Fi,n(X)、电量不足期望值EENSi,n(X)和最大增长速率
步骤14、累计第i次投运发电机后的总投运容量,并判断i<I是否成立;若成立,则将i+1赋值给i;并返回步骤9;否则,表示完成随机生产模拟。
本发明综合考虑负荷与风电时序特性的随机生产模拟方法的特点也在于:
步骤1中所述发电机组包括火电发电机组、水电发电机组和核电发电机组;所述发电机组数据包括机组的类型、台数、容量、出力上下限、平均无故障时间MTTF、平均修复时间MTTR以及耗量微增率;所述风电资源数据包括风电机组出力及其对应的时段;所述时序负荷数据为按时间顺序排列的逐小时负荷。
步骤2中是根据随机生产模拟过程中机组的不同运行状态、带基荷的状况以及分段机组容量来建立四种机组模型;
所述基荷机组模型:表示一台以额定容量在基荷运行的不分段连续运行机组;
所述启停调峰机组模型:表示一台以额定容量在基荷以上做调峰运行的不分段间歇运行机组;
所述压负荷调峰机组模型:表示一台第一分段以最小技术出力在基荷运行,第二分段在基荷以上做调峰运行的分段连续运行机组;
所述两分段启停调峰机组模型:表示一台两个分段均在基荷以上做调峰运行的分段间歇运行机组。
步骤5中所述第n次风电接纳量Wn按如下步骤获得:
步骤a、设定变量t,并初始化t=1;将初步修正时序负荷曲线Load2在t时刻的负荷值记为Load2t;
步骤b、判断n=1是否成立,若成立,则转步骤c;否则,转步骤d;
步骤c、判断t时刻是否满足:Load2t-强迫出力>=风电机组发电量;
若满足,则令t时刻的风电接纳量=风电机组发电量;否则,令t时刻风电接纳量=Load2t-强迫出力;转步骤f;
步骤d、以λ为弃风比率;判断t时刻是否满足:Load2t-(强迫出力+λ×风电装机容量×n)>=风电机组发电量;
若满足,则令t时刻风电接纳量=风电机组发电量;否则,转步骤e;
步骤e、判断Load2t-(强迫出力+λ×风电装机容量×n)>0是否成立,若成立,则令t时刻风电接纳量=Load2t-(强迫出力+λ×风电装机容量×n);否则,令t时刻风电接纳量=风电机组发电量;
步骤f、判断t=T是否成立,若成立,则表示获得第n次风电接纳量Wn;否则,将t+1赋值给t,并返回步骤b。
步骤9中是按如下步骤判断是否需要重新修正风电:
步骤(1)、判断i>1且第i次投运发电机后的总投运容量大于所述电力***的最大负荷是否成立,若成立,则选择耗量微增率最小的机组或机组分段进行投运,并转步骤(7);否则,转步骤(2);
步骤(2)、判断所有未投运的机组中是否还有强迫出力机组或机组分段,若有,则选取耗量微增率最小的强迫出力机组或机组分段进行投运,并转步骤(7);否则,转步骤(3);
步骤(3)、判断是否满足所述电力***的旋转备用约束和爬坡约束条件,若满足,则转步骤(5);否则,转步骤(4);
步骤(4)、判断所有未投运的机组中是否存在分段机组的第一分段,若存在,则选取耗量微增率最小的第一分段进行投运,并转步骤(7);否则,表示需要重新修正风电;
步骤(5)、查找所有未投运的机组中耗量微增率最小的机组或机组分段,判断该机组或机组分段是否为分段机组的第二分段;若是,则转步骤(6);否则,转步骤(7);
步骤(6)、对所述电力***的旋转备用容量和爬坡容量进行预修正,并判断是否满足预修正后的旋转备用约束和爬坡约束条件,若满足,则投入分段机组的第二分段,并转步骤(7);否则,则在排除最小比耗量机组为分段机组的第二分段的情况下,转步骤(5);
步骤(7)、表示不需要重新修正风电。
步骤9中所述电力***的约束条件包括如式(1)和式(2)所示的旋转备用约束和爬坡约束的:
式(1)和式(2)中:表示第i次投运发电机组后所述电力***的旋转备用容量;SPR0为***的旋转备用容量需求;表示第i次投运发电机组后所述电力***的增荷容量;t0为要求旋转备用容量投入的时间;为事故旋转备用容量需求;为负荷增长旋转备用容量需求;Xi为第i次投运的发电机组在等效持续负荷曲线上的投运点;表示第i次投运的发电机组在投运点Xi处最大增荷速率;Δt表示爬坡时间。
步骤10中,所述需要机组投运的转移率和不需要机组投运的转移率按利用式(3)和式(4)进行计算:
式(3)和式(4)中:Fi-1,n(Xi)表示第i-1次投运发电机组后负荷值X=Xi处的累积概率;Pi-1,n(Xi)表示第i-1次投运发电机组后负荷值X=Xi处的累积频率。
步骤11中,按如下情况修正所述旋转备用容量和爬坡容量
(a)如果第i次投运的发电机组为分段机组的第一分段,其容量为Cd i,则利用式(5)和式(6)进行修正:
(b)如果第i次投运的发电机组为分段机组的第二分段,其容量为Xd i,则利用式(7)和式(8)进行修正:
(c)如果第i次投运的发电机组为不分段机组,则利用式(9)和式(10)进行修正:
式(5)-式(10)中,Ce i表示第i次投运的发电机组的额定容量;表示第i-1次投运发电机组后电力***的旋转备用容量;表示第i-1次投运发电机组后***的爬坡容量。
步骤13中所述的累积概率Pi,n(X)、累积频率Fi,n(X)、电量不足期望值EENSi,n(X)和最大增长速率是按如下情况获得:
1)若第i次投运的发电机组为不分段机组或分段机组的第一分段,则利用式(12)、式(13)和式(14)获得累积概率Pi,n(X)、累积频率Fi,n(X)和电量不足期望值EENSi,n(X):
Pi,n(X)=Pi-1,n(X)·(1-FORi)+Pi-1,n(X-Ci)·FORi (12)
Fi,n(X)=Fi-1,n(X)·(1-FORi)+Fi-1,n(X-Ci)·FORi+ue i·FORi·[Pi-1,n(X-Ci)-Pi-1,n(X)] (13)
EENSi,n(X)=EENSi-1,n(X)·(1-FORi)+EENSi-1,n(X-Ci)·FORi (14)
式(12)、式(13)和式(14)中:X表示所述等效持续负荷曲线上的负荷值;Pi,n(X)、Fi,n(X)和EENSi,n(X)分别表示第i次机组投运后等效持续负荷曲线上负荷值为X处的累积概率、累积频率和电量不足期望值;Pi-1,n(X)、Fi-1,n(X)和EENSi-1,n(X)分别表示第i-1次机组投运后等效持续负荷曲线上负荷值为X处的累积概率、累积频率和电量不足期望值;FORi表示第i次投运的发电机组的等效两状态机组出力为零的概率;ue i表示所述等效两状态机组的等效修复率;Ci表示第i次投运的发电机组或机组分段的容量,若第i次投运的发电机组为不分段机组,则Ci=Ce i,若第i次投运的发电机组为分段机组的第一分段,则Ci=Cd i。
若第i次投运的发电机组为分段机组的第二分段,则利用式(15)、式(16)和式(17)进行反卷积运算,获得反卷积运算后的累积概率P1 i(X)、累积频率F1 i(X)和电量不足期望值EENS1 i(X):再利用式(18)、式(19)和式(20),获得累积概率Pi,n(X)、累积频率Fi,n(X)和电量不足期望值EENSi,n(X):
P1 i(X)=[Pi-1,n(X)-P1 i(X-Cd i)·FORi]/(1-FORi) (15)
F1 i(X)=[Fi-1.n(X)-F1 i(X-Cd i)·FORi-ue i·FORi·(P1 i(X-Cd i)-P1 i(X))]/(1-FORi) (16)
EENS1 i(X)=[EENSi-1,n(X)-EENS1 i(X-Cd i)]/(1-FORi) (17)
Pi,n(X)=P1 i(X)·AVi+P1 i(X-Xd i)·DFORi+P1 i(X-Ce i)·FORi (18)
Fi,n(X)=F1 i(X)·AVi+F1 i(X-Xd i)·DFORi+F1 i(X-Ce i)·FORi+AVi·λT i·[P1 i(X-Ce i)-P1 i(X)]+AVi·ρ- i,n·[P1 i(X-Xd i)-P1 i(X)]+DFORi·λD i·[P1 i(X-Ce i)-P1 i(X-Xd i)]
(19)
EENSi,n(X)=EENS1 i(X)·AVi+EENS1 i(X-Xd i)·DFORi+EENS1 i(X-Ce i)·FORi (20)
式(15)—式(20)中:AVi,DFORi,FORi分别表示所述等效三状态机组额定出力概率、降额出力概率和零出力概率;λT i和λD i分别表示所述等效三状态机组中额定状态和降额状态的强迫停运率,P1 i(X),F1 i(X),EENS1 i(X)分别表示将一段反卷积后等效持续负荷曲线上负荷值为X处的累积概率、累积频率和电量不足期望值。
与已有技术相比,本发明有益效果体现在:
1、本发明通过建立一种综合考虑负荷与风电时序特性的随机生产模拟方法,充分考虑了风电的时序特性和机组处在不同运行状态下各种运行因素对随机生产模拟的影响,是一种高效、全面的电力***随机生产模拟的方法。与传统随机生产模拟机组模型不同,本发明根据投运机组类型及带负荷位置确定机组模型,同一机组模型因带负荷位置不同而改变,有效反映出所有机组的运行情况;另外,本发明可有效考虑***备用和爬坡约束,以及机组最小开停机限制、启动失败和爬坡时间约束等因素的影响,提高了计算速度,具有重要的工程应用价值。
2、本发明采用时序负荷曲线,并提取时序负荷的相关信息,可有效反映负荷时序改变对整个随机生产模拟过程的影响。
3、本发明根据随机生产模拟过程中机组的不同运行状态(如降额、间歇、连续)、带负荷位置(是否带基荷)及分段机组,建立了四种多状态机组模型,可考虑机组启动失败率、最小开停机时间限制和爬坡速率限制,充分考虑了机组自身结构和各种运行因素的影响,提高了解析法数学模型的合理性。
4、本发明重点考虑机组发电时段中出现强迫停运的概率,采用***需要机组投入的情况下机组故障的条件强迫停运率CFOR进行卷积运算,其中条件强迫停运率CFOR不仅与已投运机组的停运特性有关,也与负荷频率转移特性有关,因此动态计算条件强迫停运率CFOR,能有效去除在***不需要机组投运时故障对卷积计算的影响。
5、本发明着重考虑风电的反调峰性和不可调度性,提出了大规模风电并网时谷荷比例弃风的方法,充分保留了风电出力的时序特性。
6、本发明提供了***的旋转备用和爬坡约束条件,并在每次投运机组或机组分段后进行校验,保证了旋转备用和爬坡约束能时刻满足,提高了***的可靠性。
7、本发明将四种多状态机组模型转化为等效两状态机组模型和等效三状态机组模型,简化了机组模型的卷积运算,提高了计算效率。
具体实施方式
如图1所示,本实施例中,综合考虑负荷与风电时序特性的随机生产模拟方法按如下步骤进行:
步骤1、获取研究周期T内电力***的发电机组数据、风电资源数据和时序负荷数据,由时序负荷数据获得原始时序负荷曲线,并记为Load1;
步骤1中的发电机组包括火电发电机组、水电发电机组和核电发电机组;发电机组数据包括机组的类型(分段情况)、台数、容量、出力上下限、平均无故障时间MTTF、平均修复时间MTTR以及耗量微增率;风电资源数据包括风电机组出力及其对应的时段;时序负荷数据为按时间顺序排列的逐小时负荷。所有含时间信息的数值计算时需要保证单位统一(如可统一以小时为单位)。
步骤2、建立四种机组模型,包括:基荷机组模型、启停调峰机组模型、压负荷调峰机组模型和两分段启停调峰机组模型;
步骤2是根据随机生产模拟过程中机组的不同运行状态(如降额、间歇、连续)、带基荷的状况(是否带基荷)以及分段机组容量(假设分段机组第一分段容量等于降额状态出力)来建立四种机组模型;
基荷机组模型(两状态模型)如图2所示:表示一台以额定容量在基荷运行的不分段连续运行机组;其停运概率与传统生产模拟的模型无异,可直接采用强迫停运率FOR(即机组零状态概率)计算。
启停调峰机组模型(五状态模型)如图3所示:表示一台以额定容量在基荷以上做调峰运行的不分段间歇运行机组;需转化为等效两状态机组模型进行计算。
压负荷调峰机组模型(四状态模型)如图4所示:表示一台第一分段以最小技术出力在基荷运行,第二分段在基荷以上做调峰运行的分段连续运行机组;在***需要机组投运的情况下,需转化为等效三状态机组模型进行计算。
两分段启停调峰机组模型(七状态模型)如图5所示:表示一台两个分段均在基荷以上做调峰运行的分段间歇运行机组。同样需转化为等效三状态机组模型进行计算。
步骤3、由发电机组数据和时序负荷数据,初步修正原始时序负荷曲线Load1,获得初步修正时序负荷曲线,记为Load2;
步骤3中的初步修正原始时序负荷曲线Load1主要为:以时间顺序排列的逐小时负荷数据,先按负的负荷模型处理承担基荷的核电和径流水电;针对水电比重较小的***,再按逐次切负荷方法,安排库容式水电进行削峰,降低负荷曲线的峰谷差;若水电所占比重较大,则将库容式水电机组在负荷曲线上寻求最佳工作位置并进行卷积运算。
逐次切负荷法:根据日负荷图得到最大负荷;用最大负荷减去该机组最大允许出力以定出该机组最初工作位置;如图6所示,对于该机组如果阴影部分的面积大于给定值,则抬高工作位置;反之则降低工作位置。让工作位置以上的负荷都由该机组带,并将大于该机组最大允许出力时段的出力取其最大允许出力值。
如图6所示,抬高或降低工作位置的步长取Δx=(E-Eg)/T,T为模拟时间长度,一般取24h,当|E-Eg|≤ε时满足迭代要求,ε为设定的精度。
步骤4、设定风电修正负荷次数为n,并初始化n=1;若风电不需要进行多次谷荷弃风,则n=1不会改变;若风电需进行多次谷荷弃风,则每进行一次谷荷弃风,n的值就会增加1,但由于多次弃风的次数无法提前确定,故n的最大值也无法提前知晓,只有完成整个随机生产过程之后才能确定(需要注意的是,只要风电修正次数改变,所有的机组需要重新安排投运,故机组参数和***约束条件参数均为第n次风电修正负荷后的值)。
步骤5、由发电机组数据获得供热机组强迫出力,根据风电资源数据和供热机组强迫出力(即所有机组最小技术出力之和),修正初步修正时序负荷曲线Load2,获得第n次风电接纳量Wn和第n次时序净负荷曲线Load3n;
步骤5中的第n次风电接纳量Wn按如下步骤获得:
步骤a、设定变量t,并初始化t=1;将初步修正时序负荷曲线Load2在t时刻的负荷值记为Load2t;
步骤b、判断n=1是否成立,若成立,则转步骤c;否则,转步骤d;
步骤c、判断t时刻是否满足:Load2t-强迫出力>=风电机组发电量;
若满足,则令t时刻的风电接纳量=风电机组发电量;否则,令t时刻风电接纳量=Load2t-强迫出力;转步骤f;
步骤d、以λ为弃风比率;判断t时刻是否满足:Load2t-(强迫出力+λ×风电装机容量×n)>=风电机组发电量;
若满足,则令t时刻风电接纳量=风电机组发电量;否则,转步骤e;
步骤e、判断Load2t-(强迫出力+λ×风电装机容量×n)>0是否成立,若成立,则令t时刻风电接纳量=Load2t-(强迫出力+λ×风电装机容量×n);否则,令t时刻风电接纳量=风电机组发电量;
步骤f、判断t=T是否成立,若成立,则表示获得第n次风电接纳量Wn;否则,将t+1赋值给t,并返回步骤b。
以T=24为例,如图7所示,当n=1时,经过初始谷荷弃风,可获得周期T内的初始弃风量和初始风电接纳接纳量W1;当n>1时,经过多次谷荷弃风直到满足约束条件,最终可获得周期T内的弃风量和风电接纳接纳量Wn;
步骤6、将第n次时序净负荷曲线Load3n转化为第n次等效持续负荷曲线,并求取第n次等效持续负荷曲线上,负荷值为X处的初始累积概率P0,n(X)、初始累积频率F0,n(X)、初始电量不足期望值EENS0,n(X)和初始最大增长速率
步骤6中的负荷值为X处的初始累积概率P0,n(X)、初始累积频率F0,n(X)、初始电量不足期望值EENS0,n(X)和初始最大增长速率是按如下步骤求取:
如图8,取从某日正午12点开始连续24h数据:
初始累积概率P0,n(X)、累积频率F0,n(X)和负荷最大增长速率的计算式如下:
式(1)-式(3)中:Pt(X),Ft(X)和为各时刻负荷值为X处的累积概率,累积频率和增长速率,其计算式如下:
式(4)-式(6)中:xt,xt+1为t和t+1时刻的负荷值。
初始电量不足期望值EENS0,n(X)为图8中阴影部分面积;本发明中负荷值X代表一系列负荷点,其值从0开始到((最大负荷+机组总容量)/ΔX+1)·ΔX为止,步长为ΔX(一般取机组容量最大公约数)。
步骤7、根据电力***的运行需求确定(第n次风电修正负荷后)电力***的初始旋转备用容量和初始爬坡容量值以及爬坡容量需求、事故备用容量需求和负荷增长备用容量需求;
步骤8、由发电机组数据获得总投运次数I;设定投运机组次数为i,并初始化i=1;
步骤8中的总投运次数I=分段机组台数×2+不分段机组台数。
步骤9、根据电力***的约束条件,以最低发电成本为目标,判断是否需要重新修正风电,若需要,则将n+1赋值给n,并返回步骤5执行,若不需要,则确定第i次投运的发电机组号;
步骤9是按如下步骤判断是否需要重新修正风电:
步骤(1)、判断i>1且第i次投运发电机后的总投运容量大于电力***的最大负荷是否成立,若成立,则选择耗量微增率最小的机组或机组分段进行投运,并转步骤(7);否则,转步骤(2);
步骤(2)、判断所有未投运的机组中是否还有强迫出力机组或机组分段,若有,则选取耗量微增率最小的强迫出力机组或机组分段进行投运,并转步骤(7);否则,转步骤(3);
步骤(3)、判断是否满足电力***的旋转备用约束和爬坡约束条件,若满足,则转步骤(5);否则,转步骤(4);
步骤(4)、判断所有未投运的机组中是否存在分段机组的第一分段,若存在,则选取耗量微增率最小的第一分段进行投运,并转步骤(7);否则,表示需要重新修正风电;
步骤(5)、查找所有未投运的机组中耗量微增率最小的机组或机组分段,判断该机组或机组分段是否为分段机组的第二分段;若是,则转步骤(6);否则,转步骤(7);
步骤(6)、对电力***的旋转备用容量和爬坡容量进行预修正,并判断是否满足预修正后的旋转备用约束和爬坡约束条件,若满足,则投入分段机组的第二分段,并转步骤(7);否则,则在排除最小比耗量机组为分段机组的第二分段的情况下,转步骤(5);
步骤(7)、表示不需要重新修正风电。
步骤9中电力***的约束条件包括如式(7)和式(8)所示的旋转备用约束和爬坡约束:
式(7)和式(8)中:表示(第n次风电修正负荷后)第i次投运发电机组后电力***的旋转备用容量,常规机组都没有投运时,其值为可快速投运的冷备用机组容量;SPR0为***的旋转备用容量需求,通常为一给定值;表示(第n次风电修正负荷后)第i次投运发电机组后电力***的增荷容量,若投入不分段机组,增荷容量不变,如投入分段机组,则需要对增荷容量进行修正;t0为要求旋转备用容量投入的时间;为事故旋转备用容量需求,即***要求其设置的事故旋转备用容量在t0小时内全部到达(若取t0=0.2h,即要求12min内***能够提供所需旋转备用容量值);为负荷增长旋转备用容量需求,负荷备用需要的平均增荷速率为式(8)不等式右边整体代表***的爬坡容量需求,由于可见其值为一变量;Xi为第i次投运的发电机组在等效持续负荷曲线上的投运点;表示(第n次风电修正负荷后)第i次投运的发电机组在投运点Xi处最大增荷速率;Δt表示考虑爬坡速率约束的爬坡时间,可取Δt=0.25h(即15min)。
除考虑***旋转备用约束和爬坡约束条件外,还考虑总投运容量与最大负荷关系,机组的随机故障、机组启动失败概率PS、机组最小开停机时间限制和机组爬坡速率的限制TR(以启停调峰机组模型为例,即从零出力状态到满额出力状态所需的爬坡时间)。
步骤10、根据第i次投运的发电机组类型及带基荷的状况,确定第i次投运的发电机组所属的机组模型;
当i=1时,由初始累积概率P0,n(X)和初始累积频率F0,n(X)获得(第n次风电修正负荷后)第i次投运发电机组时电力***需要机组投运的转移率和不需要机组投运的转移率
当i≠1时,由第i-1次投运发电机组后的累积概率Pi-1,n(X)和第i-1次累计频率Fi-1,n(X)获得(第n次风电修正负荷后)第i次投运发电机组时电力***需要机组投运的转移率和不需要机组投运的转移率
步骤10中的需要机组投运的转移率和不需要机组投运的转移率按利用式(9)和式(10)进行计算:
式(9)和式(10)中:Fi-1,n(Xi)表示(第n次风电修正负荷后)第i-1次投运发电机组后负荷值X=Xi处的累积概率;Pi-1,n(Xi)表示(第n次风电修正负荷后)第i-1次投运发电机组后负荷值X=Xi处的累积频率。
设定发电机组最小连续运行时间和最小经济停机时间假设整个***的停运持续时间和正常运行时间服从电力***需要机组投运的转移率(又称“需求率”)和不需要机组投运的转移率(又称“不需求率”)的指数分布,故:
式(11)和式(12)中:pu i是第i次投运的发电机组的运行时间大于最小连续运行时间的概率;pd i是第i次投运的发电机组的停运时间大于最小经济停机时间的概率。
机组的开停机限制约束:需将启停调峰机组和两分段启停调峰机组模型中的需求率和不需求按下式进行修正:
此时修正后的需求率和不需求在原有的基础上打了折扣,所得的机组模型为考虑机组最小开停机限制后,机组在整个随机生产模拟过程中的状态空间图,各部分概率即为机组在整个生产模拟过程中所处在该状态的概率;修正后的需求率则表示加入开停机时间限制后,机组从备用状态的向爬坡状态的实际转移率。
步骤11、根据第i次投运的发电机组号,修正获得(第n次风电修正负荷后)第i次投运发电机组后电力***的旋转备用容量和爬坡容量
步骤11中是按如下情况修正旋转备用容量和爬坡容量
(a)如果第i次投运的发电机组为分段机组的第一分段,其容量为Cd i,此时机组1段投运,第2分段处在旋转备用状态,备用容量增加量,爬坡容量也增加,增加的爬坡容量最多为2段容量,则利用式(15)和式(16)进行修正:
(b)如果第i次投运的发电机组为分段机组的第二分段,其容量为Xd i,2段退出旋转备用状态,备用容量减少,爬坡容量也随之减少,减少的爬坡容量最多为2段容量,则利用式(17)和式(18)进行修正:
(c)如果第i次投运的发电机组为不分段机组,机组投运对***的备用容量和爬坡容量没有影响,则利用式(19)和式(20)进行修正:
式(15)-式(20)中,Ce i表示第i次投运的发电机组的额定容量;表示(第n次风电修正负荷后)第i-1次投运发电机组后电力***的旋转备用容量;表示(第n次风电修正负荷后)第i-1次投运发电机组后***的爬坡容量。
步骤12、将第i次投运的发电机组所属的机组模型转化为等效两状态机组模型,如图9所示,或等效三状态机组模型,如图10所示;
基荷机组模型和启停调峰机组模型需转化为等效两状态机组模型,压负荷调峰机组模型和两分段启停调峰机组模型需转化为等效三状态机组模型;其等效两状态机组模型或等效三状态机组模型的状态概率计算如下:
1)基荷机组模型(两状态模型)可直接采用强迫停运率FOR(即机组零状态概率)计算。即:
CFOR=FOR (21)
2)启停调峰机组模型(五状态模型)需转化为等效两状态机组模型进行计算。停运概率采用条件强迫停运率CFOR计算,CFOR一般指***需要机组投运的情况下机组故障的概率;需要注意的是,图3考虑了爬坡时间约束,此时CFOR表示***需要机组投运的情况下机组出力为零的概率,图中爬坡状态为机组刚启动状态,其出力为零。故:
式(22)中:P3,P4,P5分别表示图中机组模型中状态3,4,5的概率。
3)压负荷调峰机组模型(四状态模型)在***需要机组投运的情况下,需转化为等效三状态机组模型进行计算。等效三状态的三个状态分别为:机组额定出力状态概率AV,机组降额出力状态概率DFOR和机组零状态概率FOR,对应计算式如下:
式(23)中:P1,P2,P3,P4分别表示图中机组模型中状态1,2,3,4的概率。
4)两分段启停调峰机组模型(七状态模型)同样等效为三状态机组进行计算;图4中爬坡1段为刚启动时刻,出力近似为零,爬坡3段为刚从降额开始增加出力时刻,出力近似为降额出力,其等效三状态状态概率计算如下:
式(24)中:P3,P4,P5,P6,P7分别表示图中机组模型中状态3,4,5,6,7的概率。
步骤13、将等效两状态机组模型或等效三状态机组模型进行随机生产模拟,获得(第n次风电修正负荷后)第i次投运发电机组后的累积概率Pi,n(X)、累积频率Fi,n(X)、电量不足期望值EENSi,n(X)和最大增长速率并计算获得第i次投运的发电机组的发电量、发电成本以及机组的冷、热启动次数;
步骤13中的累积概率Pi,n(X)、累积频率Fi,n(X)、电量不足期望值EENSi,n(X)和最大增长速率是按如下情况获得:
1)若第i次投运的发电机组为不分段机组或分段机组的第一分段,则利用式(26)、式(27)和式(28)获得累积概率Pi,n(X)、累积频率Fi,n(X)和电量不足期望值EENSi,n(X):
Pi,n(X)=Pi-1,n(X)·(1-FORi)+Pi-1,n(X-Ci)·FORi (26)
Fi,n(X)=Fi-1,n(X)·(1-FORi)+Fi-1,n(X-Ci)·FORi+ue i·FORi·[Pi-1,n(X-Ci)-Pi-1,n(X)] (27)
EENSi,n(X)=EENSi-1,n(X)·(1-FORi)+EENSi-1,n(X-Ci)·FORi (28)
式(26)、式(27)和式(28)中:X表示等效持续负荷曲线上的负荷值,以ΔX(一般取所有发电机组容量的最大公约数)为间隔,X=0,ΔX,2·ΔX…(最大负荷+Pi(X)/ΔX+1)·ΔX;Pi,n(X)、Fi,n(X)和EENSi,n(X)分别表示(第n次风电修正负荷后)第i次机组投运后等效持续负荷曲线上负荷值为X处的累积概率、累积频率和电量不足期望值;Pi-1,n(X)、Fi-1,n(X)和EENSi-1,n(X)分别表示(第n次风电修正负荷后)第i-1次机组投运后等效持续负荷曲线上负荷值为X处的累积概率、累积频率和电量不足期望值;FORi表示第i次投运的发电机组的等效两状态机组出力为零的概率;ue i表示等效两状态机组的等效修复率;Ci表示第i次投运的发电机组或机组分段的容量,若第i次投运的发电机组为不分段机组,则Ci=Ce i,若第i次投运的发电机组为分段机组的第一分段,则Ci=Cd i。
若第i次投运的发电机组为分段机组的第二分段,则利用式(29)、式(30)和式(31)进行反卷积运算,获得反卷积运算后的累积概率P1 i(X)、累积频率F1 i(X)和电量不足期望值EENS1 i(X):再利用式(32)、式(33)和式(34),获得累积概率Pi,n(X)、累积频率Fi,n(X)和电量不足期望值EENSi,n(X):
P1 i(X)=[Pi-1,n(X)-P1 i(X-Cd i)·FORi]/(1-FORi) (29)
F1 i(X)=[Fi-1.n(X)-F1 i(X-Cd i)·FORi-ue i·FORi·(P1 i(X-Cd i)-P1 i(X))]/(1-FORi) (30)
EENS1 i(X)=[EENSi-1,n(X)-EENS1 i(X-Cd i)]/(1-FORi) (31)
Pi,n(X)=P1 i(X)·AVi+P1 i(X-Xd i)·DFORi+P1 i(X-Ce i)·FORi (32)
Fi,n(X)=F1 i(X)·AVi+F1 i(X-Xd i)·DFORi+F1 i(X-Ce i)·FORi+AVi·λT i·[P1 i(X-Ce i)-P1 i(X)]+AVi·ρ- i,n·[P1 i(X-Xd i)-P1 i(X)]+DFORi·λD i·[P1 i(X-Ce i)-P1 i(X-Xd i)]
(33)
EENSi,n(X)=EENS1 i(X)·AVi+EENS1 i(X-Xd i)·DFORi+EENS1 i(X-Ce i)·FORi (34)
式(29)—式(34)中:AVi,DFORi,FORi分别表示等效三状态机组额定出力概率、降额出力概率和零出力概率;λT i和λD i分别表示等效三状态机组中额定状态和降额状态的强迫停运率,P1 i(X),F1 i(X),EENS1 i(X)分别表示将一段反卷积后等效持续负荷曲线上负荷值为X处的累积概率、累积频率和电量不足期望值。
步骤14、累计第i次投运发电机后的总投运容量,并判断i<I是否成立;若成立,则将i+1赋值给i;并返回步骤9;否则,表示完成随机生产模拟。
最终随机生产模拟结果中,可计算处理得到风电接纳量、弃风率,并计算得到生产模拟周期内冷、热启动次数、启动成本和冷、热启动成本以及可靠性指标EENS和LOLP;另外,还可以研究机组自身约束条件和***约束条件对随机生产模拟的影响。