CN111427265A - 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质 - Google Patents

基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质 Download PDF

Info

Publication number
CN111427265A
CN111427265A CN202010194359.2A CN202010194359A CN111427265A CN 111427265 A CN111427265 A CN 111427265A CN 202010194359 A CN202010194359 A CN 202010194359A CN 111427265 A CN111427265 A CN 111427265A
Authority
CN
China
Prior art keywords
dictionary
heavy metal
wastewater treatment
metal wastewater
sample
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
CN202010194359.2A
Other languages
English (en)
Other versions
CN111427265B (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.)
Central South University
Original Assignee
Central South 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 Central South University filed Critical Central South University
Priority to CN202010194359.2A priority Critical patent/CN111427265B/zh
Publication of CN111427265A publication Critical patent/CN111427265A/zh
Priority to PCT/CN2021/077910 priority patent/WO2021185044A1/zh
Priority to US17/908,922 priority patent/US20230089156A1/en
Application granted granted Critical
Publication of CN111427265B publication Critical patent/CN111427265B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F1/00Treatment of water, waste water, or sewage
    • C02F1/008Control or steering systems not provided for elsewhere in subclass C02F
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F1/00Treatment of water, waste water, or sewage
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F1/00Treatment of water, waste water, or sewage
    • C02F1/46Treatment of water, waste water, or sewage by electrochemical methods
    • C02F1/461Treatment of water, waste water, or sewage by electrochemical methods by electrolysis
    • C02F1/463Treatment of water, waste water, or sewage by electrochemical methods by electrolysis by electrocoagulation
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2101/00Nature of the contaminant
    • C02F2101/10Inorganic compounds
    • C02F2101/20Heavy metals or heavy metal compounds
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/005Processes using a programmable logic controller [PLC]
    • C02F2209/006Processes using a programmable logic controller [PLC] comprising a software program or a logic diagram
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/05Conductivity or salinity
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/06Controlling or monitoring parameters in water treatment pH
    • CCHEMISTRY; METALLURGY
    • C02TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02FTREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00Controlling or monitoring parameters in water treatment
    • C02F2209/40Liquid flow rate

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Hydrology & Water Resources (AREA)
  • Environmental & Geological Engineering (AREA)
  • Water Supply & Treatment (AREA)
  • Organic Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Electrochemistry (AREA)
  • General Chemical & Material Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Removal Of Specific Substances (AREA)

Abstract

本发明公开了一种基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质基于迁移学习的重金属废水处理过程异常工况智能化监测,通过对不同来源的重金属废水处理过程数据融合,能够自动的实现不同来源的重金属废水处理过程异常工况智能识别;具体为利用来源固定的重金属废水处理过程的正常样本YSD、少量来源未知的重金属废水处理过程的正常样本YTD;首先通过对YSD进行学习得到其数据表示字典DSD,然后考虑到YSD和YTD分布不同,采用迁移学习的方法,将YTD的特征融入到字典学习过程,得到泛化能力更强的字典DTD。该方法无需过程先验知识,能自适应的适应废水处理系统中的不确定性因素,能够更加准确的检测过程中相关指标的变化,实现及时地检测与预警。

Description

基于迁移学习的重金属废水处理过程异常工况智能化监测方 法、装置及存储介质
技术领域
本发明属于重金属废水处理领域,特别涉及一种基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质。
背景技术
作为地球上不可替代的自然资源与环境资源,水的储量与质量与人类的生存发展息息相关。随着近些年工业的高速发展,工业废水尤其是重金属废水对环境污染非常严重和对人类危害非常深远。同时,重金属在环境中稳定性高、难降解、迁移范围广,正逐渐成为全球性环境问题。
絮凝—电化学技术是实现重金属废水深度净化的有效途径之一。该技术用电化学技术对经传统絮凝沉淀过程处理后的废水进行深度处理,可以减少药剂的消耗,与单一的电絮凝处理废水相比,又可以节约电能。一方面,现有絮凝—电化学技术对废水的处理效果受多种过程因素(如pH值、电流密度、电导率等)的影响。另一方面,重金属废水因其来源等存在不确定性,包括气候变化、人为破坏、偶然污染、内部许多生物化学反应过程的机理不明确等,难以用明确的数学关系对重金属废水的特征进行描述。除此之外,不同来源的重金属废水的离子浓度分布差异较大,其处理效果需要经过长时间离线化验才能获取。目前,针对不同来源的重金属废水的处理工况主要依赖于工艺人员的经验进行判断,往往难以准确的识别处理过程的正常工况和异常工况。当异常工况发生时,难以有效地降低废水中的重金属离子浓度,如未对处理过程进行及时调整,其处理结果会影响下一个工序的效果。而人工处理方式操作主观盲目性大,尤其是当重金属废水来源发生变化时,其工况识别准确性较低。
发明内容
本发明为了准确识别不同来源的重金属废水处理过程的异常工况,同时克服不同来源重金属离子各浓度数据带来的数据分布畸变导致模型失配带来的问题,本发明提出一种基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质,该方法能够适应不同来源重金属废水,准确识别处理过程的工况。
字典学习是一种高效的数据表示方式,该方法假设信号能进行稀疏表示,即表示为一个字典和其稀疏编码的乘积形式,通过学习并存储一个小数据量的字典,便能还原高维信号数据,这对于高维数据的处理是卓有成效的。由于能够减少高维数据的计算与存储负荷,字典学习受到了越来越多关注。
本技术方案通过对包含处理过程中pH值、电流密度、电导率、初始重金属浓度、流量等指标在内的历史数据进行建模,在对废水处理过程进行在线监测时,考虑到不同来源的废水其初始重金属浓度的时间序列上分布不同,使用历史数据得到的模型,其检测效果可能得不到保障。因此,对来源不确定的重金属废水的处理结果进行离线化验,根据化验结果选取其正常工况下的数据(即废水中重金属离子浓度得到有效降低的数据)和确定来源下的正常工况的历史数据进行融合,通过对历史模型进行迁移学习,实现不同来源下的废水处理过程的异常工况检测。
本发明的技术方案如下:
一方面,一种基于迁移学习的重金属废水处理过程异常工况智能化监测方法,包括:
1)利用历史采集的来源固定的重金属废水处理数据样本,构建来源固定的重金属废水处理数据样本的离线字典;
2)利用历史采集的来源未知的有效重金属废水处理数据样本,对离线字典进行迁移学习,获得来源未知的有效重金属废水处理数据样本对应的扩展字典;
3)利用扩展字典,计算来源未知下的有效重金属废水处理数据样本的重构误差,并基于所述重构误差,利用核密度估计获得重金属废水处理过程工况中的控制限;
4)计算待监测数据yt在扩展字典DTD下的重构误差,若计算得到的重构误差小于控制限,则认为当前重金属废水处理过程未出现异常,否则,则认为当前重金属废水处理过程出现异常。
进一步地,所述来源固定的重金属废水处理数据样本的离线字典的构建过程如下:
步骤1.1:历史样本采集;
利用传感器采集来源固定的重金属废水处理历史样本,来源固定的样本集合为YSD
Figure BDA0002417045890000021
yi表示来源固定的第i个重金属废水处理历史样本,1≤i≤Ns,每个样本包含m个废水指标{pH值、电流密度、电导率、初始重金属浓度、流量},Ns表示YSD中包含的样本个数;
步骤1.2:依据稀疏表示原理,利用字典D1和稀疏编码X表示YSD,构建离线字典学习的目标函数,并通过求解离线字典学习的目标函数,获得YSD对应的最优初始字典DSD,以及与DSD对应的稀疏编码XSD
Figure BDA0002417045890000022
Figure BDA0002417045890000023
其中,字典D1的初始取值为从历史样本集合YSD中随机选取的K个样本按列排列形成的矩阵,K=10*m,T为稀疏编码矩阵中的每一列非零元素个数设定值,
Figure BDA0002417045890000031
和||·||0分别表示二范数和零范数;xi表示X中的第i列。
字典D1的最终取值即为DSD,D1的初始值中每一列是一个随机选取的样本,且D1具有K列;稀疏编码X的最终取值即为XSD;XSD表示YSD在DSD下的稀疏编码,DSD中每一列表示一个字典原子;
T的取值通常设置为2;
进一步地,采用K-SVD方法对所述离线字典学习的目标函数进行求解,不断更新字典D1和稀疏编码X,直到获得YSD对应的最优初始字典DSD
采用K-SVD方法进行求解,具体地,从YSD随机选取K个样本作为字典DSD的初始值,利用正交匹配追踪算法求得稀疏编码XSD后;按列对字典D1进行更新,例如更新第k列字典原子时,可写成如下形式,
Figure BDA0002417045890000032
表示X中的第k行。
Figure BDA0002417045890000033
其中,Ek等于
Figure BDA0002417045890000034
定义集合
Figure BDA0002417045890000035
表示
Figure BDA0002417045890000036
中非零项所在索引的索引集,
Figure BDA0002417045890000037
表示
Figure BDA0002417045890000038
的第i个元素,N表示
Figure BDA0002417045890000039
中的元素个数。定义Ωk为N×|ωk|矩阵,它在(ωk(i),i)处的值为1,其余值均为0。将Ωk
Figure BDA00024170458900000310
Ek相乘后即可收缩原矩阵,
Figure BDA00024170458900000311
Figure BDA00024170458900000312
Figure BDA00024170458900000313
做奇异值分解后,
Figure BDA00024170458900000314
得到dk=u(:,1),
Figure BDA00024170458900000315
当逐列更新完之后,则利用正交匹配追踪算法交替更新稀疏编码XSD。迭代更新一定次数后,得到最优的初始字典DSD
进一步地,所述利用历史采集的来源未知的有效重金属废水处理数据样本,对离线字典进行迁移学习,获得扩展字典的过程如下:
利用传感器采集来源未知的有效重金属废水处理历史样本,来源未知的有效样本集合为YTD;依据稀疏表示原理,利用初始字典DSD和对应的稀疏编码X表示YTD,构建未知来源后重金属废水处理数据样本对应的稀疏编码的目标函数,通过迁移学习,求解获得来源为未知的有效样本集合YTD对应的最优稀疏编码Xp,再利用Xp获得对应的最优字典;
Figure BDA0002417045890000041
Figure BDA0002417045890000042
其中,Dp表示迁移学习过程中的插值字典,Dp的初始值表示来源固定的重金属废水处理数据样本的离线字典对应的最优初始字典DSD,即当p=0时,D0=DSD;T为稀疏编码矩阵中的非零元素个数设定值,
Figure BDA0002417045890000043
和||·||0分别表示二范数和零范数;xi表示X中的第i列。
进一步地,通过构建来源未知的重金属废水处理数据样本对应的字典的目标函数,求解来源为未知的有效样本集合YTD对应的最优字典;
Figure BDA0002417045890000044
其中,λ表示调节参数,取值范围为[1,10];D表示待求的字典,通过迭代求得的最终字典赋值给Dp+1,Dp+1表示求得的来源为未知的有效样本集合YTD对应的最优字典;
求解过程如下:
通过对
Figure BDA0002417045890000045
进行求导,获得插值字典的更新结果,
Figure BDA0002417045890000046
再通过对扩展字典的字典原子进行缩放,
Figure BDA0002417045890000047
不断重复该过程,直到
Figure BDA0002417045890000048
δ表示停止阈值,停止阈值设置为0.01;
Xp表示第p次迭代过程中,得到的稀疏编码,
Figure BDA0002417045890000049
是Xp的转置,E是单位矩阵,
Figure BDA00024170458900000410
Figure BDA00024170458900000411
表示扩展字典DTD的第1列和第K列。
进一步地,所述重构误差的计算是指待计算重构误差的样本采集值和该样本用扩展字典与对应稀疏编码的表示值之间的二范数计算。
进一步地,利用核密度估计获得重金属废水处理过程工况中的控制限,是指按照以下公式对来源未知的历史样本的重构误差计算核密度函数,且以核密度函数在设定置信度下的取值作为对应控制限:
Figure BDA0002417045890000051
其中,e是表示待拟合的来源未知的历史样本的重构误差的分布,ei表示第i个来源未知的历史样本的重构误差,H表示带宽矩阵,n表示历史样本的总数,K[·]表示核函数;
Figure BDA0002417045890000052
是指在给定带宽矩阵H下,通过历来源未知的历史样本ei拟合出的曲线。
在本实例中,核函数采用高斯核函数,带宽矩阵采用对角线矩阵,置信度设置为0.98;
进一步地,通过设置异常废水指标定位目标函数,依次设置方向选择向量,计算每个废水指标在扩展字典下的重构误差,直到异常样本上的异常幅值收敛,确定异常废水指标;
Figure BDA0002417045890000053
Figure BDA0002417045890000054
Figure BDA0002417045890000055
其中,yf为废水处理过程存在异常的样本,yi=yfifi,yi是对yf的第i维指标进行隔离,且其他维指标的取值保持不变的待重构样本,fi是yf中第i维指标的重构幅值;xri为yi在扩展字典下的稀梳编码,xri的初始值为yf在扩展字典下的稀疏编码;Rei是yf中第i个废水指标在扩展字典上的重构误差,DTD表示扩展字典;ξi表示方向选择向量,若该向量上第i个元素为1,表示此次选择了第i个指标,其他元素全为0,ξi∈Rm
Figure BDA0002417045890000056
Figure BDA0002417045890000057
分别是通过argmin目标函数对xri和fi进行优化后得到的结果值。
不断改变方向选择向量选取的指标,例如,第1维变量表示pH值,ξ1=[1,0,0,…0,0]T表示对pH值的一个选择向量,通过找到最小的重构误差对应的i值,确定异常指标;
另一方面,一种基于迁移学习的重金属废水处理过程异常工况智能化监测装置,其特征在于,包括:
离线字典建立模块,利用历史采集的来源固定的重金属废水处理数据样本,构建离线字典;
扩展字典建立模块,利用历史采集的来源未知的有效重金属废水处理数据样本,对离线字典进行迁移学习,构建扩展字典;
控制限生成模块,用于利用扩展字典计算所有历史样本的重构误差,并采用核密度估计方法,基于所有历史样本的重构误差计算重金属废水处理过程工况中的控制限;
工况异常判断模块,用于根据扩展字典计算在线获取的待监测数据的重构误差,并将待监测数据的重构误差与控制限比较,根据比较结果判断当前重金属废水处理过程工况是否异常。
另一方面,一种计算机存储介质,用于存储程序,所述程序被执行时,用于实现上述一种基于迁移学习的重金属废水处理过程异常工况智能化监测方法。
有益效果
本发明提出一种基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质,考虑到重金属废水处理过程中,废水来源各异,不同来源的废水其数据的分布差异较大,传统的方法难以准确识别重金属废水处理过程的工况;通过对不同来源的重金属废水处理过程数据融合,能够自动的实现不同来源的重金属废水处理过程异常工况智能识别;具体为利用来源固定的重金属废水处理过程的正常样本YSD、少量来源未知的重金属废水处理过程的正常样本YTD;首先通过对YSD进行学习得到其数据表示字典DSD,然后考虑到YSD和YTD分布不同,采用迁移学习的方法,将YTD的特征融入到字典学习过程,得到泛化能力更强的字典DTD。考虑到DTD既考虑了YSD的特征,同时又将YTD的分布差异性考虑到模型中,其具有更强的数据表示能力,因此能够在线准确的识别不同来源的重金属废水的异常工况。该方法无需过程先验知识,能自适应的适应废水处理系统中的不确定性因素,能够更加准确的检测过程中相关指标的变化,实现及时地检测与预警。
附图说明
图1为本发明实例的具体流程示意图。
具体实施方式
下面将结合附图和实例对本发明做进一步的说明。
如图1所示,本实例中,一种基于迁移学习的重金属废水处理过程异常工况智能化监测方法包括:
1)利用历史采集的来源固定的重金属废水处理数据样本,构建来源固定的重金属废水处理数据样本的离线字典;
所述来源固定的重金属废水处理数据样本的离线字典的构建过程如下:
步骤1.1:历史样本采集;
利用传感器采集来源固定的重金属废水处理历史样本,来源固定的样本集合为YSD
Figure BDA0002417045890000071
yi表示来源固定的第i个重金属废水处理历史样本,1≤i≤Ns,每个样本包含m个废水指标{pH值、电流密度、电导率、初始重金属浓度、流量},Ns表示YSD中包含的样本个数;
步骤1.2:依据稀疏表示原理,利用字典D1和稀疏编码X表示YSD,构建离线字典学习的目标函数,并通过求解离线字典学习的目标函数,获得YSD对应的最优初始字典DSD,以及与DSD对应的稀疏编码XSD
Figure BDA0002417045890000072
其中,字典D1的初始取值为从历史样本集合YSD中随机选取的K个样本按列排列形成的矩阵,K=10*m,T为稀疏编码矩阵中每一列的非零元素个数设定值,
Figure BDA0002417045890000073
和||·||0分别表示二范数和零范数;xi表示X中的第i列。
字典D1的最终取值即为DSD,D1的初始值中每一列是一个随机选取的样本,且D1具有K列;稀疏编码X的最终取值即为XSD;XSD表示YSD在DSD下的稀疏编码,DSD中每一列表示一个字典原子;
T的取值通常设置为2;
采用K-SVD方法对所述离线字典学习的目标函数进行求解,不断更新字典D1和稀疏编码X,直到获得YSD对应的最优初始字典DSD
采用K-SVD方法进行求解,具体地,从YSD随机选取K个样本作为字典DSD的初始值,利用正交匹配追踪算法求得稀疏编码XSD后;按列对字典D1进行更新,例如更新第k列字典原子时,可写成如下形式,
Figure BDA0002417045890000074
表示X中的第k行。
Figure BDA0002417045890000075
其中,Ek等于
Figure BDA0002417045890000081
定义集合
Figure BDA0002417045890000082
表示
Figure BDA0002417045890000083
中非零项所在索引的索引集。定义Ωk为N×|ωk|矩阵,它在(ωk(i),i)处的值为1,其余值均为0。将Ωk
Figure BDA0002417045890000084
Ek相乘后即可收缩原矩阵,
Figure BDA0002417045890000085
Figure BDA0002417045890000086
做奇异值分解后,
Figure BDA0002417045890000087
得到dk=u(:,1),
Figure BDA0002417045890000088
当逐列更新完之后,则利用正交匹配追踪算法交替更新稀疏编码XSD。迭代更新一定次数后,得到最优的初始字典DSD
2)利用历史采集的来源未知的有效重金属废水处理数据样本,对离线字典进行迁移学习,获得来源未知的有效重金属废水处理数据样本对应的扩展字典;
所述利用历史采集的来源未知的有效重金属废水处理数据样本,对离线字典进行迁移学习,获得扩展字典的过程如下:
利用传感器采集来源未知的有效重金属废水处理历史样本,来源未知的有效样本集合为YTD;依据稀疏表示原理,利用初始字典DSD和对应的稀疏编码X表示YTD,构建未知来源后重金属废水处理数据样本对应的稀疏编码的目标函数,通过迁移学习,求解获得来源为未知的有效样本集合YTD对应的最优稀疏编码Xp,再利用Xp获得对应的最优字典;
Figure BDA0002417045890000089
其中,Dp表示迁移学习过程中的插值字典,Dp的初始值表示来源固定的重金属废水处理数据样本的离线字典对应的最优初始字典DSD,即当p=0时,D0=DSD;T为稀疏编码矩阵中每一列的非零元素个数设定值,
Figure BDA00024170458900000810
和||·||0分别表示二范数和零范数;xi表示X中的第i列。
为了不断减小新分布数据的表示误差,并且同时保持插值字典模型的连续性,通过构建来源未知的重金属废水处理数据样本对应的字典的目标函数,求解来源为未知的有效样本集合YTD对应的最优字典;
Figure BDA00024170458900000811
其中,λ表示调节参数,取值范围为[1,10];D表示待求的字典,通过迭代求得的最终字典赋值给Dp+1,Dp+1表示求得的来源为未知的有效样本集合YTD对应的最优字典;
求解过程如下:
通过对
Figure BDA00024170458900000812
进行求导,获得插值字典的更新结果:
Figure BDA0002417045890000091
为了保证表示字典原子的L2范数等于1,再通过对扩展字典的字典原子进行缩放:
Figure BDA0002417045890000092
不断重复(3)-(6),直到
Figure BDA0002417045890000093
δ表示停止阈值,停止阈值设置为0.01;
Xp表示第p次迭代过程中,得到的稀疏编码,
Figure BDA0002417045890000094
是Xp的转置,E是单位矩阵,
Figure BDA0002417045890000095
Figure BDA0002417045890000096
表示扩展字典DTD的第1列和第K列。
Figure BDA0002417045890000097
3)利用扩展字典,计算来源未知下的有效重金属废水处理数据样本的重构误差,并基于所述重构误差,利用核密度估计获得重金属废水处理过程工况中的控制限;
利用核密度估计获得重金属废水处理过程工况中的控制限,是指按照以下公式对来源未知的历史样本的重构误差计算核密度函数,且以核密度函数在设定置信度下的取值作为对应控制限:
Figure BDA0002417045890000098
其中,e是表示待拟合的来源未知的历史样本的重构误差的分布,ei表示第i个来源未知的历史样本的重构误差,H表示带宽矩阵,n表示历史样本的总数,K[·]表示核函数;
Figure BDA0002417045890000099
是指在给定带宽矩阵H下,通过历来源未知的历史样本ei拟合出的曲线。
在本实例中,核函数采用高斯核函数,带宽矩阵采用对角线矩阵,置信度设置为0.98;
4)计算待监测数据yt在扩展字典DTD下的重构误差,若计算得到的重构误差小于控制限,则认为当前重金属废水处理过程未出现异常,否则,则认为当前重金属废水处理过程出现异常。
重构误差的计算是指待计算重构误差的样本采集值和该样本用扩展字典与对应稀疏编码的表示值之间的二范数计算。
基于得到的字典模型对废水处理过程工况进行异常检测及异常指标隔离:
通过设置异常废水指标定位目标函数,依次设置方向选择向量,计算每个废水指标在扩展字典下的重构误差,直到异常样本上的异常幅值收敛,确定异常废水指标;
Figure BDA0002417045890000101
Figure BDA0002417045890000102
Figure BDA0002417045890000103
其中,yf为废水处理过程存在异常的样本,yi=yfifi,yi是对yf的第i维指标进行隔离,且其他维指标的取值保持不变的待重构样本,fi是yf中第i维指标的重构幅值;xri为yi在扩展字典下的稀梳编码,xri的初始值为yf在扩展字典下的稀疏编码;Rei是yf中第i个废水指标在扩展字典上的重构误差,DTD表示扩展字典;ξi表示方向选择向量,若该向量上第i个元素为1,表示此次选择了第i个指标,其他元素全为0,ξi∈Rm
Figure BDA0002417045890000104
Figure BDA0002417045890000105
分别是通过argmin目标函数对xri和fi进行优化后得到的结果值。
不断改变方向选择向量选取的指标,例如,第1维变量表示pH值,ξ1=[1,0,0,…0,0]T表示对pH值的一个选择向量,通过找到最小的重构误差对应的i值,确定异常指标;
通过异常隔离能够找出废水处理过程中可能的问题,例如如果发现异常是由pH导致,未来可以通过改变药剂添加量实现pH的稳定,进而保证废水处理过程的正常运行。
基于上述方法,本发明实施例还提供一种基于迁移学习的重金属废水处理过程异常工况智能化监测装置,包括:
离线字典建立模块,利用历史采集的来源固定的重金属废水处理数据样本,构建离线字典;
扩展字典建立模块,利用历史采集的来源未知的有效重金属废水处理数据样本,对离线字典进行迁移学习,构建扩展字典;
控制限生成模块,用于利用扩展字典计算所有历史样本的重构误差,并采用核密度估计方法,基于所有历史样本的重构误差计算重金属废水处理过程工况中的控制限;
工业系统异常判断模块,用于根据扩展字典计算在线获取的待监测数据的重构误差,并将待监测数据的重构误差与控制限比较,根据比较结果判断当前重金属废水处理过程工况是否异常。
应当理解,本发明各个实施例中的功能单元模块可以集中在一个处理单元中,也可以是各个单元模块单独物理存在,也可以是两个或两个以上的单元模块集成在一个单元模块中,可以采用硬件或软件的形式来实现。
本发明实施例还提供一种计算机存储介质,用于存储程序,所述程序被执行时,用于实现一种基于迁移学习的重金属废水处理过程异常工况智能化监测方法,其有益效果参见方法部分的有益效果,在此不再赘述。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。

Claims (10)

1.一种基于迁移学习的重金属废水处理过程异常工况智能化监测方法,其特征在于,包括:
1)利用历史采集的来源固定的重金属废水处理数据样本,构建来源固定的重金属废水处理数据样本的离线字典;
2)利用历史采集的来源未知的有效重金属废水处理数据样本,对离线字典进行迁移学习,获得来源未知的有效重金属废水处理数据样本对应的扩展字典;
3)利用扩展字典,计算来源未知下的有效重金属废水处理数据样本的重构误差,并基于所述重构误差,利用核密度估计获得重金属废水处理过程工况中的控制限;
4)计算待监测数据yt在扩展字典DTD下的重构误差,若计算得到的重构误差小于控制限,则认为当前重金属废水处理过程未出现异常,否则,则认为当前重金属废水处理过程出现异常。
2.根据权利要求1所述的方法,其特征在于,所述来源固定的重金属废水处理数据样本的离线字典的构建过程如下:
步骤1.1:历史样本采集;
利用传感器采集来源固定的重金属废水处理历史样本,来源固定的样本集合为YSD
Figure FDA0002417045880000011
yi表示来源固定的第i个重金属废水处理历史样本,1≤i≤Ns,每个样本包含m个废水指标{pH值、电流密度、电导率、初始重金属浓度、流量},Ns表示YSD中包含的样本个数;
步骤1.2:依据稀疏表示原理,利用字典D1和稀疏编码X表示YSD,构建离线字典学习的目标函数,并通过求解离线字典学习的目标函数,获得YSD对应的最优初始字典DSD,以及与DSD对应的稀疏编码XSD
Figure FDA0002417045880000012
Figure FDA0002417045880000013
其中,字典D1的初始取值为从历史样本集合YSD中随机选取的K个样本按列排列形成的矩阵,K=10*m,T为稀疏编码矩阵中的每一列非零元素个数设定值,
Figure FDA0002417045880000014
和||·||0分别表示二范数和零范数;xi表示X中的第i列。
3.根据权利要求2所述的方法,其特征在于,采用K-SVD方法对所述离线字典学习的目标函数进行求解,不断更新字典D1和稀疏编码X,直到获得YSD对应的最优初始字典DSD
4.根据权利要求1所述的方法,其特征在于,所述利用历史采集的来源未知的有效重金属废水处理数据样本,对离线字典进行迁移学习,获得扩展字典的过程如下:
利用传感器采集来源未知的有效重金属废水处理历史样本,来源未知的有效样本集合为YTD;依据稀疏表示原理,利用初始字典DSD和对应的稀疏编码X表示YTD,构建未知来源后重金属废水处理数据样本对应的稀疏编码的目标函数,通过迁移学习,求解获得来源为未知的有效样本集合YTD对应的最优稀疏编码Xp,再利用Xp获得对应的最优字典;
Figure FDA0002417045880000021
Figure FDA0002417045880000022
其中,Dp表示迁移学习过程中的插值字典,Dp的初始值表示来源固定的重金属废水处理数据样本的离线字典对应的最优初始字典DSD,即当p=0时,D0=DSD;T为稀疏编码矩阵中的非零元素个数设定值,
Figure FDA0002417045880000023
和||·||0分别表示二范数和零范数;xi表示X中的第i列。
5.根据权利要求4所述的方法,其特征在于,通过构建来源未知的重金属废水处理数据样本对应的字典的目标函数,求解来源为未知的有效样本集合YTD对应的最优字典;
Figure FDA0002417045880000024
其中,λ表示调节参数,取值范围为[1,10];D表示待求的字典,通过迭代求得的最终字典赋值给Dp+1,Dp+1表示求得的来源为未知的有效样本集合YTD对应的最优字典;
求解过程如下:
通过对
Figure FDA0002417045880000025
求导,获得插值字典的更新结果,
Figure FDA0002417045880000026
再通过对扩展字典的字典原子进行缩放,
Figure FDA0002417045880000027
不断重复该过程,直到
Figure FDA0002417045880000028
δ表示停止阈值,停止阈值设置为0.01;
Xp表示第p次迭代过程中,得到的稀疏编码,
Figure FDA0002417045880000029
是Xp的转置,E是单位矩阵,
Figure FDA00024170458800000210
Figure FDA00024170458800000211
表示扩展字典DTD的第1列和第K列。
6.根据权利要求1所述的方法,其特征在于,所述重构误差的计算是指待计算重构误差的样本采集值和该样本用扩展字典与对应稀疏编码的表示值之间的二范数计算。
7.根据权利要求1所述的方法,其特征在于,利用核密度估计获得重金属废水处理过程工况中的控制限,是指按照以下公式对来源未知的历史样本的重构误差计算核密度函数,且以核密度函数在设定置信度下的取值作为对应控制限:
Figure FDA0002417045880000031
其中,e是表示待拟合的来源未知的历史样本的重构误差的分布,ei表示第i个来源未知的历史样本的重构误差,H表示带宽矩阵,n表示历史样本的总数,K[·]表示核函数;
Figure FDA0002417045880000032
是指在给定带宽矩阵H下,通过历来源未知的历史样本ei拟合出的曲线。
8.根据权利要求1所述的方法,其特征在于,通过设置异常废水指标定位目标函数,依次设置方向选择向量,计算每个废水指标在扩展字典下的重构误差,直到异常样本上的异常幅值收敛,确定异常废水指标;
Figure FDA0002417045880000033
Figure FDA0002417045880000034
Figure FDA0002417045880000035
其中,yf为废水处理过程存在异常的样本,yi=yfifi,yi是对yf的第i维指标进行隔离,且其他维指标的取值保持不变的待重构样本,fi是yf中第i维指标的重构幅值;xri为yi在扩展字典下的稀梳编码,xri的初始值为yf在扩展字典下的稀疏编码;Rei是yf中第i个废水指标在扩展字典上的重构误差,DTD表示扩展字典;ξi表示方向选择向量,若该向量上第i个元素为1,表示此次选择了第i个指标,其他元素全为0,ξi∈Rm
Figure FDA0002417045880000036
Figure FDA0002417045880000037
分别是通过argmin目标函数对xri和fi进行优化后得到的结果值。
9.一种基于迁移学习的重金属废水处理过程异常工况智能化监测装置,其特征在于,包括:
离线字典建立模块,利用历史采集的来源固定的重金属废水处理数据样本,构建离线字典;
扩展字典建立模块,利用历史采集的来源未知的有效重金属废水处理数据样本,对离线字典进行迁移学习,构建扩展字典;
控制限生成模块,用于利用扩展字典计算所有历史样本的重构误差,并采用核密度估计方法,基于所有历史样本的重构误差计算重金属废水处理过程工况中的控制限;
工业系统异常判断模块,用于根据扩展字典计算在线获取的待监测数据的重构误差,并将待监测数据的重构误差与控制限比较,根据比较结果判断当前重金属废水处理过程工况是否异常。
10.一种计算机存储介质,其特征在于,用于存储程序,所述程序被执行时,用于实现如权利要求1-8任一所述的方法。
CN202010194359.2A 2020-03-19 2020-03-19 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质 Active CN111427265B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202010194359.2A CN111427265B (zh) 2020-03-19 2020-03-19 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质
PCT/CN2021/077910 WO2021185044A1 (zh) 2020-03-19 2021-02-25 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质
US17/908,922 US20230089156A1 (en) 2020-03-19 2021-02-25 Intelligent monitoring method and apparatus for abnormal working conditions in heavy metal wastewater treatment process based on transfer learning and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010194359.2A CN111427265B (zh) 2020-03-19 2020-03-19 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质

Publications (2)

Publication Number Publication Date
CN111427265A true CN111427265A (zh) 2020-07-17
CN111427265B CN111427265B (zh) 2021-03-16

Family

ID=71548055

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010194359.2A Active CN111427265B (zh) 2020-03-19 2020-03-19 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质

Country Status (3)

Country Link
US (1) US20230089156A1 (zh)
CN (1) CN111427265B (zh)
WO (1) WO2021185044A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021185044A1 (zh) * 2020-03-19 2021-09-23 中南大学 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质
CN116088307A (zh) * 2022-12-28 2023-05-09 中南大学 基于误差触发自适应稀疏辨识的多工况工业过程预测控制方法、装置、设备及介质

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114115098A (zh) * 2021-11-02 2022-03-01 浙江尔格科技股份有限公司 一种基于字典学习的冷却系统性能预警方法
CN114410984B (zh) * 2022-01-25 2023-11-17 湖南株冶有色金属有限公司 一种湿法炼锌浸出过程异常工况的控制方法
CN115129029B (zh) * 2022-07-01 2024-07-16 中南大学 基于子领域适配字典学习的工业系统故障诊断方法和系统
CN116125922B (zh) * 2023-01-09 2024-06-21 中南大学 一种基于平行式字典学习的复杂工业过程监测方法和系统
CN116125923B (zh) * 2023-01-09 2024-06-21 中南大学 基于混合变量字典学习的混杂工业过程监测方法和系统
CN116577671B (zh) * 2023-07-12 2023-09-29 中国华能集团清洁能源技术研究院有限公司 电池系统异常检测方法及装置
CN116776748B (zh) * 2023-08-18 2023-11-03 中国人民解放军国防科技大学 喉栓式变推力发动机喉栓喷管构型设计知识迁移优化方法
CN117872761A (zh) * 2024-01-12 2024-04-12 中南大学 面向多变工况的字典监测模型构建、智能监测方法及系统

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012061813A1 (en) * 2010-11-05 2012-05-10 Georgetown University Event detection, workflow analysis, and reporting system and method
CN104462019A (zh) * 2014-12-18 2015-03-25 江西理工大学 一种稀疏表示下支持向量机核函数选择方法及其应用
CN105931139A (zh) * 2016-05-05 2016-09-07 中科智水(北京)科技有限公司 水务数据标准化管理系统
CN105976028A (zh) * 2016-05-11 2016-09-28 深圳市开天源自动化工程有限公司 一种预测a2o污水处理过程中出水cod浓度的方法
US20170027168A1 (en) * 2015-07-27 2017-02-02 Stephan HEATH Methods, products, and systems relating to making, providing, and using nanocrystalline (nc) products comprising nanocrystalline cellulose (ncc), nanocrystalline (nc) polymers and/or nanocrystalline (nc) plastics or other nanocrystals of cellulose composites or structures, in combination with other materials
CN108562709A (zh) * 2018-04-25 2018-09-21 重庆工商大学 一种基于卷积自编码器极限学习机的污水处理系统水质监测预警方法
CN108821435A (zh) * 2018-07-11 2018-11-16 石家庄市桥西污水处理厂 一种污水处理中的溶氧控制方法
CN109243545A (zh) * 2018-11-21 2019-01-18 哈尔滨工业大学 一种基于机器学习的难降解有机污水处理效果预测方法
CN110580488A (zh) * 2018-06-08 2019-12-17 中南大学 基于字典学习的多工况工业监测方法、装置、设备及介质
CN110826611A (zh) * 2019-10-30 2020-02-21 华南理工大学 基于多个元分类器加权集成的stacking污水处理故障诊断方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104199441B (zh) * 2014-08-22 2017-03-01 清华大学 基于稀疏贡献图的高炉多工况故障分离方法及系统
CN104182642B (zh) * 2014-08-28 2017-06-09 清华大学 一种基于稀疏表示的故障检测方法
US20160358075A1 (en) * 2015-06-08 2016-12-08 The Regents Of The University Of Michigan System for implementing a sparse coding algorithm
CN110579967B (zh) * 2019-09-23 2020-06-30 中南大学 基于同时降维和字典学习的过程监控方法
CN111427265B (zh) * 2020-03-19 2021-03-16 中南大学 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012061813A1 (en) * 2010-11-05 2012-05-10 Georgetown University Event detection, workflow analysis, and reporting system and method
CN104462019A (zh) * 2014-12-18 2015-03-25 江西理工大学 一种稀疏表示下支持向量机核函数选择方法及其应用
US20170027168A1 (en) * 2015-07-27 2017-02-02 Stephan HEATH Methods, products, and systems relating to making, providing, and using nanocrystalline (nc) products comprising nanocrystalline cellulose (ncc), nanocrystalline (nc) polymers and/or nanocrystalline (nc) plastics or other nanocrystals of cellulose composites or structures, in combination with other materials
CN105931139A (zh) * 2016-05-05 2016-09-07 中科智水(北京)科技有限公司 水务数据标准化管理系统
CN105976028A (zh) * 2016-05-11 2016-09-28 深圳市开天源自动化工程有限公司 一种预测a2o污水处理过程中出水cod浓度的方法
CN108562709A (zh) * 2018-04-25 2018-09-21 重庆工商大学 一种基于卷积自编码器极限学习机的污水处理系统水质监测预警方法
CN110580488A (zh) * 2018-06-08 2019-12-17 中南大学 基于字典学习的多工况工业监测方法、装置、设备及介质
CN108821435A (zh) * 2018-07-11 2018-11-16 石家庄市桥西污水处理厂 一种污水处理中的溶氧控制方法
CN109243545A (zh) * 2018-11-21 2019-01-18 哈尔滨工业大学 一种基于机器学习的难降解有机污水处理效果预测方法
CN110826611A (zh) * 2019-10-30 2020-02-21 华南理工大学 基于多个元分类器加权集成的stacking污水处理故障诊断方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CHEN ZHI-MING: "Wastewater Treatment Prediction Based on Chaos-GA Optimized LS-SVM", 《2011 CHINESE CONTROL AND DECISION CONFERENCE》 *
HUI WEN: "Abnormal Event Detection via Adaptive Cascade Dictionary Learning", 《ICIP 2015》 *
秦美华: "基于STA-K 均值聚类的电化学废水处理过程离子浓度软测量", 《化工学报》 *
韩红桂: "城市污水处理过程异常工况识别和抑制研究", 《自动化学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021185044A1 (zh) * 2020-03-19 2021-09-23 中南大学 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质
CN116088307A (zh) * 2022-12-28 2023-05-09 中南大学 基于误差触发自适应稀疏辨识的多工况工业过程预测控制方法、装置、设备及介质
CN116088307B (zh) * 2022-12-28 2024-01-30 中南大学 基于误差触发自适应稀疏辨识的多工况工业过程预测控制方法、装置、设备及介质

Also Published As

Publication number Publication date
CN111427265B (zh) 2021-03-16
WO2021185044A1 (zh) 2021-09-23
US20230089156A1 (en) 2023-03-23

Similar Documents

Publication Publication Date Title
CN111427265B (zh) 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质
Chai et al. Enhanced random forest with concurrent analysis of static and dynamic nodes for industrial fault classification
WO2020049087A1 (en) Computer-implemented method, computer program product and system for anomaly detection and/or predictive maintenance
CN111982302A (zh) 一种带噪声滤波和环境温度补偿的测温方法
Tong et al. Statistical process monitoring based on a multi-manifold projection algorithm
CN112101400A (zh) 工业控制系统异常检测方法、设备和服务器、存储介质
CN112101554A (zh) 异常检测方法及装置、设备、计算机可读存储介质
Sotomayor et al. Multivariable identification of an activated sludge process with subspace-based algorithms
Ullah et al. Comparison of tree-based model with deep learning model in predicting effluent pH and concentration by capacitive deionization
Zou et al. Two‐level multi‐block operating performance optimality assessment for plant‐wide processes
Gao et al. Particle swarm optimization based RVM classifier for non-linear circuit fault diagnosis
CN114331214A (zh) 基于强化学习的域自适应轴承声纹故障诊断方法及系统
Jiang et al. A timeseries supervised learning framework for fault prediction in chiller systems
CN114266284A (zh) 开关柜绝缘缺陷类型检测方法、装置、设备和程序产品
Parvathy et al. Hybrid machine learning based false data injection attack detection and mitigation model for waste water treatment plant
Na et al. Failure detection using a fuzzy neural network with an automatic input selection algorithm
Tedesco et al. A scalable deep learning-based approach for anomaly detection in semiconductor manufacturing
Souza et al. Variable selection based on mutual information for soft sensors applications
Hu et al. Fault Diagnosis Using Deep Learning for Wastewater Treatment Processes
CN113791275B (zh) 一种单相谐波数据丢失的修复方法及系统
CN113065606B (zh) 一种基于轻量级深度学习的异常点位检测方法及系统
CN116736173B (zh) 一种储能电池模型构建和储能电池状态判定方法、装置
KR102383411B1 (ko) 동적 분별 신호 학습을 통한 산업로봇의 적응적 예지보전 장치 및 방법
CN116125922B (zh) 一种基于平行式字典学习的复杂工业过程监测方法和系统
Shiri et al. Robust switching Kalman filter for diagnostics of long-term condition monitoring data in the presence of non-Gaussian noise

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