CN105515005A - 基于开放市场环境下的多微电网系统优化方法 - Google Patents

Abstract

本发明提供基于开放市场环境下的多微电网系统优化方法，该方法包括以下步骤：收集数据；数据公布；下层优化；下层优化；上层协调；执行方案，各成员微电网的MGEMS接受并执行DMS下达的最终优化方案。本发明的效果是采用本方法在进行多微电网系统的优化时，以各微网自身利益为优化目标，微网可以自主制定自身的优化计划，增加自身所得，与不允许微电网间进行传输的多微电网系统相比，允许微电网间电能自由传输后，各微电网可以提升2％左右的所得。

Description

基于开放市场环境下的多微电网系统优化方法

技术领域

本发明属于微电网系统优化技术领域，本发明涉及一种基于开放市场环境下的多微电网系统优化方法。

背景技术

微电网是由分布式电源、储能系统、能量转换装置、监控和保护装置、负荷等汇集成的小型发、配、用电系统。将分布式电源以微电网的形式接入电网，可以充分发挥分布式电源的能源利用率，提高配电系统对分布式电源的接纳能力。随着分布式电源接入量的逐年提升，在一个局部配电系统中有可能会同时接入多个微电网，形成局部的多微电网系统。同时，我国目前正稳步推进电力体制改革，有序向社会资本开放电力市场，允许分布式电源用户和微电网系统以利益独立主体的身份参进行电能传输。可以预见，未来的微电网不仅可与所在配电系统进行电能传输，还可与临近的微电网进行电能传输，从而使含有多微电网系统的配电系统成为一个多利益主体的复杂网络，这对多微电网系统的优化提出了新的挑战。

目前，对多微电网系统优化的研究有许多，主要分为两类，一类研究较为初期，作为研究对象的微电网都被限定为只与配电系统进行电能传输，没有考虑微电网之间的电能传输，一般以微电网系统总体成本最低，或总体效益最高为优化目标，对分布式电源进行优化。另一类虽然考虑了微电网之间的电能传输，但优化过程中依然将多微电网系统视为一个整体，没有考虑各个微电网自身的利益。这与目前开放的市场环境，以及分布式电源和微电网的发展趋势并不相符，对分布式能源和微电网的有序发展缺乏有效的指导。

发明内容

针对上述问题，本发明的目的是提出一种基于开放市场环境的多微电网系统优化方法，可以解决多微电网间进行电能传输的问题，其中所关注的多微电网系统是指当多个微电网接入同一个中低压配电系统，并接受同一个配电管理机构调控时形成的复杂系统。

为实现上述目的，本发明采取的技术方案是提供基于开放市场环境下的多微电网系统优化方法，该方法建立在多个邻近的微电网接入同一个中压配电系统的多微电网系统结构上，多微电网系统中的微电网被称为成员微电网，每个微电网包含风电(WT)、光伏(PV)、微型燃气轮机(MT)、柴油发电机(DEG)和储能电池(Bat)中的一种或几种分布式电源，成员微电网通过微电网能量管理系统(MicroGridEnergyManageSystem，MGEMS)对自身内部可控发电单元的出力计划进行优化，并与所连接配电网的配电管理机构之间进行信息交互，多微电网系统中配电网的配电管理机构由配电管理系统(DistributionManagementSystem，DMS)担当，负责完成与各成员微电网之间的信息交换与共享、并协调微电网所制定的输电计划间的冲突，每个多微电网系统包含一个DMS和多个MGEMS；该方法包括以下步骤：

步骤1收集数据：多微电网系统中各成员微电网的MGEMS收集自身微电网内部未来24小时的负荷Load、光伏出力P^PV和风电出力P^WT信息，并制定自身在多微电网系统内进行电能传输费用p_sell，并将p_sell、Load、P^PV、P^WT传递给所连接配电网的DMS；

步骤2数据公布：配电网的DMS收集步骤1中各微网的负荷Load、光伏出力P^PV和风电出力P^WT信息，同时制定各微电网间进行电能传输的服务费用p_ser，并将服务费用p_ser连同步骤1中全部数据一起公布给所有成员微电网的MGEMS；

步骤3下层优化：各成员微电网的MGEMS以运行成本最低为目标函数，目标函数的数学表达式为：

$\max F_i^{MG}=\underset{t=1}{\overset{T}{\Σ}}({\mathrm{TR}}_{i,t}-F_{i,t}^{DEG}-F_{i,t}^{MT}-C_{i,t}^{OP}-{\mathrm{SC}}_{i,t})---\left(1\right)$

其中，TR_i,t表示微电网i在t时刻的电能传输成本，表示微电网i中柴油发电机在t时刻的发电燃料费用，表示微电网i中微型燃气轮机在t时刻的发电燃料费用，代表微电网i中各种分布式电源在t时刻的运行维护费用，SC_i,t表示微电网i中柴油发电机和微型燃气轮机的启动费用；

其中TR_i,t的计算公式可表示为：

${\mathrm{TR}}_{i,t}=\underset{j,j\≠i}{\Σ}{T}_t^{MG,ij}-T_{i,t}^D+p_{i,t}\·{\mathrm{Load}}_{i,t}---\left(2\right)$

