趙 茜，朱晟輝

(1.南通理工學(xué)院，江蘇 南通 226001； 2.華能國際電力江蘇能源開發(fā)有限公司南通電廠，江蘇 南通 226003)

0 引言

廠用電負(fù)荷模型的建立不僅可以精確地計(jì)算廠用電電壓，保障進(jìn)相能力仿真的精度[1]，同時(shí)根據(jù)廠用電負(fù)荷之間的關(guān)系，可以判斷系統(tǒng)內(nèi)是否存在故障[2]。通常，當(dāng)廠內(nèi)負(fù)荷所帶智能變送器發(fā)生故障或負(fù)荷尚未安裝智能變送器時(shí)，主控難以及時(shí)了解并判斷故障情況，需要通過專人值守，并結(jié)合現(xiàn)場(chǎng)異常情況、操作情況、元件動(dòng)作情況及結(jié)果進(jìn)行綜合分析，推斷故障原因[3]，不僅耗時(shí)耗力，還可能進(jìn)一步擴(kuò)大故障范圍。

本文基于最小二乘法和華能南通電廠SIS系統(tǒng)中負(fù)荷特性曲線及實(shí)際工況數(shù)據(jù)，建立準(zhǔn)確的廠用電負(fù)荷修正模型，通過其中某一負(fù)荷的有功功率計(jì)算該母線段負(fù)荷有功功率，根據(jù)信息可靠程度進(jìn)行系統(tǒng)故障診斷，并通過仿真算例來驗(yàn)證本文所提仿真方法的準(zhǔn)確性和有效性。

1 廠用電負(fù)荷靜態(tài)模型

廠用電負(fù)荷是電廠內(nèi)供電設(shè)計(jì)計(jì)算的基本依據(jù)，與廠用電母線額定有功、母線容量密切相關(guān)。要提高電廠內(nèi)供電設(shè)計(jì)的計(jì)算精度，需要建立準(zhǔn)確的廠用電負(fù)荷模型。

文獻(xiàn)[1]提出廠用電負(fù)荷靜態(tài)模型，如式(1)所示：

(1)

其中，P和Q分別為廠用電負(fù)荷的有功功率和無功功率；PN和QN分別為廠用電負(fù)荷的額定有功功率和無功功率；U為廠用電母線電壓標(biāo)幺值；Δf為頻率偏差的標(biāo)幺值。

表1為某一時(shí)間段實(shí)測(cè)華能南通電廠4A1母線用電負(fù)荷數(shù)據(jù)，代入式(1)計(jì)算出該時(shí)間段功率因數(shù)均值為0.933，該時(shí)間段4A1母線額定有功為11.921 MW，該時(shí)間段電廠4A1母線容量為12.776 8 MVA。因計(jì)算所得結(jié)果與電廠實(shí)測(cè)數(shù)據(jù)存在較大偏差，需要對(duì)式(1)進(jìn)行改進(jìn)。

1.1 廠用電負(fù)荷特性

電廠內(nèi)4A1母線上主要的廠用電中壓負(fù)荷如圖1所示，在廠用電負(fù)荷中，鍋爐和汽機(jī)系統(tǒng)的用電量占全廠的70%以上，這類系統(tǒng)中的負(fù)荷波動(dòng)嚴(yán)重影響廠用電負(fù)荷特性曲線。因此本文采用大負(fù)荷實(shí)際負(fù)載數(shù)據(jù)估算4A1母線有功功率。

表1 實(shí)測(cè)某時(shí)刻華能南通電廠4A1母線用電負(fù)荷數(shù)據(jù)

在4A1母線負(fù)載中，額定容量最大的負(fù)荷依次是吸風(fēng)機(jī)、一次風(fēng)機(jī)、循泵、凝泵、磨煤機(jī)、送風(fēng)機(jī)等。一次風(fēng)機(jī)、送風(fēng)機(jī)和吸風(fēng)機(jī)的負(fù)荷大小與機(jī)組所需的煤粉量、燃燒所需要的空氣、產(chǎn)生的煙氣量有著直接的關(guān)聯(lián)，此類負(fù)荷特性曲線的波動(dòng)會(huì)直接影響其余負(fù)荷的運(yùn)行方式，從而改變廠用電的有功功率。因此本文采用其中額定容量最大的負(fù)荷吸風(fēng)機(jī)對(duì)4A1母線的有功功率進(jìn)行預(yù)測(cè)。以實(shí)測(cè)南通華能電廠4A1母線負(fù)載數(shù)據(jù)為案例進(jìn)行分析，其中吸風(fēng)機(jī)負(fù)荷特性曲線與廠用電負(fù)荷特性曲線如圖2所示。

圖1 4A1母線負(fù)載示意圖

(a)吸風(fēng)機(jī)負(fù)荷特性曲線

(b)廠用電負(fù)荷特性曲線

由圖2可知，吸風(fēng)機(jī)負(fù)荷特性曲線與廠用電負(fù)荷特性曲線趨勢(shì)近似吻合。因此，本文將利用最小二乘法基于吸風(fēng)機(jī)負(fù)荷特性曲線預(yù)測(cè)生成廠用電負(fù)荷特性曲線，以此降低廠用電負(fù)荷數(shù)據(jù)誤差，并實(shí)時(shí)監(jiān)測(cè)廠用電負(fù)荷數(shù)據(jù)真實(shí)值的可信度。

1.2 最小二乘法

最小二乘法(又稱最小平方法)是一種數(shù)學(xué)優(yōu)化技術(shù)。它通過最小化誤差的平方和尋找數(shù)據(jù)的最佳函數(shù)匹配。利用最小二乘法可以簡便地求得未知的數(shù)據(jù)，并使這些求得的數(shù)據(jù)與實(shí)際數(shù)據(jù)之間誤差(真實(shí)目標(biāo)對(duì)象與擬合目標(biāo)對(duì)象的差)的平方和為最小[4]，使得擬合對(duì)象最大限度逼近目標(biāo)對(duì)象。

最小二乘法是解決曲線擬合問題最常用的方法。已知函數(shù)f(x)在若干點(diǎn)xi(i=1，2，…，n)處的值yi，便可根據(jù)插值原理來建立插值多項(xiàng)式作為f(x)的近似值。從給定的數(shù)據(jù)(xi，yi)出發(fā)，在某個(gè)函數(shù)類中尋求一個(gè)近似函數(shù)，來擬合這組數(shù)據(jù)。要求所得的近似曲線能最好地反映數(shù)據(jù)的基本趨勢(shì)?？傻煤瘮?shù)f(x)的表達(dá)式如式(2)所示。

f(x)=a1φ1(x)+a2φ2(x)+……+amφm(x)

(2)

其中，φk(x)是事先選定的一組線性無關(guān)的函數(shù)，ak是待定系數(shù)(k=1，2，…，m，m

1.3 基于最小二乘法的廠用電負(fù)荷預(yù)測(cè)

(3)

由于4A1母線有功功率還受其他負(fù)載影響，根據(jù)式(3)得到的結(jié)果與真實(shí)值誤差較大，因此引入噪聲以消除其他負(fù)載干擾，為簡便計(jì)算。如圖2所示，取一段時(shí)間內(nèi)吸風(fēng)機(jī)負(fù)載有功功率和4A1母線有功功率相比較，兩者比值大小如圖3所示。

圖3 吸風(fēng)機(jī)負(fù)載有功功率與4A1母線有功功率比值示意圖

由圖3可知，吸風(fēng)機(jī)負(fù)載有功功率與4A1母線有功功率比值的平均值為3.52，最大值為5.21，最小值為2.93，與額定計(jì)算的系數(shù)2.908誤差較大，因此采用額定吸風(fēng)機(jī)負(fù)載有功功率推算4A1母線負(fù)載有功功率較不合理，故采用其他算法計(jì)算4A1母線負(fù)載有功功率。

通過MATLAB 2014a平臺(tái)，采用最小二乘法對(duì)吸風(fēng)機(jī)負(fù)載有功功率與4A1母線有功功率數(shù)據(jù)進(jìn)行擬合計(jì)算，計(jì)算結(jié)果如圖4所示。

圖4 最小二乘法擬合曲線

由圖4可知，采用最小二乘法擬合成一次函數(shù)，函數(shù)系數(shù)值為[1.094 8， 4.144 3]，由此可知吸風(fēng)機(jī)負(fù)載有功功率與4A1母線有功功率的函數(shù)表達(dá)式為：

(4)

將本文最小二乘法擬合函數(shù)式(4)、式(3)中最小二乘法和額定估算得出的4A1母線負(fù)載有功功率、文獻(xiàn)[1]提出廠用電負(fù)荷靜態(tài)模型式(1)作對(duì)比，并將其與實(shí)際值相比較，采用隨機(jī)5組數(shù)據(jù)，分析對(duì)比兩者算法誤差大小，對(duì)比結(jié)果如表2所示。

由表2可知，額定估算誤差遠(yuǎn)大于最小二乘法，然而表中數(shù)據(jù)對(duì)比僅僅通過已知數(shù)據(jù)進(jìn)行比較，在實(shí)際情況中，最小二乘法精度尚不明確，由此通過采集5組隨機(jī)情況下的吸風(fēng)機(jī)負(fù)載有功功率數(shù)據(jù)，通過最小二乘法計(jì)算4A1母線有功功率數(shù)據(jù)，并以此與實(shí)際值進(jìn)行對(duì)比計(jì)算誤差，誤差結(jié)果如表3所示。

表2 額定估算與最小二乘法結(jié)果對(duì)比

表3 最小二乘法計(jì)算誤差

由表3可知，通過最小二乘法由吸風(fēng)機(jī)負(fù)載有功功率推算4A1母線負(fù)載有功功率的最大誤差為6.43%，誤差范圍小于10%，有較好的實(shí)際應(yīng)用價(jià)值。

2 結(jié)語

本文利用最小二乘法的回歸分析方法，對(duì)廠用電負(fù)荷建立修正模型。以華能南通電廠SIS系統(tǒng)中負(fù)荷特性曲線及實(shí)際工況數(shù)據(jù)為研究對(duì)象，根據(jù)負(fù)荷的有功功率計(jì)算該母線段負(fù)荷有功功率，判別系統(tǒng)內(nèi)是否存在故障，并且經(jīng)仿真算例來驗(yàn)證和分析本文所提方法。

本文僅通過最小二乘法對(duì)吸風(fēng)機(jī)負(fù)載有功功率和4A1母線負(fù)載有功功率之間的關(guān)系做一次方程，未考慮其他變量對(duì)系統(tǒng)的影響，未來擬采用卡爾曼濾波和BP神經(jīng)網(wǎng)絡(luò)加入更多變量，以消除噪聲對(duì)結(jié)果的影響。