CN106801114A - A kind of blast furnace material distribution process burden distribution matrix optimization method and system - Google Patents
A kind of blast furnace material distribution process burden distribution matrix optimization method and system Download PDFInfo
- Publication number
- CN106801114A CN106801114A CN201710003901.XA CN201710003901A CN106801114A CN 106801114 A CN106801114 A CN 106801114A CN 201710003901 A CN201710003901 A CN 201710003901A CN 106801114 A CN106801114 A CN 106801114A
- Authority
- CN
- China
- Prior art keywords
- distribution matrix
- charge
- shape
- burden distribution
- furnace charge
- 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.)
- Granted
Links
Classifications
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21B—MANUFACTURE OF IRON OR STEEL
- C21B5/00—Making pig-iron in the blast furnace
- C21B5/006—Automatically controlling the process
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21B—MANUFACTURE OF IRON OR STEEL
- C21B7/00—Blast furnaces
- C21B7/18—Bell-and-hopper arrangements
- C21B7/20—Bell-and-hopper arrangements with appliances for distributing the burden
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21B—MANUFACTURE OF IRON OR STEEL
- C21B2300/00—Process aspects
- C21B2300/04—Modeling of the process, e.g. for control purposes; CII
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Manufacturing & Machinery (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Blast Furnaces (AREA)
Abstract
The present invention provides a kind of blast furnace material distribution process burden distribution matrix optimization method and system, and the method includes:The preferable shape of charge level of setting;It is target with the shape of charge level and the deviation of preferable shape of charge level that minimize next group furnace charge, optimizes burden distribution matrix;The optimal burden distribution matrix that obtains of optimization as next group furnace charge burden distribution matrix.The system includes:Setup unit, optimization unit, output unit.The present invention provides foundation for the regulation of blast furnace material distribution matrix, the method and system that operating personnel can provide according to the present invention optimize the burden distribution matrix for obtaining, the regulation of burden distribution matrix in blast furnace material distribution operation is realized, the shape of charge level under present invention optimization gained burden distribution matrix is close to preferable shape of charge level.Make gas fluid in blast furnace reasonable layout, improve raw material availability, reach the highly efficient and productive effect of blast furnace.
Description
Technical field
The invention belongs to blast furnace process automatic control technology field, more particularly to a kind of blast furnace material distribution process burden distribution matrix is excellent
Change method and system.
Background technology
Cloth operation is one of big basic operation system of blast furnace ironmaking four, decides shape and blast furnace the inside furnace charge of charge level
Layered distribution.Directly affect gas fluid distrbution, the gas utilization rate of upper blast furnace.Charge level shape is adjusted by burden distribution system
Shape, makes blast furnace inside Gas Flow reasonable layout, increases the utilization rate of coal gas, the effect of high-yield and high-efficiency is reached, to the stabilization of blast furnace
Direct motion has profound significance.However, blast furnace be one along with high temperature, high pressure, closed large-scale reactor, due to measuring instrument
The limitation of device, is difficult to directly observe blast furnace inside in normal production now.Thus, regulation burden distribution matrix is operated by cloth
The shape of charge level for obtaining cannot also be predicted.Current practices personnel are main to pass through according to stock rod height and cross temperature temperature with reference to itself
Test and estimate shape of charge level indirectly, with extremely strong subjectivity and contingency, cause that rational cloth operation cannot be carried out often.Cause
This, it is necessary to accurate Mathematical Modeling is set up, the shape of charge level under different burden distribution matrixes is obtained, instructs operating personnel to carry out reasonably
Cloth is operated.
Patent publication No. CN106011350A proposes a kind of method of estimation of the shape of charge level of blast furnace material distribution, it is already possible to
The shape of charge level under different burden distribution matrixes is obtained, the visualization of shape of charge level is realized to a certain extent.However, the method is also
Cannot realize instructing operating personnel to be adjusted burden distribution matrix, the shape of charge level under the burden distribution matrix after adjustment is reached ideal
Shape of charge level.
The content of the invention
Regarding to the issue above, the present invention proposes a kind of blast furnace material distribution process burden distribution matrix optimization method and system, makes optimization
Shape of charge level under gained burden distribution matrix finally makes blast furnace reach stabilization, direct motion, highly efficient and productive close to preferable shape of charge level
Effect.
The technical scheme is that:
A kind of blast furnace material distribution process burden distribution matrix optimization method, including:
The preferable shape of charge level of setting;
It is target with the shape of charge level and the deviation of preferable shape of charge level that minimize next group furnace charge, optimizes burden distribution matrix;
The optimal burden distribution matrix that obtains of optimization as next group furnace charge burden distribution matrix.
Burden distribution matrix is optimized using particle cluster algorithm:The burden distribution matrix for randomly generating as particle in population initial bit
Put, the deviation between shape of charge level that the burden distribution matrix corresponding to particle current location is formed and preferable shape of charge level square
With as adaptive value, iteration optimization particle position, the population optimal location the tried to achieve i.e. optimal burden distribution matrix of next group furnace charge.
The use particle cluster algorithm optimizes burden distribution matrix, including:
The burden distribution matrix for randomly generating as particle in population initial position;
Obtain the burden distribution matrix corresponding to cloth procedure parameter, furnace charge parameter and the particle current location of next group furnace charge;
Calculate the shape of charge level that burden distribution matrix of the next group furnace charge according to corresponding to particle current location is formed;
It is inclined between the shape of charge level formed with the burden distribution matrix corresponding to particle current location and preferable shape of charge level
Poor quadratic sum is used as adaptive value, iteration optimization particle position, the population optimal location the tried to achieve i.e. optimal cloth of next group furnace charge
Material matrix.
The shape of charge level that the burden distribution matrix that next group furnace charge is calculated according to corresponding to particle current location is formed, bag
Include:
Cloth procedure parameter according to next group furnace charge calculates the initial velocity that furnace charge leaves choke valve;
The initial velocity for leaving choke valve according to furnace charge calculates speed of the furnace charge in central throat tube end;
Speed according to furnace charge in central throat tube end calculates speed when furnace charge reaches chute;
Speed when chute is reached according to furnace charge calculates speed of the furnace charge in chute end;
Each chute in burden distribution matrix according to corresponding to furnace charge in the speed calculating particle current location of chute end inclines
Furnace charge forms stockpile heap cusp abscissa under the gear of angle;
Calculate the shape of charge level under last the chute inclination angle in the burden distribution matrix corresponding to particle current location, i.e., under
Shape of charge level of a collection of furnace charge under the burden distribution matrix corresponding to particle current location.
The present invention also provides a kind of blast furnace material distribution process burden distribution matrix optimization system, including:
Setup unit:The preferable shape of charge level of setting;
Optimization unit:It is target with the shape of charge level and the deviation of preferable shape of charge level that minimize next group furnace charge, optimization
Burden distribution matrix;
Output unit:The optimal burden distribution matrix that optimization is obtained is exported as the burden distribution matrix of next group furnace charge.
The optimization unit, burden distribution matrix is optimized using particle cluster algorithm:The burden distribution matrix for randomly generating is used as in population
The initial position of particle, between shape of charge level that the burden distribution matrix corresponding to particle current location is formed and preferable shape of charge level
Deviation quadratic sum as adaptive value, iteration optimization particle position, the population optimal location tried to achieve i.e. next group furnace charge is most
Excellent burden distribution matrix.
The optimization unit, including:
Initialization module:The burden distribution matrix for randomly generating as particle in population initial position;
Acquisition module:Obtain corresponding to cloth procedure parameter, furnace charge parameter and the particle current location of next group furnace charge
Burden distribution matrix;
Shape of charge level computing module:Burden distribution matrix of the next group furnace charge according to corresponding to particle current location is calculated to be formed
Shape of charge level;
Iteration optimization module:The shape of charge level and preferable charge level formed with the burden distribution matrix corresponding to particle current location
The quadratic sum of the deviation between shape is used as adaptive value, iteration optimization particle position, the population optimal location i.e. next group tried to achieve
The optimal burden distribution matrix of furnace charge.
The shape of charge level computing module, including:
First computing module:Cloth procedure parameter according to next group furnace charge calculates the initial speed that furnace charge leaves choke valve
Degree;
Second computing module:The initial velocity for leaving choke valve according to furnace charge calculates speed of the furnace charge in central throat tube end
Degree;
3rd computing module:Speed according to furnace charge in central throat tube end calculates speed when furnace charge reaches chute;
4th computing module:Speed when chute is reached according to furnace charge calculates speed of the furnace charge in chute end;
5th computing module:Burden distribution matrix according to corresponding to furnace charge in the speed calculating particle current location of chute end
In each chute inclination angle gear under furnace charge form stockpile heap cusp abscissa;
6th computing module:Calculate under last the chute inclination angle in the burden distribution matrix corresponding to particle current location
Shape of charge level, i.e. shape of charge level of the next group furnace charge under the burden distribution matrix corresponding to particle current location.
Beneficial effect:
The present invention provides foundation for the regulation of blast furnace material distribution matrix, operating personnel can according to the method for offer of the invention and
System optimizes the burden distribution matrix for obtaining, and realizes the regulation of burden distribution matrix in blast furnace material distribution operation, present invention optimization gained cloth
Shape of charge level under matrix is close to preferable shape of charge level.Make gas fluid in blast furnace reasonable layout, improve raw material availability, reach height
The highly efficient and productive effect of stove.
Brief description of the drawings
Fig. 1 is preferable shape of charge level figure in the specific embodiment of the invention;
Fig. 2 is basic shape of charge level sectional view in the specific embodiment of the invention;
Fig. 3 is charge level shaped cross figure under first chute inclination angle in the specific embodiment of the invention;
Fig. 4 is the Growth Mechanism sectional view of polycyclic distributing charge level in the specific embodiment of the invention;
Fig. 5 is for currently shape of charge level compares figure with preferable shape of charge level in the specific embodiment of the invention;
Fig. 6 is optimization process iteration diagram in the specific embodiment of the invention;
Fig. 7 is effect of optimization figure in the specific embodiment of the invention;
Fig. 8 is blast furnace material distribution process burden distribution matrix optimization method flow chart in the specific embodiment of the invention;
Fig. 9 is blast furnace material distribution process burden distribution matrix optimization method step 2 flow chart in the specific embodiment of the invention;
Figure 10 is blast furnace material distribution process burden distribution matrix optimization method step 2-3 flow charts in the specific embodiment of the invention;
Figure 11 is blast furnace material distribution process burden distribution matrix optimization system block diagram in the specific embodiment of the invention;
Figure 12 is optimization unit block diagram in the specific embodiment of the invention;
Figure 13 is shape of charge level computing module block diagram in the specific embodiment of the invention.
Specific embodiment
Specific embodiment of the invention is elaborated below in conjunction with the accompanying drawings.
Implement the method for the present invention and system without clock-type blast furnace for string tank in present embodiment, the communicated band of furnace charge is successively
Reach material feeding jar, blanking tank, central throat tube and swivel chute, final furnace charge cloth to furnace throat inside.The blast furnace material distribution process cloth
Matrix optimizing method, as shown in figure 8, including:
Step 1, with reference to feed information and the preferable shape of charge level of current working of a furnace setting;
String furnace throat position shape of charge level of the tank without clock-type blast furnace can be to ensure good coal gas flow point under preferable charge level
Cloth and raw material availability.Generally preferable shape of charge level center journey " funnel " type, both sides are " platform " type, as shown in figure 1, setting
The preferable shape of charge level γ of radial directionbX () expression formula is:
Wherein,It is the interior angle of rest (repose) of furnace charge,It is the outer angle of rest (repose) of furnace charge, x is the distance to blast furnace center, and R is furnace throat half
Footpath;
Step 2, the shape of charge level to minimize next group furnace charge and the deviation of preferable shape of charge level are target, optimize cloth
Matrix;
Sets target function is:
Constraints is:
Wherein, ψ (x, αi, mi) it is shape of charge level function, αiIt is i-th chute inclination angle, miIt is right for i-th chute inclination angle
The rotating cycle answered;I=1, wherein 2 ... ... n, n are the cloth gear number of setting during cloth, and n generally takes 4~7.
Burden distribution matrix is optimized using particle cluster algorithm in present embodiment:The burden distribution matrix for randomly generating is used as grain in population
The initial position of son, between shape of charge level that the burden distribution matrix corresponding to particle current location is formed and preferable shape of charge level
The quadratic sum of deviation as adaptive value, iteration optimization particle position, the i.e. next group furnace charge of the population optimal location tried to achieve it is optimal
Burden distribution matrix.
The use particle cluster algorithm optimizes burden distribution matrix, as shown in figure 9, including:
Step 2-1, the burden distribution matrix for randomly generating as particle in population initial position;
Parameter setting is carried out according to Liu Gang 2# blast furnaces actual conditions:Studying factors c1=c2=2, inertia weight maximum
ωmax=0.9, minimum value ωmin=0.4, maximum iteration is 100 times, population number N=100, the dimension D=10 of particle,
I-th position x of particlei=(αi1, αi2..., αi5, mi1, mi2..., mi5), wherein αi1, αi2..., αi5Chute inclination angle is represented,
[10 °, 45 °] of span, and integer is, mi1, mi2..., mi5Represent the corresponding rotating cycle in each chute inclination angle, value
Scope is [Isosorbide-5-Nitrae], and is integer, the velocity interval [v of particlemin, vmax] it is [0,2].Particle position span and
The position of particle and speed are initialized in the span of speed.
Before in the initial position of particleEach chute inclination angle is necessary for the random integers in the range of [10 °, 45 °] in dimension,
AfterwardsRotating cycle under dimension correspondence chute inclination angle, value should be the random integers in the range of [Isosorbide-5-Nitrae].
Cloth corresponding to step 2-2, the cloth procedure parameter for obtaining next group furnace charge, furnace charge parameter and particle current location
Material matrix;
The cloth procedure parameter of next group furnace charge, including:Furnace charge batch weight M, central throat tube length H, throttle valve opening κ, group
Into the radius R of the hemispherical material valve of choke valve1, furnace throat radius R, chute length L, chute rotating speed w, chute fascinate and rubbed away from b, chute
Wipe coefficient μ, stockline depth h0The time T for passing through choke valve with furnace charge.
The furnace charge parameter of next group furnace charge includes:The averag density ρ of furnace charge, the interior angle of rest (repose) of furnace chargeWith outer angle of rest (repose)
Burden distribution matrix corresponding to particle current location, including:Chute inclination alpha=[α1, α2..., αi... αn] and
Rotating cycle m=[m corresponding to each chute inclination angle1, m2..., mi... mn], 10 °≤αi≤ 45 °, 1≤mi≤4。
The charge level shape that the burden distribution matrix of step 2-3, calculating next group furnace charge according to corresponding to particle current location is formed
Shape;
Specific steps are as shown in Figure 10, including:
Step 2-3-1, the initial velocity that choke valve is left according to the cloth procedure parameter of next group furnace charge calculating furnace charge;
According to throttle valve opening κ, the radius R of the hemispherical material valve of composition choke valve1, the batch weight M of furnace charge, the density of furnace charge
ρ, furnace charge calculate the initial velocity that furnace charge leaves choke valve by the time T of choke valveWherein S0For choke valve is horizontal
Sectional area.
Step 2-3-2, the initial velocity v that choke valve is left according to furnace charge0Calculate speed of the furnace charge in central throat tube end
v1;
Furnace charge can approximately regard furnace charge as and do the movement of falling object by central throat tube, only receive Action of Gravity Field, height of drop
Centered on length of throat H, meet following relation:G=9.8m/s2It is acceleration of gravity.
Step 2-3-3, the speed according to furnace charge in central throat tube end calculate speed v when furnace charge reaches chute2;
Furnace charge speed during being clashed into chute would detract from losing, and direction will change, and meet following relation:v2
=kcos α v1, k is furnace charge impact velocity attenuation coefficient, k=0.8.
Step 2-3-4, according to furnace charge reach chute when speed calculate furnace charge chute end speed v3;
When furnace charge is moved on chute, respectively by self gravitation, the holding power of chute, frictional force, coriolis force and by
The centrifugal force that chute rotation causes, these power are decomposed, and furnace charge will be accelerated along chute direction, meet such as ShiShimonoseki
System:
Step 2-3-5, according to furnace charge chute end speed v3Calculate the burden distribution matrix corresponding to particle current location
In each chute inclination angle gear under furnace charge form stockpile heap cusp abscissa;
Furnace charge is acted on after leaving chute by the drag of gravity and rising coal gas, inclined throw movements is done in dead zone, due to coal gas
Drag smaller negligible do not remember.Therefore furnace charge can be approximately considered and be only subject to gravity.Obtained under n chute inclination angle by kinematical equation
Furnace charge forms stockpile heap cusp abscissa xnFor:
The material under last the chute inclination angle in burden distribution matrix corresponding to step 2-3-6, calculating particle current location
Face shape, i.e. shape of charge level of the next group furnace charge under the burden distribution matrix corresponding to particle current location.
Step 2-3-6-1, the basic shape of charge level of structure;
As shown in Fig. 2 the cross section profile of shape of charge level of the furnace charge in furnace throat can approximately be regarded as be made up of straight line CB, CA,
The inside and outside angle of rest (repose) of furnace charge is the angle of this two straight lines and horizontal plane.
The function for building basic shape of charge level is:
In formula, m=1 is the rotating cycle of chute, and line segment BCA is basic shape of charge level, (xn, yh) material is formed by furnace charge
Heap heap cusp C coordinates, rLIt is straight line CB and the abscissa of horizontal line OL intersection points, rRIt is straight line CA and the abscissa of furnace wall intersection point.
Step 2-3-6-2, the relation built between the shape of charge level and basic shape of charge level of output are shape of charge level function,
Shape of charge level described by the shape of charge level function of output is the superposition of basic shape of charge level;
Shape of charge level is mainly determined that final shape of charge level is the superposition of basic shape of charge level by burden distribution matrix, builds material
Face shape function is as follows:
In formula, γ (x) is shape of charge level function, α=[α1, α2..., αi... αn] it is chute inclination angle, m=[m1,
m2..., mi... mn] it is the corresponding rotating cycle in chute inclination angle.
Cloth body in burden distribution matrix corresponding to step 2-3-6-3, structure particle current location under each chute inclination angle
Cloth volume under product constraints, i.e. each chute inclination angle is equal under the chute inclination angle under shape of charge level and upper chute inclination angle
The volume integral at shape of charge level interval;
Due to furnace charge, volume does not change in motion process, therefore the shape of charge level and initial shape of charge level (water of output
Plane) between furnace charge volume and furnace charge enter material feeding jar volume it is equal, furnace charge volume meets following constraint:
V=M/ ρ (8)
In formula, V is furnace charge volume.
Built according to formula (8), (9) under the cloth volume constraint condition under each chute inclination angle, i.e. each chute inclination angle
Cloth volume be equal to the chute inclination angle under shape of charge level with a upper chute inclination angle shape of charge level be spaced volume integral:
Cloth volume constraint condition is as follows:
In formula, miIt is the rotating cycle under i-th chute inclination angle, γiX () is the charge level shape of output under i-th chute inclination angle
Shape.
Step 2-3-6-4, according to the cloth volume constraint condition under first chute inclination angle, calculate particle current location institute
Furnace charge forms stockpile heap cusp ordinate under first chute inclination angle in corresponding burden distribution matrix;
The volume constraint that furnace charge meets under first chute inclination angle is as follows:
With the carrying out of cloth, material stack height yh1Gradually increase.When cloth is completed, by iteration furnace charge heap cusp height yh1,
First lower height y for forming stockpile heap point in chute inclination angle can be tried to achieveh1Furnace charge forms stockpile under i.e. first chute inclination angle
Heap cusp ordinate, then heap cusp C1 C1(xn1, yh1) can uniquely determine.And the interior angle of rest (repose), outer angle of rest (repose) due to furnace charge are constant, therefore
Straight line C1B1With straight line C1A1Slope it is constant, it is known that on straight line a little and its slope can uniquely determine straight line.Then can be unique
Try to achieve the intersection points B of left and right charge level and initial charge level1、A1.Can obtain its abscissa rL1、rR1。
Furnace charge institute under first chute inclination angle in step 2-3-6-5, the burden distribution matrix according to corresponding to particle current location
Stockpile heap cusp coordinate and shape of charge level function are formed, first chute in the burden distribution matrix corresponding to particle current location is solved
Shape of charge level under inclination angle.
As shown in figure 3, it is the origin of coordinates to set blast furnace center, the horizontal plane where straight line OL is initial charge level, then basis
Shape of charge level in the burden distribution matrix that formula (6), (7) can obtain corresponding to particle current location under first chute inclination angle is:
In formula, (xn1, yh1) by forming stockpile heap cusp C under first chute inclination angle1Coordinate, rL1It is straight line C1B1With
The abscissa of initial charge level intersection point, rR1It is straight line C1A1With the abscissa of initial charge level intersection point.
The body that furnace charge meets under each chute inclination angle in step 2-3-7, the burden distribution matrix according to corresponding to particle current location
The shape of charge level of the shape of charge level and basic charge level under constraint, previous chute inclination angle is accumulated, iteration obtains second chute successively
Shape of charge level under inclination angle to complete n-th each chute inclination angle of cloth, the shape of charge level under last the chute inclination angle tried to achieve, i.e. grain
The shape of charge level that burden distribution matrix corresponding to sub- current location is formed.
As shown in figure 4, with the carrying out of cloth under second chute inclination angle, material stack height yh2Gradually increase.By continuous
Change furnace charge heap cusp height yh2, it is met formula (11), the height y of furnace charge heap point under second chute inclination angle can be tried to achieveh2I.e.
Furnace charge forms stockpile heap cusp ordinate under second chute inclination angle, then heap cusp C2(xn2, yh2) can uniquely determine, and due to
The inside and outside angle of rest (repose) of furnace charge is constant, therefore straight line C2B2With straight line C2A2Slope it is constant, then can uniquely try to achieve under second chute inclination angle
The intersection points B of charge level under the left and right charge level of furnace charge formation and first chute inclination angle2、A2.Can obtain rL2, rR2.Therefore second can be obtained
Shape of charge level γ under individual chute inclination angle2(x).Similarly, the shape of charge level under the three to n-th chute inclination angle can successively be tried to achieve.
Polycyclic distributing is the superposition of institute's cloth furnace charge shape of charge level under each chute inclination angle, and current shape of charge level will be under
The initial shape of charge level of one circle cloth, i.e., the shape of charge level under first chute inclination angle will enclose the initial charge level of cloth as second
Shape, is below the shape of charge level behind the complete second chute inclination angle of cloth to the complete n-th chute inclination angle of cloth.
Shape of charge level under second and n-th chute inclination angle is respectively:
In formula, (xni, yni) by forming stockpile heap cusp coordinate, γ under i-th chute inclination angleiX () is inclined for i-th chute
Shape of charge level under angle.rLi、rRiIt is institute's cloth furnace charge under i-th chute inclination angle and upper charge level γi-1X () is in heap cusp (xni,
yni) left and right sides intersection points Bi、AiAbscissa.
Shape of charge level γ under shape of charge level γ (x) of final output and last chute inclination anglenX () is identical.
Step 2-4, the shape of charge level and preferable shape of charge level that are formed with the burden distribution matrix corresponding to particle current location
Between deviation quadratic sum as adaptive value, iteration optimization particle position, the population optimal location tried to achieve i.e. next group furnace charge
Optimal burden distribution matrix.
Step 2-4-1, adaptive value evaluation:
Adaptation value function in present embodiment be preferable shape of charge level function with particle group optimizing in optimize the cloth for obtaining
Material matrix forms the quadratic sum of the deviation of shape of charge level, i.e.,Wherein, initial charge level is current cloth
Shape of charge level under material matrix, Fig. 5 is the contrast of current shape of charge level and preferable shape of charge level.For each particle, worked as
The adaptive value of front position is compared with the adaptive value of the excellent position of the history that it is experienced, if the adaptive value of its current location is excellent
In the adaptive value of its history optimal location, then its history optimal value is replaced, with current position come optimal instead of history
Place value.The adaptive value of the individual history optimal location of each particle and the adaptive value of the history optimal location of colony are compared
Compared with, when the adaptive value of personal best particle is better than the adaptive value of history optimal location, then colony's optimal location is replaced,
Colony's optimal location is replaced with individual history optimal location.
Step 2-4-2, particle state update:
The position of the state of particle including particle and the speed of particle, using formula (15) to the speed of particle in present embodiment
Degree is updated, and the position of particle is updated using formula (16).
In formula,It is the position of k+1 moment particles,It is the position of k moment particles,It is the speed of k+1 moment particles
Degree,It is the speed of k moment particles, ξ1With ξ2The random number in the range of [0,1] is,For particle is optimal in k moment individualities
Position,It is particle in k moment population optimal locations;
Step 2-4-3, stopping iteration judging:Iterative process as shown in fig. 6, as can be seen from the figure when 20 step of iteration,
Preferable effect is reached, when iterations reaches maximum iteration 100 times, has then stopped the optimal position of iteration output current group
Put, i.e., optimal burden distribution matrix is precedingIndividual is chute inclination angle, afterBefore number is followed successively byIt is corresponding under individual chute inclination angle to slip
Groove rotating cycle, because chute is to turn to less inclination angle by larger inclination angle during cloth, it is therefore desirable to chute inclination angle
Angle it is descending sort successively, while following the chute inclination angle to carry out order change corresponding rotating cycle under the chute inclination angle
Change, obtain the burden distribution matrix of final optimization pass.Otherwise, step 2-4-1 is jumped to, proceeds iteration.Final optimization pass is obtained
Burden distribution matrix calls shape of charge level computational methods, obtains shape of charge level and preferable charge level shape under the burden distribution matrix as shown in Figure 7
The comparison diagram of shape.
Step 3, the optimal burden distribution matrix that obtains of optimization as next group furnace charge burden distribution matrix.
The present invention also provides a kind of blast furnace material distribution process burden distribution matrix optimization system, as shown in figure 11, including:
Setup unit:The preferable shape of charge level of setting;
Optimization unit:It is target with the shape of charge level and the deviation of preferable shape of charge level that minimize next group furnace charge, optimization
Burden distribution matrix;
Output unit:The optimal burden distribution matrix that optimization is obtained is exported as the burden distribution matrix of next group furnace charge.
The optimization unit, burden distribution matrix is optimized using particle cluster algorithm:The burden distribution matrix for randomly generating is used as in population
The initial position of particle, between shape of charge level that the burden distribution matrix corresponding to particle current location is formed and preferable shape of charge level
Deviation quadratic sum as adaptive value, iteration optimization particle position, the population optimal location tried to achieve i.e. next group furnace charge is most
Excellent burden distribution matrix.
The optimization unit, as shown in figure 12, including:
Initialization module:The burden distribution matrix for randomly generating as particle in population initial position;
Acquisition module:Obtain corresponding to cloth procedure parameter, furnace charge parameter and the particle current location of next group furnace charge
Burden distribution matrix;
Shape of charge level computing module:Burden distribution matrix of the next group furnace charge according to corresponding to particle current location is calculated to be formed
Shape of charge level;
Iteration optimization module:The shape of charge level and preferable charge level formed with the burden distribution matrix corresponding to particle current location
The quadratic sum of the deviation between shape is used as adaptive value, iteration optimization particle position, the population optimal location i.e. next group tried to achieve
The optimal burden distribution matrix of furnace charge.
The shape of charge level computing module, as shown in figure 13, including:
First computing module:Cloth procedure parameter according to next group furnace charge calculates the initial speed that furnace charge leaves choke valve
Degree;
Second computing module:The initial velocity for leaving choke valve according to furnace charge calculates speed of the furnace charge in central throat tube end
Degree;
3rd computing module:Speed according to furnace charge in central throat tube end calculates speed when furnace charge reaches chute;
4th computing module:Speed when chute is reached according to furnace charge calculates speed of the furnace charge in chute end;
5th computing module:Burden distribution matrix according to corresponding to furnace charge in the speed calculating particle current location of chute end
In each chute inclination angle gear under furnace charge form stockpile heap cusp abscissa;
6th computing module:Calculate under last the chute inclination angle in the burden distribution matrix corresponding to particle current location
Shape of charge level, i.e. shape of charge level of the next group furnace charge under the burden distribution matrix corresponding to particle current location.
Blast furnace material distribution process burden distribution matrix optimization method and system that the present invention is provided, can optimize to burden distribution matrix,
Make the shape of charge level under optimization gained burden distribution matrix close to preferable shape of charge level, the tune of burden distribution matrix in being operated to blast furnace material distribution
It is whole with important directive significance.
Claims (8)
1. a kind of blast furnace material distribution process burden distribution matrix optimization method, it is characterised in that including:
The preferable shape of charge level of setting;
It is target with the shape of charge level and the deviation of preferable shape of charge level that minimize next group furnace charge, optimizes burden distribution matrix;
The optimal burden distribution matrix that obtains of optimization as next group furnace charge burden distribution matrix.
2. method according to claim 1, it is characterised in that burden distribution matrix is optimized using particle cluster algorithm:Randomly generate
Burden distribution matrix as particle in population initial position, the charge level shape that the burden distribution matrix corresponding to particle current location is formed
The quadratic sum of the deviation between shape and preferable shape of charge level is used as adaptive value, and iteration optimization particle position, the population tried to achieve is optimal
Position is the optimal burden distribution matrix of next group furnace charge.
3. method according to claim 2, it is characterised in that the use particle cluster algorithm optimizes burden distribution matrix, including:
The burden distribution matrix for randomly generating as particle in population initial position;
Obtain the burden distribution matrix corresponding to cloth procedure parameter, furnace charge parameter and the particle current location of next group furnace charge;
Calculate the shape of charge level that burden distribution matrix of the next group furnace charge according to corresponding to particle current location is formed;
Deviation between the shape of charge level that is formed with the burden distribution matrix corresponding to particle current location and preferable shape of charge level
Quadratic sum is used as adaptive value, iteration optimization particle position, the optimal cloth square of the i.e. next group furnace charge of the population optimal location tried to achieve
Battle array.
4. method according to claim 3, it is characterised in that the calculating next group furnace charge is according to particle current location institute
The shape of charge level that corresponding burden distribution matrix is formed, including:
Cloth procedure parameter according to next group furnace charge calculates the initial velocity that furnace charge leaves choke valve;
The initial velocity for leaving choke valve according to furnace charge calculates speed of the furnace charge in central throat tube end;
Speed according to furnace charge in central throat tube end calculates speed when furnace charge reaches chute;
Speed when chute is reached according to furnace charge calculates speed of the furnace charge in chute end;
Each chute inclination angle shelves in burden distribution matrix according to corresponding to furnace charge in the speed calculating particle current location of chute end
The lower furnace charge in position forms stockpile heap cusp abscissa;
Calculate the shape of charge level under last the chute inclination angle in the burden distribution matrix corresponding to particle current location, i.e. next group
Shape of charge level of the furnace charge under the burden distribution matrix corresponding to particle current location.
5. a kind of blast furnace material distribution process burden distribution matrix optimizes system, it is characterised in that including:
Setup unit:The preferable shape of charge level of setting;
Optimization unit:It is target with the shape of charge level and the deviation of preferable shape of charge level that minimize next group furnace charge, optimizes cloth
Matrix;
Output unit:The optimal burden distribution matrix that optimization is obtained is exported as the burden distribution matrix of next group furnace charge.
6. system according to claim 5, it is characterised in that the optimization unit, cloth is optimized using particle cluster algorithm
Matrix:The burden distribution matrix for randomly generating as particle in population initial position, the burden distribution matrix corresponding to particle current location
The quadratic sum of the deviation between the shape of charge level that is formed and preferable shape of charge level as adaptive value, iteration optimization particle position,
The population optimal location the tried to achieve i.e. optimal burden distribution matrix of next group furnace charge.
7. system according to claim 6, it is characterised in that the optimization unit, including:
Initialization module:The burden distribution matrix for randomly generating as particle in population initial position;
Acquisition module:Obtain the cloth corresponding to cloth procedure parameter, furnace charge parameter and the particle current location of next group furnace charge
Matrix;
Shape of charge level computing module:Calculate the material that burden distribution matrix of the next group furnace charge according to corresponding to particle current location is formed
Face shape;
Iteration optimization module:The shape of charge level and preferable shape of charge level formed with the burden distribution matrix corresponding to particle current location
Between deviation quadratic sum as adaptive value, iteration optimization particle position, the population optimal location tried to achieve i.e. next group furnace charge
Optimal burden distribution matrix.
8. system according to claim 7, it is characterised in that the shape of charge level computing module, including:
First computing module:Cloth procedure parameter according to next group furnace charge calculates the initial velocity that furnace charge leaves choke valve;
Second computing module:The initial velocity for leaving choke valve according to furnace charge calculates speed of the furnace charge in central throat tube end;
3rd computing module:Speed according to furnace charge in central throat tube end calculates speed when furnace charge reaches chute;
4th computing module:Speed when chute is reached according to furnace charge calculates speed of the furnace charge in chute end;
5th computing module:In burden distribution matrix according to corresponding to furnace charge in the speed calculating particle current location of chute end
Furnace charge forms stockpile heap cusp abscissa under the gear of each chute inclination angle;
6th computing module:Calculate the charge level under last the chute inclination angle in the burden distribution matrix corresponding to particle current location
Shape, i.e. shape of charge level of the next group furnace charge under the burden distribution matrix corresponding to particle current location.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710003901.XA CN106801114B (en) | 2017-01-04 | 2017-01-04 | A kind of blast furnace material distribution process burden distribution matrix optimization method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710003901.XA CN106801114B (en) | 2017-01-04 | 2017-01-04 | A kind of blast furnace material distribution process burden distribution matrix optimization method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106801114A true CN106801114A (en) | 2017-06-06 |
CN106801114B CN106801114B (en) | 2019-02-26 |
Family
ID=58984287
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710003901.XA Active CN106801114B (en) | 2017-01-04 | 2017-01-04 | A kind of blast furnace material distribution process burden distribution matrix optimization method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106801114B (en) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107119159A (en) * | 2017-06-14 | 2017-09-01 | 内蒙古科技大学 | A kind of optimization method for the blast furnace material distribution process burden distribution matrix that there is integer programming problem |
CN107245540A (en) * | 2017-06-14 | 2017-10-13 | 内蒙古科技大学 | A kind of control strategy of blast furnace material distribution process radial direction thickness of feed layer distribution |
CN107992974A (en) * | 2017-12-13 | 2018-05-04 | 沈阳东博热工科技有限公司 | A kind of optimization method of vacuum heat burden distribution matrix |
CN108277314A (en) * | 2018-02-27 | 2018-07-13 | 首钢京唐钢铁联合有限责任公司 | A kind of distributing method and device |
CN111304388A (en) * | 2020-03-27 | 2020-06-19 | 武汉钢铁有限公司 | Method for optimizing blast furnace gas flow distribution by using scanning radar |
CN112226557A (en) * | 2020-09-27 | 2021-01-15 | 北京科技大学 | Blast furnace burden distribution process evaluation method and system, readable storage medium and application thereof |
CN114807468A (en) * | 2022-04-02 | 2022-07-29 | 武汉钢铁有限公司 | Method for improving blast furnace gas utilization rate based on blast furnace charge level monitoring |
CN114800828A (en) * | 2022-04-09 | 2022-07-29 | 沈阳建筑大学 | Intelligent optimization method for load distribution of prefabricated part concrete distribution robot |
CN113637819B (en) * | 2021-08-16 | 2022-10-18 | 中冶东方工程技术有限公司 | Blast furnace material distribution method and system based on deep reinforcement learning |
CN115584370A (en) * | 2022-09-01 | 2023-01-10 | 武安市裕华钢铁有限公司 | Method for improving air permeability of soft melting zone in blast furnace large slag quantity smelting |
CN116694837A (en) * | 2023-08-09 | 2023-09-05 | 山西建龙实业有限公司 | Distribution method for high-proportion bituminous coal injection |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102732659A (en) * | 2011-04-11 | 2012-10-17 | 宝山钢铁股份有限公司 | Burden surface profile control method of blast furnace burden distribution and control system |
CN106011350A (en) * | 2016-07-05 | 2016-10-12 | 东北大学 | Method and system for estimating charge surface shapes in blast furnace charge distribution process |
-
2017
- 2017-01-04 CN CN201710003901.XA patent/CN106801114B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102732659A (en) * | 2011-04-11 | 2012-10-17 | 宝山钢铁股份有限公司 | Burden surface profile control method of blast furnace burden distribution and control system |
CN106011350A (en) * | 2016-07-05 | 2016-10-12 | 东北大学 | Method and system for estimating charge surface shapes in blast furnace charge distribution process |
Non-Patent Citations (1)
Title |
---|
马财生等: "无钟高炉环形布料优化及布料精度", 《钢铁研究学报》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107245540A (en) * | 2017-06-14 | 2017-10-13 | 内蒙古科技大学 | A kind of control strategy of blast furnace material distribution process radial direction thickness of feed layer distribution |
CN107119159B (en) * | 2017-06-14 | 2019-02-22 | 内蒙古科技大学 | It is a kind of that there are the optimization methods of the blast furnace material distribution process burden distribution matrix of integer programming problem |
CN107119159A (en) * | 2017-06-14 | 2017-09-01 | 内蒙古科技大学 | A kind of optimization method for the blast furnace material distribution process burden distribution matrix that there is integer programming problem |
CN107992974B (en) * | 2017-12-13 | 2021-03-05 | 沈阳东博热工科技有限公司 | Optimization method of vacuum heat treatment distribution matrix |
CN107992974A (en) * | 2017-12-13 | 2018-05-04 | 沈阳东博热工科技有限公司 | A kind of optimization method of vacuum heat burden distribution matrix |
CN108277314A (en) * | 2018-02-27 | 2018-07-13 | 首钢京唐钢铁联合有限责任公司 | A kind of distributing method and device |
CN111304388B (en) * | 2020-03-27 | 2021-12-17 | 武汉钢铁有限公司 | Method for optimizing blast furnace gas flow distribution by using scanning radar |
CN111304388A (en) * | 2020-03-27 | 2020-06-19 | 武汉钢铁有限公司 | Method for optimizing blast furnace gas flow distribution by using scanning radar |
CN112226557A (en) * | 2020-09-27 | 2021-01-15 | 北京科技大学 | Blast furnace burden distribution process evaluation method and system, readable storage medium and application thereof |
CN113637819B (en) * | 2021-08-16 | 2022-10-18 | 中冶东方工程技术有限公司 | Blast furnace material distribution method and system based on deep reinforcement learning |
CN114807468A (en) * | 2022-04-02 | 2022-07-29 | 武汉钢铁有限公司 | Method for improving blast furnace gas utilization rate based on blast furnace charge level monitoring |
CN114807468B (en) * | 2022-04-02 | 2023-10-24 | 武汉钢铁有限公司 | Blast furnace charge level monitoring-based method for improving blast furnace gas utilization rate |
CN114800828A (en) * | 2022-04-09 | 2022-07-29 | 沈阳建筑大学 | Intelligent optimization method for load distribution of prefabricated part concrete distribution robot |
CN114800828B (en) * | 2022-04-09 | 2024-04-05 | 沈阳建筑大学 | Intelligent optimization method for load distribution of prefabricated part concrete distribution robot |
CN115584370A (en) * | 2022-09-01 | 2023-01-10 | 武安市裕华钢铁有限公司 | Method for improving air permeability of soft melting zone in blast furnace large slag quantity smelting |
CN115584370B (en) * | 2022-09-01 | 2023-10-13 | 武安市裕华钢铁有限公司 | Method for improving air permeability of soft melting belt in large-slag-quantity smelting of blast furnace |
CN116694837A (en) * | 2023-08-09 | 2023-09-05 | 山西建龙实业有限公司 | Distribution method for high-proportion bituminous coal injection |
CN116694837B (en) * | 2023-08-09 | 2023-10-13 | 山西建龙实业有限公司 | Distribution method for high-proportion bituminous coal injection |
Also Published As
Publication number | Publication date |
---|---|
CN106801114B (en) | 2019-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106801114B (en) | A kind of blast furnace material distribution process burden distribution matrix optimization method and system | |
CN106011350B (en) | A kind of blast furnace material distribution process shape of charge level method of estimation and system | |
CN104531924B (en) | A kind of blast furnace burden distribution Real-time Forecasting System and method | |
CN104133945B (en) | A kind of control method of blast furnace material distribution process radial direction ore coke ratio | |
CN107723399B (en) | A kind of blast furnace material distribution intelligent monitor system and method for adjustment | |
CN107169558A (en) | A kind of Modified particle swarm optimization method for realizing engineering rock mass mechanics parameter inverting | |
CN104694680B (en) | A kind of radially-arranged control method of blast furnace bed structure | |
CN110400009A (en) | Blast furnace ironmaking Multiobjective Intelligent optimization method based on self-adapted genetic algorithm | |
CN102653801B (en) | Bell-less blast furnace top burden distribution control method established on basis of genetic algorithm | |
CN103173584A (en) | Blast furnace burden-distribution control system with self-learning control function | |
CN110427715B (en) | Method for predicting furnace hearth thermal state trend based on time sequence and multiple dimensions of blast furnace | |
CN103966373A (en) | Bell-less material distributing method capable of ensuring stable running of blast furnace | |
CN106582967A (en) | Steel ball supplementing method for double-inlet and double-outlet coal mill | |
CN107245540B (en) | A kind of control strategy of blast furnace material distribution process radial direction thickness of feed layer distribution | |
Zhou et al. | Grey correlation analysis of factors influencing maldistribution in feeding device of copper flash smelting | |
CN206593486U (en) | A kind of plate belt feeder unloading type vertical cooling machine | |
CN107119159B (en) | It is a kind of that there are the optimization methods of the blast furnace material distribution process burden distribution matrix of integer programming problem | |
CN103216825A (en) | Circular-arc-shaped U valve material returning device for high-density circulating fluidized bed utilizing B-class particle | |
CN116694837A (en) | Distribution method for high-proportion bituminous coal injection | |
CN113656966A (en) | Blast furnace bell-less top on-line material distribution model simulation method | |
CN206891200U (en) | There is the vertical sinter cooler of material position adjusting means | |
JP2017160283A (en) | Temperature distribution control method and control device of vertical pyrolysis furnace | |
CN206593485U (en) | A kind of vertical sinter cooler | |
JP7077842B2 (en) | Blast furnace charge distribution prediction method, program and computer storage medium | |
CN110083936A (en) | Consider form parameter is crushed gravel particle generation method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |