CN105006844A - 一种间歇式发电并网条件下的电力系统日前鲁棒调度系统 - Google Patents

Abstract

本发明涉及电网设备技术领域，具体涉及一种间歇式发电并网条件下的电力系统日前鲁棒调度系统，包括向最优不确定集模块传输数据的采集数据模块，经过经济性最优对比和电网安全稳定运行约束条件设定后确定所述最优不确定集模块的参数，最优不确定集模块与建立的鲁棒机组组合系统通过混合整数规划求解器,如CPLEX计算出次日优化机组组合方案，次日优化机组组合方案经过安全校核后形成鲁棒调度日前计划后通过多时间尺度的电力系统鲁棒调度系统以文件形式或自动形式发布到能量管理系统。本发明适合大规模新能源并网时的机组组合系统，使生成的调度计划在保持机组状态不变的情况下应对间歇性能源的波动性，减小新能源并网后给电网带来的不利影响。

Description

一种间歇式发电并网条件下的电力系统日前鲁棒调度系统

技术领域

本发明涉及电网设备技术领域，具体涉及一种间歇式发电并网条件下的电力系统日前鲁棒调度系统。

背景技术

现有技术的传统调度应对风电波动的方法主要有两类，文献[1]WANG Yang，XIAQing.A novel security stochastic unit commitment for wind-thermal system operation[C]//2011 4th International Conference on Electric Utility Deregulation and Restructuringand Power Technologies(DRPT)(2011:386-393)和文献[2]《考虑大规模风电接入的备用容量计算[J]》(电力系统自动化，2012，36(6):1-5.王丹平，陈之栩，涂孟夫等)公开的是通过增加旋转备用来应对风电波动，这类方法简单易行，但备用容量既不易估计，又过高地增加了运行成本；另一类以随机规划技术为基础，文献[3]R.Barth，H.Brand，P.Meibom and C.Weber，“A stochastic unit-commitment model for theevaluation of the impacts of integration of large amounts of intermittent wind power，”inProc.Int.Conf.Probabilistic Methods Applied to Power Systems(PMAPS 2006)，(2006，11–15)和文献[4]E.Erdogan and G.Iyengar，“Ambiguous chance constrainedproblems and robust optimization，”(Math.Progr.2006，107(1):37–61)均提出了随机规划模型，采用随机变量描述不确定性，但这类方法局限性较多，一方面，必须提前得知精确的风功率概率分布，然而通过拟合得到的数据与实际存在一定差距，纵然获得精确数据，其形成的海量场景使计算效率降低，另一方面，若过度考虑风电功率的概率分布，有可能会破坏整个调度模型的凸规划性，除此之外，文献[5]《基于自适应粒子群优化算法的机组组合[J]》(电力系统保护与控制，2009，37(15):15-18.常文平，于海，华大鹏)和文献[6]《面向启发式调整策略和粒子群优化的机组组合问题》(电工技术学报，2009，24(12):132-141.袁晓辉，苏安俊，聂浩等)将粒子群算法应用到机组组合问题，是目前应用智能算法解决机组组合问题的典型方式。近年来鲁棒优化方法在电力系统调度中得到应用，文献[7]《一类电力系统鲁棒优化问题的博弈模型及应用实例[J]》(中国电机工程学报，2013，33(19):47-56.梅生伟，郭文涛，王莹莹等)建立了电力系统鲁棒优化的博弈模型。文献[8]《基于极限场景集的风电机组安全调度决策模型[J]》(电网技术，2013，37(6):1590-1595.高红均，刘俊勇，魏震波等)通过模拟风电概率曲线求得置信区间，在置信区间内进行鲁棒优化。文献[9]《电力系统鲁棒经济调度[J]》(电力系统自动化，2013，37(17):37-43.魏韡，刘锋，梅生伟)在此基础上考虑了风电场的时间和空间效益，更加贴合实际。在不确定集选取方面，通常采用固定不确定集或按置信区间确定不确定集的方法。前者太过依赖运行人员的经验判断，后者是一个显著的改进方案。但已有方法均没有进行不确定集的横向对比，无法说明采用的不确定集具有兼顾经济性和鲁棒性最优的综合效益。在多类型电源联合调度方面，文献[10]《基于收益最大化的水火电联合调度模型[J]》(电工技术学报，2010，25(7):150-158.吴杰康，曾建)提出了水火电联合调度时，水电在峰荷和基荷的分配原则，但依旧采用传统的确定性优化方法。文献[11]《储能系统对含风电的机组组合问题影响分析[J]》(电力系统自动化，2011，35(5):19-24.谢毓广，江晓东)将抽水蓄能机组引入风火联合调度模型，但根据蒙特卡洛法生成场景不适应风电波动较大的情况。文献[12]《多风电场并网时安全约束机组组合的混合整数规划解法[J]》(电力系统自动化，2010，34(5):29-33.叶荣，陈皓勇，王钢等)验证了极限场景对不确定集的代表性，实现了风火鲁棒调度模型，但考虑的能源种类过少。

发明内容

为了克服现有技术中的缺陷，本发明提供了一种间歇式发电并网条件下的电力系统日前鲁棒调度系统，能够应对间歇性能源大规模并网后给电网带来的冲击，提高电网消纳大规模间歇性能源的能力。通过日前计划离线计算，提前一天对次日发电机组启停及出力方案进行编排，时间分辨率为1小时。

本发明是通过如下技术方案实现的：一种间歇式发电并网条件下的电力系统日前鲁棒调度系统，所述系统包括采集数据模块，所述采集数据模块向最优不确定集模块传输数据，经过经济性最优对比和设定电网安全稳定运行约束条件后确定所述最优不确定集模块的参数，所述最优不确定集模块与建立的鲁棒机组组合系统通过混合整数规划求解器计算出次日优化机组组合方案，所述次日优化机组组合方案经过安全校核后形成鲁棒调度日前计划后通过多时间尺度的电力系统鲁棒调度系统自动形式发布到能量管理系统。

进一步地，所述采集数据模块包括间歇性电源预测数据模块、准备电网系统和机组特性数据模块，所述间歇性电源预测数据模块包括风电出力曲线预测模块及设置了参数的置信区间模块，所述准备电网系统和机组特性数据模块包括次日负荷预测曲线模块、系统安全稳定运行参数模块和各机组特性预测模块。

进一步地，所述最优不确定集模块通过综合成本运算模块采用优化方法确定最优不确定集，所述综合成本运算模块通过调度方案模块对综合成本进行统计后进行对比，所述综合成本运算模块包括运行成本模块、环境成本模块和风险成本模块，所述风险成本模块包括切负荷成本模块和弃风成本模块，所述弃风成本模块通过火电机会成本模块的计算来衡量弃风成本，所述火电机会成本模块包括火电发电成本模块和火电环境成本模块，火电机会成本为火电发电成本和火电环境成本之和。

进一步地，所述鲁棒机组组合系统包括预测场景模块和极限场景模块，所述极限场景模块包括功率平衡约束模块、线路安全约束模块、机组出力约束模块和调整时间约束模块，所述预测场景模块包括约束方程模块，所述约束方程模块包括机组侧约束模块和系统侧约束模块，所述机组侧约束模块包括火电机组约束模块、气电机组约束模块、水电机组约束模块、抽水蓄能机组约束模块和核电机组约束模块，所述系统侧约束模块包括功率平衡约束系统、旋转设备用约束系统和线路安全约束系统。

进一步地，所述环境成本模块包括发电环境成本S_e1和弃风/光环境成本S_e2，所述发电环境成本S_e1和弃风/光环境成本S_e2均包括排放物二氧化硫SO₂、氮氧化合物NO_x、粉尘TSP和二氧化碳CO₂，计算过程如下：

1)单位电量SO₂排放量计算公式为：

$Z_{\mathrm{SO}2}=Q\frac{64}{32}{S}_{\mathrm{ar}}{b}_{\mathrm{SO}2}(1-{\mathrm{\η}}_S)...\left(1\right);$

式中：Q为单位电量耗煤率，为SO₂与S分子量之比，S_ar为燃煤含硫质量分数，b_SO2为S生成SO₂的比例，η_S为脱硫效率；

2)单位电量NO_x排放量计算公式为：

$Z_{\mathrm{NOx}}=Q\frac{30}{14}{N}_{\mathrm{ar}}{b}_{\mathrm{NOx}}\left({1-\mathrm{\η}}_N\right)...\left(2\right);$

式中：为NO与N分子量之比，N_ar为燃煤含氮质量分数，b_NOx为N生成NO_x的比例，η_N为脱氮效率；

3)单位电量TSP排放量计算公式为：

Z_TSP＝Qb_TSP(1-η_TSP)……(3)；

式中：b_TSP为烟尘排放系数，η_TSP为脱尘效率；

4)单位电量CO₂排放成本计算公式为：

S＝QABP……(4),

式中：A为燃料低位发热量，B为CO₂排放系数，P为碳排放交易市场平均成交价；

所述发电环境成本模块包括计算公式：

$S_{e1}=M_{a1}(\underset{i}{\mathrm{\Σ}}0.6\frac{\mathrm{Zi}}{\mathrm{Ni}}+S)...\left(5\right),$

式中：M_a1为火电发电量，Z_i为单位电量的排放量，i为污染物种类，N_i为污染物当量值；

所述弃风/光环境成本计算公式为：

$S_{e2}=M_{a2}(\underset{i}{\mathrm{\Σ}}0.6\frac{\mathrm{Zi}}{\mathrm{Ni}}+S)...\left(6\right),$

式中：M_a2为弃新能源量，Z_i为单位电量的排放量，i为污染物种类，N_i为污染物当量值。

进一步地，所述风险成本模块包括风险成本S_r，所述风险成本S_r包括切负荷成本S_c和弃风/光成本S_a：

所述切负荷成本S_c计算公式为：

S_c＝M_cP_c……(7)，

总风险成本为弃风/光环境成本与切负荷成本之和，计算公式为：

S_f＝S_c+S_e2……(8)，

所述综合成本运算模块包括综合成本计算公式：

S_all＝S_gen+S_c+S_e+S_t……(9)，

式中，S_all为综合成本，S_gen为调度计划运行成本，S_c为切负荷成本，S_e和S_t分别为弃风成本中的火电环境成本和火电发电成本。

进一步地，所述风险成本模块包括统计风险成本，其中包括弃新能源量及切负荷量的计算公式为：

$\begin{array}{cc}{M}_{a2}=\underset{P_{\max }}{\overset{P_{\mathrm{cap}}}{\∫}}{P}_w{(P_t-P_{\max })}^d{P}_t& \∀(P_t>P_{\max })...\left(10\right),\end{array}$

$\begin{array}{cc}{M}_c=\underset{0}{\overset{P_{\min }}{\∫}}{P}_w(P_{\min }-P_t){\mathrm{dP}}_t& \∀(P_{\min }>P_t)...\left(11\right),\end{array}$

式中，M_a2、M_c分别为弃新能源量和切负荷量，P_w为新能源出力概率密度函数，P_min、P_max分别为调度计划上调极限和下调极限，P_t为风机出力，所有成本计算条件完整，进行以综合成本最优为目标的不确定集寻优计算。

进一步地，所述目标函数模块为以预测场景下系统的发电成本作为目标函数，计算公式为：

min.F(g,t)＝S_n+S_q+S_x+S_un+S_uq+S_en……(12)；

式中：S_n、S_q、S_x分别为火电、气电、核电的运行成本，S_un为火电机组启停成本，S_uq为气电机组启停成本，S_en为火电与气电的环境成本；

所述约束方程包括：

(1)所述系统侧约束模块包括：

①所述功率平衡约束系统的计算公式：

$\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{P}_g^t=D^t...\left(13\right);$

式中：P_g ^t为时段t各类电源机组出力，抽水蓄能机组抽水状态出力作负值计入，D^t为系统时段t的负荷预测值；

②所述旋转备用约束系统的计算公式：

$\underset{g\∈G_{\mathrm{aa}}}{\mathrm{\Σ}}{I}_g^t{P}_{g\max }+\underset{{g\∈G}_{\mathrm{ab}}}{\mathrm{\Σ}}{I}_g^t{P}_g^t\≥D^t+R^t...\left(14\right);$

式中：G_aa为除了风机和虚拟电动机之外的所有机组，G_ab为风机和虚拟电动机机组，R^t为系统时段t的旋转备用需求。

③所述线路安全约束系统的计算公式：

$\left|\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{gj}}{P}_g^t\right|\≤L_j...\left(15\right);$

式中：G_all为所有机组，γ_gj为机组g在线路j上的功率分布因子，L_j为线路流量限制。

(2)所述机组侧约束模块包括：

①所述火电机组约束模块包括：

机组运行成本计算公式为：

$S_n=\underset{t\∈T}{\mathrm{\Σ}}\underset{g{\∈G}_n}{\mathrm{\Σ}}{A}_n\×{\left(P_n^t\right)}^2+B_n\×P_n^t+C_n...\left(16\right),$

式中：T为调度的时段数，A_n、B_n、C_n为火电机组n成本系数，P_n ^t为火电机组n在时段t的出力；

机组提供旋转备用容量计算公式为：

$\underset{{n\∈G}_n}{\mathrm{\Σ}}{R}_n^t=\underset{{n\∈G}_n}{\mathrm{\Σ}}{I}_n^t(P_{n\max }-P_n^t)...\left(17\right),$

式中：I_n ^t为火电机组n在时段t的状态，P^h _nmax为火电机组n在时段t的最大出力；

机组出力约束计算公式为：

$P_{n\min }\≤P_n^t\≤P_{n\max }...\left(18\right),$

式中：P_nmix、P_nmax分别为火电机组n的最小出力与最大出力；

机组爬坡率约束：

${-r}_d^n\≤P_n^t-P_n^{t-1}\≤r_u^n...\left(19\right),$

式中：r_u ⁿ为火电机组n的向上爬坡率，r_d ⁿ为火电机组n的向下爬坡率；

机组开关机约束计算公式为：

$\underset{h=0}{\overset{T_u\left(n\right)-1}{\mathrm{\Σ}}}I(n,t+h)\≥(I(n,t)-I(n,t-1))\×\min (T_u\left(n\right),N_T-t+1)...\left(20\right),$

$\underset{h=0}{\overset{T_d\left(n\right)-1}{\mathrm{\Σ}}}(1-I(n,t+h))\≥(I(n,t-1)-I(n,t))\×\min (T_d\left(n\right),N_T-t+1)...\left(21\right);$

式中：T_d(n)、T_u(n)为火电机组n的最小关机、开机时间；

②所述气电机组约束模块包括：

机组运行成本计算公式为：

$S_q=\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈G_q}{\mathrm{\Σ}}{A}_q\×{\left(P_q^t\right)}^2+B_q\×P_q^t+C_q...\left(22\right),$

式中：T为调度的时段数，A_q、B_q、C_q为气电机组q成本系数，P_q ^t为气电机组q在时段t的出力；

机组提供旋转备用容量计算公式为：

$\underset{{q\∈G}_q}{\mathrm{\Σ}}{R}_q^t=\underset{{q\∈G}_q}{\mathrm{\Σ}}{I}_q^t(P_{q\max }-P_q^t)...\left(23\right),$

式中：I_q ^t为气电机组q在时段t的状态，P^t _qmax为气电机组q在时段t的最大出力；

机组出力约束计算公式为：

$P_{q\min }\≤P_q^t\≤P_{q\max }...\left(24\right),$

式中：P_qmix、P_qmax分别为气电机组q的最小出力与最大出力；

机组爬坡率约束：

${-r}_d^q\≤P_q^t-P_q^{t-1}\≤r_u^q...\left(25\right),$

式中：r_u ^q为气电机组q的向上爬坡率，r_d ^q为气电机组q的向下爬坡率；

机组开关机约束计算公式为：

$\underset{h=0}{\overset{T_u\left(q\right)-1}{\mathrm{\Σ}}}I(q,t+h)\≥(I(q,t)-I(q,t-1))\×\min (T_u\left(q\right),N_T-t+1)...\left(26\right),$

$\underset{h=0}{\overset{T_d\left(q\right)-1}{\mathrm{\Σ}}}(1-I(q,t+h))\≥(I(q,t-1)-I(q,t))\×\min (T_d\left(q\right),N_T-t+1)...\left(27\right),$

气量约束计算公式为：

$\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈\mathrm{Gq}}{\mathrm{\Σ}}{P}_q^t\≤Q_q...\left(28\right);$

式中：Q_q为因气量限制气电厂一天的最大发电量；

③所述水电机组约束模块包括：

机组提供旋转备用容量计算公式为：

$\underset{s\∈G_s}{\mathrm{\Σ}}{R}_s^t=\underset{s\∈G_s}{\mathrm{\Σ}}{I}_s^t(P_{s\max }-P_s^t)...\left(29\right),$

式中：为水电机组s在时段t的状态，为水电机组s在时段t的最大出力；

机组出力约束计算公式为：

$P_{s\min }\≤P_s^t\≤P_{s\max }...\left(30\right),$

式中：P_smin、P_smax分别为水电机组s的最小出力与最大出力；

水量约束计算公式为：

$\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈G_s}{\mathrm{\Σ}}{P}_s^t\≤Q_s...\left(31\right);$

式中：Q_s为因库容限制水电厂一天的最大发电量；

④所述核电机组约束模块包括：

机组运行成本计算公式：

$S_x=\underset{t\∈T}{\mathrm{\Σ}}\underset{g{\∈G}_x}{\mathrm{\Σ}}{B}_x\×P_x^t+C_x...\left(32\right),$

式中：T为调度的时段数，B_x、C_x为核电机组x成本系数，P_x ^t为核电机组x在时段t的出力；

机组提供旋转备用容量计算公式：

$\underset{{x\∈G}_x}{\mathrm{\Σ}}{R}_x^t=\underset{{x\∈G}_x}{\mathrm{\Σ}}{I}_x^t(P_{x\max }-P_x^t),...\left(33\right),$

式中：I_x ^t为核电机组x在时段t的状态，P^t _xmax为核电机组x在时段t的最大出力；

机组出力约束计算公式为：

$P_{x\min }\≤P_x^t\≤P_{x\max }...\left(34\right),$

式中：P_xmin、P_xmax分别为核电机组x的最小出力与最大出力；

机组爬坡率约束计算公式为：

${-r}_d^x\≤P_x^t-P_x^{t-1}\≤r_u^x...\left(35\right),$

式中：r_u ^x为核电机组x的向上爬坡率，r_d ^x为核电机组x的向下爬坡率；

机组开关机约束计算公式为：

I(x,t)＝1……(36)；

⑤所述抽水蓄能机组约束模块包括：

机组提供旋转备用容量计算公式：

$\begin{array}{cc}\underset{\mathrm{cg}{\∈G}_{\mathrm{cg}}}{\mathrm{\Σ}}{P}_{\mathrm{cg}\max }-\underset{\mathrm{cg}\∈G_{\mathrm{cg}}}{\mathrm{\Σ}}{P}_{\mathrm{cg}}^t\≤Y_0+\underset{\mathrm{cm}\∈G_{\mathrm{cm}}}{\mathrm{\Σ}}\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cm}}^t-\underset{\mathrm{cg}\∈G_{\mathrm{cg}}}{\mathrm{\Σ}}\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cg}}^t& (t\∈T)...\left(37\right),\end{array}$

$\underset{\mathrm{cg}{\∈G}_{\mathrm{cg}}}{\mathrm{\Σ}}{R}_c^t=\underset{\mathrm{cg}{\∈G}_{\mathrm{cg}}}{\mathrm{\Σ}}{I}_{\mathrm{cg}}^t(P_{\mathrm{cg}\max }-P_{\mathrm{cg}}^t)...\left(38\right),$

式中：用虚拟发电机cg和虚拟电动机cm代替抽水蓄能机组c，Y₀为水库初始所存电量，为虚拟发电机cg在时段t的状态，P_cgmax为虚拟发电机cg的最大出力；

机组出力约束计算公式：

$P_{\mathrm{cg}\min }\≤P_{\mathrm{cg}}^t\≤P_{\mathrm{cg}\max }...\left(39\right),$

$P_{\mathrm{cm}}^t=P_{\mathrm{cmr}}...\left(40\right),$

式中：P_cgmin、P_cgmax分别为虚拟发电机cg的最小出力与最大出力，P_cmr为虚拟电动机cm的额定出力；

机组状态约束计算公式：

$I_{\mathrm{cm}}^t+I_{\mathrm{cg}}^t\≤1...\left(41\right),$

式中：分别为虚拟电动机、虚拟发电机在时段t的状态；

库容约束计算公式：

$\begin{array}{cc}{W}_l\≤W_0+e\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cm}}^t-\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cg}}^t\≤W_u& (t\∈T)...\left(42\right),\end{array}$

式中：W_l、W_u分别为水库容量的下限、上限，W₀为水库初始容量，e为抽水蓄能机组的效率；

库容设定约束计算公式为：

$e\underset{t\∈T}{\mathrm{\Σ}}{P}_{\mathrm{cm}}^t-\underset{t\∈T}{\mathrm{\Σ}}{P}_{\mathrm{cg}}^t=W_1-W_0...\left(43\right),$

式中：W₁为调度结束水库容量的设定值，默认等于W₀；

所述极限场景约束系统包括：

所述功率平衡约束模块包括计算公式，

$\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{P}_{g,s}^t=D^t...\left(44\right),$

式中：P^t _g,s为极限场景s下时段t各类电源机组出力，抽水蓄能机组抽水状态出力作负值计入，D^t为系统时段t的负荷预测值；

所述线路安全约束模块包括计算公式：

$\left|\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{gj}}{P}_{g,s}^t\right|\≤L_j...\left(45\right),$

式中：γ_gj为机组g在线路j上的功率分布因子，L_j为线路流量限制；

所述机组出力约束模块包括计算公式：

