芦山地震余震发生及演化的自组织临界机制

为探讨芦山7.0级地震余震演化动力机制,应用统计地震学方法,分析了2013年4月20日—6月20日芦山地震余震序列的宏观统计分布规律;基于自组织临界理论,提出了一种新的余震模型,以期阐明相关统计地震学规律的产生动力机制,并深入讨论了该模型的自组织临界性.该余震模型的具体算法是在经典Olami-Feder-Christensen地震模型基础上,引入了应力衰减因子和应力扩散各向异性因子.研究结果表明:芦山7.0级地震余震序列的震级分布遵循Gutenberg-Richter统计规律,幂指数值约为0.766;其余震序列的时间分布遵循Omori统计规律,幂指数值约为2.52.新建模型的数值模拟能同时对芦山地震余震序列呈现出的Gutenberg-Richter和Omori统计规律给出满意的预测结果,模拟结果与实际情况高度吻合,表明龙门山断裂带处于一种自组织临界状态,芦山地震余震过程实质上是一种自组织临界现象.      

海地强震 全球救援

我国捐助财物800多万美元 当地道路、机场、港口等交通设施受损严重 救援受阻 加勒比岛国海地当地时间1月12日16时53分（北京时间1月13日5时53分），海地发生里氏7．3级地震，震中位于太子港西南约16公里处，震源深度约10公里。大地震发生不久后又连续发生多次余震。这是海地自1984年发生里氏6．7级地震以来最大的地震，也是加勒比地区数百年来遭遇的最强地震。     

2013年4月20日芦山地震最大加速度记录分析

2013年4月20日芦山地震是汶川地震后龙门山断裂带发生的又一次破坏性地震,中国数字强震动观测网络有100余台仪器触发并获得加速度记录.设置在四川省宝兴县防震减灾局的强震台(编号51BXD)在东西向获得了1.026g(-1 005.3gal)的峰值加速度记录,这是我国首次在自由地表获得峰值超过1g的记录,是我国强震动观测历史上的重要事件.经震后灾害调查,发现观测台站附近建筑物破坏并不严重,确定地震宏观烈度为Ⅶ度,这引起了社会的广泛关注与争议.通过现场考察,初步认定台站的特殊选址带来了局部地形影响,陡峭山坡放大了高频段地震动.还分析了宝兴县另外3个强震台站的观测环境,并采用全球新一代地震动预测方程进行对比分析,结果显示,该台在小于0.2 s范围内的反应谱放大明显.为了进一步说明台址受场地地形影响,对台站附近进行了背景噪声测试,采用水平/垂直谱比法计算了场地背景噪声的卓越周期、水平向放大系数,得出了随着地势由山脚向山坡上增高,场地卓越频率逐渐增大,且放大系数也随之增大的结论,这与51BXD强震动记录反应谱特征一致.最后,与汶川地震典型记录反应谱、抗震规范谱等进行对比,分析了51BXD强震动记录不同频段的特点;结合附近建筑物现场震后调查结果,初步解释了造成建筑物破坏较轻的原因,验证了局部地形影响是造成峰值加速度记录较大的主要原因.      

基于自适应匹配算法的速度脉冲型地震动模型化

速度脉冲型地震记录至今尚未有完整的波形参数识别方法，从地震动数据库中批量识别脉冲记录并搜集波形参数仍然有一定困难. 针对此问题，将脉冲量化识别、脉冲解析函数、自适应匹配算法三者结合，提出一套脉冲波形参数识别方法，可实现脉冲波形和解析函数之间的自适应匹配；该方法包括一个借助EMD (empirical mode decomposition) 的核心算法，对近场速度脉冲型地震记录实现平滑；再以Mavroeidis函数为脉冲基本模型，根据其特征参数的变化范围，分别对脉冲频率、脉冲形状系数、脉冲振荡特征系数进行了参数样本的循环枚举，并由优化算法确定模型参数取值；在此基础上，给出了当前既有的速度脉冲型记录样本的上述模型参数识别结果，并根据统计给出了地面峰值加速度分布和主脉冲周期的关系. 研究结果表明:和既有的方法利用Dickinson模型识别方法相比，本文算法对波形显著不对称的记录主脉冲周期识别错误率降低了约67%，较好的识别稳定性更有利于工程应用.      

芦山地震触发崩塌滑坡的优势方向与机理

利用芦山地震强震记录、水系图和1 754个地震触发崩塌滑坡数据,对芦山地震触发崩塌滑坡的优势方向的控制因素进行了研究.对比分析了崩塌滑坡数的优势方向与研究区自然坡体的优势坡向;基于对抗性原理,确定了芦山震区新构造应力场方向;对崩塌滑坡确定性系数分区进行了分析.认为自然坡体的优势方向是决定崩塌滑坡方向的首要因素,但应考虑不对称地形造成的坡向分异现象;新构造应力场是崩塌滑坡方向的全局性影响因素,地震台站的最大累计加速度方向是其代表区域内崩塌滑坡方向的影响因素.对于盲断层型地震,上述三因素对地震触发崩塌滑坡方向的控制效应可能是具有普适性意义的规律.      

基于随机地震动模型的黄土隧道动力响应研究

汶川地震震损表明,隧道结构在强震作用下会出现结构损伤和破坏.地震动的发生具有很强的随机性,而目前大多采用经典地震动记录对进行隧道地震动力时程分析,无法保证与震源机制和场地特性相匹配,采用确定性地震动记录,对地震动随机特性考虑尚不够完善.因此,基于杜修力地震动模型,采用脉冲函数法求解动力方程,对某黄土隧道进行了随机响应分析分析,研究了随机地震动作用下黄土隧道随机地震响应特征.计算结果表明:隧道拱顶的主应力功率谱密度函数呈单峰,其响应主要发生在低频,而拱腰、边墙、墙脚、轨道板主应力功率谱密度函数呈双峰,主要发生在低频和中频;隧道位移响应主要在发生在地震方向,其响应均方差在2.1 cm以上,轨道板位移响应很小;隧道拱腰和墙脚的主应力响应很强烈,其均方差在3.1 MPa以上;拱腰和墙脚为黄土隧道地震响应最不利的位置.该研究可以为黄土隧道抗震理论提供支撑.     

加筋土挡墙地震稳定性破裂面随机搜索法

为确定加筋土挡墙破裂面的位置，基于水平条分法提出了一种多线段破裂面表现形式. 将加筋土挡墙破裂面视为由多条线段组成，各条线段在一个平面内以不同的长度和角度相互连接；根据地震作用下水平土条的力学平衡条件，推导出与破裂面参数相关的筋材拉力计算式；将筋材总拉力作为目标函数，采用两层循环方式进行求解计算，外层为墙后填土水平表面破裂点位置循环，内层为随机角度循环；对比每次外层循环计算所得筋材总拉力，取最大值所对应的破裂面为所求加筋土挡墙临界破裂面. 通过算例对比验证了多线段破裂面计算方法的合理性，并对加筋土挡墙稳定性影响因素进行分析. 研究结果表明:以角度随机方式产生多线段破裂面的计算方法无需进行数学优化可得到合理结果，加筋土挡墙破裂面位置比对数螺旋破裂面更接近临空面；填土内摩擦角的增大使得筋材总拉力与筋材长度减小，能够增强加筋土挡墙的内部稳定性.      

地震动非参数化谱反演可靠性分析

非参数化谱反演是研究地震动震源、传播路径及场地效应的关键技术,反演的可靠性对地震震源物理的准确认识及地震危险性的有效预测具有重要影响.?为分析震源和场地项权衡的约束条件、单步法/两步法的选择、震源参数估计等关键环节对反演可靠性的影响,以2016—2017年意大利中部地震序列的地震动谱反演为例,给出确保反演可靠性的合理建...     

康定Ms5.8级地震冷竹关坡体内地震动响应特征

为了揭示坡体内不同水平深度峰值加速度变化特征及地震触发大型滑坡形成机理,在冷竹关沟两岸安置强震监测台阵,对斜坡内部不同深度地震动响应进行监测研究.2014年11月25日康定Ms5.8级地震触发了位于坡体内部不同深度的五台强震监测仪,数据揭示:自坡体表面水平向内,各监测点水平向峰值加速度逐渐减小,且0~45 m的表层向内峰值加速度下降幅度较大,在坡体内部的下降幅度变小,洞深99 m处的水平向峰值加速度约为洞口监测仪的60%;水平向加速度反应谱也表明自坡面向内各监测点加速度幅值逐渐减小,标准反应谱的动力放大系数都小于3.5;傅里叶频谱也表明自坡面向内各监测点幅值逐渐减小,且越靠近洞口傅氏谱频率成分越复杂.      

“斗山一家人”抗震救灾突击队再现玉树

2010年4月14日7点49分,7.1级强震,再一次撕裂中华民族刚刚愈合的、源自”5.12”汶川特大地震的深创巨痛。青海省玉树藏族自治州玉树县发生7.1级地震,地震波及范围3万平方公里,12个乡镇受灾。截止19日上午8时,玉树地震遇难人数为1944人,失踪216人,受伤12135人,其中重伤1434人。     

GPS资料在2015尼泊尔Mw7.8级地震中的反演应用

为了探索近场高频GPS资料在重特大地震实时发震特征反演方面的潜在作用和优势,独立使用2015年尼泊尔Mw7.8级地震震中附近的16个高频（1 Hz）GPS站点数据,基于动态数据处理技术、矩张量反演技术及跟踪方差技术,模拟实现了对该次地震的触发、实时矩心/矩张量的反演及地震事件结束判断等相关工作;结合GCMT、W-相位以及部分已公开发表成果,对反演过程及反演结果中存在的问题进行了对比分析;初步探讨了GPS同震位移场产品在节面判定中的潜在补充作用. 研究结果表明:本次地震破裂的传播方向主要为自西向东南单侧传播,平均传播速度约为3.1 km/s,地震持续时间至少为75 s;GPS数据反演得到的实时矩心及矩震级结果与GCMT、W-相位结果一致,且反演序列具有较高的稳定性和可靠性.      

THE DESIGN OF DYNAMIC SIMULATION SYSTEM ON EARTHQUAKE SURROUNDINGS

Every year ,manyearthquakesoccurallovertheworld .Earthquakealwaysbroughthumanbe ingsterribledisasterandamassofpropertyloss ,anditwillbemoredangerousalongwithmoreandmoreskyscrapers ,large scalebridgesandhugereservoirs.Becausealmostallearthquakesimulationdevicesareusedtoinvestigation ,peoplewhohavenot gonethroughtheearthquakedisasterknownlittletoit.Thispaperbringsoutanearthquakesurroundingssimulationsystembywhich peoplecanbefamiliarwithearthquake phenomena ,includeearthquakeomen ,strongshockanda…     

GPS资料获取的四川地区应变率演化状态

为了研究汶川5·12地震事件对青藏高原东南缘及四川盆地的应力再分配作用及该区域现今的地震危险性,使用1997-2015年共计355个GPS站点观测资料获取了该区域1997-2008、2009-2011、2011-2013、2013-2015年共计4期独立的速度场,基于规则格网和张力样条插值方法计算了上述速度结果对应的面膨胀率、垂直轴旋转率及最大剪应变率.以2008年汶川Ms8.0级地震为主要节点,以芦山2013年Ms7.0级地震、康定2014年Ms6.3级地震为参考节点,分析了应变率场的时空演化过程及其与地震事件的潜在关联性.研究结果显示:随着汶川震后的时间推移,沿龙门山断裂的3种应变率分布均有减弱趋势,但至今仍未恢复至震前状态,说明该区域的震后回弹受下地壳黏弹性流变物质的弛豫影响较为显著;安宁河-则木河断裂及周边的面压缩率在汶川地震之后有增强趋势,该趋势可能与区域的地震危险性有所关联;安宁河-则木河断裂带及鲜水河断裂北段是未来一段时间该研究区域内应该持续关注的主要地震危险区.      

荆岳长江大桥附近大风特征及设计风速推算研究

利用岳阳气象站1952~2000年逐年10 min平均最大风速资料序列,对大风特征进行了分析.利用极值Ⅰ型分布曲线,推算出气象站处基本风速,结合沿江观测对比记录,通过比值法把基本风速推算到设计风速.结果表明:(1)桥位附近10 min平均年最大风速极值24.7 m/s,主导风为东北风;(2)不同重现期(100,50,30,10 a)10 m高处10 min平均年最大风速(基本风速)分别为27.6,25.4,23.8,20.2 m/s;(3)考虑5%的计算误差,桥位区不同重现期(100,50,30,10 a)10 m高处10 min平均年最大风速(设计风速)分别为29.0,26.7,25.0,21.2 m/s;(4)利用近地层风速随高度的指数变化推算出300 m以下每10 m高度层的设计基准风速.     

新龙门隧道射流通风与污染物浓度的数值分析

提出了按具有运动污染源一维非恒定流模型计算隧道内气流速度和污染物浓度分布的方法，并对新龙门隧道进行了数值模拟，阐述了内燃机产生的有害气体的危害及组成成分，并按牵引功率计算出排污量；分析了影响污染浓度分布的各种因素。对列车在洞内交会地点、射流风机工作时间及风机的送风方向进行了综合比较，提出了有利工况。     

高挡墙结构的爆破地震波响应仿真分析

通过大量爆破地震波数据的采集、研究和分析,建立了相应的爆破地震模拟震荡信号模型;作为模拟爆破地表震波的信号源,对于高挡墙建筑物进行了爆破地震波激励的动态响应的系列分析;采用瞬态动力学分析技术,确定了结构的动力响应,得到结构物的控制极限震速,为爆破地震监控提供了必要的技术参数。     

Optimization of pellet induration process parameters in rotary kiln using simulation results

In order to optimize the pellet induration process parameters of rotary kiln, a mathematical model is developed and solved by using commercial software FLUENT. Orthogonal regression has been performed by using the simulated results of the pellet induration in rotary kiln. The functional relations between pellet induration quality and process parameters are obtained. The main parameters that affect the thermal process are also discussed. The results show that the high temperature zone is longer and the local high-temperature could be avoided while the velocity difference decreases between the secondary air and the carbon particle, which leads to more uniformed temperature distribution and better induration quality. The work gives a right direction on improving of pellet production.     

Static Lateral Pressure of Heavy Haul Wagon's End Wall and Side Wall Based on Soil Mechanics

A new analytical method for calculating the static lateral pressure is introduced to obtain mor accurate calculation method of the lateral pressure of end wall and side wall of heavy haul wagon. In the method the influence of mechanical characteristics of the end wall is considered. For the end wall, a basis for lightweigh design can be provided by deducing an analytical method of the lateral pressure which is closer to the experimenta results. For the side wall, calculation formula which takes the cycloid curve as projection of rupture surface i derived. The correctness of the formula is verified by comparison with the experimental results. The nonlinea distribution of lateral pressure along the height direction is achieved. Further research is done on the relationship between key parameters of the heavy haul wagon design and the lateral pressure on the side wall.