三相逆流化床的研究進展與展望

萬東玉，王勝偉，陳雙喜，王玉強

（山東玉皇化工有限公司，山東 菏澤 274512）

三相逆流化床（three-phase inverse fluidized bed，TPIFB）顆粒床層的膨脹靠向下的液流克服顆粒的靜浮力來實現，氣流逆向輔助流化[1]。三相逆流床的研究始于20世紀70年代，因其結構簡單、操作維護方便、混合傳熱傳質性能好和單位體積能耗低等一系列的優點，近年來在環境工程、食品化工、石油化工和生物化工等領域得到了廣泛的應用[2-9]。這些工業工程方面的應用引起了諸多學者對TPIFB流體力學行為及其應用方面的興趣和重視。本文主要目的是對TPIFB研究狀況及主要成果進行回顧，并對研究中存在的問題及今后應加強的方面提出看法和建議，以推動我國三相逆流化技術的持續發展。

1 TPIFB簡介

按照流體流動分布進行劃分，TPIFB主要由液體分離部分、液體分布部分、測試床體部分、氣體分布部分和氣體分離部分等五部分組成，其結構如圖1所示。氣體由塔底進氣上行，液體由塔頂分布下行，氣液逆流，共同作為反應器的能源來源。顆粒床層因其密度小于液體密度而漂浮在液體上，受下行液體的曳力作用克服液體的靜浮力作用而向下流化膨脹。氣體的引入可以促使床內液體密度下降，減小液體與顆粒間密度的差異，通過氣流擾動作用，促進了床內顆粒的流化。

圖1 三相逆流化床結構示意圖

TPIFB與傳統的三相流反應器相比具有以下突出的優勢[2-9]：① 反應器結構簡單，內部不存在機械動部件和復雜的密封件，可以防止生物反應外源污染，操作方便，維護費用低；②具有良好的流化效果，可以使較輕較小的固體顆粒完全保持懸浮狀態；③操作方式靈活，可以液體間歇操作，僅依靠 氣體鼓泡流化顆粒床層，也可以液體連續操作，顆粒床層的膨脹由氣液兩相流共同支撐；④在生物廢水處理中，液流速度低，能耗低，顆粒磨損和損耗小，能有效控制生物膜厚度；⑤氣液逆流操作，氣泡停留時間長，氣泡尺寸小，氣液傳質效率提高；⑥反應器內部軸向混合程度高，濃度、溫度分布均勻，動力學反應速率高；⑦反應器內部剪切應力小且分布均勻，熱量引入、移除容易，溫度控制方便，對剪切應力或溫度敏感的活性益菌和動植物的培養至關重要。

2 TPIFB基礎研究

隨著逆流化系統應用的不斷發展，TPIFB各方面的研究也隨之增加起來。其中，關于流型、壓力降、最小流化速度、相含率、停留時間分布、氣泡行為、傳質傳熱等流體力學特性及數值模擬方面的研究有了相關報道。

2.1 流型

三相逆流化床反應器內，一般認為液相為連續相，氣相為分散相。多年來，人們對不同的流型及流型轉變進行了大量研究，對流型特性有了一定認識[10-12]。根據TPIFB氣泡流動行為劃分出4個明確的流型區域，即氣泡分散的固定床區（伴隨或不伴隨部分流化）、鼓泡流化床區、過渡區和節涌流化床區[10]。如圖2所示，當氣液速度較低時，氣泡的湍動和液流作用于顆粒的曳力并不能克服顆粒的靜浮力，顆粒床層保持固定床狀態（圖中A點區域）。當增加液流速度或上升氣速，顆粒床層開始逐漸膨脹，直到床層全部流化。完全流化時的流化速度為最小流化速度，如圖中B點所示，此點速度下，顆粒全部運動，當小于此速度，床層變回固定床或部分流化床[11]。在鼓泡流化床區域（B與C點之間），氣泡均勻且床層均勻向下膨脹。床層膨脹隨氣速或液速增加而增加。當氣速上升超過C點，就進入了過渡區，在此區內氣泡的聚并、氣泡尺寸和氣泡頻率隨氣速增加而增加。圖中D點標明節涌流化床區域的開始，節涌氣泡上部的顆粒與氣泡一起向上運動，然后氣泡在顆粒床層表面崩破，節涌氣泡周圍的顆粒隨之降落。

曹長青等[12]通過增加表觀氣速，根據可視化觀察和床層界面平均固含率軸向分布規律確定了逆流三相湍動床內的不同流型及床層初始膨脹氣速UG1、完全膨脹氣速UG2和均勻分布氣速UG3等3個特征表觀氣速，并指出增加液體黏度，3個特征 表觀氣速都減小，增加初始床高，UG1幾乎不變，UG2和UG3減小。

圖2 三相逆流態化流型區域圖

2.2 壓力降

壓力降是計算三相逆流化床內所用能量及動力的主要參數之一，可以用來確定流化床反應器內的摩擦因子和最小流化速度。通常反應器內的壓力降隨液體流量的增加而增加，直到某一特殊的液體流量后，在非常大的液體流量范圍內，壓力降保持不變[13]。顆粒之間的碰撞頻率及其對壁面的碰撞壓力對床內壓力降具有很大的影響。Buffière等[14]利用操作參數關聯的模型考察了床內顆粒碰撞頻率和碰撞壓力，見式（1）。

式中，對于顆粒碰撞頻率模型，a=305700，b=1.01，c=1.02，d=0.99；對于顆粒壓力模型，a=423，b=1，c=1.14，d=0.75。

顆粒碰撞頻率和碰撞壓力直接影響著顆粒磨損率，也直接影響著能量損耗[14]。盡管如此，據文獻查閱了解，并沒有顆粒碰撞壓力和顆粒磨損率之間關系的報道，這也許會成為將來集中研究的一個方向。Lee[15]研究發現，增加液速可以增加壓力降，但超出最小流化液速之后，由于床層膨脹和流化的開始，壓力降隨液速增加而減小。要獲得每個縱向高度均勻間隔上的壓力梯度，可以假設整個床高上的空隙率、顆粒特性、表觀液速恒定。相比兩相系統中的壓力降正比于初始床高和液體黏度[16]，三相系統中的全床壓力梯度（靜壓梯度）分別受氣液固三相含率和密度的影響，其數學表達式為式（2）。

不同流型區域下的壓力梯度不同，而關于床內不同流型區域下局部和全床壓力梯度的研究及模型化研究的報道近乎沒有，考慮到能量的需要和動力的消耗，對各種流型下壓力降的研究還會成為研究者的研究內容之一。

2.3 最小流化速度

最小流化速度是三相流化床系統設計的重要動力學參數[17]。定義為床內壓力梯度最小時的表觀速度[18-19]，也定義為床內顆粒所受重力與液流向下作用的曳力之和平衡顆粒所受液體向上的浮力時的最小表觀速度[20]。Das等[20]利用4種不同粒徑、形狀和密度的顆粒和4種不同濃度的非牛頓液體，考察了兩種不同塔徑的三相逆流化床內的最小流化速度，結果發現固體加入量增加，床層壓力降增加，最小流化速度不變；液相黏度增加，阿基米德數減小，最小流化速度減小；塔徑減小，最小流化速度不變；可以通過雷諾數、阿基米德數和顆粒屬性（粒徑、密度和球形度）等來關聯確定最小流化速度。與Das等[20]研究結果不同的是，Han等[21]研究發現TPIFB臨界流化速度隨顆粒加入量增加而減小，并基于Comte等[22]提出的模型，修正了臨界流化速度的表達式，見式（3）。

Rengangathan等[23]研究了兩相和三相流化床系統各種變量下的最小流化速度，分別提出了TPIFB最小流化速度的經驗公式和修正的氣液擾動模型，并利用其他文獻數據比較分析發現誤差不到20%，證明了該公式和修正模型的基本可行性。Cho等[18]利用床內壓力梯度確定了最小流化液速，并發現最小流化液速依賴于床內顆粒與氣液混合物之間的密度差，表觀氣速增加，氣液混合密度減小，顆粒更易膨脹流化，最小流化速度減小。三相逆流化系統中，床層開始流化時的液速依賴于氣速，氣速越大，床層開始流化需要的液速越小[13，18，20，24]。

2.4 相含率

相含率是三相流化床的重要傳遞特性之一，決定著流化床內的流化效率。目前三相含率的研究中不乏各相含率的測量方法。氣含率的測量方法主要有體積法（或床層膨脹法）、壓差法（包括正壓差法和倒置U形壓差法）、電化學法（包括電阻法、電導法和電阻抗法）、光學法（激光探頭法、γ射線吸收法、光導纖維法和照相法）等四大類[25]。一般采用電導探頭法可以獨立測量總體液含率和局部液含率，在稀相操作時，也可以通過直接取樣法測量液含率。固含率可以通過目測法測得的床層高度結合式（4）求算。對于三相流化系統，已知其中兩相含率，利用三相含率之和為1的關系便可計算出第三相含率。

影響相含率的因素是多方面的，相互關系也錯綜復雜。Hamdad等[26]研究發現，增加表觀氣速，床層完全膨脹后，固含率保持不變，液含率增加到一最大值后減小，氣含率持續增加；添加乙醇表面活性劑可以減小床層流化需要的氣速，氣含率增加，液含率大幅減小[27]。氣含率隨表觀氣速增加而增 加[9，18-19，28-29]， 隨表觀液速增加而增加[18，28-29]。固含率隨氣速和液速增加而減小[29]。一定氣速下，增加液速，液含率增加，固含率減小，這是床層流化條件下床層膨脹的結果[18，29-30]。Shin等[29]指出液含率隨氣液表觀速度的影響很復雜，隨液體黏度增加而增加，而氣含率和固含率隨液體黏度增加而減小，并利用物性和操作條件修正了氣液相含率表達式，見式（5）、式（6）。

與Shin等[29]研究結果相矛盾的是，Kim等[27]發現氣含率隨液體黏度增加而增加，并利用操作條件和流體物性得到了氣含率的關聯式，見式（7）。

相對流體物性和操作條件對三相逆流床內相含率影響的研究，床內相含率局部分布的研究及報道較少。曹長青等[31]利用微電導探針技術測定和分析了氣液固逆流三相湍動床局部相含率的分布，研究發現顆粒軸向均勻分布時，局部固含率呈現床中心處小，床壁處大的分布規律，且在r/R=0.75～0.90范圍內有極大值，極大值出現的范圍與傳統流化床相比增大；局部氣含率呈現床中心處最大，床壁處最小的分布規律，在徑向上的非均勻性分布與傳統流化床相比愈加明顯。

2.5 停留時間分布

停留時間分布可以診斷反應器中是否存在不良流動，以便針對存在的問題改進反應器結構，而這方面的研究報道很少。Sánchez等[32]利用氯化鉀示蹤劑和瞬態響應技術考察了床內固含率對停留時間分布和液體混合的影響，結果發現平均停留時間和方差及液體混合不受固含率的影響。Renganthan 等[33]利用脈沖跟蹤技術和反卷積分析法確定了停留時間分布和液相軸向擴散系數，見式（8）、式（9）。

另外，Renganthan等[33]考察了表觀液速、靜態床層高度和阿基米德數對停留時間分布的影響，結果發現表觀液速增加，促使系統湍動性和軸向擴散系數增加；靜態床高增加，液相軸向擴散系數保持不變，貝克萊數增加；增加阿基米德數，促使最小流化速度和表觀液速增加，進一步增加液相軸向擴散系數。Kim等[34]研究發現TPIFB黏性介質系統液體軸向擴散系數隨氣液表觀速度增加而增加，隨液體黏度增加而減小，關聯式見式（10）。

Lee等[35]通過隨機模型松弛法研究發現三相逆流化床內顆粒軸向離散系數Dp隨氣液表觀速度或顆粒密度增加而增加，隨液體黏度增加而減小，并利用貝克萊數關聯式表征顆粒軸向離散系數，見 式（11）。

2.6 氣泡行為

氣泡行為是三相流化床反應器研究開發和工程放大的重要參數之一，可用于床內操作區域劃分及相流動結構的確定。氣泡行為直接決定了床內相含率的軸徑向分布、相間作用及相間傳熱傳質行為。氣泡行為參數的測定方法大致分為內置法和外置法兩類。常用的內置法有電導法、電容法和光導纖維法等，外置法有X射線法、攝像法和激光法等。外置法對床內流場無干擾，結果直觀，但多用于流化速度低的二維床，不能用于直徑較大的三維床。內置法不但適用于二維床也適用于三維床，但這類方法會不同程度地干擾流場，給測量帶來誤差[36]。

三相逆流化床內，氣泡大小隨氣液表觀速度或液體黏度增加而增加，隨固體加入量增加而減小[49]；氣泡上升速度隨表觀氣速和液體黏度增加而增加，隨表觀液速增加而減小[27，37]，隨固體加入量增加而增加[37]；氣泡頻率隨氣液表觀速度增加而增加，隨液體黏度增加而減小[27，37]。Son等[37]利用電導技術研究發現，氣泡大小和上升速度服從對數正態分布，并基于氣體飄移通量利用流體物性和操作條件關聯了氣泡大小、上升速度和頻率，見式（12）～式（14）。

Hamdad等[26]研究發現，添加乙醇表面活性劑可以減小床層流化需要的表觀氣速，抑制氣泡聚并和表面流動性，促使產生小氣泡，降低氣泡上升速度，增加氣含率。這與傳統流化床內添加表面活性劑，在氣泡表面形成表面張力梯度，抑制氣泡聚并，減小氣泡大小和上升速度相一致[38]。然而，不同的表面活性劑對氣泡的影響不同。對于大氣泡而言，具有較大的氣泡滑移速度，像氯化鈉、磷酸氫鈉及異戊醇-苯甲酸混合液等添加劑可以穩定氣泡表面，抑制氣泡破裂，減小氣含率。所以，表面活性劑具有兩種作用：要么抑制氣泡破裂減小氣含率，要么抑制氣泡聚并增加氣含率[39]。

縱觀國內外文獻，關于三相逆流化床內不同軸徑向位置上氣泡行為局部分布規律的研究報道較少，不同結構類型的氣體分布板、塔徑大小及塔內壓力等對氣泡行為影響方面的研究具有很重要的意義，而這方面的文獻報道也很少，這些也許會成為三相逆流化床未來研究的方向之一。

2.7 傳質特性

相際間傳質是反應器或接觸器設計和放大研究中重要的參數之一。三相逆流化床反應器中的宏觀反應過程主要包括氣-液相間傳質及液-固相間傳質。氣液界面的傳質是一個比較復雜的過程，到目前為止已經提出了一些描述氣液相間傳質過程的機理模型，如雙膜理論模型、溶質滲透模型以及表面更新模型等，其中最為簡單實用的是Whitman提出的雙膜理論[40]。對于難溶氣體（如氧氣）的氣液傳質過程，液膜一側的阻力遠大于氣膜一側，又由于測量手段的限制，很難準確得到氣液相界面面積，通常以液相一側的體積傳質系數kLa表示氣液兩相間的傳質系數。

目前對三相逆流化床內氣液傳質系數的研究已經取得了較大的進展，很多研究者通過對各自反應器內傳質系數的研究提出了一些經驗或半經驗的公式來描述氣液傳質系數與各影響因素之間的關聯式。三相逆流化床內，氣液體積傳質系數kLa隨 氣液表觀速度增加而增加，但受表觀氣速影響 小[8，26，41-43]，隨液體黏度增加而減小[34，41]，隨固體加入量（初始床高）增加而減小[41，43]，隨顆粒粒徑增加增加[43]，隨顆粒密度增加而增加[34，43]。Kim 等[34]和Sivasubramanian[43]分別利用操作條件和物料物性關聯了氣液傳質系數，見式（15）～式（17）。

牛頓系統

非牛頓系統

關于固體顆粒存在及固含率對kLa的影響規律，不同研究者提出了不同的結論。Hamdad等[26]研究發現在顆粒存在時，kLa隨水中乙醇表面活性劑的加入而增加；顆粒不存在時，kLa隨水中乙醇表面活性劑的加入而減小。這是由于顆粒可以破壞氣泡表面的表面活性劑，促進氣泡聚并，增加了氣泡上升速度，減小了表面活性劑的影響，促進傳質；顆粒不存在時，表面活性劑抑制氣泡聚并，盡管增加了氣含率和相界面面積，但大幅減小了氣泡表面流動性，并抑制傳質。Fahim等[8]研究發現表面活性劑生產系統中，固體顆粒存在且粒徑大于3mm時，固體顆粒可以引起氣泡的破裂，增加氣液接觸面積，提高氣液傳質效率；通過kLa關聯式（18）發現顆粒存在的參數PS指數接近于0，相對于其他參數可以忽略不計，即顆粒的存在對kLa的影響 很小。

Sánchez等[32]利用氯化鉀示蹤劑和瞬態響應技術考察了固含率對氣液傳質系數的影響，指出生物顆粒密度接近液相密度時，氣液傳質系數不受固含率的影響，并可以進行液固擬均相處理來研究TPIFB氣液傳質[44]。

近幾年，國內為關于TPIFB液固相間傳質的研究鮮有報道。Nikov等[45]利用電化學法測定電極電流來考察液固傳質速率，發現液固傳質速率受表觀液速和顆粒大小影響很小，隨液體黏度或顆粒密度增加而減小，最后提出了描述液固傳質的關聯式，見式（19）。

2.8 傳熱特性

在微生物的食品發酵和生物廢水處理等工程中，對溫度的控制尤為重要。熱量傳遞直接關系到流化床反應器的設計、反應條件的選擇和反應進程的控制等，是獲得良好反應效果的關鍵因素之一，影響著床內的流體流動特性和質量傳遞特性。Myre等[9]考察了三相逆流化床內加熱表面與床層間的平均傳熱系數，發現平均傳熱系數隨氣液表觀速度增加而增加，隨表面活性劑乙醇的加入而減小，隨固體顆粒量的增加先增加到最大值后又減小。Cho 等[18]研究了液固兩相和氣液固三相逆流化系統中加熱表面與床層間的傳熱特性，發現傳熱系數隨表觀氣速增加而增加，隨表觀液速增加先增加到最大值后又減小，并且達到傳熱系數最大值時的表觀液速隨顆粒密度或表觀氣速增加而降低；密度相對較大的聚乙烯顆粒系統下的傳熱系數高于聚丙烯顆粒系統下的傳熱系數；液體間歇系統下的傳熱系數比液體連續系統下的傳熱系數略低；隨床層空隙率增加，傳熱系數具有最大值；最后，應用量綱為1分析分別對液固兩相、液體間歇和液體連續三相系統下的傳熱系數進行了關聯，見式（20）～式（22）。

關于三相逆流化床內氣液表觀速度、顆粒密度和床層空隙率等變量對床內傳熱特性的影響，Son等[46]進行了相類似的研究和闡釋。他們還指出液體黏度增加，傳熱表面附近熱邊界層厚度增加，導致傳熱系數減小，并關聯了傳熱系數，見式（23）。

2.9 數值模擬

計算流體力學應用數值計算的方法直接求解控制方程，可以實現難以用實驗手段精確描述的局部流場特性、傳遞特性和化學反應過程，大大降低實驗次數，具有花費少、速度快、信息完整和模擬能力強等優點，可以與實驗相互比較，使理論直接面向實際，對于工程設計、開發和放大具有重要的指導意義。盡管如此，國內外關于TPIFB流體力學和傳遞特性數值模擬方面的研究鮮有報道。沈志恒等[47]基于歐拉-歐拉雙流體模型，數值模擬了倒置液固流化床內液固兩相流動行為，預測了床內顆粒的速度、濃度分布以及空隙率的變化。萬東玉等[48]基于多相流歐拉模型，應用計算流體動力學軟件，數值模擬了三相逆流湍動床流體動力學特性，考察了床內顆粒的軸徑向速度和固含率分布規律。Narayanan等[49]結合實驗，將基質溶液假設為牛頓流體，建立平推流流動模型和3種動力學模型，利用開發的軟件數值模擬并分析了自行設計的TPIFB生物反應器的特性。結果發現，不同基質溶液流率或生物酶載體加載量下生物基質轉化率的數值結果和實驗數據很接近，最大偏差為9.5%。

3 TPIFB應用研究

TPIFB反應器是在鼓泡塔反應器基礎上發展起來的一類高效多相流反應器。近年來，隨著TPIFB反應器在環境工程、食品化工、石油化工和生物化工等領域的廣泛應用，其研究開始從各種基礎性研究逐步轉向工程應用性研究，特別是在工業廢水處理方面的應用研究較多。Campos-Díaz等[6]利用乳酸桿菌作培養菌株，聚丙烯顆粒作載體，實驗測得了不同流體流量下的床層空隙率，并利用雷諾數和阿基米德數建立了適用于有/無生物膜球形顆粒的床層空隙率數學模型。Narayanan等[49]對TPIFB生物反應器進行了優化改進，并利用實驗以納米顆粒作生物酶載體考察其操作性能。研究表明，引入垂直流管的TPIFB可以實現低基質流率和低反應器體積下高的生物基質轉化率，而且所需垂直流管長度隨載體顆粒加載量增加而減小，隨生物基質轉化率增加而增加。

利用氣凝膠作為油污的吸附劑，TPIFB反應器可以實現很高的除油率。Quevedo等[50]利用疏水性處理的納米級多孔氣凝膠顆粒作固體，考察了TPIFB油水處理的除油效率和能力，結果發現，TPIFB除油效率主要受顆粒尺寸、初始床高、空隙率和流體速度影響：顆粒越小、初始床高越高、床層空隙率越小或流體速度越低，除油效果越好。利用TPIFB可以實現油水濃度2000mg/L減小到10mg/L的除油能力。Quevedo等還指出傳統的氣凝膠吸附劑使用后作焚燒或掩埋處理會造成污染和浪費，可以利用沖壓法、溶劑清洗法和過熱蒸汽法處理循環利用，提高其經濟性和使用效率。由于含油廢水的特殊性，國內外利用TPIFB反應器對含油廢水進行處理的研究仍然較少，對這方面的研究開發仍需深入。

TPIFB反應器適用于各種生物降解的無機或有機廢水處理，對各類工業廢水及生活污水均有良好的處理效果。相比其他反應器，TPIFB反應器可以實現90%以上的COD去除率，并且能耗低，生物量濃度高及卸料方便[51]。Beristain-Cardoso等[2]在TPIFB反應器內研究了不同的CH3COO-與NO3-濃度比（0.85、0.72和0.62）對廢水中硫酸鹽、硝酸鹽和乙酸鹽去除效率的影響，研究表明：兩者濃度比為0.72時，TPIFB反應器去硫、脫氮及除碳效果最佳，硫酸鹽、硝酸鹽和醋酸鹽去除率分別達到100%，產物主要為單質硫、氮氣和二氧化碳。Rajasimman等[3]利用TPIFB生物反應器對淀粉廢水進行微生物好氧化處理，考察不同水力學停留時間和初始基質濃度下淀粉廢水的處理效果，結果發現：有機質加載量為1.35kg COD/(m3·d)時，COD去除率達到最大值95.6%；有機質加載量為26.73 kg COD/(m3·d)時，COD去除率達到最小值51.8%。在TPIFB生物反應器中，Papirio等[7]利用低密度聚丙烯顆粒作生物載體，考察TPIFB生物反應器內兩種COD/硫酸鹽計量比和兩種硫酸鹽濃度進料下的COD去除率和硫酸鹽去除率，結果發現，進料基質pH值為3的床層流化超過10%時的生物量難以固定，pH值增加到5，COD/硫酸鹽計量比為4.00時的COD最高達75%，硫酸鹽去除率達97%。Sowmeyan等[52]使用粒徑1mm濕重295kg/m3的珍珠巖顆粒作生物載體，在厭氧逆流化床反應器內對酒廠廢水進行處理中發現：這種珍珠巖顆粒對生物量生長非常有效，并且在有機加載量為35kg COD/(m3·d)時，可實現COD去除率84%，沼氣產量13.22L/(L·d)；逆流厭氧流化床反應器是厭氧處理酒廠廢水的理想選擇。

TPIFB反應器廢水處理研究中，除了考察進料組分濃度、pH值和顆粒加載量等參數對除廢效果的影響外，也不乏反應器操作模式和操作參數的選擇性研究。

利用TPIFB流化床反應器處理工業廢水制甲烷的過程中，生物膜生成和達到穩態的時間較長，并且偏向于間歇模式操作，而利用連續操作模式的研究較少。Alvarado-Lassman等[4]以葡萄糖作為碳源基質考察TPIFB反應器最優操作模式，結果發現連續操作比間歇操作穩態時間短、COD去除率高、載體吸附生物量多、甲烷產率高且更接近理論值。Sokól等[5]以高強度煉廠廢水作基質、聚丙烯顆粒作生物載體，考察TPIFB生物反應器在高強度廢水處理應用中最優的操作參數，結果發現，固定床與反應器體積比（Vb/Vr）0.55、氣速0.046m/s和操作時間65h時的操作最優，COD去除率可達到96%，而將這些操作參數引用到酚類廢水和其他烴類廢水處理上時，COD去除率可分別達到95%和90%。Cresson等[53]利用TPIFB對酒水廠COD/N/P濃度比為400/7/1的廢水進行處理時，發現使用生物載體流化的最低水力學強度（操作氣速）可以促進生物膜生長，增加操作氣速會抑制生物膜生長和降低生物膜厚度，并建議生物膜生長階段使用最小動力學強度，然后增加動力學強度控制生物膜厚度。Sabarunisha等[54]利用TPIFB作生物膜反應器，用熒光假單胞菌降解酚類廢水，考察生物降解過程中氣含率、COD去除率與酚類去除率等隨反應器動力學行為的影響，結果發現，對于不同的Vb/Vr，氣含率、COD去除率與酚類去除率隨表觀氣速增加而增加，在最優表觀氣速Ugm后隨表觀氣速增加而減小；在Vb/Vr為0.20、Ugm為0.220 m/s時操作最佳，此時氣含率、COD去除率與酚類去除率分別達到0.0653、98.5%與100%。

4 結 語

TPIFB反應器因其獨特的優點，在環境工程、石油化工和廢水處理等方面已得到研究者的廣泛關注，這方面的基礎性研究和應用性研究也取得了一定的進展。然而通覽迄今已報道的TPIFB的研究狀況，容易發現：①大部分研究處在穩態和宏觀的實驗研究階段，而瞬態和微觀角度的考察研究較少；②實驗研究主要集中在流體物性和操作參數的影響上；③實驗研究為主，計算流體力學數值模擬的應用較少；④實驗結論多為一般趨勢而非定量關系，有些甚至出現矛盾；⑤主要是冷態小規模實驗裝置研究，實驗獲得的模型一般存在局限條件；⑥應用性研究仍然很少，且研究內容很分散。因此，鑒于TPIFB工業應用發展的需要，本文作者認為以下幾個方面需要做進一步的研究，以期為其提供強有力的技術支持：①由于床內流體分布均勻性的差別，在考慮載體優選、曝氣方式和操作參數優化控制的同時，進一步考察通過優化反應器結構（高徑比）、設計流體分布器和引入外場（如聲場、磁場）來提高過程優化；②利用高端的測試技術結合計算流體力學數值模擬，考察與分析床內不同軸徑向位置上流體力學行為和傳熱傳質特性；③深入研究探討TPIFB流體力學行為、傳質傳熱特性和反應過程的強化機理，建立統一性和實用性強的機理模型和反應器模型；④考慮到TPIFB小規模試驗裝置與大規模工業化裝置之間有很大差別，進行TPIFB放大方面的研究及工業應用性研究，積累設計、運行和優化經驗。

符 號 說 明

Ar—— 阿基米德數

a—— 相界面面積，m2/m3

F—— 固含率函數

FB—— 氣泡上升頻率，Hz

Fr—— 動量轉移應力項

G—— 氣泡群速度，m/s

Ga—— 伽利略數

H—— 床高，m

h——傳熱系數，W/(m2·K)

M——顆粒加入量，kg

Nu——努塞爾數

RC——密度比

Sc——施密特數

Sh——舍伍德數

U——表觀速度，m/s

ε——相含率

θ——時間，s

μ——黏度，Pa·s

μeff——有效黏度，Pa·s

ρ——密度，kg/m3

σ——表面張力，N/m

υ——空隙率

下角標

B——氣泡

G——氣體

L——液體

P——顆粒

S——固體

0——初始值

[1] Fan L S. Gas-Liquid-Solid Fluidization Engineering[M]. Boston：Butterworths Publishers，1989：2-4.

[2] Beristain-Cardoso R，Texier A C，Sierra-álvarez R，et al. Simultaneous sulfide and acetate oxidation under denitrifying conditions using an inverse fluidized bed reactor[J].Journal of Chemical Technology and Biotechnology，2008，83（9）：1197-1203.

[3] Rajasimman M，Karthikeyan C. Performance of inverse fluidized bed bioreactor in treating starch wastewater[J].Frontiers of Chemical Engineering in China，2009，3（3）：235-239.

[4] Alvarado-Lassman A，Sandoval-Ramos A，Flores-Altamirano M G. Strategies for the startup of methanogenic inverse fluidized-bed reactors using colonized particles[J].Water Environment Research，2010，5：387-391.

[5] Sokó? W，Woldeyes B. Evaluation of the inverse fluidized bed biological reactor for treating high-strength industrial wastewaters[J].Advances in Chemical Engineering and Science，2011，1：239.

[6] Campos-Díaz K E，Bandala-González E R，Limas-Ballesteros R. Fluid bed porosity mathematical model for an inverse fluidized bed bioreactor with particles growing biofilm[J].Journal of Environmental Management，2012，104：62-66.

[7] Papirio S，Esposito G，Pirozzi F. Biological inverse fluidized-bed reactors for the treatment of low pH- and sulphate-containing wastewaters under different COD/SO42-conditions[J].Environmental Technology，2013，34（9）：1141-1149.

[8] Fahim S，Dimitrov K，Vauchel P，et al. Oxygen transfer in three phase inverse fluidized bed bioreactor during biosurfactant production by Bacillus subtilis[J].Biochemical Engineering Journal，2013，76：70-76.

[9] Myre D，Macchi A. Heat transfer and bubble dynamics in a three-phase inverse fluidized bed[J].Chemical Engineering and Processing：Process Intensification，2010，49（5）：523-529.

[10] Fan L S，Muroyama K，Chern S H. Some remarks on hydrodynamics of inverse gas-liquid-solid fluidization[J].Chemical Engineering Science，1982，37：1570-1572.

[11] Buffière P，Moletta R. Some hydrodynamic characteristics of inverse three phase fluidized-bed reactors[J].Chemical Engineering Science，1999，54（9）：1233-1242.

[12] 曹長青，劉明言，王一平，等. 液體黏度對氣液固逆流三相湍動床動力學行為的影響[J]. 化工學報，2005，56（4）：627-631.

[13] Bendict R J F，Kumaresan G，Velan M. Bed expansion and pressure drop studies in a liquid-solid inverse fluidised bed reactor[J].Bioprocess Engineering，1998，19（2）：137-142.

[14] Buffière P，Moletta R. Collision frequency and collisional particle pressure in three-phase fluidized beds[J].Chemical Engineering Science，2000，55（22）：5555-5563.

[15] Lee D H. Transition velocity and bed expansion of two-phase （liquid-solid） fluidization systems[J].Korean Journal of Chemical Engineering，2001，18（3）：347-351.

[16] Bendict R J F，Kumaresan G，Velan M. Bed expansion and pressure drop studies in a liquid-solid inverse fluidised bed reactor[J].Bioprocess Engineering，1998，19（2）：137-142.

[17] Zhu Zhiping，Na Yongjie，Lu Qinggang. Effect of pressure on minimum fluidization velocity[J].Journal of Thermal Science，2007，16（3）：264-269.

[18] Cho Y J，Park H Y，Kim S W，et al. Heat transfer and hydrodynamics in two- and three- phase inverse fluidized beds[J].Industrial& Engineering Chemistry Research，2002，41（8）：2058-2063.

[19] Lee D H，Epstein N，Grace J R. Hydrodynamic transition from fixed to fully fluidized beds for three-phase inverse fluidization[J].Korean Journal of Chemical Engineering，2000，17（6）：684-690.

[20] Das B，Ganguly U P，Das S K. Inverse fluidization using non-Newtonian liquids[J].Chemical Engineering and Processing：Process Intensification，2010，49（11）：1169-1175.

[21] Han H D，Lee W，Kim Y K，et al. Phase hold-up and critical fluidization velocity in a three-phase inverse fluidized bed[J].Korean Journal of Chemical Engineering，2003，20（1）：163-168.

[22] Comte M P，Bastoul D，Hebrard G，et al. Hydrodynamics of a three-phase fluidized bed—The inverse turbulent bed[J].Chemical Engineering Science，1997，52（21）：3971-3977.

[23] Renganathan T，Krishnaiah K. Prediction of minimum fluidization velocity in two and three phase inverse fluidized beds[J].The Canadian Journal of Chemical Engineering，2003，81（3-4）：853-860.

[24] Ibrahim Y A A，Briens C L，Margaritis A，et al. Hydrodynamic characteristics of a three-phase inverse fluidized-bed column[J].AIChE Journal，1996，42（7）：1889-1900.

[25] 薛勝偉. 氣升式環流反應器流動與傳質的研究[D]. 南京：南京工業大學，2005.

[26] Hamdad I，Hashemi S，Rossi D，et al. Oxygen transfer and hydrodynamics in three-phase inverse fluidized beds[J].Chemical Engineering Science，2007，62（24）：7399-7405.

[27] Kim S D，Kang Y. Hydrodynamics，heat and mass transfer in inverse and circulating three-phase fluidized-bed reactors for wastewater treatment[J].Stud. Surf. Sci. Catal.，2006，159：103-108.

[28] Liu M，Lu C，Shi M，et al. Hydrodynamics and mass transfer in a modified three-phase airlift loop reactor[J].Petroleum Science，2007，4（3）：91-96.

[29] Shin I S，Son S M，Kang Y，et al. Phase holdup characteristics of viscous three-phase inverse fluidized beds[J].Journal of Industrial and Engineering Chemistry，2007，13（6）：971-978.

[30] Bandaru K S，Murthy D V S，Krishnaiah K. Some hydrodynamic aspects of 3-phase inverse fluidized bed[J].China Particuology，2007，5（5）：351-356.

[31] 曹長青，劉明言，郭慶杰. 氣液固逆流三相湍動床局部相含率的分布[J]. 青島科技大學學報：自然科學版，2007，2：134-137.

[32] Sánchez O，Michaud S，Escudié R，et al. Liquid mixing and gas-liquid mass transfer in a three-phase inverse turbulent bed reactor[J].Chemical Engineering Journal，2005，114（1）：1-7.

[33] Renganathan T，Krishnaiah K. Liquid phase mixing in 2-phase liquid-solid inverse fluidized bed[J].Chemical Engineering Journal，2004，98（3）：213-218.

[34] Kim S W，Kim H T，Song P S，et al. Liquid dispersion and gas-liquid mass transfer in three-phase inverse fluidized beds[J].The Canadian Journal of Chemical Engineering，2003，81（3-4）：621-625.

[35] Lee K ll，Son S Mo，Kim U Yeong，et al. Particle dispersion in viscous three-phase inverse fluidized beds[J].Chemical Engineering Science，2007，62（24）：7060-7067.

[36] 郭奮，周游，王巖，等. 氣-固流化床的氣泡參數測定——直射式光纖探頭在流化床研究中的應用[J]. 化學反應工程與工藝，1990，2（6）：71-78.

[37] Son M S，Kang S H，Kim U K，et al. Bubble properties in three-phase inverse fluidized beds with viscous liquid medium[J].Chem. Eng. Process，2007，46：736-741.

[38] Gorowara R L，Fan L S. Effect of surfactants on three-phase fluidized bed hydrodynamics[J].Industrial & Engineering Chemistry Research，1990，29（5）：882-891.

[39] Briens C L，Ibrahim Y A A，Margaritis A，et al. Effect of coalescence inhibitors on the performance of three-phase inverse fluidized-bed columns[J].Chemical Engineering Science，1999，54（21）：4975-4980.

[40] Garcia-Ochoa F，Gomez E. Bioreactor scale-up and oxygen transfer rate in microbial processes：An overview[J].Biotechnology Advances，2009，27（2）：153-176.

[41] Jin R C，Zheng P，Chen J W，et al. Gas-liquid mass transfer in three-phase inverse turbulent bed reactor[J].Chinese Journal of Biotechnology，2007，23（1）：122-126.

[42] Nikolov V，Farag I，Nikov I. Gas-liquid mass transfer in bioreactor with three-phase inverse fluidized bed[J].Bioprocess Engineering，2000，23（5）：427-429.

[43] Sivasubramanian V. Gas-liquid mass transfer in three-phase inverse fluidized bed reactor with Newtonian and non-Newtonian fluids[J].Asia-Pacific Journal of Chemical Engineering，2010，5（2）：361-368.

[44] Zhao B，Wang J，Yang W，et al. Gas-liquid mass transfer in slurry bubble systems：I. Mathematical modeling based on a single bubble mechanism[J].Chemical Engineering Journal，2003，96（1）：23-27.

[45] Nikov I，Karamanev D. Liquid-solid mass transfer in inverse fluidized bed[J].AIChE Journal，1991，37（5）：781-784.

[46] Son S M，Lee K I，Kang S H，et al. Heat transfer coefficient in viscous three-phase inverse fluidized beds[J].AIChE Journal，2007，53（11）：3011-3016.

[47] 沈志恒，劉文鐵，金記英，等. 倒置液固流化床內液固兩相流動特性的數值模擬[J]. 過程工程學報，2008（s2）：394-397.

[48] 萬東玉，曲文凱，曹長青. 氣液固三相逆流湍動床流體動力學特性的三維數值模擬[J]. 青島科技大學學報：自然科學版，2014，35（3）：277-282.

[49] Narayanan C M，Monangi T，Prasad H，et al. Studies on performance analysis and computer aided design of inverse fluidized bed bioreactors with nanosupport particles[J].International Journal of Chemical Reactor Engineering，2014，12（1）：1-13.

[50] Quevedo J A，Patel G，Pfeffer R. Removal of oil from water by inverse fluidization of aerogels[J].Industrial & Engineering Chemistry Research，2008，48（1）：191-201.

[51] Arnaiz C，Buffiere P，Elmaleh S，et al. Anaerobic digestion of dairy wastewater by inverse fluidization：The inverse fluidized bed and the inverse turbulent bed reactors[J].Environmental Technology，2013，24（11）：1431-1443.

[52] Sowmeyan R，Swaminathan G. Performance of inverse anaerobic fluidized bed reactor for treating high strength organic wastewater during start-up phase[J].Bioresource Technology，2008，99（14）：6280-6284.

[53] Cresson R，Escudié R，Carrère H，et al. Influence of hydrodynamic conditions on the start-up of methanogenic inverse turbulent bed reactors[J].Water Research，2007，41（3）：603-612.

[54] Sabarunisha B S，Radha K V. Hydrodynamic behavior of inverse fluidized bed biofilm reactor for phenol biodegradation using Pseudomonas fluorescens[J].Korean Journal of Chemical Engineering，2014，31（3）：436-445.