李存文，王在華，陳濤，馮前偉，徐克濤，張楊

（1．華電電力科學研究院有限公司，浙江省 杭州市310030；2．國網(wǎng)浙江省電力有限公司電力科學研究院，浙江省 杭州市310014；3．華電新疆發(fā)電有限公司，新疆維吾爾自治區(qū) 烏魯木齊市830063）

0 引言

為了進一步響應《煤電節(jié)能減排升級與改造行動計劃(2014—2020年)》[1]的號召，燃煤電站仍在尋求管理和技術手段，以實現(xiàn)排放達標的同時降低煙氣環(huán)保系統(tǒng)的運行維護成本，即環(huán)保系統(tǒng)的精細化運行控制顯得越發(fā)重要。而對煙氣環(huán)保系統(tǒng)性能的準確認識是其精細化控制的基礎。

獲取環(huán)保系統(tǒng)性能的途徑有性能試驗、在線監(jiān)測等。性能試驗依據(jù)于標準規(guī)范，指標全面，但試驗條件嚴苛，化驗耗時長；在線監(jiān)測主要由現(xiàn)場表計完成，指標的實時性、延續(xù)性好，但存在一些參數(shù)(如漿液品質(zhì)等)無法監(jiān)測的情況。如何準確、快速、高效地獲取環(huán)保系統(tǒng)實時性能，仍是科研技術人員需要研究的重要課題[2]。

數(shù)據(jù)挖掘技術融合人工智能、數(shù)據(jù)庫、高性能計算等多學科的成果，是將大量數(shù)據(jù)轉(zhuǎn)換為信息和知識的強大工具[3-5]。運用數(shù)據(jù)挖掘技術，研究獲取環(huán)保系統(tǒng)性能的新方法，為環(huán)保系統(tǒng)精細化控制作技術支撐，實現(xiàn)生態(tài)環(huán)境治理能力現(xiàn)代化[6]。

1 環(huán)保系統(tǒng)數(shù)據(jù)挖掘算法

1.1 數(shù)據(jù)預處理

電站運行參數(shù)來自運行實際，其數(shù)據(jù)一致性和完整性往往無法達到數(shù)據(jù)挖掘的要求。因而，數(shù)據(jù)預處理必不可少，有必要進一步研究并實現(xiàn)運行數(shù)據(jù)預處理方法的標準化。此外，鑒于運行數(shù)據(jù)在采集時已完成數(shù)量化，即自帶數(shù)值屬性，因此，數(shù)據(jù)預處理包括數(shù)據(jù)集成、數(shù)據(jù)清理和變換、統(tǒng)計分析及可視化，具體步驟如圖1所示。

數(shù)據(jù)集成將分散數(shù)據(jù)集成為數(shù)據(jù)集。

運行數(shù)據(jù)不可避免地包含異常、缺省數(shù)據(jù)和噪聲，因而，為針對性地提高數(shù)據(jù)質(zhì)量，需進行以下數(shù)據(jù)清理：一致性檢查并清理缺省值、空值和異常值；移動平均濾波法光滑去噪；剔除違背工藝邏輯的工況數(shù)據(jù)。

為規(guī)避量綱對數(shù)據(jù)認識的影響，同時為數(shù)據(jù)模型建模及迭代優(yōu)化考量，數(shù)據(jù)變換的主要內(nèi)容有：max-min 歸一化，穩(wěn)定梯度下降算法的求優(yōu)速度，以加快模型訓練速度；對有類似開或關2種狀態(tài)的參數(shù)進行二值化處理，壓縮模型變量數(shù)值域及量綱；引入新指標，如電耗量、還原劑耗量等，豐富對象性能指標。

統(tǒng)計分析及可視化包括兩部分：識別研究對象的物理邊界、數(shù)據(jù)特征及規(guī)模；將信息從數(shù)據(jù)中提煉為可解釋的知識，包括圖、表等。

1.2 環(huán)保系統(tǒng)關鍵因素分析

按所采集的測點參數(shù)計，數(shù)據(jù)集的維度通常非常大，這無疑加大了數(shù)據(jù)建模的工作負擔。因此，須對數(shù)據(jù)集進行降維，同時提取環(huán)保系統(tǒng)性能指標的關鍵因素集，而主成分分析方法在這方面十分有效[7-11]。

算法步驟如下：

1）以采樣時刻劃分數(shù)據(jù)元組，構成數(shù)據(jù)集O；

2）元組集O的標準化成矩陣S；

3）求矩陣S的相關系數(shù)矩陣R；

4）求取矩陣R的特征根λi及其特征向量ηi；

5）將標準化后的指標變量轉(zhuǎn)換為主成分ξj；

6）分析各主成分權重系數(shù)，并締結結果集U；

7）設定遴選權重系數(shù)規(guī)則，將主成分劃入關鍵因素集。

1.3 環(huán)保系統(tǒng)性能數(shù)據(jù)模型

煙氣凈化處理包含復雜的物理、化學反應，且傳熱與傳質(zhì)耦合，參數(shù)間關系復雜，使得機理模型雖然能保證邏輯趨勢的正確性，但分析的精細度面臨瓶頸，而數(shù)據(jù)模型因強擬合而具備高精度[12-13]。另外，環(huán)保系統(tǒng)內(nèi)部傳質(zhì)存在磨損、漿液中有害物質(zhì)聚集等因素，致使系統(tǒng)性能隨運行時間而衰減，因而，定期地對數(shù)據(jù)模型進行訓練更新，適應環(huán)保系統(tǒng)特性變化，既可保證模型的正確性，還可規(guī)避復雜的機理建模過程。

回歸分析是構建復雜對象特性模型的高效方法[14-18]。以關鍵因素集作為數(shù)據(jù)模型的輸入，融合多項式、Lasso、Ridge和神經(jīng)網(wǎng)絡算法[19-21]，構建環(huán)保系統(tǒng)數(shù)據(jù)模型，如圖2 所示。此外，為了均衡優(yōu)化速度，模型迭代優(yōu)化選用梯度下降算法進行。

圖2 煙氣環(huán)保系統(tǒng)數(shù)據(jù)模型搭建方法Fig.2 Method of building data model of flue gas environmental protection system

數(shù)據(jù)模型的評價指標較為豐富，結合實際并規(guī)避量綱及數(shù)值量級對模型精度的影響，綜合研究后選用均方根誤差γRMSE、決定系數(shù)γR2、平均絕對百分比誤差γMAPE這3 個指標來評價數(shù)據(jù)模型的準確性、有效性。

2 數(shù)據(jù)模型應用案例分析

某燃煤電站660 MW機組煙氣脫硫系統(tǒng)采用石灰石-石膏濕法工藝，“1爐1塔”配置，噴淋空塔并設置5層噴淋層、1層托盤。借助于廠級實時監(jiān)控信息系統(tǒng)(supervisory information system, SIS)，采集到2018年4—6月為期50天、包括107個參數(shù)的運行數(shù)據(jù)。依托Matlab平臺，將數(shù)據(jù)塊集成為數(shù)據(jù)矩陣，以采樣時刻為元組劃分原則，并去除重復數(shù)據(jù)，得到共計13983條數(shù)據(jù)元組。根據(jù)運行參數(shù)在時域上連續(xù)的特點，采用前時刻值填充后時刻缺省值，并運用移動平均濾波的方法進行去噪。

數(shù)據(jù)平滑濾波后，機組功率的瞬時波動大大減弱、數(shù)據(jù)趨勢及均值不失真，也說明數(shù)據(jù)去噪的必要性，如圖3所示。

圖3 數(shù)據(jù)平滑濾波前后結果Fig.3 Results before and after data smoothing filtering

剔除與物質(zhì)守恒及煙氣漏風等客觀規(guī)律相違背的工況數(shù)據(jù)，主要有：1）煙氣出口氧含量小于入口氧含量的數(shù)據(jù)；2）煙道兩側煙氣氧含量絕對偏差超過3.0%的數(shù)據(jù)；3）折算后出口污染物濃度大于入口的數(shù)據(jù)；4）機組有功功率小于等于零的數(shù)據(jù)。數(shù)據(jù)清理后，得到數(shù)據(jù)質(zhì)量高的8616個元組數(shù)據(jù)，即8616×111數(shù)據(jù)矩陣，占原數(shù)據(jù)集約61.62%，表明所采集運行數(shù)據(jù)的有效性高。

數(shù)據(jù)變換主要在建模中完成，而電耗和物耗指標為數(shù)據(jù)模型提供先驗知識與信息。

最后，統(tǒng)計環(huán)保系統(tǒng)參數(shù)的最值和眾數(shù)，得到該環(huán)保系統(tǒng)的物理邊界，見表1。

表1 運行數(shù)據(jù)的統(tǒng)計分析Tab.1 Statistical analysis on operation data

脫硫效率是表征煙氣脫硫系統(tǒng)性能的關鍵指標。借助于阿里云平臺，分析各參數(shù)對脫硫效率的影響，構成關鍵因素集。結果表明，影響脫硫效率的關鍵影響因素有原煙氣SO2濃度、循環(huán)泵電流、漿液pH 值、原煙氣煙塵濃度、漿液密度、煙氣量、原煙氣溫度等，如圖4 所示。該結果與機理定性分析結論一致，表明主成分分析方法具有較好的適用性和有效性，為后續(xù)脫硫系統(tǒng)的運行優(yōu)化及調(diào)整提供定量的參考。

圖4 關鍵因素集Fig.4 Key factor set

梳理關鍵因素集的結果，作為煙氣脫硫系統(tǒng)性能數(shù)據(jù)模型的輸入，采用多算法融合的回歸分析技術構建環(huán)保系統(tǒng)性能數(shù)據(jù)模型，最終通過預測凈煙氣SO2濃度、凈煙氣溫度、脫硫電耗、石灰石耗量，以達到預測煙氣脫硫系統(tǒng)性能的目的，同時完成算法調(diào)優(yōu)和模型訓練，如圖5所示。

圖5 煙氣脫硫系統(tǒng)性能預測數(shù)據(jù)模型Fig.5 Data model for performance prediction of flue gas desulphurization system

數(shù)據(jù)模型的建立成功表明了方法的可行性，而模型的準確性通過3個指標來評價界定，見表2。

表2 數(shù)據(jù)模型的主要評價指標與效果Tab.2 Main evaluation index and effect of the data model

γRMSE分別為1.09 mg/m3、0.41℃、88.82 kW，與其眾數(shù)和最大值相比，數(shù)值和量級均在工程允許范圍內(nèi)；γR2均大于0.85，模型的效果較好，能高概率地保證模型趨勢的正確性；γMAPE均小于0.10，模型能高準確度地模擬對象的輸入和輸出，間接表明算法和參數(shù)匹配較好。

為進一步驗證模型的有效性和準確性，選取不同機組負荷下運行數(shù)據(jù)，對模型進行測試，結果見表3和表4。

表3 煙氣脫硫系統(tǒng)性能預測數(shù)據(jù)模型測試案例Tab.3 Test case of performance prediction data model for flue gas desulfurization system

表4 數(shù)據(jù)模型測試結果Tab.4 Testing results of data model

石灰石制漿系統(tǒng)具有時滯延后特性，即石灰石耗量的精確確定不完全依托運行數(shù)據(jù)，因此石灰石耗量的準確性對于模型效果評價的可靠性有待深入研究。剔除石灰石耗量后，各參數(shù)的相對偏差均小于6.50%，量級與模型的γMAPE相當，且數(shù)據(jù)模型能準確復現(xiàn)負荷變化下環(huán)保系統(tǒng)狀態(tài)變化規(guī)律和趨勢，具有較強的負荷適應性，驗證了性能預測數(shù)據(jù)模型的有效性和準確性。

3 結論

突破傳統(tǒng)機理和試驗研究的方式，以燃煤電站煙氣脫硫系統(tǒng)運行數(shù)據(jù)作為切入點，運用大數(shù)據(jù)技術挖掘出影響脫硫系統(tǒng)性能的關鍵因素集，合理梳理關鍵因素并將其作為數(shù)據(jù)模型的輸入，融合回歸分析技術的多種算法，構建了煙氣脫硫系統(tǒng)性能預測數(shù)據(jù)模型。實例證明，訓練后的模型能高精度地復現(xiàn)不同負荷工況的脫硫系統(tǒng)性能，還具有較好的靈活性和可拓展性，為后續(xù)燃煤電站環(huán)保裝備的狀態(tài)監(jiān)測、趨勢分析、運行優(yōu)化提供堅實的基礎。相比于機理研究和試驗方法，該研究方法更為高效靈活、簡潔明了，豐富了燃煤電站環(huán)保研究思路，且具有啟示意義。