Blast furnace thermal and power data simulation method

Abstract

The application discloses a blast furnace thermal power and electric power data simulation method, and relates to the technical field of electric power systems. The method comprises the following steps: acquiring blast furnace electric power data and blast furnace thermal power data in the latest year, and calculating blast furnace electric power values of various sampling points; according to the blast furnace electric power values, classifying the blast furnace thermal power data, and calculating partial derivatives of the blast furnace thermal power data; calculating a blast furnace thermal power data mapping function of a simulation time according to data of the simulation time at historical time points, and calculating the blast furnace thermal power data of the simulation time according to the data of the simulation time at the historical time points and the partial derivatives of the blast furnace thermal power data. The method can simulate a blast furnace data acquisition system in a low-cost manner.

Description

Technical Field

[0001] This invention relates to power system technology, and in particular to a method for simulating blast furnace thermal and power data. Background Technology

[0002] The blast furnace environment is characterized by high temperature and high dust. Because the thermal data acquisition device of the blast furnace is in a harsh environment, the data it collects must be cleaned before it can be used.

[0003] The electrical data in blast furnace energy is characterized by large current, large high-order harmonics, and rapid amplitude changes. In actual operation, the high-order harmonic values ​​of the current change significantly every minute.

[0004] Existing blast furnace data simulations are all simulations of blast furnace electrical data, without simulation of blast furnace thermal data. Moreover, blast furnace electrical data is mostly generated using common power grid simulation data generation methods, without combining blast furnace thermal data for multi-faceted simulation, and therefore cannot reflect the actual operating conditions of the blast furnace. Summary of the Invention

[0005] In view of the deficiencies in the prior art, the technical problem to be solved by the present invention is to provide a blast furnace thermal and electrical data simulation method that can reflect the actual operating conditions of the blast furnace.

[0006] To address the aforementioned technical problems, the present invention provides a method for simulating blast furnace thermal and electrical data, characterized by the following specific steps:

[0007] 1) Obtain blast furnace power data and blast furnace thermal data from historical data of the most recent year, with a sampling point set at every m time interval;

[0008] There are 7 types of blast furnace power data: blast furnace voltage fundamental value U1, blast furnace voltage third harmonic value U3, blast furnace voltage fifth harmonic value U5, blast furnace current fundamental value I1, blast furnace current third harmonic value I3, blast furnace current fifth harmonic value I5, and blast furnace power factor value cosφ. U1, U3, U5, I1, I3, and I5 are all effective values.

[0009] There are 8 types of thermal data for blast furnaces: blast furnace air flow rate a1, blast furnace oxygen enrichment flow rate a2, blast furnace pulverized coal injection rate a3, blast furnace air pressure a4, blast furnace top pressure a5, blast furnace average top temperature over 1 minute a6, blast furnace average top temperature over 5 minutes a7, blast furnace average top temperature over 10 minutes a8.

[0010] 2) Based on U1, U3, U5, I1, I3, I5, and cosφ of each sampling point, calculate the blast furnace power value P at each sampling point. The specific calculation formula is as follows:

[0011]

[0012] 3) Select one type of blast furnace thermal data and define this type of blast furnace thermal data as the target thermal data ac;

[0013] 4) After sorting each sampling point in ascending order of the target thermal data ac, the first 50% of the sampling points are assigned to an initially empty sampling point set J1, and the remaining sampling points are assigned to an initially empty sampling point set J2.

[0014] 5) Calculate the partial derivatives of the target thermal data ac. The specific calculation formula is as follows:

[0015] F(X)`| ac =λ(P2-P1) / (ac.e-ac.b)

[0016] In the formula, F(X)`| ac Let λ be the partial derivative of the target thermal data ac, and λ be a pre-set constant value.

[0017] In the formula, P2 is the average value of the blast furnace power value of each sampling point in the sampling point set J2, P1 is the average value of the blast furnace power value of each sampling point in the sampling point set J1, ac.e is the average value of the target thermal data ac of each sampling point in the sampling point set J2, and ac.b is the average value of the target thermal data ac of each sampling point in the sampling point set J1.

[0018] 6) If the partial derivatives of all 8 types of blast furnace thermal data have been calculated, proceed to step 7); otherwise, select a type of blast furnace thermal data whose partial derivatives have not been calculated, define this type of blast furnace thermal data as the new target thermal data ac, and then proceed to step 4).

[0019] 7) Select a time point to be simulated and set that time point as time t;

[0020] 8) Generate a random perturbation s, where s is a random number with a value between 0.8 and 1.1;

[0021] The result of multiplying the blast furnace power data at time t (at the same time last year) with the random disturbance s is set as the blast furnace power simulation data at time t.

[0022] Generate eight random numbers with values ​​ranging from 0.8 to 1.1, and assign these eight random numbers to β1, β2, β3, β4, β5, β6, β7, and β8 respectively;

[0023] Set the blast furnace thermal data at time t, which was at the same time last year, as the blast furnace thermal simulation data at time t.

[0024] 9) Based on U1, U3, U5, I1, I3, I5, and cosφ in the blast furnace power simulation data at time t, calculate the blast furnace thermal data mapping function F(X) at time t. The specific calculation formula is as follows:

[0025]

[0026] 10) Based on the blast furnace thermal simulation data at time t, calculate the blast furnace thermal data at time t. The specific calculation formula is as follows:

[0027] a1.t=a1.q+(β1-1)F(X) / F(X)`| a1

[0028] a2.t=a2.q+(β2-1)F(X) / F(X)`| a2

[0029] a3.t=a3.q+(β3-1)F(X) / F(X)`| a3

[0030] a4.t=a4.q+(β4-1)F(X) / F(X)`| a4

[0031] a5.t=a5.q+(β5-1)F(X) / F(X)`| a5

[0032] a6.t=a6.q+(β6-1)F(X) / F(X)`| a6

[0033] a7.t=a7.q+(β7-1)F(X) / F(X)`| a7

[0034] a8.t=a8.q+(β8-1)F(X) / F(X)`| a8

[0035] In the formula, a1.t is the value of the blast furnace blast flow rate a1 at time t, a1.q is the simulated value of the blast furnace blast flow rate a1 at time t, and F(X)`| a1 Let F(X) be the partial derivative of the blast furnace blast flow rate a1, a2.t be the value of the blast furnace oxygen-enriched flow rate a2 at time t, and a2.q be the simulated value of the blast furnace oxygen-enriched flow rate a2 at time t. a2 Let a2 be the partial derivative of the oxygen-enriched flow rate in the blast furnace, a3.t be the value of the pulverized coal injection rate a3 in the blast furnace at time t, and a3.q be the simulated value of the pulverized coal injection rate a3 in the blast furnace at time t. F(X)`| a3Let a3 be the partial derivative of the blast furnace pulverized coal injection rate, a4.t be the value of the blast furnace blast pressure a4 at time t, and a4.q be the simulated value of the blast furnace blast pressure a4 at time t. F(X)`| a4 Let a be the partial derivative of the blast furnace blast pressure a4, a5.t be the value of the blast furnace top pressure a5 at time t, a5.q be the simulated value of the blast furnace top pressure a5 at time t, and F(X)`| a5 Let a5 be the partial derivative of the blast furnace top pressure, a6.t be the value of the blast furnace average top temperature a6 over 1 minute at time t, and a6.q be the simulated value of the blast furnace average top temperature a6 over 1 minute at time t. F(X)`| a6 Let a6 be the partial derivative of the blast furnace average top temperature over 1 minute, a7.t be the value of the blast furnace average top temperature over 5 minutes a7 at time t, and a7.q be the value of the blast furnace average top temperature over 5 minutes a7 at time t. F(X)`| a7 Let a7 be the partial derivative of the 5-minute average top furnace temperature of the blast furnace, a8.t be the value of the 10-minute average top furnace temperature of the blast furnace a8 at time t, and a8.q be the simulated value of the 10-minute average top furnace temperature of the blast furnace a8 at time t. F(X)`| a8 Let a8 be the partial derivative of the average top furnace temperature a8 over 10 minutes in the blast furnace;

[0036] 11) If all simulation times have been completed, end the simulation; otherwise, select the next time to be simulated, set that time as time t, and then go to step 8.

[0037] The blast furnace thermal and electrical data simulation method provided by this invention combines power and thermal data at historical moments using a linkage method to calculate the blast furnace thermal data mapping function and blast furnace thermal simulation data at the simulation moment. This allows for the simulation of the blast furnace data acquisition system at low cost and reflects the actual operating conditions of the blast furnace. Detailed Implementation

[0038] The technical solution of the present invention will be further described in detail below with reference to specific embodiments. However, these embodiments are not intended to limit the present invention. Any similar structures and similar variations of the present invention should be included in the protection scope of the present invention. The commas in the present invention all indicate the relationship between and. The English letters in the present invention are case-sensitive.

[0039] The blast furnace thermal and electrical data simulation method provided in this embodiment of the invention is characterized by the following specific steps:

[0040] 1) From the historical data of the most recent year, sampling points are set at intervals of m to obtain the blast furnace power data and blast furnace thermal data of each sampling point. The blast furnace power data and blast furnace thermal data can be obtained from the existing blast furnace energy management system. In this embodiment, the m time period is set to 1 minute. In other embodiments, the m time period can also be set to 2 minutes or other time periods.

[0041] There are 7 types of blast furnace power data: blast furnace voltage fundamental value U1, blast furnace voltage third harmonic value U3, blast furnace voltage fifth harmonic value U5, blast furnace current fundamental value I1, blast furnace current third harmonic value I3, blast furnace current fifth harmonic value I5, and blast furnace power factor value cosφ. U1, U3, U5, I1, I3, and I5 are all effective values.

[0042] There are 8 types of thermal data for blast furnaces: blast furnace air flow rate a1, blast furnace oxygen enrichment flow rate a2, blast furnace pulverized coal injection rate a3, blast furnace air pressure a4, blast furnace top pressure a5, blast furnace average top temperature over 1 minute a6, blast furnace average top temperature over 5 minutes a7, blast furnace average top temperature over 10 minutes a8.

[0043] 2) Based on U1, U3, U5, I1, I3, I5, and cosφ of each sampling point, calculate the blast furnace power value P at each sampling point. The specific calculation formula is as follows:

[0044]

[0045] 3) Select a type of blast furnace thermal data (e.g., blast furnace air flow rate a1) and define this type of blast furnace thermal data as the target thermal data ac;

[0046] 4) After sorting each sampling point in ascending order of the target thermal data ac, the first 50% of the sampling points are assigned to an initially empty sampling point set J1, and the remaining sampling points are assigned to an initially empty sampling point set J2.

[0047] 5) Calculate the partial derivatives of the target thermal data ac. The specific calculation formula is as follows:

[0048] F(X)`| ac =λ(P2-P1) / (ac.e-ac.b)

