稀疏台网的中小地震震源机制解反演研究——以2017年1月4日西藏仲巴4.7级地震为例

本文应用ISOLA近震全波形方法,以2017年1月4日西藏仲巴4.7级地震为例,反演稀疏台网记录的中小地震震源机制解。该地震反演所得最佳双偶机制参数为：节面Ⅰ的走向109°/倾角85°/滑动角-177°,节面Ⅱ的走向19°/倾角88°/滑动角-4°,最佳矩心位置为30.590°N、83.784°E,最佳质心深度为6km,矩震级M_W4.6。震源机制反演结果表明此次地震是一次走滑型为主的事件,其与震源区域附近历史地震震源机制解具有相同性质。本文还应用CSPS初动扫描法,利用P波初动资料和近震波形联合约束反演此次地震的震源机制,并与ISOLA近震全波形反演结果进行比较,结果表明,联合少量台站的的三分量波形数据,能够定量地判断最佳震源机制解,降低了P波初动反演结果的非唯一性,同时也约束了由于少量台站参与全波形反演引起的解的不稳定性。本文研究为中小地震震源参数测定提供了一种简单有效的方法,具有较高的稳定性和可靠性。     

2019年6月17日四川长宁6.0级地震中心震源机制解及震源区构造应力场研究

不同资料和方法给出的2019年6月17日四川长宁6.0级地震震源机制解存在较大差异,为了找到1个合适的震源机制解来研究此次地震的发震方式,通过数学方法得到了与现有震源机制解差别最小的中心震源机制解,节面I的走向、倾角、滑动角分别为194.78°、52.68°和139.16°,节面Ⅱ的走向、倾角、滑动角分别为312.44°、58.67°和45.22°,根据本次地震余震分布拟合得到的断层面的走向为312.17°,与中心震源机制的节面Ⅱ走向一致,因而推断节面Ⅱ为本次地震的发震断层面。之后,利用此次地震之前震源区地震的震源机制解,反演了震源区的震前构造应力场。结果表明,长宁6.0级地震的中心震源机制解和震源区震前应力场均为逆冲型为主兼走滑分量的类型,震前应力场压轴为NWW—SEE向,中间轴为NNE—SSW向,两轴倾角接近水平,而张轴较陡,表现为逆冲型的应力场。将反演得到的应力场投影到中心震源机制解给出的与余震分布一致的节面上,发现中心震源机制解的滑动角和应力场预测的滑动角差别仅为13.45°,表明此次地震受背景应力场控制而发生在先存的薄弱面上。     

九寨沟M_S7.0地震震源机制解及构造应力场研究

2017年8月8日四川九寨沟县发生M_S7.0地震。根据中国地震台网固定台站的记录波形,利用gCAP方法和P波初动符号反演方法求解了主震的震源机制解,其结果与哈佛CMT震源机制解、美国地质勘探局发布的震源机制以及中国地震局发布的各震源机制解有很好的一致性。此外,还求解了九寨沟地震震中附近2010—2016年间28次小震的震源机制解,利用得到的小震震源机制解反演了该区域的构造应力场。结果表明,这一区域构造应力场的最大应力轴与最小应力轴均为水平方向,其中最大应力轴方向为NWW方向,最小应力轴方向为SSW方向。     

重庆垫江MS4.4地震:一个可能发生在浅地表的地震

为更好地认识2016年8月11日重庆垫江M_S4.4地震的发震成因,本文首先分析了近场地震台站记录到的该地震产生的短周期面波（Rg）强频散现象,该现象表明本次地震具有明显的浅震特征。然后采用CAP方法反演了该地震的震源机制解和震源深度,结果显示本次地震是走向为近EW向的逆冲型地震,矩震级为M_W4.1,最佳震源深度为1—2 km。在此基础上,结合现场震害、烈度等震线长轴方位和地质剖面图等资料,分析认为在区域应力集中加载作用下,震源区沉积盖层内部的软弱地层产生挤压变形,导致地层内斜交的裂隙发生切层失稳滑动,从而诱发了垫江M_S4.4地震,同时排除了震源区NE走向的黄钦垭口地表断层为发震构造的可能。      

青藏高原区域构造应力场特征分析

利用区域测震台网的地震波形及震相数据基于CAP方法反演北纬26°~42°,东经90°~110°内的270个M_S>3.0以上地震的震源机制解;并结合GCMT目录和一些前人研究结果中该区域的共759个震源机制解数据,运用SATSI方法计算了研究区域的应力场。将研究区域按1°×1°网格化后得到了154个局部应力场分布结果,从结果上看,整个青藏高原的最大主压应力方向大致表现为顺时针且向右的旋转方式。该结果反映了青藏高原块体向NE和NNE挤压的过程中,分别在其北部和东部受到鄂尔多斯和阿拉善两个坚硬块体的阻挡,造成青藏块体增厚,块体之间物质的侧向流动。对研究区域应力型因子R值的研究显示青海祁连、甘东南区域、四川龙门山断裂带等沿青藏高原块体与阿拉善块体以及鄂尔多斯块体交界处相对应力值偏大,与近年来这些地区的地震活动性成正比。本文研究结果对比其他应力场研究结果、GPS研究结果、以及数值模拟结果具有很好的一致性,可为青藏高原地区的孕震机理、活动构造以及地震趋势判定提供可靠的参考依据。     

使用L0范数代价函数的二维相位快速解缠方法

针对最小范数解缠方法的代价函数模型以及迭代效率问题,提出一种符合L0范数准则的高效的二维相位全局最优解缠方法。在对最小范数代价函数模型特征进行分析的基础上,给出一种符合L0范数准则的代价函数,满足在相位不连续边界的切方向上相对法方向具有更强的约束作用,保证解缠迭代处理时连续相位按照代价函数的约束条件进行趋近的同时,不连续边界不会被破坏。为了解决线性方程中低频误差收敛慢的问题,根据最小范数方法独立性的特点,采用数据分块处理方法,注重高频信息以提高迭代处理效率,将低频信息处理转移到数据块间偏移量计算中。实验对比分析表明,新解缠方法在可靠性和效率方面都能够达到较好的水平。     

确定似大地水准面的Hotine-Helmert边值解算模型

空间大地测量技术的发展使大地高的观测成为可能,从而为第二大地边值问题的研究带来了新的机遇,本文对基于Helmert第二压缩法的第二边值问题(简称为Hotine-Helmert边值问题)展开研究。首先介绍了地形直接、间接影响的定义与算法,然后推导了Hotine-Helmert边值问题的解算模型。Hotine-Helmert边值理论无须计算地形压缩对重力的次要间接影响,因而较Stokes-Helmert边值理论更简单。此外,文中引入了一种低阶修正的Hotine截断核函数,该核函数较传统的截断核函数能有效地改善似大地水准面的解算精度。为了验证本文构建的Hotine-Helmert边值解算模型的有效性和实用性,本文将EIGEN-6C4模型的前360阶作为参考模型,利用Hotine-Helmert边值解算模型构建了我国中部地区6°×4°范围、1.5′×1.5′分辨率的重力似大地水准面,其精度达到±4.8 cm。     

基于多源InSAR数据的三维地表形变解算方法研究

随着不同成像参数SAR卫星不断发射升空以及新的InSAR技术不断出现,同一地区可以获取多尺度、多源、异质的InSAR地表形变监测数据。将这些多源InSAR地表形变监测数据作为观测量,实际三维地表形变作为未知参数,建立模型并解算,实现东西、南北和垂直向地表形变监测,为地质灾害预报及监测提供更全面的基础资料。本文以InSAR三维地表形变解算为核心研究内容,针对大梯度、轨道和地形、粗差等因素对解算精度的影响,利用现代测量数据处理理论和方法,消除或抑制上述因素导致的各类噪声,实现多源InSAR数据的三维地表形变高精度解算。     

CE-5器间分离的同波束干涉测量的准实时监测法

“嫦娥5号”探测器的组成包括轨道器、返回器、着陆器和上升器。在环月阶段,两个组合体(轨道器/返回器组合体和着陆器/上升器组合体)之间的分离实时监测,是飞行控制的关键检测段。本文提出利用甚长基线干涉测量(VLBI)测轨技术实现绕月探测器器间分离实时监测。特别是在器间分离前后,本文方法能够利用探测器器间两路下行信号进行同波束干涉测量(SBI),差分时延测量能够提高器间相对距离的解算精度。单基线试验分析表明,本文方法基于实测数据的“嫦娥3号”着陆器两天线的相对距离解算精度优于0.3m,平均相对距离误差约为0.15m,基于“嫦娥5号”仿真数据的双阈值判定对器间分离监测的响应时延优于30s。     

北京新机场控制点测量与净空障碍物普查工作

为保障北京市新机场建设和制定后期安全运行措施,对新机场进行控制点测量与净空障碍物普查工作。本文首先介绍了采用GNSS静态测量的方式对埋石布设的79个控制点进行观测的方案及数据处理方法,并对控制网处理的精度进行了分析。然后介绍了对机场周围障碍物进行了平面位置和高度测量的依据和方法。最终成果资料精度可靠,可供新机场建设使用。