何 曄

（1. 上海市地質(zhì)調(diào)查研究院，上海 200072，2. 自然資源部地面沉降監(jiān)測與防治重點(diǎn)實(shí)驗(yàn)室，上海 200072，3. 上海地面沉降控制工程技術(shù)研究中心，上海 200072)

地下水作為水資源的重要組成部分，既是維系區(qū)域生態(tài)環(huán)境的重要自然要素，也是保障社會(huì)穩(wěn)定不可替代的基礎(chǔ)支柱。然而全球地下水資源在氣候變化和人類活動(dòng)的影響下，消耗十分嚴(yán)重。近年來隨著城市發(fā)展及生活水平的不斷提高，地下水資源的供需矛盾日益突出，而地下水超采問題也日益嚴(yán)重，已經(jīng)引起了諸如地面沉降、地裂縫、水質(zhì)惡化和生態(tài)環(huán)境破壞等一系列與地下水超采相關(guān)的問題[1-6]。為加強(qiáng)地下水資源管理，保護(hù)優(yōu)質(zhì)的地下水資源，實(shí)現(xiàn)地下水資源可持續(xù)開發(fā)利用，保障經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展，需開展地下水超采區(qū)劃定工作[7-9]，全面掌握地下水超采區(qū)現(xiàn)狀和發(fā)展趨勢(shì)，加快推進(jìn)地下水超采區(qū)治理。

地下水超采區(qū)的劃定方法也越來越趨于多樣化并日益成熟，從最開始單一的開采系數(shù)法逐步演變?yōu)樗粍?dòng)態(tài)法、開采系數(shù)法和引發(fā)問題法等多種方法[10]。但無論水位動(dòng)態(tài)法、開采系數(shù)法還是引發(fā)問題法都無法全面反映地下水超采情況，本次利用層次分析法建立上海市地下水超采區(qū)綜合劃定框架，對(duì)水位變化、開采系數(shù)以及開采引發(fā)的問題進(jìn)行綜合分析，以期為地下水超采區(qū)劃定及科學(xué)合理制定地下水資源保護(hù)措施提供依據(jù)和參考。

1 區(qū)域概況

松散巖類孔隙水是上海地區(qū)內(nèi)主要的地下水類型，蘊(yùn)含豐富的地下水資源，是本次地下水超采區(qū)劃定的對(duì)象。按照地質(zhì)時(shí)代、水動(dòng)力條件和成因類型，將第四系松散巖類孔隙含水層劃定為全新統(tǒng)潛水含水層和微承壓含水層（僅局部分布），上更新統(tǒng)第一、二承壓含水層，中更新統(tǒng)第三承壓含水層，下更新統(tǒng)第四、五承壓含水層。

上海淺層地下水資源的開發(fā)利用主要以鄉(xiāng)村分散的土井開采為主，目前全市淺層地下水的開采程度很低，主要為以基坑工程為主的工程建設(shè)活動(dòng)的降排地下水量。全市地下水資源開發(fā)利用以深層承壓水（第二至第五承壓含水層）為主。自市區(qū)1860年開鑿第一口深井以來，深井井?dāng)?shù)及其年開采量曾一度持續(xù)增長，1965年之前，全市開采量達(dá)1.35～2.03億m3/a，達(dá)到歷史的高峰，導(dǎo)致市區(qū)嚴(yán)重的地面沉降災(zāi)害，因此1965年后進(jìn)行壓縮開采，至70年代末開采量相對(duì)穩(wěn)定，80年代后，需水量猛增的跡象加劇，而且在開采層次及地域上有所變化；隨著集約化供水能力不斷增大，目前全市地下水年度開采量持續(xù)壓減，且幅度較大（如圖1）。

圖1 上海市承壓含水層地下水開采量歷年變化Fig.1 Diachronic change of groundwater extraction from confined aquifers in Shanghai

2 傳統(tǒng)地下水超采區(qū)劃定

根據(jù)《地下水超采區(qū)評(píng)價(jià)導(dǎo)則》[11]，傳統(tǒng)地下水超采區(qū)劃定主要采用水位動(dòng)態(tài)法、開采系數(shù)法和引發(fā)問題法三種方法。數(shù)據(jù)統(tǒng)計(jì)分析以行政區(qū)域?yàn)閱卧?，其中上海市區(qū)（黃浦區(qū)、靜安區(qū)、徐匯區(qū)、長寧區(qū)、虹口區(qū)、普陀區(qū)、楊浦區(qū)）由于面積較小，合并為中心城區(qū)。

2.1 水位動(dòng)態(tài)法

水位動(dòng)態(tài)法是以地下水水位年均下降速率為評(píng)判指標(biāo)進(jìn)行超采區(qū)劃定的方法。因地下水開采造成地下水水位呈持續(xù)下降趨勢(shì)劃定為一般超采區(qū)，淺層孔隙水水位年均下降速率大于1.0m/a或深層孔隙水水位年均下降速率大于2.0m/a劃定為嚴(yán)重超采區(qū)。

以上海市為例，以2016年為基準(zhǔn)年，評(píng)價(jià)期為2001～2016年，選取了582口地下水位動(dòng)態(tài)監(jiān)測井，按各行政區(qū)統(tǒng)計(jì)水位動(dòng)態(tài)變化趨勢(shì)，統(tǒng)計(jì)結(jié)果如表1。

表1 2001～2016年地下水水位年均變化速率統(tǒng)計(jì)Table 1 Statistical table of the average annual rate of change of groundwater level from 2001 to 2016

從上表可知，上海市各含水層地下水水位總體呈現(xiàn)上升態(tài)勢(shì)，僅局部地區(qū)地下水水位有略微下降?？芍虾Ｊ懈骱畬訜o地下水位持續(xù)下降地區(qū)，根據(jù)水位動(dòng)態(tài)法，全市無地下水超采區(qū)。

2.2 開采系數(shù)法

開采系數(shù)法是以地下水開采系數(shù)為判別指標(biāo)進(jìn)行超采區(qū)劃定的方法。地下水實(shí)際開采量超過可開采量劃定為一般超采區(qū)，若年均地下水開采系數(shù)大于1.3，則劃定為嚴(yán)重超采區(qū)。

針對(duì)上海市實(shí)際情況，淺部潛水和第一承壓含水層地下水長期來未予開發(fā)利用，年度無計(jì)劃開采指標(biāo)；深部第二、三、四、五承壓含水層地下水資源豐富且水質(zhì)較好，是上海市常年開發(fā)利用的主采層。為此，根據(jù)開采系數(shù)法進(jìn)行地下水超采區(qū)劃定時(shí)，對(duì)于淺部潛水和第一承壓含水層不予考慮地下水超采區(qū)劃定。各行政區(qū)含水層實(shí)際開采量、可開采資源量以及開采系數(shù)如表2所示。

表2 2016年上海市第二 三 四 五含水層地下水開采系數(shù)統(tǒng)計(jì)表Table 2 Statistics of groundwater exploitation coefficients in the second, third, fourth and fifth aquifers in Shanghai in 2016

可知各區(qū)各含水層實(shí)際開采量均小于可開采量，故根據(jù)開采系數(shù)法，全市無地下水超采區(qū)。

2.3 引發(fā)問題法

引發(fā)問題法是以地下水開采誘發(fā)的生態(tài)與環(huán)境地質(zhì)問題為評(píng)判指標(biāo)進(jìn)行超采區(qū)劃定的方法。上海市因地下水開采造成的生態(tài)與環(huán)境地質(zhì)問題主要為地面沉降。根據(jù)引發(fā)問題法，因地下水開采引發(fā)了地面沉降則劃定為一般超采區(qū)；因地下水開采引發(fā)的年均地面沉降速率大于10mm/a劃定為嚴(yán)重超采區(qū)。

根據(jù)表3可知，2001～2016年上海市大部分區(qū)域多年平均沉降速率小于10mm/a，只有如崇明東灘新成陸區(qū)等局部區(qū)域沉降稍大，因此根據(jù)引發(fā)問題法劃定上海全市為一般超采區(qū)。

表3 上海市地面沉降速率統(tǒng)計(jì)表（2001～2016年）Table 3 Statistics of land subsidence rate in Shanghai (2001-2016)

2.4 傳統(tǒng)劃定結(jié)果分析

由上述結(jié)果可知，三種方法得出的劃定結(jié)果不盡相同，其中根據(jù)水位動(dòng)態(tài)法和開采系數(shù)法，上海市無超采區(qū)，而根據(jù)引發(fā)問題法，上海全市為一般超采區(qū)。這主要是由于水位動(dòng)態(tài)法和開采系數(shù)法更直接反映評(píng)價(jià)區(qū)域當(dāng)下的地下水位變化和地下水超采情況，而地面沉降相對(duì)地下水位變化具有滯后性及不可逆性，上海歷史地下水超采導(dǎo)致的累積地面沉降仍較大，且工程建設(shè)，特別是深基坑降水引發(fā)的局部區(qū)域差異沉降仍長期存在?？芍N傳統(tǒng)方法側(cè)重點(diǎn)不一樣，故考慮采用層次分析法建立地下水超采區(qū)綜合劃定框架。

3 基于層次分析法的地下水超采區(qū)綜合劃定

3.1 層次分析法概述

層次分析法（Analytical Hierarchy Process，簡稱AHP）是美國運(yùn)籌學(xué)家Saaty教授提出的一種多準(zhǔn)則決策方法[12]，它是對(duì)一些較為復(fù)雜、較為模糊的問題做出決策的簡易方法，特別適用于那些難于完全定量分析的問題，在決策工作中有廣泛的應(yīng)用。

其基本步驟主要為4步，第一步是將最終目標(biāo)分解成若干層次，如目標(biāo)層、準(zhǔn)測層、資源層和指標(biāo)層等；第二步是根據(jù)專家咨詢，將各層次的各元素進(jìn)行反復(fù)的兩兩比較，確定其權(quán)重；第三步是收集真實(shí)數(shù)據(jù)計(jì)算各指標(biāo)值，并得到最終評(píng)估系數(shù)（取值范圍為0～1）；第四步即在上述計(jì)算結(jié)果的基礎(chǔ)上進(jìn)行最后的評(píng)估。

3.2 地下水超采區(qū)綜合劃定

（1）構(gòu)建層次體系

以上海市地下水超采區(qū)劃定為例，主要分為目標(biāo)層、資源層和指標(biāo)層三層。本次以上海市地下水超采區(qū)劃定為最終目標(biāo)，即作為層次體系的目標(biāo)層；按地級(jí)行政區(qū)將地下水開采區(qū)域作為資源層；參照傳統(tǒng)劃定方法，選取地下水位、地下水開采和地面沉降等作為指標(biāo)層。

（2）確定權(quán)重

一般來說，確定權(quán)重是一個(gè)冗長的過程，首先需要找到合適的專家組，然后對(duì)各層次元素兩兩反復(fù)比較判斷，最終得出專家組可接受的一致性結(jié)果。這一過程需綜合考慮實(shí)際的地下水需求情況、對(duì)各項(xiàng)指標(biāo)的敏感程度、生態(tài)地質(zhì)環(huán)境問題以及政策因素等，工作量大且需經(jīng)費(fèi)支持。本次在無法滿足前述條件情況下，采取一種較簡單方式分配資源層（Aspects）和指標(biāo)層（Indicators）權(quán)重，公式如下：

其中：WA(i)為第i個(gè)資源的權(quán)重，WI(i,j)為第i個(gè)資源的第j個(gè)指標(biāo)的權(quán)重，N為資源數(shù)量，Ni為第i個(gè)資源的指標(biāo)數(shù)量。

上述分配權(quán)重方式，只要確定了資源、指標(biāo)數(shù)量，即可得出相應(yīng)權(quán)重，簡單便捷。本框架資源層共分10個(gè)資源即：中心城區(qū)、閔行區(qū)、寶山區(qū)、嘉定區(qū)、浦東新區(qū)、金山區(qū)、松江區(qū)、青浦區(qū)、奉賢區(qū)和崇明區(qū)。利用上述方式，即每區(qū)域權(quán)重為1/10，但與上海實(shí)際情況不符，故在其基礎(chǔ)上根據(jù)每個(gè)區(qū)域面積大小并結(jié)合《上海市地面沉降控制區(qū)范圍劃定方案》（上海市規(guī)劃和國土資源管理局、上海市水務(wù)局、上海市住房和城鄉(xiāng)建設(shè)管理委員會(huì)、上海市交通委員會(huì)，關(guān)于印發(fā)《上海市地面沉降控制區(qū)范圍劃定方案》的通知，滬規(guī)土資礦[2018]155號(hào)）中每個(gè)區(qū)域地面沉降防治重要程度綜合確定各區(qū)域權(quán)重。

地下水位、地下水開采和地面沉降三大指標(biāo)很難判斷哪個(gè)更重要，故采取上述公式平均分配權(quán)重。

（3）計(jì)算地下水超采系數(shù)

收集真實(shí)數(shù)據(jù)并計(jì)算各指標(biāo)值，最終計(jì)算地下水超采系數(shù)。系數(shù)的取值范圍一般為0～1，當(dāng)一個(gè)系數(shù)為0時(shí)，即表示某指標(biāo)下某資源區(qū)域地下水超采非常嚴(yán)重，為1則表示地下水超采控制的非常好，故系數(shù)值越大表示地下水超采控制的越好。最終的地下水超采系數(shù)用Ω表示，相應(yīng)的各資源、各指標(biāo)的系數(shù)為ΩA、ΩI。數(shù)據(jù)收集計(jì)算出的指標(biāo)值（x）需無量綱化且取值范圍為0～1，ΩI為x的函數(shù)。本框架計(jì)算地下水超采系數(shù)，采用線性方程，各系數(shù)計(jì)算公式如下：

本次在參照Bui等人方法[13]的基礎(chǔ)上，將地下水超采系數(shù)進(jìn)行分級(jí)，如表4所示。在傳統(tǒng)的嚴(yán)重超采區(qū)、一般超采區(qū)和非超采區(qū)之外增加了超采臨界區(qū)和開采潛力區(qū)，提高了地下水超采區(qū)劃定的精細(xì)度。

表4 地下水超采系數(shù)分級(jí)表Table 4 Scale of groundwater overexploitation Indices

（4）地下水超采綜合評(píng)估

本次三大指標(biāo)為地下水位、地下水開采和地面沉降，結(jié)合上海實(shí)際含水層情況，其中地下水位細(xì)分為潛水、第一至第五承壓含水層等6個(gè)水位變化指標(biāo)（崇明區(qū)一含缺失，只有5個(gè)水位變化指標(biāo)），地下水開采細(xì)分為第二～第五承壓含水層等4個(gè)開采系數(shù)指標(biāo)，加上1個(gè)地面沉降指標(biāo)，故共有11個(gè)小指標(biāo)（崇明區(qū)10個(gè)）。為使每個(gè)指標(biāo)與超采系數(shù)呈正相關(guān)，定義每個(gè)含水層水位變化指標(biāo)為1減去水位下降區(qū)域面積占比，每個(gè)含水層開采系數(shù)指標(biāo)為1減去開采系數(shù)（當(dāng)開采系數(shù)大于1時(shí)，取該指標(biāo)值為0），地面沉降指標(biāo)為1減去地面沉降區(qū)域面積占比。為了便于對(duì)比，指標(biāo)計(jì)算所采用數(shù)據(jù)仍參照前述傳統(tǒng)方法，以2016年為基準(zhǔn)年，評(píng)價(jià)期為2001～2016年。最終得出上海市地下水超采綜合評(píng)估結(jié)果如表5所示。

表5 上海市地下水超采區(qū)綜合劃定Table 5 Comprehensive delimitation of groundwater overexploitation zones in Shanghai

從上表可知，在綜合考慮地下水位、地下水開采及地面沉降的基礎(chǔ)上，上海市最終地下水超采系數(shù)為0.656，屬于非超采區(qū)，其結(jié)果與近年上海市壓采地下水、持續(xù)開展人工回灌，各含水層地下水位抬升且地面沉降總體趨緩可控的現(xiàn)狀較吻合。每個(gè)行政區(qū)除了金山區(qū)外都屬于非超采區(qū)，金山區(qū)雖然地下水幾乎沒有開采，但可能由于受鄰近省份歷史地下水開采影響，部分區(qū)域存在地下水位下降，故超采系數(shù)值稍低，為0.585，屬于超采臨界區(qū)。

4 結(jié)論和建議

本文在傳統(tǒng)的水位動(dòng)態(tài)法、開采系數(shù)法和引發(fā)問題法對(duì)比分析的基礎(chǔ)上，基于層次分析法建立上海市地下水超采區(qū)綜合劃定框架，統(tǒng)籌考慮地下水位、地下水開采及地面沉降，綜合劃定了上海市地下水超采區(qū)。

（1）水位動(dòng)態(tài)法、開采系數(shù)法可以比較好的反映評(píng)價(jià)區(qū)域當(dāng)下的超采情況；由于地面沉降的滯后性、不可逆性，引發(fā)問題法更反映評(píng)價(jià)區(qū)域歷史的超采情況；而本次地下水超采綜合劃定框架，統(tǒng)籌考慮了三方面因素，為后續(xù)深入研究提供了參考和借鑒。

（2）本次劃定框架資源層權(quán)重分配，綜合考慮了上海行政區(qū)域面積和各區(qū)地面沉降防治重要程度的影響，較符合上海實(shí)際情況。

（3）本次僅建立了初步的劃定框架，指標(biāo)層權(quán)重為平均分配，且地下水超采系數(shù)計(jì)算采用線性函數(shù)，后續(xù)可結(jié)合水資源管理標(biāo)準(zhǔn)、相關(guān)政策和專家意見重新調(diào)整權(quán)重，并根據(jù)監(jiān)測數(shù)據(jù)進(jìn)行計(jì)算方程擬合，建立標(biāo)準(zhǔn)的地下水超采區(qū)綜合劃定體系，使其更符合上海實(shí)際情況及政府需求。

（4）目前上海在改善地下水超采程度和控制地面沉降等方面都已取得顯著成效，但地下水位持續(xù)抬升，可能帶來一定不良效應(yīng)，需因地制宜，對(duì)地下水實(shí)行更精細(xì)化、科學(xué)化管控，加強(qiáng)地下水采灌動(dòng)態(tài)調(diào)節(jié)。