鄒德兵，王一帆，丁 剛

（1.長江勘測規(guī)劃設(shè)計研究有限責(zé)任公司，湖北 武漢 430010；2.武漢大學(xué)水資源工程與調(diào)度全國重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430072）

0 引言

已有研究指出微氣泡直徑介于1～1 000 μm 或亞微米尺度的氣泡，其物理化學(xué)特性在污水處理、生物醫(yī)藥、船舶減阻等領(lǐng)域受到了廣泛關(guān)注［1］。當(dāng)前，城市供水項(xiàng)目的興建使有壓輸水管道日益增多，而微氣泡的數(shù)量和尺寸會對水廠水處理工藝和輸水管道的安全運(yùn)行產(chǎn)生顯著影響［2］：均勻的直徑小于50 μm的微氣泡可以增加曝氣效率［3］，而較大的微氣泡則容易在管道中聚集形成大型氣泡，增加管道水頭損失，降低通水能力，危害輸水管道運(yùn)行安全，降低水廠設(shè)備工作效率。這使得對有壓輸水管道中的微氣泡進(jìn)行觀測和控制的相關(guān)研究成為必要。曹命凱等［4］研究了長距離有壓輸水管道充水過程中排氣不盡的影響，認(rèn)為管道充水時排氣不盡或管道入口摻混氣體，會影響輸水效率并增加爆管風(fēng)險，可以通過控制充水方式和設(shè)置排氣閥方式排出；楊誠志［5］則對氣泡運(yùn)動規(guī)律進(jìn)行研究，總結(jié)了氣泡起動流速與氣泡大小的量化關(guān)系。

當(dāng)前對水流中微氣泡進(jìn)行觀測的方式主要有激光成像法［6］、PIV 法［7］、探針法［8］、高速攝影圖像法等［9-12］。其中高速攝影圖像法因其測量的非接觸性以及設(shè)備維護(hù)簡單，適合有壓輸水管道模型試驗(yàn)的水流微氣泡測量。現(xiàn)有相關(guān)研究往往偏于理論研究，一般利用特定試驗(yàn)管道模型進(jìn)行測試與分析，并未與實(shí)際工程運(yùn)行直接關(guān)聯(lián)。本文以深圳市北坑水庫配套輸水工程為例，基于高速攝影圖像法和水工模型實(shí)驗(yàn)對有壓輸水管內(nèi)微氣泡在消力池消能效果、輸水流量變化時的分布特性進(jìn)行了研究，并提出了合理的有壓輸水管內(nèi)微氣泡控制措施。

1 微氣泡分布特性試驗(yàn)方案

1.1 工程概況

北坑水庫壩址位于深汕特別合作區(qū)赤石河上游峽谷出口，是深圳市水源規(guī)劃布局中“三蓄”水庫之一。其配套輸水工程由輸水隧洞段、隧洞出口消力池與有壓輸水埋管段構(gòu)成，上游起點(diǎn)為北坑水庫，下游終點(diǎn)為中心水廠。輸水隧洞段長8.93 km，隧洞洞徑2.4 m，隧洞出口中心線高程56.6 m；消力池控制水位60.3 m；消力池出口接有壓輸水埋管，該處管道中心線高程53.5 m；有壓輸水埋管段長471 m，管道直徑1.8 m，暫定埋管出口中心線高程56.0 m。工程輸水方式為全程有壓重力流輸水，設(shè)計引水流量2.5～7.64 m3/s。

1.2 物理模型設(shè)計

一般來說，正態(tài)模型設(shè)計需滿足重力相似準(zhǔn)則，且保證水流位于完全紊流區(qū)（雷諾數(shù)Re≥1 000～2 000），水深大于2 cm。在研究微氣泡的氣液兩相流正態(tài)模型試驗(yàn)中，除了滿足正態(tài)模型的各種相似條件之外，還需在調(diào)節(jié)試驗(yàn)工況過程中，保證管道壓力緩慢變化。

