污水排放与处理.

1、污水排放与处理摘要为了建立了污水处理站处理方案的优化模型，使污水源公司处理污水的投资最少，将问题归结为优化问题，运用优化理论及MATLAB软件解决各问题，在前三个问题的模型中认为自净系数恒定不变，随之将问题作为线性规划问题来处理，计算结果为：问题一的总资金最少为590万元；问题二的总资金最少为79.6万元；问题三中可以得到，总资金花费与污水源的流量，污水源的浓度，污水处理费用系数同增同减。并且当污水源的流量减少一定量时，排放污水量减少，会导致不用经过净化处理也能直接排放，并且不会影响居民生活。针对第四问，我们根据已知条件画出了和浓度C的散点图，利用MATLAB对散点图进行拟合，经检验最合适的为

2、二次多项式。利用此表达式对前三问进行新的研究。关键词：污水排放；控制；决策优化1、 问题的重述处理站1居民点1居民点2居民点3处理站3处理站4处理站201234污水源1污水源3污水源4污水源2 如图所示，某江沿岸有若干个污水源及为之配备的污水处理站和若干个居民点。污水源以的流量向污水处理站排放浓度为的污水，经污水处理站的净化处理后，污水以浓度（）通过排污口向江中排放。一般地，污水处理站净化污水的费用与处理前后的污水浓度差成正比，与所处理的污水量成正比，即，其中为比例系数，它反映了污水处理站每流量单位降低每个浓度单位所需的资金，也反映了污水处理的技术水平。设国家规定的江水水质标准为=1（mg/l

3、），与污水源、污水处理站相关的数据如下表所示。下标1234流量（l/min）5434浓度（mg/l）100705080（万元·mg·min）21.511.5以污水处理站的排水口为分界点，可以将该段江面划分为5个子段，各子段的水流量及水质分别为和。考虑到江中流水的自净作用，对每个子段引入自净系数。显然，自净系数与河流状态有着紧密地联系。在该段江水中，设=0.4，=0.5，=0.8，=0.7。设该段江水的来流的流量及污水浓度分别为=1000（l/min）和=0.8（mg/l）。1、在该段江水全面达到国家水质标准的前提下，如何使投入污水处理的总资金最少。2、在每个居民点上游水质达

4、到国家标准的条件下，如何使投入污水处理的总资金最少。3、当污水源的流量和浓度，污水处理费用系数的变化后，在上述两问中污水处理的总资金将发生什么变化？4、事实上，自净系数不仅与河流状态有关，还与江水的污水浓度有关。当江水的污水浓度较低时，流水发挥自净作用；当污水浓度逐渐增加时，流水的自净能力逐渐降低；当江水的污水浓度超过某个阈值时，流水就失去了自净功能。考虑适当的函数形式，重新研究上述三个问题。二、问题的分析前两问的目标函数相同，均为投入总资金，但前两问的分析重点在于所设环境不同，第一问求该段江水全面达到国家水质标准是的总投资最少，即是每个污水处理站的排放口处的浓度达到标准，第二问中要求每个居民

5、点上游的水质达到标准，则要求临近居民点附近的水质达标，注意居民点的位置，有的要考虑江水的自净作用，有的则不要考虑。写出前两问的限制条件，分别求解目标函数。第三问在前两问的基础上采用控制变量法，分别研究污水源的流量和浓度，污水处理费用系数的变化后总资金的变化。第四问，考虑自净系数受江水污水浓度影响根据已知条件画出了和浓度C的散点图，利用MATLAB对散点图进行拟合，可用的函数有反比例函数、一次函数、二次函数或三次函数这样的函数关系式，在对这些函数关系式进行检验找出最优函数关系式。三、基本假设1.从污水源出来的污水经过污水处理站后只有浓度的改变，不考虑体积的变化。2.本题考虑的自净系数是指江水自身

6、净化掉的杂质系数，之后的江水杂质浓度为1-。3.问题四中江面的自净系数由刚流入该江面的江水浓度决定，并且在江水净化过程中，自净系数系数保持不变。4.不经过时间变化，污水一旦排出，立刻与江中水混合均匀。4、 符号说明:第i个污水处理站净化污水的费用:第i个污水源向污水处理站排放污水的流量（单位：l/min）；:第i个污水源向污水处理站排放污水的浓度（单位：mg/L）；:经过第i个污水处理站的净化处理后通过排污口向江中排放的污水浓度（单位：mg/L）；:污水处理站净化污水的费用与处理前后的污水浓度差以及所处理的污水量的正比例系数（万元·mg·min）:江面第个子段的水流量（l/

7、min）:江面第个子段经过该子段流水净化后的水质（单位：mg/L）；:江面第i个子段未经过该子段流水净化时的水质，即排污口的污水水质（单位：mg/L）；:国家规定的江水水质标准为1（单位：mg/L）；:第j个子段的自净系数；5、 模型的建立与求解1. 问题一（1）决策变量：（2）决策目标：最小化污水处理的总资金，即：（4）约束条件：（4） 模型求解： 1）由于江水具有净化作用，所以每段江水只需考虑污水处理站排污口污水浓度的约束条件即可。 2）将上述模型输入MATLAB求解，可以得到：但由于、，故实际上结果为：即第2个，第3个，第4个污水处理点不用处理污水，也能达到国家水质标准，并使得总资金花费

