CN104180822B - 一种变形监测基准点稳定性检验方法 - Google Patents

一种变形监测基准点稳定性检验方法 Download PDF

Info

Publication number
CN104180822B
CN104180822B CN201410433361.5A CN201410433361A CN104180822B CN 104180822 B CN104180822 B CN 104180822B CN 201410433361 A CN201410433361 A CN 201410433361A CN 104180822 B CN104180822 B CN 104180822B
Authority
CN
China
Prior art keywords
statistic
observation
point
stability
datum mark
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410433361.5A
Other languages
English (en)
Other versions
CN104180822A (zh
Inventor
孙泽信
汪春桃
徐春明
刘宁
庞逸群
杨树荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huashe Testing Technology Co ltd
China Design Group Co Ltd
Original Assignee
Jiangsu Provincial Communication Planning and Design Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangsu Provincial Communication Planning and Design Institute Co Ltd filed Critical Jiangsu Provincial Communication Planning and Design Institute Co Ltd
Priority to CN201410433361.5A priority Critical patent/CN104180822B/zh
Publication of CN104180822A publication Critical patent/CN104180822A/zh
Application granted granted Critical
Publication of CN104180822B publication Critical patent/CN104180822B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种变形监测基准点稳定性检验方法,结合了多点后方交会方法、测量平差方法和数理统计方法对变形监测基准点的稳定性进行检验。本发明所提供的变形监测基准点稳定性判断的检验方法具有两方面的优点:一方面稳定性判断结果的可靠性与传统方法相比有较大的提高;另一方面本方法解决了全站仪自动化观测多点后方交会实时设站基准点稳定性判断的时效性要求问题,提高了方法的适用性。

Description

一种变形监测基准点稳定性检验方法
技术领域
本发明涉及属于测绘科学或测量数据处理技术领域,涉及一种变形监测基准点稳定性检验方法。
背景技术
变形观测与分析的前提是确保监测基准点的稳定可靠,只有在有理由说明基准点稳定的前提下,以此计算变形观测结果才是可靠的,进行基准点的稳定性检验,是变形预测分析的基础。对基准点稳定性分析研究主要是通过周期性的监测网复测,在满足测量精度的前提下,求出两期结果的差异量,以分析基准点的稳定性,分析方法包括限差分析法、t检验法、平均间隙法、单点位移分量法等。在实际的工程应用中,以限差检验法和t检验法两种方法应用较为广泛。
随着高精度全站仪的推广,全站仪自动化观测也得到了较好的应用,多点后方交会以其特有的优势得到了广泛的应用,可以极大地提高全站仪自动化变形观测的时效性及观测精度,同时后方交会基准点的稳定性也直接影响后方交会设站精度及监测精度,而传统的基准点稳定性检验方法无法满足这一时效性及多点后方交会设站模式。
发明内容
发明目的:本发明为了克服现有技术中存在的不足,本发明提供一种有效提高全站仪自动化观测多点后方交会设站基准点的稳定性、判断结果的可靠性及稳定性分析的时效性的变形监测基准点稳定性检验方法。
发明内容:为解决上述技术问题,本发明提供了一种变形监测基准点稳定性检验方法,包括以下步骤:
步骤1:通过全站仪自由设站法观测采集数据,全站仪架设在工作基点上,观测多个基准点,观测值为工作基点到基准点的水平距离和相邻两条线之间的夹角;
步骤2:利用后方交会原理建立多点后方交会模型;
步骤3:结合步骤1中采集的数据,利用步骤2中建立的多点后方交会模型计算工作基点坐标;
步骤4:计算观测值后验单位权中误差
步骤5:采用间接平差方法计算所有观测值的修正值v;
步骤6:进行基准网的稳定性判断:对观测值后验单位权中误差构造卡方检验统计量;设定卡方检验统计量的置信水平,结合卡方检验统计量的置信水平通过查表找出分位值,比较统计量与对应分位值的大小;如果统计量小于分位值,表明基准网稳定,结束基准网稳定性判断;若统计量大于分位值,进行步骤7;
步骤7:进行基准点稳定性判断:对工作基点到基准点的水平距离的观测值修正值进行t统计量构造,设置t统计量构造的置信水平,结合t统计量构造的置信水平通过查表找出分位值,比较每个统计量与分位值之间的大小,大于分位值的统计量对应的基准点为不稳定基准点,剔除不稳定点后在执行步骤3~步骤6,直至基准网符合稳定性要求。
进一步,所述步骤6中卡方检验统计量χ的计算方法为:其中,为观测值后验单位权中误差,σ0为单位权中误差,所述σ0取值为夹角度观测值中误差,f为自由度。
进一步,所述步骤7中工作基点到基准点的水平距离的观测值修正值进行t统计量构造的方法为:根据公式计算每个基准点到工作基点的水平距离的观测值修正值的t统计量,其中,表示第i个基准点到工作基点的水平距离的观测值修正值的t统计量,为第i个基准点到工作基点的水平距离的观测值,为所有基准点到工作基点的水平距离的观测值的平均值,n为基准点的个数,S为各水平距离观测值修正值的样本标准差, S 2 = Σ i = 1 n ( v s i - v ‾ s ) 2 n - 1 .
进一步,所述步骤1中基准点的个数大于5个。
进一步,所述步骤6中卡方检验统计量的置信水平为95%,所述步骤7中t统计量构造的置信水平为99%。
有益效果:与现有技术相比,本发明所提供的变形监测基准点稳定性判断的检验方法具有两方面的优点:一方面稳定性判断结果的可靠性与传统方法相比有较大的提高;另一方面本方法解决了全站仪自动化观测多点后方交会实时设站基准点稳定性判断的时效性要求问题,提高了方法的适用性。
附图说明
图1为本发明的的流程图。
具体实施方式
下面结合附图对本发明的技术方案作进一步解释。
如图1所示,一种变形监测基准点稳定性检验方法,包括以下步骤:
步骤1:通过全站仪自由设站法观测采集数据,全站仪架设在工作基点上,观测由7个基准点组成的基准网,观测值为工作基点到每个基准点的水平距离和相邻两条线之间的夹角;
步骤2:利用后方交会原理建立多点后方交会模型;
步骤3:结合步骤1中采集的数据,利用步骤2中建立的多点后方交会模型计算工作基点坐标;
步骤4:计算观测值后验单位权中误差
步骤5:采用间接平差方法计算所有观测值的修正值v;
步骤6:进行基准网的稳定性判断:对观测值后验单位权中误差构造卡方检验统计量设定卡方检验统计量的置信水平为95%,结合卡方检验统计量的置信水平计算分位值,比较统计量与对应分位值的大小;如果统计量小于分位值,表明基准网稳定,结束基准网稳定性判断;若统计量大于分位值,进行步骤7;
步骤7:进行基准点稳定性判断:对工作基点到基准点的水平距离的观测值修正值进行t统计量构造,设置t统计量构造的置信水平为99%,结合t统计量构造的置信水平计算分位值,比较每个统计量与分位值之间的大小,大于分位值的统计量对应的基准点为不稳定基准点,剔除不稳定点后在执行步骤3~步骤6,直至基准网符合稳定性要求。
实施例:
以“地铁保护区自动化监测”为应用实例,对本发明进一步阐明:
网中有1个工作基点,7个基准点作为后方交会点,7个基准点分别命名为H1、H2、H3、H4、H5、H6、H7。观测仪器采用Leica TM30型全站仪进行观测,仪器标称精度为±0.5″,±(0.6mm+10-6·D)。使用本文所述的基准点稳定性检验法进行基准点稳定性分析。各基准点起算数据如表1所示,其中P点坐标为工作基点坐标的近似值。后方交会观测值成果如表2所示。
表1
表2
表2中S1~S7为基准点H1~H7分别与工作基点构成的直线的名称,其观测值为工作基点到每个基准点的水平距离;α1表示直线S1与S2之间的夹角,即α表示相邻两条线之间的夹角。
根据利用后方交会原理推导的平差模型对观测值进行平差,具体平差成果如表3所示。基准网稳定性分析如表4所示。
表3
基准网稳定性分析如表4所示,选定置信水平α=95%,自由度为11时经查表得到卡方检验的临界值为19.675,即分位值为19.675,所构造统计量为52.360。其中,观测值个数=13,未知数=2个,所以,自由度f=13-2=11,未知数为工作基点的平面坐标,由于分位值小于构造统计量,故可判断基准网中存在显著不稳定点。
表4
从表4基准网稳定性分析表中可知基准网中存在不稳定基准点,需进行基准点稳定性检验,具体如表5所示,选定置信水平α=99%,自由度为n-1=6,经查表得统计量的临界值为t0.01,6=3.143,即分位值为3.143,按照上述方法计算的观测值统计量中,边长S4对应的t统计量为5.457,因此,认为基准点H4是不稳定的,将其剔除,对剩余的基准点重新进行后方交会平差,并再次判断基准网的稳定性。平差结果如表6所示,基准网稳定性分析如表7所示。
表5
表6
基准网稳定性分析如表7所示,选定置信水平α=95%,自由度为9时经查表得到卡方检验的临界值为16.916,所构造统计量为10.800,故可判断基准网中不存在显著不稳定点,其中,观测值个数=11,未知数=2个,自由度f=11-2=9,未知数为工作基点的平面坐标。
表7
从表7中可以看出,在剔除基准点H4后,后方交会平差模型稳定状态良好,不存在显著不稳定点,平差成果可以接受。同时对比不稳定点剔除前后的平差成果表即表3与表6,可以看出剔除不稳定点后,后验单位权中误差明显变小,进一步表明基准点稳定性检验方法具有较高的可靠性。

Claims (4)

1.一种变形监测基准点稳定性检验方法,其特征在于:包括以下步骤:
步骤1:通过全站仪自由设站法观测采集数据,全站仪架设在工作基点上,观测多个基准点,观测值为工作基点到基准点的水平距离和相邻两条线之间的夹角;
步骤2:利用后方交会原理建立多点后方交会模型;
步骤3:结合步骤1中采集的数据,利用步骤2中建立的多点后方交会模型计算工作基点坐标;
步骤4:计算观测值后验单位权中误差
步骤5:采用间接平差方法计算所有观测值的修正值v;
步骤6:进行基准网的稳定性判断:对观测值后验单位权中误差构造卡方检验统计量;设定卡方检验统计量的置信水平,结合卡方检验统计量的置信水平通过查表找出分位值,比较统计量与对应分位值的大小;如果统计量小于分位值,表明基准网稳定,结束基准网稳定性判断;若统计量大于分位值,进行步骤7;
步骤7:进行基准点稳定性判断:对工作基点到基准点的水平距离的观测值修正值进行t统计量构造,设置t统计量构造的置信水平,结合t统计量构造的置信水平通过查表找出分位值,比较每个统计量与分位值之间的大小,大于分位值的统计量对应的基准点为不稳定基准点,剔除不稳定点后在执行步骤3~步骤6,直至基准网符合稳定性要求;
所述步骤7中工作基点到基准点的水平距离的观测值修正值进行t统计量构造的方法为:根据公式计算每个基准点到工作基点的水平距离的观测值修正值的t统计量,其中,表示第i个基准点到工作基点的水平距离的观测值修正值的t统计量,为第i个基准点到工作基点的水平距离的观测值,为所有基准点到工作基点的水平距离的观测值的平均值,n为基准点的个数,S为各水平距离观测值修正值的样本标准差,
2.根据权利要求1所述的变形监测基准点稳定性检验方法,其特征在于:所述步骤6中卡方检验统计量χ的计算方法为:其中,为观测值后验单位权中误差,σ0为单位权中误差,所述σ0取值为夹角度观测值中误差,f为自由度。
3.根据权利要求1所述的变形监测基准点稳定性检验方法,其特征在于:所述步骤1中基准点的个数大于5个。
4.根据权利要求1所述的变形监测基准点稳定性检验方法,其特征在于:所述步骤6中卡方检验统计量的置信水平为95%,所述步骤7中t统计量构造的置信水平为99%。
CN201410433361.5A 2014-08-28 2014-08-28 一种变形监测基准点稳定性检验方法 Active CN104180822B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410433361.5A CN104180822B (zh) 2014-08-28 2014-08-28 一种变形监测基准点稳定性检验方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410433361.5A CN104180822B (zh) 2014-08-28 2014-08-28 一种变形监测基准点稳定性检验方法

Publications (2)

Publication Number Publication Date
CN104180822A CN104180822A (zh) 2014-12-03
CN104180822B true CN104180822B (zh) 2016-09-07

Family

ID=51962056

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410433361.5A Active CN104180822B (zh) 2014-08-28 2014-08-28 一种变形监测基准点稳定性检验方法

Country Status (1)

Country Link
CN (1) CN104180822B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108824499B (zh) * 2018-06-14 2020-11-17 合肥市测绘设计研究院 一种基于后方交会自由设站水平位移监测方法
CN110260840B (zh) * 2019-06-28 2021-06-29 中铁大桥局集团有限公司 跨海长桥桥墩绝对沉降量的测量方法及系统
CN111044022B (zh) * 2019-12-31 2020-10-02 中国科学院国家天文台 一种基于全站仪互瞄技术的基准网自动化测量方法及系统
CN113449254B (zh) * 2021-06-09 2023-09-01 华设设计集团股份有限公司 任意网型变形监测稳定性分析方法及监控点位置确定方法
CN114923459A (zh) * 2022-04-07 2022-08-19 江苏易政达信息科技有限公司 一种基于静力水准仪的沉降监测系统及测量和标定方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06331365A (ja) * 1993-05-27 1994-12-02 Hitachi Cable Ltd 回転角速度センサの校正方法
US7805269B2 (en) * 2004-11-12 2010-09-28 Philips Electronics Ltd Device and method for ensuring the accuracy of a tracking device in a volume
CN102278970B (zh) * 2011-06-14 2013-03-20 北京林业大学 一种基于全站仪角距差分的定位变形监测技术
CN103090830A (zh) * 2011-11-03 2013-05-08 中国石油天然气股份有限公司 冻土区油气管道位移监测方法及装置
CN102506812B (zh) * 2011-11-18 2014-04-23 南京市测绘勘察研究院有限公司 一种变形监测基准点稳定性判断的vt检验方法
CN103292762B (zh) * 2013-05-23 2016-08-03 青岛理工大学 判别大坝稳定性的位移监测方法
CN103352455B (zh) * 2013-07-30 2015-07-15 上海城建市政工程(集团)有限公司 开挖基坑坑底土体隆起的监测方法
CN103940364B (zh) * 2014-05-04 2016-06-22 赵鸣 地铁隧道相对变形的摄影测量方法

Also Published As

Publication number Publication date
CN104180822A (zh) 2014-12-03

Similar Documents

Publication Publication Date Title
CN104180822B (zh) 一种变形监测基准点稳定性检验方法
US7720647B2 (en) System and method for filtering a point cloud
CN105740203B (zh) 多传感器无源协同测向定位方法
CN102494940B (zh) 基于Beremin模型的材料脆性断裂评定参量的标定方法
CN107239659B (zh) 一种基于改进的k-rusle模型与成土速率计算土壤侵蚀的方法
CN104776827B (zh) Gps高程异常数据的粗差探测方法
CN104218570B (zh) 一种电能计量装置整体计量误差的在线评估方法和系统
CN103759732B (zh) 一种角度信息辅助的集中式多传感器多假设跟踪方法
CN109490072B (zh) 一种土木工程建筑用检测系统及其检测方法
CN104749251B (zh) 一种剔除水声距影响的晶粒尺寸超声评价方法
CN104280013A (zh) 基于测量坐标确定岩体结构面产状的方法
CN105484735B (zh) 一种实钻井眼轨迹与设计轨道符合率的评价方法
CN115327674B (zh) 一种大规模降雨监测雷达校验方法
CN103438872A (zh) 一种基于大坝三维前方交会测量的内外业一体化系统
CN103669319A (zh) 用于确定环境边坡危岩体稳定性的快速评价施工方法
CN102506812B (zh) 一种变形监测基准点稳定性判断的vt检验方法
CN105043390B (zh) 基于泛克里金法的重力场插值方法
CN103292773A (zh) 一种基于最小区域的对称度误差评定方法
CN103984017B (zh) 一种微震震源自动定位方法
CN103678933B (zh) 用于确定环境边坡危岩体稳定性的矩阵评价施工方法
CN105204047B (zh) 一种卫星导航系统中观测量单个粗差的探测与修复方法
CN109523129A (zh) 一种无人车多传感器信息实时融合的方法
CN104751005B (zh) 一种基于正交实验的平面度误差评定方法
CN104502889B (zh) 指纹定位中基于参考点最大距离的定位可信度计算方法
CN106558086B (zh) 一种林木应力波断层成像的速度误差校正方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP03 Change of name, title or address

Address after: 210014 Ziyun Road, Qinhuai District, Jiangsu, No. 9, No.

Patentee after: CHINA DESIGN GROUP Co.,Ltd.

Address before: 210000 Ziyun Road, Qinhuai District, Jiangsu, No. 9, No.

Patentee before: Jiangsu Provincial Communications Planning and Design Institute Co.,Ltd.

CP03 Change of name, title or address
CP01 Change in the name or title of a patent holder

Address after: 210014 No. 9, Ziyun Avenue, Qinhuai District, Jiangsu, Nanjing

Patentee after: China Design Group Co.,Ltd.

Address before: 210014 No. 9, Ziyun Avenue, Qinhuai District, Jiangsu, Nanjing

Patentee before: CHINA DESIGN GROUP Co.,Ltd.

CP01 Change in the name or title of a patent holder
TR01 Transfer of patent right

Effective date of registration: 20220429

Address after: 210000 9 Ziyun Avenue, Qinhuai District, Nanjing, Jiangsu Province

Patentee after: China Design Group Co.,Ltd.

Patentee after: Huashe Testing Technology Co.,Ltd.

Address before: 210014 No. 9 Ziyun Avenue, Qinhuai District, Nanjing City, Jiangsu Province

Patentee before: China Design Group Co.,Ltd.

TR01 Transfer of patent right