基于EEMD分解的气压改正方法研究

利用超导重力仪观测数据精确测定低于1 mHz的地球自由振荡简正模式的分裂频率,是在不与任何弹性系数发生联系的情况下改善一维密度模型的有效方法.但在该频段台站局部气压变化对重力观测数据的影响成为主要干扰来源,且具有频率依赖特性,因此精细地开展气压改正成为利用超导重力数据检测低频自由振荡信号的必要手段.本文基于EEMD方法,提出了一种具有频率依赖特性的气压改正方法.该方法将重力观测和气压变化分解成处于不同频段的本征模态函数,并在相应频段上分别进行重力-气压变化的回归分析,计算得到具有频率依赖特性的气压导纳值,精细地消除气压变化对重力观测的影响,并以此对微弱低频地球自由振荡信号开展高分辨率分析.基于本文提出的气压改正方法,利用大地震后的超导重力数据检测了频率小于1.5 mHz的低频地球自由振荡及其频谱分裂现象.研究结果表明：利用该方法进行气压改正后检测得到的各简正模具有更高的信噪比,估计的本征频率误差水平明显降低,获得的基频球型振荡₀S₂和₀S₃以及一阶球型振荡₁S₂的分裂谱峰的估计精度更高,同时还检测到了部分环型振荡在重力观测中的耦合现象.对低频地球振荡的高分辨率检测结果验证了基于EEMD分解提出的气压改正方法的有效性,同时再次证明了超导重力仪观测数据在低频地球自由振荡检测中的优势.

     

基于集合经验模态分解法的重力极潮提取与研究

重力极潮,即极移引起的重力变化,是地球由于惯性离心力变化导致地球形变的综合反映,其观测和研究有助于了解地球在长周期频段的响应,约束地球形变、地球内部构造和物理参数.区别于以往研究,本文采用了集合经验模态分解方法（Ensemble Empirical Mode Decomposition,EEMD）从全球5个超导重力仪台站连续重力观测资料中提取重力极潮、消除仪器漂移及极潮频段以外的大部分噪声信号,该方法所用的基函数是基于自身信号获得的;在此基础上,利用最小二乘匹配法分离极潮中的周年项和Chandler项,估算Chandler周期的极潮潮汐因子.结果表明,基于EEMD得到的极潮潮汐因子与前人基于其他消除仪器漂移方法（数学模型或小波分析）得到的结果相符合,精度相当,但由于这种方法是自适应的,因此本文结果更能反映实际物理过程.

     

全景PET/CT的研究进展及轴向视场的新突破

正电子发射计算机断层成像技术（PET）是一种基于功能改变的显像技术,能够通过特定的示踪剂在体内的分布特点来反映体内某方面的功能代谢水平。近年来,PET诊断得到了越来越多的临床认可,其显像系统也发展迅速,多种功能PET仪器相继推出,2017年国内推出了全球首台全景动态扫描PET/CT（轴向视场为196.6 cm）。本文将对近年来PET/CT的技术发展和全景动态扫描PET/CT研究突破进行相关综述。     

利用Excel实现快速整理CG-5型重力仪静态试验数据和计算漂移常量DRIFT

重力仪静态试验是重力勘探工作开始之前对仪器性能检查的必要环节,由于原理简单,其数据整理常不被人们重视,没有系统的方法,但整理的计算过程却又繁琐复杂,初学者在面对大量数据和多重目标时或顾此失彼,或重复计算,往往要耗费较多的工作时间.Excel电子表格具有强大的数学计算、图形显示功能,且应用广泛,易于操作,可系统整理CG-5型重力仪静态试验的原始观测记录,快速获取静态试验数据表、零点位移曲线及其拟合直线图等直观要素,同时计算得到静态零点位移曲线与直线偏差、平均零点位移率等结果参数,Excel的散点图趋势线功能也改进了静态零点位移校正参数漂移常量DRIFT的计算方法与准确度.这种静态试验数据的整理方法具有快速、准确和直观的优点.     

相关性分析在GT航空重力仪零漂改正中的应用

航空重力仪测量传感器和地面重力仪测量单元类似,往往存在零点漂移现象,通常把这种漂移视为线性关系,用前、后校数据做零漂改正处理。但是航空重力仪结构复杂,造成传感器漂移的因素很多,最为关键的因素是环境温度变化带来的扰动,实际又无法避开航空重力仪在温差剧烈变化的夏、冬两季作业。笔者以航空重力仪静态数据为例,首先进行固体潮改正并使用最小二乘法做线性零漂改正,其次使用线性相关分析方法,获得了静态数据和仪器温控温度之间具有很强的相关性特征,最后建立线性回归方程计算得到最终改正后的重力静态数据,精度有了很大的提高。这种相关性分析方法对深入理解传感器零漂改正处理方法和应用于动态测量数据中的零漂改正预研究具有实际意义。     

利用超导重力数据探测Slichter模三重分裂信号

本文基于2004苏门答腊地震后全球地球动力学计划观测网中9台超导重力仪13200 h的连续重力观测数据,联合利用整体经验模态分解(EEMD)、最优序列估计(OSE)、积谱密度分析(PSA)和自回归估计(AR)方法,探测地球固态内核平动振荡模态(Slichter模)的三重分裂信号.利用EEMD提取常规预处理后的重力残差中包含目标频段的本征模态函数(IMFs)作为最终的重力残差之后,将每个台站的重力残差平均分成无重叠的3个子数据块,再应用OSE和PSA方法获取Slichter模三重分裂谱线的积谱密度,并采用AR方法估计各积谱中弱共振信号的中心频率及其误差.结果表明,结合OSE和PSA方法探测到三个清晰的谱峰,周期分别为5.8307±7.1×10-4 h、5.2161±8.1×10-4 h和4.7536±5.2×10-4 h,分别对应Slichter模三重分裂谱线m=-1、m=0和m=+1,且三个谱峰的周期与Crossley(1992)、Rochester和Peng(1993)、Peng(1997)和Rogister(2003)基于PREM地球模型给出的Slichter模理论周期非常接近,也与Ding和Shen(2013)建议的"可能的观测结果"较为接近.参考现有的理论和观测结果,本文认为这三个信号可能是Slichter模三重分裂谱线.     

基于旋转微椭地球模型的内核平动振荡三重谱线理论模拟与实验探测

本文基于旋转微椭地球模型,采用简正模理论计算了地球内核平动振荡三重谱线的本征周期,理论上系统研究了地球内部介质（包括密度、地震波速等）分布异常对三重谱线本征周期的影响,计算了不同的内外核密度差和地核中的不同的P/S波速对应的内核平动振荡理论三重谱线周期;利用全球分布的9个超导台站,迭积每个台站长达54个月的高精度超导重力仪数据,在亚潮汐频段（0.162～0.285 cph）检测内核平动振荡三重谱线.结果发现,三重谱线本征周期对内外核边界的密度跳跃非常敏感,随着密度差的增加,以类似于双曲线的特征减小;无论是采用地球质量不变的方法还是采用浮力频率为常数的方法,计算得到的三重谱线本征周期结果相差较小,且随着内外核密度差的增大,差距逐渐减小;内、外核P波波速分布异常对三重谱线周期的影响基本相当,内核S波波速分布异常比P波波速分布异常对三重谱线周期的影响小1个量级;探测到一组信噪比较高且满足谱峰分裂特征的三重谱线的信号（0.19281,0.21456和0.24151 cph）,有极大的可能是来自于内核平动振荡.基于探测结果可以推断实际的地球模型其内外核密度差应该介于PREM模型和1066A地球模型之间,更接近于1066A模型.     

激光干涉绝对重力仪参考棱镜隔振系统仿真

在激光干涉绝对重力仪中,参考棱镜在背景振动和系统自振的作用下呈现复杂的振动模式,严重影响了整机的测量准确度.针对此问题,本文基于磁悬浮原理设计了参考棱镜的反馈式两级振动隔离系统,建立了系统的动力学模型,在不同的反馈方式下推导了相应的传递函数,进行了Matlab仿真分析以确定最优的反馈方式.系统仿真结果表明绝对速度反馈方式的隔振效果可达100 dB,与自主研发的Age-110型激光干涉绝对重力仪的联合测试结果表明,该隔振系统可将Age-110的测量精度再提高1.5倍.

     

重力梯度同步观测方法

为了克服传统陆地重力梯度测量方法的不足,提出了利用两台或多台具有电子读数的相对重力仪测量重力梯度的同步观测方法和相应的数据处理模型,该方法无需进行固体潮改正和零漂改正而直接求解重力差,垂直重力梯度实验结果验证了该方法的有效性。     

东日本大地震重力信号同震响应

对2011年3月11日日本本州东海岸附近海域Mw9.0级地震前后的观测重力资料进行了同震响应信号分析,给出了gPhone重力仪在大地震后的长周期地震面波和自由振荡的观测结果。震后12h数据经带通滤波后可以分辨出7个环绕地球传播的瑞雷波记录,在低频范围内(5mHz以下)检测到43个地球自由振荡基频振型及52个谐频振型,其分布与理论模型吻合较好。其中,在1mHz以下的低频段内观测到4个振型的谱峰分裂现象。