Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

Early Estimation of Seismic Hazard for Strong Earthquakes in Taiwan

A shakemap system providing rapid estimates of strong ground shaking could be useful for emergency response providers in a damaging earthquake. A hybrid procedure, which combines site-dependent ground motion prediction models and the limited observations of the Real-Time Digital stream output system (RTD system operated by Central Weather Bureau, CWB), was set up to provide a high-resolution shakemap in a near-real-time manner after damaging earthquakes in Taiwan. One of the main factors that affect the result of ground motion prediction analysis is the existence of site effects. The purpose of this paper is to investigate the local site effects and their influence in the ground shaking and then establish an early estimation procedure of potential hazard for damaging earthquakes. Based on the attenuation law, the site effects of each TSMIP station are discussed in terms of a bias function that is site and intensity-level dependent function. The standard deviation of the site-dependent ground motion prediction model can be significantly reduced. The nonlinear behavior of ground soil is automatically taken into account in the intensity-level dependent bias function. Both the PGA and the spectral acceleration are studied in this study. Based on the RTD data, event correctors are calculated and applied to precisely estimate the shakemap of damaging earthquakes for emergency response.     

隧道锚杆锚固质量无损检测技术

采用声波检测仪,应用声频应力波在不同波阻抗面反射的能量和相位的变化原理,对隧道的锚杆锚固质量进行了无损检测试验研究,并对检测技术进行了探讨.结果表明,作为一种工程质量管理辅助手段,采用应力波对锚杆锚固质量进行无损检测,丰富了隧道围岩锚固质量检测方法,为隧道工程建设提供更好的质量保障.     

吉林台爆破料应力应变关系的三轴试验研究

通过室内大型三轴实验，研究了吉林台水库爆破料在不同级配下的应力应变关系，得到了在一定击实功下爆破料的最大干密度随细料含量变化的规律，分析了爆破料在不同级配和不同围压下应力与应变的变化规律、轴向应变与体积应变的关系及抗剪强度变化特性．从微观的角度说明了变化规律产生的原因，得出爆破料的抗剪强度随级配的变化而变化的规律．     

青藏高原中部的东西向扩张构造运动

系统分析了1933~2003年间青藏高原及其周缘发生的745个中、强地震的震源机制解,研究了高原地壳构造运动及其动力学特征。结果表明,大量正断层型地震集中发生在青藏高原中部海拔4000m以上的地区,其中许多地震是纯正断层型地震。震源机制结果显示,该区正断层型地震的断层走向多为南北方向,断层位错矢量的水平分量均位于近东西方向,这表明青藏高原高海拔地区存在着近东西方向的扩张构造运动。地震震源应力场的研究结果表明,在高原中部高海拔地区,E-W向或WNW-ESE向的水平扩张作用控制着该区的地壳应力场。青藏高原高海拔地区近东西方向的扩张构造运动是该区引张应力场的作用结果,其动力学原因可能与持续隆升的高原自重增大引起的重力崩塌及其周边区域构造应力状况有关。而青藏高原周缘地区,除了东部边缘外,南部的喜马拉雅山前沿以及青藏高原的北部、西部边缘所发生的绝大部分地震都是逆断层型或走滑逆断层型地震。在青藏高原周缘地区,北东或者北北东方向水平挤压的构造应力场为优势应力场。在中国西部的大范围内,主压应力P轴水平分量位于NE-SW方向,形成了一个广域的NE-SW方向的挤压应力场。青藏高原及其周缘应力场特征表明,印度板块的北上运动以及它与欧亚板块之间的碰撞所形成的挤压应力场是高原强烈隆起的直接原因。在青藏高原中南部形成了近东西向引张应力场为主的区域,并以东西向扩张构造运动部分释放其应力积累。研究高原高海拔地区的引张应力场和近东西向扩张构造运动的特征,对于认识青藏高原强烈隆起的地球动力学过程与机制,有着重要的理论意义。     

中国大陆科学钻探主孔声发射法现今地应力状态的确定

对中国大陆科学钻探主孔的岩心进行了声发射测量,确定了301～1531m深度的最大主应力。并与钻孔崩落法(深度1269～1655m)测量结果进行了对比,结果表明,声发射测量所得测值基本上落在钻孔崩落法测值的趋势线上,两种方法所得结果一致,说明测量结果可信。测量结果表明科学钻探主孔地应力大小随深度增加,在浅部301m最大主应力为13.4MPa,在深部1655m为55.2MPa。随深度的增加率为0.0279MPa/m。最大主应力方向为N54°±3.3°E,且方向不随深度变化。     

山东省临邑县临盘镇地裂冒喷水油灾害初步分析

1998年7月临邑县临盘镇西十二里村发生地裂冒喷水油灾害：对此次灾害研究认为，灾害发生前地应力已有明显改变，这种改变应是自然因素造成的，油田开发难以直接引发断层的活动。因此，地裂冒喷水、油灾害是因活动性断昙——临邑断层在局部薄弱部位活动而引起的，断层活动使地下流体顺断层涌出地表，是自然地质灾害。灾害的诱发因素是多方面的，其中地面不均衡沉降耦合使得地下应力布不均习，导致在断层的薄弱部位产生应力集中，并引发灾害发生。灾害发生已使得地下不均衡应力得以释放，近期内再次发生灾害的可能性不大。