祝悅 劉家福 申琳 譚曉宇 白淑夏

(吉林師范大學(xué)旅游與地理科學(xué)學(xué)院，吉林 四平 136000)

地下水屬于自然資源中水資源重要組成成分，憑借其特有的水質(zhì)良好、水蘊(yùn)藏量豐富等特點(diǎn)，使其成為維系良好生態(tài)環(huán)境資源的要素，也是支撐我國(guó)現(xiàn)代化建設(shè)及經(jīng)濟(jì)社會(huì)高度發(fā)展的重要戰(zhàn)略物資。地下水在我國(guó)水環(huán)境中起到重要作用，特別是在干旱地區(qū)，地下水因?yàn)槠渌|(zhì)好、水量大和水溫穩(wěn)定成為很多地方的主要用水來(lái)源，地下水資源的重要性不容忽視，地下水還起到平衡地質(zhì)和生態(tài)環(huán)境的作用[1]。城市化進(jìn)程的不斷加快及經(jīng)濟(jì)水平的提升，也導(dǎo)致我國(guó)地下水資源出現(xiàn)了各種各樣的問(wèn)題，近年來(lái)地下水資源開(kāi)采利用過(guò)度，水質(zhì)被污染現(xiàn)象頻繁發(fā)生，地下水的問(wèn)題制約我國(guó)經(jīng)濟(jì)可持續(xù)發(fā)展[2]。掌握地下水資源富集性狀況，對(duì)于地下水的合理開(kāi)發(fā)使用有著重要的主導(dǎo)作用。傳統(tǒng)的地下水資源監(jiān)測(cè)方法需要消耗過(guò)多的時(shí)間與人力財(cái)力，所獲得的結(jié)果差強(qiáng)人意，并不理想，并且不適用于長(zhǎng)期觀測(cè)或?qū)崟r(shí)監(jiān)測(cè)，為克服傳統(tǒng)手段帶來(lái)的不便利性，需要通過(guò)其它的調(diào)查方法進(jìn)行研究[3]。

由于地下水的特殊性，很難用常規(guī)手段進(jìn)行富集性評(píng)估。借助遙感技術(shù)大范圍、周期短以及包含信息豐富等優(yōu)點(diǎn)，配合地理信息系統(tǒng)的空間分析處理能力，為進(jìn)行淺層地下水資源的富集性評(píng)估提供有效的解決方法。地理信息系統(tǒng)是一門高效提取并分析空間信息的工具，可以對(duì)獲得的海量信息進(jìn)行存儲(chǔ)、管理、分析及應(yīng)用[4]。遙感，即通過(guò)探測(cè)波進(jìn)行非接觸、遠(yuǎn)距離的一門新興的探測(cè)技術(shù)，憑借其探測(cè)距離遠(yuǎn)、區(qū)域廣、可以全天候多時(shí)相對(duì)目標(biāo)地進(jìn)行觀測(cè)的特點(diǎn)，近些年被廣泛應(yīng)用[5]。其中通過(guò)遙感與GIS技術(shù)相結(jié)合在淺層地下水的監(jiān)測(cè)和研究中發(fā)展迅速，許多專家學(xué)者建立不同的地下水富集性的遙感評(píng)估模型，為利用遙感與GIS技術(shù)解決地下水富集性評(píng)估問(wèn)題提供了前期理論支撐[6]。M Vasanthavigar[7]等提取地形、地貌、水網(wǎng)密度和斷裂密度以上4個(gè)指標(biāo)對(duì)盆地內(nèi)的裂隙水進(jìn)行研究。Gumma[8]等選取7個(gè)指標(biāo)進(jìn)行地下水富集性研究，結(jié)果表明，可以利用遙感和GIS技術(shù)進(jìn)行國(guó)家范圍的地下水富集性研究。許顥礫[9]等選取與地下水富集性相關(guān)的6個(gè)指標(biāo)對(duì)外伶仃島的地下水富集性進(jìn)行評(píng)估。許春華[10]等用GRSFAI方法對(duì)阿里地區(qū)進(jìn)行地下水富集性研究，證明該方法有效和適用地下水評(píng)估。本文利用遙感和GIS技術(shù)進(jìn)行數(shù)據(jù)獲取與處理，并采用GRSFAI法對(duì)吉林省淺層地下水富集性進(jìn)行評(píng)估。

1 研究區(qū)概況

吉林省位于我國(guó)東北地區(qū)中部，地處E121°38′～131°19′，N40°50′～46°19′。其轄區(qū)內(nèi)有9個(gè)地級(jí)行政區(qū)，人口眾多，總面積為18.74×104km2，約占全國(guó)總面積的2%，如圖1所示。吉林省上連黑省，下接遼城，夏季炎熱多雨，冬季白雪皚皚，全年四季分明，屬溫帶大陸性季風(fēng)氣候，全年平均溫度-2～6℃。由于自然環(huán)境影響，自東南至西北逐漸由濕轉(zhuǎn)干。年平均降水為400～600mm，受季節(jié)及區(qū)域影響較大。地面起伏自西北向東南呈明顯階梯狀變化，以中部大黑山為主要分界線，將研究區(qū)劃分為山丘區(qū)和平原區(qū)2部分。地下水資源受自然環(huán)境的限制，研究區(qū)內(nèi)地下水資源分布存在差異，地下水富集性也存在不同[11]。

圖1 研究區(qū)概況

2 數(shù)據(jù)處理與方法

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

從USGS(https://earthexplorer.usgs.gov/)上獲取2018年吉林省的Landsat 8 OLI多光譜遙感影像，此數(shù)據(jù)分辨率為30m；從NASA(https://modis.gsfc. nasa.gov/)上獲取的MODIS13Q1數(shù)據(jù)產(chǎn)品，此數(shù)據(jù)分辨率為250m；于地理空間數(shù)據(jù)云(https://www. gscloud.cn)中下載SRTM DEM數(shù)據(jù)，此數(shù)據(jù)分辨率為90m。對(duì)所下載數(shù)據(jù)進(jìn)行篩選，排除成像色彩相差大、積雪和云覆蓋過(guò)于嚴(yán)重的數(shù)據(jù)，選出云量較少且質(zhì)量較好的影像，通過(guò)對(duì)獲取的數(shù)據(jù)進(jìn)行了幾何校正，使所有的底圖之間在地理坐標(biāo)上都能相互匹配，再對(duì)數(shù)據(jù)進(jìn)行輻射校正和大氣校正，消除大氣散射和光照等因素對(duì)遙感影像的干擾，利用重采樣將所有專題圖的空間分辨率都統(tǒng)一為30m[12]。

地下水富集性評(píng)估涉及眾多指標(biāo)因素，指標(biāo)各個(gè)方面均不相同，則要對(duì)所選的指標(biāo)做歸一化處理之后才能進(jìn)行綜合加權(quán)[13]。歸一化的公式：

(1)

式中，Xi為指標(biāo)數(shù)據(jù)的統(tǒng)計(jì)值；Xmax為指標(biāo)數(shù)據(jù)值中的最大值；Xmin為指標(biāo)數(shù)據(jù)值中的最小值。

2.2 研究方法

GRSFAI法(Groundwater Remote Sensing Fuzzy Assessment Index)即地下水富集性遙感評(píng)估方法[14]，其值范圍為0～1，0代表吉林省地下水富集性最差，1代表吉林省地下水富集性最好。各個(gè)指標(biāo)對(duì)地下水富集性的影響因素均不相同，每個(gè)指標(biāo)分別對(duì)應(yīng)一個(gè)權(quán)重Wi。利用ArcGIS的空間分析工具模塊進(jìn)行加權(quán)計(jì)算，公式：

(2)

式中，Wi為各指標(biāo)的權(quán)重值；Yi為歸一化后的各指標(biāo)專題圖。

3 指標(biāo)選取

地下水富集性由多個(gè)因素決定，選取的6個(gè)指標(biāo)需要進(jìn)行歸一化處理后，利用ArcGIS進(jìn)行綜合加權(quán)計(jì)算，地表起伏度、坡度、河網(wǎng)密度、植被覆蓋度、地表溫度和土壤濕度，吉林省淺層地下水富集性的具體影響情況如表1所示[15-20]。

表1 評(píng)估指標(biāo)

4 權(quán)重確定

選取6種不同指標(biāo)因子構(gòu)建吉林省地下水富集性評(píng)估模型，由于這些指標(biāo)對(duì)地下水富集性產(chǎn)生的影響是不同的。為了使實(shí)驗(yàn)結(jié)果真實(shí)可靠，在充分利用這些指標(biāo)的基礎(chǔ)上[21]，結(jié)合分析者和專家的經(jīng)驗(yàn)知識(shí)，采用層次分析法[22]確定指標(biāo)權(quán)重，邀請(qǐng)具有豐富相關(guān)經(jīng)驗(yàn)知識(shí)的14名專家對(duì)6種指標(biāo)進(jìn)行評(píng)價(jià)打分，以量化衡定影響因子對(duì)地下水富集性的重要程度，將打分結(jié)果構(gòu)成判斷矩陣。在進(jìn)行權(quán)重計(jì)算之前，對(duì)判斷矩陣進(jìn)行一致性檢驗(yàn)，若計(jì)算結(jié)果未通過(guò)一致性檢驗(yàn)，則說(shuō)明所選取的權(quán)重值誤差過(guò)大，對(duì)指標(biāo)做出的權(quán)重值評(píng)價(jià)結(jié)果不可靠，需要重新做權(quán)重值選定直到滿足一致性檢驗(yàn)為止，經(jīng)計(jì)算一致性檢驗(yàn)結(jié)果為0.0136，符合一致性檢驗(yàn)，根據(jù)其原理和方法得到指標(biāo)權(quán)重為Wi，如表2所示。

表2 指標(biāo)權(quán)重值

5 結(jié)果分析

經(jīng)過(guò)歸一化處理的各個(gè)專題圖與計(jì)算得到的各個(gè)指標(biāo)的權(quán)重值，根據(jù)公式(2)計(jì)算得出吉林省地下水富集性情況。

通過(guò)野外實(shí)地測(cè)量數(shù)據(jù)，選取36個(gè)地下水位觀測(cè)點(diǎn)與獲得富集性評(píng)估結(jié)果進(jìn)行相關(guān)性分析[23]，驗(yàn)證結(jié)果表明，兩者之間呈指數(shù)相關(guān)，其相關(guān)系數(shù)R2為0.8011，表明可以用GRSFAI法對(duì)吉林省淺層地下水富集性進(jìn)行有效評(píng)估。相關(guān)性分析結(jié)果如圖2所示。

圖2 相關(guān)性分析結(jié)果圖

根據(jù)吉林省地下水富集性評(píng)估結(jié)果與實(shí)測(cè)數(shù)據(jù)相關(guān)性分析結(jié)果，將地下水富集性評(píng)估結(jié)果分為5個(gè)等級(jí)，具體如表3所示。再根據(jù)上述的5個(gè)等級(jí)利用ArcGIS分析得到吉林省富集性表，各個(gè)等級(jí)在各個(gè)市中的面積情況如表4所示。

由表4可以看出，吉林省地下水富集性評(píng)估結(jié)果為優(yōu)秀的區(qū)域面積為3386.43km2，占吉林省總面積的1.81%，其中松原市的優(yōu)秀區(qū)域的面積為736.48km2，是優(yōu)秀等級(jí)中面積最大的，長(zhǎng)春市次之，為500.30km2；吉林省地下水富集性評(píng)估結(jié)果為良好的區(qū)域面積為2765.01km2，占吉林省總面積的1.47%，其中四平市的良好區(qū)域面積為206.91km2，是良好等級(jí)中面積最小的；吉林省地下水富集性為一般的區(qū)域面積為2640.37km2，占吉林省總面積的1.41%，富集性為一般的區(qū)域在吉林省分布基本較為平均；富集性評(píng)估結(jié)果為及格的區(qū)域面積為46147.27km2，占吉林省總面積的24.61%，其中長(zhǎng)春市的及格區(qū)域的面積為13346.66km2，是及格等級(jí)中面積最大的，白城市次之，為9331.50km2；吉林省大部分地區(qū)的地下水富集性級(jí)別為很差，富集性評(píng)估結(jié)果為很差的區(qū)域面積為132543.48km2，占吉林省總面積的70.70%，證明吉林省的地下水情況不容樂(lè)觀，其中白城市的很差區(qū)域的面積為15535.33km2。

表3 吉林省地下水富集性評(píng)估數(shù)值表

表4 吉林省富集性表

6 結(jié)論

通過(guò)選取與地下水賦存相關(guān)的地表起伏度、坡度、河網(wǎng)密度、植被覆蓋度、地表溫度和土壤濕度6個(gè)指標(biāo)，建立吉林省地下水富集性評(píng)估模型，得到吉林省富集性表，富集性評(píng)估結(jié)果為很差、及格、一般、良好和優(yōu)秀的區(qū)域面積分別占吉林省總面積的70.70%、24.61%、1.41%、1.47%和1.81%，其中富集性評(píng)估結(jié)果在及格和及格以上的區(qū)域面積總和占吉林省總面積的29.30%。

通過(guò)富集性評(píng)估結(jié)果與實(shí)地測(cè)量36個(gè)水位觀測(cè)點(diǎn)數(shù)據(jù)進(jìn)行相關(guān)性分析，系數(shù)為0.8011，表明GRSFAI法對(duì)地下水富集性評(píng)估有效，能夠節(jié)約大量時(shí)間和成本，分析得出的吉林省地下水富集性分布表可以為今后地下水的管理工作提供一定參考。