Small Pixel 10 μm Pitch Infrared Focal Plane Array ROIC Design
-
摘要: 研制出一款小像元10μm中心距红外焦平面探测器CMOS（complementary metal oxide semiconductor）读出电路ROIC（read out integrated circuit）。读出电路设计包括积分后读出（integration then reading，ITR）和积分同时读出（integration while reading，IWR）模式，ITR模式下有2档增益，电荷满阱容量分别为4.3 Me-和1.6 Me-，其他功能包括抗晕、串口功能控制以及全芯片电注入测试功能。读出电路采用0.18μm工艺，电源电压3.3 V，测试结果表现出良好的性能：在77 K条件下，全帧频100 Hz，读出电路噪声小于0.2 mV。本文介绍了该款读出电路设计的基本架构，分析了在小的积分电容下电路抗干扰能力的设计。在测试过程中，发现了盲元拖尾现象，分析了拖尾现象产生的原因，为解决拖尾现象设计了抗晕管栅压产生电路，最后给出了整个电路的测试结果。Abstract: A kind of infrared focal plane CMOS (complementary metal oxide semiconductor) read out integrated circuit (ROIC) for small pixel applications was developed. This ROIC design includes ITR(integration then reading) and IWR(integration while reading)functions, two gains in ITR mode, charge capacities are 4.3 Me- and 1.6 Me-. Other functions include anti-blooming, series port control function and full chip current injection test function. This ROIC was fabricated in 0.18μm process, power supply voltage 3.3 V, test result show good performance of the ROIC, full frame rate is 100 Hz, noise of readout circuit is 0.2 mV. This paper introduces the basic structure of the readout circuit design, and analyzes the design of anti-interference ability of the circuit under the condition of small integral capacitance; During the test, tailing phenomenon is found, the causes of the tailing phenomenon is analyzed and In order to solve the tailing phenomenon, the Anti-blooming voltage generation circuit is designed, and the test results of the whole circuit are given at the end of the paper.
-
Keywords:
- small pixel /
- ROIC /
- infrared focal plane detector
-
0. 引言
湿法烟气脱硫利用石灰石浆液吸收烟气中的SO2，具有高效率和高可靠性等优势，已成为当前火电厂和化工厂脱硫的主力技术[1]。但是火电厂运行过程中，由于煤燃烧和SO2氧化，特别是在选择性催化还原脱销系统中催化剂作用下，SO2会更有利于氧化成生成SO3，造成烟气中SO3浓度显著增加[2-3]。烟气进入脱硫塔系统后，温度迅速冷却至酸露点以下，形成细小SO3酸雾，而单脱硫塔的SO3脱除效率仅为30%~40%[4]。
烟气中的SO3、SO2，HF及其它酸性物质会导致塔体金属发生化学腐蚀，脱硫塔内还存在电化学腐蚀、磨蚀、结晶腐蚀、垢下腐蚀和氯离子腐蚀[5]。在以上多种方式的共同作用下，受到内涂玻璃鳞片、聚烯烃共聚物、改性聚脲或纳米复合涂料等物质防护的脱硫塔仍可能发生腐蚀。
目前，对脱硫塔的检测主要有外观损伤、钢材厚度、力学性能、焊缝缺陷检测和构件变形等项目[6]。以上检测项目能够有效对停机后的脱硫塔健康状态进行评测。但是，对运行中的脱硫塔腐蚀状态进行有效检测未见文献报道。结合脱硫塔运行参数和结构参数，本文提出了采用传热学反演的方法根据表面红外热像进行运行中脱硫塔壁厚定量检测，并以某厂的烟气脱硫塔为对象进行了实验验证。
1. 脱硫塔传热模型
研究对象为图 1所示的脱硫塔，其内部环境复杂，无法布置有效的温度测量装置；脱硫过程是一个包含了传热、传质和化学反应的复杂过程，难以准确地用数学语言描述。为了建立脱硫塔的传热模型，本文进行以下简化：
1）脱硫塔为轴对称结构，内壁面热流沿周向分布均匀；
2）脱硫塔内部有玻璃鳞片防腐涂层，涂层质地均匀且热物性参数为各向同性；
3）烟气和石灰浆液对脱硫塔的传热，可等效为对脱硫塔内壁施加有沿轴向分布的加热热流；
4）忽略脱硫塔内部喷嘴和支撑结构对脱硫塔壁面温度分布的影响；
5）忽略脱硫塔内壁和防腐涂层的接触热阻；
6）忽略脱硫塔的轴向导热；
7）脱硫塔的温度场为稳态。
脱硫塔浆液区是腐蚀的重灾区，本文以浆液区段为研究对象，简化后的脱硫塔浆液区如图 2所示，具体几何参数及热物性参数见表 1。脱硫塔内壁s1存在轴向分布热流，外表面s2与环境存在对流换热，且s2面上的温度可以直接测量，其他表面s3绝热。记x＝(ϕ, θ, z)为空间坐标向量，脱硫塔温度场的控制方程为：
表 1 脱硫塔几何参数及热物性参数Table 1. Geometric parameters and thermophysical parameters ofdesulfurization towerMaterial Thickness /
mmDiameter/
mmThermal conductivity/
[W/(m·K)]Anticorrosive
coating4 5520 0.35 The tower wall 14 5528 48.85 $$\frac{1}{r}\frac{\partial }{{\partial r}}(\lambda (\mathit{\boldsymbol{x}})r\frac{{\partial T(\mathit{\boldsymbol{x}})}}{{\partial r}}) + \frac{1}{{{r^2}}}\frac{\partial }{{\partial \varphi }}(\lambda (\mathit{\boldsymbol{x}})\frac{{\partial T(\mathit{\boldsymbol{x}})}}{{\partial \varphi }}) + \frac{\partial }{{\partial z}}(\lambda (\mathit{\boldsymbol{x}})\frac{{\partial T(\mathit{\boldsymbol{x}})}}{{\partial z}}) = 0$$ (1) 边界条件为：
$$ - \lambda (\mathit{\pmb{x}})\frac{{\partial T(\mathit{\pmb{x}})}}{{\partial n}} = q(z)\;\quad \mathit{\pmb{x}} \in {s_1}$$ (2) $$ - \lambda (\mathit{\pmb{x}})\frac{{\partial T(\mathit{\pmb{x}})}}{{\partial n}} = h(T(\mathit{\pmb{x}}) - {T_f})\;\quad \mathit{\pmb{x}} \in {s_2}$$ (3) $$ - \lambda ({\mathit{\pmb{x}}})\frac{{\partial T({\mathit{\pmb{x}}})}}{{\partial n}} = 0\;\;\;\;\;{\mathit{\pmb{x}}} \in {s_3}$$ (4) 式中：q(z)为内壁沿轴向z的热流密度分布；h为表面s2的对流换热系数；Tf为环境温度；n为表面的外法线方向。
若已知脱硫塔的几何结构和热边界条件，利用有限元法（finite element method, FEM）求解公式(1)~(4)，可以确定脱硫塔的温度场T(x)，作为壁厚检测的基础。
2. 基于反问题的脱硫塔壁厚检测
运行中脱硫塔壁厚d是影响脱硫塔表面温度分布的关键因素之一，如果脱硫塔壁面腐蚀减薄，则传热热阻减小，在腐蚀部位对应的塔外表面形成局部高温区域，如图 1(b)所示。因此，可通过求解导热反问题根据脱硫塔表面红外热像进行脱硫塔壁厚检测。
2.1 壁厚检测方案
在如图 2所示的脱硫塔表面热像图中，选取高温区域的K个温度测点，依据此测量信息Tkmea（k＝1, 2, …, K）采用共轭梯度方法（conjugate gradient method，CGM）求解多变量稳态传热反问题[7-8]，进行壁厚d的反演；由于脱硫塔为薄壁结构，热扩散效应较弱，则外表面高温区域可看作与内壁腐蚀区域接近。为了便于问题讨论，在本文中，内壁腐蚀区域用圆柱近似。
然而，壁厚反演过程中正问题的计算需要已知脱硫塔内壁热边界条件如热流q(z)。而内壁热流q(z)难以直接测量。如果直接同时反演壁厚和内壁热流q(z)，可能因为测量信息不能够同时对壁厚和热流具有较大的灵敏度，造成检测系统的病态程度加剧。
因此，检测方案包括了两步：先进行内壁热流定量识别，再定量识别壁厚。内壁热流可采用CGM反演得到：在脱硫塔外表面高温区域附近同等高度的温度正常区域沿周向选取M个测点，以该测点的温度信息Tmmea（m＝1, 2, …, M）反演该位置处的脱硫塔内壁热流q；内壁热流沿高度方向变化，但是周向分布均匀，以第一步反演得到的内壁热流q作为壁厚反演中正问题的已知热边界条件，以高温区域的温度测量信息作为壁厚反演的依据，提高了测量信息对壁厚的灵敏度，有利于削弱壁厚检测问题的病态程度。
2.2 共轭梯度算法
利用共轭梯度算法求解壁厚反问题，通过迭代优化使得目标函数J(d)足够小或者达到最大迭代步imax，对应的壁厚d即为所求。目标函数J(d)可表示为：
$$J({d_i}) = \sum\limits_{k = 1}^K {{{[T_k^{{\rm{cal}}}({d_i}) - T_k^{{\rm{mea}}}]}^2}} \le \varepsilon $$ (5) 式中：Tkmea为在脱硫塔红外热像图上提取的第K个温度测量值；di为第i次迭代得到的壁厚的猜测值；Tkcal(di)是根据di进行正问题计算得到的第k个测量位置处的温度计算值。K为在红外热像图上提取的温度测量值的数目。停机标准ε可由下式表示：
$$ \varepsilon = K{\sigma ^2} $$ (6) 式中：σ为测量误差的标准差。
CGM沿着已知点处的梯度所构造出的共轭方向迭代搜索目标函数的极小点，迭代过程中对壁厚猜测值的修正可表示为：
$$ {d_i}_{ + 1} = {d_i} - {\alpha _i}{\mathit{\boldsymbol{\gamma }}_i} $$ (7) 式中：αi为搜索步长；γi为搜索方向。
搜索步长αi表示为：
$${\alpha _i}{\rm{ = }}\sum\limits_{k = 1}^K {[T_k^{{\rm{cal}}}({d_i}) - T_k^{{\rm{mea}}}]\nabla T_k^{{\rm{cal}}}({d_i}){\mathit{\boldsymbol{\gamma }}_i}} /\sum\limits_{k = 1}^K {\nabla T_k^{{\rm{cal}}}({d_i}){\mathit{\boldsymbol{\gamma }}_i}} $$ (8) 搜索方向γi可由下式表示：
$${\boldsymbol{\gamma} _i}{\rm{ = }}\nabla J({d_i}) + {\beta _i}{d_{i - 1}}$$ (9) 式中：▽J(di)为目标函数的梯度；βi为共轭系数，可根据式(10)计算：
$${\beta _i}{\rm{ = }}{\left[ {\nabla J({d_i})/\nabla J({d_{i - 1}})} \right]^2}$$ (10) 利用CGM根据红外热像图中正常区域温度反演该位置处的脱硫塔内壁热流q，其过程可参考公式(6)~(10)，在此就不一一赘述。
2.3 迭代求解流程
应用CGM根据脱硫塔红外热像图反演壁厚的计算步骤如下：
1）根据红外热像图，对异常区域进行辨识；
2）反演异常区域脱硫塔内壁热流q；
3）给出壁厚初始猜测值d0；
4）通过求解公式(1)~(4)，得到测点处的计算温度Tkmea（k＝1, 2, …, K），并代入公式(5)：
$$J({d_i}) = \sum\limits_{k = 1}^K {{{[T_k^{\rm{cal}}({d_i}) - T_k^{\rm{mea}}]}^2}} \leqslant \varepsilon $$ (11) 如果满足上述条件，di即为所求，停止迭代；否则继续；
5）按公式(8)~(10)对CGM里的参数进行更新；
6）根据公式(7)更新壁厚d的猜测值，并返回步骤4）。
3. 脱硫塔检测结果及分析
选取环境温度Tf＝20℃，对流换热系数h＝10 W/(m2·K)。异常区域温度测点数量K＝3；正常区域温度M＝2。考虑到实际测量过程中，温度测量误差是无法消除的，通过现场标定，测量误差σ＝0.055℃，ε＝0.01。
1）数值实验验证
为了验证检测系统的有效性和精确性，本文先进行脱硫塔内部缺陷检测的数值实验。在数值实验中，脱硫塔的热边界条件、几何参数和热物性参数均与实际过程相同，假设真实壁厚de＝16 mm。
数值实验中设置不同大小的测量误差，以考察测量误差对缺陷检测结果的影响。实验结果如表 2所示。
表 2 不同测量误差时的检测结果Table 2. The detection results of the different measurementerrorsMeasurement
error σDetect wall
thickness d/mmRelative
error/ %0.055 15.92 0.50 0.1 16.35 2.19 0.2 17.14 7.13 从表 2可以看出，随着测量误差的增大，壁厚检测结果的精确性下降。如σ＝0.2℃时，相对误差为7.13%，在工程上可以接受。
2）依据现场红外热像反演
① 高温区域1
如图 3所示，高温区域1的最高温度为48.2℃，区域的最大温差为1.4℃。利用基于导热反问题的脱硫塔壁厚检测方法，对高温区域1的壁厚进行计算，结果为d＝14.6 mm，即该区域玻璃鳞片厚度为0.6 mm。考虑到计算误差，可判定为防腐涂层已磨损殆尽或脱落，若不处理，塔壁金属将受到快速腐蚀。在检测后30天左右，脱硫塔停机检修，发现高温区域1的防腐涂层已脱落，证实了本方法的正确性。
![]()
图 3 高温区域1对应的红外热像图Figure 3. Infrared image of the high temperature region 1(a) Infrared image of the tower (b) Local enlargement picture② 高温区域2
如图 4所示，高温区域2的最高温度为48.3℃，区域的最大温差为1.2℃。对高温区域2的壁厚进行反演：d＝15.8 mm，即该区域防腐涂层厚度为1.8 mm，可判定为防腐涂层已减薄。
③ 高温区域3
如图 5所示，高温区域3的最高温度为48.1℃，区域的最大温差为3.7℃。对高温区域3的壁厚进行反演：d＝6.5mm。表明该区域脱硫塔金属塔壁已发生腐蚀，减薄了7.5mm，应尽快排查、检修。
4. 结论
本文采用导热反问题的方法，根据红外热像图对运行中的脱硫塔壁厚进行了检测。其中，脱硫塔壁厚和脱硫塔内壁面热流的反演均采用共轭梯度法。首先通过数值实验，验证了本方法的可行性。然后，依据红外热像进行反演，发现脱硫塔筒体腐蚀1处，防腐涂层脱落1处，防腐涂层减薄1处。在后续的停机检修时对上述部进行了复核，均验证了上述检测结果，表明了基于表面红外热像的脱硫塔壁厚定量检测方法的有效性和准确性。
防腐涂层厚度的不一致，脱硫塔内介质分布的不均匀，红外热像仪精度以及环境等因素，可能会给壁厚检测结果引入误差，造成识别精度下降。如何提高壁厚检测精度，仍是下一步研究的方向。
-
表 1 法国Sofradir公司推出的Daphnis产品信息
Table 1 Daphnis product information by Sofradir
Detector spectral
response3.7-4.8 μm FPA operating
temperatureUp to 110 K ROIC architecture Digital outputs
Direct injection input circuitROIC functionalities Programmable integration time
anti-blooming
Invertrevert
Bining
IWR or ITR
High dynamic rangeWindowing modes 320×4 minimum programmable Charge handling
capacity3main
Gains: 1.1Me-; 2.7Me-; 5.6Me-Frame rate Up to 85 Hz full frame rate NETD 20 mK(293 K, 70% well fill, 2.7Me-) 
下载: 导出CSV
表 2 读出电路主要性能参数
Table 2 Main performance parameters of readout circuit
Parameters Typical value Array 1024×768 Pixel pitch 10μm Main clock 10 MHz Charge capacity Gain1：4.3 Me- Gain0：1.6Me- Output channel 8 Output voltage swing 2 V(1-3 V) Readout mode ITR/IWR
下载: 导出CSV
表 3 沟道电荷注入效应引起的电压变化
Table 3 Voltage change caused by channel charge injection effect
C/fF ΔV/mV 800 7.925 400 15.85 200 31.7 100 63.4 50 126.8 10 634
下载: 导出CSV
表 4 列放大器对保持电容影响
Table 4 Effect of column amplifier on holding capacitance
C/fF ΔV/mV ΔV'/mV 10 629.77 21.75 40 347.35 8.32 80 204.18 2.26 100 165.81 1.91 400 43.42 0.56 800 21.88 0.29
下载: 导出CSV
表 5 读出电路测试结果
Table 5 Read out circuit test result
IWR/ITR function Normal Window mode test Normal Serial electrical interface Normal Output voltage swing 2 V Power dissipation 150 mW Output bandwidth 10 MHz Frame rate 100 Hz Noise 0.2 mV
下载: 导出CSV
-
[1] 杨超伟, 李东升, 李立华, 等. 小像元碲镉汞红外焦平面探测器的研究进展[J]. 红外技术, 2019, 41(11): 1003-1011. https://www.cnki.com.cn/Article/CJFDTOTAL-HWJS201911003.htm YANG Chaowei, LI Dongsheng, LI Lihua, et al. Review of small-pixel HgCdTe infrared focal plane detector[J]. Infrared Technology, 2019, 41(11): 1003-1011. https://www.cnki.com.cn/Article/CJFDTOTAL-HWJS201911003.htm
[2] 李俊斌, 李东升, 杨玉林, 等. 以色列SCD公司的III-Ⅴ族红外探测器研究进展[J]. 红外技术, 2018, 40(10): 936-945. http://hwjs.nvir.cn/article/id/hwjs201810003 LI Junbin, LI Dongsheng, YANG Yulin, et al. III-V semiconductor infrared detector research in SCD of Israel[J]. Infrared Technology, 2018, 40(10): 936-945. http://hwjs.nvir.cn/article/id/hwjs201810003
[3] 邓功荣, 赵鹏, 袁俊, 等. 锑基高工作温度红外探测器研究进展[J]. 红外技术, 2017, 39(9): 780-784. http://hwjs.nvir.cn/article/id/hwjs201709002 DENG Gongrong, ZHAO Peng, YUAN Jun, et al. Status of Sb-based HOT infrared detectors[J]. Infrared Technology, 2017, 39(9): 780-784. http://hwjs.nvir.cn/article/id/hwjs201709002
[4] Beletic J W, Blank R, Gulbransen D, et al. Teledyne imaging sensors: infrared imaging technologies for astronomy and civil space[C]//High Energy, Optical, and Infrared Detectors for Astronomy III, 2008, 7021: 70210H.
[5] Rogalski A. Recent progress in infrared detector technologies[J]. Infrared Physics & Technology, 2011, 54(3): 136-154. http://www.researchgate.net/profile/Antoni_Rogalski/publication/241112884_Recent_progress_in_HgCdTe_infrared_detector_technology/links/5524e12e0cf22e181e73b04e.pdf
[6] Rogalski A. Next decade in infrared detectors[C]//Electro-Optical and Infrared Systems: Technology and Applications XIV, 2017, 10433: 104330L.
[7] Caulfield J, Curzan J. Small pixel infrared sensor technology[C]//Infrared Technology and Applications XLIII, 2017, 10177: 1017725.
[8] Chen T, Catrysse P B, El Gamal A, et al. How small should pixel size be?[C]//Sensors and Camera Systems for Scientific, Industrial, and Digital Photography Applications, 2000, 3965: 451-459.
[9] Farrell J, Xiao F, Kavusi S. Resolution and light sensitivity tradeoff with pixel size[C]//Digital Photography II, 2006, 6069: 60690N.
[10] Lutz H, Breiter R, Eich D, et al. Small pixel pitch MCT IR-modules[C]//Infrared Technology and Applications XLII, 2016, 9819: 98191Y.
[11] Lutz H, Breiter R, Eich D, et al. Towards ultra-small pixel pitch cooled MW and LW IR-modules[C]//Infrared Technology and Applications XLIV, 2018, 10624: 106240B.
[12] Espuno L, Pacaud O, Reibel Y, et al. A new generation of small pixel pitch/SWaP cooled infrared detectors[C]//Electro-Optical and Infrared Systems: Technology and Applications XII, 2015: 9648: 96480H.
[13] Manissadjian A, Rubaldo L, Rebeil Y, et al. Improved IR detectors to swap heavy systems for SWaP[C]//Infrared Technology and Applications XXXVIII, 2012, 8353: 835334.
[14] Johnson John. Analysis of image forming systems[C]//Proceeding of SPIE- The International Society for Optical Engineering, 1958, 513(513) : 761.
[15] 周立庆, 宁提, 张敏, 等. 10 µm像元间距1024×1024中波红外探测器研制进展[J]. 激光与红外, 2019, 49(8): 915-920. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201908002.htm ZHOU L, NING T, ZHANG M, et al. Developments of 10 µm pixel pitch 1024×1024 MW infrared detectors[J]. Laser & Infrared, 2019, 49(8): 915-920. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201908002.htm
[16] Berthoz J, Rubaldo L, Maillard M, et al. MTF performance: measurements, modelisation, and optimization for Sofradir II-VI IR photodetectors[C]//Quantum Sensing and Nanophotonic Devices XII, 2015, 9370: 93700O.
[17] Reibel Y, Augey T, Verdet S, et al. High-performance and long-range cooled IR technologies in France[C]//Infrared Technology and Applications XXXIX, 2013, 8704: 87040B.
[18] Reibel Y, Rubaldo L, Manissadjian A, et al. High-performance MCT and QWIP IR detectors at Sofradir[C]//Electro-Optical and Infrared Systems: Technology and Applications IX, 2012, 8541: 85410A.
[19] Reibel Y, Rouvie A, Nedelcu A, et al. Large format, small pixel pitch and hot detectors at SOFRADIR[C]//Electro-Optical and Infrared Systems: Technology and Applications X, 2013, 8896: 88960B.
[20] Lefoul X, Pere-Laperne N, Augey T, et al. New SOFRADIR 10 µm pixel pitch infrared products[C]//Electro-Optical and Infrared Systems: Technology and Applications XI, 2014, 9249: 924911.
[21] Tan C L, Mohseni H. Emerging technologies for high performance infrared detectors[J]. Nanophotonics, 2018, 7(1): 169-197. http://www.onacademic.com/detail/journal_1000040103679310_e3b6.html
[22] Beletic J W, Blank R, Gulbransen D, et al. Teledyne imaging sensors: infrared imaging technologies for astronomy and civil space[C]//High Energy, Optical, and Infrared Detectors for Astronomy III, 2008: 70210H.
[23] Dorn R J, Eschbaumer S, Hall D N, et al. Evaluation of the Teledyne SIDECAR ASIC at cryogenic temperature using a visible hybrid H2RG focal plane array in 32 channel readout mode[C/]//Proc. of SPIE, 2008: DOI:10.1117/12.788717.
[24] Bai Y, Bajaj J, Beletic J W, et al. Teledyne imaging sensors: silicon CMOS imaging technologies for x-ray, UV, visible, and near infrared[C]//Proc. of SPIE, 2008, 7021: 702102.
[25] 毕查德·拉扎维. 模拟CMOS集成电路设计[M]. 西安: 西安交通大学出版社, 2003. Behzad Razavi. Design of Analog CMOS Integrated Circuits[M]. Xi'an: Xi'an JiaoTong University Press, 2003.
-
期刊类型引用(1)
1. 单幼芳. 基于深度学习算法的太赫兹人体安检图像处理技术研究. 电脑编程技巧与维护. 2024(07): 127-129 . 
百度学术
其他类型引用(1)
计量
- 文章访问数: 337
- HTML全文浏览量: 60
- PDF下载量: 125
- 被引次数: 2
