本發(fā)明涉及一種故障診斷方法，特別適用于中重型柴油機尾氣后處理scr系統(tǒng)中對輸入nox傳感器的故障診斷。

背景技術：

scr系統(tǒng)被廣泛應用于柴油機中來減少nox(x＝1,2)的排放，被認為是最具有前途的去除氮氧化物的一種尾氣后處理技術。scr(selectivecatalyticreduction)主要就是一種選擇催化還原技術，它是在系統(tǒng)入口端噴入尿素，利用其水解出的氨氣在催化劑的作用下將nox還原為對大氣無污染的氮氣和水。

在scr系統(tǒng)中，尿素的噴射是系統(tǒng)唯一的輸入控制，而nox的濃度被認為是scr尿素噴射控制中重要的狀態(tài)參數(shù)，所以nox的輸入濃度對控制算法的設計至關重要?，F(xiàn)實中輸入濃度測量一般都是借助于傳感器，但是傳感器作為一種精細的裝置，是整個scr系統(tǒng)中最易發(fā)生故障的部件。此外柴油機在工作狀態(tài)時會產(chǎn)生強烈的振動，在這種惡劣的工作環(huán)境下傳感器也極其容易發(fā)生故障。因此開發(fā)一種scr系統(tǒng)nox傳感器故障診斷方法對提高scr系統(tǒng)的經(jīng)濟效益和工作效率十分重要。

技術實現(xiàn)要素：

針對現(xiàn)有技術中存在的問題，本發(fā)明提供一種柴油機scr系統(tǒng)nox(x＝1,2)輸入傳感器故障診斷方法，所述的故障診斷方法包括如下步驟：

第一步，設置柴油機scr系統(tǒng)nox輸入傳感器故障診斷系統(tǒng)，其包括nox輸入傳感器、nox輸入濃度估計模塊、nox傳感器故障診斷模塊；所述的nox輸入傳感器和nox輸入濃度估計模塊分別通過線纜連接到nox傳感器故障診斷模塊。

第二步，根據(jù)scr系統(tǒng)參數(shù)，設置nox輸入濃度閾值δδ(y)和時間閾值t(y)；

第三步，nox輸入傳感器測量nox輸入濃度實際值δ0，并將nox輸入濃度實際值δ0的信號輸送至nox傳感器故障診斷模塊，nox輸入濃度估計模塊估算nox輸入濃度估計值并將nox輸入濃度估計值的信號輸送至nox傳感器故障診斷模塊，nox傳感器故障診斷模塊根據(jù)接收到的信號計算nox輸入濃度殘差δδ，計算公式如下：

第四步，比較nox輸入濃度殘差δδ和nox輸入濃度閾值δδ(y)，若δδ≤δδ(y)，則進入第五步；若δδ>δδ(y)，則繼續(xù)測量nox輸入濃度殘差δδ大于nox輸入濃度閾值δδ(y)持續(xù)的時間t1，若t1<t(y)，則進入第五步；否則進入第六步；

第五步，nox傳感器故障診斷模塊確定nox輸入傳感器正常，并將nox輸入濃度實際值δ0和nox輸入濃度估計值輸送至scr系統(tǒng)控制單元，控制單元輸出最終的nox輸入濃度值

第六步，nox傳感器故障診斷模塊確定nox輸入傳感器發(fā)生故障，并將nox輸入濃度實際值δ0和nox輸入濃度估計值輸送至scr系統(tǒng)控制單元，控制單元將估計值作為最終的nox輸入濃度值，并發(fā)出傳感器故障警報，停止診斷。

本發(fā)明的優(yōu)點在于：

1、本發(fā)明運用無跡卡爾曼濾波來估算nox輸入濃度的估計值，適用于scr非線性系統(tǒng)，具有較高的診斷精度。

2、采用的無跡卡爾曼濾波估算方法計算簡單，響應速度快，能夠實時在線對scr系統(tǒng)傳感器進行故障診斷。本發(fā)明可以廣泛應用于柴油機scr系統(tǒng)傳感器故障診斷中。

附圖說明

圖1是scr系統(tǒng)及相應傳感器示意圖。

圖2是本發(fā)明的nox傳感器故障診斷實時在線診斷過程示意圖。

具體實施方式

下面結合附圖對本發(fā)明的一種柴油機scr系統(tǒng)氮氧化物輸入傳感器故障診斷方法進行詳細的描述。

如圖1所示，本發(fā)明所述的scr系統(tǒng)中，scr裝置的輸入管路上設置了nox輸入傳感器和nh3輸入傳感器，輸出管路上設置了nox輸出傳感器和nh3輸出傳感器，nox輸入傳感器測量nox輸入濃度，nh3輸入傳感器測量nh3輸入濃度，nox輸出傳感器測量nox輸出濃度，nh3輸出傳感器測量nh3輸出濃度。本發(fā)明提供的故障診斷方法旨在對所述的nox輸入傳感器進行故障診斷。如圖2所示，本發(fā)明提供了一種柴油機scr系統(tǒng)nox輸入傳感器故障診斷方法，包括以下步驟：

第一步，設置柴油機scr系統(tǒng)輸入nox傳感器故障診斷系統(tǒng)，其包括nox輸入傳感器、nox輸入濃度估計模塊、nox傳感器故障診斷模塊；

第二步，根據(jù)scr系統(tǒng)參數(shù)，設置nox輸入濃度閾值δδ(y)和時間閾值t(y)；

第三步，nox輸入傳感器測量nox輸入濃度實際值δ0，并將nox輸入濃度實際值δ0的信號輸送至nox傳感器故障診斷模塊，nox輸入濃度估計模塊估算nox輸入濃度估計值并將nox輸入濃度估計值的信號輸送至nox傳感器故障診斷模塊，nox傳感器故障診斷模塊根據(jù)接收到的信號計算nox輸入濃度殘差δδ，計算公式如下：

其中，nox輸入濃度估計值的計算方法如下：

(1)根據(jù)scr系統(tǒng)內化學反應，及摩爾守恒和質量守恒建立狀態(tài)空間模型，

其中，rm表示相應反應過程的反應速率，m＝ad，de，ox，re；ad代表氨氣吸附過程反應速率；de代表氨氣解吸附過程反應速率；ox代表氨氣氧化反應速率；re代表氨氣還原nox反應速率；cno和是nox(x＝1,2)和氨氣的濃度；表示氨氣入口的濃度，cno,in表示柴油機排出尾氣中nox的濃度；km表示相應反應過程的反應系數(shù)，em表示相應反應過程的化學反應活化能；r為反應常數(shù)；t代表溫度；f是尾氣流速；v是scr裝置的體積；是scr罐內氨氣的覆蓋率；和分別為氮氧化物濃度的變化率、氨氣濃度的變化率和scr罐內氨氣的覆蓋率的變化率。

(2)將狀態(tài)空間模型與無跡卡爾曼濾波(ukf)方程結合起來。ukf估計狀態(tài)一般分為兩步，包括預測和更新。非線性系統(tǒng)的ukf方程一般表示如下：

x(k)＝f[x(k-1),u(k)]+w(k)(2)

z(k)＝h[x(k)]+v(k)，(3)

方程(2)是預測方程。x(k)代表狀態(tài)向量，u(k)是輸入向量，w(k)代表高斯過程噪聲，f()是預測函數(shù)。方程(3)是更新方程，也叫觀測方程，它包括觀測向量z(k)，觀測函數(shù)h(x)和高斯觀測噪聲v(k)。

其中預測過程包括以下幾個步驟：

a)構造sigma點；

在k-1步，根據(jù)隨機狀態(tài)變量x的統(tǒng)計量和協(xié)方差pk-1構造sigma點集。

其中λ是尺度參數(shù)，λ＝α²(nx+q)-nx。nx為狀態(tài)空間維數(shù)，在本發(fā)明中是3。

q是第二個尺度參數(shù)，一般取0或者3-nx。α設定為一個很小的常數(shù)，在本發(fā)明中取0.001。

b)對sigma點進行傳播計算

轉換公式如下，其中u代表輸入。

c)計算輸出先驗均值與誤差協(xié)方差；

輸出的先驗均值與誤差協(xié)方差pk|k-1計算公式計算如下：

其中q為噪聲協(xié)方差，和為計算先驗均值和誤差協(xié)方差的加權，定義如下：

其中β為常數(shù)，在高斯分布中，為獲得最佳估計一般取2。

更新過程通過如下步驟實現(xiàn)：

a)構造sigma點

根據(jù)上面預測過程計算出的先驗均值，再次構造sigma點。

b)計算預測輸出

傳播計算每個sigma點

則預測輸出公式如下：

c)計算卡爾曼增益

在這個過程中同時需要一個最佳的卡爾曼增益kk。計算方程如下所示：

其中，是預測出輸出誤差協(xié)方差，是均值與預測輸出交叉協(xié)方差，r是噪聲協(xié)方差。

d)計算后驗狀態(tài)估計和后驗協(xié)方差

在第k步，根據(jù)輸出的測量值，可以計算出后驗的狀態(tài)和協(xié)方差。

其中yk表示第k步的實際測量值，表示估計的觀測值。

假設nh3輸入濃度、輸出濃度和nox輸出濃度可以由傳感器測得，這樣就可以通過ukf觀測器來估計nox輸入濃度。由于在短時間內nox輸入濃度變化并不是特別劇烈，因此也可以認為nox輸入濃度的微分為零，公式如下。

在這一部分中，由于要估計nox輸入濃度，因此選擇nox輸入濃度作為預測方程的一個狀態(tài)量，另一個狀態(tài)量依然選擇scr內氨氣覆蓋率公式表達如下。

其中δt表示aukf(自適應無跡卡爾曼濾波，adaptiveunscentedkalmanfilter)計算方程中每步迭代更新時間，可以根據(jù)實際情況設定，本發(fā)明中設為0.01s。

在建立觀測方程時，選擇傳感器測量的nh3輸出濃度和nox輸出濃度作為兩個狀態(tài)量,在此認為nox輸出傳感器測得nox輸出濃度足夠精確。觀測方程如下所示。

其中定義

根據(jù)上述建立的預測方程和觀測方程，經(jīng)過ukf迭代計算，便可以估計出nox輸入濃度

第四步，比較nox輸入濃度殘差δδ和nox輸入濃度閾值δδ(y)，若δδ≤δδ(y)，則進入第五步；若δδ>δδ(y)，則繼續(xù)測量nox輸入濃度殘差δδ大于nox輸入濃度閾值δδ(y)持續(xù)的時間t1，若t1<t(y)，則進入第五步；否則進入第六步；

第五步，nox傳感器故障診斷模塊確定nox輸入傳感器正常，并將nox輸入濃度實際值δ0和nox輸入濃度估計值輸送至scr系統(tǒng)控制單元，控制單元輸出最終的nox輸入濃度值

第六步，nox傳感器故障診斷模塊確定nox輸入傳感器發(fā)生故障，并將nox輸入濃度實際值δ0和nox輸入濃度估計值輸送至scr系統(tǒng)控制單元，控制單元將估計值作為最終的nox輸入濃度值，并發(fā)出傳感器故障警報，停止診斷。