Optimization of Energy Consumption in Chemical Production Based on Descriptive Analytics and Neural Network Modeling

: Improving the energy efﬁciency of chemical industries and increasing their environmental friendliness requires an assessment of the parameters of consumption and losses of energy resources. The aim of the study is to develop and test a method for solving the problem of optimizing the use of energy resources in chemical production based on the methodology of descriptive statistics and training of neural networks. Research methods: graphic and tabular tools for descriptive data analysis to study the dynamics of the structure of energy carriers and determine possible reserves for reducing their consumption; correlation analysis with the construction of scatter diagrams to identify the dependences of the range of limit values of electricity consumption on the average rate of energy consumption; a method for training neural networks to predict the optimal values of energy consumption; methods of mathematical optimization and standardization. The authors analyzed the trends in the energy intensity of chemical industries with an assessment of the degree of transformation of the structure of the energy portfolio and possible reserves for reducing the speciﬁc weight of electrical and thermal energy; determined the dynamics of energy losses at Russian industrial enterprises; established the correlation dependence of the range of limiting values of power consumption on the average rate of power consumption; determined the optimal limiting limits of the norms for the loss of electrical energy by the example of rubbers of solution polymerization. The results of the study can be used in the development of software complexes for intelligent energy systems that allow tracking the dynamics of consumption and losses of energy resources. Using the results allows you to determine the optimal parameters of energy consumption and identify reserves for improving energy efﬁciency.


Introduction
For chemical industries, minimizing the negative impact on the environment caused by the use of fuel and electricity is a priority. Saving energy resources leads to a decrease in the load on the ecosystem and is one of the key objectives of the environmental policy. Energy Saving Programs are based on the principles of energy management, which are becoming increasingly popular in Russia and in other countries of the world. The use of the best available technologies and the achievements of energy management allows the enterprise management to build an efficient energy system in production with the lowest design and operating costs.
This research examines two interconnected topical areas-the energy component of the cost of production and the transformation of the structure of energy resources. The subject of the study is chemical industries. Chemical companies can make energy neural network to solve complex problems that do not have an explicit solution. In addition, the positive effects of using neural networks are: solving problems with unknown patterns between the input and output data; resistance to noise in the input data; adaptation to environmental changes, the ability to use neural networks in a non-stationary environment due to automatic retraining (adaptation); high performance due to parallel processing of information.
The high complexity and uncertainty of energy consumption processes make it preferable to use neural networks to study them [36][37][38][39][40]. In order to expand the mathematical tools for studying the energy efficiency of chemical industries, Jin et al. developed approaches to assessing the effectiveness of resource-saving projects based on a backpropagation neural network [41]; Xiang et al. proposed an algorithm for learning a multilayer neural network based on photonic spike timing dependent plasticity [42]; Menghi et al. systematized methods and tools for energy assessment [43]; Geng et al. assessed the energy saving of complex petrochemical industries based on complex clustering of the DEA affinity distribution [44]; Hamedi and Mokhtar applied multivariate linear regression and multilayer perceptronic artificial neural networks to build a baseline for energy consumption in a low-density polyethylene plant [45]; Shinkevich et al. showed the advantage of the network model of the value chain for analyzing the resource intensity of production [46]; Rajskaya et al. substantiated the expediency of a differentiated approach to assessing the level of innovative development of production resource-saving systems [47].
However, despite the presence of significant theoretical and methodological material and analytical data, there is a lack of research to solve the problems of efficient energy consumption by chemical industries. The lack of completeness of scientific and practical knowledge and solutions in the field under study does not allow us to objectively analyze the management processes of chemical and technological systems or to consider the specifics of production from the standpoint of optimal energy consumption.
One of the important issues that arise during energy analysis and monitoring of the implementation of energy-saving measures within the energy management system is the assessment of energy efficiency. This assessment is carried out on a few quantitative characteristics. At the same time, almost all Russian companies implementing energy efficiency programs aim to reduce the share of energy consumption in the cost of production. The urgency of this problem is high. The share of costs for electric and thermal energy in the cost of chemical products in Russia reaches an average of 13%. In comparison, in highenergy metallurgy, this figure is 9%, which confirms the relevance of energy conservation in the chemical industry. For many types of fuel and energy, the cost of resources in Russia is lower than, for example, in Western Europe. However, this does not always give Russian chemical companies an advantage due to the inefficient use of energy resources.

