CN114114443A - 一种基于水平重力势-位置转化的海洋垂线偏差测量方法 - Google Patents

一种基于水平重力势-位置转化的海洋垂线偏差测量方法 Download PDF

Info

Publication number
CN114114443A
CN114114443A CN202111406149.6A CN202111406149A CN114114443A CN 114114443 A CN114114443 A CN 114114443A CN 202111406149 A CN202111406149 A CN 202111406149A CN 114114443 A CN114114443 A CN 114114443A
Authority
CN
China
Prior art keywords
deviation
vertical
inertial
east
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
CN202111406149.6A
Other languages
English (en)
Other versions
CN114114443B (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 CN202111406149.6A priority Critical patent/CN114114443B/zh
Publication of CN114114443A publication Critical patent/CN114114443A/zh
Application granted granted Critical
Publication of CN114114443B publication Critical patent/CN114114443B/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
    • 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/16—Measuring gravitational fields or waves; Gravimetric prospecting or detecting specially adapted for use on moving platforms, e.g. ship, aircraft
    • 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)
  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Ocean & Marine Engineering (AREA)
  • Navigation (AREA)

Abstract

本发明涉及一种基于水平重力势‑位置转化的海洋垂线偏差测量方法，其技术特点是：标定不同周期变化的东西/南北方向垂线偏差与其引起的经度/纬度输出量值的转换系数函数；将东向/北向加速度计原始输出分别减去载体所在位置东西/南北方向垂线偏差图值后作为东向/北向加速度计的输出数据，消除惯性自主导航中的经度/纬度振荡发散误差；将λ海里波长的垂线偏差信息变换为周期的时间信号所需的载体速度；求得测线L上东西/南北方向垂线偏差测量值。本发明设计合理，实现了垂线偏差的精确测量，可用于支撑海洋重力场信息建设与应用。

Description

