TWI731812B - Blast furnace temperature estimating method, computer program product, and computer readable recording medium - Google Patents

Abstract

A blast furnace temperature estimating method, computer program product, and computer readable recording medium are disclosed. The method includes: driving one of a plurality of sonic wave transceivers as a sound source to emit a sonic wave; perform, on at least one of the remaining transceiver units except the sound source: sampling continuously the sonic wave to generate a sonic wave signal as a target signal; converting the target signal into a frequency-domain signal, setting a maximum amplitude peak of the frequency-domain signal to a maximum audio intensity, calculating a product of the maximum audio intensity and an azimuth proportional constant as a threshold; searching a time when the first one of amplitude peaks of the target signal exceeds the threshold being found as a rising-peak time; and calculating a sonic-wave-path temperature based on the rising-peak time, a sound-emission time, a distance of sonic wave transmission and reception, and an acoustic constant.

Description

Method for estimating temperature of blast furnace, computer program product and computer readable recording medium

The present invention relates to a temperature estimation technology, in particular to a method for estimating the temperature of a blast furnace, a computer program product and a computer readable recording medium.

In some high-temperature environments, due to the many interference sources in the environment, it is usually difficult to use existing temperature measurement components to meet the temperature measurement requirements.

Take the temperature measurement of an iron-making blast furnace as an example, by reading the value of the suction thermocouple on the temperature measuring rod on the top of the furnace, the distribution of the airflow temperature in the furnace can be understood as whether the hot air field in the furnace is uniform and whether the number of cloth loops needs to be adjusted The basis.

However, the hot air flow in the blast furnace has a high dust content, and the exhaust thermocouple on the temperature measuring rod is often blocked, which affects the operator to judge the furnace condition.

In view of this, it is necessary to provide a technical solution different from the past to solve the problems existing in the conventional technology.

An object of the present invention is to provide a method for estimating the temperature of a blast furnace. Based on the sound wave transmission characteristics, the sound wave path temperature is estimated by using the sound wave transmitting and receiving time in the blast furnace environment, which is beneficial to meet the temperature measurement requirements of the blast furnace.

The second purpose of the present invention is to provide a computer program product, based on the sound wave transmission characteristics, use the sound wave receiving and sending time in the blast furnace environment to estimate the sound wave path temperature, which is beneficial to meet the temperature measurement requirements of the blast furnace.

Another object of the present invention is to provide a computer-readable recording medium, based on the sound wave transmission characteristics, using the sound wave receiving and sending time in the blast furnace environment to estimate the sound wave path temperature, which is beneficial to meet the temperature measurement requirements of the blast furnace.

To achieve the above objective, one aspect of the present invention provides a method for estimating the temperature of a blast furnace, which is configured by a processor coupled to a memory and a plurality of transceiver units, the plurality of transceiver units in a test environment A horizontal plane is arranged equidistantly along a circular track, the processor executes instructions stored in the memory, and the method includes the steps of: driving one of the plurality of transceiver units as a sound source to emit a sound wave at a sounding time ; And perform at least one of the remaining transceiver units other than the sound source: continuous sampling of the sound wave to generate a sound wave signal as a target signal; convert the target signal into a frequency domain signal, and a part of the frequency domain signal The maximum amplitude peak value is set to an audio maximum intensity; the product of the maximum audio intensity and the azimuth proportional constant is calculated as a threshold; the time when the first amplitude peak of at least one amplitude peak of the target signal exceeds the threshold is found , As a rising peak time; and calculating a sound wave path temperature based on the rising peak time, the sounding time, a sound wave transmitting and receiving distance, and an acoustic constant.

其中，T為該音波路徑溫度；D為該音波收發距離，該音波收發距離為該目標訊號所屬的收發單元與發出該音波的收發單元之間的一距離；△t為該上昇峰值時間減去該發聲時間的一時間差值；A為該聲學常數。 In an embodiment of the present invention, the calculation method of the acoustic path temperature is as follows:

Among them, T is the sound wave path temperature; D is the sound wave transmitting and receiving distance, the sound wave transmitting and receiving distance is the distance between the transceiver unit to which the target signal belongs and the transceiver unit that emits the sound wave; △t is the rising peak time minus A time difference of the sounding time; A is the acoustic constant.

In an embodiment of the present invention, the azimuth proportional constant is less than 1, and the azimuth proportional constant has a positively correlated proportional relationship with the sound wave transmitting and receiving distance. The greater the sound wave transmitting and receiving distance, the larger the azimuth proportional constant.

In an embodiment of the present invention, the target signal is converted into the frequency domain signal after passing through a band-pass filtering process.

In an embodiment of the present invention, the target signal is converted into the frequency domain signal according to a fast Fourier algorithm.

To achieve the above objective, another aspect of the present invention provides a computer program product. After the computer program is loaded and executed, the computer can execute the method for estimating the temperature of a blast furnace as described above.

In order to achieve the above objective, another aspect of the present invention provides a computer-readable recording medium. The computer can read a program stored in the recording medium. After the computer loads and executes the program, the computer can complete the above-mentioned Method of estimating the temperature of a blast furnace.

The method for estimating the temperature of a blast furnace, the computer program product and the computer readable recording medium of the present invention convert the target signal into the frequency domain signal, and the maximum amplitude peak value of the frequency domain signal Set as the maximum intensity of the audio; calculate the product of the maximum intensity of the audio and the azimuth proportional constant as the threshold; find the time when the first amplitude peak of the target signal exceeds the threshold as the rising peak time; and basis The rising peak time, the sounding time, the sound wave transmitting and receiving distance, and the acoustic constant are used to calculate the sound wave path temperature. In this way, the sound wave path temperature can be estimated based on the sound wave transmission characteristics using the sound wave receiving and sending time in the blast furnace environment, which can improve the accuracy of the calculation of the flight time from the sound wave transmitter to the receiver, and reduce the average temperature of the medium on the sound wave path. The error is conducive to meeting the temperature measurement requirements of the blast furnace.

F2: Frequency domain signal

F7: Frequency domain signal

F8: Frequency domain signal

G: Target signal

G2: Target signal

G7: Target signal

G8: Target signal

K2: The first amplitude peak

K7: The first amplitude peak

K8: The first amplitude peak

L: Filtered signal

P2: Peak

P7: Peak

P8: Peak

P2’: Threshold

P7’: Threshold

P8’: Threshold

R1: Transceiver unit

R2: Transceiver unit

R3: Transceiver unit

R4: Transceiver unit

R5: Transceiver unit

R6: Transceiver unit

R7: Transceiver unit

R8: Transceiver unit

R9: Transceiver unit

R10: Transceiver unit

R11: Transceiver unit

R12: Transceiver unit

E: bearing

EN:Azimuth

ES: Orientation

SE: bearing

S: bearing

SW: bearing

WS: Orientation

W: bearing

WN: Orientation

NW: Orientation

N: bearing

NE: bearing

T1: voice step

T2: temperature measurement step

d12: Sending and receiving distance

d13: Sending and receiving distance

d14: Sending and receiving distance

d15: Sending and receiving distance

d16: Sending and receiving distance

d17: Sending and receiving distance

d18: Sending and receiving distance

d19: Sending and receiving distance

d1a: Sending and receiving distance

d1b: receiving and sending distance

d1c: receiving and sending distance

[Figure 1]: A schematic flow diagram of a method for estimating the temperature of a blast furnace according to an embodiment of the present invention.

[Figure 2]: A schematic diagram of the layout of the transceiver unit of the method for estimating the temperature of a blast furnace according to an embodiment of the present invention.

[Figure 3]: A schematic diagram of the signal form of the method for estimating the temperature of a blast furnace according to an embodiment of the present invention.

[Figure 4]: A schematic diagram of the rising peak time of the method for estimating the temperature of a blast furnace according to an embodiment of the present invention.

[Figure 5]: A schematic diagram of the target signal of the method for estimating the temperature of a blast furnace according to an embodiment of the present invention.

[Figure 6]: A schematic diagram of the filtered signal of the method for estimating the temperature of a blast furnace according to an embodiment of the present invention.

In order to make the above and other objectives, features, and advantages of the present invention more obvious and understandable, the preferred embodiments of the present invention will be described in detail below in conjunction with the accompanying drawings. Furthermore, the directional terms mentioned in the present invention, such as up, down, top, bottom, front, back, left, right, inside, outside, side, surrounding, center, horizontal, horizontal, vertical, vertical, axial, Radial, uppermost or lowermost layer, etc., only Refers to the direction of the attached drawings. Therefore, the directional terms used are used to describe and understand the present invention, rather than to limit the present invention.

An aspect of the present invention provides a method for estimating the temperature of a blast furnace, which can be configured as a temperature measuring device, for example, a processor coupled to a memory and several transceiver units (such as sonic transceivers). The several transceiver units are arranged equidistantly along a circular track in a horizontal plane in an environment to be tested (such as ironmaking blast furnace mouth, but not limited to this), and the processor executes the instructions stored in the memory , For performing the method for estimating the temperature of the blast furnace, as shown in Figure 1, the method includes a sounding step T1 and a temperature measurement step T2. The implementation of the method is illustrated below, but not limited to this.

Please refer to FIG. 1 again. The sounding step T1 can, for example, drive one of the transceiver units as a sound source to emit a sound wave at a sounding time.

For example, for the convenience of description, as shown in Figure 2, taking 12 transceiver units as an example, the positions of the 12 transceiver units R1~R12 can be marked as: EN, E, ES, SE, S, SW, WS, W, WN, NW, N, and NE are used to facilitate identification, but not limited to this, and the number of the transceiver unit can also be 2 or more.

As shown in Figure 2, taking the transceiver unit R1 in the position EN as the sound source to emit the sound wave as an example, in the other positions E, ES, SE, S, SW, WS, W, WN, NW, N, NE The transceiver units R2~R12 can continuously sample the sound wave to generate a sound signal each, but not limited to this. The transceiver unit R1 can also continuously sample the sound wave, but if the sound source (ie the transceiver unit) The sound wave receiving and sending distance (0) of R1) is used as a temperature measurement basis, then the transceiver unit R1 may not continuously sample the sound wave, or the processor may ignore the signal generated by the transceiver unit R1 continuously sampling the sound wave.

As shown in Figure 2, there is a separate unit between the transceiver unit R1 in the orientation EN and the transceiver units R2~R12 in the other orientations E, ES, SE, S, SW, WS, W, WN, NW, N, and NE. Sound wave transmission and reception distances d12, d13, d14, d15, d16, d17, d18, d19, d1a, d1b, d1c.

As shown in Fig. 2, the temperature measurement step T2 can be performed, for example, on at least one of the other transceiver units except the sound source: the sound wave is continuously sampled to generate a sound wave signal as a target signal; The target signal is converted into a frequency domain signal (for example, the frequency domain signals F2, F8, and F7 converted from the target signals G2, G8, and G7 shown in Figure 3 have peaks P2, P8, P7), and one of the frequency domain signals The maximum amplitude peak value is set to an audio maximum intensity; the product of the maximum audio intensity and the azimuth proportional constant is calculated as a threshold; the time when the first amplitude peak of at least one amplitude peak of the target signal exceeds the threshold is found (For example, the time when the first amplitude peaks K2, K8, and K7 of the target signals G2, G8, and G7 exceed the thresholds P2', P8', P7' shown in Figure 4) is used as a rising peak time; and A sound wave path temperature (Path_T) is calculated by rising peak time, the sounding time, a sound wave transmitting and receiving distance, and an acoustic constant. Wherein, the temperature measurement step T2 can perform the above-mentioned actions on at least one of the remaining transceiver units except the sound source to obtain the temperature of a position on a specific sound wave path; but not limited to this, the temperature measurement step T2 can also perform the above actions in sequence on each of the remaining transceiver units other than the sound source to learn the temperature distribution in the horizontal plane.

Thereby, by finding the time when the first amplitude peak of the at least one amplitude peak of the target signal exceeds the threshold as the rising peak time, it is possible to overcome the environmental noise generated by the airflow in the furnace and cause the misjudgment of the temperature measurement result to be based on Acoustic wave transmission characteristics, use the sound wave receiving and sending time in the blast furnace environment to estimate the sound wave path temperature, which can improve the accuracy of the calculation of the flight time from the sound wave transmitter to the receiving end, and reduce the error caused by the average temperature of the medium on the path of the sound wave, which is conducive to compliance Blast furnace temperature measurement requirements.

其中，T為該音波路徑溫度；D為該音波收發距離，該音波收發距離為該目標訊號所屬的收發單元與發出該音波的收發單元之間的一距離，諸如該音波收發距離d12、d13、d14、d15、d16、d17、d18、d19、d1a、d1b或d1c(如第2圖所示)；△t為該上昇峰值時間減去該發聲時間的一時間差值；A為該聲學常數，例如A=(r．R)/M，r為熱容比(heat capacity ratio)，R為理想氣體常數(gas constanst)，M為摩爾質量(molecular mass of gas)，該熱容比、理想氣體常數及摩爾質量的適用及調整係其所屬領域中具有通常知識者可以理解，在此容不贅述。 In one embodiment, the calculation method of the acoustic path temperature is as shown in the following formula:

Among them, T is the sound wave path temperature; D is the sound wave transceiver distance, the sound wave transceiver distance is the distance between the transceiver unit to which the target signal belongs and the transceiver unit that emits the sound wave, such as the sound wave transceiver distance d12, d13, d14, d15, d16, d17, d18, d19, d1a, d1b or d1c (as shown in Figure 2); Δt is the time difference of the rising peak time minus the sounding time; A is the acoustic constant, For example, A=(r.R)/M, r is the heat capacity ratio, R is the gas constanst, and M is the molecular mass of gas. The heat capacity ratio is the ideal gas The application and adjustment of constants and molar masses can be understood by those with ordinary knowledge in the field, so I won’t repeat them here.

In one embodiment, the azimuth proportional constant is less than 1, and the azimuth proportional constant is in a positively correlated proportional relationship with the sound wave transmitting and receiving distance. The greater the sound wave transmitting and receiving distance, the greater the azimuth proportional constant.

For example, as shown in Figure 2, the transceiver unit R1 at the position EN emits the sound wave, and the farther away the other transceiver units (ie, the transceiver units that receive the sound waves) from the transceiver unit R1 that emit the sound waves, such as R2~R12 ), the greater the azimuth proportional constant. For example, if the transceiver unit R1 at the position EN emits the sound wave, the azimuth proportional constant of the transceiver unit R7 at the azimuth WS is C; the azimuth proportional constant of the transceiver unit R6 at the azimuth SW and the transceiver unit R8 at the azimuth W is W1 *C; the azimuth proportional constant of the transceiver unit R5 at the azimuth S and the transceiver unit R9 at the azimuth WN is W2*C; the azimuth proportional constant of the transceiver unit R4 at the azimuth SE and the transceiver unit R10 at the azimuth NW is W3*C; The azimuth proportional constant of the transceiving unit R3 in the azimuth ES and the transceiving unit R11 in the azimuth N is W4*C; the azimuth proportional constant of the transceiving unit R2 in the azimuth E and the transceiving unit R12 in the azimuth NE is W5*C. Among them, C is a non-zero positive value less than 1; W1, W2, W3, W4, and W5 are all non-zero positive values less than 1, and W1, W2, W3, W4, and W5 are not equal, for example: W1>W2>W3>W4>W5.

In one embodiment, as shown in FIGS. 5 and 6, the target signal G may undergo a band-pass filtering process to first generate a filtered signal L, and then convert the filtered signal L into the frequency domain signal.

In one embodiment, the target signal can be converted into the frequency domain signal according to a fast Fourier (FFT) algorithm. For example, as shown in Figure 5, the target signal G2 generated by the transceiver unit R2 at the position E It can be converted into a frequency domain signal F2, the maximum amplitude peak value P2 of the frequency domain signal F2 is set to the maximum intensity of the audio; the target signal G8 generated by the transceiver unit R8 in the azimuth W can be converted into a frequency domain signal F8, the frequency domain signal F8 The maximum amplitude peak value P8 of the frequency domain signal is set to the maximum intensity of the audio; the target signal G7 generated by the transceiver unit R7 in the azimuth WS can be converted into a frequency domain signal F7, and the maximum amplitude peak value P7 of the frequency domain signal F7 is set to the maximum intensity of the audio. Not limited to this.

The following examples illustrate the practical application of the method for estimating the temperature of a blast furnace in the above embodiment of the present invention, but it is not limited thereto.

P，譬如n=1、2、...、P)之收發單元(例如在方位EN)發出一音波；接著，該處理器可進行上述測溫步驟，例如對任意收發單元(譬如其餘收發單元，也可為所有收發單元)依序進行下列測溫動作：對該音波取樣以各自取得一音波訊號；執行一帶通濾波；執行一快速傅立葉轉換取得一音頻最大強度；將該音頻最大強度各自乘上一方位比例常數後取得一閾值(如為一上昇峰值判定門檻值)；根據該上昇峰值判定門 檻值取得一上昇峰值時間，譬如尋找該目標訊號超出該閾值的至少一振幅峰值中的第一個振幅峰值發生的時間，作為該上昇峰值時間；比對各音波相對發聲源之上昇峰值時間取得各路徑音波飛行時間；根據各路徑音波飛行時間取得各路徑平均介質溫度，例如依據該上昇峰值時間、該發聲時間、一音波收發距離及一聲學常數計算一音波路徑溫度(如上式所示)；後續，該處理器可判斷該控制員是否對該測溫裝置發出一停止命令(例如按下一測溫停止按鈕)，若判斷為是，該測溫裝置停止測溫，若判斷為否，該處理器可驅使內部編號n+1(若n+1>P，則n+1=n+1-P)之收發單元做為另一發聲源，以發出另一音波作為測溫依據，再進行上述測溫動作。 For example, first, a controller can issue a temperature measurement command to the temperature measurement device (for example, press a temperature measurement start button); then, the controller can set a number of transceiver units (for example, set the number to P) , As an instruction stored in the memory, as a basis for the processor to execute; then, the processor can perform the above-mentioned sounding step T1, such as driving the internal number (for example, n, n

P, such as n=1, 2,..., P) the transceiver unit (for example in the position EN) emits a sound wave; then, the processor can perform the temperature measurement steps, for example, for any transceiver unit (such as other transceiver units) , It can also be used for all transceiver units) to perform the following temperature measurement operations in sequence: sample the sound wave to obtain a sound signal; perform a band-pass filter; perform a fast Fourier transform to obtain a maximum audio intensity; multiply the maximum audio intensity respectively A threshold is obtained after the upper azimuth proportional constant (for example, a rising peak determination threshold value); a rising peak time is obtained according to the rising peak determination threshold value, for example, searching for the first of at least one amplitude peak of the target signal exceeding the threshold The time when each amplitude peak occurs is used as the rising peak time; the sonic flight time of each path is obtained by comparing the rising peak time of each sound wave relative to the sound source; the average medium temperature of each path is obtained according to the sonic flight time of each path, for example, according to the rising peak time , The sounding time, a sound wave transmitting and receiving distance and an acoustic constant calculate a sound wave path temperature (as shown in the above formula); subsequently, the processor can determine whether the controller sends a stop command to the temperature measuring device (for example, press Temperature measurement stop button), if the judgment is yes, the temperature measurement device stops temperature measurement, if the judgment is no, the processor can drive the internal number n+1 (if n+1>P, then n+1=n+1 -P) The transceiver unit is used as another sound source, and another sound wave is used as the basis for temperature measurement, and then the above temperature measurement action is performed.

In this way, it can overcome the misjudgment of the temperature measurement results caused by the environmental noise generated by the airflow in the furnace, so that based on the sound wave transmission characteristics, the sound wave path temperature can be estimated by using the sound wave receiving and sending time in the blast furnace environment, which can improve the calculation of the flight time from the sound wave transmitter to the receiver. Accuracy, and can reduce the error caused by the average temperature estimation of the medium on the path of the sound wave, which is beneficial to meet the temperature measurement requirements of the blast furnace.

On the other hand, the present invention also provides a computer program product. After the computer program is loaded and executed, the computer can execute the method for estimating the temperature of a blast furnace as described above. For example: The computer program product may contain several program instructions, which can be implemented using existing programming languages, so as to execute the above-mentioned methods for estimating the temperature of a blast furnace, such as C or Python, but not Is limited.

On the other hand, the present invention also provides a computer-readable recording medium, such as an optical disc, a flash drive, or a hard disk, etc. The computer can read a program stored in the recording medium (such as the above-mentioned computer program), and when the computer loads the program After execution, the computer can complete the method of estimating the temperature of the blast furnace as described above.

The method for estimating the temperature of a blast furnace, the computer program product, and the computer-readable recording medium of the above-mentioned embodiment of the present invention use the time when the first peak amplitude exceeds the threshold as the rising peak time; and according to the rising peak time, The sounding time, the sound wave transmitting and receiving distance, and the acoustic constant are used to calculate the sound wave path temperature, which can overcome the environmental noise generated by the airflow in the furnace and cause the misjudgment of the temperature measurement result.

Compared with other technologies (such as using envelopes, dependency analysis, or peak-to-peak time difference of pulsed sound waves, etc.), the airflow inside the blast furnace causes environmental noise to affect the temperature measurement effect. The method for estimating the temperature of the blast furnace in the above embodiment of the present invention , Computer program products and computer-readable recording media can overcome the misjudgment of the temperature measurement results caused by the environmental noise generated by the airflow in the furnace, effectively calculate the temperature of the sound wave path, and improve the accuracy of the calculation of the flight time from the sound wave transmitter to the receiver. Reducing the error in the estimation of the average temperature of the medium on the path of the sound wave is conducive to meeting the temperature measurement requirements of the blast furnace.

In summary, the method, computer program product, and computer readable recording medium of the present invention for estimating the temperature of a blast furnace of the present invention convert the target signal into the frequency domain signal, and set the maximum amplitude peak of the frequency domain signal to The maximum intensity of the audio; the product of the maximum intensity of the audio and the azimuth proportional constant is calculated as the threshold; the time when the first amplitude peak of the target signal exceeds the threshold is found as the rising peak time; and according to the rise The peak time, the sounding time, the sound wave transmitting and receiving distance and the acoustic constant are used to calculate the sound wave path temperature, which can overcome the environmental noise generated by the airflow in the furnace, which leads to misjudgment of the temperature measurement result.

Thereby, the method for estimating the temperature of a blast furnace, the computer program product, and the computer-readable recording medium of the present invention can estimate the temperature of the sound wave path based on the sound wave transmission characteristics and use the sound wave receiving and sending time in the blast furnace environment, which can improve the sound wave transmitting end to the receiving end The accuracy of the flight time calculation can reduce the error caused by the average temperature estimation of the medium on the path of the sound wave, and improve the situation that other technologies affect the temperature measurement effect due to the environmental noise caused by the airflow inside the blast furnace, which is conducive to meeting the temperature measurement requirements of the blast furnace.

Although the present invention has been disclosed in preferred embodiments, it is not intended to limit the present invention. Anyone familiar with the art can make various changes and modifications without departing from the spirit and scope of the present invention. Therefore, the present invention The scope of protection shall be subject to the scope of the attached patent application.

T1: voice step

T2: temperature measurement step

Claims (7)

A method for estimating the temperature of a blast furnace is configured such that a processor is coupled to a memory and a plurality of transceiving units. The transceiving units are arranged equidistantly along a circular track in a horizontal plane in an environment to be measured. The processing The device executes the instructions stored in the memory, and the method includes the steps: Driving one of the several transceiver units as a sound source to emit a sound wave at a sounding time; and Perform at least one of the remaining transceiver units other than the sound source: Continuously sampling the sound wave to generate a sound wave signal as a target signal; Converting the target signal into a frequency domain signal, and setting a maximum amplitude peak of the frequency domain signal to a maximum audio intensity; Calculate a product of the maximum intensity of the audio frequency and the azimuth proportional constant as a threshold; Finding the time when the first amplitude peak of the at least one amplitude peak of the target signal exceeds the threshold as a rising peak time; and A sound path temperature is calculated according to the rising peak time, the sounding time, a sound wave transmitting and receiving distance, and an acoustic constant. The method for estimating the temperature of a blast furnace as described in claim 1, wherein the calculation method of the sonic path temperature is as follows:

Wherein, T is the sound wave path temperature; D is the sound wave transceiver distance, and the sound wave transceiver distance is the distance between the transceiver unit to which the target signal belongs and the transceiver unit that emits the sound wave;

Is the time difference of the rising peak time minus the sounding time; A is the acoustic constant. The method for estimating the temperature of a blast furnace as described in claim 1, wherein the azimuth proportional constant is less than 1, and the azimuth proportional constant has a positive correlation with the sound wave transmitting and receiving distance. The greater the sound wave transmitting and receiving distance, the azimuth proportional constant Bigger. The method for estimating the temperature of a blast furnace as described in claim 1, wherein the target signal is converted into the frequency domain signal after passing through a band-pass filtering process. The method for estimating the temperature of a blast furnace as described in claim 1, wherein the target signal is converted into the frequency domain signal according to a fast Fourier algorithm. A computer program product, when the computer program is loaded and executed, the computer can execute the method for estimating the temperature of a blast furnace as described in any one of claim items 1 to 5. A computer-readable recording medium. The computer can read a program stored in the recording medium. After the computer loads and executes the program, the computer can complete the estimation of the temperature of a blast furnace as described in any one of Claims 1 to 5 method.

Images