长江水质评价与预测的数学模型

长江水质评价与预测的数学模型

摘 要

首先,利用附录3的数据，从时间上和空间上分析长江流域水质，得出长江水质的污染程度有所增加，但不明显.大部分污染比较严重的地区都位于支流上. 在求解主要污染物高锰酸盐指数和氨氮的污染源时先计算出单独一个河段内的排污量，进而求出一个河段内包括降解的污染物总量. 计算出长江干流近一年多高锰酸盐指数和氨氮的主要污染源是位于重庆朱沱至湖北宜昌之间（每月排放约230万吨高锰酸盐和21万吨氨氮）、湖北宜昌和湖南岳阳之间（每月排放约206万吨高锰酸盐和20万吨氨氮）的工业带. 然后,借鉴马氏链模型，结合过去十年水文年全流域数据，拟合出转移矩阵，预测得到未来十年各类水占河长的百分比(如2009年长江劣V类水占河段长约26.28%，2014年长江劣V类水河段占河段长约36.93%).如果未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,2005年～2014年每年需要处理的废水分别为(单位：亿吨）：25.48，24.89，26.59，29.42，32.82，36.48，40.23，44.01，47.76，51.48.

最后,针对目前长江水质污染状况，提出了切实可行的建议，具有较强的参考价值.

关键词:马氏链模型;水质评价;水质污染;数据拟合

1

1 问题的提出

水是人类赖以生存的资源，保护水资源就是保护我们自己，对于我国大江大河水资源的保护和治理应是重中之重.专家们呼吁：“以人为本，建设文明和谐社会，改善人与自然的环境，减少污染.”

长江是我国第一、世界第三大河流，长江水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视.2004年10月，由全国政协与中国发展研究院联合组成“保护长江万里行”考察团，从长江上游宜宾到下游上海，对沿线21个重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心.为此，专家们提出“若不及时拯救，长江生态10年内将濒临崩溃”，并发出了“拿什么拯救癌变长江”的呼唤.

研究下列问题：

(1)对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况.

(2)研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在哪些地区?

(3)假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，比如研究未来10年的情况.

(4)根据预测分析，如果未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水,每年需要处理多少污水？

(5)对解决长江水质污染问题切实可行的建议和意见.

2 模型的假设

2.1 长江干流的自然净化能力均匀；

2.2水质只受高锰酸盐指数(CODMn)、氨氮(NH3-N)、溶解氧(DO) 、PH值的影响，与其

它因素无关； 2.3 污染源均匀分布在河岸两侧； 2.4 不同水质的水随机分布在全流域上； 2.5 两个监测站之间水流速均匀变化.

3 符号的约定

Wk1 表示第k监测站13个月的CODMn总量(包含上一个监测站对它的影响)

Wk2 表示第k监测站13个月NH3-N的总量(包含上一个监测站对它的影响) Wk1' 表示第k监测站13个月CODMn的总量(不包含上一个监测站对它的影响)

2

Wk2' 表示第k监测站13个月NH3-N的总量(不包含上一个监测站对它的影响) C0 表示综合降解系数

?ijk 表示在第k个监测站第i个月第j种污染物的浓度

Qo(i) 第i年长江废水排放总量

leni 第i个河段的长度

vij 第i月第j个监测点的平均水流速 Mij 第i月第j个监测点的平均水流量

4 模型的建立与求解

4.1 长江近两年多水质情况与各地区水质污染状况的定量综合评价 4.1.1 污染物含量的变化规律

要看出近两年多长江水质的变化情况，我们首先要找出主要污染物的含量随时间的变化规律.由题目我们可以知道，影响水质的主要污染物有高锰酸盐指数(CODMn)、氨氮(NH3-N)、溶解氧(DO) 、PH值.把长江全流域看成一个整体，W(i,j)为全流域第i个月第j种污染物的总量，比，即：

Pij?W(i,j) (i?1,2,??,28;j?1,2,3,4)，其中Wij??Vi*?ijk,(k?1,2,?,17),

k?117Pij为第i个月第j种污染物总量占28个月总污染物总量的百分

?W(i,j)i?128Vi表示第i个月的水流量，?ijk表示在第k个监测站的第i个月第j种污染物的浓度.

这里的

Pij都是百分率，数值上有良好的可比性，能很好的反映出三种污染物含量随

着时间(月份)的变化规律.作出28个月长江流域四项监测项目总量所占百分比随时间的变化曲线图，如图1所示：

3

图1

由图1我们可以很直观地看出污染物含量的变化规律：

（1） 污染物总量是以年为单位成周期性变化的，且相对总量逐年增加； （2） 污染物总量与水流量的变化趋势基本一致；

由附录4我们可以知道，每年长江的枯水期为1月～4月，丰水期为5月～10月平水期为11月～12月，在丰水期，污染物的总量有明显的增长，随着枯水期的到来，污染物的总量也随之减少.

（3） 四种污染物的变化趋势基本一致；

（4） 主要污染物高锰酸盐和溶解氧含量都逐年增加，氨氮的含量有所下降. 4.1.2 长江整体水质情况及各地区水质的污染状况

首先，把17个观测站按地理位置由西到东、由支流到干流的顺序重新排列(1、四川攀枝花龙洞 2、四川乐山岷江大桥 3、四川宜宾凉姜沟 4、四川泸州沱江二桥 5、重庆朱沱 6、湖北宜昌南津关 7、湖南长沙新港 8、湖南岳阳岳阳楼 9、湖南岳阳城陵矶 10、湖北丹江口胡家岭 11、湖北武汉宗关 12、江西南昌滁槎 13、江西九江蛤蟆石 14、江西九江河西水厂 15、安徽安庆皖河口 16、江苏南京林山 17、江苏扬州三江营)，根据附录3算出各观测站28个月四种污染物的平均浓度：

?(j,k)???(i,j,k)i?12828，(j?1,2,3,4;k?1,2,?,17;i?1,2,?,28)

4

其中?(i,j,k)表示第i个月第k个监测站第j种污染物的浓度.

作出四种主要监测项目(pH*、DO、CODMn、NH3-N)的平均浓度与各观测站之间的关系图，如图2所示：

图2

由图2我们可以很直观地看出各观测站4种污染物的浓度曲线，得污染物浓度的分布情况：

（1） 支流比干流污染严重，大部分重污染地区都位于支流上； （2） 四川乐山岷江大桥的氨氮和高锰酸盐的含量都很高； （3） 江西南昌滁槎的氨氮含量很高；

（4） 湖北宜昌南津关、湖南岳阳岳阳楼、湖北武汉宗关的高锰酸盐含量都很高. 由于重庆观测站朱沱位于川、渝省界，尚未进入重庆，未能反映重庆河段的污染情况.

综上所述，长江水质的污染程度虽有所增加，但不明显，这是由于长江水量目前还比较大，掩盖了问题的严重性.

4.2研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的主要污染源 4.2.1 流过各监测点的高锰酸盐和氨氮总量

由附录3的数据可以知道，2004.04～2005.04长江干流主要观测站点的平均水流量及其所对应的高锰酸盐指数和氨氮浓度，可以求出在这13个月内，每个月流过干流各

5

监测点的高锰酸盐和氨氮的总量：

13Wjk??ijk*Mik，(k?1,2,?,7)所以13个月流过干流

Qkj??Wijki?1监测点污染物为： (j?1表示CODMn含量，j?2表示NH3?N含量).得到流

经各监测点的污染物量如下表（单位：吨/月）：

观测站点 攀枝花 重庆朱沱 湖北宜昌 湖南岳阳 江西九江 安徽安庆 江苏南京 高锰酸盐 0.3078 (×10^6) 氨氮 0.1650 (×10^5) 1.2735 1.7440 2.4135 1.9937 2.0524 1.9958 1.0447 1.4189 2.1878 1.7291 1.7223 0.8494 流经4、5、6三个监测点的高锰酸盐和氨氮的总量最多，但这里没有消除上游河段的影响，

Wjk是包含上游全河段污染物的总量，不利于我们分析主要污染物高锰酸盐和

氨氮的污染源主要在哪些地区，因而求出各个河段单独的污染物排放量.以下我们考虑消除前面河段的影响.

4.2.2 单独一个河段内的排污量

dc??c0?积分解得???0e由河流污染物一维稳定衰减规律的微分方程Vdx?c0xV,其中

X表示测试河段离污染源的距离，V表示水流在该河段的平均速度，?0表示前一节点污染物的浓度，?表示前一节点对测试点的影响浓度.

2004.04～2005.04长江干流7个观测站点平均水流量、不同时间的高锰酸盐指数和氨氮含量、不同时间的水流数度，可以求出在每一个月份内，在监测点上一河段内所排放的高锰酸盐和氨氮的总量： Wi'?Wi?Wi?1e?CoxV ，把十三个月的总量相加：

Wkj??(Wikj?Wikjei?113?CoxV),k?1,2,?,7j?1,2;（程序代码见附录1）得到在长江干流不同

河段近一年多主要污染物高锰酸盐和氨氮的排放总量（单位：吨/月）：

观测点 攀枝花 重庆朱沱 湖北宜昌 湖南岳阳 江西九江 安徽安庆 江苏南京 高锰酸盐 0.3078 1.1137 1.2083 1.5759 1.0228 0.5746 1.0758 ×10^6 氨氮 0.1650 0.9669 1.0933 1.5251 0.8767 0.4081 0.0543 ×10^5

6

从表中可知重庆朱沱到湖南岳阳城陵矶之间河段排放的高锰酸盐和氨氮总量最多，其中重庆朱沱到湖北宜昌南津关之间的河段每个月排放的高锰酸盐和和氨氮分别有1.2083×10^6吨、1.0933×10^5吨，湖北宜昌南津关到湖南岳阳城陵矶之间和的河段每个月排放的高锰酸盐和氨氮分别有1.5759×10^6吨、1.5251×10^5吨.干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要在重庆朱沱到湖南岳阳城陵矶之间河段. 4.2.3 每个河段的总排放量(包括降解的污染物)

长江两岸工业分布呈工业带形态，这使对应的污染物排放量也呈带状，比如九江市岸边污染带约8公里长,重庆市在80年代就形成了81公里的岸边污染带，武汉几个大的排污口形成了几公里长的污染带.

所以我们引进“单位长度江段污染物排放量”，即每一河段一公里排放的污染物量a.设第i与第i?1监测点之间的水流速是均匀变化的，得到单位长度江段水流速变化为

??(vi?1?vi)leni，第i个河段内的有leni?1个单位的排放量，第j个单位河段排放的污染

?C0jvi?j?物经过降解流到下一个监测点“剩余”的污染物为ae余”的总污染物为

leni?1j?1leni?1j?1，leni?1个单位排放量的“剩

?ae?Cojvi?j?.当a表示高锰酸盐排放量时，这就是我们上面统计的Wk1'，

有Wk1'??aek?Cojvi?j?，解出各段对应得a值，同理可以算出氨氮排放量b值.如下所示(单

位：吨/月)：

a值： 1817.8 2994.6 b值： 157.80 270.96

5252.1

508.25

3558.3 4229.0 3545.9 305.02 310.03 17.89

攀枝花 重庆朱沱 湖北宜昌 湖南岳阳 江西九江 安徽安庆 江苏南京 对应的河段排放总量(单位：10^5吨/月):

攀枝花 重庆朱沱 湖北宜昌 湖南岳阳 江西九江 安徽安庆 江苏南京 17.251 23.268 20.693 17.756 6.893 16.418 1.4975 2.1054 2.0025 1.5221 0.5054 0.0828 重庆朱沱到湖南岳阳城陵矶之间河段排放的高锰酸盐和氨氮总量都是最多的，其中

重庆朱沱到湖北宜昌南津关之间的河段每个月排放的高锰酸盐和和氨氮分别有2.3268×10^6吨、2.1054×10^5吨，湖北宜昌南津关到湖南岳阳城陵矶之间河段每个月排放的高锰酸盐和氨氮分别有2.0693×10^6吨、2.0025×10^5吨.

7

综合5.2.2得出的数据，我们可以确定干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主要集中在在重庆朱沱到湖南岳阳城陵矶之间的河段. 4.3长江未来10年水质污染的发展趋势 4.3.1 模型的准备

马氏链模型及其基本方程 按照系统的发展，时间离散化为n?0,1,2,?，对每个n，系统的状态用随机变量Xn表示，设Xn可以取k个离散值Xn?1,2,?,k，且记

