邱志勤

(雅礱江流域水電開發(fā)有限公司，四川 成都 610051)

0 引 言

水輪發(fā)電機組是水電站的關(guān)鍵設(shè)備之一，其穩(wěn)定運行是水電站安全運行的基礎(chǔ)，也直接關(guān)系到電網(wǎng)的安全運行，決定著水電廠的經(jīng)濟效益和社會效益。二灘水電站單機容量大、調(diào)節(jié)性能好、機組運行效率高，是川渝電網(wǎng)的主力電站，承擔系統(tǒng)的調(diào)頻、調(diào)峰、調(diào)壓、斷面潮流控制及事故備用的任務(wù)。為滿足系統(tǒng)的調(diào)頻、調(diào)壓和調(diào)峰工作需要，電網(wǎng)需要二灘水電站機組頻繁變動出力以參與系統(tǒng)調(diào)節(jié)，這對二灘水電站機組的運行可靠性、調(diào)節(jié)速率與響應(yīng)時間等都提出了更高的要求，也必然會加劇機組的機械疲勞和磨損老化，另外水輪發(fā)電機組在運行中，還難免要受到泥沙磨損、氣蝕破壞等影響，客觀上需要電站生產(chǎn)運行人員能深入了解機組的運行振動特性，準確掌握機組的運行狀態(tài)，有效評估設(shè)備健康水平，從而科學的進行運行調(diào)控和檢修決策，避免對機組狀態(tài)掌握不足而導致的設(shè)備事故。

分析評估水輪發(fā)電機組運行狀態(tài)，可以從其運行振動和擺度數(shù)據(jù)著手，掌握和評估其運行振動水平，及時發(fā)現(xiàn)設(shè)備潛在缺陷和隱患，實時調(diào)整機組運行方式或及時安排機組檢修，防范于未然。

1 大數(shù)據(jù)分析方法的必要性

1.1 常規(guī)分析方式的運用

2003年，二灘電廠委托華中科技大學對二灘水電站機組振動課題做了專題研究，其研究成果揭示了二灘水電站機組運行的相關(guān)特性，對機組的運行狀態(tài)分析有一定的指導意義。其研究方法為采集機組不同負荷工況下的機組振動分析參量的穩(wěn)態(tài)運行數(shù)據(jù)，然后在二維平面里將不同的振動分析參量連點成線，繪制機組振動參量隨機組負荷的變化趨勢曲線，進而分析了機組的運行狀態(tài)[1-3]。這種研究方法由于數(shù)據(jù)具有離散性，只采集和分析了一部分的機組運行數(shù)據(jù)，且部分變量的采集數(shù)據(jù)值還有可能存在漂移偏差，僅能在一定程度上反映機組運行振動趨勢，并不能反映機組在不同水頭、不同負荷等多個變量下的完整運行特性。為全面分析機組在不同變量環(huán)境下的運行振動狀態(tài)，需要對機組不同時段、不同水頭、不同出力條件下的運行數(shù)據(jù)進行海量采集，盡最大可能從統(tǒng)計學上消除因明顯的機組故障和測量方式方法引起的數(shù)據(jù)偏差，避免導致的對機組振擺數(shù)據(jù)分析誤判的影響[4-5]。

1.2 機組穩(wěn)定性問題凸顯和工況變化

二灘水電站機組最大水頭189 m，最小水頭 135 m，變幅比1.37，具有典型水頭高、變幅大的特點。由于機組頻繁參與電網(wǎng)調(diào)峰調(diào)頻，有功出力大幅變化且調(diào)節(jié)頻繁，導致機組經(jīng)常性偏離最優(yōu)工況運行，其穩(wěn)定性問題相對突出。經(jīng)統(tǒng)計，僅2016年，二灘水電站機組出現(xiàn)振動頻繁報警的缺陷就有32例，電站自投產(chǎn)以來一共開展過12次機組A級檢修，歷年檢修原因統(tǒng)計見表1。經(jīng)分析，主要在于偏離機組最優(yōu)運行工況后，導致機組汽蝕磨損，設(shè)備老化加劇。另外，水輪發(fā)電機組進行了多次檢修調(diào)整和設(shè)備改造過后，機組原有的運行特性或可能已發(fā)生變化，尤其是機組運行診斷區(qū)間或已發(fā)生偏移，為進一步掌握機組的運行特性，有效指導電站機組生產(chǎn)運行工作，有必要對機組的歷史運行振動數(shù)據(jù)進行統(tǒng)計、分析和歸納總結(jié)。

表1 歷年二灘水電站機組A級檢修統(tǒng)計

2 R語言數(shù)據(jù)分析平臺的運用

2.1 R語言的介紹

二灘水電站機組振動運行數(shù)據(jù)龐大，且數(shù)據(jù)格式不完全一致，需要按一定的規(guī)律對振動數(shù)據(jù)進行篩選或運算，用傳統(tǒng)的數(shù)據(jù)處理方法不僅費時費力，也難以達到預期效果，但借助R語言數(shù)據(jù)分析平臺可快速繪制所需的圖形。

R語言是用于統(tǒng)計分析、繪圖的語言和操作環(huán)境。R語言作為統(tǒng)計計算和統(tǒng)計制圖的優(yōu)秀工具，能處理數(shù)據(jù)的收集、運算、探索、建模、可視化方面的工作；以能創(chuàng)建漂亮優(yōu)雅的圖形而聞名，已經(jīng)成了統(tǒng)計、預測分析和數(shù)據(jù)可視化的全球通用語言，也是現(xiàn)今最受歡迎的數(shù)據(jù)分析和可視化平臺之一[6]，其主要的特性有：

(1)R語言可以輕松地從各種類型的數(shù)據(jù)源導入數(shù)據(jù)，包括文本文件、數(shù)據(jù)庫管理系統(tǒng)、統(tǒng)計軟件，乃至專用的數(shù)據(jù)庫，并將這些數(shù)據(jù)進行篩選和批量運算。

(2)R語言是一個全面的統(tǒng)計研究平臺，提供了各式各樣的數(shù)據(jù)分析技術(shù)，幾乎任何類型的數(shù)據(jù)分析工作皆可在R語言中完成。

(3)R語言擁有頂尖水準的制圖功能，可以將復雜數(shù)據(jù)可視化。

2.2 數(shù)據(jù)處理和振動趨勢擬合

提取二灘水電站計算機監(jiān)控系統(tǒng)歷史數(shù)據(jù)庫中2015年7月～2017年8月的各機組運行振動數(shù)據(jù)，共計23.2萬余個CSV的文本格式文件，數(shù)據(jù)容量高達26 GB，約8 160萬個數(shù)據(jù)點。將每臺機組的水頭、功率、機組各部位的振動和擺度數(shù)據(jù)按時刻點做成一個R語言可以輸入的文件。將數(shù)據(jù)文件導入R語言平臺，并對所有數(shù)據(jù)進行篩選，剔除掉機組運行明顯異常的數(shù)據(jù)，有效避免對運行特性分析的影響。

二灘水電站共有6臺混流式水輪發(fā)電機組，每臺機組的負荷區(qū)間為0～550 MW，運行水頭區(qū)間為147～186 m，為更詳細的分析機組運行振動情況，將機組水頭按2 m左右分為一段，共分為15段，分別繪制機組水導、上導、下導、上機架、下機架、頂蓋等部位在各水頭段下的全負荷工況下的運行振動工況圖，共繪制990張圖，其中每一臺機組165張，然后對這些圖像進行逐一解析，可總結(jié)出其運行振動規(guī)律。以1、2號機組為例，繪制典型水頭水導擺度與機組有功功率之間的數(shù)據(jù)趨勢圖如圖1、2所示。

圖1 1號機水導+X方向擺度與機組有功功率關(guān)系(選取177 m水頭)

圖2 2號機水導+X方向擺度與有功功率關(guān)系(選取151 m水頭)

3 機組運行振動特性分析

借助R語言平臺，可以將機組運行歷史數(shù)據(jù)繪制成開展設(shè)備運行狀態(tài)分析所需要的數(shù)據(jù)趨勢圖，便于電站生產(chǎn)運行人員對機組運行特性進行分析，例如從圖1中，不難發(fā)現(xiàn)1號機組水導擺度在設(shè)計振動區(qū)間外呈現(xiàn)出的運行規(guī)律為，在低負荷工況下水導擺度隨有功增加而增加，在大負荷工況下水導擺度隨有功增加而減小，在機組設(shè)計針對區(qū)邊緣附近時水導擺度幅值最大，且機組實際運行振動區(qū)邊界較為清晰，具備以海量實際運行數(shù)據(jù)為基礎(chǔ)，復核并修正機組設(shè)計振動區(qū)間的可行性。

3.1 機組振動區(qū)邊緣運行特性分析

3.1.1 運行分析

研究發(fā)現(xiàn)，在部分水頭下，機組運行在振動區(qū)以外的邊緣區(qū)域時，如小負荷運行區(qū)靠上限區(qū)域、大負荷運行區(qū)靠下限區(qū)域，表現(xiàn)出機組振動、擺度偏大的情況，且6臺機組表現(xiàn)出一致的特點，如圖1、2所示。

機組長時間運行在振動區(qū)邊沿且振動值較大時，會引起磁拉力不均衡，造成轉(zhuǎn)子磁極匝間短路，引起機組連接部件松動、脫落，發(fā)生設(shè)備損壞事故，以及金屬部件的疲勞，從而出現(xiàn)裂紋，甚至形成裂縫直至斷裂等故障。

現(xiàn)有振動區(qū)劃分方式和AGC分配負荷邏輯會導致機組長時間運行在振動區(qū)邊緣。經(jīng)數(shù)據(jù)統(tǒng)計，二灘水電站機組運行在振動區(qū)邊沿的時間約占機組并網(wǎng)運行時間的14%。

3.1.2 解決辦法

GB/T 11348.5—2008《旋轉(zhuǎn)機械轉(zhuǎn)軸徑向振動的測量和評定 第5部分：水力發(fā)電廠和泵站機組》中將機組振擺值分為A、B、C、D共4個區(qū)間，所有振擺值都在A區(qū)范圍內(nèi)的機組稱為A區(qū)機組，以此類推。即通過A、B、C、D來評價機組的振動狀態(tài)。根據(jù)國標評價標準將機組運行狀態(tài)分為4個區(qū)域分別為：A區(qū)，機組運行狀態(tài)良好，屬于精品機組，不需要檢修調(diào)整，可以長期穩(wěn)定運行；B區(qū)，機組運行狀態(tài)合格，可以進行檢修調(diào)整，機組可長期運行；C區(qū)，機組運行狀態(tài)不好，報警，需要進行檢修調(diào)整，在該狀態(tài)下機組不可以長期運行；D區(qū)，危險，機組不可以長期運行，必須停機進行檢修，在該狀態(tài)下運行有很大危險。

水輪機組振動的測量和評定區(qū)域示意如圖3所示。

圖3 機組振動區(qū)域評價示意

根據(jù)圖3可知，二灘水電站機組振動值在150 μm以下為A區(qū)，屬于精品運行區(qū)域。為防止機組運行進入B區(qū)，可以結(jié)合大數(shù)據(jù)分析找出機組運行在B區(qū)的點，通過對各機組的振動區(qū)進行微調(diào)，使機組長期運行在A區(qū)范圍內(nèi)，以優(yōu)化機組運行工況。

3.2 機組振動異常評估策略

從振動幅值上來看三導軸承擺度幅值呈現(xiàn)如下規(guī)律：

(1)在小負荷區(qū)和大負荷區(qū)運行時，三導軸承的振擺關(guān)系，下導>水導>上導；

(2)在渦帶工況區(qū)水導的振動值最大，三導軸承的振擺關(guān)系，水導?下導>上導；

(3)振動幅值規(guī)律為三導軸承、頂蓋、機架的振動幅值隨水頭增加而增加。

結(jié)合對機組各部振動特征量運行數(shù)據(jù)的總結(jié)，歸納出機組在小負荷區(qū)、渦帶工況區(qū)、大負荷區(qū)的振動幅值范圍見表2，可以作為運行監(jiān)視過程中對機組異常或振動測量元器件是否損壞的判斷依據(jù)。

表2 穩(wěn)態(tài)運行工況下各部振動幅值 μm

3.3 機組振動運轉(zhuǎn)規(guī)律分析及運行建議

通過R語言平臺，可突破傳統(tǒng)的二維平面分析局限，將機組運行水頭、機組出力、振動擺度幅值等多個變量納入到三維空間內(nèi)，繪制出機組運行振動特性三維空間圖，便于電廠生產(chǎn)運行人員綜合分析各機組各水頭下全負荷運行特性，圖4和圖5分別為1號機組水導和頂蓋全水頭、全負荷工況下的振動特性三維圖示例。

圖4 1號機水導全水頭全負荷振動特性三維示意

圖5 1號機組頂蓋全水頭全負荷振動特性三維示意

通過對不同機組不同變量下三維運行特性圖的分析，可以深度挖掘水輪發(fā)電機組的運行規(guī)律和運行特性，也能發(fā)現(xiàn)一些設(shè)備潛在缺陷和隱患，更好的為電廠生產(chǎn)運行提供技術(shù)支持。通過對二灘水電站不同機組運行數(shù)據(jù)分析[7-10]，得出以下結(jié)論：

(1)一般情況下，二灘水電站機組在小負荷工況運行時，其振擺值大于大負荷運行工況，二灘水電站機組更適合運行在高負荷區(qū)。

(2)二灘水電站機組存在兩個振動峰值區(qū)。第一個在機組的最強渦帶振動區(qū)210～280 MW負荷區(qū)間；第二個位于機組在高水頭運行370～400 MW負荷區(qū)間。在這兩個負荷區(qū)間，機組均有點位的振動值或擺度值達到峰值，在手動調(diào)整機組負荷跨越振動區(qū)時，應(yīng)盡量避免機組在這兩個區(qū)間停留較長時間。另外，建議進一步研究監(jiān)控系統(tǒng)AGC負荷分配規(guī)則，考慮將此負荷調(diào)整邏輯加入監(jiān)控系統(tǒng)AGC程序，研究加快該負荷區(qū)域內(nèi)有功調(diào)節(jié)速率。

(3)機組在低水頭下，運行在100 MW負荷情況下，各點位均出現(xiàn)振動值增大的情況，呈現(xiàn)“小負荷振動區(qū)”特性。雖然在此運行工況下機組振擺值均在規(guī)定范圍內(nèi)，為使機組運行情況更加良好，分配機組負荷時，建議盡量規(guī)避此工況。

(4)1號機組的上導、下導和水導擺度值均偏大，6號機下導擺度值偏大，明顯大于其他機組，排除傳感器缺陷后，需要結(jié)合年度檢修工作對機組進行檢查處理。

4 結(jié) 論

水輪發(fā)電機組的振動機理比較復雜，直觀判斷和簡單的測試手段很難找到機組振動運行的規(guī)律和發(fā)現(xiàn)異常運行狀態(tài)，而運用R語言對機組運行振動數(shù)據(jù)進行分析和挖掘，無疑是一個有效探索機組運行振動規(guī)律和發(fā)現(xiàn)機組振動隱患的方法。