COULOMB STRESS CHANGE ON ACTIVE FAULTS IN SICHUAN-YUNNAN REGION AND ITS IMPLICATIONS FOR SEISMIC HAZARD

Coulomb stress change on active faults is critical for seismic hazard analysis and has been widely used at home and abroad. The Sichuan-Yunnan region is one of the most tectonically and seismically active regions in Mainland China, considering some highly-populated cities and the historical earthquake records in this region, stress evolution and seismic hazard on these active faults capture much attention. From the physical principal, the occurrence of earthquakes will not only cause stress drop and strain energy release on the seismogenic faults, but also transfer stress to the surrounding faults, hence alter the shear and normal stress on the surrounding faults that may delay, hasten or even trigger subsequent earthquakes. Previously, most studies focus on the coseismic Coulomb stress change according to the elastic dislocation model. However, the gradually plentiful observation data attest to the importance of postseismic viscoelastic relaxation effect during the analysis of seismic interactions, stress evolution along faults and the cumulative effect on the longer time scale of the surrounding fault zone. In this paper, in order to assess the seismic hazard in Sichuan-Yunnan region, based on the elastic dislocation theory and the stratified viscoelastic model, we employ the PSGRN/PSCMP program to calculate the cumulative Coulomb stress change on the main boundary faults and in inner blocks in this region, by combining the influence of coseismic dislocations of the M≥7.0 historical strong earthquakes since the Yongsheng M7.8 earthquake in 1515 in Sichuan-Yunnan region and M≥8.0 events in the neighboring area, and the postseismic viscoelastic relaxation effect of the lower crust and upper mantle. The results show that the Coulomb stress change increases significantly in the south section of the Xianshuihe Fault, the Anninghe Fault, the northern section of the Xiaojiang Fault, the southern section of the Longmen Shan Fault, the intersection of the Chuxiong-Jianshui Fault and the Xiaojiang Fault, and the Shawan section of the Litang Fault, in which the cumulative Coulomb stress change exceeds 0.1MPa. The assuming different friction coefficient has little effect on the stress change, as for the strike-slip dominated faults, the shear stress change is much larger than the normal stress change, and the shear stress change is the main factor controlling the Coulomb stress change on the fault plane. Meanwhile, we compare the Coulomb stress change in the 10km and 15km depths, and find that for most faults, the results are slightly different. Additionally, based on the existing focal mechanism solutions, we add the focal mechanism solutions of the 5 675 small-medium earthquakes(2.5≤M≤4.9)in Sichuan-Yunnan region from January 2009 to July 2019, and invert the directions of the three principal stresses and the stress shape factor in 0.1°×0.1° grid points; by combining the grid search method, we compare the inverted stress tensors with that from the actual seismic data, and further obtain the optimal stress tensors. Then, we project the stress tensors on the two inverted nodal planes separately, and select the maximum Coulomb stress change to represent the stress change at the node. The results show that the cumulative Coulomb stress change increase in the triple-junction of Sichuan-Yunnan-Tibet region is also significant, and the stress change exceeds 0.1MPa. Comprehensive analysis of the Coulomb stress change, seismic gaps and seismicity parameters suggest that more attention should be paid to the Anninghe Fault, the northern section of the Xiaojiang Fault, the south section of the Xianshuihe Fault, the southern section of the Longmen Shan Fault and the triple-junction of the Sichuan-Yunnan-Tibet region. These results provide a basis for future seismic hazard analysis in the Sichuan-Yunnan region.     

Understanding of hydraulic properties from configurations of stochastically distributed fracture networks

A systematic investigation of the effect of configurations of stochastically distributed fracture networks on hydraulic behaviour for fractured rock masses could provide either quantitative or qualitative correlation between the structural configuration of the fracture network and its corresponding hydraulic behaviour, and enhance our understanding of appropriate application of groundwater flow and contaminant transport modelling in fractured rock masses. In this study, the effect of block sizes, intersection angles of fracture sets, standard deviations of fracture orientation, and fracture densities on directional block hydraulic conductivity and representative elementary volume is systematically investigated in two dimensions by implementing a numerical discrete fracture fluid flow model and incorporating stochastically distributed fracture configurations. It is shown from this investigation that the configuration of a stochastically distributed fracture network has a significant quantitative or qualitative effect on the hydraulic behaviour of fractured rock masses. Compared with the deterministic fracture configurations that have been extensively dealt with in a previous study, this investigation is expected to be more practical and adequate, since fracture geometry parameters are inherently stochastically distributed in the field. Moreover, the methodology and approach presented in this study may be generally applied to any fracture system in investigating the hydraulic behaviours from configurations of the fracture system while establishing a ‘bridge’ from the discrete fracture network flow modelling to equivalent continuum modelling in fractured rock masses. Copyright © 2007 John Wiley & Sons, Ltd.     

Analytical formulas for the geometric and inertia quantities of the largest removable blocks around tunnels

This paper presents algorithms for determining the vertices of the maximum removable block (MB) created by a joint pyramid (JP) around a tunnel when discontinuities are fully persistent. It is shown that an MB cannot be formed by more than 4 discontinuities and this drastically limits the proliferation of rock blocks that need to be analysed. The non‐convex block obtained after the MB is tunnelled through (real maximum block, RMB) is partitioned into a set of tetrahedra, and procedures are given for determining the vertices of these tetrahedra. Geometric and inertia quantities needed for stability analysis and support/reinforcement design are determined as functions of the calculated vertices' co‐ordinates. These quantities are: RMB's volume, face areas, perimeter of the excavated surface, centroid and inertia tensor. The algorithms for their calculation are at least two times faster than other algorithms previously proposed in other applications. It is shown that the formulations presented by Goodman and Shi for translational analysis and by Tonon for rotatability analysis can be used to analyse the RMBs using the geometric quantities presented. A numerical example is presented among those used to verify these analytical procedures and their implementation. Copyright © 2006 John Wiley & Sons, Ltd.     

不同洞口位置节能砌块隐形密框复合墙体恢复力模型研究

为研究洞口位置对节能砌块隐形密框复合墙体滞回特性的影响,对6个采用低周反复加载方式的1/2缩尺试件的试验现象及数据进行分析。基于试验结果和现有恢复力模型理论,对试验数据进行拟合分析,利用退化四折线模型和回归方法,建立节能砌块隐形密框复合墙体的四折线恢复力模型。试验结果表明:当洞口位于墙体中间时,试件滞回环愈加丰满,试件初始刚度有所增加,承载力衰减速率减小,卸载刚度退化速率变慢。通过对骨架曲线特征参数的计算,并和试验骨架曲线对比,两者吻合较好,表明该恢复力模型可以较好的反映出不同洞口位置对其滞回特性的影响,可为节能砌块隐形密框复合墙体的弹塑性分析及工程应用提供参考。     

青藏高原拉萨地块新生代超钾质岩与南北向地堑成因关系

青藏高原拉萨地块广泛分布有新生代超钾质岩,岩石地球化学和Sr-Nd-Pb同位素特征表明这些超钾质岩来源于与古俯冲环境有着密切联系的含金云母的富集地幔源区,它们主要喷发于25~10 Ma。同时在拉萨地块分布有多条南北向地堑(裂谷),且它们的切割深度可能到达下地壳的深部甚至岩石圈地幔,它们主要形成于23~8 Ma。拉萨地块大多数超钾质岩沿着新生代的南北向地堑(裂谷)分布,并且它们在形成时代和空间分布上存在着明显的耦合性,结合沿着印度-雅鲁藏布江缝合带分布的中新世埃达克质岩,笔者认为这些超钾质岩很可能与中新世早期北向俯冲的印度岩石圈沿着印度-雅鲁藏布江缝合带附近发生断离,以及由此而引起拉萨地块东西向伸展构造活动产生的南北向地堑(裂谷)系统有关。     

湘赣桂地区加里东期构造变形特征及成因分析

系统解剖了华南早古生代地层出露较好且加里东运动比较典型的地区,包括广西的元宝山、越城岭、大明山、大瑶山、云开大山地区以及湘赣边境等地区。通过对其褶皱、断裂形态的描述与分析发现大明山、大瑶山地区EW向的寒武系褶皱是云开地块在晚寒武世-早奥陶世由南向北推覆挤压的结果,而桂北元宝山、越城岭地区及湘赣边境地区NE-NNE向的早古生代地层的褶皱是由于华夏地块与扬子地块在晚奥陶世-早志留世沿郴州-临武断裂收缩挤压的结果,而且这一收缩挤压是属于陆内造山事件而不是前人所说的洋陆俯冲事件和陆-陆碰撞造山事件,且加里东运动先由南向北、后由东向西逐渐拓展,变形强度由强到弱。     

应用等性块段指数法划分新场井田的构造类型

介绍了等性块段指数法定量划分井田构造类型的基本原理。利用该法对筠连矿区新场井田的构造类型进行定量评价,得出该井田构造指数为1．98,属中等构造,证明原勘探所定构造类型基本正确。     

贺兰山高压泥质麻粒岩——华北克拉通西部陆块拼合的岩石学证据

贺兰山孔兹岩系作为华北克拉通西部孔兹岩带的重要组成部分,其成因问题一直存在争议。近来,在对贺兰山富铝岩系研究过程中,发现了具有"蓝晶石+条纹长石+石榴石"特征组合的高压泥质麻粒岩。通过岩相观测与相平衡定量模拟,确定其曾经历四个阶段的变质演化过程,即早期进变质阶段(M1)、变质峰期阶段(M2)、峰期后降压阶段(M3)和晚期退变阶段(M4)。早期进变质阶段,以石榴石中包裹黑云母和白云母为特征。变质峰期阶段,以出现"蓝晶石+条纹长石+石榴石"组合为特征,变质温、压条件为850～870℃、1.4～1.5GPa。峰期后降压阶段以峰期石榴石和蓝晶石转变为堇青石和夕线石为特征,温压条件大致为840～860℃、0.6～0.8GPa。晚期退变阶段以出现新生黑云母为标志,温、压条件为780～810℃、0.45～0.55GPa。PT轨迹具有顺时针演化样式,峰期后表现为先发生近等温降压(ITD)、后发生近等压冷却(IBC)的退变过程,反映陆-陆碰撞发生高压麻粒岩相变质作用,后又快速折返到正常地壳深度的变质动力学过程,进而从岩石学角度证明华北克拉通西部孔兹岩带是由北部的阴山陆块与南部的鄂尔多斯陆块发生碰撞所致。     

华北陆块南缘中生代合峪花岗岩的地球化学特征与成因

合峪花岗岩基位于华北陆块南缘外方山地区,为豫西地区燕山期最大的岩基,出露总面积达784km2,有多期侵入的特点,可划分为6个单元。主要岩石类型为黑云母二长花岗岩。LA-ICP-MS锆石U-Pb年龄资料表明,合峪花岗岩基侵位时间为148.2~135.3Ma。合峪花岗岩基的SiO2=67.16%~75.43%,Al2O3=13.29%~15.92%,MgO=0.26%~1.12%,K2O+Na2O8%,Na2O/K2O=0.88~1.43,属于高钾钙碱性系列,ACNK=0.94~1.09,为准铝质-弱过铝质花岗岩。岩体轻稀土富集、重稀土亏损((La/Yb)N=14.5×10-6~49.9×10-6,平均值27.2×10-6),Sr含量变化较大(102×10-6~848×10-6,平均290×10-6),Y、Yb含量低(Y=3.21×10-6~17.3×10-6;Yb=0.43×10-6~2.16×10-6),Eu弱亏损(δEu=0.57~0.89),反映熔体发生过长石分离结晶作用,源岩部分熔融生成熔体时残留相组合中没有或很少有长石的存在。合峪花岗岩基ISr=0.7071~0.7090,εSr(t)=40.8~65.9,εNd(t)=-16.4~-11.2,其tDM2为1.85~2.27Ga;锆石的εHf(t)主要集中于-25.39~-5.25之间,tDM2年龄介于1301~2846Ma。同位素数据显示合峪花岗岩基的源岩主要为南秦岭及扬子地块结晶基底,还混有少量熊耳群及太华群的物质。综合区域地质演化,认为合峪花岗岩基形成于秦岭造山带中生代陆内俯冲,岩石圈增厚,幔源岩浆上涌底侵导致下地壳部分熔融形成。     

中国及邻区重力特征与块体构造单元初划

以重力资料为基础,结合其分析结果和区域地质构造现象,以及前人研究成果,在总结中国海陆地质地球物理系列图编制项目组和前人对中国岩石圈三维结构宏观特征研究成果基础上,将块体构造理论应用于解释中国海陆宏观大地构造格架,并将其划分为3个构造域和20多个块体。对3个构造域及块体构造单元特征进行了分析,指出印支运动使中国岩石圈应力场发生了根本性变化,并为现今宏观大地构造格局的形成创造了条件。