加拉白垒峰

第腮卷第2期

2021年4月

甘肃科学学报

JournalofGansu

Sciences

Vol.33No.2

Apr.2021

引用格式:DongJihong,YangChcngshcng,lGlacierDisplacementMonitoringofGy-

alaPeriBadonSAROfftT'racking'Technology[J].JournalofGansu

Sciences

,021,33(2)：1-7.3董继红，杨成

生，

张本浩，等.基于SAR偏移量跟踪技术的加拉白垒峰典型冰川位移监测JJ.甘肃科学学报,021,33(2):17」

doi:10.16/168/004-0366.2021.02.001.

基于SAR偏移量跟踪技术的加拉白垒峰

典型冰川位移监测

董继红」，杨成生」

，张本浩」

，赵超英12

(.长安大学地质工程与测绘学院，陕西西安710054;

2.地理信息工程国家重点实验室，陕西西安710054)

摘要随着淡水资源的日益紧缺，冰川作为重要的淡水资源，与人民生活的关系愈发重要，同时冰川

消融引起的地质灾害也深深影响着人们的生命和财产安全。以西藏加拉白垒峰区域典型冰川为研

究对象，选用日本的L波段ALOS-2合成孔径雷达影像数据，采用基于强度的偏移量追踪技术，对

加拉白垒峰周边冰川位移变化开展了监测研究。结果表明，加拉白垒峰北侧岗普冰川630天内方位

向最大累积形变量超过60m,其位移累积量和流速受其东侧支流的影响较大。相关研究结果可为

藏南地区类似冰崩碎屑流灾害的减灾和防灾工作提供参考。

关键词合成孔径雷达；加拉白垒峰；冰川；偏移量跟踪技术

中图分类号:P237文献标志码:A文章编号：1004-0366(2021)02-000107

在全球气候变暖背景下，冰川的加速消融和退

缩对气候的直接和间接影响也变得更加突出［

1］。

冰

川是气候变化的天然指示器，对冰川动态变化的监

测，一直备受学术界和政府关注20世纪80年

代以来全球快速增温，青藏高原气候也发生着显著

的变化，高原冰川随气候变化发生了一系列的退

缩［3］。由于冰川消融将会导致海平面上升，引发区

域淹没、冰雪灾害、冰川洪水等自然灾害

，这些都威

胁着人类居住环境的安全，并可能改变许多动植物

种群的习性和生存环境［4］，因此,对冰川变化的监

测，一安全生产自查报告 方面有助于了解高原气候的变化状况，另一方

面对保障人类生存环境的安全具有重要意义⑴.

因冰川多分布于高海拔、高陡、高寒地区，

人工

现场勘察面临工作强度大、成本高、人员无法到达等

问题。基于卫星遥感技术对冰川实施动态监测，具

有大范围、无需现场测量等优势，成为揭示冰川的变

化规律和退化状况的重要工具。然而，光学遥感技

术对冰川的动态监测主要以遥感影心律不齐什么原因 像解译与定性分

析为主归，较难对冰川的运动进行准确定量分析。

合成孔径雷达(SAR,syntheticapertureradar)因其

具有高时空分辨率、大范围监测、

不受云雾影响等技

术优势⑺，可以弥补光学遥感监测手段的不足。目

前,SAR及干涉合成孔径雷达(InSAR,

interfero

metricSAR)遥感技术已在与地表覆盖变化及地表

形变有关的许多地球物理现象的监测与解释中得到

了大量的应用［8-9］,其中偏移量跟踪技术在监测大量

级形变中有无可比拟的优势，尤其是在冰川灾害监

测中,1982年Gray等［0］将图像配准方法用于冰川

在距离和方位向的运动估计。2002年Strozzi等［11

］

利用优化的强度追踪技术对斯瓦尔巴特群岛北部的

莫纳科布林冰川进行监测，获得了有效的结果。在

国内刘国祥等［12］基于像素偏移量技术与干涉堆叠

方法，对海螺沟冰川利用地基和星基SAR数据探测

了消融区内冰川退化所引起的次生滑坡灾害

，随着

收稿日期：2020-03-18；修回日期：2020-04-23

基金项目：国家重点研发计划(2018YFC1504805)；中国地质灾害调查项目(DD20190637)；中央高校专项基金项目(CHD,)；

陕西

省自然科学基金(2019JM-245)

作者简介：董继红(1994-),男，甘肃天水人，硕士研究生，研究方向为

InSAR技术在地质灾害的应用.E-mail：*****************

通信作者：杨成生,E-mail：**************

2

甘肃科学学报

2021年第2期

高精度SAR数据的可获取程度越来越高，

偏移量追

踪技术应用愈加广泛。

研究区域加拉白垒峰位于雅鲁藏布江大拐弯

外侧，海拔7294m,是世界第85高峰。该地区冰

川属于海洋型冰川，运动较快，受气温和降水的影

响，加之地势陡峭，冰、雪崩十分频繁。2018年10

月17日凌晨，位于加拉白垒峰西侧的色东普沟冰

川发生了大规模冰川泥石流灾害，堵断了雅鲁藏布

江干流并形成堰塞湖。同年10月29日，色东普沟

因冰雪融水引发的高浓度泥石流再次堵塞雅鲁藏

布江。据调查这两次事件均与冰崩密切相关［灯，

两次堵江灾害均对雅鲁藏布江上下游一定范围内

的村镇造成了严重威胁。为此研究利用星载SAR

数据和GoogleEarth™光学数据，采用SAR影像

偏移量追踪(OfftTracking)技术对加拉白垒峰

周边典型冰川变化及分布特征进行了监测，并分析

了区域内冰川的活动特征，为今后类似灾害的监测

与防治提供参考。

1研究区概况

加拉白垒峰属于念青唐古拉山脉，走向为东西

弧形排列，主脊线上向南侧、西北侧增生着数条支

脊，多为陡壁悬崖。山谷中发育着数十条冰川，其顶

部比较平展，常年被冰雪覆盖。受气温、降水及地势

的影响，区域内冰、雪崩时有发生。色东普冰川、

岗

普冰川及则巴弄巴曲冰川为加拉白垒峰周边规模较

大的冰川，分别分布在加拉白垒峰的南、北和西侧

(见图1).

(a)

研究区地理位置分布情况

(b)色东普冰川光学影像(c)岗普冰川(左)和则巴弄巴曲冰川(右)光学影像

图1研究区域地理位置及冰川地貌

Fig.1

Geographical

locationandglaciallandformo歌颂党的征文 fthestudyarea

色东普冰川属于雅鲁藏布江大峡谷段的色东普

域，位于西藏米林市米林县派镇加拉村下游6km

处雅鲁藏布江大拐弯的左岸(见图1).沟域上游源

区地形宽阔、支沟发育，中下游主沟道狭窄。沟道内

冰川活动形成的冰磧物丰富，冰雪融水及降水提供

水流充分。沟域上游陡峭，地形纵坡降大，冰川发

育，岩土体物理风化严重，侵蚀剥蚀作用强烈"14］。

岗普冰川与色东普冰川的方向相反.区域内最高海

第33卷

董继红等：基于SAR偏移量跟踪技术的加拉白垒峰典型冰川位移监测

3

拔4169m,最低点为色东普沟口3438m,其东侧

分布着一条主要的支流，作为冰川水资源的补给.

2研究数据及方法

2.1数据介绍

为了监测加拉白垒峰周边冰川运动变化情况,

我们选取了星载SAR和光学GoogleEarth™数据

开展实验，其中SAR数据为2014年5月日本空间

局发射的ALOS-2卫星搭载L波段传感器所采集

的数据。实验收集了覆盖研究区域5景升轨SM3

模式数据，时间跨度从2017年2月一2018年10

月，研究所采用的ALOS-2SAR影像数据的详细信

息见表1考虑到该区域地形起伏较大，我们引入

30m分辨率的数字高程模型(SRTM,shuttleradar

topographymission)(https：//.

gov/)1来辅助SAR影像配准，同时用于纠正SAR

影像的局部入射角偏差引起的方位向与距离向偏移。

表1ALOS-2数据参数

Table1ALOS-2dataparameter

传感器

ALOS/PALSAR-2

轨道方向

升轨

航向角/()

—10.7

入射角/()31.4

波长

L

方位向分辨率/mX

距离向分辨率/m

3.2X4.3

影像数量

5

时间覆盖

2017年2月2日—2018年10月25日

p（x,

丫工了(厂(x,y)—u

（工x了（厂（x,y）—ur

其中：p(x,y)为归一化互相关系数；(x,y)为参考

像元位置；x—uy—u)为搜索像元位置；和s分

别为参考像元和搜索像元的像元值u和u分别为

参考和搜索窗口的像元平均值⑵].

3结果与分析

3.1冰川边界确定

冰川边界是确定冰崩灾害影响范围及碎屑流

物源的重要判断依据。我们首先利用Google

EarthTM光学影像基于目视解译的方式对加拉白垒

峰附近3处典型冰川边界进行了圈定。然而，因为

该区域常年云雾缭绕，降雨充沛，造成获取的光学

影像质量不高，所提取的研究区域的冰川并不一定

2.2偏移量追踪(OfftTracking)

技术

针对加拉白垒峰附近典型冰川碎屑流位移监

测，研究中采用了基于顾及地形起伏的SAR像素偏

移追踪技术[16].该技术是基于目标点的距离向及

方位向位置信息，利用最优解和互相关技术寻找2

幅SAR影像的同名点，并追踪雷达幅度影像中同名

特征目标点的位置变化来监测地表位移[17].基于

SAR影像的偏移追踪技术在冰川变化监测中有很

好的应用案例[18].偏移追踪技术在监测大量级地

表位移上具有较强的优势，但由于OfftTracking

技术是以像元大小来计算地表位移量，其精度比In-

SAR技术低一个数量级，一般可以达到SAR影像

像元分辨率的1/20〜1/10[19]。例如当图像分辨率

为3m时，测量运输协议合同 精度为15~30cm[20].

OfftTracking算法的基本思路是首先对地

表位移前后的2幅图像进行精配准，

并计算2幅影

像的整体偏移多项式；同时拟合并去除由于2幅图

像成像时间和空间不同所带来的系统误差偏移值;

最后根据像素偏移量和像元大小计算得出方位向和

距离向地表位移信息。

研究采用强度追踪法来计算SAR影像间的方

位向和距离向偏移量。强度跟踪法基于归一化互相

关(NNC,normalizedcrosscorrelation)准则对局部

偏移量进行准确估算。NCC依据像元进行计算，若

想获得亚像元配准精度，可以通过强度插值和相关

性插值2种不同的插值方法来实现，其中NCC的计

算准则公式为

)(s(x—u,y—u)—us)

,-----------------------------------------,(1)

(s(x—u

9y—u)—Us)2)1'2

是冰川的活动边界，不能获取与SAR影像获取时

间完全对应的冰川活动边界。为此，我们首先利用

光学影像进行初步圈定，其次利用SAR干涉影像

的相干性对基于光学影像圈定的冰川边界进行了

修正。

在InSAR相干性图中，冰川表面因雪消融、冰

川运动而呈现低相干性，可以利用该特性将其与周

围地物区别开来，进而提取冰川的边界[22].为此,

我们首先对重复轨道ALOS-2降轨数据进行粗配

准和精配准处理。因为ALOS-2数据量较少，为了

尽可能提取有效信息，我们选择了全组合形式，并计

算生成每幅干涉图的2幅SAR影像之间的相干性,

其计算公式⑺为

4

甘肃科学学报

2021年第2期

vEThTPEThTp]

其中Y代表相干系数，值域为(0,1),值越靠近1相

干性越好，相反则相干性越差；*表示复数共轭;

E[・]为数学期望;u，u为零均值的复高斯随机变

量.对进行全组合生成的干涉对计算相干系数，同

时生成相干系数图(简称“相干图”)。我们对所有相

干系数进行叠加处理，以获取稳定的相干系数结果,

并将其转换到地理坐标系下

，与获取的初始冰川边

界线进行叠加和辅助修正.最终结果如图2所示,

图2中1、2和3分别是色东普冰川

、岗普冰川和则

巴弄巴曲冰川边界.

2950'N

2945'N

2955'N

2

3

0^

2955'N

2950'N

2945'N

▲山峰

雅鲁藏布江■

□□冰川豳

05km

-,相干系数

9450'E9455'E

(b)相干图

9450'E9455'E

(a)GoogleEarthTM

