A real-time optimization method and system for operating parameters of an ethyl acetate reactive rectification column

Abstract

The application discloses an ethyl acetate reactive rectification tower operation parameter real-time optimization method and system. The method comprises the following steps: in each real-time optimization period, current data and historical data sets are obtained; the current operation data is substituted into an ethyl acetate reactive rectification mechanism model to calculate a first theoretical tower top ethyl acetate concentration; different operation offsets are applied to calculate a second theoretical tower top ethyl acetate concentration; after the correction of each second theoretical tower top ethyl acetate concentration, an operation-concentration data set is composed of different operation offsets, the current operation data and a third theoretical tower top ethyl acetate concentration corresponding to each operation offset, and the operation-concentration data set is filtered according to the calculation result of the cosine similarity; the operation optimality of each instance in the filtered operation-concentration data set is calculated, the instance with the highest operation optimality is taken as an optimized operation parameter result, and the operation parameter real-time optimization of the ethyl acetate reactive rectification tower is completed.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of industrial control, and particularly relates to a method and system for real-time optimization of operating parameters of an ethyl acetate reaction rectification tower. BACKGROUND

[0002] In the process industry, the production of ethyl acetate is mainly completed through esterification of acetic acid and ethanol with concentrated sulfuric acid as a catalyst. However, the esterification reaction is a reversible reaction, and thus there is a problem of coexistence of ethyl acetate product and acetic acid and ethanol raw materials. In order to make the reaction proceed fully and improve the conversion rate, a reaction rectification tower is generally used to realize continuous generation of the above process, that is, the chemical reaction and rectification are coupled, and the product is separated in the rectification tower while acetic acid and ethanol are reacted. This method is beneficial to the reversible reaction to proceed in a direction favorable to the reaction product, thereby improving the conversion rate.

[0003] According to the normal pressure boiling point data of acetic acid, ethanol and ethyl acetate, the constant boiling point of water-ethanol-ethyl acetate is the lowest, and ethanol is a raw material and is not expected to enter the product. Therefore, in the reaction, an excess of acetic acid is generally added to make ethanol be consumed fully. This can largely ensure the purity of ethyl acetate, and to further improve the product purity, optimization needs to be performed from the aspects of feed flow, acetic acid-ethanol ratio, reflux ratio, tower kettle temperature and other operating parameters, so as to maintain a high product purity of ethyl acetate.

[0004] For an operator, optimization of the above operating parameters requires a large amount of practical operation experience, and novice operators often lack such ability. At the same time, if a theoretical analysis is performed on the ethyl acetate reaction rectification tower, it is difficult to accurately describe the physical and chemical change processes occurring therein, and an accurate model for optimizing the operating parameters cannot be provided. Using the running data of the reaction rectification tower and combining the running mechanism of the reaction rectification tower to establish an operating parameter optimization model can overcome the difficulties faced by the above theoretical analysis method and can be closer to the actual situation of the device.

[0005] The existing operating parameter optimization technology generally establishes a strict ethyl acetate reaction rectification model, and calculates corresponding operating parameters according to actual parameters of a production device, such as the height of the rectification tower, the number of trays, the designed feed amount, etc. The above method is generally used in the design stage, that is, the operating parameters calculated by the technology are told to the operator, and the operator adjusts on this basis.

[0006] The mechanism model-based operating parameter optimization method used in the above technology has bottlenecks in actual application, that is, it can only give inaccurate operating parameter results, mainly for prompting. Moreover, it cannot be dynamically adjusted according to the running data, and cannot achieve the effect of real-time optimization of operating parameters. SUMMARY

[0007] The present application aims to solve the problems in the prior art and provide an ethyl acetate reactive rectification column operation parameter real-time optimization method and system.

[0008] In order to achieve the above-mentioned application purposes, the present application specifically adopts the following technical solutions:

[0009] In the first aspect, the present application provides an ethyl acetate reactive rectification column operation parameter real-time optimization method, which comprises the following steps:

[0010] S1. In each real-time optimization cycle, the current data is obtained, if the historical data set cannot be obtained, the current data is added to the historical data set, and after a predetermined period of time, the current data and the historical data set are reacquired, and S2 is entered; if the historical data set can be obtained, S2 is directly entered; wherein the current data includes current operation data and current column top ethyl acetate concentration data, the historical data set contains historical operation data and historical column top ethyl acetate concentration data, and the operation data includes feed flow, ethanol / acetic acid ratio, reflux ratio and column bottom temperature;