一种基于水平重力势-位置转化的海洋垂线偏差测量方法
技术领域
本发明航洋重力场测量技术领域，尤其是一种基于水平重力势-位置转化的海洋垂线偏差测量方法。
背景技术
地球重力场是反映地球物质分布特征的物理场，是地球引力场和惯性离心力场的合成场，分布于地球表面及其临近空间。重力场作为天然的物理资源，早已广泛用于各类系统中，其中比较典型的应用是惯性导航系统，但对于惯导而言，重力场这个物理场“并不完美”。“完美”的重力场，也即正常重力场模型将地球视为旋转椭球体，重力矢量指向椭球体法向。指向参考椭球法向的重力矢量在惯性导航系统中被用作与地球相关的位置参考基准。由于地球形状的不规则和内部质量分布的不均匀性，实际重力场变化非常复杂，实际观测重力场与理论上正常重力场存在着偏差，称为扰动重力场，其中二者的标量差称为重力异常，方向差称为垂线偏差。垂线偏差的存在，实际重力矢量偏离了参考椭球体法向，作为位置参考基准精度下降，必须依赖精确的重力信息予以补偿。
重力信息包含重力异常和垂线偏差信息，重力信息测量一般使用力传感器，由于力传感器难以区分运动加速度和重力加速度，必须考虑载体运动干扰的补偿。目前趋于成熟的海洋重力测量技术是指垂向重力异常的测量，其利用卫导信息对载体垂向加速度进行补偿；相比于重力异常测量，海洋运动平台水平加速度的补偿更为困难，尚未有有效的海洋垂线偏差船载测量方法。海洋垂线偏差信息主要依据卫星测高数据解算获取，空间分辨率不足且精度有限。
垂线偏差在全海域普遍存在，大部分表现为3-5角秒幅值的缓慢变化，基于卫星测高数据的垂线偏差图可以较为精确的建模这些缓慢变化区域的垂线偏差；但受限于空间分辨率，对于大洋中也较为广泛存在的小空间尺度上剧烈的垂线偏差变化，基于卫星测高数据的垂线偏差图建模能力明显下降，不能满足应用需求，如何有效提高船载高精度海洋垂线偏差测量精度是目前迫切需要解决的问题。
发明内容
本发明的目的在于克服现有技术的不足，提供一种基于水平重力势-位置转化的海洋垂线偏差测量方法，提高了船载海洋垂线偏差测量的精度。
本发明解决现有的技术问题是采取以下技术方案实现的：
一种基于水平重力势-位置转化的海洋垂线偏差测量方法，包括以下步骤：
步骤1、标定不同周期变化的东西/南北方向垂线偏差与其引起的经度/纬度输出量值的转换系数函数；
步骤2、启动惯导，将东向/北向加速度计原始输出分别减去载体所在位置东西/南北方向垂线偏差图值后作为东向/北向加速度计的输出数据，惯导完成对准后转入惯/卫组合导航解算程序M1，用于消除惯性自主导航中的经度/纬度振荡发散误差；
步骤3、载体航行进入预设的垂线偏差测量测线，通过设定航速将测线实际垂线偏差与垂线偏差图值的差值所包含波长信息变换为2倍舒勒振荡周期，将λ海里波长的垂线偏差信息变换为周期Tm的时间信号所需的载体速度vλ＝λ/Tm；
步骤4、在测线L航行时间段tL同步记录惯性自主导航经度Lon(tL)和纬度Lat(tL)、卫导经度BDLon(tL)和纬度BDLat(tL)，求取两者经度差值Lonc(tL)＝Lon(tL)-BDLon(tL)和纬度Latc(tL)＝Lat(tL)-BDLat(tL)，利用BDLon(tL)、BDLat(tL)将时间序列Lonc(tL)和Latc(tL)转换为空间序列Lonc(L)和Latc(L)，设测线L上东西/南北方向垂线偏差图值分别为Get(L)和Gnt(L)，根据步骤1得到的转换系数函数函数fe、fn获取2T周期附近的平均转换系数fe1、fn1，则测线L上东西/南北方向垂线偏差测量值分别为Ge(L)＝fe1*Lonc(L)+Get(L)，Gn(L)＝fn1*Latc(L)+Gnt(L)。
而且，所述步骤2去除惯性自主导航中的经度/纬度振荡发散误差方法包括以下步骤：
⑴设Td为半个地球周期，t0+nTd时刻以惯/卫组合导航解算程序M1获得的载体位置和姿态信息作为初始值启动惯性自主导航解算程序Jn，Jn输出经度为Lont0+nTd(t)、纬度为Latt0+nTd(t)，；
⑵按下式计算用于垂线偏差测量解算的惯导经度/纬度：(t0+nTd,t0+(n+1)Td)时间区间内t时刻惯性自主导航经度输出Lon(t)＝0.5(Lont0+(n-2)Td(t)+Lont0+(n-1)Td(t))，惯性自主导航纬度输出Lat(t)＝0.5(Latt0+(n-2)Td(t)+Latt0+(n-1)Td(t))
上式中，n＝0,1,2,3,...。
本发明的优点和积极效果是：
本发明设计合理，其利用加速度传感器无差别感知运动加速度和重力加速度的特性，把水平方向的重力势转化为惯导自主导航位置输出中的信息；采取措施有效抑制了干扰垂线偏差测量的惯导位置输出中的地球周期振荡部分；通过补偿垂线偏差图信息将待测垂线偏差变换为窄带信号，并利用惯导频响特性放大了这一窄带垂线偏差信号，从而实现了垂线偏差的精确测量，可用于支撑海洋重力场信息建设与应用。
附图说明
图1为本实施例给出的不同变化周期的东西方向垂线偏差—经度输出的转换关系图；
图2为本实施例给出的某轨迹上实际垂线偏差与垂涎偏差图信息比对图。
具体实施方式
以下结合附图对本发明实施例做进一步详述。
本发明是在高精度惯导、精密差分卫导和基于卫星测高数据解算的垂线偏差图构成海洋垂线偏差测量系统上实现的。其利用加速度计物理传感和惯导解算原理，将东西方向垂线偏差势能转化为惯导经度输出中对应信息，将南北方向垂线偏差势能转化为惯导纬度输出中对应信息，通过错时启动的惯导自主解算组合有效抑制惯导输出中地球振荡周期变化的干扰信息，通过比对惯导与卫导经纬度信息分离出惯导经纬度输出中由垂线偏差引起的部分，并按照事先标定的东西方向垂线偏差到经度输出、南北方向垂线偏差到纬度输出的映射关系，解算获取东西方向和南北方向的垂线偏差信息。
基于上述说明，本发明提出一种基于水平重力势-位置转化的海洋垂线偏差测量方法，包括以下步骤：
步骤1、标定不同周期变化的东西/南北方向垂线偏差与其引起的经度/纬度输出量值的转换系数函数。
高精度惯导工作于阻尼状态，以数据注入的方式在惯导东向/北向加速度计输出数据上叠加代表东西/南北方向垂线偏差的模拟激励数据，标定不同周期变化的东西/南北方向垂线偏差与其引起的经度/纬度输出量值的转换系数函数fe、fn，单位为米/角秒。不同变化周期的东西方向垂线偏差到经度输出的转换系数标定结果示例如图1所示，该图中T为舒勒周期，约等于84.4分钟。
步骤2、启动惯导，将东向/北向加速度计原始输出分别减去载体所在位置东西/南北方向垂线偏差图值后作为东向/北向加速度计的输出数据，惯导完成对准后转入惯/卫组合导航解算程序M1。进行垂线偏差测量时惯性自主导航由于其传感器漂移等误差源而产生的经度/纬度地球振荡发散误差将对垂线偏差测量解算造成干扰，采取下述操作去除惯性自主导航中的经度/纬度振荡发散误差：
①令Td为半个地球周期，t0+nTd(n＝0,1,2,3,...)时刻以惯/卫组合导航解算程序M1获得的载体位置和姿态信息作为初始值启动惯性自主导航解算程序Jn，Jn输出经度为Lont0+nTd(t)、纬度为Latt0+nTd(t)；
②用于垂线偏差测量解算的惯导经度/纬度输出如下：(t0+nTd,t0+(n+1)Td)(n＝2,3,...)时间区间内t时刻惯性自主导航经度输出Lon(t)＝0.5(Lont0+(n-2)Td(t)+Lont0+(n-1)Td(t))，惯性自主导航纬度输出Lat(t)＝0.5(Latt0+(n-2)Td(t)+Latt0+(n-1)Td(t))。
步骤3、载体航行进入预设的垂线偏差测量测线，保持匀速直航状态。由于步骤2中利用垂线偏差图信息对惯导加速度计输出进行了补偿，测量的对象为测线实际垂线偏差与测线垂线偏差图值的差值，实际垂线偏差与垂线偏差图值比对示例如图2所示。测线实际垂线偏差可能包含多种波长信息，但其与垂线偏差图值差值则主要表现有限波长区间信息，一般波长范围在8-12海里区间。综合考虑惯性导航系统频响特性和测量效率，通过合理设计航速将测线实际垂线偏差与垂线偏差图值的差值包含波长信息变换为2倍舒勒振荡周期即2T附近的时间信号。λ海里波长的垂线偏差信息变换为周期Tm的时间信号所需的载体速度vλ按下式计算，vλ＝λ/Tm。如采用3.5kn速度时，可将8-12海里区间波长的垂线偏差信息变换为1.6T-2.4T周期区间时间信号。
步骤4、在测线L航行时间段tL同步记录惯性自主导航经度Lon(tL)，纬度Lat(tL)，卫导经度BDLon(tL)，纬度BDLat(tL)，求取差值Lonc(tL)＝Lon(tL)-BDLon(tL)、纬度Latc(tL)＝Lat(tL)-BDLat(tL)。利用BDLon(tL)、BDLat(tL)将时间序列Lonc(tL)和Latc(tL)转换为空间序列Lonc(L)和Latc(L)，测线L上东西/南北方向垂线偏差图值分别为Get(L)，Gnt(L)。利用步骤1中的函数fe、fn获取2T周期附近的平均转换系数fe1、fn1，则测线L上东西/南北方向垂线偏差测量值分别为Ge(L)＝fe1*Lonc(L)+Get(L)，Gn(L)＝fn1*Latc(L)+Gnt(L)。
需要强调的是，本发明所述的实施例是说明性的，而不是限定性的，因此本发明包括并不限于具体实施方式中所述的实施例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，同样属于本发明保护的范围。

