刚性挡土墙后三维被动滑裂面的模型试验

针对无黏性土体,采用模型试验的方法,研究平移模式下刚性挡土墙后被动破坏滑裂面的空间形态。自主设计一种模型试验装置,重复开展6次试验,通过记录挡土墙后土体中预埋脆性玻璃条断裂的空间坐标,复原土体发生滑动的位置,绘制出挡土墙后滑裂面的三维形态图。试验结果表明：挡土墙后滑裂面具有明显的三维效应;挡土墙宽度内滑裂面纵向高度呈先缓慢增高后近似直线增高的曲面,初始破裂角度为9°,平均破裂角为26°,朗肯土压力理论的破裂角为28°;最大纵向破裂面长度为1.8倍挡土墙高度,与经典土压力理论的平面假定基本一致;滑裂面均有一定的横向扩展,主平面投影以初始扩散角约45°的斜线往外扩展,距离挡土墙最远处是宽度为0.7倍挡土墙宽度的水平线,斜线与水平线之间以半径为挡土墙宽度的圆弧过渡连接。研究结果为分析土体被动破坏的滑裂面空间形态提供了试验依据。     

Numerical investigation of the uplift performance of prestressed fiber-reinforced polymer floating piles

Finite element analyses of prestressed fiber-reinforced polymer floating piles subjected to uplift force have been conducted in this paper. First, parameters of the modified BPE model (bond–slip model at the fiber-reinforced polymer–concrete interface) were calibrated using existing pullout testing data on fiber-reinforced polymer rebars embedded in concrete. Nonlinear spring elements were used in numerical modeling to characterize the bond–slip behavior at the fiber-reinforced polymer–concrete interface. A parametric study was performed to assess the influence of rebar diameter, fiber-reinforced polymer material, embedment length, and concrete strength on the mobilized bond stress. Upon the successful modeling of the pullout performance of fiber-reinforced polymer rebars in concrete, numerical models were developed to investigate the dependence of the uplift performance of floating piles on the prestress level, uplift force, fiber-reinforced polymer type, and compressive strength of concrete.     

基于InSAR和GPS观测数据的尼泊尔地震发震断层特征参数联合反演研究

利用日本ALOS-2和欧空局Sentinel-1A卫星获得的尼泊尔地震同震形变场,结合GPS同震位移数据,联合反演了断层滑动分布特征和空间展布.结果表明:尼泊尔地震的同震形变场主要集中在150km×100km的范围内,且分为南北两个相邻的形变中心,南形变中心的视线向抬升量约为1.2m,北形变中心的视线向沉降量约为0.8m,均位于发震断层上盘.位于形变抬升区的KKN4和NAST两个GPS站,抬升量和南向运动量均达到了m级,而远离震区的其他GPS台水平和垂直观测量均在1cm以内.联合反演得到的断层位错分布主要集中在沿走向150km,沿倾向70km的范围内,最大滑动量为5.59m,平均滑动量为0.94m.断层面倾角在浅部约为7°,随着深度增加,倾角逐渐变大,到垂直深度20km时倾角接近12°;5月12日MW7.2级余震位于主震破裂区的"凹"型滑动缺损区域;主震破裂区的上边界与MBT空间位置十分吻合,主震破裂区主要集中的MBT以北50~60km处,垂直深度为8~9km,倾角为9°,继续向北时主震破裂面以10°~12°的倾角向深延伸,在18~20km可能与MHT交汇.因此,初步判定MBT为此次地震的发震断层.     

近断层应力变化快速计算方法及对强余震空间分布的指示意义:以M_w7.9汶川地震为例

构造地震一般由断层摩擦失稳所致.断层内部及周边所累积的剪切形变则通过同震滑动位移得到局部释放.因此,震后断层内部及近断层周边的静态剪切应力变化量的空间分布可通过断层面上的滑动位移分布计算得到.本文采用傅氏变换方法(FTM)计算单一有限断层同震滑移场所形成的静态剪切应力变化分布,近断层两侧的应力变化计算可由波数域内应力近似算法获得.结果表明,FTM快速有效、易于实现,有效地避免了常规应力计算中奇异值的出现.以2008年Mw7.9中国汶川大地震为例,采用前人所得有限断层滑动位移模型,得到了断层面和近断层周边准3D剪切应力分布解,并同主震后中强余震的空间分布特征作了比较.结果表明,大部分的中强余震震源位置处于剪切应力变化值为正的区域,由同震位移所产生的局部应力降峰值和均值大小同应力变化的正值大致相当,从而表明了快速且有效地计算断层内部及近断层附近的应力变化分布可以为主震后强余震发生的潜在区域提供指示意义.需要强调的是,应力变化空间分布特征的获取强烈地依赖于有限断层滑移模型解.有关滑动位移反演解的多解性对应力变化计算结果的影响,本文作了必要的讨论.     

云南鲁甸M_S6.5地震红石岩滑坡稳定性的数值模拟

2014年8月3日鲁甸MS6.5地震触发了大量的滑坡崩塌,其中,位于鲁甸县李家山村和巧家县红石岩村交界处的牛栏江干流北岸的红石岩滑坡规模巨大,与此处位于左岸的红石岩古滑坡体的前缘部分一起堵塞了牛栏江而形成高达120m、体积达1 200×104m3的大型堰塞体。通过震后开展的野外实地调查,获得了红石岩滑坡发生处的地形地貌、地质构造、岩体结构及物质组成等资料。以这些第一手资料为基础,构建了红石岩滑坡的边坡模型,并应用边坡稳定性分析软件Geo Studio中Slope/W模块分别计算了红石岩滑坡体震前坡体安全系数和地震作用下的坡体安全系数。结果表明,红石岩滑坡体发生处的坡体安全系数在地震前为1.450,处于相对稳定状态,而鲁甸地震的地震动作用则使坡体的安全系数降低至0.962,直接导致红石岩坡体的失稳。文中进一步讨论了坡体滑动面的存在与否对坡体稳定性的影响:安全系数计算的结果表明,在中强地震作用下,先存滑动面的存在是导致大型滑坡形成的重要条件;对于高陡岩质边坡,如果没有先存滑动面,只可能形成浅表性滑坡。     

向错-位错组合模型模拟渭河盆地地壳水平变形特征

基于向错-位错组合模型模拟渭河盆地内主要断层的滑动与转动运动变形,通过数值计算获得断层滑动及转动变形引起的地表水平位移,并与盆地内实测GPS水平位移进行了对比。结果表明:采用向错-位错组合模型能够完整地描述断层实际的滑动与转动运动变形状态,组合模型模拟断层滑动与转动引起的地表水平位移与实测GPS水平位移,无论在量级或是运动趋势上均具有较好的一致性特征。     

附加电离层延迟约束的实时动态PPP快速重新初始化方法

传统模式下,由于卫星信号中断等原因导致的实时动态PPP重新初始化问题将严重影响PPP的实际应用效果.针对上述问题,本文提出了一种有效的解决思路,首先设计了电离层延迟变化预报模型,并提出预报信息的定权方法,进而将上述预报结果作为虚拟观测值,联合历元间/星间差分观测数据估计得到周跳参数浮点解,并尝试固定,从而实现PPP快速重新初始化.该方法通过实现电离层预报结果权重的自适应确定,改善了参数估计结果的精度,提高了周跳修复成功率及可靠性.利用不同区域和不同时期的大量观测数据分析表明,周跳修复效果与测站位置、太阳活动及中断时长存在密切关系,当中断时间较短（90 s以内）时,均可以获得较高的成功率,平均达到90%以上,由此验证了本方法的有效性;随着中断时长进一步增大,不同地区和太阳活动情形下的修复效果存在显著差异,纬度越低,太阳活动越剧烈,成功率越低.     

汶川地震同震形变场的GPS和InSAR邻轨平滑校正与断层滑移精化反演

为克服InSAR观测汶川地震同震形变场的邻轨不连续问题,提出联合GPS观测值与邻轨平滑约束的同震位移校正方法,采用GPS观测形变去除PALSAR轨道误差引入的残留平地相位,基于形变平滑条件校正邻轨干涉相位的不连续性.ALOS/PALSAR干涉处理结果表明,校正后同震形变场的准确度与平滑性得以显著提高,InSAR高相干点残差达3.6 cm,校正后精度提高约60%,低相干点精度提高约40%,校正后形变场的邻接平滑因子标准差减小约33%,验证了轨道误差校正与邻轨平滑约束方法的准确性与可靠性.进一步基于弹性半空间位错模型的断层滑移反演结果表明,断层滑移主要分布于映秀、北川和青川地区,集中于地壳深度0~16 km范围,最大滑动量（位于北川县城）约为9.0 m,GPS反演模型残差为5.5 cm,InSAR反演模型残差达9.2 cm,InSAR反演精度约有30%的显著提高,由模型反演计算得到的地震矩为8.0469×10²⁰ N·m.     

侧向侵蚀相关的走滑断裂滑动速率计算新方法

断层滑动速率是活动构造研究中的重要内容,是反映断裂活动性和地震危险性的重要参数之一。随着测年技术不断发展和测年精度大幅度提高,全新世甚至千年尺度和百年尺度的年轻地质体的位错也越来越多地被用于断层滑动速率计算。用走滑断裂带上地质体实测年龄计算滑动速率,会受到2种因素影响:1)累积位移时间是否与所测地质体年代相符合;2)地质体位移形成过程中会受到侵蚀。在利用全新世地质体计算断层滑动速率时,应将侧向侵蚀的影响剔除。因此,文中提出1种计算走滑断层滑动速率的新方法——差值法。走滑断层上河流阶地演化与断层位错分析表明,在阶地拔河高度存在较大差异的情况下,可以利用阶地拔河高度与年龄按比例进行计算。此方法在一定程度上提高了所计算滑动速率的精度,但是需要至少有3级不同阶地的拔河高度、年龄以及位错信息。若阶地拔河高度近似呈等差排列,即各级阶地上侧向侵蚀量近似相等的情况下,利用高-低阶地累积位错量之差与对应阶地年龄差来计算滑动速率,可以在一定程度上减少上述2种因素对滑动速率的影响。应用差值法计算得到阿尔金与昆仑断裂的滑动速率为4.7~8.8mm/a,与前人获得的地质学滑动速率、测地学滑动速率、缩短速率以及强震复发周期结果一致。     

Mechanism on Moho offset induced by aseismic slip of deeply buried faults

Recent high-resolution deep seismic reflection profile across the Kunlun fault in northeastern Tibet shows clearly that the Moho is cut off by a complex thrust fault system. Moho offset is a general phenomenon, but little is known about the dynamic mechanism. In this study, contact models with Maxwell materials are used to simulate the mechanical process of Moho offset induced by the aseismic slip of deeply buried faults. Based on the seismic reflection data, we project a single fault model and a complex fault system model with two faults intersecting. The deformations of the Moho, the aseismic slips, and contact stresses on faults in different models are discussed in detail. Results show that the Moho offset might be produced by aseismic slip of deeply buried faults, and the magnitude is influenced by the friction coefficient of faults and the viscosity of the lower crust. The maximum slip occurs near the Moho on the single fault or at the crossing point of two intersecting faults system. Stress concentrates mainly on the Moho, the deep end of faults, or the crossing point. This study will throw light on understanding the mechanism of Moho offset and aseismic slip of deeply buried faults. The results of complex fault system with two faults intersecting are also useful to understand the shallow intersecting faults that may cause earthquakes.