在考虑地球曲率、成层结构、可压缩性与自重的前提下,Tanaka等[1-2]提出一套较为完备的粘弹性球体位错理论,可计算全球任意位置由地震产生的同震与震后形变(含位移、重力变化、大地水准面变化)。Gao等[3]给出了与上述理论相匹配、界面友好的计算软件,能计算30个震后时间点对应的震后形变。本研究针对Gao等的软件进行改进,可计算震后任意时间点对应的震后形变。新软件由3个部分组成:1)与32个震后时间点相对应的32套离散格林函数数值框架;2)格林函数插值计算程序,可针对上述32套格林函数数值框架进行插值运算,输出任意震后时间点对应的格林函数数值结果;3)积分计算程序,调用上述格林函数数值结果,计算任意类型地震在地表任意位置产生的同震与震后形变。一般情况下,使用者只需按要求准备辅助文件,提供发震断层模型和观测站位置信息,以及震中周围地区地幔粘滞性因子,先后运行格林函数插值计算程序和积分计算程序,即可计算出目标地震在地表任意位置产生的同震与震后形变。本文基于粘弹球体位错理论与弹性球体位错理论,分别计算2011年日本MW9.0地震引起的远场同震位移,2套结果的高度一致性证明了新程序的正确性。最后,介绍需要注意的若干事项,便于使用者掌握该软件。
鄂尔多斯地块的动力背景是青藏高原东北缘NE向推挤作用的结果。1970年以来,鄂尔多斯地块北缘现今地震活动表明,中强地震活动具有幕式活动特征,从1970~2015年共经历了4幕活动;块体西北缘与东北缘的中强地震具有迁移对跳活动特征,1970年以来共经历了4组对跳;震例统计结果揭示,地块北缘6.0级左右地震前,块体中等地震经历长时间平静具有普适特征;2005-03以来,块体所在的阴山地震带再次进入中等地震超长平静时段。其幕式活动特征、迁移特征、时间平静特征可作为区域范围强震的预测指标。
利用基于P波初动和振幅比反演震源机制解的HASH方法,得到鄂尔多斯块体周缘2008-01~2014-06共1 700个ML2.5以上地震的震源机制解,从中选取结果较为可靠的1 581个地震的震源机制解,将研究区域划分为1°×1°的网格,采用阻尼区域应力反演方法(MSATSI软件)获得了研究区域内中小地震震源机制解特征和应力场图像。鄂尔多斯块体周缘震源机制解类型以正断和走滑为主,应力场反演结果表明,鄂尔多斯周缘应力状态以拉张为主,局部区域存在剪切、挤压。块体北缘和西缘的最大水平主应力优势方位为 NE向,一致性较好;块体东缘和南缘的最大水平压应力方向与区域内主要活动构造方向基本平行,东缘最大水平主应力方向从南段的近EW向、NEE向逐渐过渡到NE向,南缘是NWW-SEE 和近 EW向。通过分析研究区域应力状态和最大水平主应力方向认为,银川-吉兰泰南端近NS向活动断裂的剪切应力状态的出现,是由于近NS向活动断裂在NE向挤压应力作用下发生了右旋剪切运动;块体西南缘NEE向的挤压和剪切的应力状态与青藏块体对鄂尔多斯块体的直接作用,以及六盘山断裂带第四纪以来强烈的挤压性质及左旋走滑运动有关;大同盆地出现局部剪切应力状态可能与盆地下方的火山有关。
利用1990~2016年水准观测资料,以GNSS点的垂向速率作为先验值约束平差,利用距水准路线上距离小于10 km的GNSS连续站的垂向速率检验,获取川滇地区的垂直运动速率场,可靠性优于0.87 mm/a,且误差分布相对均匀。研究表明,川西地区处于快速差异性隆升阶段,位于鲜水河南东段的贡嘎山地区隆升速率最快,约5~6 mm/a;四川盆地垂向运动不明显,约0~1 mm/a;滇西南地区以0~2 mm/a的速率下沉。同时,通过水准速率剖面分析了川滇地区主要断层的张压特性。
基于三维向错理论模拟倾角为60°的矩形断层在转动角不同量级和形式下引起的地表形变,分析断层转动对地表形变大小及空间分布的影响。利用向-位错模型,采用数值方法模拟鲜水河断裂带断层滑动与转动引起的地表形变,并与扣除趋势项的2009~2015年的GPS实测结果进行对比。研究表明:1)转动角大小与断层向错引起的地表位移场大小呈正比关系,位移场方向不变,转动方向不同时,地表位移场方向随之改变,大小不变,且断层附近地表位移场影响较大;2)基于向-位错组合模型对鲜水河断裂带模拟结果与GPS结果在大小和方向上均具有较好的一致性;3)鲜水河断裂带表现出一定的旋转特性,存在明显的分段差异性,具有一定的转换构造样式。
针对现有速度场模型采用数据少、现势性不强的问题,详细论述了速度场模型建立的数据获取与处理方法,利用中国大陆地壳运动监测网络中260个连续运行基准站2011~2015年的连续观测数据,采用GAMIT/GLOBK与GMT软件,获得高精度的定位及速度成果,并基于上述速度成果采用克里金插值法构建了中国大陆地壳运动速度场1°×1°格网模型。最后利用IERS公布的中国区域7个IGS站实测速度以及全国范围内均匀分布的陆态网络区域站实测速度,对模型内插结果进行检核,验证了模型的现势性、精确性和可靠性。
针对非等距沉降数据序列的建模问题,探讨了3种非等距GM(1,1)模型在沉降预测中的应用。结果显示,加权非等距GM(1,1)模型不适用于具有近似指数趋势的沉降数据,灰线性加权非等距GM(1,1)模型并非真正的非等距模型,模型预测式与时间无关,不能用于沉降预测,而对于笔者所建立的非等距GM(1,1)模型,其拟合函数等同于灰色线性回归组合模型,可用于近似指数趋势的沉降数据。通过实例分析对比了3种非等距GM(1,1)模型的应用效果,验证了上述观点的正确性。
通过对传统GM(1,1)缺陷分析和改进的基于权的PGM(1,1)建模机理描述,顾及PGM(1,1)中背景值构造时取相同的参数不能充分降低模型的预测误差,对不同的时刻引入不同的参数来改进GM(1,1)背景值序列的计算公式,将这种背景值构造方法和灰元N引入GM(1,1)建立了新的白化方程。在建立的新的白化方程基础上,用龙格-库塔法以修正的初始值计算累加值的模拟序列。针对引入的参数较多问题,采用粒子群算法寻找满足相对误差均值最优的参数,从而建立了基于粒子群优化算法和加权灰色组合的PSO-GM模型。工程实例应用表明,新模型的拟合精度高,预测效果好,相对其他两种原有模型预测精度有明显提高。
利用北京市GPS连续观测资料结合气象资料开展GPS水汽与气象要素的相关性分析。首先使用GAMIT软件解算2009-06-01~2012-04-30的北京GPS连续观测网观测数据并结合气压数据获得测站水汽序列;然后用小波变换方法对GPS水汽、温度和气压数据进行分解与重构,并对重构后的数据进行相关性分析。GPS水汽序列的变化趋势与温度呈显著正相关,与气压呈显著负相关。水汽和气压存在年周期和半年周期变化,两者存在显著负相关特性;温度有年周期变化,水汽和温度在d13重构结果的相关性最好,存在显著正相关特性。
轨道误差是InSAR数据处理中的一个重要因素,对从最初的图像配准到最后的高程值或形变值图像的生成都有着重要影响。含有误差的轨道参数造成基线误差以残差条纹的形式存在于干涉图中。本文利用4个条带的61景PALSAR升轨数据,通过InSAR干涉测量得到覆盖2011日本Tohoku地震的雷达视线向的地表同震位移场。通过联合GPS同震位移数据,采用二次多项式曲面拟合的方法去除InSAR同震形变数据中的卫星轨道误差,并分析雷达入射角变化对轨道误差去除的影响。对干涉图中拟合残差进行分析,改正后的InSAR同震形变场精度得到较大提高。
将无奇点根数和坐标旋转方法综合运用到地球静止轨道(geostationary earth orbit,GEO)卫星的18参数星历拟合,解决GEO卫星的小偏心率奇点和小倾角奇点问题。基于5颗北斗在轨GEO卫星1 a的轨道数据,分析全年星历拟合精度和稳定性,对卫星进出地影的拟合精度进行定量分析,针对Δn参数超限问题,比较绕X轴旋转5°和55°的Δn参数变化范围。结果表明,2~4 h的拟合弧长5颗GEO卫星的星历拟合精度最大为3 cm,迭代次数最多为5次;拟合弧长越长,拟合精度越低,迭代次数更少;当拟合弧段包含进出地影轨道时,拟合精度降低;旋转55°比旋转5°Δn参数的变化范围减少一个数量级,Δn参数全年无超限现象。
介绍了2013~2014年使用A10/028和FG5/214两台绝对重力仪在国家现代测绘基准体系基础设施建设一期工程中绝对重力观测的情况,并对观测结果进行分析。结果显示,两台绝对重力仪观测结果精度及较差都在5 μGal以内,满足高精度重力基准建设的技术要求;A10/028具有一定的野外观测适应性,FG5/214对重力细微变化信号的捕获具有优势,两台绝对重力仪同期观测对于提高观测结果的可靠性和准确度具有积极的作用。