鄧福成，黃斌，尹彪，李小森，易先中

(1.長(zhǎng)江大學(xué)機(jī)械工程學(xué)院，湖北荊州，434000；2.中國(guó)科學(xué)院廣州能源研究所，廣東廣州，510640；3.中國(guó)科學(xué)院廣州天然氣水合物研究中心，廣東廣州，510640；4.中國(guó)科學(xué)院天然氣水合物重點(diǎn)實(shí)驗(yàn)室，廣東廣州，510640)

天然氣水合物的全球儲(chǔ)量非常巨大，是全球常規(guī)燃料總碳量的2 倍[1]，天然氣水合物的資源量為1.0×1017～1.0×1018m3[2]，被視為可替代石油的重要能源。目前已有多個(gè)國(guó)家對(duì)天然氣水合物的開采進(jìn)行了多次試采實(shí)驗(yàn)，但天然氣水合物儲(chǔ)層出砂嚴(yán)重是其大規(guī)模商業(yè)化開采中的瓶頸問題。水合物主要分布于凍土區(qū)和海洋深水沉積物中[3]，試采中普遍使用降壓法、注熱法和氣體置換法[4]等開采方法，以降壓法最為常見，但降壓開采時(shí)儲(chǔ)層物性變化會(huì)造成大面積出砂[5-7]。開采過程中砂巖脫落的砂粒隨著流體進(jìn)入井筒，對(duì)井下設(shè)備造成破壞，從而制約天然氣水合物安全高效開采。為減小天然氣儲(chǔ)層出砂對(duì)開采的影響，人們采用機(jī)械防砂篩管來解決油氣井開發(fā)中天然氣儲(chǔ)層的出砂問題。沖縫管套作為防砂結(jié)構(gòu)中的重要元件，一般用于篩管保護(hù)套或獨(dú)立防砂篩管。當(dāng)出現(xiàn)沖蝕磨損時(shí)，管套防砂性能降低，影響油井生產(chǎn)。因此，分析沖縫管套的沖蝕壽命成為篩管研究的重要部分。

目前，國(guó)內(nèi)學(xué)者通過建立沖蝕模型，研究篩管沖蝕磨損問題，該方法已趨于成熟，并取得了豐富的研究成果。張銳等[8]建立離散顆粒篩管沖蝕模型，運(yùn)用氣固耦合模型模擬兩相流流場(chǎng)分布和沖蝕特征，為氣井防砂的沖蝕計(jì)算提供了研究思路；翟曉鵬等[9]用離散顆粒流數(shù)值模擬方法分析不同流速和體積分?jǐn)?shù)下砂粒對(duì)金屬網(wǎng)布速度分布和沖蝕速率影響，預(yù)測(cè)金屬網(wǎng)布篩管沖蝕壽命；張君濤等[10]在防砂篩管模擬中引入沖蝕模型，分析不同壓差下割縫篩管的沖蝕特征，得到?jīng)_蝕的主要集中區(qū)域，對(duì)于長(zhǎng)效防砂和提高篩管使用壽命有著重要作用；陳珊珊等[11]基于CFD 的固液兩相流沖蝕模擬分析，預(yù)測(cè)篩管壽命，并利用室內(nèi)沖蝕試驗(yàn)驗(yàn)證數(shù)值模擬預(yù)測(cè)沖蝕壽命的可行性。但這些研究都只考慮了部分篩網(wǎng)結(jié)構(gòu)下兩相流的沖蝕模擬，未見針對(duì)水氣砂三相流模擬和沖縫管套結(jié)構(gòu)等沖蝕磨損問題的研究。

首先，本文作者建立沖縫管套篩縫流道模型，基于多相流模型(Mixture)和離散相模型(DPM)來模擬不同因素下砂粒流動(dòng)對(duì)管套沖蝕的影響；其次，分析水合物分解后，氣、液和固三相流下流體流速、顆粒體積分?jǐn)?shù)、開口高度及液相體積分?jǐn)?shù)對(duì)沖縫管套沖蝕速率的影響；最后，推導(dǎo)沖縫管套壽命預(yù)測(cè)模型，為天然氣水合物儲(chǔ)層試采及商業(yè)化開采的長(zhǎng)效安全性提供理論支持。

1 水合物儲(chǔ)層物性特征

國(guó)內(nèi)學(xué)者根據(jù)神狐海域地質(zhì)情況建立了單層水合物層和水合物-游離氣層2種成藏模型[12-13]，其中第2種模型更符合現(xiàn)場(chǎng)實(shí)際，即水合物層存在于上覆層和下伏層中，且下伏層與水合物層間存在含有大量離散甲烷氣體的游離氣層。天然氣水合物常以固態(tài)形式存在，開采時(shí)主要通過降壓、注熱等方法使其分解，再將分解后的天然氣運(yùn)送至地面。水合物分解過程如下：

降壓過程中儲(chǔ)層物質(zhì)存在著相變分解與再形成，導(dǎo)致沉積物中包含氣、液和固三相物質(zhì)[14]，但水合物一般作為膠結(jié)物賦存在于沉積物中。水合物的分解導(dǎo)致砂巖應(yīng)力變化，從而致使砂巖顆粒脫落。脫落的顆粒在甲烷氣和水共同滲流的作用下發(fā)生遷移，部分未完全分解的水合物微粒也會(huì)與地層中的砂粒一起流入井筒。機(jī)械篩管是一種水合物常用的防砂手段，在生產(chǎn)過程中，甲烷氣和水將滲流進(jìn)篩管保護(hù)套，且顆粒將隨流體一起進(jìn)入保護(hù)套，造成篩管的沖蝕。

2 計(jì)算模型

水合物儲(chǔ)層生產(chǎn)過程中的沖蝕過程是氣、液和固三相流作用下的沖蝕問題。根據(jù)儲(chǔ)層中各物質(zhì)特性，計(jì)算中將氣體和液體視為連續(xù)相，顆粒為離散相。因此，水合物儲(chǔ)層中流體攜帶顆粒沖蝕計(jì)算主要分為流體控制方程、顆粒運(yùn)動(dòng)方程及沖蝕模型[15]。

2.1 氣、液流體控制方程

控制方程包括連續(xù)性方程、動(dòng)量方程和能量方程，流體分為氣體與液體，雙流體的連續(xù)性方程為[16]：

式中：α為體積分?jǐn)?shù)，%；ρ為密度，kg/m3；t為時(shí)間，s；?為拉普拉斯算子；u為速度矢量，m/s；下標(biāo)g和l分別表示氣相和液相；αl+αg= 1。

動(dòng)量方程為[17]：

式中：pl和pg分別為液相和氣相壓力，Pa；g為重力加速度，通常取9.8 m/s2；和分別為液相和氣相的應(yīng)力張量，Pa；M為相間動(dòng)量交換系數(shù)，kg/(m·s)2；μl，μs和ξ分別為液體、氣體和流體的剪切黏度，Pa·s；I為單位向量。

能量方程為[18]：

式中：σl和σg分別為液相和氣相剪應(yīng)力張量；el和eg分別為液相和氣相的比熱力學(xué)能，J/kg；K為導(dǎo)熱系數(shù)，W/(m·K)；T為熱力學(xué)溫度，K；Q為內(nèi)部熱源。

在水合物儲(chǔ)層開采中，水合物分解產(chǎn)生的氣體與液體混合，從而形成水氣兩相流。模擬中采用混合流模型，忽略水氣間的作用力。

2.2 固相顆粒運(yùn)動(dòng)方程

顆粒隨流體沖蝕沖縫管套過程中，流體(氣體和液體)視為連續(xù)相，顆粒為離散相。根據(jù)牛頓第二定律，顆粒的受力方程為[19]

式中：G為顆粒所受重力；FD，F(xiàn)Δp，F(xiàn)F，F(xiàn)V和FL分別為顆粒所受的拖曳力、壓力梯度力、升力、附加質(zhì)量力以及浮力。在顆粒運(yùn)動(dòng)中還存在許多力，如Besset 力、Magnus 力和Saffman 力等。Besset 力是瞬時(shí)的流動(dòng)阻力，與流體的穩(wěn)定性有關(guān)，當(dāng)顆粒加速度較小時(shí)，可忽略該力。Magnus和Saffman力等都只在顆粒高速旋轉(zhuǎn)的特定情況下才有意義，故受力分析中可忽略。

顆粒在流體中所受重力G為

式中：d為顆粒粒徑，m；ρs為顆粒密度，kg/m3；Vs為顆粒體積，m3。

壓力梯度力FΔp為

式中：p為壓力，Pa。

流體內(nèi)部壓強(qiáng)差形成的浮力FL為

顆粒表面流體流速差導(dǎo)致的升力FF為

式中：CL為升力系數(shù)；v為流體作用在顆粒上的速度，m/s；As為顆粒平均投影面積，m2。

流體相對(duì)于顆粒運(yùn)動(dòng)，流體產(chǎn)生的曳力(氣動(dòng)阻力)FD為

式中：CD為阻力系數(shù)；νf為顆粒速度，m/s。

顆粒加速所需要的高于顆粒本身的附加作用力FV。

顆粒在沖蝕過程中不斷撞擊壁面，顆粒對(duì)壁面的作用力是顆粒與壁面的碰撞力和摩擦力，顆粒間碰撞力忽略不計(jì)。根據(jù)GRANT 等[20]的研究，得到顆粒沖擊-反彈系數(shù)為：

式中：eN為切向反彈系數(shù)；eT為垂向反彈系數(shù)；θ為顆粒入射角度，rad。

2.3 顆粒沖蝕模型

在多相流模型計(jì)算中，歐拉-歐拉法將不同的相視為連續(xù)介質(zhì)；而歐拉-拉格朗日法則將流體視為連續(xù)相，顆粒視為離散相。本文將顆粒視為離散相，運(yùn)用動(dòng)能、質(zhì)量和能量傳遞作為橋梁來計(jì)算模型。在氣、液和固三相模型計(jì)算中，需要用歐拉-拉格朗日法進(jìn)行計(jì)算，離散相沖蝕模型計(jì)算公式[21]為

式中：Rerosion為壁面磨損速率，kg/(m2·s)；N為碰撞顆粒數(shù)；C(dp)為顆粒直徑函數(shù)，C(dp)=1.8×10-9；f(ω)為沖蝕角度函數(shù)，如式(19)所示；ω為顆粒軌跡與壁面的沖擊角度，rad；νsl為顆粒的相對(duì)速度，m/s；b(νsl)為顆粒相對(duì)速度函數(shù)，取2.6；Aface為沖蝕面積，m2。

3 沖縫管套沖蝕模型

3.1 沖縫管套篩縫模型

在天然氣水合物的開采防砂過程中，沖縫管套用于保護(hù)篩管或者獨(dú)立防砂，其與地層流體直接接觸，易受到?jīng)_蝕磨損。沖縫管套結(jié)構(gòu)如圖1所示，根據(jù)其流道結(jié)構(gòu)簡(jiǎn)化出單位流道模型結(jié)構(gòu)，建立沖縫管套篩縫的單位流體域模型。入口面為矩形截面，長(zhǎng)為22 mm，寬為3 mm；壁面高度為11 mm，出口高度根據(jù)沖縫管套的設(shè)計(jì)條件設(shè)置為不同尺寸。圖2所示為單位篩縫流道模型及網(wǎng)格模型，圖2中藍(lán)色面網(wǎng)格(網(wǎng)格模型上平面)為流場(chǎng)入口面，綠色面網(wǎng)格為篩縫壁面，黃色面網(wǎng)格為流場(chǎng)出口面(開口面)。

圖1 沖縫管套示意圖Fig.1 Schematic diagram of punched screen

圖2 單位篩縫流道模型及網(wǎng)格模型圖Fig.2 Flow channel model and grid model of unit screen

3.2 模型條件設(shè)置

