京津冀地区1992—2014年三阶段地面沉降InSAR监测

张永红，吴宏安，康永辉

中国测绘科学研究院，北京 100830



京津冀地区1992—2014年三阶段地面沉降InSAR监测

张永红，吴宏安，康永辉

中国测绘科学研究院，北京 100830

京津冀地区是我国地面沉降灾害发育最严重的地区，几乎每年都造成巨大的经济损失。本文提出了改进的时间序列InSAR技术——多主影像相干目标小基线干涉技术(MCTSB-InSAR)，利用4颗卫星摄取的3个时段的时间序列SAR影像：ERS-1/2 SAR(1992—2000年)、ENVISAT ASAR(2003—2010年)、RADARSAT-2(2012—2014年)，获取了京津冀地区1992—2014年间3个时段的地面沉降信息。经与京津两地120个以上水准测量数据进行比较，3个时段的地表监测结果的精度分别为8.7、4.7、5.4 mm/a。分析了北京和天津两市22年间地面沉降的时空变化特征，其中北京市地面沉降呈不断加重趋势；天津市地面沉降在1992—2010年间发展迅猛，在2010年以后有所减缓。同时，本文也表明MCTSB-InSAR技术是有效可靠的，在大区域地面沉降监测中具有推广应用价值。

地面沉降；InSAR；京津冀地区；地理国情监测

地面沉降是较为普遍的地质灾害，具有形成缓慢、持续时间长、影响范围广、成因机制复杂和防治难度大等特点。地面沉降已在我国20多个省(区、市)发生，每年造成数百亿元的直接经济损失，是影响经济社会发展的最重要的灾害因素之一。京津冀地区所在的华北平原是世界上最大的沉降漏斗[1]。对地面沉降进行精确、大范围、长时间监测是开展沉降控制工作的基础。以基岩标、分层标、 GNSS、水准点等为主要监测手段的地面沉降监测网络存在空间密度低、无法提供整个区域的沉降信息等问题，难以为防灾减灾提供可靠的决策依据。

始于1989年的差分合成孔径雷达干涉测量(differential interferometric SAR,D-InSAR)技术[2]曾在地表形变监测中展现了巨大潜力，因而受到广泛关注。但是，传统的两轨或三轨D-InSAR技术受到空间失相干、时间失相干和大气非均匀等因素[3-4]的限制，难以真正实用化。近10多年来得到迅速发展的时间序列InSAR技术是对传统D-InSAR技术的革新。它对时间序列干涉图进行分析，处理的对象不是整幅影像的全部像元, 而是其中具有稳定散射特性的高相干像元子集(下称“高相干点”)。时间序列InSAR技术总体上包括两类：以永久散射体干涉(permanent scatterer interferometry， PS-InSAR)[5-6]或者Persistent Scatterer Interferometry[7]为代表的单一主影像时间序列InSAR技术和以小基线集技术(small baseline subset interferometry, SBAS InSAR)[8]为代表的多主影像时间序列InSAR技术。时间序列InSAR技术已在火山、地震、滑坡、地面沉降等领域得到大量应用。

鉴于京津冀地区地面沉降的严峻局面，2014年1月国家测绘地理信息局在开展第一次全国地理国情普查工作的同时，启动了京津冀地区重要地理国情监测项目，地面沉降是其中的4项监测内容之一。目前，一些研究者已经开展过北京[9-10]、天津[11-13]、廊坊[14]、沧州[15-16]等地的地面沉降InSAR监测，但都是针对局部地区、较短时段的试验性研究。本研究利用1992年以来多颗卫星获取的SAR影像对京津冀大区域长时间跨度的地面沉降进行监测，面临着时间序列SAR数据集小、大区域点目标密度不均一、不同轨道沉降监测结果镶嵌等难题。因此，数据处理方法不能照搬现有方法。本文提出一种改进的时间序列InSAR技术，即多主影像相干目标小基线技术，对京津冀地面沉降进行遥感监测，并对地面沉降监测结果进行了精度评价，最后分析了京津两地22年间地面沉降时空变化的主要特点。

1　SAR影像数据源

根据不同合成孔径雷达卫星在京津冀地区的数据存档情况，分别选用ERS-1/2、ENVISAT、RADARSAT-2 SAR影像进行1992—2000年、2003—2010年、2012—2014年的地面沉降监测。所有影像均为单视复数据(SLC)格式。不同传感器影像的覆盖范围如图1所示。

ERS/ENVISAT SAR影像为条带扫描模式，距离向和方位向像元大小为7.90 m和3.98 m，单景影像的覆盖区域约为100 km×100 km。由于条件所限，只获取了7个Frame的ERS/ENVISAT SAR影像。轨道(Track)号分别为490、218、447、175。 这7景ERS/ENVISAT SAR影像的覆盖面积为5.03万km2，完整覆盖了北京市和天津市的平原地区，同时覆盖了河北省唐山市、廊坊市、保定市的大部分平原地区。

RADARSAT-2 SAR影像为宽模式，距离向和方位向像元大小分别为11.83 m和5.21 m，单景影像的覆盖面积约为150 km×150 km。本次监测使用了8景RADARSAT-2 SAR宽模式影像，其覆盖面积为9.36万km2，包含北京市和天津市的全部平原地区以及北京以南的河北省平原地区。为节省数据成本，每个RADARSAT-2单景覆盖区只获取了19期时间序列影像。京津冀全区域共获取了152景RADARSAT-2影像。

ERS-1/2、ENVISAT及RADARSAT-2时间序列SAR影像集的大小及获取的起始时间列于表1和表2中。可以看出，ERS-1/2 SAR影像数据集普遍较小，最小的只有11期影像。同时，ERS-1/2 卫星在其生命周期有不同的任务阶段(mission phases)，导致ERS-1/2 SAR序列影像的时间和空间基线很不均匀，造成ERS-1/2 SAR干涉失相干较严重，因此形变提取的精度不如ENVISAT ASAR影像的精度高。后面的监测结果也将印证这一点。

表1ERS-1/2 SAR及ENVISAT ASAR影像获取期数

Tab.1Image numbers and acquisition date of ERS-1/2 SAR and ENVISAT ASAR

Track/FrameERS-1/2SARENVISATASAR 期数获取时段期数获取时段218/2799211992-05-17—2000-06-09242003-06-18—2010-09-29447/2799121992-08-11—1998-07-10212004-06-18—2010-10-15175/2799131992-05-14—2000-05-21222003-06-15—2010-09-26490/2817121992-09-18—1998-08-17172003-08-05—2010-08-03218/2817151992-05-17—1999-09-02252003-06-18—2010-09-29447/2817111992-04-28—2000-10-03242003-10-17—2010-10-15175/2817181992-05-14—2000-06-25222003-06-15—2010-09-26

表2　RADARSAT-2影像获取期数

2　时间序列InSAR处理方法

本文提出了一种改进的时间序列InSAR技术，即多主影像相干目标小基线干涉技术(MCTSB-InSAR)。先利用MCTSB-InSAR技术进行单景覆盖区数据处理，获取单景影像覆盖区的地面沉降(形变)信息，包括平均沉降速率和时间序列沉降量。大区域通常需要多景SAR影像才能覆盖，因此需要将不同景，特别是不同轨道的沉降提取结果进行拼接，以得到大区域覆盖完整且一致的沉降结果。MCTSB-InSAR技术能实现多景形变结果的平差及镶嵌处理。这是MCTSB-InSAR方法区别于其他方法的一个优点。除文献[17]外，多景形变结果的连接问题很少被讨论。

2.1单景覆盖区的时间序列InSAR处理

MCTSB-InSAR技术中的单景覆盖区时间序列处理综合了PS-INSAR和SBAS InSAR方法的优点，克服了相应的缺陷。它不以单一SAR作为全部干涉像对的固定主影像，而是考虑满足小基线要求的时间序列SAR数据集任意组合的干涉对，因此不像PS-INSAR方法需要25景以上的SAR影像才能进行有效的分析。MCTSB-InSAR方法也不像SBAS方法需要先完成每个差分干涉图的相位解缠，然后再进行时间序列分析，而是通过考察空间邻近的两个高相干点的二次相位差分，同时完成形变参数解算和相位解缠。一般情况下，MCTSB-InSAR将形变分解为线性形变(主形变)加上非线性形变(次形变)。如果形变呈现强烈的非线性，则可将形变分解为高阶多项式形变加上残余形变[18]。 MCTSB-InSAR方法中的单景覆盖区时间序列处理主要包括以下5个方面。

2.1.1小基线干涉组合

对SLC格式的SAR影像完成精确配准后，根据预先设定的关于时间基线和空间基线的阈值条件自由组合生成干涉图。这种方式在处理小数据集时尤为重要，因为可以使用较小数目的SAR影像得到较多的干涉图。MCTSB-InSAR方法要求这些干涉图组成一个连通网络。如果这些干涉图形成了两个或多个小基线干涉图集，则需要手动添加几个干涉组合将这些分离的小基线集连通起来。这也是与SBAS方法的一个不同之处。因为SBAS方法要求先对干涉图进行相位解缠，然后通过奇异值分解技术解决孤立的干涉图集之间相位的一致性问题。MCTSB-InSAR方法不要求事先对干涉图进行相位解缠，孤立的干涉图容易造成形变模型参数求解不唯一。

2.1.2高相干点提取

高相干点指干涉相干性较高的像素，既包括具有优势散射体的点目标(permanent scatterer or persistent scatterer，PS)，也包括能在一段时间内保持较高相干性的分布式目标。对点目标和高相干分布式目标可以设置不同的参数及阈值条件分别提取，然后求其点集的并集作为最终的高相干点。提取PS点的参数有幅度离差[6]和平均振幅。提取高相干分布式目标的参数有平均相干系数[19]等。大区域监测中将点目标和高相干分布式目标都包括进来非常重要，因为城市区域点目标密集，但是非城区点目标很少。如果不考虑具有较高相干性的分布式目标，在非城区提取的高相干点密度不足，将难以建立高相干点有效的网络连接。

2.1.3高相干点网络连接

通常的时间序列InSAR方法[19-20]都是采用整体的Delaunay三角网将高相干点连接在一起。如果研究区仅限于城市区域，由于高相干点密度很高，整体Delaunay三角网连接比较稳健。但是，如果研究区包含很大的非城区范围(例如整景影像)，由于非城区高相干点密度不高，或者高相干点质量不高，有些连接将会在后续处理中被剔除。这样，整体Delaunay三角网连接由于冗余度不够，在剔除不满足要求的连接边后极有可能分裂成许多不连通的子网络，给后面的积分操作带来困难。为了增加网络连接的冗余度，MCTSB-InSAR技术采用了局部Delaunay三角网连接方法[21]，即将整景影像分割互相重叠的网格，在每个网格中采用Delaunay三角网连接。最终连接的冗余度依赖于网格大小和相邻网格的重叠度，网格越小，重叠度越大，整体冗余度越大。为了减弱连接边上大气相位的影响，要求连接边的长度不大于大气相关距离(根据SAR影像分辨率不同，通常设为1000～2000 m)。

2.1.4线性形变求解

(1)

在完成高相干点的网络连接后，第k个干涉图上由相邻高相干点i,j构成的连接边的干涉相位差可以表示为

(2)

(3)

式中，J表示单位复数。求解未知数Δhi,j、Δvi,j的方法有二维周期图法[5]或者空间搜索法[18,21]，二者是等价的。当对所有的边完成最大化求解后，选择γi,j≥0.7的边作为可靠的连接，那些没有被可靠边连接的高相干点将被剔除。最后，从某一参考点开始对连接边上解算出的速度和高程误差增量进行积分，将得到各相干点上形变速率和高程误差的绝对量。

2.1.5非线性形变估计

从差分干涉相位中去除高相干点的线性形变相位和高程误差相位，得到高相干点的残余相位，包括大气影响相位、非线性形变相位以及噪声相位。大气影响相位在时间序列上是不相关的，为高频信号；在空间分布是相关的，为低频信号；非线性形变相位在时间序列上是低频信号，在空间上是不相关的，为高频信号；噪声相位则是时间和空间都不相关的随机高频信号。利用残余相位3个分量的不同的时空频谱特征，可以将三者分离出来[5]。将线性形变和非线性形变相位联合起来就可得到每个影像获取时刻的累积地表形变。假定京津冀地表形变主要体现为地面沉降，基于雷达波入射角关系，将获取的雷达视线向形变量转换为垂直方向上的地面沉降值[22-23]。

2.2多景覆盖区的地面沉降结果平差与镶嵌

单景影像覆盖区的时间序列SAR影像的起始时间通常不一致，累积沉降量自然不同。 因此累积沉降量难以镶嵌。将多景覆盖区的沉降速率结果进行平差，可以获得大区域的平均沉降速率镶嵌图。

同名高相干点的识别是平差及镶嵌的基础。虽然利用影像匹配技术能通过相邻SAR振幅影像的匹配完成同名高相干点的自动识别，但耗时太长且易出错。文献[18]采用将重叠区的不同轨SAR影像转化到一个统一的坐标系统中，并将其中某轨SAR影像坐标作为统一的坐标系，在重叠区高差不大的情况下，可以通过多项式将其他轨SAR影像坐标变换到这个统一坐标系。显然，如果重叠区高差大，这种做法难以确认同名高相干点。事实上，最好的统一坐标系是地理编码后得到的大地坐标系(此时，高程导致的不同轨像元定位差异已经得到改正)，利用地理编码后同名高相干点有相同大地坐标这一原理即可有效解决这一问题。必须注意的是所有SAR影像的地理编码应该使用同种类型的DEM，如SRTM DEM，这样能保证相邻景SAR影像地理编码后有相同的平面精度。

理论上，用一个偏移值可足够说明不同轨所得形变结果的差异[18]。文献[18]采用了距离和方位向的一次多项式来说明这个差异。这种做法极易将形变场中可能存在的线性趋势分量吸收，当成是不同轨上由于轨道误差、大气误差等造成的InSAR形变差异，反而造成平差后更大的不一致。另外，文献[17]讨论的是雷达视线向形变结果的连接，真实的形变投影到不同轨的雷达视线向后可能存在沿距离和方位向线性变化的趋势性差异。本文已经将视线向形变转化为垂向沉降，消除了不同视角投影的影响，因此在进行多轨沉降结果的平差时只考虑了偏移(即零次多项式)。

平差方程由两类构成。第1类方程是根据同名高相干点的形变速率一致的准则建立的，形如式(4)

Vi(x,y)+Δi=Vj(x,y)+Δj

(4)

式中，Vi(x,y)、Vj(x,y)分别代表相邻单景覆盖区(第i景和第j景)同名高相干点(x,y)上的形变速率值；Δi、Δj为待求的第i景和第j景的形变速率偏移值。第1类方程进一步可写成

Δi-Δj=Vj(x,y)-Vi(x,y)

(5)

第2类方程是将水准(或GPS)测量结果当成控制点构成的，如式(6)

Vi(x,y)+Δi=GCP_V(x,y)

(6)

式中，GCP_V(x,y)表示在点(x,y)处经由水准测量得到的形变速率值。将如式(5)和式(6)的所有方程联立在一起，即可解得所有单景覆盖区InSAR形变速率的偏移值。在有控制点的情况下，平差求得的偏移值同时也是每个单景覆盖区形变结果的绝对定标值，这样同时完成了每景覆盖区形变速率的平差与绝对定标。如果没有控制点，则没有第2类方程。通过平差将会得到以某一个单景覆盖区形变结果为参考的相对偏移值。最后，将经偏移值改正后的同名点的形变速率取算术平均，即可实现整个大区域的形变速率的镶嵌。

3　京津冀地区22年地面沉降监测结果分析

根据上述方法，提取了京津冀地区1992—2000年、2003—2010年、2012—2014年的地面沉降速率和累积沉降量等信息，并得到了3个时段地面沉降速率的大区域镶嵌结果，如图2—图4所示。在大区域沉降结果平差之前，已经对沉降结果进行了地理编码。ERS-1/2和ENVISAT沉降结果地理编码的输出像元大小为40m，RADARSAT-2沉降结果地理编码的输出像元大小为30m。图2—图4中负值表示地面下沉，正值表示地面抬升，负值越小表示沉降越严重。从图2可以看到ERS-1/2影像得到的大区域镶嵌图上存在一些不一致，特别是T175/F2799与T175/F2817之间，这主要是由上述的ERS-1/2SAR数据集小、时空基线差异大等因素造成的。相比之下，图3和图4的大区域沉降速率镶嵌结果则有很好的一致性。这说明本文的大区域形变结果的平差与镶嵌方法是有效的。

为了标定和验证InSAR监测结果，从北京市和天津市测绘部门获取了相近时段的水准数据，包括：天津市1993—2000年184个观测数据，用于ERS-1/2结果定标(即平差控制)和精度验证；北京市2005—2010年74个观测数据、天津市2003—2010年175个观测数据，用于ENVISATASAR结果定标和精度验证；北京市2012—2014年26个观测数据、天津市2011—2013年185个观测数据，用于RADARSAT-2结果定标和精度验证。其中1/3数据用于平差中的控制，2/3数据用于精度验证。对RADARSAT-2结果进行定标和精度验证的水准点分布如图5所示，其中蓝色点为平差控制点，共70个，红色点为精度验证点，共141个。表3为各时段InSAR沉降速率的精度验证结果。可以看出，除1992—2000年时段的沉降速率精度略差外，ENVISAT和RADARSAT-2获取的沉降速率精度约在5mm/a。

表3　各时段InSAR沉降结果精度验证分析

从图4中可以清晰地看到整个京津冀三地近3年的沉降中心分布。京津冀地区地面沉降非常严重，三地沉降漏斗有连接成片的趋势。

由于3个时段的监测区基本完整地包括了京津两市的平原地区，笔者对京津两市1992—2014年的地面沉降(图6和图7)的时空变化特点进行了进一步分析。在分析沉降区面积时，首先采用克里金(Kriging)插值方法将高相干点上的沉降速率值内插至每个像元，然后进行面积统计。

图6所示，北京平原地区地面沉降1992—2000年处于初步发展阶段，零星分布一些较小的沉降中心，最大沉降速率为48 mm/a，位于门头沟区龙泉镇。经查实，该沉降主要由龙泉镇煤矿地下采煤引起(该煤矿2000年被关闭)。2003—2010年北京地区地面沉降快速发展，在东部朝阳-通州形成了沉降带，在北部海淀、昌平、顺义等地形成了5个散落分布的沉降中心，沉降速率最大达143 mm/a，位于朝阳区东坝乡(门头沟煤矿的沉降已经停止)；沉降速率超过50 mm/a的面积为265.41 km2。2012—2014年东部沉降带继续发展，沉降范围和速率均继续增大；北部散落分布的沉降中心继续扩展，连接成片，形成海淀—昌平—顺义沉降带；最大沉降速率达152 mm/a，仍然位于朝阳区东坝乡；沉降速率超过50 mm/a的面积达到433.25 km2。整体来看，在过去3个时段，北京市地面沉降呈不断加重的趋势。

从图7可以看出，天津市1992—2000年地面沉降已有一定规模，最大沉降速率为81 mm/a，位于汉沽区，沉降速率超过50 mm/a的面积达到85.73 km2。2003—2010年，天津市地面沉降的速率和范围都在迅速发展，在西北部、南部、东南部形成了3个大的沉降带；最大沉降速率为135 mm/a，位于静海区；沉降速率超过50 mm/a的面积增长到2749.1 km2。2012—2014年，整体上呈现“北增南减”的格局，北部武清、北辰和汉沽等地区沉降速率增加，南部地区沉降均有减缓；最大沉降速率达153 mm/a，位于武清区南部的王庆坨镇；沉降速率超过50 mm/a的面积缩小到1 117.55 km2。总体上，天津市地面沉降在2010年以后有所减缓。

注：较小框表示ERS/ENVISAT SAR影像覆盖范围，较大框表示RADARSAT-2 SAR影像覆盖范围，背景为SRTM DEM图1　京津冀地区地面沉降监测SAR影像覆盖范围Fig.1　Coverage of SAR data for subsidence monitoring in Beijing-Tianjin-Hebei region

图2　1992—2000年京津冀地区平均沉降速率Fig.2　Average subsidence velocity in Beijing-Tianjin-Heibei region from 1992 to 2000

图3　2003—2010年京津冀地区平均沉降速率Fig.3　Average subsidence velocity in Beijing-Tianjin-Heibei region from 2003 to 2010

从笔者收集的北京地区地下水流场监测数据、基础水文地质数据等资料来看，北京地区地面沉降分布范围与地下水降落漏斗分布在空间位置上基本一致，说明长期超量开采地下水是导致北京地面沉降发生的主要原因。

图4　2012—2014年京津冀地区平均沉降速率Fig.4　Average subsidence velocity in Beijing-Tianjin-Heibei region from 2012 to 2014

图5　用于RADARSAT-2沉降结果平差控制(蓝点)和精度验证(红点)的211个水准点分布以及141个精度验证点上的误差图Fig.5　Locations of 211 leveling points used for adjustment control and accuracy evaluation of RADARSAT-2 derived subsidence result(a). Blue points are for adjustment control, and red points for accuracy evaluation. The error of InSAR derived subsidence velocity over the 141 red points are shown in(b)

图6　北京市3个时段沉降速率对比图Fig.6　Comparison of subsidence velocity in three stages in Beijing

图7　天津市3个时段沉降速率对比图Fig.7　Comparison of subsidence velocity in three stages in Tianjin

4　结　论

本文重点介绍了京津冀大区域1992—2014年地面沉降InSAR监测的技术方法和结果。该监测结果首次提供了京津冀地区过去22年间地面沉降的完整、可靠信息，对于京津冀地区地面沉降灾害防治、规划制定、水资源合理使用、经济布局调整等具有重要的参考价值。研究成果同时表明，本文提出的大区域时间序列InSAR形变监测技术MCTSB-InSAR是有效的，可在未来全国大范围的地面沉降监测中推广应用。

致谢： 特别感谢欧空局提供的ERS-1/2 SAR及ENVISAT ASAR数据，加拿大MDA公司提供的RADARSAT-2数据，北京市测绘设计研究院、北京市水文地质工程地质大队、天津市测绘院提供的水准测量数据。

[1]何庆成, 刘文波, 李志明. 华北平原地面沉降调查与监测[J]. 高校地质学报, 2006, 12(2): 195-209. HE Qingcheng, LIU Wenbo, LI Zhiming. Land Subsidence Survey and Monitoring in the North China Plain[J]. Geological Journal of China Universities, 2006, 12(2): 195-209.

[2]GABRIEL A K, GOLDSTEIN R M, ZEBKER H A. Mapping Small Elevation Changes over Large Areas: Differential Radar Interferometry[J]. Journal of Geophysical Research: Solid Earth, 1989, 94(B7): 9183-9191.

[3]ZEBKER H A, VILLASENOR J. Decorrelation in Interferometric Radar Echoes[J]. IEEE Transactions on Geoscience and Remote Sensing, 1992, 30(5): 950-959.

[4]DING Xiaoli, LI Zhiwei, ZHU Jianjun, et al. Atmospheric Effects on InSAR Measurements and Their Mitigation[J]. Sensors, 2008, 8(9): 5426-5448.

[5]FERRETTI A, PRATI C, ROCCA F. Nonlinear Subsidence Rate Estimation Using Permanent Scatterers in Differential SAR Interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(5): 2202-2212.

[6]FERRETTI A, PRATI C, ROCCA F. Permanent Scatterers in SAR Interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(1): 8-20.

[7]HOOPER A, ZEBKER H, SEGALL P, et al. A New Method for Measuring Deformation on Volcanoes and other Natural Terrains Using InSAR Persistent Scatterers[J]. Geophysical Research Letters, 2004, 31(23): L23611.

[8]BERARDINO P, FORNARO G, LANARI R, et al. A New Algorithm for Surface Deformation Monitoring Based on Small Baseline Differential SAR Interferograms[J]. IEEE Transactions on Geoscience and Remote Sensing, 2002, 40(11): 2375-2383.

[9]宫辉力, 张有全, 李小娟, 等. 基于永久散射体雷达干涉测量技术的北京市地面沉降研究[J]. 自然科学进展, 2009, 19(11): 1261-1266. GONG Huili, ZHANG Youquan, LI Xiaojuan, et al. Research on the Land Subsidence of Beijing Based on PS-InSAR Technology[J]. Progress in Natural Science, 2009, 19(11): 1261-1266.

[10]陈蓓蓓, 宫辉力, 李小娟, 等. 基于InSAR技术北京地区地面沉降监测与风险分析[J]. 地理与地理信息科学, 2011, 27(2): 16-20. CHEN Beibei, GONG Huili, LI Xiaojuan, et al. Monitoring and Risk Analysis of Land Subsidence in Beijing Based on Interferometric Synthetic Aperture Radar (InSAR) Technique[J]. Geography and Geo-Information Science, 2011, 27(2): 16-20.

[11]范景辉, 郭华东, 郭小方, 等. 基于相干目标的干涉图叠加方法监测天津地区地面沉降[J]. 遥感学报, 2008, 12(1): 111-118. FAN Jinghui, GUO Huadong, GUO Xiaofang, et al. Monitoring Subsidence in Tianjin Area Using Interferogram Stacking Based on Coherent Targets[J]. Journal of Remote Sensing, 2008, 12(1): 111-118.

[12]LUO Qingli, PERISSIN D, LIN Hui, et al. Subsidence Monitoring of Tianjin Suburbs by TerraSAR-X Persistent Scatterers Interferometry[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 7(5): 1642-1650.

[13]LIU Guoxiang, JIA Hongguo, NIE Yunju, et al. Detecting Subsidence in Coastal Areas by Ultrashort-baseline TCPInSAR on the Time Series of High-resolution TerraSAR-X Images[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(4): 1911-1923.

[14]葛大庆, 王艳, 郭小方, 等. 利用短基线差分干涉纹图集监测地表形变场[J]. 大地测量与地球动力学, 2008, 28(2): 61-66.

GE Daqing, WANG Yan, GUO Xiaofang, et al. Surface Deformation Field Monitoring by Use of Small-baseline Differential Interferograms Stack[J]. Journal of Geodesy and Geodynamics, 2008, 28(2): 61-66.

[15]葛大庆, 王艳, 郭小方, 等. 基于相干点目标的多基线D-InSAR技术与地表形变监测[J]. 遥感学报, 2007, 11(4): 574-580.

GE Daqing, WANG Yan, GUO Xiaofang, et al. Surface Deformation Monitoring with Multi-baseline D-InSAR Based on Coherent Point Target[J]. Journal of Remote Sensing, 2007, 11(4): 574-580.

[16]张玲, 葛大庆, 郭小方, 等. 近十年来沧州地区地面沉降演化状况[J]. 上海国土资源, 2014, 35(4): 72-75, 80.

ZHANG Ling, GE Daqing, GUO Xiaofang, et al. Land Subsidence in Cangzhou over the Last Decade Based on Interferometric Time Series Analysis[J]. Shanghai Land & Resources, 2014, 35(4): 72-75, 80.

[17]KETELAAR G, VAN LEIJEN F, MARINKOVIC P, et al. Multi-Track PS-InSAR: Datum Connection and Reliability Assessment[C]∥Proceedings of Envisat Symposium 2007. Montreux, Switzerland: ESA Communication Production Office, 2007.

[18]张永红, 吴宏安, 孙广通. 时间序列InSAR技术中的形变模型研究[J]. 测绘学报, 2012, 41(6): 864-869. ZHANG Yonghong, WU Hong’an, SUN Guangtong. Deformation Model of Time Series Interferometric SAR Techniques[J]. Acta Geodaetica et Cartographica Sinica, 2012, 41(6): 864-869.

[19]MORA O, MALLORQUI J J, BROQUETAS A. Linear and Nonlinear Terrain Deformation Maps from a Reduced Set of Interferometric SAR Images[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(10): 2243-2253.

[20]HOOPER A, SEGALL P, ZEBKER H. Persistent Scatterer Interferometric Synthetic Aperture Radar for Crustal Deformation Analysis, with Application to Volcán Alcedo, Galápagos[J]. Journal of Geophysical Research: Solid Earth, 2007, 112(B7): B07407.

[21]ZHANG Lei, DING Xiaoli, LU Zhong. Modeling PSInSAR Time Series without Phase Unwrapping[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(1): 547-556.

[22]张永红, 张继贤, 龚文瑜, 等. 基于SAR干涉点目标分析技术的城市地表形变监测[J]. 测绘学报, 2009, 38(6): 482-487. DOI: 10.3321/j.issn:1001-1595.2009.06.003.ZHANG Yonghong, ZHANG Jixian, GONG Wenyu, et al. Monitoring Urban Subsidence Based on SAR Interferometric Point Target Analysis[J]. Acta Geodaetica et Cartographica Sinica, 2009, 38(6): 482-487. DOI: 10.3321/j.issn:1001-1595.2009.06.003.

[23]吴宏安, 张永红, 陈晓勇, 等. 基于小基线DInSAR技术监测太原市2003～2009年地表形变场[J]. 地球物理学报, 2011, 54(3): 673-680.

WU Hong’an, ZHANG Yonghong, CHEN Xiaoyong, et al. Ground Deformation Monitoring Using Small Baseline DInSAR Technique: A Case Study in Taiyuan City from 2003 to 2009[J]. Chinese Journal of Geophysics, 2011, 54(3): 673-680.

(责任编辑：张艳玲)

Ground Subsidence over Beijing-Tianjin-Hebei Region during Three Periods of 1992 to 2014 Monitored by Interferometric SAR

ZHANG Yonghong,WU Hongan,KANG Yonghui

Chinese Academy of Surveying and Mapping, Beijing 100830, China

The Beijing-Tianjin-Hebei region suffers the most serious ground subsidence in China, which has caused huge economic losses every year. Therefore, ground subsidence was listed as an important mission in the project of geographic conditions monitoring over Beijing-Tianjin-Hebei launched by the National Administration of Surveying, Mapping and Geoinformation in 2013. In this paper, we propose a methodology of ground subsidence monitoring over wide area, which is entitled “multiple master-image coherent target small-baseline interferometric SAR (MCTSB-InSAR)”. MCTSB-InSAR is an improved time series InSAR technique with some unique features. SAR datasets used for ground subsidence monitoring over the Beijing-Tianjin-Hebei region include ERS-1/2 SAR images acquired between 1992 to 2000, ENVISAT ASAR images acquired between 2003 to 2010 and RADARSAT-2 images acquired between 2012 to 2014. This research represents a first ever effort on mapping ground subsidence over Beijing-Tianjin-Hebei region and over such as a long time span in China. In comparison with more than 120 leveling measurements collected in Beijing and Tianjin, the derived subsidence velocity has the accuracy of 8.7 mm/year (1992—2000), 4.7 mm/year (2003—2010), and 5.4 mm/year (2012—2014) respectively. The spatial-temporal characteristics of the development of ground subsidence in Beijing and Tianjin are analyzed. In general, ground subsidence in Beijing kept continuously expanding in the period of 1992 to 2014. While, ground subsidence in Tianjin had already been serious in 1990s, had dramatically expanded during 2000s, and started to alleviate in recent years. The monitoring result is of high significance for prevention and mitigation of ground subsidence disaster, for making development plan, for efficient and effective utilization of water resource, and for adjustment of economic framework of this region. The result also indicates the effectiveness and reliability of the MCTSB-InSAR method. Thus, the MCTSB-InSAR method is applicable to monitoring ground subsidence over large areas in the future.

ground subsidence; interferometric synthetic aperture radar; Beijing-Tianjin-Hebei region; geographic conditions monitoring

The National Basic Research Program of China (973 Program) (No. 2012CB719905); The National Natural Science Foundation of China (Nos. 41271430;41304010); The National Geographic Conditions Monitoring program, NASG

ZHANG Yonghong (1973—), male, PhD, research fellow, PhD supervisor, majors in SAR interferometry and deformation monitoring, SAR polarimetry and land cover/land use mapping.

张永红，吴宏安，康永辉.京津冀地区1992—2014年三阶段地面沉降InSAR监测[J].测绘学报，2016,45(9)：1050-1058.

10.11947/j.AGCS.2016.20160072.

ZHANG Yonghong,WU Hongan,KANG Yonghui.Ground Subsidence over Beijing-Tianjin-Hebei Region during Three Periods of 1992 to 2014 Monitored by Interferometric SAR[J]. Acta Geodaetica et Cartographica Sinica,2016,45(9):1050-1058. DOI:10.11947/j.AGCS.2016.20160072.

P237

A

1001-1595(2016)09-1050-09

国家973计划课题(2012CB719905)；国家自然科学基金(41271430；41304010)；国家测绘地理信息局全国地理国情监测项目

2016-02-25

张永红(1973—)，男，博士，研究员，博士生导师，研究方向为SAR干涉测量形变监测、极化信息提取和土地利用/土地覆盖。

E-mail： yhzhang@casm.ac.cn

修回日期： 2016-07-21