基于GIS的公路橫斷面地溫時空分析方法研究——以青藏公路為例

張 興 國，劉 學 軍，陳 建 兵

（1.南京師范大學虛擬地理環境教育部重點實驗室，江蘇 南京 210046;2.信陽師范學院城市與環境科學學院，河南 信陽 464000;3.中交第一公路勘察設計研究院有限公司，陜西 西安 710075）

基于GIS的公路橫斷面地溫時空分析方法研究
——以青藏公路為例

張 興 國1，2，劉 學 軍1，陳 建 兵3

（1.南京師范大學虛擬地理環境教育部重點實驗室，江蘇 南京 210046;2.信陽師范學院城市與環境科學學院，河南 信陽 464000;3.中交第一公路勘察設計研究院有限公司，陜西 西安 710075）

研究青藏公路凍土區公路橫斷面地溫時空分析方法，獲取地溫時空變化規律，為公路結構設計、維護及保溫材料的選取提供依據。傳統的分析方法主要通過繪制地溫折線圖研究地溫隨時間或深度的變化規律，方法簡單，結果易于理解和應用，但較難發現地溫的時空變化規律。該文提出基于GIS的公路橫斷面地溫時空分析方法。選取西大灘實驗段編號為K2887＋500的橫斷面為實驗靶場，通過建立地溫場圖層、多時段融區時空動態分析等方法，得出該橫斷面地溫經歷了融區收縮、融區擴展及融區過渡3個階段，并準確計算出各觀測時間的融區范圍及融區變化速度，實現了橫斷面地溫時空量化分析，其思路對時態GIS相關開發及應用具有參考價值。

公路橫斷面;時空分析;青藏公路;地溫

0 引言

青藏公路格爾木至拉薩段位于自然環境惡劣、地質條件復雜的青藏高原腹地，穿越里程達624 km的大片連續多年凍土、島狀多年凍土區［1，2］。修筑在高原多年凍土區的路基，是以下伏多年凍土為依托，多年凍土又以負溫為基本生存條件［3］。多年凍土地溫不但代表凍土的物理學特征，而且在實際工程建設中具有重要意義。青藏公路下伏凍土的研究，特別是凍土地溫的研究，對道路病害治理、道路結構設計、保溫材料的選取等具有重要意義。

對于凍土地溫的分析，傳統的方法主要借助地溫變化曲線圖分析地溫隨路基深度及時間的變化規律［4，5］，方法簡單、易于理解，但其對于地溫在空間上的變化，乃至時空動態變化難以實現。在地溫場分析方面，研究者通常基于非GIS技術進行空間模擬［6］，效果明顯，但該方法的研究數據通常無空間參考，僅為模擬而設置，難以與已有地理數據（如DEM、公路里程數據等）集成，從而造成分析成果無法空間量測，也無法與其他地層數據關聯;而且，在時空動態分析方面缺乏有力支持。

本研究提出基于GIS的公路橫斷面地溫時空分析技術。選取西大灘實驗段編號為K2887＋500的橫斷面為實驗靶場，通過建立地溫場圖層、多時段融區時空動態分析等過程，確定了該斷面地溫經歷了融區收縮、融區擴展及融區過渡3個階段，準確計算出各階段的時間界限、融區范圍和融區變化速度，并將該方法集成至“青藏公路基礎信息平臺”。

1 公路橫斷面地溫觀測點分布

目前，在青藏公路高原多年凍土區，根據不同的觀測目的，設置有11個實驗段，分別用于不同的監測目的，如西大灘段主要用于觀測多年凍土區北界地溫，其他還有昆侖山段、楚瑪爾河段、五道梁段、北麓河段、雅瑪爾河段、沱沱河段、開心嶺段、唐古拉山北坡、唐古拉山南坡和納金河段。每個實驗段上設置1個或多個橫斷面，在每個橫斷面上設置多個地溫孔，根據地溫孔距道路中心線的相對位置，可劃分為左路肩、右路肩、路中、左天然（道路之外的一般位置）和右天然幾種類型［7］。對每個地溫孔，從地表開始向下，每間隔一段距離（如0.3 m等）安置一個測溫裝置（圖1）。這樣對于每個橫斷面對應有許多測溫點。在實際工作中，考慮到道路兩側地溫分布的相似性（除個別橫斷面道路兩側需同時監測），為了節省成本，可以僅在道路一側進行，而分析時可以假設另一側的測溫點具有相同的分布。

為了對地溫孔的觀測數據進行管理，本研究將地溫時空觀測數據以二維表的形式組織，每個地溫孔對應一張二維表格。表格的第一行除了第一列，均記錄當前地溫孔所屬的測溫點距離地表的深度，第二行以下的各行，其第一列代表時間，其他各列代表地溫觀測值（表1）。對每個橫斷面都關聯有對應的地溫孔或測溫點，對其上所有點的每個觀測時間，都可以基于表1動態生成該時段的點圖層，并且點的屬性為該時段的地溫值。

圖1 橫斷面上測溫點分布Fig.1 The distribution diagram of measurement point on cross section

表1 某地溫孔地溫Table 1 The temperature of the ground temperature hole

2 公路橫斷面時空分析方法及其實現技術

2.1 傳統的地溫時空分析方法

傳統的方法是通過折線圖分析地溫的時空變化規律，主要包括地溫孔－深度分析、測溫點－時空分析等。

地溫孔－深度分析是以地溫孔深度作為Y軸、地溫作為X軸，尋求地溫隨深度的變化規律;對任一地溫孔，從地表向下對應有多個測溫點。對任一觀測時間可以繪制出地溫隨深度的變化曲線，基于該曲線，可以得出該觀測時間地溫在不同深度的變化規律，也可以計算出0℃地溫對應的深度。本文對西大灘實驗段的K2887＋500橫斷面進行分析，該橫斷面有3個地溫孔（DWK01、DWK02及DWK03），對DWK01在2004年的4個時間點進行地溫孔－深度變化分析（圖2）。通過該曲線可以看出DWK01隨著深度的增加，地溫在距地表5.8 m處趨于穩定，為－1℃;在距地表2.5 m處，地溫達到0℃（有些月份在0℃以下）。通過多個時間的地溫孔－深度時空分析，可以得出一系列0℃地溫對應的深度，找出其中深度的最大值，即凍土上限，這對于凍土地區工程建設具有重要意義。

測溫點－時間分析是以測溫點的觀測時間為X軸、溫度為Y軸，尋求地溫隨時間的變化規律。對DWK01的多個測溫點進行分析，其深度分別是0.3 m、1.3 m、2.8 m 和3.8 m（圖3）。從圖3可以得出，深度為0.3 m和1.3 m處地溫隨時間變化明顯，并且在8月份左右地溫達到最高，為4℃以上;而深度大于2.8 m時，地溫隨時間變化較小，通常低于0℃，并且在3－4月溫度最低，為－1℃以下。

2.2 基于GIS的地溫時空分析

每個橫斷面在某一觀測時間都對應一個測溫點圖層。為了獲取地溫的時空變化規律，可以基于GIS技術，對測溫點圖層根據地溫觀測值進行內插，從而模擬出地溫場。這樣，一個橫斷面依據觀測時間對應有多個地溫場，通過對多個地溫場的對比，尋求地溫的時空變化規律。基于地溫場可以生成等溫線圖，特別是0℃等溫線的分布至關重要。0℃等溫線的深度曲線即凍土上限的分布曲線。一年四季，年間交替，0℃等溫線會隨著時間動態變化。本文以西大灘實驗段的K2887＋500橫斷面為例，進行地溫的時空分析。為了分析方便，需要對天然孔的位置進行水平或垂直移動，以求內插能夠順利進行，也滿足實際的地溫分布情況。同時，根據實際的路基橫斷面圖，獲取內插的范圍。最終觀測點的布局如圖4，通過反距離權重（IDW）進行內插，并獲取0℃等溫線（白色的曲線），如圖5。可以看出，2004年10月6日0℃等溫線總體呈拋物線形狀，最低點距地表8.3 m左右，而兩側距地表3.4 m左右。該等溫線以上均大于0℃，以下均小于0℃。

為了獲取地溫的時空變化規律，對K2887＋500橫斷面上2004年1－12月的大量觀測數據（共觀測22次，每月上旬和下旬各一次）進行地溫場及等溫線分析：1）對每一觀測時間的觀測點構建地溫場;2）基于地溫場提取融區（地溫大于0℃的區域）的分布;3）對所有觀測時間的融區進行對比分析。

實驗表明，2004年12月下旬至3月下旬，融區呈橢圓形，其外地溫均低于0℃。上側距地表距離為1.8～4.8 m，下移速度為0.5 m/半月;下側距地表距離為6.5～8.4 m，上移速度先慢后快（0.2～0.6 m/半月）;左右側距邊距3.9～6.6 m，水平內縮速度約0.2 m/半月。隨著天氣回暖，橢圓融區逐漸消失，在3月下旬幾乎消失殆盡，接近無融區狀態。該時段分析結果如圖6，圖中的封閉曲線為其對應時間的融區邊界，這里定義該階段為融區收縮階段;可以發現12月17日－3月19日，融區范圍逐漸變小，并且上邊變化較大，最后完整的橢圓被瓦解，即融區幾近消失。

2004年4月上旬開始，以道路中心為起點，融區呈半圓形向下、向左右擴展，呈輻射狀。路中向下擴展先快后慢（1.1～0.3 m/半月），至10月下旬到達極限，路中融區最大深度達8.7 m，天然地表凍土上限僅為3.2 m，定義該階段為融區擴張階段（圖7）。從圖7中可以看出，路中受公路影響比天然地表受公路影響多3 m左右。

11月上旬－12月上旬，輻射狀融區逐漸縮小，最后趨于橢圓形，定義該階段為過渡階段，即由擴張轉為收縮。11月上旬天氣開始變冷，地表氣溫降低，地表部分立即呈現非融區狀態，而地下融區的變化則是緩慢的，范圍逐漸縮小。

圖6 融區收縮階段變化Fig.6 The diagram of melt area shrinking

圖7 融區擴張階段變化Fig.7 The diagram of melt area expansion

通過以上分析，可以得出公路路基下凍土融區經歷了收縮、擴展和過渡3個階段。路中凍土影響最大深度達8.7 m，而天然孔處僅為3.2 m，凍土融區變化速度為0.3～1.1 m/半月。該成果明確顯示凍土受氣候變化和道路交通影響較大，特別是在10月下旬達到極限。為了維持凍土區青藏公路的正常運轉，減少公路病害，必須采取措施降低凍土上限值。目前，盡管采取了提高路基、鋪設保溫材料、調整路基結構、改變路基材質等措施，但效果還不理想，仍處于研究階段［8，9］。本研究基于GIS技術進行分析，可以從空間上直觀發現融區的動態變化，從定性和定量兩方面對路基橫斷面的地溫進行時空分析，比傳統的統計圖表分析更具優越性，為各種治理措施的效果評價提供了重要信息。

2.3 地溫時空分析系統實現技術

為了動態分析各實驗段的橫斷面地溫，基于Oracle、VS C＃ 2005、ArcEngine9.2、Google Earth等技術開發了青藏公路基礎信息平臺。本系統對相關的大量時空數據進行統一管理，并實現了傳統的基于折線圖的地溫時空分析和基于GIS的地溫時空分析，將地溫時空分析與各專題空間數據庫關聯，實現了公路地表、地面及地下全方位聯動，二維、三維、實景影像及視音頻等多視角觀察。同時，系統提供強大的數據輸出功能，極大地方便了科研人員的工作，提高了效率。

3 結論與討論

青藏公路橫斷面地溫時空分析涉及大量的地溫觀測數據和路基結構數據，傳統的管理技術和分析方法不便于地溫的時空模擬，難以與基礎地理數據關聯，造成相關圖件量測精度無法保障，數據多源異構，共享困難;基于GIS技術，充分發揮其數據管理、分析及強大的可視化能力，極大地方便了凍土的時空分析。本文所采用的地溫場、融區時空分析等方法，得出K2887＋500橫斷面地溫經歷了融區收縮、融區擴展及融區過渡3個階段，并對一些重要的參數（如凍土上限、收縮速度、擴展速度等）進行了量化，為青藏公路凍土分析及工程建設提供了重要依據。該方法可在所有實驗段通用，具有較強的實用價值。

該方法仍需完善和深化：1）加強時空數據模型的研究。時空數據模型是時空數據庫的核心部分，它保證時空數據庫的完整性，也決定著對實體動態信息的具體操作［10］。在后續研究中應尋求適合該領域的時空數據模型，為地溫數據管理和分析提供支撐。2）加強時空變化的可視化研究。時空變遷的可視化水平直接影響到時空規律的發現，可基于動畫技術，研究動畫過渡算法。3）本研究僅僅是對已有數據的變化分析，缺乏智能化預測功能。如為了保護路基下的凍土，其中一個辦法就是提升路基高度，那么應該提升多高，路基的詳細參數是多少，遵循什么樣的規律，這些都有待進一步研究。另外，影響地溫的因素除氣溫和深度外，其他因素如海拔、土壤、地質情況等尚未在該方法中體現，需加強多因素綜合的GIS分析技術研究。

［1］陳建兵，汪雙杰，章金釗，等.青藏公路高路基病害的形成及其機理［J］.長安大學學報（自然科學版），2008，22（6）：30－35.

［2］竇明健，胡長順，何子文，等.青藏公路多年凍土段路基病害分布規律［J］.冰川凍土，2002，24（6）：780－784.

［3］章金釗，霍明，陳建兵.多年凍土地區公路路基穩定性技術問題與對策［M］.北京：人民交通出版社，2008.24.

［4］智鵬，馮云梅，蔣秋華，等.凍土地溫自動檢測系統的設計與實現［J］.鐵道學報，2009，31（1）：59－62.

［5］潘衛東，余紹水，賈海鋒，等.青藏鐵路沿線多年凍土區地溫場變化規律［J］.冰川凍土，2002，24（6）：774－779.

［6］吳志偉，宋漢周.淺層地溫場中熱對流數值模擬［J］.巖土力學，2010，31（4）：1303－1308.

［7］汪雙杰，陳建兵，黃曉明.熱棒路基降溫效應的數值模擬［J］.交通運輸工程學報，2005，5（3）：41－46.

［8］汪雙杰，陳建兵，章金釗.保溫護道對凍土路基地溫特征的影響［J］.中國公路學報，2006，19（1）：12－16.

［9］汪雙杰，陳建兵，黃曉明.凍土路基護道地溫特征研究［J］.巖石力學與工程學報，2006，25（1）：146－151.

［10］袁一泓，高勇.面向對象的時空數據模型及其實現技術［J］.地理與地理信息科學，2008，24（3）：41－44.

Study on Spatial and Temporal Analysis Methods for the Ground Temperature of Highway Cross-Sectional Area Based on GIS：A Case Study of Qinghai-Tibet Highway

ZHANG Xing－guo1，2，LIU Xue－jun1，CHEN Jian－bing3
（1.KeyLaboratoryofVirtualGeographicEnvironment，NanjingNormalUniversity，MinistryofEducation，Nanjing210046;2.CollegeofUrbanandEnvironmentalScience，XinyangNormalUniversity，Xinyang464000;3.ChineseFirstInstituteLimitedCorporationofHighwaySurveyandDesign，Xi′an710075，China）

Ground temperature study of Qinghai-Tibet highway has been very important，which can provide key information for the structure design and the selection of insulation materials.The traditional method is based on the temperature line graph to study the ground temperature changes with time and depth，but it is more difficult to find the spatial and temporal change.On the basis of current methods，a new method which is based GIS technology is introduced.Through the establishment of geothermal field and multi-time spatial and temporal analysis，the ground temperature change of highway cross-sectional area becomes easier to obtain.The test for K2887＋500 of Xidatan test segment shows the change is divided into three stages（melt area shrinking，melt expansion and melt transition）and the melt range and speed can also be calculated accurately.The method has an important reference value for the development and application of spatial and temporal GIS.

road cross section;spatial and temporal analysis;Qinghai-Tibet highway;ground temperature

P208

A

1672－0504（2011）05－0021－04

2011－04- 21;

2011－07－10

國家自然科學基金項目（40971230）;河南省教育廳自然科學研究計劃項目（2011B170010）

張興國（1979－），男，講師，博士研究生，主要從事GIS開發與應用方面的研究。E－mail：xingguozhang＠tom.com