高效厭氧生物反應(yīng)器內(nèi)的湍流特性及結(jié)構(gòu)優(yōu)化

王　宇， 王江云， 許雙雙， 陳春茂， 馬靜園， 郭紹輝

(1.中國(guó)石油大學(xué) 北京市油氣污染防治重點(diǎn)實(shí)驗(yàn)室 重質(zhì)油國(guó)家重點(diǎn)實(shí)驗(yàn)室， 北京 102249；2.中國(guó)海油能源發(fā)展股份有限公司 北京安全環(huán)保工程技術(shù)研究院， 北京 300457； 3.中國(guó)昆侖工程公司， 北京 100037)

高效厭氧生物反應(yīng)器內(nèi)的湍流特性及結(jié)構(gòu)優(yōu)化

王宇1，2， 王江云1， 許雙雙3， 陳春茂1， 馬靜園3， 郭紹輝1

(1.中國(guó)石油大學(xué) 北京市油氣污染防治重點(diǎn)實(shí)驗(yàn)室 重質(zhì)油國(guó)家重點(diǎn)實(shí)驗(yàn)室， 北京 102249；2.中國(guó)海油能源發(fā)展股份有限公司 北京安全環(huán)保工程技術(shù)研究院， 北京 300457； 3.中國(guó)昆侖工程公司， 北京 100037)

摘要：采用數(shù)值模擬方法研究了高效厭氧生物反應(yīng)器內(nèi)的湍流及多相流傳質(zhì)過(guò)程。通過(guò)重整化群k-ε湍流模型(RNG k-ε)和歐拉多相流模型構(gòu)建反應(yīng)器內(nèi)多相湍流流動(dòng)控制方程。在此基礎(chǔ)上，應(yīng)用自定義函數(shù)(UDF)方式修正了計(jì)算流體力學(xué)軟件Fluent中的曳力模型，并植入歐拉多相流模型，建立高效厭氧生物反應(yīng)器內(nèi)氣-液-固多相流動(dòng)模型，并模擬計(jì)算了反應(yīng)器內(nèi)多相流過(guò)程，同時(shí)對(duì)其結(jié)構(gòu)進(jìn)行優(yōu)化。結(jié)果表明，修正后的多相流動(dòng)模型計(jì)算得到的流場(chǎng)特征與實(shí)驗(yàn)數(shù)據(jù)吻合，驗(yàn)證了所采用的湍流模型和計(jì)算方法的準(zhǔn)確性；優(yōu)化后的反應(yīng)器內(nèi)多相流傳質(zhì)效果明顯提高。該方法較為準(zhǔn)確地預(yù)測(cè)了反應(yīng)器內(nèi)傳質(zhì)過(guò)程，可以用于高效厭氧生物反應(yīng)器的工程應(yīng)用改進(jìn)。

關(guān)鍵詞：高效厭氧生物反應(yīng)器； 數(shù)值模擬； 多相流動(dòng)模型； 結(jié)構(gòu)優(yōu)化

隨著石油工業(yè)發(fā)展，石油開(kāi)采加工過(guò)程產(chǎn)生的污水量大幅提升，水質(zhì)日益惡化，石油化工領(lǐng)域污水處理技術(shù)升級(jí)改造已成為工業(yè)污水處理亟待解決的問(wèn)題[1-2]。高效厭氧生物處理技術(shù)兼具污水處理和能源回收、處理成本低等優(yōu)勢(shì)，在煉化污水生化處理領(lǐng)域應(yīng)用前景廣闊[3-4]。我國(guó)石油化工領(lǐng)域的污水處理技術(shù)尚處于第1代反應(yīng)器向第2代反應(yīng)器過(guò)渡階段，高效厭氧生物反應(yīng)器在石化領(lǐng)域尚處于探索階段，在石油煉制領(lǐng)域尚無(wú)成功應(yīng)用案例，因此研究高效厭氧生物反應(yīng)器內(nèi)的流態(tài)分布、水力學(xué)和傳質(zhì)特性等過(guò)程動(dòng)力學(xué)對(duì)其在石油煉化領(lǐng)域的應(yīng)用具有重要意義。當(dāng)前，國(guó)內(nèi)對(duì)高效厭氧生物反應(yīng)器的研究多集中在反應(yīng)器的工藝、三相分離器結(jié)構(gòu)、污泥馴化等方面[5-7]，對(duì)基于反應(yīng)器內(nèi)部多相流動(dòng)過(guò)程數(shù)值模擬計(jì)算的結(jié)構(gòu)優(yōu)化研究較少，通常是通過(guò)實(shí)驗(yàn)方法對(duì)反應(yīng)器進(jìn)行改進(jìn)。由于厭氧生化反應(yīng)器的啟動(dòng)時(shí)間長(zhǎng)、操作控制參數(shù)較多、影響因素復(fù)雜等原因，根據(jù)實(shí)驗(yàn)數(shù)據(jù)優(yōu)化和研究反應(yīng)器內(nèi)部結(jié)構(gòu)和多相流特點(diǎn)費(fèi)時(shí)、費(fèi)力，采用數(shù)值模擬方法對(duì)反應(yīng)器多相流動(dòng)過(guò)程進(jìn)行研究分析具有明顯優(yōu)勢(shì)[8-10]。

Fluent軟件采用目前較先進(jìn)的離散技術(shù)和計(jì)算精度控制技術(shù)，通過(guò)流場(chǎng)的數(shù)值模擬不僅可得到實(shí)驗(yàn)手段無(wú)法得到的特定信息，同時(shí)節(jié)省了科研時(shí)間和經(jīng)費(fèi)[11-13]。但是，高效厭氧生物反應(yīng)器內(nèi)部流場(chǎng)復(fù)雜，對(duì)反應(yīng)器進(jìn)行模擬計(jì)算非常困難，多數(shù)研究停留在單向流的模擬計(jì)算階段，并未真實(shí)有效地模擬反應(yīng)器內(nèi)實(shí)際流場(chǎng)特點(diǎn)。筆者基于Fluent模擬軟件，采用修正的液-固曳力模型建立了能精確反映高效厭氧生物反應(yīng)器內(nèi)液、固流動(dòng)過(guò)程的數(shù)學(xué)模型，并進(jìn)行了模擬計(jì)算，在此基礎(chǔ)上對(duì)反應(yīng)器結(jié)構(gòu)進(jìn)行優(yōu)化改造，以期為高效厭氧生物反應(yīng)器的結(jié)構(gòu)優(yōu)化提供理論基礎(chǔ)。

1高效厭氧生物反應(yīng)器(ABR)的幾何模型及網(wǎng)格劃分

高效厭氧生物反應(yīng)器(以下簡(jiǎn)稱ABR)下部直徑70 mm，高度550 mm，上部三相分離高度250 mm。采用Gambit建模軟件對(duì)反應(yīng)器進(jìn)行了完全結(jié)構(gòu)化網(wǎng)格劃分，節(jié)點(diǎn)數(shù)106004個(gè)，并使用Fluent計(jì)算軟件對(duì)其內(nèi)部流場(chǎng)進(jìn)行數(shù)值模擬計(jì)算，具體的反應(yīng)器結(jié)構(gòu)尺寸及網(wǎng)格劃分如圖1所示。

圖1　ABR的結(jié)構(gòu)尺寸及網(wǎng)格劃分

2模擬ABR內(nèi)傳質(zhì)采用的數(shù)學(xué)模型

2.1多相湍流流動(dòng)控制方程

在數(shù)值模擬計(jì)算中，不考慮污水與微生物的生物反應(yīng)過(guò)程，僅考慮反應(yīng)器內(nèi)污水和污泥的流體力學(xué)行為。假設(shè)反應(yīng)器內(nèi)污水及污泥的流動(dòng)過(guò)程是連續(xù)介質(zhì)與分散介質(zhì)的兩相流動(dòng)和分離過(guò)程，采用RNGk-ε湍流模型和歐拉多相流模型模擬反應(yīng)器內(nèi)液-固兩相流動(dòng)過(guò)程。忽略流體的密度脈動(dòng)、顆粒質(zhì)量變化率的脈動(dòng)以及非定常關(guān)聯(lián)項(xiàng)，并將各相方程中的擴(kuò)散相統(tǒng)一表示成梯度的形式，雙流體模型中各相的控制方程組可以表示為式(1)所示的統(tǒng)一形式[14]。式(1)加上一些構(gòu)成源項(xiàng)和輸運(yùn)系數(shù)的模型方程就構(gòu)成了雙流體模型的封閉控制方程組。

(1)

2.2液-固曳力模型修正

在反應(yīng)器上部污泥稀相區(qū)，污水和污泥的運(yùn)動(dòng)狀態(tài)主要由連續(xù)相流動(dòng)決定，下部污泥密相區(qū)污泥的碰撞和團(tuán)聚對(duì)污水流動(dòng)狀態(tài)產(chǎn)生較大影響。在適用于稠密液-固流動(dòng)數(shù)值模擬的雙流體模型中，液-固兩相之間的相互作用和動(dòng)量交換主要通過(guò)液-固曳力模型來(lái)實(shí)現(xiàn)。考慮固體顆粒的流化、夾帶和輸運(yùn)的情況，建立準(zhǔn)確的液-固曳力模型，確保稠密液-固反應(yīng)器內(nèi)數(shù)值模擬結(jié)果的準(zhǔn)確程度。

采用自定義函數(shù)方法(User define function, UDF)，綜合考慮污泥顆粒的團(tuán)聚及聚并過(guò)程，建立適合密相和稀相共同存在的分區(qū)液-固曳力模型，并植入歐拉多相流模型，建立用于高效厭氧生物反應(yīng)器內(nèi)氣-液-固多相流動(dòng)模型。根據(jù)反應(yīng)器內(nèi)污泥流動(dòng)特點(diǎn)和空隙率分布，綜合考慮反應(yīng)器內(nèi)三相分離系統(tǒng)和反應(yīng)器下部密相鼓泡流化床的污泥分布特點(diǎn)，將整個(gè)計(jì)算區(qū)域劃分為4個(gè)區(qū)域，各區(qū)采用不同液-固兩相曳力模型，如表1所示[15-18]。

表1　ABR流動(dòng)區(qū)域劃分及各區(qū)曳力模型

2.3反應(yīng)器內(nèi)流動(dòng)相模擬計(jì)算條件確定

根據(jù)反應(yīng)器內(nèi)污水和污泥的基本流動(dòng)過(guò)程，可以得到反應(yīng)器內(nèi)液、固兩相流動(dòng)的計(jì)算工況和邊界條件。反應(yīng)器進(jìn)水流量為2.9 L/d，反應(yīng)器下部污泥密相床層施加1.2 L藏量，污泥質(zhì)量濃度為2.57 g/L，顆粒密度為1050 kg/m3，粒徑為9 μm。反應(yīng)器從底部進(jìn)水，上升流速為0.031 m/h。根據(jù)實(shí)驗(yàn)條件僅最上部出水口打開(kāi)，故反應(yīng)器上部Outlet-water 1排水口施加常壓出口邊界條件。反應(yīng)器殼體、其他出口和內(nèi)構(gòu)件的壁面施加無(wú)滑移固壁邊界條件以及二層壁面函數(shù)模型。

3結(jié)果與討論

3.1ABR內(nèi)多相流動(dòng)數(shù)值模擬結(jié)果

3.1.1速度場(chǎng)模擬結(jié)果

圖2為ABR內(nèi)液-固兩相流速度分布云圖。由圖2可知，污水從反應(yīng)器底部進(jìn)入，經(jīng)布水板射流進(jìn)入污泥床層；在液相曳力作用下，污泥床層從團(tuán)聚狀態(tài)開(kāi)始膨脹分散；顆粒小的污泥隨污水向上運(yùn)移，射流影響消失后污泥床層重新開(kāi)始團(tuán)聚；顆粒大的污泥無(wú)法受污水曳力作用繼續(xù)上移，開(kāi)始沉降，沉降至反應(yīng)器底部，再次受污水射流作用發(fā)生膨脹分散，向上運(yùn)動(dòng)。在污泥床層膨脹分散和團(tuán)聚而上下運(yùn)動(dòng)的過(guò)程中，污水與絮狀污泥顆粒的充分混合，滿足了污泥中微生物與污水的充分接觸，強(qiáng)化了傳質(zhì)作用，實(shí)現(xiàn)污水有機(jī)質(zhì)與污泥中微生物的生化反應(yīng)過(guò)程。

圖2　ABR內(nèi)液、固兩相速度分布云圖

從圖2同時(shí)發(fā)現(xiàn)，在反應(yīng)器底部布水板作用下，污水射入污泥床層，液、固兩相在底部出現(xiàn)局部速度較高區(qū)域；沿反應(yīng)器高度向上，液、固兩相速度受底部射流影響降低，速度值減小，沿反應(yīng)器徑向分布速度逐漸均勻。三相分離器內(nèi)液相和固相的速度都非常低，由于固相密度大于液相，因此進(jìn)入三相分離器的固相較少。同時(shí)，在反應(yīng)器底部射流影響區(qū)污水與污泥傳質(zhì)作用較強(qiáng)，上部污泥濃度降低，污水僅受反應(yīng)器幾何結(jié)構(gòu)約束，速度較低。

3.1.2壓力場(chǎng)和湍動(dòng)能模擬結(jié)果

圖3為ABR內(nèi)壓力和湍動(dòng)能分布云圖。壓力分布云圖反映了污水進(jìn)入反應(yīng)器內(nèi)與污泥相互作用的流動(dòng)特性。由圖3(a)看到，由于底部污泥床層的存在，污泥和污水產(chǎn)生的流體靜壓力在反應(yīng)器底部最大；由于底部污水的沖擊和曳力的帶動(dòng)，污泥床層產(chǎn)生膨脹，使得壓力分布整體趨勢(shì)沿軸向高度方向逐漸降低，局部有壓力波動(dòng)。湍動(dòng)能是衡量湍流發(fā)展和衰退的重要指標(biāo)。由圖3(b)看到，湍動(dòng)能云圖表現(xiàn)出與流體速度分布相似的分布趨勢(shì)；在反應(yīng)器底部的污水射流區(qū)，流體湍動(dòng)劇烈，污水與污泥的傳質(zhì)作用較強(qiáng)，介質(zhì)間混合充分，污泥隨污水向上運(yùn)移，沿反應(yīng)器高度方向流體流速不斷降低，污泥在重力作用下回落，湍動(dòng)降低，傳質(zhì)作用減弱。

圖3　ABR內(nèi)壓力和湍流動(dòng)能云圖

3.1.3濃度場(chǎng)分布模擬結(jié)果

采用污泥的體積分?jǐn)?shù)表示ABR內(nèi)污泥濃度隨時(shí)間變化情況，結(jié)果示于圖4。由圖4可知，初始狀態(tài)(0 s)下，向反應(yīng)器內(nèi)施加一定體積污泥，保證床層沿高度方向濃度分布一致。計(jì)算時(shí)間為1 s時(shí)，污泥床層發(fā)生濃度變化，部分污泥團(tuán)聚下沉，部分污泥受水流作用向上運(yùn)移，污泥床層出現(xiàn)下濃上稀的分布趨勢(shì)；5 s時(shí)，底部污泥床層在污水射流作用下，出現(xiàn)類似沸騰的污泥膨脹現(xiàn)象，下部密相床層濃度分布不再均勻，同時(shí)上部稀相區(qū)范圍繼續(xù)擴(kuò)大，出現(xiàn)了明顯的床層膨脹現(xiàn)象；隨計(jì)算時(shí)間繼續(xù)延長(zhǎng)，底部密相床層膨脹現(xiàn)象更加劇烈，同時(shí)稀相區(qū)范圍向上繼續(xù)擴(kuò)大，越往上濃度越稀，并且稀相區(qū)向上運(yùn)移的速度減慢，說(shuō)明在遠(yuǎn)離污水射流作用時(shí)，污泥受重力作用逐漸回落。污泥床層增高速度減慢，最終達(dá)到平穩(wěn)狀態(tài)，出現(xiàn)底部密相區(qū)、中間過(guò)渡區(qū)、上部稀相區(qū)和極稀相區(qū)的分布趨勢(shì)。污水與污泥的傳質(zhì)作用計(jì)算需要較小時(shí)間步長(zhǎng)，用以精準(zhǔn)預(yù)測(cè)污水與污泥的真實(shí)傳質(zhì)作用。這里僅顯示計(jì)算到35 s時(shí)的濃度分布圖，但反應(yīng)器內(nèi)的基本流動(dòng)過(guò)程已經(jīng)穩(wěn)定。

3.2ABR結(jié)構(gòu)改進(jìn)及優(yōu)化

從多相流動(dòng)的數(shù)值模擬結(jié)果發(fā)現(xiàn)，反應(yīng)器結(jié)構(gòu)雖具有較好的流動(dòng)分布，但為了強(qiáng)化污水與污泥傳質(zhì)效應(yīng)，可在一些結(jié)構(gòu)參數(shù)上進(jìn)行優(yōu)化改進(jìn)。污泥沉降效果受反應(yīng)器空間以及停留時(shí)間限制。反應(yīng)器運(yùn)行過(guò)程中，上升流速較低，上部三相分離器受上升流速影響慣性分離作用減弱，固相模擬計(jì)算結(jié)果與反應(yīng)器實(shí)際運(yùn)行中出現(xiàn)的跑泥現(xiàn)象吻合，需進(jìn)一步優(yōu)化改進(jìn)反應(yīng)器出水區(qū)結(jié)構(gòu)參數(shù)。

3.2.1反應(yīng)器結(jié)構(gòu)參數(shù)改進(jìn)

圖5為ABR出水區(qū)結(jié)構(gòu)參數(shù)改進(jìn)前后對(duì)比。由圖5可見(jiàn)，改進(jìn)后反應(yīng)器出水區(qū)排水環(huán)隙截面積明顯增加，運(yùn)移到環(huán)隙附近的污泥上升流速將明顯降低，受重力作用得到充分沉降，回流到反應(yīng)區(qū)，達(dá)到降低反應(yīng)器跑泥現(xiàn)象的目的。

3.2.2反應(yīng)器結(jié)構(gòu)參數(shù)改進(jìn)結(jié)果

圖6為ABR結(jié)構(gòu)改進(jìn)前后液、固兩相速度分布云圖。由圖6看到，液、固兩相速度分布均為反應(yīng)器下部布水板上方速度較高，且沿軸向方向速度逐漸減低。改進(jìn)后環(huán)隙的液、固兩相速度明顯低于改進(jìn)前，有利于固相的重力沉降回落。

圖4　ABR內(nèi)污泥濃度場(chǎng)分布圖

圖5　ABR出水區(qū)結(jié)構(gòu)參數(shù)改進(jìn)前后對(duì)比

圖6　ABR改進(jìn)前后液、固兩相流速分布對(duì)比

圖7為ABR改進(jìn)前后出水區(qū)速度云圖。從圖7可發(fā)現(xiàn)，結(jié)構(gòu)改進(jìn)前由于環(huán)隙截面積略大于反應(yīng)器主體截面積，從量級(jí)上可知其環(huán)隙流速與反應(yīng)器主體流速基本一致，污泥難以在有限的環(huán)隙空間再次發(fā)生沉降，極易被水流帶出反應(yīng)器；結(jié)構(gòu)改進(jìn)后，增加了出水區(qū)環(huán)隙截面積，同時(shí)增加45°擋板，在保證環(huán)隙截面積不小于反應(yīng)器主體截面積基礎(chǔ)上，增加污泥沉降空間，使環(huán)隙內(nèi)的流速明顯低于底部流體上升速度，增加了污泥在環(huán)隙內(nèi)的沉降效率，此外環(huán)隙下部擋板可截流一部分上升污泥和從環(huán)隙沉降落回?fù)醢逑路降奈勰啵岣呶勰喑两敌Ч?/p>

圖7　ABR改進(jìn)前后出水區(qū)速度云圖對(duì)比

圖8為ABR改進(jìn)前后出水區(qū)內(nèi)部液相流線分布。由圖8可知，流體從反應(yīng)器底部上升流經(jīng)三相分離器與反應(yīng)器出水區(qū)時(shí)，流速略有提高，并向出水口環(huán)隙轉(zhuǎn)向。改進(jìn)前出水口環(huán)隙速度較為均勻，不利于污泥再次沉降分離；改進(jìn)后增加的環(huán)隙空間降低了污泥在環(huán)隙空間的上升流速，增加了污泥的沉降空間和時(shí)間，同時(shí)下部45°擋板結(jié)構(gòu)也利于阻擋污泥上升運(yùn)動(dòng)，緩解反應(yīng)器跑泥問(wèn)題。

圖8　ABR改進(jìn)前后出水區(qū)液相流線速度分布

圖9為ABR改進(jìn)前后器內(nèi)整體污泥濃度分布。由圖9可以看出，改進(jìn)前、后反應(yīng)器內(nèi)污泥分布趨勢(shì)一致，包括底部的污泥濃相區(qū)和上部的污泥稀相區(qū)。然而，污泥進(jìn)入改進(jìn)前的反應(yīng)器出水口環(huán)隙后，污泥濃度在環(huán)隙內(nèi)較為均勻，表明污泥會(huì)隨液體流出；反應(yīng)器結(jié)構(gòu)改進(jìn)后，出水口環(huán)隙的污泥沉降作用增強(qiáng)，污泥濃度明顯降低，出水口附近基本無(wú)跑泥現(xiàn)象，說(shuō)明反應(yīng)器結(jié)構(gòu)改進(jìn)后對(duì)污泥的截留作用明顯提升。

4結(jié)論

(1)采用分段曳力模型，并基于RNGk-ε湍流模型和歐拉多相流模型，建立了適合高效厭氧生物反應(yīng)器內(nèi)的多相流動(dòng)模型。計(jì)算得到的反應(yīng)器內(nèi)部流場(chǎng)信息與反應(yīng)器實(shí)際流動(dòng)過(guò)程相吻合，驗(yàn)證了湍流模型的準(zhǔn)確性，并能夠反映反應(yīng)器內(nèi)多相流的傳質(zhì)過(guò)程。高效厭氧生物反應(yīng)器內(nèi)的流動(dòng)特點(diǎn)是，底部污水經(jīng)布水板小孔產(chǎn)生較強(qiáng)的射流作用，沖擊污泥床層，使團(tuán)聚的污泥膨脹分散，并與污水充分混合，強(qiáng)化傳質(zhì)，增加污水與污泥中微生物的接觸頻率；同時(shí)，在污水上升流速帶動(dòng)作用下，污泥沿反應(yīng)器高度上升運(yùn)動(dòng)，上升到一定高度，污水的沖擊作用減弱，污泥重新團(tuán)聚，并在重力作用下開(kāi)始回落，少量密度低的污泥隨污水繼續(xù)運(yùn)移，有可能從出水口隨水流排出。

圖9　ABR改進(jìn)前后的污泥濃度分布

(2)反應(yīng)器出水區(qū)結(jié)構(gòu)改進(jìn)后，跑泥現(xiàn)象明顯降低，污泥沉降作用增強(qiáng)，提高了反應(yīng)器內(nèi)多相流傳質(zhì)效果，模擬計(jì)算結(jié)果與反應(yīng)器實(shí)際運(yùn)行效果一致。Fluent模擬優(yōu)化反應(yīng)器結(jié)構(gòu)可以有效提升反應(yīng)器運(yùn)行效能，可用于反應(yīng)器流體力學(xué)的性能優(yōu)化。

符號(hào)說(shuō)明：

k——湍動(dòng)能，m2/s2；

p——壓力，Pa；

Sφ——因變量φ的源項(xiàng)；

t——時(shí)間，s；

v——流體速度，m/s；

VOF——污泥體積分?jǐn)?shù)，%；

xj——通用坐標(biāo)，mm；

α——體積分?jǐn)?shù)，%；

β——孔隙率，%；

