Hermert 方差分量估计在混合水准网平差中的应用.doc

Hermert 方差分量估计在混合水准网平差中的应用1 赫尔默特方差分量估计的计算步骤赫尔默特方差分量估计的迭代计算步骤如下：（1） 将观测值按等级或不同观测来源分类，并进行验前权估计，即确定各类观测值的权的初值P1,P2,Pm；（2） 进行第一次平差，求得各类观测值的残差平方和。（3） 按（40）式进行第一次方差分量估计，求得各类观测值单位权方差的第一次估值 ，再依下式定权：式中C为任一常数，一般是选 中的某一个值。（4） 反复进行第二项和第三项，即进行：平差方差分量估计定权后再平差，直至为止，或通过必要的检验认为各类单位权方差之比等于1为止3。2 赫尔默特方差分量估计流程图下图是几何水准和三角高程混合网中的赫尔默特方差分量估计流程图输入水准点个数LM，已知点个数N，几何水准观测值个数n1,三角高程观测值个数n2，总观测值个数n=n1+n2输入已知点编号及其高程输入高差观测值及水准点间的距离输入方程的路径矩阵求解水准网的平差方程的常数项和系数矩阵对观测值定初权根据方程解出各观测值改正数根据观测值改正数及其对应的系数矩阵解出S，WV求出各类观测值的单位权方差是否等于算出各水准点的高程平差值对水准网进行精度评定结果输出是否对观测值再定权3 赫尔默特方差分量估计实例分析3-1 苏通长江公路大桥高程控制网概述一 测区概况 苏通长江公路大桥位于江苏省东南部、长江三角洲东部，北接南通市，南连苏州市，东距长江入海口约100公里，西距江阴长江公路大桥约90公里，地理位置为东经120o59，北纬31o47。桥位北岸在南通市境内的国营南通农场23大队，桥位南端在常熟市东北部的徐六泾，西距常熟发电厂约1公里。该桥是连接长江口段的重要公路通道组成部分，是我国目前跨越长江水面最宽的长江第一特大型桥梁。桥位处于冲积平原的新三角洲，地势低平，海拔高程公23米。桥位区地颀覆盖层两岸江堤之间约6.1公里。江堤高45米，宽46米。两岸江堤均可通大车和汽车。桥位西约6公里是通常汽渡。桥位区外围四周市、集镇均有公路通至桥位附近，交通十分方便。桥位区北岸居民点较少，且住房简陋，南岸村落稠密，大多民房为23层楼房。桥位区属亚热带湿润季风气候区。气候较温和，雨水充沛。受季风影响，四季分明年平均气温15.2oC，最热月平均2830oC（七月份），最冷月平均-0.23.1oC（一月份）。年高温35oC仅3.44.6天。平均降水量为1082.6mm，年最多降水日数为143天（8-9月雨季）。桥位区属长江的感潮河段，每年5月至10月为汛期。11月至次年3月份为枯水期。1月或2月水位最低。洪峰多出现在6至8月。二 测区已有高程控制点资料本项目使用的高程控制网的起算点为：长江北岸南通市境内狼山的基岩水准点II狼山乙，该点为岩层基本标石、岩层混凝土铜标志，位于南通市狼山公园，是长江流域管理办公室于1963年11月埋设。长江南岸常熟市境内虞山公园内，是江苏省工程测绘院于1998年5月埋设。三 高程控制网建立的基本原则为满足大桥在勘测、施工期的施工放样要求，特在大桥两端各约4公里范围内布设了专用的高程控制网，以国家二等水准施测；在顾及大桥施工周期长，对高程控制点的稳定性要求较高的特点，分别在大桥的两端附近各布设三个深埋水准点，其中北岸布设埋深68.5米水准点一座、埋深20米水准点两座，南岸布设埋深50米水准点一座、埋深20米水准点两座；为保障起算数据的可靠性及对深埋水准点稳定性的检测，特在北岸狼山基岩点至埋深68.5米水准点和南岸虞山基岩点至埋深50米水准点间施测国家一等水准；为保证大桥两岸施工在高程方面顺利衔接，使南、北岸高程控制网形成整体布局，特在大桥两岸高程控制网间施测测距三角高程控制网连接，从而使用苏通长江公路大桥桥位高程控制网布设成为一个统一的、整体的专用高程控制网。四 高程控制网的建立苏通长江公路大桥平面控制网由18个点组成，网形为以桥轴线为公共边的两个大地四边形，两边各扩展一个大地四边形，同时在桥轴线的延长线布设四个大地四边形，各平面控制点保证有两个以上的通视方向。平面控制网的所有点位均布设为高程控制网点。为满足对本高程控制网的稳定性要求，保障设计、施工放样的顺利进行，在上述基础上再布设6个钢管深井点。跨江高程传递是本项目中的重点，因此，专为跨江测距三角高程测量建造4个观测墩。由此，苏通大桥高程控制网共有28个高程点组成。其中：一等水准观测共计49.73公里。其中南岸虞山基岩至测区深埋水准点为30.64公里；北岸狼山基岩点至测区深埋水准点为19.09公里。二等水准观测共计31.03公里。其中南岸施测两个环线合计13.49公里；北岸施测两个环线合计17.54公里。过江水准测量用三角测距高程测量进行。本次平差计算中假定以一等水准观测从南岸虞山基岩测得深埋水准点28号高程已知。其它各27个点为待定点，各个观测值的权为各自观测点区的距离的倒数。苏通长江公路大桥高程控制网图见附图（一）。3-2 赫尔默特方差分量估计计算一 苏通长江公路大桥已知数据和观测数据表（一）为已知水准点高程表，为保证整个高程控制网的整体性，此次平差计算只用了一个已知点，即苏通长江公路大桥南岸由虞山基准点经一等水准测量引测到的28号点为已知点。表（一）已知数据表点号高程(m)283.5897表（二）为高差观测数据表，共有34个观测值，其中前面30个观测值为二等几何水准测量得到的观测值，最后4个观测值为跨江测距三角高程测量得到的观测值。表（二）观测数据表编号始点终点高差观测值(m)距离(km)1170.97473.10027810.21070.2003890.10290.1004910-9.48180.7005103-1.49430.7006321.56230.400721-1.87463.7008240.17401.200945-0.09391.3001056-0.96320.400116133.30651.700121312-0.35461.400131211-1.75491.300141110-0.38211.3001528270.91600.800162726-0.14990.1001726171.57081.3001817167.62260.200191615-0.91830.100201514-6.76450.900211428-2.27672.300222625-0.43440.6002325240.45850.7002424231.34280.500252322-2.25040.3002622211.59430.3002721200.37092.4002820190.56050.700291918-0.14441.2003018170.07311.200318160.67445.80032815-0.22885.800339160.57125.80034915-0.33985.7003.3、 Helmert具体实施步骤首先采用常规水准网平差的定权方法（）进行定权平差结果就是（表1第一行）然后采用helmert方差分量估计的方法将水准路线分为陆上、跨江两类进行处理。经过五次迭代平差。各类单位权方差和之比接近于1：1如表1所示表1 各次迭代的比较结果迭代次数系数矩阵:第1次4.376015.73760.52626.69570.079第2次0.28832.00030.05050.69020.073第3次0.21970.15540.04390.05190.846第4次0.21960.13230.04390.04410.995第5次0.21960.13090.04390.04371.005Helmert方差分量估计前后的精度分别为1.58mm和0.59mm，其对各点高程的影响见表二表2 各点高程的影响列表（单位：mm）点 号 12345678原始高程值(m)1.68403.560021.997673.734023.639132.675912.658685.8517高程值(Helmert)(m)1.68403.55861.99623.73263.63872.67552.65875.8517点 号910111213141516原始高程值(m)12.8693812.972225.747013.49093.872965.627875.98244 5.87664 高程值(Helmert)(m)12.869412.97235.74713.4905 3.87265.62755.98205.8664点 号1718192021222324原始高程值(m)7.8561612.639113.543898.528985.9259 5.853165.998145.43764 高程值(Helmert)(m)7.852912.630913.5493 8.5344 5.92665.8535 5.9979 5.4374点 号2526272829303132原始高程值(m)5.065593.471325.72174 4.378953.931724.366064.515975.59998 高程值(Helmert)(m)5.0665 3.47235.7227 4.3799 3.92154.35584.50573.5