$\begin{array}{cc}{P}_{g\min }\≤P_{g,s}^t\≤P_{g\max }& (g\∈G_{\mathrm{all}})...\left(46\right),\end{array}$

所述调整时间约束模块包括计算公式：

$\begin{array}{cc}-\frac{1}{6}{r}_d^g\≤P_{g,s}^t-P_g^t\≤\frac{1}{6}{r}_u^g& (g\∈G_{\mathrm{all}})...\left(47\right).\end{array}$

式中：为极限场景s下时段t各类电源机组出力，这里选择小时爬坡率的1/6为约束上下界，即要求机组在10min内完成场景之间的过渡，也可根据实际运行需要选取其他值。

与现有技术相比，优越效果在于：本发明提供一种间歇式发电并网条件下的电力系统日前鲁棒调度系统，能够适合大规模新能源，尤其风电并网时的机组组合，使生成的调度计划能在保持机组状态不变的情况下应对间歇性能源的波动性，减小新能源并网后给电网带来的不利影响，保证电网安全、经济运行。

附图说明

图1为本发明所述间歇式发电并网条件下的电力系统日前鲁棒调度系统模块示意图；

图2为本发明所述间歇式发电并网条件下的电力系统日前鲁棒调度系统中的采集数据模块示意图；

图3为本发明所述间歇式发电并网条件下的电力系统日前鲁棒调度系统中的综合成本运算模块的示意图；

图4为本发明所述间歇式发电并网条件下的电力系统日前鲁棒调度系统中的鲁棒机组组合系统示意图；

图5为本发明所述间歇式发电并网条件下的电力系统日前鲁棒调度系统中的间歇式发电出力概率分布曲线示意图；

图6是本发明所述间歇式发电并网条件下的电力系统日前鲁棒调度系统中的日前计划编制功能流程图；

图7是本发明所述间歇式发电并网条件下的电力系统日前鲁棒调度系统中的负荷及新能源的短期功率预测数据示意图；

图8是本发明所述间歇式发电并网条件下的电力系统日前鲁棒调度系统中的不同不确定集下的综合成本示意图。

具体实施方式

下面结合附图对本发明具体实施方式作进一步详细说明。

如图1-4所示，具体说明本发明，本发明提供的一种间歇式发电并网条件下的电力系统日前鲁棒调度系统，包括采集数据模块，所述采集数据模块向最优不确定集模块传输数据，经过经济性最优对比和设定电网安全稳定运行约束条件后确定所述最优不确定集模块的参数，所述最优不确定集模块与建立的鲁棒机组组合系统通过混合整数规划求解器计算出次日的优化机组组合方案，所述次日优化机组组合方案经过安全校核后形成鲁棒调度日前计划后通过多时间尺度的电力系统鲁棒调度系统以文件形式或自动形式发布到能量管理系统，所述采集数据模块包括间歇性电源预测数据模块、准备电网系统和机组特性数据模块，所述间歇性电源预测数据模块包括风电出力曲线预测模块及设置了相关参数的置信区间模块，所述准备电网系统和机组特性数据模块包括次日负荷预测曲线模块、系统安全稳定运行参数模块和各机组特性预测模块，所述最优不确定集模块通过综合成本运算模块采用优化方法确定最优不确定集，所述综合成本运算模块通过调度方案模块对综合成本进行统计后进行对比，所述综合成本运算模块包括运行成本模块、环境成本模块和风险成本模块，所述风险成本模块包括切负荷成本模块和弃风成本模块，所述弃风成本模块通过火电机会成本模块的计算来衡量弃风成本，所述火电机会成本模块包括火电发电成本模块和火电环境成本模块，火电机会成本为火电发电成本和火电环境成本之和。所述鲁棒机组组合系统包括预测场景模块和极限场景模块，所述极限场景模块包括功率平衡约束模块、线路安全约束模块、机组出力约束模块和调整时间约束模块，所述预测场景模块包括约束方程模块，所述约束方程模块包括机组侧约束模块和系统侧约束模块，所述机组侧约束模块包括火电机组约束模块、气电机组约束模块、水电机组约束模块、抽水蓄能机组约束模块和核电机组约束模块，所述系统侧约束模块包括功率平衡约束系统、旋转设备用约束系统和线路安全约束系统。

不确定集是包含不确定因素各种可能场景的集合。调度计划费用称为综合成本，其包括运行成本和事故成本，事故成本包括切负荷和弃风两种措施产生的成本，弃风成本用火电机会成本衡量，可分为火电发电成本和环境成本。通过调度方案统计运行成本和事故成本，求和得到综合成本，能够使综合成本最低的不确定集即为最优不确定集。以综合成本最优为目标来寻找最优不确定集，首先是几类成本的定义：

运行成本：运行成本是机组在调度周期内发电所消耗燃料的成本。

环境成本：环境成本是指为处理发电产生的燃料废渣废气而付出的成本，包括发电环境成本和弃风/光环境成本两部分。

发电环境成本：发电环境成本是调度计划消耗燃料产生的环境成本。

弃风/光环境成本：弃风/光环境成本是指弃风/光的环境效益，即弃风/光量等效火电发电产生的环境成本。

风险成本：风险成本是当新能源波动超过调度计划消纳能力而采用的弃风/光或切负荷措施带来的成本。

综合成本：综合成本是指将运行成本、环境成本和风险成本之和。

不确定集是鲁棒调度结果所必须满足的新能源出力波动集合，其范围越大，调度计划消耗燃料越多，应对新能源波动的能力越强，风险成本越低，运行成本和环境成本越高；相反，其范围越小，调度计划消耗燃料越少，应对新能源波动的能力越弱，风险成本越高，运行成本越低因此如何确定新能源最优不确定集是问题关键所在。

1)单位电量SO₂排放量计算公式为：

$Z_{\mathrm{SO}2}=Q\frac{64}{32}{S}_{\mathrm{ar}}{b}_{\mathrm{SO}2}(1-{\mathrm{\η}}_S)...\left(1\right);$

式中：Q为单位电量耗煤率，为SO₂与S分子量之比，S_ar为燃煤含硫质量分数，b_SO2为S生成SO₂的比例，η_S为脱硫效率；

2)单位电量NO_x排放量计算公式为：

$Z_{\mathrm{NOx}}=Q\frac{30}{14}{N}_{\mathrm{ar}}{b}_{\mathrm{NOx}}\left({1-\mathrm{\η}}_N\right)...\left(2\right);$

式中：为NO与N分子量之比，N_ar为燃煤含氮质量分数，b_NOx为N生成NO_x的比例，η_N为脱氮效率；

3)单位电量TSP排放量计算公式为：

Z_TSP＝Qb_TSP(1-η_TSP)……(3)；

式中：b_TSP为烟尘排放系数，η_TSP为脱尘效率；

表1污染物当量值

4)单位电量CO₂排放成本计算公式为：

S＝QABP……(4)，

式中：A为燃料低位发热量，B为CO₂排放系数，P为碳排放交易市场平均成交价；

所述发电环境成本模块包括计算公式：

$S_{e1}=M_{a1}(\underset{i}{\mathrm{\Σ}}0.6\frac{\mathrm{Zi}}{\mathrm{Ni}}+S)...\left(5\right),$

式中：M_a1为火电发电量，Z_i为单位电量的排放量，i为污染物种类，N_i为污染物当量值；

所述弃风/光环境成本计算公式为：

