CN114354502B - 基于激光声表面波的硅片加工表面损伤和残余应力表征方法 - Google Patents

基于激光声表面波的硅片加工表面损伤和残余应力表征方法 Download PDF

Info

Publication number
CN114354502B
CN114354502B CN202111574343.5A CN202111574343A CN114354502B CN 114354502 B CN114354502 B CN 114354502B CN 202111574343 A CN202111574343 A CN 202111574343A CN 114354502 B CN114354502 B CN 114354502B
Authority
CN
China
Prior art keywords
acoustic wave
surface acoustic
silicon wafer
residual stress
layer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111574343.5A
Other languages
English (en)
Other versions
CN114354502A (zh
Inventor
林滨
刘再蔚
梁小虎
杜安尧
马小康
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN202111574343.5A priority Critical patent/CN114354502B/zh
Publication of CN114354502A publication Critical patent/CN114354502A/zh
Application granted granted Critical
Publication of CN114354502B publication Critical patent/CN114354502B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明属于超声检测的技术领域,具体涉及基于激光声表面波的硅片加工表面损伤和残余应力表征方法,包括建立正演模型,基于正演模型构建声表面波频散数据的反演模型和数据误差函数,进行激光超声实验激发复频声表面波,进行频散分析得到相应的频散数据,基于实验频散数据及构建的反演模型计算得到所测区域的亚表面损伤及残余应力。本发明可准确量化硅片加工表面的亚表面损伤程度及残余应力大小,并且,其测量结果具有更高的检测精度与可信度。

Description

基于激光声表面波的硅片加工表面损伤和残余应力表征方法
技术领域
本发明属于超声检测的技术领域,具体涉及基于激光声表面波的硅片加工表面损伤和残余应力表征方法。
背景技术
单晶硅作为最受欢迎的半导体器件制造材料,现已广泛应用于电子、太阳能光伏技术和红外光学等领域。在硅片的制造中,由于磨削具有高材料去除率和低成本等特点,通常可用作精磨之前的粗加工工艺。在磨削过程中,磨粒与晶圆表面的相互作用会不可避免地导致亚表面损伤和残余应力。虽然这些缺陷可以通过蚀刻或化学机械抛光完全消除,但这些过程既费时又费钱。因此,通过控制磨削工艺参数来尽可能减少这些缺陷是具有研究意义的。相应地,也需要找到一种可靠的检测方法来评估硅片磨削过程造成的亚表面损伤和残余应力。
目前,用于检测亚表面损伤的方法可以分为两类:破坏性方法和非破坏性方法。其中,破坏性方法包括角度抛光法、横截面技术、透射电子显微镜检测和化学蚀刻法等。虽然这些破坏性方法可以提供可观察、可靠的信息,但是它们过于耗时。而非破坏性方法则使用光学或声学特性来检测,例如,基于白光干涉仪的表面粗糙度预测法、扫描声学显微镜检测法等。
然而,现有的检测方法获得的硅片的残余应力值是不准确的,并且,大多方法都忽略了亚表面损伤层的影响。为此,亟需提出一种新型的技术方案以解决上述问题。
发明内容
本发明的目的在于:针对现有技术的不足,提供基于激光声表面波的硅片加工表面损伤和残余应力表征方法,其能够将亚表面损伤程度与残余应力解耦,从而实现硅片中亚表面损伤和残余应力的同时测量。
为了实现上述目的,本发明采用如下技术方案:
基于激光声表面波的硅片加工表面损伤和残余应力表征方法,包括以下步骤:
S1、根据硅片加工表面的特性,建立声表面波在含有表层残余应力的分层各向异性介质表面传播的正演模型;
S2、根据正演模型,构建与声表面实验频散数据相关的误差函数;
S3、利用脉冲激光在硅片加工表面激发声表面波,并在预设的位置采用声表面波探头采集激发出来的声表面波垂直位移信号;
S4、使脉冲激光对应于硅片加工表面的区域为测量区域,将硅片以测量区域为中心旋转一定的角度;
S5、重复S3和S4,直至得到预设数量的声表面波信号,将所有采集到的声表面波信号进行频散分析,得到一系列声表面波实验频散数据;
S6、将S5的实验频散数据代入到S2的误差函数中并通过计算优化得到满足最小化误差函数的正演函数相关模型参数,即反演得到层介质的弹性常数和残余应力,以表征损伤程度和残余应力。
进一步地,所述S1中的正演模型至少基于以下三个声表面波边界条件获得:
条件(i):界面处的位移连续:
在(1)式中,和ui分别表示声表面波在层和基体介质中的位移,其中,i=1,2,3…n;
和ui的一般形式分别可以表示为:
在(2)式中,x1,x3αi和l3分别为待定参数,c为声表面波速度,t表示时间;
条件(ii):界面处的边界力连续:
在(3)式中,和fi分别表示声表面波在层和基体介质中的边界力,其中,i=1,2,3…n;
和fi可以由公式(4)-(7)计算得到:
其中,和sij分别表示层和基体介质中的增量应力分量,/>为层介质的增量弹性常数,/>为层介质的有效弹性常数,Cijkl为基体介质的弹性常数,/>为层介质的膨胀系数,δij表示克罗内克尔函数, Sij为层介质中的残余应力;
条件(iii):自由表面的边界力消失:
基于公式(1)-(8)可以计算出待定参数αi和l3及声表面波速度c,从而得到声表面波的正演模型,即模型参数与声表面理论频散数据之间的关系:
在(9)式中,和Cijkl分别表示层和基体的弹性常数,Sij表示层中的残余应力,h为层厚度,θ为声表面波的传播方向,/>表示正演, f为频率,/>为一系列沿着不同角度的声表面波理论频散数据。
进一步地,所述S2的误差函数为:
在(10)式中,N为实验频散数据的总数,N=∑iNi,Ni表示不同方向上的声表面波频散数据量,cmeas(f)为实验频散数据,ctheo(f)为理论频散数据,misfit的函数值的大小体现了实验频散数据与理论频散数据的匹配程度。
进一步地,所述S5还包括:将硅片测量区域的反演结果转化对应为杨氏模量分布,以杨氏模量变化量体现硅片的亚表面损伤程度。
本发明的有益效果在于:1)本发明通过基于激光声表面波的硅片加工表面损伤和残余应力表征方法,其能够将亚表面损伤程度与残余应力解耦,从而实现硅片中亚表面损伤和残余应力的同时测量,综合考虑了亚表面损伤程度和残余应力对声表面波速度的影响;2)该方法可以很好地实现硅片的亚表面损伤程度及残余应力的量化;3) 由于各向异性的假设,该方法并不局限于磨削硅片,还可以用于绝大数材料加工表面的检测。
附图说明
图1为本发明的流程图。
图2为本发明的激光超声实验的示意图。
图3为本发明的通过检测得到的一系列声表面波信号图。
图4为本发明获得的磨削硅片亚表面损伤层的示意图。
图5为本发明的声表面波实验频散数据及反演结果所对应的理论频散曲线。
图6为本发明的硅片测量区域的杨氏模量分布图。
具体实施方式
如在说明书及权利要求当中使用了某些词汇来指称特定组件,本领域技术人员应可理解,制造商可能会用不同名词来称呼同一个组件。本说明书及权利要求并不以名称的差异来作为区分组件的方式,而是以组件在功能上的差异来作为区分的准则。如在通篇说明书及权利要求当中所提及的“包含”为一开放式用语,故应解释成“包含但不限定于”。“大致”是指在可接受的误差范围内,本领域技术人员能够在一定误差范围内解决技术问题,基本达到技术效果。
在本发明的描述中,需要理解的是,术语“上”、“下”、“前”、“后”、“左”、“右”、水平”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
在发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
在现有的检测亚表面损伤层中残余应力的方法中,当依靠硅片的变形及Stoney方程时,可以基于两个假设:(1)假设硅片足够薄以使其产生可观的变形,且残余应力层足够薄;(2)残余应力在整个硅片上均匀分布。然而,由于以上的假设,不仅使得最终获得的残余应力值不准确的,并且,其还完全忽略了亚表面损伤层的影响。
为了解决上述问题,以下结合附图1~6和具体实施例对本发明作进一步详细说明,但不作为对本发明的限定。
基于激光声表面波的硅片加工表面损伤和残余应力表征方法,包括以下步骤:
S1、根据硅片加工表面的特性,在线弹性框架下建立声表面波在含有表层残余应力的分层各向异性介质表面传播的正演模型,正演模型至少基于以下三个声表面波边界条件获得:
条件(i):界面处的位移连续:
在(1)式中,和ui分别表示声表面波在层和基体介质中的位移,其中,i=1,2,3…n;
和ui的一般形式分别可以表示为:
在(2)式中,x1,x3αi和l3分别为待定参数,c为声表面波速度,t表示时间;
条件(ii):界面处的边界力连续:
在(3)式中,和fi分别表示声表面波在层和基体介质中的边界力,其中,i=1,2,3…n;
和fi可以由公式(4)-(7)计算得到:
其中,和sij分别表示层和基体介质中的增量应力分量,/>为层介质的增量弹性常数,/>为层介质的有效弹性常数,Cijkl为基体介质的弹性常数,/>为层介质的膨胀系数,δij表示克罗内克尔函数, Sij为层介质中的残余应力;
条件(iii):自由表面的边界力消失:
基于公式(1)-(8)可以计算出待定参数αi和l3及声表面波速度c,从而得到声表面波的正演模型,即模型参数与声表面理论频散数据之间的关系:
在(9)式中,和Cijkl分别表示层和基体的弹性常数,Sij表示层中的残余应力,h为层厚度,θ为声表面波的传播方向,/>表示正演, f为频率,/>为一系列沿着不同角度的声表面波理论频散数据;
S2、根据正演模型,构建与声表面实验频散数据相关的误差函数,误差函数为:
在(10)式中,N为实验频散数据的总数,N=∑iNi,Ni表示不同方向上的声表面波频散数据量,cmeas(f)为实验频散数据,ctheo(f)为理论频散数据,misfit的函数值的大小体现了实验频散数据与理论频散数据的匹配程度;
S3、利用脉冲激光在硅片加工表面激发声表面波,并在预设的位置或在一段距离外采用声表面波探头采集激发出来的声表面波垂直位移信号;
S4、使脉冲激光对应于硅片加工表面的区域为测量区域,将硅片以测量区域为中心旋转一定的角度;
S5、重复S3和S4,直至得到预设数量的声表面波信号,将所有采集到的声表面波信号进行频散分析,得到一系列声表面波实验频散数据;
S6、将S5的实验频散数据代入到S2的误差函数中并通过计算优化得到满足最小化误差函数的正演函数相关模型参数,即反演得到层介质的弹性常数和残余应力,以表征损伤程度和残余应力。
其中,在S3步骤中,采用如图2所示的实验配置,对直径为 200mm、厚度为800μm的磨削硅片表面进行多角度激光超声实验,即利用脉冲激光在硅片加工表面沿不同晶向激发复频声表面波、并在不同位置处拾取其垂直位移信号。同时,在操作过程中,在每个角度下收集两次声表面波信号,两次激发源相距10mm,从而得到如图3 所示的多个声表面波垂直位移信号。
假设亚表面层为正交各向异性介质且残余应力为等双轴残余应力,此时的正演模型含有10个待定参数,即9个正交各向异性常数和一个残余应力值;对损伤层厚度的初步了解可以得到更加准确的弹性常数及残余应力。通过截面抛光法,如图4所示,表面损伤层厚度可以确定为8μm。
在S5步骤中,对图2的实验所得的声表面波信号进行频散分析,即对垂直位移信号进行频散分析,得到对应的频散数据,如图5所示,将35°和145°方向的频散数据进行平均以对应上述的亚表面层为正交各向异性假设。接着,按照图1中的步骤进行反演,反演结果如下表一所示,同时,该结果所对应的理论频散曲线可参见图5。
表一、硅片测量区域的反演结果
此外,为了更好地理解表一中正交各向异性常数的性质,可以将其转化为平面杨氏模量分布,即将硅片测量区域的反演结果转化对应为杨氏模量分布。
如图6所示,典型角度下的杨氏模量已经被列在如下的表二中;可以看到,沿着<100>晶向的杨氏模量变化量小于沿着<110>晶向的杨氏模量变化量,此外,沿着磨削方向杨氏模量变化量小于垂直于磨削方向杨氏模量变化量。杨氏模量变化量可以看作是亚表面损伤程度的体现。
表二、硅片测量区域反演结果所对应的杨氏模量分布
可见,上述基于激光声表面波技术的硅片加工表面/亚表面损伤和残余应力综合表征方法,在线弹性框架下建立声表面波在含有表面应力的分层各向异性介质表面传播的正演模型;基于正演模型构建声表面波频散数据的反演模型;利用脉冲激光在硅片加工表面沿不同晶向激发复频声表面波,并在不同位置处拾取其垂直位移信号;对垂直位移信号进行频散分析得到相应的频散数据;基于实验频散数据及构建的反演模型得到所测区域的亚表面损伤及残余应力。该方法不仅可准确量化硅片加工表面的亚表面损伤程度及残余应力大小,并且,其测量结果具有更高的检测精度与可信度。
根据上述说明书的揭示和教导,本发明所属领域的技术人员还能够对上述实施方式进行变更和修改。因此,本发明并不局限于上述的具体实施方式,凡是本领域技术人员在本发明的基础上所作出的任何显而易见的改进、替换或变型均属于本发明的保护范围。此外,尽管本说明书中使用了一些特定的术语,但这些术语只是为了方便说明,并不对本发明构成任何限制。

Claims (3)

1.基于激光声表面波的硅片加工表面损伤和残余应力表征方法,其特征在于,包括以下步骤:
S1、根据硅片加工表面的特性,建立声表面波在含有表层残余应力的分层各向异性介质表面传播的正演模型,
所述正演模型至少基于以下三个声表面波边界条件获得:
条件(i):界面处的位移连续:
在(1)式中,和ui分别表示声表面波在层和基体介质中的位移,其中,i=1,2,3…n;
和ui的一般形式分别可以表示为:
在(2)式中,x1,x3αi和l3分别为待定参数,c为声表面波速度,t表示时间;
条件(ii):界面处的边界力连续:
在(3)式中,和fi分别表示声表面波在层和基体介质中的边界力,其中,i=1,2,3…n;
和fi可以由公式(4)-(7)计算得到:
其中,和sij分别表示层和基体介质中的增量应力分量,/>为层介质的增量弹性常数,/>为层介质的有效弹性常数,Cijkl为基体介质的弹性常数,/>为层介质的膨胀系数,δij表示克罗内克尔函数,Sij为层介质中的残余应力;
条件(iii):自由表面的边界力消失:
基于公式(1)-(8)可以计算出待定参数αi和l3及声表面波速度c,从而得到声表面波的正演模型,即模型参数与声表面理论频散数据之间的关系:
在(9)式中,和Cijkl分别表示层和基体的弹性常数,Sij表示层中的残余应力,h为层厚度,θ为声表面波的传播方向,f为频率,/>为一系列沿着不同角度的声表面波理论频散数据;
S2、根据正演模型,构建与声表面实验频散数据相关的误差函数;
S3、利用脉冲激光在硅片加工表面激发声表面波,并在预设的位置采用声表面波探头采集激发出来的声表面波垂直位移信号;
S4、使脉冲激光对应于硅片加工表面的区域为测量区域,将硅片以测量区域为中心旋转一定的角度;
S5、重复S3和S4,直至得到预设数量的声表面波信号,将所有采集到的声表面波信号进行频散分析,得到一系列声表面波实验频散数据;
S6、将S5的实验频散数据代入到S2的误差函数中并通过计算优化得到满足最小化误差函数的正演函数相关模型参数,即反演得到层介质的弹性常数和残余应力,以表征损伤程度和残余应力。
2.如权利要求1所述的基于激光声表面波的硅片加工表面损伤和残余应力表征方法,其特征在于,所述S2的误差函数为:
在(10)式中,N为实验频散数据的总数,N=∑iNi,Ni表示不同方向上的声表面波频散数据量,cmeas(f)为实验频散数据,ctheo(f)为理论频散数据,misfit的函数值的大小体现了实验频散数据与理论频散数据的匹配程度。
3.如权利要求1所述的基于激光声表面波的硅片加工表面损伤和残余应力表征方法,其特征在于,所述S5还包括:将硅片测量区域的反演结果转化对应为杨氏模量分布,以杨氏模量变化量体现硅片的亚表面损伤程度。
CN202111574343.5A 2021-12-21 2021-12-21 基于激光声表面波的硅片加工表面损伤和残余应力表征方法 Active CN114354502B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111574343.5A CN114354502B (zh) 2021-12-21 2021-12-21 基于激光声表面波的硅片加工表面损伤和残余应力表征方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111574343.5A CN114354502B (zh) 2021-12-21 2021-12-21 基于激光声表面波的硅片加工表面损伤和残余应力表征方法

Publications (2)

Publication Number Publication Date
CN114354502A CN114354502A (zh) 2022-04-15
CN114354502B true CN114354502B (zh) 2024-01-30

Family

ID=81101069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111574343.5A Active CN114354502B (zh) 2021-12-21 2021-12-21 基于激光声表面波的硅片加工表面损伤和残余应力表征方法

Country Status (1)

Country Link
CN (1) CN114354502B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116593400B (zh) * 2023-07-17 2023-10-17 国家电投集团江西电力有限公司 一种太阳能板黑斑损伤检测方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6543288B1 (en) * 1998-11-04 2003-04-08 National Research Council Of Canada Laser-ultrasonic measurement of elastic properties of a thin sheet and of tension applied thereon
JP2004077460A (ja) * 2002-06-17 2004-03-11 Nippon Steel Corp 残留応力分布測定装置及び残留応力分布測定方法
CN103674359A (zh) * 2013-12-13 2014-03-26 烟台富润实业有限公司 一种复合材料残余应力的激光超声无损检测方法及设备
CN108168747A (zh) * 2018-02-11 2018-06-15 浙江大学 一种基于激光超声的工件表面残余应力测量装置及其方法
CN109521090A (zh) * 2018-10-18 2019-03-26 天津大学 一种激光无损表征薄膜杨氏模量的优化方法
CN110057912A (zh) * 2019-03-11 2019-07-26 天津大学 基于声表面波信号处理获得弱频散条件下频散曲线的方法
CN112326786A (zh) * 2020-10-16 2021-02-05 河北工业大学 基于电磁超声Lamb波S1模态群速度的金属板应力检测方法
CN113218875A (zh) * 2021-04-30 2021-08-06 东南大学 一种金属增材制造件残余应力的激光超声测量方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6543288B1 (en) * 1998-11-04 2003-04-08 National Research Council Of Canada Laser-ultrasonic measurement of elastic properties of a thin sheet and of tension applied thereon
JP2004077460A (ja) * 2002-06-17 2004-03-11 Nippon Steel Corp 残留応力分布測定装置及び残留応力分布測定方法
CN103674359A (zh) * 2013-12-13 2014-03-26 烟台富润实业有限公司 一种复合材料残余应力的激光超声无损检测方法及设备
CN108168747A (zh) * 2018-02-11 2018-06-15 浙江大学 一种基于激光超声的工件表面残余应力测量装置及其方法
CN109521090A (zh) * 2018-10-18 2019-03-26 天津大学 一种激光无损表征薄膜杨氏模量的优化方法
CN110057912A (zh) * 2019-03-11 2019-07-26 天津大学 基于声表面波信号处理获得弱频散条件下频散曲线的方法
CN112326786A (zh) * 2020-10-16 2021-02-05 河北工业大学 基于电磁超声Lamb波S1模态群速度的金属板应力检测方法
CN113218875A (zh) * 2021-04-30 2021-08-06 东南大学 一种金属增材制造件残余应力的激光超声测量方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Inversion of surface damage and residual stress in ground silicon wafers by laser surface acoustic wave technology;Lin, Bin;ULTRASONICS;全文 *
激光超声检测技术的材料表面微小损伤检测;姬冠妮;王亚亚;史二娜;激光杂志(第009期);全文 *
激光超声表面波用于加工表面损伤检测的传播机理研究;刘雪莲;机械工程;全文 *

Also Published As

Publication number Publication date
CN114354502A (zh) 2022-04-15

Similar Documents

Publication Publication Date Title
Kushibiki et al. Development of the line-focus-beam ultrasonic material characterization system
Dallas et al. Resonance ultrasonic vibrations for crack detection in photovoltaic silicon wafers
AU2008226491B2 (en) A method and apparatus for in-line quality control of wafers
CN102414807B (zh) 确定基底完整性的非破坏性信号传播系统和方法
CN114354502B (zh) 基于激光声表面波的硅片加工表面损伤和残余应力表征方法
Liu et al. Quantifying the subsurface damage and residual stress in ground silicon wafer using laser ultrasonic technology: A Bayesian approach
CN101093212A (zh) 检测由半导体材料构成的锭件中机械缺陷的方法和设备
CN110349877B (zh) 检测晶圆键合强度的方法及晶圆键合机台
Belyaev et al. Resonance ultrasonic vibration diagnostics of elastic stress in full-size silicon wafers
CN115420203A (zh) 半球谐振子金属化镀膜均匀性的评估方法
CN110132805B (zh) 一种2219铝合金铸锭平均晶粒尺寸超声波评定方法
CN113533519B (zh) 一种非接触无损评估材料各向异性的方法和装置
US4741212A (en) Method for determining structural defects in semiconductor wafers by ultrasonic microscopy
US20050039534A1 (en) Methods for measuring strength of film and determining quality of object having the film
US6413789B2 (en) Method of detecting and monitoring stresses in a semiconductor wafer
Degertekin et al. Hertzian contact transducers for nondestructive evaluation
Nakaso et al. Measurement of acoustic reflection coefficients by an ultrasonic microspectrometer
CN110865123A (zh) 一种计算激光声表面波频散曲线的方法
KR101678872B1 (ko) 실리콘 웨이퍼 중에 존재하는 원자 공공 농도의 정량 평가 방법, 실리콘 웨이퍼의 제조 방법, 및 이 제조 방법에 의해 제조된 실리콘 웨이퍼
Liang et al. Study on the relationship between the damage of machined surface layer and the dispersion curve of laser-induced surface acoustic wave
Kushibiki et al. Application of line-focus-beam acoustic microscope to inhomogeneity detection on SAW device materials
CN111289620B (zh) 各向异性材料的弹性常数检测方法及系统
Byelyayev Stress diagnostics and crack detection in full-size silicon wafers using resonance ultrasonic vibrations
CN106153733B (zh) 高温下测固有频率的试样夹持装置及其使用方法
Schneider et al. Nondestructive testing of damage layers in semiconductor materials by surface acoustic waves

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant