曹廣斌，程啟云、2，韓世成，周煊亦、2，蔣樹義，陳忠祥

(1.中國水產(chǎn)科學(xué)研究院黑龍江水產(chǎn)研究所，黑龍江 哈爾濱150070;2.上海海洋大學(xué) 工程學(xué)院，上海201306)

在冷水魚工廠化養(yǎng)殖中，氨氮是養(yǎng)殖魚類的主要排泄物之一，鮭鱒魚類每增長1 kg 就會產(chǎn)生33 g 氮［1-2］，如不及時進行處理，水中氨氮積累到一定濃度就會對養(yǎng)殖魚類產(chǎn)生毒害作用。工廠化養(yǎng)殖水體的氨氮總量一般不應(yīng)超過1 mg/L，非離子氨不應(yīng)超過0.05 mg/L［3］。由于冷水魚養(yǎng)殖水溫較低(8 ～20 ℃)，常用的生物處理效率會大大降低［4］，很難滿足處理氨氮的要求。臭氧氧化具有氧化氨氮、凝聚懸浮物、消毒殺菌、脫色、除臭等水處理作用［5］，是處理養(yǎng)殖水體中氨氮的有效方法之一。用臭氧直接氧化氨氮其去除效率約為25.8%［6］，在加入催化劑的條件下，氧化效率可以大幅度提高［7］。但是，臭氧在水中超過一定濃度時又會對水生生物產(chǎn)生毒害，臭養(yǎng)對養(yǎng)殖魚類的安全濃度為0.008 ～0.060 mg/L［8］。在利用臭氧氧化過程中，必須使臭氧在水中能充分混合、反應(yīng)并減少殘留，以確保應(yīng)用過程的安全性。鼓泡塔能夠使氣液兩相或氣液固三相均勻地接觸流動，促進各相間的傳質(zhì)、溶解和反應(yīng)，被廣泛應(yīng)用于化學(xué)工業(yè)中的加氫、空氣氧化和費-托反應(yīng)等［9］，也用于對水體的臭氧消毒［10］。但是，鼓泡塔體積龐大，安裝、搬運困難，不適宜水產(chǎn)養(yǎng)殖中使用。

1 反應(yīng)塔的設(shè)計

為減小設(shè)備高度，增加易用性，促進臭氧與水體的氣液兩相流的傳質(zhì)與溶解，提高臭氧氧化氨氮的效率，降低水體中殘留的臭氧，本研究中，根據(jù)工廠化水產(chǎn)養(yǎng)殖水處理的實際情況，應(yīng)用氣泡流體力學(xué)的設(shè)計方法和鼓泡塔工作原理，設(shè)計了雙層逆流臭氧氧化反應(yīng)塔。

1.1 結(jié)構(gòu)原理

反應(yīng)塔結(jié)構(gòu)原理如圖1所示。工作時，養(yǎng)殖水體在循環(huán)水泵的作用下從進水管高速流入，臭氧發(fā)生器產(chǎn)生的臭氧通過臭氧入口管和微孔擴散器進入內(nèi)筒與水體混合，并在內(nèi)筒錐形部分受到水流沖擊而破碎，變成小氣泡，形成氣泡流。氣泡在內(nèi)筒體隨水流向下流動，并隨著壓力的增大，溶解度不斷升高，形成混合、溶解和反應(yīng)過程。在內(nèi)筒體底部設(shè)置4 個對稱排水槽，水氣混合體通過排水槽進入外部環(huán)形筒體內(nèi)，在外筒體內(nèi)以同樣的速度上升，并進行臭氧氧化氨氮的充分反應(yīng)。反應(yīng)后的水體從外筒體上部排出管排出，再用于水產(chǎn)養(yǎng)殖。未被溶解和反應(yīng)的臭氧從廢氣排出管排出，并進行無害化處理。在此過程中，養(yǎng)殖水體完成了臭氧氧化氨氮的處理過程，達到了去除氨氮的目的。

1.2 基于氣泡流體理論的反應(yīng)塔設(shè)計

設(shè)計計算反應(yīng)塔的各項尺寸，涉及到氣泡流體內(nèi)部的流動狀態(tài)、氣泡運動過程的質(zhì)量變化及其尺寸變化，須建立運動、質(zhì)量擴散和尺寸變化數(shù)學(xué)模型。

圖1 逆流型雙層反應(yīng)塔結(jié)構(gòu)簡圖Fig.1 Structure chart of counter－flow reaction column

1.2.1 氣泡運動狀態(tài)方程 反應(yīng)塔中的水流是確定速度的定常流，根據(jù)氣泡運動過程的受力情況，氣泡在反應(yīng)塔流體中的相對運動方程為

其中:mg為氣泡的質(zhì)量;u 為氣泡運動的速度;t為氣泡運動的時間;Fb為氣泡的浮力;Fg為氣泡的重力;Fd為阻力;Fa為附加質(zhì)量力;Fh為Basset 殘留阻力。由于氣體的密度ρg遠小于水的密度ρf，且在污染較少的水體中殘留阻力非常小，氣泡的質(zhì)量、氣泡的重力和氣泡運動的殘留阻力可忽略不計。整理后的運動方程如下［11］:

當(dāng)水流向下流動時，浮力與運動的方向相反，令u=u1，可以得到

當(dāng)水流向上流動時，浮力與運動的方向相同，令u=u2，可以得到

其中:k 為附加質(zhì)量系數(shù)(球形氣泡為0.5);V 為氣泡的體積;Ab為氣泡的截面積;CD為運動阻力系數(shù)［12］，且

式中:Re 為Reynolds 數(shù)，Re=(d·u)/ν;d 為流管直徑，ν 為水的運動黏度系數(shù)，在標準條件下ν=1.003 ×106m2/s。

1.2.2 氣泡運動過程中的溶解方程在臭氧氣泡溶解于水的過程中，任意時間氣泡的質(zhì)量變化等于氣泡界面的質(zhì)量擴散［11］。其關(guān)系表達式為

式(5)的左邊是氣泡質(zhì)量隨時間的變化率，右邊是任意半徑的氣泡界面向水中擴散的氣泡質(zhì)量。其中:R 為氣泡半徑;ρg為臭氧氣泡的密度(標準條件下為2.144 g/L);D 為臭氧在水中的擴散系數(shù);C 為擴散界面臭氧的濃度;Cs為氣泡表面臭氧的濃度(20 ℃時的飽和濃度為30 mg/L);C∞為水體中無窮遠處的臭氧濃度，即水體的臭氧濃度;r 為擴散界面的半徑;s 為與擴散半徑垂直的界面;θ 為擴散界面的方向角;sh(Sherwood Number)為質(zhì)量轉(zhuǎn)移系數(shù)，且

式中:Pe為Peclet 系數(shù)，且Pe=(2R·u)/D，D為臭氧擴散系數(shù)，標準條件下(T=293 K，P=1.013 ×105Pa，下同)，D=1.3 ×10-9m2/s［13］。

式(5)中，氣泡表面的臭氧濃度Cs可認為是臭氧在水中的飽和濃度，其計算公式可根據(jù)亨利定律表示為

其中:H 為隨著溫度變化的亨利常數(shù)(pH 為7、溫度為293 K 的標準條件下H=4.6 ×108Pa·g/L);P 為氣泡表面的壓力。H 也可按下式計算［14］:

式中，［OH-］表示OH-濃度，T表示溫度(K)。

1.2.3 氣泡運動過程中氣泡尺寸的變化方程 研究表明，當(dāng)氣泡直徑和運動的速度較小時，氣泡在水中保持球形形狀，氣泡內(nèi)外壓力應(yīng)保持平衡，其內(nèi)外壓力平衡關(guān)系如下［15］:

即

其中:Pg為氣泡內(nèi)部的壓力;P0為大氣壓;S 為氣泡離水面的距離;˙R 為氣泡直徑的變化速度，˙R=dR/dt。

當(dāng)臭氧氣泡半徑R ＞0.001 mm和水深S0＜10 m 時，氣泡內(nèi)壓力遠小于其臨界壓力Pc(Pc=5.05 ×106Pa)［16］，氣泡內(nèi)的臭氧氣體可視為理想氣體，其變化規(guī)律符合熱力學(xué)方程，即

其中，R0為氣泡的起始半徑。

1.3 水流速度的確定

在氣泡流運動過程中，其向下流動的速度應(yīng)該遠大于氣泡上升速度，才能確保氣泡在運動過程中不會聚集變成大氣泡，進而影響氣泡流體的均勻一致性。根據(jù)氣泡流體力學(xué)理論［12］，在氣流速度較大(約1 m/s)的密集鼓泡條件下，氣液混合物可形成大小相當(dāng)均勻的小氣泡，其上升的速度用下式估算:

其中:uc為氣泡上升的速度(m/s);σ 為液體的表面張力(N/m);g 為重力加速度(m/s2);μf為水的動力黏性系數(shù)(N·s/m3);ρf為水的密度(kg/m3)。

把水在20 ℃時的一些物理參數(shù):g=9.8 m/s2，ρf=1000 kg/m3，σ=0.0728 N/m，μf=1002 N·s/m3，代入式(12)可計算出氣泡由于浮力而引起的上升速度uc=0.02 m/s，取氣泡流體速度為uf=0.05 m/s。

1.4 反應(yīng)塔內(nèi)外筒直徑的計算

反應(yīng)塔內(nèi)筒直徑的計算公式為

其中:d 為內(nèi)筒直徑(m);Q 為循環(huán)泵的流量(m3/s);uf為流體的速度(m/s)。為了使氣泡流流動均勻穩(wěn)定，在設(shè)計過程中使氣泡流上、下流速相同，內(nèi)筒截面積與內(nèi)外筒之間的環(huán)形面積相等，可以得到外筒直徑。

1.5 反應(yīng)塔高度的設(shè)計計算

反應(yīng)塔選擇Ⅰ型微孔鈦板作為擴散器，其孔徑≤10 μm，氣泡平均直徑為0.33 mm［17］;在氧氣作為氣源產(chǎn)生臭氧的情況下，假設(shè)氧氣由電暈生成的100%臭氧。理想條件下的反應(yīng)塔高度應(yīng)該滿足如下要求:氣泡在內(nèi)筒向下流動和經(jīng)過內(nèi)外筒之間的空間向上流動到塔頂部的過程中，臭氧氣泡全部溶解，充分反應(yīng)沒有殘留。

利用氣泡流體理論建立的(2)、(3)、(5)和(10)式微分方程，進行反應(yīng)塔高度尺寸的設(shè)計計算。計算中運動時間t 與距離S 的計算公式為

當(dāng)R=0，S2=0 時，S=S1即為設(shè)計反應(yīng)塔高度。

1.5.1 計算的邊界條件 氣泡在水下開始上升時，其邊界條件:t0=0，R0=0.33 mm，u0=0.05 m/s。其中t0為起始時間，R0為氣泡的起始半徑;u0為氣泡流體的開始速度。計算方程中的常系數(shù)采用標準條件下的數(shù)值;由于氣泡在水中，可認為氣泡內(nèi)的溫度等于水的溫度，且水溫保持不變。

1.5.2 計算方法 運用C++語言編寫基于EULER ROMBERG 微分方程計算方法的優(yōu)化程序，應(yīng)用方程(2)、(3)、(5)和(10)進行臭氧溶解過程中氣泡的速度(u)、密度(ρg)和尺寸半徑(R)的數(shù)值計算。根據(jù)公式(14)～(16)計算氣泡上升的時間(t)和距離(S)。時間變量的步長取為0.0001 s，計算控制的相對誤差為≤0.0001。計算的初始數(shù)據(jù)和常系數(shù)如下:溫度T=293 K，氣泡表面的起始壓力P0=1.013 ×105Pa，氣泡的起始半徑R0=0.165 mm，氣泡流體的開始速度u0=0.05 m/s，水的運動黏度系數(shù)ν=1.003 ×106m2/s，臭氧擴散系數(shù)D=1.3 ×10-9m2/s，臭氧亨利常數(shù)H=4.6 ×108Pa·g/L，液體的表面張力σ=0.0728 N/m，水的動力黏性系數(shù)μf=1002 N·s/m3。計算流程圖如圖2所示。

