综合运用InSAR和机器学习技术，对四川省金阳县进行滑坡易发性评价。通过解译滑坡更新数据集，基于12个评价因子，在Python环境下使用随机森林(RF)、支持向量机(SVM)和极端梯度提升(XGBoost)3种模型进行训练，完成滑坡易发性制图，并采用ROC曲线等验证模型预测性能。对负样本进行优化，使用机器学习得到样本优化后的滑坡易发性评价结果，并利用地表LOS向形变速率对其进行更新。结果表明，3种机器学习模型均具有较好的分区效果，XGBoost模型制图效果最佳，样本优化后模型精度最高，AUC值达到0.95。通过SBAS-InSAR技术获取地表形变速率可以减少分区错误，同时赋予滑坡易发性评价结果时效性。

基于溪洛渡水电站的可移动北斗变形监测系统，利用GNSS测站水平变形模拟实验，开展静态差分、事后动态差分(PPK)和动态精密单点定位(PPP)3种方法对比分析与形变监测精度评估。结果表明，单BDS与GPS+BDS双系统的观测精度相当，在3种解算方法中，静态差分精度最高，双系统和单BDS的外符合精度分别优于0.5 mm和0.7 mm，基线差分方法通过24 h长时序观测能够更加有效地消除坝区水汽干扰和电离层及多路径误差；动态PPK双系统和单BDS结果的最大误差分别为0.5 mm和1 mm，动态PPP分别为1.7 mm和1 mm，说明PPK和PPP也能获得精度接近1 mm的形变监测结果，具有较大的应用潜力。

利用IGS发布的精密钟差数据，对已发射并在轨运行的6颗GPS_Ⅲ星钟的频率准确度、频率漂移率和频率稳定度等关键性能指标进行分析，并与前期GPS星钟进行对比。结果表明，GPS_Ⅲ搭载的铷钟整体运行稳定且数据质量较好，可能与抗干扰能力的提升有关，时频性能并无重大变化或明显不同，基本与上一代Block ⅡF铷钟一致，但频率稳定度更优，有助于提升PNT服务精度。研究成果可为其他GNSS系统空间段卫星的升级改造提供借鉴和参考。

在对滑坡监测数据粗差进行有效处理及充分顾及滑坡监测数据自身特性的基础上，提出一种结合时序分解与相似分量重组的深度学习滑坡位移组合预测模型。首先，利用孤立森林法对滑坡时序监测数据的显著粗差进行处理，再对其平稳性、自相关性、正态性进行综合分析，确定模型预测中输入特征序列的最佳长度；其次，利用集合经验模态分解(EEMD)方法，将非稳态滑坡监测数据分解为多个平稳时间序列，再结合样本熵与K-means算法将其划分为高频、中频、低频3类时间分量；最后，通过对比不同神经网络模型的预测精度，分别构建适合于3类时间分量的预测模型，再将预测结果相叠加，实现对滑坡位移的高精度预测。实验区典型滑坡体北斗/GNSS监测数据测试表明，本文组合预测模型对含有显著粗差的滑坡监测数据具有较好的适用性，相较于单一及现有组合模型可显著提高滑坡位移预测精度。

针对地下水位数据的复杂特性(包括非线性趋势、季节性波动和随机扰动)，引入Facebook开发的Prophet时间序列预测模型，旨在利用其非线性趋势捕捉、季节性波动解析及对异常值和数据缺失的灵活应对能力，显著提升地下水位异常识别的准确性。黑龙江省绥化市北林区地震台观测数据表明，Prophet模型在捕捉时间序列动态特征上表现优越，能有效识别异常。模型调整后具有高拟合精度和高预测能力，预测误差低，决定系数高。此外，模型在地震预测中能识别出与地震相关的水位异常，可为地震前兆研究提供新视角。本文结果表明Prophet模型在处理复杂时间序列数据时具有可行性，可为地震预测提供新工具。

建立BDS传递起算基准控制测量模型，比较分析BDS-3/BDS-2和BDS-3+BDS-2高精度控制测量定位性能差异。选择北京某工程9个首级控制点和中国及周边3个IGS跟踪站建立一种同时支持GPS与BDS-3/BDS-2信号的高精度控制测量定位网，选取2023年doy153-155连续观测的3个时段数据，设计4种方案进行数据处理，从起算基准传递结果和控制测量平差结果2个方面作定位性能比较分析。结果表明，相比于GPS，BDS-3和BDS-3+BDS-2在N向、E向、U向、平面和点位上较差均仅存在mm级的差异，而BDS-2在U向和点位上较差存在mm至cm级差异，在N向、E向及平面上存在mm级差异; BDS-3测量精度与BDS-3+BDS-2基本相当，BDS-2测量精度稍差，特别是在U向。

利用精密水准测量数据，结合改进的倾滑位移模型，系统分析西秦岭北缘断裂带各区段的垂向运动特征。结果表明，西秦岭北缘断裂带正经历左旋走滑与逆冲挤压的复合变形机制，且不同区段的活动性存在显著差异。垂向倾滑速率呈现从西段的0.24 mm/a向东增大到2.31 mm/a后，又减小到1.51 mm/a的空间变化趋势。通过反演获得的闭锁深度与地震精定位所揭示的小震活动深度高度吻合，说明该断裂带各区段处于不同的闭锁状态。结合地震活动性发现，天水-宝鸡段发震危险性最紧迫。

提出一种适用于复杂环境下低成本接收机精密单点定位(PPP)的自适应随机模型，该模型通过对单个历元所有GNSS卫星的伪距观测值和载波相位观测值的权重比进行自适应调节，以提升PPP性能。采用和芯星通UM982型号低成本GNSS接收机，在树荫、高楼和玻璃墙3种复杂环境下进行PPP实验。静态PPP实验结果表明，相比传统的经验定权随机模型，自适应随机模型定位精度在3种复杂环境下分别提升24%、45%和50%，收敛时间分别缩短49%、27%和24%。动态PPP实验结果表明，自适应随机模型定位精度在复杂环境下提升35%。

GNSS技术被广泛应用于滑坡监测，但在复杂和高危环境下，监测设备难以安装，而利用无人机抛投有望实现监测设备的无人化部署。针对无人机续航以及复杂环境下地形威胁等问题，抛投无人机的路径规划作为基础任务尤为关键。将复杂山区抛投无人机路径规划问题分为地质灾害区域三维实景地图构建、代价函数设计、航迹规划算法3个方面展开研究，并将基于自适应权重和Levy飞行策略的鲸鱼优化算法应用于抛投无人机航迹规划。以四川省汉源县红岩子滑坡为研究区，实现了复杂山区GNSS监测装备抛投无人机的路径规划。

为降低复杂环境下GNSS定位误差，提出一种联合高精度测站和距离交会精确估计定位点坐标的方法。该方法首先将观测方程构建为非线性高斯-赫尔默特模型，针对其中的非线性问题，引入反向传播(back-propagation，BP)神经网络进行辅助处理。与传统线性化方法相比，BP神经网络能够有效拟合复杂的非线性函数关系。仿真和实测结果表明，该方法能有效降低复杂环境对定位精度的影响，E、N、U方向定位精度分别提高78.1%、72.8%、79.2%。

针对地震前兆观测中的地电阻率法，研制出适合其观测系统的检测装置，并对装置进行理论仿真与实地检测。首先，通过对地电阻率测量系统进行剖析，指出现有检测系统的不足；然后，根据观测系统需求进行分析，设计对应的检测指标和电路，并将仪器检测精度设计至0.001级别，提升设备性能；最后，借助Simulink工具对检测装置进行理论仿真，对比研究干扰对观测结果的影响，并随机选取甘肃地区地震观测台站进行实地检测。结果表明，检测装置不仅能直观展示地电地阻率测量系统的工作状态，而且能精准评价地电阻率观测系统的实时优劣性，填补地电阻率观测系统检测装置的空白，为地震一线工作人员提供极大便利。

收集日本在其周边海域及西北太平洋通过船载调查获取的227个航次海洋重力数据集，应用统计分析法、可视成图法、相关分析法等开展日本重力调查情况及数据质量情况分析，系统分析日本重力数据调查历史、时空分布特征及质量情况。结果显示，船测重力数据与GEBCO地形数据、DTU重力模型数据具有良好的空间一致性，其中与DTU模型的相关系数为0.98，均方根误差为7.7 mGal。本文建立了航迹交点差值模型，并利用交点差误差进行数据测量精度的定量分析，总交点共73 979个，误差为11.85 mGal，内部交点24 822个，误差为9.56 mGal，外部交点49 157个，误差为12.85 mGal，数据测量准确度总体较高，数据一致性较好。

将生成对抗网络(generative adversarial networks, GAN)应用于GPS坐标时间序列的降噪处理中, 通过生成器生成时间序列数据, 并利用判别器对生成数据与干净的时间序列数据进行区分。在训练过程中, 不断迭代优化生成器, 直至判别器无法区分生成数据与真实干净数据, 从而实现降噪效果。仿真实验结果表明, 与WNNM、BiLSTM和SSA等传统方法相比, GAN方法在残差振幅、残差谱指数和速度不确定度等关键指标上与原始噪声的相关指标最接近, 展现出更高效的降噪性能。真实数据的实验结果进一步验证了本文算法在实际GPS坐标时间序列降噪中的有效性。

以美国德克萨斯大学空间研究中心提供的地心运动时间序列为实验数据，首先采用经验模态分解(empirical mode decomposition, EMD)对数据进行降噪处理；然后利用快速傅里叶变换(fast Fourier transform, FFT)和连续小波变换(continuous wavelet transform, CWT)对该数据进行频域转换、功率谱分析和周期项提取；最后采用自回归积分滑动平均(autoregressive integrated moving average, ARIMA)模型和指数平滑法对未来20个月的地心运动进行预报。结果表明，利用FFT提取的周年项的振幅和相位与以往的地心运动研究较为接近；ARIMA模型对于Y方向20个月内的地心运动预测结果较好，指数平滑法对X、Z方向的地心运动预测结果更优。

本文使用Jason-1/2/3系列卫星的VTEC数据, 从时间、空间等角度分析全球电离层地图的精度情况, 并利用太阳活动分量、日变化分量和季节变化分量建立相应的修正方程。结果表明: 1) CODE GIMs与Jason系列卫星的残差(G-J)与F_10.7呈二次函数关系, 在日变化上呈24 h的周期特征, 在季节变化上呈半年变化周期特征, 修正方程能够有效修正GIMs数据。2) G-J的平均值在太平洋的低纬度地区取得最大值, 在澳大利亚南部取得最小值; 赤道异常中的等离子输送导致G-J的中误差呈现双峰特征; G-J的相对误差在大西洋南部, 尤其是非洲南部的海洋地区取得极大值。

对2023年跨唐山断裂带北段沿北西-南东向实施的一条长约30 km的深反射剖面的初至波走时进行反演，获得唐山市北部的浅层精细速度结构剖面。速度结构剖面揭示出，测线坐标2.5~8.5 km新生代沉积厚约150 m，较为均匀；丰台-野鸡坨断裂隐伏于测线坐标3.5 km处，倾向北西，F_2-1断裂隐伏于测线坐标8.5 km处，倾向南东；这2条断裂之间存在一个北西低南东高的隐伏凸起。测线坐标8.5~15.2 km，自北西向南东新生代沉积逐渐加厚，最大厚度约250 m，显示出荆各庄凹陷为箕状新生代断陷。陡河断裂出露于15.2 km处，倾向北西。测线坐标15.2~24 km为唐山凸起，仅局部区域有很薄的新生代沉积，同等深度上，P波速度远高于两侧；巍山下方有一近直立的低速带，向下延伸至少600 m，推测为唐山-巍山-长山南坡断裂带的破碎带。测线坐标24 km至测线端部为开平凹陷的西北部，表现为一规整的盆地边缘形态，新生代沉积约130 m。陡河断裂、唐山-巍山-长山断裂为第四纪活动断裂。

在同一场地分别采用T型、直线型、正五边型和螺旋型4种不同观测台阵进行微动探测对比实验，分析互相关分段时间窗口长度、观测台阵几何形状及其展布对频率-贝塞尔(F-J)频散能量成像质量的影响。结果表明，对于相同时间长度的观测，增加互相关分段时间窗口长度能有效区分基阶和高阶频散曲线，但会减少叠加次数，降低频散谱分辨率；F-J方法对观测台阵排列无特别要求，随机分布、线性排列等均适用，在台阵口径相同的情况下，拥有更多台间距组合的观测台阵的频散成像质量和分辨率会更高；F-J方法对观测台阵的展布范围有一定要求，当展布范围过小时，F-J方法频散成像质量较差，无法分离高阶频散曲线。

基于2023年26 496个不同观测时长的北斗监测序列，构建北斗监测结果可用率、监测精度与观测时长的定量函数表达式。结果表明：1) 监测结果可用率与观测时长存在明显的正相关，6个监测点的平均相关系数为0.73；N、E、U方向的监测精度与观测时长存在明显的负相关，相关系数分别为-0.83、-0.90和-0.77。2)为保证监测结果的高时效、高可用率(>95%)和高精度(平面优于3 mm、高程优于5 mm)，建议最佳观测时长为3 h。

利用南天山地区2020—2023年流动重力观测资料，获取2024-01-23乌什M_S7.1地震前不同时间尺度下区域重力场动态演化特征，并通过功率谱分析方法，获取各阶小波重力细节对应的近似场源深度。结果表明：1)流动重力(3 a尺度)结果显示，乌什M_S7.1地震前，乌恰至巴楚地区和阿克苏地区重力变化呈现明显的四象限分布，震中位于四象限边缘及零值线附近；2)2020—2023年南天山地区重力场小波变换(4阶小波重力细节)结果显示，乌什M_S7.1地震前，乌恰至巴楚地区重力变化出现明显的四象限分布，震中位于四象限边缘及零值线附近。

利用内蒙古地区728个震源机制解和2 793次地震事件的14 460条P波极性数据，采用3种不同策略对内蒙古地区应力场进行研究，获得地壳深层应力状态。将内蒙古地区划分为1°×1°网格，分别采用基于初动极性数据的综合震源机制方法、基于震源机制的贝叶斯右手三棱锥方法(BRTM)和极性数据反演结果与震源机制解数据加权融合方法进行反演。结果表明，基于极性数据和基于震源机制的方法具有较好的一致性。由于极性数据更为丰富，可克服震源机制解数据的有限性，获得更多区域的反演结果。将基于极性数据的反演结果与震源机制解数据加权融合后，采用BRTM方法获得125个格点的反演结果，可扩充计算区域和降低结果的不确定性。研究结果显示，内蒙古地区不同区域的应力场在空间上存在较大差异，东部地区主要表现为剪切和拉张特征，大兴安岭主脊表现为剪切特征，两侧盆地表现出拉张特征，东南部的赤峰和通辽一带局部呈现挤压特征，鄂尔多斯地块北缘主要表现出拉张特征，阿拉善地块主要表现出剪切特征；最大主应力方向存在较大差异，而最小主应力方向较为一致。

基于欧洲轨道测定中心(CODE)总电子含量(TEC)数据和张衡一号卫星电子密度(N_e)数据，分析2020—2022年中国西部3次M_S≥6.0浅源地震(2020年伽师M_S6.4、2021年漾濞M_S6.4和2022年马尔康M_S6.0地震)前的电离层扰动现象。结果表明，电离层异常集中在震前1周内，并且在伽师M_S6.4地震前5 d和2 d、漾濞M_S6.4地震前7 d和马尔康M_S6.0地震前6 d出现CODE TEC和张衡一号卫星N_e同步异常。这3次强震可能通过岩石圈-大气层-电离层耦合机理模型的地球化学和电磁波通道影响电离层，震前5~15 d内震源区附近出现地球化学、水位和电磁等地震前兆异常。

联合青藏高原东南缘地区最新的水平与垂直GPS速度场，基于球面小波多尺度分析的应变求解方法，考虑区域地震矩与应变分布的权重关系，得到青藏高原东南缘三维地壳应变特征，并分析区域形变特征与地震活动的关系。结果表明：1)研究区面应变与垂直应变特征具有一致性。高面应变与低垂直应变主要集中在青藏高原腹地和金沙江北部区域，低面应变与高垂直应变主要集中在喜马拉雅东构造结与龙门山断裂。两者在空间上呈反比分布，说明GPS观测得到的水平速度与垂直速度具有耦合性，即对于该区域地壳形变与地震危险性研究，必须考虑垂向运动因素；两者在量级上也存在较好的一致性，进一步说明该区域的地壳变形可以近似认为是连续均质的弹性行为。2)应变速率类型与震源机制解具有一致性。由GPS解算得到的青藏高原东南缘应变类型可以很好地匹配该区域震源机制解。正断层震源机制主要分布在金沙江地带，逆冲断层震源机制主要分布在龙门山断裂与喜马拉雅东构造结区域，走滑断层震源机制主要集中在鲜水河-小江断裂带等走滑剪切作用显著的断裂带附近，说明由GPS观测得到的应变特征可以较好地描述该区域构造应力场特征。结合该区域应力应变空间分布特征和地震层析成像等结果，认为龙门山断裂与喜马拉雅东构造结处于挤压应变状态。

太阳光压摄动作为在轨导航卫星受到的最大非保守力，是卫星精密定轨的重要误差源，而目前多数光压模型是针对全球定位系统(GPS)卫星而建立的，对北斗三号(BDS-3)卫星精密定轨的适用性分析较少。本文借助MGEX观测数据对5种欧洲定轨中心拓展轨道模型(ECOM1-9、ECOM1-7、ECOM1-5、ECOM2、ECOMC)太阳光压模型进行定轨实验。结果表明，不同厂商研制的卫星对光压模型的适用性不同，对中国空间技术研究院(CAST)研制的卫星而言，ECOMC模型在径向、切向和法向均表现出最优适用性；对上海微小卫星工程中心(SECM)研制的卫星而言，径向、切向和法向分别为ECOM1-5、ECOMC、ECOMC模型精度最佳；对三维整体而言，CAST卫星和SECM卫星定轨均推荐采用ECOMC模型；从24 h预测结果来看，对于CAST卫星，ECOM1-9、ECOM1-7、ECOM1-5、ECOM2、ECOMC模型24 h预测精度分别为23.3 cm、20.6 cm、17.2 cm、21.8 cm、10.4 cm，对于SECM卫星，各模型预测精度分别为25.7 cm、19.1 cm、15.9 cm、12.9 cm、11.9 cm。

收集山西长治及周边地区2009年以来M_L≥2.0地震事件的波形数据, 利用Snoke初动符号加振幅比方法计算其震源机制解, 最终获得其中M_L≥2.4且RMS≤0.45的42条天然地震和6条非天然地震的震源机制解。结果显示, 研究区天然地震震源类型以走滑和正走滑型为主, 主要受NEE-SWW向挤压和NNW-SSE向拉张作用; 长治盆地区具有复杂而独特的小区域应力场特征, T轴近NS向, P轴沿文王山断裂为NEE-SWW向, 其他区域以NE-SW为优势方向, 并显示出震源节面具有较好的一致性和高倾角特征。探讨研究区内非天然地震的震源机制, 发现矿山诱发地震多为正断型, 塌陷为逆断型, 其应力轴特征表现出P轴或T轴的大幅偏转, 两轴夹角约180°, 与天然地震的P轴和T轴接近垂直存在明显不同, 这些特征可辅助判识地震类型。

以2010—2021年格陵兰地区的35个GNSS测站的小时分辨率的天顶对流层延迟(ZTD)数据为基础，结合ERA5再分析资料和T_m模型，对该区域的大气可降水量(PWV)进行求解。随后，对其周年、半周年、周日、半周日等主要季节和日变化特征进行研究。结果表明，格陵兰地区测站存在一定程度的自相关性，且夏季自相关程度弱于其他季节。周年项振幅范围为2~7 mm，半周年项振幅范围为0~3 mm，其相位在年际间变化较小，最高值出现在夏季，最低值出现在冬季。仅有少数几个测站周日项比较明显，振幅在0.08 mm左右。去除时间序列中的周年、半周年、周日项得到残差时间序列，最大的RMSE范围为1.30~3.65 mm，平均RMSE为2.08 mm，其中，北部地区RMSE偏大，南部地区较小。

利用山西及邻区2009-01—2023-06的天然地震到时数据，采用双差成像方法获得该地区重定位目录、较高分辨率的地壳三维速度结构和波速比结构，讨论该区域6级以上历史地震的孕震环境。穿过山西南部地区6级以上历史地震的速度和波速比与深度的二维剖面结果显示，1303年洪洞8级地震处在高低速异常梯度带上，地震下方存在低速异常，表明韧性的中下地壳蠕变会导致脆性上地壳的应力应变局部集中而引发地震；该地震处于波速比低值区域，且临汾盆地下方为低波速比，其两侧山区波速比略高，这与盆地下方地壳经历伸展减薄的作用过程有关。速度和波速比结构均显示，罗云山山前断裂下方呈现出浅陡深缓的产状特征，断裂切割可延伸至约30 km深度，这与已有的地震反射剖面结果较为相似，可能是该处多次发生6级以上地震的原因。山西南部多数历史地震发生在高低波速比梯度带附近，可能与这些区域位于盆山边界有关，因其地壳物质组成变化较为剧烈而易于积累应力应变，从而产生地震。

为同时提高点云平面分割的效率与可靠性，在基于多尺度超体素及区域生长算法的基础上，以超体素为生长单元，采用体素区域生长算法实现精细点云平面分割。结果表明，本文方法对室内外点云数据都有较好的分割效果，优于RANSAC、Khaloo及全局能量优化方法，且本文方法的Precision、Recall和F₁-score都超过0.90。

提出一种多故障条件下RAIM完好性风险计算方法。该方法基于最小二乘残差构建卡方检验统计量进行故障检测，并进一步根据检验统计量和定位误差的概率分布，通过搜索多个故障模式下最差情况故障向量实现完好性风险计算。实验结果表明，该方法能够有效监测多故障完好性风险，在观测几何较差时段实现风险告警，保障用户定位可靠。在航空LPV-200阶段的完好性需求下，用户定位结果的可用性达到99.47%。

基于2023年积石山M_S6.2地震震区Sentinel-1A的雷达影像数据，利用D-InSAR技术获取升降轨同震形变场，使用SDM反演程序获取发震断层滑动分布，使用Coulomb3.3程序获取区域断层库仑应力变化。结果表明，本次地震产生的升降轨同震形变场均以抬升为主，升轨最大LOS向形变量为6.8 cm，降轨最大LOS向形变量为7.6 cm。发震断层参数最优解为走向303°、倾角52°、滑动角89°，最大滑移量为0.5 m，发生在距地表 11.89 km深度处，破裂未到达地表，累计释放地震矩1.48×10¹⁸ Nm，合矩震级M_W6.1。库仑应力计算结果显示，拉脊山南缘断裂、拉脊山北缘断裂的西北段、倒淌河-临夏断裂西北段、西秦岭北缘断裂西北段处于明显的应力加载状态，以上区域在未来一段时间内发生地震的可能性较大，需要重点关注。

提出一种基于自适应噪声完备集合经验模态分解(complete ensemble empirical mode decomposition with adaptive noise, CEEMDAN)和时间卷积网络-注意力机制(temporal convolutional network-attention mechanism, TCN-Attention)算法的多尺度预测模型(简称C-TCN-A)，该模型可有效应用于GNSS高程时间序列缺失数据的插补和未来趋势的预测。该模型利用CEEMDAN对时间序列进行多尺度分解，然后基于TCN-Attention对不同尺度分量进行预测和重构得到预测结果。为验证模型的性能，选取12个观测站进行1 d与5 d预测，并与其他多种模型进行对比。结果表明，在1 d预测中，C-TCN-A的RMSE和MAE分别降低35%~40%和36%~41%，相关系数R提高25%~29%;在5 d预测中，C-TCN-A的RMSE和MAE分别降低20%~26%和20%~28%，相关系数R提高26%~33%。为验证模型的普适性，利用C-TCN-A对陆态网络99个观测站进行1 d与5 d预测。结果表明，RMSE和MAE指标总体上结果较好，误差分布集中，大多数误差小于4 mm; 预测精度存在一定的空间分布差异，西北地区效果最佳。

通过将区域地壳应力的强弱与其发震概率进行联系，提出以视应力为基础数据来计算区域发震概率的方法，主要得出以下结论：1)运用多重分形理论计算奇异值α, 揭示研究区隐含的应力富集区、贫乏区和背景区。2)推导奇异值α与b值拟合的新方法，得到奇异值背景下的b值，使地壳应力强弱与b值得到统一。3)运用产生随机数的方法计算参考区剪应力加载速率值，并运用指定公式计算其他剪应力加载速率值。4)运用对数正态分布模型计算地震发生概率时，使用定积分进行模型参数处理，结果发现，剪应力加载速率值的计算误差是影响结果不确定性的主要来源。5)与正态分布模型和BPT模型进行对比验证发现，对数正态分布模型与正态分布模型计算结果吻合较好，而与BPT模型差别较大，这可能是各模型算法不同导致的。6)研究区5级以上地震发生概率较高的地区在沂源-临朐-安丘-诸城-沂水围成的区域、沂水-五莲-莒县-沂南围成的区域和研究区北部临朐-昌乐一带。10 a发震的条件概率最高可达0.034, 累积概率可达0.861;30 a发震的条件概率最高可达0.097, 累积概率可达0.87;50 a发震的条件概率最高可达0.15, 累积概率可达0.878。

提出一种综合外部大气数据与地形数据的一体化InSAR大气延迟校正方法。在引入外部大气数据对InSAR干涉对进行大气延迟粗校正的基础上，建立区域高程与干涉相位相关的全局校正模型，抑制区域长波分布地形相关大气延迟分量；再对干涉图进行分块处理，并对每个分块建立局部最适应的地形相关大气延迟误差模型，实现对各子块的精细化大气延迟误差校正。以滑坡灾害密集发育的杂谷脑河流域为例，利用不同的大气校正方法开展InSAR干涉对的大气延迟误差校正与滑坡识别。结果显示，随着大气校正的逐步开展，InSAR相位与地形的相关性呈显著下降趋势，验证了一体化大气延迟校正方法的有效性。区域时序形变解算与滑坡灾害识别显示，一体化大气延迟校正方法可有效恢复中小型滑坡形变信号，并显著有利于滑坡变形边界的准确探测。统计发现，相比于单一外部大气校正，联合全局地形相关大气校正识别滑坡数量增加了75%，而进一步引入局部自适应大气校正后，识别滑坡数量又增加了33%。

运用小波分析对昌黎台1981—2020年地电阻率和面应变固体潮理论值进行相干谱计算。结果显示，整点值计算结果12 h、24 h周期谱相关性较显著，频次约在1~4次/月，每次持续2~3 d左右，但出现时间不规律；各周期成分在时域变化中正负相关性比较杂乱，无稳定性。日值计算结果主要存在两类周期成分，一是8~30 d内的8 d、10 d、15 d、30 d短周期，二是较长的半年、年周期。其中，短周期成分连续性差，频次约3~5次/a, 在时域变化中正负相关性较杂乱；半年和年尺度的长周期相关性相对显著，持续时间多连续且变化相对稳定，年周期成分出现时间多于半年周期。年周期相关性最稳定，呈负相关，半年周期相关性稍差，正、负相关均有。日值小波相干谱还显示了长、短周期空区的存在，与震情有关，主要反映了本辖区应力场变化和远场大震的调制作用。

