雷世雄，徐非，譚翔，彭敏俊，

1.海軍駐武漢七一九所軍事代表室，湖北武漢 430000 2.哈爾濱工程大學(xué)核科學(xué)與技術(shù)學(xué)院，黑龍江哈爾濱 150001

分布式策略的核動(dòng)力裝置故障診斷技術(shù)研究

雷世雄1，徐非2，譚翔2，彭敏俊1，2

1.海軍駐武漢七一九所軍事代表室，湖北武漢 430000 2.哈爾濱工程大學(xué)核科學(xué)與技術(shù)學(xué)院，黑龍江哈爾濱 150001

對(duì)核動(dòng)力裝置中的故障問題進(jìn)行了研究，根據(jù)系統(tǒng)結(jié)構(gòu)和運(yùn)行特點(diǎn)，基于分布式診斷策略將核動(dòng)力裝置劃分成三級(jí)子單元結(jié)構(gòu)，將故障診斷任務(wù)劃分到不同的子單元中，采用核主元分析（kernel principal component analysis）進(jìn)行故障檢測(cè)，監(jiān)測(cè)核動(dòng)力裝置的運(yùn)行狀態(tài)，采用模型法對(duì)不同的設(shè)備單元建立相應(yīng)的模型進(jìn)行故障診斷、識(shí)別，并輔以專家知識(shí)進(jìn)行補(bǔ)充，得到故障問題的最優(yōu)解，并以此為依據(jù)設(shè)計(jì)了核動(dòng)力裝置分布式故障診斷系統(tǒng)。在全范圍仿真機(jī)上的測(cè)試表明，基于分布式策略的故障診斷方法有效，核動(dòng)力裝置分布式故障診斷系統(tǒng)具有較高的可靠度，滿足故障診斷的功能要求。

分布式策略；故障診斷；診斷模型；專家知識(shí)；核主元分析

故障診斷技術(shù)是核動(dòng)力裝置運(yùn)行支持技術(shù)的重要組成部分之一，故障診斷技術(shù)的研究對(duì)核動(dòng)力裝置的安全性和經(jīng)濟(jì)性有著重要的意義，能夠增強(qiáng)核動(dòng)力裝置的安全性能，減少操作員誤操作的概率。目前，故障診斷方法大致可以分為兩類：基于數(shù)據(jù)驅(qū)動(dòng)的經(jīng)驗(yàn)?zāi)Ｐ头ê突跀?shù)學(xué)方程的物理模型法［1］。這2種方法各有優(yōu)缺點(diǎn)。物理模型法可以對(duì)診斷對(duì)象進(jìn)行本質(zhì)描述，給操縱員提供更顯著的信息，隨著對(duì)設(shè)備運(yùn)行狀況的更深的了解，可以逐漸修正監(jiān)測(cè)模型以提高精度［2］。而基于數(shù)學(xué)驅(qū)動(dòng)的診斷方法只需要被診斷對(duì)象的歷史數(shù)據(jù)和在線數(shù)據(jù)，通過對(duì)數(shù)據(jù)分析來判別對(duì)其運(yùn)行狀態(tài)。本文將數(shù)據(jù)驅(qū)動(dòng)方法與物理模型法相結(jié)合建立核動(dòng)力裝置分布式故障診斷系統(tǒng)，通過兩者的結(jié)合，使故障診斷結(jié)果可靠，可理解性更好，診斷速度快。

1 分布式故障診斷策略

基于分布式策略的診斷方法就是將整個(gè)系統(tǒng)的監(jiān)測(cè)診斷任務(wù)基于一定的原則分解為各個(gè)子單元下較為簡單的監(jiān)測(cè)診斷任務(wù)，使系統(tǒng)的診斷任務(wù)可以由若干個(gè)不同子單元共同完成，以降低故障問題的診斷難度［3］。如圖1，其設(shè)計(jì)思想主要包括以下3個(gè)方面。

圖1 分布式診斷策略的設(shè)計(jì)思想

1）對(duì)多故障診斷任務(wù)A的分配與分解。根據(jù)待診斷對(duì)象在結(jié)構(gòu)、機(jī)理和功能上的分布性和層次性特點(diǎn)，將診斷任務(wù)A分解成不同子單元下的子任務(wù)An；

2）對(duì)子任務(wù)An的求解。根據(jù)每個(gè)診斷子單元自身的特點(diǎn)，采用適合的診斷方法或領(lǐng)域知識(shí)對(duì)子任務(wù)An進(jìn)行求解，得到子任務(wù)的解Xn；

3）綜合求解。將子任務(wù)的求解結(jié)果Xn根據(jù)一定的原則進(jìn)行結(jié)論融合，排除矛盾的結(jié)論，形成系統(tǒng)診斷任務(wù)A的最終解X。

為了便于分析研究，基于核動(dòng)力裝置的系統(tǒng)功能和系統(tǒng)結(jié)構(gòu)，將分布式故障診斷系統(tǒng)的任務(wù)分解到三級(jí)結(jié)構(gòu)單元中去，本文主要針對(duì)一回路系統(tǒng)中的4個(gè)主要系統(tǒng)即主冷卻劑系統(tǒng)（reactor coolant system，RCS）、化學(xué)與容積控制系統(tǒng)（chemical and volume control system，CVCS）、設(shè)備冷卻水系統(tǒng)（component cooling water system，CCWS）及余熱排出系統(tǒng)（afterheat removal system，RHRS）進(jìn)行研究分析，如圖2所示。

圖2 分布式多故障診斷系統(tǒng)的三級(jí)結(jié)構(gòu)

2 基于KPCA的故障檢測(cè)

本文針對(duì)第一級(jí)結(jié)構(gòu)單元即四個(gè)系統(tǒng)運(yùn)用KP-CA方法進(jìn)行故障檢測(cè)，故障檢測(cè)是故障診斷的重要組成部分，是進(jìn)行故障辨識(shí)的入口條件。

2.1 KPCA建模

已知數(shù)據(jù)樣本矩陣Xn×m由n組觀測(cè)樣本、m個(gè)測(cè)量變量（m≥n）組成，即

式中：xi∈Rm( i＝1，2，…，n)。將Xn×m通過非線性映射Φ從原始數(shù)據(jù)空間U映射到高維特征空間F上［4］，得到一個(gè)n×n的對(duì)稱正定核矩陣Kn×n，核矩陣中的元素Kij為

選定核函數(shù)形式后，對(duì)K進(jìn)行高維空間上的中心化：

式中：1n為元素為1的n階常數(shù)矩陣。

非線性空間F中核矩陣Kn×n的協(xié)方差矩陣C為

協(xié)方差矩陣C的特征方程：

式中：λ為特征值，v為特征向量。由式（5）左乘(xk)得：

式中特征向量v可由特征空間的樣本映射成：

由式（4）～式（6）可得：

式（8）可簡化為

由此可得到樣本X在特征空間F中的第k個(gè)核主元tk：

2.2 主元選取

本文采用CPV方法作為主元選取的準(zhǔn)則。CPV方法反映了前a個(gè)核主元確定的模型包含原數(shù)據(jù)信息的程度：

式中：n是核矩陣K的維度，λi是協(xié)方差矩陣C的第i個(gè)特征值。當(dāng)累計(jì)貢獻(xiàn)率達(dá)到90%以上的時(shí)候，認(rèn)為核主元個(gè)數(shù)包含了原數(shù)據(jù)足夠多的信息。

2.3 故障檢測(cè)

平方預(yù)測(cè)誤差（squared prediction error，SPE）和Hotelling’s T2是在基于KPCA的檢測(cè)過程中常用的指標(biāo)［5-7］，以檢測(cè)過程是否發(fā)生異常，這2個(gè)指標(biāo)是互補(bǔ)的。

1）SPE統(tǒng)計(jì)量

SPE統(tǒng)計(jì)量描述的是實(shí)際運(yùn)行過程與樣本模型之間的偏離程度，可以衡量樣本向量在殘差空間投影的變化：

式中：I為n階單位矩陣，V∈Rn為協(xié)方差矩陣C的前a個(gè)特征向量構(gòu)成的負(fù)載矩陣，δα2表示置信水平為α?xí)r的控制限，該控制限計(jì)算公式如下：cα為標(biāo)準(zhǔn)正態(tài)分布在置信水平α為下的閾值。

2）T2統(tǒng)計(jì)量

T2統(tǒng)計(jì)量主要用于衡量主元空間中監(jiān)測(cè)參數(shù)變量的變化程度：

式中：Fa，n－a；α是帶有a和n－a個(gè)自由度、置信水平為α的F分布臨界值。

3 故障辨識(shí)

核動(dòng)力裝置中，各設(shè)備與周圍設(shè)備之間存在傳質(zhì)、傳熱，在熱力學(xué)上視為開口系統(tǒng)，滿足最基本的物理定律和原理，本文運(yùn)用基本物理方程對(duì)各個(gè)設(shè)備建立模型，在故障檢測(cè)的基礎(chǔ)上，識(shí)別出發(fā)生故障的設(shè)備單元。

3.1 單元模型

3.1.1 換熱器單元

根據(jù)換熱器的質(zhì)量能量守恒以及傳熱公式，可得如下公式：

式中：1、2、t表示分別表示熱交換器管側(cè)、殼側(cè)、管壁；cp表示熱交換器內(nèi)工質(zhì)定壓比熱容，kJ／（kg·K）；m表示熱交換器內(nèi)工質(zhì)質(zhì)量，表示熱交換器內(nèi)工質(zhì)平均密度，kg／m3；V表示熱交換器內(nèi)工質(zhì)體積，m3；Tavg表示熱交換器內(nèi)工質(zhì)平均溫度，K；Q表示熱交換器管側(cè)導(dǎo)出熱量，kJ；Qr表示熱交換器對(duì)外散熱，kJ。

當(dāng)換熱器內(nèi)發(fā)生泄漏故障時(shí)，假設(shè)泄漏是由管側(cè)向殼側(cè)發(fā)生，則式（17）應(yīng)改寫為如下形式：由式（18）可知，熱交換器在發(fā)生管側(cè)向殼側(cè)泄漏時(shí)，管側(cè)出口流量的實(shí)際值Gout1比按正常運(yùn)行計(jì)算得到的理論值Gout1小，且泄漏量可按式（18）計(jì)算。

3.1.2 離心泵單元

由伯努利方程可得，離心泵的揚(yáng)程為

與式（17）相比較可見，式（18）多出了一個(gè)含Gleak的項(xiàng)：

根據(jù)上述兩式以及相似關(guān)系式可得：

對(duì)于同一個(gè)離心泵，當(dāng)處于正常運(yùn)行時(shí)，系數(shù)K為恒定的常數(shù)；當(dāng)發(fā)生運(yùn)行故障停轉(zhuǎn)時(shí)，系數(shù)K也隨之發(fā)生改變，其進(jìn)、出口壓差迅速變?yōu)樨?fù)壓。

3.1.3 熱阱單元

在正常運(yùn)行時(shí)，熱阱中的泄漏項(xiàng)均為0，熱阱液位隨時(shí)間t積分可得：

根據(jù)泵的相似理論，對(duì)于同一臺(tái)泵，認(rèn)為其幾何相似、運(yùn)動(dòng)相似［8］。

對(duì)于任意時(shí)刻t和t＋Δt，定義一個(gè)離心泵的特征系數(shù)K，K應(yīng)為常數(shù)：

式中：′、″分別代表液相、氣相；in、out、leak表示分別代表單元進(jìn)、出口及泄漏處；A表示單元截面積，m2；H表示單元內(nèi)的液位，m；H0表示單元總高度，m。

而當(dāng)熱阱發(fā)生泄漏故障時(shí)，泄漏項(xiàng)不為0，式（22）應(yīng)改寫為可見，在泄漏故障時(shí)，熱阱液位的實(shí)際值Ht＋Δt比按正常運(yùn)行計(jì)算得到的理論值Ht＋Δt低。對(duì)熱阱子單元來說，其橫截面積A、高度H0只與結(jié)構(gòu)有關(guān)，液位H、各進(jìn)出口管線流量G可直接監(jiān)測(cè)，而密度等熱工參數(shù)可由監(jiān)測(cè)壓力、溫度查表求得。因此，可通過計(jì)算t＋Δt時(shí)刻的理論液位Ht＋Δt與實(shí)際監(jiān)測(cè)液位Ht＋Δt進(jìn)行比較，從而得到熱阱泄漏故障的診斷結(jié)果。

