李 毅，龍 瀾，吳小虎，劉一力，俞柯丞

（中國移動寧波分公司，浙江 寧波 315048）

0 引 言

閥控式鉛酸蓄 Valve-Regulated Lead Acid(VRLA)電池作為數(shù)據(jù)中心不間斷電源（Uninterruptible Power Supply，UPS）系統(tǒng)的重要組成部分，在市電突然斷電時，鉛酸蓄電池組需保障柴油發(fā)電機啟動期間的用電需求。對鉛酸蓄電池單體剩余容量的準(zhǔn)確預(yù)測可以估算后備時長，提高電池的循環(huán)壽命，是保障數(shù)據(jù)中心業(yè)務(wù)連續(xù)性和安全性的重要基礎(chǔ)[1,2]。

蓄電池容量是表征蓄電池單體儲存電能的物理量，即在特定工作條件下，蓄電池完全放電所放出的電能，是衡量其對外放電能力的重要標(biāo)志。隨著蓄電池充放電循環(huán)次數(shù)的增加以及外部環(huán)境的變化，各個電池單體的剩余容量會存在較大差異，串聯(lián)的電池組中某一蓄電池單體性能劣化，從而使得整組電池的性能亟劇下降。為保障鉛酸蓄電池組處于最優(yōu)運行狀況，至少每2年需進(jìn)行1次核對性放電試驗，以檢驗電池單體的剩余容量和性能優(yōu)劣[3]。

常用的測量蓄電池剩余容量的方法主要分為2種，第一種為核對性充放電法。在規(guī)定環(huán)境溫度下，以某一特定的放電率將蓄電池組放電至終止電壓以測量電池容量，進(jìn)而發(fā)現(xiàn)落后電池。但是該方法需要將蓄電池從直流系統(tǒng)中退出，加裝臨時負(fù)載進(jìn)行測試，重復(fù)充放電測試會加速電池老化，縮短電池的循環(huán)使用壽命，也會存在欠充、過充及短路等安全隱患[4]。第二種是安時積分法。對連續(xù)采集到的電流參數(shù)進(jìn)行時間積分得到電池組的剩余容量，但采用該方法進(jìn)行測量時將額外增加干擾誤差，導(dǎo)致測量結(jié)果不準(zhǔn)確。

考慮到蓄電池的剩余容量與電池單體的端電壓和工作環(huán)境（溫度、工作時間）等因素有關(guān)，當(dāng)電池端電壓小于截止電壓時，則認(rèn)為該電池單體報廢，因此對蓄電池剩余容量的預(yù)測可以轉(zhuǎn)化為對電池單體端電壓的預(yù)測。

根據(jù)以上分析，本研究將基于大量現(xiàn)場實測數(shù)據(jù)，搭建一種結(jié)構(gòu)簡單、運算穩(wěn)定的多元線性模型用于蓄電池剩余容量的預(yù)測，并通過對采集數(shù)據(jù)的預(yù)測分析，以驗證本文所提出模型的準(zhǔn)確性。

1 蓄電池技術(shù)指標(biāo)

1.1 蓄電池基本參數(shù)

蓄電池容量用于衡量儲存電能的多少，其定義方法是在特定的放電強度、放電電流及放電終止電壓下，完全充電的電池單體所放出的電量或放電時間長短即為蓄電池的容量值，通常用Q表示[5]。其容量計算式為

式中：K為安全系數(shù)，一般取1.25；I為負(fù)荷電流，A；T為放電小時數(shù)，h；η為蓄電池放電容量系數(shù)；t為蓄電池所處的實際環(huán)境溫度，℃；α為電池溫度系數(shù)，通常取0.008。

1.2 蓄電池數(shù)據(jù)采集及預(yù)處理

1.2.1 數(shù)據(jù)采集

本研究所采用的測試數(shù)據(jù)取自浙江省寧波市某數(shù)據(jù)中心，配置有獨立的UPS系統(tǒng)，UPS系統(tǒng)采用2N架構(gòu)。閥控式密封鉛酸蓄電池系統(tǒng)配置及性能參數(shù)如表1所示。

表1 鉛酸蓄電池系統(tǒng)配置及性能參數(shù)表

鉛酸蓄電池通過核容放電測試方法獲得現(xiàn)場實測參數(shù)，放電測試的基本步驟如下文所述。

（1）分離目標(biāo)電池組：分?jǐn)嗄繕?biāo)電池組開關(guān)箱內(nèi)的斷路器，用絕緣扳手拆除目標(biāo)電池組正負(fù)極總電纜（及120節(jié)處的中心線）；1+1 UPS系統(tǒng)的一組電池拆除后，需要將電池開關(guān)箱內(nèi)的斷路器重新合上；目標(biāo)電池組從系統(tǒng)中脫離。

（2）連接測試儀表：運用導(dǎo)線連接好放電測試儀表的電源線、正負(fù)極總線。

（3）設(shè)置測試儀表參數(shù)：標(biāo)稱容量100 Ah（150 Ah、200 Ah）、小時率0.1C3、預(yù)放電流25 A（37.5 A、50 A）、預(yù)放容量 75 Ah（112.5 Ah、150 Ah）、預(yù)放時間3 h、單體下限1.8 V、組端下限432 V、每組單體240個、下限個數(shù)6個。

（4）測試并存儲數(shù)據(jù)：確認(rèn)無誤后，閉合放電測試儀表顯示屏下方“斷路器”，開始放電（萬用表核實電壓電流），并用U盤轉(zhuǎn)存測試數(shù)據(jù)。

經(jīng)測試后得到電池單體的端電壓曲線如圖1所示。

圖1 電池測試端電壓曲線（以10節(jié)電池單體為例）

1.2.2 數(shù)據(jù)預(yù)處理

在數(shù)據(jù)采集過程中，原始數(shù)據(jù)會因傳感器故障等問題存在異常值和缺失值?；陔妷簲?shù)據(jù)具有連續(xù)性的假設(shè)，本研究運用滑動窗口方式發(fā)現(xiàn)異常值。與窗口中其他點的距離超過給定閾值的數(shù)據(jù)點被視為異常值，需要將該樣本刪除，然后通過線性插值的方式補全缺失值。

蓄電池容量預(yù)測所用的數(shù)據(jù)特征具有不同的量綱單位，如果不消除特征之間的量綱影響，將其直接輸入預(yù)測模型，將導(dǎo)致模型僅受部分量綱較大的特征變量的影響，進(jìn)而影響模型的預(yù)測性能。因此，需要對原始數(shù)據(jù)集進(jìn)行標(biāo)準(zhǔn)化處理以消除由于量綱不同產(chǎn)生的特征重要性差異。目前常用的數(shù)據(jù)標(biāo)準(zhǔn)化處理方式包括：離差標(biāo)準(zhǔn)化（min-max標(biāo)準(zhǔn)化）和標(biāo)準(zhǔn)差標(biāo)準(zhǔn)化（Z-score標(biāo)準(zhǔn)化），表2對比了2種處理方法的缺點和適用場景[6]。

表2 數(shù)據(jù)標(biāo)準(zhǔn)化處理方法比較

在數(shù)據(jù)采集過程中，數(shù)據(jù)樣本的最大值和最小值不會為某一固定值，它們將隨著采集時長以及采集頻率的增加而不斷變化，難以確定原始數(shù)據(jù)集的xmax和xmin，因此本研究采用基于標(biāo)準(zhǔn)差的方法對特征進(jìn)行標(biāo)準(zhǔn)化處理。經(jīng)過標(biāo)準(zhǔn)化處理后的特征樣本符合均值為0，標(biāo)準(zhǔn)差為1的標(biāo)準(zhǔn)正態(tài)分布。

1.2.3 特征數(shù)據(jù)集構(gòu)建

特征工程構(gòu)造得到的數(shù)據(jù)相關(guān)性通過Pearson相關(guān)系數(shù)r實現(xiàn)量化評價，具體表達(dá)式為

式中：N為數(shù)據(jù)集的長度；xi和yi代表2組特征向量。