圖2 臭氧氣泡溶解數(shù)值計算流程圖Fig.2 The flowing diagram of numerical calculation for ozone bubbles

1.6 計算結(jié)果

根據(jù)工廠化養(yǎng)殖循環(huán)量的需要，選擇系統(tǒng)水泵型號為50JYWQ25 -10，流量為25 t/h。由此計算得到反應(yīng)塔的主要尺寸:內(nèi)徑d=0.42 m，外徑Dw=0.70 m，高度S=3.00 m。

2 反應(yīng)塔應(yīng)用性能試驗

根據(jù)計算結(jié)果研制雙層逆流臭氧氧化反應(yīng)塔，并進行實際應(yīng)用效果試驗，以檢驗反應(yīng)塔氧化氨氮的能力和氧化過程中臭氧在水體中的殘留濃度。

2.1 材料與方法

2.1.1 試驗設(shè)備系統(tǒng) 試驗在冷水魚工廠化循環(huán)水養(yǎng)殖系統(tǒng)(圖3)中進行。系統(tǒng)工作過程中，從養(yǎng)魚池底部排出的養(yǎng)殖廢水經(jīng)過溢流槽后，經(jīng)微濾機過濾后進入回水槽，由水泵抽到鼓泡塔內(nèi)筒。同時，由臭氧發(fā)生器產(chǎn)生的臭氧也進入鼓泡塔內(nèi)筒，并在鼓泡塔內(nèi)外筒中混合、反應(yīng)，氧化水體中的氨氮，經(jīng)氧化處理過的水體返回養(yǎng)殖池，進行循環(huán)利用。試驗系統(tǒng)共有直徑1.8 m、高1.5 m 的圓形養(yǎng)殖池8 個，其中7 個設(shè)定為試驗池，1 個為對照池;采用XY-19 型臭氧發(fā)生器，額定氣體流量為8 m3/h，額定臭氧產(chǎn)量為100 g/h。

圖3 臭氧氧化試驗循環(huán)水養(yǎng)殖系統(tǒng)流程圖Fig.3 Flowing chart of a recirculating aquaculture system for ozone oxidation

2.1.2 材料 試驗魚為虹鱒Oncrhynchus mykiss，體質(zhì)量為(249 ±36)g，每個養(yǎng)殖池放養(yǎng)60 尾，約15 kg;按魚體質(zhì)量的1%每天9:00、14:00和19:00 投喂蛋白質(zhì)含量為45%的智利進口鮭鱒魚飼料(Salmofood);氧氣，純度不少于95%，流量為0.2 m3/h，臭氧產(chǎn)量為0.12 m3/h;臭氧氧化催化劑NaBr(分析純)，含量不少于99.0%;鼓泡塔內(nèi)水體為1.5 m3，水力停留時間為3.6 min;試驗池、沉淀池等整個循環(huán)系統(tǒng)水體為9.1 m3，對照池水體為1.3 m3。

2.1.3 方法 試驗在黑龍江水產(chǎn)研究所工廠化養(yǎng)殖試驗車間進行，試驗水溫為12 ℃。7 個試驗池在圖3所示的整個循環(huán)系統(tǒng)中進行試驗，對照池不參與系統(tǒng)循環(huán)，通過水泵在原池循環(huán)增氧。

(1)反應(yīng)塔臭氧催化氧化氨氮的效果。臭氧催化氧化氨氮由兩個化學(xué)反應(yīng)完成:一個是臭氧催化氧化的反應(yīng)，主要去除水中的離子銨NH+4，如反應(yīng)式(17)～(19)，另一個是臭氧直接氧化反應(yīng)，主要把水中的非離子氨NH3氧化為硝酸鹽，如反應(yīng)式(20)、(21)［7］。

NaBr 作為催化劑加入試驗水體，以Br/N 的比值為0.4 計算添加量，根據(jù)檢測的氨氮濃度，每天加入1 mg/L;每天8:00 第一次投喂餌料前，在鼓泡塔排水口即試驗池的入水口取2 個水樣、1 個平行樣;在對照池取1 個水樣、1 個平行樣。通過檢測養(yǎng)殖池和對照池水樣的氨氮、亞硝酸鹽、硝酸鹽，確定氧化氨氮的效果。分別采用納氏試劑比色法(GB7481 - 87)、酚二磺酸分光光度法(GB7480 -87)、N-(1 -萘基)-乙二胺分光光度法(GB7493 -87)測定氨氮、硝酸鹽、亞硝酸鹽的濃度，采用玻璃電極法(GB6920 -86)測定pH 值。

(2)養(yǎng)殖池水中臭氧濃度的控制。養(yǎng)殖池水中臭氧的安全濃度為0.01 mg/L，選用HACH 9185sc在線臭氧檢測儀檢測池水中臭氧濃度，并將檢測結(jié)果進行A/D 轉(zhuǎn)換后送給PLC 進行在線監(jiān)測和控制，調(diào)整臭氧發(fā)生器電壓，控制臭氧產(chǎn)量。在設(shè)定的臭氧發(fā)生器電壓下，檢測鼓泡塔內(nèi)筒水體的臭氧濃度。

2.2 結(jié)果與討論

2.2.1 氨氮氧化表1 為臭氧氧化過程中試驗組與對照組氨氮、硝酸鹽、亞硝酸鹽和pH 值的變化情況。由于試驗5 d 后對照組水質(zhì)惡化，魚類攝食明顯減少而停止水質(zhì)檢測。從表1可知:試驗組氨氮濃度與對照組相比，從第2天開始有明顯減少，第5天降低了54%，第5天臭氧催化氧化的效果比臭氧直接氧化(25.8%［6］)高1 倍多，表明產(chǎn)生了很好的催化氧化反應(yīng)，使水體中的大部分離子銨轉(zhuǎn)化成了氮氣，達到了去除氨氮的目的。在pH值為7 左右的環(huán)境中，大部分氨氮是以離子銨存在，在反應(yīng)過程中主要進行的是催化氧化反應(yīng)。試驗組硝酸鹽有積累現(xiàn)象，表明反應(yīng)過程也存在臭氧直接氧化的反應(yīng)，但是由于pH 值為7 左右，氨氮以非離子氨NH3存在的量非常少，直接氧化的效率很低。

試驗表明，試驗組仍存在氨氮積累的現(xiàn)象，證明臭氧還不能完全把氨氮氧化。主要原因:一是受到養(yǎng)殖水體pH 值中性的限制，不能充分發(fā)揮直接氧化的效能;二是受到養(yǎng)殖水體臭氧殘留濃度的限制，鼓泡塔內(nèi)臭氧濃度較低，反應(yīng)動力較小。試驗中亞硝酸鹽有所增加，但變化幅度較小，這是因為亞硝酸鹽是反應(yīng)的中間產(chǎn)物，殘留不會過高。試驗中pH 值呈現(xiàn)逐漸降低的趨勢，是加入催化劑NaBr后產(chǎn)生了溴酸HBrO 的原因，pH 值的降低影響了臭氧直接氧化的結(jié)果。

表1 臭氧氧化氨氮試驗的水質(zhì)變化Tab.1 The change in water quality for ammonia oxidation by ozone

NaBr 的加入和產(chǎn)生的HBrO可能在水體中產(chǎn)生溴酸鹽，從而對養(yǎng)殖魚類健康帶來不利影響。有研究表明，在Br/N 值選擇合理的條件下，Br 優(yōu)先反應(yīng)的是N，不會產(chǎn)生過多其他形式的溴化物［7］;同時，氨氮的存在可大幅減少HBrO 產(chǎn)生溴化物的可能性［18］，這也是盡管海水中存在大量Br 元素，臭氧仍然在海水養(yǎng)殖水處理中被廣泛應(yīng)用的原因［19-20］。

2.2.2 養(yǎng)殖池水中的臭氧濃度圖4 為應(yīng)用PLC在線監(jiān)測技術(shù)得到的養(yǎng)殖水體中臭氧濃度的變化情況。在檢測過程中，水體中溶解的臭氧濃度開始逐步上升，2 個多小時后基本穩(wěn)定在0.01 mg/L 附近，此時PLC 控制的臭氧發(fā)生器電壓為72 kV。這一結(jié)果表示，基于PLC 的養(yǎng)殖水體臭氧殘留濃度控制技術(shù)安全可靠，在魚類安全濃度0.008 ～0.060 mg/L［8］的范圍內(nèi)，不會對養(yǎng)殖魚類造成損傷。

圖4 5天內(nèi)養(yǎng)殖水體中臭氧殘留隨時間的變化Fig.4 Variation in dissolved ozone concentration in 5 days

2.2.3 鼓泡塔內(nèi)筒水體中的臭氧濃度 試驗表明，養(yǎng)殖池水體臭氧的安全濃度被穩(wěn)定控制時，臭氧發(fā)生器的電壓為72 kV，在該電壓條件下PLC在線檢測到的鼓泡塔內(nèi)筒水體臭氧濃度變化如圖5所示。內(nèi)筒水體的臭氧濃度在開始工作1 h 后逐漸上升，2 h 后臭氧濃度在0.09 mg/L 附近上下波動，波動趨勢比較穩(wěn)定。表明通過檢測養(yǎng)殖池水體的臭氧殘留濃度，進而用控制臭氧發(fā)生器電壓來控制臭氧產(chǎn)量的方法達到了設(shè)計要求;同時，與養(yǎng)殖池臭氧殘留濃度比較，內(nèi)筒水體臭氧濃度明顯高于安全濃度，水體中的臭氧在流經(jīng)鼓泡塔的過程中，89%的溶解臭氧與氨氮進行了氧化反應(yīng)，從而使水體到達養(yǎng)殖池時臭氧的濃度能夠符合安全濃度的要求。

圖5 鼓泡塔內(nèi)養(yǎng)殖水體臭氧濃度隨時間的變化Fig.5 Variation in dissolved ozone concentration in aquaculture water in reaction column water

3 結(jié)論

(1)根據(jù)氣泡運動、氣泡溶解和尺寸變化方程而設(shè)計的臭氧雙層反應(yīng)塔，設(shè)計合理、結(jié)構(gòu)簡單、體積小、重量輕，適應(yīng)于低溫下魚類養(yǎng)殖水體的氨氮處理。

(2)利用雙層臭氧反應(yīng)塔進行低溫氨氮催化氧化處理，可以達到去除水體中54%氨氮的效果，為冷水性魚類養(yǎng)殖水體的氨氮處理開辟了新途徑;

(3)由控制系統(tǒng)控制產(chǎn)生的臭氧和氨氮在反應(yīng)塔內(nèi)可以充分反應(yīng)，殘留在水體中的臭氧濃度不高于0.01 mg/L，符合魚安全類養(yǎng)殖的要求。

(4)由于低溫氨氮催化氧化處理還不能完全消除氨氮，實際應(yīng)用中應(yīng)結(jié)合低溫生物處理技術(shù)，進行分步處理，才能達到解決養(yǎng)殖過程氨氮積累的問題。

［1］Asgard T，Shearer K D.Dietary phosphorus requirement of juvenile Atlantic salmon，Salmo salar L.［J］.Aquacult Nutr，1997，3:17-23.

［2］Cho C Y，Hynes J D，Wood K R.Quantification of fish culture wastes by biological(nutritional)and chemical(limnological)methods［C］//University of Guelph.Proceedings of the first international symposium on nutritional strategies in management of aquaculture waste.Ontario:University of Guelph，1991:37 -50.

［3］Thomas M，Losordo，Michael P.Recirculating aquaculture tank production systems(an overview of critical considerations)［J］.SRAC Publication，1998，451:1 -6.

［4］曹廣斌，蔣樹義，韓世成.雙層浮球生物濾器設(shè)計及其水產(chǎn)養(yǎng)殖水處理性能試驗［J］.水產(chǎn)學(xué)報，2005，29(4):578 -582.

［5］Trapido M，Veressinina Y，Munter R.Ozonation advanced oxidation process of polycyclic aromatic hydrocarbons in aqueous solutions -a kinetic study［J］.Environmental Technol，1995，16(8):729-740.

［6］Lin S H，Wu C L.Removal of nitrogenous compounds from aqueous solution by ozone and ion exchange［J］.Wat Res，1996，30(8):1851 -1857.

［7］Yang M，Uesugi K，Myoga H.Ammonia removal in bubble column by ozonation in the present of bromide［J］.Wat Res，1999，33(8):1911 -1917.

［8］Bullock G L，Summerfelt S T，Noble A.Ozonation of a recirculating rainbow trout culture system:I.effects on bacterial gill disease and heterotrophic bacteria［J］.Aquaculture，1997，158:43 -55.

［9］馬躍龍，黃娟，沈春銀，等.鼓泡塔中液相混合時間的影響因素［J］.華東理工大學(xué)學(xué)報，2010，36(2):165 -172.

［10］曹廣斌，賈慧文，蔣樹義，等.循環(huán)水養(yǎng)魚系統(tǒng)中臭氧射流混合設(shè)備設(shè)計與性能測試［J］.農(nóng)業(yè)工程學(xué)報，2011，27(10):73 -78.

［11］Takemura F，Yabe A.Gas dissolution process of spherical rising gas bubbles［J］.Chemical Engineering Science，1998，53(15):2691 -2699.

［12］郭烈錦.兩相與多相流體動力學(xué)［M］.西安:西安交通大學(xué)出版社，2002:427 -449.

［13］Farines V，Baig S，Albet J，et al.Ozone transfer from gas to water in a co - current up - flow packed bed reactor containing silica gel［J］.Chemical Engineering Journal，2003，91:67 -73.

［14］Kuosa M，Laari A，Kallas J.Determination of the Henry's coefficient and mass transfer for ozone in a bubble column at different pH values of water［J］.Ozone Science and Engineering，2004，26:277 -286.

［15］Godwin R P，Chapyak E J.Resource letter BD-1:nubble dynamics［J］.Am J Phys，2000，68(3):211 -219.

［16］The Chemical Society of Japan.Handbook of Chemistry［M］.4th ed.Tokyo:Iwanami Publisher，1993:58 -63.

［17］張金松，李偉光，馬放，等.臭氧接觸裝置的傳質(zhì)與吸收試驗研究［J］.哈爾濱建筑大學(xué)學(xué)報，1997，30(2):75 -79.

［18］張書芬，王全林，沈堅，等.飲用水中臭氧消毒副產(chǎn)物溴酸鹽含量的控制技術(shù)探討［J］.水處理技術(shù)，2011，37(1):28 -32.

［19］宋奔奔，倪琦，張宇雷，等.臭氧對大菱鲆半封閉循環(huán)水養(yǎng)殖系統(tǒng)水質(zhì)凈化研究［J］.漁業(yè)現(xiàn)代化，2011，38(6):11 -15.

［20］周游，黃濱，宋協(xié)法，等.臭氧前后置工藝變化對循環(huán)水半滑舌鰨養(yǎng)殖系統(tǒng)水環(huán)境的影響［J］.漁業(yè)現(xiàn)代化，2012，39(6):5 -9.