航空重力梯度仪实时重力梯度解调方法

钱学武 ，赵立业，尹航

(1.东南大学 仪器科学与工程学院，江苏 南京 210096； 2.济南职业学院 电子工程学院，济南 250103； 3.中国自然资源航空物探遥感中心，北京 100083)

0 引言

重力梯度是重力加速度的空间微分，能够反映场源体更多的细节，具有比重力更高的分辨率。重力梯度测量对于空间科学、地球科学、地质科学的发展以及资源勘探、惯性导航等研究具有非常重要的意义[1-4]。20世纪末，澳大利亚BHP公司和美国Lockheed Martin公司合作共同研制了基于加速度计的部分张量航空重力梯度测量系统FALCON，随后Lockheed Martin公司研制了全张量重力梯度测量系统Air-FTGTM，英国ARKeX公司采用Lockheed Martin公司重力梯度专利技术开发了全张量重力梯度测量系统FTGeX，这3种重力梯度测量系统经过不断测试和改进，证实了其在能源勘探方面的巨大优势，并取得了很好效果[5]。目前，国际上正在研发的重力梯度仪主要包括旋转加速度计重力梯度仪、超导重力梯度仪、冷原子重力梯度仪、静电重力梯度仪，还有最近兴起的基于微机械加工(MEMS)技术的重力梯度仪等，其中旋转加速度计重力梯度仪是目前唯一成功用于机载/船载动基座上的重力梯度测量仪器[6]。

1 旋转加速度计重力梯度仪工作原理

将4只高精度加速度计(A)对称安装在慢速旋转圆盘上，其敏感轴方向沿着圆盘切线方向且相互正交，相对两只加速度计敏感轴方向相反，并且与圆盘旋转轴垂直，各加速度计质量中心到圆盘中心的距离相等，圆盘以固定角频率旋转，相对两只加速度计输出信号相加再相减，可以降低圆盘平动加速度和旋转角加速度对重力梯度信号的影响。对重力梯度仪输出信号进行滤波和解调，就可以得到重力梯度分量，为了实时获取圆盘旋转角度用于重力梯度解算，在圆盘中心位置处放置高精密光栅角编码器。旋转加速度计重力梯度仪测量原理如图1所示。

图1 旋转加速度计重力梯度仪测量原理

2 重力梯度解调相位角分析

重力梯度仪样机是由复杂机械结构、电气线路和控制算法等构成的高精密测量仪器。由于重力梯度仪所处的环境非常恶劣，如果直接对输出信号进行重力梯度解调，必然会引入额外的干扰，导致重力梯度解调过程中造成数据的不准确，因此，在进行重力梯度信号解调前，需要实时确定出精确的解调相位角。重力梯度解调相位角结构如图2所示。

图2 重力梯度解调相位角结构

4只加速度计输出信号经过信号调理处理，得到模拟信号Eout(t)，角度信息转换器输出同步脉冲信号给A/D转换器，A/D转换器对模拟信号进行A/D转换得到第(i-1)Δt时刻数字信号Eout(i)，其表达式为[9]：

(1)

参考信号源数字正余弦信号，它由DSP发出[9]，其表达式为：

，

(2)

式中：ϑj、ζj分别为正弦、余弦信号参考源的初始相角在第j次的修正角。

利用参考信号源数字正余弦信号对A/D转换器的输出信号进行2倍频幅值解调，通过改变其初始相角直到得到最大的幅值为止，最终的初始相角即为重力梯度解调相位角[9]，其计算表达式为：

，

(3)

式中：N为圆盘每旋转一圈角度信息转换器输出的脉冲数，ϑq、ζq分别为第q次得到的重力梯度解调相位角；对得到的重力梯度解调相位角ϑq、ζq进行组合，然后把组合的数据通过零相位低通滤波器进行滤波处理，最后对处理后的相位角进行均值处理，最终得到重力梯度解调相位角α。

3 重力梯度实时解调分析

为了实现航空重力梯度实时测量，在求解出解调相位角的基础上，添加重力梯度解调部分[10]，结构图如图3所示。

图3 航空重力梯度仪实时重力梯度解调结构

将4只加速度计输出信号经过信号调理和滤波电路进行处理，得到模拟信号Eout(t)，角度信息转换器输出同步脉冲信号给A/D转换器，A/D转换器对模拟信号Eout(t)进行A/D转换得到第(k-1)Δt时刻数字信号Eout(t)[10]，其表达式为

Eout(k)=2KIKR{(Γyy-Γxx)sin2[ωΔt(k-1)-α]+2αxycos2[ωΔt(k-1)-α]}

，

(4)

式中：k为数据序列索引；K为系统信号放大增益; Δt为采样时间;KI为加速度计标度因数；R为圆盘半径；ω为圆盘旋转角频率；α为重力梯度解调相位角； Гxx为x轴上的重力加速度分量在x轴方向上的空间导数； Гyy为y轴上的重力加速度分量在y轴方向上的空间导数； (Гyy- Гxx)为重力梯度Гyy与Гxx之差； Гxy为x(y)轴上的重力加速度分量在y(x)轴方向上的空间导数。

圆盘光栅编码器上的角度信息转换器实时输出圆盘角位置信息为θp,k，利用角位置信息θp,k对重力梯度数字信号进行实时重力梯度解调，得到具有时空信息的重力梯度数据，重力梯度解调表达式为：

，

(5)

式中：q为当前时刻的圆盘旋转圈数序列索引；M为圆盘旋转总圈数；gps为GPS提供的位置信息。

4 仿真验证

根据旋转加速度计重力梯度仪工作原理，采用高性能计算机、可编程高精度电流源、低噪声电流放大器、低噪声电压放大器、多路切换开关和高精度数字电压表等仪器组建了重力梯度半实物仿真分析系统。计算机把仿真数据发送给仿真分析系统，转换为接近真实的重力梯度仪工作环境下的重力梯度信号，通过改变各种实验参数和实验数据，可以全方位模拟重力梯度信号，有助于进一步对重力梯度仪的信号特征、误差分析、加速度计标度因子反馈调整算法以及信号处理方法等进行详细研究[11]，重力梯度半实物仿真分析系统实物图如图4所示。

图4 重力梯度仪半物理仿真系统实物

设定圆盘基线距离为1 m，加速度计标度因数为32 mA/g，重力加速度为9.81 m/s2，圆盘旋转周期为20 s，采样率为2 Hz，初始时刻，加速度计A1与x轴之间的夹角为0°，环境物体引起的方位角为0°，对重力梯度信号进行数据采集并处理，角度变动步长为0.01°，解调相位角结果如图5所示。从图5中可以看到，重力梯度解调相位角在-50.3°左右，零相位滤波处理前的解调相位角波动较大，标准差为0.14°，通过零相位滤波处理后降低了波动范围，标准差变为0.046°，最后做均值处理，得到重力梯度解调相位角均值为-50.29°。如果采用0°初始相位角作为重力梯度解调相位角，则会导致较大的重力梯度测量误差，从而影响重力梯度测量精度[10]。

图5 重力梯度解调相位角在零相位滤波前后比较波形

通过计算得到该环境物体引起的重力梯度分量(Гyy-Гxx)和Гxy的理论值分别为-335 Eu和0 Eu，如果不做重力梯度解调相位角处理，采用0°初始相位角作为重力梯度解调相位角，解调得到的重力梯度分量(Гyy-Гxx)和Гxy分别为61.89 Eu和164.62 Eu，与理论梯度值之间的误差达到了几百Eu，因此直接对原始重力梯度数据解调会产生非常大的测量误差，无法进行航空重力梯度测量；若采用计算得到的解调相位角作为重力梯度解调初始相位角对重力梯度信号进行解调，解调得到重力梯度分量(Гyy-Гxx)和Гxy分别为-335.008 5 Eu和0.000 55 Eu，与理论梯度值之间的误差小于0.01 Eu，能够满足航空重力梯度测量需求，重力梯度解调相位补偿前后的梯度信息比较如图6、图7所示。

图6 相位补偿前后的梯度信息(Гyy-Гxx)比较的波形

图7 相位补偿前后的梯度信息Гxy比较的波形

5 结论

本文针对航空重力梯度仪在测量过程中解调相位角对重力梯度测量精度的影响，提出了实时确定解调相位角和重力梯度解调求解的方法。首先提出了解调相位角的确定方法和求解算式；然后给出了实时重力梯度解调算式；接着利用解调相位角对重力梯度仪输出数据进行实时重力梯度解调，得到更高精度的重力梯度数据；最后为了验证所提方法的有效性，在重力梯度半物理仿真系统上进行半物理仿真实验测试，仿真结果表明，所给的初始相位角解算式和重力梯度实时解调算式能够明显提高重力梯度测量精度，可以为重力梯度仪的研制提供技术参考。