数学建模论文2004B.pdf

1 电力市场的输电阻塞管理电力市场的输电阻塞管理 摘要 摘要 本文首先建立了基本模型 问题 3 4 中的具体模型均由此基本模型的一部分构成 解决了一个实际的电力问题 首先 在多元线性回归模型中 对数据进行成功的回归分析 得到各线路上有功潮 流关于各发电机组出力的近似表达式 然后 设计了四种阻塞费用计算规则 其中规则 三先给出一种简明的费用计算方法 然后引入满意度的概念 将其刻画成发电商得到的 补偿的函数 对满意度用最小最大模型来找到使各发电商满意度都高的补偿规则 兼顾 简明 合理和公平 接着 在下一时段预报负荷为 MW4 982 时 运用机组出力方案预 报模型 求出清算价为 303 元 MWh 机组的出力方案如下 机组 1 2 3 4 5 6 7 8 方案 150 79 180 99 5 125 140 95 113 9 由电力市场交易规则判断出求得的出力方案将发生电流阻塞 进而运用输电阻塞管理 模型对方案进行优化 消除了输电阻塞 优化方案如下 机组 1 2 3 4 5 6 7 8 方案 153 86 66 228 69 83 152 95 80 91 117 对应阻塞费用 B 4852 69 元 当预报负荷为 MW8 1052 时 清算价是 356 元 MWh 分配预案为 机组 1 2 3 4 5 6 7 8 方案 150 81 218 2 99 5 135 150 102 1 117 判断其发生阻塞且无法消除 但不必拉闸限电 运用分层排序模型 在安全且经济 的原则下 得到最优的分配方案 机组 1 2 3 4 5 6 7 8 方案 153 88 228 99 5 152 113 102 1 117 阻塞费用85 2234 B元 在模型扩展中 我们构造出问题 4 的新模型 引入了危险系数 瘫痪损失 风险费用 等概念 巧妙的把要求既安全有经济的双目标规划转化为总费用最小的单目标规划 并 用模拟退火法求解 同时 我们又对表一表二的数据运用多元线性回归模型重新进行 3 种不同的拟合 带常数项拟合 考虑零输入零输出的拟合 剔除异常数据的拟合 得 出结论 尽管拟合方式差别很微小 但对问题结果的影响是巨大的 在灵敏度分析中 在各个机组出力变化 5 的情况下 得出超限系数 i 方案 阻 塞费用的变化很小 从而得出结论 机组出力在小范围内变化 对目标函数值没有影响 2 关键词 多元线性回归 关键词 多元线性回归 可用段容量可用段容量 n人合作对策人合作对策 满意度满意度 危险系数危险系数 瘫痪损失瘫痪损失 风险费用风险费用 问题重述 问题重述 随着我国用电紧张的缓解 电力市场化将进入新一轮的发展 这给有关产业和研究 部门带来了可预期的机遇和挑战 电力从生产到使用是瞬间完成的 电网公司在遵循电 网 安全第一 原则的前提下要制订一个电力市场交易规则 按照购电费用最小的经济 目标来运作 设某电网有若干台发电机组和若干条主要线路 每条线路上的有功潮流 输电功率 和方向 取决于电网结构和各发电机组的出力 电网每条线路上的有功潮流的绝对值有 一安全限值 限值还具有一定的相对安全裕度 即在应急情况下潮流绝对值可以超过限 值的百分比的上限 当发生输电阻塞时 需要研究如何制订既安全又经济的调度计划 已知电力市场交易规则 市场交易 调度中心在当前时段内要完成的具体操作过程 和输电阻塞管理原则 试完成如下工作 1 某电网有 8 台发电机组 6 条主要线路 表 1 和表 2 中的方案 0 给出了各机组的当前 出力和各线路上对应的有功潮流值 方案 1 32 给出了围绕方案 0 的一些实验数据 试用这些数据确定各线路上有功潮流关于各发电机组出力的近似表达式 2 设计一种简明 合理的阻塞费用计算规则 除考虑上述电力市场规则外 还需注意 在输电阻塞发生时公平地对待序内容量不能出力的部分和报价高于清算价的序外容 量出力的部分 3 假设下一个时段预报的负荷需求是 982 4MW 表 3 表 4 和表 5 分别给出了各机组 的段容量 段价和爬坡速率的数据 试按照电力市场规则给出下一个时段各机组的 出力分配预案 4 按照表 6 给出的潮流限值 检查得到的出力分配预案是否会引起输电阻塞 并在发 生输电阻塞时 根据安全且经济的原则 调整各机组出力分配方案 并给出与该方 案相应的阻塞费用 5 假设下一个时段预报的负荷需求是 1052 8MW 重复 3 4 的工作 模型假设 模型假设 1 电力从生产到使用的四大环节 发电 输电 配电和用电是瞬间完成的 3 2 每个时段的负荷预报和机组出力分配计划的参照时刻均为该时段结束时刻 3 机组当前出力是对机组在当前时段结束时刻实际出力的预测值 4 假设每台机组单位时间内能增加或减少的出力相同 该出力值称为该机组的爬坡速 率 由于机组爬坡速率的约束 可能导致选取它的某个段容量的部分 5 为了使得各机组计划出力之和等于预报的负荷需求 清算价对应的段容量可能只选 取部分 6 网方和发电商总能就阻塞费用达成协议 从而消除阻塞 保证电网的安全性 7 若 8 台机组总发电量和 6 条主要线路的总潮流有差值 则认为其为其他次要线路上 的潮流之和 那么下一时段预报的负荷需求应该等于机组总发电量 符号说明 符号说明 B 阻塞费用 i W 机组i在机组出力方案调整之前的发电功率 821K i i W 机组i在机组出力方案调整之后的发电功率 821K i i P 机组i在机组出力方案调整之前的报价 821K i i P 机组i在机组出力方案调整之后的报价 821K i Cl 方案清算价 E 各机组出力预案中网方的购电费用 E 机组出力方案调整后网方的购电费用 R 下一时段的负荷需求 ij L 机组i第 j 段的段价 ij V 机组i第 j 段的段容量 i S 机组i的爬坡速率 i W0 机组i当前发电功率 i Lim 线路i的潮流限值 i 线路i的安全裕度 4 i 线路i的潮流超过潮流限值的百分比 又名超限系数 821 WWWTiK 线路i上潮流关于各发电机组出力的近似表达式 ij LF 各机组可用段容量 if TD 线路i的危险系数 即其陷入瘫痪的概率 i TLs 线路i瘫痪给网方造成的损失 模型建立 模型建立 根据题意 对于问题3 4 基本模型建立如下 基本模型基本模型 ij PREMin 4 1 10 2 1 8 2 1 LL ji i MaxMin 6 2 1 L i BMin ts RLF i k j t kt 11 j 1 it t 1 j it t 1 j 1 it t 1 0 V 15 V 15 15V oii ijijoii oii WS F LVWS WS 当 当 其他 10 2 1 8 2 1 LL ji 其他 当 0 G 6 2 1 6 2 1 821 821 i ii ii i i i ii LimWWWT LimWWWT i Lim G i K K L L 其中 目标描述在机组处理分配预案中网方购电费用最小 目标描述每条线路上潮流绝对值超过限值的百分比最小 5 目标描述阻塞费用最小 约束描述的是可用段容量 ij LF从低到高选取直到和不小于预报的负荷 题中等 于预报的负荷是因为它允许选可用段容量的部分 两种约束对模型求解是等价的 ij LF见下表 约束描述各机组的潮流增长不应超过爬坡速率的限制 约束描述各线路潮流值超出潮流限制的百分比不超过安全裕度 并且最小值为 0 阻塞费用B的函数将在对问题2的回答中给出 我们构造各机组可用段容量 ij LF来对爬坡能力进行约束 即当线路的爬坡速率允 许线路取到段容量 ij V或是 ij V的一部分时 ij LF即为线路允许取到的段容量 否则 ij LF为0 以可用段容量 ij LF代替段容量 ij V即可对线路的爬坡能力进行限制 如下 图所示 表1 各机组可用段容量 单位 MW 机组 段 1 2 3 4 5 6 7 8 9 10 1 70 0 50 0 0 30 0 0 0 3 2 30 0 20 8 15 6 2 0 0 7 3 110 0 40 0 30 0 20 28 0 0 4 55 5 10 10 10 9 5 0 0 0 0 5 75 5 15 0 15 15 0 10 10 7 1 6 95 0 10 20 0 15 10 5 0 0 7 50 15 5 15 10 7 1 0 0 0 0 8 70 0 20 0 20 0 7 10 15 5 模型的描述见附录二 多元线性回归模型 多元线性回归模型 回归分析中最简单的形式是xy 10 yx 均为标量 称一元线性回归 1 它 的一个自然推广是x为多元变量 形如 xxxy m L 210 2 m 1 因变量y与自变量 m xx 1 L之间是否存在如 1 所示的线性关系是需要检验的 显然 如果所有的 1 mjL 都很小 y与 m xx 1 L的线性关系就不明显 所以可 令原假设为 1 0 0 mjH j L 2 当 0 H成立时 回归平方和U和残差平方和Q满足 1 1 mnmF mnQ mU F 3 在显著水平 下有 1分位数 1 1 mnmF 若 1 1 设机组j属于过发电补偿情况 此类机组共有m个 则其在低于对应报价的清算价 上出力的序外容量为 ii WW 显然有 ii WW 在 ii WW 这段长度上网方所付单价由报 价 j P 降为清算价Cl 我们无法判断在这段长度上机组j发电是否亏本 为保护发电商 的积极性 我们拟用方案调整后的报价 j P 和清算价Cl的差来作为对在低于对应报价的 清算价上出力的序外容量的补偿单价 则对机组j的过发电补偿即为 4 1 ClPWWB jjjrj jj WWjj 对于方案改变前后发电功率无变化的机组 显然对其补偿为0 则对于全体机组 阻塞费用可表示为 9 m j rj n i si BBB 11 显然此公式形式简洁 并且正如前文所解释的 它还公平地对待了产生阻塞费用的 两种情况 符合此问题的要求 是一种简明 合理的阻塞费用计算规则 规则二规则二 在目前出现的交易模式下 发电方经过网方出售电力 申报发电价格和发电容量 用 户通过网方购买电力 此处不考虑用户投标的情况 于是在不考虑网络约束的条件下 网方采用使购电费用最小的方法确定市场清算价Cl 获得各节点的发电功率分别记为 n PPP L 21 于是 网方的总购电费用为 n i i PClE 1 4 1 由以上安排的发电计划还必须经过网络约束检验 若没有阻塞 则执行发电计划 若出现传输阻塞 则必须在计及网络约束的条件下 采用使购电费用最小的目标函数 重新调度发电出力 通过优化求解 发电功率被重新安排 各节点的发电功率分别记为 n PPP L 21 阻塞消除后 发电方的结算规则如下 多发电量按发电报价结算 少发电量按Cl与发电报价之差结算 于是 支付多发电量的发电方费用为 jj m j m j jjjjr WWjjClWPWWE 4 1 11 支付少发电量的发电方费用为 ii n i i n i iiis WWiiClWPClWWE 4 1 11 因此 阻塞费用为 r s BEEE 规则三 规则三 将阻塞费用在受到影响的nm 家发电商中的合理分配看作n人合作对策 求发电 商的满意解 1 但这类问题要求分配的补偿总额是定值 因此 可考虑首先用规则一 确定某一调整后方案的阻塞费用 则对每个确定的调整方案 阻塞费用B成为定值 即 求阻塞费用时用规则一 而分配阻塞费用时用规则三 i方分配的满意度定义为 ii ii i de dx u 其中 i d 是当前得到的补偿 i e 是理想补偿 i x 是最后得到的补偿 我们认为当前补偿都为0 而理想补偿就是将其补偿到和方案调整之前一样 对m个欠发电机组 ii ii i de dx u 0 i d ClWe ii 10 对n个过发电机组 ii ii i de dx u 0 i d ClWWe iii 为追求各方的满意度都高 用最小最大模型 i i uMinMax ts Bx nm i i 1 得到 n i n i ii nm i i iiii de dB udeudx 11 1 即为各机组得到的补偿 规则四规则四 电站单位时间的运行成本与输出功率的平方成正比 2 2 C PP 即单位时间单位输出功率的成本与输出功率成正比 dc Cost PP dP 其中 Cost P 单位时间单位输出功率的成本 比例系数 P 输出功率 C P 单位时间运行成本 题目所给各机组的段价 表4 是单调不减的 实际表中是单调递增 对各机组各 段价与其对应的内容量进行拟和 我们发现满足很好的线性关系 注 图2 3 4分别为机 组报价分段数最多 10段 最少 4段 及均值 6段 拟合曲线与散点图 10 78 86 3 55 CostPP 11 80100120 150 200 250 300 350 图2 4 517 081 5 25 CostPP 100120140160180200 100 100 200 300 400 500 6 225 09 2 55CostPP 图3 150175200225250275300 100 200 300 400 500 图4 注 在第一段报价 为负值 和最后一段报价由于受市场因素影响 与拟合曲线有较大偏离 故拟和时舍弃 由于全部机组所有出力均按清算价结算 我们可近似认为发电方利润全部由清算 价与序内容量对应报价差带来 那么我们可假设报价即为单位时间单位输出功率的成 本 即对应每机组段价与发电量拟合出的曲线即为 Cost P 基于连续成本函数 我们可以使得发电方调整后利润不减少为赔偿目标 提出规则 四 12 调整规则后付给发电放的总费用 8 1 i iiiiii WWCostClWWCostClMaxE 于是阻塞费用EEB 规则的比较和评价规则的比较和评价 规则一和规则二从不同的方向上刻画了阻塞费用 形式上都比较简明 规则一更为 简明一些 同时都能适当兼顾到公平的原则 规则三形式上略显烦琐 但由于引入了满意度的概念 使得阻塞费用的分配是否公 平可以量化 规则四以发电方调整后利润不减少为赔偿目标 因此也有很好的公平性 同时显得非 常严谨 但由于其需要拟和一些曲线 这些曲线又跟实时报价和段容量有关 因此求解 时显得较为烦琐 综合以上评价 可以得出 规则三是这四个规则中对简明和公平兼顾最好的 所以在这些规则中选择规则三作为问题2的答案 3 对问题对问题3的回答的回答 基本模型中的目标和约束 组成了本问题的模型 模型模型3 1 ij PREMin 4 1 10 2 1 8 2 1 LL ji ts RLF i k j t kt 11 其他 当 当 V15 15 V 15 V 0 1 j 1t it j 1t it 1 j 1t it ioi ioiij ioi ij SW SWV SW LF 10 2 1 8 2 1 LL ji 目标和约束的意义见基本模型 算法描述算法描述 1 按段价由低到高选取各机组的可用段容量或其部分 直到它们之和等于预报的 负荷 2 这时每个机组被选入的可用段容量或其部分之和形成该时段该机组的出力分配 预案 3 最后一个被选入的段价 最高段价 称为该时段的清算价 模型模型3 1的求解的求解 13 用aMathematic软件按照算法描述对此模型进行求解 得到清算价为303元 MWh 下一个时段各机组的出力分配预案如下 表3 负荷需求为MW4 982时各机组分配预案 单位 MW 机组 1 2 3 4 5 6 7 8 方案 150 79 180 99 5 125 140 95 113 9 此外手算也不失为一种便捷方法 手算步骤按照算法描述进行 4 对问题对问题4的回答的回答 基本模型中的目标 和约束 组成了本问题的模型 模型模型4 1 i MaxMin 6 2 1 L i BMin ts RLF i k j t kt 11 其他 当 当 V15 15 V 15 V 0 1 j 1t it j 1t it 1 j 1t it ioi ioiij ioi ij SW SWV SW LF 其他 当 0 G 6 2 1 6 2 1 821 821 i ii ii i i i ii LimWWWT LimWWWT i Lim G i K K L L 目标和约束的意义见基本模型 模型模型4 1的求解的求解 按照输电阻塞的管理原则 首先要尝试调整各机组分配方案使得输电阻塞消除 借 助Lingo软件通过二分法的思想编程对方案进行调整 可找到一个负荷需求的上界 如 果预报的负荷不大于此上界 则必然能够消除输电阻塞 我们算出这个上界为 MW45 983 因此对于MW4 982的负荷需求可以消除输电阻塞 14 因此模型4 1 中的安全性目标 i MaxMin 必然得到满足 于是安全性目标退化 为约束 模型退化为求阻塞费用最小的单目标规划模型 用Lingo软件对此模型进行求 解 由于约束条件和目标条件的关系过于复杂 因此Lingo无法找到最优解 我们认为 当各线路潮流值调整到刚刚没有超过潮流限值时的阻塞费用最小 因此 我们拟用各线路潮流值与潮流限值的比值之和最大来做为目标函数 即用 6 1i i i Lim T Max 代替 BMin 函数具体程序见附录 由此我们得到了一个十分接近最优解的局部最优解 给出调 整后的各机组出力分配方案以及与之对应的阻塞费用如下 表4 负荷需求为MW4 982时调整后各机组分配方案 单位 MW 机组 1 2 3 4 5 6 7 8 方案 153 86 66 228 69 83 152 95 80 91 117 对应阻塞费用 B 4852 69元 5 对问题对问题5的回答的回答 1 将负荷需求MW8 1052代入模型3 1 得到清算价是356元 MWh 下一个时段各 机组的出力分配预案如下 表5 负荷需求为MW8 1052时各机组分配预案 单位 MW 机组 1 2 3 4 5 6 7 8 方案 150 81 218 2 99 5 135 150 102 1 117 2 在安全且经济的原则下讨论有关方案的选取 就本身而言是一个多目标的规划问 题 Max 1 xf 1 xf为安全性描述函数 Max 2 xf 2 xf为经济性描述函数 ts n RXx 由于第三问的数值982 4 可以通过调整处理方案消除阻塞 利用贪婪算法求解 但是对于本问的1052 8 是必然要超出限值的 所以就必选按照求解多目标规划的方法 来做 由于要遵循电网 安全第一 的原则 所以自然的采用分层排序法来求解 得到原 问题的一个有效解 首先将两个目标依其重要程度排序 按照 安全第一 的原则下 排序为安全性 经济性 然后依次求解两个单目标规划即可 我们定义各个线路有限潮流超过限值的值的百分比为超限系数 超限系数的最大值超限系数的最大值 越小 则整个电网的安全性越好 那么安全性描述函数可以刻画为超限系数最大值有关 的函数 在电网传输中 把超限系数降的越低 网方所付出的阻塞费用就越大 换言之 超 15 限系数越小 阻塞费用就越大 可以将经济性描述函数也刻画为与超限系数最大值的 取值有关的函数 经原规划依次求解以下两个单目标规划即可 1 p Max 1 xf ts n RXx 先求得 1 p规划的最优值为 1 p 令 11 pyfxyy 2 p min 2i xf ts yxi 则可得到多目标规划有效解 如上文分析 以及 1 xf 2 xf与超限系数的关系 则可将 1 p规划转化为以下模型 模型假设 need 下一时段的需求 i t 第i条线路的超限系数 i c 第i条线的裕度 n R 满足爬坡要求的所有可行值的集合 1 p Min Max i t ts 8 1j j xneed n i Rx i t i c 81 xxftiK 6 2 1 L i 解得最优目标为 k 然后再来求解模型 2 p 2 p Max 6 1i i t ts ktMax i 8 1j j xneed n i Rx i t i c 81 xxftiK 6 2 1 L i 用Lingo 8编程实现上述思想 求解上述模型 结果如下 16 机组 1 2 3 4 5 6 7 8 方案 153 88 228 99 5 152 113 102 1 117 并且单独运行求解 1 p的程序 得到的结果与上表相同 这说明我们定义的安全性的描述是很中肯的 然后运用分层排序法去求解第四问 与第四问的贪婪算法比较 取超限系数为0 这可由982 4直接取出来 得到一个最优值 从而说明我的安全性的描述是非常有效的 我得这样的规则是简明有效的 从而也 从一个侧面得到网络的完整性 阻塞费用85 2234 B元 6 关于超出裕度 必须拉闸限电的讨论 关于超出裕度 必须拉闸限电的讨论 当拉闸限电必然要发生的时候 考虑次序还是安全第一 但尽可能的满足需求已成 为第二位的目标 积极可能的供电给用户 可以不考虑费用的优化 我们用二分法 找出这个极限值为1094 6MW 模型的评价 模型的评价 优点 1 问题 1 我们提出多元线性回归模型来解决问题 1 改模型剔除了异常的数 据 从而得到相对更优的值 2 问题 2 我们提出 4 种符合题义的阻塞费用计算规则 从中经过比较分析得 到一个兼顾简明 合理 公平的规则 3 对于问题 3 4 5 我们提出的基本模型具有通用性 4 模型扩展中我们提出对双目标模型的改进模型 将其巧妙的化成单目标模型 进行求解 缺点 缺点 由于我们一开始就把问题 3 问题的模型刻画成一个通用模型 这对问题的求解 造成了一定困难 也使得模型描述相对困难复杂 对于问题 4 我们提出了简易模型 把双目标退化成单目标 这种退化准则在不能 调试使其小处电阻塞的情况下并最合理的 在模型扩展中 我们提出了新的退化方式 模型扩展 模型扩展 对问题对问题4中模型中模型4 1的改进的改进 模型4 1为双目标规划 我们可以将其转化为单目标规划进行求解 目标是对电 网的安全性的刻画 为了将目标的量纲应统一到费用上来 我们引入了线路的危险 系数 f D这一概念 即线路陷入瘫痪的概率 17 经过分析 函数 f D应该有如下特点 1 当线路潮流值 i T没超过潮流限值 i Lim时 0 f D 即线路没有瘫痪的可能 2 当线路潮流值 i T超过了安全裕度所允许的最大潮流值 1 ii mLim 时 1 f D 即线路一定会瘫痪 3 当线路潮流值在以上2情况之间时 f D关于潮流 i T应满足是大致下凸的函数 即随着潮流值的增加 危险系数的增长越来越快 故可将 if TD刻画成 1 1 0 2 2 其他 当 当 ii ii iii ii if Lim LimT LimT LimT TD 函数形状大致为 图5 函数形状可以很好的符合以上我们分析其应有的特点 下面量化线路i瘫痪给网方造成的损失 简称为瘫痪损失 i TLs 首先做如下假设 1 若线路i瘫痪 当前时段下经过该线路的电能全部损失 称作电能损失 i TLse 2 网方能够实时监视网络运行情况 若线路i瘫痪 网方能够及时发现并在下一时 段前修好此电路或起用备用电路 3 一条线路由好多变电站连接 这些变电站是串联关系 其中任何一个瘫痪都会 18 导致整个线路的瘫痪 所以线路瘫痪并不会导致整条线路的报废 所以对瘫痪线路的维 修或设备更换费用 i TLsr不应该是一个天文数字 i TLsr与电能损失 i TLse在数量级 上应该是可以比拟的 我们假设 5 ii TLseTLsr 由以上假设和对问题的分析 定义线路i的风险费用 6 ii iiii TDfTLse TLsrTLseTDfTRsk 因此目标求潮流值超出潮流限值的百分比最小可转化为求这6条线路总的风险 费用最小 即为 6 1i i TRskMin 则双目标规划的模型4 1可转化为单目标规划的模型4 2 模型4 2 6 1i BTRskMin i ts RLF i k j t kt 11 其他 当 当 V15 15 V 15 V 0 1 j 1t it j 1t it 1 j 1t it ioi ioiij ioi ij SW SWV SW LF 其他 当 0 G 6 2 1 6 2 1 821 821 i ii ii i i i ii LimWWWT LimWWWT i Lim G i K K L L 6 ii iiii TDfTLse TLsrTLseTDfTRsk 将模型4 1化成单目标规划后 用模拟退火法可以将其解出 19 数据拟合的其他情况数据拟合的其他情况 1 对于问题一 如果严格考虑零输入零输出 我们对不带常数项的情况同样进行 了数据拟和 得到一组新的系数 利用这组系数 我们重新计算了问题4 5 网络的 改变不影响问题2 3的回答 1 当预报负荷为982 4MW时 经过调节各机组出力分配方案 不能够消除电阻塞 按模型5 1求解得到各机组调 整后的分配方案如下 单位MW 机组 1 2 3 4 5 6 7 8 方案 135 04 58 228 95 26 152 95 102 1 117 对应阻塞费用 B 4432 13元 此时潮流 各条线路上潮流值占限制的百分比为 限值 1 2 3 4 5 6 安全裕度 102 6 98 9 97 9 87 3 96 7 102 6 2 当报价负荷为1052 8MW时 经过调节各机组方案 不能削出电阻塞 通过调节 我们发现无论怎样条就不能使 其在安全裕度内工作 所以必须采取拉闸限电 借助Lingo软件 我们用二分法的的方法求出了不拉闸限电的极限值 1043 8MW 对于这个不拉闸限电的极限值 我们按照按对问题3的求解方法求解出清算价为 356 元 MW 分配预案如下 单位 MW 机组 1 2 3 4 5 6 7 8 方案 150 81 209 2 99 5 152 150 102 1 117 按照模型5 1求解 得到各机组调整后的分配方案如下 单位 MW 机组 1 2 3 4 5 6 7 8 方案 153 58 228 99 5 152 134 2 102 1 117 对应阻塞费用 B2111元 此时潮流 各条线路上潮流值占限制的百分比为 限值 1 2 3 4 5 6 安全裕度 112 9 103 6 105 90 108 6 112 5 2 考虑到题目所给大量数据中 不准确数据对结果的影响 我们用Matlab中的回 归分析工具箱得到剔除显著性程度极小的数据的一组系数 同样 利用这组系数 我们重新计算了问题4 5 网络的改变不影响问题2 3 的回答 1 当预报负荷为982 4MW时 经过调节各机组出力分配方案 能够消除电阻塞 按模型4 1求解得到各机组调整 后的分配方案如下 单位 MW 机组 1 2 3 4 5 6 7 8 20 方案 128 76 69 24 176 54 85 03 152 151 74 102 1 117 对应阻塞费用 B 1232 11元 此时潮流 各条线路上潮流值占限制的百分比为 限值 1 2 3 4 5 6 安全裕度 100 100 100 90 2 77 9 100 2 当报价负荷为1052 8MW时 经过调节各机组方案 不能削出电阻塞 按模型5 1求解得到各机组调整后的分配 方案如下 单位 MW 机组 1 2 3 4 5 6 7 8 方案 153 79 73 197 18 96 79 152 155 102 1 117 对应阻塞费用 B 164 1元 此时潮流 各条线路上潮流值占限制的百分比为 限值 1 2 3 4 5 6 安全裕度 112 9 112 9 100 8 85 8 112 9 112 0 通过上述结果我们能够看出 不准确数据对结果的影响是很大的 在做数据拟和时 应用回归分析的方法剔除不准确数据来处理 灵敏度分析 灵敏度分析 1 安全裕度变化的灵敏度分析安全裕度变化的灵敏度分析 各