CN104375477A - Production method based on integrated optimization of production system and energy system - Google Patents

Abstract

The invention discloses a production method based on integrated optimization of a production system and an energy system. The method includes the steps of establishing a multi-period production plan module of the production system an oil refining enterprise, a multi-period operation optimization model of the energy system of the oil refining enterprise and a material flow and energy flow multi-period coupling model of the production system and the energy system, establishing a material flow and energy flow multi-period integrated optimization model according to a material cost and energy cost integrated objective function, achieving rapid optimization search of the material flow and energy flow multi-period coupling model through the model decomposition solution strategy, and conducting production according to the result of the rapid optimization search. The method can effectively solve the problems that due to the fact that the production system and the energy system of a traditional oil refining enterprise are interpedently optimized step by step, the energy medium waste is caused, and the energy supply and demand are not balanced, and by achieving the material flow and energy flow integrated optimization, the energy conservation and emission reduction can be achieved while the production profits of enterprises are improved.

Description

A kind of production method of the integrated optimization based on production system and energy system

Technical field

The present invention relates to Petrochemical Production Technologies field, be specifically related to a kind of production method of the integrated optimization based on production system and energy system.

Background technology

The production schedule is taking into full account on the basis that enterprises' production capacity and the market demand and the utilization of resources balance, to enterprise produce the kind of product, comprehensive overall arrangement that quantity, the aspect such as quality and manufacturing schedule are done, be blueprint for action and the administration base of enterprise's activity in production.The production schedule occupy key link in activity in production, can coordinate production and market, helps enterprise to make full use of resource, realizes the management objectives of enterprise.

Oil refining enterprise's (petrochemical enterprise) process units is many, the production cycle is long, product signal is many, logistics is complicated, and the Commodity flow in enterprise and energy flow influence each other.Under the condition considering changes in market demand, to the Commodity flow of oil refining enterprise with energy flow is analyzed and global optimization, formulate the production schedule and be very important.Production planning optimization is as an important step in modern oil refining enterprise production run, all the time all to optimize full factory flow distribution for target, namely under given crude oil and product supply/demand, based on plant capacity in production system, optimize the Product processing scheme (i.e. operation scheme) determining full factory device, realize full factory productive profit by resource distribution the highest.

Traditional oil refining enterprise comprises the production system based on processing materials, and in process of production for production system provides the energy system of different public work.Energy system all the time all as the subsystem of process of producing product, for materiel machining device provides energy demand.In the production run of oil refining enterprise, by production run and public work integrally, setting up the Production planning model of full factory, by the system optimization of the production schedule and public work, is the important means that oil refining enterprise saves energy and reduce the cost and increases the benefit.

Owing to realizing Commodity flow that enterprise-level produces and the integrated optimization of energy flow, need production system and the energy system integrated moulding of just enterprise, the features such as the scale obtained due to integrated moulding is large, nonlinear relationship is complicated, thus increase the difficulty that conventional linear planing method solves model optimization, existing oil refining enterprise is when formulating the production schedule, usually respectively to production system and energy system Modling model, and substep independent optimization is carried out.Although improve economic benefit to a certain extent, be limited to and do not set up oil refining enterprise's logistics and can coupling model be flowed, thus lose the effective room for promotion of full factory economic benefit, cannot accomplish that the high efficiente callback of the energy utilizes simultaneously.

Summary of the invention

For the deficiencies in the prior art, the invention provides a kind of production method of the integrated optimization based on production system and energy system.

Based on a production method for the integrated optimization of production system and energy system, comprising:

(1) the multicycle Production planning model setting up production system is combined according to oil refining enterprise's production procedure and device;

(2) set up the multicycle operation optimization model of energy system, comprise multi-state boiler model, multi-state steam turbine model, the constraint of each energy medium equilibrium of supply and demand and equipment capacity and retrain;

(3) utilize the product power consumption model of the material balance relationship of intermediate product, materiel machining device and be in harmonious proportion physical property and the public work gas discharging factor correlation model of product, set up the logistics of production system and energy system and multicycle coupling model can be flowed;

(4) model set up according to step (1), (2), (3) and the integrated objective function according to material and energy cost, set up logistics and multicycle Integrated Optimization Model can be flowed, utilize model decomposition solution strategies to logistics with can flow multicycle coupling model and carry out quick optimizing;

(5) produce according to the result of the quick optimizing in step (4).

The present invention is based on the existing production system planning model in oil refining enterprise and energy system planning model, innovation proposes a kind of modeling method being applied to oil refining enterprise's production system and energy system integrated optimization, complete logistics (Commodity flow) and obtain Integrated Optimization Model with the Model coupling that can flow (energy flow), and utilize decomposition measurement to solve Integrated Optimization Model, instruct oil refining to produce according to solving result, full factory productive profit can have been promoted and realize energy-saving and emission-reduction.

Multicycle Production planning model in described step (1) comprises the yield model of product demand constraint in production system, material balance, stock's constraint, harmonic process constraint and each materiel machining device.

Production system multicycle production planning optimization model, combines according to the production procedure of actual oil refining enterprise and device and builds.Wherein, product demand constraint representation is as follows:

${\mathrm{DP}}_{c,t}^{\mathrm{up}}\≥{\mathrm{SC}}_{c,t}\≥{\mathrm{DP}}_{c,t}^{\mathrm{lo}},---\left(1\right)$

with the upper and lower bound of the market demand of product c in difference indication cycle t, SC _c,tthe sales volume of product c in indication cycle t.

Material balance is expressed as follows:

$\underset{c\∈{\mathrm{CI}}_{u,m}}{\mathrm{\Σ}}{\mathrm{FC}}_{u,m,c,t}=\underset{c\∈{\mathrm{CO}}_{u,m}}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t},---\left(2\right)$

with distinguish input amount and the quantum of output of materiel machining assembling device u material c under operation scheme m in indication cycle t.

Stock's constraint representation is as follows:

${\mathrm{MI}}_{c,t}^{\mathrm{lo}}\≤{\mathrm{MI}}_{c,t}\≤{\mathrm{MI}}_{c,t}^{\mathrm{up}},---\left(3\right)$

${\mathrm{MI}}_{c,t}={\mathrm{MI}}_{c,t-1}+\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m\∈\mathrm{MU}}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}-\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m\∈\mathrm{MU}}{\mathrm{\Σ}}{\mathrm{FC}}_{u,m,c,t},---\left(4\right)$

with distinguish lower limit and the upper limit of the inventory limitation of product c in indication cycle t, MI _c,tthe value of inventory of product c in indication cycle t.

Harmonic process constraint representation is as follows:

$\underset{c\∈{\mathrm{CI}}_{u,m}}{\mathrm{\Σ}}{\mathrm{PI}}_{c,p}{\mathrm{FC}}_{u,m,c,t}\≥\underset{c\∈{\mathrm{CO}}_{u,m}}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}{\mathrm{PL}}_{c,p},---\left(5\right)$

$\underset{c\∈{\mathrm{CI}}_{u,m}}{\mathrm{\Σ}}{\mathrm{PI}}_{c,p}{\mathrm{FC}}_{u,m,c,t}\≤\underset{c\∈{\mathrm{CO}}_{u,m}}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}{\mathrm{PU}}_{c,p},---\left(6\right)$

PI _c,prepresent intermediate product c enter attemperation apparatus before physical property p, physical property p comprises the mass propertys such as the sulfur content of product, carbon content, nitrogen content and condensation point, PL _c,pwith PU _c,prepresent lower limit and the upper limit of the physical property p value of the final mediation product (i.e. final products) of attemperation apparatus output respectively.Wherein, intermediate product refer to the outputting material of materiel machining device, and as the input material of attemperation apparatus.

The yield model of materiel machining device comprises the yield model of each materiel machining device, inputoutput balance model, operation scheme restricted model and working ability restricted model.

The yield model of each materiel machining device refers to that each product of the side line of each processing unit (plant) output under each operation scheme is that processing capacity is multiplied by side line productive rate under scheme, is specifically expressed as follows:

FP _u,m,c,t＝α _u,m,cFF _u,m,t(7)

In formula, c is the kind (i.e. the kind of the product of materiel machining device) of the material of materiel machining device, the kind of the operation scheme that m uses for materiel machining device, u is the kind of materiel machining device, and t is production cycle (i.e. the cycle), FP _{u, m, c, t}for the output of the product c of materiel machining device u under operation scheme m, α _{u, m, c}for the productive rate of the sideline product c of device u under operation scheme m (i.e. product c), FF _{u, m, t}for the processing capacity of materiel machining device u under operation scheme m in cycle t.

Inputoutput balance model for each materiel machining device refers to that the output of production of materiel machining device under various operation scheme equals the consumption summation of a raw material, is expressed as follows:

${\mathrm{FF}}_{u,m,t}=\underset{c}{\mathrm{\Σ}}{\mathrm{FC}}_{u,m,c,t},---\left(8\right)$

Wherein, FC _{u, m, c, t}for the consumption of the spent material c of materiel machining device u under operation scheme m in cycle t.

For the operation scheme restricted model that each materiel machining device limits, be expressed as follows:

$\underset{m}{\mathrm{\Σ}}{x}_{u,m,t}=1,---\left(9\right)$

Wherein, MU is the set of operation scheme.X _{u, m, t}for 0-1 variable, represent whether materiel machining device u adopts operation scheme m within the t cycle, according to then x _{u, m, t}be 1, do not adopt then x _{u, m, t}be 0.

Working ability constraint for each materiel machining device refers to that the device output of each materiel machining device under specific operation scheme m must meet the restriction of device working ability, is expressed as follows:

$x_{u,m,t}{\mathrm{FU}}_{u,t}^{\mathrm{lo}}\≤{\mathrm{FF}}_{u,m,t}\≤x_{u,m,t}{\mathrm{FU}}_{u,t}^{\mathrm{up}},---\left(10\right)$

In formula, with represent lower limit and the upper limit of the production and processing ability of materiel machining device u in cycle t respectively.

Multi-state boiler model in described step (2) comprises the Fuel Consumption of boiler and the multi-state state constraint of steam production model and boiler, wherein:

The Fuel Consumption of boiler and steam production model, as follows:

${\mathrm{FB}}_{u,c,t}=\underset{q}{\mathrm{\Σ}}(A_{q,u,c,t}{\mathrm{XFB}}_{q-1,u,c,t}+{\mathrm{\β}}_{q,u,c,t}({\mathrm{XFB}}_{q,u,c,t}-{\mathrm{XFB}}_{q-1,u,c,t})),---\left(11\right)$

${\mathrm{SB}}_{u,c,t}=\underset{q}{\mathrm{\Σ}}(A_{q,u,c,t}{\mathrm{XSB}}_{q-1,u,c,t}+{\mathrm{\β}}_{q,u,c,t}({\mathrm{XSB}}_{q,u,c,t}-{\mathrm{XSB}}_{q-1,u,c,t})),---\left(12\right)$

FB _{u, c, t}represent the consumption of boiler u at cycle t internal consumption fuel c, SB _{u, c, t}represent the steam production that boiler u obtains at cycle t internal consumption fuel c, q is operating condition, A _{q, u, c, t}for 0-1 variable, represent whether boiler u operates (0 expression does not operate in operating mode q interval, and 1 expression operates in operating mode q interval) at cycle t in operating mode q interval, β _{q, u, c, t}for boiler u in operating mode q interval is in the linearizing continuous variable of running status of cycle t internal consumption fuel c, value is at 0 ~ 1, XFB _{q, u, c, t}and XFB _{q-1, u, c, t}represent respectively in operating mode q interval with the maximal value of boiler u in operating mode q-1 interval at cycle t internal consumption fuel c, XSB _{q, u, c, t}and XSB _{q-1, u, c, t}to be illustrated respectively in operating mode q interval and the maximal value (namely by piece wire approximation method, making the charging consumption of boiler and boiler gas production rate carry out multistage linear association process) of the steam production that boiler u obtains at cycle t internal consumption fuel c in operating mode q-1 interval;

The multi-state state constraint of boiler, is expressed as follows:

$\underset{q}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{A}_{q,u,c,t}\≤1,---\left(13\right)$

By multi-state state constraint, make boiler u can only consume a kind of fuel at the most in one-period t, and can only run in the interval q of an operating mode.

Multi-state steam turbine model in described step (2) comprises the steam consumption of steam turbine and the multi-state state constraint of generated energy model and steam turbine, wherein:

Make the steam consumption of steam turbine and steam turbine power generation amount carry out multistage linear association process by piece wire approximation method, obtain steam consumption and the generated energy model of steam turbine, embody as follows:

${\mathrm{STC}}_{u,t}=\underset{q\∈Q}{\mathrm{\Σ}}(B_{q,u,t}{\mathrm{XST}}_{q-1,u,t}+{\mathrm{\γ}}_{q,u,t}({\mathrm{XST}}_{q,u,t}-{\mathrm{XST}}_{q-1,u,t}))---\left(14\right)$

${\mathrm{ET}}_{u,t}=\underset{q\∈Q}{\mathrm{\Σ}}(B_{q,u,t}{\mathrm{XET}}_{q-1,u,t}+{\mathrm{\γ}}_{q,u,t}({\mathrm{XET}}_{q,u,t}-{\mathrm{XET}}_{q-1,u,t}))---\left(15\right)$

STC _u,trepresent the steam consumption of steam turbine u in cycle t, ET _u,trepresent the generated energy of steam turbine u in cycle t, q is the operating mode of steam turbine operation, and Q is the set of the operating mode of steam turbine operation, B _{q, u, t}for 0-1 variable, represent whether steam turbine u operates at cycle t in operating mode q interval, γ _{q, u, t}represent the linearizing continuous variable of the running status of steam turbine u in cycle t in operating mode q interval, value is at 0 ~ 1, XST _{q, u, t}and XST _{q-1, u, t}be illustrated respectively in the maximal value of steam turbine u steam consumption in cycle t in operating mode q interval and q-1 interval, XET _{q, u, t}and XET _{q-1, u, t}be illustrated respectively in the maximal value of steam turbine u generated energy in cycle t in operating mode q interval and q-1 interval;

The multi-state state constraint of steam turbine is expressed as follows:

$\underset{q\∈Q}{\mathrm{\Σ}}{B}_{q,u,t}\≤1.---\left(16\right)$

By the multi-state state constraint of steam turbine, steam turbine u can only be run in one-period t in the interval q of an operating mode at the most.

Each energy medium equilibrium of supply and demand constraint in described step (2) comprises the constraint of the steam equilibrium of supply and demand and the constraint of the electric energy equilibrium of supply and demand, wherein:

Steam equilibrium of supply and demand constraint representation is:

$\begin{array}{c}\underset{u\∈\mathrm{UB}}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{\mathrm{SB}}_{u,c,t}+\underset{u\∈\mathrm{UT}}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{\mathrm{STG}}_{u,c,t}+\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{EG}}_{u,m,c,t}-\underset{u\∈\mathrm{UT}}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{\mathrm{STC}}_{u,c,t}+{\mathrm{EP}}_{c,t}\\ +{\mathrm{LSI}}_{c,t}-{\mathrm{LSO}}_{c,t}\≥\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{ED}}_{u,m,c,t}\end{array},---\left(17\right)$

UB is boiler set, and UT is steam turbine set, and UP is materiel machining device sets; SB _{u, c, t}represent the steam production that boiler u obtains at cycle t internal consumption fuel c, STG _{u, c, t}represent the steam production that steam turbine u obtains at cycle t internal consumption fuel c, EG _{u, m, c, t}and ED _{u, m, c, t}distinguish output sum and the consumption sum of materiel machining device u each grade steam of steam c under operation scheme m in indication cycle t, STC _{u, c, t}represent that steam turbine u generates electricity in cycle t and consume the consumption of steam c, EP _c,trepresent the outsourcing amount of steam c in cycle t, LSI _c,tand LSO _c,trepresent the air inflow of the steam c of temperature-decreased pressure reducer in cycle t and air output respectively, each grade steam is respectively 3.5MPa high pressure steam, 1.5MPa middle pressure steam, 0.3MPa low-pressure steam;

Electric energy equilibrium of supply and demand constraint representation is:

${\mathrm{EP}}_{c,t}+\underset{u\∈\mathrm{UT}}{\mathrm{\Σ}}{\mathrm{ET}}_{u,t}\≥\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{ED}}_{u,m,c,t},---\left(18\right)$

EP _c,trepresent the outsourcing amount of steam c in cycle t, ET _u,trepresent the generated energy of steam turbine u in cycle t, ED _{u, m, c, t}represent materiel machining device u consumption of each grade steam of steam c under operation scheme m in cycle t.

Equipment capacity constraint representation in described step (2) is:

${\mathrm{FU}}_u^{\mathrm{lo}}\≤{\mathrm{FU}}_{u,t}\≤{\mathrm{FU}}_u^{\mathrm{up}},---\left(19\right)$

${\mathrm{LS}}_{c,t}\≤{\mathrm{LS}}_{c,t}^{\mathrm{up}},---\left(20\right)$

with the lower limit of the production capacity of difference indication equipment u and the upper limit, FU _u,tthe load of indication equipment u in cycle t, LS _c,tin indication cycle t, temperature-decreased pressure reducer is by the load of the process of steam c, in indication cycle t, temperature-decreased pressure reducer is by the upper limit of the load of the process of steam c.Equipment specifically refers to power plant, the various kinds of equipment namely in energy system, comprises boiler and steam turbine, wherein boiler for producing steam, steam turbine power generation.

Described step (3) comprises the steps:

(3-1) fuel is set up to the mass balance model comprising product inventory, sale, outsourcing, mediation production and boiler and consume, as follows:

${\mathrm{MI}}_{c,t}+\underset{u}{\mathrm{\Σ}}{\mathrm{FB}}_{u,c,t}+{\mathrm{SC}}_{c,t}={\mathrm{MI}}_{c,t-1}+\underset{u}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}+{\mathrm{PC}}_{c,t},---\left(21\right)$

Wherein, MI _c,tand MI _{c, t-1}difference indication cycle t is interior and cycle t-1 fuel is oily or the tank farm stock of gas c, SC _c,tthe sales volume of indication cycle t fuel c, FB _{u, c, t}represent that in energy system, boiler u is at the consumption of cycle t internal consumption fuel c, FP _{u, m, c, t}the output of the fuel c that attemperation apparatus u obtains under operation scheme m in indication cycle t, PC _c,tthe outsourcing amount of indication cycle t fuel c;

(3-2) for materiel machining device, produce by the product consumption of public work in related product production run and the fixing of the load of materiel machining device, the physical property of product and materiel machining device the product of setting up materiel machining device of consuming energy and to consume energy model, as follows:

ED _u,m,c,t＝λ _u,m,cFF _u,m,t+ρ _u,c,pP _c,p+CD _u,m,c， (22)

EG _u,m,c,t＝μ _u,m,cFF _u,m,t+σ _u,c,pP _c,p+CG _u,m,c， (23)

In formula, ED _{u, m, c, t}represent materiel machining device u in cycle t under operation scheme m to the demand of public work c, λ _{u, m, c}represent the materiel machining materiel machining amount of device u under operation scheme m and this consumption ratio coefficient of public work c, FF _{u, m, t}the material treatment capacity of materiel machining device u under operation scheme m in indication cycle t, ρ _{u, c, p}represent the correlation coefficient of the consumption of public work c and the physical property p of product discharge in materiel machining device u, P _c,prepresent the value of the physical property p of product c, CD _{u, m, c}represent that materiel machining device u consumes the fixed charges of public work c under operation scheme m, EG _{u, m, c, t}the output of the public work c that materiel machining device u generates under operation scheme m in indication cycle t, μ _{u, m, c}represent the materiel machining materiel machining amount of device u under operation scheme m and the production ratio coefficient of public work c, σ _{u, c, p}represent the correlation coefficient of the output of public work c and the physical property p of product discharge in materiel machining device u, CG _{u, m, c}the fixing output of indication device u public work c under operation scheme m;

(3-3) for harmonic process and boiler combustion gas generation process, by physical property and the public work gas discharging factor of association mediation product, the correlation model of fuel oil production run and Pollution of Boiler gas discharge process is set up, as follows;

$\underset{c}{\mathrm{\Σ}}{\mathrm{FC}}_{u,m,c,t}{\mathrm{PI}}_{c,p}=\underset{c}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}{P}_{c,p}---\left(24\right)$

XC _{u，c，g，t}＝ω _gf(P _c，p)FB _u，c，t(25)

In formula, PI _c,prepresent the value of the physical property p of fuel oil c before entering attemperation apparatus, XC _{u, c, g, t}represent the discharge capacity of the dusty gas g that the consume fuel oil c of boiler u in cycle t produces, ω _grepresent the emission factor of dusty gas g, f (P _c,p) represent the quality of dusty gas g.

In step (3-1), fuel oil and gas can, as the fuel input product of oil refining enterprise's energy system, be again the intermediate product that device incidentally produces simultaneously.

F (P in step (3-3) _c,p) represent by the function of the exploitation dusty gas quality of the physical property p of fuel oil c.

Described step (4) comprises the steps:

(4-1) physical property of final products is set as fixed value, the multicycle Production planning model of the production system of oil refining enterprise is relaxed into MIXED INTEGER linear model;

(4-2) the input product market demand, and turn to target with full factory material production maximum profit, rapid solving is carried out to the MIXED INTEGER linear model of production system, obtains the load of each materiel machining device in production system, operation scheme, product library storage, the physical property of product and the turnout of gas and fuel oil;

(4-3) result of described step (4-2) is substituted into the product power consumption model of materiel machining device, calculate the product consumption demand of all kinds of energy medium in production system;

(4-4) consumption demand is produced as constraint of demand using all kinds of energy medium of production system, be minimised as target with the energy cost of energy system, to energy system multicycle operation optimization model solve the product consumption of load and all kinds of energy medium obtaining each equipment in energy system;

(4-5) by the fixed value of the physical property of final products, by the turnout of the physical property of the load of materiel machining device each in production system, operation scheme, product library storage, product, gas and fuel oil, and in energy system the load of each equipment and the product consumption of all kinds of energy medium as logistics and the initial feasible solution that can flow multicycle coupling model, utilize branch-and-bound optimized algorithm, solve this logistics and multicycle coupling model can be flowed obtain the optimum material production scheme of whole oil refining enterprise and energy production consumes scheme.

Integrated objective function variable described in step (4) comprises full factory profit, profit of sales, device processing cost, purchasing of raw materials cost, stock's operation cost, fuel consumption cost, Exogenous factor medium cost, equipment operation maintenance expense and environmental emission cost.Equipment wherein in energy system refers to the energy device in energy system, comprises steam turbine and boiler.

Production planning optimization in current enterprise is usually based on production system logistic optmum, the energy requirements obtained based on logistic optmum carries out the optimization of energy system capacity planning, and such step-by-step optimization method have compressed full factory's profit improvement space and is unfavorable for that high efficiency of energy utilizes.The present invention is based on the Coupling method that plant energy consumption model, the energy medium equilibrium of supply and demand and gas discharging coupling model realize oil refining enterprise's logistics and can flow, model decomposition strategy is utilized to realize the rapid solving of integrated model, compared to traditional step-by-step optimization method, not only effectively improve the productive profit of full factory (oil refining enterprise), realize the energy-saving and emission-reduction that enterprise produces simultaneously.The present invention is portable strong, goes for different manufacturing enterprise's environment.

Accompanying drawing explanation

Fig. 1 is the production logistics system process flow diagram of corresponding oil refining enterprise in the present embodiment;

Fig. 2 is the energy system process flow diagram of corresponding oil refining enterprise in the present embodiment;

Fig. 3 is the operational efficiency curve map of energy system five boiler multistage linear;

Fig. 4 is each energy medium relation between supply and demand interaction figure between production system and energy system;

Fig. 5 is the model decomposition solution strategies schematic diagram proposed in the present invention.

Embodiment

In order to more specifically describe the present invention, below in conjunction with embodiment and accompanying drawing, the present invention is described in detail.

Based on a production method for the integrated optimization of production system and energy system, comprise the following steps:

(1) the multicycle Production planning model of described oil refining enterprise production system is set up, comprising product demand constraint, material balance, stock's constraint, harmonic process constraint and materiel machining device yield model.Fig. 1 is the flow sheet of this refinery, and wherein solid line represents liquid stream route, and dotted line represents gaseous stream route, and wherein raw material comprises crude oil and MTBE (methyl tert-butyl ether).

The multicycle Production planning model product demand constraint of oil refining enterprise's production system, is expressed as follows:

${\mathrm{DP}}_{c,t}^{\mathrm{up}}\≥{\mathrm{SC}}_{c,t}\≥{\mathrm{DP}}_{c,t}^{\mathrm{lo}},---\left(1\right)$

with the upper and lower bound of the market demand of product c in difference indication cycle t, SC _c,tthe sales volume of product c in indication cycle t.

Product (final products) in the present embodiment comprising: gas Fuel gas, 93# gasoline Gasoline, 90# gasoline Gasoline, kerosene Kerosene ,-10# diesel oil Diesel, 0# diesel oil Diesel and fuel oil Fuel oil;

Material balance constraint representation is as follows:

$\underset{c\∈{\mathrm{CI}}_{u,m}}{\mathrm{\Σ}}{\mathrm{FC}}_{u,m,c,t}=\underset{c\∈{\mathrm{CO}}_{u,m}}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t},---\left(2\right)$

FC _{u, m, c, t}and FP _{u, m, c, t}distinguish input amount and the quantum of output of materiel machining device u material c under operation scheme m in indication cycle t.

Stock's constraint representation is as follows:

${\mathrm{MI}}_{c,t}^{\mathrm{lo}}\≤{\mathrm{MI}}_{c,t}\≤{\mathrm{MI}}_{c,t}^{\mathrm{up}},---\left(3\right)$

${\mathrm{MI}}_{c,t}={\mathrm{MI}}_{c,t-1}+\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m\∈\mathrm{MU}}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}-\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m\∈\mathrm{MU}}{\mathrm{\Σ}}{\mathrm{FC}}_{u,m,c,t},---\left(4\right)$

MU is the set of operation scheme, and MP is the set of materiel machining device, with distinguish lower limit and the upper limit of the inventory limitation of product c in cycle t, MI _c,tthe value of inventory of product c in indication cycle t; Formula (4) represents that this cycle stock of product c equals cycle stock and the output sum of this cycle belongings material processing unit (plant) to this product deducts all material processing unit (plant)s to the consumption of this product.

Intermediate product refer to the outputting material of materiel machining device, and as the input material of attemperation apparatus.Intermediate product comprise: light naphthar LSR, heavy naphtha HSR, rough coal oil RKER, light gas oil LGO, atmospheric gas oil AGO, vacuum gas oil VGO, vacuum residuum RED, Reformed Gasoline RG, catalytically cracked gasoline CG, catalytic cracking diesel oil CGO, thick gas SFG, low-sulfur gas LSFG, kerosene KER, diesel oil DIE, hydrogenated fuel oil FOH, crackedoil FOF.

Harmonic process constraint representation is as follows:

$\underset{c\∈{\mathrm{CI}}_{u,m}}{\mathrm{\Σ}}{\mathrm{PI}}_{c,p}{\mathrm{FC}}_{u,m,c,t}\≥\underset{c\∈{\mathrm{CO}}_{u,m}}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}{\mathrm{PL}}_{c,p},---\left(5\right)$

$\underset{c\∈{\mathrm{CI}}_{u,m}}{\mathrm{\Σ}}{\mathrm{PI}}_{c,p}{\mathrm{FC}}_{u,m,c,t}\≤\underset{c\∈{\mathrm{CO}}_{u,m}}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}{\mathrm{PU}}_{c,p},---\left(6\right)$

PI _c,prepresent intermediate product c enter attemperation apparatus before physical property p, physical property p comprises the mass propertys such as the sulfur content of product, carbon content, nitrogen content and condensation point, PL _c,pwith PU _c,prepresent lower limit and the upper limit of the physical property p value of the final products of attemperation apparatus output respectively.

Attemperation apparatus specifically comprises: gas mixing arrangement (gas tank), 93# gasoline blending device (Gasoline blender), 90# gasoline blending device (Gasoline blender), kerosene mixing arrangement (Kerosene tank) ,-10# diesel oil blending device (Diesel blender), 0# diesel oil blending device (Diesel blender) and fuel oil attemperation apparatus (Fuel oil).

The yield model of materiel machining device comprises yield model, inputoutput balance model, operation scheme restricted model, the working ability restricted model of each materiel machining device.

Materiel machining device specifically comprises: atmospheric and vacuum distillation unit CDU, catalytic cracking unit FCC, catalytic reforming unit CRU, delayed coking unit DCU, hydro-refining unit HT, hydrogenation desulfurizing device HDS and gas desulfurizing device DS.

The yield model of each materiel machining device refers to that each product of the side line of this materiel machining device output under each operation scheme is the side line productive rate under processing capacity is multiplied by corresponding operating scheme, is specifically expressed as follows:

FP _u,m,c,t＝α _u,m,cFF _u,m,t， (7)

In formula, c is the kind (i.e. the kind of the product of materiel machining device) of the material of materiel machining device, the kind of the operation scheme that m uses for materiel machining device, and u is the kind of materiel machining device, t is production cycle (i.e. the cycle), FP _{u, m, c, t}for the output of the product c of materiel machining device u under operation scheme m, α _{u, m, c}for the productive rate of the sideline product c of device u under operation scheme m (i.e. product c), FF _{u, m, t}for the processing capacity of materiel machining device u under operation scheme m in cycle t.

Inputoutput balance model for each materiel machining device refers to that this materiel machining device output of production under each operation scheme equals, to each consumption of raw materials amount summation, to be specifically expressed as follows:

${\mathrm{FF}}_{u,m,t}=\underset{c}{\mathrm{\Σ}}{\mathrm{FC}}_{u,m,c,t},---\left(8\right)$

Wherein, FC _{u, m, c, t}for the consumption of the spent material c of materiel machining device u under operation scheme m in cycle t.

The operation scheme restricted model that each materiel machining device limits, is expressed as follows:

$\underset{m}{\mathrm{\Σ}}{x}_{u,m,t}=1,---\left(9\right)$

Wherein, x _{u, m, t}for 0-1 variable, represent whether materiel machining device u adopts operation scheme m, according to x in the t cycle _{u, m, t}be 1, do not adopt x _{u, m, t}be 0.

Working ability for each materiel machining device retrains, and is expressed as follows:

$x_{u,m,t}{\mathrm{FU}}_{u,t}^{\mathrm{lo}}\≤{\mathrm{FF}}_{u,m,t}\≤x_{u,m,t}{\mathrm{FU}}_{u,t}^{\mathrm{up}},---\left(10\right)$

Namely the materiel machining device output of each materiel machining device under operation scheme m must meet the working ability restriction of device.In formula, with represent lower limit and the upper limit of the production and processing ability of each materiel machining device u in cycle t respectively.

(2) set up the multicycle operation optimization model of the energy system of the oil refining enterprise of the present embodiment, the multicycle operation optimization model of energy system comprises multi-state boiler model, multi-state steam turbine model, each energy medium equilibrium of supply and demand and equipment capacity and retrains.The energy system process flow diagram of Tu2Shi oil refining enterprise, energy system consume fuel (comprising fuel oil, gas and rock gas) produces steam and electricity, for production system and electrical network provide the energy.The energy system of the present embodiment comprises five boiler (Bl1, Bl2, Bl3, Bl4, Bl5) with six steam turbines (T1, T2, T3, T4, T5, T6), wherein Bl1 and Bl2 produces high pressure steam, Bl3, Bl4 and Bl5 produce middle pressure steam, boiler fuel consumption can select fuel oil, gas or rock gas, the generating of steam turbine consumption gas, and T1 consumes high pressure steam and discharges low-pressure steam, T2 consumes high pressure steam and discharges middle pressure steam, and T3, T4, T5 and T6 consume middle pressure steam and discharge low-pressure steam.The steam producing different brackets is incorporated to steam pipe system unified distribution, and generating is unified to be connected to the grid.The present embodiment mesohigh steam is 3.5MPa high pressure steam, middle pressure steam is 1.5MPa, low-pressure steam is 0.3MPa.In addition, the energy system of the present embodiment is also provided with reducing-and-cooling plant and oxygen-eliminating device, by reducing-and-cooling plant steam regulation pressure grade, is also controlled the combustion position of fuel further by deaerating plant, auxiliary adjustment vapour pressure grade.

Multi-state boiler model comprises the Fuel Consumption of boiler and the multi-state state constraint model of steam production model and boiler.The operational efficiency curve map (representing the relation of boiler feed ratio and steam production ratio) of the multistage linear of five boilers in Tu3Shi oil refining enterprise energy system.

The Fuel Consumption of boiler and steam production model, be expressed as follows respectively:

${\mathrm{FB}}_{u,c,t}=\underset{q}{\mathrm{\Σ}}(A_{q,u,c,t}{\mathrm{XFB}}_{q-1,u,c,t}+{\mathrm{\β}}_{q,u,c,t}({\mathrm{XFB}}_{q,u,c,t}-{\mathrm{XFB}}_{q-1,u,c,t})),---\left(11\right)$

${\mathrm{SB}}_{u,c,t}=\underset{q}{\mathrm{\Σ}}(A_{q,u,c,t}{\mathrm{XSB}}_{q-1,u,c,t}+{\mathrm{\β}}_{q,u,c,t}({\mathrm{XSB}}_{q,u,c,t}-{\mathrm{XSB}}_{q-1,u,c,t})),---\left(12\right)$

Namely by piece wire approximation method, the charging consumption of boiler and boiler gas production rate is made to carry out multistage linear association process.Wherein, FB _{u, c, t}represent the consumption of boiler u at cycle t internal consumption fuel c, SB _{u, c, t}represent that the steam production that boiler u obtains at cycle t internal consumption fuel c, q are the operating mode run, A _{q, u, c, t}for 0-1 variable, represent whether boiler u operates (0 expression does not operate in operating mode q interval, and 1 expression operates in operating mode q interval) at cycle t in operating mode q interval, β _{q, u, c, t}for boiler u in operating mode q interval is in the linearizing continuous variable of running status of cycle t internal consumption fuel c, value is at 0 ~ 1, XFB _{q, u, c, t}and XFB _{q-1, u, c, t}represent respectively in operating mode q interval with the maximal value of boiler u in operating mode q-1 interval at cycle t internal consumption fuel c, XSB _{q, u, c, t}and XSB _{q-1, u, c, t}to be illustrated respectively in operating mode q interval and the maximal value of the steam production that boiler u obtains at cycle t internal consumption fuel c in operating mode q-1 interval.

In the present embodiment, the charging (i.e. fuel) of boiler comprises fuel oil, gas and rock gas, and wherein rock gas acquiring way is outsourcing.

The multi-state state constraint of boiler, is expressed as follows:

$\underset{q}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{A}_{q,u,c,t}\≤1,---\left(13\right)$

Represent that boiler u can only consume a kind of fuel at the most in one-period t, and can only run in the interval q of an operating mode.

Multi-state steam turbine model comprises the steam consumption of steam turbine and the multi-state state constraint model of generated energy model and steam turbine.

The steam consumption of steam turbine and generated energy model, be expressed as follows respectively:

${\mathrm{STC}}_{u,t}=\underset{q\∈Q}{\mathrm{\Σ}}(B_{q,u,t}{\mathrm{XST}}_{q-1,u,t}+{\mathrm{\γ}}_{q,u,t}({\mathrm{XST}}_{q,u,t}-{\mathrm{XST}}_{q-1,u,t})),---\left(14\right)$

${\mathrm{ET}}_{u,t}=\underset{q\∈Q}{\mathrm{\Σ}}(B_{q,u,t}{\mathrm{XET}}_{q-1,u,t}+{\mathrm{\γ}}_{q,u,t}({\mathrm{XET}}_{q,u,t}-{\mathrm{XET}}_{q-1,u,t})),---\left(15\right)$

Namely by piece wire approximation method, the steam consumption of steam turbine and steam turbine power generation amount is made to carry out multistage linear association process.STC _u,trepresent the steam consumption of steam turbine u in cycle t, ET _u,trepresent the generated energy of steam turbine u in cycle t, q is steam turbine operation operating mode, and Q is the set of steam turbine operation operating mode, B _{q, u, t}for 0-1 variable, represent whether steam turbine u operates at cycle t in operating mode q interval, γ _{q, u, t}represent the linearizing continuous variable of the running status of steam turbine u in cycle t in operating mode q interval, value is at 0 ~ 1, XST _{q, u, t}and XST _{q-1, u, t}be illustrated respectively in the maximal value of steam turbine u steam consumption in cycle t in operating mode q interval and q-1 interval, XET _{q, u, t}and XET _{q-1, u, t}be illustrated respectively in the maximal value of steam turbine u generated energy in cycle t in operating mode q interval and q-1 interval;

The multi-state state constraint of steam turbine, is expressed as follows:

$\underset{q\∈Q}{\mathrm{\Σ}}{B}_{q,u,t}\≤1,---\left(16\right)$

This multi-state state constraint represents that steam turbine u can only run at the most in one-period t in the interval q of an operating mode.

Each energy medium equilibrium of supply and demand constraint comprises the constraint of the steam equilibrium of supply and demand and retrains with the electric energy equilibrium of supply and demand.

Steam equilibrium of supply and demand constraint representation is:

$\begin{array}{c}\underset{u\∈\mathrm{UB}}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{\mathrm{SB}}_{u,c,t}+\underset{u\∈\mathrm{UT}}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{\mathrm{STG}}_{u,c,t}+\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{EG}}_{u,m,c,t}-\underset{u\∈\mathrm{UT}}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{\mathrm{STC}}_{u,c,t}+{\mathrm{EP}}_{c,t}\\ +{\mathrm{LSI}}_{c,t}-{\mathrm{LSO}}_{c,t}\≥\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{ED}}_{u,m,c,t}\end{array}---\left(17\right)$

UB is boiler set, and UT is steam turbine set, and UP is device sets; SB _{u, c, t}represent the steam production that boiler u obtains at cycle t internal consumption fuel c, STG _{u, c, t}represent the steam production that steam turbine u obtains at cycle t internal consumption fuel c, EG _{u, m, c, t}and ED _{u, m, c, t}distinguish output sum and the consumption sum of materiel machining device u each grade steam of steam c under operation scheme m in indication cycle t, STC _{u, c, t}represent that steam turbine u generates electricity in cycle t and consume the consumption of steam c, EP _c,trepresent the outsourcing amount of steam c in cycle t, LSI _c,tand LSO _c,trepresent the air inflow of the steam c of temperature-decreased pressure reducer in cycle t and air output respectively;

Electric energy equilibrium of supply and demand constraint representation is:

${\mathrm{EP}}_{c,t}+\underset{u\∈\mathrm{UT}}{\mathrm{\Σ}}{\mathrm{ET}}_{u,t}\≥\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{ED}}_{u,m,c,t},---\left(18\right)$

EP _c,trepresent the outsourcing amount of steam c in cycle t, ET _u,trepresent the generated energy of steam turbine u in cycle t, ED _{u, m, c, t}represent materiel machining device u consumption of each grade steam of steam c under operation scheme m in cycle t.

Equipment capacity constraint representation in described step (2) is:

${\mathrm{FU}}_u^{\mathrm{lo}}\≤{\mathrm{FU}}_{u,t}\≤{\mathrm{FU}}_u^{\mathrm{up}},---\left(19\right)$

${\mathrm{LS}}_{c,t}\≤{\mathrm{LS}}_{c,t}^{\mathrm{up}},---\left(20\right)$

with the lower limit of the production capacity of difference indication equipment u and the upper limit, FU _u,tthe load of indication equipment u in cycle t, LS _c,tin indication cycle t, temperature-decreased pressure reducer is by the load of the process of steam c, in indication cycle t, temperature-decreased pressure reducer is by the upper limit of the load of the process of steam c.

Fig. 4 is the relation between supply and demand interaction figure of each energy medium between production system and energy system in oil refining enterprise.Device (comprising materiel machining device and attemperation apparatus) wherein in production system is processed material (crude oil), comprises time processing process and secondary processing process etc., eventually passes harmonic process and obtain final products.Device in production system can consume different brackets steam (high pressure steam, middle pressure steam and low-pressure steam) and electricity when processing materials, also can generating portion grade steam.By consuming water, fuel (gas, fuel oil and gas in energy system, wherein outsourcing fuel mainly comprises fuel oil and rock gas) produce steam and electricity, the electricity of generation is supplied to Utilities Electric Co.'s (electrical network), steam and provides energy for production system.Boiler consume fuel producing steam (steam), steam turbine consumes steam-electric power.

(3) set up the logistics of oil refining enterprise's production system and energy system and multicycle coupling model can be flowed, comprising physical property and the public work gas discharging factor correlation model of the mass balance model of intermediate product, the product power consumption model of materiel machining device and mediation product.Specifically comprise the following steps:

(3-1) for fuel oil and gas, can as the fuel input product of oil refining enterprise's energy system, be again the intermediate product that device incidentally produces simultaneously, set up the mass balance model comprising product inventory, sale, outsourcing, mediation production and boiler and consume, as follows:

${\mathrm{MI}}_{c,t}+\underset{u}{\mathrm{\Σ}}{\mathrm{FB}}_{u,c,t}+{\mathrm{SC}}_{c,t}={\mathrm{MI}}_{c,t-1}+\underset{u}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}+{\mathrm{PC}}_{c,t},---\left(21\right)$

Wherein, MI _c,tand MI _{c, t-1}difference indication cycle t is interior and cycle t-1 fuel is oily or the tank farm stock of gas c, SC _c,tthe sales volume of indication cycle t fuel c, FB _{u, c, t}represent that in energy system, boiler u is at the consumption of cycle t internal consumption fuel c, FP _{u, m, c, t}the output of the fuel c that attemperation apparatus u obtains under operation scheme m in indication cycle t, PC _c,tthe outsourcing amount of indication cycle t fuel c;

(3-2) for materiel machining device, produce by the product consumption of public work in related product production run and the fixing of the load of materiel machining device, the physical property of product and materiel machining device the product of setting up materiel machining device of consuming energy and to consume energy model, as follows:

ED _u,m,c,t＝λ _u,m,cFF _u,m,t+ρ _u,c,pP _c,p+CD _u,m,c， (22)

EG _u,m,c,t＝μ _u,m,cFF _u,m,t+σ _u,c,pP _c,p+CG _u,m,c， (23)

In formula, ED _{u, m, c, t}represent materiel machining device u in cycle t under operation scheme m to the demand of public work c, λ _{u, m, c}represent the materiel machining materiel machining amount of device u under operation scheme m and this consumption ratio coefficient of public work c, FF _{u, m, t}the material treatment capacity of materiel machining device u under operation scheme m in indication cycle t, ρ _{u, c, p}represent the correlation coefficient of the consumption of public work c and the physical property p of product discharge in materiel machining device u, P _c,prepresent the value of the physical property p of product c, CD _{u, m, c}represent that materiel machining device u consumes the fixed charges of public work c under operation scheme m, EG _{u, m, c, t}the output of the public work c that materiel machining device u generates under operation scheme m in indication cycle t, μ _{u, m, c}represent the materiel machining materiel machining amount of device u under operation scheme m and the production ratio coefficient of public work c, σ _{u, c, p}represent the correlation coefficient of the output of public work c and the physical property p of product discharge in materiel machining device u, CG _{u, m, c}the fixing output of indication device u public work c under operation scheme m;

(3-3) for harmonic process and boiler combustion gas generation process, by physical property and the public work gas discharging factor of association mediation product, the correlation model of fuel oil production run and Pollution of Boiler gas discharge process is set up, as follows;

$\underset{c}{\mathrm{\Σ}}{\mathrm{FC}}_{u,m,c,t}{\mathrm{PI}}_{c,p}=\underset{c}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}{P}_{c,p},---\left(24\right)$

XC _u，c，g，t＝ω _gf(P _c，p)FB _u，c，t， (25)

In formula, PI _c,prepresent the value of the physical property p of fuel oil c before entering attemperation apparatus, XC _{u, c, g, t}represent the discharge capacity of the dusty gas g that the consume fuel oil c of boiler u in cycle t produces, ω _grepresent the emission factor of dusty gas g, f (P _c,p) represent the quality of dusty gas g.

(4) set up the logistics of described oil refining enterprise and multicycle Integrated Optimization Model can be flowed, utilize model decomposition solution strategies to logistics with can flow multicycle coupling model and carry out quick optimizing;

First, introduce consider material and energy cost with the integrated objective function of full factory profit maximization:

Full factory profit=profit of sales-device processing cost-purchasing of raw materials cost-stock's operation cost-fuel consumption cost-Exogenous factor medium cost-equipment operation maintenance expense-environmental emission cost.

Wherein, environmental emission cost comprises CO ₂gas discharging cost, SO _xgas discharging cost and NO _xgas discharging cost;

The all models set up in associating above-mentioned steps, set up the logistics of oil refining enterprise and can flow multicycle Integrated Optimization Model.

Then, realize the quick optimizing of integrated model by model decomposition solution strategies, Fig. 5 is that Integrated Optimization Model decomposes solution strategies schematic diagram, is divided into following steps:

(4-1) physical property of final products is set as fixed value, the multicycle Production planning model of the production system of oil refining enterprise is relaxed into MIXED INTEGER linear model (MILP);

(4-2) the input product market demand, and turn to target with full factory material production maximum profit, rapid solving is carried out to the MIXED INTEGER linear model of production system, obtains the load of each materiel machining device in production system, operation scheme, product library storage, the physical property of product and the turnout of gas and fuel oil;

(4-3) result of described step (4-2) is substituted into the product power consumption model of materiel machining device, calculate the product consumption demand of all kinds of energy medium in production system;

(4-4) consumption demand is produced as constraint of demand using all kinds of energy medium of production system, be minimised as target with the energy cost of energy system, to energy system multicycle operation optimization model (the multicycle operation optimization model of energy system is exactly MILP) solve the product consumption of load and all kinds of energy medium obtaining each equipment in energy system;

(4-5) by the fixed value of the physical property of final products, by the turnout of the physical property of the load of materiel machining device each in production system, operation scheme, product library storage, product, gas and fuel oil, and in energy system the load of each equipment and the product consumption of all kinds of energy medium as logistics and the initial feasible solution that can flow multicycle coupling model (being exactly MINLP), utilize branch-and-bound optimized algorithm, solve this logistics and multicycle coupling model can be flowed obtain the optimum material production scheme of whole oil refining enterprise and energy production consumes scheme.

(5) produce according to the result of the quick optimizing in step (4).

Setting optimizing relative tolerance is 5%, and feasible solution reaches can be held in error range, and optimizing stops automatically.Compared to tradition substep modeling optimization method, integrated moulding optimization method effect is as shown in table 1.

Table 1

Table 1 is visible, compared to existing method (step-by-step optimization model), the present embodiment (Integrated Optimization Model) not only increases the productive profit of whole enterprise, improves efficiency of energy utilization simultaneously, reduces the empty discharge capacity of gas and reduces the discharge capacity of various dusty gas.And utilize model decomposition strategy, effectively improve the Searching efficiency of integrated optimization method.

Claims (8)

1., based on a production method for the integrated optimization of production system and energy system, it is characterized in that, comprising:

(1) combine according to the production procedure of oil refining enterprise and device the multicycle Production planning model setting up production system;

(2) set up the multicycle operation optimization model of energy system, comprise multi-state boiler model, multi-state steam turbine model, the constraint of each energy medium equilibrium of supply and demand and equipment capacity and retrain;

(3) utilize the product power consumption model of the material balance relationship of intermediate product, materiel machining device and be in harmonious proportion physical property and the public work gas discharging factor correlation model of product, set up the logistics of production system and energy system and multicycle coupling model can be flowed;

(4) model set up according to step (1), (2), (3) and the integrated objective function according to material and energy cost, set up logistics and multicycle Integrated Optimization Model can be flowed, utilize model decomposition solution strategies to logistics with can flow multicycle coupling model and carry out quick optimizing;

(5) produce according to the result of the quick optimizing in step (4).

2. as claimed in claim 1 based on the production method of the integrated optimization of production system and energy system, it is characterized in that, the multicycle Production planning model in described step (1) comprises the yield model of product demand constraint in production system, material balance, stock's constraint, harmonic process constraint and each materiel machining device.

3. as claimed in claim 2 based on the production method of the integrated optimization of production system and energy system, it is characterized in that, multi-state boiler model in described step (2) comprises the Fuel Consumption of boiler and the multi-state state constraint of steam production model and boiler, wherein:

The Fuel Consumption of boiler and steam production model, as follows:

${\mathrm{FB}}_{u,c,t}=\underset{q}{\mathrm{\Σ}}(A_{q,u,c,t}{\mathrm{XFB}}_{q-1,u,c,t}+{\mathrm{\β}}_{q,u,c,t}({\mathrm{XFB}}_{q,u,c,t}-{\mathrm{XFB}}_{q-1,u,c,t})),$

${\mathrm{SB}}_{u,c,t}=\underset{q}{\mathrm{\Σ}}(A_{q,u,c,t}{\mathrm{XSB}}_{q-1,u,c,t}+{\mathrm{\β}}_{q,u,c,t}({\mathrm{XSB}}_{q,u,c,t}-{\mathrm{XSB}}_{q-1,u,c,t})),$

FB _{u, c, t}represent the consumption of boiler u at cycle t internal consumption fuel c, SB _{u, c, t}represent the steam production that boiler u obtains at cycle t internal consumption fuel c, q is operating condition, A _{q, u, c, t}for 0-1 variable, represent whether boiler u operates at cycle t in operating mode q interval, β _{q, u, c, t}for boiler u in operating mode q interval is in the linearizing continuous variable of running status of cycle t internal consumption fuel c, value is at 0 ~ 1, XFB _{q, u, c, t}and XFB _{q-1, u, c, t}represent respectively in operating mode q interval with the maximal value of boiler u in operating mode q-1 interval at cycle t internal consumption fuel c, XSB _{q, u, c, t}and XSB _{q-1, u, c, t}to be illustrated respectively in operating mode q interval and the maximal value of the steam production that boiler u obtains at cycle t internal consumption fuel c in operating mode q-1 interval;

The multi-state state constraint of boiler, is expressed as follows:

$\underset{q}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{A}_{q,u,c,t}\≤1.$

4. as claimed in claim 3 based on the production method of the integrated optimization of production system and energy system, it is characterized in that, multi-state steam turbine model in described step (2) comprises the steam consumption of steam turbine and the multi-state state constraint of generated energy model and steam turbine, wherein:

Steam consumption and the generated energy model of steam turbine are as follows:

${\mathrm{STC}}_{u,t}=\underset{q\∈Q}{\mathrm{\Σ}}(B_{q,u,t}{\mathrm{XST}}_{q-1,u,t}+{\mathrm{\γ}}_{q,u,t}({\mathrm{XST}}_{q,u,t}-{\mathrm{XST}}_{q-1,u,t})),$

${\mathrm{ET}}_{u,t}=\underset{q\∈Q}{\mathrm{\Σ}}(B_{q,u,t}{\mathrm{XET}}_{q-1,u,t}+{\mathrm{\γ}}_{q,u,t}({\mathrm{XET}}_{q,u,t}-{\mathrm{XET}}_{q-1,u,t})),$

STC _u,trepresent the steam consumption of steam turbine u in cycle t, ET _u,trepresent the generated energy of steam turbine u in cycle t, q is steam turbine operation operating mode, and Q is the set of steam turbine operation operating mode, B _{q, u, t}for 0-1 variable, represent whether steam turbine u operates at cycle t in operating mode q interval, γ _{q, u, t}represent the linearizing continuous variable of the running status of steam turbine u in cycle t in operating mode q interval, value is at 0 ~ 1, XST _{q, u, t}and XST _{q-1, u, t}be illustrated respectively in the maximal value of steam turbine u steam consumption in cycle t in operating mode q interval and q-1 interval, XET _{q, u, t}and XET _{q-1, u, t}be illustrated respectively in the maximal value of steam turbine u generated energy in cycle t in operating mode q interval and q-1 interval;

The multi-state state constraint of steam turbine is expressed as follows:

$\underset{q\∈Q}{\mathrm{\Σ}}{B}_{q,u,t}\≤1.$

5. as claimed in claim 4 based on the production method of the integrated optimization of production system and energy system, it is characterized in that, each energy medium equilibrium of supply and demand constraint in described step (2) comprises the constraint of the steam equilibrium of supply and demand and the constraint of the electric energy equilibrium of supply and demand, wherein:

Steam equilibrium of supply and demand constraint representation is:

$\underset{u\∈\mathrm{UB}}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{\mathrm{SB}}_{u,c,t}+\underset{u\∈\mathrm{UT}}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{\mathrm{STG}}_{u,c,t}+\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{EG}}_{u,m,c,t}-\underset{u\∈\mathrm{UT}}{\mathrm{\Σ}}\underset{c}{\mathrm{\Σ}}{\mathrm{STC}}_{u,c,t}+{\mathrm{EP}}_{c,t}+{\mathrm{LSI}}_{c,t}-{\mathrm{LSO}}_{c,t}\≥\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{ED}}_{u,m,c,t},$

UB is boiler set, and UT is steam turbine set, and UP is materiel machining device sets; SB _{u, c, t}represent the steam production that boiler u obtains at cycle t internal consumption fuel c, STG _{u, c, t}represent the steam production that steam turbine u obtains at cycle t internal consumption fuel c, EG _{u, m, c, t}and ED _{u, m, c, t}distinguish output sum and the consumption sum of materiel machining device u each grade steam of steam c under operation scheme m in indication cycle t, STC _{u, c, t}represent that steam turbine u generates electricity in cycle t and consume the consumption of steam c, EP _c,trepresent the outsourcing amount of steam c in cycle t, LSI _c,tand LSO _c,trepresent the air inflow of the steam c of temperature-decreased pressure reducer in cycle t and air output respectively;

Electric energy equilibrium of supply and demand constraint representation is:

${\mathrm{EP}}_{c,t}+\underset{u\∈\mathrm{UT}}{\mathrm{\Σ}}{\mathrm{ET}}_{u,t}\≥\underset{u\∈\mathrm{UP}}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{ED}}_{u,m,c,t}.$

6., as claimed in claim 5 based on the production method of the integrated optimization of production system and energy system, it is characterized in that, the equipment capacity constraint representation in described step (2) is:

${\mathrm{FU}}_u^{\mathrm{lo}}\≤{\mathrm{FU}}_{u,t}\≤{\mathrm{FU}}_u^{\mathrm{up}}$

${\mathrm{LS}}_{c,t}\≤{\mathrm{LS}}_{c,t}^{\mathrm{up}}$

with the lower limit of the production capacity of difference indication equipment u and the upper limit, FU _u,tthe load of indication equipment u in cycle t, LS _c,tin indication cycle t, temperature-decreased pressure reducer is by the load of the process of steam c, in indication cycle t, temperature-decreased pressure reducer is by the upper load limit of the process of steam c.

7., as claimed in claim 6 based on the production method of the integrated optimization of production system and energy system, it is characterized in that, described step (3) comprises the steps:

(3-1) fuel is set up to the mass balance model comprising product inventory, sale, outsourcing, mediation production and boiler and consume, as follows:

${\mathrm{MI}}_{c,t}+\underset{u}{\mathrm{\Σ}}{\mathrm{FB}}_{u,c,t}+{\mathrm{SC}}_{c,t}={\mathrm{MI}}_{c,t-1}+\underset{u}{\mathrm{\Σ}}\underset{m}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}+{\mathrm{PC}}_{c,t},$

Wherein, MI _c,tand MI _{c, t-1}with the tank farm stock of cycle t-1 fuel c in difference indication cycle t, SC _c,tthe sales volume of indication cycle t fuel c, FB _{u, c, t}represent that in energy system, boiler u is at the consumption of cycle t internal consumption fuel c, FP _{u, m, c, t}the output of the fuel c that attemperation apparatus u obtains under operation scheme m in indication cycle t, PC _c,tthe outsourcing amount of indication cycle t fuel c;

(3-2) for materiel machining device, produce by the product consumption of public work in related product production run and the fixing of the load of materiel machining device, the physical property of product and materiel machining device the product of setting up materiel machining device of consuming energy and to consume energy model, as follows:

ED _u,m,c,t＝λ _u,m,cFF _u,m,t+ρ _u,c,pP _c,p+CD _u,m,c，

EG _u,m,c,t＝μ _u,m,cFF _u,m,t+σ _u,c,pP _c,p+CG _u,m,c，

In formula, ED _{u, m, c, t}represent materiel machining device u in cycle t under operation scheme m to the demand of public work c, λ _{u, m, c}represent the materiel machining materiel machining amount of device u under operation scheme m and this consumption ratio coefficient of public work c, FF _{u, m, t}the material treatment capacity of materiel machining device u under operation scheme m in indication cycle t, ρ _{u, c, p}represent the correlation coefficient of the consumption of public work c and the physical property p of product discharge in materiel machining device u, P _c,prepresent the value of the physical property p of product c, CD _{u, m, c}represent that materiel machining device u consumes the fixed charges of public work c under operation scheme m, EG _{u, m, c, t}the output of the public work c that materiel machining device u generates under operation scheme m in indication cycle t, μ _{u, m, c}represent the materiel machining materiel machining amount of device u under operation scheme m and the production ratio coefficient of public work c, σ _{u, c, p}represent the correlation coefficient of the output of public work c and the physical property p of product discharge in materiel machining device u, CG _{u, m, c}the fixing output of indication device u public work c under operation scheme m;

(3-3) for harmonic process and boiler combustion gas generation process, by physical property and the public work gas discharging factor of association mediation product, the correlation model of fuel oil production run and Pollution of Boiler gas discharge process is set up, as follows;

$\underset{c}{\mathrm{\Σ}}{\mathrm{FC}}_{u,m,c,t}{\mathrm{PI}}_{c,p}=\underset{c}{\mathrm{\Σ}}{\mathrm{FP}}_{u,m,c,t}{P}_{c,p},$

XC _u,c,g,t＝ω _gf(P _c,p)FB _u,c,t，

In formula, PI _c,prepresent the value of the physical property p of fuel oil c, XC _{u, c, g, t}represent the discharge capacity of the dusty gas g that the consume fuel oil c of boiler u in cycle t produces, ω _grepresent the emission factor of dusty gas g, f (P _c,p) represent the quality of dusty gas g.

8., as claimed in claim 7 based on the production method of the integrated optimization of production system and energy system, it is characterized in that, described step (4) comprises the steps:

(4-1) physical property of final products is set as fixed value, the multicycle Production planning model of the production system of oil refining enterprise is relaxed into MIXED INTEGER linear model;

(4-2) the input product market demand, and turn to target with full factory material production maximum profit, rapid solving is carried out to the MIXED INTEGER linear model of production system, obtains the load of each materiel machining device in production system, operation scheme, product library storage, the physical property of product and the turnout of gas and fuel oil;

(4-3) result of described step (4-2) is substituted into the product power consumption model of materiel machining device, calculate the product consumption demand of all kinds of energy medium in production system;

(4-4) consumption demand is produced as constraint of demand using all kinds of energy medium of production system, be minimised as target with the energy cost of energy system, to energy system multicycle operation optimization model solve the product consumption of load and all kinds of energy medium obtaining each equipment in energy system;

(4-5) by the fixed value of the physical property of final products, by the turnout of the physical property of the load of materiel machining device each in production system, operation scheme, product library storage, product, gas and fuel oil, and in energy system the load of each equipment and the product consumption of all kinds of energy medium as logistics and the initial feasible solution that can flow multicycle coupling model, utilize branch-and-bound optimized algorithm, solve this logistics and multicycle coupling model can be flowed obtain the optimum material production scheme of whole oil refining enterprise and energy production consumes scheme.