一种用于修正原子干涉重力仪干涉条纹的振动补偿方法与流程

1.本发明涉及重力测量领域，尤其涉及一种用于修正原子干涉重力仪干涉条纹的振动补偿方法。


背景技术：

2.绝对重力测量，即对重力加速度绝对值的精确测量，在辅助导航、资源勘探、地球物理、计量领域都有着广泛的应用。现有的高精度绝对重力仪普遍采用自由落体法，主要包括以角锥棱镜为下落物体的激光干涉式绝对重力仪和以冷原子团为下落物体的原子干涉式绝对重力仪，以下分别简称为光学重力仪和原子重力仪。前者通过激光干涉仪测量角锥棱镜在真空中自由下落的运动轨迹，拟合求解重力加速度，后者根据原子团在自由下落的过程中经过三次拉曼激光脉冲作用形成的原子干涉条纹来计算重力加速度。由于二者的实际测量对象是下落物体或下落原子团相对于仪器内一个参考镜(角锥棱镜或平面镜)的运动加速度，因此参考镜自身的加速度必然耦合到测量结果中。理论上，忽略参考镜自身振动时，现有的绝对重力仪的测量精度可以达到微伽量级；但实际上，如果直接将参考镜放置在地面上，地面振动引起的测值离散度将超过十微伽，在复杂振动环境中甚至达到毫伽或十毫伽量级。因此，利用测振仪器获得参考镜自身的运动加速度并在计算重力加速度时对其影响进行修正是提高绝对重力仪测量精度的有效手段之一，这种修正方法一般称为振动补偿。
3.具体而言，影响原子重力仪的地面振动噪声主要为环境噪声，与测量时间和测量点的地基、坐标及环境有关，可以根据频率来粗略区分。首先是人类活动噪声，一般大于1hz，如仪器测试期间附近的人员走动，车辆等人工振源产生的振动，以及建筑物的晃动等。其次是来源于地球本身的噪声，一般在0.1hz～1hz之间，主要包括周期性地脉动和地震等地球内部运动引起的振动，其中地脉动噪声模型的加速度功率谱密度在频率为0.2hz和3hz的位置存在峰值。最后是气压波动等大气运动噪声，一般小于0.1hz，普遍变化缓慢且幅值较低，对绝对重力测量的影响可以忽略。由于绝对重力测量只依靠垂直方向的位移数据进行计算，因此对原子干涉条纹进行振动补偿时，一般忽略地面振动的水平分量，仅需使用测振仪器输出的垂直振动信号。
4.目前国际上较为成熟的原子重力仪是法国巴黎天文台研制的cag-01型原子重力仪，曾参加过多次国际重力比对，采用的振动处理方法为被动式垂直隔振系统与振动补偿的结合，其中振动补偿的具体原理为：数据采集卡同步记录原子干涉过程中的地震计输出的速度信号与原子干涉条纹，在控制软件中利用iir(infinite impulse response，无限冲激响应)滤波器和非因果低通滤波器对地震计信号进行处理，以减小地震计传递函数幅频特性不恒定、相频特性存在非线性的影响，获得更真实的参考镜振动信号，最后利用处理后的信号得到由参考镜振动引入的相位噪声，修正原子干涉条纹。在该仪器所在的测量环境中，此振动补偿方法本身可以使重力仪灵敏度提升为原来的3倍。德国汉诺威大学提出利用数字式低通和高通滤波器对振动传感器的输出信号进行前处理，以得到更为精确的参考镜
振动引入的相位噪声，从而实现振动补偿。其实验结果表明，在较安静的环境中利用高精度商用地震计进行振动补偿，或在复杂振动环境中使用精度较低的商用加速度计或其自主研制的新型光学惯性传感器进行补偿，均可以有效提高重力仪的灵敏度。
5.在我国，浙江大学自制的原子重力仪同样采用了振动补偿方法，其主要特点为利用商用反演软件修正地震计实测传递函数的影响以获得更真实的参考镜振动信号(如图1所示)。在其实验室环境中，该方法可以将重力仪的灵敏度提升至与使用被动隔振系统时相近的水平。国防科技大学对参考镜垂向振动及水平偏转引入的测量噪声进行了详细分析，利用加速度计测量原子重力仪工作时的参考镜振动，同时对加速度计与参考镜之间的传递函数进行了较为精确的标定。实测结果表明在较复杂的振动环境中其振动补偿方法可以实现显著的振动补偿效果。军事科学院国防创新院提出的振动补偿方法则采用具有高精度四通道移相探测器的迈克尔逊激光干涉仪替代传统振动传感器。该干涉仪可以直接测量出原子干涉条纹的相移，并结合信号整形、数字鉴相等方法确定拉曼激光脉冲作用前后原子干涉条纹的相位小数的变化，得到参考镜振动引起的原子干涉相位偏差。去除该相位偏差，即可对重力仪测量到的跃迁几率信号进行修正。此外，中国计量科学研究院、中科院精密测量研究院、华中科技大学、中国科技大学、国防科技大学、美国斯坦福大学、德国洪堡大学等单位也研制了不同类型的隔振系统，以减小地面振动对原子重力仪的影响。
6.然而，现有的原子重力仪适应性较差，在测量环境改变时，其振动补偿效果明显降低，从而影响了原子重力仪的测量精度，因此有必要对其进行改进。


技术实现要素：

7.本发明目的是针对上述问题，提供一种在各个测量环境下均具有最佳补偿效果的用于修正原子干涉重力仪干涉条纹的振动补偿方法。
8.为了实现上述目的，本发明的技术方案是：
9.一种用于修正原子干涉重力仪干涉条纹的振动补偿方法，包括以下步骤：
10.s1、给定扫频速率序列、重力参考值、脉冲间隔时长、有效波矢；
11.s2、导入实测数据；
12.s3、通过原始条纹计算并设定优化目标；
13.s4、根据触发信号分割振动传感器信号，得到每个跃迁几率对应的振动传感器信号；
14.s5、给定延时系数以及增益系数的起始范围；
15.s6、通过遍历搜索法或黄金分割搜索法求取最优延时系数、最优增益系数；
16.s7、通过最优延时系数、最优增益系数计算每次原子干涉时的参考镜振动速度，进而得到振动相位噪声序列、仅受振动影响的原子相位序列以及修正条纹，最终计算得到将修正条纹拟合到标准余弦信号时的rmse值。
17.进一步的，所述步骤s2中，实测数据包括与扫频速率序列对应的实测跃迁几率序列、同时段内测量参考镜振动的振动传感器垂直方向输出信号、同时段内每次原子团下落前的触发信号。
18.进一步的，所述步骤s3中，优化目标为原始条纹拟合时的rmse值、原始条纹拟合时的可决系数、干涉条纹残差的标准差、原始条纹残差与仅受振动影响的条纹残差之间的相
关系数中的一种。
19.进一步的，所述步骤s4中，用校正器传递函数对振动传感器信号进行滤波处理操作。
20.进一步的，所述步骤s6中，使用遍历搜索法或黄金分割搜索法求取最优延时系数、最优增益系数时，当优化目标为原始条纹拟合时的rmse值时，以原始条纹拟合时的rmse值达到最小值时对应的延时系数、增益系数作为最优延时系数、最优增益系数；
21.当优化目标为原始条纹拟合时的可决系数时，以原始条纹拟合时的可决系数达到最小值时对应的延时系数、增益系数作为最优延时系数、最优增益系数；
22.当优化目标为干涉条纹残差的标准差时，以干涉条纹残差的标准差达到最小值时对应的延时系数、增益系数作为最优延时系数、最优增益系数；
23.当优化目标为原始条纹残差与仅受振动影响的条纹残差之间的相关系数时，以原始条纹残差与仅受振动影响的条纹残差之间的相关系数的绝对值达到最大值时对应的延时系数作为最优延时系数；以其他优化目标达到最小值时对应的增益系数作为最优增益系数。
24.进一步的，所述步骤s7中，通过最优延时系数、最优增益系数计算每次原子干涉时的参考镜振动速度时，当振动传感器为地震计时，则计算公式为：
[0025][0026]
式中，vm(t)为参考镜的振动速度；τ为最优延时系数；k为最优增益系数；us(t)为地震计输出信号；vs为地震计敏感质量振动速度；ks为地震计的标称灵敏度；
[0027]
当振动传感器为加速度计时，则计算公式为：
[0028][0029]
式中，vm(t)为参考镜的振动速度；τ为最优延时系数；k为最优增益系数；ka为加速度计的标称灵敏度；ua(t)为加速度计输出信号。
[0030]
进一步的，所述步骤s7中，通过最优延时系数、最优增益系数计算每次原子干涉时的参考镜振动速度时，当振动传感器为加速度计时，则采用半周期积分算法计算参考镜振动速度。
[0031]
与现有技术相比，本发明具有的优点和积极效果是：
[0032]
本发明对环境具有极好的自适应性，其可以在测量环境改变时实现当时、当地最好的振动补偿效果，解决了现有振动补偿方法在测量环境改变时振动补偿效果相对下降的问题，有效提高了原子重力仪在振动环境中的振动补偿效果，从而提高了原子重力仪的测量精度，进一步给绝对重力的测量工作作出了一定的贡献。
附图说明
[0033]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可
以根据这些附图获得其他的附图。
[0034]
图1为浙江工业大学的振动补偿方法流程图；
[0035]
图2为原子干涉条纹示意图；
[0036]
图3为地震计输出信号与参考镜真实振动的框架关系图；
[0037]
图4为实施例1的振动补偿算法流程图；
[0038]
图5为实施例2的振动补偿算法流程图；
[0039]
图6为实施例10的振动补偿算法流程图。
具体实施方式
[0040]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
[0041]
一、对原子干涉重力仪干涉条纹进行振动补偿的通用原理
[0042]
原子重力仪基于冷原子物质波干涉原理，将激光脉冲作用到冷原子团，通过双光子受激拉曼跃迁或多光子布拉格衍射过程来实现原子干涉过程。以基于受激拉曼跃迁的原子重力仪为例，冷原子团在真空腔中自由下落，真空腔下方的平面参考镜将垂直向下的拉曼激光原路反射，使两束对射的激光同时作用于原子团，在此过程中三个连续的激光脉冲实现了原子波包的分束、反转和合束，分别为第1个π/2脉冲、π脉冲和第2个π/2脉冲，相邻脉冲的时间间隔为t。合束后任一能态的原子均由两条路径上的原子叠加而成，通过荧光信号可以探测脉冲作用后原子团的跃迁几率p，即干涉条纹，满足
[0043][0044]
其中，δф为原子相位，与当地的重力加速度g有关；代表条纹对比度的系数a和b可以通过对干涉条纹进行余弦拟合得到。理想情况下，原子相位δф等于其理论值δф
th
，即
[0045]
δφ
th
＝(k
eff
g-α)t2.
ꢀꢀꢀꢀꢀꢀ
(2)
[0046]
其中，k
eff
为有效波矢，α表示激光的扫频速率。重力仪的控制软件调控式(2)中的α线性增大，则原子相位δф也将线性变化，使干涉条纹呈现为标准余弦信号的形式。此时p应与由p拟合出的标准信号p
fit
完全相同，有
[0047][0048]
但实际情况下，受到各种相位噪声的影响，原子相位δф等于其实际值δфm，有
[0049][0050]
其中，δφ
vib
为参考镜振动引入的相位噪声，以下简称为振动相位；δφothers为其他噪声源引入的相位噪声的总和，以下简称为其他相位噪声。此时p由p
fit
变为
[0051][0052]
当振动相位在相位噪声中占主导时，依据上述原理可以确定修正干涉条纹p的振动补偿过程。
[0053]
第一步，对原始实测条纹p(δф
th
)进行余弦拟合(注：条纹即为曲线，括号外为曲线的纵坐标，括号内为曲线的横坐标)，得到拟合出的理论条纹p
fit
(δф
th
)和式(1)中的系数a和b。
[0054]
第二步，通过振动传感器获得一段时间t内参考镜的参考镜的真实振动速度vm，进而根据
[0055][0056]
计算出真实的振动相位δφ
vib
，其中s(t)为重力仪的灵敏度函数，t1和t3分别为第1个π/2脉冲和第2个π/2脉冲的作用时刻，均由重力仪控制软件的已知参数决定。还需要说明的是，虽然重力值g是待测的未知量，但在实际测量的初始阶段可以利用近似的粗测值g
th
代入式(2)计算理论原子相位δф
th
以获得拟合条纹p
fit
，再通过后续测量和信号处理得到精确的实测值g。
[0057]
第三步，计算出仅受振动影响时的原子相位δфv和跃迁几率pv，即
[0058][0059]
此时重新以δфv为横坐标、p为纵坐标得到的数据对p(δфv)即为修正后的干涉条纹，至此补偿过程结束。如图2所示，与原始条纹p(δф
th
)相比，修正条纹p(δфv)的横坐标有所移动，理论上应更接近标准余弦信号，即拟合出的理论条纹p
fit
(δф
th
)。也就是说，对修正后的条纹p(δфv)再次进行余弦拟合时的拟合优度应高于对原始条纹p(δф
th
)进行拟合时的拟合优度。
[0060]
二、本发明描述的振动补偿方法的具体原理
[0061]
基于传递函数简化模型以实现上述振动补偿过程的方法包括硬件和软件两方面。
[0062]
硬件方面，采用一台独立的地震计或加速度计放在参考镜旁边或下方测量其振动。以地震计为例，受机械结构及地震计自身性能的限制，地震计输出信号us(t)并不能完全真实地反应参考镜的振动速度vm(t)，二者的关系如图3所示，其中ga表示从地面振动速度到参考镜振动速度的传递函数，gb表示从地面振动速度到地震计敏感质量振动速度的传递函数，gc表示从敏感质量运动速度到地震计输出电压之间的传递函数。由于地震计内部的敏感质量与参考镜的实际位置必然存在水平和垂直差异，因此ga与gb不可能相等，且当地基材料改变时，这两个传递函数本身也可能发生变化。
[0063]
另一方面，地震计自身的传递函数gc并非幅频特性恒定、相频特性为线性的理想传递函数，且有可能随测量时间与环境的改变而发生变化(这是因为温度、湿度、气压等环境参数会对地震计内部的机电模块参数产生影响，且这些参数可能随时间发生漂移)。综上，参考镜的真实振动速度与地震计输出的电压信号之间的传递函数h满足
[0064][0065]
虽然实际情况中ga、gb、gc不易测量且可能改变，导致测量人员无法得到总传递函数h的精确值，但由于单次重力测量任务中测量环境不会有显著变化，且测量时间一般不超过24小时，因此可以利用由一个延时系数τ和一个增益系数k组成的简化模型来估算h，并认为这两个系数的数值在当前测量过程中保持不变。由此可得图3的简化模型为
[0066][0067]
其中，ks为地震计的标称灵敏度，由制造商提供。
[0068]
理论上当延时系数τ和增益系数k取最接近真实值时，从地震计输出电压us，根据式(9)计算出的参考镜振动速度vm也最接近真实值，对条纹的修正效果最好。因此振动补偿的具体思路为：给定τ和k的取值范围，选取某一表征拟合优度的参数值为优化目标，查找使该参数达到最优值时对应的最优延时系数τ
opt
和最优增益系数k
opt
。查找完毕后，设定延时为τ
opt
、增益为k
opt
的情况下，得到修正后的条纹p(δфv)，重新计算作为优化目标的参数值。
[0069]
由于实际上地震计自身的传递函数gc并非理想形式，即幅频特性非常数、相频特性非线性等特点，因此地震计的输出信号us(t)可能存在失真。所以在搜索延时系数τ和增益系数k前，可以首先用一个校正滤波器f对信号us(t)进行处理，减小其失真程度。
[0070]
硬件部分由地震计更换为加速度计时，可以直接对加速度计的输出信号ua(t)进行一次积分，结合延时系数和增益系数，得到参考镜振动信号vm(t)。
[0071]
相应的计算公式为
[0072][0073]
其中，ka为加速度计的标称灵敏度，由制造商提供；积分时间根据实际需要设定，至少包含原子团自由下落的整个时间段。也可以采用半周期积分等其他积分算法从加速度计信号ua(t)计算参考镜速度信号vm(t)。
[0074]
三、实施例
[0075]
(1)实施例1
[0076]
一种用于修正原子干涉重力仪干涉条纹的振动补偿方法，其特征在于：包括以下步骤：
[0077]
s1、给定扫频速率序列、重力参考值、脉冲间隔时长、有效波矢；
[0078]
s2、导入与扫频速率序列对应的实测跃迁几率序列、同时段内测量参考镜振动的振动传感器垂直方向输出信号、同时段内每次原子团下落前的触发信号；
[0079]
s3、通过原始条纹计算并设定原始条纹拟合时的rmse值为优化目标；
[0080]
s4、根据触发信号分割地震计信号，得到每个跃迁几率对应的地震计信号；
[0081]
s5、给定延时系数以及增益系数的起始范围；
[0082]
s6、通过遍历搜索法求取最优延时系数、最优增益系数；
[0083]
s7、通过最优延时系数、最优增益系数计算每次原子干涉时的参考镜振动速度，进而得到振动相位噪声序列、仅受振动影响的原子相位序列以及修正条纹，最终计算得到将修正条纹拟合到标准余弦信号时的rmse值。
[0084]
该实施例中，硬件传感器采用地震计，输出信号为us(t)。以拟合条纹时的rmse值(root mean square error,均方根误差)为优化目标，先遍历搜索延时系数τ，再遍历搜索增益系数k。取rmse值达到最小值时分别对应的延时系数和增益系数为τ
opt
和k
opt
。具体的算法流程如图4所示。
[0085]
(2)实施例2
[0086]
其步骤和实施例1相同，仍以拟合条纹时的rmse值为优化目标，实施例2采用黄金分割法同时搜索延时系数τ和增益系数k。具体的算法流程如图5所示。
[0087]
(3)实施例3
[0088]
对图4中“数据预处理”后的地震计信号us(t)增加一个校正步骤，即用设定好的校正器传递函数f对该信号进行滤波，具体的滤波函数可以自主编写或利用软件自身的函数库。再对校正后的信号执行实施例1中“遍历搜索最优参数”和“结果计算”的流程。
[0089]
(4)实施例4
[0090]
以拟合条纹时的r-square(可决系数)为优化目标，取r-square值达到最小值时分别对应的延时系数和增益系数为τ
opt
和k
opt
。搜索延时系数搜索和增益系数k的方法与实施例1相同，即遍历搜索。
[0091]
(5)实施例5
[0092]
以拟合条纹时的r-square(可决系数)为优化目标，取r-square值达到最小值时分别对应的延时系数和增益系数为τ
opt
和k
opt
。搜索延时系数和搜索增益系数k与实施例2相同，即采用黄金分割法。
[0093]
(6)实施例6
[0094]
rmse值具有简单直观的数学含义，但物理意义不明显。因此可以用另一种更具有物理意义的变量作为优化目标，即干涉条纹拟合残差的变化。搜索系数的步骤与实施例1相同，仍为先遍历搜索延时系数τ，再遍历搜索增益系数k。
[0095]
运行振动补偿算法前，先计算出原始条纹p(δф
th
)相较于理论条纹p
fit
(δф
th
)的残差rm(δф
th
)的标准差σm。运行振动补偿算法时，给定当前的延时系数τ，根据式(6)(7)(9)得到仅受振动影响的干涉条纹pv相较于理论条纹p
fit
的残差rv。同时，记原始条纹p与仅受振动影响的条纹pv的差值为rc。显然rc也就是原始条纹残差rm与仅受振动影响的条纹残差rv的差值，相当于去除振动影响后的条纹残差。理论上振动补偿效果越好，则振动相位δφ
vib
的计算越准确，推算出的仅由振动引入的条纹残差rv也越接近原始条纹残差rm。此时rc将趋近于仅由其他相位噪声δφ
others
引入的条纹残差，其标准差σc也应趋于最小值。因此，实施例6将σc设为优化目标，以σc达到最小值时分别对应的延时系数和增益系数为τ
opt
和k
opt
；延时系数和增益系数的搜索方法同实施例1，即分别遍历搜索。
[0096]
(7)实施例7
[0097]
以实施例6中的干涉条纹残差标准差σc为优化目标，采用实施例2中的黄金分割法，同时搜索延时系数τ和增益系数k。
[0098]
(8)实施例8
[0099]
如实施例6所述，理论上振动补偿效果越好，推算出的仅由振动引入的条纹残差rv也越接近原始条纹残差rm。因此实施例8以rm和rv的相关系数cr为优化目标，遍历搜索延时系数τ，以cr的绝对值达到最大值时对应的延时系数为τ
opt
。
[0100]
由于计算相关系数时需要对原始数据进行归一化处理，因此cr无法作为搜索增益系数时的优化目标。此处可以采用实施例1中的rmse值作为遍历搜索增益系数k的优化目标。
[0101]
(9)实施例9
[0102]
以实施例8中的相关系数cr为优化目标，遍历搜索延时系数τ。再以实施例6中的干
涉条纹残差标准差σc为优化目标，遍历搜索增益系数k。
[0103]
(10)实施例10
[0104]
振动补偿的硬件传感器更换为加速度计，输出信号为ua(t)。此时采用一次积分并结合延时系数和增益系数即可从ua(t)得到参考镜的振动速度信号vm(t)。仍与实施例1步骤相同，以拟合条纹时的rmse值为优化目标，分别遍历搜索延时系数τ和增益系数k。具体的算法流程如图6所示。
[0105]
(11)实施例11
[0106]
振动补偿的硬件传感器仍为加速度计，实施例10中从加速度信号ua(t)计算速度信号vm(t)的方法采用半周期积分法。仍以拟合条纹时的rmse值为优化目标，分别遍历搜索延时系数τ和增益系数k。
[0107]
本发明对环境具有极好的自适应性，其可以在测量环境改变时实现当时、当地最好的振动补偿效果，解决了现有振动补偿方法在测量环境改变时振动补偿效果相对下降的问题，有效提高了原子重力仪在振动环境中的振动补偿效果，从而提高了原子重力仪的测量精度，进一步给绝对重力的测量工作作出了一定的贡献。