System, method and computer program product for determining operating conditions in a chemical production plant

Abstract

A system and method are provided for determining operating conditions of a chemical production plant including at least one catalytic reactor. Operating data and catalyst aging indicators are received via a communications interface (10). At least one target operating parameter for the operating conditions of a planned or current production run is determined (14). The at least one target operating parameter for the operating conditions is used for monitoring and / or controlling the chemical production plant.

Description

[Technical field]

[0001] The present disclosure relates to systems, methods and computer program products for determining operating conditions in a planned or current production run of a chemical production plant having at least one catalytic reactor. [Background technology]

[0002] Catalytic reactors are widely used in the chemical industry to convert raw materials into valuable chemicals. Catalyst performance (such as conversion, selectivity and yield) is related to the reactor operating parameters and catalyst aging. To achieve certain performance targets, plant operators must adjust the reactor parameters based on experience and engineering knowledge. Often, reactor models based on kinetic, heat and mass transport phenomena are used to represent the catalyst performance and to allow more direct control of the reactor operation. However, such models are complex and experimental determination of kinetic and transport parameters is very tedious and expensive.

[0003] Recently, hybrid types have been presented that rely partly on process knowledge but also use machine learning methods. For example, Clough and Ramirez (AIChEJ 22, 1976, p. 1097) and Gujarati and Babu (Chem. Eng. Sci. 65, 2010, p. 2009) present knowledge-based models of styrene reactors, using them to optimize reactor operation to maximize styrene yield or process economics. In both cases, the kinetic rate equations and coefficients are taken from previous literature. These models are of limited practical use, considering that the kinetic coefficients must be calculated separately for each real catalyst (different catalysts have different reaction rates for all chemical transformations occurring in the system). Furthermore, the aging rate of the catalyst is not taken into account.

[0004] Tamsiian et al. (Comp. Chem. 40, 2012, p. 1) developed a model based on kinetic equations published in the literature and fitted rate coefficients using styrene plant data from a five-day run.

[0005] H. Lim et al. (Ind. Eng. Chem. Res. 43, 2004, p. 6441) developed a hybrid model consisting of a first-principles (knowledge-based) part describing the dependence of catalyst performance (conversion and selectivity) on reactor operating parameters, pressure, temperature, STO ratio and ethylbenzene flow rate as input variables; and a neural network model to predict the deactivation factor used in the first-principles model. The neural network predicts the updated deactivation factor based on the deactivation factor, inlet temperature, total and partial pressures of ethylbenzene and steam for the latest time step.

[0006] Shahhosseini et al. (Int. J. Chem. React. Eng. 9, 2011) developed a hybrid model to optimize the performance of an adiabatic industrial reactor system. First, a kinetic model was fitted using seven data points measured in an experimental isothermal reactor. With the kinetic coefficients thus estimated, they used a Tabu search algorithm or a genetic algorithm to efficiently optimize the operating conditions to maximize two objective functions: ethylbenzene conversion and styrene selectivity. The catalyst activity was modeled as an exponential function of time and decayed from 100% to 40% over a 48-month period based on literature values. Modeling deactivation in this way cannot reflect the effects of different operating scenarios, which is essential for operation-specific and plant-specific deactivation predictions.

[0007] Wu et al. (Lecture Notes in Computer Science, 10357, pp. 301-312, 2017) and WO2018 / 035718A1 disclose a data-driven model for real-time prediction of styrene production based on actual operating data. Thirty-three sensor-based predictors were identified from a combination of principal component analysis (PCA) and random forest variable importance analysis. Three models were compared to predict styrene monomer yield: generalized regression neural network (GRNN), GRNN after Kalman filtering, and random forest regression model. Wu et al. focus on predicting current styrene production based on available sensor data, motivated by the stated fact that the information of available sensor data is only available after analytical testing, which usually takes several hours. However, the presented method is not concerned with performance prediction based on user-defined operating scenarios, or performance forecasting including styrene monomer yield at future time points. Catalyst deactivation is not considered.

[0008] Thus, many knowledge-based or kinetic models have been developed to model styrene production in catalytic reactors, and recently, elements of machine learning (ML) have been integrated (hybrid models), in some cases for the optimization of the kinetic models rather than for performance modeling per se. In all hybrid modeling and machine learning approaches, the ML methods are trained on small data sets, usually consisting of a single experimental plant or a single run of a production plant.

[0009] Knowledge-based or kinetic hybrid reactor models are complex and often require simplifying assumptions to reduce computational effort. In most cases, not all of the physicochemical processes that contribute to catalyst performance and aging are understood. For example, in the case of styrene-based catalysts, it is not possible to accurately estimate the potassium evaporation rate, since the potassium loss is controlled by the temperature gradients occurring along the depth of the catalyst bed and within the catalyst particles, the STO ratio, the carbon dioxide content in the gaseous mixture (resulting from the gasification of coke), and the pressure gradient along the catalyst bed. It is not possible to calculate a priori the potassium loss rate at all points in the reactor bed in order to estimate the catalyst aging rate within the catalyst bed and over the catalyst's lifetime. [Prior art documents] [Patent documents]

[0010] [Patent Document 1] WO2018 / 035718A1 [Non-patent literature]

[0011] [Non-Patent Document 1] Clough and Ramirez (AIChEJ 22, 1976, p. 1097) [Non-Patent Document 2] Gujarati and Babu (Chem. Eng. Sci. 65, 2010, pp. 2009) [Non-Patent Document 3] Tamsiian et al. (Comp.Chem.40, 2012, p. 1) [Non-Patent Document 4] H. Lim et al. (Ind.Eng.Chem.Res.43, 2004, page 6441) [Non-Patent Document 5] Shahhosseini et al.(Int.J.Chem.React.Eng.9,2011) [Non-Patent Document 6] Wu et al. (Lecture Notes in Computer Science, 10357, pp. 301-312, 2017) Summary of the Invention [Problem to be solved by the invention]

[0012] It is an object of the present disclosure to provide a method for determining operating conditions in a chemical production plant including at least one catalytic reactor, which method enables robust, stable and reliable reactor operation and improves process control of the catalyst-based production plant. [Means for solving the problem]

[0013] According to a first aspect of the present invention, there is provided a system for determining operating conditions of a chemical production plant including at least one catalytic reactor, the system comprising: A system including a communications interface and a processing unit in communication with the communications interface, the system comprising: - receiving, via the communications interface, operational data indicative of predefined operating conditions for a scheduled production run or measured operational data indicative of current operating conditions for a current production run, wherein at least one operational data point includes a desired operating value indicative of a change to the current operating conditions; - receiving, via the communications interface, a catalyst aging indicator associated with a period of time that the catalyst has been used in a current or planned production run; - determining, via said processing device, at least one target operating parameter for operating conditions of a planned production run or for operating conditions of a change in a current production run based on the operating data and the catalyst aging indicator by using a data-driven model, preferably a data-driven machine learning model, where the data-driven model is parameterized according to a training data set, where the training data set is based on a historical data set having operating data, a catalyst aging indicator, and at least one target operating parameter; - providing at least one target operating parameter for a scheduled production run or a change in operating conditions in a current production run via the communication interface; It is structured as follows.

[0014] According to another example of the first aspect of the present invention, there is provided a system for determining operating conditions of a chemical production plant including at least one catalytic reactor, the system comprising a communications interface and a processing device in communication with the communications interface. (a) For a scheduled production run, the system: - receiving, via said communications interface, operational data indicative of predefined operating conditions for a scheduled production run; - receiving, via the communications interface, a catalyst aging indicator associated with a period of time during which the catalyst has been used in a scheduled production run; - determining, via the processing device, at least one target operating parameter for an operating condition of the scheduled production operation based on the operating data and the catalyst aging indicator by using a data-driven model, the data-driven model being parameterized according to a training data set, the training data set being based on a historical data set having operating data, a catalyst aging indicator, and at least one target operating parameter; - providing, via said communication interface, at least one target operating parameter for an operating condition of a scheduled production run; or (b) due to a change in a current production run, the system: - receiving, via the communications interface, measured operating data indicative of current operating conditions for a current production operation, at least one operating data point including a desired operating value indicative of a change in the current operating conditions; - receiving, via said communications interface, a catalyst aging indicator associated with a period of time that a catalyst has been used in a current production run; - determining, via the processing device, at least one target operating parameter for changing operating conditions in a current production operation based on the operating data and the catalyst aging indicator by using a data-driven model, the data-driven model being parameterized according to a training data set, the training data set being based on a historical data set having the operating data, the catalyst aging indicator, and at least one target operating parameter; - providing at least one target operating parameter for a change in operating conditions in a current production run via the communication interface; It is structured as follows.

[0015] According to a second aspect of the present invention, there is provided a system for optimizing operating conditions of a planned production run or changes in an existing production run of a chemical production plant, the system comprising the system described above and an optimization processor, the optimization processor comprising: - receiving, via the communications interface, determined target operating parameters for one or more operating conditions of a scheduled production run or a change in a current production run; - determining, via the optimization processor, based on the received target operating parameters for each operating condition, a maximum or minimum value of the target operating parameters or a maximum or minimum value of at least one optimization parameter derived from the target operating parameters; - providing, via said communication interface, maximum or minimum values ​​indicative of optimal operating conditions for a planned or current production run, such as maximum or minimum values ​​of target operating parameters or maximum or minimum values ​​of optimization parameters derived from the target operating parameters; It is structured as follows.

[0016] According to another example of the second aspect of the present invention, there is provided a system for optimizing operating conditions of a chemical production plant, the system comprising the above-mentioned system and an optimization processing device, the optimization processing device comprising: - receiving, via the communications interface, determined target operating parameters for a plurality of operating conditions determined for (a) a scheduled production run, or (b) a change in a current production run; - determining, via the optimization processor, for each operating condition based on the received target operating parameters, a minimum or maximum value of the target operating parameters or a minimum or maximum value of an optimization parameter derived from the target operating parameters; - providing, via said communication interface, minimum or maximum values ​​indicative of optimal operating conditions for either (a) a planned production run or (b) a current production run, e.g., a minimum or maximum value of a target operating parameter or a minimum or maximum value of an optimization parameter derived from the target operating parameter; It is structured as follows.

[0017] According to a third aspect of the present invention, there is provided a production monitoring and / or control system, the system including a communication interface communicatively coupled to the system, e.g., via a wired or wireless connection, for determining the above-mentioned operating conditions or for optimizing the above-mentioned operating conditions. The production monitoring and / or control system may include a display device, the display device configured to receive and display the determined operating conditions. The production monitoring and / or control system may include a control unit, the control system configured to receive the determined operating conditions and control a current or scheduled production operation of a chemical production plant based on the determined operating conditions. The determined operating conditions preferably include determined target operating parameters and operational data, optionally further including any desired operating values.

[0018] According to a fourth aspect of the present invention, there is provided a computer-implemented method for determining operating conditions for a planned production run, or for a change in an existing production run, of a chemical production plant including at least one catalytic reactor, the method comprising the steps of: - receiving, via said communications interface, operational data indicative of predefined operating conditions of a scheduled production run or measured operational data indicative of a current production run, wherein at least one operational data point includes a desired operational value indicative of a change in a current operating condition; - receiving, via said communications interface, a catalyst aging indicator relating to the duration of use of the catalyst in a current or planned production run; - determining, via a processing device, at least one target operating parameter for operating conditions of a planned production run or a change in a current production run based on the operating data and the catalyst aging indicator by using a data-driven machine learning model, the data-driven model being parameterized according to a training data set, the training data set being based on a historical data set having the operating data, the catalyst aging indicator, and at least one target operating parameter; - providing, via said communication interface, at least one target operating parameter for a scheduled production run or a change in operating condition in a current production run.

