污水处理问题的研究 16页

摘要随着社会的不断发展，排污问题备受关注，水资源密系民生，环保问题不容忽视。本文就工厂污水处理费用问题做了合理分析。针对问题一，为使江面上所有地段的水污染达到国家标准且花费最少，并结合污水排放浓度和江水的自净能力，得到此问题属于运筹学中在约朿条件下求最优解的问题。基于对江水浓度的限定与对花费最少两方面的考虑，建立线性规划模型，最终利用族go求解结果为当三个处理厂排岀的污水浓度分别为41mg/l,2l.lmg/l,50/ng/l时，得到最少费用为489.5万元。针对问题二，分别就居民点在处理站前、处理站后和处理站正对而进行分析。由江水稀释公式，可得处理费用与污水浓度Z间的关系，求得以上三种情况的最小费用分别为：居民点在处理站前吋，当处理厂1、处理厂2和处理厂3出口的污水浓度依次为63.33加g/1、60mg/1和50加g/1吋，花费为183.3333万元；居民点在对应处理站后不远处或正对面为489.5万元。针对问题三，求总费用最少可等价于求网路铺设后最短管长，其屮江中管道长度按四倍算其等效长度。先假设只有一个污水处理厂，在平而上找一点，使其与其它六点构成的最小生成树，并求其等效长度，找岀其屮等效长度最小的点即为所求。算得最小费用为273.7545万元。按上述方法计算两个污水处理厂时的情况，求得最小费用为311.0369万元。由于三个以上处理厂费用已经高于建立一个处理厂费用，于是不可取。丁•是污水处理厂个数为1,位置为居民点1与工厂1连线与江岸的交点，最小费用为273.7545万元。新城为该管道铺设区的任意位置。针对问题四，考虑到污水排放量的大小，因此将污水处理厂的位置建设在工厂区一岸，新城区的选取位置在工厂区一•附近，具污水排放到工厂区一。考虑到污水汇合的情况，首先建立只有一个污水交汇点的情况，再增加污水交汇点，逐步对模型一进行优化，最后再在河岸两边增加两点，优化管道过河的费用，最终得到污水交汇点的个数为5个的时候所用的建设费用最少，此时费用为：148.6025万元。关键词：运筹学线性规划最优解最小生成树自净能力\n一、问题重述1.1问题背景上游江水流量为IOOOxIO121/min,污水浓度为0.8加g/1,3个T厂的污水流量均为5xIO121/min,污水浓度（从上游到下游排列）分别为100,60,50（zng/l）,处理系数均为1万元/（（l012l/min）x/ng/l）,3个工厂之间的两段江而的自净系数（从上游到下游）分别为0.9,0.6o国家标准规定水的污染浓度不能超过lmt?/lc1・2具体问题（1）为了使江面上所有地段的水污染达到国家标准，最少需要花费多少费用？（2）如果只耍求三个居民点上游的水污染达到国家标准，最少需耍花费多少费用？（3）为了扩展城市发展空间，改善环境，计划在江水的上游建一个10万人口的新城，并对T厂废水及所有居民的生活污水用管线连接起来进行统一处理，已知建一个污水处理厂的成本为120万元，在江岸两边建管线的费用为3万/km,跨江管线费用为普通的4倍,居民点1,2,3的污水流量均为0.5xl0,21/min,污水浓度（从上游到下游排列）分别为25,6,4（加g/1）,新城的污水流量控制在2xl0121/min以内，且污水浓度不超过8吨/I,处理系数均为1.2万7G/（（10,2l/min）x/ng/l）ifi:根据以上条件，确定新城的位置及污水处理厂的个数和位置。使在符合国家标准规定的条件下总的费用最小。（4）在实际的施工过程中，仅用一种规格的污水采集管是不够的，如在采集管汇合后，要采用更大口径的管道。市面上有以下儿种污水采集管，并根据所给表中数据，重新确定新城及污水处理厂的位置，使总的费用最小，并给出管道铺设方案。二、问题分析通过对该污水处理所花费用最少问题的分析，可知在此问题中有多个污水浓度，江水的原始污水浓度，工厂排出的污水浓度，处理厂排出的污水浓度，以及当处理厂排出污水与江水混合后再经江水自净后的浓度，在这儿个浓度中只有经处理厂排出的污水的浓度是未知的，其关系着整个问题，耍使总费用最少，江中每段的污水浓度都达到国家标准，那么江水中污水浓度在到达下一居民点Z前须达到国家标准l〃g/lO2.1问题一的分析针对问题一，为了使江面上所有地段的水污染达到国家标准,最少需耍花费多少费用的解也就是说对丁工厂1,2,3所排出的污水经过污水处理厂处理后的污水与江水混合后的污水浓度达到国家标准。此问题属于在约束条件下求最优解的问题，我们只需要将满足经处理站处理后的污水达到国家规定的污水浓度，并且以工厂处理费用为冃标函数，求其最小值即可。\n2.2问题二的分析针对问题二，题冃没有明确指明居民点的位置，只是说明在处理站对面是居民点，所以居民区可以有多种排法，它们可以在对丿应处理站前，也可在对丿应处理站后，或正对面，因此这里只假设对应居民区与处理站前后相距不远。因此可分为三种情况考虑。一是当居民点在对应处理站前，二是当居民点在对应处理站后时，三是居民点在对应处理站的正对面的情况。2.3问题三的分析针对问题三，为使水污染达到国家标准且费用最少，对于污水处理厂的个数以及新城位置的确定，可通过设三个工厂和三个居民点为六个污水排放点，求解出连接六个污水排放点以及污水处理厂的最短距离，然后进行计算所需费用即可。2.4问题四的分析对于问题四，考虑进实际的施工条件，污水采集管需要使用不同的规格，使所花费用为最少的情况，其实就是要找到所有污水排放点距离污水处理厂最近的情况，再结合各个污水排放点的污水流量选择相应合适规格的污水采集管，从而计算所花费用。对丁•新城位置的确定，为使花费最少，新城应离处理厂越近越好。三、模型假设I•假设长江的水流量固定，不会因为加入污水或改变污水浓度而改变。2•假设污水Z间无反应，不会因为污水反应而改变污水量或污水浓度。3•假设居民区不产生污水。4.假设江水的自净作用对所有污水都有效。5•假设污水在进入长江Z后是分布均匀的。6.假设在对进行污水处理时，不改变污水流量，只改变污水浓度。7.假设管道的的最大流量大于等T15xlO,21/min。8.假设任意居民点与任意工厂Z间都可以铺设直线管道。9•假设任意河段可以修直线管道。10•假设任一点可以建设污水处理厂且其处理污水的能力大于18.5xIO121/min。II•假设居民区上游任一点可以建设新城。\n四、符号说明符号含义F总的污水处理费用4工厂i的处理系数Ci\T厂i处理前的污水浓度Ci2工厂2•处理后的污水浓度Qi工厂i的污水流量q工厂2•所排出的污水经过污水处理厂处理后的污水与江水混合后的污水浓度长江上游的江水流量coo长江上游的江水浓度工厂1,2Z间两段江面的自净系数f2工厂2,3Z间两段江面的自净系数工厂1污水处理费用工厂1,2污水处理费用Z和S,等效管长（江中管道长度按普通管道长度的4倍计算）S°污水处理厂到江边的距离%江岸两边建管线的单位费用w建污水处理厂的费用Pi污水排放点0污水处理厂s必fpj污水排放点i到污水排放点j的距离SpjfO污水排放点到污水处理厂的距离比.相应规格管道单价五、模型建立与求解5.1问题一模型的建立与求解工厂的花费受污水浓度与国家污水浓度规定浓度限制，因此，为了使江面上所有地段的水污染达到国家标准且使花费最少，可用线性规划模型来解决此问题，建立以工厂处理费用为冃标函数，如下：(1)3minF=Q\n其中，F为总的污水处理费用；q为工厂i的处理系数；5为工厂i处理前的污水浓度；％为工厂，处理后的污水浓度；Q为工厂i的污水流量。考虑到长江上游污水浓度达到国家标准以及江水的自净能力特性，此题的约朿条件为处理厂和江水自净功能的共同作用下使得丁厂江水的浓度达到国家规定的污水浓度，因此有以下约朿条件：G+0i=5i(a)+。|)+02「22?c)2最后，结合以上公式代入数值进行计算得，SM7.0369如，费用为311.0369万元,居民区污水处理厂位置为C(6.0000,1.0847),I厂区污水处理厂位置为D(6.0000,7.5119)o在平而屮搜索建立污水处理厂的两个点得到其分布情况，如图3所示:图3建立两个污水处理厂的分布格局图由图3可以得到两个污水处理厂的大致分布情况，其分别分布在居民区和工厂区，幷与工厂2、居民点2大致连成一线。5.4问题四模型的建立与求解该问题在问题三基础上提出，故知污水处理厂个数为一个。本模型为污水处理厂位置和管道铺设作分析。综合考虑实际的施工条件、各个污水排放点的污水流量大小、汇合点的污水流量大小，合理选择相应合适规格的采集管，并找到污水排放点距离污水处理厂最近的位置，从而计算所花费用。将污水处理厂与管道接头均称作交汇点：由于受到污水流量的限制，所选管道根据实际情况作变化。本模型从以下情况作分析：首先考虑所有点与污水处理厂直接用管道相连，其费用：\nF二工S；xW,+So%+W(9)Z=1其中，s,为等效管长；叱为相应规格管道单价；S。为污水处理厂到江边的距离；％为江岸两边建管线的单位费用；W为建污水处理厂的费用。由表屮管道规格屮，选择满足条件的管道。最终算出其最小费用为：359.756万元。污水处理厂(6.0921,6.3439)。然后在上述条件下加入一个交汇点记为D。该点在工厂区或者居民区，使得区中三点先交汇到该点，再一起汇入处理厂，于是分两种情况作讨论。其H标函数：36minF二工S,・xW,+工S’x叱+S°»xWd(10)•I/一4式中SgD为。到D的距离，为。到D所需管道费用。污水处理厂位置(5.9244,4.3394),居民区交汇点(5.9357,1.5676),求得其最小费用为:273.7028万元。在上述条件Z下，在厂区再取一交汇点记为E,使工厂区污水汇入该点，然后再流入污水处理厂。其冃标函数：6min尸二工SjX叱+So»xW°+SgEX%(11)式中为。到E的距离，We为。到E所需管道费用，污水处理厂位置(6.6274,4.3007)居民区上交汇点(6.2151,1.3877),工厂区上交汇点(6.20066.2722)求得其最小费用为:272.8694万元。在上述情况下，再加入一个交汇点F,即在区中求两点使得区中各点到两点Z中距离较短的一点作管道，然后两点之间作管道，最终流入污水处理厂，使其费用最低。其冃标函数：6血I工SfWq+工(12)/=|式中J表示D、E、F中任意一点，R表示D、E、F、O中任意一点，S表路程，W表费用。污水处理厂位置(6.5334,4.4357),居民点1,2的交汇点(5.4853,1.1779),居民点1,2,3的交汇点(6.6397,2.0000),工厂区1,2的交汇点(6.0658,6.7412)。最少费用为271.8127万元。在上述条件下，像工厂区中再加入一交汇点G,在工厂区屮求两点使得区屮各点到两点Z屮距离较短的一点作管道，然后两点之间作管道，最终流入污水处理厂，使其费用最低。其冃标函数：6minF二工(⑶/=1式中丿•表示D、E、F、G中任意一点，k表示D、E、F、G、O中任意一点，S表路程，W表费用。代入数值计算，污水处理厂位置为(7.0714,3.3000),居民区1与居民区2交汇点(6.0000,0.6),居民区3与居民区1,2交汇点(8.1036,2.0000),工厂区3与工厂区1,2交汇点(8.0446,4.7460),工厂1与工厂2交汇(5.7993,7.1789),由于污水处理厂与交汇点F儿乎垂直于河道，因此不需要优化，其最小费用为：268.6025万元。\n由以上分析得到，在平面中搜索污水处理厂点的位置，得到污水处理厂及管道铺设分布情况，如图4所示：图4污水处理厂及管道铺设分布图由图4可知，污水处理厂建立在丁•厂区域与长江岸线相接的位置，其中连接丁•厂、居民点、四个汇合点和污水处理厂的路线即为最优管道铺设方案，得到所用花费最少。管道铺设方案:居民点1到居民点2用B种型号的管道铺设，居民2到交汇点F用D型号管道铺设，居民点3到F用B型号管道铺设，F交汇点到污水处理厂D用型号管道铺设，工厂1到交汇点E用G管道铺设，T厂2到交汇点E用F型管道铺设，交汇点E到交汇点G用H管道铺设，工厂3到交汇点G用F管道铺设，交汇点G到污水处理厂用I管道铺设。由以上情况可知，交汇点把所有可能的管道铺设情况包含在内，于是可知，在最后一种情况得到最优解，最小花费为268.6025万元。六、模型评价与推广6.1模型评价:模型优点：（1）建立的线性规划模型简单易懂，原理清晰，依据可靠，论证有力，结论最优，\n对于不是很复杂的问题可以比较准确地计算。（2）该模型将现实中的污水处理问题用简单的线性规划问题进行分析计算，结构简单，计算方便，有利于对相似问题进行求解和对模型进行扩充，比如工厂的流水作业问题，物品运输问题，空气污染净化等问题的建模求解。（3）此问题所建立的模型是从一般问题到特殊问题的过渡，所用的数学方法为线性规划，易于用多种数学软件编程求解，例如ZJNGO,C++,matlab等。模型缺点：（1）该模型在处理此问题时有假设与理想化的思想，与实际问题的求解还有一定的距离，比如这三个污水厂排出的污水流量相等，实际中居民点是一个面，再此模型中将其看做了一个点来进行处理。（2）模型只从费用单方面考虑，忽略了处理厂与江水流量变化等的实际问题，使得模型的建立偏离一定实际，从而计算结果不太准确。6.2模型推广：通过对题冃的解读，不难发现这是一类线性规划最优解问题。我们建立了一个效益的线性规划模型，此类模型可以做到费用最低、最经济的情况，可以在现实中解决诸多效益问题。线性规划最优解是运筹学的一个重要分支。它在解决工业生产组织、经济计划、组织管理人机系统中，都发挥着重要的作用。本文模型的建立是为了解决污水处理的费用最少问题，决策者要通过概念抽彖、关系分析可将各类影响因子放入规划模型屮，可以通过相关的计算机软件得到兼顾全局的最优解。本题的求解是一个典型的线性规划最优解问题，该模型的使用范用非常广泛，规划模型在工业、商业、交通运输、工程技术、行政管理等领域有着广泛的应用。七、参考文献[1]王晓银，周保平，数学建模与数学实验，北京：科学出版社，2012.[2]刘竞，王静，通过数学建模解决污水处理问题，石家庄职业技术学院学报,第23卷第2期:74-75.[3]谢金星，薛毅，优化建模与UNDO/LING。软件，北京：清华大学出版社，2005.⑷齐微，曲线拟合的MATLAB实现和优化度检验,123000：1-6,2012。[5]司守荃孙玺簣,数学建模算法与应用,北京，国防工业出版社,2011o[6]肖华勇，数学建模竞赛优秀论文精选与点评，西安，西北工业大学出版社,2011o\n八、附件问题一的程序：min二zl+z2+z3;cl2*5-1005*cl+800二0;c1*904.5+5*c22-1010*c2=0;c2*606+5*c32-1015*c3二0;cl<=l;c2<=l;c3〈二1;zl=@if(cl2#le#100,5,0)*@abs(100-cl2);z2=@if(c22#lett60,5,0)*@abs(60-c22);z3=@if(c32#lett50,5,0)*@abs(50-c32);zl>=0;z2>=0;z3>=0;end8.2建立单个污水处理厂的程序:functionresults=dd(a,b)k=((a(l)-b(l))A2+(a(2)-b(2))A2)A(1/2);results=k;endCOline([0,15],[2,2],'Linewidth',2,'Color，，‘r')holdonline([0,15],[4.3,4.3],'Linewidth，,2,'Color'，‘r')holdonx=[l,6,12.8,1,6,12.8];y=[7.3,8.3,7.8,1,0.6,1.2];plot(x,y,'・’，‘Markersize',20,'color'，‘b');str二char('Pl','P2','P3'，‘P4‘，'P5','P6');a=[x;y];fori=l:6[aa,b]=ginput(1);holdontext(aa,b,str(i))%holdonendgridon%已知点z间的距离fori二1:6forj二1:6ifj>二i%y=ff(a(:,i),a(:,j));P=a(2,i);q=a(2,j);if(p>4.3&q>4・3)|(p<2&q<2)dddd(i,j)=dd(a(:,i),a(:,j));elseal=a(:,i);a2=a(:,j);\nk=(al(2)-a2(2))/(al(l)-a2(l));b=al(2)-k*al(1);xxlf=solve('k*x+b-2=0');xxl=eval(xxlf);xx2f=solve(，k*x+b-4.3=0,);xx2=eval(xx2f);yyl=2;yy2=4.3;xxxl=[xxl,yyl];xxx2=[xx2,yy2];ddd(i)=dd(xxxl,xxx2);dl(i)=dd(a(:,i),a(:,j));dddd(i,j)二3*ddd(i)+dl(i);ifabs(j-i)==3dddd(i,j)二dl(i)+3*2・3;endendenddddd(j,i)=dddd(i,j);endendddddpause;%污水站与各点Z间直接距离xO=l:0.5:12.8;y00=0.6:0.5:2;yll=4.3:0.5:8.3;y0=[y00,yl1];fork=l:length(xO)forr=l:length(yO)ify0(r)<=2gg=2-y0(r);elsegg=y0(r)-4.3;endaO=[xO(k);yO(r)];fori=l:6al=a(:,i);%y=ff(a(:,i),aO);ifaO(l)^=al(1)kk=(al(2)-a0(2))/(al(l)-a0(l));b=al(2)-kk*al(1);%y=k*x+bxxllf二solve('kk*x+b-2=0');xxll=eval(xxllf);xx22f二solve('kk*x+b-4・3二O');xx22=eval(xx22f);if(xxlKl)|(xxll>12.8)|(xx2212.8)dis(i)=dd(a(:,i),aO);%pauseelseyyll=2;yy22=4.3;xxxll=[xxll,yyl1];xxx22=[xx22,yy22];ddis=dd(xxxll,xxx22);disl(i)=dd(a(:,i),aO);dis(i)=3*ddis+disl(i);end\nelsee=al(2);ee=aO(2);if(e>4.3&ee>4.3)|(e<2&ee<2)dis(i)=dd(a(:,i),aO);elsedis(i)=dd(a(:,i),aO)+3*2.3;endendend%dis%判断哪段路程最短fori=l:6forj=l:6if1=Jaaa(i)=100;ifaaa(i)>dddd(i,j)aaa(i)=dddd(i,j);elseaaa(i)=aaa(i);endifaaa(i)==dddd(i,j)dddd(j,i)=100;elsedddd(j,i)=dddd(j,i);endelseifi>laaa(i)zzaaa(i-l);elseaaa(i)=100;endendifa.aa(i)>dis(i)aaa(i)=dis(i);elseaaa(i)=aaa(i);endendendaaayy(k,r)=sum(aaa)+gg;endend[c,b]=min(yy(:))[rows,cols]=find(yy==c)\nholdon;plot((rows-1)*0.5+1,(cols-1)*0.5+0.6,'・','Markersize',20,'color5,'r)；holdonline([l,6],[7.3,8.3],'Linewidth',2,'Color，，‘k')line([6,12.8],[8.3,7.8],'Linewidth',2,'Color','k')line([l,1],[2.1,7.3],'Linewidth',2,'Color','k')line([1,1],[1,2.1],'Linewidth',2,'Color'，‘k')line([l,6],[1,0.6],'Linewidth',2,'Color'，‘k')line([6,12.8],[0.6,1.2]/Linewidth',2,'Color',，k')8.3建立两个污水处理厂的程序：funlm文件:functionf=funl(x);f二(x(l厂2+x(2厂2)'(l/2)+((x(l)-5)'2+(x(2)+0・4厂2厂(l/2)+((x(l)-ll.8厂2+(x(2)-0.2厂2厂(1/2)+(l-x(2))Fun2m文件：functionf二fun2(x);f=(x(l厂2+(x(2)-6・3厂2厂(1/2)+((x⑴-5厂2+(x⑵-7・3厂2厂(l/2)+((x⑴-11.8厂2+(x(2)-6.8厂2厂(1/2)+(x(2)-3.3);执行程序：fminunc('fun2',[50.4])fminunc(，fun2,,[56.3])y=[0,5,11.8]；z=[0,-0.4,0.2];plot(y,z,';text(0,0,'居民点一')text(5,-0.4,'居民点二')text(ll.8,0.2/居民点三')holdonx二⑸；y=[0.0847];text(5,0.0845,'污水处理厂')plot(x,y,'+')holdonx=-2:1:15w二ones(1,18);plot(x,w)holdony二[0,5,11.8];z=[6.3,7.3,6.8];plot(y,z,';text(0,6.3,'工厂一')text(5,7.3,'工厂二')text(ll・8,6・8,'工厂三')holdonx二⑸；y=[6.5119];text(5,6.5119,'污水处理厂’)plot(x,y,'+')holdonx=-2:l:15w=ones(1,18);plot(x,3.3*w)axis([-215-39])