污水处理数学模型论文

污水处理数学模型论文成员：xwell随着经济的快速发展，环保问题已经成为一个不容忽视的问题，而水资源更是关系着每个居民的日常生活，因此对于污水处理这一特殊的问题我们在解决时就应该本着高效的原则去实施，在这个污水处理问题中，我们先建立了一般情况下的模型，然后将该模型应用到实际问题中从而解决了实际问题。在模型的建立中我们要考虑工厂的净化能力，江水的自净能力，在保证江水经这一系列的处理后在到达下一个居民点后要达到国家标准，还要花费最少，对该问题进行全面的分析后可知这是一个运筹学方面关于线性规划的最优解问题，在该模型的建立中我们针对江水污水浓度在每个居民点之前小于国家标准这一条件对其建立线性约束条件，然后综合考虑费用最小，在结合三个处理厂各自的情况后关于费用抽象数模型的目标函数，然后应用LINDO软件求解该问题得到当三个处理厂排出的污水浓度分别为40mg/l，20mg/l，50mg/l时，此时我们得到使江面上所有地段的水污染达到国家标准,最少需要花费费用为500万元。当从三个处理厂出来的污水浓度分别为62.222225mg/l，60mg/l，50mg/l,时，此时如果只要求三个居民点上游的水污染达到国家标准最少需要花费费用为188.8889万元。关键词：解决；污水；处理；分析问题的提出设上游江水流量为,污水浓度为0.8mg/l，3个工厂的污水流量均为,污水浓度(从上游到下游排列)分别为100,60,50(mg/l),处理系数均为1万元((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之前须达到国家标准，居民点3同理。模型假设假设与符号：假设在两个江面之间的江水流量和污水浓度在一小段范围内变化不大，即将一段江面看做一个点，我们对该点进行处理。Qi表示第i段江水的流量Si表示各工厂排出污水的流量.Ci表示第i段江水中污水的浓度Ai表示第i个污水厂的污水浓度。Xi表示第i个处理厂的污水浓度。Di表示江水与处理厂的污水混合后的污水浓度。Ri表示第i个处理厂的处理系数ti表示第i段江面的自净系数。M表示所花费用。C0表示国家规定的污水浓度其中C0=1mg/l对上述符号假设完，我们可以对上述问题进行模型的抽象，对此下建立一个一般情况下的模型。如下所示：设有n个工厂，n个处理厂与n个居民点，模型中部分相关参数在途中已进行表示如下所示工厂i+1污水浓度Ai+1流量Si+1工厂I污水浓度Ai流量Si处理厂1污水浓度X1流量S1处理厂I污水浓度Xi流量Si处理厂i+1污水浓度Xi+1流量Si+1江水流量为Qi，江水上游污水浓度为C1，各水段自净系数为ti工厂1污水浓度A1流量S1居民点1居民点i居民点i+11当处理厂江污水处理完排放到江中之后，居民点1即要取水，此时所要满足的条件是（为了解决问题方便不妨假设S1=S2=Si=S0）（Qi*C1+S1*X1）/(Qi+S0)=C0同理对居民点i其所满足的为Ci=C0,其中Ci=(Ci-1)*ti+Xi*((Qi+(i-1)*S0)*Ci-1+S0*Xi)/((Qi+(i-1)*S0)+S0)假设花费为M则有目标函数：M=∑Ri*S0*(Ai-Xi)(i=1……n)模型的建立对问题进行一般化处理后我们建立一般化的模型如下：目标函数：minM=∑Ri*S0*(Ai-Xi)(i=1……n)线性约束条件：Di=（Qi*Ci+S0*Xi）/(Qi+S0)Ci+1=ti*Dis.tDi=C0Xi=Ai模型的简化与求解在上面的一般模型中我们比较仔细的考虑了江水流量与处理厂的流量问题，但在现实生活中因污水处理厂的处理能力有限，因此其流量相对于江水流量而言较小，我们对其进行理想化的处理即整个江水的流量为一常数Q=，在求解i段江面的混合污水浓度时忽略污水厂的流量。得到的简化模型如下所示：minM=∑Ri*S0*(Ai-Xi)(i=1,2,3)Di=Ci+S0*Xi/QiCi+1=ti*Dis.tDi=C0Xi=Ai针对下面的问题：设上游江水流量为,污水浓度为0.8mg/l，3个工厂的污水流量均为,污水浓度(从上游到下游排列)分别为100,60,50(mg/l),处理系数均为1万元((mg/l)),3个工厂之间的两段江面的自净系数(从上游到下游)分别为0.9和0.6.国家标准规定水的污染浓度不能超过1mg/l.(1)为了使江面上所有地段的水污染达到国家标准,最少需要花费多少费用?(2)如果只要求三个居民点上游的水污染达到国家标准最少需要花费多少费用?建立模型：对于问题(1)求解minM=5*(100-X1)+5(60-X2)+5(50-X3)0.005X1=0.20.0045X1+0.005X2=0.280.0027X1+0.003X2+0.005X3=0.568S.tX1=100X2=60X3=50利用lindo求解可得当X1=40，X2=20，X3=50时，M=500.所以要想使江面所有地段均达到国家标准，所花最小费用为500万元。对于问题二求解：minM=5*(100-X1)+5(60-X2)+5(50-X3)0.0045X1=0.20.0027X1+0.003X2=0.568S.TX1=100X2=60X3=50利用lindo求解可得当X1=62.222225，X2=60，X3=50时，M=188.8889所以要使个居民点上游江水均达到国家标准，所花早少费用为188.8889万元。模型的优缺点优点：1)该方案简单易行,原理清晰,依据可靠,论证有力,结论最优2)该模型将现实中的污水处理问题用简单的线性规划问题进行分析计算，结构简单，计算方便，有利于对相似问题进行求解和对模型进行扩充，比如工厂的流水作业问题，物品运输问题，空气污染净化等问题的建模求解。3)此问题所建立的模型是从一般问题到特殊问题的过渡，所用的数学方法为线性规划，易于用多种数学软件编程求解，例如LINDO，C++，Matlab等。缺点：1．该模型在处理此问题时有假设与理想化的思想，与实际问题的求解还有一定的距离，比如这三个污水厂排出的污水流量相等，实际中居民点是一个面，再此模型中将其看做了一个点来进行处理2)模型只从费用单方面考虑，忽略了处理厂与江水流量变化等的实际问题，使得模型的建立偏离一定实际，从而计算结果不准确。改进此模型不甚全面，可以适当修改以适应实际的应用。参考文献胡运权《运筹学基础与应用》第四版高等教育出版社刘承平《数学建模方法》.北京.高等教育出版社.雷功炎《数学模型讲义》北京大学出版社谢金星薛毅《优化建模与LINDO/LINGO软件》清华大学出版社姜启源《数学模型》第二版高等教育出版社附件:(1)Min5A1-5X1+5A2-5X2+5A3-5X3st0.005X1=0.20.0045X1+0.005X2=0.280.0027X1+0.003X2+0.005X3=0.568X1=100X2=60X3=50A1=100A2=60A3=50LPOPTIMUMFOUNDATSTEP2OBJECTIVEFUNCTIONVALUE1)500.0000VARIABLEVALUEREDUCEDCOSTA1100.0000000.000000X140.0000000.000000A260.0000000.000000X220.0000020.000000A350.0000000.000000X350.0000000.000000ROWSLACKORSURPLUSDUALPRICES2)0.000000100.0000233)0.0000001000.0000004)0.1500000.0000005)60.0000000.0000006)40.0000000.0000007)0.0000005.0000008)0.000000-5.0000009)0.000000-5.00000010)0.000000-5.000000NO.ITERATIONS=2RANGESINWHICHTHEBASISISUNCHANGED:OBJCOEFFICIENTRANGESVARIABLECURRENTALLOWABLEALLOWABLECOEFINCREASEDECREASEA15.000000INFINITYINFINITYX1-5.0000000.500000INFINITYA25.000000INFINITYINFINITYX2-5.0000005.0000000.555556A35.000000INFINITYINFINITYX3-5.0000005.000000INFINITYRIGHTHANDSIDERANGEROWCURRENTALLOWABLEALLOWABLERHSINCREASEDECREASE20.2000000.1111110.20000030.2800000.2000000.10000040.568000INFINITY0.1500005100.000000INFINITY60.000000660.000000INFINITY40.000000750.00000030.00000650.0000008100.000000INFINITY100.000000960.000000INFINITY60.0000001050.000000INFINITY50.000000(2)Min5A1-5X1+5A2-5X2+5A3-5X3st0.0045X1=0.280.0027X1+0.003X2=0.568X1=100X2=60X3=50A1=100A2=60A3=50LPOPTIMUMFOUNDATSTEP1OBJECTIVEFUNCTIONVALUE1)188.8889VARIABLEVALUEREDUCE