凍土水熱力特性分析新方法PPT學(xué)習(xí)教案

1、會(huì)計(jì)學(xué)1凍土水熱力特性分析新方法凍土水熱力特性分析新方法格格陵陵蘭蘭加拿加拿大大俄俄羅羅斯斯Svalbard冰島冰島阿拉斯加阿拉斯加挪挪威威蒙古蒙古太平洋北冰洋大西洋一、研究背景一、研究背景北 半 球 多 年北 半 球 多 年凍土凍土占陸地面積占陸地面積的的24%24%。多 年 凍 土多 年 凍 土 和和季節(jié)凍土季節(jié)凍土大約占大約占陸地面積的陸地面積的57%57%。中 國 多 年 凍中 國 多 年 凍土土占國土面積的占國土面積的22.4%22.4%，為，為世界第世界第三三凍土大國凍土大國。中國中國中中國國第1頁/共71頁大理玉樹第2頁/共71頁一、研究背景一、研究背景內(nèi)陸鹽漬土內(nèi)陸鹽漬土第3頁

2、/共71頁一、研究背景一、研究背景凍融作用凍融作用與與鹽漬化鹽漬化作用作用重疊、相互影重疊、相互影響：響：伴隨水熱過程的伴隨水熱過程的鹽分遷移鹽分遷移影響路基工程熱影響路基工程熱、力、水文狀況、力、水文狀況路基穩(wěn)定性問題路基穩(wěn)定性問題第4頁/共71頁溫度變化溫度變化（冰強(qiáng)度、膠結(jié)力）（冰強(qiáng)度、膠結(jié)力）凍土強(qiáng)度變化凍土強(qiáng)度變化 （冰（冰- -水相變）水相變）凍、融土力學(xué)特性激變凍、融土力學(xué)特性激變0100020003000400050006000-10-8-6-4-202溫度（）強(qiáng)度（kPa）水分遷移水分遷移分分凝冰凝冰融化融化工程病害工程病害高含冰量高含冰量上上限及溫度變化限及溫度變化工程病害

3、工程病害量變：升溫量變：升溫承載力或強(qiáng)度降低承載力或強(qiáng)度降低質(zhì)變：融化質(zhì)變：融化承載力或強(qiáng)度喪失承載力或強(qiáng)度喪失一、研究背景一、研究背景第5頁/共71頁縱向裂縫縱向裂縫道路翻漿道路翻漿水熱侵蝕水熱侵蝕一、研究背景一、研究背景第6頁/共71頁未凍融（放大未凍融（放大100）9周期（放大周期（放大100）凍結(jié)凍結(jié)期期融化融化期期反復(fù)凍融的鹽漬化路反復(fù)凍融的鹽漬化路基填料孔隙度增大基填料孔隙度增大路基鹽脹及沉陷路基鹽脹及沉陷路基承載力路基承載力及結(jié)構(gòu)破壞及結(jié)構(gòu)破壞一、研究背景一、研究背景第7頁/共71頁在凍融循環(huán)和鹽分侵蝕的雙重作在凍融循環(huán)和鹽分侵蝕的雙重作用下，道路病害更加嚴(yán)重！用下，道路病害更加

4、嚴(yán)重！一、研究背景一、研究背景第8頁/共71頁道路破壞道路破壞 開展對(duì)凍融作用、干濕循環(huán)和鹽分遷移及其綜開展對(duì)凍融作用、干濕循環(huán)和鹽分遷移及其綜合作用對(duì)寒合作用對(duì)寒區(qū)工程穩(wěn)定性區(qū)工程穩(wěn)定性的影響研究是非常必要的。的影響研究是非常必要的。結(jié)結(jié)晶晶土土體體結(jié)結(jié)構(gòu)構(gòu)鹽脹鹽脹鹽分鹽分遷移遷移一、研究背景一、研究背景第9頁/共71頁匯報(bào)提綱匯報(bào)提綱一一. . 研究背景研究背景二二. . 硫酸鹽漬土凍結(jié)溫度硫酸鹽漬土凍結(jié)溫度三三. . 晶體生長的微觀過程晶體生長的微觀過程四四. . 土體凍結(jié)過程水土體凍結(jié)過程水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理第10頁/共71頁土壤中土壤中SO42-、Cl-、N

5、a+離子的含量占絕對(duì)優(yōu)勢離子的含量占絕對(duì)優(yōu)勢。鹽漬土為。鹽漬土為含氯化鈉的含氯化鈉的硫酸硫酸鹽漬土鹽漬土。青藏高原北麓河青藏高原北麓河地區(qū)地區(qū)局部鹽局部鹽漬土漬土，該地區(qū)為多年凍土區(qū)，凍該地區(qū)為多年凍土區(qū)，凍土類型為富冰凍土，土質(zhì)為亞粘土類型為富冰凍土，土質(zhì)為亞粘土，存在大面積的鹽漬化現(xiàn)象土，存在大面積的鹽漬化現(xiàn)象 。2.1 2.1 硫酸鹽漬土分布硫酸鹽漬土分布二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度離子離子分布分布第11頁/共71頁在特定含水量下（不考慮含水量對(duì)凍結(jié)溫度的影在特定含水量下（不考慮含水量對(duì)凍結(jié)溫度的影響），響），土體土體凍結(jié)溫度由兩個(gè)部分組成凍結(jié)溫度由兩個(gè)部分組成。TTTf

6、f液土溶液凍結(jié)溫度溶液凍結(jié)溫度：?。喝Q于決于溶液濃度溶液濃度和和離離子類型子類型。二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.2 2.2 基本假定基本假定土體影響土體影響：取決于：取決于孔隙大小孔隙大小和和固固- -液液表面自由能表面自由能。第12頁/共71頁),(),(papiwwiiwwddddaadpiwpwwapww)()()(,二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.3 2.3 硫酸鹽溶液凍結(jié)溫硫酸鹽溶液凍結(jié)溫度度在壓力在壓力p時(shí)，溶液的時(shí)，溶液的凝固點(diǎn)為凝固點(diǎn)為Tf，此時(shí)液相與固相平，此時(shí)液相與固相平衡。衡。在定壓下，如水分活度在定壓下，如水分活度變化變化daw，即即

7、awaw+daw，則，則TT+ dT，建立新的平衡：，建立新的平衡：則有：則有：第13頁/共71頁wwwaRT lndSdaaRTdSimwww,二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.3 2.3 硫酸鹽溶液凍結(jié)溫硫酸鹽溶液凍結(jié)溫度度溶液中溶液中水的化學(xué)勢水的化學(xué)勢為：為：代入微分方程有：代入微分方程有：純純水水和溶液和溶液的凝固點(diǎn)的凝固點(diǎn)分別分別為為Tf*和和Tf，對(duì)上式積分變，對(duì)上式積分變形后得形后得溶液的冰點(diǎn)溫度溶液的冰點(diǎn)溫度：fwwifffTaLRTTTln1第14頁/共71頁iiwwmMa1000ln二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2 2.3 .3 硫酸鹽溶液凍結(jié)

8、溫硫酸鹽溶液凍結(jié)溫度度溶液的凍結(jié)溫度取決于水分活度，溶液的凍結(jié)溫度取決于水分活度，采用采用Pitzer離子模離子模型型進(jìn)行計(jì)算。進(jìn)行計(jì)算。Mw為水的摩爾質(zhì)量，為水的摩爾質(zhì)量，m為溶質(zhì)的摩爾質(zhì)量分?jǐn)?shù)，為溶質(zhì)的摩爾質(zhì)量分?jǐn)?shù)，i為離為離子子類型，類型，為滲透系數(shù)。為滲透系數(shù)。2424242442422424221212123SOClNaSOClNaSOClSOClSONaSONaNaClClNaSOClNammmmmZCBmmZCBmmbIIAmmmSONaNaCl,/)()(第15頁/共71頁cNa2SO4=0.1958mol/lcNa2SO4=0.3127mol/lcNa2SO4=0.3915m

9、ol/lNa2SO4和和NaCl二二元元溶液凍溶液凍結(jié)溫度結(jié)溫度的的計(jì)算計(jì)算值與實(shí)驗(yàn)值值與實(shí)驗(yàn)值十分接近十分接近，說明以上，說明以上計(jì)算計(jì)算方法是可行的方法是可行的。二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.3 2.3 硫酸鹽溶液凍結(jié)溫硫酸鹽溶液凍結(jié)溫度度第16頁/共71頁sl為冰為冰液表面自由能液表面自由能，r為孔隙為孔隙半徑，半徑， r0為初始結(jié)晶半為初始結(jié)晶半 徑。徑。 冰 晶 生長pisl pl r r0冰 晶溶液冰晶生長過程中冰晶生長過程中受到受到冰和水壓力差的冰和水壓力差的作用作用，p= pi- pl，結(jié)晶壓力結(jié)晶壓力與晶體尺寸與晶體尺寸及其表面張力有關(guān)：及其表面張力有關(guān)：d

10、VdApsl二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.4 2.4 土體對(duì)凍結(jié)溫度的影響土體對(duì)凍結(jié)溫度的影響第17頁/共71頁Vs表示冰晶的摩爾體積表示冰晶的摩爾體積，T表示由表示由克克服冰液表面自由能服冰液表面自由能引起的溫度降低值。引起的溫度降低值。Tf-solution表示溶液的凍結(jié)溫度表示溶液的凍結(jié)溫度。凍結(jié)過程需要克服固液表面自由能做功，凍結(jié)過程需要克服固液表面自由能做功，做功的大小做功的大小等于等于降溫降溫所需所需的的能能量。量。因此因此，可可建立建立如下如下平衡方程平衡方程：dVdAVTTLssolutionfiwslTLTVrwisolutionfssl20二、硫酸鹽漬土凍

11、結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.4 2.4 土體對(duì)凍結(jié)溫度的影響土體對(duì)凍結(jié)溫度的影響初始結(jié)晶半徑初始結(jié)晶半徑: :第18頁/共71頁第一項(xiàng)第一項(xiàng)由由溶液性質(zhì)溶液性質(zhì)決定決定: :021rLTVTalnLRTTTwisolutionfsfwwiffslsoilfwa1soilf T02rLTVwisolutionfssl第二項(xiàng)第二項(xiàng)由由結(jié)晶半徑結(jié)晶半徑及及冰液表面自由能冰液表面自由能決定決定: :0rsoilf Tsolutionf T二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.5 2.5 硫酸鹽漬土凍結(jié)溫度硫酸鹽漬土凍結(jié)溫度通過上述分析，通過上述分析，硫酸鹽漬土凍結(jié)溫度硫酸鹽漬土凍結(jié)溫度表

12、示為表示為：當(dāng)溶液為純水時(shí)，第一項(xiàng)等于0當(dāng)冰液接觸面為平面，土體凍結(jié)溫度等于溶液凍結(jié)溫度第19頁/共71頁含鹽土的凍結(jié)溫度含鹽土的凍結(jié)溫度比相應(yīng)比相應(yīng)溶液溶液低低，兩者差值隨濃度增，兩者差值隨濃度增大大而而增大增大。鹽晶體析出前鹽晶體析出前凍結(jié)溫度凍結(jié)溫度隨隨濃度濃度呈現(xiàn)很好的呈現(xiàn)很好的線性線性關(guān)系關(guān)系。二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.5 2.5 硫酸鹽漬土凍結(jié)溫度硫酸鹽漬土凍結(jié)溫度鹽晶體析出第20頁/共71頁溫度回升點(diǎn)為溫度回升點(diǎn)為十水硫酸鈉十水硫酸鈉晶體析出的初始溫度晶體析出的初始溫度，可得，可得到不同硫酸鈉含鹽量的晶體初始析出溫度到不同硫酸鈉含鹽量的晶體初始析出溫度。高含

13、鹽量硫酸鈉鹽漬土降溫圖低含鹽量硫酸鈉鹽漬土降溫圖二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.6 2.6 鹽晶體析出對(duì)凍結(jié)溫度影響鹽晶體析出對(duì)凍結(jié)溫度影響第21頁/共71頁溶液過飽和比溶液過飽和比定義為定義為：satccS 其中，其中，c c為為溶液濃度，溶液濃度，c csatsat為為溶溶液飽和濃度。液飽和濃度。Tc0773131650sat.飽和濃度飽和濃度與與溫度溫度之間的試驗(yàn)關(guān)之間的試驗(yàn)關(guān)系為：系為：二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.6 2.6 鹽晶體析出對(duì)凍結(jié)溫度影響鹽晶體析出對(duì)凍結(jié)溫度影響Smax與與c0的關(guān)系的關(guān)系硫酸鈉鹽漬土初始濃度越高，溶硫酸鈉鹽漬土初始濃度越

14、高，溶液最大過飽和比越小，鹽晶體越液最大過飽和比越小，鹽晶體越發(fā)容易析出發(fā)容易析出。Smax與與T的關(guān)系的關(guān)系第22頁/共71頁硫酸鈉鹽漬土凍結(jié)溫度硫酸鈉鹽漬土凍結(jié)溫度對(duì)于對(duì)于含水量為含水量為18%的青藏高原粉質(zhì)的青藏高原粉質(zhì)粘土，粘土，當(dāng)含鹽量大于當(dāng)含鹽量大于1.8%時(shí)時(shí)，土體凍結(jié)溫，土體凍結(jié)溫度以上會(huì)有鹽晶體析度以上會(huì)有鹽晶體析出，造成出，造成凍結(jié)溫度先凍結(jié)溫度先回升后遞降回升后遞降。二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.6 2.6 鹽晶體析出對(duì)凍結(jié)溫度影響鹽晶體析出對(duì)凍結(jié)溫度影響(Wan & Lai*, Permafrost and Periglac. Process,

15、 2015, IF=3.0).第23頁/共71頁（1 1）鹽）鹽晶體析出引起的含水量變化晶體析出引起的含水量變化二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.7 2.7 未相變含水量與溫度關(guān)系未相變含水量與溫度關(guān)系 鹽漬土信號(hào)強(qiáng)度與溫度關(guān)系鹽漬土信號(hào)強(qiáng)度與溫度關(guān)系根據(jù)根據(jù)硫酸鈉鹽漬土信號(hào)強(qiáng)度硫酸鈉鹽漬土信號(hào)強(qiáng)度與溫度關(guān)系與溫度關(guān)系可知，可知，未凍水含量未凍水含量的的變化與變化與鹽分結(jié)晶鹽分結(jié)晶和和水分結(jié)冰水分結(jié)冰有關(guān)有關(guān)。)(.).()(001422676212676118000crrcwww其中，其中，0為初始含水量，為初始含水量，c0為為溶液溶液初始初始濃度濃度， r為溫度為溫度T時(shí)時(shí)溶

16、液的溶解度溶液的溶解度。第24頁/共71頁對(duì)溫度進(jìn)行修正之后，鹽晶體析出引起的未相變含水對(duì)溫度進(jìn)行修正之后，鹽晶體析出引起的未相變含水量變化為：量變化為：表示表示溫度溫度影響系數(shù)影響系數(shù)：)(.).()(001422676212676118000crrcwwwTec015000485123430.).(硫酸鹽漬土未相變含水量的計(jì)算硫酸鹽漬土未相變含水量的計(jì)算值與試驗(yàn)值值與試驗(yàn)值鹽結(jié)晶引起的鹽結(jié)晶引起的未相變含水未相變含水量量計(jì)算值計(jì)算值與與試驗(yàn)值試驗(yàn)值最大差值最大差值不不超過超過0.4%。二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.7 2.7 未相變含水量與溫度關(guān)系未相變含水量與溫度關(guān)系第

17、25頁/共71頁凍結(jié)溫度以下鹽晶體析出凍結(jié)溫度以下鹽晶體析出很少很少，含水量的變化主要由冰，含水量的變化主要由冰水相變決定水相變決定。此時(shí)，。此時(shí)，未凍水量未凍水量隨溫度隨溫度變化規(guī)律變化規(guī)律為：為：相同含水量時(shí)，溶液初始濃度越大，凍結(jié)溫度以下未凍相同含水量時(shí)，溶液初始濃度越大，凍結(jié)溫度以下未凍含水量越小，當(dāng)溫度低于含水量越小，當(dāng)溫度低于-10-10時(shí)，未凍含水量基本不再變化時(shí)，未凍含水量基本不再變化。kAeTf176310cAln.7787064500.ln.c02650.k（2 2）水結(jié)冰）水結(jié)冰引起的含水量變化引起的含水量變化二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度2.7 2.7 未

18、相變含水量與溫度關(guān)系未相變含水量與溫度關(guān)系其中：其中：第26頁/共71頁二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度0aaRTGln對(duì)于純組分體有：mTVpG)(則特定溫度下可變?yōu)椋?aaVRTpmln第27頁/共71頁二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度OHMMOHXMXMXMzXzM2020鹽分的活度可表示為：00,)(wXMamaXM其中表示離子的平均活度系數(shù)：/1）（XMNM第28頁/共71頁二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度通過推導(dǎo)可得到鹽結(jié)晶壓力為：)ln(ln,000wwmccVRTp式中：c為溶液的摩爾濃度，c0為溶液的飽和摩爾濃度。水分的摩爾分?jǐn)?shù)w。實(shí)驗(yàn)采用

19、青藏高原粉質(zhì)粘土。含鹽量依為1%、2.6%、3.8%，含水量為18%。干密度均控制在1.82-1.90g/cm3之內(nèi)，盡量減小土體空氣孔隙率。土將配好的土樣裝到有機(jī)玻璃罐中，使土樣直徑為10.0cm，高度為13.4cm。第29頁/共71頁二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度土樣初始孔徑分布與累計(jì)進(jìn)汞量關(guān)系第30頁/共71頁二、硫酸鹽漬土凍結(jié)溫度二、硫酸鹽漬土凍結(jié)溫度硫酸鈉含鹽量3.8%硫酸鈉含鹽量2.6%土樣降溫曲線第31頁/共71頁匯報(bào)提綱匯報(bào)提綱一一. . 研究背景研究背景二二. . 硫酸鹽漬土凍結(jié)溫度硫酸鹽漬土凍結(jié)溫度三三. . 晶體生長的微觀過程晶體生長的微觀過程四四. . 土

20、體凍結(jié)過程水土體凍結(jié)過程水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理第32頁/共71頁3.1 3.1 相變驅(qū)動(dòng)力相變驅(qū)動(dòng)力 熔體熔體- -晶體體系共存，兩相晶體體系共存，兩相自由能之差自由能之差為為：一般情況下一般情況下相變溫度區(qū)間較小相變溫度區(qū)間較小，因而近似認(rèn)為：，因而近似認(rèn)為：因此，相變時(shí)的自由能變化量為：因此，相變時(shí)的自由能變化量為：因此因此, ,冰冰- -水相變的驅(qū)動(dòng)力為過冷度水相變的驅(qū)動(dòng)力為過冷度T( )( )TSTTHG-=( )()*=THTH( )()*=TSTS()()*cl*clcl*TTLTTLLTSTTHG=-三、晶體生長的微觀過程三、晶體生長的微觀過程第33頁/共

21、71頁 晶核形成是系統(tǒng)中晶核形成是系統(tǒng)中體積自由能體積自由能和和表面自由能表面自由能共同作用的結(jié)共同作用的結(jié)果，對(duì)半徑為果，對(duì)半徑為r的球形晶核的均相成核過程，其總的自由能寫為的球形晶核的均相成核過程，其總的自由能寫為： Vm為物質(zhì)的摩爾體積；為物質(zhì)的摩爾體積； sf為晶為晶相與液相界面比表面自由能。相與液相界面比表面自由能。成核功與晶體尺寸的關(guān)系臨界晶核半徑為：臨界晶核半徑為：形成臨界球形晶核所需的功為形成臨界球形晶核所需的功為：sf*mclcrTTVLrG234+34= -TLTVrcl*msfcrit2=22223316=TLTGLVclcl*msfcr三、晶體生長的微觀過程三、晶體生長

22、的微觀過程3.2 3.2 成核臨界半徑及形成功成核臨界半徑及形成功第34頁/共71頁單位時(shí)間單位體積內(nèi)單位時(shí)間單位體積內(nèi)超臨界晶胞形成的概率超臨界晶胞形成的概率為為：分子碰撞頻率與活化能有關(guān)：分子碰撞頻率與活化能有關(guān)：超臨界分子的數(shù)目為超臨界分子的數(shù)目為：因此，因此，均勻成核的結(jié)晶率為均勻成核的結(jié)晶率為：crsnndtdNI=RTGexp00-kTGexpncrcrn0kTGexpRTGexpnndtdNIcrs000三、晶體生長的微觀過程三、晶體生長的微觀過程3.3 3.3 結(jié)晶速率結(jié)晶速率第35頁/共71頁結(jié)晶速率與溫度關(guān)系土體凍結(jié)典型溫度時(shí)間曲線結(jié)晶速率方程揭示了兩個(gè)相互競爭的因素：結(jié)晶

23、速率方程揭示了兩個(gè)相互競爭的因素：成核作用成核作用和和晶體生長晶體生長。恒 定跳 躍過 冷遞 降自 由 水凍 結(jié)束 縛 水凍 結(jié)Te大 部 分 束 縛水 凍 結(jié) ，土 體 冷 卻TscTf時(shí) 間0溫 度TTsc相 對(duì) 于 熱 力 學(xué) 平 衡 時(shí) 的 亞 穩(wěn) 態(tài) 過 度 溫 度Tf自 由 水 的 起 始 凍 結(jié) 溫 度T相 對(duì) 于0 的 冰 點(diǎn) 溫 度Te束 縛 水 幾 乎 全 部 凍 結(jié) 時(shí) 的 溫 度臨 界晶 核成 核吸 熱生 長放 熱吸 放 熱平 衡放 熱較 少冷 卻三、晶體生長的微觀過程三、晶體生長的微觀過程3.3 3.3 結(jié)晶速率結(jié)晶速率溫度較低時(shí)，分子動(dòng)能較低，分子難以克服勢壘進(jìn)入晶

24、核，因此在溫度很低情況下，晶體生長反而變得困難。第36頁/共71頁*TTT-*0=fvpsf*SrcosT2= -（1）當(dāng)當(dāng)90o時(shí)時(shí)，晶體不能在，晶體不能在T0*時(shí)進(jìn)入孔隙，并且孔隙越小時(shí)進(jìn)入孔隙，并且孔隙越小，晶體進(jìn)入孔隙的溫度越低；相反，晶體進(jìn)入孔隙的溫度越低；相反， 90o時(shí)時(shí)冰點(diǎn)溫度不會(huì)降冰點(diǎn)溫度不會(huì)降低，冰晶能夠自由的進(jìn)入孔隙。低，冰晶能夠自由的進(jìn)入孔隙。（2）當(dāng)孔隙半徑當(dāng)孔隙半徑rp10-6 m時(shí)時(shí)，孔隙結(jié)構(gòu)對(duì)相變溫度的影響較小。，孔隙結(jié)構(gòu)對(duì)相變溫度的影響較小。三、晶體生長的微觀過程三、晶體生長的微觀過程3.4 3.4 孔隙結(jié)構(gòu)對(duì)凍結(jié)溫度的影響孔隙結(jié)構(gòu)對(duì)凍結(jié)溫度的影響第37頁/共

25、71頁冰冰- -水處于平衡狀態(tài)時(shí)，水處于平衡狀態(tài)時(shí)，水活度與溫度的關(guān)系水活度與溫度的關(guān)系表示為：表示為：aw表示溫度為表示溫度為T時(shí)溶液中水的時(shí)溶液中水的活度，活度， 和和 分別表示溫度分別表示溫度純水在正常熔點(diǎn)下的潛熱和冰水純水在正常熔點(diǎn)下的潛熱和冰水比熱容之差比熱容之差水分活度與溫度關(guān)系當(dāng)當(dāng)T=T*時(shí)時(shí),aw=1,冰和水處于平衡冰和水處于平衡狀態(tài)；當(dāng)狀態(tài)；當(dāng)TT*時(shí)時(shí),aw1,冰逐漸融化成冰逐漸融化成水；當(dāng)水；當(dāng)TT*時(shí)時(shí),aw1,水逐漸凍結(jié)成冰水逐漸凍結(jié)成冰。 TTlnRCTTRTCTHaln*p*p*mw111 *mTH1pC三、晶體生長的微觀過程三、晶體生長的微觀過程3.5 3.5

26、溫度對(duì)水活度的影響溫度對(duì)水活度的影響第38頁/共71頁水凍結(jié)與融化過程的水凍結(jié)與融化過程的相變速率相變速率寫為寫為：上式中，第一項(xiàng)代表冰晶生長，第二項(xiàng)代表冰晶融化。上式中，第一項(xiàng)代表冰晶生長，第二項(xiàng)代表冰晶融化。和和 分別表示正向和反向速率常數(shù)。分別表示正向和反向速率常數(shù)。和和 分別表示水和冰的濃度分別表示水和冰的濃度(kg/m3)；n為材料參數(shù)。為材料參數(shù)。 sinwrlwnwfCaTICaTImiwwi11 fTI rTIlwCsiC三、晶體生長的微觀過程三、晶體生長的微觀過程3.6 3.6 晶體生長的宏觀動(dòng)力學(xué)模型晶體生長的宏觀動(dòng)力學(xué)模型第39頁/共71頁 物質(zhì)總是由化學(xué)勢高的位置向化學(xué)

27、勢低的位置遷移，即物質(zhì)總是由化學(xué)勢高的位置向化學(xué)勢低的位置遷移，即化化學(xué)勢梯度是物質(zhì)遷移的驅(qū)動(dòng)力學(xué)勢梯度是物質(zhì)遷移的驅(qū)動(dòng)力：fm為單位體積流體上所受的體積驅(qū)動(dòng)力。為單位體積流體上所受的體積驅(qū)動(dòng)力。單位體積流體的平均遷移速率為：單位體積流體的平均遷移速率為： 當(dāng)采用當(dāng)采用達(dá)西定律達(dá)西定律來描述水的運(yùn)動(dòng)過程時(shí)，來描述水的運(yùn)動(dòng)過程時(shí)，滲透壓力梯度是滲透壓力梯度是水分遷移的驅(qū)動(dòng)力：水分遷移的驅(qū)動(dòng)力：xVflwmm1xpflwmgxVVlwlwmlwintrlwkk1三、晶體生長的微觀過程三、晶體生長的微觀過程3.7 3.7 水分遷移水分遷移第40頁/共71頁不同尺度多孔介質(zhì)特性冰棱鏡液態(tài)水液態(tài)水凍結(jié)緣

28、基質(zhì)孔隙冰多孔介質(zhì)中冰結(jié)晶與水分遷移三、晶體生長的微觀過程三、晶體生長的微觀過程3.8 3.8 多孔介質(zhì)熱質(zhì)傳遞與結(jié)晶過程的耦合作用多孔介質(zhì)熱質(zhì)傳遞與結(jié)晶過程的耦合作用(Wu & Lai*, Inter. J. Heat and Mass Transfer, 2015, IF=2.5)第41頁/共71頁匯報(bào)提綱匯報(bào)提綱一一. . 研究背景研究背景二二. . 硫酸鹽漬土凍結(jié)溫度硫酸鹽漬土凍結(jié)溫度三三. . 晶體生長的微觀過程晶體生長的微觀過程四四. . 土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理第42頁/共71頁凍土的成分凍土的成分溫度溫度水分水分應(yīng)力應(yīng)

29、力對(duì)流對(duì)流/ /相變相變相平衡相平衡水分遷移水分遷移導(dǎo)熱性能導(dǎo)熱性能相平衡相平衡/ /凍脹力凍脹力孔壓孔壓相平衡相平衡水分遷移水分遷移固體骨架、冰、固體骨架、冰、水、氣水、氣四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.1 4.1 凍土中的物質(zhì)成分及耦合關(guān)系凍土中的物質(zhì)成分及耦合關(guān)系第43頁/共71頁 開展凍融作用下土體水分遷移機(jī)理試驗(yàn)，分析不同工開展凍融作用下土體水分遷移機(jī)理試驗(yàn)，分析不同工況條件下水分遷移過程中的凍脹機(jī)理。況條件下水分遷移過程中的凍脹機(jī)理。溫度溫度水分水分鹽分鹽分變形變形溫度溫度水分水分變形變形土水勢土水勢控溫、補(bǔ)水控溫、補(bǔ)水控溫、

30、加載控溫、加載四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4 4.2 .2 試驗(yàn)研究試驗(yàn)研究第44頁/共71頁（1 1）荷載）荷載的影響的影響受土樣頂部受土樣頂部荷載的影響荷載的影響，土體，土體凍脹凍脹和和凍結(jié)鋒面凍結(jié)鋒面隨時(shí)間的變化過程不同。條件隨時(shí)間的變化過程不同。條件中土體壓縮階段較長，凍結(jié)中土體壓縮階段較長，凍結(jié)90h，條件，條件和條件和條件的凍脹量為的凍脹量為9.18mm和和4.99mm，凍結(jié)鋒面距暖端距離分別為，凍結(jié)鋒面距暖端距離分別為48.9mm和和54.23mm。圖圖4.1 凍脹凍脹-時(shí)間時(shí)間0102030405060708090100-

31、1012345678910BDCA 凍脹 (工況) 凍脹 (工況) 時(shí)間 /h凍脹 /mm01020304050607080901000102030405060708090100DBC 凍結(jié)鋒面(工況) 凍結(jié)鋒面 (工況) 時(shí)間 /h凍結(jié)鋒面 /mm圖圖4.2 凍結(jié)鋒面凍結(jié)鋒面-時(shí)間時(shí)間0.31/cm 50kPa0.31/cm 100kPa0.31/cm 50kPa0.31/cm 100kPa四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理第45頁/共71頁0102030405060708090100-2-101234SSSCFC降溫速率 /h 20mm 4

32、0mm 60mm 80mm 降溫速率, 80mm溫度 /時(shí)間 /h-3-2-101 圖圖5.3 溫度溫度-時(shí)間時(shí)間0102030405060708090100-2-101234SSSCFC降溫速率 /h 20mm 40mm 60mm 80mm 降溫速率, 80mm溫度 /時(shí)間 /h-3-2-101 0.31/cm 50kPa0.31/cm 100kPa圖5.4 含水量-高度高度 /mm含水量 /%溫度變化均呈現(xiàn)溫度變化均呈現(xiàn)3個(gè)階段，不同的荷個(gè)階段，不同的荷載作用下，降溫速率不同；載作用下，降溫速率不同；未凍區(qū)土體壓密，未凍區(qū)土體壓密，孔隙改變，同時(shí)凍孔隙改變，同時(shí)凍結(jié)鋒面以下吸力減小，未凍區(qū)

33、水分遷結(jié)鋒面以下吸力減小，未凍區(qū)水分遷移量大于補(bǔ)給量，導(dǎo)致未凍區(qū)出現(xiàn)非移量大于補(bǔ)給量，導(dǎo)致未凍區(qū)出現(xiàn)非飽和現(xiàn)象。飽和現(xiàn)象。條件條件四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理第46頁/共71頁（2 2） 溫度梯度的影響溫度梯度的影響0102030405060708090100-1012345678910B4B3DCAB1 凍脹 (工況)凍脹 (工況) 時(shí)間 /h凍脹 /mm01020304050607080901000102030405060708090100DCB1 凍結(jié)鋒面 (工況) 凍結(jié)鋒面 (工況) 時(shí)間 /h凍結(jié)鋒面 /mm圖圖4.5 凍脹凍脹

34、-時(shí)間時(shí)間圖圖4.6 凍結(jié)鋒面凍結(jié)鋒面-時(shí)間時(shí)間較大的較大的溫度梯度決定了土體凍脹速率及鋒面移動(dòng)速率溫度梯度決定了土體凍脹速率及鋒面移動(dòng)速率，凍結(jié)，凍結(jié)90 h，條件，條件和條件和條件的凍脹量為的凍脹量為4.99 mm和和9.82 mm，凍結(jié)鋒面距暖端距離分別為，凍結(jié)鋒面距暖端距離分別為54.23 mm和和35.65 mm。0.65/cm 100kPa0.31/cm 100kPa0.31/cm 100kPa0.65/cm 100kPa四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理第47頁/共71頁圖圖4.7 凍脹速率、溫度梯度凍脹速率、溫度梯度-時(shí)間時(shí)間0

35、.31/cm 50kPa0.31/cm 100kPa圖圖4.9 含水量含水量-高度高度高度 /mm含水量 /%受荷載和溫度梯度綜合影響，初始階受荷載和溫度梯度綜合影響，初始階段凍脹速率和溫度梯度的變化呈相反段凍脹速率和溫度梯度的變化呈相反趨勢，之后溫度梯度與凍脹速率呈正趨勢，之后溫度梯度與凍脹速率呈正相關(guān)趨勢，相關(guān)趨勢，溫度梯度是凍脹主導(dǎo)因素溫度梯度是凍脹主導(dǎo)因素。01020304050607080901000.00.81.01.2B4B3B2B1dz/dt /mm/hgrad T=(T0-Tf)/Hf /cm 時(shí)間 /h凍脹速率和溫度梯度 01020304050607080

36、90100-5-4-3-2-101234SSSCFC溫度變化速率 /h 20mm 40mm 60mm 80mm 溫度改變速率, 80mm溫度 /時(shí)間 /h -3-2-101 0.65/cm 100kPa0.65/cm 100kPa溫度梯度凍脹速率圖圖4.8 溫度溫度-時(shí)間時(shí)間條件條件四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理第48頁/共71頁（3 3） 溫度邊界的影響溫度邊界的影響020406080100 120 140 160 180 200 220 240 260 280 30002468101214161820BCA 凍脹 (工況) 凍脹 (工況

37、) 時(shí)間 /h凍脹 /mm020406080100 120 140 160 180 200 220 240 260 280 300020406080100120140160180200 凍結(jié)鋒面 (工況) 凍結(jié)鋒面 (工況) 時(shí)間 /h凍結(jié)鋒面 /mm0.30/cm(+1.0 -5.0)0.30/cm(+1.8 -4.2)0.30/cm(+1.8 -4.2)0.30/cm(+1.0 -5.0)圖圖4.10 凍脹凍脹-時(shí)間時(shí)間圖圖4.11 凍結(jié)鋒面凍結(jié)鋒面-時(shí)間時(shí)間較低的冷端溫度加速了孔隙冰水相變，改變了孔隙冷生構(gòu)造，較低的冷端溫度加速了孔隙冰水相變，改變了孔隙冷生構(gòu)造，凍結(jié)凍結(jié)290h，條件，

38、條件和條件和條件的凍脹量為的凍脹量為9.8mm和和17.6mm，凍結(jié)鋒面距暖端，凍結(jié)鋒面距暖端距離分別為距離分別為50.7mm和和28.1mm。四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理第49頁/共71頁0.31/cm 50kPa圖圖5.13 含水量含水量-高度高度高度 /mm含水量 /%0.65/cm 100kPa凍脹速率圖圖4.12 溫度溫度-時(shí)間時(shí)間條件條件020406080100 120 140 160 180 200 220 240 260 280 300-5-4-3-2-101234FCSSSC溫度變化速率 /h 35mm 75mm 125

39、mm 145mm 溫度變化速率, 145mm溫度 /時(shí)間 /h -3-2-101 020406080100 120 140 160 180 200 220 240 260 280 300-5-4-3-2-101234FCSSSC溫度改變速率 /h 50mm 100mm 140mm 180mm 溫度變化速率, 140mm溫度 /時(shí)間 /h -3-2-101 溫度變化均呈現(xiàn)溫度變化均呈現(xiàn)3個(gè)階段，個(gè)階段，不同的溫不同的溫度邊界下，土體溫度變化速率不同；度邊界下，土體溫度變化速率不同；受冷端溫度邊界的影響，兩種條件下受冷端溫度邊界的影響，兩種條件下的冰透鏡體分布不同，的冰透鏡體分布不同，較低的冷端溫

40、較低的冷端溫度減少了冰透鏡體數(shù)量和厚度。度減少了冰透鏡體數(shù)量和厚度。0.30/cm(+1.8 -4.2)0.30/cm(+1.0 -5.0)四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理第50頁/共71頁基于水分基于水分遷移及凍脹試驗(yàn)，從遷移及凍脹試驗(yàn)，從多孔介質(zhì)理論多孔介質(zhì)理論出發(fā)，考慮土體單元出發(fā)，考慮土體單元質(zhì)量守恒質(zhì)量守恒、能量守恒能量守恒及及靜力平衡靜力平衡，構(gòu)建水，構(gòu)建水- -熱熱- -力耦合理論模型力耦合理論模型。土顆?？紫端哥R冰冷端暖端Tc ZcTf ZfTl ZlTm ZmTw Zw圖圖4.15 單向凍結(jié)原理圖示單向凍結(jié)原理圖示圖圖4.

41、14 正凍土體示意正凍土體示意四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究第51頁/共71頁 （1 1）土）土體特性體特性含冰量：含冰量： 導(dǎo)導(dǎo)水水率：率： 熱物理熱物理參數(shù)：參數(shù)： 901nSkkiiiwisCnSCSnCnC11iinSiSnwns11四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究00TT T T 0110TTSi第52頁/共71頁mTbaEs（2 2）壓縮）壓縮模量模量（3 3）土）土體容重體容重gSnnSngwiiis11gnnw

42、s0001 凍土強(qiáng)度對(duì)溫度有明顯的依賴性，凍土壓縮模量與凍土強(qiáng)度對(duì)溫度有明顯的依賴性，凍土壓縮模量與溫度之間存在以下經(jīng)驗(yàn)式溫度之間存在以下經(jīng)驗(yàn)式四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究第53頁/共71頁單位單位土體組成及各組分的體積含量土體組成及各組分的體積含量 （4 4）質(zhì)量守恒）質(zhì)量守恒方程方程componentsicewatersoil grains1nSin(1-Si)1-n四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究假設(shè)假設(shè)t和和t+t時(shí)

43、刻時(shí)刻各組分的各組分的質(zhì)量分別為質(zhì)量分別為 和和 ，忽略土顆粒的流入或移出，對(duì)于體積為，忽略土顆粒的流入或移出，對(duì)于體積為dV的微元體，冰、水體系質(zhì)量守恒：的微元體，冰、水體系質(zhì)量守恒：tttswimmm，ttttttswimmm，t)(t)(ttdAdzzvvdAdzzvvmmdAvdAvmmwwwiiiwiwwiiwittttttdAdzzvzvtmtmwwiiwi)(dVnSmiiidVSnmiww)1 ( dVnmss)1 ( 第54頁/共71頁 四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究dAdzzvzvtdVS

44、ntdVnSwwiiiwii)()1 ()(dVndVs 1)(1)1 ()(wwiiiwiiiwivvtnnSStTTSn第55頁/共71頁 （5 5）水分）水分遷移驅(qū)動(dòng)力遷移驅(qū)動(dòng)力四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究kvwgmsTpgPwwpzgTSkTT凍結(jié)區(qū)溫度凍結(jié)區(qū)溫度T和滲透壓力和滲透壓力Pw不是獨(dú)立的，則不是獨(dú)立的，則kT =0第56頁/共71頁 （5 5）水分）水分遷移驅(qū)動(dòng)力遷移驅(qū)動(dòng)力四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究

45、)ln(-0TTLPPiiwwawwPPPaiiPPP0TTTiwporP)(PP15 . 1)1 (iSwiporwwiaiwwPTTTLPP)1 ()(ln()1 ()(1 (00第57頁/共71頁 （6）能量守恒）能量守恒方程方程 水分遷移帶走的熱量水分遷移帶走的熱量 四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究t)dVnS(LiiTdVvCwwdV)PDT(wp-第58頁/共71頁 土體變形土體變形導(dǎo)致孔隙結(jié)構(gòu)變化，從而影響土體的熱物理性質(zhì)，導(dǎo)致孔隙結(jié)構(gòu)變化，從而影響土體的熱物理性質(zhì)，單位時(shí)間土體變形引起的熱量變化

46、，即單位時(shí)間土體變形引起的熱量變化，即變形與熱運(yùn)移耦合項(xiàng)變形與熱運(yùn)移耦合項(xiàng) dVtTnv1四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究tTnTvCPDTtnEtTEvwwwpnT)1 ()-(TSnLCEiiTnSLEiin1第59頁/共71頁 （7 7）靜）靜力平衡方程力平衡方程 由于由于固體骨架本身的壓縮性遠(yuǎn)小于孔隙的壓縮性固體骨架本身的壓縮性遠(yuǎn)小于孔隙的壓縮性 ，假設(shè)假設(shè)固體固體顆粒不可壓顆粒不可壓，土體變形由于孔隙的改變而產(chǎn)生土體變形由于孔隙的改變而產(chǎn)生。 0-dzd四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱

47、熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究porPzuv應(yīng)變應(yīng)變與與位移位移之間的關(guān)系表示為：之間的關(guān)系表示為：Tev假定凍土為假定凍土為彈性材料彈性材料，總應(yīng)變包括，總應(yīng)變包括彈性應(yīng)變彈性應(yīng)變和和溫度應(yīng)變溫度應(yīng)變：esdEd-彈性彈性應(yīng)力應(yīng)力- -應(yīng)變應(yīng)變關(guān)系為關(guān)系為：第60頁/共71頁 nnne1-00TTT靜力平衡靜力平衡方程：方程：四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究不考慮土體流失不考慮土體流失，土體顆?？傮w積為，土體顆粒總體積為常數(shù)。常數(shù)。常數(shù)Vn1nnVVndnVdV001于是

48、可得彈性應(yīng)變和溫度應(yīng)變的表達(dá)式分別為于是可得彈性應(yīng)變和溫度應(yīng)變的表達(dá)式分別為：01-00zPTTnnnEzpors）（第61頁/共71頁 用用分凝溫度分凝溫度作為冰透鏡體萌生和結(jié)束的控制溫度（作為冰透鏡體萌生和結(jié)束的控制溫度（Konrad and Morgenstern, 1980) 用用最大冰壓力最大冰壓力或中斷應(yīng)力來控制（或中斷應(yīng)力來控制（Miller, 1978; Gilpin, 1980; J.F. Nixon, 1991; Konrad and Duquennoi, 1993; Akagawa et al., 2007; H.R. THOMAS, et al., 2009;） 采用采

49、用土體裂縫形成土體裂縫形成作為冰透鏡體形成的判據(jù)（作為冰透鏡體形成的判據(jù)（Azmatch et al., 2011, 2012） （8 8）凍結(jié)）凍結(jié)緣及冰透鏡體緣及冰透鏡體四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究第62頁/共71頁 Lukas U.Arenson et al.2007Tezaera F.Azmatch et al.2012Lukas U.Arenson et al.2008四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究第63頁/共7

50、1頁cTTnnnEdzPshzpor)(-100初始條件：c=0 分凝壓力：分凝壓力： 由靜力平衡由靜力平衡方程得：方程得：tzhhlsepdzPPsepnn 分凝條件分凝條件與與初始初始孔隙率孔隙率、土體壓縮特性土體壓縮特性、溫度梯度溫度梯度及及應(yīng)力應(yīng)力有有關(guān)。關(guān)。 四、四、土體凍結(jié)過程中水土體凍結(jié)過程中水- -熱熱- -力相互作用機(jī)理力相互作用機(jī)理4.3 4.3 理論研究理論研究分凝條件：分凝條件：sepporPPtlstlssepPTTEPTTnEn)(1)(000hzseptlsepshzpornndzPnnTTnnnEdzP )(-100第64頁/共71頁 從物理現(xiàn)象分析，冰透鏡體的生成還需從物理現(xiàn)象分析，冰透鏡體的生成還需滿足：滿足：冰透鏡體冰透鏡體形成后形成后，孔隙，孔隙通道受阻通道受阻，滲透系數(shù)發(fā)生變化，根據(jù)，滲透系數(shù)發(fā)生變化，根據(jù)上述分析對(duì)滲透