4. 发电成本与运行特性建模

安全约束经济调度（SCED）的优化结果是否可信，首先取决于机组模型是否真实。若把火电机组简化为一条光滑二次曲线，优化器会给出数学上最优、工程上却无法执行的调度方案：例如让机组在禁止运行的负荷区间长时间停留，或安排它在20分钟内完成需要1小时才能实现的功率爬升。本节围绕五类最关键的运行特性——成本曲线、阀点效应、最小启停时间、爬坡约束、运行禁区——说明它们在SCED中的工程化建模方法。

一、成本曲线：从热耗率到分段线性

火电机组的燃料成本通常由热耗率试验数据拟合得到。最常用形式是二次函数：

C(P) = a·P² + b·P + c   （元/h，P为机组有功出力，单位MW）
典型示例：a = 0.0025，b = 18.5，c = 320，P ∈ [80, 400] MW
P = 300 MW 时：C = 0.0025×90000 + 18.5×300 + 320 = 225 + 5550 + 320 = 6095 元/h

二次函数在SCED中会带来非线性，工程上常将其分段线性化：把[Pmin, Pmax]划分为若干区间，每个区间用一条直线近似，并引入分段出力变量δ_k与0-1状态变量，使成本变为线性表达式，可直接交给混合整数线性规划（MILP）求解。分段数越多精度越高，但整数变量也越多；工程实践中常取3～6段，在精度与求解规模之间折中。

二、阀点效应：不可忽略的脉动波纹

汽轮机进汽阀在部分开启时会产生节流损失，使实际热耗率曲线上叠加一个周期性脉动分量，称为阀点效应。常用正弦项修正：

C(P) = a·P² + b·P + c + |e·sin(f·(Pmin − P))|
典型示例：e = 120，f = 0.035，Pmin = 80 MW
P = 200 MW 时：sin(0.035×(80−200)) = sin(−4.2) ≈ 0.8716
修正项 ≈ 120×0.8716 ≈ 104.6 元/h

阀点效应使成本曲线出现多个局部极小值，是SCED中典型的非凸来源。工程处理方式有两种：一是保留正弦项并用非线性求解器或智能算法处理；二是在分段线性化时加密采样点，把每个“波谷”附近单独划段，用足够细的线性片段逼近非凸形状，同时避免优化器误把波峰当作最优工作点。

三、最小启停时间与启停成本

机组一旦启动，必须持续运行至少T_on小时；一旦停机，必须停够T_off小时才能再启动。SCED若只做单时段优化，会把启停当作无代价的瞬时动作，导致调度方案频繁切换。工程化处理需引入跨时段状态变量u_i,t（运行/停机）和启动、停机指示变量，并施加约束：

Σ u(i,τ) ≥ T_on(i)·[u(i,t) − u(i,t−1)]，τ ∈ [t−T_on+1, t]
Σ [1−u(i,τ)] ≥ T_off(i)·[u(i,t−1) − u(i,t)]，τ ∈ [t−T_off+1, t]

同时把启动成本（热启动、温启动、冷启动分档）计入目标函数。这样，SCED从单时段静态优化扩展为多时段机组组合与经济调度的联合问题，避免“今天启、明天停”的不可执行方案。

四、爬坡约束：连接相邻时段的功率变化

机组出力不能瞬时跳变，需满足向上、向下爬坡速率限制：

P(i,t) − P(i,t−1) ≤ RU(i)·Δt   （向上爬坡）
P(i,t−1) − P(i,t) ≤ RD(i)·Δt   （向下爬坡）
典型示例：RU = 60 MW/h，RD = 60 MW/h，Δt = 15 min
则单时段最大变化量 = 60×0.25 = 15 MW

爬坡约束把各时段耦合起来，是SCED时间维度上最重要的约束之一。工程上还需区分“运行状态爬坡”与“启动爬坡”：机组从停机到并网、从并网到最小出力，其速率往往不同于正常运行爬坡，应分别设置启动爬坡速率与停机爬坡速率，否则会出现“刚并网就要求满出力”的不可行调度。

五、运行禁区：禁止停留的负荷区间

锅炉燃烧稳定性、汽轮机叶片振动、脱硝系统投运温度等原因，使机组在某些负荷区间不允许长期运行，形成运行禁区（也称禁止运行区）。设机组有k个禁区[L_j, U_j]，则出力必须落在可行域内：

P ∈ [Pmin, L1] ∪ [U1, L2] ∪ [U2, Pmax]
典型示例：Pmin=80，Pmax=400，禁区[150,190]、[250,290]
可行区间：80–150、190–250、290–400 MW

在MILP中，常用0-1变量与“大M”法把机组出力分配到某个可行子区间，并保证同一时段只落在一个区间内。禁区建模会显著增加整数变量数量，工程上应只保留真正影响运行的禁区，避免模型规模失控。

六、工程化处理要点

把上述特性放入SCED时，需注意三点。第一，量纲与时间步长统一：爬坡速率单位是MW/h，而调度时段常为15 min或30 min，必须乘以Δt换算为每时段允许的MW变化量，否则约束会松紧失当。第二，非凸与整数变量的规模控制：阀点效应和运行禁区都会引入非凸性，若对所有机组都精细建模，MILP规模会迅速膨胀；工程上可对基荷大机组精细建模，对小机组采用分段线性近似。第三，可行性优先于最优性：当负荷快速上升、多台机组同时受爬坡限制时，SCED可能无可行解，此时应检查是否预留了足够的旋转备用与快速启动机组，而不是简单放宽约束。

综上，机组成本曲线、阀点效应、最小启停时间、爬坡约束与运行禁区共同构成了SCED的机组侧模型基础。只有把这些工程特性准确转化为优化模型中的目标函数项与约束条件，调度结果才既经济又具备现场可执行性。