MCP成像探测器前端电子学增益补偿

窦双团, 付利平, 贾楠, 王天放

窦双团, 付利平, 贾楠, 王天放. MCP成像探测器前端电子学增益补偿[J]. 应用光学, 2022, 43(6): 1165-1174. DOI: 10.5768/JAO202243.0604017
引用本文: 窦双团, 付利平, 贾楠, 王天放. MCP成像探测器前端电子学增益补偿[J]. 应用光学, 2022, 43(6): 1165-1174. DOI: 10.5768/JAO202243.0604017
DOU Shuangtuan, FU Liping, JIA Nan, WANG Tianfang. Research on front-end electronics gain compensation of MCP imaging detector[J]. Journal of Applied Optics, 2022, 43(6): 1165-1174. DOI: 10.5768/JAO202243.0604017
Citation: DOU Shuangtuan, FU Liping, JIA Nan, WANG Tianfang. Research on front-end electronics gain compensation of MCP imaging detector[J]. Journal of Applied Optics, 2022, 43(6): 1165-1174. DOI: 10.5768/JAO202243.0604017

MCP成像探测器前端电子学增益补偿

基金项目: 国家自然科学基金(42174226, 41874187)
详细信息
    作者简介:

    窦双团(1992—),男,博士研究生,主要从事真空紫外光学遥感单光子成像探测技术研究。E-mail:doushuangtuan@163.com

    通讯作者:

    付利平(1972—),女,研究员,博士生导师,主要从事真空紫外光学遥感研究。E-mail:fuliping@nssc.ac.cn

  • 中图分类号: TN366

Research on front-end electronics gain compensation of MCP imaging detector

  • 摘要:

    基于微通道板电子倍增电荷分割型阳极的成像探测器常用于行星大气、气辉等微弱信号探测。针对探测器读出电路增益不相等导致探测器成像产生畸变的问题,利用理论分析结合MATLAB仿真揭示了不同畸变图像的形成机制,在该基础上提出了一种探测器读出电路增益补偿方法减小探测器成像畸变。通过MATLAB仿真和实验测试结果表明该文提出的方法能够有效减小由于读出电路增益不相等导致的探测器成像畸变。

    Abstract:

    The imaging detectors based on the electron multiplication charge splitting anode of micro-channel plate are often used for weak signal detection such as planetary atmosphere and airglow. In view of the problem of detector imaging distortion caused by unequal gain of detector readout circuit, the theoretical analysis combined with MATLAB simulations were used to reveal the formation mechanism of different distorted images, and a gain compensation method of detector readout circuit was proposed to reduce the detector imaging distortion. The MATLAB simulation and experimental results show that the proposed method can effectively reduce the detector imaging distortion due to the unequal gain of readout circuit.

  • 基于微通道板 (micro-channel plate, MCP)电子倍增电荷分割型阳极的成像探测器常用于行星大气气辉、极光等微弱信号探测,如TIMED卫星的GUVI载荷、2019年发射的ICON卫星的远紫外波段( far ultraviolet, FUV)载荷等[1-4],我国卫星计划包括FY-4电离层成像光谱仪、极光成像仪[5]、嫦娥7号的极紫外相机、包括处于规划中的太阳系边际项目、木星探测等项目均采用基于MCP电子倍增电荷分割阳极的成像探测器对紫外辐射进行探测。由此可见,该探测器在深空紫外辐射探测方面有着极广的应用前景。随着空间天气研究的深入,对成像探测器的成像质量提出了更高的要求,如更高的空间分辨率和光谱分辨率以及更小的图像畸变等。探测器成像畸变不仅影响对目标图像的分析,更有可能造成对空间天气相关事件的误判。因此,更小的图像畸变是实现提高探测器成像质量的关键指标。探测器常见的畸变有“S”型畸变、调制畸变[6]。在我们的研究中发现,探测器图像还存在另外一种容易被忽略的由于探测器读出电路增益不相等导致的畸变,此类型的畸变会降低探测器的空间分辨率。就笔者目前接触到的文献中,西安光机所的缪振华在其博士论文[7]中指出,探测器读出电路增益不相等会降低探测器的空间分辨,但是读出电路增益不相等对探测器空间分辨率影响的产生机制以及工程中如何降低或者消除此影响并未做进一步论述。本文以MCP电子倍增(wedge strip anode,WSA)电荷分割型阳极成像探测器为研究对象,以理论分析结合MATLAB仿真的研究方法,研究了读出电路增益不相等对探测器成像的影响,在此基础上提出了一种增益补偿方法,仿真和实验结果表明,该方法能够有效减小探测器读出电路增益偏差导致的图像畸变,提高探测器成像质量。

    基于MCP电子倍增WSA电荷分割型阳极的极远紫外成像探测器主要由光阴极、MCP堆、WSA电荷分割型阳极、读出电路以及数据采集系统等组成[8-10]。当目标波段的光子入射到探测器时,位于探测器前端的光阴极以一定的量子效率将入射光子转换为电子后被MCP堆倍增为106~107量级的电子云经MCP堆出射端口输出,从MCP堆出射的电子云被电荷分割型阳极的不同电极收集后由读出电路进一步处理。用于分割从MCP堆输出电子云的WSA阳极结构,如图1所示。

    图  1  WSA阳极不同电极结构示意图
    Figure  1.  Schematic diagram of different electrode structures in WSA anode

    图1所示,WSA阳极由W(wedge)、S(strip)、Z(zigzag)3个电极构成。从MCP堆出射的电子云被WSA阳极的3个电极以不同比例收集,利用这一特性可对电子云质心坐标进行编码[11-12],根据不同电极收集到的电荷量确定电子云质心坐标的公式如下所示。

    $$ \left\{ \begin{gathered} x = \frac{{2{Q_{\rm{s}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} \\ y = \frac{{2{Q_{\rm{w}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} \\ \end{gathered} \right. $$ (1)

    式中:${Q_{\rm{w}}}$、${Q_{\rm{s}}}$、${Q_{\rm{z}}}$分别代表电极$W$、$S$、$Z$收集到的电荷量;$x$和$y$分别代表电子云质心的横坐标和纵坐标。根据方程(1),测量3个电极收集的电荷量就可以反演出电子云质心坐标。实际电荷信号不便于直接测量,需要3路读出电路分别将3个电极收集到的电荷信号转换为易于处理的电压信号,通过建立电压信号和被收集电荷量的联系,间接确定电子云质心坐标。

    探测器读出电路主要由电荷灵敏前置放大器、零极相消电路、半高斯整形滤波电路组成。其中,电荷灵敏前置放大器将探测器阳极输出的电流脉冲转换为负指数衰减的电压脉冲;零极相消电路将前置放大器输出的负指数衰减脉冲的拖尾缩短,减少系统的脉冲堆积;滤波整形电路将零极相消电路输出的脉冲整形为便于后续数据采集系统采集的准高斯脉冲,同时滤除前级电路输出信号的噪声分量;探测器读出电路的信号传输模型和各级电路输出的信号波形示意图如图2所示(理论要求探测器的3路读出电路完全一致,因此,图2所示的读出电路信号传输模型对探测器的3路读出电路均适用)。

    图  2  读出电路信号传输模型和波形示意图
    Figure  2.  Signal transmission model of readout circuit and waveform diagram

    图2所示探测器读出电路信号传输模型从左到右分别表示:图2(a)探测器阳极输出电流脉冲的时域波形示意图以及时域模型,其中${Q_{{\rm{in}}}}$表示WSZ阳极任意一个电极收集的电荷量;图2(b)、(c)、(d)分别表示探测器读出电路的电荷灵敏前置放大器、零极相消电路、半高斯整形滤波器电路输出波形示意图以及相应的S域模型。图2(b)中的${C_f}$表示前置放大器用于收集阳极输出电荷的积分电容,${\tau _f}$表示前置放大器的时间常数,其决定了图2(b)所示的负指数衰减脉冲的拖尾长度;图2(c)中的${\tau _{p{\textit{z}}}}$和$ {\tau _c} $分别表示为了缩短前置放大器的脉冲拖尾,用于和${\tau _f}$相抵消的时间常数以及零极相消电路输出脉冲的时间常数;图2(d)中的${\tau _i}$表示半高斯整形滤波电路输出准高斯脉冲的时间常数,n表示滤波整形电路的整形阶数。根据图2所示的信号模型,探测器读出电路传递函数的S域模型如下所示。

    $$ \dfrac{{{V_{{\rm{out}}}}(s)}}{{{i_{\rm{D}}}(s)}} = \dfrac{1}{{{C_f}}} \times \dfrac{1}{{\left(s + \dfrac{1}{{{\tau _f}}}\right)}} \times \dfrac{{\left(s + \dfrac{1}{{{\tau _{p{\rm{z}}}}}}\right)}}{{\left(s + \dfrac{1}{{{\tau _c}}}\right)}} \times \dfrac{1}{{{{(1 + s{\tau _i})}^n}}} $$ (2)

    式中:${V_{{\rm{out}}}}(s)$和${i_{\rm{D}}}(s)$分别表示读出电路输出的准高斯电压脉冲和探测器阳极输出的窄电流脉冲的S域模型,其中电流脉冲${i_{\rm{D}}}(s)$对应的时域模型为:${i_{\rm{D}}}(s) = {Q_{{\rm{in}}}}\delta (t)$。实际电路中为方便处理,通常可令${\tau _c} = {\tau _i} = \tau $、${\tau _{p{\rm{z}}}} = {\tau _f}$,将${\tau _c} = {\tau _i} = \tau $和${\tau _{p{\rm{z}}}} = {\tau _f}$带入(2)式,对(2)式进行逆拉普拉斯变换,整理后可得读出电路输出准高斯脉冲时域表达式和准高斯脉冲的时间导数分别如(3)、(4)式所示。

    $$ {V_{{\rm{out}}}} = \frac{{{Q_{{\rm{in}}}}}}{{{C_f}}} \times \frac{1}{{n!}} \times {\left(\frac{t}{{{\tau _i}}}\right)^n} \times {{\rm{e}}^{ - \tfrac{t}{{{\tau _i}}}}} $$ (3)
    $$ \frac{{{\rm{d}}{V_{{\rm{out}}}}(t)}}{{{\rm{d}}t}} = \frac{{{Q_{{\rm{in}}}}}}{{{C_f}}} \times \frac{1}{{n!}} \times {\left(\frac{1}{{{\tau _i}}}\right)^n} \times \left(n \times {t^{n - 1}}{{\rm{e}}^{ - \tfrac{t}{{{\tau _i}}}}} - \frac{{{t^n}}}{{{\tau _i}}}{{\rm{e}}^{ - \tfrac{t}{{{\tau _i}}}}}\right) $$ (4)

    为了得到读出电路输出准高斯脉冲的峰值时间${t_p}$,令(4)式求导结果等于0,可解得峰值时间${t_p} = n\tau $,将峰值时间${t_p}$带入(3)式可以解得读出电路输出准高斯脉冲的电压峰值${V_p}$为

    $$ {V_p} = \frac{{{Q_{{\rm{in}}}}}}{{{C_f}}} \times \frac{1}{{n!}} \times {(n)^n} = g \times {Q_{{\rm{in}}}} $$ (5)

    式中:$g$表示探测器读出电路的增益,(5)式描述了探测器读出电路输出准高斯脉冲峰值电压${V_p}$和阳极电极收集的电荷量${Q_{{\rm{in}}}}$的关系,由(5)式可知,读出电路的增益$g$只和探测器前置放大器的积分电容${C_f}$以及滤波整形电路的整形阶数$n$有关。

    为了建立读出电路输出准高斯脉冲的峰值电压和阳极电极输出电荷量的联系,分别将探测器的3路读出电路的积分电容${C_{f\_w}}$、${C_{f\_s}}$、${C_{f\_{\textit{z}}}}$以及阳极的W、S、Z三个电极收集到的电荷${Q_w}$、${Q_s}$、${Q_{\textit{z}}}$分别带入(5)式整理后得到探测器3路读出电路输出准高斯脉冲电压峰值分别为

    $$ Q\left\{ \begin{gathered} {V_{\rm{w}}} = \frac{{{Q_{\rm{w}}}}}{{{C_{f\_{\rm{w}}}}}} \times \frac{1}{{n!}} \times {(n)^n} \times {{\rm e}^{ - n}} = {Q_{\rm{w}}} \times h \\ {V_{\rm{s}}} = \frac{{{Q_{\rm{s}}}}}{{{C_{f\_{\rm{s}}}}}} \times \frac{1}{{n!}} \times {(n)^n} \times {{\rm e}^{ - n}} = {Q_{\rm{s}}} \times m \\ {V_{\rm{z}}} = \frac{{{Q_{\textit{z}}}}}{{{C_{f\_{\textit{z}}}}}} \times \frac{1}{{n!}} \times {(n)^n} \times {{\rm e}^{ - n}} = {Q_{\textit{z}}} \times k \\ \end{gathered} \right. $$ (6)

    方程组(6)的$h$、$m$、$k$分别表示W、S、Z三个电极的各自读出电路的增益。假设读出电路的增益系数$h = m = k$,利用方程组(6)中3路读出电路的电压峰值和各自电极收集电荷量的关系带入公式(1),整理后可将电子云质心坐标${p_{\rm e}}(x,y)$表示为

    $$ {p_{\rm e}}(x,y) = \left\{ \begin{gathered} x = \frac{{2{Q_{\rm{s}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} = \frac{{2{V_{\rm{s}}}}}{{{V_{\rm{w}}} + {V_{\rm{s}}} + {V_{\rm{z}}}}} \\ y = \frac{{2{Q_{\rm{w}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} = \frac{{2{V_{\rm{w}}}}}{{{V_{\rm{w}}} + {V_{\rm{s}}} + {V_{\rm{z}}}}} \\ \end{gathered} \right. $$ (7)

    (7)式即为WSA阳极常用的电子云质心解码公式,由于实际电子元器件的各向异性以及温度漂移等原因导致3路读出电路的增益$h$、$m$、$k$不相等,忽略读出电路增益偏差,直接使用公式(7)对电子云质心坐标解码,必然导致解码坐标和真实坐标存在偏差,进而影响探测器成像质量,为此,将3路读出电路的增益$h$、$m$、$k$考虑进解码公式(7),整理后可以得到:

    $$ {p_{\rm e}}(x,y) = \left\{ \begin{gathered} x = \frac{{2{Q_{\rm{s}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} = \frac{{2hk{V_{\rm{s}}}}}{{mk{V_{\rm{w}}} + hk{V_{\rm{s}}} + mh{V_{\rm{z}}}}} \\ y = \frac{{2{Q_{\rm{w}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} = \frac{{2mk{V_{\rm{w}}}}}{{mk{V_{\rm{w}}} + hk{V_{\rm{s}}} + mh{V_{\rm{z}}}}} \\ \end{gathered} \right. $$ (8)

    理论上根据公式(8)可以得到电子云质心的真实坐标。但是,由于实际电路的复杂性,有时准确求出增益因子$h$、$m$、$k$难以实现,即使求出了$h$、$m$、$k$的值,由于元器件的温度漂移或者其他未知干扰也会使得增益$h$、$m$、$k$生变化,造成电子云质心解码误差。因此需要一种不依赖读出电路增益却能尽量减小电子云质心解码坐标和真实坐标偏差的解码方法。

    如第一节所述,读出电路增益不相等将导致电子云质心的解码坐标偏离真实坐标、造成解码误差,进而影响探测器成像质量。接下来将具体分析读出电路增益不相等对读出电路输出的准高斯脉冲峰值、峰值相位以及对探测器成像的影响。

    如(5)式所描述,读出电路增益$g$完全由积分电容${C_f}$、整形阶数$n$决定。为了测试探测器的3路读出电路增益不相等对准高斯脉冲峰值的影响,设W、S、Z电极输出等量电荷到各自的读出电路、读出电路适用的积分电容${C_f}$的精度为±20%(理论要求3路读出电路的前置放大器的积分电容${C_{f\_w}} = {C_{f\_s}} = {C_{f\_{\textit{z}}}} = 1\;{\rm{pF}}$,实际电子元器件都存在精度等级,±20%为常用的电子元器件精度等级,因此3路读出电路前置放大器的积分电容取值范围为:0.8 pF$ \leqslant {C_f} \leqslant $1.2 pF)、整形时间${\tau _i}$等于10 μs。根据上述条件,基于公式(5)仿真了±20%精度等级的积分电容以及不同的整形阶数($1 \leqslant n \leqslant 4$,实际电路中,整形阶数$n$通常不大于4)对读出电路输出准高斯脉冲峰值的影响。由积分电容$ {C}_{f} $的精度等级和整形阶数$n$的变化引起读出电路输出脉冲峰值变化的图像如图3(a)所示,固定$n$等于4(仿真模型使用的整形阶数为4阶),只考虑积分电容${C_f}$精度引起输出脉冲变化的图像如图3(b)所示。

    图  3  积分电容和整形阶数对输出准高斯脉冲峰值的影响
    Figure  3.  Effects of integral capacitance and shaping orders on quasi-Gaussian pulse peak

    图3(a)所示,读出电路输出准高斯脉冲的峰值和积分电容${C_f}$以及整形阶数$ n $成反比,拥有更小积分电容和更低整形阶数的读出电路输出的准高斯脉冲峰值更高。实际使用中,同一探测器不同电极读出电路使用的整形阶数是相同的。如图3(b)所示,固定整形阶数$n$等于4,读出电路输出准高斯脉冲峰值只和电容${C_f}$的大小有关。电容${C_f}$精度等级为±20%时,对于等量输入电荷,探测器的3路读出电路输出脉冲峰值最大值和最小值的差异率$dr$为

    $$\begin{split} dr = & \dfrac{{{V_{\max }} - {V_{\min }}}}{{{V_{\max }}}} \times 100{\text{%}} = \dfrac{{\dfrac{{{Q_{\rm in}}}}{{{C_{f0}}(1 - 20{\text{%}} )}} - \dfrac{{{Q_{\rm in}}}}{{{C_{f0}}(1 + 20{\text{%}} )}}}}{{\dfrac{{{Q_{\rm in}}}}{{{C_{f0}}(1 - 20{\text{%}} )}}}} =\\ &33{\text{%}}\\[-12pt] \end{split} $$ (9)

    式中:${C_{f0}}$表示积分电容的标称值,当积分电容${C_{f0}}$存在±20%的误差即${C_{f0}}(1 - 20{\text{%}} ) \leqslant {C_f} \leqslant {C_{f0}}(1 + 20{\text{%}} )$时,探测器3路读出电路输出脉冲峰值的最大值与最小值之差为最大值的33%。因此,忽略读出电路增益偏差,使用公式(7)对电子云质心坐标解码将会使解码坐标偏离真实坐标,降低探测器成像质量。

    读出电路输出准高斯脉冲峰值${V_p}$对应的峰值时间${t_p} = n\tau $,因此,准高斯脉冲的峰值时间${t_p}$及其对应的峰值相位取决于整形时间常数${\tau _i}$和整形阶数$n$。和电容的精度等级相同,假设读出电路的电阻也存在±20%的误差,则整形时间${\tau _i} = (1 \pm 20{\text{%}} ){R_0}\times (1 \pm 20{\text{%}} ){C_0}$,${\tau _i} = R\times C$、${\tau _i}$的取值范围为:$0.64\tau \leqslant {\tau _i} \leqslant 1.44\tau $($\tau = {R_0}\times {C_0}$,$\tau $表示理论整形时间)。${\tau _i}$以及$n$的变化对读出电路输出脉冲峰值相位的影响如图4(a)所示,固定$n$等于4,${\tau _i}$变化导致读出电路输出准高斯脉冲峰值相移的图像如图4(b)所示。

    图  4  整形时间$ {\tau }_{i} $和整形阶数n对读出电路输出准高斯脉冲峰值相位的影响
    Figure  4.  Effects of shaping time ${\tau _i}$ and shaping order n on peak phase of quasi-Gaussian pulse output in readout circuit

    图4展示了整形时间${\tau _i}$和整形阶数n对准高斯脉冲峰值相位差的影响即$\Delta \varphi = \Delta {t_p}$。如图4(a)所示,相位差$\Delta \varphi $和${\tau _i}$以及$n$成正比,固定$n$等于4后,整形时间常数${\tau _i}$的差值$\Delta \tau $引起的峰值电压相位差$\Delta \varphi $为

    $$ \Delta \varphi = n\Delta \tau $$ (10)

    图4(b)所示:${\tau _i}$值越大,准高斯脉冲的脉冲宽度也越宽。由整形时间常数的偏差$\Delta \tau $引起的准高斯脉冲相位差$\Delta \varphi $在低计数率的时候不会对系统的性能产生影响.计数率较高时,脉冲宽度更宽的读出电路会先发生脉冲堆积,此时可以使用脉冲堆积拒绝电路拒绝堆积的脉冲[13-14],但是会降低系统的计数率,偏离系统的预定设计指标。

    为测试读出电路增益不相等对探测器成像的影响,利用MATLAB模拟了探测器3路读出电路增益$m > k > h$等6种不同增益条件下,不同位置的电子云质心坐标解码图像如图5所示。图5中黑色圆圈组成的图像表示忽略读出电路增益不相等;使用公式(7)解码不同位置的电子云质心坐标得到的图像(畸变图像),黑色点组成的图像表示考虑读出电路增益偏差;使用公式(8)解码不同位置的电子云质心坐标得到的图像(无畸变图像)。模拟解码时,利用图5各个子图上方黑体标明的增益条件,结合上述积分电容读出电路增益的关系,3路读出电路的积分电容分别在0.8 pF、1 pF、1.2 pF中对应选择、整形阶数$n$等于4、整形时间${\tau _i}$等于10 μs。

    图  5  3路读出电路增益互不相等造成的探测器图像畸变
    Figure  5.  Detector image distortion caused by unequal gain of three readout circuits

    图5所示,不同的读出电路增益导致的探测器畸变图像各不相同。为解析图5各个子图不同畸变图像的形成机制,重新整理公式(7)和公式(8)分别如(11)式和(12)式所示。

    $$ {p_{\rm{e}}}(x,y) = \left\{ \begin{gathered} x = \dfrac{{2{Q_{\rm{s}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} = \dfrac{{2{V_{\rm{s}}}}}{{{V_{\rm{w}}} + {V_{\rm{s}}} + {V_{\rm{z}}}}} = \dfrac{2}{{\dfrac{{{V_{\rm{w}}}}}{{{V_{\rm{s}}}}} + 1 + \dfrac{{{V_{\rm{z}}}}}{{{V_{\rm{s}}}}}}} \\ y = \dfrac{{2{Q_{\rm{w}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} = \dfrac{{2{V_{\rm{w}}}}}{{{V_{\rm{w}}} + {V_{\rm{s}}} + {V_{\rm{z}}}}} = \dfrac{2}{{\dfrac{{{V_{\rm{s}}}}}{{{V_{\rm{w}}}}} + 1 + \dfrac{{{V_{\rm{z}}}}}{{{V_{\rm{w}}}}}}} \\ \end{gathered} \right. $$ (11)
    $$ {p'}_{\rm e}(x,y) = \left\{ \begin{gathered} {x'} = \dfrac{{2{Q_{\rm{s}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} = \dfrac{{2\dfrac{{{V_{\rm{s}}}}}{m}}}{{\dfrac{{{V_{\rm{w}}}}}{h} + \dfrac{{{V_{\rm{s}}}}}{m} + \dfrac{{{V_{\rm{z}}}}}{k}}} = \dfrac{2}{{\alpha \dfrac{{{V_{\rm{w}}}}}{{{V_{\rm{s}}}}} + 1 + \beta \dfrac{{{V_{\rm{z}}}}}{{{V_{\rm{s}}}}}}} \\ {y'} = \dfrac{{2{Q_{\rm{w}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} = \dfrac{{2\dfrac{{{V_{\rm{w}}}}}{h}}}{{\dfrac{{{V_{\rm{w}}}}}{h} + \dfrac{{{V_{\rm{s}}}}}{m} + \dfrac{{{V_{\rm{z}}}}}{k}}} = \dfrac{2}{{\gamma \dfrac{{{V_{\rm{s}}}}}{{{V_{\rm{w}}}}} + 1 + \upsilon \dfrac{{{V_{\rm{z}}}}}{{{V_{\rm{w}}}}}}} \\ \end{gathered} \right. $$ (12)

    (12)式中的$\alpha = {m / h}$, $\beta = {m / k}$, $\gamma = {h / {m = {1 / \alpha }}}$, $ \upsilon = {h / {k = {1 / \beta }}} $,$h$、$m$、$k$都是大于0的正数。以图5(a)为例(图5其他子图的图像畸变的形成机制分析方法相同,限于篇幅,在此不再赘述),在$m > k > h$条件下可得到$\alpha > \beta > 1$, $0 < \gamma < \upsilon < 1$。针对同一电子云质心坐标${p_e}(x,y)$,使用解码公式(11)得到的黑色圆圈图像的坐标${p_e}(x,y)$和使用解码公式(12)得到的黑色点图像的坐标${p'}_e(x,y)$在$x$和$y$方向上的差值$\Delta x$和$\Delta y$分别为

    $$ \Delta x = x - {x'} = \dfrac{2}{{\dfrac{{{V_{\rm{w}}}}}{{{V_{\rm{s}}}}} + 1 + \dfrac{{{V_{\rm{z}}}}}{{{V_{\rm{s}}}}}}} - \dfrac{2}{{\alpha \dfrac{{{V_{\rm{w}}}}}{{{V_{\rm{s}}}}} + 1 + \beta \dfrac{{{V_{\rm{z}}}}}{{{V_{\rm{s}}}}}}} > 0 $$ (13)
    $$ \Delta y = y - {y'} = \dfrac{2}{{\dfrac{{{V_{\rm{s}}}}}{{{V_{\rm{w}}}}} + 1 + \dfrac{{{V_{\rm{z}}}}}{{{V_{\rm{w}}}}}}} - \dfrac{2}{{\gamma \dfrac{{{V_{\rm{s}}}}}{{{V_{\rm{w}}}}} + 1 + \upsilon \dfrac{{{V_{\rm{z}}}}}{{{V_{\rm{w}}}}}}} < 0 $$ (14)

    (13)式的$\Delta x$大于0,表示对于同一电子云质心坐标,畸变图像的$x$坐标值大于无畸变时的坐标值${x'}$。和无畸变的黑色点组成的图像相比,黑色圆圈组成的图像在水平方向被拉伸 (如水平方向的黑色箭头所示),且$\Delta x$随着${V_{\rm{s}}}$增大($x$轴正方向)而增大,即随着$x$增大,图像在水平方向上的被拉伸量也增大(水平方向上的黑色箭头更长);同理,(14)式的$\Delta y$小于0,表示畸变图像的$y$坐标小于无畸变图像的${y'}$坐标,和无畸变的黑色点组成的图像相比,黑色圆圈组成的图像在竖直方向上的$y$坐标被压缩 (如竖直方向上的黑色箭头所示),且$\Delta y$随着${V_{\rm{w}}}$增大($y$轴正方向)而增大,表示随着$y$增大,畸变图像竖直方向上受到的压缩量也会变大(竖直方向上的箭头也越长)。无畸变图像在水平方向上受到拉伸和在竖直方向上受到压缩的共同作用下,最终形成了图5(a)所示的黑色圆圈组成的畸变图像。

    除了上述读出电路增益完全不相等外,另一种可能的情况为3路读出电路中的2路增益较为接近,另外一路读出电路增益不等于其他2路,即$h = k > m$等6种情况,其对应的图像畸变如图6所示。如图6的各个子图所示,读出电路不同增益引起的探测器图像畸变各不相同。不同畸变图像的形成机制如图6各个子图中的黑色箭头所示,图像畸变具体的分析方法同上所述。模拟解码时使用的参数如下,针对图6中各个子图上方黑体标明的增益条件,结合读出电路增益和积分电容的关系,读出电路的积分电容${C_f}$分别在1.2 pF、0.8 pF两种电容值中对应选择,整形阶数$n$等于4,整形时间${\tau _i}$等于10 μs。图6的各子图中由黑色圆圈组成的图像表示忽略读出电路增益不相等,使用公式(11)解码不同位置的电子云质心坐标得到的图像(畸变图像);黑色点组成的图像表示考虑读出电路增益偏差,使用公式(12)解码不同位置的电子云质心坐标得到的图像(无畸变图像)。

    图  6  读出单路2路增益相等,第3路增益不等于其他2路导致的图像畸变
    Figure  6.  Detector image distortion under conditions of gain of two channels of readout circuit is equal, and gain of the third channel is not equal to the other two channels

    如上文所述,探测器的读出电路增益不相等会导致探测图像产生畸变。为提高探测器成像质量,需要对读出电路进行增益补偿,矫正探测器的图像畸变。

    为减小探测器读出电路增益不相等产生的畸变,本文提出了一种读出电路增益校准补偿方法来矫正探测器图像畸变,其原理为:将等量电荷输入探测器的3路读出电路,比较3路读出电路输出准高斯脉冲的峰值,以峰值最大的读出电路的增益作为基准,其他2路读出电路峰值电压分别乘以基准电压峰值与各自电压峰值之比作为增益补偿因子补偿读出电路增益。

    探测器的功能相当于一个对外输出电荷的电荷源,可以利用小功率数控电压源和电容作为“电荷源”模拟探测器功能[15],其原理图如图7所示。

    图  7  电荷注入原理图
    Figure  7.  Schematic diagram of charge injection

    图7所示,“电荷源”由小功率数控电压源、电容${C_1}$组成。通过选择电容${C_1}$、${C_f}$以及前置放大器的放大倍数$A$,确保$(1 + A){C_f} > > {C_1}$成立时,阶跃电压源注入到电容${C_1}$中的电荷$\Delta Q$可完全转移到前置放大器中。利用阶跃电压和测试电容${C_1}$对电荷灵敏前置放大器注入的电荷量$\Delta Q$为

    $$ \Delta Q = \frac{{{C_1}(1 + A){C_f}}}{{{C_1} + (1 + A){C_f}}}V \approx {C_1}V $$ (15)

    根据(15)式,利用电压源给电容${C_1}$充电相同的时间后断开,积累在电容${C_1}$中的电荷会转移到前置放大器中。利用“电荷源”模拟探测器输出的电流脉冲,可用于WSZ阳极读出电路增益校准,读出电路增益校准方案如图8所示。

    图  8  探测器读出电路增益校准方案
    Figure  8.  Schematic diagram of gain calibration scheme of detector readout circuit

    图8所示,本文提出的增益校准方案由“电荷源”、控制逻辑(FPGA或者DSP产生用于控制图8中不同电路模块的通断)、单刀多掷开关、二选一多路选择器以及读出电路组成。校准读出电路增益前需要先关闭加载在探测器两端的直流高压电源,通过控制逻辑断开读出电路和探测器阳极电极的连接。以校准W电极读出电路的增益为例(S、Z电极校准方法和W电极相同)。通过控制逻辑控制单刀多掷开关和二选一多路选择器WMUX2_1,连通“电荷源”和W电极的读出电路;利用控制逻辑控制“电荷生成器”产生电荷注入到读出电路,读出电路将注入的电荷转换为准高斯脉冲后被探测器采集系统的ADC采集,探测器的控制单元通过寻峰算法找到被采集脉冲的峰值电压${V_{\rm pw}}$,至此完成了单次读出电路峰值采集。考虑到噪声对准高斯脉冲的干扰,重复对同一读出电路注入等量的电荷$N$次后得到$N$个脉冲峰值,单次采集的脉冲峰值记为${V_{{\rm pw} \_i}}$。利用(16)式对$N$个脉冲峰值求和取均值后得到W电极的脉冲峰值均值${V_w}$:

    $$ {V_{\rm{w}}} = \frac{{\displaystyle\sum\nolimits_{i = 1}^N {{V_{{\rm pw}\_i}}} }}{N} $$ (16)

    S电极和Z电极的校准方法和W电极相同,校准S电极和Z电极的读出电路后分别得到各自读出电路电压峰值为${V_{\rm{s}}}$和${V_{\rm{z}}}$,比较电压${V_{\rm{w}}}$、${V_{\rm{s}}}$、$ {V_{\rm{z}}} $的大小,假设3路读出电路输出的脉冲峰值最大值为${V_{\rm{w}}}$。令$r = {{{V_{\rm{w}}}} / {{V_{\rm{s}}}}}$, $ l = {{{V_{\rm{w}}}} / {{V_{\rm{z}}}}} $,$r$和$l$分别为S电极以及Z电极读出电路待求的增益补偿因子。将$r$和$l$带入公式(7),得到新的电子云质心坐标解码公式为

    $$ {p_e}(x,y)\left\{ \begin{gathered} x = \frac{{2{Q_{\rm{s}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} = \frac{{2r{V_{\rm{s}}}}}{{{V_{\rm{w}}} + {V_{\rm{s}}} + l{V_{\rm{z}}}}} \\ y = \frac{{2{Q_{\rm{w}}}}}{{{Q_{\rm{w}}} + {Q_{\rm{s}}} + {Q_{\rm{z}}}}} = \frac{{2{V_{\rm{w}}}}}{{{V_{\rm{w}}} + r{V_{\rm{s}}} + l{V_{\rm{z}}}}} \\ \end{gathered} \right. $$ (17)

    需要注意的是,公式(17)是电压${V_{\rm{w}}}$大于${V_{\rm{s}}}$和${V_{\rm{z}}}$的一种特殊形式,如果比较结果为${V_{\rm{s}}}$或${V_{\rm{z}}}$最大,则需要补偿的电压会发生改变,公式(17)也会略微改变。

    利用公式(17)对电子云质心坐标解码时和读出电路增益因子$h$、$m$、$k$的值无关,读出电路增益不相等时,利用上述增益补偿方案校准读出电路增益,得到增益补偿因子$r$和$l$后即可以对电子云质心坐标解码。

    为验证上述增益补偿方案,随机选择W、S、Z三个电极读出电路前置放大器的积分电容${C_f}$分别为0.89 pF、1.17 pF、0.98 pF、固定整形阶数$n$等于4,整形时间${\tau _i}$等于10 μs。利用公式(17)对不同位置的电子云质心坐标解码得到图像如图9(a)所示。作为对比,图9(b)给出了读出电路增益不相等,且未经增益补偿,利用公式(7)解码不同位置电子云质心的图像。

    图  9  读出电路经过增益校准补偿后的不同位置电子云质心坐标解码图
    Figure  9.  Decoded diagram of electron cloud centroid coordinates at different locations after gain calibration compensation of readout circuit

    图9(a)所示,黑色点组成的图像表示读出电路增益相等时不同位置电子云质心坐标的解码图像;黑色圆圈组成的图像为读出电路增益不相等,经过增益补偿后不同位置电子云质心坐标的解码图像。对比图9(a)的黑色圆圈组成的图像和黑色点组成的图像可知,读出电路增益不相等,经过增益补偿校准后的解码图像非常接近于读出电路增益相等时的解码图像。图9(b)为读出电路增益不相等且未经增益校准补偿的图像,其中黑色点组成的图像为读出电路增益相等时的解码图像,黑色圆圈组成的图像为读出电路增益不相等且未经增益补偿解码图像。对比图9(a)、图9(b)的黑色圆圈组成的图像可知,未经增益补偿的解码图像(图9(b)黑色圆圈组成的图像)畸变更为严重,经过增益补偿的解码图像(图9(a)黑色圆圈组成的图像)非常接近读出电路增益相等的无畸变图像。读出电路增益相等时的解码图像(图9(a)黑色点组成的图像)和经过增益补偿的图像(图9(a)黑色圆圈组成的图像)的偏差主要限制于读出电路信噪比的大小、ADC的采样精度、补偿校准次数N的大小以及增益补偿因子$r$和$l$取值的精度。

    为了验证本文提出的探测器3路读出电路增益不相等会导致探测器成像产生畸变,选取了15×15的针孔掩膜板测试探测器读出电路增益不相等引起探测器成像产生畸变,针孔掩膜板每个小孔的直径为300 μm(孔直径误差为±20 μm),孔间距为1 700 μm(孔间距误差为±40 μm),测试时将针孔掩膜板紧贴于探测器入射窗口的表面作为探测器采集的目标图像。开展畸变测试时,探测器的3路前置放大电路的3个积分电容${C_{f\_{\rm{w}}}}$、${C_{f\_{\rm{s}}}}$以及${C_{f\_{\rm{z}}}}$的电容实测值分别为0.93 pF、1.02 pF、1.08 pF($h > m > k$)。测试的畸变图像如图10所示。图10以及图11中小孔所成的像明暗不均衡是可能是由于小孔的孔径误差过大,单位时间内透过小孔的光子数量不一致最终造成小孔成的像明暗不一致。

    图  10  探测器的读出电路的增益不相等时的畸变图像
    Figure  10.  Distorted image when gain of detector readout circuit is not equal
    图  11  增益补偿后的正常图像
    Figure  11.  Normal image after gain compensation

    对比图10图9(b)可知,实验测试图像和仿真的预期图像基本相符,探测器成像的畸变成因如图9(b)的黑色箭头所示,图像在水平方向和竖直方向均被拉伸,最终在水平方向拉伸和竖直方向“拉伸”的共同作用下,形成了图10所示的畸变图像。数据采集系统将探测器采集的图像量化为512×512像素的图像,图10左上方和左下方的部分图像的像素的解码坐标超出了预先量化好的坐标的范围,因此针孔掩膜板左半部分小孔在探测器上成的图像由于畸变被“挤出”预定坐标,最终导致图10左边部分像素“缺失”。实验测试结果验证了本章提出的探测器读出电路增益不相等会导致探测器成像畸变的正确性。

    为了校正探测器的畸变图像,使用本文提出的增益补偿方法后得到的正常图像如图11所示。

    图11所示,经过增益校准后的探测器图像和图10所示的畸变图像相比,图10所示的图像在水平方向和竖直方向被“拉伸”的情况得到大幅度改善,有效降低了图10所示的探测器成像畸变,探测器的成像质量有所提高,最终得到图11所示的正常图像。实验结果验证了本章提出的增益校准补偿方案的有效性。

    仿真结果表明MCP电子倍增WSA电荷分割型阳极成像探测器读出电路增益不相等会导致探测器成像产生畸变,探测器的3路读出电路面对不同的增益条件,探测器图像产生的畸变也各不相同。仿真和实验结果验证了本文提出的读出电路增益补偿方案能有效减小由于读出电路增益不相等导致的探测器图像畸变。实验结果验证了使用本文提出的读出电路增益补偿方案,可有效降低读出电路对前置放大器积分电容精度的要求。

  • 图  1   WSA阳极不同电极结构示意图

    Figure  1.   Schematic diagram of different electrode structures in WSA anode

    图  2   读出电路信号传输模型和波形示意图

    Figure  2.   Signal transmission model of readout circuit and waveform diagram

    图  3   积分电容和整形阶数对输出准高斯脉冲峰值的影响

    Figure  3.   Effects of integral capacitance and shaping orders on quasi-Gaussian pulse peak

    图  4   整形时间$ {\tau }_{i} $和整形阶数n对读出电路输出准高斯脉冲峰值相位的影响

    Figure  4.   Effects of shaping time ${\tau _i}$ and shaping order n on peak phase of quasi-Gaussian pulse output in readout circuit

    图  5   3路读出电路增益互不相等造成的探测器图像畸变

    Figure  5.   Detector image distortion caused by unequal gain of three readout circuits

    图  6   读出单路2路增益相等,第3路增益不等于其他2路导致的图像畸变

    Figure  6.   Detector image distortion under conditions of gain of two channels of readout circuit is equal, and gain of the third channel is not equal to the other two channels

    图  7   电荷注入原理图

    Figure  7.   Schematic diagram of charge injection

    图  8   探测器读出电路增益校准方案

    Figure  8.   Schematic diagram of gain calibration scheme of detector readout circuit

    图  9   读出电路经过增益校准补偿后的不同位置电子云质心坐标解码图

    Figure  9.   Decoded diagram of electron cloud centroid coordinates at different locations after gain calibration compensation of readout circuit

    图  10   探测器的读出电路的增益不相等时的畸变图像

    Figure  10.   Distorted image when gain of detector readout circuit is not equal

    图  11   增益补偿后的正常图像

    Figure  11.   Normal image after gain compensation

  • [1]

    SCOTT A B, KENNETH F D, STEFAN E T, et al. On-orbit characterization and performance of the HIRAAS instruments aboard ARGOS: LORAAS sensor performance[J]. SPIE,2002,4485:349-360.

    [2]

    MCCOY P R, MEIER R R, WOLFRAM D K, et al. Instrumentation on the RAIDS experiment. I-Far ultraviolet imaging spectrograph and scanning grating spectrometers for the middle and near ultraviolet[J]. SPIE,1992,1745:75-88.

    [3]

    GLODSTEN O J, HUMM C D, PAXTON J L, et al. Performance of the wedge-and-strip microchannel plate detectors and electronics for the Global Ultraviolet Imager[J]. SPIE,1999,3765:408-416. doi: 10.1117/12.366521

    [4]

    KENNETH F D, NICHOLAS A C, SCOTT A B, et al. The special sensor ultraviolet limb imager instruments[J]. Journal of Geophysical Research:Space Physics,2017,122:2674-2685. doi: 10.1002/2016JA022763

    [5]

    ZHANG Xiaoxin, CHEN Bao, HE Fei, et al. Wild-field auroral imager onboard the Fengyun satellite[J]. Light:Science and Applications,2019,8(1):47-58. doi: 10.1038/s41377-019-0157-7

    [6]

    CHAUDHARI A J, JOSHI A A, WU Yibao, et al. Spatial distortion correction and crystal identification for MRI-compatible position-sensitive avalanche photodiode-based PET scanners[J]. IEEE Transaction on nuclear science, 2009, 56(3):549-556.

    [7] 缪振华. 基于楔条形阳极探测器的单光子成像系统[D]. 西安: 西安光学精密机械研究所, 2008: 148.

    MIU Zhenhua. A single photon imaging system based on wedge and strip anodes detector[D]. Xi'an: Xi'an Institute of Optics and Precision Mechanics, 2008: 148.

    [8]

    TREMSIN A S, VALLERGA J V, SIEGMUND O H W. Overview of spatial and timing resolution of event counting detectors with microchannel plates[J]. Nuclear Instruments and Methods in Physics Research A,2020,949:1-11.

    [9] 石吟馨, 戴丽英, 唐光华. 楔条形阳极紫外探测器综述[J]. 真空电子技术,2011,6:1-5. doi: 10.3969/j.issn.1002-8935.2011.01.001

    SHI Yinxin, DAI Liying, TANG Guanghua. Overview of UV detectors with wedge-and-strip anodes[J]. Vacuum Electronics,2011,6:1-5. doi: 10.3969/j.issn.1002-8935.2011.01.001

    [10]

    XU Neng, LIU Yangan, SHENG Lizhi, et, al. Novel electronic readout system for mcro-channel plate detector with wedge and strip anode[J]. Journal of Modern Optics,2019,66(6):697-701. doi: 10.1080/09500340.2019.1567844

    [11]

    HATFIELD J V, BELL S, NEAVES P I. A wedge and strip particle detector based on a current-mode approach[C]. 38th Midwest Symposium on Circuits and Systems. Proceedings, Rio de Janeiro. Brazil: IEEE, 1995: 193-196.

    [12]

    OGLETREE D Frank, BLANKMAN G C, HWANG R Q, et al. A new pulse counting low‐energy electron diffraction system based on a position sensitive detector[J]. Review of Scientific Instruments,1992,63(1):104-113. doi: 10.1063/1.1143196

    [13]

    ROZENSTEIN S. Pile up rejection cicuits[J]. Nuclear Instruments and Methods in Physics Research A,1961,11:316-320. doi: 10.1016/0029-554X(61)90031-3

    [14]

    GLENN Andrew, CHENG Qi, KAPLAN ALAN D, et al. Pulse pileup rejection using two-component Gaussian mixture model for the fast neutron detection with pulse shape discriminating scintillator[J]. Nuclear Instruments and Methods in Physics Research A,2021,998:1-8.

    [15]

    GRYBOS P. Front-end electronics for multichannel semiconductor detector systems[J]. Hep Websearch Hep,2012(8):55.

  • 期刊类型引用(1)

    1. 翟一行,王斌. 基于交叉条形阳极探测器混叠数据的目标定位方法. 激光与光电子学进展. 2024(18): 312-318 . 百度学术

    其他类型引用(0)

图(11)
计量
  • 文章访问数:  189
  • HTML全文浏览量:  60
  • PDF下载量:  34
  • 被引次数: 1
出版历程
  • 收稿日期:  2022-08-29
  • 修回日期:  2022-10-12
  • 网络出版日期:  2022-11-13
  • 刊出日期:  2022-11-14

目录

/

返回文章
返回