变形观测数据处理.ppt

《变形观测数据处理.ppt》由会员分享，可在线阅读，更多相关《变形观测数据处理.ppt（42页珍藏版）》请在装配图网上搜索。

1、变形观测与数据处理 成都理工大学地球科学学院 测绘工程系 授课 曾涛 变形观测数据处理 1 概述 2 变形观测数据预处理 3 变形分析与预报 第一节 概述 1、数据处理目的： 观测成果计算、分析时，应根据最 小二乘和统计检验原理对控制网和观测 点进行平差计算，对测量点的变形进行 几何分析与必要的物理解释。 第一节 概述 2、数据处理要求： 观测值中不应含有超限误差，观测值中的 系统误差应减弱到最小程度。 合理处理随机误差，正确区分测量误差与 变形信息。 多期观测成果的处理应建立在统一的基准 上。 按网点的不同要求，合理估计观测成果精 度，正确评定成果质量。 第一节 概述 3、数据处理前的准备工

2、作： 核对和复查外业观测成果与起算数据； 进行各项改正计算； 验算各项限差，在确认全部符合规定要求后， 方可进行计算。 第一节 概述 4、数据处理方法： 数据检验：粗差剔出；超限误差检验；稳 定性分析 平差处理：经典平差；自由网平差；秩亏 网平差；拟稳平差等 成果整理：数据整理；绘制过程线；等值 线图；变形值分布剖面图等 变形分析与预报：回归分析法 ；确定函数 模型法 第二节 变形观测数据 处理 变形观测数据预处理的目的、 内容及方法 变形观测数据整理及资料报 告的编写 第四章 变形监测资料的预处理 4.1 监测资料检核的意义与方法 复习要点 监测资料检核的意义 确保观测值的质量 (精度 、

3、可靠性 )， 为变形 分析做准备 。 监测资料检核的方法 野外检核 (限差要求 ) 室内检核 (校核 、 统计分析 、 逻辑分析 ) 大偏差问题 (粗差与突变 ) 第四章 变形监测资料的预处理 4.2 一元线性回归 分析法 进行资料的检核 复习要点 应用 一元线性回归 分析方法 进行 变形监测 资料 空间相关性 检核 的基本思路和过程 4.3 监测网观测资料的数据筛选 数据筛选 的基本 原理 与检验步骤 超限误差局部检验中 ， u检验法 、 检验法 、 t检验法等三种检验方法的本质区别 第四章 变形监测资料的预处理 4.4 监测资料的奇异值检验与插补 复习要点 自动化监测系统中数据的奇异值检验

4、方法 监测资料插补的原因和插补方法 4.5 小波变换用于信噪分离 变形分析中的 为什么要 应用小波变换 ？ 了解 小波变换 的基本概念 小波变换在变形分析中的应用体现 第四章 变形监测资料的预处理 4.6 变形监测成果的整理 复习要点 以 “ 工作基点位移对变形值的影响 ” 为例 ， 说 明变形监测成果整理的重要性 观测资料的整编 变形过程线的定义与的绘制 建筑物变形分布图 变形值的统计规律及其成因分析 以大坝为例 ， 说明引起变形的原因 第四章 变形监测资料的预处理 4.7 监测资料 的 管理 复习要点 人工管理处理 计算机辅助人工处理 数据库管理系统 监测资料的管理 监测资料的处理 监测资

5、料的解释 一、变形观测数据预处理的目的 目的：提前发现错误，剔除粗差及超 限误差，保证原始数据的正确性，为平 差处理及数据分析提供保证。 二、变形观测数据预处理的内容 观测值中的超限误差，除在观测过程 中应严格作业、认真检核随时予以排除外， 在变形分析中，还应通过检验将判定含有粗 差的观测值予以剔除。 观测值中的系统误差，除在作业中应 严格进行仪器检校、按规定观测程序操作予 以减弱或对受大气及其他影响的数据作预处 理外，对于高精度监测网还应进行系统误差 的统计检验和补偿。 三、变形观测数据预处理的方法 对于多次重复观测列中的粗差检验，可采用格 拉布斯准则或狄克逊准则；（ 3倍中误差法则） 对高

6、精度监测网的粗差检验，可采用巴尔达数 据探测法或 S稳健估计法。 稳定点的检验可采用统计检验方法。先作整体 检验，在判别有动点后再作局部检验，找出变 动点予以剔除，最后确定出稳定点组。亦可采 用按单点高程、坐标变差和观测量变差的 u、 x2、 t、 F检验法，或采用按两期平差值之差与 测量限差之比的组合排列检验法。 四、变形观测数据预处理的算例 水准点稳定性检验： 平均间隙法 超限误差检验： 、 B检验法 五、变形观测资料整理 1. 概述 变形观测除了现场进行观测取得第一手资料外，还必 须进行观测资料的整理分析。 观测资料整理分析主要包括两个方面的内容： 观测资料的整理和整编 。这一阶段的主要

7、工作是对 现场观测所取得的资料加以整理，编制成图表和说明，使 它成为便于使用的成果； 观测资料的分析。这一阶段是分析归纳建筑变形过 程、变形规律、变形幅度，分析变形的原因，变形值与引 起变形因素之间的关系，找出它们之间的函数关系，进而 判断建筑物的工作情况是否正常。在积累了大量观测数据 后，又可以进一步找出建筑物变形的内在原因和规律，从 而修正设计的理论以及所采用的经验系数。 2. 资料整理步骤： 资料的校核 变形观测资料的插补 填表 绘图 资料的校核 校核各项原始记录，检查各次变形观测值的计算有 否错误。 （ 1）原始观测记录应填写齐全，字迹清楚，不得涂 改、擦改和转抄；凡划改的数字和超限划

8、去的成果，均 应注明原因，并注明重测结果的所在页数； （ 2）平差计算成果、图表及各种检验、分析资料， 应完整、清晰、无误； （ 3）使用的图式、符号，应统一规格，描绘工整， 注记清楚； （ 4）观测成果计算和分析中的数字取位要求见规范。 变形观测资料的插补 当由于各种主、客观条件的限制，实测资料出现漏 测时，或在数据处理时需要利用等间隔观测值时，则可 利用已有的相邻测次或相邻测点的可靠资料进行插补工 作。 （ 1）按内在物理联系进行插补 按照物理意义，根据对已测资料的逻辑分析，找出主 要原因量之间的函数关系，再利用这种关系，将缺漏值 插补出来。 （ 2）按数学方法进行插补 线性内插法：由某两

9、个实测值内插此两值之间的 观测值时，可用： 为效应量； 为时间 拉格朗日内插计算：对变化情况复杂的效应 量，可按下式计算： 式中： 为效应量； 为自变量。 用多项式进行曲线拟合，即： 式中方次和拟合所用点数必须根据实际情况适 当选择。 填表 对各种变形值按时间逐点填写观测数值表。 例如，某大坝变形观测点，根据观测记录整理填 写的位移数值表见表 1和表 2。 表 1 某大坝 5观测点 1967年累计位移数值表（单位： mm） 观测 日期 1 月 10 日 2 月 11 日 3 月 10 日 4 月 11 日 5 月 10 日 6 月 10 日 7 月 11 日 8 月 11 日 9 月 10 日

10、 10 月 11 日 11 月 11 日 12 月 10 日 4.0 6.2 6.5 4.2 4.3 5.0 2.2 3.8 1.5 2.0 3.5 4.0 注： 为累计位移值。 表 2 沉降数据一览表 h1 本期 (mm) 累计 (mm) h2 本期 (mm) 累计 (mm) 基期测量 500165.30 基期测量 500129.80 第一期 500164.99 0.31 0.3100 第一期 500129.49 0.31 0.3100 第二期 500165.17 -0.18 0.1300 第二期 500130.10 -0.61 -0.3000 第三期 500164.86 0.31 0.44

11、00 第三期 500129.40 0.70 0.4000 第四期 500164.28 0.58 1.0200 第四期 500128.39 1.01 1.4100 第五期 500163.88 0.40 1.4200 第五期 500128.70 -0.31 1.1000 第六期 500164.28 -0.40 1.0200 第六期 500129.01 -0.31 0.7900 h3 本期 (mm) 累计 (mm) h4 本期 (mm) 累计 (mm) 基期测量 500209.30 基期测量 500216.30 第一期 500209.03 0.27 0.2700 第一期 500215.99 0.31

12、 0.3100 第二期 500209.43 -0.40 -0.1300 第二期 500216.48 -0.49 -0.1800 第三期 500208.91 0.52 0.3900 第三期 500215.78 0.70 0.5200 第四期 500207.81 1.10 1.4900 第四期 500214.68 1.10 1.6200 第五期 500207.90 -0.09 1.4000 第五期 500214.89 -0.21 1.4100 第六期 500208.02 -0.12 1.2800 第六期 500214.89 0.00 1.4100 h5 本期 (mm) 累计 (mm) h6 本期

13、(mm) 累计 (mm) 基期测量 500206.30 基期测量 500209.50 第一期 500206.09 0.21 0.2100 第一期 500209.29 0.21 0.2100 第二期 500206.52 -0.43 -0.2200 第二期 500209.81 -0.52 -0.3100 第三期 500205.82 0.70 0.4800 第三期 500208.89 0.92 0.6100 第四期 500204.72 1.10 1.5800 第四期 500208.10 0.79 1.4000 第五期 500204.81 -0.09 1.4900 第五期 500205.99 2.11

