MK模型的拟合方程和预测方程在土-水特征曲线中的应用及参数敏感性分析

土-水特征曲线可以预测非饱和土的各种性质（如：非饱和渗透系数、剪应力和热学性能等）。但测量土-水特征曲线耗时久且花费昂贵。为了解决这一问题,目前,很多研究都致力于从基本的岩土工程性质预测土-水特征曲线。基于此,以MK（Modified Kovács）模型的2种形式（拟合方程和预测方程）为土-水特征曲线模型,以Matlab编程语言中的cftool为拟合工具,以西宁黄土、粉砂土、红黏土和冰碛土4种细粒土为研究对象,对比拟合方程和预测方程描述细粒土土-水特征曲线的效果和差异,分析MK模型中黏附饱和度 1解 的变化规律,提出了基于MK模型的饱和度进行参数敏感性分析的计算公式。结果表明：拟合曲线和预测曲线在描述4类典型细粒土土-水特征曲线时均具有较好的效果,但拟合曲线整体上优于预测曲线;土壤质地和黏粒含量影响 值;饱和度对拟合参数 的敏感性较大,对拟合参数 的敏感性较小。     

粘性土-砂井地基单井固结模型的改进

单井固结模型的计算中通常将砂井周围土体简单划分为涂抹区和非涂抹区,不符合实际砂井周围土体的渗透系数分布复杂的事实。本文在Terzaghi固结理论的基础上提出了改进的单井固结模型,以一个待定参数流量系数Cq取代涂抹区和非涂抹区渗透系数来刻画砂井周围土体的横向渗透性特征,使单井固结问题得到高度简化又不失严密性。本文将改进模型用于非完整砂井单井固结的最终沉降量的数值计算,并将计算结果与谢康和改进法以及Hart法的解析解进行了比较,证明了改进模型数值解的可靠性。     

物理细胞自动机与岩石弹-脆-塑性性质的细观机制研究

基于能量守恒定律和岩石的基本力学性质,进一步发展了由作者提出的一种用于模拟岩石非线性破坏演化的新方法-物理细胞自动机（PCA）模型。该模型通过岩石内部（或细观）基元（或细胞）间简单的随机相互作用的综合结果来反映岩石系统整体的稳定宏观力学现象。利用PCA模型,研究了形成不同岩石本构关系的本质影响因素,揭示了岩石弹-脆-塑性性质的细观机制,为进一步认知岩石等非均质材料的力学性质提供了一种新的理论方法。同时,其研究思路和结论也可为微观和细观力学的数值模拟方法及新型复合材料的设计提供重要的借鉴。     

热-水-应力耦合影响的有限元分析

为了考虑应力拉压和压力（化学）溶解对裂隙开度的综合影响,对所建立的双重孔隙-裂隙介质热-水-应力耦合模型中裂隙开度的计算模型作了改进。通过一个假定的高放废物地质处置库算例,就岩体裂隙开度变化的3种工况,分析了岩体中的温度、孔（裂）隙水压力、地下水流速和主应力的变化、分布情况。结果显示：3种工况的计算域中温度场基本相同;孔（裂）隙水渗流场形态相似,但其量值有一定差别;工况1的裂隙开度在应力和压力（化学）溶解的共同作用下闭合量最大,负孔（裂）隙水压力增值最高;核废物的释热效应明显地改变了岩体自重应力场的水平分量,但对其垂直分量影响较小。     

水源热泵井施工中水文地质参数确定的探讨——以东煤沈阳钻探机械研制中心厂房水源热泵井施工为例

以单井抽水试验资料为基础,计算出该区的新近系砂砾岩含水层的出水量为40m^3/h。为解决水源热泵系统的热源及同层回灌问题,进行了较长时间的带有观测孔的抽水试验及回灌试验。试验表明：该区的渗透系数为32.16m/d,影响半径为166m,回灌渗透系数为14.38m/d,回灌影响半径为61m。据此计算干扰井单井抽水量为38.51m^3/h,干扰井单井回水量为11.14m^3/h,同时确定该区的抽灌井比例为1∶4。根据用水需求,该区设计抽水两口,回灌井8口,备用一口回灌井,形成该区的水源热泵系统。系统运行近两年来,单井出水量、回灌量及水位降深均符合设计要求,达到了预期效果。     

廊坊地区浅层地热能开发前景及地下水回灌试验研究

本文在对廊坊市万庄镇的浅层地热资源量进行概算的基础上,对水源热泵在该区的推广可行性做了评价;并且通过野外场地试验,对该区地下水供给能力及回灌能力进行了研究。得出如下主要结论:(1)该区域浅层地热能储量丰富,总储量约为18 075.0×1010KJ,相当于标准煤616.8×104t,具备推广水源热泵的基础条件。(2)研究区浅层含水层富水性较好,水质较好、供水能力较强,满足水源热泵的运行,也利于地下水回灌,不引起水资源的浪费。(3)在不引起次生环境地质问题前提下,建议在水源热泵系统运行时单井抽水量控制在80 m3/h左右,单井回灌量不超过40 m3/h,采用一抽两灌的原则配置抽、灌井数量,抽水井和回灌井的井间距控制在70~80 m,这样能够使得被抽取的地下水基本得到回灌补给,也可合理的利用浅层地热资源。     

考虑初始孔隙比影响的单/双峰土-水特征曲线模型研究

土样内部的孔隙结构可以分为单峰、双峰甚至多峰孔径分布。双峰孔隙结构土体的土-水特征曲线（SWCC）通常为双峰形态。现有的双峰SWCC模型大部分为单峰SWCC模型经过分段或叠加的方式所得,两段SWCC交点处的相关参数往往难以准确确定,而且很少有模型考虑孔隙比对双峰土水特性的影响。首先,详细分析和论述了广吸力范围内不同初始孔隙比土试样单/双峰SWCC的特点。其次,基于上述土水特性的分析,以及吸附水和毛细水两种不同的作用机制,结合开尔文方程推导出指数函数形式的吸附项方程,毛细项方程则以修正的Fredlund和Xing SWCC模型为基础。最后,提出一个考虑初始孔隙比影响的单/双峰SWCC模型,并利用不同类型土单/双峰SWCC的实测数据验证了模型的有效性。     

中国干旱-半干旱区表土bGDGTs与气候环境因子数据再分析

文章通过收集中国北方干旱-半干旱区域已发表的bGDGTs数据,重新计算了不同地区表土bGDGTs各组分与温度及pH等环境因子之间的相关关系。主要获得以下认识：1)不同于干旱区,半干旱区bGDGTs的异构体指标与土壤pH表现为强相关,在4个异构体指标(IBT、 CBT′、 IR′_IIa、 IR′_IIIa)中CBT′、 IR′_IIIa计算的pH结果精度较高;2)半干旱区f（5 me）与pH表现为强相关（-0.847）,指示了基于5-甲基bGDGTs的指标同样具有重建半干旱区pH的潜力。但是目前缺乏CBT′_5ME-pH之间的转换函数;3)全球温度校正方程计算得出的温度误差普遍大于中国区域校正方程（中国北方、青藏高原等区域）获取的温度误差。通过整合数据重新计算发现4个区域温度校正方程(MBT′_5ME、 Tmr、 MAATsr、 Index 1)计算的温度结果具有较小的误差,且在半干旱区重建的温度结果更为精确;另外MAATsr、 Index 1的区域温度校正方程也适用于干旱区的温度重建,...     

Green-Ampt模型渗透系数取值方法研究

格林-安姆普特（Green-Ampt）模型原理简单、使用方便,在浅层滑坡的降雨入渗分析中有很大的应用潜力。渗透系数是Green-Ampt模型中的一个主要参数,该渗透系数并不一定等于土体的饱和渗透系数。通过与理查茲（Richards）方程求解进行比较,研究了Green-Ampt模型中渗透系数的取值方法。研究发现,为获得与Richards方程相同的入渗量计算结果,需对饱和渗透系数进行修正,该修正系数的大小与入渗深度和土体类型有关。对于文中研究的土体,当Green-Ampt模型中渗透系数取为0.7倍饱和渗透系数时,由Green-Ampt模型计算的孔隙水压力分布与Richards方程计算结果最为接近,建议Green-Ampt模型中渗透系数修正系数取0.7。     

岩石损伤过程中的热-流-力耦合模型及其应用初探

岩石损伤过程热-流-力（THM）耦合问题的研究对于深部采矿等许多工程领域都具有重要的理论意义。以岩石的损伤为主线,在多场耦合分析方程中引入损伤变量,基于质量守恒和能量守恒原理,提出岩体损伤过程中的THM耦合模型。通过把均匀弹性介质THM耦合响应的模拟结果与理论分析结果进行对比,验证了程序及有限元实施的正确性。然后,用该耦合模型进行了不同地应力条件下流固耦合过程的数值模拟,探讨了水压力对于岩石损伤过程的作用机制。数值模拟表明,水压力导致了拉伸损伤范围的扩大和损伤程度的加剧,同时亦对剪切损伤具有抑制作用。     

不连续面在双重介质热-水-力三维耦合分析中的有限元数值实现

在建立双重介质热-水-力耦合微分控制方程的基础上,提出了裂隙岩体热-水-力耦合的三维力学模型,对不同介质分别建立以节点位移、水压力和温度为求解量的三维有限元格式,开发了双重介质热-水-力耦合分析的的三维有限元计算程序,在有限元数值分析中不连续面应力计算采用等厚度空间8节点节理单元进行离散,而不连续面渗流和热能计算时采用平面4节点等参单元进行离散,这样保证了不同介质之间的水量、热量交换和两类模型接触处节点水头、温度和位移相等。通过高温岩体地热开发算例,揭示了在热-水-力耦合作用下不连续面处于低应力区,其张开度随运行时间的延长呈非线性增加,非稳定渗流阶段不连续面显著地控制着渗流场的整体分布,它的水头远高于拟连续岩体介质的水头,而进入稳定渗流阶段不连续面的控渗作用不明显,由于高温岩体地热开发系统中存在大规模的热量补给,不连续面对岩体温度场分布的影响并不显著。     

A coupled model for heating and hydratation in unsaturated clays

Knowledge of transport processes of heat and moisture in soils of arid zones is vital to understanding the environmental and economic impacts of many activities: agriculture, waste disposal, geoenvironmental practices and earth sciences. Through extensive review and study on the different aspects of coupled transfer processes in swelling porous media, a general mathematical model for coupled heat, moisture, air flow and deformation problems in clayey soils is proposed in a consistent and unified manner. The model is characterized by the presence of a deformable solid matrix filled with two fluid phases (liquid water and air). In the proposed model, both pore water and air transfers are assumed to be governed by the generalized Darcy’s law. Fully coupled, non-linear partial differential equations are established and then solved by using a Galerkin weighted residual approach in space domain and an implicit integrating scheme in time domain. The obtained model has been finally validated by means of some case tests for the prediction of the thermo-hydro-mechanical behaviour of unsaturated swelling soils. The calculated relative errors between experimental and numerical results are 3% for temperature and 7% for stresses. Consequently, the developed numerical model predicts satisfactory results, compared to experimental test measures. The model is applicable to two-dimensional problems with various initial and boundary conditions; non-linear soil parameters can be easily included in this model.     

爆轰荷载作用下球空腔热流固耦合动力响应

基于饱和多孔弹性介质热流固耦合动力模型(THMD)及控制方程,研究了爆轰荷载作用下球空腔内壁面受到随时间变化热、力冲击作用下的热流固耦合动力响应。利用Laplace变换技术,对控制方程进行解耦,获得了温度、位移、孔隙水压力和应力积分形式解。采用Laplace逆变换得到数值结果,分析了热流固耦合参数的影响特性,并将热流固耦合结果与热弹性条件下的结果进行比较,以验证热流固耦合条件下结果的正确性。     

黏土岩温度-渗流-应力耦合特性试验与本构模型研究进展

高放废物处置库、垃圾填埋场等工程中常常涉及到温度场（T）、渗流场（H）和应力场（M）的耦合作用的问题。从试验和理论模型两个角度综述国内外黏土岩温度-渗流-应力耦合特性的研究进展,主要包括其传热特性、温度影响下的渗流特性、变形、强度、蠕变特性。在此基础上,重点分析了黏土岩水-热迁移模型以及热-力耦合本构模型的适应性。基于上述认识,通过试验研究了比利时Boom clay在温度作用下的强度、渗透性、蠕变性等特征。结果表明：随着温度升高,Boom clay的强度有所降低,渗透性显著增强,蠕变速率明显加快。提出了适用于Boom clay的THM耦合弹塑性损伤模型,计算结果验证了模型能合理反映温度的影响。最后,探讨了黏土岩THM耦合机理研究的不足和今后的研究方向。     

Thermo-hydro-mechanical modeling of fracturing porous media with two-phase fluid flow using X-FEM technique

In this paper, a fully coupled thermo-hydro-mechanical model is presented for two-phase fluid flow and heat transfer in fractured/fracturing porous media using the extended finite element method. In the fractured porous medium, the traction, heat, and mass transfer between the fracture space and the surrounding media are coupled. The wetting and nonwetting fluid phases are water and gas, which are assumed to be immiscible, and no phase-change is considered. The system of coupled equations consists of the linear momentum balance of solid phase, wetting and nonwetting fluid continuities, and thermal energy conservation. The main variables used to solve the system of equations are solid phase displacement, wetting fluid pressure, capillary pressure, and temperature. The fracture is assumed to impose the strong discontinuity in the displacement field and weak discontinuities in the fluid pressure, capillary pressure, and temperature fields. The mode I fracture propagation is employed using a cohesive fracture model. Finally, several numerical examples are solved to illustrate the capability of the proposed computational algorithm. It is shown that the effect of thermal expansion on the effective stress can influence the rate of fracture propagation and the injection pressure in hydraulic fracturing process. Moreover, the effect of thermal loading is investigated properly on fracture opening and fluids flow in unsaturated porous media, and the convective heat transfer within the fracture is captured successfully. It is shown how the proposed computational model is capable of modeling the fully coupled thermal fracture propagation in unsaturated porous media.     

The solar dynamo and integrated irradiance variations in the course of the 11-year cycle

Abstractthe effect of the large-scale magnetic fields generated by the solar dynamo on the radiation flux issuing from the convection zone is studied. A governing equation describing convective heat transfer is obtained in the framework of mean-field magnetohydrodynamics, with account for the influence of magnetic fields and differential rotation on the energy budget of the convection zone. The principal effects are illustrated using a one-dimensional numerical model. Calculations indicate that the influence of large-scale magnetic fields can modulate the solar irradiance with a relative amplitude of ～0.07%.     

Energy equation and stress–dilatancy relationship for sand

The energy equation is an expression of the first law of thermodynamics or the law of conservation of energy. According to the first law of thermodynamics, the externally applied work to a system is equal to the sum of dissipation energy and Helmholtz free energy of the system. However, most of the currently available stress–dilatancy relationships are based on the energy equation of Taylor-Cam Clay type, which hypothesizes that the applied plastic work is equal solely to the frictional dissipation energy. The Helmholtz free energy has been completely neglected. Recently, observed from acoustic experiments, it has been recognized that Helmholtz free energy can be caused by deformation mechanisms other than friction between particles. Thus, it is necessary to include additional terms in the energy equation in order to correctly model the stress-dilatancy behavior. This paper addresses the issue regarding the balance of this energy equation. Analyses of experimental results are presented. Specific forms of the frictional energy and Helmholtz free energy are proposed. The proposed energy equation is verified with the experimental data obtained from Silica sand, Ottawa sand, and Nevada sand.

     

Modelling of internal erosion based on mixture theory: General framework and a case study of soil suffusion

A general thermo-hydro-mechanical framework for the modelling of internal erosion is proposed based on the theory of mixtures applied to two-phase porous media. The erodible soil is partitioned in two phases: one solid phase and one fluid phase. The solid phase is composed of nonerodible grains and erodible particles. The fluid phase is composed of water and fluidized particles. Within the fluid phase, species diffuse. Across phases, species transfer. The modelling of internal erosion is contributed directly by mass transfer from the solid phase towards the fluid phase. The constitutive relations governing the thermomechanical behaviour, generalised diffusion, and transfer are structured by the dissipation inequality. The particular case of soil suffusion is investigated with a focus on constitutive laws. A new constitutive law for suffusion is constructed based on thermodynamic conditions and experimental investigations. This erosion law is linearly related to the power of seepage flow and to the erosion resistance index. Owing to its simplicity, this law tackles the overall trend of the suffusion process and permits the formulation of an analytical solution. This new model is then applied to simulate laboratory experiments, by both analytical and numerical methods. The comparison shows that the newly developed model, which is theoretically consistent, can reproduce correctly the overall trend of the cumulated eroded mass when the permeability evolution is small. In addition, the results are provided for four different materials, two different specimen sizes, and various hydraulic loading paths to demonstrate the applicability of the new proposed law.     

热-水-应力耦合弹塑性二维有限元程序的开发与应用尝试

从已有的应力平衡方程、水连续性方程、能量守恒方程和弹塑性矩阵入手,使用Galerkin方法将各控制方程分别在空间域和时间域进行离散,初步开发出了一个用于分析饱和岩土介质中热-水-应力耦合弹塑性问题的二维有限元程序,并通过2个算例进行了应用尝试。     

Finite element analyses of rate-dependent thermo-hydro-mechanical behaviors of clayey soils based on thermodynamics

Clayey soils in the vicinity of energy geostructures may be exposed to long-term periodic thermal cycles. The creep and consolidation behaviors of the clayey soils can be both rate-dependent and temperature-dependent, and the underlying physical mechanisms are merely investigated theoretically. In this study, based on the theory of thermodynamics, a fully coupled thermo-hydro-mechanical (THM) finite element (FE) program for saturated soils is developed for this purpose. The FE formulation accounts for the combined effect of rate and temperature through the novel concept of granular temperature. Simulations of THM coupled validation cases and a series of experimental observations on the soft Bangkok clay are carried out. The obtained numerical results exhibit good agreement with analytical solutions and laboratory measurements. It is found that three fundamental physical mechanisms contribute to the irreversible thermal contraction observed for normally consolidated and lightly overconsolidated clays under drained thermal cycles: (1) the thermal creep excited by mass exchange from adsorbed water to free water; (2) the mechanical creep induced by confining stresses; and (3) the increase in granular packing caused by the thermal expansion of soil particles. The thermal contraction generally stabilizes within a few thermal cycles, as a result of the noticeable reduction in the thermal creep rate. It is further demonstrated that the transient heat transfer and the heating rate can greatly influence the deformation of clays subjected to thermal cycles.