A three-dimensional temperature measurement method for metallurgical furnace based on spatial acoustic wave correction method

By arranging acoustic transducers inside the metallurgical furnace to capture the flight time of sound waves and reconstructing the three-dimensional temperature field, the problems of accuracy and stability of temperature detection in the high-temperature environment of the metallurgical furnace are solved, and efficient temperature measurement is achieved.

CN119756621BActive Publication Date: 2025-11-25CINF ENG CO LTD
View PDF 2 Cites 0 Cited by

Patent Information

Application Number
CN202411780194.1
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2024-12-05
Publication Date
2025-11-25
Estimated Expiration
2044-12-05

AI Technical Summary

Technical Problem

Existing technologies struggle to achieve accurate and stable temperature detection in high-temperature environments within metallurgical furnaces. Traditional methods are costly and produce inconsistent results, while laser vision methods are susceptible to interference, leading to errors.

Method used

A three-dimensional temperature measurement method based on spatial acoustic correction is adopted. By arranging acoustic transducers in the metallurgical furnace, the flight time of acoustic waves is captured, the three-dimensional temperature field of the furnace is reconstructed, and the temperature is measured by utilizing the positive correlation between the propagation speed of acoustic waves and temperature.

Benefits of technology

It has achieved stability and accuracy in metallurgical furnace temperature detection, meeting industrial control requirements, reducing detection costs, and improving real-time performance.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119756621B_ABST
    Figure CN119756621B_ABST
Patent Text Reader

Abstract

The application discloses a kind of based on spatial acoustic wave correction method's metallurgical furnace three-dimensional temperature measurement method, comprising the following steps: S1.in the metallurgical furnace hearth region to be measured arrangement acoustic transducer;S2.the acoustic transducer of step S1 arrangement is used, acoustic wave flight time capture is carried out, and acoustic wave signal flight time is obtained;S3.according to the acoustic wave signal flight time obtained in step S2, reconstructs space temperature field based on space coefficient matrix, and obtains metallurgical furnace hearth three-dimensional temperature field;S4.according to the three-dimensional temperature field of metallurgical furnace hearth obtained in step S3, complete the measurement of metallurgical furnace three-dimensional temperature.The three-dimensional temperature measurement method of metallurgical furnace disclosed in the application obtains three-dimensional temperature information in metallurgical furnace hearth in real time accurately by the optimization capture of acoustic wave flight time and space temperature field reconstruction, improves the stability of metallurgical furnace temperature detection.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention belongs to the field of intelligent manufacturing in non-ferrous metal smelting, and specifically relates to a three-dimensional temperature measurement method for metallurgical furnaces based on spatial acoustic correction method. Background Technology

[0002] Non-ferrous metals refer to all metals other than iron, manganese, chromium, and iron-based alloys. They can be divided into heavy metals, light metals, precious metals, and rare metals. In a broader sense, non-ferrous metals also include non-ferrous alloys, which are alloys composed of non-ferrous metals as a base and one or more other elements added. Non-ferrous metals are fundamental materials for national economic development. Most industries, including aviation, aerospace, automobiles, machinery manufacturing, power, communications, construction, and home appliances, rely on non-ferrous metals for production. With the rapid advancement of modern chemical industry, agriculture, and science and technology, the role of non-ferrous metals in human development is becoming increasingly important.

[0003] Non-ferrous metal smelting can be broadly categorized into all-pyrometallurgical smelting, all-hydrometallurgical smelting, and combined hydrometallurgical and pyrometallurgical smelting. Among these, the pyrometallurgical smelting process requires a metallurgical furnace. The metallurgical furnace is the core equipment in the pyrometallurgical smelting process, and its operational quality directly affects the direct metal recovery rate and overall recovery rate, making it a key focus of research in non-ferrous metal intelligent manufacturing. Accurate temperature detection within the metallurgical furnace has long been a challenge in the industry. On the one hand, the furnace temperature is high, containing a molten liquid of glassy and metallic phases that can interfuse with most materials. Existing contact transmitters are insufficient, with most being single-use, costly, and producing inconsistent results. On the other hand, the harsh environment inside the furnace interferes with laser vision-based detection methods, particularly with damage to the laser emitter and camera leading to errors, making the obtained results unsuitable for intelligent control of non-ferrous metal pyrometallurgical processes. Summary of the Invention

[0004] The purpose of this invention is to provide a three-dimensional temperature measurement method for metallurgical furnaces based on spatial acoustic correction, which improves the stability, accuracy and real-time performance of metallurgical furnace temperature detection and meets industrial control requirements.

[0005] This invention provides a three-dimensional temperature measurement method for metallurgical furnaces based on spatial acoustic correction, comprising the following steps:

[0006] S1. Arrange acoustic transducers in the area to be tested in the furnace chamber of a metallurgical furnace;

[0007] S2. Using the acoustic transducer arranged in step S1, the acoustic wave transit time is captured to obtain the acoustic wave signal transit time.

[0008] S3. Based on the flight time of the acoustic signal obtained in step S2, the spatial temperature field is reconstructed using the spatial coefficient matrix to obtain the three-dimensional temperature field of the metallurgical furnace chamber.

[0009] S4. Based on the three-dimensional temperature field of the metallurgical furnace obtained in step S3, complete the measurement of the three-dimensional temperature of the metallurgical furnace.

[0010] Step S1 is as follows: Several layers of acoustic transducers are arranged at equal intervals in the area to be tested in the metallurgical furnace. Each layer of acoustic transducers is installed at equal intervals on the inner wall of the metallurgical furnace, and the distance is equal to the interlayer spacing between two adjacent layers of acoustic transducers. The acoustic transducers are wrapped by cooling water jackets. A plane rectangular coordinate system is established by selecting a direction parallel to the inner wall of the furnace as the x-axis on the horizontal plane and numbering the acoustic transducers. Each layer of acoustic transducers is numbered in alphabetical order from top to bottom.

[0011] The sound wave transit time capture in step S2 specifically involves: designing a transmission signal, using a sound wave transducer to transmit the signal in a preset order, extracting the echo signal envelope, separating the peak time corresponding to the echo signal, and calculating the sound wave transit time by weighted averaging.

[0012] The specific signal transmission sequence of the acoustic transducer is as follows: Taking any vertical plane A within the furnace where the acoustic transducer is located as a reference, the acoustic transducer A is positioned at the center of that plane. ij The sound wave is emitted from the starting point. After one sound wave emission, the sound transducer B at the corresponding position on the adjacent plane B is taken in a clockwise direction. ij The sound waves are emitted. The sound wave transducers in the vertical plane of the furnace of the two metallurgical furnaces, marked A and B, complete the sound wave emission in one go according to the Fibonacci curve. When a sound wave transducer in a certain vertical plane emits a sound wave signal, the sound wave transducers in the same plane do not work, and the other sound wave transducers are all receiving ends.

[0013] The transmitted signal s(t) is expressed by the following formula:

[0014] s(t) = |sin(t)|

[0015] Where t is time and 0 ≤ t ≤ 2π;

[0016] The echo signal envelope is extracted using the Hilbert transform, expressed by the following formula:

[0017]

[0018] Where x(t) is the echo signal; H[x(t)] is the Hilbert transform of the echo signal; D is the signal transmission interval; n(τ) is the correction function, characterizing the average noise during the operation of the core equipment, the metallurgical furnace; the excitation process analytical signal is fitted using the impulse response function, expressed by the following formula:

[0019]

[0020] Where F(t) is the envelope function of the echo signal; z(t) is the projected phase of the echo signal; j is the imaginary operator; the envelope function F(t) is fitted with an impulse response function model and expressed by the following formula:

[0021]

[0022] Among them, z * The function is the envelope fitting function; F0 is the first characteristic coefficient; Δ is the second characteristic coefficient; δ is the third characteristic coefficient; α is the fourth characteristic coefficient; n * (τ) is the envelope amplitude rate correction function; the characteristic coefficient fitness function f(x) is designed and expressed by the following formula:

[0023] f(x) = s(t) * )-z *

[0024] Among them, t * For the standardized time variable and The optimal solution is obtained by using the differential evolution method to find the combination of eigenvalues ​​that minimizes the fitness function f(x), as expressed by the following formula:

[0025] X i =(X i,1 (N),X i,2 (N),X i,3 (N),X i,4 (N))

[0026] Among them, X i X is a potential optimal solution where 1 ≤ i ≤ M; M is the number of individuals in the solution space, and each individual is a 4-dimensional vector; i,1 (N) corresponds to the characteristic coefficients F0, X i,2 (N) corresponds to the characteristic coefficients Δ, X i,3 (N) corresponds to the characteristic coefficient δ, X i,4 (N) corresponds to the characteristic coefficient α; N is the number of iterations;

[0027] The optimized differential evolution method includes the following steps:

[0028] a. Adaptive mutation of the optimal solution:

[0029]

[0030] Among them, X i (N) represents the current optimal solution; This is the result after the optimal solution has been mutated; and The individual solution is randomly selected during the Nth iteration and K is the scaling factor, calculated using the following formula:

[0031]

[0032] Among them, f i (N) represents the characteristic coefficients, and the fitness function f(x) takes the value X. i (N) corresponds to fitness; f p (N) represents the characteristic coefficients, and the fitness function f(x) takes the value X. p (N) corresponds to fitness; f q (N) represents the characteristic coefficients, and the fitness function f(x) takes the value X. q (N) represents the fitness level;

[0033] b. Optimal solution adaptive crossover and selection:

[0034]

[0035] in, The individual solution after adaptive crossover; H i (N) is the average of the three solutions with the closest fitness; Ra is a random number in the interval [0,1]; cr is the crossover probability, expressed by the following formula:

[0036]

[0037] Where cr1 is the current crossover probability of the individual solution; To achieve the fitness of the individual solution after adaptive mutation; f i The average fitness of all current individual solutions; cr max This is the upper limit of the crossover probability cr; cr min Let be the lower bound of the crossover probability cr; use the following formula to represent the new individual solution X. i (N+1):

[0038]

[0039] in, To provide the individual solution after completing the mutation and crossover process fitness; f Xi The individual solution X before the mutation and crossover process i The fitness of;

[0040] c. After selection, repeat the mutation, crossover, and selection process until the number of iterations N reaches the preset limit.

[0041] Substitute the obtained optimal characteristic coefficients back into the envelope fitting function z. * And differentiate with respect to time t and let (z) *Given that z' = 0, we can solve for the time t1 and t2 when z(t) reaches its maximum and second maximum values, respectively. Then, the time t of the acoustic signal travels through a single receiver is t. f Express it using the following formula:

[0042]

[0043] In this context, Max{·} represents the maximum value operation, and Min{·} represents the minimum value operation.

[0044] The design selects vertical plane A, vertical plane B, and the bottom surface as references, dividing the measurement space into equal-sized α-planes. 3 A grid space is selected; the grid space adjacent to vertical plane A and vertical plane B in the lowest layer is selected as the initial reference number (1,1). Taking the positive direction of the x-axis, the grid is numbered according to the first digit of the number being the layer number, and the second digit being the number from left to right and from bottom to top, with the maximum number being (α,α). Taking the grid space adjacent to vertical plane A and vertical plane B in the lowest layer as the origin, the sound wave flight time in the grid is numbered from left to right and from bottom to top along the x-axis. The flight time of a single spatial grid is expressed by the following formula:

[0045]

[0046] Where, α ij Let t be the time t of the sound wave's journey. f The actual distance through the j-th spatial grid in the i-th layer; t fε α is the time of flight between the transmitter and receiver of the acoustic transducer during one measurement cycle, where ε is an integer and 0 < ε ≤ α. 3 .

[0047] Step S3 specifically involves constructing a three-dimensional temperature field of the metallurgical furnace based on the flight time of the sound wave signal in a single spatial grid obtained in step S2, since the propagation speed of sound waves is positively correlated with temperature.

[0048] This invention discloses a three-dimensional temperature measurement method for metallurgical furnaces based on spatial acoustic correction. By optimizing the capture of acoustic wave transit time and reconstructing the spatial temperature field, the method can obtain three-dimensional temperature information inside the furnace in real time and accurately, thereby improving the stability of temperature detection in metallurgical furnaces. Attached Figure Description

[0049] Figure 1 This is a flowchart illustrating the method of the present invention. Detailed Implementation

[0050] This invention provides a three-dimensional temperature measurement method for metallurgical furnaces based on spatial acoustic correction, the flowchart of which is shown below. Figure 1 As shown, it includes the following steps:

[0051] S1. An acoustic transducer is installed in the area to be tested within the metallurgical furnace chamber, specifically as follows:

[0052] Several layers of acoustic transducers are arranged at equal intervals in the test area of ​​the metallurgical furnace. Each layer of acoustic transducers is installed at equal intervals on the inner wall of the metallurgical furnace, and the distance is equal to the interlayer spacing between two adjacent layers of acoustic transducers. The acoustic transducers are wrapped by cooling water jackets. A plane rectangular coordinate system is established by selecting a direction parallel to the inner wall of the furnace as the x-axis on the horizontal plane, and the acoustic transducers are numbered according to this system. Each layer of acoustic transducers is numbered in alphabetical order from top to bottom.

[0053] S2. Using the acoustic transducer arranged in step S1, the acoustic wave transit time is captured to obtain the acoustic signal transit time, specifically:

[0054] The design involves transmitting signals using an acoustic transducer in a preset sequence. The echo signal envelope is extracted, the peak time corresponding to the echo signal is separated, and the sound wave transit time is obtained by weighted averaging.

[0055] The specific signal transmission sequence of the acoustic transducer is as follows: Taking any vertical plane A within the furnace where the acoustic transducer is located as a reference, the acoustic transducer A is positioned at the center of that plane. ij The sound wave is emitted from the starting point. After one sound wave emission, the sound transducer B at the corresponding position on the adjacent plane B is taken in a clockwise direction. ij Sound waves are emitted. Two acoustic transducers, labeled A and B, within the vertical plane of the metallurgical furnace chamber complete the sound wave emission in one pass according to the Fibonacci curve. When a certain vertical plane acoustic transducer emits a sound signal, other transducers on the same plane do not operate; all other transducers act as receivers. Taking the transducer selection sequence of the first 18 measurements as an example, the sequence of the transmitting transducers can be represented as: A ij B ij A (i+1 ) j B (i+1 ) j A (i+1)(j+1) B (i+1)(j+1) A i(j+1) B i(j+1) A (i-1)(j+1) B (i-1)(j+1) A (i-1)j B (i-1)j A (i-1)(j-1) B (i-1)(j-1) A i(j-1) B i(j-1) A (i+1)(j-1) B (i+1)(j-1) .

[0056] The transmitted signal s(t) is expressed by the following formula:

[0057] s(t) = |sin(t)|

[0058] Where t is time and 0 ≤ t ≤ 2π;

[0059] The echo signal envelope is extracted using the Hilbert transform, expressed by the following formula:

[0060]

[0061] Where x(t) is the echo signal; H[x(t)] is the Hilbert transform of the echo signal; D is the signal transmission interval; n(τ) is the correction function, characterizing the average noise during the operation of the core equipment, the metallurgical furnace; the excitation process analytical signal is fitted using the impulse response function, expressed by the following formula:

[0062]

[0063] Where F(t) is the envelope function of the echo signal; z(t) is the projected phase of the echo signal; j is the imaginary operator; the envelope function F(t) is fitted with an impulse response function model and expressed by the following formula:

[0064]

[0065] Among them, z * The function is the envelope fitting function; F0 is the first characteristic coefficient; Δ is the second characteristic coefficient; δ is the third characteristic coefficient; α is the fourth characteristic coefficient; n * (τ) is the envelope amplitude rate correction function; the characteristic coefficient fitness function f(x) is designed and expressed by the following formula:

[0066] f(x) = s(t) * )-z *

[0067] Among them, t * For the standardized time variable and The optimal solution is obtained by using the differential evolution method to find the combination of eigenvalues ​​that minimizes the fitness function f(x), as expressed by the following formula:

[0068] X i =(X i,1 (N),X i,2 (N),X i,3 (N),X i,4 (N))

[0069] Among them, X i X is a potential optimal solution where 1 ≤ i ≤ M; M is the number of individuals in the solution space, and each individual is a 4-dimensional vector; i,1(N) corresponds to the characteristic coefficients F0, X i,2 (N) corresponds to the characteristic coefficients Δ, X i,3 (N) corresponds to the characteristic coefficient δ, X i,4 (N) corresponds to the characteristic coefficient α; N is the number of iterations;

[0070] The optimized differential evolution method includes the following steps:

[0071] a. Adaptive mutation of the optimal solution:

[0072]

[0073] Among them, X i (N) represents the current optimal solution; This is the result after the optimal solution has been mutated; and The individual solution is randomly selected during the Nth iteration and K is the scaling factor, calculated using the following formula:

[0074]

[0075] Among them, f i (N) represents the characteristic coefficients, and the fitness function f(x) takes the value X. i (N) corresponds to fitness; f p (N) represents the characteristic coefficients, and the fitness function f(x) takes the value X. p (N) corresponds to fitness; f q (N) represents the characteristic coefficients, and the fitness function f(x) takes the value X. q (N) represents the fitness level;

[0076] b. Optimal solution adaptive crossover and selection:

[0077]

[0078] in, The individual solution after adaptive crossover; H i (N) is the average of the three solutions with the closest fitness; Ra is a random number in the interval [0,1]; cr is the crossover probability, expressed by the following formula:

[0079]

[0080] Where cr1 is the current crossover probability of the individual solution; To achieve the fitness of the individual solution after adaptive mutation; f i The average fitness of all current individual solutions; cr max This is the upper limit of the crossover probability cr; cr minLet be the lower bound of the crossover probability cr; use the following formula to represent the new individual solution X. i (N+1):

[0081]

[0082] in, To provide the individual solution after completing the mutation and crossover process The fitness of; The individual solution X before the mutation and crossover process i The fitness of;

[0083] c. After selection, repeat the mutation, crossover, and selection process until the number of iterations N reaches the preset limit.

[0084] Substitute the obtained optimal characteristic coefficients back into the envelope fitting function z. * And differentiate with respect to time t and let (z) * Given that z' = 0, we can solve for the time t1 and t2 when z(t) reaches its maximum and second maximum values, respectively. Then, the time t of the acoustic signal travels through a single receiver is t. f Express it using the following formula:

[0085]

[0086] In this context, Max{·} represents the maximum value operation, and Min{·} represents the minimum value operation.

[0087] The design selects vertical plane A, vertical plane B, and the bottom surface as references, dividing the measurement space into equal-sized α-planes. 3 A grid space is selected; the grid space adjacent to vertical plane A and vertical plane B in the lowest layer is selected as the initial reference number (1,1). Taking the positive direction of the x-axis, the grid is numbered according to the first digit of the number being the layer number, and the second digit being the number from left to right and from bottom to top, with the maximum number being (α,α). Taking the grid space adjacent to vertical plane A and vertical plane B in the lowest layer as the origin, the sound wave flight time in the grid is numbered from left to right and from bottom to top along the x-axis. The flight time of a single spatial grid is expressed by the following formula:

[0088]

[0089] Where, α ij Let t be the time t of the sound wave's journey. f The actual distance through the j-th spatial grid in the i-th layer; t fε α is the time of flight between the transmitter and receiver of the acoustic transducer during one measurement cycle, where ε is an integer and 0 < ε ≤ α. 3 .

[0090] S3. Based on the acoustic signal transit time obtained in step S2, and using the spatial coefficient matrix, reconstruct the spatial temperature field to obtain the three-dimensional temperature field of the metallurgical furnace chamber, specifically:

[0091] Since the speed of sound propagation is positively correlated with temperature, the three-dimensional temperature field of the metallurgical furnace is constructed based on the flight time of the sound wave signal of a single spatial grid obtained in step S2.

[0092] S4. Based on the three-dimensional temperature field of the metallurgical furnace obtained in step S3, complete the measurement of the three-dimensional temperature of the metallurgical furnace.

[0093] The method of the present invention will be described below with reference to an embodiment:

[0094] Choosing a pyrometallurgical copper smelting furnace as the application scenario, and selecting a single acoustic transducer with a power of 15kHz and 2600W, five layers of acoustic transducers with five transducers per layer are evenly spaced around the inner wall of the furnace in the area to be measured. The sequence of transducer activation for the first 18 measurements in one cycle is: A 33 B 33 A 34 B 34 A 24 B 24 A 23 B 23 A 22 B 22 A 32 B 32 A 42 B 42 A 43 B 43 A 44 B 44 The impulse response function model fitting envelope function optimization process is set with a population size of 200, a total number of iterations of 1000, an initial scaling factor of 0.1, and an initial crossover probability of 0.5. Through actual measurements using disposable cement thermocouples, the method described in this invention has an error within 15%, meeting the requirements of core equipment for non-ferrous metal pyrometallurgical processes.

Claims

1. A three-dimensional temperature measurement method for a metallurgical furnace based on spatial acoustic correction, characterized in that, Includes the following steps: S1. Arrange acoustic transducers in the area to be tested in the furnace chamber of a metallurgical furnace; S2. Using the acoustic transducer arranged in step S1, the acoustic wave transit time is captured to obtain the acoustic wave signal transit time. S3. Based on the flight time of the acoustic signal obtained in step S2, the spatial temperature field is reconstructed using the spatial coefficient matrix to obtain the three-dimensional temperature field of the metallurgical furnace chamber. S4. Based on the three-dimensional temperature field of the metallurgical furnace obtained in step S3, complete the measurement of the three-dimensional temperature of the metallurgical furnace; The sound wave transit time capture is specifically as follows: design the transmission signal, use a sound wave transducer to transmit the signal in a preset order, extract the echo signal envelope, separate the peak time corresponding to the echo signal, and calculate the sound wave transit time by weighted average. Transmit signal Express it using the following formula: in, For time and ; The echo signal envelope is extracted using the Hilbert transform, expressed by the following formula: in, This is an echo signal; The Hilbert transform of the echo signal; This is the signal transmission range; For integration variables; The correction function characterizes the average noise during the operation of the core equipment, the metallurgical furnace; the impulse response function is used to fit the analytical signal of the excitation process, expressed by the following formula: in, The envelope function of the echo signal; The projected phase of the echo signal; The imaginary operator is used; the envelope function is fitted using an impulse response function model. It can be expressed by the following formula: in, This is the function for fitting the envelope. The first characteristic coefficient; The second characteristic coefficient; It is the third characteristic coefficient; It is the fourth characteristic coefficient; Design the envelope amplitude rate correction function; design the characteristic coefficient fitness function. It can be expressed by the following formula: in, For the standardized time variable and ; express The transmitted signal at time; the fitness function with characteristic coefficients is solved using the optimized differential evolution method. The combination of characteristic coefficients that minimizes the value yields the potential optimal solution, expressed by the following formula: in, It is a potential optimal solution and ; The solution space contains the number of individuals, each of which is a 4-dimensional vector; Corresponding characteristic coefficients , Corresponding characteristic coefficients , Corresponding characteristic coefficients , Corresponding characteristic coefficients ; The number of iterations; Substitute the obtained optimal characteristic coefficients back into the envelope fitting function. and on time Differentiate and let Solve Obtain Time to obtain the maximum and second maximum values , The time of flight of the acoustic signal at a single receiver Express it using the following formula: in, This is for retrieving the maximum value. This is an operation to find the minimum value.

2. The three-dimensional temperature measurement method for metallurgical furnaces based on spatial acoustic correction method according to claim 1, characterized in that, Step S1 specifically involves: arranging several layers of acoustic transducers at equal intervals in the test area of ​​the metallurgical furnace. Each layer of acoustic transducers is installed at equal intervals on the inner wall of the metallurgical furnace, with the distance equal to the interlayer spacing between adjacent layers; the acoustic transducers are encased in cooling water jackets; and a direction parallel to the inner wall of the furnace is selected on the horizontal plane as... Establish a Cartesian coordinate system and number the acoustic transducers.

3. The three-dimensional temperature measurement method for metallurgical furnaces based on spatial acoustic correction method according to claim 2, characterized in that, The specific signal transmission sequence of the acoustic transducer is as follows: Select any vertical plane of the furnace where the acoustic transducer is located. Using this as a reference, the acoustic transducer is located at the center of the plane. The sound wave is emitted from the starting point, and after one sound wave emission, the adjacent plane is taken in a clockwise direction. Acoustic transducer at corresponding position Emit sound waves, marked as , The acoustic transducers in the vertical plane of the two metallurgical furnaces complete the acoustic wave transmission in one go according to the Fibonacci curve; when a certain vertical plane acoustic transducer transmits an acoustic wave signal, the acoustic transducers in the same plane do not work, and the other acoustic transducers are all receiving ends.

4. The three-dimensional temperature measurement method for metallurgical furnaces based on spatial acoustic correction method according to claim 3, characterized in that, The optimized differential evolution method includes the following steps: a. Adaptive mutation of the optimal solution: in, This is the current optimal solution; This is the result after the optimal solution has been mutated; and For the first The individual solution is randomly selected during the second iteration and ; The scaling factor is calculated using the following formula: in, The fitness function of characteristic coefficients Pick Corresponding fitness; The fitness function of characteristic coefficients Pick Corresponding fitness; The fitness function of characteristic coefficients Pick Corresponding fitness; Indicates a logical intersection relationship; Represents a logical union relation; b. Optimal solution adaptive crossover and selection: in, The individual solution after adaptive crossover; The average of the solutions of the three individuals with the closest fitness; For interval Random numbers; The crossover probability is expressed by the following formula: in, Solve for the current crossover probability for each individual; To achieve the fitness of the individual solution after adaptive mutation; This represents the average fitness of all current individual solutions. Crossover probability The upper limit; Crossover probability The lower limit value; use the following formula to represent the new individual solution. : in, To provide the individual solution after completing the mutation and crossover process The fitness of; Individual solutions before the mutation and crossover process The fitness of; c. After selection, repeat the mutation, crossover, and selection process until the specified number of iterations is reached. The preset parameters have been achieved.

5. The three-dimensional temperature measurement method for metallurgical furnaces based on spatial acoustic correction method according to claim 4, characterized in that, The design selects a vertical plane. Vertical plane Using the bottom surface as a reference, the measurement space is divided into equal-sized sections. Each grid space; Select the lowest space and the vertical plane and vertical plane Adjacent grid spaces are used as the initial reference numbering. ,Pick In the positive direction of the axis, the numbers are assigned according to the rule that the first digit represents the layer number, and the second digit is numbered from left to right and bottom to top, with the largest number being [number missing]. ; in the lowest space and the vertical plane and vertical plane The adjacent grid space is the origin, along... The axes are arranged from left to right and from bottom to top to encode the time of sound wave transit within the grid. The flight time of a single spatial grid is expressed by the following formula: in, Time of travel for sound waves Through the first Layer The actual distance of each spatial grid; The time of flight between the transmitter and receiver of the acoustic transducer within one measurement cycle is denoted as . Integer and .

6. The three-dimensional temperature measurement method for metallurgical furnaces based on spatial acoustic correction method according to claim 5, characterized in that, Step S3 specifically involves constructing a three-dimensional temperature field of the metallurgical furnace based on the flight time of the sound wave signal in a single spatial grid obtained in step S2, since the propagation speed of sound waves is positively correlated with temperature.

Citation Information

Patent Citations

  • Detonation flame temperature measurement system and three-dimensional temperature field reconstruction method

    CN106289568A

  • Pipeline multi-source electromagnetic acoustic emission feature extraction method based on sparse optimization

    CN118760877A