基于h2滤波的铑自给能探测器信号延迟消除方法
【技术领域】
[0001] 本发明涉及先进堆芯测量系统（核反应堆功率分布在线监测系统）所用的堆内铑 自给能中子探测器信号的处理技术，具体是基于H2滤波的铑自给能探测器信号延迟消除 方法。
【背景技术】
[0002] 用作先进堆芯测量系统堆内探测器的铑自给能中子探测器，其敏感材料铑与中子 反应产生的次生核素发生0衰变产生电流，稳态情况下该电流大小与所在位置通量成正 比，因此通过测量铑自给能探测器能够推知其所在位置中子通量。由于该类探测器电流主 要成分是由次生核素0衰变产生的，在反应堆瞬态情况（中子通量水平变化的情况）下， 该类探测器电流不能实时反映通量水平的变化，而是有一定的延迟，延迟时间参数与次生 核素的0衰变一致。因此，利用铑自给能中子探测器作中子测量装置的先进堆芯测量系 统，为了保证中子通量测量的准确性，需要对铑自给能探器的电流信号作延迟消除处理。
[0003] 由于实际的测量过程中总伴随有噪声（过程噪声和测量噪声），利用直接的数学 反演方法作延迟消除会将探测器电流信号噪声放大，最大可放大到20倍，影响测量的精 度。因此，在延迟消除处理过程中，需要有效抑制噪声的放大。
[0004] 目前应用于铑自给能探测器信号延迟的消除主要基于Kalman滤波器实现，其应 用时必须假定系统的外部扰动输入信号是一个具有已知统计特性的白噪声信号，当输入信 号的统计特性难以得到时，该方法就难以应用。
【发明内容】
[0005] 本发明的目的在于克服现有技术的不足，提供了一种基于H2滤波的铑自给能探 测器信号延迟消除方法，其应用时能对铑自给能中子探测器的电流信号进行延迟消除处 理，并能有效抑制噪声，使得铑自给能中子探测器在反应堆瞬态工况时也能正常使用，且由 于本发明采用了H2滤波器，作延迟消除时无需预先知道外部扰动输入信号的统计特性。
[0006] 本发明解决上述问题主要通过以下技术方案实现：基于H2滤波的铑自给能探测 器信号延迟消除方法，包括以下步骤：
[0007] 步骤1、建立铑与热中子的核反应模型：
[0008] 在反应堆瞬态工况下，通量的变化引起铑自给能中子探测器电流的变化并不同 步，后者较前者有一定的滞后，描述上述反应的具体公式如下：
【主权项】
1.基于H2滤波的铑自给能探测器信号延迟消除方法，其特征在于：包括以下步骤： 步骤1、建立铑与热中子的核反应模型： 在反应堆瞬态工况下，通量的变化引起铑自给能中子探测器电流的变化并不同步，后 者较前者有一定的滞后，描述上述反应的具体公式如下：
I(t) =cn(t) +Ajirij(t) (3) 其中，分别表示1(l4Rh和1(l4mRh直接引起的电荷量，n(t)表示探测器处热中 子通量对应的探测器平衡状态下的探测器电流，L、\2分别表示1MRh和1(l4niRh的衰变常 数，c表示探测器电流的瞬时响应份额，&1、&2分别表示 1(l4Rh和1(l4mRh引起的电流份额，I⑴ 表示铭自给能电流； 步骤2、采用直接变换建立核反应模型对应的离散状态方程： 令上⑴=XA(t)代入式（1)、式（2)及式（3)中，将连续时间的常微分方程直接进 行离散化，并添加噪声项，得到以下离散状态方程：
I(k) = [1 0c] ?X(k) + [1] ?V(k) (5) n(k) = [0 0 1] ?X(k) (6) 其中，
，W(k)为过程噪声项，V(k)为测量噪声项， 初始值为
步骤3、确定铑自给能探测器电流的瞬时响应份额： 步骤4、利用H2滤波器对铑自给能探测器电流信号作延迟消除： 对于一个离散控制过程系统，该系统可用一个状态方程来描述：x(k+1) =Ax(k)+Bw(k) y(k) =Cx(k)+Dw(k) (8) z(k) =Lx(k) 其中，x(k)为第k次采样点的n维状态向量，w(k)包含了系统过程噪声以及系统观测 白噪声，y(k)为第k次采样点的测量值，z(k)为1维待求向量，L为l*n维矩阵； 假定系统是渐近稳定的，则对给定的常数traceX)，要求设计一个渐近稳定的满阶线性 滤波器
系统存在一个H2滤波器，当且仅当以下的线性矩阵不等式成立
其中Y、Z为待求解的对称正定矩阵，而Q、G、F为待求解的一般矩阵； 得到上述矩阵后，H2滤波器的相关矩阵表示如下： Af = -Y、(I-Ylr1，Bf =-Y<F，Cf =G(I-YD-1 (12) 对于铑自给能探测器，由其离散状态方程可知方程（8)中的对应矩阵为：
C= [1 0c] D= [0 1] L= [0 0 1] 通过求解线性矩阵不等式（10)、（11)，可得H2滤波器矩阵Af、Bf、Cf，从而可以由如下步 骤获取消除延迟后任意时刻的探测器电流值： 由初始电流测量值j)(〇)可得1
初始〇时刻延迟消除后电流值为
对于任意k+1(k= 0, 1，...）时刻，
，而k+1时刻延迟消除后 的电流值为
2. 根据权利要求1所述的基于H2滤波的铑自给能探测器信号延迟消除方法，其特征在 于，所述步骤三确定铑自给能探测器电流的瞬时响应份额具体包括以下步骤：在反应堆启 动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆外探测器信号实测 值与铑自给能探测器信号实测值，其中，堆外探测器能够瞬时响应中子通量的变化，相应的 实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N个不同的瞬时响 应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可以得到N组铑自给能探 测器信号理论值，将理论值与铑自给能探测器信号实测值进行比较，取其中符合程度最好 的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响应份额。
3. 根据权利要求1或2所述的基于H2滤波的铑自给能探测器信号延迟消除方法，其特 征在于，在有换挡的情况下，还包括按如下的信号处理方法对原始信号进行处理：在换档区 域内，假设中子通量保持不变，然后反推中子通量密度产生的电流信号，再与探测器实际输 出电流相减，得到换挡突变分量；在换档区域外，探测器输出电流减去换挡突变分量，得到 中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除处理。
4. 根据权利要求3所述的基于H2滤波的铑自给能探测器信号延迟消除方法，其特征在 于，所述的换档区域设计结构如下： 在换档区域内k2)，假设中子通量密度不变，则有： n(k+1) =n(k) (13)
可以反推出铑自给能探测器电流信号为： I(k+1) =Ja(k+l)+cn(k+l) (16) 将反推电流（16)当成探测器实际输出电流，然后进行延迟消除； 在换档区域时间边界匕处，换档引起的电流偏置量可以由下式进行估算：
其中，夕氏）表示在匕时刻的探测器实际输出电流； 在换档区域外，将探测器实际输出电流加上式（17)表示的换档引起的电流偏置量，得 到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除。
【专利摘要】本发明公开了基于H2滤波消除铑自给能探测器信号延迟的方法，包括依次进行的以下步骤：步骤1、建立铑与热中子的核反应模型；步骤2、采用直接变换建立核反应模型对应的离散状态方程；步骤3、确定铑自给能探测器电流的瞬时响应份额；步骤4、利用H2滤波器对铑自给能探测器电流信号作延迟消除。本发明应用时能对铑自给能中子探测器的电流信号进行延迟消除处理，并能有效抑制噪声，使得铑自给能中子探测器在反应堆瞬态工况时也能正常使用，且本发明无需预先知道外部扰动输入信号的统计特性。
【IPC分类】G01T3-00, G21C17-108, G06F19-00
【公开号】CN104778984
【申请号】CN201510166196
【发明人】龚禾林, 陈长, 彭星杰, 赵文博, 刘启伟, 李向阳, 李庆, 于颖锐
【申请人】中国核动力研究设计院
【公开日】2015年7月15日
【申请日】2015年4月9日