CN115454009A - Model Predictive Control Method of Component Formulation for Chemical Production - Google Patents

Model Predictive Control Method of Component Formulation for Chemical Production Download PDF

Info

Publication number
CN115454009A
CN115454009A CN202211408678.4A CN202211408678A CN115454009A CN 115454009 A CN115454009 A CN 115454009A CN 202211408678 A CN202211408678 A CN 202211408678A CN 115454009 A CN115454009 A CN 115454009A
Authority
CN
China
Prior art keywords
distillation
temperature
distillation range
point
heating element
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.)
Pending
Application number
CN202211408678.4A
Other languages
Chinese (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.)
Changshu Institute of Technology
Pengchen New Material Technology Co Ltd
Original Assignee
Changshu Institute of Technology
Pengchen New Material 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 Changshu Institute of Technology, Pengchen New Material Technology Co Ltd filed Critical Changshu Institute of Technology
Priority to CN202211408678.4A priority Critical patent/CN115454009A/en
Publication of CN115454009A publication Critical patent/CN115454009A/en
Pending legal-status Critical Current

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
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/418Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
    • G05B19/41865Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by job scheduling, process planning, material flow
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/32Operator till task planning
    • G05B2219/32252Scheduling production, machining, job shop

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Quality & Reliability (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention relates to the technical field of chemical control, in particular to a component distribution model predictive control method for chemical production. The method obtains the reaction characteristics of each distillation range point, and matches the distillation range points in the distillation process under different control parameters according to the reaction characteristics to obtain a plurality of matched pairs. And further obtaining a switchable index according to the reaction characteristic difference between the distillation range points in the matching pair, selecting the matching pair based on the switchable index to construct an initial training set, and amplifying the initial training set according to the matching data to obtain a prediction model for training the long-short term memory production group to prepare, so as to predict real-time data and control container parameters according to the prediction parameters. The invention predicts the parameters to be adjusted of the future distillation range point by constructing the component distribution model, controls the parameters in time and ensures the production efficiency.

Description

用于化工生产的组分配制模型预测控制方法Model Predictive Control Method of Component Formulation for Chemical Production

技术领域technical field

本发明涉及化工控制技术领域,具体涉及用于化工生产的组分配制模型预测控制方法。The invention relates to the technical field of chemical control, in particular to a predictive control method for component allocation models used in chemical production.

背景技术Background technique

高沸点芳烃溶剂是以重整芳烃为原料,按国际上专用芳烃溶剂标准研制生产的,具有溶解力强、毒性低、气味小、沸点高、挥发慢、不含水和烯烃、不含氯和重金属、化学物理性能稳定及流平性好等特点。重芳烃是一种粗品,前期的性质不同,可以提取250°C轻馏份,结合馏程更高的尾油(馏程约260〜380°C),对C10芳烃进行合理利用,生产高沸点芳烃。一般的生产过程中需要经催化加氢,然后进行精馏分离,但目前这种配制参数是依靠人工经验调整的,容易出现先品质不稳定的情况,在后期也一直需要人力监控,因此在组分配制过程中需要一种智能判别是否需要改进工艺,进行加氢干预的系统,以提高蒸出效率。High-boiling aromatic solvents are developed and produced based on reformed aromatics as raw materials, according to international special aromatic solvent standards. They have strong solvency, low toxicity, low odor, high boiling point, slow volatilization, no water and olefins, no chlorine and heavy metals , stable chemical and physical properties and good leveling. Heavy aromatics are a kind of crude product with different properties in the early stage. Light fractions at 250 ° C can be extracted, combined with tail oil with a higher distillation range (distillation range about 260 ~ 380 ° C), rational utilization of C10 aromatics, production of high boiling point Aromatics. In the general production process, catalytic hydrogenation is required, followed by rectification and separation. However, at present, this preparation parameter is adjusted by manual experience, which is prone to unstable quality at first, and human monitoring is always required in the later stage. Therefore, in the assembly In the process of dispensing, a system is needed to intelligently judge whether the process needs to be improved and to intervene in hydrogenation to improve the efficiency of distillation.

发明内容Contents of the invention

为了解决上述技术问题,本发明的目的在于提供一种用于化工生产的组分配制模型预测控制方法,所采用的技术方案具体如下:In order to solve the above-mentioned technical problems, the object of the present invention is to provide a component formulation model predictive control method for chemical production, and the adopted technical scheme is as follows:

本发明提出了一种用于化工生产的组分配制模型预测控制方法,所述方法包括:The present invention proposes a method for predictive control of component formulation model for chemical production, said method comprising:

采集蒸馏过程中的每个馏程点的设定温度值、蒸汽温度和容器内加热元件温度;获得容器内反应物的外观特征混合编码;根据采样过程下的蒸汽温度、加热元件温度和容器的体积空速构建短期反应描述子;Collect the set temperature value, steam temperature and heating element temperature in the container at each distillation point in the distillation process; obtain the mixed code of the appearance characteristics of the reactants in the container; according to the steam temperature, heating element temperature and container Volumetric space velocity constructs short-term response descriptors;

根据不同控制参数下蒸馏过程中馏程点之间的外观特征混合编码差异和短期反应描述子差异对不同馏程点进行匹配,获得多个匹配对;Match different distillation range points according to the mixed coding difference of appearance characteristics and short-term response descriptor difference in the distillation process under different control parameters, and obtain multiple matching pairs;

根据匹配对中两个馏程点之间的设定温度值差异、蒸汽温度差异和加热元件温度差异获得可切换指数;选择可切换指数最大的预设数量个匹配对中的数据作为初始训练集,根据初始训练集中互相匹配的馏程点数据对初始训练集进行增广,获得训练集,根据训练集训练长短期记忆生产组分配制预测模型;The switchable index is obtained according to the set temperature value difference, steam temperature difference and heating element temperature difference between two distillation range points in the matching pair; select the data in the preset number of matching pairs with the largest switchable index as the initial training set , augment the initial training set according to the matching distillation point data in the initial training set to obtain a training set, and train the long-short-term memory production group to prepare a prediction model according to the training set;

将实时蒸馏过程中实时馏程点的实时蒸汽温度和实时加热元件温度输入长短期记忆生产组分配制预测模型中,获得下一个馏程点对应的设定温度值和体积空速,并根据下一个馏程点对应的设定温度值和体积空速对容器参数进行控制。Input the real-time steam temperature and real-time heating element temperature at the real-time distillation point in the real-time distillation process into the long-short-term memory production component allocation prediction model to obtain the set temperature value and volume space velocity corresponding to the next distillation point, and according to the following The set temperature value and volume space velocity corresponding to a distillation range point control the vessel parameters.

进一步地,所述获得容器内反应物的外观特征混合编码包括:Further, said obtaining the mixed coding of the appearance characteristics of the reactants in the container includes:

获取容器内反应物的外观特征的描写词语,并根据描写词语构建TF-IDF外观编码;对容器反应物的颜色信息进行采集,以RGB颜色空间中的分量值构建外观颜色编码;将TF-IDF外观编码和外观颜色编码合并,获得外观特征混合编码。Obtain the description words of the appearance characteristics of the reactants in the container, and construct the TF-IDF appearance code according to the description words; collect the color information of the container reactants, and construct the appearance color code with the component values in the RGB color space; Appearance coding and appearance color coding are combined to obtain a mixed coding of appearance features.

进一步地,所述获得外观特征混合编码包括:Further, said obtaining appearance feature mixed coding includes:

将外观颜色编码的维度进行扩充,获得第一扩充编码;将第一扩充编码与TF-IDF外观编码合并,获得第二扩充编码;将第二扩充编码降维至预设维数,获得外观特征混合编码。Expand the dimension of the appearance color code to obtain the first extended code; merge the first extended code with the TF-IDF appearance code to obtain the second extended code; reduce the dimensionality of the second extended code to the preset dimension to obtain the appearance features Mixed encoding.

进一步地,所述根据采样过程下的蒸汽温度、加热元件温度和容器的体积空速构建短期反应描述子包括:Further, said constructing a short-term reaction descriptor according to the steam temperature under the sampling process, the temperature of the heating element and the volume space velocity of the container includes:

所述短期反应描述子由蒸馏过程中每个馏程点下的蒸汽温度序列中值、加热元件温度序列中值、蒸汽温度序列均值、加热元件温度序列均值、蒸汽温度序列初始值和加热元件温度序列初始值构成。The short-term response descriptor consists of the median value of the steam temperature sequence, the median value of the heating element temperature sequence, the average value of the steam temperature sequence, the average value of the heating element temperature sequence, the initial value of the steam temperature sequence, and the temperature of the heating element at each distillation point in the distillation process. Sequence initial value composition.

进一步地,所述获得多个匹配对的方法包括:Further, the method for obtaining multiple matching pairs includes:

根据不同反应参数下的不同蒸馏过程之间的外观特征混合编码差异和短期反应描述子差异构建匹配距离;根据匹配距离利用KM匹配算法对不同蒸馏过程进行匹配,获得多个匹配对。The matching distance was constructed according to the mixed coding difference of appearance features and the short-term reaction descriptor difference between different distillation processes under different reaction parameters; according to the matching distance, the KM matching algorithm was used to match different distillation processes to obtain multiple matching pairs.

进一步地,所述根据不同反应参数下的不同蒸馏过程之间的外观特征混合编码差异和短期反应描述子差异构建匹配距离包括:Further, the construction of the matching distance based on the mixed encoding difference and short-term reaction descriptor difference between different distillation processes under different reaction parameters includes:

根据匹配距离公式获得匹配距离,匹配距离公式包括:The matching distance is obtained according to the matching distance formula, and the matching distance formula includes:

Figure DEST_PATH_IMAGE001
Figure DEST_PATH_IMAGE001

其中,

Figure 531474DEST_PATH_IMAGE002
为馏程点X和馏程点Y之间的匹配距离,
Figure 767283DEST_PATH_IMAGE003
为馏程点X和馏程点Y之间的短期反应描述子相似度,
Figure 840282DEST_PATH_IMAGE004
为馏程点X和馏程点Y之间的外观特征混合编码相似度,
Figure 870555DEST_PATH_IMAGE005
为馏程点X和馏程点Y在整个蒸馏过程中对应的馏程点序号的差异绝对值;若
Figure 29003DEST_PATH_IMAGE006
,则
Figure 802924DEST_PATH_IMAGE005
为无穷大。in,
Figure 531474DEST_PATH_IMAGE002
is the matching distance between distillation point X and distillation point Y,
Figure 767283DEST_PATH_IMAGE003
is the short-term response descriptor similarity between distillation point X and distillation point Y,
Figure 840282DEST_PATH_IMAGE004
Coding similarity for the mixture of appearance features between distillation point X and distillation point Y,
Figure 870555DEST_PATH_IMAGE005
is the absolute value of the difference between the distillation point numbers corresponding to the distillation point X and the distillation point Y in the whole distillation process; if
Figure 29003DEST_PATH_IMAGE006
,but
Figure 802924DEST_PATH_IMAGE005
for infinity.

进一步地,所述根据匹配对中两个蒸馏过程中的设定温度值差异、蒸汽温度差异和加热元件温度差异获得可切换指数包括:Further, said obtaining the switchable index according to the set temperature value difference, the steam temperature difference and the heating element temperature difference in the two distillation processes in the matching pair includes:

根据可切换指数公式获得可切换指数,可切换指数公式包括:The switchable index is obtained according to the switchable index formula, which includes:

Figure 730429DEST_PATH_IMAGE007
Figure 730429DEST_PATH_IMAGE007

其中,

Figure 931603DEST_PATH_IMAGE008
为馏程点X和馏程点Y之间的可切换指数,
Figure 842928DEST_PATH_IMAGE009
为馏程点X和馏程点Y之间的设定温度值的差值绝对值,
Figure 889381DEST_PATH_IMAGE002
为馏程点X和馏程点Y之间的匹配距离,
Figure 936971DEST_PATH_IMAGE010
为馏程点X和馏程点Y之间蒸汽温度序列的动态时间规整距离,
Figure 309047DEST_PATH_IMAGE011
为馏程点X和馏程点Y之间加热元件温度序列的动态时间规整距离。in,
Figure 931603DEST_PATH_IMAGE008
is the switchable index between distillation point X and distillation point Y,
Figure 842928DEST_PATH_IMAGE009
is the absolute value of the difference between the set temperature values between the distillation point X and the distillation point Y,
Figure 889381DEST_PATH_IMAGE002
is the matching distance between distillation point X and distillation point Y,
Figure 936971DEST_PATH_IMAGE010
is the dynamic time-regularized distance of the steam temperature sequence between distillation point X and distillation point Y,
Figure 309047DEST_PATH_IMAGE011
is the dynamic time warping distance of the temperature sequence of the heating element between the distillation point X and the distillation point Y.

进一步地,所述根据初始训练集中互相匹配的蒸馏过程数据对初始训练集进行增广包括:Further, said augmenting the initial training set according to the matching distillation process data in the initial training set includes:

将目标馏程点的体积空速、设定温度、蒸汽温度和加热元件温度替换为匹配馏程点所属的蒸馏过程中对应馏程点的体积空速、设定温度、蒸汽温度和加热元件温度,获得第一增广训练数据;Replace the volume space velocity, set temperature, vapor temperature, and heating element temperature of the target distillation point with the volume space velocity, set temperature, vapor temperature, and heating element temperature of the corresponding distillation point in the distillation process to which the matching distillation point belongs , to obtain the first augmented training data;

将目标馏程点下一个馏程点的体积空速、设定温度、蒸汽温度和加热元件温度替换为匹配馏程点所属的蒸馏过程中对应馏程点的体积空速、设定温度、蒸汽温度和加热元件温度,获得第二增广训练数据;Replace the volume space velocity, set temperature, steam temperature and heating element temperature of the distillation point next to the target distillation point with the volume space velocity, set temperature, steam temperature of the corresponding distillation point in the distillation process to which the matching distillation point belongs temperature and heating element temperature to obtain the second augmented training data;

分别将目标馏程点和目标馏程点下一个馏程点的体积空速、设定温度、蒸汽温度和加热元件温度同时替换为匹配馏程点所属的蒸馏过程中对应馏程点的体积空速、设定温度、蒸汽温度和容器内加热元件温度,获得第三增广训练数据;The volume space velocity, set temperature, steam temperature and heating element temperature of the target distillation point and the distillation point next to the target distillation point are respectively replaced with the volume space of the corresponding distillation point in the distillation process to which the matching distillation point belongs. Speed, set temperature, steam temperature and heating element temperature in the container to obtain the third augmented training data;

将每个馏程点的第一增广训练数据、第二增广训练数据和第三增广训练数据并获得训练集。The first augmented training data, the second augmented training data and the third augmented training data of each distillation point are used to obtain a training set.

本发明具有如下有益效果:The present invention has following beneficial effect:

本发明实施例通过构建短期反应描述子和外观特征混合编码用于表征蒸馏过程中的反应特征。根据反应特征即可对不同参数下蒸馏过程的馏程点进行匹配,进而通过选取可切换指数大的匹配对作为训练数据,并根据匹配对之间的数据对数据进行增广,使得后续的长短期记忆生产组分配制预测模型的训练过程拥有足够的训练基础,并且长短期记忆生产组分配制预测模型可根据增广数据输出一个优选的结果,保证了预测数据的可参考性,进而通过长短期记忆生产组分配制预测模型预测出未来时刻的设定温度值和体积空速,并对容器内的参数进行控制。实现了及时对容器参数的控制,保证了生产效率。In the embodiment of the present invention, a short-term reaction descriptor and a hybrid code of appearance characteristics are constructed to characterize the reaction characteristics in the distillation process. According to the reaction characteristics, the distillation range points of the distillation process under different parameters can be matched, and then by selecting the matching pair with a large switchable index as the training data, and augmenting the data according to the data between the matching pairs, the subsequent long-term The training process of the short-term memory production composition allocation prediction model has sufficient training basis, and the long-term short-term memory production composition allocation prediction model can output an optimal result according to the augmented data, ensuring the referenceability of the prediction data, and then through the long-term The short-term memory production group allocation forecast model predicts the set temperature value and volume space velocity in the future, and controls the parameters in the container. Realize timely control of container parameters and ensure production efficiency.

附图说明Description of drawings

为了更清楚地说明本发明实施例或现有技术中的技术方案和优点,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它附图。In order to more clearly illustrate the technical solutions and advantages in the embodiments of the present invention or in the prior art, the following will briefly introduce the accompanying drawings that need to be used in the description of the embodiments or the prior art. Apparently, the appended The drawings are only some embodiments of the present invention, and those skilled in the art can also obtain other drawings based on these drawings without creative work.

图1为本发明一个实施例所提供的一种用于化工生产的组分配制模型预测控制方法流程图;Fig. 1 is a flow chart of a component formulation model predictive control method for chemical production provided by an embodiment of the present invention;

图2为本发明一个实施例所提供的一种颜色信息采集过程示意图。Fig. 2 is a schematic diagram of a color information collection process provided by an embodiment of the present invention.

具体实施方式detailed description

为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例,对依据本发明提出的一种用于化工生产的组分配制模型预测控制方法,其具体实施方式、结构、特征及其功效,详细说明如下。在下述说明中,不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外,一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。In order to further explain the technical means and effects of the present invention to achieve the intended purpose of the invention, the following in conjunction with the accompanying drawings and preferred embodiments, for a component allocation model predictive control method for chemical production proposed according to the present invention, Its specific implementation, structure, feature and effect thereof are described in detail as follows. In the following description, different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. Furthermore, the particular features, structures or characteristics of one or more embodiments may be combined in any suitable manner.

除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the technical field of the invention.

下面结合附图具体的说明本发明所提供的一种用于化工生产的组分配制模型预测控制方法的具体方案。A specific scheme of a component formulation model predictive control method for chemical production provided by the present invention will be described below in conjunction with the accompanying drawings.

请参阅图1,其示出了本发明一个实施例提供的一种用于化工生产的组分配制模型预测控制方法流程图,该方法包括:Please refer to Fig. 1, which shows a flow chart of a component allocation model predictive control method for chemical production provided by an embodiment of the present invention, the method includes:

步骤S1:采集蒸馏过程中的每个馏程点的设定温度值、蒸汽温度和容器内加热元件温度;获得容器内反应物的外观特征混合编码;根据采样过程下的蒸汽温度、加热元件温度和容器的体积空速构建短期反应描述子。Step S1: Collect the set temperature value, steam temperature and heating element temperature of each distillation point in the distillation process; obtain the mixed code of the appearance characteristics of the reactants in the container; according to the steam temperature and heating element temperature in the sampling process and the volumetric space velocity of the container to construct a short-term reaction descriptor.

本发明实施例以CIO重芳烃和经切割提取轻组分后的CIO尾油按照体积比2:1混合。进行加氢轻质化的组分配制生产为例,对于组分配制,即分配加氢过程。在本发明实施例中其氢压为5MPa,以某种催化工艺为例,对于整个蒸馏过程,要实时控制蒸出的蒸汽温度,调节的节点分别在总体积的0%(最初馏程)、10%、20%、30%、40%、50%、60%、70%、80%、90%、95%以及终馏点(干点),即在本发明实施例中,一个蒸馏过程存在12个馏程点。在本发明实施例中对于最初馏程,设定温度值预设为98摄氏度,以此类推,一般是温度越来越高,在此不做限定,即每个馏程点均对应一个设定温度值。In the embodiment of the present invention, CIO heavy aromatics and CIO tail oil after cutting and extracting light components are mixed according to the volume ratio of 2:1. Take the hydrogenation and lightening component preparation production as an example, for the component preparation, that is, the distribution hydrogenation process. In the embodiment of the present invention, its hydrogen pressure is 5MPa. Taking a certain catalytic process as an example, for the entire distillation process, the temperature of the steam that is evaporated must be controlled in real time, and the nodes to be adjusted are respectively 0% of the total volume (initial distillation range), 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, 95% and final boiling point (dry point), that is, in the embodiment of the present invention, a distillation process exists 12 distillation range points. In the embodiment of the present invention, for the initial distillation range, the set temperature value is preset to 98 degrees Celsius, and so on, generally the temperature is getting higher and higher, which is not limited here, that is, each distillation range point corresponds to a set temperature value.

不同反应参数下的蒸馏过程的加热方式和控制方式不同,以及加氢程度不同,因此需要实时测量蒸馏过程下容器内的温度信息,具体容器内的温度信息包括蒸汽温度和加热元件温度。因为蒸汽温度之间有一定的数值迟滞,即反应容器加热元件温度与蒸汽温度之间的数据存在吸热、比热容的联合效应,存在一定的迟滞,即温度阻尼,这个特征与重芳烃原料特性、加氢反应的速率以及两相物质的内涵特性有关,因此在采集温度信息时需要采集蒸汽温度和加热元件温度两种温度信息。The heating and control methods of the distillation process under different reaction parameters are different, and the degree of hydrogenation is different. Therefore, it is necessary to measure the temperature information in the container under the distillation process in real time. The specific temperature information in the container includes the steam temperature and the temperature of the heating element. Because there is a certain numerical hysteresis between the steam temperature, that is, there is a joint effect of heat absorption and specific heat capacity in the data between the heating element temperature of the reaction vessel and the steam temperature, and there is a certain hysteresis, that is, temperature damping. This feature is related to the characteristics of heavy aromatics raw materials, The rate of the hydrogenation reaction is related to the internal characteristics of the two-phase substance, so when collecting temperature information, two temperature information, the steam temperature and the heating element temperature, need to be collected.

在本发明实施例中,以1Hz的速率记录蒸汽温度T_A和反应容器内加热元件温度T_B。In the embodiment of the present invention, the steam temperature T_A and the heating element temperature T_B in the reaction vessel were recorded at a rate of 1 Hz.

在本发明实施例中,对于每个馏程点之间,都有若干个蒸汽温度记录数值,对这些数值进行最近邻重采样,保证每个馏程点之间存在400个数据。即每个馏程点的数据被处理为一种定长的温度记录值,这个数据是T_A和T_B的两个数据定长序列。即每个馏程点对应的每种温度序列是由400个数据点信息构成的序列。In the embodiment of the present invention, for each distillation range point, there are several recorded steam temperature values, and nearest neighbor resampling is performed on these values to ensure that there are 400 data between each distillation range point. That is, the data of each distillation point is processed as a fixed-length temperature record value, and this data is two data fixed-length sequences of T_A and T_B. That is, each temperature sequence corresponding to each distillation point is a sequence composed of 400 data point information.

在蒸馏过程中,每个馏程点因为反应参数的影响,容器内的反应物会存在不同的形态特征,例如反应物的颜色、透明程度等信息。因此为了进一步描述蒸馏过程还需要获得蒸馏过程中每个蒸馏点下容器内的反应物外观特征混合编码,利用外观特征混合编码表征每个馏程点下反应物的形态信息,具体包括:During the distillation process, due to the influence of reaction parameters at each distillation point, the reactants in the container will have different morphological characteristics, such as the color and transparency of the reactants. Therefore, in order to further describe the distillation process, it is necessary to obtain the mixed codes of the appearance characteristics of the reactants in the container at each distillation point in the distillation process, and use the mixed codes of appearance characteristics to characterize the morphological information of the reactants at each distillation point, including:

对物性数据进行分析,获取容器内反应物外观特征的描写词语,并通过常规词语来表达:外观词语可以为微黄色、澄清等,为了保持业界词语统一,实现文字的去歧义,在本发明实施例中使用英文单词,例如Bright & Clear中对Bright和Clear进行分词,从而得到基于词库分词后的TF-IDF外观编码。需要说明的是外观编码的维度和所分词的词库大小有关,具体编码方法为本领域技术人员熟知的技术手段,在此不做赘述;获取容器内反应物外观特征的描写词语的方法可通过目测或者神经网络进行获取,在此不做限定,可根据具体实施场景自由选择获取方式。Analyze the physical property data, obtain the description words of the appearance characteristics of the reactants in the container, and express them through conventional words: the appearance words can be yellowish, clear, etc. In order to maintain the unity of words in the industry and realize the disambiguation of words, in the implementation of the present invention In the example, English words are used. For example, in Bright & Clear, Bright and Clear are segmented to obtain TF-IDF appearance codes based on word segmentation in the thesaurus. It should be noted that the dimension of the appearance coding is related to the size of the lexicon of the segmented words. The specific coding method is a technical means well known to those skilled in the art, and will not be repeated here; the method of obtaining the description words of the appearance characteristics of the reactants in the container can be obtained through Visual inspection or neural network acquisition is not limited here, and the acquisition method can be freely selected according to the specific implementation scenario.

进一步对容器反应物的颜色信息进行采集,以RGB颜色空间中的分量值构建外观颜色编码。请参阅图2,其示出了本发明实施例提供的一种本发明实施例提供的一种颜色信息采集过程示意图,本发明实施例使用ADJD-S313-QR999RGB分量色度传感器,在恒定光源亮度下测得外观RGB数值,在反应管路中构建旁路或直接安装观察窗,即可测量。其中,本实施例使用avago公司型号为adjd-s313-qr999的CMOS数字色彩传感器,光源1提供背光,得到透射的颜色参考,光源2提供漫反射颜色参考,RGB色度测量后,基于传感器的最大量程,为RGB分量归一化到0~1的三个浮点值,即每个馏程点的外观颜色编码为三个颜色通道分量的归一化数据组成的编码。The color information of the reactants in the container is further collected, and the appearance color coding is constructed with the component values in the RGB color space. Please refer to Fig. 2, which shows a schematic diagram of a color information acquisition process provided by an embodiment of the present invention. The embodiment of the present invention uses the ADJD-S313-QR999RGB component chromaticity sensor, under constant light source brightness The RGB value of the appearance is measured below, and the bypass can be built in the reaction pipeline or the observation window can be directly installed to measure. Wherein, this embodiment uses the CMOS digital color sensor of avago company model as adjd-s313-qr999, light source 1 provides backlight, obtains the color reference of transmission, light source 2 provides diffuse reflection color reference, after RGB chromaticity measurement, based on the sensor's maximum The range is three floating-point values normalized to 0~1 for the RGB components, that is, the appearance color coding of each distillation point is coded by the normalized data of the three color channel components.

进一步将TF-IDF外观编码和外观颜色编码合并,获得外观特征混合编码,具体包括:Further merge the TF-IDF appearance coding and appearance color coding to obtain a mixed coding of appearance features, including:

a.首先对外观颜色编码C内的多个归一化后的RGB分量进行维度扩增,得到12*3维的高维向量第一扩充编码H_0。a. First, dimensionally expand the multiple normalized RGB components in the appearance color code C to obtain the first extended code H_0 of a 12*3-dimensional high-dimensional vector.

b.然后与文字TF-IDF外观编码合并,得到12*3+X维的高维向量第二扩充编码H_1,其中X为TF-IDF外观编码的维度数。b. Then merge it with the text TF-IDF appearance code to obtain the second extended code H_1 of a high-dimensional vector of 12*3+X dimensions, where X is the number of dimensions of the TF-IDF appearance code.

c.进一步将第二扩充编码降维至预设维数,获得外观特征混合编码。在本发明实施例中,基于各次生产记录中的高维向量H_1,构建高维空间,基于PCA算法降维到10维,获得外观特征混合编码。即本发明实施例中预设维数为10维。c. Further reduce the dimensionality of the second extended code to the preset dimension to obtain a mixed code of appearance features. In the embodiment of the present invention, a high-dimensional space is constructed based on the high-dimensional vector H_1 in each production record, and the dimension is reduced to 10 dimensions based on the PCA algorithm to obtain a mixed code of appearance features. That is, the preset dimension in the embodiment of the present invention is 10 dimensions.

i. 具体地,从经验值上来说,10维是精度可接受的维度数,实施者可以基于该维度数继续增加维度,从而保证数据精度,降维后的向量主要用来表示未经控制时的各个蒸出过程的物性外观数据。i. Specifically, from an empirical point of view, 10 dimensions is the number of dimensions with acceptable accuracy. The implementer can continue to increase the dimension based on this number of dimensions, so as to ensure the accuracy of the data. The vector after dimension reduction is mainly used to represent the uncontrolled time Physical properties and appearance data of each evaporation process.

ii.至此,基于前期的数据得到降维后的外观特征混合编码H。ii. So far, based on the previous data, the appearance feature hybrid code H after dimensionality reduction is obtained.

外观特征混合编码H的作用是针对此次蒸馏过程,对主观信息和客观信息联合表示,从而避免人工分析的主观误差,同时提高数据在存管过程中的准确性和客观性。The role of the appearance feature mixed code H is to jointly express the subjective information and objective information for this distillation process, so as to avoid the subjective error of manual analysis, and at the same time improve the accuracy and objectivity of the data in the storage process.

根据各个生产数据,根据采样过程下的蒸汽温度、加热元件温度和容器的体积空速构建短期反应描述子,具体包括:According to each production data, short-term reaction descriptors are constructed according to the steam temperature, heating element temperature and volume space velocity of the container under the sampling process, including:

以氢压6MPa为例,在生产过程中,催化剂的各项配比会微调、且蒸出温度也有微调,对应地,以氢油体积比1000为例,体积空速V可能会在1.8到1.5之间,各个馏程点的蒸出温度会有不同的目标值。一般情况下,V越高表示可假设的催化剂活性愈高,对应的装置处理能力越大。需要说明的是,任何反应容器和反应原理都有边际效应,V不能无限提高,对于测量的容器装置,V比较大意味着单位时间里通过催化剂的原料多,原料在催化剂上的停留时间短,反应深度浅;相反,V小意味着反应时间长,一般情况下降低V对于提高反应的转化率是有利的。但是,过低的V意味着在相同处理量的情况下需要的催化剂数量较多,在馏程中,有可能变相地降低了反应效率在经济上是不合理的。因此,V代表了此时投产时,催化剂微调后的经验情况,不一定是最佳参数。因此,温度参数和体积空速V是短期反应的上下文特征的影响因素。对于V,由于已经投产,所记录的V一般是通过微调和经验调优后的。因此根据采样过程下的蒸汽温度、加热元件温度和容器的体积空速构建短期反应描述子,具体包括:Taking the hydrogen pressure of 6MPa as an example, during the production process, the proportions of the catalyst will be fine-tuned, and the steaming temperature will also be fine-tuned. Correspondingly, taking the hydrogen-oil volume ratio of 1000 as an example, the volumetric space velocity V may be between 1.8 and 1.5 Between, the distillation temperature of each distillation range point will have different target values. In general, the higher the V, the higher the assumed catalyst activity, and the greater the processing capacity of the corresponding device. It should be noted that any reaction vessel and reaction principle have marginal effects, and V cannot be increased infinitely. For the measured container device, a relatively large V means that there are more raw materials passing through the catalyst per unit time, and the residence time of raw materials on the catalyst is short. The reaction depth is shallow; on the contrary, a small V means a long reaction time, and generally lowering V is beneficial to improving the conversion rate of the reaction. However, too low V means that the amount of catalyst required under the same treatment amount is large, and in the distillation range, it is economically unreasonable to reduce the reaction efficiency in a disguised manner. Therefore, V represents the experience of fine-tuning the catalyst when it is put into production at this time, and it is not necessarily the best parameter. Therefore, temperature parameters and volumetric space velocity V are influencing factors for the contextual characteristics of short-term reactions. As for V, since it has already been put into production, the recorded V is generally adjusted through fine-tuning and experience. Therefore, a short-term response descriptor is constructed according to the steam temperature, heating element temperature, and container volume space velocity during the sampling process, including:

短期反应描述子由蒸馏过程中每个馏程点下的蒸汽温度序列中值、加热元件温度序列中值、蒸汽温度序列均值、加热元件温度序列均值、蒸汽温度序列初始值和加热元件温度序列初始值构成。即短期反应描述子为多个特征参数组成的向量,其数学表现形式为

Figure 442088DEST_PATH_IMAGE012
,其中
Figure 557811DEST_PATH_IMAGE013
为蒸汽温度序列中值,
Figure 459908DEST_PATH_IMAGE014
为加热元件温度序列中值,
Figure 2885DEST_PATH_IMAGE015
为蒸汽温度序列均值,
Figure 357643DEST_PATH_IMAGE016
为加热元件温度序列均值,
Figure 542637DEST_PATH_IMAGE017
为蒸汽温度序列初始值,
Figure 322678DEST_PATH_IMAGE018
为加热元件温度序列初始值。The short-term reaction descriptor consists of the median value of the steam temperature series, the median value of the heating element temperature series, the mean value of the steam temperature series, the mean value of the heating element temperature series, the initial value of the steam temperature series and the initial value of the heating element temperature series at each distillation point in the distillation process. value composition. That is, the short-term response descriptor is a vector composed of multiple characteristic parameters, and its mathematical expression is
Figure 442088DEST_PATH_IMAGE012
,in
Figure 557811DEST_PATH_IMAGE013
is the median value of the steam temperature series,
Figure 459908DEST_PATH_IMAGE014
is the median value of the heating element temperature sequence,
Figure 2885DEST_PATH_IMAGE015
is the mean value of the steam temperature series,
Figure 357643DEST_PATH_IMAGE016
is the mean value of the heating element temperature sequence,
Figure 542637DEST_PATH_IMAGE017
is the initial value of the steam temperature series,
Figure 322678DEST_PATH_IMAGE018
It is the initial value of the heating element temperature sequence.

步骤S2:根据不同控制参数下蒸馏过程中每个馏程点之间的外观特征混合编码差异和短期反应描述子差异对不同馏程点进行匹配,获得多个匹配对。Step S2: Match different distillation range points according to the mixed coding difference of appearance characteristics and the difference of short-term response descriptors between each distillation range point in the distillation process under different control parameters, and obtain multiple matching pairs.

首先不同体积空速V、催化剂的调整和温度的调整,会产生不同的产后数据,因此,对于短期反应模式描述子Q和馏程点的颜色参考值C,可以基于K-M算子的搜索来按照如下策略做类型分配:搜索每个短期反应模式与下一个反应馏程的反应模式的异同和匹配距离,可以为后续控制的介入程度以及温度推荐值进行模型参数的构建。即以每个馏程点的反应过程作为一个短期反应过程,将不同反应参数下蒸馏过程中馏程点之间进行匹配。例如将反应模式A下的第2个馏程点作为目标馏程点,将目标馏程点与反应模式B、反应模式D等其他所有反应模式下的所有馏程点进行匹配。First of all, different volume space velocity V, catalyst adjustment and temperature adjustment will produce different post-production data. Therefore, for the short-term reaction mode descriptor Q and the color reference value C of the distillation point, it can be searched based on the K-M operator according to The following strategy is used for type allocation: search for the similarities and differences and matching distances between each short-term reaction mode and the reaction mode of the next reaction distillation range, and can construct model parameters for the degree of subsequent control intervention and recommended temperature values. That is, the reaction process of each distillation point is regarded as a short-term reaction process, and the distillation points in the distillation process under different reaction parameters are matched. For example, the second distillation point in reaction mode A is used as the target distillation point, and the target distillation point is matched with all distillation points in reaction mode B, reaction mode D and other reaction modes.

根据匹配距离公式获得匹配距离,匹配距离公式包括:The matching distance is obtained according to the matching distance formula, and the matching distance formula includes:

Figure 36556DEST_PATH_IMAGE001
Figure 36556DEST_PATH_IMAGE001

其中,

Figure 144189DEST_PATH_IMAGE002
为馏程点X和馏程点Y之间的匹配距离,
Figure 601715DEST_PATH_IMAGE003
为馏程点X和馏程点Y之间的短期反应描述子相似度,
Figure 947246DEST_PATH_IMAGE004
为馏程点X和馏程点Y之间的外观特征混合编码相似度,
Figure 97605DEST_PATH_IMAGE005
为馏程点X和馏程点Y在整个蒸馏过程中对应的馏程点序号的差异绝对值;若
Figure 426955DEST_PATH_IMAGE006
,则
Figure 422593DEST_PATH_IMAGE005
为无穷大。in,
Figure 144189DEST_PATH_IMAGE002
is the matching distance between distillation point X and distillation point Y,
Figure 601715DEST_PATH_IMAGE003
is the short-term response descriptor similarity between distillation point X and distillation point Y,
Figure 947246DEST_PATH_IMAGE004
Coding similarity for the mixture of appearance features between distillation point X and distillation point Y,
Figure 97605DEST_PATH_IMAGE005
is the absolute value of the difference between the distillation point numbers corresponding to the distillation point X and the distillation point Y in the whole distillation process; if
Figure 426955DEST_PATH_IMAGE006
,but
Figure 422593DEST_PATH_IMAGE005
for infinity.

在匹配距离公式中,将短期反应描述子相似度和外观特征混合编码相似度相乘后再被数值1相减,即

Figure 888209DEST_PATH_IMAGE019
表示了两种数据的整体差异,进一步引入馏程点序号的差异,使得后续匹配过程能够保证将不同馏程点序号的馏程点进行匹配的同时,还保证了匹配的两个馏程点的馏程点序号差异不会过大,增加后续预测数据的参考性。In the matching distance formula, the short-term response descriptor similarity is multiplied by the appearance feature mixed coding similarity and then subtracted by the value 1, that is
Figure 888209DEST_PATH_IMAGE019
Indicates the overall difference between the two data, and further introduces the difference in the number of the distillation point, so that the subsequent matching process can ensure that the distillation points with different numbers of the distillation point are matched, and at the same time ensure the matching of the two distillation points. The difference in the serial number of the distillation procedure point will not be too large, which increases the reference of subsequent prediction data.

在本发明实施例中,

Figure 209469DEST_PATH_IMAGE003
为两个短期反应描述子之间的余弦相似度,
Figure 26115DEST_PATH_IMAGE004
为两个外观特征混合编码之间的余弦相似度。需要说明的是,余弦相似度为本领域技术人员熟知的现有技术,在此以短期反应描述子之间的余弦相似度为例,列出其余弦相似度的表达式:In the embodiment of the present invention,
Figure 209469DEST_PATH_IMAGE003
is the cosine similarity between two short-term response descriptors,
Figure 26115DEST_PATH_IMAGE004
Cosine similarity between two appearance feature mixture encodings. It should be noted that the cosine similarity is an existing technology well known to those skilled in the art. Taking the cosine similarity between short-term response descriptors as an example, the expression of the cosine similarity is listed:

Figure 825444DEST_PATH_IMAGE020
Figure 825444DEST_PATH_IMAGE020

其中,

Figure 411146DEST_PATH_IMAGE021
为短期反应描述子
Figure 903307DEST_PATH_IMAGE022
中的第
Figure 941670DEST_PATH_IMAGE023
个元素,
Figure 544690DEST_PATH_IMAGE024
为短期反应描述子
Figure 250478DEST_PATH_IMAGE025
中的第
Figure 647961DEST_PATH_IMAGE023
个元素。in,
Figure 411146DEST_PATH_IMAGE021
short-term response descriptor
Figure 903307DEST_PATH_IMAGE022
in the first
Figure 941670DEST_PATH_IMAGE023
elements,
Figure 544690DEST_PATH_IMAGE024
short-term response descriptor
Figure 250478DEST_PATH_IMAGE025
in the first
Figure 647961DEST_PATH_IMAGE023
elements.

在上述的D公式的距离约束下,K-M算法会搜索得到若干个匹配对,且匹配对之间的匹配结果大部分是馏程节点之间相差1个,例如10%和20%之间;且短期反应模式中,温度特征以及V的参数是近似的,此时的外观是相似的。Under the distance constraint of the above D formula, the K-M algorithm will search for several matching pairs, and most of the matching results between the matching pairs are 1 difference between the distillation range nodes, for example, between 10% and 20%; and In the short-term response mode, the temperature characteristics and the parameters of V are approximate, and the appearance at this time is similar.

步骤S3:根据匹配对中两个馏程点之间的设定温度值差异、蒸汽温度差异和加热元件温度差异获得可切换指数;选择可切换指数最大的预设数量个匹配对中的数据作为初始训练集,根据初始训练集中互相匹配的馏程点数据对初始训练集进行增广,获得训练集,根据训练集训练长短期记忆生产组分配制预测模型。Step S3: Obtain the switchable index according to the set temperature value difference, the steam temperature difference and the heating element temperature difference between the two distillation range points in the matching pair; select the data in the preset number of matching pairs with the largest switchable index as For the initial training set, the initial training set is augmented according to the matching distillation point data in the initial training set to obtain a training set, and the long-short-term memory production group is trained according to the training set to prepare a prediction model.

一个匹配对中的馏程点能够代表处于不同或尽可能相似的参数下,能够匹配的下一个馏程的温度值。基于这些温度值,和各个馏程点的温度信息的相似程度,筛选出匹配对中较为合适的参考匹配对:The distillation range points in a matched pair can represent the temperature values of the next distillation range that can be matched under different or as similar as possible parameters. Based on these temperature values and the similarity of the temperature information of each distillation range point, a more suitable reference matching pair is selected among the matching pairs:

对每个匹配对定义可切换指数,该指数代表在上下馏程之间,可持续延迟蒸馏的程度,若这个指数越大,则认为在此情况下可产出的馏分是较多的,即根据匹配对中两个馏程点中的设定温度值差异、蒸汽温度差异和加热元件温度差异获得可切换指数,具体包括:Define a switchable index for each matching pair, which represents the degree of sustainable delayed distillation between the upper and lower distillation ranges. If the index is larger, it is considered that there are more fractions that can be produced in this case, that is A switchable index is obtained from the difference in set temperature values, vapor temperature difference and heating element temperature difference in two distillation range points in a matched pair, including:

根据可切换指数公式获得可切换指数,可切换指数公式包括:The switchable index is obtained according to the switchable index formula, which includes:

Figure 439200DEST_PATH_IMAGE007
Figure 439200DEST_PATH_IMAGE007

其中,

Figure 580331DEST_PATH_IMAGE008
为馏程点X和馏程点Y之间的可切换指数,
Figure 875046DEST_PATH_IMAGE009
为馏程点X和馏程点Y之间的设定温度值的差值绝对值,
Figure 443431DEST_PATH_IMAGE002
为馏程点X和馏程点Y之间的匹配距离,
Figure 721965DEST_PATH_IMAGE010
为馏程点X和馏程点Y之间蒸汽温度序列的动态时间规整距离,
Figure 666788DEST_PATH_IMAGE011
为馏程点X和馏程点Y之间加热元件温度序列的动态时间规整距离。in,
Figure 580331DEST_PATH_IMAGE008
is the switchable index between distillation point X and distillation point Y,
Figure 875046DEST_PATH_IMAGE009
is the absolute value of the difference between the set temperature values between the distillation point X and the distillation point Y,
Figure 443431DEST_PATH_IMAGE002
is the matching distance between distillation point X and distillation point Y,
Figure 721965DEST_PATH_IMAGE010
is the dynamic time-regularized distance of the steam temperature sequence between distillation point X and distillation point Y,
Figure 666788DEST_PATH_IMAGE011
is the dynamic time warping distance of the temperature sequence of the heating element between the distillation point X and the distillation point Y.

在可切换指数公式中,分母为两个距离的和,即距离越大说明两个馏程点的差异越大,则可切换指数越小;分子位设定温度值的差值绝对值和匹配距离乘积的倒数,即差值绝对值越小且匹配距离越小,说明两个馏程点越接近,可切换指数越大。In the switchable index formula, the denominator is the sum of the two distances, that is, the larger the distance, the greater the difference between the two distillation points, and the smaller the switchable index; the absolute value of the difference between the numerator and the set temperature value matches The reciprocal of the distance product, that is, the smaller the absolute value of the difference and the smaller the matching distance, it means that the closer the two distillation points are, the larger the switchable index is.

选择可切换指数最大的预设数量个匹配对中的数据作为初始训练集,在本发明实施例中,选择可切换指数排名前百分之二十五的匹配对构建初始训练集。The data in the preset number of matching pairs with the largest switchable index is selected as the initial training set. In the embodiment of the present invention, the top 25% of the matching pairs with the switchable index are selected to construct the initial training set.

对于排名前百分之二十五的匹配对中,其中任一匹配对,都有馏程的先后,以较先前的那个馏程点所对应的产后数据,将初始训练集进行增广,即根据初始训练集中互相匹配的馏程点数据对初始训练集进行增广,获得训练集,具体增广方法包括:For the top 25% of the matching pairs, any matching pair has a sequence of distillation ranges, and the initial training set is augmented with the postpartum data corresponding to the previous distillation point, that is According to the matching distillation point data in the initial training set, the initial training set is augmented to obtain the training set. The specific augmentation methods include:

将目标馏程点的体积空速、设定温度、蒸汽温度和加热元件温度替换为匹配馏程点所属的蒸馏过程中对应馏程点的体积空速、设定温度、蒸汽温度和加热元件温度,获得第一增广训练数据;将目标馏程点下一个馏程点的体积空速、设定温度、蒸汽温度和加热元件温度替换为匹配馏程点所属的蒸馏过程中对应馏程点的体积空速、设定温度、蒸汽温度和加热元件温度,获得第二增广训练数据;分别将目标馏程点和目标馏程点下一个馏程点的体积空速、设定温度、蒸汽温度和加热元件温度同时替换为匹配馏程点所属的蒸馏过程中对应馏程点的体积空速、设定温度、蒸汽温度和容器内加热元件温度,获得第三增广训练数据;将每个馏程点的第一增广训练数据、第二增广训练数据和第三增广训练数据并获得训练集。Replace the volume space velocity, set temperature, vapor temperature, and heating element temperature of the target distillation point with the volume space velocity, set temperature, vapor temperature, and heating element temperature of the corresponding distillation point in the distillation process to which the matching distillation point belongs , to obtain the first augmented training data; replace the volumetric space velocity, set temperature, steam temperature and heating element temperature of the distillation point next to the target distillation point with the corresponding distillation point in the distillation process to which the matching distillation point belongs The volume space velocity, set temperature, steam temperature and heating element temperature are obtained to obtain the second augmented training data; respectively, the volume space velocity, set temperature, and steam temperature of the target distillation point and the next distillation point and the temperature of the heating element are simultaneously replaced with the volume space velocity, set temperature, steam temperature and temperature of the heating element in the vessel corresponding to the distillation process to which the distillation point belongs to obtain the third augmented training data; The first augmented training data, the second augmented training data and the third augmented training data of the process point are obtained to obtain a training set.

例如,对于目标馏程点a而言,其存在一个所属的蒸馏过程A,且在蒸馏过程A中的馏程点序号为i;目标馏程点a对应一个匹配馏程点b,匹配馏程点b所属蒸馏过程B。则目标馏程点a的第一增广训练数据是由蒸馏过程B中馏程点序号为i的馏程点数据替换而来;第二增广训练数据是将蒸馏过程A中馏程点序号为i+1的馏程点数据替换为蒸馏过程B中馏程点序号为i+1的馏程点数据;第三增广训练数据是将蒸馏过程A中馏程点序号为i的馏程点数据替换为蒸馏过程B中馏程点序号为i的馏程点数据,将蒸馏过程A中馏程点序号为i+1的馏程点数据替换为蒸馏过程B中馏程点序号为i+1的馏程点数据。For example, for the target distillation point a, there is a distillation process A to which it belongs, and the number of the distillation point in the distillation process A is i; the target distillation point a corresponds to a matching distillation point b, and the matching distillation Point b belongs to distillation process B. Then the first augmented training data of the target distillation point a is replaced by the distillation process point number i in the distillation process B; the second augmented training data is the distillation process point number i in the distillation process A The distillation point data of i+1 is replaced by the distillation process point number i+1 in the distillation process B; the third augmented training data is the distillation process point number i in the distillation process A The point data is replaced by the distillation process point number i in distillation process B, and the distillation process point number i+1 in distillation process A is replaced by the distillation process point number i in distillation process B +1 for the distillation point data.

基于K-M的匹配作用,基于短期反应模式中,温度特征以及V的参数是近似的,外观是相似的产后数据为数据源。通过增广出在此情况下可产出的馏分是较多的控制模型训练数据。长短期记忆生产组分配制预测模型可以基于Batch Size的作用下,基于增广的若干条数据做误差平差,从而得到更佳切与产后数据相匹配的预测结果。对于训练集,基于前一个馏程的T_A、T_B,标签设定为下一个馏程的T_C和空速V。Based on the matching effect of K-M, based on the short-term response mode, the temperature characteristics and the parameters of V are approximate, and the appearance is similar to the postpartum data as the data source. The fraction that can be produced in this case by augmentation is more control model training data. The long-short-term memory production group allocation prediction model can be based on Batch Size, and based on several augmented data for error adjustment, so as to obtain better prediction results that match the postpartum data. For the training set, based on T_A, T_B of the previous distillation range, the label is set to T_C and space velocity V of the next distillation range.

步骤S4:将实时蒸馏过程中实时馏程点的实时蒸汽温度和实时加热元件温度输入长短期记忆生产组分配制预测模型中,获得下一个蒸馏过程中的设定温度值和体积空速,并根据下一个蒸馏过程中的设定温度值和体积空速对容器参数进行控制。Step S4: Input the real-time steam temperature and the real-time heating element temperature at the real-time distillation point in the real-time distillation process into the long-short-term memory production component allocation prediction model to obtain the set temperature value and volume space velocity in the next distillation process, and The vessel parameters are controlled according to the set temperature value and the volume space velocity in the next distillation process.

至此,在实时的生产过程中,基于一个实时蒸馏过程中实时馏程点的实时蒸汽温度和实时加热元件温度,做重采样后输入长短期记忆生产组分配制预测模型,即可获得下一个馏程点对应的设定温度值和体积空速,通过工作人员手动设置或者蒸馏容器的自动化设置,即可达到预测性的优化控制效果,提高馏程的馏分产出量。So far, in the real-time production process, based on the real-time steam temperature and real-time heating element temperature at the real-time distillation point in a real-time distillation process, after resampling, input the long-short-term memory production component allocation prediction model to obtain the next distillation The set temperature value and volume space velocity corresponding to the process point can be set manually by the staff or automatically set by the distillation vessel, so as to achieve the predictive optimization control effect and increase the distillate output of the distillation range.

综上所述,本发明实施例获取每个馏程点的反应特征,根据反应特征将不同控制参数下蒸馏过程中的馏程点进行匹配,获得多个匹配对。进一步根据匹配对内馏程点之间的反应特征差异获得可切换指数,基于可切换指数选取匹配对构建初始训练集,并根据匹配数据对初始训练集进行增广,以获得训练集用于训练长短期记忆生产组分配制预测模型,进而实现对实时数据的预测,根据预测参数对容器参数进行控制。本发明实施例通过构建组分配制模型预测未来馏程点所需调整的参数,保证了生产效率。To sum up, the embodiments of the present invention obtain the reaction characteristics of each distillation point, and match the distillation points in the distillation process under different control parameters according to the reaction characteristics to obtain multiple matching pairs. Further, the switchable index is obtained according to the reaction characteristic difference between the internal distillation points of the matching pair, and the matching pair is selected based on the switching index to construct the initial training set, and the initial training set is augmented according to the matching data to obtain the training set for training The long-short-term memory production group configures the prediction model, and then realizes the prediction of real-time data, and controls the container parameters according to the prediction parameters. In the embodiment of the present invention, the production efficiency is ensured by constructing a component allocation model to predict the parameters that need to be adjusted at future distillation range points.

需要说明的是:上述本发明实施例先后顺序仅仅为了描述,不代表实施例的优劣。在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中,多任务处理和并行处理也是可以的或者可能是有利的。It should be noted that: the order of the above embodiments of the present invention is only for description, and does not represent the advantages and disadvantages of the embodiments. The processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. Multitasking and parallel processing are also possible or may be advantageous in certain embodiments.

本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。Each embodiment in this specification is described in a progressive manner, the same and similar parts of each embodiment can be referred to each other, and each embodiment focuses on the differences from other embodiments.

以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included in the protection of the present invention. within range.

Claims (8)

1. A method for predictive control of a component dispensing model for chemical production, the method comprising:
collecting the set temperature value, the steam temperature and the temperature of a heating element in a container of each distillation range point in the distillation process; obtaining appearance characteristic mixed codes of reactants in the container; constructing a short-term reaction descriptor according to the steam temperature, the heating element temperature and the volume space velocity of the container in the sampling process;
matching different distillation range points according to appearance characteristic mixed coding difference and short-term reaction descriptor difference between the distillation range points in the distillation process under different control parameters to obtain a plurality of matched pairs;
obtaining a switchable index based on a set temperature value difference between two boiling point values in the matched pair, a steam temperature difference, and a heating element temperature difference; selecting data in a preset number of matching pairs with the largest switchable indexes as an initial training set, amplifying the initial training set according to mutually matched distillation range data in the initial training set to obtain a training set, and training a long-term and short-term memory production component system prediction model according to the training set;
inputting the real-time steam temperature of the real-time distillation range point and the real-time heating element temperature in the real-time distillation process into the long-short term memory production component preparation prediction model to obtain the set temperature value and the volume airspeed corresponding to the next distillation range point, and controlling the container parameters according to the set temperature value and the volume airspeed corresponding to the next distillation range point.
2. The method as claimed in claim 1, wherein the obtaining of the appearance feature mixture codes of the reagents in the containers comprises:
obtaining descriptive words of appearance characteristics of reactants in the container, and constructing TF-IDF appearance codes according to the descriptive words; collecting color information of the container reactant, and constructing an appearance color code according to component values in an RGB color space; and combining the TF-IDF appearance coding and the appearance color coding to obtain the appearance characteristic mixed coding.
3. The method as claimed in claim 2, wherein the obtaining of the appearance feature hybrid coding comprises:
expanding the dimension of the appearance color code to obtain a first expanded code; merging the first extended code and the TF-IDF appearance code to obtain a second extended code; and reducing the dimension of the second extended code to a preset dimension to obtain an appearance characteristic mixed code.
4. The method as claimed in claim 1, wherein the step of constructing the short-term response descriptor according to the steam temperature, the heating element temperature and the volume space velocity of the container during sampling comprises:
the short-term reaction descriptor is composed of a steam temperature sequence median value, a heating element temperature sequence median value, a steam temperature sequence mean value, a heating element temperature sequence mean value, a steam temperature sequence initial value and a heating element temperature sequence initial value at each distillation range point in the distillation process.
5. The method of claim 1, wherein the step of obtaining a plurality of matching pairs comprises:
constructing a matching distance according to appearance characteristic mixed coding difference and short-term reaction descriptor difference between different distillation processes under different reaction parameters; and matching different distillation processes by utilizing a KM matching algorithm according to the matching distance to obtain a plurality of matching pairs.
6. The method as claimed in claim 5, wherein the step of constructing the matching distance according to the appearance feature mixed coding difference and short-term response descriptor difference between different distillation processes under different response parameters comprises:
obtaining a matching distance according to a matching distance formula, wherein the matching distance formula comprises:
Figure 438329DEST_PATH_IMAGE001
wherein,
Figure 960446DEST_PATH_IMAGE002
is the matching distance between the boiling point X and the boiling point Y,
Figure 627051DEST_PATH_IMAGE003
short term reaction descriptor similarity between boiling point X and boiling point Y,
Figure 158395DEST_PATH_IMAGE004
the similarity is coded for the appearance feature mixture between the boiling point X and the boiling point Y,
Figure 328476DEST_PATH_IMAGE005
the absolute value of the difference of the corresponding distillation range point serial numbers of the distillation range point X and the distillation range point Y in the whole distillation process is shown; if it is
Figure 21495DEST_PATH_IMAGE006
Then, then
Figure 175396DEST_PATH_IMAGE005
Is infinite.
7. The method of claim 1, wherein the obtaining the switchable index according to the set temperature value difference, the steam temperature difference and the heating element temperature difference in the two distillation processes in the matching pair comprises:
obtaining a switchable index according to a switchable index formula, the switchable index formula comprising:
Figure 792322DEST_PATH_IMAGE007
wherein,
Figure 72036DEST_PATH_IMAGE008
is a switchable index between the boiling point X and the boiling point Y,
Figure 749005DEST_PATH_IMAGE009
is the absolute value of the difference between the set temperature values between the distillation range point X and the distillation range point Y,
Figure 373891DEST_PATH_IMAGE002
is the matching distance between boiling point X and boiling point Y,
Figure 466612DEST_PATH_IMAGE010
is the dynamic time-warping distance of the steam temperature sequence between the distillation range point X and the distillation range point Y,
Figure 126132DEST_PATH_IMAGE011
is the dynamic time-warping distance of the heating element temperature sequence between boiling point X and boiling point Y.
8. The method of claim 1, wherein the augmenting the initial training set based on distillation process data that match each other in the initial training set comprises:
replacing the volume airspeed, the set temperature, the steam temperature and the heating element temperature of the target distillation range point with the volume airspeed, the set temperature, the steam temperature and the heating element temperature of the corresponding distillation range point in the distillation process to which the matched distillation range point belongs to obtain first augmentation training data;
replacing the volume airspeed, the set temperature, the steam temperature and the heating element temperature of the next distillation range point of the target distillation range point with the volume airspeed, the set temperature, the steam temperature and the heating element temperature of the corresponding distillation range point in the distillation process to which the matched distillation range point belongs to obtain second augmentation training data;
respectively replacing the volume airspeed, the set temperature, the steam temperature and the heating element temperature of the target distillation range point and the next distillation range point of the target distillation range point with the volume airspeed, the set temperature, the steam temperature and the heating element temperature of the corresponding distillation range point in the distillation process to which the matched distillation range point belongs, and obtaining third augmentation training data;
and combining the first augmented training data, the second augmented training data and the third augmented training data of each distillation range point to obtain a training set.
CN202211408678.4A 2022-11-11 2022-11-11 Model Predictive Control Method of Component Formulation for Chemical Production Pending CN115454009A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211408678.4A CN115454009A (en) 2022-11-11 2022-11-11 Model Predictive Control Method of Component Formulation for Chemical Production

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211408678.4A CN115454009A (en) 2022-11-11 2022-11-11 Model Predictive Control Method of Component Formulation for Chemical Production

Publications (1)

Publication Number Publication Date
CN115454009A true CN115454009A (en) 2022-12-09

Family

ID=84295606

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211408678.4A Pending CN115454009A (en) 2022-11-11 2022-11-11 Model Predictive Control Method of Component Formulation for Chemical Production

Country Status (1)

Country Link
CN (1) CN115454009A (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102393636A (en) * 2011-09-22 2012-03-28 任季明 Control system and method for petroleum refining process
CN104152179A (en) * 2014-09-01 2014-11-19 江苏联东化工股份有限公司 Production method of high-boiling-point aromatic solvent oil
CN105844667A (en) * 2016-03-25 2016-08-10 中国矿业大学 Structural target tracking method of compact color coding
CN109726844A (en) * 2017-10-31 2019-05-07 中国石油化工股份有限公司 Product yield method for automatically regulating, system and the storage equipment of hydrocracking unit
CN109726845A (en) * 2017-10-31 2019-05-07 中国石油化工股份有限公司 Product yield prediction technique, method for establishing model and the storage equipment being hydrocracked
CN109754113A (en) * 2018-11-29 2019-05-14 南京邮电大学 Load forecasting method based on dynamic time warping and long and short-term memory
CN110102075A (en) * 2019-03-18 2019-08-09 山东金特昂莱测控技术有限公司 The chlorination control method and device that view-based access control model identifies in bromine distillation technique
CN113764046A (en) * 2021-09-24 2021-12-07 中国石油化工股份有限公司 Method for improving yield of high value-added product in catalytic diesel oil hydroconversion

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102393636A (en) * 2011-09-22 2012-03-28 任季明 Control system and method for petroleum refining process
CN104152179A (en) * 2014-09-01 2014-11-19 江苏联东化工股份有限公司 Production method of high-boiling-point aromatic solvent oil
CN105844667A (en) * 2016-03-25 2016-08-10 中国矿业大学 Structural target tracking method of compact color coding
CN109726844A (en) * 2017-10-31 2019-05-07 中国石油化工股份有限公司 Product yield method for automatically regulating, system and the storage equipment of hydrocracking unit
CN109726845A (en) * 2017-10-31 2019-05-07 中国石油化工股份有限公司 Product yield prediction technique, method for establishing model and the storage equipment being hydrocracked
CN109754113A (en) * 2018-11-29 2019-05-14 南京邮电大学 Load forecasting method based on dynamic time warping and long and short-term memory
CN110102075A (en) * 2019-03-18 2019-08-09 山东金特昂莱测控技术有限公司 The chlorination control method and device that view-based access control model identifies in bromine distillation technique
CN113764046A (en) * 2021-09-24 2021-12-07 中国石油化工股份有限公司 Method for improving yield of high value-added product in catalytic diesel oil hydroconversion

Similar Documents

Publication Publication Date Title
CN109814513B (en) Catalytic cracking unit optimization method based on data model
Wang et al. Cross-modal food retrieval: learning a joint embedding of food images and recipes with semantic consistency and attention mechanism
Diallo et al. Auto-attention mechanism for multi-view deep embedding clustering
Yan et al. A unified pretraining framework for passage ranking and expansion
Aimin et al. Neighborhood preserving regression embedding based data regression and its applications on soft sensor modeling
Li et al. Intelligence customs declaration for cross-border e-commerce based on the multi-modal model and the optimal window mechanism
De Oliveira et al. Binary classification on French hospital data: benchmark of 7 machine learning algorithms
CN115454009A (en) Model Predictive Control Method of Component Formulation for Chemical Production
Teng et al. Machine-learned digital phase switch for sustainable chemical production
Yang et al. Segmented sequence decomposition-Informer model for deformation of arch dams
Wang et al. Synthesis and optimization of low-temperature gas separation processes
Liu et al. Hyperspectral band selection based on consistency-measure of neighborhood rough set theory
Liu et al. A learnable front-end based efficient channel attention network for heart sound classification
Han et al. Novel CNN-based transformer integrating Boruta algorithm for production prediction modeling and energy saving of industrial processes
Si et al. Prediction of brown tide algae using improved Gramian angular field and deep learning based on laser-induced fluorescence spectrum
Hu et al. Histopathology cross-modal retrieval based on dual-transformer network
Han et al. Linear optimization fusion model based on fuzzy C-means: Case study of energy efficiency evaluation in ethylene product plants
Yan et al. Application of DS evidence theory in apple’s internal quality classification
Wang et al. Industrial units modeling using self-attention network based on feature selection and pattern classification
Liu et al. Differentiable neural architecture search for domain adaptation in fault diagnosis
Krier et al. Supervised variable clustering for classification of NIR spectra.
Kumar et al. Machine learning solutions for investigating streams data using distributed frameworks: Literature review
Yuan et al. Stacked attention-based autoencoder with feature fusion and its application for quality prediction
CN109408527B (en) RDF structured query automatic construction method based on vector space
Wang et al. A Novel Reservoir Modeling Method based on Improved Hierarchical XGBoost

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