一种旋转加速度计重力梯度仪线运动误差补偿装置及方法与流程

本发明涉及一种旋转加速度计重力梯度仪线运动误差补偿装置及方法，属于精密测量技术领域。

背景技术：

动基座重力梯度勘探是一种低成本、高效率的重力梯度勘探方法；是目前世界上最先进的重力场勘探方式。重力梯度数据广泛应用于地质分析、重力场建模、高精度导航、资源勘探等。重力梯度仪具有极其重要的国防、民用价值。目前国内外在研的重力梯度仪主要有冷原子重力梯度仪、超导重力梯度仪、mems重力梯度仪等。国外已经投入商业应用的重力梯度仪主要有旋转加速度计重力梯度仪及旋转超导加速度计重力梯度仪。我国的重力梯度仪样机正处于研制中。

在动基座重力梯度勘探时，由于重力梯度仪内部的加速度计存在安装误差、加速度计一阶、高阶标度系数不匹配、电路增益不匹配等，导致重力梯度仪的加速度，角速度，角加速度传递到重力梯度仪的输出，造成测量误差。同时重力梯度仪线运动、自梯度、角运动会引起重力梯度仪前端信号调理电路过电压饱和或损坏。本发明提供一种能够实时补偿重力梯度仪线运动误差、自梯度的装置及方法，目前没有关于重力梯度仪线运动误差、自梯度实时补偿装置及技术的没有公开的报道。

技术实现要素：

技术问题：本发明提供一种能够补偿重力梯度仪线运动引起的测量误差的旋转加速度计重力梯度仪线运动误差补偿装置，既能抑制温度、电磁场等环境因素对重力梯度仪的线运动误差传递系数的影响，同时还能避免重力梯度仪加速度引起的重力梯度仪前端信号调理电路过电压饱和和过电压损坏。本发明同时提供了一种实现上述效果、解决上述问题的旋转加速度计重力梯度仪线运动误差补偿方法。

技术方案：本发明的旋转加速度计重力梯度仪线运动误差补偿装置，包括：

用于检测重力梯度仪加速度的线运动检测模块；

用于产生正交幅度调制载波的参考信号产生模块；

用于根据反馈的补偿后的重力梯度仪信号实时微调线运动误差传递系数的线运动误差传递系数处理模块；

用于根据重力梯度仪加速度信号、正交幅度调制载波、线运动误差传递系数产生线运动误差补偿信号的线运动误差补偿信号产生模块；

用于根据重力梯度仪的姿态产生自梯度补偿信号的自梯度补偿信号产生模块；

用于根据线运动误差补偿信号、自梯度补偿信号对含有线运动误差、自梯度的重力梯度仪信号进行补偿的补偿运算模块；

所述线运动检测模块的输出连接到线运动误差补偿信号产生模块的输入；所述参考信号产生模块的输出连接到线运动误差补偿信号产生模块、自梯度补偿信号产生模块的输入；所述线运动误差传递系数处理模块的输入连接到补偿运算模块的输出，所述线运动误差传递系数处理模块的输出连接到线运动误差补偿信号产生模块的输入；所述线运动误差补偿信号产生模块、自梯度补偿信号产生模块的输出连接到补偿运算模块的输入。

进一步的，本发明装置中，所述线运动检测模块包括加速度计和低通滤波器，所述加速度计安装在重力梯度仪测量坐标系的x轴，y轴，z轴，测量重力梯度仪测量坐标系的加速度ax，ay，az；所述低通滤波器则滤除加速度信号中的高频噪声。

进一步的，本发明装置中，所述参考信号产生模块包括重力梯度仪旋转圆盘轴编码器和信号发生器，所述重力梯度仪旋转圆盘轴编码器检测重力梯度仪圆盘旋转的相位角φt，所述信号发生器根据相位角φt，产生正交幅度调制载波波sinφt，sin2φt，cosφt，cos2φt。

进一步的，本发明装置中，所述线运动误差传递系数处理模块包括线运动误差传递系数输入模块和线运动误传递差系数调节模块，所述线运动误差传递系数输入模块，输入线运动误差传递系数初值：所述线运动误差传递差系数调节模块，根据反馈的补偿后的重力梯度仪信号产生调节量，微调线运动误传递差系数；线运动误差传递系数处理模块具有两种工作模式，调节模式和非条节模式，当工作在调节模式时，实时调节线运动误差传递系数；当工作在非调节模式时，线运动误差传递系保持不变。

进一步的，本发明装置中，所述线运动误差补偿信号产生模块具有两种工作模式：非补偿模式，补偿模式；非补偿模式下，产生的线运动补偿信号cl(t)为：

cl(t)＝0；

在补偿模式下，产生的线运动补偿信号cl(t)为：

式中sin2φt，cos2φt，sinφt，cosφt为t时刻输入线运动误差补偿信号产生模块的正交幅度调制载波；表示t时刻输入线运动误差补偿信号产生模块的线运动误差传递系数；ax(t)，ay(t)，az(t)表示t时刻输入线运动误差补偿信号产生模块的加速度信号。

进一步的，本发明装置中，所述自梯度补偿信号产生模块包括姿态输入模块、自梯度模型参数输入模块、正交幅值调制载波输入模块、自梯度信号产生模块；所述姿态输入模块实时获取重力梯度仪的姿态数据；所述自梯度模型参数输入模块用于设置自梯度模型参数初值；所述正交幅值调制载波输入模块用于输入自梯度信号的正交幅值调制载波；所述自梯度信号产生模块根据输入的自梯度模型参数、重力梯度仪姿态、正交幅度调制信号产生自梯度补偿信号；自梯度补偿信号产生模块具有两种工作模式：补偿模式和非补偿模式；补偿模式下，自梯度补偿信号产生模块的输出csg(t)为：

非补偿模式下的输出为：

csg(t)＝0

式中，是重力梯度仪在t时刻的姿态，p是重力梯度仪自梯度模型参数；分别是重力梯度仪的inline通道、cross通道t时刻的自梯度；sin2φt，cos2φt是t时刻的正交幅值调制载波。

进一步的，本发明装置中，所述补偿运算模块根据自梯度补偿信号、线运动误差补偿信号对重力梯度仪的输出进行补偿。

本发明的旋转加速度计重力梯度仪线运动误差补偿方法，包括以下步骤：

1)根据工作模式，计算t时刻的线运动误差传递系数：

调节模式：非调节模式：

式中表示t-1时刻的线运动误差传递系数；gc(t-1)表示t-1时刻补偿后的重力梯度仪信号；f(gc(t-1))是t时刻线运动误差传递系数的微调量，它是gc(t-1)的函数；

2)检测旋转加速度计重力梯度仪旋转圆盘在t时刻的相位角φt，计算t时刻的正交幅度调制载波sinφt，sin2φt，cosφt，cos2φt；检测旋转旋转加速度计重力梯度仪在t时刻的加速度ax(t)，ay(t)，az(t)；

3)根据下式计算t时刻的3类线运动误差补偿信号cl1(t)，cl2(t)，cl3(t)：

4)根据工作模式，计算t时刻的总的线运动误差补偿信号cl(t)：

在非补偿模式时，cl(t)＝0；

在补偿模式时，cl(t)＝cl1(t)+cl2(t)+cl3(t)；

5)根据工作模式，计算t时刻的自梯度补偿信号csg(t)：

补偿模式时，

在非补偿模式时，非补偿模式下的输出为：

csg(t)＝0

6)对t时刻含有线运动误差、自梯度的重力梯度仪信号g(t)，根据下式进行线运动误差、自梯度补偿；

gc(t)＝g(t)-cl(t)-csg(t)

式中gc(t)是t时刻补偿后的重力梯度仪信号，g(t)是t时刻含有线运动误差、自梯度的重力梯度仪信号，cl(t)是t时刻的线运动误差补偿信号，csg(t)是自梯度补偿信号；

更进一步的，本发明方法中，所述步骤1)中，t＝0时刻的线运动误差传递系数均通过标定得到。

在动基座重力梯度勘探时，由于重力梯度仪内部的加速度计存在安装误差、加速度计一阶、高阶标度系数不匹配、电路增益不匹配等，导致重力梯度仪的加速度，角速度，角加速度传递到重力梯度仪的输出，造成测量误差，此外重力梯度仪的线运动误差传递系数、角运动误差传递系数容易受温度、电磁场等环境因素的影响。同时重力梯度仪的加速度，角速度，角加速度会引起重力梯度仪前端信号调理电路过电压饱和，或过电压损坏，本发明能够解决以上问题。

