Metallurgical node temperature self-adaptive control method and system based on multi-zone coordinated regulation

Abstract

The present application relates to the technical field of temperature control, and especially relates to a metallurgical node temperature self-adaptive control method and system based on multi-zone cooperative regulation, which comprises the following steps: collecting the temperature of multiple temperature zones in a metallurgical furnace; using the linear distance of the temperature zones to perform spatial attenuation weighting on the temperature difference, and constructing a temperature gradient matrix; calculating the cosine of the included angle in combination with the direction of the main airflow velocity vector, accumulating the nonlinear term of the thermal driving gradient and the cosine weight of the included angle, and obtaining a comprehensive thermal coupling degree; using a nonlinear temperature gain term to calculate a thermal trend operator based on the temperature rise slope and the logarithmic fluctuation characteristics in the neighborhood window; and using the comprehensive thermal coupling degree and the thermal trend operator to correct the original power value to obtain a final execution power, so as to control the heating element. By comprehensively considering the spatial heat backflow and the time heat storage trend, the present application reduces the interference of airflow anisotropy and thermal inertia on temperature control, and improves the stability of metallurgical node temperature regulation.

Description

Technical Field

[0001] This invention relates to the field of temperature control technology, and in particular to a method and system for adaptive temperature control of metallurgical nodes based on multi-zone coordinated regulation. Background Technology

[0002] In the industrial preheating and sintering process of metallurgical materials, ensuring the uniformity of the temperature field at metallurgical nodes is fundamental to improving the structural performance of materials and the success rate of heat diversion. Real-time monitoring and energy distribution of multiple heating zones within the furnace body are used to cope with thermal disturbances caused by the complex and ever-changing production environment. This plays a crucial role in ensuring the consistency of the physicochemical properties of metallurgical materials and improving the safety of equipment operation.

[0003] In existing technologies, temperature regulation is typically achieved using a multi-channel controller in conjunction with a linear decoupling matrix. This method takes into account the mutual influence between temperature zones, compensating for the control output of adjacent temperature zones through a pre-set linear gain coefficient. Simultaneously, it utilizes a low-pass filter operator to process the acquired temperature signal, reducing the interference of environmental electromagnetic noise on the feedback loop, thereby improving the regulation effect of multi-temperature zone coordinated control.

[0004] However, existing technologies have limitations when dealing with temperature control tasks in metallurgical furnaces with strong thermal radiation characteristics. Because the heat transfer process inside a metallurgical furnace involves high-temperature radiation and nonlinear convection, the intensity of heat diffusion across different temperature ranges does not change linearly with temperature difference, but rather exhibits a nonlinear increasing trend as temperature rises. Linear compensation logic, based on a fixed mapping ratio, struggles to accurately capture the nonlinear heat backflow phenomenon caused by temperature fluctuations under high-temperature conditions. This deviation in thermal coupling strength leads to energy imbalances in different regions during the heating phase, increasing the probability of dynamic temperature oscillations at various metallurgical nodes, and consequently affecting the accuracy of stable control of complex temperature fields. Summary of the Invention

[0005] To address the aforementioned technical problems of imbalanced energy distribution in various temperature zones and dynamic temperature oscillations at metallurgical nodes caused by the inability of linear compensation logic to adapt to nonlinear thermal coupling and heat storage delay within metallurgical furnaces, this invention provides solutions in the following aspects.

[0006] In a first aspect, the present invention provides an adaptive control method for metallurgical node temperature based on multi-zone coordinated regulation, the method comprising the steps of: The temperatures of multiple temperature zones within the metallurgical furnace are collected at various times according to a preset sampling period. Taking any given time as the current time, the temperature difference between any two temperature zones and the straight-line distance between their physical centers are calculated. The temperature difference is then nonlinearly weighted with spatial attenuation using the straight-line distance to obtain the thermal driving gradient between any two temperature zones at the current time, thus constructing a temperature gradient matrix. The direction of the main airflow velocity vector within the metallurgical furnace is obtained. Taking any temperature zone as the reference temperature zone, the cosine of the angle between the line connecting the reference temperature zone and the other temperature zones and the direction of the main airflow velocity vector is calculated. The sign function value of the thermal driving gradient between the reference temperature zone and the other temperature zones, and the nonlinearity of the magnitude of the thermal driving gradient are then calculated. The linear power term and the weight based on the cosine of the included angle are multiplied and accumulated to obtain the comprehensive thermal coupling degree of the reference temperature zone. For any temperature zone, the temperature rise slope is calculated based on the temperature at the current time and the previous time. The logarithmic fluctuation feature of the rate of change of the temperature rise slope in the neighborhood window corresponding to the current time is extracted. The nonlinear temperature gain term formed by the temperature at the current time is used to nonlinearly weight the temperature rise slope and the logarithmic fluctuation feature to obtain the thermal trend operator of the corresponding temperature zone. The preset original power value is corrected using the comprehensive thermal coupling degree and the thermal trend operator to obtain the final execution power. The heating element corresponding to the temperature zone is controlled based on the final execution power.

[0007] This invention acquires the temperatures of multiple temperature zones within a metallurgical furnace and calculates the spatially attenuated weighted thermal driving gradient to determine the nonlinear thermal coupling strength between temperature zones. By combining the direction of the main airflow velocity vector with cosine modulation of the thermal driving gradient, the transport effect of flow field anisotropy on thermal energy migration is evaluated, obtaining a comprehensive thermal coupling degree reflecting the real-time disturbance intensity of each temperature zone. The invention utilizes the cubic term of temperature to characterize the physical law of the dramatic increase in thermal radiation intensity with increasing absolute temperature under high-temperature conditions, and extracts the accelerated process of material heat storage evolution using logarithmic fluctuation characteristics, obtaining a thermal trend operator reflecting the dynamic trend of energy accumulation. This invention uses the comprehensive thermal coupling degree and the thermal trend operator to correct the original power value, ensuring the continuity of the control output at the physical level under multi-factor interference, reducing the probability of energy distribution imbalance in each temperature zone within the metallurgical furnace and dynamic oscillations in the temperature of metallurgical nodes, and improving the accuracy of stable regulation of complex temperature fields.

[0008] Preferably, the step of using the straight-line distance to perform a nonlinear weighting of the temperature difference to achieve spatial attenuation, and obtaining the thermal driving gradient of the pairwise temperature intervals at the current moment, includes: processing the ratio of the straight-line distance to the preset thermal influence radius using a negative natural exponential function to obtain the distance attenuation amount; and multiplying the temperature difference by the distance attenuation amount to obtain the thermal driving gradient of the pairwise temperature intervals at the current moment.

[0009] This invention utilizes a negative natural exponential function to process the ratio of straight-line distance to a preset thermal influence radius to obtain the distance attenuation. This characterizes the exponential loss of heat as distance increases during heat conduction in refractory materials and non-homogeneous media. By multiplying the temperature difference by the distance attenuation, the thermal driving gradient is obtained. In the calculation process, the redundant contribution of the temperature zone at a distance is suppressed, and the interference of weak disturbances in irrelevant areas at a distance on the accuracy of signal recognition is reduced. This ensures that the system can quickly lock onto the core heat source with the highest influence weight on the controlled node.

[0010] Preferably, obtaining the direction of the main airflow velocity vector within the metallurgical furnace includes: acquiring multiple local wind speed vectors through a wind speed sensing array distributed within the furnace chamber of the metallurgical furnace; and performing vector synthesis calculations on the multiple local wind speed vectors to obtain the direction of the main airflow velocity vector.

[0011] Preferably, the step of using any temperature zone as a reference temperature zone and calculating the cosine of the angle between the line connecting the reference temperature zone and the other temperature zones with respect to the direction of the main airflow velocity vector includes: using any temperature zone as a reference temperature zone, obtaining the physical center coordinates of the reference temperature zone and the other temperature zones; calculating the relative displacement based on the physical center coordinates of the reference temperature zone and the other temperature zones to determine the geometric position vector pointing from the reference temperature zone to the other temperature zones; and performing a vector dot product operation on the geometric position vector and the direction of the main airflow velocity vector to obtain the cosine of the angle.

[0012] Preferably, the overall thermal coupling degree satisfies the following relationship: ; in, It is the first The temperature zone is in The overall thermal coupling degree at any given moment; It is a symbolic function; It is the first In the temperature gradient matrix at time t, the th Temperature zone for the first The thermal driving gradient of the temperature zone; yes The modulus length; It is the coupling dimension conversion factor; It is the weighting coefficient for the influence of airflow; It is the first Temperature zone points to the first The angle between the vector connecting the temperature zones and the direction of the main airflow velocity vector; This is the total number of sensors.

[0013] This invention obtains the comprehensive thermal coupling degree by multiplying and accumulating the nonlinear power term of the sign function and the modulus and the weight based on the cosine of the included angle. It uses the nonlinear power term to simulate the enhanced law of mixed radiation and convection heat transfer under high temperature environment, improves the system’s compensation sensitivity to strong thermal coupling interference, and realizes the accurate characterization of the multi-source nonlinear superposition process of total heat load in metallurgical furnace.

[0014] Preferably, the step of calculating the temperature rise slope based on the temperature at the current moment and the temperature at the previous moment includes: performing a difference operation on the temperature at the current moment and the temperature at the previous moment to obtain a temperature difference value; and recording the ratio of the temperature difference value to a preset sampling period as the temperature rise slope.

[0015] Preferably, the heat trend operator satisfies the following relation: ; in, It is the first The temperature zone is in The heat trend operator in real time; , They are the first The temperature zone is in , Temperature at any moment; It is the preset sampling period; It is the reference temperature; It is the first The temperature zone is in The rate of change of the temperature rise slope within the neighborhood window at any given time; It is the natural logarithm function; It is the absolute value symbol.

[0016] This invention combines the cubic power of temperature and the natural logarithm of the rate of change of the temperature rise slope to obtain a thermal trend operator. The cubic power term is used to increase the sensitivity of trend identification in the high-temperature zone, and the slope fluctuation is processed by logarithmic compression to filter out numerical glitches caused by random sensor jitter or local airflow disturbances, thereby improving the robustness of the assessment of the material's thermal storage evolution trend.

[0017] Preferably, obtaining the rate of change of the temperature rise slope within the neighborhood window includes: performing a difference operation on the temperature rise slope of each adjacent time within the neighborhood window corresponding to the current time to obtain a temperature rise slope change sequence of the neighborhood window; and performing a statistical average operation on the temperature rise slope change sequence to obtain the rate of change of the temperature rise slope within the neighborhood window.

[0018] Preferably, the final execution power satisfies the following relationship: ; in, It is the first The temperature zone is in The final execution power at any given moment; It is the first The temperature zone is in Original power value at any given time; , They are the first The temperature zone is in The overall thermal coupling degree and thermal trend operator at any given time; It is the sensitivity conversion coefficient; It is the hyperbolic tangent function; , These are functions for finding the maximum and minimum values, respectively.

[0019] This invention utilizes the hyperbolic tangent function to perform reverse suppression processing on the comprehensive thermal coupling degree and thermal trend operator, and makes the final execution power approach the original power value when the interference approaches zero, ensuring the linear transition of the adjustment process in the critical region. Combined with the limiting structure, the power modulation ratio is limited to the preset range, reducing the frequency of power output jumping beyond the physical limit, and ensuring the dynamic adaptive balance of the temperature field at the metallurgical node.

[0020] In a second aspect, the present invention provides a metallurgical node temperature adaptive control system based on multi-zone coordinated regulation. The metallurgical node temperature adaptive control system based on multi-zone coordinated regulation includes a memory and a processor. The memory stores computer program instructions, and when the computer program instructions are executed by the processor, the metallurgical node temperature adaptive control method based on multi-zone coordinated regulation of the first aspect of the present invention is implemented.

[0021] By adopting the above technical solution, the metallurgical node temperature adaptive control method based on multi-zone coordinated adjustment of the first aspect of the present invention is generated into a computer program and stored in a memory so that it can be loaded and executed by a processor. In this way, a terminal device can be made based on the memory and the processor for convenient use.

[0022] The beneficial effects of this invention are as follows: This invention addresses the temperature control oscillation problem caused by nonlinear thermal coupling and material heat storage delay within metallurgical furnaces. It utilizes a spatially attenuated weighted thermal driving gradient and the main airflow velocity vector direction for collaborative correction, obtaining a comprehensive thermal coupling degree reflecting the intensity of thermal backflow in the spatial dimension. This reduces the impact of flow field anisotropy on the identification of interference in multi-temperature zone collaborative control. By extracting a thermal trend operator based on the cubic gain of absolute temperature and the rate of change of temperature rise slope, the energy accumulation intensity of the material in different temperature ranges is evaluated, addressing the time-dimensional lag effect caused by material thermal inertia. Furthermore, the saturation characteristics of the hyperbolic tangent function are used to smooth the power, ensuring the final execution power maintains physical continuity under multi-factor interference, reducing the probability of frequent oscillations in the actuator, and improving the steady-state control accuracy and equipment operation safety of the metallurgical furnace under complex thermal conditions. Attached Figure Description

[0023] Figure 1 A flowchart of a metallurgical node temperature adaptive control method based on multi-zone coordinated regulation provided in an embodiment of the present invention; Figure 2 A comparative schematic diagram of temperature changes over time in a temperature zone provided in an embodiment of the present invention; Figure 3 A comparative schematic diagram showing the change in temperature control deviation with overall thermal coupling degree provided in an embodiment of the present invention; Figure 4 The structural block diagram of the metallurgical node temperature adaptive control system based on multi-zone coordinated regulation provided in the embodiments of the present invention is shown. Detailed Implementation

[0024] The first aspect of this invention provides an adaptive temperature control method for metallurgical nodes based on multi-zone coordinated regulation, such as... Figure 1 As shown, the method includes steps S100-S500: Step S100: Collect the temperature of multiple temperature zones inside the metallurgical furnace at various times according to the preset sampling period.

[0025] It should be noted that, considering the spatially discrete nature of the temperature distribution inside a metallurgical furnace, a global reference benchmark reflecting the furnace's thermal energy distribution structure can be obtained by establishing a unified time base and spatial vector relationship. Since the thermodynamic responses of different temperature zones differ along the time axis, synchronous sampling is performed to ensure spatiotemporal consistency in the logical calculations of subsequent adjustment commands.

[0026] Specifically, by using several temperature sensors installed at different locations within the metallurgical furnace body, the interior of the furnace is spatially divided into multiple corresponding temperature zones. During execution, the raw temperature corresponding to each sensor node is synchronously collected at a unified timestamp according to a preset sampling period, serving as the temperature feedback for each temperature zone at each moment.

[0027] Considering the thermal inertia of the metallurgical furnace, the temperature field changes have a significant time lag. If the sampling period is set too short, such as less than 0.5s, high-frequency random thermal noise can be easily introduced, causing false fluctuations in the calculated temperature rise slope, which in turn induces frequent oscillations in the actuator. If the sampling period is set too long, such as greater than 10s, it will be impossible to capture the sudden changes in thermal backflow caused by the anisotropy of the flow field in time. In this embodiment, the sampling period is preferably 1s to 5s.

[0028] Thus, the temperatures of each temperature zone at each time point and the straight-line distance between the physical centers of each temperature zone were obtained.

[0029] Step S200: Taking any time as the current time, calculate the temperature difference between any two temperature intervals and the straight-line distance between their physical centers. Use the straight-line distance to perform a nonlinear weighting of the temperature difference with spatial attenuation to obtain the thermal driving gradient between any two temperature intervals at the current time to construct a temperature gradient matrix.

[0030] It should be noted that, considering the varying physical orientations and heating loads of the heated nodes within a metallurgical furnace, heat conduction in refractory materials and non-homogeneous media exhibits significant distance loss characteristics, resulting in a non-steady-state local topological distribution rather than isotropic energy field within the furnace. If disturbance signals across all temperature zones are indiscriminately superimposed as scalars, weak redundant disturbances in distant, unrelated regions can easily be introduced, thereby reducing the signal recognition accuracy of the control logic. Such redundant disturbances are typically ignored or filtered as measurement noise in traditional PID control. However, this invention reveals that in strongly coupled multi-temperature-zone scenarios, the cumulative effect of these disturbances can induce systematic control deviations. Based on this, this invention establishes a spatial topological mapping through a matrix structure, characterizing the evolution of heat energy radiation and diffusion along different paths in discrete space at the logical level. It utilizes physical distance to nonlinearly weight the temperature difference driving force, thereby ensuring that the control system can quickly identify the core heat source with the highest impact weight on the controlled node from multi-dimensional disturbances, providing structured data support for subsequent asymmetric energy compensation.

[0031] Specifically, based on the spatial installation coordinates of each sensor, the relative positions and straight-line distances between the physical centers of each temperature zone are determined. For multiple temperature zones at each time point, a temperature gradient matrix reflecting the global spatial topology is constructed based on the temperature difference between any two temperature zones and the straight-line distance between the physical centers of any two temperature zones. Any element in the temperature gradient matrix satisfies the following relationship: ; in, It is the first In the temperature gradient matrix at time t, the th Temperature zone for the first The thermal driving gradient of the temperature zone; , They are the first , The temperature zone is in Temperature at any moment; It is the first Temperature zone and the first The straight-line distance between the physical centers of the temperature zone; It is the heat-affected radius, measured in length; It is a natural exponential function.

[0032] In this relation, The original driving force and intensity of heat conduction, when the first Temperature in the temperature zone is lower than the first In the temperature range, a negative temperature difference indicates that heat is transferred from the first... Temperature zone flow direction Temperature zone; This component is used to determine the exponential heat loss that occurs during heat conduction in a medium with increasing distance. As the linear distance between temperature zones increases, this component decreases non-linearly and tends to... This is used to suppress the contribution of distant temperature regions in the calculation results.

[0033] It should be noted that the heat-affected radius... The value needs to be set based on the thermal conductivity and physical thickness of the furnace lining material: for furnace linings using lightweight ceramic fibers or other high thermal resistance materials and with a small thermal diffusion range, the value can be reduced. The preset value, such as setting for m, to improve the sensitivity of identifying near-range coupling interference; for working conditions using dense clay bricks or silicon carbide and other high thermal conductivity materials with a wide thermally affected range, the value can be appropriately increased. The preset value, such as setting for m, to enhance the ability to capture wide-area heat backflow phenomena. In this embodiment, Preferred m.

[0034] Thus, the temperature gradient matrix at each time point was obtained.

[0035] Step S300: Obtain the direction of the main airflow velocity vector in the metallurgical furnace. Taking any temperature zone as the reference temperature zone, calculate the cosine of the angle between the line connecting the reference temperature zone and the other temperature zones with respect to the direction of the main airflow velocity vector. Multiply and sum the sign function value of the thermal driving gradient between the reference temperature zone and the other temperature zones, the nonlinear power term of the magnitude of the thermal driving gradient, and the weight based on the cosine of the angle to obtain the comprehensive thermal coupling degree of the reference temperature zone.

[0036] It should be noted that in the case of forced airflow circulation within a metallurgical furnace, heat transfer efficiency is guided by the direction of the airflow velocity vector, resulting in a physical asymmetry in heat exchange between the downstream and upstream directions. If a conventional symmetrical coupling model is used, the transport effect of the airflow on heat migration will be ignored, increasing the probability of overcompensation in the downwind temperature zone. Therefore, this invention introduces an anisotropic correction coefficient to vectorize and weight the thermally driven gradient, and uses a cosine projection operator to characterize the modulation effect of the airflow velocity field on the thermal coupling strength, thereby reconstructing the distribution law of asymmetric heat flux in logical calculations.

[0037] First, the cosine of the angle between the line connecting the temperature zones and the direction of the main airflow velocity vector is calculated. It should be noted that, considering the intensity of heat transport by airflow depends on the geometric consistency between the heat flow path and the airflow vector, obtaining the cosine of the angle between the line connecting the temperature zones and the direction of the main airflow velocity vector is to establish a spatial weighting benchmark for thermal convection gain. Since temperature zones in different orientations have different projection components relative to the airflow center, the vector angle cosine can accurately locate which temperature zones are on the airflow carrying path, thus providing a geometric basis for subsequent asymmetric physical compensation.

[0038] Specifically, the velocity vector direction of the main airflow within the metallurgical furnace is first obtained. Considering the distribution characteristics of the flow field within the metallurgical furnace, the velocity vector direction of the main airflow is acquired through a wind speed sensing array distributed within the furnace chamber. As a preferred embodiment, multiple local wind speed vectors are acquired through the wind speed sensing array distributed within the furnace chamber; then, vector synthesis operations are performed on these multiple local wind speed vectors, such as weighted average or summation, to obtain the synthesized vector direction. This synthesized vector direction is used as the velocity vector direction of the main airflow, thereby transforming the discrete flow field characteristics at multiple points into a global dominant vector. Then, taking any temperature zone as a reference temperature zone, the cosine of the angle between the line connecting this reference temperature zone and the other temperature zones and the velocity vector direction of the main airflow is calculated. Specifically, the physical center coordinates of the reference temperature zone and the other temperature zones are obtained; the relative displacement is calculated based on the physical center coordinates of the reference temperature zone and the other temperature zones to determine the geometric position vector pointing from the reference temperature zone to the other temperature zones; finally, the geometric position vector and the direction of the main airflow velocity vector are subjected to a vector dot product operation to directly calculate the cosine of the angle reflecting the spatial geometric relationship between the two, which is used to evaluate the thermal convection intensity weight in different spatial directions.

