山西省滹沱河水环境综合质量马尔柯夫变化趋势研究.docVIP

山西省滹沱河水环境综合质量马尔柯夫变化趋势研究 　　摘要 采用马尔柯夫水质动态分析法研究了山西省滹沱河7个断面的水环境综合质量的动态变化，结果表明：滹沱河7断面水环境综合质量的进步度指数的变化，从小到大依次为济胜桥＜定襄桥＜忻口桥＜乔儿沟＜代县桥＜崞阳桥＜南庄。滹沱河水质环境综合质量在空间上有明显的地域性，忻口桥至济胜桥之间有明显的峰值，是污染较严重河段，整体呈现出中间水质污染严重两头水质良好的趋势。滹沱河水环境综合质量的时间变化趋势是局部有所改善，整体趋于污染加重，Ⅴ类水质的范围在扩大，Ⅲ类水质的范围在缩小。滹沱河水质环境综合质量受产业结构和生产方式等因素的综合影响。 　　关键词 滹沱河；水环境综合质量；马尔柯夫 　　中图分类号 X824文献标识码A文章编号1007-5739（2008）20-0304-03 　　 　　马尔柯夫过程是研究某一事件的状态及状态之间转移规律的随机过程，具有描述事物随机变化的良好特性。它通过对事件初始状态（即t0时刻）的概率及状态间的概率转移关系进行研究，揭示事件发展的下一状态（即t0+Δt时刻）的变化趋势[1]。20世纪90年代中期，马尔柯夫方法被应用于河流综合水质的预报[2]，在实践中马尔柯夫预测方法与灰色预测相结合，在世界能源产量的预测方面显示出独特优势[3]。在环境质量预测方面，运用马尔柯夫过程的原理预测水环境质量转化的概率以及大气环境综合质量的动态变化研究颇多[4-9]。鉴于河流水质变化具有很强的随机性，变化态势只与当前的某种状态有关，所以马尔柯夫方法在水质动态分析方面被广泛应用。本文应用马尔柯夫方法对滹沱河水质的动态变化进行分析，目的在于了解滹沱河水环境质量的时空变化特点，为河流水资源保护和合理利用提供参考依据。 　　 　　1滹沱河基本概况 　　 　　滹沱河发源于山西省繁峙县的泰戏山，经代县、五台山北麓在原平、定襄与忻州交界处，受阻于金山，折向东流，至盂县出境，穿过太行山脉入河北省境。它在山西境内的长度为319km，流域面积18 856km2，主要支流有阳武河、云中河、牧马河和清水河等。河流两岸有平原、河谷、丘陵、山地等地貌，地形复杂多变。 　　滹沱河流域属于暖温带，温带大陆性季风气候，夏季暖热多雨，冬季寒冷干燥。据流经各县的气象数据记载，年平均气温8.0～9.0℃，最热月7月平均气温22.5～23.5℃，最冷月1月平均气温-10.0～-9.0℃，10℃积温3 100～2 400℃，年降水量400～450mm。滹沱河的生物多样性十分丰富，是湿地生物的优越栖息地[10]。 　　滹沱河水量较大，含沙量较小，水力蕴藏量为18.17万千瓦，开发利用价值较大。沿岸农牧业比较发达，20世纪80年代中期以来工业产业快速成长，目前已成为山西中部重要的工农业产业经济带。流域矿产资源丰富，以煤、铝、铁为主，矿种达50余种，形成了以煤炭、电力、建材、冶金、化工为主的地域性产业结构。20世纪90年代以来，滹沱河受工、农业污染的威胁越来越重，引起社会的广泛关注。 　　 　　2研究方法 　　 　　2.1滹沱河断面及其水质资料 　　根据忻州市环境保护局最近4年在山西省滹沱河从上游向下游设置的7个断面（图1）的监测数据，按照《地表水环境质量标准》（GB3838--2002）将滹沱河水环境综合质量划分为5个类别（即Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ类），部分断面的水环境综合质量劣于Ⅴ类水，即劣Ⅴ类水（表1）。 　　2.2马尔柯夫概率转移矩阵的构建 　　由表1可知，该矩阵有6个状态，即S={1，2，3，4，5，6}。以崞阳桥断面水环境综合质量类别结果为原始矩阵，计其他各断面的概率转移矩阵构建方法与之相同，在此不再赘述。 　　2.3进步度的确定及水环境综合质量动态评价 　　定义1：Sk=（Sijk）=（（i-j）3Pijk），i，j∈S，称矩阵Sk为概率转移矩阵Pk的进步矩阵。（i-j）3为水质改善或恶化的权重，水质改善权重为正，恶化为负，水质不变权重为0。 　　定义2：pd（Sk）=∑Sijk=（i-j）3Pijk，称pd（Sk）为该断面的水质综合质量进步度。若pd（Sk）＞0，说明第k个断面的水质综合质量状况得到改善；pd（Sk）＜0，说明第k个断面的水质综合质量状况有所恶化。若pd（S1）＜pd（S2）＜……＜pd（Sn），是指从S1断面到Sn断面中，断面的水质综合质量变化情况相对最好，S1断面的水质综合质量变化情况相对最差。 　　 　　3结果与分析 　　 　　3.1水环境综合质量的进步度特点 　　由表2可以看出，滹沱河各断面水环境综合质量的进步度指数变化趋势从小到大依次为济胜桥＜定襄桥＜忻口桥＜乔儿沟＜代县桥＜崞阳桥＜南庄，即南庄的水环境综合质量变好，济胜桥、定襄桥、忻口桥水环境