基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法

基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法
【技术领域】
[0001] 本发明设及核反应堆功率分布在线监测系统所用的堆内银自给能中子探测器信 号的处理技术，具体是基于Luenberger形式H-滤波的银自给能探测器信号延迟消除方 法。
【背景技术】
[0002] 用作先进堆巧测量系统堆内探测器的银自给能中子探测器，其敏感材料银与中子 反应产生的次生核素发生0衰变产生电流，稳态情况下该电流大小与所在位置通量成正 比，因此通过测量银自给能探测器能够推知其所在位置中子通量。由于该类探测器电流主 要成分是由次生核素0衰变产生的，在反应堆瞬态情况（中子通量水平变化的情况）下， 该类探测器电流不能实时反映通量水平的变化，而是有一定的延迟，延迟时间参数与次生 核素的0衰变一致。因此，利用银自给能中子探测器作中子测量装置的先进堆巧测量系 统，为了保证中子通量测量的准确性，需要对银自给能探器的电流信号作延迟消除处理。
[0003] 由于实际的测量过程中总伴随有噪声（过程噪声和测量噪声），利用直接的数学 反演方法作延迟消除会将探测器电流信号噪声放大，最大可放大到20倍，影响测量的精 度。因此，在延迟消除处理过程中，需要有效抑制噪声的放大。
[0004] 目前应用于银自给能探测器信号延迟的消除主要基于Kalman滤波器实现，其应 用时必须假定系统的外部扰动输入信号是一个具有已知统计特性的白噪声信号，当输入信 号是一个具有有限能量的不确定信号，其统计特性难W得到，该方法就难W应用。

【发明内容】

[0005]本发明的目的在于克服现有技术的不足，提供了一种基于Luenberger形式H-滤 波的银自给能探测器信号延迟消除方法，其应用时能对银自给能中子探测器的电流信号进 行延迟消除处理，并能有效抑制噪声，使得银自给能中子探测器在反应堆瞬态工况时也能 正常使用，且由于本发明采用了Luenberger形式的馬滤波器，作延迟消除时无需预先知道 外部扰动输入信号的统计特性。
[0006]本发明解决上述问题主要通过W下技术方案实现；基于Luenberger形式滤波 的银自给能探测器信号延迟消除方法，其特征在于，包括W下步骤：
[0007] 步骤1、建立银与热中子的核反应模型：
[0008]在反应堆瞬态工况下，通量的变化引起银自给能中子探测器电流的变化并不同 步，后者较前者有一定的滞后，描述上述反应的具体公式如下：
[0011] (3)
[001引其中，wsAg(t)表示wsAg的核密度，ii°Ag(t)表示ii°Ag的核密度，i^Ag表示i^Ag的 核密度，i^Ag表示i^Ag的核密度，4 (t)表示探测器处中子通量，iw0表示i^Ag的中子俘获 截面，W9 0表示l^Ag的中子俘获截面，W8 ^表示WSAg的P衰变常数，1K) ^表示n0Ag的0 衰变常数，i^K表示i^Ag俘获中子后产生电流的概率，i°9K表示i^Ag俘获中子后产生电流的 概率，WSK表示i^Ag发生0衰变后产生电流的概率，11化表示ii°Ag发生0衰变后产生电流 的概率，I(t)表示SPND电流；
[0013] 步骤2、采用去禪变换获取核反应模型对应的离散状态方程：
[0014] 使用Laplace变换将动态模型化为；
[0015]
[0017]I(t) = P(Xi(t)+X2(t)+X3(t))化）
[001引其中P为瞬时电流份额，q为ii°Ag(t)发生0衰变对应的延迟电流份额，r为wsAg(t)发生0衰变对应的延迟电流份额；
[001引将式（4)、妨和做进行时间离散化处理，并且添加过程噪声项和测量噪声项，可W得到离散状态方程如下：
[0020]
[0021] Ik= [P P p]Xk+Vk (8)
[002引 4k= [1 00]Xk (9)
[002引其中Xi* 4了，Wk为过程噪声，V巧测量噪声，Ts为采样时间；
[0024] 初始值为
[00巧] (10);
[0026] 步骤3、确定银自给能探测器电流的瞬时响应份额：
[0027] 在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的 堆外探测器信号实测值与银自给能探测器信号实测值；堆外探测器能够瞬时响应中子通量 的变化，相应的实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N 个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可W得到N 组银自给能探测器信号理论值，将理论值与银自给能探测器信号实测值进行比较，取其中 符合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响 应份额；
[002引步骤4、利用Luenberger形式的H-滤波器对银自给能探测器电流信号作延迟消 除：
[0029] 对于一个离散控制过程系统，该系统可用一个状态方程来描述：
[0030]X化+1)=Ax(k)+Bw化）
[0031] y(k)=CxW+DwW 做
[0032]Z(k)=Lx化）
[003引其中，x(k)为第k次采样点的n维状态向量，W似包含了系统过程噪声化及系统 观测白噪声，y化）为第k次采样点的测量值，Z化）为1维待求向量，L为l*n维矩阵； [0034] 针对离散系统巧），设计如下渐近稳定的满阶线性Luenberger滤波器
[00对
巧）
[003引式（9)为H-滤波器，对于一个给定的丫，当且仅当如下的矩阵不等式有解：
[0037]
(10)
[0038]其中Y=YTgRDXn，W eRnxr，J二jTeRmXm，H滤波器的增益K=Y-lW;
[0039] 对于银自给能探测器，由其离散状态方程可知方程巧）中的对应矩阵为：
[0045] 通过求解线性矩阵不等式（10)，可得H-滤波器矩阵K，从而可W由如下步骤获取 消除延迟后任意时刻的探测器电流值：
[0046] 由初始电流测量值j)(〇)可得
初始0时刻延迟消除后电流值为 之(〇) = (，, i(〇);
[0047] 对于任意1^+1化=0，1，...）时刻，；-(/、+1) = ,'4,种^') +公,.雌'），而1^+1时刻延迟消 除后的电流值为+1) =V'# + 0。
[0048] 本发明应用时利用Luenberger形式的H滤波器原理，在延迟消除过程中，可W 有效地抑制噪声的放大，噪声抑制效果越好，延迟效果会逐渐变差，因此，本发明应用时需 适当调节参数使延迟消除效果与噪声抑制达到最佳平衡。
[0049] 当需要探测较大动态范围的中子通量密度，相应的也需要检测大动态范围的电流 信号，而该一问题便集中在了模数转换器上。为了适应大动态范围的电流的量化，银自给能 探测器的模数转换器采样分档电阻，当电流信号在大范围变化时，模数转换器就会发生电 阻档位转换。由于各档位没有完全匹配，各档位之间的切换会造成输出信号的近似于阶跃 的突变。
[0050] 换挡引起的突变分量进入延迟消除模块后，会被严重放大，使得时域上的阶跃突 变被严重放大，影响最终信号延迟消除的质量（突变部分信号的严重失真）。在换挡时间段 里，信号的变化主要由换挡突变贡献，相对而言，由中子通量密度变化引起的电流信号变化 可W忽略。
[0051] 为了处理换档导致的信号突变问题，在有换挡的情况下，还包括按如下的信号处 理方法对原始信号进行处理；在换档区域内，假设中子通量保持不变，然后反推中子通量密 度产生的电流信号，再与探测器实际输出电流相减，得到换挡突变分量；在换档区域外，探 测器输出电流减去换挡突变分量，得到中子通量密度产生的电流信号，然后再对此电流信 号进行延迟消除处理。
[0052] 所述的换档区域设计结构如下：
[005引在换档区域内也《k《k2)，假设中子通量密度不变，则有：
[0054]n(k+1) =n(k) (11)
[0057] 可W反推出银自给能探测器电流信号为： [0058]I化+1) =p(n化+1)+Xi化+1)+X2 化+1)) (14)
[0059] 将反推电流（14)当成探测器实际输出电流，通过权利要求1中所述步骤进行延迟 消除；
[0060] 在换档区域时间边界k2处，换档引起的电流偏置量可W由下式进行估算：
[0061] D= !似-y(k,) (15)
[006引其中夕诉)表示在k2时刻的探测器实际输出电流；
[0063] 在换档区域外，需要对探测器实际输出电流进行偏置补偿W抵消换档所带来的影 响，将探测器实际输出电流加上式（15)表示的换档引起的电流偏置量，得到中子通量密度 产生的电流信号，然后再对此电流信号进行延迟消除。
[0064] 在换档区域外，需要对探测器实际输出电流进行偏置补偿W抵消换档所带来的影 响，将探测器实际输出电流加上式（15)表示的换档引起的电流偏置量，得到中子通量密度 产生的电流信号，然后再对此电流信号进行延迟消除。
[0065] 综上所述，本发明具有W下有益效果：
[0066] 1本发明整体工序简单，便于实现，能对银自给能中子探测器的电流信号进行延迟 消除处理，并能有效抑制噪声，使得银自给能中子探测器在反应堆瞬态工况时也能正常使 用；本发明基于Luenberger形式的H-滤波器实现，在输入信号是一个具有有限能量的不 确定信号时也能正常应用；本发明应用时将滤波器设计转化为相应线性矩阵不等式计算， 方便计算，可W方便地使用Matl油的LMI Too化OX进行求解。
[0067] 2本发明解决了先进堆巧测量系统（核反应堆功率分布在线监测系统）所用的堆 内银自给能中子探测器信号的延迟消除问题。利用滤波器对银自给能中子探测器信号 进行延迟消除、平滑、降噪处理，通过适当选取Luenberger形式的H-滤波器参数，能够很 好的达到信号延迟消除效果和噪声抑制效果的最佳平衡。本发明能够保证银自给能探测器 电流信号直接用于先进堆巧测量系统后续环节，而不丧失准确度。
[006引 3本发明对银自给能中子探测器的电流信号进行延迟消除处理，响应时间（阶跃 通量变化时，信号恢复到稳态电流的90%所需的时间）在2~5秒内。
[0069] 4本发明对银自给能中子探测器的电流信号延迟消除过程中，对测量电流信号进 行降噪处理，噪声放大倍数（延迟消除处理后的电流相对误差与噪声之比）抑制在1~4 倍。
[0070] 5本发明能有效处理因硬件换挡造成的阶跃对延迟消除效果的影响。
【附图说明】
[0071] 图1为本发明的银自给能中子探测器结构图
[0072] 图2为本发明一个具体实施例的处理流程图；
[0073] 图3为银与热中子核反应图。
[0074] 附图中标记及相应的零部件名称：
[007引 1-发射极，2-绝缘层，3-收集极，4-导线，5-保护壳，6-绝缘电缆，7-电流线， 8- 本底线，9-密封管，10-电流输出端。
【具体实施方式】
[0076] 下面结合实施例及附图，对本发明做进一步地的详细说明，但本发明的实施方式 不限于此。
[0077] 实施例；
[0078] 如图1所示的银自给能中子探测器结构图，其中各个序号的零部件名称对应为： 1-发射极，2-绝缘层，3-收集极，4-导线，5-保护壳，6-绝缘电缆，7-电流线，8-本底线， 9- 密封管，10-电流输出端，该银自给能中子探测器，其特性参数为：A1。8=ln2/2. 42/60S-1 =0. 0048S-1, ^。。=ln2/24. 4s_i= 0. 0284s-i，p= 0. 09，q= 0. 66，r= 0. 25 ;图 3 为银与 中子核反应原理过程图，对于图3的反应过程中，采用图1的装置进行测量。如图2所示， 基于Luenberger形式H-滤波的银自给能探测器信号延迟消除方法，包括依次进行的W下 步骤；步骤1、建立银与热中子的核反应模型；步骤2、采用去禪变换建立核反应模型对应的 离散状态方程；步骤3、确定银自给能探测器电流的瞬时响应份额；步骤4、利用Luenberger 形式的H-滤波器对银自给能探测器电流信号作延迟消除。
[0079] 本实施例建立银与热中子的核反应模型的具体实施步骤如下；如图2所示，在反 应堆瞬态工况下，通量的变化引起银自给能中子探测器电流的变化并不同步，后者较前者 有一定的滞后，描述上述反应的具体公式如下：
(3)
[0083]其中，wsAg(t)表示wsAg的核密度，ii°Ag(t)表示ii°Ag的核密度，i^Ag表示i^Ag的 核密度，i^Ag表示i^Ag的核密度，4 (t)表示探测器处中子通量，表示的中子俘获 截面，W9 0表示l^Ag的中子俘获截面，W8 ^表示WSAg的P衰变常数，1K) ^表示n〇Ag的0 衰变常数，i^K表示i^Ag俘获中子后产生电流的概率，i°9K表示i^Ag俘获中子后产生电流的 概率，WSK表示i^Ag发生0衰变后产生电流的概率，11化表示ii°Ag发生0衰变后产生电流 的概率，I(t)表示SPND电流。
[0084] 本实施例采用去禪变换建立核反应模型对应的离散状态方程的具体实施步骤如 下：
[0085] 使用Laplace变换将动态模型化为；
[0088]I(t) =P(xi(t) +X2(t) +X3(t)) (6)
[0089] 其中p为瞬时电流份额，q为ii°Ag(t)发生e衰变对应的延迟电流份额，r为 wsAg(t)发生0衰变对应的延迟电流份额。
[0090]将式（4)、妨和做进行时间离散化处理，并且添加过程噪声项和测量噪声项，可 W得到离散状态方程如下：
[0092] Ik= [P P p]Xk+Vk 巧）[009引（})k=[1 00]Xk (9)
[0091] (7)
[0094]其中=[病Xi*- 4了，Wk为过程噪声，V巧测量噪声，Ts为采样时间；
[009引初始值为「009引
vlO)
[0097] 本实施例确定银自给能探测器电流的瞬时响应份额的具体实施步骤如下：
[0098] 在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的 堆外探测器信号实测值与银自给能探测器信号实测值。堆外探测器能够瞬时响应中子通 量的变化，相应的实测值可认为是真实的中子通量。通过调整瞬时响应份额的理论值给定 N个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可W得到 N组银自给能探测器信号理论值，将理论值与银自给能探测器信号实测值进行比较，取其中 符合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响 应份额。
[0099] 本实施例利用H-滤波器对银自给能探测器电流信号作延迟消除的具体实施步 骤如下：
[0100] 对于一个离散控制过程系统，该系统可用一个状态方程来描述：
[0101] X化+1) =Ax(k)+Bw化）
[010引 y(k)=CxW+DwW (11)
[010引 Z化）=Lx化）
[0104] 其中，x(k)为第k次采样点的n维状态向量，w(k)包含了系统过程噪声W及系统 观测白噪声，y化）为第k次采样点的测量值，Z化）为1维待求向量，L为l*n维矩阵；
[0105]针对离散系统（11)，设计如下渐近稳定的满阶线性Luenberger滤波器
[0106] 壬W=作+化巧-) (12) =Lx, 一
[0107] 式（12)为滤波器，对于一个给定的丫，当且仅当如下的矩阵不等式有解： [010 引
(13)
[0109] 其中Y=YTgRDXn，WeRnxr，J二jTeRmXm，H滤波器的增益K=Y-lW;
[0110] 对于银自给能探测器，由其离散状态方程可知方程（11)中的对应矩阵为：
[0113] C= [P P P]
[0114] D= [0 1]
[0115] L= [1 0 0]
[0116] 通过求解线性矩阵不等式（13 )，可得滤波器矩阵K，从而可W由如下步骤获取 消除延迟后任意时刻的探测器电流值：
[0117] 由初始电流测量值K〇)可得
初始0时刻延迟消除后电流值为 身(())=(.、,刮〇);
[0118] 对于任意k+1化=0,1，...）时刻，种/、'+U= +公幻，而k+1时刻延迟消 除后的电流值为+ 0 =C, + 0。
[011引 实施例2;
[0120] 本实施例在实施例1的基础上做出了如下进一步限定；在有换挡的情况下，所述 步骤4采用W下方式来进行延迟消除：
[01引]在换档区域内也《k《k2)，假设中子通量密度不变，则有：
[0122]n(k+1) =n(k) (14)
[0125] 可W反推出银自给能探测器电流信号为：
[0126]I化+1) =P (n化+1)+Xi化+1)+X2化+1)) (17)
[0127] 将反推电流（17)当成探测器实际输出电流，通过步骤4进行延迟消除。
[012引在换档区域时间边界k,处，换档引起的电流偏置量可W由下式进行估算：
[0129] D = I{k^-S又k.)(18)
[0130] 其中兴&2)表示在k2时刻的探测器实际输出电流。
[0131] 在换档区域外，需要对探测器实际输出电流进行偏置补偿W抵消换档所带来的影 响，将探测器实际输出电流加上式（18)表示的换档引起的电流偏置量，得到中子通量密度 产生的电流信号，然后再对此电流信号进行延迟消除。
[0132]W上所述，仅是本发明的较佳实施例，并非对本发明做任何形式上的限制，凡是依 据本发明的技术实质上对W上实施例所作的任何简单修改、等同变化，均落入本发明的保 护范围之内。
【主权项】
1.基于Luenberger形式H C?滤波的银自给能探测器信号延迟消除方法，其特征在于： 包括以下步骤： 步骤1、建立银与热中子的核反应模型： 在反应堆瞬态工况下，通量的变化引起银自给能中子探测器电流的变化并不同步，后 者较前者有一定的滞后，描述上述反应的具体公式如下：+iioKiioAnoAga)+10? 108 λ i〇8Ag(t) 其中，ici8Ag (t)表示ici8Ag的核密度，n°Ag(t)表示n°Ag的核密度， ici7Ag表示ici7Ag的核密 度，ltl9Ag表示ltl9Ag的核密度，（Mt)表示探测器处中子通量，1(%表示 ltl7Ag的中子俘获截 面，1(ι9σ表示ltl9Ag的中子俘获截面， 1(18 λ表示ltl8Ag的β衰变常数，11(1 λ表示n°Ag的β衰 变常数，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)进行时间离散化处理，并且添加过程噪声项和测量噪声项，可以得 到离散状态方程如下：其中Xi =[戎xf ^，Wk为过程噪声，V k为测量噪声，T s为采样时间； 初始值为步骤3、确定银自给能探测器电流的瞬时响应份额： 在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆外 探测器信号实测值与银自给能探测器信号实测值；堆外探测器能够瞬时响应中子通量的变 化，相应的实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N个不 同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可以得到N组 银自给能探测器信号理论值，将理论值与银自给能探测器信号实测值进行比较，取其中符 合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响应 份额； 步骤4、利用Luenberger形式的H c?滤波器对银自给能探测器电流信号作延迟消除： 对于一个离散控制过程系统，该系统可用一个状态方程来描述： X (k+1) = Ax (k) +Bw (k) y (k) = Cx (k) +Dw (k) (8) z (k) = Lx (k) 其中，x(k)为第k次采样点的η维状态向量，w(k)包含了系统过程噪声以及系统观测 白噪声，y(k)为第k次采样点的测量值，z (k)为1维待求向量，L为l*n维矩阵； 针对离散系统（8)，设计如下渐近稳定的满阶线性Luenberger滤波器 = A^+K(yk-C^k) (9) h = 1A 式（9)为H c?滤波器，对于一个给定的γ，当且仅当如下的矩阵不等式有解：其中 Y = YTe RnXn，W e RnXr，J = jTe Rnix' H ①滤波器的增益 K = Y-1W ; 对于银自给能探测器，由其离散状态方程可知方程（8)中的对应矩阵为：通过求解线性矩阵不等式（10)，可得H C?滤波器矩阵K，从而可以由如下步骤获取消除 延迟后任意时刻的探测器电流值： 由初始电流测量值Κ〇)可彳彳，初始〇时刻延迟消除后电流值为 i(0) = CyJc(O); 对于任意k+1 (k = 0, 1，…）时刻，邱+】）=/i,i(/() + .?·(Α)，而k+i时刻延迟消除后 的电流值为辦+丨)=+ D。2. 根据权利要求1所述的基于Luenberger形式H c?滤波的银自给能探测器信号延迟 消除方法，其特征在于，在有换挡的情况下，还包括按如下的信号处理方法对原始信号进行 处理：在换档区域内，假设中子通量保持不变，然后反推中子通量密度产生的电流信号，再 与探测器实际输出电流相减，得到换挡突变分量；在换档区域外，探测器输出电流减去换挡 突变分量，得到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除处理。3. 根据权利要求1所述的基于Luenberger形式H c?滤波的银自给能探测器信号延迟 消除方法，其特征在于，所述的换档区域设计结构如下： 在换档区域内Gc1S k2)，假设中子通量密度不变，则有： n (k+1) = n (k) (11)可以反推出银自给能探测器电流信号为： I (k+1) = p (n (k+1)+X1 (k+1)+x2 (k+1)) (14) 将反推电流（14)当成探测器实际输出电流，通过权利要求1中所述步骤进行延迟消 除； 在换档区域时间边界1^2处，换档引起的电流偏置量可以由下式进行估算： D = l(k2)-y(k2) (15) 其中1 表示在1^2时刻的探测器实际输出电流； 在换档区域外，需要对探测器实际输出电流进行偏置补偿以抵消换档所带来的影响， 将探测器实际输出电流加上式（15)表示的换档引起的电流偏置量，得到中子通量密度产 生的电流信号，然后再对此电流信号进行延迟消除。
【专利摘要】本发明公开了基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法，包括依次进行的以下步骤：步骤1、建立银与热中子的核反应模型；步骤2、采用去耦变换建立核反应模型对应的离散状态方程；步骤3、确定银自给能探测器电流的瞬时响应份额；步骤4、利用Luenberger形式的H∞滤波器对银自给能探测器电流信号作延迟消除。本发明应用时能对银自给能中子探测器的电流信号进行延迟消除处理，并能有效抑制噪声，使得银自给能中子探测器在反应堆瞬态工况时也能正常使用，且由于本发明采用了Luenberger形式的H∞滤波器，作延迟消除时无需预先知道外部扰动输入信号的统计特性。
【IPC分类】G06F19/00, G21C17/108, G01T3/00
【公开号】CN104882173
【申请号】CN201510165895
【发明人】彭星杰, 李庆, 龚禾林, 陈长, 赵文博, 刘启伟, 李向阳, 于颖锐
【申请人】中国核动力研究设计院
【公开日】2015年9月2日
【申请日】2015年4月9日