[发明专利]一种风电并网的火电系统动态调度方法

权利要求书

1.一种风电并网的火电系统动态调度方法，其特征在于，具体按照如下步骤实施：

步骤1、相关参数的设置和获取

1.1）设置算法运行的相关参数

设定最大迭代次数为iter_max，iter_max为整数，迭代次数初始值iter＝1；学习因子设置为c₁、c₂，c₁和c₂取为[0,3]之间的自然数，r₁、r₂是[0,1]之间的随机数；调度时段为h，h＝1,2,...,H，智能算法的群体规模数为M，M为大于1的整数；

PID控制器的参数分别取比例增益K_P＝-1，积分项系数T_i＝iter_max/20，微分项系数T_d＝0，调节因子σ＝1；

1.2）获得风电相关参数

包括风电厂中的风力发电机数目为L，L为不小于1的整数；第l台风力发电机组的额定输出功率为P_R，l，l=1,2,,...,L；风力发电机理论设计额定风速为v_R，l，切入风速为v_CI，l，切出风速为v_CO，l；

1.3）获得火电相关参数

根据系统运行的需求，确定以下火电机组的已知参数：机组数目为N，N为不小于1的整数；第i台火电机组的最小输出功率为P_imin和最大输出功率P_imax，i＝1,2,...,N，有功输出功率下降速率ζ_idown和上升速率ζ_iup；火电机组消耗燃料费用系数α_i,β_i,γ_i，α_i,β_i,γ_i分别表示第i台发电机组耗量特性的常数项、一次项系数和二次项系数，参数值由具体系统给出；阀点效应的常系数μ_i,η_i，具体值根据机组数不同而各异，从具体系统中获得；确定计算电网有功功率网络损耗P_loss的所需要的B系数B_ii,B_0i,B₀₀，B_ii,B_0i,B₀₀分别是根据机组的性质得到；h时段火电机组的中的额定功率

步骤2、产生初始的可行调度集合

可行调度集合的产生分两步实现：首先预测调度时段h的风力发电机组的有功功率；然后再产生火电机组的可行调度集合，具体过程为：

根据电力系统中风电机组的相关参数，第h时段内第l台发电机组的有功功率由下式来预测得到：

则h时段内所有风力发电机有功功率的输出预测值为

在调度时段h，针对前N-1个火力发电机组，在各个机组的发电容量要求下，随机产生有功功率i＝1,2,...,N-1，并利用系统功率平衡约束求出第N个发电机组的有功功率若产生的满足爬坡率约束式：max{Pimin,Pih-1-ζidown}≤Pih≤min{Pimax,Pih-1+ζiup},]]>则得到该时段的可行调度值为

否则重新进行此过程生成满足约束条件的调度值；

按照上述过程进行H次，就得到一个完整的可行调度序列；

将上述过程进行M次，就得到了初始的可行调度集合；

步骤3、对可行调度集合进行评价

依据目标函数式：minf(Pih)=Σh=1HΣi=1N{αi+βi×Pih+γi×(Pih)2+|μisin[ηi(Pimin-Pih)]|}]]>来评价调度的优劣，将步骤2得到的可行调度序列代入上述的目标函数式中，即获得相应调度序列适应值，记录所有M个可行调度序列中的最优值；

步骤4、对可行调度序列进行更新

4.1）惯性权重的生成

通过反馈控制结构生成惯性权重时，将每个可行调度序列作为反馈控制结构中的被控对象，每次迭代生成的可行调度序列的适应值作为反馈信息，实现过程如下：

对每一个可行调度序列，用它第iter次迭代结果的适应值作为输出量X_iter，通过反馈通道得到反馈量其中ψ是当前所有可行调度适应值的均值，采用e_iter＝1-Y_iter作为PID控制器输入来产生新的惯性权重：将ω_iter+1用来产生可行调度序列第iter+1次迭代结果，其中变量Y_iter作为对当前可行调度序列所获得反馈信息的一种度量，

若Y_iter＜1，则表示当前可行调度在所有调度中可行度高；

若Y_iter＞1，则表示当前可行调度在所有调度中可行度低；

4.2）对可行调度序列进行更新

在第iter次迭代中，m=1,2,L,...,M，第m个可行调度序列的第i个机组在第h时刻的有功功率为有功功率的变化率为第m个可行调度前iter次迭代历史中最优经济效益时第i个机组在第h时刻有功功率为前iter次迭代中所有可行调度中获得最优经济效益时第i个机组在第h时刻的有功功率为则在迭代次数为iter+1时，第m个可行调度中第i个机组在第h时刻的有功功率参照下述的两个迭代公式而生成：

Pmih(iter+1)=Pmih(iter)+Vmih(iter+1),]]>

Vmih(iter+1)=ωiter+1·Vmih(iter)+c1·r1·(Pmbestih(iter)-Pmih(iter))+c2r2(Pgbestih(iter)-Pmih(iter)),]]>

通过以上的两个迭代公式，就能够实现对每台发电机组不同时刻的有功功率和有功功率变化率的不断更新，寻找出新的M个更符合目标函数的发电机组发电功率的可行调度序列；

步骤5、对更新后的可行调度序列进行可行性判断

若新生成的可行调度序列，第h时刻的N台发电机组的有功功率i＝1,2,…,N满足系统功率平衡约束，且每台发电机组第h时刻与第h+1时刻有功功率间满足爬坡率约束，则新生成的调度序列可行，就继续迭代更新，否则返回步骤4重新生成可行调度序列；

步骤6、对新的可行调度集合进行评价

把更新后的M个可行调度序列分别代入上述目标函数式中进行优劣比较，确定出M个可行调度序列中的最优值，将这个最优值与历史最优值进行比较，取较小者为新的全局最优值；将每个更新后的可行调度序列当前的目标函数值与自身历史最优值进行比较，确定出新的局部最优值，分别记录下每个新的最优值以便下一次迭代使用；

步骤7、判断当前最好调度是否达到精度要求，或者当前迭代次数是否达到最大迭代次数iter_max，

若不是，设置迭代次数iter＝iter+1，继续重复步骤4到步骤6；

若是，输出当前的解集中的可行调度序列，即成。