劉海嬌，仕玉治，范明元，張 欣

(1.山東省水利科學(xué)研究院，山東 濟(jì)南 250013； 2.山東省水資源與水環(huán)境重點(diǎn)實(shí)驗(yàn)室，山東 濟(jì)南 250013)

黃河三角洲是中國(guó)暖溫帶最年輕、經(jīng)濟(jì)增長(zhǎng)最快、景觀變化最劇烈的地區(qū)，也是東北亞內(nèi)陸和西太平洋鳥類遷徙重要的中轉(zhuǎn)站、越冬棲息地和繁殖地。該區(qū)域是國(guó)家、“十一五”“十二五”重點(diǎn)建設(shè)區(qū)域，同時(shí)也是山東省區(qū)域經(jīng)濟(jì)發(fā)展“一體兩翼”中北翼主體，戰(zhàn)略地位十分重要。當(dāng)?shù)亟?jīng)濟(jì)社會(huì)的快速發(fā)展給生態(tài)環(huán)境及水資源造成了很大影響，尤其是水資源短缺和不合理地開發(fā)利用使得水資源脆弱性發(fā)生改變。因此，進(jìn)行黃河三角洲水資源脆弱性評(píng)價(jià)，為解決黃河三角洲地區(qū)水資源短缺和環(huán)境脆弱性等問(wèn)題提供參考，對(duì)實(shí)現(xiàn)黃河三角洲水資源的可持續(xù)利用具有重要意義。

水資源脆弱性源于Margat在1968年提出的地下水脆弱性這一概念[1]，國(guó)外水文地質(zhì)學(xué)家開始對(duì)地下水脆弱性進(jìn)行研究，并逐步明晰地下水脆弱性概念內(nèi)涵及評(píng)價(jià)方法[2-3]，水資源脆弱性評(píng)價(jià)逐漸拓展為地下水和地表水綜合評(píng)價(jià)。我國(guó)對(duì)水資源脆弱性評(píng)價(jià)研究起步較晚，大約始于20世紀(jì)90年代中期，在水資源脆弱性定義上多是引用外文資料。而對(duì)于水資源脆弱性(vulnerability of water resources) 的研究則主要探討的是氣候變化下對(duì)水資源系統(tǒng)的影響，考慮的影響因素包括氣候變化和人類活動(dòng)的影響[4]。

當(dāng)前，水資源脆弱性評(píng)價(jià)不僅僅局限于點(diǎn)域評(píng)價(jià)，更多關(guān)注于某一區(qū)域空間水資源脆弱性分布格局，為進(jìn)一步實(shí)現(xiàn)水資源空間優(yōu)化配置和合理開發(fā)利用奠定基礎(chǔ)。本文基于GIS空間分析技術(shù)初步實(shí)現(xiàn)評(píng)價(jià)區(qū)域劃分，并建立黃河三角洲評(píng)價(jià)指標(biāo)體系，采用模糊評(píng)價(jià)方法識(shí)別黃河三角洲空間水資源脆弱性。

1 研究區(qū)概況

黃河三角洲總面積13230.8km2，行政區(qū)劃上包括山東省東營(yíng)市全部5個(gè)縣區(qū)(占93%)、濱州市的沾化縣、無(wú)棣縣和壽光北部地區(qū)(占7%)。因受基礎(chǔ)資料限制，筆者選擇占黃河三角洲面積93%的東營(yíng)市作為研究區(qū)域，對(duì)區(qū)域水資源脆弱性進(jìn)行評(píng)價(jià)。

研究區(qū)多年平均降雨量537.7mm，時(shí)空分布不均勻，年內(nèi)分配亦不均勻，汛期降水占全年的75.2%，冬、春季節(jié)降水較少。黃河是該地區(qū)最主要的可開發(fā)利用的客水資源，全區(qū)90%以上的用水需要由黃河來(lái)供給。研究區(qū)鹽堿地分布面積大，大部分地區(qū)地下水為咸水，僅在小清河以南的廣饒縣、黃河灘地及黃河故道的部分地段有地下淡水分布。在過(guò)去的一段時(shí)間內(nèi)，為了增加供水量，地下水遭受了高強(qiáng)度的開采甚至超采。

2 水資源脆弱性評(píng)價(jià)指標(biāo)體系構(gòu)建

2.1 評(píng)價(jià)指標(biāo)體系

水資源脆弱性評(píng)價(jià)指標(biāo)體系涉及范圍和內(nèi)容比較多。為便于評(píng)價(jià)分析，筆者根據(jù)指標(biāo)體系構(gòu)建的科學(xué)性、完備性、主導(dǎo)性、相互獨(dú)立性和區(qū)域性等原則，在實(shí)地調(diào)查研究區(qū)水文地質(zhì)、地形、氣候等自然條件及社會(huì)經(jīng)濟(jì)狀況，并參考文獻(xiàn)[5-7]，從自然脆弱性和人為脆弱性方面進(jìn)行評(píng)價(jià)指標(biāo)體系構(gòu)建。研究區(qū)域水資源脆弱性評(píng)價(jià)指標(biāo)體系如圖1所示。

圖1 研究區(qū)水資源脆弱性指標(biāo)體系

2.2 GIS評(píng)價(jià)分區(qū)

為實(shí)現(xiàn)水資源脆弱性空間評(píng)價(jià)，需進(jìn)行評(píng)價(jià)區(qū)域的空間劃分。依據(jù)《山東省環(huán)境地質(zhì)圖集》[8]、《東營(yíng)市水資源公報(bào)(2008)》等相關(guān)數(shù)據(jù)資料，獲得評(píng)價(jià)指標(biāo)屬性值的空間分布圖，限于本文篇幅僅給出地下水礦化度的指標(biāo)分布(圖2)。然后利用Map GIS軟件將各指標(biāo)參數(shù)圖層文件進(jìn)行疊加處理，實(shí)現(xiàn)評(píng)價(jià)區(qū)域空間劃分，研究區(qū)域疊加得到180個(gè)分區(qū)，如圖3所示。

圖2 地下水礦化度空間分布(單位：g/L)

圖3 GIS評(píng)價(jià)分區(qū)圖

3 水資源脆弱性評(píng)價(jià)模型構(gòu)建

設(shè)由n個(gè)(n=180)待優(yōu)選的樣本(GIS疊加分區(qū))，m個(gè)(m=15)樣本脆弱性指標(biāo)，構(gòu)成樣本的指標(biāo)矩陣

(i=1，2，…，m；j=1，2，…，n)

(1)

式中：xij為方案i中第j個(gè)目標(biāo)的特征值。

表2 評(píng)價(jià)分區(qū)脆弱性評(píng)價(jià)結(jié)果

水資源脆弱性評(píng)價(jià)指標(biāo)相對(duì)隸屬度計(jì)算主要有兩類，一類是目標(biāo)特征值越大方案越優(yōu)(C1、C2、C3、C4、C5、C7、C15)，采用公式(2)計(jì)算；另一類是目標(biāo)特征值越小方案越優(yōu)(C6、C8、C9、C10、C11、C12、C13、C14)，采用公式(3)計(jì)算。

越大越優(yōu)型：

(2)

越小越優(yōu)型：

(3)

式中：rij為第i個(gè)樣本的第j項(xiàng)指標(biāo)對(duì)脆弱性的相對(duì)隸屬度；ximax、ximin分別為方案i對(duì)目標(biāo)j的特征值的最大、最小值。

R為n個(gè)樣本m項(xiàng)指標(biāo)的相對(duì)隸屬度矩陣，如公式(4)所示。

(4)

采用熵權(quán)法[9]計(jì)算評(píng)價(jià)指標(biāo)的權(quán)重W

W=(wi)1×m

(5)

(6)

將指標(biāo)相對(duì)隸屬度rij及指標(biāo)權(quán)重wi代入模糊優(yōu)選模型[10]，如式(7)所示，即可對(duì)180個(gè)疊加分區(qū)進(jìn)行水資源脆弱性評(píng)價(jià)。

(7)

4 黃河三角洲水資源脆弱性評(píng)價(jià)

根據(jù)上述構(gòu)建的水資源脆弱性評(píng)價(jià)模型，按照式(5)計(jì)算出各指標(biāo)權(quán)重，如表1所示。

表1 評(píng)價(jià)指標(biāo)的權(quán)重

由表1中可知，指標(biāo)分類中人為脆弱性指標(biāo)權(quán)重比自然脆弱性指標(biāo)權(quán)重大。其中，地表水污染指數(shù)C13、單位面積平均用水量C15、地下水開采模數(shù)C9、外調(diào)水占供水量比例C14、含水層補(bǔ)給模數(shù)C4、工程調(diào)蓄能力C7、萬(wàn)元GDP用水量C12、礦化度C6占權(quán)重的前八位。

計(jì)算出各指標(biāo)的權(quán)重及相對(duì)隸屬度矩陣后，根據(jù)式(7)計(jì)算出各分區(qū)的脆弱性隸屬度，計(jì)算結(jié)果如表2所示。

