工程地質(zhì)

基于GIS的滇藏鐵路麗江—香格里拉段工程地質(zhì) 條件分區(qū)研究

  0 引 言
 
  滇藏鐵路是我國正在規(guī)劃建設(shè)的大型鐵路工程之一, 麗江—香格里拉段是繼大理—麗江段后進行的第二段規(guī)劃建設(shè)段, 該段地形地貌地質(zhì)條件非常復(fù)雜, 雖然經(jīng)過多輪論證, 線路仍難最后確定。按照初期規(guī)劃(圖1), 滇藏鐵路麗江—香格里拉段共有3個走向方案可以比選:(1)麗江—長松坪—虎跳峽上峽口—香格里拉方案(西線方案);(2)麗江—大具—白水臺—小中甸—香格里拉方案(組合方案);(3)麗江—大具—白水臺—天生橋—香格里拉方案(東線方案)。初步分析認為, 西線方案工程地質(zhì)條件相對較好, 可以作為推薦方案, 但該方案的工程地質(zhì)條件仍很復(fù)雜且存在許多重大工程地質(zhì)問題, 工程建設(shè)難度大。
 
  西南山區(qū)是我國西部開發(fā)工程建設(shè)的重點地區(qū)之一, 也是工程地質(zhì)條件最復(fù)雜的地區(qū)之一,隨著工程建設(shè)的深入, 工程地質(zhì)條件的復(fù)雜程度已經(jīng)成為制約工程建設(shè)的重要因素。工程地質(zhì)條件一直是工程地質(zhì)界研究的熱點, 丁繼新等對工程地質(zhì)條件的概念和應(yīng)用進行了相關(guān)的研究,但基于正在規(guī)劃建設(shè)的大型工程的工程地質(zhì)條件分區(qū)評價方面的研究則較少。本文以滇藏鐵路麗江—香格里拉段為例, 闡述工程地質(zhì)條件分區(qū)評價在線路比選和優(yōu)化方面的應(yīng)用;在區(qū)域地殼穩(wěn)定性評價的基礎(chǔ)上, 將基于GIS技術(shù)的層次分析法引入到麗江—香格里拉段鐵路規(guī)劃區(qū)的工程地質(zhì)條件評價中;充分利用GIS技術(shù)處理海量數(shù)據(jù)信息的優(yōu)勢, 采用層次分析法模型, 進行麗江—香格里拉段鐵路規(guī)劃區(qū)的工程地質(zhì)分區(qū)評價, 評價結(jié)果對鐵路選線具有一定指導(dǎo)意義。
 
  1 區(qū)域工程地質(zhì)概況
 
  研究區(qū)在地貌上屬構(gòu)造剝蝕高中山與深切河谷區(qū), 地勢總體西北高、東南低。地面高程多在2 500 ~ 5 000 m, 最高峰為麗江北西面的玉龍雪山, 主峰扇子陡高程5 596 m, 最低處為麗江以北白馬廠一帶金沙江河谷, 高程約1 570m, 河流切割深度最大可達2 000 m以上。
 
  研究區(qū)位于中國兩大地貌階梯青藏高原東南緣向云貴高原的過渡部位, 為揚子地臺、松潘—甘孜塊體和三江塊體的結(jié)合部位, 構(gòu)造背景甚為復(fù)雜。在印度板塊向北推擠和青藏高原南南東向擠出疊加作用下, 新構(gòu)造運動十分強烈, 表現(xiàn)為強烈的垂直差異運動和塊體側(cè)向滑移, 以及NW向斷裂右旋位移、近SN向和NNE向斷裂左旋位移的斷裂活動特征, 主要斷裂有鶴慶—洱源斷裂帶、麗江—小金河斷裂帶、金沙江斷裂帶、塔城—紅巖斷裂、中甸斷裂帶、玉龍—哈巴雪山東麓斷裂、中甸—海羅斷裂、龍蟠—喬后斷裂帶、小中甸—大具斷裂帶等。研究區(qū)新構(gòu)造運動強烈, 地震頻發(fā), 上新世以來發(fā)生多次構(gòu)造活動。研究區(qū)自前古生代至新生代的地層均有出露。
 
  2 評價方法簡介
 
  工程地質(zhì)條件分區(qū)評價的核心是確定各種影響因素對工程地質(zhì)條件影響的大小、分布與發(fā)育強度, 計算出不同地區(qū)的工程地質(zhì)條件指數(shù)作為量化指標。結(jié)合研究區(qū)的工程地質(zhì)條件和實際工程規(guī)劃現(xiàn)狀, 采用基于GIS的層次分析法作為工程地質(zhì)分區(qū)評價的研究方法。層次分析法適用于多準則、多目標復(fù)雜問題的決策分析, 可以將決策者對復(fù)雜系統(tǒng)的決策思維過程實行數(shù)量化, 為選出最優(yōu)決策提供依據(jù)[ 6] 。經(jīng)過多年的研究和應(yīng)用實踐, 不少研究者開始將GIS技術(shù)與AHP方法相結(jié)合, 大大提高了傳統(tǒng)的AHP方法在地學(xué)研究中的應(yīng)用效果。
 
  本文采用基于GIS的工程地質(zhì)條件分區(qū)評價方法, 在確定研究區(qū)、研究對象和研究目標后,對各種資料進行數(shù)據(jù)處理并建立相應(yīng)的空間數(shù)據(jù)庫;根據(jù)研究目標的特征, 分析影響目標的因素,建立目標的層次指標模型和層次結(jié)構(gòu), 構(gòu)造判斷矩陣;在專家對影響因素進行綜合評分的基礎(chǔ)上,進行層次單排序、求解權(quán)向量和一致性檢驗, 從而獲得各指標因素值, 并運用GIS空間分析功能提取分析因子。數(shù)據(jù)處理計算采用大型地理信息系統(tǒng)ArcGIS9.2 軟件, 對研究區(qū)域進行柵格化,每一個柵格作為模型評價的一個運算單元, 并將數(shù)據(jù)庫中的數(shù)據(jù)按照規(guī)則進行柵格化處理。再采用圖形疊加的模型評價方式, 將參與評價的各個因素權(quán)值分配到不同的柵格上。將各個因素進行圖形疊加, 對屬性值進行代數(shù)運算, 再將疊加后的柵格數(shù)據(jù)化生成新的圖形, 并形成最終評價結(jié)果。評價流程見圖2。
 
  
  本文將根據(jù)計算獲得的工程地質(zhì)條件指數(shù)值的分布范圍, 結(jié)合野外實際調(diào)查情況進行驗證和分析, 對各工程地質(zhì)區(qū)域的工程地質(zhì)條件進行分區(qū)評價, 對鐵路工程建設(shè)的影響及適宜性進行討論和說明。
 
  3 評價指標的確定和評價模型
 
  3.1 評價指標的確定
 
  研究區(qū)位于我國的西南山區(qū), 工程地質(zhì)條件復(fù)雜, 影響工程地質(zhì)條件的因素較多, 在充分考慮各種因素的基礎(chǔ)上, 選取地形地貌、工程地質(zhì)巖組、斜坡結(jié)構(gòu)、地質(zhì)災(zāi)害發(fā)育現(xiàn)狀、潛在震源區(qū)、活動斷裂、微地貌類型(地形與鐵路設(shè)計高程間的高差)、人類工程活動、降水量(主要考慮垂直降水量的差別)、與溝谷間的距離等10個因素作為一級評價指標。
 
  3.2 層次結(jié)構(gòu)模型
 
  根據(jù)評價指標體系組建和操作的基本原則,采用分層遞階方法, 將研究區(qū)工程地質(zhì)條件評價指標體系分為總體目標層、約束層(一級指標層)、評價指標層(二級指標層)和對象層(評價對象層)4個層次(圖3)。
 
  目標層是系統(tǒng)分析的最高層, 用以表達鐵路工程建設(shè)的適宜性。在本次分析中, 確定工程地質(zhì)條件分區(qū)是總目標。按照影響鐵路工程建設(shè)的邏輯關(guān)系, 分列地形坡度、工程地質(zhì)巖組、斜坡結(jié)構(gòu)、地質(zhì)災(zāi)害發(fā)育程度、潛在震源區(qū)、活動斷裂發(fā)育狀況、微地貌類型、人類工程活動、降水量等方面的因素, 分別考慮其對總目標的影響,這些影響因素構(gòu)成約束層。在約束層的基礎(chǔ)上,對其中各種因素具體細分, 一共有32個指標層,這些指標層的權(quán)值在上述約束層權(quán)重計算的基礎(chǔ)上進行劃分。對象層是系統(tǒng)分析的最底層, 在本次研究中, 把工程地質(zhì)條件分區(qū)作為對象層。
 
  4 評價過程及評價結(jié)果
 
  4.1 評價指標的量化途徑
 
  在上述評價指標確定后, 充分利用GIS技術(shù)強大的基礎(chǔ)數(shù)據(jù)處理和空間分析功能, 在ArcGIS9.2平臺上形成地形坡度、工程地質(zhì)巖組、斜坡結(jié)構(gòu)類型區(qū)、地質(zhì)災(zāi)害發(fā)育程度、潛在震源區(qū)、活動斷裂, 微地貌類型、人類工程活動、降水量柵格、水系距離分析等相關(guān)的柵格文件格式的專題圖(圖4(a)— (d))。基于上述專題圖層, 對于能夠直接量化的指標, 可以在矢量化的專題圖層提取相應(yīng)的數(shù)據(jù)信息, 然后對指標進行等級劃分并賦值;對于不能直接量化的指標, 采用評分比較的方法進行分區(qū)劃分等級并賦值。以上可獲得各評價指標的單因素等級量化結(jié)果。根據(jù)研究區(qū)范圍和工程地質(zhì)條件特征, 將柵格大小定為50 m×50 m, 將6 621.5 km2的研究區(qū)劃分為2 648 600個柵格單元。
 
  4.2 影響因素和指標權(quán)重的確定
 
  影響因素和評價指標的權(quán)重計算采用在專家打分法確定各影響因素影響因子大小的基礎(chǔ)上,構(gòu)造相關(guān)判斷矩陣, 求解得到影響因素和指標的權(quán)重, 表1為在匯總分析專家打分的基礎(chǔ)上建立的因素判斷比較矩陣。經(jīng)過計算, λmax =10.369,CI=0.041 0, RI=1.49, CR=0.027 5 <0.1, 符合一致性判據(jù), 各影響因素權(quán)重分配見表2。
 
  從權(quán)重分配表中可以看出, 地質(zhì)災(zāi)害發(fā)育程度和活動斷裂的權(quán)重在各因素中位于前列, 其次是工程地質(zhì)巖組, 這3 個因素是進行工程地質(zhì)條件評價的基礎(chǔ), 是影響鐵路工程選線和建設(shè)的先決條件;地形坡度和斜坡結(jié)構(gòu)是影響區(qū)域崩塌、滑坡泥石流的重要條件, 所以權(quán)重也較高;微地貌類型、人類工程活動和降水量的權(quán)值差異不大。
 
  以上分析表明, 通過專家打分法和層次分析法相結(jié)合確定的影響因素權(quán)重是符合客觀實際的。
 
  在各影響因素權(quán)值確定的基礎(chǔ)上, 對各評價指標層進行權(quán)重劃分和賦值(表3), 其中, 活動斷裂指標又按照與活動斷裂間的遠近程度進行了距離分析, 不同距離范圍的權(quán)值見表4。數(shù)值的大小表示其對鐵路工程地質(zhì)條件的影響程度, 數(shù)值越大, 表示該因素越不利于鐵路工程建設(shè)。
 
  4.3 評價計算和結(jié)果分析
 
  4.3.1 計算過程
 
  利用ArcGIS軟件空間分析的空間信息再分類功能, 按照評價指標權(quán)重對各個柵格文件的不同屬性類別進行賦值。采用ArcGIS軟件的柵格計算功能, 對經(jīng)過空間信息再分類處理的10個影響因素的柵格文件進行各柵格權(quán)重值的求和運算, 根據(jù)式(1)計算每個柵格的工程地質(zhì)條件指數(shù)值。計算結(jié)果表明研究區(qū)內(nèi)各柵格的工程地質(zhì)條件指數(shù)值分布在0.054 ~ 0.303之間, 集中分布在0.080~ 0.230之間, 平均值為0.150, 標準差為0.034,變異系數(shù)為0.227。工程地質(zhì)條件指數(shù)值分布的平穩(wěn)性較好, 反映了工程地質(zhì)條件在區(qū)域上的變化具有一定的遞變性, 與實際工程地質(zhì)情況相符。
 
  綜合考慮計算結(jié)果、野外地質(zhì)調(diào)查工程地質(zhì)勘察資料, 確定了用于進行工程地質(zhì)條件分區(qū)的指數(shù)閥值, 將研究區(qū)工程地質(zhì)條件分為較好(B≤0.130)、中等(0.130 <B≤0.155)、較差(0.1550<B≤0.200)和差(0.200 <B≤0.310)4級, 并在ArcGIS軟件平臺自動生成工程地質(zhì)條件區(qū)劃柵格圖, 經(jīng)過局部平滑和噪音處理, 獲得研究區(qū)鐵路規(guī)劃區(qū)的工程地質(zhì)條件計算成果圖。
 
  4.3.2 結(jié)果分析
 
  根據(jù)研究區(qū)工程地質(zhì)條件評價結(jié)果圖(圖4(e))和相關(guān)工程地質(zhì)資料, 將各分區(qū)中的主要工程地質(zhì)問題分析如下。
 
  (1)工程地質(zhì)條件差。該類分區(qū)主要分布于玉龍雪山東麓文筆水庫—麗江市西側(cè)—玉湖—玉龍雪山—大具—哈巴雪山(古魯巴)—俄迪條帶狀區(qū)域內(nèi), 在香格里拉縣城附近及仁河—龍蟠一帶也有小面積分布。影響本類工程地質(zhì)分區(qū)的最主要因素是活動斷裂, 突出表現(xiàn)為與強活動斷裂距離較近, 特別是在強活動斷裂1 km范圍內(nèi), 如玉龍雪山東麓和哈巴雪山中部地區(qū);在多條斷裂相交或近于相交的地區(qū)工程地質(zhì)條件亦較差, 如麗江市縣城至文筆水庫一帶處于多條中等-強活動斷裂的交匯部位,在該區(qū)域內(nèi)工程地質(zhì)條件差;地質(zhì)災(zāi)害發(fā)育地帶、軟弱松散第四系土石類、順向坡及斜坡坡度大于45°的地帶一般也位于工程地質(zhì)條件差的區(qū)域。該類分區(qū)約占研究區(qū)總面積的6.52%。
 
  (2)工程地質(zhì)條件較差。該類分區(qū)主要分布在工程地質(zhì)條件差區(qū)域的外圍, 其一方面受活動斷裂的影響, 多位于距活動斷裂1 ~ 5 km的范圍內(nèi);另一方面, 雖然部分區(qū)域與強活動斷裂之間的距離小于1 km, 如玉湖和古魯巴附近, 但由于地形較平坦或斜坡平緩, 工程地質(zhì)巖組和斜坡結(jié)構(gòu)好,遠離溝谷影響, 因而工程地質(zhì)條件屬較差級。此外, 在仁和以西的虎跳峽鎮(zhèn)至哈巴雪山之間的順向坡發(fā)育地帶工程地質(zhì)條件較差, 這些地區(qū)斜坡結(jié)構(gòu)因素起主導(dǎo)作用。該類分區(qū)約占研究區(qū)總面積的32.24%。
 
  (3)工程地質(zhì)條件中等。該類分區(qū)多呈點狀和塊狀分布, 在麗江市團山地區(qū)、拉市海一帶、月亮坪—俄迪一帶、小中甸鎮(zhèn)—碧古—阿熱一帶以及那帕海地區(qū)分布較多。該類工程地質(zhì)條件分區(qū)多位于中等活動斷裂和弱活動斷裂的影響范圍內(nèi),距離溝谷水系近, 地質(zhì)災(zāi)害較發(fā)育, 地形坡度小、斜坡結(jié)構(gòu)一般。該類分區(qū)約占評估區(qū)總面積的32.41%。
 
  (4)工程地質(zhì)條件較好。該類分區(qū)主要分布在4個地帶, 即虎跳峽鎮(zhèn)一帶、團山—文化村一帶、安南地區(qū)和阿熱北部區(qū)域, 在研究區(qū)其他區(qū)域僅零星分布。該類地區(qū)遠離活動斷裂帶, 地質(zhì)災(zāi)害不發(fā)育, 工程地質(zhì)巖組以較堅硬中厚層狀砂板巖、玄武巖巖組和堅硬塊狀碳酸鹽巖巖組為主, 地形坡度一般小于15°, 斜坡結(jié)構(gòu)好, 以橫向坡、反向坡和平坡為主。該類分區(qū)約占研究區(qū)總面積的28.83%。
 
  5  線路方案工程地質(zhì)條件評價和優(yōu)化
 
  從工程地質(zhì)分區(qū)圖(圖4(e))可以看出, 研究區(qū)內(nèi)工程地質(zhì)條件差的區(qū)域集中分布在麗江縣城至古魯巴條帶狀區(qū)域內(nèi), 麗江—香格里拉段鐵路東線方案、西線方案和組合方案都不可避免地通過這些區(qū)域。其中, 東線方案工程地質(zhì)條件較好段約占10%, 工程地質(zhì)條件中等段約占40%, 工
 
  程地質(zhì)條件較差段約占30%, 工程地質(zhì)條件差段約占20%;主要工程地質(zhì)問題是大具盆地段位于
 
  哈巴—玉龍雪山東麓斷裂5 km影響范圍內(nèi), 斜坡結(jié)構(gòu)以順向坡為主, 滑石板滑坡和金沙江深切河谷的地貌特征增加了工程地質(zhì)條件的復(fù)雜程度,阿熱到香格里拉段位于中甸斷裂的影響范圍內(nèi)。
 
  西線方案工程地質(zhì)條件較好段約占10%, 工程地質(zhì)條件中等段約占20%, 工程地質(zhì)條件較差段約占30%, 工程地質(zhì)條件差段約占40%;主要工程
 
  地質(zhì)問題是線路在麗江盆地內(nèi)與玉龍雪山東麓斷裂斜交, 尤以小玉龍附近影響最大, 從小玉龍附近開始, 鐵路設(shè)計高程和玉龍雪山山體之間的高差大, 鐵路呈NW向以隧道(玉峰寺隧道)形式通
 
  過玉龍雪山, 隧道最大埋深達1 000 m以上, 受這兩種因素的控制, 該段工程地質(zhì)條件差。組合方案工程地質(zhì)條件較好段約占15%, 工程地質(zhì)條件中等段約占25%, 工程地質(zhì)條件較差段約占30%, 工程地質(zhì)條件差段約占30%, 該方案與哈巴—玉龍雪山東麓斷裂近直交, 延伸至俄迪一帶,影響本方案的突出地質(zhì)問題是活動斷裂。
 
  在上述工程地質(zhì)條件分區(qū)和線路規(guī)劃方案分析的基礎(chǔ)上, 認為西線方案可以作為麗江—香格里拉段鐵路推薦線路方案, 但該方案南段存在較大工程地質(zhì)問題———玉龍雪山活動斷裂影響和玉峰寺深埋隧道, 需要進一步優(yōu)化。根據(jù)工程地質(zhì)分區(qū)評價結(jié)果, 建議西線方案由原來的從團山附近接入大麗線(大理—麗江鐵路)改為從麗江市南端太平村附近接入, 線路總體走向為太平村—中吉村—文筆水庫—下長水—三家村—仁河—中義—龍蟠, 在龍蟠處接原西線方案, 這樣, 可以減弱玉龍雪山東麓斷裂對鐵路的影響以及避免玉峰寺深埋隧道引起的相應(yīng)的重大工程地質(zhì)問題。
 
  6 結(jié) 語
 
  基于ArcGIS平臺和層次分析法對滇藏鐵路麗江—香格里拉段的工程地質(zhì)條件進行了分區(qū)與評價研究, 得到如下認識:
 
  (1)滇藏鐵路麗江—香格里拉段地形地貌和地質(zhì)條件非常復(fù)雜, 新構(gòu)造運動強烈, 影響鐵路工程地質(zhì)條件的因素較多, 主要有地形坡度、工程巖組、斜坡結(jié)構(gòu)類型、地質(zhì)災(zāi)害發(fā)育程度、潛在震源區(qū)、活動斷裂、微地貌類型、人類工程活動、降水量、水系距離等10個因素。
 
  (2)經(jīng)過專家打分法和層次分析法計算認為地質(zhì)災(zāi)害發(fā)育程度和活動斷裂的權(quán)重在各評價因素中位于前列, 其次是工程地質(zhì)巖組, 這3個因素是進行工程地質(zhì)條件評價的基礎(chǔ), 是影響鐵路工程選線和建設(shè)的先決條件。
 
  (3)根據(jù)計算所得的工程地質(zhì)條件指數(shù)值分布范圍和野外實際調(diào)查結(jié)果, 將研究區(qū)工程地質(zhì)條件分為差、較差、中等和較好等4 類。工程地質(zhì)條件差的地區(qū)約占總面積的6.52%, 該類地區(qū)集中在玉龍雪山東麓和哈巴雪山地區(qū);工程地質(zhì)條件較差的地區(qū)約占研究區(qū)總面積的32.24%, 在研究區(qū)內(nèi)分布面積較大;工程地質(zhì)條件中等的地區(qū)約占研究區(qū)總面積的32.41%;工程地質(zhì)條件較好的地區(qū)約占研究區(qū)總面積的28.83%。
 
  (4)綜合分析表明西線方案的工程地質(zhì)條件總體優(yōu)于其他比選方案, 基于評價結(jié)果對西線線路方案進行了優(yōu)化, 優(yōu)化后的方案減弱了玉龍雪山東麓斷裂對鐵路的影響, 可避免玉峰寺深埋隧道引起的工程地質(zhì)問題。