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


    目录

    1 主要内容

    1.1 风光场景聚类

    1.2 主模型程序结果

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

    1.4 有无储气对比

    1.5 煤价灵敏性分析

    1.6 甲烷价格灵敏性分析

    2 部分程序

    3 下载链接


    主要内容

    本程序复现《计及源荷不确定性的综合能源生产单元运行调度与容量配置两阶段随机优化》模型,采用全年光伏、风电数据通过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 甲烷价格灵敏性分析

    部分程序

        %% 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 下载链接

  • 相关阅读:
    【精选】自学网络安全的三个必经阶段(含路线图)
    微信反复读取你的相册!官方回应:为便于聊天时快速发图,最新版本将取消...
    Android OpenGL 仿自如 APP 裸眼 3D 效果
    java计算机毕业设计高校网上报销系统MyBatis+系统+LW文档+源码+调试部署
    智慧公厕技术解析 | 提高公厕管理效率的关键
    面试金典--面试题 17.21. 直方图的水量(不困难的困难题)
    【毕业设计源码】基于SSM的小程序任务调度管理信息系统设计与实现
    【ConcurrentHashMap】JDK1.7版本源码解读与分析
    ElasticSearch-查询语法(聚合查询)
    状态估计|基于 MMSE 的分析估计器的不确定电力系统分析(Matlab代码实现)
  • 原文地址:https://blog.csdn.net/superone89/article/details/136280331