楊宗政 董迎春 李芳 吳志國 張忠園

1.天津科技大學(xué)海洋與環(huán)境學(xué)院 2.天津科技大學(xué)化工與材料學(xué)院 3.天津市塘沽鑫宇環(huán)?？萍加邢薰?/p>

鉆井平臺在石油開采過程中產(chǎn)生的含油污水大多含有泥-砂等固體顆粒[1],使得往返平臺與港口間的油輪在運輸過程中產(chǎn)生的壓艙水、清洗廢水也會攜帶泥-砂。通常,港口的含油污水由就近的污水廠接收并處理,但這些含油污水產(chǎn)量大且成分復(fù)雜,所含油相的質(zhì)量濃度在10 ～ 10 000 mg/L之間,工業(yè)上多采用隔油池、氣浮破乳、重力過濾和膜分離等來處理含油污水[2-3],但上述設(shè)施投資均較高且占地面積大。水力旋流器是一種依靠兩相間的密度差,在離心力的作用下實現(xiàn)固-液或液-液兩相分離的設(shè)備,因其具有尺寸小、占地少、分離快、效率高、操作簡便、易維護等優(yōu)點,在石油化工等領(lǐng)域得到了廣泛的應(yīng)用。

CFD數(shù)值模擬作為新興的研究技術(shù),相比傳統(tǒng)的實驗研究,具有研究成本低、研發(fā)周期短的優(yōu)點,能大幅度縮減實驗的次數(shù),在水力旋流器的設(shè)計和研發(fā)領(lǐng)域得到了廣泛的應(yīng)用[4-8]。通過計算流體動力學(xué)的方法來模擬分離器的油水分離過程,可以清楚地認識到旋流器內(nèi)流場的分布[9-12],掌握其分離過程的理論依據(jù),為油水旋流器的結(jié)構(gòu)尺寸及操作參數(shù)優(yōu)化提供方便。設(shè)計了一種用于港口油、水、砂混合液分離的組合式一體化設(shè)備,通過借助CFD模擬技術(shù)對液液二級旋流器進行結(jié)構(gòu)參數(shù)及操作參數(shù)的優(yōu)化[13-14],以期設(shè)計的一體化設(shè)備能夠?qū)崿F(xiàn)油相回收資源化利用,同時使水相中油含量減少,降低后續(xù)污水生化處理段的有機負荷,以便于污水廠處理。使用了含砂的油污水進行實例驗證,可為今后設(shè)備的模擬研發(fā)提供參考。

1 研究對象及方法

1.1 研究對象

設(shè)計了一種組合式一體化油、水、砂分離設(shè)備,由收水存儲箱、提升泵、固液初級旋流器、液液二級旋流器、儲渣桶、集油桶、中間池、集水箱等組件構(gòu)成,結(jié)構(gòu)如圖1所示。油、水、砂混合液首先進入收水存儲箱,然后經(jīng)提升泵提升后進入旋流器完成三相分離,分離后的固相進入儲渣桶,液相進入集水箱,最終油相經(jīng)中間池進一步脫水收集。本研究的主要對象為一體化設(shè)備中的液液二級旋流器,其結(jié)構(gòu)如圖2所示。

1.2 研究方法

利用Fluent軟件對液液二級旋流器進行仿真模擬,分析內(nèi)部流場情況?；谡辉囼炘O(shè)計,對液液油水旋流器的結(jié)構(gòu)參數(shù)(旋流腔直徑、溢流口直徑、底流口直徑、大錐角)進行優(yōu)化,再使用優(yōu)化的旋流器進行最佳操作條件(入口流速、溢流分流比)探索。最后對模擬設(shè)計的旋流器進行油水分離能力測試,驗證模擬的結(jié)果是否準確、設(shè)計的設(shè)備尺寸是否可行。

2 旋流器的模擬研究

2.1 模型建立及網(wǎng)格劃分

表1 旋流器各部分尺寸表結(jié)構(gòu)參數(shù)尺寸結(jié)構(gòu)參數(shù)尺寸D70.00 mmLc70.00 mmDc35.00 mmLu700.00 mmDo4.00 mmα20.0°Di12.25 mmβ1.5°Du17.50 mm

液液水力旋流器內(nèi)部的流場屬于多相流模型[15-16],而且其中的流場運動屬于較為復(fù)雜的湍流流動狀態(tài),雷諾應(yīng)力模型考慮了流體運動彎曲、漩渦、旋轉(zhuǎn)等問題,但該模型對網(wǎng)格質(zhì)量要求很高。為保證計算時結(jié)果更穩(wěn)定、精確,采用六面體網(wǎng)格進行劃分,并對入口與圓柱段接口等區(qū)域采取了局部網(wǎng)格加密處理,最終總網(wǎng)格數(shù)達到572 748個,總節(jié)點數(shù)達到568 350個。其網(wǎng)格劃分示意圖如圖3所示。

2.2 求解設(shè)置及結(jié)果分析

求解器選取基于壓力的SIMPLEC 壓力-速度模型,湍流模型選擇適合模擬高速旋轉(zhuǎn)流場的雷諾應(yīng)力模型,離散格式選用二階迎風(fēng)式插值格式[16-19]。邊界條件設(shè)為:混合液為速度入口,上下出口均為自由流出口,旋流器壁面條件則設(shè)置為靜止、無移動壁面。

物性參數(shù)選取液態(tài)水為連續(xù)相,分散相為油,其油相密度為845 kg/m3,黏度為0.002 4 Pa·s,水相密度為998 kg/m3,黏度為0.001 Pa·s。

在設(shè)定入口流速為15 m/s,入口油相體積分數(shù)(Ci,下同)為1%的混合液,溢流分流比為5%的工況下進行模擬(見圖4、圖5)。如圖4油相旋轉(zhuǎn)速度矢量圖及流線圖所示,混合液沿著入口切向進入旋液分離器后,在初始速度及重力的影響下,油相在經(jīng)過旋流腔和大錐段加速后開始反旋向上溢流,但并未完全溢流出,而是停留在圓柱段內(nèi),說明該工況下從底流口涌出的混合液大部分為水相,這與圖5所示的油相體積分數(shù)分布云圖結(jié)果一致。

在設(shè)定入口流速為15 m/s的情況下,溢流分流比為5%～30%,在入口混合液中Ci從1%逐步增加到10%的過程中,其分離效率的變化如圖6所示。當(dāng)入口流速為15 m/s時,在分流比從5%增加到30%的過程中,油相體積分數(shù)為4%、7%、10%的油水分離實驗的分離效率均依次經(jīng)過急劇增加、平穩(wěn)增加、快速降低的過程,最佳分流比分別為12%、15%、20%,但分離效率均不高。因為溢流口徑較小,匯聚在溢流口附近的油滴還未排出便受到進液的不斷沖擊而分散,并隨混合液旋轉(zhuǎn)遷移至底流口排出,因此,分離效率不理想。

2.3 結(jié)構(gòu)參數(shù)優(yōu)化

液液二級旋流器的結(jié)構(gòu)參數(shù)包括入口、溢流口、底流口、旋流腔、大錐段、小錐段、平尾段等尺寸[20-21]。在不同分流比和不同混合液油相體積分數(shù)的條件下,不同尺寸的溢流口對分離效率有一定的影響。底流口直徑尺寸對分離效率的影響主要是通過影響液體在旋流器內(nèi)的停駐時間而實現(xiàn)。旋流腔的直徑尺寸決定了旋流器的處理量,其大小與分離效率相關(guān)。旋流器內(nèi)混合液的油水分離效率受離心力影響,離心力的大小與旋轉(zhuǎn)動量有關(guān),而旋轉(zhuǎn)動量由大錐角的大小決定。

因此,選取4個主要因素作為變量,分別為溢流口直徑Do、底流口直徑Du、旋流腔直徑D及大錐角α。對每個因素分別取3個水平設(shè)計正交實驗,選擇L9(34)正交表進行9次關(guān)鍵實驗,對液液二級旋流器的結(jié)構(gòu)參數(shù)進行優(yōu)化,正交試驗因素水平如表2所列。除4個主要因素外,其余的結(jié)構(gòu)尺寸不變,分別為:Dc=35 mm、Di=7 mm、Lc=70 mm、L2=89 mm、L3=535 mm、Lu=700 mm。

三、我國獸醫(yī)管理體制改革的建議機構(gòu)不健全、基層動物防疫體系薄弱，缺乏專業(yè)人才，法律不完善、等問題，嚴重影響了動物疫病防治能力和動物產(chǎn)品質(zhì)量安全水平的提高。加之目前國際上特別是我周邊國家重大動物疫病時有發(fā)生，甚至出現(xiàn)局部蔓延，對我國防控重大動物疫病形成較大壓力。因此，改革獸醫(yī)管理體制已經(jīng)成為當(dāng)前十分緊迫的任務(wù)。

表2 正交試驗因素水平表水平號因素列D/mmDo/mmDu/mmα/(°)166.5 (1.9Dc)4.20 (0.12Dc)14.0 (0.4Dc)16270.0 (2.0Dc)4.55 (0.13Dc)17.5 (0.5Dc)18373.5 (2.1Dc)4.90 (0.14Dc)21.0 (0.6Dc)20

根據(jù)正交試驗設(shè)計的條件進行數(shù)值模擬,利用Fluent軟件依次進行液液二級旋流器的幾何模型構(gòu)建和網(wǎng)格劃分(網(wǎng)格數(shù)量在57萬個左右),設(shè)定求解器、計算模型、邊界條件及物性參數(shù)等條件后,初始化計算域后求解計算。設(shè)定Y1～Y9型的旋流器均在入口流速為15 m/s、混合液中Ci為10%、溢流分流比為10%的條件下進行模擬實驗。

表3 正交實驗表試驗號因素列D/mmDo/mmDu/mmα/(°)分離效率/%Y166.54.20 14.0 1678.53Y266.54.5517.51879.49Y366.54.9021.02081.22Y470.04.2017.52078.31Y570.04.5521.01679.73Y670.04.9014.01880.40Y773.54.2021.01874.76Y873.54.5514.02079.82Y973.54.9017.51677.27K1239.24231.60238.75235.53K2238.44239.04235.07234.65K3231.85238.89235.71239.35K179.7577.2079.5878.51K279.4879.6878.3678.22K377.2879.6378.5779.78R2.462.481.231.57主次順序Do>D>α>Du最優(yōu)組合Do=4.6 mm,D=66.5 mm,α=20°,Du=14 mm

為了保證優(yōu)化模擬試驗結(jié)果不受網(wǎng)格數(shù)量的影響,對具有最大分離效率的Yb型旋流器進行網(wǎng)格無關(guān)性驗證[22]。在保證網(wǎng)格質(zhì)量的條件下,劃分出6種不同精密度的網(wǎng)格(見表4)。網(wǎng)格的無關(guān)性驗證結(jié)果如圖7所示,網(wǎng)格數(shù)量為level-3和level-4時能夠保證計算結(jié)果的準確性且計算時間較短。因此,本優(yōu)化實驗使用網(wǎng)格數(shù)量在31萬～39萬個之間的模型。

表4 網(wǎng)格劃分水平表個網(wǎng)格劃分水平網(wǎng)格數(shù)量網(wǎng)格劃分水平網(wǎng)格數(shù)量Level-1238 531Level-4381 488Level-2273 095Level-5438 674Level-3319 670Level-6497 573

優(yōu)化后的Yb型旋流器內(nèi)部油相旋轉(zhuǎn)速度流線圖、壓力圖與油相體積分數(shù)分布云圖如圖8所示,從圖8(a)中優(yōu)化后的油相運動軌跡和圖8(b)中相體積分布云圖可以看出,油水混合液從雙入口處流入旋流腔后先在圓柱段產(chǎn)生旋流,油滴在離心力和徑向壓力的作用下與水相逐漸分離向著軸線處運動并沿著軸線方向形成油核。隨著旋流場的逐漸穩(wěn)定,大部分油滴從旋流器中心位置的排油管排出。從圖8(c)中壓力分布云圖可以看出,油滴受到的徑向壓力沿著半徑方向逐漸減小。

在Yb型旋流器的A=60 mm、A=160 mm、A=210 mm處分別做截面,3種截面位置如圖9所示。

以3種截面位置的圓心為原點,其直徑上的切向、徑向、軸向速度分布如圖10所示。從圖10(a)切向速度分布圖可看出,旋流器內(nèi)的切向速度呈軸對稱分布,使流體產(chǎn)生離心力而向壁面運動,旋轉(zhuǎn)半徑沿軸心向壁面方向增加時,切向速度先增大后減小。在最大切向速度點,混合液所受離心力最大,最易發(fā)生油水分離,經(jīng)過該點后切向速度急劇減小,該位置因接近壁面而存在運動阻力。

從圖10(b)徑向速度分布圖可看出,徑向速度分布呈現(xiàn)中心對稱分布,兩側(cè)的徑向速度方向均為壁面指向中心線,且速度是逐漸減小的,說明混合液沿著壁面流動時,混合液中的油向著中心線匯集,混合液正在發(fā)生油水兩相分離過程。

從圖10(c)軸向速度分布圖可看出,旋流器內(nèi)流體的軸向速度變化趨勢相同,靠近旋流器內(nèi)壁面水相富集的環(huán)狀區(qū)域內(nèi),軸向速度方向指向底流口。流體所受的軸向速度隨半徑的減小而減小,直到在某點處變?yōu)?,成為零軸向速度點;半徑再繼續(xù)減小時,軸向速度由負變?yōu)檎曳较虬l(fā)生了改變,指向溢流口形成逆向流。此后隨著半徑的減小,軸向速度則逐步增大,同時,油相不斷向著中心軸線附近移動而聚集,最終從溢流口排出。

2.4 最佳操作參數(shù)探究

在入口流速為5 m/s、7 m/s、9 m/s、12 m/s、14 m/s、17 m/s、22 m/s,油相體積分數(shù)為10%,溢流分流比為10%的條件下進行模擬。入口流速與分離效率、壓力降的關(guān)系如圖11所示,當(dāng)入口流速從5 m/s增至12 m/s時,分離效率從21.71%快速增至79.75%。這是因為當(dāng)流速過低時,旋流器內(nèi)尚未形成逆向流,中心軸線附近未出現(xiàn)穩(wěn)定的油核,但隨著流速的增大,旋流場內(nèi)的逆向流趨于穩(wěn)定而形成油核,油相開始溢流而出。入口流速為12 ～17 m/s時,由于已經(jīng)形成了穩(wěn)定的油核,增大流速對分離效率提升不大,當(dāng)流速繼續(xù)增大時,由于進液的持續(xù)沖擊使聚集的油滴重新分散,從而導(dǎo)致分離效率有所下降,而壓降隨著入口流速的增大而增大。故Yb型旋流器分離的最佳速度參數(shù)為15 m/s。

在最佳分離速度15 m/s的條件下,分流比與分離效率、壓力降的關(guān)系如圖12所示,當(dāng)溢流分流比從0增至15%時,增大分流比能提高分離效率。當(dāng)分流比增加到15%時,分離效率開始下降,原因是分流比的增大使得混合液被迫從上出口溢流而出,且旋流場內(nèi)存在沖擊壓力,湍流程度的增加導(dǎo)致聚集的油滴粒子破裂分散成小粒徑的油滴,從而影響了分離效率。當(dāng)分流比大于30%時,會加劇未完全分離的混合液從溢流口流出,溢流口水相增多導(dǎo)致分離效率開始急劇下降。壓降不斷增大是由于分流比的增大使溢流口處的壓力增大而導(dǎo)致的。

3 實例驗證

通過CFD數(shù)值模擬得到的最佳結(jié)構(gòu)尺寸及操作參數(shù)并不代表實際工況下的分離效率,CFD模擬是通過數(shù)值方法求解一系列的微分方程組對實際物理模型在理想條件下的簡化模型進行預(yù)測,需通過實例驗證[7]。

實驗樣品取自某燃料油庫區(qū)含油污水,對液液二級旋流器的油水分離能力進行了模擬實驗驗證,混合液中的Ci為10%,經(jīng)旋流分離后取樣檢測,檢測結(jié)果見表5。

表5 水樣檢測結(jié)果序號入口流速/(m·s-1)溢流分流比/%底流口油相體積分數(shù)/%油相去除率/%15107.9021.0027106.0839.1639103.8961.09412102.7372.75514102.4375.70617102.4675.40722102.8471.6081556.5634.40915103.9860.211015152.1978.061115202.0979.091215252.1378.751315302.8271.801415406.8531.501515507.2427.62

圖13為油相分離效率模擬值與實驗值比較情況,其模擬值與實驗值大致相同,實際的分離器最佳操作參數(shù)為入口流速15 m/s,溢流分流比15%～20%,實際除油效率達到79%。經(jīng)旋流器分離后,水相中油相體積分數(shù)大幅度降低,便于后續(xù)的污水生化處理,表明優(yōu)化的設(shè)備能夠滿足油水分離預(yù)處理的要求,具有良好的應(yīng)用前景,可為今后設(shè)備的模擬研發(fā)提供參考。

4 結(jié)論

(1) 通過正交試驗對液液二級旋流器進行了結(jié)構(gòu)參數(shù)優(yōu)化,在最優(yōu)尺寸下,分離油相體積分數(shù)為10%的混合液時,模擬的最佳操作參數(shù)為入口流速15 m/s、溢流分流比在15%～20%之間,除油效率為83.10%,液液二級旋流器對油水混合液的預(yù)分離有較好的應(yīng)用前景。

(2) 入口流速一定時,適當(dāng)增加溢流分流比能夠提高油水分離效率,超過最佳分流比會導(dǎo)致未完全分離的混合液從溢流口流出,分離效率開始急劇下降。

(3) 在一定的入口流速范圍內(nèi),分離效率會隨著入口流速的增大而增大。超過該范圍時,分離效率開始下降,因為入口混合液的高速旋流沖擊會導(dǎo)致聚合的油滴破裂而分散。

(4) 對模擬研發(fā)的組合式設(shè)備進行了油水分離能力實例驗證,其最佳入口流速為15 m/s,溢流分流比在15%～20%之間,除油效率為79%,與模擬值基本一致,說明該模擬可為實際應(yīng)用提供可靠的參考。