CN113011661B - 一种用于河流生态修复的曝气控制系统及其控制方法 - Google Patents

一种用于河流生态修复的曝气控制系统及其控制方法 Download PDF

Info

Publication number
CN113011661B
CN113011661B CN202110315841.1A CN202110315841A CN113011661B CN 113011661 B CN113011661 B CN 113011661B CN 202110315841 A CN202110315841 A CN 202110315841A CN 113011661 B CN113011661 B CN 113011661B
Authority
CN
China
Prior art keywords
river
aeration
data
value
parameters
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
CN202110315841.1A
Other languages
English (en)
Other versions
CN113011661A (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN202110315841.1A priority Critical patent/CN113011661B/zh
Publication of CN113011661A publication Critical patent/CN113011661A/zh
Application granted granted Critical
Publication of CN113011661B publication Critical patent/CN113011661B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • C—CHEMISTRY; METALLURGY
    • C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F7/00—Aeration of stretches of water
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/21—Design, administration or maintenance of databases
    • G06F16/215—Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067—Enterprise or organisation modelling
    • C—CHEMISTRY; METALLURGY
    • C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2103/00—Nature of the water, waste water, sewage or sludge to be treated
    • C02F2103/007—Contaminated open waterways, rivers, lakes or ponds
    • C—CHEMISTRY; METALLURGY
    • C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
    • C02F2209/00—Controlling or monitoring parameters in water treatment
    • C02F2209/005—Processes using a programmable logic controller [PLC]
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Educational Administration (AREA)
  • General Engineering & Computer Science (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • Biomedical Technology (AREA)
  • Computational Linguistics (AREA)
  • Artificial Intelligence (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Databases & Information Systems (AREA)
  • Hydrology & Water Resources (AREA)
  • Environmental & Geological Engineering (AREA)
  • Water Supply & Treatment (AREA)
  • Chemical & Material Sciences (AREA)
  • Organic Chemistry (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明涉及一种用于河流生态修复的曝气控制系统及其控制方法，该系统包括感知模块、控制模块和曝气模块，感知模块用于采集河流概况、河流水质、河流水动力和河流生物四个维度的河流基础数据，控制模块用于分析处理感知模块采集到的数据，将其换算成控制曝气模块的参数，并给曝气模块发送控制信号。该方法包括获取四个维度的数据，控制模块对输入数据经过预处理后，输入神经网络模型，计算出通过曝气修复河流的目标值；启动曝气模块，实时监控河流曝气情况，维持河流的溶解氧浓度在预期水平；直至河流的现状达到神经网络模型得到的曝气修复河流的目标值，停止曝气并通知管理人员。本发明在保证曝气效果的前提下，降低曝气设备的使用能耗。

Description

一种用于河流生态修复的曝气控制系统及其控制方法
技术领域
本发明涉及一种用于河流生态修复的曝气控制系统及其控制方法，属于河流生态修复的技术领域。
背景技术
河流水体中的溶解氧的浓度，是反应河流污染情况的重要指标。在河流污染治理过程中，首先需要对河流水体的溶解氧含量进行改变。河流的溶解氧浓度是河流耗氧过程和复氧过程共同作用的结果。生态健康的河流，耗氧速率与复氧速率处于平衡状态，河流污染之后，水体中过量的有机物导致河流的耗氧速率增加，溶解氧浓度不断下降，河流的氧平衡被打破，河流出现黑臭等环境问题。随着有机物的减少，耗氧速率逐渐下降，当复氧速率大于耗氧速率时，溶解氧浓度开始回升，河流恢复污染前的状态，这是河流自然的自净过程。然而，由于外源有机物持续输入超出了河流本身的自净能力，需要人为增加河流的复氧速率，曝气技术成为河流污染治理与生态修复工程中提高复氧速率的常用技术。
在河流污染治理与生态修复的工程实践中，采用的曝气控制方式，多是连续曝气或间歇曝气，曝气量多是凭操作人员的经验人为设定，存在曝气不足或过量的问题，会造成能源浪费和运行成本增加。近年，出现了基于实时溶解氧浓度的曝气控制模型，以期解决这个问题。但是，由于难以定义健康河流的溶解氧浓度，这类控制方式取得的实际效果也并不理想。
发明内容
为了解决上述技术问题，本发明提供一种用于河流生态修复的曝气控制系统及其控制方法，其具体技术方案如下：
一种用于河流生态修复的曝气控制系统，包括感知模块、控制模块和曝气模块，所述感知模块用于采集河流概况、河流水质、河流水动力和河流生物四个维度的河流基础数据，所述控制模块用于分析处理感知模块采集到的数据，将其换算成控制曝气模块的参数，并给曝气模块发送控制信号。
进一步的，所述河流概况的数据包括河流宽度、河流深度、水面坡降、河流中弘线长度、河流起止断面距离、城市连接廊道数量、最大的廊道数量、水系节点个数、污染源情况、周围土地利用类型中的一种或多种；
所述河流水质的数据包括水温、溶解氧浓度、氧化还原电位、浊度、氨氮浓度、硝酸盐氮浓度、总氮、总磷、耗氧速率、生物需氧量、化学需氧量中的一种或多种；
所述河流水动力的数据包括河流的流速、流量、剪切力中的一种或多种；
所述河流生物的数据包括微生物、原生动物、后生动物、鱼类、藻类、浮游植物、沉水植物中的一种或多种。
一种用于河流生态修复的曝气系统的控制方法，包括以下步骤：
步骤一：获取河流概况、河流水质、河流水动力和河流生物四个维度的数据，通过人机交互的方式输入控制模块，控制模块对输入数据经过预处理后，输入神经网络模型，计算出通过曝气修复河流的目标值；
步骤二：启动曝气模块，感知模块和控制模块同步开始工作，根据曝气模块启动前输入控制模块的数据，以及感知模块实时监测的数据，每隔六小时计算一次河流最佳曝气量，按公式(1)(2)计算曝气量，向曝气模块发送控制信号，控制曝气的强度和时间；根据感知模块传输给控制模块的数据，实时监控河流曝气情况，维持河流的溶解氧浓度在预期水平；
Figure BDA0002991113060000021
e(k)＝DOmeasure(k)-DOset/default(k) (2)
其中，Q(k)表示曝气量，Kp表示e(k)项的系数，Ki表示
Figure BDA0002991113060000022
项的系数，Kd表示[e(k)-e(k-1)]项的系数，k表示k时刻，e(k)表示k时刻测量值与设置/预设溶解氧浓度值的偏差，e(k-1)表示k-1时刻测量值与设置/预设值的偏差，DOmeasure(k)表示测量得到的k时刻溶解氧浓度，DOset/default(k)表示k时刻曝气器控制模块溶解氧的输入值或预设值；
步骤三：循环步骤二，直至感知模块监测的河流的现状达到神经网络模型得到的曝气修复河流的目标值，停止曝气并通知管理人员。
进一步的，所述步骤一中的神经网络模型，建立方法为：
步骤1：收集基础数据，选取有代表性的河流，收集包括河流曝气修复前后常规理化指标、生物指标，以及河流宽度、深度、功能、周围土地利用类型、污染源情况，至少80组数据；
步骤2：数据预处理，首先剔除明显异常和错误的数据，然后基于公式(3)做共线性分析，剔除共线性高的数据，再用排序算法做进一步降维，将输入指标降低到10个以内；最后，将处理好的数据划分为训练样本集和测试样本集；(随机按7：3划分训练样本集和测试样本集)
Figure BDA0002991113060000023
其中，R2为决定系数；VIF≥20剔除，VIF＜20保留；
步骤3：训练神经网络模型，将预处理后的结果作为输入层，将溶解氧浓度、氨氮浓度、化学需氧量、河流耗氧速率、自养生物占比作为输出层，输入层神经元数量与预处理后的指标个数相同，输出层神经元数量为5个；隐含层数量与神经元个数根据神经网络模型的质量进行调整，默认为两层；完成建模后，训练神经网络模型，根据输出结果，评估模型的损失函数，按公式(4)计算平均平方误差MSE表示神经网络模型的损失，直至图谱表现为水平直线，判定为符合要求；
Figure BDA0002991113060000031
其中，h(xi)代表预测的第i个值，yi代表实际的第i个值；m1为测试集样本数量；
步骤4：使用测试样本集测试模型的准确性，保证神经网络模型拟合良好，不存在欠拟合和过拟合的现象，准确率需达到90％以上。
进一步的，四个维度的河流基础数据处理过程如下：
(1)基于河流基础数据，计算河流复氧速率、弯曲度和连通度，
(1.1)复氧速率计算方法：对于水深较深、流速缓慢的河流，使用公式(5)(7)计算；对于水深较浅、流速较大(靠实验和经验确定)的河流，使用公式(6)(7)计算，
Figure BDA0002991113060000032
Figure BDA0002991113060000033
式中，k2为复氧系数；D为氧在水体中的分子扩散系数；u为平均流速；H为河流水深；J为水面坡降；κ为卡门常数，一般取0.4；g为重力加速度；
Figure BDA0002991113060000034
式中，O为溶解氧浓度；Os为饱和溶解氧浓度，t为时间；
使用公式(8)计算河流弯曲度：
Ie＝S/L (8)
式中，Ie为河流弯曲度；S为河流中泓线长度；L为河流起止断面距离；
(1.2)河流连通度计算方法：水系的环通度α使用公式(9)计算；水系节点率β使用公式(10)计算；水系网络连通度γ使用公式(11)计算，
Figure BDA0002991113060000041
Figure BDA0002991113060000042
Figure BDA0002991113060000043
式中，α为水系的环通度；M为城市连接的廊道数量；V为水系节点个数；β为水系节点率；Mmax为最大的廊道数量；
(2)层次分析法结合熵权法，对水质指标和水动力指标降维；
(2.1)层次分析法得到主观赋权的权重值：
2.1.a.建立层次结构模型，根据收集的水质和水动力参数，分目标层、准则层和方案层建立层次结构模型；
2.1.b.建立判断矩阵，根据步骤a建立的层次结构模型构建成对比较矩阵(判断矩阵)，包括分别基于目标层、准则层和方案层的判断矩阵，并计算特征向量和特征值；
2.1.c.进行层次单排序及其一次性检验，分别计算准则层对目标层和方案层对准则层的权重排序，分别做一次性检验；
2.1.d.进行层次总排序及其一次性检验，计算方案层对目标的权重排序及其一致性检验；
(2.2)熵权法得到客观赋权的权重值：
2.2.a.构建原始数据矩阵X＝(xij)m×n，并归一化，
2.2.b.使用公式(12)计算熵值
Figure BDA0002991113060000044
式中，ej为j指标的熵值；m2为曝气河流样本数量；n为指标数量；Pij为xij归一化后的矩阵；
2.2.c.使用公式(13)计算权重
Figure BDA0002991113060000045
式中，jj为j指标的权重；
2.2.d.使用公式(14)归一化处理
Figure BDA0002991113060000046
式中，βj为j指标的客观权重值；
利用公式(15)计算综合权重值
Figure BDA0002991113060000051
式中，w为由乘法归一组合赋权法得到的权重值；αj为主观赋权得到的权重值；βj为客观赋权得到的权重值；
(3)将生物数据划分为自养生物和异养生物，分别计算生物完整性指数；计算方法为：
3.a.选择受损点和参照点，选择生态修复前的数据为受损点，生态修复后的数据为参照点；
3.b.依据敏感性分析结果剔除参数，对生物参数做敏感性分析，其中IQ小于2的参数，不再做进一步分析，对于IQ值大于或等于2的参数，保留，并进入步骤3.c；
3.c.对剔除共线性参数，敏感性分析筛选后的参数进一步做Pearson相关分析和共线性
分析，按照|r|>0.75或VIF≥20的标准，剔除共线性的参数；
Figure BDA0002991113060000052
Figure BDA0002991113060000053
式中，r为相关系数；VIF为方差膨胀因子；R2为决定系数；m2为曝气河流样本数量；xi,yi为成对的样本；
3.d.按照三分法，首先计算构成生物完整性指数值(IBI值)的参数在参照点的统计分布，拟定其评价指标的评分标准，据此将每个样点的参数的评价值实测值分别进行评分，然后将选出的参数得分进行累加，最后得到每个样点的IBI总分；
3.e.按照比值法，首先计算生物参数在全部样点中的95％或5％分位数的值；最后将选出的参数的分值进行累加，即得出各样点的IBI总分值。
进一步的，所述步骤二所需的曝气量，即维持河流最佳需氧量的计算方法为：
步骤(1)：构建河流食物网模型，以生物与营养物质之间的相互作用、生物间相互作用为基础，构建河流食物网模型；
步骤(2)：基于动力学参数和实际收集的数据，模拟河流生态系统的生物量变化和营养物质变化，通过比对实际数据和模拟数据，率定动力学参数，使河流食物网模型的准确度符合需求；
步骤(3)：对河流食物网模型的初始输入值做敏感性分析，找到控制河流生态系统的关键变量；
步骤(4)：以当前溶解氧浓度为基准，以当前温度下的饱和溶解氧浓度为上限，以0.2mg/L为间隔，生成溶解氧序列；将生成的溶解氧序列依次输入河流食物网模型，其余初始值使用当前时间的在线监测值，动力学参数使用经率定的参数；根据步骤(3)敏感性分析识别出的关键变量，以关键变量最优为目标，得到河流的最佳需氧量。
进一步的，所述率定动力学参数的过程是：先使用文献或实验测得的动力学参数，输入模型进行模拟，得到预测值与实际值比较，|预测值-实际值|/实际值可以表示模型的准确度。根据模型的需求不同，选取不同的准确度阈值。(例如要求准确度为80％，则根据准确度低于80％的预测值，修改该指标对应的动力学参数，直至所有指标的预测值达到80％。)
进一步的，找到控制河流生态系统的关键变量的具体过程为：根据公式(18)计算不同变量的标准局部敏感系数
Figure BDA0002991113060000061
式中，Si,j为标准局部敏感系数；αi为模型参数相对变化量；Vj为模型输出变量的相对变化量；
Figure BDA0002991113060000062
为参数默认值改变1％；
标准局部敏感系数Si,j大于1表示变量敏感，即关键变量。
本发明的有益效果是：
1、能根据河流现状，科学设定河流曝气修复的目标和策略，准确判断河流是否需要曝气，避免出现停止曝气后，河流水质出现“反弹”。
2、能精细化控制河流修复的整个曝气过程，通过计算最佳需氧量，在保证曝气效果的前提下，降低曝气设备的使用能耗，达到节约能源，降低运行成本的目的。
3、本发明通过神经网络模型建模和预训练处理，在应用与河流时，只需采集河道的思维数据，并输入神经网络模型，便可计算出河流的最佳需氧量以及曝气量，做到科学治理河道水体。
附图说明
图1是本发明的控制流程图，
图2是神经网络平均平方误差图，
图3是动力学模型状态变量结构图，
图4是本发明数据流向概图，
图5是本发明曝气目标模型(基于神经网络算法)的建模过程图，
图6是本发明曝气目标模型模型使用图，
图7是本发明最佳曝气量模型的运算过程概图。
具体实施方式
现在结合附图对本发明作进一步详细的说明。这些附图均为简化的示意图，仅以示意方式说明本发明的基本结构，因此其仅显示与本发明有关的构成。
如图1所示，本发明系统包括三大板块，分别为感知模块、控制模块、曝气模块。
感知模块主要是指用于采集河流基础数据的各种设备，采集的数据为：河流概况、河流水质、河流水动力和河流生物四个维度的河流基础数据，包括理化指标监测、生物监测和河流耗氧速率在线监测，四个维度的数据，有两个用途，分别是计算最佳曝气量和确定能否停止曝气。采集的数据输入控制模块，数据输入方式可以是人机交互界面输入、传感器中的任意一种或多种。
其中生物监测的监测方式可以是在线监测、现场实测、实验室监测中的任意一种或多种。监测指标可以是微生物、原生动物、后生动物、鱼类、藻类、浮游植物、沉水植物中的一种或多种。
控制模块包括建模、预处理、训练，形成稳定的神经网络模型，然后使用神经网络模型，处理河流的四个维度数据，得出河流的最佳需氧量和曝气量。
曝气模块可选用现有的曝气设备，或者可选专利号为201910914913.7公开的一种用于河流修复的曝气系统，对河流进行曝气，控制模块与其控制连接，控制其曝气流速、时间、曝气量等参数。
下面介绍本发明系统的控制方法，包括以下步骤：
步骤一：获取河流概况、河流水质、河流水动力和河流生物四个维度的数据，通过人机交互的方式输入控制模块，控制模块对输入数据经过预处理后，输入神经网络模型，计算出通过曝气修复河流的目标值；
步骤二：启动曝气模块，感知模块和控制模块同步开始工作，根据曝气模块启动前输入控制模块的数据，以及感知模块实时监测的数据，每隔六小时计算一次河流最佳曝气量，按公式(1)(2)计算曝气量，向曝气模块发送控制信号，控制曝气的强度和时间；根据感知模块传输给控制模块的数据，实时监控河流曝气情况，维持河流的溶解氧浓度在预期水平；
Figure BDA0002991113060000081
e(k)＝DOmeasure(k)-DOset/default(k) (2)
式中，Q(k)表示曝气量，Kp表示e(k)项的系数，Ki表示
Figure BDA0002991113060000082
项的系数，Kd表示[e(k)-e(k-1)]项的系数，k表示k时刻，e(k)表示k时刻测量值与设置/预设溶解氧浓度值的偏差，e(k-1)表示k-1时刻测量值与设置/预设值的偏差，DOmeasure(k)表示测量得到的k时刻溶解氧浓度，DOset/default(k)表示k时刻曝气器控制模块溶解氧的输入值或预设值；
步骤三：循环步骤二，直至感知模块监测的河流的现状达到神经网络模型得到的曝气修复河流的目标值，停止曝气并通知管理人员。
其中，建立神经网络模型及其训练过程如下：
步骤1：收集基础数据，选取有代表性的河流，收集包括河流曝气修复前后常规理化指标、生物指标，以及河流宽度、深度、功能、周围土地利用类型、污染源情况，至少80组数据；
步骤2：数据预处理(该数据预处理过程与下面模型使用过程中的数据预处理是一致的。建立神经网络模型时，首先对数据预处理，使用神经网络模型测试数据时，也要首先进行同样的数据预处理，剔除明显异常数据。)，首先剔除明显异常和错误的数据，然后基于公式(3)做共线性分析，剔除共线性高的数据，再用排序算法做进一步降维，将输入指标降低到10个以内；最后，将处理好的数据划分为训练样本集和测试样本集；(随机按7：3划分训练样本集和测试样本集)
Figure BDA0002991113060000083
式中，R2为决定系数；VIF≥20剔除，VIF＜20保留；
步骤3：训练神经网络模型，将预处理后的结果作为输入层，将溶解氧浓度、氨氮浓度、化学需氧量、河流耗氧速率、自养生物占比作为输出层，输入层神经元数量与预处理后的指标个数相同，输出层神经元数量为5个；隐含层数量与神经元个数根据神经网络模型的质量进行调整，默认为两层；完成建模后，训练神经网络模型，根据输出结果，评估模型的损失函数，按公式(4)计算平均平方误差MSE表示神经网络模型的损失，直至图谱表现为水平直线，判定为符合要求；
Figure BDA0002991113060000091
h(xi)代表预测的第i个值，yi代表实际的第i个值；m1为测试集样本数量；
步骤4：使用测试样本集测试模型的准确性，保证神经网络模型拟合良好，不存在欠拟合和过拟合的现象，准确率需达到90％以上(参见图2)。
上述神经网络模型建立以及训练后，可采集待曝气河流的基本数据，输入上述建模以及训练后的模型中，进行如下步骤，即运行神经网络模型：
(1)基于河流基础数据，计算河流复氧速率、弯曲度和连通度，
(1.1)复氧速率计算方法：对于水深较深、流速缓慢的河流，使用公式(5)(7)计算；对于水深较浅、流速较大(靠实验和经验确定)的河流，使用公式(6)(7)计算，
Figure BDA0002991113060000092
Figure BDA0002991113060000093
式中，k2为复氧系数；D为氧在水体中的分子扩散系数；u为平均流速；H为河流水深；J为水面坡降；κ为卡门常数，一般取0.4；g为重力加速度；
Figure BDA0002991113060000094
式中，O为溶解氧浓度；Os为饱和溶解氧浓度，t为时间；
使用公式(8)计算河流弯曲度：
Ie＝S/L (8)
式中，Ie为河流弯曲度；S为河流中泓线长度；L为河流起止断面距离；
(1.2)河流连通度计算方法：水系的环通度α使用公式(9)计算；水系节点率β使用公式(10)计算；水系网络连通度γ使用公式(11)计算，
Figure BDA0002991113060000095
Figure BDA0002991113060000096
Figure BDA0002991113060000097
式中，α为水系的环通度；M为城市连接的廊道数量；V为水系节点个数；β为水系节点率；Mmax为最大的廊道数量；
(2)层次分析法结合熵权法，对水质指标和水动力指标降维；
(2.1)层次分析法得到主观赋权的权重值：
2.1.a.建立层次结构模型，根据收集的水质和水动力参数，分目标层、准则层和方案层建立层次结构模型；
2.1.b.建立判断矩阵，根据步骤a建立的层次结构模型构建成对比较矩阵(判断矩阵)，包括分别基于目标层、准则层和方案层的判断矩阵，并计算特征向量和特征值；
2.1.c.进行层次单排序及其一次性检验，分别计算准则层对目标层和方案层对准则层的权重排序，分别做一次性检验；
2.1.d.进行层次总排序及其一次性检验，计算方案层对目标的权重排序及其一致性检验；
(2.2)熵权法得到客观赋权的权重值：
2.2.a.构建原始数据矩阵X＝(xij)m×n，并归一化，
2.2.b.使用公式(12)计算熵值
Figure BDA0002991113060000101
式中，ej为j指标的熵值；m2为曝气河流样本数量；n为指标数量；Pij为xij归一化后的矩阵；
2.2.c.使用公式(13)计算权重
Figure BDA0002991113060000102
式中，Uj为j指标的权重；
2.2.d.使用公式(14)归一化处理
Figure BDA0002991113060000103
式中，βj为j指标的客观权重值；
利用公式(15)计算综合权重值
Figure BDA0002991113060000104
式中，w为由乘法归一组合赋权法得到的权重值；αj为主观赋权得到的权重值；βj为客观赋权得到的权重值；
(3)将生物数据划分为自养生物和异养生物，分别计算生物完整性指数；计算方法为：
3.a.选择受损点和参照点，选择生态修复前的数据为受损点，生态修复后的数据为参照点；
3.b.依据敏感性分析结果剔除参数，对生物参数做敏感性分析，其中IQ小于2的参数，不再做进一步分析，对于IQ值大于或等于2的参数，保留，并进入步骤3.c；
3.c.对剔除共线性参数，敏感性分析筛选后的参数进一步做Pearson相关分析和共线性分析，按照|r|>0.75或VIF≥20的标准，剔除共线性的参数；
Figure BDA0002991113060000111
Figure BDA0002991113060000112
式中，r为相关系数；VIF为方差膨胀因子；R2为决定系数；m2为曝气河流样本数量；xi,yi为成对的样本；
3.d.按照三分法，首先计算构成生物完整性指数值(IBI值)的参数在参照点的统计分布，拟定其评价指标的评分标准，据此将每个样点的参数的评价值实测值分别进行评分，然后将选出的参数得分进行累加，最后得到每个样点的IBI总分；
3.e.按照比值法，首先计算生物参数在全部样点中的95％或5％分位数的值；最后将选出的参数的分值进行累加，即得出各样点的IBI总分值。
(4)将上述计算结果输入训练好的神经网络模型，即可得到曝气修复河流的目标值。
步骤二所需的曝气量，即维持河流最佳需氧量的计算方法为：
步骤(1)：构建河流食物网模型，以生物与营养物质之间的相互作用、生物间相互作用为基础，构建河流食物网模型(参见图3)；
步骤(2)：基于动力学参数和实际收集的数据，模拟河流生态系统的生物量变化和营养物质变化，通过比对实际数据和模拟数据，率定动力学参数，使河流食物网模型的准确度符合需求；
率定动力学参数的过程是：先使用文献或实验测得的动力学参数，输入模型进行模拟，得到预测值与实际值比较，|预测值-实际值|/实际值可以表示模型的准确度。根据模型的需求不同，选取不同的准确度阈值。(例如要求准确度为80％，则根据准确度低于80％的预测值，修改该指标对应的动力学参数，直至所有指标的预测值达到80％。)
步骤(3)：对河流食物网模型的初始输入值做敏感性分析，找到控制河流生态系统的关键变量；
找到控制河流生态系统的关键变量的具体过程为：根据公式(18)计算不同变量的标准局部敏感系数
Figure BDA0002991113060000121
式中，Si,j为标准局部敏感系数；αi为模型参数相对变化量；Vj为模型输出变量的相对变化量；
Figure BDA0002991113060000122
为参数默认值改变1％；
标准局部敏感系数Si,j大于1表示变量敏感，即关键变量。
步骤(4)：以当前溶解氧浓度为基准，以当前温度下的饱和溶解氧浓度为上限，以0.2mg/L为间隔，生成溶解氧序列；将生成的溶解氧序列依次输入河流食物网模型，其余初始值使用当前时间的在线监测值，动力学参数使用经率定的参数；根据步骤(3)敏感性分析识别出的关键变量，以关键变量最优为目标，得到河流的最佳需氧量。
下面具体阐述河流最佳需氧量(曝气量)与曝气修复河流的目标值的关系：
参见图4，感知模块获取数据，传输给控制模块，控制模块将数据转化成控制信号，发送给曝气模块；控制信号有两种，一种是控制是否曝气，此时，需要曝气目标模型做计算，另一种是调节曝气的强弱，此时，需要最佳曝气量模型做计算。
曝气目标模型是基于神经网络算法的，包括训练模型和使用模型，训练模型仅训练一次即可，训练完成后，可一直使用；训练模型和使用模型使用的数据需要做相同的预处理(降维、去共线性、特征工程等)。模型训练过程如图5所示，其中，步骤8对模型准确度评估后，符合要求则结束；否则，重复步骤3，直至模型准确度符合要求。模型的使用过程如图6所示，在已训练模型的基础上，输入数据，预处理，经过神经网络模型运算，输出曝气修复河流的目标值。
最佳曝气量模型即上述步骤二的详细算法，可参见图7，通过输入数据，构建食物网模型，识别关键变量，计算最佳曝气量。
以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项发明技术思想的范围内，进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。

Claims (7)

1.一种用于河流生态修复的曝气系统的控制方法，其特征在于：包括以下步骤：
步骤一：获取河流概况、河流水质、河流水动力和河流生物四个维度的数据，通过人机交互的方式输入控制模块，控制模块对输入数据经过预处理后，输入神经网络模型，计算出通过曝气修复河流的目标值，
神经网络模型的建立方法为：
步骤1：收集四个维度的河流基础数据：选取有代表性的河流，收集包括河流曝气修复前后常规理化指标、生物指标，以及河流宽度、深度、功能、周围土地利用类型、污染源情况，至少80组数据；
步骤2：数据预处理：首先剔除明显异常和错误的数据，然后基于公式(3)做共线性分析，剔除共线性高的数据，再用排序算法做进一步降维，将输入指标降低到10个以内；最后，将处理好的数据划分为训练样本集和测试样本集；
Figure FDA0003722892030000011
R2为决定系数；VIF≥20剔除，VIF＜20保留；
步骤3：训练神经网络模型：将预处理得到的结果作为输入层，将溶解氧浓度、氨氮浓度、化学需氧量、河流耗氧速率、自养生物占比作为输出层，输入层神经元数量与预处理后的指标个数相同，输出层神经元数量为5个；隐含层数量与神经元个数根据神经网络模型的质量进行调整，默认为两层；完成建模后，训练神经网络模型，根据输出结果，评估模型的损失函数，按公式(4)计算平均平方误差MSE表示神经网络模型的损失，直至图谱表现为水平直线，判定为符合要求；
Figure FDA0003722892030000012
h(xi)代表预测的第i个值，yi代表实际的第i个值；m1为测试集样本数量；
步骤4：使用测试样本集测试模型的准确性，保证神经网络模型拟合良好，不存在欠拟合和过拟合的现象，准确率需达到90％以上；
步骤二：启动曝气模块，感知模块和控制模块同步开始工作，根据曝气模块启动前输入控制模块的数据，以及感知模块实时监测的数据，每隔六小时计算一次河流最佳曝气量，按公式(1)(2)计算曝气量，向曝气模块发送控制信号，控制曝气的强度和时间；根据感知模块传输给控制模块的数据，实时监控河流曝气情况，维持河流的溶解氧浓度在预期水平；
Figure FDA0003722892030000021
e(k)＝DOmeasure(k)-DOset/default(k) (2)
其中，Q(k)表示曝气量，Kp表示e(k)项的系数，Ki表示
Figure FDA0003722892030000022
项的系数，Kd表示[e(k)-e(k-1)]项的系数，k表示k时刻，e(k)表示k时刻测量值与设置/预设溶解氧浓度值的偏差，e(k-1)表示k-1时刻测量值与设置/预设值的偏差，DOmeasure(k)表示测量得到的k时刻溶解氧浓度，DOset/default(k)表示k时刻曝气器控制模块溶解氧的输入值或预设值；
步骤三：循环步骤二，直至感知模块监测的河流的现状达到神经网络模型得到的曝气修复河流的目标值，停止曝气并通知管理人员。
2.根据权利要求1所述的用于河流生态修复的曝气系统的控制方法，其特征在于：四个维度的河流基础数据处理过程如下：
(1)基于河流基础数据，计算河流复氧速率、弯曲度和连通度，
(1.1)复氧速率计算方法：对于水深较深、流速缓慢的河流，使用公式(5)(7)计算；对于水深较浅、流速较大的河流，使用公式(6)(7)计算，
Figure FDA0003722892030000023
Figure FDA0003722892030000024
式中，k2为复氧系数；D为氧在水体中的分子扩散系数；u为平均流速；H为河流水深；J为水面坡降；κ为卡门常数，取0.4；g为重力加速度；
Figure FDA0003722892030000025
式中，O为溶解氧浓度；Os为饱和溶解氧浓度，t为时间；
使用公式(8)计算河流弯曲度：
Ie＝S/L (8)
式中，Ie为河流弯曲度；S为河流中泓线长度；L为河流起止断面距离；
(1.2)河流连通度计算方法：水系的环通度α使用公式(9)计算；水系节点率β使用公式(10)计算；水系网络连通度γ使用公式(11)计算，
Figure FDA0003722892030000031
Figure FDA0003722892030000032
Figure FDA0003722892030000033
式中，α为水系的环通度；M为城市连接的廊道数量；V为水系节点个数；β为水系节点率；Mmax为最大的廊道数量；
(2)层次分析法结合熵权法，对水质指标和水动力指标降维；
(2.1)层次分析法得到主观赋权的权重值：
2.1.a.建立层次结构模型，根据收集的水质和水动力参数，分目标层、准则层和方案层建立层次结构模型；
2.1.b.建立判断矩阵，根据步骤a建立的层次结构模型构建成对比较矩阵，即判断矩阵，包括分别基于目标层、准则层和方案层的判断矩阵，并计算特征向量和特征值；
2.1.c.进行层次单排序及其一次性检验，分别计算准则层对目标层和方案层对准则层的权重排序，分别做一次性检验；
2.1.d.进行层次总排序及其一次性检验，计算方案层对目标的权重排序及其一致性检验；
(2.2)熵权法得到客观赋权的权重值：
2.2.a.构建原始数据矩阵X＝(xij)m×n，并归一化，
2.2.b.使用公式(12)计算熵值
Figure FDA0003722892030000034
其中i＝1，2，...，m；j＝1，2，...，n (12)
式中，ej为j指标的熵值；m2为曝气河流样本数量；n为指标数量；Pij为xij归一化后的矩阵；
2.2.c.使用公式(13)计算权重
Figure FDA0003722892030000035
其中j＝1，2，3，...，n (13)
式中，Uj为j指标的权重；
2.2.d.使用公式(14)归一化处理
Figure FDA0003722892030000036
式中，βj为j指标的客观权重值；
利用公式(15)计算综合权重值
Figure FDA0003722892030000041
式中，w为由乘法归一组合赋权法得到的权重值；αj为主观赋权得到的权重值；βj为客观赋权得到的权重值；
(3)将生物数据划分为自养生物和异养生物，分别计算生物完整性指数；计算方法为：
3.a.选择受损点和参照点，选择生态修复前的数据为受损点，生态修复后的数据为参照点；
3.b.依据敏感性分析结果剔除参数，对生物参数做敏感性分析，其中IQ小于2的参数，不再做进一步分析，对于IQ值大于或等于2的参数，保留，并进入步骤3.c；
3.c.对剔除共线性参数，敏感性分析筛选后的参数进一步做Pearson相关分析和共线性分析，
按照|r|＞0.75或VIF≥20的标准，剔除共线性的参数；
Figure FDA0003722892030000042
Figure FDA0003722892030000043
式中，r为相关系数；VIF为方差膨胀因子；R2为决定系数；m2为曝气河流样本数量；xi，yi为成对的样本；
3.d.按照三分法，首先计算构成生物完整性指数值(IBI值)的参数在参照点的统计分布，拟定其评价指标的评分标准，据此将每个样点的参数的评价值实测值分别进行评分，然后将选出的参数得分进行累加，最后得到每个样点的IBI总分；
3.e.按照比值法，首先计算生物参数在全部样点中的95％或5％分位数的值；最后将选出的参数的分值进行累加，即得出各样点的IBI总分值。
3.根据权利要求1所述的用于河流生态修复的曝气系统的控制方法，其特征在于：所述步骤二所需的曝气量，即维持河流最佳需氧量的计算方法为：
步骤(1)：构建河流食物网模型，以生物与营养物质之间的相互作用、生物间相互作用为基础，构建河流食物网模型；
步骤(2)：基于动力学参数和实际收集的数据，模拟河流生态系统的生物量变化和营养物质变化，通过比对实际数据和模拟数据，率定动力学参数，使河流食物网模型的准确度符合需求；
步骤(3)：对河流食物网模型的初始输入值做敏感性分析，找到控制河流生态系统的关键变量；
步骤(4)：以当前溶解氧浓度为基准，以当前温度下的饱和溶解氧浓度为上限，以0.2mg/L为间隔，生成溶解氧序列；将生成的溶解氧序列依次输入河流食物网模型，其余初始值使用当前时间的在线监测值，动力学参数使用经率定的参数；根据步骤(3)敏感性分析识别出的关键变量，以关键变量最优为目标，得到河流的最佳需氧量。
4.根据权利要求3所述的用于河流生态修复的曝气系统的控制方法，其特征在于：所述率定动力学参数的过程是：先使用文献或实验测得的动力学参数，输入模型进行模拟，得到预测值与实际值的比较，(预测值-实际值)/实际值表示模型的准确度。
5.根据权利要求4所述的用于河流生态修复的曝气系统的控制方法，其特征在于：找到控制河流生态系统的关键变量的具体过程为：根据公式(18)计算不同变量的标准局部敏感系数
Figure FDA0003722892030000051
式中，Si，j为标准局部敏感系数；αi为模型参数相对变化量；Vj为模型输出变量的相对变化量；
Figure FDA0003722892030000052
为参数默认值改变1％；
标准局部敏感系数Si，j大于1表示变量敏感，即关键变量。
6.权利要求1-5任一所述的用于河流生态修复的曝气系统的控制方法中的用于河流生态修复的曝气控制系统，其特征在于：包括感知模块、控制模块和曝气模块，所述感知模块用于采集河流概况、河流水质、河流水动力和河流生物四个维度的河流基础数据，所述控制模块用于分析处理感知模块采集到的数据，将其换算成控制曝气模块的参数，并给曝气模块发送控制信号。
7.根据权利要求6所述的用于河流生态修复的曝气控制系统，其特征在于：所述河流概况的数据包括河流宽度、河流深度、水面坡降、河流中弘线长度、河流起止断面距离、城市连接廊道数量、最大的廊道数量、水系节点个数、污染源情况、周围土地利用类型中的一种或多种；
所述河流水质的数据包括水温、溶解氧浓度、氧化还原电位、浊度、氨氮浓度、硝酸盐氮浓度、总氮、总磷、耗氧速率、生物需氧量、化学需氧量中的一种或多种；
所述河流水动力的数据包括河流的流速、流量、剪切力中的一种或多种；
所述河流生物的数据包括微生物、原生动物、后生动物、鱼类、藻类、浮游植物、沉水植物中的一种或多种。
CN202110315841.1A 2021-03-24 2021-03-24 一种用于河流生态修复的曝气控制系统及其控制方法 Active CN113011661B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110315841.1A CN113011661B (zh) 2021-03-24 2021-03-24 一种用于河流生态修复的曝气控制系统及其控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110315841.1A CN113011661B (zh) 2021-03-24 2021-03-24 一种用于河流生态修复的曝气控制系统及其控制方法

Publications (2)

Publication Number Publication Date
CN113011661A CN113011661A (zh) 2021-06-22
CN113011661B true CN113011661B (zh) 2022-08-23

Family

ID=76406585

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110315841.1A Active CN113011661B (zh) 2021-03-24 2021-03-24 一种用于河流生态修复的曝气控制系统及其控制方法

Country Status (1)

Country Link
CN (1) CN113011661B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113698029B (zh) * 2021-07-14 2022-07-22 河海大学 一种闸控非恒定水位河口水质净化系统及其运行方法
CN113830896B (zh) * 2021-09-27 2022-12-27 湖南工商大学 基于图神经网络的好氧段氧气预警方法、系统及设备
CN114169242A (zh) * 2021-12-10 2022-03-11 浙江嘉科新能源科技有限公司 一种基于参数不确定性分析河道生态增氧的智能控制算法
CN117541078B (zh) * 2023-11-21 2024-05-28 交通运输部规划研究院 一种基于人工运河开发的生态保护策略定制方法
CN117585794B (zh) * 2024-01-19 2024-03-26 四川绿境科兴环境科技有限公司 污水处理曝气控制方法、装置、设备及存储介质
CN117699999B (zh) * 2024-02-06 2024-04-26 深圳市深水龙岗水务集团有限公司 一种用于供水厂处理工艺的溶氧曝气监控系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101607762A (zh) * 2009-07-09 2009-12-23 东莞圣源环保科技有限公司 河流水体生态修复方法
CN104915805A (zh) * 2015-07-07 2015-09-16 郑州大学 一种河流生态修复效果评估方法
CN107480775A (zh) * 2017-08-14 2017-12-15 中国水产科学研究院淡水渔业研究中心 一种基于数据修复的池塘溶解氧预测方法
AU2020101936A4 (en) * 2018-05-04 2020-10-01 Jinan University Method of eutrophic water body restoration based on ecosystem model and biomanipulation technology

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101607762A (zh) * 2009-07-09 2009-12-23 东莞圣源环保科技有限公司 河流水体生态修复方法
CN104915805A (zh) * 2015-07-07 2015-09-16 郑州大学 一种河流生态修复效果评估方法
CN107480775A (zh) * 2017-08-14 2017-12-15 中国水产科学研究院淡水渔业研究中心 一种基于数据修复的池塘溶解氧预测方法
AU2020101936A4 (en) * 2018-05-04 2020-10-01 Jinan University Method of eutrophic water body restoration based on ecosystem model and biomanipulation technology

Also Published As

Publication number Publication date
CN113011661A (zh) 2021-06-22

Similar Documents

Publication Publication Date Title
CN113011661B (zh) 一种用于河流生态修复的曝气控制系统及其控制方法
CN108564229A (zh) 一种基于arima模型的工业污水进水量预测的方法
CN108898215B (zh) 一种基于二型模糊神经网络的污泥膨胀智能辨识方法
CN111291937A (zh) 基于支持向量分类与gru神经网络联合的处理污水水质预测方法
CN115577018B (zh) 一种水质监测数据的智能处理方法及系统
CN102313796B (zh) 一种污水处理生化需氧量软测量方法
CN110675036B (zh) 一种基于随机森林优化微生物指数的高原河流生态健康评价方法
CN103606006B (zh) 基于自组织t‑s模糊神经网络的污泥沉降指数软测量方法
CN107480775A (zh) 一种基于数据修复的池塘溶解氧预测方法
CN108562709A (zh) 一种基于卷积自编码器极限学习机的污水处理系统水质监测预警方法
CN111598754B (zh) 基于水质模型的河网区水质断面污染贡献率核算方法
CN112966891A (zh) 一种河流水环境质量预测方法
CN110132629B (zh) 一种利用支持向量机预测农村生活污水处理设施运行有效性的方法
CN105092493B (zh) 一种基于水样类型识别的水体cod光学测量方法
CN117170221A (zh) 一种污水处理的人工智能控制系统
CN104063609B (zh) 一种利用神经网络辅助判定污染源监测数据有效性的方法
CN109657287B (zh) 基于综合得分法水文模型精度识别的方法
CN105160422B (zh) 基于自组织级联神经网络的污水处理出水总磷预测方法
CN107665288A (zh) 一种化学需氧量的水质软测量预测方法
CN117196883A (zh) 一种基于人工智能的污水处理决策优化方法及系统
CN116011222A (zh) 一种基于主客观赋权云模型的隧道结构安全评价方法
CN112561307B (zh) 流域水环境大数据画像系统及方法
CN113269384B (zh) 一种对河流系统健康状态进行预警的方法
CN115422850A (zh) 基于efdc和边缘计算判断河流入河排污口污染贡献率的方法
CN112634347B (zh) 一种活性污泥形态学和污泥体积指数svi软测量方法

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