74文章解读与程序——《计及源荷不确定性的综合能源生产单元运行调度与容量配置两阶段随机优化程序》已提供下载资源

最新推荐文章于 2024-11-06 17:51:55 发布

电网论文源程序

最新推荐文章于 2024-11-06 17:51:55 发布

阅读量1.2k

点赞数 34

分类专栏：论文与完整源程序文章标签：能源

本文链接：https://blog.csdn.net/LIANG674027206/article/details/140678688

版权

论文与完整源程序专栏收录该内容

268 篇文章 277 订阅

订阅专栏

👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆

《《《《《《《《更多资源还请持续关注本专栏》》》》》》》

论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html

电网论文源程序-CSDN博客电网论文源程序擅长文章解读,论文与完整源程序,等方面的知识,电网论文源程序关注python,机器学习,计算机视觉,深度学习,神经网络,数据挖掘领域.https://blog.csdn.net/LIANG674027206?type=download

本程序采用全年光伏、风电数据通过kmeans聚类得到6种场景，构建了随机优化模型，在研究融合P2G与CCS的IEPU系统框架基础上，建立了各关键设备及生产环节数学模型，基于混合整数线性规划(mixed integer linear programming，MILP)算法，以全生命周期内经济成本最低为优化目标，考虑物料及能量平衡约束，实现典型周内各设备功率的最优逐时调度优化，并得到最佳综合能源系统中碳捕集+电制氢+甲烷化+氢存储+CO2存储的容量配置结果。程序和原文献的差别:原文献采用的是双层模型，将智能算法和MILP算法进行结合，但是考虑到智能算法由于容易陷入局优导致结果不稳定的问题，本程序采用单层MILP算法直接对模型进行复现;原文献算例采用的是典型周的方式，本文采用全年聚类形成6个典型日作为分析基础。

部分代码展示：

%% 3.2.1节有储气（碳-氢）模式程序如下：
%%  
clc
clear
close all
%%  初始化设备参数及运行变量
addpath('..\光照强度与电负荷生成');
load('IPV1');  load('IPV2'); load('IPV3'); load('IPV4'); load('IPV5'); load('IPV6');   
load('Eload1');  load('Eload2'); load('Eload3'); load('Eload4'); load('Eload5'); load('Eload6');
It = [IPV1,IPV2,IPV3,IPV4,IPV5,IPV6];
Edemand = [Eload1,Eload2,Eload3,Eload4,Eload5,Eload6]; %电负荷
load('Nday'); %各典型日频次   
nday = [Nday(1)*ones(1,24),Nday(2)*ones(1,24),Nday(3)*ones(1,24),Nday(4)*ones(1,24),Nday(5)*ones(1,24),Nday(6)*ones(1,24)];
%% 以下，注意是把六个典型日的约束一起写，维度是24*6
T = 24*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
E_CCSmax = 4.4551e+03; %kW

效果展示：