[go: up one dir, main page]

WO2004050920A1 - Procede de regulation du fonctionnement d'une cuve a reaction, dispositif, programme, et support d'enregistrement lisible par ordinateur - Google Patents

Procede de regulation du fonctionnement d'une cuve a reaction, dispositif, programme, et support d'enregistrement lisible par ordinateur Download PDF

Info

Publication number
WO2004050920A1
WO2004050920A1 PCT/JP2003/012997 JP0312997W WO2004050920A1 WO 2004050920 A1 WO2004050920 A1 WO 2004050920A1 JP 0312997 W JP0312997 W JP 0312997W WO 2004050920 A1 WO2004050920 A1 WO 2004050920A1
Authority
WO
WIPO (PCT)
Prior art keywords
temperature
reaction vessel
heat flux
inverse problem
problem analysis
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.)
Ceased
Application number
PCT/JP2003/012997
Other languages
English (en)
Japanese (ja)
Inventor
Hiroyuki Yoshino
Yasushi Kurita
Masatoshi Miyawaki
Hiroyuki Furuta
Junichi Nakagawa
Masahiro Yamamoto
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.)
Nippon Steel Corp
Original Assignee
Nippon Steel Corp
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
Priority claimed from JP2002348429A external-priority patent/JP2004003800A/ja
Priority claimed from JP2003044672A external-priority patent/JP4299554B2/ja
Priority claimed from JP2003085817A external-priority patent/JP4559708B2/ja
Application filed by Nippon Steel Corp filed Critical Nippon Steel Corp
Publication of WO2004050920A1 publication Critical patent/WO2004050920A1/fr
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • C—CHEMISTRY; METALLURGY
    • C21—METALLURGY OF IRON
    • C21B—MANUFACTURE OF IRON OR STEEL
    • C21B7/00—Blast furnaces
    • C21B7/24—Test rods or other checking devices
    • C—CHEMISTRY; METALLURGY
    • C21—METALLURGY OF IRON
    • C21C—PROCESSING OF PIG-IRON, e.g. REFINING, MANUFACTURE OF WROUGHT-IRON OR STEEL; TREATMENT IN MOLTEN STATE OF FERROUS ALLOYS
    • C21C5/00—Manufacture of carbon-steel, e.g. plain mild steel, medium carbon steel or cast steel or stainless steel
    • C21C5/28—Manufacture of steel in the converter
    • C21C5/42—Constructional features of converters
    • C21C5/46—Details or accessories
    • C21C5/4673—Measuring and sampling devices
    • C—CHEMISTRY; METALLURGY
    • C21—METALLURGY OF IRON
    • C21B—MANUFACTURE OF IRON OR STEEL
    • C21B2300/00—Process aspects
    • C21B2300/04—Modeling of the process, e.g. for control purposes; CII
    • C—CHEMISTRY; METALLURGY
    • C21—METALLURGY OF IRON
    • C21C—PROCESSING OF PIG-IRON, e.g. REFINING, MANUFACTURE OF WROUGHT-IRON OR STEEL; TREATMENT IN MOLTEN STATE OF FERROUS ALLOYS
    • C21C2300/00—Process aspects
    • C21C2300/06—Modeling of the process, e.g. for control purposes; CII
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P10/00—Technologies related to metal processing
    • Y02P10/25—Process efficiency

Definitions

  • the present invention relates to a reaction vessel operation management method suitable for use in managing the operation of a reaction vessel involving a high-temperature gas reaction or a liquid reaction such as a blast furnace, a steel material heating furnace by combustion, a coal gasification reaction furnace, etc.
  • the present invention relates to an apparatus, a computer program, and a computer-readable recording medium. Background art
  • reaction vessel involving a high-temperature gas reaction or a liquid reaction such as a blast furnace, a steel heating furnace by combustion, or a coal gasification reaction furnace
  • a reaction vessel for example, combustion behavior
  • the state inside the reaction vessel has been estimated based on the temperature measured by a thermocouple embedded inside the wall of the reaction vessel. For example, if there is a sudden rise in temperature, it is estimated that abnormal heat is generated in the reaction vessel around the thermocouple, and if there is an extreme drop in temperature, heat is generated in the reaction vessel around the thermocouple. This is an empirical method such as estimating that a decrease in the calorific value has occurred such as a reduction in the reaction area.
  • Figure 35 shows a two-dimensional cross section near the furnace wall of a reaction vessel (heating furnace) in which multiple thermocouples “X” are embedded.
  • the boundary is indicated by a broken line in the furnace wall, one-dimensional means that only the heat flow in the direction along the broken line is considered. That is, for example, assuming heat conduction in la ⁇ lb ⁇ lc or ld ⁇ le directions, the heat flux on the furnace inner surface is estimated.
  • it is general to determine the heat flux on the unknown inner furnace surface, assuming that the cooling condition of the outer furnace surface is known. Of course, it is also possible to reverse the known and unknown boundary conditions.
  • Japanese Patent Application Laid-Open Publication No. 2000-1 — 324 217 discloses an inverse problem analysis of an unsteady one-dimensional heat conduction equation from a thermocouple embedded in a blast furnace hearth. It describes a method for estimating the heat flux.
  • One of these methods is to estimate the heat flux at the opposite end point from the thermocouple temperature change at one point and the cooling condition (assumed to be known) at the end point.
  • Such cooling conditions are given by the heat transfer coefficient and the temperature of the cooling water.
  • uncertain estimated values This may have a negative effect on the accuracy of the heat flux estimate at the opposite end point, estimated using the value in the inverse problem.
  • thermocouple temperature change is described.However, since it is a method to solve by giving one of the two points as a fixed temperature boundary condition However, it is difficult to estimate by capturing the relative temperature change between the two points, and it is possible to estimate the heat flux under the fixed temperature boundary condition, but outside of the side selected as the fixed temperature boundary condition The end point heat flux on the extension line cannot be estimated.
  • the change in thickness due to the in-furnace melt attached to the surface of the refractory and the change in heat flux are simultaneously estimated, rather than the method of finding the heat flux at both ends while fixing the analysis length.
  • Method. Introducing a mouthpiece that increases or decreases the amount of adherence by coagulation and dissolution phenomena into the inverse problem analysis firstly complicates the calculation procedure and makes the calculation unstable. Second, if a calculation procedure that changes the analysis length is entered at each time step, uncertainties may be mixed in the method of estimating the temperature distribution before and after the change in the length. Estimation accuracy may be worse I can't deny it.
  • the conventional inverse problem analysis method is inadequate in many respects, and a new method for simultaneously estimating the heat flux at both ends by fixing the analysis length from multiple thermocouple information has been established. Therefore, a technique for accurately and stably estimating unsteady changes in heat flux is important.
  • the original unsteady one-dimensional inverse one-dimensional heat conduction problem is to estimate the boundary conditions on the inner and outer surfaces of the furnace at the same time.
  • temperature fluctuations measured by a certain thermocouple may be caused by a change in heat flux on the inner surface of the reaction vessel as described above, or may be caused by poor contact of a cooling device installed outside the reaction vessel. It cannot be distinguished whether the change is due to a change in heat flux on the outer surface of the reaction vessel.
  • the heat conduction phenomenon should also occur in the vertical direction across the broken line shown in Fig. 35, and it is necessary to solve the inverse heat conduction problem in two dimensions. In this case, even if the upper and lower boundaries in Fig. 35 are assumed to be adiabatic, it is necessary to construct a two-dimensional inverse problem for estimating the fine heat flux distribution on the left and right boundaries.
  • thermocouple information a method of estimating the temperature distribution by trial and error and simultaneously estimating the temperature distribution at both ends is also conceivable, so that the temperature change can be sufficiently expressed from multiple pieces of thermocouple information.
  • the calculation becomes complicated as the number of thermocouples increases, and it is extremely difficult to obtain a temperature distribution solution that satisfies the measured temperature changes of all thermocouples.
  • the unknown heat flux is calculated by alternately changing the two end points, and apparently the heat flux of the end point can be estimated at the same time.
  • the measurement temperature as the fixed temperature boundary condition is alternately calculated repeatedly, and when the measured temperature and the calculated temperature of both thermocouples agree to some extent, both ends in the time step Heat flux solution.
  • thermocouple positions when minimizing the absolute value of the difference between the measured temperature and the calculated temperature at two thermocouple positions, a criterion for how to balance the minimization at two independent thermocouple positions is required. It is difficult to set properly. It goes without saying that if this method is extended to the case of multiple thermocouples, it will be extremely difficult to determine the solution.
  • the time course of the heat flux obtained by inverse problem analysis has been independently plotted and evaluated.However, if the change of the heat flux is large or small, it is very vague, It is difficult to judge how non-stationary changes (rapid changes) are in the state inside the container.
  • the present invention has been made in view of the points mentioned above, and c to that of the heat flux changes and temperature changes in the inner and outer surfaces of the reaction vessel and can be simultaneously estimated the main purpose is the inverse problem The purpose is to stabilize the analysis and to evaluate the degree of unsteadiness of the state in the reaction vessel. Disclosure of the invention
  • the operation management method for a reaction container is a method for managing the operation of a reaction container for controlling the operation of a reaction container accompanied by a temperature change reaction, wherein a plurality of the reaction containers are arranged at least in a thickness direction inside a wall of the reaction container. Based on the temperature measured at the measured temperature measurement point, the inverse problem analysis using an unsteady heat conduction equation is performed, thereby obtaining a temperature or a heat flux on the inner surface and the outer surface of the reaction vessel. It has features.
  • another feature of the operation management method for the reaction vessel of the present invention is that the temperature measured at each of the temperature measurement points, the temperature on the inner surface and the outer surface of the reaction vessel, Or the minimum value of the sum of the squares of the difference between the temperature at each of the temperature measurement points calculated from the assumed value of the heat flux and the unsteady heat conduction equation is the inner surface and outer surface of the reaction vessel. In that the temperature or the heat flux is determined.
  • the method includes a step of increasing the number of digits after the decimal point of the temperature data measured at the temperature measurement point, and The point is that the increased number of temperature data is used for the inverse problem analysis.
  • the method includes a step of performing a filtering process on the temperature data measured at the temperature measurement point, and using the temperature data after the filtering process in the inverse problem analysis.
  • the method includes a step of increasing the number of digits after the decimal point of the temperature data measured at the temperature measurement point, and a step of filtering the temperature data with the number of digits after the decimal point increased. Is used in the above inverse problem analysis.
  • Another feature of the operation management method for a reaction vessel of the present invention is that, based on the temperature data measured at the temperature measurement points, analysis is performed by a steady-state method, whereby the inner surface of the reaction vessel is analyzed.
  • the analysis procedure by the steady-state method for obtaining the temperature or the heat flux, the temperature or the heat flux on the inner surface of the reaction vessel obtained by the inverse problem analysis, and the temperature or the heat flux obtained by the analysis procedure by the steady-state method And a comparison procedure for comparison.
  • the inverse problem analysis is an inverse problem analysis using an interpolation function or an external ⁇ function that satisfies an unsteady heat conduction equation. On the point.
  • An operation management apparatus for a reaction vessel is an operation management apparatus for a reaction vessel for managing an operation of a reaction vessel accompanied by a temperature change reaction, and a plurality of the operation management apparatuses are arranged at least in a thickness direction inside a wall of the reaction vessel.
  • the computer program of the present invention is a computer program for managing the operation of a reaction vessel accompanied by a temperature change reaction, wherein the computer program is measured at least at a plurality of temperature measurement points arranged in the thickness direction inside the wall of the reaction vessel.
  • the computer program is measured at least at a plurality of temperature measurement points arranged in the thickness direction inside the wall of the reaction vessel.
  • FIG. 1 is a diagram showing a schematic configuration of an operation management device for a reaction vessel according to the first embodiment.
  • FIG. 2 is a flowchart for explaining the inverse problem analysis processing.
  • FIG. 3 is a diagram for explaining an arrangement relationship of thermocouples in an example of the first embodiment.
  • 4A to 4E are diagrams for explaining the analysis results in the example of the first embodiment.
  • FIG. 5 is a diagram showing a schematic configuration of an operation management device for a reaction vessel according to the second embodiment.
  • FIG. 6 is a diagram showing a schematic configuration of another operation management device for a reaction vessel in the second embodiment.
  • FIG. 7 is a diagram for explaining an arrangement relationship of thermocouples in an example of the second embodiment.
  • 8A to 8D are diagrams for explaining the analysis result in the example of the second embodiment.
  • FIG. 9 is a diagram for explaining the arrangement relationship of thermocouples in another example of the second embodiment.
  • FIGS. 10A to 10D are diagrams for explaining analysis results in another example of the second embodiment.
  • FIG. 11 is a diagram showing a schematic configuration of a reaction vessel operation management device according to the third embodiment.
  • FIG. 12 is a diagram for explaining the analysis by the stationary method.
  • FIG. 13 is a diagram for explaining the analysis by the stationary method.
  • FIG. 14 is a diagram for explaining a state in which the amount of heat input suddenly changes.
  • FIG. 15 is a diagram for explaining the problem setting.
  • FIG. 16 is a characteristic diagram showing a change in heat flux at each cross-sectional position when a heat flux that changes stepwise at a position of 4 m is given.
  • FIG. 17 is a characteristic diagram showing a temperature change at each cross-sectional position when a heat flux that changes stepwise at a position of 4 m is given.
  • FIG. 18 is a diagram showing an example in which the horizontal axis is mapped as the non-stationary index and the vertical axis is mapped as the steady heat flux.
  • FIGS. 19A to 19E are diagrams for describing the analysis results in the example of the third embodiment.
  • FIG. 20 is a diagram showing an example in which the horizontal axis represents the unsteadiness index and the vertical axis represents the constant heat flux in the example of the third embodiment.
  • FIG. 21 is a flowchart for explaining the inverse problem analysis processing.
  • FIG. 22 is a diagram for explaining an arrangement relationship of thermocouples in Example 1 of the third embodiment.
  • FIG. 23 is a diagram for explaining an arrangement relationship of thermocouples in Example 1 of the third embodiment.
  • FIGS. 24A and 24B are diagrams showing the results obtained by performing a forward problem analysis with artificial boundary conditions and determining the change over time in thermocouple temperature.
  • Fig. 25A and Fig. 25B are diagrams showing the results of performing a forward problem analysis with artificial boundary conditions and determining the change over time in thermocouple temperature.
  • FIG. 26 is a diagram schematically showing a relationship among a reference point, a known temperature point, and an estimated point.
  • FIG. 27 is a diagram schematically showing a relationship among a reference point, a known temperature point, and an estimated point.
  • Figure 28 shows the results of comparing the heat flux given by the forward problem analysis with the heat flux found by the inverse problem analysis.
  • Figure 29 shows the results of comparing the heat flux given by the forward problem analysis with the heat flux found by the inverse problem analysis.
  • FIGS. 30A to 30D are diagrams for explaining the analysis results in Example 1 of the fourth embodiment.
  • FIGS. 31A to 31D illustrate the analysis results in Example 1 of the fourth embodiment.
  • FIG. 31A to 31D illustrate the analysis results in Example 1 of the fourth embodiment.
  • FIG. 32 is a diagram schematically showing a relationship among a reference point, a known temperature point, and an estimated point.
  • FIG. 33 is a diagram for explaining an analysis result in Example 2 of the fourth embodiment.
  • FIG. 34 is a diagram for explaining the analysis result in Example 2 of the fourth embodiment.
  • FIG. 35 is a diagram showing a two-dimensional cross section of a reaction vessel in which a plurality of thermocouples are embedded.
  • FIG. 1 shows a schematic configuration of the operation management device for a reaction vessel of the present embodiment.
  • the operation management device for the reaction vessel has an input section 101 for inputting temperature data measured by a thermocouple embedded in the wall of the reaction vessel (see Fig. 35).
  • Inverse problem analysis to determine the temperature or heat flux on the inner and outer surfaces of the reaction vessel by performing an inverse problem analysis using the transient heat conduction equation based on the temperature data input to the input unit 101 Unit 102 and an output unit 103 for displaying the temperature or heat flux on the inner surface and the outer surface of the reaction vessel calculated by the inverse problem analysis unit 102, for example, on a display (not shown).
  • Equation (1) The unsteady heat conduction equation used for the inverse problem analysis is expressed as shown in the following equation (1).
  • p is the density of the material in the reaction vessel
  • c P is the specific heat of the material in the reaction vessel
  • ⁇ is the calculated value of the temperature inside the reaction vessel
  • t is time
  • k is the thermal conductivity of the material in the reaction vessel.
  • Heat conduction inverse problem analysis is based on the unsteady heat conduction equation that governs the calculation domain. Then, assuming the temperature inside the region as known, estimating boundary conditions such as temperature and heat flux at the region boundary. On the other hand, the heat conduction order problem analysis refers to estimating the temperature inside a region from known boundary conditions such as temperature and heat flux at the region boundary.
  • Examples of the two-dimensional inverse problem analysis method include, for example, a method disclosed by the present applicant in Japanese Patent Application Laid-Open No. 2002-2066958. It can also be applied to title analysis.
  • As an example of one-dimensional inverse problem analysis an analysis method proposed by Beck et al. Is known (Beck and others, Inverse Heat Conduction, 1985, Wiley, New Yo: rk).
  • the concept shown in Japanese Patent Application Laid-Open No. 2002-209658 is used as a method of analyzing the inverse problem. That is, as shown in the following equation (2), the temperature Y measured by each thermocouple arranged in a certain one-dimensional direction (such as la ⁇ lb ⁇ lc or ld ⁇ le shown in Fig. 35) and the reaction Y
  • the assumed value that minimizes the sum of the squares of the difference from the temperature T at each thermocouple position calculated by the transient heat conduction equation from the assumed value of the heat flux on the inner and outer surfaces of the vessel is the inner surface of the reaction vessel.
  • t Incidentally obtained as heat flux at the outer surface, J is the number of thermocouples.
  • Equation (3) shows an example of a zero-order regularization term.
  • p is the number of divisions of the estimated heat flux, ⁇ .
  • I a regularization parameter obtained from empirical values. 0 ⁇ () (3)
  • Equation (5) represents the objective function representing the deviation of the measured temperature Y and calculated temperature T.
  • Equation (6) is an objective function added to stabilize the calculation, and has the function of suppressing rapid changes in values in the space division direction.
  • ⁇ in equation (6) are regularization parameters obtained from certain empirical values.
  • the objective function is set so that the square of the difference between the temperature Y measured by a certain thermocouple and the temperature ⁇ calculated from the assumed value of the heat flux by the heat conduction equation model is minimized. ing.
  • an objective function for performing regularization in the spatial direction is set so that the solution is stable even if there is a temperature measurement error.
  • equation (4) Is used as the overall objective function, as shown in the following equation (7), a minimum point is searched for the unknown heat flux division region.
  • Equation (8) indicates the heat flux, and it is assumed that the heat flux in the future time m + r _ 1 time step is constant from the heat flux q m force in the m time step.
  • i which indicates the number of divisions in the time direction, changes up to the maximum M time steps
  • j which indicates the number of thermocouples, changes up to J
  • k which indicates the number of divisions of the heat flux distribution
  • T * is a temperature reference value
  • q * is a heat flux reference value.
  • Equation (9) The number of divisions of the heat flux distribution Equation (9) is a simultaneous equation for estimating the change in heat flux when a temperature change occurs.
  • this equation (9) Is used to determine the heat flux q at both ends.
  • the calculated temperature at the thermocouple position in the previous time step is set as the initial T *, and q is obtained by Eq. (9).
  • This q is given as the boundary condition of the forward problem heat conduction equation model calculated in parallel, and the temperature distribution is calculated.
  • the calculated temperature value obtained here is used as the next temperature reference value T *, and q is revised again (substituting into equation (9) to find q again).
  • Equation (10) represents a kind of sensitivity matrix. In short, it shows the ratio of the magnitude of the change of the calculated temperature T at the thermocouple position to the unit change of the heat flux q at the boundary endpoint. are doing.
  • the value per unit time step can be calculated at each time step by the forward problem calculation performed at the same time as the inverse analysis.
  • a more desirable solution will be described using a one-dimensional inverse problem analysis as an example. As described above, even if a one-dimensional inverse problem with the heat flux of the two end faces (the inner surface and outer surface of the reaction vessel) as unknown boundary conditions is constructed (formulated), in principle, a solution must be obtained. Can be done.
  • thermocouples there may be multiple solutions depending on the number of thermocouples and the thermophysical conditions of the materials, and the calculation may be unstable.
  • One of the reasons is that if the combination of “heat flux difference between unknown end faces” can be properly selected, there may be countless combinations of heat flux that express the temperature change at discrete temperature measurement points. Because there is. In particular, heat transfer In the case of a substance with low conductivity, even when estimating the boundary condition that causes the surface temperature to become extremely large or small, if only the temperature change at discrete measurement points is reproduced, one Recognition as a solution can also occur. This not only is impossible as a real phenomenon, but also makes the inverse problem calculation very unstable.
  • the temperature of the thermocouple at the start of the inverse problem analysis (the temperature at the discrete measurement points) is given as a known value, but the initial conditions of the temperature distribution in other analysis domains are unknown. Is common. For this reason, the calculation is started from the tentative initial temperature distribution arbitrarily given, and while the calculation steps are proceeding, the actual temperature distribution is searched and estimated, and the temperature distribution is corrected gradually to an appropriate temperature distribution, and stably It is necessary to use a calculation logic that allows the calculation to proceed (for example, the temperature distribution referred to here means, in the calculation procedure of inverse problem analysis, the solution of the above equation (9) must be modified. This is the calculated value of the forward problem heat conduction equation model calculated in parallel).
  • a temporary heat flux q is given as a heat flux on one of the inner surface and the outer surface of the reaction vessel, here, on the outer surface.
  • T surf indicates the temperature at the unknown boundary, here the outer surface of the reaction vessel.
  • This surface temperature T surf is usually also used to perform a forward problem analysis in order to correct the heat flux value during the inverse problem analysis, but it is equivalent to the surface temperature obtained in the forward problem analysis.
  • the reference temperature Tb is a temperature outside the inside and outside and inside surfaces of the reaction vessel.
  • the cooling condition of the reaction vessel is determined based on the water temperature and the like in the case of water cooling.
  • the heat flux q on the left side of the above equation can be given as if it were known heat flux information.
  • the heat transfer coefficient h and / or the external reference temperature T b are changed by several points, that is, the value of the temporary heat flux q on the outer surface of the reaction vessel is shaken by several points (point K), and the reaction is performed.
  • the temporary heat flux q on the outer surface of the vessel and the heat flux on the inner surface of the reaction vessel that minimizes the sum of the squares of the differences in temperature ⁇ and ⁇ when given each temporary heat flux information q Are obtained (Step S203).
  • the one-dimensional inverse problem analysis for simultaneously obtaining the heat flux changes on the inner surface and the outer surface of the reaction vessel can be stably executed. If the temperature change and the heat flux change on the inner surface and the outer surface of the reaction vessel can be simultaneously estimated, for example, the temperature fluctuation at a certain temperature measurement point is caused by the heat flux change on the inner surface of the reaction vessel. It is possible to distinguish whether the heat flux is caused by a change in heat flux on the outer surface of the reaction vessel caused by poor contact of a cooling device installed outside the reaction vessel or the like.
  • the above method is simple to apply to one-dimensional inverse problem analysis, and is often effective as a practical problem.
  • the reason is that the upper and lower ends of the reaction vessel are generally insulated (symmetric) in many cases, and there is no practical problem.
  • inverse problem analysis can be performed assuming one dimension in the thickness direction within the area demarcated by the broken line in Fig. 35, and the results can be combined into two dimensions by combining them vertically.
  • a two-dimensional inverse problem analysis is required.
  • Such two-dimensional analysis consists of making the heat flux divisions at the left and right ends in Fig. 35 finer upward, and finding the minimum heat flux distribution by least squares at the temperature at these thermocouple positions. This is equivalent, and the present invention may be applied in accordance with the same method as the inverse problem formulation disclosed in Japanese Patent Application Laid-Open No. 2002-209658.
  • the heat flux at the upper and lower ends in Fig. 35 may be unknown or known, but considering the stability of calculation, an appropriate heat flux (for example, adiabatic Etc.) and it is desirable to make it known. Based on the same idea, extension to three-dimensional analysis can be easily performed.
  • thermocouples 1 t and 1 s are embedded at 0.1 [m] and 0.2 [m] from the outer end of the furnace. ing.
  • the position of the thermocouple is arranged to be deviated to the outside of the furnace.
  • FIGS. 4A to 4E show measured temperature data and analysis results of the thermocouples 1 s and 1 t.
  • the horizontal axis is the number of days.
  • Figure 4A shows the changes over time of the two thermocouples 1 s and It. According to the results, two places where the temperature is noticeably higher can be observed as enclosed by the ellipse (high temperature 1 and high temperature 2).
  • Figures 4B and 4C show the results of calculating the temperature at the furnace inner end (the inner surface of the reactor) and the temperature at the furnace outer end (the outer surface of the reactor) using the inverse problem analysis method described above. .
  • the high-temperature phenomenon 1 shows that the inner and outer ends of the furnace were simultaneously heated, and the effect of the reaction activation high temperature in the furnace affected the outer end of the furnace. It can be inferred.
  • the high temperature phenomenon 2 almost no temperature fluctuation was observed at the inner end of the furnace, and it can be inferred that the temperature was raised simply because the cooling capacity outside the furnace was reduced.
  • Fig. 4'D shows the change of the heat flux at the outer end and the inner end of the furnace similarly obtained by the inverse problem analysis method.
  • Fig. 4E shows a characteristic diagram converted to the heat transfer coefficient h on the outer end of the furnace. From these results, the heat transfer coefficient h also fluctuated greatly near the high temperature phenomenon 2, and the cooling capacity gradually increased after about 300 days, and dropped sharply near the high temperature phenomenon 2. It can be understood that there is an abnormality in cooling due to some factor.
  • FIG. 5 shows a schematic configuration of the operation management device for the reaction vessel of the present embodiment.
  • reference numeral 201 denotes an arithmetic unit, which performs an operation to increase the number of digits after the decimal point using the temperature data input to the input unit 101, as described later.
  • the inverse problem analysis unit 103 performs an inverse problem analysis using an unsteady heat conduction equation based on the temperature data in which the number of digits after the decimal point has been increased by the arithmetic unit 201, thereby obtaining a reaction vessel.
  • the unit time step of the inverse problem calculation is long, it is necessary to estimate a change that goes far back in the past at least for that time step.
  • the inverse problem calculation is a calculation result obtained by time averaging in the calculation time step, the restored heat flux change also captures a dull change, and a rapid change below the time step cannot be captured.
  • shortening the unit time step of the inverse problem calculation means that during that short time, the small amount of heat transmitted from the position where the heat flux changes (for example, the inner surface of the blast furnace in contact with the hot metal) to the thermocouple position This means that temperature changes must be captured.
  • the position of the thermocouple is far from the position where the heat flux changes, and the thermal conductivity of the material between them is small, the temperature movement in a short time becomes very small.
  • thermocouple Therefore, the present inventors have conducted intensive studies in order to be able to estimate the past heat flux distribution as close as possible in a short time step while utilizing temperature data measured by a thermocouple.
  • the global temperature When looking at the flow of data, it is common to evaluate the value by rounding down the temperature data after the decimal point.
  • the manufacturer's guaranteed range of the temperature measurement device that converts the thermocouple contact voltage change into temperature may be up to one decimal place in most cases. Often temperature data is used. In other words, the collected and recorded temperature data is coarse, with zero or one decimal place.
  • the time transition of temperature-the flow of temperature change is more important than the accuracy of temperature itself. If the temperature changes in a staircase, pursuing only the accuracy of the temperature itself, the inverse problem analysis would be unstable. In particular, when aiming for inverse problem analysis in a short time step, an unpredictable step-like temperature change may occur in that time step. Therefore, even if an inverse problem analysis is performed using this stepwise temperature change to obtain a heat flux solution, there is a possibility that a physically impossible solution may be reached. In fact, in most cases the solution diverges, making it impossible to continue the inverse problem analysis.
  • Generating temperature data with an increased number of decimal places from temperature data with a small number of decimal places in a rational manner is an effective method for smoothing step-like temporal changes. It is.
  • the arithmetic unit 201 collects temperature data with zero or one digit below the decimal point with a sampling time shorter than the time step in inverse problem analysis, averages them over time, and performs inverse problem analysis.
  • the representative value of the temperature data used in the calculation time step in. For example, a small digit number of temperatures sampled at 5-minute intervals shorter than the time step (1 hour) in the inverse problem analysis is simply averaged over a 1-hour range, and a representative value of the 1-hour-step temperature data is obtained. I do. This makes it possible to apparently reduce the number of decimal points in the one-hour temperature data.
  • the temperature data used for the inverse problem calculation is based on the temperature data (measured by the thermocouple) with zero decimal places and one decimal place that guarantees the accuracy measured by the thermocouple.
  • the number of decimal places of It is possible to increase.
  • the number of digits after the decimal point of the temperature data used for the inverse problem calculation may be increased by another method. For example, in a thermocouple, even if the assurance range of the temperature data itself is one decimal place, it is actually calculated based on the conversion formula that converts the voltage at the thermocouple junction into temperature, or the calibration is performed. A common method is to convert the voltage-temperature conversion table to temperature.
  • the calculated temperature value of many digits is rounded off to one digit after the decimal point, and it is actually possible to obtain the number of decimal digits according to the capacity of the computer.
  • a conversion formula for converting the voltage of a K-type thermocouple into temperature there is the following formula.
  • T a + b ⁇ X + c ⁇ ⁇ 2 + d ⁇ X 3 + e ⁇ X + f ⁇ X 5 + g ⁇ X 6 + h ⁇ X 7 + i ⁇ X 8
  • the contact voltage is a true contact voltage after cold junction temperature compensation.
  • a filter 101a is provided in the input section 101 as shown in FIG. To be applied. That is, using the filter 10 la, Correct the temperature data measured in pairs and use the corrected temperature data for inverse problem analysis. If a filter is used, even if the original temperature data has a small number of decimal places, the decimal places can be apparently generated in the temperature data. Calculation of problem analysis is very stable.
  • the mouth-pass filter corrects the measured data in order to remove the influence of fluctuation due to high-frequency noise mixed in the measured data sampled at equal time intervals.
  • thermocouple position caused by a change inside the furnace far away is a low frequency change in time.
  • high-frequency noise due to other causes is superimposed on the actual thermocouple contact voltage change (original signal before temperature conversion). Therefore, using temperature data that has been processed to remove ⁇ -frequency noise from actual measured temperature changes in an inverse problem analysis is a very physically meaningful method. Therefore, this method not only can improve the calculation stability of the inverse problem analysis, but also can improve the accuracy of estimating the heat flux.
  • the number of digits after the decimal point was increased by some other method, as well as the temperature data whose number of digits was increased by taking the above-mentioned time average or using a conversion formula.
  • Equations (12) to (14) show the equations of the three-term method, the five-term method, and the seven-term method, respectively, and the performance as a low-pass filter increases as the number of measured temperature data used for correction increases. Is improved. That is, only lower frequency signals can be extracted. Basically, it can be expressed by the coefficient of each term. 5) It becomes like (17).
  • FIG. 7 shows a one-dimensional transient heat conduction model using a thermocouple embedded in a carbon brick installed on the bottom of the blast furnace bottom.
  • TC 1 indicates a high-temperature thermocouple
  • TC 2 indicates a low-temperature thermocouple. From the temperature changes measured by these thermocouples TC1 and TC2, the unsteady heat flux q1 on the high-temperature heat flux surface and the unsteady heat flux q2 on the cooling surface are calculated by the inverse problem analysis described above. Estimate at the same time.
  • the position of the high-temperature heat flux surface is a fixed point 4.0 m deep from the cooling surface.
  • FIG. 8A shows the estimation result of the heat flux q1 on the high temperature side obtained in the present example.
  • the horizontal axis of 8D represents the date, and shows the results from September 1 to October 1. Actually, the calculation starts before September 1, and ends on October 1st.
  • the time step of the inverse problem analysis is 8 time steps.
  • Case 1 rounds off the temperature data measured by the thermocouple to one decimal place. This is the calculation result when filter processing is performed on the data with 0 decimal places.
  • Case 2 is the calculation result when filtering the temperature data with two decimal places rounded to one decimal place.
  • FIG. 8B and 8C show plots of the temperature data (temperature data used for inverse problem analysis) after filtering the high-temperature-side thermocouple TC1 and low-temperature-side thermocouple TC2.
  • FIG. 8D shows the temperature difference between case 1 and case 2 (case 1—case 2) at the respective positions of thermocouples T C1 and T C2.
  • the Spencer-type equation (Equation (2 1)) of up to 21 terms is used.
  • a technique is used in which the number of terms is gradually reduced as approaching.
  • the expressions (15) to (21) are used in a Windows format notation. That is, the number of terms decreases as the date approaches October 1, and finally equation (15) is used, and no filtering is performed on the last temperature data.
  • Fig. 9 shows a schematic diagram similar to Fig. 7 for cases where the thermocouples are installed at different depths.
  • Figure 10A shows the estimation results of the high-temperature side heat flux q1 when the calculation time step of the inverse problem analysis is set to 8 hours and when it is set to 6 hours. The results of these calculations are extracted from September 1 to October 1 of the long-term calculation results.
  • FIGS. 10B and 10C show characteristic diagrams obtained by plotting the temperature data of the high-temperature-side thermocouple TC1 and the low-temperature-side thermocouple TC2 after the smoothing. Further, FIG. 10D shows the difference between the high-temperature side thermocouple temperature and the low-temperature side thermocouple temperature (T C 1 -T C 2) corresponding to the steady state heat flux.
  • the quasi-stationary heat flux in the characteristic diagram in Fig. 10D can be interpreted as the average heat flux passing around the thermocouple position, and the value of the unsteady heat flux estimated at the hot surface 4 m deep
  • the peaks do not always correspond to each other, but in this case, a peak lag of about 2-3 days is observed.
  • the unsteady heat flux based on the inverse problem analysis detects the movement in the furnace approximately 2-3 days earlier than the conventional steady-state method.
  • the 6-hour step calculation shows a movement with a sharper contour. If the time step is lengthened, only the averaged motion within that time can be captured, so shortening the time step as much as possible is more likely to express the actual internal motion.
  • FIG. 11 shows a schematic configuration of the operation management apparatus for a reaction vessel of the present embodiment.
  • reference numeral 301 denotes an analysis unit based on the steady-state method. By performing analysis based on the temperature data obtained through the input unit 101 by the steady-state method, the temperature on the inner surface of the reaction vessel is calculated. Or heat flux.
  • Reference numeral 302 denotes a comparison unit which compares the temperature or heat flux obtained by the inverse problem analysis unit 102 with the temperature or heat flux obtained by the analysis unit 301 by the steady-state method.
  • the output unit 103 displays the result compared by the comparison unit 302 on, for example, a display (not shown). Since the inverse problem analysis used in the present embodiment has been described in the first embodiment, the description is omitted here.
  • the temperature at the inner surface of the reaction vessel is obtained from the temperature measurement values at two points in the depth direction (thickness direction of the wall (brick)) in the one-dimensional stationary method.
  • the temperatures at the two points are connected by a straight line, and the temperature at the outside point is the temperature on the inner surface of the reaction vessel.
  • the heat flux is the slope of the straight line, it can be obtained by multiplying the slope by the thermal conductivity k (the intermediate term in the following equation (22)).
  • a straight line that satisfies the temperature at the plurality of points is obtained by the linear least squares method.
  • the temperature at the outside point is defined as the temperature on the inner surface of the reaction vessel.
  • the heat flux can be obtained by multiplying the slope of the straight line by the thermal conductivity k (the rightmost term in the following equation (22)).
  • the heat flux (steady heat flux) obtained as a result of the analysis by the steady-state method can be expressed by a straight line, so that it takes the same value regardless of the cross-section, which is the heat flux (unsteady heat flux) obtained as a result of the inverse problem analysis. Bunch).
  • the temperature gradient changes sensitively due to the effect of the outside, so it is easy to largely separate from the straight line of the steady method.
  • Fig. 14 for example, if the inflow calorific value increases suddenly, the temperature gradient will increase locally, and if the inflow calorie decreases suddenly, the temperature gradient will decrease locally. If such a change is not instantaneous but extends for a long time (steady state), the distribution of the temperature gradient is gradually eliminated and converges on a certain straight line (steady state).
  • the implication of estimating the heat flux on the inner surface of the reaction vessel by inverse problem analysis using the unsteady heat conduction equation is to accurately capture this local change in the temperature gradient. It is important to be able to make decisions.
  • the unsteady heat flux is lower than the steady heat flux, it can be estimated that the steady heat flux (average heat flux) will tend to decrease in the future. Conversely, if the unsteady heat flux is higher than the steady heat flux, it can be estimated that the steady heat flux (average heat flux) will tend to increase in the future. Also, if the heat flux difference is large, the unsteadiness is large, so it can be judged that this suggests that the average heat flux will rapidly rise and fall in the future.
  • the same evaluation can be performed by comparing the solution of the inverse problem analysis and the analysis by the steady-state method for the estimated temperature at the fixed position, instead of comparing the heat flux at the fixed position. Also, by examining the starting point where the heat flux difference between the two begins to change rapidly, it is possible to identify what caused the sudden rise or fall.
  • the difference between the unsteady heat flux (or temperature) at a fixed position and the steady heat flux (or temperature) is used as an indicator of unsteadiness, and There is also a method of evaluating The non-stationarity index is calculated by the following formulas (23), (24) It is possible to define the heat flux difference or the temperature difference at a fixed position as a numerical value, so that a quantitative grasp can be obtained. It is important to clarify at which fixed position the value of the non-stationarity index is defined, and the value changes depending on the defined position, as described above.
  • the fixed position reference is fixed at a very shallow position, the unsteadiness may become dull as shown in Fig. 16, and the desired comparison result may not be obtained. In such a case, it is desirable to set the standard for the reactor with the largest brick thickness, even if it is virtually.
  • the non-stationary index only estimates the magnitude and directionality of the subsequent changes based on the current steady-state heat flux (steady-state temperature), and therefore depends on the magnitude of the heat flux at that time. The meaning changes. That is, the relationship between both the unsteadiness index and the absolute value of the heat flux is important.
  • Fig. 18 schematically shows an example of mapping that concept.
  • the horizontal axis shows the unsteadiness index ⁇ ⁇ ⁇ vine ⁇
  • the vertical axis shows the steady heat flux, and this operation requires that the steady heat flux be kept above a certain reference value Q.
  • the critical value of the unsteadyness index ⁇ T tran is set to 120 ° C.
  • the value depends on the concept of stabilizing the operation, the allowable heat flux change, and the thermophysical properties of the reaction vessel wall. Of course, it changes depending on the value and the estimated position (thickness position). Also, here, the value on the vertical axis is the steady heat flux, but there is no problem even if it is unsteady heat flux.
  • FIG. 9 shows a case where a thermocouple embedded in a carbon brick installed in the bottom of the blast furnace bottom is used.
  • An example of an inverse problem analysis is shown, assuming a one-dimensional transient heat conduction model. Two thermocouples are biased toward the cooling surface side. TC 1 indicates a high temperature side thermocouple, and TC 2 indicates a low temperature side thermocouple. From the temperature changes measured by these thermocouples TC1 and TC2, the unsteady heat flux q1 on the high-temperature heat flux surface and the unsteady heat flux q2 on the cooling surface were simultaneously determined by the inverse problem analysis described above. presume. The position of the high-temperature heat flux surface is a fixed point 4.
  • Figures 19C and 19D show the transition of the temperature at the high-temperature thermocouple TC1 and the temperature at the low-temperature thermocouple TC2 (temperature data used for the inverse problem analysis).
  • the horizontal axis is the date Yes, showing results from September 1 to January 25.
  • Figure 19A shows the results of the heat flux by the inverse problem analysis on the high-temperature heat flux surface (4 m position) estimated using the above temperature data, and the results of the heat flux by the steady-state method.
  • the transition of the unsteady heat flux by the inverse problem analysis shows larger fluctuations, while the steady heat flux by the steady method changes slowly. Comparing these changes, it can be inferred that from around January 8, the difference in heat flux between the two methods increased, and that there was an unsteady change in the hot metal flow in the furnace.
  • Figure 19B shows the results of the temperature estimation by the inverse problem analysis on the high-temperature heat flux surface (at a position of 4 m) estimated using the above temperature data, and the temperature estimation by the steady-state method. From this result, it is possible to read the same tendency as in the case of heat flux.
  • Figure 19E shows the change in the unsteadiness index ⁇ tran tran at the high-temperature heat flux surface (4 m position) as a function of temperature. Expressing with such indicators makes it easier to see when the change originated.
  • FIG. 20 shows the result of mapping this calculation result between the unsteadiness index ⁇ T tran and the steady heat flux.
  • the results show that the critical value of the unsteadyness index ⁇ T tra ⁇ was set to 130 ° C. with the lowest reference value of the steady-state heat flux set at 350 ° W / m 2 , whereas the critical value of 130 ° C. In this case, the steady-state heat flux decreases rapidly after this value is cut, and eventually falls significantly below the reference value. Did not take any action.
  • Unsteady index delta T t ra n are, so show a deviation from stationarity inside fixed point of the furnace. It can be grasped signs of early changes than before. Therefore, if an appropriate action standard can be set by utilizing the characteristics of this indicator, it will lead to a reduction in production costs, but it will be possible to promote marginal operation, which is a bad operating condition. It enables agile operation design to detect signs of deterioration in the furnace at an early stage and take immediate recovery action.
  • the schematic configuration of the reaction vessel operation management device of the present embodiment is the same as that of FIG. 1 described in the first embodiment, but is input to the input unit 101 in the inverse problem analysis unit 102.
  • the temperature or heat flux on the inner and outer surfaces of the reaction vessel by performing inverse problem analysis using an internal function or an external function that satisfies the transient heat conduction equation based on the temperature data It is characterized by the following.
  • the boundary condition on the inner surface and outer surface of the reaction vessel as a simpler method is obtained.
  • the change is estimated.
  • the interpolation function or the outer ⁇ function is a function that connects the temperature at the measurement point and expresses a region other than that point, for example, the whole or a part of the analysis region.
  • Linear functions such as linear function approximation and spline interpolation are known as internal functions that cannot be extrapolated. However, functions that can be extrapolated while satisfying the transient heat conduction equation are not known.
  • Interpolation refers to estimating unknown points inside known points
  • extrapolation refers to estimating outside and around known points.
  • the inverse problem analysis unit 102 first expresses the solution of the unsteady heat conduction equation using a predetermined internal function or external function and parameters (step S 2 101) .
  • t represents time
  • X, y, and z represent position vector elements, and can be applied to a general three-dimensional coordinate system.
  • ⁇ ⁇ , ⁇ " ⁇ ⁇ , ⁇ ⁇ , A y , ⁇ ⁇ , X, ⁇ , and ⁇ represent appropriate arbitrary constants, and the optimum value varies depending on the target system. The values of these arbitrary constants Care must be taken in choosing
  • This function F ( X , y, z, t) automatically satisfies the transient heat conduction equation (25).
  • this function F (x, y, z, t) is used to express the solution of the unsteady heat conduction equation in general, it is expressed as the following equation (27).
  • x j , y j , and z j represent each element of an arbitrary reference position vector, represents an arbitrary reference time, and x , y, z, and t are points at which temperature is to be estimated.
  • N j is the number of reference position vectors and the number of reference times in the time direction, respectively.
  • a is a parameter, but if this value is determined, it is possible to determine an arbitrary position vector (X, y, z) and a temperature distribution T (X, y, z, t) at time t. You can.
  • the value of the parameter in the solution of the transient heat conduction equation expressed by the above equation (27) is determined using the temperature information measured by the thermocouple (step S2102). The value of this parameter a can be determined by solving the following simultaneous equation (28).
  • a ki l is the temperature T measured by the thermocouple. (x k , y k , z k , t 1 ), where the subscript k indicates the measurement position and the subscript 1 indicates the sampling time.
  • the heat conduction inverse problem is based on the unsteady heat conduction equation governing the calculation domain, assuming that the temperature inside the domain is known, and the boundary conditions such as temperature and heat flux at the domain boundary or the initial conditions.
  • the heat conduction order problem refers to the problem of estimating temperature information inside a region based on known boundary conditions.
  • the temperature distribution at the wall boundary of the reaction vessel is also estimated at the same time, and although indirect, the boundary between the inner surface and the outer surface of the reaction vessel is obtained from the temperature information measured by the thermocouple. This is the inverse problem of deciding the conditions.
  • the temperature gradient at the boundary can be estimated not only from the temperature distribution at the wall boundary of the reaction vessel but also from the temperature distribution in the vicinity, the heat flux at the wall boundary position of the reaction vessel is consequently obtained.
  • the change can also be estimated.
  • Figs. 22 and 23 show one-dimensional transient heat using thermocouples embedded in carbon bricks installed on the bottom (Fig. 2 2) and side walls (Fig. 2 3) of the blast furnace bottom.
  • a model that attempted an inverse problem analysis assuming conduction is schematically shown.
  • two thermocouples are embedded toward the cooling surface side
  • TC1 is a hot-side thermocouple
  • TC2 is a low-temperature-side thermocouple.
  • the inverse problem analysis was used to determine the unsteady heat flux q1 on the hot heat flux surface and the unsteady heat flux q2 on the cooling surface. Estimate at the same time.
  • Fig. 22 and 23 show one-dimensional transient heat using thermocouples embedded in carbon bricks installed on the bottom (Fig. 2 2) and side walls (Fig. 2 3) of the blast furnace bottom.
  • TC1 is a hot-side thermocouple
  • TC2 is a low-temperature-side thermocouple.
  • the inverse problem analysis was used to determine the unstead
  • the high-temperature heat flux surface position is a fixed point (4.0 m) at a depth of 4.0 ⁇ 1 from the cooling surface, and in Fig. 23, the high-temperature heat flux surface position is 2.0 m from the cooling surface.
  • the fixed point (x 2.0 m).
  • ⁇ x it is also an important factor to set the reference position vector and the reference time.
  • Figures 26 and 27 show the reference point (reference position vector and reference time), the known temperature point (the thermocouple position vector whose temperature is being measured and the known temperature time), and the estimated point (temperature estimation position position). The relationship between the distance and the estimated time is schematically shown. In the examples of Figs. 26 and 27, the reference point and the known temperature point are matched, and Fig. 26 shows that three temperature points are used in the time direction and Fig.
  • FIG. 27 shows that five temperature points are used in the time direction. I have.
  • the estimation point selects the point closest to the present in the time direction (shown as "present"), and sets the position vector to the q1 and q2 positions did.
  • the estimated points are indicated by "X”.
  • the heat fluxes q1 and q2 at the estimation point are obtained with the settings in Fig. 26.
  • the heat flux The heat fluxes q 1 and q 2 are obtained, and ⁇ ⁇ ⁇ ⁇ that can reproduce the boundary conditions relatively well compared to the heat fluxes 1 and q 2 in Figs.
  • the value was determined by solving the simultaneous equations in Eq. (28).
  • Eq. (27) the temperature distribution of the entire region is estimated, so the values of ql and q2 cannot be directly obtained. Therefore, the temperature (T p ) at the position vector (x p ) of the estimated point and the temperature (T p at 3. O mm inside of X p , c ) near the estimated point pole ( ⁇ ⁇ ) p , c ) was estimated and calculated by the following equation (31), and it was assumed that the heat flux q was at the estimated point.
  • Fig. 28 shows the results of comparing the heat fluxes q1 and q2 given by the forward problem analysis with the heat fluxes q1 and q2 found by the inverse problem analysis for the model in Fig. 22.
  • the solid lines are q1 (thick line) and q2 (thin line) given by the forward problem, and the plots are the estimated values of q1 (command) and q2 (country) obtained by the inverse analysis.
  • Time step of inverse analysis is 28800 seconds, the value of tau chi at this time was 1800000 seconds.
  • the value of ⁇ depends on the unit system of the analysis, but in this analysis, the MKS unit system is used.
  • Fig. 29 shows the results of comparing q1 and q2 given in the forward problem analysis with q1 and q2 found in the inverse problem analysis for the model in Fig. 23.
  • the solid lines are q 1 (thick line) and q 2 (fine line) given by the forward problem, and the plots are the estimated values of q 1 ( ⁇ ) and q 2 (garden) obtained by the inverse analysis.
  • thermocouple TC1 temperature data in the characteristic diagrams of Figs. 30C and 31C thermocouple TC2 temperature data in the characteristic diagrams of Figs. 30D and 31D
  • the characteristics of Fig. 30A and 31A The heat flux ql on the high temperature side in the figure and the heat flux q2 on the low temperature side in the characteristic diagrams in Figs. 30B and 31B are simultaneously estimated.
  • the ⁇ ordinary method '' is a method of estimating the once-through heat flux that has been conventionally used, and multiplies the absolute value of the temperature difference between TG 1 and TC 2 by the thermal conductivity k x to calculate the distance between the TC 1 position and the TC 2 position. It is the value re-divided by the absolute value of In particular, it is noticeable in Figs. 31A to 31D, but it can be seen that ql of the steady-state method changes in the same shape several days later than the ql estimation. This is the effect of the lag time due to the effect of the heat capacity (p CP ) of the blast furnace carbon brick.
  • Figure 32 schematically shows the relationship between the reference point, the known temperature point, and the estimated point.
  • the future temperature change is estimated in the time direction with the position of the estimation point being the position of TC1 and the position of TC2.
  • the unit time step of the temperature known point (past) and the estimated point (future) is
  • the reference point is set by shifting the half-hour step by 14400 seconds, and the time step of the reference point is doubled to 57600 seconds.
  • Figure 33 shows an example of the analysis results.
  • the usual implicit difference approximation method time step 28,800 seconds was used.
  • the plot ( ⁇ , garden) in Fig. 33 shows the estimation results by this method, using the six known temperature points (three TC1 and three TC2) immediately before the plot. This is the result of estimating six estimated points (future points) (see Figure 32). Overall, it has very good estimation performance. However, it can be seen that some points show estimates in completely different directions.
  • the operation management device for a reaction vessel is configured by a computer CPU or MPU, RAM, ROM, or the like, and is realized by the operation of a computer program stored in RAM or ROM. Therefore, the computer program itself for realizing the functions of the above-described embodiment on a computer implements the functions of the above-described embodiment, and constitute the present invention.
  • means for supplying the computer program to a computer for example, a recording medium storing the computer program constitutes the present invention.
  • a recording medium for example, a flexible disk, a hard disk, an optical disk, a magneto-optical disk, a CD-ROM, a magnetic tape, a nonvolatile memory card, a ROM, or the like can be used.
  • the unsteadiness of the state in the reaction vessel is obtained. Strength can be evaluated.

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Manufacturing & Machinery (AREA)
  • Materials Engineering (AREA)
  • Metallurgy (AREA)
  • Organic Chemistry (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

L'invention concerne un dispositif de régulation du fonctionnement d'une cuve à réaction. Ce dispositif comprend une unité (101) de saisie dans laquelle sont introduites les données relatives à la température mesurée par des thermocouples situés sur les surfaces interne et externe de la cuve à réaction, et une unité (102) d'analyse de problème inverse permettant de déterminer les températures ou les flux thermiques au niveau des surfaces interne et externe de la cuve à réaction par une analyse de problème inverse faisant appel à une équation de la chaleur en transitoire utilisant les données de température introduites dans l'unité (101) de saisie. L'unité (102) d'analyse de problème inverse calcule les températures T au niveau des positions des thermocouples en utilisant l'équation de la chaleur en transitoire, à partir des valeurs temporaires des flux de chaleur mesurées au niveau des surfaces interne et externe des cuves à réaction. L'unité (102) d'analyse de problème inverse définit ensuite, en tant que flux de chaleur, les valeurs temporaires des flux de chaleur pour lesquelles la somme du carré des différences entre les températures Y mesurées par les thermocouples disposés de manière linéaire, et les températures T calculées, est minimale.
PCT/JP2003/012997 2002-11-29 2003-10-09 Procede de regulation du fonctionnement d'une cuve a reaction, dispositif, programme, et support d'enregistrement lisible par ordinateur Ceased WO2004050920A1 (fr)

Applications Claiming Priority (6)

Application Number Priority Date Filing Date Title
JP2002348429A JP2004003800A (ja) 2002-04-11 2002-11-29 反応容器の操業管理方法、装置、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体
JP2002-348429 2002-11-29
JP2003-044672 2003-02-21
JP2003044672A JP4299554B2 (ja) 2003-02-21 2003-02-21 逆問題解析方法、装置、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体
JP2003085817A JP4559708B2 (ja) 2003-03-26 2003-03-26 反応容器内の状態を評価する評価方法、操業管理方法、評価装置、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体
JP2003-085817 2003-03-26

Publications (1)

Publication Number Publication Date
WO2004050920A1 true WO2004050920A1 (fr) 2004-06-17

Family

ID=32475218

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2003/012997 Ceased WO2004050920A1 (fr) 2002-11-29 2003-10-09 Procede de regulation du fonctionnement d'une cuve a reaction, dispositif, programme, et support d'enregistrement lisible par ordinateur

Country Status (2)

Country Link
TW (1) TW200415343A (fr)
WO (1) WO2004050920A1 (fr)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104540968A (zh) * 2012-07-25 2015-04-22 技术资源有限公司 启动一种熔炼工艺
CN111723498A (zh) * 2020-05-07 2020-09-29 西安电子科技大学 一种基于微处理器的输油管道内层温度外推监测方法
CN118670553A (zh) * 2024-06-17 2024-09-20 西北大学 一种管道测温温度修正方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06264153A (ja) * 1993-03-09 1994-09-20 Sumitomo Metal Ind Ltd 連続式加熱炉内のスラブ温度予測方法
JPH06265415A (ja) * 1993-03-12 1994-09-22 Kobe Steel Ltd 非定常熱伝導熱流束測定方法
JP2000317594A (ja) * 1999-05-10 2000-11-21 Nippon Steel Corp 溶融金属鋳型内の凝固シェル厚み及びパウダー流入厚み予測方法
JP2001234217A (ja) * 2000-02-28 2001-08-28 Nippon Steel Corp 高炉炉底状況の推定および予測方法
JP2002206958A (ja) * 2001-01-10 2002-07-26 Nippon Steel Corp 湯面高さ検知装置、方法、及びコンピュータ読み取り可能な記憶媒体

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06264153A (ja) * 1993-03-09 1994-09-20 Sumitomo Metal Ind Ltd 連続式加熱炉内のスラブ温度予測方法
JPH06265415A (ja) * 1993-03-12 1994-09-22 Kobe Steel Ltd 非定常熱伝導熱流束測定方法
JP2000317594A (ja) * 1999-05-10 2000-11-21 Nippon Steel Corp 溶融金属鋳型内の凝固シェル厚み及びパウダー流入厚み予測方法
JP2001234217A (ja) * 2000-02-28 2001-08-28 Nippon Steel Corp 高炉炉底状況の推定および予測方法
JP2002206958A (ja) * 2001-01-10 2002-07-26 Nippon Steel Corp 湯面高さ検知装置、方法、及びコンピュータ読み取り可能な記憶媒体

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104540968A (zh) * 2012-07-25 2015-04-22 技术资源有限公司 启动一种熔炼工艺
US20150259758A1 (en) * 2012-07-25 2015-09-17 Technological Resources Pty. Limited Starting a smelting process
US9771626B2 (en) * 2012-07-25 2017-09-26 Technological Resources Pty. Limited Starting a smelting process
CN104540968B (zh) * 2012-07-25 2018-09-14 塔塔钢铁公司 启动一种熔炼工艺
CN111723498A (zh) * 2020-05-07 2020-09-29 西安电子科技大学 一种基于微处理器的输油管道内层温度外推监测方法
CN111723498B (zh) * 2020-05-07 2024-04-09 西安电子科技大学 一种基于微处理器的输油管道内层温度外推监测方法
CN118670553A (zh) * 2024-06-17 2024-09-20 西北大学 一种管道测温温度修正方法

Also Published As

Publication number Publication date
TW200415343A (en) 2004-08-16

Similar Documents

Publication Publication Date Title
EP2113552B1 (fr) Dispositif évaluant la surface des parois d'un four à coke, dispositif d'assistance pour réparer les surfaces d'un four à coke, procédé évaluant la surface des parois d'un four à coke, procédé d'assistance pour réparer les surfaces d'un four à coke et programme informatique
Britsch et al. Natural circulation FLiBe loop overview
WO2004050920A1 (fr) Procede de regulation du fonctionnement d'une cuve a reaction, dispositif, programme, et support d'enregistrement lisible par ordinateur
JP4299554B2 (ja) 逆問題解析方法、装置、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体
JP4695376B2 (ja) 加熱又は冷却特性評価方法及び装置、反応容器の操業管理方法及び装置、コンピュータプログラム、並びにコンピュータ読み取り可能な記録媒体
JP2007071686A (ja) 容器壁の温度又は熱流束の推定方法、装置、コンピュータプログラム、及びコンピュータ読み取り可能な記録媒体、並びに容器壁厚みの推定方法
JP2011257351A (ja) 熱処理シミュレーション方法
JP3978090B2 (ja) 湯面位置検知方法、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体
JP2005082862A (ja) 反応容器の内表面位置推定方法、装置、及びコンピュータプログラム
TWI391676B (zh) 動態截止頻率改變之濾波器之方法及裝置
JP4559708B2 (ja) 反応容器内の状態を評価する評価方法、操業管理方法、評価装置、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体
JP4681127B2 (ja) 湯面高さ検知装置、方法、及びコンピュータ読み取り可能な記憶媒体
JP2002514749A (ja) 冶金炉の温度測定方法
JP4753374B2 (ja) 容器壁の厚み推定方法、装置、コンピュータプログラム
JP5064433B2 (ja) 容器内表面の熱流束の推定方法、装置及びプログラム
WO2014030118A2 (fr) Procédé et système de détermination du profil d'usure réfractaire dans un haut fourneau
JP4833621B2 (ja) 反応容器の温度又は熱流束の推定方法、装置、コンピュータプログラム、及びコンピュータ読み取り可能な記録媒体
JP2002266011A (ja) 高炉の炉内状況推定方法
JP2004003800A (ja) 反応容器の操業管理方法、装置、コンピュータプログラム、及びコンピュータ読み取り可能な記憶媒体
JP4743781B2 (ja) 容器の内壁面の温度及び熱流束の推定方法、装置、並びにコンピュータプログラム
JP4850803B2 (ja) 材料内部の温度推定方法、熱流束の推定方法、装置、及びコンピュータプログラム
JP2008070340A (ja) 超音波を用いた温度測定方法
JP2021130860A (ja) 炉熱予測装置および該方法ならびに炉熱制御案内装置および該方法
JP2002363619A (ja) 高炉下部の制御方法
Hashim et al. Wetting Front Studies on Heated Plate with Orthogonal Water Jet Impingement

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): BR CN KR

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PT RO SE SI SK TR

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
122 Ep: pct application non-entry in european phase