基于像差补偿的近红外显微干涉硅通孔测量

基于像差补偿的近红外显微干涉硅通孔测量 Measurement of through silicon via by near-infrared micro interferometry based on aberration compensation first-author 吴春霞（1997-），女，江苏扬州人，硕士研究生，2020年于金陵科技学院获得学士学位，主要研究方向为光学干涉精密测量。E-mail： 15651828978@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927356&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927358&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927357&type= 22.01332855 32.00400162 吴春霞 （1997-），女，江苏扬州人，硕士研究生，2020年于金陵科技学院获得学士学位，主要研究方向为光学干涉精密测量。E-mail： 15651828978@163.com 15651828978@163.com 高志山（1966-），男，江苏淮安人，教授，1992年于中国科学院长春光学精密机械研究所获得硕士学位，2000年于南京理工大学获得博士学位，主要研究方向为精密光学测试、先进光学设计。E-mail： zhishgao@njust.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927359&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927361&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927360&type= 22.01334000 32.00400162 高志山 （1966-），男，江苏淮安人，教授，1992年于中国科学院长春光学精密机械研究所获得硕士学位，2000年于南京理工大学获得博士学位，主要研究方向为精密光学测试、先进光学设计。E-mail： zhishgao@njust.edu.cn zhishgao@njust.edu.cn corresp E-mail： guozy15@njust.edu.cn E-mail： guozy15@njust.edu.cn 郭珍艳（1986-），女，陕西宝鸡人，副教授，2014年于南京理工大学获得博士学位，主要研究方向光学干涉精密测量等。E-mail： guozy15@njust.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927362&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927365&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927363&type= 22.01332855 32.00400162 郭珍艳 （1986-），女，陕西宝鸡人，副教授，2014年于南京理工大学获得博士学位，主要研究方向光学干涉精密测量等。E-mail： guozy15@njust.edu.cn guozy15@njust.edu.cn 南京理工大学 电子工程与光电技术学院，江苏 南京 210094 南京理工大学 电子工程与光电技术学院，江苏 南京 210094 School of Electronic and Optical Engineering， Nanjing University of Science and Technology， Nanjing 210094， China School of Electronic and Optical Engineering， Nanjing University of Science and Technology， Nanjing 210094， China 为了实现深宽比大于6∶1的硅通孔的形貌测量，提出了一种基于像差补偿的近红外显微干涉检测方法。该方法采用近红外宽带光作为光源，能够穿透硅通孔，检测系统中内置变形镜自适应像差补偿模块，主动补偿硅通孔引入的调制像差。在检测硅通孔三维形貌时，依据COMSOL Multiphysics有限元仿真软件得到的三维硅通孔高深宽比结构对探测光的调制像差规律，指导设置需变形镜补偿的像差种类和量值，用基于频域的评价函数指标阈值，判定硅通孔底部图像的聚焦状态，获得待测硅通孔清晰的底部图像，本质上提升探测光的重聚焦能力。在此基础上，使用垂直扫描干涉法得到待测硅通孔的深度与其三维形貌分布。实验测量了直径为10 μm、深度为65 μm、深宽比为6.5∶1和直径为10 μm、深度为103 μm、深宽比为10.3∶1的硅通孔深孔，并与高精度SEM的测量结果对比，深度测量的相对误差为1%。与白光显微干涉测量结果对比表明，本文所提出的方法可以获得清晰的高深宽比硅通孔的底部图像，有效增强底部的宽谱干涉信号和对比度，能够准确测量更高深宽比硅通孔的三维形貌。 Through silicon via （TSV） technology refers to the silicon micro-hole processing technology used to achieve vertical conduction and interconnection between chips. The three-dimensional （3D） morphology of TSV samples is commonly measured using destructive scanning electron microscope （SEM） profile imaging technology. In addition， white light interference technology has the advantage of non-destructive measurement； however， it is difficult to probe the light reaching the bottom of TSV samples with an aspect ratio higher than 6：1， and the topography is also distorted. To address this problem， a near-infrared micro-interference detection method based on aberration compensation is proposed in this paper. The method uses near-infrared broadband light as the probe light source， which can penetrate through the TSV. Moreover， an adaptive aberration compensation module in the form of a deformable mirror is introduced into the detection system to compensate for the modulation aberration that is synchronously generated by the TSV. The first step of detecting the TSV 3D topography is to set the type and quantity of the aberration that needs to be compensated by the deformable mirror. This is obtained using the finite element simulation software COMSOL Multiphysics based on the aberration modulation law for a 3D TSV high aspect ratio structure under the probe light. Subsequently， the evaluation function index threshold based on the frequency domain is used to estimate the focus state of the TSV bottom image. Furthermore， the refocusing ability of the probe light is essentially improved and a clear bottom image of the tested TSV is obtained. On this basis， the depth value of the tested TSV and its 3D topography distribution is calculated using vertical scanning interferometry. Finally， two types of TSV samples with a diameter of 10 μm and depth of 65 μm （aspect ratio 6.5）， and diameter of 10 μm and depth of 103 μm （aspect ratio 10.3）， are measured using the proposed method. The results show that the relative error of depth measurement is 1% compared with that of the high-precision SEM measurement results. Moreover， the proposed method can obtain a clear image of the TSV bottom with a high aspect ratio， and can also effectively enhance the wide spectrum interference signal and fringe contrast at the bottom compared with white light microscopic interference. Therefore， the proposed method can be used to precisely measure the 3D topography of TSVs with a higher aspect ratio. 无损测量 光学显微干涉 自适应像差补偿 COMSOL仿真 硅通孔形貌 nondestructive measurement optical microscopic interference adaptive aberration compensation COMSOL simulation through silicon via morphology 1 引 言 1 随着电子元器件向着高性能、低功耗、尺寸更小的方向不断发展，硅通孔（Through Silicon Via， TSV）工艺成为半导体器件三维封装的技术途径之一。TSV是指在芯片和芯片之间制作垂直导通，实现芯片之间的堆叠和互连的技术 ［ 1 1 - 2 2 1-2 ］ 。然而，半导体器件三维集成对TSV的尺度提出了更高的要求，而且半导体制造的高效率对TSV三维形貌的在线无损测量提出更快的要求。对于高深宽比TSV三维形貌的测量，非光学测量主要包括扫描电子显微术（SEM） ［ 3 3 ］ 、原子力显微术（AFM） ［ 4 4 ］ 等方法。这种方法具有破坏性，仅对陪片进行测量，且制样和SEM扫描拼接测量全过程耗时数小时，测量效率低，不适用于在线测量。光学测量主要包括白光显微干涉法 ［ 5 5 - 7 7 5-7 ］ 、共焦显微成像法 ［ 8 8 - 9 9 8-9 ］ 、光谱反射法 ［ 10 10 - 15 15 10-15 ］ 等，具有非接触、无破坏、测量效率高等优点。韩国科技大学的研究人员提出了一种结合低相干近红外干涉光谱技术、共聚焦技术和光学显微技术的新型混合光学探头，实现了直径为5.954 μm、深度为40.420 μm、深宽比约为6.8∶1的TSV样品的测量 ［ 9 9 ］ 。韩国标准科学研究院的研究小组采用近红外飞秒脉冲激光的光梳作为光源，基于光谱干涉法实现了深宽比为7∶1的TSV的深度测量 ［ 10 10 ］ 。共焦显微成像法测量的样品深宽比有限，光谱反射法无法直接获得待测样品的三维形貌分布。白光显微干涉法利用宽带光的低相干特性，与显微成像技术相结合，可一次性复原出待测样品的三维形貌 ［ 16 16 ］ 。国内中北大学的研究人员利用白光显微干涉法测量了宽度为22.6 μm、深度为89.78 μm的沟槽，深宽比为4∶1，但对于宽度小于5 μm的高深宽比沟槽，探测光则无法到达沟槽的底部 ［ 7 7 ］ 。白光显微干涉法的工作波长为可见光，对于TSV硅通孔材料不透明，经过干涉显微物镜的聚焦光束探测高深宽比TSV底部时，仅中央小数值孔径（Numerical Aperture， NA）的光束为有效探测光束，边缘NA的光束经过深孔的侧壁多次反射形成串扰光束。因此，韩国首尔大学的研究人员通过在照明光路中引入孔径光阑，限制物镜NA，排除多次反射光的串扰，得到了直径为4.27 μm、深度为47.9 μm、深宽比为11.2∶1的TSV的顶部直径、孔深度等参数，但底部形貌失真，横向分辨率低，仅为数微米 ［ 5 5 - 6 6 5-6 ］ 。 由于硅材料在近红外波段具有较高的透过率，采用近红外宽带光源穿透TSV的高深宽比结构，有望突破硅基材料对大NA探测光的遮挡，增大探测光的光通量，增强近红外宽带光干涉信号，用于提取TSV的形貌信息。然而，大NA光束会被TSV结构调制，降低光束聚焦性，产生像差，严重影响成像质量和干涉条纹的对比度。自适应光学是以成像波前像差校正为目标的波前控制技术，能实时补偿由光束传播、系统光学误差、成像环境和运动扰动等因素引入的波前畸变，从而改善系统成像性能，获取高分辨率目标图像，目前广泛应用于天文成像或者生物成像中的动态像差补偿 ［ 17 17 ］ 。本文引入变形镜自适应像差补偿技术，消除近红外光源测量TSV高深宽比结构所引入的调制像差，提高探测光的聚焦能力，增强TSV的底部成像和干涉信号。 为了实现对高深宽比TSV尺寸的高精度无损测量并实现TSV的三维形貌复原，本文提出了基于像差补偿的近红外显微干涉法。依据COMSOL Multiphysics有限元仿真软件得到的三维TSV高深宽比结构对探测光的调制像差规律，设置有待变形镜补偿的像差种类和量值，引入基于频域的评价函数指标阈值，判定TSV底部图像的聚焦状态，获得待测TSV清晰的底部像，提升探测光的重聚焦能力。在此基础上，利用构建的近红外显微干涉系统对多种深宽比的TSV样品进行垂直扫描干涉测量，解算出TSV的深度值并得到了其三维形貌分布。 2 原 理 1 2.1　基于像差补偿的近红外显微干涉光学测量系统 2 图1 图1 为Linnik型近红外反射式显微干涉测量系统的光路。整个光学系统主要由4部分组成：（1）近红外宽光谱低相干科勒照明系统，工作波长为1.25~1.45 μm；（2）由参考臂和测试臂构成的Linnik型显微干涉系统；（3）由变形镜（Deformable Mirror, DM）构成的光瞳像差主动补偿模块；（4）由压电陶瓷（PZT）和红外探测器CCD构成的干涉图同步扫描采集与处理系统。近红外短相干光源LS发出的光束经过由透镜L 1 ，L 2 ，L 3 组成的柯勒照明系统产生多视场均匀照明光后，经立方分光棱镜BS 1 分为测试光和参考光；打开光阑L，参考光经过平面反射镜M 2 ，M 3 和显微物镜MO 2 ，照射到平面反射镜M 4 上，并原路返回经立方分光棱镜BS 1 和透镜L 4 被红外探测器CCD接收；测试光经DM反射，依次经反射镜M 1 和显微物镜MO 1 ，照射放置于压电陶瓷PZT上的待测样品S；待测样品S被照明后，带有调制像差的反射光原路返回至变形镜DM，光阑L遮挡住参考臂后变形镜对像差进行补偿，补偿结束后，打开光阑L使得参考光在红外探测器CCD上与测试光形成干涉。 测量系统的工作波长为1.25~1.45 μm，它是硅基材料的透射窗口波段，该波段的大NA探测光可穿透硅基材料到达深孔底部，有利于提高高深宽比TSV深孔底部的探测光通量。现有的近红外显微物镜主要集中于生物医学领域，其波段一般是0.8~1.1 μm，不是硅基最有效的近红外透过窗口，长工作距、大NA近红外物镜的设计与研制关键技术没有先例可参考。该系统配备了自主研发的大NA（0.5~0.9）近红外宽光谱（1.25~1.45 μm）显微物镜模组MO 1 和MO 2 。另一方面，大NA近红外探测光入射高深宽比TSV深孔时，必然引起衍生的调制像差，从而严重影响探测光在深孔底部的能量聚焦。为此，该系统配备了由变形镜DM构成的光瞳像差主动补偿模块，用于补偿调制像差，提高探测光的聚焦能力，增强TSV的底部成像和干涉信号。 2.2　像差补偿与三维形貌提取原理 2 宽带光源发出经显微物镜会聚的大NA探测光束，穿过具有一定厚度的待测硅晶圆或硅基TSV样品后，必然在测试臂引入像差，影响宽带光干涉显微检测系统的成像性能和干涉信号对比度。 应用自适应光学或主动光学进行像差补偿，首要问题是能够获取和表征系统光瞳面上包含像差的光瞳函数，不失一般性，将出瞳处的波前相位看作是一个标准球面波和像差波前的叠加。标准球面波可用理想透镜的相位变换因子表示，即： <math display="block" id="M1"><mi>t</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced><mo>=</mo><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced open="[" close="]" separators="|"><mrow><mo>-</mo><mi mathvariant="normal">i</mi><mfrac><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi>f</mi></mrow></mfrac><mfenced separators="|"><mrow><msup><mrow><mi>x</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>y</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927375&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927373&type= 48.68333054 9.90600014 . 需要求解球面波经过样品调制后产生的像差 <math id="M2"><mi>φ</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927379&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927377&type= 12.02266598 4.82600021 。因而，出瞳处的波前复振幅表示为： <math display="block" id="M3"><msub><mrow><mi>g</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><mi>A</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced open="{" close="}" separators="|"><mrow><mi mathvariant="normal">i</mi><mfenced open="[" close="]" separators="|"><mrow><mi>t</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced><mo>+</mo><mi>φ</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></mrow></mfenced></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927383&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927382&type= 62.39933777 6.34999990 ， 式中 <math id="M4"><mi>A</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927389&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927387&type= 12.61533356 4.82600021 是出瞳处振幅分布。光瞳函数为： <math display="block" id="M5"><msub><mrow><mi>P</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><mi>A</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced open="[" close="]" separators="|"><mrow><mi mathvariant="normal">i</mi><mi>φ</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927393&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927391&type= 43.26466751 5.58799982 ， 式中： <math id="M6"><mi>φ</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927397&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927395&type= 12.02266598 4.82600021 是球面波经过样品调制后产生的光瞳像差。本文选用标准Zernike多项式来拟合像差，即 <math id="M7"><mi>φ</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>j</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>a</mi></mrow><mrow><mi>j</mi></mrow></msub></mrow><msub><mrow><mi>Z</mi></mrow><mrow><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927400&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927398&type= 27.85533524 8.97466660 ， <math id="M8"><msub><mrow><mi>a</mi></mrow><mrow><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927404&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927401&type= 2.62466669 3.72533321 是具体的Zernike项系数， <math id="M9"><msub><mrow><mi>Z</mi></mrow><mrow><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927410&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927407&type= 3.04800010 3.72533321 是Zernike项， N 为整数。 如 图1 图1 所示，DM用来补偿像差，该像差补偿系统由多个驱动器单元组成，采用机械变形法来匹配期望的共轭波前，每个驱动器对表面 <math id="M10"><msub><mrow><mi>Z</mi></mrow><mrow><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927416&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927412&type= 12.53066730 4.82600021 的作用效果即为变形镜的影响函数 ［ 18 18 - 19 19 18-19 ］ 。这些影响函数均为线性函数，通过施加合适的幅度 <math id="M11"><msub><mrow><mi>A</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927422&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927419&type= 3.80999994 3.72533321 至 M 个驱动器上产生期望的表面 <math id="M12"><mi>S</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927428&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927425&type= 11.93799973 4.82600021 。影响函数的线性方程组表达式如下： <math display="block" id="M13"><mi>S</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>M</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>A</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><msub><mrow><mi>Z</mi></mrow><mrow><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927434&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927431&type= 38.26933289 8.55133343 . 操控变形镜进行像差补偿就是将变形镜每个驱动器的影响函数 <math id="M14"><msub><mrow><mi>Z</mi></mrow><mrow><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927439&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927436&type= 12.53066730 4.82600021 用泽尼克多项式拟合，即可将需要的各项泽尼克像差系数转变为施加给变形镜各个驱动器的电压值，即： <math display="block" id="M15"><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>M</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>A</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><msub><mrow><mi>Z</mi></mrow><mrow><mi>i</mi></mrow></msub><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>j</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>a</mi></mrow><mrow><mi>j</mi></mrow></msub></mrow><msub><mrow><mi>Z</mi></mrow><mrow><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927445&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927442&type= 37.42266846 8.97466660 . 变形镜对光瞳像差补偿达到预期目标后，即可在CCD上得到干涉信号对比度增强的干涉条纹和清晰的底部成像。在此基础上，利用垂直扫描干涉法 ［ 20 20 ］ 对TSV深孔展开三维形貌测量。近红外光为宽带光源，其显微干涉条纹图是不同波长光的干涉叠加结果，光强分布为： <math display="block" id="M16"><mi>I</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>I</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><msubsup><mo>∫</mo><mrow><msub><mrow><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi>c</mi></mrow></msub><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi>b</mi></mrow></msub></mrow><mrow><msub><mrow><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi>c</mi></mrow></msub><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi>b</mi></mrow></msub></mrow></msubsup><mrow><mi>ψ</mi><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo></mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mfrac><mrow><mn mathvariant="normal">4</mn><mi mathvariant="normal">π</mi><mo stretchy="false">(</mo><mi>z</mi><mo>-</mo><msub><mrow><mi>z</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mi>λ</mi></mrow></mfrac><mi mathvariant="normal">d</mi><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927452&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927448&type= 70.01933289 8.63599968 ， 其中： <math id="M17"><msub><mrow><mi>I</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927458&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927455&type= 2.53999996 3.72533321 为背景光强， <math id="M18"><msub><mrow><mi>λ</mi></mrow><mrow><mi mathvariant="bold">c</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927463&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927461&type= 2.87866688 3.72533321 为光源中心波长， <math id="M19"><mn mathvariant="normal">2</mn><msub><mrow><mi>λ</mi></mrow><mrow><mi mathvariant="normal">b</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927469&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927466&type= 4.74133301 3.72533321 为光谱的带宽， <math id="M20"><mi>ψ</mi><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927475&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927472&type= 7.45066643 3.80999994 为干涉图在CCD上关于波长 <math id="M21"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927582&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927579&type= 1.77800000 2.79399991 的能量分布， <math id="M22"><mi>z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927488&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927484&type= 1.77800000 2.79399991 为PZT的相对位置， <math id="M23"><msub><mrow><mi>z</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927493&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927490&type= 2.96333337 3.72533321 为零光程差位置。因此，当像面像素点在零光程差位置即 <math id="M24"><mi>z</mi><mo>=</mo><msub><mrow><mi>z</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927499&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927495&type= 9.39799976 3.72533321 时，干涉强度最大。当PZT在垂直方向上移动时，测试臂与参考臂的光程差为零，像面上获取的干涉信号最大，形成一系列相干峰。此时，一系列相干峰位置对应于样品面的一系列相对高度，这些随视场点变化的相对高度信息，形成了样品表面的三维形貌。 2.3　TSV底部成像清晰度评价函数 2 变形镜对光瞳像差补偿时，由于补偿到位的图像有更尖锐清晰的边缘，相较于未补偿到位的模糊图像具有更多的高频分量，因此可使用基于频域的评价函数来判定TSV底部的成像状态，判断变形镜是否达到预期的补偿效果 ［ 21 21 ］ 。分析图像频域特性的常用手段是傅里叶变换，对于连续图像 <math id="M25"><mi>g</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927507&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927501&type= 11.76866722 3.97933316 ，当 <math id="M26"><mo>∫</mo><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mo>+</mo><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>g</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mi mathvariant="normal">d</mi><mi>x</mi><mi mathvariant="normal">d</mi><mi>y</mi><mo><</mo><mi mathvariant="normal">∞</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927513&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927510&type= 39.03133011 7.87400007 时，存在频域变换，即： <math display="block" id="M27"><mi>F</mi><mfenced separators="|"><mrow><mi>μ</mi><mo>,</mo><mi>ν</mi></mrow></mfenced><mo>=</mo><mo>∫</mo><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mo>+</mo><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><mi>g</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mo stretchy="false">[</mo><mo>-</mo><mi>j</mi><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><mo stretchy="false">(</mo><mi>μ</mi><mi>x</mi><mo>+</mo><mi>ν</mi><mi>y</mi><mo stretchy="false">)</mo><mo stretchy="false">]</mo><mi mathvariant="normal">d</mi><mi>x</mi><mi mathvariant="normal">d</mi><mi>y</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927519&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927516&type= 74.42199707 7.87400007 . 对于 <math id="M28"><mi>M</mi><mo>×</mo><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927526&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927523&type= 11.59933376 2.87866688 个像素的图像进行二维离散傅里叶变换，得到： <math display="block" id="M29"><mi>G</mi><mfenced separators="|"><mrow><mi>k</mi><mo>,</mo><mi>l</mi></mrow></mfenced><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>x</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>y</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></munderover></mstyle><mrow><mi>g</mi><mfenced separators="|"><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow></mfenced></mrow></mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced open="[" close="]" separators="|"><mrow><mo>-</mo><mi mathvariant="normal">j</mi><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><mfenced separators="|"><mrow><mfrac><mrow><mi>k</mi><mi>x</mi></mrow><mrow><mi>M</mi></mrow></mfrac><mo>+</mo><mfrac><mrow><mi>l</mi><mi>y</mi></mrow><mrow><mi>N</mi></mrow></mfrac></mrow></mfenced></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927532&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927529&type= 72.13600159 9.90600014 . 图像的功率谱函数表示为： <math display="block" id="M30"><mi>P</mi><mfenced separators="|"><mrow><mi>k</mi><mo>,</mo><mi>l</mi></mrow></mfenced><mo>=</mo><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>G</mi><mfenced separators="|"><mrow><mi>k</mi><mo>,</mo><mi>l</mi></mrow></mfenced></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927539&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927536&type= 31.07266808 5.41866684 . 二维离散傅里叶变换后的图像清晰度评价函数为： <math display="block" id="M31"><mi>f</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>M</mi><mi>N</mi></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>k</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></munderover></mstyle><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>l</mi><mo>=</mo><mn mathvariant="normal">0</mn></mrow><mrow><mi>N</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow></munderover></mstyle><mrow><mroot><mrow><msup><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>l</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot></mrow></mrow><mi>P</mi><mfenced separators="|"><mrow><mi>k</mi><mo>,</mo><mi>l</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927544&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927542&type= 51.13866425 8.63599968 . 3 仿 真 1 TSV对大NA探测光标准球面波引起的调制像差与TSV的孔直径、深度和边界条件等表征参数有关，随待测TSV样品表征参数的变化，补偿像差的量值和种类会相应发生变化。因此，为了让变形镜能高效补偿由TSV引起的调制像差，使近红外宽带光显微干涉系统检测TSV底部时能获取清晰的底部图像和明显的相干峰干涉信号，仿真研究了TSV对探测光产生的像差调制规律。 3.1　TSV模型建立 2 COMSOL Multiphysics仿真软件以有限元法为基础，通过求解偏微分方程（单物理场）或偏微分方程组（多物理场）来模拟真实的物理现象。本文仿真所用的功能组件为COMSOL的波动光学模块，该模块从麦克斯韦方程出发，将仿真区域划分为多个离散网格，通过定义边界条件设定入射场，以设置初始值，通过不断的迭代计算出近场的复振幅。如 图2 图2 所示，这里以直径为10 μm、深度为65 μm的TSV阵列为例，阐述仿真中TSV模型的建立方法（彩图见期刊电子版）。图中灰色区域为TSV，常规标准大气压下，硅折射率为3.5，蓝色区域为空气，折射率为1.0。红色线条为入射边界，绿色线条为反射电场的接收边界，模型周围设置10 μm的完美匹配层，用于吸收到达边界的电场。设置入射电场为平面波，定义光源的工作波长为1.325 μm，焦距为 <math id="M32"><mrow><mrow><mroot><mrow><mn mathvariant="normal">3</mn></mrow><mrow/></mroot></mrow><mo>/</mo><mrow><mn mathvariant="normal">10</mn></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927550&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927547&type= 10.15999985 4.57200003 μm，使探测光恰好聚焦到TSV的底部，且NA=0.5。设置最大网格单元尺寸为 <math id="M33"><mrow><mrow><mi>λ</mi></mrow><mo>/</mo><mrow><mn mathvariant="normal">10</mn></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927555&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927552&type= 6.60400009 3.30200005 ，最小网格单元尺寸为 <math id="M34"><mrow><mrow><mi>λ</mi></mrow><mo>/</mo><mrow><mn mathvariant="normal">20</mn></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927561&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927558&type= 6.60400009 3.30200005 ，构建网格。模型求解的自由度与网格划分粗细有关，所有模型求解的自由度数在 <math id="M35"><mn mathvariant="normal">1</mn><mo>×</mo><msup><mrow><mn mathvariant="normal">10</mn></mrow><mrow><mn mathvariant="normal">7</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927567&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927563&type= 11.51466751 3.04800010 数量级，达到网格解析波长的所需标准，在内存允许的情况下保证了计算的准确性。 3.2　仿真结果及分析 2 由于不能保证探测光的光轴与样品底部的中心轴完全重合，因此需要研究样品底部不同位置采样点的像差规律。本文仿真了探测光聚集在TSV底部 A ， B 两点的像差， A 点是TSV底部的中心点， B 点是TSV底部偏离中心4 μm的一个点。平面波从左侧边界入射，聚焦在TSV的底部位置并反射回去，电场传输过程如 图3 图3 所示。COMSOL仿真得到光场近场复振幅后，用近轴透镜构建显微物镜模型，模型中样品位于显微物镜的前焦面上，光瞳面位于显微物镜的后焦面上，构成具有傅里叶变换运算的光场传播光路模型，在MATLAB中计算光瞳函数。 图4 图4 和 图5 图5 分别为光瞳面的波面相位和标准泽尼克多项式拟合结果。可以看出，点的像差主要是系数为3.69 <math id="M36"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927582&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927579&type= 1.77800000 2.79399991 的离焦和系数为0.2 <math id="M37"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927603&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927601&type= 1.77800000 2.79399991 的球差， B 点的像差主要是系数为3.87 <math id="M38"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927603&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927601&type= 1.77800000 2.79399991 的离焦、系数为 - 2.5 <math id="M39"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927603&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927601&type= 1.77800000 2.79399991 的 x 方向彗差和系数为1.01 <math id="M40"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927603&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927601&type= 1.77800000 2.79399991 的球差。由此可知，TSV底部中心点（ A ）的像差主要类型为离焦和球差，TSV底部偏中心点（ B ）的像差主要类型为离焦、球差和彗差。 4 实验与结果分析 1 4.1　实验装置与补偿效果 2 根据 图1 图1 系统光路构建的Linnik型近红外显微干涉测量系统如 图6 图6 所示，实验使用NA为0.5、放大倍率为20×、工作波段在1.25~1.45 μm的显微物镜。整个测试过程分为三步：第一步，根据3.2小节的仿真结果，设置需变形镜补偿的像差种类和量值，获得待测TSV底部清晰的像和干涉条纹；第二步，通过PZT驱动待测样品进行垂直扫描，在红外探测器CCD上同步接收待测样品不同深度表面的干涉条纹图；第三步，采用垂直扫描干涉算法对干涉图进行处理，得到待测样品的深度和宽度测量结果。实验过程中，将宽度为10 μm、深度为65 μm的TSV样品放在PZT上方的测试平台上，调试光路使顶部清晰成像且存在高对比度干涉条纹，随后启动PZT对TSV孔进行垂直扫描成像。该实验以 <math id="M41"><msub><mrow><mi>λ</mi></mrow><mrow><mi mathvariant="normal">c</mi></mrow></msub><mo>/</mo><mn mathvariant="normal">8</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927636&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927633&type= 6.68866634 4.14866638 的步长对TSV样品进行垂直扫描，共采集TSV孔顶部和底部426幅图像，前213幅为顶部附近的干涉图像，后213幅为底部附近的干涉图像。 最终，采集到的未进行像差补偿的TSV孔顶底部的干涉图像如 图7 图7 所示，图中所示为间隔 <math id="M42"><msub><mrow><mi mathvariant="bold-italic">λ</mi></mrow><mrow><mi mathvariant="bold">c</mi></mrow></msub><mo>/</mo><mn mathvariant="normal">8</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927649&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927645&type= 6.85799980 4.14866638 步长的连续4幅图。样品顶部的干涉条纹清晰可见，但由于探测底部的入射光和反射光经过TSV高深宽比结构的调制，样品底部的干涉条纹变得微弱，难以分辨，孔底部图像也极为模糊。 为使样品底部图像清晰且增强底部的宽谱干涉信号和对比度，需要进行光瞳像差的补偿。实验用于像差补偿的装置是ALPAO变形镜，由69个驱动器单元组成，光瞳直径为10.5 mm，最大倾斜补偿量为60 μm，最大球差补偿量为40 μm。根据3.2节仿真所得的各项像差系数和占比，转换成相应的电压值驱动变形镜，变形镜将补偿电压分成20次对光瞳像差进行迭代补偿。利用2.3节所述的基于频域的评价函数作为图像清晰度评定指标，结合爬山搜索法确定最佳的像差补偿位置。如 图8 图8 所示，当像差补偿迭代到第12次时， 式（10） 式（10） 中的图像清晰度评价函数 <math id="M43"><mi>f</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927663&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=39927660&type= 2.45533323 3.64066648 最大，即此时TSV的底部像最清晰，如 图9 图9 所示。 补偿达到预期目标后，在测试光路中加入参考光，对样品进行垂直扫描干涉测量。TSV孔顶底部的干涉图如 图10 图10 所示，与补偿前的 图7 图7 比较可发现，样品的底部图像明显更清晰。比较补偿前后底部一个像素点的干涉信号，如 图11 图11 和 图12 图12 所示，补偿前混叠相干信号较多，不利于相干信号的识别和解调；而补偿后的相干信号包络变得清晰，干涉信号对比度得到有效增强。 4.2　近红外显微干涉与白光显微干涉测量对比 2 基于4.1节的测量结果，采用垂直扫描干涉算法对所获得的426幅顶底部干涉图像进行处理，得到待测TSV的三维形貌复原图和一个方向的一维轮廓曲线，分别如 图13 图13 （a）和 图13 图13 （b）所示。 1,2;3,4; 参照现有的国际标准［ISO 5436-1：2000（E）］和W/3准则，在TSV样品测量中，TSV底部的理论宽度为 W ，则高度测量结果的示值读取区域为距离TSV理论边缘位置 W /3以外的区域。重复测量10次并计算TSV深度的均值、标准差，以及均值相对SEM检测结果的绝对误差和相对误差，视场范围内其中一个孔（ 图13 图13 （a）方框）的10次测量结果如 表1 表1 所示。深度平均值为65.95 μm，顶部直径为9.60 μm，底部直径为9.20 μm，与 图13 图13 （c）的扫描电镜结果（65.27 μm）对比，深度测量相对误差为1%。从10次检测结果的标准差可以看出，本文方法具有良好的重复性和鲁棒性。 Number Result 1 65.97 μm 2 67.06 μm 3 66.64 μm 4 65.33 μm 5 66.44 μm 6 65.73 μm 7 65.78 μm 8 65.39 μm 9 66.07 μm 10 66.20 μm Mean value 65.95 μm Standard deviation 0.42 μm Absolute error 0.68 μm Relative error 1% 图13 图13 （d）是白光显微干涉法的测量结果，依次为TSV的一维轮廓、三维形貌和底部干涉条纹图。由于TSV的深宽比较大，白光显微干涉法的探测光难以到达样品底部，底部成像极为模糊且无干涉信号，10次白光显微干涉实验均无法测量出样品的三维形貌。 为进一步研究像差补偿近红外显微干涉法相对于白光显微干涉法测量TSV的优越性，本文分别用这两种方法测量了直径为20 μm、深度为77 μm、深宽比为3.85∶1和直径为10 μm、深度为103 μm、深宽比为10.3∶1的TSV样品。 图14 图14 （a）和 图14 图14 （b）分别是白光显微干涉法和像差补偿近红外显微干涉法测量深宽比为3.85∶1的TSV孔顶底部的干涉图，每种方法测量10次，深度平均值分别为76.62 μm和77.12 μm，但对于TSV的底部成像，后者明显比前者更加清晰，孔底部的直径平均值分别为11.79 μm和20.00 μm。因此，与白光显微干涉技术相比，本文方法可以在保证深度测量精度的情况下提高TSV的底部成像清晰度。对于深宽比为10.3∶1的TSV，本文方法也可准确测量出其深度值，如 图15 图15 所示，10次深度测量结果的平均值为103.85 μm，顶底部直径平均值分别为10.02 μm和10.72 μm。 1,2; 1,2; 5 结 论 1 本文提出了基于像差补偿的近红外显微干涉法，用于测量深宽比大于6∶1的TSV。采用能够穿透硅通孔的近红外宽带光作为光源，利用COMSOL Multiphysics有限元仿真软件研究了TSV高深宽比结构对探测光的调制规律，用于指导变形镜补偿的像差种类和量值，分析发现TSV底部中心点位置的像差主要为离焦和球差，底部偏离中心点位置的像差增加了彗差。随后，在实验中用变形镜自适应像差补偿模块对TSV引入的调制像差进行补偿，用基于频域的评价函数评定底部图像的清晰状态，获得待测TSV清晰的底部图像。最后在此基础上，使用垂直扫描干涉法得到待测TSV的深度与其三维形貌分布。实验测量了直径为10 μm、深度为65 μm、深宽比为6.5∶1和直径为10 μm、深度为103 μm、深宽比为10.3∶1两种TSV深孔，与高精度的SEM测量结果对比，深度测量的相对误差为1%。最后，与白光显微干涉技术进行对比，结果表明本文方法可以获得清晰的高深宽比TSV的底部图像，有效增强底部的宽谱干涉信号和对比度，能够准确测量更大高深宽比TSV的三维形貌。 参考文献 PAK J S ， RYU C ， KIM J . Electrical characterization of trough silicon via （TSV） depending on structural and material parameters based on 3D full wave simulation ［C］. 2007 International Conference on Electronic Materials and Packaging . 1922，2007 ， Daejeon ， Korea （South） . IEEE ， 2008 ： 1 - 6 . doi: 10.1109/emap.2007.4510331 http://dx.doi.org/10.1109/emap.2007.4510331 Electrical characterization of trough silicon via （TSV） depending on structural and material parameters based on 3D full wave simulation KIM H ， CHO J ， KIM M ， et al . Measurement and analysis of a high-speed TSV channel ［J］. IEEE Transactions on Components， Packaging and Manufacturing Technology ， 2012 ， 2 （ 10 ）： 1672 - 1685 . doi: 10.1109/tcpmt.2012.2207900 http://dx.doi.org/10.1109/tcpmt.2012.2207900 Measurement and analysis of a high-speed TSV channel 李智 ， 王向军 . MEMS中几何量的测试方法 ［J］. 微细加工技术 ， 2003 （ 1 ）： 51 - 56 . doi: 10.3321/j.issn:1004-924X.2003.01.007 http://dx.doi.org/10.3321/j.issn:1004-924X.2003.01.007 MEMS中几何量的测试方法 LI ZH ， WANG X J . Measurement methods of geometric characteristics in MEMS ［J］. Microfabrication Technology ， 2003 （ 1 ）： 51 - 56 . （in Chinese） . doi: 10.3321/j.issn:1004-924X.2003.01.007 http://dx.doi.org/10.3321/j.issn:1004-924X.2003.01.007 Measurement methods of geometric characteristics in MEMS GARCı́A R ， PÉREZ R . Dynamic atomic force microscopy methods ［J］. Surface Science Reports ， 2002 ， 47 （ 6/7/8 ）： 197 - 301 . doi: 10.1016/s0167-5729(02)00077-8 http://dx.doi.org/10.1016/s0167-5729(02)00077-8 Dynamic atomic force microscopy methods JO T ， KIM S ， PAHK H . 3D measurement of TSVs using low numerical aperture white-light scanning interferometry ［J］. Journal of the Optical Society of Korea ， 2013 ， 17 （ 4 ）： 317 - 322 . doi: 10.3807/josk.2013.17.4.317 http://dx.doi.org/10.3807/josk.2013.17.4.317 3D measurement of TSVs using low numerical aperture white-light scanning interferometry HYUN C ， KIM S ， PAHK H . Methods to measure the critical dimension of the bottoms of through-silicon vias using white-light scanning interferometry ［J］. Journal of the Optical Society of Korea ， 2014 ， 18 （ 5 ）： 531 - 537 . doi: 10.3807/josk.2014.18.5.531 http://dx.doi.org/10.3807/josk.2014.18.5.531 Methods to measure the critical dimension of the bottoms of through-silicon vias using white-light scanning interferometry 薛晨阳 ， 孔繁华 ， 张文栋 ， 等 . 利用白光干涉原理测量MEMS深槽结构 ［J］. 传感技术学报 ， 2006 ， 19 （ 5 ）： 1516 - 1518， 1522 . doi: 10.3969/j.issn.1004-1699.2006.05.054 http://dx.doi.org/10.3969/j.issn.1004-1699.2006.05.054 利用白光干涉原理测量MEMS深槽结构 XUE CH Y ， KONG F H ， ZHANG W D ， et al . Measuring MEMS deep-trench structures with white light interference ［J］. Chinese Journal of Sensors and Actuators ， 2006 ， 19 （ 5 ）： 1516 - 1518， 1522 . （in Chinese） . doi: 10.3969/j.issn.1004-1699.2006.05.054 http://dx.doi.org/10.3969/j.issn.1004-1699.2006.05.054 Measuring MEMS deep-trench structures with white light interference AHN H ， PARK J ， KIM J A ， et al . Optical fiber-based confocal and interferometric system for measuring the depth and diameter of through silicon vias ［J］. Journal of Lightwave Technology ， 2016 ， 34 （ 23 ）： 5462 - 5466 . doi: 10.1109/jlt.2016.2618419 http://dx.doi.org/10.1109/jlt.2016.2618419 Optical fiber-based confocal and interferometric system for measuring the depth and diameter of through silicon vias AHN H ， BAE J ， PARK J ， et al . A hybrid non-destructive measuring method of three-dimensional profile of through silicon vias for realization of smart devices ［J］. Scientific Reports ， 2018 ， 8 ： 15342 . doi: 10.1038/s41598-018-33728-w http://dx.doi.org/10.1038/s41598-018-33728-w A hybrid non-destructive measuring method of three-dimensional profile of through silicon vias for realization of smart devices JIN J ， KIM J W ， KANG C S ， et al . Precision depth measurement of through silicon vias （TSVs） on 3D semiconductor packaging process ［J］. Optics Express ， 2012 ， 20 （ 5 ）： 5011 - 5016 . doi: 10.1364/oe.20.005011 http://dx.doi.org/10.1364/oe.20.005011 Precision depth measurement of through silicon vias （TSVs） on 3D semiconductor packaging process KU Y S ， HUANG K C ， HSU W . Characterization of high density through silicon vias with spectral reflectometry ［J］. Optics Express ， 2011 ， 19 （ 7 ）： 5993 - 6006 . doi: 10.1364/oe.19.005993 http://dx.doi.org/10.1364/oe.19.005993 Characterization of high density through silicon vias with spectral reflectometry KU Y S ， YANG F S . Reflectometer-based metrology for high-aspect ratio via measurement ［J］. Optics Express ， 2010 ， 18 （ 7 ）： 7269 - 7280 . doi: 10.1364/oe.18.007269 http://dx.doi.org/10.1364/oe.18.007269 Reflectometer-based metrology for high-aspect ratio via measurement 杨楠卓 . 基于光谱反射原理的高深宽比沟槽深度测量技术研究 ［D］. 长春 ： 长春理工大学 ， 2020 . 基于光谱反射原理的高深宽比沟槽深度测量技术研究 YANG N ZH . Research on Depth Measurement of the High Aspect Ratio Trenches Based on Spectral Reflectometry ［D］. Changchun ： Changchun University of Science and Technology ， 2020 . （in Chinese） Research on Depth Measurement of the High Aspect Ratio Trenches Based on Spectral Reflectometry 张传维 . 基于模型的微纳深沟槽结构红外反射谱测量方法研究 ［D］. 武汉 ： 华中科技大学 ， 2009 . 基于模型的微纳深沟槽结构红外反射谱测量方法研究 ZHANG CH W . Research on Method of Model-based Infrared Reflectometry for Micro/nano Deep Trench Structures ［D］. Wuhan ： Huazhong University of Science and Technology ， 2009 . （in Chinese） Research on Method of Model-based Infrared Reflectometry for Micro/nano Deep Trench Structures 高志山 ， 袁群 ， 孙一峰 ， 等 . 微结构显微光学无损检测方法 ［J］. 光子学报 ， ， 2022 ， 51 （ 8 ）： 17 - 34 . doi: 10.3788/gzxb20225108.0851501 http://dx.doi.org/10.3788/gzxb20225108.0851501 微结构显微光学无损检测方法 GAO ZH SH ， YUAN Q ， SUN Y F ， et al . Non-destructive test methods of microstructures by optical microscopy ［J］. Acta Photonica Sinica ， 2022 ， 51 （ 8 ）： 17 - 34 . （in Chinese） . doi: 10.3788/gzxb20225108.0851501 http://dx.doi.org/10.3788/gzxb20225108.0851501 Non-destructive test methods of microstructures by optical microscopy 马龙 . 白光扫描干涉测量方法与系统的研究 ［D］. 天津 ： 天津大学 ， 2011 . doi: 10.3969/j.issn.1672-6030.2011.01.008 http://dx.doi.org/10.3969/j.issn.1672-6030.2011.01.008 白光扫描干涉测量方法与系统的研究 MA L . Study on White Light Scanning Interferometry ： Measurement Method and System ［D］. Tianjin ： Tianjin University ， 2011 . （in Chinese） . doi: 10.3969/j.issn.1672-6030.2011.01.008 http://dx.doi.org/10.3969/j.issn.1672-6030.2011.01.008 Study on White Light Scanning Interferometry： Measurement Method and System 王玉坤 ， 曹召良 ， 李大禹 ， 等 . 液晶-变形镜自适应光学系统的数据采集与处理软件设计 ［J］. 光学 精密工程 ， 2018 ， 26 （ 6 ）： 1507 - 1516 . doi: 10.3788/ope.20182606.1507 http://dx.doi.org/10.3788/ope.20182606.1507 液晶-变形镜自适应光学系统的数据采集与处理软件设计 WANG Y K ， CAO ZH L ， LI D Y ， et al . Design of liquid crystal-deformable mirror adaptive optical system data acquisition and process software ［J］. Opt. Precision Eng. ， 2018 ， 26 （ 6 ）： 1507 - 1516 . （in Chinese） . doi: 10.3788/ope.20182606.1507 http://dx.doi.org/10.3788/ope.20182606.1507 Design of liquid crystal-deformable mirror adaptive optical system data acquisition and process software 刘维 ， 徐珺楠 ， 金玳冉 ， 等 . 基于CNN-SPGD算法的相干光通信像差校正方法研究 ［J］. 光学 精密工程 ， 2022 ， 30 （ 6 ）： 743 - 754 . doi: 10.37188/OPE.20223006.0743 http://dx.doi.org/10.37188/OPE.20223006.0743 基于CNN-SPGD算法的相干光通信像差校正方法研究 LIU W ， XU J N ， JIN D R ， et al . Research on aberration correction method of coherent optical communication based on CNN-SPGD algorithm ［J］. Opt. Precision Eng. ， 2022 ， 30 （ 6 ）： 743 - 754 . （in Chinese） . doi: 10.37188/OPE.20223006.0743 http://dx.doi.org/10.37188/OPE.20223006.0743 Research on aberration correction method of coherent optical communication based on CNN-SPGD algorithm HIDDLESTON H R ， LYMAN D D ， SCHAFER E L . Comparisons of deformable-mirror models and influence functions ［C］. San Diego ， ' 91 . Proc SPIE 1542， Active and Adaptive Optical Systems ， DiegoSan ， CA ， USA . 1991， 1542 ： 20 - 33 . doi: 10.1117/12.48791 http://dx.doi.org/10.1117/12.48791 Comparisons of deformable-mirror models and influence functions 孙一峰 ， 高志山 ， 范筱昕 ， 等 . 低相干干涉高深宽比结构高度检测信号处理技术 ［J］. 光子学报 ， 2022 ， 51 （ 9 ）： 331 - 341 . doi: 10.3788/gzxb20225109.0912001 http://dx.doi.org/10.3788/gzxb20225109.0912001 低相干干涉高深宽比结构高度检测信号处理技术 SUN Y F ， GAO ZH SH ， FAN X X ， et al . Signal processing for height measurement of high-aspect-ratio structure based on low-coherence interferometry ［J］. Acta Photonica Sinica ， 2022 ， 51 （ 9 ）： 331 - 341 . （in Chinese） . doi: 10.3788/gzxb20225109.0912001 http://dx.doi.org/10.3788/gzxb20225109.0912001 Signal processing for height measurement of high-aspect-ratio structure based on low-coherence interferometry 王昕 . 含噪声图像的多聚焦融合算法 ［J］. 光学 精密工程 ， 2011 ， 19 （ 12 ）： 2977 - 2984 . doi: 10.3788/ope.20111912.2977 http://dx.doi.org/10.3788/ope.20111912.2977 含噪声图像的多聚焦融合算法 WANG X . Multi-focus fusion algorithm for noisy images ［J］. Optics and Precision Engineering ， 2011 ， 19 （ 12 ）： 2977 - 2984 . （in Chinese） . doi: 10.3788/ope.20111912.2977 http://dx.doi.org/10.3788/ope.20111912.2977 Multi-focus fusion algorithm for noisy images 10.37188/OPE.20233103.0301 O436.1 TH744 A 1004-924X（2023）03-0301-12 吴春霞,马剑秋,高志山等.基于像差补偿的近红外显微干涉硅通孔测量[J].光学精密工程,2023,31(03):301-312. 吴春霞,马剑秋,高志山等.基于像差补偿的近红外显微干涉硅通孔测量[J].光学精密工程,2023,31(03):301-312. WU Chunxia,MA Jianqiu,GAO Zhishan,et al.Measurement of through silicon via by near-infrared micro interferometry based on aberration compensation[J].Optics and Precision Engineering,2023,31(03):301-312. WU Chunxia,MA Jianqiu,GAO Zhishan,et al.Measurement of through silicon via by near-infrared micro interferometry based on aberration compensation[J].Optics and Precision Engineering,2023,31(03):301-312. 国家重点研发计划资助项目（No.2019YFB2005500）；国家自然科学基金资助项目（No.62175107，No.U1931120）；江苏省六大人才高峰项目（No.RJFW-019）；中国科学院光学系统先进制造技术重点实验室基金资助项目（No.KLOMT190201）；上海在线检测与控制技术重点实验室基金资助项目（No.ZX2021102） 2022-09-08 2022-10-11 2023-02-10 2023-02-15 光学精密工程 3 31