污水处理模型(最终版).doc

#*污水处理模型污水处理模型摘要摘要随着经济的快速发展，环保问题已经成为一个不容忽视的问题，而水资源更是关系着每个居民的日常生活，因此对于污水处理这一特殊的问题我们在解决时就应该本着高效的原则去实施，在这个污水处理问题中，我们先建立了一般情况下的模型，然后将该模型应用到实际问题中从而解决了实际问题。在模型的建立中我们要考虑工厂的净化能力，江水的自净能力，在保证江水经这一系列的处理后在到达下一个居民点后要达到国家标准，还要花费最少，对该问题进行全面的分析后可知这是一个运筹学方面关于线性规划的最优解问题，在该模型的建立中我们针对江水污水浓度在每个居民点之前小于国家标准这一条件对其建立线性约束条件，然后综合考虑费用最小，在结合三个处理厂各自的情况后，关于费用抽象数模型的目标函数，运用 LINGO9.0 规划软件求解，最后求得使江面上所有地段的水污染浓度达到国家标准时的最小费用为 5 万元。关键词： 污水处理 自净系数 污水流量 处理系数 污水浓度#*一、一、问题重述问题重述如下图，由若干工厂的污水经排污口流入某江，各口有污水处理站，处理 站对面是居民点。工厂 1 上游江水流量和污水浓度，国家标准规定的水的污染 浓度，以及各个工厂的污水流量和污水浓度都已知道。设污水处理费用与污水 处理前后的浓度差和污水流量成正比，使每单位流量的污水下降一个浓度单位 需要的处理费用（称处理系数）为已知，处理后的污水与江水混合，流到下一 个排污口之前，自然状态下江水也会使污水浓度降低一个比例系数（称自净系 数）该系数可以估计。试确定各污水处理站出口的污水浓度，使在符合国家标 准规定的条件下总的处理费用最小。先建立一般情况下的数学模型，再求解以下的具体问题：设上游江水流量为，污水浓度为,三个工厂的污水min/10100012llmg /8 . 0流量均为，污水浓度（从上游到下游排列）分别为 100，60，50（min/10512l） ，处理系数均为 1 万元/，3 个工厂之间的两段江面lmg /)/(min)/10(12lmgl的自净系数（从上游到下游）分别为 0.9，0.6。国家规定的污水浓度不能超过 1。lmg /（1）为了使江面上所有地段的水污染达到国家标准，最少需要花费多少 费用？ （2）如果只要求三个居民点上游的水污染达到国家标准，最少需要花费工厂 1工厂 2工厂 3处理站 1处理站 2处理站 3 江水居民点 1居民点 2居民点 3#*多少费用？二、二、问题分析问题分析通过对该污水处理所花费用最少问题的分析，我们可知在此问题中有多个 污水浓度，江水的原始污水浓度，工厂排出的污水浓度，处理厂排出的污水浓 度，以及当处理厂排出污水与江水混合后再经江水自净后的浓度，在这几个浓 度中只有经处理厂排出的污水的浓度是未知的，其关系着整个问题，要使总费 用最少，江中每段的污水浓度都达到国家标准，江水中污水浓度在到达下一居民点之前须达到国家标准 1（），那么问题的重点就在于对污水浓度的认lmg /识。在问题中有三个工厂以及对应的三个污水处理厂，那么这三个污水处理厂 各向江中投放的污水浓度就要有一个界值，又因当处理厂将污水排到江中之后 污水会随着江水不断向下游移动，因此下游污水的浓度与上游污水的浓度是紧 密相关的，即江面中每段污水的浓度都是有联系的，在模型的建立过程中我们 就要考虑应用递推的方法进行相邻两端之间污水浓度的联系，在问题的求解中 因所花费用都是用来对污水的处理，因此对个处理厂排出的污水浓度的确定就 显得至关重要，只有确定了这三个未知数即这三个界值后，我们才能建立目标 函数从而进一步得到最小花费。基于对江水浓度的限定与对花费最少两方面的 考虑，我们建立了线性规划模型。具体问题分析如下：对于第一个问题(1)为了使江面上所有地段的水污染达到国家标准,最少需要花费多少费用的解也就是说对于工厂 1 所排出的污水经过污水处理厂处理后的污水与江水混合 后的污水浓度就得达到国家标准。同时工厂 2，3 排出的经过处理的污水与江水 经过自净的水混合后也要达到国家标准。这样在求解具体问题的时候每个限制 条件在江水与工厂排出的水混合时进行设定。对于第二个问题(2)如果只要求三个居民点上游的水污染达到国家标准最少需要花费多少费用,对居民点 1 来说其上游的江水污水浓度为 0.8(),低于国家的标准污水lmg /浓度，无需考虑。也就是说在第二，三个居民点之前，污水浓度必须达到国家 标准，此时处理问题的限制条件发生在第二三个居民点处。这时工厂 1 排出的 污水经过污水厂的处理之后与江水混合，再经过江水自净到达居民点 2 之前须 达到国家标准，居民点 3 同理。#*三、三、模型假设模型假设(1) 河水的水流量和污水浓度短时间内不受天气与居民用水影响，只与工厂的 排放有关； (2) 河水的自我净化能力在短时间内不会发生改变；既自净系数不变； (3) 工厂排出的污水能在很短的时间内很好地与江水均匀融合； (4) 各污染物之间不会发生化学反应，也没有物理沉淀； (5) 工厂均能正常运作，不发生任何事故； (6) 河水和工厂的水流量均衡，污染物浓度平均；四、四、符号定义及模型假设符号定义及模型假设符号定义：符号定义：表示第 段江水的流量iQi表示各工厂排出污水的流量iS表示第 段江水中污水的浓度iCi表示第 个污水厂的污水浓度iAi表示第 个处理厂的污水浓度iXi表示江水与处理厂的污水混合后的污水浓度iD表示第 个处理厂的处理系数iRi表示第 段江面的自净系数iti表示所花费用M表示国家规定的污水浓度，其中 =1OCOClmg /#*模型假设模型假设：设有个工厂，个处理厂与个居民点，模型中部分相关参数在途中已进行表nnn 示如下所示：工厂， 污水浓度， 流量；1i1iA1iS工厂， 污水浓度， 流量；iiAiS处理厂 1， 污水浓度， 流量；1X1S处理厂， 污水浓度， 流量；iiXiS处理厂， 污水浓度， 流量；1i1iX1iS江水流量为，江水上游污水浓度为，各水段自净系数为； iQiCit工厂 1， 污水浓度， 流量；1A1S居民点 居民点， 居民点 。1i1i当处理厂将污水处理完排放到江中之后，居民点 1 即要取水，此时所要满足的条件是（为了解决问题方便不妨假设）021SSSSIOIICSQXSCQ)/(0111同理对居民点其所满足的为,其中i01CC ) 1(/()1) 1() 1(0000SSiQXSCSiQXtCCiiiiiiii假设花费为则有M目标函数： )(0iiiXASRM)1(ni#*五、五、模型的建立及求解模型的建立及求解模型的建立：模型的建立：对问题进行一般化处理后我们建立一般化的模型如下：目标函数：min )(0iiiXASRM)1(ni线性约束条件： )/()(00SQXSCQDiiiiiiiiDtC1ts0CDiiiAX 模型求解：模型求解： 在上面的一般模型中我们比较仔细的考虑了江水流量与处理厂的流量问题， 但在现实生活中因污水处理厂的处理能力有限，因此其流量相对于江水流量而言较小，我们对其进行理想化的处理即整个江水的流量为一常数，在求解 段Qi江面的混合污水浓度时忽略污水厂的流量。得到的简化模型如下所示：min )(0iiiXASRM)3 , 2 , 1( iiiiiQXSCD/0iiiDtC1ts0CDiiiAX 对于问题对于问题(1)(1)求解：求解：min )50(5)60(5)100(5321XXXM#*2 . 0005. 01X28. 0005. 00045. 021XX568. 0005. 0003. 00027. 0321XXXts1001X602X503X利用 lingo 求解可得当，时，.401X202X503X500M所以要想使江面所有地段均达到国家标准，所花最小费用为 500 万元。对于问题二求解：对于问题二求解：min )50(5)60(5)100(5321XXXM2 . 00045. 01X568. 0003. 00027. 021XXts1001X602X503X利用 lingo 求解可得当，时，,所222225.621X602X503X8889.188M以要使个居民点上游江水均达到国家标准，所花最少费用为 188.8889 万元。#*六、六、模型的评价模型的评价优点：优点：1) 该方案简单易行,原理清晰,依据可靠,论证有力,结论最优2 ) 该模型将现实中的污水处理问题用简单的线性规划问题进行分析计算，结 构简单，计算方便，有利于对相似问题进行求解和对模型进行扩充，比如工厂 的流水作业问题，物品运输问题，空气污染净化等问题的建模求解。3) 此问题所建立的模型是从一般问题到特殊问题的过渡，所用的数学方法为 线性规划，易于用多种数学软件编程求解，例如 LINDO，C+，MATLAB 等。缺点：缺点：1该模型在处理此问题时有假设与理想化的思想，与实际问题的求解还有一定 的距离，比如这三个污水厂排出的污水流量相等，实际中居民点是一个面，再 此模型中将其看作了一个点来进行处理2) 模型只从费用单方面考虑，忽略了处理厂与江水流量变化等的实际问题，使 得模型的建立偏离一定实际，从而计算结果不准确。七、七、参考文献参考文献1、谭永基，蔡志杰. 数学模型M.上海：复旦大学出版社. 2005 2、薛定全，陈阳泉. 高等应用数学问题的 MATLAB 求解M.北京：清华大学 出版社.2004 3、郑汉鼎，刁在筠编著 数学规划M.,济南：山东教育出版社，1997 4、谢金星，薛毅编著 优化建模与 LINDO/LINGO 软件M. 北京：清华大学 出版社 2005附录：附录：(1)Min 5A1-5X1+5A2-5X2+5A3-5X3s.t0.005X1<=0.20.0045X1+0.005X2<=0.28#*0.0027X1+0.003X2+0.005X3<=0.568X1<=100X2<=60X3<=50A1=100A2=60A3=50LP OPTIMUM FOUND AT STEP 2OBJECTIVE FUNCTION VALUE1) 500.0000VARIABLE VALUE REDUCED COSTA1 100.000000 0.000000X1 40.000000 0.000000A2 60.000000 0.000000X2 20.000002 0.000000A3 50.000000 0.000000X3 50.000000 0.000000ROW SLACK OR SURPLUS DUAL PRICES2) 0.000000 100.0000233) 0.000000 1000.0000004) 0.150000 0.0000005) 60.000000 0.0000006) 40.000000 0.0000007) 0.000000 5.0000008) 0.000000 -5.0000009) 0.000000 -5.00000010) 0.000000 -5.000000NO. ITERATIONS= 2RANGES IN WHICH THE BASIS IS UNCHANGED:#*OBJ COEFFICIENT RANGESVARIABLE CURRENT ALLOWABLE ALLOWABLECOEF INCREASE DECREASEA1 5.000000 INFINITY INFINITYX1 -5.000000 0.500000 INFINITYA2 5.000000 INFINITY INFINITYX2 -5.000000 5.000000 0.555556A3 5.000000 INFINITY INFINITYX3 -5.000000 5.000000 INFINITYRIGHTHAND SIDE RANGESROW CURRENT ALLOWABLE ALLOWABLERHS INCREASE DECREASE2 0.200000 0.111111 0.2000003 0.280000 0.200000 0.1000004 0.568000 INFINITY 0.1500005 100.000000 INFINITY 60.0000006 60.000000 INFINITY 40.0000007 50.000000 30.000006 50.0000008 100.000000 INFINITY 100.0000009 60.000000 INFINITY 60.00000010 50.000000 INFINITY 50.000000(2)Min 5A1-5X1+5A2-5X2+5A3-5X3st0.0045X1<=0.280.0027X1+0.003X2<=0.568X1<=100X2<=60X3<=50A1=100A2=60A3=50#*LP OPTIMUM FOUND AT STEP 1OBJECTIVE FUNCTION VALUE1) 188.8889VARIABLE VALUE REDUCED COSTA1 100.000000 0.000000X1 62.222225 0.000000A2 60.000000 0.000000X2 60.000000 0.000000A3 50.000000 0.000000X3 50.000000 0.000000ROW SLACK OR SURPLUS DUAL PRICES2) 0.000000 1111.1112063) 0.220000 0.0000004) 37.777775 0.0000005) 0.000000 5.0000006) 0.000000 5.0000007) 0.000000 -5.0000008) 0.000000 -5.0000009) 0.000000 -5.000000NO. ITERATIONS= 1RANGES IN WHICH THE BASIS IS UNCHANGED:OBJ COEFFICIENT RANGESVARIABLE CURRENT ALLOWABLE ALLOWABLECOEF INCREASE DECREASEA1 5.000000 INFINITY INFINITYX1 -5.000000 5.000000 INFINITYA2 5.000000 INFINITY INFINITYX2 -5.000000 5.000000 INFINITY#*A3 5.000000 INFINITY INFINITYX3 -5.000000 5.000000 INFINITYRIGHTHAND SIDE RANGESROW CURRENT ALLOWABLE ALLOWABLERHS INCREASE DECREASE2 0.280000 0.170000 0.2800003 0.568000 INFINITY 0.2200004 100.000000 INFINITY 37.7777755 60.000000 73.333336 60.0000006 50.000000 INFINITY 50.0000007 100.000000 INFINITY 100.0000008 60.000000 INFINITY 60.0000009 50.000000 INFINITY 50.000000