留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于SBAS-InSAR的矿区地表沉降监测与分析

李达 邓喀中 高晓雄 牛海鹏

李达, 邓喀中, 高晓雄, 牛海鹏. 基于SBAS-InSAR的矿区地表沉降监测与分析[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1531-1537. doi: 10.13203/j.whugis20160566
引用本文: 李达, 邓喀中, 高晓雄, 牛海鹏. 基于SBAS-InSAR的矿区地表沉降监测与分析[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1531-1537. doi: 10.13203/j.whugis20160566
LI Da, DENG Kazhong, GAO Xiaoxiong, NIU Haipeng. Monitoring and Analysis of Surface Subsidence in Mining Area Based on SBAS-InSAR[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1531-1537. doi: 10.13203/j.whugis20160566
Citation: LI Da, DENG Kazhong, GAO Xiaoxiong, NIU Haipeng. Monitoring and Analysis of Surface Subsidence in Mining Area Based on SBAS-InSAR[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1531-1537. doi: 10.13203/j.whugis20160566

基于SBAS-InSAR的矿区地表沉降监测与分析

doi: 10.13203/j.whugis20160566
基金项目: 

国家自然科学基金 41272389

江苏高校优势学科建设工程资助项目 SZBF2011-6-B35

详细信息
    作者简介:

    李达, 硕士生, 主要从事InSAR技术的理论和应用研究。cumt_ld@163.com

    通讯作者: 邓喀中, 教授。kzdeng@cumt.edu.cn
  • 中图分类号: P258

Monitoring and Analysis of Surface Subsidence in Mining Area Based on SBAS-InSAR

Funds: 

The National Natural Science Foundation of China 41272389

the Priority Academic Program Development of Jiangsu Higher Education Institutions SZBF2011-6-B35

More Information
    Author Bio:

    LI Da, postgraduate, specializes in the theory and application of InSAR. E-mail: cumt_ld@163.com

    Corresponding author: DENG Kazhong, professor. E-mail:kzdeng@cumt.edu.cn
图(9)
计量
  • 文章访问数:  826
  • HTML全文浏览量:  91
  • PDF下载量:  289
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-05-08
  • 刊出日期:  2018-10-05

基于SBAS-InSAR的矿区地表沉降监测与分析

doi: 10.13203/j.whugis20160566
    基金项目:

    国家自然科学基金 41272389

    江苏高校优势学科建设工程资助项目 SZBF2011-6-B35

    作者简介:

    李达, 硕士生, 主要从事InSAR技术的理论和应用研究。cumt_ld@163.com

    通讯作者: 邓喀中, 教授。kzdeng@cumt.edu.cn
  • 中图分类号: P258

摘要: 针对雷达差分干涉测量技术(differential interferometry synthetic aperture radar,DInSAR)易受时空失相关、大气相位延迟等因素的影响,采用小基线集技术(small baseline subset,SBAS)对13景TerraSAR-X数据进行时序处理,估算并去除了残余数字高程模型误差、大气延迟误差和轨道误差,分析了研究区域2012-2013年沉降速率,发现2310和1301两工作面最大下沉速率分别为40 mm/a和50 mm/a;分析时序累计沉降值发现,2306、2308、2310这3个工作面在2012年11月15日之前地面沉降并不明显;在2310和1301工作面分别提取3个失相干现象表现较为缓慢的候选点进行时序分析,发现沉降值和时间成线性变化关系,且开采时间越早,其沉降特性越符合线性变化;将SBAS和DInSAR两种方法获得的累计沉降值进行差值分析,发现两种方法的差值在5 mm以内;在2310工作面走向方向和倾向方向选取若干观测点,并提取各观测点的时序沉降值对工作面时序沉降进行量化分析,实验表明SBAS-InSAR技术在矿区地表沉降监测与分析方面具有良好的应用前景。

English Abstract

李达, 邓喀中, 高晓雄, 牛海鹏. 基于SBAS-InSAR的矿区地表沉降监测与分析[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1531-1537. doi: 10.13203/j.whugis20160566
引用本文: 李达, 邓喀中, 高晓雄, 牛海鹏. 基于SBAS-InSAR的矿区地表沉降监测与分析[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1531-1537. doi: 10.13203/j.whugis20160566
LI Da, DENG Kazhong, GAO Xiaoxiong, NIU Haipeng. Monitoring and Analysis of Surface Subsidence in Mining Area Based on SBAS-InSAR[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1531-1537. doi: 10.13203/j.whugis20160566
Citation: LI Da, DENG Kazhong, GAO Xiaoxiong, NIU Haipeng. Monitoring and Analysis of Surface Subsidence in Mining Area Based on SBAS-InSAR[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1531-1537. doi: 10.13203/j.whugis20160566
  • 煤炭是我国的主体能源,长期占据着我国能源消费结构中65%以上,煤炭资源合理开发和利用对于中国经济发展和生态环境具有重要意义[1]。在地下煤炭资源被采出后,开采区周围原始应力平衡被打破,引起应力重新分布和覆岩移动,随着采空区扩大到一定的范围后,岩层移动就会向上波及地表,使地表发生沉降和变形[2]。在过去的几十年中,由于煤炭资源高强度、大面积的开采,使得开采所引起的地表沉陷及生态环境灾害问题日益突出,因此对于矿区开采沉陷的监测与分析具有重要的研究意义[3]

    雷达差分干涉测量技术(differential interferometry synthetic aperture radar, DInSAR)是利用对同一地区观测的多景SAR影像进行相干处理,获取高精度、高空间微小形变的信息技术[4]。DInSAR技术经过十几年的发展,已经在城市地面沉降监测、山体滑坡监测及矿区沉降监测等方面有了广泛的应用研究[5-7]。文献[8]探索DInSAR技术在矿区开采沉陷监测中的应用,证明DInSAR技术监测矿区地面沉降的可行性,但时空基线、大气效应等误差影响了DInSAR技术的测量精度[9],而小基线集技术(small baseline subset, SBAS)可以减轻时空失相关、大气延迟的影响,应用数量有限的影像得到毫米级的时序沉降量[10-12]。文献[13-15]已经利用SBAS技术对城市地面、高原冻土和长三角洲进行沉降监测研究,并取得较好的实验结果,但SBAS技术在煤矿地区的研究应用较少。煤矿地区地表散射特性不稳定,容易造成失相干,很难用传统的DInSAR方法实现长时间的形变监测,而用SBAS技术可以实现对矿区的长时间序列的形变监测。本文应用SBAS技术对13景TerraSAR-X数据进行处理,获得并分析了研究区域的沉降速率和时序累计沉降值,得出SBAS技术在矿区地表沉降监测与分析方面具有良好的应用前景。

    • 实验区域位于陕西省榆林市某矿区西北部,内有蒜皮滩、吕汗界等村庄,地表主要由沙漠和少量植被组成。实验区域包含4个工作面,编号分别为2310、2308、2306、1301(由于井下资料不全,为方便后续分析,本文将东边工作面编号为1301)。1301、2306、2308、2310工作面彼此相邻,采用综合机械化长壁垮落法开采,采深170 m,采厚2.5~3 m,各工作面间距40 m;1301工作面长965 m,宽240 m,开采时间为2012-04-2012-09;2306工作面长1 100 m,宽205 m,开采时间为2010-02-2010-10;2308工作面长1 126 m,宽245 m,开采时间为2010-12-2011-05;2310工作面长1 120 m,宽180 m,开采时间为2011-06-2012-04。根据开采沉陷规律,采深在100~200 m时,地表点移动的总时间为1~2 a,由此可知,1301工作面和2310工作面的地表会继续发生沉降。

    • 实验选用13景TerraSAR-X数据,影像时间跨度为2012-10-02-2013-03-05,数据为X波段,波长是3.1 cm, 影像的中心入射角为30.66°,方位向和距离向的分辨率分别为0.8 m和0.9 m,实验裁切影像大小为距离向为1 440像素,方位向为2 000像素,不进行多视处理,差分干涉测量中所用的DEM是30 m分辨率的SRTM数据。

      实验基于gamma软件实现DInSAR处理,基于StaMPS/MTI软件实现SBAS处理。DInSAR处理主要包括影像配准、去平地效应、地理编码、DEM相位差分等步骤。StaMPS/MTI首先根据幅度信息筛选出经过斜距向谱滤波后相位在短时间内失相干现象表现较为缓慢的SDFP(slowly decorrelation filter phase)候选点,然后根据SDFP候选点相位信息进行迭代运算,从中筛选出最终的SDFP。由于SDFP的精化选取是根据干涉相位进行的,干涉对的相干质量对SDFP的选取具有较大影响。StaMPS/MTI软件在组合干涉对时认为两景SAR影像的相干质量仅仅与时间/空间基线有关,并用简单的线性函数进行描述,而现实中影响相干质量的因素众多。为降低噪声等因素影响,在提取相干点目标之前需要检查每个干涉对的相干质量,去除相干性差的干涉对。

      基于13景影像,分别设置时间基线和垂直基线阈值为90 a和800 m,去掉其中相干性差的干涉对,最终挑选出33个干涉对用于时序分析,其中最大时间基线为88 a,最大空间基线为264 m,时间基线和空间基线详细信息见图 1。StaMPS/MTI处理首先设置振幅离差阈值为0.6,初步选取33万个SDFP点,然后根据相关系数迭代计算出相位误差,再根据相位误差筛选出31万个SDFP点,去掉那些由于相邻像元误差错误估计的SDFP点,最终筛选出14万个稳定的SDFP点,对这些SDFP点进行三维相位解缠,并在时间域和空间域上进行高通、低通滤波处理,最后分离出沉降相位、大气延迟相位、DEM误差相位。

      图  1  时空基线图

      Figure 1.  The Baseline Map of Time and Space

    • 通过上述数据处理后,获得2012-10-02-2013-03-05期间榆林某矿的雷达视线方向年平均沉降速率,根据公式Δhr/cosθr为卫星视线方向沉降,θ为入射角,Δh为垂直沉降值),将雷达视线方向的沉降速率换算到垂直方向。从图 2可以看出研究区域的房屋及道路等建筑物上面选取到较多的SDFP点,同时在4个工作面上也选取到大量的SDFP点,而其他时序方法很难在植被和沙漠地区选取到密集的SDFP点。结果显示该矿区4个工作面都有沉降,1301工作面沉降速率最大,沉降速率为30~50 mm/a,这是因为该工作面刚结束开采工作。其他3个工作面也均有沉降,2310工作面沉降速率为10~40 mm/a,2308和2306两个工作面沉降速率均为10~25 mm/a。

      图  2  研究区域垂直沉降速率

      Figure 2.  Vertical Subsidence Rate of the Study Area

    • 图 3是以2012-10-02为起始时间,其他时间相对起始时间的时序累计沉降值。从图 3中可以看出,在2012年10月-11月期间,1301工作面有5~8 mm的沉降量,其他3个工作面没有明显的沉降发生;从2012-11-15开始,2310工作面开始出现沉降值,并随着时间推移,沉降值不断变大,沉降范围也不断变大,到2013-03-05沉降值达到20 mm,说明2310工作面并没有达到稳定,沉降还在继续,1301工作面的累计沉降值达到30 mm,以后还会不断增大。

      图  3  研究区域时序累计沉降

      Figure 3.  The Timing Cumulative Subsidence of the Study Area

    • 为了分析研究区域的沉降性质,在2310工作面和1301工作面分别以100 m为间隔提取出6个观测点的时序沉降值,位置如图 4(b)所示,利用简单的线性拟合模型拟合出下沉值和时间的关系见图 5。从图 5可以看出,下沉值和时间呈线性关系,并发现开采时间越早的区域,对应SDFP点的时序沉降值越符合线性关系,如3号SDFP点。

      图  4  DInSAR和SBAS累计沉降值

      Figure 4.  The Cumulative Subsidence of DInSAR and SBAS

      图  5  研究区域沉降特征

      Figure 5.  The Settlement Characteristics of the Study Area

    • 由于研究区域没有实测值,为了评估SBAS结果的可靠性,对13景TerraSAR影像进行DInSAR处理,组成12个干涉对,对其进行影像配准、地理编码、去平地效应、DEM差分、相位解缠等处理,最终获得12个差分干涉图,在ArcGIS中将其进行累加处理,获得DInSAR处理的累计沉降图,如图 4(a)所示。在DInSAR、SBAS获得的累计沉降图走向方向提取观测点的沉降值,并去掉受到植被和噪声影响的观测点,采用折线图进行分析。

    • 图 6中可见,SBAS方法和DInSAR方法得到的累计沉降值在变化趋势上大体相同。2310工作面走向方向沉降值是一个由小变大再到由大变小的过程,两种方法在监测量级方面也较为接近,SBAS方法监测到最大沉降值是20 mm,DInSAR方法监测到最大沉降值是23 mm。为了更准确比较两种方法结果,将两种方法的结果进行作差处理,差值统计结果如图 7所示。虽然没有和实测值进行比较,但基于以往DInSAR方法的监测精度,可认为SBAS方法得到的结果是可靠的。

      图  6  DInSAR、SBAS沉降值

      Figure 6.  The Settlement of DInSAR and SBAS

      图  7  SBAS和DInSAR差值统计

      Figure 7.  Difference Statistics of SBAS and DInSAR

    • 图 4(b)是2012-10-02-2013-03-05期间研究区域的累计沉降值,1301工作面累计沉降值达到18~30 mm,2310工作面累计沉降值也达到10~20 mm。为了量化分析2310工作面的时序累计沉降值,在走向方向和倾向方向分别提取一系列点进行分析,去掉受植被和噪声影响的观测点,将走向和倾向方向的观测点画成折线图,如图 8图 9所示。

      图  8  走向时序累计沉降值

      Figure 8.  Timing Cumulative Subsidence of Trending

      图  9  倾向时序累计沉降值

      Figure 9.  Timing Cumulative Subsidence of Orientation

      图 8可知,随着时间的增加,2310工作面的沉降值不断变大,大的沉降区域沉降值由开始的3 mm,经过5个月时间沉降值增大到20 mm;小的沉降区域沉降值由开始的2 mm增加到12 mm。观察累计沉降变化曲线,可以看到2310工作面有一个快速下沉盆地,并以40 mm/a的速度沉降,其他下沉缓慢区域则以20 mm/a的速度沉降,这是因为2310工作面是由北向南开采,工作面南边区域比北边区域活跃,下沉速度较大,相同时间内累计沉降值也大,这一点在1301工作面也有很好验证,由图 4(a)可以看出1301工作面最南端区域累计沉降值比北边区域大。

      图 9也可以看出,2310工作面在倾向方向上表现出下沉盆地的特性,累计沉降值是先增大再减小的过程,下沉盆地中心沉降值从开始的1.5 mm以40 mm/a的平均速度增加到22 mm,说明2310工作面并没有稳定,未来沉降值还会不断地增大。

    • 针对传统的DInSAR方法时空失相关等问题,利用SABS方法对13景TerraSAR影像进行处理,得到研究区域的年平均沉降速率、时序累计沉降值等高空间分辨率的时间序列形变信息。通过比较分析SBAS和DInSAR两种方法获得的累计沉降值,表明两种方法监测结果接近,SBAS方法获得结果是可靠的。研究表明,2310和1301工作面时序沉降是线性沉降,最大沉降速度分别为40 mm/a和50 mm/a;在累计沉降图的走向和倾向方向提取一些观测点,通过这些观测点分析了2310工作面的时序沉降变化,观察到该工作面开采时间晚的区域比开采时间早的区域沉降速度大的情况。通过以上分析表明,SBAS技术可以为矿区地表的变形监测与分析提供新的监测手段。

参考文献 (15)

目录

    /

    返回文章
    返回