$T_t^{MG,ij}=\left\{\begin{array}{c}{p_{sell}}^j\&CenterDot;P_t^{MG,ij}+p_{ser,ij}\&CenterDot;P_t^{MG,ij}(P_t^{MG,ij}<0)\\ {p_{sell}}^i\&CenterDot;P_t^{MG,ij}-p_{ser,ij}\&CenterDot;P_t^{MG,ij}(P_t^{MG,ij}\&GreaterEqual;0)\end{array}\right.---\left(3\right)$

$T_{i,t}^D=\left\{\begin{array}{cc}{p}_{SP}\&CenterDot;P_t^{D,i}& (P_t^{D,i}\&GreaterEqual;0)\\ p_{BP}\&CenterDot;P_t^{D,i}& (P_t^{D,i}<0)\end{array}\right.---\left(4\right)$

其中，式(2)中的表示微电网i在t时刻与在同一多微电网系统下除微电网i以外的成员微电网之间传输电能所需的传输费用和服务费用之和，表示微电网i在t时刻与配电系统的电能传输费用，p_i,t·Load_i,t表示微电网向自身负荷的供电所得，p_i,t为微网i在t时刻向用户供电的价格，Load_i,t为微网i在t时刻的负荷大小，式(3)中的p_sell ⁱ表示微网i的送电价格，p_ser,ij表示微网i与微网j间进行传输电能的服务费用价格，P_t ^MG,ij表示微电网i与微电网j之间的传输的电量，大于零时表示微电网i向微电网j输送的电量，小于零时表示微电网j向微电网i输送的电量，式(4)中p_SP表示配电网向微电网的送电价格，p_BP表示微电网向配电网的送电价格，P_t ^D,i表示微网i与配电网之间的传输的电量，大于零时表示配电系统向微电网i输送的电量，小于零时表示微电网i向配电系统输送的电量；

DEG的燃料费用表示为：

$F^{DEG}=\underset{m=1}{\overset{ND}{\&Sigma;}}(d^{DEG}+e^{DEG}\&CenterDot;P_m^{DEG}+f^{DEG}\&CenterDot;{\left(P_m^{DEG}\right)}^2)---\left(5\right)$

其中，F^DEG表示柴油发电机的发电成本，ND表示柴油发电机数目，d^DEG，e^DEG，f^DEG表示发电成本系数，由发电机性能和柴油价格决定，表示柴油发电机的有功出力；

MT的燃料费用表示为：

$F^{MT}=\underset{m=1}{\overset{NT}{\&Sigma;}}(\frac{C_{nl}}{L}\&times;\frac{{P_m}^{MT}}{\mathrm{\&eta;}})---\left(6\right)$

其中，F^MT表示微型燃气轮机的发电成本，NT表示微型燃气轮机数目，C_nl为天然气价格($/m3)，L为天然气热值(kWh/m³)，η表示燃气轮机的发电效率(％)，表示微型燃气轮机的有功出力；

分布式电源的运行维护费用表示为：

$C_{i,t}^{OP}=\underset{n}{\overset{ND+NT+NW+NP}{\&Sigma;}}(P_{i,t}^n\&times;{K_{OP}}^n)---\left(7\right)$

其中，NW表示风机个数，NP表示光伏个数，表示第n台机组的出力，K_OP ⁿ表示第n台机组的运行维护费用系数，光伏风电机组燃料成本为零；

DEG和MT的启动费用SC_i,t科表示为：

${\mathrm{SC}}_{i,t}=\underset{n}{\overset{ND+NT}{\&Sigma;}}\&lsqb;{\mathrm{hst}}_n+{\mathrm{cst}}_n\&CenterDot;\&lsqb;1-\exp \left(\frac{-{\mathrm{TF}}_t^n}{{\mathrm{CLT}}^n}\right)\&rsqb;\&rsqb;---\left(8\right)$

其中，hst_n表示机组热启动费用，cst_n表示机组冷启动费用，CLTⁿ表示机组冷却时间，TF_t ⁿ表示第n台机组到t时刻时处于关停状态下的持续时间；

下层优化过程中的约束条件包括负荷平衡约束、电量传输约束、爬坡约束和储能运行约束，具体表达式为：

1)微电网负荷平衡约束

$P_{i,t}^{Load}=\underset{m=1}{\overset{ND}{\&Sigma;}}{P}_{m,t}^{DEG}+\underset{n=1}{\overset{NT}{\&Sigma;}}{P}_{n,t}^{MT}+P_t^{WT}+P_t^{PV}+P_t^{Bat}+P_t^{D,i}+\underset{j,j\&NotEqual;i}{\&Sigma;}{P}_t^{MG,ji}---\left(9\right)$

其中，表示微电网i在t时刻的负荷；

2)电量传输约束

P_t ^MG,ji＝-P_t ^MG,ij(10)

$P_{t,min}^{MG,ij}\&le;P_t^{MG,ij}\&le;P_{t,max}^{MG,ij}$

$P_{t,min}^{MG,ij}=max\{P_{\min ,i,t}^{DEG}+P_{min,i,t}^{MT}+P_{i,t}^{WT}+P_{i,t}^{PV}-P_{i,t}^{Load},P_{j,t}^{Load}-P_{max,j,t}^{DEG}-P_{max,j,t}^{MT}-P_{j,t}^{WT}-P_{j,t}^{PV}\}$

$P_{t,\max }^{MG,ij}=\min \{P_{\max ,i,t}^{DEG}+P_{\max ,i,t}^{MT}+P_{i,t}^{WT}+P_{i,t}^{PV}-P_{i,t}^{Load},P_{j,t}^{Load}-P_{\min ,j,t}^{DEG}-P_{\min ,j,t}^{MT}-P_{j,t}^{WT}-P_{j,t}^{PV}\}---\left(11\right)$

其中，式(10)表示微电网输电矩阵是一个反对称矩阵，式(11)表示微电网i和j之间传输电量的上下限，和分别表示柴油发电机出力的上下限，和分别表示微型燃气轮机出力的上下限；

3)爬坡约束

$U_{m,t}^{DEG}{\mathrm{PL}}_{m,t}^{DEG}\&le;P_{m,t}^{DEG}\&le;U_{m,t}^{DEG}{\mathrm{PU}}_{m,t}^{DEG}---\left(12\right)$

