寒区大坝心墙土料冬季冻融与防控监测

两河口水电站是雅砻江干流规划开发中的中游控制性龙头电站,大坝为300米级砾石土心墙堆石坝。由于地处川西高原气候区,冬季气候寒冷干燥,大坝心墙土料填筑过程中面临着冻融问题的困扰。基于一个完整冬季的现场监测,系统分析了这一寒区大坝心墙建筑过程中砾石土、接触黏土土料温度变化规律、冻融特征与影响因素以及现有保温措施防冻效果。结果表明,无保温措施条件下,砾石土、接触黏土均出现了负温冻结现象,其中砾石土最大冻结深度达20 cm,接触黏土达14 cm,土料冻结持续时间不超过1个昼夜,为短时冻土。土料降温冻结过程以与外界大气对流换热过程为主,受气温、风速条件影响显著,波动范围较大,而升温融化过程以太阳辐射增温过程为主,与有效辐射起始时间密切相关,因此波动范围较小。监测期内,采用三布两膜保温材料覆盖可有效防止心墙土料冻结的发生,有、无保温材料覆盖条件下浅层土料最低温度可相差约5℃。结合监测结果和现场实际,建议将现有心墙防渗土料半幅填筑方式转换为全幅填筑,并研发轻质、憎水保温材料及其快速收放机械设备,以提高填筑作业进度、强化防冻效果。系统完整的现场实测结果,可为未来寒区水电大坝建设提供基础数据和参考。     

Tectonics and Chemical Transport: The Rules of the Resource Game

The most basic solid earth resources for use by man include construction materials, energy resources, metals, ceramics, soils. With the exception of construction materials, the formation of economic mineral deposits and soils involves the thermal or gravity driven flow of various types of geofluids. The search for mineral deposits thus involves the understanding of the exact nature of fluids, their mass and flow patterns. Given the giant thermal disturbances associated with the formation and destruction of Gondwana, it must be no surprise that these regions contain very significant mineral resources on a large scale.     

块石路基对流特性实验研究

利用块石层调控冻土路基地温是冻土工程主动冷却路基的一种重要措施,并在青藏铁路重大工程建设中得到广泛应用.通过块石层中空气流动特性的深入研究,对其结构优化和效能提升具有重要意义,而不同温差条件下块石层中空气的流动特性,以及微风速的实际测定一直是其中亟待解决的问题.通过高精度微风速探测,首次获得封闭条件下块石层的空气流动特征.结果表明:在边界温度波动条件下,块石层内的降温过程与其内部自然对流过程密切相关,在对流作用下快速完成降温过程;升温过程则在相对较长时段内主要由内部的热传导换热完成,由此导致温度变化曲线非对称这一特殊现象的出现,并随深度增加而愈加明显.块石层上下端面的温差是对流发生及强度的关键控制因素,当温差达到一定量值时流速随温差的增加而加大,同时空气温度的变化也会对空气流动过程产生影响;块石层内整个空气流动过程是一个传递过程,也是热量传递的过程.该实验研究结果将对该种工程措施机理的进一步认识和改进、相关模拟计算参数的选取,以及在冻土区高等级公路等重大工程建设中的进一步有效应用均具有重要意义.     

非饱和膨润土的有效热传导特性模型

膨润土缓冲材料热传导特性的研究,对于高放废物深地质处置系统的安全评价至关重要。基于串、并联原理,通过将土体孔隙划分为与固相基质并联和串联两部分,提出了考虑矿物成分、颗粒亲水性、孔隙率及饱和度等因素的非饱和膨润土有效热传导系数的4种预测形式,建立了基于4种形式线性组合的有效热传导特性预测模型。详细讨论了模型参数的确定方法,并讨论了孔隙率、饱和度和孔隙结构、颗粒亲水性等因素对土体有效热传导特性的影响。基于MX-80膨润土和高庙子膨润土热传导特性试验成果,对模型的预测性能进行了验证。结果表明,由于膨润土颗粒尺寸较小且具有亲水特性,孔隙内的空气与水宜采用并联描述。研究成果对于非饱和膨润土的导热性能以及工程屏障系统的THM耦合数值模拟研究具有一定的参考价值。     

考虑残余含气量的非饱和土的水力耦合本构模型

大量的非饱和土干湿循环试验表明,当土体处于吸湿过程直至吸力降低为0 kPa时,土体并不能达到完全饱和状态,还存在一定的残余气体。在高饱和度时,由于残余气体以封闭气泡的形式分布在土体中,土体呈现较大的压缩性,使其与饱和土的性质不同。在这种状态下,现有的非饱和土本构模型预测到的土饱和度为100%,与试验结果存在一定的偏差。为了使本构模型在高饱和度状态时具有较高的精度,对非饱和土的毛细滞回和塑性变形耦合本构模型进行了修正,使其能够考虑残余含气量的影响。通过预测与实测结果比较,证明了新模型能够有效地模拟残余含气量对非饱和土力学特性的影响。     

Sources of dissolved oxygen in monitoring and pumping wells

Groundwater monitoring and pumping wells set in anoxic aquifers require attention to keep the groundwater free of dissolved oxygen (DO). In properly constructed monitoring or pumping wells, two processes can however still introduce oxygen to anoxic groundwater: (1) permeation of oxygen through polymer materials such as silicone, PVC, HDPE or Teflon, and (2) thermally driven convection, which can occur in all types of piezometers or wells, regardless of construction material, when the water table or pressure head is close (<10 m) to the land surface. Here, field measurements (temperature and DO well loggings) from a monitoring well in Bilthoven, the Netherlands, are combined with analytical and numerical modelling to investigate the role of both processes on oxygenation of anoxic groundwater in wells. The results of numerical and analytical modeling show that both permeation and convection can introduce oxygen into anoxic wells to near saturation concentrations. In the field data gathered, convection is primarily responsible for oxygen intrusion up to a depth of around 12 m. Oxygen intrusion through convection and permeation in monitoring and pumping wells may influence groundwater sampling and analyses, and may contribute to well clogging, depending on site conditions. The combination of field and modelling provides new insights into these processes, which can be used for both groundwater sampling and pumping well design.     

Application of Relevance Vector Machine for Prediction of Ultimate Capacity of Driven Piles in Cohesionless Soils

This paper examines the potential of relevance vector machine (RVM) in prediction of ultimate capacity of driven piles in cohesionless soils. RVM is a Bayesian framework for regression and classification with analogous sparsity properties to the support vector machine (SVM). In this study, RVM has been used as a regression tool. It can be seen as a probabilistic version of SVM. In this study, RVM model outperforms the artificial neural network (ANN) model based on root-mean-square-error (RMSE) and mean-absolute-error (MAE) performance criteria. It also estimates the prediction variance. An equation has been developed for the prediction of ultimate capacity of driven piles in cohesionless soils based on the RVM model. The results show that the RVM model has the potential to be a practical tool for the prediction of ultimate capacity of driven piles in cohesionless soils.     

土壤温室气体昼夜变化及其环境影响因素研究

通过对北京东灵山草地和桦树林土壤气体CO₂,N₂O和CH₄浓度及其排放通量的昼夜连续观测,探讨了生长季节草地和森林土壤温室气体昼夜变化及其环境影响因素。研究表明:1)土壤CO₂排放通量昼高夜低,N₂O排放通量有明显小时尺度波动,但昼夜变化不突出;土壤CO₂和N₂O浓度昼夜变化不明显,且与排放通量波动不一致;土壤是大气CH₄的一个汇,相对厌氧的环境可能有利于土壤吸收CH₄。2)无雨时气温昼夜变化通过影响土壤表层的气体扩散和CO₂产生过程,来影响土壤CO₂和N₂O的地表排放通量,而对土壤10cm以下CO₂和N₂O的产生影响不大。小时尺度的土壤CO₂和N₂O浓度波动则可能还有其他影响因素或机制。3)降雨时土壤渗水引起的土壤空气对流取代气体浓度扩散成为土壤与大气空气交换的主要方式,导致土壤CO₂和N₂O排放通量的同步波动。降雨渗水较多时,较多的溶解氧随着雨水进入土壤内,会促进土壤CO₂的生成和抑制N₂O的产生。4)土壤CO₂与N₂O浓度存在显著的正相关关系,反映出土壤CO₂和N₂O有相对稳定的产率比。土壤有效碳可能是造成土壤CO₂与N₂O浓度正相关的主要原因,土壤空气的氧分压则可能是造成土壤CO₂和N₂O浓度波动不一致的重要因素。     

古玉器的化学风化机理初探之一——粉末态模拟实验

矿物的溶解和分解是矿物风化的重要过程,因此探明矿物的溶解和分解过程将有助于揭示其风化机理。古玉在长期的埋藏过程中,处于土壤及土壤水的溶解反应体系中。本文初步模拟古玉埋藏的酸、碱性土壤环境,通过对不同材质玉料的粉末形态在酸、碱缓冲溶液中的浸泡实验,初步探讨了不同材质玉料在酸、碱环境下的风化机理。     

透壁式通风管-块石复合气冷路基室内模型试验研究

为了研究透壁式通风管-块石复合气冷路基的降温效果,针对年均气温约-3.5℃,平均风速2.5 m·s^-1,主导风向为西北方向的高原环境条件开展室内模型试验,分析了单一块石气冷路基和透壁式通风管-块石复合气冷路基的孔隙空气对流速度、特征点地温变化过程以及模型整体温度场变化过程.结果表明：在透壁式通风管的疏导作用下,透壁式通风管与块石层的复合结构能够起到强化路基体对流的效果,复合路基块石孔隙中的空气流速比单一块石路基提高约20%,由此导致复合路基模型底部的降温幅度是单一块石路基模型的2.2倍.模型整体温度场表明,复合路基能够起到储存冷量、降低下伏多年冻土地温的作用.     

Experimental simulation of magma–carbonate interaction beneath Mt. Vesuvius, Italy

We simulated the process of magma–carbonate interaction beneath Mt. Vesuvius in short duration piston-cylinder experiments under controlled magmatic conditions (from 0 to 300 s at 0.5 GPa and 1,200 °C), using a Vesuvius shoshonite composition and upper crustal limestone and dolostone as starting materials. Backscattered electron images and chemical analysis (major and trace elements and Sr isotopes) of sequential experimental products allow us to identify the textural and chemical evolution of carbonated products during the assimilation process. We demonstrate that melt–carbonate interaction can be extremely fast (minutes), and results in dynamic contamination of the host melt with respect to Ca, Mg and ⁸⁷Sr/⁸⁶Sr, coupled with intense CO₂ vesiculation at the melt–carbonate interface. Binary mixing between carbonate and uncontaminated melt cannot explain the geochemical variations of the experimental charges in full and convection and diffusion likely also operated in the charges. Physical mixing and mingling driven by exsolving volatiles seems to be a key process to promote melt homogenisation. Our results reinforce hypotheses that magma–carbonate interaction is a relevant and ongoing process at Mt. Vesuvius and one that may operate not only on a geological, but on a human timescale.     

地幔对流的数值模拟方法

地幔对流的数值模拟研究为理解地壳运动、地幔内部的性质和地球的动力学机制提供了手段。笔者回顾了地幔对流数值模拟的发展,总结了各种常用的描述地幔对流的数学模型,边界条件以及数值解法。分析了级数展开方法,有限差分方法,有限元方法和多重网格方法的性质和特点,针对某些方法编写了二维地幔对流的计算机软件,使用了实际的地质参数,无穷大的Ｐｒａｎｄｔｌ数,在１０４到１０７之间的Ｒａｙｌｅｉｇｈ数,并给出了数值结果和地质解释。     

三相非饱和松散颗粒介质中弹性波的弥散性分析

由已建立了的非饱和土这种三相松散介质中的动力控制方程,得到波动方程的解,表明在非饱和土中存在3种P波（P1、P2和P3波）和1种S波。根据推导出的波速以及衰减的解析表达式,以实际工程中的尾砂土这种非饱和多孔介质为例,通过Matlab计算工具对非饱和土中的4种体波的弥散特性进一步数值分析,给出了典型的弥散和衰减曲线。结果表明,在一定的频率范围内,4种波的波速及衰减随频率的升高而增大,其中P2、P3波受频率影响强烈,而P1、S波受频率影响则相对较小。此外,饱和度的增加也将导致P1波波速的增大,但对P2、P3波和S波影响不大。     

透壁通风管路堤的对流换热和蒸发散热

针对透壁通风管路堤中透壁通风管管壁与空气之间的对流换热和土体水分通过管壁小孔的蒸发散热机制,分析了开孔率、风速及含水率等因素对透壁通风管管壁对流换热和水分蒸发散热的影响,并具体给出了管壁对流换热系数和蒸发散热系数的计算公式。冬季路堤由于冻土层未冻水含量较小而使管壁小孔的水分蒸发散热较弱,路堤总的降温效果主要由管壁对流换热效应控制,而暖季通风管内空气与管壁的对流换热效应可使路堤土体增温,同时,由于通风管周围融土的未冻水含量较大,而使得通过管壁小孔的水分蒸发散热较强,可部分或全部抵消对流换热引起的增温效应,而有利于路堤的稳定。     

非饱和土二维固结简化计算的研究

针对较高饱和度的非饱和土二维固结问题展开研究。对于较高饱和度的非饱和土,如饱和度 时,将孔隙中气、水近似地看作可压缩的气、水混合流体后,非饱和土近似为土骨架和混合流体的二相土。考虑混合流体的压缩性,建立混合流体的连续方程。联立平衡方程和混合流体的连续方程,求出应力-应变和混合流体压力;再建立水连续方程求解水压力,继而求出气压力、吸力等。算例表明：加载和消散过程中,混合流体与水压力变化基本一致,气压的作用并不大;地基变形过程与高速公路填筑过程中地基变形发生规律一致。说明该简化方法是合理的,并促进了非饱和土固结变形计算走向实用化。     

On the Validity of the Boussinesq Approximation for the Elder Problem

The Elder problem is an example of a density driven flow, motivated by experiments of a thermally driven convection in porous media. It is a mathematical benchmark problem used for code verification of density driven flow simulators and comparison of issues related to its numerical treatment. Its bifurcation diagram with respect to the Rayleigh number is investigated on a hierarchy of uniformly refined grids. Eleven stationary solutions are shown to exist for the Elder problem. Similar solutions can be found using the Boussinesq approximation. Despite this similarity the corresponding bifurcation diagrams are shown to be topologically not equivalent. This gives rise to serious doubts on the validity of Boussinesq approximation for this model problem. Grid convergence is investigated for the numerically obtained solutions.     

软土中已打入桩对土体位移的遮帘作用研究

采用侧向移动土中桩的理论分析结果以及有限元数值模拟,较为系统地分析了沉桩工程中已打入桩对土体侧向位移及地表隆起的遮帘作用,结果表明已打入桩对土体侧向位移的总体遮帘作用并不显著。采用不同切向劲度系数的接触面单元以模拟已打入桩间隔时间长短对地表隆起的遮帘效应,结果表明：打入桩间隔期较短时,已打入桩对地表隆起的遮拦作用主要体现在对土体位移的“水平阻隔”上,土体沿桩土界面发生流动向上隆起;而间隔期较长时,桩土间存在有较大的摩阻力而限制了土体的位移,此时已打入桩对地表隆起的遮拦作用体现在竖向摩阻力对土体隆起的“抑制”上。     

碎(块)石路堤孔隙空气对流运动的Darcy定律适用性

针对青藏高原多年冻土区路基工程中采用的碎石路堤地温调控技术,分析了碎石层孔隙空气对流运动符合Darcy定律的适用范围及其可能偏离Darcy定律的非线性效应,并给出了考虑非线性效应的Darcy定律修正方程.同时,讨论了适合于路堤碎石层的空气渗透系数及其表征空气流动非线性效应的Ergun常数的估算公式.     

The potential for free and mixed convection in sedimentary basins

 Free thermal convection and mixed convection are considered as potential mechanisms for mass and heat transport in sedimentary basins. Mixed convection occurs when horizontal flows (forced convection) are superimposed on thermally driven flows. In cross section, mixed convection is characterized by convection cells that migrate laterally in the direction of forced convective flow. Two-dimensional finite-element simulations of variable-density groundwater flow and heat transport in a horizontal porous layer were performed to determine critical mean Rayleigh numbers for the onset of free convection, using both isothermal and semi-conductive boundaries. Additional simulations imposed a varying lateral fluid flux on the free-convection pattern. Results from these experiments indicate that forced convection becomes dominant, completely eliminating buoyancy-driven circulation, when the total forced-convection fluid flux exceeds the total flux possible due to free convection. Calculations of the thermal rock alteration index (RAI=q·∇T) delineate the patterns of potential diagenesis produced by fluid movement through temperature gradients. Free convection produces a distinct pattern of alternating positive and negative RAIs, whereas mixed convection produces a simpler layering of positive and negative values and in general less diagenetic alteration. Received, January 1999/Revised, June 1999/Accepted, July 1999     

Desiccation shrinkage of non‐clayey soils: multiphysics mechanisms and a microstructural model

Analysis of macroscopic desiccation shrinkage experiments indicates that most, but not all of the shrinkage during drying occurs while soil is still saturated. Shrinkage practically ceases and air starts to penetrate the soil, when the water content is still quite high, for example, above 20% for the tested soils. The remaining, unsaturated drying process occurs with a much‐reduced shrinkage rate. In this context, we examine data of the pore system evolution as represented by the mercury porosimetry experimental results. The process is then modeled as a two‐stage process of deformation and evacuation of a two‐tube vessel system driven by the external evaporation flux. In the first stage, Poiseuille flow occurs through the vessels. The amount of water evaporated in this stage equals to the reduction of volume of the vessel through the deformation of its walls. This stage ends when a negative water pressure (suction) required to further deform the vessel reaches a critical value at which air enters the pore space. Two physical interpretation of such threshold are discussed. In the subsequent stage, evaporation proceeds with a receding liquid/vapor interface starting from the open end, incrementally emptying the vessel but with a marginal water flow and vessel deformation. The leading variables of the process are identified, and a quantifiable multiphysics meso‐scale scenario of models is established. Copyright © 2012 John Wiley & Sons, Ltd.