硅光电二极管伽马射线探测器
李宣成
黑体辐射2023/05/17原创 核技术与核仪器仪器仪表 IP:湖北
中文摘要
研究与设计基于硅光电二极管的伽马射线探测器,对硅光电二极管的伽马探测器进行了简单建模,制作了样机实物,并对样机进行简单测试。
关键词
光电二极管伽马射线探测器

该研究与设计为本人自主完成,现将内容发布到科创论坛为科研爱者提供参考。任何人均可使用(本人能力有限😂,有可能会出错,还请大佬指正🙏)。


硅光电二极管是可用于探测伽马射线的探测器件中最便宜的一种,理论上可以实现低分辨率的能谱探测,但因其原本设计并不是用来探测伽马射线,在实际应用中需要特殊的偏压结构与灵敏的信号提取电路才能从噪声中提取出核脉冲信号。本文将使用硅光电二极管实现伽马射线的探测,并提取核脉冲信号。

1 .探测电路框架结构

基于硅光电二极管的伽马射线探测电路主要包括反偏电路和信号提取与处理电路两个主要模块,如图1所示。伽马射线传感部分主要有两种方案:1、使用闪烁体将射线转换成可见光,再利用硅光电二极管将微光转换成电信号。2、直接探测,让射线直接进入硅光电二极管耗尽层产生电信号。信号提取电路输入端均使用JFET运放进行提取处理,以获得有用的核脉冲信号。

                  image.png         

图1. 硅光电二极管探测设计方案

方案一设计较为复杂,设计架构可拓展性更强,其反偏电压可在2.5V-30V之间调节,通过修改电路参数可以适配硅光电倍增管(SiPM)器件。方案二在模型仿真的指导下通过器件选型优化可获得较好的输出脉冲波形,在经过大量的实验后伽马射线探测器成本可降至60元以内。

2.硅光电二极管电路建模

硅光电二极管反偏模式如图2所示,其中二极管以反向偏置接入电路,由电源VCC为其提供反偏电压;电阻RA为限流电阻防止二极管击穿,同时防止二极管产生的信号进入电源发生衰减;C1为隔直电容,用于提取二极管上的核脉冲信号;RL为输出负载电阻。

image.png

2 硅光电二极管反偏模式

在反偏模式下光电二极管内部简化模型如图3所示,其中I2为光信号产生电流;Rs为光电二极管的反偏等效电阻,该值与探测器中耗尽层自激发相关,温度越低自激发产生的载流子越少,其阻值越大;Ct为探测器的极间电容[1],该值与两级间耗尽层厚度以及平板面积成正比,耗尽层越厚、平面面积越小极间电容也越小。图中省去了内部理想二极管,因为在反偏下理想二极管不会有任何电流,且会增加计算负担。

image.png

3 反偏模式下的光电二极管内部简化模型

光电二极管的反偏电流(也称为暗电流)I1的大小与反偏电压和工作温度相关,温度越高、反偏电压越大,I1的值也越大[2]。以S2387为例,暗电流与反偏电压的关系如图4所示。

image.png

4 S2387反偏电压与暗电流关系图

由图4可见,暗电流随反偏电压的增长是非线性的。采用恒定的反偏电压,在25摄氏度室温参数下搭建具体模型如如图5所示。其中,RsC1I1为实际测量值。测试条件为:将RA设为1MΩRL以示波器输入阻抗代替,C1取值为1uF

image.png

5 反偏模式下的硅光电二极管模型

在测量时,需要使用铝箔对光电二极管做避光处理,通过多次测量后取平均值的方法得到实际测量值。

二极管两端的反偏电压测量需要利用示波器的触发模式采集,在示波器探头接触光电二极管两端的瞬间其电压值就会下降。这个现象的原因是示波器探头的输入阻抗相比光电二极管要小得多,通常示波器的输入阻抗为1MΩ-10MΩ,但探测器的内阻通常为GΩ单位 [3]

在反向偏压下的结电容Ct可通过测量C1开路端对地的电容得出。因为C1的容值远大于Ct的容值,此时测试的电容量可以近似的看做Ct的容值[4]

在估算时可以对暗电流与内阻进行简单化简。在静态模式下考虑温度等变量时,暗电流恒流源会发生变化,而内电阻会随信号传输过程中分压的变化而变化。在实际推导过程中发现,结电容Ct的大小对输出信号幅值影响相当大,而内电阻与恒流源对输出信号幅值的影响相对较小。在温度和偏置电压均恒定不变时,可将暗电流现象用内电阻电流简化处理,以降低建模难度。简化后的电路如下图6所示。

image.png

 

6 反偏模式下的光电二极管简化模型

3.闪烁体光脉冲仿真建模

要想求出该模型的信号输出,首先要知道输入光的波形,具体来说是闪烁体产生的光波波形。常见的无机闪烁体光衰中有两种成分,其光强度I(t)随时间衰变规律为[5] 

image.png   (1)

其中,I_f 、I_s分别为快慢光衰的组成比例,入射粒子的能量成正比关系;τ_f 、τ_s分别为快慢光衰的衰减时间量,同种材料中衰减时间量为常数。参考文献[6]其中碘化铯闪烁体光波波形如图7所示。

  image.png

图7 参考文献[6]中碘化铯闪烁体光波波形

图7 的数学表达式结构与式(1)相似。根据该文献中的数据和指数衰减拟合函数,其光衰部分波形N(t)的表达式为:

    image.png        (2)

其中,快光衰时间参数τ_f≈2.48∙10^(-9),慢光衰时间参数τ_s≈1.912∙10^(-8),快光衰比例系数A_f≈-0.397,慢光衰比例系数A_s≈-0.0704,基线修正参数y_0≈0.00047。

式 (3-2)可较为准确的拟合表达闪烁体光衰,但该式中有两个指数函数,导致推导工作量增大。通过分析可以发现,慢光衰成分要比快光衰小得多,可以考虑使用快光衰成分近似表达光衰期间的波形。将慢光衰的成分组合到快光衰中,对快光衰时间参数乘上比例系数1.2,并将基线降低0.04,可得到与原拟合函数波 式(2)形相似的波形N_sf (t):

    image.png          (3)

简化后的光衰拟合函数N_sf (t)与原拟合函数N(t)的差异如图8所示。

image.png

图8 N_sf (t)与N(t)波形的比较

通过测试得出一个与真实波形高度相似的平滑函数:

image.png      (4)

其中,A=5,α=4×10^6,β=4×10^5,B=-8×10^(-8)。

完整的拟合函数波形I_Photon (t)如图9所示,可近似的代表闪烁体输出光波波形。                                           image.png

9 拟合函数与I_Photon (t)参考波形对比

拟合波形所参考的光电倍增管采集波形是负向的,而实际光波波形是的正向的,将拟合的函数进行重新定义,并省略整体波形贡献很小的参数B,得到:

    image.png     (5)

其中,A=5,α=10,β=1。

image.png

10 重新定义的I_Photon (t)函数

式(5)拟合出的闪烁体发出光波波形形状如图10所示,波形外观与原拟合波形相似,但时间与峰值还未进行约束。理论上射线的能量越大,闪烁体发出的光能越多,光波波形基本保持不变。将上述波形进行积分,利用发出光子总能量来约束波形,具体包含两个约束关系:1、波形下降沿的时间,对应光波光衰特性。2、波形积分面积,对应光子输出总能量。

约束一:通过查阅相光资料可以得知光衰平均时间为40nS~80nS[7],这里取中间值60nS衰减到峰值的5%

I_Photon (t)峰值点为(0.1,0.184),波形相对下降到峰值5%时点为(0.567,0.009),下降时长为0.467s,与60ns做比得到时间缩放参数m_1=0.467s/60ns=7.783×10^6。则修正时间函数为:

image.png (6)

约束二:这里假设入射的伽马射线最高能量为300KeV,伽马射线入射碘化铯(铊)闪烁体后产生的光子平均能量为:2.25eV≈3.6×10^(-19) J。

产生的光子数量约为:(( 52000)⁄(1MeV ))×300KeV=15600 (个)。

总光子能量为:15600×2.25eV=35100eV=35.1KeV≈5.616×10^(-15)  J。

假设1mm厚碘化铯(铊)闪烁体产生的光子能量以2.25eV平均值进行计算。设光子能量参数为m_2,则对光波波形m_2×I_Photon (m_1×t)的积分为光子的总能量。则产生的光能为:

image.png

已知积分公式:

image.png    (7)

代入后计算得产生的光子总能量为:

image.png

则m_2=  (5.616×10^(-15) J)/(6.36×10^(-9)J )=8.83×10^(-7)。。

闪烁体在300KeV的伽马射线下产生的能量约束与时间约束的拟合函数为:

image.png            (8)

将式(8)代入数据后得到:

image.png     (9)

在硅光电二极管中描述光强与电流转换的参数为光敏感度S_k,量纲为A/W,是指输出电流变化(A安培)与输入光功率(W瓦)变化的比值[8]。所以为了后续求解过程还需要将光功率波形转换成光电流信号。

通常晶体内光自吸收损耗可以高达40%,在折射率突变的传播界面损耗一般在10%以内[9]。假设光传输效率:η_(Photon-m)=54%。

结合硅光电二极管在520nm波长下的光敏感度S_k,可以将P_Photon (t)转换成光电流:

image.png              (10)

例如S2387-1010520nm波长下的S_k≈0.33 (  A/W  )。则该光信号在S2387-1010内部光电流I_Photoelectron (t)的波形如图11所示。

image.png

11 光电二极管内部光电流波形

当t=12.85ns时,光电电流达到最大值14.81μA。由于该光电流是硅光电二极管内部信号,该信号还需要代入具体电路中计算才能得到真实的输出信号。

I_Photoelectron (t)加入u(t)阶跃信号限制其定义域,代入参数后的普拉斯变换为:

  image.png

因选用的拟合函数在(-∞,+∞)绝对可积,则其变换后的收敛域包含极轴。

image.png  

定义域:(s≥-7.783×10^7)。

其中光电信号的拉式变换简化表达如下:

image.png  (11)


定义域:(s≥-B_Photoelectron),当伽马射线能量为300KeV时各参数为 

A_Photoelectron=1.57×10^(-7)、B_Photoelectron=7.78×10^7、

C_Photoelectron=7.78×10^6。


4.闪烁体光脉冲在硅光电二极管中的响应估算

在第3节通过近似计算得到了1mm厚的碘化铯(铊)搭配S2387-1010光电二极管,单个300KeV伽马射线在其内部产生光电电流信号I_Photoelectron (t)。带入测试电路中求解实际输出信号V_out (t)。测试电路中如图12所示。

     image.png        

12 光电二极管光脉冲估算模型

首先设输入电流信号为I_IN=I_2=I_Photoelectron (t),光电二极管端电压为U_2(即R_s两端电压)。其中,R_S、C_t是在反偏模式下实际测量出来的,R_L为后级JFET运放电路的输入阻抗,R_A 、C_1为外围电路元器件。

外围器件参数与偏置电压为:R_A=1MΩ,C_1=1uF,Vcc=30V。

实际测量得出下列参数:R_S≈374MΩ,C_t≈620pF,R_L≈100MΩ,U_2≈29.92V。

利用拉普拉斯变换后求解Vout:定义区域:(s≥-B_Photoelectron)。

image.png (12)

将Vout(s)简化为零极点模式:

image.png (13)

定义区域:(s≥-B_Photoelectron)。其中:

image.png

代入输入光电信号P_Photon (s)的表达式得:  

image.png (14)

此时利用留数法将上式做拉式逆变换,求解出时域信号。其中,对应的共轭极点求解如下:

image.png

式(14)求解出的各留数对应的信号:

image.png

其中:

image.png

输出信号为各留数所对应的信号:

image.png (15)

最终得到输出电压信号V_OUT (t)的波形如图13所示,可见输出信号的最大值在800nV左右,产生的信号非常小。

image.png

13 V_OUT (t)估算波形

降低结电容C_t可以获得更好的信号输出。如果使用结电容C_t=2pF的MAQ87-CQ硅光电二极管,在其他参数不变的条件下波形输出如图14所示,可见输出信号的峰值在248uV左右,与结电容为620pF的输出相比信号峰值增加了310倍。

image.png

14 将MAQ87-CQ参数代入模型后的输出

 综上所述,选用硅光电二极管器件做闪烁型探测器产生的信号幅度非常小,脉冲提取难度较大;但脉冲宽度较大,使用微伏信号带宽放大器可以实现脉冲信号的提取。


5. 硅光电二极管探测器方案一制作:

方案一利用碘化铯(铊)闪烁体搭配大面积S2387-1010BR硅光电二极管实现伽马射线的探测。该设计所采用的碘化铯(铊)闪烁体(CsI(Ti))具体参数如表3.1所示。对于低能伽马射线,2毫米的CsI(Ti)晶体俘获探测效率接近100%[10] ,对于能量较高伽马射线能射线则需要厚的更厚的晶体才能实现高捕捉效率。

        表1 CsI(Ti)晶体参数[11]

密度(g/cm3)密度(g/cm3)光子产额(ph/MeV)光子产额(ph/MeV)峰值波长(nm)吸收长度(cm)成分
4.511.7852000550100037CsI(Ti)

3mm厚的碘化铯(铊)闪烁体光输出率为56.3%[12]。设计采用3mm厚10*10mm的方形薄片晶体作为实验晶体,如图15(a)所示。晶体的五个面镀有2um厚的银全反射膜,以提高光子输出效率。

image.png

15 闪烁体外形设计图与闪烁体实物图。

硅光电二极管选型:

2 常见的大面积硅光二极管的部分参数

image.png

无论是否有光照,硅光电二极管内部都存在着各类噪声,例如热噪声、散粒噪声、白噪声、暗电流噪声等。这些噪声的功率总和定义为光电二极管内部的等效噪声功率(NEP)[13],在选取时需要NEP参数尽量小。硅光电二极管的反偏电压越大,内部耗尽层厚度增加也就越大,结电容也会下降,有利于提高输出信号的幅值。

为了使闪烁体发出的光尽可能地被硅光电二极管吸收,选用感光面为10*10mm大面积硅光电二极管,另外对于晶体发出的560nm光敏感性要高。根据查询产品资料给出了一些常见的大面积硅光二极管的部分参数如表2所示。

通过参数对比S8650、S3590型号性能最适用于该方案,但设计最终使用S2387-1010BR。S2387-1010BR在560nm波长的光敏感度为0.33A/W,理论上也可以产生可以观测的信号。经过实际测量,S2387-1010BR静态结电容为3nF,在29V反偏电压下实测结电容减小为620pF。

探测器组装步骤如下:

1、如图16所示,通过3D打印的结构将闪烁体与光电二极管固定,在闪烁体晶体与光电二极管添加PMX-200光学硅油以增加光耦合,压紧后采用环氧树脂胶水固定整体结构。

image.png

图3.16 闪烁体与硅光电二极管安装结构

2、利用铅环做电磁屏蔽与辐射屏蔽。为了保证探测器有良好的噪声屏蔽和探测方向性,设计铅环做外围屏蔽。如图17所示,利用翻砂陶铸工艺制作屏蔽铅环。首先将3D打印器件反扣在红砂内压实后取出3D打印器件,在红纱内得到铅环结构的倒模,将融化的铅倒入红砂倒模中即得到该部件。

image.png

17 铅屏蔽环3D模型与实物

3、将闪烁体与光电二极管安装在铅环内环中,如图18所示。将纳米四氧化三铁和环氧树脂按体积比1:2混合,得到不透光且导磁良好的功能性固定胶,将该胶涂抹在探测器的背面,并使用铝箔覆盖,在使用时将铅环接地处理。

image.png

18 探头实物图

探测器电路设计:

探测器电路原理图如图19所示,采用Tl072 JEFT输入运放二级级联方式进行信号放大。其中R1对于光电二极管为输出负载,对应建模中的负载RL。运放使用单节9V电池以单电源方式供电,探测器端偏压电路使用三节12V电池串联供电,利用稳压二极管为探测器提供30V的反向压。

image.png

19 探测器电路原理图设计

探测器实物制作如图20所示。

image.png

(a.探测器电路实物,b.探测器与屏蔽结构组装实物,c.探测器整体实物外观)

20 探测器实物

测试使用232Th(独居石标本)射线激发,激发的波形如图21所示。

image.png

图21 闪烁体-硅光电二极管探测232Th输出脉冲

为了验证该脉冲为232Th释放的伽马射线激发,设计了如下对比实验:

1.232Th放置时可观测到明显脉冲信号,待稳定工作一段时间(20min),期间可观察到大量且明显的脉冲信号;撤出232Th后脉冲消失,说明该脉冲确实为232Th所激发的。

2.为了进一步验证是伽马射线所激发的脉冲,在探测器工作时使用铅板、铝板分别阻挡232Th,实验发现当铝板阻挡时探测脉冲数量有一定的下降,从每平均秒8.62个下降到每秒8.37个;当使用5mm铅板阻挡时,探测脉冲几乎为零,测试结果证明输出脉冲确实为伽马射线激发所产生的。

232Th射线测试结果可以看出,与改进前相比输出信号的幅度明显增加。同时,在测试结果中输出核脉冲的幅值有明显的阶梯性,表明该探测器可能具有一定的能谱分析能力。

为了验证该探测器是否有能谱分析能力,采用人工统计的方法创建能谱信息,具体方法为:设置示波器触发并读取脉冲的峰值,进行人工统计,共统计了309个核脉冲数据如图22(a),实现能谱图形的构建,如图22(b)所示。

image.png

图22 手工统计能谱数据图。

通过实际测量,该方案输出能谱分辨率较差,几乎不能辨认特征峰信息,无法探测能量较低的241Am(烟雾报警器)射线。该方案性能较差的主要原因是选用的光感器件结电容过大,造成输出脉冲信号幅值过小。通常硅光电二极管结电容越小,其感光面积也就越小。由于工艺水平等条件限制,很难在1mm*1mm的感光面上安装闪烁体颗粒,但可以尝试利用PN结耗尽层对射线敏感效应实现探测,使用结电容尽量小、耐压足够高的硅光电二极管,在反偏模式下进行直接探测。


6.硅光电二极管探测器方案二

在方案一的研究过程中发现:结电容参数越小的器件输出信号越大。例如常见的MAQ87-CQ感光面积只有2mm2,无偏压时结电容为2pF,在增加偏压后结电容会进一步的减小。伽马射线在PN结中也会直接激发产生载流子,从而实现伽马射线-电信号转换。

信号放大电路:


方案二电路如下图23所示,相比方案一结构简化许多。在本次实验中使用MAQ87-CQ的暗电流为0.1nA,使用24V电池直接为其提供反偏电压。

image.png

图23 方案二实际应用电路

选用运放OPA656NB作为核心器件,信号放大电路仿真结果如图24所示。通过仿真分析,电路没有发生自激振荡。

image.png

图24 方案二微小信号放大电路输入输出仿真

电路的频率与增益倍数仿真如图25所示。Y轴为log(增益倍数),X轴为频率线性。

image.png

图25 方案二微小信号放大电路频率响应仿真

电路实物如下图26所示。运放供电与反偏电压供电同样使用12V23A)电池供电。电路采用铁盒屏蔽以降低输出噪声,采用大面积双面覆铜板以提供更好的信号地平面,信号地平面与外壳地通过42mm铜螺栓连接。在多种屏蔽抗干扰的优化下,输出信号的到很好的噪声抑制。

image.png         

(a.探测器正面,b.探测器电路板背面)

图26 硅光电二极管伽马射线探测器方案二实物

该探测器不仅可以探测钍-232射线,还可以探测能量较低的镅-241射线,输出波形相比方案一有着很大改善,如下图27所示,(a)为Am-241激发波形、(b)为钍-232激发波形。

image.png

图27硅光电二极管伽马射线探测器方案二输出信号

从图27中可以明显看出,不同能量的伽马射线激发的脉冲形状有着明显的不同,为了更好的提取脉冲信号,还需要对信号进行滤波处理。

在设计滤波器电路前,对示波器采样到的模拟信号进行拟合,如图28所示。

image.png

图28 探测器输出信号的数学拟合

图28的拟合函数主表达式为:

image.png

由于滤波器生成的频率响应需要用硬件实现,所以滤波优化的框架采用30阶布拉克曼窗构建的FIR低通滤波,理论上使用该框架生成的滤波器频率响应可以使用硬件电路近似实现。最终得到滤波后较为纯净的信号如图29所示。

image.png

图29  FIR低通滤波输出

在设计FIR低通滤波时主要考虑了输出波形的纯净度,最后得出有较好输出波形的滤波器与常见滤波器的频谱响应不太一样,其中通带内信号放大了10倍,过渡带比较平摊。如图30(a)所示。之后参考该频谱响应设计硬件滤波器电路,最终实现的硬件滤波器频谱响应如图30(b)所示。

image.png

图30  滤器频谱响应

通过参考图30(a)频谱设计的硬件滤波器有少许偏差,经过多次调节后均未能实现两者频率响应完全相同,但频响大体形状相似,最后构建的实际滤波器电路频率响应如图30(b)所示。滤波器采用多通道反馈结构设计,电路如图31所示。主运放采用OPA656NB,单电源驱动设计。

image.png

图31 硬件滤波器原理图

硬件滤波器在实际电路中搭建如图32所示。

image.png

图32硬件滤波器实物(款内为LPF电路)


信号通过硬件滤波器电路后,输出脉冲波形如图33所示。

image.png

图33 硬件滤波器输出脉冲波形



7. 方案二探测器测试:

7.1使用烟雾报警器进行测试:

image.png

图34 硅光电二极管伽马射线探测器对241Am能谱测试结果

image.png

图35 Theremino_MCA 核素查询结果

硅光电二极管伽马射线探测器在该项测试中有脉冲输出,对应采集的能谱如图34所示。能谱的波形和已知能谱相反,低能区峰值高而高能区峰值低且高能区明显可见两个峰。

这可能是由于光电二极管对能量高的伽马射线捕捉能力差所导致的,分析其工作原理不难看出该探测器敏感区域为半导体材料中的耗尽层,未经过半导体工艺改造的耗尽层通常很薄,能量高的伽马射线可能还未与耗尽层中的物质进行能量交换就穿透过去,造成其对能量高的伽马射线捕捉能力差。对于所捕捉的主峰能量分辨率为29.7%

7.2 低能射线探测效率测试:

为了获得该探测器的射线捕捉效率,这里使用了自制的低能伽马射线能谱仪,该仪器的之后在后期排版整理后会发布到科创论坛上。供大家参考设计,相关的设计源文件也会全部开源到论坛中。

实验具体内容:通过将241Am源分别放置8cm、9cm、12cm、17cm处进行200秒的计数统计,平均每秒的计数结果如下表所示。其中高能伽马射线能谱仪在该测试中无法探测241Am所释放的低能伽马射线。

                探测器在不同位置的计数

探测器

8cm

9cm

12cm

17cm

Si-PIN探测器

0.765

0.655

0.315

0.255

低能伽马射线能谱仪

515.2

448.7

323.7

193.5

低能伽马射线能谱仪的有效探测面积为706mm2,Si-PIN探测器的有效探测面积为4mm2。在该实验中以低能伽马射线能谱仪在8cm处的测试数据可知,球冠面积与整个球面的比值为1:116.77,则可以估算出放射源放射性活度为60159.4Bq。Si-PIN探测器与低能伽马射线能谱仪的4mm2等效探测面积计数做对比,结果如下表所示。

                低能伽马射线测试中4mm2探测面积等效计数率

探测器(4mm2探测面积)

8cm

9cm

12cm

17cm

低能伽马射线能谱仪

0.75

0.593

0.333

0.166

Si-PIN探测器实测值

2.919

2.542

1.834

1.096

从上表中的对比数据可以看出。在探测241Am低能射线时,Si-PIN探测器相比低能伽马射线能谱仪探测效率低约4倍。其中低能伽马射线能谱仪是使用碘化铯(铊)+PMT+特制窗制作,对低能伽马射线(x射线能端)探测效率极高。

7.3 使用独居石进行测试:


 经过多次调节与改进,硅光电二极管伽马射线探测器未能在232Th能谱测试中体现任何峰如图36所示,但可以对232Th产生脉冲响应。

image.png

36 硅光电二极管伽马射线探测器对232Th能谱测试结果

为了验证不是源的问题,使用CH282闪烁探测器对同一独居石样品进行能谱采集,结果如图37所示。

image.png

37 CH282探测器对232Th能谱测试结果

综上该探测器对高能射线的能谱捕捉很糟糕。

7.4 高能射线探测效率测试

探测器对高能伽马射线源不同距离的计数率对比实验,通过将232Th源与分别放置10cm、12cm、15cm、20cm、处进行200秒的计数统计,平均每秒的计数结果如下表所示。

                高能伽马射线测试中探测器在不同位置的计数

探测器

10cm

12cm

15cm

20cm

Si-PIN探测器

0.11

0.085

0.07

0.045

低能伽马射线能谱仪

40.1

37.5

24.3

16.3

CH282-1闪烁探测器

327.8

293.7

204.5

140.5

其中高能伽马射线能谱仪的有效探测面积为1963mm2。在该实验中高能伽马射线能谱仪探测效率最高,以高能伽马射线能谱仪在10cm处球冠面积与整个球面的比值为1:66.98,则可以估算出放射源放射性活度为21957.6Bq。三种探测器4mm2探测面积的计数率如下表所示。

                高能伽马射线测试中4mm2探测面积计数率


探测器(4mm2探测面积)

10cm

12cm

15cm

20cm

Si-PIN探测器

0.11

0.085

0.07

0.045

低能伽马射线能谱仪

0.227

0.212

0.138

0.092

CH282-1闪烁探测器

0.656

0.598

0.417

0.286



从上表中的对比数据可以看出。在探测高能射线时Si-PIN探测器探测效率低的离谱,相比CH282低6.3倍。

总结:

该设计基本上是手焊没有PCB为大家提供,但对该类型的探测器进行了较为详细的建模,建模过程运用了大量的工程简化,由于本人能力有限其中的推导可能会出现不恰当的简化,但可以为大家提供一个数量级的参考,知道在设计该类型的探测器时所处理的信号的是极微小的。

该探测的能谱性能在整体上是低于闪烁体-光电倍增管方案的,在探测效率上也略低于闪烁体方案,特别是探测高能射线时。

该类型探测器适用于探测低能射线,在探测Am241射线时与闪烁体方案相比探测效率低约4倍,但成本却低57倍。

参考文献:

[1] 张文在 ,姚历农 ,綦杰. 光电二极管探测器直接测量种子源剂量场分布[J]. 中国原子能科学研究院年报,2006(00):269.

[2] 蔡霆力 ,张启兴 ,曾怡 ,等. 光电二极管在辐射探测系统中的线性应用研究[J]. 核电子学与探测技术,2010,30(03):405-409.

[3] 韩婷婷. 频率综合器关键模块的研究与设计[D]. 东南大学, 2016.

[4] 贾牧霖 ,曾国强 ,马雄楠. 一款低成本硅PIN光电二极管偏置电路的设计及应用[J]. 现代电子技术,2014,37(13):159-161.

[5] 郭亚军 ,杨春晖. CsI(Tl)闪烁晶体研究进展[J]. 人工晶体学报,2004(01):77-80.

[6] 任国浩 ,宋朝晖 ,张子川 ,等.纯碘化铯(CsI)晶体的发光与光衰减特性研究[J]. 无机材料学报,2017,32(02):169-174.

[7] Grim J Q , Qi L , Ucer K B , et al. High excitation density, quenching, and radiative kinetics in CsI:Tl scintillator[J]. American Physical Society, 2010.

[8] 钟华强 ,骆志平 ,刘森林 ,等.基于PIN光电二极管的个人辐射探测器研制[J]. 原子能科学技术,2020,54(02):360-365.

[9] 林俏露 ,李公平 ,潘小东 ,等. CsI(Tl)光电二极管探测器效率和输出电流的模拟[J]. 甘肃科学学报,2013,25(01):107-111.

[10] 任国浩 ,宋朝晖 ,张子川 ,等.纯碘化铯(CsI)晶体的发光与光衰减特性研究[J]. 无机材料学报,2017,32(02):169-174.

[11] XXXXXXXXXXXXXXXXXXXXXXXXXXXXXX/halide-scintillators/XXXXXXXXml

[12] 王静 ,苏明 ,杨祖华 ,等. CsI(Tl)闪烁体软X光能量响应的模拟研究[J]. 强激光与粒子束,2015,27(12):98-102.

[13] 张催 ,李公平. CsI(Tl)光电二极管探测器模拟研究[J]. 无损探伤,2014,38(06):12-16.




                                               



[修改于 1年0个月前 - 2023/05/18 02:21:11]

+1  学术分    RodTech    2023/05/18 理论部分充足,过程科学,结果尚可
来自:物理高能技术 / 核技术与核仪器仪器与装备 / 仪器仪表动手实践:实验报导严肃内容:预印本
10
 
22
已屏蔽 原因:{{ notice.reason }}已屏蔽
{{notice.noticeContent}}
~~空空如也
黑体辐射 作者
1年0个月前 修改于 1年0个月前 IP:湖北
920682
硅光电二极管伽马探测器方案一的改进尝试

在上文中对该探测器整体进行了简单描述,里面不乏会有些欠缺的地方,该实践在理论上属于微小信号放大,为了优化方案一中噪声问题,又做了一下工作,但是实验结果却十分让人失望。

在您观看本结内容时,首先说明该过程为失败的案例,本次对硅光电二级管伽马射线探测器的超低噪改进没有成功。主要原因如下:

  1. 在没有对电路噪声源分析时就开始设计,造成电路主要噪声源没有被排出。

  2. 超低噪运放的选型失败,ADA4817的输入端虽然为JFET,但与输出端实测有1兆的内部电阻,没能实现ADA4817的开环增益倍数放大,输出信号较小。

=======正文=======

针对方案一的超低噪优化主要分为一下几个方面:

1、降低电路噪声,使用基准电压芯片为探测器提供更加低噪声的反向偏置电压。

2、设计正负5V超低噪电源为电路供电。

3、采用超低噪声高增益JFET输入运放设计信号提取电路。

1.探测器偏压电路

image.png

改进版探测器反偏电路

改进的探测器偏压电路如图1所示。采用3节12V电池供电,稳压器件调整为LM431基准电压芯片。LM431是TL431电压基准芯片的改进版,与TL431相比具有更低的噪声、更低的功耗、更小的温漂。在LM431输出的基础上增加一阶RC滤波,以获得噪声更低的基准电压。同时,为了抑制R13电阻热噪声,RC中的电容取值调整到uF级别,并搭配两个100nF低ESR电容。


2.超低噪声电源设计

    设计采用9V干电池供电,使用LT3471 DC-DC电路产生两路正负7V输出。LT3471电路设置有自动升降压功能,设计输入电压范围为3V~10V,可以充分利用9V干电池电能。电路工作频率为1.2MHz,通过在输出端增加Pi型滤波电路可以将输出纹波减小到uV级。这两路正负7V电源分别送入精密LDO线性稳压器件,产生两路超低噪声正负5V电源为电路供电。

    采用LTspice仿真软件对LT3471 DC-DC电路仿真设计,如图2所示。其中,+7V采用Boost拓扑结构设计,-7V输出采用SEPIC(Single Ended Primary Inductor Converter)拓扑结构设计,每路分别增加有独立的反馈环路。电路设计输出电流为500mA,+7V输出效率大于80%,-7V输出效率大于60%。

  image.png

图2 LT3471 DC-DC电路仿真设计图

LT3471的反馈环路设计的较滞后,因负载为两个运放,电流波动非常小,以牺牲快速响应获得尽量减小输出纹波。上电波形仿真如图3所示,+7V电源在上电1.2mS后开始稳定,-7V因采用SEPIC结构输出噪声偏大,响应速度比+7V更慢,上电后1.8mS开始进入稳定状态。

image.png

图3 LT3471电路上电波形仿真

    对电路输出噪声进行仿真测试,结果如图5所示。图中黄色为+7V输出噪声,粉色为-7V输出噪声。其中,图4-a为LT3471未加Pi型滤波的输出噪声的频谱,图4-b为LT3471加Pi型滤波的输出噪声的频谱,可见输出噪声中的高频分量显著降低。经过优化后的输出噪声中最高频率分量为-80dB。

image.png

图4 LT3471 DC-DC输出噪声分析

将正负7伏电源输入两颗LDO进一步降低噪声。+7V送入TPS7A4901DGNR 超低噪LDO进一步降噪得到超低噪声+5V输出,-7V送入TPS7A3401DGNR 超低噪LDO进一步降噪得到超低噪声-5V输出。完整的超低噪声正负电源原理图设计如图5所示。

image.png

图5 超低噪电源设计电路

在实际设计PCB时器件布局、器件选型以及线路规划等都需要十分注意。器件布局中要把功率器件规划到一起,开关线路尽量短,避免1.2M的开关信号从PCB线路中发射出去,可以使用滤波器件尽量包围功率器件。器件选型对于噪声控制至关重要,特别是功率电感应尽量选择漏磁极小磁损耗低的器件。线路规划时尽量将输出的正负5V电源线走在PCB的内层,并在传输区域保持完整的地平面,电源线每隔5mm增加一颗100nF的电容接地以保持超低噪特性。


3.微小信号放大电路

微小信号放大电路仿真电路如图3.28所示,探测器输出的微弱信号由初级超低噪运放ADA4817提取。其低频带宽噪声为 2.1nV/√Hz 、 1.4pA/√Hz,粉红噪声(1/f噪声)为7nV/√Hz 、 13pA/√Hz。防止失调电压被放大输出,ADA4817采用C3电容耦合输出。电路中增加C1电容极点,防止ADA4817运放产生震荡。

提取后的信号送入ADA4851-4WYRUZ-R7进行二次放大,ADA4851-4WYRUZ-R7运放芯片中有四个GBW=80Mhz的JFET运放。每级放大倍数为2,四级个同相比例运算放大电路级联,将信号放大16倍输出。仿真电路如图6所示。

image.png

6 微小信号放大仿真电路

仿真结果如图7所示,可以看出ADA4817有没有发生自激振荡。在仿真中没有添加噪声源,通过分析可知电路中主要噪声源为100M电阻的热噪声为1.26uV/√Hz,代入仿真增益中可知输出噪声为126mV/√Hz在四级联放大过程中,C4、R12等构成高通滤波可以对低频噪声滤波,进一步提高输出信噪比。

image.png

图7 微小信号放大电路仿真输入输出

微小信号放大电路PCB设计原理图如图8所示。在仿真电路的基础上增加了很多电容以降低噪声,在每路运放的反馈电阻上并联2pF负反馈电容以提高电路稳定性。

image.png

图8 微小信号放大电路原理图

超低噪电源电路、微小信号放大电路、改进版探测器反偏电路各模块线路连接如图9所示,其中9V、36V、GND为全局通用顶层中部分省略。

image.png

图9 各模块线路连接图

PCB电路设计如图10所示,采用4层板设计,电源与信号线均包埋在内层。PCB增加了两个屏蔽罩,分别用来屏蔽功率器件区域和运放电路区域。

image.png

(XXXXp_layer层设计图,XXXXyer_1层设计图,XXXXyer_2层设计图,XXXXttom layer层设计图,e.元器件布局图,f.实物图。)

图10 PCB设计与实物

探测器通过3D打印结构件与PCB固定,入射面采用0.4mm厚的石墨片作为“窗材料”,如图11所示。石墨具有导电、遮光的特性,不同于与金属箔,石墨内为碳元素,相对原子序数低有利于射线穿透。

image.png

图11 探测器石墨窗安装

安装完成后的整体实物如图12所示,测试时需要将其置于屏蔽铁管内实现稳定工作。

image.png

图12 闪烁体-硅光电二极管伽马射线探测器实物

测试结果:

钍232  (77).png

输出变得更差了。与原方案对比:信号放大倍数减小了,信噪比变得更低了。

原因分析:

  1. ADA4817的内部存在消除Vos的1兆欧电阻,使运放放大无法达到开环增益倍数。

  2. 电路主体的采样架构还是没有改变,所以以前有的固有噪声现在还存在,只是换了噪声更低的器件器件,做了更多的噪声更低的PCB而已,并没有理性思考、运用超低噪设计的精髓。


工程源文件:

1、PCB文件:

attachment icon PCB文件.zip 13.48MB ZIP 26次下载

2、电源仿真文件:

attachment icon LT3471-500mA.zip 70.02MB ZIP 21次下载

3、

引用
评论
5
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
虎哥
1年0个月前 修改于 1年0个月前 IP:四川
920683

强烈赞扬,罕见的、系统的研发记录。在761设计中发现PIN管加偏压之后噪音剧烈增加,且噪音的增加比低能脉冲幅度的增加更快,会导致截止能量抬高,241Am都能被掩盖,于是最后只好不加偏压了。但信号实在太小,为了赶快收场,用了一支静电运放max4477,此物似乎没有同级别的国产替代,留下了软肋。。。

引用
评论
2
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
黑体辐射作者
1年0个月前 IP:湖北
920684
引用虎哥发表于2楼的内容
强烈赞扬,罕见的、系统的研发记录。在761设计中发现PIN管加偏压之后噪音剧烈增加,且噪音的增加比低...

感谢感谢,这阵会把相关研究都慢慢发到论坛里的。🙏

引用
评论
2
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
qiuzheru1
1年0个月前 修改于 1年0个月前 IP:瑞士
920688

有个疑问,你这设计中图23的TIA(或者说电荷放大器)用了100MOhm反馈电阻,这个值很大,带宽怕是对杂散电容也比较敏感。为什么不考虑选择更大的反馈电容来控制带宽,而是另外设计一级低通滤波?

如果系统分辨率受宽带噪声影响,而且不追求高计数率,在PD接近能买到最好的极限的情况下,把前级电荷放大的时间常数改大,是否对减小噪声能有帮助?当然,参考各大模拟芯片公司的Application note (搜索TIA noise,很多),仔细分析噪声来源是最好的。

引用
评论
1
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
黑体辐射作者
1年0个月前 修改于 1年0个月前 IP:湖北
920689
引用qiuzheru1发表于4楼的内容
有个疑问,你这设计中图23的TIA(或者说电荷放大器)用了100MOhm反馈电阻,这个值很大,为什么...

这个是因为JFET运放的特性,在几乎开环时,输入端不能浮空,因为有Vos的存在,在该电路中彻底浮空输入端会造成输出进入+-Vin。

关于噪声抑制,我后面还会发个帖子,是这个探测器的超低噪版本,到时候会把PCB和仿真一并开源,到时候欢迎您指点指点超低噪技术🙏


引用
评论(5)
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
巨大橘子
1年0个月前 IP:上海
920696

大佬就是大佬,一出手就是重量级,期待后续。

我有个小问题,关于7.1中硅光电管测得的Am-241光子能谱,为何69.8keV峰几乎和59.5keV一样大,按理说把Am-241从65-75keV区域里所有的峰加一块产额也只有0.0034%,和59keV处的35.9%相去甚远,在以往文献的SDD、HPGe以及闪烁体/正比管测试的谱图中,似乎也很少见到这么明显的峰,有没有可能是信号处理导致的(直球怀疑theremino和声卡部分)

另外估测活度可以考虑一下特征全能峰面积+效率+产额补偿嘛,看大佬的探头应该是2''NaI,这种尺寸的晶体的效率应该有超级多的资料可以参考,这个大小本底也有三位数cps了,还蛮大的。

引用
评论
1
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
黑体辐射作者
1年0个月前 修改于 1年0个月前 IP:湖北
920702
引用巨大橘子发表于6楼的内容
大佬就是大佬,一出手就是重量级,期待后续。我有个小问题,关于7.1中硅光电管测得的Am-241光子能...

大佬您好,我在B站关注您了,还参与您的抽奖活动了😁

1.这个完全是有可能的,有可能是开了FRI后的幻觉峰。这个是经过很长时间采样出来的数据。这个能谱的问题之前在论坛上求助过:AM241能谱求助。这个是原始数据:(没有标定的)

高能区修正-iir0.png

修正MQ87-AM241.png

   所以很有可能是幻觉,感谢提醒。

软件配置:

attachment icon 配置.txt 12.16KB TXT 2次下载

attachment icon 配置2.txt 12.16KB TXT 2次下载

采样数据:

attachment icon 数据1.txt 81.21KB TXT 2次下载

attachment icon 数据2.txt 82.27KB TXT 1次下载

attachment icon 数据3.txt 82.36KB TXT 1次下载

XXXX241能谱测试所采用的方案,是不加闪烁体直接测量的那个。背景计数非常低mac上显示的是0,背景计数是手算出来的,所以没能提供有效的背景能谱。感谢您提出的全峰面积法,这个我会改进的。


引用
评论
3
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
WHUlyq
11个月25天前 IP:甘肃
921250

看到参考导师的文献很神奇的感觉😂



引用
评论
1
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
巨大橘子
11个月24天前 IP:上海
921298
引用黑体辐射发表于7楼的内容
大佬您好,我在B站关注您了,还参与您的抽奖活动了😁。1.这个完全是有可能的,有可能是开了FRI后的...

好耶,感谢参与,祝大佬抽到好东西😂

这个怪怪的峰看起来是和其他峰同步冒出来的,那不妨试试用一些物体当做滤光片,在逐步改变射线不同能量的比例,看看各个峰的高度/面积有没有随之变化,效果应该蛮明显的。

大佬如果感觉thermino软件用起来有局限,推荐一发interspec,github上免费下载的,可以看到各个峰的参考比例什么的,功能很多很方便,打开的时候文件名应该是不能含有中文的,可以试试看。

引用
评论(1)
1
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
老母鸡火箭筒
11个月2天前 IP:江苏
922152

单一直流偏置不一定能满足需求,楼主选的pin管的即使在较高偏置电压下的结电容还是比较大的,会导致脉冲的响应没有那么快,一般用PIN管的话给的偏压会比较极限,接近雪崩电压点,而且会设计猝灭电路以保证每个计数的灵敏识别。


引用
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论

想参与大家的讨论?现在就 登录 或者 注册

文件下载
加载中...
{{errorInfo}}
{{downloadWarning}}
你在 {{downloadTime}} 下载过当前文件。
文件名称:{{resource.defaultFile.name}}
下载次数:{{resource.hits}}
上传用户:{{uploader.username}}
所需积分:{{costScores}},{{holdScores}}下载当前附件免费{{description}}
积分不足,去充值
文件已丢失

当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
支持的图片格式:jpg, jpeg, png
插入公式
评论控制
加载中...
文号:{{pid}}
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
收藏
取消收藏
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
管理提醒
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}