CN112731545B - 一种海洋不规则测网数据平差方法 - Google Patents
一种海洋不规则测网数据平差方法 Download PDFInfo
- Publication number
- CN112731545B CN112731545B CN202011509770.0A CN202011509770A CN112731545B CN 112731545 B CN112731545 B CN 112731545B CN 202011509770 A CN202011509770 A CN 202011509770A CN 112731545 B CN112731545 B CN 112731545B
- Authority
- CN
- China
- Prior art keywords
- adjustment
- line
- survey
- supplementary
- measuring
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
Abstract
本发明涉及一种海洋不规则测网数据的平差方法,属于船载重力仪观测数据系统差校准及数据处理的应用领域,所述方法基于补充测线与主测线、联络线交点分布及大小等因素,建立6种平差模型,在此基础上,选取测线精度、交点差均方根及平差结果整体均一性在内的三个指标进行分析评估确定最佳平差方式。本发明方法在进行平差时,考虑到补充测线与主测线和联络测线斜交产生的交点,避免因舍弃这部分交点而造成平差精度的降低,同时也解决了直接遍历全部测线而造成交点差重复计算的问题。因此本方法充分利用了现场观测资料并防止在平差过程中进一步造成新的系统误差。
Description
技术领域
本发明属于船载重力仪观测数据平差处理及精度分析的应用领域,具体地涉及一种海洋不规则测网数据平差方法。
背景技术
海洋重力测量是指在海洋范围内进行的重力加速度的测定工作,为研究地球形状和地球内部构造、勘探海洋矿产资源、保障航天和远程武器发射等提供重力资料。海洋重力测量与陆地重力测量具有相同的研究目的和作用,是地球重力测量学研究的重要组成部分。海洋调查过程中,受仪器设备、海况、定位精度等多种因素影响,重力观测值内包含一定大小的偶然误差和系统误差。海洋重力测网平差作为内业数据处理、质量控制的最后一步,其关键性不言而喻。传统的测网仅依据主测线和联络线的交点个数进行划分,以m条主测线、n条联络线组成的测网为例,具有m×n个交点时为规则测网,因走向、长度差异造成交点数少于m×n个交点的测网则为不规则测网。该分类方法并未考虑随着同一区域内多个航次的往复测量、国际共享资料的逐年累加,乃至航次内的补充测线等问题。
发明内容
本发明提供了一种海洋不规则测网数据的平差方法,所述方法基于补充测线与主测线、联络线交点分布及大小等因素,建立6种平差模型,在此基础上,选取测线精度、交点差均方根及平差结果整体均一性在内的三个指标进行分析评估确定最佳平差方式。
本发明是通过如下技术方案来实现的:
一种海洋不规则测网数据平差方法,所述方法包括以下步骤:
1)按照测线航迹分布、航速在内的因素,剔除弯折段、停航段及无交点在内的无效数据,将测线分为主测线、联络线、补充测线三类;补充测网为实际开展作业过程中受现场实际因素制约影响,产生的与主测线、联络线斜交的测线;
2)针对该补充测线数据采用6种方式进行分析:①构建系数矩阵时,先暂时舍弃补充测线与主测线的交点,即将补充测线设定为主测线进行平差计算;②构建系数矩阵时,先暂时舍弃补充测线与联络线的交点,即将补充测线设定为联络线进行平差计算;③以与主测线、联络线何种类型测线交点个数少且误差相对更小则将补充测线设定为何种测线为原则,逐条设定补充测线的类型并平差;④将全部补充测线视为一个系统,经平差后的规则测网视为另一个高精度系统,将两个系统做差,将补充测线系统调差至已平差的高精度规则测网系统中;⑤经方式④处理后,再对全部平差结果测线遍历进行平差;⑥遍历全部原始测线完成数据平差;
3)针对上述6种类型平差方式,采用测线精度、交点差均方根及平差结果整体均一性在内的三个指标选定最终平差方式。
进一步,所述的步骤3)选定平差方式的具体方法:首先对比各平差方式的均方根,舍弃较大均方根的平差方式,然后针对平差结果一致的情况,依据测网系统精度去除测网系统精度较低的平差方式,再结合平差后交点差的均一性确定最终平差方式。
本发明与现有技术相比的有益效果:
现有平差方法只是针对主测线和联络线开展测网数据平差工作,并未考虑随着同一区域内多个航次的往复测量、国际共享资料的逐年累加,乃至航次内的补充测线等问题。本发明方法在进行平差时,引入了补充测线概念,并详细考虑了该类型测线与主测线和联络测线斜交产生的交点,将实际调查测网数据划分为6种基础平差模型;避免因舍弃这部分交点而造成平差精度的降低,同时也解决了直接遍历全部测线的方式而造成的交点差重复计算的问题。综上所述,本方法充分利用了现场观测资料并防止进一步造成新的系统误差,可以广泛应用于不规则复杂测网下测线数据的平差工作中。
附图说明
图1本发明方法的流程图;
图2MW9006测线分割图;
图3为各种方式平差结果图:(a)平差前交点误差分布图,(b)方式①交点误差分布图,(c)方式②和③交点误差分布图,(d)方式④交点误差分布图, (e)方式⑤和⑥交点误差分布图。
图4方式⑤和⑥平差后交点差对比图。
具体实施方式
下面通过实施例来对本发明的技术方案做进一步解释,但本发明的保护范围不受实施例任何形式上的限制。
实施例1
1)以共享航次MW9006为例。方法流程如图1所示,按照测线航迹分布、航速等因素,剔除弯折段、停航段及无交点等无效数据,将测线分为50条主测线(M),28条联络线(L),3条补充测线(S)三类(图2)。
2)针对该补充测线数据可采用4种方式进行分析:①构建系数矩阵时,先暂时舍弃补充测线与主测线的交点,即将补充测线设定为主测线进行平差计算;②构建系数矩阵时,先暂时舍弃补充测线与联络线的交点,即将补充测线设定为联络线进行平差计算;③以与主测线、联络线何种类型测线交点个数少且误差相对更小则将补充测线设定为何种测线为原则,逐条设定补充测线的类型并平差;④将全部补充测线视为一个系统,经平差后的规则测网视为另一个高精度系统,将两个系统做差,将补充测线系统调差至已平差的高精度规则测网系统中;⑤经方式④处理后,再对全部平差结果测线遍历进行平差;⑥遍历全部原始测线完成数据平差;
具体平差过程为:
海洋重力测网交点处的主测线(gij)、联络线(gji)交点观测值可表示为:
gij=gtij+εij+Δij (1)
gji=gtji+βji+Δ′ij (2)
式中i,j为测线编号,gt为交点真实值,εij、βji分别为第i条主测线和第j 条联络线的系统误差,△ij为随机误差,交点差为:
dij=gij-gji=εij-βji+Δij-Δ′ij (3)
数据处理过程中将各个测点视为等精度观测;采用半系统差调整方式计算时,分别对主测线和联络线进行一次半系统差调整后,交点差则变为:
迭代计算至交点差代数和为零。
最终在计算数据精度时以均方根作为评价参数,其计算公式为:
6种类型的计算结果为:
方式1:均方根=±2.14×10-5m/s2
方式2:均方根=±1.84×10-5m/s2
方式3:均方根=±1.84×10-5m/s2
方式4:均方根=±2.01×10-5m/s2
方式5:均方根=±1.81×10-5m/s2
方式6:均方根=±1.81×10-5m/s2
经对比,其中方式⑤和⑥经平差后,测线交点差相等,但并非每条测线的平差值均相同,方式⑥测线平差值普遍比方式⑤大,经计算方式⑤的调差值总和为-18.86,方式⑥的调差值总和为71.42。分析认为主要原因为方式⑤判断原主测线、联络线组成的测网系统精度更高,在完成补偿测线的独立系统平差时,先完成了整体调差(消除相对系统误差)。因此方式⑤和⑥判断依据是如果主测线、联络线的测网系统精度更高,采用方式⑤,如果补充测线精度更高,采用方式⑥。
本实施例在对6种结果进行评价时,选取了测线精度、交点差均方根及平差结果整体均一性等三个要素进行分析评估。6种方式的交点误差分布如图3和图4所示。
经对比:
1、方式①和④的交点差均方根较大,代表平差结果较差,所以不能采用。
2、方式⑤和⑥的交点均方根一致,但依据现场实例分析发现方式⑤的主测线、联络线的测网系统精度更高,所以不能采用方式⑥;
3、本实例属于方式②和③平差结果一致的特殊实例,故将两种类型合并为 1种;
4、比对方式②、③和⑤,虽然平差后均方根仅相差0.03×10-5m/s2,但采用方式②、③所获得结果集中在-8到8,而方式⑤则集中在-10到10,所以方式②、③所获得的结果均一性更佳。
综上所述,最终确定实例中方式②和③为最佳的平差方式,最终航次内114 个交点均方根降为±1.84×10-5m/s2。
传统计算时,或者直接舍弃补充测线开展平差,往往只能导致损失部分交点数据不参与平差而降低精度;或者直接遍历测线的方式计算全部交点并参与平差,但存在导致系统误差的问题。本实施例克服了上述两个缺点,充分利用了现场观测资料并防止进一步造成新的系统误差。
Claims (1)
1.一种海洋不规则测网数据平差方法,其特征在于所述方法包括以下步骤:
1)按照测线航迹分布、航速在内的因素,剔除弯折段、停航段及无交点在内的无效数据,将测线分为主测线、联络线、补充测线三类;补充测网为实际开展作业过程中受现场实际因素制约影响,产生的与主测线、联络线斜交的测线;
2)针对该补充测线数据采用6种方式进行分析:①构建系数矩阵时,先暂时舍弃补充测线与主测线的交点,即将补充测线设定为主测线进行平差计算;②构建系数矩阵时,先暂时舍弃补充测线与联络线的交点,即将补充测线设定为联络线进行平差计算;③以与主测线、联络线何种类型测线交点个数少且误差小则将补充测线设定为何种测线为原则,逐条设定补充测线的类型并平差;④将全部补充测线视为一个系统,经平差后的规则测网视为另一个高精度系统,将两个系统做差,将补充测线系统调差至已平差的高精度规则测网系统中;⑤经方式④处理后,再对全部平差结果测线遍历进行平差;⑥遍历全部原始测线完成数据平差;
3)针对上述6种类型的平差方式,采用测线精度、交点差均方根及平差结果整体均一性在内的三个指标选定最终平差方式;
所述选定最终平差方式的具体方法:首先对比各平差方式的均方根,舍弃大均方根的平差方式,然后针对平差结果一致的情况,依据测网系统精度去除测网系统精度低的平差方式,再结合平差后交点差的均一性确定最终平差方式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011509770.0A CN112731545B (zh) | 2020-12-18 | 2020-12-18 | 一种海洋不规则测网数据平差方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011509770.0A CN112731545B (zh) | 2020-12-18 | 2020-12-18 | 一种海洋不规则测网数据平差方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112731545A CN112731545A (zh) | 2021-04-30 |
CN112731545B true CN112731545B (zh) | 2022-07-22 |
Family
ID=75603938
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011509770.0A Active CN112731545B (zh) | 2020-12-18 | 2020-12-18 | 一种海洋不规则测网数据平差方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112731545B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102168806A (zh) * | 2011-01-30 | 2011-08-31 | 中国海洋石油总公司 | 一种天然气长输管道位置信息采集方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10474767B2 (en) * | 2016-01-26 | 2019-11-12 | Saudi Arabian Oil Company | Gravity modeling a rifted continental margin |
CN108828669B (zh) * | 2018-04-02 | 2020-07-10 | 中国石油天然气集团有限公司 | 一种二维相交测线静校正处理方法、装置及系统 |
CN109870733A (zh) * | 2019-03-08 | 2019-06-11 | 中国石油天然气集团有限公司 | 一种陆地勘探数据的测量装置以及勘探数据的处理方法 |
CN110927823B (zh) * | 2019-12-04 | 2021-03-16 | 自然资源部第一海洋研究所 | 一种应用于船载重力仪的格值校正方法 |
-
2020
- 2020-12-18 CN CN202011509770.0A patent/CN112731545B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102168806A (zh) * | 2011-01-30 | 2011-08-31 | 中国海洋石油总公司 | 一种天然气长输管道位置信息采集方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112731545A (zh) | 2021-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109141426B (zh) | 一种水下重力匹配导航适配区的方法 | |
CN106126961B (zh) | 灰色粒子群卫星钟差自适应预测方法及系统 | |
CN106646645B (zh) | 一种重力正演加速方法 | |
CN106354942B (zh) | 一种基于熵理论的装配接触应力分布评估方法 | |
CN105654483A (zh) | 三维点云全自动配准方法 | |
CN108344409B (zh) | 提高卫星姿态确定精度的方法 | |
CN109883406B (zh) | 基于更少点的全站仪概略位姿的计算方法及系统 | |
CN109446726A (zh) | 基于大数据分析的获取边坡变形三维数据的边坡监测系统 | |
CN108846200B (zh) | 一种基于迭代法的准静态桥梁影响线识别方法 | |
CN112950008A (zh) | 一种飞行器的航迹规划评价方法 | |
CN111339615A (zh) | 一种基于两步插值的飞机油箱重心计算与补偿方法 | |
CN109031224A (zh) | 一种雷达系统可靠性模糊分配方法 | |
CN113486591B (zh) | 一种卷积神经网络结果的重力多参量数据密度加权反演方法 | |
CN112731545B (zh) | 一种海洋不规则测网数据平差方法 | |
CN115342814A (zh) | 一种基于多传感器数据融合的无人船定位方法 | |
CN113419280B (zh) | 基于改进椭圆拟合的叠前裂缝密度估算方法 | |
CN112131752B (zh) | 一种基于拟准检定的超强崩溃污染率抗差估计算法 | |
CN113221278A (zh) | 一种车载光电平台照准架轻量化方法 | |
CN104517121A (zh) | 一种基于粒子群优化的空间大数据字典学习方法 | |
CN108595792B (zh) | 一种飞机油箱燃油测量传感器布局优化方法 | |
CN109633781B (zh) | 地质属性获取方法及装置、电子设备、存储介质 | |
CN108304649B (zh) | 一种高层建筑形变预测方法 | |
CN108731616B (zh) | 一种基于云模型的螺旋锥齿轮齿面测量点自适应分布方法 | |
CN114993604B (zh) | 一种基于深度学习的风洞天平静态校测方法 | |
CN111768435B (zh) | 一种应用于零件自动找正的自适应步长点云配正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |