何朋立 郭 力 王在泉

(1.洛陽理工學院土木工程系，河南 洛陽 471023； 2.青島理工大學理學院，山東 青島 266033)

深厚表土層井壁溫度豎向附加應力分析

何朋立1郭 力1王在泉2

(1.洛陽理工學院土木工程系，河南 洛陽 471023； 2.青島理工大學理學院，山東 青島 266033)

深厚表土層豎井井壁破裂是嚴重危害煤礦生產(chǎn)的地質災害。溫度作用是導致井壁內豎向應力增大而誘發(fā)井壁破裂的一種重要因素。為了分析深厚表土層立井井壁因溫度作用產(chǎn)生縱向膨脹對井壁受力狀態(tài)的影響，在分析井壁與土體相互作用的基礎上，根據(jù)熱力學和彈性理論建立了井壁和土體剪切作用的彈塑性模型，推導了井壁溫度豎向附加應力的計算公式。進一步分析表明：井壁內溫度豎向附加應力隨深度呈不斷增大的趨勢，并在基巖附近達到最大值；溫度豎向附加應力隨著彈塑性剪切段分界點的深度近似呈線性增加。這些結論有助于加深對深厚表土層井壁破裂機理的認識，并對井壁破裂預防治理具有理論指導作用。

縱向膨脹 溫度作用 土體與井壁相互作用 彈塑性剪切模型 溫度豎向附加應力

多年來，我國華東礦區(qū)的井壁破裂事故給煤炭生產(chǎn)造成重大損失[1]。學界經(jīng)過多年持續(xù)不斷的理論分析、實驗模擬和現(xiàn)場實測等大量研究，認為底部含水層疏水給井壁造成的豎向附加力是井壁破裂的主要原因[2-3]。這一重要結論為井壁設計以及破裂預防和治理起到了重要作用。但深厚表土層井壁受到土體、地下水和溫度等多種作用，除了底部含水層疏水形成的豎向附加力外，溫度對井壁影響也不可忽視。文獻[4]通過彈性理論和熱力學分析，得到了井壁中溫度應力的解析解；文獻[5]對井壁溫度應力特征進行了分析，研究了井壁彈性模量、泊松比、井壁溫差等因素與溫度應力的影響；文獻[6-7]在考慮溫度應力影響的基礎上，探討了各種荷載對井壁破裂的影響。但這些關于井壁溫度應力的分析都是把井筒看作兩端自由的厚壁圓筒，按照平面應變模型進行彈性力學求解。事實上，溫度升高還會導致深厚表土層井壁因縱向膨脹在井壁內產(chǎn)生溫度豎向附加應力，與其他荷載形成的豎向應力分量疊加后超過井壁混凝土抗壓強度而造成井壁破裂。

本研究根據(jù)深厚表土層井壁實際所處的工作環(huán)境，考慮在溫度升高時井壁與周圍土體的相互剪切作用，基于彈性力學和熱力學原理，建立了井壁與土體相互作用的彈塑性模型，并推導了溫度豎向附加力應力公式，在此基礎上分析了溫度豎向附加應力對于井壁影響與特征。

1 溫度升高時井壁與土體相互作用模型

溫度對深厚表土層地區(qū)井壁的影響可分為3個方面：第1方面，溫差導致在井壁內產(chǎn)生溫度自應力；第2方面，井筒徑向膨脹變形受到周圍土體的阻止而產(chǎn)生溫度應力；第3方面，井筒縱向膨脹造成井壁與土體發(fā)生相互剪切作用受到土體阻止而產(chǎn)生的溫度豎向附加應力?，F(xiàn)場監(jiān)測表明[8-9]：井口標高和井壁應變均隨季節(jié)呈周期性波動，這是井筒縱向變形的反映。

在高溫季節(jié)，井筒產(chǎn)生縱向膨脹，井壁受到土體阻止其縱向變形的溫度豎向附加力。由于表土層井壁底部固接，井壁縱向伸長向上發(fā)展，井壁和土體間相對位移變形由基巖向地面呈逐漸增大的線性分布(見圖1)。當井壁與周圍土體的剪切相對位移值小于剪切強度對應的極限位移值u0時，井壁與周圍土體之間產(chǎn)生彈性剪切力；當井壁與周圍土體間的相對位移超過剪切強度的位移值時，井壁與土體之間產(chǎn)生的豎向剪切力為塑性剪切力；在彈性剪切段和塑性剪切段的交界面處，彈(塑)性剪切力達到最大值，而交界面的位置則隨溫度變化而變動。

圖1 溫度升高時井壁與土體剪切作用帶示意

2 井壁溫度豎向附加應力計算

為簡化起見，可近似認為彈性段和塑性段溫度升高引起的井壁豎向附加力沿井筒軸向成直線分布(如圖1所示)。塑性段由溫度引起的豎向附加力fOA為

(1)

從圖1可以看出，只要確定了彈性剪切段和塑性剪切段的交界面位置，溫度豎向附加力即可得到求解。鑒于彈性剪切段和塑性剪切段的分界標準時彈性極限剪切位移值u0，因此可以利用井筒OB段由于溫度升高產(chǎn)生的自由伸長量ΔLOB和井筒由于溫度豎向附加力作用的縮短量ΔLt之差等于u0的關系來確定，即

(2)

根據(jù)相似三角形原理，由圖1可以得到

(3)

要計算井筒OB段的伸長量ΔLOB就需要計算該段井筒的軸力Nz，根據(jù)靜力平衡條件，可以得到井筒OB段軸力Nz1為

(4)

根據(jù)胡克定律，可以得到ΔLt為

(5)

由熱力學原理，對于OB段井筒由于溫度升高ΔT產(chǎn)生的伸長量可以通過式(6)求得:

(6)

立井壁溫度升高和內外壁溫差引起的井壁伸長量和根據(jù)溫度應力場理論計算。由于井壁內溫度沿縱向梯度變化很小，可以按照常數(shù)來考慮。將井壁看作內外壁存在溫差、兩端自由的厚壁圓筒，根據(jù)彈性力學[5-7]，可得到OB段井壁因為溫度變化以及內外壁溫度差異作用引起的伸長量

(7)

分別將ΔLOB表達式和ΔLt代入式(2)并進行整理可得

(8)

對式(8)進行整理和等量代換，可得

(9)

進一步可以代換為

d1zO3+d2zO2+d3zO+d4=0.

(10)

式(10)為一個關于zO的一元三次方程，利用求根公式可以求得其解析解為

(11)

式中，

d4=3HD1-3u0,D1=

由于溫度升高井筒和土體之間相互作用產(chǎn)生的塑性剪切帶和彈性剪切帶分界面位置zO確定以后，井壁豎向伸長在井壁內產(chǎn)生的溫度應力即可得到。在塑性剪切段(OA段)，根據(jù)式(1)，井壁橫截面上由于溫度產(chǎn)生的豎向附加力在井壁內形成的軸力和豎向附加應力分別為

(12)

(13)

對于彈性剪切段(OB段)橫截面上由溫度產(chǎn)生的軸力豎向附加應力分別為

(14)

(15)

3 實例計算

現(xiàn)以一個井筒為例進行分析，說明溫度豎向附加應力的計算，并通過實例分析各種因素對豎向溫度附加應力的影響。某礦井井筒穿過的表土層厚度為235 m，井壁為C33鋼筋混凝土結構，井筒內半徑a=6.0m，外半徑b=7.0m，泊松比μ=0.2，井壁重度γc=25kN/m3,表土層重度為γs=19kN/m3,井壁外表面溫度Tb=15 ℃，井壁內表面溫度Ta=35 ℃，鋼筋混凝土井壁彈性模量E=3.0×107kPa，井壁熱膨脹系數(shù)αf=10×10-6/℃;u0=0.007m;λ=0.333;k0=0.3。

從前文的公式中可以看出，井壁溫度豎向附加應力的求解關鍵在于確定井壁與土體剪切帶的彈塑性分界深度zO。為此，將相關的數(shù)據(jù)代入到式(11)得到zO=10.5m，然后將其他相關數(shù)據(jù)代入到式(15)即可得到表土與基巖交界處由溫度豎向附加力在井壁內產(chǎn)生的溫度豎向附加應力最大值為σzt2=2.50MPa。溫度豎向附加應力沿井壁軸向的變化規(guī)律如圖2所示。為了分析不同深度彈塑性交界點位置對井壁中溫度豎向附加應力的影響，選取zO分別為2、4、6、8m進行計算，得到了不同的彈塑性分界深度下井壁內溫度豎向附加應力的分布規(guī)律。從圖2可以看出，溫度豎向附加應力都是隨著深度呈增大趨勢，在表土與基巖交界面處達到最大值，這與大部分井壁在此處破裂的工程實踐一致；從每一條曲線來看，溫度豎向附加應力在井筒上部沿深度增加的幅度要大于下部，這與井壁外表面的溫度豎向附加力的分布有關；從圖2還可以看出，隨著彈塑性分界點深度增加，井壁內的溫度豎向附加應力也在增加，并且增加幅度逐漸加大，這說明溫度對井壁的影響在增大，這與溫度升越高，井壁縱向膨脹變形越大，彈塑性交界點的深度就越大的結論一致。圖3為井壁彈塑性交界點深度與井壁最大溫度豎向附加應力關系曲線，從中可以看出，井壁表土層底部與基巖交界處的最大溫度豎向附加應力隨彈塑性交界點的深度近似呈線性增加。

4 結 論

圖2 溫度豎向附加力隨井壁深度變化規(guī)律

圖3 溫度豎向附加應力與彈塑性分界深度關系

(1)在溫度升高時，深厚表土層井壁的變形除了徑向膨脹外，還因縱向伸長與周圍土層形成豎向剪切帶而在井壁內產(chǎn)生溫度豎向附加應力。

(2)表土段井壁只能由基巖向地面伸長，上段井壁和土體間的相對位移大于極限剪切位移，二者之間的剪切帶處于塑性狀態(tài)；下段井壁與土體間的相對位移小于彈性極限位移，剪切帶處于彈性狀態(tài)。

(3)溫度豎向附加應力沿井壁縱向呈不斷增大趨勢，在基巖處達到最大值；彈塑性分界點深度越大，井壁內溫度豎向附加應力的增大幅度也越大。

[1] 崔廣心.特殊地層條件下豎井井壁破壞機理及防治技術[J].建井技術，1998，19(1)：28-32. Cui Guangxin.The failure mechanism of shaft lining under special strata and the technique of preventing the shaft from fracturing[J].Well Construction Technology,1998，19(1)：28-32.

[2] 楊維好，崔廣心，周國慶，等.特殊地層條件下井壁破裂機理與防治技術的研究(之一)[J].中國礦業(yè)大學學報，1996,25(4)：1-5. Yang Weihao，Cui Guangxin,Zhou Guoqing,et al.Fracture mechanism of shaft lining under special strata condition and the technique preventing the shaft from fracturing (part one)[J].Journal of China University of Mining and Technology，1996,25(4)：1-5.

[3] 崔廣心，楊維好，呂恒林.深厚表土層中的凍結壁和井壁[M].徐州：中國礦業(yè)大學出版社，1998:125-129. Cui Guangxin，Yang Weihao，Lu Henglin.Frozen Wall and Shaft Lining in Thick Alluvium[M].Xuzhou：China University of Mining and Technology Press，1998:125-129.

[4] 經(jīng)來旺，高全臣，徐輝東，等.凍結壁融化階段井壁溫度應力研究[J].巖土力學，2004，25(9)：1358-1362. Jing Laiwang，Gao Quanchen，Xu Huidong，et al.Thermal stress analysis of shaft in melting stage of frozen wall[J].Rock and Soil Mechanics，2004，25(9)：1358-1362.

[5] 劉金龍，陳陸望，王吉利.立井井壁溫度應力特征分析[J].巖土力學，2011,32(8)：2386-2390. Liu Jinlong，Chen Luwang，Wang Jili.Characteristic analysis of temperature stresses of shaft wall[J].Rock and Soil Mechanics，2011,32(8)：2386-2390.

[6] 劉金龍，陳陸望，王吉利.考慮溫度應力影響的立井井壁強度設計方法[J].巖石力學與工程學報，2011,30(8)：1557-1563. Liu Jinlong，Chen Luwang，Wang Jili.Design method for shaft walls strength of mine considering influence of temperature stresses[J].Chinese Journal of Rock Mechanics and Engineering，2011,30(8)：1557-1563.

[7] 何朋立.考慮溫度應力的深厚表土層井壁破裂分析[J].工業(yè)建筑，2013，43(3)：75-78. He Pengli.Rupture Analysis of shaft lining in the deep alluvium considering the temperature stress[J].Industrial Construction，2013，43(3)：75-78.

[8] 梁恒昌，周國慶，劉志強.地層沉降和井壁附加應變實測耦合分析[J].煤炭學報，2000，25(1)：40 -45. Liang Hengchang，Zhou Guoqing，Liu Zhiqiang.In-situ tests on ground surface settlement and additional strain of shaft linings[J].Journal of China Coal Society，2000，25(1)：40-45.

[9] 黃定華，俞懷君，馬嘉榮，等.臨渙礦區(qū)地表沉降對井壁破裂的影響[J].煤炭科學技術，1991(7)：42-47. Huang Dinghua，Yu Huaijun，Ma Jiarong，et al.Effects of the surface subsidence on shaft wall fracture in Linhuan mine area[J].Coal Science Technology，1991(7)：42-47.

[10] 李世平.巖石力學簡明教程[M].徐州：中國礦業(yè)大學出版社，1986. Li Shiping.Rock Mechanics Introductory Tutorial[M].Xuzhou：China Mining University Press，1986.

[11] 經(jīng)來旺，張 浩，徐輝東，等.井筒變形機理分析及高強高性能混凝土井壁材料研究[M].合肥：中國科學技術大學出版社，2011. Jing Laiwang，Zhang Hao，Xu Huidong，et al.Shaft Deformation Mechanism Analysis and High Strength High Performance Concrete Shaft Materials Research[M].Hefei：Science and Technology University of China Press，2011.

(責任編輯 徐志宏)

Analysis of Shaft Lining Temperature Vertical Stress in Deep Topsoil

He Pengli1Guo Li1Wang Zaiquan2

(1.SchoolofCivilEngineering，LuoyangInstituteofScienceandTechnology，Luoyang471023，China；2.CollegeofScience，QingdaoTechnologicalUniversity，Qingdao266033，China)

The shaft lining rupture in deep topsoil is a geological disaster，which seriously harm the coal mine production.Temperature action is an important factor that results in the increasing of vertical wellbore stress and induces the fracture of the shaft lining.In order to analyze the temperature effect of longitudinal expansion on shaft lining in deep alluvium，on the basis of analyzing the interaction between shaft lining and soil，the el-plastic shearing model of the soil and soil shearing was built according to thermodynamics and elastic theory，and the calculation formula of temperature vertical additional stress is induced.Further analysis shows that the temperature vertical additional stress in shaft lining increases with depth continuously，and reached the maximum value near the bedrock.The temperature vertical additional stress displays as an approximate linear increase with the depth of el-plastic shear section cut-off point.These findings contributes to deepening the understanding on the mechanism of shaft lining rupture in deep topsoil，and provides the theoretical basis for the prevention and governance of shaft lining fracture accident.

Longitudinal expansion，Temperature action，Interaction between shaft lining and soil mass，El-plastic shearing model，Temperature vertical additional stress

2015-01-26

何朋立(1976—)，男，講師，博士。

TD263

A

1001-1250(2015)-04-119-04