有益效果：本发明与现有技术相比，具有以下优点：

本发明提供一种能够对重力梯度仪线运动误差、自梯度进行实时补偿的装置和方法，它根据反馈的补偿后的重力梯度仪信号实时微调线运动误差传递系数，能够抵消重力梯度仪环境因素比如温度、磁场等对重力梯度仪线运动误差传递系数的影响，能够补偿重力梯度仪线运动引起的测量误差，提高重力梯度仪的测量精度。由于重力梯度仪线运动误差信号，自梯度信号均进行了正交幅值调制，直接对重力梯度仪解调前的信号进行补偿，因此还能够解决由于重力梯度仪加速度、自梯度引起的前端信号调理电路过电压饱和，过电压损坏问题。

附图说明

图1旋转加速度计重力梯度仪线运动误差、自梯度实时补偿装置结构图。

图2旋转加速度计重力梯度仪线运动传感器安装示意图

具体实施方式

下面结合实施例和说明书附图对本发明作进一步地说明。

由图1可知，本发明一种旋转加速度计重力梯度仪线运动误差补偿装置，它包括线运动检测模块，参考信号产生模块，线运动误差传递系数处理模块，线运动误差补偿信号产生模块，自梯度补偿信号产生模块，补偿运算模块，重力梯度仪加速度计信号处理模块，重力梯度信号恢复模块。

线运动检测模块，为线运动误差补偿信号产生模块提供重力梯度仪加速度信号；它由加速度计，低通滤波器组成；如图2所示，加速度计安装在重力梯度仪测量坐标系的x轴，y轴，z轴，测量重力梯度仪测量坐标系的加速度ax，ay，az；低通滤波器则滤除加速度信号中的高频噪声。

参考信号产生模块，用于产生正交幅度调制载波为线运动误差补偿信号产生模块提供正交幅度调制信号；参考信号产生模块，由重力梯度仪旋转圆盘轴编码器及信号发生器组成；重力梯度仪旋转圆盘轴编码器检测重力梯度仪圆盘旋转相位角φt，信号发生器根据相位角φt，产生正交幅度调制载波波sinφt，sin2φt，cosφt，cos2φt。

线运动误差传递系数处理模块，为线运动误差补偿信号产生模块提供线运动误差系数；线运动误差传递系数处理模块，由线运动误差传递系数输入模块，线运动误传递差系数调节模块组成；线运动误差传递系数输入模块，输入线运动误差传递系数初值：线运动误差传递差系数调节模块，根据反馈的补偿后的重力梯度仪信号gc(t)及工作模式，微调线运动误传递差系数线运动误差传递系数处理模块的两种工作模式为：调节模式和非调节模式；这两种工作模式下的线运动误差传递系数的调节如下：

调节模式：非调节模式：

线运动误差补偿信号产生模块，根据重力梯度仪线运动信号ax(t)，ay(t)，az(t)，正交幅度调制载波sin2φt，cos2φt，sinφt，cosφt，线运动误差传递系数产生线运动误差补偿信号cl(t)；线运动误差补偿信号产生模块具有2种工作模式，非补偿模式和补偿模式；非补偿模式下，产生的线运动补偿信号cl(t)为：

cl(t)＝0；

在补偿模式下，产生的线运动补偿信号cl(t)为：

自梯度补偿信号产生模块，具有两种工作模式，补偿模式和非补偿模式，在补偿模式下，产生的自梯度补偿信号csg(t)为：

在非补偿模式时，产生的自梯度补偿信号csg(t)为：csg(t)＝0。

重力梯度仪加速度计信号处理模块，对重力梯度仪旋转圆盘上安装的四只加速度计输出信号进行处理，得到含有线运动误差、自梯度的重力梯度仪信号g(t)。

补偿运算模块，根据线运动误差补偿信号cl(t)，自梯度补偿信号csg(t)，对含有线运动误差的重力梯度仪信号g(t)，根据下式进行线运动误差补偿：

gc(t)＝g(t)-cl(t)-csg(t)

补偿后的重力梯度信号gc(t)被反馈到线运动误差传递系数处理模块，用于实时调节重力梯度仪线运动误差传递系数，同时补偿后的重力梯度信号gc(t)输入重力梯度信号恢复模块，解调输出重力梯度信号。