武沧D-INSAR.docVIP

D-InSAR技术在矿山塌陷和地面沉降监测中的应用 张振生 （河北省遥感中心，石家庄 050021） 摘要：差分干涉测量（D-InSAR）技术，是合成孔径雷达（SAR）卫星应用的一个拓展。雷达图像的差分干涉图可用于监测cm级或更微小的地面形变，具有全天候、大面积监测地面沉降和矿山塌陷的优势，本文以武安矿山塌陷和沧州市地面沉降监测为例，介绍了这一新技术在灾害监测领域中的实际应用。 关键词：矿山塌陷；地面沉降；雷达图像；差分干涉测量 0 引言 上世纪90年代以来，欧美一些发达国家对机载和星载（包括航天飞机）合成孔径雷达理论和应用做了深入研究，获取了大量的商用SAR图像，D-InSAR技术在大地测量方面的研究和应用得到极大发展，一方面D-InSAR技术可以在大面积范围内（100km×100km）监测地面的微小形变，具有不需要人员进入灾害地区测量的特点，相对于GPSGLONASS全球定位系统和经典的大地测量手段（高精度水准测量）而言，有着巨大的优势和广阔的应用潜力。另一方面，高精度水准和GPS技术只能监测有限的、离散的控制点，而D-InSAR一幅图像就可以控制1万平方公里的地表形变监测，其空间分辨率可达到5m×20m。近年来基于D-InSAR的相应研究课题和项目多达上千项，在地震形变、山体滑坡、火山活动、矿区地面塌陷、城市地面沉降等人为因素造成的自然灾害监测方面，取得了一系列重要的研究成果。目前在太 空运行的雷达卫星有欧洲空间局的ERS－1/2、ENVISAT-1、日本的JERS-1和加拿大的RADARSAT等。其中ERS－1/2生成间隔至多相差一天的SAR数据,为InSAR的研究提供了全 球陆地面积75%的地形数据。本文以武安矿山塌陷和沧州地面沉降监测为例，介绍D-InSAR技术在该领域的应用。 雷达差分干涉测量处理方法 SAR数据处理的一般流程 要获取地表形变信息，则必须消除区域的地形信息。目前常采用的方法主要是以Massonnet为代表的DEME法（又称“二轨法”）和以Zebker为代表的三轨法。二轨法是利用实验区地表变化前后二幅影像生成干涉纹图，再利用事先获取的DEM数据模拟纹图，从干涉纹图中去除地形信息以得到地表变化信息。这种方法的优点是无须进行相位解缠（Massonnet，1993B），减少处理工作量。缺点是对于无DEM数据的地区无法采用上述方法，另外在引入DEM数据的同时有可能带入新的误差。三轨法是利用三景影像生成二幅干涉纹图，一幅反映地形信息（对于ERS-1/2 SAR数据来说，一般采用相隔一天的重复轨道数据），一幅反映地表形变信息，进行平地效应消除后，分别进行相位解缠，最后利用差分干涉原理计算得到地表形变信息。三轨法的主要优点是无需地面信息，数据间的配准较易实现。缺点是相位解缠的好坏将影响最终结果。 二轨法的主要步骤 ①根据SAR数据的时间间隔、季节等因素选择合适的图像对。 ②采用多级配准的方法对干涉相对精确配准。 ③基于局部地形坡度对干涉图像对进行滤波处理，然后生成干涉图。 ④计算DEM的点间隔与干涉纹图像元间隔之间的比值，对DEM点进行过采样，基于多普勒方程、斜距方程和椭球方程，利用轨道参数将DEM转换到雷达坐标系统，将DEM数据转换为相位值。 ⑤从干涉图中减去利用DEM模拟的干涉图得到形变图。 ⑥将形变图投影至地理坐标系中。 三轨法的一般处理流程 三次过境法差分InSAR处理流程见图1。 图1 三次过境法差分InSAR处理的工作流程 矿山地面塌陷监测 河北武安矿山地面塌陷监测采用的是2次和3次过境法差分InSAR处理，表3是用JERS-1和Envisat卫星SAR数据对武安矿区地面塌陷差分干涉测量综合结果。选择JERS-1 SAR 14景存档数据，用于测量92-98年时间段的地表沉降，同时购置了2003年以来的6景Envisat SAR数据，以分析现今地表沉陷。图形面积均为100km×100km。综合图反映了武安煤矿区近10年中地面沉降的分布范围和沉降量。 时相及数据组合 根据SAR数据的时相和相对基线距，将JERS-1 SAR数据分成10组2次差分分组和6组3次差分分组（见表1、2），以形成几种时间间隔。 表1 JERS-1数据2次差分分组 开始时间 结束时间 时间间隔（月） 垂直基线距（米） 930227 930412 1.5 1245 930412 940214 10 -44 950201 960416 14 -740 960416 960503 1.5 -1030 960826 970105 4.5 -720 970105 970517 4.5 -163 970517 980504 12 335 980504 980731 2.5 -393 980