基于云模型的東平湖枯水期地下水水質(zhì)評價

胡尊芳,孫建峰,宋印勝,王仕昌,韓玉珍

(1.山東省物化探勘查院,山東 濟(jì)南　250013; 2.山東大學(xué)土建與水利學(xué)院,山東 濟(jì)南　250061)

?

基于云模型的東平湖枯水期地下水水質(zhì)評價

胡尊芳1,2,孫建峰1,宋印勝1,王仕昌1,韓玉珍1

(1.山東省物化探勘查院,山東 濟(jì)南250013; 2.山東大學(xué)土建與水利學(xué)院,山東 濟(jì)南250061)

水質(zhì)評價；南水北調(diào)東線工程;東平湖;枯水期;地下水;云模型

東平湖位于山東省泰安市,是南水北調(diào)東線工程的調(diào)蓄水庫,同時也擔(dān)負(fù)著黃河下游地區(qū)和山東省西水東調(diào)水資源調(diào)蓄調(diào)配的功能,地位重要。地下水是東平湖水環(huán)境的重要組成部分,一方面地下水是當(dāng)?shù)鼐用裰饕纳钏?另一方面地下水通過與地表水的相互交換影響地表水水質(zhì),進(jìn)而影響南水北調(diào)工程受水區(qū)域的用水安全,因此有必要對東平湖地下水水質(zhì)進(jìn)行調(diào)查評價。

有關(guān)地下水水質(zhì)評價的方法有:單因子評價法、綜合指數(shù)法、人工神經(jīng)網(wǎng)絡(luò)模型、模糊綜合評價法、灰色聚類法等[1],但每種方法都有其局限性,鑒于此,筆者引入云模型地下水水質(zhì)評價方法。云模型是基于概率論和模糊數(shù)學(xué)理論而構(gòu)建出來的一種定性概念與其定性數(shù)值表示之間的不確定性轉(zhuǎn)換模型[2]。云模型能綜合考慮評價過程中評價標(biāo)準(zhǔn)的非線性、評價因子的多樣性以及評價過程之中的模糊性與隨機(jī)性等特點,在空氣質(zhì)量評價、水資源再生能力評價、土地利用風(fēng)險評價、地震風(fēng)險評價等多方面得到廣泛應(yīng)用[3-7]。

云模型在水環(huán)境評價領(lǐng)域也有應(yīng)用,如楊文[8]首次采用云模型與Spearman秩相關(guān)系數(shù)相結(jié)合的方法對寧?？h地表水環(huán)境質(zhì)量進(jìn)行了評價,結(jié)果表明模型分析結(jié)果與實際情況基本符合,云模型可以很好彌補(bǔ)Spearman秩相關(guān)系數(shù)中存在的不均勻性量化值缺陷;丁昊等[9]利用云模型對我國14個代表性湖庫的富營養(yǎng)化情況進(jìn)行了評價,并與國內(nèi)外廣泛應(yīng)用的營養(yǎng)狀態(tài)指數(shù)法、評分公式法及模糊綜合評價法所得的結(jié)果進(jìn)行了對比,驗證了云模型方法的可行性和有效性;白曉燕等[10]采用基于Critic權(quán)重的云模型對東江干流水質(zhì)進(jìn)行了評價,結(jié)果表明云模型方法合理可行,且更為直觀地反映了水質(zhì)評價中的模糊性和隨機(jī)性;王昭洋等[11]采用改進(jìn)賦權(quán)的多維云模型對國內(nèi)12個水庫的富營養(yǎng)化進(jìn)行了評價,準(zhǔn)確獲得了湖庫水質(zhì)的富營養(yǎng)化級別。

目前,東平湖地下水水質(zhì)評價工作開展的相對較少,在2013年10月19日南水北調(diào)東線一期工程全線通水后,相關(guān)的工作更是沒有展開。本文采用云模型對南水北調(diào)東線工程通水前后的東平湖地下水水質(zhì)進(jìn)行評價,以期為南水北調(diào)東線工程的運行管理提供參考。

1　材料與方法

1.1樣品采集與分析

圖1　東平湖水系圖及采樣點位置

1.2云模型評價

表1　評價因子的分級標(biāo)準(zhǔn)　mg/L

鑒于各評價因子對水質(zhì)評價結(jié)果的貢獻(xiàn)不同,權(quán)重也不相同,筆者根據(jù)評價因子的實測質(zhì)量濃度和各級別標(biāo)準(zhǔn)濃度計算權(quán)重。設(shè)評價因子合集X={x1,x2,x3,x4,x5,x6},它的權(quán)重矩陣為A=[a1,a2,a3,a4,a5,a6],權(quán)重的計算公式為

(1)

(2)

b. 特征參數(shù)確定。在確定評價因子和評價標(biāo)準(zhǔn)的基礎(chǔ)上,根據(jù)文獻(xiàn)提出的方法確定云模型的3個數(shù)字特征:期望Ex,熵En,超熵He[12]。假設(shè)對于某變量VQa,具有上下邊界[Bmin,Bmax],則變量VQa的云模型的3個數(shù)字特征可以通過下式計算:

(3)

式中:Ex為云滴在論域空間分布的期望,最能代表定性概念的中心點,是概念量化的最典型樣本點;En為定性概念的不確定性度量,一方面是定性概念重合度量,反映論域空間可被定性概念接受的云滴取值范圍,另一方面又是定性概念隨機(jī)性度量,反映了能夠代表這個定性概念云滴的離散度;He為熵的不確定性度量,反映了在論域空間代表該語言值的所有點的不確定度的凝聚性,其取值范圍間接反映了云滴的厚度,由熵的隨機(jī)性和模糊性共同確定;Bmin、Bmax分別為變量VQa的最小和最大邊界;k為常數(shù),本文中k=0.01。得到不同等級云模型數(shù)字特征計算方法見表2。

表2中a～d分別代表評價因子xi的分級標(biāo)準(zhǔn)中的各個邊界值,以總硬度為例,其各邊界值分別為:a=150,b=300,c=450,d=550。

圖2　6個評價因子對應(yīng)各水質(zhì)級別的云模型

c. 云模型的生成。根據(jù)確定的云模型數(shù)字特征,運用matlab R2007b編寫正態(tài)云發(fā)生器及半云發(fā)生器程序,輸入相應(yīng)的云參數(shù),分別生成6個評價因子的云模型,如圖2所示。橫坐標(biāo)代表評價因子的質(zhì)量濃度,縱坐標(biāo)代表評價因子的某一實測值通

表2　不同等級云模型數(shù)字特征計算方法

過云模型映射所對應(yīng)的確定度U。橫坐標(biāo)從左到右分別代表地下水質(zhì)量標(biāo)準(zhǔn)的Ⅰ級到Ⅴ級所對應(yīng)的云。以圖2(a)總硬度為例,某實測點的總硬度為280 mg/L時,由正向云發(fā)生器可以得到總硬度隸屬于各地下水水質(zhì)級別的確定度分別為:UⅠ(總硬度)=0,UⅡ(總硬度)=0.545 8,UⅢ(總硬度)=0.164 6,UⅣ(總硬度)=0,UⅤ(總硬度)=0,其中UⅡ值最大,UⅢ次之,其他確定度為0,說明總硬度為280 mg/L時,隸屬于水質(zhì)級別Ⅱ類的可能性最大,隸屬于Ⅲ類有一定的可能,但是不隸屬于其他的水質(zhì)級別。同理,也可以求出其他評價因子的確定度,得到相對應(yīng)的確定度矩陣。

d. 評價綜合確定度的計算。在得到各評價因子隸屬于各水質(zhì)級別的確定度矩陣[X][Y]后,乘以相應(yīng)的權(quán)重A,得到各評價因子的綜合確定度A[X][Y],計算其各列值的綜合確定度,最大確定度所在的級別即為監(jiān)測點所屬的水質(zhì)級別。

2　結(jié)果與討論

南水北調(diào)東線工程通水前后東平湖地下水水質(zhì)評價結(jié)果見表3。從表3可以看出,南水北調(diào)東線工程通水前,東平湖地下水水質(zhì)級別為Ⅰ～Ⅴ類均有分布,而在南水北調(diào)東線工程通水后,東平湖地下水水質(zhì)為Ⅰ類或Ⅴ類,其中S1、S3、S5點的水質(zhì)級別在通水前后沒有差別,分別為Ⅰ、Ⅴ、Ⅴ類,S2、S6、S7點在通水后的水質(zhì)級別優(yōu)于通水前,分別由Ⅱ、Ⅲ、Ⅴ類提升至Ⅰ類,而S4點的水質(zhì)級別在通水后下降一個級別,由Ⅳ類下降到Ⅴ類。另外還發(fā)現(xiàn)無論通水前后,東平湖東部地下水水質(zhì)均劣于西部地下水水質(zhì)。Wang等[13]的研究表明，東平湖東部地區(qū)的底泥重金屬污染更加嚴(yán)重,說明東部地區(qū)存在污染源,推測可能是工礦業(yè)污染或農(nóng)業(yè)污染。

表3　通水前后各采樣點水質(zhì)級別確定度及水質(zhì)評價結(jié)果

另外,即使通水前后的水質(zhì)級別一致,其確定度的不同又代表了不同意義,如S1點在通水前后的地下水水質(zhì)級別均為Ⅰ類,但通水后S1點的UⅠ更大,故通水后S1點的水質(zhì)級別更傾向于Ⅰ類,水質(zhì)更好;對于S3點來說,通水前后水質(zhì)級別均為Ⅴ類,且兩者的確定度相差不大,因此可以認(rèn)為通水前后S3點的地下水水質(zhì)沒有明顯變化;對于S5點來說,通水前后的地下水水質(zhì)級別均為Ⅴ類,但通水后S5點的UⅤ更大,更傾向于Ⅴ類,水質(zhì)也更差。

3　結(jié)　論

a. 云模型評價方法可以克服單一評價指標(biāo)的隨機(jī)性和模糊性,能夠客觀合理地評價東平湖地下水水質(zhì)。

b. 南水北調(diào)東線工程通水后,在7個采樣點所采集的地下水水樣的水質(zhì)級別并未變差,原因可能是因為南水北調(diào)東線工程的引水為長江地表水,水中重金屬質(zhì)量濃度較低,因此客水的引入降低了東平湖地下水中重金屬質(zhì)量濃度。

c. 東平湖東部地區(qū)的地下水水質(zhì)劣于西部地區(qū),是因為東部地區(qū)存在污染源,推測可能是工礦業(yè)污染或農(nóng)業(yè)污染。

在后續(xù)的研究工作中,建議增加地下水采樣點數(shù)目并對客水水質(zhì)進(jìn)行分析,以期進(jìn)一步了解南水北調(diào)東線工程對東平湖地下水水質(zhì)的影響。

[ 1 ] 王文強(qiáng).地下水水質(zhì)評價方法淺析[J].地下水,2007,29(6): 37-39.(WANG Wenqiang.A brief analysis on groundwater quality evaluation methods[J].Ground Water,2007,29(6):37-39.(in Chinese))

[ 2 ] 李德毅,杜鹢.不確定性人工智能[M].北京: 國防工業(yè)出版社,2005.

[ 3 ] 張峰,張鵬林,呂志勇,等.云模型在城鎮(zhèn)空氣質(zhì)量評價中的應(yīng)用[J].環(huán)境科學(xué)與技術(shù),2009,32(6):160-164.(ZHANG Feng,ZHANG Penglin,LYU Zhiyong,et al.Assessment of urban air quality based on cloud models[J].Environmental Science & Technology,2009,32(6):160-164.(in Chinese))

[ 4 ] 賈琦,段春青,陳曉楠.黃河流域水資源可再生能力評價的云模型[J].中國人口·資源與環(huán)境,2010,20(9): 48-52.(JIA Qi,DUAN Chunqing,CHEN Xiaonan.Cloud model of water resources reproducible ability in the Yellow River Basin[J].China Population,Resources and Environment,2010,20(9):48-52.(in Chinese))

[ 5 ] 周啟剛,張曉媛,王兆林.基于正態(tài)云模型的三峽庫區(qū)土地利用生態(tài)風(fēng)險評價[J].農(nóng)業(yè)工程學(xué)報,2014,30(23): 289-297.(ZHOU Qigang,ZHANG Xiaoyuan,WANG Zhaolin.Land use ecological risk evaluation in Three Gorges Reservoir area based on normal cloud model[J].Transactions of the Chinese Society of Agricultural Engineering,2014,30(23):289-297.(in Chinese))

[ 6 ] 張秋文,章永志,鐘鳴.基于云模型的水庫誘發(fā)地震風(fēng)險多級模糊綜合評價[J].水利學(xué)報,2014,45(1): 87-95.(ZHANG Qiuwen,ZHANG Yongzhi,ZHONG Ming.A cloud model based approach for multi-hierarchy fuzzy comprehensive evaluation of reservoir-induced seismic risk[J].Journal of Hydranlic Engineering,2014,45(1):87-95.(in Chinese))

[ 7 ]王賀，劉高峰，王慧敏.基于云模型的城市極端雨洪災(zāi)害風(fēng)險評價[J] .水利經(jīng)濟(jì)，2014，32(2)：15-18.(WANG He,LIU Gaofeng,WANG Huimin.Urban extreme rainstorm flood loss based on cloud model[J] .Journal of Economics of Water Resources，2014，32(2)：15-18.(in Chinese))

[ 8 ] 楊文.基于云模型的寧?？h地表水環(huán)境質(zhì)量評價[D].金華: 浙江師范大學(xué),2013.

[ 9 ] 丁昊,王棟.基于云模型的水體富營養(yǎng)化程度評價方法[J].環(huán)境科學(xué)學(xué)報,2013,33(1): 251-257.(DING Hao,WANG Dong.The evaluation method of water eutrophication based on cloud model[J].Acta Scientiae Circumstantiae,2013,33(1):251-257.(in Chinese))

[10] 白曉燕,賴成光,陳曉宏,等.基于Critic-cloud模型的東江干流水質(zhì)評價[J].水資源保護(hù),2014,30(5): 26-31,59.(BAI Xiaoyan,LAI Chengguang,CHEN Xiaohong,et al.Evaluation of water quality of trunk stream of Dongjiang River based on Critic-cloud model[J].Water Resources Protection,2014,30(5):26-31,59.(in Chinese))

[11] 王昭洋,王立,王小藝,等.改進(jìn)賦權(quán)的水體富營養(yǎng)化多維云模型評價[J].計算機(jī)與應(yīng)用化學(xué),2015,32(3): 346-350.(WANG Zhaoyang,WANG Li,WANG Xiaoyi,et al.Multi-dimensional cloud model evaluation of water eutrophication based on improved weight method[J].Computers and Applied Chemistry,2015,32(3):346-350.(in Chinese))

[12] 丁昊.云模型在水體富營養(yǎng)化評價及水文序列展延中的應(yīng)用[D].南京: 南京大學(xué),2013.

[13] WANG Yunqiang,YANG Liyuan,KONG Linghao,et al.Spatial distribution,ecological risk assessment and source identification for heavy metals in surface sediments from Dongping Lake,Shandong,East China[J].Catena,2015,125: 200-205.

Evaluation of groundwater quality of Dongping Lake in dry season based on cloud model

HU Zunfang1, 2, SUN Jianfeng1, SONG Yinsheng1, WANG Shichang1, HAN Yuzhen1

(1.ShandongInstituteofGeophysicalandGeochemicalExploration,Jinan250013,China;2.SchoolofCivilEngineering,ShandongUniversity,Jinan250061,China)

water quality evaluation; South-to-North Water Diversion Eastern Route Project; Dongping Lake; dry season; groundwater; cloud model

10.3880/j.issn.1004-6933.2016.05.015

山東省物化探勘查院科研基金(2015)；山東省地質(zhì)勘查項目(魯勘字(2011)121);山東省地礦局地質(zhì)科技攻關(guān)項目(KY201516)

胡尊芳(1986—)，女,博士,主要從事污水處理與水環(huán)境評價研究。E-mail:shiwaixianshu@126.com

X824

A

1004-6933(2016)05-0074-05

2015-10-22編輯：王芳)