[0019] According to another example of the fourth aspect of the present invention, there is provided a computer-implemented method for determining operating conditions of a chemical production plant including at least one catalytic reactor, the method comprising the steps of: (a) For a scheduled production run, the method includes the steps of: - receiving, via a communication interface, operating data indicative of predefined operating conditions for a scheduled production run; - receiving, via said communications interface, a catalyst aging indicator associated with a period of time that the catalyst has been used in a scheduled production run; - determining, via a processing device, at least one target operating parameter for an operating condition of a planned production operation based on the operating data and the catalyst aging indicator by using a data-driven model, the data-driven model being parameterized according to a training data set, the training data set being based on a historical data set having the operating data, the catalyst aging indicator, and at least one target operating parameter; - providing, via said communication interface, at least one target operating parameter for an operating condition of a scheduled production run; or (b) a change in a current production operation, the method comprising the steps of: - receiving, via a communications interface, measured operating data indicative of current operating conditions for a current production run, wherein at least one operating data point includes a desired operating value indicative of a change to the current operating conditions; - receiving, via said communications interface, a catalyst aging indicator associated with a period of time that a catalyst has been used in a current production run; - determining, via a processing device, at least one target operating parameter for changing operating conditions in a current production operation based on the operating data and the catalyst aging indicator by using a data-driven model, the data-driven model being parameterized according to a training data set, the training data set being based on a historical data set having the operating data, the catalyst aging indicator, and at least one target operating parameter; - providing, via said communications interface, at least one target operating parameter for a changed operating condition in a current production operation.

[0020] According to a fifth aspect of the present invention, there is provided a method for optimizing operating conditions of a chemical production plant for a planned production run or changes in an existing production run, the method comprising the steps of: - receiving, via a communication interface, determined target operating parameters for one or more operating conditions for a scheduled production run or a change in a current production run; - determining, via an optimization processor, based on the received target operating parameters for each operating condition, a maximum or minimum value of the target operating parameters or a maximum or minimum value of an optimization parameter derived from the target operating parameters; - providing, via said communication interface, a maximum or minimum value indicative of an optimal operating condition for a planned or current production run, such as a maximum or minimum value of a target operating parameter or a maximum or minimum value of an optimization parameter derived from the target operating parameter.

[0021] According to another example of the fifth aspect of the present invention, there is provided a method for optimizing operating conditions of a chemical production plant for a planned production run or a change in an existing production run, the method comprising the steps of: - receiving, via a communications interface, determined target operating parameters for a plurality of operating conditions determined for (a) a scheduled production run or (b) a change in a current production run; - determining, via an optimization processor, based on the received target operating parameters for each operating condition, a maximum or minimum value of the target operating parameters or a maximum or minimum value of an optimization parameter derived from the target operating parameters; - providing, via said communications interface, minimum or maximum values ​​indicative of optimal operating conditions for either (a) a planned production run, or (b) a current production run, e.g., a minimum or maximum value of a target operating parameter, or a minimum or maximum value of an optimization parameter derived from the target operating parameter.

[0022] According to a sixth aspect of the present invention, there is provided a method for monitoring and / or controlling a chemical production plant, the method comprising the step of determining or optimizing the operating conditions as described above. The method may comprise displaying the determined operating conditions on a display device and / or controlling a current or planned production operation in the chemical production plant based on the determined operating conditions. The determined operating conditions preferably comprise determined target operating parameters, and / or operating data optionally including desired operating values.

[0023] According to a seventh aspect of the present invention, there is provided a method for training a data-driven model to determine operating conditions for a planned production run, or for changes in a current production run, of a chemical production plant including at least one catalytic reactor, the method comprising the steps of: - receiving, via a communications interface, a training data set based on a historical data set including the operating data, the catalyst aging indicator, at least one target operating parameter, and optionally plant metadata; - training, via a processing device, a driven model by adjusting parameter settings according to a training data set; - providing the trained data-driven model via a communication interface.

[0024] According to an eighth aspect of the present invention, a computer program or computer program product, or a computer readable instruction non-volatile storage medium containing computer readable instructions, when loaded and executed by a processing device, performs the methods disclosed herein.

[0025] According to a ninth aspect of the present invention, there is provided a catalyst comprising a catalyst type identifier associated with a data-driven model trained according to the methods described herein, the model being trained for a catalyst type indicated by the catalyst type identifier. In other words, there is provided a catalyst system comprising a catalyst and a catalyst type identifier associated with a data-driven model trained according to the methods described herein, the model being trained for a catalyst type indicated by the catalyst type identifier.

[0026] According to a tenth aspect of the present invention, there is provided a catalyst comprising a catalyst type identifier associated with a computer program as described herein. In other words, there is provided a catalyst system comprising a catalyst and a catalyst type identifier associated with a computer program as described herein.

[0027] According to an eleventh aspect of the present invention there is provided a chemical reaction process using a catalyst associated with a data-driven model trained in accordance with the present disclosure, said data-driven model being used to design plant components and optimise operation of the chemical reaction process to achieve target performance, particularly for operating conditions of a chemical production plant including at least one catalytic reactor.

[0028] The following disclosure applies to the systems, methods, computer programs, computer readable non-volatile storage media, catalysts, chemical reaction processes and computer program products similarly disclosed herein. Therefore, no distinction is made between the systems, methods, computer programs, computer readable non-volatile storage media, or computer program products. All features are disclosed in relation to the systems, methods, computer programs, computer readable non-volatile storage media, catalysts, chemical reaction processes and computer program products disclosed herein.

[0029] The present invention provides a system or method for determining operating conditions or parameters solely based on data-driven models, which allows for more robust, stable and reliable reactor operation and enhanced process control in catalyst-based production plants, and is more accurate than known approaches based on kinetic or hybrid models, in particular because there is no a priori information on reaction kinetics and therefore no estimates or simplifying assumptions about the underlying process are required.

[0030] Specifically, the determination takes into account catalyst deactivation by providing a catalyst aging indicator related to the period of time the catalyst has been used in the current production run. By including the catalyst aging indicator in the input parameters of the model, the operating conditions can be determined more accurately because the model inherently accounts for catalyst deactivation or aging. By accounting for catalyst deactivation, the determination of the operating parameters can be applied to a wide range of chemical production plants and operating parameters. In contrast, driven models that do not account for catalyst aging have limited applicability. In particular, such models are limited to the chemical production plant on which the model was trained, sometimes only applicable to a single run and a narrow range of operating conditions at the production plant.

[0031] The disclosed system and method performs short-term discrete forecasting and time series forecasting, the latter of which can cover the catalyst deactivation process during the entire production run. The method and system allow the plant operator to improve and optimize the daily operating policy according to the energy cost, the market supply of raw materials or the demand for the plant's products, and other constraints that may occur within the plant, such as the shutdown of various plant parts or equipment.

[0032] Additionally, the system and method allows for predicting short-term and / or forecasting long-term behavior of a catalyst-based production plant using plant operating data and catalyst aging indicators, thereby enhancing process control in a catalyst-based production plant, for example, by allowing a plant operator to easily evaluate planned or current operations based on existing operating conditions in the production plant. Thus, the system provides a powerful tool for planning, monitoring and control of production processes.

[0033] In one embodiment of the present invention, the chemical production plant comprises one or more catalytic reactors. In the context of the present invention, the term "catalytic reactor" means a chemical reactor, typically containing a catalyst, in which catalytic chemical reactions take place. The catalytic reactor may be a fixed bed catalytic reactor. The chemical production plant may be a styrene production plant.

[0034] In this context, a data-driven model, preferably a data-driven machine learning model or simply a data-driven model, refers to a trained mathematical model that is parameterized according to a training data set to reflect the reaction kinetics or physicochemical process of a chemical production plant or a catalytic reactor. An untrained mathematical model refers to a model that does not reflect the reaction kinetics or physicochemical process, e.g., an untrained mathematical model is not derived from physical laws that provide scientific generalizations based on empirical observations. Thus, the reaction kinetics or physicochemical properties may not be inherent to the untrained mathematical model. An untrained model does not reflect such properties. The untrained mathematical model can be parameterized by feature engineering and training with the respective training data set. The result of such training is a simply data-driven model, preferably a data-driven machine learning model, that reflects the reaction kinetics or physicochemical properties as a result of the training process, preferably only as a result of the training process.

[0035] Historical data in this context refers to a data set including at least operating data, catalyst aging indicators, and at least one target operating parameter, each data set being associated with one production run. Thus, each data set includes data related to a production run of one catalyst from start to finish. Such data may be measured and recorded during the production run over the life of the catalyst, for example, from start to finish after catalyst replacement to finish when the catalyst again requires replacement.

[0036] The chemical reaction catalyzed is not limiting to the scope of the present invention. As an example, the catalyst may be a dehydrogenation catalyst. In particular, the catalyst of the catalytic reactor may be an iron oxide-based catalyst for dehydrogenating aliphatic or alkyl aromatic hydrocarbons to form the corresponding unsaturated hydrocarbons. Examples of such dehydrogenation processes include the dehydrogenation of ethylbenzene to styrene, the dehydrogenation of isopropylbenzene to alpha-methylstyrene, the dehydrogenation of butene to butadiene, or the dehydrogenation of isoamylene to isoprene. The method and system of the present invention are particularly useful in styrene production plants. Preferably, the chemical production plant is a styrene production plant, which converts ethylbenzene to styrene using a styrene catalyst. The preparation of styrene typically involves a heterogeneously catalyzed dehydrogenation of ethylbenzene in the presence of steam. The catalytic dehydrogenation of aliphatic or alkyl aromatic hydrocarbons is usually carried out industrially at temperatures ranging from 500 to 700° C. in the presence of steam. In these processes, the hydrocarbon and steam are typically mixed and passed over an iron oxide dehydrogenation catalyst at high temperature and low pressure.

[0037] The term operational data refers to quantities indicative of the operational status of a chemical production plant. In particular, such quantities relate to measurement data collected during production operations of the chemical production plant and may be derived directly or indirectly from such measurement data. In a preferred embodiment, the operational data comprises sensor data measured via sensors installed in the chemical production plant, quantities derived directly or indirectly from such sensor data, analytical data measured on samples taken from the chemical production plant, quantities derived directly or indirectly from such analytical data, or combinations thereof.

[0038] The sensor data may include measured quantities available in the chemical production plant by installed sensors, e.g. temperature sensors, pressure sensors, flow rate sensors, etc. The analytical data may include quantities provided from analytical measurements of samples extracted at any point in the process or over time from the chemical production plant. In particular, such analytical data may include compositions of reactants, starting materials, products and / or by-products, determined e.g. via gas chromatography from samples extracted during the production process, at different stages of the production process, e.g. before or after a catalytic reactor. The analytical data preferably form the basis for determining the performance characteristics of the catalyst.

[0039] The set of operational data may include "raw data," which refers to basic unprocessed analytical and / or sensor data, or process or derived parameters derived directly or indirectly from the raw data. In the case of a chemical production plant, derived parameters may include average inlet temperature across multiple catalytic reactors derived from respective temperature sensors; steam to oil ratio derived from raw data of steam and reactant flow rates; conversion and selectivity derived from analytical data before and after the reactor; any type of normalization data, such as production values ​​normalized by catalyst amount or catalyst mass; any data derived from time series data, such as cumulative production, maximum load to date, or combinations thereof.