在模型與原型壓強(qiáng)相似的情況下，模型中的氣泡由于本身尺寸更小，受到壓強(qiáng)影響相對更小，從而使得氣泡直徑在部分長度比尺下并不滿足幾何相似條件。根據(jù)郭永鑫、楊開林等的研究成果［14，15］，使用正態(tài)模型進(jìn)行管道氣泡相關(guān)試驗(yàn)時，考慮壓強(qiáng)影響后氣泡在原型與模型中直徑比值λq滿足式（1）：

式中：λl為長度比尺；Pa、Pm分別為標(biāo)準(zhǔn)大氣壓、任意時刻模型氣泡所在位置水壓；

整體模型包括水循環(huán)系統(tǒng)、上游水庫、進(jìn)水塔、輸水隧洞段、消力池、有壓輸水埋管段、下游水廠。其中，輸水隧洞段模型通過阻力環(huán)替代部分輸水隧洞，有壓輸水埋管段模型簡化由進(jìn)口水平段、轉(zhuǎn)彎及下彎段、出口水平段、出口段組成。進(jìn)水塔、輸水隧洞段、消力池及有壓輸水埋管段采用透明有機(jī)玻璃制作，以滿足管道糙率相似要求及觀測需要。模型完成后的模型過流能力測試結(jié)果表明，模型的泄流能力與水頭損失與原型相似，能滿足設(shè)計要求。模型布置如圖1所示。

圖1 物理模型布置圖（單位：m）Fig.1 Layout of the design of the model test

1.3 高速攝影裝置及測試原理

模型試驗(yàn)采用非接觸式的高速攝影圖像法研究氣泡在有壓輸水埋管內(nèi)的運(yùn)動過程（包括微氣泡大小、數(shù)量等）。將非頻閃LED強(qiáng)光源置于有壓輸水埋管下方，采用VEO高速攝影機(jī)對典型管段進(jìn)行拍攝，利用VEO 攝影機(jī)自帶的PCC圖像采集軟件對拍攝圖像進(jìn)行采集和初步切分，接著對切分后的圖像進(jìn)行灰度處理，再進(jìn)一步通過調(diào)整霍夫變換算法中圓心的累加器閾值以及邊緣梯度，可以查找不同尺寸微氣泡并計算出其所占像素。最后檢測每幀微氣泡像素與總像素占比以及每幀氣泡數(shù)量，并依據(jù)圖像中標(biāo)尺長度計算采樣塊面積，從而計算模型中氣泡實(shí)際大小。通過統(tǒng)計模塊給出上述測量結(jié)果的統(tǒng)計參數(shù)，即可得到該管段的微氣泡數(shù)量、直徑等數(shù)據(jù)。

試驗(yàn)中采用的高速攝影裝置示意圖如圖2所示。

圖2 高速攝影測試裝置示意圖Fig.2 Schematic of high-speed photography device

1.4 試驗(yàn)工況

經(jīng)過前期多次調(diào)試，試驗(yàn)最終選擇典型采樣管段圖像寬度（沿管軸線）為130 mm（模型值），VEO 高速攝影機(jī)拍攝參數(shù)如下：圖像分辨率1 280×800、攝像速度400 fps，曝光時間1 600 μS，分析采用圖像共2 400幀（6 s）。

通過對輸水埋管運(yùn)行工作條件的分析，選取如表1所示的4個代表性工況進(jìn)行試驗(yàn)。

表1 模型試驗(yàn)典型工況及控制參數(shù)Tab.1 Typical working conditions and control parameters of model test

2 微氣泡分布特性試驗(yàn)成果

2.1 微氣泡分布特性

模型試驗(yàn)中，在各工況下調(diào)試上下游水位及流量，保證出水口中心水廠水位穩(wěn)定，使得輸水埋管段穩(wěn)定運(yùn)行工作，進(jìn)而對輸水埋管段微氣泡分布特性進(jìn)行觀測。各工況進(jìn)口水平段圖像對比如圖3所示。

圖3 試驗(yàn)工況下進(jìn)口水平段圖像對比Fig.3 Comparison of images under test conditions

由圖3 可知，試驗(yàn)工況下輸水埋管進(jìn)口水平段內(nèi)氣泡數(shù)量很少，難見直徑大于1 mm 的氣泡，僅工況3［如圖3（c）所示］在圖像中發(fā)現(xiàn)了少數(shù)直徑較大的氣泡。采用圖像分析軟件進(jìn)行圖像分析后，得到各工況在拍攝歷時內(nèi)觀測到的微氣泡大小與總數(shù)量對應(yīng)關(guān)系如圖4所示，微氣泡分布特性關(guān)鍵參數(shù)統(tǒng)計如表2所示。

表2 各工況微氣泡分布特性統(tǒng)計Tab.2 Statistics of the distribution characteristics of microbubbles

圖4 試驗(yàn)工況下微氣泡分析結(jié)果Fig.4 Analysis results of microbubbles under test conditions

通過圖4 和表2 可以發(fā)現(xiàn)，隨著測試工況發(fā)生變化，微氣泡的分布特性發(fā)生了顯著改變，工況1、工況2、工況4 微氣泡面積基本在0.05 mm2以下，工況3 下觀測到的微氣泡面積多大于0.05 mm2。每幀最多同時存在氣泡數(shù)量和氣泡總數(shù)量也隨著測試工況變化也發(fā)生了較大改變。需進(jìn)一步分析輸水流量和與消力池消能率對微氣泡分布特性的影響。

2.2 消力池消能率對微氣泡分布特性的影響

加入示蹤粒子（示蹤劑）后的工況2 與工況3 消力池流態(tài)如圖5所示。

圖5 消力池流態(tài)Fig.5 Flow pattern of stilling basin

通過圖5（a）、（b）可以發(fā)現(xiàn)，兩工況下示蹤粒子在消力池水體表面分布較均勻，沒有發(fā)現(xiàn)表面渦旋。通過圖5（c）、（d）可以發(fā)現(xiàn)，兩工況下消力池內(nèi)均未出現(xiàn)吸氣漩渦，不會導(dǎo)致氣體從消力池表面被吸入輸水埋管。這一現(xiàn)象也表明，有壓輸水埋管內(nèi)微氣泡主要來源于上游水體進(jìn)入消力池消能后所形成的殘余氣泡。

通過圖5（c）、（d）中示蹤劑分布可以發(fā)現(xiàn)，消力池消能過程中，主流從輸水隧洞進(jìn)入消力池后直接沖擊消力池內(nèi)隔板后，分別從隔板過流孔和隔板上緣進(jìn)入消力池后室再進(jìn)入輸水埋管。工況2下主流明顯被隔板分為兩股進(jìn)入隔板后消力池后室發(fā)生摻混；工況3 下隔板的分流效果不如工況2 明顯，主流很少穿過隔板過流孔，而是翻越隔板上緣進(jìn)入消力池后室。同時，工況2 下表面示蹤粒子更加集中分布于靠近下游側(cè)，而工況3下表面示蹤粒子則分布相對更均勻，計算兩種工況消能率發(fā)現(xiàn)，工況2 消能率為95.69%，工況3 消能率僅為87.01%?？梢哉J(rèn)為，工況2中由輸水隧洞進(jìn)入消力池水流能量更大，消能過程中水流擾動更強(qiáng)，而工況3 下消力池內(nèi)流態(tài)平穩(wěn)，水流擾動很小，主流中包含的微氣泡等更易集中進(jìn)入消力池下游輸水埋管。

通過以上分析以及圖4 和表2 數(shù)據(jù)可知：在流量不變時上游庫水位更低的工況，水流在消力池內(nèi)受到擾動越強(qiáng)，上游水體所含氣泡在消力池內(nèi)更易散成體積更小的微氣泡，進(jìn)而使得進(jìn)入輸水埋管的水流中微氣泡分布密度更低，尺寸更小。

2.3 輸水流量對微氣泡分布特性的影響

