地球物理学报  2018, Vol. 61 Issue (7): 3037-3043   PDF    
宽频带超导重力仪的研制进展
刘习凯, 马东, 陈亮, 刘向东     
华中科技大学物理学院基本物理量测量教育部重点实验室, 武汉 430074
摘要:本文报道一种构建宽频带高分辨率超导重力仪(SG)的方法.与已商业化的GWR超导重力仪类似,新型超导重力仪同样使用磁悬浮超导检验质量来构建弹簧振子结构,不同的是弹簧振子的固有频率更高,可以有效覆盖地脉动频带.同时,利用基于超导量子干涉器件(SQUID)的位移传感技术,大幅提高了位移传感灵敏度,以保证仪器在测量带宽扩展之后仍具有与GWR仪器相当的测量分辨率,从而实现对时变重力信号与地球背景噪声信号的同时高分辨观测.宽频带超导重力仪核心部件垂向超导加速度计(VSA)的设计、组装和初步测试已经完成.结果表明,加速度计的噪声本底为8 μGal/Hz1/2@0.1 Hz,实现了对地脉动信号的高分辨率测量.而为了将仪器噪声本底降低到GWR仪器水平,加速度计的固有频率需要进一步减小,该部分参数优化工作正在进行.本文将对垂向超导加速度计的工作原理、结构和测试结果进行详细介绍.
关键词: 宽频带      垂向超导加速度计      分析      测试     
Progress in the development of a wide-band superconducting gravimeter
LIU XiKai, MA Dong, CHEN Liang, LIU XiangDong     
MOE Key Laboratory of Fundamental Physical Quantities Measurement, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China
Abstract: In this paper, a method to construct superconducting gravimeter (SG) has been described. Similar to the commercial GWR superconducting gravimeter (SG), a superconducting test mass is levitated to form the spring-oscillator. However, the spring-oscillator is designed to have higher natural frequency, so that the frequency band can cover the microseism signals. At the same time, more sensitive technology based on superconducting quantum interference devices (SQUID) is used to detect the displacement of the spring-oscillator. In this way, high sensitivity can be attained even if the instrument bandwidth has been broadened. Therefore, the constructed SG can record both the time-varying gravity signals and the earth ambient noise simultaneously. A prototype of vertical superconducting accelerometer (VSA) has been fabricated with a nature frequency of about 20 Hz. The accelerometer shows a noise floor of 8 μGal/Hz1/2@ 0.1 Hz, and the microseism signals have been measured by the VSA with high sensitivity. In order to improve the instrumental noise level to be better than the GWR SG, the natural frequency of the VSA must be further reduced. The works to optimize the accelerometer parameters is on the way. The working principle, structure and test results of the superconducting accelerometer will be present in detail.
Key words: Wide band    Vertical superconducting accelerometer    Theoretical analysis    Test    
0 引言

地球时变重力场包含地球系统物质运动的丰富信息,反映了地球各圈层基本的变化特性,因而时变重力场的观测是地球物理学与地球动力学研究的重要手段(孙和平等,2006).长期以来,超导重力仪以其低噪声(0.1~0.3 μGal/Hz1/2@1~20 mHz)和高稳定性(μGal/年水平),被公认为是地面检测微小重力变化的最佳仪器(Prothero and Goodking, 1968; Goodking,1999),在地球时变重力场的测量中扮演着重要角色(Hinderer et al., 2009Sun et al., 2002).目前,超导重力仪由美国GWR仪器公司独家生产.GWR型超导重力仪利用超导载流线圈将超导检验质量悬浮起来构成弹簧振子,使用电容位移传感技术对弹簧振子位移进行检测(Goodking,1999).为了达到高分辨测量的目的,弹簧振子结构的固有频率被设置在0.1 Hz,因而仪器的带宽也被限制在0.1 Hz以下.

