在综合分析国内外地形变前兆研究现状的基础上,结合汶川8.0级地震前GPS、应变、重力等研究结果,提出临震前震源区域可能出现变形不动点现象。地壳变形过程中的不动点现象是临震前区域应力场进入临界状态的重要标志,是地震潜在危险区域划定的重要依据。借助不动点理论,给出汶川地震前变形不动点集合的演化,结合大区域地震活动空区与远场显著形变异常分布,探讨变形不动点现象发生的动力学背景,验证了汶川地震前地壳变形的不动点现象是大区域应力场有序运动的结果。
基于19个连续GNSS观测站的多种采样率(10 Hz、1 Hz和30 s)数据,分析2017-09-03朝鲜核爆对中国东北地区地壳形变场的影响。静态位移分析结果显示,距核爆现场最近的GNSS观测站JLCB没有产生明显的同震形变,与基于均匀半无限弹性空间断裂位错模型正演模拟的结果一致;GNSS应变结果显示,此次核爆在东北区域引起主应变的略微调整,但总体调整量在10-8以内,没有引起区域范围内面应变的变化;高频GNSS动态形变分析结果显示,核爆在JLCB站处没有引起明显的波形变化。另外,高频GNSS受到多路径效应的影响,可能会引起cm级的扰动,利用恒星日滤波方法可以抑制多路径的影响,得到精度更高、更可靠的动态形变结果。
基于1999~2007期、2009~2013期及2013~2016期GPS速度场数据,利用最小二乘配置方法解算南北地震带附近区域应变率参数,对比应变率变化并分析变化原因。结合研究时段内MW6.0以上地震的震中分布,分析应变率参数与地震分布的相关关系;再结合地震地质,分析研究区内主要断层的地震危险性。结果显示,青藏高原近期存在加速扩展特征,其东北缘存在明显的介质性质横向不均匀特征,MW6.0以上地震与应变率的大小存在相关关系;研究区内的地震危险区分别为祁连山中段危险区、玛沁-玛曲危险区、安宁河-小江断裂系危险区及红河断裂带中南段危险区。
基于2018-01以来云南地区GNSS连续观测时间序列,采用克里金插值方法对形变场进行网格化、计算应变,并获取面应变,结合基线数据动态分析两次通海5.0级地震前区域应变场特征和震后云南地区的应变变化情况。结果表明,通海地震前,云南整体区域应变积累存在显著增强的背景,并在震前短期内出现快速反向调整;通海地震后,短期内以腾冲为代表的小滇西地区、大寨-会泽等滇东北地区挤压应变积累出现明显的减弱趋势,滇南地区的勐海-勐腊区域存在压性增强现象。
通过完备总体经验模态分解(ICEEMDAN)与最小二乘法开展GNSS垂向序列的噪声与非构造形变去除研究。首先,采用ICEEMDAN方法对GPS垂向序列进行分解,分析GNSS垂向序列的年周期项、半年周期项和趋势项;然后,通过改进的小波阈值和排列熵值法对GNSS垂向序列进行降噪处理;最后,利用最小二乘法拟合非构造形变序列,找出与其最为接近的分解层数并进行改正,获得去除噪声和非构造形变的GNSS垂向形变信息。结果表明,利用ICEEMDAN方法可实现GNSS垂向序列噪声的有效去除;采用ICEEMDAN与最小二乘法可提高GNSS垂向序列形变测量精度,其中RMS值平均减小1.28 mm,WRMS差值平均增大0.33 mm。
将一种基于数据驱动的RegEM算法引入GPS坐标时间序列插值中,分别采用模拟不同比例连续缺失数据与实测含缺失数据,比较RegEM与拉格朗日方法、三次样条方法、正交多项式方法的插值效果与性能。结果表明,对于模拟不同比例连续缺失数据,RegEM算法插值效果均优于传统方法,且在大量数据连续缺失的情况下效果最优;对于实测含缺失数据,RegEM方法插值所得序列保留方差最大化效果最好,约为正交多项式方法的1.17倍、三次样条方法的1.38倍。
在双天线GNSS/MEMS组合测姿中,低精度、低成本的MEMS加速度计输出通常包含比较大的噪声,导致由加速度计调平计算得到的水平角观测值精度较低,影响组合测姿的精度和可靠性。不同于传统的时域降噪方法,从频域的角度出发,根据实际车载运动中加速度计输出具有因汽车固有结构产生的高频振动噪声的特点,利用滑动窗口的Butterworth低通滤波器平滑加速度计输出。动态车载实验结果表明,该滤波器降噪效果明显。与滤波前相比,水平角观测值精度提升了约4倍,达到(0.643°,0.546°)。在此基础上,滤波后组合测姿精度也从(0.124°,0.738°,0.532°)提升到(0.122°,0.074°,0.052°),其中水平角精度提升了大约一个数量级,滤波后的测姿结果更加准确可靠。
对单/双频GPS/BDS/Galileo混合双差相对定位模型算法进行研究,推导了GPS/BDS/Galileo混合双差相对定位模型。选择短基线数据进行实验,对GPS、GPS/BDS标准双差、GPS/BDS混合双差、GPS/Galileo混合双差以及GPS/BDS/Galileo混合双差解算模式在定位精度、模糊度固定速度方面进行对比分析。结果表明,相对于GPS单系统或GPS/BDS双系统标准双差解算模式,双系统或三系统混合双差解算模式能有效提高模糊度首次固定时间,定位精度也有一定程度的提升。
通过分析由ERA-Interim气象再分析资料积分方法得到的天顶对流层总延迟随高程变化的规律,提出一种基于垂直剖面函数的天顶对流层延迟(ZTD)插值算法。该算法以ZTD的垂直分布规律为基础,通过垂直剖面函数实现ZTD在高程方向上的精准投影延拓,可以避免因高差较大造成的空间内插结构畸形。采用IGS站提供的高精度对流层产品进行实验验证表明,该算法相对于传统算法能够有效提高ZTD改正值的精度,尤其在高差超过1 km的情况下,相对于反距离加权法精度提升了96%,相对于空间回归法精度提升了79%。
利用香港卫星定位参考站网GNSS观测数据,提取强热带风暴“塔拉斯”与热带风暴“洛克”影响期间各测站天顶方向对流层延迟,反演香港区域大气可降水量;根据香港区域49个天文台气象站提供的实测降雨量数据,分析大气可降水量与实际降雨量的相关性,以及两次台风对香港区域水汽时空分布的不同影响。结果表明,大气可降水量在台风影响前期均上升,在大量降雨后回落,但在连续台风的间歇期间,仍高于台风来临前的水平;水汽累积是大量降雨的前提条件,当水汽累积量相近时,水汽累积时长与累积降雨量呈正相关;台风期间大气可降水量值超过65 mm的区域面积与台风等级相关,台风路径对局部水汽分布有一定的影响。
首先介绍小波变换的基本原理,分析影响小波多分辨率分解效果的因素,提出选择coif5小波基函数、交叉验证确定分解阶次的小波多尺度分解方法;然后在模型实验中利用该方法对模型局部和区域异常实现有效分离,并结合径向对数功率谱方法对分解重构异常进行定性、半定量分析;最后对川滇地区201309~201404期流动重力差分动态变化异常进行多尺度分解。在分解得到的D4~D6异常、重构的局部和区域异常中,2014年鲁甸MS6.5和景谷MS6.6地震主震位置均存在明显的异常变化特征。结合地质构造背景信息对该异常特征进行分析解释。
基于动力学方法比较分析了双精度与四精度模式下重力场模型的解算精度,主要包括缔合勒让德函数计算、数值积分器及重力场反演结果。结果显示,在勒让德函数计算方面,部分角度在双精度模式下计算至1 900阶以后会出现溢出问题,而在四精度模式下任何角度都满足精度要求,并且计算结果比双精度模式高8个量级。数值积分器Adams预测校正法积分1 d的位置和速度误差,在四精度模式下比在双精度模式下高4个量级。在精密轨道反演重力场计算方面,动力学方法在双精度及四精度模式下反演结果一致,统计其计算至60阶的累计大地水准面误差为1.29×10-5 m,这是因为动力学方法的线性误差相对计算误差而言是主要误差;非线性动力学方法在四精度模式下比在双精度模式下高7个量级,其大地水准面误差分别为8.92×10-15 m和8.16×10-8 m。
在江汉盆地的同一场地分别采用三重圆形、T形及L形3种不同台阵形式进行微动探测观测对比实验,结果表明,3种台阵形式提取到的频散曲线基本一致,但频散图谱能量集中度具有明显的差异。三重圆形台阵频散图谱能量集中度最高,L形、T形台阵在2.0~4.6 Hz频段范围内集中度较高,而在1.2~2.0 Hz频段范围内较分散;相比T形和L形台阵,三重圆形台阵更易于探测到地层分界面附近的波速突变特征,即三重圆形台阵比其他两种台阵形式准确性及可靠性更高。