Claims (2)

1.一种基于水平重力势-位置转化的海洋垂线偏差测量方法，其特征在于：包括以下步骤：
步骤1、标定不同周期变化的东西/南北方向垂线偏差与其引起的经度/纬度输出量值的转换系数函数；
步骤2、启动惯导，将东向/北向加速度计原始输出分别减去载体所在位置东西/南北方向垂线偏差图值后作为东向/北向加速度计的输出数据，惯导完成对准后转入惯/卫组合导航解算程序M1，用于消除惯性自主导航中的经度/纬度振荡发散误差；
步骤3、载体航行进入预设的垂线偏差测量测线，通过设定航速将测线实际垂线偏差与垂线偏差图值的差值所包含波长信息变换为2倍舒勒振荡周期，将λ海里波长的垂线偏差信息变换为周期Tm的时间信号所需的载体速度vλ＝λ/Tm；
步骤4、在测线L航行时间段tL同步记录惯性自主导航经度Lon(tL)和纬度Lat(tL)、卫导经度BDLon(tL)和纬度BDLat(tL)，求取两者经度差值Lonc(tL)＝Lon(tL)-BDLon(tL)和纬度Latc(tL)＝Lat(tL)-BDLat(tL)，利用BDLon(tL)、BDLat(tL)将时间序列Lonc(tL)和Latc(tL)转换为空间序列Lonc(L)和Latc(L)，设测线L上东西/南北方向垂线偏差图值分别为Get(L)和Gnt(L)，根据步骤1得到的转换系数函数函数fe、fn获取2T周期附近的平均转换系数fe1、fn1，则测线L上东西/南北方向垂线偏差测量值分别为Ge(L)＝fe1*Lonc(L)+Get(L)，Gn(L)＝fn1*Latc(L)+Gnt(L)。
2.根据权利要求1所述的一种基于水平重力势-位置转化的海洋垂线偏差测量方法，其特征在于：所述步骤2去除惯性自主导航中的经度/纬度振荡发散误差方法包括以下步骤：
⑴设Td为半个地球周期，t0+nTd时刻以惯/卫组合导航解算程序M1获得的载体位置和姿态信息作为初始值启动惯性自主导航解算程序Jn，Jn输出经度为Lont0+nTd(t)、纬度为Latt0+nTd(t)，；
⑵按下式计算用于垂线偏差测量解算的惯导经度/纬度：(t0+nTd,t0+(n+1)Td)时间区间内t时刻惯性自主导航经度输出Lon(t)＝0.5(Lont0+(n-2)Td(t)+Lont0+(n-1)Td(t))，惯性自主导航纬度输出Lat(t)＝0.5(Latt0+(n-2)Td(t)+Latt0+(n-1)Td(t))
上式中，n＝0,1,2,3,...。
CN202111406149.6A 2021-11-24 2021-11-24 一种基于水平重力势-位置转化的海洋垂线偏差测量方法 Active CN114114443B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111406149.6A CN114114443B (zh) 2021-11-24 2021-11-24 一种基于水平重力势-位置转化的海洋垂线偏差测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111406149.6A CN114114443B (zh) 2021-11-24 2021-11-24 一种基于水平重力势-位置转化的海洋垂线偏差测量方法