Γφ——輸運(yùn)系數(shù)，m2/s；

ρ——流體密度，kg/m3；

φ——通用變量；

下標(biāo)

g——?dú)庀唷?/p>

j——笛卡爾坐標(biāo)系；

k——第k相。

參考文獻(xiàn)

[1] 羅毅. 中國(guó)環(huán)境統(tǒng)計(jì)年報(bào)2010[M].第一版. 北京：中國(guó)環(huán)境科學(xué)出版社, 2011:1.

[2] 2013—2017年中國(guó)環(huán)保行業(yè)研究及市場(chǎng)投資決策報(bào)告[M/CD].http://www.askci.com/reports/2011/11/1594826 142682.shtml, 2012, 37.

[3] VAN LIER J B, MAHMOUD N, ZEEMAN G. Anaerobic wastewater treatment[J].Biological Wastewater Treatment Principles Modeling & Design, 2008, 3(4):346-362.

[4] 賀延齡. 廢水的厭氧生物處理[M].北京：中國(guó)輕工業(yè)出版社，1998.

[5] 陳小光, 鄭平. 超高效螺旋式厭氧生物反應(yīng)器流態(tài)研究[J]. 環(huán)境科學(xué)學(xué)報(bào)，2010, 30(5): 941-946. (CHEN Xiaoguang, ZHENG Ping. Flow patterns of a super-high-rate spiral anaerobic bioreactor[J]. Acta Scientiae Circumstantiae, 2010, 30(5): 941-946.)

[6] 朱曉玲, 劉永紅, 劉宗寬, 等. 流態(tài)研究推動(dòng)下的高效厭氧生物反應(yīng)器的發(fā)展[J].榆林學(xué)院學(xué)報(bào)， 2008, 18(6):63-66.(ZHU Xiaoling, LIU Yonghong, LIU Zongkuan, et al. The development of high-rate anaerobic apparatus promoted by the study of flow pattern of such reactors[J].Journal of Yulin College, 2008, 18(6):63-66.)

[7] 楊積德, 陳曉娟. 厭氧工藝在低濃度廢水處理中的應(yīng)用[J].環(huán)境保護(hù)與循環(huán)經(jīng)濟(jì)， 2012, (3):51-54.(YANG Jide, CHEN Xiaojuan. Anaerobic process in the application of the low concentration wastewater treatment[J].Environmental Protection and Circular Economy, 2012, (3):51-54.)

[8] 何連生, 朱迎波, 席北斗, 等. 高效厭氧生物反應(yīng)器研究動(dòng)態(tài)及趨勢(shì)[J]. 環(huán)境工程，2004, 22(1): 7-11. (HE Liansheng, ZHU Yingbo, XI Beidou, et al. Research development and trends of high efficient anaerobic reactors[J]. Environmental Engineering, 2004, 22(1): 7-11.)

[9] 劉鋒平， 李薇， 李繼強(qiáng)， 等. HRT對(duì)UASB-SMBR(PTFE)組合工藝處理某油田含油廢水性能的影響[J].石油學(xué)報(bào)(石油加工)， 2012, 28(6):1053-1060.(LIU Fengping, LI Wei, LI Jiqiang, et al. The effect of HRT on UASB-SMBR (PTFE) process in treating oily wastewater of one oil field[J].Acta Petrolei Sinica (Petroleum Processing Section), 2012, 28(6):1053-1060.)

[10] 李波, 張慶文, 洪厚勝, 等. 攪拌反應(yīng)器中計(jì)算流體力學(xué)數(shù)值模擬的影響因素研究進(jìn)展[J].化工進(jìn)展， 2009, 28(1):7-11.(LI Bo, ZHANG Qingwen, HONG Housheng, et al. Several factors of CFD numerical simulation in stirred tank[J].Chemical Industry and Engineering Progress, 2009, 28(1):7-11.)

[11] 朱紅鈞, 林元華, 謝龍漢. FLUENT流體分析及仿真實(shí)用教程[M].北京：人民郵電出版社， 2010.

[12] 韓占忠, 王敬, 蘭小平. FLUENT—流體工程仿真計(jì)算實(shí)例與應(yīng)用[M].北京：北京理工大學(xué)出版社， 2010.

[13] 張凱, 王瑞金, 王剛. Fluent技術(shù)基礎(chǔ)與應(yīng)用實(shí)例[M].北京：清華大學(xué)出版社， 2010.

[14] 馮留海， 王江云， 毛羽， 等. 突擴(kuò)突縮管內(nèi)液-固沖蝕的數(shù)值模擬[J].石油學(xué)報(bào)(石油加工)， 2014, 30(6):1080-1085.(FENG Liuhai, WANG Jiangyun, MAO Yu, et al. Numerical simulation of liquid-solid erosion in the sudden expansion and contraction tube[J].Acta Petrolei Sinica (Petroleum Processing Section), 2014, 30(6):1080-1085.)

[15] ERGUN S. Fluid flow through packed columns[J].Chemical Engineering Progress, 1952, 48(2):89-94.

[16] GIDASPOW D. Multiphase Flow and Fluidization: Continuum and Kinetic Theory Description[M].San Diego: Academic Press， 1994.

[17] WEN C Y, YU Y H. A generalized method for predicting the minimum fluidization velocity[J].AIChE Journal, 1966, 12(3):610-612.

[18] 曹斌. 大差異多元顆粒氣固流化床流動(dòng)規(guī)律的研究[D]. 北京：中國(guó)石油大學(xué)， 2006.

Structure Optimization and Turbulent Flow Characteristics inHigh Efficient Anaerobic Biological Reactor

WANG Yu1,2, WANG Jiangyun1, XU Shuangshuang3, CHEN Chunmao1, MA Jingyuan3, GUO Shaohui1

(1.StateKeyLaboratoryofHeavyOilProcessing,BeijingKeyLaboratoryofOil&GasPollutionControl,ChinaUniversityofPetroleum,Beijing102249,China; 2.BeijingSafetyEnvironmentalProtectionEngineeringTechnologyResearchInstitute,CNOOCEnergyTechnology&ServicesLimited，Beijing300457,China;3.ChinaKunlunContracting&EngineeringCorporation,Beijing100037,China)

Abstract：The turbulent flow characteristics and multiphase flow mass transfer process in the high efficient anaerobic biological reactor were studied by numerical simulation. Multiphase turbulent flow control equation was built based on RNG k-ε (renormalization group, RNG, k-ε model) turbulent model and Euler multi-phase model, on the basis of which, the gas-liquid-solid multiphase flow model was built by loading the user defined function (UDF) to modify the drag model in Fluent solver and implanting the Euler multi-phase model, and the flow characteristics were calculated and the structure in high efficient anaerobic biological reactor was optimized. The results confirmed that the flow characteristics calculated with modified multiphase model was in reasonable agreement with experimental data, which verified the modified model and calculating method. And the multiphase mass transfer effect was improved obviously in modified reactor. An accurate prediction during the mass transfer process in the reactor could be obtained by the modified multiphase flow model, which could be used for the engineering prediction of high efficient anaerobic biological reactor.

Key words：high efficient anaerobic biological reactor; numerical simulation; multiphase flow model; structure optimization

收稿日期：2015-03-10

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目(21306229)和中國(guó)昆侖工程公司項(xiàng)目(2013GJTC-06-03)資助

文章編號(hào)：1001-8719(2016)03-0614-08

中圖分類號(hào)：X703.3

文獻(xiàn)標(biāo)識(shí)碼：A

doi:10.3969/j.issn.1001-8719.2016.03.024

第一作者： 王宇，女，博士，從事稠油加工污水的高效厭氧處理工藝方面的研究

通訊聯(lián)系人： 郭紹輝，男，教授，博士，從事石油工業(yè)環(huán)境污染治理與修復(fù)領(lǐng)域方面的研究；Tel:010-89732278;E-mail:cupgsh@163.com； 王江云，男，助理研究員，博士，從事多相流動(dòng)的數(shù)值模擬與實(shí)驗(yàn)方面的研究；Tel:010-89733293;E-mail:wangjy@cup.edu.cn