深井巷道圍巖溫度場分布規律研究

黎明鏡

隨著我國煤礦開采深度的不斷增加，其高溫環境引起的熱害問題越來越受到國內外學者和工程技術人員的關注，深井巷道巖體溫度是熱害礦井的最大熱源。巷道圍巖的溫度場分布研究是基礎問題之一，本文基于地質學和熱力學理論，給出巷道圍巖溫度場受不同的風流冷卻時間的無因次解析式，并導出調熱圈外半徑的近似計算式，針對淮南礦區深部開采的高地溫(42℃左右)，采用有限元軟件ANSYS9.0對深井巷道圍巖的溫度分布進行了數值模擬分析。

1 圍巖溫度場分布的解析解

開掘巷道瞬間，巖體內各處溫度都等于原始巖溫 Ty，當恒定風溫為 Tf(Tf＜Ty)的空氣通過后，巖熱就開始流向巷道。在t1時刻內，圍巖壁溫度為 Tb，隨著時刻的繼續增長，圍巖壁溫也不斷接近風溫(如圖1所示)，即limTb→Tk，但嚴格地說，Tb永遠不會絕對等于Tk。當圍巖壁溫極為接近風溫后，調熱圈抵達相對穩定的最大半徑[1]。

新掘進的巷道，圍巖中受風流冷卻的調熱圈不大，這時圍巖可當作半無限體來求出溫度場。當巷道掘通100 d以上時，圍巖受風流冷卻時間較長，調熱圈就會變大到相當程度，這時圍巖的溫度場將受到巷道形狀的影響，就不能把圍巖單純看作是半無限體，要采用其他方法求解[2]。

1.1 受風流冷卻時間短的圍巖溫度場

受風流冷卻時間短的圍巖可當作半無限體來求出溫度場。設圍巖均質、干燥，根據傳熱學有[3]:

1.2 受風流冷卻時間長的圍巖溫度場的確定

因為受風溫冷卻時間長的巷道圍巖與風流熱交換是一個復雜的過程[4]，在熱交換過程中又受到各種因素的影響，為了使條件簡化，便于討論其本質，作出如下假設:1)巷道圍巖均質且各向同性;2)初始溫度等于原始巖溫;3)巷道圍巖壁面換熱條件一致;4)巷道內流動空氣恒溫;5)圍巖傳遞的熱量完全傳給風流。

在以上假設的基礎上，圍巖與風流開始進行熱交換后，圍巖溫度場開始隨時間變化，圍巖的導熱過程是非穩態導熱過程。由熱量守恒定律和傅里葉導熱定律，建立三維熱傳導方程，求得無因次溫度的解析解為(求解過程詳見文獻[5][6]):

根據式(2)即可確定受風流冷卻時間長的巷道圍巖溫度場。

2 計算實例

以淮南朱集礦為實例，用式(1)，式(2)計算其-906 m井底車場巷道溫度場的解析解，用有限元軟件ANSYS9.0計算其數值解。

朱集煤礦-906 m井底車場位置的主要巖性為鋁質泥巖、粉砂巖、中砂巖。該礦區構造復雜、斷層發育，特別是深部斷層的交匯，形成深部地熱上行的良好通道，加之上部的泥巖和泥質砂巖構成熱儲構造的復合蓋層，形成該礦區高溫達42℃左右，具有相當大程度的熱害。數據采集巷道為通風巷道，斷面形狀為半圓拱形，直墻高度1.7 m，半圓拱半徑3.7 m，錨噴支護。

2.1 實測熱物性參數

溫度場計算參數見表1。

表1 溫度場計算參數表

2.2 圍巖溫度場的數值計算

模型的建立取巷道周邊圍巖的溫度場分布，因此，模型的計算范圍根據對稱性取巷道上半部分圓心的上、下各30 m，左、右各40 m(見圖2)，深為-906 m，在二維平面內計算，單元類型為4節點四邊形實體單元(Plane55)，共劃分2 964個單元，計算幾何模型及網格劃分如圖2所示。

建立模型后，進入求解器，設置分析類型為瞬態分析，設定圍巖模型初始溫度分布分別為42℃的均勻溫度場。選取巷道的輪廓線，分別施加20℃的對流荷載，巖壁參考溫度為25℃，將實體模型轉換為有限元模型進行求解。

將同一條件下，按公式計算出的溫度場各離散點上溫度隨時間變化而變化的計算值與數值計算值做對比分析可知，解析解和模擬值基本吻合，解析解的溫度場分布范圍略小于模擬值。

2.3 圍巖溫度場的分布規律

從ANSYS數值計算結果中調出不同時間、不同徑向深度的溫度分布數據制成圖。由圖3，圖4可知，圍巖溫度場一年內最大的變溫區域位于距巷道壁面2.3 m～6 m之間，該處圍巖溫度的變化產生較大的溫度應力和變形，現場松動圈實測結果表明該范圍的圍巖破碎嚴重，形成了相當大的熱害。

3 結語

1)本文建立了熱害礦井巷道溫度場的數學模型，給出巷道圍巖溫度場受不同的風流冷卻時間的無因次解析式，并導出調熱圈外半徑的近似計算式。2)以淮南地區某高溫礦井的實測數據作為參數，用有限元軟件ANSYS計算出圍巖溫度場的分布，解析解與數值模擬解得以相互印證。3)朱集礦-906 m井底車場圍巖溫度場一年內最大的變溫區域位于距巷道壁面2.3 m～6 m之間。該處圍巖溫度的變化產生較大的溫度應力和變形，現場松動圈實測結果表明該范圍的圍巖破碎嚴重，形成了相當大的熱害，計算結果可為支護設計提供參考。

[1] 岑衍強，侯棋棕.礦內熱環境工程[M].武漢:武漢工業大學出版社，1989:50-101.

[2] 張素芬.礦井通風降溫的有效性與局限性[J].煤炭安全，1992(7):78-79.

[3] 趙以蕙.礦井通風與空氣調節[M].北京:中國礦業大學出版社，1996:193-230.

[4] 趙 利，楊德源.礦內風溫預測計算及程序編制[J].煤礦安全，1996(11):25-28.

[5] 周心權，吳 兵.礦井火災救災理論與實踐[M].北京:煤碳工業出版社，1996:300-320.

[6] Scerban，A.N.The Scientific Basis of Calculations and Control of Climatic Conditions in Deep Mines[M].AN Ukraine R.S.R.，Kiev，1959:164-250.

[7] 田錦州.誤差函數erf(x)近似計算及其在開采沉陷預計中的應用[J].煤礦開采，2009，14(2):33-35.