绝对重力仪的技术发展: 光学干涉和原子干涉 下载: 1906次特邀综述
1 引言
重力加速度随着时间和空间而不断变化,重力加速度测量按照测量结果分为绝对重力测量和相对重力测量,绝对重力测量通常为相对重力测量提供参考标准,是保证所有重力加速度测量结果具有溯源性和准确性的必要手段。在国际单位制(SI)中,重力加速度的单位是 m·s-2;在实际应用中,通常用“伽”(1 Gal=1×10-2 m·s-2)、“毫伽”(1 mGal=1×10-5 m·s-2)和“微伽”(1 μGal=1×10-8 m/s2)来表征重力加速度的测量结果。重力加速度测量在计量科学、资源勘探、海洋监测等领域有着广泛的应用,它既是航天器飞行等无源导航的主要方法,又是开展地球科学研究、揭示地球物理现象的关键手段。目前,激光干涉绝对重力仪和原子干涉绝对重力仪是开展绝对重力测量的主要手段;国际上激光干涉绝对重力仪的合成标准不确定度最优可达到1.8 μGal,原子干涉绝对重力仪的合成标准不确定度最优可达到4.5 μGal。与原子干涉绝对重力仪相比,激光干涉绝对重力仪起步时间早,发展更为成熟,以FG5/FG5X为代表的商用产品已成为相关行业主要使用的仪器,在历次绝对重力仪国际比对中也占据绝对主导地位。与激光干涉绝对重力仪相比,原子干涉绝对重力仪无机械磨损、测量效率和灵敏度高,具有非常好的发展前景。本文简要阐述了激光干涉方法和原子干涉方法开展绝对重力测量的原理,概括性介绍了激光干涉绝对重力仪和原子干涉绝对重力仪的研究现状和技术进展。最后,对绝对重力测量的未来发展方向进行了讨论与展望。
2 激光干涉绝对重力仪
2.1 基本原理
在地球表面及其附近,物体在真空环境中自由下落的位移随时间变化的关系可以表示为
式中:z0为第一个测量点的位置;v0为落体在此处的初始速度;g0为此处的重力加速度;gz为此处的重力梯度垂直分量。实际情况中,落体的下落距离一般小于0.5 m,该段内重力梯度垂直分量gz近似为定值,此时(1)式可以改写为
因此,在gz未知或已知但精度不高时,通过测量物体在真空中的自由下落轨迹z(t)并对其进行二次拟合,可以求解出有效高度处的重力加速度值,理论上表示为
这里的有效高度即与前面的z1相关联,指的是在下落起点的正下方、与起点距离为z1的高度。
目前典型的激光干涉绝对重力仪分别以激光波长和原子钟作为长度基准和时间基准,内部装有角锥棱镜的落体作为敏感元件[1],落体在真空中自由下落的运动加速度即为当地的重力加速度。角锥棱镜是一种可以保证反射光与入射光绝对平行、仅传播方向相反的特殊光学器件,以它为敏感元件可以自然消除水平扰动对垂直方向上重力测量的影响,大幅提高测量精度。绝对重力仪的主要光学测量结构为迈克耳孙干涉仪,如
2.2 技术实现
激光干涉绝对重力仪的结构如
2.2.1 仪器布局
早期的激光干涉绝对重力仪的结构布局与现在不太相同,其中最典型的是美国JILA实验室于20世纪80年代研制的JILA型绝对重力仪,如
之后,JILA实验室进一步研制出FG5型和FG5X型绝对重力仪,如
图 4. Micro-g Lacoste公司的FG5和FG5X型绝对重力仪[5-6]
Fig. 4. FG5 and FG5X absolute gravimeters from Micro-g Lacoste company[5-6]
2.2.2 落体运动方向
目前,绝大多数激光干涉绝对重力仪中,落体在释放前位于真空腔内顶部,自由下落后在底部被传动机构中的托盘承接,再被运输至顶部,从而实现仪器的反复测量,如上述FG5型绝对重力仪。但也有部分仪器中的落体在释放前始终位于真空腔内底部,被弹射器射出,具有先上抛后下落的运动轨迹。这种方案有两大优势,一是运动路径对称,可以减少残余空气的干扰;二是最少只需测量两个落体经过的位置点便可在一次下落后计算出重力加速度。因此,早期真空腔的真空度受限时,采用这种方案可以用比同类的自由下落式绝对重力仪更为紧凑的实验结构来实现测量。此外,上抛式绝对重力仪采用激光干涉来实现多点测量后,相比具有相同高度的自由落体型绝对重力仪,落体的飞行时间更长,对应更大的数据量。从统计角度来看,地面微振动造成的影响也会更小。不过此类仪器中落体在被弹射时初速度中可能具有较大的水平速度分量,会导致科里奥利力对测量结果产生影响,这可能是一项影响最大的误差源。上抛式的绝对重力仪中最典型且目前精度最高的是意大利计量院研制的IMGC-02型绝对重力仪,如
2.2.3 传动机构
激光干涉绝对重力仪普遍采用自由下落式落体释放方案,但不同型号的仪器中传动机构也存在差别,共有钢带、钢丝绳、凸轮、齿轮、齿条等多种形式。FG5型、FG5X型绝对重力仪以及中国计量科学研究院(NIM)NIM-3A型绝对重力仪采用的传动机构是钢带;清华大学T-1型绝对重力仪则采用钢丝绳进行机械传动;中国科学院测量与地球物理研究所IGG-02型绝对重力仪采用的传动机构是齿轮齿条。2005年,美国JILA实验室Faller等[9]在JILA型绝对重力仪的基础上研制出FGC型凸轮式绝对重力仪,如
2.2.4 振动处理方式
上述绝对重力仪普遍采用被动式或主动式垂直隔振系统来抑制地面振动对参考棱镜的影响。除此以外,还有一种抑制地面振动影响的方法是直接测量参考棱镜的振动信号并用该信号补偿测得的干涉信号以修正测量值,中国计量科学研究院研制的NIM系列绝对重力仪即属于此类。1975年,中国计量科学研究院研制出我国第一台固定的自由落体式绝对重力仪,测量准确度为100 μGal。1982年,几乎与国外同行同步,中国计量科学研究院完成了我国第一台可移动的绝对重力仪的研制,即NIM-1型绝对重力仪[10]。该仪器曾参与第一届绝对重力仪国际比对(ICAG),测量结果的合成标准不确定度为20 μGal,达到当时的国际先进水平。1985年,中国计量科学研究院成功研制了第二代的NIM-2型绝对重力仪,该仪器分别以碘稳频激光器和铷原子钟为长度基准和时间基准,测量落体轨迹中的大量位置点,从而计算出绝对重力值。与美国JILA实验室有所区别的是,该仪器首次采用地震计获得地面的微振动信号,并利用该信号对测得的落体运动轨迹进行补偿,在不使用隔振系统的情况下使测值的离散度得以大幅降低。NIM-2型绝对重力仪参与了第二次绝对重力仪国际比对[11],并在国内和区域内完成了大量验证性测量。在此基础上,中国计量科学研究院进一步研发完成了新一代NIM-3A型绝对重力仪,并于2014年被国家质检总局批准为重力加速度社会公用计量标准装置,如
除上述典型仪器外,法国乔治S.A.公司制造的GA60型上抛式绝对重力仪是第一台可移动的商用绝对重力仪[12];日本东京大学地震研究所、俄罗斯计量院研制的绝对重力仪使用压电致动器来释放下落物体[13-14];德国马普研究所的Rothleitner[15]也研制出了MPG-1和MPG-2型激光干涉式绝对重力仪。
2.3 重力测量不确定度评估
激光干涉绝对重力仪利用评估来分析测量偏差,对测量值进行修正,同时给出评估修正的不确定度。激光干涉绝对重力仪不确定度的B类评估一般由以下几个部分组成。
2.3.1 来源于仪器自身的误差
1) 时间与位移测量的准确程度
根据落体的运动方程即(2)式,可以得到时间、位移的测量引入的重力加速度的不确定度分量为
式中:
2) 空气阻尼
通过下式可以计算出真空腔内残余空气的阻尼引起的重力测量偏差[5,16],
式中:Fd为空气阻力;m为落体质量;A为落体表面积;ρ为平均气体密度;V ≈ 476 m/s为氮气分子在T = 300 K时的平均速度;v为信号采集过程中下落物体相对于残余空气的最大速度。平均气体密度ρ又可以在已知真空腔真空度(气压)P、玻尔兹曼常数kB、室温T、氮气分子质量mr的情况下用理想气体状态方程推导得到。该项的修正值即可设定为此偏差值的相反数,后续分析同理。该项引起的重力测量不确定度一般在0.1 μGal左右。
3)温度梯度
真空腔内的温度梯度导致残余空气形成气压梯度,从而引起重力测量偏差和不确定度。偏差值可以用下式估计[5]。
式中:Acros为落体在与垂直方向正交的平面内的截面积;T为真空腔内的温度。利用常用参数估计,落体实际行程两端的温度梯度一般小于0.1 K,引起的不确定度约为0.1 μGal。
4) 光束垂直度和光束发散
现有激光器发出的光束一般为高斯光束,具有一定的发散角,即使经过准直也不可能使发散角降为0。但使用激光干涉法测量落体轨迹时必须保证激光具有较小的发散角,且波矢方向沿垂直方向,否则会引入一般为0.5 μGal左右的测量误差[15]。假设光束偏角θ为矩形分布,此时可以按照下式来计算光束垂直度引起的修正值和不确定度[17],
式中:重力参考值g0取9.8 m·s-2即可。还可以按照下式来计算光束具有ϕ的发散角时对应引起的重力测量偏差和不确定度,
式中:系数K可以根据实际经验和所用激光器资料来选取,如10%。
5) 光速有限
虽然光速远大于落体的运动速度,但毕竟光速有限,因此干涉条纹实际形成的时刻将略微滞后于测量臂光束的波前与反射镜相遇的时刻。这两个时刻的时间差称为延迟时间,会导致重力测值略大于真值,称为有限光速效应,在精度达到微伽量级的绝对重力测量中必须予以考虑。直接利用光速和理论公式来修正落体轨迹中的测量点时刻,并且修正后的时间位移对数据进行二次拟合的修正方法比较复杂和繁琐,研究人员一般不采用,而是在单次结束后直接修正此项测量引入的偏差[18-19]。干涉条纹信号给出的落体轨迹测量实际上是一组等位移间隔的时间位移对,根据如下公式可以计算利用最小二乘法来拟合时的修正值和不确定度[20]。
式中:v0和σv0分别为落体的初速度及其不确定度;c为光速,由此得到的不确定度一般也在0.5 μGal左右。
6) 光电探测器的非线性效应
激光干涉仪中的光电探测器采集激光干涉条纹信号时会对其进行放大,但此过程也会引入相移。由于实际条纹信号为含有频率变化的信号,当该相移为常量或在信号的频带内具有线性相频曲线时,引入的重力测量偏差为0;但如果该相频曲线具有非线性,则一定会引入重力测量偏差[5]。目前常用的光电探测器造成的该项不确定度一般在1 μGal以内[5]。
7) 仪器本身的自吸引效应
绝对重力仪自身对落体的万有引力会引起重力测量不确定度,这种现象称为自吸引效应。通过有限元分析等方法可以对该效应进行精确评估,偏差一般在2 μGal以内,不确定度约为0.1 μGal[21]。
8) 落体旋转
自由落体式的绝对重力仪中,在落体被释放的一瞬间,受释放机构转矩作用的影响,落体将具有微小的旋转角速度,在下落过程中具有微小旋转。落体自由下落期间只受重力作用,因此它将绕自身质心旋转,且角速度基本保持不变;但由于落体内的部件在加工和装配过程中存在误差,落体的质心不可能与其内部角锥棱镜的光心完全重合。激光干涉仪实际测量的重力值是角锥棱镜光心的运动加速度,其中包含了光心具有的向心加速度。落体质心与光心的距离矢量可以分解为水平方向与垂直方向两部分,水平矢量引起的向心加速度沿水平方向,与重力方向垂直,不会对测量造成干扰;但垂直分量会引入重力测量偏差和不确定度。由于每次下落的释放状态随机,落体的真实角速度并不统一,因此一般不讨论该项引入的偏差。相应的不确定度可以表示为
式中:ωmax为落体可能的最大旋转速度;Rmax为落体光心和质心间距的最大值。通过精确设计和调节真空腔内的传动机构的控制参数以及落体的加工和装配过程,可以保证ωmax和Rmax在较小的范围内。一般而言,该项引入的测量不确定度为0.5 μGal左右。
9) 地面反弹效应
在落体被释放的一瞬间,绝对重力仪的静态质量突然减小,导致仪器受到的地面支持力大于其自身的重力,称为地面反弹效应。这种效应本质上是一种发生在垂直方向上的受力不平衡,等效为一个瞬时冲击,可能导致隔振系统内的参考棱镜或激光干涉仪内的光学器件(包括分光镜等)产生微振动,从而对落体轨迹的测量引入误差。实际上该效应可能是高精度绝对重力测量不确定度的主要来源之一[5]。地面反弹效应同样具有随机性,因此也仅对该项的不确定度进行讨论。根据相关文献,利用下式可以从计算重力测值时的二次拟合残差(一般近似于正弦波形)中推导出地面反弹效应引入的不确定度[22],一般在0.1 μGal左右[5]。
式中:T为干涉条纹信号的持续采集时间;A0和f0分别为拟合残差波形的幅值和频率。
10) 有效高度的计算误差
如前所述,落体的释放状态具有随机性,因此落体每次的下落轨迹都不可能完全相同,使得测量开始时刻(第一个采集到的数据点)对应的落体初速度v0和重力加速度真值g0都可能有所改变,从而在有效高度的计算中引入误差。通过实验测定测量开始时刻的落体速度不确定度和重力不确定度,将其代入有效高度的计算公式,配合测量点所在地的测量不确定度重力梯度,可以得到该项一般也为0.1 μGal左右。
总体而言,来源于上述大多数误差源的不确定度都可以通过使用精度更高的配套设备来降低,包括铷原子钟、激光器、光电探测器、光束准直器等;部分误差源造成的影响需要依靠研究人员的设计、调试及操作来降低,如光束垂直度引起的不确定度可以通过调节干涉仪的可动器件来减小,利用分子泵和离子泵进一步提高真空腔内的真空度可以减小空气阻尼,通过更为合理的机械设计可以减小温度梯度、落体旋转、仪器自吸引、地面反弹效应、有效高度计算等因素导致的不确定度;还有部分因素,如光速有限,可能需要更为复杂的软件算法才能降低其影响。目前正在使用中的各类绝对重力仪在配套设备精度的选择上基本一致,仪器自身不确定度的差异主要来源于研究人员在光机电设计方面的不同,由此导致每种类型重力仪的主导不确定度因素也有所不同,但基本上各类仪器的合成标准不确定度保持在10 μGal以内。
2.3.2 来源于测量环境的误差
1) 大气压强
测量点所在地大气压强的增大等效于测量点上方大气质量的增大,使得大气对落体的万有引力增大,导致重力测值小于重力真值。可以通过引入气压影响因子fB = 0.3 μGal/hPa来修正此项影响,
式中:p0表示重力测量点的实测气压,由具有不确定度的气压计测量得到,因此会引起重力测值的不确定度,一般为1~3 μGal。定期校准气压计可有效减小此项不确定度。
式中:pn为标准大气压强,满足其中psea=1013.25 hPa为海平面标准大气压;hm为测量点所在地的海拔高度;L=0.0065 K/m为垂直温度梯度;T0=288.15 K为海平面标准温度;g=9.80665 m/s2为模型中所取地表重力加速度平均值;M=0.0289644 kg/mol为干空气的分子质量;R=8.31447 J/(mol ·K)为理想气体常数。
2) 固体潮汐
受太阳和月球对地球周期性变化的引力影响,地球重力场也具有周期性的变化。这种变化称为重力潮汐,其中造成地球形变的潮汐力称为固体潮汐,由此引起的重力加速度的变化一般在±150 μGal之间。因此进行绝对重力测量时通常会进行潮汐修正,具体的修正值可以通过专业软件如Tsoft等来计算。在一定的精度范围内,固体潮汐可以作为评估重力测值精密度的参考标准,即潮汐修正前的重力测值越靠近潮汐的理论变化曲线,则绝对重力仪的测量精密度越高。目前此项引入的测量不确定度仅为0.1 μGal左右。
3) 海洋负荷
太阳与月球的引力除引起固体潮汐外还会引发海潮。海潮意味着地球表面的质量分布存在明显变化,因此也会引起重力值的变化,这种现象称为海洋负荷。海洋负荷引起的不确定度一般为0.3 μGal左右,实际的重力变化幅度从近海至内陆逐渐递减。一般而言,现有的绝对重力仪进行24 h的重复测量就可以将海洋负荷的影响降至最小。
4) 电磁力
电磁力对绝对重力测量的影响可以根据其物理性质来分别讨论。首先,如果制作落体时采用了磁性材料,磁场效应将引起较大的重力测量偏差,因此现有绝对重力仪中落体的零部件一般均采用非磁性材料来制作,此时由外界磁场对落体的吸引和排斥非常小,该项测量偏差可以忽略不计。其次,制作落体的材料主要为各类金属,均属于导电材料,因此在磁场环境下落体的自由下落将产生涡流效应,形成与运动方向相反的阻力。构成该磁场环境的磁场源包括真空腔外的电机、离子泵中的磁铁等,形成的磁场在落体的实际运动范围内通常与地磁场在同一量级,所以目前主要的激光干涉式绝对重力仪在评估时通常忽略该项不确定度。然后,主要由金属材料构成的真空腔外壳将在落体外形成法拉第笼,使得真空腔外的电场不会对落体产生静电力,真空腔内也没有电场。最后,在自由落体式绝对重力仪中,只要采用同种材料制作落体与托盘上相互接触的部位(如前所述,托盘是带动落体上升至释放点并承接下落后的落体的传动部件),就可以保证两者之间几乎不存在接触电压,因此相应的静电力也可以忽略不计。
5) 极移
地球的自转轴方向在最大惯性轴附近几米范围内进行周期性的运动,从而导致的最大重力值变化可达13 μGal。极移包含两个频率分量,周期分别为12个月和14个月,前者主要来源于地球作为一个刚体在大气环流影响下的受迫摆动,后者相当于地球作为一个弹性体的自由摆动[23-24]。重力测量的极移修正公式(单位:nm·s-2)由国际绝对重力基准网的绝对观测数据处理标准给出。
式中:ω为地球自转角速度(单位:rad/s);a=6378136.6(1)为水准椭球半长轴;φ和λ分别为地理纬度和地理经度(单位:rad);xp和yp是极坐标(单位:rad),由国际地球自转服务(IERS)网站提供[25]。该项改正一般用于计算重力测值的修正值,引入的测量不确定度很小,可取为0.1 μGal以内。
6) 科里奥利力
受地球自转的影响,如果落体的速度在自由下落的过程中具有东西方向的水平度分量,则落体受到的科里奥利力会有垂直方向上的分量,由此引起的重力测量偏差为
式中:ΩN为地球的标称角速度;vEW为落体速度在水平面东西方向上的分量;φ为重力测点的地理纬度。该项引起的测量不确定度一般约为0.3 μGal。
3 原子干涉绝对重力仪
3.1 基本原理
原子干涉绝对重力仪基于冷原子物质波干涉原理,利用激光脉冲作用到冷原子团,通过双光子受激拉曼跃迁或多光子布拉格衍射过程[26-27],实现原子波包的分束合束操作,完成原子干涉过程。本文以受激拉曼跃迁过程为例,如
图 8. 基本原理。(a)受激拉曼跃迁;(b)原子干涉示意图
Fig. 8. Basic principle. (a) Stimulated Raman transition; (b) schematic of atomic interference
原子干涉过程一般采用三个拉曼激光脉冲作用,如
拉曼激光脉冲与原子作用时,两束激光的相位差将映射到原子物质波函数,则经过三个脉冲后,两条路径上原子波包的相位差为
式中:Φi表示第i个脉冲时刻激光相位差。另一方面,原子波包沿着两条路径自由演化的相位差,利用路径积分方法计算的结果为0。因此原子干涉仪的相位即
式中:g表示待测重力加速度;T表示脉冲间隔时间;α表示拉曼激光频率扫描速率,为了补偿原子下落过程的多普勒效应。因此,通过改变α值,可以改变干涉相位,得到原子干涉条纹。
原子干涉绝对重力仪可以通过拉曼脉冲作用到冷原子实现重力加速度的测量,整个系统的工作原理如
图 9. 原子干涉绝对重力仪的工作原理示意图
Fig. 9. Schematic of working principle of atom interferometry absolute gravimeter
按照原子运动的路径,可以将原子干涉绝对重力仪分为喷泉式和自由下落式两种类型。本文以自由下落式原子干涉绝对重力仪为例,如
3.2 技术实现
原子干涉绝对重力仪按照功能可以分为真空物理系统、光学系统、控制系统和隔振等部分[30]。真空物理系统用于制备冷原子样品,激光系统用于实现激光冷却、原子干涉操控和探测等,控制系统用于完成原子干涉绝对重力仪整机工作流程,隔振及其他辅助部件用于实现系统稳定工作。各部分的组成和要求如
图 11. 原子干涉绝对重力仪系统结构框图[30]
Fig. 11. System structure diagram of atom interferometry absolute gravimeter[30]
3.2.1 真空物理系统
原子干涉绝对重力仪的真空物理系统包括真空腔、磁屏蔽、磁场线圈等部分。
图 12. NIM-AGRb-1型原子干涉绝对重力仪的真空系统[29]
Fig. 12. Vacuum system of NIM-AGRb-1 atom interferometry absolute gravimeter [29]
自由下落式原子干涉绝对重力仪真空物理系统结构紧凑,容易实现小型化[31-32]。喷泉式原子干涉绝对重力仪的真空物理部分稍有不同,原子干涉区位于磁光阱的上部,可以实现更长的拉曼脉冲间隔时间T,提高重力测量灵敏度[33-34]。
法国巴黎天文台提出一种简化真空物理系统的方案,采用金字塔形状的角锥反射镜,放置于真空系统中,可以实现单束激光制备冷原子、拉曼干涉和探测等过程,省去多个真空窗口,增强真空系统的稳定性[35]。法国MuQuans公司采用了这种方案,研制实现商品原子干涉绝对重力仪[36]。美国加州理工学院伯克利分校Müller小组也采用了类似的方案,有利于实现可移动型原子干涉绝对重力仪[37]。
3.2.2 光学系统
原子干涉绝对重力仪所需激光一般由单台或多台半导体激光器产生,以铷原子为例,对应的激光波长为780 nm,需要5种以上的频率,分别为冷却激光、重泵激光、探测激光、拉曼主激光、拉曼从激光等。通常采用饱和吸收谱或调制转移谱等方案,稳定参考激光频率,再通过声光调制器移频产生各种频率的激光。拉曼激光可以采用电光调制器或光学锁相实现。
图 13. NIM-AGRb-1型原子干涉绝对重力仪的光学系统示意图
Fig. 13. Schematic of optical system of NIM-AGRb-1 atom interferometry absolute gravimeter
近年来,各研究单位相继开展了小型化光学系统的研制[38-42]。法国航天局采用一台窄线宽1560 nm光纤激光器,经过相位调制(1 GHz)后利用周期极化铌酸锂晶体(PPLN)倍频产生780 nm激光,通过一级边带的倍频激光锁定到85Rb饱和吸收谱线上,而载波的倍频光可用于产生冷却激光和探测激光,再利用另一个相位调制器(6.5 GHz),一级边带倍频产生重泵激光,如
3.2.3 隔振
振动噪声会通过拉曼激光反射镜影响激光与原子相互作用时的相位,引起原子干涉相位的变化,是原子干涉绝对重力仪最重要的噪声来源之一。常用的隔振方案包括被动隔振装置、主动隔振装置、数据修正等。NIM-AGRb-1型原子干涉绝对重力仪采用基于LaCoste弹簧的被动隔振装置,时间常数达到20 s以上,具有结构简单,维护方便等优点[44]。美国斯坦福大学[45]、德国洪堡大学[46]、中国科学院精密测量科学与技术创新研究院(APM)[47]、华中科技大学(HUST)[48]等单位采用主动隔振方案,一般通过被动隔振平台上放置测振计测量剩余振动噪声,并通过滤波环路反馈到执行元件,实现振动噪声的主动伺服抑制。该方案可以将低频噪声(1~10 Hz)降低1个量级左右,效果良好。法国巴黎天文台[49]、浙江工业大学(ZJUT)[50]等单位采用了数据修正的方案,利用测振计或加速度计直接测量平台剩余振动噪声,经过数据处理后直接用于原子干涉相位的修正,通过软件方式降低振动噪声对重力测量的影响。
图 14. 小型光学系统示意图。(a)法国航天局小组研制的小型光学系统示意图[39];(b)激光频率分布
Fig. 14. Schematic of a small optical system. (a) Schematic of a small optical system developed by the French space agency group[39]; (b) laser frequency distribution
3.3 重力测量性能评估
3.3.1 测量灵敏度
原子干涉绝对重力仪可用于长期连续测量重力,如
图 15. 长期连续测量数据。(a) NIM-AGRb-1原子干涉绝对重力仪长期重力测量数据; (b)实验数据与理论差值[28]
Fig. 15. Long term gravity measurement data. (a) Long term gravity measurement data of NIM-AGRb-1 atom interferometry absolute gravimeter; (b) experimental data and theoretical difference [28]
影响原子干涉绝对重力仪测量灵敏度的因素包括振动噪声、拉曼激光的频率、相位和强度噪声、探测噪声等[49]。
表 1. 国内外部分研究小组的原子重力仪灵敏度
Table 1. Sensitivity of atomic gravimeters of some research groups at home and abroad
|
3.3.2 系统误差评估
原子干涉绝对重力仪的另一个重要参数是准确度,各种物理因素会导致重力测量绝对值的偏移,因此需要进行系统误差评估和不确定度分析。一般将系统误差分为两类,一类与拉曼激光的波矢方向无关,包括二阶塞曼效应、单光子相移、射频相移等,可以通过改变波矢方向降低影响;另一类与波矢方向有关,包括双光子相移、自引力效应、科里奥利力效应、拉曼激光垂直偏差、激光波长偏差、波前畸变等效应,需要进行理论或实验研究。此外,还需要考虑各种环境因素,如气压变化、极移等效应。
表 2. NIM-AGRb-1系统误差评估表[28]
Table 2. NIM-AGRb-1 system error evaluation table[28]
|
波前畸变效应来源于拉曼激光受到光学元件的影响,波前并非完美的平面,与原子团作用时引起重力测量值的偏移。为了评估该效应引起的系统误差,可以通过改变原子温度[52]或选择不同运动轨迹的原子[53],不确定度一般为3~4 μGal量级。2018年,法国巴黎天文台小组报道采用偶极阱蒸发冷却实现原子温度冷却至50 μK,通过改变原子温度,测量重力加速度数值,外推至原子温度为0,评估正常原子温度的偏移,实现波前畸变效应不确定度降至1.3 μGal[52]。采用该方案,有望将原子干涉绝对重力仪的准确度进一步提高。
4 结束语
目前,激光干涉绝对重力仪还是主流应用的仪器类型,占据着全球大部分市场份额。激光干涉绝对重力仪测量结果的合成标准不确定度最优达到1.8 μGal,相关的不确定度分析表明,这一水平已逼近测量极限精度。未来,激光干涉绝对重力仪将向着进一步提高其可靠性,并扩展其在海洋、航空、卫星等多领域多平台的应用发展。原子干涉绝对重力仪是一种基于冷原子物质波干涉原理的新型绝对重力仪,具有准确度高、稳定度好、测量速度快、无机械磨损等特点,是绝对重力测量非常有潜力的技术方案和发展方向,其标准不确定度已经接近激光干涉绝对重力仪,随着波前畸变效应评估水平的提高,有望进一步降低至1 μGal水平。自2009年开始,原子干涉绝对重力仪连续参与了绝对重力仪国际比对,特别是 2017年在北京由中国计量科学研究院主办的第10届绝对重力仪国际比对中,来自国内6家单位研制的原子干涉绝对重力仪参加了研究比对,实现了与激光干涉绝对重力仪的相互验证[54]。近年来,可移动型原子干涉绝对重力仪发展迅速,实现了船载和机载的重力测量实验[55-56],国外多家公司已实现商品化,这些进展为原子干涉绝对重力仪向着地球物理、资源勘探、地震监测、惯性导航、****等领域的应用奠定基础。最后,不论是激光干涉绝对重力仪还是原子干涉绝对重力仪,其发展都离不开激光器的发展,也离不开激光干涉、激光稳频、激光调控、激光冷却等激光技术的发展。应该说,激光技术六十多年的发展进步,为绝对重力仪的技术发展奠定了基础,提供了机遇。
致谢 本文参考了中国计量科学研究院时间频率计量科学研究所光学重力组和原子重力组的工作,本文的撰写得到了庄伟博士和冯金扬博士的大力帮助,一并表示感谢。
[1] Peck E R. Theory of the corner-cube interferometer[J]. Journal of the Optical Society of America, 1948, 38(12): 1015-1024.
[2] Klopping F. FG5X absolute gravimeter user's manual[EB/OL]. ( 2015-10-26)[2020-07-24]. http:∥microglacoste.com/wp-content/uploads/2018/01/FG5X-Manual-115060001.pdf.
[3] Zumberge M A, Rinker R L, Faller J E, et al. A portable apparatus for absolute measurements of the earth's gravity[J]. Metrologia, 1982, 18(3): 145-152.
[4] Rinker RL. Super spring- a new type of low frequency vibration isolator[D]. Boulder: University of Colorado at Boulder, 1983.
[6] Klopping F. FG5X absolute gravity meters[EB/OL]. ( 2014-12-01)[2020-07-24]. http:∥microglacoste.com/product/fg5-X-absolute-gravimeter/.
[8] D'Agostino G. Desogus S , Germak A , et al. The new IMGC-02 transportable absolute gravimeter: measurement apparatus and applications in geophysics and volcanology[J]. Annals of Geophysics, 2008, 51: 39-49.
[9] Faller J E, Vitouchkine A L. A new small cam-driven absolute gravimeter[J]. International Association of Geodesy Symposia, 2005, 129: 276-279.
[10] Feng Y Y, Zhang G Y, Li D X, et al. A transportable absolute gravimeter for determining the acceleration due to the earth's gravity[J]. Metrologia, 1982, 18(3): 139-143.
[13] Tsubokawa T, Svitlov S. New method of digital fringe signal processing in an absolute gravimeter[J]. IEEE Transactions on Instrumentation and Measurement, 1999, 48(2): 488-491.
[14] Vitushkin L F, Orlov O A. Absolute ballistic gravimeter ABG-VNIIM-1 by D.I. Mendeleyev research institute for metrology[J]. Gyroscopy and Navigation, 2014, 5(4): 283-287.
[15] RothleitnerC. Ultra-high precision, absolute, earth gravity measurements [D]. Erlangen: University ofErlangen-Nuremberg, 2008.
[16] Kennard EH. Kinetic theory of gases, with an introduction to statistical mechanics[M]. New York:McGraw-Hill, 1938.
[17] D'AgostinoG, Desogus S, Germak A, et al. The assessment of the measurement error due to a non-vertical laser beam path in absolute gravimeters[C]∥Cahier du Centre Europeen dé Geodynamique et de Séismologie, 2006: 26.
[22] Svetlov S M. An absolute gravimeter and vibration disturbances: a frequency responses method[J]. Springer Berlin Heidelberg, 1997: 47-54.
[23] Wahr J M. Deformation induced by polar motion[J]. Journal of Geophysical Research Atmospheres, 1985, 90(B11): 9363-9368.
[24] Chandler S C. On the variation of latitude, I[J]. The AstronomicalJournal, 1891, 11: 59- 61.
[25] McCarthy DD, Petit G. IERSconventions ( 2003)[ R]. Frankfurt am Main: IERS, 2004.
[26] Kasevich M, Chu S. Atomic interferometry using stimulated Raman transitions[J]. Physical Review Letters, 1991, 67(2): 181-184.
[27] Giltner D M. McGowan R W, Lee S A. Atom interferometer based on Bragg scattering from standing light waves[J]. Physical Review Letters, 1995, 75(14): 2638-2641.
[28] Dalibard J, Cohen-Tannoudji C. Laser cooling below the Doppler limit by polarization gradients: simple theoretical models[J]. Journal of the Optical Society of America B, 1989, 6(11): 2023-2045.
[29] Wang S K, Zhao Y, Zhuang W, et al. Shift evaluation of the atomic gravimeter NIM-AGRb-1 and its comparison with FG5X[J]. Metrologia, 2018, 55(3): 360-365.
[32] Bidel Y, Carraz O, Charrière R, et al. Compact cold atom gravimeter for field applications[J]. Applied Physics Letters, 2013, 102(14): 144107.
[33] Hu Z K, Sun B L, Duan X C, et al. Demonstration of an ultrahigh-sensitivity atom-interferometry absolute gravimeter[J]. Physical Review A, 2013, 88(4): 043610.
[35] Bodart Q, Merlet S, Malossi N, et al. A cold atom pyramidal gravimeter with a single laser beam[J]. Applied Physics Letters, 2010, 96(13): 134101.
[37] Wu XJ, PagelZ, Malek BS, et al., 2019, 5(9): eaax0800.
[38] Schmidt M, Prevedelli M, Giorgini A, et al. A portable laser system for high-precision atom interferometry experiments[J]. Applied Physics B, 2011, 102(1): 11-18.
[39] Theron F, Carraz O, Renon G, et al. Narrow linewidth single laser source system for onboard atom interferometry[J]. Applied Physics B, 2015, 118(1): 1-5.
[40] Wang Q Y, Wang Z Y, Fu Z J, et al. A compact laser system for the cold atom gravimeter[J]. Optics Communications, 2016, 358: 82-87.
[43] Caldani R, Merlet S. Pereira dos Santos F, et al. A prototype industrial laser system for cold atom inertial sensing in space[J]. The European Physical Journal D, 2019, 73(12): 248-256.
[46] Hauth M, Freier C, Schkolnik V, et al. First gravity measurements using the mobile atom interferometer GAIN[J]. Applied Physics B, 2013, 113: 49-55.
[47] Tang B, Zhou L, Xiong Z Y, et al. A programmable broadband low frequency active vibration isolation system for atom interferometry[J]. Review of Scientific Instruments, 2014, 85(9): 093109.
[48] Zhou M K, Xiong X, Chen L L, et al. Note: a three-dimension active vibration isolator for precision atom gravimeters[J]. Review of Scientific Instruments, 2015, 86(4): 046108.
[51] Gillot P, Francis O, Landragin A, et al. Stability comparison of two absolute gravimeters: optical versus atomic interferometers[J]. Physics, 2014, 51(5): L9-L11.
[53] Zhou M K, Luo Q, Chen L L, et al. Observing the effect of wave-front aberrations in an atom interferometer by modulating the diameter of Raman beams[J]. Physical Review A, 2016, 93: 043610.
Article Outline
吴书清, 李天初. 绝对重力仪的技术发展: 光学干涉和原子干涉[J]. 光学学报, 2021, 41(1): 0102002. Shuqing Wu, Tianchu Li. Technical Development of Absolute Gravimeter: Laser Interferometry and Atom Interferometry[J]. Acta Optica Sinica, 2021, 41(1): 0102002.