Publications (2)

Publication Number Publication Date
CN114114443A true CN114114443A (zh) 2022-03-01
CN114114443B CN114114443B (zh) 2023-05-19

Family

ID=80371981

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111406149.6A Active CN114114443B (zh) 2021-11-24 2021-11-24 一种基于水平重力势-位置转化的海洋垂线偏差测量方法

Country Status (1)

Country Link
CN (1) CN114114443B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5673191A (en) * 1995-04-10 1997-09-30 Atlantic Richfield Company Method and apparatus for identifying geological structures using wavelet analysis of potential fields
CN102175213A (zh) * 2011-01-30 2011-09-07 中国测绘科学研究院 Gps水准高程传递方法
CN106526700A (zh) * 2016-10-10 2017-03-22 北京航天控制仪器研究所 一种利用地球重力场评估重力传感器分辨率的方法
CN109085655A (zh) * 2018-09-19 2018-12-25 中国船舶重工集团公司第七0七研究所 一种水下平台重力测量方案与验证方法
CN110703355A (zh) * 2019-09-19 2020-01-17 东华理工大学 一种星载加速度计的校准方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5673191A (en) * 1995-04-10 1997-09-30 Atlantic Richfield Company Method and apparatus for identifying geological structures using wavelet analysis of potential fields
CN102175213A (zh) * 2011-01-30 2011-09-07 中国测绘科学研究院 Gps水准高程传递方法
CN106526700A (zh) * 2016-10-10 2017-03-22 北京航天控制仪器研究所 一种利用地球重力场评估重力传感器分辨率的方法
CN109085655A (zh) * 2018-09-19 2018-12-25 中国船舶重工集团公司第七0七研究所 一种水下平台重力测量方案与验证方法
CN110703355A (zh) * 2019-09-19 2020-01-17 东华理工大学 一种星载加速度计的校准方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
毛经伦，等: "利用EGM2008地球重力场模型 研究川西地区水平构造应力场" *

