污水处理模型最终稿

污水处理数学模型论文
摘要：随着经济的快速发展，环保问题已经成为一个不容忽视的问题，而水资源更是关系着每个居民的日常生活，因此对于污水处理这一特殊的问题我们在解决时就应该本着高效的原则去实施，在这个污水处理问题中，我们先建立了一般情况下的模型，然后将该模型应用到实际问题中从而解决了实际问题。

在模型的建立中我们要考虑工厂的净化能力，江水的自净能力，在保证江水经这一系列的处理后在到达下一个居民点后要达到国家标准，还要花费最少，对该问题进行全面的分析后可知这是一个运筹学方面关于线性规划的最优解问题，在该模型的建立中我们针对江水污水浓度在每个居民点之前小于国家标准这一条件对其建立线性约束条件，然后综合考虑费用最小，在结合三个处理厂各自的情况后关于费用抽象数模型的目标函数，然后应用LINGO软件求解
问题的提出
设上游江水流量为1000*10^12(l/min) ,污水浓度为0.8 mg/l，3个工厂的污水流量均为5*10^12(l/min) ,污水浓度(从上游到下游排列)分别为
100,60,50(mg/l),处理系数均为1万元( (10^12(l/min)*(mg/l)),3个工厂之间的两段江面的自净系数从上游到下游)分别为0.9和0.6.国家标准规定水的污染浓度不能超过1mg/l.
(1)为了使江面上所有地段的水污染达到国家标准,最少需要花费多少费用?
(2)如果只要求三个居民点上游的水污染达到国家标准最少需要花费多少费用?
问题的分析
通过对该污水处理所花费用最少问题的分析，我们可知在此问题中有多个污水浓度，江水的原始污水浓度，工厂排出的污水浓度，处理厂排出的污水浓度，以及当处理厂排出污水与江水混合后再经江水自净后的浓度，在这几个浓度中只有经处理厂排出的污水的浓度是未知的，其关系着整个问题，要使总费用最少，江中每段的污水浓度都达到国家标准，江水中污水浓度在到达下一居民点之前须达到国家标准1（mg/l），那么问题的重点就在于对污水浓度的认识。

在问题中有三个工厂以及对应的三个污水处理厂，那么这三个污水处理厂各向江中投放的污水浓度就要有一个界值，又因当处理厂将污水排到江中之后污水会随着江水不断向下游移动，因此下游污水的浓度与上游污水的浓度是紧密相关的，即江面中每段污水的浓度都是有联系的，在模型的建立过程中我们就要考虑应用递推的方法进行相邻两端之间污水浓度的联系，在问题的求解中因所花费用都是用来对污水的处理，因此对个处理厂排出的污水浓度的确定就显得至关重要，只有确定了这三个未知数即这三个界值后，我们才能建立目标函数从而进一步得到最小花费。

基于对江水浓度的限定与对花费最少两方面的考虑，我们建立了线性规划模型。

具体问题分析如下：
对于第一个问题
(1)为了使江面上所有地段的水污染达到国家标准,最少需要花费多少费用的
解也就是说对于工厂1所排出的污水经过污水处理厂处理后的污水与江水混合后的污水浓度就得达到国家标准。

同时工厂2，3排出的经过处理的污水与江水经过自净的水混合后也要达到国家标准。

这样在求解具体问题的时候每个限制条件在江水与工厂排出的水混合时进行设定。

对于第二个问题
(2)如果只要求三个居民点上游的水污染达到国家标准最少需要花费多少费用,对居民点1来说其上游的江水污水浓度为0.8(mg/l),低于国家的标准污水
浓度，无需考虑。

工厂1排出的污水经过污水厂的处理之后与江水混合再经过江水自净到达居民点2 之前须达到国家标准。

工厂2排出的污水经过污水厂的处理之后与江水混合再经过江水自净到达居民点3 之前须达到国家标准。

理论上，溶液浓度的计算为：
含污浓度=（处理厂排污浓度×工厂污水流量+江水流量×江中污水浓度）/ （江水流量+工厂污水流量）
实际上，因为江水流量远大于工厂污水流量，所以
含污浓度=（处理厂排污浓度×工厂污水流量+江水流量×江中污水浓度）/ 江水流量
要使得水的污染浓度符合国家标准规定的浓度，合理考虑江水的自净能
力，可知当工厂各污水处理站排出污水与江水混合的那一处浓度最高。

所以，
只要使得A，B，C三点处污水浓度符合国家标准即可。

为了分析问题方便，设
出每个工厂处理后的污水浓度分别为X，Y，Z，并忽略各工厂的排污流量大
小，从而列出符合国家标准条件下浓度的不等式。

根据已知条件，又可得到各
工厂处理费用与其流量，处理污水浓度之间的关系，进而通过等量代换，得到
各个工厂处理费用与各处理厂处理后的污水浓度的不等式关系。

利用所建立的
数学模型，通过可行域，求得最优解，即在符合国家标准的排污浓度的条件
下，总的处理费用最小。

模型假设：
1. 假设江水的水流量固定，不会因为加入污水或改变污水浓度而改变。

2. 假设污水之间无反应，不会因为污水反应而改变污水量或污水浓度。

3. 假设居民区不产生污水。

4. 假设江水的自净作用对所有污水都有效。

5. 假设污水在进入江水之后是分布均匀的。

6. 假设污水在进入江水之后不会流入上游。

7. 假设江水进行自净作用时，不改变江水本身流量。

8. 假设在对进行污水处理时，不改变污水流量，只改变污水浓度
模型求解：
设有N个工厂和N个污水处理站，
用Ci 表示第i 段江水未和处理厂的污水混合前的污水浓度
Ai 表示第i个污水厂的污水浓度。

Xi 表示第i个处理厂后的污水浓度。

Di 表示江水与处理厂处理好后的污水混合后的污水浓度。

Ri 表示第i个处理厂的处理系数
ti 表示第i段江面的自净系数。

S 表示各工厂排出污水的流量也即经处理厂处理后排出的流量
Q 表示江水的流量
M 表示所花费用。

C0 表示国家规定的污水浓度
其中 C0=1mg/l
第一种情况：江面上所有地段的水污染达到国家标准
（1）对于居民点1
当处理厂江污水处理完排放到江中之后，居民点1即要取水，此时所要满足的条件是D1=（Q×C1+S×X1）/(Q+S)<=C0
（2）同理对居民点i
其所满足的为 Di<=C0,其中
Di=（(Q+(i-1)×S)×Ci+S×Xi)/((Q+(i-1)×S)+S)
目标函数：
M=∑Ri×S×(Ai-Xi)(i=1……n)
在上面的一般模型中我们比较仔细的考虑了江水流量与处理厂的流量问题，但在现实生活中因污水处理厂的处理能力有限，因此其流量相对于江水流量而言较小，我们对其进行理想化的处理即整个江水的流量为一常数Q ，在求解i段江面的混合污水浓度时忽略污水厂的流量。

得到的简化模型如下所示：
目标函数 min M=∑Ri×S×(Ai-Xi) (i=1,2,3)
Di=Ci+S×Xi/Q
Ci+1=ti×Di
约束条件 Di<=C0
第二种情况：只要求三个居民点上游的水污染达到国家标准
（1）对于居民点1
当处理厂江污水未排放到江中之后，居民点1即可取水，此时所要满足的条件是 C1<=C0
（2）同理对居民点i
其所满足的为 Ci<=C0,其中
Di=（(Q+(i-1)×S)×Ci+S×Xi)/((Q+(i-1)×S)+S)
Ci=Ti×Di-1
Xi<=Ai
目标函数：
M=∑Ri×S×(Ai-Xi)(i=1……n)
在上面的一般模型中我们比较仔细的考虑了江水流量与处理厂的流量问题，但在现实生活中因污水处理厂的处理能力有限，因此其流量相对于江水流量而言较小，我们对其进行理想化的处理即整个江水的流量为一常数Q ，在求解i段江面的混合污水浓度时忽略污水厂的流量。

得到的简化模型如下所示：
目标函数 min M=∑Ri×S×(Ai-Xi) (i=1,2,3)
Di=Ci+S×Xi/Qi
Ci+1=ti×Di
约束条件 Di<=C0
Xi<=Ai
针对下面的问题：
设上游江水流量为1000 *10^12(l/min),污水浓度为0.8mg/l ，3个工厂的污水流量均为 ,污水浓度(从上游到下游排列)分别为100,60,50(mg/l),处理系数均为1万元( 10^12(l/min)(mg/l)),3个工厂之间的两段江面的自净系数(从上游到下游)分别为0.9和0.6.国家标准规定水的污染浓度不能超过1mg/l.
(1)为了使江面上所有地段的水污染达到国家标准,最少需要花费多少费用?
(2)如果只要求三个居民点上游的水污染达到国家标准最少需要花费多少费用?
对于问题(1)求解
model:
min=5*(100-x1)+5*(60-x2)+5*(50-x3);
c1=0.8;
d1=0.8+5*x1/1000;
c2=0.9*d1;
d2=c2+5*x2/1000;
c3=0.6*d2;
d3=c3+5*x3/1000;
d1<=1;
d2<=1;
d3<=1;
x1<=100;
x2<=60;
x3<=50;
End
求解得
Global optimal solution found.
Objective value: 500.0000
Total solver iterations: 0
Variable Value Reduced Cost
X1 40.00000
0.000000
X2 20.00000
0.000000
X3 50.00000
0.000000
C1 0.8000000
0.000000
D1 1.000000
0.000000
C2 0.9000000
0.000000
D2 1.000000
0.000000
C3 0.6000000
0.000000
D3 0.8500000
0.000000
Row Slack or Surplus Dual Price
1 500.0000 -
1.000000
2 0.000000
0.000000
3 0.000000 -1000.000
4 0.000000 -1000.000
5 0.000000 -1000.000
6 0.000000
0.000000
7 0.000000
0.000000
8 0.000000 100.0000
9 0.000000 1000.000
10 0.1500000 0.000000
11 60.00000 0.000000
12 40.00000 0.000000
13 0.000000 5.000000
利用lingo 求解可得当X1=40，X2=20，X3=50时，M=500.
所以要想使江面所有地段均达到国家标准，所花最小费用为500万元。

对于问题二求解：
model:
min=5*(100-x1)+5*(60-x2)+5*(50-x3);
c1=0.8;
d1=0.8+5*x1/1000;
c2=0.9*d1;
d2=c2+5*x2/1000;
c3=0.6*d2;
d3=c3+5*x3/1000;
c1<=1;
c2<=1;
c3<=1;
x1<=100;
x2<=60;
x3<=50;
End
求解得
Global optimal solution found.
Objective value: 188.8889
Total solver iterations: 0
Variable Value Reduced Cost
X1 62.22222
0.000000
X2 60.00000
0.000000
X3 50.00000
0.000000
C1 0.8000000
0.000000
D1 1.111111
0.000000
C2 1.000000
0.000000
D2 1.300000
0.000000
C3 0.7800000
0.000000
D3 1.030000
0.000000
Row Slack or Surplus Dual Price
1 188.8889 -
1.000000
2 0.000000
0.000000
3 0.000000 -1000.000
4 0.000000 -1111.111
5 0.000000
0.000000
6 0.000000
0.000000
7 0.000000
0.000000
8 0.2000000
0.000000
9 0.000000 1111.111
10 0.2200000
0.000000
11 37.77778
0.000000
12 0.000000
5.000000
13 0.000000
5.000000
利用lingo 求解可得当X1=62.222225，X2=60，X3=50时，M=188.8889
所以要使个居民点上游江水均达到国家标准，所花早少费用为188.8889万元。

模型的优缺点
优点：
1)该方案简单易行,原理清晰,依据可靠,论证有力,结论最优
2 )该模型将现实中的污水处理问题用简单的线性规划问题进行分析计算，结构简单，计算方便，有利于对相似问题进行求解和对模型进行扩充，比如工厂的流水作业问题，物品运输问题，空气污染净化等问题的建模求解。

缺点：
1．该模型在处理此问题时有假设与理想化的思想，与实际问题的求解还有一定的距离，比如这三个污水厂排出的污水流量相等，实际中居民点是一个面，再此模型中将其看做了一个点来进行处理
2) 模型只从费用单方面考虑，忽略了处理厂与江水流量变化等的实际问题，使得模型的建立偏离一定实际，从而计算结果不准确。

3.工厂只考虑汇入河流主干道的情况，但实际情况中会有许多支流汇入，从而又降低了污水浓度，所以我们可以考虑有多个支流汇入江中的情况。

4.当计算污水处理后汇入江中的浓度Di时，忽略了工厂的污水流量Q从而使计算过程得以简化，但在实际生活中，若工厂流量不能被忽略时，此计算结果将会使得所得混合浓度偏大，从而影响处理费用的最优求解。