[0049] In the formula, F(X)`| ac Let λ be the partial derivative of the target thermal data ac, and let λ be a pre-set constant value, typically 0.8.

[0050] In the formula, P2 is the average value of the blast furnace power value of each sampling point in the sampling point set J2, P1 is the average value of the blast furnace power value of each sampling point in the sampling point set J1, ac.e is the average value of the target thermal data ac of each sampling point in the sampling point set J2, and ac.b is the average value of the target thermal data ac of each sampling point in the sampling point set J1.

[0051] 6) If the partial derivatives of all 8 types of blast furnace thermal data have been calculated, proceed to step 7); otherwise, select a type of blast furnace thermal data whose partial derivatives have not been calculated (e.g., blast furnace oxygen enrichment flow rate a2), define this type of blast furnace thermal data as the new target thermal data ac, and then proceed to step 4).

[0052] 7) Select a time point to be simulated and set that time point as time t;

[0053] 8) Generate a random perturbation s, where s is a random number with a value between 0.8 and 1.1;

[0054] The result of multiplying the blast furnace power data at time t (at the same time last year) with the random disturbance s is set as the blast furnace power simulation data at time t.

[0055] Generate eight random numbers with values ​​ranging from 0.8 to 1.1, and assign these eight random numbers to β1, β2, β3, β4, β5, β6, β7, and β8 respectively;

[0056] Set the blast furnace thermal data at time t, which was at the same time last year, as the blast furnace thermal simulation data at time t.

[0057] 9) Based on U1, U3, U5, I1, I3, I5, and cosφ in the blast furnace power simulation data at time t, calculate the blast furnace thermal data mapping function F(X) at time t. The specific calculation formula is as follows:

[0058]

[0059] 10) Based on the blast furnace thermal simulation data at time t, calculate the blast furnace thermal data at time t. The specific calculation formula is as follows:

[0060] a1.t=a1.q+(β1-1)F(X) / F(X)`| a1

[0061] a2.t=a2.q+(β2-1)F(X) / F(X)`| a2

[0062] a3.t=a3.q+(β3-1)F(X) / F(X)`| a3

[0063] a4.t=a4.q+(β4-1)F(X) / F(X)`| a4

[0064] a5.t=a5.q+(β5-1)F(X) / F(X)`| a5

[0065] a6.t=a6.q+(β6-1)F(X) / F(X)`| a6

[0066] a7.t=a7.q+(β7-1)F(X) / F(X)`| a7

[0067] a8.t=a8.q+(β8-1)F(X) / F(X)`| a8

[0068] In the formula, a1.t is the value of the blast furnace blast flow rate a1 at time t, a1.q is the simulated value of the blast furnace blast flow rate a1 at time t, and F(X)`| a1 Let F(X) be the partial derivative of the blast furnace blast flow rate a1, a2.t be the value of the blast furnace oxygen-enriched flow rate a2 at time t, and a2.q be the simulated value of the blast furnace oxygen-enriched flow rate a2 at time t. a2 Let a2 be the partial derivative of the oxygen-enriched flow rate in the blast furnace, a3.t be the value of the pulverized coal injection rate a3 in the blast furnace at time t, and a3.q be the simulated value of the pulverized coal injection rate a3 in the blast furnace at time t. F(X)`| a3 Let a3 be the partial derivative of the blast furnace pulverized coal injection rate, a4.t be the value of the blast furnace blast pressure a4 at time t, and a4.q be the simulated value of the blast furnace blast pressure a4 at time t. F(X)`| a4 Let a be the partial derivative of the blast furnace blast pressure a4, a5.t be the value of the blast furnace top pressure a5 at time t, a5.q be the simulated value of the blast furnace top pressure a5 at time t, and F(X)`| a5 Let a5 be the partial derivative of the blast furnace top pressure, a6.t be the value of the blast furnace average top temperature a6 over 1 minute at time t, and a6.q be the simulated value of the blast furnace average top temperature a6 over 1 minute at time t. F(X)`| a6 Let a6 be the partial derivative of the blast furnace average top temperature over 1 minute, a7.t be the value of the blast furnace average top temperature over 5 minutes a7 at time t, and a7.q be the value of the blast furnace average top temperature over 5 minutes a7 at time t. F(X)`| a7 Let a7 be the partial derivative of the 5-minute average top furnace temperature of the blast furnace, a8.t be the value of the 10-minute average top furnace temperature of the blast furnace a8 at time t, and a8.q be the simulated value of the 10-minute average top furnace temperature of the blast furnace a8 at time t. F(X)`| a8 Let a8 be the partial derivative of the average top furnace temperature a8 over 10 minutes in the blast furnace;

[0069] 11) If all simulation times have been completed, end the simulation; otherwise, select the next time to be simulated, set that time as time t, and then go to step 8.

Claims