離散相采用DPM 模型，對(duì)壁面沖蝕進(jìn)行氣、液和固三相三維湍流穩(wěn)態(tài)計(jì)算。計(jì)算中湍流模型采用RNGk-ε模型，運(yùn)用DPM 模型對(duì)顆粒運(yùn)動(dòng)進(jìn)行追蹤，固體球形顆粒由入口面注入。流體為氣相和液相混合，壓力和速度耦合用Coupled 算法。流體模型為多相流模型，液相與氣相的入口速度相同。當(dāng)液相體積分?jǐn)?shù)較小時(shí)，液相均勻分散在氣相中，忽略水氣間的相間作用力，將流道區(qū)域內(nèi)流體擬定為氣液混合流體(等效為在氣液混合均勻的單相流體)，兩者共享同一壓力場(chǎng)和速度場(chǎng)，且忽略相對(duì)滑移速度。

氣相為甲烷，考慮低溫、高壓對(duì)氣體的影響，密度為1.225 kg/m3；液相為水，密度為998.2 kg/m3；固相為砂粒顆粒，密度為1 500 kg/m3，沉積物中大部分顆粒是粉砂，其顆粒粒徑為4～63 μm[22]，取顆粒平均粒徑為40 μm，粒徑正態(tài)分布。由于分析的流道模型是管套內(nèi)實(shí)際流道模型的基本單元，并沿軸線軸向分布，且分析模型尺寸較小，故不考慮重力影響。入口處邊界條件設(shè)置速度入口、流體和顆粒的速度及方向。流動(dòng)中不涉及熱力學(xué)，即不考慮流動(dòng)過程中水合物二次生成等現(xiàn)象，流動(dòng)過程視為等溫絕熱。管壁條件采用無滑移壁面，出口邊界設(shè)置outflow。

2017年我國(guó)在南海神狐海域第一輪試采中累計(jì)產(chǎn)氣量超過3.0×104m3，日平均產(chǎn)量為5.0×103m3，最高日產(chǎn)量為3.5×104m3；第二次試采創(chuàng)造了產(chǎn)氣總量為8.614×105m3，日平均產(chǎn)量為2.87×104m3的世界紀(jì)錄[23]。假定該井的水平段長(zhǎng)度為100 m，按日平均產(chǎn)量2.87×104m3可計(jì)算得到流體流速為0.125 m/s，計(jì)算時(shí)初始流體流速取0.100 m/s。根據(jù)申志聰?shù)萚24]建立的模型，水合物儲(chǔ)層下方存在游離氣層。若水平井處于水合物層中，開采初始階段壓差使得水合物分解和游離氣層逸散出的混合氣體流入管道。此時(shí)，氣體流量將遠(yuǎn)大于初始流量?？紤]到各因素對(duì)管道破壞的影響，合理選取不同因素的取值，如表1所示。

表1 模擬因素水平表Table 1 Simulation factor level table

在沖蝕計(jì)算中，由于流體流速及組成比例不同，顆粒數(shù)量根據(jù)體積來確定更適宜，顆粒體積分?jǐn)?shù)(含砂體積分?jǐn)?shù))按在入口面處所占面積比進(jìn)行計(jì)算，因流體流速發(fā)生變化，單位時(shí)間內(nèi)進(jìn)入篩縫流道的顆?？傮w積隨之改變。計(jì)算得到不同顆粒體積分?jǐn)?shù)和不同流體流速下的顆?？傎|(zhì)量流量：

式中：Qm為顆?？傎|(zhì)量流量，kg/s；φs為顆粒體積分?jǐn)?shù)，%；vs為顆粒在入口處速度，m/s；Aintel為入口面積，m2。

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

4.1 流場(chǎng)特征

分析計(jì)算在入口流速為0.1 m/s、顆粒體積分?jǐn)?shù)為5%、液相體積分?jǐn)?shù)為0.25%和開口高度為0.30 mm 條件下顆粒對(duì)篩縫結(jié)構(gòu)的沖蝕磨損情況。圖3所示為上述條件下沖蝕速率分布云圖。從圖3可見，顆粒進(jìn)入流道對(duì)底部壁面和斜端壁面造成沖蝕磨損，且底部壁面靠近開口處的沖蝕磨損速率最大。

圖3 模型沖蝕速率分布云圖Fig.3 Cloud map of model erosion rate distribution

沖蝕過程中顆粒會(huì)在底部與壁面發(fā)生碰撞[25]后經(jīng)出口面流出，造成篩縫開口處的沖蝕磨損。圖4所示為顆粒軌跡圖。從圖4可看出，部分顆粒撞擊壁面后，顆粒隨流體向出口運(yùn)動(dòng)。分析顆粒運(yùn)動(dòng)軌跡，發(fā)現(xiàn)顆粒在篩縫壁面和底部中線區(qū)域的停留時(shí)間長(zhǎng)、速度小、軌跡復(fù)雜。

圖4 顆粒軌跡圖Fig.4 particle trajectory

4.2 沖蝕速率的影響因素

4.2.1 流體速度

由式(20)可知：當(dāng)流速改變時(shí)，進(jìn)入流道內(nèi)的顆粒質(zhì)量流量增大，單位時(shí)間內(nèi)撞擊壁面的顆粒數(shù)量增加、顆粒撞擊頻率升高，沖蝕速率增大。根據(jù)表1中數(shù)據(jù)，分析計(jì)算得到?jīng)_縫管套沖蝕速率。圖5所示為不同顆粒體積分?jǐn)?shù)下沖蝕速率隨流速變化折線圖，對(duì)數(shù)據(jù)進(jìn)行擬合，沖蝕速率y和流速x的擬合關(guān)系為

式中：ɑ和b均為冪函數(shù)參數(shù)。圖5中各組擬合線與相應(yīng)數(shù)據(jù)的重合程度較好。篩縫上沖蝕速率隨流速增加呈現(xiàn)指數(shù)增長(zhǎng)趨勢(shì)[26]。當(dāng)流速增加時(shí)，顆粒獲得更大動(dòng)能，對(duì)壁面的沖擊磨損程度增大，沖蝕速率增大。

圖5 不同顆粒體積分?jǐn)?shù)下沖蝕速率隨流速變化折線圖Fig.5 Line chart of erosion rate change with velocity at different sand volume fractions

4.2.2 顆粒體積分?jǐn)?shù)

圖6所示為不同流速下沖蝕速率隨顆粒體積分?jǐn)?shù)變化。從圖6可得：沖蝕速率隨顆粒體積分?jǐn)?shù)增加，呈現(xiàn)線性增長(zhǎng)趨勢(shì)。由式(20)可知，隨顆粒體積分?jǐn)?shù)增大，顆粒質(zhì)量流量增加。顆粒數(shù)量增多，顆粒與壁面碰撞頻率上升，造成更嚴(yán)重的沖蝕磨損。SUN 等[13]研究發(fā)現(xiàn)：當(dāng)顆粒質(zhì)量流量達(dá)到一定值后，沖蝕速率隨質(zhì)量流量增加出現(xiàn)增長(zhǎng)緩慢的趨勢(shì)。顆粒增多導(dǎo)致顆粒間相互碰撞更加劇烈，動(dòng)能減小，沖蝕速率降低。在本次分析中，未考慮顆粒間碰撞的因素。

4.2.3 液相體積分?jǐn)?shù)

圖7所示為不同液相體積分?jǐn)?shù)下沖蝕速率隨流速變化折線圖。由圖7可知：沖蝕速率隨著液相體積分?jǐn)?shù)增加而增大，呈現(xiàn)出線性增長(zhǎng)的趨勢(shì)。由于流道開口面積遠(yuǎn)小于入口面積，相同流量流體通過開口面的流速將遠(yuǎn)大于入口面的速度。流體中液相黏性大，影響顆粒運(yùn)動(dòng)，增加顆粒動(dòng)能；而氣相黏性小，顆粒速度的增加不明顯[27]。因此，液相體積分?jǐn)?shù)增加時(shí)，出口處顆粒動(dòng)能增加，沖蝕速率增大。在計(jì)算中，未考慮重力的影響，故氣液兩相的密度對(duì)沖蝕速率影響不大，而兩者黏度的差異造成氣相更易從流道中逸出，殘余的液相對(duì)顆粒運(yùn)動(dòng)產(chǎn)生影響。

圖7 不同液相分?jǐn)?shù)下沖蝕速率隨流速變化折線圖Fig.7 Line chart of erosion rate change with flow velocity at different liquid phase volume fractions

在流速為0.2 m/s 時(shí)，沖蝕速率隨液相體積分?jǐn)?shù)增加出現(xiàn)先減小后增大的現(xiàn)象。由于模型計(jì)算中2種流體為均勻混合，但兩者黏度的差異導(dǎo)致混合流體在流經(jīng)篩縫底部時(shí)，液相更容易附著于底部。氣體逸出篩縫后，在底部沖蝕區(qū)域液體隨著體積分?jǐn)?shù)增加更易形成液膜。液相體積分?jǐn)?shù)增加對(duì)沖蝕速率有著增大和抑制2個(gè)方面作用：

1)由于液相對(duì)顆粒動(dòng)能變化起主導(dǎo)作用，隨著液相體積分?jǐn)?shù)增加，沖蝕速率增大；

2)當(dāng)液相體積分?jǐn)?shù)增加時(shí)，由于液相更易附著于壁面，部分沖蝕區(qū)域?qū)⑿纬梢耗?，減小顆粒動(dòng)能，從而減弱顆粒對(duì)管套壁面的沖蝕作用。

4.2.4 出口面高度

圖8所示為不同開口高度下沖蝕速率隨流速變化折線圖。由圖8可見：當(dāng)流速一定時(shí)，開口高度減小，壁面沖蝕速率增大。開口高度減小為0.10 mm時(shí)，隨流體流速增加，沖蝕速率的增長(zhǎng)速率增加得更為明顯。當(dāng)開口高度變小時(shí)，出口面積減小使得出口處流速增加、動(dòng)能增大，碰撞壁面的顆粒數(shù)量增多。顆粒動(dòng)能和碰撞數(shù)量增加造成沖蝕區(qū)域的沖蝕程度加劇。

圖8 不同開口高度下沖蝕速率隨流速變化折線圖Fig.8 Line diagram of erosion rate change with flow velocity at different opening heights

開口高度分別為0.10 mm和0.15 mm時(shí)，堵塞現(xiàn)象導(dǎo)致沖蝕速率增長(zhǎng)幅度明顯。分析計(jì)算中顆粒粒徑為0.02～0.06 mm，當(dāng)粒徑大于開口的1/3時(shí)，易形成砂橋，發(fā)生堵塞現(xiàn)象[28]。開口高度再次減小或顆粒數(shù)量增加時(shí)，該現(xiàn)象發(fā)生的概率將增大。在現(xiàn)場(chǎng)應(yīng)用中，若部分篩管發(fā)生堵塞，由于生產(chǎn)井產(chǎn)量一定，將增加未堵塞篩管的過流速度，增大沖蝕速率。

4.3 壽命預(yù)測(cè)分析

4.3.1 預(yù)測(cè)模型計(jì)算

由圖3可知，最大沖蝕速率出現(xiàn)在沖縫管套篩縫底部邊緣區(qū)域。分析計(jì)算入口流速為0.10 m/s、顆粒體積分?jǐn)?shù)為5%、液相體積分?jǐn)?shù)為1%和開口高度0.30 mm條件下篩管結(jié)構(gòu)的沖蝕速率，得出篩縫底部的沖蝕速率云圖如圖9所示。從圖9可知，底部沖蝕區(qū)域面積占比約為2/3。

圖9 篩縫底部沖蝕云圖Fig 9 Cloud map of bottom of screen seam

擬用最大沖蝕速率ER表示底部區(qū)域的平均沖蝕速率E(kg·m-2·s-1)，但由于沖蝕速率會(huì)隨時(shí)間而變化，故需引入修正系數(shù)k。從圖9可看出：部分沖蝕區(qū)域沖蝕速率與最大沖蝕速率相差較大，沖蝕區(qū)域的沖蝕速率的最小值與最大值的比值為0.1，綜合考慮后取k=0.1。

沖蝕區(qū)域的平均沖蝕速率E為

顆粒沖蝕篩管將造成篩縫開口高度增大，導(dǎo)致篩管過濾精度發(fā)生改變，防砂失效。沖蝕深度按壁面材料的質(zhì)量損失進(jìn)行計(jì)算。該沖蝕模型的沖蝕主要發(fā)生在開口底部一側(cè)，計(jì)算沖縫管套的質(zhì)量損失MR為

式中：S為篩縫處沖蝕的面積，m2。

根據(jù)體積公式，將質(zhì)量損失轉(zhuǎn)化為體積損失后，可計(jì)算沖蝕深度為

式中：L為沖蝕深度，mm；VR為體積損失，m3；ρb為管套材料的密度，kg/m3。

沖縫管套預(yù)測(cè)使用壽命T為

式中：D為篩縫允許沖蝕深度，mm。

管套常用材料304 不銹鋼密度為7 930 kg/m3。根據(jù)GILLESPIE 等[27]進(jìn)行的沖蝕磨損試驗(yàn)，縫寬增加0.05 mm時(shí)，出砂量會(huì)嚴(yán)重增多?？紤]到?jīng)_蝕發(fā)生在底部，故D=0.05 mm。

4.3.2 壽命預(yù)測(cè)結(jié)果

根據(jù)式(25)，當(dāng)允許沖蝕深度為0.05 mm 時(shí)，預(yù)測(cè)使用壽命與沖蝕速率呈反比例關(guān)系，分析計(jì)算后可得到管套壽命預(yù)測(cè)數(shù)據(jù)。在天然氣產(chǎn)出的初始階段，流速大于初始條件值。按最大日產(chǎn)量3.5×104m3計(jì)算，流速取0.2 m/s。將此條件下的數(shù)據(jù)繪制為沖縫管套預(yù)測(cè)壽命折線圖如圖10所示。

圖10 不同因素下沖縫管套預(yù)測(cè)壽命折線圖Fig.10 Line diagram of predicted life of punched screen for different factors

從圖10可以得出：沖縫管套的預(yù)測(cè)使用壽命隨流體速度增加而減小，在速度增加的初始段，預(yù)測(cè)使用壽命下降趨勢(shì)明顯，該情況與實(shí)際符合。預(yù)測(cè)使用壽命隨著開口高度增加而增大，在開口高度為0.15～0.30 mm 間減少趨勢(shì)為線性關(guān)系。開口高度為0.15 mm和0.10 mm時(shí)有所偏差的原因是堵塞現(xiàn)象，導(dǎo)致沖蝕速率驟增、預(yù)測(cè)壽命減小幅度大。預(yù)測(cè)使用壽命隨流體中液體體積分?jǐn)?shù)和顆粒體積分?jǐn)?shù)上升而減小。在液相體積分?jǐn)?shù)為0.5%時(shí)，由于底部區(qū)域殘余的液相較少，無法形成液膜，導(dǎo)致沖蝕速率大、預(yù)測(cè)壽命小。沖蝕速率與顆粒體積分?jǐn)?shù)之間的關(guān)系為線性關(guān)系，由式(25)和圖10(c)可知，預(yù)測(cè)使用壽命與顆粒體積分?jǐn)?shù)之間的關(guān)系為反比關(guān)系。

5 結(jié)論

1)利用CFD 方法分析了沖縫管套篩縫的流場(chǎng)分布情況，管套篩縫受到的沖蝕磨損主要集中在底部邊緣區(qū)域。

2)在水合物儲(chǔ)層開采中，流體流速、液相體積分?jǐn)?shù)增加和開口高度減小都是通過直接或間接方式增大顆粒動(dòng)能，進(jìn)而增大沖蝕速率，而顆粒體積分?jǐn)?shù)增加通過增加顆粒碰撞壁面的數(shù)量來增大沖蝕速率。

3)流體中液相體積分?jǐn)?shù)增加對(duì)沖蝕有著增大和抑制2 種作用。當(dāng)液相體積分?jǐn)?shù)增大時(shí)，一方面，會(huì)增加顆粒動(dòng)能，增大沖蝕速率；另一方面，在沖蝕區(qū)域形成液膜，削弱顆粒動(dòng)能、減少摩擦磨損。

4)沖縫管套預(yù)測(cè)使用壽命在流速增大的初始階段下降較為明顯。當(dāng)液相體積分?jǐn)?shù)增大時(shí)，形成的液膜將起到一定的減小沖蝕、延長(zhǎng)管套使用壽命的作用。