[0011] S2. The current operation data is substituted into the ethyl acetate reactive rectification mechanism model to calculate the first theoretical column top ethyl acetate concentration;

[0012] S3. After different operation offsets are respectively applied to the current operation data, each operation data added with an operation offset is substituted into the ethyl acetate reactive rectification mechanism model to calculate the second theoretical column top ethyl acetate concentration corresponding to each operation offset, the first theoretical column top ethyl acetate concentration is corrected according to the deviation of the current column top ethyl acetate concentration data, the third theoretical column top ethyl acetate concentration corresponding to each operation offset is obtained, and the operation-concentration data set is composed of different operation offsets, the current operation data and the third theoretical column top ethyl acetate concentration corresponding to each operation offset;

[0013] S4. The column top ethyl acetate concentration data in the historical data set is sorted from high to low, the concentration rising optimization direction is calculated from the sorted column top ethyl acetate concentration data, and the operation-concentration data set in S3 is filtered according to the calculation result of the cosine similarity, to obtain the filtered operation-concentration data set;

[0014] S5. The operation optimality of each instance in the filtered operation-concentration data set is calculated, the instance with the highest operation optimality is taken as the optimized operation parameter result, and the operation parameter real-time optimization of the ethyl acetate reactive rectification column is completed.

[0015] On the basis of the above-mentioned scheme, each step can be realized in the following preferred specific manner.

[0016] As a preferred embodiment of the first aspect, in step S2, the function form of the ethyl acetate reactive distillation mechanism model is:

[0017] c = h(f, r f , t)

[0018] wherein c represents the theoretical top ethyl acetate concentration; f represents the feed flow rate; r f represents the ethanol / acetic acid ratio; r represents the reflux ratio; and t represents the column bottom temperature.

[0019] As a preferred embodiment of the first aspect, in step S3, the method for generating the operation offset is: calculating the minimum positive integer p that satisfies the following condition:

[0020]

[0021] The operation offset of the jth group is represented as follows:

[0022]

[0023] wherein X represents the base number, and the value range is X ∈ {2, 3, 5, 7}; X i represents the i-th power of the base number X; N represents the preset number of groups of operation offsets; represents the inverted decimal number after representing j in binary; represents the inverted decimal number after representing j in ternary; represents the inverted decimal number after representing j in quinary; represents the inverted decimal number after representing j in septenary.

[0024] As a preferred embodiment of the first aspect, in step S3, the formula for correcting the second theoretical top ethyl acetate concentration is:

[0025] c 3 = c 2 + c 0 - c 1

[0026] wherein c 0 represents the current top ethyl acetate concentration data; c 1 represents the first theoretical top ethyl acetate concentration; c 2 represents the second theoretical top ethyl acetate concentration; and c 3 represents the third theoretical top ethyl acetate concentration.

[0027] As a preferred embodiment of the first aspect, in step S4, the specific process of filtering the operation-concentration dataset is as follows: obtaining the operation offset in the historical dataset from the concentration increase optimization direction, calculating the cosine similarity of the operation offset of each instance in the operation-concentration dataset and the operation offset in the historical dataset, and deleting the instance with a cosine similarity less than 0 from the operation-concentration dataset, thereby obtaining the filtered operation-concentration dataset.

[0028] As a preferred embodiment of the first aspect, in step S5, the operation optimality OpOp is calculated as follows:

[0029]

[0030] wherein F represents the operation offset of the feed flow; F0 represents the operation offset of the feed flow in each instance in the filtered operation-concentration dataset relative to the current feed flow; R F represents the operation offset of the ethanol / acetic acid ratio; R F0 represents the operation offset of the ethanol / acetic acid ratio in each instance in the filtered operation-concentration dataset relative to the current feed flow; R represents the operation offset of the reflux ratio; R0 represents the operation offset of the reflux ratio in each instance in the filtered operation-concentration dataset relative to the current feed flow; T represents the operation offset of the column bottom temperature; T0 represents the operation offset of the column bottom temperature in each instance in the filtered operation-concentration dataset relative to the current feed flow; F 2 ,R F 2 ,R 2 ,T 2 , respectively represent the square of F, R F ,R,T,F0,R F0 ,R0,T0; C represents the ethyl acetate concentration offset.

