CN116861224B - 基于间歇过程软测量建模方法的间歇过程软测量建模系统 - Google Patents

基于间歇过程软测量建模方法的间歇过程软测量建模系统 Download PDF

Info

Publication number
CN116861224B
CN116861224B CN202311127076.6A CN202311127076A CN116861224B CN 116861224 B CN116861224 B CN 116861224B CN 202311127076 A CN202311127076 A CN 202311127076A CN 116861224 B CN116861224 B CN 116861224B
Authority
CN
China
Prior art keywords
data
matrix
representing
batch
intermittent
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
CN202311127076.6A
Other languages
English (en)
Other versions
CN116861224A (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.)
Beijing Paike Shenghong Electronic Technology Co ltd
Original Assignee
Ludong 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 Ludong University filed Critical Ludong University
Priority to CN202311127076.6A priority Critical patent/CN116861224B/zh
Publication of CN116861224A publication Critical patent/CN116861224A/zh
Application granted granted Critical
Publication of CN116861224B publication Critical patent/CN116861224B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Operations Research (AREA)
  • Computing Systems (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明属于生产过程数据处理及产品质量在线测量技术领域，公开了基于间歇过程软测量建模方法的间歇过程软测量建模系统。该系统搭载在青霉素间歇生产线，利用间歇过程软测量建模方法获取的在线质量预测值输入青霉素仿真平台，对青霉素生产工艺参数进行调控，获取青霉素最终生产工艺参数数据；该方法包括：对间歇过程的数据进行预处理；对间歇过程的数据进行特征提取；利用最小二乘支持向量回归方法对间歇过程进行预测；利用工业生产过程数据进行产品质量在线预测。本发明通过将间歇过程数据提取出来与最小二乘支持向量机有效结合的方法来建立模型，可以针对数据进行有效的预测，避免不良因素对生产质量的影响，从而提高生产效率。

Description

基于间歇过程软测量建模方法的间歇过程软测量建模系统
技术领域
本发明属于生产过程数据处理及产品质量在线测量技术领域，尤其涉及一种基于间歇过程软测量建模方法的间歇过程软测量建模系统。
背景技术
现在的工业生产过程分为连续生产过程和间歇生产过程两类。连续生产就是运行状态长期稳定运行，不改变生产条件。通常这种生产方式主要是生产一些生产量巨大的单一产品。间歇过程就是在一定的时间里，按照制定好的工艺流程进行工业生产，使得原材料被加工成为预想的产品。但是整个过程有许多因素具有不稳定性，一旦生产不稳定，生产出来的产品与想要的产品之间产生误差，不利于工业生产过程中要求度较高的产品的生产。间歇过程在工业生产中运用时间长久并且是传统的生产方式，具有生产成本低、工作流程简单等特点。因此应当建立该生产过程的软测量模型，模型会使对该过程的生产质量有一定的预测，这样就可以避免各种不良因素对生产质量产生不良影响，有利于提高生产质量。
通过上述分析，现有技术存在的问题及缺陷为：由于生产条件的变化或运行条件的不稳定性，过程数据往往呈现出不同的数据特征，如非线性、重复性和慢时变特性等问题。这些数据存在的问题，对拟合模型的效果表现并不理想，影响了建模的有效性，不能够实时准确预报生产产品质量指标，从而降低泛化性能，影响预测质量且影响了生产效率。
发明内容
为克服相关技术中存在的问题，本发明公开实施例提供了一种基于间歇过程软测量建模方法的间歇过程软测量建模系统，具体属于工业生产过程产品质量预测的软测量建模技术领域，具体涉及一种基于慢特征分析与最小二乘支持向量回归的间歇过程软测量建模方法。
所述技术方案如下：一种基于间歇过程软测量建模方法的间歇过程软测量建模系统，该系统搭载在青霉素间歇生产线，利用间歇过程软测量建模方法获取的在线质量预测值输入青霉素仿真平台，对青霉素生产工艺参数进行调控，获取青霉素最终生产工艺参数数据，所述青霉素生产工艺参数包括：间歇过程中的底物浓度、溶解氧浓度、生物量浓度、通气体积、二氧化碳浓度、通气速率、搅拌器功率、基材进料流量、发酵罐温度、基材进料温度；
所述间歇过程软测量建模系统包括：批次运行状态的离线模型建立模块，用于利用LSFA算法，进行间歇过程的数据动态-非线性特性的信息挖掘，建立批次运行状态的离线模型；
在线质量预测值计算模块，用于根据批次运行状态的离线模型，计算间歇过程在线数据的得分向量，并计算回归模型的在线质量预测值；
所述间歇过程软测量建模方法包括以下步骤：
S1，利用LSFA算法进行间歇过程的数据动态-非线性特性的信息挖掘，建立批次运行状态的离线模型；
S2，根据批次运行状态的离线模型，计算间歇过程在线数据的得分向量，并计算回归模型的在线质量预测值。
在步骤S1中，建立批次运行状态的离线模型，具体包括以下步骤：
S11：收集批次生产过程的数据，采用马氏距离对间歇过程的异常数据进行预处理；
S12：对间歇过程的数据进行动态矩阵扩展；
S13：对动态扩展矩阵进行主成分分析降维，提取间歇过程数据中的主元信息；
S14：对降维数据进行慢特征提取；
S15：利用最小二乘支持向量回归方法对间歇过程进行预测，构建质量与慢特征之间的回归模型。
在步骤S11中，采用马氏距离对间歇过程的异常数据进行预处理，具体过程如下：
收集多个批次生产过程的数据，第个批次过程输入矩阵为/>，其中，/>表示变量个数，/>表示样本采样次数，/>表示实数集，/>为第/>个批次过程输入矩阵；
计算同一批次过程输入数据的均值，其中，/>表示均值个数，表示同一批次过程输入数据的均值，/>表示第一批数据的样本均值，/>表示矩阵的转置；
变量与/>间的协方差矩阵为/>：
（1）
式中，表示协方差矩阵，/>表示样本个数，/>表示第/>行第/>列变量，/>表示第/>列均值，/>表示第/>行第/>列变量，/>表示第/>列均值；
分别计算第个批次过程输入矩阵/>中每一个样本点到均值的马氏距离：
（2）
式中，表示样本马氏距离，/>表示样本变量，/>表示均值，/>表示矩阵转置，/>表示协方差矩阵；
对马氏距离进行升序排序，通过设置阈值，遍历选取前/>个马氏距离最小的数据进行保留，剔除掉不符合设定阈值的异常值。
在步骤S12中，对间歇过程的数据进行动态矩阵扩展，具体包括：
批次测量数据的样本为，其中/>分别表示测量的批次数、工艺变量和采样次数，批次数据按三维数据形式构造；将数据按照时间序列分批展开，展开后的矩阵：
（3）
式中，表示扩展后/>个批次的数据矩阵，/>表示第/>个批次中第/>时刻采样的第/>个变量值；
引入时移重复作为附加变量来扩展原始输入矩阵，时刻的增广输入矩阵：
（4）
式中，表示增广矩阵，/>表示第/>个批次的样本矩阵；
通过上述方法对间歇过程的自变量进行动态扩展，得到扩展后结果/>。
在步骤S13中，对动态扩展矩阵进行主成分分析降维，提取间歇过程数据中的主元信息，具体包括：
对数据进行动态矩阵扩展后，获得的数据有个样本/>，其中，/>是一个/>维向量，对矩阵进行线性变换，降维后的矩阵为：
（5）
其中，是一个/>维向量并且/>；/>是一个/>的矩阵；
令；/>表示负载矩阵，/>表示负载向量，/>分别表示降维前后样本的维数，/>表示第/>维负载向量的转置，/>表示矩阵转置；
根据样本求均值：
（6）
样本的投影均值计算为：
（7）
根据样本均值计算样本的协方差矩阵：
（8）
计算样本投影后的方差：
（9）
根据拉格朗日函数构造目标函数为：
（10）
对目标函数求偏导：
（11）
令偏导等于0，得到结果：
（12）
式中，是/>的特征值，/>是对应的特征向量，只对/>进行排序，将前/>个最大特征值所对应的特征向量提取出来得到矩阵/>。
在步骤S14中，对降维数据进行慢特征提取，具体包括：
通过对间歇过程的自变量进行PCA降维，得到降维后结果/>，根据公式（7）、（9）计算降维后数据/>的均值/>和标准差/>，对降维后数据/>进行标准化得到：；/>表示标准化后的数据；
求数据的协方差矩阵，并进行SVD分解得到酉矩阵：
（13）
式中，表示标准化数据矩阵的转置，/>表示SVD分解后的酉矩阵，/>表示矩阵转置；
利用分解后的矩阵对进行白化变换：
（14）
（15）
式中，表示白化矩阵，/>表示奇异值矩阵，/>表示酉矩阵的转置，/>表示白化变量的一阶导数矩阵，/>表示连续两个样本之间的差；
求解白化矩阵的协方差矩阵，进行二次SVD分解：
（16）
式中，表示矩阵/>的一阶导数协方差矩阵，/>表示酉矩阵的转置，/>表示奇异值矩阵，/>表示酉矩阵；
按照升序进行排序，提取前列作为数据的特征向量/>，最终得到负载向量/>：
（17）
式中，表示负载矩阵；
该步骤保证提取的特征均值为0，方差为1且各向量之间不相关；
根据得到的权重向量，求解慢特征：
（18）
式中，表示慢特征。
在步骤S15中，构建质量与慢特征之间的回归模型，具体包括：
将获得的数据采用LSSVR进行回归预测，得到优化的目标函数为：
；
（19）
式中，为权向量，/>为偏置向量，/>为模型的误差变量，/>为惩罚系数，/>为目标函数，/>为负载向量的转置，/>为非线性映射函数；
根据拉格朗日函数：
（20）
其中，为拉格朗日乘子，/>为构造的拉格朗日函数；
根据最优化KKT条件求解，得到模型公式为：
（21）
其中，表示核函数。
在步骤S2中，根据批次运行状态的离线模型，计算间歇过程在线数据的得分向量，并计算回归模型的在线质量预测值，具体包括：
S20：采集新过程样本数据；
S21：计算中样本点到均值的马氏距离：
（22）
其中，均值采用离线建模阶段训练得到的同一批次过程输入数据的均值，根据求解的马氏距离值与阈值比较，若值大于阈值，则对测试数据的异常值进行剔除；
S22：构建时序增广矩阵，对新过程样本进行动态矩阵扩展：
（23）
S23：对动态扩展矩阵进行主成分分析降维，根据步骤S13中训练的权重向量，提取数据中的主元信息：
（24）
S24：步骤S14中计算出负载向量，将获得降维后的数据经过慢特征分析方法进行慢特征提取：
（25）
S25：利用最小二乘支持向量回归方法对间歇过程进行预测，构建质量与慢特征之间的回归模型，最终得到预测值：
（26）
其中，是拉格朗日乘子，/>为偏置向量，/>为核函数。
本发明的另一目的在于提供一种间歇过程软测量建模系统，实施所述的间歇过程软测量建模方法，该系统包括：
批次运行状态的离线模型建立模块，用于利用LSFA算法，进行间歇过程的数据动态-非线性特性的信息挖掘，建立批次运行状态的离线模型；
在线质量预测值计算模块，用于根据批次运行状态的离线模型，计算间歇过程在线数据的得分向量，并计算回归模型的在线质量预测值。
进一步，利用所述间歇过程软测量建模系统搭建精细化工间歇生产线，进行生产工艺参数调控。
结合上述的所有技术方案，本发明所具备的优点及积极效果为：本发明提供一种间歇过程软测量建模方法包括以下步骤：对间歇过程的数据进行预处理：对间歇过程的数据进行特征提取；利用最小二乘支持向量回归方法对间歇过程进行离线建模；利用工业生产过程数据进行在线预测。本发明中基于线性慢特征分析与最小二乘支持向量回归的建模方法，相较于传统的数据驱动建模方法，能够在批处理数据的情况下，有效提高新过程关键变量预测精度，避免不良因素对生产质量的影响，从而提高生产效率。通过对间歇过程中数据的提取，可以提高模型的生产质量预测精度。利用慢特征方法核函数跟踪过程的变化趋势，从而提高模型的预测精度。通过青霉素发酵过程仿真平台Pensim2.0模型产生数据，利用MATLAB仿真验证了所用方法的有效性。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本公开的实施例，并与说明书一起用于解释本公开的原理：
图1是本发明实施例提供的间歇过程软测量建模系统示意图；
图2是本发明实施例提供的间歇过程软测量建模方法流程图；
图3是本发明实施例提供的间歇过程软测量建模方法框架图；
图4是本发明实施例提供的间歇过程软测量建模方法原理图；
图5是本发明实施例提供的根据本发明所述的实施例的训练集质量预测结果图；
图6是本发明实施例提供的根据本发明所述的实施例的测试集质量预测结果图；
图中：1、批次运行状态的离线模型建立模块；2、在线质量预测值计算模块。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其他方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似改进，因此本发明不受下面公开的具体实施的限制。
实施例1，本发明实施例提供一种基于间歇过程软测量建模方法的间歇过程软测量建模系统，该系统搭载在青霉素间歇生产线，利用间歇过程软测量建模方法获取的在线质量预测值输入青霉素仿真平台，对青霉素生产工艺参数进行调控，获取青霉素最终生产工艺参数数据，所述青霉素生产工艺参数包括：间歇过程中的底物浓度、溶解氧浓度、生物量浓度、通气体积、二氧化碳浓度、通气速率、搅拌器功率、基材进料流量、发酵罐温度、基材进料温度。
如图1所示，本发明实施例提供的间歇过程软测量建模系统包括以下步骤：
批次运行状态的离线模型建立模块1，用于利用LSFA算法，进行间歇过程的数据动态-非线性特性的信息挖掘，建立批次运行状态的离线模型；
在线质量预测值计算模块2，用于根据批次运行状态的离线模型，计算间歇过程在线数据的得分向量，并计算回归模型的在线质量预测值。
如图2所示，本发明实施例提供的间歇过程软测量建模方法包括以下步骤：
S1，利用LSFA算法，进行间歇过程的数据动态-非线性特性的信息挖掘，建立批次运行状态的离线模型；
S2，根据批次运行状态的离线模型，计算间歇过程在线数据的得分向量，并计算回归模型的在线质量预测值。
在本发明实施例中，步骤S1利用LSFA算法，进行间歇过程的数据动态-非线性特性的信息挖掘，建立批次运行状态的离线模型具体包括以下步骤：
S11：收集批次生产过程的数据，采用马氏距离对间歇过程的异常数据进行预处理；
S12：对间歇过程的数据进行动态矩阵扩展；
S13：对动态扩展矩阵进行主成分分析降维，提取间歇过程的数据中的主元信息；
S14：对降维数据进行慢特征提取；
S15：利用最小二乘支持向量回归方法对间歇过程进行预测，构建质量与慢特征之间的回归模型。
数据预处理旨在采用马氏距离剔除异常值，步骤S11中，具体过程如下：
收集多个批次生产过程的数据，第个批次过程输入矩阵为/>，其中，/>表示变量个数，/>表示样本采样次数，/>表示实数集，/>为第/>个批次过程输入矩阵；
计算同一批次过程输入数据的均值，其中，/>表示均值个数，表示同一批次过程输入数据的均值，/>表示第一批数据的样本均值，/>表示矩阵的转置；
变量与/>间的协方差矩阵为/>：
（1）
式中，表示协方差矩阵，/>表示样本个数，/>表示第/>行第/>列变量，/>表示第/>列均值，/>表示第/>行第/>列变量，/>表示第/>列均值；
分别计算第个批次过程输入矩阵/>中每一个样本点到均值的马氏距离：
（2）
式中，表示样本马氏距离，/>表示样本变量，/>表示均值，/>表示矩阵转置，/>表示协方差矩阵；
对马氏距离进行升序排序，通过设置阈值，遍历选取前/>个马氏距离最小的数据进行保留，剔除掉不符合设定阈值的异常值。
动态矩阵扩展旨在解决系统中批与批之间的动态特性，步骤S12中，具体过程如下：
批次测量数据的样本为，其中/>分别表示测量的批次数、工艺变量和采样次数，批次数据按三维数据形式构造；将数据按照时间序列分批展开，扩展后的矩阵：
（3）
式中，表示扩展后/>个批次的数据矩阵，/>表示第/>个批次中第/>时刻采样的第/>个变量值；
引入时移重复作为附加变量来扩展原始输入矩阵，时刻的增强输入矩阵：
（4）
式中，表示增广矩阵，/>表示第/>个批次的样本矩阵；
通过上述方法对间歇过程的自变量进行动态扩展，得到扩展后结果/>。
主成分分析旨在对数据进行降维，将高维空间中无法进行比较的原始数据映射到低维空间，提取数据主元信息，步骤S13中，具体过程如下：
对数据进行动态矩阵扩展后，获得的数据有个样本/>，其中，/>是一个/>维向量，对矩阵进行线性变换，降维后的矩阵为：
（5）
其中，是一个/>维向量并且/>；/>是一个/>的矩阵；
令；/>表示负载矩阵，/>表示负载向量，/>分别表示降维前后样本的维数，/>表示第/>维负载向量的转置，/>表示矩阵转置；
根据样本求均值：
（6）
样本的投影均值计算为：
（7）
根据样本均值计算样本的协方差矩阵：
（8）
计算样本投影后的方差：
（9）
根据拉格朗日函数构造目标函数为：
（10）
对目标函数求偏导：
（11）
令偏导等于0，得到结果：
（12）
式中，是/>的特征值，/>是对应的特征向量，只对/>进行排序，将前/>个最大特征值所对应的特征向量提取出来得到矩阵/>。
慢特征分析旨在从输入数据中提取出变化尽可能缓慢的潜变量，步骤S14中，具体过程如下：
通过对间歇过程的自变量进行PCA降维，得到降维后结果/>，根据公式（7）、（9）计算降维后数据/>的均值/>和标准差/>，对降维后数据/>进行标准化得到：；/>表示标准化后的数据；
求数据的协方差矩阵，并进行SVD分解得到酉矩阵：
（13）
式中，表示标准化数据矩阵的转置，/>表示SVD分解后的酉矩阵，/>表示矩阵转置；利用分解后的矩阵对/>进行白化变换：
（14）
（15）
式中，表示白化矩阵，/>表示奇异值矩阵，/>表示酉矩阵的转置，/>表示白化变量的一阶导数矩阵，/>表示连续两个样本之间的差；
求解白化矩阵的协方差矩阵，进行二次SVD分解：
（16）
式中，表示矩阵/>的一阶导数协方差矩阵，/>表示酉矩阵的转置，/>表示奇异值矩阵，/>表示酉矩阵；
按照升序进行排序，提取前列作为数据的特征向量/>，最终得到负载向量/>：
（17）
式中，表示负载矩阵；
该步骤保证提取的特征均值为0，方差为1且各向量之间不相关；
根据得到的权重向量，求解慢特征：
（18）
式中，表示慢特征。
最小二乘支持向量回归旨在构建质量变量与慢特征变量之间的模型，对变量进行回归预测，步骤S15中，具体过程如下：
将获得的数据采用LSSVR进行回归预测，得到优化的目标函数为：
，
（19）
式中，为权向量，/>为偏置向量，/>为模型的误差变量，/>为惩罚系数，/>为目标函数，/>为负载向量的转置，/>为非线性映射函数；
根据拉格朗日函数：
（20）
其中，为拉格朗日乘子，/>为构造的拉格朗日函数；
根据最优化KKT条件求解，得到模型公式为：
（21）
其中，表示核函数。
采用试凑法对LSSVR的参数和/>进行寻优：保持选取建模的数据个数、数据的值和高斯核函数中核参数/>的值相同的情况下，设置LSSVR中/>的值在0.01~10000之间变化，为了减少计算量，采用逐渐增大值区间的方法形成/>的遍历向量，同样将其应用于/>的值向量。
通过计算LSSVR模型的RMSE（RootMeanSquareError），可以得到均方根误差矩阵，其中，/>表示/>的指标个数，/>表示/>的指标个数，其中设置/>。
通过选择LSSVR模型均方根误差矩阵的最小值，输出与索引值对应的两个参数值，并通过迭代参数值的上下限来找到最优参数值。
在本发明实施例中，步骤S2，根据批次运行状态的离线模型，计算间歇过程在线数据的得分向量，并计算回归模型的在线质量预测值具体包括：
S20：采集新过程样本数据；
S21：计算中样本点到均值的马氏距离：/>
（22）
其中，均值采用离线建模阶段训练得到的同一批次过程输入数据的均值，根据求解的马氏距离值与阈值比较，若值大于阈值，则对测试数据的异常值进行剔除；
S22：构建时序增广矩阵，对新过程样本进行动态矩阵扩展：
（23）
S23：对动态扩展矩阵进行主成分分析降维，根据步骤S13中训练的权重向量，提取数据中的主元信息：
（24）
S24：步骤S14中计算出负载向量，将获得降维后的数据经过慢特征分析方法进行慢特征提取：
（25）
S25：利用最小二乘支持向量回归方法对间歇过程进行预测，构建质量与慢特征之间的回归模型，最终得到预测值：
（26）
其中，是拉格朗日乘子，/>为偏置向量，/>为核函数。
LSSVR中的参数和/>应用模型训练选取的最优参数值。
本发明实施例解决了间歇过程的软测量模型因数据不足，与传统的模型相比，本发明的优点在于本方法可以在大量批处理数据的情况下，通过将间歇过程数据提取出来与最小二乘支持向量机有效结合的方法来建立模型，可以针对数据进行有效的预测，避免不良因素对生产质量的影响，从而提高生产效率。
实施例2，作为本发明的另一种实施方式，以测量的批次数，工艺变量和采样次数分别为50，9和800为例，如图3所示，本发明实施例提供的一种间歇过程软测量建模方法，包括：输入样本X，主成分分析，提取K个主成分、慢特征分析，选取S维慢特征，最小二乘支持向量回归，获得预测值Y。
如图4所示，具体为一种基于慢特征分析与最小二乘支持向量回归的间歇过程软测量建模方法，包括：
第一步：收集多个批次生产过程的数据，第个批次过程输入矩阵为/>，其中，/>表示变量个数，/>表示样本采样次数，/>表示实数集。计算同一批次过程输入数据的均值/>，其中，/>表示均值个数。
变量与/>间的协方差矩阵为/>：/>
（1）
分别计算第个批次过程输入矩阵/>中每一个样本点到均值的马氏距离：
（2）
对马氏距离进行升序排序，通过设置阈值，遍历选取前/>个马氏距离最小的数据进行保留，剔除掉不符合设定阈值的异常值。
第二步：选取批次测量数据的样本为，其中/>分别表示测量的批次数、工艺变量和采样次数，批次数据按三维数据形式构造；将数据按照时间序列分批展开，扩展后的矩阵：
（3）
引入时移重复作为附加变量来扩展原始输入矩阵，时刻的增强输入矩阵：
（4）
通过上述方法对间歇过程的自变量进行动态扩展，得到扩展后结果/>。
第三步：对数据进行动态矩阵扩展后，获得的数据有个样本/>，其中，/>是一个/>维向量，对矩阵进行线性变换，降维后的矩阵为：
（5）
其中，是一个/>维向量并且/>；/>是一个/>的矩阵；
令；/>表示负载矩阵，/>表示负载向量，/>分别表示降维前后样本的维数，/>表示第/>维负载向量的转置，/>表示矩阵转置；
根据样本求均值：
（6）
样本的投影均值计算为：
（7）
根据样本均值计算样本的协方差矩阵：
（8）
计算样本投影后的方差：
（9）
根据拉格朗日函数构造目标函数为：
（10）
对目标函数求偏导：
（11）
令偏导等于0，得到结果：
（12）
式中，是/>的特征值，/>是对应的特征向量，只对/>进行排序，将前/>个最大特征值所对应的特征向量提取出来得到矩阵/>。
第四步：通过对间歇过程的自变量进行PCA降维，得到降维后结果/>，根据公式（7）、（9）计算降维后数据/>的均值/>和标准差/>，对降维后数据/>进行标准化得到：；
求数据的协方差矩阵，并进行SVD分解得到酉矩阵：
（13）
利用分解后的矩阵对进行白化变换：
（14）
（15）
求解白化矩阵的协方差矩阵，进行二次SVD分解：
（16）
按照升序进行排序，提取前列作为数据的特征向量/>，最终得到负载向量/>：
（17）
该步骤保证提取的特征均值为0，方差为1且各向量之间不相关；
根据得到的权重向量，求解慢特征：
（18）
第五步：将获得的数据采用LSSVR进行回归预测，得到优化的目标函数为：
，
（19）
式中，为权向量，/>为偏置向量，/>为模型的误差变量，/>为惩罚系数，/>为目标函数，/>为负载向量的转置，/>为非线性映射函数；
根据拉格朗日函数：
（20）
其中，为拉格朗日乘子，/>为构造的拉格朗日函数；
根据最优化KKT条件求解，得到模型公式为：
（21）
其中，表示核函数。
采用试凑法对LSSVR的参数和/>进行寻优：保持选取建模的数据个数、数据的值和高斯核函数中核参数/>的值相同的情况下，设置LSSVR中/>的值在0.01 ~ 10000之间变化，为了减少计算量，采用逐渐增大值区间的方法形成/>的遍历向量，同样将其应用于/>的值向量。
通过计算LSSVR模型的RMSE（Root Mean Square Error），可以得到均方根误差矩阵 ，其中，/>表示/>的指标个数，/>表示/>的指标个数，其中设置/>。
通过选择LSSVR模型均方根误差矩阵的最小值，输出与索引值对应的两个参数值，并通过迭代参数值的上下限来找到最优参数值。
第六步：通过控制青霉素仿真平台Pensim2.0模型中的参数范围来模仿工业生产，来产生所需数据。具体过程如下：
设置间歇过程中的底物浓度取值范围为12-18 g/L, 溶解氧浓度为1.16 mmol/L，生物量浓度取值范围为0.05-0.15 g/L，通气体积取值范围为95-105 L,二氧化碳浓度设置为0.5 mmol/L，通气速率的取值范围为6-10 L/hr，搅拌器功率取值范围为25-35 W，基材进料流量取值范围为0.035-0.045 L/hr，发酵罐温度取值范围为298-300 K，基材进料温度取值范围为296-298 K；
通过改变参数的取值来输入到仿真平台中，训练生成不同的数据，本发明选取了50批次的数据构成模型的样本，其中样本的维度为选取的9个变量。挑选前30批数据返回步骤S11作为输入进行模型训练，直至步骤S15。在模型训练完成以后，输入剩余批次数据作为测试集返回步骤S20，直至步骤S25结束。
在本实施例中，本发明利用上述方法，通过上述参数设置，可得到青霉素仿真过程产品质量预测结果。为了验证本发明方法的有效性，本实例对比了该方法与经过动态内主成分分析方法处理数据的对比结果，以此来检验本发明的预测效果。
由图5和图6结果对比可以看出，本发明方法的回归预测结果更接近于真实值，相比另一种方法有所改进，提升了模型的预测精度，拟合效果更好。
需要注意的是，除非另有说明，本申请使用的技术术语或者科学术语应当为本发明所属领域技术人员所理解的通常意义。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。
以上所述，仅为本发明较优的具体的实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，都应涵盖在本发明的保护范围之内。

Claims (4)

1.一种基于间歇过程软测量建模方法的间歇过程软测量建模系统，其特征在于，该系统搭载在青霉素间歇生产线，利用间歇过程软测量建模方法获取的在线质量预测值输入青霉素仿真平台，对青霉素生产工艺参数进行调控，获取青霉素最终生产工艺参数数据，所述青霉素生产工艺参数包括：间歇过程中的底物浓度、溶解氧浓度、生物量浓度、通气体积、二氧化碳浓度、通气速率、搅拌器功率、基材进料流量、发酵罐温度、基材进料温度；
所述间歇过程软测量建模系统包括：批次运行状态的离线模型建立模块（1），用于利用LSFA算法，进行间歇过程的数据动态-非线性特性的信息挖掘，建立批次运行状态的离线模型；
在线质量预测值计算模块（2），用于根据批次运行状态的离线模型，计算间歇过程在线数据的得分向量，并计算回归模型的在线质量预测值；
所述间歇过程软测量建模方法包括以下步骤：
S1，利用LSFA算法进行间歇过程的数据动态-非线性特性的信息挖掘，建立批次运行状态的离线模型；
S2，根据批次运行状态的离线模型，计算间歇过程在线数据的得分向量，并计算回归模型的在线质量预测值；
在步骤S1中，建立批次运行状态的离线模型，具体包括以下步骤：
S11：收集批次生产过程的数据，采用马氏距离对间歇过程的异常数据进行预处理；
S12：对间歇过程的数据进行动态矩阵扩展；
S13：对动态扩展矩阵进行主成分分析降维，提取间歇过程数据中的主元信息；
S14：对降维数据进行慢特征提取；
S15：利用最小二乘支持向量回归方法对间歇过程进行预测，构建质量与慢特征之间的回归模型；
在步骤S13中，对动态扩展矩阵进行主成分分析降维，提取间歇过程数据中的主元信息，具体包括：
对数据进行动态矩阵扩展后，获得的数据有个样本/>，其中，/>是一个/>维向量，对矩阵进行线性变换，降维后的矩阵为：
（5）
其中，是一个/>维向量并且/>；/>是一个/>的矩阵；
令；/>表示负载矩阵，/>表示负载向量，/>分别表示降维前后样本的维数，/>表示第/>维负载向量的转置，/>表示矩阵转置；
根据样本求均值：
（6）
样本的投影均值计算为：
（7）
根据样本均值计算样本的协方差矩阵：
（8）
计算样本投影后的方差：
（9）
根据拉格朗日函数构造目标函数为：
（10）
对目标函数求偏导：
（11）
令偏导等于0，得到结果：
（12）
式中，是/>的特征值，/>是对应的特征向量，只对/>进行排序，将前/>个最大特征值所对应的特征向量提取出来得到矩阵/>；
在步骤S14中，对降维数据进行慢特征提取，具体包括：
通过对间歇过程的自变量进行PCA降维，得到降维后结果/>，根据公式（7）、（9）计算降维后数据/>的均值/>和标准差/>，对降维后数据/>进行标准化得到：/>；/>表示标准化后的数据；
求数据的协方差矩阵，并进行SVD分解得到酉矩阵：
（13）
式中，表示标准化数据矩阵的转置，/>表示SVD分解后的酉矩阵，/>表示矩阵转置；
利用分解后的矩阵对进行白化变换：
（14）
（15）
式中，表示白化矩阵，/>表示奇异值矩阵，/>表示酉矩阵的转置，/>表示白化变量的一阶导数矩阵，/>表示连续两个样本之间的差；
求解白化矩阵的协方差矩阵，进行二次SVD分解：
（16）
式中，表示矩阵/>的一阶导数协方差矩阵，/>表示酉矩阵的转置，/>表示奇异值矩阵，/>表示酉矩阵；
按照升序进行排序，提取前列作为数据的特征向量/>，最终得到负载向量/>：
（17）
式中，表示负载矩阵；
该步骤保证提取的特征均值为0，方差为1且各向量之间不相关；
根据得到的权重向量，求解慢特征：
（18）
式中，表示慢特征；
在步骤S15中，构建质量与慢特征之间的回归模型，具体包括：
将获得的数据采用LSSVR进行回归预测，得到优化的目标函数为：
；
（19）
式中，为负载向量，/>为偏置向量，/>为模型的误差变量，/>为惩罚系数，/>为目标函数，/>为负载向量的转置，/>为非线性映射函数；
根据拉格朗日函数：
（20）
其中，为拉格朗日乘子，/>为构造的拉格朗日函数；
根据最优化KKT条件求解，得到模型公式为：
（21）
其中，表示核函数。
2.根据权利要求1所述的基于间歇过程软测量建模方法的间歇过程软测量建模系统，其特征在于，在步骤S11中，采用马氏距离对间歇过程的异常数据进行预处理，具体过程如下：
收集多个批次生产过程的数据，第个批次过程输入矩阵为/>，其中，/>表示变量个数，/>表示样本采样次数，/>表示实数集，/>为第/>个批次过程输入矩阵；
计算同一批次过程输入数据的均值，其中，/>表示均值个数，/>表示同一批次过程输入数据的均值，/>表示第/>批数据的样本均值，/>表示矩阵的转置；
变量与/>间的协方差矩阵为/>：
（1）
式中，表示协方差矩阵，/>表示样本个数，/>表示第/>行第/>列变量，/>表示第/>列均值，/>表示第/>行第/>列变量，/>表示第/>列均值；
分别计算第个批次过程输入矩阵/>中每一个样本点到均值的马氏距离：
（2）
式中，表示样本马氏距离，/>表示样本变量，/>表示均值，/>表示矩阵转置，/>表示协方差矩阵；
对马氏距离进行升序排序，通过设置阈值，遍历选取前/>个马氏距离最小的数据进行保留，剔除掉不符合设定阈值的异常值。
3.根据权利要求1所述的基于间歇过程软测量建模方法的间歇过程软测量建模系统，其特征在于，在步骤S12中，对间歇过程的数据进行动态矩阵扩展，具体包括：
批次测量数据的样本为，其中/>分别表示测量的批次数、工艺变量和采样次数，批次数据按三维数据形式构造；将数据按照时间序列分批展开，展开后的矩阵：
（3）
式中，表示扩展后/>个批次的数据矩阵，/>表示第/>个批次中第/>时刻采样的第/>个变量值；
引入时移重复作为附加变量来扩展原始输入矩阵，时刻的增广输入矩阵：
（4）
式中，表示增广矩阵，/>表示第/>个批次的样本矩阵；
通过上述方法对间歇过程的自变量进行动态扩展，得到扩展后结果/>。
4.根据权利要求1所述的基于间歇过程软测量建模方法的间歇过程软测量建模系统，其特征在于，在步骤S2中，根据批次运行状态的离线模型，计算间歇过程在线数据的得分向量，并计算回归模型的在线质量预测值，具体包括：
S20：采集新过程样本数据；
S21：计算中样本点到均值的马氏距离：
（22）
其中，均值采用离线建模阶段训练得到的同一批次过程输入数据的均值，根据求解的马氏距离值与阈值比较，若值大于阈值，则对测试数据的异常值进行剔除；
S22：构建时序增广矩阵，对新过程样本进行动态矩阵扩展：
（23）
S23：对动态扩展矩阵进行主成分分析降维，根据步骤S13中训练的权重向量，提取数据中的主元信息：
（24）
S24：步骤S14中计算出负载向量，将获得降维后的数据经过慢特征分析方法进行慢特征提取：
（25）
S25：利用最小二乘支持向量回归方法对间歇过程进行预测，构建质量与慢特征之间的回归模型，最终得到预测值：
（26）
其中，是拉格朗日乘子，/>为偏置向量，/>为核函数。
CN202311127076.6A 2023-09-04 2023-09-04 基于间歇过程软测量建模方法的间歇过程软测量建模系统 Active CN116861224B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311127076.6A CN116861224B (zh) 2023-09-04 2023-09-04 基于间歇过程软测量建模方法的间歇过程软测量建模系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311127076.6A CN116861224B (zh) 2023-09-04 2023-09-04 基于间歇过程软测量建模方法的间歇过程软测量建模系统

Publications (2)

Publication Number Publication Date
CN116861224A CN116861224A (zh) 2023-10-10
CN116861224B true CN116861224B (zh) 2023-12-01

Family

ID=88223769

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311127076.6A Active CN116861224B (zh) 2023-09-04 2023-09-04 基于间歇过程软测量建模方法的间歇过程软测量建模系统

Country Status (1)

Country Link
CN (1) CN116861224B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117708691B (zh) * 2024-02-06 2024-05-10 东北大学 间歇过程监控方法及存储介质、计算机设备

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008048752A1 (en) * 2006-10-17 2008-04-24 International Business Machines Corporation Method and apparatus for improved regression modeling
CN102340620A (zh) * 2011-10-25 2012-02-01 重庆大学 基于马氏距离的视频图像背景检测方法
CN105138972A (zh) * 2015-08-11 2015-12-09 北京天诚盛业科技有限公司 人脸认证方法和装置
CN109858190A (zh) * 2019-03-13 2019-06-07 江南大学 一种基于果蝇算法优化梯度提升回归树的青霉素发酵过程软测量建模方法
CN110210687A (zh) * 2019-06-13 2019-09-06 中南大学 一种基于局部加权慢特征回归的非线性动态生产过程产品质量预测方法
CN110942258A (zh) * 2019-12-10 2020-03-31 山东科技大学 一种性能驱动的工业过程异常监测方法
CN111144017A (zh) * 2019-12-30 2020-05-12 北京化工大学 一种基于ff-rvm的多时段间歇过程软测量建模方法
CN114707424A (zh) * 2022-04-27 2022-07-05 浙江科技学院 基于质量相关慢特征分析算法的化工过程软测量方法
WO2022251701A1 (en) * 2021-05-27 2022-12-01 Sail Sv Llc X-ray imaging system and method
WO2023024210A1 (zh) * 2021-08-23 2023-03-02 江苏大学 一种基于傅里叶域主成分分析的数据降维方法
CN115985411A (zh) * 2022-10-17 2023-04-18 浙江科技学院 基于高斯过程回归模型的高分子熔融指数软测量方法
CN116300774A (zh) * 2023-05-23 2023-06-23 蓝星智云(山东)智能科技有限公司 基于主元分析与核密度估计的间歇过程可视化监控方法
WO2023160922A1 (en) * 2022-02-25 2023-08-31 Siemens Aktiengesellschaft Automatically quantifying an uncertainty of a prediction provided by a trained regression model

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105373098B (zh) * 2015-11-26 2017-10-24 江南大学 一种基于变量时段分解的间歇过程无模型在线滚动优化方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008048752A1 (en) * 2006-10-17 2008-04-24 International Business Machines Corporation Method and apparatus for improved regression modeling
CN102340620A (zh) * 2011-10-25 2012-02-01 重庆大学 基于马氏距离的视频图像背景检测方法
CN105138972A (zh) * 2015-08-11 2015-12-09 北京天诚盛业科技有限公司 人脸认证方法和装置
CN109858190A (zh) * 2019-03-13 2019-06-07 江南大学 一种基于果蝇算法优化梯度提升回归树的青霉素发酵过程软测量建模方法
CN110210687A (zh) * 2019-06-13 2019-09-06 中南大学 一种基于局部加权慢特征回归的非线性动态生产过程产品质量预测方法
CN110942258A (zh) * 2019-12-10 2020-03-31 山东科技大学 一种性能驱动的工业过程异常监测方法
CN111144017A (zh) * 2019-12-30 2020-05-12 北京化工大学 一种基于ff-rvm的多时段间歇过程软测量建模方法
WO2022251701A1 (en) * 2021-05-27 2022-12-01 Sail Sv Llc X-ray imaging system and method
WO2023024210A1 (zh) * 2021-08-23 2023-03-02 江苏大学 一种基于傅里叶域主成分分析的数据降维方法
WO2023160922A1 (en) * 2022-02-25 2023-08-31 Siemens Aktiengesellschaft Automatically quantifying an uncertainty of a prediction provided by a trained regression model
CN114707424A (zh) * 2022-04-27 2022-07-05 浙江科技学院 基于质量相关慢特征分析算法的化工过程软测量方法
CN115985411A (zh) * 2022-10-17 2023-04-18 浙江科技学院 基于高斯过程回归模型的高分子熔融指数软测量方法
CN116300774A (zh) * 2023-05-23 2023-06-23 蓝星智云(山东)智能科技有限公司 基于主元分析与核密度估计的间歇过程可视化监控方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
1种基于遗失数据重构的软测量方法;李庆华;;计算机与应用化学(第05期);全文 *
Soft-measuring approach to on-line predict BOD based on PCA time-delay neural network;Weili Ran;《 Fifth World Congress on Intelligent Control and Automation (IEEE Cat. No.04EX788)》;全文 *
基于多向核主元分析的青霉素生产过程在线监测;刘世成;王海清;李平;;浙江大学学报(工学版)(第02期);全文 *
局部最小二乘支持向量机回归在线建模方法及其在间歇过程的应用;刘毅;王海清;李平;;化工学报(第11期);全文 *

Also Published As

Publication number Publication date
CN116861224A (zh) 2023-10-10

Similar Documents

Publication Publication Date Title
CN116861224B (zh) 基于间歇过程软测量建模方法的间歇过程软测量建模系统
CN108667069B (zh) 一种基于偏最小二乘法回归的短期风电功率预测方法
CN110794782A (zh) 基于jy-mkpls的间歇工业过程在线质量预测方法
CN110263488B (zh) 一种基于集成即时学习的工业混炼胶门尼粘度软测量方法
CN107403196B (zh) 基于谱聚类分析的即时学习建模对脱丁烷塔底丁烷浓度进行预测的方法
CN102601881B (zh) 一种橡胶硬度的在线质量监控和预测模型更新的方法
CN110046377B (zh) 一种基于异构相似度的选择性集成即时学习软测量建模方法
CN112199829A (zh) 基于密度峰值聚类和即时学习的多工况多阶段批次过程监测方法
CN113031553B (zh) 一种基于自动聚类结合偏最小二乘的间歇过程质量预测方法
CN109616161B (zh) 一种基于孪生支持向量回归机的发酵过程软测量方法
CN111352408A (zh) 一种基于证据k近邻的多工况流程工业过程故障检测方法
CN102621953B (zh) 一种橡胶硬度的在线自动质量监控和预测模型更新的方法
CN117334271B (zh) 一种基于指定属性生成分子的方法
CN100370453C (zh) 稀土串级萃取分离组分含量软测量方法
CN112434739A (zh) 一种基于多核学习的支持向量机的化工过程故障诊断方法
CN116776252A (zh) 一种改进Mallow's Cp变量选择的工业过程软测量方法和系统
CN112506050B (zh) 一种基于潜变量过程迁移模型的间歇过程集成优化方法
CN110879873B (zh) 一种基于哈密顿蒙特卡洛采样的vine copula相关性描述的软测量方法及系统
CN114320266A (zh) 一种基于支持向量机的致密油藏常规井产量预测方法
Chen et al. Dynamic modeling of biotechnical process based on online support vector machine
CN111797574A (zh) 聚合物分子量分布的集成高斯过程回归模型方法
CN111948155B (zh) 精细全光谱结合改进gs-svr的复杂水体硝酸盐定量分析方法
CN111145838A (zh) 基于多向数据模型的青霉素发酵过程迭代学习卡尔曼滤波方法
CN111178627A (zh) 一种基于spca的神经网络混合优化预测方法
Ke et al. Data-driven real-time control method for process equipment in flow shop towards product quality improvement

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
TR01 Transfer of patent right

Effective date of registration: 20240310

Address after: 103, 1st floor, building 2, yard 16, Changxing Road, Fengtai District, Beijing

Patentee after: BEIJING PAIKE SHENGHONG ELECTRONIC TECHNOLOGY Co.,Ltd.

Country or region after: China

Address before: 264000 No. 186 Hongqi Middle Road, Zhifu District, Shandong, Yantai

Patentee before: LUDONG University

Country or region before: China

TR01 Transfer of patent right