用随机有限断层方法模拟2016-11-13新西兰Mw7.8地震

在未来大地震发震构造已知的前提下,提出用经验公式确定大地震断层长度、宽度、地震矩和滑动分布等参数,并用随机有限断层方法预测未来大地震。用随机有限断层方法模拟2016-11-13新西兰M_W7.8地震的12个基岩台站地震动时程和反应谱,并用模拟的频谱幅值与记录频谱幅值平均比值确定模拟误差。结果表明,在周期为0~10 s的范围内,模拟误差在0.92~1.08之间,不同频率模拟误差的标准差变化范围不超过1。95%的置信区间宽度随频率无明显变化,模拟结果反映了地震动记录的平均效果。用准随机方法重新获取新西兰大地震断层初始破裂点和位错滑动分布,模拟地震动并计算模拟误差随频率的变化,得到与最初模型相似的结论,从而进一步证实了本文提出的获取震源参数并用随机有限断层方法预测未来大地震的可靠性,且其特别适用于大地震远场模拟。     

基于共轭梯度搜索的病态问题处理方法

在最小二乘平差准则基础上,把病态平差问题转化为无约束的二次规划问题,并利用优化理论分析病态对平差解的影响。通过共轭梯度搜索算法在可行域中寻找最优步长因子,自动寻找最速下降方向,并给出迭代初值的设置方法。分析近似计算中病态问题与局部最优解的关系,讨论局部最优解的快速迭代方法,并通过实例验证算法的有效性,计算迭代的速度。由于整个过程没有对法方程系数矩阵进行求逆计算,该算法可用于处理大规模系数矩阵高病态的平差问题。     

基于全球分布的地震观测数据估计0S0模态的特征频率及品质因子

选取2011-03-11日本MW9.3地震后全球地震台站记录的约220 h波形数据,基于一定的判断准则,筛选出质量较高的448个垂向观测数据,采用频域AR方法对0S0模态的特征频率和品质因子进行估值,并以信噪比为权重进行加权平均,得到特征频率f=0.814 658 4±5.3×10-7,品质因子Q=5 586±12。该结果与已有研究成果相比精度更高,有助于约束现有地球介质密度和衰减模型,使其更接近真实地球。     

青藏高原东南缘构造应力场重构

采用有限元方法,针对青藏高原东南缘建立更细致、更精确的三维有限元弹性模型。选取9种不同的应力边界条件,分别进行优化分析后处理,将对应台站形变模拟值与GPS实测值进行误差分析,最终选取最佳方案作为古构造应力场。结果表明,青藏高原东南缘4 Ma BP的古应力场主要起源于中国大陆周围板块的相互作用,特别是印度板块NNE向强烈碰撞作用,成为中国大陆尤其是西南部青藏高原地区构造应力场最主要的动力来源,控制各个块体相互作用的方式和运动格局。青藏高原东南缘古应力场主要包括几个力源：西北部青藏高原侧向挤压造成的WE向应力约105 MPa;西南部直接来自于印度板块的NE-WS向应力约70 MPa;南部NS向作用力33 MPa;东南部扬子块体侧向NW-SE阻挡力56 MPa;北东部受扬子块体强烈EW向阻挡力90 MPa。这些力源共同作用于青藏高原东南缘,形成现今复杂应力场。
     

基于TCP-InSAR的采动区铁路动态沉降监测与时序分析

为高效获取采动区长时间序列形变,监控煤炭开采对矿区铁路的影响,研究一种基于TCP-InSAR(temporarily coherent point interferometric synthetic aperture radar)的采动区铁路形变监测方法。该方法根据时间序列SAR影像间的相干性,选取临时相干点构建Delaunay三角网,并通过离群值探测去除具有相位模糊度的TCP间的弧段,最后采用最小二乘解算得到区域地表变形。实验使用2016-10～2017-04时间段内的15景Sentinel-1A数据,利用TCP-InSAR技术获得某矿区的铁路形变。结果表明,受采动影响,该时间段内铁路最大下沉值为95mm,最大倾斜坡度为0.37‰。利用TCP-InSAR技术可实现采动区内铁路长时间的动态形变监测。     

GRACE时变信号恢复方法对比分析

针对利用GRACE重力数据反演地表质量变化过程中产生的时变信号削弱和泄露现象,分别采用尺度因子法和迭代恢复法对信号进行恢复,并选用CSR Mascon产品从时间序列变化和空间分布两个角度对两种恢复方法的效果进行对比分析。结果表明:1)尺度因子法虽实现上较为简单,但其计算结果完全依赖于所选的先验模型,当所选模型的可靠性较差时,会影响到GRACE反演结果的可靠性。2)迭代恢复法的恢复过程仅受GRACE原始观测值的约束,能较好地恢复研究区域的长期趋势信号和周期信号,但在部分区域可能存在信号过量恢复的现象。建议采用迭代恢复法进行GRACE时变信号恢复。     

非等间距GM(1,1)模型的总体最小二乘算法及其病态问题

在非等间距GM(1,1)模型中,系数矩阵中有无误差的常数项和有误差的随机项,并且系数矩阵与观测向量误差同源,即系数矩阵与观测向量中有相同的元素存在,这些相同元素应该有相同的改正数,为此本文推导了一种适合非等间距GM(1,1)模型求解的总体最小二乘算法。同时,考虑到非等间距GM(1,1)模型中存在病态问题时影响总体最小二乘计算结果的稳定性,提出对系数矩阵常数列乘以某一常数的方法,以改善病态问题。     

大连长兴岛海域疏浚区GF-1 WFV影像悬浮物遥感反演研究

利用2017年8-9月期间大连长兴岛海域实测遥感反射率、悬浮物浓度及GF-1 WFV数据,研究了该区域悬浮物浓度的遥感反演算法,发现遥感反射率与悬浮物浓度的散点图存在两种不同的变化趋势,将2008-2015年期间辽东湾内其它海域的数据与长兴岛海域数据进行叠加后,可初步推断产生两种变化趋势的主要原因在于疏浚区与非疏浚区悬浮物的粒径大小不同,导致其光学特性和反演算法也存在差异。将两种算法分别应用于2017年9月20日辽东湾及长兴岛近岸海域的GF-1 WFV4影像,发现采用非疏浚区算法反演长兴岛海域悬浮物浓度时,出现明显的低估现象,而疏浚区算法则能较好地反映出疏浚施工过程中悬浮物的分布规律,具有较好的工程应用价值。     

深海滑坡灾害的物质点法模拟

深海滑坡可以沿着斜坡搬运大量沉积物,是危害巨大的海洋地质灾害之一。评估深海滑坡的流滑机制并量化分析其对海底结构物的冲击对海洋资源的开发利用至关重要。采用物质点法对滑坡的滑动特性进行了总应力分析。设计了三组不同的土体、海床参数,展现了三种不同的滑动形态:延展、块体和扩张。随后复现了流滑体对部分暴露管线的冲击过程,并对冲击力的稳定值进行了分析,分析中考虑了流滑体的惯性、抗剪强度和静压力的影响。研究表明,物质点法能够用于模拟深海滑坡的流滑过程及其对海底结构物的冲击,可以为海洋工程中的实际设计工作提供服务。