CN113539382A - 一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统 - Google Patents

一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统 Download PDF

Info

Publication number
CN113539382A
CN113539382A CN202110591042.7A CN202110591042A CN113539382A CN 113539382 A CN113539382 A CN 113539382A CN 202110591042 A CN202110591042 A CN 202110591042A CN 113539382 A CN113539382 A CN 113539382A
Authority
CN
China
Prior art keywords
data
esterification
machine learning
learning model
reaction temperature
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
CN202110591042.7A
Other languages
English (en)
Other versions
CN113539382B (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.)
Zhejiang Supcon Software Co ltd
Zhejiang Supcon Technology Co Ltd
Original Assignee
Zhejiang Supcon Software Co ltd
Zhejiang Supcon 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 Zhejiang Supcon Software Co ltd, Zhejiang Supcon Technology Co Ltd filed Critical Zhejiang Supcon Software Co ltd
Priority to CN202110591042.7A priority Critical patent/CN113539382B/zh
Publication of CN113539382A publication Critical patent/CN113539382A/zh
Application granted granted Critical
Publication of CN113539382B publication Critical patent/CN113539382B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10Analysis or design of chemical reactions, syntheses or processes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/21Design, administration or maintenance of databases
    • G06F16/215Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/044Recurrent networks, e.g. Hopfield networks
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/70Machine learning, data mining or chemometrics
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/90Programming languages; Computing architectures; Database systems; Data warehousing
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Software Systems (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Medical Informatics (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Molecular Biology (AREA)
  • Mathematical Physics (AREA)
  • Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
  • General Factory Administration (AREA)

Abstract

本发明提出了一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统,包括:选取影响酯化反应温度、酯化真空度两个关键工艺参数的设备,基于选取的设备在亚磷酸二甲酯生产线上的关联位号采集历史数据;对采集到的历史数据进行数据清洗;根据清洗后的历史数据划分工况,将不同工况下的历史数据保存到样本库中;建立时序性机器学习模型,根据样本库对时序性机器学习模型进行训练;通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测,根据预测结果对酯化反应温度、酯化真空度进行分析,生成预警提示。本发明能够快速利用历史数据构建具有学习能力的时序性机器学习模型,克服机理建模过程耗时较久、建模过程复杂等问题。

Description

一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统
技术领域
本发明属于化工生产控制领域,尤其涉及一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统。
背景技术
亚磷酸二甲酯是一种重要的有机化工合成原料及中间体,广泛应用在有机膦酸型缓蚀剂、合成塑料助剂、染料添加剂、助燃剂以及农药的制备工业中。亚磷酸二甲酯生产过程中,酯化反应温度与酯化真空度两个影响生产效率的关键工艺参数很容易因工况影响而发生波动,当波动发生时,现阶段通常通过机理建模或人工判断的方式分析定位异常原因,均存在耗时耗力、建模复杂的问题,影响了关键工艺参数的监控效果,进而降低了亚磷酸二甲酯的生产效率。
发明内容
为了解决现有技术中存在的缺点和不足,本发明提出了一种亚磷酸二甲酯关键工艺参数的预警定位方法,包括:
选取影响酯化反应温度、酯化真空度两个关键工艺参数的设备,基于选取的设备在亚磷酸二甲酯生产线上的关联位号采集历史数据;
对采集到的历史数据进行数据清洗;
根据清洗后的历史数据划分工况,将不同工况下的历史数据保存到样本库中;
建立时序性机器学习模型,根据样本库对时序性机器学习模型进行训练;
通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测,根据预测结果对酯化反应温度、酯化真空度进行分析,生成预警提示。
可选的,所述对采集到的历史数据进行数据清洗,包括:
判断历史数据中是否存在缺失点,若存在则对缺失点进行填充,所述填充方法包括均值插补法、中位数插补法、众数插补法、k近邻法、回归法以及插值函数法;
确定历史数据中的异常数据,对异常数据进行平滑处理,所述平滑处理的方法包括滑动平均滤波法、一阶滞后滤波法、中位值平均滤波法以及加权递推平均滤波法。
可选的,所述根据清洗后的历史数据划分工况,将不同工况下的历史数据保存到样本库中,包括:
根据历史数据中的进料量、进料量配比、产量以及关键工艺参数划分不同的工况;
在每种工况下将历史数据按照预设比例保存到样本库中。
可选的,所述建立时序性机器学习模型,根据样本库对时序性机器学习模型进行训练,包括:
步骤一:基于LSTM的架构建立时序性机器学习模型;
步骤二:将样本库中的数据作为训练数据输入时序性机器学习模型,得到时序性机器学习模型的输出结果,记录从输入训练数据至得到输出结果的预测时长yTP
步骤三:将输出结果中的酯化反应温度预测值与同时序的实际酯化反应温度进行第一均方根误差的计算,将输出结果中的酯化真空度预测值与同时序的实际酯化真空度进行第二均方根误差的计算,将第一均方根误差和第二均方根误差相加得到预测精确度yACC
步骤四:基于评价公式EC对时序性机器学习模型的训练效果进行评价,所述评价公式为:
EC=w1yACC+w2yTP
其中,w1、w2分别为预设的精确度权重和预测时长权重,w1、w2的取值范围均为小于1的正数,yACC、yTP的取值范围均为正数;
步骤五:若EC的值大于预设训练阈值,则调整时序性机器学习模型中激活函数、神经网络层数、数据长度以及正则化的参数,重复步骤二至步骤四,直至EC的值不大于预设训练阈值时停止训练。
可选的,所述预警定位方法还包括,在通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测的同时,对样本库进行更新,更新过程包括:
当酯化反应温度、酯化真空度的预测结果与实际结果的误差超过预设值时,采集设备在此时的运行数据,将所述运行数据保存到样本库中。
可选的,所述通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测,根据预测结果对酯化反应温度、酯化真空度进行分析,生成预警提示,包括:
获取酯化反应温度、酯化真空度的预测结果,当预测结果超过预设预警值时获取当前关联位号对应的实时数据;
对实时数据进行相关性分析,根据实时数据与酯化反应温度、酯化真空度的相关程度进行排序,将排序结果对应的关联位号作为预警提示输出。
本发明还基于同样的思路提出了一种亚磷酸二甲酯关键工艺参数的预警定位系统,其特征在于,所述预警定位系统包括:
采集模块:用于选取影响酯化反应温度、酯化真空度两个关键工艺参数的设备,基于选取的设备在亚磷酸二甲酯生产线上的关联位号采集历史数据;
预处理模块:用于对采集到的历史数据进行数据清洗;
样本库模块:用于根据清洗后的历史数据划分工况,将不同工况下的历史数据保存到样本库中;
训练模块:用于建立时序性机器学习模型,根据样本库对时序性机器学习模型进行训练;
预警模块:用于通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测,根据预测结果对酯化反应温度、酯化真空度进行分析,生成预警提示。
可选的,所述样本库模块具体用于:
根据历史数据中的进料量、进料量配比、产量以及关键工艺参数划分不同的工况;
在每种工况下将历史数据按照预设比例保存到样本库中。
可选的,所述训练模块具体用于:
步骤一:基于LSTM的架构建立时序性机器学习模型;
步骤二:将样本库中的数据作为训练数据输入时序性机器学习模型,得到时序性机器学习模型的输出结果,记录从输入训练数据至得到输出结果的预测时长yTP
步骤三:将输出结果中的酯化反应温度预测值与同时序的实际酯化反应温度进行第一均方根误差的计算,将输出结果中的酯化真空度预测值与同时序的实际酯化真空度进行第二均方根误差的计算,将第一均方根误差和第二均方根误差相加得到预测精确度yACC
步骤四:基于评价公式EC对时序性机器学习模型的训练效果进行评价,所述评价公式为:
EC=w1yACC+w2yTP
其中,w1、w2分别为预设的精确度权重和预测时长权重,w1、w2的取值范围均为小于1的正数,yACC、yTP的取值范围均为正数;
步骤五:若EC的值大于预设训练阈值,则调整时序性机器学习模型中激活函数、神经网络层数、数据长度以及正则化的参数,重复步骤二至步骤四,直至EC的值不大于预设训练阈值时停止训练。
可选的,所述预警模块具体用于:
获取酯化反应温度、酯化真空度的预测结果,当预测结果超过预设预警值时获取当前关联位号对应的实时数据;
对实时数据进行相关性分析,根据实时数据与酯化反应温度、酯化真空度的相关程度进行排序,将排序结果对应的关联位号作为预警提示输出。
本发明提供的技术方案带来的有益效果是:
充分考虑工业企业亚磷酸二甲酯生产系统的实际运行状况,在不改变原有工艺和装置的前提下,能够快速利用历史数据构建具有学习能力的时序性机器学习模型,对影响亚磷酸二甲酯生产效果的酯化反应温度、酯化真空度进行预测,帮助工艺操作人员提前预判酯化反应温度、酯化真空度的发展趋势,克服机理建模过程耗时较久、建模过程复杂等问题,并通过关联位号实现生产线设备的快速定位,提前控制来稳定关键工艺参数,保证生产的安全性,延长生产周期,提高生产效率。
除此之外,本发明提出的技术方案结合了亚磷酸二甲酯生产过程中的数据类型完成数据清洗过程,提高了数据采集质量,进而提高预测精确度。
附图说明
为了更清楚地说明本发明的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提出的一种亚磷酸二甲酯关键工艺参数的预警定位方法的流程示意;
图2为本发明提出的一种亚磷酸二甲酯关键工艺参数的预警定位系统的结构框图。
具体实施方式
为使本发明的结构和优点更加清楚,下面将结合附图对本发明的结构作进一步地描述。
实施例一
如图1所示,本发明提出了一种亚磷酸二甲酯关键工艺参数的预警定位方法,包括:
S1:选取影响酯化反应温度、酯化真空度两个关键工艺参数的设备,基于选取的设备在亚磷酸二甲酯生产线上的关联位号采集历史数据;
S2:对采集到的历史数据进行数据清洗;
S3:根据清洗后的历史数据划分工况,将不同工况下的历史数据保存到样本库中;
S4:建立时序性机器学习模型,根据样本库对时序性机器学习模型进行训练;
S5:通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测,根据预测结果对酯化反应温度、酯化真空度进行分析,生成预警提示。
本实施例充分考虑工业企业亚磷酸二甲酯生产系统的实际运行状况,在不改变原有工艺和装置的前提下,能够快速利用历史数据构建具有学习能力的时序性机器学习模型,对影响亚磷酸二甲酯生产效果的酯化反应温度、酯化真空度进行预测,帮助工艺操作人员提前预判酯化反应温度、酯化真空度的发展趋势,克服机理建模过程耗时较久、建模过程复杂等问题,并通过关联位号实现生产线设备的快速定位,提前控制来稳定关键工艺参数,保证生产的安全性,延长生产周期,提高生产效率。
在本实施例中,首先需要梳理整个亚磷酸二甲酯生产系统的工艺路线,寻找可能影响酯化反应温度、酯化真空度的变量及其对应设备。同时,在采集历史数据之前可通过实际调研确定影响变量是否可测、仪表是否正常以及数据通信是否正常等问题,并根据采集数据中的缺失点、异常数据的分布比例评估数据质量,如果因仪表等故障导致数据质量差的问题,可通过校准、改造、更换仪表等方式解决,对于缺失关键测点等问题,可通过加装测点等方式保证获取关键信息数据。最终确认实际生产过程中设备对应的可测量、数据可靠的关键变量,例如酯化反应温度、三氯化磷与甲醇流量、冷冻水温度及各级反应釜真空度等。
为了提高采集到的历史数据的质量,本实施例对采集到的历史数据进行数据清洗的过程包括:
判断历史数据中是否存在缺失点,若存在则对缺失点进行填充,所述填充方法包括均值插补法、中位数插补法、众数插补法、k近邻法、回归法以及插值函数法。
所述均值插补法、中位数插补法、众数插补法即根据缺失点的属性,选择其他数据的平均值、中位数或众数进行插补。所述k近邻法为根据欧式距离或相关分析来确定与缺失点距离最近的K个数据,将这K个数据加权平均来估计该缺失点处的填补值。所述回归法为根据已有数据与其他相关变量的数据建立拟合模型来预测缺失点的填补值。所述插值函数法为利用已知数据建立合适的插值函数f(x),缺失点的填补值由对应点xi求出的函数值f(xi)近似代替。
确定历史数据中的异常数据,所述异常数据包括离群值以及流程行业噪声。对异常数据进行平滑处理,所述平滑处理的方法包括滑动平均滤波法、一阶滞后滤波法、中位值平均滤波法以及加权递推平均滤波法。
所述滑动平均滤波法包括:建立N个数据的滑窗,滑窗移动,每采进一个新数据,就释放最早采集的数据,而后求包括新数据在内的N个数据的算数平均值y,即有
Figure BDA0003089566620000081
其中,y为N个数据的平均值,N为滑窗长度,x为滑窗内的数据。滑动平均滤波平滑度高,但灵敏度低。
所述一阶滞后滤波法包括:
y(k)=αx(k)+(1-α)y(k-1);
其中y(k)平均滤波值,α为滤波系数,取值范围为[0,1],y(k-1)为上一次滤波输出值。一阶滞后滤波法对于周期性噪音具有良好的抑制作用,但数据滞后,灵敏度低。
所述中位值平均滤波法包括:建立N个数据的滑窗,滑窗移动,每采进一个新数据,就释放最早采集的数据,去除滑窗内的一个最大值和一个最小值,然后计算滑窗内的N-2个数据的算数平均值。能够有效的抑制随机干扰,但计算速度较慢。
所述加权递推平均滤波法包括:对滑动平均滤波法的改进,即不同时刻的数据加以不同的权通常是,越接近现时刻的数据,权取得越大。给予新采样值的权系数越大,则灵敏度越高,但信号平滑度低。
通过上述数据清洗的过程,进一步提高获取的历史数据的准确性,从而有利于提高后续时序性机器学习模型的训练效果以及预测准确性,并且能够及时对设备仪表出现的问题进行排查,尽最大限度保证数据采集的有效性,避免由于设备仪表问题反复出现数据采集异常的情况,从而提高数据清洗的效率。
在本实施例中,所述根据清洗后的历史数据划分工况,将不同工况下的历史数据保存到样本库中,包括:
根据历史数据中的进料量、进料量配比、产量以及关键工艺参数划分不同的工况。例如,根据关键工艺参数中的酯化反应温度划分为高温、正常反应温度、低温三种工况,根据产量划分为高产、低产、正常产量三种工况。
在每种工况下将历史数据按照预设比例保存到样本库中,用于后期时序性机器学习模型的建立。在一示例中对每种工况下的历史数据折算预设比例后保存到样本库中;在另一示例中对每种工况下的历史数据进行聚类,再将每个类的聚类中心保存到样本库中。
为了对样本库中的工况与数据能够及时适应实际生产情况,本实施例还包括在通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测的同时,对样本库进行更新,更新过程包括:
当酯化反应温度、酯化真空度的预测结果与实际结果的误差超过预设值时,采集设备在此时的运行数据,将所述运行数据保存到样本库中。
在本实施例中,所述建立时序性机器学习模型,根据样本库对时序性机器学习模型进行训练,包括:
步骤一:基于LSTM的架构建立时序性机器学习模型;
步骤二:将样本库中的数据作为训练数据输入时序性机器学习模型,得到时序性机器学习模型的输出结果,记录从输入训练数据至得到输出结果的预测时长yTP
步骤三:将输出结果中的酯化反应温度预测值与同时序的实际酯化反应温度进行第一均方根误差的计算,将输出结果中的酯化真空度预测值与同时序的实际酯化真空度进行第二均方根误差的计算,将第一均方根误差和第二均方根误差相加得到预测精确度yACC
步骤四:基于评价公式EC对时序性机器学习模型的训练效果进行评价,所述评价公式为:
EC=w1yACC+w2yTP
其中,w1、w2分别为预设的精确度权重和预测时长权重,w1、w2的取值范围均为小于1的正数,yACC、yTP的取值范围均为正数;
步骤五:若EC的值大于预设训练阈值,则调整时序性机器学习模型中激活函数、神经网络层数、数据长度以及正则化的参数,重复步骤二至步骤四,直至EC的值不大于预设训练阈值时停止训练。
LSTM是一种用于处理序列数据的神经网络,主要是为了解决长序列训练过程中的梯度消失和梯度爆炸问题,相比普通的RNN,LSTM能够在更长的序列中有更好的表现。根据时序性机器学习模型建立预测模型的优点在于不需要再每个变量对目标变量的时滞时间进行单独的计算,只需对数据长度进行调整。本领域技术人员应当知道如何根据LSTM架构建立所述时序性机器学习模型,此处不再赘述。
通过本实施例提出的评价公式对时序性机器学习模型的训练效果进行判断,根据亚磷酸二甲酯的生产工艺需要,训练好的时序性机器学习模型应该兼具准确度及较短的预测时长,保证现场技术人员能够有时间处理未来发生的异常情况。以上述评价公式为目标,实现在线模型参数寻优实现最优目标,完成在线学习和自我优化。同时通过不断的训练迭代调整模型中的参数,本实施例中可通过网格搜索、随机搜索、贝叶斯优化算法、粒子群寻优等方法结合实现参数的优化更新。在特定的工况下,在基础的LSTM架构模型和样本库数据的基础上,以上述评价公式为目标,实现在线模型的激活函数、神经网络层数、数据长度、正则化方法的超参寻优,寻优实现最优目标,完成在线学习和自我优化。
在本实施例中,所述通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测,根据预测结果对酯化反应温度、酯化真空度进行分析,生成预警提示,包括:
获取酯化反应温度、酯化真空度的预测结果,当预测结果超过预设预警值时获取当前关联位号对应的实时数据;
将获取的实时数据与预先存储的标准数据范围进行比对,筛选出超出标准数据范围的实时数据;
获取酯化反应温度、酯化真空度的预测结果,当预测结果超过预设预警值时获取当前关联位号对应的实时数据;
对实时数据进行相关性分析,根据实时数据与酯化反应温度、酯化真空度的相关程度进行排序,将排序结果对应的关联位号作为预警提示输出。
本实施例中的相关性分析包括:监测实时数据,当酯化反应温度、酯化真空度等关键工艺参数出现在预设时段内持续上升或下降,或者超过控制上下限的情况下,采集最新的实时数据;对实时数据进行预处理以及数据对齐,预处理方法参考上述S2的过程;基于spearman系数或最大互信息系数分别分析实时数据与酯化反应温度、酯化真空度的相关程度。
本实施例通过训练时序性机器学习模型实现对影响酯化反应温度、酯化真空度的因素进行趋势分析以及超限分析,并根据相关程度表现的影响因素重要性调整相应的设备,例如若时序性机器学习模型对关键工艺指标酯化反应温度、酯化真空度的预测趋势超过控制人员设定的阈值,会发出预警提示,控制人员可根据预警提示中提供的排序结果找到排名靠前的关联位号对应的设备,通过控制手段如三氯化磷进料流量、甲醇进料流量、循环冷却水量进行提前调节,避免关键因素未来发生波动,将其控稳在一定范围。当关键工艺参数实际发生波动时,控制人员也可根据异常原因排序排查原因,缩短排查时间,缩短关键因素异常波动的时间。
实施例二
如图2所示,本实施例提出了一种亚磷酸二甲酯关键工艺参数的预警定位系统6,包括:
采集模块61:用于选取影响酯化反应温度、酯化真空度两个关键工艺参数的设备,基于选取的设备在亚磷酸二甲酯生产线上的关联位号采集历史数据;
预处理模块62:用于对采集到的历史数据进行数据清洗;
样本库模块63:用于根据清洗后的历史数据划分工况,将不同工况下的历史数据保存到样本库中;
训练模块64:用于建立时序性机器学习模型,根据样本库对时序性机器学习模型进行训练;
预警模块65:用于通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测,根据预测结果对酯化反应温度、酯化真空度进行分析,生成预警提示。
本实施例充分考虑工业企业亚磷酸二甲酯生产系统的实际运行状况,在不改变原有工艺和装置的前提下,能够快速利用历史数据构建具有学习能力的时序性机器学习模型,对影响亚磷酸二甲酯生产效果的酯化反应温度、酯化真空度进行预测,帮助工艺操作人员提前预判酯化反应温度、酯化真空度的发展趋势,克服机理建模过程耗时较久、建模过程复杂等问题,并通过关联位号实现生产线设备的快速定位,提前控制来稳定关键工艺参数,保证生产的安全性,延长生产周期,提高生产效率。
在本实施例中,首先需要梳理整个亚磷酸二甲酯生产系统的工艺路线,寻找可能影响酯化反应温度、酯化真空度的变量及其对应设备。同时,在采集历史数据之前可通过实际调研确定影响变量是否可测、仪表是否正常以及数据通信是否正常等问题,并根据采集数据中的缺失点、异常数据的分布比例评估数据质量,如果因仪表等故障导致数据质量差的问题,可通过校准、改造、更换仪表等方式解决,对于缺失关键测点等问题,可通过加装测点等方式保证获取关键信息数据。最终确认实际生产过程中设备对应的可测量、数据可靠的关键变量,例如酯化反应温度、三氯化磷与甲醇流量、冷冻水温度及各级反应釜真空度等。
为了提高采集到的历史数据的质量,本实施例中预处理模块62具体用于:
判断历史数据中是否存在缺失点,若存在则对缺失点进行填充,所述填充方法包括均值插补法、中位数插补法、众数插补法、k近邻法、回归法以及插值函数法。
所述均值插补法、中位数插补法、众数插补法即根据缺失点的属性,选择其他数据的平均值、中位数或众数进行插补。所述k近邻法为根据欧式距离或相关分析来确定与缺失点距离最近的K个数据,将这K个数据加权平均来估计该缺失点处的填补值。所述回归法为根据已有数据与其他相关变量的数据建立拟合模型来预测缺失点的填补值。所述插值函数法为利用已知数据建立合适的插值函数f(x),缺失点的填补值由对应点xi求出的函数值f(xi)近似代替。
确定历史数据中的异常数据,所述异常数据包括离群值以及流程行业噪声。对异常数据进行平滑处理,所述平滑处理的方法包括滑动平均滤波法、一阶滞后滤波法、中位值平均滤波法以及加权递推平均滤波法。
所述滑动平均滤波法包括:建立N个数据的滑窗,滑窗移动,每采进一个新数据,就释放最早采集的数据,而后求包括新数据在内的N个数据的算数平均值y,即有
Figure BDA0003089566620000131
其中,y为N个数据的平均值,N为滑窗长度,x为滑窗内的数据。滑动平均滤波平滑度高,但灵敏度低。
所述一阶滞后滤波法包括:
y(k)=αx(k)+(1-α)y(k-1);
其中y(k)平均滤波值,α为滤波系数,取值范围为[0,1],y(k-1)为上一次滤波输出值。一阶滞后滤波法对于周期性噪音具有良好的抑制作用,但数据滞后,灵敏度低。
所述中位值平均滤波法包括:建立N个数据的滑窗,滑窗移动,每采进一个新数据,就释放最早采集的数据,去除滑窗内的一个最大值和一个最小值,然后计算滑窗内的N-2个数据的算数平均值。能够有效的抑制随机干扰,但计算速度较慢。
所述加权递推平均滤波法包括:对滑动平均滤波法的改进,即不同时刻的数据加以不同的权通常是,越接近现时刻的数据,权取得越大。给予新采样值的权系数越大,则灵敏度越高,但信号平滑度低。
通过上述数据清洗的过程,进一步提高获取的历史数据的准确性,从而有利于提高后续时序性机器学习模型的训练效果以及预测准确性,并且能够及时对设备仪表出现的问题进行排查,尽最大限度保证数据采集的有效性,避免由于设备仪表问题反复出现数据采集异常的情况,从而提高数据清洗的效率。
在本实施例中,所述样本库模块63具体用于:
根据历史数据中的进料量、进料量配比、产量以及关键工艺参数划分不同的工况。例如,根据关键工艺参数中的酯化反应温度划分为高温、正常反应温度、低温三种工况,根据产量划分为高产、低产、正常产量三种工况。
在每种工况下将历史数据按照预设比例保存到样本库中,用于后期时序性机器学习模型的建立。在一示例中对每种工况下的历史数据折算预设比例后保存到样本库中;在另一示例中对每种工况下的历史数据进行聚类,再将每个类的聚类中心保存到样本库中。
为了对样本库中的工况与数据能够及时适应实际生产情况,本实施例还包括更新模块,用于在通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测的同时,对样本库进行更新,更新过程包括:
当酯化反应温度、酯化真空度的预测结果与实际结果的误差超过预设值时,采集设备在此时的运行数据,将所述运行数据保存到样本库中。
在本实施例中,所述训练模块64具体用于:
步骤一:基于LSTM的架构建立时序性机器学习模型;
步骤二:将样本库中的数据作为训练数据输入时序性机器学习模型,得到时序性机器学习模型的输出结果,记录从输入训练数据至得到输出结果的预测时长yTP
步骤三:将输出结果中的酯化反应温度预测值与同时序的实际酯化反应温度进行第一均方根误差的计算,将输出结果中的酯化真空度预测值与同时序的实际酯化真空度进行第二均方根误差的计算,将第一均方根误差和第二均方根误差相加得到预测精确度yACC
步骤四:基于评价公式EC对时序性机器学习模型的训练效果进行评价,所述评价公式为:
EC=w1yACC+w2yTP
其中,w1、w2分别为预设的精确度权重和预测时长权重,w1、w2的取值范围均为小于1的正数,yACC、yTP的取值范围均为正数;
步骤五:若EC的值大于预设训练阈值,则调整时序性机器学习模型中激活函数、神经网络层数、数据长度以及正则化的参数,重复步骤二至步骤四,直至EC的值不大于预设训练阈值时停止训练。
LSTM是一种用于处理序列数据的神经网络,主要是为了解决长序列训练过程中的梯度消失和梯度爆炸问题,相比普通的RNN,LSTM能够在更长的序列中有更好的表现。根据时序性机器学习模型建立预测模型的优点在于不需要再每个变量对目标变量的时滞时间进行单独的计算,只需对数据长度进行调整。本领域技术人员应当知道如何根据LSTM架构建立所述时序性机器学习模型,此处不再赘述。
通过本实施例提出的评价公式对时序性机器学习模型的训练效果进行判断,根据亚磷酸二甲酯的生产工艺需要,训练好的时序性机器学习模型应该兼具准确度及较短的预测时长,保证现场技术人员能够有时间处理未来发生的异常情况。以上述评价公式为目标,实现在线模型参数寻优实现最优目标,完成在线学习和自我优化。同时通过不断的训练迭代调整模型中的参数,本实施例中可通过网格搜索、随机搜索、贝叶斯优化算法、粒子群寻优等方法结合实现参数的优化更新。在特定的工况下,在基础的LSTM架构模型和样本库数据的基础上,以上述评价公式为目标,实现在线模型的激活函数、神经网络层数、数据长度、正则化方法的超参寻优,寻优实现最优目标,完成在线学习和自我优化。
在本实施例中,所述预警模块65具体用于:
获取酯化反应温度、酯化真空度的预测结果,当预测结果超过预设预警值时获取当前关联位号对应的实时数据;
将获取的实时数据与预先存储的标准数据范围进行比对,筛选出超出标准数据范围的实时数据;
获取酯化反应温度、酯化真空度的预测结果,当预测结果超过预设预警值时获取当前关联位号对应的实时数据;
对实时数据进行相关性分析,根据实时数据与酯化反应温度、酯化真空度的相关程度进行排序,将排序结果对应的关联位号作为预警提示输出。
本实施例中的相关性分析包括:监测实时数据,当酯化反应温度、酯化真空度等关键工艺参数出现在预设时段内持续上升或下降,或者超过控制上下限的情况下,采集最新的实时数据;对实时数据进行预处理以及数据对齐,预处理方法参考上述预处理模块62执行的功能;基于spearman系数或最大互信息系数分别分析实时数据与酯化反应温度、酯化真空度的相关程度。
本实施例通过训练时序性机器学习模型实现对影响酯化反应温度、酯化真空度的因素进行趋势分析以及超限分析,并根据相关程度表现的影响因素重要性调整相应的设备,例如若时序性机器学习模型对关键工艺指标酯化反应温度、酯化真空度的预测趋势超过控制人员设定的阈值,会发出预警提示,控制人员可根据预警提示中提供的排序结果找到排名靠前的关联位号对应的设备,通过控制手段如三氯化磷进料流量、甲醇进料流量、循环冷却水量进行提前调节,避免关键因素未来发生波动,将其控稳在一定范围。当关键工艺参数实际发生波动时,控制人员也可根据异常原因排序排查原因,缩短排查时间,缩短关键因素异常波动的时间。
上述实施例中的各个序号仅仅为了描述,不代表各部件的组装或使用过程中的先后顺序。
以上所述仅为本发明的实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种亚磷酸二甲酯关键工艺参数的预警定位方法,其特征在于,所述预警定位方法包括:
选取影响酯化反应温度、酯化真空度两个关键工艺参数的设备,基于选取的设备在亚磷酸二甲酯生产线上的关联位号采集历史数据;
对采集到的历史数据进行数据清洗;
根据清洗后的历史数据划分工况,将不同工况下的历史数据保存到样本库中;
建立时序性机器学习模型,根据样本库对时序性机器学习模型进行训练;
通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测,根据预测结果对酯化反应温度、酯化真空度进行分析,生成预警提示。
2.根据权利要求1所述的一种亚磷酸二甲酯关键工艺参数的预警定位方法,其特征在于,所述对采集到的历史数据进行数据清洗,包括:
判断历史数据中是否存在缺失点,若存在则对缺失点进行填充,所述填充方法包括均值插补法、中位数插补法、众数插补法、k近邻法、回归法以及插值函数法;
确定历史数据中的异常数据,对异常数据进行平滑处理,所述平滑处理的方法包括滑动平均滤波法、一阶滞后滤波法、中位值平均滤波法以及加权递推平均滤波法。
3.根据权利要求1所述的一种亚磷酸二甲酯关键工艺参数的预警定位方法,其特征在于,所述根据清洗后的历史数据划分工况,将不同工况下的历史数据保存到样本库中,包括:
根据历史数据中的进料量、进料量配比、产量以及关键工艺参数划分不同的工况;
在每种工况下将历史数据按照预设比例保存到样本库中。
4.根据权利要求1所述的一种亚磷酸二甲酯关键工艺参数的预警定位方法,其特征在于,所述建立时序性机器学习模型,根据样本库对时序性机器学习模型进行训练,包括:
步骤一:基于LSTM的架构建立时序性机器学习模型;
步骤二:将样本库中的数据作为训练数据输入时序性机器学习模型,得到时序性机器学习模型的输出结果,记录从输入训练数据至得到输出结果的预测时长yTP
步骤三:将输出结果中的酯化反应温度预测值与同时序的实际酯化反应温度进行第一均方根误差的计算,将输出结果中的酯化真空度预测值与同时序的实际酯化真空度进行第二均方根误差的计算,将第一均方根误差和第二均方根误差相加得到预测精确度yACC
步骤四:基于评价公式EC对时序性机器学习模型的训练效果进行评价,所述评价公式为:
EC=w1yACC+w2yTP
其中,w1、w2分别为预设的精确度权重和预测时长权重,w1、w2的取值范围均为小于1的正数,yACC、yTP的取值范围均为正数;
步骤五:若EC的值大于预设训练阈值,则调整时序性机器学习模型中激活函数、神经网络层数、数据长度以及正则化的参数,重复步骤二至步骤四,直至EC的值不大于预设训练阈值时停止训练。
5.根据权利要求1所述的一种亚磷酸二甲酯关键工艺参数的预警定位方法,其特征在于,所述预警定位方法还包括在通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测的同时,对样本库进行更新,更新过程包括:
当酯化反应温度、酯化真空度的预测结果与实际结果的误差超过预设值时,采集设备在此时的运行数据,将所述运行数据保存到样本库中。
6.根据权利要求1所述的一种亚磷酸二甲酯关键工艺参数的预警定位方法,其特征在于,所述通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测,根据预测结果对酯化反应温度、酯化真空度进行分析,生成预警提示,包括:
获取酯化反应温度、酯化真空度的预测结果,当预测结果超过预设预警值时获取当前关联位号对应的实时数据;
对实时数据进行相关性分析,根据实时数据与酯化反应温度、酯化真空度的相关程度进行排序,将排序结果对应的关联位号作为预警提示输出。
7.一种亚磷酸二甲酯关键工艺参数的预警定位系统,其特征在于,所述预警定位系统包括:
采集模块:用于选取影响酯化反应温度、酯化真空度两个关键工艺参数的设备,基于选取的设备在亚磷酸二甲酯生产线上的关联位号采集历史数据;
预处理模块:用于对采集到的历史数据进行数据清洗;
样本库模块:用于根据清洗后的历史数据划分工况,将不同工况下的历史数据保存到样本库中;
训练模块:用于建立时序性机器学习模型,根据样本库对时序性机器学习模型进行训练;
预警模块:用于通过训练好的时序性机器学习模型对酯化反应温度、酯化真空度进行预测,根据预测结果对酯化反应温度、酯化真空度进行分析,生成预警提示。
8.根据权利要7所述的一种亚磷酸二甲酯关键工艺参数的预警定位系统,其特征在于,所述样本库模块具体用于:
根据历史数据中的进料量、进料量配比、产量以及关键工艺参数划分不同的工况;
在每种工况下将历史数据按照预设比例保存到样本库中。
9.根据权利要7所述的一种亚磷酸二甲酯关键工艺参数的预警定位系统,其特征在于,所述训练模块具体用于:
步骤一:基于LSTM的架构建立时序性机器学习模型;
步骤二:将样本库中的数据作为训练数据输入时序性机器学习模型,得到时序性机器学习模型的输出结果,记录从输入训练数据至得到输出结果的预测时长yTP
步骤三:将输出结果中的酯化反应温度预测值与同时序的实际酯化反应温度进行第一均方根误差的计算,将输出结果中的酯化真空度预测值与同时序的实际酯化真空度进行第二均方根误差的计算,将第一均方根误差和第二均方根误差相加得到预测精确度yACC
步骤四:基于评价公式EC对时序性机器学习模型的训练效果进行评价,所述评价公式为:
EC=w1yACC+w2yTP
其中,w1、w2分别为预设的精确度权重和预测时长权重,w1、w2的取值范围均为小于1的正数,yACC、yTP的取值范围均为正数;
步骤五:若EC的值大于预设训练阈值,则调整时序性机器学习模型中激活函数、神经网络层数、数据长度以及正则化的参数,重复步骤二至步骤四,直至EC的值不大于预设训练阈值时停止训练。
10.根据权利要求7所述的一种亚磷酸二甲酯关键工艺参数的预警定位系统,其特征在于,所述预警模块具体用于:
获取酯化反应温度、酯化真空度的预测结果,当预测结果超过预设预警值时获取当前关联位号对应的实时数据;
对实时数据进行相关性分析,根据实时数据与酯化反应温度、酯化真空度的相关程度进行排序,将排序结果对应的关联位号作为预警提示输出。
CN202110591042.7A 2021-05-28 2021-05-28 一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统 Active CN113539382B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110591042.7A CN113539382B (zh) 2021-05-28 2021-05-28 一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110591042.7A CN113539382B (zh) 2021-05-28 2021-05-28 一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统

Publications (2)

Publication Number Publication Date
CN113539382A true CN113539382A (zh) 2021-10-22
CN113539382B CN113539382B (zh) 2024-04-05

Family

ID=78124460

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110591042.7A Active CN113539382B (zh) 2021-05-28 2021-05-28 一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统

Country Status (1)

Country Link
CN (1) CN113539382B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114202124A (zh) * 2021-12-14 2022-03-18 浙江中控技术股份有限公司 一种环氧树脂合成反应停止时间的预测方法
CN114462309A (zh) * 2022-01-19 2022-05-10 苏州浪潮智能科技有限公司 一种物理仿真模型训练方法、装置、设备和存储介质
CN115494196A (zh) * 2022-09-19 2022-12-20 山东二叶制药有限公司 用于新舒巴坦酸提取工艺的品质保持方法
WO2024016552A1 (zh) * 2022-07-21 2024-01-25 浙江中控技术股份有限公司 流化床反应温度的确定方法、相关装置及计算机存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109814513A (zh) * 2019-03-20 2019-05-28 杭州辛孚能源科技有限公司 一种基于数据模型的催化裂化装置优化方法
US20200184131A1 (en) * 2018-06-27 2020-06-11 Dalian University Of Technology A method for prediction of key performance parameter of an aero-engine transition state acceleration process based on space reconstruction
CN111353631A (zh) * 2019-11-26 2020-06-30 国网山东省电力公司电力科学研究院 基于多层lstm的火电厂凝汽器真空度预测方法
CN111563706A (zh) * 2020-03-05 2020-08-21 河海大学 一种基于lstm网络的多变量物流货运量预测方法
CN111999657A (zh) * 2020-10-29 2020-11-27 北京航空航天大学 一种电动汽车的锂离子电池剩余寿命的行驶里程评估方法
CN112101480A (zh) * 2020-09-27 2020-12-18 西安交通大学 一种多变量聚类与融合的时间序列组合预测方法
CN112529234A (zh) * 2019-09-18 2021-03-19 上海交通大学 基于深度学习的地表水质预测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200184131A1 (en) * 2018-06-27 2020-06-11 Dalian University Of Technology A method for prediction of key performance parameter of an aero-engine transition state acceleration process based on space reconstruction
CN109814513A (zh) * 2019-03-20 2019-05-28 杭州辛孚能源科技有限公司 一种基于数据模型的催化裂化装置优化方法
CN112529234A (zh) * 2019-09-18 2021-03-19 上海交通大学 基于深度学习的地表水质预测方法
CN111353631A (zh) * 2019-11-26 2020-06-30 国网山东省电力公司电力科学研究院 基于多层lstm的火电厂凝汽器真空度预测方法
CN111563706A (zh) * 2020-03-05 2020-08-21 河海大学 一种基于lstm网络的多变量物流货运量预测方法
CN112101480A (zh) * 2020-09-27 2020-12-18 西安交通大学 一种多变量聚类与融合的时间序列组合预测方法
CN111999657A (zh) * 2020-10-29 2020-11-27 北京航空航天大学 一种电动汽车的锂离子电池剩余寿命的行驶里程评估方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李爱莲;张帅;: "基于LSTM网络的焦炉火道温度预测建模仿真", 计算机仿真, no. 06, pages 1 - 8 *
翁国庆: "基于聚类LSTM深度学习模型的主动配电网电能质量预测", 高技术通讯, no. 07, pages 1 - 13 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114202124A (zh) * 2021-12-14 2022-03-18 浙江中控技术股份有限公司 一种环氧树脂合成反应停止时间的预测方法
CN114462309A (zh) * 2022-01-19 2022-05-10 苏州浪潮智能科技有限公司 一种物理仿真模型训练方法、装置、设备和存储介质
CN114462309B (zh) * 2022-01-19 2023-11-17 苏州浪潮智能科技有限公司 一种物理仿真模型训练方法、装置、设备和存储介质
WO2024016552A1 (zh) * 2022-07-21 2024-01-25 浙江中控技术股份有限公司 流化床反应温度的确定方法、相关装置及计算机存储介质
CN115494196A (zh) * 2022-09-19 2022-12-20 山东二叶制药有限公司 用于新舒巴坦酸提取工艺的品质保持方法

Also Published As

Publication number Publication date
CN113539382B (zh) 2024-04-05

Similar Documents

Publication Publication Date Title
CN113539382B (zh) 一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统
CN111444953B (zh) 一种基于改进粒子群优化算法的传感器故障监测方法
CN105302848A (zh) 一种设备智能预警系统的评估值校准方法
CN109919356B (zh) 一种基于bp神经网络区间需水预测方法
US20240144043A1 (en) Prediction model for predicting product quality parameter values
CN109155018A (zh) 用于kpi性能分析的具有无约束因变量的mpc
CN116380445B (zh) 基于振动波形的设备状态诊断方法及相关装置
CN113065703A (zh) 一种结合多种模型的时间序列预测方法
CN109240276B (zh) 基于故障敏感主元选择的多块pca故障监测方法
CN115878603A (zh) 基于k近邻算法与gan网络的水质缺失数据插补算法
CN113408659A (zh) 一种基于数据挖掘的建筑能耗集成分析方法
CN114329347B (zh) 一种电能表计量误差预测方法、装置及存储介质
CN114962239A (zh) 一种基于智能物联网的设备故障检测方法
CN114139446A (zh) 一种用于特殊精馏过程组分在线检测软测量方法
CN114819178A (zh) 一种铁路施工进度指标预测和在线更新方法
CN112541296A (zh) 一种基于pso-lssvm的so2预测方法
CN116662925A (zh) 一种基于加权稀疏神经网络工业过程软测量方法
KR102110319B1 (ko) 학습 데이터 생성 시스템
CN116720080A (zh) 同源气象要素融合检验方法
CN116363843A (zh) 试验室设备预警系统
CN115828744A (zh) 一种白光发光二极管故障在线诊断与寿命预测方法
CN113887571A (zh) 改进svr算法的电子设备故障预测方法
CN115906602A (zh) 一种基于深度学习Transformer自编码的风电机组塔筒倾倒监测评价方法
CN117829822B (zh) 一种电力变压器故障预警方法及系统
CN117076260A (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