光学影像数据

注：图例中0.9为相干系数值，代表相干性较好，在相干图中表现为黄色，具体见图例.

图中红色二角为加拉白垒峰，海拔

7294m,属于图幅内海拔最咼点.

图2研究区域冰川边界分布情况

Fig.2Distributionmap

ofglacierboundary

ofthestudyarea

3.2偏移量结果及分析

为了量化研究区域内冰川的位移量，

我们利用

瑞士GAMMA软件对ALOS-2数据根据2.2中的

原理与参数设置开展了偏移量追踪技术处理.首先

对SAR影像进行校正，并基于DEM进行精配准,

然后利用SAR影像强度基于强度互相关函数进行

估值,获取精确的距离向和方位向像素偏移量,并对

获取的像素偏移量结果进行了中值滤波，最终将所

获取的像素偏移量转换为距离向和方位向位移。

为了将SAR影像坐标系下的位移结果转换到

地理坐标系下，我们利用外部DEM对SAR强度图

进行了地理编码处理，获取地理坐标与SAR影像坐

标之间的转换关系，并将此转换关系应用到距离向

和方位向位移结果中，

从而获取了地理坐标系下距

离向和方位向的冰川表面位移时间序列结果.研究

区域典型冰川2017年2月2日—2018年10月25

日时间段在距离向、方位向上位移形变时间序列结

果分别如图3、图4所示.从图3中我们可以看到

岗普冰川随着时间的推移，冰川的变化主要集中在

冰川上端和尾端

，同时我们也发现岗普冰川在黑色

50

25

0

距离向

位移/m

(a)2017年2月2

日一2017年

7月

6

日在距离向位移

(b)2017年2月2日一2017年

11月9日在距离向位移

(c)2017年2月2

日一2018年

5月10

日在距离向位移

(d)2017

年2月2

日一2018年

10

月25

日在距离向位移

图3研究区2017年2月2日—2018年10月25日距离向位移时间序列

Fig.3Timeriesofthestudyareadistancedisplacementfrom2ndofFeb.2017to25thOct.2018

第

33卷董继红等：基于SAR偏移量跟踪技术的加拉白垒峰典型冰川位移监测

5

(a)2017年2月

2日一2017

年

(b)2017年2月2日一2017

年(c)2017年2月

2日一2018

年(d)2017

年2月2日一2018年

7月

6

日在方位向位移11月9

日在方位向位移5月10日在方位向位移

10月25日在方位向位移

图4研究区2017年2月2日至2018年10月25日方位向位移时间序列

Fig.4Timeries

ofthestudyarea

azimuthdisplacementfrom2ndFeb.2017to25thOct.2018

圈内变化复杂，经分析原因为该区域位于岗普冰川的

物质积累区域，形变量较为复杂

。

红线圈定区域位于

色东普冰川，

我们发现监测到的色东普冰川的形变主

要集中在该区域，且2018年5月10日—2018年10

月25日之间形变量最大，形变特征较复杂，原因可牡丹鹦鹉怎么养 能

与2018年10月17日发生在该区域的泥石流灾害有

关，导致位移特征复杂。同时比较图3和图4发现,

方位向的形变量更大，经分析原因是冰川流向近似于

方位向一致，与地形坡度

、坡向等联合导致。

同时，为了更好比较冰川随时间变化的位移特

征，我们选取岗普冰川提取了冰川位移的剖线(剖线

位置为图3(a)、图4(a)中AA，),该冰川在监测时间

段内方位向和距离向的累积位移时间序列如图5

所示。

剖线距离/km

(a)距离向

・2017年2月2日一2017年7月

6日・2017年2月2日一2017年11月9日

人2017举2月2日一

2018举5月10日v2017^2

月2

日一2018举10

月行文思路 25

日

60

40

20

0

-20

80

6400

4

000

06

4

800桅

5600g

1.7

244.8

剖线距离/km

(b)方位向

▼

2017年2月2日—2017年7月6日♦2017年2月

2日—2017

年11月

9日

2017举2月2日—2018^5月100►2017^2月2日—

2018^10月

25日

图5岗普冰川距离向和方位向剖线AA'时间序列

Fig.5SectionlineAAoftimeriesofdistanceand

azimuthdisplacementofthestudyarea

通过提取岗普冰川的剖线，结合图5我们可以

获得以下结果：

(1)距离向上：冰川位移具有较好的累积性，累

积位移在距离向呈现2个峰值，分别位于剖线AA

的1.7km和4.8km处。我们推测距离A点0〜

1.7km范围内的冰磧物随着冰川的位移在AA

1.7km处受局部地形的影响在此堆积，并形成第1

个峰值。而在4.8km处出现峰值(图5(a)距离A

点4.8km处虚线)的原因，一方面是上游物质在此

处运移堆积，另一方面是岗普冰川东侧有一处支流

冰川(香菇英文 图3(a)B点处)，该支流冰川汇集到冰川主流

上，使得在冰川主流上AA4.8km处上游冰磧物

源增加，并在AA4.8km处汇集形成第2个峰值。

(2)方位向上：我们发现在距离A点1km处,

2017年2月2日—2018年10月25日的累积形变

较前3个时间段的形变复杂，主要原因可能是时间

跨度较大，形变过大，造成地表特征爱的教育读后感 破坏，进而影响

获取结果的精度。同时我们也发现岗普冰川在方位

向距离A点4.8km处出现最大累积形变量，主要是

受其东侧支流冰川的影响，使得在距离A点4.8km

处上游冰磧物源累积，在630天的时间内最大累积

形变量超过了60mo

4结语

加拉白垒峰地区是冰、雪崩灾害发生频繁的区

域之一，开展该地区典型冰川的变化监测研究，对该

6

甘肃科学学报

2021年第2期

地区类似灾害的监测与防治具有重要的参考意义。

研究利用ALOS-2影像对加拉白垒峰附近的2个典

型冰川岗普冰川和色东普冰川的冰川位移量变化特

征进行监测。研究结果表明，基于InSAR相干性和

光学影像可以很好地确定冰川边界。位于加拉白垒

峰北侧的岗普冰川位移受地形及冰川支流汇集的影

响较大，该冰川在630天的时间内方位向(近南北)

最大累积位移量超60m.同时，加拉白垒峰南侧的

色东普冰川的局部位移受2018年发生的冰川泥石

流事件影响较大.

冰川研究内容丰富，如流速、体积变化以及冰物

质平衡等,研究中仅采用SAR数据对加拉白垒峰地

区的典型冰川位移开展了初步研究，下一步还需结

合光学偏移量技术和其他地形资料进行更深层次分

析和解译。

参考文献

：

:1]李志斐.气候变化对青藏高原水资源安全的影响[J].国际安全

研究,2018,36(3):42-63,157.

[2]

管伟瑾.西昆仑主峰地区现代冰川变化[D].兰州：兰州大学,

2020.

[3]李明月，孙学军，李胜楠，等.青藏高原及其周边地区冰川融水

径流无机水化学特征研究进展[].冰川冻土，2020,42(2)：

562-574.

:4]冀琴.1990—2015年喜马拉雅山冰川变化及其对气候波动的

响应[D].兰州:兰州大学，2018.

[5]吴珊珊，姚治君，姜丽光，等.现代冰川体积变化研究方法综述

[].地球科学进展,015,0(2):237-246.

[6]AltenaB.

ScambosT,

tingRecent

Short-termGlacierVelocityEvolutionOverSouthernAlaska

fromaLargeCollectionofLandsatData]J].CryosphereDis

cussions.2019.13(3)

:795-814.

:7]刘国祥

，

陈强，罗小军，等.InSAR原理与应用[M].北京：科学

出版社2019.

:8

]张波，张瑞，刘国祥，等.基于SAR影像的贡巴冰川末端冰湖

年际变化监测及溃决规律分析[J].

武汉大学学报：

信息科学

版,2019,44(7):1054-1064.

:9]周华云.基于SBAS-InSAR技术对五道梁多年冻土区地面形

变监测与分析[D].兰州：兰州交通大学,018.

[10]Gray免费铺货 AL

,MattarKE,VachonPW,esults

fromtheRADARSATAntarcticMapping

MissionData：

Esti

mation

ofGlacier

MotionUsingaSimpleRegistrationProce-

dure[C]/'/‘IEEEInternationalGeoscience&RemoteSensing

,19

98.

[11]StrozziT

,

acierMotionEsti

mationUsingSAROfft-trackingProcedures[J

].IEEE

TransactionsonGeoscienceand

RemoteSensing,2002,40

(11):2384-2391.

[12]刘国祥，张波，

张瑞，等.联合卫星SAR和地基SAR

的海螺沟

冰川动态变化及次生滑坡灾害监测[J].武汉大学学报：信息

科学版,2019,44(7):980-995.

[3]童立强，涂杰楠，裴丽鑫，等.雅鲁藏布江加拉白垒峰色东普流

域频繁发生碎屑流事件初步探讨[]•工程地质学报，2018,26

(6):1552-1561.

[4]刘传正，吕杰堂，童立强，等.雅鲁藏布江色东普沟崩滑-碎屑

流堵江灾害初步研究[J].中国地质

,019,6(2):219-234.

[15]Yang

Chengsheng,

Han

Bingquan

,ZhaoChaoying,ealCo-

and

Post-ismicDeformationMechanismsoftheMw7.3

Iran

Earthquake(2017)RevealedbySentinel-1InSARObrva

tions^

J].

RemoteSens.2019.11(4)：418.

[16]LiuX

.Zhao

C,

Zhang

ationoftheBaigeLand-

slide,Tibet,China,RevealedThroughtheIntegrationofCross-

PlatformALOS/PALSAR-1andALOS/PALSAR-2

SAROb-

rvations[J

].GeophysicalRearchLeters,2020,47(3).

[7]林昊.基于D-InSAR和OfftTracking技术的滑坡形变场提

取研究[D].北京：中国地质大学,014.

[18]ZhouY.

Li

Z,G-lacierMass

LossintheKarakoram

RegionDuringthe1970sto2000RevealedbyKH-9

Images

andSrtmDem[J].JournalofGlaciology,2017,63(238):331-

342.

[19]Haemmig(Asssmentof

GlacialLakeOutburstFloodsfromKyagarGlacier,Karako-

ramMountains,China[J].AnnalsofGlaciology,2014,55

(66):34-44.

[0]王群，张蕴灵，范景辉，等.利用高分三号影像监测依嘎冰川表

面运动

[J].武汉大学学报：信息科学版,020,5(3):460466.

[1]张晓博，赵学胜，葛大庆，等.利用新型C波段雷达卫星研究南

伊内里切克冰川运动特征[].武汉大学学报：信息科学版,

2019,4(3)：429-435.

[2]郭唯娜

，

可长生活习惯英语 青，范宇宾.基于SAR干涉数据的东帕米尔高原

冰川变化[J].地球信息科学学报,019,1(11):1790-1801.

第33卷

董继红等：基于SAR偏移量跟踪技术的加拉白垒峰典型冰川位移监测

7

TypicalGlacierDisplacementMonitoringofGyalaPeri

BadonSAROfftTrackingTechnology

DongJihong1,YangChengsheng1,ZhangBenhao1?ZhaoChaoying1'2

(ofGeologyEngineering

and

Geomatics,ChanganUniversity,Xian

710054,China;

eyLaboratoryof

Geo-information

Engineering,Xi'an710054,China)

AbstractWiththeincreasingshortageoffreshwaterresources?

glaciers?asanimportant,freshwaterre

source,ametime,thegeologicaldis

asters

caud

bythe

glaciermeltingalsodeeplyaffect,thesafetyofpeople'slifeandproperty.

Thispaper

takes

thetypicalglaciersin

theGyalaPeriareaofTibet,astherearchobject?lecting

theL-bandALOS-

2syntheticaperture

radarimagedataofJapan,andadoptingtheintensity-badofft,trackingtechnology

tomonitorandrearch

onthe

changesofglaciersdisplacement,ultsshow

that,themaximum

accumulationofGangpu

GlacieronthenorthsideofGyalaPeriexceeds60m

in630

days,andits

displacement,accumulation

andflow

velocity

aregreatlyaffectedbyitast,e-

vantrearchresultsofthispapercouldprovidereferencesfordisastermitigationandpreventionofsimilar

ice

avalanchedebris

flowdisastersinsouthernTibet.

KeywordsSyntheticapertureradar；GyalaPeri；Glacier；ng

（本文责编：冯婷）