充填料漿自流輸送管道磨損機理研究
2018-05-02 11:13:15毛明發王炳文朱家銳常偉華
金屬礦山 2018年4期
毛明發 王炳文 朱家銳 常偉華
(1.云南馳宏鋅鍺股份有限公司，云南 曲靖 655011；2.中國礦業大學(北京)資源與安全工程學院，北京 100083)
隨著我國淺部資源不斷枯竭，開發深部礦產資源將成為一種趨勢。在這種趨勢的驅動下，集安全性能高、環境保護好、礦石回收率高以及控制地壓強等諸多優點為一體的充填采礦法，尤其是高濃度膠結充填技術，必將成為未來深部礦床開發的主要開采方法之一。但是，隨著礦山開采深度的增加，充填系統的安全性、穩定性以及可靠性等問題成為了以充填采礦技術為主要開采方法礦山的技術瓶頸，特別是膏體充填或高濃度料漿充填系統中出現的高頻率堵管、爆管問題[1],因此，對充填管道的磨損機理進行研究具有重要的意義。韓文亮等[2]分析了物料在管道輸送過程產生真空不滿管流的原因。張欽禮、劉同有、王賢來等[3-6]從動量、能量以及磨蝕理論2個角度闡述了充填鉆孔內管道磨損的原因，引入了管壁單位面積耗能量的概念。付林[7]利用 Fluent數值模擬軟件對質量濃度為47%的油煤漿在彎管部位的流場進行仿真模擬計算。但是，相關充填鉆孔內空氣高速流動加速充填管道氧化腐蝕的報道相對較少；同時，能夠定量給出確定充填管道特殊位置磨損區域大小的計算方法也十分罕見，相關利用ANSYS CFD數值模擬軟件仿真模擬高濃度充填料漿在管道內的流動情況來解釋充填管道磨損的研究也較少。
基于上述分析，本研究在 Finnie[8-9]的微切削理論的基礎上，結合某礦充填管道的磨損情況，從動量、能量、料漿相變理論以及ANSYS CFD軟件數值模擬結果等多角度深入分析充填管道磨損機理。同時，研究各因素對管道磨損的影響程度，以便為充填料漿管道自流輸送系統的設計提供參考。
1 管道自流輸送系統
自流充填系統中，料漿流動的動力來自料漿在垂直管段產生的靜壓頭[10]。如圖1所示，如果充填系統中管徑為恒定值，豎直管高差為H0，水平管長度為L0，料漿體重為γp，摩擦阻力損失ip，則料漿輸送有3種情況：當γp·H0
2 管道磨損機理研究
2.1 充填管道磨損情況
某金屬礦山充填鉆孔深度145 m，管徑146 mm，水平管長度450 m。經過充填鉆孔管道的料漿充填量達到20萬m3后，鉆孔的磨損位置主要表現在3處：①距進料口900 mm處磨損沖蝕嚴重，耐磨層及外管均已被磨穿，如圖2(a)所示；②距鉆孔入口40～90 m段充填管磨損逐漸增大，出現絲扣連接下端耐磨層脫落以及局部耐磨層磨損等現象，如圖2(b)所示；③距進料彎管90 m往后，耐磨層已經磨薄甚至消失，鉆孔直徑變大，充填料將直接沖刷外層，見圖2(c)。另鉆孔與水平管交界處的直角彎管磨損嚴重，出現了窄長的沖刷槽，見圖2(d)。因此，本研究將鉆孔分為入口段、空氣—料漿交界面段和鉆孔底部直角彎管段3個區段來深入研究鉆孔的磨損機理。
圖1 充填料漿管道自流輸送模型Fig.1 Pipeline gravity transportationmodel of filling slurry
圖2 充填管道不同管段的磨損情況Fig.2 Wear patterns of filling pipe on different section
2.2 鉆孔入料口處磨損機理研究
鉆孔入口處充填料漿的運動軌跡如圖3所示。料漿通過彎管進入鉆孔入口時具有一定的初速度，因此料漿在管道內將做拋物線運動。當料漿碰撞到管壁時，將對管壁進行沖擊和摩擦雙重作用，同時，速度方向急劇變化，產生了巨大的沖量和能量損失，這部分能量損失主要用于對管道的磨損作用。鉆孔入口處的沖蝕磨損長度為C，管道直徑為D。根據能量、動量、微切削理論[8-9]以及充填鉆孔入口的三角幾何關系推導出C=D，鉆孔入口的最長磨損距離h的計算式、料漿在B點的入射角β和入射速度vB以及體積沖蝕率VB的數學表達式分別為
圖3 鉆孔入口處充填料漿的運動軌跡Fig.3 Movement path of filling slurryat the entrance of drilling
(1)
(2)
(3)
(4)
式中，M為顆粒的質量；K為顆粒分數；p為顆粒與耙材間的彈性流動壓力；g為自由落體加速度。
其中，
由以上各式可知，最長磨損距離的位置與斜管的坡度、質量料漿流量(Q，kg/s)和管道直徑等因素有關，且具有隨著管徑、斜管坡度的增大而增大，隨著流量的升高而降低的特性。
2.3 鉆孔內管道磨損機理的研究
2.3.1 沖擊磨損位置的確定
空氣～料漿交界面處的沖擊磨損是鉆孔管道磨損最為嚴重的位置之一，可以通過推導出交界面距鉆孔入口的高度計算表達式來準確預測出沖擊磨破位置。若理想地認為自由段料漿的剩余能量僅用管道磨損，那么交界面距鉆孔入口的高度數學表達式為
(5)
式中，H0為垂直管高度，m；L0為水平管長度，m；γp為料漿體重，t/m3；ip為料漿的水力坡度，kPa/m(假設系統管徑相同)。
但是，在充填料漿的非滿管流輸送系統中，料漿自由下落高度要高于式(5)所計算的結果。因為單位質量的料漿在自由下落段內所具有的位能并不全部用于管道磨損。因此，引入能量分配系數f，它表示單位質量的料漿用于克服管道磨損所消耗的能量與其在自由下落段所具有的最大位能的比值。以空氣—料漿交界面為基準面，單位質量的料漿在自由段所具有的勢能Ez主要消耗在2個方面：一部分能量用于管道磨損，記為Em；另一部分能量用于克服料漿在滿管流段的相應距離Hp的水力坡降，記為Ed。
Ez=mgH1/mg=H1,
(6)
Em=f·Ez,
(7)
Ed=(1-f)·Ez.
(8)
根據水力學理論原理推導出充填料漿非滿管流輸送過程中自由下落段的表達式為
(9)
由式(9)可以看出，Hs的大小由水平管道的長度L0、水力坡度ip、鉆孔高度H0以及鉆孔管道能量磨損分配系數f共同決定。自由下落段高度H1隨水平管道的長度增大而減小，同時沖擊磨破的位置距離鉆孔入口越近，實際自由下落高度的幾何模型如圖4。
圖4 空氣—料漿交界面的實際位置Fig.4 Actual location of air-slurry interface
2.3.2 基于動量理論研究管道磨損
2.3.2.1料漿對空氣—料漿界面的沖擊力
由于充填鉆孔管道內空氣—交界面是具有一定速度的動態面。因此，根據相關理論推導出料漿對交界面的沖擊力
(10)
式中，k為試驗修正系數；其他符號意義同上。
2.3.2.2 充填料漿動量分析
通過建立如圖5所示的物理模型，假設充填料漿進入豎直管口的初始速度為v1，料漿恰好到達空氣—料漿交界面處的最大瞬時速度為v2，滿管流段的料漿流速為v，管口至交界面之間的距離為H1。
圖5 充填鉆孔料漿內沖擊物理模型Fig.5 Impact physical model offilling slurry in the drilling
(1)充填料漿速度的變化。假設充填料漿在豎直管道內的自由下落段做加速直線運動。根據牛頓定律以及能量守恒方程等理論推導出料漿達到交界面的最大速度為
(11)
(2)充填料漿動量變化。當充填料漿下落到交界面速度達到最大瞬時速度v2時，其動量達到最大值。由于料漿在交界面發生碰撞，料漿速度在極短時間由v2減小到v，因此，充填料漿將會對管壁產生巨大的沖量和沖擊力，加速管道的局部磨損，甚至可導致管道的磨破、破裂。
充填料漿滿管流段的流速為v，那么在沖擊時間Δt內，流過管道斷面的料漿總質量
Δm=ρj·v2·Δt·S,
(12)
式中，ρj為料漿密度，g/cm3；S為充填管道斷面面積，m2。
根據動量定理，充填料漿對管壁的沖力
(13)
將式(12)代入式(13)得：
F=ρj·S·v2(v2-v).
(14)
從式(14)可以看出充填料漿在空氣—料漿交界面處對管壁的沖擊力F與速度v2和速度差v2-v分別成正比。因此，充填料漿的速度是影響沖擊力最主要因素，同時也說明料漿自由下落帶高度越大，流速越大，其對管壁的沖擊磨損也就越嚴重。
2.3.3 基于能量理論研究管道磨損
2.3.3.1 局部管道磨損
充填料漿從管口運動到空氣—料漿交界面消耗的總能量等于進入管口的初始動能加上重力所做的功。假設單位時間Δt內進入管道的料漿的單位質量Δm=ρj·Δt·Q1，Q1為料漿體積流量，m3/s。這部分料漿從管口運動到交界面消耗的能量分為2部分：充填料漿由v1加速到v2過程中的能量損失為Ek1；由v2減小到v過程中能量損失為Ek2。
(15)
式(15)能夠充分解釋解管道偏斜率對管道磨損影響。在一個相對穩定的充填系統中，式(15)中唯一的變量就是v2。管道的偏斜率越高，自由下落段料漿與管壁的接觸面積就越大，造成管壁對料漿向下運動的阻力越大，因此料漿對管壁所做的功就越多，即Ek1越大。因此，局部的管道磨損就會越嚴重。
(16)
這部分能量主要消耗在空氣—料漿交界面附近的管壁上。因此可用管壁單位面積的耗能量Eh作為反映管壁磨損重要指標，Eh越大，料漿對管壁的磨損越嚴重。
(17)
式中，Sh為料漿和管壁有效的磨損接觸面積，m2。Sh=π·D·v·Δt，代入(17)得：
(18)
式(18)能夠充分詮釋為什么空氣—料漿交界面附近的管壁磨損最嚴重。在一個既定的充填系統中，式(18)僅有v2為變量，管壁單位面積的耗能量Eh與v2的二次方成正比例關系。由此可見，非滿管流的存在是造成豎直管道或鉆孔壽命短的首要原因。因此，在充填料漿自流輸送管網設計以及礦山應用過程中，我們都要以最大限度地實現滿管流輸送為宗旨，從而將管道磨損降到最低。
2.3.3.2 整體管道磨損
對于充填管網的整個管道系統來說，管壁單位面積的耗能量Eh可以表示為
(19)
式中，Sq為整個管網的管壁面積，m2;Sq=π·D(H0+L0);α′為修正系數，0<α′≤1。
從式(19)可以看出單位質量料漿與管壁有效磨損的接觸面積越小，單位面積的耗能量就越大，導致管壁磨損越嚴重；反之亦然。因此，從整個充填管網系統來看，充填料漿只有在滿管輸送的條件下才能實現單位面積耗能量的最小，這樣才能將管道的磨損量降到最低。
2.3.4基于鉆孔管道內氣體流動形式研究管道磨損
鉆孔內料漿自由下落區分為3個區域，即空化區、空氣區和水躍區[2]，如圖6。在垂直豎管內料漿首先由連續的漿體柱分離成小段體，然后離散成團狀和顆粒狀，呈散射流形式向下加速運動[13]。相鄰2個段體之間產生真空度。在空氣區由于離散漿體的高速下落，離散漿體頂部產生負壓，底部產生高壓，加速空氣在管孔內的循環流動，導致管道的氧化腐蝕嚴重。在這不斷循環的正負壓強的作用下，將造成鉆孔內自由下落段絲扣連接處下端的耐磨層大面積脫落。再者，微氣泡重新凝結過程中將放出一定熱量，進一步管壁的化學腐蝕創造條件，從而加速管道損壞。
圖6 深井充填管道空氣柱分區示意Fig.6 Partition map of air column indeep mine filling pipeline
2.4 鉆孔下端直角彎管磨損機理研究
2.4.1 鉆孔下端直角彎管壓力
單段自流輸送系統中料漿對垂直管壓力可表示為
(20)
式中，β1為壓力修正系數；其他符號意義同前。
由式(20)可見，料漿對垂直管內的壓力隨著H的增大而增大，當H=H0時，壓力取最大值Pmax。
(21)
2.4.2 料漿對鉆孔下端直角彎管的沖擊力
在充填系統開機啟動階段，料漿在鉆孔內做自由落體運動。料漿運動到直角彎管處時，其速度達到最大值，法向速度將在極短時間內減小至零。因此，將對直角彎管管壁產生很大的沖擊力，彎管處料漿的速度變化如圖7所示。
圖7 啟動階段料漿在彎管處的料漿流速Fig.7 Slurry flow velocity at the bend in startup stage
彎管處沖擊面的法向沖擊力大小以及距彎管入口沖擊點位置可表示為
(22)
L=ξ·Rsinα1.
(23)
(24)
(25)
從上式不難看出，最大沖擊力隨著垂直管高度的增大而增大，隨著曲率半徑的增大而減小。
2.4.3 基于ANSYS CFD數值模擬分析磨損機理
通過應用 ANSYS FLUENT的前處理軟件 DM 按照充填管道滿管流段的直角彎管進行幾何建模。充填料漿質量濃度76%，灰砂比1∶6，料漿密度1.90 t/m3，體積流量90 m3/h。采用速度進口，壓力出口來對管道充填料漿流動模型設置邊界條件，壁面采用無滑移壁面。從而模擬得出充填料漿在彎管處的速度分布矢量圖和剪切應力分布圖。
(1)速度分布。從充填管道XY切平面的局部速度矢量圖8可以看出，料漿在彎管處流速分布發生了變化，即邊界層變薄，彎管底部更加接近核流區(高速區)，從而加速了管道的磨損。
(2)剪切應力分布。充填管壁處的剪切應力直觀地反映了管道容易磨損位置。由圖9可知，彎管處的最易磨損位置為彎管的底部。這主要是由于料漿
圖8 彎管處XY切面的速度分布Fig.8 Velocity distribution of XY plane in bend pipe
在靠近彎管入口附近的底部流速方向發生了改變，從而對管底產生了一定的力，導致該處的剪切應力增大。
圖9 彎管壁面的剪切應力分布Fig.9 Shear stress distribution of the bend pipe wall
2.5 不同因素對管道磨損的影響程度
首先選擇空氣—料漿交界面的沖擊力、滿管率、啟動時垂直管地底部彎管處的最大沖擊力以及輸送過程中垂直管與水平管交接處的壓力來作為表征充填管道磨損的重要指標。然后，利用控制變量法來研究管徑、體積濃度、充填倍線、垂直管高度以及固體顆粒的平均粒度等因素分別對磨損指標的影響程度。其中不同影響因素的變化情況如表1所示。從圖10～圖13可以看出，對充填管道磨損顯著的因素主要為管道直徑、垂直管高度以及充填倍線，其他因素對管磨的影響相對較小。
表1 不同影響因素的變化情況Table 1 Variety of different influence factors
圖10 各因素對滿管率的影響程度Fig.10 Influence of various factorson the full pipe ratio■—密度；▲—管徑；?—粒級組成；◆—體積濃度；▼—充填倍線；●—垂直管高度
圖11 各因素對彎管處沖擊力的影響程度Fig.11 Influence of various factorson the impact on bend pipe■—密度；▲—管徑；?—粒級組成；◆—體積濃度；▼—充填倍線；●—垂直管高度
圖12 各因素對彎管處壓強的影響程度Fig.12 Influence of various factorson the pressure on bend pipe■—密度；▲—管徑；?—粒級組成；◆—體積濃度；▼—充填倍線；●—垂直管高度
3 結 論
通過詳細分析鉆孔損傷探測錄像資料以及現場記錄管道破損數據，結合充填料漿在充填系統中的流態分析結果，提出了計算垂直鉆孔管道中非滿管流段高度的新方法。建立了料漿對鉆孔入口以及底部直角彎管的沖擊力與磨損量的物理力學模型，確定了相關磨損位置的磨損區域大小。從動量、能量以及料漿相變理論等多角度探討了充填管道的磨損機理，優化了管壁單位面積耗能量的計算式，提出了垂直鉆孔自由下落段空氣區內氣體循環流動概念，并建立氣體流動物理模型。利用控制變量法，定量分析了不同影響因素對管道磨損的影響程度。結果表明，對充填管道磨損顯著的因素主要為管道直徑、垂直管高度以及充填倍線，其他因素對管磨的影響相對較小。
圖13 各因素對空氣—料漿交界面沖擊力的影響程度Fig.13 Influence of various factorson impact on interface of air-slurry■—密度；▲—管徑；?—粒級組成；◆—體積濃度；▼—充填倍線；●—垂直管高度
[1] 孫恒虎,段鴻杰,劉文永,等.膏體充填自流輸送系統特征參數分析[J].有色金屬,2000，52(3)：7-10.
Sun Henghu,Duan Hongjie，Liu Wenyong,et al.Analysis of characteristic parameters of self-flowing pastefill systems[J].Nonferrous Metals，2000，52(3)：7-10.
[2] 韓文亮,張志平.長距離輸送管道中的真空不滿流及其預防[J].金屬礦山,1994(11):48-52.
Han Wenliang，Zhang Zhiping.Vacuum-unfull flow in the long-distance transmission pipe line and its prevention[J].Metel Mine,1994(11):48-52.
[3] 王賢來,鄭晶晶,張欽禮,等.充填鉆孔內管道磨損的影響因素及保護措施[J].礦冶工程,2009(5):9-12.
Wang Xianlai，Zheng Jingjing，Zhang Qinli,et al.Influencing factors on the abrasion of pipelines in backfill drilling and protective measures[J].Mining and Metallurgical Engineering,2009(5):9-12.
[4] 劉同有.充填采礦技術與應用[M].北京:冶金工業出版社,2001.
Liu Tongyou.Filling Mining Technology and Application[M].Beijing：Metallurgical Industry Press，2001.
[5] 張德明,王新民,鄭晶晶,等.深井充填鉆孔內管道磨損機理及成因分析[J].武漢理工大學學報,2010(13):100-105.
Zhang Deming，Wang Xinmin，Zheng Jingjing,et al.Wear mechanism and causes of backfilling drill-holes pipelines in deep mine[J].Journal of Wuhan University of Technology,2010(13):100-105.
[6] 張欽禮,鄭晶晶,王新民,等.充填管道磨損形式及機理分析研究[J].金屬礦山,2009(5):115-118.
Zhang Qinli，Zheng Jingjing，Wang Xinmin,et al.A study of the abrasion form and mechanism of backfilling pipelines[J].Metel Mine,2009(5):115-118.
[7] 付 林,高炳軍.油煤漿輸送管道彎管部位流場的數值模擬與磨損預測[J].化工機械,2009(5):463-466.
Fu Lin，Gao Bingjun.Numerical simulation of the flow field and wear prediction of the elbows of coal-oil slurry transporting pipelines[J].Chemical Engineering & Machinery,2009(5):463-466.
[8] Finnie J.Erosion of surfaces by solid particles oberflachen erosion durch feste teilchen[J].Wear,1960,3(2)：87-103.
[9] Finnie J,Kabil Y H.Some observations on the erosion of ductile metals[J].Wear,1965,8(l):60-69.
[10] 何哲祥,古德生.礦山充填管道水力輸送研究進展[J].有色金屬,2008,60(3):116-120.
He Zhexiang,Gu Desheng.Development of research on hydraulic transport of mine filling in pipelines[J].Nonferrous Metals，2008,60(3):116-120.
[11] 王新民.基于深井開采的充填材料與管輸系統的研究[D].長沙:中南大學,2005.
Wang Xinmin.A Study of Filling Materials and Pipeline Transportation Systems in Deep Mines[D].Changsha：Central South University，2005.
[12] Zhang Qinli,Hu Guanyu,Wang Xinmin.Hydraulic calculation of gravity transportation pipeline system for backfill slurry[J].Journal of Central South University of Technology,2008,15(5):645-649.
[13] 王海瑞.全尾砂膠結充填自流輸送管路改造及優化[D].長沙：中南大學,2010.
Wang Hairui.The Alteration and Optimization of Cemented Whole Tailings Gravity Transport Backfilling System[D].Changsha：Central South University，2010.
