(1.天津大學 水利工程仿真與安全國家重點實驗室，天津 300072；2.華能瀾滄江水電股份有限公司，昆明 650214)

1 研究背景

我國水能資源豐富，理論蘊藏量6.94億kW，技術可開發(fā)量5.42億kW。其中西部高海拔地區(qū)理論蘊藏量2.055億kW，位居全國第一，技術可開發(fā)量1.16億kW，是未來水電開發(fā)的重點。但西部高海拔地區(qū)具有氣壓低、氧氣稀薄的特點，平均海拔超過4 000 m，平均氣壓在65.25 kPa以下，不足海平面的2/3，空氣密度0.57～0.89 kg/m3，含氧量只有平原地區(qū)的40%～50%。高壩泄洪的高速水流涉及水氣兩相流復雜機理，會受到環(huán)境氣壓的影響。

國外水電工程海拔普遍較低，除空泡動力學機理性研究外，幾乎沒有高海拔水電工程水力學特性的研究成果。國內(nèi)已建工程海拔大多在2 000 m以下，少量高海拔地區(qū)的水電工程在設計時對環(huán)境氣壓的考慮簡單粗糙。Plesset[1]在考慮了空泡表面張力的情形下，用 Rayleigh方法計算了外壓變化時蒸汽泡直徑的變化。吳建華[2]分析水中的空氣含量隨周圍壓力和溫度的變化，當壓力降低時，空氣含量將減少，如果原型和模型含氣量不同，擴散空化將不同，所得的初生空化數(shù)也將不同。朱宗銘等[3]通過研究發(fā)現(xiàn)環(huán)境壓力越小，水輪機沖蝕和空蝕交互磨損程度越嚴重。因此，環(huán)境氣壓對高速水流的流態(tài)空化等特性具有一定程度的影響[4-8]。楊吉健等[9]指出處于高海拔地區(qū)的工程必須考慮空蝕破壞的可能性。夏葉青等[10]在設計西藏拉洛水利樞紐多功能泄洪隧洞時，將洞身設計為有壓隧洞并采用圓形斷面，確?？紤]高原地區(qū)低壓的影響后洞身壓力仍滿足要求。蒙富強等[11]進行藏木水電站泄洪閘水工模型試驗時明渠水流空化數(shù)為0.22～0.31，考慮到海拔較高，預計極有可能發(fā)生空化空蝕破壞，建議采取合理的摻氣措施保護底板。

理論上，環(huán)境氣壓會影響水體中氣泡的分布和尺寸，從而影響水流空化空蝕及其他相關特性。但有些水電工作者認為環(huán)境氣壓對高速水流的水力學特性影響微乎其微，不需要特別關注。因此，考慮環(huán)境氣壓因素，定量研究其對高海拔地區(qū)高速水流空化空蝕特性的影響程度，對西部高海拔地區(qū)擬建大批水電工程的摻氣減蝕設計具有重要的參考價值。

2 環(huán)境氣壓對初生空化數(shù)的影響

2.1 單氣核臨界壓強受環(huán)境氣壓影響程度

空化與空蝕的根本原因是氣核在低壓區(qū)生長和失穩(wěn)轉(zhuǎn)化成空泡后在高壓區(qū)潰滅，沒有氣核就不會有空化。因此空化發(fā)生的準則是：空化核由穩(wěn)定狀態(tài)生長到慣性失穩(wěn)狀態(tài)。在常溫且溫度不變的條件下，流場中最小瞬時壓強pmin小于給定氣核相應的臨界壓強pc，且該氣核在壓強為pmin的區(qū)域中運動至某一點時，它的半徑恰好生長到臨界半徑Rc，則給定的氣核在該點失穩(wěn)而形成空泡[12]。臨界壓強是水流空化發(fā)生的關鍵參數(shù)。

W表示氣核的韋伯數(shù)，則有

W=2S/[R0(p0-pv)] 。

(1)

式中：S為液體的表面張力系數(shù)；R0為氣泡的初始半徑；p0為泡外的液體壓強；pv為飽和蒸汽壓強。

Rc為氣泡臨界半徑，則有

(2)

pc為臨界壓強，則有

(3)

假設泡內(nèi)的蒸汽壓強pv=2 500 N/m2，水下1 m深處平衡的氣核半徑R0=10-5m，水的表面張力系數(shù)S=0.075 N/m，則氣核失穩(wěn)的臨界壓強隨大氣壓變化如表1。由表1可知，在水壓一定的情況下，氣泡失穩(wěn)的臨界壓強隨著大氣壓的減小而減?。缓０蚊可仙?00 m，氣壓每降低5 600 Pa，氣泡失穩(wěn)的臨界壓強降低9%～16%。因此在水深1 m處，半徑10-5m的穩(wěn)定氣核，泡內(nèi)蒸汽壓強2 500 N/m2的假設條件下，大氣壓對單個氣泡失穩(wěn)的臨界壓強影響較大，大氣壓越低，單個氣泡臨界壓強越小，越不容易發(fā)展為空泡。

表1 大氣壓對氣核失穩(wěn)臨界壓強的影響Table 1 Effect of atmospheric pressure on criticalpressure of gas nucleus instability

圖1 不同水深處氣泡失穩(wěn)的臨界壓強隨大氣壓的變化趨勢Fig.1 Variation tendency of critical pressure of bubble instability in different depths of water with atmospheric pressure

圖1為不同水深處氣泡失穩(wěn)的臨界壓強隨大氣壓的變化趨勢。假設位于1,10,20 m處的同等大小的穩(wěn)定氣核，當大氣壓變化時，氣核失穩(wěn)的臨界壓強也隨之變化；且水深越深，大氣壓對氣核失穩(wěn)臨界壓強的影響程度越小，水深10 m處的穩(wěn)定氣核，海拔每上升500 m，氣壓每降低5 600 Pa，氣泡失穩(wěn)的臨界壓強降低2%左右；水深20 m處的穩(wěn)定氣核，海拔每上升500 m，氣壓每降低5 600 Pa，氣泡失穩(wěn)的臨界壓強降低1%左右。水深越深，同等大小氣核的失穩(wěn)臨界壓強越大，氣核越容易發(fā)生失穩(wěn)。

2.2 單氣核臨界壓強變化對初生空化數(shù)的影響

初生空化數(shù)σi計算式為

(4)

式中：p∞,i，v∞,i分別為剛剛察覺到水流中有空化發(fā)生時參考點的壓強和相應的斷面平均流速；ρ為水的密度?？栈跎桥c氣核失穩(wěn)聯(lián)系在一起的，氣核一旦失穩(wěn)，立即開始“爆發(fā)性”生長，同時氣核變成空泡。初生空化數(shù)的實測過程表明，當參考點的壓強和流速分別為p∞,i，v∞,i時，系統(tǒng)中空化初生部位的瞬時壓強pmin已低于與水溫相應的飽和蒸汽壓強pv。因此單氣核失穩(wěn)時的臨界壓強與宏觀空化現(xiàn)象發(fā)生時的參考點壓強之間存在著量級之間的關系。因此假設海拔0 m處流速為30 m/s的水流初生空化數(shù)為0.2，根據(jù)式(4)推算出空化初生時參考點壓強p∞,i為92 500 Pa，而此條件下單氣核失穩(wěn)臨界壓強pc為506.37 Pa。

根據(jù)上述pc與p∞,i之間的量級關系，若海拔0 m時30 m/s流速的水流初生空化數(shù)為0.2，則海拔每升高500 m，初生空化數(shù)減小程度約0.05%～0.09%；若海拔0 m時30 m/s流速的水流初生空化數(shù)為2，則海拔每升高500 m，初生空化數(shù)減小程度約0.005%～0.009%，見表2。由此可知，若水流初生空化數(shù)越小，環(huán)境氣壓對初生空化數(shù)的影響程度越明顯；海拔對初生空化數(shù)的影響程度很小，海拔每升高500 m，初生空化數(shù)減小程度基本都在<0.1%的量級水平。也就是說，環(huán)境氣壓對單個氣核的失穩(wěn)有較大影響，但單個氣核的失穩(wěn)臨界壓強變化程度對于整個系統(tǒng)空化初生現(xiàn)象的影響已極小。

表2 海拔對初生空化數(shù)的影響程度Table 2 Influence of altitude on incipient cavitation number

3 環(huán)境氣壓對水流空化數(shù)的影響

氣核的存在是形成空化的內(nèi)在因素，壓強降低是形成空化的外在動力。工程上，采用空化數(shù)作為衡量實際水流是否發(fā)生空化的指標，空化數(shù)是一個無量綱數(shù)，以K表示。

(5)

式中：p為邊界的絕對壓強；v為相應的斷面平均流速；pv為當?shù)氐恼羝麎簭?；hw為相應斷面水流邊界壓強水頭；ha為相應斷面大氣壓強水頭；hv為相應水溫下的水的汽化壓強水頭。

假設hw=10 m，水的溫度20 ℃，hv=0.233 8 m，斷面水流流速為30 m/s。根據(jù)式(5)得到水流空化數(shù)隨海拔高度的變化，見圖2。在水流壓強、水流流速、水體溫度均不變的情況下，海拔越高，氣壓越低，水流空化數(shù)越??；海拔每升高500 m，氣壓每下降5 600 Pa，水流空化數(shù)減小2.5%～3.5%的量級水平。

圖2 空化數(shù)K隨海拔的變化規(guī)律(v=30 m/s)Fig.2 Variation regularity of water cavitation numberK with elevation(v=30 m/s)

高速水流發(fā)生空化的條件即水流空化數(shù)小于初生空化數(shù)(K

4 結(jié) 論

作為未來水電開發(fā)熱點的西藏地區(qū)具有高海拔、低氣壓的環(huán)境特點，而目前國內(nèi)外泄洪消能技術并未系統(tǒng)考慮大氣壓對高速水流水力學特性的影響。本文首次考慮環(huán)境氣壓因素影響，采用理論分析方法，分析論證環(huán)境氣壓對高速水流空化特性的影響。主要結(jié)論如下：

(1)大氣壓對高速水流空化特性是有影響的。高海拔地區(qū)的低氣壓環(huán)境會導致水流初生空化數(shù)和水流空化數(shù)同時減小，但水流空化數(shù)減小的程度遠遠大于初生空化數(shù)減小程度，因此高海拔地區(qū)的高速水流空化空蝕風險增大。

(2)建議高海拔地區(qū)泄洪系統(tǒng)水力設計考慮低氣壓環(huán)境的影響，海拔每升高500 m，水流空化數(shù)宜較常壓計算值減小3%左右的水平進行空化空蝕風險評價。

(3)高速水流空化特性涉及氣核與水體之間復雜的微觀的作用機理，本文主要分析氣壓因素的影響，因此假設在不同環(huán)境氣壓、不同水深條件下穩(wěn)定存在的氣核尺寸是不變的，且氣核的失穩(wěn)臨界壓強與空化初生的參考點壓強之間存在量級關系。

需指出：文中的案例數(shù)據(jù)是為分析初生空化數(shù)及水流空化數(shù)隨海拔變化的影響程度而設定，結(jié)論可供設計參考，具體工程問題還需具體分析。由于目前還未有高海拔地區(qū)高壩泄洪原型觀測資料，減壓模型試驗技術具有縮尺效應，本文成果還有待更多案例尤其是原型觀測結(jié)果論證。