经济调度问题的建模及求解

您所在的位置:网站首页 爬坡爬坡怎么写扒 经济调度问题的建模及求解

经济调度问题的建模及求解

2024-07-12 21:00:00| 来源: 网络整理| 查看: 265

前言:在上期内容火电机组经济调度建模及求解——基础篇中,我们介绍了基础的火电机组经济调度模型,本期内容我们将利用动态规划解决考虑爬坡约束的火电机组经济调度问题,在此问题中我们将涉及爬坡率的概念:

电厂灵活性被认为是管理电力负荷变化和提供电网支持服务的重要工具。衡量这种灵活性的一个指标是爬坡率,即发电厂燃气轮机发电增加或减少产出的速率。(来源:Google)

当负荷变化时,发电机组需要以一定的斜率从一个出力值变为另外一个出力值。由于火电机组存在机械/热力方面的约束,每分钟的出力不能超过规定的最大变化率限值。当电力系统中负荷出现一个小的增量时,为实现最优的经济调度(发电成本最低),需要在满足爬坡率约束的前提下,降低部分机组的发电量的同时增加其他发电机的出力。因此该问题可以视为一个在时间上前后关联的多阶段决策问题,并可以利用动态规划的思想实现问题求解。(来源:《电力系统发电、运行和控制》)

本期推文内容介绍

爬坡约束的建模以及编程实现、排污成本的计算model.addRange()的用法介绍动态规划中两个连续时间点的约束表述——用于爬坡约束的描述x.sum()对变量求和得到最优解后,利用最优解计算其他目标值——getValue()的用法介绍 Part 1 考虑爬坡约束的火电机组经济调度模型

基础的火电机组经济调度模型为: min ⁡ P g , t T C = ∑ g , t a g P g , t 2 + b g P g , t + c g P g min ⁡ ≤ P g , t ≤ P g max ⁡ ∑ g P g , t ≥ L t \begin{aligned} \min _{P_{g, t}} \mathrm{TC}=& \sum_{g, t} a_{g} P_{g, t}^{2}+b_{g} P_{g, t}+c_{g} \\ & P_{g}^{\min } \leq P_{g, t} \leq P_{g}^{\max } \\ & \sum_{g} P_{g, t} \geq L_{t} \end{aligned} Pg,t​min​TC=​g,t∑​ag​Pg,t2​+bg​Pg,t​+cg​Pgmin​≤Pg,t​≤Pgmax​g∑​Pg,t​≥Lt​​ 在此基础上加入爬坡约束,以避免因为负荷扰动,出现火电机组违反爬坡率的情况: P g , t − P g , t − 1 ≤ R U g P g , t − 1 − P g , t ≤ R D g \begin{aligned} & P_{g, t}-P_{g, t-1} \leq \mathrm{RU}_{g} \\ & P_{g, t-1}-P_{g, t} \leq \mathrm{RD}_{g} \end{aligned} ​Pg,t​−Pg,t−1​≤RUg​Pg,t−1​−Pg,t​≤RDg​​ 式中, P g , t P_{g, t} Pg,t​、 P g , t − 1 P_{g, t-1} Pg,t−1​分别为 t 时刻与 t-1 时刻的发电功率, R U g \mathrm{RU}_{g} RUg​、 R U g \mathrm{RU}_{g} RUg​分别表示火电机组单位时间内增加或减少的出力,即:爬坡率上限。由于火电机组发电需要考虑排污费用,在得到最优调度结果以后,需要计算总的排污成本: E M = ∑ g , t d g P g , t 2 + e g P g , t + f g \begin{aligned} \mathrm{EM}=& \sum_{g, t} d_{g} P_{g, t}^{2}+e_{g} P_{g, t}+f_{g} \end{aligned} EM=​g,t∑​dg​Pg,t2​+eg​Pg,t​+fg​​

Part 2 Python+Gurobi 代码实现 import gurobipy as gp from gurobipy import GRB '参数设定' a_g = [0.12, 0.17, 0.15, 0.19] # 发电二次项系数 b_g = [14.8, 16.57, 15.55, 16.21] # 发电一次项系数 c_g = [89, 83, 100, 70] # 发电常数项 d_g = [1.2, 2.3, 1.1, 1.1] # 排污二次项系数 e_g = [-5, -4.24, -2.15, -3.99] # 排污一次项系数 f_g = [3, 6.09, 5.69, 6.2] # 排污常数项 P_g_min = [28, 20, 30, 20] # 发电量最小值 P_g_max = [200, 290, 190, 260] # 发电量最大值 RU_g = [40, 30, 30, 50] # 爬坡率——UP RD_g = [40, 30, 30, 50] # 爬坡率——Down L_t = [510, 530, 516, 510, 515, 544, 646, 686, 741, 734, 748, 760, 754, 700, 686, 720, 714, 761, 727, 714, 618, 584, 578, 544] # 24小时的负载 '定义模型及名称' model = gp.Model('Cost-Based-DED') '变量' P_g = model.addVars(4, 24, vtype=GRB.CONTINUOUS, name='P_g') # # 边界约束描述——方法1 # model.addConstrs((P_g[i, t]>=P_g_min[i] for i in range(4) for t in range(23)), name="Con_Pg_min") # model.addConstrs((P_g[i, t]


【本文地址】

公司简介

联系我们

今日新闻


点击排行

实验室常用的仪器、试剂和
说到实验室常用到的东西,主要就分为仪器、试剂和耗
不用再找了,全球10大实验
01、赛默飞世尔科技(热电)Thermo Fisher Scientif
三代水柜的量产巅峰T-72坦
作者:寞寒最近,西边闹腾挺大,本来小寞以为忙完这
通风柜跟实验室通风系统有
说到通风柜跟实验室通风,不少人都纠结二者到底是不
集消毒杀菌、烘干收纳为一
厨房是家里细菌较多的地方,潮湿的环境、没有完全密
实验室设备之全钢实验台如
全钢实验台是实验室家具中较为重要的家具之一,很多

推荐新闻


图片新闻

实验室药品柜的特性有哪些
实验室药品柜是实验室家具的重要组成部分之一,主要
小学科学实验中有哪些教学
计算机 计算器 一般 打孔器 打气筒 仪器车 显微镜
实验室各种仪器原理动图讲
1.紫外分光光谱UV分析原理:吸收紫外光能量,引起分
高中化学常见仪器及实验装
1、可加热仪器:2、计量仪器:(1)仪器A的名称:量
微生物操作主要设备和器具
今天盘点一下微生物操作主要设备和器具,别嫌我啰嗦
浅谈通风柜使用基本常识
 众所周知,通风柜功能中最主要的就是排气功能。在

专题文章

    CopyRight 2018-2019 实验室设备网 版权所有 win10的实时保护怎么永久关闭