CN114137626B - 一种通过调制拉曼激光确定原子重力仪等效高度的方法 - Google Patents
一种通过调制拉曼激光确定原子重力仪等效高度的方法 Download PDFInfo
- Publication number
- CN114137626B CN114137626B CN202111439466.8A CN202111439466A CN114137626B CN 114137626 B CN114137626 B CN 114137626B CN 202111439466 A CN202111439466 A CN 202111439466A CN 114137626 B CN114137626 B CN 114137626B
- Authority
- CN
- China
- Prior art keywords
- atomic
- coordinate system
- raman
- formula
- height
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000001069 Raman spectroscopy Methods 0.000 title claims abstract description 92
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000007704 transition Effects 0.000 claims abstract description 30
- 230000010355 oscillation Effects 0.000 claims abstract description 16
- 230000005484 gravity Effects 0.000 claims description 33
- 238000005259 measurement Methods 0.000 claims description 17
- 238000005305 interferometry Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本申请涉及一种通过调制拉曼激光确定原子重力仪等效高度的方法。所述方法包括:通过由载波和经电光调制器调制后的产生的±1阶边带组成的拉曼激光对驱动原子双光子受激拉曼跃迁,在原子不同的自由下落时刻扫描拉曼光的脉冲持续时间,测量记录每个拉曼脉冲持续时间和每个时刻对应的原子跃迁概率;在每个时刻拟合拉比振荡概率曲线,得到每个时刻对应的有效拉比频率,再对每个时刻对应的有效拉比频率进行拟合,得到原子运动坐标系水平面和仪器坐标系水平面之间的高度差,进而得到原子重力仪相对于参考水平的等效点高度。本发明可以简单而精确地计算原子重力仪相对于特定参考水平的等效高度。
Description
技术领域
本申请涉及原子重力仪技术领域,特别是涉及一种通过调制拉曼激光确定原子重力仪等效高度的方法。
背景技术
近二十年来,基于物质波干涉仪的双光子激发拉曼跃迁原子重力仪发展迅速,已成为一种应用广泛的强大仪器。随着1μGal(10-8m/s2)级的绝对测量原子重力仪逐渐成为重力的基准。近年来,在车、船、飞机上搭载的原子重力仪在野外应用方面表现出很强的适应性和广阔的前景。
在地球上进行测量时,原子重力仪存在一个棘手的问题。由于垂直重力梯度的存在,由具有自由落体原子的重力仪测量的绝对重力值不能被视为常数。因此测量的沿原子运动轨迹的平均重力,可以等于某一点的绝对重力,称为等效点,其相对于参考水平面的距离称为等效高度。在不同绝对重力仪的比较中,所有仪器的参考水平面是统一的,这意味着在考虑约300μGal/m的垂直重力梯度时,1cm的等效高度误差将导致3μGal的重力误差。因此,准确确定原子重力仪的等效高度将显着改善其测量***误差。
然而,原子自由落体坐标系与仪器坐标系是相对独立,因而难以确定等效点在仪器坐标系或其他实体坐标系中的位置。如何在两个坐标系之间建立正确的联系来计算等效高度变得至关重要。在典型的原子重力仪中,有以下几种解决方案:
一是使用一个物理定义而不是具体的对象来建立连接。再次根据其他具体对象间接估计物理定义对象的位置时,不可避免地会产生误差。一般来说,制造和安装几乎不可能与设计完全相同。特别是,由于磁场和激光的不匹配等原因,实际的MOT中心会发生偏离。
二是考虑通过具体对象直接确定等效高度。现有技术公开了一种使用原子干涉直接测量原子团和底部反射镜之间距离的方法。通过在π脉冲处对拉曼激光器施加频率跳跃,会引起与距离线性相关的额外干涉相位,获得2mm的测量不确定度。这种方法精确但难以实现。首先,不适当的频率跳跃会导致不必要的额外相移和拉曼激光功率的变化。此外,较长的跳频响应时间不易缩短,这与由干涉区域长度所限制的较短干涉时间相冲突。此外,等效高度的测量精度取决于原子干涉仪本身的优异性能,这并非所有原子重力仪都具备。
发明内容
基于此,有必要针对上述技术问题,提供一种能够提高确定原子重力仪等效高度的精确度和适应性方法,即本专利所述的的通过调制拉曼激光确定原子重力仪等效高度的方法。
一种通过调制拉曼激光确定原子重力仪等效高度的方法,所述方法包括:
在拉曼脉冲原子干涉重力仪中,以MOT中心为原点建立原子运动坐标系,同时以反射镜上表面为XY平面,Z轴与所述原子运动坐标系一致,建立仪器坐标系;
通过由载波和经电光调制器调制后的产生的±1阶边带组成的拉曼激光对驱动原子双光子受激拉曼跃迁,在原子不同的自由下落时刻扫描拉曼光的脉冲持续时间,测量记录每个拉曼脉冲持续时间和每个时刻对应的原子跃迁概率;其中,每个原子下落时刻扫描一系列拉曼光脉冲持续时间所获得的原子跃迁概率对应该时刻原子拉比振荡曲线;
根据预设的第一公式,在每个时刻拟合拉比振荡概率曲线,得到每个时刻对应的有效拉比频率;
根据预设的第二公式,对所述每个时刻对应的有效拉比频率进行拟合,得到所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差;
根据计算得到原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标、所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差,以及预设的参考水平面与仪器坐标系的高度,得到原子重力仪相对于所述参考水平面的等效点高度。
在其中一个实施例中,还包括:所述第一公式为:
其中,dam为中间变量,dam=zm,A-za,A,zm,A为所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差,za,A为ta时刻原子团在所述原子运动坐标系中的z轴坐标,τ为拉曼脉冲持续时间,P(dam,τ)为参数dam和参数τ对应的原子跃迁概率。
在其中一个实施例中,还包括:将所述第一公式进一步优化,得到第一公式优化后的拟合公式:
P(τ)=b1[1-exp(-b2τ)cos(Ωτ)]+b4
其中,b1、b2、Ω、b4均为待拟合的参数;
通过所述第一公式优化后的拟合公式在每个时刻ta拟合拉比振荡概率曲线P(τ),得到每个时刻ta对应的有效拉比频率Ω(ta)。
在其中一个实施例中,还包括:所述第二公式为:
其中,Ω0为拉曼激光组中只有E0和E+1时的有效拉比频率,Δk=kn+1-kn为拉曼激光对的波矢差,Δ0,+1为E0与E+1组成的拉曼激光对驱动原子双光子拉曼跃迁时的单光子失谐,Δ-1,0为E-1与E0组成的拉曼激光对驱动原子双光子拉曼跃迁时的单光子失谐。
在其中一个实施例中,还包括:将所述第二公式进一步优化,得到第二公式优化后的拟合公式:
其中,ΩA、zm,A、υ0、ΩB均为待拟合参数,g0为测得的当地绝对重力值;
根据所述第二公式优化后的拟合公式,以ta自变量,以有效拉比频率Ω(ta)为因变量进行拟合,得到所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差拟合值
在其中一个实施例中,还包括:根据所述原子运动坐标系和所述仪器坐标系之间的高度差拟合值,通过公式/>计算zm,A的真实值,在此需要提前预估zm,A的值以设定调节参数i,其中zm,A的预估值与真实值之间的误差不大于Td。
在其中一个实施例中,还包括:根据计算得到原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标、所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差,以及预设的参考水平面与仪器坐标系的高度,得到原子重力仪相对于所述参考水平面的等效点高度为:
heq=zm,A-zeq,A-dmr
其中,zeq,A为计算得到原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标,计算式为v1表示原子团在自由下落至t1时刻且获得第一个拉曼脉冲的反冲速度vrec转移后的速度,T表示拉曼脉冲之间的间隔时间,z1,A表示t1时刻原子团的位置,dmr为所述预设的参考水平面与仪器坐标系水平面的高度差。
在其中一个实施例中,还包括:所述拉曼激光对驱动的原子为87Rb原子。
在其中一个实施例中,还包括:所述拉曼脉冲中的脉冲为π/2-π-π/2序列。
上述通过调制拉曼激光确定原子重力仪等效高度的方法,通过由载波和经电光调制器调制后的产生的±1阶边带组成的拉曼激光对驱动原子双光子受激拉曼跃迁,在原子不同的自由下落时刻扫描拉曼光的脉冲持续时间,测量记录每个拉曼脉冲持续时间和每个时刻对应的原子跃迁概率;在每个时刻拟合拉比振荡概率曲线,得到每个时刻对应的有效拉比频率,再对每个时刻对应的有效拉比频率进行拟合,得到原子运动坐标系水平面和仪器坐标系水平面之间的高度差;根据计算得到原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标、原子运动坐标系水平面和仪器坐标系水平面之间的高度差,以及预设的参考水平面与仪器坐标系的高度,得到原子重力仪相对于参考水平面的等效点高度。本发明通过调制拉曼激光器,与受激拉曼跃迁相关的有效拉比频率不再是常数,而是垂直方向上的周期性余弦函数,使得可以简单而精确地计算原子重力仪相对于特定参考水平的等效高度。
附图说明
图1为一个实施例中通过调制拉曼激光确定原子重力仪等效高度的方法的流程示意图;
图2为一个实施例中典型拉曼脉冲原子干涉重力仪的基本原理图;
图3为一个实施例中在某一ta时刻拟合的拉比振荡曲线P(τ);
图4为一个实施例中对Ω(ta)进行拟合得到的曲线。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
在一个实施例中,如图1所示,提供了一种通过调制拉曼激光确定原子重力仪等效高度的方法,包括以下步骤:
步骤102,在拉曼脉冲原子干涉重力仪中,以MOT中心为原点建立原子运动坐标系,同时以反射镜上表面为XY平面,Z轴与原子运动坐标系一致,建立仪器坐标系。
本专利的理论推导具体如下:
典型拉曼脉冲原子干涉重力仪的基本原理图如图2所示。在这个模型中,建立了两个坐标系,两个坐标系的Z轴是相同的,其中一个是以MOT中心为原点的原子运动坐标系O-XAYAZA,另一个是以反射镜上表面为XIYI平面的仪器坐标系O-XIYIZI。通过中间变量dam=zm,A-za,A建立两个坐标系之间的连接,其中zm,A是镜子在O-XAYAZA坐标系中的位置,也就是原子运动坐标系水平面和仪器坐标系水平面之间的高度差,而表示在时间ta时原子团的坐标,g0是zA=0处的重力,v0表示原子团在t0时刻的初速度。
用π/2-π-π/2的序列分析拉曼脉冲原子重力仪,其中,在O-XAYAZA坐标系中测得的绝对重力可以表示为
gm=g0+γ(7g0T2/12+v1T+z1,A)=g0+γzeq,A
其中原子自由落体从t=0开始,γ和T分别表示线性垂直重力梯度和拉曼脉冲之间的间隔时间,v1表示原子团在自由下落至t1时刻且获得第一个拉曼脉冲的反冲速度vrec转移后的速度,此时原子团的位置是z1,A。所以,原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标是所以原子干涉重力仪的测量点相对于特定参考水平面的等效高度的计算公式为
heq=zm,A-zeq,A-dmr
其中,dmr为所述预设的参考水平面与仪器坐标系水平面的高度差。其中关键的未知项是原子运动坐标系和仪器坐标系之间沿Z轴的垂直距离,即zm,A。
下面论述确定zm,A的具体步骤。
步骤104,通过由载波和经电光调制器调制后的产生的±1阶边带组成的拉曼激光对驱动原子双光子受激拉曼跃迁,在原子不同的自由下落时刻扫描拉曼光的脉冲持续时间,测量记录每个拉曼脉冲持续时间和每个时刻对应的原子跃迁概率。其中,每个原子下落时刻扫描一系列拉曼光脉冲持续时间所获得的原子跃迁概率对应该时刻原子拉比振荡曲线。
具体地,通过由载波和经电光调制器(EOM)调制后产生的边带组成的拉曼激光对驱动原子双光子拉曼跃迁。期望的频率差为ωm的拉曼激光对由调制后产生的En和En+1边带组成,其中单光子失谐为Δn,n+1。由于高阶边带的幅度显着降低,可以只考虑了载波E0和E±1的影响。因此,不同位置的有效拉比频率Ω可以表示为第二公式,为:
其中,Ω0为拉曼激光组中只有E0和E+1时的有效拉比频率,Δk=kn+1-kn为拉曼激光对的波矢差,Δ0,+1为E0与E+1组成的拉曼激光对驱动原子双光子拉曼跃迁时的单光子失谐,Δ-1,0为E-1与E0组成的拉曼激光对驱动原子双光子拉曼跃迁时的单光子失谐;显然,Ω以余弦形式随dam变化,位置上的周期Td=π/Δk。在本具体实施例中,使用的87Rb原子,ωm=6.834GHz,Td的值为2.19cm。
在原子自由下落时间ta处对原子施加脉冲拉曼光,拉曼光的持续时间为τ,同时测量记录该条件的原子跃迁概率P。在实验中,扫描参数ta和τ来获取数据P(ta,τ),即在不同的时刻ta的原子拉比振荡曲线P(τ)。
下一步是根据预设的第一公式,在每个时刻ta拟合拉比振荡曲线P(τ)获得Ω(ta)的值。
步骤106,根据预设的第一公式,在每个时刻拟合拉比振荡概率曲线,得到每个时刻对应的有效拉比频率。
具体地,第一公式为:
其中,dam为中间变量,dam=zm,A-za,A,zm,A为所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差,za,A为ta时刻原子团在所述原子运动坐标系中的z轴坐标,τ为拉曼脉冲持续时间,P(dam,τ)为参数dam和参数τ对应的原子跃迁概率。
在实际中,由于振荡幅值随拉曼光持续时间衰减的效应,实际拟合公式可进一步优化如下
P(τ)=b1[1-exp(-b2τ)cos(Ωτ)]+b4
其中b1,b2,b4,和Ω均为待拟合的参数。
步骤108,根据预设的第二公式,对每个时刻对应的有效拉比频率进行拟合,得到原子运动坐标系水平面和仪器坐标系水平面之间的高度差。
根据预设第二公式,对Ω(ta)进行拟合,以得到参数zm,A。
具体地,第二公式为:
其中,Ω0为拉曼激光组中只有E0和E+1时的有效拉比频率,Δk=kn+1-kn为拉曼激光对的波矢差,Δ0,+1为E0与E+1组成的拉曼激光对驱动原子双光子拉曼跃迁时的单光子失谐,Δ-1,0为E-1与E0组成的拉曼激光对驱动原子双光子拉曼跃迁时的单光子失谐。
在实验中,以中间变量dam=zm,A-za,A对该拟合公式进行优化如下
其中,ΩA、zm,A、υ0、ΩB均为待拟合参数,g0为测得的当地绝对重力值,需要指明,该公式中对绝对重力值g0的精度要求较低。
由于在预设第二公式中,Z轴上位置上具有空间周期性,周期为Td,因此在上一步骤中,拟合得到的并非zm,A的真值,二者通过如下公式进行计算
其中,i为整数,它的值通过提前预估zm,A的值来确定,预估的误差不超过Td即可。
在上述两次拟合中,本实施例的拟合结果图3,4。
步骤110,根据计算得到原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标、原子运动坐标系水平面和仪器坐标系水平面之间的高度差,以及预设的参考水平面与仪器坐标系的高度,得到原子重力仪相对于参考水平面的等效点高度。
具体地,计算原子干涉重力仪的等效高度前面已经给出,计算公式如下
heq=zm,A-zeq,A-dmr
在本专利具体实施例中,预设dmr=0,计算得到等效高度测量的不确定度为2.1×10-4m。
上述通过调制拉曼激光确定原子重力仪等效高度的方法中,通过由载波和经电光调制器调制后的产生的±1阶边带组成的拉曼激光对驱动原子双光子受激拉曼跃迁,在原子不同的自由下落时刻扫描拉曼光的脉冲持续时间,测量记录每个拉曼脉冲持续时间和每个时刻对应的原子跃迁概率;在每个时刻拟合拉比振荡概率曲线,得到每个时刻对应的有效拉比频率,再对每个时刻对应的有效拉比频率进行拟合,得到原子运动坐标系水平面和仪器坐标系水平面之间的高度差;根据计算得到原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标、原子运动坐标系水平面和仪器坐标系水平面之间的高度差,以及预设的参考水平面与仪器坐标系的高度,得到原子重力仪相对于参考水平面的等效点高度。本发明通过调制拉曼激光器,与受激拉曼跃迁相关的有效拉比频率不再是常数,而是垂直方向上的周期性余弦函数,使得可以简单而精确地计算原子重力仪相对于特定参考水平的等效高度。
应该理解的是,虽然图1的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图1中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (9)
1.一种通过调制拉曼激光确定原子重力仪等效高度的方法,其特征在于,所述方法包括:
在拉曼脉冲原子干涉重力仪中,以MOT中心为原点建立原子运动坐标系,同时以反射镜上表面为XY平面,Z轴与所述原子运动坐标系一致,建立仪器坐标系;
通过由载波和经电光调制器调制后的产生的±1阶边带组成的拉曼激光对驱动原子双光子受激拉曼跃迁,在原子不同的自由下落时刻扫描拉曼光的脉冲持续时间,测量记录每个拉曼脉冲持续时间和每个时刻对应的原子跃迁概率;其中,每个原子下落时刻扫描一系列拉曼光脉冲持续时间所获得的原子跃迁概率对应该时刻原子拉比振荡曲线;
根据预设的第一公式,在每个时刻拟合拉比振荡概率曲线,得到每个时刻对应的有效拉比频率;
根据预设的第二公式,对所述每个时刻对应的有效拉比频率进行拟合,得到所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差;
根据计算得到原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标、所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差,以及预设的参考水平面与仪器坐标系的高度,得到原子重力仪相对于所述参考水平面的等效点高度。
2.根据权利要求1所述的方法,其特征在于,所述第一公式为:
其中,dam为中间变量,dam=zm,A-za,A,zm,A为所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差,za,A为ta时刻原子团在所述原子运动坐标系中的z轴坐标,τ为拉曼脉冲持续时间,P(dam,τ)为参数dam和参数τ对应的原子跃迁概率,Ω表示有效拉比频率。
3.根据权利要求2所述的方法,其特征在于,在每个时刻ta拟合拉比振荡概率曲线P(τ),得到每个时刻ta对应的有效拉比频率Ω(ta),还包括:
将所述第一公式进一步优化,得到第一公式优化后的拟合公式:
P(τ)=b1[1-exp(-b2τ)cos(Ωτ)]+b4其中,b1、b2、Ω、b4均为待拟合的参数;
通过所述第一公式优化后的拟合公式在每个时刻ta拟合拉比振荡概率曲线P(τ),得到每个时刻ta对应的有效拉比频率Ω(ta)。
4.根据权利要求1所述的方法,其特征在于,所述第二公式为:
其中,Ω表示有效拉比频率,Ω0为拉曼激光组中只有E0和E+1时的有效拉比频率,Δk=kn+1-kn为拉曼激光对的波矢差,Δ0,+1为E0与E+1组成的拉曼激光对驱动原子双光子拉曼跃迁时的单光子失谐,Δ-1,0为E-1与E0组成的拉曼激光对驱动原子双光子拉曼跃迁时的单光子失谐。
5.根据权利要求4所述的方法,其特征在于,对所述每个时刻对应的有效拉比频率进行拟合,得到所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差,包括:
将所述第二公式进一步优化,得到第二公式优化后的拟合公式:
其中,ΩA、zm,A、υ0、ΩB均为待拟合参数,g0为测得的当地绝对重力值;
根据所述第二公式优化后的拟合公式,以ta自变量,以有效拉比频率Ω(ta)为因变量进行拟合,得到所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差拟合值
6.根据权利要求1所述的方法,其特征在于,在得到所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差拟合值之后,还包括:
根据所述原子运动坐标系和所述仪器坐标系之间的高度差拟合值通过公式计算zm,A的真实值,在此需要提前预估zm,A的值以设定调节参数i,其中zm,A的预估值与真实值之间的误差不大于Td。
7.根据权利要求1所述的方法,其特征在于,根据计算得到原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标、所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差,以及预设的参考水平面与仪器坐标系的高度,得到原子重力仪相对于所述参考水平面的等效点高度,包括:
根据计算得到原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标、所述原子运动坐标系水平面和所述仪器坐标系水平面之间的高度差,以及预设的参考水平面与仪器坐标系的高度,得到原子重力仪相对于所述参考水平面的等效点高度为:
heq=zm,A-zeq,A-dmr
其中,zeq,A为计算得到原子重力仪的绝对重力测量值在原子运动坐标系中对应的等效点坐标,计算式为v1表示原子团在自由下落至t1时刻且获得第一个拉曼脉冲的反冲速度vrec转移后的速度,T表示拉曼脉冲之间的间隔时间,z1,A表示t1时刻原子团的位置,dmr为所述预设的参考水平面与仪器坐标系水平面的高度差。
8.根据权利要求1所述的方法,其特征在于,所述拉曼激光对驱动的原子为87Rb原子。
9.根据权利要求1所述的方法,其特征在于,所述拉曼脉冲中的脉冲为π/2-π-π/2序列。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111439466.8A CN114137626B (zh) | 2021-11-30 | 2021-11-30 | 一种通过调制拉曼激光确定原子重力仪等效高度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111439466.8A CN114137626B (zh) | 2021-11-30 | 2021-11-30 | 一种通过调制拉曼激光确定原子重力仪等效高度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114137626A CN114137626A (zh) | 2022-03-04 |
CN114137626B true CN114137626B (zh) | 2023-11-21 |
Family
ID=80389421
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111439466.8A Active CN114137626B (zh) | 2021-11-30 | 2021-11-30 | 一种通过调制拉曼激光确定原子重力仪等效高度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114137626B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117347928A (zh) * | 2023-12-06 | 2024-01-05 | 湖南国科通导时空科技有限公司 | 一种激光功率比参数自动标定方法、***及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013144503A1 (fr) * | 2012-03-28 | 2013-10-03 | ONERA (Office National d'Etudes et de Recherches Aérospatiales) | Appareil d ' interferometrie atomique de mesure inertielle differentielle |
CN110596785A (zh) * | 2019-10-23 | 2019-12-20 | 中国人民解放军军事科学院国防科技创新研究院 | 一种适用于原子干涉重力仪的便携式振动噪声修正补偿方法及装置 |
CN111610572A (zh) * | 2020-05-09 | 2020-09-01 | 中国人民解放军军事科学院国防科技创新研究院 | 一种原子干涉重力仪波前相移快速评估方法及*** |
CN111897024A (zh) * | 2020-06-08 | 2020-11-06 | 北京无线电计量测试研究所 | 一种冷原子重力仪及检测方法 |
CN113219545A (zh) * | 2021-04-26 | 2021-08-06 | 中国人民解放军军事科学院国防科技创新研究院 | 一种原子干涉重力仪振动噪声及波前畸变误差双补偿方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IT1403617B1 (it) * | 2010-12-29 | 2013-10-31 | Eni Spa | Dispositivo di misura gravimetrica assoluta a interferometria atomica per applicazioni geofisiche particolarmente per il monitoraggio di giacimenti di idrocarburi |
US8879059B2 (en) * | 2012-02-16 | 2014-11-04 | The Charles Stark Draper Laboratory, Inc. | Methods and apparatus for controlled generation of hyperfine polarizations and coherences |
-
2021
- 2021-11-30 CN CN202111439466.8A patent/CN114137626B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013144503A1 (fr) * | 2012-03-28 | 2013-10-03 | ONERA (Office National d'Etudes et de Recherches Aérospatiales) | Appareil d ' interferometrie atomique de mesure inertielle differentielle |
CN110596785A (zh) * | 2019-10-23 | 2019-12-20 | 中国人民解放军军事科学院国防科技创新研究院 | 一种适用于原子干涉重力仪的便携式振动噪声修正补偿方法及装置 |
CN111610572A (zh) * | 2020-05-09 | 2020-09-01 | 中国人民解放军军事科学院国防科技创新研究院 | 一种原子干涉重力仪波前相移快速评估方法及*** |
CN111897024A (zh) * | 2020-06-08 | 2020-11-06 | 北京无线电计量测试研究所 | 一种冷原子重力仪及检测方法 |
CN113219545A (zh) * | 2021-04-26 | 2021-08-06 | 中国人民解放军军事科学院国防科技创新研究院 | 一种原子干涉重力仪振动噪声及波前畸变误差双补偿方法及装置 |
Non-Patent Citations (2)
Title |
---|
原子干涉重力仪集成光源***综述;吕梦洁;王光明;颜树华;王国超;朱凌晓;王亚宁;贾爱爱;张旭;电子测量与仪器学报(007);全文 * |
基于原子干涉技术的精密重力测量研究;陈乐乐;罗覃;邓小兵;谈玉杰;毛德凯;张恒;周敏康;段小春;邵成刚;胡忠坤;;中国科学:物理学 力学 天文学(07);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114137626A (zh) | 2022-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11175139B2 (en) | Hybrid inertial measurement system and method using a light pulse cold atom interferometer | |
CN114137626B (zh) | 一种通过调制拉曼激光确定原子重力仪等效高度的方法 | |
JP2009533666A (ja) | 光ファイバ測定方法及び装置、及び電動ジャイロスコープ | |
CN103961101A (zh) | 磁共振成像设备及其控制方法 | |
CN108646203A (zh) | 一种纳米尺度的微波磁场测量方法 | |
NL9002842A (nl) | Magnetische resonantie werkwijze en inrichting ter reductie van beeldfouten in een magnetisch resonantie beeld. | |
CN109443337B (zh) | 一种基于金刚石内nv色心的定位导航***与方法 | |
CN105091844B (zh) | 一种高精度动态测角装置及方法 | |
CN114119596A (zh) | 多相位图像任意反射面速度干涉仪数据处理方法 | |
CN108759708A (zh) | 一种对尖峰噪声免疫的高精度相干峰定位方法 | |
CN114415262A (zh) | 一种基于等效零偏的重力仪测量误差补偿方法 | |
Lee et al. | Zero lock-in implementation by phase wrapping/unwrapping in a ring laser gyroscope | |
CN111964658B (zh) | 一种旋转场驱动的核磁共振陀螺闭环磁共振方法 | |
CN113639737A (zh) | 一种光纤陀螺零偏确定方法 | |
Shavshin | Development of automatic control of optical signal gain in a atomic frequency standard based on rubidium-87 atoms | |
Cheiney et al. | Demonstration of a robust hybrid classical/quantum accelerometer | |
Celikel et al. | Establishment of all digital closed-loop interferometric fiber-optic gyroscope and scale factor comparison for open-loop and all digital closed-loop configurations | |
Arinushkina et al. | Improvement of the frequency standard on cesium atoms used in spacecraft for remote sensing of the Earth | |
CN101144720A (zh) | 光纤陀螺渡越时间在线精密测量方法 | |
CN110411434A (zh) | 一种原子陀螺仪三维原位快速磁补偿方法 | |
Zhang et al. | Detection method for the singular angular velocity intervals of the interferometric fiber optic gyroscope scale factor | |
CA1329251C (en) | Optical fibre measuring device, gyrometer, central navigation and stabilizing system | |
CN116413706A (zh) | 一种移动载体上激光雷达同时建图和内参标定方法 | |
CN111947705B (zh) | 一种对磁光阱垂向位置的漂移标定及反馈锁定方法 | |
CN104296740A (zh) | 一种基于反熔丝fpga的光纤陀螺主控板晶振选取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |