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.     

The First Appearance of Cattle in Denmark Occurred 6000 Years Ago: An Effect of Cultural or Climate and Environmental Changes

Stable carbon and nitrogen isotope ratios from bones of contemporaneous Late Atlantic aurochs and early cattle in eastern Denmark are significantly different and provide information on the origin and feeding strategies of the earliest domestic cattle. The data show that the early cattle were feeding on grass right from the beginning 4000 cal. yr BC. In contrast, the youngest aurochs population primarily browsed and grazed from the dense forest floor resulting in rather negative δ¹³C values measured on bone collagen. The oldest aurochs have similar isotope values to the earlier cattle, whereas the youngest aurochs have similar values to Late Atlantic red deer from the same locality. As eastern Denmark was largely covered by forest, speculations on the origin of the grazing areas are many. The grass may have grown in openings in the forest, at the forest fringe, or more likely on the newly reclaimed coastal land areas exposed by the decreasing rate of eustatic sea‐level rise contemporaneously with isostatic uplift, during the Littorina transgressions. The stable isotope values do not indicate that leaf foddering of the early cattle was of importance.     

On global forces of nature driving the Earth’s climate. Are humans involved?

The authors identify and describe the following global forces of nature driving the Earth’s climate: (1) solar radiation as a dominant external energy supplier to the Earth, (2) outgassing as a major supplier of gases to the World Ocean and the atmosphere, and, possibly, (3) microbial activities generating and consuming atmospheric gases at the interface of lithosphere and atmosphere. The writers provide quantitative estimates of the scope and extent of their corresponding effects on the Earth’s climate. Quantitative comparison of the scope and extent of the forces of nature and anthropogenic influences on the Earth’s climate is especially important at the time of broad-scale public debates on current global warming. The writers show that the human-induced climatic changes are negligible.     

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

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

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

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

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

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

滇黔桂盆地及其邻区石炭纪至二叠纪层序地层格架及三级海平面变化的全球对比

属于被动大陆边缘裂谷盆地的滇黔桂盆地，自加里东运动之后，泥盆纪开始拉开，整个晚古生代延续发展。受同生断裂的控制，在深水盆地中发育大小不等的孤立碳酸盐台地，因而在滇黔桂盆地及其邻区形成特殊的“台-盆-丘-槽”的古地理格局。对不同古地理背景下的典型剖面进行三级层序划分并对其进行空间追索和对比，在滇黔桂盆地及其邻区的石炭纪和二叠纪地层中可以识别出12个三级层序，相当于晚古生代25个三级沉积层序(2-5Ma；SQ1至SQ25)中的SQ24至SQ25；以地层记录中的的两种相变面和两种穿时性为基本要素，可以建立研究区石炭系和二叠系的层序地层格架；层序地层格架反映了三级层序的基本特征：空间上相序的有序性和时间上环境变化的同步性。研究区的石炭系和二叠系，主要为一套碳酸盐岩地层，其中在连陆台地上发育3套煤系地层，它们分别组成石炭纪三级层序SQ15跨系的三级层序SQ19的HST以及二叠纪吴家坪期的三级层序SQ24；而且在连陆台地边缘和孤立台地上，于阳新世的茅口亚世以及乐平世的长兴期发育海绵生物礁，海绵生物礁分别构成了以下二叠系三级层序即阳新世的SQ22和SQ23以及长兴期的SQ25。3套煤系地层和两套海绵生物礁的发育，使研究区的石炭系和二叠系形成了一个与欧美地区完全不同的层序地层序列。在研究区石炭系和二叠系所识别出的12个三级层序中，石炭纪一二叠纪船山世地层可以划分出6个三级层序。这种划分与Busch等对北美相同层位的划分相似，也就是说该6个三级层序的形成时限大于10Ma；较长的形成时限可能反映了形成在联合古陆汇聚时期即全球构造相对稳定时期三级海平面变化的一个特性。因此，研究区的石炭系和二叠系所识别和划分出的12个三级层序，其数量比Ross和Ross(1985)对欧美地区石炭系和二叠系所划分出的50多个三级层序要少得多。     

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

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

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

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