发明名称 一种半导体制造过程的多性能预测方法
摘要 本发明涉及一种半导体制造过程的多性能预测方法。首先选取半导体生产线上表征其状态的在制品水平参数、设备参数和工件参数为性能指标的影响因素;采集生产线相关数据,利用主成分分析法进行预处理,剔除冗余信息,并利用贝叶斯神经网络构建多性能预测模型,引入贝叶斯方法控制网络模型的复杂度;然后通过模型精度验证法分析模型性能是否符合性能评价准则,并对不符合标准的网络模型结构进行在线修正;最后确定影响工件平均加工周期和设备利用率的关键因素。本发明改善了半导体领域性能预测模型受约束条件限制、泛化性能较差等缺点,解决了半导体领域单一性能预测模型不能适用于多性能预测建模问题,是对半导体制造过程多性能预测方法的改进。
申请公布号 CN103745273A 申请公布日期 2014.04.23
申请号 CN201410005327.8 申请日期 2014.01.06
申请人 北京化工大学 发明人 曹政才;刘雪莲;刘民;李博;王炅;邱明辉
分类号 G06Q10/04(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 北京思海天达知识产权代理有限公司 11203 代理人 张慧
主权项 1.一种半导体制造过程的多性能预测方法,其特征在于包括以下步骤:步骤1,确定需预测的性能指标及半导体制造过程中影响性能指标的因素;需预测的性能指标包括工件平均加工周期和设备利用率;影响性能指标的因素包括在制品水平、设备参数和工件参数;步骤2,从半导体生产线采集步骤1中所述预测性能指标及影响因素的相关数据,然后对采集的样本数据进行预处理;预处理方法如下:步骤2.1,采用主成分分析方法剔除样本数据中的冗余信息,降低输入数据维数;用向量x=(x<sub>1</sub>,x<sub>2</sub>,...,x<sub>n</sub>)<sup>T</sup>表示性能指标影响因素,n为向量x中变量的个数;若x=(x<sub>1</sub>,x<sub>2</sub>,...,x<sub>n</sub>)<sup>T</sup>的协方差阵的n个特征根为λ<sub>1</sub>,λ<sub>2</sub>,...,λ<sub>n</sub>,λ<sub>1</sub>≥λ<sub>2</sub>≥...≥λ<sub>n</sub>>0,计算特征根对应性能指标影响参数的累计贡献率d:<maths num="0001"><![CDATA[<math><mrow><mi>d</mi><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&lambda;</mi><mi>j</mi></msub></mrow></mfrac></mrow></math>]]></maths>累计贡献率达到95%时的m即为剔除冗余信息后数据的维数;求特征根λ<sub>i</sub>对应的m个正交化单位特征向量a<sub>i</sub>=(a<sub>i1</sub>,a<sub>i2</sub>,...,a<sub>in</sub>),由u<sub>i</sub>=a<sub>i</sub>·x得到x降维后的m个主元u<sub>i</sub>,i=1,2,...,m,u<sub>i</sub>表达式为:u<sub>1</sub>=a<sub>11</sub>x<sub>1</sub>+a<sub>12</sub>x<sub>2</sub>+...+a<sub>1n</sub>x<sub>n</sub>u<sub>2</sub>=a<sub>21</sub>x<sub>1</sub>+a<sub>22</sub>x<sub>2</sub>+...+a<sub>2n</sub>x<sub>n</sub>............................................u<sub>m=</sub>a<sub>m1</sub>x<sub>1</sub>+a<sub>m2</sub>x<sub>2</sub>+...+a<sub>mn</sub>x<sub>n</sub>步骤2.2,利用“最大最小法”对筛选后的数据进行归一化;步骤3,采用贝叶斯神经网络方法建立半导体制造过程多性能预测模型;步骤3.1,将步骤1所确定的影响预测性能指标的参数作为贝叶斯神经网络的输入,预测性能指标作为输出;步骤3.2,选择网络结构,初始化网络参数;网络模型采用包括输入层I、隐含层H和输出层O的I-H-O三层结构;由输入层节点数I和输出层节点数O,根据隐含层神经元个数经验公式<img file="FDA0000453626120000012.GIF" wi="246" he="82" />计算隐含层节点数H,然后在H值周围用实验法确定最优的隐层节点数;设权值先验概率P(w|α,H)服从高斯分布,w为网络权值,α为控制网络的复杂程度的超参数,完成网络参数的初始化;步骤3.3,利用步骤2所得数据对贝叶斯神经网络进行训练;步骤4,对步骤3建立的预测模型进行在线修正;步骤4.1,按照步骤2的方法获取输入、输出实时数据并完成预处理;步骤4.2,通过模型精度验证法分析模型性能是否符合性能评价准则;模型预测精度为:<maths num="0002"><![CDATA[<math><mrow><mi>pre</mi><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>|</mo><mfrac><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>y</mi><mi>i</mi><mo>'</mo></msubsup></mrow><msub><mi>y</mi><mi>i</mi></msub></mfrac><mo>|</mo></mrow></math>]]></maths>如果模型预测精度达到了预测精度期望值,表明模型符合性能评价准则,不需要进行修正;否则,不符合性能评价准则,需要对贝叶斯神经网络模型结构进行修正,修正过程同步骤3;步骤5,应用网络权值分析法确定影响工件平均加工周期和设备利用率的关键因素;步骤5.1,获取步骤4建立的预测模型的权值,对每一个输入计算其对输出的贡献值;输入x<sub>i</sub>对输出的贡献值c<sub>i</sub>:<maths num="0003"><![CDATA[<math><mrow><msub><mi>c</mi><mi>i</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>h</mi></munderover><msub><mi>c</mi><mi>ij</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>h</mi></munderover><mfrac><msub><mi>w</mi><mi>ij</mi></msub><mrow><msub><mi>w</mi><mrow><mn>1</mn><mi>j</mi></mrow></msub><mo>+</mo><msub><mi>w</mi><mrow><mn>2</mn><mi>j</mi></mrow></msub><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>w</mi><mi>mj</mi></msub></mrow></mfrac></mrow></math>]]></maths>式中,i=1,2,...,m,h为网络隐层节点数,c<sub>ij</sub>为输入x<sub>i</sub>对隐含层第j个节点的贡献值,w<sub>ij</sub>为网络权值;步骤5.2,计算重要影响因素阈值;选取贡献值的平均值为阈值:<maths num="0004"><![CDATA[<math><mrow><mi>t</mi><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>c</mi><mi>i</mi></msub></mrow><mi>m</mi></mfrac></mrow></math>]]></maths>式中,t为阈值;步骤5.3,确定对输出平均加工周期和设备利用率影响较大的重要因素;将步骤5.1得到的每个输入的贡献值与步骤5.2得到的阈值比较,贡献值大于阈值的输入为工件平均加工周期和设备利用率的重要影响因素。
地址 100029 北京市朝阳区北三环东路15号