CN102289566A - 独立运行模式下的微电网多时间尺度能量优化调度方法 - Google Patents

Abstract

本发明公开了一种独立运行模式下的微电网多时间尺度能量优化调度方法。将微电网的经济运行分为日前计划和实时调度两个阶段，在日前计划阶段，将一个调度周期分为24个时段，基于日前预测数据，建立日前机组启停优化计划模型；在实时调度阶段，将遵循日前计划的开停机结果，基于实时超短期预测数据，并实时监测储能单元的能量状态，根据净负荷大小及储能单元所处的不同能量状态区间采取不同的能量优化，以确定各可控型微电源的有功功率调度指令、卸荷功率指令及切负荷指令。本发明适用于微电网在独立运行模式下的能量优化调度，既可提高微电网运行的经济性和可靠性，又可使储能设备的能量状态处于安全的工作范围内，利于延长储能单元的使用寿命。

Description

独立运行模式下的微电网多时间尺度能量优化调度方法

技术领域

本发明涉及微电网规划、运行、调度领域，尤其涉及一种独立运行模式下的微电网多时间尺度能量优化调度方法。

背景技术

在能源需求和环境保护的双重压力下，以光伏和风力发电为代表的分布式发电(DG)技术获得了越来越多的重视与应用，上述DG***通常接入配电网实现并网运行，接入位置和运行策略大多由用户自由确定，且各DG之间处于分散无联络的自由状态。随着可再生能源发电***的增加，DG***日前增多，这种运行模式不公对原有电网的安全性构成威胁，而且也不得于整体能源综合利用效率的提高。

微电网概念的提出为DG运行提供了一个新的模式，微电网是指由分布式电源、储能装置、能量变换装置、相关负荷和监控***、保护装置汇集而成的小型发配电***，既可以与外部电网并网运行，也可以独立运行。通常情况下，微电网以联网模式运行，以增强微电网运行的灵活性和可靠性。然而在某些情况下微电网却不得不独立自治运行，比如，大电网故障时微电网与大电网断开而独立运行，偏远孤岛、牧场、边防等电网建设难度大、设备投资成本高的地区由于无大电网存在微电网只能独立自治运行。

微电网中可再生能源发电占据一定的比例，如光伏电池、风力发电机等，它们对环境和气候的变化比较敏感，其发电所依赖的光照资源和风资源具有随机波动特性，使得它们的出力是不稳定的，尤其当天气发生剧变时，它们的出力随之发生更大的改变。独立运行时，由于没有大电网的支撑，为平衡可再生能源发电的随机波动，改善电能质量，维持***稳定，一般会配备一定容量的储能设备。铅酸蓄电池因其较为低廉的价格、能满足应用于微电网时的功率密度需求和响应速度需求，被认为是最合适的储能设备。

目前，有关微电网经济运行方面的研究主要集中于并网运行模式下的研究，而在独立运行模式下的微电网经济运行研究却很少，目前尚无公认成熟的解决方案；现有少量的研究只是集中于微网实时经济运行优化调度，忽视了大时间尺度规划(如日前)对***运行经济性的影响；同时，现有的少量研究也都是对独立运行微电网进行统一化建模，而没有考虑独立运行微电网在通过其内部压频控制单元吸收间歇式能源发电的出力波动性与负荷的波动功率时会引起储能单元能量状态的频繁改变，甚至超出其安全能量状态范围，以至缩减储能单元使用寿命，增加微电网运行维护成本；另外，现有的少量研究也没考虑间歇式能源发电过剩与负荷过重等极端情况，而在现实微电网中又不得不考虑微电网的各种运行情况，以保证微电网安全、可靠及经济运行。

发明内容

针对现有技术的不足，本发明的目的在于提供了一种独立运行模式下的微电网多时间尺度能量优化调度方法。

本发明的目的是通过以下技术方案来实现的，它包括如下步骤：：

1)统计微电网运行历史数据，建立微电网内所有可控型微电源的成本-出力曲线的非线性函数，并将其分段线性化；

2)采集微电网负荷信息数据、气象信息数据，综合微电网运行的历史数据，对负荷/风能/太阳能进行未来一天的预测，得到未来一天内微电网的负荷/风能/太阳能预测数据；

3)将微电网未来一天内的经济运行分为24个时段，以微电网全天运行成本最小为目标函数，其中所有可控型微电源使用分段线性化模型，考虑微电网内部的各时段能量平衡、各设备元件的出力限制/爬坡率限制/开停机成本，基于步骤2)中的日前负荷/风能/太阳能预测数据，将此微电网日前计划问题构成一个混合整数线性规划问题的数学模型进行求解，得到各时段可控型微电源的机组日前启停优化计划方案；

4)在微电网实时运行过程中，以每15分钟为一调度周期，即将每小时划分为4个调度时段，全天划分为nT＝24*4＝96个调度时段，在每次调度时刻监测储能单元的能量状态SOS，采集微电网负荷信息数据、气象信息数据以，对负荷/风能/太阳能进行超短期预测，得到本调度时段内微电网的负荷/风能/太阳能预测数据；

5)根据步骤3)的机组日前启停优化计划方案得到当前时段处于开机状态的可控型微电源集合，确定处于开机状态的各可控型微电源的基点运行功率的上下限根据步骤4)得到的本调度时段内微电网的负荷/风能/太阳能预测数据确定净负荷功率大小；

6)根据步骤4)监测到的该调度时段储能单元的能量状态所处不同状态区间，以及步骤5)确定的不同净负荷功率大小，为独立运行模式下的微电网制定不同的能量优化策略，并建立相应的能量优化模型，通过模型求解得到该时段的微电网经济运行调度方案；

7)由步骤6)得到的微电网经济运行调度方案形成微电网调度指令，发布给微电网中的可控型微电源、可再生能源发电微电源、卸荷装置以及负荷的控制器，使得微电网在下一时段按照指定方式安全经济运行；

8)在下一调度时刻，判断是否达到第nT个时段，如果不是，则重复步骤4)，如果是，则重复步骤2)。

与背景技术相比，本发明具有的有益效果是：

(1)传统的独立运行微电网能量优化调度没有进行长时间的计划，本发明方法将微电网的经济运行分为日前计划和实时调度两个阶段，长时间尺度的日前计划能确保微电网运行的整体经济性，短时间尺度的实时调度考虑了微电网运行的实时运行情况，能兼顾微电网运行的安全性与经济性。

(2)传统的独立运行微电网能量优化调度都是对微电网建立统一化的经济模型，其只考虑了微电网运行的多数情况，而没考虑微电网运行的少数极端情况，本发明方法将储能单元的能量状态划分为多个区间，根据微电网实时运行时所处的不同储能状态与净负荷状态采取不同的能量优化策略，考虑了微电网运行的所有可能情况，提高了微电网运行的安全性与储能单元的使用寿命。

(3)在涉及的能量优化模型中引入了需求侧负荷响应与过剩能量的切除控制，确保微电网运行的安全稳定性。本发明方法将储能单元的能量状态划分为多个区间，通常情况下不对储能单元进行充放电功率调度，即指定其功率运行基点为0，而只是用作为压频控制单元吸纳微电网内的不平衡功率，减少了微电网对储能单元容量的需求，提高了微电网的投资成本与维护成本，搞高了经济性。

附图说明

图1是独立运行模式下的微电网多时间尺度能量优化调度方法流程图。

图2是独立运行模式下的微电网实时能量优化调度流程图a。

图3是独立运行模式下的微电网实时能量优化调度流程图b。

图4是独立运行模式下的微电网实时能量优化调度流程图c。

图5是实施例微电网结构图。

图6是可控型微电源成本-出力曲线分段线性化示意图。

图7是日前计划得到的可控型微电源日前启停优化计划安排。

图8是采用本发明得到的微电网运行结果。

具体实施方式

以下结合附图和实施例对本发明作进一步的说明。

如图1所示，本发明一种独立运行模式下的微电网多时间尺度能量优化调度方法，包括如下步骤：

1)统计微电网运行历史数据，建立微电网内所有可控型微电源的成本-出力曲线的非线性函数，并将其分段线性化；

2)采集微电网负荷信息数据、气象信息数据，综合微电网运行的历史数据，对负荷/风能/太阳能进行未来一天的预测，得到未来一天内微电网的负荷/风能/太阳能预测数据；

3)将微电网未来一天内的经济运行分为24个时段，以微电网全天运行成本最小为目标函数，其中所有可控型微电源使用分段线性化模型，考虑微电网内部的各时段能量平衡、各设备元件的出力限制/爬坡率限制/开停机成本，基于步骤2)中的日前负荷/风能/太阳能预测数据，将此微电网日前计划问题构成一个混合整数线性规划问题的数学模型进行求解，得到各时段可控型微电源的机组日前启停优化计划方案，

上述混合整数线性规划问题模型的数学模型为：

min f(x，u)

