赖远明 黄文熙讲座-冻土水热力特性分析新方法-2015-4-10

合集下载

冻结砂土的损伤试验研究_金龙

冻结砂土的损伤试验研究_金龙

文章编号:1000-0240(2008)02-0306-07冻结砂土的损伤试验研究 收稿日期:2007-06-11;修订日期:2007-11-12 基金项目:国家自然科学基金重点项目(40730736);中国科学院知识创新工程重要方向项目(KZCX3-SW -351)资助 作者简介:金龙(1982—),男,甘肃榆中人,2005年毕业于西安交通大学,现为硕士研究生,主要从事冻土力学的研究. E -mail :jlcoolmail @金龙1,赖远明1,高志华2,李双洋1(1.中国科学院寒区旱区环境与工程研究所冻土工程国家重点实验室,甘肃兰州 730000;2.长安大学建筑工程学院,陕西西安 710061)摘 要:任何数值的荷载,都会导致冻土内部冰体塑性流动和冰晶体在土体内的重新定向,以及土颗粒的重新排列.这几方面的因素会使冻土产生微结构的损伤,同时表现出明显的各向异性.基于损伤力学理论,推导出冻土的各向异性损伤变量的表达式,利用-6℃冻结砂土的常规三轴试验数据,获得损伤变量的变化曲线。

结果表明:损伤变量随应变的增大而增大,呈现出双曲线的变化规律;径向损伤大于轴向损伤;围压有助于冻土结构的强化,减少冻土结构的损伤量.这些规律为利用损伤力学理论建立冻土的损伤本构关系提供了依据.关键词:冻结砂土;损伤;横观各向同性;变化规律中图分类号:P642.14文献标识码:A0 引言 冻土是由岩土颗粒、冰、未冻水和空气组成复合介质,它与融土的本质区别是冻土中冰的存在,土冻结后强度将提高很多倍,其中主要依赖于冻土粘聚力的大幅度提高,而冻土的粘聚力主要受冰控制[1-2].对冻土施加荷载后,将导致内部冰体的塑性流动,冰晶体的重新定向,以及土颗粒的重新排列,从而使冰的粘聚力下降,在土体内将会产生微裂纹,随着荷载的增大,最终导致土体的破坏.这一过程,可以看作是冻土体内部的损伤产生和发展的过程. 损伤理论,是将固体物理学、材料强度理论和连续介质力学统一起来研究材料的本构关系,解释材料的破坏机理,预估其剩余寿命[3-4].国外一些学者不断完善发展了损伤力学理论[5-7],而目前国内外对于土的损伤理论研究较少,沈珠江[8]和赵锡宏等[9]做了一些这方面的研究.维亚洛夫等[10]在试验基础上引入了定性表征冻土微结构损伤的初步概念;苗天德等[11]在蠕变试验和微结构观测研究的基础上,选择了描述冻土微结构变化的损伤内变量,建立了内变量演化方程及冻土蠕变的一般损伤理论.H e Ping 等[12]基于各向同性损伤假设,依据热力学原理,并引用损伤力学理论,分析了冻土所具有的粘弹塑力学性质以及围压的强化弱化双重性质;李栋伟等[13-14]以冻土试验为基础,提出了冻粘土服从DP 屈服准则的粘塑性损伤变量,研究了冻土粘弹塑性损伤耦合本构关系. 在荷载的作用下,土体将会发生微结构的损伤[11],同时,会使冰晶体产生明显的各向异性,特别是在塑性状态[2].因此,如何正确的选用损伤变量,成为建立损伤本构模型的关键.本文试图通过分析冻结砂土的常规三轴试验数据,利用损伤力学理论推导出冻结砂土的各向异性损伤变量的表达式,做出其变化规律曲线,为以后利用损伤力学理论建立冻土的损伤本构关系提供依据.1 试验描述及分析 本试验选用的土样为青藏线砂土,其颗粒成分如表1.试验时将砂土击实到一定的干重度,制备成直径为61.8m m ,高度125mm 的试样,抽气饱和后置于低温环境下冻结48h .然后将冻结土样取第30卷 第2期2008年4月冰 川 冻 土JO U RN A L OF G L ACIO L OG Y A N D GEOC RYO LO G YV ol .30 N o .2Apr .2008出,并装到低温试验机上,在特定温度下恒温24h ,加围压固结5min ,然后施加轴向荷载.试样的平均含水量为12.7%,平均干容重为19.7kN ·m -3,温度控制精度为±0.1℃.围压从0.3M Pa 一直加到18M Pa ,剪切速率为2.08×10-2m m ·s-1.表1 试样的颗粒组成T able 1 P article fractio ns o f the soil粒径范围和百分比/%<0.075mm 0.075~0.1mm 0.1~0.25mm 0.25~0.5mm 0.5~1m m >1m m 6.24131.08154.414.6922.2410.869 对-6℃下的冻结砂土进行了不同围压下的常规三轴试验,部分试验结果如图1.图1 -6℃冻结砂土的应力应变图F ig .1 Stress -strain curves of thefro zen sandy soil a t -6℃ 由图1可以看出,当围压<3M Pa 时,具有明显的应变软化现象,当围压>3M Pa 时,应力应变关系表现为应变硬化现象;在应力应变的全过程中,弹性应变占很小部分,变形过程主要表现为塑性变形过程;在σ3<10MPa 时,应力峰值随围压的增大而增大,而在σ3>10M Pa 时,应力峰值随围压的增大而减小,这是因为,冻土在高围压下发生压融现象[12,15],即随着围压的增大,冻土中孔隙冰的压融和微裂隙的发育,造成了冻土强度的弱化.2 损伤变量的推导 Kachanov 在研究金属的蠕变断裂时,第一次提出用连续性变量Χ描述材料的损伤状态:Χ=A~A(1)式中:A 为试样初始面积;A ~为受损后试样的净面积或者有效面积. Rabotnov 引入一个与连续性变量相对应的变量D ,称为损伤变量:D =A -A~A(2)式中:D =0对应于无损伤状态;D =1对应于完全损伤状态;0<D <1对应于不同程度的损伤状态. 令σ=F /A 为横截面上的名义应力;σ~=F /A ~为有效截面上的应力,称为有效应力,由σ·A =σ~·A ~,得:σ~=σ1-D(3)根据Lemaitre 等效应力原理:应力σ作用在受损材料上引起的应变与等效应力σ~作用在无损材料上引起的应变等价,可得到一维问题中受损材料的应力应变关系:ε=σE~=σ~E =σE (1-D )(4)即σ=E ε(1-D )(5)式中:E ~为受损材料的弹性模量,即有效弹性模量;E 为无损材料的弹性模量,即初始弹性模量;D 为损伤变量. 由式(4)和式(5)可知,材料的损伤变量可以通过材料弹性模量的变化来测定,即:D =1-E~E(6) 式(6)得出来的损伤变量是由一维问题中推得的,描述的是各向同性损伤.而冻土在加载时,由于冰晶体的各向异性,所以实际的损伤应该是三维的,采用各向同性损伤往往不能正确描述冻土的损伤过程.在三维问题中,有效应力为:{σ~}=[M (D )]{σ}(7)式中:[M (D )]为损伤有效矩阵,损伤有效矩阵为:[M (D )]=diag {1/W 11,1/W 22,1/W 33,1/W 12,1/W 23,1/W 31}(8)式中:W ij =[(1-D i )(1-D j )]1/2,i ,j =1,2,3;D 1,D 2,D 3分别为主轴的损伤变量.冻土在加载过程中,横向平面内是各向同性的,但在垂直于平面的方向性质不同,表现为横观各向同性损伤,即D 2=D 3.3072期金龙等:冻结砂土的损伤试验研究 假设冻土在初始状态(无损状态)下为各向同性,无损土的柔度矩阵[C]-1为:[C]-1=1E-vE-vE0001E-vE0001E0002(1+v)E00S2(1+v)E02(1+v)E(9)式中:E为弹性模量;v为泊松比. 损伤土的有效柔度矩阵为:[C]-1=[M(D)]T[C]-1[M(D)](10)式中:有效柔度矩阵是一个6×6的对称矩阵,其各个分量为:C-1ii=1E(1-D i)2,C-1ij=-vE(1-D i)(1-D j),i,j=1,2,3;C-144=2(1+v)E(1-D1)(1-D2),C-155=2(1+v)E(1-D2)(1-D3),C-166=2(1+v)E(1-D3)(1-D1);其余分量均为0.又因为:{εe}=[C]-1{σ}(11) 由式(10)和式(11),得到横观各向同性损伤变量的表达式:D1=1-E~E12(12)D3=1-vV~(1-D1)(13)式中:V~为受损材料的有效泊松比.3 损伤试验分析 试验证明,加载初期冻土体就已存在微裂隙发育,如何量化的描述损伤变量和损伤门槛值是研究冻土损伤力学的关键.一些学者利用C T技术,对冻土的加载过程进行了动态监测,从而得到了冻土的损伤规律.刘增利等[16]对冻土单轴压缩进行了CT测试,给出了冻土单轴压缩过程中不同承载阶段细观结构损伤演化特征;Zhang Shujuan等[17]借助C T技术对动态条件下泥岩的冻融损伤及强度和变形特征进行了研究.由式(12)和式(13),我们设想从冻土弹性模量的变化,得到损伤变量和损伤门槛值.这样,通过常规三轴试验就可以实现,不需要借助其它辅助设备. 从图1可以看出,应力-应变曲线存在明显的转折点.在加载初期,应变变化较缓,随着应力增大,到达转折点后,应变变化较快。

冻土水热力特性分析新方法

冻土水热力特性分析新方法
6000 5000 4000 3000 2000 1000
0 -10 -8 -6 -4 -2 0 2
温度(℃)
水分迁移—分 凝冰—融化— 工程病害
高含冰量—上 限及温度变 化—工程病害
量变:升温—承载力或强度降低 质变:融化—承载力或强度丧失
一、研究背景
沉降变形
水热侵蚀
工程的修建改变了原有冻土环境的水热平衡状态,随着全球气 候变暖和多年冻土的退化,冻融灾害对寒区道路的破坏日益严 重,危及交通安全。
(1.2676r ' 1)w 2.2676 0(wr ' 142c0 )
盐渍土信号强度与温度关系
其中,θ0为初始含水量,c0为溶液初始浓度, r’为温度T’时 溶液的溶解度。
二、硫酸盐渍土冻结温度
2.7 未相变含水量与温度关系
对温度进行修正之后,盐晶体析出引起的未相变含水量
变化为:
(180c0 w ) 0
内陆盐渍土 沿海盐渍土
高原盐渍土
我国盐渍土面积约20多万平方公里,居世界第四。
一、研究背景
高原盐渍土
冻融作用与盐渍化 作用重叠、相互影 响:
➢伴随水热过程的 盐分迁移
➢影响路基工程热、 力、水文状况
➢路基稳定性问题
一、研究背景
强度(kPa)
温度变化—(冰强度、胶结力)—冻土强度变化 —(冰-水相变)—冻、融土力学特性激变
T*
0.0 200 210 220 230 240 250 260 270 280
Liw T f solution
T
Vs sl
dA dV
Vs表示冰晶的摩尔体积,ΔTγ表示由克服冰液表面自由能引 起的温度降低值。Tf-solution表示溶液的冻结温度。

开放边界条件下青藏铁路抛石路基的降温效果分析

开放边界条件下青藏铁路抛石路基的降温效果分析

1 控制方程与数学模型
目前, 在青藏铁路建设过程中, 已大量采用抛石
路基结构, 抛石层多为上下封闭、左右开放, 外界风 和气压变化直接影响抛石路基内的空气对流方式. 因此, 计算中建立了路基区与空气区的传热传质耦 合模型.
由于铁路道碴铺层和抛石层可以看作多孔介质,
其内部热对流方式为非稳态的非等温渗流, 方程组 为连续性方程、动量方程和能量方程[14,15]:
将图 1 所示路基模型按平面四节点单元划分, 采 用自适应网格技术, 即在计算中根据精度需要程序 对网格密度进行自动调整, 网格数目变化在 10~15 万 之间.
热初始条件按(11)式, 不考虑升温的天然地表温
度方程作为Ⅲ、Ⅳ区的上边界条件进行反复多年计 算, 直至得到稳定的温度场为止, 将此时获得的 7 月 15 日温度场作为该区计算的温度初始条件; 同时, 为 确保计算结果的有效性和实用性, 路堤区域(Ⅰ, Ⅱ, Ⅴ区)温度取该时间浅地表土层温度. 速度初始条件 取不考虑能量方程((4)和(8)式), 定温条件下以(15)式 为边界条件, 计算得到的 7 月 15 日稳定速度场作为 计算流体区速度初始条件.
考虑全球气候变暖的影响, 根据文献[18]取青藏 高原未来 50 年年平均气温上升 2.6℃, 设初始年平均 气温为−4.0℃, 由于路基模型上表面温度变化不仅与 环境空气有关, 而且受到太阳辐射等复杂因素的综 合 作 用 , 故 根 据 附 面 层 理 论 [19]及 长 期 现 场 观 测 资 料[20], 对计算域的热边界条件进行如下设定:
在计算中应用显热容法假设介质相变发生在wwwscichinacom580中国科学地球科学第35当建立等效体积热容时应考虑温度间隔的影响同时假设介质在正冻未冻时的体积热容c不取决于温度因此简化构造出9b由于求解模型中的温度场压力场和速度场是一个强非线性问题无法求得解析解故计算中在交叉网格上采用控制容积积分法来对方程进行离散后利用simple算法来求解离散方程1517收敛判别条件为各变量误差小于10通过以上方法本文对开放边界抛石路基的温度场和速度场进行了数值计算分析

寒区冻土层水理性质研究

寒区冻土层水理性质研究

寒区冻土层水理性质研究戴长雷;常龙艳;孙思淼;吕雅洁;刘月;黄集华【摘要】寒区是我国重要的地理分区,在我国75%的国土覆盖着多年或季节性冻层土壤.冻层的不透水性、抑制蒸发作用、蓄水调节作用及调节产汇流作用,使寒区的水文循环机理与过程比非寒区复杂的多.从冻土层的水理性质的相关概念及特征参数、冻层水分迁移、融雪水入渗、冻土保墒和地下水溢流冰模型5个方面梳理国内外冻层水理性质的相关研究成果,并对国际领先研究成果一俄科院冻土所论著《寒区冻结层上水》进行内容分析,以期为农业生产、春季融雪径流、越冬期冻土层持水量、越冬期地表水对潜水补给量计算等研究提供有效的理论依据.【期刊名称】《黑龙江大学工程学报》【年(卷),期】2013(004)002【总页数】9页(P1-8,18)【关键词】冻土层;水理性质;冻土保墒;融雪入渗;溢流积冰;寒区【作者】戴长雷;常龙艳;孙思淼;吕雅洁;刘月;黄集华【作者单位】黑龙江大学水利电力学院,哈尔滨150080;黑龙江大学寒区地下水研究所,哈尔滨150080;黑龙江大学水利电力学院,哈尔滨150080;黑龙江大学寒区地下水研究所,哈尔滨150080;伯明翰大学地理地球环境科学学院,英国伯明翰;黑龙江大学水利电力学院,哈尔滨150080;黑龙江大学寒区地下水研究所,哈尔滨150080;中冶焦耐工程技术有限公司,辽宁大连116065;中冶焦耐工程技术有限公司,辽宁大连116065【正文语种】中文【中图分类】P6410 引言冻土层是季节性冻土区土壤中特有的土层,一般位于土壤包气带中。

其水理性质(容水性质、持水性质、给水性质、透水性质)所表现出的不透水作用、蓄水调节作用和抑制蒸发作用,使冻融期降水下渗、土壤含水率的变化状态和降雨、径流、蒸发的“三水”转换关系均具有不同于无冻条件下的动态规律和特点[1-3]。

在寒区,低温与冻融作用使水文循环机理与过程更加复杂,将非寒区理论和方法应用于寒区水文循环存在诸多问题[4-6],其中突出问题之一就是春季融雪水在包气带冻土层的入渗过程与机理[7-8]。

人工多圈管冻结水热耦合数值模拟研究

人工多圈管冻结水热耦合数值模拟研究

第26卷第2期岩石力学与工程学报V ol.26 No.2 2007年2月Chinese Journal of Rock Mechanics and Engineering Feb.,2007人工多圈管冻结水热耦合数值模拟研究汪仁和1,2,3,李栋伟1,2,3(1. 安徽理工大学土木工程系,安徽淮南 232001;2. 现代矿业工程安徽省重点实验室,安徽淮南 232001;3. 矿山建设工程安徽省高校重点实验室;安徽,淮南 232001)摘要:人工多圈管冻结过程中温度场和水分场的计算,是一个多物理场耦合和非线性数学问题,其影响因素多、相互之间关系复杂。

基于相似理论原理和人工多圈管冻结模型试验,提出渗流方程中的导水系数是温度梯度的函数,建立了多圈管正冻土中水热耦合数学模型,并采用有限元方法实现了对多圈管冻结温度场和水分场耦合的数值分析。

数值模拟温度场的发展趋势、水分场迁移特征与模型试验吻合良好,表明用提出的水热耦合数学模型计算多圈管冻结温度场和水分场是可行的,对多圈管冻结壁设计有参考意义。

关键词:岩土工程;模型试验;温度梯度;水热耦合;数值分析中图分类号:TU 47 文献标识码:A 文章编号:1000–6915(2007)02–0355–05 RESEARCH ON HYDRO-THERMAL COUPLING NUMERICAL SIMULATION WITH ARTIFICIAL MULTI-FREEZING-TUBE CYCLESWANG Renhe1,2,3,LI Dongwei1,2,3(1. Department of Civil Engineering,Anhui University of Science and Technology,Huainan,Anhui232001,China;2. The key laboratory of Modern mining Engineering of Anhui Province,Huainan,Anhui232001,China;3. Mine Building Key ofLaboratory of Colleges and Universities in AnHui Province,Huainan,Anhui232001,China)Abstract:During the coupling analysis of moisture field and temperature field in freezing soil by using artificial multi-freezing-tube cycles,a lot of factors should be considered,including multi-fields and nonlinear problems,etc.. Their relations are complicated. The model tests of freezing process with multi-freezing-tube cycles are conducted. The regulation of variation of temperature field and moisture field with freezing time are obtained. It is proved that the temperature is the source of variation of moisture field;and the hydro-thermal coupling dominant equation is proposed. The finite element method is used to carry out hydro-thermal coupling movement with the hydro-thermal coupling parameters being solved simultaneously. It can simulate the characters of temperature forming and development,freezing moisture field transferring,as well as the coupling relation of the two fields. The results of numerical simulations and model tests have good coincidence. The hydro-thermal coupling analysis in this study has important reference to designing and optimizing freezing wall. The results have promotion role for both technological progress and theoretical development of shaft sinking by artificial freezing process with multi-freezing-tube-cycles in thick overburden.Key words:geotechnical engineering;model testing;temperature gradient;hydro-thermal coupling;numerical analysis收稿日期:2006–04–24;修回日期:2006–05–30基金项目:安徽省自然科学基金资助项目(050440501);现代矿山建设工程安徽省高校重点实验室基金项目;安徽理工大学博士、硕士基金项目作者简介:汪仁和(1956–),男,博士,1982年毕业于淮南矿业学院矿井建设专业,现任教授、博士生导师,主要从事岩土力学方面的教学与研究工作。

柴油机用DPF技术及其再生系统的匹配与优化

柴油机用DPF技术及其再生系统的匹配与优化
大气污染防治
综采工作面旋转风幕隔尘实验研究 王鹏飞"刘荣华"汤"梦"张"文"桂"哲(#!' 柴油机用 I5)技术及其再生系统的匹配与优化 胡乃涛"张学敏"林超群"余"鹏(!%$ 电解铝用滤袋缝合工艺对 5B$J' 过滤效率的影响 杨"勇"沈恒根"李纪伟"周刘轲"薛婷婷(!%! 基于 B+KL= 的轻型车颗粒物排放来源和特征分析 郝艳召"邓顺熙"邱兆文"陈瑞华(!&' 溶胶0凝胶法制备颗粒状催化剂低温 =D?脱硝性能 闫东杰"玉"亚"黄学敏"徐"颖(!$( 华北地区温度和风速对冬季灰霾形成的潜在影响 苏"盛"宋明月"吴长江"许"武"刘光远"陈兴宇(!$# 湍流场中颗粒破碎的数值模拟 刘含笑"姚宇平"郦建国(!(8 兰州市冬季细颗粒中微量金属元素及无机可溶性离子来源分析 王丹璐"张有贤"牛武江(!,, 普通醋纤香烟滤嘴中烟碱分布模拟 李珊红"唐"奇"李彩亭"曾光明"王垭曼"丁倩倩(!'' 小波分析在郑州市供暖期 5B$J' 浓度相关性分析中的应用 马跃先"孙明明"张"杰"殷巍杰"武汉清(!9% 沿面放电活化空气用于亚硫酸铵氧化 韩长民"张"轶"商克峰"吴"森"李"杰"吴"彦(!9'
固体废物处置
添加垃圾焚烧飞灰对烧结矿冶金性能的影响 郑慧敏"刘清才"孟"飞"朱博洪"兰苑培"丁"健(!8% 过氧乙酸和亚铁联用调质强化活性污泥过滤脱水性能 王彩霞"张伟军 "王东升"王庆飞(!8' 堆肥污泥施入黄土对植物生长及 DM 吸收的影响

蒙古高原土壤湿度时空变化格局及其对气候变化的响应

蒙古高原土壤湿度时空变化格局及其对气候变化的响应

第33卷,第1期国土资源遥感Vol.33,No.1 2021年3月REMOTE SENSING FOR LAND&RESOURCES M v-,2021dol:10.6046/atzyyy.2020122引用格式:王佳新,萨楚拉,毛克彪,等.蒙古高原土壤湿度时空变化格局及其对气候变化的响应[J].国土资源遥感,2021,33(1):231-239.(Wang J X, Sa C L,Mao KB,ct al.Temporal and spatial variation of soil moisture in the Mongolian Plateau and its re­sponse to climate change[J].Remote Sensing for Land and Resources,2021,33(1):231-239.)蒙古高原土壤湿度时空变化格局及其对气候变化的响应王佳新2,萨楚拉2,毛克彪2,孟凡浩2,罗敏2,王牧兰2(1.内蒙古师范大学地理科学学院,呼和浩特010022;2.中国农业科学院农业资源与农业区划研究所,北京100081;3.内蒙古自治区遥感与地理信息系统重点实验室,呼和浩特010022)摘要:土壤水作为一种重要的水资源,其时空分布和动态变化对植被的分布和生长状况具有显著影响。

蒙古高原是典型的干旱-半干旱气候区,也是亚欧大陆温带草原的主要构成部分,气候变化导致的土壤含水量变化无疑会对草原生态系统健康及稳定产生直接影响;厘清蒙古高原土壤湿度时空特征及其对气候变化的响应有助于为生态保护相关政策的制定提供科学支撑。

本研究基于GLDAS-Noah土壤湿度数据,运用线性回归分析、相关分析和Mann-KendV l M-K)检验等方法,分析1982—2018年不同深度土壤湿度的时空格局、变化趋势及突变性等特征,并与CRU温度、降水数据相结合探讨土壤湿度对气象因子变化的响应。

黄河未来气候情景下冰情特征演变分析

黄河未来气候情景下冰情特征演变分析

2024年3月水 利 学 报SHUILI XUEBAO第55卷 第3期文章编号:0559-9350(2024)03-0355-12收稿日期:2023-06-30;网络首发日期:2024-01-17网络首发地址:https:??kns.cnki.net?kcms?detail?11.1882.TV.20240116.1723.002.html基金项目:第二次青藏高原综合科学考察研究任务(2019QZKK0207-02);国家自然科学基金项目(U2243221,U2243239,51979291,52009144);中国水科院科研专项(HY0145B032021)作者简介:周中元(1999-),硕士生,主要从事河冰水动力学研究。

Email:zhongyuan123ab@163.com通信作者:王涛(1975-),博士,教授级高级工程师,主要从事河冰水动力学及冰情预报研究。

Email:wangtao@iwhr.com黄河未来气候情景下冰情特征演变分析周中元1,王 涛1,孙亚翡1,陈玉状1,路锦枝1,2(1.中国水利水电科学研究院流域水循环与调控国家重点实验室,北京 100038;2.清华大学水利水电工程系,北京 100084)摘要:了解未来气候变化如何影响河流冰情特征对于研究冰凌洪水灾害、水电生产以及大坝管理运行等问题至关重要。

基于黄河流域气象观测数据以及第六次国际耦合模式比较计划(CMIP6)中8种全球气候模式(GCMs)的日均气温数据,评估了各GCMs在分位数增量映射(QDM)偏差校正前后对于黄河流域凌汛期日平均气温的模拟能力,预估了黄河流域未来凌汛期气温变化趋势。

建立了最大冰厚以及封冻历时预测模型,并预估了黄河防凌重点区域黄河宁蒙段未来最大冰厚以及封冻历时的变化趋势。

研究表明,在SSP1-2.6、SSP2-4.5和SSP5-8.5三种气候情景下预计2015—2100年期间黄河流域凌汛期平均气温升温速率分别为0.014、0.031和0.067℃?a,黄河巴彦高勒断面21世纪内最大冰厚值将会分别下降8.5、19.5和39.5cm。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

二、硫酸盐渍土冻结温度
2.1 硫酸盐渍土分布
青藏高原北麓河地区局部盐 渍土,该地区为多年冻土区,冻
土类型为富冰冻土,土质为亚粘
土,存在大面积的盐渍化现象 。
离子分布
土壤中SO42-、Cl-、Na+离子的含量占绝对优势。盐渍土为含 氯化钠的硫酸盐渍土。
二、硫酸盐渍土冻结温度
2.2 基本假定
在特定含水量下(不考虑含水量对冻结温度的影
M M X X 0 H2O M M z M X M z X 0 H2O
盐分的活度可表示为:
a ( M M X
X
0 )(m ) a w, 0
其中γ±表示离子的平均活度系数:
M X 1/ ( M N )
二、硫酸盐渍土冻结温度
1
2 m Na m Cl m SO 2
4
二、硫酸盐渍土冻结温度
2.3 硫酸盐溶液冻结温度
cNa2SO4=0.1958mol/l cNa2SO4=0.3127mol/l
cNa2SO4=0.3915mol/l
Na2SO4 和 NaCl 二元溶液冻 结温度的计算值与实验值 十分接近,说明以上计算
二、硫酸盐渍土冻结温度
2.6 盐晶体析出对冻结温度影响
溶液过饱和比定义为:
c S csat
Smax与T的关系
其中,c为溶液浓度,csat为溶 液饱和浓度。 饱和浓度与温度之间的试验关 系为:
csat 0.3165 1.0773
T
Smax与c0的关系
硫酸钠盐渍土初始浓度越高,溶 液最大过饱和比越小,盐晶体越 发容易析出。
在 定 压下 , 如 水分活 度 变化 daw ,即 aw→aw+daw , 则 T→T+ dT,建立新的平衡:
w d w i d i
则有:
w w i ( ) p , aw dΤ ( ) Τ , p daw ( ) p dΤ Τ a w Τ
二、硫酸盐渍土冻结温度
俄-中输油管线 新藏公路
大连
黑大公路、哈大客专
玉树
西线南水北调
青藏铁路
青藏高速公路
大理
青康公路(高速)
一、研究背景
内陆盐渍土 沿海盐渍土
高原盐渍土
我国盐渍土面积约20多万平方公里,居世界第四。
一、研究背景
冻融作用与盐渍化 作用重叠、相互影 响:
伴随水热过程的 盐分迁移
影响路基工程热、 力、水文状况 路基稳定性问题
响),土体冻结温度由两个部分组成。
T f 土 T f 液 T
溶液冻结温度:取 决于溶液浓度和离 子类型。
土体影响:取决于 孔隙大小和固-液表 面自由能。
二、硫酸盐渍土冻结温度
2.3 硫酸盐溶液冻结温度
在压力p时,溶液的凝固点为Tf,此时液相与固相平衡。
w (Τ , p, aw ) i (Τ , p)
2.3 硫酸盐溶液冻结温度
溶液中水的化学势为:
w w RT ln aw
代入微分方程有:
RT S w dΤ daw S dΤ m ,i aw
纯水和溶液的凝固点分别为Tf*和Tf,对上式积分变
形后得溶液的冰点温度:
Tf 1 T f T f R L wi ln a w T f
过饱和溶液中生长1mol晶体时,其体系的吉布斯自由能降低为:
G RT ln a a0
对于纯组分体有:
( G )T Vm p
则特定温度下可变为:
p RT a ln Vm a0
二、硫酸盐渍土冻结温度
对于晶体若其组成为固定形式MυMXυX·υ0H2O,其由υM个带
zM电荷的正离子和υX个带zX电荷的负离子和υ0个水分子组成。 其晶体溶解反应为:
二、硫酸盐渍土冻结温度
土样初始孔径分布与累计进汞量关系
二、硫酸盐渍土冻结温度
土样降温曲线
硫酸钠含盐量3.8%
硫酸钠含盐量2.6%
汇报提纲
一. 研究背景
二. 硫酸盐渍土冻结温度 三. 晶体生长的微观过程 四. 土体冻结过程水-热-力相互作用机理
三、晶体生长的微观过程
3.1 相变驱动力
熔体-晶体体系共存,两相自由能之差为:
二、硫酸盐渍土冻结温度
2.6 盐晶体析出对冻结温度影响
硫酸钠盐渍土冻结温度
对于含水量为18%
的青藏高原粉质粘土,
当含盐量大于1.8%时, 土体冻结温度以上会有 盐晶体析出,造成冻结 温度先回升后递降。
(Wan & Lai*, Permafrost and Periglac. Process, 2015, IF=3.0).
融化期
一、研究背景
在冻融循环和盐分侵蚀的双重作 用下,道路病害更加严重!
一、研究背景
水分 迁移 盐分 迁移
冻胀
结晶
融化
水分聚集 及消散
土体 结构
盐胀
道路破坏
开展对冻融作用、干湿循环和盐分迁移及其综 合作用对寒区工程稳定性的影响研究是非常必要的。
汇报提纲
一. 研究背景
二. 硫酸盐渍土冻结温度 三. 晶体生长的微观过程 四. 土体冻结过程水-热-力相互作用机理
因此,冰-水相变的驱动力为过冷度ΔT
三、晶体生长的微观过程
3.2 成核临界半径及形成功
对半径为r的球形晶核的均相成核过程,其总的自由能写为:
4 πr 3 Lcl ΔT 2 ΔGc = + 4 π r γsf * 3 Vm T
自由能 (J)
10
-18
晶核形成是系统中体积自由能和表面自由能共同作用的结果,
2 slVsT f solution Lwi r0
Tf soil
第一项由溶液性质决定:
aw
1
Tf soil
当溶液为纯水时,第一项等于0
第二项由结晶半径及冰液表面自由能决定: Tf solution Tf soil r0
当冰液接触面为平面,土体冻结温度等于溶液冻结温度
二、硫酸盐渍土冻结温度
方法是可行的。
二、硫酸盐渍土冻结温度
2.4 土体对冻结温度的影响
冰晶生长过程中
r
受到冰和水压力差的
溶液
冰 晶 生长
pl

作用,△p= pi- pl,结 晶压力与晶体尺寸及 其表面张力有关:
sl
冰 晶 pi
r0
γsl 为冰液表面自由能, r 为孔隙半 径, r0为初始结晶半 径。
p sl
f Ae k
其中: A 1.763ln c0 1 0.645ln c0 0.7787
k 0.0265
T
相同含水量时,溶液初始浓度越大,冻结温度以下未冻含 水量越小,当温度低于-10℃时,未冻含水量基本不再变化。
二、硫酸盐渍土冻结温度
2.8 硫酸盐渍土盐胀力及盐胀对土体影响
强度(kPa)
一、研究背景
沉降变形
水热侵蚀
工程的修建改变了原有冻土环境的水热平衡状态,随着全球气 候变暖和多年冻土的退化,冻融灾害对寒区道路的破坏日益严 重,危及交通安全。
纵向裂缝
道路翻浆
一、研究背景
未冻融(放大100)
冻结期
路基承载力 及结构破坏
9周期(放大100)
反复冻融的盐渍化路 基填料孔隙度增大 路基盐胀及沉陷
η表示温度影响系数:
(0.2343c0 1.0485)e 0.015T
盐结晶引起的未相变含水量 计算值与试验值最大差值不超 过0.4%。
硫酸盐渍土未相变含水量的计算 值与试验值
二、硫酸盐渍土冻结温度
2.7 未相变含水量与温度关系
(2)水结冰引起的含水量变化 冻结温度以下盐晶体析出很少,含水量的变化主要由冰水 相变决定。此时,未冻水量随温度变化规律为:
二、硫酸盐渍土冻结温度
2.7 未相变含水量与温度关系
根据硫酸钠盐渍土信号强度与 温度关系可知,未冻水含量的变
化与盐分结晶和水分结冰有关。
(1)盐晶体析出引起的含水量变化

(180c0 w ) 0 (1.2676r ' 1) w 2.2676 0 ( w r ' 142c0 )
dA dV
二、硫酸盐渍土冻结温度
2.4 土体对冻结温度的影响
冻结过程需要克服固液表面自由能做功,做功的大小等 于降温所需的能量。因此,可建立如下平衡方程:
Liw T f solution
dA T Vs sl dV
Vs表示冰晶的摩尔体积,ΔTγ表示由克服冰液表面自由能引 起的温度降低值。Tf-solution表示溶液的冻结温度。
通过推导可得到盐结晶压力为:
p
RT
Vm
(ln
c 0 ln w ) c0 w, 0
式中:c为溶液的摩尔浓度,c0为溶液的饱和摩尔浓度。水分 的摩尔分数χw。
实验采用青藏高原粉质粘土。含盐量依为1%、2.6%、3.8%, 含水量为18%。干密度均控制在1.82-1.90g/cm3之内,尽量 减小土体空气孔隙率。土将配好的土样装到有机玻璃罐中, 使土样直径为10.0cm,高度为13.4cm。
初始结晶半径:
r0
2 slVsT f solution Lwi T
二、硫酸盐渍土冻结温度
2.5 硫酸盐渍土冻结温度
通过上述分析,硫酸盐渍土冻结温度表示为:
Tf 2 slVsT f solution Tf T R Lwi r0 1 f ln aw Lwi
盐渍土信号强度与温度关系
其中,θ0为初始含水量,c0为溶液初始浓度, r ’为温度T’时
相关文档
最新文档