IBIS-L系统与TM30全站仪在边坡变形监测中的对比分析
- 格式:pdf
- 大小:251.08 KB
- 文档页数:4
第35卷第4期 2015年8月 大地测量与地球动力学 Journal of Geodesy and Geodynamics VoI.35 No.4
Aug.,2015
文章编号:1671—5942(2015)04—0718—04
IBIS—L系统与TM30全站仪在 边坡变形监测中的对比分析
余红楚 左小清 字陈波 徐小坤 王利启 刘海波 1 昆明理工大学国土资源工程学院,昆明市文昌巷68号,650093 2华能澜沧江水电有限公司,昆明市世纪城中路1号,650214
摘要:使用IBIS-L系统对大型水电站边坡长时间的监测数据进行分析和处理,并与TM30全站仪进行变 形监测精度对比。试验表明,IBIS系统可得到在雷达视线方向优于mm级精度的形变结果,优于TM30全站 仪,可以实现高分辨率、高精度、实时的大型水电站边坡长时间变形监测。 关键词:IBIS-I 系统;TM30全站仪;边坡变形监测;对比分析 中图分类号:P258 文献标识码:A
传统水电站边坡变形监测技术主要包括水准 仪、全站仪及测量机器人等。水准仪精度高,但在 陡峭的边坡受到观测高差的限制,难以施测。采 用全站仪和测量机器人需要在边坡上埋设大量观 测墩,人力资源成本高,前期准备周期长,后期需 要对每个观测墩的数据进行计算整理 ]。IBIS—L 系统获取监测区域的二维影像,利用合成孔径技 术和步进频率连续波技术实现雷达影像方位向和 距离向的高分辨率,通过比较不同时间影像中目 标点的电磁波相位信息,采用干涉技术求取监测 区域的变形量,实现优于mm级的微变形监 测¨2 ]。本文对IBIS-L系统为期一个月的大型水 电站边坡监测数据进行分析,并与TM30全站仪 进行变形监测精度对比。 1试验方案与数据处理 1.1试验方案 试验在某水电站边坡进行。选取位于不同边 坡上的3个监测点进行对比分析,其中包含两个 自动观测点和一个人工观测点。 TM30全站仪监测频率为4 d测量一次,监 测周期30 d。利用水电站提供的GPS基准网中 的控制点作为监测基准点。数据采集完成后先进 行仪器常数改正、气象改正和倾斜改正等距离改 正,再计算监测点坐标。 IBIS—L系统安置在边坡对岸的简易棚内。 试验中IBIS—L系统的测点无法直接与边坡已有 的TM30全站仪监测点相重合,所以为了方便对 比分析,在TM30全站仪监测点旁安置角形反射 器,以便准确识别全站仪监测点的位置。IBIS—L 系统24 h实时监测,一次采样时问为6 min。 1.2监测数据处理 1.2.1 TM30全站仪数据处理 TM3O全站仪测量采用08—12、O8—15、O8—21 和08—25观测值的均值作为初值,计算之后每期 的三维位移向量。由于IBIS—L系统测量的是仪 器视线方向位移,必须进行投影变换。采用 TM30全站仪测出IBIS—L系统中心的三维坐标, 计算视线向量。设向量A一(X ,Y ,Z ),B一 (Xz,Y ,Z。),则A向量在B向量上的投影为 X1X2+yly2+Z1Z2 __ 。 据此,计算TM30监测 的三维位移向量在IBIS—L系统视向的投影,并计 算累积位移(表1、2)。 1.2.2 IBIS-L系统监测数据处理 1)数据定标和聚焦处理 IBIS-L系统获得一维信号数据,对数据进行 定标和聚焦处理,得到距离向分辨率0.5 m、方位 向分辨率4.4 mrad的二维扇形图像 。 2)干涉处理 通过设定阈值(强度值或相干系数)去掉质量
收稿日期:2014 12 28 项目来源:糯扎渡电厂边坡多波束雷达观测研究项目(NZDDC2014/P16)。 第一作者简介:余红楚,硕士生,主要研究方向为InSAR数据处理,E—mail:1063716089@qq.corn。 通讯作者:左小清,博士,教授,主要从事数据挖掘和遥感图像处理研究,E—mail:zuoxq@163.COFI1。 第35卷第4期 余红楚等:IBIS-L系统与TM30全站仪在边坡变形监测中的对比分析 719 O8—15 08—17 08—21 08—25 08—28 09一O2 O9一O6 O9—1 0.0O1 3 o.oo2 o o.oo1 2 —0.002 4 o.ooo 2 o.oo2 1 o.oo1 9 0.O03 3 0.O00 4 o.002 8 o.oo1 7 o.oo1 6 0.003 6 o.002 1 0.000 2 o.oo1 4 o.oo2 1 o.002 1 0.0O1 2 o.OOO 7 o.OOO 8 0.000 5 0.OO1 4 0.OO1 6 o.oo1 4 o.oo3 6 o.o0O 8 o.OO1 2 o.OO2 1 o.oo2 o 0.0O1 6 o.oo3 4 o.oo1 2 o.oo1 3 o.oo1 9 o.oo3 2 o.oo2 o 0.001 7 o.oo2 2 o.ooo 5 o.ooo 4 0.00o 5 o.oo1 1 O.0Ol 9 0.003 2 0.001 4 O.O03 8 o.oo3 4 o.oo1 3 0.O00 6 0.002 9 o.OO1 6 0.003 1 0.002 4 0.001 6 o.OO1 5 o.ool 6 346.809 2 637.435 6 457.192 5 346.806 8 637.de35 4 457.19O 4 346.8O1 9 637.428 3 457.19o 8 346.804 5 637.425 3 457.192 4 346.808 1 637.424 4 457.192 2 346.8O9 5 637.426 5 457.194 3 346.81O 7 887.524 5 685.796 3 318.621 O 887.525 9 685.797 9 318.6l9 6 887.529 5 685.798 7 318.619 8 887.527 4 685.796 7 318.621 4 887.530 o 685.802 6 318.623 1 887.528 1 685.802 2 318.616 4 887.529 8 21.926 6 51.977 2 97.584 9 21.927 1 51.976 1 97.583 0 21.934 3 51.974 7 97.586 8 21.940 8 51.976 o 97.586 2 21.933 7 —51.977 6 97.583 1 2l_931 3 51.979 2 97.581 6 21.932 9 o.oo1 1 0.O03 2 0.004 0 637.425 4 685.809 4 —51.975 2 o.002 8 —0.000 8 —0.003 5 457.191 5 318.620 6 —97.585 1 o.002 7 一o.001 5 0.001 1 346.803 7 887.528 3 —21.931 8 o.002 3 一o.002 8 o.oo3 9 637.425 2 685.800 3 —51.971 3 一o.ooo 7 一o.002 1 一o.oo1 7 457.190 8 318.618 5 —97.586 8 不好的点,保证干涉图像的可靠性_5_。对数据进 行干涉处理,首先选择一个主影像,然后分别和乘0 余的影像进行干涉。干涉后,可以设定合适的滤 波窗口对于涉图进行滤波处理,除去噪声的影响, 保证相位的连续性,加快解缠的速度 ]。滤波结 束后可以对滤波后的影像进行相位解缠,获得真 实的距离相位[7]。 3)差分处理 对干涉图进行两两差分处理,获得不同时刻 距离向相位变化值。把差分获得的相位转换为距 离值,得到位移量。 4)误差校正 地面雷达观测的误差主要为视向形变误差, 包括干涉相位误差和系统工作频率偏移,干涉相 位误差包括大气干扰相位误差、干涉相位噪声和 散射体去相干。试验中,监测周期为30 d,系统频 率能维持稳定,系统的频率偏移极其微小,其误差 在0.O1 mm内,可以忽略不计。散射体具有高相 关性,散射体去相干可以忽略 。 1.175 5 o.8o2 2 1.186 7 0.475 2 0.972 2 2.8l3 3 2.714 5 1.735 7 1.654 5 o.876 6 2.546 2 2.090 6 1.784 4 2.405 7 —1.428 5 1.315 o o.812 2 o.340 1 2.O56 4 1.347 o 2.O81 4 2.403 8 0.722 O 1.4=58 2
大气水汽含量在时间空间上的变化引起的雷 达信号延迟是雷达干涉测量中不可忽视的误 差_9]。根据雷达图像的热信噪比、估计信噪比、相 关性值和位移曲线,结合实地情况,选择热信噪比 高、估计信噪比高、相关性值接近1,且位移量在 误差范围内的点作为GCP点(地面控制点)。 GCP点相位变化只受环境变化影响。采用单个 GCP点或多个GCP点加权平均的方式对观测墩 所在的像素点进行校正,从而得到像素点的真实 位移。例如,点位1经过GCP点校正后的位移变 化如图l所示,从图上看,位移变化在1 mm之内 波动。这是因为GCP点和观测点位在雷达图像 上处在不同的像元,GCP点受环境的影响不能完 全反映观测点位微波路径上的情况,校正后存在 一定的误差。 IBIS—L系统是连续观测,TM30全站仪是非 连续的,监测周期为4 d一次,位移量是观测当天
1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 72O 大地测量与地球动力学 2015年8月 。 。 0 萋。
_o -0
●● … I● -
_ ‘ i t ’ i 6 一 I 泣 l
。 I i l ,.. ● 量 !。 t- - ,‘J
, I1 j I c. I
..
霹 『
I== 啊! __ : 国阻l ● 1 : 一 r’ :
l1 I I , I 。: 。 I _ . .j l I一 , 1
1 I 口 ’
-- -_ j 1 , ' I 孵 { 『
图1点位1经GCP点校正后的位移变化图 Fig.1 The displacement change graph of point one after GCP points correction
某一时段或某一时刻的观测值,故选取IBIS—L系 统相应观测日期的平均位移量与TM30全站仪 监测值进行比较。从IBIS—L系统观测数据中提 取出O8—15、08—17、08—21等的位移量,计算它们每 天的平均位移量和累积位移量(表3、4)。 2精度对比分析 根据TM30全站仪和IBIS_L系统在雷达视 线方向监测的累积位移值,绘制图2、3。 由图2看出,TM30全站仪测量每个点的累 积位移都在4-3 mlTl内波动,呈不规则分布。这 是因为水电站边坡都比较稳定,30 d内形变值在 1~2 mlTl左右,而由于观测距离较远,约1 km, TM30全站仪观测时受环境和视差等影响较大, 误差在Film级,无法监测出2~3 ITIITI的位移值。 由图3看出,IBIS—L系统监测误差在_4-1 iTlm 之内波动,但能监测出1Tim级的形变趋势,所以 点2不存在形变,点l在雷达视线方向呈现远离 形变趋势,点3在雷达视线方向呈靠近形变趋势。 点位1和点位3处在不同的边坡,实际位移方向 与雷达视线方向的夹角不同,实际位移在雷达视 线方向的投影有正有负,呈靠近IBIS—I 系统或远 离IBIS—L系统趋势。 采用高次曲线拟合和多项式拟合,对点位1 和点位3的30 d雷达观测数据进行分析,拟合曲 线如图4所示。点位1呈倒指数曲线图形,点位 3呈双曲线图形。