现代大陆动力学及构造学研究认为,壳内软弱带控制着地壳与地幔间相互作用的强弱,制约岩石圈层间的耦合关系,同时也影响地块的相互作用方式。这种耦合作用可以通过比较不同深度的形变和应力的相关特征来分析。青藏高原东南缘20~35 km的深度广泛存在低速带,有可能在现今地震及构造活动区内由壳内部分熔融承担了潜在的解耦功能。基于1999~2015年GPS坐标时间序列数据得到速度场,结合前人给出的SKS剪切波分裂数据,分析该区水平方向形变耦合特征;同时搜集莫霍面和岩石圈底面形变数据,给出垂向上的形变耦合特征。利用EGM2008重力模型的11~36阶球谐系数计算青藏高原东南缘岩石圈底面地幔对流应力场;搜集并整理2000年至今的1 131个震源机制解数据,采用区域应力张量阻尼反演法得到该区多震层应力场;进一步分析两种应力场的相关性。在ABAQUS软件中建立青藏高原东南缘岩石圈三维粘弹性有限元模型,同时加入上地壳的边界位移约束条件和底面地幔对流拖曳力加载条件,分别对水平方向和垂直方向的应力应变作分析,并给出强震剖面结果,探讨壳幔耦合强震孕育机制,据此分析青藏高原东南缘各个断层的地震危险性。研究认为,应力的不均匀发展是应变能积累的关键,同时垂向上圈层间的耦合关系极大地影响块体间的作用方式;壳内软弱带广泛发育的块体应力出现分层现象,边缘区域横向纵向应力变化明显,与之接触的上伏地壳地震危险性较高。
以河北省为例,利用GPS水汽进行MODIS水汽的校正模型研究。通过GPS水汽与MODIS水汽的比较发现,两者的相关系数超过0.879,均方根误差达到2.31 mm,说明两者存在较好的相关性和一定的偏差;利用线性回归方法分别构建城市和区域MODIS水汽校正模型,实验证明,利用实测GPS水汽对城市模型和区域模型进行可靠性验证,9个测站的模型检验的均方根误差小于1 mm。两类模型均可有效提高MODIS水汽精度,满足气象预报应用的要求。
对常用的3种天顶对流层延迟改正模型(Saastamoninen模型、Hopfield模型和EGNOS模型)进行误差分析,代入气象元素及测站位置误差,得出各模型的ZTD估值受误差影响的程度。使用C++语言实现以上3种模型,选取一系列不同纬度和高程的IGS站,利用IGS分析中心提供的气象文件,结合测站时空信息,导入程序进行模拟计算,并与IGS对流层产品进行比较,对改正模型进行质量评价。结果表明,Saastamoninen模型和Hopfield模型能够较准确地根据地面气象资料反映对流层延迟的日变化,Saastamoninen模型的改正精度略高于Hopfield模型;同时,无需实测气象资料的EGNOS模型RMS小于0.1 m,也可满足GNSS m级定位需求。
针对北斗广域差分系统实时电离层延迟格网改正信息是基于伪距观测计算,而相位观测值具有更高的观测精度,提出伪距数据结合相位历元间差分数据综合解算电离层延迟的算法。基于陆态网GPS实测数据对该伪距相位综合算法进行实验验证,结果表明,综合伪距相位算法相比单纯伪距法可以获得平滑、稳定、精度更高的站星斜路径方向电离层延迟值。在此基础上,全天所有实时电离层延迟格网的内符合精度平均提高30.2%,外符合精度平均提高27.8%。
采用非组合精密单点定位技术估算VTEC,并选取一系列与测站共址的COSMIC掩星电离层产品共同作为约束条件,利用函数模型近实时地探测GNSS测站上空电离层电子密度剖面特性。为了验证新方法的可行性,选取东亚11个IGS跟踪站2016年4个季节中3月、6月、9月和12月最后一周的观测数据及与测站共址的COSMIC掩星ionPrf产品进行实验。将每一个反演结果与同一时刻COSMIC掩星探测的电离层电子密度(作为真值)进行比较。结果表明,新方法实时反演的电子密度剖面与掩星电离层产品相比具有很好的一致性,所获得的f0F2与hmF2参数精度皆优于IRI2016模型。
引入一种改进的集合经验模式分解方法(improved ensemble empirical mode decomposition,IEEMD)进行多路径误差的建模。与EEMD方法相比较,IEEMD方法可以有效地克服模态分量数目不一致和分解不彻底等问题,并将原始序列分解为不同尺度的模态分量。同时,考虑到不同模态分量中高斯白噪声的能量密度与平均周期之积为常数,设计一种自动选择尺度与重构的方法,用于模态分量的选择与重构,进而构建GPS多路径误差模型。在此基础上,采用恒星日滤波技术,进行邻近周日GPS坐标序列中多路径误差的实时削弱。实验结果表明,采用同样的尺度选择方法和恒星日滤波技术,使用IEEMD方法可以得到比EEMD方法精度更高的GPS坐标序列。
基于多GNSS实验跟踪网(the multi-GNSS experiment,MGEX)的实测数据,对当前Galileo单系统标准单点定位(SPP)和精密单点定位(PPP)性能进行评估与分析。结果表明:1)Galileo双频标准单点定位的水平精度优于2 m,垂向精度优于4 m;2)GAL-201、GAL-202可用于SPP解算,引入二者之后,平均定位可用率提升10%,达到80%以上;3)当观测时长为3 h时,Galileo静态PPP可达到水平方向优于5 cm、垂向优于7 cm的定位精度;而随着观测时长的增加,其精度逐步提高,当观测时长为24 h时,水平精度优于1 cm,垂向精度优于3 cm。
利用资源三号01星(ZY3-01)实测GPS数据,基于简化动力学定轨方法和残差法估计该星GPS天线的在轨相位中心变化(phase center variation,PCV)模型,并分析该星PCV对精密定轨的影响。实验表明,通过考虑PCV模型,ZY3-01星定轨结果在重叠弧段对比上,三维位置RMS值有4.5 mm的精度提升,SLR检核RMS值有1.2 mm的精度提升,且各SLR测站检核结果也均有不同程度的精度提升。
针对空间直线拟合不宜直接采用总体最小二乘算法和混合总体最小二乘算法的问题,提出一种基于PEIV模型的总体最小二乘法的空间直线拟合算法。首先将空间直线的标准方程进行变换,改写为总体最小二乘的EIV模型;然后针对系数矩阵的特点,将模型转换为更加合理的PEIV模型,线性化成类似于最小二乘间接平差形式,采用迭代的方法求解拟合参数。平差过程保证了系数矩阵重复元素的改正数一致,常数元素的改正数为零,符合实际理论;最后,通过算例比较验证了该方法的可行性和优越性。
为分析网形结构对重力网平差结果的影响,利用北天山测网改造后的最新4期流动重力观测数据,分析测网中存在长支线对整网数据处理结果的影响。结果显示,长支线的存在使北天山最大异常变化达到260 μGal。将支线连成闭合环再进行平差计算,得到的最大重力变化量减小到60 μGal,表明测线上的误差积累导致支线末端测点重力值可靠性差,严重影响异常重力变化的识别。由于观测原因形成长支线时,为正确识别异常信息,采取适当的数据处理方法,将支线连接成闭合环,可以有效提高测网整体平差结果的可靠性。
针对捷联式GNSS/INS组合航空重力测量系统,比较在当地水平坐标系和惯性坐标系下重力扰动的估算方法及精度。结果表明,惯性坐标系下基于加速度更新的卡尔曼滤波方法比当地水平坐标系下基于位置、速度更新的卡尔曼滤波方法模型简捷且重力扰动的估算精度相当,在半波长分辨率为7.5 km时,6条重复测线重力扰动差值的标准差平均分别为1.21 mGal和1.27 mGal;对于相同测线,两种方法计算得到的重力扰动差值的标准差平均为0.65 mGal;两种方法求得的重力扰动中都存在明显的系统误差,将其按时间进行线性补偿后明显减小,测线间的最大差值分别由8.20 mGal、8.19 mGal减小到2.70 mGal、2.54 mGal,平均内符合精度分别从1.21 mGal、1.27 mGal提高到1.06 mGal和1.10 mGal。
用采集时钟方波信号加以分析的方法评估地震数据采集器的授时、守时服务精度、稳定性及系统特性。将具有时钟误差日志记录的数采的时间精度测试结果与日志记录相比较,二者完全一致;用时钟方波信号实测的数采经验传递函数与用厂家提供的滤波器系数并校正滤波器时间延迟计算的理论频率特性相符。实测结果显示,由于数采守时误差变化呈非线性,即使同一型号的不同数采间,其守时能力也可以差别很大;对于不同的数采及抽取滤波器,由于滤波延迟时间校正的参考点不一致,可导致数采记录数据的时间服务存在1个采样点左右的时间差,这对时间精度要求高的研究工作是不能接受的。
为揭示陕西石泉井不同频带水位对气压和固体潮的响应特征,采用频谱分析、相干函数和时移分析方法,对石泉井2015-11-01~2016-02-29的井水位、气压和理论重力固体潮等进行深入分析。结果表明,低频带(f<0.5 cpd)井水位对气压响应较为一般且波动较大,对固体潮响应很差;中频带(0.5~8 cpd)井水位对固体潮响应很好,同时,对气压成分中的K1、S2和S3频点有较好的响应;高频带(f>8 cpd)井水位与气压和固体潮的相干性均较差,这可能由于该频带井水位的信噪比较低和气压能量较弱等因素所致。此外,在全频带内,随着气压周期的增大,井水位的滞后时间也相应从1 min增至720 min;在中低频带的某些频点或频段,井水位对气压响应的时移存在超前和异常波动现象。
关键词: