在介绍经验映射函数NMF和动态映射函数VMF1、GMF模型的基础上，利用ECMWF提供的格网数据，具体研究和分析了这3种映射函数1 a中的时间变化特征及其随高程、纬度、高度角变化的特点，并比较了3种映射函数随高度角的变化对斜路径延迟估计的影响。

在分析观测数据解码以及观测数据质量指标的基础上编制软件，对采集到的原始GPS观测数据进行解码与多路径效应、周跳等质量指标的分析计算，并对质量指标进行直观、实时的图形显示。介绍所写软件在工程实践中的应用，并与TEQC软件进行比较，证明了软件的正确性和可靠性。

In GPS regional network coordinate time series, common mode errors (CME) are widespread. In order to extract CME, a sub-region is divided to get them superimposed and filtered. However, there is no principle to rely on for dividing sub-regions. According to the CME mentioned above, regarding three stations in southern California as centers, we divide the space into nine different spaces within different scales, carrying out the PCA space filtering method, to get the CME in each space scale. After analyzing and comparing the CME in each space scale, we can get the spatial response characteristics in CME, which provides a principle for the dividing sub-regions in the scale field.

介绍了广播星历精度评估的基本原理和方法，并利用2013-01~2015-04共28个月的广播星历数据，分析比较了北斗不同类型卫星的轨道精度、钟差精度及整体精度的短期和长期变化趋势。结果表明，GEO卫星广播星历的轨道精度约为1.8 m，钟差精度优于6 ns；IGSO和MEO广播星历的轨道精度优于1.2 m，钟差精度优于4 ns，整体上优于GEO卫星。从长期变化趋势来看，北斗广播星历的精度有逐渐提高的趋势。

对MGEX发布的北斗差分码偏差（differential code biases, DCB）产品的精度与稳定性进行质量分析，比较其日解产品与周解产品，并利用交叠式Allan方差初步计算分析北斗差分码偏差的随机误差特性。结果表明，北斗系统IGSO卫星DCB产品精度与稳定性最优，MEO次之，GEO最差。3种产品中，DCBB2B3的精度与稳定性最优，DCBB1B2次之，DCBB1B3最差。实验分析得出，北斗系统DCB随机误差组合主要为WN+GM，DCBB1B2、DCBB1B3与DCBB2B3的高斯白噪声分别为0.16 ns、0.15 ns、0.09 ns；一阶马尔可夫过程噪声分别为0.093 ns/d、0.049 ns/d、0.058 ns/d，相关时间分别为42.63 d、85.79 d、42.63 d。

根据GM(1，1) 、BP神经网络、卡尔曼滤波的特点，建立以灰色关联度最大为准则的最优非负变权组合预测模型，在一定程度上抑制了误差“放大”的效应。与以误差平方和最小为准则的变权组合模型、各最优加权组合模型、各单一模型进行对比分析，结果表明，本文模型预测精度最高。

将最小二乘配置方法的高斯型经验协方差函数模型应用于区域跨断层形变曲线的异常辅助识别。通过川滇地区跨断层形变对芦山7.0级、康定6.3级、鲁甸6.5级等地震的前兆异常判定曲线，说明最小二乘配置内插拟合可以进行区域跨断层水准形变曲线的趋势模拟和辅助查找区间异常特征，对于跨断层形变曲线异常识别有一定意义。

从观测值奇异值、观测值分布及反演模型的先验信息3个方面入手，研究观测值对反演结果的影响。以观测点到微分断层的距离倒数为权，构建顾及观测值权值影响的反演模型。以位错理论为基础，结合遗传算法，对龙门山断裂带进行三维滑动速率反演研究。结果表明，经过奇异值剔除的数据应用于顾及观测值权值的反演模型，给定合理的“先验信息”，其结果比已知观测值直接参与反演模型更稳定、可靠。

利用2004～2013年GPS速度场数据，分析青藏块体东北缘的应变场和剖面应变积累特征。结果显示：1）稳定的阿拉善地块及鄂尔多斯地块分别阻挡了祁连山断裂带与贺兰山构造区分别向NE、E方向的物质迁移；2）主压应变率、最大剪应变率高值区均分布在祁连山构造带东段、海原断裂东段至六盘山断裂西段，面膨胀率压性高值区主要在祁连山构造带东段，3）祁连山断裂东段显示为以压性为主、兼具左旋走滑，海原断裂东段至六盘山断裂西段则表现为以左旋走滑为主、压性为辅的相对运动特征。

根据渤海湾区域1999~2013年多期GPS站点速率资料，以营潍断裂带为界分为东西两个地块计算运动参数，求解主应变率和断裂带的应变率。基于应变率计算结果研究了区域应变-应力场。结果表明，1999~2011年渤海湾区域地壳整体受北东东-南西西向压应力和北北西-南南东向张应力作用，2011~2013年区域地壳压应力主方向为近南北向，张应力主方向呈近东西向，致使这一变化的主要原因是2011年日本3.11大地震。渤海湾区域不同地段的应力也存在差异，局部性变化明显。营潍断裂带不同时段的正应变率和剪应变率亦不同，1999~2011年为右旋剪应变-剪应力作用，2011~2013年为左旋剪应变-剪应力作用。

为研究汶川和芦山地震的孕育过程，利用单台单震法，将地震丛发区内发生的地震粗略视为“广义重复地震”，系统分析龙门山断裂带及附近区域1990～2015年波速比时空变化过程。结果发现，波速比在汶川地震前降低，芦山地震后明显恢复。对于龙门山断裂带北段的地震丛发区，大部分台站波速比的同步时间变化过程可分为两个阶段：一是在汶川地震前，2000年～2008-05波速比出现异常降低-回返-发震过程，低值异常持续时间约8 a。汶川地震后，波速比快速回升。二是2010年开始再次显现波速比持续性稳定下降，在2013-04降至最低点，发生芦山地震，震后波速比再次快速回升到正常值。对于龙门山断裂带南段及附近的地震丛发区，大部分台站波速比的同步时间变化从2000年～2013-04波速比维持长时间的异常低值状态，低值持续时间约13 a，直至2013-04芦山地震后波速比快速回升。不论是龙门山南段还是北段，均可看到波速比基本一致的变化，2000～2013年的异常降低-恢复过程，显示了汶川8级地震与芦山7级地震的整体性和差异性，孕震过程长达14 a。由此可见，汶川地震和芦山地震很可能是同一动力学过程孕育的两次先后发生的强震。

采用有限元方法模拟隧道地震波场，采用波场快照与时间记录相结合的方法，研究空洞对隧道地震波场传播的影响，并对含空洞模型的时间记录进行处理，得到数值模型的速度云图和反射层位图。经数据处理表明，采用TSPwin设定默认值处理得到的速度云图与模型设定的空洞位置比较一致，在提取反射层位图上，空洞反射层呈现条带状特征，需要结合速度云图来确定空洞位置，且P波预报的准确性相对较高。对TSP系统的抗噪性进行研究，表明其具有良好的抗噪性能。最后通过对工程实例的处理，验证了数值模拟所得的结论。

分别利用Helmert和KTH法模拟试算了不同地形起伏区域内横向密度扰动对于大地水准面的影响。结果表明，在地形起伏较小的平原或丘陵地带，横向密度扰动的影响一般不会超过cm级，在精度允许的范围内可以忽略其影响；在地形起伏较大的山区，横向密度扰动的影响可达cm甚至dm级，对于cm级大地水准面精化而言，需要考虑其影响。此外，基于Helmert与KTH法解算的结果在地形起伏剧烈的山区差异较大，由于Helmert法计算地形效应时忽略了外区地形质量的影响，其解算结果可能存在偏差。总体而言，在地形起伏剧烈的山区建模时，需利用多种计算方法基于实测数据分别进行试算，通过计算的变密度大地水准面与实测的GPS水准数据进行比较分析，得到适合该区域的大地水准面建模方法。

根据均匀分布在全球各大洲的IGS测站信息，选取105个数据质量良好的近海岸测站以及29个内陆测站，利用ECCO发布的海底压力数据计算得到由非潮汐海洋负载引起的IGS测站位置N、E、U方向上的位移。分析了由地理位置差异导致的测站受非潮汐海洋负载效应的影响存在的差异，并修正了SOPAC提供的IGS测站坐标时间序列，探讨了非潮汐海洋负载效应引起的测站位移对近海岸测站坐标时间序列的影响。结果表明，近海区域（<50 km）测站受到非潮汐海洋负载影响U方向平均量级为5 mm，远大于内陆（>500 km）测站U方向的影响；经过非潮汐海洋负载修正，71%的测站高程坐标时间序列RMS值得以减小，其中最大可使RMS值降低13%；大西洋、太平洋西海岸受其影响比东海岸大，考虑到其受非潮汐海洋负载效应影响显著，使用大西洋与太平洋西海岸地区测站坐标时间序列前必须进行非潮汐海洋负载修正。

对四平台、榆树台地电场北分量Ex和东分量Ey不同时段分钟值与长春台分量地磁、分量应变资料进行日变形态、日变幅等方面的对比分析，运用FFT将各观测资料的优势周期分解出来，以期进一步得到各物理量之间的相互关系，同时对地电场在地震前的对应性进行相关总结。结果表明：1）两台地电场日变化都以双峰单谷为主；两测向优势周期都是以12 h、8 h为主。与分量磁场相比，东分量By与Ex比较吻合，北分量Bx与Ey的对应性很高；与分量应变相比，12 h的优势周期为共有周期。2）两台地电场都能够记录到明显的地电暴同步变化。3）Ey与Bx的日变幅比较同步，Ex与By的日变幅比较同步，且呈现出较明显的季节变化。4）两台地电场都能够较为清晰地记录到一定的短临异常，但两台震前异常表现形式不尽相同，出现异常的时间也不同：震中距越近，出现异常越早，异常幅度也越大。

在保留RZB型电容式分量钻孔应变仪传感器功能不变的前提下，充分吸取了原RZB型电容式多分量钻孔应变仪及其他前兆仪器的优点，对其主机部分进行升级改造，实现了RZB多分量的宽频带响应（采样率提升为秒采样）；采用了低功耗、高稳定性的ARM920T高速处理器内核、固化WinCE操作系统，为仪器的稳定性、连续性提供了技术保障；网络接口板外嵌触摸式显示屏，极大地方便了仪器的参数配置以及现场检修；对原供电方式进行改进，实现了集防雷、供电、充电、交直流切换于一体以及软件的远程更新、软件工作的详细日志查询等功能。通过对新仪器数据分析发现，其连续率得到大幅提升，并记录到分钟采样所记录不到的详细地震波响应过程。

分析丹江口井网3口井水位、水温2013～2015年的观测数据，建立各个测项的正常动态，总结井网对2015-04-25尼泊尔M8.1地震的响应特征，探讨同震响应差异的可能原因，初步得到水位同震响应幅度与震级、井震距的关系。结果表明，井网水位年动态受降雨影响有一定规律，月正常动态普遍显示出固体潮效应与气压效应；水温年动态平稳，年变幅度在0.1 ℃左右。对于尼泊尔地震同震效应，唐扒井水位阶升-阶降式脉冲变化和水温陡降-缓升式脉冲变化是由水位和水温的变化机理不同导致的；雷庄井与唐扒井水位表现出不同的同震响应特征，判断与它们所处断裂的不同走向有关；盛湾井水温阶降型同震响应幅度较唐扒井小得多，可能是观测层岩性及地下水类型等差异导致的。

对地震预警信息发布时延作调研后认为，预警接收终端的发布时延小于1 s；广播、电视发布时延为几s；互联网发布（新媒体、手机APP、第三方平台推送）时效性受并发量影响，在千万级用户时发布时延是30～40 s。因此，应采用分区域分级发布策略，以提高发布时效。以汶川地震为例，针对社会公众发布预警信息的时效分析结果表明，在首台触发3～4 s内开始发布预警信息，震中20～50 km范围内95%以上的网民有1～7 s的预警时间，而破裂最大区域（北川）95%以上的网民将获得30 s以上的预警时间。