CN114970375B - Rectification process monitoring method based on real-time sampling data - Google Patents

Rectification process monitoring method based on real-time sampling data Download PDF

Info

Publication number
CN114970375B
CN114970375B CN202210901926.2A CN202210901926A CN114970375B CN 114970375 B CN114970375 B CN 114970375B CN 202210901926 A CN202210901926 A CN 202210901926A CN 114970375 B CN114970375 B CN 114970375B
Authority
CN
China
Prior art keywords
rectification process
state
monitoring
data
conservation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210901926.2A
Other languages
Chinese (zh)
Other versions
CN114970375A (en
Inventor
郑鹏
韦兴鹏
史新玉
巩克乐
郑伟
徐�明
杨昭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong Feiyang Chemical Co ltd
Original Assignee
Shandong Feiyang Chemical Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong Feiyang Chemical Co ltd filed Critical Shandong Feiyang Chemical Co ltd
Priority to CN202210901926.2A priority Critical patent/CN114970375B/en
Publication of CN114970375A publication Critical patent/CN114970375A/en
Application granted granted Critical
Publication of CN114970375B publication Critical patent/CN114970375B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B01PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01DSEPARATION
    • B01D3/00Distillation or related exchange processes in which liquids are contacted with gaseous media, e.g. stripping
    • B01D3/14Fractional distillation or use of a fractionation or rectification column
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24Querying
    • G06F16/245Query processing
    • G06F16/2458Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2462Approximate or statistical queries
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24Querying
    • G06F16/245Query processing
    • G06F16/2458Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2477Temporal data queries
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/04Manufacturing
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Business, Economics & Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Mathematical Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Strategic Management (AREA)
  • Probability & Statistics with Applications (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Computational Linguistics (AREA)
  • Fuzzy Systems (AREA)
  • Algebra (AREA)
  • Evolutionary Computation (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Quality & Reliability (AREA)
  • Chemical & Material Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Computing Systems (AREA)
  • Development Economics (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Medical Informatics (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Manufacturing & Machinery (AREA)
  • Health & Medical Sciences (AREA)

Abstract

The invention relates to the technical field of rectification process monitoring, and provides a rectification process monitoring method based on real-time sampling data, which comprises the following steps: sampling the rectification process in real time, and performing correlation calculation on the sampled data to obtain actual state parameters of the rectification process; constructing a state prediction network to train under each working condition state in the rectification process to obtain network output under each working condition state to form a knowledge base, comparing the output of the real-time sampling data after training with the knowledge base for calculation, and further monitoring the rectification process; and constructing a conservation relation model according to each conservation parameter in the rectification process, and monitoring the rectification process to obtain a monitoring result. The invention solves the technical problems of inaccurate monitoring and overhigh cost in the prior art, and realizes the technical effects of high-precision monitoring and low cost.

Description

Rectification process monitoring method based on real-time sampling data
Technical Field
The invention relates to the technical field of rectification process monitoring, in particular to a rectification process monitoring method based on real-time sampling data.
Background
The rectification operation is an important link in the oil refining and chemical production process. The control process of the rectifying tower directly affects the product quality, yield and energy consumption, so the monitoring of the rectifying tower is highly valued by people.
The rectifying tower is an object with multiple inputs and outputs, and consists of multiple stages of tower plates, the internal mechanism is complex, the control action is correspondingly slow, the incidence relation among parameters is complex, the control requirement is higher, the rectifying process is more effectively controlled by monitoring the data of the rectifying process in real time, the industrial production safety and the production efficiency are greatly improved, and the methods are numerous, the invention patent application number of China is 202111414908.3, and the invention provides a soft measurement method for the online detection of components in a special rectifying process, which mainly comprises the following steps: selecting soft measurement model input variables by adopting principal component analysis random forest combination variables, introducing generalized robust loss functions into a limit gradient algorithm, optimizing loss function hyperparameters by adopting a Bayes method, training by combining the soft measurement model input variables with the limit gradient algorithm to obtain a trained ARXGboost model, and substituting real-time data into the ARXGboost model for online detection.
However, in the process of implementing the technical solution of the invention in the above application embodiments, the inventor of the present application finds that the above technology has at least the technical problems of inaccurate measurement and high cost.
Disclosure of Invention
The invention provides a rectification process monitoring method based on real-time sampling data, solves the technical problems of inaccurate monitoring and high cost in the prior art, and achieves the technical effects of high-precision monitoring and low cost.
The invention provides a rectification process monitoring method based on real-time sampling data, which specifically comprises the following technical scheme:
the rectification process monitoring method based on real-time sampling data comprises the following steps:
s1, sampling in real time in a rectification process, and performing correlation calculation on sampling data to obtain actual state parameters in the rectification process;
s2, constructing a state prediction network to train under each working condition state in the rectification process to obtain network output under each working condition state to form a knowledge base, comparing the output of the real-time sampling data after training with the knowledge base for calculation, and further monitoring the rectification process;
and S3, constructing a conservation relation model according to each conservation parameter in the rectification process, and monitoring the rectification process to obtain a monitoring result.
Further, the step S1 includes:
sampling in real time under the normal operation state of the rectifying tower to obtain sampling dataX
Figure 271149DEST_PATH_IMAGE002
NWhich represents the number of samples to be sampled,
Figure 669901DEST_PATH_IMAGE004
to representtThe data vector of the rectifying tower collected at the moment,
Figure 535963DEST_PATH_IMAGE006
and each data vector comprises a rectification process state parameter obtained by directly sampling and consulting data, a rectification process state parameter obtained by directly obtaining parameters through a calculation formula, and a more accurate state parameter obtained by calculating the concentration of each component in the rectification process through an optimization formula, and further samples a data matrix to provide a more accurate parameter basis for monitoring the rectification process.
Further, the step S2 includes:
by defining the set of conditions occurring during the rectificationConAnd concentrating the working conditions into a state parameter set under each working condition
Figure 246430DEST_PATH_IMAGE008
Training the state prediction network sequentially by defining the optimized connection weight vector to obtain the final state prediction network training result
Figure 131341DEST_PATH_IMAGE010
(ii) a And all training results are stored in a knowledge base and used for providing reference basis with the actual working state of the rectification process.
Further, the step S2 includes:
a set of rectification process state vector data is sampled in real time
Figure 742189DEST_PATH_IMAGE012
Obtaining an output state vector through a state prediction network
Figure 687142DEST_PATH_IMAGE014
(ii) a Calculating output vectors of all working condition state sets stored in a knowledge base by defining an angle approximation function Sa
Figure 55544DEST_PATH_IMAGE016
Approximate angle set from actual sampled data
Figure 212987DEST_PATH_IMAGE018
Further, if
Figure 914227DEST_PATH_IMAGE020
First in the setsIndividual element calculation result
Figure 122092DEST_PATH_IMAGE022
Less than 15 DEG, the sampled data is considered as the secondsAnd the data state under different working conditions realizes the monitoring of the rectification process.
Further, the step S3 includes:
by constructing a balance model, calculating each large conservation law followed in the rectification process, and ensuring the normal operation of the rectification process; the specific model is constructed as follows:
Figure 541572DEST_PATH_IMAGE024
wherein,Xa set of sampled data is represented as,Fa set of sampled data extracts is represented,Wthe material conservation relation matrix is expressed,Hotthe heat conservation relation matrix is expressed,Qthe matrix is expressed in the conservation of energy relationship,Othrepresent the other array of conservation relationships,Ysa constraint matrix is shown to define conditions,Outrepresenting the output result of the model; aggregating sampled data using a sampled data extraction setXExtracting the state parameters of rectification process related to conservation, and then making them be related to conservation relationWHotQOthCalculating and comparing, constraining under the limited condition to obtain whether each conservation relation is satisfied, and obtaining the output resultOut
The invention has at least the following technical effects or advantages:
1. according to the invention, the state parameter of the rectification process is acquired more accurately by defining the calculation of the state parameter of the rectification process, so that a data basis is provided for accurate monitoring, the monitoring accuracy is further improved, and the monitoring cost is reduced.
2. The method obtains the output vectors of the output layers under different working conditions by carrying out network training on the state parameters of the rectification process under various working conditions through the state prediction network, basically comprises all the working conditions possibly occurring in the rectification process, has reference comprehensiveness, and more accurately monitors the rectification process.
3. According to the invention, by optimizing the connection weight vector array in the state prediction network, each working condition in the rectification process is trained more quickly and accurately, and the monitoring cost is reduced.
4. The invention calculates the approximate angles of the real-time sampling data and the elements in the state sets of various working conditions of the knowledge base by defining the angle approximate function, and compares the approximate angles to obtain the working condition states corresponding to the sampling data, thereby more accurately realizing the monitoring of the rectification process.
5. According to the invention, by constructing the conservation relation model, the rectification process is monitored according to the multi-aspect conservation law and specific limiting conditions, the accuracy of the rectification process monitoring is ensured, the waste of materials in the rectification process is avoided, and the monitoring cost is further reduced.
Drawings
FIG. 1 is a flow diagram of a rectification process according to the present invention;
FIG. 2 is a step diagram of a rectification process monitoring method based on real-time sampling data according to the present invention.
In the figure: 1-feed preheater, 2-tower bottom product cooler, 3-condenser, 4-tower top product cooler, 5-tower plate, 6-rectifying tower and 7-reboiler.
Detailed Description
The embodiment of the application provides a rectification process monitoring method based on real-time sampling data, solves the technical problems of inaccurate monitoring and overhigh cost in the prior art, and comprises the following specific steps:
firstly, sampling in real time in a rectification process, and carrying out correlation calculation on sampling data to obtain actual state parameters in the rectification process; constructing a state prediction network for defining an optimized connection weight vector to train under each working condition state in the rectification process to obtain network output under each working condition state to form a knowledge base, comparing the output of the real-time sampling data after training with the knowledge base for calculation, and further monitoring the rectification process; and constructing a conservation relation model according to each conservation parameter in the rectification process, and monitoring the rectification process to obtain a monitoring result. The state parameters of the rectification process are acquired more accurately by defining the calculation of the state parameters of the rectification process, so that a data basis is provided for accurate monitoring, the monitoring accuracy is further improved, and the monitoring cost is reduced; output layer output vectors under different working conditions are obtained by carrying out network training on the rectification process state parameters under various working conditions through a state prediction network, all the working conditions possibly occurring in the rectification process are basically contained, the method has reference comprehensiveness, and the rectification process is more accurately monitored; by optimizing the connection weight vector array in the state prediction network, all working conditions in the rectification process are trained more quickly and accurately, and the monitoring cost is reduced; the method comprises the steps of calculating the approximate angles of real-time sampling data and elements in various working condition state sets of a knowledge base by defining an angle approximate function, comparing the approximate angles to obtain the working condition states corresponding to the sampling data, and more accurately monitoring the rectification process; by constructing a conservation relation model, the rectification process is monitored according to the multi-aspect conservation law and specific limiting conditions, the monitoring accuracy of the rectification process is ensured, the waste of materials in the rectification process is avoided, and the monitoring cost is further reduced.
In order to better understand the technical solution, the technical solution will be described in detail with reference to the drawings and the specific embodiments.
Referring to FIG. 1, the rectification process is generally carried out in a rectification column, wherein the feed is fed from the middle of the column, and under the action of a feed preheater, the column section above the feed plate is a rectification section for increasing the concentration of volatile components in the rising gas phase plate by plate, and the column plates inside and below the feed plate are stripping sections for extracting volatile components in the falling liquid phase plate by plate; wherein the tower plates are the places for gas-liquid two-phase mass transfer and heat transfer, and the gas-liquid two-phase mass transfer is carried out on each tower plate; leading the steam led out from the tower top to a condenser, taking part of condensate as reflux liquid to return to the rectifying tower from the tower top, and taking the rest distillate as a tower top product; the liquid extracted from the tower bottom is partially gasified by a reboiler, the steam rises along the tower, and the rest liquid is used as a tower bottom product.
Referring to the attached figure 2, the rectification process monitoring method based on real-time sampling data comprises the following steps:
s1, sampling in real time in a rectification process, and performing correlation calculation on sampling data to obtain actual state parameters in the rectification process;
s11, storing and calculating the acquired data of the rectification process to obtain relevant state parameters of the rectification process;
under the running state of the rectifying tower, the process is carried outSampling in real time to obtain sampled dataX
Figure 768285DEST_PATH_IMAGE026
NWhich represents the number of samples to be sampled,
Figure 494671DEST_PATH_IMAGE028
to representtThe data vector of the rectifying tower collected at the moment,
Figure 47006DEST_PATH_IMAGE030
wherein each data vector comprises a rectification process state parameter obtained by direct sampling and obtained by utilizing a direct acquisition parameter through a calculation formula:
inlet flow rate
Figure 921159DEST_PATH_IMAGE032
Temperature at inlet
Figure 545038DEST_PATH_IMAGE034
Pressure at the top of the column
Figure 486449DEST_PATH_IMAGE036
Outlet flow rate
Figure 177062DEST_PATH_IMAGE038
Amount of reflux
Figure 243239DEST_PATH_IMAGE040
Component 1 outlet flow
Figure 405230DEST_PATH_IMAGE042
Component
2 outlet flow
Figure 434103DEST_PATH_IMAGE044
… …, component (a)mOutlet flow rate
Figure 797083DEST_PATH_IMAGE046
Liquid outlet flow from the bottom of the column
Figure 944030DEST_PATH_IMAGE048
Amount of liquid refluxed from the bottom of the column
Figure 572630DEST_PATH_IMAGE050
mTemperature of the column trays
Figure 957475DEST_PATH_IMAGE052
mPressure of the column plate
Figure 84831DEST_PATH_IMAGE054
Flow rate of condenser cooling water
Figure 420872DEST_PATH_IMAGE056
Steam flow of reboiler
Figure 924666DEST_PATH_IMAGE058
Concentration of light component output at tower top and tower bottom, and heating quantity of reboiler
Figure 164017DEST_PATH_IMAGE060
Cooling capacity of condenser
Figure 196695DEST_PATH_IMAGE062
Concentration of light component output from each column plateCEnthalpy of feed product
Figure 285612DEST_PATH_IMAGE064
Enthalpy of the overhead product
Figure 593096DEST_PATH_IMAGE066
Enthalpy of bottom product
Figure 421375DEST_PATH_IMAGE068
And other related calculated state parameters to obtain a set of state parameters
Figure 421692DEST_PATH_IMAGE070
Figure 201167DEST_PATH_IMAGE072
MTo representtThe number of the state parameters of the rectification process contained in the data vector of the rectification tower collected at any moment.
In particular, the volatility of each component is obtained by looking up relevant data according to the feeding property
Figure 312343DEST_PATH_IMAGE074
And other related component properties, and further calculating to obtain relative volatility and other related component characteristic parameters;
from the above, it can be known that the data matrix is sampledXCan be expressed as:
Figure 791866DEST_PATH_IMAGE076
s12, further limiting the distillation process state parameters in the step S11;
and calculating to obtain the concentration of the light component output according to an optimized An Tuoyin equation and a gas-liquid phase equilibrium equation:
Figure DEST_PATH_IMAGE078
wherein,
Figure DEST_PATH_IMAGE080
Figure DEST_PATH_IMAGE082
is shown astAt the sampling time ofiThe concentration of the effluent of the tray components in the liquid and vapor phases,
Figure DEST_PATH_IMAGE084
is shown astAt the sampling time ofiThe pressure of the tower plates on the layer,
Figure DEST_PATH_IMAGE086
indicates the temperature of the ith tray at the tth sampling time,
Figure DEST_PATH_IMAGE088
Figure DEST_PATH_IMAGE090
Figure DEST_PATH_IMAGE092
Figure DEST_PATH_IMAGE094
is An Tuo factor constant.
According to the method, the state parameters of the rectification process are acquired more accurately by defining the calculation of the state parameters of the rectification process, so that a data basis is provided for accurate monitoring, the monitoring accuracy is further improved, and the monitoring cost is reduced.
S2, constructing a state prediction network to train under each working condition state in the rectification process to obtain network output under each working condition state to form a knowledge base, comparing the output of the real-time sampling data after training with the knowledge base for calculation, and further monitoring the rectification process;
the operating conditions include normality, feed flow fluctuation, tower bottom temperature sensor failure and other operational problems that can occur during the rectification process.
S21, collecting state parameters under each working condition state in the rectification process;
the invention defines the working condition set occurring in the rectification process as Con,
Figure DEST_PATH_IMAGE096
qthe total number of the working conditions is shown,
Figure DEST_PATH_IMAGE098
is shown askUnder the condition of the various working conditions,
Figure DEST_PATH_IMAGE100
wherein the elements are in each working conditionLAnd (4) setting state parameters of the rectification process.
Defining parameters of the state prediction network:
collecting the state parameters of each working condition in the working condition set
Figure DEST_PATH_IMAGE102
Sequentially used as the input of the state prediction network, the number of the obtained input training vectors is L, and the number of the corresponding obtained network input nodes (the dimension of the input mode vector) is LM(ii) a Defining a network connection weight matrixW
Figure DEST_PATH_IMAGE104
RExpressing the number of neurons in an output layer;
network learning rate:
Figure DEST_PATH_IMAGE106
in general, are
Figure DEST_PATH_IMAGE108
Figure DEST_PATH_IMAGE110
Figure DEST_PATH_IMAGE112
Which represents the initial learning rate of the initial learning,nrepresents the number of learning times, num represents the total number of learning times, and decreases as the learning time increases;
defining a neighborhood
Figure DEST_PATH_IMAGE114
Is an output layer neuronjIs related to the number of nodes it contains, and is followed bynThe increase of (2) is gradually reduced, and the neighborhood is determined by the determined winning neuron nodejA central, area containing several neurons, which area may be of any shape, but is generally uniformly symmetrical, generally a square or circular area,
Figure DEST_PATH_IMAGE116
Figure DEST_PATH_IMAGE118
the expression is to be taken to the whole,
Figure DEST_PATH_IMAGE120
represent
Figure DEST_PATH_IMAGE122
The initial value of (c).
An input layer: for input training vector
Figure DEST_PATH_IMAGE124
Carrying out normalization processing to obtain normalized data vector
Figure DEST_PATH_IMAGE126
To network connection weight matrix
Figure DEST_PATH_IMAGE128
To (1)jNetwork connection weight vector of individual neurons
Figure DEST_PATH_IMAGE130
Normalization is carried out to obtain a normalized network connection weight vector
Figure DEST_PATH_IMAGE132
Compute input layer
Figure DEST_PATH_IMAGE134
Connecting weight vectors to each output neuron node
Figure DEST_PATH_IMAGE136
European distance of
Figure DEST_PATH_IMAGE138
. Note book
Figure DEST_PATH_IMAGE140
Is n time input layeriThe output of each of the plurality of neurons,
Figure DEST_PATH_IMAGE142
representing input neuronsiTo output layer neuronsjThe invention defines the Euclidean distance
Figure DEST_PATH_IMAGE144
Comprises the following steps:
Figure DEST_PATH_IMAGE146
further, calculating a minimum distance output nodejAs the best node, note
Figure DEST_PATH_IMAGE148
The invention adjusts and combines neurons
Figure DEST_PATH_IMAGE150
And the connection weight vectors of all neuron nodes in its neighborhood,
Figure DEST_PATH_IMAGE152
taking the next input vector as a training vector set, normalizing the input vector, and calculating the Euclidean distance again until the training is finished
Figure DEST_PATH_IMAGE154
Obtaining the state prediction network training result
Figure DEST_PATH_IMAGE156
Training the elements in the formed working condition set item by item to obtain the final state prediction network training result
Figure DEST_PATH_IMAGE158
(ii) a And all training results are stored in a knowledge base and used for providing reference basis with the actual working state of the rectification process.
The invention obtains the output layer output vectors under different working conditions by carrying out network training on the rectification process state parameters under various working conditions through the state prediction network, basically comprises all the working conditions possibly occurring in the rectification process, has reference comprehensiveness, and more accurately monitors the rectification process.
According to the invention, by optimizing the connection weight vector array in the state prediction network, each working condition in the rectification process is trained more quickly and accurately, and the monitoring cost is reduced.
S22, calculating and monitoring the real-time sampling data according to the output result of the state prediction network;
real-time sampling of a set of distillation process state vector data
Figure DEST_PATH_IMAGE160
Obtaining an output state vector through a state prediction network
Figure DEST_PATH_IMAGE162
. The invention calculates the output vector of each working condition state set stored in a knowledge base by defining an angle approximate function Sa
Figure DEST_PATH_IMAGE164
The set of approximate angles to the actual sampled data is:
Figure DEST_PATH_IMAGE166
then
Figure DEST_PATH_IMAGE168
Further, if
Figure DEST_PATH_IMAGE170
First in a setsIndividual element calculation result
Figure DEST_PATH_IMAGE172
Less than 15 DEG, the sampled data is considered as the secondsThe state of the data under the seed condition,
Figure DEST_PATH_IMAGE174
the invention calculates the approximate angles of the real-time sampling data and the elements in the state sets of various working conditions of the knowledge base by defining the angle approximate function, and compares the approximate angles to obtain the working condition states corresponding to the sampling data, thereby more accurately realizing the monitoring of the rectification process.
And S3, comparing all conservation parameters in the rectification process to construct a conservation relation model, and monitoring the rectification process to obtain a monitoring result.
The invention calculates each large conservation law followed in the rectification process by constructing a conservation relation model, and ensures the normal operation of the rectification process. The specific model is constructed as follows:
Figure DEST_PATH_IMAGE176
wherein,Xa set of sampled data is represented as,Frepresenting a set of sampled data extractions,Wthe material conservation relation matrix is expressed,Hotthe heat conservation relation matrix is expressed,Qthe matrix is expressed in the conservation of energy relationship,Othrepresent the other array of conservation relationships,Ysa constraint matrix of a defined condition is represented,Outrepresenting the output result of the model;
in particular, the qualifying constraint matrixYsThe settings are made by the staff according to the specific rectification.
In particular, the conservation relation matrix comprises conservation relations existing at any stage in the rectification process, such as element inclusion, total material conservation relation, material conservation relation of light components and other related material conservation relations in the material conservation relation matrix;
aggregating sampled data using a sampled data extraction setXExtracting the state parameters of rectification process related to conservation, and then making them be related to conservation relation matrixWHotQOthCalculating and comparing, constraining under the limited condition to obtain whether each conservation relation is satisfied, and obtaining the output resultOut. The specific process is as follows:
Figure DEST_PATH_IMAGE178
according to the invention, by constructing the conservation relation model, the rectification process is monitored according to the multi-aspect conservation law and specific limiting conditions, the accuracy of the rectification process monitoring is ensured, the waste of materials in the rectification process is avoided, and the monitoring cost is further reduced.
In conclusion, the rectification process monitoring method based on real-time sampling data is completed.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While preferred embodiments of the present invention have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. Therefore, it is intended that the appended claims be interpreted as including preferred embodiments and all such alterations and modifications as fall within the scope of the invention.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.

Claims (3)

1. A rectification process monitoring method based on real-time sampling data is characterized by comprising the following steps:
s1, sampling in real time in a rectification process, and performing correlation calculation on sampling data to obtain actual state parameters in the rectification process;
s2, constructing a state prediction network to train under each working condition state in the rectification process to obtain network output under each working condition state to form a knowledge base, comparing the output of the real-time sampling data after training with the knowledge base for calculation, and further monitoring the rectification process;
s3, constructing a conservation relation model according to each conservation parameter in the rectification process, and monitoring the rectification process to obtain a monitoring result;
the step S2 includes:
by defining the set of conditions occurring during the rectificationConAnd concentrating the working conditions into a state parameter set under each working condition
Figure 785985DEST_PATH_IMAGE002
Training the state prediction network sequentially by defining the optimized connection weight vector to obtain the final state prediction network training result
Figure 86254DEST_PATH_IMAGE004
(ii) a All training results are stored in a knowledge base and used for providing reference basis with the actual working state of the rectification process;
a set of rectification process state vector data is sampled in real time
Figure 41572DEST_PATH_IMAGE006
Obtaining an output state vector through a state prediction network
Figure 801718DEST_PATH_IMAGE008
(ii) a Calculating and storing in a knowledge base by defining an angle approximation function SaFinal state prediction network training results of each working condition state set
Figure 239390DEST_PATH_IMAGE010
Approximate angle set from actual sampled data
Figure 375973DEST_PATH_IMAGE012
(ii) a If it is
Figure 982535DEST_PATH_IMAGE014
First in a setsAn element
Figure DEST_PATH_IMAGE015
Is less than 15 deg., the sampled data is considered to be the firstsAnd the data state under different working conditions realizes the monitoring of the rectification process.
2. The rectification process monitoring method based on real-time sampling data according to claim 1, wherein the step S1 comprises the following steps:
sampling in real time under the normal operation state of the rectifying tower to obtain sampling dataX
Figure DEST_PATH_IMAGE017
NWhich represents the number of samples to be sampled,
Figure DEST_PATH_IMAGE019
to representtThe data vector of the rectifying tower collected at the moment,
Figure DEST_PATH_IMAGE021
wherein each data vector comprises a direct sampling acquisition, a consulting data acquisition, a rectification process state parameter obtained by utilizing a direct acquisition parameter through a calculation formula, and a more accurate state parameter obtained by utilizing an optimization formula to calculate the concentration of each component in the rectification process.
3. The method for monitoring the rectification process based on the real-time sampling data as claimed in claim 1, wherein the step S3 comprises:
by constructing an equilibrium model, calculating each large conservation law followed in the rectification process, and ensuring the normal operation of the rectification process; the specific model is constructed as follows:
Figure DEST_PATH_IMAGE023
wherein,Xa set of sampled data is represented as,Frepresenting a set of sampled data extractions,Wthe material conservation relation matrix is expressed,Hotthe heat conservation relation matrix is expressed,Qthe matrix is expressed in the conservation of energy relationship,Othrepresent the other array of conservation relationships,Ysa constraint matrix is shown to define conditions,Outrepresenting the output result of the model; aggregating sampled data using a sampled data extraction setXExtracting the state parameters of rectification process related to conservation, and then making them be related to conservation relation matrixWHotQOthCalculating and comparing, constraining under the limited condition to obtain whether each conservation relation is satisfied, and obtaining the output resultOut
CN202210901926.2A 2022-07-29 2022-07-29 Rectification process monitoring method based on real-time sampling data Active CN114970375B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210901926.2A CN114970375B (en) 2022-07-29 2022-07-29 Rectification process monitoring method based on real-time sampling data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210901926.2A CN114970375B (en) 2022-07-29 2022-07-29 Rectification process monitoring method based on real-time sampling data

Publications (2)

Publication Number Publication Date
CN114970375A CN114970375A (en) 2022-08-30
CN114970375B true CN114970375B (en) 2022-11-04

Family

ID=82968541

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210901926.2A Active CN114970375B (en) 2022-07-29 2022-07-29 Rectification process monitoring method based on real-time sampling data

Country Status (1)

Country Link
CN (1) CN114970375B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117385106A (en) * 2023-12-12 2024-01-12 山东飞扬化工有限公司 Rapid calculation method of feeding amount, feeding method and feeding system

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1693297A (en) * 2005-03-16 2005-11-09 扬子石油化工股份有限公司 Process for advanced controlling rectifying apparatus of butadiene
CN102053613A (en) * 2010-12-31 2011-05-11 北京三博中自科技有限公司 Energy real-time monitoring system of industrial process equipment and monitoring method thereof
CN102183892A (en) * 2011-05-10 2011-09-14 上海交通大学 Load change energy consumption optimizing control method of three-column methanol distillation system
CN107096252A (en) * 2017-05-04 2017-08-29 万华化学集团股份有限公司 The method that rectifying column tower top cold is automatically controlled
CN108287474A (en) * 2017-12-27 2018-07-17 上海交通大学 Based on the probabilistic catalytic reforming reactor robust operation optimization method of raw material
CN110222372A (en) * 2019-05-08 2019-09-10 中国水利水电科学研究院 A kind of Flow of River water quality real-time predicting method and device based on data assimilation
CN111914889A (en) * 2020-06-13 2020-11-10 宁波大学 Rectifying tower abnormal state identification method based on brief kernel principal component analysis
CN111905396A (en) * 2020-06-13 2020-11-10 宁波大学 Rectification process real-time monitoring method based on online sampling data driving
CN112966399A (en) * 2021-04-15 2021-06-15 苏州大学张家港工业技术研究院 Pulse tube refrigerator working condition prediction method and system based on machine learning
CN112989920A (en) * 2020-12-28 2021-06-18 华东理工大学 Electroencephalogram emotion classification system based on frame-level feature distillation neural network
CN113051828A (en) * 2021-03-30 2021-06-29 重庆大学 Online prediction method for natural gas water dew point driven by technological parameters
CN113433906A (en) * 2021-06-24 2021-09-24 之江实验室 Method for product prediction and distillation operation parameter optimization of distillation device
CN113537400A (en) * 2021-09-14 2021-10-22 浙江捷瑞电力科技有限公司 Branch neural network-based edge computing node allocation and exit method
CN113723686A (en) * 2021-08-31 2021-11-30 江南大学 Multitask ash box prediction method and system for energy consumption in organic silicon monomer fractionation process
CN114240243A (en) * 2021-12-30 2022-03-25 无锡雪浪数制科技有限公司 Rectifying tower product quality prediction method and device based on dynamic system identification
CN114239430A (en) * 2021-12-03 2022-03-25 东南大学 Method and system for predicting NOx at furnace outlet based on numerical simulation
CN114460848A (en) * 2022-02-08 2022-05-10 万华化学集团股份有限公司 Method and device for controlling sensitive plate of rectifying tower and computer readable storage medium
CN114741772A (en) * 2022-05-09 2022-07-12 北京化工大学 Method for solving design problem of partition wall rectifying tower based on heuristic algorithm

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7292899B2 (en) * 2005-08-15 2007-11-06 Praxair Technology, Inc. Model predictive control having application to distillation
WO2008055209A2 (en) * 2006-10-31 2008-05-08 Pavilion Technologies, Inc. Integrated model predictive control of distillation and dehydration sub-processes in a biofuel production process
CN100490930C (en) * 2006-12-26 2009-05-27 浙江大学 Method and system for controlling high-purity rectification of rectifying tower based on generalized prediction control
US9014858B2 (en) * 2008-09-30 2015-04-21 Rockwell Automation Technologies, Inc. Energy optimizer for dehydrating biofuels through distillation towers and molecular sieves
CN109669017B (en) * 2017-10-17 2021-04-27 中国石油化工股份有限公司 Refinery distillation tower top cut water ion concentration prediction method based on deep learning
CN116249991A (en) * 2020-07-24 2023-06-09 华为技术有限公司 Neural network distillation method and device
CN112633406A (en) * 2020-12-31 2021-04-09 天津大学 Knowledge distillation-based few-sample target detection method

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1693297A (en) * 2005-03-16 2005-11-09 扬子石油化工股份有限公司 Process for advanced controlling rectifying apparatus of butadiene
CN102053613A (en) * 2010-12-31 2011-05-11 北京三博中自科技有限公司 Energy real-time monitoring system of industrial process equipment and monitoring method thereof
CN102183892A (en) * 2011-05-10 2011-09-14 上海交通大学 Load change energy consumption optimizing control method of three-column methanol distillation system
CN107096252A (en) * 2017-05-04 2017-08-29 万华化学集团股份有限公司 The method that rectifying column tower top cold is automatically controlled
CN108287474A (en) * 2017-12-27 2018-07-17 上海交通大学 Based on the probabilistic catalytic reforming reactor robust operation optimization method of raw material
CN110222372A (en) * 2019-05-08 2019-09-10 中国水利水电科学研究院 A kind of Flow of River water quality real-time predicting method and device based on data assimilation
CN111914889A (en) * 2020-06-13 2020-11-10 宁波大学 Rectifying tower abnormal state identification method based on brief kernel principal component analysis
CN111905396A (en) * 2020-06-13 2020-11-10 宁波大学 Rectification process real-time monitoring method based on online sampling data driving
CN112989920A (en) * 2020-12-28 2021-06-18 华东理工大学 Electroencephalogram emotion classification system based on frame-level feature distillation neural network
CN113051828A (en) * 2021-03-30 2021-06-29 重庆大学 Online prediction method for natural gas water dew point driven by technological parameters
CN112966399A (en) * 2021-04-15 2021-06-15 苏州大学张家港工业技术研究院 Pulse tube refrigerator working condition prediction method and system based on machine learning
CN113433906A (en) * 2021-06-24 2021-09-24 之江实验室 Method for product prediction and distillation operation parameter optimization of distillation device
CN113723686A (en) * 2021-08-31 2021-11-30 江南大学 Multitask ash box prediction method and system for energy consumption in organic silicon monomer fractionation process
CN113537400A (en) * 2021-09-14 2021-10-22 浙江捷瑞电力科技有限公司 Branch neural network-based edge computing node allocation and exit method
CN114239430A (en) * 2021-12-03 2022-03-25 东南大学 Method and system for predicting NOx at furnace outlet based on numerical simulation
CN114240243A (en) * 2021-12-30 2022-03-25 无锡雪浪数制科技有限公司 Rectifying tower product quality prediction method and device based on dynamic system identification
CN114460848A (en) * 2022-02-08 2022-05-10 万华化学集团股份有限公司 Method and device for controlling sensitive plate of rectifying tower and computer readable storage medium
CN114741772A (en) * 2022-05-09 2022-07-12 北京化工大学 Method for solving design problem of partition wall rectifying tower based on heuristic algorithm

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于深度学习的精馏过程质量参数软测量方法研究;卢政印;《中国优秀硕士学位论文全文数据库》;20220115;全文 *
工程化精馏实训装置操作条件优化研究;王卫霞等;《黑龙江科技信息》;20141125(第33期);全文 *
精馏塔的故障诊断方案设计;由宏君;《贵州化工》;20040430(第02期);全文 *

Also Published As

Publication number Publication date
CN114970375A (en) 2022-08-30

Similar Documents

Publication Publication Date Title
CN104965967B (en) A kind of yield real-time predicting method of atmospheric and vacuum distillation unit
CN114970375B (en) Rectification process monitoring method based on real-time sampling data
CN108803520B (en) Dynamic process monitoring method based on variable nonlinear autocorrelation rejection
Jorge et al. Simulation and analysis of a sugarcane juice evaporation system
Xie et al. A novel robust data reconciliation method for industrial processes
CN113420500B (en) Intelligent atmospheric and vacuum system
CN109669415A (en) A kind of dynamic process monitoring method based on the analysis of structuring canonical variable
CN109255477A (en) A kind of wind speed forecasting method and its system and unit based on depth limit learning machine
Aguel et al. Parametric study and modeling of cross-flow heat exchanger fouling in phosphoric acid concentration plant using artificial neural network
CN110869728A (en) Method for determining the stress level in a material of a process engineering device
CN112597430A (en) Method for optimizing operation parameters of complex rectifying tower
Abdul-Wahab et al. Optimization of multistage flash desalination process by using a two-level factorial design
CN115014454A (en) Soft measurement method, system, equipment and medium for main steam flow of thermal power generating unit
Al-Shayji et al. Predictive modeling of large-scale commercial water desalination plants: data-based neural network and model-based process simulation
Yao et al. Fault diagnosis based on RseNet-LSTM for industrial process
CN116415394A (en) Two-loop equipment model simulation method integrating machine learning and mechanism model
CN114492551A (en) Soft measurement strategy based on improved GWO-SVM
Xie et al. TSC prediction and dynamic control of BOF steelmaking with state-of-the-art machine learning and deep learning methods
CN116821695A (en) Semi-supervised neural network soft measurement modeling method
Chen et al. Performance prediction and optimization of the air‐cooled condenser in a large‐scale power plant using machine learning
Sánchez Escalona et al. Artificial neural network modeling of hydrogen sulphide gas coolers ensuring extrapolation capability
CN115631804A (en) Method for predicting outlet concentration of sodium aluminate solution in evaporation process based on data coordination
CN105259790A (en) Time registration method for multi-parameter data during alumina production and evaporation process
CN116662925A (en) Industrial process soft measurement method based on weighted sparse neural network
CN107694139B (en) A kind of plate distillation column coefficient of performance on-line analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant