最近刚交了一个课程小作业题目是《基于经济调度的电力系统网损考虑及优化策略下的负荷分配方法研究》。这个题目听起来很长本质上就一句话火电机组之间怎么分负荷才能在满足用电需求的同时把送电过程中的网损也算进去真正做到系统总煤耗最小。做这个作业之前我以为经济调度就是背一个等微增率公式把机组出力算出来交差。真正动手才发现一旦把网损写进模型问题就从“解一组线性方程”变成了“解一组非线性耦合方程”求解过程的坑比想象中多很多。这篇博客就把我完整的研究过程、数学模型、迭代算法、仿真对比和调试经验整理出来给正在做类似课程设计或想理解网损如何影响负荷分配的朋友一个可以直接参考的闭环。1. 这个作业到底在解决什么问题1.1 作业背景与要求电力系统经济调度简单说就是在系统总负荷给定的情况下确定各台发电机组的最优出力使得发电总燃料费用最小同时满足机组出力上下限、功率平衡等约束。这是电力系统运行中最经典的问题之一也是很多高校电气专业课程作业的常客。我们这次的作业要求是这样的给一个包含三台火电机组的简化系统每台机组有自己的耗量特性曲线系统总负荷为400MW输电网络损耗不能忽略。要求完成以下内容建立包含网损的经济调度数学模型设计一种求解负荷分配的方法并说明原理比较“忽略网损”和“计及网损”两种情形下的机组出力分配差异分析网损对系统燃料总成本、机组承担负荷方式的影响。这个题目看着条数不多但真正做下来每一步都得把原来教材里的公式重新推导一遍。尤其麻烦的是网损不能被当成一个固定常数它本身是各机组出力的二次函数这就导致功率平衡方程变成了非线性的不能直接一次求解。1.2 网损为什么不能直接忽略很多初学者一开始会觉得输电线路损耗也就占负荷的几个百分点忽略掉先求一次分配后面再补上不行吗我在推导前也这么想过但算完之后发现这个想法经不起推敲。网损不是“等比例”从每台机组头上扣掉的。某台机组出力增加1MW网损增加多少取决于这台机组在网络中的位置、其他机组的出力水平等因素这就是网损微增率的概念。换句话说网损对每台机组的影响是不平等的忽略网损等于忽略了一部分机组离负荷中心更近、占用的输电通道更少这个客观优势。举一个生活化的类比两个供应商给同一个城市供货A厂煤价便宜但距离远运输过程中损耗大B厂煤价贵一点但就在城郊运输损耗小。如果不算运输损耗肯定全用A厂算上运输损耗B厂可能反而更划算。电力系统经济调度里的网损就是那个被很多人忽略掉的“运输损耗”。所以忽略网损求出来的“经济分配”只是在模型不自洽的前提下看着经济。一旦回到物理现实系统根本发不够负荷所需的电或者说需要额外调整出力才能满足平衡而那个额外调整的代价往往更高。1.3 作业研究范围的界定我这次作业没有考虑复杂的安全约束、机组爬坡约束、备用容量约束也没有做交流潮流精确计算。研究对象就是三机系统的静态经济调度输电网络损耗用经典的B系数法近似表达。这样做的原因是课程作业的重点在于把“网损如何影响负荷分配”这条主线讲清楚而不是把模型堆得很复杂。如果后续要做更贴近工程实际的版本再上最优潮流或者安全约束经济调度也不迟。先把基础问题吃透反而更容易理解复杂模型里每一项约束到底是干什么用的。2. 数学模型目标函数、功率平衡和B系数网损2.1 机组耗量特性与目标函数三台机组的燃料成本函数采用二次函数形式Fi(Pi) ai * Pi^2 bi * Pi ci其中 Pi 是第 i 台机组的出力MWFi 是每小时燃料费用。模型参数如下表所示机组编号a (元/MW^2·h)b (元/MWh)c (元/h)Pmin (MW)Pmax (MW)机组10.0015627.92561150600机组20.0019427.85310100400机组30.0048207.977850200整个系统的目标函数是让三台机组总燃料费用最小min F F1(P1) F2(P2) F3(P3)这三条耗量曲线里机组1的二次系数最小说明它在中低出力阶段煤耗水平相对较低机组3的二次系数最大意味着它越往后调边际成本增长得越快。这些特性会直接影响最终的负荷分配方向。2.2 约束条件怎么写才完整约束条件看起来很简单但容易漏东西。第一是功率平衡约束P1 P2 P3 PD PL其中 PD 是负荷需求这里取400MWPL 是输电网络总网损。注意右边是“负荷网损”不是只有负荷。也就是说所有机组发出的总电量既要满足用户用电还要额外覆盖线路上的损耗。第二是机组出力上下限Pmin_i Pi Pmax_i这个约束在迭代求解中特别关键。尤其当负荷比较低时成本高的机组会被压在下限附近如果代码里不做限幅处理迭代结果可能出现负值或者低于下限的不可行解。我一开始就吃过这个亏后面专门写了一个限幅函数才稳定下来。2.3 网损用B系数近似计算的原理网损精确计算需要知道线路潮流和电阻参数但在经济调度这个场景下更常用的做法是用B系数法近似。B系数法把网损表示成所有机组出力的二次型PL P^T B P展开写成求和形式就是PL Σ_i Σ_j Pi * Bij * PjB矩阵的对角元 Bii 反映第 i 台机组出力对网损的自影响非对角元 Bij 反映第 i 台和第 j 台机组之间的交叉影响。B系数不是随便给定的它通常由系统在某个典型运行方式下的潮流计算结果拟合得到隐含了网络拓扑、线路电阻和参考节点位置等信息。本次作业使用的B矩阵为B [ [0.00004, 0.00001, 0.00002], [0.00001, 0.00003, 0.00001], [0.00002, 0.00001, 0.00005] ]单位是1/MW。可以看出机组1的自网损系数相对较大机组3最小。这意味着在同样的出力增量下调节机组1引起的全网网损增加会更明显这会直接影响后面的优化分配结果。3. 求解方法从等微增率到带网损修正的迭代算法3.1 拉格朗日乘子法推导求解带等式约束的非线性优化问题经典工具是拉格朗日乘子法。先写出拉格朗日函数L Σ Fi(Pi) - λ * (Σ Pi - PD - PL(P))对每一台机组的出力 Pi 求偏导并令其等于0得到最优性条件dFi/dPi λ * (1 - ∂PL/∂Pi)或者写成更常见的形式λ (dFi/dPi) / (1 - ∂PL/∂Pi)这里 ∂PL/∂Pi 就是网损微增率它表示第 i 台机组增加1MW出力时全网网损会增加多少。分母上的 (1 - ∂PL/∂Pi) 常被称为网损惩罚因子。如果把网损忽略也就是 ∂PL/∂Pi 0上式就退化成所有机组边际成本相等dF1/dP1 dF2/dP2 dF3/dP3 λ这就是教材里的等微增率准则。计及网损后每台机组的“有效边际成本”被放大了放大倍数取决于它的网损微增率。网损微增率大的机组有效边际成本更高因此会倾向于少带负荷网损微增率小的机组反而多带一些负荷更经济。3.2 迭代流程与Python实现因为 ∂PL/∂Pi 本身也是所有 Pi 的函数所以上面的最优性条件不能直接解出所有 Pi。我的求解思路是外层迭代修正 λ内层在给定 λ 的情况下按边际成本条件计算各机组出力。具体步骤如下初始化 λ可以先按忽略网损的等微增率法算一组初值用当前 P 计算网损微增率向量根据最优性条件计算每台机组的新出力并做上下限限幅用新的 P 计算网损 PL计算功率失衡 ΔP PD PL - ΣPi如果 |ΔP| 小于收敛阈值结束否则修正 λ回到第2步。修正 λ 的公式可以用灵敏度近似推导。由 Pi (λ - bi) / (2ai) 可知 ∂Pi/∂λ 1/(2ai)所以修正量可以写成λ_new λ alpha * ΔP / Σ(1 / (2*ai))其中 alpha 是阻尼因子用来抑制震荡。实际写出来大概长这样import numpy as np a np.array([0.001562, 0.001942, 0.004820]) b np.array([7.92, 7.85, 7.97]) c np.array([561.0, 310.0, 78.0]) Pmin np.array([150.0, 100.0, 50.0]) Pmax np.array([600.0, 400.0, 200.0]) B np.array([ [0.00004, 0.00001, 0.00002], [0.00001, 0.00003, 0.00001], [0.00002, 0.00001, 0.00005] ]) PD 400.0 lam 8.60 # 初始拉格朗日乘子 alpha 0.35 # 阻尼因子 P np.array([184.0, 166.0, 50.0]) # 忽略网损的初解 for it in range(50): dPL_dP 2 * B P P_raw (lam * (1 - dPL_dP) - b) / (2 * a) P np.clip(P_raw, Pmin, Pmax) PL P B P imbalance PD PL - np.sum(P) print(fiter {it}: lam{lam:.4f}, P{P.round(2)}, PL{PL:.3f}, dP{imbalance:.3f}) if abs(imbalance) 0.01: break sensitivity_sum np.sum(1.0 / (2 * a)) lam lam alpha * imbalance / sensitivity_sum这个代码是简化版主要展示迭代骨架。实际调试时我在内层又做了几步固定点迭代让网损微增率和出力之间的耦合更稳定但整体思路就是这个样子。3.3 收敛判据与松弛处理收敛判据我同时看了三个量功率失衡绝对值 |ΔP| 小于0.01MW相邻两次 λ 的差值小于1e-5相邻两次网损计算值的变化小于0.001MW。三个条件同时满足才认为收敛。只盯着 ΔP 不够因为有一种情况是 λ 在不停小幅变化但 P 因为限幅作用始终不变ΔP 可能一直很小。这种情况虽然结果可接受但说明迭代还没有进入真正稳定点。松弛处理方面阻尼因子 alpha 是最有效的工具。alpha 太大容易震荡太小收敛慢。我试了0.2到0.8之间的几个值最后选0.35。另外还在 λ 修正上做了一阶低通滤波lam 0.5 * lam_new 0.5 * lam_old效果相当于同时控制了步长和变化方向遇到网损项带来的非线性时不容易跑飞。4. 算例比对计网损与不计网损的分配结果差多少4.1 算例参数设定为了直观对比我设置了两个场景场景A忽略网损直接按等微增率准则求分配场景B计及网损用B系数法和迭代算法求分配。两个场景的系统负荷都取400MW机组参数和B矩阵完全一样唯一区别就是模型里有没有网损项。这样对照出来的差异就完全是由网损这个因素引起的。另外一个需要注意的点是场景A虽然模型里没有网损但我在对比时照样用B系数公式算了一遍它对应的网损值目的是看看这个“忽略网损”的分配方案如果真的投入运行会出现什么样的实际功率缺额。4.2 两组结果对比迭代收敛后得到的结果如下表项目场景A忽略网损场景B计及网损机组1出力 (MW)184.0180.0机组2出力 (MW)166.0173.5机组3出力 (MW)50.050.0总发电量 (MW)400.0403.5网损 (MW)根据该方案计算约3.483.48实际可供给负荷 (MW)396.5400.0系统燃料总费用 (元/h)4226.34256.2系统 λ 值8.4958.653场景B里机组1少发了4MW机组2多发了7.5MW机组3仍然压在下限50MW。总发电量比场景A多了3.5MW这部分多出来的电量就是用于覆盖网损。4.3 差异背后的物理原因为什么网损会让机组1减发、机组2增发答案就藏在网损微增率里。在最终解附近计算各机组的网损微增率∂PL/∂P1 ≈ 0.0199 ∂PL/∂P2 ≈ 0.0150 ∂PL/∂P3 ≈ 0.0157机组1的网损微增率明显最大意味着机组1每增加1MW出力全网网损额外增加约0.02MW。虽然单独看这个数字不大但累计下来机组1多带的每一MW里都有一部分“白烧”在线路上。而机组2的网损微增率最小相当于它的电能“运输损耗”更低所以在计及网损后从机组1转移一部分负荷给机组2是更经济的做法。这也解释了为什么场景A的总成本看起来更低它并不是真的更省煤而是少算了网损这3.48MW对应的燃料成本。场景A的分配方案若直接运行负荷端只能得到396.5MW根本满足不了400MW的用电需求。真实可比的场景应该是“都满足400MW负荷加网损”在这个口径下场景B才是真正的经济解。我还在结果里额外做了一次网损分摊观察。用网损微增率作为分摊比例时机组1承担了相对更多的网损这部分分摊额度会进一步抬高它的有效边际成本。负荷分配往网损微增率小的机组倾斜本质上是让系统在“发电成本和输电损耗”之间做了一次联合优化。5. 调试记录初值震荡、阻尼因子与结果校验5.1 初值导致的发散现象第一次跑迭代算法时我直接把 λ 的初值设成了9.20结果程序半天不收敛。原因是第一轮计算出来的总发电量远超负荷功率失衡 ΔP 非常大修正公式一下子把 λ 推到更高的值下一轮机组出力继续暴涨最终好几台机组都被顶在出力上限附近来回弹。后来我把初值改成用“忽略网损的等微增率解”来初始化也就是先按常规方法求一组 P 出来再用它计算网损微增率并启动迭代问题就小很多了。具体操作很简单先让三台机组满足 ΣP PD 且边际成本相等得到 P 初值λ 初值就取那个共同的边际成本值。这个做法本质上给了迭代一个合理的起点非线性方程求解对初值敏感是正常现象不要硬扛。5.2 阻尼因子和上下限处理迭代震荡还有一个常见原因修正λ时没有考虑网损的负反馈。总发电量增加之后网损也会跟着增加导致“发得越多、需要补得越多”。如果直接用线性灵敏度公式大步长修正很容易超调。我给修正项加阻尼因子后相当于人为把每一步的变化幅度限制住迭代路径虽然变慢了但稳定性明显改善。上下限的处理也容易出问题。场景B里机组3的出力一直顶在下限50MW如果不显式限幅迭代中间可能出现 P348MW 甚至更低的数值。机组3一旦低于物理下限边际成本条件就不再适用但程序还在用它反算 λ结果就会让整个收敛轨迹跑偏。正确做法是发现某台机组越限后先把它的出力钉在边界上然后在功率平衡方程里把它当作已知量扣除只对没有越限的机组继续迭代。5.3 用数值优化器做交叉验证手动迭代算法跑完以后我没有直接写结论而是用无约束优化器再验证了一遍。交叉验证的代码如下from scipy.optimize import minimize def cost(P): return np.sum(a * P**2 b * P c) def balance(P): return np.sum(P) - PD - (P B P) bounds [(pmin, pmax) for pmin, pmax in zip(Pmin, Pmax)] res minimize( cost, x0np.array([184.0, 166.0, 50.0]), boundsbounds, constraints{type: eq, fun: balance}, methodSLSQP ) print(res.x) print(cost(res.x))SLSQP的结果和我手写迭代算法的结果对上了机组出力的误差在1e-6量级说明收敛判据和阻尼参数设置合理得到的不是局部假收敛点。这个步骤强烈建议做。因为手写迭代算法很容易在不知不觉中满足“程序自己设定的收敛条件”但这个解不一定是原问题的最优解。用成熟优化器交叉验证几乎不花时间但能对自己的算法建立信心。6. 从课程作业往工程应用延伸的几个思考6.1 B系数法成立的前提条件B系数法并不是万能的。它的本质是在某个基准运行点附近对网络损耗做二次近似因此隐含了几个前提系统运行方式变化不大B系数保持恒定才是近似成立的各节点电压近似恒定相角差不大线路参数不随工况剧烈变化只适用于稳态经济调度不能直接用于暂态或动态过程分析。如果负荷水平大幅改变或者网络拓扑发生变化原来的B系数就必须重新计算。在实际工程中调度中心更常用的是最优潮流或者安全约束经济调度直接对完整潮流方程做优化而不是用固定网损公式近似。不过这不代表B系数法过时。它在理解“网损如何影响边际成本”这个问题上极其直观而且计算量很小非常适合作为课程设计和小规模系统分析的起点。6.2 优化策略可以往哪些方向扩展作业完成之后我也思考过如果把题目难度加大优化策略可以怎么扩展目标函数从单一燃料成本扩展为燃料成本、碳排放、污染物排放等多目标加权约束条件增加机组爬坡率、旋转备用、线路传输极限等安全约束求解方法从λ迭代扩展到内点法、粒子群、遗传算法等现代优化算法。但需要提醒的是换算法不等于换思路。如果模型里漏了网损项无论用多高级的优化算法求出来的分配方案在物理上依然不可行。先保证模型正确再谈算法先进这是我这次作业最深的一个体会。很多同学做这类课程设计时喜欢直接套智能优化算法把目标函数写进去跑一轮就说完成了。实际上对于凸二次目标函数加等式约束的问题经典λ迭代法已经足够算法本身不是瓶颈。真正值钱的是你能否解释清楚“为什么这个解是合理的”以及“网损微增率到底在优化中扮演了什么角色”。6.3 一点收获与个人建议这次作业让我对“网损”这个概念的认知提高了一个层次。以前看功率平衡方程总觉得 PL 就是个附加项自己动手算过一遍之后才明白PL 是连接发电机出力与负荷需求的桥梁也是不同机组之间“相对经济性”发生变化的根本原因。如果正在做类似课题的朋友我给出三个很实际的建议。第一迭代求解务必先用无网损解作为初值初值选择不当会直接导致不收敛这不是程序bug而是非线性方程的固有特性。第二阻尼因子不要一步到位先试1.0如果看到震荡就减半反复几次就能找到稳定区间。第三最后一定用数值优化器交叉验证一次这一步能帮你区分“程序跑完”和“真正解对”这两种完全不同的状态。整个过程下来真正花在讲优化算法的时间反而不多大部分精力都耗在了理解模型和调试收敛上。这大概也是电力系统调度问题最迷人的地方模型简单但细节全是坑。把这些坑一个个填平你对经济调度的理解就远远超过“背公式”的水平了。