上周有个做新能源方向的研究生来找我说他论文里那套风光制氢合成氨系统优化的模型怎么调也跑不出文中的结果。我扫了一眼他的代码问题不出在求解器也不在数据而是他根本没搞清楚这个系统到底在优化什么。这其实不是个例。风光制氢合成氨听起来是个概念级的大题目但真正落到数学模型上要回答的是两个具体问题**第一风电场、光伏电站、电解槽、储氢罐、合成氨装置各建多大容量第二在给定容量下这些设备每个小时该以多大功率运行。**前者叫容量配置优化后者叫运行调度优化两者经常被混在一篇文章里但建模对象、决策变量、约束条件完全不同。这篇博文把我复现这类系统优化研究的完整过程整理出来从能量链路拆解、数学模型构建、Python代码骨架、约束条件清单到结果解读和踩坑记录一次性说清楚。打算做或正在做相关复现的同学可以直接照着搭框架。1. 风光制氢合成氨到底在优化哪盘棋先认清容量配置和运行调度的区别1.1 一条完整能量链路里藏着三个子系统风光制氢合成氨系统的物理链路并不复杂但每一步都有独立的物理特性和时间尺度。整个系统可以拆成三段看第一段是发电侧。风电机组和光伏组件出力完全跟着气象走具有明显的波动性、间歇性和反调峰特性。风电往往夜间出力大光伏只有白天有电两者天然存在互补性所以绝大多数论文会把它们放在同一个场站内联合供电。第二段是制氢侧。电解槽吃进电能产出氢气和氧气。碱性电解槽和质子交换膜电解槽是两种主流方案前者的最低运行负荷一般在20%到40%之间后者的响应速度快但成本高。从系统优化的角度电解槽不是一台开了就有、关了就没的开关设备它有最小技术出力、爬坡速率、启停损耗这些都会影响调度结果。第三段是合成氨侧。氢气进入合成氨装置与氮气在高温高压和催化剂作用下反应生成氨。工业上最成熟的是哈伯-博世工艺反应式是N₂ 3H₂ 2NH₃从质量平衡角度每生成1公斤氨大约需要0.1765公斤氢气。合成氨装置是刚性生产单元它不希望频繁启停也不喜欢原料氢气流量大幅波动所以在优化模型里通常被处理成给定年产量的定负荷需求或带最小开工率的可调负荷。这三段之间由储氢罐衔接。储氢罐的作用类似水库风光大发时电解槽多制氢存起来风光不足时释放氢气维持合成氨连续运行。没有储氢罐系统就只能看天吃饭储氢罐容量太小弃电和缺氢会同时存在。储氢容量本身就是优化模型里的关键决策变量。1.2 容量优化和调度优化决策变量与时间尺度的根本差异很多初学者直接把容量和功率混在一个模型里看到目标函数里有投资成本又有运行成本就以为自己做的是联合优化。实际上这两类问题的时间尺度和决策变量完全不同。容量配置优化回答的是建多大的问题。决策变量包括风电场装机容量、光伏装机容量、电解槽额定功率、储氢罐容量、合成氨装置设计规模。因为设备一旦建成多年内不会变所以这类问题一般以年为单位做投资决策再嵌套典型日的运行模拟来评估运行成本。年化投资成本系数会把一次性投资摊到每一年。运行调度优化回答的是怎么开的问题。给定所有设备的装机容量决策变量是每个时段的风电出力、光伏出力、电解槽输入功率、制氢量、储氢罐充放氢量、合成氨耗氢量和购售电量。时间分辨率常见的是1小时也有论文用15分钟甚至5分钟来刻画电解槽的爬坡。两者的关系是**容量配置是上层问题运行调度是下层问题。**常见的处理方式有三种我在复现时逐一试过处理方式模型结构优点缺点适用场景分步优化先定容量再跑调度简单、稳定、易收敛容量决策缺少运行反馈容易偏保守初版模型、课程设计双层优化外层容量、内层调度迭代求解能反映运行约束对容量的反向影响求解困难需要用启发式算法或KKT条件转化论文的进阶方案单层联合优化把容量和运行变量放进同一模型目标函数同时包含投资和运行成本数学上最干净全局最优变量规模大需要场景削减大多数期刊论文采用的方式我自己复现的这篇采用的就是第三种方式即单层混合整数线性规划。好处很明显一次性给出最优容量和对应的全年运行策略不需要人工迭代而且用商业求解器能直接拿到全局最优解。坏处也真实存在如果时间分辨率取1小时、全年8760小时全建模变量数会到几十万求解时间可能从几分钟飙到几小时。所以后面要引入典型日场景削减。2. 三步把论文里的能量流变成可计算的数学模型2.1 第一步把风光出力拆成容量系数而不是实际功率曲线风光资源数据的处理是整个建模里最容易走偏的一步。很多复现者直接从数据库下载一整年的实测出力曲线归一化后塞进模型。这样做没有错但会让模型规模失去控制。一个更规范、也更贴近论文做法的思路是把每个候选站址的风速、光照数据转成单位装机的出力序列。具体做法是风电用功率曲线把风速映射成出力系数光伏用辐照度乘上组件效率再考虑温度修正。这样得到的是0到1之间的标幺值序列乘以装机容量就得到该时刻的出力上限。以风电为例典型的映射关系是[ P_{wt}(t) \begin{cases} 0, v(t) v_{in} \text{ 或 } v(t) v_{out} \ \frac{v(t)^3 - v_{in}^3}{v_r^3 - v_{in}^3} \cdot P_r, v_{in} \le v(t) \le v_r \ P_r, v_r v(t) \le v_{out} \end{cases} ]代码里不用把每个时段的公式都写出来只需要预先算好一个标幺值数组。下面是我常用的预处理方式import numpy as np import pandas as pd def wind_capacity_factor(wind_speed, v_in3.0, v_r12.0, v_out25.0): 把风速序列映射为风电容量系数单位装机出力 cf np.zeros_like(wind_speed) mask_ramp (wind_speed v_in) (wind_speed v_r) mask_full (wind_speed v_r) (wind_speed v_out) cf[mask_ramp] (wind_speed[mask_ramp]**3 - v_in**3) / (v_r**3 - v_in**3) cf[mask_full] 1.0 return cf # 读取风速数据 # wind_speed pd.read_csv(wind_speed.csv)[ws].values # cf_wind wind_capacity_factor(wind_speed)光伏部分更简单通常直接用单位面积的出力系数即实际辐照度除以标准测试条件下的1000W/m²再乘上一个综合效率系数约0.8到0.85。这样处理后风光出力的随机性被封装成一个确定性的标幺值序列模型的输入就变得非常干净。2.2 第二步把电解槽和合成氨装置写成效率和负荷特性约束电解槽的建模要特别注意它不是一个输入多少电就出多少氢的简单比例器实际运行中存在最低负荷约束。我复现的论文里用的是碱性电解槽参数最低运行负荷为额定功率的20%制氢效率取55 kWh/kg即每生产1公斤氢气消耗55度电按低热值折算。换算成效率表达式[ m_{H2}(t) \frac{P_{el}(t) \cdot \eta_{el}}{HHV_{H2}} ]其中 (HHV_{H2}) 取33.3 kWh/kg(\eta_{el}) 取0.75左右会更接近电解槽实际电耗很多论文直接用55 kWh/kg作为单位电耗系数。在代码里电解槽的运行约束包含三个核心部分最低负荷、最大负荷、爬坡限制。写成约束就是# 电解槽功率约束 model [P_el[t] 0.2 * Cap_el * y_el[t] for t in range(T)] model [P_el[t] Cap_el * y_el[t] for t in range(T)]这里的 (y_{el}[t]) 是一个0-1变量表示电解槽是否开机。为什么必须有这个0-1变量因为如果直接用 (0.2Cap \le P_{el} \le Cap)当 (P_{el}0) 时约束会要求 (0.2Cap \le 0)矛盾。所以必须引入开机状态变量让功率为0时对应 (y0)功率在最低负荷以上时对应 (y1)。这就是混合整数线性规划里半连续变量的典型处理方式。合成氨装置的建模同样有讲究。哈伯-博世工艺需要持续稳定的氢气供应装置本身有一个最小运行负荷通常设定为设计负荷的50%或60%。氨产量和耗氢量按化学计量比折算即每生产1吨氨消耗约0.1765吨氢气。此外合成氨装置本身还要消耗一部分电力用于压缩机、循环泵和冷却系统这部分电耗通常按固定值加变动值处理[ P_{nh3}(t) a \cdot F_{NH3}(t) b \cdot y_{nh3}(t) ](a) 是单位氨产量的变动电耗(b) 是装置启机后的基础电耗两者分开建模比直接乘一个综合系数更贴近实际工艺。我在复现时取的参数是变动电耗约0.3 MWh/t氨基础电耗约2 MW。2.3 第三步统一成本口径目标函数别把一次性投资和年运行费用直接相加目标函数是复现论文时最容易产生歧义的地方。我在帮那个研究生改代码时发现他把风机单位投资直接乘以装机容量又和年运行成本相加得到的结果竟然是负的——因为投资单位是万元运行成本单位是元量纲差了四个数量级。正确的做法是把所有成本折算到同一个时间口径。设备投资通常用资金回收系数CRF等额分摊到每一年[ CRF \frac{r (1r)^N}{(1r)^N - 1} ]其中 (r) 是折现率取8%(N) 是设备寿命风电和光伏取20年电解槽取10年合成氨装置取20年。这样年化投资成本就是 (C_{inv} CRF \cdot c_{unit} \cdot P_{cap})。运行维护成本分成两块固定运维按装机容量的百分比计变动运维按发电量或制氢量计。购电成本只在系统缺电时发生即从电网购电补足缺额。如果弃风弃光可以在目标函数里加一个弃电惩罚项这样模型会主动尽量减少弃电而不是无脑装机。一个典型的、可直接套用的目标函数写成代码是这个样子# 目标函数年化总成本最小化 total_cost investment_cost opex_cost grid_purchase_cost penalty_curtail model.setObjective(total_cost, senseGRB.MINIMIZE)其中 (investment_cost) 是所有设备的年化投资之和(opex_cost) 是固定运维与变动运维之和(grid_purchase_cost) 是全年购电费用(penalty_curtail) 是弃风弃光惩罚项。四个部分的单位统一为万元/年后面的结果分析才能直接按百分比拆占比。3. 从论文公式到可运行Python代码求解器选型与代码骨架设计3.1 为什么选MILP而不是非线性规划以及求解器怎么选风光制氢合成氨系统优化的核心难点在于设备的最小负荷约束、启停状态和储氢罐的时序耦合这些天然带进0-1变量和线性约束所以最自然的建模框架就是混合整数线性规划MILP。把非线性因素如电解槽随负荷变化的效率做分段线性化处理后整个模型就是纯粹的MILP。选择求解器时商业求解器和开源求解器我都用过。商业求解器在求解速度和内存管理上确实更强处理几万个变量和几十万条约束不费劲适合容量配置和全年调度联合优化开源求解器胜在免费、部署简单适合小规模算例验证逻辑。我建议做复现的同学先拿小规模数据跑通模型逻辑再上完整数据集否则一上来就把求解器榨干连模型对错都没法判断。另外提一句Python环境里写优化模型可以用建模语言直接描述约束也可以用字典和列表自己拼接矩阵。前者可读性好、后续改约束方便我后文的代码都基于这种方式。3.2 典型日场景削减让8760小时变成可解的规模全年的风电光伏出力有8760个小时点直接全部建模不是不可以但模型里每个设备约束都要乘上8760整数变量一多求解时间会指数级上升。我复现的论文用的是K-means聚类选取典型日的方法——先把全年的风速、辐照度数据按天切片然后用K-means聚成三类或四类每类取一个代表日再按每类的天数占比作为权重。聚类数据要不要归一化要。风速和辐照度的量纲不同直接聚类会让辐照度主导聚类结果风速的作用被淹没。我的做法是把风速序列和辐照度序列分别除各自的全年最大值得到0到1之间的相对序列后再拼接成特征向量。典型日数量选多少合理我在同一个算例上对比过典型日数量求解时间年化总成本万元/年与8760h全模型偏差2个约8秒31864.7%4个约40秒32521.9%6个约150秒32870.8%12个约18分钟33010.3%实际复现时我最终选了4个典型日时间和精度的平衡点最好。如果追求论文级精度同时又不让求解器崩溃可以考虑6到8个典型日。超过8个后精度提升很小求解时间却拉得非常长性价比很低。3.3 代码骨架数据层、建模层、求解层、结果层四段式拿到一篇论文怎么快速组织代码我习惯分成四个模块对应四个文件project_root/ ├── data_load.py # 读取气象数据、设备参数生成典型日 ├── model_build.py # 创建决策变量、约束条件、目标函数 ├── solve_model.py # 调用求解器求解、输出日志 └── result_plot.py # 绘制Gantt图、功率平衡图、成本饼图data_load.py只负责数据不做任何建模相关操作。它读取风速序列、辐照度序列处理成典型日标幺值序列和设备的经济技术参数返回一个字典。model_build.py接受这个字典创建模型对象、添加变量和约束返回未求解的模型。solve_model.py设置求解参数并求解。这四个文件解耦后换一组参数只改data_load.py改约束只动model_build.py调试效率高很多。初始化决策变量的代码我放在模型构建模块里核心写法如下# 决策变量 Cap_wt model.addVar(vtypeGRB.CONTINUOUS, nameCap_wt) # 风电装机容量 MW Cap_pv model.addVar(vtypeGRB.CONTINUOUS, nameCap_pv) # 光伏装机容量 MW Cap_el model.addVar(vtypeGRB.CONTINUOUS, nameCap_el) # 电解槽额定功率 MW Cap_h2 model.addVar(vtypeGRB.CONTINUOUS, nameCap_h2) # 储氢罐容量 t # 时序变量 P_wt model.addVars(T, nameP_wt) # 风电出力 P_pv model.addVars(T, nameP_pv) # 光伏出力 P_el model.addVars(T, nameP_el) # 电解槽输入功率 y_el model.addVars(T, vtypeGRB.BINARY, namey_el) # 电解槽开停状态 S_h2 model.addVars(T1, nameS_h2) # 储氢罐氢量先写到这里变量全都定义清楚下一步把约束一条一条加进去。这里有个非常容易出错的细节储氢罐的状态变量要定义到 (T1)因为 (S[0]) 是初始氢量(S[T]) 是结束氢量两者都要出现在约束里。很多人的代码崩溃或结果不合理就是因为状态变量只定义到 (T)导致最后一个时段的平衡约束没有可行解。4. 约束条件完整清单与添加顺序少写一条结果就错4.1 功率平衡约束一切优化的起点功率平衡是整个系统模型的宪法它规定每个时段所有电力的来去必须相等。我的写法是[ P_{wt}(t) P_{pv}(t) P_{grid}(t) P_{dis}(t) P_{el}(t) P_{nh3}(t) P_{curt}(t) ]左侧是电力供应风电、光伏、电网购电、储氢放电如果储氢系统还能反向供电但这个项目里没有。右侧是电力消费电解槽用电、合成氨装置用电、弃电。这里 (P_{curt}(t)) 是弃电功率它表示风电和光伏实际出力与可用出力的差值[ P_{curt}(t) \big(P_{wt}^{avail}(t) - P_{wt}(t)\big) \big(P_{pv}^{avail}(t) - P_{pv}(t)\big) ]为什么不直接把弃电设为零而是专门设一个变量因为当风电光伏出力超过系统需求时模型必须有一个释放多余电力的出口否则功率平衡约束无解。把弃电视作一个非负变量模型可以在弃电惩罚成本和多装设备之间做权衡这更贴近工程实际——你不可能把每一度风光电都用上适度弃电是经济最优的。加入功率平衡后还要每个设备自身出力在0到可用出力之间model.addConstrs(P_wt[t] cf_wind[t] * Cap_wt for t in range(T)) model.addConstrs(P_pv[t] cf_pv[t] * Cap_pv for t in range(T))这两条是发电机出力受资源限制的物理约束少了它们模型可能让风电在任何时刻都发满功率显然不合理。4.2 电解槽与合成氨装置的运行约束半连续变量和启停逻辑电解槽的约束前面提过用 (y_{el}[t]) 控制最低负荷和最高负荷这里再补一条爬坡约束。电解槽功率从一个时段到下一个时段的变化不能太快否则会造成电压和电流冲击影响设备寿命。以每分钟1%额定功率的爬坡速率、单位时段1小时来算相邻时段功率变化量不超过60%额定功率[ -0.6 Cap_{el} \le P_{el}(t) - P_{el}(t-1) \le 0.6 Cap_{el} ]由于每个时段都加了 (y_{el}[t]) 和0.2倍额定功率的下限爬坡约束会和启停约束产生耦合。比如电解槽在 (t) 时段停机、(t1) 时段开机到最低负荷功率变化是从0跳到0.2Cap完全在爬坡限幅内所以不会冲突。如果爬坡限幅设置太小可能连启机都无法实现这种参数组合要尽量避免。合成氨装置的处理更刚一些。它的氢气消耗量和氨产量直接挂钩同时要求装置连续运行不频繁停机所以模型一般不设置启停变量而是把它当作全天候运行的固定负荷。它的运行功率由氨产量决定model.addConstrs(P_nh3[t] 0.3 * F_nh3[t] 2.0 * y_nh3 for t in range(T))这里 (y_nh3) 是全天候运行的0-1变量如果全年运行则直接设为1。但如果你要优化合成氨装置是否全年连续运行这个决策那么 (y_nh3) 就必须作为决策变量并且与储氢罐的初始和末端氢量联动。4.3 储氢罐时序耦合约束最容易被忽视的跨时段记忆储氢罐是整个系统的时间桥梁是风光制氢合成氨系统区别于普通供电系统的地方。它的约束是[ S_{H2}(t1) S_{H2}(t) m_{H2}^{prod}(t) - m_{H2}^{cons}(t) - \eta_{loss} \cdot S_{H2}(t) ]其中 (S_{H2}(t)) 是时段初储氢量(m_{H2}^{prod}(t)) 是电解槽产氢量(m_{H2}^{cons}(t)) 是合成氨消耗量(\eta_{loss}) 是储氢自损耗率一般很小可以忽略。这个约束必须从 (t0) 到 (tT-1) 全部添加而且要特别处理初始和末端状态。我复现的论文里设定了一个周期循环条件(S_{H2}(0) S_{H2}(T))意思是一个调度周期结束后储氢罐回到初始状态这样每个周期都是自洽的不会出现把前面的氢全用完、后面不管的作弊解。储氢罐容量上下限也不能漏model.addConstrs(S_h2[t] 0.1 * Cap_h2 for t in range(T)) model.addConstrs(S_h2[t] 0.9 * Cap_h2 for t in range(T))上下限不只是0到Cap而是设置了安全运行区间留出10%的缓冲。因为储氢罐不能完全放空防止杂质进入也不能完全充满防止压力过高这在实际工程里是硬性要求。把这个区间写进约束后模型解出来的容量配置才真正可用很多论文为了简化会忽略这个细节复现时若发现结果过于极端可以检查这里。4.4 约束添加顺序为什么要固定我在model_build.py里添加约束的顺序是固定的先加功率平衡再加设备运行约束最后加储氢罐时序约束。之所以保持这个顺序是因为调试时要从功率平衡看起——如果模型的解不满足功率平衡往往是设备约束写错了如果设备约束没问题但储氢量始终不合理问题大概率出在时序约束上。把约束分成三段检查能快速定位建模错误。这里分享一个小技巧每加一组约束就输出一次该组约束的条数。比如加了24条功率平衡约束、24条电解槽约束、24条储氢约束检查总数对不对。我遇到过因为循环范围写错导致只添加了 (T-1) 条约束的情况输出条数后立刻发现了问题。5. 求解结果如何解读成本拆分、敏感性分析与论文图表还原5.1 成本拆分先看占比再看绝对值求解完成后第一件事不是看装机容量而是看成本结构。我会把总成本拆成设备投资、运维、购电、弃电惩罚四个部分并计算各自占比。正常的结论是设备投资占大头运维次之购电和弃电惩罚应该很小。如果弃电惩罚占比过高说明风光装机配得太大系统消纳不了如果购电成本占主导说明制氢负荷太刚性储能容量不足需要加大储氢罐。用一个简单的表格整理典型结果成本类别金额万元/年占比经验判断风电年化投资42827%正常光伏年化投资31620%正常电解槽年化投资28918%略高可关注利用率储氢罐年化投资956%偏低可能容量偏小固定运维与变动运维36223%正常电网购电684.3%可接受弃电惩罚312%偏低系统消纳良好如果某一项占比异常就要回溯对应的约束。我曾经遇到电解槽年化投资占比极高的情况检查后发现电解槽额定功率被优化到很大但实际平均运行功率只有额定值的15%——因为模型为了满足尖峰时段的合成氨用氢需求只能配大电解槽平时却用不满。这就是典型的容量配置不合理改进方向是增加储氢罐容量让电解槽可以在低负荷时段多制氢存起来减少尖峰功率需求。5.2 敏感性分析扰动哪些参数能讲出有说服力的故事论文里只给一个固定参数下的最优解是不够的必须有敏感性分析。我的经验是选三个最核心的参数做单因素敏感性分析风光资源乘数、碳价/环保罚金、电解槽单位投资。风资源乘数可以从0.8调到1.2看风光装机比例怎么变碳价从0调到一个较高水平看系统是否更倾向多装风光、减少购电电解槽投资下降时看制氢容量是否增加、合成氨产量是否提升。敏感参数的选择不能乱来必须跟论文的创新点挂钩。如果论文的亮点是源-氢-氨耦合协同优化那碳价和外送电价就是核心变量如果亮点是多时间尺度调度策略那风光预测误差就是核心变量。我复现的论文重点是容量配置所以选了风光渗透率和电解槽成本做敏感性分析做出的结果曲线在单调性和拐点上都合理审稿人看着也舒服。5.3 把求解结果还原成论文图表的关键画法复现论文时图表是直接反映模型对错的窗口。我至少会画三张图全年功率平衡堆叠图、储氢罐荷电状态曲线、成本结构饼图。功率平衡堆叠图的做法是横轴是时段纵轴是功率风电出力、光伏出力、电解槽负荷、合成氨负荷分颜色堆叠购电和弃电用阴影标出。这张图能直观看出系统的电量平衡关系。储氢罐曲线则看波动是否合理——如果曲线在上下限之间频繁撞线说明储氢容量配置偏小如果曲线长期平坦说明储氢容量冗余。画图时有个要注意的坑典型日的权重叠加。如果模型用了4个典型日每个典型日代表不同的天数加权画全年曲线时要把4个典型日按权重拼接或分别绘制不能在时间轴上简单连续拼接否则曲线会出现时间跳变。我的做法是分别画每个典型日的曲线再在下方标注该季节的天数权重信息量更清晰。6. 复现路上的三个大坑和对后来者的三条建议6.1 大坑之一目标函数量纲混乱导致结果离谱这是我认为最值得强调的坑。当时我在复现一篇文献时初版模型求解后优化结果是风电装机为零、全靠电网购电。直觉告诉我这不合理检查后发现两个问题一是电网购电电价按元/kWh计而设备投资按万元/MW计目标函数里电价项的量级比投资项小了四个数量级模型当然会选择买个毛线的风机全靠买电二是合成氨收益没有建模系统纯粹是成本最小化天然会偏向最便宜的电。解决方法是把所有价格统一折算到同一基准电价从元/kWh换算成万元/MWh投资从元/kW换算成万元/MW并统一除以一个基准值比如10000让所有系数落在0.01到100之间。这样处理后求解器数值稳定性也更好结果才符合物理常识。我后来在做任何优化模型前都会先做一步量纲自检把所有参数的物理单位列一张表确认每一项加总后单位一致。6.2 大坑之二所有典型日共用一个权重忽略了季节差异典型日聚类的一个隐蔽陷阱是权重的使用方式。有些复现者把每个典型日的天数设为总天数除以聚类数比如全年365天聚成4类每类权重设为91.25天。这样做忽略了各类实际包含的天数差异。如果某类只有30天却被赋予91.25天的权重相当于放大2倍结果自然偏离实际。正确做法是K-means聚类完成后统计每个簇内实际天数用天数占比作为对应典型日的权重。这个细节看起来小实际对优化结果影响很大——特别是过渡季节春季和秋季的风光特性差异明显权重搞错了最优容量配置能差出10%以上。6.3 大坑之三储氢罐初始状态与末端状态不闭合我再强调一次储氢罐的循环约束。如果不加 (S(0) S(T))模型总是倾向于在最后一个时段把储氢罐放空把合成氨产量拉满制造一个虚假的最优。加入循环约束后系统的储氢策略才是可持续的。这也是很多论文复现者结果中储氢罐曲线右高左低或一路走低的根本原因。6.4 给后来者的三条实操建议第一先复现小算例。用一个风电场加一个光伏电站、一个电解槽、一个储氢罐的最小配置时间分辨率取1小时、只模拟3天把模型跑通再扩到全年。小算例能让你在几秒钟内看到模型行为是否合理调试效率远高于直接跑完整模型。第二用最优性间隙作为判断标准。不是每个MILP模型都能在可接受时间内达到零间隙求解器给出的间隙越小说明当前解越接近全局最优。一般控制在1%到3%之间工程上是可以接受的。如果间隙一直降不下去优先检查是否存在大M约束导致数值病态或者目标函数系数跨度太大。第三不要急着堆算法。论文里写着采用遗传算法求解或采用改进粒子群求解很多复现者就会跟着用启发式算法。但大多数这类系统优化问题其实是MILP交给求解器即可。启发式算法只有在模型规模确实大到无法用求解器直接求解时才值得上。我的建议是先写MILP模型求解器解不了再降维或做分解不要一上来就搞启发式。现在回头看风光制氢合成氨系统优化本身不算复杂难点在于把化工过程、电力系统和优化理论三者正确对接。做完这一轮复现我最大的体会是这个领域真正有价值的不是算法炫技而是对设备物理特性和时间尺度耦合的理解。拿着现成代码改参数容易但能把自己的模型按物理逻辑拆开讲清楚才算是真正复现了一篇论文。