李 鋒,李 峰,楊榮杰,王洪祥
(貴州北盤(pán)江電力股份有限公司董箐分公司,貴州 黔西南布依族苗族自治州 562200)
在水輪機(jī)控制系統(tǒng)中需診斷采集信號(hào),信號(hào)正常時(shí)方可使用?,F(xiàn)有的信號(hào)診斷方法(斷線、超量程、階躍或突變)效果單一,由此提出一種從測(cè)量值倒推函數(shù)的信號(hào)診斷算法,從倒推中求取信號(hào)的一階系數(shù)、二階系數(shù)及高階系數(shù),以此作為信號(hào)某階的特征值。對(duì)水輪機(jī)控制系統(tǒng)中的重要信號(hào)進(jìn)行特性分析,診斷此信號(hào)是否正常。此方法適用范圍廣,可互補(bǔ)采用。
根據(jù)信號(hào)量屬性不同,可對(duì)診斷方式進(jìn)行分類。本研究提出的算法原理是從測(cè)量值來(lái)倒推信號(hào)變化特征的函數(shù)方程,對(duì)其做微分處理,求取其一階項(xiàng)、二階項(xiàng)及高階項(xiàng)系數(shù),從而得到特征值,主要應(yīng)用于低頻信號(hào)。
假設(shè)一個(gè)信號(hào)量的一般性函數(shù)為:
y=a·t+b·t2+c·t3
(1)
y為信號(hào);t為時(shí)間;a、b、c為需要得到的一階系數(shù)、二階系數(shù)及三階系數(shù)。
分別做多階微分處理,得到一、二、三階項(xiàng)的系數(shù)。如果直接反推是有難度的。如果做多階差分,那么根據(jù)多階差分原理:當(dāng)時(shí)間t的增量足夠小時(shí),一階差分與數(shù)值增量Δt為線性關(guān)系;二階差分與數(shù)值增量Δt的平方為線性關(guān)系;三階差分與數(shù)值增量Δt的三次方為線性關(guān)系;N階差分與數(shù)值增量Δt的n次方為線性關(guān)系。求取信號(hào)的多階差分,即增量Δt的n階線性關(guān)系,或者差商,分別取一階差商-a、二階差商-b、三階差商-c、N階差商-n。差商的定義為:n階差分與Δt增量的n次方的比值,即得到信號(hào)y的N階特征值:
(2)
取Δt->0趨近于極限時(shí),某階的差商就是某階導(dǎo)數(shù),變成了上述一般表達(dá)式的求導(dǎo),即微分關(guān)系?,F(xiàn)實(shí)應(yīng)用中通常不會(huì)選擇Δt->0。以Δt=1為例,取得的特征系數(shù)不是嚴(yán)格相等的,而是相似的,但依然能夠體現(xiàn)出某一階項(xiàng)系數(shù)的特性。
在水輪機(jī)控制系統(tǒng)工程量中,除機(jī)組頻率信號(hào)、主配開(kāi)度信號(hào)及水頭信號(hào)(或者進(jìn)水口壓力)可能出現(xiàn)二階項(xiàng)系數(shù)以外,其余信號(hào)正常時(shí)大部分時(shí)間都只有一階項(xiàng)。信號(hào)的一階項(xiàng)系數(shù)a的最大值基本穩(wěn)定或在設(shè)計(jì)中能找出與之對(duì)應(yīng)的理論最大值,對(duì)可能出現(xiàn)二階項(xiàng)特征的工程量,其二階項(xiàng)的系數(shù)b變化幅值不大,而高階項(xiàng)的系數(shù)都為零或幾乎可以忽略為零。由于求取無(wú)限多階差分的現(xiàn)實(shí)意義不大,當(dāng)信號(hào)主要由一階和二階信號(hào)組成時(shí),在做處理時(shí)求取一階、二階、三階項(xiàng)系數(shù)及四階項(xiàng)系數(shù),四階項(xiàng)即為高階項(xiàng)[1-2](此處不在超小周期做應(yīng)用,如在毫秒級(jí)或微秒級(jí),只取相似值,不選擇Δt->0)。
當(dāng)一、二階項(xiàng)系數(shù)發(fā)生較大改變或高階項(xiàng)系數(shù)不再為零,就是信號(hào)異常在特征值上的體現(xiàn)或此時(shí)系統(tǒng)行為不正常。當(dāng)信號(hào)特征體現(xiàn)出震蕩特性時(shí),經(jīng)過(guò)一階微分處理后會(huì)出現(xiàn)0值(極點(diǎn)),應(yīng)考慮是否判定為異常。
一階項(xiàng)的函數(shù)特征,當(dāng)信號(hào)正常時(shí),特征函數(shù)為:
(3)
一旦穩(wěn)定,c即為常量,對(duì)時(shí)間t分別求微分后得到:
(4)
其變化率即一階項(xiàng)系數(shù)為定值或當(dāng)穩(wěn)態(tài)時(shí)為0。
當(dāng)測(cè)量到的信號(hào)如圖1所示:
圖1 信號(hào)變化的斜率改變Fig.1 Slope change of signal change
其函數(shù)特征為:
(5)
對(duì)此信號(hào)特征做一階微分處理后,在其拐點(diǎn)處倒推函數(shù)的一階項(xiàng)系數(shù)會(huì)發(fā)生突然改變a->b,此時(shí)即可判定為信號(hào)特征發(fā)生改變[3]。
含有二階項(xiàng)函數(shù)特征,當(dāng)信號(hào)正常時(shí),特征函數(shù)為:
(6)
對(duì)時(shí)間t分別求二階微分后得到:
(7)
即二階項(xiàng)系數(shù)為定值或變化較小。那么當(dāng)其二次項(xiàng)系數(shù)變化較大且高階項(xiàng)系數(shù)出現(xiàn)不可忽略的數(shù)值,不再為零時(shí),即應(yīng)定義為異常狀態(tài)。由于水輪機(jī)控制是一個(gè)非線性的、非最小相位系統(tǒng),單獨(dú)研究二階項(xiàng)很復(fù)雜,故需要配合高階項(xiàng)系數(shù)判斷[4]。
震蕩特征的傳遞函數(shù)為:
(8)
ωn為自然角頻率;ξ為阻尼比;s為算子。
震蕩特征的阻尼系數(shù)較大時(shí),波形會(huì)很快收斂,在水輪機(jī)調(diào)節(jié)系統(tǒng)中屬于正常的調(diào)節(jié)形態(tài),對(duì)于信號(hào)分析的意義不大,如圖2。
圖2 大阻尼系數(shù)響應(yīng)Fig.2 Large damping coefficient response
只有當(dāng)阻尼系數(shù)較小時(shí)才會(huì)體現(xiàn)出震蕩特性,此種情況有分析價(jià)值,如圖3。
當(dāng)符合一階項(xiàng)函數(shù)特性的信號(hào)值出現(xiàn)震蕩特性時(shí),從診斷角度并不關(guān)心其弧度多少,也不關(guān)心其在交替之間的直線是否平滑,故震蕩特性中的一階項(xiàng)系數(shù)與二階項(xiàng)系數(shù)診斷意義不大,而高階項(xiàng)系數(shù)可以參考。
震蕩收斂還是震蕩發(fā)散對(duì)于水輪機(jī)控制系統(tǒng)很重要,如果對(duì)震蕩特性做一階微分處理后,微分值=0時(shí)會(huì)出現(xiàn)上、下極點(diǎn)。當(dāng)一個(gè)方向只出現(xiàn)一個(gè)極點(diǎn),在調(diào)速器系統(tǒng)的工程量中意義不大。當(dāng)同一方向出現(xiàn)一個(gè)以上的極點(diǎn)時(shí),應(yīng)判定為異常。需要指出的是,當(dāng)震蕩特征出現(xiàn)后,有可能會(huì)伴隨二階項(xiàng)和高階項(xiàng)系數(shù)不為零。但二階項(xiàng)和高階項(xiàng)系數(shù)不為零出現(xiàn)時(shí),不一定是震蕩。
正常的導(dǎo)葉開(kāi)度信號(hào)應(yīng)符合一階項(xiàng)的函數(shù)特征,即:
(9)
穩(wěn)態(tài)時(shí),c應(yīng)為常量。這里-a和b分別取接力器的最大開(kāi)啟速度(一般收主配開(kāi)度限制)和最大關(guān)閉速度(節(jié)流器)。但由于接力器本身及工況不同,其函數(shù)特征具有一些特殊工況,應(yīng)分類以下幾種情況分別處理:
1)二階項(xiàng)和高階項(xiàng)系數(shù)不可忽略為零,此種特征不符合接力器運(yùn)動(dòng)的本質(zhì)規(guī)律,當(dāng)判定為異常。
2)接力器本身帶有容性特性,傳遞函數(shù)為:
(10)
Ty為接力器時(shí)間常數(shù),s為算子。當(dāng)接力器開(kāi)始動(dòng)作時(shí),微分處理后的二階項(xiàng)系數(shù)會(huì)出現(xiàn)數(shù)值隨后消失。故較小的二階項(xiàng)系數(shù)應(yīng)當(dāng)視為正常現(xiàn)象。當(dāng)二階項(xiàng)系數(shù)大于0.3或出現(xiàn)高次項(xiàng)不為0時(shí),應(yīng)判定為異常。判斷條件0.3作為閾值在軟件中的設(shè)計(jì)參數(shù),便于調(diào)試期間修改。
3)接力器緩沖裝置。對(duì)于帶有緩沖裝置的接力器,在進(jìn)入緩沖段時(shí),微分處理后的二階項(xiàng)系數(shù)有可能出現(xiàn)較大的數(shù)值。這是緩沖裝置的固有特性,沒(méi)有意義。故將緩沖段范圍在軟件中設(shè)計(jì)為參數(shù),忽略其緩沖段的二階項(xiàng)系數(shù)。
4)一次調(diào)頻。當(dāng)轉(zhuǎn)速超出人工失靈區(qū)時(shí),調(diào)速器會(huì)進(jìn)入一次調(diào)頻狀態(tài)。在一次調(diào)頻接近預(yù)定值時(shí),PID波形特征會(huì)出現(xiàn)慣性特征。此時(shí)會(huì)出現(xiàn)二次項(xiàng)系數(shù),而波形的二次項(xiàng)系數(shù)會(huì)隨著調(diào)速器PI參數(shù)的變化而體現(xiàn)出不同特性,故應(yīng)當(dāng)忽略。當(dāng)一次調(diào)頻動(dòng)作時(shí),應(yīng)將頻率實(shí)時(shí)值帶入增量PID算式中求出理論值,并校驗(yàn)其水輪機(jī)控制系統(tǒng)的PID輸出。
5)空載和發(fā)電穩(wěn)態(tài)??蛰d并網(wǎng)和發(fā)電穩(wěn)態(tài)時(shí),調(diào)速器處于PID的閉環(huán)反饋控制,此時(shí)的接力器開(kāi)度會(huì)出現(xiàn)震蕩特征,同方向會(huì)出現(xiàn)一個(gè)以上的極點(diǎn)。當(dāng)進(jìn)入穩(wěn)態(tài)后,極點(diǎn)絕對(duì)差值A(chǔ)-B小于參數(shù)設(shè)定時(shí),應(yīng)視為正常現(xiàn)象。當(dāng)同方向極點(diǎn)大于2個(gè)且極點(diǎn)絕對(duì)差值A(chǔ)-B大于2%(在軟件中設(shè)計(jì)為參數(shù))時(shí),表示空載PID參數(shù)不好或出力出現(xiàn)震蕩。當(dāng)同方向極點(diǎn)大于2個(gè)且相同方向的極點(diǎn)差值A(chǔ)′-A為正,且大于2次(在軟件中設(shè)計(jì)為參數(shù)),表示PID參數(shù)的PI環(huán)節(jié)增益太大,有可能無(wú)法收斂,此時(shí)應(yīng)當(dāng)預(yù)警。
6)分段關(guān)閉。當(dāng)主機(jī)的過(guò)度過(guò)程有分段關(guān)閉需求時(shí),應(yīng)關(guān)注其拐點(diǎn)的接力器位置是否有較大偏差。
主配開(kāi)度特征主要受伺服控制中的閉環(huán)反饋K0和K1影響,如圖4所示。
圖4 伺服控制邏輯Fig.4 Servo control logic
閉環(huán)反饋差值的特征變化復(fù)雜,一階項(xiàng)和二階項(xiàng)系數(shù)沒(méi)有診斷意義。研究主配開(kāi)度特征形態(tài)主要有以下幾種情況:①震蕩特征。主配開(kāi)度不能出現(xiàn)震蕩特征,且主配抽動(dòng)對(duì)于機(jī)組的傷害較大,故體現(xiàn)出震蕩特征時(shí)應(yīng)預(yù)警。②頻繁出現(xiàn)高階項(xiàng)不為零。此種情況是主配閥芯卡澀的特征之一,應(yīng)預(yù)警并告知維護(hù)人員檢查。③隨動(dòng)特征。應(yīng)實(shí)時(shí)監(jiān)測(cè)主配開(kāi)度是否長(zhǎng)期偏離給定大于5%(在軟件中設(shè)計(jì)為參數(shù)),如果滿足則表示伺服隨動(dòng)的特性較差,應(yīng)當(dāng)預(yù)警。主配開(kāi)度-[(接力器給定-接力器反饋)xK0-主配開(kāi)度]>=5%。④主配拒動(dòng)。當(dāng)過(guò)速115%且主配不在全關(guān)時(shí),1 s后觸發(fā)主配拒動(dòng)預(yù)警。
水頭信號(hào)受壓力鋼管壓力的影響,特征滿足剛性水擊模型特征,體現(xiàn)為導(dǎo)葉快速變化出現(xiàn)功率反調(diào)及甩負(fù)荷時(shí)引水鋼管升壓。剛性水擊傳遞函數(shù)如下:
(11)
Y為接力器行程,Tw為水擊時(shí)間常數(shù),s為算子。
除以上情況外,水頭的正常特征是在正常水頭范圍內(nèi)緩慢變化。研究水頭信號(hào)主要有以下兩種異常特征形態(tài):①震蕩特征。當(dāng)水頭出現(xiàn)震蕩特征時(shí)為異?,F(xiàn)象,應(yīng)預(yù)警。在一些有較長(zhǎng)引水鋼管的機(jī)組,水頭的震蕩特征會(huì)有其他體現(xiàn),如球閥自激振現(xiàn)象,從而導(dǎo)致水壓呈震蕩波動(dòng)。二階項(xiàng)和高階項(xiàng)系數(shù)不可忽略,在水頭信號(hào)上此種現(xiàn)象可能會(huì)伴隨震蕩特征,也可能是傳感器問(wèn)題。②階躍。一旦出現(xiàn)階躍,表示信號(hào)本身品質(zhì)差,應(yīng)重點(diǎn)檢查傳感器。
正常的系統(tǒng)壓力只符合一階函數(shù)特征,故對(duì)其特征分析應(yīng)包含二次項(xiàng)和高次項(xiàng)系數(shù)。此種現(xiàn)象一般為傳感器問(wèn)題或存在信號(hào)干擾。震蕩特征與時(shí)間間隔在軟件中設(shè)計(jì)為參數(shù),大周期震蕩時(shí)應(yīng)忽略。如果周期小于10 min,應(yīng)檢查是否壓力管路存在漏油。階躍表示信號(hào)本身品質(zhì)差,可能存在干擾。對(duì)于系統(tǒng)壓力來(lái)說(shuō),一階函數(shù)特征以外的其他特征一旦是真實(shí)反應(yīng)測(cè)量值,故障程度較嚴(yán)重,應(yīng)及時(shí)排查。
根據(jù)多階差分原理提出倒推多項(xiàng)式各階項(xiàng)系數(shù),以此為特征值作為信號(hào)診斷的判據(jù)。對(duì)水輪機(jī)控制系統(tǒng)的重要信號(hào)進(jìn)行分析,根據(jù)工程量的自身特性來(lái)定義某信號(hào)屬于幾階信號(hào),結(jié)合特征值的算法識(shí)別信號(hào)特征及形態(tài),從而實(shí)現(xiàn)信號(hào)診斷。在生產(chǎn)過(guò)程中,很多控制系統(tǒng)無(wú)法自動(dòng)識(shí)別異常,而維護(hù)人員則能判斷出有問(wèn)題,這是因?yàn)闇y(cè)量值在允許范圍內(nèi),但形態(tài)不對(duì),故研究典型的異常特征形態(tài)并不斷完善,具有現(xiàn)實(shí)意義。