洪屏抽蓄電站地下廠房對(duì)不同頻域壓力脈動(dòng)的振動(dòng)響應(yīng)

張 軍，伍鶴皋 ，劉 建，石長(zhǎng)征

(1.江西洪屏抽水蓄能有限公司，江西 靖安 330603；2.武漢大學(xué)水資源與水電工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，武漢 430072)

洪屏抽水蓄能電站一期裝機(jī)容量4×300 MW，機(jī)組正常轉(zhuǎn)速500 r/min，對(duì)應(yīng)轉(zhuǎn)頻8.33 Hz，轉(zhuǎn)輪葉片數(shù)為9，電站額定水頭高達(dá)540 m，具有高轉(zhuǎn)速高水頭的特征。水輪機(jī)供應(yīng)商提供了水泵水輪機(jī)模型試驗(yàn)報(bào)告，其中包含540 m水頭附近開(kāi)度γ從6°～26°，相鄰測(cè)量工況間隔2°情況下的脈動(dòng)壓力實(shí)測(cè)資料。目前，國(guó)內(nèi)外抽水蓄能電站水泵水輪機(jī)組脈動(dòng)壓力模型驗(yàn)收試驗(yàn)主要是根據(jù)IEC60193的相關(guān)標(biāo)準(zhǔn)進(jìn)行驗(yàn)收，一般是規(guī)定97%置信度的實(shí)測(cè)壓力脈動(dòng)峰-峰雙振幅值不超過(guò)保證值[1-4]，而不考慮脈動(dòng)壓力的頻率特性。目前涉及尾水管內(nèi)壓力脈動(dòng)的絕大部分文獻(xiàn)認(rèn)為尾水管內(nèi)的壓力脈動(dòng)以低頻為主[3,5,6]，對(duì)尾水管壓力脈動(dòng)進(jìn)行頻譜分析時(shí)將截?cái)囝l率取為5倍的額定轉(zhuǎn)頻甚至更小[7]，還有文獻(xiàn)提到對(duì)尾水管壓力脈動(dòng)進(jìn)行低通濾波[5]。但從尾水管脈動(dòng)壓力實(shí)測(cè)資料來(lái)看，其壓力脈動(dòng)的頻率成分是比較豐富的[7-9]。本文開(kāi)篇提到水輪機(jī)供應(yīng)商提供的模型實(shí)測(cè)數(shù)據(jù)也表明尾水管內(nèi)壓力脈動(dòng)的頻率成分是比較豐富的。本文依據(jù)壓力脈動(dòng)實(shí)測(cè)資料，分析其頻譜特性，發(fā)現(xiàn)尾水管中出現(xiàn)了約為9倍轉(zhuǎn)頻(相當(dāng)于一倍轉(zhuǎn)輪葉片頻率)的較高頻率，而且是尾水管內(nèi)壓力脈動(dòng)的主頻。考慮到抽水蓄能電站的地下廠房嵌固于圍巖中，且廠房混凝土結(jié)構(gòu)厚實(shí)，廠房結(jié)構(gòu)的自振頻率較高，中高頻的壓力脈動(dòng)對(duì)廠房振動(dòng)的影響程度反而更大。因此，有必要研究分析尾水管內(nèi)不同頻域脈動(dòng)壓力對(duì)廠房振動(dòng)的影響。

1 尾水管脈動(dòng)壓力時(shí)程信號(hào)

圖1為尾水管脈動(dòng)壓力傳感器布置示意圖，肘管段有a和b兩個(gè)測(cè)點(diǎn)，錐管段有c和d兩個(gè)測(cè)點(diǎn)，傳感器采樣頻率為4 800 Hz，時(shí)間持續(xù)10 s。選擇最接近額定出力的測(cè)量工況，該工況的詳細(xì)說(shuō)明見(jiàn)表1。采用FFT方法分析4個(gè)測(cè)點(diǎn)的脈動(dòng)壓力頻譜特性，對(duì)比情況見(jiàn)表2。

圖1 尾水管脈動(dòng)壓力傳感器布置示意圖

項(xiàng) 目數(shù)值導(dǎo)葉開(kāi)度γ/(°)20模型轉(zhuǎn)速/(r·min-1)1500模型水頭/m82．527對(duì)應(yīng)原型轉(zhuǎn)速/(r·min-1)500對(duì)應(yīng)原型水頭/m538．780對(duì)應(yīng)原型出力/MW298比例系數(shù)DP/DM7．65986

注：DP和DM分別對(duì)應(yīng)水輪機(jī)原型和模型轉(zhuǎn)輪直徑。

表2 FFT變換后的幅值和主頻

注：fn為原型或模型水輪機(jī)的轉(zhuǎn)頻，本文中原型轉(zhuǎn)頻為8.33 Hz，模型轉(zhuǎn)頻為25 Hz。

圖2(a)為尾水肘管下游側(cè)測(cè)點(diǎn)a前0.5 s的脈動(dòng)壓力時(shí)程曲線，圖2(b)為測(cè)點(diǎn)a的FFT振幅譜，從振幅譜中可以看出，壓力脈動(dòng)幅值較小，均在0.1%以下，其主頻f/fn=9。如若按照常規(guī)的處理辦法，設(shè)計(jì)低通濾波器將高頻部分濾掉，或者截?cái)囝l率取為5fn或6fn，那么仍會(huì)忽略尾水管脈動(dòng)壓力中9fn頻率的部分，與實(shí)際不符。本文分析了水輪機(jī)廠家提供的540 m水頭附近開(kāi)度γ從6°～26°的11種測(cè)量工況，發(fā)現(xiàn)當(dāng)導(dǎo)葉開(kāi)度γ逐漸接近20°時(shí)(對(duì)應(yīng)原型水輪機(jī)額定出力工況)，9fn頻率的脈動(dòng)壓力成分便越明顯，并逐漸成為尾水管脈動(dòng)壓力的主頻。由于本工程中轉(zhuǎn)輪葉片數(shù)為9，則轉(zhuǎn)輪葉片數(shù)頻率或轉(zhuǎn)輪葉片的葉柵頻率9fn將是對(duì)廠房振動(dòng)產(chǎn)生重要影響的頻率[10,11]，因此，有必要研究9fn附近頻率的尾水管壓力脈動(dòng)對(duì)廠房振動(dòng)的影響。

圖2 脈動(dòng)壓力時(shí)程曲線及FFT頻譜

在Matlab中設(shè)計(jì)FFT濾波器，將圖2(a)中的脈動(dòng)壓力信號(hào)進(jìn)行濾波處理，分離出7fn～9fn頻率范圍的信號(hào)，以及剩余部分的信號(hào)，分別見(jiàn)圖2(c)和圖2(e)，圖2(d)和圖2(f)為對(duì)應(yīng)的FFT振幅譜。圖2(c)中的脈動(dòng)壓力信號(hào)主頻為低頻，且幅值較小。圖2(c)和圖2(e)中的脈動(dòng)壓力信號(hào)疊加后幾乎可以得到圖2(a)中的原始脈動(dòng)壓力信號(hào)，這也是采用FFT法進(jìn)行濾波的好處，即在時(shí)域和頻域上均能達(dá)到較好的濾波效果[12]。

2 計(jì)算荷載和方案

圖3為洪屏抽水蓄能電站地下廠房結(jié)構(gòu)的三維有限元實(shí)體模型，其中圖3(a)為廠房整體三維模型，圖3(b)為沿機(jī)組中心線的橫剖面圖，有限元模型中廠房結(jié)構(gòu)按實(shí)際尺寸模擬，機(jī)組段長(zhǎng)度為22.5 m，廠房寬度22.0 m，廠房結(jié)構(gòu)柱尺寸一般為1.0 m×1.0 m，發(fā)電機(jī)層樓板厚度為0.6 m，中間層和水輪機(jī)層樓板厚度為0.5 m，并對(duì)球閥吊物孔附近的樓板局部加厚0.2 m，各層主梁尺寸一般為0.9 m×1.6 m或0.8 m×1.5 m，次梁尺寸一般為0.8 m×1.3 m或0.6 m×1.2 m。蝸殼外圍混凝土和機(jī)墩風(fēng)罩采用C30混凝土，其他結(jié)構(gòu)采用C25混凝土。該電站所處地理位置地質(zhì)條件較好，以Ⅱ類(lèi)圍巖為主。地下廠房的邊界條件采用法向彈簧模擬蝸殼層以上高程的上下游圍巖的約束作用，而蝸殼層以下部分混凝土與圍巖之間作為剛性連接，廠房底部直接施加固定約束。

圖3 地下廠房結(jié)構(gòu)三維模型

有關(guān)文獻(xiàn)討論了壓力脈動(dòng)原型與模型的轉(zhuǎn)換關(guān)系[8,13-15]，普遍認(rèn)為原型壓力脈動(dòng)幅值與水頭的比值與模型相同。對(duì)于壓力脈動(dòng)頻率的轉(zhuǎn)換關(guān)系，目前甚至沒(méi)有總結(jié)出半經(jīng)驗(yàn)公式[9]，因此，在沒(méi)有更準(zhǔn)確的轉(zhuǎn)換關(guān)系的情況下，本文使得原型和模型的f/fn保持一致。文獻(xiàn)[13]提到的尾水管內(nèi)壓力脈動(dòng)頻率換算關(guān)系見(jiàn)公式(1)。

