《科学计算与数学建模第九章(精品).ppt》由会员分享，可在线阅读，更多相关《科学计算与数学建模第九章(精品).ppt（72页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、科学计算与数学建模科学计算与数学建模中南大学数学科学与计算技术学院中南大学数学科学与计算技术学院 综合评价方法及其应用综合评价方法及其应用第九章 长江水质的综合评价 综合评价方法及其应用长江水质的综合评价模型长江水质的综合评价模型1综合评价方法的基本概念综合评价方法的基本概念2评价指标的规范化处理评价指标的规范化处理3456污染源的确定污染源的确定水质综合评价模型水质综合评价模型 1 长江水质的综合评价模型9.1.19.1.1长江水质的评价和预测（长江水质的评价和预测（20052005年大学数学建模年大学数学建模A A题）题）1、问题、问题 水是人类赖以生存的资源，保护水资源就是保护我们自己，
2、对于水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重。专家们呼吁：我国大江大河水资源的保护和治理应是重中之重。专家们呼吁：“以以人为本，建设文明和谐社会，改善人与自然的环境，减少污染。人为本，建设文明和谐社会，改善人与自然的环境，减少污染。”长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。已引起了相关政府部门和专家们的高度重视。2004年年10月，由全国政月，由全国政协与中国发展研究院联合组成协与中国发展研究院联合组成“保护长江万里行
3、保护长江万里行”考察团，从长江上考察团，从长江上游宜宾到下游上海，对沿线游宜宾到下游上海，对沿线21个重点城市做了实地考察，揭示了一幅个重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心。为此，专家们提出长江污染的真实画面，其污染程度让人触目惊心。为此，专家们提出“若不及时拯救，长江生态若不及时拯救，长江生态10年内将濒临崩溃年内将濒临崩溃”（附件），并发出（附件），并发出了了“拿什么拯救癌变长江拿什么拯救癌变长江”的呼唤（附件的呼唤（附件2）。）。附件附件3给出了长江沿线给出了长江沿线17个观测站（地区）近两年多主要水质指个观测站（地区）近两年多主要水质指标的检测数据
4、，以及干流上个观测站近一年多的基本数据（站点标的检测数据，以及干流上个观测站近一年多的基本数据（站点 距离、水流量和水流速）。通常认为一个观测站（地区）距离、水流量和水流速）。通常认为一个观测站（地区）的水质污染主要来自于本地区的排污和上游的污水。一般的水质污染主要来自于本地区的排污和上游的污水。一般说来，江河自身对污染物都有一定的自然净化能力，即污说来，江河自身对污染物都有一定的自然净化能力，即污染物在水环境中通过物理降解、化学降解和生物降解等使染物在水环境中通过物理降解、化学降解和生物降解等使水中污染物的浓度降低。反映江河自然净化能力的指标称水中污染物的浓度降低。反映江河自然净化能力的指标
5、称为降解系数。事实上，长江干流的自然净化能力可以认为为降解系数。事实上，长江干流的自然净化能力可以认为是近似均匀的，根据检测可知，主要污染物高锰酸盐指数是近似均匀的，根据检测可知，主要污染物高锰酸盐指数和氨氮的降解系数通常介于和氨氮的降解系数通常介于0.10.50.10.5之间，比如可以考虑之间，比如可以考虑取取0.20.2(单位：单位：1/1/天天)。附件。附件4 4是是“1995200419952004年长江流域年长江流域水质报告水质报告”给出的主要统计数据。下面的附表是国标给出的主要统计数据。下面的附表是国标(GB3838-2002)(GB3838-2002)给出的给出的地表水环境质量标
6、准地表水环境质量标准中中4 4个主个主要项目标准限值，其中要项目标准限值，其中、类为可饮用水。类为可饮用水。请你们研究下列问题：请你们研究下列问题：（1）对长江近两年多的水质情况做出定量的综合）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况评价，并分析各地区水质的污染状况（2）研究、分析长江干流近一年多主要污染物高）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区锰酸盐指数和氨氮的污染源主要在哪些地区?序序号号类类 类类 类类 类类 类类劣劣类类 17.5（或饱（或饱和率和率90%）65 3 2 0 2高锰酸盐指高锰酸盐指数数(CODMn)2
7、 4610 15 3氨氮氨氮（NH3-N）0.15 0.5 1.0 1.5 2.0 4PH值（无值（无量纲量纲）6-9 分类分类标准值标准值项目项目溶解溶解(DO)附表附表:地表水环境质量标准地表水环境质量标准（GB3838GB383820022002）中）中4 4个主要项目标准限值个主要项目标准限值 单位：单位：mg/Lmg/L2、问题说明、问题说明 针对长江水质的综合评价这一问题，采用动态加权综合评价方法来解针对长江水质的综合评价这一问题，采用动态加权综合评价方法来解决。假设决。假设1717个城市为被评价对象个城市为被评价对象 ，共有四项评价指标（或，共有四项评价指标（或属性）属性）DOD
8、O、CODMnCODMn、NH3-N NH3-N 和和PHPH值，分别记为值，分别记为 和和 ，前三项指标都，前三项指标都有有6 6个等级个等级 ，相应的分类区间值如表（，相应的分类区间值如表（1 1）所示，而）所示，而PHPH值没有值没有等级之分。等级之分。地表水环境质量标准地表水环境质量标准(GB3838(GB38382002)2002)中中4 4个主要项目标准限值个主要项目标准限值 mg/Lmg/L指 标 类类类类 类类 类类 类 劣劣类类 溶解氧溶解氧(DO)7.5,)6,7.5)5,6)3,5)2,3)0,2 高锰酸盐指数高锰酸盐指数(CODMn)(0,2(2,4(4,6(6,10(
9、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 2 综合评价方法的基本概念综合评价的问题综合评价的问题:对被评价对象所进行的客观、公正、合理的对被评价对象所进行的客观、公正、合理的全面评价。通常的综合评价问题都是有若干个同类的全面评价。通常的综合评价问题都是有若干个同类的被评价对象被评价对象(或系统或系统)，每个被评价对象往往都涉及到多个属性（或指标）。，每个被评价对象往往都涉及到多个属性（或指标）。综合评价的目的综合评价的目的:根据系统的属性判断确定这些系统的运行根据系统的属性判断确定这些
10、系统的运行（或发展）状况哪个优，哪个劣，即按优劣对各被评价对象进排序（或发展）状况哪个优，哪个劣，即按优劣对各被评价对象进排序或分类或分类。这类问题又称为这类问题又称为多属性（或多指标）的综合评价问题多属性（或多指标）的综合评价问题。综合评价的应用综合评价的应用:研究多目标决策问题的前提，因此研究解决研究多目标决策问题的前提，因此研究解决这类问题在实际中是很有意义的，特别是在政治、经济、社会及军这类问题在实际中是很有意义的，特别是在政治、经济、社会及军事管理、工程技术及科学决策等领域都有重要的应用价值。事管理、工程技术及科学决策等领域都有重要的应用价值。构成综合评价问题的五个要素分别为构成综合
11、评价问题的五个要素分别为:被评价对象、评价指标、被评价对象、评价指标、权重系数、综合评价模型和评价者。权重系数、综合评价模型和评价者。被评价对象就是综合评价问题中所研究的对象，或称为系统。通被评价对象就是综合评价问题中所研究的对象，或称为系统。通常情况下，在一个问题中被评价对象是属于同一类的，且个数要大于常情况下，在一个问题中被评价对象是属于同一类的，且个数要大于1 1，不妨假设一个综合评价问题中有，不妨假设一个综合评价问题中有 个被评价对象（或系统），分别个被评价对象（或系统），分别记为记为。9.2.29.2.2评价指标评价指标评价指标评价指标是反映被评价对象是反映被评价对象(或系统或系统)
12、的运行的运行(或发展或发展)状况的基本状况的基本要素。通常的问题都是有要素。通常的问题都是有多项指标多项指标构成，每一项指标都是从不同的侧构成，每一项指标都是从不同的侧面刻画系统所具有某种特征大小的一个度量。面刻画系统所具有某种特征大小的一个度量。一个综合评价问题的评价指标一般可用一个综合评价问题的评价指标一般可用一个向量一个向量表示，其中每一表示，其中每一个分量就是从一个侧面反映系统的状态，即称为个分量就是从一个侧面反映系统的状态，即称为综合评价的指标体系。综合评价的指标体系。9.2.19.2.1被评价对象被评价对象9.2.39.2.3权重系数权重系数 每一综合评价的问题都有相应的评价目的，
13、针对某种评价目每一综合评价的问题都有相应的评价目的，针对某种评价目的，各评价指标之间的相对重要性是不同的，评价指标之间的这种的，各评价指标之间的相对重要性是不同的，评价指标之间的这种相对重要性的大小可以用权重系数来刻画。如果用相对重要性的大小可以用权重系数来刻画。如果用 来表示评价来表示评价指标指标 的的权重系数权重系数，则应有则应有 且且 。评价指标体系应遵守的原则：系统性、科学性、可比性、可测性评价指标体系应遵守的原则：系统性、科学性、可比性、可测性（即可观测性）和独立性。这里不妨设系统有（即可观测性）和独立性。这里不妨设系统有 个评价指标（或属性）个评价指标（或属性），分别记为，分别记为
14、 ，即评价指标向量为，即评价指标向量为 。注意到注意到:当各被评价对象和评价指标值都确定以后，问题的综当各被评价对象和评价指标值都确定以后，问题的综合评价结果就完全依赖于权重系数的取值了，即权重系数确定的合合评价结果就完全依赖于权重系数的取值了，即权重系数确定的合理与否，直接关系到综合评价结果的可信度，甚至影响到最后决策理与否，直接关系到综合评价结果的可信度，甚至影响到最后决策的正确性。的正确性。9.2.49.2.4综合评价模型综合评价模型对于多指标（或多因素）的综合评价问题对于多指标（或多因素）的综合评价问题,就是要通过建立合就是要通过建立合适的综合评价数学模型将多个评价指标综合成为一个适的
15、综合评价数学模型将多个评价指标综合成为一个整体的综合评整体的综合评价指标价指标，作为综合评价的依据，从而得到相应的评价结果。，作为综合评价的依据，从而得到相应的评价结果。不妨假设不妨假设n n个被评价对象的个被评价对象的m m个评价指标向量为个评价指标向量为 指标权重向量为指标权重向量为 ，由此构造综合评价函数为，由此构造综合评价函数为 。如果已知各评价指标的如果已知各评价指标的 个观测值为个观测值为 5 5 ，其中，其中 ，则可以计算出各系统的综合，则可以计算出各系统的综合评价值评价值 ，。根据。根据 值的大小将这值的大小将这 个系统进行排序或分类，即得到综合评价结果。个系统进行排序或分类，
16、即得到综合评价结果。9.2.59.2.5评价者评价者评价者是直接参与评价的人，可以是某一个人，也可以是一个团评价者是直接参与评价的人，可以是某一个人，也可以是一个团体。对于评价目的选择、评价指标体系确定、评价模型的建立和权重体。对于评价目的选择、评价指标体系确定、评价模型的建立和权重系数的确定都与评价者有关。系数的确定都与评价者有关。综合评价的一般步骤：综合评价的一般步骤：明确评价目的；确定被评价对象；建立评价指明确评价目的；确定被评价对象；建立评价指标体系（包括评价指标的原始值、评价指标的若标体系（包括评价指标的原始值、评价指标的若干预处理等）；确定与各项评价指标相对应的权干预处理等）；确定
17、与各项评价指标相对应的权重系数；选择或构造综合评价模型；计算各系统重系数；选择或构造综合评价模型；计算各系统的综合评价值的综合评价值,并给出综合评价结果。并给出综合评价结果。9.3.1 9.3.1 评价指标类型的一致化评价指标类型的一致化一般说来，在评价指标一般说来，在评价指标中中可能包含有可能包含有“极大型极大型”指标、指标、“极小型极小型”指标、指标、“中间型中间型”指标和指标和“区间型区间型”指标。指标。极大型指标极大型指标:总是期望指标的取值越大越好；总是期望指标的取值越大越好；极小型指标极小型指标:总是期望指标的取值越小越好；总是期望指标的取值越小越好；中间型指标中间型指标:总是期望
18、指标的取值既不要太大，也不要太小为好，即取总是期望指标的取值既不要太大，也不要太小为好，即取适当的中间值为最好适当的中间值为最好;区间型指标区间型指标:总是期望指标的取值最好是落在某一个确定的区间内为最总是期望指标的取值最好是落在某一个确定的区间内为最好。好。1.1.极小型指标极小型指标:对于某个极小型指标对于某个极小型指标，则通过变换，则通过变换,或变换或变换,其中其中为指标为指标的可能取值的可能取值的最大值，即可将指标的最大值，即可将指标极大化。极大化。3 评价指标的规范化处理 2 2.中间型指标中间型指标:对于某个中间型指标，则通过变换对于某个中间型指标，则通过变换其中其中和和分别为指标
19、分别为指标的可能取值的最大值和最的可能取值的最大值和最小值，即可将中间型小值，即可将中间型指标极大化。指标极大化。3 3.区间型指标区间型指标：对于某个区间型指标对于某个区间型指标，则通过变换，则通过变换其中其中为指标为指标的最佳稳定的区间，的最佳稳定的区间，和和分别为指标分别为指标的可能取值的最的可能取值的最大值和最小值。即可将区间型指标大值和最小值。即可将区间型指标极大化。极大化。9.3.2 评价指标的无量纲化 在实际中的评价指标在实际中的评价指标之间，往往之间，往往都存在着各自不同的单位和数量级，使得这些指标之间存在着不可公度都存在着各自不同的单位和数量级，使得这些指标之间存在着不可公度
20、性，这就为综合评价带来了困难，尤其是为综合评价指标建立和依据这性，这就为综合评价带来了困难，尤其是为综合评价指标建立和依据这个指标的大小排序产生不合理性。个指标的大小排序产生不合理性。如果不对这些指标作相应的无量纲处理，则在综合评价过程中就会如果不对这些指标作相应的无量纲处理，则在综合评价过程中就会出出“大数吃小数大数吃小数”的错误结果，从而导致最后得到错误的评价结论。的错误结果，从而导致最后得到错误的评价结论。无量纲化处理又称为指标数据的标准化无量纲化处理又称为指标数据的标准化,或规范化处理或规范化处理常用方法常用方法:标准差方法、极值差方法和功效系数方法等。标准差方法、极值差方法和功效系数
21、方法等。假设假设个评价指标个评价指标，在此不妨，在此不妨假设已进行了类型的一致化处理，并都有假设已进行了类型的一致化处理，并都有组样本观测值组样本观测值，则将其作无量纲化处理。，则将其作无量纲化处理。其中其中。显然显然指标的均值和均指标的均值和均方差分别为方差分别为0和和1，即即是无量纲的指标，称之为是无量纲的指标，称之为的标准观测值。的标准观测值。2.极值差方法极值差方法:令令，其中其中。则。则是无量纲的指标观测值。是无量纲的指标观测值。3.功效系数法功效系数法:令令，其中其中均为确定的常数。均为确定的常数。表示表示“平移量平移量”，表示表示“旋转量旋转量”，即表示即表示“放大放大”或或“缩
22、小缩小”倍数，则倍数，则。譬如若取譬如若取，则，则。1.标准差方法标准差方法:令令，4 水质综合评价模型9.4.1基本假设基本假设（1）本文只以长江流域中的）本文只以长江流域中的17个观测点为研究对象，不考虑长江个观测点为研究对象，不考虑长江流域的其它部分和未提到的其它支流。流域的其它部分和未提到的其它支流。（2）假设高锰酸盐和氨氮的降解系数都为）假设高锰酸盐和氨氮的降解系数都为0.2。（3）各年的水质情况的检测数据互不影响。）各年的水质情况的检测数据互不影响。（4）各个污染指标不相关，互不影响。）各个污染指标不相关，互不影响。（5）评价和预测水质时忽略其他环境因素。）评价和预测水质时忽略其他
23、环境因素。（6）各监测站的监测数据代表该地区的水质情况）各监测站的监测数据代表该地区的水质情况（7）长江干流的水流速度均匀变化。）长江干流的水流速度均匀变化。9.4.2符号说明符号说明（1）DO表示水中溶解氧表示水中溶解氧（2）CODMn表示水中高锰酸盐指数表示水中高锰酸盐指数（3）NH3-N表示水中氨氮表示水中氨氮9.4.3问题分析问题分析整个水质评价工作应包括五个环节：整个水质评价工作应包括五个环节：1、确定调查范围，设计检测站点；、确定调查范围，设计检测站点；2、水污染调查与监测，得到各站点的监测值；、水污染调查与监测，得到各站点的监测值；3、确定评价指标与水质分级以及各指标在各等级的标
24、准值；、确定评价指标与水质分级以及各指标在各等级的标准值；4、建立数学模型，进行综合评价；、建立数学模型，进行综合评价；5、划分水环境质量等级，并作出评价结论。、划分水环境质量等级，并作出评价结论。可用下面的流程图（图可用下面的流程图（图1）来表示，其中环节）来表示，其中环节4和环节和环节5是本文要是本文要做的工作，即建立数学模型进行综合评价，并划分水环境质量等级，做的工作，即建立数学模型进行综合评价，并划分水环境质量等级，作出评价结论。作出评价结论。图图1 1 水质评价流程图水质评价流程图9.4.4各观测站点分布图各观测站点分布图水质评价工作的第一步就是确定调查范围，设计检测站点。由题水质评
25、价工作的第一步就是确定调查范围，设计检测站点。由题目中的附件目中的附件3可知此次水质评价工作的范围为长江中下游地区（从四可知此次水质评价工作的范围为长江中下游地区（从四川攀枝花到江苏扬州），在这之间设计了川攀枝花到江苏扬州），在这之间设计了17个监测站点。个监测站点。参照中国地图册，画出长江中下游流域与各个观测站点（地区）参照中国地图册，画出长江中下游流域与各个观测站点（地区）的大致分布，如下面图的大致分布，如下面图3所示：所示：9.4.5监测数据的采取监测数据的采取水质评价工作的第二步是每隔一段时间对水质调查与水质评价工作的第二步是每隔一段时间对水质调查与监测，得到各个指标的监测值序列。题目
26、中的附件监测，得到各个指标的监测值序列。题目中的附件3已给已给出了这出了这17个检测站近两年多主要水质指标的监测数据。个检测站近两年多主要水质指标的监测数据。9.4.6标准的选取标准的选取水质评价工作的第三步是确定评价指标、水质分级水质评价工作的第三步是确定评价指标、水质分级以及质量标准值。这里质量标准值采用以及质量标准值。这里质量标准值采用地表水环境质量地表水环境质量标准标准(GB3838-2002)中的标准值。评价指标为水质的中的标准值。评价指标为水质的3个个主要指标：主要指标：DO（溶解氧）、（溶解氧）、CODMn（高锰酸盐指数）（高锰酸盐指数）和和NH3-N（氨氮）。评价等级设为（氨氮
27、）。评价等级设为6个等级：个等级：类、类、类、类、类、类、类、类、类和劣类和劣类。类。9.4.7模型建立模型建立水质评价工作的第三步是建立数学模型，进行综合评水质评价工作的第三步是建立数学模型，进行综合评价。由于水质污染物浓度受水文流量及污染物排放因素的价。由于水质污染物浓度受水文流量及污染物排放因素的影响较大，存在随机性，而水质综合评价又存在模糊性。影响较大，存在随机性，而水质综合评价又存在模糊性。因此，本文提出了水质评价的模糊概率综合评价模型，把因此，本文提出了水质评价的模糊概率综合评价模型，把概率统计与模糊数学有机地结合起来，它能较全面地评价概率统计与模糊数学有机地结合起来，它能较全面地
28、评价水质状况。水质状况。（1）确定评价指标、水质分级与各指标标准值确定评价指标、水质分级与各指标标准值设评价指标有设评价指标有m个，水域水质分个，水域水质分n级，则级，则评价指标集合评价指标集合U=（1）水质分级集合水质分级集合V=（2）i指标（指标（i=1，2，m），），j级水质（级水质（j=1，2，n）的指标标准值的指标标准值参考参考地表水环境质量标准地表水环境质量标准(GB3838-2002)中的相关数值。中的相关数值。（2）污染物监测值统计概率分析污染物监测值统计概率分析设设i指标污染物监测值共有指标污染物监测值共有个，其中介于个，其中介于至至之之间的监测值有间的监测值有个，那么，对于
29、个，那么，对于i指标而言，介于至之间的监测值指标而言，介于至之间的监测值发生的统计概率为发生的统计概率为（3）其中其中i=1，2，m；j=1，2，n。（3）隶属度刻划水质分级界限隶属度刻划水质分级界限已知水质等级标准值为已知水质等级标准值为，i指标污染物介于指标污染物介于至至之间的之间的个监测值的平均值为个监测值的平均值为，则，则（4）式中，式中，为第为第i指标污染物监测值中介于指标污染物监测值中介于至至之之间的间的个监测值的第个监测值的第k个。个。则对则对i指标而言，分别属于第指标而言，分别属于第j级水质和第级水质和第j-1级水质的程度，即级水质的程度，即对第对第j级水质和第级水质和第j-1
30、级水质的隶属度可由下列各式推求，即级水质的隶属度可由下列各式推求，即对对DO，有，有（5）对其他指标（对其他指标（CODMn，NH3-N），正好与上述方法相），正好与上述方法相反，即反，即（6）（4）相对距离法计算权重相对距离法计算权重根据各评价参数超标情况进行加权，超标越多，加权越大，权重根据各评价参数超标情况进行加权，超标越多，加权越大，权重可采用下式计算，即可采用下式计算，即对对DO：（7）对其他指标（对其他指标（CODMn，NH3-N）：）：（8）式中式中第第i种污染物的监测值算术平均值种污染物的监测值算术平均值DO的饱和浓度的饱和浓度第第i种污染物各级标准值算术平均值种污染物各级标准
31、值算术平均值为了避免这一现象的发生，可用相对距离法计算单项为了避免这一现象的发生，可用相对距离法计算单项指标的权重。具体计算如下：指标的权重。具体计算如下：对对DO：（9）其中其中i=1,2.m,j=1,2.n.对其他指标（对其他指标（CODMn，NH3-N）：）：(10)式中其中式中其中i=1,2.m,j=1,2.n第第i种污染物第种污染物第j级标准值级标准值用（用（9）、（）、（10）式计算权重，方法简单客观，同时也肯定了污染）式计算权重，方法简单客观，同时也肯定了污染超标大者其权重亦大的基本思想，但也存在如下的问题：超标大者其权重亦大的基本思想，但也存在如下的问题：按照国际（按照国际（G
32、B38382002）给出的）给出的地表水环境质量标准地表水环境质量标准，同，同一污染物的五类标准一污染物的五类标准是单向递增的，因此，如果两种指标分别是单向递增的，因此，如果两种指标分别属于属于k类和类和r类水质，设类水质，设kr，按照污染超标大者其权重亦大，则应有，按照污染超标大者其权重亦大，则应有。但由于不同的污染指标的相邻两级标准值增长不同步，可能造成上。但由于不同的污染指标的相邻两级标准值增长不同步，可能造成上述次序的混乱。述次序的混乱。为了进行模糊运算，各单项指标权重进行归一化处理，即为了进行模糊运算，各单项指标权重进行归一化处理，即(i=1,2,.m)(11)m个指标权重构成权重矩
33、阵，即个指标权重构成权重矩阵，即(12)（5）单指标模糊概率评价单指标模糊概率评价根据上述方法，分别求出根据上述方法，分别求出i指标，并分别统计指标，并分别统计j级水和级水和j-1级水的隶级水的隶属度和，以及介于至之间的属度和，以及介于至之间的i指标监测值发生的统计概率，就可以应用指标监测值发生的统计概率，就可以应用模糊概率方法对模糊概率方法对i指标进行单指标评价。指标进行单指标评价。i指标发生各级水质的模糊概率分别为指标发生各级水质的模糊概率分别为(13)其中其中i=1，2，m。i指标分别发生指标分别发生n级水质的概率就构成级水质的概率就构成i指标单指标模指标单指标模糊概率评价矩阵，即糊概率
34、评价矩阵，即m个指标的单指标模糊概率矩阵就构成模糊概率关系个指标的单指标模糊概率矩阵就构成模糊概率关系矩阵，即矩阵，即(14)（6）多指标模糊概率综合评价多指标模糊概率综合评价已求出权重矩阵及模糊概率关系矩阵，则对已求出权重矩阵及模糊概率关系矩阵，则对m个评价个评价指标综合评价，得到多指标模糊概率综合评价矩阵，即指标综合评价，得到多指标模糊概率综合评价矩阵，即(15)其中，其中，(16)即为水域水体出现即为水域水体出现j级水质的模糊综合概率。级水质的模糊综合概率。（7）综合评价指数综合评价指数定义水域水质综合评价指数定义水域水质综合评价指数P为各个等级水质出现的为各个等级水质出现的概率与等级序
35、号相乘在求和（类似于概率论中的全概率公概率与等级序号相乘在求和（类似于概率论中的全概率公式），即为：式），即为：(17)9.4.8模型求解模型求解把水质评价工作中的第二步工作所采集的监测数据代入模型，把水质评价工作中的第二步工作所采集的监测数据代入模型，经过计算可求得结果如下表：经过计算可求得结果如下表：等级等级等级等级等级等级综合评价指数四川攀枝花0.84040.09860.05550.0055001.2261重庆朱沱0.74040.24590.01370001.2737湖北宜昌南津关0.70870.26630.02500001.3163湖南岳阳城陵矶0.48400.46620.049900
36、01.5659江西九江河西水厂0.81360.186400001.1864安徽安庆皖河口0.72200.27150.00650001.2845江苏南京林山0.76800.23060.00140001.2335四川乐山岷江大桥0.15920.37740.26580.14200.055602.5573四川宜宾凉姜沟0.70390.20210.041860.02870.023401.4657四川泸州沱江二桥0.50270.27590.12880.04160.05100.04922.1575湖北丹江口胡家岭0.98760.012400001.0124湖南长沙新港0.30800.29950.29080.
37、1018002.1864湖南岳阳岳阳楼0.42670.42030.14410.0090001.7355湖北武汉宗关0.58570.36270.04850.0030001.4688江西南昌滁槎0.17580.12990.26780.11640.31020.30515.0859江西九江蛤蟆石0.57530.31520.08360.0259001.5601江苏扬州三江营0.58270.36640.04860.0023001.4706结果分析：由上表可以发现，长江流域干流上监测站点水质都在结果分析：由上表可以发现，长江流域干流上监测站点水质都在第一等级和第二等级之间，水质情况较好。长江流域支流上部分监
38、测第一等级和第二等级之间，水质情况较好。长江流域支流上部分监测点水质处于第一等级和第二等级之间，水质情况较好；部分监测点点水质处于第一等级和第二等级之间，水质情况较好；部分监测点（四川乐山岷江大桥、四川泸州沱江二桥（四川乐山岷江大桥、四川泸州沱江二桥和湖南长沙新港）和湖南长沙新港）处于第处于第二与第三等级之间，水质受到轻度污染；但在监测站江西南昌滁槎二与第三等级之间，水质受到轻度污染；但在监测站江西南昌滁槎（赣江）水质处于劣五级，污染十分严重。通过查阅资料发现（赣江）水质处于劣五级，污染十分严重。通过查阅资料发现：由于：由于目前赣江的主要流经城市南昌市内大量工业废水和生活废水不经处理目前赣江的
39、主要流经城市南昌市内大量工业废水和生活废水不经处理就排入江内，造成氨氮、粪大肠菌群浓度过高，而南昌市只建有一个就排入江内，造成氨氮、粪大肠菌群浓度过高，而南昌市只建有一个日处理污水能力日处理污水能力8万吨的朝阳污水处理厂，但对于目前市内每日万吨的朝阳污水处理厂，但对于目前市内每日100万万吨的污水排放量，城市处理污水能力远远跟不上工业的的发展，造成吨的污水排放量，城市处理污水能力远远跟不上工业的的发展，造成现在赣江的水体质量严重超标，达到重度污染级，这里是改善长江水现在赣江的水体质量严重超标，达到重度污染级，这里是改善长江水质的重点处理区之一。质的重点处理区之一。总的说来，长江流域（干流及支流
40、）的大部分水质情况还是比较总的说来，长江流域（干流及支流）的大部分水质情况还是比较理想的，但是个别地区污染比较严重，不容忽视。理想的，但是个别地区污染比较严重，不容忽视。河流的水质污染程度随径污比变化。河流的径流量和排入河流的河流的水质污染程度随径污比变化。河流的径流量和排入河流的污水、污物量决定了径污比。排污量相同，河流的径流量大则污染轻，污水、污物量决定了径污比。排污量相同，河流的径流量大则污染轻，反之就重。河流的径流量随时间变化，因此河水污染程度也随时间变反之就重。河流的径流量随时间变化，因此河水污染程度也随时间变化。这点在模型求解过程中没有考虑到，将在下面的模型的进一步求化。这点在模型
41、求解过程中没有考虑到，将在下面的模型的进一步求解中给出。解中给出。9.4.9模型的进一步求解：模型的进一步求解：（1）时间的分段处理）时间的分段处理现在先对时间进行分段处理，即分为丰水期（现在先对时间进行分段处理，即分为丰水期（5月到月到10月）、枯水期（月）、枯水期（1月到月到4月）和平水期（月）和平水期（11月到月到12月）月）三个时期。然后对各个时期的监测数据进行模糊概率综合三个时期。然后对各个时期的监测数据进行模糊概率综合评价，计算得到的结果如下表：评价，计算得到的结果如下表：表表2各监测点各时期污染等级分布各监测点各时期污染等级分布监测点丰水期平水期枯水期水文年四川攀枝花1.3218
42、1.05591.19711.2261重庆朱沱1.17391.35721.51021.2737湖北宜昌南津关1.44821.12481.19011.3163湖南岳阳城陵矶1.56851.51781.60291.5659江西九江河西水厂1.24221.1061.14381.1864安徽安庆皖河口1.37671.15241.19721.2845江苏南京林山1.38461.1681.1241.2335四川乐山岷江大桥2.21792.37153.40492.5573四川宜宾凉姜沟1.28831.78121.96451.4657四川泸州沱江二桥1.66731.82943.80492.1575湖北丹江口胡家
43、岭1.01641.02171.00211.0124湖南长沙新港2.00022.66422.432.1864湖南岳阳岳阳楼1.66131.70741.92981.7355湖北武汉宗关1.57231.5191.31451.4688江西南昌滁槎4.28587.4335.63315.0859江西九江蛤蟆石1.35382.43181.66691.5601江苏扬州三江营1.65881.28781.26291.4706结果分析：由表结果分析：由表2可以发现，三个时期（丰水期、平水期和枯水期）可以发现，三个时期（丰水期、平水期和枯水期）各监测站的污染等级分布与前面模型求解得到的全年的综合评价等级各监测站的污染
44、等级分布与前面模型求解得到的全年的综合评价等级基本一致，三个时期之间的等级分布也基本相同，说明该河流的水质基本一致，三个时期之间的等级分布也基本相同，说明该河流的水质情况不会随季节发生非常显著性的变化。当是还是存在一定的变化，情况不会随季节发生非常显著性的变化。当是还是存在一定的变化，17个监测站大致的可以分为丰水期污染程度小于枯水期的污染程度和个监测站大致的可以分为丰水期污染程度小于枯水期的污染程度和丰水期污染程度大于枯水期的污染程度两个类别：丰水期污染程度大于枯水期的污染程度两个类别：1、丰水期污染、丰水期污染枯水期污染枯水期污染原因：丰水期雨量充足，水流量大，容易从上游或河流周边地区原因
45、：丰水期雨量充足，水流量大，容易从上游或河流周边地区带来大量污染物而污染下游水质，如果这个负面作用过强就会超过因带来大量污染物而污染下游水质，如果这个负面作用过强就会超过因径流量大而使径污比减少的正面作用，因此丰水期污染反而更加严重。径流量大而使径污比减少的正面作用，因此丰水期污染反而更加严重。特别是如果雨量过大在河流流经地区发生洪涝灾害，则污染将更加严特别是如果雨量过大在河流流经地区发生洪涝灾害，则污染将更加严重。重。（2）污染范围与污染历时评价）污染范围与污染历时评价全面的水质状况，本来应从三个方面来评定：第一是污染强度，全面的水质状况，本来应从三个方面来评定：第一是污染强度，即各种污染物
46、在水体中的浓度；第二是污染范围，即在水域中各种污即各种污染物在水体中的浓度；第二是污染范围，即在水域中各种污染物所影响的空间大小；第三是污染历时，即各种污染物在水域中污染物所影响的空间大小；第三是污染历时，即各种污染物在水域中污染所持续的时间。染所持续的时间。上面只对污染强度一个方面进行了评价，其实模糊概率综合模型上面只对污染强度一个方面进行了评价，其实模糊概率综合模型也可以容易地得到另外两个方面的评价。下面举例来说明用模糊概率也可以容易地得到另外两个方面的评价。下面举例来说明用模糊概率综合模型是如何来评价污染物的污染范围的（对污染历时的评价也类综合模型是如何来评价污染物的污染范围的（对污染历
47、时的评价也类似可得）。似可得）。以污染物以污染物CODMn为例来对其污染范围进行评价：为例来对其污染范围进行评价：用单指标模糊概率评价法对长江流域用单指标模糊概率评价法对长江流域17个观测站的个观测站的CODMn污染污染状况进行评价，结果如下表：状况进行评价，结果如下表：表表3各监测点各监测点CODMn污染状况污染状况 等级等级等级等级等级等级综合等级四川攀枝花0.68040.21790.10090001.4223重庆朱沱0.82500.16790.00710001.1821湖北宜昌南津关0.59460.35180.05360001.4589湖南岳阳城陵矶0.22320.66070.11610
48、001.8929江西九江河西水厂0.77320.226800001.2268安徽安庆皖河口0.70180.28040.01790001.3161江苏南京林山0.85000.150000001.15四川乐山岷江大桥0.15360.41790.25540.1732002.4482四川宜宾凉姜沟0.64110.31250.02860.0179001.4232四川泸州沱江二桥0.43390.40710.15360.0054001.7304湖北丹江口胡家岭0.9750.02500001.025湖南长沙新港0.717960.24110.04110001.3232湖南岳阳岳阳楼0.21790.49820.2
49、6340.0205002.0866湖北武汉宗关0.39290.55180.05540001.6625江西南昌滁槎0.76430.150.04910.0366001.358江西九江蛤蟆石0.46430.30710.16960.0589001.8232江苏扬州三江营0.51790.40360.07860001.5607结果分析：由上表可以发现，长江中下游流域的结果分析：由上表可以发现，长江中下游流域的CODMn污染主污染主要分布在四川岷江上游（与大渡河汇合前）、湖南洞庭湖和江西鄱阳要分布在四川岷江上游（与大渡河汇合前）、湖南洞庭湖和江西鄱阳湖地区。由于岷江上游距长江干流较远，所以它对长江干流的湖地
50、区。由于岷江上游距长江干流较远，所以它对长江干流的CODMn污染影响较小，而洞庭湖和鄱阳湖都与长江相接，它们的污污染影响较小，而洞庭湖和鄱阳湖都与长江相接，它们的污染直接影响到长江干流的污染，例如湖南岳阳城陵矶的染直接影响到长江干流的污染，例如湖南岳阳城陵矶的CODMn污染污染就比较严重。就比较严重。因此，总的来讲，长江干流的因此，总的来讲，长江干流的CODMn污染范围不是很大，主要污染范围不是很大，主要是分布在长江与洞庭湖、鄱阳湖相接的河段及其附近河段。是分布在长江与洞庭湖、鄱阳湖相接的河段及其附近河段。9.4.10模型的优点模型的优点（1）把水污染监测浓度看成一个离散的随机变量，用概率统计