2004年全国大学生数学建模竞赛B题.pdf
《2004年全国大学生数学建模竞赛B题.pdf》由会员分享,可在线阅读,更多相关《2004年全国大学生数学建模竞赛B题.pdf(12页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、欢迎光临中国数学建模网www.shumo.org 感谢您对网站建设一如既往的支持1 输电阻塞管理的优化调整摘要本问题是一个输电阻塞管理的优化调整问题,即是在发生输电阻塞时,制订既安全又经济的方案对预案进行调整,安排机组的出力分配问题。首先我们用多元线性回归的方法分析了题设给出的实验数据,得到各线路上有功潮流关于各发电机组出力的近似表达式,同时对表达式作了相应的显著性检验,认为所得近似表达式是可信赖的;其次我们用序内容量不出力部分的损失和序外容量出力部分的损失进行加权补偿的办法,给出了阻塞费用的计算规则;然后我们以电力市场交易原则为前提,在 0 方案发生时段为当前时段等一系列合理假设的基础上,以
2、各机组爬坡速率和段容为限制条件,设计了按段价从低到高排序搜索的算法,对给定的下一时段的的负荷需求作出各机组出力分配预案,并利用上述近似表达式求出每条线路上的有功潮流值,将之与线路潮流限值相比,判断该方案是否引起输电阻塞。对发生输电阻塞的出力方案需进行优化调整,使调整后的出力分配方案既安全又经济。根据输电阻塞管理原则,我们建立了一系列的求解模型:可消除输电阻塞的调整方案的非线性规划模型:以阻塞费用最小为目标,以序外容量出力部分和序内容量不出力部分为决策变量在上述模型解空间内没有最优解的情况下,建立使用安全裕度输电的非线性规划模型:以线路上潮流值超过限值百分比之和最小为目标,以序内容量不出力部分和
3、序外容量出力部分为决策变量,同时给出安全与经济约束下的调整模型在模型仍没有最优解的情况下,建立拉闸限电的线性规划模型因此,只要给定下一时段的负荷预报,顺次用这几个模型求解,直至可以求出最优解则认为已经得到调整后的优化出力分配方案。最后我们对上述模型进行了分析并给出了改进方向。对问题三和问题五,对下一时段的预报负荷需求作出各机组出力分配预案如下:负荷机组 1 机组 2 机组 3 机组 4 机组 5 机组 6 机组 7 机组 8 982.4 150 79 180 99.5 125 140 95 113.9 1058.2 150 81 218.2 99.5 135 150 102.1 117 对于不
4、同的阻塞费用补偿权重,可以得到不同的调整方案和不同的阻塞费用,这里我们仅以经济补偿权重1122,为例,得到问题四调整后的出力分配方案,相应的阻塞最小费用为 3997.074 元,具体的调整后的出力分配方案如下:机组1 2 3 4 5 6 7 8 方案150.2613 88 228 90.1786 152 96.86 60.1 117 欢迎光临中国数学建模网www.shumo.org 感谢您对网站建设一如既往的支持2 一、问题重述我国电力市场初期是发电侧电力市场,采取交易与调度一体化的模式。电网公司在组织交易、调度和配送时,必须遵循电网“安全第一”的原则,同时制订一个电力市场交易规则,按照购电费
5、用最小的经济目标来运作。市场交易调度中心根据负荷预报和交易规则制订满足电网安全运行的调度计划各发电机组的出力(发电功率)分配方案;在执行调度计划的过程中,还需实时调度承担AGC(自动发电控制)辅助服务的机组出力,以跟踪电网中实时变化的负荷。电网中每条线路的有功潮流(取决于电网结构和各发电机组的出力)的绝对值有一安全限值,限值还有一定的相对安全裕度(即在应急情况下潮流绝对值可超过限值的百分比上限)。如果某条线路上的有功潮流的绝对值超出限值,称为输电阻塞。当发生输电阻塞时,需要制订既安全又经济的调度计划来安排各机组的出力分配。电力市场交易规则:1.以 15 分钟为一个时段组织交易,各机组在当前时段
6、开始时刻前将可用出力由低到高分成至多10 段给出下一个时段的报价,每个段的长度称为段容量,每个段容量的报价称为段价,段价按段序数单调不减。2.当前时段内,市场交易调度中心根据下一时段的负荷预报,按段价从低到高选取各机组的段容量或其部分,直到它们之和等于预报的负荷,这时每个机组被选入的段容量或其部分之和形成该时段该机组的出力分配预案(初始交易结果)。最后一个被选入的段价(最高段价)称为该时段的清算价,该时段全部机组的所有出力均按清算价结算。市场交易调度中心需要计算当执行各机组出力分配预案时电网各主要线路上的有功潮流,判断是否会出现输电阻塞。如果不出现,接受各机组出力分配预案;否则,按如下原则实施
7、阻塞管理:输电阻塞管理原则:(1)调整各机组出力分配方案使得输电阻塞消除。(2)如果(1)做不到,还可以使用线路的安全裕度输电,以避免拉闸限电(强制减少负荷需求),但要使每条线路上潮流的绝对值超过限值的百分比尽量小。(3)如果无论怎样分配机组出力都无法使每条线路上的潮流绝对值超过限值的百分比小于相对安全裕度,则必须在用电侧拉闸限电。(4)当改变各机组出力分配预案时,一些通过竞价取得发电权的发电容量(称序内容量)不能出力;而一些在竞价中未取得发电权的发电容量(称序外容量)要在低于对应报价的清算价上出力。因此,网方应该为因输电阻塞而不能执行初始交易结果付出代价,网方在结算时应适当地给发电商以经济补
8、偿,由此引起的费用称为阻塞费用。网方在电网安全运行的保证下应考虑尽量减少阻塞费用。我们要做的是:1.根据围绕当前出力分配方案的各机组的当前出力和各线路对应的有功潮流值的实验数据,确定各线路上有功潮流关于各发电机组出力的近似表达式。2.设计一种简明、合理的阻塞费用计算规则,须考虑电力市场规则,还应在输电阻塞时公平地对待序内容量不能出力的部分和报价高于清算价的序外容量出力的部分。3.假设下一时段预报的负荷需求是982.4MW,根据给出的各机组的段容量、段价和爬坡速率的数据,按照电力市场规则给出下一个时段各机组的出力分配预案。4.按潮流限值,检查得到的出力分配预案是否引起输电阻塞,在发生输电阻塞时,
9、根据安全且经济的原则,调整各机组出力分配方案,并给出与该方案相应的阻塞费用。5.假设下一时段预报的负荷需求是1052.8MW,重复 34的工作。欢迎光临中国数学建模网www.shumo.org 感谢您对网站建设一如既往的支持3 二、问题分析本题是一个输电阻塞管理的优化调度问题,目的是在根据电力市场交易规则制订的各机组出力分配预案引起输电阻塞时,对该分配预案进行优化调整,使得因调整引起的阻塞费用在电网安全运行的保证下尽量的少,即实现购电总费用最小的经济目标。在电力市场交易规则和市场交易-调度中心具体操作下,我们要解决的是如何调整各机组出力分配方案,使其在市场交易过程中既安全又经济。问题一:根据给
10、出的32 组数据,求出各线路的有功潮流关于各机组当前出力的近似表达式,每一条线路的有功潮流值是关于8 个发电机组出力的多元函数,这样我们首先考虑其是一个多元的线性关系,于是我们想到用多元线性回归分析来寻求变量间关系的近似表达式,这样可以根据多元线性回归给出有功潮流关于各机组出力的近似表达式,然后对其检验看是否可以信赖。问题二:给出一个简明,合理的计算阻塞费用的规则,不只是要考虑电力市场规则,还要注意序内不出力和序外高于清算价的出力。这样我们可以对序内容量不出力部分和序外容量出力部分的损失给以相应的补偿,因为它有很主观因素在里面,所以我们借用层次分析的思想,给两者的损失前均加个权重。这样可以根据
11、具体情况,确定权重的大小。给出阻塞费用的计算公式。问题三:0 方案作为当前方案时,要求给出下一时段预报负荷需求为982.4MW 的各机组的分配预案。我们就要根据电力市场规则设计一个算法,用计算机编程来解决这个问题。问题四:利用问题三的结果,把得到的各个机组的分配预案代到第一问得到的各线路潮流值关于各机组出力的表达式,求出其有功潮流值,与各个线路的潮流限值比较,判断是否发生阻塞,若发生阻塞则就要按市场调度中心给出的输电阻塞管理规则进行方案调整。首先要考虑调整各个机组的分配预案,从而消除阻塞,但在消除阻塞的过程中要考虑相应的阻塞费用,应使其尽量小。若不能消除阻塞,则考虑使用安全裕度输电,但这时要求
12、使每条线路的潮流的绝对值超过限值的百分比尽量小。如果用以上两种方案都不能解决输电阻塞,那只有采取拉闸限电,但要使拉闸限电消除阻塞时所限电量最小,这时的阻塞费用最小。我们可以根据上面调整顺序建立相应的处理模型。问题五:表面上看它是第三和第四问的重复,但又可能有不同的意义,即三、四问是一种阻塞情况,而第五问又是一种阻塞情况,所以我们要把第四问的三个处理各种情况阻塞的模型建立完善。三、符号约定TT为一个时段的时间R下一时段预报的负荷需求if第i条线路的有功潮流jx第j机组的出力jv第j机组爬坡速率欢迎光临中国数学建模网www.shumo.org 感谢您对网站建设一如既往的支持4 id第i条线路有功潮
13、流的限值iw第i线路的相对安全裕度jkN第j机组第k段的段容量jkQ第j机组第k段报价jD第j机组在出力分配预案中的出力P出力分配预案的清算价0 jk预案中第j机组最后被选入的段序数0 jM第j机组当前出力jB第j机组序内不能出力的部分jkC第j机组第k段序外出力(0 jkk)S阻塞费用jM发生发电阻塞时,调整出力分配方案后的第j机组出力iH第i线路上有功潮流绝对值超过限值的百分比1z,2z取 0 或 1 约定:i表示线路序号,i=1,2,,,6 j表示机组序号,j=1,2,,,8 k表示机组段序号,k=1,2,,,10 四、模型假设对本模型的假设如下:每个时段的负荷预报和机组出力分配计划的参
14、照时刻均为该时段结束时刻 当前时段即是以各机组以0 方案出力的时段 机组当前出力是对机组在当前时段结束时刻实际出力的预测值 每台机组单位时间内能增加或减少的出力相同,该出力值称为该机组的爬坡速率 对于所有机组,若两段的段价相同,段容量大的优先选择,段容量相同的按机组序数选择欢迎光临中国数学建模网www.shumo.org 感谢您对网站建设一如既往的支持5 五、模型建立与求解51各线路上有功潮流关于各机组出力的近似表达式综合以上分析,根据围绕当前方案0 的各机组出力方案以及线路上对应的有功潮流值的 32 组实验数据,由于每条线路上的有功潮流值是由各机组的出力分配也即多个变量决定的,这种利用实验数
15、据来寻求变量间关系的近似表达式的问题,我们很自然的想到测试技术中解决多个变量的回归问题,于是我们利用Matlab 中提供的多元线性回归函数 regress()对之求解,得到各线路上有功潮流值关于各机组出力的近似表达式如下:11234567821234567831=1 10.2 9 65+0.0 8 2 8+0.0 48 3+0.0 53 0+0.1 1 99-0.0 2 54+0.1 22 0+0.1 21 6-0.0 0 12=1 3 1.2 28 9-0.0 546+0.1 27 9-0.0 000+0.0 3 33+0.0 86 8-0.1 12 4-0.0 1 89+0.0 9 8 7
16、=-1 0 8.8 7 3 2-0.0 69 5+0.0fxxxxxxxxfxxxxxxxxfx23456784123456785123616-0.1 5 66-0.0 0 9 9+0.1 24 5+0.0 021-0.0 0 25-0.2 0 1 4=7 7.4 81 7-0.0 345-0.1 02 4+0.2 052-0.0 2 08-0.0 11 8+0.0 060+0.1 4 49+0.0 7 6 5=1 3 2.9 7 4 5+0.0 00 5+0.2 433-0.0 6 46-0.0 4 11xxxxxxxfxxxxxxxxfxxx45678612345678-0.0 652+0
17、.0 70 3-0.0 0 43-0.0 0 8 9=1 2 0.6 63 3+0.2 378-0.0 60 2-0.0 7 79+0.0 9 30+0.0 469+0.0 00 1+0.1 659+0.0 0 07xxxxxfxxxxxxxx利用 MATLAB 可以得到六个方程的相关系数,回归方程的F统计量,拒绝无效假设的概率,将其值分别列表如下:回归方程相关系数回归方程统计量拒绝无效假设概率1f0.9995 5861.5 0 2f0.9996 7228.7 0 3f0.9999 22351.7 0 4f0.9999 25582.5 0 5f0.9996 6971.8 0 6f0.9998
18、17454.5 0 另外我们利用所得到的方程对方案0 的有功潮流进行了验证,并给出了相对误差,其结果如下:线路1 2 3 4 5 6 0 方案的有功潮流值164.78 140.87-144.25 119.09 135.44 157.69 表达式求出的潮流值164.7192 140.8278-144.20 119.050 135.367 157.615 相对误差(10-4)3.6922 2.995 3.4558 3.368 5.3773 4.7441 同时我们给出了六个方程的残差图形如下图所示,显而易见,这个回归方程是极显著的。从残差图形可以看出,各个数据的残差离零点均较近,且残差的置信区间均包
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2004 全国大学生 数学 建模 竞赛
限制150内