滤波的银自给能探测器信号延迟消除方法
【技术领域】
[0001] 本发明设及核反应堆功率分布在线监测系统所用的堆内银自给能中子探测器信 号的处理技术，具体是基于滤波的银自给能探测器信号延迟消除方法。
【背景技术】
[0002] 用作先进堆巧测量系统堆内探测器的银自给能中子探测器，其敏感材料银与中子 反应产生的次生核素发生0衰变产生电流，稳态情况下该电流大小与所在位置通量成正 比，因此通过测量银自给能探测器能够推知其所在位置中子通量。由于该类探测器电流主 要成分是由次生核素0衰变产生的，在反应堆瞬态情况（中子通量水平变化的情况）下， 该类探测器电流不能实时反映通量水平的变化，而是有一定的延迟，延迟时间参数与次生 核素的0衰变一致。因此，利用银自给能中子探测器作中子测量装置的先进堆巧测量系 统，为了保证中子通量测量的准确性，需要对银自给能探器的电流信号作延迟消除处理。
[0003] 由于实际的测量过程中总伴随有噪声（过程噪声和测量噪声），利用直接的数学 反演方法作延迟消除会将探测器电流信号噪声放大，最大可放大到20倍，影响测量的精 度。因此，在延迟消除处理过程中，需要有效抑制噪声的放大。
[0004] 目前应用于银自给能探测器信号延迟的消除主要基于Kalman滤波器实现，其应 用时必须假定系统的外部扰动输入信号是一个具有已知统计特性的白噪声信号，当输入信 号的统计特性难W得到时，该方法就难W应用。
【发明内容】
[0005] 本发明的目的在于克服现有技术的不足，提供了一种基于&滤波的银自给能探测 器信号延迟消除方法，其应用时能对银自给能中子探测器的电流信号进行延迟消除处理， 并能有效抑制噪声，使得银自给能中子探测器在反应堆瞬态工况时也能正常使用，且由于 &滤波器不要求了解输入信号的统计特性，且由于本发明采用了H2滤波器，作延迟消除时 无需预先知道外部扰动输入信号的统计特性。
[0006] 本发明解决上述问题主要通过W下技术方案实现；基于&滤波的银自给能探测器 信号延迟消除方法，其特征在于，包括W下步骤：
[0007] 步骤1、建立银与热中子的核反应模型：
[0008] 在反应堆瞬态工况下，通量的变化引起银自给能中子探测器电流的变化并不同 步，后者较前者有一定的滞后，描述上述反应的具体公式如下：
[0009]
[0011]I(t) = (1。化1。9O iwAg+iwRiwoiwAg)Mt)
[001引 做
[0013] +1…护°入1 …Ag(t) + …8k…8AwsAg(t)
[0014]其中，wsAg(t)表示wsAg的核密度，ii°Ag(t)表示ii°Ag的核密度，i^Ag表示i^Ag的 核密度，i^Ag表示i^Ag的核密度，d) (t)表示探测器处中子通量，iw〇表示i^Ag的中子俘获 截面，W90表示i^Ag的中子俘获截面，W8^表示wsAg的P衰变常数，^表示nc>Ag的0 衰变常数，i^K表示i^Ag俘获中子后产生电流的概率，i°9K表示i^Ag俘获中子后产生电流的 概率，WSK表示i^Ag发生0衰变后产生电流的概率，11化表示ii°Ag发生0衰变后产生电流 的概率，I(t)表示SPND电流；本步骤的目的在于由第一性原理出发，导出银自给能探测器 产生信号该一物理过程所对应的连续时间变量数学模型，而核反应模型是应用滤波器进行 延迟消除的基础。
[0015] 步骤2、采用去禪变换获取核反应模型对应的离散状态方程：
[0016] 使用Laplace变换将动态模型化为；
[0019]I(t) =P(xi(t) +X2(t) +X3(t)) (6)
[0020] 其中p为瞬时电流份额，q为ii°Ag(t)发生e衰变对应的延迟电流份额，r为 i°8Ag(t)发生0衰变对应的延迟电流份额；
[002。 将式（4)、妨和做进行时间离散化处理，并且添加过程噪声项和测量噪声项，可 W得到离散状态方程如下：
[0022]
(7)
[0023] Ik=[PPp]Xk+Vk巧）
[0024]<!),=[1 0 0]x, (9)
[00幼其中4了，Wk为过程噪声，Vk为测量噪声，Ts为采样时间；
[002引初始值为[0027] (10)；
[0028] 步骤3、确定银自给能探测器电流的瞬时响应份额：
[0029] 在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的 堆外探测器信号实测值与银自给能探测器信号实测值，堆外探测器能够瞬时响应中子通量 的变化，相应的实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N 个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可W得到N 组银自给能探测器信号理论值，将理论值与银自给能探测器信号实测值进行比较，取其中 符合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响 应份额。
[0030] 步骤4、利用肥滤波器对银自给能探测器电流信号作延迟消除：
[0031] 对于一个离散控制过程系统，该系统可用一个状态方程来描述：
[0032] X化+1) = Ax化)+Bw化）
[0033] y似=Cx似+Dw化）（11)
[0034] Z化）=Lx化）
[003引其中，x(k)为第k次采样点的n维状态向量，W似包含了系统过程噪声化及系统 观测白噪声，y化）为第k次采样点的测量值，Z化）为1维待求向量，L为l*n维矩阵；
[0036] 假定系统是渐近稳定的，则对给定的常数trace〉0,要求设计一个渐近稳定的满阶 线性滤波器
[0037]
(12)
[0038] 系统存在一个肥滤波器，当且仅当W下的线性矩阵不等式成立
[0041] 其中Y、Z为待求解的对称正定矩阵，而Q、G、F为待求解的一般矩阵；
[0042] 得到上述矩阵后，肥滤波器的相关矩阵表示如下：
[0043] Af=-Y-iQ(I-Y-1幻-1，Bf=-Y-中，Cf=G(I-Y-1幻-1 (14)
[0044] 对于银自给能探测器，由其离散状态方程可知方程巧）中的对应矩阵为：
[0050] 通过求解线性矩阵不等式（11)、（12)，可得肥滤波器矩阵Af、Bf、Cf，从而可W由如 下步骤获取消除延迟后任意时刻的探测器电流值：
[0051] 由初始电流测量值j)(0)可福
初始0时刻延迟消除后电流值为 i(o) = c.VW));
[005引对于任意k+1化=0, 1，...）时刻，雖' +1) =A-WO+Av(/0，而k+1时刻延迟消除后的电流值为+1) = ^ +1)。
[0053] 本发明应用时利用&滤波器原理，在延迟消除过程中，可W有效地抑制噪声的放 大，噪声抑制效果越好，延迟效果会逐渐变差，因此，本发明应用时需适当调节参数使延迟 消除效果与噪声抑制达到最佳平衡。
[0054] 当需要探测较大动态范围的中子通量密度，相应的也需要检测大动态范围的电流 信号，而该一问题便集中在了模数转换器上。为了适应大动态范围的电流的量化，银自给能 探测器的模数转换器采样分档电阻，当电流信号在大范围变化时，模数转换器就会发生电 阻档位转换。由于各档位没有完全匹配，各档位之间的切换会造成输出信号的近似于阶跃 的突变。
[00巧]换挡引起的突变分量进入延迟消除模块后，会被严重放大，使得时域上的阶跃突 变被严重放大，影响最终信号延迟消除的质量（突变部分信号的严重失真）。在换挡时间段 里，信号的变化主要由换挡突变贡献，相对而言，由中子通量密度变化引起的电流信号变化 可W忽略。
[0056] 在有换挡的情况下，还包括按如下的信号处理方法对原始信号进行处理；在换档 区域内，假设中子通量保持不变，然后反推中子通量密度产生的电流信号，再与探测器实际 输出电流相减，得到换挡突变分量；在换档区域外，探测器输出电流减去换挡突变分量，得 到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除处理。