基于h∞濾波的銀自給能探測器信號延遲消除方法
【技術領域】
[0001] 本發明涉及核反應堆功率分布在線監測系統所用的堆內銀自給能中子探測器信 號的處理技術，具體是基于H濾波的銀自給能探測器信號延遲消除方法。
【背景技術】
[0002] 用作先進堆芯測量系統堆內探測器的銀自給能中子探測器，其敏感材料銀與中子 反應產生的次生核素發生0衰變產生電流，穩態情況下該電流大小與所在位置通量成正 比，因此通過測量銀自給能探測器能夠推知其所在位置中子通量。由于該類探測器電流主 要成分是由次生核素0衰變產生的，在反應堆瞬態情況（中子通量水平變化的情況）下， 該類探測器電流不能實時反映通量水平的變化，而是有一定的延遲，延遲時間參數與次生 核素的0衰變一致。因此，利用銀自給能中子探測器作中子測量裝置的先進堆芯測量系 統，為了保證中子通量測量的準確性，需要對銀自給能探器的電流信號作延遲消除處理。
[0003] 由于實際的測量過程中總伴隨有噪聲（過程噪聲和測量噪聲），利用直接的數學 反演方法作延遲消除會將探測器電流信號噪聲放大，最大可放大到20倍，影響測量的精 度。因此，在延遲消除處理過程中，需要有效抑制噪聲的放大。
[0004] 目前應用于銀自給能探測器信號延遲的消除主要基于Kalman濾波器實現，其應 用時必須假定系統的外部擾動輸入信號是一個具有已知統計特性的白噪聲信號，當輸入信 號的統計特性難以得到時，該方法就難以應用。
【發明內容】
[0005] 本發明的目的在于克服現有技術的不足，提供了一種基于H濾波的銀自給能探 測器信號延遲消除方法，其應用時能對銀自給能中子探測器的電流信號進行延遲消除處 理，并能有效抑制噪聲，使得銀自給能中子探測器在反應堆瞬態工況時也能正常使用，且由 于H--濾波器不要求了解輸入信號的統計特性，本發明無需預先知道外部擾動輸入信號的 統計特性。
[0006] 本發明解決上述問題主要通過以下技術方案實現：基于H濾波的銀自給能探測 器信號延遲消除方法，其特征在于，包括以下步驟：
[0007] 步驟1、建立銀與熱中子的核反應模型：
[0008] 在反應堆瞬態工況下，通量的變化引起銀自給能中子探測器電流的變化并不同 步，后者較前者有一定的滯后，描述上述反應的具體公式如下：
[0010]
[0011]I(t) = (109K109 0 109+Ag+107K107 〇 107Ag) (i) (t)
[0012] (3)
[0013] +110K110 入110Ag(t) +108K108 入108Ag(t)
[0014]其中，1(l8Ag(t)表示1(l8Ag的核密度，n°Ag(t)表示n°Ag的核密度，1OTAg表示1OTAg的 核密度，1(l9Ag表示1(l9Ag的核密度，巾（t)表示探測器處中子通量，1OT〇表示1(l7Ag的中子俘獲 截面，1°9〇表示1°1 8的中子俘獲截面，1°!^表示1、8的0衰變常數， 11(^表示11°48的0 衰變常數，1OTK表示1OTAg俘獲中子后產生電流的概率，1(I9K表示1(l9Ag俘獲中子后產生電流的 概率，1(I8K表示1(l8Ag發生0衰變后產生電流的概率，n°K表示n°Ag發生0衰變后產生電流 的概率，I(t)表示SPND電流；
[0015] 步驟2、采用去耦變換獲取核反應模型對應的離散狀態方程：
[0016] 使用Laplace變換將動態模型化為：
[0019]I(t) =p(x: (t) +x2 (t) +x3 (t)) (6)
[0020] 其中p為瞬時電流份額，q為1^AgU)發生0衰變對應的延遲電流份額，r為 1CI8Ag(t)發生0衰變對應的延遲電流份額；
[0021] 將式⑷、（5)和（6)進行時間離散化處理，并且添加過程噪聲項和測量噪聲項，可 以得到離散狀態方程如下：
[0023]Ik=[ppp]xk+vk (8)
[0024] <i>k= [1 0 0]xk (9)
[0025] 其中A=|^- < ，wk為過程噪聲，vk為測量噪聲，Ts為采樣時間。
[0026] 初始值為
[0027]
[0028] 步驟3、確定銀自給能探測器電流的瞬時響應份額：
[0029] 在反應堆啟動物理實驗階段，通過升/降反應堆功率形成功率臺階，記錄相應的 堆外探測器信號實測值與銀自給能探測器信號實測值，堆外探測器能夠瞬時響應中子通量 的變化，相應的實測值可認為是真實的中子通量，通過調整瞬時響應份額的理論值給定N 個不同的瞬時響應份額預測值，再將堆外探測器信號實測值代入離散狀態方程，可以得到N 組銀自給能探測器信號理論值，將理論值與銀自給能探測器信號實測值進行比較，取其中 符合程度最好的某組理論值相應的瞬時響應份額預測值為后續延遲消除所采用的瞬時響 應份額；
[0030] 步驟4、利用H濾波器對銀自給能探測器電流信號作延遲消除：
[0031] 對于一個離散控制過程系統，該系統可用一個狀態方程來描述：
[0032] x(k+l) =Ax(k)+Bw(k)
[0033] y(k) =Cx(k)+Dw(k) (11)
[0034] z(k) =Lx(k)
[0035] 其中，x(k)為第k次采樣點的n維狀態向量，w(k)包含了系統過程噪聲以及系統 觀測白噪聲，y(k)為第k次采樣點的測量值，z(k)為1維待求向量，L為l*n維矩陣；
[0036] 假定系統是漸近穩定的，則對給定的常數丫>〇,要求設計一個漸近穩定的滿階線 性濾波器
[0038] 系統存在一個H濾波器，當且僅當以下的線性矩陣不等式成立
[0040] 其中Y、Z為待求解的對稱正定矩陣，而Q、G、F為待求解的一般矩陣；
[0041] 得到上述矩陣后，H濾波器的相關矩陣表示如下：
[0042] Af =-Y-1Q(I-Ylr1，Bf =-Y-1F，Cf =G(I-YD-1 (14)
[0043] 對于銀自給能探測器，由其離散狀態方程可知方程（11)中的對應矩陣為：
[0046] C=[ppp]
[0047] D= [0 1]
[0048] L= [1 0 0]
[0049] 通過求解線性矩陣不等式（13)，可得濾波器矩陣Af、Bf、Cf，從而可以由如下步 驟獲取消除延遲后任意時刻的探測器電流值：
[0050] 由初始電流測量值j>(〇)可得
，初始〇時刻延遲消除后電流值為 z(0) =Cfx(Q);
[0051] 對于任意k+1(k= 0, 1，? ? ?）時刻，
，而k+1時刻延遲消除 后的電流值為
[0052] 本發明應用時利用H濾波器原理，在延遲消除過程中，可以有效地抑制噪聲的 放大，噪聲抑制效果越好，延遲效果會逐漸變差，因此，本發明應用時需適當調節參數使延 遲消除效果與噪聲抑制達到最佳平衡。
[0053] 當需要探測較大動態范圍的中子通量密度，相應的也需要檢測大動態范圍的電流 信號，而這一問題便集中在了模數轉換器上。為了適應大動態范圍的電流的量化，銠自給能 探測器的模數轉換器采樣分檔電阻，當電流信號在大范圍變化時，模數轉換器就會發生電 阻檔位轉換。由于各檔位沒有完全匹配，各檔位之間的切換會造成輸出信號的近似于階躍 的突變。
[0054] 換擋引起的突變分量進入延遲消除模塊后，會被嚴重放大，使得時域上的階躍突 變被嚴重放大，影響最終信號延遲消除的質量（突變部分信號的嚴重失真）。在換擋時間段 里，信號的變化主要由換擋突變貢獻，相對而言，由中子通量密度變化引起的電流信號變化 可以忽略。
[0055] 在有換擋的情況下，還包括按如下的信號處理方法對原始信號進行處理：在換檔 區域內，假設中子通量保持不變，然后反推中子通量密度產生的電流信號，再與探測器實際 輸出電流相減，得到換擋突變分量；在換檔區域外，探測器輸出電流減去換擋突變分量，得 到中子通量密度產生的電流信號，然后再對此電流信號進行延遲消除處理。
[0056] 3、根據權