数模基于粒子群算法的水库群联合优化调度问题_数学建模竞赛论文(36页).docx
《数模基于粒子群算法的水库群联合优化调度问题_数学建模竞赛论文(36页).docx》由会员分享,可在线阅读,更多相关《数模基于粒子群算法的水库群联合优化调度问题_数学建模竞赛论文(36页).docx(35页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、-数模基于粒子群算法的水库群联合优化调度问题_数学建模竞赛论文-第 33 页数学建模校内竞赛论文论文题目:基于粒子群算法的水库群联合优化调度问题2013-07-18基于粒子群算法的水库群联合优化调度问题摘要由于水库在水资源的优化配置中发挥着重要的作用,本文基于粒子群算法建立了水库群联合调度的优化模型,使得它们在保证安全运行的基础上发挥最大的“群体”效应。针对问题一,由于本题是关于各级枢纽电站的统一优化调度问题,具有非线性、离散性等特点。本文采用基于群体智能的启发式全局搜索算法粒子群(Particle Swarm Optimization)算法,通过更新粒子的个体最优解和全局最优逐渐向全局最优靠
2、近。而在枢纽电站的寻求全局最优解的时候,容易陷入局部最优(pBest)的陷阱当中,而本文运用的PSO算法能够有很多措施可以避免它,而找到全局最优(gBest)本文在构造库容水位的关系时,运用SPSS进行了回归分析,得出了关于水位的回归函数。在本文构建水库联合调度优化模型,以嘉陵江上的利泽,渭沱,草街和井口作为水库阶梯分析其在各个约束条件下的发电量最大值优化。发现在仅考虑本段优化的条件下,得到一年之内的最大发电电量为35.63亿千瓦时,相对常年平均发电电量29.016亿千瓦时增长近21.8%。得到了最优化调度的水位和流量关系。针对问题二,运用第一题中的模型,考虑上游枢纽对第一级阶梯的入库流量的影
3、响,考察枯水期时,分析得出上游枢纽电站对于枢纽工程的影响程度。得到最优化调度的流量水位关系。针对问题三,考虑梯级水库的正常蓄水位的提高和水轮机效率对联合调度的影响,运用第一题中得到的模型,考虑正常蓄水位提高0.5米,水轮机效率提高10%的工作效率情况下的最优化水位流量调度关系,分析每个水电站对于正常蓄水位和水轮机效率的灵敏度,并且分析得出其中的影响因素。关键词:粒子群算法(PSO) SPSS回归分析 水库群联合调度优化 一 问题重述随着水资源的不断开发利用,往往在一条河流上或一个流域内建成一批水库,形成了一个水库群。原有的单库分散调度的方式不利于流域内水利综合效益的发挥。水库群的形成,改变了原
4、来单库或少库的水力条件,各水库之间存在相互影响,这就需要站在全流域的高度,采取联合调度的方式,开展水库群优化调度,让它们在保证安全的基础上发挥最大的“群体”效益。重庆市已建成投产的水电站:涪江渭沱航电枢纽工程、涪江富金坝航电枢纽工程、嘉陵江草街航电枢纽工程;以及在建的水电站: 嘉陵江利泽航运枢纽工程、嘉陵江井口航电枢纽工程等。综合考虑水资源利用,在站内优化调度的基础上实施嘉陵江(重庆段)水库(发电)联合优化调度,将草街枢纽及嘉陵江上下游梯级电站纳入统一优化调度,根据来水水情和负荷情况,考虑航运用水要求,优化协调各级枢纽电站问的负荷分配,合理化控制水位,尽可能使流域电站在最优流量区运行,有效解决
5、不同来水条件下的航运和发电等问题,最大限度节约一次能源,优化航运和发电效益。现需要解决以下问题:1、综合考虑各种约束条件,以总出力最大为目标函数建立井口、草街、渭沱、利泽航电枢纽工程联合优化调度模型,设计并实现你的算法,参考附件数据计算最大出力和年发电量,并同实际发电量进行比较。2、利用模型定量分析富金坝和在利泽上游修建大型枢纽工程对联合调度的影响。3、利用模型定量分析提高正常蓄水位和水轮机功率对联合调度的影响。二 模型的假设1、假设水头损失忽略不计,即净水头等于毛水头。2、假设各水库一月份的初始水位为正常蓄水位。3、假设下游水库的入库流量为上游直接汇入该水库的各水库出库流量之和。4、假设水电
6、站的引用流量为其出库流量。三 符号说明符号符号的意义N水电站的出力水轮机的额定流量水电站的入库流量水电站的出库流量H毛水头X水库水位V水库库容时间间隔水轮发电机组的总效率W水电站发电量E发电量提高率k水库编号四 模型的建立与求解4.1针对问题一4.1.1问题一的分析题目要求以总出力最大为目标函数,考虑各种约束条件建立对四个梯级水电站的联合优化调度模型,这是一个具有非线性、离散性特点的复杂优化模型,不适合用动态规划、逐步优化等经典的算法。因此本文采用模拟进化算法中的粒子群(PSO)算法进行求解,它利用群体优势,在没有集中控制、不提供全局模型的前提下来寻找复杂问题的解决方案,具有健壮性、灵活性、高
7、效性等优点。因此,本题的求解思路是:建立总出力最大的目标函数设定各类约束条件应用粒子群算法求解编写matlab程序并得出结果对模型进行评价分析4.1.2 建立目标函数水电站的出力公式为,其中Q为引用流量;H为净水头,根据假设不考虑水头损失,则H为毛水头;为水轮发电机组的总效率,一般为0.700.90,在本题中设定为0.8。考虑到引用流量Q受到水轮机额定流量和出库流量的限制,因此: (4-1-1)则水库群总出力为: (4-1-2)其中k=1,2,3,4分别指代利泽水电站,渭沱水电站,草街水电站和井口水电站。建立的目标函数为: (4-1-3)本文以水位x为决策变量。4.1.3 设定约束条件根据相关
8、的水利资料,一般水电站应满足以下几个约束条件:(1)水库水量平衡约束: (4-1-4)其中,代表k号水库在t时刻的库容;为入库流量;为出库流量。首先,为了得到库容V与水位x之间的关系,本文利用SPSS软件对二者进行拟合,以井口水电站为例:由附件材料的数据表4.1:表4.1井口水电站水位库容关系水位(米)库容(万立方米)17511349176.0412095177.5149521802000019044930196.76000020090037204.7120000210152093213.4180000215.7200000220233023图4-1 井口水电站库容-水位拟合图由图4-1可知,
9、当拟合是该二次函数关系时,能够非常好地反映二者的关系。表 4.2 井口水电站库容-水位拟合程度表因变量:井口水电站库容方程模型汇总R 方Fdf1df2Sig.二次.9981866.31129.000自变量为 井口水电站水位表 4.3 井口水电站库容-水位拟合参数估计因变量:井口水电站库容方程参数估计值常数b1b2二次2916184.769-33754.54198.060自变量为 井口水位。因此,井口水电站库容与水位的函数关系式为: (4-1-5)取为1个月的时间,t为月份,则井口水库水量平衡条件为:(4-1-6)同理可得其他三个水电站的库容-水位关系式。利泽水电站的拟合模型参数如表4.4所示:
10、表 4.4 利泽水电站库容-水位拟合参数估计因变量:利泽库容方程参数估计值常数b1b2二次1636753.449-16390.36541.026自变量为 利泽水位。所以,利泽水电站库容与水位的函数关系式为: (4-1-7)则其水库水量平衡条件为:(4-1-8)渭沱水电站的拟合模型参数如表4.5 所示:表 4.5 渭沱水电站库容-水位拟合参数估计因变量:渭沱库容 方程参数估计值常数b1b2二次-54398.56101.327 自变量:渭沱水位所以,渭沱水电站库容与水位的函数关系式为: (4-1-9)则其水库水量平衡条件为: (4-1-10)草街水电站的拟合模型参数如表4.6所示:表 4.6 草街
11、水电站库容-水位拟合参数估计因变量:草街库容方程参数估计值常数b1b2b3三次2910277.368-24215.780.000.249自变量:草街水位。所以,草街水电站库容与水位的函数关系式为: (4-1-11)则其水库水量平衡条件为:(4-1-12)(2)水库蓄水水位约束: (4-1-13)其中,代表k号水库在t时刻的水位,单位为米(m);为水库死水位,用于保证航运以及正常用水;为正常蓄水位。 具体约束如下所示: (4-1-14)(3)出库流量约束: (4-1-15)其中,代表k号水库在t时刻的出库流量;最小的出库流量在实际中可为0;出库流量的最大值根据每个水电站设计的运行方式决定,具体约
12、束如下所示: (4-1-16)(4)电站出力约束: (4-1-17)其中, 代表k号水库在t时刻的电站出力大小,单位为兆瓦(Mw);由于本题的目标函数是总出力最大,所以只需设定上限,具体约束如下所示: (4-1-18)4.1.3 粒子群算法及其matlab工具箱简介(1)粒子群算法(PSO)简介粒子群优化算法是一类基于群体的演化算法,基本思想是通过群体中个体之间的协作和信息共享来寻找最优解。应用 PSO 算法求解优化问题时,问题的解对应于搜索空间中 “粒子”(particle),每个粒子都有自己的位置和速度(决定飞行的方向和距离),还有一个由优化问题对应函数决定的适应度值(fitness)。粒
13、子寻优的过程就是粒子记忆、追随当前最优粒子,在解空间不断搜索的过程。粒子的每次迭代过程不是完全随机的,如果找到较好的解,就将其保存下来,并以此为依据来寻找下一个解。PSO 求解优化问题时,首先要初始化生成一群随机粒子(对应于问题的一组随机解),在每一次迭代过程中,粒子通过跟踪两个“极值”来更新自己。一个是到目前为止粒子本身所找到的最优解,即个体极值( pbest),另一个是到目前为止整个种群找到的最优解,即全局极值(gbest)。找到这个解后,粒子根据式(4-1-19)和式(4-1-20)来更新自己的速度和位置。设粒子的速度信息、位置信息等可以用一个N 维向量表示。位置向量,速度向量,粒子目前
14、最优位置(对应于问题的局部最优解) ,粒子目前最优位置(对应于问题的全局最优解) 。则速度和位置更新公式如下:(4-1-19) (4-1-20)式中,w为惯性权值;、为粒子在第 k +1次迭代中第 n 维数的速度和更新后的位置;c1和c2为加速系数,分别调节向全局最好粒子和个体最好粒子方向飞行的最大步长,若太小,则粒子可能远离目标区域,若过大可能飞过目标区域,合适的c1,c2可加速收敛且不宜陷入局部最优,通常取;,是(0,1)之间的随机数。为了防止粒子远离搜索空间,粒子的每一维速度都被限制在。(2)粒子群算法求解最优化问题步骤根据上述的PSO的基本原理,PSO求解优化问题的步骤如下:Step1
15、: 初始化粒子种群,设,根据具体研究问题,随机生成L个粒子,其中,、分别为第l个粒子的位置向量和速度向量,;Stept2: 计算每个粒子当前的适应度;Stept3: 更新局部最佳个体和全局最佳个体。找出目前为止,每个个体的最佳适应度值,记为,并记其相应的位置向量为,;找出目前为止,群体的最佳适应度值记为,并记其相应的位置向量为,;Stept4: 更新个体的位置向量和速度向量,产生下一代粒子群。对所有个体中的每一个速度元素和位置元素,采用下式进行更新:(4-1-21) (4-1-22) Stept5: 检验是否满足算法迭代终止条件,若满足,则停止迭代,输出最优结果;否则令,转向b。迭代终止条件(
16、为预先设置的最大迭代次数);或则 (为事先给定的精度)(4-1-23)(3)Matlab粒子群算法工具箱MATLAB粒子群算法工具箱(Psot),是把粒子群算法的核心部分封装起来,只将算法的可调参数提供给用户。该工具箱结构简单,易于掌握,是解决数值优化问题的有力工具。工具箱中包含pso_Trelea_vectorized函数,该函数实现了整个粒子群的初始化、个体最优值计算及更新、全局最有值计算及更新、个体速度更新、个体位置更新等功能。在实际的计算过程中,只需要编写好目标函数并调用pso_Trelea_vectorized.m函数,即可实现使用粒子群算法来寻优的过程。4.1.4 粒子群算法在水库
17、群联合优化调度中的应用根据式(4-1-2)和(4-1-3),可以列写目标函数(也称为适应度函数)的M文件,其中x1、x2、x3、x4分别代表利泽、渭沱、草街、井口的水位,y1,y2,y3,y4分别代表它们的出库流量,文件命名为Max_N.M,具体程序代码见附件(1)。再根据式(4-1-4)(4-1-18),列写包含约束条件的主程序,命名为Find_Max。在主程序中调用pso_Trelea_vectorized.m函数时,根据变量个数以及求解需要,本文设定维数为n=8,粒子群规模为600,最大迭代次数为500次,学习因子C1=C2=2,由于粒子数目较大,因此惯性权重w设定随迭代次数的增加从1.
18、9线性降至0.5。具体的主程序代码和pso_Trelea_vectorized.m函数代码分别见附件(2)和附件(3),最后的运行结果如下所示:以一月份为例,水库群联合优化调度的求解结果为:N1 =1.895853822625391e+004N2 =2.849020775056230e+003N3 =8.946814799251813e+004N4 =5.806147015470056e+004PSO: 251/600 iterations, GBest = 169337.17714852883.ans = 1.0e+005 * 0.00212059585505 0.0020516517333
19、0 0.00202469029754 0.00177451808131 0.00230936244650 0.00124735496805 0.00414264474790 0.00834451898195 1.69337177148529从图(4-2)中也可以看到,PSOt能快速找到满足约束条件下目标函数的最优值,搜索速度和收敛速度都很快,效果十分理想。图4-2 最优个体适应度值程序的运行结果即如表4.7所示:表4.7 一月份各水库的优化调度模型水库编号水库名称入库流量(m3/s)初始水位(m)出库流量(m3/s)末水位(m)出力(Mw)发电量(亿kwh)1利泽230212.3230.942
20、12.0618.960.1372渭沱123206124.74205.172.850.0213草街429.1203414.26202.4789.470.6444井口452177.5834.45177.4558.060.418合计169.341.220需要指出的是,由于一月份的初始水位未知并且难以获得,所以本文将其设定为正常蓄水位,各水库的末水位和出库流量由程序的求解结果ans可得;各水库的出力由运行结果N值得到;而上表的发电量为一个月时间的发电量,即: (4-1-24)注意到出库流量大的水电站其出力不一定大,这是因为出力还与水头有关,同理,可得各水库在每个月份的优化调度模型的最优解:表4.8 利
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数模 基于 粒子 算法 水库 联合 优化 调度 问题 数学 建模 竞赛 论文 36
限制150内