下载文档原格式
摘要问题一主要是对长江水质的综合评价和对各地水质的分析，在对长江水质的综合评价中，我们运用了标准值法，分别算出各类污染物的标准值，然后加权求和，以此判断水质类别，通过对长江水质近两年各地水质类别和个污染物浓度的处理，我们得到，长江水质处于II 类水质类别，但是特别接近III 类水质类别。
在对各地分析中，我们做出了近两年每月各地水质类别变化的折线图，图见正文，并算出各地水质平均值，得出结论湖北丹江口胡家岭属于I 水质类别，四川攀枝花，重庆朱沱，等地属于II 类水质，四川宜宾凉姜沟，四川泸州沱江二桥，等地属于III 类水质类别，四川乐山岷江大桥属于IV 水质类别，江西南昌滁槎属于V 水质类别。
问题二中主要统计了高锰酸盐和氨氮化物引起重大污染的次数以及发生的地区，并计算出各地区发生次数占总次数的比重，图见正文，得到氨氮化物污染源在江西南昌,湖南长沙，和四川攀枝花。
高锰酸盐污染源在湖南岳阳，江西九江。
问题三中主要解决的是未来十年长江水质的预测，我们用了回归模型，利用过去十年数据，分别利用Matlab 拟合了废水量--时间，废水浓度--Ⅰ类水百分比，Ⅱ类水百分比，Ⅲ类水百分比，Ⅳ类水百分比，Ⅴ类水百分比，劣Ⅴ类水百分比的曲线，预测了未来十年各类别水所占百分比，结果见正文表（3-3）问题四中主要解决要使Ⅳ类和Ⅴ类水控制在20%之内，消除劣Ⅴ类水，即使非可饮用水控制在20%之内，我们利用第三问的回归模型得到非可饮用水在20%时，各年的废水浓度从而得到每年需处理的废水量 问题五中我们根据前四问的模型对提高长江水质提出了相关建议。
关键词：标准值法 回归模型 长江水质年份2005 2006 2007 2008 2009 2010 2011 2012 2013 2014需处理废水单位（亿吨）143.39 166.29 190.79 216.99 244.89 274.49 305.69 338.69 373.19 409.49问题重述水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重。
2005高教社杯全国大学生数学建模竞赛题目A题: 长江水质的评价和预测水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重。
专家们呼吁：“以人为本，建设文明和谐社会，改善人与自然的环境，减少污染。
”长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。
2004年10月，由全国政协与中国发展研究院联合组成“保护长江万里行”考察团，从长江上游宜宾到下游上海，对沿线21个重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心。
为此，专家们提出“若不及时拯救，长江生态10年内将濒临崩溃”（附件１），并发出了“拿什么拯救癌变长江”的呼唤（附件2）。
附件3给出了长江沿线17个观测站（地区）近两年多主要水质指标的检测数据，以及干流上７个观测站近一年多的基本数据（站点距离、水流量和水流速）。
通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。
一般说来，江河自身对污染物都有一定的自然净化能力，即污染物在水环境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低。
反映江河自然净化能力的指标称为降解系数。
事实上，长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于0.1~0.5之间，比如可以考虑取0.2(单位：1/天)。
附件4是“1995~2004年长江流域水质报告”给出的主要统计数据。
下面的附表是国标(GB3838-2002)给出的《地表水环境质量标准》中4个主要项目标准限值，其中Ⅰ、Ⅱ、Ⅲ类为可饮用水。
请你们研究下列问题：（1）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。
（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?（3）假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况。
长江水质的评价和预测李云锋王勇...本文利用长江流域近两年多主要城市水质检测数据，通过对原始数据进行归一化综合处理，确定了水质新的综合评判指标函数ψ。
在对整个长江流域所有观测站的位置关系作一定的简化假设后，得到长江综合评定函数值ψ=0．4331，水质为良好。
主要污染物为氨氮。
通过建立污染浓度的反应扩散方程，本文用三种方法反演出未知的污染源强迫函数f（x，t），并对，（x，t）的三种数据加以综合分析，分别给出了高锰酸钾盐和氨氮污染源的主要分布地区。
为了对长江未来水质污染发展趋势进行预测，本文建立了回归分析模型并对回归系数进行了F检验，结果是如果不采取有效的治理措施。
长江可饮用水将逐年下降，且10年后可饮用水所占长江水总量的比例将不到50％。
根据这一预测结果，我们进而使用二元线性回归模型。
通过对各种不可饮用水进行综合考虑，得到如下结果：要在未来10年内使长江干流的不可饮用水（IV类和V类水）的比例控制在20％以内，且没有劣V 类水，那么每年污水处理量至少为75．195亿吨长江水质的评价和预测.pdf (370.52 KB)水质的评价和预测模型张震张超...本文首先考虑到水质类别的差异和相同类别水质在数量上的差异对综合评价的影响。
构造“S”形的变权函数，对属于不同水质类别的同种污染指标进行“动态加权”，建立基于逼近理想点排序法的评价模型和利用灰色关联度的分析方法。
对长江水质状况做出了综合评价：其次，根据7个观测站的位置将干流分成8段，把每段河道内所有污染源都等效为一个段中央的连续稳定源，分别利用稳态条件下的一维水质模型及质量守恒定律。
得出中间6段每个月的排污量，综合比较各河段一年多来的总排污量得到主要污染源的分布区域：然后，用每年不可饮用类水的百分比之和刻画水质状况。
综合利用灰色GM（1，1）模型和时间序列分析方法，对变化趋势进行了预测：最后，建立不可饮用类水的百分比与长江水总流量和废水排放量的线性回归模型，计算在满足约束条件下排污量的极限值，用排污量的预测值减去极限值，得到未来10年的污水处理量水质的评价和预测模型.pdf (283.07 KB)长江水质的评价预测模型谯程骏张东辉...本问题是一个对长江的水质进行综合评价、预测和控制的问题。
承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们将受到严肃处理。
我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名) ：1.2.3.指导教师或指导教师组负责人(打印并签名)：日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：长江水质的评价和预测摘要针对第一问首先对本题所给的数据进行分析，发现长江水质中各个污染物指标较多为了简化解题过程，比较简明地分析长江流域各个监测站水质情况，对长江流域整体水质作分析我们运用主成分分析法来减少变量，提取主成分，使其更容易得出这十七个观测站的总排名得出了江西、四川的一些地区污染比较严重。
对于第二问我们得出长江干流主要污染物CODMn和NH3-N的平均排放量第三问先运用灰色预测GM（1，1）来预测未来十年水污染的趋势，这是因为水质问题是一个复杂的非线性系统，但是由于数据样本少，需要预测的时间长。
由此我们得出未来十年长江的排污量。
第四问我们根据第三问得出的数据结合BP神经网络可以进一步算出每年需要处理的排污量。
关键词：主成分分析GM（1，1）BP神经网络问题重述作为中华民族最具象征性的的河流，长江哺育了灿烂悠久的中华文明，她是我国赖以生存的生命线，北方的干旱需要他的乳汁哺育因此有了南水北调，经济的发展所需的电力由她提供、因此有了三峡工程，可是就是这样一条对我们生存和发展都有重大影响的母亲河正在一步步走向没落---长江的污染程度日益加重，我们要积极对长江进行调查研究，通过科学的研究分析来拯救我们的生命线。
长江水质的评价和预测摘 要:河流污染物浓度一维稳态衰减规律:uKxeC C -⋅=0对于点源，河水和污水的稀释混合方程为：Ep EE p p Q Q Q C Q C C +⋅+⋅=排污口允许纳污量（水环境容量）计算公式为：pp Ei ni p C C Q Q Q S W ⋅-+⋅=∑=)(1四川攀枝花、湖南岳阳和江苏扬州，这三个地区高锰酸盐的浓度常年比较高，是污染物高锰酸盐的主要排放地。
重庆朱沱、岳阳城陵矶、江苏扬州，这3个地区氨氮（NH3-N ）的浓度常年比较高，是污染物氨氮（NH3-N ）的主要排放地。
近几年来，长江流域水质的组成为Ⅰ类占总流域长度的2%、Ⅱ类水占总流域长度的27%、Ⅲ类水占总流域长度的39%、Ⅳ类水占总流域长度的15%、Ⅴ类水占总流域长度的7%、劣Ⅴ类水占总流域长度的10%。
其中Ⅰ类水由十年前的15%变化为目前的2%，呈减小趋势。
劣V 类水由3%变化为10%，呈增大趋势。
可饮用水在长江中的比例在逐年的变小，而Ⅳ类水以后的水在逐年的变大。
全干流域IV 类水、V 类水和劣V 类水所占百分比的预测模型分别为：IVF= 4.3-2.18* t + 0.06*t 2 +10.5* Ln(t) ， VIF= 2.3 +0.5* t + 0.005* t 2 -1.14*Ln(t) ， VIF= 2.3 +0.5* t + 0.005* t 2 -1.14*Ln(t) ，式中t 为时间（年）。
从1995年到2004年这10年间废水排放总量用最小二乘法数据模拟，可得变化规律为：QF ＝167.375＋3.68* t +0.835* t 2 ，式中为t 时间（年）。
可预测今后10废水排放总量为：308.9,331.8,356.3,382.5,410.4,440.0,471.3,504.2,538.7,575.0亿吨。
如果未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,每年需要在现有处理的基础增加污水处理数量：24，47，71，98，125，155，186，219，254，295亿吨。
一：问题的提出A题: 长江水质的评价和预测水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重。
专家们呼吁：“以人为本，建设文明和谐社会，改善人与自然的环境，减少污染。
”长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。
2004年10月，由全国政协与中国发展研究院联合组成“保护长江万里行”考察团，从长江上游宜宾到下游上海，对沿线21个重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心。
为此，专家们提出“若不及时拯救，长江生态10年内将濒临崩溃”（附件１），并发出了“拿什么拯救癌变长江”的呼唤（附件2）。
附件3给出了长江沿线17个观测站（地区）近两年多主要水质指标的检测数据，以及干流上７个观测站近一年多的基本数据（站点距离、水流量和水流速）。
通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。
一般说来，江河自身对污染物都有一定的自然净化能力，即污染物在水环境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低。
反映江河自然净化能力的指标称为降解系数。
事实上，长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于0.1~0.5之间，比如可以考虑取0.2(单位：1/天)。
附件4是“1995~2004年长江流域水质报告”给出的主要统计数据。
下面的附表是国标(GB3838-2002)给出的《地表水环境质量标准》中4个主要项目标准限值，其中Ⅰ、Ⅱ、Ⅲ类为可饮用水。
请你们研究下列问题：（1）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。
（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?（3）假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况。
（4）根据你的预测分析，如果未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,那么每年需要处理多少污水？ （5）你对解决长江水质污染问题有什么切实可行的建议和意见。
二、模型假设1．主要污染物高锰酸盐指数和氨氮的降解系数通常介于0.1~0.5之间，不考虑时期、温度等影响，取0.2 (单位：1/天)； 2． 水流量和水流速以恒定值流动；3． 污染物在较短的时间内基本能混合均匀，即在排污口断面瞬时完成均匀混合； 4． 假定污染物浓度仅在河流纵向上发生变化，在横向和竖向的浓度梯度为零 5．增加和减少的水量只在河流上排污口、取水口、干支流汇合口等节点处发生；三、一维水质数学模型对于河流而言，一维模型假定污染物浓度仅在河流纵向上发生变化，主要适用于同时满足以下条件的河段：1）宽浅河段；2）污染物在较短的时间内基本能混合均匀；3）污染物浓度在断面横向方向变化不大，横向和垂向的污染物浓度梯度可以忽略。
如果污染物进入水域后，在一定范围内经过平流输移、纵向离散和横向混合后达到充分混合，或者根据水质管理的精度要求允许不考虑混合过程而假定在排污口断面瞬时完成均匀混合，即假定水体内在某一断面处或某一区域之外实现均匀混合，则不论水体属于江、河、湖、库的任一类，均可按一维问题概化计算条件。
在忽略离散作用时，描述河流污染物浓度一维稳态衰减规律的微分方程为：Kc dx dcu-=将dt dxu =代入，得到 Kcdt dc -=积分解得 uKxeC C -⋅=0式中：u 是河流断面平均流速，单位为m/s ； x 是沿程距离，单位为km ；K 是综合降解系数，单位为1/d ；C 是沿程污染物浓度，单位为mg/L ；C 0是前一个节点后污染物浓度，单位为mg/L 。
四、定常设计条件下河流稀释混合模型对于点源，河水和污水的稀释混合方程为：Ep EE p p Q Q Q C Q C C +⋅+⋅=式中：C ——完全混合的水质浓度（mg/L ）；Q p 、 C p ——上游来水设计水量（m 3/s ）与设计水质浓度（mg/L ）；Q E 、 C E ——污水设计流量（m 3/s ）与设计排放浓度（mg/L ）； 由于污染源作用可线性迭加，多个污染源排放对控制点或控制断面的影响，等于各个污染源单个影响作用之和，符合线性迭加关系。
单点源计算可迭加使用，计算多点源条件。
单断面或单点约束条件，可根据节点平衡，递推多断面或多点约束条件。
对于可概化为完全均匀混合类的排污情况，排污口与控制断面之间水域的允许纳污量计算公式为：单点源排放：pp E p C C Q Q Q S W ⋅-+⋅=)(式中：W C ——水域允许纳污量(g/L)；S ——控制断面水质标准(mg/L)多点源排放： pp Ei ni p C C Q Q Q S W ⋅-+⋅=∑=)(1式中：Q Ei ——第i 个排污口污水设计排放流量(m 3/s)； n ——排污口个数。
五、水环境容量（1）河流概化污染物进入河流后，在一定范围内经过平流输移、纵向离散和横向混合后达到充分混合，或者根据水质管理的精度要求允许不考虑混合过程而假定在排污口断面瞬时完成均匀混合，可按一维问题概化计算条件，建立水质模型。
河流一维水质模型由河段和节点两部分组成，节点指河流上排污口、取水口、干支流汇合口等造成河道流量发生突变的点，水量与污染物在节点前后满足物质平衡规律（忽略混合过程中物质变化的化学和生物影响）。
河段指河流被节点分成的若干段，每个河段内污染物的自净规律符合一阶反应规律。
图 河流一维模型概化示意图上界j#k#如图4所示，假定功能区内有i 个节点，则将河流分成i ＋1个河段。
在节点处，要利用节点均匀混合模型进行节点前后的物质守恒分析，确定节点后的河段流量和污染物浓度。
节点后的河段要以节点平衡后的流量和污染物浓度为初始条件，按照一维降解规律计算到下一个节点前的污染物浓度。
（2）节点平衡方程。
考虑干流、支流、取水口、排污口均在同一节点的最复杂情况，水量平衡方程为：Q 干流混合后＝Q 干流混合前＋Q 支流＋Q 排污口－Q 取水口污染物平衡方程为（忽略混合过程的不均匀性）：取水口排污口支流干流混和前取水口取水口排污口排污口支流支流干流混和前干流混和前干流混合后＝Q Q Q Q Q *C Q *C Q *C Q *C C ++++++（3）环境容量计算将 ji i Q Q Wi C C +=54.31/＋代入模型，得到一维模型水环境容量的计算公式为：)(*)*(*54.31*4.86j i i uKx Q Q C eC Wi +-=式中：Wi ——第i 个排污口允许排放量，t/a ；C i ——河段第i 个节点处的水质本底浓度，mg/l ；C ——沿程浓度，mg/l ；Q i ——河道节点后流量，m 3/s ；Q j ——第i 节点处废水入河量，m 3/s ；u ——第i 个河段的设计流速，m/s ；x ——计算点到第i 节点的距离，m 。
若以下界处作为功能区考核的断面，按照上述方法沿程计算整个功能区的沿程污染物浓度变化规律，如图5所示。
若功能区要求是Ⅲ类，COD Ⅲ类标准是20mg/l ，则计算结果模拟结果显然超过Ⅲ类水质标准要求，就要通过削减每一个排污口的排污量来重新计算，直到计算结果满足水质标准要求为止。
这时各个排污口的排污量之和，即∑=ni WiW 就是此环境功能区内的一个水环境容量值。
六、模型计算6.1、对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况 主要污染物指标变化特征注：1、各地区各值的折线图形变化情况分为以下五种情况：①平稳，②上升，③下降，④波动(振幅较平缓)，⑤振荡（振幅较剧烈）;2、表中所列数为最小值和最大值；3、水质类别为本地区28个月汇总数据，数值越小水质越好；其中将Ⅰ类变为数字1，将Ⅱ类变为数字2，将Ⅲ类变为数字3，将Ⅳ类变为数字4，将Ⅴ类变为数字5，将劣Ⅴ类变为数字6，然后逐个相加而得。
由上述分析可得，长江水所受的污染在随时间而加剧，水值随时间而变坏。
所以，治理长江迫在眉捷！6.2 长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要地区 （1）、长江流域高锰酸盐（CODMn ）浓度变化情况对长江沿线的十一个城市进行研究和分析，得到如上高锰酸盐（CODMn ）浓度的变化情况图，它们分别表示在枯水期、平水期和丰水期高锰酸盐（CODMn ）浓度随地域的不同而发生变化的情况。
由图可见，四川攀枝花、湖南岳阳和江苏扬州，这三个地区高锰酸盐的浓度常年比较高，是污染物高锰酸盐的主要排放地。
