滇池高原内陆湖泊水环境现状与治理对策.docxVIP

滇池高原内陆湖泊水环境现状与治理对策 云南省最大的内陆湖泊，水量约15.7亿3m。它是云南省重要的水源。近年来，污染逐渐加剧。目前，整个湖泊的水质属于较差的v类。虽然已经投入大量资金进行管理，但效果并不明显。 本文应用时空地质统计学模拟滇池水体水质的时空分布, 强调水质分布的时空位置和适度外推预测, 用时空变异函数揭示水质指标的时空变异规律, 并据此用时空普通克立格方法估计其水质指标浓度在不同时间的空间分布。 1 草海海域水质指标的频率直方图 滇池位于昆明市的西南边, 有29条河道汇入其中, 是整个昆明生产和生活污水的接纳体。滇池整体呈弓形凸向东侧, 以海埂为界分为南北两个部分, 北部为草海, 南部为外海, 整个水域约300km2, 分布有10个观测站 (草海2个, 外海8个) , 见图1。由于草海水域狭长且与外海仅有一条很窄的航道连通, 其所得的测量具有很大的特殊性, 从水质指标的频率直方图来看, 其明显与外海的测量数据不属同一分布, 因此本文仅选择使用了外海8个观测站的数据, 包括总氮 (TN) 、总磷 (TP) 、高锰酸盐指数 (CODMn) 、溶解氧 (DO) 、生化需氧量 (BOD) 、叶绿素a (Chl-a) 等6个指标在2005～2006年每个月的观测数据, 共192×6个数据。 从6个指标与时间或空间的位置关系来看, 未见其有明显的依赖于时间或空间的规律性联系, 同时各个指标在经自然对数转换后多服从单峰正态分布 (TP除外) , 因此这些水质指标可以看作时空域内的随机变量。TP的分布不具有单峰正态分布特征, 而是呈双峰特征, 表明磷至少有两个不同的来源, 这与磷主要来源于滇池东岸的农业面源污染和滇池南岸的昆阳磷矿的事实相符。 2 时空变异函数的应用 实际上, 自然界中的多数变量, 尤其是环境变量, 均一定程度地存在依赖时间、空间变化的规律, 可以用时空变异函数(Space-time Semivariogram) 来描述。 时空变异函数γ (hs, ht) 是描述时空变量Z (s, t) 与时空变化 (hs, ht) 之间规律的一种工具, 在Z (s, t) 满足二阶平稳(Second-order stationary) 假设的条件下, 其定义为: 2γ(hs,ht) =E[Z(s+hs,t+ht) -Z(s,t) ]2(1) 其中:Z (s, t) 为时空点 (s, t) 处的值, hs、ht分别为空间、时间间隔变量, s= (s1, s2, …, sd) ∈D?Rd(一般地d=3) 代表空间坐标, t∈T?R为时间坐标, 且D×T=Rd+1。 据此可计算实验时空变异函数, 而理论时空变异函数则一般用乘积模型、乘积-和模型、张量积三次光滑曲面法等及其修正方法来模拟, 其中的乘积-和模型为本文使用的模拟方法, 其定义如下: γ(hs,ht) =γ(hs, 0) +γ(0,ht) -K×γ(hs, 0) ×γ(0,ht) (2) 其中?K=sill[γ(hs,0)]+sill[γ(0,ht)]?sill[γ(hs,ht)]sill[γ(hs,0)]×sill[γ(0,ht)]其中?Κ=sill[γ(hs,0)]+sill[γ(0,ht)]-sill[γ(hs,ht)]sill[γ(hs,0)]×sill[γ(0,ht)] 且K必需满足充分必要条件0K≤1max{sill[γ(hs,0)];sill[γ(0,ht)]}?sill0Κ≤1max{sill[γ(hs,0)];sill[γ(0,ht)]}?sill为变异函数的基台值。 应用上述乘积-和模型, 可以求得时空变异函数的理论模型, 据此可以计算任意时空点上 (一般是变程以内) 的时空变异函数值。图2列出了其中3种水质指标的实验时空变异函数及其积-和理论模型。从计算结果来看, 时间边界变异函数的拟合度比空间边界变异函数更好, 其拱高明显较大, 具有更好的连续性。表明水质指标随时间变化的变异规律较明显, 而随空间变化的变异规律则较差。 根据滇池水体的空间展布形态, 需考虑其空间变异函数的各项异性, 即分别计算NS和EW向的实验变异函数。从计算结果来看, NS方向和EW方向的空间变异函数的差异不大, 且由于本文中所用数据仅限于样本数很少的8空间点, 因此本文未考虑空间变异函数的各向异性。 为了检验时空变异函数在应用时空克立格插值方法模拟滇池水质时空分布中的有效性, 需要计算其交叉验证系数。应用上述理论时空变异函数计算了6种指标各观测时空点 (192×6) 的估计值 (其中参与克立格估计的样品点为距待估点最近的时间、空间点各3个, 但不包括待估点本身) , 求出了各样品时空点的实测值与估计值之间的相关系数。从结果来看, TN (0.7) 、Mn (0.76)