ai(n)?P(Xn?i)，即状态概率，从Xn?i到Xn?1?j的概率记pij?P(Xn?1?j|Xn?i)，即转移概率.如果Xn?1的取值只取决于Xn的取值及转移概率，而与Xn?1,Xn?2,?的取值无关，那么这种离散状态按照离散时间的随机转移过程称为马氏链.由状态转移的无后效性和全概率公式可以写出马氏链的基本方程为

ai(n?1)??aj(n)pjij?1k i?1,2,?,k

并且ai(n)和

pij应满足：

?a?n??1ii?1k， n?0,1,2,??

pij?0， i,j?1,2,?,k

?pj?1kij?1， i?1,2,?,k

4.3.2 水质污染的发展趋势

根据水环境质量的不同，把水环境分为6类(Ⅰ类、Ⅱ类、Ⅲ类、Ⅳ类、Ⅴ类、劣Ⅴ类)，其中Ⅰ、Ⅱ、Ⅲ类为可饮用水.各类水质间可以相互转化，且下一个状态只取决于上一状态，与以前的状态无关，这正符合马氏链模型.

把6类水质分为四种状态：Ⅰ、Ⅱ类同属于状态1；Ⅲ类属于状态2；Ⅳ、Ⅴ类同属于状态3；劣Ⅴ类属于状态4.

用随机变量Xn表示第n年的状态，Xn?1表示水质属于Ⅰ类或Ⅱ类，Xn?2表示水质属于Ⅲ类，Xn?3表示水质属于Ⅳ类或Ⅴ类，Xn?4表示水质属于劣Ⅴ类，n=0,1,2,?.用ai(n)表示第n年处于状态i的概率，i=1,2,3，4，即ai(n)?P(Xn?i).用

pij表示已知今年处于状态i，来年处于状态j的概率，i,j=1,2,3，4，即

8

pij?P(Xn?1?j|Xn?i).第n+1年的状态Xn?1只取决于第n年的状态，Xn和转移概率

pij，而与以前的状态Xn?1,Xn?2,?无关，即状态转移具有无后效性.第n+1年的状态概

率可由概率公式得到:

?a1(n?1)?a1(n)p11?a2(n)p21?a3(n)p31?a4(n)p41?a(n?1)?a(n)p?a(n)p?a(n)p?a(n)p?2112222332442??a3(n?1)?a1(n)p13?a2(n)p23?a3(n)p33?a4(n)p43??a4(n?1)?a1(n)p14?a2(n)p24?a3(n)p34?a4(n)p44

这里的ai(n)为第n年处于第i种状态的水域长度占河长的百分比.对于不同的转移矩阵P，来年有不同的百分比，结合附录4十年不同水质变化数据，可以建立以下模型：

min

?xi?1921i+

?xi?1922i+

?xi?1923i+

?xi?1924i , 其中

x1i?a1(i?1)?a1(i)p11?a2(i)p21?a3(i)p31?a4(i)p41； x2i?a2(i?1)?a1(i)p12?a2(i)p22?a3(i)p32?a4(i)p42； x3i?a3(i?1)?a1(i)p13?a2(i)p23?a3(i)p33?a4(i)p43； x4i?a4(i?1)?a1(i)p14?a2(i)p24?a3(i)p34?a4(i)p44；

s.t.

?pj?14ij?1, i=1,2,3,4； ??(1) ； ??(2)

0?pij?1由模型的实际意义，不同状态水质之间的转化应满足关系： (1)向自身状态转化的概率不小于向其它状态转化的概率； (2)向相邻状态转化的概率不小于向相隔状态转化的概率；

p11-p12?0;p12-p13?0;p13-p14?0; ??(3)

p22-p21?0;p22-p23?0;p21-p24?0; p23-p24?0; ??(4)

9

p34-p31?0;p32-p31?0;p33-p32?0; p33-p34?0; ??(5)

p42-p41?0;p43-p42?0; p44-p43?0; ??(6)

对于不同的水期，不同的流域可以求出对应的转移矩阵P，然后再递推出长江未来十年各类水质的变化情况.现对水文年的全流域十年数据进行分析，用数学软件lingo8

?0.451?0.445进行求解，得到转移矩阵P???0.215??00.4510.4490.21500.0980.1060.35500??0?,其对应可以求出长江未来十0.215??1?年水文年的各类水质a(i)?a(i?1)*P, i?11...20; 有：

第11年 第12年 第13年 第14年 第15年 第16年 第17年 第18年 第19年 第20年 状态1(％) 状态2(％) 状态3(％) 状态4(％) 34.8610 35.0265 14.3620 15.7505 34.3794 34.5246 12.2577 18.8384 33.4866 33.6298 11.4098 21.4738 32.5040 32.6434 10.9257 23.9269 31.5183 31.6536 10.5521 26.2760 30.5534 30.6847 10.2172 28.5447 29.6156 29.7429 9.9002 30.7414 28.7058 28.8292 9.5951 32.8699 27.8238 27.9433 9.3000 34.9329 26.9688 27.0847 9.0141 36.9324 用枯水期全流域的数据求出干流未来十年水质的变化情况：

状态1(％) 第11年 31.3932 第12年 30.4556 第13年 29.1473 第14年 27.8598 第15年 26.6261 第16年 25.4468 第17年 24.3196 第18年 23.2424 第19年 22.2129 第20年 21.2290

状态2(％) 35.8743 34.4546 32.9880 31.5371 30.1416 28.8067 27.5308 26.3113 25.1459 24.0321 10

状态3(％) 15.7575 14.1754 13.4065 12.7932 12.2241 11.6823 11.1648 10.6703 10.1976 9.7459 状态4(％) 16.9750 20.9144 24.4582 27.8099 31.0082 34.0642 36.9848 39.7760 42.4435 44.9930

用丰水期全流域的数据求出干流未来十年水质的变化情况：

状态1(％) 第11年 29.9735 第12年 29.6774 第13年 29.2404 第14年 28.7597 第15年 28.2531 第16年 27.7306 第17年 27.1997 第18年 26.6656 第19年 26.1322 第20年 25.6023

用水文年干流和支流的数据可以求出各自未来十年水质的变化情况（见附录2、3）. 通过对5个表格的分析对比，我们可以知道长江水质变化呈现统一趋势，劣Ⅴ类水河长百分比逐年增加，支流污染比干流污染严重(如水文年劣Ⅴ类水十年变化：干流10.5％～23.6％，支流15.4％～37.6％)，枯水期劣Ⅴ类水河长的变化比丰水期快(枯水期17％～45％，丰水期12.7％～28.8％).在不采取有效治理措施的情况下，劣Ⅴ类水河长百分比会以一个很快的增长速度不断增大，致使长江陷入深度危机，若不及时拯救，10年之内，长江水系生态将濒临崩溃. 4.4在满足要求情况下的污水处理方法

第i年排放的污染物为 Qo(i)，水文年干流第i年四个状态水量百分比为

a1(i),a2(i),a3(i),a4(i)，有关系Qo(i)?a*a1(i)?b*a2(i)?c*a3(i)?d*a4(i)?e i?1,2,?,10，其中a,b,c,d,e为常系数.对过去十年数据，用最小二乘法拟和得到Qo(i)和a1(i),a2(i),a3(i),a4(i)的关系为：

状态2(％) 37.0704 36.3473 35.6879 35.0172 34.3389 33.6590 32.9815 32.3096 31.6455 30.9907 状态3(％) 20.2299 19.1445 18.2493 17.5023 16.8664 16.3143 15.8255 15.3851 14.9821 14.6082 状态4(％) 12.7262 14.8307 16.8223 18.7208 20.5416 22.2962 23.9934 25.6397 27.2402 28.7988 Q0(i)=1.612*a2(i)+2.706*a3(i)+3.983*a4(i)+104.312;

用5.3计算得到的未来十个水文年干流四个状态水量百分比，而我们的废水处理目标要求要把a4(i)?0而且a3(i)?20，即Q0(i)'?1.1612*a2(i)?2.706*0.2?104.312这样可以得到每年要处理的污水量为：Q0(i)?Q0(i)'，即（单位：亿吨）：

25.48，24.89，26.59，29.42，32.82，36.48，40.23，44.01，47.76，51.48 . 4.5对解决长江水质污染问题的建议和意见

11

意见：目前长江造成局部污染严重而总体水质良好的原因不在水量和水源缺乏，而在水

源浪费，废水处理率低甚至出现“直排”现象，人们环保意识薄弱. 建议：（1）抓好对人民环保意识的教育，进一步提高水资源保护意识，增强对水资源保

护重要性、紧迫性的认识唤醒民众的环保意识. （2）强化流域水资源保护机构，制定流域管理法规，统一规划和监督管理水环境

质量，防治水环境污染. （3）依据水环境功能用途要求和水体稀释自净能力，建立不同类型的保护区，优

先保护好生活用水水源，使其不受污染.

（4）积极引进和开发无废或少废，不用水或少用水的工业技术，研究适合流域工

业特色和自然环境特点的废污水处理利用和资源化技术，加快建设城市污水处理设施和资源化工程，降低产业水耗，提高废水利用率，使有限的水资源发挥更大的经济、社会和环境效益. （5）以水环境质量目标和污染物总量控制目标为导向，合理规划工业布局和调整

工业结构，严格执行环保“三同时”制度和乡镇企业污染防治法规，使新增废污水及污染物排放量得到有效控制. （6）从生态环境和社会发展对水环境的需求出发，确定社会和环境可承受的水资

源开发程度，合理规划水利工程布局、调节调度水资源水量，使水资源的环境功能得以充分的发挥和利用.

5 模型的评价与推广

5.1 模型的评价 5.1.1 模型的优点：

（1） 模型能抓住影响水质的主要因素（高锰酸盐指数(CODMn)、氨氮(NH3-N)、溶

解氧(DO)），能正确预测长江水质变化，过程清晰明了，结果科学合理; （2） 模型具有较好的通用性，实用性强，对现实有很强的指导意义;

（3） 在求解主要污染物高锰酸盐指数和氨氮的污染源时先计算出单独一个河段内的排污量，进而求出一个河段内包括降解的污染物总量. 5.1.2 模型的不足以及需要改进的地方：

（1）监测点太少，不能全面反映长江流域水质; （2）没有考虑长江沿岸地形对流速、流量的影响. 5.2 模型的推广

我们建立模型的方法和思想对其他类似的问题也很适用，适用于其它河流、湖泊水质的评价和预测.

12

参考文献

[1] 姜启源等,数学模型[M],北京：高等教育出版社，2004.

[2] 李强等,Maple基础应用教程[M],北京：中国水利水电出版社,2004. [3] 宋兆基等，MATLAB6.5在科学计算中的应用[M]：清华大学出版社.2005. [4] 总参谋部测绘局编制，中华人民共和国地图集[M]：星球地图出版社.2000. [5] 董哲仁等.中国江河1000问[M]：黄河水利出版社.2001.

附录

附录1

单独一个河段内的排污量

da 是221*6 的矩阵，记录了题目所给附录3的04年4月～05年4月的数据 v 是13*7矩阵，记录了题目所给附录3的水流量数据 len=[ 0 950 1728 2123 2623 2787 3251];

w_Mn=zeros(1,7); % 临时记录7个监测点的高锰酸盐含量 w_NH=zeros(1,7); % 临时记录7个监测点的氨氮含量

pure1=zeros(1,7); % 记录7个监测点13个月的高锰酸盐含量 pure2=zeros(1,7); % 记录7个监测点13个月的氨氮含量

pure1(1)=116820 ; pure2(1) =6261; % 四川攀枝花监测点为第一个河段的上点 flag=0; flag2=0; for i=1:221

temp=mod(i,17);

if temp~=0 && temp<=7 t=temp; flag=flag+1; m=mod(flag,7); if m==0 m=7; end

n=(flag-m)/7;

w_Mn(t)=da(i,3)*V(n,m); w_NH(t)=da(i,4)*V(n,m); end

if mod(i,17)==0 flag2=flag2+1; w1=w_Mn; w2=w_NH;

vv=v(flag2,:);

for j=7:-1:2 ％ 第j个监测点每一个月的污染物量 Co=0.2; ff=1; ％ 减去了上一河段的影响

time=2*(len(j)-len(j-1))/((vv(j)+vv(j-1))*3.6*24); wo=w1(j-1); temp1=wo*exp(-Co*time); wo2=w2(j-1); temp2=wo2*exp(-Co*time); w1(j)=w1(j)-temp1; w2(j)=w2(j)-temp2; pure1(j)=pure1(j)+w1(j);

13

pure2(j)=pure2(j)+w2(j); end end end

plot(pure1,'-p') hold on

plot(pure2,'-*') grid on

附录2

用水文年干流的数据可以求出干流未来十年水质的变化情况：

状态1(％) 状态2(％) 状态3(％) 状态4(％) 第11年 29.5362 40.2088 19.6769 10.5780 第12年 30.8505 39.3196 17.6890 12.1409 第13年 31.2546 38.5548 16.5215 13.6691 第14年 31.2012 37.8428 15.7883 15.1677 第15年 30.9142 37.1645 15.2828 16.6386 第16年 30.5104 36.5090 14.8976 18.0831 第17年 30.0504 35.8705 14.5770 19.5021 第18年 29.5655 35.2462 14.2920 20.8963 第19年 29.0719 34.6341 14.0276 22.2663 第20年 28.5781 34.0335 13.7760 23.6124

附录3

用水文年支流的数据可以求出支流未来十年水质的变化情况：

状态1(％) 状态2(％) 状态3(％) 状态4(％) 第11年 34.5207 34.9669 14.9693 15.4430 第12年 33.7418 34.0741 13.6514 18.4326 第13年 32.7668 33.0698 12.8704 21.1930 第14年 31.7463 32.0320 12.3119 23.8098 第15年 30.7276 31.0009 11.8524 26.3191 第16年 29.7294 29.9925 11.4410 28.7372 第17年 28.7585 29.0125 11.0566 31.0724 第18年 27.8174 28.0628 10.6903 33.3295 第19年 26.9061 27.1435 10.3383 35.5121 第20年 26.0244 26.2539 9.9988 37.6229

14

本文来源：https://www.bwwdw.com/article/nm5a.html