本发明专利技术公开了基于Luenberger形式H2滤波的铑自给能探测器信号延迟消除方法，包括依次进行的以下步骤：步骤1、建立铑与热中子的核反应模型；步骤2、采用直接变换建立核反应模型对应的离散状态方程；步骤3、确定铑自给能探测器电流的瞬时响应份额；步骤4、利用Luenberger形式的H2滤波器对铑自给能探测器电流信号作延迟消除。本发明专利技术应用时能对铑自给能中子探测器的电流信号进行延迟消除处理，并能有效抑制噪声，使得铑自给能中子探测器在反应堆瞬态工况时也能正常使用，且由于本发明专利技术采用了Luenberger形式的H2滤波器，作延迟消除时无需预先知道外部扰动输入信号的统计特性。
【技术实现步骤摘要】
本专利技术设及核反应堆功率分布在线监测系统所用的堆内锭自给能中子探测器信 号的处理技术，具体是基于Luenberger形式&滤波的锭自给能探测器信号延迟消除方法。
技术介绍
用作先进堆巧测量系统堆内探测器的锭自给能中子探测器，其敏感材料锭与中子 反应产生的次生核素发生0衰变产生电流，稳态情况下该电流大小与所在位置通量成正 比，因此通过测量锭自给能探测器能够推知其所在位置中子通量。由于该类探测器电流主 要成分是由次生核素0衰变产生的，在反应堆瞬态情况（中子通量水平变化的情况）下， 该类探测器电流不能实时反映通量水平的变化，而是有一定的延迟，延迟时间参数与次生 核素的0衰变一致。因此，利用锭自给能中子探测器作中子测量装置的先进堆巧测量系 统，为了保证中子通量测量的准确性，需要对锭自给能探器的电流信号作延迟消除处理。 由于实际的测量过程中总伴随有噪声（过程噪声和测量噪声），利用直接的数学 反演方法作延迟消除会将探测器电流信号噪声放大，最大可放大到20倍，影响测量的精 度。因此，在延迟消除处理过程中，需要有效抑制噪声的放大。 目前应用于锭自给能探测器信号延迟的消除主要基于Kalman滤波器实现，其应 用时必须假定系统的外部扰动输入信号是一个具有已知统计特性的白噪声信号，当输入信 号是一个具有有限能量的不确定信号，其统计特性难W得到，该方法就难W应用。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供了一种基于Luenberger形式馬滤 波的锭自给能探测器信号延迟消除方法，其应用时能对锭自给能中子探测器的电流信号进 行延迟消除处理，并能有效抑制噪声，使得锭自给能中子探测器在反应堆瞬态工况时也能 正常使用，且由于本专利技术采用了Luenberger形式的馬滤波器，作延迟消除时无需预先知道 外部扰动输入信号的统计特性。 本专利技术解决上述问题主要通过W下技术方案实现；基于Luenberger形式馬滤波 的锭自给能探测器信号延迟消除方法，其特征在于，包括W下步骤： 步骤1、建立锭与热中子的核反应模型： 在反应堆瞬态工况下，通量的变化引起锭自给能中子探测器电流的变化并不同 步，后者较前者有一定的滞后，描述上述反应的具体公式如下：I(t) =cn(t) +AjIHi(t) (3) 其中，mi(t)、m2(t)分别表示IM化和iMm化直接引起的电荷量，n(t)表示探测器处 热中子通量对应的探测器平衡状态下的探测器电流，A1、A2分别表示化和IMm化的衰 变常数，C表示探测器电流的瞬时响应份额，ai、a2分别表示化和IMm化引起的电流份额， I(t)表示锭自给能电流； 步骤2、采用直接变换获取核反应模型对应的离散状态方程： 令Ja(t)=Aimi(t)代入式（1)、式（2)及式（3)中，将连续时间的常微分方程直 接进行离散化，并添加噪声项，得到W下离散状态方程：[0015：(4)?X(k) + ?V(k) (5) n(k) = ? X(k) (6) 其中，W(k)为过程噪声项，V(k)为测量噪声项，[001引初始值为(7)； 步骤3、确定锭自给能探测器电流的瞬时响应份额： 在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的 堆外探测器信号实测值与锭自给能探测器信号实测值；堆外探测器能够瞬时响应中子通量 的变化，相应的实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N 个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可W得到N 组锭自给能探测器信号理论值，将理论值与锭自给能探测器信号实测值进行比较，取其中 符合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响 应份额； 步骤4、利用Luenberger形式的肥滤波器对锭自给能探测器电流信号作延迟消 除： 对于一个离散控制过程系统，该系统可用一个状态方程来描述：X化+1) =Ax化)+Bw化）y似=Cx似+Dw化） 做 Z似=Lx化） 针对离散系统巧），设计如下渐近稳定的满阶线性Luenberger滤波器巧）[003。 式（9)为最优&滤波器，当且仅当如下的优化问题有解：[003引（；10) 其中J满足如下矩阵不等式： 其中Y=YTgRDXn，WGRnxr，J=jTeR-Xm，肥滤波器的增益K二Y-lW; 对于锭自给能探测器，由其离散状态方程可知方程巧）中的对应矩阵为： 通过求解线性矩阵不等式（10)、（11)、（12)，可得肥滤波器矩阵K，从而可W由如 下步骤获取消除延迟后任意时刻的探测器电流值： 由初始电流测量值W〇)可捐，初始0时刻延迟消除后电流值 为S(0) =C,刮 0);[004引对于任意k+1 化=0,1，...）时刻，种A+1)=.如(/0 + /q.va'）-巧(A))，而k+1 时刻 延迟消除后的电流值为^'(/( +。=。(^ + 0。 本专利技术应用时利用Luenberger形式的肥滤波器原理，在延迟消除过程中，可W有 效地抑制噪声的放大，噪声抑制效果越好，延迟效果会逐渐变差，因此，本专利技术应用时需适 当调节参数使延迟消除效果与噪声抑制达到最佳平衡。 当需要探测较大动态范围的中子通量密度，相应的也需要检测大动态范围的电流 信号，而该一问题便集中在了模数转换器上。为了适应大动态范围的电流的量化，锭自给能 探测器的模数转换器采样分档电阻，当电流信号在大范围变化时，模数转换器就会发生电 阻档位转换。由于各档位没有完全匹配，各档位之间的切换会造成输出信号的近似于阶跃 的突变。 在有换挡的情况下，还包括按如下的信号处理方法对原始信号进行处理；在换档 区域内，假设中子通量保持不变，然后反推中子通量密度产生的电流信号，再与探测器实际 输出电流相减，得到换挡突变分量；在换档区域外，探测器输出电流减去换挡突变分量，得 到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除处理。 换挡引起的突变分量进入延迟消除模块后，会被严重放大，使得时域上的阶跃突 变被严重放大，影响最终信号延迟消除的质量即突变部分信号的严重失真。在换挡时间段 里，信号的变化主要由换挡突变贡献，相对而言，由中子通量密度变化引起的电流信号变化 可W忽略。 为了处理换档导致的信号突变问题，进一步的，所述的换档区域设计结构如下： [005。 在换档区域内（ki《k《k2)，假设中子通量密度不变，则有： n (k+1) = n (k) (13) 可W反推出锭自给能探测器电流信号为： I 化+1) = Ja 化+1)+cn 化+1) (16) 将反推电流（16)当成探测器实际输出电流，通过权利要求1中所述步骤进行延迟 消除；[005引在换档区域时间边界k,处，换档引起的电流偏置量可W由下式进行估算： D=/(L_)-y(L_)(17) 其中.P(年)表示在k2时刻的探测器实际输出电流；在换档区域外，需要对探测器实 际输出电流进行偏置补偿W抵消换档所带来的影响，将探测器实际输出电流加上式（17) 表示的换档引起的电流偏置量，得到中子通量密度产生的电流信号，然后再对此电流信号 进行延迟消除。 综上所述，本专利技术具有W下有益效果： 1本专利技术整体工序简单，便于实现，能对锭自给能中子探测器的电流信号进行延迟本文档来自技高网...
【技术保护点】
基于Luenberger形式H2滤波的铑自给能探测器信号延迟消除方法，其特征在于：包括以下步骤：步骤1、建立铑与热中子的核反应模型：在反应堆瞬态工况下，通量的变化引起铑自给能中子探测器电流的变化并不同步，后者较前者有一定的滞后，描述上述反应的具体公式如下：∂m2(t)∂t=a2n(t)-λ2m2(t)---(1)]]>∂m1(t)∂t=a1n(t)+λ2m2(t)-λ1m1(t)---(2)]]>I(t)＝cn(t)+λ1m1(t) (3)其中，m1(t)、m2(t)分别表示104Rh和104mRh直接引起的电荷量，n(t)表示探测器处热中子通量对应的探测器平衡状态下的探测器电流，λ1、λ2分别表示104Rh和104mRh的衰变常数，c表示探测器电流的瞬时响应份额，a1、a2分别表示104Rh和104mRh引起的电流份额，I(t)表示铑自给能电流；步骤2、采用直接变换获取核反应模型对应的离散状态方程：令Ja(t)＝λ1m1(t)代入式(1)、式(2)及式(3)中，将连续时间的常微分方程直接进行离散化，并添加噪声项，得到以下离散状态方程：X(k+1)=e-λ1·Tsλ2·(1-e-λ1·Ts)a1·(1-e-λ1·Ts)0e-λ2·Tsa2λ2·(1-e-λ2·Ts)001·X(k)+001W(k)---(4)]]>I(k)＝[1 0 c]·X(k)+[1]·V(k) (5)n(k)＝[0 0 1]·X(k) (6)其中，X(k)=Ja(k)m2(k)n(k),]]>W(k)为过程噪声项，V(k)为测量噪声项，初始值为X(0)=Ja(0)m2(0)n(0)=(a1+a2)·I(0)a2λ2·I(0)I(0)---(7);]]>步骤3、确定铑自给能探测器电流的瞬时响应份额：在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆外探测器信号实测值与铑自给能探测器信号实测值；堆外探测器能够瞬时响应中子通量的变化，相应的实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可以得到N组铑自给能探测器信号理论值，将理论值与铑自给能探测器信号实测值进行比较，取其中符合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响应份额；步骤4、利用Luenberger形式的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维矩阵；针对离散系统(8)，设计如下渐近稳定的满阶线性Luenberger滤波器x^k+1=Ax^k+K(yk-Cx^k)z^k=Lx^k---(9)]]>式(9)为最优H2滤波器，当且仅当如下的优化问题有解：minJ,Y,WTr{J}---(10)]]>其中J满足如下矩阵不等式：JBTY-DTWTYB-WDY≥0---(11)]]>YYA-WC0ATY-CTWTYLT0LI≥0---(12)]]>其中Y＝YT∈Rn×n，W∈Rn×r，J＝JT∈Rm×m，H2滤波器的增益K＝Y‑1W；对于铑自给能探测器，由其离散状态方程可知方程(8)中的对应矩阵为：A=e-λ1·Tsλ2·(1-e-λ1·Ts)a1·(1-e-λ1·Ts)0e-λ2·Tsa2λ2·(1-e-λ2·Ts)001]]>B=000010]]>C＝[1 0 c]D＝[0 1]L＝[0 0 1]通过...
【技术特征摘要】
【专利技术属性】
技术研发人员：彭星杰，李庆，龚禾林，陈长，赵文博，刘启伟，李向阳，于颖锐，
申请(专利权)人：中国核动力研究设计院，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。