近年来,地球背景噪声由于包含丰富的地下结构信息,被人们广泛关注和研究(Ide et al., 2011Sun et al., 2002).其中1 Hz以上的成分被认为是人类活动产生,而0.005~0.3 Hz的噪声主要归于海浪与大陆架的相互作用、海底洋流作用、地震等,具有较高的研究价值.GWR型超导重力仪无法对0.1 Hz以上的地球背景噪声进行有效测量.如果研制一种测量分辨率与GWR仪器相当、仪器带宽达到1 Hz的超导重力仪,就能利用同一台仪器完成对时变重力信号和地球背景噪声的同时测量.如果使用这样的超导重力仪组网观测(Imanish et al., 2004)地球低频背景噪声(Xu et al., 2004),可以获得丰富的地下结构特征(孙和平等,2004),因而在地球壳幔结构的研究、资源勘探技术和地震预测方法研究中具有广泛的应用前景(陈运泰,2009刘子维等,2011赵根模等,2001Kerr,1998Shen et al., 2011).

为此,我们提出一种构建宽频带高分辨超导重力仪的方法,其要点是采用具有更高灵敏度的基于SQUID的位移探测方法(Chan and Paik, 1987)替代电容位移传感,使仪器在固有频率增大,加速度到位移的传递函数减小的情况下,仍然具有满足需求的高分辨率.

我们已经完成宽频带超导重力仪中的核心单元——垂向超导加速度计的设计和制作,并进行了初步的测试,在实验室环境下,噪声水平为8 μGal/Hz1/2.本文将对该加速度计的工作原理、结构和测试结果进行介绍.

1 垂向超导加速度计动力学与噪声分析

垂向超导加速度计是超导重力仪的核心部件,本节首先对垂向超导加速度计进行了完整的理论分析.结合具体的超导电路,从检验质量动力学方程出发,给出加速度计的固有频率和传递函数,进而计算加速度计的固有噪声水平.

1.1 垂向超导加速度计工作原理

垂向超导加速度计的工作原理如图 1所示.在超导载流线圈所在回路中注入持久超导电流后,超导检验质量表面会形成具有确定分布的超导屏蔽电流,从而给检验质量一个方向向上的磁斥力,用来平衡检验质量重力,实现检验质量的磁悬浮,构成磁力弹簧振子结构.当重力加速度发生变化时,悬浮的检验质量与固定的超导线圈之间位置关系发生相应的变化,根据Messina效应,此时检验质量的超导屏蔽电流会随之改变,屏蔽电流被超导线圈探测的磁通亦随之改变,这样,检验质量的位移会被转化成超导探测线圈的磁通变化量.由于超导回路具有类磁通守恒性质,超导回路中电流会自发改变,将检验质量的位移进一步转化为超导回路中的电流变化,然后利用SQUID对回路中的变化电流进行高灵敏度检测,最终以电压形式进行信号输出.

图 1 垂向超导加速度计原理图 Fig. 1 Principle of the vertical superconducting accelerometer
1.2 垂向超导加速度计理论分析

超导重力仪器的理论分析方法已经由马里兰大学Paik和Chan等给出(Chan et al., 1987),我们结合自主设计的超导电路(如图 1所示),对超导加速度计进行了详细的分析.

以检验质量偏离平衡位置向上为正方向,当检验质量位移为x时,线圈的电感变化可以用(1)式表示(Chan et al., 1987):

(1)

其中L10是检验质量处于平衡位置时的有效电感值(有效电感指超导线圈在自身通电电流产生的磁通与探测到的超导屏蔽电流产生的磁通共同作用下的电感值),λ是线圈相关参数,可以通过实验测量得到.

图 1中的超导电路结构含有2个超导回路,在与检验质量相互作用的超导悬浮回路中通入恒定电流I0后,两个回路的磁通必须守恒(Chan et al., 1987):

(2)

其中I为超导悬浮回路变化电流,i为通过互感耦合至SQUID输入线圈的电流,为有效互感.联立式(1)和(2)可以解出超导回路变化电流对检验质量位移的依赖关系:

(3)

对(3)式进行Taylor展开得到位移到电流的传递函数:

(4)

外界加速度所引起的检验质量位移变化由弹簧振子的固有频率决定,在超导加速度计中,固有频率由超导载流线圈所产生的电磁作用力决定.回路的电流变化已由式(3)给出,那么系统的电磁能为:

(5)

电磁能对位移的微分给出超导线圈电磁作用力:

(6)

式(6)中右边第一项为静态作用力,用来抵消检验质量的重力,实现检验质量的磁悬浮:

(7)

其中m为检验质量的质量,g0为当地重力加速度.式(6)右边的第二项是与检验质量位移成比例、方向相反的动态作用力,具有弹簧振子回复力的性质.

检验质量的动力学方程可以表示如下:

(8)

其中γ为振子阻尼系数,在高真空环境中,加速度计的Q值很高,阻尼的影响较小.Δg为时变重力信号,aP为仪器工作平台的加速度.加速度计不能将平台的运动加速度与时变重力信号区分开来,当平台运动噪声过大时,会影响重力加速度的测量分辨率.

(4) 和(8)两式可以给出加速度到探测电流的传递函数:

(9)

其中ω0为电磁弹簧振子的固有频率,

(10)

(10)式表明,当悬浮线圈有效电感L10和耦合器互感M确定时,通过改变L2线圈电感值,可以对固有频率进行调节.

1.3 仪器固有噪声

垂向超导加速度计的固有噪声包括弹簧振子热噪声以及SQUID探测系统噪声(Chan et al., 1987).

在工作温度为T时,弹簧振子的热噪声功率谱密度可由(11)式表示:

(11)

其中kB为玻尔兹曼常数;T为仪器工作温度;m为检验质量的质量;γ为弹簧振子的阻尼系数.

SQUID可以看作是电流到电压的转换放大器,SQUID所引入的固有噪声由加速度计的传递函数和SQUID的能量阈值共同决定.我们使用的SQUID是从美国Star Cryoelectronics公司购置的直流偏置SQUID,其能量阈值可由式(12)近似表示(Moody,2011):

(12)

是由白噪声和低频噪声组成.

SQUID探测系统的固有噪声通过(13)式转换为加速度测量噪声:

(13)

式中Li为SQUID输入线圈电感;Hgi为加速度到SQUID输入电流的传递函数,已经在1.2节进行了详细计算.

仪器的固有噪声是弹簧振子热噪声和SQUID放大器噪声之和:

(14)

根据表 1所设计的超导电路参数,弹簧振子固有频率为17.4 Hz,低频时加速度到SQUID输入线圈电流传递函数为0.01 A/m/s2.垂向超导加速度计固有噪声曲线如图 2所示.

表 1 超导电路参数表 Table 1 Parameters of the superconducting circuit
图 2 垂向超导加速度计固有噪声谱密度图 Fig. 2 Intrinsic noise of the vertical superconducting accelerometer

当固有频率为17.4 Hz时,在低频处,加速度计的固有噪声主要受限于SQUID放大器的噪声.在1 mHz处,固有噪声为0.1 μGal/Hz1/2,在1 Hz处接近0.01 μGal/Hz1/2.式(10)和(15)表明,当固有频率减小时,低频固有噪声会进一步降低,分辨率会进一步提高.当固有频率为1 Hz时,在1 Hz以下频段,加速度计固有噪声会进一步降低1个量级(参见图 2).因此,宽频带的高分辨超导加速度计方案具有理论可行性.

2 垂向超导加速度计的结构设计

加速度计的结构主要包含三部分:检验质量、超导线圈陶瓷骨架、探头封装结构.其整体示意图如图 3所示.

图 3 加速度计结构示意图 Fig. 3 Structure schematic of the accelerometer

检验质量采用中间有隔板的圆筒,质量比GWR的空心超导铌球重1~2个量级,为0.1 kg,大检验质量的使用有助于降低系统的机械热噪声.同时在其侧壁柱面安装超导线圈,对检验质量的水平运动进行限制,减小地面水平振动耦合干扰.

