基于Luenberger形式H∞滤波的铑自给能探测器信号延迟消除方法
【技术领域】
[0001] 本发明设及核反应堆功率分布在线监测系统所用的堆内锭自给能中子探测器信 号的处理技术，具体是基于Luenberger形式H-滤波的锭自给能探测器信号延迟消除方 法。
【背景技术】
[0002] 用作先进堆巧测量系统堆内探测器的锭自给能中子探测器，其敏感材料锭与中子 反应产生的次生核素发生0衰变产生电流，稳态情况下该电流大小与所在位置通量成正 比，因此通过测量锭自给能探测器能够推知其所在位置中子通量。由于该类探测器电流主 要成分是由次生核素0衰变产生的，在反应堆瞬态情况（中子通量水平变化的情况）下， 该类探测器电流不能实时反映通量水平的变化，而是有一定的延迟，延迟时间参数与次生 核素的0衰变一致。因此，利用锭自给能中子探测器作中子测量装置的先进堆巧测量系 统，为了保证中子通量测量的准确性，需要对锭自给能探器的电流信号作延迟消除处理。
[0003] 由于实际的测量过程中总伴随有噪声（过程噪声和测量噪声），利用直接的数学 反演方法作延迟消除会将探测器电流信号噪声放大，最大可放大到20倍，影响测量的精 度。因此，在延迟消除处理过程中，需要有效抑制噪声的放大。
[0004] 目前应用于锭自给能探测器信号延迟的消除主要基于Kalman滤波器实现，其应 用时必须假定系统的外部扰动输入信号是一个具有已知统计特性的白噪声信号，当输入信 号是一个具有有限能量的不确定信号，其统计特性难W得到，该方法就难W应用。
【发明内容】
[0005] 本发明的目的在于克服现有技术的不足，提供了一种基于Luenberger形式H-滤 波的锭自给能探测器信号延迟消除方法，其应用时能对锭自给能中子探测器的电流信号进 行延迟消除处理，并能有效抑制噪声，使得锭自给能中子探测器在反应堆瞬态工况时也能 正常使用，且由于本发明采用了Luenberger形式的H-滤波器，作延迟消除时无需预先知 道外部扰动输入信号的统计特性。
[0006] 本发明解决上述问题主要通过W下技术方案实现；基于Luenberger形式滤波 的锭自给能探测器信号延迟消除方法，其特征在于，包括W下步骤：
[0007] 步骤1、建立锭与热中子的核反应模型：
[0008] 在反应堆瞬态工况下，通量的变化引起锭自给能中子探测器电流的变化并不同 步，后者较前者有一定的滞后，描述上述反应的具体公式如下：
[0009] =cun(f) - /U1U(/) ( 1 ) dt一一-
[0010] (2)
[0011]I(t) =cn(t) +Ajiiii(t) (3)
[0012] 其中，mi(t)、m2(t)分别表示IM化和iMm化直接引起的电荷量，n(t)表示探测器处 热中子通量对应的探测器平衡状态下的探测器电流，A1、A2分别表示化和IMm化的衰 变常数，C表示探测器电流的瞬时响应份额，ai、a2分别表示化和IMm化引起的电流份额， I(t)表示锭自给能电流；本步骤的目的在于由第一性原理出发，导出锭自给能探测器产生 信号该一物理过程所对应的连续时间变量数学模型，而核反应模型是应用滤波器进行延迟 消除的基础；
[0013] 步骤2、采用直接变换建立核反应模型对应的离散状态方程：
[0014] 由于探测器的电流信号都是通过离散采样获取，步骤1建立的连续状态方程需要 转换为离散状态方程；令Ja(t) =Aimi(t)代入式（1)、式（2)及式（3)中，将连续时间的常 微分方程直接进行离散化，并添加噪声项，得到W下离散状态方程：
[0015]
(4)
[0016]I(k) = [1 0c] ?X(k) + [1] ?V(k) (5)
[0017]n(k) = [0 0 1] .X〇〇 化）
[0018] 其中，
IW(k)为过程噪声项，V(k)为测量噪声项，
[001引初始值为 [0020]
(7)，
[0021] 步骤3、确定锭自给能探测器电流的瞬时响应份额：
[0022] 锭自给能探测器瞬时响应份额C可W由理论进行大致地估算，但理论估算值与真 实值之间的不匹配将会导致滤波效果的下降，可W设计如下实验较准确地确定出瞬时响应 份额C，在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆 外探测器信号实测值与锭自给能探测器信号实测值。堆外探测器能够瞬时响应中子通量 的变化，相应的实测值可认为是真实的中子通量。通过调整瞬时响应份额的理论值给定N 个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可W得到N 组锭自给能探测器信号理论值，将理论值与锭自给能探测器信号实测值进行比较，取其中 符合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响 应份额；
[002引步骤4、利用Luenberger形式的H-滤波器对锭自给能探测器电流信号作延迟消 除：
[0024] 对于一个离散控制过程系统，该系统可用一个状态方程来描述：
[00巧]X化+1)= Ax(k)+Bw化）
[0026]y似=Cx似+Dw化） 做
[0027] Z化）=Lx化）
[002引其中，x(k)为第k次采样点的n维状态向量，W似包含了系统过程噪声化及系统 观测白噪声，y化）为第k次采样点的测量值，Z化）为1维待求向量，L为l*n维矩阵；
[0029] 针对离散系统巧），设计如下渐近稳定的满阶线性Luenberger滤波器
[0030] 私"=少*+巧於-巧i) (9) 乏k二L文k
[0031] 式巧）为滤波器，对于一个给定的丫，当且仅当如下的矩阵不等式有解：
[0032]
(10)
[0033]其中Y=YTgRDXn，WeRDXr，J二jTeRmXm，H滤波器的增益K=Y-lW;
[0034] 对于锭自给能探测器，由其离散状态方程可知方程巧）中的对应矩阵为：
[0040] 通过求解线性矩阵不等式（10)，可得滤波器矩阵K，从而可W由如下步骤获取 消除延迟后任意时刻的探测器电流值：
[0041] 由初始电流测量值K〇)可得
初始0时刻延迟消除后电流 值为i(o) =c,i(〇);
[004引对于任意k+1 化=0, 1，...）时刻，-;(/、+1)=,如(/0 + 巧(/0)，而k+1 时刻 延迟消除后的电流值为+1)=+])。
[0043] 本发明应用时利用Luenberger形式的H滤波器原理，在延迟消除过程中，可W 有效地抑制噪声的放大，噪声抑制效果越好，延迟效果会逐渐变差，因此，本发明应用时需 适当调节参数使延迟消除效果与噪声抑制达到最佳平衡。
[0044] 当需要探测较大动态范围的中子通量密度，相应的也需要检测大动态范围的电流 信号，而该一问题便集中在了模数转换器上。为了适应大动态范围的电流的量化，锭自给能 探测器的模数转换器采样分培电阻，当电流信号在大范围变化时，模数转换器就会发生电 阻档位转换。由于各档位没有完全匹配，各档位之间的切换会造成输出信号的近似于阶跃 的突变。
[0045] 在有换挡的情况下，还包括按如下的信号处理方法对原始信号进行处理；在换档 区域内，假设中子通量保持不变，然后反推中子通量密度产生的电流信号，再与探测器实际 输出电流相减，得到换挡突变分量；在换档区域外，探测器输出电流减去换挡突变分量，得 到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除处理。
[0046] 换挡引起的突变分量进入延迟消除模块后，会被严重放大，使得时域上的阶跃突 变被严重放大，影响最终信号延迟消除的质量（突变部分信号的严重失真）。在换挡时间段 里，信号的变化主要由换挡突变贡献，相对而言，由中子通量密度变化引起的电流信号变化 可W忽略。
[0047] 为了处理换档导致的信号突变问题，进一步的，在有换挡的情况下，所述步骤4采 用W下方式来进行延迟消除：
[004引在换档区域内也《k《k2)，假设中子通量密度不变，则有：
[0049] n(k+l) = n(k) (11)
[0050]
[0052] 可W反推出锭自给能探测器电流信号为：
[0053] I(k+1) =Ja化+l)+cn化+1) (14)
[0054] 将反推电流（14)当成探测器实际输出电流，通过权利要求1中所述步骤进行延迟 消除；
[0055] 在换档区域时间边界k,处，换档引起的电流偏置量可W由下式进行估算：
[0056] 公= /(A'2)-v(A'2) (15)
[0057] 其中夕(^)表示在k2时刻的探测器实际输出电流；在换档区域外，需要对探测器 实际输出电流进行偏置补偿W抵消换档所带来的影响，将探测器实际输出电流加上式（15) 表示的换档引起的电流偏置量，得到中子通量密度产生的电流信号，然后再对此电流信号 进行延迟消除。
[0058] 综上所述，本发明具有W下有益效果：
[0059] 1本发明整体工序简单，便于实现，能对锭自给能中子探测器的电流信号进行延迟 消除处理，并能有效抑制噪声，使得锭自给能中子探测器在反应堆瞬态工况时也能正常使 用；本发明基于Luenberger形式的H-滤波器实现，在输入信号是一个具有有限能量的不 确定信号时也能正常应用；本发明应用时将滤波器设计转化为相应线性矩阵不等式计算， 方便计算，可W方便地使用Matl油的LMIToo化OX进行求解；
[0060] 2本发明解决了先进堆巧测量系统（核反应堆功率分布在线监测系统）所用的堆 内锭自给能中子探测器信号的延迟消除问题。利用滤波器对锭自给能中子探测器信号 进行延迟消除、平滑、降噪处理，通过适当选取Luenberger形式的H-滤波器参数，能够很 好的达到信号延迟消除效果和噪声抑制效果的最佳平衡，本发明能够保证锭自给能探测器 电流信号直接用于先进堆巧测量系统后续环节，而不丧失准确度；
[0061] 3本发明对锭自给能中子探测器的电流信号进行延迟消除处理，响应时间即阶跃 通量变化时，信号恢复到稳态电流的90%所需的时间在2~10秒内；
[0062] 4本发明对锭自给能中子探测器的电流信号延迟消除过程中，对测量电流信号进 行降噪处理，噪声放大倍数即延迟消除处理后的电流相对误差与噪声之比抑制在1~8 倍；
[0063] 5本发明能有效处理因硬件换挡造成的阶跃对延迟消除效果的影响。
【附图说明】
[0064] 图1为本发明的锭自给能中子探测器结构图
[0065] 图2为本发明一个具体实施例的处理流程图；
[0066] 图3为锭与热中子核反应图。
[0067] 附图中标记及相应的零部件名称：
[006引1-发射极，2-绝缘层，3-收集极，4-导线，5-保护壳，6-绝缘电缆，7-电流线， 8- 本底线，9-密封管，10-电流输出端。
【具体实施方式】
[0069] 下面结合实施例及附图，对本发明做进一步地的详细说明，但本发明的实施方式 不限于此。
[0070] 实施例：
[0071] 如图1所示的锭自给能中子探测器结构图，其中各个序号的零部件名称对应为： 1-发射极，2-绝缘层，3-收集极，4-导线，5-保护壳，6-绝缘电缆，7-电流线，8-本底线， 9- 密封管，10-电流输出端，该锭自给能中子探测器，其特性参数为：A1=