[0039] Then, the overall thermal coupling degree, reflecting the comprehensive disturbance level of each temperature zone, is calculated. It should be noted that the total heat load within the metallurgical furnace is the result of the nonlinear superposition of multiple heat sources, and the heat transfer rate under high-temperature conditions exhibits a nonlinear evolution characteristic with temperature difference. Calculating the overall coupling degree, which includes power-law correction and vector gain, is to transform the discrete spatial gradient into a trend operator that can be directly used for control law modulation. Introducing a dimension conversion factor is to harmonize the thermally driven gradient in the spatial dimension with the temperature variation rate in the time dimension, thereby ensuring that the energy compensation amount and the actuator's action amplitude are accurately aligned in physical and logical terms.

[0040] Based on the above logic, for any given moment, the overall thermal coupling degree of any temperature region satisfies the following relationship: ; in, It is the first The temperature zone is in The overall thermal coupling degree at any given moment; It is a sign function used to determine the direction of heat flow. When the time is positive, the output is 1. Output -1 when it is negative. Output 0 when the value is zero, used to maintain the algebraic properties of heat inflow and outflow during the summation process; It is the first In the temperature gradient matrix at time t, the th Temperature zone for the first The thermal driving gradient of the temperature zone; yes The modulus length; It is the coupling dimension conversion factor, with units of ; It is the weighting coefficient for the influence of airflow; It is the first Temperature zone points to the first The angle between the vector connecting the temperature zones and the direction of the main airflow velocity vector; This is the total number of sensors.

[0041] In this relation, The vector direction used to determine thermal energy migration ensures that the heat can be identified in the first... Temperature accumulation or towards the first Temperature loss; This is to characterize the nonlinear enhancement law of mixed radiation and convection heat transfer under high temperature environment. This nonlinear amplification improves the compensation sensitivity of the present invention to strong interference. The value of this power term is not limited to 1.5, and can be adjusted between 1.2 and 1.8 according to the actual heat transfer characteristics. The coupling strength is directionally increased based on the airflow direction, when the first... Temperature zone is in the first When the temperature zone is upwind, the cosine projection value is positive. Increasing this value characterizes the physical fact that airflow accelerates heat transfer, thus transforming the discrete spatial thermal driving force into a reflection of the first... The comprehensive thermal coupling degree of the real-time disturbance intensity in the temperature range.

[0042] It should be noted that the weighting coefficient for airflow influence and dimensional conversion factor The value needs to be set according to the power specifications of the furnace circulating fan and the furnace volume: for scenarios equipped with high-power circulating fans, the value needs to be increased. And decrease If we assume It is 0.5. The value is 0.08; for scenarios with relatively stable flow fields, it needs to be appropriately reduced. And increase If we assume It is 0.2. It is 0.12. In this embodiment, The preferred value is 0.35. The preferred value is 0.1.

[0043] Thus, the comprehensive thermal coupling degree of each temperature zone at each time point was obtained.

[0044] Step S400: For any temperature zone, calculate the temperature rise slope based on the temperature at the current time and the previous time, extract the logarithmic fluctuation feature of the rate of change of the temperature rise slope in the neighborhood window corresponding to the current time, and use the nonlinear temperature gain term formed by the temperature at the current time to perform nonlinear weighting on the temperature rise slope and the logarithmic fluctuation feature to obtain the thermal trend operator of the corresponding temperature zone.

[0045] It should be noted that during the high-temperature heating process of metallurgical materials, the furnace lining and the heated material exhibit significant heat storage characteristics, resulting in a time lag between the real-time signal fed back by the temperature sensor and the actual applied control energy. Adjusting solely based on the current instantaneous temperature deviation is insufficient to address the continuous heating or cooling trends caused by the material's thermal inertia. Therefore, this invention introduces a dynamic trend assessment mechanism reflecting the rate of energy change. By constructing a time evolution operator incorporating an absolute temperature nonlinear gain, it is used to assess the energy accumulation intensity of the material in different temperature ranges.

[0046] First, for any temperature zone, the temperature rise slope is calculated based on the current temperature and the previous temperature. It should be noted that, considering the temperature change trajectory over time directly reflects the energy supply and demand balance within the furnace, calculating the temperature difference between adjacent sampling periods is to obtain the instantaneous thermal momentum of the heated nodes. Since environmental noise in the metallurgical process may interfere with the accuracy of a single sample, slope normalization using the sampling period is used to establish a logical benchmark that can objectively characterize the rate of temperature rise.

[0047] Specifically, taking any given moment as the current moment, the temperature of each temperature zone at the current moment and the temperature at the previous moment are obtained; the temperature is then used to perform a differential operation with a preset sampling period to obtain the temperature rise slope at the corresponding moment.

[0048] Then, a neighborhood window is constructed for each time step, and the rate of change of the temperature rise slope within the neighborhood window is obtained. It should be noted that, to accurately assess the abruptness and sustained trend of temperature changes at the current sampling moment, it is not sufficient to rely solely on the instantaneous temperature difference; it is also necessary to rely on the historical sampling sequence within a certain range prior to that moment to statistically determine the baseline level and fluctuation intensity of local thermal energy changes. Constructing a time-dimensional neighborhood window is to filter out random thermal noise at the logical level and capture the energy accumulation characteristics with physical inertia.

[0049] Specifically, a neighborhood window is constructed with the current sampling time as the endpoint. As a preferred implementation, the size of the neighborhood window is preferably set to 10 sampling periods. This window size is positively correlated with the characteristic time constant of heat propagation within the metallurgical furnace. Considering the electrothermal response of heating elements in conventional metallurgical furnaces and their conduction delay in refractory materials, at conventional sampling frequencies, the significant evolution of thermal momentum typically occurs within a span of 3 to 5 sampling periods. A neighborhood of 10 sampling periods can fully encompass the entire thermal response process caused by a single power adjustment, ensuring the accuracy of the local temperature rise slope statistics, while avoiding interference from excessively large windows due to distant historical operating conditions or irrelevant process switching. Implementers can adjust the window size according to the actual furnace volume, heating power reserve, and the thermal inertia of the heated material: for large refining furnaces with huge volumes and slow thermal responses, the window can be appropriately increased to 20 to 30 sampling periods; for fast-response scenarios such as small samples and high-frequency induction heating, the window can be appropriately reduced to 5 to 7 sampling periods.

[0050] After constructing the neighborhood window, the temperature rise slope of each adjacent time point is differentially calculated within the neighborhood window corresponding to each time point to obtain the temperature rise slope change sequence of the neighborhood window; the temperature rise slope change sequence is statistically averaged to obtain the temperature rise slope change rate within the neighborhood window.

[0051] Finally, a thermal trend operator reflecting the material's heat storage evolution trend is calculated. It should be noted that the specific heat capacity and thermal radiation capacity of metallurgical materials increase non-linearly with increasing temperature, resulting in a much higher energy accumulation rate at high temperatures compared to low temperatures. Introducing a cubic term of absolute temperature to weight the temperature rise rate is intended to characterize the accelerated process of heat storage evolution under high-temperature conditions, thus providing feedforward features with physical a priori support for predicting potential temperature overshoot.

[0052] Based on the above logic, the heat trend operator satisfies the following relation: ; in, It is the first The temperature zone is in The time-in-time thermal trend operator is used to characterize the dynamic trend of energy accumulation; , They are the first The temperature zone is in , Temperature at any moment; It is the preset sampling period; It is a reference temperature for material properties, used to establish an evaluation benchmark for thermal energy excitation; It is the first The temperature zone is in The rate of change of the temperature rise slope within the neighborhood window at time t is calculated by taking its dimensionless pure value when performing logarithmic operations. It is the natural logarithm function; It is the absolute value symbol.

[0053] In this relation, It is used to characterize the basic heat exchange rate at the current moment, reflecting the instantaneous derivative of temperature fluctuations. By using the cube of temperature as a nonlinear gain term, the physical law of the dramatic increase in thermal radiation intensity with increasing absolute temperature under high-temperature conditions is characterized, enabling the present invention to have higher trend recognition sensitivity in high-heat regions. By examining the first Logarithmically compressing the slope fluctuations within the neighborhood window corresponding to a given time allows for the extraction of thermal inertial fluctuation intensity while filtering out numerical glitches caused by random sensor jitter or local airflow disturbances.

[0054] It should be noted that the reference temperature The value needs to be set according to the rated process temperature of the metallurgical material: for scenarios using low-melting-point alloys or low-temperature baking processes, the value can be reduced. The preset value, such as setting for To improve the sensitivity of this invention to changes in heat storage in low-temperature regions; for scenarios using high-alumina refractory materials or high-temperature refining conditions, the sensitivity can be appropriately increased. The preset value, such as setting for This ensures that the value of the thermal trend operator remains within a stable growth range under ultra-high temperature conditions. In this embodiment, Preferred .

[0055] At this point, the thermal trend operators for each temperature zone have been obtained.

[0056] Step S500: Correct the preset original power value using the integrated thermal coupling degree and thermal trend operator to obtain the final execution power, and control the heating element corresponding to the temperature zone based on the final execution power.

[0057] It should be noted that, in order to achieve unified compensation for spatial thermal interference and temporal heat storage trends, the calculated comprehensive thermal coupling degree and thermal trend operator need to be logically integrated. Directly adding the interference term to the basic control quantity can easily cause abrupt changes in control gain, leading to frequent jitter in the actuator and thus shortening the lifespan of the heating element. Therefore, this invention introduces a composite sensitivity coefficient with dimensional conversion capabilities, combined with a nonlinear mapping function with saturation characteristics, to smoothly correct the basic control power. This ensures that the control output maintains physical continuity under multi-factor interference, thereby restoring the true temperature control target requirements in complex thermal coupling environments.

[0058] Specifically, the final execution power is calculated based on the linear summation results of the comprehensive thermal coupling degree and the thermal trend operator. It should be noted that, considering that both the thermal backflow interference in the spatial dimension and the thermal inertia in the temporal dimension are physically manifested as a disruption of the thermal balance of the controlled node, the composite interference intensity is constructed by algebraically summing the two, and the hyperbolic tangent function is introduced as a nonlinear mapping operator. This is to provide a natural saturation buffer when the interference fluctuates violently, and to avoid power output jumps that exceed physical limits.

[0059] Based on the above logic, for any temperature range, its final execution power at any given moment satisfies the following relationship: ; in, It is the first The temperature zone is in The final execution power at any given moment; It is the first The temperature zone is in The original power value at any given time is usually calculated by a conventional PID controller or feedforward controller based on the deviation between the set temperature and the feedback temperature; , They are the first The temperature zone is in The overall thermal coupling degree and thermal trend operator at any given time; It is the sensitivity conversion factor, with units of ; It is the hyperbolic tangent function; , These are functions for finding the maximum and minimum values, respectively.

[0060] In this relation, Achieved the first The temperature zone is in A comprehensive evaluation of spatial interference and temporal trends at any given moment; this relationship is expressed through... The structure establishes reverse suppression logic; when enhanced external heat backflow is detected, it means... Increase or accelerate the internal heat storage rate As the interference increases, the modulation coefficient automatically decreases, representing the physical process of reducing power output to suppress temperature rise overshoot. Compared to traditional subtractive compensation, this reverse suppression structure has a gradual characteristic where the modulation coefficient approaches 1 when the interference approaches zero, ensuring a linear transition of the adjustment process in the critical region and avoiding power jumps under interference-free conditions. , The constructed limiting structure ensures that the power modulation ratio is within the range of 100%. The range is used to avoid power overflow or negative power during the adjustment process, ensuring the dynamic adaptive balance of the temperature field at the metallurgical node.

[0061] It should be noted that the sensitivity conversion coefficient The value of needs to be set based on the standard deviation of the thermal inertia of the controlled system and the allowable process deviation: for precision metallurgical conditions with high thermal sensitivity and narrow allowable fluctuation range, it can be increased. The preset value, such as setting The value is set to 0.4 to enhance the sensitivity to suppressing small temperature rise trends; for large-volume, thermally inertial, and insensitive to power changes in extensive heating scenarios, the value can be appropriately reduced. The preset value, such as setting The value is set to 0.2 to maintain the system's regulating power. In this embodiment, The preferred value is 0.3.

[0062] The final execution power of each temperature zone at each time is converted into a standard analog signal or pulse width modulation signal by the controller and sent to the thyristor voltage regulator or solid-state relay connected to the heating element of the corresponding temperature zone. By adjusting the output effective voltage or conduction duty cycle of the heating power supply in real time, the energy compensation is applied to the heating physical field inside the furnace without loss, thereby physically offsetting the temperature fluctuations caused by spatial heat backflow and material heat storage inertia, and ensuring the constant temperature adaptive balance of the multi-temperature zone field inside the metallurgical furnace.

[0063] like Figure 2 As shown in the figure, this is a comparative diagram of temperature changes over time in different temperature zones. The horizontal axis represents time, and the vertical axis represents the temperature of any given zone. The dashed line corresponds to traditional PID control, while the solid line corresponds to the method of this invention. From the figure, we can see that during the initial heating and transition phases, the dashed line exhibits several large-amplitude oscillations, demonstrating significant overshoot characteristics. This reflects the lag in the adjustment capability of the traditional PID control method when facing heat storage inertia. The solid line converges rapidly as it approaches the target value, and the entire transition process is smooth without significant fluctuations. This indicates that this invention, through advance prediction using a thermal trend operator, can effectively counteract heat storage inertia and improve the stability of the thermal field regulation.

[0064] For example Figure 3As shown in the figure, this is a comparative schematic diagram of temperature control deviation changing with the overall thermal coupling degree. To visually demonstrate the anti-interference performance of this invention under different loads, the overall thermal coupling degree in the figure has been normalized. The horizontal axis represents the normalized value of the overall thermal coupling degree, and the vertical axis represents the temperature control deviation, that is, the absolute deviation of the temperature of any temperature zone from the preset target temperature at each sampling time. The dashed line with dots corresponds to traditional PID control, and the solid line with diamonds corresponds to the method of this invention. From the deviation trend under different operating conditions: in the range of low overall thermal coupling degree, the deviation of both methods remains at a low level. As the overall thermal coupling degree increases, the dashed line shows a sharp upward trend, indicating that the anti-interference ability of the traditional PID control method is severely degraded under strong interference environment. The solid line only shows a gradual trend with the increase of interference intensity, and the value remains in the low range. This shows that this invention, by vectorizing and weighting the spatial thermal coupling, can still maintain accurate temperature control effect under harsh interference conditions and has stronger adaptive capability.

[0065] Thus, the adaptive temperature control of metallurgical nodes based on multi-zone coordinated regulation was completed.

[0066] The second aspect of this embodiment provides an adaptive temperature control system for metallurgical nodes based on multi-zone coordinated regulation, such as... Figure 4 As shown, the metallurgical node temperature adaptive control system based on multi-zone coordinated regulation includes a memory and a processor. The memory stores computer program instructions, and when the computer program instructions are executed by the processor, the metallurgical node temperature adaptive control method based on multi-zone coordinated regulation of the first aspect of the present invention is implemented.

[0067] The metallurgical node temperature adaptive control system based on multi-zone coordinated regulation also includes other components well known to those skilled in the art, such as communication buses and communication interfaces. Their settings and functions are known in the art and will not be described in detail here.

[0068] In this invention, the aforementioned memory can be any tangible medium containing or storing a program that can be used or combined with an instruction execution system, apparatus, or device. For example, a computer-readable storage medium can be any suitable magnetic or magneto-optical storage medium, such as resistive random access memory (DRAM), dynamic random access memory (DRAM), static random access memory (SRAM), enhanced dynamic random access memory (DRAM), high-bandwidth memory, hybrid memory cube, etc., or any other medium that can be used to store desired information and can be accessed by an application, module, or both. Any such computer storage medium can be part of a device or accessible to or connected to a device.

[0069] The above are all preferred embodiments of the present invention and are not intended to limit the scope of protection of the present invention. Therefore, all equivalent changes made in accordance with the structure, shape and principle of the present invention should be covered within the scope of protection of the present invention.

Claims

1. A metallurgical node temperature adaptive control method based on multi-zone coordinated regulation, characterized in that, include: The temperature of multiple temperature zones inside the metallurgical furnace at various times is collected according to the preset sampling period. Taking any moment as the current moment, calculate the temperature difference between any two temperature intervals and the straight-line distance between their physical centers. Use the straight-line distance to perform a nonlinear weighting of the temperature difference with spatial decay to obtain the thermal driving gradient between any two temperature intervals at the current moment and construct the temperature gradient matrix. Obtain the direction of the main airflow velocity vector in the metallurgical furnace. Take any temperature zone as the reference temperature zone and calculate the cosine of the angle between the line connecting the reference temperature zone and the other temperature zones with respect to the direction of the main airflow velocity vector. Multiply and sum the sign function value of the thermal driving gradient between the reference temperature zone and the other temperature zones, the nonlinear power term of the modulus of the thermal driving gradient, and the weight based on the cosine of the angle to obtain the comprehensive thermal coupling degree of the reference temperature zone. For any temperature zone, the temperature rise slope is calculated based on the temperature at the current time and the previous time. The logarithmic fluctuation feature of the rate of change of the temperature rise slope in the neighborhood window corresponding to the current time is extracted. The nonlinear temperature gain term formed by the temperature at the current time is used to nonlinearly weight the temperature rise slope and the logarithmic fluctuation feature to obtain the thermal trend operator of the corresponding temperature zone. The preset original power value is corrected by the comprehensive thermal coupling degree and thermal trend operator to obtain the final execution power, and the heating element corresponding to the temperature zone is controlled based on the final execution power.

2. The adaptive temperature control method for metallurgical nodes based on multi-zone coordinated regulation according to claim 1, characterized in that, The step of using the linear distance to spatially attenuate the temperature difference through nonlinear weighting to obtain the thermal driving gradient of each pair of temperature intervals at the current moment includes: The distance attenuation is obtained by processing the ratio of the straight-line distance to the preset thermally affected radius using a negative natural exponential function. Multiplying the temperature difference by the distance attenuation, we obtain the thermal driving gradient between the two temperature intervals at the current moment.

3. The adaptive temperature control method for metallurgical nodes based on multi-zone coordinated regulation according to claim 1, characterized in that, The acquisition of the main gas flow velocity vector direction within the metallurgical furnace includes: Multiple local wind speed vectors are collected by a wind speed sensing array distributed within the furnace chamber of the metallurgical furnace. The direction of the main airflow velocity vector is obtained by performing vector synthesis on the multiple local wind speed vectors.

4. The adaptive temperature control method for metallurgical nodes based on multi-zone coordinated regulation according to claim 1, characterized in that, The step of calculating the cosine of the angle between the line connecting any temperature zone to the other temperature zones and the direction of the main airflow velocity vector includes: Using any temperature zone as a reference temperature zone, obtain the physical center coordinates of the reference temperature zone and the other temperature zones; The relative displacement between the reference temperature zone and the physical center coordinates of the other temperature zones is calculated to determine the geometric position vector from the reference temperature zone to the other temperature zones. The geometric position vector and the direction of the main airflow velocity vector are used to perform a vector dot product operation to obtain the cosine of the included angle.

5. The adaptive temperature control method for metallurgical nodes based on multi-zone coordinated regulation according to claim 1, characterized in that, The overall thermal coupling degree satisfies the following relationship: ； in, It is the first The temperature zone is in The overall thermal coupling degree at any given moment; It is a symbolic function; It is the first In the temperature gradient matrix at time t, the th Temperature zone for the first The thermal driving gradient of the temperature zone; yes The modulus length; It is the coupling dimension conversion factor; It is the weighting coefficient for the influence of airflow; It is the first Temperature zone points to the first The angle between the vector connecting the temperature zones and the direction of the main airflow velocity vector; This is the total number of sensors.

6. The adaptive temperature control method for metallurgical nodes based on multi-zone coordinated regulation according to claim 1, characterized in that, The calculation of the temperature rise slope based on the temperature at the current time and the previous time includes: The temperature difference is obtained by performing a difference calculation between the current temperature and the previous temperature. The ratio of the temperature difference to the preset sampling period is denoted as the temperature rise slope.

7. The adaptive temperature control method for metallurgical nodes based on multi-zone coordinated regulation according to claim 1, characterized in that, The heat trend operator satisfies the following relation: ； in, It is the first The temperature zone is in The heat trend operator in real time; , They are the first The temperature zone is in , Temperature at any moment; It is the preset sampling period; It is the reference temperature; It is the first The temperature zone is in The rate of change of the temperature rise slope within the neighborhood window at any given time; It is the natural logarithm function; It is the absolute value symbol.

8. The adaptive temperature control method for metallurgical nodes based on multi-zone coordinated regulation according to claim 7, characterized in that, The acquisition of the rate of change of the temperature rise slope within the neighborhood window includes: Within the neighborhood window corresponding to the current time, the temperature rise slope of each adjacent time is differentially calculated to obtain the temperature rise slope change sequence of the neighborhood window. The temperature rise slope change sequence is statistically averaged to obtain the temperature rise slope change rate within the neighborhood window.

9. The adaptive temperature control method for metallurgical nodes based on multi-zone coordinated regulation according to claim 1, characterized in that, The final execution power satisfies the following relationship: ； in, It is the first The temperature zone is in The final execution power at any given moment; It is the first The temperature zone is in Original power value at any given time; , They are the first The temperature zone is in The overall thermal coupling degree and thermal trend operator at any given time; It is the sensitivity conversion coefficient; It is the hyperbolic tangent function; , These are functions for finding the maximum and minimum values, respectively.

10. A metallurgical node temperature adaptive control system based on multi-zone coordinated regulation, characterized in that, The metallurgical node temperature adaptive control system based on multi-zone coordinated regulation includes a processor and a memory. The memory stores computer program instructions, which, when executed by the processor, implement the metallurgical node temperature adaptive control method based on multi-zone coordinated regulation according to any one of claims 1-9.

Images