根據(jù)最大隸屬度原則及研究區(qū)實(shí)際概況，將各分區(qū)水資源脆弱性評(píng)價(jià)隸屬度級(jí)別分為5個(gè)等級(jí)：①?gòu)?qiáng)脆弱性：相對(duì)隸屬度0.1～0.3；②較強(qiáng)脆弱性：相對(duì)隸屬度0.3～0.4；③中度脆弱性：相對(duì)隸屬度0.4～0.6；④低脆弱性：相對(duì)隸屬度0.6～0.7；⑤弱脆弱性：相對(duì)隸屬度0.7～1.0。黃河三角洲水資源脆弱性評(píng)價(jià)結(jié)果如圖4所示。

圖4 研究區(qū)水資源脆弱性評(píng)價(jià)結(jié)果

評(píng)價(jià)結(jié)果顯示，研究區(qū)整體水資源脆弱性程度較高，其中中度脆弱以上地區(qū)占研究區(qū)總面積的93.7%，沒(méi)有弱度脆弱的區(qū)域。水資源脆弱性空間變化較大，沿海及小清河區(qū)域水資源脆弱性較高，沿黃河區(qū)域脆弱性相對(duì)較低。分析原因是沿海地區(qū)受到海水入侵，礦化度較高，對(duì)跨流域外調(diào)水依賴性較大；小清河地區(qū)受到地下水過(guò)度開采及地表水污染嚴(yán)重及農(nóng)業(yè)面源污染等影響，水資源脆弱性較高；而地處沿黃河區(qū)域水資源相對(duì)較為豐富，引、輸水干渠配套工程較為完善，其水資源脆弱性相對(duì)較低，處于中等水平。

5 建 議

水資源脆弱性評(píng)價(jià)是識(shí)別區(qū)域空間水資源受到人類活動(dòng)干擾后恢復(fù)原來(lái)狀態(tài)的能力，為進(jìn)一步實(shí)現(xiàn)水資源空間優(yōu)化配置和合理開發(fā)利用奠定基礎(chǔ)。研究表明，黃河三角洲脆弱性較高，空間區(qū)域變化較大，水資源脆弱性影響因素較多且不同。因此需要在現(xiàn)有黃河三角洲水系基礎(chǔ)上，規(guī)劃建設(shè)濱海骨干生態(tài)河道，綜合整治區(qū)域排水河道，完善灌區(qū)灌排干支渠，增設(shè)必要蓄水設(shè)施，把黃河、區(qū)域河道和平原水庫(kù)等水源工程串聯(lián)起來(lái)，構(gòu)成一個(gè)集輸水、蓄水、灌溉、供水、排水、改堿、咸水微咸水利用及自然保護(hù)區(qū)濕地補(bǔ)水于一體的工程水網(wǎng)和生態(tài)水系，做到多水聯(lián)供和循環(huán)利用，切實(shí)轉(zhuǎn)變水資源利用方式，減弱黃河三角洲水資源脆弱性，建立黃河三角洲水資源高效生態(tài)利用模式。

[1] DOERFLIGER N，JEANNIN P Y，ZWAHLEN F.Water vulnerability assessment in karst environments a new method of defining protection areas using a multi-attribute approach and GIS tools[J].Environmental Geology，1999，39(2)：165-176.

[2] 李劍穎.官?gòu)d水庫(kù)流域水資源脆弱性評(píng)價(jià)研究[D].北京：北京師范大學(xué)，2007.

[3] 張保祥.黃水河流域地下水脆弱性評(píng)價(jià)與水源保護(hù)區(qū)劃分研究[D].北京：中國(guó)地質(zhì)大學(xué)，2006.

[4] 沈珍瑤，楊志峰，曹瑜.環(huán)境脆弱性研究評(píng)述[J].地質(zhì)科技情報(bào)，2003，22(3)：91-94.

[5] 劉綠柳.水資源脆弱性及其定量評(píng)價(jià)[J].水土保持通報(bào)，2002，22(2)：41-44.

[6] 鄒君，劉蘭芳，田亞平,等.地表水資源的脆弱性及其評(píng)價(jià)初探[J].資源科學(xué)，2007，29(1)：92-98.

[7] 馮少輝，李靖，朱振峰，等.云南省滇中地區(qū)水資源脆弱性評(píng)價(jià)[J].水資源保護(hù)，2010，26(1)：13-16.

[8] 山東省地質(zhì)礦產(chǎn)局.山東省環(huán)境地質(zhì)圖集[M].濟(jì)南：山東省地圖出版社，1996：66.

[9] 張先起，梁川.基于熵權(quán)的模糊物元模型在水質(zhì)綜合評(píng)價(jià)中的應(yīng)用[J].水利學(xué)報(bào)，2005，36(9)：1057-1061.

[10] 陳守煜.復(fù)雜水資源系統(tǒng)優(yōu)化模糊識(shí)別理論與應(yīng)用[M].長(zhǎng)春：吉林大學(xué)出版社，2002：10-13.