$s.t\left\{\begin{array}{cc}\begin{array}{c}h(x,u)=0\\ \underset{\‾}{g}\≤g(x,u)\≤\stackrel{\‾}{g}\end{array}& x\∈R,u\∈\left\{\mathrm{0,1}\right\}\end{array}\right.$

其中：

优化变量x，u定义为：

$x=\left\{\begin{array}{cc}\begin{array}{c}{P}_{\mathrm{Gi}}^t,D_{\mathrm{Gi}}^{t,k}\\ \mathrm{\Δ}{P}_{\mathrm{no}+}^t,{\mathrm{\ΔP}}_{\mathrm{no}-}^t\end{array}& i\∈S_G,t\∈S_T\end{array}\right.$

$u=u_{\mathrm{Gi}}^t,u_{\mathrm{Gi}*}^t,v_{\mathrm{Gi}}^{t,k},i\∈S_G,t\∈S_T$

目标函数f(x，u)定义为：

$f(x,u)=\underset{t\∈S_T}{\mathrm{\Σ}}(\underset{i\∈S_G}{\mathrm{\Σ}}(u_{\mathrm{Gi}}^t{A}_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}\left(F_{\mathrm{Gi}}^k{D}_{\mathrm{Gi}}^{t,k}\right)+K_{\mathrm{OMi}}{P}_{\mathrm{Gi}}^t+S_{\mathrm{Gi}}^{\mathrm{on}}{u}_{\mathrm{Gi}*}^t)+\mathrm{\σ}({\mathrm{\ΔP}}_{\mathrm{no}+}^t+{\mathrm{\ΔP}}_{\mathrm{no}-}^t))$

等式约束h(x，u)包括：

(1).功率平衡约束：

$\underset{i\∈S_G}{\mathrm{\Σ}}{P}_{\mathrm{Gi}}^t+\underset{i\∈S_I}{\mathrm{\Σ}}{P}_{\mathrm{Ii}}^t+{\mathrm{\ΔP}}_{\mathrm{no}+}^t-{\mathrm{\ΔP}}_{\mathrm{no}-}^t=\underset{i\∈S_L}{\mathrm{\Σ}}{P}_{\mathrm{Li}}^t,t\∈S_T$

(2).可控型微电源出力定义：

$P_{\mathrm{Gi}}^t=u_{\mathrm{Gi}}^t{B}_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{D}_{\mathrm{Gi}}^{t,k},i\∈S_G,t\∈S_T$

(3).可控型微电源分段运行归属标记位互斥条件：

$\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,k}=u_{\mathrm{Gi}}^t,i\∈S_G,t\∈S_T$

(4).可控型微电源最小开停机时间约束中的等式约束

$\underset{t=1}{\overset{G_i^{\mathrm{on}}}{\mathrm{\Σ}}}(1-u_{\mathrm{Gi}}^t)=0$ $G_i^{\mathrm{on}}=\min \{N_T,(\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{on}}}-T_{\mathrm{Gi}}^{\mathrm{on}})u_{\mathrm{Gi}}^0\}$

$\underset{t=1}{\overset{G_i^{\mathrm{off}}}{\mathrm{\Σ}}}{u}_{\mathrm{Gi}}^t=0$ $G_i^{\mathrm{off}}=\min \{N_T,(\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{off}}}-T_{\mathrm{Gi}}^{\mathrm{off}})(1-u_{\mathrm{Gi}}^0)\}$

不等式约束g(x，u)包括：

(1).可控型微电源分段出力值定义

$\underset{j=k+1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j}\≤\frac{D_{\mathrm{Gi}}^{t,k}}{B_{\mathrm{Gi}}^{k+1}-B_{\mathrm{Gi}}^k}\≤\underset{j=k}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j},i\∈S_G,t\∈S_T,k=1...L_{\mathrm{Gi}}$

(2).可控型微电源开始开机标记位定义

$u_{\mathrm{Gi}*}^t\≥u_{\mathrm{Gi}}^t-u_{\mathrm{Gi}}^{t-1},i\∈S_G,t\∈S_T$

(3).可控型微电源爬坡率约束

$\left\{\begin{array}{cc}\begin{array}{c}{P}_{\mathrm{Gi}}^t-P_{\mathrm{Gi}}^{t-1}\≤\mathrm{\ΔT}\stackrel{\‾}{{\mathrm{\ΔP}}_{\mathrm{Gi}}}\\ P_{\mathrm{Gi}}^{t-1}-P_{\mathrm{Gi}}^t\≤\mathrm{\ΔT}\underset{\‾}{\mathrm{\Δ}{P}_{\mathrm{Gi}}}\end{array}& i\∈S_G,t\∈S_T\end{array}\right.$

(4).可控型微电源最小开停机时间约束

$\left\{\begin{array}{cc}\begin{array}{cc}\underset{t=k}{\overset{k+\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{on}}}-1}{\mathrm{\Σ}}}{u}_{\mathrm{Gi}}^t\≥\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{on}}}(u_{\mathrm{Gi}}^k-u_{\mathrm{Gi}}^{k-1})& k\∈[G_i^{\mathrm{on}}+1,N_T-\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{on}}}+1]\\ \underset{t=k}{\overset{N_T}{\mathrm{\Σ}}}\left(u_{\mathrm{Gi}}^t(u_{\mathrm{Gi}}^k-u_{\mathrm{Gi}}^{k-1})\right)\≥0& k\∈[N_T-\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{on}}}+2,N_T\end{array}& i\∈S_G\end{array}\right.$

$\left\{\begin{array}{cc}\begin{array}{cc}\underset{t=k}{\overset{k+\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{off}}}-1}{\mathrm{\Σ}}}(1-u_{\mathrm{Gi}}^t)\≥\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{off}}}(u_{\mathrm{Gi}}^{k-1}-u_{\mathrm{Gi}}^k)& k\∈[G_i^{\mathrm{off}}+1,N_T-\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{off}}}+1]\\ \underset{t=k}{\overset{N_T}{\mathrm{\Σ}}}(1-u_{\mathrm{Gi}}^t-(u_{\mathrm{Gi}}^{k-1}-u_{\mathrm{Gi}}^k))\≥0& k\∈[N_T-\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{off}}}+2,N_T\end{array}& i\∈S_G\end{array}\right.$

(5).可控型微电源最大开停机次数约束

$\underset{t\∈S_T}{\mathrm{\Σ}}{u}_{\mathrm{Gi}*}^t\≤N_{\mathrm{on}}^{\max },i\∈S_G$

(6).参与压频控制的微电源或储能单元运行状态约束

$\underset{i\∈S_{\mathrm{Vf}}}{\mathrm{\Σ}}{u}_{\mathrm{Gi}}^t\≥1$

其中，各符号定义如下：S_T为时段集合、S_G为可控型微电源集合；S_I为不可控型微电源集合、S_L为内部负荷集合、S_Vf为参与压频控制的微电源或储能单元集合、f(x，u)为目标函数、N_T为总时段数、

为可控型微电源成本曲线参数、为可控型微电源分段曲线归属状态、

为可控型微电源分段曲线取值状态、L_Gi为可控型微电源成本曲线分段数、

为可控型微电源开机成本、K_OMi为可控型微电源运行维护成本、

ΔP_Gi为可控型微电源出力变化率界限、

为可控型微电源最短连续运行/停运时间、

为可控型微电源初始时刻已连续运行/停运时间、

为可控型微电源最大开关机次数、

为可控型微电源有功出力、

为可控型微电源工作状态(0关1开)、

为可控型微电源开始开机标记位、

为不可控型微电源出力、

为负荷功率、

为***总发电功率与总负荷功率之间的差值

为***总负荷功率与总发电功率之间的差值

其中，不等式约束g(x，u)中的(6).参与压频控制的微电源或储能单元运行状态约束，其压频控制是指，微电网在独立运行时，必须要有至少一个微电源或储能单元参与压频控制，以为微电网提供稳定的电压和频率，如果多个微电源或储能单元同时参与压频控制，则它们将通过下垂控制分摊微电网内的功率波动，其中，储能单元始终参与压频控制，部分可控型微电源也可参与压频控制，其余可控型微电源均为定有功无功控制，即PQ控制；

4)在微电网实时运行过程中，以每15分钟为一调度周期，即将每小时划分为4个调度时段，全天划分为nT＝24*4＝96个调度时段，在每次调度时刻监测储能单元的能量状态SOS，采集微电网负荷信息数据、气象信息数据以，对负荷/风能/太阳能进行超短期预测，得到本调度时段内微电网的负荷/风能/太阳能预测数据，

其中，储能单元的能量状态SOS是反映其剩余存储能量占其总容量比例的技术参数，定义为：

$\mathrm{SOS}=\frac{C_{\mathrm{net}}}{C}=1-\frac{3.6\×{10}^{-6}\×\∫\mathrm{Pdt}}{C}$

式中：C_net-储能单元剩余存储能量，kWh；

C-储能单元总容量，kWh；

P-单元的放电功率，W。

5)根据第3)步的机组日前启停优化计划方案得到当前时段处于开机状态的可控型微电源集合，确定处于开机状态的各可控型微电源的基点运行功率的上下限

根据第4)步得到的本调度时段内微电网的负荷/风能/太阳能预测数据确定净负荷功率大小，

处于开机状态的各可控型微电源的基点运行功率的上下限

的确定是按如下3个步骤进行：

A)确定该时段参与压频控制的单元总共需提供的上调或下调功率裕量为：

${\mathrm{\ΔP}}_{\mathrm{\Σ}}^t=e_I\·\underset{i\∈S_I}{\mathrm{\Σ}}{P}_{\mathrm{Ii}}^t+e_L\·\underset{i\∈S_L}{\mathrm{\Σ}}{P}_{\mathrm{Li}}^t$

-在该调度时段内参与压频控制的单元总共需提供的上调或下调功率裕量；

e_I-不可控微电源出力功率预测的最大误差；

e_L-负荷功率预测的最大误差；

-不可控型微电源出力；

-负荷功率；

B)确定处于开机状态的各压频控制单元需要提供的上调或下调功率裕量：

