CN116246724A - Zinc rotary kiln key monitoring variable soft measurement method, device, terminal and medium - Google Patents

Zinc rotary kiln key monitoring variable soft measurement method, device, terminal and medium Download PDF

Info

Publication number
CN116246724A
CN116246724A CN202310150189.1A CN202310150189A CN116246724A CN 116246724 A CN116246724 A CN 116246724A CN 202310150189 A CN202310150189 A CN 202310150189A CN 116246724 A CN116246724 A CN 116246724A
Authority
CN
China
Prior art keywords
rotary kiln
kiln
zinc
digital twin
temperature
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.)
Pending
Application number
CN202310150189.1A
Other languages
Chinese (zh)
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.)
Peng Cheng Laboratory
Original Assignee
Peng Cheng Laboratory
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 Peng Cheng Laboratory filed Critical Peng Cheng Laboratory
Priority to CN202310150189.1A priority Critical patent/CN116246724A/en
Publication of CN116246724A publication Critical patent/CN116246724A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/20Identification of molecular entities, parts thereof or of chemical compositions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10Analysis or design of chemical reactions, syntheses or processes

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Manufacture And Refinement Of Metals (AREA)

Abstract

The invention discloses a method, a device, a terminal and a medium for soft measurement of key monitoring variables of a zinc rotary kiln, which comprise the following steps: analyzing a temperature field conservation theory and a chemical reaction process in the zinc rotary kiln, and constructing a mechanism model of coupling of the temperature and the chemical reaction of the zinc rotary kiln; analyzing and mining through preset observation data, and identifying key parameters of the model to obtain digital twin parameters of the rotary kiln consistent with equipment working conditions; aiming at the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln, adopting a variable step length bidirectional finite difference to carry out digital twin simulation solution; according to the change of the equipment regulation and control parameters and the external disturbance data, dynamically updating the digital twin parameters and timely simulating to obtain the temperature in the rotary kiln and the chemical composition in the kiln, so that the real-time feedback of the regulation and control operation is realized. The invention provides a zinc rotary kiln digital twin model construction and virtual-real interaction method, and key monitoring variable soft measurement of a zinc rotary kiln is realized through digital twin.

Description

Zinc rotary kiln key monitoring variable soft measurement method, device, terminal and medium
Technical Field
The invention relates to the technical field of metallurgical equipment modeling and application, in particular to a method, a device, a terminal and a medium for soft measurement of key monitoring variables of a zinc rotary kiln.
Background
The nonferrous metallurgy industry has the characteristics of high energy consumption and high carbon emission. The zinc rotary kiln is used as core equipment of a recovery link in a zinc smelting process, and realizes separation and recovery of zinc in leaching residues through high-temperature combustion and oxidation-reduction reaction after mixing of zinc leaching residues and coke, so that the characteristics of high energy consumption and high carbon emission are more obvious. However, fluctuation of factors such as raw materials, coke grade, environment, slag-to-carbon ratio and the like causes random change of reaction conditions in the rotary kiln. If the reaction condition in the kiln is insufficient, the recovery rate of zinc is low, and resource waste is caused; if the reaction in the kiln is excessive, the consumption of coke is excessive, so that the waste of coke use and the increase of carbon emission are caused, and meanwhile, the potential safety hazard is caused by the burning through of the kiln wall. The optimal reaction condition in the kiln is achieved by workers through regulating and controlling the rotation speed, the feeding amount and the blasting state of the kiln head, so that the minimization of the coke consumption is realized on the premise of ensuring the process index that the zinc recovery rate is more than or equal to 95%. Therefore, to realize the optimal operation of energy saving and carbon reduction of the zinc rotary kiln, key variables reflecting the reaction conditions in the kiln must be monitored in real time to ensure that the interior of the zinc rotary kiln is always in an optimal stable reaction state.
Because of the unique characteristics of long axial dimension, high-temperature sealing and continuous rotation, the key monitoring variables reflecting the reaction conditions in the kiln are difficult to directly obtain, so that workers can only observe the appearance of kiln head flame, kiln tail temperature and kiln slag components by means of eyes to carry out equipment regulation and control, and the zinc rotary kiln has high blindness and randomness. In order to realize the simulation of the reaction condition in the kiln, the traditional zinc rotary kiln modeling simulation method is mostly carried out offline based on theoretical and empirical assumptions, the error of a simulation result is larger due to the imperfection of the data fusion of the established model and equipment, and meanwhile, the simulation result is inconsistent with the field condition due to the fact that the model fails to dynamically track the working condition change of the actual equipment, so that the application value is limited.
For example, the existing method mainly adopts one-dimensional finite difference to simulate the internal components of the rotary kiln through temperature field analysis and reaction analysis based on mass conservation and energy conservation, but is difficult to apply to field environments with changeable actual working conditions due to the fact that the influence of time-varying working condition data is not considered. For another example, a numerical simulation model and a sintering temperature mechanism model are established by adopting a numerical simulation technology based on thermal balance analysis, and error compensation is carried out on the established model through a BP neural network, but complex multi-field coupling and mutual influence in the rotary kiln are not considered, and the established model can only reflect a certain single physical characteristic in the rotary kiln.
The digital twin is based on the update of physical model, history data and real-time data, integrates the simulation process of multiple disciplines, multiple physical quantities and multiple scales, and synchronously maps the virtual reality of the equipment in the digital space, thereby reflecting the actual running condition of the equipment. The effective fusion of digital twinning mechanism and data, the real-time mapping of virtual model and real equipment is the biggest difference from traditional offline modeling simulation. The digital twin model with high precision and dynamic consistency is a necessary means for effectively monitoring the internal reaction condition of the rotary kiln. By constructing a digital twin model for reflecting the reaction condition of the zinc rotary kiln, synchronous mapping is maintained with the on-site rotary kiln in real time, and the defect that the mechanism data is difficult to fuse, virtual and real are difficult to synchronize in the traditional modeling simulation is overcome. The monitoring results of key variables in the kiln, including the temperature of the high-temperature reaction zone, the reaction rate, the zinc recovery progress and the carbon consumption progress, are obtained through digital twin real-time simulation, so that the problem that the reaction condition of the high-temperature zone can only be inferred by indirectly relying on limited observation data is solved, and the regulation and control effect of the rotary kiln is fed back in time.
However, the research shows that no report on a method for realizing soft measurement of key variables of a rotary kiln by performing model construction and verification on the zinc rotary kiln based on a digital twin concept exists yet.
Accordingly, there is a need in the art for improvement.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a method, a device, a terminal and a medium for soft measurement of key monitoring variables of a zinc rotary kiln, so as to realize soft measurement of the key monitoring variables of the zinc rotary kiln through digital twinning.
The technical scheme adopted for solving the technical problems is as follows:
in a first aspect, the invention provides a method for soft measurement of key monitoring variables of a zinc rotary kiln, which comprises the following steps:
analyzing a temperature field conservation theory and a chemical reaction process in the zinc rotary kiln, and constructing a mechanism model of coupling of the temperature and the chemical reaction of the zinc rotary kiln;
analyzing and mining through preset observation data, and identifying key parameters of the model to obtain digital twin parameters of the rotary kiln consistent with equipment working conditions;
aiming at the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln, adopting a variable step length bidirectional finite difference to carry out digital twin simulation solution;
according to the change of the equipment regulation and control parameters and the external disturbance data, dynamically updating the digital twin parameters and timely simulating to obtain the temperature in the rotary kiln and the chemical composition in the kiln, so that the real-time feedback of the regulation and control operation is realized.
In one implementation, the analysis of temperature field conservation theory and chemical reaction processes in the zinc rotary kiln comprises the following steps:
The rotary kiln is divided into a kiln head area, a high-temperature reaction area and a kiln tail area, wherein a complex heat transfer process between materials, gas and kiln walls occurs in the sectional area, and a chemical reaction process between the materials and the gas occurs in the high-temperature reaction area in sections.
In one implementation, the analysis of the conservation theory of the temperature field and the chemical reaction process in the zinc rotary kiln builds a mechanism model of the coupling of the temperature and the chemical reaction of the zinc rotary kiln, and the analysis comprises the following steps:
analyzing convection heat transfer and radiation heat transfer among the flue gas, the kiln wall and the materials, and constructing a heat conservation equation of the flue gas, the materials and the kiln wall;
the high-temperature decomposition of the preheating section and the coke combustion of the high-temperature reaction zone, the reduction of zinc oxide and the reaction of iron compounds mainly occur in the volatilizing kiln;
and (3) based on the principle analysis of chemical reaction by mass conservation and energy conservation, constructing a principle model of the coupling of the temperature of the zinc rotary kiln and the chemical reaction.
In one implementation, the constructing a mechanism model of the zinc rotary kiln temperature and the chemical reaction coupling includes:
temperature field analysis and modeling are independently carried out, and a preliminary result of the temperature field is used as a main calculation basis of the reaction rate in the chemical reaction;
And correcting the temperature field and modeling the coupling through analysis and calculation of enthalpy change in the chemical reaction process.
In one implementation manner, the analyzing and mining through preset observation data, identifying key parameters of a model, and obtaining digital twin parameters of the rotary kiln consistent with equipment working conditions includes:
acquiring and storing temperature, concentration and pressure data in the running process of the rotary kiln in a discrete time sequence data form through a DCS system;
according to the operation mode experience of on-site observation and regulation of the rotary kiln, analyzing and guiding rotary kiln working conditions, and analyzing and processing DCS data and flame image observation data;
and extracting the complement data corresponding to the rotary kiln state through time sequence data analysis and image processing.
In one implementation manner, the extracting the complement data corresponding to the rotary kiln state through time sequence data analysis and image processing includes:
and extracting time domain features, frequency domain features, static features of the flame image and dynamic features of the flame image corresponding to the rotary kiln state through the time sequence data analysis and the image processing.
In one implementation manner, the extracting the complement data corresponding to the rotary kiln state through time sequence data analysis and image processing includes:
And (3) recursively calculating a parameter estimation value once per observation of data, and updating and verifying the digital twin model according to the prediction result.
In one implementation, the digital twin simulation solution is performed by adopting a variable step bidirectional finite difference according to the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln, and the method comprises the following steps:
dividing the rotary kiln into N units along the axial direction, and determining a mass conservation formula existing in each unit;
deducing in a differential mode according to a mass conservation formula in each unit to form a recursion formula;
and carrying out digital twin simulation solving by adopting the variable step length bidirectional finite difference.
In one implementation, the digital twin simulation solution using the variable step bidirectional finite difference includes:
aiming at the problem that the material and gas countercurrent and boundary conditions in the zinc rotary kiln are not in the same position, a kiln head material composition matrix MH1 is defined at the kiln head of the rotary kiln, and a kiln tail material composition matrix MT1 is defined at the kiln tail of the rotary kiln;
obtaining initial conditions of a gas phase in a kiln head material component matrix MH1 according to compressed air and oxygen-enriched drum charging amount of the kiln head, and obtaining initial conditions of a material phase in a kiln tail material component matrix MT1 according to raw material component assay of kiln tail feeding;
And respectively solving a kiln head material component matrix MH1 and a kiln tail material component matrix MT1, updating the material phase in the kiln head material component matrix MHi by the material phase result in the kiln tail material component matrix MTi after one-time solving, and updating the gas phase in the kiln tail material component matrix MTi by the gas phase result in the kiln head material component matrix MHi.
In one implementation manner, the dynamically updating digital twin parameters and timely simulating to obtain the temperature and chemical composition in the rotary kiln according to the change of the equipment regulation parameters and the external disturbance data comprises the following steps:
and according to the equipment regulation and control parameters and the change of the external disturbance data, obtaining variables of the temperature, the reaction rate, the zinc generation progress and the carbon consumption progress in the rotary kiln through dynamically updated digital twin simulation, and realizing real-time feedback of regulation and control operation.
In a second aspect, the invention provides a zinc rotary kiln key monitoring variable soft measurement device, comprising:
the model construction module is used for analyzing the conservation theory of the temperature field and the chemical reaction process in the zinc rotary kiln and constructing a mechanism model of coupling of the temperature and the chemical reaction of the zinc rotary kiln;
the digital twin parameter module is used for analyzing and mining through preset observation data, identifying key parameters of the model and obtaining digital twin parameters of the rotary kiln consistent with the working conditions of equipment;
The simulation solving module is used for carrying out digital twin simulation solving by adopting a variable step length bidirectional finite difference aiming at the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln;
and the dynamic updating module is used for dynamically updating the digital twin parameters according to the change of the equipment regulation and control parameters and the external disturbance data, timely simulating to obtain the temperature and chemical components in the rotary kiln, and realizing the real-time feedback of the regulation and control operation.
In a third aspect, the present invention provides a terminal comprising: the system comprises a processor and a memory, wherein the memory stores a zinc rotary kiln key monitoring variable soft measurement program, and the zinc rotary kiln key monitoring variable soft measurement program is used for realizing the operation of the zinc rotary kiln key monitoring variable soft measurement method according to the first aspect when being executed by the processor.
In a fourth aspect, the present invention further provides a medium, where the medium is a computer readable storage medium, where the medium stores a zinc rotary kiln key monitoring variable soft measurement program, where the zinc rotary kiln key monitoring variable soft measurement program is used to implement the operation of the zinc rotary kiln key monitoring variable soft measurement method according to the first aspect when the zinc rotary kiln key monitoring variable soft measurement program is executed by a processor.
The technical scheme adopted by the invention has the following effects:
According to the method, a digital twin model for reflecting the reaction condition of the zinc rotary kiln is constructed, synchronous mapping is kept with the on-site rotary kiln in real time, and the defect that mechanism data are difficult to fuse, virtual and real are difficult to synchronize in traditional modeling simulation is overcome; the monitoring results of key variables in the kiln, including the temperature of a high-temperature reaction zone, the reaction rate, the zinc recovery progress and the carbon consumption progress, are obtained through digital twin real-time simulation, so that the problem that the reaction condition of the high-temperature zone can only be inferred by indirectly relying on limited observation data is solved, the regulation and control effect of the rotary kiln is fed back in time, and the green intelligent operation of the zinc rotary kiln is serviced; the invention provides a zinc rotary kiln digital twin model construction and virtual-real interaction method, and key monitoring variable soft measurement of a zinc rotary kiln is realized through digital twin.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to the structures shown in these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for soft measurement of key monitored variables of a zinc rotary kiln in one implementation of the invention.
FIG. 2 is a schematic diagram of modeling of temperature and chemical reaction coupling in one implementation of the invention.
FIG. 3 is a schematic diagram of dynamic update of digital twin parameters in one implementation of the invention.
FIG. 4 is a schematic diagram of digital twin key parameter prediction in one implementation of the invention.
FIG. 5 is a schematic diagram of a two-way finite difference solution gas-solid two-phase countercurrent process in one implementation of the present invention.
Fig. 6 is a schematic diagram of the digital twin platform effect of a zinc rotary kiln in one implementation of the invention.
FIG. 7 is a schematic representation of digital twin simulation results in one implementation of the invention.
Fig. 8 is a functional schematic of a terminal in one implementation of the invention.
The achievement of the objects, functional features and advantages of the present invention will be further described with reference to the accompanying drawings, in conjunction with the embodiments.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more clear and clear, the present invention will be further described in detail below with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
Exemplary method
The digital twin is based on the update of physical model, history data and real-time data, integrates the simulation process of multiple disciplines, multiple physical quantities and multiple scales, and synchronously maps the virtual reality of the equipment in the digital space, thereby reflecting the actual running condition of the equipment. The effective fusion of digital twinning mechanism and data, the real-time mapping of virtual model and real equipment is the biggest difference from traditional offline modeling simulation. The digital twin model with high precision and dynamic consistency is a necessary means for effectively monitoring the internal reaction condition of the rotary kiln. By constructing a digital twin model for reflecting the reaction condition of the zinc rotary kiln, synchronous mapping is maintained with the on-site rotary kiln in real time, and the defect that the mechanism data is difficult to fuse, virtual and real are difficult to synchronize in the traditional modeling simulation is overcome. The monitoring results of key variables in the kiln, including the temperature of the high-temperature reaction zone, the reaction rate, the zinc recovery progress and the carbon consumption progress, are obtained through digital twin real-time simulation, so that the problem that the reaction condition of the high-temperature zone can only be inferred by indirectly relying on limited observation data is solved, and the regulation and control effect of the rotary kiln is fed back in time.
However, the method for realizing soft measurement of key variables of the rotary kiln by performing model construction and verification on the zinc rotary kiln based on the digital twin concept is found through retrieval.
Aiming at the technical problems, the embodiment provides a key monitoring variable soft measurement method of a zinc rotary kiln, which is used for solving the defects of difficult mechanism data fusion and difficult virtual-real synchronization in traditional modeling simulation by constructing a digital twin model for reflecting the reaction condition of the zinc rotary kiln and maintaining synchronous mapping with the on-site rotary kiln in real time; and the monitoring results of key variables in the kiln, including the temperature of a high-temperature reaction zone, the reaction rate, the zinc recovery progress and the carbon consumption progress, are obtained through digital twin real-time simulation, so that the key monitoring variables of the zinc rotary kiln are measured in a soft mode.
As shown in fig. 1, the embodiment of the invention provides a key monitoring variable soft measurement method for a zinc rotary kiln, which comprises the following steps:
and S100, analyzing a temperature field conservation theory and a chemical reaction process in the zinc rotary kiln, and constructing a mechanism model of coupling of the temperature and the chemical reaction of the zinc rotary kiln.
In this embodiment, the method for soft measurement of key monitoring variables of a zinc rotary kiln is applied to a terminal, which includes but is not limited to: a computer, etc.
The zinc rotary kiln is important energy consumption and carbon emission equipment in zinc smelting, and the control depending on the experience of workers is blind and underoptimized due to the lack of key monitoring variables in the kiln. The digital twin of virtual-real fusion and virtual-real synchronization is a necessary means for realizing real-time accurate monitoring of key variables of the zinc rotary kiln. In order to overcome the defects of the prior art, the embodiment of the invention provides a method for constructing a digital twin model of a zinc rotary kiln and interacting virtual and real, and realizes soft measurement of key monitoring variables of the zinc rotary kiln through digital twin.
Firstly, a mechanism model of zinc rotary kiln temperature and chemical reaction coupling is constructed by analyzing a temperature field conservation theory and a chemical reaction process in the zinc rotary kiln, key parameters of the model are identified by analyzing and mining limited observation data, and rotary kiln digital twin consistent with equipment working conditions is realized. On the basis, a digital twin simulation solving method based on a variable-step bidirectional Euler method is provided for the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln, and the simulation solving precision of key variables in the kiln is improved. Compared with the data actually collected on the industrial site, the result shows that the simulation solving effect of the zinc rotary kiln digital twin gas-solid two-phase countercurrent process is better.
Therefore, according to the changes of equipment regulation and control parameters (blanking speed, kiln body rotating speed and blasting state) and external disturbance data (environment temperature, coke grade and coke ratio), key variables such as temperature, reaction rate, zinc generation progress, carbon consumption progress and the like in the rotary kiln can be obtained through digital twin timely simulation of dynamic update, real-time feedback of regulation and control operation is realized, so that regulation and control reference is provided for workers, and the problems of blindness and randomness in regulation caused by the loss of the key variables can be effectively solved. Finally, the optimal operation of energy conservation and carbon reduction of the zinc rotary kiln is realized, and the green and intelligent zinc rotary kiln is served for nonferrous metallurgical equipment.
Specifically, in one implementation of the present embodiment, the step S100 includes the following steps before:
and S001, dividing the rotary kiln into a kiln head region, a high-temperature reaction region and a kiln tail region, wherein a complex heat transfer process between materials, gas and kiln walls occurs in the sectional region, and a chemical reaction process between the materials and the gas occurs in the high-temperature reaction region in a sectional manner.
In the embodiment, complicated heat and mass transfer processes such as material transfer, heat transfer, chemical reaction and the like exist in the zinc rotary kiln, wherein a temperature field is an important influencing factor of a plurality of reactions and fields, and a chemical reaction field is a direct embodiment of a final process index of zinc recovery of the zinc rotary kiln, so that digital twin modeling of the temperature field and the chemical reaction field is mainly performed. The rotary kiln is divided into a kiln head area, a high-temperature reaction area and a kiln tail area, wherein a complex heat transfer process among materials, gas and kiln walls occurs in the sectional area, and a chemical reaction process among the materials and the gas occurs in the high-temperature reaction area. The complex chemical reaction in the kiln is accompanied by heat transfer, and mainly comprises convection heat transfer and radiation heat transfer among the smoke, kiln walls and materials. The flue gas exchanges heat with heat radiation and kiln walls and materials through thermal convection, the materials exchange heat with the kiln walls through thermal radiation, and the kiln outer wall exchanges heat with the external environment through radiation and convection.
Specifically, in one implementation of the present embodiment, step S100 includes the steps of:
s101, analyzing convection heat transfer and radiation heat transfer among smoke, kiln walls and materials, and constructing a heat conservation equation of the smoke, the materials and the kiln walls;
step S102, mainly generating pyrolysis of a preheating section and coke combustion of a high-temperature reaction zone, reduction of zinc oxide and reaction of iron compounds in a volatilizing kiln;
and step S103, carrying out mechanism analysis of chemical reaction based on mass conservation and energy conservation, and constructing a mechanism model of the coupling of the temperature of the zinc rotary kiln and the chemical reaction.
In the embodiment, on the basis of analyzing heat transfer of complex chemical reaction in the kiln, a heat conservation equation of smoke, materials and kiln walls is constructed. In addition, reactions such as pyrolysis in a preheating section, coke combustion in a high-temperature reaction zone, reduction of zinc oxide, and reaction of iron compounds mainly occur in the volatilizing kiln, and chemical reaction mechanism analysis is performed based on mass conservation and energy conservation.
Under steady-state conditions, assuming that the heat transfer process at any radial section of the rotary kiln does not take into account time-varying effects, each cell satisfies the energy conservation equation:
Figure SMS_1
Wherein M represents mass flow, cp represents specific heat, Q represents heat change, T represents temperature, x represents kiln axial length, subscripts s and g represent material and gas, subscripts ew and es represent bare kiln wall and bare material, subscripts cw and cs represent blanket kiln wall and blanket material, R and CV represent heat radiation and heat conduction, ΔH represents enthalpy change of chemical reaction, Q c Indicating the coke combustion exotherm. Since there is no energy accumulation in the kiln wall, the energy conservation equation of the kiln wall is:
Q sh =Q g-ew +Q es-ew +Q cs-cw (2)
wherein Q is sh Indicating heat transfer from the kiln walls to the environment. Because the heat release of carbon combustion is the main source of heat in the rotary kiln, the chemical reaction enthalpy change delta H is ignored, and the heat conduction and heat convection coefficients are substituted into the formulas (1) and (2), the method can be further obtainedTo the formula:
Figure SMS_2
wherein A represents heat exchange area, h represents heat exchange coefficient, and i and j represent components of materials and gases. And further carrying out mechanism modeling on the material chemical reaction process in the kiln on the basis of the mechanism analysis modeling of the temperature field of the rotary kiln. The main component in the zinc leaching slag is ZnFe 2 O 4 、ZnSO 4 And ZnO, the pyrolysis reaction which occurs first in the kiln head region mainly includes:
6ZnFe 2 O 4 +C=6ZnO+4Fe 3 O 4 +CO 2
Fe 3 O 4 +CO=3FeO+CO 2
2ZnSO 4 =2ZnO+2SO 2 +O 2 (4)
the main chemical reaction of the high temperature reaction zone in the kiln is divided into two parts, namely the interior of the material layer and the upper part of the material layer. Coke combustion, reduction of zinc oxide, reaction of iron compounds, etc. occur inside the material layer, and specific reactions are as follows formula 6:
C+O 2 =CO 2 (1) CO 2 +C=2CO(2)
ZnO+CO=Zn↑+CO 2 (3) ZnO+C=Zn↑+CO(4)
FeO+CO=Fe+CO 2 (5) ZnO+Fe=Zn↑+FeO(6) (5)
The CO combustion and the oxidation process of Zn simple substance mainly occur above the material layer, and the specific reaction formula is as follows:
2CO+O 2 =2CO 2 (7) 2Zn↑+O 2 =2ZnO(8) (6)
the rate of the 8 chemical reactions is an important factor affecting the zinc oxide production efficiency of a high-temperature reaction zone in the kiln and is also a key parameter of digital twin modeling. Rate modeling of the main reaction process in the kiln was performed using the arrhenius equation for multi-step finite reaction rates:
Figure SMS_3
wherein x and y represent the reactant or the product, c is the concentration of the substance, a and b are the reaction stages, A represents the factor before the reaction, ea represents the activation energy, R represents the molar gas constant, and T represents the absolute temperature. According to the Arrhenius equation of the chemical reaction and the multi-step reaction, the rate formula of 9 reactions is constructed:
Figure SMS_4
according to the material generation and consumption conditions of the reaction in the kiln, constructing a mass conservation formula of solid and gas components in the kiln:
Figure SMS_5
the main factors influencing the recovery efficiency and the carbon consumption of the zinc rotary kiln are an internal temperature field and a chemical reaction process.
Specifically, in one implementation of the present embodiment, step S103 includes the steps of:
step S103a, temperature field analysis and modeling are independently carried out, and a preliminary result of the temperature field is used as a main calculation basis of the reaction rate in the chemical reaction;
Step S103b, correction of a temperature field and coupling modeling are performed through analysis and calculation of enthalpy change in the chemical reaction process.
In this embodiment, in order to implement the coupling analysis of the temperature field and the chemical reaction field in the kiln, based on the premise that the heat generated by burning coke in the kiln accounts for more than 90% of the total heat, the temperature field analysis and modeling are performed independently, the preliminary result of the temperature field is used as the main calculation basis of the reaction rate in the chemical reaction, and then the correction and coupling modeling of the temperature field are performed through the analysis calculation of the enthalpy change Δh in the chemical reaction process, and the specific flow is shown in fig. 2:
modifying the formula (1) and the formula (2) by using the chemical reaction enthalpy, wherein:
ΔH s =r 3 ΔH 3 +r 4 ΔH 4 +r 5 ΔH 5 +r 6 ΔH 6 (10)
ΔH g =r 8 ΔH 8 (11)
according to the embodiment, a mechanism model of zinc rotary kiln temperature and chemical reaction coupling is constructed by analyzing a temperature field conservation theory and a chemical reaction process in the zinc rotary kiln, key parameters of the model are identified by analyzing and mining limited observation data, and rotary kiln digital twin consistent with equipment working conditions is realized.
As shown in fig. 1, in an implementation manner of the embodiment of the invention, the method for soft measurement of key monitoring variables of the zinc rotary kiln further comprises the following steps:
and step S200, analyzing and mining through preset observation data, and identifying key parameters of the model to obtain digital twin parameters of the rotary kiln consistent with the working conditions of equipment.
In this embodiment, because the regulation parameters, the observation data and the disturbance data change in real time in the running process of the zinc rotary kiln, the digital twin model built from the analysis of the multi-field coupling mechanism is only consistent with the physical characteristics of the actual rotary kiln in the model structure, and the key parameters in the model are required to keep synchronous with the continuously changing working conditions of the actual equipment. The digital twin model unknown parameters are obtained through analysis and processing based on limited observation data, and the regulation parameters and the disturbance data are used as input parameters of twin model simulation. The on-site available observation data are fewer, main observation data (namely preset observation data) comprise worker experience, kiln head flame/material morphology, kiln body/kiln tail temperature and kiln slag components, the digital twin unknown parameters of the zinc rotary kiln are required to be effectively analyzed and calculated based on the limited observation data, and dynamic tracking of a model to a real rotary kiln is realized, so that dynamic consistency of the twin model and actual equipment is ensured.
Specifically, in one implementation of the present embodiment, step S200 includes the steps of:
step S201, acquiring and storing temperature, concentration and pressure data in the running process of the rotary kiln in a discrete time sequence data form through a DCS system;
Step S202, analyzing and guiding rotary kiln working conditions according to on-site observation and operation mode experience of regulating and controlling a rotary kiln, and analyzing and processing DCS data and flame image observation data;
step S203, extracting the complement data corresponding to the rotary kiln state through time sequence data analysis and image processing.
In the embodiment, data such as temperature, concentration and pressure in the running process of the rotary kiln are collected and stored through the DCS system in the form of discrete time sequence data. As shown in fig. 3, the operation mode experience of the rotary kiln is regulated and controlled by combining the site observation of workers, the analysis and the guidance of the rotary kiln working condition are carried out, and the finite observation data such as DCS data, flame images and the like are analyzed and processed. And extracting the complement data related to the rotary kiln state through time sequence data analysis and image processing.
Specifically, in one implementation of the present embodiment, step S203 includes the steps of:
step S203a, extracting a time domain feature, a frequency domain feature, a static feature of the flame image and a dynamic feature of the flame image corresponding to the rotary kiln state through the time sequence data analysis and the image processing.
In this embodiment, the extracted complement data related to the rotary kiln state mainly includes a time domain feature, a frequency domain feature, a static feature of a flame image, and a dynamic feature of the flame image; in extracting data, it is necessary to pre-process the data.
Wherein the data preprocessing includes numerical transformation and supplementation of missing values, trend term elimination, smoothing and noise reduction, and encoding and transformation of data attributes. The trend term elimination is used for eliminating zero drift of the sensor signal caused by temperature change, and is realized by a polynomial least square method based curve fitting method. Let the rotary kiln sensor collect one item of data as Y= (Y) 1 ,y 2 ,...,y m ) Meanwhile, the fitting curve to the sampled data is assumed to be:
Figure SMS_6
to fit the sampled data as closely as possible, the sum of squares of the errors needs to be satisfied
Figure SMS_7
The least, i.e. the solution satisfying its least squares method. Find a= [ a ] n ,a n-1 ,...a 0 ] T And obtaining a fitting equation of the curve. A satisfies the following conditions:
Figure SMS_8
and calculating the above equation to obtain a fitting equation of the time series data. Smoothing and noise reduction of the sensing data are realized by adopting an average method, and the formula is as follows:
Figure SMS_9
feature extraction of the zinc rotary kiln sensing data is mainly carried out from a time domain and a frequency domain. Analysis is carried out in the time domain according to the time domain amplitude waveform, and the extracted characteristics are an effective value Xrms and a peak value X max The calculation formulas of the skewness alpha and the kurtosis beta are as follows:
Figure SMS_10
the extraction of the frequency domain features adopts Fourier transformation to convert the time domain signals into frequency domain signals for analysis. Let the discrete signal collected by the rotary kiln be (x) 0 ,x 1 ,x 2 ,...,x N-1 ) The fourier transform of its spectrum calculation is:
Figure SMS_11
to reduce the computational effort, a series of iterative operations, i.e., a fast fourier transform algorithm, is employed to reduce the computational effort. The effective characteristics of the rotary kiln DCS data are used for on-line identification and calibration of digital twin model parameters.
Specifically, in one implementation of the present embodiment, step S203 includes the following steps:
step S204, a parameter estimation value is calculated once in each observation of data recursion, and the digital twin model is updated and verified according to the prediction result.
In this embodiment, for nonlinear characteristics inside the rotary kiln and for implementing virtual-real synchronization of digital twin, a support vector regression method is adopted to implement model parameter online prediction, that is, parameter estimation values are recursively calculated once every time data is observed, and update and verification of a digital twin model are implemented according to a prediction result. If the verification result meets the requirement, completing the digital twin model of virtual-real synchronous mapping; and if the verification result is not satisfied, the identification of the model parameters is carried out again, so that the digital twin model of the zinc rotary kiln is recursively realized.
In one implementation manner of the embodiment, the analysis of the working condition and the reaction mechanism of the rotary kiln is used for obtaining the mass ratio of the temperature in the kiln and the related substances in the kiln, which has a great influence on the process index of the rotary kiln. Therefore, data acquisition is carried out by installing sensing equipment on a zinc smelting site, and the distribution results of the reaction rates and the mass ratio of the existing kilns with other sizes are obtained based on the existing literature, so that SVR model construction and prediction are carried out on different reaction rates in the rotary kiln. Based on the field data and the related data, the related results of the content of substances in the kiln and the reaction rate of 85 groups can be obtained, and SVR model training is carried out after data analysis processing and normalization, and is particularly shown in fig. 4.
Input x= [ m1 ] consisting of mass and temperature of components in kiln i ,m2 i ,...,mk i ,Ts i ,Tg i ,Tsh i ]To train different chemical reaction rates y= [ r1, r2, ], r9 within the kiln]. Where i=1, 2..n is the number of axial segments of the kiln, mk i For the mass of the ith segment of the mass k, ts (Ts) i Is the material temperature, tg i Tsh is the gas temperature i Is kiln wall temperature. And then training to obtain an SVR model for realizing reaction rate prediction. The input data for reaction rate prediction comes from a real volatilizing kilnCollecting data on site, including equipment parameters, material quality at initial time, gas and solid temperature, etc. to realize recursion prediction of reaction rate in kiln, therefore, key parameters of the rotary kiln digital twin model are continuously updated, and dynamic tracking of real equipment is realized. The pseudo code of the rotary kiln digital twin dynamic update based on the SVR prediction model is as follows:
Figure SMS_12
according to the embodiment, the maximum likelihood parameter estimation of the online recursion of the digital twin model of the zinc rotary kiln can be obtained by solving and taking the extremum, so that the online recursion method is used for updating the twin model coefficient, the coupling coefficient and the boundary condition, realizing the dynamic tracking of the zinc rotary kiln on site and ensuring the virtual-real synchronization of the digital twin.
As shown in fig. 1, in an implementation manner of the embodiment of the invention, the method for soft measurement of key monitoring variables of the zinc rotary kiln further comprises the following steps:
And step S300, aiming at the problem of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln, adopting a variable step length bidirectional finite difference to carry out digital twin simulation solution.
In the embodiment, on the basis of the digital twin model simulation, a digital twin simulation solving method based on a variable-step-length bidirectional Euler method is adopted for solving the problems of gas-solid countercurrent and gas-solid boundary conditions in the kiln, so that the simulation solving precision of key variables in the kiln is improved, the error in the material simulation solving process in the zinc rotary kiln is continuously reduced, and the method meets the test and actual conditions. Compared with the data actually collected on the industrial site, the result shows that the simulation solving effect of the zinc rotary kiln digital twin gas-solid two-phase countercurrent process is better.
Specifically, in one implementation of the present embodiment, step S300 includes the steps of:
step S301, dividing the rotary kiln into N units along the axial direction, and determining a mass conservation formula existing in each unit;
step S302, deducing in a differential mode according to a mass conservation formula in each unit to form a recursion formula;
and step S303, carrying out digital twin simulation solving by adopting the variable step length bidirectional finite difference.
At present, the simulation of key variables in a rotary kiln is mainly calculated by adopting a finite difference method, however, due to the unique characteristics of bidirectional countercurrent of materials and gases of the rotary kiln, the problem that the boundary conditions of the materials and gases are precisely defined is difficult to solve by the existing method, and the solving process is divergent. The modeling process needs to consider the problem that the initial conditions are difficult to solve due to the fact that the gas-solid countercurrent of the zinc rotary kiln is not located at the same boundary, so that the mechanism modeling and simulation of a material component field in the convection reaction process in the kiln are realized, the coupling solution of a temperature field and a chemical reaction field is realized, and the accuracy of the simulation result in the zinc rotary kiln is improved.
In this embodiment, as shown in fig. 5, in order to construct and simulate a reaction field model in the rotary kiln, the rotary kiln is divided into N units along the axial direction, and a mass conservation formula exists in each unit:
Figure SMS_13
wherein the method comprises the steps of
Figure SMS_14
For the mass of a substance j at node i, < +.>
Figure SMS_15
For the rate of reaction of substance j at node i, M j The amount of the substance j, v is the movement speed of the material or the gas, wherein the axial dimension of the zinc rotary kiln is 68m, and the inner diameter is 4.5m. Deriving by adopting a differential form, thereby forming a recursive formula:
Figure SMS_16
the meaning of this formula is that the mass of the substance j at node i+1 is equal to the sum of the mass of the substance j at the previous node i and the mass of the reaction product substance. The material enters the zinc rotary kiln from the kiln tail, and the initial condition of the material exists at the kiln head; the gas is blown into the rotary kiln from the kiln tail, and the initial condition of the gas exists at the kiln tail. If the rotary kiln digital twin model established by the analysis of the chemical reaction and the mass conservation mechanism is solved by adopting a finite difference method, the simulation result error is large due to the fact that initial conditions are defined on the same boundary, and the solution is not in accordance with the actual situation due to the fact that simulation divergence of certain reaction processes is caused. Meanwhile, because the temperature change speed of different sections in the kiln is different, the finite difference with a fixed step length is adopted, so that larger errors exist in transition sections of certain nodes such as a preheating zone and a high-temperature reaction zone. Therefore, a variable-step bidirectional finite difference method is provided for carrying out digital twin solution simulation of the zinc rotary kiln, and the simulation is particularly shown in fig. 5.
Specifically, in one implementation of the present embodiment, step S303 includes the steps of:
step S303a, aiming at the problem that the material and gas phase countercurrent and boundary conditions in the zinc rotary kiln are not in the same position, defining a kiln head material composition matrix MH1 at the kiln head of the rotary kiln, and defining a kiln tail material composition matrix MT1 at the kiln tail of the rotary kiln;
step S303b, obtaining initial conditions of a gas phase in a kiln head material composition matrix MH1 according to compressed air and oxygen-enriched drum intake calculation of the kiln head, according to the raw material composition of kiln tail feeding, obtaining initial conditions of material phases in a kiln tail material composition matrix MT1;
step S303c, respectively solving a kiln head material component matrix MH1 and a kiln tail material component matrix MT1, updating the material phase in the kiln head material component matrix MHi by the material phase result in the kiln tail material component matrix MTi after one-time solving, and updating the gas phase in the kiln tail material component matrix MTi by the gas phase result in the kiln head material component matrix MHi.
In this embodiment, grid cells are divided along the length direction of the kiln, and the step dx of each cell is defined as follows based on the change slope of the temperature in the cell of the preliminary analysis, so that densification of cell division is realized in a region with a severe temperature change range, solving errors caused by rapid temperature change are eliminated, and twin simulation precision is improved.
Figure SMS_17
Aiming at the problem that the material and gas countercurrent in the zinc rotary kiln and the boundary conditions are not in the same position, 1 material component matrix MH1 (kiln head) and MT1 (kiln head) with the same size are respectively defined at the kiln head and the kiln tail of the rotary kiln, the number of rows in the matrix represents the number N of nodes, and the number of columns represents the material type.
ZnFe mainly taking solid state items in kiln head area reaction into consideration in zinc rotary kiln 2 O 4 、ZnSO 4 、Fe 3 O 4 SO of FeO, fe, znO, C and gaseous terms 2 、CO、CO 2 、O 2 Composition change, and C, fe, feO, znO (solid) of high temperature reaction zone material phase and CO, CO of gas phase 2 、O 2 Reaction process of 9 substances of Zn (gas state) and ZnO (gas state). The initial conditions of the material phases in the kiln tail component matrix MT1 are obtained by assaying the raw material components fed from the kiln tail, the initial conditions of the gas phase in the kiln head component matrix MH1 are calculated by compressed air and oxygen-enriched air blowing amount of the kiln head. In the digital twin model solving process of the rotary kiln, each iteration process carries out the respective solving of two matrixes, after the one-time solving is completed, the material phase in the kiln head component matrix MHi is updated by the material phase result in the kiln head component matrix MTi, and the gas phase in the kiln head component matrix MHi is updated by the gas phase result in the kiln head component matrix MHi.
After the primary solution is finished, calculating the algorithm difference of the two matrixes MTi and MHi, and if the error does not meet the requirement that the simulation calculation does not reach the kiln stable state, continuing to carry out iterative solution; until the two matrix MHK and MTK errors meet the requirements, the solving process is finished, at the moment, the two component matrices in the kiln are considered to be approximate, and the simulation result also reaches the stable state in the kiln; at the moment, the content of each material component in the kiln obtained by calculation and solution is considered as the working condition state of the rotary kiln after the rotary kiln is regulated and controlled and stably operated. When the new working conditions such as the rotary kiln rotating speed, the blanking speed, the blast volume and the like change, simulation solution under the new working conditions is carried out to obtain key monitoring variables in the rotary kiln, so that the method is used for guiding workers to carry out optimization regulation and control on the rotary kiln.
As shown in fig. 1, in an implementation manner of the embodiment of the invention, the method for soft measurement of key monitoring variables of the zinc rotary kiln further comprises the following steps:
step S400, dynamically updating digital twin parameters according to the changes of equipment regulation and control parameters and external disturbance data, and timely simulating to obtain the temperature and chemical components in the rotary kiln, so that real-time feedback of regulation and control operation is realized.
In this embodiment, on the basis of the above-mentioned digital twin simulation solving method, dynamic update is performed according to the changes of the on-site equipment regulation parameters and external disturbance data, so as to guide a worker to perform optimized regulation and control of the rotary kiln.
Specifically, in one implementation of the present embodiment, step S400 includes the following steps:
and S401, obtaining variables of the temperature, the reaction rate, the zinc generation progress and the carbon consumption progress in the rotary kiln through dynamically updated digital twin simulation according to the equipment regulation and control parameters and the change of the external disturbance data, and realizing real-time feedback of regulation and control operation.
In the embodiment, the equipment regulation parameters comprise parameters such as blanking speed, kiln body rotating speed, blasting state and the like, and the external disturbance data comprise ambient temperature, coke grade and coke duty ratio; according to the change of the parameters, key variables such as the temperature, the reaction rate, the zinc generation progress, the carbon consumption progress and the like in the rotary kiln can be obtained through dynamically updated digital twin in-time simulation, and the real-time feedback of the regulation and control operation is realized for the regulation and control reference of workers.
In an actual application scenario of the embodiment, a digital twin simulation kernel in the kiln is constructed through analysis of the mechanical coupling of a temperature field and a chemical reaction field of the zinc rotary kiln, and based on the analysis of working condition data and the updating of twin parameters, a digital twin software platform for monitoring the key state of the zinc rotary kiln is realized as shown in fig. 6. The software platform adopts an MVVM architecture (Model-View Model) to realize the separation of rich visual information and the front end and the rear end of the twin mechanism kernel. The platform function mainly comprises functions of rotary kiln running state monitoring, working condition identification, data analysis, parameter updating, key variable simulation and the like.
Wherein the temperature field and the chemical reaction process of the rotary kiln digital twin kernel part adopt the bidirectional finite difference with variable step length to carry out solution simulation on the basis of mechanism analysis modeling and parameter dynamic updating, the improvement of the simulation result precision of the digital twin model of the zinc rotary kiln is realized, and the key monitoring variables reflecting the working conditions in the kiln are finally obtained through simulation, wherein the key monitoring variables comprise an axial temperature field model and an axial component reaction progress model, and the result is shown in figure 7.
The rotary kiln digital twin simulation solution is carried out by adopting a traditional method, the gas phase calculation process diverges, and the solution result deviates from the actual condition far, for example, the mass ratio of the kiln head is more than 1, the CO2 component is less than 0 and the like obtained by calculating the gas phase O2 component. Adopting a variable-step bidirectional finite difference solution digital twin model of the zinc rotary kiln, wherein the error in the material simulation solution process in the zinc rotary kiln is continuously reduced, and the method accords with the test and actual conditions; the digital twin simulation result comparison verification is carried out through the CO concentration component of the kiln tail and the Zn content in the kiln slag at present under the limitation of the current data detection conditions, and the kiln slag Zn content and the kiln tail CO content in the simulation result of the zinc rotary kiln are more consistent with the actual DCS monitoring result through multiple comparison verification. The simulation method provided by the embodiment improves the precision of digital twin simulation solving of the rotary kiln and provides feedback guidance for rotary kiln parameter regulation and control.
In another implementation manner of this embodiment, the corresponding modification scheme includes one or a combination of the following:
(1) The embodiment provides a coupling modeling thought of chemical reaction of a temperature field, is not limited to coupling of the two reaction processes, and can be popularized to coupling modeling of the temperature field and a flow field, and coupling modeling of the temperature field and a material field.
(2) The embodiment provides a variable-step bidirectional finite difference digital twin simulation solving method of the zinc rotary kiln, and can be popularized to further methods belonging to finite differences, such as Euler-Euler method, bidirectional Dragon-Gregorian tower method and bidirectional improved Euler method.
(3) The embodiment provides a digital twin virtual-real synchronization method and concept of a zinc rotary kiln, which mainly comprises operation mode analysis, time sequence data analysis flame characteristic extraction and twin parameter generation, and the actual implementation can comprise all or a plurality of methods.
(4) The embodiment provides a method and a concept of digital twin virtual-real synchronization of a zinc rotary kiln, wherein data preprocessing comprises numerical conversion and supplement of missing values, trend term elimination, smoothing and noise reduction and coding and conversion of data attributes. Implementations may include all or some of the methods.
(5) Aiming at the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions of the zinc rotary kiln, the embodiment provides a digital twin simulation solving method based on a variable step length bidirectional finite difference, and the simulation solving precision of key variables in the kiln is improved. The method is also applicable to other rotary kilns such as cement rotary kilns, garbage incineration rotary kilns and lime rotary kilns which have gas-solid two-phase countercurrent scenes.
(6) The embodiment provides a soft measurement method of key variables in a zinc rotary kiln based on digital twinning, which is not limited to soft measurement of temperature and chemical reaction components, and can be also applied to soft measurement of variables such as material movement, gas flow and the like in the rotary kiln.
The following technical effects are achieved through the technical scheme:
according to the embodiment, by constructing a digital twin model for reflecting the reaction condition of the zinc rotary kiln, synchronous mapping is kept with the on-site rotary kiln in real time, and the defect that the mechanism data is difficult to fuse, virtual and real are difficult to synchronize in the traditional modeling simulation is overcome; the monitoring results of key variables in the kiln, including the temperature of a high-temperature reaction zone, the reaction rate, the zinc recovery progress and the carbon consumption progress, are obtained through digital twin real-time simulation, so that the problem that the reaction condition of the high-temperature zone can only be inferred by indirectly relying on limited observation data is solved, the regulation and control effect of the rotary kiln is fed back in time, and the green intelligent operation of the zinc rotary kiln is serviced; the embodiment provides a zinc rotary kiln digital twin model construction and virtual-real interaction method, and key monitoring variable soft measurement of the zinc rotary kiln is realized through digital twin.
Exemplary apparatus
Based on the above embodiment, the invention also provides a key monitoring variable soft measurement device of a zinc rotary kiln, comprising:
the model construction module is used for analyzing the conservation theory of the temperature field and the chemical reaction process in the zinc rotary kiln and constructing a mechanism model of coupling of the temperature and the chemical reaction of the zinc rotary kiln;
the digital twin parameter module is used for analyzing and mining through preset observation data, identifying key parameters of the model and obtaining digital twin parameters of the rotary kiln consistent with the working conditions of equipment;
the simulation solving module is used for carrying out digital twin simulation solving by adopting a variable step length bidirectional finite difference aiming at the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln;
and the dynamic updating module is used for dynamically updating the digital twin parameters according to the change of the equipment regulation and control parameters and the external disturbance data, timely simulating to obtain the temperature and chemical components in the rotary kiln, and realizing the real-time feedback of the regulation and control operation.
Based on the above embodiment, the present invention also provides a terminal, and a functional block diagram thereof may be shown in fig. 8.
The terminal comprises: the system comprises a processor, a memory, an interface, a display screen and a communication module which are connected through a system bus; wherein the processor of the terminal is configured to provide computing and control capabilities; the memory of the terminal comprises a storage medium and an internal memory; the storage medium stores an operating system and a computer program; the internal memory provides an environment for the operation of the operating system and computer programs in the storage medium; the interface is used for connecting external equipment such as mobile terminals, computers and other equipment; the display screen is used for displaying corresponding information; the communication module is used for communicating with a cloud server or a mobile terminal.
The computer program is used for realizing the operation of the zinc rotary kiln key monitoring variable soft measurement method when being executed by a processor.
It will be appreciated by those skilled in the art that the functional block diagram shown in fig. 8 is merely a block diagram of some of the structures associated with the present inventive arrangements and is not limiting of the terminal to which the present inventive arrangements may be applied, and that a particular terminal may include more or less components than those shown, or may combine some of the components, or have a different arrangement of components.
In one embodiment, a terminal is provided, including: the system comprises a processor and a memory, wherein the memory stores a zinc rotary kiln key monitoring variable soft measurement program which is used for realizing the operation of the zinc rotary kiln key monitoring variable soft measurement method when being executed by the processor.
In one embodiment, a storage medium is provided, wherein the storage medium stores a zinc rotary kiln key monitored variable soft measurement program, which when executed by a processor is used to implement the operations of the zinc rotary kiln key monitored variable soft measurement method as above.
Those skilled in the art will appreciate that implementing all or part of the above-described methods may be accomplished by way of a computer program comprising instructions for the relevant hardware, the computer program being stored on a non-volatile storage medium, the computer program when executed comprising the steps of the embodiments of the methods described above. Any reference to memory, database, or other medium used in embodiments provided herein may include non-volatile and/or volatile memory.
In summary, the invention provides a method, a device, a terminal and a medium for soft measurement of key monitoring variables of a zinc rotary kiln, wherein the method comprises the following steps: analyzing a temperature field conservation theory and a chemical reaction process in the zinc rotary kiln, and constructing a mechanism model of coupling of the temperature and the chemical reaction of the zinc rotary kiln; analyzing and mining through preset observation data, and identifying key parameters of the model to obtain digital twin parameters of the rotary kiln consistent with equipment working conditions; aiming at the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln, adopting a variable step length bidirectional finite difference to carry out digital twin simulation solution; according to the change of the equipment regulation and control parameters and the external disturbance data, dynamically updating the digital twin parameters and timely simulating to obtain the temperature in the rotary kiln and the chemical composition in the kiln, so that the real-time feedback of the regulation and control operation is realized. The invention provides a zinc rotary kiln digital twin model construction and virtual-real interaction method, and key monitoring variable soft measurement of a zinc rotary kiln is realized through digital twin.
It is to be understood that the invention is not limited in its application to the examples described above, but is capable of modification and variation in light of the above teachings by those skilled in the art, and that all such modifications and variations are intended to be included within the scope of the appended claims.

Claims (13)

1. The key monitoring variable soft measurement method for the zinc rotary kiln is characterized by comprising the following steps of:
analyzing a temperature field conservation theory and a chemical reaction process in the zinc rotary kiln, and constructing a mechanism model of coupling of the temperature and the chemical reaction of the zinc rotary kiln;
analyzing and mining through preset observation data, and identifying key parameters of the model to obtain digital twin parameters of the rotary kiln consistent with equipment working conditions;
aiming at the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln, adopting a variable step length bidirectional finite difference to carry out digital twin simulation solution;
according to the change of the equipment regulation and control parameters and the external disturbance data, dynamically updating the digital twin parameters and timely simulating to obtain the temperature in the rotary kiln and the chemical composition in the kiln, so that the real-time feedback of the regulation and control operation is realized.
2. The method for soft measurement of key monitoring variables of a zinc rotary kiln according to claim 1, wherein the analysis of temperature field conservation theory and chemical reaction process in the zinc rotary kiln comprises the following steps:
the rotary kiln is divided into a kiln head area, a high-temperature reaction area and a kiln tail area, wherein a complex heat transfer process between materials, gas and kiln walls occurs in the sectional area, and a chemical reaction process between the materials and the gas occurs in the high-temperature reaction area in sections.
3. The method for soft measurement of key monitoring variables of a zinc rotary kiln according to claim 1, wherein the analysis of the temperature field conservation theory and the chemical reaction process in the zinc rotary kiln to construct a mechanism model of coupling of the temperature and the chemical reaction of the zinc rotary kiln comprises the following steps:
analyzing convection heat transfer and radiation heat transfer among the flue gas, the kiln wall and the materials, and constructing a heat conservation equation of the flue gas, the materials and the kiln wall;
the high-temperature decomposition of the preheating section and the coke combustion of the high-temperature reaction zone, the reduction of zinc oxide and the reaction of iron compounds mainly occur in the volatilizing kiln;
and (3) based on the principle analysis of chemical reaction by mass conservation and energy conservation, constructing a principle model of the coupling of the temperature of the zinc rotary kiln and the chemical reaction.
4. The method for soft measurement of key monitored variables of a zinc rotary kiln according to claim 3, wherein the constructing a mechanism model of the zinc rotary kiln temperature and the chemical reaction coupling comprises:
temperature field analysis and modeling are independently carried out, and a preliminary result of the temperature field is used as a main calculation basis of the reaction rate in the chemical reaction;
and correcting the temperature field and modeling the coupling through analysis and calculation of enthalpy change in the chemical reaction process.
5. The method for soft measurement of key monitoring variables of a zinc rotary kiln according to claim 1, wherein the analyzing and mining by preset observation data, identifying key parameters of a model, and obtaining digital twin parameters of the rotary kiln consistent with equipment working conditions comprises the following steps:
acquiring and storing temperature, concentration and pressure data in the running process of the rotary kiln in a discrete time sequence data form through a DCS system;
according to the operation mode experience of on-site observation and regulation of the rotary kiln, analyzing and guiding rotary kiln working conditions, and analyzing and processing DCS data and flame image observation data;
and extracting the complement data corresponding to the rotary kiln state through time sequence data analysis and image processing.
6. The method for soft measurement of key monitoring variables of a zinc rotary kiln according to claim 5, wherein the extracting of the complement data corresponding to the rotary kiln state through time series data analysis and image processing comprises the following steps:
and extracting time domain features, frequency domain features, static features of the flame image and dynamic features of the flame image corresponding to the rotary kiln state through the time sequence data analysis and the image processing.
7. The method for soft measurement of key monitoring variables of a zinc rotary kiln according to claim 5, wherein the extracting of the complement data corresponding to the rotary kiln state through time series data analysis and image processing comprises the following steps:
And (3) recursively calculating a parameter estimation value once per observation of data, and updating and verifying the digital twin model according to the prediction result.
8. The method for soft measurement of key monitoring variables of a zinc rotary kiln according to claim 1, wherein aiming at the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln, digital twin simulation solution is performed by adopting variable step-length bidirectional finite difference, and the method comprises the following steps:
dividing the rotary kiln into N units along the axial direction, and determining a mass conservation formula existing in each unit;
deducing in a differential mode according to a mass conservation formula in each unit to form a recursion formula;
and carrying out digital twin simulation solving by adopting the variable step length bidirectional finite difference.
9. The method for soft measurement of key monitoring variables of a zinc rotary kiln according to claim 8, wherein the step-variable bidirectional finite difference is adopted for digital twin simulation solution, and the method comprises the following steps:
aiming at the problem that the material and gas countercurrent and boundary conditions in the zinc rotary kiln are not in the same position, a kiln head material composition matrix MH1 is defined at the kiln head of the rotary kiln, and a kiln tail material composition matrix MT1 is defined at the kiln tail of the rotary kiln;
Obtaining initial conditions of a gas phase in a kiln head material component matrix MH1 according to compressed air and oxygen-enriched drum charging amount of the kiln head, and obtaining initial conditions of a material phase in a kiln tail material component matrix MT1 according to raw material component assay of kiln tail feeding;
and respectively solving a kiln head material component matrix MH1 and a kiln tail material component matrix MT1, updating the material phase in the kiln head material component matrix MHi by the material phase result in the kiln tail material component matrix MTi after one-time solving, and updating the gas phase in the kiln tail material component matrix MTi by the gas phase result in the kiln head material component matrix MHi.
10. The method for soft measurement of key monitoring variables of a zinc rotary kiln according to claim 1, wherein dynamically updating digital twin parameters and timely simulating to obtain the temperature and chemical composition in the rotary kiln according to the change of equipment regulation parameters and external disturbance data comprises the following steps:
and according to the equipment regulation and control parameters and the change of the external disturbance data, obtaining variables of the temperature, the reaction rate, the zinc generation progress and the carbon consumption progress in the rotary kiln through dynamically updated digital twin simulation, and realizing real-time feedback of regulation and control operation.
11. The utility model provides a zinc rotary kiln key monitoring variable soft measurement device which characterized in that includes:
the model construction module is used for analyzing the conservation theory of the temperature field and the chemical reaction process in the zinc rotary kiln and constructing a mechanism model of coupling of the temperature and the chemical reaction of the zinc rotary kiln;
the digital twin parameter module is used for analyzing and mining through preset observation data, identifying key parameters of the model and obtaining digital twin parameters of the rotary kiln consistent with the working conditions of equipment;
the simulation solving module is used for carrying out digital twin simulation solving by adopting a variable step length bidirectional finite difference aiming at the problems of gas-solid two-phase countercurrent and gas-solid boundary conditions in the kiln;
and the dynamic updating module is used for dynamically updating the digital twin parameters according to the change of the equipment regulation and control parameters and the external disturbance data, timely simulating to obtain the temperature and chemical components in the rotary kiln, and realizing the real-time feedback of the regulation and control operation.
12. A terminal, comprising: the system comprises a processor and a memory, wherein the memory stores a zinc rotary kiln key monitoring variable soft measurement program which is used for realizing the operation of the zinc rotary kiln key monitoring variable soft measurement method according to any one of claims 1-10 when being executed by the processor.
13. A medium, characterized in that the medium is a computer readable storage medium, the medium stores a zinc rotary kiln key monitoring variable soft measurement program, and the zinc rotary kiln key monitoring variable soft measurement program is used for implementing the operation of the zinc rotary kiln key monitoring variable soft measurement method according to any one of claims 1-10 when being executed by a processor.
CN202310150189.1A 2023-02-06 2023-02-06 Zinc rotary kiln key monitoring variable soft measurement method, device, terminal and medium Pending CN116246724A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310150189.1A CN116246724A (en) 2023-02-06 2023-02-06 Zinc rotary kiln key monitoring variable soft measurement method, device, terminal and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310150189.1A CN116246724A (en) 2023-02-06 2023-02-06 Zinc rotary kiln key monitoring variable soft measurement method, device, terminal and medium

Publications (1)

Publication Number Publication Date
CN116246724A true CN116246724A (en) 2023-06-09

Family

ID=86634625

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310150189.1A Pending CN116246724A (en) 2023-02-06 2023-02-06 Zinc rotary kiln key monitoring variable soft measurement method, device, terminal and medium

Country Status (1)

Country Link
CN (1) CN116246724A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116520712A (en) * 2023-07-03 2023-08-01 贵阳职业技术学院 Self-adaptive regulation and control method, system, terminal and medium for feeding amount of reaction container

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116520712A (en) * 2023-07-03 2023-08-01 贵阳职业技术学院 Self-adaptive regulation and control method, system, terminal and medium for feeding amount of reaction container
CN116520712B (en) * 2023-07-03 2023-08-22 贵阳职业技术学院 Self-adaptive regulation and control method, system, terminal and medium for feeding amount of reaction container

Similar Documents

Publication Publication Date Title
CN108764517B (en) Method, equipment and storage medium for predicting change trend of silicon content in molten iron of blast furnace
CN112100916B (en) Method, device, electronic equipment and medium for constructing reinforcement learning model
CN107299170B (en) A kind of blast-melted quality robust flexible measurement method
CN111831719A (en) Intelligent control method and system for blast furnace ironmaking production process
CN107038307B (en) The Roller Conveying Kiln for Temperature that mechanism is combined with data predicts integrated modelling approach
CN107368125B (en) A kind of blast furnace temperature control system and method based on CBR Yu the parallel mixed inference of RBR
CN111915080B (en) Raw fuel cost optimal proportioning method based on molten iron quality constraint
CN116246724A (en) Zinc rotary kiln key monitoring variable soft measurement method, device, terminal and medium
CN106709197A (en) Molten iron silicon content predicting method based on slide window T-S fuzzy neural network model
CN105974896B (en) A kind of pneumatic steelmaking Optimal Control System and method based on information physical fusion
CN110598958B (en) Ladle hierarchical management analysis method and system
CN105574297B (en) Self adaptation blast furnace molten iron silicon content trend prediction method
Jiang et al. Real-time moisture control in sintering process using offline–online NARX neural networks
CN113192568B (en) Method and system for forecasting desulfurization end point of refining furnace
CN114968997B (en) Sintering end point forecasting method based on space-time feature extraction
CN111666683A (en) Method and system for predicting distribution state of furnace burden in blast furnace
CN116738863B (en) Method for building digital twin-based external refining CO2 digital management platform
CN104865367A (en) Dry quenching coke burn-out rate real-time measurement method
CN118039015A (en) Converter endpoint element content prediction method based on dynamic time sequence characteristics
CN113512620B (en) Dynamic control method for endpoint carbon in whole converter smelting process of gas analysis and sublance
Cui et al. [Retracted] Prediction Model of Hot Metal Silicon Content Based on Improved GA‐BPNN
CN114525372B (en) Multi-mode fusion based blast furnace state monitoring method and device
Mitra et al. Evolutionary algorithms in ironmaking applications
Niu et al. Prediction of sinter chemical composition based on ensemble learning algorithms
An et al. Applications of evolutionary computation and artificial intelligence in metallurgical industry

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