CN103439926A - Gas optimization scheduling device of iron and steel enterprise - Google Patents

Abstract

The invention relates to a gas optimization scheduling device of an iron and steel enterprise. The scheduling device is connected with an energy management system and comprises a gas managing device, an optimizing device and a display device, wherein the gas managing device, the optimizing device and the display device are sequentially connected. The gas managing device acquires gas consumption amount historical data, gas generation amount historical data and storage amount data of a gas cabinet in each procedure through the energy management system of the enterprise to predict the production amount and the consumption amount of gas at present or in future; on the basis of the predicted data, the optimizing device adopts improved distribution estimation algorithms to optimize distribution of the gas, and the optimization principle aims to make the scattering amount smallest and the purchasing amount largest; the display device displays the gas distribution scheme obtained by the optimizing device for reference of users, so that the gas can be reasonably distributed. Compared with the prior art, the gas optimization scheduling device has the advantages of reducing the scattering amount and the purchasing amount of the gas and the like.

Description

A kind of iron and steel enterprise coal gas Optimized Operation device

Technical field

The present invention relates to a kind of gas dispatching device, especially relate to a kind of iron and steel enterprise coal gas Optimized Operation device.

Background technology

Iron and steel is the industry of high energy consumption, high pollution, maximum discharge.Energy cost accounts for 30% of its manufacturing cost, takes full advantage of coal gas, and reducing bleeding rate is the important channel of energy-saving and emission-reduction.Iron and steel enterprise relates to a lot of operations, and coal gas demand type and the demand of each operation are different, operation in coal gas demand in the same time not also in dynamic change.Due to the complex relationship between coal gas and operation, need a kind of work of coordinating the two of installing badly.This matching requirements can be planned as a whole whole gas pipe, possesses complete gas allocation and controls ability.

In Chinese patent " iron and steel enterprise's coke-oven gas cabinet position prediction balance method " (numbering: CN102109837A), utilize factor association analysis theory with data-driven modeling Regression Forecasting Technology, iron and steel enterprise's coke-oven coal gas holder position to be predicted.And wish that the result dope can make the dispatcher understand in advance the variation situation of gas chamber, regulate in time user's gas using quantity, reduce the amount of diffusing of coal gas, improve the utilization factor of coal gas.But this patent does not relate to concrete dispatching method, there is no the prediction of truly realizing and the combination of dispatching.

At Chinese patent " a kind of iron and steel enterprise energy simulating scenes preparation method ", (numbering: CN102073951A), be by Construction of A Model, apparatus for establishing model and product model, carry out the scene establishment.This emulation tries hard to reproduce the production run of iron and steel enterprise.Wish can be used in method general, that can reflect steel product thing mass flow change, Product Process, unusual service condition, overhaul of the equipments.But this patent just rests on visual degree, and be to pay close attention to technological parameter more, the situation of production status, do not reach the purpose to the scheduling of coal gas intelligent optimization.

In Chinese patent " a kind of smelter by-product gas off-line analysis on-line scheduling method " (numbering: CN102915006A).This dispatching method has a data collection server to be connected with coal gas production and marketing management system, and this data server can carry out on-line optimizing scheduling according to the mathematical forecasting model established and scheduling model.This is a kind of off-line analysis on-line scheduling method.But this dispatching method just considers how to make coal gas diffusion amount minimum, and does not consider the coal gas purchase cost, and simplex side's computing method of adopting of this patent, easily be absorbed in without the quagmire of separating.The most important thing is, this patent just predicted value based on current time is dispatched coal gas, and does not consider next coal gas surplus situation constantly.If gas chamber corresponding to the coal gas of the surplus of current time deposit in like this, and the also surplus and other coal gas Shortcomings phenomenons will be had to again this coal gas diffusion in atmosphere of next this coal gas constantly, other coal gas are again in not enough state simultaneously.

Summary of the invention

Purpose of the present invention is exactly to provide a kind of iron and steel enterprise coal gas Optimized Operation device in order to overcome the defect that above-mentioned prior art exists.

Purpose of the present invention can be achieved through the following technical solutions:

A kind of iron and steel enterprise coal gas Optimized Operation device, this dispatching device is connected with energy management system, it is characterized in that, and described dispatching device comprises gas pipe device, optimizer and the display device connected successively; The gas pipe device utilizes the energy management system of enterprise to obtain each operation gas consumption amount and growing amount historical data and gas chamber data on stock, and coal gas current time and turnout and the consumption in the following moment are predicted; Optimizer, based on above-mentioned predicted data, adopts improved Estimation of Distribution Algorithm to be optimized gas allocation, and the principle of optimization is to make the amount of diffusing of coal gas minimum minimum with purchase volume; The gas allocation scheme that display device display optimization device obtains, for reference, carry out reasonably gas allocation.

Described optimizer, based on above-mentioned predicted data, adopts improved Estimation of Distribution Algorithm to be optimized and to be specially gas allocation:

Take coal gas diffusion amount and purchase volume minimum be the final optimization pass target, meeting under the prerequisite of process constraint, operation is carried out to priority definition, then first coal gas is carried out to original allocation, the principle of original allocation is the priority according to operation, take and consume coke-oven gas and coal gas of converter is principle as far as possible, the minimum heat that meets operation of take requires as benchmark, finally coal gas is carried out to secondary distribution, this distribution belongs to distribution according to need, iteration is distributed, and utilizes improved Estimation of Distribution Algorithm to regulate the gas allocation amount of each operation.

Described operation is carried out priority definition and is specially:

At first will be by the nonadjustable operation priority ordering of gas consumption structure during Operation Sequencing, it is less important preferential with the operation that consumes coke-oven gas, again preferential with the operation that consumes coal gas of converter.

Describedly first coal gas is carried out to original allocation and is specially:

The minimum heat demand that meets operation of take is principle, at first distribute coke-oven gas, travel through the operation relevant to coke-oven gas, until coke-oven gas distributes, once coke-oven gas distributes the distribution that has just entered coal gas of converter, the same traversal operation relevant to coal gas of converter, until, after fully coal gas of converter being distributed, finally do same operation to blast furnace gas.

Describedly coal gas carried out to secondary distribution be specially:

1) determine coal gas side's formula case of current time and the relation of next gas chamber storage constantly;

2) one group of variable data[to be solved is set];

3) determine optimization aim, comprise that the amount of diffusing is minimum and the coal gas purchase volume is minimum;

4) adopt improved Estimation of Distribution Algorithm to be distributed, be about to searching process and be divided into two levels, evenly divide whole solution interval, with scatter searching, guide the Evolution of Population direction, the interval produced according to scatter search produces filial generation, with histogram model statistic mass's probability again.

Coal gas side's formula case of described definite current time is specially with the relation of next gas chamber storage constantly:

1), if be carved with the coal gas surplus when next, go up a moment and will gas chamber be put to sky as far as possible and go to meet the next storage demand of coal gas constantly;

2) if the next situation of coal gas Shortcomings constantly, current time will reserve corresponding coal gas in gas chamber as far as possible, to meet the next coal gas deficiency situation that constantly may exist.

Described variable to be solved specifically comprises the amount of blast furnace gas and coke-oven gas equivalent exchange, the amount of blast furnace gas and coal gas of converter equivalent exchange, the amount of coke-oven gas and coal gas of converter equivalent exchange, three kinds of coal gas offer respectively the consumption again of the operation with corresponding mode of communicating, three kinds of coal gas offer respectively the disappearance amount of the operation with corresponding mode of communicating and the amount of next blast furnace gas constantly and coke-oven gas equivalent exchange, the amount of blast furnace gas and coal gas of converter equivalent exchange, the amount of coke-oven gas and coal gas of converter equivalent exchange.

The improved Estimation of Distribution Algorithm of described employing distributes to be specially sets up population space, divide equally spaced interval in space, there is own corresponding probability in each interval, this probability is characterized by the probability that excellent individual occurs, by each, interval probability is evaluated the member in reference set, this algorithm is provided with two reference sets, excellent solution reference set and the good reference set of diversity, pass through the calculating of genetic operator by the interval of selecting at random respectively respectively from two reference sets, generate an interval number, from this interval, the random number that generates forms the next generation again, of future generation through calculating, upgrade again probability, upgrade the reference set member, so constantly repeat until meet termination condition, termination condition is for arriving the number of iterations of setting,

1) initialization population

Set up space and interval original state, mathematical description according to relation and the system of operation and coal gas, set the variable data[of property value for solving], whole region of search is divided into to N equally spaced interval, each interval probability of initialization is 1/N, setting population scale is size, the random first generation that generates in search volume;

2) parameter is upgraded, and according to the solution generated, calculates corresponding fitness value g (x _k)

$g\left(x_k\right)=\frac{1}{\underset{i=0}{\overset{2}{\mathrm{\Σ}}}[{\mathrm{\α}}_i{f}_i\left(t\right)+\underset{m=0}{\overset{M-1}{\mathrm{\Σ}}}{\mathrm{con}}_{\mathrm{mi}}\left(t\right)]}$

α _ifor the i class coal gas of correspondence shared weight, f _i(t) the total class that is the mode of communicating between operation and coal gas for the amount of diffusing, the M of i class coal gas is counted summation, con _mi(t) for to have the disappearance amount of the operation of m kind mode of communicating to i class coal gas, and disappearance is measured this part and be can be understood as purchase volume;

According to fitness value calculation distance, renewal interval probability, renewal reference set, define each interval probability update rule: get m the individuality that rank is forward, judge which interval is this m excellent solution drop on respectively on, suppose to have g excellent solution to drop on i interval, i interval probability update rule:

$p\left(i\right)=(1-\mathrm{\α})\×p\left(i\right)+\mathrm{\α}\×\underset{j=1}{\overset{g}{\mathrm{\Σ}}}\frac{g\left(x_g\right)}{\underset{k=1}{\overset{m}{\mathrm{\Σ}}}g\left(x_k\right)}$

P (i) is interval probability.Wherein α is learning rate, and this probability Renewal model can guarantee

$\underset{i=0}{\overset{N-1}{\mathrm{\Σ}}}(1-\mathrm{\α})\×p\left(i\right)+\mathrm{\α}\×\underset{j=1}{\overset{g}{\mathrm{\Σ}}}\frac{g\left(x_g\right)}{\underset{k=1}{\overset{m}{\mathrm{\Σ}}}g\left(x_k\right)}=1;$

3) upgrade reference set

After upgrading each interval probability, get the individual preferably interval of b1 and be combined into excellent disaggregation W, interval distance and definition criterion:

$D\left(i\right)=\underset{j=0}{\overset{N-1}{\mathrm{\Σ}}}|j-i|\×{\mathrm{num}}_j;$

$D\left(j\right)=\underset{j=0}{\overset{N-1}{\mathrm{\Σ}}}|j-i|\×{\mathrm{num}}_i;$

Num _jby j to i process interval number, num _iby i to j the interval number of process, get b2 maximum interval of D (i) and be combined into the set that diversity is good, as the subset D with reference to collection, in order more accurately to search out optimum solution, set an integer M, after the i interval of j the attribute of population is listed in first of excellent solution set W and surpasses M time, for this attribute, the hunting zone of algorithm starts to narrow down to the i interval, initialization operation is again carried out in search to this attribute, use the same method and searched for, and the hunting zone of other dimensions and strategy remain unchanged;

4) generate filial generation

After having upgraded reference set, for choosing between two sub-concentration zones of reference set, to turn method by roulette selected, for the choosing probability and will obtain through conversion of interval, take excellent disaggregation as example here:

$\mathrm{pr}\left(i\right)=\frac{p\left(i\right)}{\underset{j=0}{\overset{b1-1}{\mathrm{\Σ}}}p\left(j\right)}$

The interval choosing method of diversity in gathering well is similar, merges reference set C=D ∪ W;

The definite of filial generation is to realize by the union operation of two subsets, to merge be that intersection, mutation operator by genetic algorithm realized to subset herein, the subset produced is an interval, the numerical range interval according to this generates filial generation, subset merges per generation carries out size time, generates size filial generation.

Compared with prior art, the present invention has the following advantages:

1) dispatching method provided by the invention can reduce the amount of diffusing of coal gas, reduces the coal gas purchase volume.

2) Estimation of Distribution Algorithm that the present invention proposes can converge to optimization solution rapidly, in scheduling process, can not produce hysteresis because of complicated data, has greatly improved the feasibility of this scheduling scheme.

