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.     

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

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

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

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

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

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

邯峰矿区古构造应力场与构造演化

采用共轭剪节理应力反演方法,恢复了邯郸-峰峰矿区晚古生代以来的3期古构造应力场,进而探讨了煤田构造的演化历史,将其分为4大阶段：①中生代早期近NS向挤压,煤系后期改造初动期;②中生代晚期SE-NW向挤压,奠定煤田构造格架的基础;③中生代末至古近纪NW-SE向拉张,煤田构造格架定型;④新近纪以来近东西向拉张,煤田构造的现代活动.     

Laboratory Studies Of Wind Stress Over Surface Waves

Simultaneous laboratory observations of wind speed, wind stress, and surfacewind-wave spectra are made under a variety of wind forcing patterns using cleanwater as well as water containing an artificial surfactant. Under typical experimentalconditions, more than half of the total stress is supported by the wave-induced stressrather than by the surface viscous stress. When the surfactant reduces the shortwind-wave spectra, the wind stress also decreases by as much as 20–30% at agiven wind speed. When the wind forcing is modulated in time, the wind stresstends to be higher under decreasing wind than under increasing wind at a givenwind speed, mainly because the response of short wind-wave spectra to varyingwind forcing is delayed in time. These examples clearly demonstrate that therelationship between the wind speed and the wind stress can be significantlymodified if the surface wave field is not in equilibrium with the wind forcing.Next, we examine whether the wind stress is estimated accurately if the wave-inducedstress by all surface wave components is explicitly evaluated by linear superpositionand is added to the surface viscous stress. It is assumed that the surface viscous stressis uniquely related to the wind speed, and that the wind input rate is determined by thelocal, reduced turbulent stress rather than the total stress. Our wind stress estimatesincluding the wave contributions agree well with observed wind stress values, evenif the surface wave field is away from its equilibrium with the wind in the presenceof surface films and/or under time-transient wind forcing. These observations stronglysuggest that the wind stress is accurately evaluated as a sum of the wave-induced stressand the surface viscous stress. At very high winds, our stress estimates tend to be lowerthan the observations. We suspect that this is because of the enhancement of wind stressover very steep (or breaking) short wind-waves.     

阿尔金地区构造应力场及其对金属矿产分布的控制作用

构造应力场研究可以为区域矿产预测提供依据。主要根据节理、擦痕的测试分析，结合区域构造解析，确定阿尔金山东西向拉配泉－－红柳沟构造带经历了三期不同方向的应力作用：印支期以前的南北向挤压作用，印支－燕山期的北西－南东向的挤压和中生代晚期－新生代北东东－南南西向的挤压；并利用有限元数值模拟前两期构造应力场的演变过程，进而讨论了区域构造应力场对内生金属矿产分布规律的控制作用指出了区域寻找大型内生金属矿产的有利地段。     

A method to predict the group fissuring and faulting caused by regional groundwater decline

Ground fissuring is a recurrent problem in many countries where water extraction surpasses the natural recharge of aquifers. Due to differential settlement, the soil layer undergoes deformation and cracks with serious consequences for civil infrastructure. Here, we propose an approximate analysis of the fissuring process that can be used to predict the location of cracks, which increasingly affect some middle- and large-sized cities in the world. For that purpose, the ground loss theory is applied to sediments overlying a sinusoidal-shaped graben. This analysis shows the existence of a tensile zone at the border of the graben with maximal values on its shoulder where tension cracks are more likely to appear. It also shows that soil deformation under differential settlement may evolve into ground faulting if water withdrawal continues. Finally, when a crack has completely developed, the tensile zone shifts towards the center of the graben, creating a new area for potential cracking and faulting.     

Crack damage stress as a composite function of porosity and elastic matrix stiffness in dolomites and limestones

Twenty-five uniaxial compression tests were performed to determine stress at onset of dilation, referred to herein as “the crack damage stress,” in heterogeneous dolomites and limestones. A simplified model for crack damage stress (σ_cd) is developed here using porosity, elastic modulus, Poisson's ratio and three empirical coefficients. The model shows that when porosity decreases and elastic modulus increases, σ_cd rapidly increases and approaches its maximum value. On the other hand, when porosity increases and elastic modulus decreases, σ_cd rapidly decreases and approaches its minimum value. The proposed model is validated for six heterogeneous limestone and dolomite formations which are widely distributed in Israel.