太湖流域水质的预测与防治策略
- 格式:doc
- 大小:625.50 KB
- 文档页数:19
太湖流域水质的预测与防治策略摘要本文在参考一些数据和环境资料的基础上用“移动平均数预测法”对太湖流域的水质情况进行了详尽的分析研究,且对前几个周的水质情况给出了适当的预测,让大家更加深刻地认识到目前太湖流域严峻的环境问题。
对于问题一,我们根据“移动平均数预测法”预测出第1—6周时间内断面7斜路港河(苏-浙省界)的主要监测项目(PH 、DO 、Mn COD 、3NH N )的具体数值及水质类别;然后我们根据各组数据拟合得到各监测项目的变化曲线;最后我们用线性加权对各个监测项目加权,最后我们画出断面7的水质变化曲线;各表及曲线图见附件二。
对于问题二,我们根据移动平均数预测法利用第16—10周的数据推测出第9周各断面的各监测项目的数值;同理用16—8周的数据可以得到第7周的水质状况表;各表见附件三。
对于问题三,与问题二的方法一样我们可以预测出第18、20周的水质预测项目数据,再根据《地表水环境质量标准》得到其水质状况表见附件四,最后我们对预测精确度进行了灵敏度分析,分析知道e s 为68%、2e s 为95%。
关键字:移动平均数预测法 拟合 线性加权 预测精确度一、问题重述1.1 太湖流域简介太湖流域面积236500km,其以平原为主,占总面积的4/6,水面占1/6,其余为丘陵和山地。
三面滨江临海,西部自北而南分别以茅山山脉、界岭和天目山与秦淮河、水阳江、钱塘江流域为界。
地形特点为周边高、中间低。
中间为平原、洼地,包括太湖及湖东中小湖群、湖西洮鬲湖及南部杭嘉湖平原,西部为天目山、茅山及山麓丘陵。
北、东、南三边受长江口及杭州湾泥沙淤积的影响,形成沿江及沿海高地,整个地形呈碟状。
1.2 相关情况为了能够及时全面的掌握主要流域重点断面水体的水质状况,预警或预报重大(流域性)水质污染事故,需要对水质进行实时连续监测和远程监控。
为此,从1999年9月开始至2003年12月,国家环保总局在松花江、辽河、海河、黄河、淮河、长江、珠江、太湖、巢湖、滇池等流域建设了82个水质自动监测站。
监测项目为五参数(水温、pH、DO、TB及EC)、高锰酸盐指数、氨氮和TOC八项指标。
附表一所给的数据是太湖流域7个断面2005年第1周至第16周的具体数据。
水质的分类见附表二1.3 问题的提出1. 在第1-6周断面7:斜路港河(苏-浙省界),由于某种原因未获得数据,请根据有关数据推测该断时间内断面7的水质变化曲线;2. 因数据处理过程中的失误,将第7周、第9周的数据遗失,试通过数学知识模拟该时间段的数据;3. 根据所给数据预测2005年第18周、第20周的断面水质状况。
二、问题分析问题一,要得到水质变化曲线就要知道1—6周断面7的主要监测项目的数据。
我们运用移动平均数预测法来推测出各主要监测项目的数据。
将各个监测项目所得的数据拟合,得到拟合曲线;再用线性加权得到各个监测项目的权重,我们可以画出断面7的水质变化曲线。
问题二,分析附表一给出的第16—10周的水质状况表中各断面的PH 的数据,用移动平均数预测法推测出第9周各断面PH 的数值;同理可得到各断面DO 、Mn COD 、3NH N -的数值;进而得到第9周的水质状况表。
运用同样的方法我们可以推测出第7周各断面的PH 、DO 、Mn COD 、3NH N - 的数值,从而得到第7周的水质状况表。
问题三,在问题二的基础上用同样的方法可以预测出第18、20周的水质状况。
三、模型假设1、题中所给的数据是正确的2、水质变化没有明显的季节性四、符号约定1、i x :时间序列的观察量(1,2,,,i t =⋅⋅⋅⋅⋅⋅)2、t x :在周期t 内预测的时间序列值3、N :观察的次数4、t e :t x 的预测值误差5、MAD :所有t e 的平均绝对值6、e s :预测误差的标准差五、问题求解5.1移动平均数预测法我们定义,1t f 为观察t x 后制定的周期1t +的预测周期。
在移动平均法中,1t f =最后N 次观察的平均值=t x ,1t x -,2t x -…1t N x -+的平均值,其中N 是指定的参数。
我们取N=3,针对前6周的数据预测各个断面各主要监测项目的变化状况。
现在对断面1江苏无锡沙渚的监测项目PH 的前6个月的数据移动平均分析:在第16—14周中,我们还没有观察到3个月的数据,所以当N=3时,我们无法计算出这3个月的移动平均指标预测。
对于第13周求得预测值13,17.927.897.97.9033f ++==,对于第12周,求得预测值12,17.897.97.857.883f ++==,对于第11周求得预测值12,17.97.857.887.883f ++==。
表1PH 断面1江苏无锡沙渚的移动平均预测值(N=3)注意我们的预测从一个周期到下一周期的“移动”,用最近观察的数据代替“最老”的平均观察数据。
N 的选择我们应当如何选择N ,使预测误差最小,即如何计算移动平均数的周期数量?为了更好的说明问题,我们定义二个数学名词:量度:平均绝对偏差(MAD )作为预测精确度的量度 预测误差t e :()t t t e x x =-的预测值 由表1可知:137.857.9030.053e =-=-,127.887.880e =-=,117.957.880.13e =-=,1312110.05300.130.06133e e e MAD ++-++===因此平均来说监测项目PH 的每周偏差为0.061。
我们试图把下周的预测值作为最后N 个周实际数值的平均数,即N 取什么值才会使我们的平均绝对误差最小呢?我们取1,2,,12N =⋅⋅⋅.下面用 Excel OFFSET 函数进行统计说明:1)步骤一、B0复制到C6:C11,D 列为C 与B 列的差的绝对值,求D6:D11和然后除1,得到的值为0.23;2)步骤二、C7为B1与B2的和,列以此类推,D 列为C 与B 列的差的绝对值,求D7:D11和然后除1,得到的值为0.23;3)同理得数据如H 列所示,其中最小值为0.025,对应的MDA 为2 我们用MATLAB 可求各断面各监测项目的N 对应最小MAD ,程序见附件一。
结果如下表:5.2 问题一由移动平均数预测法得断面7第1—6周各监测项目的数据如下表:PH DO COD NH3-N 第6周 7.03 2.168511 3.782324 1.648518 第5周 7.03 1.454713 3.344683 1.77076 第4周 7.03 1.266773 2.509401 1.061161 第3周 7.03 0.978 1.927282 1.029687 第2周 7.03 1.173599 2.312738 0.918354 第1周7.031.4083192.7752861.071413断面1 PH DO CODNH3-N 1 0.038 2 0.296 6 0.15 1 0.257 断面2 5 0.079 6 0.047 1 0.3 6 0.272 断面3 5 0.222 1 0.412 1 0.6 4 0.065 断面4 2 0.189 1 0.397 1 0.8 3 0.158 断面5 6 0.023 2 0.57 6 0.133 4 0.414 断面6 6 0.093 5 0.056 4 0.842 1 0.253 断面710.27551.245 0.71 60.17由PH的1—6周的数据拟合曲线如下图:周数17.03y=由DO的1—6周的数据拟合曲线如下图:周数20.0017 6.5373y x=+由MnCOD的1—6周的数据拟合曲线如下图:周数230.00480.047.147y x x =-+由3NH N -的1—6周的数据拟合曲线如下图:周数240.00070.0004 1.656y x x =++我们用线性加权对各个监测项目加权,即212340.10.20.40.30.002130.156 3.53y y y y y x x =+++=-+。
最后我们画出断面7的水质变化曲线如下图:周数2=-+y x x0.002130.156 3.535.3 问题二由移动平均数预测法模拟得第7、9周的数据,再根据水质类别定义知道各断面的水质类别,水质状况表如下:第7周(2月7日-2月13日)水质类别PH DO COD NH3-N 本周上周江苏无锡沙渚7.955 9.628 5.197 0.88 ⅢⅢ江苏宜兴山嘴 6.654 10.789 5.3 0.787 ⅤⅢ江苏苏州西山7.16 11.6 2.6 0.42 ⅡⅡ浙江湖州新塘港 6.907 9.17 5.5 0.694583 ⅣⅣ上海青浦急水港 4.77933 10.385 3.0995 0.694583 劣Ⅴ劣Ⅴ浙江嘉兴王江泾 5.57589 10.3144 4.124875 1.82 ⅤⅣ浙江嘉兴斜路港7.03 10.45168 4.124875 0.882153 ⅢⅤ第9周(2月21日-2月27日)水质类别PH DO COD NH3-N 本周上周江苏无锡沙渚7.97 9.415 5.083 1.49 ⅢⅢ江苏宜兴山嘴 6.602 11.067 1.8 0.62 ⅢⅢ江苏苏州西山 2.9144 11 2.2 0.5275 ⅡⅡ浙江湖州新塘港 4.7582 9.22 3 0.879167 ⅢⅢ上海青浦急水港 3.70743 10.11 2.013833 0.879167 劣Ⅴ劣Ⅴ浙江嘉兴王江泾 4.32534 10.1624 2.253458 2.33 ⅤⅤ浙江嘉兴斜路港 6.51 10.31188 2.253458 1.120972 ⅣⅤ5.4 问题三由移动平均数预测法推测得第18、20周的数据,再根据水质类别定义知道各断面的水质类别,水质状况表如下:第18周(4月25日-5月1日)水质类别 PH DO COD NH3-N本周 上周 江苏无锡沙渚 7.9125 8.7 4.9917 0.23 Ⅲ Ⅱ 江苏宜兴山嘴 6.59 10.1022 3.2 0.4403 Ⅱ Ⅲ 江苏苏州西山 7.51 9.31 4.8 2.38316 劣Ⅴ Ⅱ 浙江湖州新塘港 8.5275 7.48 5.8 0.892222 Ⅲ Ⅲ 上海青浦急水港 7.18333 5.6525 5.872222 5.64375 劣Ⅴ 劣Ⅴ 浙江嘉兴王江泾 6.50111 2.3864 7.0375 1.25 Ⅳ Ⅳ 浙江嘉兴斜路港 7.06 2.516 6.2 1.746667Ⅴ Ⅳ第20周(5月2日-5月8日) 水质类别 PH DO COD NH3-N本周 上周 江苏无锡沙渚 7.91063 8.705 4.7803 0.23 Ⅲ Ⅲ 江苏宜兴山嘴 6.6488 10.12173 3.2 0.377323 Ⅱ Ⅱ 江苏苏州西山 7.52258 9.31 4.8 0.224375 Ⅲ Ⅲ 浙江湖州新塘港 8.55188 7.48 5.8 0.915062 Ⅲ Ⅳ 上海青浦急水港 7.20898 5.7375 5.970525 6.028359劣Ⅴ 劣Ⅴ 浙江嘉兴王江泾 6.57401 2.394416 7.521094 1.25 Ⅳ Ⅳ 浙江嘉兴斜路港 7.06 2.71064 5.42 1.697685Ⅴ Ⅳ六、预测精确度的灵敏度分析通过观察可知预测误差服从正态分布的模型,我们可以使用MAD 来估计e s =预测误差的标准差。