1. A method for simulating the thermal and electrical data of a blast furnace, characterized in that, The specific steps are as follows: 1) Obtain blast furnace power data and blast furnace thermal data from historical data of the most recent year, with a sampling point set at every m time interval; There are seven types of blast furnace power data: blast furnace voltage fundamental value U1, blast furnace voltage third harmonic value U3, blast furnace voltage fifth harmonic value U5, blast furnace current fundamental value I1, blast furnace current third harmonic value I3, blast furnace current fifth harmonic value I5, and blast furnace power factor. U1, U3, U5, I1, I3, and I5 are all valid values; There are 8 types of thermal data for blast furnaces: blast furnace air flow rate a1, blast furnace oxygen enrichment flow rate a2, blast furnace pulverized coal injection rate a3, blast furnace air pressure a4, blast furnace top pressure a5, blast furnace average top temperature over 1 minute a6, blast furnace average top temperature over 5 minutes a7, blast furnace average top temperature over 10 minutes a8. 2) Based on the U1, U3, U5, I1, I3, I5, ... of each sampling point The blast furnace power value P at each sampling point is calculated using the following formula: 3) Select one type of blast furnace thermal data and define this type of blast furnace thermal data as the target thermal data ac; 4) After sorting each sampling point in ascending order of the target thermal data ac, the first 50% of the sampling points are assigned to an initially empty sampling point set J1, and the remaining sampling points are assigned to an initially empty sampling point set J2. 5) Calculate the partial derivatives of the target thermal data ac. The specific calculation formula is as follows: F(X)`| ac =λ(P2-P1) / (ac.e-ac.b) In the formula, F(X)`| ac Let λ be the partial derivative of the target thermal data ac, and λ be a pre-set constant value. In the formula, P2 is the average value of the blast furnace power value of each sampling point in the sampling point set J2, P1 is the average value of the blast furnace power value of each sampling point in the sampling point set J1, ac.e is the average value of the target thermal data ac of each sampling point in the sampling point set J2, and ac.b is the average value of the target thermal data ac of each sampling point in the sampling point set J1. 6) If the partial derivatives of all 8 types of blast furnace thermal data have been calculated, proceed to step 7); otherwise, select a type of blast furnace thermal data whose partial derivatives have not been calculated, define this type of blast furnace thermal data as the new target thermal data ac, and then proceed to step 4). 7) Select a time point to be simulated and set that time point as time t; 8) Generate a random perturbation s, where s is a random number with a value ranging from 0.8 to 1.1; The result of multiplying the blast furnace power data at time t at the same time last year with the random disturbance s is set as the blast furnace power simulation data at time t; Generate eight random numbers with values ​​ranging from 0.8 to 1.1, and assign these eight random numbers to β1, β2, β3, β4, β5, β6, β7, and β8 respectively; Set the blast furnace thermal data at time t, which was at the same time last year, as the blast furnace thermal simulation data at time t. 9) Based on the blast furnace power simulation data at time t, U1, U3, U5, I1, I3, I5, The specific formula for calculating the blast furnace thermal data mapping function F(X) at time t is as follows: 10) Based on the blast furnace thermal simulation data at time t, calculate the blast furnace thermal data at time t. The specific calculation formula is as follows: a1.t=a1.q+(β1-1)F(X) / F(X)`| a1 a2.t=a2.q+(β2-1)F(X) / F(X)`| a2 a3.t=a3.q+(β3-1)F(X) / F(X)`| a3 a4.t=a4.q+(β4-1)F(X) / F(X)`| a4 a5.t=a5.q+(β5-1)F(X) / F(X)`| a5 a6.t=a6.q+(β6-1)F(X) / F(X)`| a6 a7.t=a7.q+(β7-1)F(X) / F(X)`| a7 a8.t=a8.q+(β8-1)F(X) / F(X)`| a8 In the formula, a1.t is the value of the blast furnace blast flow rate a1 at time t, a1.q is the simulated value of the blast furnace blast flow rate a1 at time t, and F(X)`| a1 Let F(X) be the partial derivative of the blast furnace blast flow rate a1, a2.t be the value of the blast furnace oxygen-enriched flow rate a2 at time t, and a2.q be the simulated value of the blast furnace oxygen-enriched flow rate a2 at time t. a2 Let a2 be the partial derivative of the oxygen-enriched flow rate in the blast furnace, a3.t be the value of the pulverized coal injection rate a3 in the blast furnace at time t, and a3.q be the simulated value of the pulverized coal injection rate a3 in the blast furnace at time t. F(X)`| a3 Let a3 be the partial derivative of the blast furnace pulverized coal injection rate, a4.t be the value of the blast furnace blast pressure a4 at time t, and a4.q be the simulated value of the blast furnace blast pressure a4 at time t. F(X)`| a4 Let F(X) be the partial derivative of the blast furnace blast pressure a4, a5.t be the value of the blast furnace top pressure a5 at time t, a5.q be the simulated value of the blast furnace top pressure a5 at time t, and F(X)`| a5 Let F(X) be the partial derivative of the blast furnace top pressure a5, a6.t be the value of the blast furnace average top temperature a6 over 1 minute at time t, a6.q be the simulated value of the blast furnace average top temperature a6 over 1 minute at time t, and F(X)`| a6 Let a6 be the partial derivative of the blast furnace average top temperature a6 over 1 minute, a7.t be the value of the blast furnace average top temperature a7 over 5 minutes at time t, and a7.q be the value of the blast furnace average top temperature a7 over 5 minutes at time t. F(X)`| a7 Let a7 be the partial derivative of the 5-minute average top furnace temperature of the blast furnace, a8.t be the value of the 10-minute average top furnace temperature of the blast furnace a8 at time t, and a8.q be the simulated value of the 10-minute average top furnace temperature of the blast furnace a8 at time t. F(X)`| a8 Let a8 be the partial derivative of the average top furnace temperature a8 over 10 minutes in the blast furnace; 11) If all simulation times have been completed, end the simulation; otherwise, select the next time to be simulated, set that time as time t, and then go to step 8.

Images