張曉光 趙立新 徐保蕊 蔣明虎 鄧 鑫

(1. 東北石油大學(xué)機(jī)械科學(xué)與工程學(xué)院；2. 中國(guó)石油天然氣管道工程有限公司)

錐角變化對(duì)旋流除砂性能影響的數(shù)值模擬*

張曉光**1趙立新1徐保蕊1蔣明虎1鄧 鑫2

(1. 東北石油大學(xué)機(jī)械科學(xué)與工程學(xué)院；2. 中國(guó)石油天然氣管道工程有限公司)

針對(duì)在油田開(kāi)展采出液除砂的必要性進(jìn)行了介紹，對(duì)所研究的除砂型旋流器結(jié)構(gòu)和基本原理進(jìn)行了說(shuō)明。利用數(shù)值模擬軟件Fluent，將某油田現(xiàn)場(chǎng)采出液作為分離對(duì)象，對(duì)主直徑為56mm的除砂旋流器進(jìn)行變錐角結(jié)構(gòu)數(shù)值模擬研究。得出不同錐角結(jié)構(gòu)旋流器的速度場(chǎng)、壓力降和砂相體積分?jǐn)?shù)分布的變化規(guī)律，并對(duì)不同錐角結(jié)構(gòu)旋流器的除砂分離效率進(jìn)行了對(duì)比。研究表明：當(dāng)錐角為5°時(shí)，旋流器內(nèi)部產(chǎn)生的切向速度最大，底流出口砂相體積分?jǐn)?shù)分布最高，但也伴隨著較大的壓力降，所完成的砂相分離效率最高為96.30%。

旋流器 錐角 除砂 分離性能 數(shù)值模擬

隨著油井開(kāi)采時(shí)間的延長(zhǎng)，油井采出液中含砂量越來(lái)越高，甚至一些油井從投產(chǎn)時(shí)就存在嚴(yán)重的采出液含砂問(wèn)題[1]?，F(xiàn)有的沉降清砂方法周期長(zhǎng)、占地面積大、所耗費(fèi)成本較高，已不能適應(yīng)油井采出液產(chǎn)液量逐年上升、含砂量逐年增加的現(xiàn)狀。近年來(lái)不進(jìn)罐清砂技術(shù)(如旋流分離)在石油工業(yè)中應(yīng)用廣泛。除砂用水力旋流器是利用固液兩相介質(zhì)間的密度差將分散相——砂相從連續(xù)相——液相中分離出來(lái)的一種技術(shù)方法，具有裝置緊湊、占地面積小、分離周期短、處理量調(diào)節(jié)范圍大及設(shè)備投資成本小等優(yōu)點(diǎn)。

除砂用水力旋流器是一種應(yīng)用廣泛的分離設(shè)備，不僅可用于油田除砂、脫泥，還可以用于礦場(chǎng)分級(jí)分離、洗滌等。研究表明旋流器錐角結(jié)構(gòu)參數(shù)對(duì)旋流器的分離性能影響較大，但由于旋流器處理介質(zhì)不同、主要結(jié)構(gòu)形式不同，錐角變化對(duì)旋流器分離性能的影響規(guī)律也不盡相同[2，3]。越來(lái)越多的研究表明，利用數(shù)值模擬方法研究旋流器內(nèi)分離流場(chǎng)的方法具有成本低、速度快的優(yōu)點(diǎn)，模擬結(jié)果確實(shí)可反映旋流器內(nèi)流體的流動(dòng)規(guī)律[4，5]。

油井采出液除砂對(duì)于旋流器的要求主要是提高旋流器底流口排砂濃度，筆者結(jié)合油田采出液含砂的實(shí)際情況，借助于計(jì)算流體動(dòng)力學(xué)CFD數(shù)值模擬方法，研究分析了一種單錐形水力旋流器的錐角變化對(duì)旋流器速度場(chǎng)、壓力降、底流含砂濃度和除砂分離效率的影響規(guī)律，給出了適用于油井采出液除砂旋流器的最佳錐角。

1 旋流分離器結(jié)構(gòu)分析

旋流器的工作原理是運(yùn)用離心作用進(jìn)行液固兩相分離，水力旋流器本身無(wú)任何運(yùn)動(dòng)部件，由待分離的固液兩相混合液以一定的壓力從水力旋流器入口切向進(jìn)入，產(chǎn)生強(qiáng)烈旋轉(zhuǎn)運(yùn)動(dòng)，由于固液兩相的密度差，兩相所受的離心力、向心浮力和流體曳力的大小不同，大部分的重質(zhì)相(固相)經(jīng)由旋流器底流出口排出，而輕質(zhì)相(液相)則由溢流出口排出，從而達(dá)到固液兩相分離的目的[6]。

旋流器通常由入口、旋流腔、錐段、溢流出口和底流出口組成(圖1)[7]，主要結(jié)構(gòu)參數(shù)包括主直徑D1、錐角α、入口當(dāng)量直徑di、旋流腔長(zhǎng)度L1、溢流管直徑Do、溢流管伸入長(zhǎng)度Lo及底流管直徑Du等，其中，主直徑與錐角兩個(gè)參數(shù)最為重要，這是因?yàn)樾黧w長(zhǎng)度由D1與α決定，其他結(jié)構(gòu)參數(shù)也均與D1成一定比例關(guān)系。入口形式選擇切向雙入口矩形結(jié)構(gòu)，其當(dāng)量水力直徑計(jì)算值為主直徑的0.22倍，矩形入口斷面與旋流腔壁面相切，能較好地消除來(lái)液短路的死區(qū)，減弱可能出現(xiàn)的循環(huán)流，使入口來(lái)液液流更穩(wěn)定。其他主要參數(shù)設(shè)置分別為：Lo=0.714D1、Do=0.357D1、L1=D1。

圖1 旋流器結(jié)構(gòu)示意圖

圖2為不同錐角結(jié)構(gòu)旋流器的結(jié)構(gòu)形式對(duì)比圖，變換錐角主要對(duì)旋流器錐段的長(zhǎng)度產(chǎn)生影響，而旋流器旋流腔、溢流管段、底流管段的結(jié)構(gòu)參數(shù)不變。模擬結(jié)構(gòu)確定旋流器主直徑為D1=56mm，用Hy-α5、Hy-α8和Hy-α10分別代表旋流器錐角為5、8、10°的結(jié)構(gòu)，其錐段計(jì)算長(zhǎng)度分別為5.71D1、3.58D1和2.86D1。

圖2 不同錐角水力旋流器結(jié)構(gòu)形式對(duì)比

2 參數(shù)設(shè)置與網(wǎng)格劃分

2.1介質(zhì)物性參數(shù)與初始條件設(shè)置

進(jìn)行模擬計(jì)算時(shí)，對(duì)應(yīng)介質(zhì)的物性參數(shù)為：油相密度870kg/m3，動(dòng)力粘度0.046 1kg/(m·s)，油相粒徑0.09mm；水的密度998.2kg/m3，動(dòng)力粘度1.003g/(m·s)；砂的密度2 500kg/m3，動(dòng)力粘度1.72×10-5kg/(m·s)，砂平均粒徑0.09mm。

根據(jù)油田某站實(shí)際情況，初始模擬計(jì)算邊界條件設(shè)置為：入口流量4m3/h，入口油相體積分?jǐn)?shù)為7.9%，入口含砂1.5%，溢流分流比為80%；入口法向速度分量為9.26m/s，其他兩個(gè)方向的速度分量為零；溢流出口和底流出口均設(shè)定為自由出口。

2.2網(wǎng)格劃分

