多目标电力系统环境经济调度问题牛奔，王红摘要：随着我国经济的飞速发展，电力工业成为支持我国经济发展的基础工业。
随着电力工业的快速发展，废气、废水的排放等环境污染问题引起各界的广泛关注。
近年来，许多国家限制了火电厂对有害气体的排放量，因此，同时考虑经济因素和污染排放量的多目标电力系统环境经济调度问题，就成为了电力工业十分关注的优化问题。
本案例深入展示了电力系统环境经济调度的多目标问题，适用于运筹学中多目标规划、非线性规划、启发式算法等模块的教学。
关键词：电力系统环境经济调度问题；多目标规划；非线性规划；启发式算法1多目标电力系统环境经济调度问题的提出电力工业是能源工业的重要组成部分，是推动人类文明及支撑社会经济发展的重要基础。
近年来，随着中国经济的持续快速发展，对电力的需求十分强劲。
为了有效缓解电力供需矛盾，国家加快了电力建设步伐。
电力项目建设不仅有力地缓解了各地电力供应紧张的局面，而且对电力工业结构调整与合理布局发挥了重要作用。
我国的主要发电方式为火力发电，这种方式以煤炭消耗为主。
但是，发电用煤的平均灰份高达28%左右，基本上是没有经过洗选的动力煤，外加污染控制和治理技术落后，致使火力发电行业成为二氧化硫、氮氧化物、烟尘等大气污染物的主要排放源，同时也是废水、粉煤灰和炉渣等固体废弃物的主要排放源。
近年来，电力行业的环境污染问题受到广泛关注，许多国家制定了限制火电厂有害气体排放的法规。
火力发电行业控制污染气体、液体、固体排放量的压力日趋上升。
因此，在保证可靠供电的前提下，如何以最低的成本和最少的污染使电力系统正常运行，即电力系统环境经济调度优化，这个多目标优化问题成为电力行业至关重要的优化问题。
2 IEEE-30总线测试系统IEEE-30总线的电力系统有6个发电机，41条线，其单线结构如图1所示。
这是一个标准的测试系统，调度的目的是使得经济成本最低，同时环境污染最小，因此这是一个多目标优化问题。
发电机的燃料消耗成本、固定损耗率及氮氧化物排放量相关数据如表1和表2所示。
表3和表4是系统相关的详细数据。
图1. IEEE-30节点电力系统的单线图表2. 发电机氮氧化合物排放量注：线流量（$lP ）是标准值的110%。
3多目标电力系统环境经济调度问题的数学模型多目标环境经济调度问题即为解一个带有多个等式约束和不等式约束的多个非线性目标函数的最小值，且这些目标函数之间有着相互制约的关系（燃料费用和污染排放）。
目标函数在考虑环境经济调度的情况下，多目标经济调度可以采用以下两个目标函数： （1）电力系统发电燃料总耗量或发电燃料总费用，可以用发电机有功出力的二阶多项式表示：)/($21h P c P b a F Gi i Ni Gi i i G ++=∑= （1）其中，G F 为发电燃料总费用，Gi P 为系统内第i 台发电机的有功出力，N 为系统内发电机组的数目；i i b a ,和i c 分别表示第i 台发电机组耗量特性的常数项、一次项系数和二次项系数。
（2）最小化污染排放量。
考虑到环境污染对生态平衡的影响（如酸雨及臭氧层的破坏），一些法律规定各电厂必须控制氮氧化物和硫氧化物的排放量，以减小空气污染。
另外，有些法规对热辐射也有限制，为了不失一般性，下面仅给出考虑氮氧化物排放限制的情况。
氮氧化物的排放量分别表示为：)/]()(10[122MWh ton e x P c P b a EGi i xP l i Ni Gi i Gi i i GNO +++=∑=- （2）)/]()(10[,2,12,,,2MWh ton ex P c P b a EGii P l i Ni Gii Gi i i GSO +++=∑=- （3）其中，i i i i i l x c b a ,,,,是表示发电机组i 的x NO 或2SO 费用系数。
约束条件（1）功率平衡约束：这是一个等式约束，系统发电机总出力必须满足系统总负荷与传输线路网损之和：01=--∑=LOSS D Ni GiP P P（4）其中，Gi P 为发电机i 的出力；D p 为系统总负荷；loss P 为系统网损。
网络损耗可表示为：∑∑===N i Nj Gjij Gi loss P B P P 11（5）其中，ij B 为网络损耗系数。
（2）机组发电容量约束：这是一个不等式约束，发电机输出功率必须维持在系统稳定运行要求的范围之内：max min i i i P P P ≤≤ （6）本案例所研究的电力经济调度问题的目标函数及约束条件，可描述为一个带有等式与不等式约束的非线性多目标优化问题，其数学表述如下：{}G G E F ,min （7）()()⎭⎬⎫⎩⎨⎧≤=00..x g x h t s （8）上式中，g 与h 分别为上节提及的不等式与等式约束。
最优折衷解多目标问题的解不是唯一的，而是一组Pareto 解集，其中的每个解都是满足条件的。
然而，在实际运行中，调度人员必须从该解集中做出最优选择，最终选择的解便称之为“最优折衷解”。
这里，我们引用了模糊隶属度函数来表示每个Pareto 解中各个目标函数对应的满意度，定义模糊隶属度函数如下：⎪⎪⎩⎪⎪⎨⎧≥--≤=max min max maxmin 0,, 1,i i ii ii i i f f f f f f f f f u i m a xm i n ii i f f f << (9) 其中当0=if u 时表示对某个目标函数值完全不满意，而当1=i f u 时则表示对某个目标函数值完全满意。
对于Pareto 集中的每个解，应用下式求解其标准化满意度值：/,111N N objobj M k k ki i i i k μμμ∑∑∑==== （10） 其中M 为Pareto 集中解的个数，obj N 为待优化目标函数的个数。
最优折衷解即为具有最大标准化满意度值kμ的解。
4基于多目标菌群算法的环境经济调度问题求解多目标菌群算法细菌有一个特别的感知、行动和决策的机制，每个细菌在移动的过程中，当它发现有利环境时，会释放一种引诱剂，以便于其它的细菌能朝它移动的方向移动。
当它发现不利环境时，会释放一种忌避剂，以提醒其它细菌远离。
一般来说，多目标优化算法的目标值通过反映个体健康状况的适应生存机制获得。
考虑函数1212(,,...)(,,...,)n n J J J J f f f ==，多目标菌群优化模拟这种个体间相互吸引和排斥的社会行为，这种吸引排斥的行为用下式表示： 1(,(,,))(,(,,))Si i cc cci J P j k l Jj k l θθθ==∑∑∑∑∑====--+---=Si Si Pm i mm repelent repelent Pm i i mm attract attract w h w d 11121)])(exp([)])(exp([θθθθ（11）其中，(,((,,))cc J p j k l θ是成本函数值，把它加在实际成本函数中，对其进行最小化处理，它表示随时间变化的成本函数。
S 是细菌的总数，p 代表在每个细菌中需要优化的参数数量，attract d , attract w , repelent h , repelent w 是可供选择的不同的系数。
（1）健康排序方法对于有两个目标的MBFO ，需要计算每一个目标函数的值t (1, 2): (,1,,)t J i j k l +，再令)),,1(),,,1((),,1,(),,1,(l k j P l k j J l k j i J l k j i J i cc t t ++++=+θ(12)对于给定的k 和l , 每一个细菌个体i （i =1,…,S ）的健康值是其趋药过程中所有目标函数值的总和：1(,,,)cN it health t j J J i j k l ==∑(13)i health J 用来计算第i 个细菌在它的生命周期内获得营养和躲避有害物的成本。
i health J 值越高，表示第i 个细菌的健康值越低。
按照i health J 的值从小到大排序，i health J 值最高的细菌死亡的可能性最大，而i health J 值最小的细菌将会进行繁殖。
（2）Pareto 占优机制假设第i 个细菌的健康排序位于S 个细菌的前50%（具有较低的健康成本i health J ）。
按照i health J 值大小顺序，如果第j 个细菌受第i 个细菌的支配，那么第j 个细菌将会消亡。
为了保持细菌个体（解）的多样性，因这种机制而死亡的细菌的数量不能超过群体数量的10%。
为了保持群体的规模不变，选择较好的细菌进行繁殖，这些繁殖出的新个体跟它们的父代具有相同的属性（它们所在位置跟父代的位置一致）。
在计算、排序、分离和繁殖的整个过程中，产生一个具有更好解的菌落。
细菌在生命周期内可以直行、翻转或是两者交替进行。
细菌在某个方向()C i 上翻转的步长及其在该方向基础上调整的角度()i 共同决定了细菌的位置。
在这里调整函数用于提高多目标搜寻的效率。
对于一个多目标最优化问题来说，边界控制非常重要。
如果个体离开了可行域，为了使这些个体仍然有效，有两种基于经验的策略可以执行：其一，产生新的个体代替移出可行域的个体；其二，把移出可行域的个体放置在边界上，并且改变它的前进方向。
考虑到要保持个体的多样性，并且需要在整个区域内搜寻最优解，在进入下一个优化周期时，以某一概率给出一个非可行解的边界，以代替之前的非可行解边界。
下面我们简要概括多目标细菌觅食优化算法的步骤： 各个参数所代表的意义如下p :搜寻空间的维度； S :细菌数量； c N : 趋药过程的步数； sN : 直行的步数； re N : 繁殖的次数； edN :迁徙过程的步数； ed P :消亡的概率；()C i :每次直行或翻转的行程长度。
1：初始化参数,,,,,,,()(1,2,,),ic s re ed edp s N N N N P C i i S θ=…。
2：迁徙循环：1l l =+。
3：繁殖循环：1k k =+。
4：趋药过程循环：1j j =+a ：当1,2,...,i l S ==时，按照下面的步骤对细菌i 执行趋药行为。
b ：计算两个适应值函数),,,(),,,,(21l k j i J l k j i J 。
c ：保留函数值),,,(11l k j i J J last =,),,,(22l k j i J J last =。
d ：翻转：产生一个随机向量()ni R ∆∈，其中的每个元素()m i ∆, 1,2,...,m S =, 是[]1,1-之间的一个随机数。
e ：移动：按照式（3）更新位置。