${\mathrm{\ΔP}}_{\mathrm{Gi}}^t={\mathrm{\ΔP}}_{\mathrm{\Σ}}^t\·\frac{\stackrel{\‾}{P_{\mathrm{Gi}}}}{\stackrel{\‾}{P_s}+\underset{i\∈S_{\mathrm{GVf}}^t}{\mathrm{\Σ}}\stackrel{\‾}{P_{\mathrm{Gi}}}}$

-该调度时刻内第i台处于开机状态的可控型微电源需要提供的上调或下调功率裕量；

-在该调度时段内参与压频控制的单元总共需提供的上调或下调功率裕量；

-第i台处于开机状态的可控型微电源的最大输出功率；

-储能单元的最大输出功率；

-在时刻处于开机状态的参与压频控制的可控型微电源集合；

C)确定处于开机状态的各可控型微电源的基点运行功率的上下限

对于参与压频控制的可控型微电源：

$\left\{\begin{array}{c}{P}_{\mathrm{Gi}}^{\min }=\underset{\‾}{P_{\mathrm{Gi}}}+{\mathrm{\ΔP}}_{\mathrm{Gi}}^t\\ P_{\mathrm{Gi}}^{\max }=\stackrel{\‾}{P_{\mathrm{Gi}}}-{\mathrm{\ΔP}}_{\mathrm{Gi}}^t\end{array}\right.$

对于不参与压频控制的可控型微电源：

$\left\{\begin{array}{c}{P}_{\mathrm{Gi}}^{\min }=\underset{\‾}{P_{\mathrm{Gi}}}\\ P_{\mathrm{Gi}}^{\max }=\stackrel{\‾}{P_{\mathrm{Gi}}}\end{array}\right.....;$

其中，第5)步中的净负荷功率指总负荷超短期预测功率减去不可控微电源总超短期预测输出功率，即

$P_{\mathrm{net}}=\underset{t\∈S_L}{\mathrm{\Σ}}{P}_{\mathrm{Li}}^t-\underset{i\∈S_I}{\mathrm{\Σ}}{P}_{\mathrm{Ii}}^t$

式中：P_net-净负荷功率；

-不可控型微电源出力；

-负荷功率；

S_I-不可控型微电源集合；

S_L-内部负荷集合。

6)根据第4)步监测到的该调度时段储能单元的能量状态所处不同状态区间，以及第5)步确定的不同净负荷功率大小，为独立运行模式下的微电网制定不同的能量优化策略，并建立相应的能量优化模型，通过模型求解得到该时段的微电网经济运行调度方案，

如图2～图4所示，制定不同的能量优化策略，是按如下步骤进行的：

a)判断净负荷功率P_net是否满足P_net≥0，如满足，则进行步骤b)～g)，如不满足，则转到步骤h)；

b)如满足P_net≥0，则判断净负荷功率P_net与第5)步确定的处于开机状态的各可控型微电源的基点运行功率的下限之和

的关系是否满足

如满足P_net≥0，不满足

则转到步骤d)，如满足P_net≥0，且满足则进一步判断是否满足SOS＜SOSx，SOS_max1为设定的储能单元的最大允许储能状态，如满足SOS＜SOS_max1，则确定处于开机状态的各可控型微电源的输出功率指令均取为

同时卸荷功率指令为

无切负荷指令，得到该时段的微电网经济运行调度方案；如满足P_net≥0、不满足SOS＜SOS_max1，则进行步骤c)；

c)当步骤b)中不满足SOS＜SOS_max1时，进一步判断是否满足SOS＞SOS_min，如果满足SOS＞SOS_min，则计算储能单元的允许充电功率P_chmax，并判断是否满足

如满足则确定处于开机状态的各可控型微电源的输出功率指令均取为

同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，如不满足则确定处于开机状态的各可控型微电源的输出功率指令均取为

同时卸荷功率指令为

无切负荷指令，得到该时段的微电网经济运行调度方案；如果不满足SOS＞SOS_min，则需要对储能单元充电，确定其充电功率P_ch1，

$P_{\mathrm{ch}1}=\min (\frac{(\frac{{\mathrm{SOS}}_{\max 2}+{\mathrm{SOS}}_{\min }}{2}-\mathrm{SOS})\·C_{\mathrm{stor}}}{\mathrm{\Δt}},P_{\mathrm{ch}\_\max })$

P_ch1-储能单元充电功率；

SOS-储能单元当前的储能状态；

C_stor-储能单元容量，kWh；

P_{ch_max}-储能单元的最大可充电功率，kW；

此时储能单元相当于负荷，进一步判断是否满足如满足

则确定处于开机状态的各可控型微电源的输出功率指令均取为

同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，如不满足

则建立优化模型，优化分配处于开机状态的各可控型微电源的输出功率指令，同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，其中，建立的优化模型为

min f(x)

$s.t\left\{\begin{array}{cc}\begin{array}{c}h(x,u)=0\\ \underset{\‾}{g}\≤g(x,u)\≤\stackrel{\‾}{g}\end{array}& x\∈R,u\∈\left\{\mathrm{0,1}\right\}\end{array}\right.$

其中：

优化变量x，u定义为：

$x=P_{\mathrm{Gi}}^t,D_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

$u=v_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

目标函数f(x)定义为：

$f\left(x\right)=\underset{i\∈S_G^t}{\mathrm{\Σ}}(A_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}\left(F_{\mathrm{Gi}}^k{D}_{\mathrm{Gi}}^{t,k}\right)+K_{\mathrm{OMi}}{P}_{\mathrm{Gi}}^t)$

等式约束h(x，u)包括：

(1).功率平衡约束：

$\underset{i\∈S_G^t}{\mathrm{\Σ}}{P}_{\mathrm{Gi}}^t=\underset{i\∈S_L}{\mathrm{\Σ}}{P}_{\mathrm{Li}}^t+P_{\mathrm{ch}1}-\underset{i\∈S_I}{\mathrm{\Σ}}{P}_{\mathrm{Ii}}^t,i\∈S_G^t,t\∈S_T$

(2).可控型微电源出力定义：

$P_{\mathrm{Gi}}^t=B_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{D}_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

(3).可控型微电源分段运行归属标记位互斥条件：

$\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,k}=1,i\∈S_G^t,t\∈S_T$

不等式约束g(x，u)包括：

(1).可控型微电源分段出力值定义

$\underset{j=k+1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j}\≤\frac{D_{\mathrm{Gi}}^{t,k}}{B_{\mathrm{Gi}}^{k+1}-B_{\mathrm{Gi}}^k}\≤\underset{j=k}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j},i\∈S_G^t,t\∈S_T,k=1...L_{\mathrm{Gi}}$

其中，各符号定义如下：

为本时段处于开机状态的可控型微电源集合；S_I为不可控型微电源集合、S_L为内部负荷集合、f(x)为目标函数、

为可控型微电源成本曲线参数、为可控型微电源分段曲线归属状态、

为可控型微电源分段曲线取值状态、L_Gi为可控型微电源成本曲线分段数、K_OMi为可控型微电源运行维护成本、

为可控型微电源有功出力、

为不可控型微电源出力、

为负荷功率；

d)当步骤b)中满足P_net≥0，不满足

时，进一步判断是否满足SOS＞SOS_max2，如果满足SOS＞SOS_max2，则计算储能单元至少可提供的放电功率P_dh1，

$P_{\mathrm{dh}1}=\min (\frac{(\mathrm{SOS}-\frac{{\mathrm{SOS}}_{\max 2}+{\mathrm{SOS}}_{\min }}{2})\·C_{\mathrm{stor}}}{\mathrm{\Δt}},P_{\mathrm{dh}\_\max })$

P_dh1-储能单元放电功率；

SOS-储能单元当前的储能状态；

C_stor-储能单元容量，kWh；

P_{dh_max}-储能单元的最大可放电功率，kW；

并转到步骤e)；如果不满足SOS＞SOS_max2，则进一步判断是否满足SOS＞SOS_min，如果满足SOS＞SOS_min，则令P_dh1＝0，同时转到步骤e)，如果不满足SOS＞SOS_min，则对储能单元以功率P_ch1充电，并转到步骤g)，其为P_ch1定义为：

$P_{\mathrm{dh}1}=\min (\frac{(\frac{{\mathrm{SOS}}_{\max 2}+{\mathrm{SOS}}_{\min }}{2}-\mathrm{SOS})\·C_{\mathrm{stor}}}{\mathrm{\Δt}},P_{\mathrm{dh}\_\max })$

P_ch1-储能单元充电功率；

SOS-储能单元当前的储能状态；

C_stor-储能单元容量，kWh；

P_{ch_max}-储能单元的最大可充电功率，kW；

e)判断是否满足

如满足

则确定处于开机状态的各可控型微电源的输出功率指令均取为

储能单元放电功率指令为同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，如不满足则进一步判断是否满足

如不满足

则转到步骤f)，如满足，则储能单元功率指令为放电P_dh1，并由各可控型微电源优化分配P_net-P_dh1，同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，其中，各可控型微电源优化分配P_net-P_dh1对应的优化模型如下：

min f(x，u)

$s.t\left\{\begin{array}{cc}\begin{array}{c}h(x,u)=0\\ \underset{\‾}{g}\≤g(x,u)\≤\stackrel{\‾}{g}\end{array}& x\∈R,u\∈\left\{\mathrm{0,1}\right\}\end{array}\right.$

其中：

优化变量x，u定义为：

$x=P_{\mathrm{Gi}}^t,D_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

$u=v_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

目标函数f(x，u)定义为：

$f\left(x\right)=\underset{i\∈S_G^t}{\mathrm{\Σ}}(A_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}\left(F_{\mathrm{Gi}}^k{D}_{\mathrm{Gi}}^{t,k}\right)+K_{\mathrm{OMi}}{P}_{\mathrm{Gi}}^t)$

等式约束h(x，u)包括：

(1).功率平衡约束：

$\underset{i\∈S_G^t}{\mathrm{\Σ}}{P}_{\mathrm{Gi}}^t=P_{\mathrm{net}}-P_{\mathrm{dh}1},i\∈S_G^t,t\∈S_T$

(2).可控型微电源出力定义：

$P_{\mathrm{Gi}}^t=B_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{D}_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

(3).可控型微电源分段运行归属标记位互斥条件：

$\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,k}=1,i\∈S_G^t,t\∈S_T$

不等式约束g(x，u)包括：

(1).可控型微电源分段出力值定义

$\underset{j=k+1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j}\≤\frac{D_{\mathrm{Gi}}^{t,k}}{B_{\mathrm{Gi}}^{k+1}-B_{\mathrm{Gi}}^k}\≤\underset{j=k}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j},i\∈S_G^t,t\∈S_T,k=1...L_{\mathrm{Gi}}$

其中，各符号定义如下：

为本时段处于开机状态的可控型微电源集合；S_I为不可控型微电源集合、S_L为内部负荷集合、f(x)为目标函数、

为可控型微电源成本曲线参数、

为可控型微电源分段曲线归属状态、

为可控型微电源分段曲线取值状态、L_Gi为可控型微电源成本曲线分段数、K_OMi为可控型微电源运行维护成本、

为可控型微电源有功出力、P_net为净负荷功率；

f)如步骤e)中不满足

则计算储能单元可提供的最大放电功率P_dh max，其计算式为：

$P_{\mathrm{dh}\max }=\min (\frac{\left(\frac{\mathrm{SOS}-{\mathrm{SOS}}_{\min }}{2}\right)\·C_{\mathrm{stor}}}{\mathrm{\Δt}},P_{\mathrm{dh}\_\max })$

P_dh max-储能单元放电功率，

SOS-储能单元当前的储能状态，

C_stor-储能单元容量，kWh，

P_{dh_max}-储能单元的最大允许放电功率，kW，

并进一步判断是否满足

如果满足则各可控型微电源输出功率指令均为

同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，如果不满足

则转到步骤g)；

g)建立负荷可中断优化模型，根据模型求解结果确定各可控型微电源的输出功率指令及负荷切除指令，同时卸荷功率指令为0得到该时段的微电网经济运行调度方案，其中负荷可中断优化模型如下：

max f(x，u)

$s.t\left\{\begin{array}{cc}\begin{array}{c}h(x,u)=0\\ \underset{\‾}{g}\≤g(x,u)\≤\stackrel{\‾}{g}\end{array}& x\∈R,u\∈\{0,1\}\end{array}\right.$

其中：

优化变量x，u定义为：

$\left\{\begin{array}{cc}\begin{array}{c}x=P_{\mathrm{Gi}}^t,D_{\mathrm{Gi}}^{t,k},P_{\mathrm{stor}}^t,D_{\mathrm{Stor}}^k\\ u=x_i^t,v_{\mathrm{Gi}}^{t,k},u_{\mathrm{stor}}^t\end{array}& i\∈S_G^t,t\∈S_T\end{array}\right.$

目标函数f(x，u)定义为：

$f(x,u)=\underset{i\∈S_L^t}{\mathrm{\Σ}}(p_i^t{x}_i^t-b_i^t\stackrel{\‾}{x_i^t})\·P_{\mathrm{Li}}^t+\underset{t\∈S_L-S_L^t}{\mathrm{\Σ}}\left(c{P}_{\mathrm{Li}}^t\right)$

$-(\underset{i\∈S_G^t}{\mathrm{\Σ}}(A_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}\left(F_{\mathrm{Gi}}^k{D}_{\mathrm{Gi}}^{t,k}\right)+K_{\mathrm{OMi}}{P}_{\mathrm{Gi}}^t)+f(u_{\mathrm{stor}}^t,P_{\mathrm{stor}}^t))$

其中，

为储能单元放电罚函数的线性化表示，定义为：

$f(u_{\mathrm{stor}}^t,P_{\mathrm{stor}}^t)=u_{\mathrm{sotr}}^t{A}_{\mathrm{Stor}}^1+\underset{k=1}{\overset{L_S}{\mathrm{\Σ}}}\left(F_{\mathrm{Stor}}^k{D}_{\mathrm{Stor}}^k\right)$

储能单元放电罚函数

设计为

$C\left(P_{\mathrm{stor}}^t\right)=\mathrm{\δ}\·P_{\mathrm{stor}}^t\·\mathrm{\Δt}$

δ＝a₁+a₂·dSOS+a₃·P_dh+a₄·dSOS·P_dh+a₅·dSOS²

dSOS＝SOS-SOS^min

等式约束条件h(x，u)包括：

(1).功率平衡约束：

$\underset{i\∈S_G^t}{\mathrm{\Σ}}{P}_{\mathrm{Gi}}^t+\underset{i\∈S_I}{\mathrm{\Σ}}{P}_{\mathrm{Ii}}^t+P_{\mathrm{Stor}}^t=\underset{i\∈s_L^t}{\mathrm{\Σ}}{x}_i^t{P}_{\mathrm{Li}}^t+\underset{i\∈S_L-S_L^t}{\mathrm{\Σ}}{P}_{\mathrm{Li}}^t,i\∈S_G^t,t\∈S_T$

(2).可控型微电源出力定义：

$P_{\mathrm{Gi}}^t=B_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{D}_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

(3).可控型微电源分段运行归属标记位互斥条件：

$\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,k}=1,i\∈S_G^t,t\∈S_T$

(4).储能单元出力定义：

$P_{\mathrm{Stor}}^t=u_{\mathrm{sotr}}^t{B}_{\mathrm{Stor}}^1+\underset{k=1}{\overset{L_S}{\mathrm{\Σ}}}{D}_S^k$

(5).储能单元分段运行归属标记位互斥条件：

$\underset{k=1}{\overset{L_S}{\mathrm{\Σ}}}{v}_{\mathrm{stor}}^k=u_{\mathrm{stor}}^t$

不等式约束g(x，u)包括：

(1).可控型微电源分段出力值定义：

$\underset{j=k+1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j}\≤\frac{D_{\mathrm{Gi}}^{t,k}}{B_{\mathrm{Gi}}^{k+1}-B_{\mathrm{Gi}}^k}\≤\underset{j=k}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j},i\∈S_G^t,t\∈S_T,k=1...L_{\mathrm{Gi}}$

(2).储能单元分段出力值定义：

$\underset{j=k+1}{\overset{L_S}{\mathrm{\Σ}}}{v}_{\mathrm{Stor}}^j\≤\frac{D_{\mathrm{Stor}}^k}{B_{\mathrm{Stor}}^{k+1}-B_{\mathrm{Stor}}^k}\≤\underset{j=k}{\overset{L_S}{\mathrm{\Σ}}}{v}_{\mathrm{Stor}}^j,k=1...L_{\mathrm{Stor}}$

其中，各符号定义如下：

为本时段处于开机状态的可控型微电源集合；S_I为不可控型微电源集合、S_L为内部负荷集合、

为内部可中断负荷集合、f(x，u)为目标函数、

是第i个可中断负荷与微电网运营商签订的合同电价(元/kWh)，

α_i是可中断负荷的电价系数，对于折扣式可中断负荷，α_i≤1，对于高赔偿可中断负荷，α_i＝1、p₀是售电电价(元/kWh)、

是可中断负荷被微电网中断后的单位赔偿费用(元/kWh)，b_i＝β_ip₀、β_i是中断赔偿系数，对于折扣式可中断负荷，β_i＝0，即没有中断后的赔偿费用、

为负荷、

是第i个可中断负荷的开断状态，1-未断开，0-断开、

表示取反；、

为可控型微电源成本曲线参数、为储能单元罚函数曲线参数、

为可控型微电源分段曲线归属状态、

为储能单元放电罚函数分段曲线归属状态、

为可控型微电源分段曲线取值状态、

为储能单元放电罚函数分段曲线取值状态、L_Gi为可控型微电源成本曲线分段数、L_S为储能单元放电罚函数曲线分段数、K_OMi为可控型微电源运行维护成本、

为可控型微电源有功出力；

h)当净负荷功率不满足P_net≥0时，进一步判断是否满足SOS＜SOS_max1，如不满足SOS＜SOS_max1，则各可控型微电源输出功率指令均为

同时卸荷功率指令为无切负荷指令，得到该时段的微电网经济运行调度方案，如满足SOS＜SOS_max1，则计算储能单元允许的最大充电功率P_ch max，并转到步骤i)，其中，P_ch max的计算式为：

$P_{\mathrm{ch}\max }=\min (\frac{({\mathrm{SOS}}_{\max 1}-\mathrm{SOS})\·C_{\mathrm{Stor}}}{\mathrm{\Δt}},P_{\mathrm{ch}\_\max })$

P_ch max-储能单元充电功率；

SOS-储能单元当前的储能状态；

C_stor-储能单元容量，kWh；

P_{ch_max}-储能单元的最大允许充电功率，kW；

i)判断是否满足

如满足

则各可控型微电源输出功率指令均为

同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案；如不满足

则各可控型微电源输出功率指令均为

同时卸荷功率指令为

无切负荷指令，得到该时段的微电网经济运行调度方案。

7)由第6)步得到的微电网经济运行调度方案形成微电网调度指令，发布给微电网中的可控型微电源、可再生能源发电微电源、卸荷装置以及负荷的控制器，使得微电网在下一时段按照指定方式安全经济运行；

8)在下一调度时刻，判断是否达到第nT个时段，如果不是，则重复步骤4)，如果是，则重复步骤2)。

以下结合附图，对本发明的实施例作详细说明，该发明的总体流程图如图1所示，该发明的第6)步所涉及到的流程图如图2～图4所示。

实施例：

考虑如图5所示的微电网，微电网内含柴油发电机(DE)、微型燃气轮机(MT)、燃料电池(FC)、风力发电机(WT)、光伏电池(PV)及蓄电池(Storage)，微电网与大电网连接的公共耦合点(PCC)保持断开，即微电网独立运行，其中，微电网中的柴油发电机、微型燃气轮机、燃料电池为可控型微电源，风力发电机、光伏电池为不可控型微电源，微电网由蓄电池与燃料电池共同参与压频控制，柴油发电机与微型燃气轮机采用PQ控制，蓄电池能量状态区间划分为SOS_amx1＝0.9，SOS_max2＝0.8，SOS_min＝0.4，采用本发明对独立运行模式下的微电网进行实时能量优化调度。

1)统计微电网运行历史数据，建立微电网内所有可控型微电源的成本-出力曲线的非线性函数，并将其分段线性化，分段线性化的形式如图6所示。以某型号的燃料电池为例，线性化后的参数如表1所示，使用分段线性化模型，即可以使用混合整数线性规划模型来建模微网的日前计划问题，保证问题的求解便捷。

表1某型号燃料电池成本-出力曲线分段线性化参数

2)采集微电网负荷信息数据、气象信息数据，综合微电网运行的历史数据，对负荷/风能/太阳能进行未来一天的预测，得到未来一天内微电网的负荷/风能/太阳能预测数据。

3)以1小时为一时段，将微电网未来一天内的经济运行分为24个时段，以微电网全天运行成本最小为目标函数，其中所有可控型微电源使用分段线性化模型，考虑微电网内部的各时段能量平衡、各设备元件的出力限制/爬坡率限制/开停机成本，基于步骤2)中的日前负荷/风能/太阳能预测数据，将此微电网日前计划问题构成一个混合整数线性规划问题进行求解，得到各时段可控型微电源的机组日前启停优化计划方案，如图7所示。

4)在微电网实时运行过程中，以每15分钟为一调度周期，即将每小时划分为4个调度时段，全天划分为nT＝24*4＝96个调度时段，在每次调度时刻监测储能单元的能量状态SOS，采集微电网负荷信息数据、气象信息数据以，对负荷/风能/太阳能进行超短期预测，得到本调度时段内微电网的负荷/风能/太阳能预测数据；

5)根据步骤3)的机组日前启停优化计划方案得到当前时段处于开机状态的可控型微电源集合，确定处于开机状态的各可控型微电源的基点运行功率的上下限

根据步骤4)得到的本调度时段内微电网的负荷/风能/太阳能预测数据确定净负荷功率大小；

6)根据步骤4)监测到的该调度时段储能单元的能量状态SOS所处不同状态区间，以及步骤5)确定的不同净负荷功率大小，为独立运行模式下的微电网制定不同的能量优化策略，并建立相应的能量优化模型，通过模型求解得到该时段的微电网经济运行调度方案；

7)由步骤6)得到的微电网经济运行调度方案形成微电网调度指令，发布给微电网中的可控型微电源、可再生能源发电微电源、卸荷装置以及负荷的控制器，使得微电网在下一时段按照指定方式安全经济运行；

8)在下一调度时刻，判断是否达到第24个时段，如果不是，则重复步骤4)，如果是，则重复步骤2)。

图8为采用本发明对独立运行模式下的微电网算例实时能量优化调度得到的结果，其中，图(a)超短期功率预测；(b)实际功率输出；(c)切负荷与卸荷；(d)可控型微电源功率输出指令；(e)压频控制单元的实际功率输出；(f)储能单元的能量状态SOS。

分析图7可知，在0-2点，由于微网内部负荷较少，故在日前计划只安排了一台柴油发电机(DE)，而在4-23点期间，用电负荷攀升，为尽可能满足负荷需求，安排三台可控型微电源全部开机，在24点，微网内部负荷有所下降，此段时间内安排了柴油发电机(DE)和燃料电池(FC)两台可控型微电源。

分析图8可知，图8(a)显示8-11点期间为微电网负荷最大时段，图8(c)显示在8-11点期间三台可控型微电源均按基点运行功率的上限值运行，由于仍不能全部满足负荷需求，因此该段时间内增加了蓄电池的输出功率(图8(e)可见)并通过需求侧负荷竞价策略切除了部分负荷(图8(c)P_Cut可见)，以维持微电网内部电能供需平衡；由于本微电网算例内部不可控微电源功率比例较小，所以全天都没出现功率过剩的情况，过剩功率卸荷指令一直为0(图8(c)中P_Shed)；图8(e)为压频控制单元的实际功率输出，因为压频控制单元将分摊微电网内部的随机性波动功率，所以燃料电池的实际功率输出(图8(e)P_FC)是在其功率运行指令值(图8(d)P_FC)基础上叠加了实时运行时的随机性波动功率分量；图8(e)P_Storage为蓄电池充放电功率，由于蓄电池参与压频控制，且在多数情况下不对其进行功率调度，即指定其功率运行基点为0，所以多数情况下其只分摊微电网内的波动功率，充放电功率较小，但在某些情况下，比如负荷过大时，若各微电源都按基点功率上限值运行仍不能完全满足所有负荷，则需要调度蓄电池放电，对应本算例在8-10点，体现为蓄电池有较大放电功率，但蓄电池在全天的充放电功率仍然相对较小，说明减少了微电网对储能单元容量的需求，提高了微电网的投资成本与维护成本，提高了微电网运行经济性。图8(f)为蓄电池能量状态SOS，可见其始终在安全能量状态范围0.5～0.9之间，防止过充电或过放电，能延长蓄电池使用寿命，保证微电网安全可靠的运行。

综上所述，通过本实施例的测试结果，说明本发明提出的一种独立运行模式下的微电网多时间尺度能量优化调度方法，能兼顾微电网运行的安全性与经济性，同时减小了微电网对储能单元容量的需求，且适应微电网负荷过重或新能源发电过剩等各种运行情况。

Claims (7)

1.一种独立运行模式下的微电网多时间尺度能量优化调度方法，其特征在于，包括如下步骤：

1)统计微电网运行历史数据，建立微电网内所有可控型微电源的成本-出力曲线的非线性函数，并将其分段线性化；

2)采集微电网负荷信息数据、气象信息数据，综合微电网运行的历史数据，对负荷/风能/太阳能进行未来一天的预测，得到未来一天内微电网的负荷/风能/太阳能预测数据；

3)将微电网未来一天内的经济运行分为24个时段，以微电网全天运行成本最小为目标函数，其中所有可控型微电源使用分段线性化模型，考虑微电网内部的各时段能量平衡、各设备元件的出力限制/爬坡率限制/开停机成本，基于步骤2)中的日前负荷/风能/太阳能预测数据，将此微电网日前计划问题构成一个混合整数线性规划问题的数学模型进行求解，得到各时段可控型微电源的机组日前启停优化计划方案；

4)在微电网实时运行过程中，以每15分钟为一调度周期，即将每小时划分为4个调度时段，全天划分为nT＝24*4＝96个调度时段，在每次调度时刻监测储能单元的能量状态SOS，采集微电网负荷信息数据、气象信息数据以，对负荷/风能/太阳能进行超短期预测，得到本调度时段内微电网的负荷/风能/太阳能预测数据；

5)根据步骤3)的机组日前启停优化计划方案得到当前时段处于开机状态的可控型微电源集合，确定处于开机状态的各可控型微电源的基点运行功率的上下限

根据步骤4)得到的本调度时段内微电网的负荷/风能/太阳能预测数据确定净负荷功率大小；

6)根据步骤4)监测到的该调度时段储能单元的能量状态所处不同状态区间，以及步骤5)确定的不同净负荷功率大小，为独立运行模式下的微电网制定不同的能量优化策略，并建立相应的能量优化模型，通过模型求解得到该时段的微电网经济运行调度方案；

7)由步骤6)得到的微电网经济运行调度方案形成微电网调度指令，发布给微电网中的可控型微电源、可再生能源发电微电源、卸荷装置以及负荷的控制器，使得微电网在下一时段按照指定方式安全经济运行；

8)在下一调度时刻，判断是否达到第nT个时段，如果不是，则重复步骤4)，如果是，则重复步骤2)。

2.根据权利要求1所述的一种独立运行模式下的微电网多时间尺度能量优化调度方法，其特征在于：所述步骤3)中的混合整数线性规划问题的数学模型为：

min f(x，u)

$s.t\left\{\begin{array}{cc}\begin{array}{c}h(x,u)=0\\ \underset{\‾}{g}\≤g(x,u)\≤\stackrel{\‾}{g}\end{array}& x\∈R,u\∈\left\{\mathrm{0,1}\right\}\end{array}\right.$

其中：

优化变量x，u定义为：

$x=\left\{\begin{array}{cc}\begin{array}{c}{P}_{\mathrm{Gi}}^t,D_{\mathrm{Gi}}^{t,k}\\ \mathrm{\Δ}{P}_{\mathrm{no}+}^t,{\mathrm{\ΔP}}_{\mathrm{no}-}^t\end{array}& i\∈S_G,t\∈S_T\end{array}\right.$

$u=u_{\mathrm{Gi}}^t,u_{\mathrm{Gi}*}^t,v_{\mathrm{Gi}}^{t,k},i\∈S_G,t\∈S_T$

目标函数f(x，u)定义为：

$f(x,u)=\underset{t\∈S_T}{\mathrm{\Σ}}(\underset{i\∈S_G}{\mathrm{\Σ}}(u_{\mathrm{Gi}}^t{A}_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}\left(F_{\mathrm{Gi}}^k{D}_{\mathrm{Gi}}^{t,k}\right)+K_{\mathrm{OMi}}{P}_{\mathrm{Gi}}^t+S_{\mathrm{Gi}}^{\mathrm{on}}{u}_{\mathrm{Gi}*}^t)+\mathrm{\σ}({\mathrm{\ΔP}}_{\mathrm{no}+}^t+{\mathrm{\ΔP}}_{\mathrm{no}-}^t))$

等式约束h(x，u)包括：

(1)功率平衡约束：

$\underset{i\∈S_G}{\mathrm{\Σ}}{P}_{\mathrm{Gi}}^t+\underset{i\∈S_I}{\mathrm{\Σ}}{P}_{\mathrm{Ii}}^t+{\mathrm{\ΔP}}_{\mathrm{no}+}^t-{\mathrm{\ΔP}}_{\mathrm{no}-}^t=\underset{i\∈S_L}{\mathrm{\Σ}}{P}_{\mathrm{Li}}^t,t\∈S_T$

(2)可控型微电源出力定义：

$P_{\mathrm{Gi}}^t=u_{\mathrm{Gi}}^t{B}_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{D}_{\mathrm{Gi}}^{t,k},i\∈S_G,t\∈S_T$

(3)可控型微电源分段运行归属标记位互斥条件：

$\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,k}=u_{\mathrm{Gi}}^t,i\∈S_G,t\∈S_T$

(4)可控型微电源最小开停机时间约束中的等式约束：

$\underset{t=1}{\overset{G_i^{\mathrm{on}}}{\mathrm{\Σ}}}(1-u_{\mathrm{Gi}}^t)=0,G_i^{\mathrm{on}}=\min \{N_T,(\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{on}}}-T_{\mathrm{Gi}}^{\mathrm{on}})u_{\mathrm{Gi}}^0\}$

$\underset{t=1}{\overset{G_i^{\mathrm{off}}}{\mathrm{\Σ}}}{u}_{\mathrm{Gi}}^t=0{,G}_i^{\mathrm{off}}=\min \{N_T,(\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{off}}}-T_{\mathrm{Gi}}^{\mathrm{off}})(1-u_{\mathrm{Gi}}^0)\}$

不等式约束g(x，u)包括：

(1)可控型微电源分段出力值定义：

$\underset{j=k+1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j}\≤\frac{D_{\mathrm{Gi}}^{t,k}}{B_{\mathrm{Gi}}^{k+1}-B_{\mathrm{Gi}}^k}\≤\underset{j=k}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j},i\∈S_G,t\∈S_T,k=1...L_{\mathrm{Gi}}$

(2)可控型微电源开始开机标记位定义：

$u_{\mathrm{Gi}*}^t\≥u_{\mathrm{Gi}}^t-u_{\mathrm{Gi}}^{t-1},i\∈S_G,t\∈S_T$

(3)可控型微电源爬坡率约束：

$\left\{\begin{array}{cc}\begin{array}{c}{P}_{\mathrm{Gi}}^t-P_{\mathrm{Gi}}^{t-1}\≤\mathrm{\ΔT}\stackrel{\‾}{{\mathrm{\ΔP}}_{\mathrm{Gi}}}\\ P_{\mathrm{Gi}}^{t-1}-P_{\mathrm{Gi}}^t\≤\mathrm{\ΔT}\underset{\‾}{\mathrm{\Δ}{P}_{\mathrm{Gi}}}\end{array}& i\∈S_G,t\∈S_T\end{array}\right.$

(4)可控型微电源最小开停机时间约束：

$\left\{\begin{array}{cc}\begin{array}{cc}\underset{t=k}{\overset{k+\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{on}}}-1}{\mathrm{\Σ}}}{u}_{\mathrm{Gi}}^t\≥\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{on}}}(u_{\mathrm{Gi}}^k-u_{\mathrm{Gi}}^{k-1})& k\∈[G_i^{\mathrm{on}}+1,N_T-\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{on}}}+1]\\ \underset{t=k}{\overset{N_T}{\mathrm{\Σ}}}\left(u_{\mathrm{Gi}}^t(u_{\mathrm{Gi}}^k-u_{\mathrm{Gi}}^{k-1})\right)\≥0& k\∈[N_T-\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{on}}}+2,N_T\end{array}& i\∈S_G\end{array}\right.$

$\left\{\begin{array}{cc}\begin{array}{cc}\underset{t=k}{\overset{k+\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{off}}}-1}{\mathrm{\Σ}}}(1-u_{\mathrm{Gi}}^t)\≥\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{off}}}(u_{\mathrm{Gi}}^{k-1}-u_{\mathrm{Gi}}^k)& k\∈[G_i^{\mathrm{off}}+1,N_T-\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{off}}}+1]\\ \underset{t=k}{\overset{N_T}{\mathrm{\Σ}}}(1-u_{\mathrm{Gi}}^t-(u_{\mathrm{Gi}}^{k-1}-u_{\mathrm{Gi}}^k))\≥0& k\∈[N_T-\underset{\‾}{T_{\mathrm{Gi}}^{\mathrm{off}}}+2,N_T\end{array}& i\∈S_G\end{array}\right.$

(5)可控型微电源最大开停机次数约束：

$\underset{t\∈S_T}{\mathrm{\Σ}}{u}_{\mathrm{Gi}*}^t\≤N_{\mathrm{on}}^{\max },i\∈S_G$

(6)参与压频控制的微电源或储能单元运行状态约束：

$\underset{i\∈S_{\mathrm{Vf}}}{\mathrm{\Σ}}{u}_{\mathrm{Gi}}^t\≥1$

其中，各符号定义如下：S_T为时段集合、S_G为可控型微电源集合；S_I为不可控型微电源集合、S_L为内部负荷集合、S_Vf为参与压频控制的微电源或储能单元集合、f(x，u)为目标函数、N_T为总时段数、

为可控型微电源成本曲线参数、为可控型微电源分段曲线归属状态、

为可控型微电源分段曲线取值状态、L_Gi为可控型微电源成本曲线分段数、

为可控型微电源开机成本、K_OMi为可控型微电源运行维护成本、

ΔP_Gi为可控型微电源出力变化率界限、为可控型微电源最短连续运行/停运时间、

为可控型微电源初始时刻已连续运行/停运时间、

为可控型微电源最大开关机次数、

为可控型微电源有功出力、为可控型微电源工作状态(0关1开)、

为可控型微电源开始开机标记位、

为不可控型微电源出力、

为负荷功率、

为***总发电功率与总负荷功率之间的差值

为***总负荷功率与总发电功率之间的差值

3.根据权利要求2所述的一种独立运行模式下的微电网多时间尺度能量优化调度方法，其特征在于：所述的不等式约束中(6)参与压频控制的微电源或储能单元运行状态约束，其压频控制是指，微电网在独立运行时，必须要有至少一个微电源或储能单元参与压频控制，以为微电网提供稳定的电压和频率，如果多个微电源或储能单元同时参与压频控制，则它们将通过下垂控制分摊微电网内的功率波动，其中，储能单元始终参与压频控制，部分可控型微电源也可参与压频控制，其余可控型微电源均为定有功无功控制，即PQ控制。

4.根据权利要求1所述的一种独立运行模式下的微电网多时间尺度能量优化调度方法，其特征在于：所述步骤4)的储能单元的能量状态SOS是反映其剩余存储能量占其总容量比例的技术参数，定义为：

$\mathrm{SOS}=\frac{C_{\mathrm{net}}}{C}=1-\frac{3.6\×{10}^{-6}\×\∫\mathrm{Pdt}}{C}$

式中：C_net-储能单元剩余存储能量，kWh；

C-储能单元总容量，kWh；

P-单元的放电功率，W。

5.根据权利要求1所述的一种独立运行模式下的微电网多时间尺度能量优化调度方法，其特征在于：所述步骤5)的确定处于开机状态的各可控型微电源的基点运行功率的上下限

是指包括如下步骤：

A)确定该时段参与压频控制的单元总共需提供的上调或下调功率裕量为：

${\mathrm{\ΔP}}_{\mathrm{\Σ}}^t=e_I\·\underset{i\∈S_I}{\mathrm{\Σ}}{P}_{I_i}^t+e_L\·\underset{i\∈S_L}{\mathrm{\Σ}}{P}_{\mathrm{Li}}^t$

-在该调度时段内参与压频控制的单元总共需提供的上调或下调功率裕量；

e_I-不可控微电源出力功率预测的最大误差；

e_L-负荷功率预测的最大误差；

-不可控型微电源出力；

-负荷功率；

B)确定处于开机状态的各压频控制单元需要提供的上调或下调功率裕量：

${\mathrm{\ΔP}}_{\mathrm{Gi}}^t={\mathrm{\ΔP}}_{\mathrm{\Σ}}^t\·\frac{\stackrel{\‾}{P_{\mathrm{Gi}}}}{\stackrel{\‾}{P_s}+\underset{i\∈S_{\mathrm{GVf}}^t}{\mathrm{\Σ}}\stackrel{\‾}{P_{\mathrm{Gi}}}}$

