CN101463407B - Method for calculating converter steel melting lime adding amount - Google Patents
Method for calculating converter steel melting lime adding amount Download PDFInfo
- Publication number
- CN101463407B CN101463407B CN2008102290126A CN200810229012A CN101463407B CN 101463407 B CN101463407 B CN 101463407B CN 2008102290126 A CN2008102290126 A CN 2008102290126A CN 200810229012 A CN200810229012 A CN 200810229012A CN 101463407 B CN101463407 B CN 101463407B
- Authority
- CN
- China
- Prior art keywords
- lime
- basicity
- input
- heat
- adding amount
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 239000004571 lime Substances 0.000 title claims abstract description 108
- 235000008733 Citrus aurantifolia Nutrition 0.000 title claims abstract description 105
- 235000011941 Tilia x europaea Nutrition 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 44
- 229910000831 Steel Inorganic materials 0.000 title claims description 39
- 239000010959 steel Substances 0.000 title claims description 39
- 238000002844 melting Methods 0.000 title claims description 7
- 230000008018 melting Effects 0.000 title claims description 7
- 239000002893 slag Substances 0.000 claims abstract description 56
- 238000012706 support-vector machine Methods 0.000 claims abstract description 50
- 239000010459 dolomite Substances 0.000 claims abstract description 23
- 229910000514 dolomite Inorganic materials 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 230000000694 effects Effects 0.000 claims abstract description 12
- 101100399296 Mus musculus Lime1 gene Proteins 0.000 claims description 98
- 239000000292 calcium oxide Substances 0.000 claims description 41
- ODINCKMPIJJUCX-UHFFFAOYSA-N calcium oxide Inorganic materials [Ca]=O ODINCKMPIJJUCX-UHFFFAOYSA-N 0.000 claims description 41
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 29
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 26
- 239000000377 silicon dioxide Substances 0.000 claims description 16
- 229910052742 iron Inorganic materials 0.000 claims description 13
- 230000036284 oxygen consumption Effects 0.000 claims description 13
- 229910004298 SiO 2 Inorganic materials 0.000 claims description 12
- 238000007664 blowing Methods 0.000 claims description 12
- XWHPIFXRKKHEKR-UHFFFAOYSA-N iron silicon Chemical compound [Si].[Fe] XWHPIFXRKKHEKR-UHFFFAOYSA-N 0.000 claims description 12
- 238000010606 normalization Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 11
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 claims description 10
- 230000003647 oxidation Effects 0.000 claims description 10
- 238000007254 oxidation reaction Methods 0.000 claims description 10
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 9
- 229910052799 carbon Inorganic materials 0.000 claims description 9
- DALUDRGQOYMVLD-UHFFFAOYSA-N iron manganese Chemical compound [Mn].[Fe] DALUDRGQOYMVLD-UHFFFAOYSA-N 0.000 claims description 9
- DPTATFGPDCLUTF-UHFFFAOYSA-N phosphanylidyneiron Chemical compound [Fe]#P DPTATFGPDCLUTF-UHFFFAOYSA-N 0.000 claims description 9
- 229910052717 sulfur Inorganic materials 0.000 claims description 9
- 239000011593 sulfur Substances 0.000 claims description 9
- BRPQOXSCLDDYGP-UHFFFAOYSA-N calcium oxide Chemical compound [O-2].[Ca+2] BRPQOXSCLDDYGP-UHFFFAOYSA-N 0.000 claims description 8
- 239000000203 mixture Substances 0.000 claims description 7
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 6
- 229910052760 oxygen Inorganic materials 0.000 claims description 6
- 239000001301 oxygen Substances 0.000 claims description 6
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 claims description 5
- 239000011575 calcium Substances 0.000 claims description 5
- 229910052791 calcium Inorganic materials 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 229910052710 silicon Inorganic materials 0.000 claims description 5
- 239000010703 silicon Substances 0.000 claims description 5
- 101100373011 Drosophila melanogaster wapl gene Proteins 0.000 claims description 3
- 238000012821 model calculation Methods 0.000 claims description 3
- 210000004483 pasc Anatomy 0.000 claims description 3
- -1 bath temperature Substances 0.000 claims description 2
- 238000009628 steelmaking Methods 0.000 abstract description 9
- 230000003068 static effect Effects 0.000 abstract description 5
- 238000004458 analytical method Methods 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000010276 construction Methods 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 238000004090 dissolution Methods 0.000 description 6
- 239000000463 material Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000012535 impurity Substances 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 235000019739 Dicalciumphosphate Nutrition 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 2
- 238000003723 Smelting Methods 0.000 description 2
- 239000001506 calcium phosphate Substances 0.000 description 2
- NEFBYIFKOOEVPA-UHFFFAOYSA-K dicalcium phosphate Chemical compound [Ca+2].[Ca+2].[O-]P([O-])([O-])=O NEFBYIFKOOEVPA-UHFFFAOYSA-K 0.000 description 2
- 229940038472 dicalcium phosphate Drugs 0.000 description 2
- 229910000390 dicalcium phosphate Inorganic materials 0.000 description 2
- AMWRITDGCCNYAT-UHFFFAOYSA-L hydroxy(oxo)manganese;manganese Chemical compound [Mn].O[Mn]=O.O[Mn]=O AMWRITDGCCNYAT-UHFFFAOYSA-L 0.000 description 2
- 229910052698 phosphorus Inorganic materials 0.000 description 2
- 239000011574 phosphorus Substances 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 108010074864 Factor XI Proteins 0.000 description 1
- 239000005864 Sulphur Substances 0.000 description 1
- QCJQWJKKTGJDCM-UHFFFAOYSA-N [P].[S] Chemical compound [P].[S] QCJQWJKKTGJDCM-UHFFFAOYSA-N 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000036571 hydration Effects 0.000 description 1
- 238000006703 hydration reaction Methods 0.000 description 1
- UQSXHKLRYXJYBZ-UHFFFAOYSA-N iron oxide Inorganic materials [Fe]=O UQSXHKLRYXJYBZ-UHFFFAOYSA-N 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000000395 magnesium oxide Substances 0.000 description 1
- CPLXHLVBOLITMK-UHFFFAOYSA-N magnesium oxide Inorganic materials [Mg]=O CPLXHLVBOLITMK-UHFFFAOYSA-N 0.000 description 1
- AXZKOIWUVFPNLO-UHFFFAOYSA-N magnesium;oxygen(2-) Chemical compound [O-2].[Mg+2] AXZKOIWUVFPNLO-UHFFFAOYSA-N 0.000 description 1
- WPBNNNQJVZRUHP-UHFFFAOYSA-L manganese(2+);methyl n-[[2-(methoxycarbonylcarbamothioylamino)phenyl]carbamothioyl]carbamate;n-[2-(sulfidocarbothioylamino)ethyl]carbamodithioate Chemical compound [Mn+2].[S-]C(=S)NCCNC([S-])=S.COC(=O)NC(=S)NC1=CC=CC=C1NC(=S)NC(=O)OC WPBNNNQJVZRUHP-UHFFFAOYSA-L 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- QMQXDJATSGGYDR-UHFFFAOYSA-N methylidyneiron Chemical compound [C].[Fe] QMQXDJATSGGYDR-UHFFFAOYSA-N 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- NDLPOXTZKUMGOV-UHFFFAOYSA-N oxo(oxoferriooxy)iron hydrate Chemical compound O.O=[Fe]O[Fe]=O NDLPOXTZKUMGOV-UHFFFAOYSA-N 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
Landscapes
- Carbon Steel Or Casting Steel Manufacturing (AREA)
Abstract
The invention belongs to the automatic control technical field and relates to the construction of a converter steelmaking static model, in particular to a calculating method of the adding amount of lime during a converter steelmaking process. In the aspect of constructing an alkalinity deviation estimation model, a proper variable is selected as the input of the alkalinity deviation estimation model through the analysis of reasons causing deviations between the input alkalinity and the slag terminal alkalinity. Then, qualified heat data in a historical database are used for constructing the alkalinity deviation estimation model of a support vector machine and forecasting the deviation between the input alkalinity and the slag terminal alkalinity of the current heat. In the aspect of the calculation of the lime adding amount, a predictive value of the alkalinity deviation estimation model is used for rectifying alkalinity parameters in an empirical formula and eliminating effects of the lime adding amount on the adding amount of dolomite. Finally, a calculation formula of the lime adding amount is obtained. The calculating method of the adding amount of lime during the converter steelmaking process has the beneficial effects of effectively improving the calculation accuracy of the lime adding amount and simultaneously ensuring the slag terminal alkalinity to satisfy technological requirements.
Description
Technical field
The invention belongs to technical field of automatic control, relate to the foundation that static model are produced in converter steelmaking, the method for calculation of lime adding amount in particularly a kind of converter steelmaking production process.
Background technology
It is with impurity contents such as molten iron carbon drop, intensification, reduction phosphorus sulphur that converter steelmaking is produced, and obtains the commercial run of qualified molten steel., add slag making materials impurity is removed from molten steel the impurity element oxidation in the molten iron by top blast oxygen.Slagging regime is one of process system important during converter steelmaking is produced, and the basicity of slag is to weigh the important indicator of slag making quality, also be influence phosphorus, sulphur removes one of principal element of ratio.Basicity is normally defined CaO content and SiO in the slag
2The ratio of content: R=w (CaO)/w (SiO
2).Control suitable lime adding amount, guarantee that the slag end-point alkalinity is most important for producing qualified molten steel.
The lime adding amount computation model belongs to the part of converter static model.At present, the modeling method of calculating at lime adding amount mainly contains:
Experimental formula method based on material balance, the experimental formula method with silicone content in the molten iron and process goal basicity as the reference (Chen Zhongwei that determines lime adding amount, Yuan keeps the static mathematical model and the realization [J] of modest .LD converter smelting. steel-making, 2000,16 (5): 31-34);
The increment Return Law based on statistical theory, the increment homing method is the lime adding amount (Zhu Guangjun of the current heat of basic calculation with the increment size between current heat molten iron condition and the calibration furnace time, this river of beam. Converter static control is optimized model [J]. steel-making, 1999,15 (4): 25-28).
Yet these two kinds of lime adding amount method of calculation can't guarantee well that the slag end-point alkalinity satisfies processing requirement, sometimes even make the basicity of slag fluctuation between heat very big.
Summary of the invention
The technical problem to be solved in the present invention provides a kind of method for calculating converter steel melting lime adding amount, this method is by the deviation between basicity estimation of deviation model prediction input basicity and the slag end-point alkalinity, revise the basicity parameter in the lime adding amount empirical Calculation formula, satisfy processing requirement to guarantee the slag end-point alkalinity.
Technical scheme of the present invention is:
Analysis is bessemerized influences lime dissolved factor in the process, select suitable input variable according to analytical results, and set up SVMs basicity estimation of deviation model, forecasts the deviation between current heat input basicity and the terminal point basicity of slag.Utilize the basicity value in the estimated value correction experimental formula of basicity deviation afterwards, the lime adding amount computation model that is improved.
The step of SVMs basicity estimation of deviation modelling is as follows:
Step 1: the principal element of analyzing influence slag end-point alkalinity comprises slag composition, bath temperature, lime quality and blowing oxygen quantity.
Step 2: determine the input variable of SVMs basicity estimation of deviation model, comprise molten iron silicon content, molten iron manganese content, molten iron phosphorus content, molten steel sulfur content, rhombspar add-on, Molten Steel End Point, lime activity and the total oxygen-consumption of heat.
(1) data logging that selection has the terminal point slag composition to chemically examine in historical data base;
(2) data logging of selecting endpoint carbon content and temperature to hit simultaneously;
(3) select molten iron silicon content between 0.2~0.6%, and the data logging of slag end-point alkalinity between 2.8~3.5, standard heat storehouse set up.
(1) the input basicity R of each heat in the base of calculation heat storehouse
Input:
Wherein, CaO
LimeThe calcium oxide content of bringing into for lime; CaO
DolomiteThe calcium oxide content of bringing into for rhombspar; SiO
2limeThe silica volume of bringing into for lime; SiO
2dolomiteBe the silica volume of bringing into by rhombspar; K=2.14 is SiO
2The ratio of the molecular mass of/Si; K * w[Si] * W
HMSilica volume for pasc reaction generation in the molten iron; SiO
2scrapBe silica volume of bringing into by steel scrap and the silica volume that reacts generation thereof;
(2) calculate deviation delta R between input basicity and the slag end-point alkalinity, as the output of SVMs basicity estimation of deviation model:
ΔR=R
input-R
aim
Wherein, R
InputBe input basicity; R
AimBe target basicity;
(3) normalization method input/output variable, and the kernel function and the loss function of selection SVMs;
(4) use the crosscheck method to determine the parameter of SVMs basicity estimation of deviation model.
The step that current heat lime adding amount is calculated is as follows:
(1) search and this heat molten iron silicon content, molten iron manganese content, molten iron phosphorus content, molten steel sulfur content, rhombspar add-on, lime activity and the most close heat of terminal point target temperature in standard heat storehouse, with the reference of total oxygen-consumption of this heat as the total oxygen-consumption of current heat, use the process goal temperature of current heat of molten steel to replace Molten Steel End Point simultaneously, as the input of this heat SVMs basicity estimation of deviation model;
(2) the current heat molten iron silicon content that will determine, molten iron manganese content, molten iron phosphorus content, molten steel sulfur content, rhombspar add-on, Molten Steel End Point, lime activity and the total oxygen-consumption normalization method of heat are brought normalized each input variable value into SVMs basicity estimation of deviation Model Calculation;
(3) with the anti-normalization method of calculation result of SVMs basicity estimation of deviation model, obtain predicted value to current heat basicity deviation.
Step 2: use the basicity deviate of forecast that the basicity parameter in the experimental formula is revised, and calculate lime adding amount:
Wherein, K=2.14 is SiO
2The ratio of the molecular mass of/Si; W[Si] be the massfraction of silicon in the molten iron; W
HMBe molten steel quality; R
AimBe target basicity; Δ R is the estimated value of basicity deviation; W (CaO
Lime_eff) be the calcareous amount mark of the efficient oxidation in the lime; W
DolomiteBe the rhombspar add-on; W (CaO
Dolomite_eff) be the massfraction of the efficient oxidation calcium in the rhombspar.
After the current heat finishing blowing, if carry out slag chemical examination, and basicity of slag and endpoint carbon content and temperature satisfy processing requirement, then this heat information added standard heat storehouse, uses for setting up SVMs basicity estimation of deviation model.
The gordian technique of aforesaid method is to have set up SVMs basicity estimation of deviation model, forecasts the deviation between current heat input basicity and the slag end-point alkalinity, and basicity deviate is according to weather report adjusted lime adding amount.
Effect of the present invention and benefit are to remedy deviation between input basicity and the end-point alkalinity preferably, improve the accuracy that lime adding amount calculates, make the slag end-point alkalinity more near the target basicity value of processing requirement, help controlling the stable in properties of slag between heat.
Description of drawings
Fig. 1 is a method of calculation synoptic diagram of the present invention.
Wherein, R
AimBe terminal point target basicity; x
1, x
2... x
nInput for SVMs basicity estimation of deviation model; Δ R is input basicity and slag end-point alkalinity predicted value; R
InputBe input basicity; W
LimeBe lime adding amount; R
EndBe the slag end-point alkalinity; ε
1Actual value for input basicity and slag end-point alkalinity deviation; ε
2Poor for basicity deviation predicted value and actual value.
Fig. 2 is a SVMs basicity estimation of deviation model structure synoptic diagram.
Wherein, x
1Be molten iron silicon content; x
2Be molten iron manganese content; x
3Be the molten iron phosphorus content; x
4Be molten steel sulfur content; x
5Be the rhombspar add-on; x
6Be terminal temperature; x
7Be lime activity; x
8Be the stove oxygen-consumption; Y is the basicity deviation; S is final support vector number.
The extensive error curve diagram of test samples of correspondence when Fig. 3 gets different value for parameter beta and C.
Fig. 4 is the predict the outcome figure of SVMs basicity estimation of deviation model to test sample book.
Fig. 5 is the synoptic diagram that concerns between the calculated value of test sample book lime adding amount and the actual value.
Wherein, long and short dash line represents that calculated value equates with actual value, and two interior points of solid line all satisfy Error Absolute Value in 1 ton scope, and two interior points of dotted line all satisfy Error Absolute Value in 0.7 ton scope.
Embodiment
Be described in detail specific embodiments of the invention below in conjunction with technical scheme and accompanying drawing 1.
Analyze the influencing factor that produces the basicity deviation, determine the input of SVMs basicity estimation of deviation model.
In actual production, deviation is that incomplete dissolving by lime causes between input basicity and the blow end point basicity of slag.At the blowing initial stage, the lime surface that adds converter generates fusing point up to 2130 ℃ and fine and close hard Dicalcium Phosphate (Feed Grade) layer, hinders the dissolving of lime.The dissolving of lime continues whole converting process, in addition during finishing blowing lime still in dissolving.By the principal element of analyzing influence lime dissolution rate, determine the input of SVMs basicity estimation of deviation model:
A. the composition of slag: it has a significant impact the lime dissolution rate, and the content of compositions such as the calcium oxide in the slag, magnesium oxide, manganese oxide and ferric oxide all affects the dissolution rate of lime.Slag is mainly by the auxiliary material dissolving of the oxidation of elements such as the silicon in the molten iron, phosphorus, manganese, iron and adding and generate.
B. bath temperature: the bath temperature height helps the reduction of slag viscosity, quickens slag and permeates in lime block, impels the rapid fusion of the Dicalcium Phosphate (Feed Grade) formation slag that comes off, thereby improves the dissolution rate of lime.
C. the quality of lime: surface porosity, void content height, the quickened lime that response capacity is strong, help slag and enter in the lime block, enlarge reaction area, quicken the dissolution process of lime.The activity of lime is the physical quantity that characterizes the lime hydration speed of response, also is an important parameter that influences the lime dissolution rate.
D. blowing oxygen quantity: the dissolving of lime can continue whole process, and the carbon drop amount that how much depends primarily on of blowing oxygen quantity and the requirement of intensification amount increase the abundant dissolving that blowing oxygen quantity helps lime, improves basicity of slag.
To sum up, eight input variables determining SVMs basicity estimation of deviation model are respectively: molten iron silicon content, molten iron manganese content, molten iron phosphorus content, molten steel sulfur content, rhombspar add-on, Molten Steel End Point, lime activity and the total oxygen-consumption of heat.
After the input of determining SVMs basicity estimation of deviation model, need to select suitable sample to carry out modeling.In actual production, the slag of not all heat all can be chemically examined, and therefore selects to have the data of terminal point slag composition chemical examination; Endpoint carbon content and temperature are to weigh the leading indicator of bessemerizing quality, therefore select endpoint carbon content and temperature to hit the data logging of heat simultaneously; For guaranteeing not carry out the deslagging operation in the converting process, select the data logging of molten iron silicon content between 0.2~0.6%; According to the requirement of on-the-spot technology, the slag end-point alkalinity preferably is controlled between 2.8~3.5, thereby selects the data logging of slag end-point alkalinity between 2.8~3.5.So selecting data is in order to guarantee the accuracy of SVMs basicity estimation of deviation model.
Use the above-mentioned data of selecting to set up SVMs basicity estimation of deviation model, model is output as the difference of input basicity and slag end-point alkalinity.Input basicity R
InputBe defined as: contain in the calcium oxide content that contains in all auxiliary materials of the converter of packing into and the main auxiliary material and the ratio of the silica volume that may generate.Calculation formula is
Wherein:
CaO
LimeThe calcium oxide content of bringing into for lime;
CaO
DolomiteThe calcium oxide content of bringing into for rhombspar;
SiO
2limeThe silica volume of bringing into for lime;
SiO
2dolomiteBe the silica volume of bringing into by rhombspar;
K=2.14 is SiO
2The ratio of the molecular mass of/Si;
K * w[Si] * W
HMSilica volume for pasc reaction generation in the molten iron;
SiO
2scrapBe silica volume of bringing into by steel scrap and the silica volume that reacts generation thereof.
SVMs basicity estimation of deviation model is output as the deviation delta R of input basicity and blow end point basicity, has:
ΔR=R
input-R
aim (2)
Wherein:
R
InputBe input basicity;
R
AimBe target basicity.
Because there is different dimensions in selected inputoutput data, and different dimensional data numerically differ bigger, some in addition differ several magnitude.For avoiding variable bigger on the numerical value to fall into oblivion less variable, each input and output variable is done normalized:
In the formula (4):
I=1,2, L, 8 corresponding i variablees;
K=1,2, L, k sample in the corresponding every group of variable of n, n is a total sample number;
x
I (min)It is the minimum value in the i group variable;
x
I (max)It is the maximum value in the i group variable;
x
i(k) be k original value in the i group variable;
The computation process of SVMs is as follows:
If sampled data is { (x
1, y
1), (x
2, y
2) ..., (x
N, y
N), x wherein
i∈ R
m, y
i∈ R, i=1 ..., N.SVMs at first passes through nonlinear mapping Φ () the importation { x of sample
1, x
2..., x
NBe mapped to high-dimensional feature space F, carry out linear regression then
f(x)=w
TΦ(x)+b (4)
Wherein, w ∈ F, b represent biasing, w
TThe dot product of Φ (x) expression w and Φ (x) vector.(w, asking for b) can draw by the optimization to following formula unknown quantity
Wherein ‖ ‖ is an Euclidean distance, and C is the regularization coefficient, E
ε(g) be loss function.In order to obtain sparse solution, improve computing velocity to test sample book, the error loss function is chosen as and has the insensitive loss function of ε among the present invention:
For each sample point is introduced two relaxation factor ξ
nAnd ξ
n, optimization problem is converted into the convex quadratic programming shown in the following formula:
s.t.y
i-w
TΦ(x
i)-b≤εt
i+ξ
i (7)
For separating this optimization problem, introduce Lagrangian multiplier and this optimization problem be converted into lagrange duality problem:
0≤a
n≤C
Wherein, a
nWith
Be the Lagrange multiplier, k (x
n, x
m) be the kernel function equation, kernel function is chosen as radially basic kernel function among the present invention:
Wherein, β is the RBF width.
RBF width beta and regular terms coefficient C use the crosscheck method to determine.
Because the application of an insensitive loss function only has minority
Value is not 0, and its pairing vector is a support vector.The value of biasing b can calculate by Karush-Kuhn-Tucker (KKT) condition.
Finally set up SVMs basicity estimation of deviation model, its structure iron as shown in Figure 2:
In order to obtain actual basicity deviate, need do following anti-normalization method conversion to the output result of SVMs basicity estimation of deviation model:
ΔR=y′(k)×(y
max-y
min)+y
min (11)
Wherein:
Δ R is the estimated value of basicity deviation;
The normalized value that y ' (k) exports for model;
y
MaxBe the maximum deflection difference value in the learning sample;
y
MinBe the minimum deviation value in the learning sample.
What use when setting up SVMs basicity estimation of deviation model is the data logging that the blowing result meets processing requirement in the historical data.When current heat was predicted, some input variables can't directly obtain, and need determine these variablees.The method of determining is as follows:
For Molten Steel End Point, only after finishing blowing, could obtain, and bessemerize terminal temperature is required relatively strictness, the target temperature of terminal temperature and processing requirement is more approaching, therefore, select the target temperature of current heat processing requirement to replace the input of Molten Steel End Point as SVMs basicity estimation of deviation model.
For rhombspar add-on and the total oxygen-consumption of heat, search and the historical heat of each smelting condition of current heat success the most similar and that target call is identical (the end point carbon temperature is hit and basicity satisfies processing requirement) are with the input as current heat SVMs basicity estimation of deviation model of the rhombspar add-on of this history heat and the total oxygen-consumption of heat.
Molten iron silicon content, molten iron manganese content, molten iron phosphorus content, molten steel sulfur content, rhombspar add-on, Molten Steel End Point, lime activity and the total oxygen-consumption of heat of definite current heat are carried out normalization method with setting up the employed data of SVMs basicity estimation of deviation model.
Bring each input variable value of the current heat after the normalization method into SVMs basicity estimation of deviation model and calculate,, obtain predicted value Δ R current heat basicity deviation again to the anti-normalization method of output valve of SVMs basicity estimation of deviation model.
The experimental formula of lime adding amount is relevant with the technology basicity requirement of slag, is expressed as:
Wherein:
R
AimTarget basicity for slag;
W (CaO
Lime_eff) be the calcareous amount mark of the efficient oxidation in the lime, and w (CaO
Lime_eff)=w (CaO
Lime)-R
Aim* w (SiO
2lime);
K=2.14 is SiO
2The ratio of the molecular mass of/Si;
W[Si] be the massfraction of silicon in the molten iron;
W
HMBe molten steel quality.
Owing to also contain CaO and SiO in the rhombspar
2, can the basicity value of terminal point be exerted an influence, for eliminating the influence of rhombspar, experimental formula is revised as lime adding amount:
Wherein:
K=2.14 is SiO
2The ratio of the molecular mass of/Si;
W
DolomiteBe the rhombspar add-on;
W (CaO
Dolomite_eff) be the massfraction of the efficient oxidation calcium in the rhombspar;
And w (CaO
Dolomite_eff)=w (CaO
Dolomite)-R
Aim* w (SiO
2dolomite).
Further, use the basicity deviate of SVMs basicity estimation of deviation model prediction that the basicity parameter in the experimental formula is revised, the calculation formula that finally obtains lime adding amount is:
The lime adding amount that uses following formula to calculate has taken into full account influences the reason that basicity produces deviation, and it is compensated, and guarantees that at utmost the slag end-point alkalinity satisfies processing requirement.
After the current heat finishing blowing, if carry out slag chemical examination, and basicity of slag and end point carbon temperature satisfy processing requirement, then this heat information added standard heat storehouse, uses when later heat is set up SVMs basicity estimation of deviation model.
Example: the validity for checking institute of the present invention extracting method, adopt the actual production data of certain 150 tons of converter of steel mill to test.Select 230 groups of data according to the requirement of setting up SVMs basicity estimation of deviation model step 3, preceding 150 groups are used as modeling, and the back is used as test for 80 groups.When using definite radially sound stage width degree parameter beta of crosscheck method and regular terms coefficient C, use 150 preceding 100 groups of setting up in the modulus certificate to carry out modeling, test in back 50, definite process of parameter is used the grid search method, its Search Results as shown in Figure 3, selecting one group of optimum parameter is β=0.62, C=24000.
Use the SVMs basicity estimation of deviation model of setting up that 80 groups of test datas are carried out the basicity estimation of deviation, the relation between predicted value and the calculated value as shown in Figure 4.The accurate calculating that accurately is established as lime adding amount of basicity estimation of deviation model provides strong guarantee.With the predicted value substitution formula (14) of basicity deviation, the relation between lime adding amount that obtains and the actual add-on as shown in Figure 5.As can be seen, be that the data point of coordinate is evenly distributed on dashdotted both sides with calculated value and actual value, the absolute value major part of the error between calculated value and the actual value is all in 0.7 ton scope.
Use same data, method of the present invention and experimental formula method and increment homing method are compared, the result is as shown in table 1.
Lime adding amount Model Calculation value and the square error between the actual value based on SVMs basicity estimation of deviation are 0.2042 ton, Error Absolute Value is 96.34% less than the shared ratio of 1 ton test sample book point, and Error Absolute Value is 89.02% less than the shared ratio of 0.7 ton test sample book point.The method of this paper all is better than existing method on every index, and brings up to 0.7 ton hour when error precision by 1 ton, and the degree that accuracy rate descends is significantly less than existing method.
The comparison of table 1 and traditional method
Claims (4)
1. method for calculating converter steel melting lime adding amount, set up SVMs basicity estimation of deviation model, utilize the deviation between current heat input basicity of basicity estimation of deviation model prediction and slag end-point alkalinity, use the basicity parameter in this deviate compensation experimental formula; Described experimental formula is as follows:
Wherein, R
AimTarget basicity for slag; W (CaO
Lime_eff) be the massfraction of the efficient oxidation calcium in the lime; 2.14 be SiO
2The ratio of the molecular mass of/Si; W[Si] be the massfraction of silicon in the molten iron; W
HMBe molten steel quality;
The massfraction calculation formula of the efficient oxidation calcium is as follows in the lime:
w(CaO
lime_eff)=w(CaO
lime)-R
aim×w(SiO
2lime)
Calculate lime adding amount after the compensation basicity parameter, its calculation formula is as follows;
Wherein, K=2.14 is SiO
2The ratio of the molecular mass of/Si; W[Si] be the massfraction of silicon in the molten iron; W
HMBe molten steel quality; R
AimBe target basicity; Δ R is the deviation between input basicity and the slag end-point alkalinity; W (CaO
Lime_eff) be the calcareous amount mark of the efficient oxidation in the lime; W
DolomiteBe the rhombspar add-on; W (CaO
Dolomite_eff) be the massfraction of the efficient oxidation calcium in the rhombspar, its calculation formula is as follows:
w(CaO
dolomite_eff)=w(CaO
dolomite)-R
aim×w(SiO
2dolomite)
It is characterized in that:
The step of SVMs basicity estimation of deviation modelling is as follows:
(1) principal element of analyzing influence slag end-point alkalinity comprises slag composition, bath temperature, lime quality and blowing oxygen quantity;
(2) determine the input variable of SVMs basicity estimation of deviation model, comprise molten iron silicon content, molten iron manganese content, molten iron phosphorus content, molten steel sulfur content, rhombspar add-on, Molten Steel End Point, lime activity and the total oxygen-consumption of heat;
(3) make up the standard heat storehouse that is applicable to modeling;
(4) set up SVMs basicity estimation of deviation model;
The step that current heat lime adding amount is calculated is as follows:
(1) deviate of current heat input basicity of forecast and slag end-point alkalinity;
(2) use the basicity deviate of forecast that the basicity parameter in the experimental formula is revised;
(3) calculate current heat lime adding amount;
After the current heat finishing blowing,, then this heat information is added standard heat storehouse, use for setting up SVMs basicity estimation of deviation model if basicity of slag, endpoint carbon content and temperature satisfy processing requirement.
2. a kind of method for calculating converter steel melting lime adding amount according to claim 1 is characterized in that: described step (3) makes up and is applicable to that the standard heat storehouse of modeling comprises:
(1) data logging that selection has the terminal point slag composition to chemically examine in historical data base;
(2) data logging of selecting endpoint carbon content and temperature to hit simultaneously;
(3) select molten iron silicon content between 0.2~0.6%, and the data logging of slag end-point alkalinity between 2.8~3.5, standard heat storehouse set up.
3. a kind of method for calculating converter steel melting lime adding amount according to claim 1 is characterized in that: described step (4) is set up SVMs basicity estimation of deviation model and is comprised:
(1) the input basicity R of each heat in the base of calculation heat storehouse
Input:
Wherein, CaO
LimeThe calcium oxide content of bringing into for lime; CaO
DolomiteThe calcium oxide content of bringing into for rhombspar; SiO
2limeThe silica volume of bringing into for lime; SiO
2dolomiteBe the silica volume of bringing into by rhombspar; K=2.14 is SiO
2The ratio of the molecular mass of/Si; K * w[Si] * W
HMSilica volume for pasc reaction generation in the molten iron; SiO
2scrapBe silica volume of bringing into by steel scrap and the silica volume that reacts generation thereof;
(2) calculate deviation delta R between input basicity and the slag end-point alkalinity, as the output of SVMs basicity estimation of deviation model:
ΔR=R
input-R
aim
Wherein, R
InputBe input basicity; R
AimBe target basicity;
(3) normalization method input/output variable, and the kernel function and the loss function of selection SVMs;
(4) use the crosscheck method to determine the parameter of SVMs basicity estimation of deviation model.
4. a kind of method for calculating converter steel melting lime adding amount according to claim 1 is characterized in that: comprise in the step 1 to current heat lime adding amount calculating:
(1) search and this heat molten iron silicon content, molten iron manganese content, molten iron phosphorus content, molten steel sulfur content, rhombspar add-on, lime activity and the most close heat of terminal point target temperature in standard heat storehouse, with the reference of total oxygen-consumption of this heat as the total oxygen-consumption of current heat, use the process goal temperature of current heat of molten steel to replace Molten Steel End Point simultaneously, as the input of this heat SVMs basicity estimation of deviation model;
(2) the current heat molten iron silicon content that will determine, molten iron manganese content, molten iron phosphorus content, molten steel sulfur content, rhombspar add-on, Molten Steel End Point, lime activity and the total oxygen-consumption normalization method of heat are brought normalized each input variable value into SVMs basicity estimation of deviation Model Calculation;
(3) with the anti-normalization method of calculation result of SVMs basicity estimation of deviation model, obtain predicted value to current heat basicity deviation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102290126A CN101463407B (en) | 2008-11-22 | 2008-11-22 | Method for calculating converter steel melting lime adding amount |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102290126A CN101463407B (en) | 2008-11-22 | 2008-11-22 | Method for calculating converter steel melting lime adding amount |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101463407A CN101463407A (en) | 2009-06-24 |
CN101463407B true CN101463407B (en) | 2011-02-02 |
Family
ID=40804198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102290126A Expired - Fee Related CN101463407B (en) | 2008-11-22 | 2008-11-22 | Method for calculating converter steel melting lime adding amount |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101463407B (en) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103194563A (en) * | 2013-04-26 | 2013-07-10 | 北京科技大学 | Method for rapidly and fully deslagging converter based on physical property control of slag |
CN106503413B (en) * | 2015-08-31 | 2019-06-18 | 上海梅山钢铁股份有限公司 | A method of accurately calculating desulfurizing iron magnesium powder amount |
CN106987676B (en) * | 2017-02-13 | 2018-11-13 | 唐山不锈钢有限责任公司 | A kind of converter basicity dynamic control method |
CN107092714A (en) * | 2017-03-13 | 2017-08-25 | 昆明理工大学 | A kind of method of copper smelting by pyrometallurgy clinker slag type quantitative optimization |
CN111479935B (en) * | 2017-12-15 | 2021-11-05 | 杰富意钢铁株式会社 | Method for refining molten iron |
TWI658143B (en) * | 2018-03-30 | 2019-05-01 | 中國鋼鐵股份有限公司 | Method of reducing lime in basic oxygen furnace |
CN108982576A (en) * | 2018-06-28 | 2018-12-11 | 浙江大学 | A kind of high-precision propylene polymerization production process optimal soft survey instrument of chaos |
CN109946993A (en) * | 2019-03-18 | 2019-06-28 | 浙江大学 | A kind of thermal coupling air separation plant energy-consumption monitoring system based on hybrid modeling |
CN110904300A (en) * | 2019-12-27 | 2020-03-24 | 安徽工业大学 | Efficient dephosphorization and furnace protection method for converter slag melting point control model |
JP7469646B2 (en) * | 2020-06-05 | 2024-04-17 | 日本製鉄株式会社 | Converter blowing control device, statistical model building device, converter blowing control method, statistical model building method and program |
CN113362903B (en) * | 2021-06-02 | 2022-05-20 | 邯郸钢铁集团有限责任公司 | Method for intelligently adding lime in TSC (thyristor switched capacitor) stage of large converter |
CN115595396A (en) * | 2022-11-03 | 2023-01-13 | 山东莱钢永锋钢铁有限公司(Cn) | Method for controlling converter process and end point temperature |
CN116030900B (en) * | 2023-03-24 | 2023-06-16 | 安徽瑞邦数科科技服务有限公司 | Method, device, equipment and storage medium for controlling component content of chemical product |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1196757A (en) * | 1995-10-31 | 1998-10-21 | 伯利恒钢铁公司 | Method and apparatus to determine and control carbon content of steel in BOF vessel |
CN1757758A (en) * | 2005-11-17 | 2006-04-12 | 钢铁研究总院 | Rotary furnace steelmaking process and end point control system |
-
2008
- 2008-11-22 CN CN2008102290126A patent/CN101463407B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1196757A (en) * | 1995-10-31 | 1998-10-21 | 伯利恒钢铁公司 | Method and apparatus to determine and control carbon content of steel in BOF vessel |
CN1757758A (en) * | 2005-11-17 | 2006-04-12 | 钢铁研究总院 | Rotary furnace steelmaking process and end point control system |
Non-Patent Citations (4)
Title |
---|
郭亚芬等.转炉静态模型探索研究与应用.《工业自动化应用实践——全国(第五届)炼钢、连铸和轧钢自动化学术会议论文集》.2002,149-154. * |
陈俊东等.转炉静态机理模型与节能降耗.《河北理工学院学报》.2007,第29卷(第1期),32-35. * |
陈忠伟等.LD转炉冶炼的静态数学模型及实现.《炼钢》.2000,第16卷(第5期),31-34. * |
韩敏等.贪婪核主元模糊神经网络在转炉炼钢终点预报中的应用.《信息与控制》.2008,第37卷(第4期),494-499. * |
Also Published As
Publication number | Publication date |
---|---|
CN101463407A (en) | 2009-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101463407B (en) | Method for calculating converter steel melting lime adding amount | |
CN109447346B (en) | Converter oxygen consumption prediction method based on gray prediction and neural network combined model | |
Liu et al. | A dynamic analytics method based on multistage modeling for a BOF steelmaking process | |
CN103397140B (en) | The system and method for required refining quantity of slag during a kind of on-line prediction LF refining desulfuration | |
CN107368125B (en) | A kind of blast furnace temperature control system and method based on CBR Yu the parallel mixed inference of RBR | |
CN103645694A (en) | PS copper converter converting process intelligent decision and terminal forecast method and device | |
CN105974896B (en) | A kind of pneumatic steelmaking Optimal Control System and method based on information physical fusion | |
EP3989013A1 (en) | Method for controlling process, operation guidance method, method for operating blast furnace, method for manufacturing molten iron, and device for controlling process | |
CN113192568B (en) | Method and system for forecasting desulfurization end point of refining furnace | |
CN114091871A (en) | Blast furnace iron-making ore blending method and system | |
CN102601881A (en) | Method for monitoring on-line quality and updating prediction model of rubber hardness | |
CN104680012A (en) | Calculating model for sintering and burdening | |
CN104419799A (en) | Method for predicting carbon content of high-carbon steel on line during converter smelting | |
Ruuska et al. | Mass-balance based multivariate modelling of basic oxygen furnace used in steel industry | |
JP2012167366A (en) | Phosphorus concentration prediction apparatus, and blowing control method | |
Zhang et al. | Prediction of lime utilization ratio of dephosphorization in BOF steelmaking based on online sequential extreme learning machine with forgetting mechanism | |
CN111518980B (en) | Correction method and system for converter end point carbon content prediction model | |
CN100583136C (en) | Blast furnace iron-smelting expert system decision-making synthesis method | |
Pavlov et al. | Decision-making support in blast-furnace operation | |
CN116127345B (en) | Converter steelmaking process mode design method based on deep clustering generation countermeasure network | |
Nurkkala et al. | Blast furnace dynamics using multiple autoregressive models with exogenous inputs | |
JP2012167365A (en) | Quicklime concentration prediction apparatus, and blowing control method | |
KR101246436B1 (en) | Prediction method for product measuring of pig iron | |
KR102579633B1 (en) | Operation result prediction method, learning model learning method, operation result prediction device, and learning model learning device | |
Zhang et al. | Establishment of neural network prediction model for terminative temperature based on grey theory in hot metal pretreatment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110202 Termination date: 20131122 |