李恒超,鄧英爾,2,史振環(huán),宋 剛,鐘貴莉
(1.成都理工大學(xué) 環(huán)境與土木工程學(xué)院,四川 成都 610000;2.成都理工大學(xué) 地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護國家重點實驗室,四川 成都 610059)
自我國20世紀70年代末期大規(guī)模經(jīng)濟建設(shè)以后,經(jīng)濟的快速發(fā)展,一并帶來了水資源緊缺、污水排放不達標生態(tài)系統(tǒng)破壞等社會性問題,迫使人們不得不認真考慮水資源的承載能力問題。水資源承載力是在此之上引申出來的一個概念,強調(diào)水資源的最大開發(fā)規(guī)?;蛘咚Y源對經(jīng)濟社會發(fā)展的支撐能力[1]。隨后,我國學(xué)者對此展開了一系列的研究,到21世紀初已經(jīng)形成了以模糊綜合評價法,系統(tǒng)動力學(xué)法,多目標模型分析法為代表相對完善的應(yīng)用體系。其中模糊綜合評價法是應(yīng)用最多的一種之一。該方法是基于模糊數(shù)學(xué)的綜合評標方法。該綜合評價法根據(jù)模糊數(shù)學(xué)的隸屬度理論把定性評價轉(zhuǎn)化為定量評價,即用模糊數(shù)學(xué)對受到多種因素制約的事物或?qū)ο笞龀鲆粋€總體的評價。它具有結(jié)果清晰,系統(tǒng)性強的特點,能較好地解決模糊的、難以量化的問題[2]。
貢力等[3]應(yīng)用模糊綜合評判方法對蘭州市水資源承載力進行了分析評價,研究了蘭州市未來政策施后水資源承載力的動態(tài)變化過程。張忠學(xué)[2]等運用模糊綜合評價法對綏化市北林區(qū)1996年到2012年的農(nóng)業(yè)水資源進行承載力研究;高俊[4]從主體支撐力、調(diào)控力、客體壓力3個角度選擇15項典型指標構(gòu)建水資源承載力評價體系,科學(xué)評價了呼和浩特市2015-2019年的水資源承載狀況。吳凡等[5]為定量評價新疆各地州水資源承載力狀況,在構(gòu)建了模糊集對評價模型基礎(chǔ)之上,結(jié)合五元減法集對勢對新疆15個地州水資源承載力進行評價并診斷其發(fā)展態(tài)勢。馮書倉等[6]側(cè)重于持續(xù)發(fā)展質(zhì)量評價體系,選用水資源條件、供水需水系統(tǒng)、社會經(jīng)濟發(fā)展水平和生態(tài)環(huán)境狀況四個一級指標和22個二級指標,建立了基于模糊評價法的水資源可持續(xù)發(fā)展質(zhì)量評價體系。國外學(xué)者直接以水資源承載力為主題的專項研究幾乎沒有,而是在探討水資源與生態(tài)承載力及人口承載力的歷程中提出了一些相似內(nèi)涵的概念,如礦區(qū)承載力(MAWRCC)[7],水資源緊缺指數(shù)(water stress index)[8]等。
本次研究區(qū)位于河北省張家口市壩上地區(qū),目前相關(guān)該地區(qū)的研究大多停留在大范圍上,對于該地區(qū)的精細化研究尚顯缺乏,本文將根據(jù)相關(guān)資料,根據(jù)模糊判別法建立水資源承載力評價體系,對該地區(qū)的水資源現(xiàn)狀及開發(fā)利用情況、當(dāng)前社會經(jīng)濟條件下的水資源人口承載力進行探討,并對未來的水資源承載狀況做出預(yù)測,旨在對該地區(qū)未來水資源合理利用和布置及產(chǎn)業(yè)調(diào)整提供一定的理論和數(shù)據(jù)支撐。
本區(qū)地處河北西北部,北京市西北方向,位于內(nèi)蒙古高原的“南緣”,俗稱“壩上”。地理坐標東經(jīng)113°50′~116°04′,北緯40°44′~41°32′,占張家口市總面積的37.4%。包括張北縣、康??h、沽源縣和尚義縣北部壩上部分及察北、塞北兩個管理區(qū),總面積12 863 km2??傆^其貌,北部及東南部為低山丘陵,最高峰冰山梁位于沽源縣東南部,海拔標高2 211 m。壩上地區(qū)氣候明顯受內(nèi)蒙古高壓中心影響,夏季涼爽短促,雨量集中,冬季寒冷漫長、干燥;春秋兩季多風(fēng)、稍降雨,氣溫變化劇烈。氣候總體特征是干旱、少降水、多風(fēng)沙,屬于典型的寒溫帶半干旱大陸季風(fēng)氣候。
壩上地區(qū)以內(nèi)陸河水系為主,無境外水量補給,水資源的形成主要靠流域內(nèi)有限的大氣降水補給,形成途徑單一。據(jù)多年降雨量變化分析,各縣的降水量趨勢性變化不明顯。八十年代尚義縣、沽源縣降水量明顯偏枯,其他縣與多年平均值接近;九十年代至本世紀,降水量接近平均值,但近幾年豐枯交替更加頻繁。綜合來看,該地區(qū)多年平均降水量380 mm,降水量年際變化大,且由東南向西北逐漸減少;多年平均蒸發(fā)量1 655 mm,是多年平均降水量的5倍。區(qū)內(nèi)河流較為發(fā)育,沽源縣內(nèi)灤河、潮白河、內(nèi)陸河三大水系貫穿四方,張北縣內(nèi)共有25條河流,總長度793 km,尚義縣為內(nèi)陸河流域,主要有大青溝和、二龍河,流域面積247 km2,唯康??h境內(nèi)僅分布十幾條季節(jié)性小河流。1994年,地表水供水量1 615萬 m3,地下水供水量4 760萬 m3,今年內(nèi)由于壩上地區(qū)降水徑流減少,地表供水工程老化失修,供水保證率低,對地下水的開采量逐年上升,其在整個區(qū)域內(nèi)水資源利用量的比例也不斷提高。
2000年以來,由于壩上地區(qū)大力發(fā)展農(nóng)業(yè),以發(fā)展錯季蔬菜為支柱產(chǎn)業(yè)的水澆地面積迅速擴大,地下水開采量逐年加大,機井?dāng)?shù)量猛增,地下水位連年下降,深層水的開采利用率也逐年提高。全區(qū)現(xiàn)存數(shù)十座中小型水庫,機井?dāng)?shù)量已接近6萬眼。
該區(qū)地下水資源采補失調(diào),在水資源利用效率上仍有較大提升空間,主要原因是農(nóng)業(yè)需求的擴大,機井群遍布全區(qū),在大片蔬菜種植區(qū)更是出現(xiàn)扎堆打井的現(xiàn)象,機井布局沒有合理規(guī)劃,再加上大部分節(jié)水工程年久失修或超過設(shè)計使用年限失去使用功能,灌溉效率低,水資源浪費較為嚴重。另一方面由于多數(shù)農(nóng)供機沒有采取止水措施,造成多個含水層的水力連通,致使深、淺層水混用和混采,水質(zhì)和生態(tài)環(huán)境問題逐步凸顯。
設(shè)評判指標集合為U={u1,u2,u3…,um},評語集合為V={v1,v2,v3…,vn},u、v代表的是具體的變量意義,其中U代表的是綜合評判因素的集合,評判因素的選取因人而異,大多是根據(jù)研究區(qū)水資源和社會經(jīng)濟特定選定;V代表的是評語組成的集合。首先對評判指標集合中的某一項指標ui進行單因素評判,確定其所對應(yīng)評價vj的隸屬度rij,這樣就得到第i個因素的單因素評判集合ri={ri1,ri2,ri3},m個評判指標的評價集構(gòu)造出評價矩陣R。
(1)
(2)
式中a為綜合評判結(jié)果矩陣B的水資源承載力的綜合評分值,a值越高,代表該區(qū)域水資源承載力越大。為了突出占優(yōu)勢等級的作用,進行區(qū)域水資源承載力的綜合評價與分析K值,本次研究區(qū)為干旱區(qū),結(jié)合前人研究[11],K取1。
與水資源承載力相關(guān)聯(lián)的評價指標很多,為了研究的準確性和針對性,在查閱相關(guān)文獻的基礎(chǔ)上[2,12],結(jié)合研究區(qū)自身特點和經(jīng)濟發(fā)展情況,此次研究選取水資源開發(fā)率(u1):75%頻率可供水總量與水資源總量之比,反映水資源現(xiàn)階段開發(fā)程度;水資源利用率(u2):75%頻率可供水總量與可利用水資源總量之比;農(nóng)業(yè)水資源利用率(u3):農(nóng)業(yè)水資源用量與可利用水資源用量之比,反映水資源對農(nóng)業(yè)的支持程度;人均水量(u4):75%頻率可供水總量與當(dāng)?shù)乜側(cè)丝谥?;生活需水定額(u5);耕地灌溉率(u6)6個因素作為評價指標。各評價指標及其分級見表1。
此次將評價指標分為3級。其中V1級表示水資源承載力接近飽和值,進一步開發(fā)發(fā)展?jié)摿π?,?yīng)及時調(diào)整農(nóng)牧業(yè)用水,提高水資源利用率,減少不必要浪費。V3級表示水資源承載力良好,后續(xù)開發(fā)潛力大,可以適當(dāng)加大灌溉規(guī)模,增加農(nóng)業(yè)投入。V2級處于兩者之間,表示水資源開發(fā)已達到一定規(guī)模,但仍舊擁有一定的開發(fā)潛力,可酌情適當(dāng)開發(fā)。為了更定量反映研究區(qū)水資源承載力情況,對上述3個評語等級進行1分制評分,即a1=0.05,a2=0.5,a3=0.95,這樣就將定性S描述水資源承載力的評判因素轉(zhuǎn)換成可以定量評價該地區(qū)的水資源情況的數(shù)值,數(shù)值越大表示水 資源承載力越強,后續(xù)可開發(fā)利用性更強。
表1 各評價指標及其分級
根據(jù)上述6個評價因素的實際值,對照其分級指標來推求評判矩陣的隸屬函數(shù)Uvi(ui)的值。操作過程中采用模糊化處理并構(gòu)建隸屬函數(shù),使其值在各級間平滑過渡。對于評價因素 u1、u2、u3、u5、u6,各評語級隸屬函數(shù)計算公式如下:
對于評價因素u4,其評語及隸屬度計算公式如為:
上式中:評價等級V1與V2的臨界值為k1;評價等級V2與V3的臨界值為k3;評價等級V2的中點值為k2,即k2=(k1+k3)/2。
根據(jù)張家口市壩上地區(qū)社會經(jīng)濟條件和水資源利用狀態(tài),分別統(tǒng)計了壩上地區(qū)1993年、2003年、2013年6個評價指標。2023年為預(yù)測數(shù)據(jù),是根據(jù)該地區(qū)用水現(xiàn)狀并結(jié)合有關(guān)規(guī)劃資料,按供需量進行預(yù)測得到的基礎(chǔ)數(shù)據(jù)計算而成。由此計算得到6個評價因素的數(shù)值,具體見表2。按照上述隸屬度計算公式,計算出評價因素對應(yīng)不同評價等級的隸屬度rij,從而構(gòu)造出整個評價矩陣。
在得到評價矩陣B和權(quán)重矩陣A的基礎(chǔ)上,按照公式B=AR計算得出張家口市水資源承載力評價結(jié)果矩陣B,各年各指標隸屬度矩陣R結(jié)果如下:
表2 壩上地區(qū)評價因素指標統(tǒng)計
綜上,1993年壩上地區(qū)水資源承載力最終評價矩陣B計算過程如下:
表3 張家口市壩上地區(qū)多年水資源承載力綜合評價結(jié)果
從表3可以看出,1993-2023年這40 a間張家口市壩上地區(qū)水資源承載力綜合評分值對V1的隸屬度顯著高于V3和V2的隸屬度,且對V1的隸屬度逐年上升,而對V2及V3的隸屬度則不斷下降,說明該地區(qū)水資源開發(fā)利用率已接近極限,繼續(xù)開發(fā)潛力較小。在不同年份的綜合評分中,承載力評分分值均小于0.5,評價結(jié)果對V1的隸屬度遠遠大于V3的隸屬度,預(yù)測年份2013年對V1的隸屬度竟高達0.581,說明在不斷地社會發(fā)展中,該地區(qū)水資源承載已接近飽和,后期應(yīng)注重調(diào)整農(nóng)牧業(yè)結(jié)構(gòu),提高水資源利用率。
對隸屬度矩陣R的計算結(jié)果分析可以得出,壩上地區(qū)在1993-2023年這40 a間,水資源利用率u2、人均水量u4、生活需水定額u5對V1的隸屬度都在逐年增加,而對V2的隸屬度大體上都有逐年下降的趨勢,說明隨著社會、生活的不斷向前進步,人們對水資源的需求量也翌日劇增 ,這一切都迫使水資源的開發(fā)進入一個快速發(fā)展期,大量的水資源被開發(fā),被利用;水資源開發(fā)率u1、農(nóng)業(yè)水資源利用率u3、耕地灌溉率u6對V3的隸屬度高于對V1和V2的隸屬度,且大體上對V3的隸屬度有不斷下降的走勢,這說明隨著農(nóng)牧業(yè)的發(fā)展,尤其是農(nóng)業(yè)的發(fā)展需求,灌溉面積的加大對水資源的需求量不斷增加,該地區(qū)的水資源開發(fā)利用率正逐步迫近極限,這同樣側(cè)面佐證了上述推論。
(1)通過以上計算,河北省張家口市壩上地區(qū),1993-2013年這30 a間水資源承載力綜合得分數(shù)值一直處于不斷下降的趨勢,從1993年的0.437下降到2013年的0.226,且對V1的隸屬度遠大于對V2和V3的隸屬度,都說明隨著該地區(qū)大力發(fā)展農(nóng)牧業(yè),對水資源的需求量不斷增大,導(dǎo)致地下水超采,水資源承載已接近飽和,后續(xù)開發(fā)潛力較小。
(2)通過2023年預(yù)估數(shù)據(jù)對壩上地區(qū)進行水資源承載力評價,綜合評分數(shù)值為0.188,對V1和V2的隸屬度大致相似,表明若照目前這個用水趨勢,后期有極大可能遭遇水資源供給不足,水資源無法承載當(dāng)?shù)厣a(chǎn)的情形。
(3)應(yīng)采取適當(dāng)措施緩解水資源壓力,從而提高水資源承載力,實現(xiàn)可持續(xù)發(fā)展。筆者認為應(yīng)從用水方式,取水手段兩方面來考慮:
壩上地區(qū)以農(nóng)業(yè)為支柱產(chǎn)業(yè),農(nóng)業(yè)用水量在全年地下水開采量中的占比高達80%以上,故應(yīng)創(chuàng)新農(nóng)業(yè)發(fā)展態(tài)勢,推廣耐旱節(jié)水品種,提倡高效節(jié)水種植和澆灌方式,提高水資源利用率。
大力推進河北省南水北調(diào)工程、引黃工程建設(shè)及其配套水廠及輸送網(wǎng)管工程建設(shè),引入長江、黃河水進行灌溉,從而大幅度緩解本地水資源壓力。