旋塞球閥是鉆柱內(nèi)防噴系統(tǒng)中的關(guān)鍵設(shè)備，在氣固兩相流下球閥易受磨損而失效，并造成嚴(yán)重的井噴事故。為此，將計(jì)算流體動(dòng)力學(xué)理論與沖蝕磨損理論相結(jié)合，運(yùn)用FLUENT軟件對(duì)球閥壁面在氣固兩相流下的磨損分布情況進(jìn)行研究，并進(jìn)一步分析了球閥結(jié)構(gòu)參數(shù)對(duì)于球閥壁面磨損的影響規(guī)律。

　　旋塞球閥是鉆柱內(nèi)防噴系統(tǒng)中的關(guān)鍵設(shè)備，在氣固兩相流下球閥易受磨損而失效，并造成嚴(yán)重的井噴事故。為此，將計(jì)算流體動(dòng)力學(xué)理論與沖蝕磨損理論相結(jié)合，運(yùn)用FLUENT軟件對(duì)球閥壁面在氣固兩相流下的磨損分布情況進(jìn)行研究，并進(jìn)一步分析了球閥結(jié)構(gòu)參數(shù)對(duì)于球閥壁面磨損的影響規(guī)律。結(jié)果表明：當(dāng)氣固兩相流流經(jīng)球閥時(shí)，固體顆粒會(huì)與氣流分離，并在壁面上產(chǎn)生三處磨損集中區(qū);隨著球閥開度的減小，球閥壁面磨損量會(huì)急劇增大，且閥球內(nèi)通道壁面上的磨損集中區(qū)由塊狀逐步轉(zhuǎn)化為帶狀，而球閥出口處的磨損集中區(qū)則會(huì)逐漸向下移動(dòng);球閥流道直徑的減小也會(huì)使得壁面磨損量增加，但磨損集中區(qū)的分布基本不變。真空技術(shù)網(wǎng)(http://www.mp99x.cn/)認(rèn)為研究結(jié)果可為進(jìn)一步優(yōu)化球閥流道結(jié)構(gòu)以減輕其壁面磨損提供理論依據(jù)。

0、引言

　　在石油天然氣工業(yè)中，旋塞球閥常被用于鉆柱內(nèi)防噴系統(tǒng)，以預(yù)防和處理鉆柱內(nèi)井噴，但近年鉆井技術(shù)的充分發(fā)展使得鉆具的工作環(huán)境越來越惡劣，鉆柱常受到含有固體雜質(zhì)的氣流或液流對(duì)其的沖蝕而發(fā)生磨損，而旋塞球閥作為鉆柱的關(guān)鍵設(shè)備，其在氣固兩相流下受磨損而失效的現(xiàn)象尤為嚴(yán)重。

　　閥門的磨損問題廣泛存在于含固相顆粒輸送的工程中，因此許多學(xué)者針對(duì)閥門的這一問題進(jìn)行了試驗(yàn)研究。Nkleberg等對(duì)石油天然氣工業(yè)上常用的針形節(jié)流閥進(jìn)行了關(guān)于沖蝕磨損問題的試驗(yàn)研究，觀察并分析了閥門各處在氣固兩相流下的磨損情況。Wheeler等對(duì)近海閘閥在不同內(nèi)涂層下的壁面磨損情況進(jìn)行了試驗(yàn)研究，并觀察了受到磨損后閘閥壁面的微觀形態(tài)。Fang等針對(duì)電液伺服閥的沖蝕磨損問題建立了物理失效模型，并對(duì)模型的準(zhǔn)確性進(jìn)行了實(shí)驗(yàn)驗(yàn)證。但是由于實(shí)驗(yàn)研究常受到現(xiàn)場(chǎng)條件的限制，且通過實(shí)驗(yàn)只能得到閥門磨損的總體質(zhì)量損失，但得不到具體的磨損集中區(qū)單位面積上的磨損速率，而實(shí)際生產(chǎn)中閥門的失效形式常為局部的壁面減薄和穿透，因此越來越多的專家學(xué)者嘗試?yán)�用�?shù)值模擬方法對(duì)各類閥門在含固相顆粒輸送過程中的沖蝕磨損問題進(jìn)行研究。Forder等在CFD軟件中建立沖蝕磨損量計(jì)算模型，并據(jù)此對(duì)控制閥的磨損情況進(jìn)行了預(yù)測(cè)。Atkinson等自行設(shè)計(jì)了應(yīng)用于石油天然氣工業(yè)的抗蝕閥流體通道，并利用CFD軟件分別對(duì)其在含砂氣流下的磨損情況進(jìn)行研究。Zhu等應(yīng)用Fluent軟件對(duì)針形閥在氣固流下的磨損特性進(jìn)行了研究，并對(duì)閥芯進(jìn)行了流固耦合條件下的位移模擬。

　　由于鉆柱內(nèi)球閥隨鉆桿在地層中工作，因此對(duì)其開展現(xiàn)場(chǎng)實(shí)驗(yàn)較為困難。而為了防止發(fā)生因球閥受到磨損而失效所導(dǎo)致的井噴事故，利用數(shù)值模擬方法對(duì)球閥進(jìn)行氣固兩相流下的磨損特性研究很有必要。

1、球閥物理模型及模擬條件確定

　　流道直徑為80mm的球閥在任意開度下的結(jié)構(gòu)及其尺寸如圖1所示，為使湍流充分發(fā)展，在球閥前后增加了一定長(zhǎng)度的直管段。并以開度50%，流道直徑80mm作為球閥的標(biāo)準(zhǔn)狀態(tài)。

圖1 任意開度下球閥流道結(jié)構(gòu)示意圖

　　根據(jù)生產(chǎn)中的實(shí)際情況確定模擬條件，本研究中所選用的連續(xù)相介質(zhì)為甲烷，考慮到流域進(jìn)出口溫度變化較小且氣體流速不高，故將連續(xù)相介質(zhì)設(shè)置為不可壓縮流體，其密度為0.6679kg/m3，動(dòng)力粘度為1.087×10-5Pa·s;固體顆粒簡(jiǎn)化為球形，其直徑為100μm，密度為2700kg/m3;而球閥為碳鋼球閥。

2、數(shù)學(xué)模型的建立

　　2.1、氣體控制方程及RNGk-ε模型

　　氣體的連續(xù)性方程及動(dòng)量方程分別為：

　　(1)

　　(2)

　　式中，ρ為流體密度，kg/m3;t為時(shí)間，s;ui、uj是流體時(shí)均速度分量，m/s;p為流體微元體上的壓力，Pa;μ為流體動(dòng)力粘度，Pa·s;Si為動(dòng)量方程的廣義源項(xiàng)。

　　為使方程組封閉，必須選用合適的湍流模型。RNGk-ε模型作為標(biāo)準(zhǔn)k-ε模型的改進(jìn)模型，可以更好地處理流線彎曲程度較大的流動(dòng)，包括旋轉(zhuǎn)均勻剪切流、包含有射流和混合流的自由流動(dòng)、邊界層流動(dòng)等，因此本研究中選用RNGk-ε湍流模型。

　　2.2、離散相顆粒軌道計(jì)算模型

　　Fluent中的離散相模型(DPM)只考慮氣體與顆粒之間的相互作用，而不考慮顆粒間的相互作用，故可用來模擬湍流中稀疏顆粒相的運(yùn)動(dòng)。運(yùn)用離散相模型進(jìn)行模擬時(shí)要控制顆粒體積分?jǐn)?shù)小于10%，但其質(zhì)量承載率可以大于10%。

　　基于拉格朗日法的離散相顆粒作用力平衡微分方程在x方向上的形式為：

　　(3)

　　式中，up為離散相顆粒的速度，m/s;u為流體流速，m/s;t為時(shí)間，s;FD(u-up)為顆粒在單位質(zhì)量下的曳力;gx為重力加速度在x方向上的分量，m/s2;ρp為顆粒密度，kg/m3;ρ為流體密度，kg/m3;Fx為單位質(zhì)量下其它作用力在x方向上的分量。

　　2.3、磨損量計(jì)算模型

　　Edwards等研究了固體顆粒對(duì)碳鋼和鋁的沖刷磨損，并在綜合考慮了多種因素的基礎(chǔ)上提出了一種磨損速率計(jì)算模型，而考慮到本研究中所選用的球閥材料為碳鋼，故選用Edwards模型作為磨損量計(jì)算模型，其表達(dá)式為：

　　(4)

　　式中，Rerosion為壁面的磨損速率，kg/(m2·s);Nparticles為與壁面碰撞的顆粒數(shù);mp為顆粒的質(zhì)量流量，kg/s;α為顆粒與壁面的碰撞角度，rad;f(α)為關(guān)于顆粒碰撞角度的函數(shù);υ為顆粒相對(duì)于壁面的速度，m/s;b(υ)為關(guān)于此相對(duì)速度的函數(shù);A為與壁面材料有關(guān)的系數(shù);Fs為顆粒的形狀系數(shù)，尖銳顆粒取1，半圓形顆粒取0.53，圓形顆粒取0.2;Aface為壁面上計(jì)算單元的面積，m2。

3、網(wǎng)格劃分及邊界條件

　　氣體入口采用速度邊界條件，大小為20m/s，出口采用自由流邊界條件，壁面無滑移。固體顆粒由入口面釋放，且其進(jìn)入速度等于入口處的流體速度。顆粒相在壁面處的邊界類型為反彈(reflect)，在出口處的邊界類型則為逃逸(escape)。另外，顆粒質(zhì)量流量設(shè)置為16.286kg/s，此條件下固體顆粒體積濃度約為6%，符合Fluent中離散相模型的相關(guān)要求。

　　進(jìn)行網(wǎng)格劃分時(shí)，對(duì)流體域設(shè)置全局單元尺寸，并對(duì)球閥內(nèi)部通道處的流域進(jìn)行局部加密。而通過對(duì)網(wǎng)格及釋放顆粒數(shù)無關(guān)性的研究，最終確定將全局單元尺寸設(shè)置為6mm，加密處單元尺寸設(shè)置為1.5mm，而由入口釋放的顆粒數(shù)為15000。

4、結(jié)果討論與分析

　　4.1、顆粒運(yùn)動(dòng)規(guī)律及壁面磨損分布情況分析

　　顆粒對(duì)壁面的碰撞是造成其受到磨損的最主要原因，因此其在整個(gè)流體通道內(nèi)的運(yùn)動(dòng)軌跡是研究壁面磨損分布的主要依據(jù)。而由于在上、下游直管段大部分區(qū)域隨流運(yùn)動(dòng)的顆粒與壁面產(chǎn)生的碰撞很少，因此該區(qū)域壁面的磨損量可忽略不計(jì)�？梢酝ㄟ^對(duì)圖3中顆粒運(yùn)動(dòng)軌跡和圖4中垂直于Z軸的截面顆粒質(zhì)量濃度分布的分析，來研究圖2中磨損集中區(qū)的成因。以標(biāo)準(zhǔn)狀態(tài)下的球閥為例，當(dāng)氣固流流經(jīng)球閥時(shí)，球閥壁面會(huì)產(chǎn)生三處磨損集中區(qū)(圖2中的1、2、3)，其中，圖2中1區(qū)域處位于球閥入口處的閥球外壁，顆粒在此處發(fā)生碰撞后的軌跡受球形凸面影響較大，且絕大多數(shù)顆粒在碰撞過程中動(dòng)量損失較大，因此這些顆粒在1區(qū)域發(fā)生碰撞后會(huì)被流速較大的氣流直接帶入球閥內(nèi)部，或在上游直管段處與壁面發(fā)生二次碰撞后隨流進(jìn)入球閥內(nèi)部(圖3)，而由圖4(a)可知，球閥入口附近流域的右側(cè)因顆粒碰撞閥球外壁后反彈而濃度較大;2區(qū)域處位于球閥內(nèi)通道的一側(cè)壁面，該處產(chǎn)生磨損的原因主要是閥球內(nèi)通道與顆粒運(yùn)動(dòng)軌跡形成了一定角度，而無論與閥球外壁碰撞與否，大部分顆粒都在慣性力的作用下于球閥入口處與氣流分離，而由圖4(b)可知，在球閥入口處與氣流分離后的顆粒主要偏向一側(cè)運(yùn)動(dòng)，并在該區(qū)域壁面上發(fā)生碰撞;3區(qū)域處位于靠近球閥出口的下游直管段，該處壁面因受到從球閥出口處流出的顆粒的碰撞而產(chǎn)生磨損，這些顆粒同樣是在自身慣性的作用下于球閥出口處與氣流分離，并在隨后與下游段壁面發(fā)生碰撞的，而由圖4(c)也可知，在球閥出口附近的流域內(nèi)顆粒偏向一側(cè)聚集，并于下游直管段壁面發(fā)生碰撞。

圖2 球閥壁面磨損分布圖

圖3 固體顆粒典型運(yùn)動(dòng)軌跡

圖4 流道截面顆粒質(zhì)量濃度分布圖

　　4.2、球閥結(jié)構(gòu)參數(shù)對(duì)壁面磨損的影響規(guī)律

