基于卫星雷达数字影像矿区地表沉降研究.docVIP

基于卫星雷达数字影像矿区地表沉降研究 　　摘要：利用卫星遥感技术完成矿区地表沉降情况的动态监测，卫星遥感技术在矿区地表沉降监测中的应用研究，不仅要利用卫星勘查矿区地表的特征表现，更要解决卫星遥感监测技术的精度问题，为此，本文将重点阐述基于卫星雷达数字影像的矿区地表沉降监测技术。 　　关键词：卫星雷达；数字影像；矿区；沉降 　　中图分类号： [TU196.2] 文献标识码： A 文章编号： 　　矿产资源的过度开采容易破坏矿区地下地质结构，引发地表沉降区、滑坡、崩塌、泥石流、地裂缝、地面沉降和地面积水等事故，甚至导致矿区透水事故，造成严重的安全隐患甚至人员伤亡。通过对矿区长期、动态的沉降监测，对矿区沉降规律进行分析，可及时掌握矿区地质环境的破坏程度，从而针对实际情况制定相应应急处理方案。同时，矿区沉降监测获取到的地表形变数据信息可为矿区的规划、治理和发展提供数据支撑，有利于促进矿山地质灾害环境的一体化管理。因此，矿区沉降监测是矿区开采和可持续发展的重要组成部分，对预防潜在地质灾害具有重要意义。以下将探讨一种卫星雷达数字影像应用领域中的地面沉降监测方法[1]。 　　1卫星雷达数字影像的基本原理 　　基本原理是利用同一矿区范围的多幅卫星雷达数字影像组成干涉图，干涉图中不仅包含因地形起伏引起的干涉相位，还包含由地表位移引起的形变相位，通过对干涉图进行差分处理(去除地形影响)来获取地表微量形变。对于一幅干涉图而言，形变相位的干涉影像是通过已配准的主辅卫星雷达数字影像的共轭相乘所得到的： 　　 　　式中：(r， n)表示像素坐标；M和S分别为复数形式的主辅卫星雷达数字影像；IF为复数形式的干涉影像，而其相位Δφ可以写成以下的形式： 　　 　　式中：λ 表示雷达波长；R 为斜距；θ 为雷达入射角；k 为整周模糊度；B⊥ 和B//分别为干涉对的垂直基线和平行基线；φnoise代表相位噪声；φatm为大气相位。 　　等号右侧第一项表示形变相位，是指在主辅影像所间隔的时间范围内，地表在雷达视线方向上的位移D所引起的干涉相位。第二项表示由地表高程H引起的相位变化，即地形相位。去除地形相位的方法有两种：一种是利用仅包含地形信息的卫星雷达数字影像对生成干涉图进而得到地形相位；另一种是利用已有的外部DEM模拟地形信息从而实现地形相位的去除。因此差分干涉方法可归结为两种模式：使用外部DEM的二轨方法；使用多图像干涉生成地形影像对的三轨或四轨方法。第三项表示平地相位，是由不同像素在参考椭球面的斜距差所引起的。B//可通过卫星精密轨道数据求得。大气相位(φatm)可以使用干涉图堆叠的方法来去除。相位噪声(φnoise)包括热噪声、时空失相关噪声、多普勒失相关噪声以及数据处理过程产生的噪声等，多视和滤波都可以达到去噪的效果[2]。 　　 　　图1 卫星雷达数字影像技术实施程序 　　2卫星雷达数字影像在矿区地表沉降监测中的应用 　　卫星遥感影像具有覆盖范围大、更新周期短、费用少及多光谱信息量大等特点，被广泛应用于大范国、多用途、多时域的城市综合调查和研究。由于遥感卫星的轨道高、分辨率较低、光谱波段较多等原因，因此，卫星影像上记录矿区地表信息也有其特点。 　　地下煤炭资源的大规模开采，破坏了上覆岩层的稳定性，从而破碎、垮落，导致地表的下沉、塌陷。这些地表沉降区造成大面积土地资源的破坏，在高潜水位地区甚至出现大量积水，导致区域生态环境的变化，因此，对煤矿沉降区因地制宜进行复垦、治理，对恢复被破坏的土地资源，维护区域的生态环境有重要意义。为此，需要掌握沉降区地随开采活动的扩展规律，沉降区面积、破坏程度，积水区范围和深度，水污染状况等数据资料。采用常规的技术手段，如地形测量是难以胜任的，特别是对于积水沉降区来说更是如此。遥感图像则能快速、周期性地提供沉降区的大量有关信息，为矿区土地复垦工程服务。卫星雷达影像的空间配准是干涉影像处理中关键的一步。由于地形对相位的影响非常敏感，如果两幅图像所对应的同名点有所错开，将会产生很大的测量误差。因此，要求主辅图像空间配准精度一定要达到亚像元级的水平。配准的关键在于计算两幅复图像所对应的同名点的相对偏移量，一旦确定了两幅图像同名点之间的相对偏移量，就可以对一幅图像进行插值重采样，完成两幅图像的配准。根据所选取的两幅卫星雷达影像数据，采用基于几何模型的同名点预测自动配准算法，两幅图像经过该算法处理后，得到匹配点。记下相关系数最大时窗口中心像元的行列偏移量，将其作为控制点进行二阶多项式拟合来确定两幅图像中任何一个匹配位置的相对偏移量，即可以确定几何变换模型，最后根据相对偏移量用双线性插值法对辅图像作插值重采样，从而完成配准操作。配准后控制点残差的方均根误差为0.1～0.2像元且有效控制点一般至少在20以上，基