技术领域
[0001] 本
发明属于电
力系统技术领域，具体涉及一种冷热电联供微电网系统的双层 优化方法。
背景技术
[0002] 工业技术的快速发展，对
能源供应提出越来越高的要求。传统的石油、
煤炭 等化石
燃料，总量虽然丰富，但是不可再生。而在我国，化石能源虽然总量大， 但人均占有量远低于世界平均
水平，因此提高能源利用效率成为缓解能源危机的 重要方法。冷热电联供系统的供能形式多样，能够通过
能量的
梯级利用，使一次 能源利用率高达90％，成为国内外关注的重点。
[0003] 但是现有的冷热电联供型多微网主动配电系统其能源利用率的提升收到了
瓶颈，所以需要一个新的技术方案来解决这个问题。
发明内容
[0004] 本发明所要解决的技术问题在于：针对背景技术的
缺陷，本发明提供一种冷 热电联供微电网系统的双层优化方法，可得出冷热电联供型多微网和主动配电网 的双层优化结果，有效提高能源利用效率，降低系统运行成本。
[0005] 本发明了解决以上技术问题采用以下技术方案：
[0006] 本发明提供一种冷热电联供微电网系统的双层优化方法，包括以下步骤：
[0007] S1：对冷热电联供型多微网内设备进行数学建模；
[0008] S2：建立双层优化模型的上层主动配电网层优化模型；
[0009] S3：建立双层优化模型的下层冷热电联供型多微网层优化模型；
[0010] S4：采用遗传
算法和混合整数线性规划进行求解，得出上层模型和下层模型 的日前优化调度计划。
[0011] 进一步的，本发明提出的一种冷热电联供微电网系统的双层优化方法，所述 步骤S1中的冷热电联供型多微网内设备包括微型
燃气轮机、燃气
锅炉、余热锅 炉、
吸收式制冷机组、
蒸汽换热装置、电制冷机、储能装置及
可再生能源发电装 置，所述数学建模的具体步骤如下：
[0012] S1-1：建立燃气轮机的数学模型：
[0013] ηc＝(8.935+33.157β-27.081β2+17.989β3)/100×100％
[0014] ηr＝(82.869-30.173β+24.644β2-16.371β3)/100×100％
[0015]
[0016]
[0017] 其中，β为机组电负荷率；ηc为燃气轮机发电效率，ηr为燃气轮机热回收效 率，QGT为燃气轮机排气余热量，PGT为燃气轮机发电功率，ηl为燃气轮机
散热损 失系数，VGT为运行时间内燃气轮机所消耗的
天然气量，LHVNG为天然气热值；t 为运行时间；
[0018] S1-2：建立储能装置的数学模型：
[0019]
[0020] 其中，E(t)为储能装置在t时段储存的能量，Δt为t时段到t+1时段的时间 间隔，Pabs(t)为t时段储能功率，Prelea(t)为t时段放能功率，μ为储能装置自身向 环境散能损失或自损耗的能量系数，ηabs为储能装置的蓄能效率，ηrelea为储能 装置放能效率；
[0021] S1-3：建立燃气锅炉的数学模型：
[0022]
[0023] 其中，PGB为燃气锅炉的热功率， 为燃气锅炉在Δt时间段消耗的燃气量，ηGB为燃气锅炉的热效率；LNG为燃气热值；
[0024] S1-4：建立换热装置的数学模型：
[0025] PHX,out＝PWH,heatηHX
[0026] 其中，PHX,out为蒸汽热水换热装置输出制热功率，PWH,heat为
余热锅炉的输出 蒸汽中用于制热的热功率，ηHX为蒸汽热水换热装置的转换效率；
[0027] S1-5：建立吸收式制冷机的数学模型：
[0028] PAC,out＝PWH,coolηAC
[0029] 其中，PAC,out为蒸汽型吸收式制冷机输出制冷功率，PWH,cool为余热锅炉的输出 蒸汽中用于制冷的输入功率，ηAC为蒸汽型吸收式制冷机的制冷效率；
[0030] S1-6：建立电制冷机的数学模型：
[0031] PEC,out＝PEC,inηEC
[0032] 其中，PEC,out为电制冷机的输出制冷功率，PEC,in为电制冷机的输入电功率，ηEC为电制冷机的能效比。
[0033] 进一步的，本发明提出的一种冷热电联供微电网系统的双层优化方法，所述 步骤S2中建立上层主动配电网层的优化模型包括以下步骤：
[0034] S2-1：建立上层优化目标函数，即主动配电网层的日运行经济成本最小：
[0035] minFDN＝FG-Fs
[0036]
[0037]
[0038] 其中，FDN为配电网的总生产成本，FG为配电网中发
电机组的发电成本，Fs为主动配电网与各个混合能源系统之间交互电功率的
费用，PGi(t)为第i台常规 机组在时段t的出力，ai、bi、ci为对应的成本系数，n表示
发电机组的数量， 为时段t主动配电网与冷热电联供型微网j交互的电功率，正值表示电网向 多微网售电，负值表示电网从多微网购电，τj(t)为t时刻电网与微网j的实时交 易电价，m为冷热电联供型微网个数；
[0039] S2-2：建立主动配电网的约束条件，具体包括电功率平衡约束、常规机组出 力上下限约束、主动配电网旋转备用约束、联络线传输功率约束：
[0040] 电功率平衡约束：
[0041]
[0042] 常规机组出力上下限约束：
[0043] PGi,min≤PGi(t)≤PGi,max
[0044] 配电网旋转备用约束：
[0045]
[0046] 联络线传输功率约束：
[0047]
[0048] 其中， 为主动配电网在t时段的电负荷预测值，PGi,max和PGi,min为机组i 的有功出力上下限，RDN(t)为主动配电网在t时段的备用需求， 和 为主动配 电网向各个冷热电联供型微网传输功率的上下限。
[0049] 进一步的，本发明提出的一种冷热电联供微电网系统的双层优化方法，所述 步骤S3中建立下层冷热电联供型多微网层的优化模型包括以下步骤：
[0050] S3-1：冷热电联供型多微网经济优化调度模型中的优化目标函数：
[0051]
[0052] 其中，j为冷热电联供型微网的编号，m为冷热电联供型微网的个数，Ffuel,j为 第j个冷热电联供型微网的燃料费用，Fdisnet,j为第j个冷热电联供型微网与主动配 电网交互的功率费用；
[0053] S3-2：建立冷热电联供型多微网经济优化调度模型的约束条件，包括功率平 衡约束、设备容量约束、设备运行约束，其中，设备容量及运行约束条件是满足 各个设备的出力功率上下限限制，所述功率平衡约束具体如下：
[0054] 建立冷功率平衡约束方程：
[0055]
[0056] 其中， 为电制冷机的制冷功率， 为吸收式制冷机制冷功率，其来源 为余热锅炉收集的燃气轮机废热量， 为冷热电联供型微网中冷负荷的需求量；
[0057] 建立热功率平衡约束方程：
[0058]
[0059] 其中， 为燃气锅炉的输出热功率， 为蒸汽热水换热装置的制热功率， 为冷热电联供型微网中热负荷的需求量；
[0060] 建立电功率平衡约束方程：
[0061]
[0062] 其中， 为燃气轮机的发电功率， 为冷热电联供型微网和主动配电网 的逐时电功率交换值， 为
风力发电功率， 为
光伏发电功率， 为冷热电 联供型微网电负荷功率， 为冷热电联供型微网中电制冷机耗电功率； 为
蓄电池放电功率； 为
蓄电池充电功率；
[0063] 建立设备出力功率上下限约束：
[0064]
[0065] 其中， 和 为燃气轮机出力功率的最小值和最大值， 和 为燃 气锅炉出力功率的最小值和最大值， 和 为换热装置出力功率的最小值和 最大值，和 为吸收式制冷机出力功率的最小值和最大值， 和 为 电制冷机出力功率的最小值和最大值， 和 为蓄电池充电功率的最小值和 最大值， 和 为蓄电池放电功率的最小值和最大值。
[0066] 进一步的，本发明提出的一种冷热电联供微电网系统的双层优化方法，所述 步骤S4具体为：
[0067] 采用
遗传算法与混合整数线性规划
软件相结合的方法进行求解，上层采用遗 传算法，对主动配电网的发电机组出力和与各个冷热电联供型微网间购售电功率 进行寻优，下层调用混合整数线性规划软件对调度周期内各个冷热电联供型微网 的设备出力进行计算，并返回各个冷热电联供型微网与主动配电网购售电的功率 值，上层计算适应度值，通过遗传操作寻优，得到设定目标下的最优解。
[0068] 进一步的，本发明提出的一种冷热电联供微电网系统的双层优化方法，所述 步骤S4中对于双层优化模型的求解算法流程步骤具体如下：
[0069] S4-1：读取基本数据；
[0070] S4-2：随机生成初始种群数据；
[0071] S4-3：调用混合整数线性规划软件求解下层模型冷热电联供型多微网系统的 优化模型；
[0072] S4-4：计算上层模型的个体适应度值；
[0073] S4-5：判断
迭代次数是否大于30；
[0074] S4-6：如果判断结果为是，则输出上层和下层的最优解；
[0075] S4-7：如果判断结果为否，则重新形成遗传种群，重复步骤S4-2，直到判断 结果为是。
[0076] 本发明采用以上技术方案，与
现有技术相比的有益效果：
[0077] 本发明与现有技术相比，通过建立双层优化模型，分析主动配电系统和冷热 电联供型多微网系统的经济成本和设备出力等情况，不但具有重要的工程指导意 义，而且能够有效的提高冷热电联供型多微网主动配电系统的能源利用效率，降 低系统运行成本。
附图说明
[0079] 图2为本发明冷热电联供型多微网主动配电系统结构图。
[0080] 图3为本发明双层优化模型的求解
流程图。
具体实施方式
[0081] 为了使本发明的目的、技术方案及优点更加清楚明白，以下结合
实施例和说 明书附图对本发明的技术方案进行清楚、完整的描述，显然，此处所描述的具体 实施例仅用以解释本发明，并不用于限定本发明。
[0082] 本技术领域技术人员可以理解的是，除非另外定义，这里使用的所有术语(包 括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相 同的意义。还应该理解的是，诸如通用字典中定义的那些术语应该被理解为具有 与现有技术的上下文中的意义一致的意义，并且除非像这里一样定义，不会用理 想化或过于正式的含义来解释。
[0083] 如图1和图2所示，本发明提供一种冷热电联供微电网系统的双层优化方法， 包括以下步骤：
[0084] S1：对冷热电联供型微网的设备进行数学建模；
[0085] S2：建立双层优化模型的上层主动配电网层优化模型；
[0086] S3：建立双层优化模型的下层冷热电联供型多微网层优化模型；
[0087] S4：采用遗传算法和混合整数线性规划进行求解，得出上层模型和下层模型 的日前优化调度计划。
[0088] 本实施例步骤S1中的冷热电联供型多微网内设备包括微型燃气轮机、燃气 锅炉、余热锅炉、吸收式制冷机组、蒸汽换热装置、电制冷机、储能装置及可再 生能源发电装置，数学建模的具体步骤如下：
[0089] S1-1：建立燃气轮机的数学模型：
[0090] ηc＝(8.935+33.157β-27.081β2+17.989β3)/100×100％
[0091] ηr＝(82.869-30.173β+24.644β2-16.371β3)/100×100％
[0092]
[0093]
[0094] 其中，ηc为燃气轮机发电效率，ηr为燃气轮机热回收效率，QGT为燃气轮机 排气余热量，PGT为燃气轮机发电功率，ηl为燃气轮机散
热损失系数，VGT为运行 时间内燃气轮机所消耗的天然气量，LHVNG为天然气热值；
[0095] S1-2：建立储能装置的数学模型：
[0096]
[0097] 其中，E(t)为储能装置在t时段储存的能量，Δt为t时段到t+1时段的时间 间隔，Pabs(t)为t时段储能功率，Prelea(t)为t时段放能功率，μ为储能装置自身向 环境散能损失或自损耗的能量系数，ηabs为储能装置的蓄能效率，ηrelea为储能 装置放能效率；
[0098] S1-3：建立燃气锅炉的数学模型：
[0099]
[0100] 其中，PGB为燃气锅炉的热功率， 为燃气锅炉在Δt时间段消耗的燃气量，ηGB为燃气锅炉的热效率；
[0101] S1-4：建立换热装置的数学模型：
[0102] PHX,out＝PWH,heatηHX
[0103] 其中，PHX,out为蒸汽热水换热装置输出制热功率，PWH,heat为余热锅炉的输出 蒸汽中用于制热的热功率，ηHX为蒸汽热水换热装置的转换效率；
[0104] S1-5：建立吸收式制冷机的数学模型：
[0105] PAC,out＝PWH,coolηAC
[0106] 其中，PAC,out为蒸汽型吸收式制冷机输出制冷功率，PWH,cool为余热锅炉的输出 蒸汽中用于制冷的输入功率，ηAC为蒸汽型吸收式制冷机的制冷效率；
[0107] S1-6：建立电制冷机的数学模型：
[0108] PEC,out＝PEC,inηEC
[0109] 其中，PEC,out为电制冷机的输出制冷功率，PEC,in为电制冷机的输入电功率，ηEC为电制冷机的能效比。
[0110] 本实施例步骤S2中建立上层主动配电网层的优化模型包括以下步骤：
[0111] S2-1：建立上层优化目标函数，即主动配电网层的日运行经济成本最小：
[0112] minFDN＝FG-Fs
[0113]
[0114]
[0115] 其中，FDN为配电网的总生产成本，FG为配电网中发电机组的发电成本，Fs为主动配电网与各个混合能源系统之间交互电功率的费用，PGi(t)为第i台常规 机组在时段t的出力，ai、bi、ci为对应的成本系数，n表示发电机组的数量， 为时段t主动配电网与冷热电联供型微网j交互的电功率，正值表示电网向 多微网售电，负值表示电网从多微网购电，τj(t)为t时刻电网与微网j的实时交 易电价，m为冷热电联供型微网个数；
[0116] S2-2：建立主动配电网的约束条件，具体包括电功率平衡约束、常规机组出 力上下限约束、主动配电网旋转备用约束、联络线传输功率约束：
[0117] 电功率平衡约束：
[0118]
[0119] 常规机组出力上下限约束：
[0120] PGi,min≤PGi(t)≤PGi,max
[0121] 配电网旋转备用约束：
[0122]
[0123] 联络线传输功率约束：
[0124]
[0125] 其中， 为主动配电网在t时段的电负荷预测值，PGi,max和PGi,min为机组i 的有DN功出力上下限，R (t)为主动配电网在t时段的备用需求， 和 为主动配 电网向各个冷热电联供型微网传输功率的上下限。
[0126] 本实施例步骤S3中建立下层冷热电联供型多微网层的优化模型包括以下步 骤：
[0127] S3-1：冷热电联供型多微网经济优化调度模型中的优化目标函数：
[0128]
[0129] 其中，j为冷热电联供型微网的编号，m为冷热电联供型微网的个数，Ffuel,j为 第j个冷热电联供型微网的燃料费用，Fdisnet,j为第j个冷热电联供型微网与主动配 电网交互的功率费用；
[0130] S3-2：建立冷热电联供型多微网经济优化调度模型的约束条件，包括功率平 衡约束、设备容量约束、设备运行约束，其中，设备容量及运行约束条件是满足 各个设备的出力功率上下限限制，所述功率平衡约束具体如下：
[0131] 建立冷功率平衡约束方程：
[0132]
[0133] 其中， 为电制冷机的制冷功率， 为吸收式制冷机制冷功率，其来源 为余热锅炉收集的燃气轮机废热量， 为冷热电联供型微网中冷负荷的需求量；
[0134] 建立热功率平衡约束方程：
[0135]
[0136] 其中， 为燃气锅炉的输出热功率， 为蒸汽热水换热装置的制热功率， 为冷热电联供型微网中热负荷的需求量；
[0137] 建立电功率平衡约束方程：
[0138]
[0139] 其中， 为燃气轮机的发电功率， 为冷热电联供型微网和主动配电网 的逐时电功率交换值， 为
风力发电功率， 为光伏发电功率， 为冷热电 联供型微网电负荷功率， 为冷热电联供型微网中电制冷机耗电功率；
[0140] 建立设备出力功率上下限约束：
[0141]
[0142] 其中， 和 为燃气轮机出力功率的最小值和最大值， 和 为燃 气锅炉出力功率的最小值和最大值， 和 为换热装置出力功率的最小值和 最大值，和 为吸收式制冷机出力功率的最小值和最大值， 和 为 电制冷机出力功率的最小值和最大值， 和 为蓄电池充电功率的最小值和 最大值， 和 为蓄电池放电功率的最小值和最大值。
[0143] 本实施例步骤S4具体为：
[0144] 采用遗传算法(GA)与混合整数线性规划软件(Cplex)相结合的方法进行 求解，上层采用GA算法，对主动配电网的发电机组出力和与各个冷热电联供型 微网间购售电功率进行寻优，下层调用Cplex对调度周期内各个冷热电联供型微 网的设备出力进行计算，并返回各个冷热电联供型微网与主动配电网购售电的功 率值，上层计算适应度值，通过遗传操作寻优，得到设定目标下的最优解。
[0145] 如图3所示，步骤S4中对于双层优化模型的求解算法流程步骤具体如下：
[0146] S4-1：读取基本的设备参数、调度周期、成本系数等基本数据；
[0147] S4-2：随机生成初始种群数据；
[0148] S4-3：调用混合整数线性规划软件求解下层模型冷热电联供型多微网系统的 优化模型；
[0149] S4-4：计算上层模型的个体适应度值；
[0150] S4-5：判断迭代次数是否大于30；
[0151] S4-6：如果判断结果为是，则输出上层和下层的最优解；
[0152] S4-7：如果判断结果为否，则重新形成遗传种群，重复步骤S4-2，直到判断 结果为是；
[0153] S4-8：结束。
[0154] 本实施例采用天津中新生态城的典型冷热电联供型多微网系统算例，分别计 算两种传统运行方式下的多微网系统总成本，并与本发明提出的双层优化调度方 法进行对比，结果如表1所示：
[0155]运行方式 多微网系统总运行成本(元)
以热定电 39265.3
以电定热 40894.2
双层优化 37541.9
[0156] 表1不同运行方式下多微网系统总运行成本
[0157] 与传统“以热定电”和“以电定热”运行方式相比，本发明提出的双层优化 调度方法可分别降低5.33％和9.58％的总运行成本。通过本发明提出的优化调 度方法，调度冷热电联供型多微网系统中各微网设备出力和各微网与电网、相邻 微网间交互电功率值，可以显著降低多微网系统的总运行成本。
[0158] 本技术领域技术人员可以理解的是，本发明中已经讨论过的各种操作、方法、 流程中的步骤、措施、方案可以被交替、更改、组合或删除。进一步地，具有本 发明中已经讨论过的各种操作、方法、流程中的其他步骤、措施、方案也可以被 交替、更改、重排、分解、组合或删除。进一步地，现有技术中的具有与本发明 中公开的各种操作、方法、流程中的步骤、措施、方案也可以被交替、更改、重 排、分解、组合或删除。
[0159] 上面结合附图对本发明的实施方式作了详细地说明，但是本发明并不局限于 上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本 发明宗旨的前提下做出各种变化。
