长江水质预测和评价.doc

长江水质预测和评价 李杰,,张峻,杨晓青 （西安电子科技大学 西安） 指导老师:朱强 刘丽霞 摘要:利用对评价分析问题有着普遍适用性的集对分析法对长江总体和各地区进行了评价，并预测了未来10年水质的变化情况，得到了长江水质上游好于下游，总体逐年变坏的结论。对于第二个问题建立了基于质量守恒的平衡方程，忽略次要因素，应用数值积分和微分方程，对长江各地排污量进行计算，找到了主要污染物的来源。应用拟合的方法得出了未来10年每年至少需要处理的污水量。最后本文给出了评价和预测结果，通过这些结果给出了符合实际并且可行的意见。 一.问题提出 水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重。专家们呼吁：“以人为本，建设文明和谐社会，改善人与自然的环境，减少污染。” 长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。2004年10月，由全国政协与中国发展研究院联合组成“保护长江万里行”考察团，从长江上游宜宾到下游上海，对沿线21个重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心. 为此，专家们提出“若不及时拯救，长江生态10年内将濒临崩溃”（附件１），并发出了“拿什么拯救癌变长江”的呼唤（附件2）。 附件3给出了长江沿线17个观测站（地区）近两年多主要水质指标的检测数据，以及干流上７个观测站近一年多的基本数据（站点距离、水流量和水流速）。通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。一般说来，江河自身对污染物都有一定的自然净化能力，即污染物在水环境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低。反映江河自然净化能力的指标称为降解系数。事实上，长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于0.1~0.5之间，比如可以考虑取0.2　(单位：1/天)。附件4是“1995~2004年长江流域水质报告”给出的主要统计数据。下面的附表是国标(GB3838-2002) 给出的《地表水环境质量标准项目标准限值。NH3-N）NH3-N）：长江总体水质联系度； N：长江上观测点总数（N＝17）； S：N个观测点中属于I级的个数； P：N个观测点中属于劣V：N个观测点中属于i＋1级的个数； M：观测的总月数（M＝28）； ：长江上第i个观测点水质的联系度； S’：M个月中水质属于I级的月数； P’：M个月中水质属于劣V：M个月中水质属于i－1级的月数； K：1995～2004年年份的编号（K＝1～20）； i:第i年 ：第i年长江的总的排污量； ：第i年长江的总流量； ：分别为第四、五及劣五类水第i年的河长百分比； ：表示点单位时间内污染物总量 ：表示前方(-1点以前)积累的污染物单位时间内到达点的总量； ：第点单位时间内排放的污染物质量 ：表示第点的污染物浓度 ：表示第点的水流量 ： ：降解系数 五、模型建立与求解 1、问题一模型与求解 1.1 集对分析方法简介 集对分析的基本思路是：在一定的问题背景下对所论的两个集合所具有的特性作同、异、反分析并加以度量刻划，得出这两个集合在所论问题背景下的同、异、反联系度表达式μ=a＋bi＋cj，μ这里称为联系数，对于一个具体问题即是联系度的概念，a为同一度，b为差异度，c为对立度，i为差异度系数，在[－1,1]中取值，j为对立系数，运算时恒取－1，在－1→1之间变化，体现了确定性与不确定性之间的相互转换，随着i→0，不确定性明显增加，而i取－1与1，都是确定性的，这是问题的实质。将集对分析方法用于水评价，可视待评价地区的水的某项指标和其标准分别为两个集合，将这两个集合构成一个集对，若该指标处于评价级别中，则认为是同一；若处于相隔的评价级别中，则认为是对立，若指标在相邻的评价级别中，则认为是异；取差异系数i在[－1,1]之间变化，越接近所要评价的级别，i越接近l，越接近相隔的评价级别，i越接近－l。 其中，，， ，j＝－1。 利用集对分析法对总体长江水质进行评价。首先将长江各段构建一个集对，对长江上N个观测点，设有S个属于I级，P个属于劣V、、 、 分别表示属于II、III、Ⅳ、Ⅴ级的个数。则联系度表达式为： 利用集对分析法对长江各观测点水质进行评价。对M个月，设S’表示属于I级的月数，P’表示属于劣V表示属于j级的月数（j＝II、III、Ⅳ、Ⅴ）（i表示第i各观测点）。则各观测点的联系度表达式为： 根据实际情况我们把长江水质划分为5个等级： 水质分别为：优、良、中、差、非常差 1.3模型求解 在求解过程中i我们分别取范围的中间值，即： ，，，。 经统计计算可得下表 长江个月总体关联度