Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

平衡剖面的制作流程及其地质意义

平衡剖面技术是地质思维和计算机技术的结晶,使对断层构造的研究提高到定量阶段,其依据是在垂直构造走向的剖面上,地层长度和面积(2D)或体积(3D)是均衡的。在此原理基础上利用数学手段对盆地的构造发育史进行正演和反演模拟,直观地再现地下构造的原始几何形态,迅速提供地震剖面的构造解释方案,并对解释结果进行检验(不平衡的剖面其解释一般有问题),为深刻认识构造发育史、分析油气运移及聚集规律提供依据,提高了工作效率。其结果也为盆地模拟、油藏模拟、定量计算构造伸缩量等地质研究打下了坚实的基础[1]。     

用双三次样条函数和GPS资料反演现今中国大陆构造形变场

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.     

Bacterial productivity in the Prydz Bay and its adjacent waters,Antarctic

1　IntroductionBacteriaandtheiractivitiesplayanimportantroleintheelementalbiogeochemicalcyclesandenergytransformingintheocean (Zhenetal.1 997) .DortchandPackard(1 989) proposedthatfoodwebsintheeutrophicwatersaredominatedbythebiomassofprimaryproducerswhilefoodwebsintheoligotrophicwatersaredominatedbythebiomassofmicrobes.Heterotrophicbacteriahadbeenshowntoplayanimportantroleinthedecompositionoflarge ,rapidlysinkingorganicparticleswithinandbelowtheeuphot iczone ,andfurthertoaffecttheelementaldyn…     

莱州湾南岸潍河下游地区咸水入侵灾害成因及特征

通过对莱州湾南岸咸水入侵较严重地潍河下游地区晚更新世以来沉积特征及现代自然环境条件变化的分析，探讨了沉积相对咸水入侵产生及空间范围特征的环境机理。晚更新世以来的三次海平面升降变化造成了潍河下游地区海陆沉积环境交替，形成了巨厚的海陆交互相沉积层。海进时期，大面积的滨海平原被淹没，在近海平原洼地滞留的海水经过蒸发、浓缩变为卤水，成为咸水入侵的物源；海退后陆源碎屑在滨海地区沉积形成了巨厚的古河道砂层。20世纪70年代末期以来，随着对地下淡水的过度开采，淡咸水水头压力差减小．卤水通过古河道砂层快速南侵。通过对潍河下游地区100余个地质钻孔水化学连续监测资料分析，阐明了咸水入侵的特征，有针对性地提出了咸水入侵的防治措施。     

西南矿区山体崩塌成因机制分析及防治对策

西南地区的矿区具有相似的地理、地质环境,通过四川芙蓉煤矿白岩崩塌的典型实例,对西南矿区山体崩塌的成灾环境进行调查研究,揭示了西南矿区山体崩塌的成因机制,并简要提出了对该煤矿开发与综合治理建议,指出了加强防治工作的紧迫性,强调开展系统工程防灾的重要性,并为该煤矿解决正常采掘接替及工作面的布置问题作出了科学的指导.     

黄土高原古植被与再造山川秀美

要想恢复秀美的山川,主要在于恢复被人类严重破坏的原始植被.黄土高原的原始植被到底是什么样子,学者们说法不一.孢粉分析是最重要的方法.目前,研究结果已揭示出黄土高原古植被在第四纪以来显示着时间上、空间上的变化.古土壤代表间冰期夏季风盛行的湿热气候,植被比较繁盛;黄土代表冰期冬季风盛行的干冷气候,植被凋凌.全新世中期5000～6000年前,现代间冰期气候最佳时期,半坡人生存之时,森林茂密并有亚热带植物入侵,但兰州一带1万年来一直保持着荒漠草原的面貌.因此,黄土高原植被的恢复应参照各地的原始植被并考虑现代气候条件,因地制宜地实施.在东南部可以恢复森林,中西部可以恢复森林草原,西北部只能种草恢复草原.     

多条太阳光球FeI线在磁场中形成及形成深度的综合分析(Ⅰ)

针对多通道滤光器太阳磁场望远镜的磁场观测定标及掌握(所采用的)谱线特征之需要,取VAL宁静太阳大气模型计算了7条Fel光球线的Stokes轮廓、形成深度、贡献函数分布,从而较为系统地对多条反常及正常Zeeman线的特征及性质作出了分析与总结,解释了在磁光效应影响下Stokes Q,U参量的形成深度曲线在近线心区域处出现陡峭峰值的原因。     

长江中下游地区地下水卤素元素的形成及其分布规律

本文以丰富的实际资料，论证了地下水的卤素元素（F、Cl、Br、I）的形成、含量及其分布规律与含水介质成分、上覆岩土性质、地下水退流条件、氧化还原环境、地下水矿化度之间的关系。根据江汉平原东部区和鄱阳湖区地下水中Br、I元素的调查研究结果及它们形成的控制因素与分布规律，结合长江三角洲南部区水文地球化学环境条件分析对比，指出该区是一个形成Br、I矿泉水的有利地区。