元謀干熱河谷強侵蝕區陷穴空間格局分析

何湘 舒成強 蔣良群

摘 要：陷穴是云南元謀干熱河谷強侵蝕區主要的潛蝕地貌之一，是一種可以導致水土流失的一種地質災害類型。選取元謀干熱河谷具有典型性的強侵蝕區——沙地村范圍內的陷穴作為研究目標，采取紙質版影像圖與電子圖相結合的方式進行野外數據采集，基于ArcGIS10.2軟件的空間分析工具，對陷穴的最鄰近指數、二階性質的K函數以及核密度進行定量研究，并對占據重要比重的串珠狀類型陷穴進行統計分析。研究結果顯示，沙地村餡穴總體上呈現為以233.72 m為半徑的范圍內的聚集分布模型，在此范圍內，陷穴的聚集模式呈現為先隨著半徑的增加而增加，當半徑為102.47 m時陷穴聚集程度達到峰值；之后隨著半徑進一步增加陷穴的聚集程度逐漸降低；當半徑超出233.72 m時，陷穴體現出分散的特征；總體而言，沙地村陷穴分為3個高密度區與3個低密度區，此類區域均為地形破碎的沖溝發育區；串珠狀陷穴占總陷穴數的58.7%，每串陷穴的陷穴數有3～6、9～11、14～15個，且多以3～5個為主，大多分布于沖溝內部。研究結果揭示了干熱河谷強侵蝕區陷穴在一定的范圍內具有較強的聚集性，往往與沖溝相伴生，多分布于沖溝中上游位置，特別是沖溝溝頭分布較為普遍。

關鍵詞：干熱河谷；潛蝕；陷穴；空間格局

DOI： 10.13651/j.cnki.fjnykj.2018.07.015

Abstract： Sink hole is one of the main underground erosion landforms in the intensive erosion zone of Yuanmou dryhot valley in Yunnan， and a type of geological disaster that can cause water and soil loss. The sink holes of Shadi Village in the typical erosion area in Yuanmou dryhot valley， were selected as the research objectives， the field data was collected by both paper data image and electronic image. Based on the spatial analysis tools in ArcGIS10.2 software， we quantitatively studied the nearest neighbor index， the K function of the secondorder property and the nuclear density of the sink holes， and statistically analyzed the beadlike sinking holes which occupying an important proportion. The results showed that the sink holes in Shadi Village generally presented aggregated distribution within the radius of233.72 m. Within this range， the aggregation degree of sinking holes appeared to increase firstly as the radius increasing， and the aggregated degree peaked at102.47 m； then decreased with the increase of the radius， and showed the characteristics of dispersion when the radius exceeded233.72 m. In general， the sink holes in Shadi Village could be divided into 3 highdensity areas and 3 lowdensity areas， all of which were gully development areas with broken terrain. The beadlike sinking holes accounted for58.7% of the total， and the number of holes per string was 3-6， 9-11 or 14-15， largely 3-5， which mostly distributed inside the gully. The results revealed that the sink holes in the intensive erosion zone of the dryhot valley was with strong clustering character in a certain range， often associated with the gullies， and distributed largely in the middle and upper parts of the gully， especially the gully head.

Key words： Dryhot valley； underground erosion； sink hole； spatial pattern

潛蝕地貌是由于地表水局部集中，沿垂直節理、劈理、裂隙或洞穴等進入地下進行潛蝕作用所形成的一類地貌[1-2]。由潛蝕作用所形成的土壤管道是一種與土壤有關的自然災害，影響著地下徑流系統[3]。土壤的潛蝕能破壞各類工程設施，且破壞速度達到驚人的程度[4]。潛蝕作用在不同類型的氣候、地貌和土壤環境中都有分布，在山坡、沖溝以及耕地等微地貌上都存在[5-7]，并導致嚴重的水土流失。在某些地區潛蝕對于泥沙的貢獻可以達到3%～61%[8]，如在英國有近30%的土地是受到了潛蝕作用的影響[9]。