3) two_phase assignment method provided by the invention, minimum to meeting system coal gas diffusion rate from the minimum requirements that meets the operation heat demand, the requirement that purchase volume is minimum.Whole assigning process is incremental, and managerial personnel can carry out decision-making according to the result of each process in the meantime, because the source of model hypothesis coal gas has two kinds, a kind of is to produce by self, and a kind of is to pass through means of purchase.In fact, coal gas can also be by carrying and come from other producing region, but this does not affect the feasibility of this dispatching method, managerial personnel can be according to the purchase volume shown in display, suitably input appropriate coal gas, be that purchase volume in result of calculation of the present invention can be subdivided into two parts, operational throughput and purchase volume.

4) three kinds of component devices provided by the invention complement each other, and predicted data derives from the gas pipe device, utilize core allocation strategy and the algorithm dispatched in optimizer, calculate allocation strategy, finally are shown in display device.Three kinds of equipment, Each performs its own functions, greatly facilitates staff's management.

5) the present invention is through certain iron and steel enterprise's checking, and this device can be realized coke oven, coal gas of converter zero are diffused to (system is actual diffuses), and the amount of diffusing of blast furnace gas is reduced to 87% of the actual amount of diffusing of system.

The accompanying drawing explanation

Fig. 1 is iron and steel enterprise's coal gas network diagram, and this network comprises blast furnace, converter, three kinds of schematic diagram of coke-oven gas, and related 13 kinds of operations in the coal gas network.

Fig. 2 is the calcspar of iron and steel enterprise's coal gas Optimized Operation device, and wherein optimizer is built-in two stage dispatching methods and the improved Estimation of Distribution Algorithm provided in the present invention;

Fig. 3 is improved Estimation of Distribution Algorithm flow process, and this algorithm mainly is comprised of two parts, histogram model and scatter search;

Fig. 4 is iron and steel enterprise's two stages gas dispatching process flow diagram, and the beginning condition is established t=0, and end condition can be the system operation time of setting or the unconditional termination of system.

Embodiment

Below in conjunction with the drawings and specific embodiments, the present invention is described in detail.

Embodiment

At first, need to be to operation in iron and steel enterprise and coal gas be communicated with situation investigated, these operations comprise the hydrogen station mix bypass, coal tar company, hot rolling mill, cold rolling mill, to old liberated area output, hydrogen station, steel-making, heat generator, pelletizing, sintering and pool trouble cogeneration plant, and their operation numbering is started from scratch and is once increased progressively, obtain be communicated with as shown in table 1.

Wherein, the mode of communicating one of operation and coal gas has 4 kinds of UNICOM's modes, it is respectively No. 0 UNICOM's mode, be that operation only has coke-oven gas that heat can be provided, No. 1 UNICOM's mode is that operation can provide heat jointly by coke-oven gas and blast furnace gas, No. 2 UNICOM's modes are that operation can provide heat jointly by coke-oven gas and coal gas of converter, and No. 3 UNICOM's modes are that operation can provide heat demand jointly by three kinds of coal gas.

And, according to operation priority dicision rules, the operation related in the coal gas network is sorted.

Deposit the priority of this operation and operation in the gas pipe device with the situation that is communicated with of coal gas, simultaneously, deposit the historical data of the heat demand of the turnout of three kinds of coal gas and each operation in the gas pipe device.Dope current time and the Gas Production amount in next moment and the heat demand of operation according to these historical datas.Export these predicted values by the gas pipe device and be transported to optimizer.At first optimizer carries out respectively the initialization in two moment according to these predicted values and distributes, and the concrete operations that initialization distributes are as Fig. 3, and initialization is once judged after distributing:

If No. 0 operation of current time is not enough, directly outsourcing, outsourcing amount 0, if No. 6 operations are inadequate, judge whether not through fine setting again.No. 1-5 and the directly outsourcing of 7-10 deficiency.

Then, once finely tune, fine setting is mainly for No. 6 operations, reason is that, except No. 0 and No. 6 operations, all the other operations all can be by blast furnace gas and coke-oven gas service, first to distribute coke-oven gas in principle, the reallocation coal gas of converter, even with distributing blast furnace gas, the heat summation that they provide under the prerequisite that just meets heat demand is consistent.That is,, if not foot phenomenon of coal gas occurs other operations, current time is to go to meet by adjusting gas proportion the coal gas disappearance amount of current operation.And if the situation of coal gas quantity not sufficients occurs in No. 6 operations, probably on the operation high in No. 6 operations of priority ratio, to consume too much, and blast furnace gas just has surplus, can be by carrying out the equivalent exchange of blast furnace gas and coke-oven gas on operation in front, principle of equal value is that calorie value equates, the selected suitable higher operation of priority consumes blast furnace gas of equal value, and this part coke-oven gas is supplied with to No. 6 operations.

The result of distributing according to initialization, coal gas is carried out to secondary distribution, next coal gas is constantly put the coal gas that sky goes to store the surplus in a moment as far as possible, the surplus of will trying one's best of upper one coal gas constantly goes to meet the next deficiency of corresponding coal gas constantly, and the gas chamber in a upper moment will be put the remaining coal gas that sky goes to store next correspondence constantly as far as possible.To meet the explained hereafter constraint of iron and steel enterprise simultaneously in the secondary distribution process.By these transformations, it is the variable relation in optimizer.If variable 22 dimensions, variable is defined as follows:

Data[0], data[1], data[2] and mean respectively the amount of current time blast furnace gas and coke-oven gas equivalent exchange, the amount of blast furnace gas and coal gas of converter equivalent exchange, the amount of coke-oven gas and coal gas of converter equivalent exchange.

Data[3], data[4], data[5], data[6] mean respectively coke-oven gas offer have 0,1,2, the part of consumption again of the operation of 3 kind of UNICOM's mode.

Data[7], data[8], data[9], data[10] mean respectively coke-oven gas offer have 0,1,2, the operation of 3 kind of UNICOM's mode, meet next disappearance amount part constantly.

Data[11], data[12] mean respectively coal gas of converter offer have 2, the part of consumption again of the operation of 3 kind of UNICOM's mode.

Data[13], data[14] mean respectively coal gas of converter offer have 1, the operation of 3 kind of UNICOM's mode, meet next disappearance amount part constantly.

Data[15], data[16] mean respectively blast furnace gas offer have 2, the part of consumption again of the operation of 3 kind of UNICOM's mode.

Data[17], data[18] mean respectively blast furnace gas offer have 1, the operation of 3 kind of UNICOM's mode, meet next disappearance amount part constantly.

Data[19], data[20], data[21] and mean respectively the amount of next blast furnace gas constantly and coke-oven gas equivalent exchange, the amount of blast furnace gas and coal gas of converter equivalent exchange, the amount of coke-oven gas and coal gas of converter equivalent exchange.

According to scheduling model, these variablees have following relation

$\min :{\mathrm{\α}}_i\underset{i=0}{\overset{2}{\mathrm{\Σ}}}{f}_i+\underset{m=0}{\overset{3}{\mathrm{\Σ}}}\underset{k=0}{\overset{\mathrm{end}}{\mathrm{\Σ}}}{{\mathrm{con}}^{*}}_m\left(t\right)$

s.t.data[3]＜＝ma _i(0)

data[4]+data[11]＜＝ma _k(1)

data[5]+data[15]＜＝ma _k(2)

data[6]+data[12]+data[16]＜＝ma _i(3)

data[7]＜＝con _k+1(0)

data[8]+data[17]＜＝con _k+1(1)

data[9]+data[13]＜＝con _k+1(2)

data[10]+data[14]+data[18]＜＝con _k+1(3)

And these data[t], the value that t is 3-18 all is less than the surplus value through coal gas after equivalent exchange.Con ^* _m(t+1) the coal gas amount that while meaning t+1, etching system also lacks after the coal gas of gas chamber supplements through the reserved son of current time, m is 0,1,2,3.

Each variable is constraint mutually, according to operation priority, distributed, often carry out a step distribution, the numerical value restriction relation of 22 dimension variablees will change accordingly, for example coke-oven gas has distributed a part to certain operation, the transformed coke-oven coal of this operation tolerance will become many, and the amount of the assignable coke-oven gas of operation of back will tail off.There is complicated relation between 22 dimension variablees, utilize mathematical method to be solved, can only utilize intelligent algorithm to solve.

The result of trying to achieve is presented on display device, and the staff can be dispatched according to the data that show on display device.

Claims (8)

