GF-1图像在中印边境楚鲁松杰村地质灾害调查中的应用
张焜1,2, 李晓民1,2, 马世斌1,2, 刘世英1,2, 李生辉1,2
1.青海省青藏高原北部地质过程与矿产资源重点实验室,西宁 810012
2.青海省地质调查院,西宁 810012

第一作者简介: 张 焜(1973-),男,高级工程师,主要从事地质矿产遥感技术应用研究。Email:zhangkun0623@sina.com

摘要

为了解高分一号卫星(GF-1)图像的应用性能和查明中印边境楚鲁松杰村一带的地质灾害发育规律,利用DEM和GF-1图像对该地区的岩性、构造和地质灾害进行详细的遥感解译。研究结果表明,GF-1数据的融合图像对灾害体的微地貌形态反映较明显,能够用于识别地质灾害体的位置、边界范围等特征信息; 区内地质灾害类型主要为滑坡和寒冻风化碎屑流等,滑坡、崩塌、泥石流多发生在[25°,35°]坡度段内,而寒冻风化碎屑流的发生则与坡度的增大成正相关; 地质灾害运移方向与主构造线方向、基岩层面倾向基本一致。研究还揭示出该区内分布多级河流阶地,较大的垂直落差及丰富的第四系松散堆积物是地质灾害发育的重要因素。

关键词: 高分一号卫星(GF-1); 中印边境; 地质灾害; 寒冻风化碎屑流; 遥感调查
中图分类号:TP 79 文献标志码:A 文章编号:1001-070X(2016)02-0139-10 doi: 10.6046/gtzyyg.2016.02.22
Application of GF-1 image to geological disaster survey in Cosibsumgy village on Sino-India border area
ZHANG Kun1,2, LI Xiaomin1,2, MA Shibin1,2, LIU Shiying1,2, Li Shenghui1,2
1. Qinghai Province North of the Qinghai-Tibet Plateau Geological Process and Mineral Resources Key Laboratories, Xining 810012,China
2. Institute of Geological Survey of Qinghai Province, Xining 810012, China
Abstract

Using data obtained from DEM and GF-1 Satellite(called GF-1) concerning lithological characteristics in Cosibsumgy village on Sino-India border area, the authors conducted detailed remote sensing interpretation of structures and geological disasters to understand the application performances of GF-1 and find out the patterns of geological disasters that have happened in this area. The results show that the fused image of GF-1 data appears to be evident for the micro morphology of the disasters. It can identify the locations and outlines of the geological disasters. Landslide and frost weathering debris flow are dominant in the geological disasters. The landslide, collapse and mudslide mostly occur in the gradients ranging [25°, 35°], whereas the frost weathering debris flow has positive correlation with the gradients. The migration directions of the geological disasters are basically consistent with the directions of the main tectonic lines and the tendency of the bedrocks. Meanwhile, the multilevel river terraces existing in the area and the huge vertical fall from Quaternary system are the main causes for geological disasters, as revealed in this study.

Keyword: GF-1 satellite; Sino-India border; geological disaster; frost weathering debris flow; remote sensing interpretation
0 引言

西藏自治区札达县楚鲁松杰村一带地处中印边境地区。在该区开展基础地质、矿产、地貌、水文地质、工程地质、交通设施与重要地物等综合遥感调查, 对该区的社会经济建设、生态环境保护等有着重要意义, 尤其对边境地区的国防建设具有重要的战略意义。由于特殊的地理、地质环境及气候条件, 该区域内的滑坡、泥石流、崩塌及滚石等地质灾害频发, 常常阻断交通; 曾有地质灾害波及到印度境内帕里河流域并影响到当地的安全, 致使印度反映强烈[1, 2]。如何减少地质灾害对该区域内交通和国防建设带来的危害, 是亟待研究的一个课题。

遥感技术能快速、高效地获取地表信息, 在滑坡、崩塌、泥石流等地质灾害调查、监测和研究中发挥了重要作用 [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]。随着我国国产卫星技术的快速发展, 国产卫星在资源、环境、测绘、海洋等领域的应用越来越重要[13, 14, 15, 16]。2013年4月26日在我国酒泉卫星发射中心成功发射的高分一号卫星(GF-1)是中国高分辨率对地观测系统的首发星, 对提高我国高分辨率卫星数据(以下简称高分数据)自给率具有重大意义。本文以西藏自治区中印边境楚鲁松杰村一带的地质灾害为研究对象, 采用GF-1卫星数据, 对其全色、多光谱数据进行融合, 制作正射遥感影像图; 并将遥感数据与DEM数据叠加生成三维遥感影像图, 进行地质灾害解译。从应用角度了解GF-1卫星数据特色及其在地质灾害调查方面的能力; 在分析研究区的自然与地质环境、水文地质条件的基础上, 查明研究区地质灾害的发生地点、种类、规模和范围, 探讨该区域地质灾害形成和发育的规律。

1 研究区概况

研究区位于西藏阿里地区西部与印度接壤地带, 属喜马拉雅山脉, 面积约800 km2。其中, 境内行政区划属西藏自治区阿里地区扎达县, 面积约400 km2; 境外面积约400 km2, 属印度西姆拉地区。区内为典型的高山峡谷地貌, 丘状高原及构造侵蚀形成的深切峡谷为其主要地貌特征。大多数山峰海拔超过5 000 m, 现代河流发育的高程为3 000~3 800 m; 相对高差多在1 500 m以上, “ V” 形沟谷发育。山峰上发育典型的山岳冰川, 具有各种类型的冰川地貌。

研究区位于帕米尔构造结东翼, 属北喜马拉雅特提斯沉积带, 主要由元古宇基底岩系聂拉木岩群(AnZN)、古生界盖层浅变质岩系色龙群(PS)和中生界才里群(JC)台地碳酸盐沉积组成, 并伴有喜山早期含电气石二云母花岗岩的侵位。聂拉木岩群是构成喜马拉雅地区结晶基底的主体, 整体为一套火山-沉积建造; 色龙群为一套浅变质岩系, 由灰绿色砂板岩和变砂岩组成, 夹少量碳酸盐沉积, 主要为浅海-滨海相陆源碎屑沉积建造; 才里群则为一套台地相碳酸盐建造。构造变形以面理改造和小型褶皱发育为特征, 叠加有后期脆性断裂改造; 断裂构造以NW向多期大型压剪为主要特征, NE向和近SN向断裂多属主压断裂所派生 [17, 18] 。从所处的特殊地质构造、气候背景条件分析, 研究区属滑坡、崩塌、泥石流等地质灾害易发区。

2 遥感技术方法
2.1 GF-1卫星数据

GF-1卫星搭载了2台2 m分辨率全色/8 m分辨率多光谱相机和4台16 m分辨率多光谱相机[19], 其技术指标见表1

表1 GF-1卫星有效载荷技术指标 Tab.1 Payload parameters from GF-1 satellite

本文使用的GF-1原始数据的获取时间为2014年6月26日, 包括8 m分辨率的多光谱数据和2 m分辨率的全色数据, 均自带有理函数模型(rational function model, RFM)的元数据文件。图像中云雪覆盖量较少(约占1%), 辐射畸变、几何畸变和噪声均较小, 图像清晰, 整体质量较好, 适用于进行遥感地质解译。

2.2 图像处理

在遥感专业处理软件Erdas中的LSP模块下采用 “ DEM+RFM+GCP” (“ 数字高程模型+有理函数模型+地面控制点” )的方法分别对GF-1多光谱和全色原始数据进行正射纠正。其中, 读取RPC.XML文件中的内、外定向参数作为GF-1的RFM模型参数; 选用由IRS-P5卫星数据制作的数字高程模型(DEM)数据作为本次正射纠正处理的辅助数据; 选择GF-1多光谱图像的B4(R)B3(G)B2(B)波段组合并采用Pansharp算法与全色波段进行融合处理; 最后根据本文研究对象对融合后的图像进行局部拉伸处理, 制作遥感正射影像图(图1)。

图1 GF-1图像处理流程Fig.1 Flow chart of GF-1 image processing

2.3 三维可视化

利用ArcInfo-ArcScene软件平台, 以制作好的遥感正射影像图为基础, 按照一定的原则, 通过叠加研究区DEM数据, 实现对研究区三维地形特征与地理要素的立体观察。

2.4 遥感解译

在充分收集已有地质环境与地质灾害资料的基础上, 建立地质灾害及其发育地质环境的遥感解译标志; 利用Erdas, PhotoShop和ArcGIS等软件平台, 以人机交互解译为主要手段, 获取相关地质背景和地质灾害信息。

地质背景解译主要参考研究区已有基础地质资料, 以正射影像图为基础图像进行岩性识别和构造解译, 着重解译第四系、断层、微地貌及阶地等的特征、位置和范围。境外地区的地质解译则是在掌握我国境内地质情况及遥感影像特征基础上, 利用对比分析方法, 按照从已知到未知的原则, 采取推演、扩展、外延的方法进行。

3 地质背景遥感解译
3.1 地层解译

应用GF-1遥感数据, 结合前人1∶ 25万区域地质调查成果资料, 对各类地质体进行了较详细的遥感解译, 尤其是对不同成因的第四系及其空间展布进行了识别与圈定。基岩区可见友谊桥岩组(Pt1y)变粒岩、石英片岩和石英岩, 达巴劳组(O1d)变质长石石英砂岩, 石器坡组(S1s)粉砂质板岩夹薄层变质石英砂岩, 普鲁组(S2-3p)厚层状细晶灰岩, 凉泉组(D1l)白云岩及白云质板岩, 色龙群(P2-3S)石英砂岩夹长石砂岩, 聂聂雄拉组下段(J2n1)泥晶灰岩夹粉砂岩、石英砂岩和聂聂雄拉组上段(J2n2)亮晶颗粒灰岩等岩性出露。上述地层在GF-1图像中多以“ V” 形沟谷、均一色调为主要特征, 灰岩出露区表现为典型的尖棱状“ 灰岩陡崖” 微地貌景观; 片岩在图像中以发育与片理走向基本一致的平行状纹形、色带为典型特征; 变质砂岩、石英岩则呈似层状影像特征, 常形成较低缓的光秃圆滑垅脊或丘岗。

下更新统香孜组(Qp1x)地层仅残存于现今朗钦藏布两岸, 主要沿楚鲁松杰村— 昌格里赞格— 纳姆吉亚一带展布, 出露的海拔高度为3 300~4 200 m; 在GF-1图像中呈灰褐色或浅黄褐色色调, 为较平坦的层状地貌, 影像结构细腻光滑, 局部地段被后期冲洪积所改造, 发育明显的扇状、树枝状水系; 具有较明显的河流相沉积特征, 其中昌格里赞格地区河流相沉积保留比较完整, ⑦⑧⑨级河流阶地构成了区内的最高级阶地, 与邻区札达盆地附近的沉积对比[20, 21], 其时代应为早更新世。④⑤⑥级河流阶地为中更新世冲积, ②③级河流阶地为晚更新世冲积(图2)。

图2 昌格里赞格地区河流阶地Fig.2 River terrace in Changgelizange rigion

GF-1图像显示, 古河流一直向西(现今朗钦藏布方向)侧蚀迁移, 且河流阶地的阶面逐渐变窄, 前缘阶坡垂直高度增大, 反映该区经历了强烈的差异性隆升与河流下切作用, 也表现出青藏高原间歇性抬升愈演愈烈的频繁性和剧烈性。

中更新统冲积(Q pal2)沿河岸山坡零星分布, 表现为山前坡麓、山前倾斜台地, 出露海拔高度多在3 900~4 000 m。图像结构光滑, 形成平坦的高位阶地, 边部冲沟密布, 水系不发育(图3)。该套地层与全新统河流沉积相伴, 对应于狮泉河镇西部的中更新统河流相沉积。该时期为青藏高原40~30 ka的“ 泛湖期” [22], 邻区的班公错古湖达到最高湖面[23]。中更新统冰碛(Q pgl2)构成高侧碛垄, 在地貌上呈突起的冰碛小丘或长条状冰碛垅。

图3 中更新统河流相沉积Fig.3 River deposition of middle Pleistocene

上更新统冰水堆积(Q pfgl3)分布于沟谷的山前坡脚地带, 呈平缓扇状地形。上更新统冲积(Q pal3)出露的海拔高度为3 200~3 600 m, 沿河岸两侧断续分布, 表现为宽窄不一的阶地。上更新统冲洪积(Q ppal3)沿山麓或河谷呈带状、扇状分布。

全新统冰碛(Qhgl)位于现代冰川冰舌前沿及冰川退缩不久的冰斗之中, 呈长条形、弧形, 舌形前缘清晰可见, “ U” 形冰槽谷明显, 具终碛堤或终碛垅地貌特征。全新统坡积(Qhdl)分布于喜马拉雅山山脉南北陡峻山地的岩石坡面上, 是基岩经长期寒冻作用发生崩解的岩石沿坡面堆积的产物。全新统冲洪积(Qhpal)以面状展布于现代河谷宽谷地带。全新统洪积(Qhpl)主要分布于山间沟谷出口处, 呈山间宽谷及山前扇形地貌。全新统冲积(Qhal)沿现代河流呈带状、条带状或条状分布, 影像结构细腻光滑。

3.2 构造解译

区内断裂和褶皱构造发育, 其中以NW向断裂为主, NE(或NNE)向和近SN向断裂多为后期断裂, 常将NW向断裂错断(图4)。

图4 断裂构造影像特征Fig.4 Remote sensing image feature of fault structure

3.3 岩浆岩解译

区内岩浆活动较强烈, 以中酸性岩浆岩为主, 主要分布于研究区东北部, 岩性为中细粒二长花岗岩(η γ E3), 在遥感图像中呈浅灰白色— 浅乳黄色色调, 形态多为不规则块状, 典型树枝状水系发育, 普遍具较为细密的斑点状影纹, 影纹结构较为粗糙。

4 地质灾害发育特征遥感分析
4.1 灾害体解译效果

崩塌在GF-1图像中边界完整, 清晰可辨; 陡坎状后壁连续齐整, 色调较周围地物浅, 阴影明显; 崩塌堆积区呈灰白色或灰褐色, 其线理结构和图案与周围明显不同, 呈团块状或流线状纹理; 崩塌体在坡脚形成锥状堆积体或倒石锥, 边界呈弧形锯齿状, 崩积物多而不均, 故影像粗糙, 常显斑状纹理。因此, 在GF-1图像中能够识别出崩塌区和崩塌堆积区的空间形态和展布特征, 但对发育在后壁的次级崩塌或小型崩塌体难以识别(图5)。

图5 崩塌影像特征Fig.5 Remote sensing image feature of collapse

滑坡在GF-1多光谱图像中一般显示簸箕形、舌形、弧形和环形等地貌特征, 圈状弧形突出, 后壁较陡, 坡度可达[30° , 40° ], 后壁之上冲沟不发育。滑坡体多呈台阶状, 故形成一系列近于平行展布的陡坎, 滑坡体的整体纹理和色调与背景环境明显不一致[10]。本文利用GF-1融合图像解译出滑坡的空间展布形态, 滑坡堆积体的长度、宽度等特征部位信息(图6), 也能够识别出古滑坡、新滑坡等类型。

图6 滑坡影像特征Fig.6 Remote sensing image feature of landslide

区内泥石流在GF-1多光谱图像中的轮廓比较清晰, 物源区、流通区和堆积区容易辨别。物源区山体较破碎, 呈勺状、漏斗状、椭圆状等三面环山之冰蚀围谷, 第四系松散堆积物(冰碛、残坡积)较丰富, 坡面及细沟侵蚀严重, 支沟普遍发育; 流通区内水系呈线状或树枝状, 沟谷曲直不一、切割明显, 沟谷一般较窄、断面形态呈“ V” 形; 堆积区多为沟谷下游出口, 地形突变而趋于平缓, 常形成扇形堆积, 扇形大小不一。堆积扇上水流不固定, 多呈漫流, 影像结构粗细间杂。因GF-1多光谱融合图像中的色调不够丰富, 局部堆积区的色调、纹理特征不甚明显, 对植被的反映也不够清晰, 泥石流的发育阶段不易判识(图7)。

图7 泥石流影像特征Fig.7 Remote sensing image feature of debris flow

寒冻风化碎屑流主要呈散体结构分布于我国西部高寒山区, 成分以碎石及少量岩屑为主, 外貌特征与大型采石场加工碎石料后的堆石特征类似, 为典型的由寒冻风化引起的基岩剥蚀产物, 受高山冰缘环境中的海拔、坡向、温度和降水、岩性构造及地形条件等条件制约[24, 25, 26]。区内寒冻风化碎屑流多分布于陡坡之上, 呈面状展布, 影像粗糙, 灰白或灰褐色调, 具斑点状或斑状纹理, 在GF-1多光谱图像中可以判识其位置和范围等特征信息(图8)。

图8 寒冻风化碎屑流影像特征Fig.8 Remote sensing image feature of frost weathering debris flow

4.2 灾害体分布规律

经详细综合解译, 发现区内共发育地质灾害86处, 主要类型为滑坡和寒冻风化碎屑流, 崩塌和泥石流次之。其中: 崩塌5处, 滑坡41处, 泥石流4条, 寒冻风化碎屑流36处。参照中国地质调查局地质调查技术标准《滑坡崩塌泥石流灾害详细调查规范( 1∶ 50 000)》中的分类原则, 区内发生的地质灾害按规模划分为特大型1处、大型9处、中型39处、小型37处。

4.2.1 崩塌

研究区内的崩塌主要发生于朗钦藏布(印度称之为萨特莱杰河)两岸, 均位于高山峡谷地段, 共5处, 海拔高程范围为3 400~3 600 m, 均发育在友谊桥岩组(Pt1y)石英片岩和石英片岩夹变粒岩组合中。

利用DEM数据将全区坡度分成4个级别: [5° , 15° )为平缓, [15° , 25° )为缓, [25° , 35° )为中陡, ≥ 35° 为高陡。统计结果表明, 崩塌主要发育在中陡和高陡的坡度段内, 缓坡度地段发育崩塌较少(仅占20%), 均为中、小型崩塌(表2)。

表2 崩塌基本特征统计 Tab.2 Statistical basic characters of collapses

4.2.2 滑坡

区内滑坡主要发生在第四系松散堆积物和海拔3 200~4 800 m并以沉积岩为主的冰蚀发育区, 共41处。其中发育在第四系松散堆积物(主要为下更统世香孜组)中的有27处, 约占滑坡总数的65.85%; 其余主要发生在色龙群(P2-3S)石英砂岩夹粉砂岩组合和友谊桥岩组(Pt1y)石英片岩夹变粒岩组合中。

统计结果表明, 滑坡主要发生在[5° , 15° )坡度段内, 约占滑坡总数的53.66%; 发生在[15° , 25° )坡度段内的滑坡占39.02%; ≥ 25° 的地段发育滑坡较少, 所占比例仅为7.32%; 区内滑坡均为中、小型滑坡(表3)。

表3 滑坡基本特征统计 Tab.3 Statistical basic characters of landslides

4.2.3 泥石流

区内泥石流共发育4条, 以沟口堆积型为主, 规模均较小, 所处高程多在2 900~5 900 m之间, 流域坡度均在[16° , 27° ]坡度段内。研究发现, 区内波石流发育受控于友谊桥岩组(Pt1y)石英片岩、石英片岩夹变粒岩, 色龙群(P2-3S)石英砂岩夹粉砂岩以及第四系洪积、冰碛、冰水堆积、坡积等沉积物。局部泥石流沟两侧可见滑坡和崩塌, 细沟和坡面侵蚀严重(表4)。

表4 泥石流基本特征统计 Tab.4 Statistical basic characters of mudslides

4.2.4 寒冻风化碎屑流

区内寒冻风化碎屑流多发生在海拔高程3 400~5 800 m之间, 其物质组成主要为全新统坡积(Qhdl)砂土、碎石以及少量岩屑, 原岩为聂聂雄拉组二段(J2n2)灰岩夹砂岩和板岩, 凉泉组(D1l)白云岩及白云质板岩和友谊桥岩组(Pt1y)石英片岩及石英片岩夹变粒岩等。碎屑流相对高度多在200~300 m之间, 最高可达1 600 m, 主要为中型(占发育总数的72.22%)。大部分碎屑流发生在> 35° 坡度段内, 所占比例约为66.67%; 其次发生在[25° , 35° ]坡度段内, 约占30.56%(表5)。

表5 寒冻风化碎屑流基本特征统计表 Tab.5 Statistical basic characters of frost weathering debris flows
5 地质灾害形成条件分析
5.1 地形地貌条件

研究区内高山峡谷发育, 地形高差悬殊, 说明曾经发生过强烈的外动力地质作用, 为地质灾害的发生提供了足够的势能和巨大的空间。不同坡度段内发育的地质灾害统计见图9

图9 不同坡度段内发育的地质灾害统计Fig.9 Statistical diagram of geological disasters in different slope sections

对解译结果的研究发现, 滑坡、崩塌和泥石流在不同坡度段内具较好的相关性, 多集中发生在中陡的[25° , 35° )坡度段内; 而寒冻风化碎屑流的发生则与坡度的增大成正相关— — 在一定的区间内, 坡度越大, 寒冻风化碎屑流越易发生, 表明区内地质灾害的发生明显受地形地貌条件的控制。

5.2 地质构造条件

地质灾害的形成和分布主要受地形地貌、地层岩性、地质构造、新构造活动、水文地质条件、地震以及气候、人类活动等因素制约[27, 28]。其中, 地层岩性、地质构造和新构造活动是地质灾害形成的基本条件。

5.2.1 地层岩性

区内发生地质灾害共计86处(图10)。按不同地层中地质灾害发育的个数来统计, 与第四系松散堆积物有关的地质灾害有67处(占地质灾害总数的77.91%), 说明中陡和高陡坡度段内发育的大量松散堆积物为地质灾害的形成奠定了丰富的物质基础。从岩性及水文地质条件来看, 松散堆积物孔隙度大, 透水性强, 而下伏基岩透水性相对较差, 易形成隔水层。地表降水的入渗使松散堆积物呈饱和状态, 下渗后的地下水软化土体, 增大土体重量和孔隙水压力, 使土体力学强度降低, 因而易在土体中或基岩面上形成滑带或滑移面, 导致松散堆积物滑动, 形成滑坡。

图10 研究区地质简图Fig.10 Simplified geological map of study area

基岩中发育的地质灾害有19处(占地质灾害总数的22.09%), 主要为滑坡和崩塌。上述地质灾害主要发生在友谊桥岩组(Pt1y)石英片岩和石英片岩夹变粒岩组合及聂聂雄拉组下段(J2n1)的灰岩夹砂岩组合中, 这些地层中的岩性软— 硬、厚— 薄相间, 多具软弱夹层, 且岩层间多以规模不等的片理化带为界, 各岩性层片理发育, 在寒冻风化、降水、雪融作用以及上覆重物的作用下极易发育成滑动面, 从而发生地质灾害[29]

5.2.2 地质构造

研究区由于青藏高原的强烈隆升, 其周缘地区构造变形十分强烈, 加之气候环境特殊, 成为地质灾害事件集中发生的区域[30]。地质构造对区内地质灾害的形成和发育具有明显的影响: ①控制了地貌的形成和发育。区内的地貌特征反映了地壳上升、河流下切的运动过程, 以本区河流相沉积为例, 多级河流阶地所组成的层状地貌就是新构造运动最近时期的产物, 反映了地壳的间歇性垂向运动特征, 是高原地面抬升的直接证据[31] 。区内广泛分布多级河流阶地, 且早更新世河流相沉积与现代河流的海拔高程相差近1 000 m, 标志着现代地壳构造运动仍存在显著的继承性垂直运动, 也表明了0.15 Ma (共和运动[32])以来该区又进入了强烈的隆升时期。差异性隆升造就了区内山势险峻、沟谷深切、临空面发育的地貌特征(如在朗钦藏布两岸形成300~1 500 m的高陡边坡), 为区内地质灾害的发育提供了临空条件。②限制了地质灾害的分布范围。由图10可以看出, 区内主构造线方向为NW向, NE向和近SN向次之。地质灾害的条带状分布基本上受控于这3组断裂带, 大部分地质灾害发育在断裂带附近。同时, 对区内基础地质资料分析和地质灾害运移方向的统计(图11)表明: 基岩层面倾向为110° ~160° , 地质灾害的运移方向主要为220° 和150° , 分别与主构造线方向、基岩层面倾向基本一致。这一事实说明, 区内断裂、大型节理、面理和层理等成灾构造是导致在这一地区频发崩塌和滑坡等地质灾害的直接原因。③改变了岩土体的结构、物理性质和力学强度。以寒冻风化碎屑流为例, 影响寒冻风化速率的主要因素是节理的密集程度, 特别是有断裂带通过的地段会加速寒冻风化[19]。而区内构造变形十分强烈, 岩石中发育大量节理和面理构造, 且在压力、剪切力、重力等作用下, 岩石易发生破碎和位移, 也就更易发生寒冻风化碎屑流。

图11 地质灾害运移方向玫瑰花图Fig.11 Rose pattern for migration direction of geological disasters

6 结论

1)高分一号卫星(GF-1)数据的融合图像对地质体(或灾害体)的微地貌形态特征反映较明显, 能够用于识别灾害体的位置、边界范围等特征信息。

2)本文依据微地貌特征, 在GF-1融合图像中识别并圈定出不同类型的第四系。在楚鲁松杰村— 昌格里赞格— 纳姆吉亚一带新解译出下更新统香孜组地层、中更新统冲积和上更新统冲积, 完善了区内第四系地层系统, 为研究区内新构造运动及青藏高原隆升补充了与第四系有关的资料。

3)区内共发育地质灾害86处, 主要灾害类型为滑坡和寒冻风化碎屑流, 崩塌和泥石流次之。其中: 滑坡41处, 崩塌5处, 泥石流4条, 寒冻风化碎屑流36处。

4)滑坡、崩塌和泥石流在不同坡度段内具较好的相关性, 多集中发生在中陡的[25° , 35° ]坡度段内; 而寒冻风化碎屑流的发生则与坡度的增大成正相关, 在一定的区间内, 坡度越大, 寒冻风化碎屑流越易发生。

5)地质灾害的运移方向分别与主构造线方向和基岩层面倾向基本一致。

6)区内第四系松散堆积物广泛分布(大多数分布在海拔高程4 000 m左右), 较大的垂直落差及组成物质的松散性成为诱发滑坡、泥石流和寒冻风化碎屑流的重要因素。

The authors have declared that no competing interests exist.

参考文献
[1] 李加洪, 蒋卫国, 张松梅, . 基于遥感与GIS的西藏帕里河滑坡动态监测分析[J]. 自然灾害学报, 2006, 15(4): 24-27.
Li J H, Jiang W G, Zhang S M, et al. RS and GIS based land slide monitoring of Pali River in Tibet[J]. Journal of Natural Disasters, 2006, 15(4): 24-27. [本文引用:2]
[2] 王治华, 徐起德, 杨日红, . 中印边界附近帕里河上的滑坡灾害遥感调查[J]. 科技导报, 2007, 25(6): 27-31.
Wang Z H, Xu Q D, Yang R H, et al. Remote sensing survey for a land slide disaster in the boundary river between China and India[J]. Science & Technology Review, 2007, 25(6): 27-31. [本文引用:2]
[3] 王治华. 大型个体滑坡遥感调查[J]. 地学前缘, 2006, 13(5): 516-523.
Wang Z H. Large scale individual land slide remote sensing[J]. Earth Science Frontiers, 2006, 13(5): 516-523. [本文引用:1]
[4] 王治华. 中国滑坡遥感及新进展[J]. 国土资源遥感, 2007, 19(4): 7-1032. doi: DOI: 106046/gtzyyg. 2007. 04. 02.
Wang Z H. Remote sensing for land slides in China and its recent progress[J]. Remote Sensing for Land and Resources, 2007, 19(4): 7-1032. doi: DOI:10.6046/gtzyyg.2007.04.02. [本文引用:1]
[5] 杨金中, 杨日红. 遥感技术在三峡库区千将坪滑坡研究中的应用[J]. 国土资源遥感, 2007, 19(4): 85-89. doi: DOI: 106046/gtzyyg. 2007. 04. 19.
Yang J Z, Yang R H. The application of remote sensing digital technique to the investigation of the Qianjiangping land slide in the Three Gorges Area[J]. Remote Sensing for Land and Resources, 2007, 19(4): 85-89. doi: DOI:10.6046/gtzyyg.2007.04.19. [本文引用:1]
[6] 童立强, 祁生文, 刘春玲. 喜马拉雅山东南地区地质灾害发育规律初步研究[J]. 工程地质学报, 2007, 15(6): 721-729.
Tong L Q, Qi S W, Liu C L. Preliminary study of geohazard development patterns in the southeast region of Himalaya mountains[J]. Journal of Engineering Geology, 2007, 15(6): 721-729. [本文引用:1]
[7] 秦绪文, 张志, 杨军杰, . 5. 12汶川地震四川省平武县次生地质灾害遥感分析[J]. 地质科技情报, 2009, 28(2): 12-15.
Qin X W, Zhang Z, Yang J J, et al. Analysis on the secondary geological disaster in Pingwu of5·12 Wenchuan Earthquake, Sichuan based on remote sensing information[J]. Geological Science and Technology Information, 2009, 28(2): 12-15. [本文引用:1]
[8] 王润生, 熊盛青, 聂洪峰, . 遥感地质勘查技术与应用研究[J]. 地质学报, 2011, 85(11): 1699-1743.
Wang R S, Xiong S Q, Nie H F, et al. Remote sensing technology and its application in geological exploration[J]. Acta Geologica Sinica, 2011, 85(11): 1699-1743. [本文引用:1]
[9] 朱静, 唐川. 遥感技术在我国滑坡研究中的应用综述[J]. 遥感技术与应用, 2012, 27(3): 458-464.
Zhu J, Tang C. An overview of remote sensing applications for land slide research in China[J]. Remote Sensing Technology and Application, 2012, 27(3): 458-464. [本文引用:1]
[10] 童立强, 郭兆成. 典型滑坡遥感影像特征研究[J]. 国土资源遥感, 2013, 25(1): 86-92. doi: DOI: 106046/gtzyyg. 2013. 01. 16.
Tong L Q, Guo Z C. A study of remote sensing image features of typical land slides[J]. Remote Sensing for Land and Resources, 2013, 25(1): 86-92. doi: DOI:10.6046/gtzyyg.2013.01.16. [本文引用:2]
[11] 童立强, 聂洪峰, 李建存, . 喜马拉雅山地区大型泥石流遥感调查与发育特征研究[J]. 国土资源遥感, 2013, 25(4): 104-112. doi: DOI: 106046/gtzyyg. 2013. 04. 17.
Tong L Q, Nie H F, Li J C, et al. Survey of large-scale debris flow and study of its development characteristics using remote sensing technology in the Himalayas[J]. Remote Sensing for Land and Resources, 2013, 25(4): 104-112. doi: DOI:10.6046/gtzyyg.2013.04.17. [本文引用:1]
[12] 魏永明, 魏显虎, 陈玉. 岷江流域映秀——茂县段地震次生地质灾害分布规律及发展趋势分析[J]. 国土资源遥感, 2014, 26(4): 179-186. doi: DOI: 106046/gtzyyg. 2014. 04. 28.
Wei Y M, Wei X H, Chen Y. Analysis of distribution regularity and development tendency of earthquake secondary geohazards in Yingxiu-Maoxian section along the Minjiang River[J]. Remote Sensing for Land and Resources, 2014, 26(4): 179-186. doi: 0.6046/gtzyyg.2014.04.28. [本文引用:1]
[13] 张祖勋, 张永军. 利用国产卫星影像构建我国地理空间信息[J]. 测绘地理信息, 2012, 37(5): 7-9.
Zhang Z X, Zhang Y J. Building geospatial information of China with domestic satellite imagery[J]. Journal of Geomatics, 2012, 37(5): 7-9. [本文引用:1]
[14] 路云阁, 刘采, 王姣. 基于国产卫星数据的矿山遥感监测一体化解决方案——以西藏自治区为例[J]. 国土资源遥感, 2014, 26(4): 85-90. doi: DOI: 106046/gtzyyg. 2014. 04. 14.
Lu Y G, Liu C, Wang J. Integrated solutions for mine remote sensing monitoring based on domestic satellite images: A case study of Tibet[J]. Remote Sensing for Land and Resources, 2014, 26(4): 85-90. doi: DOI:10.6046/gtzyyg.2014.04.14. [本文引用:1]
[15] 陈玲, 梁树能, 周艳, . 国产高分卫星数据在高海拔地区地质调查中的应用潜力分析[J]. 国土资源遥感, 2015, 27(1): 140-145. doi: DOI: 106046/gtzyyg. 2015. 01. 22.
Chen L, Liang S N, Zhou Y, et al. Potential of applying domestic high-resolution remote sensing data to geological survey in high altitudes[J]. Remote Sensing for Land and Resources, 2015, 27(1): 140-145. doi: DOI:10.6046/gtzyyg.2015.01.22. [本文引用:1]
[16] 朱利, 李云梅, 赵少华, . 基于GF-1号卫星WFV数据的太湖水质遥感监测[J]. 国土资源遥感, 2015, 27(1): 113-120. doi: DOI: 106046/gtzyyg. 2015. 01. 18.
Zhu L, Li Y M, Zhao S H, et al. Remote sensing monitoring of Taihu Lake water quality by using GF-1 satellite WFV data[J]. Remote Sensing for Land and Resources, 2015, 27(1): 113-120. doi: DOI:10.6046/gtzyyg.2015.01.18. [本文引用:1]
[17] 西藏地质调查院. 1∶25万斯诺乌山幅、狮泉河幅地质调查成果与进展[J]. 沉积与特提斯地质, 2005, 25(1/2): 57-61.
Xizang Institute of Geological Survey. 1∶250 000 Sinuowushan and Shiquanhe Sheets in Xizang[J]. Sedimentary Geology and Tethyan Geology, 2005, 25(1/2): 57-61. [本文引用:1]
[18] 河北省地质调查院区调所. 1∶25万日新幅、札达县幅、姜叶马幅地质调查成果与进展[J]. 沉积与特提斯地质, 2005, 25(1/2): 71-79.
Hebei Institute of Geological Survey. 1∶250 000 Rixin, Zand a and Gya'nyima Sheets in Xizang[J]. Sedimentary Geology and Tethyan Geology, 2005, 25(1/2): 71-79. [本文引用:1]
[19] 邱学雷. 高分一号卫星工程首批影像图发布将为国土、环境、农业等领域提供精准服务[J]. 国防科技工业, 2013(6): 14-16.
Qiu X L. Gaofen-1 satellite first released image map which will provide accurate service to land resources, environment and agriculture[J]. Defence Science and Technology Industry, 2013(6): 14-16. [本文引用:2]
[20] 朱大岗, 孟宪刚, 邵兆刚, . 西藏阿里札达盆地上新世-早更新世的古植被、古环境与古气候演化[J]. 地质学报, 2007, 81(3): 295-306.
Zhu D G, Meng X G, Shao Z G, et al. Evolution of the paleovegetation, paleoenvironment and paleoclimate during pliocene-early pleistocene in Zhada Basin, Ali, Tibet[J]. Acta Geologica Sinica, 2007, 81(3): 295-306. [本文引用:1]
[21] 朱大岗, 孟宪刚, 邵兆刚, . 西藏札达盆地上新统-早更新统河湖相地层多重地层划分与对比[J]. 地球学报, 2008, 29(2): 137-144.
Zhu D G, Meng X G, Shao Z G, et al. Division and correlation of pliocene and early pleistocene fluvio-lacustrine strata in the Zhada Basin, Tibet, China[J]. Acta Geoscientica Sinica, 2008, 29(2): 137-144. [本文引用:1]
[22] 郑绵平, 袁鹤然, 赵希涛, . 青藏高原第四纪泛湖期与古气候[J]. 地质学报, 2006, 80(2): 169-180.
Zheng M P, Yuan H R, Zhao X T, et al. The Quaternary Pan-lake(Overflow)Period and Paleoclimate on the Qinghai-Tibet Plateau[J]. Acta Geologica Sinica, 2006, 80(2): 169-180. [本文引用:1]
[23] 李炳元. 青藏高原大湖期[J]. 地理学报, 2000, 55(2): 174-182.
Li B Y. The Last greatest lakes on the Xizang(Tibetan)plateau[J]. Acta Geographica Sinica, 2000, 55(2): 174-182. [本文引用:1]
[24] 李树德. 天山乌鲁木齐河源胜利达坂岩石风化剥蚀速率初探[J]. 冰川冻土, 1981, 3(s1): 114-119.
Li S D. Study on denudation rate of Heyuanshenglidaban rock of Tianshan mountains in Urumchi City[J]. Journal of Glaciology and Geocryology, 1981, 3(s1): 114-119. [本文引用:1]
[25] 刘耕年, 熊黑钢. 中国天山高山冰缘环境中的寒冻风化剥蚀作用及其影响因素[J]. 冰川冻土, 1992, 14(4): 332-341.
Liu G N, Xiong H G. The denudation rate of gelifraction and its influential factors in the periglacial environment of Tianshan mountains(in China)[J]. Journal of Glaciology and Geocryology, 1992, 14(4): 332-341. [本文引用:1]
[26] 王启国, 张良平. 西藏中波河段大型碎屑流成因机制及防治对策[J]. 岩石力学与工程学报, 2011, 30(11): 2293-2302.
Wang Q G, Zhang L P. Formation mechanism and preventive countermeasure of large-scale frost weathering debris flow on Zhongbo reach of Yuqu river in Tibet[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(11): 2293-2302. [本文引用:1]
[27] 段永侯. 我国地质灾害的基本特征与发展趋势[J]. 第四纪研究, 1999(3): 208-216.
Duan Y H. Basic characters of geo-hazards and its development trend in China[J]. Quaternary Sciences, 1999(3): 208-216. [本文引用:1]
[28] 张春山, 张业成, 胡景江, . 中国地质灾害时空分布特征与形成条件[J]. 第四纪研究, 2000, 20(6): 559-564.
Zhang C S, Zhang Y C, Hu J J, et al. Spatial and temporal distribution characteristics and forming conditions of Chinese geological disasters[J]. Quaternary Sciences, 2000, 20(6): 559-564. [本文引用:1]
[29] 刘春玲, 祁生文, 童立强, . 喜马拉雅山地区重大滑坡灾害及其与地层岩性的关系研究[J]. 工程地质学报, 2010, 18(5): 669-676.
Liu C L, Qi S W, Tong L Q, et al. Great land slides in Himalaya mountain area and their occurrence with lithology[J]. Journal of Engineering Geology, 2010, 18(5): 669-676. [本文引用:1]
[30] 彭建兵, 马润勇, 卢全中, . 青藏高原隆升的地质灾害效应[J]. 地球科学进展, 2004, 19(3): 457-466.
Peng J B, Ma R Y, Lu Q Z, et al. Geological hazards effects of uplift of Qinghai-Tibet plateau[J]. Advance in Earth Sciences, 2004, 19(3): 457-466. [本文引用:1]
[31] 潘保田, 李吉均, 李炳元. 青藏高原地面抬升证据讨论[J]. 兰州大学学报: 自然科学版, 2000, 36(3): 100-111.
Pan B T, Li J J, Li B Y. Discussion on evidences for surface uplift of the Qinghai Tibetan Plateau[J]. Journal of Lanzhou University: Natural Sciences, 2000, 36(3): 100-111. [本文引用:1]
[32] 李吉均, 方小敏, 马海洲, . 晚新生代黄河上游地貌演化与青藏高原隆起[J]. 中国科学(D辑), 1996, 26(4): 316-322.
Li J J, Fang X M, Ma H Z, et al. Geomorphological and environmental evolution in the upper reaches of the Yellow River during the Late Cenozoic[J]. Science in China(Series D), 1996, 26(4): 316-322. [本文引用:1]