孙肖南，孙 陶，李江坤，张文峰，吴 雪
(1.核工业航测遥感中心，石家庄 050002；
2.铀资源地球物理勘查技术重点实验室，石家庄 050002；
3.河北省航空探测与遥感技术重点实验室，石家庄 050002)
天然放射性元素发出的伽玛射线特征峰能量固定，各特征峰数据在能谱测量中的记数道位置也固定；
但是在航空伽玛能谱测量中，对伽玛射线的接收主要采用NaI(Tl)晶体与光电倍增管电路，碘化钠晶体发光效率和光电倍增管的放大性能受温度变化的影响，在能谱测量工作中，测量数据会随着环境温度的变化而产生能谱数据记录道数的非线性漂移[1-2]，影响测量数据的准确性。因此在放射性伽玛能谱测量中，设计自动稳谱电路，减小能谱测量数据的漂移，提高测量数据的质量有着非常重要的意义；
伽玛能谱测量自动稳谱技术是在测量中利用软件寻找天然放射性元素具有的伽玛射线能量特征峰数据，计算中心峰位与相应的特征峰标准峰位道数偏差，跟据偏差数据控制稳谱电路调整能谱测量信号脉冲的大小，减小特征峰记录道数的漂移；
常用于稳谱的天然放射性元素特征峰有：40K 1.46 MeV、214Bi 1.76 MeV、208Th 2.615 MeV以及测试源137Cs 0.662 MeV[3]；
航空伽玛测量稳谱电路设计采用软件寻峰控制，根据测量数据特征峰位的漂移变化，利用软件程序调整控制能谱信号增益、调整信号电平相对基线的偏移变化，减小环境温度变化引起的能谱数据非线性漂移；
设计内容包括：1)高压调整控制电路：在稳谱软件控制下利用单片机调整光电倍增管高压，改变信号增益，稳定信号输出幅度；
2)信号脉冲电平叠加与消减电路：在稳谱软件控制下利用单片机操作信号脉冲电平叠加与消减，平移能谱信号脉冲幅度；
3)单片机控制电路：执行单片机与数据收录系统的数据通信和数据传输，跟据接收到的稳谱命令控制稳谱电路的工作状态[4]。
图1 电路系统结构框图
在能谱仪器测量工作中，伽玛射线照射到能谱探测器晶体上，在晶体内部产生电离，电离产生的荧光在光电倍增管阴极打出电子，电子在光电倍增管内部经多级倍增放大，在阳极产生电荷信号脉冲，其阳极信号的放大与光电倍增管的高压存在指数增长关系[5]，本设计电路利用stm32f407调整光电倍增管的高压控制光电倍增管信号脉冲输出的幅度；
电荷脉冲信号经放大电路进入信号脉冲幅度调整电路，脉冲幅度调整电路可以在信号脉冲幅度上进行电平电压的叠加，调整信号的增加或减小，利用stm32f407控制脉冲幅度调整电路调整信号对基线的偏移；
多道分析器电路将接收到的脉冲信号分析、处理，输出1 024道全谱数据；
并通过单片机电路发送到能谱数据收录系统；
在数据收录系统中运行的稳谱软件，循环不断的对接收到的1 024道全谱数据进行设定时间的计数累加，当设定的计数累加时间结束，大量的能谱测量累加数据，特征峰计数曲线符合高斯分布[6]，稳谱寻峰程序利用逐道比较法在1 024道能谱测量累加数据中寻找可用的特征峰以及特征峰数据的左右边界，计算特征峰位计数的中心峰位和中心峰位与标准峰位的偏移大小，以及特征峰位之间的相对记录数据位置的偏移量，计算出调整控制参数并通过串口发送到stm32f407控制单元，利用stm32f407控制稳谱电路执行相关功能的稳峰操作，稳定能谱特征锋位的数据记录，稳谱程序工作如图2所示[7-8]。
图2 稳谱程序流程框图
在连续的伽玛能谱测量数据计数累加中，各道能谱数据计数经一定时间的测量数据累加，累加数据中放射性元素能谱特征峰计数曲线符合正态分布统计规律，特征峰中心峰位计算如下：
nl为特征峰的左边界道值,nr为特征峰的右边界道值,Ci代表各道计数，i代表能谱道数。
稳谱常用的几种放射性元素能量及特征峰中心峰位记录道数如表1所示。
表1 放射性元素特征峰位
放射性能谱测量自动稳谱电路设计包括：1)stm32f407单片机数据传输和稳谱控制电路：负责接收多道能谱分析器发送的1 024道能谱数据，并通过串口将接收的能谱数据发送到数据收录系统，通过串口接收稳谱控制命令并执行稳谱命令的相关操作;2)由stm32f407控制的光电倍增管高压调整控制电路;3)由stm32f407控制的能谱脉冲信号加减电平电路;4)系统电源供电电路。
电路设计中，能谱探测器光电倍增管采用进口型号，高压工作范围为1 000～1 400 V；
高压电路设计采用国产高压电源模块HG-P152-1M66，12 V电压输入，最高输出1 500 V可调整高压。
2.1 stm32F407单片机数据传输和控制电路设计
稳谱控制处理电路采用ST公司Cortex M4内核stm32f407单片机及外围电路构成；
利用stm32f407单片机接收能谱多道分析器电路发送的1 024道全谱数据，并通过串口将接收到的全谱数据发送到数据收录系统，同时stm32F4单片机通过串口接收数据收录系统发来的稳谱控制命令，根据稳谱命令内容控制光电倍增管的高压调整和控制信号脉冲幅度调整电路的电平调整数值；
STM32F407单片机设计电路如图3所示。
图3 stm32f407单片机控制电路电路设计图
stm32f407单片机通过串口PA10实时发送1 024道全谱数据到收录系统，同时不断监测串口的接收状态，通过串口PA9、接收数据收录系统稳谱程序发送来的稳谱高压调整指令和峰位调整指令；
单片机根据指令中调整数据的大小 通过PA4和PA5口输出相应的模拟控制电压，输出的模拟信号DAC2调整光电倍增管工作高压，DAC1调整能谱接收信号脉冲幅值增量，单片机通过PD0和PD1接口分别操纵模拟开关SGM3001_A和SGM3001_B控制能谱数据信号脉冲幅度基值的增加或减小[9-10]。
2.2 光电倍增管高压调整控制电路设计
光电倍增管高压调整由程序命令进行控制，高压调整电路设计如图4所示，高压模块选用HG-P152-1M66，高压输出为0～1 500 V，高压调整控制电压为0～5 V；
设计电路选用光电倍增管高压工作范围在1 000～1 400 V，考虑到所使用光电倍增管的型号更换，增强电路的兼容性，设计电路高压调整范围为700～1 400 V；
光电倍增管高压调整控制是由stm32f407单片机D/A模拟输出来实现，设置单片机控制电压输出范围为：0～2.5 V，由J2插口送入高压调整控制电路，与MAX6062输出的基准电压叠加，控制高压模块的调整电压输入端，电路中放大电路芯片采用AD8065，AD8065输出端电压设计符合下式规律：VHV-ADJ=VD/A+2，输出的调整控制电压2～4.5 V，对应高压输出范围为650～1 400 V。
图4 高压调整电路设计图
2.3 能谱信号脉冲幅度加减电平电路
能谱信号脉冲幅度加减电平调整由程序命令控制，设计电路如图5所示。信号脉冲幅度调整设计为：信号脉冲到来时允许进行幅度调整，信号脉冲消失后保持基线状态，信号脉冲幅度调整功能不起作用；
电路设计主要包括放大电路、信号比较电路和高速电子开关电路构成，所采用的芯片主要有：AD8065、MAX913、SGM7222和SGM3001。电路利用接收到的信号脉冲通过MAX913比较器产生与接收到的信号脉冲同步的方波信号，控制SGM7222电子开关的导通与关闭，保证只有信号脉冲到来时，才起调节作用；
SGM7222模拟开关控制信号为高电平时，将单片机控制输出的模拟电压接入信号脉冲放大电路，模拟电压由稳谱控制控制命令数据转换得到，模拟电压作为信号脉冲幅度增加或减小的调整值；
SGM7222模拟开关控制信号为低电平时，将地电平接入信号电路，使得电路基线电平保持不变。能谱信号脉冲幅度的调节大小及方向由运行在数据收录系统的自动稳谱程进行控制，STM32f407单片机通过串口接收发送来的稳谱命令，操纵脉冲幅度调整电路执行相应的命令动作，单片机通过对两个电子开关SGM3001的控制操作来实现信号脉冲幅度增加或减小的操作功能，单片机PD0口操纵的SGM3001电子开关控制信号脉冲幅度的增加操作，PD1口连接的SGM3001电子开关控制信号基脉冲幅度的减小操作[11-12]。
图5 信号脉冲幅度调整电路设计图
2.4 电源电路
设计电路供电采用12 V电压输入，12 V电压分别供给高压电源模块HG-P152-1M66和开关电源芯片MP1482，通过电源芯片MP1482及相关元件电路得到5 V电路工作电压，其输出电流可达2 A，效率93%，工作温度范围-40～85 ℃。
自动稳谱程序在收录系统中不断的对接收到的放射性全谱数据进行累加，在设定的数据累加时间结束，寻峰程序采用逐道比较法寻找适合于稳谱控制的特征峰，计算特征峰中心峰位的数据偏移，以及各特征峰位的相对偏移数据，通过串口向stm32f407单片机控制电路发送相应的能谱特征峰数据计数偏离调整命令，stm32f407单片机控制和操纵稳谱电路执行相应命令的数据操作；
自动稳谱电路控制采用stm32f407单片机，stm32f407软件编程包括单片机系统的初始化程序，主要有系统时钟初始化、GPIO接口初始化、DMA通道初始化、中断初始化、以及串口初始化和模拟端口的初始化，单片机初始化程序完毕后自动进入数据传输状态，将从能谱多道分析器电路接收到的1 024道全谱数据，通过串口发送到数据收录系统，同时不断的监测串口输入，通过串口接收数据收录系统发出的稳谱操作指令，根据稳谱控制命令执行相应的电路操作。stm32f407单片机接收稳谱数据指令，调整控制程序设计如图6所示。
图6 stm32f407单片机控制程序设计框图
stm32f407单片机通过PC11引脚接收多道能谱分析器发出的1 024道全谱数据，通过DMA通道传输到串口发送缓冲区，按照相应的数据格式通过串口发送引脚PA9实时向收录系统发送所接收的1 024道全谱数据，同时串口接收引脚PA10实时的等待接收收录系统稳谱程序发来的稳谱控制命令，命令格式设计为：高压调整命令为：Vdddd 加回车换行符结束，其中“V”为英文大写字母，后跟的“dddd”为4位要调整的十进制数据，0 系统测试利用串口软件程序通过串口分别发送高压调整命令和脉冲信号幅度调整命令，检查调整电路对调整参数的反应状态。 光电倍增管工作点设定高压为1 207 V，对应的调整数据为2 178，以2 178为调整数据中心点，连续的增加和减小数值，测试高压的变化； 表2 调整参数与高压测试表 由测量数据得到调整数据与高压输出线性趋势线：y=0.263 3x+633.41，线性相关系数R2=0.999 7。 由高压调整测试数据表显示，光电倍增管高压由单片机控制可以进行0.2～0.3伏的变化的调整，对于钾道1.76 MeV的特征峰有0.4道左右变化量的调整，且测试数据的线性度比较好，满足高压调整设计的要求。 1 024道能谱数据特征峰道数偏移与调整参数对应关系测试，利用信号幅度调整命令通过串口发送调整指令，采用改变调整数据参数测量调整输出电压进行测试，电路设计中：每道信号电压=2 000 mV/1 024=1.953 mV/道，利用测试电压与道数偏移的线性关系计算出偏移道数：道数偏移=V/1.953 V为测试电压，单位为毫伏。脉冲幅度调整参数与道偏移测试如表3所示。 表3 脉冲幅度调整参数与能谱偏移道数测试表 由测量数据得到峰位调整参数与能谱道数偏移趋势线：y=0.406 1x+406.12，线性相关系数R2=0.999 9。 由脉冲幅度调整测试数据可看出脉冲幅度调整，一个调整数字的增量改变可使能谱峰位有0.4道左右的变化，调整作用明显且测试数据的线性度很好，满足信号脉冲峰位调整的设计要求，设计电路可以实现航空伽玛能谱测量自动稳谱控制电路峰位调整功能。 在航空伽玛能谱测量中，能谱自动稳谱电路对于减小测量数据的漂移，提高伽玛能谱测量数据的质量、起着非常重要的作用；4.1 光电倍增管高压与调整参数对应关系测试
高压调整参数测试数据如表2所示。
4.2 信号脉冲增加消减电平电路测试
对于温差较大时，能谱数据的漂移往往是非线性的[1]，本稳谱电路设计采用软件控制，利用改变光电倍增管高压的方式调整信号增益，在能谱测量信号脉冲上叠加控制电平、调整能谱信号相对基线的漂移变化，结合两种调整方式可以控制能谱峰位的漂移变化。在稳谱电路的测量试验中，对于放射性元素137Cs、40K、214Bi、208Th特征峰位测量数据线性拟合R2=0.999 9，1 024道能谱测量数据特征峰最大峰漂1.5道，达到了稳谱电路设计的要求[10]。稳谱软件程序对光电倍增管的高压控制可以达到0.2～0.3 V调整变化；
软件程序对信号幅度增量的调整控制，由于受到单片机D/A输出精度限制，信号幅度调整电路最小调整增量电平为0.9 mV，为提高调整精度设计可以采用更高精度的外部D/A电路。