多目标优化模型v3.docxVIP

1、简介 此为水资源的多目标优化配置模型。 水流沿箭头方向流动，从水库1，一直到生态断面2。其中，riv2和riv4为支流入流，其余河道为干流，模型供水区有5个，每个供水区都通过引水渠（绿色channel）和退水渠（红色channel）与河道相连。汇流点起到平衡、传承流量的作用，同时能保证水不会倒流。例如有引、退水渠连接的河道，如果不设置汇流点，很可能退的水会被其相应的引水渠引走，就有问题了。 2、约束 模型中最重要的约束就是水量平衡约束。 对于水库res，需要考虑其蓄水的变化（自身变化量），流入量+蓄水变化量=流出量+损失量+供水量； 对于河道riv，流入量=流出量+损失量+供水量； 对于渠道cha，流入量=流出量+损失量； 对于供水区dis，进入的水为供水，其中一部分全部消耗（耗水），一部分通过退水渠排出（退水）。耗水=供水*耗水系数。流入量=耗水量+退水量 对于汇流点con，流入量=流出量； 另外，变量还有上下限约束。 具体见5、约束的数学表达。 3、目标 模型目标是一方面尽量保证社会经济供水，即给供水区的供水与其需水差的最小（F1），另一方面尽量保证不损害生态，即满足生态断面的水流特征要求（F2）。具体见6、目标函数。 4、变量 模型中的变量为res1的出流、cha1的入流、cha3的入流、cha5的入流、cha7的入流、cha9的入流，其余量可利用以上六个变量求出。因为模型与时间有关（30年连续系列），变量每天都发生变化，因此每天都有六个变量，模型变量总数大致为6*365*30（其实还要考虑闰年的，变量数会比这个多一点点）。 5、约束的数学表达 （1）所有的Qin（t）为t时段入流[in]，Qout（t）为t时段出流[out]，L（t）为t时段损失（蒸发+渗漏）[loss]，Ws（t）为t时段供水[water supply]。 （2）t：以日为步长，共有30年 （3）已知连续30年t时段的水库1入流Qres1_in(t)，河道2和4的入流Qriv2_in(t)、Qriv4_in(t)，5个供水区各自的需水量Qdis_dem(t) 与 耗水系数d（t），每个单元每时段的损失系数，变量要求的上下限。 （4）所有变量均为正值 （5）黄底标出的参数为已知值 水库RES Qres_in(t) + V(t-1) =Qres_out(t) + V(t) + Lres(t) + Ws_res(t) Lres(t) = ((V(t-1) + V(t))/2 –Vmin)*a(t) Vmin≤V(t)≤Vmax Qres_out(t) ≤Qmax 7月-10月（汛期） V(t) ≤ Vlimit+△V1+△V2+△V3 △V1 ≤ Vuse – Vlimit △V2 ≤ Vflood – Vuse △V3 ≤ Vmax – Vflood 1-6、11-12月（非汛期） V(t)≤Vuse RES1 Qres1_in(t) Ws_res1(t) = 0 河道RIV Qriv_in(t) = Qriv_out(t) + Lriv(t )+ Ws_riv(t) Qriv_min≤Qriv_in(t)≤Qriv_max Qriv_min≤Qriv_out(t)≤Qriv_max RIV1 Qriv1_in(t) = Qres_out(t) Ws_riv1(t) = 0 RIV2 Qriv2_in(t) Ws_riv2(t) = Qcha1_in(t) RIV3 Qriv3_in(t) = Qcon1_out(t) Ws_riv3(t) = 0 RIV4 Qriv4_in(t) Ws_riv4(t) = Qcha3_in(t) RIV5 Qriv5_in(t) = Qcon2_out(t) Ws_riv5(t) = Qcha5_in(t) RIV6 Qriv6_in(t) = Qcon3_out(t) Ws_riv6(t) = Qcha7_in(t)+ Qcha9_in(t) 汇流点CON Qcon_in(t) = Qcon_out(t) CON1 Qcon1_in(t) = Qriv1_out(t) + Qriv2_out(t) + Qcha2_out(t) CON2 Qcon2_in(t) = Qriv3_out(t) + Qriv4_out(t) + Qcha4_out(t) CON3 Qcon3_in(t) = Qriv5_out(t) + Qcha6_out(t) CON4 Qcon4_in(t) = Qriv6_out(t) + Qcha8_out(t) + Qcha10_out(t) 渠道CHA Qcha_in(t) = Qcha_out(t) + Lcha(t) Lcha(t)