本发明专利技术公开了基于H∞滤波消除钒自给能探测器信号延迟的方法，包括依次进行的以下步骤：步骤1、建立钒与热中子的核反应模型；步骤2、采用直接变换建立核反应模型对应的离散状态方程；步骤3、确定钒自给能探测器电流的瞬时响应份额；步骤4、利用H∞滤波器对钒自给能探测器电流信号作延迟消除。本发明专利技术应用时能对钒自给能中子探测器的电流信号进行延迟消除处理，并能有效抑制噪声，使得钒自给能中子探测器在反应堆瞬态工况时也能正常使用，且由于本发明专利技术采用了H∞滤波器，作延迟消除时无需预先知道外部扰动输入信号的统计特性。
【技术实现步骤摘要】
本专利技术设及先进堆巧测量系统（核反应堆功率分布在线监测系统）所用的堆内饥 自给能中子探测器信号的处理技术，具体是基于H-滤波的饥自给能探测器信号延迟消除 方法。
技术介绍
用作先进堆巧测量系统堆内探测器的饥自给能中子探测器，其敏感材料饥与中子 反应产生的次生核素发生0衰变产生电流，稳态情况下该电流大小与所在位置通量成正 比，因此通过测量饥自给能探测器能够推知其所在位置中子通量。由于该类探测器电流主 要成分是由次生核素0衰变产生的，在反应堆瞬态情况（中子通量水平变化的情况）下， 该类探测器电流不能实时反映通量水平的变化，而是有一定的延迟，延迟时间参数与次生 核素的0衰变一致。因此，利用饥自给能中子探测器作中子测量装置的先进堆巧测量系 统，为了保证中子通量测量的准确性，需要对饥自给能探器的电流信号作延迟消除处理。 由于实际的测量过程中总伴随有噪声（过程噪声和测量噪声），利用直接的数学 反演方法作延迟消除会将探测器电流信号噪声放大，最大可放大到20倍，影响测量的精 度。因此，在延迟消除处理过程中，需要有效抑制噪声的放大。 目前应用于饥自给能探测器信号延迟的消除主要基于Kalman滤波器实现，其应 用时必须假定系统的外部扰动输入信号是一个具有已知统计特性的白噪声信号，当输入信 号的统计特性难W得到时，该方法就难W应用。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供了一种基于H-滤波的饥自给能探 测器信号延迟消除方法，其应用时能对饥自给能中子探测器的电流信号进行延迟消除处 理，并能有效抑制噪声，使得饥自给能中子探测器在反应堆瞬态工况时也能正常使用，且由 于本专利技术采用了H-滤波器，作延迟消除时无需预先知道外部扰动输入信号的统计特性。 本专利技术解决上述问题主要通过W下技术方案实现；基于H-滤波的饥自给能探测 器信号延迟消除方法，包括W下步骤： 步骤1、建立饥与热中子的核反应模型： 在反应堆瞬态工况下，通量的变化引起饥自给能中子探测器电流的变化并不同 步，后者较前者有一定的滞后，描述上述反应的具体公式如下： (1) 其中，巫（t)为中子通量，Nsi(t)为Sly的核密度，Ns2(t)为52y的核密度，〇51为Sly 的中子吸收截面，人52为52V的衰变常量，I(t)为饥自给能探测器输出电流，Kpy为瞬时成分 的灵敏度，Kgy为瞬时成分的灵敏度； 将式（1)变形得到如下等式：I(t) =cn(t) +Anig(t) (4)[001引其中，mi(t)、m2(t)分别代表51v，52v的存量，c为电流的瞬时份额，a为电流的e衰变份额； 步骤2、采用直接变换获取核反应模型对应的离散状态方程： 令L(t) =Am2(t)代入式（2)、式做及式（4)中，将连续时间的常微分方程直 接进行离散化，并添加噪声项，得到W下离散状态方程：[00过初始值为C8) 步骤3、确定饥自给能探测器电流的瞬时响应份额； 步骤4、利用H-滤波器对饥自给能探测器电流信号作延迟消除： 对于一个离散控制过程系统，该系统可用一个状态方程来描述： (9)[002引其中，x(k)为第k次采样点的n维状态向量，w似包含了系统过程噪声化及系统 观测白噪声，y化）为第k次采样点的测量值，Z化）为1维待求向量，L为l*n维矩阵； 假定系统是渐近稳定的，则对给定的常数丫〉0,要求设计一个渐近稳定的满阶线 性滤波器(10) 系统存在一个H-滤波器，当且仅当W下的线性矩阵不等式成立 「003引(11) 其中Y、Z为待求解的对称正定矩阵，而Q、G、F为待求解的一般矩阵； 得到上述矩阵后，H-滤波器的相关矩阵表示如下：Af=-Y-iQ(I-Y-i幻-1，Bf=-Y-中，Cf=G(I-Y-1 幻-1 (12) 对于饥自给能探测器，由其离散状态方程可知方程巧）中的对应矩阵为： D = [00川 L = [0 U 通过求解线性矩阵不等式（11)，可得滤波器矩阵4,、8,八,，从而可^由如下步 骤获取消除延迟后任意时刻的探测器电流值： 由初始电流测量值W〇)可得的始0时刻延迟消除后电流值为 i(〇) =C,i-(0)； 对于任意k+1化=0, 1，...）时刻，；a' + 1) = /!,刮/。+公,.v''(A')，而k+1时刻延迟消 除后的电流值为维+1) =c, -W: +U。 核反应模型是应用滤波器进行延迟消除的基础，本专利技术由第一性原理出发，步骤 一导出饥自给能探测器产生信号该一物理过程所对应的连续时间变量数学模型。由于探测 器的电流信号都是通过离散采样获取，步骤2将步骤1建立的连续状态方程转换为离散状 态方程。 本专利技术应用时利用H-滤波器原理，在延迟消除过程中，可W有效地抑制噪声的 放大，噪声抑制效果越好，延迟效果会逐渐变差，因此，本专利技术应用时需适当调节参数使延 迟消除效果与噪声抑制达到最佳平衡。 饥自给能探测器瞬时响应份额C可W由理论进行大致地估算，但理论估算值与真 实值之间的不匹配将会导致滤波效果的下降，为了准确地确定出瞬时响应份额C，进一步 的，所述步骤3确定饥自给能探测器电流的瞬时响应份额具体包括如下步骤；在反应堆启 动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆外探测器信号实测 值与饥自给能探测器信号实测值，其中，堆外探测器能够瞬时响应中子通量的变化，相应的 实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N个不同的瞬时响 应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可W得到N组饥自给能探 测器信号理论值，将理论值与饥自给能探测器信号实测值进行比较，取其中符合程度最好 的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响应份额。 当需要探测较大动态范围的中子通量密度，相应的也需要检测大动态范围的电流 信号，而该一问题便集中在了模数转换器上。为了适应大动态范围的电流的量化，饥自给能 探测器的模数转换器采样分档电阻，当电流信号在大范围变化时，模数转换器就会发生电 阻档位转换。由于各档位没有完全匹配，各档位之间的切换会造成输出信号的近似于阶跃 的突变。 换挡引起的突变分量进入延迟消除模块后，会被严重放大，使得时域上的阶跃突 变被严重放大，影响最终信号延迟消除的质量（突变部分信号的严重失真）。在换挡时间段 里，信号的变化主要由换挡突变贡献，相对而言，由中子通量密度变化引起的电流信号变化 可W忽略。 为了处理换档导致的信号突变问题，进一步的，在有换挡的情况下，还包括按如下 的信号处理方法对原始信号进行处理；在换档区域内，假设中子通量保持不变，然后反推中 子通量密度产生的电流信号，再与探测器实际输出电流相减，得到换挡突变分量；在换档区 域外，探测器输出电流减去换挡突变分量，得到中子通当前第1页1 2 3 本文档来自技高网...
【技术保护点】
基于H∞滤波的钒自给能探测器信号延迟消除方法，其特征在于：包括以下步骤：步骤1、建立钒与热中子的核反应模型：在反应堆瞬态工况下，通量的变化引起钒自给能中子探测器电流的变化并不同步，后者较前者有一定的滞后，描述上述反应的具体公式如下：ddtN51(t)=-N51(t)σ51Φ(t)ddtN52(t)=N51(t)σ51Φ(t)-λ52N52(t)I(t)=KpvN51(t)σ51Φ(t)+Kgvλ52N52(t)---(1)]]>其中，Φ(t)为中子通量，N51(t)为51V的核密度，N52(t)为52V的核密度，σ51为51V的中子吸收截面，λ52为52V的衰变常量，I(t)为钒自给能探测器输出电流，Kpv为瞬时成分的灵敏度，Kgv为瞬时成分的灵敏度；将式(1)变形得到如下等式：∂m1(t)∂t=-an(t)---(2)]]>∂m2(t)∂t=an(t)-λm2(t)---(3)]]>I(t)＝cn(t)+λm2(t) (4)其中，m1(t)、m2(t)分别代表51V，52V的存量，c为电流的瞬时份额，a为电流的β衰变份额；步骤2、采用直接变换获取核反应模型对应的离散状态方程：令Ja(t)＝λm2(t)代入式(2)、式(3)及式(4)中，将连续时间的常微分方程直接进行离散化，并添加噪声项，得到以下离散状态方程：X(k+1)=e-λ·Ts(1-e-λ·Ts)01·X(k)+01W(k)---(5)]]>I(k)＝[1 c]·X(k)+[1]·V(k) (6)n(k)＝[0 1]·X(k) (7)其中，X(k)=Ja(k)n(k)]]>初始值为X(0)=Ja(0)n(0)=I(0)I(0)---(8)]]>步骤3、确定钒自给能探测器电流的瞬时响应份额；步骤4、利用H∞滤波器对钒自给能探测器电流信号作延迟消除：对于一个离散控制过程系统，该系统可用一个状态方程来描述：x(k+1)＝Ax(k)+Bw(k)y(k)＝Cx(k)+Dw(k) (9)z(k)＝Lx(k)其中，x(k)为第k次采样点的n维状态向量，w(k)包含了系统过程噪声以及系统观测白噪声，y(k)为第k次采样点的测量值，z(k)为1维待求向量，L为l*n维矩阵；假定系统是渐近稳定的，则对给定的常数γ>0，要求设计一个渐近稳定的满阶线性滤波器x^(k+1)=Afx^(k)+Bfy^(k)z^(k)=Cfx^(k)---(10)]]>系统存在一个H∞滤波器，当且仅当以下的线性矩阵不等式成立ZZZAZAZB0*YYA+FC+QYA+FCYB+FD0**ZZ0L′-G′***Y0L′****I0*****γ2I>0---(11)]]>其中Y、Z为待求解的对称正定矩阵，而Q、G、F为待求解的一般矩阵；得到上述矩阵后，H∞滤波器的相关矩阵表示如下：Af＝‑Y‑1Q(I‑Y‑1Z)‑1,Bf＝‑Y‑1F,Cf＝G(I‑Y‑1Z)‑1 (12)对于钒自给能探测器，由其离散状态方程可知方程(9)中的对应矩阵为：A=e-λ·Ts(1-e-λ·Ts)01]]>B=0010]]>C＝[1 c]D＝[0 1]L＝[0 1]通过求解线性矩阵不等式(11)，可得H∞滤波器矩阵Af、Bf、Cf，从而可以由如下步骤获取消除延迟后任意时刻的探测器电流值：由初始电流测量值可得x^(0)=y^(0)y^(0),]]>初始0时刻延迟消除后电流值为z^(0)=Cfx^(0);]]>对于任意k+1(k＝0,1,...)时刻，而k+1时刻延迟消除后的电流值为z^(k+1)=Cfx^(k+1).]]>...
【技术特征摘要】
【专利技术属性】
技术研发人员：龚禾林，陈长，彭星杰，赵文博，刘启伟，李向阳，李庆，于颖锐，
申请(专利权)人：中国核动力研究设计院，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。