Contained density currents with high volume of release

Contained density currents with high volume of release reflect against the boundaries of the reception environment commonly leading to oscillatory flow. These flows exist in sediment clarifiers, compromising their operations, and deposited signatures of contained turbidity currents are found as part of the infill of sedimentary basins; for operation of the former and interpretation of the latter it is essential to understand the dynamic processes of these flows. Six lock‐exchange experiments with different initial densities were made in a horizontal flume, where the volume of the saline mixture in the lock was equivalent to the volume of the ambient fluid. A further two tests, with a repeated initial density, were made: one with high volume of release and very long duration; and another with low volume of release. Firstly, the movement of the current is discussed, including the oscillations within the experimental tank involving the density current and an upper layer counter‐current. It is shown that the cyclic behaviour is self‐similar with the reduced gravity of the initial density in the lock. Secondly, entrainment and water mixing processes are characterized. The time evolution of mixing is characterized qualitatively by analysing the background potential energy of the density distributions to show that mixing occurs even in the earlier stages of the current, and mainly within the first cycle of the oscillation. Quantified analysis reveals that, in currents with high volume of release, entrainment discharge is one order of magnitude higher, mainly due to the larger interface between the ambient fluid and the current. A model for the evolution of the mixing process is proposed for density currents with high volume of release. Finally, the dynamics of the head of the current is analysed. The entrainment in the head, when compared to the entrainment in the remainder body of the flows, is less important for the configuration with a larger lock.     

Bounding surface model for soil resistance to cyclic lateral pile displacements with arbitrary direction

The development of a two-surface elastic–plastic bounding surface P–Y model for cyclic lateral pile motions is described. The kinematic-hardening model is applicable to the analysis of pile foundations subjected to loading with arbitrary azimuths relative to the pile axis. The model realistically captures the hysteretic energy damping associated with dynamic loading of subsea foundations through physically correct plastic mechanisms and provides results consistent with those observed in physical tests including cyclic loading. Its performance is demonstrated in element states of stress and in pile foundation analyses. The development based on the incremental theory of plasticity results in more robust solutions than may be obtained using alternative elastic, variable moduli and deformation plasticity formulations.     

体素模型的战场环境建模研究

针对基于面模型的战场环境建模难以展现战场环境要素的内部属性和规律等的不足,该文提出了利用体素模型实现战场环境建模的基本方法。该方法能够记录战场环境各要素的不同属性,反映其内部属性和规律,准确地模拟战场地理环境的变换过程,可实现战场环境各要素的统一建模。文章利用体素模型实现了对战场环境的认知与抽象,对战场环境体素数据模型进行了分类。最后以战场地形环境要素为例,利用战场环境基本体元实现了对战场地形环境的体素建模。     

Coupled Model of Two-phase Debris Flow, Sediment Transport and Morphological Evolution

The volume fraction of the solid and liquid phase of debris flows,which evolves simultaneously across terrains,largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper,we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth,the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress,and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition(Iverson and Ouyang,2015) is presented. In this formula,the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe's Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.     

Time-domain simulation of berthing problem between FPSO and shuttle tanker in waves

A time-domain simulation method based on potential flow model has been developed to investigate the berthing problem between two floating bodies in wave. The boundary value problem is formulated with respect to an earth-fixed coordinate system because the relative positions of the two vessels continuously change during the berthing operation. The classical finite element method is used to solve the Laplace equation in the fluid domain with moving boundary. The linearized free-surface boundary conditions are integrated in time by applying 4th-order Adams–Bashforth–Moulton method. A simple re-mesh algorithm with local and global mesh systems is introduced to update mesh by considering large horizontal movement of the berthing vessel. The developed numerical method is used to investigate the berthing problem between a FPSO and shuttle tanker in waves. The focus is on the wave-induced motion response during the berthing process. The characteristics of the motion responses in berthing operation are examined with various wave frequencies, berthing speeds and wave headings.     

Incremental elastoplastic FEM for simulating the deformation process of suction caissons subjected to cyclic loads in soft clays

This paper presents an incremental elastoplastic finite element method (FEM) to simulate the undrained deformation process of suction caisson foundations subjected to cyclic loads in soft clays. The method is developed by encoding the total-stress-based bounding surface model proposed by the authors in the ABAQUS software package. According to the model characteristics, elastoplastic stress states associated with the incremental strains of each iteration are determined using the sub-incremental explicit Euler algorithm, and the state parameters describing the cyclic accumulative rates of strains are updated by setting state variables during the calculations. The radial fallback method is also proposed to modify the stress states outside the bounding surface to the surface during determination of the elastoplastic stress states. The stress reversals of soil elements are judged by the angle between the incremental deviatoric stress and the exterior normal vector at the image stress point on the bounding surface to update the mapping centre and state variables during cyclic loading. To assess the general validity of the method, the reduced scale model tests and centrifuge tests of suction caissons subjected to cyclic loads are simulated using the method. Predictions are in relative good agreement with test results. Compared with the limit equilibrium and quasi-static methods, the method can not only determine the cyclic bearing capacity, but can also analyse the deformation process and the failure mechanisms of suction caisson under cyclic loads in soft clays.     

近60年来长江河口河势变化及其对水动力和盐水入侵的影响I.河势变化

河势是影响河口水动力和盐水入侵基本因子。本文利用20世纪50和70年代长江河口海图,数值化岸线和水深,结合2012年长江河口实测水深资料,分析长江河口自50年代以来的河势变化。长江河口为分汊河口,50年代仅为二级分汊,至70年代才形成三级分汊,四口入海的河势格局。70年代相比于50年代,北支淤浅严重,其上、中、下段容积变化分别为-64.13×106、-306.60×106和-639.27×106 m3,对应的变化率分别为-16.30%、-22.74%和-25.69%,均显著减小;南支的上、中、下段容积变化分别为-28.61×106、-35.69×106和126.43×106 m3,相应的变化率分别为-1.30%、-2.12%和4.36%;北港由于崇明浅滩和横沙浅滩的淤浅,下段容积明显减小,其上段和下段容积变化分别为109.21×106和-797.14×106 m3,对应的变化率分别为5.01%和-15.25%;南港上段由于河道淤浅容积减小,下段北由于铜沙浅滩被冲开形成北槽,导致水深变深、容积增加,其上段、下段北和下段南容积变化分别为-238.95×106、203.58×106和153.34×106 m3,对应的变化率分别为-8.96%、6.85%和3.26%。2012年相比于70年代,北支由于大量淤浅和围垦容积大幅减小,其上、中、下段容积变化分别为-199.06×106、-504.61×106和-654.12×106 m3,对应的变化率分别为-60.45%、-48.44%和-35.38%;南支的上、中、下段容积变化分别为92.34×106、193.01×106和-163.62×106 m3,相应的变化率分别为4.24%、11.73%和-5.40%;北港上段青草沙水库的围垦和下段横沙东滩的围垦造成面积和容积减小,其上段和下段容积变化分别为-154.64×106和-511.79×106 m3,对应的变化率分别为-6.75%和-11.55%;南港由于上段河道刷深而下段九段沙以及南汇边滩淤浅、围垦,导致其容积上段增加,下段减小,上段、下段北和下段南容积变化分别为136.39×106、-658.28×106和-1266.11×106 m3,对应的变化率分别为5.62%、-20.73%和-26.06%。     

塔里木溢流玄武岩火山通道的三维结构及其热成因气体释放

塔里木上奥陶统-志留系沉积地层中广泛发育早二叠世溢流玄武岩的火山通道相岩床-岩墙网络。三维地震数据解释结果显示,这些火山通道以平行围岩地层的岩床和斜切围岩地层的"碟状岩床"为主要特征。在玄武岩喷发过程中,火山通道岩浆的热量可以导致沉积围岩发生热接触变质并将沉积围岩中的有机质转化为"热成因气体"。在塔北英买2井区火山通道烘烤沉积围岩模型基础上,利用有限元热模拟方法确定了该区早二叠世玄武岩喷发时火山通道热烘烤影响范围随时间的变化。基于沉积围岩有机质丰度估算,该区热烘烤成因甲烷释放量可达11.3Gt(即113亿吨)。如果整个塔里木溢流玄武岩省具有与英买2地区相同的释放强度,则塔里木溢流玄武岩省活动期间释放的甲烷总量可达7062.5Gt,必然导致非常显著的环境效应。同时,玄武岩火山通道岩浆引起的热接触变质作用对已存在的油藏具有明显的破坏作用,塔里木盆地古生界总量约8~10Gt的油藏破坏和大量沥青的形成可能与此有关。     

准噶尔盆地吉木萨尔凹陷梧桐沟组层序地层划分

准噶尔盆地吉木萨尔凹陷二叠系梧桐沟组作为一套重要的储集岩,具有良好的油气成藏前景,但对其层序划分及地层沉积样式的认识一直存在争议。在构造背景分析的基础上,结合岩芯、测井及地震资料,对研究区梧桐沟组层序地层划分及地层沉积样式研究表明,梧桐沟组沉积时期构造强度较弱,地形较为平缓,地层在全区稳定分布;中晚二叠世盆地发生造山运动,吉木萨尔凹陷东南边缘区域经过构造抬升,上部地层遭受不同程度的剥蚀,形成了现今"底平顶削"的地层样式。以凹陷中部少数地层保存较全的井的地层叠加样式分析为基础,通过井-震结合识别不同类型界面,特别是最大湖泛面,建立了区域层序地层格架,即梧桐沟组为一个完整的长期基准面旋回(三级层序),并在其内部识别出5个中期旋回。梧桐沟组地层在不同的层序发育时期表现出不同的旋回叠加样式:下段沉积时期,即最大湖泛面以下,随着可容纳空间的增大,地层表现为明显的退积叠加样式;上段沉积时期,即最大湖泛面以上,随着可容纳空间的减小,地层表现为明显的进积叠加样式,符合沉积物体积分配原理。     

岩溶坝基帷幕灌浆高注浆量控制措施与效果分析

帷幕灌浆是坝基防渗处理的一种常用技术手段,而岩溶地区因其岩溶发育致使岩溶裂隙和通道纵横交错,且难以获取其规律性,灌浆过程中往往会出现注浆量偏高的现象。为了能够既有效控制帷幕灌浆成本,又确保帷幕灌浆施工质量,就需要采用科学、合理、可靠的施工技术措施和手段,对高注浆量的情况进行重点控制和改进。本文结合某水电站岩溶坝基帷幕灌浆中突出的高注浆量现象,通过有针对性地从技术措施和控制手段等方面的综合探索与应用,最终达到了坝基防渗的目的,并最大限度地节约了帷幕灌浆的时间,有效控制了高材料消耗所带来的高成本问题,为类似项目的帷幕灌浆积累了一定的经验。