第26卷第4期2001年8月武汉大学学报・信息科学版
G eomatics and Information Science of Wuhan University Vol. 26No. 4Aug. 2001
文章编号:1000-050X(2001) 04-0354-07文献标识码:A
一种基于可靠性的工程控制网优化设计新方法
张正禄1 罗年学1 黄全义1 梅文胜1
(1 武汉大学测绘科学与技术学院, 武汉市珞喻路129号,430079)
摘 要:在论述工程控制网优化设计方法及网的精度、可靠性、灵敏度以及费用等准则的基础上, 提出了一种基于观测值内部可靠性指标的工程控制网模拟法优化设计的新思想和算法, 指出了该方法的优点和特点; 介绍了该优化设计方法所采用的“科傻”软件, 并用实例说明了用该方法进行工程控制网优化设计的计算步骤、优化效益以及优化设计的必要性。
关键词:工程控制网; 优化设计; 可靠性; 多余观测分量中图法分类号:P221. 2; P258
按用途划分, 施工控制网, () 和约束网() , 独立网只固定一个点和一个方向, 即已知一点的坐标和一条边的方位角, 网的尺度由测边确定; 约束网有两个或两个以上已知点。一般工程控制网以采用独立网为宜, 已知点和已知方位角可以在假定坐标系中给定, 也可以采用国家或地方坐标系的一个已知点和一个已知方位角。
网的优化设计有解析法和模拟法两种。解析法是以最优化理论为基础的严密方法, 其数学模型一般表示为:
f (X ) →min φ(X ) ≥0(1) ψ(X ) =0式中, 第一式为目标函数; 第二、三式为约束条件。优化设计的实质是通过给定的约束条件求目标函数的极值而得到最优解。目标函数可以是精度、可靠性、灵敏度或费用等指标, 约束条件也可以是上述指标。有人总以为解析法的结果是最优的, 实际上并非如此。因为目标函数和约束条件都是一种经过简化的数学模型, 不能完全反映客观实际。最有代表性的例子是以费用为目标函数或约束条件时, 都是以观测值的权的总和来代表建网费用, 实际上这是很粗糙的。而且, 绝大多数情况
收稿日期:2001-02-23。
项目来源:国家测绘局教学研究基金资助项目(97GH01) 。
, 优化后还存在取整的问题, 同一测站上还要划为便于操
作的等权观测问题, 经过取整和调整后的观测方案不再是最优的了。解析法不但存在上述缺点, 同时还因优化任务不同, 目标函数和约束条件的数学模型也不尽相同, 故研制通用的解析法优化设计软件的难度较大。尽管解析法优化设计方面的研究已很多, 却难以推广, 在实际生产中也鲜有应用。
模拟法实际上是一种试算法。一般作法是:根据优化的任务和设计者的知识和经验, 制定初始设计方案, 用模拟的观测值对该方案作平差计算, 对平差结果作评价, 对初始方案进行修改, 再计算、再修改, 如此多次重复, 直至认为满意为止。模拟法需要一个好的平差程序, 并且该程序还能作观测值模拟, 显示网图及各种精度指标和可靠性数值指标, 便于人机对话式操作, 算法好(如采用序贯算法) , 速度快。但模拟法即使有一个好的优化设计程序, 仍依赖于设计者的知识和经验, 同样的要求, 甚至同样的初始设计方案, 不同人设计出的最终方案也有差别, 甚至差别还很大。由于解析法和模拟法的缺陷, 造成在实际工作中人们很少对工程控制网作优化设计, 大多根据测量负责人的知识、经验和习惯进行布网和观测。对于高精度施工控制网或变形监测网, 往往精度和费用偏高; 而对于测图控制网, 大多用多级
第4期 张正禄等:一种基于可靠性的工程控制网优化设计新方法 355
附合导线进行加密。这些作法虽然都符合规范要求, 但谈不上网的优化设计, 有不少弊病。
针对上述情况, 本文提出了一种基于观测值内部可靠性指标的工程控制网优化设计方法。该方法本质上属于模拟法, 但加进可靠性指标后, 具有量化的优化设计准则, 不以人的知识和经验为转移, 优化结果既具有一致性, 也不失严密性。该方法使用笔者研制的现代测量控制网数据处理通用软件包———科傻系统系列软件之二———Cosa2(或称科达普施Codaps ) 。
低, 当r i 等于1时, 则该观测值完全成为多余, 即使删除之网平差结果也不变。
2) r i 愈小, 该观测值的粗差愈难被发现, r i =0时, 即使有大的粗差或错误也无法发现。根据内部可靠性与外部可靠性具有一致性的性质[3], 平差结果受粗差的影响随r i 的减小而增大。r i 愈大, 则较小的粗差也能发现, r i =1时, 观测值的粗差能完全确定, 粗差对平差结果的影响随r i 的增大而减小。
3) 对于一个确定的网和设计方案, 即在网形和网的观测值总数n (或多余观测数r ) 确定的情况下, 观测值之间的精度相差愈大, 则内部可靠性r i 的值相差也愈大。观测值的精度愈高, 则相应的r i 愈小; 观测值的精度愈低, 其r i 愈大, 即观。
4) () , 多余观测r i 也愈大, 反之, 若r i 大, 则
r 也必定大(建网费用将随r 的增大而增加) 。
1 观测值的内部可靠性及其性质
网的可靠性可分为网的总体可靠性、观测值的内部可靠性(也称局部可靠性) 和外部可靠性。在此, 主要讨论观测值的内部可靠性。
对于一个工程控制网来说, 2-1
) l 差模型(l , A x , σ0P 性r i 为[1]:
=(V V P )
ii
(2)
5) 对于独立网来说, 与边长和角度是不变量
且满足
n
∑r
1
一样, 观测值的内部可靠性r i 是与基准的位置无
关的不变量。
i
=r =n -t
r i
(3) (4)
ω0/ 0l i =σi ・
2 几个有关问题
2. 1 多余观测数和平均多余观测分量
式中, ω0为非中心参数, 对于单个观测值粗差而
言, 其取值与显著水平α和检验功效γ有关, 常取2. 79(α=0. 05, γ=0. 80) 或4. 13(α=0. 001, γ=0. 80) [2]; r 为网的多余观测数, r i 为观测值l i 的多余观测分量, r i 可以反映控制网发现观测值l i (中误差为σ通过统计检验, r i i ) 中粗差的能力。愈大, l i 中能被发现粗差的下界值 l i 愈小; 或对同一个粗差, 检验功率愈大。因此, r i 被定义为观测值的内部可靠性。
观测值相互独立(本文只讨论这种情况) 时, 有:
2
(5) r i =1-σ^2i /σi 式中, σ^i 为l i 的平差值的中误差。若观测值l i 的
精度很高, 则平差后的精度提高很小, 有σ^i ≈σi , 这时r i →0。若无多余观测, 有σ^i =σi , r i =0; 反之, 若l i 的精度很低, 则平差后精度将显著提高, 即σ^i
1) 0≤r i ≤1。r i 愈小, 该观测值在网中的地
设控制网的网点数为m , 已知点数为m k , 未
知点数为m u , 已知边和已知方位角数分别为k s 和k a , 进行有方向观测的设站数为m l , 只进行有边长观测的设站数为m s , 其中有方向或边长观测的已知点设站数为m 1, 总的观测值数为n (其中, 方向观测值数为n l , 边长观测值数为n s ) , 必要观测值个数为t , 多余观测数为r , 平均多余观测分量为 r i , 其中, 方向和边长的平均多余观测分量分别为 r li 、 r si , 方向、边长的最大最小多余观测分量分别为r l max 、r l min 、r s max 、r s min 。某观测值的多余观测分量用r i 、r li 、r si 表示。对于按方向的间接观测平差模型, 网的多余观测数r 和平均多余观测分量可按下式计算:
r =n -t =n -(3m l +2m s -2m 1-k s -k a )
(6) (7)
r =
n
位愈高, 若r i 等于零, 则该观测值不可缺少, 否则将产生形亏。r i 愈大, 该观测值在网中的地位愈
例如, 对于附合导线, 有m 1=2, m l =m u +2, 方向观测值数n l =2(m u +2) , 边长观测值数n s =m u +1, 未知数个数t =3(m u +2) -2×2, 故r =n l +n s
武汉大学学报・信息科学版 2001年356
-t =2m u +4+m u +1-(3m u +6-4) =3。
表1列出了附合导线含不同未知点数时观测值个数和平均多余观测分量。当未知点多于4个
时, r i 小于0. 2。从可靠性指标来看, 附合导线不能作为施工控制网和变形监测网。 在优化设计时, 为了使网有足够的可靠性, 在
表1 附合导线的观测值个数和平均多余观测分量
Tab. 1 Observation Number and Average Ebservation Redundant for Traverses 未知点数
n
180. 375
2110. 27
3140. 21
4170. 18
5200. 15
6230. 13
7260. 115
8290. 10
9320. 09
r i
此提出平均多余观测分量设计值的概念, 用 r i 表示。对于不同用途、不同要求的网, 可以取 r i 的00设计值 r i 为0. 3～0. 6。由 r i 值, 可按下式计算
验值。若要求提高, 还可取小一些, 如k =1～
2. 0。
2. 3 网的质量准则问题
观测值数的设计值n 0, 有:
n
=
1- r i
(8)
令
n i 0=1- r i 0
t =n i
网的质量准则主要是精度、可靠性、建网费
用, 对于变形监测网还包括灵敏度和可区分性。
、相对点位精(或主分量) , 点位精度与基准的位置有关, 对于独立网来说, 最好是将最靠近网的重心的点作为已知点, 以通过该点的最接近中心线的方向作为起始方向。这样可保证点位精度在数值上达到最小。实际上, 应将相对点位精度或最弱边精度作为精度准则, 因为它们是与基准的位置无关的不变量。在模拟法优化设计中, 应取先验单位权中误差计算各精度指标。
网的可靠性指标建议采用平均多余观测分量 r i 衡量, 而且 r si 和 r li 相差不应太大。对于相同的观测方案, r si 和 r li 随边角精度的改变而改变, 例如若提高方向观测的精度, 则 r li 将减小。同时, (r l max -r l min ) 与(r s max -r s min ) 也不能太大。
由§1中观测值内部可靠性的性质4) 可知, 由于r i 的总和反映了多余观测数, 因此, 也在很大程度反映了建网费用, 即r 愈大, 建网费用愈高。
网的灵敏度表现为网点在特定方向上的精度。只有在特殊情况下才需要作灵敏度计算。一般可以采用误差椭圆近似地进行评价, 即网点在要求的方向上误差应较小, 如尽量靠近误差椭圆的短轴方向。
n i 0为观测值的平均必要观测分量, :2. 2 , 纯测角网已逐渐消亡, 纯测边网也很少采用, 广泛应用的边角网或导线网存在边角精度的匹配问题。众所周知, 测边引起(或控制) 纵向误差, 方向引起(或控制) 横向误差, 它们都与边长密切相关。设方向中误差为m ″r , 测边的固定误差和比例误差分别为a 和b , 边长为s , 则由方向中误差引起的横向误差为:
m (9) m q =・S
ρ″
由边长中误差引起的纵向误差为:
m L =
22
a +(bs ) 或 m L =a +bs
(10)
所谓边角精度完全匹配, 是指应满足m q =m L 。由于网的边长变化和仪器的限制, 边角精度匹配是相对的, 不匹配是绝对的。一般认为, 当
m q 和m L 满足下述关系
k
m L ≤m q ≤km L , (k =1～2. 5) (11)
时, 边角精度是基本匹配的。因为观测值的精度与其可靠性成反比, 边角精度是否匹配也可通过可靠性来判别。对于初始方案即边角全测的情况(这时n l ≈2n s ) , 通过计算方向和边长的多余观测数和平均多余观测分量r l 、r s 、 r li 、 r si , 可根据 r li 、 r si 的大小来判断边角精度是否匹配, 一般应满足 r li ≤ r si ≤k r li , k 一般在1～2. 5之间。这
k
里的k 和式(11) 中的k 值无直接联系, 且都是经
3 优化设计思想、计算步骤和特点
观测值的内部可靠性与观测值的精度有密切关系, 而观测值的精度又与建网费用有关。变形监测网的灵敏度实际是网点在特定方向上的精度, 它也取决于网的观测方案设计和观测值的精
第4期 张正禄等:一种基于可靠性的工程控制网优化设计新方法 357
度, 此外, 变形与粗差的可区分性也必然涉及到观测值的精度。因此, 观测值的内部可靠性与观测值的精度、建网费用、监测网的灵敏度和可区分性存在密切的关系, 有的关系可以简化为如式(3) 、(4) 、(5) 或文献[4]所描述的数学公式, 有的关系却难以进行数学表达。由观测值内部可靠性的性质和前述的讨论, 基于可靠性的优化设计思想和计算步骤如下。
1) 一个网必须要有一定的多余观测, 多余观测数r 愈大, 则网的可靠性愈好, 但建网费用也愈高。
2) 在多余观测数一定的情况下, 观测值之间的精度相差不要太大, 边角观测值之间的精度应基本匹配。对于边角全测的初始方案, 也可根据边角观测值的平均多余观测分量来判断边角精度匹配的情况。
3) 根据网的设计要求、所使用的仪器、设计和实地踏勘, , 故有最大的多余观测数, 是一个“肥网”或“密网”。
4) 模拟初始观测方案, 进行平差计算, 对精度、可靠性乃至灵敏度计算结果并进行分析。首先, 检验确定的观测精度是否合理, 若不合理, 则需作适当调整。在观测值精度基本合理的基础上, 基于观测值内部可靠性指标按从“肥”到“瘦”、从“密”到“疏”的策略[5]进行网的优化设计。具
体地说, 先确定一个恰当的 r i 值, 再按式(8) 计算
00
n , 根据差值(n -n ) 确定应删去的观测值个数
4 优化设计软件
本文所提出的优化设计方法可使用任何一种网平差软件进行, 但要求软件还具有计算观测值的多余观测分量和模拟观测方案的功能。
科傻二[6]能按严密的间接观测模型对一、二维网进行从一个到数千个未知点, 包括上万个观测值的等权或不等权平差。除平差功能外, 还可自动推算网点近似坐标, 自动生成最小独立闭合环[7], 进行贯通误差影响值估算、坐标换算、粗差探测与定值定位、方差分量估算等[8]。该软件包与网的优化设计有关的功能是:能作网的模拟计算, 计算观测值的多余观测分量, 根据需要输出各种精度评定结果, 。
:网。l , 边长固定误差l (mm ) ,
比例误差l (ppm ) , 精度号l
……
点名, 点类型, X 坐标(m ) , Y 坐标(m ) ……
测站点, 照准点, A , 方位角值测站点点号
照准点点号, 观测值类型, 精度号
……其中, Ⅰ为观测精度部分; Ⅱ为控制点坐标和已知方位角部分; Ⅲ为观测方案部分。
某模拟桥梁边角网的OB2文件为:1. 8, 3, 2(, m 1)
(若只有一组观测精度, 则仅有第一行, 且m 1可省去。第2行如:2. 5, 5, 3, m 2)
n q 。对计算的观测值多余观测分量按从大到小
的顺序排列, 删去n d 个多余观测分量较大的那
些观测值, 然后重新作观测值模拟计算。如果观测值精度选择恰当, 仅作一二次迭代计算即可得到网的优化设计方案。
该方法的特点是:初始方案是一个观测精度和观测值个数都有富裕的全边角网, 整个优化设计过程中的关键是如何删除多余观测和调整观测精度。按此法删除的多余观测具有确定性且不致于引起形亏(这是该法的最大优点和特点) 。初始方案中每测站上方向观测的零方向选择要恰当, 最好选距离偏短的边。对于同一测站上的边和方向, 还可根据通视条件和外界环境的影响给不同的精度, 使结果更符合实际。
1, 0(表示已知点) , 5000, 50002, 1(表示未知点) , 6000, 5000
……
1, 2, A , 0. 0(表示1至2点的已知方位角为
) 0°00′00. 0″
1(下面按顺时针方向输照准点)
S W , L (表示方向观测) , 1(与m 1对应, 若只
有一组观测精度, 可省去)
S W , S (表示边长观测) , 1N W , L , 1N W , S , 1
2, L , 12, S , 1
N E , L , 1
武汉大学学报・信息科学版 2001年358
N E , S , 1M E , L , 1M E , L , 1S E , L , 1S E , S , 1度m p max ≤4. 5mm , 最弱边精度m s /s 优于
1/(12万) 。初始方案为观测50个方向值, 25条对向边的全边角网(如图1(a ) 所示) 。观测值精度为:m r =1. 8″, a =3mm , b =2×10-6。初始方案的计算结果为:n =75, t =21, r =54, r i =0. 72, r li =0. 76, r si =0. 64(边长观测值较方向观测值的多余观测分量要小, 说明边长的精度要高于方向的精度) 。点位精度均值m p =2. 0mm , 最大点位精度m p max =2. 3mm , m s /s =1/(17. 5万) , 显然质量太高。
按前述步骤进行如下优化设计与计算。
00
1) 取 r i =0. 5, 由式(11) 计算得n =42, 应删去n d =75-42=33个观测值。按观测值的多余观测分量从大到小排序, 确定删除r li ≥0. 80的观测值22个, r si ≥0. 6611个, 重新作, : r . si 0, m p =3. 1mm , p . , m s (. 2万) , 网的质量仍0
2) 取 r 0i =0. 42进行计算, 得n =36, n d =39。在前述的基础上, 再删去3个r li >0. 70和3
2(此测站下面的输入与第1站完全相同)
…
程序可根据网点坐标和观测方案自动生成平差所需要的输入文件(网名・IN2) , 文件中的观测值系根据坐标反算后加上模拟的观测值精度。为此, 在程序中, 提供了生成正态标准随机数的功能, 对这些正态标准随机数还进行有各种检验(如正态性、周期误差、偏度、峰度、方差、均值检验等) 。
模拟计算得到的观测值文件(网名・IN2) 与实测得到的完全一致, 因不含粗差和系统误差, 故更符合最小二乘平差要求。种精度和可靠性指标, 的优化设计5 算 例
5. 1 算例1 模拟桥梁边角网
个r si >0. 59共6个方向和边长观测值, 计算得
r li =0. 47, r si =0. 29, m p =3. 2mm , m p max =3. 7mm , m s /s =1/(13. 9万) 。网的质量仍偏高。
3) 由于 r i 已达到设计的极限值, 为此, 可根
该网共8个点, 其中点1、2位于桥轴线上, 北
岸4个点, 南岸3个点, 江中岛上1个点。以1为已知点, 1至2的方位角为已知方位角, 按独立网
进行设计。优化设计要求为: r 1≥0. 4, 最弱点精
据测量仪器和作业习惯, 改变观测值的精度, 例如
让测边的精度保持不变, 而m r 降为2. 5″, 仍按前面优化计算所得到的观测方案进行计算
(a ) 初始方案 (b ) 优化方案
图1 模拟桥梁边角网初始方案和优化方案
Fig. 1 The Simulated Triangu-trilateration Network for Bridge Project Primary Design and Optrmization Design
第4期 张正禄等:一种基于可靠性的工程控制网优化设计新方法 359
得 r li =0. 5, r si =0. 23, m p =3. 6mm , m p max =4. 3mm , m s /s =1/(12. 1万) 。该结果满足设计的要求, 所对应的为网的优化设计方案(参见图1(b ) ) 。优化设计方案相对于初始方案, 共减少了39个观测值, 方向观测精度也从1. 8″降低到2. 5″, 达到综合最优。
观测, 6个点上只有边长观测, 共有98个方向观
测值, 88个边长观测值。观测精度为:m r =1. 0″,
(7) 计算得:。按公式(6) 、t =55, r =131, r i =0. 70, r li =0. 69, r si =0. 72(有多条边的多余观测分量为1) 。由于该网是约束网, 可在相同基准下用点位精度和最弱边相对精度来衡量网的质量。优化方案取 r 0i =0. 54, 由a =1mm , b =1×10
-6
需要说明的是, 改变起始点和起始方位角, 例
如以S W 为已知点, S W-1为已知方位角, 仍按网的优化设计方案计算, 所有的可靠性指标不变, 边长相对精度也不变, 但点位精度改变了, m p 从3. 6mm 增大到8. 7mm , m p max 从4. 3mm 变大到12. 8mm 。
5. 2 算例2 某大型实测施工控制网
式(11) 计算得n 0=119, 应去掉67个观测值。将
计算出的r li 和r si 按大小顺序排列, 结合经验和作业习惯, 去掉31个方向观测值, 36个边长观测值, 作为优化方案, 该方案的观测精度不变, r li =0. 54, r si =0. 53。我们将实测方案和优化方案进行比较结果列于表2, 表中的精度系根据验后单位数中误差计算。
1, 本文
该网总点数27个, 其中已知点10个。在每
个点上都设站观测, 其中在21个点上进行了方向
表2 Tab. 2 The Real Observation Network 案
况
k
n L
n s
r
精 度/mm
m p
m p max
m s
s/m
s/m s
r i
0. 700. 54
优化方案
21
1010
186119
9867
8852
13164
1. 31. 5
1. 71. 9
0. 80. 8
11. 3211. 32
1400014000
所提出的方法简单易行, 优化结果很少依赖人们的主观认识, 一致性较好。从初始方案到优化设计方案只需作一二次迭代计算即可。根据本文所提出的算法和公式, 容易编写网优化设计的自动计算程序。从算例2可见, 该方法可用于对现有网的评价, 优化方案与实测方案比较, 在精度降低非常小的情况下, 其工作量和建网费用的减小却是十分显著的。
译. 北京:测绘出版社,1989
2 李德仁. 误差处理和可靠性理论. 北京:测绘出版社,
1988
3 陈永奇, 张正禄, 吴子安, 等. 高等应用测量. 武汉:武
汉测绘科技大学出版社,1996
4 张正禄, 李晓东. 观测值的灵敏度影响系数及其在监
测网优化设计中的应用. 测绘学报,1989,18(4) :249～
257
5 张正禄. 工程平面控制网模拟法优化设计策略. 武测
科技,1993(1) :1～4
6 结 语
工程控制网特别是施工控制网、变形监测网的优化设计仍然是一个值得研究和应用的问题。本文在深入研究观测值内部可靠性指标的基础上, 归纳了观测值多余观测分量与精度、建网费用乃至灵敏度等质量指标的关系, 提出了基于观测值可靠性的模拟法优化设计的思想和算法, 指出了该方法的特点和优点, 并在通用平差软件基础上研制了相应的优化设计程序。通过算例说明了网的优化设计过程, 用该思想和方法, 还可在对现有网平差的同时, 评价网在设计上的好坏。
参 考 文 献
1 佩尔策. 现代工程测量控制网的理论和应用. 张正禄
6 张正禄, 黄全义, 罗年学, 等. 全站式地面测量工程一
体化自动化系统研究. 武汉测绘科技大学学报,1999,
24(1) :79～83
7 冯 琰, 张正禄, 罗年学. 最小独立闭合环与附合导线
的自动生成算法. 武汉测绘科技大学学报, 1998, 23
(3) :255～259
8 罗年学, 张正禄, 黄全义, 等. WIN95版现代测量控制
网数据处理通用软件的若干技术和算法. 见:姚楚光, 张正禄. 水利水电测绘科技论文集, 武汉:武汉测绘科技大学出版社,1998. 169～174
作者简介:张正禄, 教授, 博士生导师。现主要从事精密工程测量、变形监测分析与预报、测量数据处理和工程信息系统方面的科研和教学工作。主要代表成果有:地面测量工程内外业一体化和数据处理自动化系统; 超站式GPS 和全站仪集成测绘系统; 测量机器人、GPS 和多传感器专用仪器应用研究。已发表论文120余篇。
E -mail :[email protected]. cn
武汉大学学报・信息科学版 2001年360
A N e w Method B ased on R eliability of Engineering Control
N et work for N etsoptimization
ZHA N G Zhengl u L U O N ianx ue HUA N G Q uanyi M EI Wensheng
(1 School of G eodesy and G eomatics ,Wuhan University ,129Luoyu Road ,Wuhan ,China ,430079)
1
1
1
1
Abstract :The optimization of engineering control network , especially that of the deformation monitoring network and precise control network for large construction is a very important work in network design. But in practice it is still difficult for the users to realize the network optimization with analytical or analogical method. In this paper ,after interpretation on the optimization design methods of engineering control network and the quality standards such as precision , reliability , sensitivity and costs ,a new algorithm of optimization design with simulative method for engineer 2ing control network based on the inner reliability of observations is presented. The be 2tween inner reliability and other quality standards is and characteristics of this new method are that the and not on the knowledge and ,such as aver 2age and matching criterion for the preci 2sion of ,are proposed and defined. A general software packet for
data and adjustment of surveying network (cosa ) used for the network optimization de 2sign with this method is introduced. With this software ,a primary network observation plan is simulated by 3. OB2file. Then the corespondent observation file 3. IN2will be created automat 2ically according to the above 3. OB2file. After adjustment ,the redundant of every observation will be calculated. Based on the average redundant and observation design number ,surveyors may make a decision that which observations could be deleted. Finally ,two examples are given. One is a simulated triangu-trilateration network for bridge project. The work load of the final optimal mea 2suring plan is reduced to 36from the initial 75observations and the precision and reliability are still satisfied for the design. The other is a real large construction network. For this network ,there are 186observations at the beginning. After the optimal processing the observations are reduced to 119and the precision of the point with the lowest precision is not degraded significantly. The opti 2mization benefits and necessity can be seen clearly through the above examples. Furthermore ,this method can not only make the optimization of network ,but also be used to interpret the rationality of an existing network.
K ey w ords :engineering control network ;optimization design ;reliability ;observations redundant
About the author :ZH ANG Zhenglu , professor , Ph. D supervisor. He is concentrated on the research and education in precise engineering geodesy ,deformation monitoring analysis and forecast ,measurement data processing and engineering geo-in formation system. He has made contribution in the areas of autormatic system of terrestrial control surveying data gathering and processing integrated super GPS and total sta 2tion measurement sytem ,the application of measurement robot , GPS and special instrument with multi-sensors. H is published papers are more than 120.
E -mail :[email protected]. cn
第26卷第4期2001年8月武汉大学学报・信息科学版
G eomatics and Information Science of Wuhan University Vol. 26No. 4Aug. 2001
文章编号:1000-050X(2001) 04-0354-07文献标识码:A
一种基于可靠性的工程控制网优化设计新方法
张正禄1 罗年学1 黄全义1 梅文胜1
(1 武汉大学测绘科学与技术学院, 武汉市珞喻路129号,430079)
摘 要:在论述工程控制网优化设计方法及网的精度、可靠性、灵敏度以及费用等准则的基础上, 提出了一种基于观测值内部可靠性指标的工程控制网模拟法优化设计的新思想和算法, 指出了该方法的优点和特点; 介绍了该优化设计方法所采用的“科傻”软件, 并用实例说明了用该方法进行工程控制网优化设计的计算步骤、优化效益以及优化设计的必要性。
关键词:工程控制网; 优化设计; 可靠性; 多余观测分量中图法分类号:P221. 2; P258
按用途划分, 施工控制网, () 和约束网() , 独立网只固定一个点和一个方向, 即已知一点的坐标和一条边的方位角, 网的尺度由测边确定; 约束网有两个或两个以上已知点。一般工程控制网以采用独立网为宜, 已知点和已知方位角可以在假定坐标系中给定, 也可以采用国家或地方坐标系的一个已知点和一个已知方位角。
网的优化设计有解析法和模拟法两种。解析法是以最优化理论为基础的严密方法, 其数学模型一般表示为:
f (X ) →min φ(X ) ≥0(1) ψ(X ) =0式中, 第一式为目标函数; 第二、三式为约束条件。优化设计的实质是通过给定的约束条件求目标函数的极值而得到最优解。目标函数可以是精度、可靠性、灵敏度或费用等指标, 约束条件也可以是上述指标。有人总以为解析法的结果是最优的, 实际上并非如此。因为目标函数和约束条件都是一种经过简化的数学模型, 不能完全反映客观实际。最有代表性的例子是以费用为目标函数或约束条件时, 都是以观测值的权的总和来代表建网费用, 实际上这是很粗糙的。而且, 绝大多数情况
收稿日期:2001-02-23。
项目来源:国家测绘局教学研究基金资助项目(97GH01) 。
, 优化后还存在取整的问题, 同一测站上还要划为便于操
作的等权观测问题, 经过取整和调整后的观测方案不再是最优的了。解析法不但存在上述缺点, 同时还因优化任务不同, 目标函数和约束条件的数学模型也不尽相同, 故研制通用的解析法优化设计软件的难度较大。尽管解析法优化设计方面的研究已很多, 却难以推广, 在实际生产中也鲜有应用。
模拟法实际上是一种试算法。一般作法是:根据优化的任务和设计者的知识和经验, 制定初始设计方案, 用模拟的观测值对该方案作平差计算, 对平差结果作评价, 对初始方案进行修改, 再计算、再修改, 如此多次重复, 直至认为满意为止。模拟法需要一个好的平差程序, 并且该程序还能作观测值模拟, 显示网图及各种精度指标和可靠性数值指标, 便于人机对话式操作, 算法好(如采用序贯算法) , 速度快。但模拟法即使有一个好的优化设计程序, 仍依赖于设计者的知识和经验, 同样的要求, 甚至同样的初始设计方案, 不同人设计出的最终方案也有差别, 甚至差别还很大。由于解析法和模拟法的缺陷, 造成在实际工作中人们很少对工程控制网作优化设计, 大多根据测量负责人的知识、经验和习惯进行布网和观测。对于高精度施工控制网或变形监测网, 往往精度和费用偏高; 而对于测图控制网, 大多用多级
第4期 张正禄等:一种基于可靠性的工程控制网优化设计新方法 355
附合导线进行加密。这些作法虽然都符合规范要求, 但谈不上网的优化设计, 有不少弊病。
针对上述情况, 本文提出了一种基于观测值内部可靠性指标的工程控制网优化设计方法。该方法本质上属于模拟法, 但加进可靠性指标后, 具有量化的优化设计准则, 不以人的知识和经验为转移, 优化结果既具有一致性, 也不失严密性。该方法使用笔者研制的现代测量控制网数据处理通用软件包———科傻系统系列软件之二———Cosa2(或称科达普施Codaps ) 。
低, 当r i 等于1时, 则该观测值完全成为多余, 即使删除之网平差结果也不变。
2) r i 愈小, 该观测值的粗差愈难被发现, r i =0时, 即使有大的粗差或错误也无法发现。根据内部可靠性与外部可靠性具有一致性的性质[3], 平差结果受粗差的影响随r i 的减小而增大。r i 愈大, 则较小的粗差也能发现, r i =1时, 观测值的粗差能完全确定, 粗差对平差结果的影响随r i 的增大而减小。
3) 对于一个确定的网和设计方案, 即在网形和网的观测值总数n (或多余观测数r ) 确定的情况下, 观测值之间的精度相差愈大, 则内部可靠性r i 的值相差也愈大。观测值的精度愈高, 则相应的r i 愈小; 观测值的精度愈低, 其r i 愈大, 即观。
4) () , 多余观测r i 也愈大, 反之, 若r i 大, 则
r 也必定大(建网费用将随r 的增大而增加) 。
1 观测值的内部可靠性及其性质
网的可靠性可分为网的总体可靠性、观测值的内部可靠性(也称局部可靠性) 和外部可靠性。在此, 主要讨论观测值的内部可靠性。
对于一个工程控制网来说, 2-1
) l 差模型(l , A x , σ0P 性r i 为[1]:
=(V V P )
ii
(2)
5) 对于独立网来说, 与边长和角度是不变量
且满足
n
∑r
1
一样, 观测值的内部可靠性r i 是与基准的位置无
关的不变量。
i
=r =n -t
r i
(3) (4)
ω0/ 0l i =σi ・
2 几个有关问题
2. 1 多余观测数和平均多余观测分量
式中, ω0为非中心参数, 对于单个观测值粗差而
言, 其取值与显著水平α和检验功效γ有关, 常取2. 79(α=0. 05, γ=0. 80) 或4. 13(α=0. 001, γ=0. 80) [2]; r 为网的多余观测数, r i 为观测值l i 的多余观测分量, r i 可以反映控制网发现观测值l i (中误差为σ通过统计检验, r i i ) 中粗差的能力。愈大, l i 中能被发现粗差的下界值 l i 愈小; 或对同一个粗差, 检验功率愈大。因此, r i 被定义为观测值的内部可靠性。
观测值相互独立(本文只讨论这种情况) 时, 有:
2
(5) r i =1-σ^2i /σi 式中, σ^i 为l i 的平差值的中误差。若观测值l i 的
精度很高, 则平差后的精度提高很小, 有σ^i ≈σi , 这时r i →0。若无多余观测, 有σ^i =σi , r i =0; 反之, 若l i 的精度很低, 则平差后精度将显著提高, 即σ^i
1) 0≤r i ≤1。r i 愈小, 该观测值在网中的地
设控制网的网点数为m , 已知点数为m k , 未
知点数为m u , 已知边和已知方位角数分别为k s 和k a , 进行有方向观测的设站数为m l , 只进行有边长观测的设站数为m s , 其中有方向或边长观测的已知点设站数为m 1, 总的观测值数为n (其中, 方向观测值数为n l , 边长观测值数为n s ) , 必要观测值个数为t , 多余观测数为r , 平均多余观测分量为 r i , 其中, 方向和边长的平均多余观测分量分别为 r li 、 r si , 方向、边长的最大最小多余观测分量分别为r l max 、r l min 、r s max 、r s min 。某观测值的多余观测分量用r i 、r li 、r si 表示。对于按方向的间接观测平差模型, 网的多余观测数r 和平均多余观测分量可按下式计算:
r =n -t =n -(3m l +2m s -2m 1-k s -k a )
(6) (7)
r =
n
位愈高, 若r i 等于零, 则该观测值不可缺少, 否则将产生形亏。r i 愈大, 该观测值在网中的地位愈
例如, 对于附合导线, 有m 1=2, m l =m u +2, 方向观测值数n l =2(m u +2) , 边长观测值数n s =m u +1, 未知数个数t =3(m u +2) -2×2, 故r =n l +n s
武汉大学学报・信息科学版 2001年356
-t =2m u +4+m u +1-(3m u +6-4) =3。
表1列出了附合导线含不同未知点数时观测值个数和平均多余观测分量。当未知点多于4个
时, r i 小于0. 2。从可靠性指标来看, 附合导线不能作为施工控制网和变形监测网。 在优化设计时, 为了使网有足够的可靠性, 在
表1 附合导线的观测值个数和平均多余观测分量
Tab. 1 Observation Number and Average Ebservation Redundant for Traverses 未知点数
n
180. 375
2110. 27
3140. 21
4170. 18
5200. 15
6230. 13
7260. 115
8290. 10
9320. 09
r i
此提出平均多余观测分量设计值的概念, 用 r i 表示。对于不同用途、不同要求的网, 可以取 r i 的00设计值 r i 为0. 3～0. 6。由 r i 值, 可按下式计算
验值。若要求提高, 还可取小一些, 如k =1～
2. 0。
2. 3 网的质量准则问题
观测值数的设计值n 0, 有:
n
=
1- r i
(8)
令
n i 0=1- r i 0
t =n i
网的质量准则主要是精度、可靠性、建网费
用, 对于变形监测网还包括灵敏度和可区分性。
、相对点位精(或主分量) , 点位精度与基准的位置有关, 对于独立网来说, 最好是将最靠近网的重心的点作为已知点, 以通过该点的最接近中心线的方向作为起始方向。这样可保证点位精度在数值上达到最小。实际上, 应将相对点位精度或最弱边精度作为精度准则, 因为它们是与基准的位置无关的不变量。在模拟法优化设计中, 应取先验单位权中误差计算各精度指标。
网的可靠性指标建议采用平均多余观测分量 r i 衡量, 而且 r si 和 r li 相差不应太大。对于相同的观测方案, r si 和 r li 随边角精度的改变而改变, 例如若提高方向观测的精度, 则 r li 将减小。同时, (r l max -r l min ) 与(r s max -r s min ) 也不能太大。
由§1中观测值内部可靠性的性质4) 可知, 由于r i 的总和反映了多余观测数, 因此, 也在很大程度反映了建网费用, 即r 愈大, 建网费用愈高。
网的灵敏度表现为网点在特定方向上的精度。只有在特殊情况下才需要作灵敏度计算。一般可以采用误差椭圆近似地进行评价, 即网点在要求的方向上误差应较小, 如尽量靠近误差椭圆的短轴方向。
n i 0为观测值的平均必要观测分量, :2. 2 , 纯测角网已逐渐消亡, 纯测边网也很少采用, 广泛应用的边角网或导线网存在边角精度的匹配问题。众所周知, 测边引起(或控制) 纵向误差, 方向引起(或控制) 横向误差, 它们都与边长密切相关。设方向中误差为m ″r , 测边的固定误差和比例误差分别为a 和b , 边长为s , 则由方向中误差引起的横向误差为:
m (9) m q =・S
ρ″
由边长中误差引起的纵向误差为:
m L =
22
a +(bs ) 或 m L =a +bs
(10)
所谓边角精度完全匹配, 是指应满足m q =m L 。由于网的边长变化和仪器的限制, 边角精度匹配是相对的, 不匹配是绝对的。一般认为, 当
m q 和m L 满足下述关系
k
m L ≤m q ≤km L , (k =1～2. 5) (11)
时, 边角精度是基本匹配的。因为观测值的精度与其可靠性成反比, 边角精度是否匹配也可通过可靠性来判别。对于初始方案即边角全测的情况(这时n l ≈2n s ) , 通过计算方向和边长的多余观测数和平均多余观测分量r l 、r s 、 r li 、 r si , 可根据 r li 、 r si 的大小来判断边角精度是否匹配, 一般应满足 r li ≤ r si ≤k r li , k 一般在1～2. 5之间。这
k
里的k 和式(11) 中的k 值无直接联系, 且都是经
3 优化设计思想、计算步骤和特点
观测值的内部可靠性与观测值的精度有密切关系, 而观测值的精度又与建网费用有关。变形监测网的灵敏度实际是网点在特定方向上的精度, 它也取决于网的观测方案设计和观测值的精
第4期 张正禄等:一种基于可靠性的工程控制网优化设计新方法 357
度, 此外, 变形与粗差的可区分性也必然涉及到观测值的精度。因此, 观测值的内部可靠性与观测值的精度、建网费用、监测网的灵敏度和可区分性存在密切的关系, 有的关系可以简化为如式(3) 、(4) 、(5) 或文献[4]所描述的数学公式, 有的关系却难以进行数学表达。由观测值内部可靠性的性质和前述的讨论, 基于可靠性的优化设计思想和计算步骤如下。
1) 一个网必须要有一定的多余观测, 多余观测数r 愈大, 则网的可靠性愈好, 但建网费用也愈高。
2) 在多余观测数一定的情况下, 观测值之间的精度相差不要太大, 边角观测值之间的精度应基本匹配。对于边角全测的初始方案, 也可根据边角观测值的平均多余观测分量来判断边角精度匹配的情况。
3) 根据网的设计要求、所使用的仪器、设计和实地踏勘, , 故有最大的多余观测数, 是一个“肥网”或“密网”。
4) 模拟初始观测方案, 进行平差计算, 对精度、可靠性乃至灵敏度计算结果并进行分析。首先, 检验确定的观测精度是否合理, 若不合理, 则需作适当调整。在观测值精度基本合理的基础上, 基于观测值内部可靠性指标按从“肥”到“瘦”、从“密”到“疏”的策略[5]进行网的优化设计。具
体地说, 先确定一个恰当的 r i 值, 再按式(8) 计算
00
n , 根据差值(n -n ) 确定应删去的观测值个数
4 优化设计软件
本文所提出的优化设计方法可使用任何一种网平差软件进行, 但要求软件还具有计算观测值的多余观测分量和模拟观测方案的功能。
科傻二[6]能按严密的间接观测模型对一、二维网进行从一个到数千个未知点, 包括上万个观测值的等权或不等权平差。除平差功能外, 还可自动推算网点近似坐标, 自动生成最小独立闭合环[7], 进行贯通误差影响值估算、坐标换算、粗差探测与定值定位、方差分量估算等[8]。该软件包与网的优化设计有关的功能是:能作网的模拟计算, 计算观测值的多余观测分量, 根据需要输出各种精度评定结果, 。
:网。l , 边长固定误差l (mm ) ,
比例误差l (ppm ) , 精度号l
……
点名, 点类型, X 坐标(m ) , Y 坐标(m ) ……
测站点, 照准点, A , 方位角值测站点点号
照准点点号, 观测值类型, 精度号
……其中, Ⅰ为观测精度部分; Ⅱ为控制点坐标和已知方位角部分; Ⅲ为观测方案部分。
某模拟桥梁边角网的OB2文件为:1. 8, 3, 2(, m 1)
(若只有一组观测精度, 则仅有第一行, 且m 1可省去。第2行如:2. 5, 5, 3, m 2)
n q 。对计算的观测值多余观测分量按从大到小
的顺序排列, 删去n d 个多余观测分量较大的那
些观测值, 然后重新作观测值模拟计算。如果观测值精度选择恰当, 仅作一二次迭代计算即可得到网的优化设计方案。
该方法的特点是:初始方案是一个观测精度和观测值个数都有富裕的全边角网, 整个优化设计过程中的关键是如何删除多余观测和调整观测精度。按此法删除的多余观测具有确定性且不致于引起形亏(这是该法的最大优点和特点) 。初始方案中每测站上方向观测的零方向选择要恰当, 最好选距离偏短的边。对于同一测站上的边和方向, 还可根据通视条件和外界环境的影响给不同的精度, 使结果更符合实际。
1, 0(表示已知点) , 5000, 50002, 1(表示未知点) , 6000, 5000
……
1, 2, A , 0. 0(表示1至2点的已知方位角为
) 0°00′00. 0″
1(下面按顺时针方向输照准点)
S W , L (表示方向观测) , 1(与m 1对应, 若只
有一组观测精度, 可省去)
S W , S (表示边长观测) , 1N W , L , 1N W , S , 1
2, L , 12, S , 1
N E , L , 1
武汉大学学报・信息科学版 2001年358
N E , S , 1M E , L , 1M E , L , 1S E , L , 1S E , S , 1度m p max ≤4. 5mm , 最弱边精度m s /s 优于
1/(12万) 。初始方案为观测50个方向值, 25条对向边的全边角网(如图1(a ) 所示) 。观测值精度为:m r =1. 8″, a =3mm , b =2×10-6。初始方案的计算结果为:n =75, t =21, r =54, r i =0. 72, r li =0. 76, r si =0. 64(边长观测值较方向观测值的多余观测分量要小, 说明边长的精度要高于方向的精度) 。点位精度均值m p =2. 0mm , 最大点位精度m p max =2. 3mm , m s /s =1/(17. 5万) , 显然质量太高。
按前述步骤进行如下优化设计与计算。
00
1) 取 r i =0. 5, 由式(11) 计算得n =42, 应删去n d =75-42=33个观测值。按观测值的多余观测分量从大到小排序, 确定删除r li ≥0. 80的观测值22个, r si ≥0. 6611个, 重新作, : r . si 0, m p =3. 1mm , p . , m s (. 2万) , 网的质量仍0
2) 取 r 0i =0. 42进行计算, 得n =36, n d =39。在前述的基础上, 再删去3个r li >0. 70和3
2(此测站下面的输入与第1站完全相同)
…
程序可根据网点坐标和观测方案自动生成平差所需要的输入文件(网名・IN2) , 文件中的观测值系根据坐标反算后加上模拟的观测值精度。为此, 在程序中, 提供了生成正态标准随机数的功能, 对这些正态标准随机数还进行有各种检验(如正态性、周期误差、偏度、峰度、方差、均值检验等) 。
模拟计算得到的观测值文件(网名・IN2) 与实测得到的完全一致, 因不含粗差和系统误差, 故更符合最小二乘平差要求。种精度和可靠性指标, 的优化设计5 算 例
5. 1 算例1 模拟桥梁边角网
个r si >0. 59共6个方向和边长观测值, 计算得
r li =0. 47, r si =0. 29, m p =3. 2mm , m p max =3. 7mm , m s /s =1/(13. 9万) 。网的质量仍偏高。
3) 由于 r i 已达到设计的极限值, 为此, 可根
该网共8个点, 其中点1、2位于桥轴线上, 北
岸4个点, 南岸3个点, 江中岛上1个点。以1为已知点, 1至2的方位角为已知方位角, 按独立网
进行设计。优化设计要求为: r 1≥0. 4, 最弱点精
据测量仪器和作业习惯, 改变观测值的精度, 例如
让测边的精度保持不变, 而m r 降为2. 5″, 仍按前面优化计算所得到的观测方案进行计算
(a ) 初始方案 (b ) 优化方案
图1 模拟桥梁边角网初始方案和优化方案
Fig. 1 The Simulated Triangu-trilateration Network for Bridge Project Primary Design and Optrmization Design
第4期 张正禄等:一种基于可靠性的工程控制网优化设计新方法 359
得 r li =0. 5, r si =0. 23, m p =3. 6mm , m p max =4. 3mm , m s /s =1/(12. 1万) 。该结果满足设计的要求, 所对应的为网的优化设计方案(参见图1(b ) ) 。优化设计方案相对于初始方案, 共减少了39个观测值, 方向观测精度也从1. 8″降低到2. 5″, 达到综合最优。
观测, 6个点上只有边长观测, 共有98个方向观
测值, 88个边长观测值。观测精度为:m r =1. 0″,
(7) 计算得:。按公式(6) 、t =55, r =131, r i =0. 70, r li =0. 69, r si =0. 72(有多条边的多余观测分量为1) 。由于该网是约束网, 可在相同基准下用点位精度和最弱边相对精度来衡量网的质量。优化方案取 r 0i =0. 54, 由a =1mm , b =1×10
-6
需要说明的是, 改变起始点和起始方位角, 例
如以S W 为已知点, S W-1为已知方位角, 仍按网的优化设计方案计算, 所有的可靠性指标不变, 边长相对精度也不变, 但点位精度改变了, m p 从3. 6mm 增大到8. 7mm , m p max 从4. 3mm 变大到12. 8mm 。
5. 2 算例2 某大型实测施工控制网
式(11) 计算得n 0=119, 应去掉67个观测值。将
计算出的r li 和r si 按大小顺序排列, 结合经验和作业习惯, 去掉31个方向观测值, 36个边长观测值, 作为优化方案, 该方案的观测精度不变, r li =0. 54, r si =0. 53。我们将实测方案和优化方案进行比较结果列于表2, 表中的精度系根据验后单位数中误差计算。
1, 本文
该网总点数27个, 其中已知点10个。在每
个点上都设站观测, 其中在21个点上进行了方向
表2 Tab. 2 The Real Observation Network 案
况
k
n L
n s
r
精 度/mm
m p
m p max
m s
s/m
s/m s
r i
0. 700. 54
优化方案
21
1010
186119
9867
8852
13164
1. 31. 5
1. 71. 9
0. 80. 8
11. 3211. 32
1400014000
所提出的方法简单易行, 优化结果很少依赖人们的主观认识, 一致性较好。从初始方案到优化设计方案只需作一二次迭代计算即可。根据本文所提出的算法和公式, 容易编写网优化设计的自动计算程序。从算例2可见, 该方法可用于对现有网的评价, 优化方案与实测方案比较, 在精度降低非常小的情况下, 其工作量和建网费用的减小却是十分显著的。
译. 北京:测绘出版社,1989
2 李德仁. 误差处理和可靠性理论. 北京:测绘出版社,
1988
3 陈永奇, 张正禄, 吴子安, 等. 高等应用测量. 武汉:武
汉测绘科技大学出版社,1996
4 张正禄, 李晓东. 观测值的灵敏度影响系数及其在监
测网优化设计中的应用. 测绘学报,1989,18(4) :249～
257
5 张正禄. 工程平面控制网模拟法优化设计策略. 武测
科技,1993(1) :1～4
6 结 语
工程控制网特别是施工控制网、变形监测网的优化设计仍然是一个值得研究和应用的问题。本文在深入研究观测值内部可靠性指标的基础上, 归纳了观测值多余观测分量与精度、建网费用乃至灵敏度等质量指标的关系, 提出了基于观测值可靠性的模拟法优化设计的思想和算法, 指出了该方法的特点和优点, 并在通用平差软件基础上研制了相应的优化设计程序。通过算例说明了网的优化设计过程, 用该思想和方法, 还可在对现有网平差的同时, 评价网在设计上的好坏。
参 考 文 献
1 佩尔策. 现代工程测量控制网的理论和应用. 张正禄
6 张正禄, 黄全义, 罗年学, 等. 全站式地面测量工程一
体化自动化系统研究. 武汉测绘科技大学学报,1999,
24(1) :79～83
7 冯 琰, 张正禄, 罗年学. 最小独立闭合环与附合导线
的自动生成算法. 武汉测绘科技大学学报, 1998, 23
(3) :255～259
8 罗年学, 张正禄, 黄全义, 等. WIN95版现代测量控制
网数据处理通用软件的若干技术和算法. 见:姚楚光, 张正禄. 水利水电测绘科技论文集, 武汉:武汉测绘科技大学出版社,1998. 169～174
作者简介:张正禄, 教授, 博士生导师。现主要从事精密工程测量、变形监测分析与预报、测量数据处理和工程信息系统方面的科研和教学工作。主要代表成果有:地面测量工程内外业一体化和数据处理自动化系统; 超站式GPS 和全站仪集成测绘系统; 测量机器人、GPS 和多传感器专用仪器应用研究。已发表论文120余篇。
E -mail :[email protected]. cn
武汉大学学报・信息科学版 2001年360
A N e w Method B ased on R eliability of Engineering Control
N et work for N etsoptimization
ZHA N G Zhengl u L U O N ianx ue HUA N G Q uanyi M EI Wensheng
(1 School of G eodesy and G eomatics ,Wuhan University ,129Luoyu Road ,Wuhan ,China ,430079)
1
1
1
1
Abstract :The optimization of engineering control network , especially that of the deformation monitoring network and precise control network for large construction is a very important work in network design. But in practice it is still difficult for the users to realize the network optimization with analytical or analogical method. In this paper ,after interpretation on the optimization design methods of engineering control network and the quality standards such as precision , reliability , sensitivity and costs ,a new algorithm of optimization design with simulative method for engineer 2ing control network based on the inner reliability of observations is presented. The be 2tween inner reliability and other quality standards is and characteristics of this new method are that the and not on the knowledge and ,such as aver 2age and matching criterion for the preci 2sion of ,are proposed and defined. A general software packet for
data and adjustment of surveying network (cosa ) used for the network optimization de 2sign with this method is introduced. With this software ,a primary network observation plan is simulated by 3. OB2file. Then the corespondent observation file 3. IN2will be created automat 2ically according to the above 3. OB2file. After adjustment ,the redundant of every observation will be calculated. Based on the average redundant and observation design number ,surveyors may make a decision that which observations could be deleted. Finally ,two examples are given. One is a simulated triangu-trilateration network for bridge project. The work load of the final optimal mea 2suring plan is reduced to 36from the initial 75observations and the precision and reliability are still satisfied for the design. The other is a real large construction network. For this network ,there are 186observations at the beginning. After the optimal processing the observations are reduced to 119and the precision of the point with the lowest precision is not degraded significantly. The opti 2mization benefits and necessity can be seen clearly through the above examples. Furthermore ,this method can not only make the optimization of network ,but also be used to interpret the rationality of an existing network.
K ey w ords :engineering control network ;optimization design ;reliability ;observations redundant
About the author :ZH ANG Zhenglu , professor , Ph. D supervisor. He is concentrated on the research and education in precise engineering geodesy ,deformation monitoring analysis and forecast ,measurement data processing and engineering geo-in formation system. He has made contribution in the areas of autormatic system of terrestrial control surveying data gathering and processing integrated super GPS and total sta 2tion measurement sytem ,the application of measurement robot , GPS and special instrument with multi-sensors. H is published papers are more than 120.
E -mail :[email protected]. cn