Basic Research Methods
As the main research methods, the article uses graphic and tabular tools for descriptive data analysis to study the dynamics of the structure of energy consumption and determine possible reserves for reducing its consumption. To determine the measure of the spread of the values of the specific weight of the resource relative to its average value, the variance of the values (D), which make up the data array on the structure of energy resources, was calculated using the formula: where x is the value of the specific weight of the energy resource, M is the average value of the specific weight of the energy resource in the energy portfolio, n is the number of periods.
The standard deviation of the population of data on the structure of energy resources (S), calculated as a percentage and more accessible for data interpretation, is defined as the positive square root of the variance of the population: To determine the possible reserve for reducing the share of energy in the consumption structure, the calculation of the interquartile distribution, namely the first lower or 0.25 quartile and the third upper or 0.75 quartile, is proposed. The reserve for reducing the specific weight of the resource (R) is defined as the difference between the actual value and the value of the lower-order quartile according to the formula: where F is the actual value of the specific weight of the resource for the period, Q 1 is the lower-order quartile, ∑f i is the sum of the frequencies of the series, n is the period, i is the type of energy resource.
To predict the maximum possible values of energy consumption, a statistical learning method is proposed-a method for training neural networks based on a regression predictive model. In contrast to traditional computing systems, optimization of parameters using neural networks allows solving problems with unknown patterns of development of the situation, an unstable relationship between input and output data. In the chemical industry, in relation to which the research was carried out, the probability of emergency situations, unstable operation of technological lines can determine the non-stationarity of energy consumption. This requires optimization systems to be adaptive for robust modeling. Neural networks have the ability to adapt to changes in external parameters and uncertainty of input data. Therefore, this method was chosen for modeling.
In this study, we used a type of multilayer neural network MLP, where one backpropagation algorithm trains all layers. Multilayer perceptrons MLP can be interpreted as an input-output model in which weights and constraints on value bias are free parameters of the model. Activation functions for hidden and output neurons can be both smoothed (logistic σ(t) = 1 1+e −t , hyperbolic th(t) = e t −e −t e t +t −t ) and the identity g(t) = t and exponent exp(t) = e t .
The iterative method of numerical optimization BFGS (Broyden-Fletcher-Goldfarb-Shanno algorithm) was used as a network learning algorithm, aimed in this case at finding the local minimum energy consumption at the upper limit of the range. The analysis methods were implemented using the Microsoft Office Excel software from the MS Office package and the STATISTICA software package.

Information Basis of the Research
The main informational study of the structure of fuel and energy resources in the chemical industry and forecasting the optimal consumption of energy resources is the dynamics of the energy component of the cost of production (Table 1) and the structure of fuel and energy resources in the production of chemicals from 2012 to 2019 (8 years) ( Table 2). The calculations used data from the Federal State Statistics Service (Russia), as well as data from Russian chemical enterprises. An important parameter of energy efficiency in industry is the energy component of the production cost or the share of fuel and energy resources in the production cost (Table 1). In 2019, according to the Federal State Statistics Service, in the production of chemicals and chemical products, the share of all fuel and energy resources was 13.25%. At the same time, the share of thermal energy in the cost of production was 6.13%, electric energy was 4.93% and water for technical needs was 0.18%. The share of fuel in the cost of production (gas, fuel oil, coal, gasoline, diesel fuel) amounted to 2.01% in total.
The averaged structure of fuel and energy resources in the production of chemicals and chemical products according to the Federal State Statistics Service included, in 2019, 37.2% of electric energy and 46.3% of thermal energy (Table 2). Traditionally, 1-2% falls on water. Among the types of fuel, combustible natural gas (12.7%) occupies a predominant position, and the other types account for no more than 2%. The movement of the structure of fuel and energy resources in dynamics for 2012-2019 for certain items is quite significant. This may be due to the change in the cost of energy resources, the transformation of the product range, the introduction of new technological equipment, as well as the effect of the implementation of corporate energy-saving programs.
Thus, to increase the competitiveness of Russian products, the priority direction of the energy conservation policy for the production of chemicals and chemical products should be to reduce the share of fuel and energy resources in the cost of production. Modernization of the energy management system at the enterprises of the industry will make it possible to implement a complex of interrelated elements in the form of organizational measures, technical means and software, and methodological support in the direction of energysaving and sustainable consumption of energy resources.

