CN114137625A - 一种基于两套惯导互观测的海洋垂线偏差测量方法 - Google Patents

一种基于两套惯导互观测的海洋垂线偏差测量方法 Download PDF

Info

Publication number
CN114137625A
CN114137625A CN202111404563.3A CN202111404563A CN114137625A CN 114137625 A CN114137625 A CN 114137625A CN 202111404563 A CN202111404563 A CN 202111404563A CN 114137625 A CN114137625 A CN 114137625A
Authority
CN
China
Prior art keywords
information
sets
inertial navigation
latitude
longitude
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.)
Granted
Application number
CN202111404563.3A
Other languages
English (en)
Other versions
CN114137625B (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.)
707th Research Institute of CSIC
Original Assignee
707th Research Institute of CSIC
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 707th Research Institute of CSIC filed Critical 707th Research Institute of CSIC
Priority to CN202111404563.3A priority Critical patent/CN114137625B/zh
Publication of CN114137625A publication Critical patent/CN114137625A/zh
Application granted granted Critical
Publication of CN114137625B publication Critical patent/CN114137625B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • G01V7/02—Details
    • G01V7/06—Analysis or interpretation of gravimetric records
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Navigation (AREA)

Abstract

本发明涉及一种基于两套惯导互观测的海洋垂线偏差测量方法，其技术特点是：在垂线偏差基准点，校准两套惯导位置输出；进入测线前以卫导位置为基准观测惯导经纬度误差，获取经纬度误差中地球周期振荡项的相位；L测线航行期间以卫导位置为基准获得惯导经纬度误差，去除惯导自身误差源引起的误差项后，得到L测线东西向、南北向垂线偏差测量值。本发明设计合理，在惯导经纬度误差中分离出垂线偏差，给出了完整的垂线偏差测量方案，具有精度高、实用性强等特点。

Description

一种基于两套惯导互观测的海洋垂线偏差测量方法
技术领域
本发明海洋重力信息测量技术领域，涉及海洋平台的垂线偏差测量，尤其是一种基于两套惯导互观测的海洋垂线偏差测量方法。
背景技术
海洋是国防屏障、资源宝藏和重要通道，与国家安全与权益、可持续发展以及战略性资源保障息息相关。海洋重力场信息是一种重要的海洋环境信息，重力信息包括重力异常和垂线偏差。目前重力异常测量方法已经趋于成熟，但是尚不具备海洋垂线偏差测量能力。
发明内容
本发明的目的在于填补现有技术的空白，提供一种设计合理且精度高的基于两套惯导互观测的水下平台海洋垂线偏差测量方法。
本发明解决现有的技术问题是采取以下技术方案实现的：
一种基于两套惯导互观测的海洋垂线偏差测量方法，包括以下步骤：
步骤1、载体航行至垂线偏差基准点，将两套惯导位置输出校准为基准点准确位置，记录基准点东西方向、南北方向垂线偏差值ge0、gn0；
步骤2、沿东西方向设计垂线偏差测线，载体进入测线前在测线同纬度且垂线偏差变化平缓区域游弋一段时间，以卫导经纬度为基准观测两套惯导经纬度误差，提取两套惯导纬度误差中地球周期振荡相位
Figure BDA0003371841520000011
经度误差中地球周期振荡相位
Figure BDA0003371841520000012
步骤3、载体进入垂线偏差测线L航行，期间以卫导经纬度为基准，得到两套惯导的经度误差lonE1(t)、lonE2(t)和纬度误差latE1(t)、latE2(t)，基于两套惯导经度差值lon12(t)、纬度差值lat12(t)估计两套惯导自身误差源引起的经度误差项lonE1a(t)、lonE2a(t)和纬度误差项latE1a(t)、latE2a(t)，则L测线南北向垂线偏差为：
gn(t)＝latE1(t)-latEla(t)+gn0或gn(t)＝latE2(t)-latE2a(t)+gn0
东西向垂线偏差为：
Figure BDA0003371841520000013
其中
Figure BDA0003371841520000014
为纬度。
而且，所述步骤3中基于两套惯导经度差值lon12(t)估计2套惯导自身误差源引起的经度误差项lonE1a(t)、lonE2a(t)的具体方法为：
(1)提取lon12(t)中地球周期振荡项幅值A12e、相位
Figure BDA0003371841520000015
依据
Figure BDA0003371841520000016
和A12e、
Figure BDA0003371841520000017
构建三角形求得两套惯导经度误差中地球周期振荡项幅值A1e、A2e；
(2)Ω为地球周期对应的角频率，计算
Figure BDA0003371841520000018
Figure BDA0003371841520000021
Figure BDA0003371841520000022
将lon12R(t)中信息与lonE1R(t)、-lonE2R(t)中信息进行特征比对，lon12R(t)中与lonE1R(t)中信息特征符合的特征信息为D1(t)，与-lonE2R(t)中信息特征符合的特征信息为D2(t)，则
Figure BDA0003371841520000023
而且，所述步骤3中基于2套惯导纬度差值lat12(t)估计2套惯导自身误差源引起的纬度误差项latE1a(t)、latE2a(t)的具体方法为：
(1)提取lat12(t)中地球周期振荡项幅值A12n、相位
Figure BDA0003371841520000024
依据
Figure BDA0003371841520000025
和A12n、
Figure BDA0003371841520000026
构建三角形求得两套惯导纬度误差中地球周期振荡项幅值A1n、A2n；
(2)Ω为地球周期对应的角频率，计算
Figure BDA0003371841520000027
Figure BDA0003371841520000028
Figure BDA0003371841520000029
将lat12R(t)中信息与latE1R(t)、-latE2R(t)中信息进行特征比对，lat12R(t)中与latE1R(t)中信息特征符合的特征信息为F1(t)，与-latE2R(t)中信息特征符合的特征信息为F2(t)，则
Figure BDA00033718415200000210
本发明的优点和积极效果是：
本发明设计合理，利用两套惯导经纬度差值观测还原惯导经纬度误差中的地球周期振荡项，并通过特征匹配识别惯导经纬度误差中其它由惯导自身因素引起的部分，支撑实现从惯导经纬度误差中分离出垂线偏差，给出了完整的垂线偏差测量方案，具有精度高、实用性强等特点。
附图说明
图1为本发明的基于两套惯导经纬度差值中地球周期振荡项得到两套惯导各自经纬度误差中地球周期振荡幅值的示意图。
具体实施方式
以下结合附图对本发明实施例做进一步详述。
本发明是在两套中高精度惯导和一套卫导组成海洋垂线偏差测量系统上实现的。两套惯导工作于阻尼状态，其阻尼网络参数设置一致，以保证两套惯导对垂线偏差输出频率响应的一致性。
基于上述说明，本发明提出一种基于两套惯导互观测的海洋垂线偏差测量方法，包括以下步骤：
一种基于两套惯导互观测的海洋垂线偏差测量方法，包括以下步骤：
步骤1、载体航行至垂线偏差基准点，将两套惯导位置输出校准为基准点准确位置，记录基准点东西方向、南北方向垂线偏差值ge0、gn0。
步骤2、沿东西方向设计垂线偏差测线，载体进入测线前在测线同纬度且垂线偏差变化平缓区域游弋一段时间，以卫导经纬度为基准观测两套惯导经纬度误差，提取两套惯导纬度误差中地球周期振荡相位
Figure BDA0003371841520000031
经度误差中地球周期振荡相位
Figure BDA0003371841520000032
步骤3、载体进入垂线偏差测线L航行，期间以卫导经纬度为基准，得到两套惯导的经度误差lonE1(t)、lonE2(t)和纬度误差latE1(t)、latE2(t)，则L测线南北向、东西向垂线偏差具体求解过程为：
(1)惯导经度误差中包含垂线偏差引起的部分和惯导自身误差源引起的部分，求取垂线偏差必须分离这两部分。两套惯导经度差值lon12(t)中不包含垂线偏差引起的部分，主要是两套惯导自身误差源引起部分的合成，因而基于lon12(t)可估计2套惯导自身误差源引起的经度误差项lonE1a(t)、lonE2a(t)，具体方法为：
①提取lon12(t)中地球周期振荡项幅值A12e、相位
Figure BDA0003371841520000033
依据
Figure BDA0003371841520000034
和A12e、
Figure BDA0003371841520000035
构建三角形求得两套惯导经度误差中地球周期振荡项幅值A1e、A2e。在xoy平面作长度为A12e、与x夹角为
Figure BDA0003371841520000036
的线段EF，过E点作与x夹角为
Figure BDA0003371841520000037
直线a，过F点作与x夹角为
Figure BDA0003371841520000038
直线b，a和b交相交于点，量取线段OE长度即为A1e，量取线段OF长度即为A2e。具体过程如图1所示，图中
Figure BDA0003371841520000039
②Ω为地球周期对应的角频率，计算
Figure BDA00033718415200000310
Figure BDA00033718415200000311
Figure BDA00033718415200000312
将lon12R(t)中信息与lonE1R(t)、-lonE2R(t)中信息进行特征比对，lon12R(t)中与lonE1R(t)中信息特征符合的特征信息为D1(t)，与-lonE2R(t)中信息特征符合的特征信息为D2(t)，则
Figure BDA00033718415200000313
(2)惯导纬度误差中包含垂线偏差引起的部分和惯导自身误差源引起的部分，求取垂线偏差必须分离这两部分。两套惯导纬度差值lat12(t)中不包含垂线偏差引起的部分，主要是两套惯导自身误差源引起部分的合成，因而基于lat12(t)可估计2套惯导自身误差源引起的纬度误差项latE1a(t)、latE2a(t)，具体方法为：
①提取lat12(t)中地球周期振荡项幅值A12n、相位
Figure BDA00033718415200000314
依据
Figure BDA00033718415200000315
和A12n、
Figure BDA00033718415200000316
构建三角形求得两套惯导纬度误差中地球周期振荡项幅值A1n、A2n。在xoy平面作长度为A12n、与x夹角为
Figure BDA00033718415200000317
的线段EF，过E点作与x夹角为
Figure BDA00033718415200000318
直线a，过F点作与x夹角为
Figure BDA00033718415200000319
直线b，a和b交相交于点，量取线段OE长度即为A1n，量取线段OF长度即为A2n。具体过程如图1所示，图中
Figure BDA00033718415200000320
②Ω为地球周期对应的角频率，计算
Figure BDA00033718415200000321
Figure BDA00033718415200000322
Figure BDA00033718415200000323
将lat12R(t)中信息与latE1R(t)、-latE2R(t)中信息进行特征比对，lat12R(t)中与latE1R(t)中信息特征符合的特征信息为F1(t)，与-latE2R(t)中信息特征符合的特征信息为F2(t)，则
Figure BDA0003371841520000041
(3)则L测线南北向、东西向垂线偏差分别为：
gn(t)＝latE1(t)-latE1a(t)+gn0或gn(t)＝latE2(t)-latE2a(t)+gn0
Figure BDA0003371841520000042
或
Figure BDA0003371841520000043
其中
Figure BDA0003371841520000044
为纬度。
需要强调的是，本发明所述的实施例是说明性的，而不是限定性的，因此本发明包括并不限于具体实施方式中所述的实施例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，同样属于本发明保护的范围。

Claims (3)

1.一种基于两套惯导互观测的海洋垂线偏差测量方法，其特征在于：包括以下步骤：
步骤1、载体航行至垂线偏差基准点，将两套惯导位置输出校准为基准点准确位置，记录基准点东西方向、南北方向垂线偏差值ge0、gn0；
步骤2、沿东西方向设计垂线偏差测线，载体进入测线前在测线同纬度且垂线偏差变化平缓区域游弋一段时间，以卫导经纬度为基准观测两套惯导经纬度误差，提取两套惯导纬度误差中地球周期振荡相位
Figure FDA0003371841510000011
经度误差中地球周期振荡相位
Figure FDA0003371841510000012
步骤3、载体进入垂线偏差测线L航行，期间以卫导经纬度为基准，得到两套惯导的经度误差lonE1(t)、lonE2(t)和纬度误差latE1(t)、latE2(t)，基于两套惯导经度差值lon12(t)、纬度差值lat12(t)估计两套惯导自身误差源引起的经度误差项lonE1a(t)、lonE2a(t)和纬度误差项latE1a(t)、latE2a(t)，则L测线南北向垂线偏差为：
gn(t)＝latE1(t)-latE1a(t)+gn0或gn(t)＝latE2(t)-latE2a(t)+gn0
东西向垂线偏差为：
Figure FDA0003371841510000013
或
Figure FDA0003371841510000014
其中
Figure FDA0003371841510000015
为纬度。
2.根据权利要求1所述的一种基于两套惯导互观测的海洋垂线偏差测量方法，其特征在于：所述步骤3中基于两套惯导经度差值lon12(t)估计2套惯导自身误差源引起的经度误差项lonE1a(t)、lonE2a(t)的具体方法为：
(1)提取lon12(t)中地球周期振荡项幅值A12e、相位
Figure FDA0003371841510000016
依据
Figure FDA0003371841510000017
和A12e、
Figure FDA0003371841510000018
构建三角形求得两套惯导经度误差中地球周期振荡项幅值A1e、A2e；
(2)Ω为地球周期对应的角频率，计算
Figure FDA0003371841510000019
Figure FDA00033718415100000110
Figure FDA00033718415100000111
将lon12R(t)中信息与lonE1R(t)、-lonE2R(t)中信息进行特征比对，lon12R(t)中与lonE1R(t)中信息特征符合的特征信息为D1(t)，与-lonE2R(t)中信息特征符合的特征信息为D2(t)，则
Figure FDA00033718415100000112
3.根据权利要求1所述的一种基于两套惯导互观测的海洋垂线偏差测量方法，其特征在于：所述步骤3中基于2套惯导纬度差值lat12(t)估计2套惯导自身误差源引起的纬度误差项latE1a(t)、latE2a(t)的具体方法为：
(1)提取lat12(t)中地球周期振荡项幅值A12n、相位
Figure FDA00033718415100000113
依据
Figure FDA00033718415100000114
和A12n、
Figure FDA00033718415100000115
构建三角形求得两套惯导纬度误差中地球周期振荡项幅值A1n、A2n；
(2)Ω为地球周期对应的角频率，计算
Figure FDA00033718415100000116
Figure FDA00033718415100000117
Figure FDA00033718415100000118
将lat12R(t)中信息与latE1R(t)、-latE2R(t)中信息进行特征比对，lat12R(t)中与latE1R(t)中信息特征符合的特征信息为F1(t)，与-latE2R(t)中信息特征符合的特征信息为F2(t)，则
Figure FDA0003371841510000021
CN202111404563.3A 2021-11-24 2021-11-24 一种基于两套惯导互观测的海洋垂线偏差测量方法 Active CN114137625B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111404563.3A CN114137625B (zh) 2021-11-24 2021-11-24 一种基于两套惯导互观测的海洋垂线偏差测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111404563.3A CN114137625B (zh) 2021-11-24 2021-11-24 一种基于两套惯导互观测的海洋垂线偏差测量方法

