摘要：風沙地區鐵路選線主要受制于地質因素、環境因素，實現沿線區域地質環境建模是減災選線及線路優化設計的基礎。為了提高沙害區線路方案設計的效率與水平，借助遙感技術實現鐵路沿線風沙災害信息識別與提取，設計多源異構地質數據的入庫規則與存儲方式，建立基于 GIS 的風沙地區鐵路選線不良地質信息庫，實現災害數據的統一管理;利用獲取的選線區域矢量、柵格、專題數據，基于 GIS 建模功能構建風沙區選線地質環境模型，實現不良地質信息的多維表達與空間分析，通過可視化的立體選線方法進行線路方案的優化設計。以格庫鐵路某沙害區選線為例，建立線路三維地質環境實體模型，通過分析風沙流對平面穿越和立體繞避 2 個線路方案的危害程度，得到了較為合理的設計結果。研究表明本方法能夠為今后風沙地區鐵路選線設計提供一定的技術支撐。

　　關鍵詞：風沙地區;鐵路選線;不良地質建模;GIS;優化設計

　　引言

　　地質問題一直是影響鐵路選線設計的重要因素，不良地質的空間分布及危害程度決定了線路走向、構造物布設和工程投資等內容。風沙危害是一種由風沙流運動而引起的典型的不良地質災害，鐵路沙害的主要形式為風蝕路基、鋼軌磨蝕和沙埋鐵路，其貫穿沙區鐵路建設、運營全過程[1-3]。

　　風沙地區鐵路選線設計需要著重考慮沙害問題，采取合理的線路方案或沙害整治措施，盡可能減弱、降低風沙流對線路工程的影響。根據《國家綜合立體交通網規劃綱要》(2021)、《國家“十四五”發展規劃》(2021)和《西部陸海新通道總體規劃》(2019)，西北地區鐵路在未來一段時期建設量仍然很大，而這些區域戈壁、沙漠廣泛分布，擬修建的很多條線路都穿行于風沙危害頻發區域，風沙地區鐵路選線的技術水平、作業效率面臨更大的挑戰，如何利用先進的計算機信息建模技術為選線服務是風沙地區鐵路選線急需解決的問題。

　　風沙災害是一種發生于地表的地質危害，通過利用現代化勘測設計技術對選線范圍內的沙害區進行識別、提取與表達，建立線路沿線的帶狀地理環境模型，進而分析線路與不良地質的空間關系來實現方案優化設計，是當前風沙地區鐵路減災選線設計的必然要求與發展趨勢。易思蓉等[4]最早提出了虛擬環境選線系統的智能環境模型，通過研發虛擬選線系統模擬真實地理環境來輔助鐵路選線設計。蒲浩等[5]采用 OSG 技術二次開發了人機交互式設計的鐵路三維可視化設計平臺，實現了鐵路三維場景的快速構建。

　　呂希奎等[6]利用遙感技術提取線路沿線區域地質信息，對基于地質建模的復雜艱險山區鐵路三維選線技術進行了研究。陳情等[7]利用Google Earth 影像數據和遙感技術動態解譯了滇藏鐵路選線范圍內的不良地質，實現了高海拔大高差地區地質選線的方案比選。RS、GIS 技術在鐵路選線領域中多被用來提取沿線的地質信息[8-9]，在風沙鐵路中的應用主要是沙害區域的識別，且應用到風沙區鐵路選線中主要以二維影像和文本為主，信息的利用率和直觀程度不夠[10]，適用于選線分析的信息多維表達研究也不深入，與地質環境可視化的空間立體選線目標還有一定的差距，已不能適應于新時期風沙地區鐵路減災選線的復雜性。

　　鑒于此，本文提出矢量數據、柵格數據、專題數據相集成的風沙地區鐵路選線地質環境建模方法，首先通過遙感技術實現地物及不良地質的信息提取，利用多源數據和 GIS 技術建立風沙地區鐵路選線地質環境數據庫，基于 GIS 平臺將數字地形數據和影像數據集成構建三維地形環境，然后疊加不良地質形成可用于風沙地區鐵路選線的三維可視化地質環境，通過在三維地質環境中分析不良地質與線路的空間關系及風沙危害程度，實現線路方案的優化設計，可為提高風沙地區鐵路選線設計水平和效率提供技術支持。

　　1 基于 GIS 的風沙地區鐵路選線地質環境數據庫構建

　　1.1 地形數據獲取

　　三維地形環境是數字化鐵路選線設計的環境基礎，通過構建地形環境可對線路沿線地表信息特征進行三維可視化表達[11]。高程數據 DEM(DigitalElevation Model)和影像數據是建立三維虛擬地質環境的基礎，高程數據的獲取可利用互聯網提供的免費數據服務，從國家地理空間數據云平臺網站或美國地質調查局 UCGS網站進行數據下載，操作時通過上傳矢量邊界或選定行政區劃來制定下載的范圍。影像數據獲取常用的方法是通過無人機、航測飛機等航測方式來獲取影像或利用互聯網獲取免費影像資源。

