[发明专利]基于卡尔曼滤波的铑自给能探测器信号延迟消除方法

摘要

本发明涉及核反应堆芯测量系统探测器信号处理技术领域，具体公开了一种基于卡尔曼滤波的铑自给能探测器信号延迟消除方法。该方法的具体步骤为：1、建立铑与中子的核反应模型；2、建立卡尔曼滤波模型；3、利用卡尔曼滤波对铑自给能中子探测器电流信号作延迟消除；3.1、获得卡尔曼滤波算法中的系统过程白噪声方差矩阵Q和系统观测白噪声方差矩阵为R；3.2、采集铑自给能探测器电流值，进行模数转换后，利用卡尔曼滤波对铑自给能中子探测器电流信号作延迟消除；本发明所述的一种基于卡尔曼滤波的铑自给能探测器信号延迟消除方法，可以对测量电流信号进行降噪处理，可以保证响应时间足够小的情况下，噪声放大倍数抑制在1~8倍。

主权项

1.一种基于卡尔曼滤波的铑自给能探测器信号延迟消除方法，其特征在于：该方法具体包括如下步骤：步骤1、建立铑与中子的核反应模型；在反应堆瞬态工况下，通量的变化引起铑自给能中子探测器电流的变化并不同步，后者较前者有一定的滞后，描述上述反应的具体公式如下：∂m2(t)∂t=a2n(t)-λ2m2(t)---(1)]]>∂m1(t)∂t=a1n(t)+λ2m2(t)-λ1m1(t)---(2)]]>I(t)＝cn(t)+λ₁m₁(t)                        （3）其中，m₁(t)，m₂(t)代表探测器¹⁰⁴Rh和^104mRh直接引起的电流成分，n(t)表示探测器处中子通量对应的探测器平衡状态下的探测器电流；λ₁，λ₂表示¹⁰⁴Rh和^104mRh的衰变常数；c表示探测器电流的瞬时响应成分；a₁、a₂分别表示¹⁰⁴Rh和^104mRh引起的电流份额；I(t)表示铑自给能电流；步骤2、建立卡尔曼滤波模型；对于一个离散控制过程的系统，该系统可用一个状态方程来描述：X(k+1)=F(k+1|k)X(k)+GW(k)                    （4）系统的测量值：Z(k)=H(k)X(k)+IV(k)                          （5）系统待求量：Y(k)=LX(k)                                   （6）其中，X(k)为第k次采样点的n维状态向量，W(k)系统过程白噪声，其方差矩阵为Q，V(k)为系统观测白噪声，其方差矩阵为R，F(k+1|k)为n*n状态转移矩阵，Z(k)为第k次采样点的m维测量值，H(k)为m*n观测矩阵，G为n*n控制系统矩阵，I为m*n噪声控制矩阵。Y(k)为l维待求向量，L为l*n维矩阵；步骤3、利用卡尔曼滤波对铑自给能中子探测器电流信号作延迟消除；步骤3.1、获得卡尔曼滤波算法中的系统过程白噪声方差矩阵Q和系统观测白噪声方差矩阵为R；获取自给能中子探测器测量电流的白噪声方差R和系统过程白噪声方差矩阵Q；步骤3.2、采集铑自给能探测器电流值，进行模数转换后，利用卡尔曼滤波对铑自给能中子探测器电流信号作延迟消除；对公式（1）、公式（2）、公式（3）作拉普拉斯变换，得到如下等式：I(s)n(s)=c+a1·λ1s+λ1+a2·λ1·λ2s2+(λ1+λ2)·s+λ1·λ2---(7)]]>平衡态时，等式变为：I0n0=c+a1+a2=1---(8)]]>于是公式（7）变为：I(s)=n(s)·I0n0(c+a1·λ1s+λ1+a2·λ1·λ2s2+(λ1+λ2)·s+λ1·λ2))---(9)]]>对公式（9）进行拉普拉斯逆变换，得到如下状态方程：∂x1(t)∂t=1c(a1·λ1-a2·g)·n(t)-λ1x1(t)---(10)]]>∂x2(t)∂t=1ca2·g·n(t)-λ2x2(t)---(11)]]>I(t)＝[c,c,c]·X(t)                        （12）其中，g=λ1·λ2λ1-λ2---(13)]]>X(t)=n(t)x1(t)x2(t)---(14)]]>初始值：X(0)=n(0)1c(a1-a2·g/λ1)·n(0)1c(a2·g)·n(0)---(15)]]>离散状态方程为：X(k+1)=1001c(a1-a2·g/λ1)·(1-eλ1·Ts)e-λ1·Ts01ca2·g·(1-e-λ2·Ts)/λ20e-λ2·Ts·X(k)+100010001·W(k)---(16)]]>I(k)＝[c c c]·X(k)                             （17）n(k)＝[1 0 0]·X(k)                             （18）其中，X(k)=n(k)x1(k)x2(k)]]>初始值：X(0)=I(0)1c(a1-a2·g/λ1)·I(0)1c(a2·g)·I(0)---(19)]]>P(0)=100000000---(20)]]>公式（16）、（17）、（18）分别对应卡尔曼模型中的公式（4）、（5）、（6）。则延迟消除步骤如下：采集k时刻铑自给能探测器电流值I(k)，对应公式（5）的Z(k)，继续下面步骤：进行预估计：XP(k+1)=F(k+1|k)X(k)                            （21）获得协方差矩阵P(k+1)=F(k+1|k)P(k)*F(k+1|k)'+G*Q*G'            （22）获得卡尔曼增益K(k+1)=P(k+1)*H(k)′*(H(k)*P(k+1)*H(k)'+R)-1    （23）获得下一步状态值X(k+1)=XP(k+1)+K(k+1|k)*(Z(k+1)-H(k)*XP(k+1))    （24）更新协方差矩阵P(k+2)=P(k+1)-K(k+1)*H(k)*P(k+1)                 （25）计算获得待求量Y(k)=LX(k)                                       （26）其中，Y(k)即为k时刻铑自给能探测器电流值I（k）对应的延迟消除后的值。