Publications (2)

Publication Number Publication Date
CN114137625A true CN114137625A (zh) 2022-03-04
CN114137625B CN114137625B (zh) 2023-04-28

Family

ID=80391781

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111404563.3A Active CN114137625B (zh) 2021-11-24 2021-11-24 一种基于两套惯导互观测的海洋垂线偏差测量方法

Country Status (1)

Country Link
CN (1) CN114137625B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117346772A (zh) * 2023-12-04 2024-01-05 华中光电技术研究所(中国船舶集团有限公司第七一七研究所) 一种适用于水下平台惯性导航系统的误差校正方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1118663A (en) * 1965-06-03 1968-07-03 North American Aviation Inc Inertial navigation system error correcting methods
CN107677292A (zh) * 2017-09-28 2018-02-09 中国人民解放军国防科技大学 基于重力场模型的垂线偏差补偿方法
US20210293543A1 (en) * 2005-06-22 2021-09-23 James R. Huddle Apparatus and Method for Integrating Continuous and Discontinuous Inertial Instrument

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1118663A (en) * 1965-06-03 1968-07-03 North American Aviation Inc Inertial navigation system error correcting methods
US20210293543A1 (en) * 2005-06-22 2021-09-23 James R. Huddle Apparatus and Method for Integrating Continuous and Discontinuous Inertial Instrument
CN107677292A (zh) * 2017-09-28 2018-02-09 中国人民解放军国防科技大学 基于重力场模型的垂线偏差补偿方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DONGKAI DAI等: "Dynamic measurement of high-frequency deflections of the vertical based on the observation of INS/GNSS integration attitude error" *
安文;许江宁;吴苗;江鹏飞;: "垂线偏差对捷联惯导姿态对准的影响和补偿" *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117346772A (zh) * 2023-12-04 2024-01-05 华中光电技术研究所(中国船舶集团有限公司第七一七研究所) 一种适用于水下平台惯性导航系统的误差校正方法及系统
CN117346772B (zh) * 2023-12-04 2024-02-20 华中光电技术研究所(中国船舶集团有限公司第七一七研究所) 一种适用于水下平台惯性导航系统的误差校正方法及系统

Also Published As

Publication number Publication date
CN114137625B (zh) 2023-04-28

Similar Documents

Publication Publication Date Title
CN110487301B (zh) 一种雷达辅助机载捷联惯性导航系统初始对准方法
CN103630137B (zh) 一种用于导航系统的姿态及航向角的校正方法
CN104597471B (zh) 面向时钟同步多天线gnss接收机的定向测姿方法
CN103217159B (zh) 一种sins/gps/偏振光组合导航系统建模及动基座初始对准方法
CN103323007B (zh) 一种基于时变量测噪声的鲁棒联邦滤波方法
CN103674030B (zh) 基于天文姿态基准保持的垂线偏差动态测量装置和方法
CN102486377B (zh) 一种光纤陀螺捷联惯导系统初始航向的姿态获取方法
CN103591949B (zh) 三轴姿态测量系统非正交性误差的正交补偿方法
CN105091907B (zh) Sins/dvl组合中dvl方位安装误差估计方法
CN101949703A (zh) 一种捷联惯性/卫星组合导航滤波方法
CN103900565A (zh) 一种基于差分gps的惯导系统姿态获取方法
CN109282804B (zh) 一种单轴光纤陀螺寻北算法
CN109085655B (zh) 一种水下平台重力测量方案与验证方法
CN103344260A (zh) 基于rbckf的捷联惯导系统大方位失准角初始对准方法
CN103674059A (zh) 一种基于外测速度信息的sins水平姿态误差修正方法
CN104697520A (zh) 一体化无陀螺捷联惯导系统与gps系统组合导航方法
CN111722295A (zh) 一种水下捷联式重力测量数据处理方法
CN103791918A (zh) 一种舰船捷联惯导系统极区动基座对准方法
CN114137625A (zh) 一种基于两套惯导互观测的海洋垂线偏差测量方法
CN110514200B (zh) 一种惯性导航系统及高转速旋转体姿态测量方法
CN113900069B (zh) 一种基于干涉成像高度计的垂线偏差计算方法及其系统
CN109085656B (zh) 一种面向特征的高精度重力图构建与插值方法
CN103901459B (zh) 一种mems/gps组合导航系统中量测滞后的滤波方法
CN108007457B (zh) 一种基于细分时间片的监控导航系统异步数据融合方法
CN112292578B (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