　　1.2 不良地質數據

　　獲取風沙危害作為一種發生于地表的地質災害現象，具有明顯的平面特征。解譯遙感影像是識別地質信息的一種先進可靠技術，通過建立解譯標志和參照物，對擬選線區域內的遙感影像進行解譯可得到地物、水系、滑坡、泥石流、風沙區等信息[12]，還可以利用現有的地質調查數據對其進行補充。風速、風向決定著風沙運動的強度、路徑，通過國家氣象科學數據中心和架設移動式氣象站獲取區域氣象數據。

　　風沙區選線的不良地質信息組成，基于 GIS 技術可進一步分析不良地質的分布規律、發展趨勢和危害程度等。將線路沿線帶狀區域作為柵格數據操作的研究區，根據土地利用現狀圖、Landsat-5TM 遙感影像、中國沙漠分布圖等數據，借助軟件解譯功能實現監督分類，判別選線區域的不良地質、地物、植被等，沿線地物識別結果。著重提取風沙區時，可對影像進行 7、4、2 波段的假彩色合成，這種波段組合可有效區分沙丘、沙地與其他地物類型。

　　(2)矢量地質信息獲取不良地質分布及范圍識別以后，使用 ENVI 軟件的 AOI 模塊提取地質對象的空間位置、邊界點坐標、周長和范圍等信息，還可用按腌膜提取工具實現不良地質柵格格式數據的提取。

　　為了更直觀有效地進行三維可視化分析，將柵格數據、矢量數據疊加到一起實現不良地質的矢量柵格一體化表達。由于不同類型沙丘對線路設計的影響作用不同，在風沙區空間位置及范圍確定后，需要在遙感解譯的基礎上借助 Google Earth 軟件和現場監測、實地調研、既有研究成果數據對沙丘類型、移動特征、移動方向等信息進行補充完善。根據風沙地貌學相關知識[13-15]，并結合實地考察調研，分別建立流動沙丘/沙地、半固定沙丘/沙地、固定沙丘/沙地、戈壁風沙流、風蝕的解譯標志，其中建立的流動沙丘解譯標志。

　　利用 Google Earth 的添加路徑工具，選取風沙區界線的起點，依次選取下一點直至沙害區地質界線繪制完成，設置不同的文件名后保存為相應的KML 文件，將 GE 的 WGS84 地理坐標系轉換為需要的坐標系統后，基于 GIS 的點轉柵格功能、Excel轉表功能、TXT 文件構建矢量風沙區范圍，轉化為相應的面圖層，在面 shp 的 Table 中將沙害類型、移動特征、危害程度等信息添加到進去，利用分類統計結果可研究各類型沙害的空間分布格局及組合關系。

　　格庫鐵路沿線主要的風沙危害類型空間分布。由矢量數據可知，沿線范圍內風積沙地段 369.5km，戈壁風沙流地段 235.9km，風蝕地段 146.3km。根據風沙流對鐵路的危害類型及作用機理，本文將風沙類型主要分為流動沙丘、流動沙地、半固定沙丘、半固定沙地、戈壁風沙流、風蝕、固定沙丘 7 種類型，其中，流動沙地主要分布在羅布莊至通古斯巴、瓊吐爾至羅布莊、東柴山至大烏斯、茫崖湖東至茫崖湖等地段，是線路與防沙工程設計的重點區域。

　　(3)地形數據地面是風沙流運動的環境基礎，地表的粗糙度、起伏度、坡度、坡向等影響著沙源的分布，對風沙流的形成以及運動發展有重要的決定作用。這些數據雖然也是柵格形式，卻是不良地質區域平面表達分析的重要內容，需通過二次分析后得到，利用 GIS 的 Spatial Analyst 空間分析模塊處理 DEM 可獲得這類數據。

　　(4)風況數據風沙流是一種典型的氣固兩相流，只有當風速大于起沙風速時才會發生風沙流運動[16]，因此風況決定著區域內風沙流的危害程度，是風沙災害分析的重要內容。這類數據獲取的途徑主要是利用國家氣象站的監測數據和實地架設觀測儀器實現數據獲取，經規范化處理后可得到主導風向、平均風速、起沙風頻率、風玫瑰、輸沙勢等。基于 GIS 的非空間數據空間化方法，利用沿線氣象站點數據對格庫鐵路沿線風速數據空間化結果。

　　1.3 地質環境數據存儲設計

　　鐵路選線過程中使用的數據具有多源異構特性，基于 GIS 技術將選線區域不良地質體坐標位置作為各類型地質信息的組織框架，通過 GIS 建立的地質信息庫可實現地形、地質、地物、正射影像、沙害類型、地下水深等空間數據、非空間數據的存儲，將各種要素數據的特征信息以屬性表的形式存儲，并與矢量數據鏈接，進而利用分析模塊實現風沙地區選線不良地質信息的提取與表達。

　　Arc GIS主要是用 Geo database 實現數據存儲[17]，根據風沙地區鐵路選線數據要素需求，對地貌類型、沙丘疏密度、植被類型及覆蓋度和土壤類型等按要素特征分別存儲在不同的 Feature dataset 中，沙丘、沙地、林地及草地、自然保護區等以 polygon 格式存儲，褶皺、斷層、線路設計方案、風沙區邊界等以 line格式存儲，城鎮、工礦企業、居民點等以 point 格式存儲。

　　2 風沙區選線地質環境建模

　　2.1 建立三維地形環境數字高程模型 DEM 是在一個區域內以密集地形模型點的坐標表達地面形態和描述地表起伏情況，也可利用其提取各種地形參數，被廣泛應用于公路、鐵路選線設計領域[18]。

　　DEM 最主要的三種表示模型是不規則三角網模型、規則格網模型和等高線模型，其中 TIN 模型在線路設計方面應用最多。根據獲取的影像數據和高程數據，利用 GIS的表面模型將柵格數據轉化為 TIN 網，然后疊加該區域的影像數據，即可快速建立選線區域的三維地形環境。不同顏色表示不同的高程范圍，能比較直觀地觀察、分析風沙區的位置、范圍、面積等信息，結合區域內不良地質的分布可實現線路方案的初步選定。

　　2.2 平面表達分析

　　風沙地區鐵路選線設計更加強調了地質選線、環境選線、減災選線的理念，對各種不良地質的準確表達有助于避免線路設計修建及運營過程中的一些重大地質問題，特別是局部不良地質的有效表達更為重要。沙丘是在地表沙粒在風力的作用下搬運、堆積形成的地貌形態，根據沙丘流動程度可將其劃分為固定沙丘、半固定沙丘和流動沙丘[19]，其中流動沙丘會隨著時間發展向一定方向移動，且流動沙丘區域內的沙粒極易被風吹起，故這類沙丘會淹沒道路、村莊、農田等，對基礎設施、社會經濟的影響最大，而半固定、固定沙丘由于表面沙流相對穩定，發生風沙流危害現象的可能性較小，造成的破壞也較輕[20]。不同類型沙丘都具有特有的表面形態。

　　風沙區是一種具有典型平面特性的不良地質，利用遙感解譯得到的影像數據將不良地質矢量信息在地理環境中進行疊加，可以得到不良地質的表面區域。利用 ENVI 軟件獲取的風沙區邊界坐標文件和 GIS 地理數據庫新建 polygon 要素，實現不同類型沙丘的空間范圍勾繪。

　　沙丘地貌提取完成以后，按照沙丘的流動特征、形態特征對其分類，并將數字地貌類型轉化為相應的矢量圖。起伏度是指比較范圍內最大高程與最小高程之間的差值，經反復試驗，將起伏度的閾值設定為 16m 可較好地實現沙丘、沙地的類型區分。用 Arc Map 的點線面編輯工具勾畫沙丘單元，并添加相應屬性信息。由于部分區域存在非平滑的邊界折線，還需對其進行平滑處理。

　　利用 GIS 空間分析可實現面積統計、空間分布特征分析，選線設計中當線路穿過風沙區時，對有交集的不良地質進行高亮顯示，并根據 ID 調入已存儲的不良地質屬性信息，結合三維空間分析功能實現不良地質對線路影響程度的判別，另外將 2 期影像聯合分析還可得到沙丘的移動規律、空間范圍變化等。

　　2.3 立面表達分析

　　不良地質信息的立體表達主要是通過線路沿線的三維地形環境和鉆孔數據，采用面向對象的表達方法，利用 GIS 平臺將地質信息集成到三維環境中進行有效表達。將帶有位置信息的鉆孔數據轉化為三維空間點，分層次提取每一層鉆孔點，將獨立保存的鉆孔數據依次分別生成每一層的地質體，完成地質體的三維建模，立體表達結果見圖 13(a)。線路作為一個帶狀的三維空間實體，在不良地質區域往往會穿越斷層，而斷層的存在破壞了地層連續性，其對橋梁、隧道的局部方案設計有很大影響。而斷層區域的三維地質體建模就是通過分別存儲斷層兩側的鉆孔數據，在建模后將地層分界點進行連接即可。

　　2.4 不良地質體與地形融合

　　鐵路選線設計時，為了更好地分析對比線路方案，避免部分優秀方案的缺失，通過將生成的不良地質實體疊加到線路三維地形環境中，實現兩者的融合[21]。首先要確定融合區域，不良地質體與地形融合影響區域確定的實質就是將不良地質體的邊界投影到地形網格上，形成一個多邊形區域。兩者間的對應關系可以為一對多或者多對一，即多個融合區域對應一個地形網格。地質實體模型讀入后，可通過分解模型的包圍盒進行融合區域的確定。進行地質體的坐標轉換后，通過求解融合區域與地形塊網格，利用約束三角網刪除內部的三角形，以融合區域的邊界為特征邊對剩余的三角形進行重構網，不良地質體循環構網后實現地質體與地形的融合。

　　3 風沙區線路方案優化設計

　　3.1 工程概況格庫鐵路是連接青海省格爾木市和新疆自治區庫爾勒市的一條國鐵Ⅰ級客貨共線電氣化鐵路[22]，全長 1213.7km，設計速度 120km/h。某沙害段線路走行于山前沖洪積平原及低山丘陵區，地勢略起伏，為線路穿越東柴山小越嶺地段，附近發育有風蝕殘丘、風蝕洼地、流動沙地、半固定沙地沙丘及戈壁風沙流等不良地質，局部段落風積沙嚴重。風沙地區鐵路選線設計時，通過建立沿線三維地形環境和不良地質體模型，能夠較為全面地考慮風沙危害對線路平面走向的影響。根據沙害類型及程度，不僅可以采用平面繞避、橫穿的方法，當平面繞避對線路走向影響較大時，還可以對不良地質體和構造物進行空間立體化分析，通過設置橋隧采用立面繞避的方式，進而確定線路最優設計方案。

　　3.2 線路方案優化設計為繞避流動沙地、沙丘等不良地質，在保證工程設置安全性的前提下盡可能節約投資，利用獲取的數據基于 GIS 建立選線區域三維地理環境及不良地質表達模型，在此基礎上進行平、立面一體化的地質選線，多角度地尋求最優方案。設計的兩個工程方案如圖 15 所示，隧道方案對風沙危害采用立面避開設計，線路以短直方向設隧道穿越該埡口，出隧道后線路繼續足坡而下至比較終點;路基方案則以路基結構形式通過該沙害地段，從比較起點引出后向西以路基繞東柴山足坡而下至比較終點，結合沿線地理環境條件，在線路兩側的沙害區域采用“遠阻近固”的工程防沙體系。

　　4 結論

　　(1)根據鐵路數字化減災選線設計的目標，研究并建立了適用于風沙地區鐵路選線不良地質信息的內容組成及其識別、獲取方法，設計多源異構地質信息的組織框架與入庫規則，構建了基于GIS 的風沙地區鐵路選線不良地質信息庫。(2)采用矢量-柵格-專題一體化的風沙區地質環境建模方法，基于 GIS 的三維建模、空間分析功能，實現了不良地質信息在選線環境中的多維表達與空間分析，進而能夠在三維地質環境中進行線路方案的快速比選與優化設計。

　　(3)以格庫鐵路某沙害區選線為例，利用獲取的數據構建地質環境模型，通過空間分析風沙危害對不同線路方案的影響程度，得到較為合理的設計方案，表明基于 GIS 的風沙地區不良地質環境建模方法能比較準確地實現線路方案的優化設計，可為今后風沙地區鐵路選線提供一定的技術支持。

　　參考文獻：

　　[1] XIE Shengbo, QU Jianjun, HAN Qingjie, etal. Experimental definition and its significance on theminimum safe distance of blown sand between theproposed Qinghai-Tibet Expressway and theexisting Qinghai-Tibet Railway[J]. ScienceChina(Technological Sciences), 2020, 63(12): 2664-2676.

　　[2] 賈燕.荒漠區鐵路風沙防治關鍵技術[J].干旱區資源與環境，2021，35(12)：86-92.

　　[3] CHEN Boyu, CHENG Jianjun, XIN Lingui, et al.Effectiveness of hole plate-type sand barriers in reducingAeolian sediment flux: Evaluation of effect of hole size[J].Aeolian Research,2019,38：1-12.

　　[4] 易思蓉，聶良濤.基于虛擬地理環境的鐵路數字化選線設計系統[J].西南交通大學學報，2016，51(2)：373-380.

　　[5] 蒲浩，羅詩瀟，李偉，等.基于 OSG 的鐵路站場三維場景層次細節建模研究[J].鐵道工程學報，2017，34(3)：63-67.

　　[6] 呂希奎，陳進杰.鐵路數字化選線三維地質環境建模方法[J].鐵道學報，2015，37(8)：91-97

　　作者：高玉祥 1，程建軍 2，董曉峰 3，李澤宇 4