根據(jù)專家知識，蓄電池剩余容量主要和放電電流、溫度、電池內(nèi)阻等眾多因素相關(guān)，考慮到數(shù)據(jù)測試平臺的限制，電池溫度和內(nèi)阻難以準(zhǔn)確測量。當(dāng)電池端電壓小于截止電壓時，則認(rèn)為該電池單體報廢，因此對蓄電池剩余容量的預(yù)測可以轉(zhuǎn)化為對電池單體端電壓的預(yù)測。由現(xiàn)場測試的端電壓波動曲線可以發(fā)現(xiàn)，在初始放電階段，電池單體以較大速率放電迅速達(dá)到谷電壓，然后以較小的幅度回升至電壓極大值，最后以較小的速率放電至截止電壓。因此，電池單體的端電壓與放電時間t、初始放電電壓Ustart、峰電壓Upeak、總電壓U、總電流I等因素有關(guān)。圖2為根據(jù)Pearson相關(guān)系數(shù)繪制的特征熱力圖。

圖2 特征向量相關(guān)系數(shù)矩陣熱圖

從圖中可以看出，電池單體端電壓與峰電壓Upeak、放電時間t、總電壓U等特征變量存在較強的相關(guān)性，因此本研究選擇這些變量作為模型的輸入特征。

2 蓄電池剩余容量預(yù)測模型

2.1 模型介紹

根據(jù)現(xiàn)場測試的電池單體端電壓波動曲線，可以發(fā)現(xiàn)端電壓近似為線性變化，因此本研究擬采用多元線性回歸模型對電池單體的端電壓進(jìn)行預(yù)測，進(jìn)而預(yù)測蓄電池的剩余容量。

多元線性回歸模型是基于對2個或2個以上的特征變量與標(biāo)簽值的相關(guān)分析，建立線性回歸模型進(jìn)行預(yù)測的方法[7]。多元線性回歸模型的一般形式為

式中：表示觀測值y的估計值；表示p個回歸系數(shù)；x1，x2…xp代表p個解釋變量；ε為誤差項，用于解釋由隨機變量引起觀測值y的變化。

2.2 評價指標(biāo)

為了評估模型的預(yù)測性能，本研究選用均方根誤差RMSE和擬合優(yōu)度R2來檢驗和衡量算法模型的效果和精確度。其中RMSE用于表示實際值與預(yù)測結(jié)果之間的偏差，R2用于衡量預(yù)測模型的精度，該值越接近1，意味著模型的預(yù)測性能越好[8]。RMSE和R2的表達(dá)式為

式中：m為預(yù)測的總次數(shù)；yt和t分別代表t時刻端電壓的真實值和預(yù)測值；t為電池單體實測端電壓的均值。

3 實驗驗證

3.1 模型訓(xùn)練

為了充分利用數(shù)據(jù)集，本研究采用80%的樣本作為訓(xùn)練集，20%的樣本作為測試集，以保證模型預(yù)測的準(zhǔn)確性。簡單多元線性回歸模型基于Python的機器學(xué)習(xí)庫scikit_learn庫實現(xiàn)。

本研究選用第21組電池的實測數(shù)據(jù)集為例，用經(jīng)過訓(xùn)練的模型進(jìn)行測試，得到該電池單體端電壓預(yù)測值與實測值的對比結(jié)果如圖3所示。

圖3 第21節(jié)電池單體端電壓預(yù)測值與實測值對比結(jié)果

3.2 模型評估

為了評估多元線性回歸模型的預(yù)測性能，本研究進(jìn)一步對比了該模型的均方根誤差RMSE和擬合優(yōu)度R2值，結(jié)果如表3所示。

表3 基于多元線性回歸的電池單體端電壓預(yù)測性能

RMSE越小，R2越大，表示模型預(yù)測精度越高。結(jié)果表明，該模型的預(yù)測均方根誤差（RMSE）為0.0526 V，擬合優(yōu)度R2為78.4%，基本滿足工程預(yù)測需求。

4 結(jié) 論

本研究針對數(shù)據(jù)機房鉛酸蓄電池剩余容量難以預(yù)測性的問題，將剩余容量預(yù)測轉(zhuǎn)化為對電池單體兩端端電壓的預(yù)測，進(jìn)一步提出了基于簡單多元線性回歸的蓄電池剩余容量預(yù)測方法，并運用浙江省寧波市某數(shù)據(jù)中心的大量實測數(shù)據(jù)作為算例進(jìn)行了模型仿真訓(xùn)練和測試。結(jié)果表明基于多元線性回歸方法建立的蓄電池容量預(yù)測模型具有較優(yōu)的預(yù)測性能，基本滿足工程需求。