摘要：文章以孕災地質條件為基礎選取評價指標，採用信息量模型法對全區地質災害進行易發性評價。在易發性評價的基礎上，採用歷史月累計降雨量（大於5年）開展地質災害危險性綜合評價，分別評估人員和基礎設施易損性，設定權重綜合確定承載體易損性；將危險性和易損性評價結果疊加運算，利用ArcGIS軟體的自然間斷點分級方法，根據單元信息量累積頻度分布曲線特徵，將本區的地質災害易發性劃分為極高、高、中、低四個等級，形成風險評價與區劃結果。為支撐國土空間規劃、地質災害防治規劃及重大工程選址等相關工作提供科學依據。
關鍵詞：GIS；信息量模型法；自然間斷點分級法
地質災害包括地裂縫、地面沉降、崩塌、滑坡、地震、火山、泥石流等，是指人為或者自然因素的作用下形成的對人類的生命和財產、環境造成損失和破壞的地質作用現象。我國是地質災害嚴重的國家之一，每年的地質災害都威脅著人民的生命和財產安全，也制約著國民經濟的發展，基於GIS中信息量模型法和自然間斷點分級法對地質災害易發性和危險性分區進行評價和分析不僅符合黃土高原地區地質災害點的實際情況，對模型參數也有很好的擬合。
1GIS的概念
GIS是一種新型的技術，它是在計算機硬、軟體系統支持下，對整個或部分地球表層（包括大氣層）空間中的有關地理分布數據進行採集、儲存、管理、運算、分析、顯示和描述的技術系統。它可以按照空間位置或地理坐標對空間數據進行處理、研究空間的實體、對數據進行有效的管理。GIS是一門信息產業，廣泛應用於各個行業。世界上GIS軟體就有四百多種，它們風格不同，大小也不一樣。在國外有ARCGIS、MGE、GENAMAP；國內有CITYSTAR、MAPGIS和Geostar等。
2GIS在國內外的應用現狀
雖然有很多GIS軟體，但是大體上分為在GIS的基礎上，應用函數庫二次開發出特有的地理信息系統軟體；應用GIS系統來處理用戶數據這兩種情況。現今已經應用到設施管理、資源管理、城市和區域的規劃、教育、石油和天然氣自動製圖等方面。目前GIS已經成功應用在政府管理、資源保護、環境保護、城市規劃建設等很多領域。隨著我國經濟的發展，GIS在測繪、交通運輸、軍事等領域發揮著舉足輕重的作用，GIS的應用主要包括綜合分析評價與模擬預測；在地理空間數據管理中的應用；建立專題信息系統和區域信息系統；在地圖製圖中的應用等。
3GIS信息量模型法在地質災害風險評價中的應用
3.1地質災害易發性評價方法確定
地質災害易發性定量評價方法較多，主要模型有信息量模型、專家系統模型、灰色系統模型、非線性模型及模式識別模型等。考慮到各模型的可行性、數據的易獲取性，本次沁源縣地質災害易發性評價採用信息量模型法。信息量模型反映了一定地質環境下最易致災因素及其細分區間的組合；具體是通過特定評價單元內某種因素作用下地質災害發生頻率與區域地質災害發生頻率相比較實現的。對應某種因素特定狀態下的地質災害信息量公式可表示為：式中：I——對應因素A、j狀態（或區間）下地質災Aj→B害B發生的信息量；N——對應因素A、j狀態（或區間）下地質災害分布的j單元數；N——調查區已知有地質災害分布的單元總數；S——因素A、j狀態（或區間）分布的單元數；S為調j查區單元總數。當I>0時，反映了對應因素A、j狀態（或區間）下地Aj→B質災害發生傾向的信息量較大，地質災害發生的可能性較大，或者說利於地質災害發生；當I<0時，表明因素Aj→BA、j狀態（或區間）條件下，不利於地質災害發生；當I=0時，表明因素A、j狀態（或區間）不提供有關地質災Aj→B害發生與否的任何信息，即因素A、j狀態（或區間）可以剔除掉，排除其作為地質災害預測因子。
3.2評價單元及評價指標選取
本次該縣一般調查區採用10m×10m大小的柵格單元進行，經統計共劃分為254.9萬個評價單元。本次地質災害易發性評價側重的是崩塌、滑坡等地質災害發育的數量多少及其活躍程度。評價指標體系主要為地質災害形成條件，包括地質環境條件（地貌、坡度、坡高、工程地質岩組、河流、植被）和地質災害誘發因素（人類工程活動）兩大類。
3.3各單因素評價指標信息量的確定
本次信息量模型中求出每一因子圖層中各類型地質災害的分布密度。為統計地質災害在每種因素的各個類別中的分布情況，需要利用空間分析功能來實現，在ArcGIS平台下計算各因子圖層各類別的信息量值。利用空間分析模塊中的區域統計功能，將地質災害點文件導入各單因素柵格文件後，得到在相應區域內地質災害點數目，用以上計算結果及已獲調查數據，得到地質災害信息量值，最終得到各評價指標的信息量表，見表1。
4信息量指標權重分配
根據各單因素信息量計算結果，將各種狀態的信息量值按大小排序，較前的狀態分別為：1）在距離道路30-60m區間內；2）坡高>30m；3）坡度>60°；4）距離居民地<10m區間內；5）距離河流在100-200m範圍內；6）植被指數0-0.4（覆蓋率低）；7）粘性土岩土體類型；8）梁峁狀黃土丘陵類型；9）坡型>0.5（權重分配見表2）。上述單因素分布區間屬本次一般調查區地質災害易發的區域，這些因素對崩塌等地質災害的形成發育起決定性的作用，是地質災害形成的主要控制條件。
5應用自然間斷點分級法評價易發分區
利用GIS柵格計算功能對各評價因子進行柵格疊加運算，得到各評價單元信息權值結果，數值越大反映各因素對地質災害發生的貢獻越大，發生地質災害的可能性就越大。經過ArcGIS軟體計算得到本次地質災害易發性評價信息量值區間為[-0.5937，1.7608]，並根據信息量值結果，利用ArcGIS軟體的自然間斷點分級方法，根據單元信息量累積頻度分布曲線特徵和曲線明顯分布的拐點，將本區的地質災害易發性劃分為極高、高、中、低四個等級。其中極高易發區信息量值區間為[1.7608，0.2154]，高易發區信息量值區間為[0.2154，0.0255]，中易發區信息量值區間為[0.0255，-0.2330]，低易發區信息量值區間為[-0.2330，-0.5937]。
5.1地質災害危險性評價
本次一般調查區地質災害危險性評價是在地質災害易發性評價基礎上，採用歷史月累計降雨量（大於5年）進行評價。本次調查收集沁源縣近40年（1980—2019年）歷史逐月降雨量數據，按照克里金插值法生產降雨量等值線。同時考慮到地質災害易發性和降雨影響的重要性不同，本次按照專家打分法確定了調查區地質災害易發性和降雨影響指標的權重，分別給予權重值為0.8和0.2，利用GIS柵格計算功能對地質災害易發性和降雨量兩個評價因子進行歸一化處理，之後柵格疊加運算，得到危險性歸一化結果，數值越大反映地質災害危險性越大。利用ArcGIS軟體的自然間斷點分級方法，根據單元信息量累積頻度分布曲線特徵，最終將危險性評價結果分為四級：極高危險性、高危險性、中危險性和低危險性。
5.2地質災害易損性評價
易損性是指承災體可能遭受地質災害破壞的嚴重程度。易損性評價流程主要包括三部分：1）評價因子的確定；2）易損性評價權重的確定；3）在完成易損性評價因子量化的前提上，利用ArcGIS軟體疊加各評價指標得出評價單元易損性最終值，再將值進行重分類即等級劃分，完成沁源縣地質災害易損性劃分。將承災體易損性分為建築物易損性、人員易損性、交通設施易損性三大類，分別評價各自易損性，而後將不同類型承災體易損性按照一定的權重疊加，獲得綜合易損性評價圖。本次評價因子的權重採用專家打分法，確定人員、建築物、交通設施的權重值分別為0.6、0.3、0.1。將不同類型承災體易損性進行疊加，獲得綜合易損性評價圖，並採用自然斷點法將易損性分為四級：極高易損性、高易損性、中易損性、低易損性。
5.3地質災害風險評價
地質災害風險是指在一定區域和時期內，各類承災體因地質災害而造成的損失的可能性。地質災害風險評價即對發生地質災害的可能性及預計造成的人身安全和財產損失財進行的定量化的分析與評價。本次一般調查區風險評價是以10×10m柵格單元為表達形式的區域地質災害風險評價。考慮到承災體本身對災害的抵禦能力和相同結構的承災體在不同危險區的危險程度，充分結合地質災害危險性評價及易損性評價模型，共同建立地質災害定量風險評價體系，地質災害的危險性和易損性評價結果疊加運算採用矩陣分析方法得到了研究區各評價單元的風險性值，為保證計算在同一尺度空間，地質災害危險性和易損性評價指標均應當是歸一化處理後的數值。風險性採用矩陣計算公式如下：R=H×V式中，R——地質災害風險指數，H——區域地質災害危險性值，V——區域地質災害易損性值。根據上述得到的危險性評價和易損性評價結果得出區域地質災害風險性值，根據自然斷點法將風險性值分為四級，分別為極高風險、高風險、中風險、低風險。對風險性信息量圖進行概化和部分手動調整後，將整個沁源縣地質災害風險性分為極高風險區、高風險區、中風險區和低風險區4個大區，並結合承災體種類及其分布的區域，進一步分為9個亞區，其中極高風險區1個，高風險亞區2個，中風險亞區4個，低風險亞區2個。
5.4成果資料庫建設
GIS軟體在柵格矢量轉換過程中，會產生許多假的面集和許多面積很小或不協調面集單元，再次通過GIS的消除、融合功能，消除不合理因素，最終得到評價單元面數據集，並建立空間資料庫。
6結語
利用GIS強大的分析能力和空間信息管理能力，在地質災害的研究中給予了強大的研究方法和解決方案，它不僅能夠快速管理各種外動力因素，還能有效地應對影響地質災害的各種內在因素。但信息量模型還有待進一步研究並在實踐中完善，例如影響因素的確定、影響因子的確定、臨界值的確定。其中，專家的實地調查和經驗尤為重要，調查和經驗是任何技術和模型所無法替代的，如何將專家意見和模型相結合有待繼續研究。
參考文獻：
[1]吳信才.MAPGIS地理信息系統[M].北京:電子工程出版社,2016.
[2]吳樹仁,董誠,石菊松,等.地質災害信息系統研究——以重慶市豐都縣為例[J].第四紀研究,2015,23(6):683-692.
[3]張曉博,蔡恩琪,等.GIS技術在地質災害預測中的應用[J].科技與企業,2013(16):234-235.
作者:何濤 單位:山西省地質勘查局二一二地質隊
關鍵詞：GIS；信息量模型法；自然間斷點分級法
地質災害包括地裂縫、地面沉降、崩塌、滑坡、地震、火山、泥石流等，是指人為或者自然因素的作用下形成的對人類的生命和財產、環境造成損失和破壞的地質作用現象。我國是地質災害嚴重的國家之一，每年的地質災害都威脅著人民的生命和財產安全，也制約著國民經濟的發展，基於GIS中信息量模型法和自然間斷點分級法對地質災害易發性和危險性分區進行評價和分析不僅符合黃土高原地區地質災害點的實際情況，對模型參數也有很好的擬合。
1GIS的概念
GIS是一種新型的技術，它是在計算機硬、軟體系統支持下，對整個或部分地球表層（包括大氣層）空間中的有關地理分布數據進行採集、儲存、管理、運算、分析、顯示和描述的技術系統。它可以按照空間位置或地理坐標對空間數據進行處理、研究空間的實體、對數據進行有效的管理。GIS是一門信息產業，廣泛應用於各個行業。世界上GIS軟體就有四百多種，它們風格不同，大小也不一樣。在國外有ARCGIS、MGE、GENAMAP；國內有CITYSTAR、MAPGIS和Geostar等。
2GIS在國內外的應用現狀
雖然有很多GIS軟體，但是大體上分為在GIS的基礎上，應用函數庫二次開發出特有的地理信息系統軟體；應用GIS系統來處理用戶數據這兩種情況。現今已經應用到設施管理、資源管理、城市和區域的規劃、教育、石油和天然氣自動製圖等方面。目前GIS已經成功應用在政府管理、資源保護、環境保護、城市規劃建設等很多領域。隨著我國經濟的發展，GIS在測繪、交通運輸、軍事等領域發揮著舉足輕重的作用，GIS的應用主要包括綜合分析評價與模擬預測；在地理空間數據管理中的應用；建立專題信息系統和區域信息系統；在地圖製圖中的應用等。
3GIS信息量模型法在地質災害風險評價中的應用
3.1地質災害易發性評價方法確定
地質災害易發性定量評價方法較多，主要模型有信息量模型、專家系統模型、灰色系統模型、非線性模型及模式識別模型等。考慮到各模型的可行性、數據的易獲取性，本次沁源縣地質災害易發性評價採用信息量模型法。信息量模型反映了一定地質環境下最易致災因素及其細分區間的組合；具體是通過特定評價單元內某種因素作用下地質災害發生頻率與區域地質災害發生頻率相比較實現的。對應某種因素特定狀態下的地質災害信息量公式可表示為：式中：I——對應因素A、j狀態（或區間）下地質災Aj→B害B發生的信息量；N——對應因素A、j狀態（或區間）下地質災害分布的j單元數；N——調查區已知有地質災害分布的單元總數；S——因素A、j狀態（或區間）分布的單元數；S為調j查區單元總數。當I>0時，反映了對應因素A、j狀態（或區間）下地Aj→B質災害發生傾向的信息量較大，地質災害發生的可能性較大，或者說利於地質災害發生；當I<0時，表明因素Aj→BA、j狀態（或區間）條件下，不利於地質災害發生；當I=0時，表明因素A、j狀態（或區間）不提供有關地質災Aj→B害發生與否的任何信息，即因素A、j狀態（或區間）可以剔除掉，排除其作為地質災害預測因子。
3.2評價單元及評價指標選取
本次該縣一般調查區採用10m×10m大小的柵格單元進行，經統計共劃分為254.9萬個評價單元。本次地質災害易發性評價側重的是崩塌、滑坡等地質災害發育的數量多少及其活躍程度。評價指標體系主要為地質災害形成條件，包括地質環境條件（地貌、坡度、坡高、工程地質岩組、河流、植被）和地質災害誘發因素（人類工程活動）兩大類。
3.3各單因素評價指標信息量的確定
本次信息量模型中求出每一因子圖層中各類型地質災害的分布密度。為統計地質災害在每種因素的各個類別中的分布情況，需要利用空間分析功能來實現，在ArcGIS平台下計算各因子圖層各類別的信息量值。利用空間分析模塊中的區域統計功能，將地質災害點文件導入各單因素柵格文件後，得到在相應區域內地質災害點數目，用以上計算結果及已獲調查數據，得到地質災害信息量值，最終得到各評價指標的信息量表，見表1。
4信息量指標權重分配
根據各單因素信息量計算結果，將各種狀態的信息量值按大小排序，較前的狀態分別為：1）在距離道路30-60m區間內；2）坡高>30m；3）坡度>60°；4）距離居民地<10m區間內；5）距離河流在100-200m範圍內；6）植被指數0-0.4（覆蓋率低）；7）粘性土岩土體類型；8）梁峁狀黃土丘陵類型；9）坡型>0.5（權重分配見表2）。上述單因素分布區間屬本次一般調查區地質災害易發的區域，這些因素對崩塌等地質災害的形成發育起決定性的作用，是地質災害形成的主要控制條件。
5應用自然間斷點分級法評價易發分區
利用GIS柵格計算功能對各評價因子進行柵格疊加運算，得到各評價單元信息權值結果，數值越大反映各因素對地質災害發生的貢獻越大，發生地質災害的可能性就越大。經過ArcGIS軟體計算得到本次地質災害易發性評價信息量值區間為[-0.5937，1.7608]，並根據信息量值結果，利用ArcGIS軟體的自然間斷點分級方法，根據單元信息量累積頻度分布曲線特徵和曲線明顯分布的拐點，將本區的地質災害易發性劃分為極高、高、中、低四個等級。其中極高易發區信息量值區間為[1.7608，0.2154]，高易發區信息量值區間為[0.2154，0.0255]，中易發區信息量值區間為[0.0255，-0.2330]，低易發區信息量值區間為[-0.2330，-0.5937]。
5.1地質災害危險性評價
本次一般調查區地質災害危險性評價是在地質災害易發性評價基礎上，採用歷史月累計降雨量（大於5年）進行評價。本次調查收集沁源縣近40年（1980—2019年）歷史逐月降雨量數據，按照克里金插值法生產降雨量等值線。同時考慮到地質災害易發性和降雨影響的重要性不同，本次按照專家打分法確定了調查區地質災害易發性和降雨影響指標的權重，分別給予權重值為0.8和0.2，利用GIS柵格計算功能對地質災害易發性和降雨量兩個評價因子進行歸一化處理，之後柵格疊加運算，得到危險性歸一化結果，數值越大反映地質災害危險性越大。利用ArcGIS軟體的自然間斷點分級方法，根據單元信息量累積頻度分布曲線特徵，最終將危險性評價結果分為四級：極高危險性、高危險性、中危險性和低危險性。
5.2地質災害易損性評價
易損性是指承災體可能遭受地質災害破壞的嚴重程度。易損性評價流程主要包括三部分：1）評價因子的確定；2）易損性評價權重的確定；3）在完成易損性評價因子量化的前提上，利用ArcGIS軟體疊加各評價指標得出評價單元易損性最終值，再將值進行重分類即等級劃分，完成沁源縣地質災害易損性劃分。將承災體易損性分為建築物易損性、人員易損性、交通設施易損性三大類，分別評價各自易損性，而後將不同類型承災體易損性按照一定的權重疊加，獲得綜合易損性評價圖。本次評價因子的權重採用專家打分法，確定人員、建築物、交通設施的權重值分別為0.6、0.3、0.1。將不同類型承災體易損性進行疊加，獲得綜合易損性評價圖，並採用自然斷點法將易損性分為四級：極高易損性、高易損性、中易損性、低易損性。
5.3地質災害風險評價
地質災害風險是指在一定區域和時期內，各類承災體因地質災害而造成的損失的可能性。地質災害風險評價即對發生地質災害的可能性及預計造成的人身安全和財產損失財進行的定量化的分析與評價。本次一般調查區風險評價是以10×10m柵格單元為表達形式的區域地質災害風險評價。考慮到承災體本身對災害的抵禦能力和相同結構的承災體在不同危險區的危險程度，充分結合地質災害危險性評價及易損性評價模型，共同建立地質災害定量風險評價體系，地質災害的危險性和易損性評價結果疊加運算採用矩陣分析方法得到了研究區各評價單元的風險性值，為保證計算在同一尺度空間，地質災害危險性和易損性評價指標均應當是歸一化處理後的數值。風險性採用矩陣計算公式如下：R=H×V式中，R——地質災害風險指數，H——區域地質災害危險性值，V——區域地質災害易損性值。根據上述得到的危險性評價和易損性評價結果得出區域地質災害風險性值，根據自然斷點法將風險性值分為四級，分別為極高風險、高風險、中風險、低風險。對風險性信息量圖進行概化和部分手動調整後，將整個沁源縣地質災害風險性分為極高風險區、高風險區、中風險區和低風險區4個大區，並結合承災體種類及其分布的區域，進一步分為9個亞區，其中極高風險區1個，高風險亞區2個，中風險亞區4個，低風險亞區2個。
5.4成果資料庫建設
GIS軟體在柵格矢量轉換過程中，會產生許多假的面集和許多面積很小或不協調面集單元，再次通過GIS的消除、融合功能，消除不合理因素，最終得到評價單元面數據集，並建立空間資料庫。
6結語
利用GIS強大的分析能力和空間信息管理能力，在地質災害的研究中給予了強大的研究方法和解決方案，它不僅能夠快速管理各種外動力因素，還能有效地應對影響地質災害的各種內在因素。但信息量模型還有待進一步研究並在實踐中完善，例如影響因素的確定、影響因子的確定、臨界值的確定。其中，專家的實地調查和經驗尤為重要，調查和經驗是任何技術和模型所無法替代的，如何將專家意見和模型相結合有待繼續研究。
參考文獻：
[1]吳信才.MAPGIS地理信息系統[M].北京:電子工程出版社,2016.
[2]吳樹仁,董誠,石菊松,等.地質災害信息系統研究——以重慶市豐都縣為例[J].第四紀研究,2015,23(6):683-692.
[3]張曉博,蔡恩琪,等.GIS技術在地質災害預測中的應用[J].科技與企業,2013(16):234-235.
作者:何濤 單位:山西省地質勘查局二一二地質隊
收藏