$S_{e2}=M_{a2}(\underset{i}{\mathrm{\Σ}}0.6\frac{\mathrm{Zi}}{\mathrm{Ni}}+S)...\left(6\right),$

式中：M_a2为弃新能源量，Z_i为单位电量的排放量，i为污染物种类，N_i为污染物当量值。

进一步地，所述风险成本模块包括风险成本S_r，所述风险成本S_r包括切负荷成本S_c和弃风/光成本S_a：

所述切负荷成本S_c计算公式为：

S_c＝M_cP_c……(7)，

总风险成本为弃风/光环境成本与切负荷成本之和：

S_f＝S_c+S_e2……(8),

所述综合成本运算模块包括综合成本计算公式：

S_all＝S_gen+S_c+S_e+S_t……(9),

式中，S_all为综合成本，S_gen为调度计划运行成本，S_c为切负荷成本，S_e和S_t分别为弃风成本中的火电环境成本和火电发电成本。

进一步地，所述风险成本模块包括统计风险成本，其中包括弃新能源量及切负荷量的计算公式为：

$\begin{array}{cc}{M}_{a2}=\underset{P_{\max }}{\overset{P_{\mathrm{cap}}}{\∫}}{P}_w{(P_t-P_{\max })}^d{P}_t& \∀(P_t>P_{\max })...\left(10\right),\end{array}$

$\begin{array}{cc}{M}_c=\underset{0}{\overset{P_{\min }}{\∫}}{P}_w(P_{\min }-P_t){\mathrm{dP}}_t& \∀(P_{\min }>P_t)...\left(11\right),\end{array}$

式中，M_a2、M_c分别为弃新能源量和切负荷量，P_w为新能源出力概率密度函数，P_min、P_max分别为调度计划上调极限和下调极限，P_t为风机出力，所有成本计算条件完整，进行以综合成本最优为目标的不确定集寻优计算。

进一步地，所述目标函数模块为以预测场景下系统的发电成本作为目标函数，计算公式为：

min.F(g,t)＝S_n+S_q+S_x+S_un+S_uq+S_en……(12)；

式中：S_n、S_q、S_x分别为火电、气电、核电的运行成本，S_un为火电机组启停成本，S_uq为气电机组启停成本，S_en为火电与气电的环境成本；

所述约束方程包括：

(1)所述系统侧约束模块包括：

①所述功率平衡约束系统的计算公式：

$\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{P}_g^t=D^t...\left(13\right);$

式中：P_g ^t为时段t各类电源机组出力，抽水蓄能机组抽水状态出力作负值计入，D^t为系统时段t的负荷预测值；

②所述旋转备用约束系统的计算公式：

$\underset{g\∈G_{\mathrm{aa}}}{\mathrm{\Σ}}{I}_g^t{P}_{g\max }+\underset{{g\∈G}_{\mathrm{ab}}}{\mathrm{\Σ}}{I}_g^t{P}_g^t\≥D^t+R^t...\left(14\right);$

式中：G_aa为除了风机和虚拟电动机之外的所有机组，G_ab为风机和虚拟电动机机组，R^t为系统时段t的旋转备用需求。

③所述线路安全约束系统的计算公式：

$\left|\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{gj}}{P}_g^t\right|\≤L_j...\left(15\right);$

式中：G_all为所有机组，γ_gj为机组g在线路j上的功率分布因子，L_j为线路流量限制。

(2)所述机组侧约束模块包括：

①所述火电机组约束模块包括：

机组运行成本计算公式为：

$S_n=\underset{t\∈T}{\mathrm{\Σ}}\underset{g{\∈G}_n}{\mathrm{\Σ}}{A}_n\×{\left(P_n^t\right)}^2+B_n\×P_n^t+C_n...\left(16\right),$

式中：T为调度的时段数，A_n、B_n、C_n为火电机组n成本系数，P_n ^t为火电机组n在时段t的出力；

机组提供旋转备用容量计算公式为：

$\underset{{n\∈G}_n}{\mathrm{\Σ}}{R}_n^t=\underset{{n\∈G}_n}{\mathrm{\Σ}}{I}_n^t(P_{n\max }-P_n^t)...\left(17\right),$

式中：I_n ^t为火电机组n在时段t的状态，P^h _nmax为火电机组n在时段t的最大出力；

机组出力约束计算公式为：

$P_{n\min }\≤P_n^t\≤P_{n\max }...\left(18\right),$

式中：P_nmix、P_nmax分别为火电机组n的最小出力与最大出力；

机组爬坡率约束：

${-r}_d^n\≤P_n^t-P_n^{t-1}\≤r_u^n...\left(19\right),$

式中：r_u ⁿ为火电机组n的向上爬坡率，r_d ⁿ为火电机组n的向下爬坡率；

机组开关机约束计算公式为：

$\underset{h=0}{\overset{T_u\left(n\right)-1}{\mathrm{\Σ}}}I(n,t+h)\≥(I(n,t)-I(n,t-1))\×\min (T_u\left(n\right),N_T-t+1)...\left(20\right),$

$\underset{h=0}{\overset{T_d\left(n\right)-1}{\mathrm{\Σ}}}(1-I(n,t+h))\≥(I(n,t-1)-I(n,t))\×\min (T_d\left(n\right),N_T-t+1)...21),$

式中：T_d(n)、T_u(n)为火电机组n的最小关机、开机时间；

②所述气电机组约束模块包括：

机组运行成本计算公式为：

$S_q=\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈G_q}{\mathrm{\Σ}}{A}_q\×{\left(P_q^t\right)}^2+B_q\×P_q^t+C_q...\left(22\right),$

式中：T为调度的时段数，A_q、B_q、C_q为气电机组q成本系数，P_q ^t为气电机组q在时段t的出力；

机组提供旋转备用容量计算公式为：

$\underset{{q\∈G}_q}{\mathrm{\Σ}}{R}_q^t=\underset{{q\∈G}_q}{\mathrm{\Σ}}{I}_q^t(P_{q\max }-P_q^t)...\left(23\right),$

式中：I_q ^t为气电机组q在时段t的状态，P^t _qmax为气电机组q在时段t的最大出力；

机组出力约束计算公式为：

$P_{q\min }\≤P_q^t\≤P_{q\max }...\left(24\right),$

式中：P_qmix、P_qmax分别为气电机组q的最小出力与最大出力；

机组爬坡率约束：

${-r}_d^q\≤P_q^t-P_q^{t-1}\≤r_u^q...\left(25\right),$

式中：r_u ^q为气电机组q的向上爬坡率，r_d ^q为气电机组q的向下爬坡率；

机组开关机约束计算公式为：

$\underset{h=0}{\overset{T_u\left(q\right)-1}{\mathrm{\Σ}}}I(q,t+h)\≥(I(q,t)-I(q,t-1))\×\min (T_u\left(q\right),N_T-t+1)...\left(26\right),$

$\underset{h=0}{\overset{T_d\left(q\right)-1}{\mathrm{\Σ}}}(1-I(q,t+h))\≥(I(q,t-1)-I(q,t))\×\min (T_d\left(q\right),N_T-t+1)...\left(27\right),$

气量约束计算公式为：

$\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈\mathrm{Gq}}{\mathrm{\Σ}}{P}_q^t\≤Q_q...\left(28\right);$

式中：Q_q为因气量限制气电厂一天的最大发电量；

③所述水电机组约束模块包括：

机组提供旋转备用容量计算公式为：

$\underset{s\∈G_s}{\mathrm{\Σ}}{R}_s^t=\underset{s\∈G_s}{\mathrm{\Σ}}{I}_s^t(P_{s\max }-P_s^t)...\left(29\right),$

式中：为水电机组s在时段t的状态，为水电机组s在时段t的最大出力；

机组出力约束计算公式为：

$P_{s\min }\≤P_s^t\≤P_{s\max }...\left(30\right),$

式中：P_smin、P_smax分别为水电机组s的最小出力与最大出力；

水量约束计算公式为：

$\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈G_s}{\mathrm{\Σ}}{P}_s^t\≤Q_s...\left(31\right);$

式中：Q_s为因库容限制水电厂一天的最大发电量；

④所述核电机组约束模块包括：

机组运行成本计算公式：

$S_x=\underset{t\∈T}{\mathrm{\Σ}}\underset{g{\∈G}_x}{\mathrm{\Σ}}{B}_x\×P_x^t+C_x...\left(32\right),$

式中：T为调度的时段数，B_x、C_x为核电机组x成本系数，P_x ^t为核电机组x在时段t的出力；

机组提供旋转备用容量计算公式：

$\underset{{x\∈G}_x}{\mathrm{\Σ}}{R}_x^t=\underset{{x\∈G}_x}{\mathrm{\Σ}}{I}_x^t(P_{x\max }-P_x^t)...\left(33\right),$

式中：I_x ^t为核电机组x在时段t的状态，P^t _xmax为核电机组x在时段t的最大出力；

机组出力约束计算公式为：

$P_{x\min }\≤P_x^t\≤P_{x\max }...\left(34\right),$

式中：P_xmin、P_xmax分别为核电机组x的最小出力与最大出力；

机组爬坡率约束计算公式为：

${-r}_d^x\≤P_x^t-P_x^{t-1}\≤r_u^x...\left(35\right),$

式中：r_u ^x为核电机组x的向上爬坡率，r_d ^x为核电机组x的向下爬坡率；

机组开关机约束计算公式为：

I(x,t)＝1……(36)；

⑤所述抽水蓄能机组约束模块包括：

机组提供旋转备用容量计算公式：

$\begin{array}{cc}\underset{\mathrm{cg}{\∈G}_{\mathrm{cg}}}{\mathrm{\Σ}}{P}_{\mathrm{cg}\max }-\underset{\mathrm{cg}\∈G_{\mathrm{cg}}}{\mathrm{\Σ}}{P}_{\mathrm{cg}}^t\≤Y_0+\underset{\mathrm{cm}\∈G_{\mathrm{cm}}}{\mathrm{\Σ}}\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cm}}^t-\underset{\mathrm{cg}\∈G_{\mathrm{cg}}}{\mathrm{\Σ}}\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cg}}^t& (t\∈T)...\left(37\right),\end{array}$

$\underset{\mathrm{cg}{\∈G}_{\mathrm{cg}}}{\mathrm{\Σ}}{R}_c^t=\underset{\mathrm{cg}{\∈G}_{\mathrm{cg}}}{\mathrm{\Σ}}{I}_{\mathrm{cg}}^t(P_{\mathrm{cg}\max }-P_{\mathrm{cg}}^t)...\left(38\right),$

式中：用虚拟发电机cg和虚拟电动机cm代替抽水蓄能机组c，Y₀为水库初始所存电量，为虚拟发电机cg在时段t的状态，P_cgmax为虚拟发电机cg的最大出力；

机组出力约束计算公式：

$P_{\mathrm{cg}\min }\≤P_{\mathrm{cg}}^t\≤P_{\mathrm{cg}\max }...\left(39\right),$

$P_{\mathrm{cm}}^t=P_{\mathrm{cmr}}...\left(40\right),$

式中：P_cgmin、P_cgmax分别为虚拟发电机cg的最小出力与最大出力，P_cmr为虚拟电动机cm的额定出力；

机组状态约束计算公式：

$I_{\mathrm{cm}}^t+I_{\mathrm{cg}}^t\≤1...\left(41\right),$

式中：分别为虚拟电动机、虚拟发电机在时段t的状态；

库容约束计算公式：

$\begin{array}{cc}{W}_l\≤W_0+e\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cm}}^t-\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cg}}^t\≤W_u& (t\∈T)...\left(42\right),\end{array}$

式中：W_l、W_u分别为水库容量的下限、上限，W₀为水库初始容量，e为抽水蓄能机组的效率；

库容设定约束计算公式为：

$e\underset{t\∈T}{\mathrm{\Σ}}{P}_{\mathrm{cm}}^t-\underset{t\∈T}{\mathrm{\Σ}}{P}_{\mathrm{cg}}^t=W_1-W_0...\left(43\right),$

式中：W₁为调度结束水库容量的设定值，默认等于W₀。

所述极限场景约束系统包括：

所述功率平衡约束模块包括计算公式，

$\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{P}_{g,s}^t=D^t...\left(44\right),$

式中：P^t _g,s为极限场景s下时段t各类电源机组出力，抽水蓄能机组抽水状态出力作负值计入，D^t为系统时段t的负荷预测值；

所述线路安全约束模块包括计算公式：

$\left|\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{gj}}{P}_{g,s}^t\right|\≤L_j...\left(45\right),$

式中：γ_gj为机组g在线路j上的功率分布因子，L_j为线路流量限制；

所述机组出力约束模块包括计算公式：

$\begin{array}{cc}{P}_{g\min }\≤P_{g,s}^t\≤P_{g\max }& (g\∈G_{\mathrm{all}})...\left(46\right),\end{array}$

所述调整时间约束模块包括计算公式：

$\begin{array}{cc}-\frac{1}{6}{r}_d^g\≤P_{g,s}^t-P_g^t\≤\frac{1}{6}{r}_u^g& (g\∈G_{\mathrm{all}})...\left(47\right),\end{array}$

式中：为极限场景s下时段t各类电源机组出力，这里选择小时爬坡率的1/6为约束上下界，即要求机组在10min内完成场景之间的过渡，也可根据实际运行需要选取其他值。

本实施例中，采用混合整数规划CPLEX求解器求解计算；其步骤具体为：首先在预测场景下得到最优解，然后校验其能否满足极限场景约束和安全约束，若能则输出结果，若不能则返回输入端调整发电计划重新计算。本系统中，每日计算得到次日优化机组组合方案，由多时间尺度的电力系统鲁棒调度系统将经过安全校核的鲁棒调度日前计划以文件方式或者自动发布的方式发布到能量管理系统。

本发明的工作原理：本发明综合考虑了水、火、风、气、核、抽水蓄能的多种能源结构，提出一种含多类型能源的鲁棒机组组合系统，根据风电功率分布及一定置信水平来找到风功率预测误差带，通过引入极限场景来应对大规模风电的预测误差给电网安全运行带来的影响。同时，系统考虑环境保护成本，适合当前构建生态文明的发展战略。

如图6所示，首先，采集间歇性电源的预测数据，准备电网系统及机组特性数据；其次，以经济性最优为目标函数，以保证电网安全稳定运行为约束条件，运用极限场景法应对间歇性能源波动，建立鲁棒机组组合系统；再次，通过调用多时间尺度的电力系统鲁棒调度系统的静态安全校核服务，实现鲁棒调度日前计划的静态安全校核，当出现支路或断面潮流越限时，采用安全校正算法调整发电计划以消除支路及断面的潮流约束；最后，每日计算得到次日优化机组组合方案，由多时间尺度的电力系统鲁棒调度系统将经过安全校核的鲁棒调度日前计划以文件方式或者自动发布的方式发布到能量管理系统。

