区域沉降观测稳健迭代计算与点位稳定性分析-精.docVIP

区域沉降观测稳健迭代计算与点位稳定性分析 赵显富(唐山大学　063000) 　　【摘　要】　本文介绍变形监测平差方法及应用环境，提出一种新的定权方法用于区域沉降观测中基准网点稳健迭代平差，计算合理的位移量，分析点位稳定性。测绘信息网　　一、变形监测平差方法及其应用环境　　平差方法［1］总体上分为外坐标平差法与内坐标平差法两种。变形监测网数据处理多采用内坐标平差法，亦即自由网平差。依据网中相对稳定点多少的不同，自由网平差分为经典、普通秩亏、拟稳秩亏及加权秩亏等。这些平差方法分别属于固定基准、重心基准、拟稳基准和加权重心基准平差。它们在最小二乘、最小范数条件约束下，可以用统一的数学模型来表达，即 　　　　　　　　　　　　　　(1) 　　其中，第1式为误差方程，第2式为最小二乘条件，第3式为最小范数条件。GT为基准条件方程系数，Px为未知数的权，它代表点位的稳定程度。对式(1)求解得　　未知数： 　　　　　　　　　　　　　(2) 　　未知数协因数：测绘信息网 　　　　　　　　　　　　　　(3) 　　上述模型对于同一监测网，通过对PX的取值不同，即可实现前述的各种平差。当知道网中一点绝对稳定时，取该点PjX=1，其余各点PjX=0，这时即为经典自由网平差；当网中各点无先验条件，即各点具有同等稳定程度时，取各点的PX均等于1，这时平差结果即为普通秩亏自由网平差；当网中一部分点相对另一部分点稳定时，取相对稳定点组的PjX=1，其余点组PjX=0，平差结果即为拟稳秩亏自由网平差；当网中各点没有明显的稳定与非稳定界限，即各点均有一定的稳定性指标时，取PX＝［0，1］区间内的数值，其平差结果即为加权秩亏自由网平差［2］。　　由此可见，代表点位稳定程度的权PX，在平差中起着重要的作用，它决定了平差基准，从而直接影响平差结果及点位位移量的计算。因此，在变形监测数据处理时，应根据项目的具体情况，合理确定PX，从而选取适当的平差模型。PX一般是难以合理确定的，而稳健估计［3］原理正是利用点位权PX与位移值d之间的函数关系，通过加权秩亏平差迭代计算出较为切合实际的与PX。测绘信息网　　二、稳健迭代平差计算模型与迭代权的确定　　1. 稳健迭代平差计算模型　　对于同一沉降监测网，当网形及观测纲要两周期不变时，利用式(2)计算结果只随L而变。设两周期高程未知数平差值为，误差方程常数为LⅠ、LⅡ，则有位移量平差值，即 　　　　　　　　　(4) 设 　　　　　　　　　　　　　　(5) 　　则式(4)简写为测绘信息网 　　　　　　　　　　　　　　　(6) 　　由式(6)可以看出，为了求得与实际尽可能吻合的位移量，其关键是PX的选取，即选择合理的基准。而PX是与点位的稳定性(位移量)有关的，即PX是的函数，这一函数总体成反比关系，在文献［3］中给出了几种迭代权函数。而问题最终解决办法只能用迭代计算法，即给出一个初始值（0）（可用普通秩亏自由网平差得到），然后利用（0）计算迭代权，再利用式(6)计算位移量，直到迭代计算收敛为止。　　2. 区域沉降迭代权的确定　　若监测网点的稳定性具有一定的先验信息，则利用文献［3］中迭代权函数进行迭代，均能得到较为满意的结果。但当无任何先验信息时，其迭代结果与实际出入较大。例如图1所示的沉降观测网，假设实际上有两点各沉降4 mm，其余两点稳定不动，利用式(6)进行普通秩亏平差时，为了满足基准条件Σ=0(最小范数条件)，其平差结果必然得到两非稳定点各下沉2 mm，两稳定点各上升2 mm的错误结论，使稳定与非稳定点无法区分。利用这一初始结果，无论用文献［3］中哪种迭代权函数进行迭代，均无法得到正确位移量。 图 1 　　尽管无法知道沉降量大小，但利用下述定权方法进行稳健迭代可以得到切合实际的相对沉降量，找出网中相对最为稳定的点组。　　设首次利用普通秩亏平差得两周期高程位移值为（0），其中为伪上升点组，为伪下沉点组。分如下三种情况确定迭代模型。　　(1) 完全沉降区域测绘信息网　　对于区域沉降观测而言，该地区处于下沉状态是可知的［4］(例如受荷载而下沉、港区沉降等观测)。设伪上升点组加权均值位移量为 　　　　　　　　　　　　　　　　(7) 剔除伪上升量后的位移量为 　　　　　　　　　　　　　　　　(8) 下一次迭代权取测绘信息网 　　　　　　　　　　　　　　　　　(9) 　　利用式(6)得加权秩亏平差解为（1），重复式(7)～(9)，进行反复迭代定权、平差。式(7)～(9)的通式为 　　　　　　　　　　　　　　　(10) 　　当两次迭代结果在一定精度内，即　　当|(k+1)-（k）|≤ε　(一般取ε为0.001～0.01 mm)时，迭代结束。通过多个模拟算例［5］，这一迭代定权法均能取得与实际相符合的沉