[发明专利]基于H∞滤波的银自给能探测器信号延迟消除方法

权利要求书

1.基于H∞滤波的银自给能探测器信号延迟消除方法，其特征在于：包括以下步骤：

步骤1、建立银与热中子的核反应模型：

在反应堆瞬态工况下，通量的变化引起银自给能中子探测器电流的变化并不同步，后者较前者有一定的滞后，描述上述反应的具体公式如下：

d108Ag(t)dt=σ107107Agφ(t)-λ108108Ag(t)---(1)]]>

d110Ag(t)dt=σ109109Agφ(t)-λ110110Ag(t)---(2)]]>

I(t)=(K109σ109109Ag+K107σ107107Ag)φ(t)+K110λ110110Ag(t)+K108λ108108Ag(t)---(3)]]>

其中，¹⁰⁸Ag(t)表示¹⁰⁸Ag的核密度，¹¹⁰Ag(t)表示¹¹⁰Ag的核密度，¹⁰⁷Ag表示¹⁰⁷Ag的核密度，¹⁰⁹Ag表示¹⁰⁹Ag的核密度，φ(t)表示探测器处中子通量，¹⁰⁷σ表示¹⁰⁷Ag的中子俘获截面，¹⁰⁹σ表示¹⁰⁹Ag的中子俘获截面，¹⁰⁸λ表示¹⁰⁸Ag的β衰变常数，¹¹⁰λ表示¹¹⁰Ag的β衰变常数，¹⁰⁷K表示¹⁰⁷Ag俘获中子后产生电流的概率，¹⁰⁹K表示¹⁰⁹Ag俘获中子后产生电流的概率，¹⁰⁸K表示¹⁰⁸Ag发生β衰变后产生电流的概率，¹¹⁰K表示¹¹⁰Ag发生β衰变后产生电流的概率，I(t)表示SPND电流；

步骤2、采用去耦变换获取核反应模型对应的离散状态方程：

使用Laplace变换将动态模型化为：

dx1(t)dt=-λ110x1(t)+λ110qpx3(t)---(4)]]>

dx2(t)dt=-λ108x2(t)+λ108rpx3(t)---(5)]]>

I(t)＝p(x₁(t)+x₂(t)+x₃(t))             (6)

其中p为瞬时电流份额，q为¹¹⁰Ag(t)发生β衰变对应的延迟电流份额，r为¹⁰⁸Ag(t)发生β衰变对应的延迟电流份额；

将式(4)、(5)和(6)进行时间离散化处理，并且添加过程噪声项和测量噪声项，可以得到离散状态方程如下：

xk+1=100qp·(1-e-λ108Ts)e-λ108Ts0rp·(1-e-λ110Ts)0e-λ110Ts·xk+100·wk---(7)]]>

I_k＝[p p p]x_k+v_k             (8)

φ_k＝[1 0 0]x_k                  (9)

其中xk=φkx1kx2kT,]]>w_k为过程噪声，v_k为测量噪声，T_S为采样时间。

初始值为

x0=φ0x10x20=I0prI0p2rqI0---(10);]]>

步骤3、确定银自给能探测器电流的瞬时响应份额：

在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆外探测器信号实测值与银自给能探测器信号实测值，堆外探测器能够瞬时响应中子通量的变化，相应的实测值可认为是真实的中子通量，通过调整瞬时响应份额的理论值给定N个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可以得到N组银自给能探测器信号理论值，将理论值与银自给能探测器信号实测值进行比较，取其中符合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响应份额；

步骤4、利用H∞滤波器对银自给能探测器电流信号作延迟消除：

对于一个离散控制过程系统，该系统可用一个状态方程来描述：

x(k+1)＝Ax(k)+Bw(k)

y(k)＝Cx(k)+Dw(k)                 (11)

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)---(12)]]>

系统存在一个H∞滤波器，当且仅当以下的线性矩阵不等式成立

ZZZAZAZB10*YYA+FC+QYA+FCYB1+FD10**ZZ0L′-G′***Y0L′****I0*****γ2I>0---(13)]]>

其中Y、Z为待求解的对称正定矩阵，而Q、G、F为待求解的一般矩阵；

得到上述矩阵后，H∞滤波器的相关矩阵表示如下：

A_f＝-Y^-1Q(I-Y^-1Z)^-1,B_f＝-Y^-1F,C_f＝G(I-Y^-1Z)^-1       (14)

对于银自给能探测器，由其离散状态方程可知方程(11)中的对应矩阵为：

A=100qp·(1-e-λ108Ts)e-λ108Ts0rp·(1-e-λ110Ts)0e-λ110Ts]]>

B=100000]]>

C＝[p p p]

D＝[0 1]

L＝[1 0 0]

通过求解线性矩阵不等式(13)，可得H∞滤波器矩阵A_f、B_f、C_f，从而可以由如下步骤获取消除延迟后任意时刻的探测器电流值：

由初始电流测量值可得x^(0)=y^(0)pr·y^(0)p2rq·y^(0),]]>初始0时刻延迟消除后电流值为z^(0)=Cfx^(0);]]>

对于任意k+1(k＝0,1,...)时刻，而k+1时刻延迟消除后的电流值为z^(k+1)=Cfx^(k+1).]]>