一种航空重力梯度仪实时重力梯度解调方法及装置制造方法及图纸

本发明专利技术公开了一种航空重力梯度仪实时重力梯度解调方法及装置，属于精密测量领域。方法包括：首先对四只加速度计输出信号进行信号调理和滤波；然后A/D转换器根据角度信息转换器发出的同步脉冲信号对滤波后的重力梯度模拟信号进行A/D转换，得到重力梯度数字信号；接下来利用参考信号源对重力梯度数字信号进行相位角解调，得到解调相位角；最后利用解调相位角、角度信息和参考信号源对重力梯度数字序列进行实时重力梯度解调，得到具有时空信息的重力梯度数据。使用本发明专利技术提供的航空重力梯度仪实时重力梯度解调方法，能够进一步提高重力梯度测量精度，操作简单方便，同时可为后续的离线处理和地质分析提供丰富的数据材料。

【技术实现步骤摘要】
一种航空重力梯度仪实时重力梯度解调方法及装置
本专利技术属于精密测量
，具体涉及一种航空重力梯度仪实时重力梯度解调方法及装置。
技术介绍
高精度航空重力梯度测量技术对于空间科学、地球科学、地质科学、能源勘探、惯性导航等领域的研究具有重要意义。目前，重力梯度仪主要有旋转加速度计重力梯度仪、超导重力梯度仪、冷原子重力梯度仪、静电重力梯度仪，以及基于微机械加工(MEMS)技术的重力梯度仪等，其中旋转加速度计重力梯度仪是目前唯一能够商用，并且成功用在机载/船载动基座上的重力梯度测量仪器。航空重力梯度测量以其经济、高效及对小尺度地质目标特有的高分辨率而备受青睐。旋转加速度计重力梯度仪是一种高精密测量仪器，电子元件稳定性、电子线路电磁兼容性问题、敏感器工作性能、旋转机构稳定性、温控性能等因素都会对重力梯度测量精度造成严重影响，这些因素会使重力梯度仪输出信号包含各种噪声，从而降低了测量精度。现有技术中为提高重力梯度测量精度做了很多工作，取得了一些成果，例如，专利CN108873093B公开了一种航空重力梯度仪自梯度补偿方法，通过对运载体运行过程中由于姿态改变引起的自梯度进行补偿，提高重力梯度测量精度；专利CN110471123A公开了一种旋转加速度计重力梯度仪数据诊断及处理方法，通过对重力梯度数据进行合理的数据诊断和处理，提高测量精度；专利CN109212619B公开了一种旋转加速度计重力梯度仪线运动误差补偿装置及方法，专利CN109212629B公开了一种旋转加速度计重力梯度仪角运动误差补偿装置及方法，专利CN109212620B公开了一种动基座旋转加速度计重力梯度仪误差补偿装置及方法，专利CN110068876A公开了一种基于载体自振动航空重力梯度仪运动误差补偿方法，这些误差补偿方法对重力梯度仪工作过程中的角运动误差、线运动误差、自梯度、垂向运动误差进行研究并通过特定的方法予以补偿，提高了机器测量性能。然而，除了以上影响因素外，重力梯度解调相位角的精确确定也是一个重要的影响因子，其直接影响了重力梯度测量性能。解调相位角的误差主要来自旋转机构稳定性、电机控制稳定性、电路噪声以及信号处理导致的传输延时等方面，但目前针对航空重力梯度仪解调相位角的误差补偿并没有相关文献发布或发表。
技术实现思路
专利技术目的：为了克服现有技术的不足，本专利技术提供一种操作简单、实施方便，并且具有实用性强的一种航空重力梯度仪实时重力梯度解调方法及装置，解决了重力梯度解调相位角不精确而导致的重力梯度测量性能下降问题，能够为高精度航空重力梯度测量提供重要的应用参考。技术方案：为了实现以上目的，本专利技术采用以下技术方案：第一方面，提供一种航空重力梯度仪实时重力梯度解调方法，包括以下步骤：对四只加速度计输出的信号进行信号调理和滤波，得到重力梯度模拟信号；A/D转换器根据角度信息转换器发出的同步脉冲信号对滤波后的重力梯度模拟信号进行A/D转换，得到重力梯度数字信号；利用参考信号源对重力梯度数字信号进行相位角解调，得到解调相位角；利用解调相位角、角度信息转换器输出的角度信息和参考信号源对重力梯度数字序列进行实时重力梯度解调，得到具有时空信息的重力梯度数据。其中，所述重力梯度数字信号表达形式为：Eout(k)＝2KIKR{(Γyy-Γxx)sin2[ωΔt(k-1)-α]+2Γxycos2[ωΔt(k-1)-α]}式中，k为数据序列索引，K为系统信号放大增益，Δt为采样时间，KI为加速度计标度因数，R为圆盘半径，ω为圆盘旋转角频率，α为重力梯度解调相位角，Eout(k)为第(k-1)Δt时刻数字信号，Гxx为x轴上的重力加速度分量在x轴方向上的空间导数，Гyy为y轴上的重力加速度分量在y轴方向上的空间导数，(Гyy-Гxx)为重力梯度Гyy与Гxx之差，Гxy为x(y)轴上的重力加速度分量在y(x)轴方向上的空间导数，所述x轴和y轴所在重力梯度仪坐标系OXYZ是东北天地理坐标系，原点O为重力梯度仪测量中心。进一步地，所述利用参考信号源对重力梯度数字信号进行相位角解调包括：发出数字正弦信号Ss(k,i)和数字余弦信号Sc(k,i)作为参考信号；利用参考信号对A/D转换器的输出信号进行相位角解调，计算公式如下：式中，ζi分别为正弦、余弦信号参考源的初始相角在第i次的修正角，等式右侧|max表示寻找最大值对应的正弦、余弦信号参考源的初始相角，N为圆盘每旋转一圈角度信息转换器输出的脉冲数，ζq分别为第q次得到的重力梯度解调相位角。作为优选的实施方式，得到重力梯度解调相位角和ζq后，对和ζq进行组合，把组合后的数据通过零相位低通滤波器进行滤波处理，对处理后的相位角进行均值处理，得到最终的重力梯度解调相位角α。所述重力梯度解调表达式为：式中，p为当前时刻的圆盘旋转圈数序列索引，M为圆盘旋转总圈数，gps为GPS提供的位置信息，θp,k为角度信息转换器实时输出的圆盘角位置信息。第二方面，提供一种航空重力梯度仪实时重力梯度解调装置，包括：信号调理模块，用于对四只加速度计输出的信号进行信号调理和滤波，得到重力梯度模拟信号；A/D转换模块，用于根据角度信息模块发出的同步脉冲信号对滤波后的重力梯度模拟信号进行A/D转换，得到重力梯度数字信号；角度信息模块，用于发出同步脉冲信号以及实时输出圆盘角度信息；相位角解调模块，利用参考信号源对重力梯度数字信号进行相位角解调，得到解调相位角；以及重力梯度解调模块，利用解调相位角、角度信息和参考信号源对重力梯度数字序列进行实时重力梯度解调，得到具有时空信息的重力梯度数据。有益效果：本专利技术从实际测量系统出发，首次给出了航空重力梯度仪实时重力梯度解调的处理方法，可以实现重力梯度解调相位角补偿和高精度重力梯度解调，提高重力梯度测量精度，且该方法具有操作简便、便于实施、应用性强等特点，对提高动基座航空重力梯度测量性能具有重要参考价值。附图说明图1为本专利技术的航空重力梯度仪实时重力梯度解调方法的工作原理示意图；图2为根据本专利技术实施例的重力梯度解调相位补偿前后的重力梯度信息(Гyy-Гxx)比较的波形图；图3为根据本专利技术实施例的重力梯度解调相位补偿前后的重力梯度信息Гxy比较的波形图。具体实施方式下面结合附图对本专利技术的技术方案作进一步说明。图1为本专利技术航空重力梯度仪实时重力梯度解调方法的原理示意图，结合图1，对本专利技术的航空重力梯度仪实时重力梯度解调方法进行说明，该方法包括以下步骤：1)重力梯度仪放置在稳定平台上，以重力梯度仪测量中心为原点O，东北天地理坐标系为重力梯度仪坐标系，记为oxyz；2)四只加速度计输出信号经过信号调理和滤波电路进行处理，得到模拟信号Eout(t)，角度信息转换器输出同步脉冲信号给A/D转换器，A/D转换器对模拟信号Eout(t)进行A/D转换得到第(k-1)Δt时刻数字本文档来自技高网...

【技术保护点】
1.一种航空重力梯度仪实时重力梯度解调方法，其特征在于，所述方法包括以下步骤：/n对四只加速度计输出的信号进行信号调理和滤波，得到重力梯度模拟信号；/nA/D转换器根据角度信息转换器发出的同步脉冲信号对滤波后的重力梯度模拟信号进行A/D转换，得到重力梯度数字信号；/n利用参考信号源对重力梯度数字信号进行相位角解调，得到解调相位角；/n利用解调相位角、角度信息转换器输出的角度信息和参考信号源对重力梯度数字序列进行实时重力梯度解调，得到具有时空信息的重力梯度数据。/n

【技术特征摘要】
1.一种航空重力梯度仪实时重力梯度解调方法，其特征在于，所述方法包括以下步骤：
对四只加速度计输出的信号进行信号调理和滤波，得到重力梯度模拟信号；
A/D转换器根据角度信息转换器发出的同步脉冲信号对滤波后的重力梯度模拟信号进行A/D转换，得到重力梯度数字信号；
利用参考信号源对重力梯度数字信号进行相位角解调，得到解调相位角；
利用解调相位角、角度信息转换器输出的角度信息和参考信号源对重力梯度数字序列进行实时重力梯度解调，得到具有时空信息的重力梯度数据。


2.根据权利要求1所述的航空重力梯度仪实时重力梯度解调方法，其特征在于，所述A/D转换器得到的重力梯度数字信号表达形式为：
Eout(k)＝2KIKR{(Γyy-Γxx)sin2[ωΔt(k-1)-α]+2Γxycos2[ωΔt(k-1)-α]}
式中，k为数据序列索引，K为系统信号放大增益，Δt为采样时间，KI为加速度计标度因数，R为圆盘半径，ω为圆盘旋转角频率，α为重力梯度解调相位角，Eout(k)为第(k-1)Δt时刻数字信号，Гxx为x轴上的重力加速度分量在x轴方向上的空间导数，Гyy为y轴上的重力加速度分量在y轴方向上的空间导数，(Гyy-Гxx)为重力梯度Гyy与Гxx之差，Гxy为x(y)轴上的重力加速度分量在y(x)轴方向上的空间导数，所述x轴和y轴所在重力梯度仪坐标系OXYZ是东北天地理坐标系，原点O为重力梯度仪测量中心。


3.根据权利要求1所述的航空重力梯度仪实时重力梯度解调方法，其特征在于，所述利用参考信号源对重力梯度数字信号进行相位角解调包括：
发出数字正弦信号Ss(k,i)和数字余弦信号Sc(k,i)作为参考信号；
利用参考信号对A/D转换器的输出信号进行相位角解调，计算公式如下：



式中，ζi分别为正弦、余弦信号参考源的初始相角在第i次的修正角，等式右侧|max表示寻找最大值对应的正弦、余弦信号参考源的初始相角，N为圆盘每旋转一圈角度信息转换器输出的脉冲数，ζq分别为第q次得到的重力梯度解调相位角。


4.根据权利要求3所述的航空重力梯度仪实时重力梯度解调方法，其特征在于，所述方法还包括：对得到的重力梯度解调相位角和ζq进行组合，把组合后的数据通过零相位低通滤波器进行滤波处理，对处理后的相位角进行均值处理，得到最终的重力梯度解调相位角，作为重力梯度解调的解调相位角α。


5.根据权利要求2所述的航空重力梯度仪实时重力梯度解调方法，其特征在于，所述重力梯度解调表达式为：



式中，p为当前时刻的圆盘旋转圈数序列索引，M为圆盘旋转总圈数，gps为GPS提供的位置信息，θp,k为角度信息转换器实时输出的圆盘角位置信息。


6.一种航空重力梯度仪实时重力梯度解调装置，其特...

【专利技术属性】
技术研发人员：钱学武，赵立业，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏;32