GIS 在地下水污染分析評價中的應用與技巧

黃金玉，姜月華，蘇晶文

(南京地質調查中心，江蘇南京 210016)

GIS 在地下水污染分析評價中的應用與技巧

黃金玉?，姜月華，蘇晶文

(南京地質調查中心，江蘇南京 210016)

以長江三角洲地區地下水為研究對象，對地下水污染進行調查和分析評價，將獲得的野外調查數據和樣品測試數據通過GIS技術進行信息處理，完成地下水檢出與超標評價、單指標質量評價、分類指標質量評價以及單指標污染評價，并將分析評價結果通過GIS圖件形式顯示輸出。本文詳細介紹GIS在地下水污染分析評價中的應用過程和方法，并根據實際情況對GIS使用技巧進行闡述。

地理信息系統；地下水；質量評價；等值線

1 引 言

2 GIS在地下水質量分析評價中應用

水資源是支撐長江三角洲地區經濟社會可持續發展的戰略資源，地下水是水資源的重要組成部分。然而，隨著人類活動的加劇，石油類污染，城市垃圾和生產生活污水的不合理處置以及農業生產中農藥、化肥的大量使用，導致地下水有機污染狀況日趨加重。地下水污染引起水中有害、有毒元素不斷增多，水質下降甚至惡化，給長江三角洲地區飲用水安全帶來嚴重威脅。

將GIS技術應用到地下水質量分析評價中，不僅僅是幾個簡單的操作步驟，而是一個系統而復雜的操作過程，通過不同功能模塊結合不同操作完成。

2.1 GIS操作過程

對地下水污染分析評價，首先需要將野外空間數據進行投影轉換，地下水采樣點通過經緯度值記錄其空間地理位置，讀取對應地理坐標值，通過GIS投影技術轉換為矢量圖元，然后運用GIS參數設定和等值線分析對采樣點(GIS點圖元)屬性值進行分析，過程如圖1所示。

圖1 地下水污染分析評價GIS運用流程圖

2.2 運用GIS技術采樣點矢量化

通過野外獲取各個研究區域的分析點的相關地理信息，如地理坐標(經緯度或公里網值)、地理位置等，通過測試獲得每個分析點的指標值。對于地下水質量分析評價首要操作，則是將分析點進行投影轉換，再根據轉換的區域地理分布結果進一步進行質量評價。運用GIS投影變換功能，GIS在進行投影轉換時需要獲取兩方面的投影參數，一個是用戶投影參數，實際是指野外獲取的分析點的實際地理坐標，在長三角地區獲取的野外坐標單位是經緯度，另外一個參數是結果投影參數，是指作為底圖的矢量電子地圖的地圖參數，根據實際參數值設置完成后，GIS則通過內部運算法則進行轉換，并將投影結果通過GIS點文件投影顯示于底圖中。

2.3 運用GIS進行地下水質量評價

在通過GIS投影轉換的過程中，測試指標值同時轉換為點圖元的屬性值，在評價中就是對這些指標測試數據進行分析，運用GIS分析評價即是對點圖元屬性的分析。

確定參評指標共56項，其中，常規指標25項，非常規指標31項。

(1)常規指標(25項)

無機指標(23項):pH、鐵、錳、銅、鋅、鋁、鈉、氯化物、硫酸鹽、總硬度、溶解性總固體、耗氧量、砷、鎘、六價鉻、鉛、汞、硒、氟化物、硝酸鹽、亞硝酸鹽、氨氮、碘化物；

有機指標(2項):三氯甲烷、四氯化碳。

(2)非常規指標(31項)

揮發性有機物指標(21項):1，1，1-三氯乙烷、三氯乙烯、四氯乙烯、二氯甲烷、1，2-二氯乙烷、1，1，2-三氯乙烷、1，2-二氯丙烷、三溴甲烷、氯乙烯、溴二氯甲烷、一氯二溴甲烷、1，1-二氯乙烯、1，2-二氯乙烯、氯苯、鄰二氯苯、對二氯苯、苯、甲苯、乙苯、二甲苯、苯乙烯；

半揮發性有機物指標(10項):總六六六、γ-BHC (林丹)、總滴滴涕、六氯苯、七氯、苯并(a)芘、萘、蒽、熒蒽、苯并〔b〕熒蒽。

對于單指標質量評價的56項指標分為無機和有機兩大類，需要運用GIS兩類功能，一個是對于有機單指標質量評價通過GIS點圖元參數進行分析區分，在點圖元屬性中存在各點的評價指標測試值，選擇GIS中根據屬性賦參數的操作功能，通過圖元顏色參數對各類水進行分析顯示；另外一個是對于無機單指標質量評價則是通過GIS等值線分析進行，這個分析過程相對較復雜，在MapGIS軟件中進入DTM分析模塊，裝入前面已轉換好的點圖元文件，對此點圖元中需要進行分析的屬性進行點屬性值的提取操作，屬性值提取后進行離散數據網格化得到后綴為Grd網格化文件，針對Grd文件進行平面等值線分析，得到的分析結果通過GIS區域、點以及線文件綜合顯示。

3 GIS應用關鍵步驟

3.1 采樣點空間數據轉換[1～4]

對于地下水野外空間數據點(采樣點)的投影，最關鍵的是投影參數的設置以及屬性的賦值。長江三角洲地區地下水污染評價項目選擇使用MapGIS軟件進行空間數據的轉換，MapGIS讀取的空間數據文件格式是文本格式，將原始數據文件以文本格式存儲后，進入GIS投影變換功能模塊，選擇用戶文件投影菜單，此時系統提供了投影參數設置的交互界面，在此界面中對投影參數進行設置，而投影參數的設置是整個投影操作中較為核心的步驟，需要了解設置的兩個投影參數的具體含義:

用戶投影參數——指野外實際獲得的地理坐標參數，通常是指兩類單位坐標:經緯度和公里網，針對經緯度坐標系類型應選擇地理坐標系，坐標單位是DDDMMSS.SS，對于公里網坐標系單位應選擇大地坐標系，比例尺為1∶1，坐標單位是m；

結果投影參數——指作為投影底圖的地圖參數，通常坐標系類型是投影直角坐標系，區別于用戶投影參數是在結果投影參數中需要輸入更多的參數值，如投影類型(包括有通用橫向墨卡托投影坐標系(UTM)、蘭伯特等角圓錐投影坐標系LAMBERT CONFORMAL CONIC、墨卡托(正軸等角圓柱)投影坐標系MERCATOR、高斯-克呂格(橫切橢圓柱等角))，每個投影類型都對應各自的投影參數值，如高斯-克呂格(橫切橢圓柱等角)對應著投影中心點經度等，蘭伯特等角圓錐投影坐標系對應第一標準緯度、第二標準緯度、中央子午線經度等。

3.2 運用GIS技術進行有機單指標評價

使用GIS參數進行有機單指標評價，通過MapGIS主界面進入輸入編輯功能模塊，打開已完成采樣點(分析點)空間數據轉換的點文件，選擇點屬性編輯菜單中的根據屬性統改參數功能，在屬性條件選擇界面上，選擇對應單指標屬性名稱，完成條件語句”單指標值＞=判定標準值＆＆單指標值＜=判定標準值”，此時GIS通過此判定語句選擇符合條件的點圖元后提供參數設置對話框，將圖元參數中的顏色設置為對應顏色。

3.3 運用GIS技術進行無機單指標評價[7]

相對于有機單指標評價方法，無機單指標評價更復雜，在進行等值線分析中有兩個關鍵參數設置界面比較重要，一個是離散數據網格化中網格參數設置和網格化方法的選擇，在這個界面中網格參數設置開始GIS系統會選擇默認參數值，網格線數X方向的通常默認為101，網格化方法選擇Kriging(克立金)法網格化，并將網格化處理好的文件?.Grd保存；另外一個參數設置界面是等值線參數設置界面，這其中等值層值的設定很重要，根據無機單指標評價標準對起始等值線層值、終止等值線層值以及步長增量(即等值間距)進行選擇設置。

4 運用GIS技術的思考與技巧

通過長江三角洲地區地下水采樣點分布特點發現，該地區采樣點分布不均勻呈不規則形，點間距多樣性，而運用GIS中Kriging等值線分析針對分布均勻有規律等距的分析點做出的等值線更美觀標準，因此，通過多克立金等值線分析參數分析發現，在離散數據網格化時需要根據采樣點分布情況對一些參數進行修訂。在GIS完成等值線分析時會使用其內部運算公式，如X和Y方向的等值線間距就有如下公式:

X方向等之間間距=(X方向最大值-X方向最小值)/X等值線數量

Y方向等之間間距=(Y方向最大值-Y方向最小值)/Y等值線數量

通過公式可以知道X或Y方向間距和等值線數量是互相推算出來了，只要根據實際情況設定其中一個參數值即可，可根據采樣點實際間距將X和Y方向值設定為符合實際分析的數值。

在網格化方法選定后，有兩個重要的參數需要設置，其中最大搜索圓半徑的運算公式如下:

最大搜索圓半徑=X方向等之間間距×Y方向等之間間距

另外一個關鍵參數是變差函數類型中的變程值，這個參數的缺省值為數據點最小包圍矩形對角線的一半長度，但也存在其運算公式，如下:

變程數據值=開根號((X方向最大值-X方向最小值)×(X方向最大值-X方向最小值)+(Y方向最大值-Y方向最小值×(Y方向最大值-Y方向最小值))/2

在實際操作時發現將最大搜索圓半徑調整大小對最后的分析結果有很大的影響，當半徑大時分析結果圖件與實際分析點分布不太符合，而將搜索半徑調整至合適值成圖效果不但美觀準確而且分析區域與分析點分布區域對應。

5 實際完成效果

長江三角洲地區地下水采樣點屬性結構如表1，通過GIS分析結果如圖2所示。

地下水污染屬性結構 表1

圖2 通過GIS參數分析結果(左)和通過GIS等值線分析結果(右)

通過GIS在長江三角洲地區地下水污染調查項目中的應用，總結出應用方法和技巧，也從中明白將GIS技術運用到各個領域都需要根據實際情況對常規應用方法進行改進和分析。

[1] 尹朝陽，尹鵬偉.關于MapGIS中比例變換、格式轉換等相關問題的處理方法與技巧[J].城市勘測，2012，10 (5):112～113.

[2] 陳勇，劉輝，史瑞芝等.數字圖數據到MapGIS數據的格式轉換[J].解放軍測繪學院學報，2004，21(2):154～157.

[3] 李恒民，陳永波，李娜等.基于ArcGISEngine的地圖投影與坐標轉換的研究與實現[J].城市勘測，2012，10(5): 105～108.

[4] 王虎文.WGS-84坐標系與任意坐標系坐標轉換[J].內蒙古石油化工，2004，5(30):18～19.

[5] 李永兵，陳旭瑞，胡俊峰等.基于GIS的地質數據庫系統:研究現狀和發展趨勢[J].地球物理學進展，2002，17 (3):532～539.

[6] 鄧昭明，王俊.GIS技術在區域地下水資源信息系統中的應用研究[J].水文地質工程地質，2004(5):106～108.

[7] 劉明柱，陳鴻漢，葉念軍等.GIS在區域地下水資源評價中的應用[J].水利學報，2002(1):142.

[8] 吳信才.地理信息系統原理與方法[M].北京:電子工業出版社，2002:202～212.

GIS Application in the Analysis and Evaluation of Groundwater Pollution and Skills

Huang Jinyu，Jiang YueHua，Su Jingwen
(Nanjing Center of Geological Survey，Nanjing 210016，China)

With the groundwater in the Yangtze River Delta as the research object，carries on the investigation and the analysis of groundwater pollution，the field investigation data and sample data will get the information processing by GIS technology，complete detection and evaluation of groundwater exceed the standard，the single index evaluation，classification index of quality evaluation and the single pollution index of evaluation，and the evaluation results through the GIS map display output.This paper introduces application of GIS in the process and method of evaluation of groundwater pollution analysis，and according to the actual situation of GIS using the technique described.

GIS；groundwater；quality evaluation；contour

1672-8262(2013)04-161-03

P208.2，P641.69

A

2013—03—01

黃金玉(1982—)，女，助理研究員，主要從事地理信息系統和水文地質環境方面的研究。

中國地質調查局項目(1212010634400)