3.2 基于知識(shí)的模型

對(duì)于核動(dòng)力裝置來說，有些單元難以根據(jù)物理原理建立相應(yīng)的數(shù)學(xué)過程模型；有些單元雖然可以建模，但存在更為簡便的診斷方法。對(duì)于這2類單元，本文中采用基于專家知識(shí)的模型進(jìn)行補(bǔ)充診斷。專家知識(shí)的表示方法采用產(chǎn)生式規(guī)則，一般形式為

Rule：IF＜條件＞THEN＜動(dòng)作＞（24）例如，反應(yīng)堆運(yùn)行過程中，冷卻劑和慢化劑經(jīng)高能中子照射后將發(fā)生16O（n，p）16N反應(yīng)，該反應(yīng)只存在與反應(yīng)堆中，并且17N的半衰期約為4.11 s，很快就衰變?yōu)?6N。16N是堆內(nèi)水放射性的重要來源，基本不存在于二回路中。因此，可以通過監(jiān)測(cè)二次側(cè)工質(zhì)16N放射性來作為蒸汽發(fā)生器傳熱管是否發(fā)生泄漏的判斷依據(jù)。

4 結(jié)果測(cè)試

核動(dòng)力裝置分布式故障診斷系統(tǒng)（NPPDFDS）的工作流程如圖3所示。NPPDMFDS主要包括：數(shù)據(jù)處理、工況識(shí)別、故障監(jiān)測(cè)、故障識(shí)別4個(gè)模塊。其中，數(shù)據(jù)處理模塊先將采集而來的數(shù)據(jù)進(jìn)行預(yù)處理，以去除噪聲，有利于后續(xù)故障診斷。工況識(shí)別模塊式識(shí)別當(dāng)前的運(yùn)行工況，本文在100%功率運(yùn)行工況下進(jìn)行分析。故障檢測(cè)模塊利用KPCA作為故障檢測(cè)的方法，對(duì)當(dāng)前工況下的運(yùn)行狀態(tài)進(jìn)行監(jiān)測(cè)，以判斷核動(dòng)力裝置是否正常運(yùn)行，是否有異常發(fā)生。故障識(shí)別運(yùn)用機(jī)理模型對(duì)各個(gè)設(shè)備子單元進(jìn)行診斷分析，從而得出導(dǎo)致異常的設(shè)備。其中KPCA故障檢測(cè)模型，根據(jù)CPV方法確定第一級(jí)結(jié)構(gòu)單元的核主元個(gè)數(shù)，當(dāng)前a個(gè)主元的累積貢獻(xiàn)率達(dá)到90%以上的時(shí)候，即認(rèn)為選取的核主元包含足夠多的原數(shù)據(jù)信息。如圖4所示。其中，主冷卻劑系統(tǒng)（RCS）監(jiān)測(cè)參數(shù)有140個(gè)，選取的核主元個(gè)數(shù)為7?；瘜W(xué)容積控制系統(tǒng)（CVCS）監(jiān)測(cè)參數(shù)有68個(gè)，取3個(gè)核主元。設(shè)備冷卻水（CCWS）系統(tǒng)的監(jiān)測(cè)參數(shù)有15個(gè)，取2個(gè)主元；余熱排出系統(tǒng)（RHRS）監(jiān)測(cè)參數(shù)有19個(gè)，該系統(tǒng)在功率穩(wěn)態(tài)運(yùn)行時(shí)處于被隔離狀態(tài)，系統(tǒng)參數(shù)不發(fā)生變化，取1個(gè)核主元。

圖3 NPPDMFDS的工作過程

圖4 累計(jì)貢獻(xiàn)率法確定主元個(gè)數(shù)

取置信水位α為0.05，各系統(tǒng)在100%功率穩(wěn)態(tài)運(yùn)行時(shí)的控制限如表1所示。

表1 PCA模型中各監(jiān)測(cè)子單元的控制限

4.1 單一故障測(cè)試

在全范圍仿真機(jī)上100%功率穩(wěn)態(tài)運(yùn)行工況下，運(yùn)行核動(dòng)力裝置分布式故障診斷系統(tǒng)。20 s時(shí)在1＃蒸汽發(fā)生器中插入故障度為1%的蒸汽發(fā)生器傳熱管破裂事故（SGTR），SGTR是核動(dòng)力裝置中最常見的故障之一。SPE統(tǒng)計(jì)量示意和T2統(tǒng)計(jì)量示意如圖5、6。

圖5 SPE統(tǒng)計(jì)量示意圖

圖6 T2統(tǒng)計(jì)量示意圖

如圖5、6所示，插入SGTR故障后，KPCA故障檢測(cè)模塊對(duì)數(shù)據(jù)處理模塊處理后的數(shù)據(jù)進(jìn)行檢測(cè)，22 s時(shí)主冷卻劑系統(tǒng)的SPE統(tǒng)計(jì)量超出控制限，23 s時(shí)T2統(tǒng)計(jì)量超出控制限。當(dāng)SPE和T2同時(shí)超出控制限時(shí)，說明核動(dòng)力裝置運(yùn)行異常，NPPDMFDS 在23 s時(shí)診斷出系統(tǒng)異常，發(fā)出報(bào)警。此時(shí)，將故障檢測(cè)模塊得出的結(jié)果作為故障識(shí)別模塊的入口觸發(fā)條件，進(jìn)入故障辨識(shí)模塊。

各個(gè)設(shè)備的機(jī)理模型接收其所需的數(shù)據(jù)進(jìn)行故障識(shí)別，若數(shù)據(jù)代入機(jī)理模型仍滿足物理守恒定律或功能要求，那么該設(shè)備并沒有發(fā)生故障，若數(shù)據(jù)代入機(jī)理模型不滿組物理守恒定律或功能要求，則說明該設(shè)備發(fā)生故障。如圖7，蒸汽發(fā)生器的知識(shí)模型得出1＃SG蒸汽出口的16N放射性突然呈幾何增加，然后略有下降。在約7.13 s（即16N的半衰期）后，維持在一個(gè)穩(wěn)定值；而2＃SG蒸汽出口的16N放射性則基本不變。根據(jù)3.2節(jié)中的蒸汽發(fā)生器子單元模型，可得到“1＃SG中發(fā)生SGTR”的診斷結(jié)果。

圖7 1%SGTR狀態(tài)下蒸汽發(fā)生器出口16N放射性

4.2 故障組

為了證明NPPDFDS對(duì)多故障同樣具有診斷辨識(shí)功能，本文選取以下故障組作為實(shí)例。

故障組由再生熱交換器傳熱管泄漏和下泄熱交換器傳熱管泄漏2個(gè)故障組成，每個(gè)故障的故障度分別取5%。在100%功率穩(wěn)態(tài)運(yùn)行工況下，運(yùn)行NPPDFDS。20 s時(shí)先插入5%的再生熱交換器傳熱管泄漏，22 s時(shí)插入5%的下泄熱交換器傳熱管泄漏故障。如圖8、9所示，插入SGTR故障后，22 s時(shí)化學(xué)和容積控制系統(tǒng)的SPE統(tǒng)計(jì)量超出控制限，24 s 時(shí)T2統(tǒng)計(jì)量超出控制限。當(dāng)SPE和T2統(tǒng)計(jì)量同時(shí)超出控制限時(shí)，說明故障檢測(cè)模塊檢測(cè)出核動(dòng)力裝置運(yùn)行異常。NPPDFDS在24 s時(shí)診斷出系統(tǒng)異常，發(fā)出報(bào)警。此時(shí)，觸發(fā)故障辨識(shí)模塊進(jìn)入工作狀態(tài)。將數(shù)據(jù)傳輸?shù)礁鱾€(gè)設(shè)備機(jī)理模型中，其中下泄熱交換器和再生熱交換器的模型計(jì)算結(jié)果如圖10所示，以下泄熱交換器為對(duì)象進(jìn)行計(jì)算，下泄子單元出口量實(shí)際監(jiān)測(cè)值小于理論值。根據(jù)上述模型可知，下泄熱交換器傳熱管發(fā)生泄漏故障。

圖8 SPE統(tǒng)計(jì)量示意圖

圖9 T2統(tǒng)計(jì)量示意圖

圖10 下泄子單元出口流量

下泄子單元出口流量如圖11。

圖11 下泄子單元出口流量