　　球閥的結(jié)構(gòu)參數(shù)包括開度和流道直徑兩方面。為探究在其它條件不變的情況下開度對(duì)于球閥磨損情況的影響，共設(shè)置了15%、25%、50%、75%共4種不同的球閥開度。如圖5、6所示，球閥開度變化對(duì)其壁面磨損的分布和大小均有影響。隨著球閥開度的減小，其壁面處的最大磨損速率和總體質(zhì)量損失均會(huì)相應(yīng)增大，且增大速率隨著開度的減小而加快;另外，當(dāng)開度逐漸減小時(shí)，閥球通道內(nèi)壁上的磨損集中區(qū)(圖2中2區(qū)域)會(huì)由塊狀逐漸轉(zhuǎn)至帶狀，且下游段磨損集中區(qū)(圖2中3區(qū)域)會(huì)逐步向下移動(dòng)。出現(xiàn)上述的變化是因?yàn)楫?dāng)球閥開度減小后，球閥出入口處的過流截面更加狹窄，氣流在該處的速度變化更為劇烈，且對(duì)固體顆粒的作用更加明顯，使其在分離后能以更大的速度沖擊閥球通道內(nèi)壁面，雖然球閥開度的變化會(huì)改變顆粒沖擊壁面時(shí)的攻角并進(jìn)而影響磨損量，但與氣流對(duì)顆粒的加速作用相比，攻角變化對(duì)于磨損量的影響要小得多。因此，隨著球閥開度的減小，壁面磨損速率的最大值和總體質(zhì)量損失均會(huì)增大，且磨損速率最大的區(qū)域會(huì)出現(xiàn)在閥球內(nèi)壁上;球閥開度的減小同時(shí)意味著閥球內(nèi)通道和直管段之間的夾角更大，這使得閥球通道內(nèi)壁沿氣流方向上可與固體顆粒發(fā)生碰撞的區(qū)域更大，故該處集中區(qū)會(huì)由塊狀逐漸轉(zhuǎn)化為帶狀，除此之外，夾角的增大也會(huì)使固體顆粒在閥球內(nèi)壁碰撞后以一個(gè)較大的角度發(fā)生反彈，而這會(huì)使得下一次與顆粒發(fā)生碰撞的區(qū)域更加靠近球閥出口處，故下游段磨損集中區(qū)會(huì)逐步向下方移動(dòng)。

圖5 不同開度下球閥壁面磨損分布圖

圖6 不同開度下球閥壁面磨損量變化曲線

　　流體通道直徑的變化也會(huì)對(duì)球閥壁面的磨損產(chǎn)生一定的影響。為探究在其它條件不變的情況下流道直徑對(duì)于球閥磨損情況的影響，共設(shè)置了65mm、80mm、100mm、125mm共四種不同的球閥流道直徑尺寸。如圖7、8所示，球閥壁面的最大磨損速率和總體質(zhì)量損失均隨流道直徑的增大而減小，但磨損的分布情況變化不大。這是因?yàn)楫?dāng)閥門開度保持不變，而流道直徑增大時(shí)，閥腔的內(nèi)部空間會(huì)變大，而流體在流經(jīng)閥球時(shí)，其速度劇烈增大的趨勢(shì)將會(huì)減緩，固體顆粒對(duì)壁面的沖擊強(qiáng)度也會(huì)隨之減小，故球閥壁面的最大磨損速率和總體質(zhì)量損失均隨流道直徑的增大而減小，但流道直徑的增加并未改變流體通道的總體構(gòu)造，故真空技術(shù)網(wǎng)(http://www.mp99x.cn/)認(rèn)為并不會(huì)對(duì)球閥壁面整體的磨損分布情況產(chǎn)生大的影響。

圖7 不同流體通道直徑下球閥壁面磨損分布圖

圖8 不同流體通道直徑下球閥壁面磨損量變化曲線

5、結(jié)論

　　1)當(dāng)氣固兩相流流經(jīng)球閥時(shí)，其壁面上共存在三處磨損集中區(qū)，第一處位于球閥入口處的閥球外壁上，此處顆粒與壁面碰撞后的軌跡受球形凸面影響較大;第二處位于通道內(nèi)低速區(qū)一側(cè)壁面上，這個(gè)區(qū)域內(nèi)的磨損主要是由固體顆粒在球閥入口處與氣流分離后沖擊閥球通道內(nèi)壁而產(chǎn)生的;第三處磨損集中區(qū)在球閥出口處附近，固體顆粒在閥球內(nèi)壁處發(fā)生碰撞并反彈，在自身慣性作用下于球閥出口處與氣流分離并沖擊壁面，最終造成該區(qū)域的磨損。

　　2)重點(diǎn)研究了球閥結(jié)構(gòu)參數(shù)(開度、流道直徑)的變化對(duì)球閥壁面磨損的影響，球閥開度及流道直徑的減小均會(huì)使壁面上最大磨損速率和總體質(zhì)量損失增大，但只有開度的變化才會(huì)對(duì)磨損分布產(chǎn)生較大的影響，且當(dāng)開度減小時(shí)，位于通道內(nèi)壁的集中區(qū)(圖2中2區(qū)域)由塊狀逐漸轉(zhuǎn)變?yōu)閹�狀，而球閥出口附近的集中區(qū)(圖2中3區(qū)域)則逐漸向下移動(dòng)。

　　3)以本次研究結(jié)果為基礎(chǔ)，后續(xù)將進(jìn)一步研究通過改進(jìn)球閥內(nèi)部流道結(jié)構(gòu)的方法以減輕其壁面的磨損，從而達(dá)到提高球閥的使用壽命的目的。