污水排放模型
第三次作业
张峰华 材料学院 材料成型及控制工程04班 20123631 刘泽 材料学院 材料成型及控制工程04班 20123627 杨海鹏 材料学院 冶金工程03班 20123203
一、问题重述
如下图，有若干工厂的污水经排污口流入某江，各口有污水处理站，处理站对面是居民点。工厂1上游江水流量和污水浓度，国家标准规定的水的污染程度，以及各个工厂的污水流量和污水浓度均已知道。设污水处理费用与污水处理前后的浓度差和污水流量成正比，使每单位流量的污水下降一个浓度单位需要的处理费用（称处理系数）为已知。处理后的污水与江水混合，流到下一个排污口之前，自然状态下的江水也会使污水浓度下降一个比例系数（称自净系数），该系数可以估计。试确定各污水处理站出口的污水浓度，使在符合国家标准规定的条件下总的处理费用最少。
工厂1 2 工厂 工厂3 处理站1 处理站 2 处理站3
江水 居民点1 居民点2 居民点3 先建立一般情况下的数学模型，再求解以下的具体问题： 设上游江水流量为1000的污水流量均为5?1012L/min，污水浓度为
0.8mg/L，3
个工厂
?1012L/min，污水浓度（从上游到下游排列）分别为
，
100,60,50(mg/L)12((L/min)?(mg/L))，处理系数均为1万元/103个工厂之间的两段江面的自净系数（从上游到下游）分别为0.9，0.6。国家标准规定水的污染浓度不能超过
1mg/L。
（1）为了使江水上所有地段的水污染达到国家标准，最少需要花费多少费
用？
（2）如果只要求三个居民点上游的水污染达到国家标准，最少需要花费多
少费用？
二、问题分析
1
问题（1）为了使江面上所有地段的水污染达到国家标准，求最少需要的花费。此问题属于在约束条件下求最优解的问题，我们只需要将满足经处理站处理后的污水达到国家规定的污水浓度，并且以工厂处理费用为目标函数，求其最小值即可。
问题（2）若只要求三个居民点上游的水污染达到国家标准，求最少需要的花费。考虑到江水有自净的功能，这样处理厂就不必将污水处理成符合国家规定的污水浓度，只要使每个居民点上游的江水再其自净作用过后能达到国家规定的污水浓度就可以。那么此题的约束条件就变成了在处理厂和江水自净功能的共同作用下使得工厂江水的浓度达到国家规定的污水浓度。目标函数依然是工厂处理费用。
又因为题目中没有明确的说明居民点的具体位置，只是说在处理站的对面，
那么根据居民点位置的不同再利用上述思路进行分析和求解。
三、模型的假设和符号说明
1.假设长江的水流速固定，不会因为加入污水或改变污水浓度而改变。 2.假设污水之间无反应，不会因为污水反应而改变污水量或污水浓度。 3.假设居民区不产生污水。
4.假设江水的自净作用对所有污水都有效。 5.假设污水在进入长江之后是分布均匀的。 6.假设污水在进入长江之后不会流入上游。
7.假设江水进行自净作用时，不改变江水本身流量。
8.假设在对进行污水处理时，不改变污水流量，只改变污水浓度。 9.假设3个工厂之间的两段江面，各自单位距离的自净能力相同。
c0：上游江水的污染浓度
c11：工厂1产生的污水浓度 c12：工厂1经处理后的污水浓度 c21：工厂2产生的污水浓度 c22：工厂2经处理后的污水浓度 c31：工厂3产生的污水浓度 c32：工厂3经处理后的污水浓度 c1：从处理站1流出的污水与江水混合之后的污水浓度 c2：从处理站2流出的污水与江水混合之后的污水浓度 c3：从处理站3流出的污水与江水混合之后的污水浓度
vvv
0：江水上游水流量
v：工厂1产生的污水流量 b1：工厂1的处理系数
12：
工厂2产生的污水流量
b2：工厂2的处理系数 b3：工厂3的处理系数
3：工厂3产生的污水流量
2
aazz12：
工厂1、2之间江面的自净系数 工厂2、3之间江面的自净系数
23：
z：工厂1的污水处理费
12：工厂2的污水处理费
3：
工厂3的污水处理费
Z：3个工厂总费用 i=1,2,3分别代表工厂1，工厂2，工厂3
五、模型的建立与求解
题目中已知的条件有：c0=0.8mg/L
c11=100mg/L c21=60mg/L c31=50mg/L
12v0=1000?10L/min
12v=v=v123=5?10L/min
a12=0.9
a23=0.6
江水本有的污水量=江水流量（v1）×江水污水浓度（c0）
工厂i处理后的污水量=工厂i的污水量（vi）×处理后污水浓度（ci2）
ci=（江水本有的污水量+工厂i处理后的污水量）/（江水流量+工厂i的污水
量）
污水下降浓度=原污水浓度（ci1）-处理后的污水浓度（ci2）
工厂i处理污水的费用（zi）=污水下降浓度×工厂i的污水量（vi）×工厂i的污水处理系数（bi）
按照上述条件即可列出以下： 问题（1）的约束条件为：
?vc?c1?vc1v?v0011201
222c1?1000?0.8?5c12?11000?5 (1000?5)?0.9c1?5c22?11000?5?5 (100?0?5?5)c20.?6c325?11000??5?55
?v??vc?acvc2?v?v?v12101012232012?1332
c2?c3??v?v??vc?acvc3?v?v?v?v0123?1
z?0 （i=1,2,3）
i 3
目标函数为：
33?bminZ??zi??i?ci1ci2?i?1i?1 3minZ??zi?5?(100?c12)?5?(60?c22)?5?(50?c32)i?1
用lingo9.0求解后得：
Local optimal solution found.
Objective value: 489.5000 Extended solver steps: 0 Total solver iterations: 13
Variable Value Reduced Cost
Z1 295.0000 0.000000
Z2 194.5000 0.000000
Z3 0.2555717E-07 0.000000
C12 41.00000 0.000000
C1 1.000000 0.000000
C22 21.10000 0.000000
C2 1.000000 0.000000
C32 50.00000 0.000000
C3 0.8433498 0.000000
污水处理厂1、处理厂2和处理厂3出口的浓度依次为41.00 mg/l、21.10 mg/l和50.00 mg/l时，工厂1花费295万元，工厂2花费194.5万元，工厂3不用处理污水。所以为了使江水上所有地段的水污染达到国家标准，最少需要花费489.5万元。
对于问题（2）题目没有明确指明居民点的位置，只是说明在处理站对面是居民点，所以居民区可以有多种排法，它们可以在对应处理站前，也可在对应处理站后，或正对面，因此我们这里只假设对应居民区与处理站前后相距不远。 （1）当居民点在对应处理站前时如下图所示：
4