陷穴屬于潛蝕地貌類型之一，是特定地層條件下因浸水、積水而發生局部濕陷或各種管道在水流等侵蝕作用下進一步擴容、塌陷而形成的露出于地面的豎向洞穴，可以影響土壤的侵蝕，導致土地破壞，加速溝谷系統的發育[10-11]，是一種嚴重的地質災害類型[12]。云南元謀干熱河谷位于長江中上游重要的生態屏障保護區，但干熱河谷內溝壑縱橫，水土流失極為嚴重，特別是在強侵蝕區，潛蝕地貌分布極為廣泛，是區域溝谷發育、地形破壞的重要影響因素之一，但目前還不清楚此類地貌在侵蝕區的分布規律。因此，有必要對潛蝕地貌在侵蝕區的分布格局進行深入研究。本研究選取元謀干熱河谷具有代表性的強侵蝕區——沙地村的陷穴作為研究對象，通過野外調查結合遙感影像資料，利用ArcGIS10.2軟件空間分析工具分析了沙地村陷穴空間分布格局，以期能找出陷穴在強侵蝕區的分布規律及其原因，為干熱河谷水土流失治理及生態恢復提供理論依據。

1 研究區概況

元謀干熱河谷位于滇中高原北面的龍川江河谷盆地內（E 101°35′至E 102°25′，N 25°25′至N 26°07′），全縣國土面積2021.69 km2。元謀干熱河谷四周高、中間低，四周高山圍繞所形成的一個典型河谷盆地，境內龍川江自南往北縱穿盆地，盆地內高差懸殊，整個盆地自南向北呈傾斜的“筲箕凹”形[13]。元謀區內最高處海拔2835.9 m，最低處海拔898 m，高差達到1937.9 m，相對高差大造成自然環境要素在垂直方向上分異明顯，從海拔最低處到海拔最高處可分為：寬谷平壩區、壩周低山丘陵區、中低山區和高中山區4個垂直地帶[14]。由于橫斷山脈的深度切割，加之河谷焚風效應顯著，氣候干熱，具有明顯干濕季交替的亞熱帶季風氣候型。年均降水量較少，僅615.1 mm，并且時間分配不均，一半以上的降水集中在6月至10月，蒸發量高達3640.5 mm，年平均氣溫21.8℃，年均日照時數2647 h，正常年無霜，被譽為“天然溫室”[15-16]。

本研究區域位于干熱河谷中部，元謀縣城附近的南東方向，該區域靠近人類聚居區，四周為人類活動較為頻繁的農業種植區。研究區地形崎嶇，沖溝系統發育不完善，是元謀干熱河谷水土流失較為嚴重的強侵蝕區之一。

2 材料與方法

2.1 數據來源

2.1.1 影像資料 遙感影像主要來自于P星數據和Bigemap下載的衛圖數據兩部分。兩者可以形成互補，有利于野外調查時和實際情況對比。

2.1.2 陷穴空間數據采集 以Bigemap下載的19級衛圖（分辨率為0.51 m）為基礎圖件，采用華測LT600手持GPS平板電腦（亞米級）采集陷穴空間數據。在實際操作中，為了提高野外調查工作效率，并保證調查數據的精度，野外數據采集采用紙質版影像圖與電子圖相結合的方法。首先在野外找到至少3處具有明顯參考性的地物標志，在紙質圖中予以標注，同時在電子圖中找到該類標志物予以標注坐標，另外利用GPS平板自帶坐標系統進行該標志物的坐標標識，在Mapcloud office軟件中對兩組坐標進行配準校正，最后將校正后的影像圖導入平板電腦。處理后的圖像坐標系統與實際地理坐標具有極好的吻合度，在提高野外工作效率的同時也保證了數據采集的精度。

2.2 研究方法

研究區域選擇在元謀縣附近的強侵蝕區——沙地村范圍內。自2010-2018年對該地侵蝕地貌進行連續調查，為獲得最新的陷穴空間分布格局，本研究采用的數據為2017-2018年最新調查的數據，包括野外數據調查和室內數據處理兩部分：野外重點調查陷穴的空間分布特征，同時記錄陷穴之間的距離以及陷穴的平面形態特征。數據處理以ArcGIS10.2為工具，將陷穴看成是空間上一系列的點源數據，采用最鄰近指數、二階K函數以及核密度分析來分析判別陷穴的空間格局。

3 結果與分析

3.1 陷穴空間分布格局

3.1.1 陷穴的聚集分布 本研究以野外采集的252組陷穴坐標數據為基礎，利用Arc GIS10.2的空間統計中最鄰近指數工具計算其最鄰近指數。從陷穴最鄰近指數計算結果（圖1）可知，陷穴間最鄰距離最小值為4.67，最大值為140.85，平均值為36.34。最鄰近指數為0.24<1，在空間上為聚集模式。對顯著性進行檢驗，得到的Z值為-23.51倍標準差。根據王遠飛的描述[17]，若Z值的符號為負，則模式趨向于聚集；若符號為正，則趨向于均勻分布。可知陷穴總體上是呈現聚集分布的模式。

