王金哲,張光輝,聶振龍,嚴(yán)明疆,王瑩

(中國地質(zhì)科學(xué)院水文地質(zhì)環(huán)境地質(zhì)研究所,河北石家莊050061)

淺層地下水系統(tǒng)的變化不僅與當(dāng)?shù)氐淖匀坏乩憝h(huán)境有關(guān),也和人類活動強(qiáng)度密切相關(guān)。目前,大規(guī)模人類活動無時不在悄然改變著天然水循環(huán)的大氣、地表、土壤和地下各個過程,致使現(xiàn)代環(huán)境下水循環(huán)呈現(xiàn)出明顯的“天然—人工”二元特性[1-2]。但由于人類活動是一個抽象的綜合性概念,許多方面量化比較困難,因而以往的研究多為定性或集中于單一因素的分析[3-5]。如何將人類活動對水循環(huán)的干擾影響分離出來從而直觀認(rèn)識其作用程度已成為近幾年的研究熱點(diǎn)?；诖?近幾年也有學(xué)者在人類活動對水循環(huán)的影響量化研究進(jìn)行了有益的探索,采用的方法歸納起來有地球動力學(xué)方法[6]、人工神經(jīng)網(wǎng)絡(luò)缺省因子檢驗(yàn)法[7]、變異系數(shù)法[8]和投影尋蹤回歸模型[9]。但不同的研究方法有不同的適用范圍和優(yōu)缺點(diǎn),從計(jì)算精度、步驟繁瑣程度、指標(biāo)容量等方面考慮,本研究嘗試采用加權(quán)綜合指數(shù)法,選取研究區(qū)影響淺層地下水系統(tǒng)的典型人類活動作為指標(biāo),以MAPGIS為數(shù)據(jù)處理平臺,建立了人類活動對淺層地下水干擾程度量化研究的指標(biāo)體系,并以滹沱河流域平原區(qū)為例,進(jìn)行了實(shí)例驗(yàn)證。

1 評價方法

1.1 方法原理

加權(quán)綜合指數(shù)法是通過加權(quán)來測定一組項(xiàng)目的綜合變動狀況。該方法將與評價目標(biāo)有關(guān)的所有要素綜合整理在一起,確定各要素的相對重要性,給出定量指標(biāo),然后通過數(shù)學(xué)方法求解,作為綜合分析、評價的基礎(chǔ)[10]。在本項(xiàng)研究中,綜合指數(shù)法就是將影響淺層地下水系統(tǒng)的分指標(biāo)指數(shù)和權(quán)重進(jìn)行迭加形成綜合指數(shù),根據(jù)綜合指數(shù)量化評價人類活動對淺層地下水系統(tǒng)造成的干擾程度。

其評價模型如下:

式中:R——綜合評價指數(shù);ai——評價參數(shù)的權(quán)重;Xi——評價參數(shù);n——評價參數(shù)的個數(shù)。

1.2 權(quán)重計(jì)算

權(quán)重計(jì)算時選用灰色關(guān)聯(lián)度法。其計(jì)算模型為:

設(shè)有m個與母因素(x0)有一定關(guān)聯(lián)作用的子因素(x1,x2,…,xm),每個評價因子都有n個統(tǒng)計(jì)值,構(gòu)成母序列和子序列。

母、子序列在t=1到t=n的關(guān)聯(lián)性用關(guān)聯(lián)系數(shù)表示:

式中:Δmin,Δmax——m條子線在區(qū)間[1,n]母線的距離 Δok(i)的最大值與最小值;ζ——分辨系數(shù),一般取0.5。

于是有第k條子線與母線在[1,n]間的關(guān)聯(lián)度為:

采用公式(4)使關(guān)聯(lián)度和為“1”對關(guān)聯(lián)度進(jìn)行標(biāo)準(zhǔn)化,標(biāo)準(zhǔn)化后的關(guān)聯(lián)度即可做為每個評價因子的權(quán)重。

1.3 綜合評價

根據(jù)計(jì)算出來的一系列R值,對照相應(yīng)的評價分級標(biāo)準(zhǔn),然后根據(jù)評價區(qū)的具體條件進(jìn)行綜合分析,判別評價對象的真實(shí)狀況,判別依據(jù)見表1。一般認(rèn)為,人類活動對地下水干擾程度大于0.5時,人類活動對地下水動態(tài)起主導(dǎo)作用[6]。

表1 人類活動對淺層地下水干擾程度綜合指數(shù)等級劃分及其含義

2 研究區(qū)概況

滹沱河流域平原區(qū)位于河北省中南部,地處滹沱河以南、滏陽河以北、太行山以東山前傾斜平原,主要包括石家莊地區(qū)和衡水地區(qū)的大部分縣市,邢臺地區(qū)的柏鄉(xiāng)、新河兩縣小部分區(qū)域包括在內(nèi)?？偯娣e8 204 km2。

該區(qū)屬溫帶半濕潤、半干旱大陸性季風(fēng)氣候,降水量年際變化較大,多年平均降雨量493 mm。滹沱河曾是本區(qū)最大河流,屬海河的子牙河水系,自西向東貫穿該區(qū),在獻(xiàn)縣與滏陽河交匯,流域面積2.48×104km2,目前已干涸。黃壁莊水庫和崗南水庫是位于滹沱河出山口的兩個串聯(lián)水庫,庫容分別是1.50×109m3和 1.20×109m3。

3 評價體系構(gòu)建及數(shù)據(jù)處理

3.1 評價體系構(gòu)建

體系的構(gòu)建首先需選取指標(biāo),指標(biāo)選取依據(jù)主導(dǎo)型、可度量、可操作性和覆蓋面廣的原則,并考慮可直觀反映研究區(qū)地下水系統(tǒng)的現(xiàn)狀因子,影響和反映人類活動的簡單易解釋因子,人類活動與地下水系統(tǒng)在內(nèi)的相互影響、相互制約的表征因子,并能通過更新數(shù)據(jù)就可以調(diào)整為新的指標(biāo)[11]。從上述因素考慮,選取了8個指標(biāo)作為滹沱河流域平原區(qū)人類活動對淺層地下水干擾程度量化研究的評價因子,分別為開采量、人口、機(jī)井?dāng)?shù)量、糧食產(chǎn)量、農(nóng)作物播種面積、蔬菜產(chǎn)量、水澆地面積和工業(yè)總產(chǎn)值。

3.2 數(shù)據(jù)處理

3.2.1 指標(biāo)變量數(shù)據(jù)整理 滹沱河流域平原區(qū)包含19個縣市。指標(biāo)的原始數(shù)據(jù)變量序列搜集整理于1958—2008年石家莊市、衡水和邢臺市統(tǒng)計(jì)年鑒,通過核算8個指標(biāo)變量分別在20世紀(jì)50—90年代和21世紀(jì)初期的平均值,得到8個指標(biāo)在不同時代的變量值。

3.2.2 變量數(shù)據(jù)標(biāo)準(zhǔn)化 為使8個指標(biāo)之間具有可比性,從而完整地組合到一起,首先統(tǒng)一這些指標(biāo)的量綱。采用極值標(biāo)準(zhǔn)化方法對指標(biāo)進(jìn)行標(biāo)準(zhǔn)化。

3.2.3 變量數(shù)據(jù)提取 根據(jù)研究區(qū)現(xiàn)有資料局限、研究程度、面積范圍大小和評價精度要求,運(yùn)用MAPGIS軟件對研究區(qū)進(jìn)行網(wǎng)格化,選取的剖分網(wǎng)格為正方形。由于該區(qū)地下水資源評價、水文地質(zhì)條件、開發(fā)利用程度等方面研究程度較高,根據(jù)基礎(chǔ)數(shù)據(jù)的最小單元和示范區(qū)面積,確定剖分網(wǎng)格間距為1.25 km,共剖分單元128×76個,其中有效計(jì)算單元共計(jì)5 645個。

3.3 確定權(quán)重

采用灰色關(guān)聯(lián)度法確定權(quán)重。通常情況下,淺層地下水的變化指標(biāo)用地下水水位或地下水位埋深表示,在本研究中,用地下水開采量作為淺層地下水系統(tǒng)變化的特征值。原因有兩個方面:一是根據(jù)收集資料的實(shí)際情況,指標(biāo)變量都是以縣為單位的多年動態(tài)變化值,每個年代人類活動對淺層地下水的干擾都是在前期各項(xiàng)人類活動多年累計(jì)基礎(chǔ)上疊加的結(jié)果。即各項(xiàng)指標(biāo)反映的是以縣為單位的區(qū)域值,而通常意義上的地下水位或地下水埋深是反映某一點(diǎn)的值,不能把人類活動對區(qū)域地下水系統(tǒng)的干擾很好地反映出來。二是地下水開采量是與地下水互動相關(guān)的變量,能較準(zhǔn)確地表征淺層地下水系統(tǒng)的變化,而且一個地區(qū)的地下水開采量是用這一地區(qū)內(nèi)不同位置點(diǎn)的開采量值進(jìn)行累加得到的變量。所以,用地下水開采量作為淺層地下水系統(tǒng)變化的特征值,也就是母序列。子序列為7個指標(biāo)變量,經(jīng)過計(jì)算,得到人類活動對淺層地下水干擾程度評價7個指標(biāo)的權(quán)重值(表2)。

表2 人類活動各指標(biāo)對淺層地下水的干擾權(quán)重

3.4 評價指數(shù)計(jì)算

在取得各指標(biāo)的基本參數(shù)和權(quán)重后,運(yùn)用公式(1)就可得到人類活動對淺層地下水的干擾程度這一綜合指數(shù)。

把不同時代的綜合指數(shù)與相對應(yīng)的坐標(biāo)建立空間對應(yīng)關(guān)系,借助MAPGIS工作平臺,使用空間分析的DTM功能就可繪制工作區(qū)人類活動對淺層地下水干擾程度的分布圖。

4 研究區(qū)人類活動對淺層地下水干擾程度定量評價

本研究中繪制了20世紀(jì)50—90年代和21世紀(jì)初期的人類活動對淺層地下水干擾程度分布(圖1),通過對5個時段進(jìn)行評價,可縱觀20世紀(jì)50年代以來人類活動對滹沱河流域平原區(qū)淺層地下水系統(tǒng)干擾程度的演變特征。

20世紀(jì)50年代,滹沱河流域平原區(qū)農(nóng)業(yè)落后,工業(yè)不發(fā)達(dá),河道自然溢流,濕地遍布,淺層地下水的開采量為1.35×109m3。這一時期人類活動為低基數(shù)發(fā)展階段,對淺層地下水干擾程度綜合指數(shù)變化區(qū)間為0～0.4,都小于0.5,也就是,淺層地下水系統(tǒng)保持自然平衡狀態(tài),以自然因素調(diào)節(jié)為主,人類活動沒有占據(jù)主導(dǎo)地位。人類活動對淺層地下水干擾程度0～0.2之間的分布面積為2 461.5 km2,主要分布在研究區(qū)的中東部;0.2～0.4之間的分布面積為6 343.5 km2,主要分布研究區(qū)的西部地區(qū)。這一時期,人類活動對淺層地下水系統(tǒng)的干擾程度較弱,大部分地區(qū)的淺層地下水系統(tǒng)自身修復(fù)能力較強(qiáng)。

20世紀(jì)60年代,“以蓄為主”推動的第一次水利建設(shè)高潮和圍繞建成“上蓄、中疏、下排、適當(dāng)?shù)販狈篮轶w系開始的第二次水利建設(shè)高潮,使研究區(qū)的地下水系統(tǒng)不再處于平靜的自然平衡狀態(tài),自然狀態(tài)下的補(bǔ)、徑、排關(guān)系被破壞。但受當(dāng)時經(jīng)濟(jì)技術(shù)條件的限制,人類活動對淺層地下水的干擾程度較弱,淺層地下水的開采量為1.71×109m3,地下水系統(tǒng)可以修復(fù),之后建立新的平衡。這一時期,人類活動對淺層地下水的干擾程度綜合指數(shù)變化區(qū)間為0～0.4,0～0.2之間的分布面積為1 750.5 km2,主要分布在中西部;0.2～0.4之間的分布面積為6 343.5 km2,主要分布在東部。與20世紀(jì)50年代相比,干擾程度低的區(qū)域在減小,干擾程度較低的區(qū)域在逐漸增大,也就是說,人類活動對淺層地下水的干擾程度增強(qiáng),淺層地下水系統(tǒng)對遭到的破壞可以自身修復(fù)。

圖1 滹沱河流域平原區(qū)不同時代人類活動對淺層地下水干擾程度分布圖

20世紀(jì)70年代,圍湖造田和大規(guī)模的打井高潮,以及工業(yè)的復(fù)蘇使人類活動出現(xiàn)高潮,對地下水的影響明顯增強(qiáng),淺層地下水開采量增加為2.55×109m3,表現(xiàn)為水澆地面積的增加使得農(nóng)業(yè)對地下水的開采量明顯增加,20世紀(jì)60年代中期水澆地面積近為3.00×104hm2,20世紀(jì)70年代末為2.03×105hm2。這一時期,滹沱河流域平原區(qū)人類活動對淺層地下水干擾程度綜合指數(shù)變化區(qū)間為0～0.6,出現(xiàn)了人類活動對淺層地下水干擾占主導(dǎo)作用的大面積區(qū)域,面積有5 321.4 km2,位于研究區(qū)的西部和中部大部分地區(qū)。干擾程度低的區(qū)域急劇縮減,僅為79.6 km2,幾乎消失殆盡。與20世紀(jì)60年代相比,人類活動對淺層地下水的干擾程度仍在增強(qiáng),大部分地區(qū)的淺層地下水系統(tǒng)修復(fù)需要較長時間。

20世紀(jì)80年代,研究區(qū)水利化程度進(jìn)一步發(fā)展,機(jī)井?dāng)?shù)量大幅增加,糧食產(chǎn)量增長效果顯著,與之相應(yīng),淺層地下水系統(tǒng)受到強(qiáng)烈的擾動,淺層地下水開采量為3.03×109m3,是研究時段的最大消耗時期。這一時期,滹沱河流域平原區(qū)人類活動對淺層地下水干擾程度綜合指數(shù)變化區(qū)間為0.2～0.7,人類活動對淺層地下水的干擾強(qiáng)度向更強(qiáng)方向發(fā)展,出現(xiàn)干擾程度較強(qiáng)的區(qū)域,面積為1 414.4 km2,主要分布在中部和西部。而且,干擾程度低的區(qū)域完全消失,干擾程度較低至一般的區(qū)域在20世紀(jì)70年代基礎(chǔ)上有所減小,減小面積為1 414.4 km2,即這部分縮減的區(qū)域轉(zhuǎn)化為干擾程度較強(qiáng)的區(qū)域。與20世紀(jì)70年代相比,人類活動對淺層地下水的干擾程度不斷增強(qiáng),地下水系統(tǒng)遭到的破壞不斷增強(qiáng),可以自身修復(fù)的區(qū)域不斷縮小。

20世紀(jì)90年代,隨著地下水的資源功能、生態(tài)功能和地質(zhì)環(huán)境功能等諸多問題的出現(xiàn),人類活動被迫進(jìn)入理性高基數(shù)發(fā)展階段,主動探求更好的水資源利用配置,如1995年的黃壁莊水庫水引入石家莊市、合理井灌和渠灌等,但受區(qū)內(nèi)地理位置和缺少地表水的現(xiàn)實(shí),造成人們?nèi)匀徊坏貌灰蚤_采地下水為主,淺層地下水的開采量為2.82×109m3。這一時期,人類活動對淺層地下水干擾程度綜合指數(shù)變化區(qū)間為0.2～0.8,干擾程度繼續(xù)增強(qiáng),干擾程度占主導(dǎo)作用的區(qū)域繼續(xù)迅速擴(kuò)展,出現(xiàn)干擾程度強(qiáng)的區(qū)域,位于石家莊市。干擾程度較強(qiáng)的區(qū)域面積從20世紀(jì)70年代的0 km2增至90年代的4 899.1 km2,主要位于研究區(qū)的西部和中部地區(qū),干擾程度低的區(qū)域萎縮至研究區(qū)東部的小部分區(qū)域,干擾程度較低的區(qū)域萎縮的速度加快,與20世紀(jì) 80年代相比,從 3 404 km2縮減到1 297.7 km2,減少了2 106.3 km2。干擾程度低的區(qū)域完全消失。說明人類活動對淺層地下水的干擾程度持續(xù)增強(qiáng),地下水系統(tǒng)遭到的破壞持續(xù)增強(qiáng),絕大部分地區(qū)的淺層地下水系統(tǒng)修復(fù)需要較長時間,甚至不能修復(fù)。

進(jìn)入21世紀(jì),滹沱河流域平原區(qū)實(shí)施了農(nóng)業(yè)上的節(jié)水降耗措施和工業(yè)生產(chǎn)用水的循環(huán)率等技術(shù)推廣應(yīng)用,石家莊市作為地下水的用水大戶也從20世紀(jì)90年代開始加大地表水入市力度,這些措施的實(shí)施,使人類活動對淺層地下水的干擾程度有所降低,淺層地下水的開采量為2.77×109m3,相比20世紀(jì)90年代有所減小。但由于前期的發(fā)展和積累,人類活動強(qiáng)度對淺層地下水的干擾強(qiáng)度在20世紀(jì)90年代的基礎(chǔ)上繼續(xù)增大。目前,綜合指數(shù)變化區(qū)間為0.2～1.0,0.4～1.0之間分布的區(qū)域幾乎占據(jù)整個研究區(qū),面積為 7 374.6 km2,占研究區(qū)總面積的83.8%。干擾程度強(qiáng)的區(qū)域?yàn)橹猩嫌蔚氖仪f市和正定縣、辛集市、寧晉縣和深州市的縣城位置,呈斑塊狀分布,面積有1 624 km2。說明人類活動對淺層地下水的干擾程度持續(xù)增強(qiáng),地下水系統(tǒng)遭到的破壞更為嚴(yán)重,不能修復(fù)的區(qū)域逐漸擴(kuò)展。

從20世紀(jì)50年代到21世紀(jì)初期,人類活動對滹沱河流域平原淺層地下水的干擾程度在時空分布上呈現(xiàn)出明顯的特征?？臻g上,不同時代都表現(xiàn)出從西部向東部干擾程度逐漸減小;從整個研究區(qū)看,干擾程度增大的趨勢從西部向東部逐漸推進(jìn);干擾程度從分級簡單的0～0.4發(fā)展到目前0.2～1.0不同級別共存的復(fù)雜分布狀態(tài)。時間上,隨著時間的推移,整個研究區(qū)的干擾程度都表現(xiàn)出增強(qiáng)態(tài)勢,在不同地區(qū)增強(qiáng)幅度不一致,人類活動對淺層地下水干擾影響經(jīng)歷了不占據(jù)主導(dǎo)作用→逐漸占據(jù)主導(dǎo)作用→幾乎完全占據(jù)主導(dǎo)作用的發(fā)展過程。

5 結(jié)論

采用加權(quán)綜合指數(shù)法量化評價人類活動對區(qū)域淺層地下水干擾程度的方法,根據(jù)主導(dǎo)性、可度量性、可操作和覆蓋面廣的原則,選取了量化評價指標(biāo),在此基礎(chǔ)上,對近50 a來滹沱河流域平原區(qū)人類活動對淺層地下水的干擾程度進(jìn)行了量化評價?？偟膩砜?人類活動對滹沱河流域平原區(qū)淺層地下水干擾程度時空分布上呈現(xiàn)出明顯的特征。

在空間上,不同時代都表現(xiàn)出從西部向東部干擾程度逐漸減小;從整個研究區(qū)發(fā)展形勢看,干擾程度增大的趨勢從西部向東部逐漸推進(jìn);干擾程度從分級簡單的0～0.4發(fā)展到目前0.2～0.4,0.4～0.6,0.6～0.8和0.8～1.0不同級別共存的復(fù)雜分布狀態(tài)。

在時間上,隨著時間的推移,整個研究區(qū)的干擾程度都表現(xiàn)出增強(qiáng)態(tài)勢,在不同地區(qū)增強(qiáng)幅度不一致;人類活動對淺層地下水干擾影響經(jīng)歷了從不占據(jù)主導(dǎo)作用到逐漸占據(jù)主導(dǎo)作用再到幾乎完全占據(jù)主導(dǎo)作用的發(fā)展過程。

[1] 王浩,王建華,秦大庸,等.基于二元水循環(huán)模式的水資源評價理論方法[J].水利學(xué)報(bào),2006,37(12):1496-1502.

[2] 王建華,江東.黃河流域二元水循環(huán)要素反演研究[M].北京:科學(xué)出版社,2005:89-90.

[3] 陳志清.歷史時期黃河下游的淤積決口改道與人類活動的關(guān)系[J].地理科學(xué)進(jìn)展,2001,20(1):44-50.

[4] 李新,周宏飛.人類活動干預(yù)后的塔里木河水資源可持續(xù)利用問題[J].地理研究,1998,17(2):171-177.

[5] 劉昌明,成立.黃河干流下游斷流的徑流序列分析[J].地理學(xué)報(bào),2000,55(3):257-265.

[6] 楊建強(qiáng),羅先香,高振會.GIS支持下人類活動對地下水動態(tài)影響的定量分析[J].水科學(xué)進(jìn)展,2003,14(3):358-362.

[7] 馮紹元,霍再林,康少忠,等.干旱內(nèi)陸區(qū)自然—人工條件下地下水動態(tài)的ANN模型[J].水利學(xué)報(bào),2007,38(7):873-878.

[8] 許月卿.土地利用對地下水位的影響[J].地理研究,2005,24(2):222-228.

[9] 張翠云,王昭.黑河流域人類活動強(qiáng)度的定量評價[J].地球科學(xué)進(jìn)展,2004(S):387-390.

[10] 賈俊平.統(tǒng)計(jì)學(xué)[M].北京:清華大學(xué)出版社,2004:395-396.

[11] 夏軍,黃國和,龐進(jìn)武,等.可持續(xù)水資源管理:理論?方法?應(yīng)用[M].北京:化學(xué)工業(yè)出版社,2005:152-156.