Descriptive Statistics of Changes in the Structure of Fuel and Energy Resources in the Production of Chemicals and Chemical Products
Analysis of the database on the energy component of the cost of chemical products showed the importance of studying the dynamics of the structure of fuel and energy resources used by enterprises in this industry for production. Complex chemical-technological systems have more complex reserves of energy conservation and have additional effects from the optimization of energy flows. At the enterprise level, the possibilities for modernization and development of energy systems are not limited to technical measures but are complemented by solutions and analytical methods for identifying reserves, as well as methods for improving management systems. Moreover, the introduction of technical tools is not exhaustive. The power supply of the chemical-technological system constantly requires adjustment due to the development of technical progress. The energy management system assumes the organization of system monitoring and analysis, where the energy component of the production cost is considered among the energy efficiency indicators. Monitoring and energy analysis are necessary in order to provide information support for the continuous process of improving the results of the use of energy resources. The database on the structure of energy resources used in the production of chemical products (Table 2) shows a significant variation in the dynamics of fuel and energy.
To simplify the interpretation of the results of graphical data analysis, Table 3 shows descriptive statistics of the structure of fuel and energy resources in the production of chemicals and chemical products in 2012-2019. Table 3. Descriptive statistics of changes in the structure of fuel and energy resources in the production of chemicals and chemical products in 2012-2019. A given space with random values of the consumption of fuel and energy resources X with distribution P X can be described by the distribution function:

Statistical Parameter
The distribution function is given by the formula: thus, the distribution function of the specific gravity of resource consumption X is the function F(x), the value of which at the point x is equal to the probability of the event {X ≤ x}, consisting only of the space of elementary events ω, expressed as X (ω) ≤ x. When evaluating the range of variables, the formulas for the lower boundary of the whisker X 1 and the upper boundary of the whisker X 2 are: where Q 1 is the first quartile (lower-order quartile), Q 3 is the third quartile, k is a coefficient of 1.5. Statistical parameters make it possible to generalize the primary results, present the distribution of their frequencies and identify the central trends in the change in the capacity of the energy portfolio relative to a resource. The reserve for reducing the specific weight of the resource is defined as the difference between the actual value of 2019 and the value of the lower order quartile Q 1 .
Analysis of the change in the share of electric energy in the energy portfolio of chemical industries produced some results. A positive trend is the greater number of cases with the share of electricity in the volume of resources of 31-33%, included in the quartile of the lower or first order (Q 1 = 32.13), which can be visually observed on the swing diagram. The variance of the analyzed resource relative to other resources has an average value D = 10.6. Accordingly, the standard deviation of the dynamic range of electrical energy S = 3.2%. The asymmetry coefficient A = 0.212 indicates a greater deviation of the values in the positive direction, i.e., the tendency for the predominance of values above the average value (M = 34.99%) remains. A negative kurtosis coefficient characterizes the absence of peaked data distribution. In general, the analysis shows a low spread in the share of electricity in the energy portfolio and a high frequency of minimum possible values, which is a positive point and shows the reserve for a decrease in the indicator at 5.07 percentage points.
An analysis of the change in the share of thermal energy is characterized by the largest magnitude of the swing (D = 27.71) and a deviation from the average value of 5.2%. At the same time, according to the coefficient of asymmetry (A = −0.108), high values of the specific weight of thermal energy, equal to 53-55%, prevail in the dynamic series. The distribution of the indicator values in dynamics is flat-topped and has no peak values (E = −2.418). The analysis of descriptive statistics shows the instability of the share of thermal energy in the total volume of energy resources and the predominance of high values of the indicator in the sample. In this case, the possible reserve for the decrease in the indicator will be less significant and amount to 0.95 percentage points.
An analysis of the change in the share of fuel (combustible gas, fuel oil, coal, gasoline, diesel fuel) in the energy portfolio of chemical industries shows the predominance in the dynamics of high values of the specific gravity of fuel (15-16%). The variance of the analyzed resource has a low value (D = 4.71). Accordingly, the standard deviation of the dynamic range of the fuel is insignificant (S = 2.2%). In general, the analysis shows a stable dynamic of the share of fuel resources with a bias of data towards higher values. With the indicator level in 2019 of 15.15%, the possible reserve for its reduction can be estimated at 3.4 percentage points.
The dynamic series of the share of water for technical needs in the energy profile of chemical industries shows a relatively even distribution of values with the established dynamics of the indicator of 1.7-1.4% in 2017-2019. Accordingly, there is a minimum variance (D = 0.18) and a series deviation (S = 0.42%). For this energy resource, the possible reserve for reducing the specific weight of water in the volume of resources is not more than 0.04 percentage points.
To determine the possible reserve for reducing the share of energy in the structure of resources, the interquartile distribution was calculated, namely the first lower or 0.25 quartile and the third upper 0.75 quartile. The reserve for reducing the specific weight of the resource (R) is defined as the difference between the actual value of 2019 and the value of the lower-order quartile by Formula 3.
Thus, the analysis showed the largest possible reserve for reducing the specific weight of electrical energy with existing production technologies. Given the high cost of this resource, the reduction in the cost of production due to energy conservation can be tangible. At the end of 2019, the share of electric energy in the cost of chemical production amounted to 4.93%. At the same time, the minimum value for the study period was noted at 4.5%. At the same time, thermal energy cannot be excluded from the list of priorities, the share of which in the cost of production in 2019 reached 6.13%.

Neural Network Models for Predicting Optimal Energy Consumption at the Upper Limit of the Range of Values
The method of statistical training of neural networks predicted the maximum possible values of energy consumption. As the initial data, the quarterly dynamics of the structure of energy resources consumption in the production of chemicals for 32 time periods was used (quarterly data for 2012-2019). The learning algorithm is based on a regression predictive model, where the objective function is: where C opt max is the optimal value of energy consumption at the upper limit of the range (maximum value), C 1−n min is the value of energy consumption at the lower limit of the range (minimum value) for the study period, C 1−n med is the average value of energy consumption for the study period, S n is the standard deviation of the value in the variational series, D n is the variance of the data series, n is the number of periods in the data series.
The categorical input variable is expressed in terms of C 1−n med , and the continuous input variable is expressed as C 1−n min . A linear regression relationship is determined between the continuous target variable and the input variables. Accordingly, the continuous target variable is the most optimal energy consumption value at the upper limit of the C opt max range. Since the structure of the energy portfolio is limited by a certain set of energy resources with limiting values of consumption according to chemical technology and the level of reliability (95.0%) of the parameters of the time series, the following conditions must be satisfied: A block diagram of finding the optimal energy consumption by the method of training neural networks is shown in Figure 1. The procedure includes five main stages:

•
Step 1. Primary statistical processing of an array of data on the structure of energy resources consumption with obtaining the necessary parameters for training neural networks.

•
Step 2. Statement of the optimization problem-determination of the minimum possible value of the objective function-the parameter of energy resource consumption at the upper limit of the range, considering the restrictions.

•
Step 3. Selection of neural network training parameters: BFGS iterative numerical optimization method as a network training algorithm and MLP type of neural network.

•
Step 4. The procedure for optimizing the target parameter under the condition of maintaining the integrity of the energy portfolio structure ∑ C opt max E 1−9 = 100. The algorithm is executed if the inequalities under the constraints lim min ≤ ∑ C opt max E n ≤ lim max are true.

•
Step 5. Determination of the target optimal parameter by choosing the smallest value ∑ C opt max E n from the variants of training cycles-iterations N with the true inequality lim min ≤ ∑ C opt max E n ≤ lim max . Input parameters were set at 20 networks for training and 5 networks for saving. The results of the neural network models are presented in Table 4.  Input parameters were set at 20 networks for training and 5 networks for saving. The results of the neural network models are presented in Table 4. According to the resulting architecture, each of the saved MLP networks has five input variables and one variable at the output C opt max according to the objective function, and the number of hidden neurons varies from 3 to 8. The highest performance or specific weight of correct classification on the training and control samples is possessed by 3-5 neural networks (96-99%). The iterative method of numerical optimization BFGS (Broyden-Fletcher-Goldfarb-Shanno algorithm) was used as a network learning algorithm, aimed in this case at finding the local maximum consumption of energy resources. The optimal value at the highest level of control performance (network No. 4) was found in 24 training cycles.
In the study, the object of statistical training is the structure of energy resources of chemical industries, which includes nine components: thermal energy (E 1 ), electrical energy (E 2 ), water (E 3 ), natural combustible gas (E 4 ), diesel fuel (E 5 ), fuel oil (E 6 ), coal (E 7 ), gasoline (E 8 ) and other fuels (E 9 ). The specific weight of each resource E i in the general structure of energy resources is as follows: Thus, as a result of training the neural network (network No. 4 MLP 5-3-1), the optimal values of energy consumption were obtained at the upper limit of the range of values (maximum values), which, as can be seen, in most cases do not have a significant deviation from the objective function (Table 5). During the statistical processing of the data, the resources of natural combustible gas (E 4 ) and gasoline (E 8 ) were removed from the sample due to the presence of values with a high level of variation. Table 5. Neural network training results: optimal consumption in the structure of energy resources at the upper limit of the range of values.

Energy Types
Objective Function: The optimal Value of C The optimal structure of the energy resources of chemical industries contains E 1 = 55.48%, E 2 = 39.59%, E 3 = 2.35%, E 5 = 1.71%, E 6 = 0.21%, E 7 = 0.19%, E 9 = 0.48%. Undoubtedly, the absence of the primary specified components E 4 and E 8 in the structure of energy resources requires additional analysis with the participation of engineers in the field of chemical technology.
The use of the neural network training method made it possible to simulate the optimal structure of energy resources of chemical production, which is an important component of the strategy of resource conservation and increases the competitiveness of the Russian industry. In addition, a significant parameter of resource efficiency of production is the loss of fuel and energy resources. Losses of energy resources in distribution networks are inevitable; therefore, it is important that they do not exceed a reasonable level. An increase in the rates of technological consumption of fuel and energy resources indicates the problems that have arisen, and to eliminate them it is necessary to identify the causes of inappropriate costs and determine the directions of their reduction.  (Figure 2). However, the amount of electricity losses of 15.2% of the total volume of transmitted energy remains relatively high in Russia. In comparison, in industrially developed countries it is as follows: Germany-7.6%, Finland-6.7%, USA-6.5%, Canada-6.3%, Japan-5%. component of the strategy of resource conservation and increases the competitiveness of the Russian industry. In addition, a significant parameter of resource efficiency of production is the loss of fuel and energy resources. Losses of energy resources in distribution networks are inevitable; therefore, it is important that they do not exceed a reasonable level. An increase in the rates of technological consumption of fuel and energy resources indicates the problems that have arisen, and to eliminate them it is necessary to identify the causes of inappropriate costs and determine the directions of their reduction.

