杜紅普，申芙榮，李 敏

（1.河北工業(yè)大學(xué) 能源與環(huán)境工程學(xué)院，天津 300401；2.河北工業(yè)大學(xué) 土木工程學(xué)院，天津 300401；3.河北省土木工程技術(shù)研究中心，天津 300401）

0 前言

溫度上升，將引起土中固、液相組分的體積膨脹，誘發(fā)孔隙液向混合物邊界的排出趨勢(shì)，改變混合物中的孔隙壓力、應(yīng)力與應(yīng)變場(chǎng).力與熱的相互作用在地下建筑、道路、土壩等巖土工程領(lǐng)域以及儲(chǔ)熱系統(tǒng)中均有體現(xiàn)，尤其是隨著儲(chǔ)熱品質(zhì)的提高，熱作用不僅制約著濕份遷移，而且對(duì)巖土介質(zhì)的強(qiáng)度、變形均有一定的影響.伴隨著熱濕遷移研究成為主流的同時(shí)，針對(duì)收縮、形變引起孔隙結(jié)構(gòu)變化的力學(xué)特性也逐步引起人們的重視.

Reuss[1]研究認(rèn)為，當(dāng)儲(chǔ)熱溫度高于60℃后，埋地?fù)Q熱器附近的土層在克服不同深度地應(yīng)力后引發(fā)局部孔隙結(jié)構(gòu)的改變，進(jìn)而大大削弱其傳熱速率.Sanavia和Hu[2-3]指出水分遷移對(duì)干燥收縮和土壤熱導(dǎo)率的影響顯著.Morit[4]借助室內(nèi)和現(xiàn)場(chǎng)大尺度的軟粘土儲(chǔ)能試驗(yàn)，研究了熱作用（最高排熱溫度70℃）引起的粘土變形模型.Saxe[5]通過(guò)長(zhǎng)期對(duì)間歇性取、排熱監(jiān)測(cè)，認(rèn)為土壤高溫區(qū)會(huì)發(fā)生熱應(yīng)變，導(dǎo)致結(jié)構(gòu)穩(wěn)定性受到影響.作者 [6]在開(kāi)展儲(chǔ)熱試驗(yàn)過(guò)程中也發(fā)現(xiàn)熱源面出現(xiàn)干涸裂隙，需進(jìn)一步考慮儲(chǔ)熱中的力學(xué)問(wèn)題.

然而，目前土層儲(chǔ)熱的熱力分析多傾向定性提出而未能形成系統(tǒng)研究體系，直接針對(duì)土壤儲(chǔ)熱過(guò)程中熱濕作用下變形位移場(chǎng)的研究尚十分有限.具體應(yīng)用工程中，Wang[7]考慮了地源熱泵與太陽(yáng)能—地下高溫儲(chǔ)熱的聯(lián)合應(yīng)用，對(duì)該系統(tǒng)的熱濕參數(shù)進(jìn)行實(shí)時(shí)采集，綜合分析了不同儲(chǔ)熱策略下的運(yùn)行特性.而關(guān)于熱作用下形變研究多集中在農(nóng)業(yè)地表干裂和巖土工程開(kāi)裂問(wèn)題上.Hu[3]采用多尺度多物理場(chǎng)方法，模擬研究了收縮過(guò)程中的孔隙率變化.Sabnis[8]考慮外界環(huán)境中不同濕度條件下的粘土力學(xué)特性，對(duì)收縮應(yīng)變、彈性模量及拉伸強(qiáng)度等參數(shù)進(jìn)行分析.目前針對(duì)儲(chǔ)熱系統(tǒng)中的回填材料熱縮特性，專門(mén)描述較少.圖1給出了常用回填材料在60℃排熱條件下3 d后的干縮現(xiàn)象，可以看出，最左側(cè)的粘土發(fā)生變形量最大，尤其是靠近埋管的區(qū)域，存在不同方向裂縫，相比而言，砂加膨潤(rùn)土及砂出現(xiàn)干裂的程度較小，尤其是砂無(wú)明顯的裂隙.這主要是由于高飽和度下砂呈流態(tài)，使得回填較為密實(shí)，而砂加膨潤(rùn)土中，2種材料的粘聚程度存在差異，膨潤(rùn)土自身的濕變性很強(qiáng)，在一定程度上加劇了混合材料的離散性，這也是在工程實(shí)踐中常遇到的問(wèn)題.

本文針對(duì)上述3種材料，結(jié)合室內(nèi)土體收縮試驗(yàn)及極差分析方法，研究土樣類型、干密度及溫度3因素影響下的形變差異特性.

圖1 不同回填材料失水收縮現(xiàn)象Fig.1 Shrinkageof backfillmaterials

1 試驗(yàn)材料及試驗(yàn)方法

1.1 試驗(yàn)材料

試驗(yàn)土樣分砂、粘土和砂加膨潤(rùn)土，粘土風(fēng)干過(guò)2mm篩，砂子粒徑范圍為0.1～0.5 mm.鈉基膨潤(rùn)土在干燥器中110℃烘干過(guò)2mm篩，3種材料的基本性質(zhì)見(jiàn)表1.

砂加膨潤(rùn)土的導(dǎo)熱性試驗(yàn)表明，膨潤(rùn)土的添加量存在一個(gè)最優(yōu)比范圍：10%～12%，該數(shù)值同時(shí)也在ASHRAE所推薦的范圍之內(nèi)[9].因此，試驗(yàn)過(guò)程中膨潤(rùn)土的添加量均設(shè)置為11%.

1.2 試驗(yàn)方法

鑒于實(shí)際工程利用過(guò)程中，砂加膨潤(rùn)土及砂回填的條件均為高含水率土層.同時(shí)，天津地區(qū)送樣平均質(zhì)量含水率在20%～25%.3種回填材料均按23%質(zhì)量含水率配制，并放入保濕器中24h.采用分層夯實(shí)法制備土樣，3次壓實(shí)，層間打毛.試樣尺寸直徑50 mm，高110mm，干密度為1.45、1.55和1.65 g/cm3.

利用SS-1型收縮儀測(cè)定收縮過(guò)程，烘箱控制恒溫條件（依據(jù)現(xiàn)場(chǎng)的高溫儲(chǔ)熱情況，選取測(cè)試溫度為30、60和80℃），待烘箱溫度穩(wěn)定后，放入已固定有土樣的收縮儀，調(diào)整百分表讀數(shù)，開(kāi)始記錄.記錄時(shí)間開(kāi)始以1h為間隔，隨后逐步延長(zhǎng)，待讀數(shù)穩(wěn)定后停止試驗(yàn).同期測(cè)定土樣的質(zhì)量變化，測(cè)定結(jié)束后，將土樣放在105℃烘箱中烘至恒溫，反算各階段的含水率.

1.3 正交試驗(yàn)設(shè)計(jì)及試驗(yàn)結(jié)果

選取土樣類型、干密度和溫度為因素，每個(gè)因素包含3個(gè)水平（見(jiàn)表2），以含水率、線縮率和體縮率為評(píng)價(jià)指標(biāo)，進(jìn)行正交試驗(yàn)設(shè)計(jì).

表1 3種材料的基本性質(zhì)Tab.1 Characteristic ofmaterials

表2 收縮試驗(yàn)涉及的因素及水平Tab.2 Factorsand levelof the shrinkage test

表3 正交試驗(yàn)設(shè)計(jì)及測(cè)定結(jié)果Tab.3 Orthogonalexperimentaldesign andmeasurement results

正交試驗(yàn)設(shè)計(jì)條件下線索率、體縮率、收縮系數(shù)及失水率結(jié)果見(jiàn)表3.可以發(fā)現(xiàn)：粘土的收縮性較大，遠(yuǎn)大于砂和砂加膨潤(rùn)土的.在不考慮溫度和干密度的影響，粘土的線縮率在4%～6%、砂土的線縮率在0.36%～0.58%、砂加膨潤(rùn)土的在0.21%～0.42%.對(duì)于粘土而言不僅存在縱向收縮，而且存在橫向收縮，粘土的體縮率可達(dá)15.22%.砂土和砂加膨潤(rùn)土在橫向變形幾乎沒(méi)有變形，縱向的變形較小，兩種土樣的體縮率僅為3%～5%.

2 試驗(yàn)結(jié)果分析

2.1 3種土樣的收縮規(guī)律

隨溫度的升高，3種土樣的收縮性均增強(qiáng)（圖2）.粘土的線縮率隨時(shí)間呈先快速增長(zhǎng)隨后區(qū)趨于穩(wěn)定的變化過(guò)程，曲線近似為拋物線型.溫度越高，坡度越陡，初期的變形越劇烈，達(dá)到變形穩(wěn)定的時(shí)間越快.在30℃溫度下，變形穩(wěn)定時(shí)間為300 h；在60℃溫度下，變形穩(wěn)定時(shí)間為160 h；在80℃溫度下，變形穩(wěn)定時(shí)間不到60 h.

砂土的線縮率隨烘干時(shí)間的變化較小，曲線僅在開(kāi)始的十幾個(gè)小時(shí)內(nèi)存在一定的弧度，隨后變化曲線趨于一條直線.砂加膨潤(rùn)土的線縮率隨烘干時(shí)間的變化趨勢(shì)和砂土的相似，整體略緩于砂.

圖2 不同溫度下線縮率隨時(shí)間的變化Fig.2 Changesof linear shrinkageatdifferent temperatures

30℃溫度條件下，粘土、砂和砂加膨潤(rùn)土的含水率隨時(shí)間的變化趨勢(shì)較為相近（圖3）.隨時(shí)間延長(zhǎng)，土樣中的含水率逐漸降低.總體上，砂的變化趨勢(shì)最陡，其次是粘土和砂加膨潤(rùn)土.這主要與土中水的賦存狀態(tài)有關(guān)，砂的孔隙大，賦存與土顆粒孔隙間水多為自由水，易于散失；而對(duì)于粘土和膨潤(rùn)土而言，由于雙電子層的作用，在一定程度上，緩和了水分的散失速度.當(dāng)溫度達(dá)60℃以上，砂的含水率隨時(shí)間的變化曲線與粘土和砂加膨潤(rùn)土不再一致.砂的含水率在初期出現(xiàn)一個(gè)明顯的降低，隨后趨于穩(wěn)定.60℃條件下，砂在150h內(nèi)含水率從23%下降到3%后穩(wěn)定；80℃條件下，砂在30 h內(nèi)含水率從23%下降到2.5%后穩(wěn)定.在這兩個(gè)溫度條件下，粘土和砂加膨潤(rùn)土的含水率隨時(shí)間的變化曲線接近且緩于砂，60℃條件下，粘土和砂加膨潤(rùn)土的含水率穩(wěn)定時(shí)間為400 h左右；80℃條件下，在110 h左右.總體上，溫度越高，曲線變化均越陡，達(dá)到穩(wěn)定的時(shí)間也越快，3種土樣中，溫度對(duì)于砂土的含水率變化的影響最敏感.

圖3 不同溫度下含水率隨時(shí)間的變化Fig.3 Changes of moistureat different temperatures

對(duì)比3種類型土的線縮率隨含水率的變化曲線（圖4），可以發(fā)現(xiàn)，隨土樣含水率的降低，土的線縮率逐步增大.含水率的變化對(duì)于粘土的線縮率影響較大，而對(duì)于砂土和砂加膨潤(rùn)土的影響較小.溫度越高，線縮率與含水率的關(guān)系曲線越陡.

圖4 不同溫度下線縮率隨含水率的變化Fig.4 Relation between linear shrinkage and moistureat different temperatures

2.2 極差分析

依據(jù)式 (1)、式 (2)可得到線縮率和體縮率．表4、表5中的為因素對(duì)線縮率、體縮率、收縮系數(shù)及失水率總影響程度；為i因素下的線縮率、體縮率、收縮系數(shù)及失水率； 為 平均值.

表4 線縮率和體縮率極差分析結(jié)果Tab.4 Range analysisof line and volume shrinkage線縮率

從表4和表5中可以看出，3因素對(duì)土樣收縮性影響的主次關(guān)系為：土的類型（X1）gt;溫度（X2）gt;干密度（X3）.土樣的收縮性主要取決于土的類型.在線縮率、體縮率和收縮系數(shù)指標(biāo)方面的最優(yōu)組合為：砂加膨潤(rùn)土—30℃—1.65 g/cm3（X13X21X31）.即砂加膨潤(rùn)土的抗縮性較小，密度越大，土樣的收縮性越小，溫度越低，土樣的收縮性越小.

表5 收縮系數(shù)和失水率極差分析結(jié)果Tab.5 Range analysisof shrinkage coefficientand water loss rate

3 結(jié)論

1）粘土的收縮形變遠(yuǎn)大于砂和砂加膨潤(rùn)土，砂的形變略大于砂加膨潤(rùn)土.粘土在高溫條件下可達(dá)下的變形量可達(dá)6%，而砂及砂加膨潤(rùn)土的形變僅發(fā)生在開(kāi)始的幾個(gè)小時(shí)內(nèi)，隨后迅速趨于穩(wěn)定.溫度高低對(duì)初始變形影響顯著，溫度越高，土樣達(dá)到變形穩(wěn)定的時(shí)間越短.

2）在3種溫度的影響作用下，粘土、砂及砂加膨潤(rùn)土的含水率隨時(shí)間分布與自身的持水能力相吻合，遞減程度依次是砂加膨潤(rùn)土、粘土和砂.溫度越高，含水率變化趨勢(shì)越陡，材料內(nèi)水份達(dá)到平衡的時(shí)間越快.

3）隨溫度提高，粘土線收縮率與含水率的分布逐步由近似線性趨于指數(shù)，而砂和砂加膨潤(rùn)土由于自身較小的收縮特性，受含水率影響程度較小.隨溫度增加，砂土土樣的含水率變化較大，而砂加膨潤(rùn)土土樣的含水率變化較為緩慢，膨潤(rùn)土的添加有利于緩解砂土的失水速率，改善砂的失水收縮性.

4）3因素對(duì)土樣收縮性影響的主次關(guān)系為：土的類型（X1）gt;溫度（X2）gt;干密度（X3）.土樣的收縮性主要取決于土的類型.綜合線縮率、體縮率和收縮系數(shù)指標(biāo)，最優(yōu)組合為：砂加膨潤(rùn)土—30℃—1.65g/cm3.

[1]ReussM，Beck M，Müller JP.Design of a seasonal thermalenergy storage in theground[J].Solar Energy，1997，59（4-6）：247-257.

[2]Sanavia L，Bonifetto G，LalouiL.Thermo-elasto-plastic consolidation analysisw ithwater phase change[C]//PoromechanicsV@sProceedingsof the Fifth BiotConferenceon Poromechanics.ASCE，2013-2022.

[3]Hu LB，Péron H，HueckelT，etal.Dryingshrinkageofde for mableporous media：mechanis msinduced by the fluid removal[J].ASCEGeotechnical SpecialPublication，2007，157：1-10.

[4]Moritz L.Geotechnicalpropertiesof clay atelevated temperatures[R].Swedish Geotechnical Institute，ReportSGI47，1995.

[5]Saxe SD.Ground sourceheatstorageand thermo-physical responseof softclay[D].USA：Massachusetts Instituteof Technology，2009.

[6]杜紅普.土壤高溫?zé)釢襁w移過(guò)程中的溫/濕度峰值現(xiàn)象研究 [D].天津：河北工業(yè)大學(xué)，2010.

[7]Wang EY，F(xiàn)ung A S，QiCY，LeongW H.Performancepredictionofahybrid solarground-sourceheatpump system[J].Energy and Buildings，2012，47：600-611.

[8]Sabnis A，Abdallah I，Nazarian S，etal.Developmentof shrinkagemodel for claysbased onmoisture variation[C]//Experimental and Applied Modeling of Unsaturated Soils.ASCE，2010：166-172.

[9]Wang Huajun，Lu Junchao，Qi Chengying.Thermal conductivities of sand-bentonitemixtures as a backfillmaterial of geothermal borehole[J].TransactionsGeothermalResourcesCouncil，2011，35：1135-1138.