以下是本发明系统的一个实际算例，以2014年国内某省级电网为例进行仿真计算。其中，火电占47.1％，气电占11.4％，水电占7.7％，核电占15.4％，抽水蓄能占3％，风电占13.2％，光伏发电占2.5％。

如图7所示，其中新能源的容量渗透率为24.5％，近似为逆调峰特性。

如图8所示，随着不确定集的增大，“运行+环境成本”逐渐升高，而风险成本逐渐降低，其两者之和综合成本在20％-30％范围内达到最低。在20％-30％范围内以1％为误差范围，采用二分搜索法寻找综合成本最低值，经搜索，确定最优不确定集比例为25％，因此按25％缩放比例形成的不确定集即为最优不确定集。

表2是传统法与鲁棒法的经济性比较。鲁棒法运行成本高，但能有效避免弃风及切负荷带来的风险损失，综合成本较传统法低。

定义综合效益为鲁棒法相对传统法的综合成本减少量。表3统计了不同运行情况下的综合效益。由此可知，1)新能源为逆调峰时综合效益更高；2)新能源渗透率越高，综合效益越高。

表2传统法与鲁棒法的经济性比较(元)

系统	风险成本	运行成本	综合成本
				传统法	78668	12002210	12080878
鲁棒法	44374	12019087	12063461

表3不同情况下综合效益比较(元)

新能源渗透率	新能源特性	综合效益
			12.7％	逆调峰	153330
12.7％	顺调峰	17416
			19.1％	顺调峰	490245

本发明并不限于上述实施方式，在不背离本发明的实质内容的情况下，本领域技术人员可以想到的任何变形、改进、替换均落入本发明的保护范围。

Claims (5)

1.一种间歇式发电并网条件下的电力系统日前鲁棒调度系统，其特征在于，所述系统包括采集数据模块，所述采集数据模块向最优不确定集模块传输数据，经过经济性最优对比和设定电网安全稳定运行约束条件后确定所述最优不确定集模块的参数，所述最优不确定集模块与建立的鲁棒机组组合系统通过混合整数规划求解器计算出次日优化机组组合方案，所述次日优化机组组合方案经过安全校核后形成鲁棒调度日前计划后通过多时间尺度的电力系统鲁棒调度系统自动形式发布到能量管理系统。

2.根据权利要求1所述间歇式发电并网条件下的电力系统日前鲁棒调度系统，其特征在于，所述采集数据模块包括间歇性电源预测数据模块、准备电网系统和机组特性数据模块，所述间歇性电源预测数据模块包括风电出力曲线预测模块及设置了相关参数的置信区间模块，所述准备电网系统和机组特性数据模块包括次日负荷预测曲线模块、系统安全稳定运行参数模块和各机组特性预测模块。

3.根据权利要求1所述间歇式发电并网条件下的电力系统日前鲁棒调度系统，其特征在于，所述最优不确定集模块通过综合成本运算模块采用优化方法确定最优不确定集，所述综合成本运算模块通过调度方案模块对综合成本进行统计后进行对比，所述综合成本运算模块包括运行成本模块、环境成本模块和风险成本模块，所述风险成本模块包括切负荷成本模块和弃风成本模块，所述弃风成本模块通过火电机会成本模块的计算来衡量弃风成本，所述火电机会成本模块包括火电发电成本模块和火电环境成本模块，火电机会成本为火电发电成本和火电环境成本之和。

4.根据权利要求1所述间歇式发电并网条件下的电力系统日前鲁棒调度系统，其特征在于，所述鲁棒机组组合系统包括预测场景模块和极限场景模块，所述极限场景模块包括功率平衡约束模块、线路安全约束模块、机组出力约束模块和调整时间约束模块，所述预测场景模块包括约束方程模块，所述约束方程模块包括机组侧约束模块和系统侧约束模块，所述机组侧约束模块包括火电机组约束模块、气电机组约束模块、水电机组约束模块、抽水蓄能机组约束模块和核电机组约束模块，所述系统侧约束模块包括功率平衡约束系统、旋转设备用约束系统和线路安全约束系统。

5.根据权利要求3所述间歇式发电并网条件下的电力系统日前鲁棒调度系统，其特征在于，所述环境成本模块包括发电环境成本S_e1和弃风/光环境成本S_e2，所述发电环境成本S_e1和弃风/光环境成本S_e2均包括排放物二氧化硫SO₂、氮氧化合物NO_x、粉尘TSP和二氧化碳CO₂，计算过程如下：

1)单位电量SO₂排放量计算公式为：

$Z_{{\mathrm{SO}}_2}=Q\frac{64}{32}{S}_{\mathrm{ar}}{b}_{{\mathrm{SO}}_2}(1-{\mathrm{\η}}_S)......\left(1\right);$

2)单位电量NO_x排放量计算公式为：

$Z_{\mathrm{NOx}}=Q\frac{30}{14}{N}_{\mathrm{ar}}{b}_{\mathrm{NOx}}(1-{\mathrm{\η}}_N)......\left(2\right);$

3)单位电量TSP排放量计算公式为：

Z_TSP＝Qb_TSP(1-η_TSP)……(3)；

4)单位电量CO₂排放成本计算公式为：

S＝QABP……(4)，

所述发电环境成本模块包括计算公式为：

$S_{e1}=M_{a1}(\underset{i}{\mathrm{\Σ}}0.6\frac{\mathrm{Zi}}{\mathrm{Ni}}+S)......\left(5\right),$

所述弃风/光环境成本计算公式为：

$S_{e2}=M_{a2}(\underset{i}{\mathrm{\Σ}}0.6\frac{\mathrm{Zi}}{\mathrm{Ni}}+S)......\left(6\right),$

所述风险成本模块包括风险成本S_r，所述风险成本S_r包括切负荷成本S_c和弃风/光成本S_a：

所述切负荷成本S_c计算公式为：

S_c＝M_cP_c……(7)，

总风险成本为弃风/光环境成本与切负荷成本之和,计算公式为：

S_f＝S_c+S_e2……(8)，

所述综合成本运算模块包括综合成本计算公式：

S_all＝S_gen+S_c+S_e+S_t……(9)，

所述风险成本模块包括统计风险成本，其中包括弃新能源量及切负荷量的计算公式为：

$M_{a2}=\underset{P_{\max }}{\overset{P_{\mathrm{cap}}}{\∫}}{P}_w(P_t-P_{\max })d{P}_t\∀(P_t>P_{\max })......\left(10\right),$

$M_c=\underset{0}{\overset{P_{\min }}{\∫}}{P}_w(P_{\min }-P_t)d{P}_t\∀(P_{\min }>P_t)......\left(11\right),$

所述目标函数模块为以预测场景下系统的发电成本作为目标函数，计算公式为：

min.F(g,t)＝S_n+S_q+S_x+S_un+S_uq+S_en……(12)；

所述约束方程包括：

(1)所述系统侧约束模块包括：

①所述功率平衡约束系统的计算公式为：

$\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{P}_g^t=D^t......\left(13\right);$

②所述旋转备用约束系统的计算公式为：

$\underset{g\∈G_{\mathrm{aa}}}{\mathrm{\Σ}}{I}_g^t{P}_{g\max }+\underset{g\∈G_{\mathrm{ab}}}{\mathrm{\Σ}}{I}_g^t{P}_g^t\≥D^t+R^t......\left(14\right);$