利用CFD前處理軟件Gambit進(jìn)行模型構(gòu)建和網(wǎng)格劃分，劃分網(wǎng)格如圖3所示，所劃分網(wǎng)格基本采用六面體結(jié)構(gòu)性網(wǎng)格形式，不同錐角結(jié)構(gòu)形式旋流器的入口段、旋流腔段、溢流出口管段和底流出口管段網(wǎng)格劃分一致，對(duì)于錐段部分依照等間距劃分網(wǎng)格方法，使得不同錐角旋流器的錐段部分網(wǎng)格劃分密度一致。同時(shí)進(jìn)行網(wǎng)格獨(dú)立性檢驗(yàn)[8]，分別劃分不同數(shù)量的網(wǎng)格水平，對(duì)比不同網(wǎng)格劃分水平對(duì)底流出口含砂體積分?jǐn)?shù)變化的影響，結(jié)果顯示網(wǎng)格數(shù)為281 200時(shí)，底流出口砂相體積分?jǐn)?shù)分布基本不隨網(wǎng)格數(shù)的增多而變化。

圖3 旋流器模型建立和網(wǎng)格劃分

3 速度場(chǎng)影響分析

主要對(duì)旋流器內(nèi)速度場(chǎng)、壓力降和固相體積分?jǐn)?shù)分布進(jìn)行分析。

旋流器內(nèi)液流為三維湍流流動(dòng)，其速度場(chǎng)用圓柱坐標(biāo)系來(lái)表示，分為3個(gè)分速度：切向速度vt、軸向速度va和徑向速度vr。由于切向速度在數(shù)值上遠(yuǎn)大于軸向速度和徑向速度分量，因此分析中暫忽略軸向和徑向兩個(gè)速度分量，主要針對(duì)切向速度進(jìn)行重點(diǎn)研究[9～12]。

vt是旋流器內(nèi)速度場(chǎng)中最重要的一個(gè)速度分量，vt的大小決定著旋流器中產(chǎn)生的離心加速度和離心力的大小，也是固液分離中固相顆粒分離的先決條件。由于旋流器為軸對(duì)稱結(jié)構(gòu)，且入口采用對(duì)稱雙入口結(jié)構(gòu)形式，因此選取旋流腔與錐段連接處橫截面I作為不同錐角結(jié)構(gòu)旋流器內(nèi)部速度場(chǎng)的對(duì)比截面，該截面上的速度變化直接受到旋流器錐段變化的影響，更能體現(xiàn)旋流器錐角變化對(duì)分離器內(nèi)部流體流動(dòng)的影響規(guī)律。

對(duì)比截面I上，不同錐角結(jié)構(gòu)旋流器的vt分布對(duì)比曲線圖如圖4所示，由圖可見(jiàn)，截面I上切向速度呈周向?qū)ΨQ分布，3種錐角結(jié)構(gòu)旋流器的vt變化趨勢(shì)基本一致：中心處速度均為零；隨半徑的增加，切向速度也逐漸增加，且存在一個(gè)最大切向速度點(diǎn)vtmax，在該點(diǎn)處隨著半徑的進(jìn)一步增加，旋流器vt逐漸減小，在邊壁處的切向速度減小為零，因此以最大切向速度為界，使得旋流器內(nèi)部分為內(nèi)渦流區(qū)和外渦流區(qū)兩個(gè)區(qū)域，由于溢流分流比設(shè)置為80%，明顯高于底流分流比，可以看出內(nèi)渦流區(qū)也明顯大于外渦流區(qū)。

圖4 切向速度分布對(duì)比

由圖4也可看出3種錐角結(jié)構(gòu)旋流器的切向速度變化規(guī)律的不同：在徑向位置(0～3mm)范圍內(nèi)，Hy-α5旋流器變化梯度最大，隨錐角的增加，該區(qū)域的切向速度變化梯度逐漸減小，且Hy-α8與Hy-α10兩種錐角結(jié)構(gòu)旋流器的切向速度梯度變化差異不大；徑向位置(3～20mm)范圍內(nèi)，Hy-α5旋流器的切向速度在一定范圍內(nèi)基本不變，然后繼續(xù)增大到最大值，Hy-α8與Hy-α10兩種旋流器切向速度變化呈線性增加趨勢(shì)也逐漸增大到最大切向速度，隨著錐角的增加，旋流器內(nèi)部所能達(dá)到的最大切向速度點(diǎn)值逐漸減小，Hy-α5旋流器的最大切向速度值最高。因此，其他設(shè)置參數(shù)一致情況下，旋流器錐角增加對(duì)內(nèi)外渦流的分界面影響很小，但會(huì)改變其旋轉(zhuǎn)強(qiáng)度，Hy-α5結(jié)構(gòu)旋流器內(nèi)部流場(chǎng)的旋轉(zhuǎn)強(qiáng)度最大。

4 壓力降影響分析

水力旋流器的壓力降為進(jìn)口處壓力與溢流出口管(或底流出口)處的壓力之差。圖5為不同錐角結(jié)構(gòu)旋流器內(nèi)部壓力降的分布對(duì)比圖，對(duì)比截面為旋流器的縱向截面z=0，由圖可見(jiàn)在截面z=0上，旋流器內(nèi)部壓力降由旋流器器壁處向中心逐漸增大，在旋流器中心處達(dá)到最大值，沿旋流器器壁向底流出口方向壓力降呈逐漸增大趨勢(shì)，在旋流腔內(nèi)部，旋流器中心與溢流出口底部相交位置處出現(xiàn)最大壓力降值，且由該位置沿旋流器中心軸向向上到溢流出口和向下到底流出口截面，壓力降呈現(xiàn)出逐漸降低趨勢(shì)，即3種不同錐角旋流器內(nèi)部的最大壓力降均發(fā)生在溢流出口伸入部分底部。3種不同錐角旋流器所達(dá)到的最大壓力降不同，錐角越小壓力降越大，因Hy-α5旋流器錐段最長(zhǎng)，底流壓力降與溢流出口相比其他錐角結(jié)構(gòu)旋流器壓力降則更高；Hy-α8和Hy-α10旋流器，由于錐段減小，底流出口壓力降明顯低于溢流出口壓力降。因此當(dāng)旋流器其他設(shè)置參數(shù)不變時(shí)，僅錐角變化時(shí)對(duì)旋流器內(nèi)部壓力降的影響較為明顯。

圖5 不同錐角旋流器內(nèi)部z=0截面上的壓力降分布對(duì)比

由圖5可見(jiàn)旋流器內(nèi)部主要壓力降分布包括溢流出口壓力降Δpo，底流出口壓力降Δpu和溢流伸入管處最大壓力降Δpmax，對(duì)比不同錐角結(jié)構(gòu)旋流器的主要壓力降分布如圖6所示。由圖可見(jiàn)，隨錐角增加壓力降呈逐漸降低趨勢(shì)，錐角越小壓力降越大，且溢流壓力降比底流壓力降增加幅度也最大，錐角由8°變化到10°時(shí)二者壓力降變化較小；底流壓力降最低，因?yàn)樽罡邏毫捣植荚谝缌鞴苌烊氩糠值撞?，與溢流出口較近，溢流壓力降比最高壓力降值降低幅度不大。

5 砂相體積分?jǐn)?shù)分布影響分析

砂相經(jīng)旋流器離心分離后，主要分布在旋流器器壁臨近區(qū)域，因此主要對(duì)比不同錐角結(jié)構(gòu)旋流器器壁臨近區(qū)域內(nèi)的砂相體積分?jǐn)?shù)Fs分布。在相同操作參數(shù)下，不同錐角旋流器固相體積分?jǐn)?shù)梯度分布如圖7所示，對(duì)比可見(jiàn)，砂相經(jīng)旋流腔和錐段分離后基本沿旋流器內(nèi)壁流向底流出口，整體上底流含砂濃度最高。Hy-α5、Hy-α8和Hy-α10旋流器所達(dá)到的底流砂相體積分?jǐn)?shù)最高值分別為29.9%、20.0%和17.5%。

圖7 不同錐角旋流器砂相體積分?jǐn)?shù)梯度分布

為進(jìn)一步對(duì)比底流出口含砂濃度受錐角結(jié)構(gòu)變化的影響，針對(duì)底流出口截面上的砂相體積分?jǐn)?shù)分布曲線進(jìn)行對(duì)比，如圖8所示。由圖可見(jiàn)，底流出口截面不同徑向位置處的砂相體積分?jǐn)?shù)是變化的，基本呈軸向?qū)ΨQ分布，且由邊壁到底流出口中心處隨半徑減小而逐漸降低，在臨近邊壁7mm范圍內(nèi)，也是旋流器底流出口砂相含量最集中的部分，含砂體積分?jǐn)?shù)隨錐角的增加而逐漸降低，Hy-α5結(jié)構(gòu)旋流器的底流出口含砂體積分?jǐn)?shù)分布最高；在中心處和臨近部分，含砂相體積分?jǐn)?shù)很小，小于0.1%，其分布對(duì)旋流器底流出口含砂濃度影響較小。整體上Hy-α5結(jié)構(gòu)旋流器的底流出口含砂濃度最高。

圖8 不同錐角旋流器底流出口截面上砂相體積分?jǐn)?shù)分布曲線對(duì)比

6 固液分離效果影響分析

除砂分離效率Es為底流出口固相砂的質(zhì)量流率與入口來(lái)液中固相砂質(zhì)量流率的比：

(1)

式中Mi——入口砂的質(zhì)量流率，kg/s；

Mu——底流出口砂的質(zhì)量流率，kg/s。

計(jì)算得出不同錐角結(jié)構(gòu)旋流器除砂分離效率變化如圖9所示。由圖可見(jiàn)，旋流器錐角變化對(duì)油田采出液砂相分離效率影響較大，在研究范圍內(nèi)，錐角越小旋流器除砂分離效率越高，當(dāng)錐角為5°時(shí)，旋流器除砂分離效率最高可達(dá)96.30%。

7 結(jié)論

7.1不同錐角結(jié)構(gòu)旋流器的切向速度變化趨勢(shì)基本一致，截面I上切向速度呈周向?qū)ΨQ分布，以最大切向速度為界，旋流器內(nèi)部速度場(chǎng)被分為內(nèi)渦流區(qū)和外渦流區(qū)現(xiàn)象明顯，隨著錐角的增加，旋流器內(nèi)部所能達(dá)到的最大切向速度點(diǎn)值逐漸減小，旋流器錐角變化對(duì)速度場(chǎng)內(nèi)外渦流的分界面徑向位置影響較小，但對(duì)其旋轉(zhuǎn)強(qiáng)度影響較大，錐角為5°旋流器內(nèi)部流場(chǎng)的旋轉(zhuǎn)強(qiáng)度最大。

7.2旋流器內(nèi)溢流管伸入管底部所產(chǎn)生的壓力降值最大，溢流出口壓力降次之，底流出口的壓力降最低，錐角越小壓力降越大。

7.35°錐角結(jié)構(gòu)旋流器的底流出口含砂體積分?jǐn)?shù)最大，且砂相主要分布在旋流器器壁臨近區(qū)域，該結(jié)構(gòu)計(jì)算得出的砂相分離效率也最高。

7.4綜上可見(jiàn)，在一定范圍內(nèi)減小旋流器錐角可以提高除砂效率，但其壓力降也會(huì)增大，實(shí)際當(dāng)中應(yīng)綜合考量除砂分離效率與生產(chǎn)能耗等進(jìn)行參數(shù)選擇。

[1] 呂瑞典，李君裕，王遠(yuǎn)明，等.油井產(chǎn)液除砂旋流器試驗(yàn)研究[J].石油機(jī)械，1995，23(6)：18～23.

