为从现今庞大的GNSS观测网络中快速检测断层蠕滑形变,精细反演蠕滑时空分布及演变特征,提出集GNSS网络滤波、地表形变信息提取和地下断层蠕滑时空分布反演三者于一体的方法。该方法同时采用整个GNSS网络的时空观测阵列,利用断层形变高空间相关的特点,对覆盖断裂带地表的GNSS位移时空序列进行主成分分析,利用主成分信息快速检测并反演蠕滑断层的时空分布与演变过程。以2005年苏门答腊MW8.6地震震后余滑和2006年墨西哥Guerrero州慢滑移为例,本文方法可成功检测并反演蠕滑断层的时空分布及演变特征,结果与相关研究成果吻合。
借助Sentinel-1A卫星升降轨数据,利用合成孔径雷达差分干涉测量(DInSAR)和Okada均匀滑动模型反演获得四川荣县地区三维同震形变场及荣县地震震源参数。结果表明,此次地震造成的断裂现象大体沿NS向,两侧沉降区域中间伴随抬升现象,最大垂直形变量约18 mm;荣县地震发震断层符合逆断层特性,断层长约9.5 km,宽约2.1 km。采用时间序列分析技术获取该地区2018-11~2019-03各时间段的累计形变,结果发现,该区域在震前形变波动较小,地震是造成累计形变较大的主要原因,震后地表变化趋于稳定。
运用专业地震解释及构造建模软件,综合利用焦作地区地震数据、钻孔数据、地质信息对活断层二维浅层地震剖面进行层位与断层解释,同时对断层平面组合与空间展布规律进行研究,实现地质资料与地震数据的有机结合,构建近地表Q3地层底界面、Q地层底界面、N地层底界面分布模型和第四纪断层模型。以第四纪断层模型为基本骨架,以近地表地层分布为主要分层,模拟Q-Q3地质体、N-Q地质体的构造形态和断层发育情况,系统反映近地表N地层底界面至Q3地层底界面的空间结构特征。
基于首都圈地震台网观测数据,利用Hypo2000方法对2016-09-10河北唐山M4.1地震主震震中进行准确测定。采用近震波形反演方法,确定主震的最佳双力偶震源机制解为:节面Ⅰ走向122°、倾角60°、滑动角-42°;节面Ⅱ走向236°、倾角54°、滑动角-142°;矩震级约为MW4.3。结合近震深度震相sPL确定主震的震源深度约为6 km。使用双差定位法对余震进行重新精确定位,结果显示,余震深度集中分布在4~9 km,整体形态近于铅直,表明发震断层面具有倾向北西、倾角较陡的特点,与节面Ⅱ的性质较为一致,据此推测节面Ⅱ为发震断层面,节面Ⅰ为辅助面。将发震断层面参数与震源区附近断层性质进行对比分析认为,发震断层应为唐山断裂带中唐山-丰南断裂。
基于2008年汶川MW7.9地震精细的同震滑动分布模型,计算此次地震引起的周边地区库仑应力变化;利用变化系数(CV)评价库仑应力变化的可靠性,并评估接收断层走向、倾角、滑动角、摩擦系数、Skempton系数对库仑应力变化可靠性的影响;根据CV可靠阈值剔除不可靠库仑应力变化,分析可靠库仑应力变化与余震分布之间的关系。结果表明,库仑应力变化在正值区或负值区内部的可靠性较高,在接近正负值交界处库仑应力变化的可靠性急剧下降;接收断层倾角对库仑应力变化可靠性的影响最大;79.3%的余震位于可靠库仑应力变化分布区域,剔除不可靠库仑应力变化及不可靠库仑应力变化分布区域内余震后,余震触发率下降3.1%。
采用甘肃省CORS网和中国大陆构造环境监测网络中共48个台站的GPS观测数据,解算得到观测台站的垂直位移,并与GRACE时变重力场Mascon模型解CSR RL05M数据计算得到的垂直形变进行比较,分析区域地表垂直形变特征。结果表明,研究区内台站垂直形变存在局部特征,甘肃庆阳和平凉地区垂直形变与其他地区存在明显差异,相关系数、均方根减少量和周年信号减少量均高于其他地区;扣除趋势项后,观测台站GPS垂直位移与GRACE垂直形变时间序列相关系数均值为0.72,GPS和GRACE周年信号振幅均值分别为6.00 mm和3.70 mm,周年信号减少量和均方根误差减少量均值分别为0.51和0.29;研究区内GPS垂直位移和GRACE垂直形变时间序列一致性较强,GRACE垂直形变能有效解释50%以上的GPS垂直位移周年信号,GPS垂直位移时间序列包含的非构造形变中平均约29%来源于环境负载变化所引起的负荷形变。
基于云南丽江台超导重力仪OSG-066超过5 a的连续重力潮汐观测资料,获取高精度潮汐参数并构建地方重力潮汐模型,为丽江地区精确重力潮汐改正提供参考。研究9个全球海潮模型在丽江台的重力影响改正有效性,用时域内回归分析法和频域内小波分析法分别研究大气重力导纳值的变化特征,并利用观测潮汐数据估算自由核章动本征参数。结果表明,9个海潮模型在丽江台的重力影响改正效果基本一致,NAO.99b模型的改正精度较高;大气重力导纳值在时域和频域内均显示出较为明显的季节性变化特征,且与降雨量变化存在负相关关系。自由核章动本征周期为431.5(411.2, 454.0)恒星日,品质因子为 -2 796。
针对经验模态分解(empirical mode decomposition, EMD)降噪过程中采用相关系数(p)准则确定分界本征模态函数(intrinsic mode function, IMF)分量K值存在不准确性的问题,选用复合评价指标(T),综合考虑曲线平滑度(r)与均方根误差(RMSE)2个指标值,提出一种改进的EMD降噪方法。利用9个模拟数据和2个陆态网基准站的实测GPS高程时间序列数据进行验证,结果表明,复合评价指标比单一的相关系数指标所确定的分界IMF分量K值更准确,能够更可靠地识别噪声与信号的分界点,使降噪效果更佳。
针对目前常用的GPS实时钟差确定算法主要基于实时观测数据流实现,使产品可靠性和完整性受限于网络质量,及广播星历和IGU-P精度较低的问题,提出并实现基于小时观测文件拼接的实时钟差确定算法。该算法基于小时观测文件数据,通过将钟差估计与自适应超短期钟差预报相结合以确定实时钟差。30 d实时在线结果显示,基于小时观测文件确定的实时钟差精度约为0.25 ns,优于广播星历和IGU-P,与IGS RTS提供的实时钟差精度相当,且能保证产品的可靠性和完整性。
为提高Bernese GNSS software数据处理效率,将英特尔数学核心函数库(math kernel library,MKL)应用于Bernese 精密定轨数据处理,对比分析多个MKL矩阵求逆函数与Bernese 原有程序的计算效率。使用2019-03全球200个测站北斗/GNSS数据进行实验分析,结果表明,采用参数预消除策略时,参数预消除步骤消耗时间明显大于矩阵求逆,使用MKL处理数据效率提升不明显;而未采用参数预消除策略时,使用MKL矩阵求逆函数可显著提高矩阵求逆效率,其中dpotri函数矩阵求逆计算效率最高,消耗时间平均值为133 s,相比Bernese 原有程序计算速度可提高13倍。
利用MGEX(multi-GNSS experiment)发布的BDS卫星差分码偏差(differential code bias, DCB)产品,比较分析不同太阳活动水平下BDS卫星DCB产品的稳定性变化特性,并采用差分自回归移动平均(auto-regressive integrated moving average, ARIMA)时间序列预测模型,实现不同太阳活动水平下BDS卫星DCB的短期预报。结果表明,在太阳活动高年,BDS卫星DCB日解值稳定度、月稳定性均明显低于太阳活动低年,且不同卫星星座类型的BDS卫星DCB稳定性也存在差异;ARIMA时间序列预报结果与MGEX发布值符合程度较好,优于多项式拟合法预测结果。
计算位于不同纬度的3个MGEX测站在不同太阳活动及地磁场活动情况下GPS、GLONASS、Beidou和Galileo卫星导航定位系统的电离层高阶项(HOI)延迟,分析高阶项延迟对多系统PPP结果的影响,总结不同太阳活动和地磁活动对PPP的影响。结果表明,电离层延迟随纬度的增加而减小;太阳活动是影响电离层高阶项延迟的主要因素,在太阳活动指数较低时,二阶项延迟不超过2 cm,三阶项延迟约为1 mm;当太阳指数较高时,二阶项延迟超过2.5 cm,三阶项延迟可达到5 mm;当磁暴发生时,高阶项延迟会存在一定程度的增长。将电离层高阶项延迟改正后,太阳活动指数较低时点位精度提升效果明显。
针对北斗混合系统中不同轨道类型、BDS-2/BDS-3卫星间观测值的差异,提出利用基于北斗UERE的Helmert方差分量估计模型对北斗混合系统定位的随机模型进行优化。对基于等权模型、普通Helmert方差分量估计模型及本文提出的优化模型进行单点定位实验,从定位精度与定位成功率上对3种模型的解算结果进行对比分析。结果表明,相比于等权模型和普通Helmert方差分量估计模型,优化模型能提高在E、N、U方向上的定位精度,同时减少因卫星星座空间几何构型较差而导致的定位精度较低的情况,从而提高定位成功率。