灰色聚类在邢台山前倾斜平原区孔隙水脆弱性评价中的应用
灰色聚类在邢台山前倾斜平原区
孔隙水脆弱性评价中的应用
王占辉~姜正义~徐素娟~谭俊
(河北省地矿局第十一地质大队水环处,河北 邢台 054000) 摘 要:以DRASTIC模型为基础~采用层次分析法依据当地的实际情况确定权重~结合MAPGIS的叠加技术~形成了评价地下水污染难易程度的新体系~并将其应用于邢台市山前倾斜平原区孔隙水脆弱性评价。评价结果从趋势上反映了评价区域内孔隙水受污染的难易程度~与实际情况吻合较好。
关键词:地下水 脆弱性 DRASTIC模型 MAPGIS 0 引言
地下水脆弱性是1968年Margat提出的。他初步把地下水脆弱性定义为污染源从地
表
关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf
渗透与扩散到地下水面的可能性,并且通过图件来描述地下水对污染的脆弱程度。此后,随着对地下水脆弱性复杂影响因素的深入认识和研究水平的逐渐提高,地下水脆弱性概念也在不断的丰富、完善和发展。80年代中期以前,对地下水脆弱性的定义主要是从含水层的地质、水文地质特性等内部因素角度考虑,把地下水的脆弱性理解为含水层的一种内在的自然属性。也就是把地下水脆弱性单纯地理解为固有脆弱性。而80年代末以后,对地下水脆弱性的定义方式有所突破,除考虑地下水系统内部因素外,还考虑了土地利用和污染源等外部人为因素对脆弱性的影响。
1灰色聚类方法简介
灰色聚类是以灰数的白化权函数生成为基础,将一些观测指标或观测对象聚类成若干个可定义的类别的方法。
1.1灰色聚类样本的构成
ikk记=1,2,…,n为聚类对象;j=1,2,…,m为聚类指标;=1,2,…,为聚类灰数,即灰类。d为第i个聚类对象对于第j个聚类指标的样本值,D是以d为元素的样本ijij矩阵。
dd?d,,11121m,,dd?d21222m,,Dd ,,()ij,,????
,,dd?dn1n2nm,,
为便于对各样本指标进行综合分析并使聚类结果具有可比性,须对环境质量聚类样本各个指标的白化数和灰类进行
标准
excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载
化处理。数据标准化处理就是对水质聚类样本各个指标的白化值和灰类进行无量纲化处理,以便于对各样本指标进行 综合分析并使聚类结果具有可比性。常用的方法为污染指数法,这样即可达到无量纲化的目的,又可反映污染物的相对污染状况。
1.2白化函数的构成
jk白化函数反应了聚类指标对灰类的亲疏关系。第个指标的灰类1、灰类(k=2,3,…,h-1)和灰类h的白花函数分别为:
1 , ;xx,,m, fxxxxxxxx()()/() , ,,,,,;,1jhhmmh
,0 , ;xx,h,
0 , ;xx,,0,()/() , ;xxxxxxx,,,,000mm,,fxxxxxxxx()()/() , ;,,,,, ,jkhhmmh
, 1 , ;xx,m,
,0 , ;xx,h,
1 , ;xx,,m, fxxxxxxxx()()/() , ;,,,,,,000jhmm
,0 , ;xx,0,
f设为第j个聚类指标属于k灰类的白化函数,则白化函数有三种形式,其jk
,f中为白化函数的阈值 jkjk
1
fffj1j2j3
11
00λλj2j3ddd
0λj1
1.3聚类权与聚类系数的计算
聚类权是衡量各指标对同一灰类的权重。在孔隙水含水层脆弱性分析中,由于聚类指标的意义、量纲相同,绝对值相差不大,因此需按式(1-2)进行处理,即按
,jk (1,1) η,jkm
,,jkj,1
若聚类指标的意义、量纲不同,且不同指标的样本值在数量上相关很大时,不能
η按上式计算,否则将会引起评估偏差。 jk
,ik聚类系数反映了聚类对象对灰类的隶属程度。 ik
n
,,f(d),,(i,1,2,?,n;k,1,2,3),ikjkijjk,j1
jf(d)kd式中为指标,灰类的白化权函数,可由样本值代入相应的白化jkijij
权函数计算式中求得:
σ,max{σ}ikik*1,k,5ik若,则称聚类对象属于灰类。从各行找出最大聚类系数,
i该最大聚类系数所属灰类即为聚类对象(各测点)所属灰类(级别)。 2实例分析
2.1含水层环境背景
邢台市山前倾斜平原主要分布在东部,其孔隙水含水层水平方向岩性由西向东,由卵石、砾石逐渐变为粗砂、中砂,局部还有中细砂、细粉砂。其含水层由厚变薄,层数变多。富水性也由冲洪积扇轴部向其边缘减弱的规律。而垂直方向岩性由上到下,全新统较薄,一般透水不含水,上更新统为本区的主要含水层(局部疏干),又称为浅层水,岩性以底部卵砾石、粗砂为主,局部为中细砂。中更新统、下更新统岩性以中粗砂为主,又称为深层水,局部为中细砂,且砂多呈风化状、含粘土质较多,富水性较差。
2.2评价指标确定
结合邢台市孔隙水含水层具体情况,本次选取七个代表性指标:地下水埋深、含水层岩性、降水入渗补给量、土壤类型、地形坡度、含水层渗透系数和包气带
介质类型,由于评价因子指标意义、量纲不同,利用进行评分。
孔隙水脆弱性分类指标表
序评价因子 I类 II类 III类 IV类 V类 号
地下水埋深(m) ?0 10,15 15,20 20,25 >25 1 评分 6 4 3 2 1
以卵砾石以粗砂含砾以粗砂为以中砂为以细砂为含水层岩性 为主 为主 主 主 主 2
评分 10 8 6 5 4 降水入渗补给量>175 150,175 125,150 100,125 ?100 ,mm/a, 3
评分 7 6 5 4 3
冰积、残破土壤类型 裸露、洪积 冲积 冲洪积 风积 积 4
评分 10 9 6 3 1 地形坡度,?, 0,2 2,3 3,6 6,9 9,12
5
评分 10 8 7 5 3 含水层渗透系数>400 300,400 200,300 100,200 ?100 ,m/d, 6
评分 10 8 6 4 2 包气带介质类型 卵砾石 砂类 粉土 粉质粘土 粘土 7 评分 10 7 4 2 1 灰色聚类权η的计算结果 jk
灰类
I类 II类 III类 IV类 V类 指标
0.375 0.25 0.187 0.125 0.0625 1
0.242 0.18 0.15 0.12 0.3032
0.24 0.2 0.16 0.12 0.283
0.31 0.20 0.10 0.034 0.3454
0.242 0.21 0.15 0.09 0.3035
0.26 0.2 0.13 0.066 0.3336
0.29 0.16 0.083 0.041 0.4167