Graphical Statistics of the Dynamics and Interdependence of Losses of Electrical and Thermal Energy in the Distribution Networks of Industrial Production
For the studied years 2012-2019, there is a positive trend in reducing the amount of energy losses, both in electrical (reduction index in 2012-2019-93.8%) and thermal energy (reduction index in 2012-2019-93.3%) (Figure 2). However, the amount of electricity losses of 15.2% of the total volume of transmitted energy remains relatively high in Russia. In comparison, in industrially developed countries it is as follows: Germany-7.6%, Finland-6.7%, USA-6.5%, Canada-6.3%, Japan-5%.  Correlation coefficient of heat and electric energy losses Kor = 0.98. This indicates a close relationship between the parameters. A stable reduction in the amount of energy losses in distribution networks is probably not a random phenomenon but is due to purposeful management actions of executive authorities and enterprise management [48].  Correlation coefficient of heat and electric energy losses Kor = 0.98. This indicates a close relationship between the parameters. A stable reduction in the amount of energy losses in distribution networks is probably not a random phenomenon but is due to purposeful management actions of executive authorities and enterprise management [48]. Russian chemical producers have indicators of the consumption of fuel and energy resources. These indicators are published in the industry's best available technology reference books. These guides give the minimum and maximum possible consumption of energy resources for specific types of products. The range between the minimum and maximum values of resource consumption is justified by existing traditional and innovative production technologies, the degree of deterioration of energy networks and the specificity of the chemical technology of the product.
For example, the Russian polymer manufacturer PJSC "Nizhnekamskneftekhim" produces five grades of solution polymerized rubbers with a total volume of about 1.9 million tons per year. Each type of rubber has a standardized minimum and maximum value of the specific consumption of energy resources. The amount of excess of the max value over the min value of the resource consumption (r) is calculated by the formula: Pr = Cmax(r) / Cmin(r) × 100. (11) The calculation results are shown in Table 6. The value of the excess of the max flow rate over the min value differs significantly. Synthetic butadiene-styrene rubber DSSK-2565F and synthetic cis-butadiene lithium rubber SKD-L have the highest specific energy consumption, for which the production of one ton, respectively, requires 667-693 kWh and 750-790 kWh of electric energy and 6.1-8.3 Gcal and 6.5-6.9 Gcal of thermal energy. As can be seen, for these types of rubber, there is a smaller corridor between the limiting values of resource consumption: the excess of min over max Pr = 103.9% and Pr = 105.3% for electrical energy. For thermal energy, the range of limit values is somewhat wider than for electrical energy. The reason for this phenomenon is the technological features of heat exchange processes and equipment.  Russian chemical producers have indicators of the consumption of fuel and energy resources. These indicators are published in the industry's best available technology reference books. These guides give the minimum and maximum possible consumption of energy resources for specific types of products. The range between the minimum and maximum values of resource consumption is justified by existing traditional and innovative production technologies, the degree of deterioration of energy networks and the specificity of the chemical technology of the product.
For example, the Russian polymer manufacturer PJSC "Nizhnekamskneftekhim" produces five grades of solution polymerized rubbers with a total volume of about 1.9 million tons per year. Each type of rubber has a standardized minimum and maximum value of the specific consumption of energy resources. The amount of excess of the max value over the min value of the resource consumption (r) is calculated by the formula: P r = C max(r) / C min(r) × 100.
The calculation results are shown in Table 6. The value of the excess of the max flow rate over the min value differs significantly. Synthetic butadiene-styrene rubber DSSK-2565F and synthetic cis-butadiene lithium rubber SKD-L have the highest specific energy consumption, for which the production of one ton, respectively, requires 667-693 kWh and 750-790 kWh of electric energy and 6.1-8.3 Gcal and 6.5-6.9 Gcal of thermal energy. As can be seen, for these types of rubber, there is a smaller corridor between the limiting values of resource consumption: the excess of min over max Pr = 103.9% and Pr = 105.3% for electrical energy. For thermal energy, the range of limit values is somewhat wider than for electrical energy. The reason for this phenomenon is the technological features of heat exchange processes and equipment. Table 6. Limit values of the specific consumption of energy resources to produce rubbers for solution polymerization. The range of limit values of specific energy consumption as a percentage of Dcr is expressed mathematically through the excess of the max value over the min value: Dcr = Pr-100. The dependence of the range of limit values of specific energy consumption on the average consumption rate is confirmed in the scatter diagrams. Figure 4 shows an absolute inverse relationship between the values of the electrical capacity of rubber and the range of the min-max values corridor (correlation coefficient Kor = −0.77). Therefore, it can be assumed that in the production of high-power products, the rate of loss of electrical energy should be significantly lower. Attention should be paid to technical losses of electricity due to physical processes in electrical equipment, avoidance of idle time and downtime of high-power installations. The range of limit values of specific energy consumption as a percentage of Dcr is expressed mathematically through the excess of the max value over the min value: Dcr = Pr-100. The dependence of the range of limit values of specific energy consumption on the average consumption rate is confirmed in the scatter diagrams. Figure 4 shows an absolute inverse relationship between the values of the electrical capacity of rubber and the range of the min-max values corridor (correlation coefficient Kor = −0.77). Therefore, it can be assumed that in the production of high-power products, the rate of loss of electrical energy should be significantly lower. Attention should be paid to technical losses of electricity due to physical processes in electrical equipment, avoidance of idle time and downtime of high-power installations. The dependence of the range of limiting values of the specific consumption of heat energy on the average consumption rate is weakly negative, with small tightness of connection (Kor = −0.35) ( Figure 5). The lower correlation of the parameters is due to the use of secondary heat resources in sufficiently large volumes due to the integration of cold and hot flows in chemical processes. It should also be noted that in addition to purchasing heat energy on the market, enterprises can produce it independently, spending electric energy on it. Nevertheless, we consider it expedient to use the criterion of the average specific resource consumption to produce products when calculating the rates of heat loss. The dependence of the range of limiting values of the specific consumption of heat energy on the average consumption rate is weakly negative, with small tightness of connection (Kor = −0.35) ( Figure 5). The lower correlation of the parameters is due to the use of secondary heat resources in sufficiently large volumes due to the integration of cold and hot flows in chemical processes. It should also be noted that in addition to purchasing heat energy on the market, enterprises can produce it independently, spending electric energy on it. Nevertheless, we consider it expedient to use the criterion of the average specific resource consumption to produce products when calculating the rates of heat loss.

Development of Optimal Values of the Rate of Losses of Electrical Energy for Rubbers of Solution Polymerization Based on Training a Neural Network
To determine the rate of losses of electrical energy Ni in the production of rubbers of solution polymerization, depending on the range of limit values of energy consumption proposed in the handbook on the best available technologies, a regression method was used for training neural networks. The objective function is: ( ) p r is the maximum value of electricity consumption for polymer, p is polymer (solution polymerization rubber).
The optimization problem is to find the minimum possible value of the rate of electrical energy losses for the polymer (solution polymerization rubber), considering the established restrictions. The limits, according to the proposed approach, are determined based on the minimum and maximum values of electricity consumption for polymer in accordance with the Standard for Best Available Techniques [49]: The block diagram of finding the optimal energy consumption by the method of training neural networks is shown in Figure 6. The sequence of operations is according to the block diagram for determining the optimal energy consumption (see Figure 1). As a network learning algorithm, the BFGS numerical optimization method was also used to calculate the inverse Hessian based on changes in gradient values and changes in weights. The gradient vector of the error function is computed using the usual backpropagation procedure. The results of the neural network models are presented in Table 7. All five neural networks have a high proportion of correct classification in the training and control samples. The highest training performance is observed in the second neural network MLP 5-9-1, representing five input variables, one variable at the output Ni p opt according to the objective function and nine hidden neurons. At the same time, the optimal value at the highest level of control performance was found in 32 training steps.
As a result of network training, synthetic cis-isoprene rubber SKI-3 was removed from the sample. The optimal values of the rate of electrical energy losses for the four types of solution polymerization rubbers in percentage are presented in Table 8. Table 8. Results of neural network training: the optimal value of the rate of electrical energy losses for solution polymerization rubber as a percentage.

Rubber
Objective Function: Thus, as a result of training the neural network, the optimal values of the rate of electrical energy losses were obtained. Note that the background information for statistical training is based on the best available technology reference material. When designing or modernizing the production of chemicals and chemical products, it is necessary to revise the norms for the loss of energy resources, considering, along with other parameters, the absolute value of the level of energy intensity of products. Using the example of solution polymerization rubbers using the neural network training method, the boundary values of the rate of electrical energy losses have been developed (Table 9): • At a level of electricity consumption of more than 500 kWh per ton of product, the rate of electricity losses should not exceed 7% of the total consumed resource on average; • At a level of electricity consumption of 300-500 kWh per ton of product, possible resource losses should not be more than 10%; • At a power consumption level of up to 300 kWh per ton of product, the rate of resource loss is no more than 13%.

Min Electricity Consumption Per Ton of Product Electricity Loss Rate Rubber
More than 500 kWh No more than 7% Synthetic cis-butadiene lithium rubber SKD-L Synthetic butadiene-styrene rubber DSSK-2565F 300-500 kWh No more than 10% Synthetic cis-isoprene rubber SKI-3 Synthetic cis-butadiene neodymium rubber SKD-N Up to 300 kWh No more than 13% Halobutyl rubber GBK These boundary values were proposed by the authors on the basis of the above results of descriptive analytics and training of neural networks with respect to the time series of energy resources and the dependence of the range of limiting values of energy consumption on the average rate established in reference books on the best available technologies.

Discussion
The study of the dynamics of the energy component of the production cost, structure and maximum consumption of energy resources by chemical enterprises allows us to draw the following conclusions:

1.
Various methods can be used to solve the research problems posed by the authors. Comparative analysis showed that it is advisable to use Artificial Neural Networks for modeling and optimization. This tool makes it possible to solve problems with unknown patterns of development of the situation and an unstable relationship between input and output data. In chemical production, the likelihood of emergency situations and the instability of the operation of technological lines determines the non-stationary consumption of energy resources. This requires optimization systems to be adaptive for robust modeling. Neural networks have the ability to adapt to changes in external parameters and the uncertainty of the input data. Therefore, Artificial Neural Networks were chosen as the main modeling and optimization tool.

2.
In Russian chemical industries, at least 13% of the cost of production is made up of fuel and energy resources, which includes 37.2% of electrical energy and 46.3% of thermal energy in their structure. A retrospective analysis of the structure of the energy portfolio of the industry shows a visible transformation of the proportion of resources due to a change in their cost, the introduction of new technological lines, as well as the effect of the implementation of corporate energy conservation programs. The descriptive method for analyzing data on the structure of energy resources of chemical industries made it possible to estimate the potential reserves for reducing the specific weight of electric energy with existing production technologies at 5.07 percentage points, heat energy at 0.95 percentage points, fuel at 3.4 percentage points and water for technical needs at 0.04 percentage points.

3.
A positive trend has been established to reduce the amount of energy losses at Russian industrial enterprises in 2012-2019 by 6.2% for electrical energy and 6.7% for thermal energy. On average, the value of electricity losses in 15.2% of the total volume of transmitted energy remains high in comparison with industrialized countries, where the maximum values of losses are in the range of 7-8%. Based on the revealed correlation dependences of the range of limiting values of electricity consumption from the average rate of energy consumption, the necessity of taking into account the absolute value of the energy consumption of products when determining the rates of losses of energy resources at the design stage or modernization of chemical production is justified.

4.
By the method of training neural networks using a regression predictive model, the optimal values of energy consumption at the upper limit of the range of values are determined, which in most cases do not have a significant deviation from the objective function. Using the BFGS numerical optimization method, the optimal structure of the energy portfolio of chemical industries is obtained, containing 55.4% of thermal energy, 39.6% of electric energy, 2.3% of water for technical needs and 2.7% of fuel resources.

5.
Approbation of the proposed approach based on statistical training to determine the norms of losses of electrical energy on the example of rubbers of solution polymerization showed that at a level of resource consumption of more than 500 kWh per ton of product, the rate of losses of electrical energy should not exceed 7%; at a consumption level of 300-500 kWh, no more than 10%; at a flow rate of up to 300 kWh, no more than 13%. The research materials can be used in the development of software complexes for intelligent energy systems that allow tracking the dynamics of consumption and losses of energy resources, determine the cause-and-effect relationships of deviations in resource consumption from a given trajectory, and establish the optimal vector of sustainable energy consumption and parameters of reserves for increasing energy efficiency. Funding: This research received no external funding.