CN111207854A - 一种基于分布式光纤温度传感器的数据处理算法 - Google Patents

一种基于分布式光纤温度传感器的数据处理算法 Download PDF

Info

Publication number
CN111207854A
CN111207854A CN202010032159.7A CN202010032159A CN111207854A CN 111207854 A CN111207854 A CN 111207854A CN 202010032159 A CN202010032159 A CN 202010032159A CN 111207854 A CN111207854 A CN 111207854A
Authority
CN
China
Prior art keywords
optical fiber
signal
stokes
stokes signal
temperature sensor
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
CN202010032159.7A
Other languages
English (en)
Other versions
CN111207854B (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.)
Guoxing Huijin Shenzhen Technology Co Ltd
Original Assignee
Guoxing Huijin Shenzhen Technology 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 Guoxing Huijin Shenzhen Technology Co Ltd filed Critical Guoxing Huijin Shenzhen Technology Co Ltd
Priority to CN202010032159.7A priority Critical patent/CN111207854B/zh
Publication of CN111207854A publication Critical patent/CN111207854A/zh
Application granted granted Critical
Publication of CN111207854B publication Critical patent/CN111207854B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/32Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K15/00Testing or calibrating of thermometers
    • G01K15/005Calibration

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)

Abstract

本发明涉及一种基于分布式光纤温度传感器的数据处理算法,包括设置系统参数并创建测量任务;打开激光器发送脉冲,数据采集模块采集光纤的斯托克斯信号和反斯托克斯信号并分别求取其平均值;关闭激光器,数据采集模块采集背景噪声信号并求取其平均值;分别将斯托克斯信号和反斯托克斯信号的平均值减去背景噪声信号的平均值,并进行滑窗平滑处理,然后利用插值算法将经滑窗平滑处理的斯托克斯信号和反斯托克斯信号拉伸至固定的采样空间间隔内,计算光纤同一点位置上的斯托克斯信号和反斯托克斯信号的比值序列以得到光纤的温度序列。本发明提高了光纤温度分辨率的一致性,提高了长距离的测量精度,进而提高了整个分布式光纤温度传感器的应用范围。

Description

一种基于分布式光纤温度传感器的数据处理算法
技术领域
本发明涉及光纤传感技术领域,特别涉及一种基于分布式光纤温度传感器的数据处理算法。
背景技术
分布式光纤温度传感器主要依据光纤的光时域反射和光纤的背向拉曼散射温度效应,测量沿光纤每一点的温度。现在市场上的分布式光纤温度传感器采用的都是单端口光纤配置模式,无法保证其温度测量精度,同时温度分辨率随着光纤长度的增加而变差,测量不准确。因此,现有技术有待发展。
发明内容
本发明的目的在于克服现有技术中存在的缺点,提供一种提高光纤各区域温度分辨率稳定性的基于分布式光纤温度传感器的数据处理算法。
为了实现上述目的,本发明提供如下技术方案:
本发明提供了一种基于分布式光纤温度传感器的数据处理算法,包括以下步骤:
S100,设置系统参数并创建测量任务;
S200,反射信号采集:打开激光器发送脉冲,数据采集模块采集光纤的斯托克斯信号和反斯托克斯信号,并分别求取斯托克斯信号和反斯托克斯信号的平均值;
S300,噪声信号采集:关闭激光器,数据采集模块采集背景噪声信号,并求取背景噪声信号的平均值;
S400,数据处理:分别将斯托克斯信号和反斯托克斯信号的平均值减去背景噪声信号的平均值,得到去噪后的斯托克斯信号和反斯托克斯信号并对其进行滑窗平滑处理,利用插值算法将经过滑窗平滑处理的斯托克斯信号和反斯托克斯信号拉伸至固定的采样空间间隔内,计算光纤同一点位置上的斯托克斯信号和反斯托克斯信号的比值序列,以得到光纤的温度序列。
进一步地,噪声信号采集还包括对数据采集模块采集到的背景噪声信号进行滤波处理,以去除背景噪声信号中的高频部分。
进一步地,所述系统参数包括背景噪声采样次数、后向散射信号采样次数、采样空间间隔以及单、双端口模式等相关参数。
进一步地,对光纤的温度序列计算方差。
进一步地,当采用双端口模式时,将光纤两端的数据对齐处理,得到光纤同一点位置上的斯托克斯信号和反斯托克斯信号的比值序列,以得到光纤的温度序列。
进一步地,通过两个端口的斯托克斯信号和反斯托克斯信号的比值序列计算损失校正参数,得到损失校正参数序列。
进一步地,分别计算两个端口的斯托克斯信号和反斯托克斯信号在每一点温度的权重,并根据两者的权重计算得到温度对数曲线η序列,最后根据温度对数曲线η序列得到优化后的温度序列。
进一步地,还包括S500,数据打包:将原始数据打包,便于后续软件的解析处理。
本发明技术方案具有的有益效果:
本发明的基于分布式光纤温度传感器的数据处理算法,通过反射信号的采集、噪声信号的采集以及数据处理得到光纤的温度序列,提高了光纤每个测量点上温度分辨率的一致性,提高了长距离的测量精度,进而提高了整个分布式光纤温度传感器的应用范围。
附图说明
图1是本发明的流程图;
图2是本发明反射信号采集的流程图;
图3是本发明噪声信号采集的流程图;
图4是本发明数据处理的流程图;
图5是本发明总的流程图;
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明,本发明实施例中所有方向性指示(诸如上、下、左、右、前、后……)仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等,如果该特定姿态发生改变时,则该方向性指示也相应地随之改变。
在本发明中,除非另有明确的规定和限定,术语“连接”、“固定”等应做广义理解,例如,“连接”可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
另外,在本发明中如涉及“第一”、“第二”等的描述仅用于描述目的,而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。
请参考图1至图5,本发明提供了一种基于分布式光纤温度传感器的数据处理算法,包括以下步骤:
S100,设置系统参数并创建测量任务,PC端配置软件通过网络登录到分布式光纤温度传感器后,发送该设备对应的配置系统参数,并根据接收到的配置系统参数创建测量任务,每个测量任务中包含其对应测量使用得到的光开关端口、光纤长度、背景噪声采样次数、真实信号采样次数以及单双端模式等信息;
S200,反射信号采集:测量开始后,光开关切到对应端口,打开激光器发送脉冲,数据采集模块同时开始采集光纤的斯托克斯信号和反斯托克斯信号,并分别求取斯托克斯信号和反斯托克斯信号的平均值,即进行累加求平均值;
S300,噪声信号采集:在完成上述真实信号的数据采集之后,关闭激光器,光开关切换到断开位置,按照测量任务配置的次数测量背景噪声信号,数据采集模块采集背景噪声信号,并求取背景噪声信号的平均值,即进行累加求平均值;其中,在系统稳定工作后,不过激光器是否开启,背景噪声一直存在。
S400,数据处理:在反射信号采集中采集到的斯托克斯信号和反斯托克斯信号是包含背景噪声在内的,因此分别将斯托克斯信号和反斯托克斯信号的平均值减去背景噪声信号的平均值,得到去噪后的斯托克斯信号和反斯托克斯信号并对其进行滑窗平滑处理,利用插值算法将经过滑窗平滑处理的斯托克斯信号和反斯托克斯信号拉伸至固定的采样空间间隔内,计算光纤同一点位置上的斯托克斯信号和反斯托克斯信号的比值序列,以得到光纤的温度序列。
在本实施方式中,对去噪后的斯托克斯信号和反斯托克斯信号进行滑窗平滑处理是根据该数据左右相邻N个点的不同权重重新计算该位置的值,左右N个点再加上自身,即为奇数个点参与整个算法,也就是说滑窗的大小必须为奇数。滑窗平滑处理的目的是尽量消除或降低信号过冲带来的影响,而过冲效果是指在温度突变位置或者光纤接头位置散射信号突然发生变化导致测量结果出现突起或下陷的尖峰,这种物理现象对测量准确性有干扰,所以需要通过该步骤进行消除。
在本实施方式中,由于斯托克斯信号和反斯托克斯信号传播的速率不一样,而数据采集模块采集数据是同时进行的,因此采集到1000米远位置散射回来的斯托克斯信号与同时采集到的反斯托克斯信号并不是来自同样位置上的。而计算温度需要同样位置上的斯托克斯信号和反斯托克斯信号,因此需要将两个信号每次采集的结果都对齐到相同位置,插值算法就是将斯托克斯信号和反斯托克斯信号的值都缩放到相同距离的位置,以保证在计算温度时不会出现位置偏移。
在本实施方式中,斯托克斯信号和反斯托克斯信号的折射率存在细微差别,两者在光纤中的传输速度不一致,这意味着斯托克斯信号和反斯托克斯信号两次采样之间的距离是不相等的。而分布式光纤温度传感器最终输出某个位置上的温度值是通过该位置上的斯托克斯信号和反斯托克斯信号的比值关系得到的,因此斯托克斯信号和反斯托克斯信号在空间上的不一致,会导致两者的位置偏移随着距离的增加而越来越大,测量得到的温度数据也会越来越不准确。通过插值算法,将斯托克斯信号和反斯托克斯信号都拉伸到固定的采样空间间隔内,以保证在计算温度时不会出现位置偏移,从而提高光纤每个测量点上温度分辨率的一致性,提高了长距离的测量精度,进而提高了整个分布式光纤温度传感器的应用范围。
本发明的基于分布式光纤温度传感器的数据处理算法,通过反射信号的采集、噪声信号的采集以及数据处理得到光纤的温度序列,提高了光纤每个测量点上温度分辨率的一致性,提高了长距离的测量精度,进而提高了整个分布式光纤温度传感器的应用范围。
优选地,噪声信号采集还包括对数据采集模块采集到的背景噪声信号进行滤波处理,以去除背景噪声信号中的高频部分,保证背景噪声信号的平稳,有利于提高测量结果的精度。
在本实施方式中,所述系统参数包括背景噪声采样次数、后向散射信号采样次数、采样空间间隔以及单、双端口模式等相关参数,有助于之后的步骤顺利进行,具体的系统参数如表1所示:
Figure BDA0002364713690000061
表1
进一步地,对光纤的温度序列计算方差,可以知道斯托克斯信号和反斯托克斯信号等相关信号的稳定性,对于检测光纤上的信号是否异常有帮助,也参与了信号处理中各个数据的补偿处理。
作为一种实施例,当采用双端口模式时,将光纤两端的数据对齐处理,得到光纤同一点位置上的斯托克斯信号和反斯托克斯信号的比值序列,以得到光纤的温度序列,通过光纤两端的数据对齐,得到的温度分辨率一致性稳定。
进一步地,通过两个端口的斯托克斯信号和反斯托克斯信号的比值序列计算损失校正参数,得到损失校正参数序列。光在传播的过程中会逐渐衰减,散射回来的信号也会随着距离的增加而越来越弱,通过光纤自身的衰减特性参数计算的损失校正参数,可以根据光纤长度对采集到的信号进行补偿。
优选地,分别计算两个端口的斯托克斯信号和反斯托克斯信号在每一点温度的权重,并根据两者的权重计算得到温度对数曲线η序列,最后根据温度对数曲线η序列得到优化后的温度序列。在双端口的温度计算中,信噪比好的那端数据所占权重更大,而温度对数曲线η序列就是根据两个端口各自的曲线计算出一条更合理的温度对数曲线,对温度对数曲线η序列取反对数得到优化后的温度序列。
具体地,所述的基于分布式光纤温度传感器的数据处理算法还包括S500,数据打包:由于原始数据没有数据结构,因此需要将原始数据打包,便于后续软件的解析处理。数据打包后还会经过Zone数据计算以及Alarm数据计算,其中Zone是光纤上用户指定的区域,是根据客户的应用场景指定的,对用户感兴趣的区域温度进行算法处理也是根据用户自身需求来定义的;Alarm数据计算是在Zone计算结果的基础上来进行的检测,判断相关结果是否满足特定的范围,如果超出预设范围,需要根据Alarm定义执行其他操作,比如软件弹出告警提示等。
本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无法对所有的实施方式予以穷举。凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。

Claims (8)

1.一种基于分布式光纤温度传感器的数据处理算法,其特征在于,包括以下步骤:
S100,设置系统参数并创建测量任务;
S200,反射信号采集:打开激光器发送脉冲,数据采集模块采集光纤的斯托克斯信号和反斯托克斯信号,并分别求取斯托克斯信号和反斯托克斯信号的平均值;
S300,噪声信号采集:关闭激光器,数据采集模块采集背景噪声信号,并求取背景噪声信号的平均值;
S400,数据处理:分别将斯托克斯信号和反斯托克斯信号的平均值减去背景噪声信号的平均值,得到去噪后的斯托克斯信号和反斯托克斯信号并对其进行滑窗平滑处理,利用插值算法将经过滑窗平滑处理的斯托克斯信号和反斯托克斯信号拉伸至固定的采样空间间隔内,计算光纤同一点位置上的斯托克斯信号和反斯托克斯信号的比值序列,以得到光纤的温度序列。
2.根据权利要求1所述的基于分布式光纤温度传感器的数据处理算法,其特征在于,噪声信号采集还包括对数据采集模块采集到的背景噪声信号进行滤波处理,以去除背景噪声信号中的高频部分。
3.根据权利要求1所述的基于分布式光纤温度传感器的数据处理算法,其特征在于,所述系统参数包括背景噪声采样次数、后向散射信号采样次数、采样空间间隔以及单、双端口模式等相关参数。
4.根据权利要求1所述的基于分布式光纤温度传感器的数据处理算法,其特征在于,对光纤的温度序列计算方差。
5.根据权利要求3所述的基于分布式光纤温度传感器的数据处理算法,其特征在于,当采用双端口模式时,将光纤两端的数据对齐处理,得到光纤同一点位置上的斯托克斯信号和反斯托克斯信号的比值序列,以得到光纤的温度序列。
6.根据权利要求5所述的基于分布式光纤温度传感器的数据处理算法,其特征在于,通过两个端口的斯托克斯信号和反斯托克斯信号的比值序列计算损失校正参数,得到损失校正参数序列。
7.根据权利要求5所述的基于分布式光纤温度传感器的数据处理算法,其特征在于,分别计算两个端口的斯托克斯信号和反斯托克斯信号在每一点温度的权重,并根据两者的权重计算得到温度对数曲线η序列,最后根据温度对数曲线η序列得到优化后的温度序列。
8.根据权利要求1所述的基于分布式光纤温度传感器的数据处理算法,其特征在于,还包括S500,数据打包:将原始数据打包,便于后续软件的解析处理。
CN202010032159.7A 2020-01-13 2020-01-13 一种基于分布式光纤温度传感器的数据处理算法 Active CN111207854B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010032159.7A CN111207854B (zh) 2020-01-13 2020-01-13 一种基于分布式光纤温度传感器的数据处理算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010032159.7A CN111207854B (zh) 2020-01-13 2020-01-13 一种基于分布式光纤温度传感器的数据处理算法

Publications (2)

Publication Number Publication Date
CN111207854A true CN111207854A (zh) 2020-05-29
CN111207854B CN111207854B (zh) 2021-10-29

Family

ID=70784259

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010032159.7A Active CN111207854B (zh) 2020-01-13 2020-01-13 一种基于分布式光纤温度传感器的数据处理算法

Country Status (1)

Country Link
CN (1) CN111207854B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113008407A (zh) * 2021-02-26 2021-06-22 中煤科工集团重庆研究院有限公司 多参数分布式光纤信号解调仪
CN113155309A (zh) * 2021-03-09 2021-07-23 电子科技大学 一种分布式光纤拉曼测温系统中光纤色散的补偿方法
CN116105891A (zh) * 2023-02-10 2023-05-12 之江实验室 高温传感装置及高温传感装置的标定方法
CN116678506A (zh) * 2023-08-02 2023-09-01 国检测试控股集团南京国材检测有限公司 一种无线传输热损检测装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5618108A (en) * 1992-06-16 1997-04-08 Kabushiki Kaisha Toshiba Temperature distribution measuring apparatus using an optical fiber
CN101256098A (zh) * 2008-03-14 2008-09-03 秦一涛 一种分布式光纤温度传感器系统的温度测量方法
CN103033283A (zh) * 2011-09-30 2013-04-10 上海华魏光纤传感技术有限公司 一种分布式光纤测温系统降噪方法
CN103115693A (zh) * 2013-01-17 2013-05-22 长飞光纤光缆有限公司 一种分布式光纤拉曼测温系统
RU2583060C1 (ru) * 2015-02-09 2016-05-10 Общество с ограниченной ответственностью "ПетроФайбер" Способ измерения температурного распределения в объекте и устройство для его осуществления
CN105953941A (zh) * 2016-04-29 2016-09-21 深圳艾瑞斯通技术有限公司 一种基于拉曼散射的分布式光纤测温方法及装置
CN106323500A (zh) * 2015-07-08 2017-01-11 中国电力科学研究院 一种实现温度自校准的光纤拉曼测温系统及其校准方法
CN110375880A (zh) * 2019-08-15 2019-10-25 泰华智慧产业集团股份有限公司 分布式光纤测温系统及温度动态测量方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5618108A (en) * 1992-06-16 1997-04-08 Kabushiki Kaisha Toshiba Temperature distribution measuring apparatus using an optical fiber
CN101256098A (zh) * 2008-03-14 2008-09-03 秦一涛 一种分布式光纤温度传感器系统的温度测量方法
CN103033283A (zh) * 2011-09-30 2013-04-10 上海华魏光纤传感技术有限公司 一种分布式光纤测温系统降噪方法
CN103115693A (zh) * 2013-01-17 2013-05-22 长飞光纤光缆有限公司 一种分布式光纤拉曼测温系统
RU2583060C1 (ru) * 2015-02-09 2016-05-10 Общество с ограниченной ответственностью "ПетроФайбер" Способ измерения температурного распределения в объекте и устройство для его осуществления
CN106323500A (zh) * 2015-07-08 2017-01-11 中国电力科学研究院 一种实现温度自校准的光纤拉曼测温系统及其校准方法
CN105953941A (zh) * 2016-04-29 2016-09-21 深圳艾瑞斯通技术有限公司 一种基于拉曼散射的分布式光纤测温方法及装置
CN110375880A (zh) * 2019-08-15 2019-10-25 泰华智慧产业集团股份有限公司 分布式光纤测温系统及温度动态测量方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113008407A (zh) * 2021-02-26 2021-06-22 中煤科工集团重庆研究院有限公司 多参数分布式光纤信号解调仪
CN113008407B (zh) * 2021-02-26 2024-04-09 中煤科工集团重庆研究院有限公司 多参数分布式光纤信号解调仪
CN113155309A (zh) * 2021-03-09 2021-07-23 电子科技大学 一种分布式光纤拉曼测温系统中光纤色散的补偿方法
CN116105891A (zh) * 2023-02-10 2023-05-12 之江实验室 高温传感装置及高温传感装置的标定方法
CN116105891B (zh) * 2023-02-10 2024-01-09 之江实验室 高温传感装置及高温传感装置的标定方法
CN116678506A (zh) * 2023-08-02 2023-09-01 国检测试控股集团南京国材检测有限公司 一种无线传输热损检测装置
CN116678506B (zh) * 2023-08-02 2023-10-10 国检测试控股集团南京国材检测有限公司 一种无线传输热损检测装置

Also Published As

Publication number Publication date
CN111207854B (zh) 2021-10-29

Similar Documents

Publication Publication Date Title
CN111207854B (zh) 一种基于分布式光纤温度传感器的数据处理算法
EP2290320A1 (en) Film thickness measurement apparatus
CA2469817A1 (en) System and method for measuring optical distance
US20060285105A1 (en) Polarized lightwave reflectometry method (potdr)
JPH0499943A (ja) 伝送媒体の試験方法
WO2004081524B1 (en) Measurement site dependent data preprocessing method for robust calibration and prediction
US5442434A (en) Method for finding and measuring optical features using an optical time domain reflectometer
EP0926479A1 (en) Optical loss measurement
CN109752720B (zh) 基于分布式振动传感器振动信息定位方法及系统
CN114696893B (zh) 基于otdr的事件定位方法、系统、装置及存储介质
CN112285742A (zh) 一种在频域内估算相干测风激光雷达载噪比的方法
JPH06105217B2 (ja) 分光測定法
CN112033572B (zh) 一种光时域反射计中补偿弯曲损耗及斯托克斯光与反斯托克斯光衰减的方法
CN111289142B (zh) 一种分布式光纤拉曼测温系统信号自适应干扰抵消方法
CN112229501A (zh) 一种基于ofdr的振动信号自动检测装置和方法
CN111207857B (zh) 一种利用分布式光纤温度传感器测量光纤长度的方法
CN112697303B (zh) 面向智能电网的分布式光纤传感系统和检测方法
CN106353009B (zh) Dts系统衰减系数不同时标定方法
CN210327579U (zh) 一种光缆振动的定位装置
CN108303196B (zh) 全同光纤光栅数字化解调及区域温度监测方法及其系统
CN112268608A (zh) 一种强干扰dvs的高精度解调算法
CN105030193A (zh) 一种背景多普勒消除系统及方法
CN205449326U (zh) 一种基于光纤fp腔实现激光拍频测温的装置
CN115684076B (zh) 多通道红外气体传感器数据处理方法
CN112050967B (zh) 一种光纤温度分布测试仪的光纤温度自动校准及补偿方法

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