以再生熱交換器為對(duì)象進(jìn)行計(jì)算，下泄子單元出口量實(shí)際監(jiān)測(cè)值小于理論值，根據(jù)熱交換器物理模型，再生熱交換器未發(fā)生故障，但考慮到下泄熱交換器泄漏的影響，對(duì)理論值進(jìn)行修正后，監(jiān)測(cè)值大于修正理論值，則再生熱交換器也發(fā)生傳熱管泄漏故障。因此該分布式故障診斷系統(tǒng)可以準(zhǔn)確的診斷出故障，識(shí)別出發(fā)生故障的設(shè)備單元。

5 結(jié)論

1）測(cè)試結(jié)果表明，采用基于分布式的故障診斷策略，基于KPCA的故障檢測(cè)方法可以迅速的檢測(cè)出系統(tǒng)的異常狀態(tài)，為故障診斷提供入口條件；基于物理模型和知識(shí)的模型可以充分發(fā)揮分布式診斷策略的優(yōu)勢(shì)，從而消除故障征兆影響，得到可靠地故障診斷結(jié)果

2）在進(jìn)一步完善的基礎(chǔ)上，基于分布式的故障診斷策略和基于物理模型及知識(shí)的診斷方法的研究可為核動(dòng)力裝置故障診斷系統(tǒng)提供一定的參考或借鑒作用。

［1］MA Jianping，JIANG Jin.Application of fault detection and diagnosismethods in nuclear power plants：a review［J］.Progress in Nuclear Energy，2011，53：255-266.

［2］劉繼偉.基于大數(shù)據(jù)的多尺度狀態(tài)監(jiān)測(cè)方法及其應(yīng)用［D］.北京：華北電力大學(xué)，2013：5-10.

［3］閆修平，彭敏俊，成守宇.分布式故障診斷策略及其在核電站中的應(yīng)用［J］.核動(dòng)力工程，2009，30（2）：99-103.

［4］倪國文，劉愛倫.基于核主元分析的非線性故障檢測(cè)研究［J］.儀器儀表學(xué)報(bào)，2009，30（6）：443-447.

［5］LEE JM，YOO C K，CHOISW，et al.Nonlinear process monitoring using kernel principal component analysis［J］.Chemical Engineering Science，2004，59（1）：223-234.

［6］LEE JM，YOO C K，LEE I B.Fault detection of batch processes usingmultiway kernel principal component analy-sis［J］.Computers and Chemical Engineering，2004，28 （9）：1837-1840.

［7］WANG Haiqing，SONG Zhihuan，WANG Hui.Fault detec-tion behavior and performance analysis of PCA-based processmonitoringmethods［J］.Journal of Chemical Indus-try and Engineering，2002，53（3）：297-300.

［8］高璞珍.核動(dòng)力裝置用泵［M］.哈爾濱：哈爾濱工程大學(xué)出版社，2004：20-25.

Research on fault diagnosis technology for a nuclear power unit based on the distributed strategy

LEIShixiong1，XU Fei2，TAN Xiang2，PENGMinjun1，2
1.Military Representative Office of Navy in 719 Institute，Wuhan 430000 2.College of Nuclear Science and Technology，Harbin Engineering University，Harbin 150001

This paper studies the fault problem in a nuclear power unit.According to the structure and operation characteristics of the system，it is divided into three-level subunits on the basis of the distributed strategy.Kernel principle component analysis is used as the fault detection method，and physicalmodels and expert knowledge that corresponds to each unit are used to diagnose faults so as to eliminate the influence of symptom coupling on diagno-sis results and get optimal solution of fault problem under abnormal conditions.A NPPDMFDS is developed on the basis of themethod mentioned above.The test on the full scope of simulator shows that the fault diagnosismethod based on the distributed strategy is effective and NPPDMFDS has high reliability，meeting the requirement of fault diagnosis.

distributed strategy；fault diagnosis；diagnosismodel；expert knowledge；kernel principal componenta-nalysis

10.3969／j.issn.1009-671X.201411004

TK05

A

1009-671X（2015）02-068-07

2014-11-09.

日期：2015-03-25.

雷世雄（1968-），男，高工.

彭敏俊，E-mail：heupmj＠hrbeu.edu.cn.

http：／／www.cnki.net／kcms／detail／23.1191.U.20150326.0846.002.html