基于StaMPS和SBAS-InSAR技术获取2017-03—2022-10河北隆尧地裂缝及周边的地表形变结果。结果显示，2种技术监测结果较为一致，说明结果可靠；巨鹿县周边的地面沉降较严重，已形成较大规模沉降区，最大形变速率超过-30 mm/a; 提取的时序结果均呈周期性变化，且与地下水水位和降水量变化一致性较高；隆尧地裂缝两侧地表的最大形变速率差超过40 mm/a, 最大形变量差超过200 mm, 实地调查发现，地裂缝引发的地表差异形变对当地造成巨大影响。

高程基准是测绘基准的重要组成部分，可为国家重大工程提供数据支撑。针对高程基准因受地表形变等因素影响而难以精准监测的问题，本文基于时序InSAR地表形变监测和水准数据时间序列分析地表形变结果，对形变区开展高程基准周期性复测，实现以高程基准成果更新、点位可用性分析、点位沉降情况分析为主要内容的动态监测，并以山东省为例进行实验。结果表明，该方法能有效降低监测成本和监测难度，成果可靠性强，可有效提升测绘基准保障能力和服务质量，为区域高程基准动态监测提供一种新思路。

为了解决传统硬阈值方法在处理异常测站时难以应对模糊边界的问题, 提出一种新方法, 将构造运动和观测误差引起的测站异常均视为粗差, 并利用隶属度函数构建IGGIII等价权函数, 以减小异常测站的影响。新方法引入模糊集理论, 是一种基于标准化残差属于粗差的模糊子集隶属度函数改进的IGGIII抗差估计方法(membership function IGGIII, MF-IGGIII)。首先, 通过模糊统计确定标准化残差受粗差污染程度最大的模糊子集的隶属度函数; 然后, 依据这些隶属度函数值构建等价权函数, 并进行选权迭代。模拟实验和对环渤海区域GNSS速度场的拟合分析结果表明, 该方法在异常测站检测方面的表现优于假设检验法、均值漂移法和IGGIII法, 显示出较强的异常测站识别能力。

针对无线节点地震仪在物探与地震监测过程中实时交互能力不足、功能匮乏的现状，设计基于Android的节点地震仪监控与野外布设综合平台。平台采用低功耗蓝牙(bluetooth low energy, BLE)通信技术，通过自定义加密无线通信协议实现移动终端与节点式仪器的数据通信与设备控制，完成对单点设备的单/三分量地震动速度信号实时显示、状态查询与单台/多台同步设置；应用超高频射频识别技术、QR扫描和电子地图实现野外快速布设节点地震仪阵列测线；应用SQLite数据库和本地文件系统实现测线可靠记录、参数设置和系统日志的管理。经测试，平台可提高节点地震仪器系统野外物探勘测和流动监测工作效率。

地球内部圈层存在各向异性, 且各向异性成因机制不同。剪切波分裂可揭示不同构造、不同深度各向异性介质特征, 为地球深部动力学演化研究和地震预测提供重要信息。本文介绍地球内部各向异性成因机制, 总结剪切波分裂方法发展情况, 以及利用剪切波分裂方法研究地壳和地幔介质各向异性的成果与进展, 并在此基础上综合分析壳幔耦合模式。

分析在不同重力场模型下的定轨精度与轨道外推精度发现, 重力场模型阶数为100与150时的定轨精度和轨道外推精度均高于50阶重力场模型的结果。在GRACE-C和GRACE-D简化动力学定轨中, EGM2008与Tongji-Grace02k的3D方向定轨精度最高可达22.6 mm和24.0 mm, 这2种重力场模型定轨精度均优于ITG-Grace03。在轨道外推中, Tongji-Grace02k精度最高, 能达到亚dm级别; 其次是EGM2008, ITG-Grace03精度最低, 二者均在dm级。通过优化随机脉冲参数, 可以提升3种重力场的轨道外推精度, 模型阶数为100与150时的优化效果比阶数为50时明显。

利用音频大地电磁测深(AMT)与高密度电阻率相结合的方法，在辽宁新城子深井进行地电阻率观测，探测地下电性结构分布情况。通过精细探测，圈定了含水区的空间位置、范围和形态等关键参数，为多极距多层地震地电阻率场地选择和设计提供了重要参考。