3.1.2 不同尺度有不同的聚集模式 采用二階性質的K函數方對陷穴的分布格局進行深入分析。研究結果表明，陷穴在以233.72 m為半徑的范圍內呈聚集分布。在此范圍內，陷穴的聚集模式呈現為先隨著半徑的增加聚集程度逐漸增加，在半徑為102.47 m時陷穴聚集程度達到峰值，隨著半徑進一步增加陷穴的聚集程度逐漸降低；當半徑超出233.72 m時，陷穴體現出分散的特征（圖2）。

3.1.3 陷穴的核密度分析 為了進一步分析陷穴分布的位置，采用核密度分析法計算了沙地村陷穴密集區域。研究結果顯示，沙地村陷穴分為3處高分布區和3處次低分布區（圖3）。其中，A1高密度區屬于風化殘余臺地，該臺地上植被稀疏，零星分布著小葉桉，地表土壤為燥紅土，開裂現象明顯，該臺地邊緣發育著處于幼年期的沖溝系統，溝內陷穴多呈密集的串珠狀分布。A2陷穴高分布區地形較為破碎，臺地多為處于幼年期和青年期的沖溝所切割，沖溝下切和側蝕十分強烈。該處的沖溝有一個共同的特點：一側溝岸坡度較緩，另一側大致為直立形式。陷穴基本發育在直立溝岸的正下方，呈密集的串珠式分布。沖溝四周的植被稀少，零散分布有一些禾草類植物。此處的串珠陷穴之間大都已經貫通，中間僅剩1個土橋，土橋橋面寬度10～40 cm。A3陷穴高分布區中間為臺地，臺地四周為發育短小沖溝所環繞。陷穴主要散落在短小的沖溝中及部分集中在一個小平臺上，平臺上的陷穴多為種植的灌木枯死遺留的樹坑發育而來。

此外還有3個比較集中的次低分布區，分別為圖3中的B1、B2、B3三處。B1和B2處分布陷穴數量不多，相對較為零散；B3處及其周圍區域較之B1和B2處，陷穴數量多，分布范圍廣，在核心區相對集中。陷穴次低分布區地貌類型主要為發展處于穩定期的沖溝。這類沖溝形態上體現為溝道寬闊，以淺 “U” 型溝為主， 溝底多為沖積物所覆蓋，溝底多呈平坦狀，低矮草本植被生長良好。

3.2 串珠狀陷穴空間分布格局

串珠狀陷穴一般由幾個連續的陷穴組成（圖4），在研究區占有較大的比重。沙地村的串珠狀陷穴主要有以下幾種類型：3個1串、4個1串、5個1串、6個1串、9個1串、10個1串、11個1串、14個1串、15個1串，每種類型都有一定的數量并占總陷穴數量一定比例。從調查結果（表1）可知，串珠狀陷穴總數共有148個，占總陷穴數量（252個）的58.7 %。其中以3個、4個、5個1串為主要串珠組，共有20個，占總串珠組（27個）的74.0%；3～5個串珠狀陷穴共有77個，占串珠陷穴總數（148個）的52.0%。3個、4個1串的陷穴主要分布于元馬鎮翠峰小學后方的沖溝群中，即次低分布區的B3區域及其周圍；而5個以上成1串的陷穴多分布于A1區和A2區周圍區域（圖4）。這類區域的沖溝多數比較短小，溝內植被覆蓋度低，多為裸地，側蝕和下切強烈，沖溝兩岸的直線距離小。

4 結論與討論

（1）陷穴總體上呈現為以233.72 m為半徑的范圍內的聚集分布模型。在此范圍內，陷穴的聚集模式呈現為先隨著半徑的增加而增加，當半徑為102.47 m時陷穴聚集程度達到峰值；之后，隨著半徑進一步增加陷穴的聚集程度逐漸降低；當半徑超出233.72 m時，陷穴體現出分散的特征。