Also Published As

Publication number Publication date
CN114114443B (zh) 2023-05-19

Similar Documents

Publication Publication Date Title
Wang et al. A novel SINS/DVL tightly integrated navigation method for complex environment
CN111323050B (zh) 一种捷联惯导和多普勒组合系统标定方法
CN101949703B (zh) 一种捷联惯性/卫星组合导航滤波方法
CN102279002B (zh) 星敏感器测量坐标系与载体坐标系转换矩阵的标定方法
CN102252677A (zh) 一种基于时间序列分析的变比例自适应联邦滤波方法
CN106500693A (zh) 一种基于自适应扩展卡尔曼滤波的ahrs算法
CN109085655B (zh) 一种水下平台重力测量方案与验证方法
Vasconcelos et al. INS/GPS aided by frequency contents of vector observations with application to autonomous surface crafts
CN110763872A (zh) 一种多普勒测速仪多参数在线标定方法
CA2256964C (en) Method of locating hydrophones
CN109507706A (zh) 一种gps信号丢失的预测定位方法
CN111812737B (zh) 水下导航与重力测量一体化系统
CN110631573B (zh) 一种惯性/里程计/全站仪多信息融合方法
CN111982105B (zh) 一种基于sins/lbl紧组合的水下导航定位方法及系统
Kiselev et al. Autonomous underwater robot as an ideal platform for marine gravity surveys
CN111504278A (zh) 基于自适应频域积分的海浪检测方法
CN113900069B (zh) 一种基于干涉成像高度计的垂线偏差计算方法及其系统
Song et al. Novel underwater glider-based absolute oceanic current observation solutions
CN109084755B (zh) 一种基于重力视速度与参数辨识的加速度计零偏估计方法
CN111220146A (zh) 一种基于高斯过程回归学习的水下地形匹配定位方法
CN114323065A (zh) 基于多手段融合的水下自主导航系统误差监测与估计方法
CN111982126B (zh) 一种全源BeiDou/SINS弹性状态观测器模型设计方法
CN116559966B (zh) 基于sins/ldv组合的重力测量方法及系统
CN114114443B (zh) 一种基于水平重力势-位置转化的海洋垂线偏差测量方法
Thompson Altimeter data and geoid error in mesoscale ocean prediction: Some results from a primitive equation model

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