详情

原文传递一种全聚焦成像的声束校准方法

专利名称：	一种全聚焦成像的声束校准方法
摘要：	本发明是一种全聚焦成像的声束校准方法，该方法通过实际测量标准试块中横通孔的幅值信息来建立声束校准矩阵以补偿声波能量衰减，校准的步骤如下：制作标准试块；全矩阵数据的采集及存储；子全矩阵数据的提取；对每一个子全矩阵数据进行全聚焦成像；提取横通孔幅值信息并生成原始校准数据；对原始校准数据做线性插值计算生成声束校准矩阵；对全聚焦成像进行声束校准。本发明方法可有效减小缺陷定量检测误差和检测漏检率，改善全聚焦图像能量均匀性。
专利类型：	发明专利
国家地区组织代码：	北京;11
申请人：	中国航发北京航空材料研究院
发明人：	徐娜;许路路;王东升;梁菁;何方成
专利状态：	有效
申请号：	CN201910010408.X
公开号：	CN109490419A
代理机构：	中国航空专利中心 11008
代理人：	陈宏林
分类号：	G01N29/06(2006.01)I;G;G01;G01N;G01N29
申请人地址：	100095 北京市海淀区北京市81号信箱科技发展部
主权项：	1.一种全聚焦成像的声束校准方法，其特征在于：该方法的步骤如下：步骤一、制作标准试块该标准试块(4)由实际被检测材料制成，具有长方体外形，长度不小于所采用超声阵列探头长度的2倍，高度不小于拟实际检测的最大深度，宽度不小于超声阵列探头宽度，在标准试块(4)长度方向正中心位置、高度方向自上而下等间距加工一组沿宽度方向的横通孔(3)，横通孔(3)的直径为拟实际检测定量评价时的标准缺陷尺寸，最上端横通孔中心距标准试块上表面位置为拟实际检测近表面分辨力的1/2，最下端横通孔中心距标准试块下表面位置为拟实际检测近底面分辨力的1/2，相邻的横通孔(3)的中心间距为横通孔直径的3‑5倍之间；步骤二、全矩阵数据采集及存储将包含N个阵元晶片的超声阵列探头(1)放置在标准试块(4)上表面的中心位置上，N不小于64，在声束校准过程中不移动超声阵列探头(1)的位置，设置相控阵超声探伤仪，从超声阵列探头(1)中第一号阵元晶片起依次使每一个阵元晶片发射超声波，同时，每个阵元晶片发射超声波时，依次采集超声阵列探头(1)中每一个阵元晶片接收到的超声回波信号，并将超声回波数据进行存储，作为一个全矩阵数据，记为Sij，其中，i＝1,2,…,N；j＝1,2,…,N，其中，i代表发射阵元晶片的序号，j代表接收阵元晶片的序号；步骤三、子全矩阵数据提取依次提取全矩阵数据Sij中从第k个阵元晶片起的连续M个阵元晶片发射和接收的超声回波信号，M为实际检测时所使用阵元晶片个数，M不大于32个，k＝1,2,…,N‑M+1，获得一系列子全矩阵数据Vij，其中，i＝k,k+1,…,k+M‑1；j＝k,k+1,…,k+M‑1，定义每一个子全矩阵数据中所包含的阵元晶片为一个子阵列，对提取的每一个子全矩阵数据Vij进行希尔伯特变换，得到一系列子全矩阵包络数据Gij；步骤四、对每一个子全矩阵数据进行全聚焦成像，其过程如下：建立二维直角坐标系oxz，其中，坐标原点o(5)位于超声阵列探头的中心位置，x轴(7)与标准试块(4)的长度方向一致，z轴(8)与标准试块(4)的高度方向一致，将整个标准试块(4)作为成像区域，沿x轴(7)和z轴(8)的方向进行网格划分，网格划分间隔不大于0.2mm，网格划分形成的交点定义为聚焦点(9)，数量计为W个，再依次定义每一个子阵列的中心位置为子阵列原点ok(6)；依次根据每个子全矩阵包络数据Gij中对应的发射、接收阵元晶片到聚焦点(9)的距离计算声波传播时间，再在每个聚焦点位置将信号幅值进行叠加，即可通过公式Ⅰ和公式Ⅱ计算出在子全矩阵包络数据Gij下的任意聚焦点(x,z)的聚焦点幅值Ip(x，z，k)其中，(xi,0)和(xj,0)分别为发射和接收阵元晶片的坐标值，c为超声波在标准试块中的传播速度；将一个子全矩阵包络数据下所计算出的整个成像区域中所有聚焦点幅值的集合记为Fk，则：Fk＝{Ip1(x，z，k),Ip2(x，z，k),…,Ipn(x，z，k),…,IpW(x，z，k)}，至此，根据每个子全矩阵包络数据Gij便获得一系列全聚焦成像数据Fk；步骤五、提取横通孔(3)幅值信息并生成原始校准数据根据所获得的全聚焦成像数据生成一系列全聚焦图像，依次在每一个全聚焦图像中提取出全部横通孔(3)的最大横通孔幅值的信息及相对于子矩阵原点ok的相对位置信息，获得全部横通孔幅值数据D＝{dmn}和相应位置信息(Xmn,Zmn)，其中，m代表一系列全聚焦图像的数量，n代表每个全聚焦图像中横通孔的个数，横通孔幅值数据共计m*n个，同时，根据相应位置信息获得横通孔幅值数据D沿x轴和z轴的位置范围(X1,X2)和(Z1,Z2)，取出横通孔幅值数据D中的最大值Dmax，并将Dmax分别除以横通孔幅值数据D生成原始校准数据D1＝{Dmax/dmn}；步骤六、原始校准数据线性插值计算首先，在x轴的位置范围(X1,X2)和z轴的位置范围(Z1,Z2)内，沿x轴和z轴进行网格划分，网格划分间隔不大于0.2mm，网格划分形成的交点定义为新聚焦点，然后，根据原始校准数据D1，通过线性插值运算得到在新聚焦点的校准数据，即为声束校准矩阵C＝{c1,c2,…cq}，q代表一系列新聚焦点的个数；步骤七、对全聚焦成像进行声束校准：将包含M个阵元晶片的超声阵列探头放置在被检测工件的上表面，按照步骤二采集全矩阵数据并存储，获得一组全矩阵数据Vhl(h＝1,2,…,M；l＝1,2,…,M)，并对采集的全矩阵数据进行希尔伯特变换，获得一组全矩阵包络数据Ghl；以超声阵列探头中心位置作为坐标原点，确定该全矩阵数据Vhl中全部阵元晶片所对应的坐标位置，取该全矩阵数据Vhl中全部发射、接收阵元晶片到每个新聚焦点的距离，计算声波传播时间，再在每个新聚焦点位置将信号幅值进行叠加并乘以对应的声束校准值cp，cp为声束校准矩阵C中的一个，即可计算出每个新聚焦点声束校准后的新聚焦点幅值I′p(x，z)其中，tij(x,z)由公式Ⅱ计算得出，c为超声波在被检测工件中的传播速度，至此，根据公式Ⅲ便获得在位置范围(X1,X2)和(Z1,Z2)内全部新聚焦点声束校准后的幅值，即获得了声束校准后的全聚焦成像结果。
所属类别：	发明专利