超导线圈陶瓷骨架有两种类型.一种是与检验质量中间隔板相对,用来安装平面超导线圈,中间隔板的两侧分别安装1个单层密绕平面超导线圈,下部超导线圈对检验质量进行悬浮和位移探测,上部超导线圈进行辅助探测;另一种是环绕检验质量的陶瓷骨架,用来安装侧壁线圈,在检验质量侧壁均匀地安装有8组柱面线圈,其安装示意图如图 4所示.正对的两个线圈为一组,两组线圈相互正交对称安装,实现对检验质量两个水平方向自由度运动的限制.

图 4 超导加速度计中线圈安装示意图 Fig. 4 Installation of the superconducting coils

探头封装结构实现对内部陶瓷骨架的固定以及对超导线圈等元件的电磁屏蔽.组装完成的垂向超导加速度计如图 5所示.

图 5 组装好的超导加速度计 Fig. 5 Assembled accelerometer
3 初步测试结果与分析 3.1 实验环境

仪器测试在实验室液氦实验杜瓦中进行.实验过程中,采用高精度CMG-3T型微震仪对地面震动进行实时监控,同时对加速度计的分辨率进行标定.

3.2 加速度计测试

图 6为检验质量未起浮时SQUID输出信号,直接反映了仪器的本底噪声.噪声主要来源于电磁干扰,高频电磁信号会经由电流注入线等直接传输至超导回路中.

图 6 检验质量未起浮,加速度计本底信号 Fig. 6 Background noise of the accelerometer

在将6.8 A悬浮电流注入到超导电路后,检验质量进入磁悬浮状态,仪器的输出信号如图 7所示.

图 7 检验质量起浮后自由振动信号 Fig. 7 Oscillation signal of the levitated test mass

我们利用该加速度计监测实验室地面震动,并与CMG-3T型微震仪测量数据进行对比,结果如图 8所示.加速度计已经实现了对地脉动信号的高分辨率测量,噪声水平为8 μGal/Hz1/2@0.1 Hz.在0.1 Hz以上频段加速度计的测量明显受限于地面震动.

图 8 垂向超导加速度计噪声谱密度(红色)与微震仪测量(蓝色)对比 Fig. 8 PSD (Power Spectral Density) of the accelerometer (red) and the seismometer (blue)

加速度计的固有噪声水平远小于目前的测量水平,说明现阶段加速度计噪声主要受电磁干扰、温度波动等外部噪声的限制.外部噪声对应的仪器输出电压一定时,固有频率越小, 加速度计的加速度到位移的传递函数越大,等效的加速度噪声越小.当固有频率为17.4 Hz时,加速度噪声为8 μGal/Hz1/2@0.1Hz,如果将加速度计的固有频率降低到1 Hz附近时,外部噪声所等效的加速度噪声将降低1/300,至10-2μGal/Hz1/2的水平.

4 结论

本文提出了一种新型宽频带超导重力仪的构建方案,并对其核心部件垂向超导加速度计进行了详细的分析和介绍.自主研制的固有频率为17.4 Hz的垂向超导加速度计实现了检验质量的稳定磁悬浮,噪声本底为8 μGal/Hz1/2@0.1 Hz,主要受限于温度波动等外部噪声.

为了突破外部噪声对垂向超导加速度计噪声水平的限制,仪器的固有频率需进一步降低.如果固有频率可以降低到1 Hz,那么外部噪声所等效的加速度噪声将减小至10-2μGal/Hz1/2的水平,仪器的噪声水平将会显著降低.固有频率的调节优化工作已经在进行之中.

致谢

感谢课题组中方晖、张燕在加速度计结构设计、超导线圈绕制方面的贡献,同时感谢胡远旺、胡咸军、肖再林、陶文兵、任小龙、赵坤雷在长期的低温实验中给予的帮助.

References
Chan H A, Paik H J. 1987. Superconducting gravity gradiometer for sensitive gravity measurements. Ⅰ. Theory. Physical Review D, 35(12): 3551-3571. DOI:10.1103/PhysRevD.35.3551
Chan H A, Moody M V, Paik H J. 1987. Superconducting gravity gradiometer for sensitive gravity measurements. Ⅱ. Experiment. Physical Review D, 35(12): 3572-3597. DOI:10.1103/PhysRevD.35.3572
Chen Y T. 2009. Earthquake prediction:Review and prospect. Science China Series D:Earth Sciences (in Chinese), 39(12): 1633-1658.
Goodking J M. 1999. The superconducting gravimeter. Review of Scientific Instruments, 70(11): 4131-4152. DOI:10.1063/1.1150092
Hinderer J, Crossley D, Warburton R J. 2009. Gravimetric Methods Superconducting gravity meters. //Schubert G ed. Treatise on Geophysics. Amsterdam: Elsevier.
Ide S, Baltay A, Beroza G C. 2011. Shallow dynamic overshoot and energetic deep rupture in the 2011 MW9.0 Tohuko-Oki earthquake. Science, 32(6036): 1426-1429.
Imanish Y, Sato T, Higashi T, et al. 2004. A network of superconducting gravimeters detects submicrogal coseismic gravity changes. Science, 306(5695): 476-478. DOI:10.1126/science.1101875
Kerr R A. 1998. Seismology:A slow start for earthquakes. Science, 279(5353): 985. DOI:10.1126/science.279.5353.985
Liu Z W, Wei J, Hao H T, et al. 2011. Pre-seismic gravity disturbance of high frequency before MS9.0 earthquake in Japan. Journal of Geodesy and Geodynamics (in Chinese), 32(2): 4-8.
Moody M V. 2011. A superconducting gravity gradiometer for measurements from a moving vehicle. Review of Scientific Instruments, 82(9): 094501. DOI:10.1063/1.3632114
Prothero W A, Goodking J M. 1968. A superconducting gravimeter. Review of Scientific Instruments, 39(9): 1257-1262. DOI:10.1063/1.1683645
Shen W B, Wang D J, Hwang C. 2011. Anomalous signals prior to Wenchuan earthquake detected by superconducting gravimeter and broadband seismometers records. Journal of Earth Science, 22(5): 640-651. DOI:10.1007/s12583-011-0215-4
Sun H P, Hu H T, Jentzsch G, et al. 2002. Tidal gravity observations obtained with a superconducting gravimeter at Wuhan/China and its application to geodynamics. Journal of Geodynamics, 33(1): 187-198.
Sun H P, Xu J Q, Ducarme B. 2004. Detection of the translational oscillations of the earth's solid inner core based on the international superconducting gravimeter observations. Chinese Science Bulletin, 49(11): 1165-1176. DOI:10.1360/03wd0242
Sun H P, Xu J Q, Li Q. 2006. Detail spectral structure of earth's gravity field and its application. Progress in Geophysics (in Chinese), 21(2): 345-352.
Xu J Q, Sun H P, Ducarme B. 2004. A global experimental model for gravity tides of the earth. Journal of Geodynamics, 38(3-5): 293-306. DOI:10.1016/j.jog.2004.07.003
Zhao G M, Yang G S, Chen H R. 2001. Quiescent foreshock and earthquake prediction. Earthquake (in Chinese), 21(1): 69-77.
陈运泰. 2009. 地震预测:回顾与展望. 中国科学D辑:地球科学, 39(12): 1633-1658.
刘子维, 韦进, 郝洪涛, 等. 2011. 日本MS9.0地震前的重力高频扰动. 大地测量与地球动力学, 32(2): 4-8.
孙和平, 徐建桥, DucarmeB. 2004. 基于国际超导重力仪观测资料检测地球固态内核的平动振荡. 科学通报, 49(8): 803-813.
孙和平, 徐建桥, 黎琼. 2006. 地球重力场的精细频谱结构及其应用. 地球物理学进展, 21(2): 345-352.
赵根模, 杨港生, 陈化然. 2001. 寂静的前震与地震预测. 地震, 21(1): 69-77.