一种基于可靠性的综合项目工程控制网优化设计新方法.doc

一种基于可靠性工程控制网优化设计新办法* 收稿日期 -2-23国家测绘局教学研究项目（编号97GH01）张正禄 罗年学 黄全义 梅文胜 巢佰崇 （武汉大学测绘科学与技术学院，武汉市珞瑜路129号，430079）摘要：在阐述工程控制网优化设计办法及网精度、可靠性、敏捷度以及费用等准则基本上，提出了一种基于观测值内部可靠性指标工程控制网模仿法优化设计新思想和算法，指出了该办法长处和特点；简介了该优化设计办法所采用“科傻”软件，并用实例阐明了用该办法进行工程控制网优化设计计算环节、优化效益以及优化设计必要性。核心词：工程控制网；优化设计；可靠性；多余观测分量1 概述按用途划分，工程控制网可分为测图控制网、施工控制网、安装控制网和变形监测网；按基准或已知数据划分，可分为独立网（或称典型自由网）和约束网（或称强制网），独立网只固定一种点和一种方向，即已知一点坐标和一条边方位角，网尺度由测边拟定；约束网有两个或两个以上已知点。普通工程控制网以采用独立网为宜，已知点和已知方位角可以在假定坐标系给定，也可以采用国家或地方坐标系一种已知点和一种已知方位角。工程控制网可以采用地面测量技术或空间测量技术（重要是GPS技术）建立，它优化设计是一种古老而又时新问题，所涉及内容非常广。本文重要讨论用地面测量技术建立平面施工控制网和变形监测网优化设计问题，且又以独立网为重点。网优化设计有解析法和模仿法两种，它们都需要依托程序在计算机上作大量和复杂计算，显然都是机助法。解析法是以最优化理论为基本严密办法，其数学模型普通表达为 （1）第一式称目的函数，第二、三式称约束条件。优化设计实质是在给定约束条件通过求目的函数极值而得到最优解。目的函数可以是精度、可靠性、敏捷度或费用等指标，约束条件也可以是上述指标。人们总觉得解析法成果是最优，事实上并非如此。由于目的函数和约束条件都是一种通过简化数学模型，总不能完全反映客观实际。最有代表性例子是以费用为目的函数或约束条件时，都是以观测值权总和来代表建网费用。事实上这是很粗糙。并且，绝大多数状况下，解析法是对观测值权进行最佳分派，优化后还存在取整问题，同一测站上还要划为便于操作等权观测问题。通过取整和调节后观测方案不再是最优了。解析法不但存在上述缺陷，同步还因优化任务不同，目的函数和约束条件数学模型也不尽相似，故研制通用解析法优化设计通用软件难度较大。尽管解析法优化设计方面研究已诸多，却难于推广，在实际生产中也鲜有应用。模仿法事实上是一种试算法。普通作法是：依照优化任务和设计者知识和经验，制定初始设计方案，用模仿观测值对该方案作平差计算，对平差成果作评价，对初始方案进行修改，再计算、再修改，如此多次重复，直至以为满意为止。从上可见，模仿法需要一种好平差程序，该程序还能作观测值模仿，能显示网图及各种精度指标和可靠性数值指标，能便于人机对话式操作，算法好（如采用序贯算法），速度快。模仿法缺陷是：虽然有一种好优化设计程序，但仍依赖于设计者知识和经验，同样规定，甚至同样初始设计方案，不同人设计出最后方案也有差别，甚至差别还很大。由于解析法和模仿法缺陷，导致在实际工作中，人们很少对工程控制网作优化设计，大多依照测量负责人知识、经验和习惯布网和观测，对于高精度施工控制网或变形监测网，往往精度和费用偏高，而对于测图控制网，大多用多级附合导线进行加密，这些作法虽然都符合规范规定，但谈不上网优化设计。事实上有不少弊病。针对上述状况，本文提出了一种基于观测值内部可靠性指标工程控制网优化设计办法。该办法在本质上属于模仿法，但加进可靠性指标后，具备量化优化设计准则，不以人知识和经验为转移，优化成果既具备一致性，也不失严密性。该法所使用软件是咱们所研制当代测量控制网数据解决通用软件包科傻系统系列软件之二Cosa2（或称科达普施Codaps）。2 观测值内部可靠性及其性质网可靠性可分网总体可靠性，观测值内部可靠性（也称局部可靠性）和外部可靠性。在此，重要讨论观测值内部可靠性。对于一种工程控制网来说，由其间接观测平差模型（）可得观测值li内部可靠性ri为1： （2）且满足 （3） （4）上述公式中，0为非中心参数，对于单个观测值粗差而言，其取值与明显水平和检查功能关于，常取2.79（=0.05 =0.80）或4.13（=0.001 =0.80）2。ri为观测值li多余观测分量，r为网多余观测数，ri可以反映控制网发现观测值li（中误差为i）中粗差能力。ri愈大，通过记录检查，能发现li中粗差下界值愈小；或对同一种粗差，检查功率愈大。因而，ri被定义为观测值内部可靠性。观测值互相独立（本文只讨论这种状况）时，有 （5）上式中i为li平差值中误差。若观测值li精度很高,则平差后精度提高很小，有ii，这时ri0。特例，若无多余观测，有i=i，ri=0；反之，若li精度很低，则平差后精度将明显提高，即ii，这时ri1。特例，对已知边作边长观测，由于已知边误差等于零，即si=0，则rsi=1。观测值内部可靠性具备如下性质：1）0ri1。ri愈小，该观测值在网中地位愈高，若ri等于零，则该观测值不可缺少，否则将产生形亏。ri愈大，该观测值在网中地位愈低，当ri等于1，则该观测值完全成为多余，虽然删除其网平差成果也不变。2）ri愈小，该观测值粗差愈难被发现，ri=0时，虽然有大粗差或错误也无法发现，依照内部可靠性与外部可靠性具备一致性性质3，平差成果受粗差影响随ri减小而增大。ri愈大，则较小粗差也能发现，ri=1时，观测值粗差能完全拟定，粗差对平差成果影响随ri增大而减小。3）对于一种拟定网和设计方案，即在网形和网观测值总数n（或多余观测数r）拟定状况下，观测值之间精度相差愈大，则内部可靠性ri值相差也愈大。观测值精度愈高，则相应ri愈小，观测值精度愈低，其ri愈大。即观测值内部可靠性与观测值精度成反比。4）对于拟定网和观测精度（方向与边长观测精度重要由仪器和测回数决定）状况下，多余观测数r愈大，则观测值ri也愈大，反之，若ri大，则r也必然大（建网费用将随r增大而增长）。5）对于独立网来说，与边长和角度是不变量同样，观测值内部可靠性ri是与基准位置无关不变量。该性质可通过算例验证，其严格数学证明有待读者去作。3 几种关于问题31 多余观测数和平均多余观测分量设控制网网点数为m，已知点数为mk，未知点数为mu，已知边和已知方位角数分别为ks和ka（不是由已知点坐标反算得到），进行有方向观测设站数为ml，只进行有边长观测设站数为ms，其中有方向或边长观测已知点设站数为m1，总观测值数为n（其中，方向观测值数为nl，边长观测值数为ns），必要观测值个数为t，多余观测数为r，平均多余观测分量为，其中，方向和边长平均多余观测分量分别为，方向、边长最大最小多余观测分量分别为rlmax、rlmin、rsmax、rsmin，某观测值多余观测分量用ri、rli、rsi表达。对于按方向间接观测平差模型，网多余观测数r和平均多余观测分量可按下式计算r=n-t=n-(3ml+2ms-2m1-ks-ka) (6)= (7)例如，对于附合导线，有m1=2，ml=mu+2，方向观测值数nl=2（mu+2），边长观测值数ns=mu+1，未知数个数t=3（mu+2）-22，故r=nl+ns t=2mu+4+mu+1-（3mu+6-4）=3上述各种值均由程序进行自动计算。表1列出了附合导线含不同未知点数时观测值个数和平均多余观测分量。当未知点多于4个时，不大于0.2，从可靠性指标来看，附合导线不能作为施工控制网和变形监测网。在测图控制网中，采用附合导线进行逐级加密是值得注意。表1 附合导线观测值个数和平均多余观测分量未知点数123456789n811141720232629320.3750.270.210.180.150.130.1150.100.09在优化设计时，为了使网有足够可靠性，在此提出平均多余观测分量设计值概念，用表达。对于不同用途不同规定网，咱们可以取设计值为0.30.6。由值，可按下式计算观测值数设计值n0，有： （8）32 边角精度匹配问题及可靠性鉴别法 随着全站仪普遍应用，纯测角网已逐渐消灭，纯测边网也很少采用，广泛应用边角网或导线网存在边角精度匹配问题。咱们懂得，测边引起（或控制）纵向误差，方向引起（或控制）横向误差，它们都与边长密切有关。设方向中误差为mr，测边固定误差和比例误差分别为a和b，边长为s，则由方向中误差引起横向误差为 mrmq= S （9） q由边长中误差引起纵向误差为： 或mL=a+bs （10）所谓边角精度完全匹配，是指应满足mq=mL。由于网边长变化和仪器限制，边角精度匹配是相对，不匹配是绝对，普通以为当mq和mL满足下述关系： （11）式中，k=13时，都以为边角精度是基本匹配。由于观测值精度与其可靠性有成反比关系，边角精度与否匹配也可通过可靠性来鉴别，对于初始方案即边角全测状况（这时nl2ns），通过计算方向和边长多余观测数和平均多余观测分量rL、rs、，可依照、大小来判断边角精度与否匹配，普通应满足k，k普通在12.0之间。33 网质量准则问题网质量准则重要是精度、可靠性、建网费用，对于变形监测网还涉及敏捷度和可区别度。精度准则常采用有点位精度，相对点位精度（涉及误差椭圆）以及特性值以及主元（或主分量）等指标。实际应用中，计算点位精度和相对点位精度就足够了。值得注意是，点位精度与基准位置关于，对于独立网来说，最佳是将最接近网重心点作为已知点，以通过该点最接近中心线方向作为起始方向。这样可保证点位精度在数值上达到最小。事实上，应将相对点位精度或最弱边精度作为精度准则，由于它们是与基准位置无关不变量。记住：无论是独立网还是约束网，只有在相似基准下进行精度比较才故意义。并且，在模仿法优化设计中，应取先验单位数中误差计算各精度指标。网可靠性指标建议采用平均多余观测分量衡量，并且，和相差不应太大。对于相似观测方案，和随边角精度变化而变化，例如，若提高方向观测精度，则将减小。同步，（rlmax-rlmin）与（rsmax-rsmin）也不能太大。由观测值内部可靠性性质4）可知，由于ri总和反映了多余观测数，因而，也在很大限度反映了建网费用，即r愈大，建网费用愈高。网敏捷度，体现为网点在特定方向上精度，只有在特殊状况下才需要作敏捷度计算。普通，可以采用误差椭圆近似地进行评价，即网点在规定方向上误差应较小，如尽量接近误差椭圆短轴方向。4 基于可靠性优化设计思想、计算环节和特点观测值内部可靠性与观测值精度有密切关系，而观测值精度又与建网费用关于，并且，变形监测值网敏捷度实际是网点在特定方向上精度，它也取决于网观测方案设计和观测值精度，此外，变形与粗差可区别性也必然涉及到观测值精度。因而，观测值内部可靠性与观测值精度、建网费用、监测网敏捷度和可区别性存在密切关系，关于系可以简化为如（3）、（4）、（5）式或文献4所描述数学公式，关于系却难于进行数学表达。由观测值内部可靠性性质和前述讨论，咱们可以归纳为如下基于可靠性优化设计思想和计算环节如下：1）一种网必要要有一定多余观测，多余观测数r愈大，则网可靠性愈好，但建网费用也愈高。2）在多余观测数一定状况下，观测值之间精度相差不要太大，边角观测值之间精度应基本匹配，对于边角全测初始方案，也可依照边角观测值平均多余观测分量来判断边角精度匹配状况。3）依照对于网设计规定，所使用仪器，图上设计和实地踏勘，拟定观测精度和初始观测方案。观测精度应选用仪器所能达到最高精度，使优化时有减少余地；初始观测方案应对所有也许观测边和方向进行全测，故有最大多余观测数，是一种“肥网”或“密网”。4）模仿初始观测方案，进行平差计算，对精度、可靠性乃至敏捷度计算成果进行分析，一方面拟定观测精度拟定与否合理，若不合理，则需作恰当调节，在观测值精度基本合理基本上，基于观测值内部可靠性指标按从“肥”到“瘦”，从“密”到“疏”方略5进行网优化设计。详细地说，先拟定一种恰当值，再按（8）式计算n，依照差值（nn）拟定应删去观测值个数nq。对计算观测值多余观测分量按从大到小顺序排列，删去nd个多余观测分量较大那些观测值。然后重新作观测值模仿计算。如果观测值精度选取恰当，仅作一、二次迭代计算即可得到网优化设计方案。该办法特点是：初始方案方案总是一种观测精度和观测值个数均有富裕全边角网方案，如果该方案还达不到设计规定话，则阐明或者是设计规定太高，或者是所拥有仪器设备精度不够高。整个优化设计过程是如何删除多余观测和调节观测精度。按此法删除多余观测具备拟定性且不致于引起形亏（这是该法最大长处和特点）。初始方案中每测站上方向观测零方向选取要恰当，最佳选距离偏短边。对于同一测站上边和方向，还可依照通视条件和外界环境影响给定不同精度，使成果更符合实际。这些涉及优化设计中某些技巧，读者还可以在网优化设计实践中进一步总结。5 优化设计软件本文所提出优化设计办法可使用任何一种网平差软件进行，但规定软件还具备计算观测值多余观测分量和模仿观测方案功能。这里重要推荐咱们所研制“当代测量控制网数据解决通用软件包”即科傻系统系列软件之二科傻二（cosa2）6或称科达普施（Codaps）。科傻二既可单独使用，也可与科傻一联合使用（这时能自动进行改正计算、概算并生成平差输入文献），能按严密间接观测模型对一、二维网进行从一种到数千个未知点，涉及上万个观测值等权或不等权平差。除平差功能外，还可自动推算网点近似坐标，自动生成最小独立闭合环7，进行贯通误差影响值估算、坐标换算、粗差探测与定值定位、方差分量估算等8。该软件包与网优化设计关于功能是：能作网模仿计算，计算观测值多余观测分量，依照需要输出各种精度评估成果，网图显绘和生成各种报表等。下面重要简介模仿法优化设计功能。软件规定先生成一种原始观测方案文献：网名OB2，该文献构造如下：方向中误差1,边长固定误差1（mm）,比例误差1（ppm）,精度号1 点名，点类型，X坐标（m），Y坐标（m） 测站点，照准点，A，方位角值测站点点号 照准点点号，观测值类型，精度号为观测精度某些，为控制点坐标和已知方位角某些；为观测方案某些。6.1中算例1OB2文献为：1.8，3，2（，m1）（若只有一组观测精度，则仅有第一行，且m1可省去。第2行如：2.5，5，3，m2）1，0（表达已知点），5000，5000（为已知坐标）2，1（表达未知点），6000，5000（为近似坐标）1,2,A，0.0(表达1至2点已知方位角为00000.0)1（下面按顺时针方向输照准点SW，L（表达方向观测），1（与m1相应，若只有一组观测精度，可省去）SW，S（表达边长观测），1NW，L，1NW，S，12，L，12，S，1NE，L，1NE，S，1ME，L，1ME，L，1SE，L，1SE，S，12（此测站输入与第1站完全相似）程序可依照网点坐标和观测方案自动生成平差所需要输入文献（文献名为网名IN2），文献中观测值系依照坐标反换后加上模仿观测值精度。为此，在程序中，提供了生成正态原则随机数功能，对这些正态原则随机数还进行有各种检查（如正态性、周期误差、偏度、峰度、方差、均值检查等）。模仿计算得到观测值文献（网名1N2）与实测得到完全一致，因不含粗差和系统误差，故更符合最小二乘平差规定。通过平差计算网各种精度和可靠性指标，即可按前述环节进行网优化设计。6 算例和分析61 算例1-模仿桥梁边角网该网共8个点，其中1、2位于桥轴线上，北岸4个点，南岸3个点，江中岛上1个点。以1为已知点，1至2方位角为已知方位角，按独立网进行设计。优化设计规定：0.4，最弱点精度mpmax4.5mm，最弱边精度ms/s优于1/12万。初始方案为观测50个方向值，25条对向边全边角网（如图1（a）所示）。观测值精度为：mr=1.8，a=3mm，b=2ppm。初始方案计算成果：n=75，t=21，r=54，=0.72，=0.76，=0.64（边长观测值较方向观测值多余观测分量要小，阐明边长精度要高于方向精度）。点位精度均值=2.0mm，最大点位精度mpmax=2.3mm，ms/s=1/17.5万，显然质量太高。按前述环节进行如下优化设计与计算：1）取=0.5，由式（11）计算得n=42，应删去nd=75-42=33个观测值，按观测值多余观测分量从大到小排序，拟定删除rli0.80观测值22个，rsi0.66观测值11个，重新作平差计算，得：=0.53，=0.40，=3.1mm，mpmax= 3.6mm，ms/s= 1/14.2万，网质量依然偏高。2）取=0.42进行计算，得n=36，nd=39，在前述基本上，再删去3个rli0.70和3个rsi0.59共6个方向和边长观测值，计算得：=0.47，=0.29，=3.2mm，mpmax= 3.7mm，ms/s= 1/13.9万。网质量仍偏高。3）由于已达到设计极限值，为此，可依照测量仪器和作业习惯，变化观测值精度，例如，让测边精度保持不变，而mr降为2.5，仍按前面优化计算所得到观测方案进行计算，成果为：=0.5，=0.23，=3.6mm，mpmax= 4.3mm，ms/s= 1/12.1万，该成果满足设计规定，所相应为网优化设计方案（参见图1（b）。优化设计方案相对于初始方案，共减少了39个观测值，方向观测精度也从1.8减少到2.5，在精度、可靠性、经费等方面达到综合最优。图1 模仿桥梁边角网初始方案(a)和优化方案(b)阐明：变化起始点和起始方位角，例如以SW为已知点，SW-1为已知方位角，仍按网优化设计方案计算，所有可靠性指标不变，边长相对精度也不变，但点位精度变化了，从3.6mm增大到8.7mm，mpmax从4.3mm变大到12.8mm。62 算例2-某大型实测施工控制网该网总点数27个，其中已知点10个。在每个点上都设站观测，其中在21个点上进行了方向观测，6个点上只有边长观测，共有98个方向观测值，88个边长观测值。观测精度为：mr=1.0，a=1mm，b=1 ppm。按公式（6）、（7）计算得：t=321+26-102=55 r=131，=0.70，=0.69，=0.72（有多条边多余观测分量为1）。由于该网是约束网,可在相似基准下用点位精度和最弱边相对精度来衡量网质量。优化方案取=0.54，由式（11）计算得n=119，应去掉67个观测值。将计算出rli和rsi按大小顺序排列，结合经验和作业习惯，去掉31个方向观测值，36个边长观测值，作为优化方案，该方案观测精度不变，=0.54， =0.53。咱们将实测方案和优化方案进行比较成果列于表2，表中精度系依照验后单位数中误差计算。表2 某大型施工控制网实测方案和优化设计方案比较表 数 项 值 目方案基 本 情 况可 靠 性精度（mm）mmknnLnsrmpmaxmss（m）s/ms实测方案211018698881310.701.31.70.811.3214000优化方案21101196752640.541.51.90.811.32140006.3 分析 从算例1优化过程和优化成果可见，本文所提出办法既有理论根据，又简易行，优化成果很少依赖人们主观结识，一致性较好。从初始方案到优化设计方案只需作一、二次迭代计算即可。依照本文所提出算法和公式，容易编写网优化设计自动计算程序。从算例2可见，该办法可用于对既有网评价，优化方案与实测方案比较，在精度减少非常小状况下，其工作量和建网费用减小却是十分明显。相信这种状况在实际生产中还是较普遍，也阐明了工程控制网优化设计依然是一种需要研究和应用问题。7 结语工程控制网特别是施工控制网、变形监测网优化设计依然是一种值得研究和应用问题。本文在进一步研究观测值内部可靠性指标基本上，归纳了观测值多余观测分量与精度、建网费用乃至敏捷度等质量指标关系，提出了基于观测值可靠性模仿法优化设计思想和算法，指出了该办法特点和长处，并在通用平差软件基本上研制了相应优化设计程序，通过算例阐明了网优化设计过程，用该思想和办法，还可在对既有网平差同步，评价网在设计上好坏。参照文献：1佩尔策主编,张正禄编译,当代工程测量控制网理论和应用.北京:测绘出版社,19892李德仁.误差解决和可靠性理论.北京:测绘出版社,19883陈永奇等.高等应用测量.武汉:武汉测绘科技大学出版社,19964张正禄,李晓东.观测值敏捷度影响系数及其在监测网优化设计中应用.测绘学报,1989,18(4):2492575张正禄.工程平面控制网模仿法优化设计方略.武测科技,1993(1):146张正禄,黄全义等.全站式地面测量工程一体化自动化系统研究.武汉测绘科技大学学报,1999,24(1):79837冯琰,张正禄等.最小独立闭合环与附合导线自动生成算法.武汉测绘科技大学学报,1998,23(3):2552598罗年学,张正禄等.WIN95版当代测量控制网数据解决通用软件若干技术和算法.水利水电测绘科技论文集,1998,169174A new Method based on Reliability of Engineering ContralNetwork for NetsoptimizationZHANG Zhenglu，LUO Nianxue，HUANG Quangyi，MEI Wensheng，CAO Baichong (School of Surreying Scince and Technology，Wuhan University，129 Luoyu Road，Wuhan，China，430079)In this paper，after discussing 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 engineering control network which based on the inner reliability of observation is presented. It points out the advantages and characteristics of this new method. The general software packed for data processing of surveying network ( COSA) used for the method is introduced. Finally ，the calculating steps，optimization benefits and necessity are given by two examples.