CN108873093B - 一种航空重力梯度仪自梯度补偿方法 - Google Patents
一种航空重力梯度仪自梯度补偿方法 Download PDFInfo
- Publication number
- CN108873093B CN108873093B CN201810761713.8A CN201810761713A CN108873093B CN 108873093 B CN108873093 B CN 108873093B CN 201810761713 A CN201810761713 A CN 201810761713A CN 108873093 B CN108873093 B CN 108873093B
- Authority
- CN
- China
- Prior art keywords
- carrier
- under
- gravity
- gradient
- attitude angle
- 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
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
- G01V7/16—Measuring gravitational fields or waves; Gravimetric prospecting or detecting specially adapted for use on moving platforms, e.g. ship, aircraft
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Manufacturing & Machinery (AREA)
- Aviation & Aerospace Engineering (AREA)
- Ocean & Marine Engineering (AREA)
- Testing Of Balance (AREA)
Abstract
本发明公开了一种航空重力梯度仪自梯度补偿方法,首先把重力梯度仪放置在高精度稳定平台上,调整运载体姿态并实时记录运载体姿态、燃油质量以及在该状态下对应的重力梯度测量信息;然后将初始重力梯度值与不同状态下的重力梯度测量值相减得到相应状态下的自梯度值;最后对含有运载体姿态和燃油质量信息的航空重力梯度测量值与相同状态下的自梯度值进行相减,得到重力梯度仪真实测量到的重力梯度信息。本发明可以实现航空重力梯度的自梯度补偿,能够进一步提高航空重力梯度的测量精度。
Description
技术领域
本发明涉及一种航空重力梯度仪自梯度补偿方法,尤其涉及一种能够显著提高航空重力梯度测量精度的补偿方法,属精密测量技术领域。
背景技术
高精度重力梯度测量对能源勘探、空间科学、地球科学、地质科学的发展以及在惯性导航等方面具有非常重要的应用价值和应用前景。重力梯度测量相对于重力测量有明显的优势,对于地质结构的形状、尺寸和埋藏深度等信息,重力梯度测量能够提供更高的解析度,特别适合于探测和研究局部地质结构及其细节。
自梯度是由重力梯度仪运载体(如:船舶、飞机、飞艇等)姿态改变导致其分布质量相对于重力梯度仪的空间存在位置变化而产生的梯度变化。在进行航空重力梯度测量时,重力梯度仪被放置在高精度惯性稳定平台上,重力梯度仪相对于惯性空间保持方位不变,但是运载体在运行过程中会受到气流、风向等因素的影响,不可避免地引起运载体姿态的改变,运载体姿态的改变会直接导致运载体分布质量相对于重力梯度仪在空间上的位置变化;同时,运载体在运行过程中,燃油质量逐渐减小,变化的燃油质量相对于重力梯度仪所处空间引起引力变化。因此,为了提高航空重力梯度测量精度,必须对航空重力梯度仪进行自梯度补偿。
发明内容
技术问题:本发明提供一种简单、高效、实现方便,具有重要实用价值的航空重力梯度仪自梯度补偿方法。
技术方案:本发明的一种航空重力梯度仪自梯度补偿方法,包括以下步骤:
1)以重力梯度仪中心为原点,东北天方向为坐标轴建立重力梯度仪坐标系oxgygzg和运载体坐标系oxbybzb,所述重力梯度仪坐标系oxgygzg表示为g系,运载体坐标系oxbybzb表示为b系;初始时刻g系与b系重合;
2)调整运载体姿态,采集并记录运载体姿态信息、燃油质量信息和重力梯度仪输出信息,然后进行自梯度标定;
3)根据下式计算g系到b系的旋转矩阵
其中,ψ、λ和分别为运载体的航向角、俯仰角和横滚角;
根据下式计算在g系下运载体姿态变化前后的重力梯度张量元素变化量,即运载体姿态变化引起的自梯度分量:
其中,为在g系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在x轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在y轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在z轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在y轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在z轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在z轴方向上的重力加速度分量在z轴方向上的空间导数变化量,I6×6为6行6列单位矩阵,为在b系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在x轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在y轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在z轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在y轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在z轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在z轴方向上的重力加速度分量在z轴方向上的空间导数, 为在b系下运载体姿态角为ψ、λ和时重力梯度仪输出的重力梯度分量信息,为重力梯度向量矩阵;
4)根据下式计算在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量:
其中,t为时间,To为燃油总消耗时间,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γxx,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γxy,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γxz,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γyy,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γyz,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γzz,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γxx,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γxy,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γxz,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γyy,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γyz,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γzz, 为在b系下燃油无消耗时运载体姿态角为ψ、λ和时运载体姿态改变引起的自梯度分量;
5)根据下式计算在g系下重力梯度仪测量到的真实重力梯度分量:
其中,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γxx,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γxy,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γxz,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γyy,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γyz,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γzz,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γxx,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γxy,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γxz,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γyy,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γyz,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γzz。
进一步的,本发明方法中,所述步骤2)中的运载体姿态信息为俯仰角、航向角和横滚角信息,由重力梯度仪稳定平台的角传感器获得。
进一步的,本发明方法中,所述步骤2)的内容在地面上进行,所述自梯度标定是在采集并记录运载体姿态信息和燃油质量信息过程中,根据状态调整前后重力梯度仪输出的重力梯度变化量确定不同状态与重力梯度仪实际输出之间的对应关系。
进一步的,本发明方法中,所述步骤3)中的重力梯度向量矩阵根据下式计算:
其中,C11为中的第1行、第1列元素,C12为中的第1行、第2列元素,C13为中的第1行、第3列元素,C21为中的第2行、第1列元素,C22为中的第2行、第2列元素,C23为中的第2行、第3列元素,C31为中的第3行、第1列元素,C32为中的第3行、第2列元素,C33为中的第3行、第3列元素。
本发明中,在地面上确定出运载体姿态、燃油质量与重力梯度仪输出之间的对应关系,根据该对应关系实时对重力梯度仪测量到的重力梯度信息进行补偿处理,实现自梯度补偿,从而提高航空重力梯度的测量精度。
有益效果:本发明与现有技术相比,具有以下优点:
目前,有关航空重力梯度测量过程中的自梯度补偿的具体方法很少有相关文献公布或发表。本发明给出了航空重力梯度仪的自梯度补偿方法,综合考虑了运载体姿态和燃油消耗对重力梯度的测量影响。首先通过调整运载体姿态和燃油质量,确定出运载体姿态、燃油消耗与航空重力梯度仪输出信息之间的对应关系;然后在进行航空重力梯度测量过程中,根据该对应关系实时对重力梯度信息进行自梯度补偿,最后获得真实重力梯度信息。采用所提方法可以对重力梯度仪输出信息进行实时在线重力梯度补偿,能够显著提高航空重力梯度测量精度,且该方法具有操作简便、易于实施、应用性强等特点,对提高我国航空重力测量技术具有重要参考价值。
附图说明
图1为全张量航空重力梯度仪自梯度补偿处理流程示意图。
图2为运载体航向角与自梯度分量间的关系。
图3为运载体俯仰角与自梯度分量间的关系。
图4为运载体横滚角与自梯度分量间的关系。
图5a为运载体航向角和燃油消耗引起的自梯度(Γxx分量)。
图5b为运载体航向角和燃油消耗引起的自梯度(Γxy分量)。
图5c为运载体航向角和燃油消耗引起的自梯度(Γxz分量)。
图5d为运载体航向角和燃油消耗引起的自梯度(Γyy分量)。
图5e为运载体航向角和燃油消耗引起的自梯度(Γyz分量)。
图5f为运载体航向角和燃油消耗引起的自梯度(Γzz分量)。
图6a为运载体俯仰角和燃油消耗引起的自梯度(Γxx分量)。
图6b为运载体俯仰角和燃油消耗引起的自梯度(Γxy分量)。
图6c为运载体俯仰角和燃油消耗引起的自梯度(Γxz分量)。
图6d为运载体俯仰角和燃油消耗引起的自梯度(Γyy分量)。
图6e为运载体俯仰角和燃油消耗引起的自梯度(Γyz分量)。
图6f为运载体俯仰角和燃油消耗引起的自梯度(Γzz分量)。
图7a为运载体横滚角和燃油消耗引起的自梯度(Γxx分量)。
图7b为运载体横滚角和燃油消耗引起的自梯度(Γxy分量)。
图7c为运载体横滚角和燃油消耗引起的自梯度(Γxz分量)。
图7d为运载体横滚角和燃油消耗引起的自梯度(Γyy分量)。
图7e为运载体横滚角和燃油消耗引起的自梯度(Γyz分量)。
图7f为运载体横滚角和燃油消耗引起的自梯度(Γzz分量)。
图8为重力梯度分量补偿误差。
具体实施方式
下面结合实施例和说明书附图对本发明作进一步地说明。
图1为全张量航空重力梯度仪自梯度补偿处理流程示意图,具体实施步骤如下:
1)以重力梯度仪中心为原点,东北天方向为坐标轴建立重力梯度仪坐标系oxgygzg和运载体坐标系oxbybzb,所述重力梯度仪坐标系oxgygzg表示为g系,运载体坐标系oxbybzb表示为b系;初始时刻g系与b系重合;
2)调整运载体姿态,采集并记录运载体姿态信息、燃油质量信息和重力梯度仪输出信息,然后进行自梯度标定;
3)根据下式计算g系到b系的旋转矩阵
其中,ψ、λ和分别为运载体的航向角、俯仰角和横滚角;
根据下式计算在g系下运载体姿态变化前后的重力梯度张量元素变化量,即运载体姿态变化引起的自梯度分量:
其中,为在g系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在x轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在y轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在z轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在y轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在z轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在z轴方向上的重力加速度分量在z轴方向上的空间导数变化量,I6×6为6行6列单位矩阵,为在b系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在x轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在y轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在z轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在y轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在z轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在z轴方向上的重力加速度分量在z轴方向上的空间导数, 为在b系下运载体姿态角为ψ、λ和时重力梯度仪输出的重力梯度分量信息,为重力梯度向量矩阵;
4)根据下式计算在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量:
其中,t为时间,To为燃油总消耗时间,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γxx,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γxy,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γxz,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γyy,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γyz,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γzz,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γxx,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γxy,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γxz,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γyy,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γyz,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γzz, 为在b系下燃油无消耗时运载体姿态角为ψ、λ和时运载体姿态改变引起的自梯度分量;
5)根据下式计算在g系下重力梯度仪测量到的真实重力梯度分量:
其中,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γxx,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γxy,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γxz,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γyy,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γyz,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γzz,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γxx,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γxy,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γxz,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γyy,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γyz,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γzz。
本发明方法的一种实施例中,所述步骤2)中的运载体姿态信息为俯仰角、航向角和横滚角信息,由重力梯度仪稳定平台的角传感器获得。
本发明方法的一种实施例中,所述步骤2)的内容在地面上进行,所述自梯度标定是在采集并记录运载体姿态信息和燃油质量信息过程中,根据状态调整前后重力梯度仪输出的重力梯度变化量确定不同状态与重力梯度仪实际输出之间的对应关系。
本发明方法的一种实施例中,所述步骤3)中的重力梯度向量矩阵根据下式计算:
其中,C11为中的第1行、第1列元素,C12为中的第1行、第2列元素,C13为中的第1行、第3列元素,C21为中的第2行、第1列元素,C22为中的第2行、第2列元素,C23为中的第2行、第3列元素,C31为中的第3行、第1列元素,C32为中的第3行、第2列元素,C33为中的第3行、第3列元素。
为了便于对本发明实施例中的方案进行验证,进行如下参数设定和仿真:
首先对运载体姿态引起的自梯度进行仿真,设定运载体姿态角度调整范围为0°~180°,角度变化量为0.001°,运载体分布质量、分布密度和在运载体坐标系下的中心位置坐标信息如下表:
质量分布 | 质量分布体积(m×m×m) | 分布密度(kg/m<sup>3</sup>) | 位置坐标(m) |
前部 | 1.0×1.0×0.8 | 1200 | (0,2.5,0.5) |
中部 | 1.0×1.0×0.5 | 1000 | (0,0.5,0.5) |
左翼 | 2.0×1.0×0.3 | 500 | (-4.0,0.5,1) |
右翼 | 2.0×1.0×0.3 | 500 | (4.0,0.5,1) |
后部 | 1.8×1.2×0.5 | 800 | (0,-3.0,0.5) |
在该参数下进行自梯度仿真,运载体在不同姿态下引起的自梯度如图2、图3、图4所示,从图中可以看出,运载体姿态引起的自梯度最大达到了10Eu。
然后,对运载体燃油消耗引起的自梯度进行仿真,设定燃油密度780kg/m3,初始时刻燃油体积为1.8m×2.8m×0.6m,燃油消耗时间为6h,在该参数下进行自梯度仿真,不同姿态下对应的燃油消耗所引起的自梯度结果如图5a至图5f、图6a至图6f、图7a至图7f所示,从图中可以看出,燃油消耗引起的自梯度最大达到了30Eu。
最后,对所提出的自梯度补偿方法进行仿真,设定运载体运行高度为100m,梯度输出空间间距为50m,异常体密度为2000kg/m3,异常体体积为20000m×1000m×1000m,异常体中心位置坐标为(0,-2000,-1100)m,运载体运行的起始坐标为(-60000,0,0)m,终止坐标为(60000,0,0)m,在该参数下进行自梯度补偿效果仿真,仿真结果如图8所示,从图中可以看出,5个重力梯度分量误差均小于0.1Eu,能够满足航空重力梯度测量要求。
应理解上述实施例仅用于说明本发明技术方案的优选实施方式,而不用于限制本发明的范围。在阅读了本发明之后,本领域技术人员对本发明的各种等同形式的修改和替换均落于本申请权利要求所限定的保护范围。
Claims (3)
1.一种航空重力梯度仪自梯度补偿方法,其特征在于,该方法包括以下步骤:
1)以重力梯度仪中心为原点,东北天方向为坐标轴建立重力梯度仪坐标系oxgygzg和运载体坐标系oxbybzb,所述重力梯度仪坐标系oxgygzg表示为g系,运载体坐标系oxbybzb表示为b系;初始时刻g系与b系重合;
2)调整运载体姿态,采集并记录运载体姿态信息、燃油质量信息和重力梯度仪输出信息,然后进行自梯度标定;
3)根据下式计算g系到b系的旋转矩阵
其中,ψ、λ和分别为运载体的航向角、俯仰角和横滚角;
根据下式计算在g系下运载体姿态变化前后的重力梯度张量元素变化量,即运载体姿态变化引起的自梯度分量:
其中,为在g系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在x轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在y轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在z轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在y轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在z轴方向上的空间导数变化量,为在g系下运载体姿态角为ψ、λ和时对应的在z轴方向上的重力加速度分量在z轴方向上的空间导数变化量,I6×6为6行6列单位矩阵,为在b系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在x轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在y轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在x轴方向上的重力加速度分量在z轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在y轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在y轴方向上的重力加速度分量在z轴方向上的空间导数,为在b系下运载体姿态角为ψ、λ和时对应的在z轴方向上的重力加速度分量在z轴方向上的空间导数, 为在b系下运载体姿态角为ψ、λ和时重力梯度仪输出的重力梯度分量信息,为重力梯度向量矩阵,所述重力梯度向量矩阵根据下式计算:
其中,C11为中的第1行、第1列元素,C12为中的第1行、第2列元素,C13为中的第1行、第3列元素,C21为中的第2行、第1列元素,C22为中的第2行、第2列元素,C23为中的第2行、第3列元素,C31为中的第3行、第1列元素,C32为中的第3行、第2列元素,C33为中的第3行、第3列元素;
4)根据下式计算在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量:
其中,t为时间,To为燃油总消耗时间,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γxx,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γxy,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γxz,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γyy,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γyz,为在t时刻,在g系下运载体姿态角为ψ、λ和时对应的燃油消耗引起的自梯度分量Γzz,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γxx,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γxy,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γxz,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γyy,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γyz,为在b系下燃油无消耗、运载体姿态角为ψ、λ和时运载体姿态变化引起的自梯度分量Γzz, 为在b系下燃油无消耗时运载体姿态角为ψ、λ和时运载体姿态改变引起的自梯度分量;
5)根据下式计算在g系下重力梯度仪测量到的真实重力梯度分量:
其中,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γxx,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γxy,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γxz,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γyy,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γyz,为在t时刻,在g系下运载体姿态角为ψ、λ和时重力梯度仪测量到的真实重力梯度分量Γzz,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γxx,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γxy,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γxz,)为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γyy,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γyz,为在t时刻,在b系下运载体姿态角为ψ、λ和时重力梯度仪实际输出的重力梯度分量Γzz。
2.根据权利要求1所述的航空重力梯度仪自梯度补偿方法,其特征在于,所述步骤2)中的运载体姿态信息为俯仰角、航向角和横滚角信息,由重力梯度仪稳定平台的角传感器获得。
3.根据权利要求1所述的航空重力梯度仪自梯度补偿方法,其特征在于,所述步骤2)的内容在地面上进行,所述自梯度标定是在采集并记录运载体姿态信息和燃油质量信息过程中,根据状态调整前后重力梯度仪输出的重力梯度变化量确定不同状态与重力梯度仪实际输出之间的对应关系。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810761713.8A CN108873093B (zh) | 2018-07-12 | 2018-07-12 | 一种航空重力梯度仪自梯度补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810761713.8A CN108873093B (zh) | 2018-07-12 | 2018-07-12 | 一种航空重力梯度仪自梯度补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108873093A CN108873093A (zh) | 2018-11-23 |
CN108873093B true CN108873093B (zh) | 2019-11-05 |
Family
ID=64301156
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810761713.8A Active CN108873093B (zh) | 2018-07-12 | 2018-07-12 | 一种航空重力梯度仪自梯度补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108873093B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109766812B (zh) | 2018-12-31 | 2021-03-09 | 东南大学 | 一种旋转加速度计重力梯度仪运动误差事后补偿方法 |
CN109709628B (zh) * | 2019-02-15 | 2020-08-14 | 东南大学 | 一种旋转加速度计重力梯度仪标定方法 |
CN113885100B (zh) * | 2021-09-28 | 2023-06-20 | 中国船舶重工集团公司第七0七研究所 | 一种重力梯度仪平台质量自梯度补偿方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106507770B (zh) * | 2012-04-05 | 2014-09-10 | 上海新跃仪表厂 | 一种用无拖曳技术实现重力梯度卫星姿轨控的方法 |
CN105044795A (zh) * | 2015-05-29 | 2015-11-11 | 中国船舶重工集团公司第七0七研究所 | 利用光纤陀螺减小旋转加速度计重力梯度仪测量误差的方法 |
CN105044798A (zh) * | 2015-06-29 | 2015-11-11 | 东南大学 | 旋转加速度计重力梯度仪加速度计标度因子反馈调整方法 |
KR101739390B1 (ko) * | 2015-12-11 | 2017-05-24 | 국방과학연구소 | 중력오차보상을 통한 관성항법장치의 자체정렬 정확도 향상기법 |
CN107576992A (zh) * | 2017-09-04 | 2018-01-12 | 东南大学 | 一种重力梯度仪自标定方法及离心梯度补偿方法 |
CN107870371A (zh) * | 2017-12-05 | 2018-04-03 | 东南大学 | 一种动基座重力梯度仪自梯度补偿方法 |
-
2018
- 2018-07-12 CN CN201810761713.8A patent/CN108873093B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106507770B (zh) * | 2012-04-05 | 2014-09-10 | 上海新跃仪表厂 | 一种用无拖曳技术实现重力梯度卫星姿轨控的方法 |
CN105044795A (zh) * | 2015-05-29 | 2015-11-11 | 中国船舶重工集团公司第七0七研究所 | 利用光纤陀螺减小旋转加速度计重力梯度仪测量误差的方法 |
CN105044798A (zh) * | 2015-06-29 | 2015-11-11 | 东南大学 | 旋转加速度计重力梯度仪加速度计标度因子反馈调整方法 |
KR101739390B1 (ko) * | 2015-12-11 | 2017-05-24 | 국방과학연구소 | 중력오차보상을 통한 관성항법장치의 자체정렬 정확도 향상기법 |
CN107576992A (zh) * | 2017-09-04 | 2018-01-12 | 东南大学 | 一种重力梯度仪自标定方法及离心梯度补偿方法 |
CN107870371A (zh) * | 2017-12-05 | 2018-04-03 | 东南大学 | 一种动基座重力梯度仪自梯度补偿方法 |
Non-Patent Citations (1)
Title |
---|
旋转式重力梯度仪自身梯度补偿研究;丁昊 等;《导航与控制》;20150430;第14卷(第2期);第29-32页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108873093A (zh) | 2018-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108873093B (zh) | 一种航空重力梯度仪自梯度补偿方法 | |
Bögel et al. | Test and calibration of the DLR Falcon wind measuring system by maneuvers | |
CN102736118B (zh) | 一种用于全球重力场测量的综合型卫星系统 | |
Cai et al. | The first airborne scalar gravimetry system based on SINS/DGPS in China | |
CN106353824A (zh) | 航空磁通门磁梯度张量仪的系统校正及磁干扰补偿融合方法 | |
US20160131794A1 (en) | Systems and methods for a gravity survey using a free-fall gravity sensor | |
CN106052694B (zh) | 基于重力矢量及其梯度张量对单个运动物体进行定位跟踪的方法 | |
CN106767671B (zh) | 基于三维电子罗盘的地质结构面产状计算方法 | |
CN109470241A (zh) | 一种具备重力扰动自主补偿功能的惯性导航系统及方法 | |
CN109709628A (zh) | 一种旋转加速度计重力梯度仪标定方法 | |
CN102906601B (zh) | 重力梯度仪 | |
Metzger et al. | Measuring the 3-D wind vector with a weight-shift microlight aircraft | |
Rodi et al. | Correction of static pressure on a research aircraft in accelerated flight using differential pressure measurements | |
CN109856691A (zh) | 一种基于梯度法的航空重力矢量向下延拓方法及系统 | |
CN110244367B (zh) | 一种基于地面多基站的ztem系统姿态补偿方法 | |
Rogers et al. | An investigation into the feasibility of using a modern gravity gradiometer instrument for passive aircraft navigation and terrain avoidance | |
CN104809352A (zh) | 基于正演的拖曳式航磁全张量梯度数据软补偿方法 | |
CN110736484A (zh) | 基于陀螺仪及磁传感器融合的背景磁场标定方法 | |
CN106290969B (zh) | 一种考虑减速伞气动力影响的风速风向探测方法 | |
Luo et al. | First unmanned aerial vehicle airborne gravimetry based on the CH-4 UAV in China | |
CN107804487A (zh) | 一种基于自适应偏差控制的跳跃式再入返回落点预报方法 | |
CN104297525A (zh) | 基于火箭橇试验的惯性测量系统加速度计标定方法 | |
Gnadt et al. | Signal enhancement for magnetic navigation challenge problem | |
Jekeli | Accuracy requirements in position and attitude for airborne vector gravimetry and gradiometry | |
CN206470416U (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 |