基于h∞濾波的銠自給能探測(cè)器信號(hào)延遲消除方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及先進(jìn)堆芯測(cè)量系統(tǒng)（核反應(yīng)堆功率分布在線監(jiān)測(cè)系統(tǒng)）所用的堆內(nèi)銠 自給能中子探測(cè)器信號(hào)的處理技術(shù)，具體是基于濾波的銠自給能探測(cè)器信號(hào)延遲消除 方法。
【背景技術(shù)】
[0002] 用作先進(jìn)堆芯測(cè)量系統(tǒng)堆內(nèi)探測(cè)器的銠自給能中子探測(cè)器，其敏感材料銠與中子 反應(yīng)產(chǎn)生的次生核素發(fā)生0衰變產(chǎn)生電流，穩(wěn)態(tài)情況下該電流大小與所在位置通量成正 比，因此通過(guò)測(cè)量銠自給能探測(cè)器能夠推知其所在位置中子通量。由于該類(lèi)探測(cè)器電流主 要成分是由次生核素0衰變產(chǎn)生的，在反應(yīng)堆瞬態(tài)情況（中子通量水平變化的情況）下， 該類(lèi)探測(cè)器電流不能實(shí)時(shí)反映通量水平的變化，而是有一定的延遲，延遲時(shí)間參數(shù)與次生 核素的0衰變一致。因此，利用銠自給能中子探測(cè)器作中子測(cè)量裝置的先進(jìn)堆芯測(cè)量系 統(tǒng)，為了保證中子通量測(cè)量的準(zhǔn)確性，需要對(duì)銠自給能探器的電流信號(hào)作延遲消除處理。
[0003] 由于實(shí)際的測(cè)量過(guò)程中總伴隨有噪聲（過(guò)程噪聲和測(cè)量噪聲），利用直接的數(shù)學(xué) 反演方法作延遲消除會(huì)將探測(cè)器電流信號(hào)噪聲放大，最大可放大到20倍，影響測(cè)量的精 度。因此，在延遲消除處理過(guò)程中，需要有效抑制噪聲的放大。
[0004] 目前應(yīng)用于銠自給能探測(cè)器信號(hào)延遲的消除主要基于Kalman濾波器實(shí)現(xiàn)，其應(yīng) 用時(shí)必須假定系統(tǒng)的外部擾動(dòng)輸入信號(hào)是一個(gè)具有已知統(tǒng)計(jì)特性的白噪聲信號(hào)，當(dāng)輸入信 號(hào)的統(tǒng)計(jì)特性難以得到時(shí)，該方法就難以應(yīng)用。

【發(fā)明內(nèi)容】

[0005] 本發(fā)明的目的在于克服現(xiàn)有技術(shù)的不足，提供了一種基于H濾波的銠自給能探 測(cè)器信號(hào)延遲消除方法，其應(yīng)用時(shí)能對(duì)銠自給能中子探測(cè)器的電流信號(hào)進(jìn)行延遲消除處 理，并能有效抑制噪聲，使得銠自給能中子探測(cè)器在反應(yīng)堆瞬態(tài)工況時(shí)也能正常使用，且由 于本發(fā)明采用了H濾波器，作延遲消除時(shí)無(wú)需預(yù)先知道外部擾動(dòng)輸入信號(hào)的統(tǒng)計(jì)特性。
[0006] 本發(fā)明解決上述問(wèn)題主要通過(guò)以下技術(shù)方案實(shí)現(xiàn)：基于H濾波的銠自給能探測(cè) 器信號(hào)延遲消除方法，包括以下步驟：
[0007] 步驟1、建立銠與熱中子的核反應(yīng)模型：
[0008] 在反應(yīng)堆瞬態(tài)工況下，通量的變化引起銠自給能中子探測(cè)器電流的變化并不同 步，后者較前者有一定的滯后，描述上述反應(yīng)的具體公式如下：
[0009]
【主權(quán)項(xiàng)】
1.基于H濾波的銠自給能探測(cè)器信號(hào)延遲消除方法，其特征在于：包括以下步驟： 步驟1、建立銠與熱中子的核反應(yīng)模型： 在反應(yīng)堆瞬態(tài)工況下，通量的變化引起銠自給能中子探測(cè)器電流的變化并不同步，后 者較前者有一定的滯后，描述上述反應(yīng)的具體公式如下：
I(t) =cn(t) +Ajirij(t) (3) 其中，分別表示1(l4Rh和1(l4mRh直接引起的電荷量，n(t)表示探測(cè)器處熱中 子通量對(duì)應(yīng)的探測(cè)器平衡狀態(tài)下的探測(cè)器電流，L、\2分別表示1MRh和1(l4niRh的衰變常 數(shù)，c表示探測(cè)器電流的瞬時(shí)響應(yīng)份額，&1、&2分別表示 1(l4Rh和1(l4mRh引起的電流份額，I⑴ 表示銘自給能電流； 步驟2、采用去耦變換建立核反應(yīng)模型對(duì)應(yīng)的離散狀態(tài)方程： 對(duì)式（1)、式（2)及式（3)作拉普拉斯變換，得到如下等式：
平衡態(tài)時(shí)，等式變?yōu)?br>于是式（4)變?yōu)?br>對(duì)式（6)進(jìn)行拉普拉斯逆變換，得到如下?tīng)顟B(tài)方程
I(t) = [c,c,c] ?X(t) (9) 其中
初始值
式（7)、式（8)、式（9)對(duì)應(yīng)的離散狀態(tài)方程為
I(k) =[ccc] *X(k) + [l] ?V(k) (12) n(k) = [1 0 0] ?X(k) (13) 其中，
初始值為
步驟3、確定銠自給能探測(cè)器電流的瞬時(shí)響應(yīng)份額； 步驟4、利用H濾波器對(duì)銠自給能探測(cè)器電流信號(hào)作延遲消除： 對(duì)于一個(gè)離散控制過(guò)程系統(tǒng)，該系統(tǒng)可用一個(gè)狀態(tài)方程來(lái)描述：x(k+1) =Ax(k)+Bw(k) y(k) =Cx(k)+Dw(k) (15) z(k) =Lx(k) 其中，x(k)為第k次采樣點(diǎn)的n維狀態(tài)向量，w(k)包含了系統(tǒng)過(guò)程噪聲以及系統(tǒng)觀測(cè) 白噪聲，y(k)為第k次采樣點(diǎn)的測(cè)量值，z(k)為1維待求向量，L為l*n維矩陣； 假定系統(tǒng)是漸近穩(wěn)定的，則對(duì)給定的常數(shù)y>〇,要求設(shè)計(jì)一個(gè)漸近穩(wěn)定的滿(mǎn)階線性濾 波器
系統(tǒng)存在一個(gè)H濾波器，當(dāng)且僅當(dāng)以下的線性矩陣不等式成立
其中Y、Z為待求解的對(duì)稱(chēng)正定矩陣，而Q、G、F為待求解的一般矩陣； 得到上述矩陣后，H濾波器的相關(guān)矩陣表示如下： Af = -Y、(I-Ylr1，Bf =-Y<F，Cf =G(I-YD-1 (18) 對(duì)于銠自給能探測(cè)器，由其離散狀態(tài)方程可知方程（15)中的對(duì)應(yīng)矩陣為：
C=[ccc] D= [0 1] L= [1 0 0] 通過(guò)求解線性矩陣不等式（17)，可得H濾波器矩陣Af、Bf、Cf，從而可以由如下步驟獲 取消除延遲后任意時(shí)刻的探測(cè)器電流值： 由初始電流測(cè)量值j>(〇)可得
，初始〇時(shí)刻延遲消除后電 流值為
對(duì)于任意k+l(k= 0, 1，...）時(shí)刻，
而k+1時(shí)刻延遲消除后 的電流值為
2.根據(jù)權(quán)利要求1所述的基于H濾波的銠自給能探測(cè)器信號(hào)延遲消除方法，其特征 在于，所述步驟3確定銠自給能探測(cè)器電流的瞬時(shí)響應(yīng)份額的具體步驟如下：在反應(yīng)堆啟 動(dòng)物理實(shí)驗(yàn)階段，通過(guò)升/降反應(yīng)堆功率形成功率臺(tái)階，記錄相應(yīng)的堆外探測(cè)器信號(hào)實(shí)測(cè) 值與銠自給能探測(cè)器信號(hào)實(shí)測(cè)值，其中，堆外探測(cè)器能夠瞬時(shí)響應(yīng)中子通量的變化，相應(yīng)的 實(shí)測(cè)值可認(rèn)為是真實(shí)的中子通量；通過(guò)調(diào)整瞬時(shí)響應(yīng)份額的理論值給定N個(gè)不同的瞬時(shí)響 應(yīng)份額預(yù)測(cè)值，再將堆外探測(cè)器信號(hào)實(shí)測(cè)值代入離散狀態(tài)方程，可以得到N組銠自給能探 測(cè)器信號(hào)理論值，將理論值與銠自給能探測(cè)器信號(hào)實(shí)測(cè)值進(jìn)行比較，取其中符合程度最好 的某組理論值相應(yīng)的瞬時(shí)響應(yīng)份額預(yù)測(cè)值為后續(xù)延遲消除所采用的瞬時(shí)響應(yīng)份額。
3. 根據(jù)權(quán)利要求1或2所述的基于H濾波的銠自給能探測(cè)器信號(hào)延遲消除方法，其 特征在于，在有換擋的情況下，還包括按如下的信號(hào)處理方法對(duì)原始信號(hào)進(jìn)行處理：在換檔 區(qū)域內(nèi)，假設(shè)中子通量保持不變，然后反推中子通量密度產(chǎn)生的電流信號(hào)，再與探測(cè)器實(shí)際 輸出電流相減，得到換擋突變分量；在換檔區(qū)域外，探測(cè)器輸出電流減去換擋突變分量，得 到中子通量密度產(chǎn)生的電流信號(hào)，然后再對(duì)此電流信號(hào)進(jìn)行延遲消除處理。
4. 根據(jù)權(quán)利要求3所述的基于H濾波的銠自給能探測(cè)器信號(hào)延遲消除方法，其特征 在于，所述的換檔區(qū)域設(shè)計(jì)結(jié)構(gòu)如下： 在換檔區(qū)域內(nèi)k2)，假設(shè)中子通量密度不變，則有： n(k+1) =n(k) (19)
可以反推出銠自給能探測(cè)器電流信號(hào)為： I(k+1) =c(n(k+1) +x: (k+1) +x2 (k+1)) (22) 將反推電流（22)當(dāng)成探測(cè)器實(shí)際輸出電流，然后進(jìn)行延遲消除； 在換檔區(qū)域時(shí)間邊界匕處，換檔引起的電流偏置量可以由下式進(jìn)行估算：
其中，J'(t)表示在匕時(shí)刻的探測(cè)器實(shí)際輸出電流； 在換檔區(qū)域外，將探測(cè)器實(shí)際輸出電流加上式（23)表示的換檔引起的電流偏置量，得 到中子通量密度產(chǎn)生的電流信號(hào)，然后再對(duì)此電流信號(hào)進(jìn)行延遲消除。
【專(zhuān)利摘要】本發(fā)明公開(kāi)了基于H∞濾波消除銠自給能探測(cè)器信號(hào)延遲的方法，包括依次進(jìn)行的以下步驟：步驟1、建立銠與熱中子的核反應(yīng)模型；步驟2、采用去耦變換建立核反應(yīng)模型對(duì)應(yīng)的離散狀態(tài)方程；步驟3、確定銠自給能探測(cè)器電流的瞬時(shí)響應(yīng)份額；步驟4、利用H∞濾波器對(duì)銠自給能探測(cè)器電流信號(hào)作延遲消除。本發(fā)明應(yīng)用時(shí)能對(duì)銠自給能中子探測(cè)器的電流信號(hào)進(jìn)行延遲消除處理，并能有效抑制噪聲，使得銠自給能中子探測(cè)器在反應(yīng)堆瞬態(tài)工況時(shí)也能正常使用，且本發(fā)明無(wú)需預(yù)先知道外部擾動(dòng)輸入信號(hào)的統(tǒng)計(jì)特性。
【IPC分類(lèi)】G06F19-00, G21C17-108
【公開(kāi)號(hào)】CN104778985
【申請(qǐng)?zhí)枴緾N201510166372
【發(fā)明人】龔禾林, 陳長(zhǎng), 彭星杰, 趙文博, 劉啟偉, 李向陽(yáng), 李慶, 于穎銳
【申請(qǐng)人】中國(guó)核動(dòng)力研究設(shè)計(jì)院
【公開(kāi)日】2015年7月15日
【申請(qǐng)日】2015年4月9日