贵阳市地下水水质评价

2024-12-15

贵阳市地下水水质评价（精选8篇）

贵阳市地下水水质评价篇1

地下水水质现状评价是地下水资源评价和保护的重要内容,通过对地下水水化学资料的分析,可以科学地评价区域地下水体的质量状况。可为防治评价区内水质恶化和制定水资源管理决策方案提供科学依据,对评价区内地下水资源的可持续开发利用和综合管理有很好的指导作用[1]。目前地下水水质评价的方法有很多,主要有属性识别、综合指数法、人工神经网络模型、模糊综合评判法、灰色聚类法等。由于受到诸多因素的影响和制约及方法自身存在的不足,至今仍没有一个广泛认可的评价模型。如综合指数法中权重不好确定和忽略了水质分级界线的模糊性;人工神经网络需要的样本数多, 网络结构的优劣因人而异;模糊综合评判中隶属函数和权重矩阵的构造;灰色聚类法中白化函数的选择和聚类权的确定往往因人而异,造成评价模式难以通用。本文在陈守煜教授可变模糊集理论基础上[2,3,4,5],提出地下水水质评价的可变模糊评价方法。该方法能够科学、合理地处理水质分级界线的模糊性,确定样本指标对各级指标标准区间的相对隶属度,为了增强二元比较分析法确定权重的可操作性,使其应用更为便捷,提出二元比较模糊决策分析简捷方法,以此确定各指标权重;此外,该方法能够通过变化模型及其参数,全面地给出样本的评价等级,提高对样本等级评价的可信度。

1 相对差异函数模型

设X0=[a,b]为实轴上模糊可变集合 $\underset{˜}{V}$ 的吸引域,即 $0 < D_{\underset{˜}{A}} \leq 1$ 区间,X=[c,d]为包含X0(X0⊂X)的某一上、下界范围域区间。如图1所示。

根据模糊可变集合 $\underset{˜}{V}$ 定义可知[c,a]与[b,d]均为 $\underset{˜}{V}$ 的排斥域,即 $- 1 \leq D_{\underset{˜}{A}} (u) ＜ 0$ 区间。设M为吸引域区间[a,b]中 $D_{\underset{˜}{A}} (u) = 1$ 的点值,M不一定为区间[a,b]的中点值,需按物理分析确定。x为X区间内的任意点的量值,则当x落入M点左侧时,相对差异函数模型可为:

