三微网优化matlab+yalmip 采用matlab+yalmip编程，实现三个微网的优化调度

三微网优化matlab+yalmip 采用matlab+yalmip编程，实现三个微网的优化调度，分别包括微网内燃气轮机、燃料电池、储能等主体约束，还包括微网间互供以及和电网间购售电约束，程序以成本最低作为目标，有对应的参考资料，出图共17张

一、系统概述

三微网优化调度系统基于MATLAB与YALMIP工具构建，以实现三个微网的低碳经济优化调度为核心目标。系统整合了微网内燃气轮机、燃料电池、光伏、风机等分布式电源，以及储能装置，同时支持微网间、微网与大电网的能量交互。通过建立多维度约束条件与成本优化模型，在保障能源供应稳定性的前提下，最小化系统运行总成本，兼顾经济收益与环保要求，适用于多微网互联场景下的高效能源管理。

二、核心功能模块

（一）数据导入与预处理模块

该模块负责读取并初步处理三微网系统运行所需的基础数据，为后续优化计算提供数据支撑，具体功能包括：

1. 可再生能源数据导入：读取三个微网各自的光伏出力数据（Ppv1、Ppv2、Ppv3）与风机出力数据（Pwt1、Pwt2、Pwt3），这些数据反映了不同时段可再生能源的发电能力，是微网能源供应的重要组成部分。
2. 负荷数据处理：导入三个微网的基础负荷数据（PL1、PL2、PL3），并根据实际需求将负荷值乘以1.5进行调整，得到最终的微网负荷需求，确保负荷数据与实际运行场景匹配。
3. 电价数据导入：获取三类关键电价数据，分别是微网与大电网间的购电电价（Cp1）、微网与微网间的分时电价（Cp2）、微网与大电网间的售电电价（Cp3），电价数据是成本计算与能量交互决策的核心依据。
4. 供需差额计算：通过遍历一天24个时刻，计算每个微网在各时刻的可再生能源出力与负荷的差额（pd1、pd2、pd3），以及三个微网整体的供需差额（pdz），快速判断各时刻微网的能源盈余或短缺状态，为调度策略制定提供初步参考。

（二）决策变量定义模块

该模块定义了三微网优化调度过程中所有需要求解的决策变量，涵盖能源生产、存储、交互等全环节，具体变量分类如下：

1. 分布式电源出力变量：包括三个微网各自的燃气轮机出力（PG1、PG2、PG3）与燃料电池出力（FC1、FC2、FC3），用于描述可控分布式电源的发电功率调节范围。
2. 储能运行变量：包含储能充电功率（Psch1、Psch2、Psch3）、储能放电功率（Psdis1、Psdis2、Psdis3）、储能容量（Es1、Es2、Es3），以及储能充放电状态二进制变量（Uch1、Udis1等），全面刻画储能系统的运行状态与能量变化。
3. 能量交互变量：分为微网与大电网交互变量（PMbuy1、PMsell1等，代表微网向大电网购电、售电功率）和微网间交互变量（Pwbuy12、Pwsell12等，代表不同微网间的购电、售电功率），同时定义了对应的能量交互状态二进制变量（Um1、Uw12等），用于约束能量交互的方向与启停状态。

（三）参数配置模块

该模块为三微网系统各设备与运行规则配置关键参数，确保优化模型符合实际设备特性与运行要求，主要参数类别包括：

1. 分布式电源参数：设定每个微网燃气轮机的出力上下限（如PG1min=80kW、PG1max=800kW）、爬坡功率限制（如Rd1=PG1min×0.5、Rup1=PG1max×0.5），以及燃料电池的出力上下限（如FC1min=50kW、FC1max=500kW），约束可控电源的出力调节速度与范围。
2. 储能系统参数：配置储能初始容量（如Es01=800kWh）、容量上下限（如Esmax1=Es01×0.9、Esmin1=Es01×0.2）、充放电功率上限（如Ps1max=200kW）及充放电效率（如eta1=0.95、eta11=0.97），同时设定储能放电损耗系数（u=0.02）与单位充放电成本系数（如Ks1=0.08），精准反映储能系统的运行特性与成本。
3. 能量交互参数：定义微网与大电网联络线功率上限（PMmax=1500kW）、微网间联络线功率上限（PWmax=500kW），限制能量交互的最大规模，保障系统稳定。
4. 成本与环保参数：设置分布式电源成本系数（如a1=0.7939为微网1燃气轮机成本系数、b1=0.6126为微网1燃料电池成本系数）、可再生能源维护成本系数（Kpv=0.08、Kwt=0.11），以及环保相关参数（如v1=0.04为燃气轮机单位电量CO₂排放系数、v3=0.031为单位CO₂治理费用），为成本与环保目标计算提供依据。

