【全部更新完毕】2024电工杯A题数学建模详细思路代码文章分享

本文主要是介绍【全部更新完毕】2024电工杯A题数学建模详细思路代码文章分享,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

A 题:园区微电网风光储协调优化配置
摘要
在全球范围内,气候变化和环境污染问题日益严重,减少碳排放和实现可持续发展成为各国的共同目标。新能源,尤其是风能和光伏发电,因其清洁、可再生的特性,正在全球范围内得到广泛应用。然而,风能和光伏发电的间歇性和波动性导致与用电负荷的时序不匹配,产生了弃风弃光问题。为解决这一问题,储能技术的应用显得尤为重要。通过储能设备,可以存储多余的风光电能,在用电高峰时段释放,提高风光发电的利用率,减少对主电网的依赖。

通过合理配置储能系统和优化运行策略,三个园区在可再生能源利用和经济性方面都有显著提升。特别是园区B,通过有效的储能系统配置和管理,实现了最低的单位电量平均供电成本。园区A尽管在经济性上表现不如园区B,但通过进一步优化,仍有较大提升空间。此次研究为实现经济高效的园区微电网系统提供了参考和指导。

在分析联合园区储能配置及其经济性时,我们通过线性规划优化方法,比较了联合园区在未配置储能和最优配置储能两种情况下的经济性。联合园区的购电量和总供电成本均显著低于独立园区的总和,单位电量平均供电成本也较低。同时,联合园区的电池功率上限和容量明显减少,表明通过共享储能资源,可以有效降低投资成本并提高能源利用效率。通过联合运行和优化储能配置,可以实现更高效的能源利用和显著的成本节约,具有重要的实际应用价值和推广前景。

随着园区电负荷需求的增加,如何在满足更高负荷的同时保持经济性和可持续性成为关键问题。本文通过模拟分析了在电负荷增加50%的条件下,各园区分别在独立运营和联合运营模式下的风光储协调配置方案及其经济性。结果表明在电负荷增加50%的情况下,联合运营模式展现了更高的经济性和资源利用效率。

关键词 :风能发电,光伏发电,混合整数线性规划 (MILP),Gurobi,磷酸铁锂电池

目录
A 题:园区微电网风光储协调优化配置 1
摘要 1
一、背景 4
二、问题重述与分析 6
三、模型假设与参数说明 8
四、问题求解 10
4.1问题一 各园区独立运营储能配置方案及其经济性分析 11
4.1.1分析未配置储能时各园区运行的经济性 11
4.1.2各园区分别配置50kW/100kWh储能,制定储能最优运行策略及购电计划 14
4.1.3是否最优及最优储能配置方案 19
4.2问题二 联合园区储能配置方案及其经济性分析 24
4.2.1若未配置储能,分析联合园区运行经济性 24
4.2.2制定联合园区的总储能最优配置方案 28
4.2.3与各园区独立运营相比,园区联合运营的经济收益 31
4.3问题三 园区风、光、储能的协调配置方案及其经济性分析 34
4.3.1分别按各园区独立运营、联合运营制定风光储协调配置方案 34
4.3.2按全年12个月典型日风光发电功率数据,制定各园区独立运营的风光储协调配置方案 47
五、展望 55
六、 参考文献 56
附录【自行删减】 57


### https://docs.qq.com/doc/DVWRIQUlKaVNqcWFr

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

四、问题求解
求解工具选择Gurobi搭配Matlab
Gurobi Optimizer 是一个用于数学优化(数值优化)问题求解的商业软件,它能够高效地求解线性规划(LP)、二次规划(QP)、二次约束规划(QCP)、整数规划(IP)、混合整数规划(MIP)、非线性规划(NLP)等问题。Gurobi 被广泛应用于工业、金融、能源、物流等多个领域,是优化求解的首选工具之一。

4.1问题一 各园区独立运营储能配置方案及其经济性分析
4.1.1分析未配置储能时各园区运行的经济性
对于ABC三个园区,都采用线性规划的方式来解决电网的优化调配问题优化方法:线性规划问题
优化变量:,
目标函数:
限制变量:
电平衡:
其他:
求解结果:
表1给出了三个园区未配置储能时各园区运行的经济性的购电量,弃风弃光电量,总供电成本和单位电量平均供电成本信息,可以看出:
1)园区A的购电量和供电成本最高,说明其对主电网的依赖程度较高,自发电能力不足。相反,园区B购电量和供电成本最低,表明其自发电能力较强,对主电网依赖较小。
2)园区C的弃风弃光电量最高,说明其可再生能源发电能力较强,但由于负荷时序不匹配,未能充分利用这些能源。减少弃风弃光电量是提高可再生能源利用率的关键。
3)园区A的单位电量平均供电成本最高,这不仅是因为购电量大,还因为可再生能源利用率低,导致整体供电成本较高。园区B的单位电量平均供电成本最低,说明其可再生能源利用率较高,经济性较好。

园区 A B C
购电量/kWh 4874.125 2432.3 2699.39
弃风弃光电量/kWh 951.2 897.5 1128.02
总供电成本/¥ 6863.1625 4902.38 5510.87
单位电量平均供电成本/(¥/kWh) 0.8686448 0.63585 0.7087
表1 三个园区未配置储能时各园区运行的经济性

在这里插入图片描述

图4 园区A未配置储能各时段电力调配信息

4.1.2各园区分别配置50kW/100kWh储能,制定储能最优运行策略及购电计划
优化方法:混合整数线性规划问题
优化变量:,,, ,
目标函数:

限制变量:
电平衡:
电池状态:
电池容量变换范围限制:
电池充放电状态限制:
不能同时充放电:
其他:

在这里插入图片描述

求解结果:
表2给出了三个园区未配置储能时各园区运行的经济性的购电量,弃风弃光电量,总供电成本和单位电量平均供电成本信息,可以看出:
1)园区A的购电量和供电成本最高,说明其对主电网的依赖程度较高,自发电能力不足。相反,园区B购电量和供电成本最低,表明其自发电能力较强,对主电网依赖较小。

4.1.3是否最优及最优储能配置方案
显然50kW/100kWh的 方案不是最优的,因为可以明显看出(2)的各个结果中仍然出现了大量的弃电现象,这会导致大量的浪费我们需要增加电池的电量来获取更好的电量‘去峰填谷’效果
优化方法:混合整数线性规划问题
优化变量:,,, ,,,
目标函数:

限制变量:
电平衡:
电池状态:
电池容量变换范围限制:
电池充放电状态限制:
不能同时充放电:
其他:
优化结果:
表3给出了三个园区配置最优储能时各园区运行的经济性的购电量,弃风

4.2.3与各园区独立运营相比,园区联合运营的经济收益
经济收益比较:
计算各园区独立运营和联合运营情况下的总购电量、总弃风弃光电量、总供电成本和单位电量平均供电成本。
园区 A B C 总和 联合园区
购电量/kWh 4015.667 1622.30625 1681.35195 7319.325 7149.719563
弃风弃光电量/kWh 0 0 0 0 0
总供电成本/¥ 6604.40039 4447.857303 5176.151573 16228.41 14480.59104
单位电量平均供电成本/(¥/kWh) 0.83589424 0.576894592 0.665657353 0.739382 0.619172662
电池功率上限/kW 194.625 212.9 260.69 668.215 448.115
电池容量/kWh 1129.55 626.1940789 1269.758125 3025.502 1469.145313
表7 独立园区与联合园区最优配置储能时运行经济性比较

在这里插入图片描述

表7展示了独立园区与联合园区最优配置储能时运行经济性比较,从表中可以看出

  1. 购电量
    独立园区:A、B、C三个园区的购电量分别为4015.667 kWh、1622.30625 kWh和1681.35195 kWh,总和为7319.325 kWh。
    联合园区:购电量为7149.719563 kWh,比独立园区的总和减少169.605437 kWh。
    分析:联合园区的购电量减少表明,通过联合运行和共享储能资源,可以有效降低购电需求,提高能源利用效率。
  2. 弃风弃光电量:
    所有园区的弃风弃光电量均为0。

4.3问题三 园区风、光、储能的协调配置方案及其经济性分析
4.3.1分别按各园区独立运营、联合运营制定风光储协调配置方案
条件:各园区的最大负荷增长50%,且负荷波动特性不变。风电、光伏电源的配置成本分别为3000元/kW、2500元/kW。投资回报期为5年。
分时电价如表8所示(需要提供表8的具体分时电价)。
时段 电价(元/kWh)
7:00-22:00 1
其余时段 0.4
表8 分时电价表
计算负荷增长后的最大负荷:
P_load.A_new = PLmax.A * 1.5
P_load.B_new = PLmax.B * 1.5
P_load.C_new = PLmax.C * 1.5
独立运营的风光储协调配置方案:
确定各园区风电和光伏的装机容量,使得满足未来的负荷需求。计算储能的配置容量和功率,使得能在高峰负荷时段提供足够的电力,并在风光发电高峰时段存储电能。
联合运营的风光储协调配置方案:
计算联合园区的总负荷需求和总风光装机容量。计算联合园区的储能配置容量和功率。
经济性分析:计算各方案的总投资成本(风电、光伏、储能)。计算年化成本(考虑5年投资回报期)。
比较不同方案的经济性。

优化变量:,,, ,,
目标函数:

限制变量:
电平衡:
电池状态:
电池容量变换范围限制:
电池充放电状态限制:
不能同时充放电:
其他:
优化结果:
表9为增长负荷之后的风光储协调配置方案经济性,从该表中可以看出:
购电量
在这里插入图片描述

在这里插入图片描述

图16联合园区增长负荷之后的风光储协调配置信息
时间段 购电量/kWh 弃风弃光电量/kWh 充电量/kWh 放电量/kWh
0 288.7293216 0 0 0
1 0 0 290.6295812 0
2 0 0 161.0289092 0
3 0 0 280.9827567 0
4 0 0 978.6979893 0
5 0 0 855.2761331 0
6 0 0 0 125.7532545
7 0 0 0 1090.747681
8 0 0 0 458.2943
9 0 0 309.81864 0
10 0 0 0 325.4529367
11 0 0 387.4384707 0
12 0 0 913.0820729 0
13 0 0 253.4907529 0
14 0 0 0 457.116122
15 0 0 0 488.5146185
16 0 0 193.054425 0
17 0 0 212.939244 0
18 0 0 415.5817153 0
19 0 0 726.5380996 0
20 0 0 0 322.9343346
21 0 0 0 623.4313861
22 0 0 0 527.7193439
23 0 0 0 975.6853304
表13联合增长负荷之后的风光储协调运行各时段参数
图13~16分别展示了园区A,B,C和联合园区在增长负荷之后的风光储协调配置信息,

在这里插入图片描述

问题代码展示:A_1_1_A

%A_1_1_A
% 清除工作区
clear;clc;close all;warning off;
%读取参数%正常读取
% P_LOAD=xlsread('附件1:各园区典型日负荷数据.xlsx');
% P_PVWT=xlsread('附件2:各园区典型日风光发电数据.xlsx');
% P_load=P_LOAD(:,2)';
% P_PV=750*P_PVWT(:,2)';%若无法读取直接使用该数据
P_load=[275,275,277,310,310,293,293,380,375,281,447,447,447,405,404,403,268,313,287,288,284,287,277,275];%电负荷/kW
P_PV=750*[0,0,0,0,0,0,0,.0058,.3026,.6020,.7711,.8555,.8531,.7842,.6437,.4242,.0619,0,0,0,0,0,0,0,];%光伏输出/kW%常数
N_T=24;%24h
C_buy=1;%购电价格:元/kWh
C_WT=.5;    %风电价格:元/kWh
C_PV=.4;    %光伏价格:元/kWh% 决策变量
P_buy = sdpvar(1,N_T);%购电功率
P_waste = sdpvar(1,N_T);%弃电功率% 设置约束条件
C=[P_PV+P_buy==P_load+P_waste      %电平衡P_buy >= 0                      % 购电功率非负P_waste >= 0                    % 弃电功率非负];ops = sdpsettings('solver', 'gurobi', 'verbose', 2, 'debug', 1);F = sum(P_PV*C_PV+C_buy*P_buy);
result = optimize(C,F,ops);
if result.problem == 0 % problem =0 代表求解成功disp(value(F));
elsedisp('求解出错');
end% 定义要删除的文件名
filename = 'A_1_1_A.xlsx';
% 使用delete函数删除文件
if exist(filename, 'file')delete(filename);disp(['已成功删除文件:' filename]);
elsedisp(['文件 ' filename ' 不存在。']);
end
data = table({'购电量/kWh', '弃风弃光电量/kWh', '总供电成本/¥', '单位电量平均供电成本/(¥/kWh)'}', ...[sum(value(P_buy)), sum(value(P_waste)),value(F),value(F)/sum(P_load)]', 'VariableNames', {'经济参数', 'A_1_1_A'});
writetable(data, 'A_1_1_A.xlsx');figure('Position', [880, 60, 800, 500]);  % [左下角 x 坐标, 左下角 y 坐标, 宽度, 高度]
b=bar([value(P_buy);value(P_PV);-value(P_waste)]','stacked');hold on;
b(1).FaceColor = [156/255 168/255 184/255];
b(2).FaceColor = [122/255 103/255 71/255];
b(3).FaceColor = [201/255 156/255 78/255];
plot(P_load,'k','LineWidth',2);hold on;
legend('购电','光伏','弃电','电负荷','Orientation', 'horizontal')%横着放图例
legend('Location', 'NorthWest');  % 根据需要选择位置
title('A\_1\_1\_A电平衡');saveas(gcf, 'A_1_1_A电平衡.jpg');

这篇关于【全部更新完毕】2024电工杯A题数学建模详细思路代码文章分享的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1000222

相关文章

Ilya-AI分享的他在OpenAI学习到的15个提示工程技巧

Ilya(不是本人,claude AI)在社交媒体上分享了他在OpenAI学习到的15个Prompt撰写技巧。 以下是详细的内容: 提示精确化:在编写提示时,力求表达清晰准确。清楚地阐述任务需求和概念定义至关重要。例:不用"分析文本",而用"判断这段话的情感倾向:积极、消极还是中性"。 快速迭代:善于快速连续调整提示。熟练的提示工程师能够灵活地进行多轮优化。例:从"总结文章"到"用

2024年流动式起重机司机证模拟考试题库及流动式起重机司机理论考试试题

题库来源:安全生产模拟考试一点通公众号小程序 2024年流动式起重机司机证模拟考试题库及流动式起重机司机理论考试试题是由安全生产模拟考试一点通提供,流动式起重机司机证模拟考试题库是根据流动式起重机司机最新版教材,流动式起重机司机大纲整理而成(含2024年流动式起重机司机证模拟考试题库及流动式起重机司机理论考试试题参考答案和部分工种参考解析),掌握本资料和学校方法,考试容易。流动式起重机司机考试技

poj3468(线段树成段更新模板题)

题意:包括两个操作:1、将[a.b]上的数字加上v;2、查询区间[a,b]上的和 下面的介绍是下解题思路: 首先介绍  lazy-tag思想:用一个变量记录每一个线段树节点的变化值,当这部分线段的一致性被破坏我们就将这个变化值传递给子区间,大大增加了线段树的效率。 比如现在需要对[a,b]区间值进行加c操作,那么就从根节点[1,n]开始调用update函数进行操作,如果刚好执行到一个子节点,

hdu1394(线段树点更新的应用)

题意:求一个序列经过一定的操作得到的序列的最小逆序数 这题会用到逆序数的一个性质,在0到n-1这些数字组成的乱序排列,将第一个数字A移到最后一位,得到的逆序数为res-a+(n-a-1) 知道上面的知识点后,可以用暴力来解 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#in

hdu1689(线段树成段更新)

两种操作:1、set区间[a,b]上数字为v;2、查询[ 1 , n ]上的sum 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdl

【专题】2024飞行汽车技术全景报告合集PDF分享(附原数据表)

原文链接: https://tecdat.cn/?p=37628 6月16日,小鹏汇天旅航者X2在北京大兴国际机场临空经济区完成首飞,这也是小鹏汇天的产品在京津冀地区进行的首次飞行。小鹏汇天方面还表示,公司准备量产,并计划今年四季度开启预售小鹏汇天分体式飞行汽车,探索分体式飞行汽车城际通勤。阅读原文,获取专题报告合集全文,解锁文末271份飞行汽车相关行业研究报告。 据悉,业内人士对飞行汽车行业

活用c4d官方开发文档查询代码

当你问AI助手比如豆包,如何用python禁止掉xpresso标签时候,它会提示到 这时候要用到两个东西。https://developers.maxon.net/论坛搜索和开发文档 比如这里我就在官方找到正确的id描述 然后我就把参数标签换过来

高效录音转文字:2024年四大工具精选!

在快节奏的工作生活中,能够快速将录音转换成文字是一项非常实用的能力。特别是在需要记录会议纪要、讲座内容或者是采访素材的时候,一款优秀的在线录音转文字工具能派上大用场。以下推荐几个好用的录音转文字工具! 365在线转文字 直达链接:https://www.pdf365.cn/ 365在线转文字是一款提供在线录音转文字服务的工具,它以其高效、便捷的特点受到用户的青睐。用户无需下载安装任何软件,只

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n

uva 10014 Simple calculations(数学推导)

直接按照题意来推导最后的结果就行了。 开始的时候只做到了第一个推导,第二次没有继续下去。 代码: #include<stdio.h>int main(){int T, n, i;double a, aa, sum, temp, ans;scanf("%d", &T);while(T--){scanf("%d", &n);scanf("%lf", &first);scanf