${\mathrm{PL}}_{m,t}^{DEG}=max(P_{min,m}^{DEG},P_{m,t-1}^{DEG}-{\mathrm{RD}}_m^{DEG})---\left(13\right)$

${\mathrm{PU}}_{m,t}^{DEG}=min(P_{max,m}^{DEG},P_{m,t-1}^{DEG}+{\mathrm{RU}}_m^{DEG})---\left(14\right)$

其中，式(12)、(13)、(14)表示柴油发电机的爬坡约束，分别表示柴油发电机在t时刻的出力、出力下限和出力上限，分别表示柴油发电机额定出力的上下限，分别表示柴油发电机爬坡下限和上限，微型燃气轮机的约束以柴油发电机类推；

4)储能系统运行约束

$E_t^{Bat}=E_{t-1}^{Bat}-P_t^{Bat}---\left(15\right)$

$DPc\&le;P_t^{Bat}\&le;DPd---\left(16\right)$

$E_{\min }^{Bat}\&le;E_t^{Bat}\&le;E_{\max }^{Bat}---\left(17\right)$

其中，式(16)表示储能系统的充放电约束，DPc、DPd分别表示最大充电深度和最大放电深度，表示t时刻电池的剩余电量，式(17)是蓄电池容量约束，分别表示电池容量的上下限；

通过以上优化目标函数和约束条件得出自身内部柴油发电机出力计划P^DEG、微型燃气轮机出力计划P^MT、储能电池出力计划P^bat和与同一多微电网系统中除自身以外的成员微电网之间的电能传输计划P^MG,ij即表示微电网i向微电网j输送的电量，并传递给DMS；

步骤4上层协调：配电网的DMS检验步骤3中所有成员微电网的出力计划和电能传输计划，按照多微电网系统电能传输规则，解决各成员微电网电能传输计划中的冲突，并将传输计划P^MG,ij通过DMS与MGEMS之间的信息通道返回给MGEMS；

步骤5执行方案：各成员微电网的MGEMS接受并执行DMS下达的最终优化方案。

本发明的效果是采用本方法在进行多微电网系统的优化时，以各微网自身利益为优化目标，微网可以自主制定自身的优化计划，增加自身所得，与不允许微电网间进行传输的多微电网系统相比，允许微电网间电能自由传输后，各微电网可以提升2％左右的所得。

附图说明

图1是本发明中所用具体实例的多微电网结构示意图

图2是本发明中所用到的各微网日负荷曲线；

图3是本发明中所用到的各微网风电和光伏出力预测曲线；

图4是本发明中所用到的配电系统协调流程图；

图5是本发明的多微电网系统协调流程图。

具体实施方式

下面结合附图和实施例对本发明的基于开放市场环境下的多微电网系统优化方法进行详细的描述。

本发明的基于开放市场环境下的多微电网系统优化方法是建立在多个邻近的微电网接入同一个中压配电系统的多微电网系统结构上，多微电网系统中的微电网被称为成员微电网，每个微电网包含风电(WT)、光伏(PV)、微型燃气轮机(MT)、柴油发电机(DEG)和储能电池(Bat)中的一种或几种分布式电源，成员微电网通过微电网能量管理系统(MicroGridEnergyManageSystem，MGEMS)对自身内部可控发电单元的出力计划进行优化，并与所连接配电网的配电管理机构之间进行信息交互，多微电网系统中配电网的配电管理机构由配电管理系统(DistributionManagementSystem，DMS)担当，负责完成与各成员微电网之间的信息交换与共享、并协调微电网所制定的输电计划间的冲突，每个多微电网系统包含一个DMS和多个MGEMS；该方法包括以下步骤：

步骤1收集数据：多微电网系统中各成员微电网的MGEMS收集自身微电网内部未来24小时的负荷Load、光伏出力P^PV和风电出力P^WT信息，并制定自身在多微电网系统内进行电能传输费用p_sell，并将p_sell、Load、P^PV、P^WT传递给所连接配电网的DMS。

步骤2数据公布：配电网的DMS收集步骤1中各微网的负荷Load、光伏出力P^PV和风电出力P^WT信息，同时制定各微电网间进行电能传输的服务费用p_ser，并将服务费用p_ser连同步骤1中全部数据一起公布给所有成员微电网的MGEMS。

步骤3下层优化：各成员微电网的MGEMS以运行成本最低为目标函数，目标函数的数学表达式为：

$\max F_i^{MG}=\underset{t=1}{\overset{T}{\&Sigma;}}({\mathrm{TR}}_{i,t}-F_{i,t}^{DEG}-F_{i,t}^{MT}-C_{i,t}^{OP}-{\mathrm{SC}}_{i,t})---\left(1\right)$

其中，TR_i,t表示微电网i在t时刻的电能传输成本，表示微电网i中柴油发电机在t时刻的发电燃料费用，表示微电网i中微型燃气轮机在t时刻的发电燃料费用，代表微电网i中各种分布式电源在t时刻的运行维护费用，SC_i,t表示微电网i中柴油发电机和微型燃气轮机的启动费用。

其中TR_i,t的计算公式可表示为：

${\mathrm{TR}}_{i,t}=\underset{j,j\&NotEqual;i}{\&Sigma;}{T}_t^{MG,ij}-T_{i,t}^D+p_{i,t}\&CenterDot;{\mathrm{Load}}_{i,t}---\left(2\right)$

$T_t^{MG,ij}=\left\{\begin{array}{c}{p_{sell}}^j\&CenterDot;P_t^{MG,ij}+p_{ser,ij}\&CenterDot;P_t^{MG,ij}(P_t^{MG,ij}<0)\\ {p_{sell}}^i\&CenterDot;P_t^{MG,ij}-p_{ser,ij}\&CenterDot;P_t^{MG,ij}(P_t^{MG,ij}\&GreaterEqual;0)\end{array}\right.---\left(3\right)$

$T_{i,t}^D=\left\{\begin{array}{cc}{p}_{SP}\&CenterDot;P_t^{D,i}& (P_t^{D,i}\&GreaterEqual;0)\\ p_{BP}\&CenterDot;P_t^{D,i}& (P_t^{D,i}<0)\end{array}\right.---\left(4\right)$

其中，式(2)中的表示微电网i在t时刻与在同一多微电网系统下除微电网i以外的成员微电网之间传输电能所需的传输费用和服务费用之和，表示微电网i在t时刻与配电系统的电能传输费用，p_i,t·Load_i,t表示微电网向自身负荷的供电所得，p_i,t为微网i在t时刻向用户供电的价格，Load_i,t为微网i在t时刻的负荷大小，式(3)中的p_sell ⁱ表示微网i的送电价格，p_ser,ij表示微网i与微网j间进行传输电能的服务费用价格，P_t ^MG,ij表示微电网i与微电网j之间的传输的电量，大于零时表示微电网i向微电网j输送的电量，小于零时表示微电网j向微电网i输送的电量，式(4)中p_SP表示配电网向微电网的送电价格，p_BP表示微电网向配电网的送电价格，P_t ^D,i表示微网i与配电网之间的传输的电量，大于零时表示配电系统向微电网i输送的电量，小于零时表示微电网i向配电系统输送的电量。

DEG的燃料费用表示为：

$F^{DEG}=\underset{m=1}{\overset{ND}{\&Sigma;}}(d^{DEG}+e^{DEG}\&CenterDot;P_m^{DEG}+f^{DEG}\&CenterDot;{\left(P_m^{DEG}\right)}^2)---\left(5\right)$

其中，F^DEG表示柴油发电机的发电成本，ND表示柴油发电机数目，d^DEG，e^DEG，f^DEG表示发电成本系数，由发电机性能和柴油价格决定，表示柴油发电机的有功出力；

MT的燃料费用表示为：

$F^{MT}=\underset{m=1}{\overset{NT}{\&Sigma;}}(\frac{C_{nl}}{L}\&times;\frac{{P_m}^{MT}}{\mathrm{\&eta;}})---\left(6\right)$

其中，F^MT表示微型燃气轮机的发电成本，NT表示微型燃气轮机数目，C_nl为天然气价格($/m3)，L为天然气热值(kWh/m³)，η表示燃气轮机的发电效率(％)，表示微型燃气轮机的有功出力；

分布式电源的运行维护费用表示为：

$C_{i,t}^{OP}=\underset{n}{\overset{ND+NT+NW+NP}{\&Sigma;}}(P_{i,t}^n\&times;{K_{OP}}^n)---\left(7\right)$

其中，NW表示风机个数，NP表示光伏个数，表示第n台机组的出力，K_OP ⁿ表示第n台机组的运行维护费用系数，光伏风电机组燃料成本为零；

DEG和MT的启动费用SC_i,t科表示为：

${\mathrm{SC}}_{i,t}=\underset{n}{\overset{ND+NT}{\&Sigma;}}\&lsqb;{\mathrm{hst}}_n+{\mathrm{cst}}_n\&CenterDot;\&lsqb;1-\exp \left(\frac{-{\mathrm{TF}}_t^n}{{\mathrm{CLT}}^n}\right)\&rsqb;\&rsqb;---\left(8\right)$

其中，hst_n表示机组热启动费用，cst_n表示机组冷启动费用，CLTⁿ表示机组冷却时间，TF_t ⁿ表示第n台机组到t时刻时处于关停状态下的持续时间；

下层优化过程中的约束条件包括负荷平衡约束、电量传输约束、爬坡约束和储能运行约束，具体表达式为：

1)微电网负荷平衡约束

$P_{i,t}^{Load}=\underset{m=1}{\overset{ND}{\&Sigma;}}{P}_{m,t}^{DEG}+\underset{n=1}{\overset{NT}{\&Sigma;}}{P}_{n,t}^{MT}+P_t^{WT}+P_t^{PV}+P_t^{Bat}+P_t^{D,i}+\underset{j,j\&NotEqual;i}{\&Sigma;}{P}_t^{MG,ji}---\left(9\right)$

其中，表示微电网i在t时刻的负荷；

2)电量传输约束

P_t ^MG,ji＝-P_t ^MG,ij(10)

$P_{t,min}^{MG,ij}\&le;P_t^{MG,ij}\&le;P_{t,max}^{MG,ij}$

$P_{t,min}^{MG,ij}=max\{P_{\min ,i,t}^{DEG}+P_{min,i,t}^{MT}+P_{i,t}^{WT}+P_{i,t}^{PV}-P_{i,t}^{Load},P_{j,t}^{Load}-P_{max,j,t}^{DEG}-P_{max,j,t}^{MT}-P_{j,t}^{WT}-P_{j,t}^{PV}\}$

$P_{t,\max }^{MG,ij}=\min \{P_{\max ,i,t}^{DEG}+P_{\max ,i,t}^{MT}+P_{i,t}^{WT}+P_{i,t}^{PV}-P_{i,t}^{Load},P_{j,t}^{Load}-P_{\min ,j,t}^{DEG}-P_{\min ,j,t}^{MT}-P_{j,t}^{WT}-P_{j,t}^{PV}\}---\left(11\right)$

其中，式(10)表示微电网输电矩阵是一个反对称矩阵，式(11)表示微电网i和j之间传输电量的上下限，和分别表示柴油发电机出力的上下限，和分别表示微型燃气轮机出力的上下限。

3)爬坡约束

$U_{m,t}^{DEG}{\mathrm{PL}}_{m,t}^{DEG}\&le;P_{m,t}^{DEG}\&le;U_{m,t}^{DEG}{\mathrm{PU}}_{m,t}^{DEG}---\left(12\right)$

${\mathrm{PL}}_{m,t}^{DEG}=max(P_{min,m}^{DEG},P_{m,t-1}^{DEG}-{\mathrm{RD}}_m^{DEG})---\left(13\right)$

${\mathrm{PU}}_{m,t}^{DEG}=min(P_{max,m}^{DEG},P_{m,t-1}^{DEG}+{\mathrm{RU}}_m^{DEG})---\left(14\right)$

其中，式(12)、(13)、(14)表示柴油发电机的爬坡约束，分别表示柴油发电机在t时刻的出力、出力下限和出力上限，分别表示柴油发电机额定出力的上下限，分别表示柴油发电机爬坡下限和上限，微型燃气轮机的约束以柴油发电机类推。

4)储能系统运行约束

$E_t^{Bat}=E_{t-1}^{Bat}-P_t^{Bat}---\left(15\right)$

$DPc\&le;P_t^{Bat}\&le;DPd---\left(16\right)$

$E_{\min }^{Bat}\&le;E_t^{Bat}\&le;E_{\max }^{Bat}---\left(17\right)$

其中，式(16)表示储能系统的充放电约束，DPc、DPd分别表示最大充电深度和最大放电深度，表示t时刻电池的剩余电量，式(17)是蓄电池容量约束，分别表示电池容量的上下限。

通过以上优化目标函数和约束条件得出自身内部柴油发电机出力计划P^DEG、微型燃气轮机出力计划P^MT、储能电池出力计划P^bat和与同一多微电网系统中除自身以外的成员微电网之间的电能传输计划P^MG,ij即表示微电网i向微电网j输送的电量，并传递给DMS。

步骤4上层协调：配电网的DMS检验步骤3中所有成员微电网的出力计划和电能传输计划，按照多微电网系统电能传输规则，解决各成员微电网电能传输计划中的冲突，并将传输计划P^MG,ij通过DMS与MGEMS之间的信息通道返回给MGEMS。

步骤5执行方案：各成员微电网的MGEMS接受并执行DMS下达的最终优化方案。

上述步骤4中所述的多微电网系统电能传输规则具体为：

1)微电网需要首先满足自身的负荷需求，只有在自身负荷都满足的情况下，多余的电量才能够向外输送；

2)多微电网系统内电能传输费用p_sell应满足

p_BP<p_sell≤p_SP(18)

3)微电网间进行电能传输时需要配电系统提供线路容量、信息交换等配套服务，因此传输电能会产生服务成本，并由传输电能的双方承担；

4)当各微电网制定的与其他成员微电网的电能传输计划中出现冲突的情况时，缺电量最大的微电网拥有优先购电权。

3、根据权利要求1所述的基于开放市场环境下的多微电网系统优化方法，其特征是：步骤4中所述的电能传输计划冲突的解决方法，其具体步骤为：

步骤一：DMS读取各成员微电网的电量传输计划，对每个时段t生成各微电网间输电矩阵A，其中

A_ij＝P_t ^MG,ij(i≠j)(19)

A_ii＝0(20)

步骤二：比较矩阵的对应元素A_ij和A_ji，将绝对值较小的一个作为实际的输送电量，再将A中各元素去除实际的输送电量部分形成新的矩阵A’；

步骤三：求A’各行的和，若第i行和为正，则将微电网i放入多电微电网序列中，若第i行和为负，则将微电网i放入缺电微电网序列中，按照电量缺额从大到小对缺电微电网进行排序，按照送电电价从低到高对多电微电网进行排序；

步骤四：电量缺额排在第一位的缺电微电网优先按步骤三中对多电微电网排序的顺序向多电微电网引入电量，电量缺额完全满足后，排在下一位的缺电微电网才有资格按顺序向剩余的多电微电网引入电量，直到所有多电微电网的电量余额用尽或所有缺电微电网的电量缺额补足；

步骤五：若所有多电微电网的电量余额都平衡后，微电网依然还有剩余的电量缺额，则剩余的缺电量由配电系统补足，若所有缺电微电网的电量缺额都平衡后，微电网依然还有剩余的电量余额，则剩余的多电量全都输送给配电系统。

选取一个含三个微电网的配电系统，其结构如图1所示，各微电网中的分布式发电设备情况如表1所示。

表1各微电网分布式发电设备信息表

如图5所示，对系统进行如下处理：

第一步：收集数据：多微电网系统中各成员微电网的MGEMS收集自身微电网内部未来24小时的负荷Load、光伏出力P^PV和风电出力P^WT信息，其中各微电网内部的日负荷预测曲线如图2所示，风电和光伏预测出力情况如图3所示。各MGEMS制定自身在多微电网系统内进行电能传输的费用p_sell，其中MG1、MG2、MG3的电能传输费用分别为0.635元/kW·h、0.61元/kW·h、0.66元/kW·h。接下来各成员微电网的MGEMS将p_sell、Load、P^PV、P^WT传递给所连接配电网的DMS。

第二步：数据公布：配电网的DMS收集第一步中各微网的负荷Load、光伏出力P^PV和风电出力P^WT信息，同时制定各微电网间进行电能传输的服务费用p_ser为0.01元/kW·h。接下来DMS将服务费用p_ser连同第一步中全部数据一起公布给所有成员微电网的MGEMS。

第三步：下层优化：各成员微电网的MGEMS以运行成本最低为目标函数，以负荷平衡约束、电量传输约束、爬坡约束和储能运行约束为约束条件，其中，柴油发电机参数d＝0.4333，e＝0.2333，f＝0.0074，微型燃气轮机参数C_nl＝0.76元/m3，L＝9.7kW·h/m³，η的取值如下式所示：

$\mathrm{\&eta;}=0.0753{\left(\frac{P^{MT}}{65}\right)}^3-0.3095{\left(\frac{P^{MT}}{65}\right)}^2+0.4174\left(\frac{P^{MT}}{65}\right)+0.1068$

各类机组的运行维护费用系数如表3所示：

表3各类机组运行维护费用系数

MG1的MGEMS通过以上优化目标函数和约束条件得出自身内部柴油发电机出力计划deg1、微型燃气轮机出力计划mt1、储能电池出力计划bat1和与MG2、MG3之间的电能传输计划MG1_2、MG1_3，结果如表4所示。MG2的MGEMS通过以上优化目标函数和约束条件得出自身内部微型燃气轮机出力计划mt2、储能电池出力计划bat2和与MG1、MG3之间的电能传输计划MG2_1、MG2_3，结果如表5所示。MG3的MGEMS通过以上优化目标函数和约束条件得出自身内部柴油发电机出力计划deg3、储能电池出力计划bat3和与MG1、MG2之间的电能传输计划MG1_2、MG1_3，结果如表6所示。

表4MG1优化结果

表5MG2优化结果

表6MG3优化结果

接下来MG1将MG1_2、MG1_3，MG2将MG2_1、MG2_3，MG3将MG3_1、MG3_2传递给配电网的DMS。

第四步：上层协调：过程如图4所示。配电网的DMS检验第三步中的MG1_2、MG1_3、MG2_1、MG2_3、MG3_1、MG3_2，按照多微电网系统电能传输规则，解决各成员微电网电能传输计划中的冲突，对每个时段t，具体步骤如下：

步骤一：DMS读取t时段的MG1_2、MG1_3、MG2_1、MG2_3、MG3_1、MG3_2，生成各微电网间输电矩阵A，其中

$A=\left[\begin{array}{ccc}0& MG1\_2& MG1\_3\\ MG2\_1& 0& MG2\_3\\ MG3\_1& MG3\_2& 0\end{array}\right]$

步骤二：比较矩阵的对应元素A_ij和A_ji，将绝对值较小的一个作为实际的输送电量，再将A中各元素去除实际的输送电量部分形成新的矩阵A’；

步骤三：求A’各行的和，若第i行和为正，则将微电网i放入多电微电网序列中，若第i行和为负，则将微电网i放入缺电微电网序列中，按照电量缺额从大到小对缺电微电网进行排序，按照送电电价从低到高对多电微电网进行排序；

步骤四：电量缺额排在第一位的缺电微电网优先按步骤三中对多电微电网排序的顺序向多电微电网引入电量，电量缺额完全满足后，排在下一位的缺电微电网才有资格按顺序向剩余的多电微电网引入电量，直到所有多电微电网的电量余额用尽或所有缺电微电网的电量缺额补足；

步骤五：若所有多电微电网的电量余额都平衡后，微电网依然还有剩余的电量缺额，则剩余的缺电量由配电系统补足，若所有缺电微电网的电量缺额都平衡后，微电网依然还有剩余的电量余额，则剩余的多电量全都输送给配电系统。

对每个时段t重复以上5个步骤，得到最终的电能传输计划P^MG,ij，如表7所示：

表7最终电能传输计划

将最终的电能传输计划P^MG,ij通过DMS与MGEMS之间的信息通道返回给各成员微电网的MGEMS。

第五步：执行方案：各成员微电网的MGEMS接受并执行DMS下达的最终优化方案。

Claims (3)

1.一种基于开放市场环境下的多微电网系统优化方法，该方法建立在多个邻近的微电网接入同一个中压配电系统的多微电网系统结构上，多微电网系统中的微电网被称为成员微电网，每个微电网包含风电(WT)、光伏(PV)、微型燃气轮机(MT)、柴油发电机(DEG)和储能电池(Bat)中的一种或几种分布式电源，成员微电网通过微电网能量管理系统(MicroGridEnergyManageSystem，MGEMS)对自身内部可控发电单元的出力计划进行优化，并与所连接配电网的配电管理机构之间进行信息交互，多微电网系统中配电网的配电管理机构由配电管理系统(DistributionManagementSystem，DMS)担当，负责完成与各成员微电网之间的信息交换与共享、并协调微电网所制定的输电计划间的冲突，每个多微电网系统包含一个DMS和多个MGEMS；该方法包括以下步骤：

步骤1收集数据：多微电网系统中各成员微电网的MGEMS收集自身微电网内部未来24小时的负荷Load、光伏出力P^PV和风电出力P^WT信息，并制定自身在多微电网系统内进行电能传输费用p_sell，并将p_sell、Load、P^PV、P^WT传递给所连接配电网的DMS；

步骤2数据公布：配电网的DMS收集步骤1中各微网的负荷Load、光伏出力P^PV和风电出力P^WT信息，同时制定各微电网间进行电能传输的服务费用p_ser，并将服务费用p_ser连同步骤1中全部数据一起公布给所有成员微电网的MGEMS；

步骤3下层优化：各成员微电网的MGEMS以运行成本最低为目标函数，目标函数的数学表达式为：

${\mathrm{maxF}}_i^{MG}=\underset{t=1}{\overset{T}{\&Sigma;}}({\mathrm{TR}}_{i,t}-F_{i,t}^{DEG}-F_{i,t}^{MT}-C_{i,t}^{OP}-{\mathrm{SC}}_{i,t})---\left(1\right)$

其中，TR_i,t表示微电网i在t时刻的电能传输成本，表示微电网i中柴油发电机在t时刻的发电燃料费用，表示微电网i中微型燃气轮机在t时刻的发电燃料费用，代表微电网i中各种分布式电源在t时刻的运行维护费用，SC_i,t表示微电网i中柴油发电机和微型燃气轮机的启动费用；

其中TR_i,t的计算公式可表示为：

${\mathrm{TR}}_{i,t}=\underset{j,j\&NotEqual;i}{\&Sigma;}{T}_t^{MG,ij}-T_{i,t}^D+p_{i,t}\&CenterDot;{\mathrm{Load}}_{i,t}---\left(2\right)$

$T_t^{MG,ij}=\left\{\begin{array}{c}{p_{sell}}^j\&CenterDot;P_t^{MG,ij}+p_{ser,ij}\&CenterDot;P_t^{MG,ij}(P_t^{MG,ij}<0)\\ {p_{sell}}^j\&CenterDot;P_t^{MG,ij}+p_{ser,ij}\&CenterDot;P_t^{MG,ij}(P_t^{MG,ij}\&GreaterEqual;0)\end{array}\right.---\left(3\right)$

$T_{i,t}^D=\left\{\begin{array}{c}{p}_{SP}\&CenterDot;P_t^{D,i}(P_t^{D,i}\&GreaterEqual;0)\\ p_{BP}\&CenterDot;P_t^{D,i}(P_t^{D,i}<0)\end{array}\right.---\left(4\right)$

其中，式(2)中的表示微电网i在t时刻与在同一多微电网系统下除微电网i以外的成员微电网之间传输电能所需的传输费用和服务费用之和，表示微电网i在t时刻与配电系统的电能传输费用，p_i,t·Load_i,t表示微电网向自身负荷的供电所得，p_i,t为微网i在t时刻向用户供电的价格，Load_i,t为微网i在t时刻的负荷大小，式(3)中的p_sell ⁱ表示微网i的送电价格，p_ser,ij表示微网i与微网j间进行传输电能的服务费用价格，表示微电网i与微电网j之间的传输的电量，大于零时表示微电网i向微电网j输送的电量，小于零时表示微电网j向微电网i输送的电量，式(4)中p_SP表示配电网向微电网的送电价格，p_BP表示微电网向配电网的送电价格，表示微网i与配电网之间的传输的电量，大于零时表示配电系统向微电网i输送的电量，小于零时表示微电网i向配电系统输送的电量；

DEG的燃料费用表示为：

$F^{DEG}=\underset{m=1}{\overset{ND}{\&Sigma;}}(d^{DEG}+e^{DEG}\&CenterDot;P_m^{DEG}+f^{DEG}\&CenterDot;{\left(P_m^{DEG}\right)}^2)---\left(5\right)$

其中，F^DEG表示柴油发电机的发电成本，ND表示柴油发电机数目，d^DEG，e^DEG，f^DEG表示发电成本系数，由发电机性能和柴油价格决定，表示柴油发电机的有功出力；

MT的燃料费用表示为：

$F^{MT}=\underset{m=1}{\overset{NT}{\&Sigma;}}(\frac{C_{nl}}{L}\&times;\frac{{P_m}^{MT}}{\mathrm{\&eta;}})---\left(6\right)$

其中，F^MT表示微型燃气轮机的发电成本，NT表示微型燃气轮机数目，C_nl为天然气价格($/m3)，L为天然气热值(kWh/m³)，η表示燃气轮机的发电效率(％)，表示微型燃气轮机的有功出力；

分布式电源的运行维护费用表示为：

$C_{i,t}^{OP}=\underset{n}{\overset{ND+NT+NW+NP}{\&Sigma;}}(P_{i,t}^n\&times;{K_{OP}}^n)---\left(7\right)$

其中，NW表示风机个数，NP表示光伏个数，表示第n台机组的出力，K_OP ⁿ表示第n台机组的运行维护费用系数，光伏风电机组燃料成本为零；

DEG和MT的启动费用SC_i,t科表示为：

${\mathrm{SC}}_{i,t}=\underset{n}{\overset{ND+NT}{\&Sigma;}}\&lsqb;{\mathrm{hst}}_n+{\mathrm{cst}}_n\&CenterDot;\&lsqb;1-\exp \left(\frac{-{\mathrm{TF}}_t^n}{{\mathrm{CLT}}^n}\right)\&rsqb;\&rsqb;---\left(8\right)$

其中，hst_n表示机组热启动费用，cst_n表示机组冷启动费用，CLTⁿ表示机组冷却时间，表示第n台机组到t时刻时处于关停状态下的持续时间；

下层优化过程中的约束条件包括负荷平衡约束、电量传输约束、爬坡约束和储能运行约束，具体表达式为：

1)微电网负荷平衡约束

$P_{i,t}^{Load}=\underset{m=1}{\overset{ND}{\&Sigma;}}{P}_{m,t}^{DEG}+\underset{n=1}{\overset{NT}{\&Sigma;}}{P}_{n,t}^{MT}+P_t^{WT}+P_t^{PV}+P_t^{Bat}+P_t^{D,i}+\underset{j,j\&NotEqual;i}{\&Sigma;}{P}_t^{MG,ji}---\left(9\right)$

其中，表示微电网i在t时刻的负荷；

2)电量传输约束

$P_t^{MG,ji}=-P_t^{MG,ij}---\left(10\right)$

$\begin{array}{c}{P}_{t,\min }^{MG,ij}\&le;P_t^{MG,ij}\&le;P_{t,\max }^{MG,ij}\\ P_{t,\min }^{MG,ij}=\max \left\{P_{\min ,i,t}^{DEG}+P_{\min ,i,t}^{MT}+P_{i,t}^{WT}+P_{i,t}^{PV}-P_{i,t}^{Load},P_{j,t}^{Load}-P_{\max ,j,t}^{DEG}-P_{\max ,j,t}^{MT}-P_{j,t}^{WT}-P_{j,t}^{PV}\right\}\\ P_{t,\max }^{MG,ij}=\min \left\{P_{\max ,i,t}^{DEG}+P_{\max ,i,t}^{MT}+P_{i,t}^{WT}+P_{i,t}^{PV}-P_{i,t}^{Loda},P_{j,t}^{Load}-P_{\max ,j,t}^{DEG}-P_{\max ,j,t}^{MT}-P_{j,t}^{WT}-P_{j,t}^{PV}\right\}\end{array}---\left(11\right)$

其中，式(10)表示微电网输电矩阵是一个反对称矩阵，式(11)表示微电网i和j之间传输电量的上下限，和分别表示柴油发电机出力的上下限，和分别表示微型燃气轮机出力的上下限；

3)爬坡约束

$U_{m,t}^{DEG}{\mathrm{PL}}_{m,t}^{DEG}\&le;P_{m,t}^{DEG}\&le;U_{m,t}^{DEG}{\mathrm{PU}}_{m,t}^{DEG}---\left(12\right)$

${\mathrm{PL}}_{m,t}^{DEG}=max(P_{min,m}^{DEG},P_{m,t-1}^{DEG}-{\mathrm{RD}}_m^{DEG})---\left(13\right)$

${\mathrm{PU}}_{m,t}^{DEG}=min(P_{max,m}^{DEG},P_{m,t-1}^{DEG}+{\mathrm{RU}}_m^{DEG})---\left(14\right)$

其中，式(12)、(13)、(14)表示柴油发电机的爬坡约束，分别表示柴油发电机在t时刻的出力、出力下限和出力上限，分别表示柴油发电机额定出力的上下限，分别表示柴油发电机爬坡下限和上限，微型燃气轮机的约束以柴油发电机类推；

4)储能系统运行约束

$E_t^{Bat}=E_{t-1}^{Bat}-P_t^{Bat}---\left(15\right)$

$\mathrm{DPc}\&le;P_t^{\mathrm{Bat}}\&le;\mathrm{DPd}---\left(16\right)$

$E_{min}^{Bat}\&le;E_t^{Bat}\&le;E_{\max }^{Bat}---\left(17\right)$

其中，式(16)表示储能系统的充放电约束，DPc、DPd分别表示最大充电深度和最大放电深度，表示t时刻电池的剩余电量，式(17)是蓄电池容量约束，分别表示电池容量的上下限；

通过以上优化目标函数和约束条件得出自身内部柴油发电机出力计划P^DEG、微型燃气轮机出力计划P^MT、储能电池出力计划P^bat和与同一多微电网系统中除自身以外的成员微电网之间的电能传输计划P^MG,ij即表示微电网i向微电网j输送的电量，并传递给DMS；

步骤4上层协调：配电网的DMS检验步骤3中所有成员微电网的出力计划和电能传输计划，按照多微电网系统电能传输规则，解决各成员微电网电能传输计划中的冲突，并将传输计划P^MG,ij通过DMS与MGEMS之间的信息通道返回给MGEMS；

步骤5执行方案：各成员微电网的MGEMS接受并执行DMS下达的最终优化方案。

2.根据权利要求1所述的基于开放市场环境下的多微电网系统优化方法，其特征是：步骤4中所述的多微电网系统电能传输规则具体为：

1)微电网需要首先满足自身的负荷需求，只有在自身负荷都满足的情况下，多余的电量才能够向外输送；

2)多微电网系统内电能传输费用p_sell应满足

p_BP<p_sell≤p_SP(18)

3)微电网间进行电能传输时需要配电系统提供线路容量、信息交换等配套服务，因此传输电能会产生服务成本，并由传输电能的双方承担；

4)当各微电网制定的与其他成员微电网的电能传输计划中出现冲突的情况时，缺电量最大的微电网拥有优先购电权。

3.根据权利要求1所述的基于开放市场环境下的多微电网系统优化方法，其特征是：步骤4中所述的电能传输计划冲突的解决方法，其具体步骤为：

步骤一：DMS读取各成员微电网的电量传输计划，对每个时段t生成各微电网间输电矩阵A，其中

$A_{\mathrm{ij}}=P_t^{\mathrm{MG},\mathrm{ij}}(i\&NotEqual;j)---\left(19\right)$

A_ii＝0(20)

步骤二：比较矩阵的对应元素A_ij和A_ji，将绝对值较小的一个作为实际的输送电量，再将A中各元素去除实际的输送电量部分形成新的矩阵A’；

步骤三：求A’各行的和，若第i行和为正，则将微电网i放入多电微电网序列中，若第i行和为负，则将微电网i放入缺电微电网序列中，按照电量缺额从大到小对缺电微电网进行排序，按照送电电价从低到高对多电微电网进行排序；

步骤四：电量缺额排在第一位的缺电微电网优先按步骤三中对多电微电网排序的顺序向多电微电网引入电量，电量缺额完全满足后，排在下一位的缺电微电网才有资格按顺序向剩余的多电微电网引入电量，直到所有多电微电网的电量余额用尽或所有缺电微电网的电量缺额补足；

步骤五：若所有多电微电网的电量余额都平衡后，微电网依然还有剩余的电量缺额，则剩余的缺电量由配电系统补足，若所有缺电微电网的电量缺额都平衡后，微电网依然还有剩余的电量余额，则剩余的多电量全都输送给配电系统。