${\begin{matrix} D_{\underset{˜}{A}} (u) = (\frac{x - a}{Μ - a})^{β} & x \in [a, Μ] \\ D_{\underset{˜}{A}} (u) = - (\frac{x - a}{c - a})^{β} & x \in [c, a] \end{matrix} (1)$

x落入M点右侧时,其相对差异函数模型为:

${\begin{matrix} D_{\underset{˜}{A}} (u) = (\frac{x - b}{Μ - b})^{β} & x \in [Μ, b] \\ D_{\underset{˜}{A}} (u) = - (\frac{x - b}{d - b})^{β} & x \in [b, d] \end{matrix} (2)$

x落入X区间外时

$D_{\underset{˜}{A}} (u) = - 1 u \notin [c, d] (3)$

公式(1)、(2)中β为非负指数,常取β=1,即相对差异函数模型为线性函数。

$D_{\underset{˜}{A}} (u)$ 确定以后,根据公式(4)求解相对隶属度 $μ_{\underset{˜}{A}} (u)$ 。

$μ_{\underset{˜}{A}} (u) = \frac{1 + D_{\underset{˜}{A}} (u)}{2} (4)$

从以上分析可知,当x=a、x=b时, $μ_{\underset{˜}{A}} (u) = 0.5$ ; 当x=M时, $μ_{\underset{˜}{A}} (u) = 1$ ;当x=c、x=d时, $μ_{\underset{˜}{A}} (u) = 0$ ;当x∉[c,d]时,满足 $μ_{\underset{˜}{A}} = 0$ 。

2 二元对比确定指标权重方法

二元对比模糊决策分析法是以我国语言与思维习惯为基础,以模糊集的余集定义为求解手段,严格遵循一致性检验条件,针对无结构决策问题而提出的一种有效简便的求解方法。该方法数学理论严谨,多次试验数据表明确定的权重有代表性,依据专业知识确定指标重要性二元对比后即可进行权重的确定。为了增强二元对比分析法确定权重的可操作性,使其应用更为便捷,本文在系统介绍二元比较模糊决策分析基础上,引申出其简捷方法。

2.1 一致性检验

二元比较判断矩阵的一致性问题是非结构决策合理与科学性的基础,被研究的对象组,若在属性排序明确的前提下不需要进行一致性检验;在最重要指标明确时,为简化操作,可以不进行一致性检验,直接按简化方法求解;若各研究对象属性排序模糊不定的情况下,必须进行一致性检验,具体的检验原则可参考文献[6]。

一致性检验后,所得模糊标度矩阵各行和数由大到小排列,此排序即为指标集在满足一致性条件下关于重要性的排序,具体证明可见[7]。

2.2 定性指标权重的确定

指标集D就模糊概念 $\underset{˜}{A}$ 重要性所给的排序作二元比较,若二元比较矩阵为:

$\begin{array}{l} B = (\begin{matrix} b_{11} & b_{12} & \dots & b_{1 m} \\ b_{21} & b_{22} & \dots & b_{2 m} \\ ⋮ & ⋮ & ⋱ & ⋮ \\ b_{m 1} & b_{m 2} & \dots & b_{m m} \end{matrix}) = (b_{g k}) (5) \\ g = 1, 2, \dots, m; k = 1, 2, \dots, m \end{array}$

B称为指标集D关于A对重要性二元比较定量排序一致性模糊标度矩阵,bkg表示指标g与指标k相对于重要性的二元比较。

由于对重要性二元比较定量排序矩阵一致性模糊标度矩阵的对角线元素表示自身的比较,其元素值为0.5,它将元素分为上下两个三角,上三角中的元素值从对角线元素0.5开始,每行元素值自左向右递增,每列元素自下向上递增,因此上三角矩阵自对角线元素0.5开始的每一行元素满足:

$0.5 = b_{g k} \leq b_{g (k + 1)} \leq b_{g (k + 2)} \leq \dots \leq b_{g m} \leq \dots \leq 1 (6)$

由互补性条件[8]得到下三角矩阵,见式(7)。

$b_{k g} = 1 - b_{g k} (7)$

根据矩阵B构造相及矩阵Φ:

$Φ = (\begin{matrix} 1 & b_{12} / b_{21} & \dots & b_{1 m} / b_{m 1} \\ b_{21} / b_{12} & 1 & \dots & b_{2 m / m 2} \\ ⋮ & ⋮ & ⋱ & ⋮ \\ b_{m 1} / b_{1 m} & b_{m 2} / b_{2 m} & \dots & 1 \end{matrix}) (8)$

重要性的相对比较是权重的一种测度,按照模糊集的特点测度上限为1,故有序相及矩阵Φ的上三角元素均定义为1,对矩阵Φ每一行取小,即Φ的第一列元素。

$\begin{array}{l} ω^{´} = (ω_{1}^{´}, ω_{2}^{´}, \dots, ω_{m}^{´}) = (1, b_{21} / b_{12}, \dots, b_{m 1} / b_{1 m}) = \\ (1, \frac{1 - b_{12}}{b_{12}}, \dots, \frac{1 - b_{1 j}}{b_{1 j}}, \dots, \frac{1 - b_{1 m}}{b_{1 m}}) (9) \end{array}$

式(9)给出了指标集对重要性的相对隶属度向量即非规一化权向量。对其规一化处理,得指标权向量:

$ω^{´} = (ω_{1}^{´}, ω_{2}^{´}, \dots, ω_{m}^{´}) = (ω_{1}^{´} / \sum_{i = 1}^{m} ω_{i}^{´}, ω_{2}^{´} / \sum_{i = 1}^{m} ω_{i}^{´}, \dots, ω_{m}^{´} / \sum_{i = 1}^{m} ω_{i}^{´}) (10)$

按照我国的思维习惯确定语气算子,根据公式(6)的两个边界值0.5与1,在同样与无可比拟之间插入语义逐渐加重的语气算子,相应模糊标度在0.5与1.0之间以线形增值0.05插入模糊标度值。已知模糊标度根据式(9)计算可得相应相对隶属度,为便于应用模糊概念语气算子与模糊标度值、相对隶属度间的关系归入表1。也可根据实际情况,将表1进一步细化。相对隶属度即非规一化权重已知的前提下依据式(10)确定指标权重。

3 地下水水质模糊可变评价方法

3.1 方法的建立

设有n个待识别的水质评价对象

$X = {x_{1}, x_{2}, \dots, x_{n}} (11)$

评价对象j的特性用m个指标特征值来表示

$x_{j} = (x_{1 j}, x_{2 j}, \dots, x_{m j})^{Τ} (12)$

则n个评价对象的全体指标特征可用m×n阶矩阵表示

$X = (x_{i j}) (13)$

式中:xij为第j个评价对象第i个指标的特征值;i=1,2,…,m;j=1,2,…,n。评价对象依据m个指标按c个级别的指标标准特征值进行识别,则有m×c阶指标标准特征值矩阵

$Y = (y_{i h}) (14)$

式中:yih为级别h指标i的标准特征值,h=1,2,…,c。

参照指标标准值矩阵和实际情况确定水质评价可变集合的吸引(为主)域矩阵与范围域矩阵

$Ι_{a b} = ([a_{i h}, b_{i h}]), Ι_{c d} = ([c_{i h}, c_{i h}]) (15)$

根据各级别的实际情况确定水质评价吸引域[aih,bih]中 $D_{\underset{˜}{A}} (x_{i j})_{h} = 1$ 的点值Mih的矩阵

$Μ = (Μ_{i h}) (16)$

根据式(1～3)计算相对差异度,再由式(4)计算指标对h级的相对隶属度矩阵

$U_{h} = [μ_{\underset{˜}{A}} (x_{i j})_{h}] (17)$

设m个指标的权向量为

$\hat{w} = \begin{matrix} (w_{1} & w_{2} & \dots & w_{m}) \end{matrix} = w_{i} (18)$

满足 $\sum_{i = 1}^{m} w_{i} = 1$ 。

参考连续统上任一点x指标特征值i的相对隶属度 $μ_{\underset{˜}{A}} (x_{i j})_{h}$ 和 $μ_{\underset{˜}{Ζ}^{c}} (x_{i j})_{h}$ 关于左、右极点的广义权距离分别为:

$\begin{array}{l} d_{g} = {\sum_{j = 1}^{m} [ω_{j} (1 - μ_{\underset{˜}{A}} (x_{i j})_{h})]^{p}}^{1 / p} (19) \\ d_{b} = {\sum_{j = 1}^{m} [ω_{j} (1 - μ_{\underset{˜}{A}^{c}} (x_{i j})_{h})]^{p}}^{1 / p} = {\sum_{j = 1}^{m} [ω_{j} μ_{\underset{˜}{A}} (x_{i j})_{h})]^{p}}^{1 / p} (20) \end{array}$

模糊可变评价模型为:

$_{i} u_{h}^{´} = \frac{1}{1 + (d_{g} / d_{b})^{a}} (21)$

式中:iu′h为评价对象关于等级概念A对吸引性质 $\underset{˜}{A}$ 的非归一化相对隶属度;α为模型优化准则参数,α=1为最小一乘方准则,α=2为最小二乘方准则;p为距离参数,p=1为海明距离,p=2为欧式距离;其中j为评价指标序号,j=1,2,…,m。

通常情况下模型(21)中α和p有4种搭配:

$α = 1, p = {\begin{cases} 2 \\ 1 \end{cases} ‚ α = 2, p = {\begin{cases} 2 \\ 1 \end{cases} (22)$

(1)当α=1,p=1时,式(21)变为:

$_{i} u_{h}^{´} = \sum_{j = 1}^{m} [ω_{j} μ_{\underset{˜}{A}} (x_{i j})_{h}] (23)$

式(21)为一模糊综合评价模型,是一个线性模型,属于模糊可变评价模型的一个特例。

(2)当α=1,p=2时,式(21)变为:

$_{i} u_{h}^{´} = \frac{d_{b}}{d_{b} + d_{g}} (24)$

此时式(21)为理想点模型,属于模糊可变评价模型的又一个特例。

(3)当α=2,p=1时,式(21)变为:

$\begin{array}{l} _{i} u_{h}^{´} = \frac{1}{1 + [(1 - d_{b}) / d_{b}]^{2}} (25) \\ d_{b} = \sum_{j = 1}^{m} [ω_{j} μ_{\underset{˜}{A}} (x_{i j})_{h}] (26) \end{array}$

式(21)为Sigmoid型函数,可用以描述神经网络系统中神经元的非线性或激励函数。

(4)当α=2,p=2时,式(21)变为:

$\begin{array}{l} _{i} u_{h}^{´} = \frac{1}{1 + (d_{g} / d_{b})^{2}} (27) \\ d_{g} = \sqrt{\sum_{j = 1}^{m} {ω_{j} [1 - μ_{\underset{˜}{A}} (x_{i j})_{h}]}^{2}} (28) \\ d_{b} = \sqrt{\sum_{j = 1}^{m} [ω_{j} μ_{\underset{˜}{A}} (x_{i j})_{h})]^{2}} (29) \end{array}$

此时模糊可变评价模型变为模糊优选模型。

由此可见,模糊可变评价模型是一个变化模型,在可变模糊集理论中是一个十分重要的模型,可广泛应用于模糊概念的识别问题。通过不同参数组合,可对评价成果的可靠性进行验证。

由可变评价模型可得非归一化的综合相对隶属度矩阵

$U^{´} = (_{i} u_{h}^{´}) (30)$

将式(30)归一化处理得到综合相对隶属度矩阵

$U = (_{i} u_{h}) (31)$

式中:

$_{i} u_{h} =_{i} u_{h}^{´} / \sum_{h = 1}^{c}_{i} u_{h}^{´} (32)$

级别特征值公式为:

$Η = (1, 2, \dots, c) \cdot U (33)$

3.2 求解步骤

(1)根据已有观测资料及评价标准(依据),确定待评价对象现状指标特征值与指标标准值矩阵X、Y;

(2)参照指标标准值矩阵Y结合实际情况确定水质评价可变集合的吸引(为主)域矩阵与范围域矩阵Iab、Icd;

(3)结合实际情况确定水质评价吸引域Iab中 $D_{\underset{˜}{A}} (x_{i j})_{h} = 1$ 的点值的矩阵M;

(4)根据式(1)～(3)计算相对差异度矩阵,再由式(4)计算指标对h级的相对隶属度矩阵Uh;

(5)通过二元比较法确定指标的权向量;

(6)采取不同的参数准则,利用模糊可变评价模型(21),求得综合相对隶属度矩阵并作归一化处理;

(7)求级别特征值矩阵,并据此做出评价。

4 应用实例

4.1 实例应用1

以邯郸市化工区地下水水质资料为例[9],该区主要有镀锌厂、树脂厂和磷肥厂等,共有7个水井测点,每一测点测量水的硬度、SO $_{4}^{- 2}$ 、CL-、F和有机磷5项评价指标,依据水质所受污染的情况,将其分为,I级为污染水,II级为已污染,III级为严重污染。各指标现状指标特征值与指标标准值见表2。

mg/L

根据表2可得邯郸市化工区地下水水质的现状指标特征值与指标标准值矩阵X、Y,结合各指标物理意义确定吸引(为主)域矩阵Iab、Icd,其中指标为开区间的吸引域与范围域的确定以相邻指标区间距离作为参考距离。点值矩阵M为吸引域矩阵Iab中相对差异度等于1的点组成的矩阵,按物理意义分析可得。

五项指标重要程度与研究区域厂区结构有关,综合考虑各项因素确定CL-最为重要,以其作为标准与其他各指标进行二元比较分析,得到二元比较判断矩阵并作一致性检验,氟与磷均处于“同样”与“稍稍”边缘,硬度与SO $_{4}^{2 -}$ 在“稍稍”与“略微”之间,确定权重为ω=(0.179,0.161,0.242,0.219,0.198)。

$\begin{array}{l} X = \\ [\begin{matrix} 1 289 & 1 804 & 1 176 & 689 & 1 422 & 1 382 & 388 \\ 192.05 & 153.21 & 277.80 & 142.20 & 1 120.05 & 217.02 & 279.02 \\ 1 856.05 & 2 573.50 & 2 094.07 & 782.32 & 726.88 & 1 694.78 & 216.04 \\ 1.20 & 1.35 & 1.35 & 0.73 & 0.99 & 0.67 & 0.63 \\ 0.050 & 0.171 & 0.384 & 0.019 & 0.028 & 0.051 & 0.022 \end{matrix}] \\ Y = [\begin{matrix} < 250 & [250, 400] & > 400 \\ < 250 & [250, 400] & > 500 \\ < 250 & [250, 350] & > 350 \\ < 1.0 & [1.0, 1.5] & > 1.5 \\ < 0.05 & [0.05, 0.1] & > 0.1 \end{matrix}] \\ Ι_{a b} = [\begin{matrix} [0, 250] & [250, 400] & [400, 550] \\ [0, 250] & [250, 500] & [500, 750] \\ [0, 250] & [250, 350] & [350, 450] \\ [0, 1.0] & [1.0, 1.5] & [1.5, 2.0] \\ [0, 0.05] & [0.05, 0.1] & [0.1, 0.15] \end{matrix}] \\ Ι_{c d} = [\begin{matrix} [0, 400] & [0, 550] & [250, 550] \\ [0, 500] & [0, 750] & [250, 750] \\ [0, 350] & [0, 450] & [250, 450] \\ [0, 1.5] & [0, 2.0] & [1.0, 2.0] \\ [0, 0.1] & [0, 0.15] & [0.05, 0.15] \end{matrix}] \\ Μ = [\begin{matrix} 0 & 250 & 550 \\ 0 & 250 & 750 \\ 0 & 250 & 450 \\ 0 & 1.0 & 2.0 \\ 0 & 0.05 & 0.15 \end{matrix}] \end{array}$

判断矩阵Iab、Icd与M判断样本特征值xij在Mih点的左侧还是右侧,据此选用不同公式计算差异度矩阵及相对隶属度矩阵,采用不同参数组合根据模糊可变评价模型计算级别特征值,为增强结果的可比性,本次计算考虑了二元比较确定的权重及文献[9]中指标权重两种情况,分别为权重1权重2,取各次评价结果平均情况进行等级评价,等级判断标准H<1.67为I级,1.67≤H<2.5为II级,H>2.5为III级。评价结果见表3。将计算所得级别特征值从大到小排序,邯郸市化工区地下水水质污染程度排序依次为:测点3、测点2、测点5、测点1、测点6、测点4和测点7,可见测点3污染情况最为严重,测点7污染程度较轻。从表中可以看出,两组权重计算所得4组评价结果基本一致,级别特征值稳定在一个范围内,波动较小。不同参数组合下的评价结果稳定,从而说明了可变模糊评价法具有较高的可靠性。

将上述评价结果与文献[10]所采用的属性识别法及综合评价法,文献[9]所采用的改进TOPSIS评价法的计算结果进行比较,比较结果见表4。可变模糊评价法采用4种类型的结果均值作为评判标准,采用该法评价地下水水质是可行的、可靠的。另外,可变模糊评价法用级别特征值体现污染情况,有利于确切掌握水质污染的具体情况,以量化概念表征模糊事物进展程度,更加直观表征水质污染动态。

4.2 实例应用2

基于德州市水利局2004 年对德州市德城区地下水的监测数据[11],选取代表性的3个监测井观测数据进行水质评价,为增强可比性,水质标准分级延用原文中GB/T 14848-93标准,水质相关评价数据见表5。

利用地下水可变模糊评价法进行综合评判,评价结果及结果合理性比较见表6。可变模糊评价法采用4种类型的结果均值作为评判标准,该方法能够通过变化模型及指标权重,合理地确定出样本的评价等级,提高了样本等级评价的可信度。

mg/L

5 结语

本文利用相对差异函数表征地下水水质样本与评价等级间的吸引特性与排斥特性,系统分析二元对比模糊决策分析法特点的基础上引申出二元对比确定指标权重的简捷方法,使权重确定更简便易行,增强了可操作性,将以上理论与可变模糊集相结合提出了地下水水质评价的可变模糊评价方法,该方法集成理论的严谨性及操作的简捷性于一体。实例应用结果表明,可变模糊评价法在地下水水质评价中的评价是合理的,方法灵活性较强,评价结果可信度高。最后以量化的形式表征水质污染这一模糊事物的进展,有助于管理人员及时掌握水质动态,做出科学决策。该方法计算简便、思路清晰、计算结果合理,可推广到项目决策、识别和各种资源评价等方面,具有广泛的应用前景。

摘要：在可变模糊集理论基础上,结合二元比较确定指标权重方法,提出地下水质评价可变模糊评价方法,该方法能够合理地确定样本指标对各级指标标准区间的相对隶属度,并且能够通过变化模型及指标权重,确定出样本的评价等级,提高对样本等级评价的可信度。将该方法应用于地下水水质综合评价,评价结果表明,由于考虑了模型参数与权重的可变性,评价结果可信度高。

关键词：可变模糊评价,地下水,水质评价

参考文献

[1]苏耀明,苏小四.地下水水质评价的现状与展望[J].水资源保护.2007,23(2):4-9.

[2]陈守煜.工程模糊集理论与应用[M].北京:国防工业出版社,1998.

[3]陈守煜.工程可变模糊集理论与模型-模糊水文水资源学数学基础[J].大连理工大学学报,2005,45(2):308-312.

[4]陈守煜.水资源与防洪系统可变模糊集理论与方法[M].大连:大连理工大学出版社,2005.

[5]陈守煜.复杂水资源系统优化模糊识别理论与应用[M].长春:吉林大学出版社,2002.

[6]陈守煜.求解系统无结构决策问题的新途径[J].大连理工大学学报,1993,33(6):705-710.

[7]陈守煜.复杂水资源系统优化模糊识别理论与应用[M].长春:吉林大学出版社,2002.

[8]陈守煜.符合我国语言习惯的决策思维模式[J].大连理工大学学报(社会科学版),1999,20(2):8-10.

[9]张先起,梁川,刘慧卿.基于熵权的改进TOPSIS法在水质评价中的应用[J].哈尔滨工业大学学报,2007,39(10):1 670-1 672.

[10]舒军龙,潘仲麟.属性识别理论模型在环境质量评价中的应用[J].干旱环境监测,1998,14(4):224-227.

[11]王文强.综合指数法在地下水水质评价中的应用[J].水利科技与经济,2008,14(1):54-55.

贵阳市地下水水质评价篇2

萍乡市地下水水质监测与初步评价

随着工农业经济的发展和城市人口增加,水资源的保护和开发利用日趋重要,通过对萍乡市地下水成因类型的分析,并对地下水水质进行监测与评价发现,在火乐冲、制药厂、铝厂、萍钢福利厂等地段地下水阶段性地会受到轻度污染,为地下水易污染区,地下水质量属较差级,其它地段地下水质量良好.

作者：李锦才张炔郭伟作者单位：江西省九○一地质工程勘察院,江西,萍乡,337000刊名：西部探矿工程英文刊名：WEST-CHINA EXPLORATION ENGINEERING年，卷(期)：200921(9)分类号：P641关键词：地下水水质监测初步评价

葫芦岛市地下水水质评价篇3

1 评价依据与参数

地下水资源质量评价依据中华人民共和国国家标准《地下水质量标准》 (GB/T14848-93) 和中华人民共和国国家标准《生活饮用水卫生标准》 (GB5749-85) 。

地下水资源质量评价参数选用:pH、硫酸盐、氯化物、溶解性总固体、总硬度、硝酸盐氮、亚硝酸盐氮、高锰酸盐指数、氨氮、铁、锰、铜、锌、铅、镉、挥发酚、氰化物、六价铬、砷、汞、氟化物、阴离子合成洗涤剂等共22项。

2 评价方法

地下水资源质量评价:采用单指标评价方法, 按《地下水质量标准》所列分类指标将水质划分为五类, 不同类别标准值相同时, 从优不从劣, 以《地下水质量标准》III类水标准值的上限值确定为地下水控制标准, 单项指标水质最差为综合水质类别。

3 评价现状

3.1 按面积比

3.1.1 全市地下水水质

从评价结果看, 全市地下水水质较差, 超标率59.8%, Ⅲ、Ⅳ、Ⅴ类水质分别占全市总面积的40.2%、24.5%、35.3%。

从全市水资源各流域分区看, 大凌河、青龙河地下水质较好均达到Ⅲ类以上标准;其次为六股河, 超标率仅为13.7%;小凌河、小凌河~六股河、六股河~山海关水质最差, 超标率为100%。

3.1.2 山丘区、平原区地下水水质

山丘区:地下水水质较差, 超标率56.5%, Ⅲ、Ⅳ、Ⅴ类水分别占山丘区总面积的43.5%、18.0%、38.5%, 主要超标项目为氯化物、总硬度、溶解性总固体、铁、锰、氨氮、硝酸盐。其中大凌河、青龙河地下水质较好均达到Ⅲ类以上标准;其次为六股河, 超标率仅为13.6%;小凌河、小凌河~六股河、六股河~山海关水质最差, 超标率为100%。葫芦岛市山丘区地下水各类水质比重。

平原区:地下水水质较山丘区差, 超标率高达96.4%, Ⅲ、Ⅳ类水分别占平原区总面积的3.6%、96.4%, 主要超标项目为铁、硝酸盐氮、亚硝酸盐氮, 超标率分别为41.1%、9.9%、45.5%。从全市各水资源分区来看, 六股河超标率为23.1%, 小凌河~六股河、六股河~山海关超标率均为100%。葫芦岛市平原区地下水各类水质比重。

3.2 按水资源量

3.2.1 全市不同水质地下水资源量

经评价全市地下水水质较差超标率为77.5%。在全市地下水资源量4.18亿m3中, Ⅲ、Ⅳ、Ⅴ类水资源量分别为0.94亿m3、2.10亿m3、1.14亿m3, 分别占全市地下水资源总量的2 2.5%、5 0.2%、2 7.3%, 见表1、表2。

流域分区:大凌河、青龙河水质最好均达Ⅲ类以上标准;小凌河、小凌河~六股河、六股河~山海关水质最差、超标率为100%;六股河水质较好, 超标率仅为17.9%。

行政分区:建昌县水质最好均达Ⅲ类以上标准;南票区、连山区、龙港区、兴城市水质最差、超标率为100%;绥中县水质较差、超标率为8 0.0%。

3.2.2 山丘区、平原区不同水质地下水资源量

全市山丘区地下水资源量2.64亿m3, 其中Ⅲ、Ⅳ、Ⅴ类水资源量分别为0.90亿m3、0.60亿m 3、1.14亿m3。全市平原区地下水资源量1.70亿m3, 其中Ⅲ、Ⅳ类水资源量分别为0.08亿m3、1.62亿m3。

流域分区:大凌河地下水资源量0.26亿m3均为山丘区地下水资源量、水质为Ⅲ类。小凌河地下水资源量0.54亿m3均为山丘区地下水资源量, Ⅳ、Ⅴ类水分别占7.4%、92.6%。小凌河~六股河山丘区地下水资源量为0.76亿m3, Ⅳ、Ⅴ类水分别占57.9%、42.1%;平原区地下水资源量0.82亿m3均为Ⅳ类水。六股河山丘区地下水资源量为0.72亿m3, Ⅲ、Ⅳ类水分别占83.3%、16.7%;平原区地下水资源量0.11亿m3, Ⅲ、Ⅳ类水分别占72.7%、27.3%。六股河~山海关山丘区地下水资源量0.32亿m3均为Ⅴ类水;平原区地下水资源量0.77亿m3均为Ⅳ类水。青龙河地下水资源量0.04亿m3均为山丘区地下水资源量、水质为Ⅲ类。

行政分区:建昌县地下水资源量0.67亿m3均为山丘区地下水资源量、水质为Ⅲ类。南票区地下水资源量0.09亿m3均为山丘区地下水资源量、水质为Ⅴ类。连山区山丘区地下水资源量为0.67亿m3、水质为Ⅴ类;平原区地下水资源量0.10亿m3、水质为Ⅳ类。龙港区山丘区地下水资源量为0.06亿m3、水质为Ⅴ类;平原区地下水资源量0.06亿m3、水质为Ⅴ类。兴城市山丘区地下水资源量0.60亿m3、水质为Ⅲ类;平原区地下水资源量0.69亿m3、水质为Ⅳ类。绥中县山丘区地下水资源量为0.55亿m3, Ⅲ、Ⅴ类水分别占41.8%、58.2%;平原区地下水资源量0.85亿m3, Ⅲ、Ⅳ类水分别占9.4%、9 0.6%。

4 结语

贵阳市地下水水质评价篇4

加权优序法在地下水质评价中的应用

加权优序法是多目标评价与决策的方法之一,对加权优序法的`决策步骤作了详细地介绍,并据此对黑龙洞泉域水质进行了评价,结果表明用加权优序法进行水质评价,思路清晰,计算方法简捷,实用性好.

作者：任传栋王志真高佳刘成良 REN Chuan-dong WANG Zhi-zhen GAO Jia LIU Cheng-liang 作者单位：河北工程大学,河北,邯郸,056038刊名：水科学与工程技术英文刊名：WATER SCIENCES AND ENGINEERING TECHNOLOGY年，卷(期)：“”(5)分类号：X143关键词：加权优序法多指标评价水质评价 AHP

贵阳市地下水水质评价篇5

关键词：模糊模式识别法,地下水,相对隶属度,水质评价

水质评价就是按照评价目标,选择相应的水质参数、水质标准和评价方法,对水体的质量利用价值及水的处理要求作出评定。水质评价是合理开发利用和保护水资源的一项基本工作。地下水水质评价的主要内容包括pH值、总硬度分布特征、矿化度、水化学等类型分布。模糊模式识别法是由被誉为我国模糊水文水资源学科“工程模糊集理论、模型与应用”的创始人、大连理工大学陈守煜[2]教授提出的一个理论,在地表水水量评价中取得了良好的效果,现在该识别法已应用于水库水质评价[3]、船舶配载评价、地震岩性预测,甚至房地产估价等领域,所取得的结果比较令人满意。但是在地下水水质评价方面的应用还不多见,为此本文尝试用该法来评价地下水水质,并结合咸阳市实地观测资料来说明此种方法的可行性。

1 模式的原理

模糊模式识别法的中心思想是依据模糊数学提出了相对隶属度理论,并通过“距离平方和最小”的聚类准则以及模糊算子等评价程序的运用得出评价结果。其所需要引入和使用的参数有指标特征值矩阵和指标标准特征值矩阵,在此基础上建立起评判模型。

1.1 指标特征值矩阵和指标标准特征值矩阵

X为n个井位组成的样本集合,每个井有m项评价指标。则有实测指标矩阵:

其中,xij为第j个井第i个因子的实测值,i=1,2,…,m;j=1,2,…,n。

m项指标按5级环境标准评价,依据GB/T 14848-93Ⅲ地下水质量标准值则有指标标准矩阵:

其中,yih为h级标准指标i的标准值,h=1,2,…,5。

用相对隶属度(rij,sih)来描述环境污染的模糊性,规定:对实测值浓度大污染越严重的指标i,xij≤yi1的,对于“污染”的相对隶属度rij=0;xij≥yi5的,rij=1。实测值浓度大污染越轻的指标i,xij≥yi1的,对于“污染”的相对隶属度rij=0;xij≤yi5的,rij=1。介于yi1与yi5之间的,实测指标相对隶属度rij用式(3)表示:

$r_{i j} = \frac{x_{i j} - y_{i 1}}{y_{i 5} - y_{i 1}}$ (3)

类似我们可规定:指标i的1级标准值对“污染”的相对隶属度sih=0;指标i的5级标准值对“污染”的相对隶属度sih=1;介于1级与5级之间的指标i的h级标准值的指标标准相对隶属度sih用式(4)表示:

$s_{i h} = \frac{y_{i h} - y_{i 1}}{y_{i 5} - y_{i 1}}$ (4)

用上述两公式可以将矩阵Xm×n,Ym×5变为相应的实测指标相对隶属度矩阵和指标标准相对隶属度矩阵。

1.2 模糊模式识别评判模型的建立

隶属度可定义为权重,因此矩阵R5×4不仅描述了超标,且表示了权重,其按列归一化矩阵(超标权重矩阵):

其中,wij为超标权重, $w_{i j} = \frac{r_{i j} v_{i}}{\sum_{i = 1}^{m} r_{i j} v_{i}} ‚ v_{i}$ 为各指标的影响权重。

$\sum_{i = 1}^{m} w_{i j} = 1, \forall j$ (8)

井位j以相对隶属度uhj隶属于h级环境标准,则有井位相对隶属度矩阵:

满足的约束条件:

$\sum_{h = 1}^{5} u_{h j} = 1 ‚ \forall j$ ;0≤uhj≤1 (10)

则样本j与级别h间的差别用权距离向量来表示:

dj=(d1j,d2j,d3j,d4j,d5j)。

其中,

$d_{h j} = {\sum_{i = 1}^{m} (w_{i j} | r_{i j} - s_{i h} |)^{p}}^{p^{- 1}}$ (11)

其中,dhj为样本j与级别h间的广义权距离;p为距离参数,取欧拉距离2。

我们可建立样本j级别h对A的最优相对隶属度目标函数:

min ${F (u_{h j}) = \sum_{h = a_{j}}^{b_{j}} u_{h j}^{2} d_{h j}^{α}}$ (12)

其中,α为优化准则参数,取海明距离1。

由条件式(10)和式(12)来构造拉格朗日函数,求得样本j级别h对A的最优相对隶属度函数公式:

$u_{h j} = (d_{h j}^{α} \sum_{k = a_{j}}^{b_{j}} d_{k j}^{- α})^{- 1}, d_{h j} \neq 0, a_{j} \leq h \leq b_{j}$ (13)

样本j以相对隶属度uhj隶属于h级环境标准,则有最优相对隶属度矩阵:

其中,uij为井位最优相对隶属度;wij为超标权重;rij为实测指标相对隶属度;sih为指标标准相对度;h为级别代码;amin,amax分别与前aj,bj意义一致。

2 实例分析

咸阳市位于关中平原的中部,位于东经107°38′～109°10′,北纬34°11′～35°32′之间,南与西安市隔水相望,北与甘肃相连,东与渭南、铜川市毗邻,西与宝鸡市接壤。咸阳市是一个典型的依靠地下水进行工农业生产和居民生活的城市,地下水开采量占总用水量的80%～90%,近年来地下水受到了持续污染,水质恶化,对生产生活产生了较大威胁。因此本文以咸阳市为例进行分析(见表1)。

依据GB/T 14848-93Ⅲ地下水质量标准值,并结合咸阳市的实际情况,划分表1中五种指标的分类标准如表2所示。

则实测指标矩阵和指标标准矩阵分别为:

用上述公式可以将矩阵Xm×n,Ym×5变为相应的实测指标相对隶属度矩阵和指标标准相对隶属度矩阵:

计算超标权重矩阵,取vi=1/5得:

求得最优相对隶属度矩阵并考虑最大隶属度大于0.5的情况,得地下水水质的分类级别及变化趋势:

从计算结果来看,2号井的水质最差,处于5级;1号井从计算的相对隶属度来看,属于4级标准,但由于u51>0.5,即向5级水质过渡的趋势明显,因此,我们把它归为4级～5级;3号井,4号井的水质较好,达到2级标准;此外,从结果看,3号井和4号井虽同为2级水,但3号井的水质要优于4号井,这与实际情况中分析表现出来的结果相符。从整体上说,这四口井水质情况从好到坏依次是:3号井,4号井,1号井,2号井。通过与当地环保部门的分析资料对比,可以得出这个结果是可信的,是比较可观的。

3结语

通过以上的理论分析和实例验证,在目前评价地下水水质及其变化趋势的表现上,模糊模式识别法具有其独特的优点,主要表现在以下几方面:

1)该方法提出相对隶属度的概念,并通过加权距离平方和最小值来作为判定污染等级的依据。在多因子共同影响环境时,由式(14)可知,综合考虑了样本所处级别的所有信息,使结果更综合,更全面,也更客观,比以往的模糊综合评判法更能反映因子之间的作用。2)从计算结果分析对比来看,该法可以兼顾最大隶属度大于0.5的情况,并能得出水质的变化趋势,推溯污染的主要原因,这就可以为以后防止地下水进一步恶化提供依据,找到合理快捷的解决方案。3)该法的计算结果量化程度比较高,可以准确的提供有关评判参数,这就给对比分析提供了有力的武器和可信的证据,因此操作性较强;也可以编写相应的程序,实现在计算机上的自动化演算,从而提高水质分析预测的速度。

参考文献

[1]尹国勋,李振山.地下水污染与防治[M].北京:中国环境科学出版社,2005:44-51.

[2]陈守煜,赵瑛琪.模糊模式识别理论模型与水质评价[J].水利学报,1991(6):35-40.

[3]王国利,陈守煜,李成林.模糊模式识别模型在碧流河水库水质评价中的应用[J].大连理工大学学报,1997,11(6):700-703.

[4]郝聚民,林焰,纪卓尚,等.一种新的模糊模式识别理论模型及其在船舶配载评价上的应用研究[J].Ship Engineering,1998(3):20-22.

贵阳市地下水水质评价篇6

1 区内地质概况

1.1 气象、水文

黑河市属于中温带大陆性季风气候, 冬季严寒漫长, 春季风大雨少, 夏季湿热多雨, 秋季温低早霜。年降雨量在400-600mm之间, 主要集中在7-9月份, 年水面蒸发量在550-700mm之间, 陆地蒸发量在350-400mm之间。

市内有大小河流 60 条, 其中流域面积大于 200km2 以上的河流 26 条, 这些河流分属于黑龙江、嫩江水系。黑龙江在市区东侧由南向北流过, 市内流域面积 1.44 万 km2, 当地产水量 20.87×108km3, 多年平均进出水量 457×108km3。

1.2 地形地貌

黑河市区位于小兴安岭北部, 地形总的趋势是西北高, 东南部较低。海拨高度由500~800m递减到120m。依据地貌形态特征划分为低山、丘陵、冰蚀冰碛台地、二级阶地、一级阶地、漫滩及熔岩台地七个形态单元。

1.3 地质概况

黑河市区横跨罕达气优地槽褶皱带与小兴安岭-松嫩地块, 地层发育较齐全。上奥陶统至上志留统连续沉积了半深海、浅海、潮间带中酸性火山岩、陆源碎屑岩、碳酸盐地层。在此期间有深源超基性-基性-中基性岩浆岩沿断裂喷发及侵入。到中生界构造运动开始活跃及老的基底拉张断陷, 沉积了侏罗纪至上第三纪地层, 新生界地层主要分布在现代河谷及其两侧。

1.4 水文地质条件

根据地下水的赋存条件和水力性质, 将本区地下水划分为:松散岩类孔隙水、碎屑岩类裂隙、孔隙水和基岩裂隙水三大类型。

1.4.1 松散岩类孔隙水

主要蕴藏于第四系砂砾石层中。含水层厚度7~10m, 单井涌水量500~1000m3/d, 水位埋深一般在3~6m, 为区内主要供水目的层。主要以接受黑龙江的侧向补给为主, 大气降水的渗入补给为辅。另外, 该层与基岩裂隙水具有较好的水力联系, 可接受基岩裂隙水的侧向补给, 还可接受农田间灌溉水的回渗补给。河流对地下水主要起排泄作用, 人工取水是其主要的排泄途径之一。

1.4.2 碎屑岩类裂隙、孔隙水

1.4.2. 1 碎屑岩类孔隙潜水

蕴藏于第三系孙吴组砂岩、砂砾岩中。含水层厚1.5~3.5 m, 单井涌水量10~100m3/d, 无供水意义。主要以大气降水垂直渗入及基岩裂隙水的侧向补给, 以地下迳流形式排泄于下游。

1.4.2. 2 碎屑岩类裂隙、孔隙微承压水

分布于向斜盆地的砂岩、砂砾岩及含泥质砂砾岩中。含水层厚一般为20~42 m, 顶板埋深37~81 m, 为本区主要供水目的层之一。可接受上层水的越流补给, 同时还可接受盆地边缘基岩裂隙水的侧向补给, 以侧向迳入渗方式排泄于区外, 人工取水及自流也是其排泄途径之一。

1.4.3 基岩裂隙水

1.4.3. 1 风化带网状裂隙水

广泛分布于岩浆岩、变质岩及火山岩、火山碎屑岩等风化裂隙中, 风化带厚度15~55m, 一般厚度30~40m。

1.4.3. 2 孔洞裂隙水

分布于熔岩台地一带, 含水层厚一般为30~50 m, 单井涌水量100~1000m3/d, 地下水位埋深12~37m。

1.4.3. 3 构造裂隙水

分布于向斜盆地中, 岩性主要为凝灰质砂砾岩与砂质泥岩、凝灰岩互层, 夹较厚层玄武岩。含水带分布一般较稳定, 含水层厚一般为18~73 m, 地下水位埋深1~19m。单井涌水量一般100~1000m3/d。

风化带网状裂隙水及孔洞裂隙水均可接受大气降水的渗入补给, 次为地下水侧向迳流补给。多以较短的地下水迳流排泄于附近沟谷中或转化成地表水, 地表植被蒸腾作用也是其排泄途径之一, 还可通过泉水溢出形式进行排泄。

构造裂隙水主要接受风化网状裂隙水的补给, 只在局部地段含水层出露地带, 可直接接受大气降水的渗入补给, 通过泉或地下水溢出及人工取水进行排泄。

2 地下水水质现状

2.1 地下水水质特征

黑河市地下水一般多为无色、无味、无嗅、透明、无肉眼可见物, 化学类型主要为HCO3-Ca-Na型、HCO3-Ca型、HCO3-Na型水。PH值一般6.5~8, 铁含量一般为0~8.4mg/L, 锰含量一般为0~0.76mg/L。按照国家颁布的 (生活饮用水卫生标准) , 黑河市地下水做为饮用水, 其感官性指标, 一般化学指标, 毒理学指标及细菌学指标较为理想, 化学指标中仅亿别项目超过标准。其中75%的水点总铁含量超过0.3mg/L标准, 大部分地方氟含量偏低, 一般0.22mg/L, 个别村、屯偏高, 超过标准1.0mg/L。碎屑岩类裂隙孔隙水中含氟量较高, 部分符合饮用水标准, 多数超过饮用水标准。个别地方铅含量超过标准, 占分析水点近1%。黑龙江沿岸及黑河镇, 硝酸根离子及亚硝酸根离子含量偏高, 黑河镇最高达90mg/L, 氯和硫根离子含量增高, 地下水变为氯化物硫酸型水。主要原因是人为污染所致。总之, 区内地下水除铁、锰等指标超过饮用水标准外, 其它项目一般符合标准要求, 为比较适合饮用的地下水。

2.2 地下水污染程度评价

人为因素是本区浅层地下水污染的主要根源, 深层地下水有稳定隔水顶板, 渗透能力弱, 隔污能力强, 地下水没有受到污染。在工业及居民集中的城区内虽然粉质粘土层自净能力较强, 但因排污水沟渠开挖较深, 接近砂砾石层, 工业“三废”及生活污水、垃圾均都排放到洼地及沟渠中, 使污水经包气带下渗至地下, 污染地下水, 未经处理的固体废渣, 通过大气降水对固体物的淋漓, 使有害物质渗入地下, 污染地下水, 工业废气则随降雨渗入地下水中, 间接对地下水污染。工矿企业污水部分监测数据见表1。

近郊及远郊由于农田使用速效化肥, 对地下水造成一定的污染, 特别是近郊地区 (菜地) 更为突出。

3 结论

3.1 黑河市地下水类型属中性低矿化重碳酸型淡水。水化学类型主要为 HCO3- Ca- Na 型、HCO3- Ca 型、HCO3- Na 型水, 除铁、锰等指标超过饮用水标准外, 其它指标均符合标准要求, 为比较适合饮用的地下水。

3.2 在工业及居民集中的城区内, 局部由于受污染为氯化物硫酸型水, 不宜饮用。

4 建议

4.1 建立地下水动态监测网点, 实现水资源管理的信息化和科学化。

4.2 在工业及居民集中的城区内修建污水处理厂和垃圾处理厂, 完善配套排污管线。

4.3 对有能力产生污染的企业, 要配备污水净化设备, 严禁采用渗坑 (窖) 处理、排放污水。

贵阳市地下水水质评价篇7

1.1 模糊物元及复合模糊物元

给定事物的名称N, 它关于特征c有量值为v, 以有序三元R= (N, c, v) 组作为描述事物的基本元, 简称物元。其中, 量值v具有模糊性, 便称为模糊物元。如果事物N有n个特征c1, c2, …, cn和相应的模糊量值v1, v2, …, vn, 称R为n维模糊物元, 简记为R= (n, c, v) 。如果m个事物的n维物元组合在一起, 使构成m个事物n维复合物元Rmn。若将Rmn的量值改写为模糊物元量值, 称为m个事物n维复合模糊物元, 记作

$R_{m n} = [\begin{matrix} Μ_{1} & \dots & Μ_{m} \\ C_{1} & v_{11} & \dots & v_{m 1} \\ ⋮ & ⋮ & ⋮ \\ C_{n} & v_{1 n} & \dots & v_{m n} \end{matrix}] (1)$

式中:Rmn为m个事物的n个评价指标的复合物元;Mi为第i个事物 (i =1, 2, …, m) ;Ck为第k个特征 (k=1, 2, …, n) ;vik为第i个事物第k个特征对应的模糊量值。

1.2 从优隶属度原则

各单项评价指标相应的模糊量值, 从属于标准方案各对应评价指标相应的模糊量值隶属程度, 称为从优隶属度。由此建立的原则, 称为从优隶属原则。由于从优隶属度一般为正值, 可采用下面类型的指标:

越大越优型 μik=Vik/maxVik (2)

越小越优型 μik=minVik/Vik (3)

式中:μik为从优隶属度;Vik为第i个事物第k项特征对应的量值;maxVik、minVik分别为各事物中每一项特征所有量值Vik中的最大值和最小值。由此可以构建从优隶属度模糊物元 $\tilde{R}$ mn:

$\tilde{R}_{m n} = [\begin{matrix} Μ_{1} & \dots & Μ_{m} \\ C_{1} & μ_{11} & \dots & μ_{m 1} \\ ⋮ & ⋮ & ⋮ \\ C_{n} & μ_{1 n} & \dots & μ_{m n} \end{matrix}] (4)$

1.3 标准模糊物元与差平方复合模糊物元

标准模糊物元Ron是指从优隶属度模糊物元 $\tilde{R}$ mn中各评价指标的从优隶属度的最大值或最小值。本文以最大值作为最优, 也就是各指标从优隶属度均为1。

若以Δij (i =1, 2, …, n;j =1, 2, …, m) 表示标准模糊物元Ron与复合从优隶属度模糊物元R～mn中的各项差的平方, 则组成差平方复合模糊物元RΔ, 即Δji= (μ0i-μji) 2 (i =1, 2, …, n;j =1, 2, …, m) , 可表示为:

$R_{Δ} = [\begin{matrix} Μ_{1} & \dots & Μ_{m} \\ C_{1} & Δ_{11} & \dots & Δ_{m 1} \\ ⋮ & ⋮ & ⋮ \\ C_{n} & Δ_{1 n} & \dots & Δ_{m n} \end{matrix}] (5)$

2 熵值法确定权重系数

“熵权”理论是一种客观赋权方法, 熵是系统无序程度的度量, 它可以度量数据所提供的有效信息量, 因此, 可以用熵来确定权重。

信息论创始人申农 (C.E.Shannon) 将熵引人信息论, 用以表示系统的不确定性、稳定程度和信息量, 或者说是信息源信号状态多少——即其信息量的度量。信息熵的概念来源于信号通讯理论, 它是通过分析组成通讯信号的数字或符号的统计结构特征来定量表示信号通讯的能力——即信息量的大小。本文引用信息科学中信息熵的概念来确定水资源系统中各个评价指标的权重, 方法及求解步骤如下:

(1) 构建m个事物的n个评价指标的判断矩阵R= (vij) mn, (i =1, 2, …, n;j =1, 2, …, m) 。

(2) 将判断矩阵归一化处理, 得到归一化判断矩阵A:

对大者为优的指标而言, 有

$a_{i j} = \frac{v_{i j} - v_{\min}}{v_{\max} - v_{\min}} (6)$

对小者为优的指标而言, 有

$a_{i j} = \frac{v_{\max} - v_{i j}}{v_{\max} - v_{\min}} (7)$

式中:vmax、vmin分别为同指标下不同事物中最满意者或最不满意者 (越小越满意或越大越满意) 。

(3) 定义熵。在有n个评价指标, m个被评价对象的评估问题中, 可以确定第i个评价指标的熵为:

$\begin{array}{l} Η_{i} = - k \sum_{j = 1}^{m} f_{i j} l n f_{i j} (8) \\ f_{i j} = \frac{a_{i j}}{\sum_{j = 1}^{m} a_{i j}} (i = 1 ‚ 2 ‚ \dots ‚ n ； j = 1 ‚ 2 ‚ \dots ‚ m) \\ k = 1 / l n m \end{array}$

为使lnfij有意义, 当fij=0时, 根据水质评价的实际意义, 可以理解lnfij为一较大的数值, 与fij相乘趋于0, 故可认为fijlnfij=0.但当fij=1, fijlnfij也等于0, 这显然与熵所反映的信息无序化程度相悖, 不切合实际, 故需对fij加以修正, 将其定义为:

$f_{i j} = \frac{1 + a_{i j}}{\sum_{j = 1}^{m} (1 + a_{i j})} (9)$

(4) 定义熵权w。定义了第i个评价指标的熵之后, 可得到第i个评价指标的熵权定义, 即:

$w_{i} = \frac{1 - Η_{i}}{n - \sum_{i = 1}^{n} Η_{i}} (其中 0 \leq w_{i} \leq 1 ‚ \sum_{i = 1}^{n} w_{i} = 1) (10)$

3 贴近度和综合评价

考虑到本文具有综合评价的意义, 采用M (·, +) 算法, 先乘后加运算欧氏贴近度ρHj, 则

$ρ Η_{j} = 1 - \sqrt{\sum_{i = 1}^{n} w_{i} Δ_{j i}} (j = 1 ‚ 2 ‚ \dots ‚ m) (11)$

式中:ρHj为第n个方案与标准方案之间的相互接近程度, 其值越大表示两者越接近, 反之则相离较远, 以此来构造欧氏贴近度复合模糊物元RρH, 则

$R_{ρ Η} = [\begin{matrix} Μ_{1} & Μ_{2} & \dots & Μ_{m} \\ ρ Η_{j} & ρ Η_{1} & ρ Η_{2} & \dots & ρ Η_{m} \end{matrix}] (12)$

欧氏贴近度是表示各方案与标准方案 (最优方案) 之间的贴近程度, 可以根据欧氏贴近度的大小对各方案进行优劣排序, 并可进行分类。

4 实例应用

本文选取河南省叶县具有代表性的7个监测点进行水质基于熵权的模糊物元综合评价研究, 监测点分别为曹庄、李村、堰口、小河赵、姜庄、娘娘庙和刘店, 选取总硬度、锰、硝酸盐、亚硝酸盐、氨氮和溶解性总固体六项评价指标, 各指标在7个监测点的监测数据见表1。根据国家地下水质量标准 (GB/T14848-93) 的规定, 将水质级别划分为5级, 分级指标见表1。

4.1 建立评价模型

(1) 构建复合模糊物元。根据表1数据, 对7个监测点和分级标准建立12个样品6个指标的复合模糊物元。

(2) 根据式 (2) , b1～b6以越大越优原则构建从优隶属度模糊物元R～mn:

$\begin{array}{l} \tilde{R}_{m n} = \\ [\begin{matrix} 曹庄 & 李村 & 堰口 & 小河赵 & 姜庄 & 娘娘庙 & 刘店 & Ⅰ 级 & Ⅱ 级 & Ⅲ 级 & Ⅳ 级 & Ⅴ 级 \\ b_{1} & 0.687 3 & 0.756 4 & 0.612 7 & 0.703 6 & 0.765 5 & 0.598 2 & 0.600 0 & 0.272 7 & 0.545 5 & 0.818 2 & 1.000 0 & 1.000 0 \\ b_{2} & 0.080 0 & 0.090 0 & 0.120 0 & 0.020 0 & 0.020 0 & 0.140 0 & 0.110 0 & 0.050 0 & 0.050 0 & 0.100 0 & 1.000 0 & 1.000 0 \\ b_{3} & 0.090 5 & 0.359 4 & 1.000 0 & 0.678 9 & 0.309 5 & 0.007 5 & 0.066 4 & 0.033 3 & 0.083 2 & 0.332 8 & 0.499 2 & 0.4992 \\ b_{4} & 0.578 9 & 0.842 1 & 0.263 2 & 1.000 0 & 0.684 2 & 0.263 2 & 0.315 8 & 0.005 3 & 0.052 6 & 0.105 3 & 0.526 3 & 0.526 3 \\ b_{5} & 0.147 2 & 0.142 8 & 0.157 2 & 0.255 6 & 1.000 0 & 0.240 0 & 0.594 4 & 0.011 1 & 0.011 1 & 0.11 1 & 0.277 8 & 0.277 8 \\ b_{6} & 0.523 0 & 0.840 0 & 0.603 0 & 0.535 0 & 0.718 0 & 0.195 0 & 0.330 0 & 0.150 0 & 0.250 0 & 0.500 0 & 1.000 0 & 1.000 0 \end{matrix}] \end{array}$

(3) 根据标准模糊物元和 $\tilde{R}$ mn构建差平方模糊复合物元RΔ:

$\begin{array}{l} R_{Δ} = \\ [\begin{matrix} 曹庄 & 李村 & 堰口 & 小河赵 & 姜庄 & 娘娘庙 & 刘店 & Ⅰ 级 & Ⅱ 级 & Ⅲ 级 & Ⅳ 级 & Ⅴ 级 \\ b_{1} & 0.097 8 & 0.059 4 & 0.150 0 & 0.087 8 & 0.055 0 & 0.161 5 & 0.160 0 & 0.528 9 & 0.206 6 & 0.033 1 & 0.000 0 & 0.000 0 \\ b_{2} & 0.846 4 & 0.828 1 & 0.774 4 & 0.960 4 & 0.960 4 & 0.739 6 & 0.792 1 & 0.902 5 & 0.902 5 & 0.810 0 & 0.000 0 & 0.000 0 \\ b_{3} & 0.827 2 & 0.410 4 & 0.000 0 & 0.103 1 & 0.476 8 & 0.985 1 & 0.871 6 & 0.934 6 & 0.840 5 & 0.445 2 & 0.250 8 & 0.250 8 \\ b_{4} & 0.177 3 & 0.024 9 & 0.542 9 & 0.000 0 & 0.099 7 & 0.542 9 & 0.468 1 & 0.989 5 & 0.897 5 & 0.800 6 & 0.224 4 & 0.224 4 \\ b_{5} & 0.727 2 & 0.734 8 & 0.710 3 & 0.554 2 & 0.000 0 & 0.577 6 & 0.164 5 & 0.977 9 & 0.977 9 & 0.790 1 & 0.521 6 & 0.521 6 \\ b_{6} & 0.227 5 & 0.025 6 & 0.157 6 & 0.216 2 & 0.079 5 & 0.648 0 & 0.448 9 & 0.722 5 & 0.562 5 & 0.250 0 & 0.000 0 & 0.000 0 \end{matrix}] \end{array}$

(4) 用熵值法计算各指标权重。各指标的实际值按式 (6) 进行归一化处理得到判断矩阵Aij:

$A_{i j} = [\begin{matrix} 0.532 6 & 0.945 7 & 0.087 0 & 0.630 4 & 1.000 0 & 0.000 0 & 0.010 9 \\ 0.500 0 & 0.583 3 & 0.833 3 & 0.000 0 & 0.000 0 & 1.000 0 & 0.750 0 \\ 0.083 7 & 0.354 6 & 1.000 0 & 0.676 4 & 0.304 3 & 0.000 0 & 0.059 3 \\ 0.428 6 & 0.785 7 & 0.000 0 & 1.000 0 & 0.571 4 & 0.000 0 & 0.071 4 \\ 0.005 2 & 0.000 0 & 0.016 9 & 0.131 6 & 1.000 0 & 0.113 4 & 0.526 9 \\ 0.508 5 & 1.000 0 & 0.632 6 & 0.527 1 & 0.810 9 & 0.000 0 & 0.209 3 \end{matrix}]$

由式 (8) 、 (9) 计算各指标的熵Hi:

$Η_{i} = (0.980 6 ‚ 0.984 6 ‚ 0.984 5 ‚ 0.982 2 ‚ 0.982 0 ‚ 0.988 7)$

由式 (10) 计算可得各指标的权重wi:

$w_{i} = (0.199 3 ‚ 0.158 3 ‚ 0.159 0 ‚ 0.182 9 ‚ 0.184 9 ‚ 0.115 9)$

(5) 计算贴近度。由式 (11) 、 (12) 可以得到各样品的贴进度RpH:

$\begin{array}{l} R_{p Η} = [\begin{matrix} 曹庄 & 李村 & 堰口 \\ p Η j & 0.308 5 & 0.407 1 & 0.366 5 \end{matrix} \\ \begin{matrix} 小河赵 & 姜庄 & 娘娘庙 & 刘店 \\ 0.440 1 & 0.484 0 & 0.233 8 & 0.318 9 \end{matrix} \\ \begin{matrix} Ⅰ 级 & Ⅱ 级 & Ⅲ 级 & Ⅳ 级 & Ⅴ 级 \\ 0.082 2 & 0.146 9 & 0.274 0 & 0.578 9 & 0.578 9 \end{matrix}] \end{array}$

4.2 结果分析

经过计算, 得出各水质监测点的评价结果。另外采用了综合评价法对各监测点进行评价, 在确定各监测点的最终水质级别时, 按照置信度准则进行评价 (置信度λ 取值在0.7左右) 。两种方法评价结果对比见表2所示。

评价结果表明, 各取水点的水质级别基本为Ⅲ级, 满足国家标准规定的生活饮用水水质标准。但也有一些取水点污染相当严重, 如小河赵紧靠县城, 西边有盐厂, 北边有蓝光电厂等, 附近污水排放比较严重, 地下水水质受到很大程度的污染;姜庄则紧挨一些养猪场、养鸡场, 附近污水排放也是相当严重, 地下水水质受到不同程度的污染。而其他地区也有不同程度的污染但污染程度相对较轻。

与模糊综合评价方法进行对比 (表2) , 结果表明, 模糊物元评价方法和模糊综合评价方法得到的结果基本一致, 但是结合当地的实际情况模糊物元比模糊综合评价更准确一些。如曹庄、娘娘庙附近有一些养猪场、养鸡场, 地下水水质受到了不同程度的污染, 更符合Ⅲ级标准, 模糊综合评价结果存在一些误差。同时表明了基于熵权的模糊物元地下水水质评价模型应用在地下水水质评价中是合理可行的, 且计算简便实用。

5 结语

(1) 本文将欧氏贴近度与模糊物元结合起来, 建立了基于欧氏贴近度的模糊物元模型, 该模型以解决不相容问题为核心, 适用于多因子评价问题。并将该模型应用到地下水水质综合评价中, 进行评价时将分类等级作为物元的事物来构造复合模糊物元矩阵, 通过计算与理想物元的贴近程度来实现水样水质的等级评价。

(2) 引入熵值理论, 从数据本身所反映的信息无序化效用值来计算权重系数, 避免了专家赋权的偏好性, 能减少人为的主观因素对评价结果的影响, 可提高水质评价模型结果的精度。

(3) 实际计算结果表明, 该模型理论简单, 使用简便, 解决了地下水水质评价问题中的模糊性和不相容性, 评价结果更准确、合理。

(4) 与模糊综合评价结果相比, 模糊物元评价结果更加符合实际情况, 表明该模型用于解决该问题是可行的。同时该模型还可应用于城市环境质量评价、区域水资源承载力评价等领域, 具有广泛的推广前景。

参考文献

[1]谢季坚, 刘承平.模糊数学方法及其应用[M].武汉:华中科技大学出版社, 2004.

[2]蔡文.物元模型及应用[M].北京:科学技术文献出版社, 1994.

[3]张斌, 雍歧东, 肖芳淳.模糊物元分析[M].北京:石油工业出版社, 1997.

[4]闫文周, 顾连胜.熵权决策法在工程评价中的应用[J].西安建筑科学大学学报, 2004, 36 (1) :98-100.

[5]肖芳淳.输气管道土壤腐蚀性模糊物元贴近度聚类分析[J].油气储运, 1999, 18 (6) :35-37.

[6]邹同庆.地下水水质的物元评价方法及应用[J].常德师范学院学报 (自然科学版) , 2001, 13 (4) :31-34.

[7]杨立君.用模糊数学方法综合评价地下水水质[J].矿山设计研究, 1988 (3) :28-33.

贵阳市地下水水质评价篇8

关键词：地下水质量评价,单因子评价法,综合评价法,影响因素

1 研究区概况

凤阳县XX矿区位于江淮丘陵区北缘, 地形为丘陵岗地, 属北亚热带季风气候, 年平均气温15.01℃, 年平均降水量923.53 mm, 年平均蒸发量1559.82 mm。本区出露地层属华北地层大区晋冀鲁豫地层区徐淮地层分区, 主要有新太古界庄子里组、峰山李组。矿区大地构造位置处于中朝准地台华北陆块徐淮地块蚌埠隆起, 区内断层有近东西向及北北东向两组。本次工作通过对矿区矿坑道水、5个抽水孔及周边民井等共计采水样15个进而对地下水水质进行分析。

2 水文地质条件

2.1 含水岩组

(1) 松散岩类孔隙弱含水岩组:岩性主要为粉质粘土、粘土, 分布于矿区的北部、东部及西南部, 厚0~21.79 m。就整个矿区看, 由北东向西南厚度有所增大。单井涌水量<10 m3/d, 水位埋深0.84~8.40 m。地下水水质较差, 水化学类型以HCO3-Ca、Na为主。 (2) 裂隙弱含水岩组:岩性主要为黑云角闪斜长片麻岩、角闪斜长片麻岩、斜长片麻岩、闪长岩及大理岩、榴辉岩等, 分布于矿区的中部。该含水岩组埋深0~21.79 m, 厚度大于300 m。单井涌水量<50 m3/d, 其顶部0.81~21.65 m有一风化带, 岩性多为闪长岩, 岩石风化强烈, 多呈砂粒状、碎块状。

2.2 地下水补给、迳流、排泄

松散岩类孔隙含水岩组主要接受大气降水垂直入渗补给, 蒸发为其主要排泄途径、并向下补给基岩裂隙含水岩组。

基岩裂隙含水岩组主要接受上覆松散岩类孔隙含水岩组的补给, 地下水由西南向东北缓慢迳流, 排泄以侧向迳流和人工开采为主。

3 地下水质情况分析评价

在本区内共采集了水样15个进行化学污染成份分析, 测试内容主要是水的p H、总硬度、溶解性总固体、SO42-、Cl-、NH4+、NO3-、NO2-、F-、Cu、Pb、Cd、Hg、As、Zn、Ni、Cr等组分。测试结果中, 以阴离子划分, HCO3型水占总水样60%, HCO3·Cl型水占总26%, Cl·HCO3型水占总7%, HCO3·SO4型水占总7%, 依据地下水质量标准 (GB/T14848-93) , 对矿区地下水质量和生活饮用水水质分析评价, 其中5个抽水孔及矿坑道水样代表矿区基岩地下水, 矿山生活用水及8个民井水样代表浅层地下水。

3.1 生活饮用水水质分析评价

采用单因子评价方法:

(1) 对于评价标准为定值的水质因子, 其标准指数计算公式:

式中:Pi—第i个水质因子的标准指数, 无量纲;

Ci—第i个水质因子的监测浓度值, mg/L;

Csi—第i个水质因子的标准浓度值, m g/L。

(2) 对于评价标准为区间值的水质因子 (如ph值) , 其标准指数计算公式:

式中:Pph—p H的标准指数, 无量纲;

PH—p H的检测值;

PHsu—标准中p H的上限值;

PHsd—标准中p H的下限值。

根据测试结果, 矿区钻孔、矿坑道水、天津铺、杜涧子、张家洼及小王府民井地下水均符合生活饮用水国家卫生标准, 而矿山生活用水井、大王府、许港家、栾家及后付家民井地下水各有不同指数超标, 其中矿山生活用水铁含量0.40 mg/L, 超标0.33倍, 锰含量0.18 mg/L, 超标0.8倍;大王府总硬度587.03 mg/L, 超标0.3倍, 溶解性总固体1156.94 mg/L, 超标0.16倍;许港家溶解性总固体1177.23 mg/L, 超标0.18倍;栾家亚硝酸盐含量0.14 mg/L, 超标1.12倍;后付家总硬度539.99 mg/L, 超标0.2倍, 氯化物含量286.79 mg/L, 超标0.15倍。

3.2 地下水质量的分析评价

4 水质状况影响因素分析

通过评价对比结果, 可知Z K 3 0 4、ZK1108、ZK908、ZK1320钻孔、矿坑道水, 天津铺、杜涧子及小王府民井符合生活饮用水国家卫生标准且水质良好;ZK920钻孔及张家洼民井虽符合生活饮用水标准但水质较差;矿山生活用水井、大王府、许港家、栾家及后付家民井地下水各有不同指数超标且水质较差, 其中矿山生活用水铁、锰超标, 主要由于Qp3地层富含铁锰结核, 经过水、岩 (土) 长期的溶滤、交代作用, 地层中的铁锰则以离子形式进入地下水中, 导致含量超标;大王府总硬度及溶解性总固体超标, 许港家溶解性总固体超标, 而大王府水样总硬度几乎全为暂时硬度, 故不进行分析, 两家民井溶解性总固体超标, 主要因其地形为丘陵岗地, 地势高差不大, 水交替强度弱, 且以蒸发排泄为主, 从而引起浅层地下水矿化度升高;栾家亚硝酸盐超标, 经调查其主要原因是水井成井工艺较简单, 无过滤砾层和止水工序, 加之当地居民环保意识差, 生活垃圾随意堆放造成环境污染, 地下水交替时亚硝酸盐进入地下水中, 导致亚硝酸盐超标;后付家主要是总硬度及氯化物超标, 主要原因为地下径流缓慢, 又以蒸发排泄为主, 致使地下水中氯及钙镁浓缩导致超标。

5 结语

通过对凤阳县XX矿区地下水水质情况的分析, 得出该区基岩地下水质属于良好, 适用于生活饮用及工、农业用水, 但矿区浅层地下水水质较差。建议采矿的同时加强环境保护意识, 采取必要的环保处理措施, 避免水土污染对矿区周围居民身体健康造成危害。

参考文献

[1]《安徽省凤阳县江山矿区铅锌金矿勘探地质报告》.安徽省地质矿产勘查局三一二地质队, 2014.

【贵阳市地下水水质评价】推荐阅读：

贵阳市建设10-09

贵阳市绿化带12-07