[0040] In particular in chemical production plants, conversion, selectivity and yield are derived from analytical data. Conversion here refers to the percentage of reactants that are converted fractionally, preferably completely, in the reactor. For example, in the case of styrene production, this corresponds to converting the starting material ethylbenzene to a desired product. The selectivity of the desired product refers to the amount of the converted reactants that are converted to the desired product. For example, in the case of styrene production, this corresponds to the selectivity of the reaction from ethylbenzene to styrene. The yield of the desired product refers to the mathematical product of the conversion and the product-specific selectivity. The yield can be expressed by the percentage of the reactants charged to the reactor that are converted to the desired product.

[0041] In a further embodiment, plant metadata indicative of a physical plant layout is received via the communication interface. The plant metadata may include reactor characteristics, e.g., predefined by the physical plant layout, describing reactor characteristics, etc., that may be relevant to plant or reactor performance. The plant metadata may include, e.g., the number of reactors through which the reactant mixture subsequently passes (e.g., 2 or 3 reactors), the total amount of catalyst, the catalyst volume per reactor, the dimensions of each reactor (length, diameter, height...), the type of catalyst used in the plant, or a combination thereof. In a further embodiment, the determination of at least one target operating parameter via the processing unit is additionally based on the plant metadata using a data-driven model, and the learning data set is additionally based on a historical data set including the plant metadata. Including the plant metadata in the determination of the operating conditions allows for building a data-driven model that is applicable to different plants, thereby increasing the number of data points available for training the data-driven model. As a result, the data-driven model may broadly capture the operating conditions of different plants operating under different operating conditions with different physical plant layouts, allowing for more accurate determination.

[0042] In further embodiments, the historical data set includes data from multiple production runs, multiple plants, and / or multiple catalyst batches of the same type of catalyst. Including multiple production runs in the learning can cover different operating conditions of the same or different plants. Including data from multiple catalyst batches can account for differences between catalyst batches. Data from multiple catalyst batches can include data from multiple production runs, at least one production run using a different production batch of the same type of catalyst. Including multiple plants includes data from one or more production runs from different plants. In this way, covering operating conditions of different plants can provide a wider scope of applicability for the model. In this context, the same type of catalyst refers to the same type of catalyst formulation. Multiple catalyst batches include the same type of catalyst provided from different production batches or different delivery dates.

[0043] In further embodiments, the catalyst aging indicator is based on a time, period, amount derived from the time-dependent operating data, and / or an amount derived cumulatively from the time-dependent operating data. The catalyst aging indicator can be defined by the time spent by the catalyst in the reactor at the reaction conditions since first contacting the reactant mixture. Additionally or alternatively, the cumulative loading or cumulative production of the catalyst may be used as the catalyst aging indicator, preferably defined by the reactant feed amount or the total amount of reactant converted from the start of the operation to the predicted start time. Apart from the indicators mentioned here, any other amount that can function as a catalyst aging indicator may be utilized. The catalyst aging indicator may be provided via a client device, where the plant operator inputs the time of the start of the production run or the time period from the start of the production run to the predicted start point. Alternatively or additionally, the catalyst aging indicator may be determined based on a time series of operating data, preferably since the start of the production run, which may be determined based on an operating profile in the operating data. Such an operating profile may include a particular temperature, pressure, flow rate profile, or a combination of such profiles.

[0044] Projected or changed operating conditions in a current production run refer to operating conditions under which a chemical production plant may operate in the future or after a forecast start time. Such operating conditions may include operating conditions for one or more of a discrete time point, a plurality of discrete time points, a period of time, or a plurality of time intervals. In the latter case, the operating conditions may include different time intervals during which at least one predefined or desired operating value assumes different values ​​at the different intervals. The plurality of discrete time points, periods, or time intervals may extend from the forecast start time through the remainder of the production run to the end of the production run, at which point replacement of the catalyst is required.

[0045] The operating conditions for the planned production run may also refer to operating data of the operating conditions of the chemical production plant that are predefined before the start of the production run. The predefined operating conditions may include a set of predefined operating data that defines the operating conditions. Determining the operating conditions for the planned production run is particularly useful for production design and planning before the start of the production run. The operating conditions for the current production run may refer to changes in the operating conditions currently set in the chemical production plant. The operating conditions may include a set of operating data that defines the current operating conditions, the operating data being at least one operating data point that differs from the operating parameters as currently set in the chemical production plant. Determining the operating conditions for the current production run is particularly useful for monitoring and controlling production during the current production run.

[0046] Determining operating conditions involves predicting or forecasting the performance / behavior of a catalyst or reactor and is commonly used to illustrate the application of a data-driven model to an appropriate set of input parameters. In the example of a styrene production plant, determining refers to determining the reactor inlet temperature and associated selectivity of the ethylbenzene to styrene reaction given certain operating conditions, for example, via operational data. The selection of output parameters for such a model may be driven by the fact that most plants are operated by operators adjusting the inlet temperature to achieve a particular desired conversion. However, given the same data set, models can be developed in a similar manner with different allocations of inputs and outputs, such as using the reactor temperature as an input to predict the expected conversion under the operating conditions, as further detailed below.

[0047] In one embodiment, the determination of at least one target operating parameter for the operating conditions is based on a short-term model that determines target operating parameters for discrete time points, or a long-term model that determines target operating parameters, in particular for a certain period of time in the future. In the case of a short-term model, the at least one target operating parameter for the operating conditions is determined based on a short time frame, where a short time frame can refer to a discrete or single time point. The short-term model may be based on different machine learning techniques, including regression models, such as linear regression models, non-linear regression models, Bayesian linear regression, random forest models, neural networks, or combinations thereof. Other approaches may also be applied. Preferably, the short-term model does not have an inherent time dependency. More preferably, the short-term model is capable of predicting performance or operating conditions at time t based on model inputs at the same time t.

[0048] In the case of a planned production run, the operating data may define predefined operating conditions. The catalyst aging indicators may be estimated for the planned production run, for example, assuming that the predefined operating conditions are constant over time or change over a predefined time period, or based on a previous production run. The short-term model determination may be applied at multiple discrete time points, for example, for predefined operating conditions that are constant over the entire production run, and the catalyst aging indicators may be estimated for each of the discrete time points. Such an implementation is particularly advantageous for designing and planning future production runs.

[0049] In the case of a current production operation, the operating data may define the currently set operating conditions in the chemical production with at least one desired operating value indicating a change or deviation from the current operating conditions. The determination via the short-term model may be applied at one or more discrete points in time to determine the effect of at least one desired operating value indicating a deviation from the currently set operating conditions in the chemical production, for example. Such an implementation is particularly advantageous for monitoring and controlling the current production operation.

[0050] For long-term models, at least one target operating parameter for the operating conditions is determined based on a long time frame, where a long time frame can refer to multiple time points in the future. The number of points, and therefore the forecast horizon, depends on the time scale of the time dynamics in the production plant. For catalyst-based production processes, such dynamics can be determined by the aging dynamics of the catalyst and its time axis. For heterogeneous catalytic reactions such as styrene production, such time scales can be in the range of weeks, months, or years. Long-term models can be based on time series forecasting methods. Such methods include, for example, known regression methods such as autoregressive models, in particular autoregression (AR), moving average (MA), autoregressive moving average (ARMA), autoregressive integrated moving average (ARIMA), vector autoregression (VAR), vector autoregressive moving average (VARMA), vector autoregressive moving average with exogenous regressors (VARMAX), random forest models, neural networks, convolutional neural networks, recurrent neural networks, or combinations thereof. Other approaches are also applicable.

[0051] In one embodiment, a time series of at least one target operating parameter measured, predicted, or derived during a current or previous production run up to the prediction start time is received via the communication interface. For example, for a planned operation, a target operating parameter measured, predicted, or derived during a previous production run can be used. Preferably, the at least one target operating parameter is determined at one or more time points following the prediction start time based on the operating data, optionally including a desired operating value, the time series of the at least one target operating parameter, the catalyst aging indicator, and optionally the plant metadata using a data-driven model. More preferably, the data-driven model includes an inherent time dependency. In this embodiment, the at least one target operating parameter may include an uncontrollable or endogenous parameter that cannot be controlled in the chemical production plant via machine settings. In contrast, the operating data, including the at least one desired operating value, may include a controlled or exogenous parameter in the chemical production plant via machine settings. Determining the target operating parameters at the prediction start time t may include, for example, determining them for times t, ..., t+N (N>0) based on the non-controlled parameters or a subset of the non-controlled parameters up to time t-1 and optionally based on the controlled parameters or a subset of the controlled parameters, for example at times t, ..., t+N, optionally further including, for example, t-1, t-2, ... or other times suitable for making the prediction with the particular structure of the selected model.

[0052] Preferably, the long-term models are time-series-based models that include an inherent time dependency and predict target operating parameters at times t, ..., t+N based on model inputs at times up to time t-1, where t refers to the start of the prediction in time. In this way, the long-term models are able to determine reactor or catalyst performance and changes in reactor or catalyst performance, preferably over long periods of time, even if there is no available information of uncontrolled parameters beyond the start of the prediction. These models have at least some inherent time dependency and predict performance at times t, ..., t+N based on uncontrolled parameters up to time t-1 and optionally controlled parameters, for example at times t, ..., t+N or other times suitable for making the prediction with the particular structure of the selected model.

[0053] In one embodiment, the processing device is further configured to pre-process the operational data prior to determination via the data-driven model. Preferably, the pre-processing includes a transformation into quantities independent of the physical plant layout. The pre-processing allows taking into account systematic and non-systematic differences between different plants. Thus, a broad applicability of the determination of the target operational parameters is possible, even if data from the particular plant in question was not used to train the data-driven model. In particular, the transformation includes systematic factors as input parameters of the data-driven model and / or normalization of the operational data.

[0054] In a further embodiment, the data-driven model is validated before determining the at least one target operating parameter. Such validation increases the interpretability and reliability of the determined target operating parameter. For validation, the at least one target operating parameter and the operating data measured or derived during the current production operation at one or more time points can be received via the communication interface. The operating data can be used to determine the at least one target operating parameter corresponding to the time point at which the measured target operating parameter is available. The result of the determined target operating parameter at a particular time point can be compared with the measured operating parameter corresponding to the same time point. If the comparison leads to a valid model operation, e.g., if the difference is less than a threshold value or if no systematic error is identified, the determination of the at least one target operating parameter for the operating conditions can continue based on a set of operating data including at least one desired operating value. If the comparison results in an invalid model operation, a warning can be issued, e.g., via a display or audio device, indicating that the model operation is not suitable for monitoring and / or controlling the target chemical plant at the current operating conditions.

[0055] In a further embodiment, the data-driven model is selected based on the type of catalyst before determining the at least one target operating parameter. The type of catalyst may be received via metadata indicating the type of catalyst used, for example via a catalyst type identifier. The type of catalyst may specify a catalyst formulation. Such an option allows the system, computer program, computer program product, and method to be used very flexibly, covering different catalyst types as well as different chemical production plants. Further flexibility can be added by selecting the data-driven model based on input and output parameter identifiers before determining the at least one target operating parameter. The input and output parameter identifiers can specify which parameters are used as operating data and which are used as target operating parameters in the system or method. In such an implementation, different data-driven models may be stored in the memory of the system, including the catalyst type and / or parameter identifiers of each model, based on which the selection step is performed.

[0056] In one embodiment, the training of the data-driven model is performed based on training data including a catalyst type identifier. The trained model may be stored with the catalyst type identifier and preferably provided with the catalyst type including the catalyst type identifier. Such identification may be performed electronically, for example via a catalyst ID stored in a database or a mobile storage medium attached to a catalyst supply vessel in association with the respective chemical production plant. In such a case, the model is preferably trained based on a training data set measured in a production run using a catalyst of the same catalyst type as indicated by the catalyst type identifier. Thus, the catalyst can be bundled into the data-driven model to assist in a more robust control of production plants using catalysts of a particular catalyst type.

[0057] A method is provided for predicting the short-term behavior and forecasting the long-term behavior of a catalyst, for example in a fixed-bed catalytic reactor, as a function of reactor operating parameters or data and catalyst life or catalyst aging indicators. The method uses a data-driven model, preferably a data-driven machine learning model, and does not contain a priori information on reaction kinetics. The model is capable of predicting both the short-term and long-term behavior of the catalyst as a function of input parameters, such as typical reactor operating parameters or data and parameters derived from sensor and analytical raw data available in the production plant. A software product for implementing the method is also provided. As an application example, the method is used to predict and forecast the behavior of catalysts and technical reactors for the conversion of ethylbenzene to styrene.

[0058] The present disclosure provides a computer-implemented method for predicting short-term performance and forecasting long-term performance of a catalyst in a chemical production plant including a catalytic reactor, including catalyst aging effects. The method includes using a mathematical model of the chemical production plant, particularly the catalytic reactor, that is based on machine learning and does not include a priori information on reaction kinetics, and uses input parameters selected from sensor raw data, derived parameters, reactor operating parameters or data, plant metadata, and parameters indicative of catalyst life or catalyst aging indicators.

[0059] In one embodiment, the operational data is selected from sensor data available from a chemical production plant, particularly a catalytic reactor, such as analytical data from a gas chromatography (GC) analysis, and the derived parameters described above.

[0060] In one embodiment, the operating data includes the inlet and outlet temperatures of the one or more catalytic reactors, preferably the inlet and outlet temperatures of each of the one or more catalytic reactors, the inlet and outlet pressures of the one or more catalytic reactors, preferably the inlet and outlet pressures of each of the one or more catalytic reactors, and the composition of the reaction mixture at the inlet and outlet of the one or more catalytic reactors, preferably the composition of the reaction mixture at the inlet and outlet of each of the one or more catalytic reactors.

[0061] In one embodiment, specifically applied to a styrene production plant, the operational data includes steam-to-oil ratio (STO), liquid hourly space velocity (LHSV), total styrene production normalized by catalyst load, target ethylbenzene conversion, styrene selectivity, average inlet temperature, normalized pressure after the last catalytic reactor, normalized pressure drop in one or more catalytic reactors, temperature loss in one or more catalytic reactors, normalized deviation of temperature loss compared to expected (calculated based on conversion) value. The advantage of such a parameter set is that it remains interpretable, since many of the parameters correspond to or can be easily interpreted in the context of actual operating parameters. Other approaches to reduce the dimensionality of the problem (such as PCA or RFA) may result in parameters that are useful for the modeling process and prediction accuracy, but are often poorly interpretable.

[0062] The data-driven model, preferably performing time series forecasting, can fully cover the catalyst deactivation process in the entire production run in the reactor system, allowing the plant operator to improve and optimize the daily operating policy according to the energy cost, market supply of raw materials or demand for the plant products, and other constraints that may occur in the plant, such as the shutdown of various plant parts or equipment.

[0063] In one embodiment, the forecast period covers the remaining life of the catalyst as determined from the operating conditions constraints, preferably the maximum reactor temperature. The model allows the remaining life of the catalyst in use to be determined based on the expected operating conditions and the constraints of these operating conditions, e.g., the maximum reactor temperature or pressure. This allows the operator to plan with confidence the remaining time until the next catalyst change, and also allows the simulation of different operating policies to extend the catalyst life, if necessary, or to obtain maximum production from the remaining time until the next planned plant shutdown.

[0064] In one embodiment, the data-driven model is used to predict the operating conditions of a catalyst in a production plant that does not provide historical data to train the data-driven model.

[0065] In one embodiment, the output of the data-driven model is used to optimize operating conditions for a planned or current production run of the chemical production plant. In such an embodiment, target operating parameters determined via the data-driven model for one or more operating conditions are received and provided to an optimization processor. The optimization may include one or more optimization goals. The optimization goal may be identified, for example, by an optimization parameter or a target parameter to be optimized. In this context, the optimization goal may further include finding a minimum or maximum value for a particular optimization parameter or target parameter. In addition to the target parameters for the multiple operating conditions, the optimization goal may be received and provided to the optimization processor. For example, an operator may be provided with a selection of possible optimization goals, and a selected optimization goal based on a user selection may be received.

[0066] In the case of one or more optimization objectives, one optimal solution may exist and be provided as a result of the optimization. Such an optimal solution may be provided to a monitoring and / or control system. The optimal solution may be displayed on a display device and may be used to control the chemical production process. In the case of multiple optimization objectives or multi-objective optimization, one or more optimal solutions may exist and be provided as a result of the optimization. The multiple optimal solutions may be provided to a monitoring and / or control system and displayed on a display device. In such cases, a plant operator may choose between multiple optimal solutions, simplifying the decision-making process in the complex context of operating a chemical production plant.

[0067] In a further embodiment, the optimization processor determines a minimum or maximum value of the target operating parameter or a minimum or maximum value of at least one optimization parameter derived from the target operating parameter based on the target operating parameter output by the data-driven model. In one exemplary scenario, the remaining life of the catalyst may be an optimization parameter derived from the target operating parameter, and finding its maximum value is the goal of the optimization process. Furthermore, constraints can also be included in the optimization problem. Constraints on the target operating parameter side can include a maximum reactor temperature, a minimum selectivity, or a minimum production volume per day. In another exemplary scenario, the remaining life of the catalyst or the production volume until the time scheduled for catalyst replacement may be an optimization parameter derived from the target operating parameter, and finding its maximum value is the goal of the optimization process. Constraints on the target operating parameter side can include a maximum reactor temperature, a minimum selectivity, or a minimum production volume per day, and optionally a maximum production volume.

[0068] In another exemplary scenario, multiple objectives can be combined. For example, the remaining life of the catalyst at production beyond the catalyst life becomes an optimization parameter derived from the target operating parameters, and finding the combination of maximum remaining life and maximum production becomes the goal of the optimization process. Constraints on the target operating parameters side can include maximum reactor temperature, minimum selectivity, or minimum production per day. In a further example, the multiple objectives further include optimal timing of catalyst replacement, further considering cost aspects such as remaining production, energy requirements, catalyst replacement cost, or a combination thereof. For multi-objective optimization, known Pareto optimization techniques can be used.

[0069] In one embodiment, the data-driven model is used to simulate the expected catalyst performance, production rate (e.g., total amount of styrene produced in a period of time), energy requirements (e.g., related to heating of the reactor or steam reactor), and costs or benefits (e.g., related to reactant costs, market price of products, energy costs, etc.) for a selected set of operating data or operating parameters in a catalytic reactor. Here, production yields, energy consumption, CO2 emissions, costs, by-products, catalyst replacement intervals, remaining life, or combinations thereof can be target operating parameters or optimization parameters derived from target operating parameters. In the optimization, different constraints can be defined, such as, for example, constraints on the operating data on which the determination of the target parameters is based.

[0070] In one embodiment, the data-driven model is used to simulate expected catalyst performance, production rates, and energy requirements for a set of operating data that cannot actually be achieved at a given plant setup (e.g., lower pressure levels (deeper vacuum), lower STO ratios, or additional reactors). This helps plant managers better evaluate the economics of potential upgrades to their plants.

[0071] In one embodiment, the data-driven model is used to simulate and / or estimate the expected performance and operating conditions of a new plant where the catalyst has never been used before.

[0072] The present disclosure also provides a computer program product configured to carry out the method of the present disclosure. In one embodiment, the computer program product is a computer program implemented in a chemical production plant or a styrene production plant, in particular in a computing unit (computer) integrated therein and / or connected thereto. In one embodiment, the computer program product is integrated into a dashboard of the chemical production plant or the styrene production plant.

[0073] In one embodiment, the computer program product comprises an interface for inputting operating conditions of the chemical production plant. The operating conditions can be historical data from a previous production run or actual operating conditions of the chemical production plant in a current run.

[0074] In one embodiment of the computer program product, input parameters or operating data for the forecasting are manually uploaded to the computer program product and / or read from the process control system by the computer program product. In one embodiment, the input parameters or operating data are provided in a formatted data table. In one embodiment, a graphical user interface is provided for manually uploading data to the computer program product. In another embodiment, an application programming interface is provided for uploading data to the computer program product.

[0075] In one embodiment, the user of the computer program product provides all the information specified above (e.g. raw data, corresponding units, plant metadata) to the computer program product. This information can be entered manually, uploaded via a structured data file, or provided via an application programming interface (manual or automatic). For long-term forecasting, the input time series of target operating parameters preferably cover at least the range required by the time lag structure used in the model. Control parameters may also be included if time lags are also used in the model.

[0076] In one embodiment, exogenous operating parameters are input by a user, for example for multiple operating scenarios for which performance is desired to be predicted. All additional data processing, including formatting, aggregation and prediction, is performed by the software product.

[0077] In one embodiment, the predictions generated by the computer program product are presented to the user via a graphical user interface, as structured data files, or via an application programming interface.

[0078] According to a further embodiment, a computer program product is a computer program product which, when loaded into a memory of a computing device and executed by at least one processor of the computing device, performs the computer-implemented method steps described above.

[0079] The computer program product may be a stand-alone unit or may be used in conjunction with or incorporated in a computer system that includes one or more remote terminals or devices that communicate with a central computer via a network, such as the Internet or an intranet. Thus, the computer or processor and related components described herein may be part of a local computer system, a remote computer or an online system, or a combination thereof. Any databases and computer program products described herein may be stored in the computer's internal memory or on a non-transitory computer readable medium.

[0080] Another aspect of the present disclosure is a computer system for predicting the performance of a catalyst in a chemical production plant or for determining the operating conditions of the chemical production plant. The computer system comprises at least an interface component configured to access and read operational parameters or operational data and catalyst specific parameters, in particular catalyst aging indicators, and a processor unit configured to implement a data-driven model and predict the performance of the catalyst by feeding the reactor operational parameters or operational data and catalyst specific parameters provided via the interface component to the data-driven model. In one embodiment, the computer system is coupled to a chemical production plant including a catalytic reactor via a wired and / or wireless communication connection and is configured to access and read the reactor operational parameters or operational data and / or catalyst specific parameters at least partially automatically from a process control system of the chemical production plant including the catalytic reactor via the interface component.

[0081] A further aspect of the present disclosure is a computer-implemented method for training a data-driven model of a chemical production plant including one or more catalytic reactors based on machine learning to predict or forecast the performance of a catalyst in the catalytic reactor. The method includes providing a mathematical model as an initial basis, providing historical data, for example from multiple production runs of the same type of catalyst and / or from multiple production runs of multiple chemical production plants or catalytic reactors including the same type of catalyst, and optionally previously determined target operating parameters or operating and performance parameters of the chemical production plant, in particular the catalytic reactor, accessing and importing the provided historical data into the mathematical model, adapting parameter settings of the data-driven mathematical model to the provided historical data, providing an update of the data-driven mathematical model based on the adapted parameter settings, and setting the updated data-driven mathematical model as an initial basis and iteratively repeating the method steps.

[0082] Another aspect of the disclosure is a computer-implemented method for determining operational and / or performance parameters or target operational parameters of a chemical production plant including one or more catalytic reactors, the method including accessing sensor data indicative of operating conditions present in a reactor system and catalyst-specific sensor data indicative of a catalyst currently being used in the reactor system, determining operational and / or performance parameters or target operational parameters of the chemical production plant using a data-driven model, where the data-driven model is parameterized according to a training data set including historical data, e.g., historical data from multiple production runs of the same type of catalyst, and / or historical data from multiple production runs in multiple catalytic reactors including the same type of catalyst, and previously determined operational and performance parameters or target operational parameters, and providing the determined operational and / or performance parameters or target operational parameters of the chemical production plant including one or more catalytic reactors.

[0083] A further aspect of the present disclosure is a computer-implemented method for determining operational and / or performance parameters of a chemical production plant including one or more catalytic reactors, the method including accessing sensor data indicative of operating conditions present in the reactor system and catalyst-specific sensor data indicative of a catalyst currently being used in the reactor system, determining operational and / or performance parameters or target operational parameters of the chemical production plant using a data-driven model, where the data-driven model is parameterized according to a training data set including historical data, e.g., historical data from multiple production runs of the same type of catalyst and / or historical data from multiple production runs in multiple catalyst systems including the same type of catalyst, and previously determined operational and performance parameters or target operational parameters of the chemical production plant, and providing the determined operational and / or performance parameters or target operational parameters of the chemical production plant.

[0084] Another aspect of the present disclosure is a control system for controlling a chemical production plant including one or more catalytic reactors, comprising the computer system described above and a control unit configured to control actual and / or planned production operations in the chemical production plant based on the provided operational and / or performance parameters or target operational parameters of the chemical production plant.

[0085] The computer program may be stored and distributed on a suitable medium, such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless communication systems.

[0086] However, the computer program may also be presented over a network, such as the World Wide Web, and may also be downloaded into the working memory of a data processor from such a network.

[0087] According to a further exemplary embodiment of the invention, a data carrier or data storage medium for making a computer program available for download is provided, which computer program is arranged to perform a method according to one of the previously described embodiments of the invention.

[0088] It should be understood that the embodiments described herein are not mutually exclusive, and that one or more of the described embodiments can be combined in various ways, as would be understood by one of ordinary skill in the art.

[0089] A computer program for executing any of the methods of the present invention may be stored on a computer-readable storage medium (e.g., a non-transitory computer-readable storage medium). The computer-readable storage medium may be a floppy disk, a hard disk, a CD (compact disk), a DVD (digital versatile disk), a USB (universal serial bus) drive, a RAM (random access memory), a ROM (read only memory), an EPROM (erasable programmable read only memory). The computer-readable medium may be a data communication network, such as the Internet, that allows downloading of the program code. The methods, systems, and apparatuses described herein may be implemented as software in a digital signal processor, a DSP, a microcontroller, or any other side processor, or as hardware circuitry in an application specific integrated circuit, an ASIC, a CPLD, an FPGA, or other suitable device. The present invention may be implemented in digital electronic circuitry, computer hardware, firmware, software, or a combination thereof, for example in available hardware of conventional mobile devices, or in new hardware dedicated to processing the methods described herein, as described in more detail below. [Brief description of the drawings]

[0090] Exemplary embodiments of the present invention are illustrated in the accompanying drawings. It should be noted, however, that the accompanying drawings depict only certain embodiments of the present invention and therefore should not be considered as limiting its scope. The present invention may encompass other equally effective embodiments. [Figure 1] FIG. 1 illustrates an exemplary embodiment of a method for determining operating conditions for a planned production run, or for changes in a current production run, of a chemical production plant that includes at least one catalytic reactor. [Diagram 2] FIG. 1 illustrates a high-level workflow for pre-processing raw data from individual plants into a dataset ready for model training or prediction. [Diagram 3] FIG. 1 illustrates an exemplary workflow for identifying a selection of raw datasets as the basis for a model. [Figure 4] FIG. 1 illustrates an example implementation of a system for determining operating conditions for a planned production run or for changes in a current production run of a chemical production plant. Detailed Description of the Invention

[0091] The present disclosure provides a computer-implemented method for predicting short-term performance of a catalyst, including the effects of catalyst aging, and / or forecasting long-term performance in a chemical production plant including a catalytic reactor. The method includes using a mathematical model of the catalytic reactor, which is based on machine learning and does not include a priori information on reaction kinetics, and uses input parameters selected from sensor raw data, derived parameters, reactor operating parameters or data, plant metadata, and parameters indicative of catalyst aging. The method, system, computer program, and computer program product disclosed herein are further described with respect to a styrene production plant serving as an example. The method, system, computer program, and computer program product disclosed herein are applicable to other production plants having at least one catalytic reactor, particularly a fixed bed reactor.

[0092] In the case of a styrene production plant, a feed stream containing ethylbenzene is mixed with steam in a mixer. This mixed stream is fed to a catalytic reactor containing a potassium promoted iron oxide based catalyst for dehydrogenation to styrene monomer. The styrene production plant further includes temperature sensors, pressure sensors, flow sensors, etc. at various locations to monitor the operating conditions of the plant.

[0093] In the production of styrene monomer, ethylbenzene can be dehydrogenated in an adiabatic radial flow reactor. Ethylbenzene is mixed with steam in a specific ratio, known as the steam-to-oil (STO) ratio, to provide heat for the endothermic dehydrogenation process and to prevent reduction and coking of the potassium-promoted iron oxide-based catalyst used in the process. The reaction is usually carried out at elevated temperatures and below atmospheric pressure in a configuration with at least two reactors in series. Intermediate reheating can be used to compensate for the energy consumed in the reaction. The low pressure, dilution with steam, and high temperature favor the dehydrogenation of ethylbenzene, resulting in a higher equilibrium conversion.

[0094] Over the life of the catalyst, potassium, a coke gasification promoter, evaporates from the catalyst and is carried downstream of the catalyst bed to the cooling outlet, causing activity loss due to catalyst coking. To compensate for catalyst aging, plant operators increase the inlet temperature over the life of the catalyst to keep the ethylbenzene conversion constant. Higher temperatures adversely affect selectivity to styrene by promoting cracking and the production of by-products such as benzene and toluene. In addition, potassium deposited in the catalyst bed and fines generated from the catalyst cause an increase in the reactor inlet pressure, which is thermodynamically unfavorable. If the inlet temperature or pressure increases beyond the plant's operating limits, operation must be stopped and the catalyst replaced.

[0095] The performance of the catalyst and its aging rate depend on the reactor operating parameters, such as, for example, inlet temperature, STO ratio, inlet and outlet pressures, ethylbenzene flow rate, number of reactors used, etc. Optimizing the reactor operation can significantly improve the process economics. This has been attempted in the past by using different types of reactor models fitted to catalyst operating data derived from experimental and industrial reactors. These models were based on knowledge or assumptions about the reaction rates of the main and side reactions, mass, heat and impulse transport phenomena, adsorption / desorption of the various chemical species present in the system, rates of coke gasification and potassium loss, etc. The methods, systems, computer programs and computer program products described herein enable more robust and reliable process control.

[0096] FIG. 1 illustrates an exemplary embodiment of a method for determining operating conditions for a planned production run, or for changes in a current production run, of a chemical production plant that includes at least one catalytic reactor.

[0097] If the decision relates to a scheduled production run, operational data indicative of predefined operating conditions for the scheduled production run is received in a first step 10 via the communication interface. The predefined operating data may result from a previous production run. If the decision relates to a change in a current production run, measured operational data indicative of current operating conditions for the current production run is received in a first step 10 via the communication interface. Furthermore, at least one operational data point may be adjusted to include a desired operational value indicative of a change in the current operating conditions. The operational data may include sensor data measured via sensors installed in the chemical production plant, quantities derived directly or indirectly from said sensor data, analytical data measured on samples taken from the chemical production plant, and / or quantities derived directly or indirectly from said analytical data.

[0098] Additionally, catalyst aging indicators relating to the duration the catalyst has been used in a current or planned production run are received via the communications interface. The catalyst aging indicators may be based on a point in time, a period of time, amounts obtained from the time-dependent operating data, and / or amounts cumulatively derived from the time-dependent operating data. Additionally, plant metadata indicative of a physical plant layout can be received via the communications interface.

[0099] In a second step 12, the operational data and plant metadata may be pre-processed via a processing device prior to determining the at least one target operational parameter. Preferably, the pre-processing includes conversion to quantities that are independent of the physical plant layout.

[0100] In a third step 14, at least one target operating parameter of an operating condition of a planned production run or a change in a current production run is determined via the processing unit based on the operating data and the catalyst aging indicator using a data-driven model. The determination of the at least one target operating parameter may further be based on plant metadata. The data-driven model is parameterized according to a training data set. The training data set may be based on a historical data set including the operating data, the catalyst aging indicator, the at least one target operating parameter, and optionally the plant metadata. The historical data set may include data from multiple operations, multiple plants, and / or multiple catalyst batches. The determination of the at least one target operating parameter may be based on a short-term model that determines the target operating parameters at discrete time points, or a long-term model that determines the target operating parameters for a period of time.

[0101] For long term models, a time series of at least one target operating parameter, measured, predicted or derived for a current or previous production run up to a prediction start time point, may be received via the communications interface in step 10. Determination (14) via the processing unit of the at least one target operating parameter at one or more time points following the prediction start time point may be based on operational data including the desired operating values, the time series of the at least one target operating parameter, catalyst aging indicators, and optionally plant metadata using a data driven model. Preferably, the data driven model includes an inherent time dependency.

[0102] In a fourth step 16, at least one target operating parameter for the operating conditions of the planned production run or a change in the current production run may be provided via the communication interface.

[0103] In a fifth step 18, the determined target operating conditions may be provided to an optimization processor for optimizing operating conditions of a planned production run or a change in a current production run of the chemical production plant, where the determined target operating parameters for a plurality of operating conditions of the planned production run or a plurality of operating conditions of a change in a current production run may be received via a communication interface between the optimization processor and the processor, and based on the received target operating parameters, a minimum or maximum value of the target operating parameters or a minimum or maximum value of an optimization parameter derived from the target operating parameters may be determined via the optimization processor, and a minimum or maximum value indicative of an optimal operating condition for the planned production run or for the change in the current production run may be provided via the communication interface.

[0104] Pre-processing of driving data FIG. 2 illustrates an exemplary workflow for pre-processing data into a form suitable for model training, short-term model prediction, or long-term model forecasting.

[0105] In a first step, measured operational data is received, and optionally subsequently pre-processed. If data is being prepared for learning, such operational data may include historical data sets from multiple runs, multiple plants, and / or multiple catalyst batches. If data is being prepared for prediction or forecasting, such operational data may include measured operational data indicative of current operating conditions for a current production run. The operational data preferably includes sensor data measured through sensors installed in the styrene production plant, and / or analytical data measured on samples taken from the styrene production plant.

[0106] Each production plant has many sensors available, usually hundreds or even thousands, providing raw data at their own sampling rate. In addition, analytical data, e.g. results from gas chromatography, are available for the specific time that a sample was taken from the plant. How frequently these data are available varies from plant to plant, but typical frequencies range from once a day to once a week.

[0107] In a first pre-processing step, driving parameters can be selected from the driving data that form the input parameters of the data-driven model. These input parameters may be derived from raw parameters such as sensor data or analytical data. An exemplary process on how to select these raw parameters from all available parameters when setting up the data-driven model is described in further detail below.

[0108] In a second pre-processing step, the analytical and selected sensors may be combined based on their timestamps and specifically pre-processed to a common time scale. In many plants, only daily collections of raw sensor data may be available, rather than high-frequency raw data, which is also the typical frequency at which analytical data is available, and a daily average may be used to bring all data to the same time base. Other fusion techniques are equally applicable and well known, such as interpolation of daily data (both analytical and sensor), and more or less frequent sampling of the subjects to create, for example, hourly data.

[0109] Additionally, plant metadata may be received, such as catalyst activity and number of reactors, which may be added to the dataset as numeric or categorical variables to complete the input parameter set for the model, including the derived parameters.

[0110] In a third pre-processing step, the received and selected operating data can be filtered and smoothed. Here, for example, time points and periods of maintenance intervals, start-up phases, irregularities, outliers, etc. can be identified and optionally filtered. To achieve this, many options are known to the skilled person, and any combination of such methods can be used. For example, this procedure can include applying feasible absolute thresholds based on knowledge of the catalyst area, such as minimum reactor temperature, maximum steam-to-oil ratio, maximum pressure after the reactor, etc.; identifying outliers by comparing each value or set of values ​​with the distribution of all other values ​​of the same parameter or set of parameters available from the respective production plant, for example using absolute thresholds based on a distribution metric such as six times the interquartile range (a conservative threshold used in the example application), or alternatively identifying irregularities based on the estimated likelihood of a data point arising from the overall distribution and / or based on large jumps in a parameter compared to the monthly coefficient of variation of this parameter.

[0111] In a fourth pre-processing step, missing data points on a common time scale can be detected and replaced with statistically determined values. Such potentially missing parameters can be imputed. In particular, if analytical data are sampled less frequently than the selected time base, these data can be imputed. Imputation can be determined from various methods such as simple mean value imputation, forward or backward filling, weighted average, or estimates from a Kalman filter or equivalent estimator. The same method may be applied to replace outliers. The start of a run may be identified based on expert-defined criteria, for example, in the styrene catalyst example, as the first data point with a space velocity greater than 0.2 / h per hour (>0.2 / h). All derived parameters are then calculated, including cumulative plant production, etc.

[0112] Depending on the nature of the derived parameters, especially the cumulative parameters, identified start-up phases, downtimes, etc. can be removed from the data set before or after calculating the respective derived parameters. In an exemplary embodiment, the cumulative production since start-up on which the implemented aging indicators are based is calculated before removing start-up phases from the data set, since they contribute to catalyst aging even though they are not part of the operating conditions covered by the model.

[0113] One further pre-processing step, not shown in Figure 2, could be the conversion of data from different plants to common units, preferably directly before the workflow of Figure 2 begins, but could also be done before or after any step of the process.

[0114] At this point, the data can be used for learning, forecasting with long-term models, or prediction with short-term models. However, because the filtering procedure may result in some gaps in the data, and because the time scale for forecasting is usually much longer (e.g. months) compared to the available time bases (e.g. days), an additional aggregation step (e.g. weekly mean or median aggregation) may be performed on the data before using it in the long-term models.

[0115] Parameter Selection Process The operating data or raw parameters of interest can be determined based on the workflow outlined in Figure 3. The parameters available to the catalyst expert are based on a combination of raw sensor data, analytical data, and typical operating parameters derived therefrom (e.g., conversion, space velocity, steam-to-oil ratio, selectivity, etc.). In a first step, when using data from multiple plants, it is possible to select only parameters available from all plants of interest (this criterion is not used when modeling only one plant) and to filter out redundant parameters (e.g., rescaled parameters) and zero-variance parameters.

[0116] A correlation matrix for the remaining parameters can then be calculated and a clustering, for example a hierarchical clustering, of the parameters can be performed to identify pairs or clusters of parameters with similar information. Instead of a clustering algorithm, a simple filter for high (anti)correlation values ​​can be applied, for example (>0.90 or <-0.90) or (>0.95 or <-0.95). However, clustering the correlation values ​​with all other parameters makes it even possible to identify parameters that have a low direct correlation but very similar correlation values ​​with all other parameters. From each cluster of two or more parameters, only one parameter can be retained, based on certain selection criteria such as: the parameters must be available in all plants; typical parameters used routinely by the operators must be retained (this is an exception that could lead to retaining one or more parameters) and the parameters must represent interpretable quantities, which is also the reason not to reduce the dimensionality by principal component transformation or similar methods that generate features that cannot be directly interpreted by the operator.

[0117] Based on these criteria, the parameter set can be iteratively reduced by adjusting the cluster threshold as necessary. The number of parameters can be iteratively reduced until a parameter set (1, ..., M) remains, where only low correlations between parameters survive. The remaining relatively high correlation values ​​come from parameters that are important to the operator and therefore may not be removed.

[0118] Once the parameter set is reduced, the raw sensor and analytical data required to derive all of these remaining parameters is identified.

[0119] Normalization of parameters and plant metadata Different plants have different production levels and their typical operating conditions often show systematic differences. To take such differences between different plants into account, one of the two strategies shown below, or a combination of the two strategies, can be applied:

[0120] In a further pre-processing step, normalization may occur. In one embodiment, plant metadata may be received that indicates a physical plant layout. Such plant metadata may include reactor layout, such as number of reactors, active catalyst amount, reactor type, dimensions, or a combination thereof.

[0121] Normalization of the operational data can be done to make the model input parameters (except for the plant metadata) as independent as possible from the specific plant layout. Many parameters, such as steam to oil ratio, conversion, selectivity, etc., do not need to be normalized. Furthermore, parameters such as liquid hourly space velocity (LHSV) are inherently normalized by catalyst volume. Additionally or alternatively, the amount of total product in a catalyst production run per catalyst, or cumulative total production, may be normalized by active catalyst volume, as this is a more comparable measure of the "aging" of each unit volume of catalyst than flow time (which does not capture variations in production levels) or unnormalized cumulative production (which has different meanings for different size reactors).

[0122] Additionally or alternatively, the pressure may be normalized to an initial value at start-up, e.g., the median value for the first 90 days of flow, to focus on the effects of aging rather than plant-to-plant differences, which may be more pronounced in some cases. Additionally or alternatively, the pressure drop in each unit, such as a reactor or heat exchanger, may be normalized to the space velocity or total flow rate, since it is known that this varies with the superficial linear velocity of the mixed gas. Additionally or alternatively, the average inlet temperature of multiple, e.g., two or three, reactors may be determined.

[0123] Ultimately, methods for making operational data, analytical data, and quantities derived therefrom comparable between plants exist beyond the examples given in this section; there may be other useful reactor / plant metadata, particularly relating to reactor geometry, that could be used in a similar fashion.

[0124] The reason for finding such normalized parameters is that using data from multiple different plants offers a number of important advantages: 1) the parameter space covered by the collected data set is much larger than that of the individual plants, which usually operate around a relatively narrow set of operating parameters that operators are usually familiar with, from which erroneous deviations can result in significant financial losses. Thus, models trained on such collected data sets may provide predictions outside the operating range of a particular plant, since they contain information that is not available from its own historical data. 2) Due to the long life span of the exemplary styrene catalyst (2-3 years), certain catalysts can only be utilized 1-4 times per plant, severely limiting the number of deactivation processes that can be observed per plant (each run can only independently observe the deactivation of the catalyst once). The collected data set allows a larger number of deactivation processes to be included in the training data. 3) Combined with a parameter selection process that emphasizes the interpretability and availability of commonly used operating parameters, the learned models can be applied to plants where no prior real data is available (a situation often encountered when technical proposals for new plants are required), which would not be possible if the models were strictly linked to the specific set of sensors available at a particular plant.

[0125] Short-term model Any regression model can be used for short-term prediction of catalyst behavior, and those skilled in the art are aware of various typical model candidates. Depending on the nature of the data of a chemical plant, typically, several potential plants with a small number of independent runs and (even after normalizing many parameters) relatively simple in nature and less flexible models (with inherent biases) can be selected. In scenarios with a small number of independent runs, flexible models such as random forest regression may fit the training data set better but extrapolate poorly to new data. Therefore, depending on the number of runs available, an appropriate type or combination of regression-based models can be selected.

[0126] One possible model is an ensemble of linear models trained on a subset of the training data, for example to predict the mean reactor inlet temperature and the selectivity of the reaction. There are two advantages to using such an ensemble: first, using a model ensemble for prediction can lead to more accurate predictions, for example by taking the median prediction of all models (see, e.g., Ensemble Methods - Foundations and Algorithms, Zhi-Hua Zhou; CRC Press 2012); second, training an ensemble in this way provides an estimate of the model uncertainty, by using a range of prediction values ​​(or the 10% and 90% percentiles, or a range of other prediction quantiles, ...).

[0127] When training data from multiple catalyst batches, multiple runs, and / or multiple plants are used, the training and test data sets may be split by individual runs. Additionally or alternatively, the training and validation data sets may be split by individual plants or catalyst batches. For example, a random set of about 75% of the training set may be selected and parameters may optionally be normalized (zero mean and unit variance) during preprocessing. The remaining data may be used as a validation data set to test the trained model.

[0128] Long-term model Numerous mathematical models can be used for time series forecasting, ranging from autoregressive models to recurrent neural networks. Requirements for a model to apply to the present exemplary problem include: 1) applicability to multivariate time series, i.e., the prediction of multiple endogenous (uncontrolled) parameters that are subject to expected long-term trends, e.g., pressure, inlet temperature, or selectivity, and 2) integration of exogenous (controlled) parameters, i.e., integration of influential parameters such as steam-to-oil ratio, LHSV, target conversion, etc., that are known or externally controlled and therefore do not require prediction by the model.

[0129] A preferred embodiment of the method includes a mathematical model that allows for regularization to avoid overfitting, which is likely to occur when multiple time lags are included in the mathematical model. In one embodiment for generating the mathematical model used in the method of the present disclosure, an ensemble of mathematical models is implemented that is based primarily on a combination of penalized linear models and penalized vector autoregressive models with exogenous variables (VARX). For an overview of penalized VAR(X) models, including different structural regularization methods, see, for example, [arXiv:1508.07497v1 (Nicholson et al. 2018, VARX-L: Structural Regularization for Large Vector Autoregressions with Exogenous Variables)] and references therein.

[0130] The training / test split is done in the same way as for the short-term models during runs: All candidate models can be trained to predict the endogenous variables at time t based on the history of the exogenous variables at time t and the endogenous variables up to a maximum time lag m (tm,...t-1). Applying such a model repeatedly step by step, using at each new step the predictions of the endogenous variables as input for the next step, allows predictions to be made for any number of steps ahead.

[0131] The training procedure may be performed using leave-one-out cross-validation at the plant level, where for N plants in the dataset, N sets of training data are generated (consisting of data from all other plants), whose parameters can optionally be standardized to have mean 0 and a comparable range of values, and the trained model can be evaluated on the plants left out (validation set). Finally, the model hyperparameters (e.g., regularization parameters) can be selected that give the most robust performance on the N validation sets, measured, for example, by the average root mean square of one-step-ahead predictions, and the model can be trained on the full training set.

[0132] Since the mathematical model used in the disclosed method is based on machine learning, it must be trained with historical data from at least one production run in a chemical production plant including at least one catalytic reactor before it can be used to predict short-term or long-term performance including degradation of the catalyst used in such reactor system. In one embodiment of the disclosed method, the model is trained with historical data from multiple production runs of the same type of catalyst. In a further embodiment of the disclosed method, the model is trained with historical data of production runs in multiple reactor systems including the same type of catalyst. In both of these embodiments, the historical data may be provided from different manufacturing batches of the same type of catalyst. In one embodiment of the disclosed method, the operational data of the multiple production runs and the potential catalyst aging indicators are normalized, for example as described in the above section, before being used to train the model.

[0133] It has been found that using operational data from multiple production runs, whether from one and the same plant or from different plants, improves the quality of the prediction and expands the range of operational parameters covered by the prediction or forecast model. The inclusion of data originating from multiple runs and plants further improves the versatility of the prediction or forecast for application to production plants for which no data was available during model training. However, since different catalysts vary widely in their catalytic properties (reaction rates) and morphological properties (transport properties), it is preferred that all data used to train the model are provided from plants that used the same catalyst formulation. This includes the use of the same type of catalyst provided from different production batches or different delivery dates.

[0134] In one embodiment, the input parameters for the model are selected from sensor data available from the reactor system, analytical data such as gas chromatography (GC) analysis, and derived parameters as described in the respective sections above.

[0135] In one embodiment, the operating data or reactor operating parameters include the inlet and outlet temperatures of each reactor, the inlet and outlet pressures of each reactor, and the composition of the reaction mixture at the inlet and outlet of each reactor.

[0136] In one embodiment, the operating data or input parameters for the model include steam to oil (STO) ratio, liquid hourly space velocity (LHSV), total styrene production normalized by catalyst loading, target ethylbenzene conversion, styrene selectivity, average inlet temperature, normalized pressure after the last reactor, normalized pressure drop across the reactors, temperature loss across the reactors, and normalized deviation of temperature loss compared to expected value (calculated based on conversion).

[0137] An advantage of the parameter sets used in the disclosed methods or systems is that many of the parameters correspond to or can be easily interpreted in the context of actual operating parameters and therefore remain interpretable. Other approaches to reduce the dimensionality of the problem (e.g., PCA or RFA) can be useful for the modeling process and predictive accuracy, but often result in parameters that lack interpretability.

[0138] The mathematical model used in the disclosed method is capable of time series forecasting and full coverage of the catalyst deactivation process in the entire production run of the reactor system, allowing the plant operator to improve and optimize the daily operating policy depending on the energy costs, market supply of raw materials or demand for the plant products, and other limitations that may occur within the plant, such as the shutdown of various plant parts or equipment.

[0139] The subject matter of the present disclosure is further described and illustrated in the following examples. EXAMPLES

[0140] Generate a dataset for training a mathematical model To build the mathematical model, data from several production runs of the same type of catalyst (BASF S6-42) were used. The dataset covers information on 11 industrial plants with two or three reactors and one to four production runs at each plant. For each individual plant, a complete set of parameters (including analytical and sensor data) was collected. Sensor data were typically available with a daily resolution, whereas analytical data were available with a daily to weekly resolution.

[0141] The process of model parameter selection has been described above (Figure 3).

[0142] Table 1 shows the operational data from the raw sensors and analytical parameters selected at each time point to derive all parameters relevant to the model. These unit and format choices are only a few examples of what can be used; temperature could be specified in degrees Fahrenheit, pressure in mmHg, and alternative date and time formats could be used.

[0143] [Table 1]

[0144] The number of reactors and the total catalyst amount were additionally used as metadata. Table 2 below lists the parameter sets (derived parameters and plant metadata) used to train the different models.

[0145] [Table 2]

[0146] Example use: Short-term models

[0147] Model Development As mentioned above, an ensemble of 50 linear regression models is trained on a subset of the training data set, with each subset split by run. Each of the 50 reduced training sets is randomly selected to represent up to seventy-four percent (~74%) of the available runs, improving predictions and providing estimates of the local uncertainty of the model predictions (given that each model is trained on a different subset of operating conditions, etc.). It is important to note that the subsets of data need to be split across runs (or across plants) rather than a random sampling of training data points. Otherwise, all models will be trained on roughly the same distribution of operating conditions, resulting in little inter-model variability.

[0148] Specifically, in this example, the parameters "Temperature" (average reactor inlet temperature) and "Selectivity" (of the reaction to the desired styrene) were predicted using all other parameters in Table 2.

[0149] In this specific implementation, the reaction temperature itself is one of the major influencing factors on the observed selectivity. Therefore, the prediction of both parameters is done in two steps. First, an ensemble of 50 models is trained to predict temperature based on all parameters in Table 2 (except selectivity and temperature). A second ensemble of 50 models is trained to predict selectivity using all parameters in Table 2 (except only "selectivity") as input. This second ensemble is then used to predict selectivity using the predicted temperature from the first model ensemble as one of the input parameters.

[0150] In practical use of the developed model, the workflow of predicting both parameters consecutively was implemented in a single prediction function that receives input parameters and predicts both temperature and selectivity, and therefore the combination of both model ensembles can be externally viewed as a single instance of the "short-term" model in all applications of the model.

[0151] In addition to directly using the parameters in Table 2, interaction terms were also considered and finally the "steam to oil ratio-conversion" interaction was integrated into the short-term model for predicting temperature. Such interactions or higher order polynomial terms, such as a quadratic term in temperature for selectivity prediction, can be easily implemented in the predictive model using a statistical programming language without the need to extend the basic parameter set provided in the model.

[0152] Which type of model to choose, whether to use ensemble techniques, and which higher-order terms, transformations, and interactions of input parameters to use depend on the specific question and dataset and are typical model development procedures of a data scientist.

[0153] Short-term model usage example 1 One exemplary use case is that a plant operator may want to use the short-term model to estimate the reactor temperature required to achieve a desired target conversion at a particular operating condition not previously used in the plant. For a styrene production plant, such operating conditions could be an increase in steam-to-oil ratio, a decrease in LHSV, or a change in feed composition. Ideally, the starting point for this prediction should be the current plant state, including the normalized aging parameter "totalProduction".

[0154] Model Implementation FIG. 4 illustrates a client-server setup of a production monitoring and / or control system including a client side with a processing device including a user application and a server side with a processing device for a service provider. The client side and the server side may be communicatively coupled via a communication interface, e.g., wired or wireless. The client side may include a display device. Preferably, the user application on the client side is configured to receive and display the operating data, the desired values, the catalyst aging indicators, the plant metadata, the target operating parameters, or the determined operating conditions. More preferably, the client side user application is configured to receive the target operating parameters or the determined operating conditions and control a current or scheduled production operation of the chemical production plant based on the target operating parameters or the determined operating conditions. The client side user application may be integrated as part of a process monitoring and / or control system of the chemical plant.

[0155] Raw sensor and analytical data are recorded in the production plant together with the necessary plant metadata since the start of the production run. These data are pre-processed according to a data pre-processing workflow, for example as shown in Fig. 2, in particular using the same workflow (filter steps, imputation, threshold, collection, ...) that was used to prepare the training dataset of the respective model. Preferably, this pre-processing is implemented by the party that developed the model and can be provided to the user application by direct integration of the pre-processing functionality or via an application programming interface (API).

[0156] After the data is pre-processed, the user can adjust the operating parameters of interest to provide desired operating values ​​indicative of changes in the current operating conditions. The adjusted set of input parameters (also called operating scenarios, or forecast scenarios) is then transferred to a forecast function. This forecast function, for example implemented locally in the user application or addressable via an API, performs all operations as described above. The results are reported to the user application, for example to compare and select different scenarios. The user can further associate the input parameters and / or forecast parameters with, for example, costs or other quantities that affect the decision process regarding the operating parameters to use in the plant.

[0157] Figure 4 visualizes an exemplary implementation concept. Raw data is collected automatically or manually at a production plant and transferred via an API to a service provider, which processes the raw data into the correct format according to the workflow as described above. The transformed data set is provided to a user, such as a plant operator, who can define different scenarios based on current values. These scenarios can be forwarded to the same or a second API running a model, which can provide the user with the corresponding forecasts or predictions.

[0158] Short-term model use case 2 In the second use case, an expert may want to provide an estimate of the reactor temperature and the corresponding selectivity development before the catalyst is installed in the plant. This scenario arises often, for example, during the preparation of a technology proposal, which provides the customer with one or more what-if operating scenarios and their implications before the catalyst purchase decision. Thus, an accurate prediction is desirable for a plant where no data is provided for model learning.

[0159] By making some simplifying assumptions, for example estimating the "aging" parameter totalProduction based on hypothetical operating data, which usually includes constant operating parameters for the entire run, replacing the reactor temperature losses with estimates based on operating conditions, etc., the short-term model described above can be used for this purpose, provided a set of typical operating parameters is provided by the customer. What the model lends itself to is a direct consequence of the selection criteria used in the parameter selection workflow (Figure 3), and in particular, all parameters that are usually used by the operator to monitor and control the plant can be kept as inputs to the model.

[0160] Presenting the model to the user The user application interface can include blocks of input parameters to the model that can be entered automatically or manually and used to interactively define operating scenarios; blocks of model outputs for the defined scenarios (e.g., average reactor inlet temperature and selectivity with prediction ranges); and, optionally as further outputs, local response of the model predictions to parameters of interest (e.g., predicted reactor inlet temperature at target conversion ranges).

[0161] The input parameter blocks can be filled in automatically if current plant operating data is available (e.g., Use Case 1) or can be completely manually defined if no actual plant data is available (e.g., Use Case 2). Based on the user-tuned input parameters, the parameters predicted by the model (in this example, reactor inlet temperature and styrene selectivity) can be displayed to the user, for example, in text (for single point predictions) or graphical form (local model response to variation of a single parameter). These are provided merely as examples, as there are many more use cases in which such models are used, implemented, and presented to the model user.

[0162] Exemplary application: Long-term model

[0163] Model Development Starting from the processed dataset used to develop the short-term model (Table 2), several additional steps are performed to prepare the training dataset for the long-term model: First, a smaller number of parameters is selected from the list, and then each parameter is collected on a weekly basis.

[0164] The input parameters of the model may be constant (reactors, CatVolume), controlled (conversion, SOR, LHSV), calculated based on these constant or controlled values ​​(totalProduction), or uncontrolled for the operating scenario of interest (temperature, selectivity, pressureOut, deltaP, deltaT, dTdev_norm).

[0165] The latter parameters are labeled endogenous or uncontrolled parameters, while the former parameters (which are known over the entire forecast horizon, assuming the operating scenario is carried out as planned) are labeled exogenous or controlled parameters. Table 3 provides an overview of the parameters used to develop an exemplary forecasting model and the allocations to both types of parameters.

[0166] [Table 3]

[0167] In developing the example model, several different VARX-type candidate models mentioned above are trained, which differ, for example, in their regularization methods (elastic net, ridge) and maximum time lags (4 to 10 weeks).

[0168] For all of these different model combinations, the training procedure was performed in the same way using leave-one-out cross-validation at the plant level. All trained models were then candidates for the final ensemble that included the long-term model.

[0169] Although all learning was done for one step ahead of prediction (it can be applied iteratively and can forecast any number of steps), the final selection of the model requires that it also performs well for longer forecasts. To evaluate this, forecasts were performed from multiple time points for each available drive in the training and test datasets, using the actual exogenous variables of those drives as driving scenarios, and determining each step ahead of the error distribution of the forecasts.

[0170] Models were selected that had clear long-term biases and narrow error distributions throughout the entire range (even though the error distributions for long forecasts are expected to be moderately wide) for reactor inlet temperature and selectivity in both the training and test data. Model selection was based on the error distributions of forecasts performed on the training and test sets from 30 candidate models.

[0171] In this example development of a long-term model for styrene catalysis, the final models selected were three two-stage models with time lags up to 10 weeks trained with local gradients using different types of regularization. For the ensemble forecast, all models are run repeatedly independent of each other and only the complete individual forecasts are averaged to provide the ensemble forecast. In the presented example, this was a more efficient implementation compared to collecting forecasts after each individual step, which is one possible alternative implementation as an alternative to the ensemble forecast.

[0172] As in the single-model case, the final ensemble is evaluated on the test data to obtain an N-day-ahead error distribution, which serves as the error estimate for the forecast in the illustrated implementation.

[0173] As with the short-term example above, the ensemble forecasts and expected error distributions can be viewed as a single instance of the “long-term” model when the user applies the model, regardless of the details of the underlying procedure (ensemble averaging, multi-step models, etc.).

[0174] Long-term model usage example In an exemplary use case, a plant operator may want to estimate the remaining life of a catalyst depending on different scenarios of plant operation and plant-specific limits on operating parameters. Such scenarios may include changes in LHSV, steam-to-oil ratio, or target conversion levels; while there are many motivations for considering different scenarios, a typical question is whether it is ultimately worthwhile to gain catalyst life by operating the plant at a lower production level during a period when styrene prices are low.

[0175] The end of catalyst life may depend on many conditions that are locally defined in each production plant, but one limit of catalyst life is in all cases the reactor temperature or pressure, which must not exceed a plant-specific threshold. With predictions from the long-term model, the end of catalyst life can be estimated based on this threshold and compared to different user-defined scenarios.

[0176] Model Implementation The model can be implemented conceptually similar to that described in FIG. 4 with some minor modifications. For example, the processing function implemented locally or via an API needs to be adapted to provide the data format used for training the long-term model, in this particular example after the procedure performed on the short-term model data, the weekly collection step. For example, the forecasting function implemented locally or via an API now not only receives a set of driving conditions to forecast a single point, but also receives driving scenario values ​​for the exogenous parameters as well, instead of the delayed endogenous data required by the model. Similarly, the model output is a full forecast of all or only a portion of the endogenous variables of the model.

[0177] Presenting the model to the user In implementing a long-term model, the endogenous parameters used in the model will need to span at least the last L weeks, where L is the maximum time lag used for any model element included in the "long-term" model. This data may be obtained automatically from raw plant data or data in a suitable format as described in the data pre-processing and preparation section, or it may be entered or uploaded manually.

[0178] The input blocks of control parameters may be used by the operator to generate operating scenarios over extended periods of time; for example, in the example application, the steam-to-oil ratio, target conversion, and LHSV can be pre-planned up to three independent segments to simulate future changes as well. These scenarios can be much more complex than those described here. The planned operating scenarios can be displayed with the history of these parameters, preferably to facilitate continuity of some of the control. The control parameters selected for user manipulation may consist of any subset of all exogenous parameters used in the model. In this particular implementation, the other three exogenous parameters are constant or can be directly derived from other inputs or time.

[0179] The forecasts provided by the long-term model can be displayed in tabular form and subject to data export and further analysis, for example based on predicted trends, or the results can be displayed graphically to the user for visual inspection of one or more different projected scenarios.

[0180] The interface may include a part where selected control (exogenous) parameters can be defined for a user-defined operating scenario (e.g., target conversion, LHSV, steam-to-oil ratio; plots showing actual data up to the start of the prediction and for future operating scenarios). In another part of the interface, selected endogenous parameters can be presented to the user (e.g., reactor inlet temperature and styrene selectivity; the start of the prediction is indicated by a horizontal dashed line, everything from that point on is predicted by the model, including estimates of the prediction error). Additionally, predictions of all or some of the endogenous parameters can be displayed in a data frame that can be exported by the user.

Claims

1. 1. A system for determining operating conditions of a chemical production plant including at least one catalytic reactor, the system comprising: a communications interface; and a processing device in communication with the communications interface; (a) For a scheduled production run, the system: receiving (10) operational data indicative of predefined operating conditions for the planned production run via said communications interface; receiving (10) via said communications interface a catalyst aging indicator relating to the period of time that the catalyst has been used in said scheduled production run; via said processing device, determining (14) at least one target operating parameter for an operating condition of said planned production operation based on said operating data and said catalyst aging indicator by using a data-driven model, said data-driven model being parameterized according to a training data set, said training data set being based on a historical data set having said operating data, said catalyst aging indicator, and at least one target operating parameter; providing (16) said at least one target operating parameter for an operating condition of said scheduled production run via said communication interface; or (b) for a change in a current production run, the system: receiving (10) measured operating data indicative of current operating conditions for a current production run via a communication interface, where the operating data includes data indicative of operating conditions at at least one point in time, the data including desired operating values ​​indicative of a change from currently set operating conditions; receiving (10) via said communications interface a catalyst aging indicator associated with the period during which a catalyst has been used in said current production run; via a processing device, determining (14) at least one target operating parameter for a change in operating conditions in the current production operation based on the operating data and the catalyst aging indicator by using a data-driven model, the data-driven model being parameterized according to a training data set, the training data set being based on a historical data set having the operating data, the catalyst aging indicator, and at least one target operating parameter; providing (16) via said communication interface said at least one target operating parameter for a changed operating condition in said current production operation; A system configured as follows: The system comprises: receiving (10) plant metadata indicative of a physical plant layout via said communication interface; via the processing device, determining (14) at least one target operational parameter further based on the plant metadata by using a data-driven model, where the training data set is based on a historical data set further including plant metadata; the system is configured to pre-process (12) the operational data and the plant metadata via the processing device prior to the determination of the at least one target operational parameter, the pre-processing including conversion to quantities independent of the physical plant layout; The pretreatment is - pre-processing selected data from the received operational data to a common time scale; and - detecting missing data points in said common time scale and replacing them with statistically determined values; Including, the system.

2. 2. The system of claim 1, wherein the operational data includes sensor data measured via sensors installed in the chemical production plant, quantities derived directly or indirectly from such sensor data, analytical data measured on samples taken from the chemical production plant, and / or quantities derived directly or indirectly from such analytical data.

3. The system of claim 1 or 2, wherein the historical data set includes data from multiple runs, multiple plants, and / or multiple catalyst batches.

4. The system of any one of claims 1 to 3, wherein the catalyst aging indicator is based on a point in time, a period of time, an amount derived from time-dependent operating data, and / or an amount cumulatively derived from time-dependent operating data.

5. 5. The system of claim 1, wherein the system, via the processing device, determines (14) at least one target operating parameter based on a short-term model that determines a target operating parameter for a discrete time point or a long-term model that determines a target operating parameter for a period of time.

6. The system of any one of claims 1 to 5, wherein the system is configured to receive (10) a time series of the at least one target operating parameter measured, predicted, or derived during a current production run or from a previous production run to a predicted time point via the communication interface, and the system is further configured to determine (14) via the processing device at least one target operating parameter for one or more time points subsequent to the predicted time point based on the operating data, the time series of the at least one target operating parameter, and the catalyst aging indicator using the data-driven model having an inherent time dependency.

7. A system for optimizing operating conditions of a chemical production plant, the system comprising a system according to any one of claims 1 to 6 and an optimization processing device in communication with the communication interface, the optimization processing device comprising: receiving (18) via said communication interface determined target operating parameters for one or more operating conditions determined for either (a) a planned production run or (b) a change to a current production run; determining (18) via the optimization processor, based on the received target operating parameters for each operating condition, a maximum or minimum value of the target operating parameters or a maximum or minimum value of an optimization parameter derived from the target operating parameters; - providing (18) via said communication interface a maximum or minimum value indicative of optimal operating conditions of either (a) a planned production run or (b) a current production run.

8. A production monitoring and / or control system comprising a communication interface communicatively coupled to a system for determining operating conditions according to any one of claims 1 to 6 or a system for optimizing operating conditions according to claim 7.

9. 9. The production monitoring and / or control system of claim 8, comprising a display device configured to receive and display the determined operating conditions, or comprising a control unit configured to receive the determined operating conditions and control current or scheduled production operations in the chemical production plant based on the determined operating conditions.

10. 1. A computer-implemented method for determining operating conditions for a chemical production plant including at least one catalytic reactor, the method comprising the steps of: (a) For a scheduled production run, the method includes the steps of: receiving (10) operational data indicative of predefined operating conditions for a planned production run via a communication interface, - receiving (10) via said communications interface a catalyst ageing indicator associated with the period during which the catalyst has been used in said scheduled production run; determining (14) via a processing device at least one target operating parameter for operating conditions of the planned production operation based on the operating data and the catalyst aging indicator by using a data-driven model, the data-driven model being parameterized according to a training data set, the training data set being based on a historical data set having the operating data, the catalyst aging indicator and at least one target operating parameter; providing (16) via said communication interface said at least one target operating parameter for an operating condition of said scheduled production run; or (b) for a change in a current production run, the method comprising the steps of: - receiving (10) via a communication interface measured operating data indicative of current operating conditions for a current production run, said operating data including data indicative of operating conditions at at least one point in time, said data including desired operating values ​​indicative of a change from the currently set operating conditions; - receiving (10) via said communications interface a catalyst ageing indicator associated with the period of time that a catalyst has been used in said current production run; determining (14) via the processing device at least one target operating parameter for a change in operating conditions in a current production operation based on the operating data and the catalyst aging indicator by using a data-driven model, the data-driven model being parameterized according to a training data set, the training data set being based on a historical data set having the operating data, the catalyst aging indicator and at least one target operating parameter; - providing (16) said at least one target operating parameter for a changed operating condition in a current production run via said communication interface, Furthermore, the method comprises the steps of: - receiving (10) via said communication interface plant metadata indicative of a physical plant layout; and determining (14), via the processing device, at least one target operating parameter further based on the plant metadata by using a data-driven model, where the training data set is based on a historical data set further including plant metadata; Including, and pre-processing (12) the operational data and the plant metadata via the processing device prior to determining the at least one target operational parameter, the pre-processing including converting the operational data and the plant metadata into quantities that are independent of the physical plant layout; The pretreatment is pre-processing selected data from the received operational data to a common time scale; and detecting missing data points in said common time scale and replacing them with statistically determined values.

11. The method of claim 10 , wherein the historical data set includes data from multiple runs, multiple plants, and / or multiple catalyst batches.

12. 12. The method of claim 10 or 11, wherein the determination of the at least one target operating parameter is based on a short-term model that determines a target operating parameter for a discrete time point or a long-term model that determines a target operating parameter for a period of time.

13. 13. The method of any one of claims 10 to 12, wherein a time series of the at least one target operating parameter measured, predicted, or derived during a current production run or from a previous production run to a prediction time point is received via the communications interface (10), and the data-driven model is used to determine (14) at least one target operating parameter for one or more time points subsequent to the prediction time point based on the operating data, the time series of the at least one target operating parameter, and the catalyst aging indicator via the processing device, wherein the data-driven model includes an inherent time dependency.

14. 1. A method for optimizing operating conditions of a planned production run or changes in an existing production run of a chemical production plant, the method comprising the steps of: receiving (18) via the communication interface determined target operating parameters for one or more operating conditions determined by the method of claim 10 for either (a) a planned production run or (b) a change in a current production run; - determining (18) via an optimization processor, based on the received target operating parameters for each operating condition, a maximum or minimum value of said target operating parameters or a maximum or minimum value of an optimization parameter derived from said target operating parameters; - providing (18) via said communication interface maximum or minimum values ​​indicative of optimal operating conditions of either (a) a planned production run or (b) a current production run; The method according to claim 1,

15. A method for monitoring and / or controlling a chemical production plant, the method comprising the step of performing a method for determining operating conditions according to any one of claims 10 to 13 or a method for optimizing operating conditions according to claim 14.

16. A computer program which, when read and executed by a processing device, performs the method of any one of claims 10 to 15.

17. A computer readable non-volatile storage medium having computer readable instructions which, when read and executed by a processor, perform the method of any one of claims 10 to 15.

18. a catalyst; A catalyst system comprising a computer program according to claim 16 and an associated catalyst type identifier, A memory of the catalyst system stores a plurality of data-driven models including a catalyst type identifier; A catalyst system in which, in processing executed by the computer program, the data-driven model including the received catalyst type identifier is selected based on the received catalyst type identifier.

19. a catalyst; A catalyst system comprising a computer-readable non-volatile storage medium according to claim 17 and an associated catalyst type identifier, The computer-readable non-volatile storage medium stores a plurality of data-driven models including catalyst type identifiers; A catalyst system in which, during processing executed by computer-readable instructions contained in the computer-readable non-volatile storage medium, the data-driven model including the catalyst type identifier is selected based on a received catalyst type identifier.

Images