量值K〇)可I
，初始0时刻延迟消除后电流值为 i(〇) =r,i(〇)；
[0126] 对于任意k+1化=0, 1，...）时刻，；t/i+1) = ,4，刮/。+公,;W)，而k+1时刻延迟消 除后的电流值为+1)=(', +])。
[0127] 实施例2;
[0128] 本实施例在实施例1的基础上做出了如下进一步限定；在有换挡的情况下，所述 步骤4采用W下方式来进行延迟消除：
[0129] 在换档区域内也《k《k2)，假设中子通量密度不变，则有：
[0130] n(k+1) =n(k) (16)
[0131]
[0133] 可W反推出银自给能探测器电流信号为：
[0134] I化+1) =p(n化+1)+Xi化+1)+X2 化+1)) (19)
[0135] 将反推电流（19)当成探测器实际输出电流，通过步骤4进行延迟消除。
[0136] 在换档区域时间边界k2处，换档引起的电流偏置量可W由下式进行估算：
[0137] D = !似-交似 (20)
[013引其中少化)表示在k2时刻的探测器实际输出电流；
[0139] 在换档区域外，需要对探测器实际输出电流进行偏置补偿W抵消换档所带来的影 响，将探测器实际输出电流加上式（20)表示的换档引起的电流偏置量，得到中子通量密度 产生的电流信号，然后再对此电流信号进行延迟消除。
[0140] W上所述，仅是本发明的较佳实施例，并非对本发明做任何形式上的限制，凡是依 据本发明的技术实质上对W上实施例所作的任何简单修改、等同变化，均落入本发明的保 护范围之内。
【主权项】
1.基于Luenberger形式H2/H C?混合滤波的银探测器信号延迟消除方法，其特征在于： 包括以下步骤： 步骤1、建立银与热中子的核反应模型： 在反应堆瞬态工况下，通量的变化引起银自给能中子探测器电流的变化并不同步，后 者较前者有一定的滞后，描述上述反应的具体公式如下：其中，1(l8Ag(t)表示ltl8Ag的核密度，n°Ag(t)表示 n°Ag的核密度，ltl7Ag表示ltl7Ag的核密 度，ltl9Ag表示ltl9Ag的核密度，（Mt)表示探测器处中子通量，1(17 ?表示ltl7Ag的中子俘获截 面，1(ι9σ表示ltl9Ag的中子俘获截面， 1(18 λ表示ltl8Ag的β衰变常数，11(1 λ表示1Mg的β衰 变常数，ici7K表示ltl7Ag俘获中子后产生电流的概率， ltl9K表示ltl9Ag俘获中子后产生电流的概 率，ltl8K表示 ltl8Ag发生β衰变后产生电流的概率，lltlK表示n°Ag发生β衰变后产生电流的 概率，I⑴表示SPND电流； 步骤2、采用去耦变换获取核反应模型对应的离散状态方程： 使用Laplace变换将动态模型化为：I(t) = P (X1 (t)+X2 (t)+X3 (t)) (6) 其中P为瞬时电流份额，q为1uiAgU)发生β衰变对应的延迟电流份额，r为lcl8Ag(t) 发生β衰变对应的延迟电流份额； 将式（4)、（5)和（6)进行时间离散化处理，并且添加过程噪声项和测量噪声项，可以得 到离散状态方程如下：Ik= [p P p]xk+vk (8) <i>k= [I O 0]x k (9) 其中= χ?_ χ?，wk为过程噪声，v k为测量噪声，T s为采样时间， 初始值为(10); 步骤3、确定银自给能探测器电流的瞬时响应份额： 在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆外 探测器信号实测值与银自给能探测器信号实测值；堆外探测器能够瞬时响应中子通量的变 化，相应的实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N个不 同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可以得到N组 银自给能探测器信号理论值，将理论值与银自给能探测器信号实测值进行比较，取其中符 合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响应 份额； 步骤4、利用Luenberger形式的H2/H c?混合滤波器对银自给能探测器电流信号作延迟 消除： 对于一个离散控制过程系统，该系统可用一个状态方程来描述： X (k+1) = Ax (k) +Bw (k) y (k) = Cx (k) +Dw (k) (11) z (k) = Lx (k) 其中，x(k)为第k次采样点的η维状态向量，w(k)包含了系统过程噪声以及系统观测 白噪声，y(k)为第k次采样点的测量值，z (k)为1维待求向量，L为l*n维矩阵； 针对离散系统（11)，设计如下渐近稳定的满阶线性Luenberger滤波器 4+i -C^-) (12) zk = Lxk 式（12)为H2/H c?混合滤波器，对于一个给定的γ，当且仅当如下的优化问题有解： (13) 其中J满足如下矩阵不等式：(14)(15) 其中 Y = YTe RnXn，W e RnXr，J = jTe Rnix' H2/H ①混合滤波器的增益 K = Yl; 对于银自给能探测器，由其离散状态方程可知方程（11)中的对应矩阵为：通过求解线性矩阵不等式（13)、（14)、（15)，可得H2/H c?混合滤波器矩阵K，从而可以 由如下步骤获取消除延迟后任意时刻的探测器电流值： 由初始电流测量值MO)可初始〇时刻延迟消除后电流值为 Z(O) = CVv(O); 对于任意k+l(k = 0, 1，...）时刻，$认+ 1) = 4对幻+ jgJ?⑷，而k+Ι时刻延迟消除后 的电流值为沿+】）=G Λ·(/( +1)。2. 根据权利要求1所述的基于Luenberger形式H 2/H c?混合滤波的银探测器信号延迟 消除方法，其特征在于，在有换挡的情况下，还包括按如下的信号处理方法对原始信号进行 处理：在换档区域内，假设中子通量保持不变，然后反推中子通量密度产生的电流信号，再 与探测器实际输出电流相减，得到换挡突变分量；在换档区域外，探测器输出电流减去换挡 突变分量，得到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除处理。3. 根据权利要求1所述的基于Luenberger形式H 2/H c?混合滤波的银探测器信号延迟 消除方法，其特征在于，所述的换档区域设计结构如下： 在换档区域内Gc1S k2)，假设中子通量密度不变，则有： n (k+1) = n (k) (16)可以反推出银自给能探测器电流信号为： I (k+1) = p (n (k+1) +X1 (k+1) +x2 (k+1)) (19) 将反推电流（19)当成探测器实际输出电流，通过权利要求1中所述步骤进行延迟消 除； 在换档区域时间边界1^2处，换档引起的电流偏置量可以由下式进行估算： D = /(L)- y(L) (20) 其中At)表示在1^2时刻的探测器实际输出电流；在换档区域外，需要对探测器实际输 出电流进行偏置补偿以抵消换档所带来的影响，将探测器实际输出电流加上式（20)表示 的换档引起的电流偏置量，得到中子通量密度产生的电流信号，然后再对此电流信号进行 延迟消除。
【专利摘要】本发明公开了基于Luenberger形式H2/H∞混合滤波的银探测器信号延迟消除方法，包括依次进行的以下步骤：步骤1、建立银与热中子的核反应模型；步骤2、采用去耦变换建立核反应模型对应的离散状态方程；步骤3、确定银自给能探测器电流的瞬时响应份额；步骤4、利用Luenberger形式的H2/H∞混合滤波器对银自给能探测器电流信号作延迟消除。本发明应用时能对银自给能中子探测器的电流信号进行延迟消除处理，并能有效抑制噪声，使得银自给能中子探测器在反应堆瞬态工况时也能正常使用，且由于该方法仅要求测量误差通道所对应的滤波误差方差有一个上界，从而当输入信号是一个具有有限能量的不确定信号时，银自给能中子探测器也能正常应用。
【IPC分类】G21C17/108, G01T3/00, G06F19/00
【公开号】CN104882177
【申请号】CN201510166117
【发明人】彭星杰, 李庆, 龚禾林, 陈长, 赵文博, 刘启伟, 李向阳, 于颖锐
【申请人】中国核动力研究设计院
【公开日】2015年9月2日
【申请日】2015年4月9日