（2）陷穴核密度分析結果顯示，沙地村陷穴分為3處高密度區和3處次低分布區。高密度分布區地形總體特征為風化殘余臺地、地形破碎嚴重，陷穴主要分布于沖溝系統內部，另一部分由于發育于枯死灌木遺留洞坑而分布與殘留臺地。次低分布區陷穴主要分布于發育穩定的沖溝系統內，該類溝谷以淺“U”型溝為主，溝底多為沖積物所覆蓋，溝底多呈平坦狀，低矮草本植被生長良好。

（3）陷穴在元謀干熱河谷分布范圍廣，在沖溝發育的強侵蝕區有著廣泛的分布；在非強侵蝕區，諸如工礦分布區、農業耕作區等也都有分布。研究陷穴在更大范圍內的空間分布格局，弄清陷穴空間分布格局的影響因素對于在干熱河谷有效開展社會經濟活動及對于金沙江流域生態環境保護具有較強的理論指導意義。關于陷穴在更大范圍內的空間分布格局有待于下一步研究。

（4）快速獲取陷穴的分布數據是高效開展研究的必須手段。本研究陷穴形態尺度較小，難以采用遙感影像獲取精確的位置信息，而地面調查又受限于人力及地形可達性影響，不能高效、全面地完成野外數據調查。因此，未來采用小型無人機結合地面測量是解決目前難以高效采集野外數據的一個有效手段。

參考文獻：

[1]鄧青春，張斌，羅君，等.元謀干熱河谷潛蝕地貌的類型及形成條件[J].干旱區資源與環境， 2014，28（8）：138-144.

[2]王景明，王君.冀中南黃土潛蝕地貌與黃土構造節理[J].地理研究，1994，13（1）：90-93.

[3]BHAGYALEKSHMI S，MOHAN M，SREEDHARAN K，et al.Evaluation of the factors controlling the spatial distribution of soil piping：a case study from the southern Western Ghats，India[J].Arabian Journal of Geosciences，2015，8（10）：8055-8067.

[4]BAILLIE I C，FAULKNER P H，ESPIN G D，et al.Problems of Protection Against Piping and Surface Erosion in Central Tunisia [J].Environmental Conservation，1986，13（1）： 27-32.

[5]郭玉紅，郎南軍，和麗萍，等.元謀干熱河谷8種植被類型的林地土壤特性研究[J].西部林業科學，2007，36（2）：56-64.

[6]孫長安，王煒煒，董磊，等.我國植被恢復對土壤性狀影響研究綜述[J].長江科學院院報，2008，25（3）：6-8，22.

[7]JONES J A A.Soil piping and its hydrogeomorphic function[J].cuaternarioy Geomorfologia，1994，8（3）：77-102.

[8]SAYER A M，WALSH R P D，CLARKE M A，et al.The role of pipe erosion and slopewash in sediment redistribution in small rainforest catchments，Sabah，Malaysia[J].IAHSAISH Publication，2004，288：29-36.

[9]JONES J A A.Implications of natural soil piping for basin management in upland Britain [J].Land Degradation & Development，2004，15（3）：325-349.

[10]ZHANG T，WILSON G V.Spatial distribution of pipecollapses in Goodwin Creek Watershed，Mississippi[J].Hydrological Processes，2013，27（14）：2032-2040.

[11]VERACHTERT E，EECKHAUT M V D，POESEN J，et al.Factors controlling the spatial distribution of soil piping erosion on loessderived soils：A case study from central Belgium[J].Geomorphology，2010，118（3）：339-348.

[12]ARGENTIERI A，CARLUCCIO R，CECCHINI F，et al.Early stage sinkhole formation in the Acque Albule basin of central Italy from geophysical and geochemical observations[J].Engineering Geology，2015，191：36-47.

[13]起樹華，王建彬.元謀干熱河谷氣候生態環境變化的初步分析[J].氣象研究與應用，2007，28（S2）：125-127.

[14]拜得珍，紀中華，廖承飛，等.淺析干熱河谷水土保持型生態農業土壤工程技術——以云南元謀為例[J].西南農業學報，2004，17（S1）：272-275.

[15]王小丹，鐘祥浩，范建容，等.金沙江干熱河谷元謀盆地沖溝溝頭形態學特征研究[J].地理科學，2005，25（1）：63-67.

[16]史凱，張斌，艾南山，等.元謀干熱河谷近50a降水量時間序列的DFA分析[J].山地學報，2008，26（5）：553-559.

[17]王遠飛.空間數據分析方法[M].北京：科學出版社，2007.

（責任編輯：林玲娜）