14、 3.5100 第六期 500204.90 -0.09 1.4000 第六期 500208.10 -2.11 1.4000 绘图 绘制各种变形过程线，建筑物变形分布图等。 观测点变形过程线可明显地反映出变形的趋势、 规律和幅度，对于初步判断建筑物的工作情况是否正 常是非常有用的。 图 1是根据表 1绘制的某大坝 5观测点的位移过程 线。图中横坐标表示时间，纵坐标为观测点的累计位 移值。 图 1 某大坝 5观测点的位移过程线 图 2 珠江骏景沉降观测 T P曲线图 建筑物等沉降曲线图 回弹量纵、横断面图 地基土深层侧向位移图 滑坡观测点的位移与沉降综合曲线图 第三节变形分析与预报 变形分析与预报

15、的内容 变形分析与预报的方法 几何分析在于确定变形量的大小、方向及其变化 , 即 变形体形态的动态变化。物理解释在于确定引起变形的原 因 (例如是由某种荷载为主引起的周期性变形 )和确定变形 的模式 (属于弹性变形还是塑性变形 , 是 相对或绝对变形 量的大小、几何分布和变化规律。自身内部变形还是整体 变形等 )。一般来说 , 几何分析是基础 , 主要是确定相对 和绝对位移量 , 物理解释则是从本质上 认识变形。变形 的物理解释和变形预报可根据确定函数法如 动力学方程 进行 , 也可根据大量的监测资料用统计分析法进 行 , 同 时还可将两种方法结合起来进行综合分析预报。对于 为 基准进行约束平

16、差为宜。周期观测点场稳定性的统计检验 一 、 变形分析与预报的内容 二 、 变形分析与预报的方法 回归分析法 时间序列分析法 频谱分析法 卡尔曼滤波法 有限元法 人工神经网络法 小波分析法 系统论方法 时间序列是按时间顺序的一组数字序列。 一 组观测 值，若沿 着时间先后有顺序产生， 则称此组观测值为一时间 序列，而正整 数 N被 称 为时间 序列的 长度。 时间序列分析就是利用这组数列，应用数理统计 方法加以处理，以预测未来事物的发展。 时间序列分析是定量预测方法之一。 该方法方法简单易行，便于掌握，但准确性差， 一般只适用于短期预测。 时间序列分析法简介 它的基本原理： 一是承认事物发展的

17、延续性。应用过去数 据，就能推测事物的发展趋势。 二是考虑到事物发展的随机性。任何事物发 展都可能受偶然因素影响，为此要利用统计分析 中加权平均法对历史数据进行处理。 时间序列分析法的基本原理 时间序列分析法的基本计算流程 以平稳时间序列为例 平稳随机序列或平稳随机过程 . 它具有如下 两个特点 : 第一 ,它的数学期望、方差不随时间变化 ,亦即 E( Xt) = ( 1) E( Xt - ) 2 = 2 (2) 第二 ,协方差函数只与时间间隔 t-s 有关： Cov(Xt,Xs) = E(Xt-)( Xs -)= R(t - s) (3) 对某一随机序列的进行样本平均值、方差、自协方 差、自

18、相关函数 ,可作如下的一个估计： 如果平稳序列 Zt 的数学期望为 0 ,自协方差函数有 如下性质 Rz () = Cov( zt , zt +) = z2 = 0 或 = 0 0 则称 Zt 为白噪声序列或纯随机序列 假定 zt 为白噪声 ,且序列 Xt 满足 则称 Xt 为 P 阶自回归序列 ,这种模型简记为 AR(p) ,在 上式中参数 i i = 1 ,2 , P 称为自回归系数。 对上式两边乘因子 Xt- 再取数学期望后除以方差 2x , 得： 确定 AR( P) 模型的阶数 ,即选取 AR( P) 的最佳阶数 P： 选择 P 阶数的标准自然是预报时的残差方差最小为原则 P阶数标准通

19、常使用最小信息准则 经验表明 ,对样本不大的长期预报来说 ,按 r() 的绝对 值选择，一般选择绝对值较大 3个前期变量 算例 某站以 1980-1994年平均位移量 xt ( t = 1 ,2 , ,15) (表 1) , 作为各态历经的平稳时间序列的一个现实 ,试建立 位移的 Yule- walker 自回归方程组 ,并对第 16 年 (即 1995年 )的位移量进行预报。 (7) 回报检验 ,将历史个例因子代入上式求出历年的 。例如 ,计算 (n=9 ,即第 9年位移量距平的计算值 ) 应将 xd(9 - 4) = xd5 = 9 , xd(9 - 7) = xd2 = - 10 , xd(9 - 8) = xd1 = - 5 代入自回归方程 ,即 xd9 = - 0.38 9+0.34 (-10)+0.18 (-5)=-7.72 - 8 类似地算出 xd10 , xd11 , , xd15计算值与实测值列表 3。 (8) 对第 16（ 1995）年位移量进行预报 于是第 16年位移量计算值为 : 即 1995年位移量的预报值为 505.7mm.