[0031] In a second aspect, the present application provides an ethyl acetate reactive rectification column operation parameter real-time optimization system, which comprises:

[0032] a data acquisition module, configured to acquire current data in each real-time optimization cycle, if the historical dataset cannot be acquired, the current data is added to the historical dataset, and after a preset period of time, the current data and the historical dataset are reacquired, and S2 is entered; if the historical dataset can be acquired, S2 is directly entered; wherein the current data comprises current operation data and current column top ethyl acetate concentration data, the historical dataset contains historical operation data and historical column top ethyl acetate concentration data, and the operation data comprises feed flow, ethanol / acetic acid ratio, reflux ratio, and column bottom temperature;

[0033] A calculation module is configured to substitute the current operation data into an ethyl acetate reactive distillation mechanism model to calculate a first theoretical column top ethyl acetate concentration;

[0034] A data optimization module is configured to substitute the current operation data into the ethyl acetate reactive distillation mechanism model after different operation offsets are respectively applied to the current operation data to calculate a second theoretical column top ethyl acetate concentration corresponding to each operation offset, correct each second theoretical column top ethyl acetate concentration according to a deviation between the first theoretical column top ethyl acetate concentration and the current column top ethyl acetate concentration data to obtain a third theoretical column top ethyl acetate concentration corresponding to each operation offset, and form an operation-concentration data set composed of the different operation offsets, the current operation data and the third theoretical column top ethyl acetate concentration corresponding to each operation offset.

[0035] A filtering module is configured to sort the column top ethyl acetate concentration data in the historical data set from high to low, calculate a concentration rising optimization direction from the sorted column top ethyl acetate concentration data, and filter the operation-concentration data set in S3 according to a calculation result of the cosine similarity to obtain a filtered operation-concentration data set.

[0036] A result acquisition module is configured to calculate an operation optimality of each instance in the filtered operation-concentration data set, take an instance with the highest operation optimality as an optimized operation parameter result, and complete real-time optimization of the operation parameters of the ethyl acetate reactive distillation column.

[0037] Compared with the prior art, the present application has the following beneficial effects:

[0038] The present application solves the technical problem of giving an operation parameter combination capable of gradually reaching an optimal product purity by analyzing the operation mechanism of the reactive distillation column and combining the operation data of the ethyl acetate reactive distillation.

[0039] The key point of the present application is how to use the ethyl acetate operation data and combine the mechanism model of the ethyl acetate reactive distillation to realize optimization of the operation parameters. The present application has the beneficial effect of being able to mine operation parameter optimization knowledge from the operation data of the device, thereby reducing the requirement for accuracy. BRIEF DESCRIPTION OF DRAWINGS

[0040] Figure 1 The step flowchart of the present application;

[0041] Figure 2 The process flow diagram of the ethyl acetate reactive distillation column in the embodiment of the present application;

[0042] Figure 3 The system block diagram of the present application. DETAILED DESCRIPTION

[0043] In order to make the above objectives, characteristics and advantages of the present application more obvious and easy to understand, the specific embodiments of the present application are described in detail below with reference to the drawings. In the following description, a large number of specific details are set forth in order to facilitate a thorough understanding of the present application. However, the present application can be implemented in many other ways different from those described herein, and those skilled in the art can make similar improvements without departing from the spirit of the present application, so the present application is not limited to the specific embodiments disclosed below. The technical features in each embodiment of the present application can be combined accordingly without conflict.

[0044] In the description of the present application, it should be understood that the terms "first", "second" are only used for distinguishing description purposes and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the technical features indicated. Therefore, the features defined with "first", "second" can explicitly or implicitly include at least one of the features.

[0045] The ethyl acetate reactive rectification column process involved in the present embodiment is shown in Figure 2 , the feed is a mixture of ethanol / acetic acid in a certain ratio, and the overhead product is ethyl acetate product.

[0046] As shown in Figure 1 , in a preferred implementation of the present application, the above-mentioned real-time optimization method of ethyl acetate reactive rectification column operating parameters includes the following S1-S5 steps. The specific implementation process is described below.

[0047] S1. In each real-time optimization cycle, the current data is obtained, if the historical data set cannot be obtained, the current data is added to the historical data set, and after a predetermined period of time, the current data and the historical data set are reacquired, and S2 is entered; if the historical data set can be obtained, S2 is directly entered; wherein the current data includes the current operating data and the current overhead ethyl acetate concentration data c, the historical data set contains historical operating data and historical overhead ethyl acetate concentration data c, and the operating data includes feed flow rate f, ethanol / acetic acid ratio r f , reflux ratio r and column bottom temperature t.

[0048] It should be noted that in step S1 of the present application, the above-mentioned predetermined period is generally 30 to 60 minutes. In step S1 of the present embodiment, the predetermined period is set to 30 minutes.

[0049] In the present embodiment, the current operating data is first obtained, specifically the feed flow rate f = 100 kg / h, the ethanol / acetic acid ratio r f= 0.9, reflux ratio r = 4 and column bottom temperature t = 80 degrees Celsius, overhead ethyl acetate concentration c = 78%. The historical data set size is 1, i.e. one set of historical data is stored, specifically feed flow rate 110 kg / h, ethanol / acetic acid ratio 0.93, reflux ratio 4.2, column bottom temperature 79 degrees Celsius, overhead ethyl acetate concentration 77%.

[0050] S2. Substitute the current operating data into the ethyl acetate reactive distillation mechanism model to calculate a first theoretical overhead ethyl acetate concentration.

[0051] It should be noted that in step S2 of the present application, the implementation method of the ethyl acetate reactive distillation mechanism model belongs to the prior art, and its functional form is:

[0052] c = h(f, r f , t)

[0053] wherein c represents the first theoretical overhead ethyl acetate concentration, dimensionless; f represents the feed flow rate, unit kg / h; r f represents the ethanol / acetic acid ratio, dimensionless; r represents the reflux ratio, dimensionless; and t represents the column bottom temperature, unit degrees Celsius.

[0054] In this embodiment, the current operating data of feed flow rate f = 100 kg / h, ethanol / acetic acid ratio r f = 0.9, reflux ratio r = 4 and column bottom temperature t = 80 degrees Celsius are substituted into the ethyl acetate reactive distillation mechanism model to calculate the first theoretical overhead ethyl acetate concentration as 78.5%.

[0055] S3. After applying different operating offsets to the current operating data respectively, substitute each operating data with the added operating offset into the ethyl acetate reactive distillation mechanism model to calculate a second theoretical overhead ethyl acetate concentration corresponding to each operating offset, correct each second theoretical overhead ethyl acetate concentration according to the deviation of the first theoretical overhead ethyl acetate concentration from the current overhead ethyl acetate concentration data to obtain a third theoretical overhead ethyl acetate concentration corresponding to each operating offset, and form an operation-concentration data set composed of different operating offsets, the current operating data and the third theoretical overhead ethyl acetate concentration corresponding to each operating offset.

[0056] It should be noted that in step S3 of the present application, the generation method of the above-mentioned operating offset is as follows: calculate the smallest positive integer p that satisfies the following condition:

[0057]

[0058] The operating offset of the jth group is represented as follows:

[0059]

[0060] wherein X represents a base, and the value range of X is X∈{2, 3, 5, 7}; X i represents the i-th power of the base X; N represents the preset group number of operation offset; represents the inverted decimal of j represented by 2; represents the inverted decimal of j represented by 3; represents the inverted decimal of j represented by 5; represents the inverted decimal of j represented by 7.

[0061] In the present embodiment, different operation offsets are applied to the current operation data of feed flow 100 kg / h, ethanol / acetic acid ratio 0.9, reflux ratio 4, and column bottom temperature 80 degrees Celsius, and in the present embodiment, 100 groups of operation offset combinations are generated, and the minimum positive integer p=6 satisfying the above conditions is calculated. Taking the second group of operation offsets as an example, the calculation result is as follows:

[0062]

[0063] Here, taking the in the second group of operation offsets as an example, the calculation method of the result is described. First, the group number index j=2 is represented by 10 using 2, and then the result is represented as the inverted decimal 0.01, and the inverted decimal is expanded using 2 to obtain

[0064] It should be noted that in step S3 of the present application, the formula for correcting the second theoretical column top ethyl acetate concentration is as follows:

[0065] c 3 =c 2 +c 0 -c 1

[0066] wherein c 0 represents the current column top ethyl acetate concentration data; c 1 represents the first theoretical column top ethyl acetate concentration; c 2 represents the second theoretical column top ethyl acetate concentration; and c 3 represents the third theoretical column top ethyl acetate concentration.

[0067] In the embodiment, the N=100 sets of operation offset are substituted into the ethyl acetate reactive rectification mechanism model together with the current operation data to obtain N=100 sets of corresponding second theoretical column top ethyl acetate concentration, then the second theoretical column top ethyl acetate concentration is corrected according to the deviation of the first theoretical column top ethyl acetate concentration obtained in S2 from the current column top ethyl acetate concentration data to obtain third theoretical column top ethyl acetate concentration. Finally, the 100 sets of operation offset, the current operation data and the corresponding 100 sets of third theoretical column top ethyl acetate concentration constitute the operation-concentration data set.

[0068] S4. The column top ethyl acetate concentration data in the historical data set is sorted in descending order, the concentration increase optimization direction is calculated from the sorted column top ethyl acetate concentration data, and the operation-concentration data set in S3 is filtered according to the calculation result of the cosine similarity to obtain the filtered operation-concentration data set.

[0069] It should be noted that in step S4 of the present application, the specific process of filtering the operation-concentration data set is: obtaining the operation offset in the historical data set from the concentration increase optimization direction, calculating the cosine similarity of the operation offset of each instance in the operation-concentration data set and the operation offset in the historical data set, and deleting the instances with cosine similarity less than 0 from the operation-concentration data set to finally obtain the filtered operation-concentration data set.

[0070] In the embodiment, the cosine similarity of 100 sets of operation offset and the historical operation offset [-10, -0.3, 0.2, 1] needs to be calculated, and according to this, 58 sets of operation offset are deleted, that is, the size of the filtered operation-concentration data set is 42.

[0071] S5. The operation optimality of each instance in the filtered operation-concentration data set is calculated, the instance with the highest operation optimality is taken as the optimized operation parameter result, and the real-time optimization of the operation parameters of the ethyl acetate reactive rectification column is completed.

[0072] It should be noted that in step S5 of the present application, the operation optimality is calculated as follows:

[0073]

[0074] Wherein, F represents the operation offset of the feed flow; F0 represents the operation offset of the feed flow in each instance in the filtered operation-concentration data set relative to the current feed flow; R F represents the operation offset of the ethanol / acetic acid ratio; R F0represents the operation offset of the ethanol / acetic acid ratio in each instance in the filtered operation-concentration data set relative to the current feed flow rate; R represents the operation offset of the reflux ratio; R0 represents the operation offset of the reflux ratio in each instance in the filtered operation-concentration data set relative to the current feed flow rate; T represents the operation offset of the column bottom temperature; T0 represents the operation offset of the column bottom temperature in each instance in the filtered operation-concentration data set relative to the current feed flow rate; F 2 , R 2 ,T 2 , represent the square of F, R F , R, T, F0, R F0 , R0, T0, respectively; C represents the ethyl acetate concentration offset.

[0075] In the present embodiment, the above calculation formula of the operation optimality can select operation offsets with a wider variation direction. The instance with the highest operation optimality can be selected as the result of real-time optimization, i.e., new operation parameters (feed flow rate, ethanol / acetic acid ratio, reflux ratio, column bottom temperature) are obtained.

[0076] It should be further noted that the above-mentioned method for real-time optimization of the operation parameters of the ethyl acetate reactive rectification column in the embodiments can be essentially executed by a computer program or module. Therefore, based on the same inventive concept, another preferred embodiment of the present application also provides an ethyl acetate reactive rectification column operation parameter real-time optimization system corresponding to the above-mentioned method for real-time optimization of the operation parameters of the ethyl acetate reactive rectification column, as shown in Figure 3 which comprises:

[0077] a data acquisition module, configured to acquire current data in each real-time optimization period, if the historical data set cannot be acquired, the current data is added to the historical data set, and after a preset period, the current data and the historical data set are reacquired, and S2 is entered; if the historical data set can be acquired, S2 is directly entered; wherein the current data includes current operation data and current column top ethyl acetate concentration data, the historical data set contains historical operation data and historical column top ethyl acetate concentration data, and the operation data includes feed flow rate, ethanol / acetic acid ratio, reflux ratio, and column bottom temperature;

[0078] a calculation module, configured to substitute the current operation data into an ethyl acetate reactive rectification mechanism model to calculate a first theoretical column top ethyl acetate concentration;

[0079] The data optimization module is configured to: substitute each operation data after adding a different operation offset into the ethyl acetate reaction rectification mechanism model to calculate a second theoretical tower top ethyl acetate concentration corresponding to each operation offset; correct each second theoretical tower top ethyl acetate concentration according to a deviation between the first theoretical tower top ethyl acetate concentration and the current tower top ethyl acetate concentration data to obtain a third theoretical tower top ethyl acetate concentration corresponding to each operation offset; and generate an operation-concentration data set composed of the different operation offsets, the current operation data, and the third theoretical tower top ethyl acetate concentration corresponding to each operation offset.

[0080] The filtering module is configured to: sort the tower top ethyl acetate concentration data in the historical data set in descending order; calculate a concentration increase optimization direction from the sorted tower top ethyl acetate concentration data; and filter the operation-concentration data set in S3 according to the calculation result of the cosine similarity to obtain a filtered operation-concentration data set.

[0081] The result obtaining module is configured to: calculate an operation optimality of each instance in the filtered operation-concentration data set; and take an instance with the highest operation optimality as an optimized operation parameter result to complete real-time optimization of the operation parameters of the ethyl acetate reaction rectification tower.

[0082] The above-described embodiments are only a preferred scheme of the present application, and are not intended to limit the present application. Those skilled in the related art can make various changes and modifications without departing from the spirit and scope of the present application. Therefore, any technical scheme obtained by equivalent replacement or equivalent transformation falls within the protection scope of the present application.

Claims

1. A method for real-time optimization of operating parameters of an ethyl acetate reactive distillation column, characterized in that, Includes the following steps: S1. In each real-time optimization cycle, acquire the current data. If the historical dataset cannot be acquired, add the current data to the historical dataset. After a preset cycle, acquire the current data and the historical dataset again, and proceed to S2. If the historical dataset can be acquired, proceed directly to S2. The current data includes the current operating data and the current ethyl acetate concentration at the top of the column. The historical dataset contains historical operating data and historical ethyl acetate concentration at the top of the column. The operating data includes feed flow rate, ethanol / acetic acid ratio, reflux ratio, and column bottom temperature. S2. Substitute the current operating data into the ethyl acetate reactive distillation mechanism model to calculate the ethyl acetate concentration at the top of the first theoretical column; S3. After applying different operation offsets to the current operation data, substitute each operation data after adding the operation offset into the ethyl acetate reactive distillation mechanism model, calculate the second theoretical column top ethyl acetate concentration corresponding to each operation offset, and correct each second theoretical column top ethyl acetate concentration according to the deviation between the first theoretical column top ethyl acetate concentration and the current column top ethyl acetate concentration data to obtain the third theoretical column top ethyl acetate concentration corresponding to each operation offset. The operation-concentration dataset is composed of different operation offsets, the current operation data, and the third theoretical column top ethyl acetate concentration corresponding to each operation offset. S4. Sort the ethyl acetate concentration data at the top of the column in the historical dataset from high to low, calculate the optimal direction for concentration increase from the sorted ethyl acetate concentration data at the top of the column, and filter the operation-concentration dataset in S3 according to the calculation result of cosine similarity to obtain the filtered operation-concentration dataset. S5. Calculate the operational optimality of each instance in the filtered operational-concentration dataset, and take the instance with the highest operational optimality as the optimized operational parameter result to complete the real-time optimization of the operational parameters of the ethyl acetate reactive distillation column.

2. The method for real-time optimization of operating parameters of an ethyl acetate reactive distillation column as described in claim 1, characterized in that, In step S2, the functional form of the reactive distillation mechanism model for ethyl acetate is: c=h(f,r f ,r,t) Where c represents the theoretical concentration of ethyl acetate at the top of the column; f represents the feed flow rate; r f represents the ethanol / acetic acid ratio; r represents the reflux ratio; t represents the reboiler temperature.

3. The method for real-time optimization of operating parameters of an ethyl acetate reactive distillation column as described in claim 1, characterized in that, In step S3, the method for generating the operation offset is as follows: calculate the smallest positive integer p that satisfies the following conditions: The operation offset for the j-th group is represented as follows: Where X represents the cardinality, and its value ranges from X∈{2, 3, 5, 7}; X i Represents the base X raised to the power of i; N represents the number of preset operation offset groups; This represents the inverted decimal of j when represented in binary. This represents the inverted decimal of j when represented in base 3; This represents the inverted decimal of j when represented in base 5; This represents the inverted decimal of j when represented in base 7.

4. The method for real-time optimization of operating parameters of an ethyl acetate reactive distillation column as described in claim 1, characterized in that, In step S3, the formula for correcting the concentration of ethyl acetate at the top of the second theoretical column is as follows: c 3 =c 2 +c 0 -c 1 Among them, c 0 This indicates the current concentration of ethyl acetate at the top of the column; c 1 c represents the concentration of ethyl acetate at the top of the first theoretical column. 2 c represents the concentration of ethyl acetate at the top of the second theoretical column. 3 This indicates the concentration of ethyl acetate at the top of the third theoretical column.

5. The method for real-time optimization of operating parameters of an ethyl acetate reactive distillation column as described in claim 1, characterized in that, In step S4, the specific process of filtering the operation-concentration dataset is as follows: obtain the operation offset in the historical dataset from the direction of concentration increase optimization, calculate the cosine similarity between the operation offset of each instance in the operation-concentration dataset and the operation offset in the historical dataset, and delete instances with a cosine similarity less than 0 from the operation-concentration dataset, and finally obtain the filtered operation-concentration dataset.

6. The method for real-time optimization of operating parameters of an ethyl acetate reactive distillation column as described in claim 1, characterized in that, In step S5, the operation optimality OpOp is calculated as follows: Where F represents the operational offset of the feed flow rate; F0 represents the operational offset of the feed flow rate relative to the current feed flow rate in each instance of the filtered operational-concentration dataset; R F R represents the operational offset of the ethanol / acetic acid ratio. F0 R represents the operational offset of the ethanol / acetic acid ratio relative to the current feed flow rate for each instance in the filtered operational-concentration dataset; R represents the operational offset of the reflux ratio; R0 represents the operational offset of the reflux ratio relative to the current feed flow rate for each instance in the filtered operational-concentration dataset; T represents the operational offset of the reboiler temperature; T0 represents the operational offset of the reboiler temperature relative to the current feed flow rate for each instance in the filtered operational-concentration dataset; F 2 , R 2 ,T 2 , F and R respectively F ,R,T,F0,R F0 R0,T0 squared; C represents the ethyl acetate concentration offset.

7. A real-time optimization system for operating parameters of an ethyl acetate reactive distillation column, characterized in that, include: The data acquisition module is used to acquire the current data in each real-time optimization cycle. If the historical dataset cannot be acquired, the current data is added to the historical dataset, and the current data and historical dataset are acquired again after a preset cycle, and then the process proceeds to S2. If the historical dataset can be acquired, the process proceeds directly to S2. The current data includes the current operating data and the current ethyl acetate concentration data at the top of the column. The historical dataset contains historical operating data and historical ethyl acetate concentration data at the top of the column. The operating data includes feed flow rate, ethanol / acetic acid ratio, reflux ratio, and column bottom temperature. The calculation module is used to substitute the current operating data into the ethyl acetate reactive distillation mechanism model to calculate the ethyl acetate concentration at the top of the first theoretical column; The data optimization module is used to apply different operation offsets to the current operation data, substitute each operation data after adding the operation offset into the ethyl acetate reactive distillation mechanism model, calculate the second theoretical column top ethyl acetate concentration corresponding to each operation offset, correct each second theoretical column top ethyl acetate concentration according to the deviation between the first theoretical column top ethyl acetate concentration and the current column top ethyl acetate concentration data, and obtain the third theoretical column top ethyl acetate concentration corresponding to each operation offset. The operation-concentration dataset is composed of different operation offsets, the current operation data, and the third theoretical column top ethyl acetate concentration corresponding to each operation offset. The filtering module is used to sort the ethyl acetate concentration data at the top of the column in the historical dataset from high to low, calculate the optimal direction for concentration increase from the sorted ethyl acetate concentration data at the top of the column, and filter the operation-concentration dataset in S3 according to the calculation result of cosine similarity to obtain the filtered operation-concentration dataset. The results acquisition module is used to calculate the operational optimality of each instance in the filtered operation-concentration dataset, and take the instance with the highest operational optimality as the optimized operation parameter result to complete the real-time optimization of the operation parameters of the ethyl acetate reactive distillation column.

Images