由圖4（a）、（b）、（d）可知，增大輸水流量使得管道內(nèi)微氣泡直徑減小的同時，拍攝歷時內(nèi)捕捉到的微氣泡總數(shù)量先增加后減少。根據(jù)表2數(shù)據(jù)可知，隨著輸水流量的增大，微氣泡占采樣塊面積最大百分比逐漸增加，說明微氣泡分布密度更大。

在有壓輸水埋管進(jìn)口測得工況1、工況2、工況4 管道中心線流速分別為0.80、1.63、2.35 m/s。分析認(rèn)為小流量時管內(nèi)流速很小，使得管內(nèi)氣泡難以被水流帶出或隨水流運(yùn)動很慢，使得觀測到的氣泡直徑較大、數(shù)量較少；流量增大時，進(jìn)入管道的微氣泡數(shù)量更多，但管內(nèi)流速也逐漸增大，使得微氣泡不再能附著在管壁，聚集過程也縮短，觀測到的微氣泡直徑減小；流量進(jìn)一步增大時，流速增大使得微氣泡移動更快，觀測到的微氣泡尺寸變化不明顯，但分布密度變小。

可以認(rèn)為適當(dāng)增大管道輸水流量可以有效減少管內(nèi)微氣泡數(shù)量和分布密度。對于本工程而言，在上游水位為98 m 時，輸水流量從2.50 m3/s 增加至5.35 m3/s 時，管道內(nèi)微氣泡直徑減小，分布密度增加；輸水流量從5.35 m3/s 繼續(xù)增加至7.64 m3/s時，微氣泡直徑變化不明顯，但分布密度更小。

3 微氣泡控制措施試驗(yàn)成果

以微氣泡分布特性試驗(yàn)成果為基礎(chǔ)，從充水排氣方式、排氣閥設(shè)置及改變管道直徑等方面進(jìn)行了管內(nèi)微氣泡的控制措施試驗(yàn)。

3.1 充水排氣方式

根據(jù)《給排水設(shè)計手冊》［16］，在管道系統(tǒng)初次啟用或檢修重啟時，需要遵循小流量充水原則對管道進(jìn)行充水排氣。即滿足“盡可能地將充水時的滿管平均流速規(guī)定為0.3 m/s，最大不超過0.6 m/s”的情況下，同時也要滿足流速大于不淤流速0.5 m/s，（實(shí)際采用流速原型值為0.6 m/s）。為驗(yàn)證小流量預(yù)充水效果，選取微氣泡特性試驗(yàn)中微氣泡尺寸和分布密度較大的工況3進(jìn)行對比，對不同排氣方式下（通水流量直接充水和小流量預(yù)充水）管道內(nèi)微氣泡狀態(tài)進(jìn)行觀測。

在上游水庫、消力池、下游水廠水位均穩(wěn)定后，工況3 運(yùn)行條件下，不同排氣方式下管道內(nèi)微氣泡堆積狀態(tài)如圖6所示。

圖6 不同排氣方式下管內(nèi)微氣泡堆積情況Fig.6 Microbubble accumulation with different exhaust methods

充水排氣試驗(yàn)結(jié)果表明，采取通水流量直接充水會產(chǎn)生明顯的微氣泡堆積，雖未形成大型氣泡或氣囊，但對有壓輸水管道的微氣泡控制不利。采取小流量預(yù)充水后，管壁沒有觀察到明顯微氣泡附著和堆積，微氣泡控制效果良好。

3.2 設(shè)置排氣閥

為研究有壓輸水管道沿程微氣泡狀態(tài)，驗(yàn)證沿程排氣閥對微氣泡控制的有效性，用模型排氣孔對排氣閥進(jìn)行模擬，根據(jù)相關(guān)規(guī)程及文獻(xiàn)［17，18］，考慮管道起伏轉(zhuǎn)彎情況，在有壓輸水埋管段進(jìn)口水平段、轉(zhuǎn)彎及下彎段彎道前后、出口水平段、出口段控制閥前共設(shè)置排氣孔5 處。工況3 管段沿程氣泡狀態(tài)對比如圖7所示。對圖8（a）、（c）、（d）的典型測試部位進(jìn)行微氣分布特性統(tǒng)計，結(jié)果如表4所示。

圖7 管道沿程氣泡堆積狀態(tài)Fig.7 Microbubble accumulation along the pipe

圖8 不同管徑微氣泡圖像對比Fig.8 High-speed photographic images of microbubbles in different diameter pipes

由圖8 及表3 可知，在沿程設(shè)置排氣閥后，水流中微氣泡尺寸和數(shù)量均呈現(xiàn)比較明顯的沿程下降趨勢，微氣泡占管道面積總百分比也明顯下降。在水廠內(nèi)的輸水管道出口，最大微氣泡面積僅0.006 mm2（模型值），說明沿程設(shè)置排氣閥可以有效排出直徑較大的微氣泡，對有壓輸水管道內(nèi)微氣泡數(shù)量及尺寸控制有利。

表3 設(shè)置排氣閥后各管段微氣泡分布特性關(guān)鍵參數(shù)Tab.3 Parameters for the distribution of microbubbles in pipeline with exhaust valves

3.3 改變管道直徑

為研究管徑變化下有壓輸水埋管內(nèi)水流微氣泡的分布特性，將2.4 m 與1.8 m 管徑下管內(nèi)微氣泡分布特性試驗(yàn)數(shù)據(jù)進(jìn)行了對比。在采取小流量預(yù)充水直至設(shè)計通水流量后，工況3 進(jìn)口水平段不同管徑下微氣泡觀測結(jié)果對比如圖8所示。微氣泡分布特性關(guān)鍵統(tǒng)計參數(shù)對比如表3所示。

通過表4 的統(tǒng)計參數(shù)可知，當(dāng)管道直徑為1.8 m 時，有壓輸水管內(nèi)微氣泡直徑和分布密度均明顯低于2.4 m 直徑管道。在工況3 試驗(yàn)條件下，2.4 m 直徑管內(nèi)流速為1.18 m/s，小于1.8 m直徑管內(nèi)流速1.63 m/s，說明縮小管徑提高流速使得微氣泡在未堆積形成尺寸較大氣泡或氣囊時就被水流帶至下游，使得管內(nèi)微氣泡尺寸和分布密度減小。

表4 不同管徑微氣泡分布特性統(tǒng)計Tab.4 Statistical characteristics of microbubbles in different diameter pipes

上述分析表明，適當(dāng)縮小管徑提高管道流速有利于對氣泡尺寸和分布密度進(jìn)行控制。但由于管道結(jié)構(gòu)和過流能力的限制，輸水管道直徑不能過小，需保證管內(nèi)流速不大于管道抗沖流速，且滿足過流能力要求。

4 結(jié)論

以北坑水庫輸水配套輸水工程水工模型試驗(yàn)為依托，采用高速攝影圖像法，對有壓輸水管道內(nèi)微氣泡分布特性及微氣泡控制措施進(jìn)行了研究，研究結(jié)果表明：

（1）在輸水隧洞段出口消力池內(nèi)未觀測到表面渦旋和吸氣漩渦，有壓輸水埋管內(nèi)微氣泡主要來源于上游水體進(jìn)入消力池消能后形成的殘余氣泡。在消能率更高的工況下，較大氣泡在消力池內(nèi)更易受到擾動被打散成更小的微氣泡，使得進(jìn)入下游埋管的微氣泡直徑更小，分布密度也更低。

（2）在上游水位98 m，輸水流量從2.50 m3/s增加至5.35 m3/s時，管道內(nèi)微氣泡直徑減小，分布密度增加；流量進(jìn)一步增加至7.64 m3/s時，微氣泡直徑變化不明顯，但分布密度下降。

（3）管道初次啟用或檢修后重啟時保證小流量充水，沿程合理設(shè)置排氣閥，在滿足防沖要求和過流能力要求的前提下適當(dāng)減小管徑，均可對有壓輸水管內(nèi)微氣泡直徑和分布密度進(jìn)行有效控制。