（四）约束条件构建模块

该模块构建了覆盖设备运行、能量平衡、交互规则等多维度的约束条件，确保优化结果的可行性与合理性，核心约束包括：

1. 分布式电源运行约束
   -出力约束：限制燃气轮机与燃料电池的出力在设定的上下限范围内，如PG1≥PG1min、PG1≤PG1max，FC1≥FC1min、FC1≤FC1max，避免设备超额定工况运行。
   -爬坡约束：约束燃气轮机相邻时刻的出力变化量不超过爬坡功率限制，如[PG1(2:24);PG1(1)]-PG1≤Rup1、PG1-[PG1(2:24);PG1(1)]≤Rd1，防止出力波动过大影响系统稳定。
2. 储能系统运行约束
   -充放电功率约束：根据储能充放电状态二进制变量，限制充电功率（如0≤Psch1≤Uch1×Ps1max）与放电功率（如0≤Psdis1≤Udis1×Ps1max），确保充放电功率不超过设备上限，且同一时刻储能不同时处于充放电状态（Uch1+Udis1≤1）。
   -容量约束：约束储能容量在设定的上下限范围内（如Es1≥Esmin1、Es1≤Esmax1），同时根据充放电功率、效率与损耗系数，建立储能容量随时间变化的动态模型，如t=1时Es1(1)=(1-u)×0.5×Esmax1+(eta1×Psch1(1)-Psdis1(1)/eta11)，t>1时Es1(t)=(1-u)×Es1(t-1)+(eta1×Psch1(t)-Psdis1(t)/eta11)。
3. 能量平衡约束
   -微网内功率平衡：对每个微网，在每个时刻建立功率平衡方程，确保能源供应与需求匹配。以微网1为例，PMbuy1+Pwbuy12+Pwbuy13-Pwsell12-Pwsell13-PMsell1-Psch1-PL1+PG1+Psdis1+Ppv1+Pwt1+FC1=0，涵盖了外购电、微网间购电、售电、储能充放电、分布式电源出力与负荷等所有能量流。
   -微网间功率平衡：约束微网间购售电功率相等，如Pwbuy12=Pwsell21、Pwbuy13=Pwsell31，确保能量在微网间的流动守恒，无功率损耗或盈余偏差。
4. 能量交互约束
   -微网与大电网交互约束：根据交互状态二进制变量，限制购电功率（如0≤PMbuy1≤(1-Um1)×PMmax）与售电功率（如0≤PMsell1≤Um1×PMmax），确保交互功率不超过联络线上限，且同一时刻不同时进行购电与售电。
   -微网间交互约束：类似微网与大电网交互约束，根据微网间交互状态二进制变量，限制微网间购售电功率（如0≤Pwbuy12≤(1-Uw12)×PWmax、0≤Pwsell12≤Uw12×PWmax），保障微网间交互的合理性。

（五）目标函数构建与求解模块

1. 目标函数构建：以三微网系统24小时总运行成本最小为目标，总成本涵盖分布式电源发电成本（燃气轮机、燃料电池）、储能运行成本、微网与大电网交互成本、微网间交互成本、可再生能源维护成本及环保成本。例如，微网1燃气轮机成本CG1=a1×PG1，微网1与大电网交互成本CM1=Cp1×PMbuy1-Cp3×PMsell1，环保成本C1=v3×(v1×(PG1+PG2+PG3)+v2×(PMbuy1+PMbuy2+PMbuy3))，最终目标函数为各成本项之和的累加。
2. 优化求解：配置优化求解器（选用Gurobi solver），通过调用YALMIP工具的求解函数，在构建的约束条件下求解目标函数，得到各决策变量的最优值。求解过程中会判断求解是否成功，若成功则输出总运行成本及各决策变量（如燃气轮机出力、储能充放电功率、能量交互功率等）的最优结果；若失败则提示“求解出错”。

（六）结果分析与可视化模块

该模块对优化求解结果进行处理，并通过图表形式直观展示，便于用户分析系统运行特性，主要功能包括：

1. 关键参数计算：根据求解得到的储能容量，计算储能荷电状态（SOC），如soc1=Es1/Es01，反映储能系统的能量存储水平。
2. 图表生成
   -储能SOC图表：以柱状图展示三个微网24小时内储能SOC变化，对比不同微网储能的运行状态。
   -微网功率流图表：以堆叠柱状图展示每个微网24小时内各类功率流（如购电、售电、分布式电源出力、储能充放电等），并叠加负荷曲线，清晰呈现微网内能量供需平衡情况。
   -电价图表：以柱状图展示24小时内配网购电电价、配网售电电价与微网间电价的变化，为电价影响分析提供依据。
   -可再生能源与负荷图表：以折线图展示每个微网24小时内光伏、风机出力与负荷的变化，分析可再生能源出力与负荷的匹配度。
   -设备运行状态图表：包括储能充放电状态柱状图、燃气轮机与燃料电池出力柱状图、微网间能量交互功率柱状图等，全面展示设备运行与能量交互的最优结果。

三、系统特色与优势

1. 多目标协同优化：系统同时考虑经济成本与环保要求，在最小化运行成本的同时，通过引入CO₂排放系数与治理费用，降低系统碳排放，实现低碳经济运行目标。
2. 灵活的能量交互机制：支持微网间点对点能量交互与微网与大电网的双向交互，优先通过微网间能量互补实现可再生能源就地消纳，减少对大电网的依赖，降低大电网运行压力。
3. 精准的设备建模：对燃气轮机、燃料电池、储能等设备进行精细化建模，考虑设备的出力限制、爬坡特性、充放电效率等实际参数，确保优化结果符合设备运行规律，具备实际可操作性。
4. 全面的约束覆盖：构建了涵盖设备运行、能量平衡、交互规则等多维度的约束体系，保障优化结果的可行性与系统运行的稳定性，避免出现设备过载、功率失衡等问题。
5. 直观的结果展示：通过丰富的图表展示优化结果，能够清晰呈现系统24小时内的运行特性，为用户分析调度策略的合理性、优化系统参数提供有力支撑。

四、应用场景与价值

该系统适用于多微网互联的园区、区域能源系统等场景，能够实现以下应用价值：

1. 提升能源利用效率：通过优化分布式电源出力与能量交互，最大化可再生能源消纳率，减少能源浪费，提升整体能源利用效率。
2. 降低运行成本：通过合理制定能量交互策略与设备运行计划，最小化系统运行成本，为微网运营方带来经济收益。
3. 保障系统稳定运行：通过严格的约束条件与精细化的设备建模，确保系统在优化运行过程中保持稳定，避免设备故障与功率失衡等问题。
4. 助力低碳转型：通过低碳目标融入优化模型，减少系统碳排放，符合当前低碳经济发展要求，为能源系统绿色转型提供技术支持。

五、使用建议

1. 数据准备：使用前需确保光伏、风机、负荷、电价等基础数据的准确性与完整性，数据格式需与代码导入要求一致（如Excel格式），避免因数据问题导致求解结果偏差。
2. 参数调整：根据实际设备型号与运行场景，调整分布式电源、储能、能量交互等相关参数（如出力上下限、效率、电价等），确保模型与实际系统匹配。
3. 求解器配置：确保MATLAB环境中已正确安装Gurobi求解器与YALMIP工具，且求解器许可证有效，避免因求解器问题导致求解失败。
4. 结果分析：结合生成的各类图表，深入分析系统运行特性，如可再生能源消纳率、储能SOC变化、能量交互成本等，根据分析结果进一步优化系统参数与调度策略。

三微网优化matlab+yalmip 采用matlab+yalmip编程，实现三个微网的优化调度，分别包括微网内燃气轮机、燃料电池、储能等主体约束，还包括微网间互供以及和电网间购售电约束，程序以成本最低作为目标，有对应的参考资料，出图共17张