1. iron and steel enterprise's coal gas Optimized Operation device, this dispatching device is connected with energy management system, it is characterized in that, and described dispatching device comprises gas pipe device, optimizer and the display device connected successively; The gas pipe device utilizes the energy management system of enterprise to obtain each operation gas consumption amount and growing amount historical data and gas chamber data on stock, and coal gas current time and turnout and the consumption in the following moment are predicted; Optimizer, based on above-mentioned predicted data, adopts improved Estimation of Distribution Algorithm to be optimized gas allocation, and the principle of optimization is to make the amount of diffusing of coal gas minimum minimum with purchase volume; The gas allocation scheme that display device display optimization device obtains, for reference, carry out reasonably gas allocation.

2. a kind of iron and steel enterprise according to claim 1 coal gas Optimized Operation device, is characterized in that, described optimizer, based on above-mentioned predicted data, adopts improved Estimation of Distribution Algorithm to be optimized and to be specially gas allocation:

Take coal gas diffusion amount and purchase volume minimum be the final optimization pass target, meeting under the prerequisite of process constraint, operation is carried out to priority definition, then first coal gas is carried out to original allocation, the principle of original allocation is the priority according to operation, take and consume coke-oven gas and coal gas of converter is principle as far as possible, the minimum heat that meets operation of take requires as benchmark, finally coal gas is carried out to secondary distribution, this distribution belongs to distribution according to need, iteration is distributed, and utilizes improved Estimation of Distribution Algorithm to regulate the gas allocation amount of each operation.

3. a kind of iron and steel enterprise according to claim 2 coal gas Optimized Operation device, is characterized in that, described operation is carried out priority definition and is specially:

At first will be by the nonadjustable operation priority ordering of gas consumption structure during Operation Sequencing, it is less important preferential with the operation that consumes coke-oven gas, again preferential with the operation that consumes coal gas of converter.

4. a kind of iron and steel enterprise according to claim 2 coal gas Optimized Operation device, is characterized in that, describedly first coal gas carried out to original allocation and be specially:

The minimum heat demand that meets operation of take is principle, at first distribute coke-oven gas, travel through the operation relevant to coke-oven gas, until coke-oven gas distributes, once coke-oven gas distributes the distribution that has just entered coal gas of converter, the same traversal operation relevant to coal gas of converter, until, after fully coal gas of converter being distributed, finally do same operation to blast furnace gas.

5. a kind of iron and steel enterprise according to claim 2 coal gas Optimized Operation device, is characterized in that, describedly coal gas is carried out to secondary distribution is specially:

1) determine coal gas side's formula case of current time and the relation of next gas chamber storage constantly;

2) one group of variable data[to be solved is set];

3) determine optimization aim, comprise that the amount of diffusing is minimum and the coal gas purchase volume is minimum;

4) adopt improved Estimation of Distribution Algorithm to be distributed, be about to searching process and be divided into two levels, evenly divide whole solution interval, with scatter searching, guide the Evolution of Population direction, the interval produced according to scatter search produces filial generation, with histogram model statistic mass's probability again.

6. a kind of iron and steel enterprise according to claim 5 coal gas Optimized Operation device, is characterized in that, coal gas side's formula case of described definite current time is specially with the relation of next gas chamber storage constantly:

1), if be carved with the coal gas surplus when next, go up a moment and will gas chamber be put to sky as far as possible and go to meet the next storage demand of coal gas constantly;

2) if the next situation of coal gas Shortcomings constantly, current time will reserve corresponding coal gas in gas chamber as far as possible, to meet the next coal gas deficiency situation that constantly may exist.

7. a kind of iron and steel enterprise according to claim 5 coal gas Optimized Operation device, it is characterized in that, described variable to be solved specifically comprises the amount of blast furnace gas and coke-oven gas equivalent exchange, the amount of blast furnace gas and coal gas of converter equivalent exchange, the amount of coke-oven gas and coal gas of converter equivalent exchange, three kinds of coal gas offer respectively the consumption again of the operation with corresponding mode of communicating, three kinds of coal gas offer respectively the disappearance amount of the operation with corresponding mode of communicating and the amount of next blast furnace gas constantly and coke-oven gas equivalent exchange, the amount of blast furnace gas and coal gas of converter equivalent exchange, the amount of coke-oven gas and coal gas of converter equivalent exchange.

8. a kind of iron and steel enterprise according to claim 5 coal gas Optimized Operation device, it is characterized in that, the improved Estimation of Distribution Algorithm of described employing distributes to be specially sets up population space, divide equally spaced interval in space, there is own corresponding probability in each interval, this probability is characterized by the probability that excellent individual occurs, by each, interval probability is evaluated the member in reference set, this algorithm is provided with two reference sets, excellent solution reference set and the good reference set of diversity, pass through the calculating of genetic operator by the interval of selecting at random respectively respectively from two reference sets, generate an interval number, from this interval, the random number that generates forms the next generation again, of future generation through calculating, upgrade again probability, upgrade the reference set member, so constantly repeat until meet termination condition, termination condition is for arriving the number of iterations of setting,

1) initialization population

Set up space and interval original state, mathematical description according to relation and the system of operation and coal gas, set the variable data[of property value for solving], whole region of search is divided into to N equally spaced interval, each interval probability of initialization is 1/N, setting population scale is size, the random first generation that generates in search volume;

2) parameter is upgraded, and according to the solution generated, calculates corresponding fitness value g (x _k)

$g\left(x_k\right)=\frac{1}{\underset{i=0}{\overset{2}{\mathrm{\Σ}}}[{\mathrm{\α}}_i{f}_i\left(t\right)+\underset{m=0}{\overset{M-1}{\mathrm{\Σ}}}{\mathrm{con}}_{\mathrm{mi}}\left(t\right)]}$

α _ifor the i class coal gas of correspondence shared weight, f _i(t) the total class that is the mode of communicating between operation and coal gas for the amount of diffusing, the M of i class coal gas is counted summation, con _mi(t) for to have the disappearance amount of the operation of m kind mode of communicating to i class coal gas, and disappearance is measured this part and be can be understood as purchase volume;

According to fitness value calculation distance, renewal interval probability, renewal reference set, define each interval probability update rule: get m the individuality that rank is forward, judge which interval is this m excellent solution drop on respectively on, suppose to have g excellent solution to drop on i interval, i interval probability update rule:

$p\left(i\right)=(1-\mathrm{\α})\×p\left(i\right)+\mathrm{\α}\×\underset{j=1}{\overset{g}{\mathrm{\Σ}}}\frac{g\left(x_g\right)}{\underset{k=1}{\overset{m}{\mathrm{\Σ}}}g\left(x_k\right)}$

P (i) is interval probability.Wherein α is learning rate, and this probability Renewal model can guarantee

$\underset{i=0}{\overset{N-1}{\mathrm{\Σ}}}(1-\mathrm{\α})\×p\left(i\right)+\mathrm{\α}\×\underset{j=1}{\overset{g}{\mathrm{\Σ}}}\frac{g\left(x_g\right)}{\underset{k=1}{\overset{m}{\mathrm{\Σ}}}g\left(x_k\right)}=1;$

3) upgrade reference set

After upgrading each interval probability, get the individual preferably interval of b1 and be combined into excellent disaggregation W, interval distance and definition criterion:

$D\left(i\right)=\underset{j=0}{\overset{N-1}{\mathrm{\Σ}}}|j-i|\×{\mathrm{num}}_j;$

$D\left(j\right)=\underset{j=0}{\overset{N-1}{\mathrm{\Σ}}}|j-i|\×{\mathrm{num}}_i;$

Num _jby j to i process interval number, num _iby i to j the interval number of process, get b2 maximum interval of D (i) and be combined into the set that diversity is good, as the subset D with reference to collection, in order more accurately to search out optimum solution, set an integer M, after the i interval of j the attribute of population is listed in first of excellent solution set W and surpasses M time, for this attribute, the hunting zone of algorithm starts to narrow down to the i interval, initialization operation is again carried out in search to this attribute, use the same method and searched for, and the hunting zone of other dimensions and strategy remain unchanged;

4) generate filial generation

After having upgraded reference set, for choosing between two sub-concentration zones of reference set, to turn method by roulette selected, for the choosing probability and will obtain through conversion of interval, take excellent disaggregation as example here:

$\mathrm{pr}\left(i\right)=\frac{p\left(i\right)}{\underset{j=0}{\overset{b1-1}{\mathrm{\Σ}}}p\left(j\right)}$

The interval choosing method of diversity in gathering well is similar, merges reference set C=D ∪ W;

The definite of filial generation is to realize by the union operation of two subsets, to merge be that intersection, mutation operator by genetic algorithm realized to subset herein, the subset produced is an interval, the numerical range interval according to this generates filial generation, subset merges per generation carries out size time, generates size filial generation.

Images