Disclosure of Invention
Aiming at the problems, the invention provides a circulating cooling water minimum pressure difference energy-saving control system and method based on a process medium multi-temperature target. The cooling supply quantity of the circulating cooling water is automatically changed according to the cooling requirement of the industrial field process medium, and the energy-saving control is realized. In order to achieve the purpose, the invention adopts the following specific technical scheme:
a circulating cooling water minimum pressure difference energy-saving control system based on a process medium multi-temperature target comprises a circulating cooling water system and a control energy-saving system; the circulating cooling water system comprises N cooling towers, a water suction pool, a water supply pump unit, a water outlet pipe group, M 'heat exchangers and a water return pipe group, wherein the N cooling towers, the water suction pool, the water supply pump unit, the water outlet pipe group, the M' heat exchangers and the water return pipe group are arranged along a circulating water path; a cooling pool is arranged in the cooling tower; the water outlet pipe group comprises N water inlet pipes of the water suction tanks, L water outlet pipes of the water suction tanks, a water supply main pipe and M' water supply branch pipes; the water return pipe group comprises M' water return branch pipes, a water return main pipe and N upper tower water return pipes; any cooling tower is correspondingly connected with the water suction pool through one water suction pool inlet pipe, the water suction pools are connected with the water supply main pipe through L water suction pool outlet pipes in parallel, and the water supply main pipe supplies water to M 'heat exchangers in a one-to-one correspondence manner through M' water supply branch pipes; a water feeding pump is arranged on the water outlet pipe of the water suction pool in parallel; any heat exchanger is correspondingly connected with the water return main pipe through a water return branch pipe, and the water return main pipe is correspondingly connected with N cooling towers one by one through N upper tower water return pipes which are arranged in parallel; process medium temperature sensors are arranged in the M' heat exchangers and used for acquiring real-time temperature detection values of various process media;
a cooling water supply temperature sensor and a cooling water supply pressure sensor are arranged on the water supply header pipe, the cooling water supply temperature sensor is used for acquiring a cooling water supply temperature detection value, and the cooling water supply pressure sensor is used for acquiring a cooling water supply pressure detection value; a cooling return water temperature sensor and a cooling return water pressure sensor are arranged on the return water main pipe, the cooling return water temperature sensor is used for acquiring a cooling return water temperature detection value, and the cooling return water pressure sensor is used for acquiring a cooling return water pressure detection value; an upper tower valve is respectively arranged on the N upper tower return water pipes;
the key technology is as follows: the control energy-saving system comprises a process medium temperature least-unfavorable point selector, a process medium temperature controller and a water supply and return pressure difference inner ring PID controller; the process medium temperature minimum point selector obtains a process medium temperature minimum point according to a process medium temperature deviation value sequence collected by the circulating cooling water system; the process medium temperature controller is used for obtaining a cooling water supply and return water pressure difference set value by combining field operation data according to the least unfavorable point of the process medium temperature;
obtaining a cooling water supply pressure detection value and a cooling water return pressure detection value which are obtained by the detection of the circulating cooling water system, and then obtaining a cooling water supply water pressure difference detection value; the difference between the cooling water supply and return pressure difference detection value and the cooling water supply and return pressure difference set value is obtained to obtain a cooling water supply and return pressure difference deviation value; and the water supply and return pressure difference inner ring PID controller adjusts the opening of all the upper tower valves according to the cooling water supply and return pressure difference deviation value, so that the outlet water flow of the water supply pump unit is changed.
Through the design, on the basis of a circulating cooling water system, a circulating cooling water minimum differential pressure energy-saving control system based on a process medium multi-temperature target is established by combining water pump variable-frequency variable-pressure water supply control. The difference value of the process medium online temperature detection value and the temperature set value is called process medium temperature deviation, and the minimum point of the temperature deviation of each process medium in the whole heat exchanger group is called the process medium temperature least-beneficial point. The circulating cooling water system controls the temperature of all process media in an industrial field to be maintained in a normal, safe and energy-saving production temperature range according to the temperature deviation and deviation change rate of the most unfavorable point of the process media temperature and the real-time detection value of the cooling water supply temperature, determines the required amount of cooling water flow in real time, adjusts an upper tower valve and changes the outlet amount of cooling water of a water pump. Finally, on the basis of meeting the core requirement on the production process medium temperature control, the energy conservation of a water pump motor in the circulating cooling water system is further realized fundamentally.
The further technical scheme is as follows: the process medium temperature deviation value sequence comprises M' process medium temperature deviation values, and any process medium temperature deviation value is equal to the difference value between the corresponding process medium real-time temperature detection value and the corresponding process medium temperature set value; the process medium temperature deviation change rate sequence comprises M' process medium temperature deviation change rates, wherein any one of the process medium temperature deviation change rates is the ratio of the temperature change value of two adjacent detection time periods corresponding to the process medium to the last detection time period.
By adopting the scheme, the temperature acquisition of the process media is realized by combining with a circulating cooling water system, the difference values in all the process media are obtained by combining with the detection value and the set value, and the difference values are used for searching the heat exchanger and the process media which are closest to the temperature set threshold value to control the upper tower valve of the whole system, so that each process medium is ensured to be kept in a safe and energy-saving temperature range.
A control method of a circulating cooling water minimum differential pressure energy-saving control system based on multiple temperature targets of process media is characterized by comprising the following specific steps:
s1: setting a sampling period, and carrying out field operation data acquisition on the circulating cooling water system;
the circulating cooling water system collects real-time temperature detection values of M 'process media in M' heat exchangers on site through the process media temperature sensor; the circulating cooling water system collects a cooling water supply temperature detection value on site through the cooling water supply temperature sensor; the circulating cooling water system collects a cooling water supply pressure detection value on site through the cooling water supply pressure sensor; the circulating cooling water system collects a cooling return water temperature detection value on site through the cooling return water temperature sensor; the circulating cooling water system collects a cooling return water pressure detection value on site through the cooling return water pressure sensor;
s2: obtaining M 'process medium temperature deviation values and M' corresponding process medium temperature deviation change rates according to M 'process medium real-time temperature detection values acquired on site and M' corresponding process medium temperature set values; m' process medium temperature deviation values form the process medium temperature deviation value sequence; m' process medium temperature deviation change rate values form the process medium temperature deviation change rate sequence;
s3: the process medium temperature least-benefit point selector selects the temperature deviation minimum value in the process medium temperature deviation value sequence as the process medium temperature least-benefit point, and obtains the process medium temperature deviation change rate corresponding to the process medium temperature least-benefit point and the corresponding heat exchanger;
s4: the process medium temperature controller acquires a process medium temperature deviation value and a process medium temperature deviation change rate of a corresponding heat exchanger according to the least favorable point of the process medium temperature, and inputs field operation data of the corresponding heat exchanger into a corresponding process medium temperature prediction model to obtain a cooling water supply and return water differential pressure set value and a process medium temperature predicted value;
s5: the cooling water supply pressure detection value and the cooling water return pressure detection value are subjected to difference to obtain a cooling water supply and return pressure difference detection value; and the difference between the set value of the cooling water supply and return pressure difference and the detected value of the cooling water supply and return pressure difference is used for obtaining a cooling water supply and return pressure difference deviation value, and the cooling water supply and return pressure difference deviation value is sent to the water supply and return pressure difference inner ring PID controller to adjust the opening of the upper tower valve, so that the outlet water flow of the water supply and supply pump unit is changed.
According to the method, the closed loop for controlling the water temperature of the heat exchanger in the production area monitors the cooling condition of each production area to be cooled in real time, and the opening degree of the upper tower valve is adjusted to accurately control the required cooling water amount in the heat exchanger of each production area to be cooled, so that the waste of cooling resources is reduced while the temperature of the process medium is controlled within a safe production range, energy is saved, and the cooling effect is better. When the opening of the upper tower valve reaches the maximum and the temperature of the process medium cannot be effectively controlled, a measure for increasing the refrigeration effect of the cooling tower can be taken.
Still further, the specific step of finding the minimum value in the sequence of the process medium temperature deviation values as the process medium temperature least favorable point by using the process medium temperature least favorable point selector in the step S3 includes:
s31: initializing, setting M 'technological medium temperature deviation values to form a difference group, totaling M' technological medium temperature deviation values, and enabling W to bek=M’;k=1
S32: let Wk+1W is equal to Wk + Xk+1Can be divided by M' evenly, X is a gap with large filling difference; and X is equal to 0-M' -1;
s33: calculating Wk+2=Wk+1/M’
S34: from Wk+2In the group, the minimum value is found out from the M' process medium temperature deviation values of each group by adopting a cross comparison method to obtain Wk+2A process medium temperature deviation value;
s35: judgment of Wk+2Whether it is equal to 1; if so, taking the process medium temperature deviation value as the least unfavorable point of the process medium temperature; otherwise, let k be k + 2; return is made to step S32.
And a most unfavorable point selector is additionally arranged, and a large number of minimum values of temperature deviations of the process medium are searched to determine the most unfavorable point of the process medium temperature and pay attention to the real-time performance and accuracy of the selection of the most unfavorable point of the process medium.
Still further, the building steps of the M' process medium temperature prediction model deep learning neural networks in the step S4 are as follows:
s411: numbering M' heat exchangers, acquiring historical data generated by running of a circulating cooling water system in X sampling periods, and screening the acquired historical data according to screening conditions to obtain training data of a process medium temperature prediction model in the heat exchangers;
s412: determining characteristic variables from historical data, taking process medium temperature deviation, process medium temperature deviation change rate, cooling water supply temperature detection value and cooling water supply and return pressure difference detection value in the historical data of the heat exchanger I as input data, performing data normalization processing to obtain a normalized data set, and dividing the normalized data set into a training sample set and a test sample set;
s413: performing greedy unsupervised pre-training on the training sample set layer by layer based on a stacked automatic encoder to obtain a weight matrix W of an input layer and a hidden layer initialized by a deep learning neural network and a threshold matrix B of the input layer and the hidden layer;
s414: and (3) fine adjustment of parameters: finely adjusting a weight matrix W of an input layer and a hidden layer initialized by a deep learning neural network and a threshold matrix B of the input layer and the hidden layer until the iteration times reach the maximum value of the iteration times to obtain an initial process medium temperature prediction model based on a stacking automatic encoder;
s415: and (5) evaluating the initial process medium temperature prediction model obtained in the step (S414) by using the test sample data set to obtain a process medium temperature prediction model based on the stacking automatic encoder.
Compared with the traditional shallow neural network, the process medium temperature prediction model based on the stacked automatic encoder effectively solves a series of problems caused by the random initialization of the parameters of the traditional neural network, can effectively mine the implicit relation of each data, and greatly improves the accuracy of process medium temperature prediction control in the industrial production field; compared with the traditional stack automatic encoder method, the method can improve the prediction control precision of the temperature of the process medium in the industrial production field. The method is applied to the energy-saving control of the circulating cooling water system in the industrial production field, can quickly and accurately give the set value of the cooling water supply and return water pressure difference, is beneficial to industrial production field management personnel to master the temperature change trend of the process medium in each heat exchanger in real time, and most importantly, reduces the energy consumption of the circulating cooling water system.
Still further, the characteristic variables comprise a process medium temperature deviation value, a process medium temperature deviation change rate, a process medium real-time temperature detection value, a cooling feed water temperature detection value and a cooling feed water return pressure difference detection value in the circulating cooling water system in any heat exchanger.
Still further, the screening condition is historical data that the process medium temperature is within a safe and energy-saving temperature value interval, and the safe and energy-saving temperature value interval is within a process medium temperature threshold value interval.
The difference between the safe and energy-saving temperature value interval and the process medium temperature threshold value interval is set after induction according to the historical experience technology of technicians.
Further, in step S413, when greedy unsupervised pre-training is performed on the training sample set layer by layer, the training sample set is divided into P groups of small batch training samples, training is performed in sequence, a Dropout technique is adopted, a part of neurons are randomly selected to pause working, iteration is performed in sequence, training is performed layer by layer, and weight matrix W, input layer threshold matrix B and hidden layer threshold matrix B of the input layer and hidden layer initialized by the neural network are deeply learned.
In step S413, the stacked automatic encoder sae (stacked automatic encoder) is a typical deep learning neural network, and its basic constituent unit is an Automatic Encoder (AE), and its network structure is composed of an encoder and a decoder: the input vector is mapped to a feature vector in the hidden layer by the encoder, and then the feature phasor is reconstructed to the original input vector by the decoder.
When given an input sample set X ═ XiI is more than or equal to 1 and less than or equal to N, wherein N is the total number of samples, and x isiThe ith training sample in the sample set has dimension n. Let H ═ HiI is more than or equal to 1 and less than or equal to N is a hidden layer characteristic vector set, hiThe feature vector corresponding to the ith sample has a dimension M', and the coding relationship between X and H is:
H=sf(WX+B)
in the formula: w is a weight matrix of the input layer and the hidden layer; b is an input layer and hidden layer threshold matrix; sfFor the neuron activation function of the encoder, a Sigmoid function is usually used, which has good feature recognition:
sf(z)=1/(1+exp(-z))
in the formula: z is the input vector.
The decoder is the inverse operation of the encoder, and takes the characteristic vector of the hidden layer as the input vector
In order to output the set of vectors,
and an output vector corresponding to the ith sample is obtained, and the dimension is n, so that the expression of the decoder is as follows:
in the formula: w' is a weight matrix of the hidden layer and the output layer; b' is a hidden layer and a transport layerA threshold matrix of the layer; sg is the neuron activation function of the decoder. The automatic encoder achieves the purpose of feature extraction by minimizing the reconstruction error between the output vector and the input vector, and the formula of the reconstruction error is as follows:
and continuously adjusting the network weight and the threshold by using a gradient descent algorithm to reduce the reconstruction error, wherein the formula is as follows:
in the formula: l is the learning rate;
to represent
Calculating the deviation of the weight W;
to represent
The bias is calculated for the threshold B.
The SAE is a deep learning neural network model composed of stacked AE, and the output of the lower AE is used as the input of the upper AE. A gradual abstraction of features is achieved by stacking of AEs, ultimately resulting in more compact, useful features. WnIs a weight matrix of the n-1 st hidden layer and the n-th hidden layer, BnThe threshold matrix is the n-1 st hidden layer and the nth hidden layer. The training process comprises two steps of greedy layer-by-layer unsupervised pre-training and supervised fine tuning. Greedy non-supervision pre-training layer by layer obtains the initial weight and the threshold value of the network through layer by layer training, the input of the bottom layer AE is original data, and the output data of the hidden layer is used as the input data of the upper layer AE.
After the layered pre-training is completed, the hidden layers are stacked, and the relationship between the input data and the output data is expressed as follows:
in the formula: f is an activation function, x
iW, B are respectively the network initialization weight and threshold value obtained by pre-training layer by layer,
is the predicted value of the ith sample. Constructing an error loss function of the actual value and the predicted value, wherein the formula is as follows:
in the formula: n is the total number of samples, yiIs the actual value of the ith sample. The whole network weight and the threshold value are finely adjusted through the back propagation from top to bottom, and the error between the predicted value and the actual value is reduced.
The SAE can effectively extract high-order features of data through greedy unsupervised pre-training layer by layer, better approaches complex functions, reduces parameter optimization space, can quickly obtain network parameters, and improves deep feature learning capability of a neural network.
Compared with a shallow machine learning algorithm, when the deep network processes high-dimensional data, the overfitting problem is more easily generated due to the complex network structure of the deep network, so that the generalization capability of the model is limited.
Dropout is a mainstream anti-overfitting technique, and its basic idea is: during model training, a part of nodes are randomly selected to be out of operation, the nodes store the weight of the last iteration, and the output is set to be 0. And the selected nodes restore the weight values reserved before in the process of the next iteration, and part of the nodes are randomly selected again to repeat the process. The network structure is changed to a certain extent in each iteration process, and by adopting a Dropout technology, part of neurons selected randomly do not work temporarily, so that the combined action among specific nodes is reduced, the dependence of network output on the states of the specific nodes is reduced, and overfitting is prevented.
Further, when the initial process medium temperature prediction model is evaluated, the improved stack automatic encoder after test training is tested by using a test sample data set, and the average percent error (MAPE) is used as a standard for measuring the evaluation performance of the improved stack automatic encoder, wherein the expression is as follows:
in the formula: y is
i、
And respectively obtaining an actual value of the pressure difference of the circulating cooling water supply and return water of the ith sample and a predicted value obtained through SAE evaluation. In step S414, when the parameters are fine-tuned, a top-down small-batch RM' SProp optimization method is used to fine-tune the weight matrix W threshold matrix B initialized by the deep learning neural network.
The method comprises the following specific steps: the global learning rate l and the decay rate ρ are set, and the initial cumulative variable r1 is 0 and r2 is 0.
Selecting a small batch of data sets containing M' samples from the training set, and calculating the gradient according to an error loss function:
calculating the cumulative squared gradient, as shown in equation (9):
in the formula: as an element-by-element product symbol.
Updating the weight and threshold parameters, respectively:
and when the iteration times reach the requirements, stopping operation, otherwise, returning to the step 2 to continue to execute the calculation.
Further describing, in step S5, the specific control content of the supply and return water pressure difference inner ring PID controller on the opening of the upper tower valve and the outlet water flow of the water supply pump unit is as follows:
if the deviation value P of the cooling supply and return water differential pressure is greater than 0, the supply and return water differential pressure inner ring PID controller sends an upper tower valve amplification control signal for controlling the upper tower valve to increase the valve opening; if the cooling water supply and return pressure difference deviation value P is equal to 0, the water supply and return pressure difference inner ring PID controller keeps the original control state;
if the deviation value P of the cooling supply and return water pressure difference is less than 0, the supply and return water pressure difference inner ring PID controller sends a control signal for adjusting the upper tower valve to be small, and the control signal is used for controlling the upper tower valve to reduce the opening of the valve; through the design, the opening degree of the upper tower valve is adjusted by combining the differential pressure deviation value of cooling water supply and return water, so that the temperature of all process media is ensured to be within the threshold value. And ensures the environmental protection performance.
The invention has the beneficial effects that: the cooling water supply quantity of the circulating cooling water can be automatically changed according to the cooling requirement of the industrial field process medium without manual adjustment, the energy utilization rate of the circulating cooling water system is improved, and energy-saving control is realized.
Detailed Description
The following provides a more detailed description of the embodiments and the operation of the present invention with reference to the accompanying drawings.
A minimum differential pressure energy-saving control system of circulating cooling water based on multiple temperature targets of process media comprises a circulating cooling water system and an energy-saving control system.
In this embodiment, the chemical plant a is taken as an example, and the circulating cooling water energy saving control is performed. As can be seen from fig. 1, a process diagram of a circulating cooling water system of a chemical plant a, which can be seen in combination with fig. 1, includes N cooling towers, a water suction pool, a water supply pump unit, a water outlet pipe group, M' heat exchangers and a water return pipe group, which are arranged along a circulating water path; a cooling pool is arranged in the cooling tower; in this embodiment, M' is a positive integer, there are three heat exchanger groups: a synthesis zone heat exchanger group, a urea zone heat exchanger group and a power zone heat exchanger group.
The water outlet pipe group comprises N water inlet pipes of the water suction tanks, L water outlet pipes of the water suction tanks, a water supply main pipe and M' water supply branch pipes; in the present embodiment, the number of heat exchangers M' is 50. Wherein, 50 heat exchangers are divided into three heat exchanger groups, namely a synthesis area heat exchanger group, a urea heat exchanger group and a power heat exchanger group. The three heat exchanger groups have 50 heat exchangers in total. In this embodiment, the water return pipe group includes M' water return branch pipes, a water return main pipe, and N upper tower water return pipes;
in this embodiment, any one of the cooling towers is connected to the water suction pool through one water inlet pipe of the water suction pool, the water suction pool is connected to the water supply main pipe through L water outlet pipes of the water suction pool in parallel, and the water supply main pipe supplies water to M 'heat exchangers through M' water supply branch pipes in a one-to-one correspondence manner; a water feeding pump is arranged on the water outlet pipe of the water suction pool in parallel;
in this embodiment, any one of the heat exchangers is correspondingly connected with the water return main pipe through one water return branch pipe, and the water return main pipe is correspondingly connected with the N cooling towers one by one through N upper tower water return pipes which are parallel;
in this embodiment, process medium temperature sensors are arranged in each of the M' heat exchangers, and the process medium temperature sensors are used for acquiring real-time temperature detection values of various process media;
in the embodiment, a cooling feed water temperature sensor and a cooling feed water pressure sensor are arranged on the feed water header pipe, the cooling feed water temperature sensor is used for acquiring a cooling feed water temperature detection value, and the cooling feed water pressure sensor is used for acquiring a cooling feed water pressure detection value;
in this embodiment, a cooling return water temperature sensor and a cooling return water pressure sensor are arranged on the return water main pipe, the cooling return water temperature sensor is used for acquiring a cooling return water temperature detection value, and the cooling return water pressure sensor is used for acquiring a cooling return water pressure detection value;
in this embodiment, N upper tower return pipes are respectively provided with an upper tower valve;
in this embodiment, as can be seen from fig. 2, 3 and 5, the control energy-saving system includes a process medium temperature least-benefit selector, a process medium temperature controller and a supply-return water pressure difference inner loop PID controller; the process medium temperature minimum point selector obtains a process medium temperature minimum point according to a process medium temperature deviation value sequence collected by the circulating cooling water system;
the process medium temperature controller is used for obtaining a cooling water supply and return water pressure difference set value by combining field operation data according to the least unfavorable point of the process medium temperature;
obtaining a cooling water supply pressure detection value and a cooling water return pressure detection value which are obtained by the detection of the circulating cooling water system, and then obtaining a cooling water supply water pressure difference detection value; the difference between the cooling water supply and return pressure difference detection value and the cooling water supply and return pressure difference set value is obtained to obtain a cooling water supply and return pressure difference deviation value; and the water supply and return pressure difference inner ring PID controller cools the water supply and return pressure difference deviation value to adjust and control the opening of all the upper tower valves, so that the outlet water flow of the water supply pump unit is changed.
Referring to fig. 3, a schematic diagram of a heat exchanger group in a synthesis area is shown, cooling feed water enters each heat exchanger in an industrial production area through a water outlet pipe group, a temperature difference exists between a process medium in the heat exchanger and the cooling feed water, the process medium continuously transfers heat to the cooling feed water to realize self temperature reduction, and the cooling feed water flows out of the heat exchanger group carrying heat exchange quantity and is sent to a water return pipe group to become cooling return water. The whole heat exchange process must ensure that the temperatures of different process media in each heat exchanger always meet the respective production requirements, so that all process media in an industrial production field need to be subjected to real-time temperature detection and control.
Preferably, the process medium temperature deviation value sequence comprises M' process medium temperature deviation values, the process medium temperature deviation values are sequentially calculated according to the serial numbers of the heat exchangers, and the process medium temperature deviation value sequence is obtained according to the serial numbers. The temperature deviation value of any process medium is equal to the difference value between the real-time temperature detection value corresponding to the process medium and the temperature set value corresponding to the process medium;
the process medium temperature deviation change rate sequence comprises M' process medium temperature deviation change rates, wherein any process medium temperature deviation change rate is the ratio of the temperature change value of two adjacent detection time periods corresponding to the process medium to the last detection time period.
In this embodiment, the process media temperature deviation value corresponds to a process media temperature deviation change rate.
Referring to fig. 6, a control method of a circulating cooling water minimum differential pressure energy-saving control system based on multiple temperature targets of a process medium is characterized by comprising the following specific steps:
s1: setting a sampling period, and carrying out field operation data acquisition on the circulating cooling water system;
the circulating cooling water system collects real-time temperature detection values of M 'process media in M' heat exchangers on site through the process media temperature sensor;
the circulating cooling water system collects a cooling water supply temperature detection value on site through a cooling water supply temperature sensor;
the circulating cooling water system collects a cooling water supply pressure detection value on site through a cooling water supply pressure sensor;
the circulating cooling water system collects a cooling return water temperature detection value on site through a cooling return water temperature sensor;
the circulating cooling water system collects a cooling backwater pressure detection value on site through a cooling backwater pressure sensor;
s2: obtaining M 'process medium temperature deviation values and corresponding M' process medium temperature deviation change rates according to M 'process medium real-time temperature detection values acquired in the step S1 on site and M' process medium temperature set values; forming a process medium temperature deviation value sequence by the M' process medium temperature deviation values; m' process medium temperature deviation change rate values form a process medium temperature deviation change rate sequence;
s3: the process medium temperature least-benefit point selector selects the temperature deviation minimum value in the process medium temperature deviation value sequence as the process medium temperature least-benefit point, and obtains the process medium temperature deviation change rate corresponding to the process medium temperature least-benefit point and the corresponding heat exchanger;
s4: the process medium temperature controller acquires a process medium temperature deviation value and a process medium temperature deviation change rate of a corresponding heat exchanger according to the least favorable point of the process medium temperature, and inputs field operation data of the corresponding heat exchanger into a corresponding process medium temperature prediction model to obtain a cooling water supply and return water differential pressure set value and a process medium temperature predicted value;
s5: the cooling water supply pressure detection value and the cooling water return pressure detection value are subjected to difference to obtain a cooling water supply and return pressure difference detection value; and the difference between the set value of the cooling water supply and return pressure difference and the detected value of the cooling water supply and return pressure difference is used for obtaining a cooling water supply and return pressure difference deviation value, and the cooling water supply and return pressure difference deviation value is sent to the water supply and return pressure difference inner ring PID controller to adjust the opening of the upper tower valve, so that the outlet water flow of the water supply pump unit is changed. In the present embodiment, the feed pump is a centrifugal pump, and the centrifugal pump is connected to the frequency converter. In step S5, the supply and return water pressure difference inner ring PID controller adjusts the opening of the upper tower valve, so that the specific control content of changing the outlet water flow of the water supply pump set is as follows:
if the deviation value P of the cooling supply and return water differential pressure is greater than 0, the supply and return water differential pressure inner ring PID controller sends an upper tower valve amplification control signal for controlling the upper tower valve to increase the valve opening; the water supply and return pressure difference inner ring PID controller also sends a water supply pump unit outlet water flow reduction control signal for reducing the water pump unit outlet water flow; if the cooling water supply and return pressure difference deviation value P is equal to 0, the water supply and return pressure difference inner ring PID controller keeps the original control state; if the deviation value P of the cooling supply and return water pressure difference is less than 0, the supply and return water pressure difference inner ring PID controller sends a control signal for adjusting the upper tower valve to be small, and the control signal is used for controlling the upper tower valve to reduce the opening of the valve; and the supply and return water pressure difference inner ring PID controller also sends a water supply pump unit outlet water flow increasing control signal for increasing the water pump unit outlet water flow. As can be seen from the first table, the heat exchangers of the circulating cooling water system of the chemical plant A are listed, wherein DSC computer screen temperature display points are arranged in the first table after heat exchange of the process medium.
According to the requirements of the heat exchangers listed in the table, the design indexes and the control ranges of the process media in each heat exchanger are different, so that the outlet temperature of each process medium is detected one by one to realize the overall control of the temperature of the process media in order to ensure that all the process media are maintained in the respective normal production temperature range. The heat exchangers on the industrial production site are huge in quantity, the process media are various in types, and the difference value between the temperature detection value of each process media and the corresponding temperature set value is calculated on line to form a real-time updated process media temperature deviation sequence. If the method of calculating the cooling demand of each process medium and then determining the adjustment amount of the upper tower valve is adopted, the control process is complex and has large redundancy, the adjustment delay of the cooling water flow due to large calculation amount can be caused, and the control energy-saving effect is influenced.
Therefore, the invention provides the most unfavorable point for selecting and controlling the temperature of the process medium in real time, wherein the most unfavorable point of the temperature of the process medium refers to the process medium corresponding to the minimum value in the difference sequence of the detected temperature and the set temperature of each process medium, namely the process medium with the highest possibility of exceeding the set temperature range in each heat exchanger is controlled in real time, and the temperature of the whole process medium is equivalently controlled. With the continuous change of the temperature of the process medium along with the proceeding of industrial production, the temperature deviation and the deviation change rate of the process medium at the most unfavorable point are updated in real time, so that the temperature deviation and the deviation change rate of the process medium at the most unfavorable point are updated in real time as the input of the process medium temperature controller.
As can be seen from fig. 4, the specific steps of finding the minimum value in the sequence of the process medium temperature deviation values as the process medium temperature least favorable point by using the process medium temperature least favorable point selector in step S3 are as follows:
s31: initializing, setting M 'technological medium temperature deviation values to form a difference group, totaling M' technological medium temperature deviation values, and enabling W to bek=M’;k=1;
S32: let Wk+1W is equal to Wk + Xk+1Can be divided by M' evenly, X is a gap with large filling difference; and X is equal to 0-M' -1; s33: calculating Wk+2=Wk+1/M’
S34: from Wk+2In the group, the minimum value is found out from the M' process medium temperature deviation values of each group by adopting a cross comparison method to obtain Wk+2A process medium temperature deviation value;
s35: judgment of Wk+2Whether it is equal to 1; if so, taking the process medium temperature deviation value as the least unfavorable point of the process medium temperature; otherwise, let k be k + 2; return is made to step S32.
Further, as can be seen from fig. 7, in step S4, the step of establishing the deep learning neural network of any heat exchanger is: s411: historical data generated by the operation of the circulating cooling water system in X sampling periods are used as training data of a process medium temperature prediction model in the heat exchanger after the obtained historical data are screened according to screening conditions; wherein X is a positive integer, and in this embodiment, X is equal to 1000.
S412: taking the process medium temperature deviation, the process medium temperature deviation change rate, the cooling water supply temperature detection value and the cooling water supply and return water pressure difference detection value in the historical data of the heat exchanger I as input data according to the determined characteristic variables from the historical data, carrying out data normalization processing to obtain a normalized data set, and dividing the normalized data set into a training sample set and a test sample set;
s413: performing greedy unsupervised pre-training on the training sample set layer by layer based on a stacked automatic encoder to obtain a weight matrix W of an input layer and a hidden layer initialized by a deep learning neural network and a threshold matrix B of the input layer and the hidden layer;
s414: and (3) fine adjustment of parameters: finely adjusting a weight matrix W of an input layer and a hidden layer initialized by a deep learning neural network and a threshold matrix B of the input layer and the hidden layer until the iteration times reach the maximum value of the iteration times to obtain an initial process medium temperature prediction model based on a stacking automatic encoder;
s415: and (5) evaluating the initial process medium temperature prediction model obtained in the step (S414) by using the test sample data set to obtain a process medium temperature prediction model based on the stacking automatic encoder.
As can be seen from fig. 7, the characteristic variables include a process medium temperature deviation value ei and a process medium temperature deviation change rate Δ e in any one of the heat exchangersiA real-time temperature detection value Tci of the process medium, a cooling water supply temperature detection value Tgs and a cooling water supply and return pressure difference detection value P in the circulating cooling water systemΔj。
In this embodiment, the screening condition is historical data that the process medium temperature is within a safe and energy-saving temperature value interval, and the safe and energy-saving temperature value interval is within a process medium temperature threshold value interval.
In this embodiment, the difference between the safe and energy-saving temperature value interval and the process medium temperature threshold value interval is 4 ℃.
In step S413, when greedy unsupervised pre-training is performed on the training sample set layer by layer, the training sample set is divided into P groups of small batch training samples, training is performed in sequence, a Dropout technique is adopted, a part of neurons are randomly selected to suspend working, iteration is performed in sequence, training is performed layer by layer, and weight matrix W, input layer threshold matrix B and hidden layer threshold matrix B of the input layer and hidden layer initialized by the neural network are deeply learned. In step S414, when the parameters are fine-tuned, a top-down small-batch RMSProp optimization method is used to perform fine-tuning on the weight matrix W of the input layer and the hidden layer initialized by the deep learning neural network, and the threshold matrix B of the input layer and the hidden layer.
The pre-control algorithm is directly oriented to field measurement data of an industrial circulating cooling water system, and nonlinear mapping relations among the measurement data of each heat exchanger, including temperature deviation and temperature deviation change rate of process media in each heat exchanger, a cooling water supply temperature detection value, a cooling water supply and return water differential pressure set value and control quantity prediction quantity, namely the cooling water supply and return water differential pressure set value are respectively established through a deep-layer framework. A two-stage off-line training learning method of pre-training-parameter fine tuning is adopted, and a Dropout technology and a RMSPROP technology are introduced to optimize the process medium temperature model parameters in each heat exchanger. The trained model can extract high-order characteristics beneficial to the process medium temperature prediction control effect by means of a hidden mode of deep structure mining data. In addition, the method can improve the generalization capability of the model through unsupervised training of a large number of unlabeled samples.
The number of hidden layers of the deep learning neural network trained by the process medium in the constructed SAE heat exchanger and the number of neurons of each hidden layer have certain influence on the evaluation precision and the off-line training time. Taking the field detection data of the circulating cooling water system as sample input data, and setting the number of neurons in the hidden layer by layer: the optimal number of the neurons of the hidden layer at the layer 1 is determined and fixed, then a layer is added to determine the optimal number of the neurons of the hidden layer at the layer 2, and the like is repeated until the average percentage error (MAPE) is not increased any more.
In order to make the technical solution of the present invention clearer, the principle of the stack automatic encoder used in the present invention is explained. SAE is a typical deep learning neural network, and its basic constituent unit is an Automatic Encoder (AE), and its network structure is shown in fig. 8. An Automatic Encoder (AE) network structure is shown in fig. 8, and is composed of an encoder and a decoder: the input vector is mapped to a feature vector in the hidden layer by the encoder, and then the feature phasor is reconstructed to the original input vector by the decoder.
When given an input sample set X ═ XiI is more than or equal to 1 and less than or equal to N, wherein N is the total number of samples, and x isiThe ith training sample in the sample set has dimension n. Let H ═ HiI is more than or equal to 1 and less than or equal to N is a hidden layer characteristic vector set, hiThe feature vector corresponding to the ith sample has a dimension M', and the coding relationship between X and H is:
H=sf(WX + B) wherein: w is a weight matrix of the input layer and the hidden layer; b is an input layer and hidden layer threshold matrix; sf is the neuron activation function of the encoder, and a sigM' oid function is usually adopted, which has good feature identification:
sf(z)=1/(1+exp(-z))
in the formula: z is the input vector.
The decoder is the inverse operation of the encoder, and takes the characteristic vector of the hidden layer as the input vector
In order to output the set of vectors,
and an output vector corresponding to the ith sample is obtained, and the dimension is n, so that the expression of the decoder is as follows:
in the formula: w' is a weight matrix of the hidden layer and the output layer; b' is a threshold matrix of a hidden layer and an output layer; sg is the neuron activation function of the decoder.
The automatic encoder achieves the purpose of feature extraction by minimizing the reconstruction error between the output vector and the input vector, and the formula of the reconstruction error is as follows:
and continuously adjusting the network weight and the threshold by using a gradient descent algorithm to reduce the reconstruction error, wherein the formula is as follows:
in the formula: l is the learning rate;
to represent
Calculating the deviation of the weight W;
to represent
The bias is calculated for the threshold B.
The SAE is a deep learning neural network model composed of stacked AE, and the output of the lower AE is used as the input of the upper AE. Gradual abstraction of features is achieved by stacking of AEs, eventually leading to more compact, useful features, as shown in fig. 9 for a stacked self-encoder architecture. WnIs a weight matrix of the n-1 st hidden layer and the n-th hidden layer, BnThe threshold matrix is the n-1 st hidden layer and the nth hidden layer. The training process comprises two steps of greedy layer-by-layer unsupervised pre-training and supervised fine tuning. Greedy non-supervision pre-training layer by layer obtains the initial weight and the threshold value of the network through layer by layer training, the input of the bottom layer AE is original data, and the output data of the hidden layer is used as the input data of the upper layer AE.
After the layered pre-training is completed, the hidden layers are stacked, and the relationship between the input data and the output data is expressed as follows:
in the formula: f is an activation function, x
iW, B are respectively the network initialization weight and threshold value obtained by pre-training layer by layer,
is the predicted value of the ith sample. Constructing an error loss function of the actual value and the predicted value, wherein the formula is as follows:
in the formula: n is the total number of samples, yiIs the actual value of the ith sample. The whole network weight and the threshold value are finely adjusted through the back propagation from top to bottom, and the error between the predicted value and the actual value is reduced.
The SAE can effectively extract high-order features of data through greedy unsupervised pre-training layer by layer, better approaches complex functions, reduces parameter optimization space, can quickly obtain network parameters, and improves deep feature learning capability of a neural network.
Compared with a shallow machine learning algorithm, when the deep network processes high-dimensional data, the overfitting problem is more easily generated due to the complex network structure of the deep network, so that the generalization capability of the model is limited. Dropout is a mainstream anti-overfitting technique, and its basic idea is: during model training, a part of nodes are randomly selected to be out of operation, the nodes store the weight of the last iteration, and the output is set to be 0. And the selected nodes restore the weight values reserved before in the process of the next iteration, and part of the nodes are randomly selected again to repeat the process. The network structure will change to some extent in each iteration process, and by adopting the Dropout technology, part of neurons selected randomly do not work temporarily, as shown in fig. 10, the coaction among specific nodes is reduced, the dependence of the network output on the states of the specific nodes is reduced, and overfitting is prevented.
In the parameter fine-tuning stage, a top-down small-batch RMSProp optimization method is adopted to fine-tune the weight and the threshold of the network until the iteration number reaches a set value, and the method specifically comprises the following steps:
in step 1, a global learning rate l and a decay rate ρ are set, and an initialization cumulative variable r1 is 0 and r2 is 0.
Step 2, selecting a small batch of data sets containing M' samples from the training set, and calculating the gradient according to an error loss function:
and 3, calculating the accumulated square gradient as shown in the formula (9):
in the formula: as an element-by-element product symbol.
And 4, respectively updating the weight and the threshold parameters:
and 5, stopping operation when the iteration times meet the requirements, and otherwise, returning to the step 2 to continue to execute the calculation.
The process medium temperature prediction control algorithm can effectively represent the process medium temperature deviation and the deviation change rate, and the complex function between the industrial circulating cooling water supply temperature and the cooling water supply pressure set value, quickly and accurately carry out prediction control on the process medium temperature on the industrial production site, and meanwhile, the algorithm has good generalization capability and strong prediction adaptability to process media with different temperature change characteristics. Compared with manual experience calculation and adjustment, the accuracy of process medium temperature control is greatly improved, the calculation time overhead is saved, and industrial production field management personnel can master the process medium temperature change trend in each heat exchanger in real time; the process medium temperature prediction control algorithm comprises the steps of off-line training and on-line application of a process medium temperature prediction model in a heat exchanger based on a stack automatic encoder to give a cooling water supply and return water differential pressure set value of a circulating cooling water system: the stack automatic encoder belongs to a deep learning neural network, and compared with the traditional shallow neural network, the stack automatic encoder effectively solves a series of problems caused by the random initialization of the parameters of the traditional neural network, can effectively mine the implicit relation of each data, and greatly improves the accuracy of the prediction control of the temperature of the process medium in the industrial production field; compared with the traditional stack automatic encoder method, the improved stack automatic encoder algorithm for parameter fine adjustment by utilizing the RMSProp optimization method has better generalization capability and can improve the prediction control precision of the industrial production field process medium temperature. The method is applied to the energy-saving control of the circulating cooling water system in the industrial production field, can quickly and accurately give the set value of the cooling water supply and return water pressure difference, is beneficial to industrial production field management personnel to master the temperature change trend of the process medium in each heat exchanger in real time, and most importantly, reduces the energy consumption of the circulating cooling water system.
It should be noted that the above description is not intended to limit the present invention, and the present invention is not limited to the above examples, and those skilled in the art may make variations, modifications, additions or substitutions within the spirit and scope of the present invention.