新疆精河M_W6.3地震产生的静态应力变化研究

利用精河M_W6.3地震有限断层破裂模型,计算了精河地震产生的位移场、应力场、周围主要断层上的静态库仑应力变化以及主震对余震的触发作用。结果表明:(1)精河地震产生的地表隆升最大值约为6.6cm,沉降最大值约为1.8cm;水平位移方向呈现震中南北侧向震中汇聚、震中东西侧向外"流出"的特点。(2)精河地震产生的水平面应力场展布南北侧物质主要受到指向震中的拉张力作用,东西两侧物质主要受到因震中过剩物质东西向排出而导致的东西向挤压力作用。(3)震中西侧距震中约20km的库松木契克山前断裂中段和震中东北部距震中约50km的四棵树-古尔图南断裂西段的库仑应力加载均大于0.01MPa,即2处为地震危险区。(4)在震源深度为8~12km的余震事件中,约有85.5%处于库仑应力加载区,即受到主震的的触发作用;在深度为4~8km的余震事件中,约有87%受到主震的应力触发作用。     

REPRESENTATIONS OF COULOMB FRICTION FOR DYNAMIC ANALYSIS

In many engineering problems involving friction, the friction is treated as Coulomb friction, where the magnitude of the friction force is constant but its direction is always opposite to that of the sliding velocity. In dynamic problems, the direction of the sliding velocity can change quite often. The many changes of the velocity direction cause many discontinuities in the friction force, complicating the process of evaluating the response of systems involving friction. In this paper, it is shown that the discontinuous Coulomb friction force can be represented by at least four different continuous functions. Each of these functions involves one constant that controls the level of accuracy of that function's representation of the friction force. The accuracy of the various representations is verified by comparing the response of a single degree freedom system, obtained through numerical solutions utilizing these representations, with an exact analytical solution. © 1997 by John Wiley & Sons, Ltd.     

动态应力触发研究进展

本文从库仑破裂应力与动态应力触发原理、观测事实与研究现状以及存在的问题与展望等3个方面介绍了动态应力触发的研究进展。其中,在主震对余震活动、强震对后续远场地震活动、不同破裂类型的地震之间及强震对火山(泥火山)的动态应力触发作用、一次大地震中子事件之间的动态应力触发关系以及动态应力触发的时间延迟等方面进行了详细讨论。分析认为,动态应力触发理论在解释震后余震分布、远程触发以及对火山触发性喷发等方面得到了较好的结果;动态应力触发可能会受断层破裂方式和方向的影响;大地震中子事件之间的动态应力触发关系对震源研究有重要意义;动态应力触发现象普遍存在时间延迟。     

2010智利莫尔M_W8.8地震对其破裂区及附近地区地震活动性的影响

强震的发生与地壳中应力增加有关。库仑破裂应力被用于研究由于某次强震发生所产生的应力变化。2010年智利莫尔M_W8.8地震发生之后,其破裂区以南地区没有发生M≥6.0地震,而在其以北地区于2015年在智利发生了伊亚佩尔M_W8.3级地震。在这两个地区地震活动性悬殊。本文利用库仑破裂应力研究了2010年智利莫尔M_W8.8地震在这两个地区产生的应力变化。结果表明,2010年智利莫尔M_W8.8地震在其破裂区以北地区产生的应力变化为正值,有利于断层错动发震,而在其破裂区以南地区产生的应力变化为负值,具有抑制断层错动发震的作用。     

2015年四川乐山金口河M5.0地震孕震构造条件分析

2015年1月14日乐山金口河M5.0地震发生在历史地震强度较低的川南山区与四川盆地交界一带。基于四川区域地震台网的震相报告与波形资料,采用双差定位法对地震序列进行重新定位,同时,采用CAP波形反演方法及HASH方法反演了主震及序列中8次M_L≥2.0地震的震源机制解。另外,利用Coulomb3计算了主震发生后库仑应力改变量,得到的结果如下：①重新定位结果显示,金口河M5.0地震位于（103.18°E,29.32°N）,震源深度16.6km,略深于波形反演结果（12km）。序列分布在NNW向天全-荥经断裂和NE向西河-美姑断裂的交汇部位,余震序列在空间上呈NE向展布。②M5.0主震的机制解为节面Ⅰ：走向350°/倾角46°/滑动角107°,节面Ⅱ：走向146°/倾角47°/滑动角73°,表现为走向NW（NNW）、中等倾角的逆冲型运动方式。序列中其余8次M_L≥2.0余震大多以走向NE的逆冲型地震为主,个别为走滑或正断层类型。主震和大部分余震的节面方向不一致,主震节面方向与余震长轴方向也不一致。③主震后库仑应力改变量显示,余震主要发生在主震引起的库仑破裂应力增加的区域。综合分析推测,NNW向天全-荥经断裂为本次地震主震的发震构造,倾向NE的机制解节面Ⅰ指出了该断裂的几何产状;M5.0主震发生后,立即触发了其旁侧的NE向西河-美姑断裂,并激发了多次余震。     

主动与被动状态下墙体侧向位移近似计算

极限状态下墙体侧向位移对土压力计算和支挡结构设计影响显著。根据Rankine变形体和Coulomb刚塑体模型,将墙后土体变形分别当作单剪和直剪试验中试样的剪切过程,以达到极限剪切变形（剪应变或单位长度剪切位移）作为进入主被动状态标准,构建了土体变形与墙体位移的几何关系,提出了反映土体变形与强度特性,同时考虑静止时初始应力状态影响的墙体极限侧向位移近似计算模型。分析表明：土体极限剪切变形、滑移区范围、初始应力状态是影响墙体极限位移的核心要素,其中极限剪切变形占据主导作用,是导致不同颗粒组成及密实程度土体进入极限状态所需墙体位移差异显著的主要原因,而主被动区范围不同和因静止土压力系数 1引起的初始剪切变形,则是被动状态墙体位移远大于主动的关键因素;算例中主动与被动状态下墙体位移与墙高之比分别介于0.5‰～13.2‰和?0.4%～?5.2%,且主动状态下细粒土墙体位移大于粗粒土,计算结果与工程经验及相关文献模型试验基本一致。     

2014年鲁甸6.5级地震成因、破裂特征和余震分布特征的库仑应力作用

根据2014年鲁甸M_S6.5地震的区域构造和余震共轭分布特征,本文首先计算了1733年小江断裂带北段M7.75地震,1850年则木河断裂带M7.5地震和1974年马边M_S7.1地震对鲁甸M_S6.5地震震源机制解两个节面的黏弹性库仑应力作用,结果显示NNW向发生主破裂的包谷垴-小河断裂受到这3次地震,尤其是1850年M7.5地震明显的库仑应力作用,我们认为则木河断裂的高速左旋走滑运动以及7级以上强震的重复发生对于包谷垴-小河断裂的强震孕育和断裂演化方面具有促进作用;然后分析了鲁甸6.5级地震的共轭破裂与余震分布特征,并分别计算了两个共轭破裂面单独破裂对另一破裂面的库仑应力作用,结果显示NEE向破裂促进NNW向破裂的发生,而NNW向破裂后则阻碍了NEE向破裂的进一步发展,最终发展成以NNW向破裂为主的共轭破裂事件;最后计算了共轭破裂产生的库仑应力变化对余震的影响,认为位于NEE向破裂面西侧的余震集中分布主要是由于应力触发作用而形成.     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

断层地震孕育和发生的不稳定性模型

在以前的工作中,考虑直立走滑型断层地震,假设断层面微元破裂强度遵循Weibull概率分布,由细观力学方法推导出断层面的宏观本构关系是一个非线性函数,表现为弹性-软化塑性特征,在此基础上用稳定性理论研究了地震稳定性问题.而实际断层大多是倾斜的,为此,本文首先建立了由围岩和倾斜断层构成的平面地震力学模型,采用宏观的断层载荷-变形的全过程曲线,详细讨论了倾斜断层地震的不稳定性问题.结果表明,远场一旦施加位移,断层也同时错动,这可能与实际情况不符合.为了更好的模拟断层的初始能量累计过程,进一步对断层本构模型进行改进.考虑断层面破裂强度,采用Coulomb破裂准则,则断层表现为刚塑性本构关系,只有当断层面剪应力达到一个临界值时,断层才开始错动.研究表明,对于倾斜断层地震,与直立走滑型断层地震一样,系统刚度比β（围岩切线刚度与断层刚度最大值之比）是决定地震失稳的重要参数,只有当β<1时才会出现地震失稳,且伴随应力突跳和围岩应变能释放.当β≥1时,仅仅是断层无震滑动,不会发生地震.在远场应施以位移形式边界条件,以致地震失稳发生在平衡路径的位移转向点并伴有应力突跳.     

断层横向构造在逆冲型地震破裂中的作用——以汶川地震小鱼洞断层为例

大地震破裂大多由横向构造(如阶区、弯曲和分叉)所分割的多个段落组成.2008年5·12汶川地震破裂沿北东走向上穿过了多个横向构造部位,特别在震中北东45 km的位置,小鱼洞断层、北川断层和彭灌断层三者之间呈现复杂的断裂切割相交关系.复杂断层几何结构对破裂的扩展是有抑制还是促进的作用？在相交的断裂段之间是否存在最优的破裂顺序？本文以库仑应力分析为手段,探讨在汶川同震破裂初始30 s内,破裂在多分支断裂中选择扩展路径时的可能应力相互作用.库仑应力分析显示:如果北川断层先发生破裂,其滑动对小鱼洞断层和彭灌断层均产生强烈负应力的抑制作用,而彭灌断层的滑动却反而对小鱼洞断层和北川断层浅部有强烈正应力的促进作用.因此,从准静态应力分析角度,彭灌断层先于北川断层发生破裂的可能性较大,这一破裂顺序与小鱼洞断层参与同震破裂过程的事实相符.此外,小鱼洞断层在链接北川和彭灌断层的同震位移中可能起到桥梁作用,但非静态应力的影响.横向构造在逆冲型地震破裂扩展过程中起到的牵引作用使得逆冲型地震破裂能够比走滑型地震跨越更宽的阶区.横向构造是逆冲断裂带内广泛发育的构成单元,因此在地震危险性分析的最大潜在震级测算中应该考虑其作用.