目录
1 概述
2 案例
3 Matlab实现
3.1 题目及分析
3.2 Matlab代码
4 七个题全部代码(值得学习)
1 概述
“碳中和”目标驱动下未来电力系统必将是高比例可再生能源电力系统,可再生能源输出功率强随机波动性导致系统运行中功率实时平衡困难;储能被认为是保障系统功率实时平衡的有效手段,由于储能成本相对昂贵,利用储能平衡系统功率将增加系统运行成本;下面以高比例风电电力系统为例,探究“供给侧” 低碳化转型对电力系统运行经济性、可靠性影响。待研究系统包含火电、风电、储能和负荷,火电机组 3 台、装机容量 1050MW; 某日风电、负荷归一化功率(1.0p.u.风电对应其装机容量,1.0p.u.负荷对应最大负荷功率)数据见附件 1,风电渗透率(最大风电功率与最大负荷功率之比)递增可能造成系统弃风、失负荷,影响系统功率平衡。
定义:系统单位供电成本=系统发电总成本/系统总负荷电量,发电总成本= 火电成本+风电成本+储能成本+弃风损失+失负荷损失,其中:
⚫ 火电成本包括运行成本、碳捕集成本,其中火电运行成本由运行维护成本和发电煤耗成本构成,发电煤耗与其出力关系:
,F 为机组煤耗量(kg/h),P 为机组出力/MW;运行维护成本按照 0.5 倍煤耗成本考虑,碳捕集成本取决于碳排放量及碳捕集单价,火电机组相关参数如附表 1 所示,电煤价格为 700 元/t。
⚫ 风电成本仅考虑运维成本,相关参数如附表 2。
⚫ 储能成本由投资成本、运维成本构成,相关参数如附表 3。注:在计算每天成本时,需将投资成本平摊至每天,即平均每天投资成本=总投资成本/运行年限/365 天。
⚫ 弃风损失按 0.3 元/kWh 计算,失负荷损失按 8 元/kWh 计算。
2 案例
假设系统日负荷功率最大值 900MW,单位碳捕集成本分别为 0 元/t、60 元 /t、80 元/t、100 元/t,摄动风电渗透率,分析计算以下问题:
1. 无风电接入,火电以最小成本运行,绘制机组日发电计划曲线,计算系统单位供电成本,将结果填于表 1 相应栏(保留三位有效数字)。
2. 风电装机 300MW、替代机组 3 时,系统功率平衡发生什么变化?弃风电量多少?在此场景下,为减少弃风又不失负荷,风电接入装机容量可以降低多少?
3. 风电装机 600MW、替代机组 2 时,系统功率平衡又发生什么变化?在此场景下,为不失负荷,风电接入容量可增加多少?
4. 针对上述 2-3 风电替代场景,考虑上述 4 种碳捕捉成本,系统按最低发电成本供电,计算系统单位供电成本,并将相关计算结果填入表 2、3。
5. 风电装机 900MW、替代机组 2、3 时,失负荷电量多少?为不失负荷,需要配置的最小储能容量是多少(储能充放电效率 90%)?考虑储能成本、单位碳捕捉成本(取 60 元/t),此时系统单位供电成本多大?
6. 当负荷功率不变,试分析风电替代容量递增給系统可靠供电带来哪些挑战?为保障可靠供电,系统单位供电成本发生了什么变化?结合上述计算结果进行定量分析。
拓展题:
7. 针对附件 2 所示的十五天负荷功率(最大值 1200MW)、风电功率(装机
容量 1200MW),在风电替代火电机组 2、3 场景下,系统功率平衡存在什么问题?
请设计可能的功率平衡解决方案,并论述方案的可行性和有效性。
如果需要更多真题及论文+代码, 电工杯历年优秀论文
3 Matlab实现
下面我们先以第二题为例,然后给出Matlab代码:
3.1 题目及分析
先考虑风电装机300MW直接替代机组3的发电功率曲线,观察系统功率平衡发生的变化情况。然后在此场景下,以最小发电成本为目标函数,将负荷分配给两台火电机组1、2和风电机组,绘制机组日发电计划图。再将风电装机容量作为一个决策变量,以最小发电成本为目标函数,求取最小发电成本对应的风电装机容量,绘制机组日发电计划图。
1. 风电装机300MW、直接替代机组3时,机组日发电计划曲线如图2.1(1)所示:
系统功率平衡发生变化如下,
功率为正的部分为弃风的功率时段,功率为负的部分为失负荷的时段。有47个时段存在弃风的情况,弃风电量达到698.404MWh;有49个时段存在失负荷的情况,失负荷电量达到279.799MWh。
2. 此时的系统发电成本=火电成本+风电成本+弃风损失+失负荷损失,以最小成本为目标函数,绘制此时的机组日发电计划曲线,如下图所示:
有5个时段存在弃风的情况,弃风电量达到19.988MWh;不存在失负荷的情况。
3. 在此场景下,为减少弃风,则需要又不失负荷。再以最小成本为目标函数,此时求得的风电装机容量仍然为300MW,即不考虑降低风电的接入装机容量。风电的运维成本低,这也是我们提倡不断推进风电并网促进消纳风电的原因。此时的机组日发电计划曲线没有改变,还是图2.2(1)。
3.2 Matlab代码
(1)main1
(2)main2
(3)Yalmip_Cplex
(4) 数据
4 七个题全部代码(值得学习)
点开第一题:
然后再看看第七题: