王喬磊

(河南理工大學(xué)能源科學(xué)與工程學(xué)院，河南 焦作 454003;霍州煤電呂臨能化煤電綜合項(xiàng)目部，山西 霍州 033200)

國(guó)外采用凍結(jié)法鑿井較多的國(guó)家有英國(guó)、德國(guó)、波蘭、前蘇聯(lián)等，這些國(guó)家技術(shù)相對(duì)成熟，鑿井深度均超450 m，但它們多為含水量多、地壓大的基巖地區(qū)，然而深井凍結(jié)施工工程，特別是對(duì)深厚沖積層進(jìn)行施工不多見(jiàn)［1］。

本文采用ANSYS軟件對(duì)霍州煤電龐龐塔礦井凍結(jié)法施工進(jìn)行模擬，其中，施工過(guò)程中每個(gè)節(jié)點(diǎn)的溫度及熱物理參數(shù)，利用有限元分析法便可進(jìn)行計(jì)算。

1 多排管凍結(jié)情況下溫度場(chǎng)規(guī)律的研究

對(duì)以單、雙排管凍結(jié)溫度場(chǎng)的研究為基礎(chǔ)發(fā)展起來(lái)的多排管凍結(jié)溫度場(chǎng)進(jìn)行研究時(shí)［2］，首先需要將管道在不同地點(diǎn)布置，然后分析中排與外1排之間、外1排與最外排之間等位置的溫度場(chǎng)與鹽水溫度、間距、排距間等影響因素之間的關(guān)系，得到影響凍結(jié)壁形成的上述各因素的規(guī)律，進(jìn)而引導(dǎo)設(shè)計(jì)。

1.1 ANSYS環(huán)境下凍結(jié)區(qū)域的網(wǎng)格劃分

網(wǎng)格劃分以平面三角形為單元的形式進(jìn)行，將網(wǎng)格最密集的地方布置在凍結(jié)管區(qū)域，預(yù)測(cè)的凍結(jié)區(qū)較密，而較稀疏的地方布置在外邊界處。凍結(jié)溫度場(chǎng)計(jì)算網(wǎng)格劃分圖見(jiàn)圖1。

1.2 多排管凍結(jié)溫度場(chǎng)相互影響規(guī)律研究

凍結(jié)管周圍的土性參數(shù)、布置方案對(duì)多排管凍結(jié)溫度的計(jì)算有著重要影響，對(duì)于各自的排距、間距等各取5個(gè)水平，其所取的排距、間距及其對(duì)應(yīng)的鹽水溫度、導(dǎo)熱系數(shù)對(duì)照表見(jiàn)表1。

圖1 凍結(jié)溫度場(chǎng)計(jì)算網(wǎng)格劃分圖

表1 計(jì)算中各因素水平取值表

2 斜井底部?jī)鼋Y(jié)溫度場(chǎng)的模擬

斜井凍結(jié)段的溫度變化規(guī)律可總結(jié)為兩點(diǎn):一是深度方向向下發(fā)展的凍結(jié)壁厚度遠(yuǎn)小于向上發(fā)展的凍結(jié)壁厚度;二是排外遠(yuǎn)小于排內(nèi)的發(fā)展速度［3］。

凍結(jié)管的布置方式、含水率、土層的熱物理參數(shù)、鹽水溫度、初始地溫等因素明顯影響著底部?jī)鼋Y(jié)的溫度場(chǎng)。為研究典型參數(shù)下底部?jī)鼋Y(jié)溫度場(chǎng)的發(fā)展情況，下面以管道下部的土層與凍結(jié)段的一半為基礎(chǔ)，在ANSYS環(huán)境中建立模型來(lái)分析。

2.1 計(jì)算模型的建立

以布置管道所在區(qū)域的下半部分和下部土層為基礎(chǔ)建立模型，建好的幾何模型見(jiàn)圖2［4］。

圖2 幾何模型圖

在ANSYS環(huán)境中運(yùn)用掃掠命令對(duì)模型進(jìn)行劃分，注意應(yīng)先對(duì)模型的頂面以平面單元?jiǎng)澐郑钡綕M足要求;然后將完成網(wǎng)格劃分的頂面作為“源”面，沿垂線向下劃分成體單元［5］。為提高計(jì)算精度應(yīng)將布置凍結(jié)管的地方加密，劃分的整體網(wǎng)格圖見(jiàn)圖3，其局部放大圖見(jiàn)圖4。

圖3 數(shù)值模擬網(wǎng)格整體圖

圖4 數(shù)值模擬網(wǎng)格局部圖

2.2 計(jì)算結(jié)果的分析

為了研究底部?jī)鼋Y(jié)溫度場(chǎng)的發(fā)展規(guī)律，在ANSYS軟件的通用后處理命令中定義了9條路徑。路徑具體位置示意圖見(jiàn)圖5。

圖5 路徑的具體位置圖

對(duì)得到的結(jié)果從時(shí)空兩方面進(jìn)行分析，發(fā)現(xiàn)凍結(jié)壁的發(fā)展是參差不齊的，具體情況如下:

1)在凍結(jié)前期時(shí)，凍結(jié)管同一排軸、界面交點(diǎn)下部?jī)鼋Y(jié)壁厚度小于下部?jī)鼋Y(jié)壁厚度，這是因?yàn)閮鼋Y(jié)初期凍結(jié)壁還來(lái)不及交圈，所以軸、界面交點(diǎn)下部?jī)鼋Y(jié)壁厚度為0。

2)在凍結(jié)中期時(shí)，同排凍結(jié)管下部?jī)鼋Y(jié)壁發(fā)展較慢，其發(fā)展速度遠(yuǎn)小于軸、界面交點(diǎn)下部?jī)鼋Y(jié)壁厚度。

3)由圖5可知，在凍結(jié)中后期中排管軸、界面交點(diǎn)下部?jī)鼋Y(jié)壁厚度小于中排凍結(jié)管下部?jī)鼋Y(jié)壁厚度;同排管軸、界面交點(diǎn)下部?jī)鼋Y(jié)壁厚度大于外1排和最外排管下部?jī)鼋Y(jié)壁厚度。

3 不同凍結(jié)壁厚度的力學(xué)與變形特征模擬

3.1 力學(xué)計(jì)算模型的建立

本文根據(jù)凍結(jié)施工的實(shí)際條件建立數(shù)學(xué)計(jì)算模型，其周圍邊界設(shè)定的條件如下［6］:

1)將煤層上覆巖層的重量μ定義為模型的上邊界，其計(jì)算公式如下:

式中:

η—上覆巖層每立方所受力，kN/m3;

L—煤層的厚度，m。

2)將底板定義為模型的下邊界，并將底板簡(jiǎn)化為x方向能夠移動(dòng);y方向用固定鉸支座固定的形式，即y方向速度為0。

3)將實(shí)體煤巖體定義為模型的兩側(cè)邊界，并將實(shí)體煤巖體簡(jiǎn)化為y方向能夠移動(dòng);x方向用固定鉸支座固定的形式，即x方向速度為0。建立的力學(xué)計(jì)算模型見(jiàn)圖6。

圖6 力學(xué)計(jì)算模型圖

3.2 計(jì)算方案分析

本文采用以下三種方案分析了厚度對(duì)溫度場(chǎng)的影響，其方案如下［7］:

方案1:無(wú)凍結(jié)施工;方案2:凍結(jié)壁的平均凍結(jié)溫度為－9℃，凍結(jié)厚度為3 m;方案3:凍結(jié)壁的平均凍結(jié)溫度為－9℃，凍結(jié)厚度為5 m。

3.3 計(jì)算結(jié)果分析

在無(wú)凍結(jié)施工情況下，其頂板下沉量為1.4 m左右;當(dāng)凍結(jié)厚度為3 m，其頂板下沉量為0.94 m左右;當(dāng)凍結(jié)厚度為5 m，其頂板下沉量為0.73 m左右，凍結(jié)壁厚度與頂板下沉關(guān)系曲線見(jiàn)圖7。

圖7 凍結(jié)壁厚度與頂板下沉關(guān)系曲線示意圖

在無(wú)凍結(jié)施工情況下，其圍巖垂直應(yīng)力為0.77 MPa;當(dāng)凍結(jié)厚度為3 m，其圍巖垂直應(yīng)力為1.52 MPa;當(dāng)凍結(jié)厚度為5 m，其圍巖垂直應(yīng)力為1.22 MPa，凍結(jié)壁厚度與圍巖垂直應(yīng)力關(guān)系曲線見(jiàn)圖8。

圖8 凍結(jié)壁厚度與圍巖垂直應(yīng)力關(guān)系曲線圖

ANSYS環(huán)境下三種工況下的井筒變形云圖見(jiàn)圖9～11。

圖9 無(wú)凍結(jié)壁時(shí)井筒變形云圖

圖10 凍結(jié)壁厚度為3 m時(shí)井筒變形云圖

圖11 凍結(jié)壁厚度為5 m時(shí)井筒變形云圖

4 結(jié)語(yǔ)

本文遵循由簡(jiǎn)入繁的研究思路，對(duì)影響凍結(jié)溫度場(chǎng)的因素進(jìn)行了分析，得出了以下主要結(jié)論:

1)底部?jī)鼋Y(jié)溫度場(chǎng)從時(shí)空兩方面上來(lái)說(shuō)都是不均勻的，各個(gè)位置的凍結(jié)壁平均溫度、凍結(jié)壁向下發(fā)展的厚度都不同。

2)影響凍結(jié)溫度場(chǎng)的三個(gè)主要因素分別為:鹽水溫度、凍結(jié)管間距與凍結(jié)時(shí)間。

上述針對(duì)斜井凍結(jié)施工過(guò)程的模擬分析，過(guò)程較為詳細(xì)，可為今后凍結(jié)施工領(lǐng)域的研究與探索提供一定的幫助。

［1］蘇樂(lè)斌.斜井凍結(jié)技術(shù)在龐龐塔礦千萬(wàn)噸項(xiàng)目部斜井施工中的應(yīng)用［J］.中國(guó)煤炭工業(yè)，2012(7):56－59.

［2］李 林，徐兵狀，趙根全，等.凍結(jié)法鑿井中局部?jī)鼋Y(jié)技術(shù)對(duì)已成井壁保護(hù)的溫度場(chǎng)分析［J］.煤炭工程，2012(2):27－29.

［3］王曉健，程 樺，姚直書(shū)，等.特厚沖積層礦井豎向局部保溫凍結(jié)技術(shù)研究與應(yīng)用［J］.煤炭工程，2009(7):70－80.

［4］薛玉庭，李向陽(yáng)，陳陽(yáng)光.井筒凍結(jié)法施工的常見(jiàn)問(wèn)題及防治措施［J］.中國(guó)西部科技，2009(10):56－59.

［5］許舒榮，李長(zhǎng)忠，趙根全，等.快速高效凍結(jié)技術(shù)在立井井筒凍結(jié)施工中的應(yīng)用［J］.建井技術(shù)，2007，28(5):17－19.

［6］靳巍巍，陳有亮.隧道凍結(jié)法施工三維有限元溫度場(chǎng)及性狀分析［J］.地下空間與工程學(xué)報(bào)，2007(05):35－38.

［7］王祥玉.凍結(jié)法施工的立井井壁結(jié)構(gòu)探討分析［J］.煤炭技術(shù)，2008(07):40－43.