(1)

式中:f為頻率；D為轉(zhuǎn)輪直徑；H為水頭；下標(biāo)P和M分別表示原型和模型。

代入表1中的數(shù)據(jù)可得fP/fM=1/3。由于模型轉(zhuǎn)頻是原型轉(zhuǎn)頻的3倍，因此本文保持原型和模型的f/fn一致，實(shí)際上與式(1)的換算關(guān)系效果相同。

圖1中4個(gè)測(cè)點(diǎn)的FFT振幅譜規(guī)律幾乎一致，但尾水錐管的兩個(gè)測(cè)點(diǎn)的脈動(dòng)壓力幅值相對(duì)較大，文獻(xiàn)[5]中也提到了類(lèi)似規(guī)律。而文獻(xiàn)[13]則認(rèn)為尾水肘管的壓力脈動(dòng)幅值較大。綜合考慮，本文采用尾水肘管下游側(cè)測(cè)點(diǎn)a的脈動(dòng)壓力實(shí)測(cè)值換算后對(duì)尾水管進(jìn)行加載，加載區(qū)域?yàn)閳D1中陰影部分。計(jì)算時(shí)廠房結(jié)構(gòu)阻尼比參考相關(guān)文獻(xiàn)取為0.02[13,15]。設(shè)計(jì)了表3所示的3種方案，除壓力脈動(dòng)荷載不同外，3種方案的其他計(jì)算條件均保持一致，對(duì)廠房結(jié)構(gòu)在尾水管脈動(dòng)壓力作用下的瞬態(tài)動(dòng)力響應(yīng)進(jìn)行分析。

表3 計(jì)算方案

3 廠房結(jié)構(gòu)的動(dòng)力響應(yīng)

3.1 樓板結(jié)構(gòu)的動(dòng)力響應(yīng)

考慮到發(fā)電機(jī)層樓板布置有各種電氣設(shè)備，加之工作人員在發(fā)電機(jī)層樓板上開(kāi)展工作等，發(fā)電機(jī)層樓板的振動(dòng)最容易被感知到。樓板長(zhǎng)度和寬度方向的尺寸與厚度方向相比要大得多，因此，樓板的豎直向動(dòng)力響應(yīng)一般較水平向要大。圖4給出了3種計(jì)算方案下發(fā)電機(jī)層樓板豎直向振動(dòng)速度均方根等值線；圖5為發(fā)電機(jī)層樓板豎直向振動(dòng)加速度均方根等值線。

圖4 發(fā)電機(jī)層樓板豎直向振動(dòng)速度均方根等值線(單位：mm/s)

圖5 發(fā)電機(jī)層樓板豎直向振動(dòng)加速度均方根等值線(單位：m/s2)

表4中給出了樓板第四象限特征點(diǎn)F1的振幅、振動(dòng)速度和振動(dòng)加速度的峰值和均方根值。

從圖4和圖5中可以看出，3種方案下振動(dòng)速度和振動(dòng)加速度的均方根分布一般是樓板邊緣處較大，越靠近風(fēng)罩的樓板振動(dòng)響應(yīng)越小，說(shuō)明尾水管脈動(dòng)壓力首先引起剛度較小的立柱和邊墻的振動(dòng)，然后再通過(guò)立柱和邊墻傳遞到樓板上。球閥吊物孔、樓梯孔和附屬吊物孔附近的樓板振動(dòng)響應(yīng)相對(duì)于風(fēng)罩附近有所增大，但增加幅度不大。

3種方案比較來(lái)看，方案三的豎直向振動(dòng)速度均方根、振動(dòng)加速度均方根要明顯大于方案二的相應(yīng)值，說(shuō)明頻率在7fn～9fn范圍內(nèi)的脈動(dòng)壓力對(duì)樓板豎直向振動(dòng)的貢獻(xiàn)程度更大。以特征點(diǎn)F1(見(jiàn)圖4)為例，方案三的豎直向振動(dòng)速度均方根為0.96 mm/s，方案二為0.45 mm/s；從能量角度看，發(fā)電機(jī)層樓板豎直向振動(dòng)能量的80%以上來(lái)源于7fn～9fn范圍內(nèi)脈動(dòng)壓力的貢獻(xiàn)。若尾水管脈動(dòng)壓力主頻不是以低頻為主，那么樓板豎直向振動(dòng)能量的絕大部分將來(lái)自于尾水管內(nèi)高頻脈動(dòng)壓力的貢獻(xiàn)。此外，越靠近尾水管的樓板區(qū)域(X>0,Z>0，即第四象限)振動(dòng)響應(yīng)程度越大，但沿樓板長(zhǎng)度和寬度方向的衰減很快，樓板振動(dòng)幅度較大的區(qū)域并不大。方案二情況下，樓板豎直向振動(dòng)速度和加速度的絕對(duì)數(shù)值不大，且整個(gè)樓板的速度和加速度分布比較均勻，說(shuō)明若尾水管內(nèi)壓力脈動(dòng)幅值較小且以低頻為主，其引起廠房樓板的振動(dòng)幅度將不大。

表4 樓板特征點(diǎn)F1豎直向振幅、振動(dòng)速度和振動(dòng)加速度

另外，討論3種計(jì)算方案下樓板特征點(diǎn)F1豎直向振動(dòng)速度的時(shí)程曲線及其頻譜特性，見(jiàn)圖6。7fn～9fn頻率范圍內(nèi)的尾水管脈動(dòng)壓力引起的樓板豎直向振動(dòng)速度要明顯大過(guò)以低頻為主的壓力脈動(dòng)引起的振動(dòng)速度，進(jìn)一步說(shuō)明尾水管中高頻壓力脈動(dòng)對(duì)廠房樓板振動(dòng)的影響程度較大。參照?qǐng)D2，可以發(fā)現(xiàn)樓板豎直向振動(dòng)速度的FFT頻譜與脈動(dòng)壓力頻譜整體看較為接近，但圖6中方案二振動(dòng)速度的主頻不明顯，說(shuō)明低頻脈動(dòng)壓力對(duì)樓板振動(dòng)速度的影響較小。

圖6 特征點(diǎn)F1豎直向振動(dòng)速度時(shí)程曲線及FFT頻譜

3.2 尾水錐管段的動(dòng)力響應(yīng)

考慮到尾水管進(jìn)人孔處的尾水錐管段鋼襯沒(méi)有混凝土包裹，加之脈動(dòng)壓力直接作用在尾水錐管上，該部分尾水管的動(dòng)力響應(yīng)值得關(guān)注。因此，選取尾水錐管段的一個(gè)特征點(diǎn)W1，如圖7所示。圖7中粗線框出的部分即為尾水管鋼襯無(wú)混凝土包裹范圍，尺寸達(dá)到了1.5 m×2.0 m(寬×高)，而W1點(diǎn)為尾水管鋼襯上最接近進(jìn)人孔中心的一個(gè)節(jié)點(diǎn)，即W1點(diǎn)處于尾水管在脈動(dòng)壓力作用下振動(dòng)響應(yīng)最大的區(qū)域。

W1點(diǎn)的徑向振幅、振動(dòng)速度和振動(dòng)加速度的峰值和均方根見(jiàn)表5。分析3種方案尾水管直錐段的振幅、振動(dòng)速度和振動(dòng)加速度，發(fā)現(xiàn)7fn～9fn頻率范圍的壓力脈動(dòng)對(duì)振幅的影響程度相對(duì)更大，對(duì)振動(dòng)速度的影響程度次之，對(duì)振動(dòng)加速度的影響程度很小。從特征點(diǎn)W1的徑向振幅峰值看，7fn～9fn頻率范圍的壓力脈動(dòng)引起的振幅峰值為2.65 μm，此時(shí)整體的振幅峰值為3.03 μm；振幅均方根值分別為1.00和1.27 μm。說(shuō)明尾水管振幅受壓力脈動(dòng)的影響更加敏感。圖8為W1點(diǎn)徑向振幅時(shí)程曲線及FFT頻譜。

表5 尾水錐管特征點(diǎn)W1 徑向振幅、振動(dòng)速度和振動(dòng)加速度

圖7 尾水錐管段特征點(diǎn)W1位置示意圖

圖8 特征點(diǎn)W1 徑向振幅時(shí)程曲線及FFT頻譜

4 結(jié) 語(yǔ)

通過(guò)以上分析，初步得到以下結(jié)論。

(1)尾水管脈動(dòng)壓力作用下樓板第四象限邊緣處振幅、振動(dòng)速度和振動(dòng)加速度較大，越靠近風(fēng)罩的樓板振動(dòng)響應(yīng)越小；孔洞附近的樓板振動(dòng)響應(yīng)有所增大，但增加 幅度不大。

(2)尾水管壓力脈動(dòng)中幅值較大的高頻成分貢獻(xiàn)了樓板豎直向振動(dòng)能量的絕大部分，但振動(dòng)速度和加速度沿樓板長(zhǎng)度和寬度方向的衰減速度較快。如果尾水管中出現(xiàn)了主頻為高頻的壓力脈動(dòng)，即使幅值較小，也需引起重視，關(guān)注其對(duì)水電站廠房樓板結(jié)構(gòu)的振動(dòng)影響。

(3)無(wú)混凝土包裹的尾水管進(jìn)人孔處鋼襯的振幅對(duì)壓力脈動(dòng)的幅值更加敏感，若尾水管內(nèi)出現(xiàn)幅值較大的壓力脈動(dòng)，將會(huì)引起較大的振幅；其振動(dòng)速度和加速度受整個(gè)頻域范圍內(nèi)的壓力脈動(dòng)影響程度大，單獨(dú)某一段頻域內(nèi)的脈動(dòng)壓力，引起尾水錐管段較大振動(dòng)速度和振動(dòng)加速度的可能性也不大。

□

[1] 王勝軍，胡清娟．寶泉抽水蓄能電站水泵水輪機(jī)模型驗(yàn)收試驗(yàn)[J]．水力發(fā)電，2008，34(10)：59-62．

[2] 章存建，高從闖，陳忠賓．溧陽(yáng)抽水蓄能電站水泵水輪機(jī)模型驗(yàn)收試驗(yàn)[J]．水力發(fā)電，2013，39(3)：79-81．

[3] 楊翼蜀，楊 波，龔 莉．黃金坪電站水輪機(jī)模型試驗(yàn)研究[J]．東方電機(jī)，2012，(2)：1-2，42．

[4] 劉海輝，李冬陽(yáng)，田樹(shù)鵬．蒲石河抽水蓄能電站水泵水輪機(jī)模型驗(yàn)收試驗(yàn)及性能分析[J]．東北水利水電，2011,(11)：46-48．

[5] 何成連，王正偉，邱 華．水輪機(jī)尾水管內(nèi)部壓力脈動(dòng)試驗(yàn)研究[J]．機(jī)械工程學(xué)報(bào)，2002，38(11)：62-65．

[6] 桂中華，唐 澍，潘羅平．混流式水輪機(jī)尾水管非定常流動(dòng)模擬及不規(guī)則壓力脈動(dòng)預(yù)測(cè)[J]．中國(guó)水利水電科學(xué)研究院學(xué)報(bào)，2006，4(1)：68-73．

[7] 王正偉，周凌九，何成連．尾水管壓力脈動(dòng)的模擬與現(xiàn)場(chǎng)實(shí)測(cè)[J]．清華大學(xué)學(xué)報(bào)(自然科學(xué)版)，2005，45(8)：1 138-1 141．

[8] 邵 杰，陳鐵軍，劉樹(shù)紅，等．軸流式原型和模型水輪機(jī)的壓力脈動(dòng)相似性分析[J]．水力發(fā)電學(xué)報(bào)，2009，28(4)：166-170．

[9] 鄭 源，汪寶羅，屈 波．混流式水輪機(jī)尾水管壓力脈動(dòng)研究綜述[J]．水力發(fā)電，2007，33(2)：66-69．

[10] 鄭東飛，易忠有．張河灣抽水蓄能電站水泵水輪機(jī)主要參數(shù)研究[C]∥抽水蓄能電站工程建設(shè)文集2008．北京：中國(guó)電力出版社，2011．

[11] 唐培甲．巖灘水電站水輪機(jī)振動(dòng)問(wèn)題的研究[J]．紅水河，2000，19(3)：59-62．

[12] 萬(wàn)永革．?dāng)?shù)字信號(hào)處理的MATLAB實(shí)現(xiàn)[M]．北京：科學(xué)出版社，2012．

[13] 歐陽(yáng)金惠，陳厚群，李德玉．三峽電站廠房結(jié)構(gòu)振動(dòng)計(jì)算與試驗(yàn)研究[J]．水利學(xué)報(bào)，2005，36(4)：484-490．

[14] 陳 婧，馬震岳，劉志明，等．水輪機(jī)壓力脈動(dòng)誘發(fā)廠房振動(dòng)分析[J]．水力發(fā)電，2004，30(5)：24-27．

[15] 蔣逵超，伍鶴皋，申 艷，等．脈動(dòng)壓力作用下直埋式蝸殼非線性動(dòng)力響應(yīng)[J]．水力發(fā)電學(xué)報(bào)，2007，26(4)：103-109．