污水处理厂选址问题的数学模型.docVIP

污水处理厂选址问题的数学模型 1 引言 六七十年代以来，工业化国家为保护生态环境，已把发展污水处理系统作为重要问题提上议事日程．一种发展方向是在地区范围内统一考虑，而不是每个村镇各自为政．这样做是有充分理由的．一方面，在区域范围内建立一个污水处理系统，可以更好的考虑利用地形和河流；另一方面，建造与运行大型污水处理厂，所需的费用相对比较低，即工厂规模越大，处理单位数量污水的费用越低．以字母C表示处理一定量污水所需的费用，显然它是被净化的污水量x的函数，按照已发表的资料，在适当单位下，C正比于xa，0a1；一个为十万人口服务的处理厂所需费用仅是六个为一万人服务的小厂费用之和．污水处理厂所需费用与规模间的这种关系对本文所进行的讨论有重要意义． 2 污水处理厂选址模型 2.1 问题的表述与模型 一个地区的污水处理系统可由一个网络来描述，涉及所有的村镇，所有可能选择的处理厂位置，以及所有可能铺设的管网弧，弧的方向取决于污水的流向．以hi记结点i每年产生的污水量，显然hi0．如果结点i是可能建厂的结点，以非负的量xi表示建立后的处理厂每年净化的污水量．显然xi是未知的，它们必须按照某种优化原则解出；如果解得的xi取零值，则表示不在该处建厂．从结点i到结点j的弧记为（i，j），以非负的量yij表示每年通过该段弧的污水流量．如果在最后的解中yij=0，则说明这段管道无须修建．与结点i相关的一段管道以及一座污水处理厂的年运行与维修费用依次是yij与xi函数，分别记为bij(yij)和ai(xi)．这些函数的形式是已知的．此处不用考虑建设费用，因为在这里面我们可以把这笔费用分摊到每年的维护费用中去计算．把以上描述的网络记为G＝（N,A），其中N为结点集合，A为有向弧集合，那么在网络G上寻求最佳污水处理系统建设方案的问题表示为如下数学模型，即使目标函数 （1） 在约束条件 　　　　　　　　　 　　　　　　　　　　 （2） 下达到最小．这组约束表示在每一个结点处水量守恒，由于每一个结点都有一个相应的约束条件，所以约束总数恰恰与结点数相等，这与对目标函数的最优分析有重要的意义。 2.2 最优方案的定性分析 由上文可知，一个污水处理厂所需的费用是其基本处理能力的凹（或向上凸）函数．由此，目标函数C是其所有变量的凹函数．另一方面，所有满足约束条件的点构成一个凸集合即可行解集是凸集．由此利用凹函数的jensen不等式与凸集性质，不难知道目标函数C的极小值在可行解集的一个顶点处达到．这样，上述问题的求解可以类似于线性规划，由逐一检查凸多面体的顶点来完成． 而如上凸多面体的任何一个顶点，对应于约束方程组的一个可行解．因而取非零值变元xi与yij的数目最多等于约束方程个数，也就是说和网络结点数相等．这是因为每个结点处只有一个约束条件，即水量守恒的条件．这一结果表明：对任何结点i，xi与yij之中，只可能有一个量不为零．理由是：若解在某一个结点上有两个或两个以上的非零变量，为使所有非零变量数等于结点数，由鸽笼定理可知，则至少存在一个结点i，使xi，yij都是零，这意味着既不能在该处建立污水处理厂，又不允许污水流出，这样的结点显然是不可能或不必考虑的．对任何一个结点i，xi，yij中只有一个非零，说明或者该点是污水处理厂所在地，且工厂处理流到该点的全部污水；或者该点不设厂，且全部污水通过一条管道流出． 在一个有向图中，从一个结点出发的有向弧的数目称为该结点的出度．有上述可知，最佳污水处理方案所对应的网络，只包含出度为0或1的结点，而且显然不包含环路，除了一个结点外，所有结点出度为1，而唯一的例外结点出度为0，那么这一网络只能是树，例外结点是树的根．前面的分析表明，描述最佳污水处理系统的网络一般是由一组彼此不相连接的树组成，每个树的根结点数目就是污水处理厂的数目． 问题求解 3.1 模型求解 上述的数学模型有各种不同的求解方法，以下介绍其中的一个。由于污水处理厂选址的模型自然成树型结构．另外因为地形、河流、城镇等等的原因，使得水流向的固定，因此模型里面连接结点的弧都是有向弧．所以上述模型构成的就是有向图，也是树型结构． 因此，我们首先来考虑原来网络所有可能的生成树，即包括了网络所有的个结点但只利用了条弧的树，对每一生成树进行求解，然后再选择最优解． 为此，我们令图T=（N，A）是树，以S表示T的某个部分子树的结点集，同一子树弧集记为G（S），显然SN，G（S）A．所有的S构成的集合记为，对任意的S，(S)表示该子树的根，H（S）表示N的一个子集，他的任一元素不属于S，但有A中的一条弧，将其连接到S的某一结点；对任意N，K（）表示所有N中有路径连接到的结点的集合． 利用以上记号，重