8、最少，将值代入公式得到最小总资金为590万元。2. 问题二（1）决策变量：（2）决策目标：最小化污水处理的总资金，即：（3）约束条件：（4） 模型求解：1） 居民点1紧靠在污水处理站2的下游，居民点2和居民点3分别在污水处理站3和4的下游（非紧靠）。2） 要求临近居民点上游的水质达标，即污水处理站1不需考虑。3） 将上述模型输入MATLAB求解，可以得到：但由于、，故实际上结果为：即第1个，第3个，第4个污水处理点不用处理污水，也能达到国家水质标准，并使得总资金花费最少，为79.6万元。3. 问题三（1） 决策变量：、（2）决策目标：污水处理的总资金，即：（3）约束条件：同问题一及问题二（4）

9、模型求解：将上述模型输入MATLAB求解，可以得到：题号qurT处理费变化量第一问100790+33.9%第二问100217.2+183.9%第一问-100390-33.9%第二问-1000/第一问0100690+16.9%第二问0100136.5+78.4%第一问0-100490-16.9%第二问0-1000/第一问000.5737.5+25.0%第二问000.5102+33.3%第一问00-0.5442.5-25.0%第二问00-0.551-33.3%由上面的表格可知：（1）当固定两个变化量，改变其中一个时：q、u、r增加或减少时，处理费T增加或减少；（2）当三个量减小到一定程度时，处理费

10、会达到0。这个可以理解为，当污水的浓度、污水的排放量小到一定程度时，根本就无需建造污水处理站，这显然是我们所希望发生的。4.问题四 （1）函数拟合本题的问题中值在变化，可以用MATLAB先画出（C,）的散点图如下图所示： 通过一次、二次和三次拟合可以得到如图所示的拟合图线，并得到对应的如下函数表达式：一次拟合：；二次拟合：；三次拟合：，然后结合散点图比较上面三条曲线，一次函数显然没有二次函数和三次函数拟合的程度好，所以我们选择用二次或三次拟合，又由于二次和三次拟合结果相当，为了实际当中计算简便，我们选择采用二次函数来拟合自净系数与污水浓度的关系。（2）利用二次函数对上述三个问题重新进行研究。利

11、用迭代法得到问题一的结果如下：=0.800；= 0.959；= 0.934； =0.959=0.2402；=0.0764；=0.0962；=0.0766从结果可以看出浓度越小自净系数越大，符合题目要求。6、 模型的评价1、利用最优模型，充分考虑约束条件找到决策目标的最优解，具有实际意义。2、考虑自净系数受污水浓度的影响，更贴近实际，误差减小。3、污水处理厂排放污水后，混合需要一定时间，模型没有考虑这种不均匀性，会有误差。4、模型中将江水看作理想流体为了与实际结合更加紧密，模型有待进一步改进7、 模型的推广 据题目中提供的数据，建立了污水处理站处理方案的优化模型，并利用优化模型的相关知识，确定了

12、最小总投资，对污水处理方面有一定的参考价值。模型中分析问题、解决问题的一些独到方法，优化模型的总体思想，对其他数学问题及一般模型仍可使用。另外，针对江水污染等自然灾害，我们的方法对于环保部门可以作为分析解决问题的一种参考。8、 引用1洪 锋.沿江污水排放控制.湖南 长沙：国防科技大学9、 附录说明：代码中为；为问题一syms v1;syms v2;syms v3;syms v4;f1=a1/q1*v1+q0*c0/q1-c;f2=a2/q2*v2+a1*(1-b1)/q2*v1+q0*c0*(1-b1)/q2-c;f3=a3/q3*v3+a2*(1-b2)/q3*v2+a1*(1-b1)*(1

13、-b2)/q3*v1+q0*c0*(1-b1)*(1-b2)/q3-c;f4=a4/q4*v4+a3*(1-b3)/q4*v3+a2*(1-b2)*(1-b3)/q4*v2+a1*(1-b1)*(1-b2)*(1-b3)/q4*v1+q0*c0*(1-b1)*(1-b2)*(1-b3)/q4-c;v1,v2,v3,v4=solve(f1,f2,f3,f4,'v1','v2','v3','v4')v1 =41v2 =203/2v3 =1015/6v4 =1017/5问题二v1=100;f2=a2/q2*v2+a1*(1-b1)/q2

14、*100+q0*c0*(1-b1)/q2-c;f3=a3/q3*v3+a2*(1-b2)*(1-b3)/q3*v2+a1*(1-b1)*(1-b2)*(1-b3)/q3*v1+q0*c0*(1-b1)*(1-b2)*(1-b3)/q3-c;f4=a4/q4*v4+a3*(1-b3)*(1-b4)/q4*v3+a2*(1-b2)*(1-b3)*(1-b4)/q4*v2+a1*(1-b1)*(1-b2)*(1-b3)*(1-b4)/q4*v1+q0*c0*(1-b1)*(1-b2)*(1-b3)*(1-b4)/q4-c;v2,v3,v4=solve(f2,f3,f4,'v2','v3','v4')v2 =229/4v3 =3037/10v4 =116383/500问题四syms v1;syms v2;syms v3;syms v4;v1=41c1=a1/q1*v1+q0*c0/q1b2=1.8028*f12-4.2006*f1+2.4469v2=70c2=a2/q2*v2+a1*(1-b1)/q2*v1+q0*c0*(1-b1)/q2b3=1.8028*f22-4.200