CN110068876B - 基于载体自振动航空重力梯度仪运动误差补偿方法 - Google Patents

基于载体自振动航空重力梯度仪运动误差补偿方法 Download PDF

Info

Publication number
CN110068876B
CN110068876B CN201910464038.7A CN201910464038A CN110068876B CN 110068876 B CN110068876 B CN 110068876B CN 201910464038 A CN201910464038 A CN 201910464038A CN 110068876 B CN110068876 B CN 110068876B
Authority
CN
China
Prior art keywords
accelerometer
vertical
signals
frequency
signal
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
CN201910464038.7A
Other languages
English (en)
Other versions
CN110068876A (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 CN201910464038.7A priority Critical patent/CN110068876B/zh
Publication of CN110068876A publication Critical patent/CN110068876A/zh
Application granted granted Critical
Publication of CN110068876B publication Critical patent/CN110068876B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V13/00Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Navigation (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种基于载体自振动航空重力梯度仪运动误差补偿方法,其技术特点在于:包括以下步骤:步骤1、对加速度计组合信号在飞机螺旋桨旋转频率上进行正余弦解调;步骤2、对垂向加速度计输出信号在飞机螺旋桨旋转频率上进行正余弦解调;步骤3、求得两组加速度计敏感轴垂向安装误差实时信号;步骤4、对垂向加速度计输出信号在二倍频频点上进行正余弦解调;步骤5、得到重力梯度信号中加速度计敏感轴俯仰角安装误差予以补偿。本发明合理补偿航空重力梯度仪的垂向加速度误差,提高了系统的测量精度。

Description

基于载体自振动航空重力梯度仪运动误差补偿方法
技术领域
本发明属于旋转加速度计式航空重力梯度仪技术领域,涉及旋转加速度计式航空重力梯度仪垂向运动误差补偿方法,尤其是一种基于载体自振动航空重力梯度仪运动误差补偿方法。
背景技术
目前,旋转加速度计式重力梯度仪是对地球表面微小重力梯度变化进行连续测量的仪器。如图2所示,作为核心敏感器的重力梯度测量组件基于加速度计位置差分测量原理,通过机械旋转的方式将重力梯度张量分量调制到系统旋转频率的二倍频处,加速度计输出与重力梯度张量分量之间的关系可表示为:
(a1+a3)-(a2+a4)=2R(Γxxyy)sin2ωt+4RΓxycos2ωt (1)
式中a1、a2、a3、a4是四个加速度计的输出,R是加速度计检测质心到旋转中心的距离,Γxx、Γyy、Γxy是重力梯度张量分量,ω是旋转装置的旋转角速度。将系统输出的加速度计组合信号进行2ω频率上的同步解调得到最终的重力梯度张量信号Γxxyy和Γxy。同时定义Γuv=Γxxyy。进行航空测量时,由惯性稳定平台承载重力梯度敏感器,隔离载体角运动对重力梯度测量的影响,同时为敏感器提供基准坐标系。
利用飞机等机动平台开展重力梯度勘查具有效率高、成本低、环境适应性强的优点,可在人员难以到达的复杂地理环境中部署使用。但在机载测量条件下,飞机在飞行过程中垂向振动的二倍频分量与四只加速度计的敏感轴垂向安装误差角相耦合,进入加速度计组合信号的二倍频分量上,对重力梯度测量信号造成干扰。经测量,机载垂向运动加速度的二倍频分量约为2×10-3g量级,在每只加速度计敏感轴垂向安装误差精度优于5″的条件下,对重力梯度信号的干扰为5000E。因此,机载条件下对载体垂向加速度误差的补偿是影响仪器测量性能的重要环节,需通过特定的方法予以补偿。
发明内容
本发明的目的在于克服现有技术的不足,提出一种基于载体自振动航空重力梯度仪运动误差补偿方法,能够合理补偿旋转加速度计式重力梯度仪在机载测量时,输出的重力梯度信号中存在由载体垂向加速度与加速度计敏感轴垂向安装误差相耦合引起的测量误差。
本发明解决其现实问题是采取以下技术方案实现的:
一种旋转加速度计式航空重力梯度仪垂向线运动测量误差的补偿方法,包括以下步骤:
步骤1、对加速度计组合信号在飞机螺旋桨旋转频率上进行正余弦解调,得到两路信号
Figure BDA0002078908370000021
Figure BDA0002078908370000022
其中
Figure BDA0002078908370000023
为加速度计组合信号在飞机螺旋桨旋转频率上的正弦分量,
Figure BDA0002078908370000024
为加速度计组合信号在飞机螺旋桨旋转频率上的余弦分量;
步骤2、对垂向加速度计输出信号在飞机螺旋桨旋转频率上上进行正余弦解调,得到两路信号
Figure BDA0002078908370000025
Figure BDA0002078908370000026
其中
Figure BDA0002078908370000027
为垂向加速度在飞机螺旋桨旋转频率上的正弦分量,
Figure BDA0002078908370000028
为垂向加速度在飞机螺旋桨旋转频率上的余弦分量;
步骤3、将由从步骤1的加速度计组合信号得到的
Figure BDA0002078908370000031
Figure BDA0002078908370000032
信号与从步骤2的垂向加速度计输出信号得到的
Figure BDA0002078908370000033
Figure BDA0002078908370000034
信号进行相除处理,得到两组加速度计敏感轴垂向安装误差实时信号Δβs和Δβc,将两组信号Δβs和Δβc求平均得到更加真实的Δβ实时信号。其中Δβs是由飞机螺旋桨旋转频率正弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβc是由飞机螺旋桨旋转频率余弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβ是两个安装误差实时信号的平均值;
步骤4、对垂向加速度计输出信号在二倍频频点上进行正余弦解调,得到两路信号
Figure BDA0002078908370000035
Figure BDA0002078908370000036
Figure BDA0002078908370000037
是垂向加速度在旋转二倍频上的正弦分量,
Figure BDA0002078908370000038
是垂向加速度在旋转二倍频上的余弦分量;
步骤5、将由从垂向加速度计输出信号得到的
Figure BDA0002078908370000039
Figure BDA00020789083700000310
信号与Δβ信号进行相乘处理,得到重力梯度信号中加速度计敏感轴俯仰角安装误差ΔΓuv(az)和ΔΓxy(az),予以补偿。
本发明的优点和有益效果:
1、本发明能够依据固定翼飞机飞行测量时的运动特点和重力梯度测量的误差机理,利用本是误差源的飞机螺旋桨振动作为加速度计敏感轴垂向安装误差的观测激励项,对其进行实时监测,并依据该项安装误差实时监测值对重力梯度测量信号中的垂向加速度误差进行补偿。
2、本发明针对动态飞行测量的特点和误差原理,提出利用本是误差源的飞机自振动作为观测激励项,合理补偿航空重力梯度仪的垂向加速度误差,提高了系统的测量精度;
3、本发明在不依靠外界信息、不增加系统的复杂程度的基础上,仅使用系统稳定平台的惯性导航数据对重力梯度信号中的垂向加速度误差进行补偿,提高了航空重力梯度仪系统的可靠性。
附图说明
图1为本发明的机载试验中飞机平飞过程中的垂向振动谱示意图;
图2为本发明的旋转加速度计式重力梯度敏感器测量原理示意图。
具体实施方式
以下结合附图对本发明实施例作进一步详述:
一种基于载体自振动航空重力梯度仪运动误差补偿方法,可以用于补偿机载测量时由固定翼飞机垂向运动引起的重力梯度测量误差。
在本实施例中,机载条件下重力梯度仪由加速度计敏感轴垂向安装误差带来的重力梯度敏感器输出信号误差为:
ΔAout(Δβi)=a0z[(Δβ1+Δβ3)-(Δβ2+Δβ4)]
+[Γxz(Δβ2-Δβ4)-Γyz(Δβ1-Δβ3)]sinωt
+[Γxz(Δβ1-Δβ3)+Γyz(Δβ2-Δβ4)]cosωt (2)
式中,ΔAout(Δβi)是由加速度计敏感轴垂向安装误差引起的重力梯度敏感器输出加速度计组合信号的误差,βi(i=1,2,3,4)是第i个加速度计敏感轴的垂向角度安装误差,Γxz、Γyz是对应方向的重力梯度张量,a0z是载体的垂向加速度。
从中可以看出,误差第二项、第三项其载波频率为一倍频,可由同步解调技术除去,误差第一项表示载体垂向加速度的二倍频分量与敏感轴垂向安装误差角相耦合,进入重力梯度测量通道。因此,重力梯度测量信号中由载体垂向加速度引起的误差可表示为:
Figure BDA0002078908370000041
式中:
Figure BDA0002078908370000042
为垂向加速度的正弦二倍频分量,
Figure BDA0002078908370000043
为垂向加速度的余弦二倍频分量。经测量,机载垂向运动加速度的二倍频分量约为2×10-3g量级,在安装精度|Δβi|<5″的条件下对解调后重力梯度信号造成的影响量级为5000E,在系统动态使用时需对该项误差进行补偿。
考虑在测量过程中重力梯度仪处于高动态的飞行条件下,重力梯度敏感器中加速度计敏感轴垂向安装误差(Δβ=(Δβ1+Δβ3)-(Δβ2+Δβ4))在机载高机动过程中变化0.1″,对重力梯度信号的影响为20E量级,不可忽略,因此需对加速度计敏感轴垂向安装误差进行实时监测并采取后处理方式进行补偿。
图1为机载试验中固定翼飞机平飞运动中的垂向振动谱,可以观察到在113.2Hz处有一个0.01g的分量,为固定翼飞机螺旋桨的工作频率。可以利用此频点上的振动激励对加速度计敏感轴垂向安装误差进行实时观测并在后处理的过程中予以补偿。
一种旋转加速度计式航空重力梯度仪垂向线运动测量误差的补偿方法,如图1所示,包括以下步骤:
步骤1、对加速度计组合信号在飞机螺旋桨旋转频率上进行正余弦解调,得到两路信号
Figure BDA0002078908370000051
Figure BDA0002078908370000052
其中
Figure BDA0002078908370000053
为加速度计组合信号在飞机螺旋桨旋转频率上的正弦分量,
Figure BDA0002078908370000054
为加速度计组合信号在飞机螺旋桨旋转频率上的余弦分量;
步骤2、对垂向加速度计输出信号在飞机螺旋桨旋转频率上上进行正余弦解调,得到两路信号
Figure BDA0002078908370000055
Figure BDA0002078908370000056
其中
Figure BDA0002078908370000057
为垂向加速度在飞机螺旋桨旋转频率上的正弦分量,
Figure BDA0002078908370000058
为垂向加速度在飞机螺旋桨旋转频率上的余弦分量;
步骤3、将由从步骤1的加速度计组合信号得到的
Figure BDA0002078908370000059
Figure BDA00020789083700000510
信号与从步骤2的垂向加速度计输出信号得到的
Figure BDA00020789083700000511
Figure BDA00020789083700000512
信号进行相除处理,得到两组加速度计敏感轴垂向安装误差实时信号Δβs和Δβc,将两组信号Δβs和Δβc求平均得到更加真实的Δβ实时信号。其中Δβs是由飞机螺旋桨旋转频率正弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβc是由飞机螺旋桨旋转频率余弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβ是两个安装误差实时信号的平均值;
步骤4、对垂向加速度计输出信号在二倍频频点上进行正余弦解调,得到两路信号
Figure BDA0002078908370000061
Figure BDA0002078908370000062
Figure BDA0002078908370000063
是垂向加速度在旋转二倍频上的正弦分量,
Figure BDA0002078908370000064
是垂向加速度在旋转二倍频上的余弦分量;
步骤5、将由从垂向加速度计输出信号得到的
Figure BDA0002078908370000065
Figure BDA0002078908370000066
信号与Δβ信号进行相乘处理,得到重力梯度信号中加速度计敏感轴俯仰角安装误差ΔΓuv(az)和ΔΓxy(az),予以补偿。
需要强调的是,本发明所述实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。

Claims (1)

1.一种基于载体自振动航空重力梯度仪运动误差补偿方法,其特征在于:包括以下步骤:
步骤1、对加速度计组合信号在飞机螺旋桨旋转频率上进行正余弦解调,得到两路信号
Figure FDA0002688302670000011
Figure FDA0002688302670000012
其中
Figure FDA0002688302670000013
为加速度计组合信号在飞机螺旋桨旋转频率上的正弦分量,
Figure FDA0002688302670000014
为加速度计组合信号在飞机螺旋桨旋转频率上的余弦分量;
步骤2、对垂向加速度计输出信号在飞机螺旋桨旋转频率上上进行正余弦解调,得到两路信号
Figure FDA0002688302670000015
Figure FDA0002688302670000016
其中
Figure FDA0002688302670000017
为垂向加速度在飞机螺旋桨旋转频率上的正弦分量,
Figure FDA0002688302670000018
为垂向加速度在飞机螺旋桨旋转频率上的余弦分量;
步骤3、将由从步骤1的加速度计组合信号得到的
Figure FDA0002688302670000019
Figure FDA00026883026700000110
信号与从步骤2的垂向加速度计输出信号得到的
Figure FDA00026883026700000111
Figure FDA00026883026700000112
信号进行相除处理,得到两组加速度计敏感轴垂向安装误差实时信号Δβs和Δβc,将两组信号Δβs和Δβc求平均得到更加真实的Δβ实时信号;其中Δβs是由飞机螺旋桨旋转频率正弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβc是由飞机螺旋桨旋转频率余弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβ是两个安装误差实时信号的平均值;
步骤4、对垂向加速度计输出信号在二倍频频点上进行正余弦解调,得到两路信号
Figure FDA00026883026700000113
Figure FDA00026883026700000114
Figure FDA00026883026700000115
是垂向加速度在旋转二倍频上的正弦分量,
Figure FDA00026883026700000116
是垂向加速度在旋转二倍频上的余弦分量;
步骤5、将由从垂向加速度计输出信号得到的
Figure FDA0002688302670000021
Figure FDA0002688302670000022
信号与Δβ信号进行相乘处理,得到重力梯度信号中加速度计敏感轴俯仰角安装误差ΔΓuv(az)和ΔΓxy(az),予以补偿。
CN201910464038.7A 2019-05-30 2019-05-30 基于载体自振动航空重力梯度仪运动误差补偿方法 Active CN110068876B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910464038.7A CN110068876B (zh) 2019-05-30 2019-05-30 基于载体自振动航空重力梯度仪运动误差补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910464038.7A CN110068876B (zh) 2019-05-30 2019-05-30 基于载体自振动航空重力梯度仪运动误差补偿方法

Publications (2)

Publication Number Publication Date
CN110068876A CN110068876A (zh) 2019-07-30
CN110068876B true CN110068876B (zh) 2021-01-26

Family

ID=67372135

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910464038.7A Active CN110068876B (zh) 2019-05-30 2019-05-30 基于载体自振动航空重力梯度仪运动误差补偿方法

Country Status (1)

Country Link
CN (1) CN110068876B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111650664B (zh) * 2020-06-30 2022-08-26 东南大学 一种航空重力梯度仪实时重力梯度解调方法及装置
CN113433600B (zh) * 2021-06-23 2022-04-12 中国船舶重工集团公司第七0七研究所 一种重力仪安装误差角标定方法
CN113447994B (zh) * 2021-06-25 2022-04-12 中国船舶重工集团公司第七0七研究所 一种基于旋转加速度计原理的重力梯度全张量测量方法
CN113885098B (zh) * 2021-09-07 2023-04-28 中国船舶重工集团公司第七0七研究所 一种重力敏感器低频频率响应误差在线建模及补偿方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2612677A1 (en) * 2006-11-23 2008-05-23 Technological Resources Pty. Limited Compensation for unwanted accelerations in a gravity gradiometer
DE102008054739A1 (de) * 2008-12-16 2010-06-17 Robert Bosch Gmbh Navigationsverfahren
CN103743395A (zh) * 2014-01-17 2014-04-23 哈尔滨工程大学 一种惯性重力匹配组合导航系统中时间延迟的补偿方法
CN105204083A (zh) * 2015-09-14 2015-12-30 武汉大学 基于绝对重力仪的垂直重力梯度提取方法及系统
CN108931824A (zh) * 2018-04-27 2018-12-04 东南大学 一种动基座旋转加速度计重力梯度仪误差增益系数标定方法
CN109470241A (zh) * 2018-11-23 2019-03-15 中国船舶重工集团公司第七0七研究所 一种具备重力扰动自主补偿功能的惯性导航系统及方法
CN109709628A (zh) * 2019-02-15 2019-05-03 东南大学 一种旋转加速度计重力梯度仪标定方法
CN109766812A (zh) * 2018-12-31 2019-05-17 东南大学 一种旋转加速度计重力梯度仪运动误差事后补偿方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AUPS114702A0 (en) * 2002-03-18 2002-04-18 Bhp Billiton Innovation Pty Ltd Enhancement of sensors for airborne operation

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2612677A1 (en) * 2006-11-23 2008-05-23 Technological Resources Pty. Limited Compensation for unwanted accelerations in a gravity gradiometer
DE102008054739A1 (de) * 2008-12-16 2010-06-17 Robert Bosch Gmbh Navigationsverfahren
CN103743395A (zh) * 2014-01-17 2014-04-23 哈尔滨工程大学 一种惯性重力匹配组合导航系统中时间延迟的补偿方法
CN105204083A (zh) * 2015-09-14 2015-12-30 武汉大学 基于绝对重力仪的垂直重力梯度提取方法及系统
CN108931824A (zh) * 2018-04-27 2018-12-04 东南大学 一种动基座旋转加速度计重力梯度仪误差增益系数标定方法
CN109470241A (zh) * 2018-11-23 2019-03-15 中国船舶重工集团公司第七0七研究所 一种具备重力扰动自主补偿功能的惯性导航系统及方法
CN109766812A (zh) * 2018-12-31 2019-05-17 东南大学 一种旋转加速度计重力梯度仪运动误差事后补偿方法
CN109709628A (zh) * 2019-02-15 2019-05-03 东南大学 一种旋转加速度计重力梯度仪标定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于旋转加速度计原理的重力梯度测量技术研究与试验;杨晔 等;《导航定位与授时》;20170731;第4卷(第4期);第20-28页 *
旋转加速度计重力梯度仪重力梯度信号仿真;蔡体菁 等;《物探与化探》;20151231;第39卷;第76-79页 *

Also Published As

Publication number Publication date
CN110068876A (zh) 2019-07-30

Similar Documents

Publication Publication Date Title
CN110068876B (zh) 基于载体自振动航空重力梯度仪运动误差补偿方法
CN110207697B (zh) 基于角加速度计/陀螺/加速度计的惯性导航解算方法
Crawford et al. A sensitive fast-response probe to measure turbulence and heat flux from any airplane
CN102023051B (zh) 一种测量星上有效载荷三轴角位移高频微振动的方法
JP5474938B2 (ja) 加速度計および誤差補償方法
CN108759845A (zh) 一种基于低成本多传感器组合导航的优化方法
JPH06288771A (ja) 飛行機の慣性航行計器の校正方法
CN201561759U (zh) 惯性姿态方位测量装置
US10654564B2 (en) Aircraft assembly including deflection sensor
CN102620892A (zh) 一种转动部件的动平衡测试方法
CN103323625A (zh) 一种mems-imu中加速度计动态环境下的误差标定补偿方法
CN102322873A (zh) 一种分布式pos地面演示验证系统
CN109766812A (zh) 一种旋转加速度计重力梯度仪运动误差事后补偿方法
JP6629691B2 (ja) センサパッケージおよび自動運転車両
CN111650664B (zh) 一种航空重力梯度仪实时重力梯度解调方法及装置
Rodi et al. Correction of static pressure on a research aircraft in accelerated flight using differential pressure measurements
CN102072729A (zh) 用于测量旋转飞行器姿态的测量装置及其测量方法
CN106885587A (zh) 旋翼扰动下惯性/gps组合导航外杆臂效应误差补偿方法
KR101658473B1 (ko) Mems자이로스코프의 가속도 민감도 보정 방법
US6810739B1 (en) Accelerometer augmented leveling device
CN109931959B (zh) 硅微陀螺仪正交误差校正方法
US7461548B2 (en) Method and a device for measuring the speed of an aircraft, in particular a rotorcraft at low speed
Pérez-D'Arpino et al. Development of a low cost inertial measurement unit for uav applications with kalman filter based attitude determination
CN114370930A (zh) 利用mhd微角振动传感器测量卫星结构微角振动的方法
EA037452B1 (ru) Гравиметрический узел

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