[2] 張丹，陳曄.錐角對(duì)固-液水力旋流器流場(chǎng)及其分離性能的影響[J].流體機(jī)械，2009，37(8)：11～16.

[3] 高淑玲，魏德洲，韓聰，等.底流口直徑和錐角對(duì)旋流器流場(chǎng)的影響[J].東北大學(xué)學(xué)報(bào)(自然科學(xué)版)，2010，31(5)：728～732.

[4] 盧濤.礦物顆粒在旋流復(fù)合力場(chǎng)中運(yùn)動(dòng)行為的數(shù)值模擬[D].沈陽(yáng)：東北大學(xué)，2011.

[5] 蔣明虎，劉道友，趙立新，等.錐角對(duì)水力旋流器壓力場(chǎng)和速度場(chǎng)的影響[J].化工機(jī)械，2011，38(5)：572～576.

[6] 趙立新，李楓.離心分離技術(shù)[M].哈爾濱：東北林業(yè)大學(xué)出版社，2006.

[7] Zhao L X，Jiang M H，Xu B R，et al.Development of a New Type High-efficient Inner-cone Hydrocyclone[J].Chemical Engineering Research and Design，2012，90(12)：2129～2134.

[8] 趙立新，韓麗艷，鄭國(guó)興，等.同向出流倒錐式旋流器的結(jié)構(gòu)模擬分析[J].流體機(jī)械，2013，41(2)：19～24.

[9] 徐保蕊.T形管油水分離器流場(chǎng)分布特性及分離性能研究[D].大慶：東北石油大學(xué)，2012.

[10] 趙立新，王尊策，李楓，等.液液水力旋流器流場(chǎng)特性與分離特性研究(一)——錐角變化對(duì)切向速度場(chǎng)的影響[J].化工裝備技術(shù)，1999，20(4)：1～3.

[11] 趙立新，李楓，王尊策，等.液液水力旋流器流場(chǎng)特性與分離特性研究(二)——錐角變化對(duì)軸向速度場(chǎng)的影響[J].化工裝備技術(shù)，1999，20(5)：1～3.

[12] 趙立新，王尊策，蔣明虎.液液水力旋流器流場(chǎng)特性與分離特性研究(三)——水力旋流器徑向速度測(cè)試方法[J].化工裝備技術(shù)，1999，20(5)：4～6.

NumericalSimulationofConeAngleEffectonDe-sandingPerformanceofHydrocyclone

ZHANG Xiao-guang1, ZHAO Li-xin1,XU Bao-rui1, JIANG Ming-hu1, DENG Xin2

(1.CollegeofMechanicalScienceandEngineering,NortheastPetroleumUniversity,Daqing163318,China;2.ChinaPetroleumPipelineEngineeringCorporation,Langfang065000,China)

The necessities of de-sanding produced liquid in oilfields were introduced, including both structure and fundamental principle of the cyclone-desander. Through making use of Fluent software and taking the produced fluid from an oilfield as a separation object, the research on numerically-simulating variable cone angles’ structure of the 56mm-diameter-cyclone-desander was implemented to obtain the distribution and the variation law of the velocity, pressure drop and sand phase volume fraction. Comparing the sand separation efficiency of the hydrocyclone with different cone angles shows that, the tangential velocity (vt) in the cyclone is maximum when the cone angle is 5° and the sand volume fraction distribution in the bottom outlet is the highest along with a larger pressure drop and a sand separation efficiency up to 96.30%.

hydrocyclone, cone angle, desanding, separation performance, numerical simulation

*國(guó)家“863”計(jì)劃課題(2012AA061303)，國(guó)家教育部高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金課題(20132322110002)。

**張曉光，女，1984年1月生，講師。黑龍江省大慶市，163318。

TQ051.8

A

0254-6094(2016)06-0798-05

2015-12-25，

2016-10-28)