CN113030243A - 面向成像的地上钢结构件腐蚀缺陷漏磁检测方法及系统 - Google Patents

面向成像的地上钢结构件腐蚀缺陷漏磁检测方法及系统 Download PDF

Info

Publication number
CN113030243A
CN113030243A CN202110301694.2A CN202110301694A CN113030243A CN 113030243 A CN113030243 A CN 113030243A CN 202110301694 A CN202110301694 A CN 202110301694A CN 113030243 A CN113030243 A CN 113030243A
Authority
CN
China
Prior art keywords
leakage detection
magnetic
steel structural
flux leakage
structural member
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.)
Granted
Application number
CN202110301694.2A
Other languages
English (en)
Other versions
CN113030243B (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.)
Tsinghua University
Wuxi Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Tsinghua University
Wuxi Power Supply Co of State Grid Jiangsu Electric Power Co Ltd
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 Tsinghua University, Wuxi Power Supply Co of State Grid Jiangsu Electric Power Co Ltd filed Critical Tsinghua University
Priority to CN202110301694.2A priority Critical patent/CN113030243B/zh
Publication of CN113030243A publication Critical patent/CN113030243A/zh
Application granted granted Critical
Publication of CN113030243B publication Critical patent/CN113030243B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/72Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
    • G01N27/82Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws
    • G01N27/83Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws by investigating stray magnetic fields
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)

Abstract

面向成像的地上钢结构件腐蚀缺陷漏磁检测方法及系统,通过高精度漏磁检测装置调节被测地上钢结构件中的磁化方向,采集钢结构件在不同磁化方向下的漏磁检测信号,再根据漏磁检测信号的不同类型,对经直流滤波处理和归一化处理后的漏磁检测信号进行信号合成和边缘提取处理,得到面向反演成像所需的钢结构件腐蚀缺陷的轮廓参数,便于腐蚀缺陷的高精度成像。通过调节两对磁化线圈中的电流大小实现对被测地上钢结构件中的磁化方向进行控制,结合对在不同磁化方向下得到的漏磁检测信号的综合处理,最终实现对腐蚀缺陷的高精度成像检测。该方法具有装置结构简单、操作方便、缺陷成像检测精度高等优点,有利于腐蚀缺陷的高精度成像。

Description

面向成像的地上钢结构件腐蚀缺陷漏磁检测方法及系统
技术领域
本发明涉及无损检测领域,更具体地,涉及面向成像的地上钢结构件腐蚀缺陷漏磁检测方法及系统。
背景技术
漏磁检测由于其具有检测速度快、对检测环境要求低等特点已被广泛应用与地上钢结构件腐蚀缺陷检测中。通过采集地上钢结构件的漏磁检测信号,对检测到的漏磁信号进行分析处理,实现腐蚀缺陷的反演成像及评估是地上钢结构检测漏磁检测的最终目的,高精度的腐蚀缺陷反演结果可为地上钢结构的安全评估与维护提供有效的指导依据,因此,地上钢结构件腐蚀缺陷的高精度漏磁成像检测非常重要。
在现有的相关技术中,多数漏磁检测方法仅能实现对被测钢结构件单一方向的磁化,极大地限制了对缺陷的检测效果,如发明专利(201110350081.4)公开了一种储罐底板检测装置,采用两个永磁体组提供单一方向的磁化场,实现对储罐底板的缺陷检测,但是这个单一方向的磁化检测方法,仅能对单一走向的缺陷有较好的检测效果,对于其它方向走向的缺陷或缺陷轮廓不规则的缺陷,检测效果不佳。发明专利(201610614526.8)公开了一种基于旋转磁化场的钢板漏磁检测装置,基于初级线圈和次级线圈分别采用U形铁芯和铜棒为被测钢板提供不同方向的磁化,但是这种方法需要通过初级线圈中的磁化电路和次级线圈中的感应电流相互配合,装置结构设计复杂,且磁化场的方向无法任意调节,仅能对垂直于磁化场分布的裂纹缺陷有较为灵敏的检测效果,而对于没有明显方向性的腐蚀缺陷检测效果不佳。
发明内容
为解决现有技术中存在的不足,本发明的目的在于,提供面向成像的地上钢结构件腐蚀缺陷漏磁检测,通过高精度漏磁检测装置调节被测地上钢结构件中的磁化方向,采集钢结构件在不同磁化方向下的漏磁检测信号,再根据漏磁检测信号的不同类型,对经直流滤波处理和归一化处理后的漏磁检测信号进行信号合成和边缘提取处理,得到面向反演成像所需的钢结构件腐蚀缺陷的轮廓参数,便于腐蚀缺陷的高精度成像。采用的高精度漏磁检测装置结构简单、漏磁检测操作方便、腐蚀缺陷检测精度高,有利于腐蚀缺陷的高精度成像。
本发明采用如下的技术方案。
腐蚀缺陷漏磁检测方法的步骤如下:
步骤1,通过成对控制高精度漏磁检测装置中磁化线圈中的施加电流,实现对被测地上钢结构件中的磁化方向的任意调节,采集地上钢结构件在不同磁化方向下的漏磁检测信号;
步骤2,对步骤1采集的漏磁检测信号先直流滤波处理再归一化处理;
步骤3:根据漏磁检测信号的不同类型,对步骤2处理后的漏磁检测信号进行信号合成;
步骤4:对步骤3的合成信号进行边缘提取处理,得到地上钢结构件腐蚀缺陷的轮廓参数,实现缺陷高精度成像;
优选地,步骤1包括:
步骤1.1,高精度漏磁检测装置水平放置在待检测的地上钢结构件上;
步骤1.2,成对控制磁化线圈内施加的电流大小,实现磁轭中心的磁场方向沿逆时针方向旋转扫描;
步骤1.3,在步骤1.2进行的同时,由磁传感器阵列采集地上钢结构件腐蚀缺陷处的漏磁检测信号;
优选地,在步骤1.1中,
高精度漏磁检测装置的励磁部分包括:第一铁芯、第二铁芯、第三铁芯、第四铁芯、第一磁化线圈、第二磁化线圈、第三磁化线圈、第四磁化线圈、磁轭、磁传感器阵列;
其中,磁传感器阵列水平布置在地上钢结构件上;第一铁芯、第二铁芯、第三铁芯、第四铁芯垂直于地上钢结构件,并按照十字形顶点位置呈逆时针方向依次分布在磁传感器阵列四周;缠绕在第一铁芯上的第一磁化线圈、缠绕在第二铁芯上的第二磁化线圈、缠绕在第三铁芯上的第三磁化线圈和缠绕在第四铁芯上的第四磁化线圈,缠绕匝数相同、缠绕方向也相同;铁芯上方放置磁轭,磁传感器阵列位于磁轭正下方,并且磁传感器阵列中心与磁轭中心的连线垂直于地上钢结构件;
地上钢结构件、第一铁芯、第二铁芯、第三铁芯、第四铁芯和磁轭组成高精度漏磁检测装置的闭合磁回路;
磁轭是轴对称结构的,且至少有两条对称轴是相互垂直的;
优选地,对于磁轭的形状可以是但不限制于十字形、圆心、正方形、八边形,对于磁轭的材料可以是但不限制于软铁板、A3钢板、软磁合金;
优选地,在步骤1.2中,
四组磁化线圈中每两组为一对,第二磁化线圈和第四磁化线圈为第一对磁化线圈,第一磁化线圈和第三磁化线圈为第二对磁化线圈;
第一对磁化线圈和第二对磁化线圈施加的电流大小分别为A1和A2,则A1和A2及磁轭中心处的磁场偏转角度θ满足如下关系:
Figure BDA0002986512640000031
式中,
A1为第一对磁化线圈的施加电流,A2为第二对磁化线圈的施加电流;
θ为磁轭中心的磁场方向的偏转角度,该偏转角度以仅对第二对磁化线圈施加电流时磁轭中心的磁场方向为基准,按照逆时针方向旋转一周计算,即从θ为0°沿着逆时针方向旋转至360°;
当A1为零且A2不为零时,磁轭中心的磁场方向与第二对磁化线圈的磁化方向平行时,θ为0°;当A1不为零且A2为零时,磁轭中心的磁场方向与第一对磁化线圈的磁化方向平行时,θ为90°;
当A1和A2大小相同、方向相反时,第一对磁化线圈的磁场方向与第二对磁化线圈的磁场方向相互垂直,并且K为恒定值,因此磁轭中心的磁场大小保持不变,并且磁场偏转角度θ为45°或135°,此时磁轭中心处的磁场沿固定方向不变;
因此,通过控制两对磁化线圈中的施加电流,实现磁轭中心的磁场偏转方向沿逆时针依次偏转0°、45°、90°和135°四个扫描角度;
实际检测过程中,根据地上钢结构件的检测需求优化磁场的扫描方向及路径;考虑漏磁场的对称特性,仅将磁轭中心的磁场方向沿逆时针扫描半周,即仅从θ为0°沿着逆时针方向旋转至180°;还考虑到检测数据存储容量的限制,仅扫描θ为0°、45°、90°和135°有限个数的磁场方向;
优选地,在步骤2中,
对于有腐蚀缺陷的地上钢结构件,对其漏磁检测信号矩阵的直流滤波处理具体包括:
(1)在相同环境条件下,采用相同的高精度漏磁检测装置采集无腐蚀缺陷的地上钢结构件的漏磁检测信号,得到无腐蚀缺陷的地上钢结构件的漏磁检测信号矩阵B0
(2)在相同的磁场方向扫描条件下,从有腐蚀缺陷的地上钢结构件的漏磁检测信号Bi中减去无腐蚀缺陷的地上钢结构件的漏磁检测信号
Figure BDA0002986512640000032
即实现了直流滤波处理;
对于有腐蚀缺陷的地上钢结构件,将其直流滤波处理后的漏磁检测信号矩阵转换到[-1,1]范围内,实现归一化处理,满足如下关系式:
Figure BDA0002986512640000041
式中,
Norm(·)表示一种归一化处理函数;
Bi表示在不同磁场方向扫描时采集到的有腐蚀缺陷的地上钢结构件的漏磁检测信号,Bi=[B1,B2,...,Bn],i=1,2,...,n,定义n为不同磁场方向扫描时所记录的漏磁信号数量;
Figure BDA0002986512640000042
表示在不同磁场方向扫描时采集到的无腐蚀缺陷的地上钢结构件的漏磁检测信号,
Figure BDA0002986512640000043
Figure BDA0002986512640000044
表示有腐蚀缺陷的地上钢结构件的经直流滤波处理和归一化处理后的在不同磁场方向扫描时采集到的漏磁检测信号,
Figure BDA0002986512640000045
优选地,在步骤3中,
漏磁检测信号包括漏磁检测信号的水平分量信号、漏磁检测信号的垂直分量信号或漏磁检测信号法向分量信号中的至少一种;根据漏磁检测信号的不同类型,对步骤2处理后的漏磁检测信号进行信号合成,具体如下:
(1)对于漏磁检测信号的水平分量信号,信号合成矩阵满足如下关系式:
Figure BDA0002986512640000046
式中,
Bh_combine表示步骤2处理后的漏磁检测信号的水平分量信号的合成信号矩阵,
Figure BDA0002986512640000047
表示在不同磁场方向扫描时经步骤2处理后的漏磁检测信号的水平分量信号矩阵,
Figure BDA0002986512640000048
Ghi表示
Figure BDA0002986512640000049
所对应的处理模块;
*表示卷积运算;
abs(·)表示对矩阵数据取绝对值的函数;
(2)对于漏磁检测信号的垂直分量信号,信号合成矩阵满足如下关系式:
Figure BDA00029865126400000410
式中,
Bv_combine表示步骤2处理后的漏磁检测信号的垂直分量信号的合成信号矩阵,
Figure BDA0002986512640000051
表示在不同磁场方向扫描时经步骤2处理后的漏磁检测信号的垂直分量信号矩阵,
Figure BDA0002986512640000052
Gvi表示
Figure BDA0002986512640000053
所对应的处理模块;
优选地,
Ghi、Gvi与磁轭中心的磁场偏转角度θ有关,具体满足如下关系式:
当θ∈[0°,22.5°)∪[157.5°,202.5°)∪[337.5°,360°)时,
Figure BDA0002986512640000054
当θ∈[22.5°,67.5°)∪[202.5°,247.5°)时,
Figure BDA0002986512640000055
当θ∈[67.5°,112.5°)∪[247.5°,292.5°)时,
Figure BDA0002986512640000056
当θ∈[112.5°,157.5°)∪[292.5°,337.5°)时,
Figure BDA0002986512640000057
(3)对于漏磁检测信号的法向分量信号,信号合成矩阵满足如下关系式:
Figure BDA0002986512640000058
式中,
Bf_combine表示步骤2处理后的漏磁检测信号的法向分量信号的合成信号矩阵,
Figure BDA0002986512640000059
表示在不同磁场方向扫描时经步骤2处理后的漏磁检测信号的法向分量信号矩阵,
Figure BDA00029865126400000510
ki表示
Figure BDA00029865126400000511
的合成系数,满足0≤ki≤1;
面向成像的地上钢结构件腐蚀缺陷漏磁检测系统,包括漏磁检测信号采集模块、漏磁检测信号综合处理模块、漏磁检测信号反演成像模块;
漏磁检测信号采集模块,用于在不同磁化方向下采集漏磁检测信号;该模块内置高精度漏磁检测装置,通过成对控制高精度漏磁检测装置中磁化线圈中的施加电流,实现对被测地上钢结构件中的磁化方向的任意调节;
漏磁检测信号综合处理模块,用于对采集到的漏磁检测信号依次进行直流滤波处理、归一化处理、信号合成处理以及边缘提取处理;
漏磁检测信号反演成像模块,用于显示地上钢结构件腐蚀缺陷的三维图像;
优选地,漏磁检测信号综合处理模块包括直流滤波处理单元、归一化处理单元、水平分量和垂直分量信号合成处理单元、法向分量信号合成处理单元、边缘提取处理单元。
本发明的有益效果在于,与现有技术相比,通过调节两对磁化线圈中的电流大小实现对被测地上钢结构件中的磁化方向进行控制,结合对在不同磁化方向下得到的漏磁检测信号的综合处理,最终实现对腐蚀缺陷的高精度成像检测。该方法具有装置结构简单、操作方便、缺陷成像检测精度高等优点,有利于腐蚀缺陷的高精度成像。
附图说明
图1是本发明的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法的流程图;
图2是本发明的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法的高精度漏磁检测装置励磁部分的结构示意图;
图3是本发明的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法的磁轭中心处磁场偏转角度θ的位置示意图;
图4(a)、(b)、(c)、(d)分别是本发明的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法在磁轭中心的磁场偏转角度依次为0°、45°、90°和135°时,采集的地上钢结构件腐蚀缺陷处的漏磁检测信号的水平分量经过直流滤波处理和归一化处理后的二维图像;
图5(a)、(b)、(c)、(d)分别是本发明的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法在磁轭中心的磁场偏转角度依次为0°、45°、90°和135°时,采集的地上钢结构件腐蚀缺陷处的漏磁检测信号的水平分量经过直流滤波处理和归一化处理后,再经过卷积运算后得到的二维图像;
图6是本发明的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法得到的漏磁检测信号的水平分量信号的合成信号矩阵的二维图像;
图7(a)、(b)、(c)、(d)分别是本发明的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法在磁轭中心的磁场偏转角度依次为0°、45°、90°和135°时,采集的地上钢结构件腐蚀缺陷处的漏磁检测信号的法向分量经过直流滤波处理和归一化处理后的二维图像;
图8是本发明的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法得到的漏磁检测信号的水平分量信号的合成信号矩阵的二维图像。
具体实施方式
下面结合附图对本申请作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本申请的保护范围。
如图1,面向成像的地上钢结构件腐蚀缺陷漏磁检测方法的步骤如下:
步骤1,通过成对控制高精度漏磁检测装置中磁化线圈中的施加电流,实现对被测地上钢结构件中的磁化方向的任意调节,采集地上钢结构件在不同磁化方向下的漏磁检测信号;
步骤1包括:
步骤1.1,高精度漏磁检测装置水平放置在待检测的地上钢结构件上。
具体地,
本发明优选实施例中,高精度漏磁检测装置的励磁部分如图2所示,包括:第一铁芯、第二铁芯、第三铁芯、第四铁芯、第一磁化线圈、第二磁化线圈、第三磁化线圈、第四磁化线圈、磁轭、磁传感器阵列。
其中,磁传感器阵列水平布置在地上钢结构件上;第一铁芯、第二铁芯、第三铁芯、第四铁芯垂直于地上钢结构件,并按照十字形顶点位置呈逆时针方向依次分布在磁传感器阵列四周;缠绕在第一铁芯上的第一磁化线圈、缠绕在第二铁芯上的第二磁化线圈、缠绕在第三铁芯上的第三磁化线圈和缠绕在第四铁芯上的第四磁化线圈,缠绕匝数相同、缠绕方向也相同;铁芯上方放置磁轭,磁传感器阵列位于磁轭正下方,并且磁传感器阵列中心与磁轭中心的连线垂直于地上钢结构件。
磁轭是轴对称结构的,且至少有两条对称轴是相互垂直的。
值得注意的是,所属领域技术人员可以任意选择磁轭的形状和材料,对于磁轭的形状可以是但不限制于十字形、圆心、正方形、八边形等,对于磁轭的材料可以是但不限制于软铁板、A3钢板、软磁合金等。本发明的优选实施例磁轭采用十字形的软铁板是仅用于地上钢结构件腐蚀缺陷高精度漏磁成像检测非限制性的较佳选择。所属领域技术人员可以采用其它形状和材料的磁轭进行检测。
地上钢结构件、第一铁芯、第二铁芯、第三铁芯、第四铁芯和磁轭组成高精度漏磁检测装置的闭合磁回路。
步骤1.2,成对控制磁化线圈内施加的电流大小,实现磁轭中心的磁场方向沿逆时针扫描。
具体地,
在本发明优选实施例中,四组磁化线圈中每两组为一对,第二磁化线圈和第四磁化线圈为第一对磁化线圈,第一磁化线圈和第三磁化线圈为第二对磁化线圈。
第一对磁化线圈和第二对磁化线圈施加的电流大小分别为A1和A2,则A1和A2及磁轭中心处的磁场偏转角度θ满足如下关系:
Figure BDA0002986512640000081
式中,
A1为第一对磁化线圈的施加电流,A2为第二对磁化线圈的施加电流;
θ为磁轭中心的磁场方向的偏转角度,如图3所示,该偏转角度以仅对第二对磁化线圈施加电流时磁轭中心的磁场方向为基准,按照逆时针方向旋转一周计算,即从θ为0°沿着逆时针方向旋转至360°。
当A1为零且A2不为零时,磁轭中心的磁场方向与第二对磁化线圈的磁化方向平行时,θ为0°;当A1不为零且A2为零时,磁轭中心的磁场方向与第一对磁化线圈的磁化方向平行时,θ为90°。
当A1和A2大小相同、方向相反时,第一对磁化线圈的磁场方向与第二对磁化线圈的磁场方向相互垂直,并且K为恒定值,因此磁轭中心的磁场大小保持不变,并且磁场偏转角度θ为45°或135°,此时磁轭中心处的磁场沿固定方向不变。
因此,通过控制两对磁化线圈中的施加电流,实现磁轭中心的磁场偏转方向沿逆时针依次偏转0°、45°、90°和135°四个扫描角度。
实际检测过程中,根据地上钢结构件的检测需求优化磁场的扫描方向及路径;本发明优选实施例中,考虑漏磁场的对称特性,仅将磁轭中心的磁场方向沿逆时针扫描半周,即仅从θ为0°沿着逆时针方向旋转至180°;本发明优选实施例中,还考虑到检测数据存储容量的限制,仅扫描θ为0°、45°、90°和135°有限个数的磁场方向。
步骤1.3,在步骤1.2进行的同时,由磁传感器阵列采集地上钢结构件腐蚀缺陷处的漏磁检测信号。
步骤2,对步骤1采集的漏磁检测信号矩阵先直流滤波处理再归一化处理。
对于有腐蚀缺陷的地上钢结构件,对其漏磁检测信号矩阵的直流滤波处理具体包括:
(1)在相同环境条件下,采用相同的高精度漏磁检测装置采集无腐蚀缺陷的地上钢结构件的漏磁检测信号,得到无腐蚀缺陷的地上钢结构件的漏磁检测信号矩阵B0
本优选实施例中,无腐蚀缺陷的地上钢结构件的结构、形状、材质等尽量要与具有腐蚀缺陷的钢结构件一致。鉴于实际工程中,钢结构件的设计、制造以及安装均采用标准化流程,因此上述要求具备实现的可能性。
(2)在相同的磁场方向扫描条件下,从有腐蚀缺陷的地上钢结构件的漏磁检测信号Bi中减去无腐蚀缺陷的地上钢结构件的漏磁检测信号
Figure BDA0002986512640000082
即实现了直流滤波处理。
对于有腐蚀缺陷的地上钢结构件,将其直流滤波处理后的漏磁检测信号矩阵转换到[-1,1]范围内,实现归一化处理,满足如下关系式:
Figure BDA0002986512640000091
式中,
Norm(·)表示一种归一化处理函数,满足如下关系式:
Figure BDA0002986512640000092
式中,
max(·)表示取矩阵中所有数据的最大值,
min(·)表示取矩阵中所有数据的最小值;
Bi表示在不同磁场方向扫描时采集到的有腐蚀缺陷的地上钢结构件的漏磁检测信号,Bi=[B1,B2,...,Bn],i=1,2,...,n,定义n为不同磁场方向扫描时所记录的漏磁信号数量;
Figure BDA0002986512640000093
表示在不同磁场方向扫描时采集到的无腐蚀缺陷的地上钢结构件的漏磁检测信号,
Figure BDA0002986512640000094
Figure BDA0002986512640000095
表示有腐蚀缺陷的地上钢结构件的经直流滤波处理和归一化处理后的在不同磁场方向扫描时采集到的漏磁检测信号,
Figure BDA0002986512640000096
值得注意的是,所属领域技术人员可以任意选择不同数据处理软件中不同的归一化函数,本发明优选实施例所采用的归一化函数仅是非限制性的优选。
步骤3:根据漏磁检测信号的不同类型,对步骤2处理后的漏磁检测信号进行信号合成。
具体地,
漏磁检测信号包括漏磁检测信号的水平分量信号、漏磁检测信号的垂直分量信号或漏磁检测信号法向分量信号中的至少一种。
Bh表示漏磁检测信号的水平分量信号矩阵,Bhi表示在不同磁场方向扫描时采集到的漏磁检测信号的水平分量信号矩阵,Bh=[Bh1,Bh2,...,Bhn]。
Figure BDA0002986512640000097
表示经步骤2处理后的漏磁检测信号的水平分量信号矩阵,
Figure BDA0002986512640000098
表示在不同磁场方向扫描时经步骤2处理后的漏磁检测信号的水平分量信号矩阵,
Figure BDA0002986512640000099
Bv表示漏磁检测信号的垂直分量信号矩阵,Bvi表示在不同磁场方向扫描时采集到的漏磁检测信号的垂直分量信号矩阵,Bv=[Bv1,Bv2,...,Bvn]。
Figure BDA00029865126400000910
表示经步骤2处理后的漏磁检测信号的垂直分量信号矩阵,
Figure BDA00029865126400000911
表示在不同磁场方向扫描时经步骤2处理后的漏磁检测信号的垂直分量信号矩阵,
Figure BDA0002986512640000101
Bf表示漏磁检测信号的法向分量信号矩阵,Bfi表示在不同磁场方向扫描时采集到的漏磁检测信号的法向分量信号矩阵,Bf=[Bf1,Bf2,...,Bfn]。
Figure BDA0002986512640000102
表示经步骤2处理后的漏磁检测信号的法向分量信号矩阵,
Figure BDA0002986512640000103
表示在不同磁场方向扫描时经步骤2处理后的漏磁检测信号的法向分量信号矩阵,
Figure BDA0002986512640000105
根据漏磁检测信号的不同类型,对步骤2处理后的漏磁检测信号进行信号合成,具体如下:
(1)对于漏磁检测信号的水平分量信号,信号合成矩阵满足如下关系式:
Figure BDA0002986512640000106
式中,
Bh_combine表示步骤2处理后的漏磁检测信号的水平分量信号的合成信号矩阵,
Ghi表示
Figure BDA0002986512640000107
所对应的处理模块;
*表示卷积运算;
abs(·)表示对矩阵数据取绝对值的函数。
(2)对于漏磁检测信号的垂直分量信号,信号合成矩阵满足如下关系式:
Figure BDA0002986512640000108
式中,
Bv_combine表示步骤2处理后的漏磁检测信号的垂直分量信号的合成信号矩阵,
Gvi表示
Figure BDA0002986512640000109
所对应的处理模块。
具体地,
Ghi、Gvi与磁轭中心的磁场偏转角度θ有关,具体满足如下关系式:
当θ∈[0°,22.5°)∪[157.5°,202.5°)∪[337.5°,360°)时,
Figure BDA00029865126400001010
当θ∈[22.5°,67.5°)∪[202.5°,247.5°)时,
Figure BDA00029865126400001011
当θ∈[67.5°,112.5°)∪[247.5°,292.5°)时,
Figure BDA00029865126400001012
当θ∈[112.5°,157.5°)∪[292.5°,337.5°)时,
Figure BDA00029865126400001110
(3)对于漏磁检测信号的法向分量信号,信号合成矩阵满足如下关系式:
Figure BDA0002986512640000111
式中,
Bf_combine表示步骤2处理后的漏磁检测信号的法向分量信号的合成信号矩阵,
ki表示
Figure BDA0002986512640000112
的合成系数,满足0≤ki≤1。
步骤4:对步骤3的合成信号矩阵进行边缘提取处理,即可得到地上钢结构件腐蚀缺陷的轮廓,实现缺陷高精度成像。
值得注意的是,所属领域技术人员可以任意采用各类边缘检测算法实现边缘提取处理,这些算法包括但不限于Canny边缘检测算法、Sobel边缘检测算法、Laplacian边缘检测算法等。
实施例1。按照步骤1对有腐蚀缺陷的地上钢结构件的漏磁检测信号进行采集,磁轭中心的磁场偏转角度θ依次为0°、45°、90°和135°时,采集到漏磁检测信号的水平分量信号矩阵依次为Bh1、Bh2、Bh3和Bh4,经步骤2处理后的漏磁检测信号的水平分量信号矩阵依次为
Figure BDA0002986512640000113
Figure BDA0002986512640000114
这些矩阵的二维图像依次如图4中(a)、(b)、(c)和(d)所示;再对这些矩阵进行信号合成得到漏磁检测信号的水平分量信号的合成信号矩阵Bh_combine,该矩阵满足如下关系式:
Figure BDA0002986512640000115
式中,
Figure BDA0002986512640000116
Figure BDA0002986512640000117
Figure BDA0002986512640000118
Figure BDA0002986512640000119
的二维图像依次如图5中(a)、(b)、(c)和(d)所示。合成信号矩阵Bh_combine的二维图像如图6所示。
从图6中可以清晰地看出地上钢构件腐蚀缺陷的轮廓形状,对合成信号矩阵Bh_combine进行边缘提取处理,本发明优选实施例采用Canny边缘检测算法进行边缘提取处理,最终得到缺陷轮廓的高精度成像,该缺陷轮廓平均相对误差仅为2.18%,可见缺陷轮廓计算精度高。
实施例2。按照步骤1对有腐蚀缺陷的地上钢结构件的漏磁检测信号进行采集,磁轭中心的磁场偏转角度θ依次为0°、45°、90°和135°时,采集到漏磁检测信号的法向分量信号矩阵依次为Bf1、Bf2、Bf3和Bf4,经步骤2处理后的漏磁检测信号的法向分量信号矩阵依次为
Figure BDA0002986512640000121
Figure BDA0002986512640000122
这些矩阵的二维图像依次如图7中(a)、(b)、(c)和(d)所示;再对这些矩阵进行信号合成得到漏磁检测信号的法向分量信号的合成信号矩阵Bf_combine,该矩阵满足如下关系式:
Figure BDA0002986512640000123
式中,
本发明优选实施例中,ki均取1。
合成信号矩阵Bf_combine的二维图像如图8所示。
从图8中可以清晰地看出地上钢构件腐蚀缺陷的轮廓形状,对合成信号矩阵Bf_combine进行边缘提取处理,本发明优选实施例采用Canny边缘检测算法进行边缘提取处理,最终得到缺陷轮廓的高精度成像,该缺陷轮廓平均相对误差仅为2.12%,可见缺陷轮廓计算精度高。
地上钢结构件腐蚀缺陷高精度漏磁成像检测系统,包括漏磁检测信号采集模块、漏磁检测信号综合处理模块、漏磁检测信号反演成像模块。
漏磁检测信号采集模块,用于在不同磁化方向下采集漏磁检测信号;该模块内置高精度漏磁检测装置,通过成对控制高精度漏磁检测装置中磁化线圈中的施加电流,实现对被测地上钢结构件中的磁化方向的任意调节。
漏磁检测信号综合处理模块,用于对采集到的漏磁检测信号依次进行直流滤波处理、归一化处理、信号合成处理以及边缘提取处理。
漏磁检测信号反演成像模块,用于显示地上钢结构件腐蚀缺陷的三维图像。
优选地,漏磁检测信号综合处理模块包括直流滤波处理单元、归一化处理单元、水平分量和垂直分量信号合成处理单元、法向分量信号合成处理单元、边缘提取处理单元。
本发明的有益效果在于,与现有技术相比,通过调节两对磁化线圈中的电流大小实现对被测地上钢结构件中的磁化方向进行控制,结合对在不同磁化方向下得到的漏磁检测信号的综合处理,最终实现对腐蚀缺陷的高精度成像检测。该方法具有装置结构简单、操作方便、缺陷成像检测精度高等优点,有利于腐蚀缺陷的高精度成像。
本发明申请人结合说明书附图对本发明的实施示例做了详细的说明与描述,但是本领域技术人员应该理解,以上实施示例仅为本发明的优选实施方案,详尽的说明只是为了帮助读者更好地理解本发明精神,而并非对本发明保护范围的限制,相反,任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。

Claims (8)

1.面向成像的地上钢结构件腐蚀缺陷漏磁检测方法,其特征在于,
所述腐蚀缺陷漏磁检测方法的步骤如下:
步骤1,通过成对控制高精度漏磁检测装置中磁化线圈中的施加电流,实现对被测地上钢结构件中的磁化方向的任意调节,采集地上钢结构件在不同磁化方向下的漏磁检测信号;
步骤2,对步骤1采集的漏磁检测信号先直流滤波处理再归一化处理;
步骤3:根据漏磁检测信号的不同类型,对步骤2处理后的漏磁检测信号进行信号合成;
步骤4:对步骤3的合成信号进行边缘提取处理,得到地上钢结构件腐蚀缺陷的轮廓参数,实现缺陷高精度成像。
2.根据权利要求1所述的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法,其特征在于,
步骤1包括:
步骤1.1,高精度漏磁检测装置水平放置在待检测的地上钢结构件上;
所述高精度漏磁检测装置的励磁部分包括:第一铁芯、第二铁芯、第三铁芯、第四铁芯、第一磁化线圈、第二磁化线圈、第三磁化线圈、第四磁化线圈、磁轭、磁传感器阵列;
其中,磁传感器阵列水平布置在地上钢结构件上;第一铁芯、第二铁芯、第三铁芯、第四铁芯垂直于地上钢结构件,并按照十字形顶点位置呈逆时针方向依次分布在磁传感器阵列四周;缠绕在第一铁芯上的第一磁化线圈、缠绕在第二铁芯上的第二磁化线圈、缠绕在第三铁芯上的第三磁化线圈和缠绕在第四铁芯上的第四磁化线圈,缠绕匝数相同、缠绕方向也相同;铁芯上方放置磁轭,磁传感器阵列位于磁轭正下方,并且磁传感器阵列中心与磁轭中心的连线垂直于地上钢结构件;
地上钢结构件、第一铁芯、第二铁芯、第三铁芯、第四铁芯和磁轭组成高精度漏磁检测装置的闭合磁回路;
步骤1.2,成对控制磁化线圈内施加的电流大小,实现磁轭中心的磁场方向沿逆时针方向旋转扫描;
四组磁化线圈中每两组为一对,第二磁化线圈和第四磁化线圈为第一对磁化线圈,第一磁化线圈和第三磁化线圈为第二对磁化线圈;
第一对磁化线圈和第二对磁化线圈施加的电流大小分别为A1和A2,则A1和A2及磁轭中心处的磁场偏转角度θ满足如下关系:
Figure FDA0002986512630000021
式中,
A1为第一对磁化线圈的施加电流,A2为第二对磁化线圈的施加电流;
θ为磁轭中心的磁场方向的偏转角度,该偏转角度以仅对第二对磁化线圈施加电流时磁轭中心的磁场方向为基准,按照逆时针方向旋转一周计算,即从θ为0°沿着逆时针方向旋转至360°;
步骤1.3,在步骤1.2进行的同时,由磁传感器阵列采集地上钢结构件腐蚀缺陷处的漏磁检测信号。
3.根据权利要求2所述的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法,其特征在于,
在步骤1.1中,磁轭是轴对称结构的,且至少有两条对称轴是相互垂直的。
4.根据权利要求1所述的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法,其特征在于,
在步骤2中,
对于有腐蚀缺陷的地上钢结构件,对其漏磁检测信号矩阵的直流滤波处理具体包括:
(1)在相同环境条件下,采用相同的高精度漏磁检测装置采集无腐蚀缺陷的地上钢结构件的漏磁检测信号,得到无腐蚀缺陷的地上钢结构件的漏磁检测信号矩阵;
(2)在相同的磁场方向扫描条件下,从有腐蚀缺陷的地上钢结构件的漏磁检测信号中减去无腐蚀缺陷的地上钢结构件的漏磁检测信号,即实现了直流滤波处理;
对于有腐蚀缺陷的地上钢结构件,将其直流滤波处理后的漏磁检测信号矩阵转换到[-1,1]范围内,实现归一化处理,满足如下关系式:
Figure FDA0002986512630000022
式中,
Norm(·)表示一种归一化处理函数,
Bi表示在不同磁场方向扫描时采集到的有腐蚀缺陷的地上钢结构件的漏磁检测信号,Bi=[B1,B2,…,Bn],i=1,2,…,n,定义n为不同磁场方向扫描时所记录的漏磁信号数量;
Figure FDA0002986512630000023
表示在不同磁场方向扫描时采集到的无腐蚀缺陷的地上钢结构件的漏磁检测信号,
Figure FDA0002986512630000024
Figure FDA0002986512630000031
表示有腐蚀缺陷的地上钢结构件的经直流滤波处理和归一化处理后的在不同磁场方向扫描时采集到的漏磁检测信号,
Figure FDA0002986512630000032
5.根据权利要求1所述的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法,其特征在于,
在步骤3中,
漏磁检测信号包括漏磁检测信号的水平分量信号、漏磁检测信号的垂直分量信号或漏磁检测信号法向分量信号中的至少一种;根据漏磁检测信号的不同类型,对步骤2处理后的漏磁检测信号进行信号合成,具体如下:
(1)对于漏磁检测信号的水平分量信号,信号合成矩阵满足如下关系式:
Figure FDA0002986512630000033
式中,
Bh_combine表示步骤2处理后的漏磁检测信号的水平分量信号的合成信号矩阵,
Figure FDA0002986512630000034
表示在不同磁场方向扫描时经步骤2处理后的漏磁检测信号的水平分量信号矩阵,
Figure FDA0002986512630000035
Ghi表示
Figure FDA0002986512630000036
所对应的处理模块;
*表示卷积运算;
abs(·)表示对矩阵数据取绝对值的函数;
(2)对于漏磁检测信号的垂直分量信号,信号合成矩阵满足如下关系式:
Figure FDA0002986512630000037
式中,
Bv_combine表示步骤2处理后的漏磁检测信号的垂直分量信号的合成信号矩阵,
Figure FDA0002986512630000038
表示在不同磁场方向扫描时经步骤2处理后的漏磁检测信号的垂直分量信号矩阵,
Figure FDA0002986512630000039
Gvi表示
Figure FDA00029865126300000310
所对应的处理模块;
(3)对于漏磁检测信号的法向分量信号,信号合成矩阵满足如下关系式:
Figure FDA00029865126300000311
式中,
Bf_combine表示步骤2处理后的漏磁检测信号的法向分量信号的合成信号矩阵,
Figure FDA0002986512630000041
表示在不同磁场方向扫描时经步骤2处理后的漏磁检测信号的法向分量信号矩阵,
Figure FDA0002986512630000042
ki表示
Figure FDA0002986512630000043
的合成系数,满足0≤ki≤1。
6.根据权利要求5所述的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法,其特征在于,
所述处理模块Ghi、Gvi与磁轭中心的磁场偏转角度θ有关,满足如下关系式:
当θ∈[0°,22.5°)∪[157.5°,202.5°)∪[337.5°,360°)时,
Figure FDA0002986512630000044
当θ∈[22.5°,67.5°)∪[202.5°,247.5°)时,
Figure FDA0002986512630000045
当θ∈[67.5°,112.5°)∪[247.5°,292.5°)时,
Figure FDA0002986512630000046
当θ∈[112.5°,157.5°)∪[292.5°,337.5°)时,
Figure FDA0002986512630000047
7.一种基于权利要求1至6所述的面向成像的地上钢结构件腐蚀缺陷漏磁检测方法的面向成像的地上钢结构件腐蚀缺陷漏磁检测系统,包括漏磁检测信号采集模块、漏磁检测信号综合处理模块、漏磁检测信号反演成像模块,其特征在于,
漏磁检测信号采集模块,用于在不同磁化方向下采集漏磁检测信号;该模块内置高精度漏磁检测装置,通过成对控制高精度漏磁检测装置中磁化线圈中的施加电流,实现对被测地上钢结构件中的磁化方向的任意调节;
漏磁检测信号综合处理模块,用于对采集到的漏磁检测信号依次进行直流滤波处理、归一化处理、信号合成处理以及边缘提取处理;
漏磁检测信号反演成像模块,用于显示地上钢结构件腐蚀缺陷的三维图像。
8.根据权利要求7所述的地上钢结构件腐蚀缺陷高精度漏磁成像检测系统,其特征在于:所述漏磁检测信号综合处理模块包括直流滤波处理单元、归一化处理单元、水平分量和垂直分量信号合成处理单元、法向分量信号合成处理单元、边缘提取处理单元。
CN202110301694.2A 2021-03-22 2021-03-22 面向成像的地上钢结构件腐蚀缺陷漏磁检测方法及系统 Active CN113030243B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110301694.2A CN113030243B (zh) 2021-03-22 2021-03-22 面向成像的地上钢结构件腐蚀缺陷漏磁检测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110301694.2A CN113030243B (zh) 2021-03-22 2021-03-22 面向成像的地上钢结构件腐蚀缺陷漏磁检测方法及系统

Publications (2)

Publication Number Publication Date
CN113030243A true CN113030243A (zh) 2021-06-25
CN113030243B CN113030243B (zh) 2023-08-01

Family

ID=76472181

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110301694.2A Active CN113030243B (zh) 2021-03-22 2021-03-22 面向成像的地上钢结构件腐蚀缺陷漏磁检测方法及系统

Country Status (1)

Country Link
CN (1) CN113030243B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114674914A (zh) * 2022-04-07 2022-06-28 中国矿业大学 一种消除管道周向励磁漏磁检测时背影磁场影响的方法
KR102607277B1 (ko) * 2023-06-07 2023-11-29 한전케이피에스 주식회사 결함 감지 시스템 및 그 방법

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59217159A (ja) * 1983-05-26 1984-12-07 Mitsubishi Electric Corp 欠陥検出装置
JPH04313057A (ja) * 1991-01-19 1992-11-05 Japan Aircraft Mfg Co Ltd 探傷装置
CN202083672U (zh) * 2011-05-19 2011-12-21 中国人民解放军国防科学技术大学 基于场量检测的漏磁探头
CN105021694A (zh) * 2015-07-02 2015-11-04 清华大学 不完整信号下的漏磁检测缺陷量化与显示方法
CN106290552A (zh) * 2016-08-01 2017-01-04 四川大学 一种基于旋转磁化场的钢板漏磁检测装置
JP2017150904A (ja) * 2016-02-23 2017-08-31 東京ガスエンジニアリングソリューションズ株式会社 探傷装置および探傷方法
CN107607612A (zh) * 2017-08-11 2018-01-19 清华大学 基于漏磁信号的垂直分量的缺陷轮廓识别方法及装置
CN107741454A (zh) * 2017-09-04 2018-02-27 清华大学 基于法向分量的伪三维漏磁信号缺陷轮廓识别方法
CN111929356A (zh) * 2020-07-07 2020-11-13 清华大学 钢材缺陷磁成像装置及方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59217159A (ja) * 1983-05-26 1984-12-07 Mitsubishi Electric Corp 欠陥検出装置
JPH04313057A (ja) * 1991-01-19 1992-11-05 Japan Aircraft Mfg Co Ltd 探傷装置
CN202083672U (zh) * 2011-05-19 2011-12-21 中国人民解放军国防科学技术大学 基于场量检测的漏磁探头
CN105021694A (zh) * 2015-07-02 2015-11-04 清华大学 不完整信号下的漏磁检测缺陷量化与显示方法
JP2017150904A (ja) * 2016-02-23 2017-08-31 東京ガスエンジニアリングソリューションズ株式会社 探傷装置および探傷方法
CN106290552A (zh) * 2016-08-01 2017-01-04 四川大学 一种基于旋转磁化场的钢板漏磁检测装置
CN107607612A (zh) * 2017-08-11 2018-01-19 清华大学 基于漏磁信号的垂直分量的缺陷轮廓识别方法及装置
CN107741454A (zh) * 2017-09-04 2018-02-27 清华大学 基于法向分量的伪三维漏磁信号缺陷轮廓识别方法
CN111929356A (zh) * 2020-07-07 2020-11-13 清华大学 钢材缺陷磁成像装置及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HONGYU SUN等: ""3D focusing acoustic lens optimization method using multi-factor and multi-level orthogonal test designing theory "", 《APPLIED ACOUSTICS》, pages 1 - 9 *
ZHANG JIAYING等: ""A comparative study between magnetice field distortion and magnetic flux leakage techniques for surface defect shape reconstruction in steel plates"", 《SENSORS AND ACTUATORS A:PHYSICAL》, pages 10 - 20 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114674914A (zh) * 2022-04-07 2022-06-28 中国矿业大学 一种消除管道周向励磁漏磁检测时背影磁场影响的方法
KR102607277B1 (ko) * 2023-06-07 2023-11-29 한전케이피에스 주식회사 결함 감지 시스템 및 그 방법

Also Published As

Publication number Publication date
CN113030243B (zh) 2023-08-01

Similar Documents

Publication Publication Date Title
CN113030243B (zh) 面向成像的地上钢结构件腐蚀缺陷漏磁检测方法及系统
CN106657789A (zh) 一种螺纹全景图像合成方法
CN104200442A (zh) 基于改进的canny边缘检测的非局部均值MRI图像去噪方法
CN110992363B (zh) 一种基于生成对抗网络的三维漏磁检测缺陷重构方法
CN106856003A (zh) 轴类工件侧表面缺陷检测图像的展开校正方法
TWI521448B (zh) Vehicle identification system and method
DE102014212499A1 (de) Verfahren und Vorrichtung zur Streuflussprüfung
CN111383194B (zh) 一种基于极坐标的相机畸变图像校正方法
CN109682824B (zh) 基于图像融合的钢丝绳无损检测方法及其装置
CN104794692A (zh) 一种图像去锯齿的系统
CN111462198A (zh) 一种尺度、旋转和辐射不变性的多模态影像配准方法
CN113702384A (zh) 一种回转构件表面缺陷检测装置、检测方法及标定方法
CN115963397A (zh) 一种电机定子内轮廓表面缺陷快速在线检测方法及装置
CN115601407A (zh) 一种红外与可见光图像配准方法
US10429471B2 (en) MRI with variable density sampling
CN112396651B (zh) 一种通过二角度图像实现设备定位的方法
Angani et al. Magnetic leakage testing using linearly integrated hall and GMR sensor arrays to inspect inclusions in cold-rolled strips
CN108537831B (zh) 对增材制造工件进行ct成像的方法及装置
CN111028215A (zh) 一种基于机器视觉的检测钢卷端面缺陷方法
CN106767566A (zh) 一种工件质量在线监测评估方法及监测系统
CN109345501A (zh) 基于频域图像增强的螺旋缠绕式物料缺陷检测方法
CN109443220A (zh) 基于测行程传感器和机器视觉的绕线机线圈长度测量方法
CN115598177A (zh) 基于微分法的车载电缆内部热缺陷轮廓检测方法
CN113240752B (zh) 一种内参和外参协同标定方法和装置
Vervaeke Inline magnet inspection using fast high resolution MagCam magnetic field mapping and analysis

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