👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html
论文可参考:
1、考虑气电联合需求响应的气电综合能源配网系统协调优化运行
2、Robust Network Hardening Strategy for Enhancing Resilience of Integrated Electricity and NaturalGas Distribution Systems Against Natural Disasters
在未来多能互补、综合能源系统的背景下,传统配电网和配气网独立调度运行的模式已经不能满足多种能源互补的运行要求。为此,该文提出气电综合能源配网系统最优潮流的凸优化方法,即利用二阶锥规划方法对配电网潮流方程约束进行处理,并提出运用增强二阶锥规划与泰勒级数展开相结合的方法对天然气潮流方程约束进行处理,进而将非线性的气电综合能源配网系统优化调度问题转化为混合整数二阶锥规划模型,为气电综合能源配网的气/电协调优化运行和规划设计提供支撑。
部分代码展示:
clc; clear ;
close all;
NT = 24;
eval('PDN_IEEE33');% 载入配电网数据
eval('GDN_Data'); % 载入配气网数据
%% PDN 配电网 数据 提取和处理
NL = size(PDN_data,1);
NB = NL + 1;
Li = PDN_data(:,2);
Lj = PDN_data(:,3);
R = PDN_data(:,6)/ZB;
X = PDN_data(:,7)/ZB;
% 负荷数据 可根据自己需要调增 数据可根据文献或老师的自行调整
Pd_level = [0;PDN_data(:,4)]/sum(PDN_data(:,4)); % 各节点有功、无功负荷占总额定负荷比
Qd_level = [0;PDN_data(:,5)]/sum(PDN_data(:,5));
% one day中各小时的总负荷功率 数据可根据文献或老师的自行调整
Pd_day = [1.5 1.5 1.5 1.5 1.5 1.6 1.9 2.1 2.1 2.0 2.0 2.1...
2.3 2.3 2.3 2.3 2.3 2.3 2.4 2.5 2.4 2.3 2.0 1.8 ]; % MW
Qd_day = [1.4 1.4 1.3 1.3 1.3 1.3 1.3 1.4 1.6 1.7 1.8 1.8 ...
1.8 1.8 1.8 1.7 1.7 1.7 1.5 1.6 1.5 1.5 1.6 1.6 ]; % MVar
Pd = Pd_level*Pd_day(1:24);
Qd = Qd_level*Qd_day(1:24);
Pd = Pd/SB;
Qd = Qd/SB; % p.u
U2_l = 0.95^2;
U2_u = 1.05^2; % p.u.
NT = 24;
NL = size(PDN_data,1);
NB = NL+1;
A0 = zeros(NB, NL); % 节点-支路关联矩阵
for i = 1:NL
A0(Li(i,1),i) = -1;
A0(Lj(i,1),i) = 1;
end
A = A0(2:end,1:end); % 降阶节点支路关联矩阵 P_injection = A * P_flow
M = inv(A); % P_flow = M * P_injection , 表征 线路潮流(首端) 与 所有非平衡节点注入 的关系
%% GDN 配气网 数据 提取和处理
% 生成配气网相关参数
n_Node = length(Gas_Node_Pressure(1,:)); % 节点数
n_Pipe = length(Gas_Pipe(:,1)); % 管道数
n_GW = length(Gas_Source(:,1)); % 气井,产气源
n_Comp = length(Gas_Compressor(:,1)); % 压缩机数量
n_GS = length(Gas_Storge(:,1)); % 燃气存储库
n_GT = length(Gas_GenUnit(:,1)); % 燃气轮机
% 管道参数
P_m = Gas_Pipe(:,2); % 管道首、末节点
P_n = Gas_Pipe(:,3);
K_Pipe = Gas_Pipe(:,4); % 管道参数/C_mn (kcf/bar)
% 节点参数
Ind_Node = Gas_Node_Pressure(1,:)';
Pre_min = Gas_Node_Pressure(2,:)';
Pre_max = Gas_Node_Pressure(3,:)';
% 气井参数 Gas Well
Ind_GW = Gas_Source(:,2);
G_GW_max = Gas_Source(:,3); % 气源最大出力/kcf
% 负荷参数
Ind_GLoad = Gas_Load(:,2); % 气负荷-节点集
GLoad_factor = Gas_Load(:,3); % 负荷分配因子
Gas_Load_24h = [ % 单位 kcf/h
46.784 41.7 39.666 40.174 43.226 47.802 47.294 51.362 52.378 52.888 51.87 51.87 ...
50.854 48.31 48.82 49.836 51.362 52.888 52.888 57.464 56.448 55.43 51.362 46.784];
Gas_Demand = zeros(n_Node,NT); % 所有节点 气负荷24h 矩阵
Gas_Demand(Ind_GLoad,:) = GLoad_factor * Gas_Load_24h; % 24小时各节点气负荷矩阵,kcf/h 暂取1小时
% 压缩机参数
Comp_m = Gas_Compressor(:,2);
Comp_n = Gas_Compressor(:,3);
gama_Cmn = Gas_Compressor(:,4); % 压缩因子
alpha_Cmn = Gas_Compressor(:,5); % 压缩机消耗燃气消耗系数
gama_Cmn = Gas_Compressor(:,4); % 压缩因子
alpha_Cmn = Gas_Compressor(:,5); % 压缩机消耗燃气消耗系数
% 燃气轮机参数
Ind_GT_PDN = Gas_GenUnit(:,2);
Ind_GT_GDN = Gas_GenUnit(:,3);
eff_GT = Gas_GenUnit(:,5); % 气转电效率
P_min_GT = Gas_GenUnit(:,6); % 产电有功值上下界
P_max_GT = Gas_GenUnit(:,7); %
Q_min_GT = Gas_GenUnit(:,8); % 无功功值上下界
Q_max_GT = Gas_GenUnit(:,9); %
% GDN 元素连接关系
A_GP = zeros(n_Node, n_Pipe); % 气网的节点-管道关联矩阵
for i = 1:n_Pipe
A_GP(P_m(i),i) = 1;
A_GP(P_n(i),i) = -1;
end
A_GC = zeros(n_Node, n_Comp); % 压缩机的节点-管道关联矩阵
for i = 1:n_Comp
A_GC(Comp_m(i),i) = 1;
A_GC(Comp_n(i),i) = -1;
end
A_GW = zeros(n_Node, n_GW); % 燃气供应GW 所处于节点的索引矩阵
for i = 1: n_GW
A_GW(Ind_GW(i),i) = 1;
end
A_GT_GDN = zeros(n_Node, n_GT); % 燃气轮机所处于GDN节点的索引矩阵
for i = 1:n_GT
A_GT_GDN(Ind_GT_GDN(i),i) = 1;
end
效果展示:
101号资源-源程序:论文可在知网下载《考虑气电联合需求响应的气电综合能源配网系统协调优化运行》本人博客有解读资源-CSDN文库https://download.csdn.net/download/LIANG674027206/89271347👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html