数学建模长江水质的评价和预测.doc

摘要本文在给定数据的基础上，建立了水质综合评价模型；污染源依靠流量、流速和降解系数的模型；灰色预测模型，对未来十年污水治理做了预测。针对问题一，做出标准化的参数与相应权值，建立合理的综合评价函数，得出了各地各时间内的综合评价值，得到湖北丹江口水质最好、江西南昌谁知最差的结论。针对问题二，根据流量、流速和降解系数建立了各地段排污量的模型，得到高锰酸盐与氨氮排污量最大的地段都是湖北宜昌到湖南岳阳段。针对问题三、四，建立了灰色预测模型，并给出了污水处理方案。针对问题五，提出了整治长江污染的几点建议：加强宣传力度、加强有关部门监督、整治沿江工业。模型较全面的运用了所给数据，建模方法比较科学，但还存在具体数值设立上主观性的问题。关键词：综合评价、灰色预测1.问题重述1.1问题背景长江是我国第一、世界第三大河流，是我国唯一具有全国意义的战略水源地，是我国水资源供需平衡的最后防线。但是近几年的统计数据表明，长江水质污染日益严重，正面临着前所未有的六大危机：森林覆盖率严重下降，泥沙含量增加，生态环境急剧恶化；枯水期不断提前，长江断流日益逼近；水质严重恶化，重金属含量非常高，危及沿江许多城市的饮用水，癌症肆虐沿江城乡，长江两岸有些地方已经成为癌症高发区；物种受到威胁，珍稀水生物日益灭绝；固体废物污染严重，威胁水闸与电厂；湿地面积日益缩减，水的天然自洁功能日益丧失。综观上述：长江危机已经达到令人触目惊心的地步，因此治理保护长江的任务迫在眉睫。1.2问题提出进行长江水质评价和预测是致力保护长江的一个重要步骤。所谓的长江水质评价和预测是指通过物理或化学手段获取长江水环境检测数据，通过信息技术将这些检测数据转换为确定长江水环境状况的信息，获取长江水环境现状及其水质分布状况，分析长江现在存在的问题，抓主要矛盾，再预测其以后的发展趋势，制定综合防治措施与方案。现给出了统计出的关于长江流域的一系列检测数据以及国际水质标准的标限值，要求我们研究如下几个问题并对解决长江水质污染问题提出可行性建议。问题一：对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染情况。问题二：研究、分析长江干流近一年多主要污染物高锰酸钾指数和氨氮的污染源主要分布在哪些区域。问题三：假想如果不采取措施治理长江，根据所给出的过去10年内长江流域水质报告给出的统计数据，对长江未来10年的水质污染发展趋势做出预测分析。问题四：要求基于问题三的分析，在满足未来10年内年内江干流的IV类和V类水的比例都控制在20%以内，且没有劣V类水，求出每年需要处理污水的吨数。1.3 研究意义我们现在看到的情况是这样的：长江好像患了早期癌症，如果我们不及时治理，很快就会发展为晚期癌症，等到它真的重蹈黄河、淮河覆辙再公之于众，就晚了。虽说网上公布了许多关于长江现状的数据，但那些数据都是零散的，抽象的，普通人在其中不能得到有用的信息，对长江的现状还是很漠然，就要我们通过有效地数据处理，运用适当的数学方法，将零散的数据转化为具体的文字和图片，让更多的人产生危机意识，让更多的大老板能适当的停下手中的机器，呼吁更多的人参与到保护长江的行动中来，这也是一件造福我们子孙后代的有意义的事。2.模型假设1. 主要污染物高锰酸盐指数和氨氮的降解系数取0.2(单位：1/天)。2. 两观测点之间江水的流速是所测的两点流速的平均值。3. 不考虑洪水干旱等特殊气候对水质的影响。4.3.符号说明 评价对象（17个城市） 评价参数 (DO、CODMn、NH3-H、PH) 评价等级 参数权重 降解系数：第i个观测点的监测污染浓度。第i个观测点污水流至第i+1个观测点时的污染浓度：第i个观测点的监测水流量：第i个观测点的排污量：第i个观测点与第i+1个观测点之间的河长。：第i个观测点与第i+1个观测点水流速度的平均值。：第i个观测点到第i+1个观测点水流所需的时间 每年总流量 每年干流占总河长的比例 四、五类水占干流的比例 六类水占干流的比例 四、五类污水处理量 六类污水处理量 污水处理总量 干流流量4.模型建立与求解4.1.1问题一分析问题一要求对长江水质做出定量综合评价，并分析各地水质情况。从题中可以看出，溶解氧（OD）、高锰酸盐指数（CODMn）、铵盐（NH3-N）、只要有一个为高类别，则水质等级为高类别，所以评价指标存在“质差”、与“量差”的关系，在确定综合评价指标时，既要体现不同类型指标的差异，也要体现同类型指标的数量差异。采用动态加权综合评价方法可以有效地解决这一问题。4.1.2问题一模型建立模型建立分三步走：数据标准化确定权值建立综合评价函数。（一）数据标准化假设17个城市为评价对象S1、S2、S17，共四项评价指标：DO、CODMn、NH3-H和PH值，分别记做,前三项指标有6个等级相应分类区间如下表所示:指标I类II类III类IV类V类劣V类溶解氧（DO)7.5,)6,7.5)5,6)3,5)2,3)0,2高锰酸盐指数(CODMn)(0,2(2,4(4,6(6,10(10,15(15,)氨氮（NH3-N）(0,0.15（0.15，0.5(0.5,1(1,1.5(1.5,2(2,)PH值（无量纲）6,9（1）溶解氧（DO）的标准化注意到溶解氧(DO)为极大型指标，首先将数据指标做极小化处理，即令倒数变换，相应的分类标准区间变为，然后通过极差变换将其数据标准化，对应的分类区间随之变为（2）高锰酸盐指数的标准化高锰酸盐指数本身就是极小型指标，即由极差变换将其数据标准化，即令，对应的分类区间随之变为（3）氨氮的标准化氨氮也是极小型指标，对指标数据做极差变换将其数据标准化，即令，对应的区间随之变为（4）PH值的处理酸碱度的大小反映出水质呈酸碱性的程度，通常的水生物都适应于中性水质，即酸碱度的平衡值（PH值略大于7），在这里不妨取正常值的中值7.5.当PH7.5时偏酸性，而偏离值越大水质就越坏，PH值属于中间型指标。为此，对所有的PH值指标数据做均值差处理，即令，则将其数据标准化。（二）确定权值取动态加权函数为偏大型正态分布函数，即其中在这里取指标的I类水标准区间的中值，即，由确定。代入数据计算得：（三）建立综合评价函数前三个参数分配权值为0.75，由第四个参数PH值的特殊性取定权值为0.25.则某城市某一时间的水质综合评价指标为：经计算可得各城市水质综合评价指标值，即可得到一个阶的综合评价矩阵。4.1.3模型一求解：根据计算结果做出月均综合评价走势图表，如下：月份306307308309310311312401402403404405406407平均评价值0.32 0.34 0.350.32 0.32 0.55 0.64 0.67 1.05 0.53 0.44 0.32 0.32 0.30 月份408409410411412501502503504505506507508509平均评价值0.30 0.24 0.36 0.33 0.48 0.37 0.31 0.23 0.25 0.26 0.27 0.28 0.27 0.27 从图中可看出第9个月份即2004年2月污染最严重，整体污染呈下降趋势。排序结果如下表：城市/排序S1S2S3S4S5S5S7S8S9S10S11S12S13S14S15S16S17平均评价值0.20 0.15 0.17 0.21 0.15 0.17 0.15 0.75 0.28 0.52 0.11 0.50 0.31 0.24 2.07 0.26 0.23 Borda数263320300217343297302672201813919316218210233224总排序7351126416101311514121789从表中可以看出,水质最差的是S15江西南昌滁槎 赣江（鄱阳湖入口）,其次是S8四川乐山岷江大桥 岷江（与大渡河汇合前）。水质最好的是S11湖北丹江口 胡家岭丹江口水库（库体），排在第二的是S5江西九江河西水厂 干流（鄂-赣省界）。4.2.1问题二分析：下一个地方的污染物量由上游排污量与本地排污量组成，所以本地排污量（一个地方到下一个地方的排污量）为下一个地区的污染量减去上游排污量（本地区污染量到下一地区净化后的量）。所以可以建立水质依靠流量、流速和降解系数的数学模型，从而算出长江干流沿岸各个地段的排污量。4.2.2问题二建模：查找河水自净能力的参考文献得：其中k=0.2/天=长江干流各段排污量为：4.2.3模型求解求解过程中注意单位的换算，结果为：2004年4月2005年4月各地CODMn排出量m（单位：kg/s）地点/时间四川攀枝花重庆朱沱湖北宜昌湖南岳阳江西九江安徽安庆江苏南京2004.04 44.324242.844084.283681.348482.425144.5811/2004.05 21.526757.911073.616892.2097108.537961.9434/2004.06 37.783573.368077.586975.665750.320365.8000/2004.07 48.494766.412198.769568.432753.136757.2131/2004.08 10.906772.2439101.376272.256862.150558.9114/2004.09 184.8731120.7294197.9590128.2419105.4752173.3202/2004.10 24.765067.512675.083245.365064.155671.9649/2004.11 14.725619.386731.171332.095337.375235.2664/2004.12 7.930117.535637.424730.341123.960429.7966/2005.01 4.97428.587341.767926.156331.977435.9986/2005.02 6.35838.998625.535531.927836.986834.7284/2005.03 8.863510.858928.863425.735157.771938.6917/2005.04 7.113712.910820.995830.199246.454533.1462/平均值32.510744.561568.802656.921158.517557.0279/2004年4月2005年4月各地NH3-N排出量m（单位：kg/s）地点/时间四川攀枝花重庆朱沱湖北宜昌湖南岳阳江西九江安徽安庆江苏南京2004.04 2.63875.46866.63197.57586.47370.5867/2004.05 3.41835.75736.54176.54098.48143.0928/2004.06 2.47896.67577.77464.32165.39701.9435/2004.07 2.85823.29238.55354.16856.95071.5316/2004.08 0.02734.37678.41077.57425.20741.6911/2004.09 7.253513.290615.128317.42488.64535.3443/2004.10 3.61053.93917.84364.53155.17561.4554/2004.11 2.17132.44234.31643.20862.03750.9793/2004.12 3.31131.09603.31531.58231.97351.4897/2005.01 2.15610.96563.43972.28872.45983.1679/2005.02 1.96050.80533.43141.33872.46594.6809/2005.03 2.52590.97882.95671.57254.92150.4293/2005.04 2.28260.68653.04082.11393.02964.4198/平均值2.82253.82886.26034.94174.86302.3702/从表中可以看出高锰酸盐污染源主要在湖北宜昌到湖南岳阳段内；氨氮污染源主要也在湖北宜昌到湖南岳阳段内。5.3.1模型的准备灰色预测模型：灰色预测是一种对含有不确定因素的系统进行预测的方法。灰色预测通过鉴别系统因素之间发展趋势的相异程度，即进行关联分析，并对原始数据进行生成处理来寻找系统变动的规律，生成有较强规律性的数据序列，然后建立相应的微分方程，从而预测事物未来发展趋势的一种方法。本小题运用灰色模型中累加生成的方式生成了数列，构建了模型，依照过去10年关于长江流域的统计数据，对长江未来10年内的发展趋势做了预测分析。5.3.2数据的初始化处理1.对水质等级作出标准化处理得：2.设长江的污染程度为，各个等级水质的河长占评价河长的比例为3.对所给的年份作出标准化处理得：5.3.3 模型的建立1.分析可知，长江总的污染程度为：2.附件4给出了长江流域近10年来各个时期不同流域的水质情况，利用MATLAB软件编程绘图，得出如下的水质污染折线图：枯水期丰水期水文年注：红线全流域 绿线干流 蓝线支流单纯图上也可看出长江水质在逐年恶化。5.3.3.建立灰色模型预测水质污染的发展趋势：已知参考数据列为，将其做一次累加（AGO）,将得到的生成数列记为，则有=。其中取的均值数列：，则有，于是我们建立起了灰微分方程如下：，相应的白化微分方程为。（1）5.3.4 模型的求解记,则由最小二乘法可知，要求得使达到最小值的的条件为。求解方程（1）得：为了确保建模方法的可行性，还需要对已知的数据做必要的检验处理。下面假设计算数列的级比： ，k=2，3，n均落在可容覆盖内。其次，再建立模型GM（1，1），则可得到预测值如下：且有，经数据验证，得出未来在未来10年内，污染会越来越严重。5.4模型四的建立与求解5.4.1 模型四的建立问题四要求基于问题三中分析的关于长江未来水质污染发展的趋势，并在确保未来10年内每年干流的IV类和V类水的比例在20%以内，且没有劣V类水的条件下，求解出每年需要处理污水的吨数。1.利用MATLAB软件编程绘图，得出如下折线图：c1与c2干流的流量干流占总河长的