③所述线路安全约束系统的计算公式为：

$\left|\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{gj}}{P}_g^t\right|\≤L_j......\left(15\right);$

(2)所述机组侧约束模块包括：

①所述火电机组约束模块包括：

机组运行成本计算公式为：

$S_n=\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈G_n}{\mathrm{\Σ}}{A}_n\×{\left(P_n^t\right)}^2+B_n\×P_n^t+C_n......\left(16\right),$

机组提供旋转备用容量计算公式为：

$\underset{n\∈G_n}{\mathrm{\Σ}}{R}_n^t=\underset{n\∈G_n}{\mathrm{\Σ}}{I}_n^t(P_{n\max }-P_n^t)......\left(17\right),$

机组出力约束计算公式为：

$P_{n\min }\≤P_n^t\≤P_{n\max }......\left(18\right),$

机组爬坡率约束公式为：

$-r_d^n\≤P_n^t-P_n^{t-1}\≤r_u^n......\left(19\right),$

机组开关机约束计算公式为：：

$\underset{h=0}{\overset{T_u\left(n\right)-1}{\mathrm{\Σ}}}I(n,t+h)\≥(I(n,t)-I(n,t-1))\×\min (T_u\left(n\right),N_T-t+1)......\left(20\right),$

$\underset{h=0}{\overset{T_d\left(n\right)-1}{\mathrm{\Σ}}}(1-I(n,t+h))\≥(I(n,t-1)-I(n,t))\×\min (T_d\left(n\right),N_T-t+1)......\left(21\right);$ ②所述气电机组约束模块包括：

机组运行成本计算公式为：

$S_q=\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈G_q}{\mathrm{\Σ}}{A}_q\×{\left(P_q^t\right)}^2+B_q\×P_q^t+C_q......\left(22\right);$

机组提供旋转备用容量计算公式为：

$\underset{q\∈G_q}{\mathrm{\Σ}}{R}_q^t=\underset{q\∈G_q}{\mathrm{\Σ}}{I}_q^t(P_{q\max }-P_q^t)......\left(23\right),$

机组出力约束计算公式为：

$P_{q\min }\≤P_q^t\≤P_{q\max }......\left(24\right),$

机组爬坡率约束公式：

$-r_d^q\≤P_q^t-P_q^{t-1}\≤r_u^q......\left(25\right),$

机组开关机约束计算公式为：

$\underset{h=0}{\overset{T_u\left(q\right)-1}{\mathrm{\Σ}}}I(q,t+h)\≥(I(q,t)-I(q,t-1))\×\min (T_u\left(q\right),N_T-t+1)......\left(26\right),$

$\underset{h=0}{\overset{T_d\left(q\right)-1}{\mathrm{\Σ}}}(1-I(q,t+h))\≥(I(q,t-1)-I(q,t))\×\min (T_d\left(\text{q}\right),N_T-t+1)......\left(27\right);$

气量约束计算公式为：

$\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈\mathrm{Gq}}{\mathrm{\Σ}}{P}_q^t\≤Q_q......\left(28\right);$

③所述水电机组约束模块包括：

机组提供旋转备用容量计算公式为：

$\underset{s\∈G_s}{\mathrm{\Σ}}{R}_s^t=\underset{s\∈G_s}{\mathrm{\Σ}}{I}_s^t(P_{s\max }-P_s^t)......\left(29\right),$

机组出力约束计算公式为：

$P_{s\min }\≤P_s^t\≤P_{s\max }......\left(30\right),$

水量约束计算公式为：

$\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈G_s}{\mathrm{\Σ}}{P}_s^t\≤Q_s......\left(31\right);$

④所述核电机组约束模块包括：

机组运行成本计算公式为：

$S_x=\underset{t\∈T}{\mathrm{\Σ}}\underset{g\∈G_x}{\mathrm{\Σ}}{B}_x\×P_x^t+C_x......\left(32\right),$

机组提供旋转备用容量计算公式为：

$\underset{x\∈G_x}{\mathrm{\Σ}}{R}_x^t=\underset{x\∈G_x}{\mathrm{\Σ}}{I}_x^t(P_{x\max }-P_x^t)......\left(33\right),$

机组出力约束计算公式为：

$P_{x\min }\≤P_x^t\≤P_{x\max }......\left(34\right),$

机组爬坡率约束计算公式为：

$-r_d^x\≤P_x^t-P_x^{t-1}\≤r_u^x......\left(35\right),$

机组开关机约束计算公式为：

I(x,t)＝1……(36)；

⑤所述抽水蓄能机组约束模块包括：

机组提供旋转备用容量计算公式为：

$\underset{\mathrm{cg}\∈G_{\mathrm{cg}}}{\mathrm{\Σ}}{P}_{\mathrm{cg}\max }-\underset{\mathrm{cg}\∈G_{\mathrm{cg}}}{\mathrm{\Σ}}{P}_{\mathrm{cg}}^t\≤Y_0+\underset{\mathrm{cm}\∈G_{\mathrm{cm}}}{\mathrm{\Σ}}\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cm}}^t-\underset{\mathrm{cg}\∈G_{\mathrm{cg}}}{\mathrm{\Σ}}\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cg}}^t(t\∈T)......\left(37\right),$

$\underset{\mathrm{cg}\∈G_{\mathrm{cg}}}{\mathrm{\Σ}}{R}_c^t=\underset{\mathrm{cg}\∈G_{\mathrm{cg}}}{\mathrm{\Σ}}{I}_{\mathrm{cg}}^t(P_{\mathrm{cg}\max }-P_{\mathrm{cg}}^t)......\left(38\right),$

机组出力约束计算公式为：

$P_{\mathrm{cg}\min }\≤P_{\mathrm{cg}}^t\≤{+P}_{\mathrm{cg}\max }......\left(39\right),$

$P_{\mathrm{cm}}^t=P_{\mathrm{cmr}}......\left(40\right),$

机组状态约束计算公式为：

$I_{\mathrm{cm}}^t+I_{\mathrm{cg}}^t\≤1......\left(41\right),$

库容约束计算公式为：

$W_l\≤W_0+e\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cm}}^t-\underset{t=0}{\overset{t}{\mathrm{\Σ}}}{P}_{\mathrm{cg}}^t\≤W_u(t\∈T)......\left(42\right),$

库容设定约束计算公式为：

$e\underset{t\∈T}{\mathrm{\Σ}}{P}_{\mathrm{cm}}^t-\underset{t\∈T}{\mathrm{\Σ}}{P}_{\mathrm{cg}}^t=W_1-W_0......\left(43\right),$

所述极限场景约束系统为：

所述功率平衡约束模块计算公式为：

$\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{P}_{g,s}^t=D^t......\left(44\right),$

所述线路安全约束模块计算公式为：

$\left|\underset{g\∈G_{\mathrm{all}}}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{gj}}{P}_{g,s}^t\right|\≤L_j......\left(45\right),$

所述机组出力约束模块计算公式为：

$P_{g\min }\≤P_{g,s}^t\≤P_{g\max }(g\∈G_{\mathrm{all}})......\left(46\right),$

所述调整时间约束模块计算公式为：

$-\frac{1}{6}{r}_d^g\≤P_{g,s}^t-P_g^t\≤\frac{1}{6}{r}_u^g(g\∈G_{\mathrm{all}})......\left(47\right).$