matlab|计及源荷不确定性的综合能源生产单元运行调度与容量配置随机优化模型

目录

[1 主要内容](#1 主要内容)

[1.1 风光场景聚类](#1.1 风光场景聚类)

[1.2 主模型程序结果](#1.2 主模型程序结果)

[1.3 随机模型和确定性模型对比](#1.3 随机模型和确定性模型对比)

[1.4 有无储气对比](#1.4 有无储气对比)

[1.5 煤价灵敏性分析](#1.5 煤价灵敏性分析)

[1.6 甲烷价格灵敏性分析](#1.6 甲烷价格灵敏性分析)

[2 部分程序](#2 部分程序)

[3 下载链接](#3 下载链接)


1 主要内容

本程序复现《计及源荷不确定性的综合能源生产单元运行调度与容量配置两阶段随机优化》模型,采用全年光伏、风电数据通过kmeans聚类得到6种场景,构建了随机优化模型,在研究融合P2G与CCS的IEPU系统框架基础上,建立了各关键设备及生产环节数学模型,基于混合整数线性规划(mixed integer linear programming, MILP)算法,以全生命周期内经济成本最低为优化目标,考虑物料及能量平衡约束,实现典型周内各设备功率的最优逐时调度优化,并得到最佳综合能源系统中碳捕集+电制氢+甲烷化+氢存储+CO2存储的容量配置结果。

  • 程序和原文献的差别:原文献采用的是双层模型,将智能算法和MILP算法进行结合,但是考虑到智能算法由于容易陷入局优导致结果不稳定的问题,本程序采用单层MILP算法直接对模型进行复现;原文献算例采用的是典型周的方式,本文采用全年聚类形成6个典型日作为分析基础。
  • 程序优势:本程序不仅复现了主模型,而且实现了文献的对比算例和敏感性分析部分,分别是确定性模型和随机模型对比、无储气与有储气对比、甲烷价格灵敏度分析、煤价灵敏度分析。
  • 程序采用matlab+yalmip(求解器为gurobi)进行优化,由于非线性问题的存在,该程序采用cplex求解会报错,需要学习的同学可以下载安装gurobi并配置好环境后运行使用。

1.1 风光场景聚类

通过全年风光数据聚类得到6个典型日并计算得到每个典型日的频次。

1.2 主模型程序结果

1.3 随机模型和确定性模型对比

原文结果:因原始数据偏差,结果不尽相同,但趋势是一致的。

1.4 有无储气对比

1.5 煤价灵敏性分析

1.6 甲烷价格灵敏性分析

2 部分程序

复制代码
    %% 1.1.1光伏设备模型 
    E_PVmppt = sdpvar(1,T);  %光伏板mppt发电功率
    A_PV = sdpvar(1,1);  %光伏板面积/m2
    k = 0.200;  %1平方米的光伏板1000w/m2的标准电功率为200w
    E_PVr = sdpvar(1,1);  %光伏板额定发电功率
    ita_PV = 0.200/1000;
    %文章内写了两个E_PV,有错位,本代码将其改为E_PVmppt与E_PV
    E_PV = sdpvar(1,T);  %光伏板有效发电功率
    E_PV_cur = sdpvar(1,T);  %弃光功率
    %后文算例中出现135MW的光伏容量配置结果,那么这里的限值就算用300MW吧,即300 000kW.
    E_PVr_max = 300000;  %光伏板额定发电功率.kW
    %之后,这里直接将约束也写上,省的再回头来写约束了。
    C=[];
    C=[C, E_PVr ==  A_PV*k,
          E_PVmppt == E_PVr*ita_PV/k*It,
          E_PVmppt == E_PV + E_PV_cur,
          0<=E_PVr,E_PVr<=E_PVr_max,
          %补充
          E_PV >= 0,
          E_PV_cur >= 0,
          A_PV >= 0, 
    ];
    %% 1.1.2 CCS 模型   
    V_CO2_PGU = sdpvar(1,T);  %火电机组的二氧化碳排放量
    E_PGU = sdpvar(1,T);  %火电机组发电功率
    e_PGU = 0.46;   %火电机组的二氧化碳排放强度,见表1的 0.46 N.m3CO2/kW.h
    ita_CCS_max = 0.65;%碳捕集效率最大值  0.65 
    V_CO2_CCSmax = sdpvar(1,T); %碳捕集最大功率(体积)
    V_CO2_CCS = sdpvar(1,T); %实际碳捕集功率(体积)
    V_CO2_cur = sdpvar(1,T); %碳捕集功率耗散部分功率(体积)
    lamdaCO2 = 0.1937;   %碳捕集功率耗电系数 kW.h/N.m3CO2
    E_CCS = sdpvar(1,T); %碳捕集耗电功率
    %从图5可以找出火电机组的最大出力功率180MW,最大爬坡常出现在119时刻与162时刻的正负50MW
    %火电机组的最小出力功率90MW,
    %表1中给出的火电机组容量为300000kW
    E_PGUmax = 300000;  %kW
    E_PGUmin = 90000;   %kW
    dita_E_PGUmax = 50000;%kW
    dita_E_PGUmin = -50000;%kW
    %从图7可知CO2捕集的最大功率是23000m3每小时
    %由此计算碳捕集的最大电功率为 0.1937*23000 = 4.4551e+03 kW

3 下载链接

相关推荐
程高兴1 小时前
基于Matlab的车牌识别系统
开发语言·matlab
XuX034 小时前
手搓雷达图(MATLAB)
matlab·贴图
freexyn8 小时前
Matlab自学笔记五十一:(推荐)输入参数的数量和可变数量的输入
笔记·算法·matlab
不吃酸的柠檬8 小时前
MATLAB 中的图形绘制
人工智能·机器学习·matlab
studyer_domi8 小时前
Matlab 复合模糊PID
开发语言·matlab
DarrenPig1 天前
【新能源科学与技术】MATALB/Simulink小白教程(一)实验文档【新能源电力转换与控制仿真】
matlab·开源·github·simulink·交流
简简单单做算法1 天前
基于GA遗传优化TCN-BiGRU注意力机制网络模型的时间序列预测算法matlab仿真
matlab·tcn-bigru·时间序列预测·注意力机制·ga遗传优化
Cc小跟班2 天前
MATLAB小技巧记录(特殊符号、图例位置...)
开发语言·算法·matlab
随风飘摇的土木狗2 天前
【MATLAB第116期】基于MATLAB的NBRO-XGBoost的SHAP可解释回归模型(敏感性分析方法)
matlab·gsa·敏感性分析·特征排序·灵敏度分析·shap·可解释
海天一色y2 天前
matlab设置不同颜色的柱状图
matlab·信息可视化