-该调度时刻内第i台处于开机状态的可控型微电源需要提供的上调或下调功率裕量；

-在该调度时段内参与压频控制的单元总共需提供的上调或下调功率裕量；

-第i台处于开机状态的可控型微电源的最大输出功率；

-储能单元的最大输出功率；

-在时刻处于开机状态的参与压频控制的可控型微电源集合；

C)确定处于开机状态的各可控型微电源的基点运行功率的上下限

对于参与压频控制的可控型微电源：

$\left\{\begin{array}{c}{P}_{\mathrm{Gi}}^{\min }=\underset{\‾}{P_{\mathrm{Gi}}}+{\mathrm{\ΔP}}_{\mathrm{Gi}}^t\\ P_{\mathrm{Gi}}^{\max }=\stackrel{\‾}{P_{\mathrm{Gi}}}-{\mathrm{\ΔP}}_{\mathrm{Gi}}^t\end{array}\right.$

对于不参与压频控制的可控型微电源：

$\left\{\begin{array}{c}{P}_{\mathrm{Gi}}^{\min }=\underset{\‾}{P_{\mathrm{Gi}}}\\ P_{\mathrm{Gi}}^{\max }=\stackrel{\‾}{P_{\mathrm{Gi}}}\end{array}\right..$

6.根据权利要求1所述的一种独立运行模式下的微电网多时间尺度能量优化调度方法，其特征在于：所述步骤5)的净负荷功率指总负荷超短期预测功率减去不可控微电源总超短期预测输出功率，即

$P_{\mathrm{net}}=\underset{t\∈S_L}{\mathrm{\Σ}}{P}_{\mathrm{Li}}^t-\underset{i\∈S_I}{\mathrm{\Σ}}{P}_{\mathrm{Ii}}^t$

P_net-净负荷功率；

-不可控型微电源出力；

-负荷功率；

S_I-不可控型微电源集合；

S_L-内部负荷集合。

7.根据权利要求1所述的一种独立运行模式下的微电网多时间尺度能量优化调度方法，其特征在于：所述步骤6)的为独立运行模式下的微电网制定不同的能量优化策略，并建立相应的能量优化模型，是指包括如下步骤：

a)判断净负荷功率P_net是否满足P_net≥0，如满足，则进行步骤b)～g)，如不满足，则转到步骤h)；

b)如满足P_net≥0，则判断净负荷功率P_net与由权利要求4所确定的处于开机状态的各可控型微电源的基点运行功率的下限之和

的关系是否满足

如满足P_net≥0，不满足

则转到步骤d)，如满足P_net≥0，且满足

则进一步判断是否满足SOS＜SOS_max1，SOS_max1为设定的储能单元的最大允许储能状态，如满足SOS＜SOS_max1，则确定处于开机状态的各可控型微电源的输出功率指令均取为

同时卸荷功率指令为

无切负荷指令，得到该时段的微电网经济运行调度方案；如满足P_net≥0、

不满足SOS＜SOS_max1，则进行步骤c)；

c)当步骤b)中不满足SOS＜SOS_max1时，进一步判断是否满足SOS＞SOS_min，如果满足SOS＞SOS_min，则计算储能单元的允许充电功率P_ch max，并判断是否满足

如满足则确定处于开机状态的各可控型微电源的输出功率指令均取为

同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，如不满足

则确定处于开机状态的各可控型微电源的输出功率指令均取为

同时卸荷功率指令为

无切负荷指令，得到该时段的微电网经济运行调度方案；如果不满足SOS＞SOS_min，则需要对储能单元充电，确定其充电功率P_ch1，

$P_{\mathrm{ch}1}=\min (\frac{(\frac{{\mathrm{SOS}}_{\max 2}+{\mathrm{SOS}}_{\min }}{2}-\mathrm{SOS})\·C_{\mathrm{stor}}}{\mathrm{\Δt}},P_{\mathrm{ch}\_\max })$

P_ch1-储能单元充电功率；

SOS-储能单元当前的储能状态；

C_stor-储能单元容量，kWh；

P_{ch_max}-储能单元的最大可充电功率，kW；

此时储能单元相当于负荷，进一步判断是否满足

如满足

则确定处于开机状态的各可控型微电源的输出功率指令均取为

同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，如不满足

则建立优化模型，优化分配处于开机状态的各可控型微电源的输出功率指令，同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，其中，建立的优化模型为

min f(x)

$s.t\left\{\begin{array}{cc}\begin{array}{c}h(x,u)=0\\ \underset{\‾}{g}\≤g(x,u)\≤\stackrel{\‾}{g}\end{array}& x\∈R,u\∈\left\{\mathrm{0,1}\right\}\end{array}\right.$

其中：

优化变量x，u定义为：

$x=P_{\mathrm{Gi}}^t,D_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

$u=v_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

目标函数定义为：

$f\left(x\right)=\underset{i\∈S_G^t}{\mathrm{\Σ}}(A_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}\left(F_{\mathrm{Gi}}^k{D}_{\mathrm{Gi}}^{t,k}\right)+K_{\mathrm{OMi}}{P}_{\mathrm{Gi}}^t)$

等式约束h(x，u)包括：

(1)功率平衡约束：

$\underset{i\∈S_G^t}{\mathrm{\Σ}}{P}_{\mathrm{Gi}}^t=\underset{i\∈S_L}{\mathrm{\Σ}}{P}_{\mathrm{Li}}^t+P_{\mathrm{ch}1}-\underset{i\∈S_I}{\mathrm{\Σ}}{P}_{\mathrm{Ii}}^t,i\∈S_G^t,t\∈S_T$

(2)可控型微电源出力定义：

$P_{\mathrm{Gi}}^t=B_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{D}_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

(3)可控型微电源分段运行归属标记位互斥条件：

$\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,k}=1,i\∈S_G^t,t\∈S_T$

不等式约束g(x，u)包括：

(1)可控型微电源分段出力值定义

$\underset{j=k+1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j}\≤\frac{D_{\mathrm{Gi}}^{t,k}}{B_{\mathrm{Gi}}^{k+1}-B_{\mathrm{Gi}}^k}\≤\underset{j=k}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j},i\∈S_G^t,t\∈S_T,k=1...L_{\mathrm{Gi}}$

其中，各符号定义如下：

为本时段处于开机状态的可控型微电源集合；S_I为不可控型微电源集合、S_L为内部负荷集合、f(x)为目标函数、

为可控型微电源成本曲线参数、

为可控型微电源分段曲线归属状态、

为可控型微电源分段曲线取值状态、L_Gi为可控型微电源成本曲线分段数、K_OMi为可控型微电源运行维护成本、

为可控型微电源有功出力、

为不可控型微电源出力、

为负荷功率；

d)当步骤b)中满足P_net≥0，不满足

时，进一步判断是否满足SOS＞SOS_max2，如果满足SOS＞SOS_max2，则计算储能单元至少可提供的放电功率P_dh1，

$P_{\mathrm{dh}1}=\min (\frac{(\mathrm{SOS}-\frac{{\mathrm{SOS}}_{\max 2}+{\mathrm{SOS}}_{\min }}{2})\·C_{\mathrm{stor}}}{\mathrm{\Δt}},P_{\mathrm{dh}\_\max })$

P_dh1-储能单元放电功率；

SOS-储能单元当前的储能状态；

C_stor-储能单元容量，kWh；

P_{dh_max}-储能单元的最大可放电功率，kW；

并转到步骤e)；如果不满足SOS＞SOS_max2，则进一步判断是否满足SOS＞SOS_min，如果满足SOS＞SOS_min，则令P_dh1＝0，同时转到步骤e)，如果不满足SOS＞SOS_min，则对储能单元以功率P_ch1充电，并转到步骤g)，其为P_ch1定义为：

$P_{\mathrm{dh}1}=\min (\frac{(\frac{{\mathrm{SOS}}_{\max 2}+{\mathrm{SOS}}_{\min }}{2}-\mathrm{SOS})\·C_{\mathrm{stor}}}{\mathrm{\Δt}},P_{\mathrm{dh}\_\max })$

P_ch1-储能单元充电功率；

SOS-储能单元当前的储能状态；

C_stor-储能单元容量，kWh；

P_{ch_max}-储能单元的最大可充电功率，kW；

e)判断是否满足

如满足

则确定处于开机状态的各可控型微电源的输出功率指令均取为

储能单元放电功率指令为

同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，如不满足

则进一步判断是否满足

如不满足

则转到步骤f)，如满足，则储能单元功率指令为放电P_dh1，并由各可控型微电源优化分配P_net-P_dh1，同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，其中，各可控型微电源优化分配P_net-P_dh1对应的优化模型如下：

min f(x，u)

$s.t\left\{\begin{array}{cc}\begin{array}{c}h(x,u)=0\\ \underset{\‾}{g}\≤g(x,u)\≤\stackrel{\‾}{g}\end{array}& x\∈R,u\∈\left\{\mathrm{0,1}\right\}\end{array}\right.$

其中：

优化变量x，u定义为：

$x=P_{\mathrm{Gi}}^t,D_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

$u=v_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

目标函数f(x，u)定义为：

$f\left(x\right)=\underset{i\∈S_G^t}{\mathrm{\Σ}}(A_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}\left(F_{\mathrm{Gi}}^k{D}_{\mathrm{Gi}}^{t,k}\right)+K_{\mathrm{OMi}}{P}_{\mathrm{Gi}}^t)$

等式约束h(x)包括：

(1).功率平衡约束：

$\underset{i\∈S_G^t}{\mathrm{\Σ}}{P}_{\mathrm{Gi}}^t=P_{\mathrm{net}}-P_{\mathrm{dh}1},i\∈S_G^t,t\∈S_T$

(2).可控型微电源出力定义：

$P_{\mathrm{Gi}}^t=B_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{D}_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

(3).可控型微电源分段运行归属标记位互斥条件：

$\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,k}=1,i\∈S_G^t,t\∈S_T$

不等式约束g(x，u)包括：

(1).可控型微电源分段出力值定义

$\underset{j=k+1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j}\≤\frac{D_{\mathrm{Gi}}^{t,k}}{B_{\mathrm{Gi}}^{k+1}-B_{\mathrm{Gi}}^k}\≤\underset{j=k}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j},i\∈S_G^t,t\∈S_T,k=1...L_{\mathrm{Gi}}$

其中，各符号定义如下：

为本时段处于开机状态的可控型微电源集合；S_I为不可控型微电源集合、S_L为内部负荷集合、f(x)为目标函数、

为可控型微电源成本曲线参数、为可控型微电源分段曲线归属状态、

为可控型微电源分段曲线取值状态、L_Gi为可控型微电源成本曲线分段数、K_OMi为可控型微电源运行维护成本、

为可控型微电源有功出力、

为净负荷功率；

f)如步骤e)中不满足

则计算储能单元可提供的最大放电功率P_dh max，其计算式为：

$P_{\mathrm{dh}\max }=\min (\frac{\left(\frac{\mathrm{SOS}-{\mathrm{SOS}}_{\min }}{2}\right)\·C_{\mathrm{stor}}}{\mathrm{\Δt}},P_{\mathrm{dh}\_\max })$

P_dh max-储能单元放电功率，

SOS-储能单元当前的储能状态，

C_stor-储能单元容量，kWh，

P_{dh_max}-储能单元的最大允许放电功率，kW，

并进一步判断是否满足

如果满足

则各可控型微电源输出功率指令均为

同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案，如果不满足

则转到步骤g)；

g)建立负荷可中断优化模型，根据模型求解结果确定各可控型微电源的输出功率指令及负荷切除指令，同时卸荷功率指令为0得到该时段的微电网经济运行调度方案，其中负荷可中断优化模型如下：

max f(x，u)

$s.t\left\{\begin{array}{cc}\begin{array}{c}h(x,u)=0\\ \underset{\‾}{g}\≤g(x,u)\≤\stackrel{\‾}{g}\end{array}& x\∈R,u\∈\{0,1\}\end{array}\right.$

其中：

优化变量x，u定义为：

$\left\{\begin{array}{cc}\begin{array}{c}x=P_{\mathrm{Gi}}^t,D_{\mathrm{Gi}}^{t,k},P_{\mathrm{stor}}^t,D_{\mathrm{Stor}}^k\\ u=x_i^t,v_{\mathrm{Gi}}^{t,k},u_{\mathrm{stor}}^t\end{array}& i\∈S_G^t,t\∈S_T\end{array}\right.$

目标函数f(x，u)定义为：

$f(x,u)=\underset{i\∈S_L^t}{\mathrm{\Σ}}(p_i^t{x}_i^t-b_i^t\stackrel{\‾}{x_i^t})\·P_{\mathrm{Li}}^t+\underset{t\∈S_L-S_L^t}{\mathrm{\Σ}}\left(c{P}_{\mathrm{Li}}^t\right)$

$-(\underset{i\∈S_G^t}{\mathrm{\Σ}}(A_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}\left(F_{\mathrm{Gi}}^k{D}_{\mathrm{Gi}}^{t,k}\right)+K_{\mathrm{OMi}}{P}_{\mathrm{Gi}}^t)+f(u_{\mathrm{stor}}^t,P_{\mathrm{stor}}^t))$

其中，

为储能单元放电罚函数的线性化表示，定义为：

$f(u_{\mathrm{stor}}^t,P_{\mathrm{stor}}^t)=u_{\mathrm{sotr}}^t{A}_{\mathrm{Stor}}^1+\underset{k=1}{\overset{L_S}{\mathrm{\Σ}}}\left(F_{\mathrm{Stor}}^k{D}_{\mathrm{Stor}}^k\right)$

储能单元放电罚函数

设计为

$C\left(P_{\mathrm{stor}}^t\right)=\mathrm{\δ}\·P_{\mathrm{stor}}^t\·\mathrm{\Δt}$

δ＝a₁+a₂·dSOS+a₃·P_dh+a₄·dSOS·P_dh+a₅·dSOS²

dSOS＝SOS-SOS^min

等式约束条件h(x，u)包括：

(1).功率平衡约束：

$\underset{i\∈S_G^t}{\mathrm{\Σ}}{P}_{\mathrm{Gi}}^t+\underset{i\∈S_I}{\mathrm{\Σ}}{P}_{\mathrm{Ii}}^t+P_{\mathrm{Stor}}^t=\underset{i\∈s_L^t}{\mathrm{\Σ}}{x}_i^t{P}_{\mathrm{Li}}^t+\underset{i\∈S_L-S_L^t}{\mathrm{\Σ}}{P}_{\mathrm{Li}}^t,i\∈S_G^t,t\∈S_T$

(2).可控型微电源出力定义：

$P_{\mathrm{Gi}}^t=B_{\mathrm{Gi}}^1+\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{D}_{\mathrm{Gi}}^{t,k},i\∈S_G^t,t\∈S_T$

(3).可控型微电源分段运行归属标记位互斥条件：

$\underset{k=1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,k}=1,i\∈S_G^t,t\∈S_T$

(4).储能单元出力定义：

$P_{\mathrm{Stor}}^t=u_{\mathrm{sotr}}^t{B}_{\mathrm{Stor}}^1+\underset{k=1}{\overset{L_S}{\mathrm{\Σ}}}{D}_S^k$

(5).储能单元分段运行归属标记位互斥条件：

$\underset{k=1}{\overset{L_S}{\mathrm{\Σ}}}{v}_{\mathrm{stor}}^k=u_{\mathrm{stor}}^t$

不等式约束g(x，u)包括：

(1).可控型微电源分段出力值定义：

$\underset{j=k+1}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j}\≤\frac{D_{\mathrm{Gi}}^{t,k}}{B_{\mathrm{Gi}}^{k+1}-B_{\mathrm{Gi}}^k}\≤\underset{j=k}{\overset{L_{\mathrm{Gi}}}{\mathrm{\Σ}}}{v}_{\mathrm{Gi}}^{t,j},i\∈S_G^t,t\∈S_T,k=1...L_{\mathrm{Gi}}$

(2).储能单元分段出力值定义：

$\underset{j=k+1}{\overset{L_S}{\mathrm{\Σ}}}{v}_{\mathrm{Stor}}^j\≤\frac{D_{\mathrm{Stor}}^k}{B_{\mathrm{Stor}}^{k+1}-B_{\mathrm{Stor}}^k}\≤\underset{j=k}{\overset{L_S}{\mathrm{\Σ}}}{v}_{\mathrm{Stor}}^j,k=1...L_{\mathrm{Stor}}$

其中，各符号定义如下：

为本时段处于开机状态的可控型微电源集合；S_I为不可控型微电源集合、SL为内部负荷集合、为内部可中断负荷集合、f(x，u)为目标函数、

是第i个可中断负荷与微电网运营商签订的合同电价(元/kWh)，

α_i是可中断负荷的电价系数，对于折扣式可中断负荷，α_i≤1，对于高赔偿可中断负荷，α_i＝1、p₀是售电电价(元/kWh)、

是可中断负荷被微电网中断后的单位赔偿费用(元/kWh)，b_i＝β_ip₀、β_i是中断赔偿系数，对于折扣式可中断负荷，β_i＝0，即没有中断后的赔偿费用、

为负荷、

是第i个可中断负荷的开断状态，1-未断开，0-断开、

表示取反；、

为可控型微电源成本曲线参数、

为储能单元罚函数曲线参数、

为可控型微电源分段曲线归属状态、

为储能单元放电罚函数分段曲线归属状态、

为可控型微电源分段曲线取值状态、

为储能单元放电罚函数分段曲线取值状态、L_Gi为可控型微电源成本曲线分段数、L_S为储能单元放电罚函数曲线分段数、K_OMi为可控型微电源运行维护成本、

为可控型微电源有功出力；

h)当净负荷功率不满足P_net≥0时，进一步判断是否满足SOS＜SOS_max1，如不满足SOS＜SOS_max1，则各可控型微电源输出功率指令均为

同时卸荷功率指令为

无切负荷指令，得到该时段的微电网经济运行调度方案，如满足SOS＜SOS_max1，则计算储能单元允许的最大充电功率P_ch max，并转到步骤i)，其中，P_ch max的计算式为：

$P_{\mathrm{ch}\max }=\min (\frac{({\mathrm{SOS}}_{\max 1}-\mathrm{SOS})\·C_{\mathrm{Stor}}}{\mathrm{\Δt}},P_{\mathrm{ch}\_\max })$

P_ch max-储能单元充电功率；

SOS-储能单元当前的储能状态；

C_stor-储能单元容量，kWh；

P_{ch_max}-储能单元的最大允许充电功率，kW；

i)判断是否满足

如满足

则各可控型微电源输出功率指令均为同时卸荷功率指令为0，无切负荷指令，得到该时段的微电网经济运行调度方案；如不满足

则各可控型微电源输出功率指令均为

同时卸荷功率指令为

无切负荷指令，得到该时段的微电网经济运行调度方案。

Images